CN108678738A - 运用分形理论确定基质-高渗条带功能型聚合物驱剩余油分布及运移规律的方法 - Google Patents
运用分形理论确定基质-高渗条带功能型聚合物驱剩余油分布及运移规律的方法 Download PDFInfo
- Publication number
- CN108678738A CN108678738A CN201810559947.4A CN201810559947A CN108678738A CN 108678738 A CN108678738 A CN 108678738A CN 201810559947 A CN201810559947 A CN 201810559947A CN 108678738 A CN108678738 A CN 108678738A
- Authority
- CN
- China
- Prior art keywords
- remaining oil
- image
- box
- counting dimension
- area
- 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.)
- Granted
Links
- 229920000642 polymer Polymers 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000009826 distribution Methods 0.000 title claims abstract description 38
- 238000013508 migration Methods 0.000 title claims abstract description 33
- 230000005012 migration Effects 0.000 title claims abstract description 33
- 230000008595 infiltration Effects 0.000 title claims abstract description 18
- 238000001764 infiltration Methods 0.000 title claims abstract description 18
- 230000008859 change Effects 0.000 claims abstract description 31
- 238000006073 displacement reaction Methods 0.000 claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims abstract description 20
- 238000002474 experimental method Methods 0.000 claims abstract description 8
- 238000002347 injection Methods 0.000 claims description 29
- 239000007924 injection Substances 0.000 claims description 29
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 15
- 210000004534 cecum Anatomy 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 12
- 238000004519 manufacturing process Methods 0.000 claims description 11
- 238000003709 image segmentation Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 8
- 230000000007 visual effect Effects 0.000 claims description 7
- 241001212149 Cathetus Species 0.000 claims description 6
- 238000004513 sizing Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 239000007788 liquid Substances 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000002360 preparation method Methods 0.000 claims description 3
- 238000004445 quantitative analysis Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 abstract description 9
- 239000011148 porous material Substances 0.000 abstract description 6
- 239000000284 extract Substances 0.000 abstract description 3
- 238000005192 partition Methods 0.000 abstract description 2
- 239000003921 oil Substances 0.000 description 111
- 239000011521 glass Substances 0.000 description 9
- 238000011084 recovery Methods 0.000 description 8
- 238000001259 photo etching Methods 0.000 description 7
- 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
- 230000007423 decrease Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 239000003292 glue Substances 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 229920006395 saturated elastomer Polymers 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000012153 distilled water Substances 0.000 description 2
- 235000019441 ethanol Nutrition 0.000 description 2
- 238000007687 exposure technique Methods 0.000 description 2
- 239000005357 flat glass Substances 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000032683 aging Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000004452 microanalysis Methods 0.000 description 1
- 229920002120 photoresistant polymer Polymers 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000035882 stress Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
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
- 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)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
一种运用分形理论确定基质‑高渗条带功能型聚合物驱剩余油分布及运移规律的方法。主要目的在于克服剩余油运移规律研究方法过于单一的缺点,本发明提出采用分形理论确定基质‑高渗条带功能型聚合物驱剩余油分布及运移规律的方法。具体解决方案为:制作微观物理模型,通过微观驱油装置进行基质‑高渗条带功能型聚合物驱驱油实验,收集渗流过程的连续图像,根据实际情况将驱替范围分区,通过图像处理,去除实际孔隙骨架,提取完整剩余油分布形态,按分区网格计算网格内剩余油计盒维数,总结相同网格及相邻网格计盒维数变化曲线,对比相邻网格计盒维数变化规律,确定基质‑高渗条带功能型聚合物驱剩余油分布及运移规律。
Description
技术领域
本发明隶属于油气田开发领域,尤其涉及到一种分形理论研究剩余油分布规律的装置与方法。
背景技术
随着油田开发的不断深入,国内大多数油田皆已进入高含水、高采出程度的“双高”阶段,针对二次采油未能采出的未波及区的剩余油和波及区的残余油,认识剩余油为油田二次采油及三次采油提供重要依据尤为重要。剩余油又分为宏观剩余油和微观剩余油两大类。目前微观方面关于剩余油的研究更侧重微观驱替机理及剩余油分布形态。更多的是采用微观可视模型,虽然直观,但是驱替过程中多区域同时流动,运移规律显得抽象,缺少对确定功能型聚合物驱剩余油分布及运移规律的方法。此外,剩余油研究方法较为单一,缺少用直观的规律曲线来分析剩余油分布及运移规律的方法。分形几何理论研究的是具有统计自相似性的对象,自诞生以来,在理论研究和实际应用方面都得到了快速发展,已有研究将分形理论用在孔隙结构、裂缝评价、储层非均质性描述,但关于剩余油方面目前只有计盒维数的计算公式,尚无具体运用分形理论研究剩余油分布及运移规律的方法。
发明内容
为了解决背景技术中所提到的技术问题,本发明提供一种运用分形理论确定基质-高渗条带功能型聚合物驱剩余油分布及运移规律的方法,该种方法将分形理论运用到剩余油分布研究领域,采用分形计盒维数定量化确定基质-高渗条带功能型聚合物驱剩余油分布及运移规律,从而克服现有技术中的缺陷。
本发明的技术方案是:该种运用分形理论确定化学驱剩余油分布及运移规律的方法,其特征在于包括如下步骤:
发明步骤一:布有基质-高渗条带的微观可视模型的制备
微观模型制备的基本原理是利用曝光技术将设计好的储层孔隙结构转移到涂有感光胶的玻璃片上并显影,然后利用氢氟酸腐蚀将显影部分刻蚀出具有一定深度的孔隙结构,而其他部分因分别被光刻胶和蜡保护而不被腐蚀,最后盖上盖玻片,这样便得到可视化驱替用的光刻玻璃岩心微观模型。根据上述原理,设制长4.5mm宽4.5mm的微观可视模型,并在对角线处布置基质高渗条带,其最大宽度不超过对角线长度的5分之1。基质高渗条带两侧的喉道半径为0.5-1.5μm,高渗条带处喉道半径为2-3μm。
发明步骤二:将步骤一中获得的微观物理模型置入微观驱油装置中进行微观驱油实验,采集功能型聚合物驱替过程图像
根据上述所制模型,实验步骤如下:
(1)用乙醇、石油醚、蒸馏水清洗光刻玻璃模型;
(2)将模型抽真空,饱和所模拟地层水;
(3)45℃条件下模型饱和油,并在此温度下静置24h;
(4)以0.01ml/min的驱替速度向模型中注入水直至含水率98%,后注入功能型聚合物至残余油饱和度;
(5)通过数字显微摄像系统采集驱替全程图像。
微量泵通过管线和注水活塞容器的下部阀门连接,活塞容器的上部阀门通过管线和光刻玻璃岩心模型入口的控制阀门连接,光刻玻璃岩心模型出口的控制阀门通过管线和量筒连接。
所述微量泵给整个驱替装置提供动力。
所述注水活塞容器为注入水的容器。
所述注功能型聚合物活塞容器为注入功能型聚合物的容器。
所述光刻玻璃为与实际储层接近的模型。
所述量筒承接采出液。
所述六通提供多个通路开关控制。
所述恒温箱使整个实验流程保持在地层温度下。
发明步骤三:将步骤二中采集到的驱替全程图像进行剩余油图像提取.
(1)将步骤二中采集到的驱替全程图像的RGB格式图像通过photoshop软件转化为灰度图像并对图像进行预处理,所述预处理的步骤包括进行亮度、对比度和锐化的处理以使得剩余油和孔隙能够得到清晰显示;
(2)采用迭代阀值法通过photoshop软件,对于预处理后的图像分割图像中的剩余油,提出剩余油分布完整图像;
进行图像分割时,最重要的部分就是如何选取阈值,如果阈值选择不合理,很容易将目标对象理解成为背景图像,或者将背景图像理解为目标对象。以步骤2中所获得的驱替图像中的一张图像,渗流图像为例,来说明具体路径。
具体路径如下:
①求出渗流图像中的最大灰度值Zmax和最小灰度值Zmin,并设置初始值,T0=(Zmin+Zmax)/2;
②根据阈值Tk(k=0,1,2,3…)将渗流图像分割成目标图像和背景图像两部分,再按照公式(1)和公式(2)求出这两部分的平均灰度值ZO和ZB;
(1)
(2)
式中,Z(i,j)为图像中像素点(i,j)的灰度值;N(I,j)为像素点(i,j)的加权系数,一般为1。
③求出新的阈值,Tk+1=(ZO+ZB)/2;
④如果Tk= Tk+1,则迭代结束,否则k=k+1,重复②—④。
上述的KTZ是通过photoshop软件按照公式算的,在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 数变化曲线图按照如下规律确定基质功能型聚合物驱条件下主流线剩余油分布及运移规律,设定△A=A1﹣A2,其中,A1为孔隙维数曲线与x坐标轴间的面积;A2 为化学驱剩余油计盒维数曲线与x坐标轴间的面积,△A反映了区域的驱替程度。
(2)利用对应两翼区域的计盒维数随 PV 数变化曲线图按照如下规律确定基质功能型聚合物驱条件下两翼剩余油分布及运移规律,
(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。基质高渗条带两侧的喉道半径为0.5-1.5μm,高渗条带处喉道半径为2-3μm。
步骤二:进行微观驱替实验,采集功能型聚合物驱替过程图像。
实验步骤如下:
(1)用石油醚、乙醇、蒸馏水清洗微观仿真玻璃刻蚀模型;
(2)将清洗干净后的玻璃可视模型抽真空后,饱和模拟地层水;
(3)20℃条件下在模型中饱和油,并在此温度条件下静置老化24h;
(4)以0.01ml/min的驱替速度向模型中注入水直至含水率98%,后注入功能型聚合物至残余油饱和度;
(5)分析微观驱替图像。
发明步骤三:将步骤二中采集到的驱替全程图像进行剩余油图像提取。
(1)将步骤二中采集到的驱替全程图像的RGB格式图像通过photoshop软件转化为灰度图像并对图像进行预处理,所述预处理的步骤包括进行亮度、对比度和锐化的处理以使得剩余油和孔隙能够得到清晰显示;
(2)采用迭代阀值法通过photoshop软件,对于预处理后的图像分割图像中的剩余油,提出剩余油分布完整图像;
进行图像分割时,最重要的部分就是如何选取阈值,如果阈值选择不合理,很容易将目标对象理解成为背景图像,或者将背景图像理解为目标对象。以步骤2中所获得的驱替图像中的一张图像,渗流图像为例,来说明具体路径。
具体路径如下:
①求出渗流图像中的最大灰度值Zmax和最小灰度值Zmin,并设置初始值,T0=(Zmin+Zmax)/2;
②根据阈值Tk(k=0,1,2,3…)将渗流图像分割成目标图像和背景图像两部分,再按照公式(1)和公式(2)求出这两部分的平均灰度值ZO和ZB;如图4所示。
迭代阈值法的实现步骤如下:
①求出渗流图像的最大灰度值Zmax和最小灰度值Zmin,并设置初始值,T0=(Zmin+Zmax)/2;
②根据阈值Tk(k=0,1,2,3…)将渗流图像分割成目标图像和背景图像两部分,再按照公式(1)和公式(2)求出这两部分的平均灰度值ZO和ZB;
(1)
(2)
式中,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所示,图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所示,图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 数变化曲线图按照如下规律确定基质功能型聚合物驱条件下主流线剩余油分布及运移规律,设定△A=A1﹣A2 ,其中,A1为孔隙维数曲线与x坐标轴间的面积;A2 为化学驱剩余油计盒维数曲线与x坐标轴间的面积,△A反映了区域的驱替程度。如图6所示为主流线上区域计盒维数的变化曲线,△A1>△A2 >△A3 ,说明注入端采出程度最大,注采中部较小,采出端采出程度最小。随着功能型聚合物的进行,注入端计盒维数在注入量0至0.1PV区间略微下降,说明注入端在水驱过程中绝大部分剩余油已被驱替,注入端少量的剩余油受功能型聚合物驱替导致剩余油不断运移,注入端剩余油不断减少,在注入量为0.1PV之后,计盒维数趋于平缓,说明注入端已无功能型聚合物可动剩余油,功能型聚合物已无动用潜力。而注采中部在注入端维数下降中计盒维数上下波动,说明注采中部的剩余油受到注入聚合物的波及以及注入端剩余油的运移,计盒维数在注入量为0.1PV至0.25PV区间快速下降,说明功能型聚合物已进入注采中部驱替剩余油,注入量在0.25PV后,计盒维数趋于平缓,说明注采中部已无可动剩余油,聚驱再无动用潜力。采出端计盒维数在0.2PV前上下波动,说明采出端受功能型聚合物的波及以及注采中部剩余油的运移,在0.2PV至0.3PV维数大幅下降,说明功能型聚合物进入并驱替了采出端的剩余油,在注入量为0.3PV后1、2、3区域计盒维数趋于平缓,说明功能型聚合物已无动用潜力。
(2)利用对应两翼区域的计盒维数随 PV 数变化曲线图按照如下规律确定基质功能型聚合物驱条件下两翼剩余油分布及运移规律,如图7所示为侧翼区域计盒维数的变化曲线,△A4>△A6 >△A5 ,说明近井端采出程度最大,远井端较小,盲端端采出程度最小。随着聚驱的进行,近井端在注入量为0至0.2PV维数大幅下降,说明功能型聚合物进入近井端驱替剩余油导致剩余油发生运移,在注入量为0.2PV后,计盒维数趋于平缓,说明近井端已无聚驱可动剩余油。盲端的计盒维数在注入量为0至0.25PV区间上下波动,说明盲端受聚合物的波及以及近井端剩余油的运移,在注入量0.25PV至0.3PV区间计盒维数大幅下降,说明聚合物已注入到盲端,驱替剩余油,后在0.3PV时盲端已无聚驱可动剩余油。远井端在注入量为0至0.3PV区间时,计盒维数上下波动,说明受聚合物波及以及盲端剩余油的运移,在注入量为0.3PV至0.35PV区间计盒维数快速下降,说明剩余油被大幅动用,在0.35PV时各区计盒维数趋于平缓,说明聚驱再无动用潜力。
(3)利用对应主流线与两翼之间的区域的计盒维数随PV数变化曲线图,定量分析主流线与两翼间剩余油分布及运移规律;如图8所示,两翼较注采中部计盒维数变化幅度大且起始值高于注采中部。说明水驱结束以后,侧翼残余油高于注采中部,功能型聚合物有效进入侧翼驱替剩余油,两翼剩余油动用程度较注采中部高。
(4)利用对应特征区的计盒维数随 PV 数变化曲线图,按照如下规律确定基质功能型聚合物驱条件下特征区剩余油分布及运移规律,即:当计盒维数减小时,说明水驱波及到侧翼。
如图9所示特征区10在注入量0至0.1PV区间的计盒维数为0,说明主流线区域没有剩余油运移至特征区10,而在注入量为0.1PV至0.2PV区间和0.2PV至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 数变化曲线图按照如下规律确定基质功能型聚合物驱条件下主流线剩余油分布及运移规律:设定△A=A1﹣A2 ,其中,A1为孔隙维数曲线与x坐标轴间的面积;A2 为化学驱剩余油计盒维数曲线与x坐标轴间的面积,△A反映了区域的驱替程度;
(2)利用对应两翼区域的计盒维数随 PV 数变化曲线图按照如下规律确定基质功能型聚合物驱条件下两翼剩余油分布及运移规律;
(3)利用对应主流线与两翼之间的区域的计盒维数随PV数变化曲线图,定量分析主流线与两翼间剩余油分布及运移规律;
(4)利用对应特征区的计盒维数随 PV 数变化曲线图,按照如下规律确定基质功能型聚合物驱条件下特征区剩余油分布及运移规律,即:当计盒维数减小时,说明水驱波及到侧翼。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810559947.4A CN108678738B (zh) | 2018-06-02 | 2018-06-02 | 运用分形理论确定基质-高渗条带功能型聚合物驱剩余油分布及运移规律的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810559947.4A CN108678738B (zh) | 2018-06-02 | 2018-06-02 | 运用分形理论确定基质-高渗条带功能型聚合物驱剩余油分布及运移规律的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108678738A true CN108678738A (zh) | 2018-10-19 |
CN108678738B CN108678738B (zh) | 2021-11-12 |
Family
ID=63809653
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810559947.4A Expired - Fee Related CN108678738B (zh) | 2018-06-02 | 2018-06-02 | 运用分形理论确定基质-高渗条带功能型聚合物驱剩余油分布及运移规律的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108678738B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004095259A1 (en) * | 2003-03-26 | 2004-11-04 | Exxonmobil Upstream Research Company | Performance prediction method for hydrocarbon recovery processes |
US20070251389A1 (en) * | 2006-04-26 | 2007-11-01 | Dina Katsir | Composite inorganic membrane for separation in fluid systems |
CN102339339A (zh) * | 2010-07-23 | 2012-02-01 | 中国石油化工股份有限公司 | 一种分析缝洞型油藏剩余油分布的方法 |
CN103339488A (zh) * | 2011-01-27 | 2013-10-02 | 普拉德研究及开发股份有限公司 | 非常规岩石样品的气体吸附分析 |
CN104076046A (zh) * | 2013-03-28 | 2014-10-01 | 中国石油化工股份有限公司 | 多孔介质中剩余油微观分布图像采集与定量表征方法 |
CN105551004A (zh) * | 2015-12-29 | 2016-05-04 | 中国石油大学(华东) | 一种基于岩心ct图像处理的剩余油微观赋存表示方法 |
CN107153074A (zh) * | 2017-03-21 | 2017-09-12 | 山东省科学院海洋仪器仪表研究所 | 一种基于高精度成像的微观剩余油定量评价方法 |
CN108062789A (zh) * | 2017-12-20 | 2018-05-22 | 中国石油天然气股份有限公司 | 岩心样品选取方法和装置 |
-
2018
- 2018-06-02 CN CN201810559947.4A patent/CN108678738B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004095259A1 (en) * | 2003-03-26 | 2004-11-04 | Exxonmobil Upstream Research Company | Performance prediction method for hydrocarbon recovery processes |
US20060224369A1 (en) * | 2003-03-26 | 2006-10-05 | Yang Shan H | Performance prediction method for hydrocarbon recovery processes |
US20070251389A1 (en) * | 2006-04-26 | 2007-11-01 | Dina Katsir | Composite inorganic membrane for separation in fluid systems |
CN102339339A (zh) * | 2010-07-23 | 2012-02-01 | 中国石油化工股份有限公司 | 一种分析缝洞型油藏剩余油分布的方法 |
CN103339488A (zh) * | 2011-01-27 | 2013-10-02 | 普拉德研究及开发股份有限公司 | 非常规岩石样品的气体吸附分析 |
CN104076046A (zh) * | 2013-03-28 | 2014-10-01 | 中国石油化工股份有限公司 | 多孔介质中剩余油微观分布图像采集与定量表征方法 |
CN105551004A (zh) * | 2015-12-29 | 2016-05-04 | 中国石油大学(华东) | 一种基于岩心ct图像处理的剩余油微观赋存表示方法 |
CN107153074A (zh) * | 2017-03-21 | 2017-09-12 | 山东省科学院海洋仪器仪表研究所 | 一种基于高精度成像的微观剩余油定量评价方法 |
CN108062789A (zh) * | 2017-12-20 | 2018-05-22 | 中国石油天然气股份有限公司 | 岩心样品选取方法和装置 |
Non-Patent Citations (5)
Title |
---|
AKHLAGHI AMIRI,H.A.等: "Pore-scale modeling of non-isothermal two phase flow in 2D porous media: Influences of viscosity, capillarity, wettability and heterogeneity", 《INTERNATIONAL JOURNAL OF MULTIPHASE FLOW》 * |
MOGENSEN,K.等: "A Dynamic Two-Phase Pore-Scale Model of Imbibition", 《TRANSPORT IN POROUS MEDIA》 * |
孙焕泉等: "胜坨油田特高含水期剩余油分布仿真模型", 《石油勘探与开发》 * |
朱九成等: "指进、剩余油形成与分布的物理模拟", 《新疆石油地质》 * |
李中锋等: "微观物理模拟水驱油实验及残余油分布分形特征研究", 《中国石油大学学报( 自然科学版)》 * |
Also Published As
Publication number | Publication date |
---|---|
CN108678738B (zh) | 2021-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105041307A (zh) | 一种碎屑岩油气储层优势渗流通道识别的工业化流程 | |
CN108868715A (zh) | 运用分形理论确定基质-高渗条带化学驱剩余油分布及运移规律的方法 | |
CN106651158B (zh) | 一种超低渗致密油藏水平井注水开发见效程度定量评价方法 | |
CN111810101B (zh) | 一种水驱油藏动态分析方法及装置 | |
CN109209312A (zh) | 一种适用聚合物驱的资源开发方法 | |
CN108843310A (zh) | 运用分形理论确定基质-高渗条带水驱剩余油动用规律的方法 | |
CN109522619B (zh) | 一种低渗油田周期采油方法 | |
CN108678738A (zh) | 运用分形理论确定基质-高渗条带功能型聚合物驱剩余油分布及运移规律的方法 | |
CN108843311A (zh) | 运用分形理论确定水驱剩余油分布及运移规律的方法 | |
CN107288595B (zh) | 一种注水利用率的评价方法 | |
CN111350498B (zh) | 一种中高渗油藏特高含水开发期弱驱分布特征描述的方法 | |
CN108952696A (zh) | 运用分形理论确定化学驱剩余油分布及运移规律的方法 | |
CN110909477B (zh) | 一种酸蚀裂缝粗糙程度定量表征方法 | |
CN109344508B (zh) | 一种缝洞型碳酸盐岩油藏模型建立方法 | |
CN107169227A (zh) | 一种分段压裂水平井的粗网格模拟方法及系统 | |
CN106703765A (zh) | 一种调整注采液位的可视化人工填砂驱油装置及其应用 | |
CN108825222A (zh) | 运用分形理论确定功能型聚合物驱剩余油分布及运移规律的方法 | |
CN107227946A (zh) | 驱替实验的实验数据获取方法及系统 | |
CN112502677B (zh) | 一种基于多元线性回归的注水开发效果评价方法 | |
RU2290501C1 (ru) | Способ разработки нефтяной залежи | |
CN112392456B (zh) | 一种井网布局合理性判断方法及井网布局优化方法 | |
CN107808374A (zh) | 一种用于评价油气藏压裂效果的数值化方法 | |
CN110069850A (zh) | 一种基于bim创建堤坝横断面的方法 | |
CN1517510A (zh) | 油田注水井分层流量调配方法 | |
CN110778314A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20211112 |