它三口采 油井，向岩心中注水，至压力平稳，记录平稳压力，计算岩心沿此采出井方向基质渗透率，另 外三口井也重复上面操作测定渗透率；
[0045] ④进行饱和油操作，在注入井用恒流栗注入模拟油，待任意一口采出井不再出水 并恒定出油之后关闭该井，直至全部井位饱和完全结束，计算饱和油体积，并在81°C条件下 老化超过72小时。
[0046] (三）实验部分
[0047] ①待充分老化后，以lmL/min的流量进行一次水驱，各井用试管接采出液，计算含 水率变化，同时水驱过程中每隔半个小时把所有采出井出口端放在同一试管测定岩心综合 含水率，至综合含水率达到98%，一次水驱结束；
[0048] ②以0. 5mL/min的流量注入0. 1PV交联聚合物溶液，充分候凝，实验过程中各井用 试管接采出液，计算含水率变化，同时监测各采压点压力变化，并密集取样，每半个小时测 定岩心综合含水率；
[0049] ③以lmL/min的流量进行后续水驱，水驱过程中监测各采压点压力变化，并密集 取样，同时每隔半小时测定综合岩心综合含水率，至达到98%，实验结束；
[0050] ④测定各测压点和各井取样中聚合物浓度，拟合曲线，如图7-9所示，浓度值的确 定用的是淀粉-碘化镉法测量取样点采出液聚合物浓度值。
[0051] 分别在（10 :00-10 :15,14 :00-14 :15,18 :00-18 :15, 22 :00-22 :15, 2 :00-2 :15) 打开1号、3号、5号、7号内环测压点，取液15分钟；（10 :15-10 :30,14 :15-14 :30,18 : 15-18 :30, 22 :15-22 :30, 2 :15-2 :30)打开2号、4号、6号、8号外环测压点，取液15分钟； 汇总绘制注胶压力分布图，如图2-6所示。
[0052] 2.根据各取样点上聚合物溶液值采用克里金插值算法反演出模型其他部位的聚 合物浓度值。
[0053] 1)以得到各取样点的聚合物浓度值作为原始数据，对该原始数据进行特高值处 理，检验是否服从正态分布；
[0054] 2)计算实验半变异函数，η个方向；
[0055] 变异函数是地质统计学所具有的基本工具，它既可以用来描述区域化变量的空间 结构变化，又可以描述其随机变化。假设空间点X只在一维的X轴上变化，把区域化变量 Ζ(χ)在x、x+h两点处的值之差的一半定义为Ζ(χ)在X轴方向上的半变异函数，即：
[0056] γ(x,h) =E(Z(x)-Z(x+h) )/2
[0057] 在满足平稳假设时，单个方向上求实验半变异函数yYh)的计算公式为：
[0058]
[0059] 式中Z(Xl)-Xl处的观测值；
[0060] N(h) 一计算数值点的个数
[0061] 通过计算yYh)值作图，选取理论模型（球状、指数模型等），并求出基台，变程。
[0062] 3)根据待估计区域的信息情况进行超块和次块的划分，并就每行次块的起始与终 止号建立文件，并对对已知信息点进行归块。
[0063] 4)以待估计超块为中心，向外围增加若干个超块作为一个邻域，本实施例中为8 个，进行滑动估计；
[0064] 5)根据参数基台，变程及选取的模型，求γ(Xl，x])或r(Xl，V)，建立方程组。
[0065] 线性估计量为：
[0066]
[0067] 其中4为被估计区域真值Zv的估计值；λi为克里金权系数；Z 估计值邻域内 的观测值；要使估计值霉最佳，必须满足下列两个条件：
[0068] 无偏性
[0069]
[0070] 估计方差σ/最小[0071] 通过方程组：
[0076] 其中λi为克里金权系数；γ(XXj)为观测值xgX名间的半变异函数值；
[0072]
[0073]
[0074]
[0075] r(Xl，V)为观测值与被估计区V内所有点之间的平均变异函数；r(V，V)为被估计区V内所 有点的平均半变异函数；μ为拉格朗日常数。
【主权项】
1. 一种注入聚合物推进的前缘模拟监测方法，其特征在于，该监测方法包括W下步 骤： 1) 在=维物理诱注实验模型取样点上采集聚合物溶液浓度值； 2) 针对各取样点上聚合物溶液浓度值采用克里金插值算法预测模型其他部位的聚合 物浓度值； 3)根据所预测出模型其他部位的聚合物浓度值实现对聚合物推进前缘模拟的监测。2. 根据权利要求1所述的注入聚合物推进的前缘模拟监测方法，其特征在于，所述步 骤2)中克里金算法的计算过程如下： A. W取样点上采集的聚合物溶液浓度值为原始数据，在单个方向上求解实验半变异函 数丫"^化)，通过计算丫^^化）值作图，选取理论模型，求出基台和变程； B.根据待估计区域的信息情况进行超块和次块的划分，对已知信息点进行归块； C. W待估计超块为中心向外围增加若干个超块作为一个邻域，进行滑动估计； 化根据参数基台，变程及选取的模型，建立克里金方程； E.W无偏性和最小估计方差为约束调节，求解克里金方程，使估计值《最佳，求解得 到结果即为待估计区域3.根据权利要求1所述的注入聚合物推进的前缘模拟监测方法，其特征在于，所述半 变异函数丫^^化）的公式为：其中Z(Xi)为Xi处的观测值，即在取样点X1处采集的聚合物溶液浓度值；NOi)为计算 数值点的个数，即取样点的个数。4.根据权利要求3所述的注入聚合物推进的前缘模拟监测方法，其特征在于，所述步 骤D中所采用的方程为：其中^为被估计区域真值Zv的估计值；A1为克里金权系数；Zi为估计值邻域内的观 测值。5.根据权利要求1所述的注入聚合物推进的前缘模拟监测方法，其特征在于，所述步 骤1)是采用平面非均质岩屯、来模拟储层平面非均质性，在采出井方向设置高渗条带，注入 强度符合要求的调驱体系聚合物进行驱油，沿程在各取样点采集聚合物浓度值并监测压 力。
【专利摘要】本发明涉及一种注入聚合物推进的前缘模拟监测方法，属于原油采收技术领域。本发明首先通过在实验模型取样点上采集聚合物溶液浓度值；然后针对各取样点上聚合物溶液浓度值采用克里金插值算法预测模型其他部位的聚合物浓度值；根据所预测出模型其他部位的聚合物浓度值实现对聚合物推进前缘模拟的监测。本发明利用数学地质的方式，运用克里金插值方法，根据一个区域内外若干信息样品的浓度值，对该区域作出一种线性无偏和最小估计方差的估计，从而实现对整个区域注入聚合物推进前缘模拟的监测，并且可以适用于交联聚合物推进前缘模拟的监测。
【IPC分类】E21B47/10, G06F17/50
【公开号】CN105373643
【申请号】CN201510537423
【发明人】林杨, 肖磊, 张丽庆, 郭艳, 权国胜, 王正欣, 孔昭柯, 刘艳华
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司河南油田分公司勘探开发研究院
【公开日】2016年3月2日
【申请日】2015年8月27日