CN108868715B - 运用分形理论确定基质-高渗条带化学驱剩余油分布及运移规律的方法 - Google Patents
运用分形理论确定基质-高渗条带化学驱剩余油分布及运移规律的方法 Download PDFInfo
- Publication number
- CN108868715B CN108868715B CN201810559938.5A CN201810559938A CN108868715B CN 108868715 B CN108868715 B CN 108868715B CN 201810559938 A CN201810559938 A CN 201810559938A CN 108868715 B CN108868715 B CN 108868715B
- Authority
- CN
- China
- Prior art keywords
- residual oil
- image
- box
- grid
- counting dimension
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000013508 migration Methods 0.000 title claims abstract description 35
- 230000005012 migration Effects 0.000 title claims abstract description 35
- 239000000126 substance Substances 0.000 title claims abstract description 30
- 230000035699 permeability Effects 0.000 title description 2
- 238000006073 displacement reaction Methods 0.000 claims abstract description 41
- 239000011148 porous material Substances 0.000 claims abstract description 35
- 230000008859 change Effects 0.000 claims abstract description 21
- 230000008569 process Effects 0.000 claims abstract description 10
- 238000002474 experimental method Methods 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims abstract description 3
- 238000002347 injection Methods 0.000 claims description 55
- 239000007924 injection Substances 0.000 claims description 55
- 150000001875 compounds Chemical class 0.000 claims description 34
- 238000004519 manufacturing process Methods 0.000 claims description 24
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 17
- 239000011159 matrix material Substances 0.000 claims description 13
- 230000000007 visual effect Effects 0.000 claims description 11
- 238000000605 extraction Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 239000007788 liquid Substances 0.000 claims description 4
- 229920006395 saturated elastomer Polymers 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000003709 image segmentation Methods 0.000 claims description 2
- 238000004445 quantitative analysis Methods 0.000 claims 1
- 239000000758 substrate Substances 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 8
- 230000007547 defect Effects 0.000 abstract description 2
- 238000000638 solvent extraction Methods 0.000 abstract 1
- 239000011521 glass Substances 0.000 description 11
- 238000001259 photo etching Methods 0.000 description 6
- RTZKZFJDLAIYFH-UHFFFAOYSA-N Diethyl ether Chemical compound CCOCC RTZKZFJDLAIYFH-UHFFFAOYSA-N 0.000 description 4
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 4
- KRHYYFGTRYWZRS-UHFFFAOYSA-N Fluorane Chemical compound F KRHYYFGTRYWZRS-UHFFFAOYSA-N 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 238000002360 preparation method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 238000009738 saturating Methods 0.000 description 3
- 238000004140 cleaning Methods 0.000 description 2
- 230000007797 corrosion Effects 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 2
- 239000006059 cover glass Substances 0.000 description 2
- 239000012153 distilled water Substances 0.000 description 2
- 238000002296 dynamic light scattering Methods 0.000 description 2
- 239000008398 formation water Substances 0.000 description 2
- 239000003292 glue Substances 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 229920002120 photoresistant polymer Polymers 0.000 description 2
- 229920000642 polymer Polymers 0.000 description 2
- 230000032683 aging Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/194—Segmentation; Edge detection involving foreground-background segmentation
Landscapes
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种运用分形理论确定基质‑高渗条带化学驱剩余油分布及运移规律的方法。主要目的在于克服剩余油运移规律研究方法过于单一的缺点，本发明提出采用分形理论确定基质‑高渗条带化学驱剩余油分布及运移规律的方法。具体解决方案为：制作微观物理模型，通过微观驱油装置进行基质‑高渗条带化学驱驱油实验，收集渗流过程的连续图像，根据实际情况将驱替范围分区，通过图像处理，去除实际孔隙骨架，提取完整剩余油分布形态，按分区网格计算网格内剩余油计盒维数，总结相同网格及相邻网格计盒维数变化曲线，对比相邻网格计盒维数变化规律，确定基质‑高渗条带化学驱剩余油分布及运移规律。
Description
技术领域
本发明隶属于油气田开发领域，尤其涉及到一种分形理论研究剩余油分布规律的装置与方法。
背景技术
随着油田开发的不断深入，国内大多数油田皆已进入高含水、高采出程度的“双高”阶段，针对二次采油未能采出的未波及区的剩余油和波及区的残余油，认识剩余油为油田二次采油及三次采油提供重要依据尤为重要。剩余油又分为宏观剩余油和微观剩余油两大类。目前微观方面关于剩余油的研究更侧重微观驱替机理及剩余油分布形态。更多的是采用微观可视模型，虽然直观，但是驱替过程中多区域同时流动，运移规律显得抽象，缺少对确定基质-高渗条带化学驱剩余油分布及运移规律的方法。此外，剩余油研究方法较为单一，缺少用直观的规律曲线来分析剩余油分布及运移规律的方法。分形几何理论研究的是具有统计自相似性的对象，自诞生以来，在理论研究和实际应用方面都得到了快速发展，已有研究将分形理论用在孔隙结构、裂缝评价、储层非均质性描述，但关于剩余油方面目前只有计盒维数的计算公式，尚无具体运用分形理论研究剩余油分布及运移规律的方法。
发明内容
为了解决背景技术中所提到的技术问题，本发明提供一种运用分形理论确定基质-高渗条带化学驱剩余油分布及运移规律的方法，该种方法将分形理论运用到剩余油分布研究领域，采用分形计盒维数定量化确定基质-高渗条带化学驱剩余油分布及运移规律，从而克服现有技术中的缺陷。
本发明的技术方案是：该种运用分形理论确定化学驱剩余油分布及运移规律的方法，其特征在于包括如下步骤：
发明步骤一：布有基质-高渗条带的微观可视模型的制备
微观模型制备的基本原理是利用曝光技术将设计好的储层孔隙结构转移到涂有感光胶的玻璃片上并显影，然后利用氢氟酸腐蚀将显影部分刻蚀出具有一定深度的孔隙结构，而其他部分因分别被光刻胶和蜡保护而不被腐蚀，最后盖上盖玻片，这样便得到可视化驱替用的光刻玻璃岩心微观模型。根据上述原理，设制长4.5mm宽4.5mm的微观可视模型，并在对角线处布置基质高渗条带，其最大宽度不超过对角线长度的5分之1。基质高渗条带两侧的喉道半径为1-2μm，高渗条带处喉道半径为4-5μm。
发明步骤二：将步骤一中获得的微观物理模型置入微观驱油装置中进行微观驱油实验，采集化合物驱替过程图像
根据上述所制模型，实验步骤如下：
（1）用乙醇、石油醚、蒸馏水清洗光刻玻璃模型；
（2）将模型抽真空，饱和所模拟地层水；
（3）45℃条件下模型饱和油，并在此温度下静置24h；
（4）以0.01ml/min的驱替速度向模型中注入水直至含水率98%，后注入化合物至残余油饱和度；
（5）通过数字显微摄像系统采集驱替全程图像。实验装置如图1所示。
微量泵通过管线和注水活塞容器的下部阀门连接，活塞容器的上部阀门通过管线和光刻玻璃岩心模型入口的控制阀门连接，光刻玻璃岩心模型出口的控制阀门通过管线和量筒连接。所述微量泵给整个驱替装置提供动力。所述注水活塞容器为注入水的容器。所述注化合物活塞容器为注入化合物的容器。所述光刻玻璃为与实际储层接近的模型。所述量筒承接采出液。所述六通提供多个通路开关控制。所述恒温箱使整个实验流程保持在地层温度下。
发明步骤三：将步骤二中采集到的化合物驱替过程图像进行剩余油图像提取，具体实现路径：（1）将步骤二中采集到的化学驱阶段RGB图像通过photoshop软件转化为灰度图像，并对图像进行亮度、对比度、锐化等预处理；（2）图像分割当中最重要的部分就是如何选取阈值,如果阈值的选择不合理,很容易将目标对象（剩余油）理解成为背景图像,或者将背景图像理解为目标对象。
当孔隙喉道半径与化合物分子回旋半径之比大于5倍时，注入低相对分子质量聚合物不会发生堵塞，小于5倍时会出现堵塞。通过动态光散射仪测量出化合物分子的回旋半径并与设计出的喉道半径比较，将化合物无法通过的孔喉通过photoshop软件去除微观孔喉并采用迭代阀值法通过photoshop软件分割图像中的剩余油，提出剩余油分布完整图像。
迭代阈值法的实现步骤如下:
①求出渗流图像中的最大灰度值Zmax和最小灰度值Zmin，并设置初始值，T0=(Zmin+Zmax)/2；
②根据阈值Tk(k=0,1,2,3…)将渗流图像分割成目标图像和背景图像两部分，再按照公式(1)和公式(2)求出这两部分的平均灰度值ZO和ZB；
式中，Z(i,j)为图像中像素点(i,j)的灰度值；N(I,j)为像素点(i,j)的加权系数，一般为1。
③求出新的阈值，Tk+1=(ZO+ZB)/2；
④如果Tk= Tk+1，则迭代结束，否则k=k+1,重复②—④。
选取上述阀值即可分割剩余油。
发明步骤四：通过MATLAB软件对第三步获取的剩余油图像网格化处理并计算初始饱和剩余油图像的计盒维数，具体路径如下：
（1）设置像素点，图像大小为A×B个像素点；
（2）划分网格：用大小为n×n网格去覆盖剩余油全部图像。n是整数，则分割比率为r=n/A，原图被分成了A/n个网格块，其中n序列为区域图像A和B的最大公约数的全部因数；
（3）被占用的判定规则：任何一个网格块只要包含剩余油就被认为是被占用的网格块，统计总共被占的网格块数并记为Nr，此时就获得了一组(r，Nr)；
（4）计算剩余油占有网格数：改变网格大小，重新获取r，Nr，具体步骤重复（2）、（3）过程，就能得到一系列(r，Nr)；
（5）计盒维数的计算：对（ln1/r，lnNr）进行线性拟合，运用最小二乘法，其中直线的斜率值D就是该图像中剩余油的计盒维数。
发明步骤五：通过Photoshop软件对步骤四中已进行完网格化处理的剩余油图像九等分；
其中，通过Photoshop软件添加图像长和宽的三等分点所在直线将图像分成九个区域；靠近注入介质端所在区域为注入端，靠近采出液体端所在区域为采出端，图像中心所在区域为注采中部，这三个区域为主流线区；所述主流线区的两侧区域为两翼区，其中，靠近注入端的区域为近井端，靠近采出端的区域为远井端，靠近注采中部的区域为盲端；在所述盲端与注采中部的边界上任意选取两处作为特征区，用于观察主流线区域是否波及到盲端；
发明步骤六：对步骤五中得到的主流线区域、两翼区域、主流线与两翼之间的区域以及特征区，计算这些区域的计盒维数并绘制计盒维数随PV数变化曲线；
计算该区域的计盒维数具体路径如下：
（1）设置像素点，图像大小为A×B个像素点；
（2）划分网格：用大小为n×n网格去覆盖剩余油全部图像。n是整数，则分割比率为r=n/A，原图被分成了A/n个网格块，其中n序列为区域图像A和B的最大公约数的全部因数；
（3）被占用的判定规则：任何一个网格块只要包含剩余油就被认为是被占用的网格块，统计总共被占的网格块数并记为Nr，此时就获得了一组(r，Nr)；
（4）计算剩余油占有网格数：改变网格大小，重新获取r，Nr，具体步骤重复（2）、（3）过程，就能得到一系列(r，Nr)；
（5）计盒维数的计算：对（ln1/r，lnNr）进行线性拟合，运用最小二乘法，其中直线的斜率值D就是该图像的计盒维数
发明步骤七：利用经过步骤六所获得的四组曲线分析基质化合物驱剩余油分布及运移规律，确定的具体路径如下：
（1）利用对应主流线区域的计盒维数随 PV 数变化曲线图按照如下规律确定基质化学驱条件下主流线剩余油分布及运移规律，
设定△R=R孔隙 —R维数 ，其中R孔隙为孔隙的分形维数与x轴之间的长度；R维数为剩余油维数与x轴之间的长度。△R反映了区域的驱替程度。当△R介于1.1-1.3之间时，剩余油呈孤立状，零散的分布在孔喉中，说明驱替程度完全。当△R处于1-1.1时，剩余油呈条带状形态分布在少数连通的孔喉中，此时剩余油多数被驱替。
（2）利用对应两翼区域的计盒维数随 PV 数变化曲线图按照如下规律确定基质化学驱条件下两翼剩余油分布及运移规律，即：当△R处于1-1.1时，部分孔喉空间充满剩余油,呈斑块状形态分布，此时剩余油绝大多数被驱替。当△R处于0.8-1时，剩余油呈条带状形态分布在少数连通的孔喉中，此时剩余油多数被驱替。当△R处于1-1.1时，部分孔喉空间充满剩余油,呈斑块状形态分布
（3）利用对应主流线与两翼之间的区域的计盒维数随PV数变化曲线图，定量分析主流线与两翼间剩余油分布及运移规律；
（4）利用对应特征区的计盒维数随 PV 数变化曲线图，按照如下规律确定基质化学驱条件下特征区剩余油分布及运移规律，即:当计盒维数减小时，说明水驱波及到侧翼。
本发明具有如下有益效果：为了克服剩余油运移规律研究方法过于单一的缺点，本发明提出采用分形理论研究剩余油分布规律。在本发明中，制作微观物理模型，通过微观驱油装置进行驱油实验，收集渗流过程的连续图像，根据实际情况将驱替范围分区，通过图像处理，去除实际孔隙骨架，提取完整剩余油分布形态，按分区网格计算网格内剩余油计盒维数，总结相同网格及相邻网格计盒维数变化曲线，对比相邻网格计盒维数变化规律，获取剩余油的分布及运移规律。本发明通过计算剩余油分布的分形维数，能够准确描述不同区域剩余油分布的复杂程度，分形维数越大，剩余油分布越复杂，驱替效果越不好，通过分形维数曲线的变化趋势，来分析驱替过程中剩余油的分布及运移规律，对油田开发中的开发方案具有重要意义。
附图说明：
图1是本发明具体实施时所采用的驱替实验装置图。
图2是本发明具体实施时采集到的原始图像。
图3是本发明具体实施时锐化处理后的图像。
图4是本发明具体实施时，去骨架后完整剩余油图像的提取图。
图5是本发明具体实施时剩余油图像九等分的区域划分示意图。
图6是本发明具体实施时主流线区维数随注入PV数变化曲线图。
图7是本发明具体实施时侧翼区域维数随注入PV数变化曲线图。
图8是本发明具体实施时主流线及侧翼相邻区域维数随注入PV数变化曲线图。
图9是本发明具体实施时特征区维数随注入PV数变化曲线图。
图中1-注入端；2-注采中部；3-采出端；4、7分别为近井两翼；5、8分别为盲端；6、9分别为远井两翼；10、11分别为特征区；12-微量泵；13-注水活塞容器；14-摄像头；15-光刻玻璃岩心模型；16-计算机；17-量筒；18-六通；19-恒温箱；20-注化合物活塞容器。
具体实施方式：
下面结合附图对本发明作进一步说明：
发明步骤一：布有基质-高渗条带的微观可视模型的制备
微观模型制备的基本原理是利用曝光技术将设计好的储层孔隙结构转移到涂有感光胶的玻璃片上并显影，然后利用氢氟酸腐蚀将显影部分刻蚀出具有一定深度的孔隙结构，而其他部分因分别被光刻胶和蜡保护而不被腐蚀，最后盖上盖玻片，这样便得到可视化驱替用的光刻玻璃岩心微观模型。根据上述原理，设制长4.5mm宽4.5mm的微观可视模型，并在对角线处布置基质高渗条带，其最大宽度不超过对角线长度的5分之1。基质高渗条带两侧的喉道半径为1-2μm，高渗条带处喉道半径为4-5μm。
步骤二：将步骤一中获得的微观物理模型置入微观驱油装置中进行微观驱油实验， 采集化合物驱替过程图像
实验步骤如下：
（1）用石油醚、乙醇、蒸馏水清洗微观仿真玻璃刻蚀模型；
（2）将清洗干净后的玻璃可视模型抽真空后，饱和模拟地层水；
（3）20℃条件下在模型中饱和油，并在此温度条件下静置老化24h；
（4）以0.01ml/min的驱替速度向模型中注入水直至含水率98%，后注入化合物至残余油饱和度；
（5）分析微观驱替图像。
步骤三：将采集到的化合物驱替全过程图像进行剩余油图像提取.
（1）将步骤2中采集到的化学驱阶段的RGB图像通过photoshop软件转化为灰度图像，并对图像进行亮度、对比度、锐化等预处理如图2、图3所示；
（2）通过分割当中最重要的部分就是如何选取阈值,如果阈值的选择不合理,很容易将目标对象（剩余油）理解成为背景图像,或者将背景图像理解为目标对象。
当孔隙喉道半径与化合物分子回旋半径之比大于5倍时，注入低相对分子质量聚合物不会发生堵塞，小于5倍时会出现堵塞。通过动态光散射仪测量出化合物分子的回旋半径并与设计出的喉道半径比较，将化合物未能通过的孔喉通过photoshop软件去除微观孔喉并采用迭代阀值法通过photoshop软件分割图像中的剩余油，提出剩余油分布完整图像。
迭代阈值法的实现步骤如下:
①求出渗流图像的最大灰度值Zmax和最小灰度值Zmin，并设置初始值，T0=(Zmin+Zmax)/2；
②根据阈值Tk(k=0,1,2,3…)将渗流图像分割成目标图像和背景图像两部分，再按照公式(1)和公式(2)求出这两部分的平均灰度值ZO和ZB；
式中，Z(i,j)为图像中像素点(i,j)的灰度值；N(l,j)为像素点(I,j)的加权系数,一般为1。
③求出新的阈值，Tk+1=(ZO+ZB)/2；
④如果Tk= Tk+1，则迭代结束,否则k=k+1，重复②—④。
经上述计算图像阀值T为28，通过photoshop软件选取上述阀值，即可得到分割的剩余油图像。如图4所示。
发明步骤四：通过MATLAB软件对第三步获取的剩余油图像网格化处理并计算初始饱和剩余油图像的计盒维数。
（1）设置像素点，图像大小为240×180个像素点；
（2）划分网格：用大小为n×n网格去覆盖剩余油全部图像。n是整数，则分割比率为r=n/A原图被分成了A/n个网格块，其中n序列为区域图像240和180的最大公约数的全部因数即（1，2，3，4，5，6，10，15，20，30，60）；
（3）被占用的判定规则：任何一个网格块只要包含剩余油就被认为是被占用的网格块，统计总共被占的网格块数并记为Nr，此时就获得了一组(r，Nr)；
（4）计算剩余油占有网格数：改变网格大小，重新获取r，Nr，具体步骤重复（2）、（3）过程，就能得到一系列(r，Nr)；
（5）计盒维数的计算：对（ln1/r，lnNr）进行线性拟合，运用最小二乘法，其中直线的斜率值D就是图像中剩余油的计盒维数。
发明步骤五：通过Photoshop软件对步骤四中已进行完网格化处理的剩余油图像九等分；
其中，通过Photoshop软件添加图像长和宽的三等分点所在直线将图像分成九个区域；靠近注入介质端所在区域为注入端，靠近采出液体端所在区域为采出端，图像中心所在区域为注采中部，这三个区域为主流线区；所述主流线区的两侧区域为两翼区，其中，靠近注入端的区域为近井端，靠近采出端的区域为远井端，靠近注采中部的区域为盲端；在所述盲端与注采中部的边界上任意选取两处作为特征区，用于观察主流线区域是否波及到盲端；上述区域如图5所示。
发明步骤六：对步骤五中得到的主流线区域、两翼区域、主流线与两翼之间的区域以及特征区，计算这些区域的计盒维数并绘制计盒维数随PV数变化曲线；
计算该区域的计盒维数具体路径如下：
（1）设置像素点，图像大小为A×B个像素点；
（2）划分网格：用大小为n×n网格去覆盖剩余油全部图像。n是整数，则分割比率为r=n/A，原图被分成了A/n个网格块，其中n序列为区域图像A和B的最大公约数的全部因数；
（3）被占用的判定规则：任何一个网格块只要包含剩余油就被认为是被占用的网格块，统计总共被占的网格块数并记为Nr，此时就获得了一组(r，Nr)；
（4）计算剩余油占有网格数：改变网格大小，重新获取r，Nr，具体步骤重复（2）、（3）过程，就能得到一系列(r，Nr)；
（5）计盒维数的计算：对（ln1/r，lnNr）进行线性拟合，运用最小二乘法，其中直线的斜率值D就是该图像的计盒维数。主流线区域的计盒维数并随PV数变化曲线如图6所示。两翼区域的计盒维数随PV数变化曲线如图7。主流线与两翼的相邻区域的计盒维数随PV数变化曲线如图8。特征区的计盒维数随PV数变化曲线如图9。
发明步骤七：利用经过步骤六所获得的四组曲线分析基质化合物驱剩余油分布及运移规律，确定的具体路径如下：
（1）利用对应主流线区域的计盒维数随 PV 数变化曲线图按照如下规律确定基质化学驱条件下主流线剩余油分布及运移规律，设定△R=R孔隙 —R维数 ，其中R孔隙为孔隙的分形维数与x轴之间的长度；R维数为剩余油维数与x轴之间的长度。△R反映了区域的驱替程度。当△R介于1.1-1.3之间时，剩余油呈孤立状，零散的分布在孔喉中，说明驱替程度完全。当△R处于1-1.1时，剩余油呈条带状形态分布在少数连通的孔喉中，此时剩余油多数被驱替。如图6所示为主流线上区域计盒维数的变化曲线，化学驱结束后，注入端与注采中部△R1、△R2 介于1.1-1.3之间，说明剩余油呈孤立状，化学驱替完全。采出端△R3 介于1-1.1之间，说明剩余油斑块状，绝大多数剩余油被驱替。随着化学驱的进行，注入端计盒维数在注入量0至0.1PV区间略微下降，说明注入端在水驱过程中绝大部分剩余油已被驱替，注入端少量的剩余油受化合物驱替导致剩余油不断运移，注入端剩余油不断减少，在注入量为0.1PV之后，计盒维数趋于平缓，说明注入端已无化合物可动剩余油，化合物已无动用潜力。而注采中部在注入端维数下降中计盒维数上下波动，说明注采中部的剩余油受到化合物的波及以及注入端剩余油的运移，计盒维数在注入量为0.1PV至0.2PV区间快速下降，说明化合物已进入注采中部驱替剩余油，注入量在0.2PV后，计盒维数趋于平缓，说明注采中部已无可动剩余油，化学驱再无动用潜力。采出端计盒维数在0.2PV前上下波动，说明采出端受化合物的波及以及注采中部剩余油的运移，在0.2PV至0.4PV维数大幅下降，说明化合物进入并驱替了采出端的剩余油，在注入量为0.4PV后1、2、3区域计盒维数趋于平缓，说明化合物已无动用潜力。
（2）利用对应两翼区域的计盒维数随 PV 数变化曲线图按照如下规律确定基质化学驱条件下两翼剩余油分布及运移规律，即：当△R处于1-1.1时，部分孔喉空间充满剩余油,呈斑块状形态分布，此时剩余油绝大多数被驱替。当△R处于0.8-1时，剩余油呈条带状形态分布在少数连通的孔喉中，此时剩余油多数被驱替。当△R处于1-1.1时，部分孔喉空间充满剩余油,呈斑块状形态分布。如图7所示为侧翼区域计盒维数的变化曲线，化学驱结束后，近井端与远井端△R4、△R6 介于1 -1.1之间，说明剩余油呈斑块状，绝大多数剩余油被驱替，而盲端△R5 介于0.8-1之间，说明剩余油呈条带状，部分剩余油未被驱替。随着化学驱的进行，近井端在注入量为0至0.2PV维数大幅下降，说明化合物进入近井端驱替剩余油导致剩余油发生运移，在注入量为0.2PV后，计盒维数趋于平缓，说明近井端已无化学驱可动剩余油。盲端的计盒维数在注入量为0至0.2PV区间上下波动，说明盲端受化合物的波及以及近井端剩余油的运移，在注入量0.2PV至0.35PV区间计盒维数大幅下降，说明化合物已注入到盲端，驱替剩余油，后在0.35PV时盲端已无化学驱可动剩余油。远井端在注入量为0至0.3PV区间时，计盒维数上下波动，说明受化合物波及以及盲端剩余油的运移，在注入量为0.3PV至0.35PV区间计盒维数快速下降，说明剩余油被大幅动用，在0.35PV时各区计盒维数趋于平缓，说明化学驱再无动用潜力。
（3）利用对应主流线与两翼之间的区域的计盒维数随PV数变化曲线图，定量分析主流线与两翼间剩余油分布及运移规律；如图8所示，两翼较注采中部计盒维数变化幅度大且起始值高于注采中部。说明水驱结束以后，侧翼残余油高于注采中部，化合物有效进入侧翼驱替剩余油，两翼剩余油动用程度较注采中部高。
（4）利用对应特征区的计盒维数随 PV 数变化曲线图，按照如下规律确定基质化学驱条件下特征区剩余油分布及运移规律，即:当计盒维数减小时，说明水驱波及到侧翼。如图9所示特征区10在注入量0至0.05PV区间的计盒维数为0，说明主流线区域没有剩余油运移至特征区10，而在注入量为0.05PV至0.15PV区间和0.15PV至0.35PV区间计盒维数先增后降至为0，说明主流线剩余油运移至特征点10，后又被驱替至特征区无剩余油。特征区11同理。
Claims (1)
1.一种运用分形理论确定基质-高渗条带化学驱剩余油分布及运移规律的方法，其特征在于包括如下步骤：
步骤一，布有基质-高渗条带的微观可视模型的制备；
步骤二，将步骤一中获得的微观物理模型置入微观驱油装置中进行微观驱油实验，采集化合物驱替驱替过程图像；
步骤三，将步骤二中采集到的化合物驱替过程图像进行剩余油图像提取；具体路径如下：
(1)将步骤二中采集到的驱替全程图像的RGB格式图像通过photoshop软件转化为灰度图像并对图像进行预处理,所述预处理的步骤包括进行亮度、对比度和锐化的处理以使得剩余油和孔隙能够得到清晰显示；
(2)采用迭代阀值法通过photoshop软件，对于预处理后的图像分割图像中的剩余油，提出剩余油分布完整图像；
步骤四，通过MATLAB软件对第三步获取的剩余油图像网格化处理并计算初始饱和剩余油图像的计盒维数，具体路径如下：
(1)设置像素点，图像大小为A×B个像素点；
(2)划分网格：用大小为n×n网格去覆盖剩余油全部图像；n是整数，则分割比率为r＝n/A，原图被分成了A/n个网格块，其中n序列为区域图像A和B的最大公约数的全部因数；
(3)被占用的判定规则：任何一个网格块只要包含剩余油就被认为是被占用的网格块，统计总共被占的网格块数并记为Nr，此时就获得了一组(r，Nr)；
(4)计算剩余油占有网格数：改变网格大小，重新获取r，Nr，具体步骤重复(2)、(3)过程，就能得到一系列(r，Nr)；
(5)计盒维数的计算：对(ln1/r，lnNr)进行线性拟合，运用最小二乘法，其中直线的斜率值D就是该图像中剩余油的计盒维数；
步骤五：通过Photoshop软件对步骤四中已进行完网格化处理的剩余油图像九等分；其中，通过Photoshop软件添加图像长和宽的三等分点所在直线将图像分成九个区域；靠近注入介质端所在区域为注入端，靠近采出液体端所在区域为采出端，图像中心所在区域为注采中部，这三个区域为主流线区；所述主流线区的两侧区域为两翼区，其中，靠近注入端的区域为近井端，靠近采出端的区域为远井端，靠近注采中部的区域为盲端；在所述盲端与注采中部的边界上任意选取两处作为特征区，用于观察主流线区域是否波及到盲端；
步骤六，对步骤五中得到的主流线区域、两翼区域、主流线与两翼之间的区域以及特征区，计算这些区域的计盒维数并绘制计盒维数随PV数变化曲线；计算这些区域的计盒维数具体路径如下：
(1)设置像素点，图像大小为A×B个像素点；
(2)划分网格：用大小为n×n网格去覆盖剩余油全部图像；n是整数，则分割比率为r＝n/A，原图被分成了A/n个网格块，其中n序列为区域图像A和B的最大公约数的全部因数；
(3)被占用的判定规则：任何一个网格块只要包含剩余油就被认为是被占用的网格块，统计总共被占的网格块数并记为Nr，此时就获得了一组(r，Nr)；
(4)计算剩余油占有网格数：改变网格大小，重新获取r，Nr，具体步骤重复(2)、(3)过程，就能得到一系列(r，Nr)；
(5)计盒维数的计算：对(ln1/r，lnNr)进行线性拟合，运用最小二乘法，其中直线的斜率值D就是该图像的计盒维数；
步骤七，利用经过步骤六所获得的四组曲线分析基质化合物驱剩余油分布及运移规律，确定的具体路径如下：
(1)利用对应主流线区域的计盒维数随PV数变化曲线图按照如下规律确定基质化学驱条件下主流线剩余油分布及运移规律，设定△R＝R孔隙—R维数，其中R孔隙为孔隙的分形维数与x轴之间的长度；R维数为剩余油维数与x轴之间的长度；△R反映了区域的驱替程度；当△R介于1.1-1.3之间时，剩余油呈孤立状，零散的分布在孔喉中，说明驱替程度完全；当△R处于1-1.1时，剩余油呈条带状形态分布在少数连通的孔喉中，此时剩余油多数被驱替；
(2)利用对应两翼区域的计盒维数随PV数变化曲线图按照如下规律确定基质化学驱条件下两翼剩余油分布及运移规律，即：当△R处于1-1.1时，部分孔喉空间充满剩余油,呈斑块状形态分布，此时剩余油绝大多数被驱替；当△R处于0.8-1时，剩余油呈条带状形态分布在少数连通的孔喉中，此时剩余油多数被驱替；当△R处于1-1.1时，部分孔喉空间充满剩余油,呈斑块状形态分布；
(3)利用对应主流线与两翼之间的区域的计盒维数随PV数变化曲线图，定量分析主流线与两翼间剩余油分布及运移规律；所述定量分析路径如下：两翼较注采中部计盒维数变化幅度大且起始值高于注采中部，说明水驱结束以后，侧翼残余油高于注采中部，化合物有效进入侧翼驱替剩余油，两翼剩余油动用程度较注采中部高；
(4)利用对应特征区的计盒维数随PV数变化曲线图，按照如下规律确定基质化学驱条件下特征区剩余油分布及运移规律，即:当计盒维数减小时，说明水驱波及到侧翼。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810559938.5A CN108868715B (zh) | 2018-06-02 | 2018-06-02 | 运用分形理论确定基质-高渗条带化学驱剩余油分布及运移规律的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810559938.5A CN108868715B (zh) | 2018-06-02 | 2018-06-02 | 运用分形理论确定基质-高渗条带化学驱剩余油分布及运移规律的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108868715A CN108868715A (zh) | 2018-11-23 |
| CN108868715B true CN108868715B (zh) | 2020-05-01 |
Family
ID=64335770
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810559938.5A Expired - Fee Related CN108868715B (zh) | 2018-06-02 | 2018-06-02 | 运用分形理论确定基质-高渗条带化学驱剩余油分布及运移规律的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108868715B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112200689B (zh) * | 2019-07-08 | 2024-03-29 | 中国石油天然气股份有限公司 | 确定油藏渗流场潜力分散度的方法和装置 |
| CN112554853A (zh) * | 2020-12-11 | 2021-03-26 | 成都北方石油勘探开发技术有限公司 | 一种控制注水或注气的方法及系统 |
| CN113982573B (zh) * | 2021-11-22 | 2023-07-07 | 中国石油大学(北京) | 用于抑制驱替前缘指进的方法及处理器 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7575681B2 (en) * | 2004-07-06 | 2009-08-18 | Schlumberger Technology Corporation | Microfluidic separator |
| CN101424185B (zh) * | 2008-11-24 | 2012-09-05 | 大庆油田有限责任公司 | 注采主流线与分流线油藏剩余油参数特征及关系的地球化学测定方法 |
| CN102367731A (zh) * | 2011-09-23 | 2012-03-07 | 中国石油天然气股份有限公司 | 一种描述水驱后微观剩余油分布规律的方法 |
| CN102852516A (zh) * | 2012-04-19 | 2013-01-02 | 北京大学 | 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置 |
| CN103556994B (zh) * | 2013-11-19 | 2016-11-23 | 中国石油大学(华东) | 缝洞型储层剩余油分布的实验检测系统及检测方法 |
| CN104573339B (zh) * | 2014-12-24 | 2018-03-09 | 中国石油大学(北京) | 页岩气储层的地质参数确定方法和装置 |
| CN106437639A (zh) * | 2016-10-19 | 2017-02-22 | 中国石油化工股份有限公司 | 一种co2混相驱油效率和剩余油分布评价方法和专用装置 |
| CN107153074B (zh) * | 2017-03-21 | 2019-10-29 | 山东省科学院海洋仪器仪表研究所 | 一种基于高精度成像的微观剩余油定量评价方法 |
| CN206681717U (zh) * | 2017-05-07 | 2017-11-28 | 西南石油大学 | 一种模拟水和聚合物对剩余油在砂岩内分布影响的实验装置 |
-
2018
- 2018-06-02 CN CN201810559938.5A patent/CN108868715B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN108868715A (zh) | 2018-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108868715B (zh) | 运用分形理论确定基质-高渗条带化学驱剩余油分布及运移规律的方法 | |
| Zhang et al. | An integrated algorithm to evaluate flow direction and flow accumulation in flat regions of hydrologically corrected DEMs | |
| CN102913233B (zh) | 一种基于无因次对比图版的优势流动通道识别方法 | |
| CN108843310B (zh) | 运用分形理论确定基质-高渗条带水驱剩余油动用规律的方法 | |
| CN106703776A (zh) | 压裂参数优选方法 | |
| CN104331579A (zh) | 一种低渗透储层原油边界层的模拟方法 | |
| CN106600440B (zh) | 一种低渗透油藏调剖堵水动态指标选井的方法 | |
| CN106779211A (zh) | 一种用于油田注采井网驱替开发的产能预测方法 | |
| CN109142875B (zh) | 一种利用数字岩心获取致密砂岩岩石电学特性的方法 | |
| CN114091163B (zh) | 一种基于有限体积法的河道构筑物过流数值模拟方法 | |
| CN108843311B (zh) | 运用分形理论确定水驱剩余油分布及运移规律的方法 | |
| CN106530397B (zh) | 一种基于稀疏剖面地质轮廓线的地质面三维重建方法 | |
| CN104712328B (zh) | 快速评价复杂油藏中单个流动单元动用状况的方法 | |
| CN108952696B (zh) | 运用分形理论确定化学驱剩余油分布及运移规律的方法 | |
| CN112267873A (zh) | 一种模拟地层条件的单裂缝调驱可视化实验装置及方法 | |
| CN108678738B (zh) | 运用分形理论确定基质-高渗条带功能型聚合物驱剩余油分布及运移规律的方法 | |
| CN113158594B (zh) | 基于铸体薄片微观驱替模拟的驱油效率分析方法 | |
| CN108825222B (zh) | 运用分形理论确定功能型聚合物驱剩余油分布及运移规律的方法 | |
| CN107975366B (zh) | 一种堵水剂在多孔介质中封堵调剖可视化评价方法 | |
| CN202440786U (zh) | 河流入渗模拟装置 | |
| CN113643300A (zh) | 一种基于Seg-CapsNet算法的路面裂缝像素级别检测方法 | |
| CN110909466A (zh) | 提高差孔隙分选储层平均毛管压力曲线计算精度的方法 | |
| CN110286027A (zh) | 考虑红柳根系影响的河岸土体冲刷参数的量化方法 | |
| CN107808374A (zh) | 一种用于评价油气藏压裂效果的数值化方法 | |
| CN111091476B (zh) | 针对强非均质地层的油藏数值模拟方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200501 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |