CN105676309A - 裂缝性油藏逆向渗吸采收率预测方法 - Google Patents
裂缝性油藏逆向渗吸采收率预测方法 Download PDFInfo
- Publication number
- CN105676309A CN105676309A CN201610060046.1A CN201610060046A CN105676309A CN 105676309 A CN105676309 A CN 105676309A CN 201610060046 A CN201610060046 A CN 201610060046A CN 105676309 A CN105676309 A CN 105676309A
- Authority
- CN
- China
- Prior art keywords
- crack
- matrix
- finite elements
- imbibition
- fracture
- 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
- 238000005213 imbibition Methods 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000011084 recovery Methods 0.000 title claims abstract description 46
- 239000011159 matrix material Substances 0.000 claims abstract description 116
- 238000010168 coupling process Methods 0.000 claims abstract description 22
- 230000008569 process Effects 0.000 claims abstract description 20
- 230000008859 change Effects 0.000 claims abstract description 18
- 238000004088 simulation Methods 0.000 claims abstract description 18
- 230000008878 coupling Effects 0.000 claims abstract description 13
- 238000005859 coupling reaction Methods 0.000 claims abstract description 13
- 238000005520 cutting process Methods 0.000 claims abstract description 5
- 206010017076 Fracture Diseases 0.000 claims description 58
- 208000010392 Bone Fractures Diseases 0.000 claims description 48
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 37
- 239000011435 rock Substances 0.000 claims description 29
- 239000012530 fluid Substances 0.000 claims description 25
- 238000011161 development Methods 0.000 claims description 21
- 230000035699 permeability Effects 0.000 claims description 16
- 239000008346 aqueous phase Substances 0.000 claims description 12
- 239000012071 phase Substances 0.000 claims description 12
- 230000000694 effects Effects 0.000 claims description 10
- 238000009792 diffusion process Methods 0.000 claims description 6
- 238000003384 imaging method Methods 0.000 claims description 6
- 239000011148 porous material Substances 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims description 4
- 238000005481 NMR spectroscopy Methods 0.000 claims description 3
- 241001282153 Scopelogadus mizolepis Species 0.000 claims description 3
- 238000002591 computed tomography Methods 0.000 claims description 3
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 3
- 229910052753 mercury Inorganic materials 0.000 claims description 3
- 230000007547 defect Effects 0.000 abstract 1
- 239000003921 oil Substances 0.000 description 11
- 238000002474 experimental method Methods 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 4
- 230000001788 irregular Effects 0.000 description 4
- 239000010779 crude oil Substances 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000009533 lab test Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 210000004534 cecum Anatomy 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明涉及一种裂缝性油藏逆向渗吸采收率预测方法,该方法考虑孤立裂缝、基质岩块切割的不规则性以及裂缝与基质之间接触的多样性,建立裂缝性介质模型,并基于离散裂缝模型的数值模拟方法及裂缝与基质之间的耦合,数值模拟裂缝性介质模型的渗吸过程与规律,预测采收率变化。本发明全面模拟裂缝性油藏的渗吸过程和渗吸规律,研究预测渗吸作用下裂缝性介质内剩余油变化,实现了对裂缝的精确刻画,能够模拟出孤立裂缝的影响,预测结果准确,弥补常规数值模拟方法的不足。
Description
技术领域
本发明属于石油开采技术领域,具体地说,涉及一种预测裂缝性油藏逆向渗吸采收率的方法,特别是针对存在孤立裂缝的基质岩块。
背景技术
裂缝性油藏中原油的采收率低、见水早,基质内原油的驱替难度大。对于裂缝性油藏,特别是低渗透油藏,通常是见高含水后,采取闷井停产一段时间后再生产,依靠毛管力作用,即渗吸作用来驱替基质内部的原油。闷井时间的长短、渗吸作用下的采收率变化等方面的研究对裂缝性油藏提高开采率具有重要的意义。
现有的基质内渗吸采收率预测方法主要有三种:室内实验、解析解分析和数值解分析。大多集中在室内小尺度岩样的实验分析和确定采收率与无因次时间tD之间的变化曲线上(Mattax、Kyte(1962);Bourbiaux、Kalaydjian(1990);Zhang(1996))。Zhang等人(1996)的实验表明:不同的边界条件下得到的采收率与tD之间的变化曲线,均逼近统一表达式。Cilet等人(1998)实验分析了逆向渗吸过程中的影响因素,提出了预测采收率的关系式其中R为采收率,R∞为最终采收率,常数ωD和n取0.0033和0.6,tD为无量纲时间。还有研究者通过数学计算的方法来研究渗吸过程中的采收率,给出了逆向渗吸过程中预测采收率的解析表达式 其中为无量纲时间(Zimmerman和Bodvasson(1989);Zimmerman(1990);Chen(1995);Tavassoil(2005);Li(2011))。Ma等人(1997)通过不同亲水类型的岩石实验结果分析,认为不同亲水类型的岩心,其采收率与无因次时间tD的变化曲线相似,最后均可通过采收率与无因次时间tD来表示,并给出关系式其中,R为采收率,R∞为最终采收率,α为常数,当取α=0.05时适用于大多数实验结果。Ma等人给出的关系式适用于强亲水介质,且忽略了润湿性的影响。Gupta和Civan(1994)和Cil(1998)引入了油水接触角θ的余弦项,代表介质的润湿性影响。当介质亲水时,则cosθ=1,当忽略润湿性的影响,则α是接触角分布的函数。Zhou(2002)通过高孔隙度低渗透率的岩心的物理实验,研究了较大流度比范围的渗吸过程,并给出了新的表达式其中tD为无量纲时间,t为时间,k为渗透率,φ为孔隙度,σ为表面张力,θ为接触角,μnw为流体粘度,Lc特征长度。采用数值模拟的方法来研究渗吸过程(Pooladi-Darvish、Firoozabadi(2000);Behbahani(2006);Bagherinezhad、Pishvaie(2014)),大多是建立在连续介质模型基础上,Behbahani(2006)应用模拟软件,采用细网格模拟渗吸过程。上述三种方法各有优势和不足:
(1)室内物理实验是地下基质岩心在室内条件下物理过程的真实再现。但室内岩心大多是圆柱形或者长方体等规则形状;鉴于物理过程的缓慢性和实验条件的有限性,大多实验岩样尺寸为小岩心岩样;在室内实验过程中,借助于测量仪器仪表实现实验数据测量,一方面测量数据具有局限性,另一方面测量数据容易受到人为干扰和存在测量误差;物理实验岩心的加工制作时间和物理实验过程耗时比较长,需要实验周期和人力时间较长。
(2)渗吸作用的解析解是对渗吸过程采用数学方程进行描述,求解控制方程的解析解。但要求岩心模型为均匀各向同性,难以适用于非均质和各向异性情况,对于岩样内部发育裂缝的情况更无法获得其解析解;方程是建立在圆柱坐标或者笛卡尔坐标下,解析解可用于分析圆柱或者长方体等规则形状的渗吸规律,难以适用于不规则形状的岩样模型。
(3)数值模拟是采用数值法求解控制方程来预测渗吸过程。数值模拟可以适用于不规则形状的岩样模型;数值模拟可以方便的模拟不同边界条件、不同岩样模型大小、不同时刻的渗吸状况,可以获得模型中任意网格点处的数值,具有省时边界、获得数据量丰富的优点。当前数值模拟主要基于连续介质模型。
目前,在逆向渗吸研究中,研究者认为,无论哪一种方法获得的采收率R和无因次时间tD之间的变化关系均可表示成:
其中,R为采收率,R∞为最终采收率,α为常数,当取α=0.05时适用于大多数实验结果。
对于裂缝性油藏而言,不仅存在连通性较好的裂缝网络,还存在裂缝网络连通性差的孤立裂缝或裂缝盲端,由于具有一定的开度,其毛管力要小于基质孔隙的毛管力,使得基质岩块表相处强裂缝具有非均质性,进而影响了基岩逆向渗吸采收率变化。因此,基质岩块内有无裂缝、裂缝的多少和形态对整个岩块的采收率变化影响不同,目前还未见如何预测发育孤立裂缝的基质内的渗吸采收率的相关报道。
发明内容
本发明针对现有裂缝性油藏预测逆向渗吸采收率时存在的结果不准确等上述不足,提供一种裂缝性油藏逆向渗吸采收率预测方法,该方法考虑孤立裂缝、基质岩块切割的不规则性以及裂缝与基质之间接触的多样性,全面模拟裂缝性油藏的渗吸过程和渗吸规律,研究预测渗吸作用下裂缝性介质内剩余油变化,实现了对裂缝的精确刻画,能够模拟出孤立裂缝的影响,预测结果准确。
本发明的技术方案是:一种裂缝性油藏逆向渗吸采收率预测方法,含有以下步骤:
(一)考虑裂缝对基质岩块切割的不规则性、裂缝与基质之间接触的多样性以及基质中存在的孤立裂缝,建立裂缝性介质模型,该裂缝性介质模型为基质岩块切割后任一岩块边缘裂缝首尾相接组成;
(二)采用非结构化网格对裂缝性介质模型进行离散化处理,获得离散裂缝模型;
(三)根据离散裂缝模型确定基质的渗吸方程和裂缝的渗吸方程,并根据含水饱和度和流量关系进行基质有限单元之间、裂缝有限单元之间、基质有限单元与裂缝有限单元之间的耦合,耦合后获得整体的有限元方程;
(四)数值模拟裂缝性介质模型的渗吸过程与规律,预测渗吸作用下饱和度变化和采收率变化。
在根据本发明上述预测方法中,步骤(一)中,建立裂缝性介质模型的步骤为:
(一)搜集反应裂缝发育状况的岩心、地震、成像测井的资料;
(二)根据岩心统计结果、地震预测结果、成像测井检测结果,统计裂缝发育状况、裂缝发育密度、裂缝长度、裂缝开度、裂缝倾角的信息,获得大尺度裂缝网络的连通性信息、裂缝切割基质岩块形状信息以及裂缝与基质的接触关系信息;
(三)通过CT扫描、扫描电镜、高压压汞和核磁共振分析小尺度裂缝的连通性,确定基质内部连通性差、或孤立裂缝发育状况、或连通性差和孤立裂缝发育状况;
(四)根据裂缝网络对基质岩块的切割尺寸和接触关系,以及连通性差、或孤立裂缝发育状况、或连通性差和孤立裂缝发育状况,确定裂缝性介质模型。
在根据本发明上述预测方法中,步骤(三)中,基质的渗吸方程表达式为:
式中,vm,w为基质内水相流体流速,Pm,c为基质中毛管力,Km为基质的绝对渗透率,km,rw为基质中水相流体相对渗透率,km,ro为基质中油相流体相对渗透率,μw为水相流体年度,μo为油相流体粘度;
进一步表示为:
其中,Sm,w为基岩含水饱和度,基岩内分流量fm,w为:
定义基岩内毛管力扩散系数Dm,e为:
则有:
同时满足质量守恒,则有:
采用混合有限元法求解基质的渗吸方程式(6)和(7)获得基质有限单元上的渗吸方程表达式为:
式中,SmwA,e为基质有限单元e上面含水饱和度,为矩阵Ae的逆矩阵,为基质有限单元e在边li上边含水饱和度,i=1,2,3,φm为基质孔隙度。
裂缝的渗吸方程表达式为:
式中,vf,w为裂缝内水相流体流速,Pf,c为裂缝中毛管力,Kf为裂缝的绝对渗透率,kf,ro为裂缝中油相流体相对渗透率,μo为油相流体粘度,Sf,w为裂缝含水饱和度,裂缝内分流量ff,w为:
定义裂缝内毛管力扩散系数Df,e为:
其数值的大小与裂缝开度有关,当裂缝开度大于0.1mm,因裂缝内部毛管力较小,Df,e数值很小,当裂缝开度小于等于0.1mm时,其内部毛管力不可忽略。则有:
同时满足质量守恒,则有:
采用混合有限元法求解裂缝的渗吸方程式(13)和(14)获得裂缝有限单元上的渗吸方程表达式为:
式中,Sfwl,l为裂缝有限单元l上边含水饱和度,为矩阵Ll的逆矩阵,为裂缝有限单元l在端点di上含水饱和度,i=1,2。
在根据本发明上述预测方法中,步骤(三)中,基质有限单元之间的耦合过程为:定义在两个相邻有限单元之间的边界上,通过边界流出第一有限单元的流量等同于通过边界流入第二有限单元的流量,第一有限单元在边界上的含水饱和度等于第二有限单元在边界上的含水饱和度,如此完成基质有限单元之间的耦合;基质有限单元与裂缝有限单元之间的耦合过程为:定义在裂缝有限单元与相邻的两个基质有限单元之间,流出裂缝有限单元的流量等于通过边界流入到两个基质有限单元的流量之和,如此完成基质有限单元与裂缝有限单元之间的耦合;裂缝有限单元之间的耦合过程为:定义在裂缝有限单元的节点处,经过该节点流出裂缝有限单元的流量等于经过该节点流入到裂缝有限单元的流量,在裂缝有限单元节点处,各个裂缝有限单元在节点上的含水饱和度相等,如此完成裂缝有限单元之间的耦合。
在根据本发明上述预测方法中,步骤(三)中,耦合后获得的整个有限元方程为:
式中,m为基质,f为裂缝,Dm、Bm,m、Bm,f、Bf,f、Mm,m、Mm,f、Mf,f、Df、Mf均为系数矩阵。
本发明的有益效果是:本发明考虑了孤立裂缝与基质孔隙之间的差异,以及裂缝对基质岩块切割的不规则性和裂缝与基质之间接触的多样性,建立裂缝性介质模型,并基于离散裂缝模型全面模拟裂缝性油藏的渗吸过程和渗吸规律,对裂缝进行精确的刻画。与现有预测方法相比,本发明可适用于被裂缝切割后形成的不规则形状的基质岩块;本发明不仅可以考虑基质的非均质性,还可以考虑基质岩块内孤立裂缝与基质之间毛管力大小的差异,以及孤立裂缝的不同发育形态;本发明可利用CT成像和扫描电镜对基质内小尺度裂缝的识别,通过非结构化网格的剖分来刻画小尺度裂缝的形态。通过本发明预测方法预测的渗吸采收率更加准确,本发明属于数值模拟,能很好的模拟出孤立裂缝的影响,快速精确,弥补常规数值模拟方法的不足。
附图说明
图1为本发明实施例1中裂缝性介质模型。
图2为本发明实施例裂缝性介质模型非结构化网格剖分图。
图3为本发明实施例1中离散裂缝模型。
图4为本发明实施例2中裂缝性介质模型。
图5为本发明实施例2中离散裂缝模型。
图6为本发明渗吸采收率预测方法与现有渗吸采收率方法的数值模拟数据与实验结果对比图。
图7为本发明具体实施例含裂缝的裂缝性介质模型示意图。
图8为本发明具体实施例无因次时间tD=1.03时的无裂缝和含裂缝的含水饱和度分布图。
图9为本发明具体实施例无因次时间tD=11.1时的无裂缝和含裂缝的含水饱和度分布图。
图10为本发明具体实施例无因次时间tD=708时的无裂缝和含裂缝的含水饱和度分布图。
图11为本发明具体实施例不同长度裂缝及对应于裂缝性介质模型中的位置图。
图12为本发明具体实施例不同长度裂缝对应的采收率与无因次时间的变化曲线。
具体实施方式
以下结合附图和具体实施例对本发明作出进一步说明。
实施例1:一种裂缝性油藏逆向渗吸采收率预测方法,含有以下步骤:
步骤一:考虑裂缝对基质岩块切割的不规则性、裂缝与基质之间接触的多样性以及基质中存在的孤立裂缝,建立裂缝性介质模型,如图1所示,该裂缝性介质模型为基质岩块切割后不规则五边形岩块边缘裂缝首尾相接组成。其具体步骤为:
(一)搜集反应裂缝发育状况的岩心、地震、成像测井的资料;
(二)根据岩心统计结果、地震预测结果、成像测井检测结果,统计裂缝发育状况、裂缝发育密度、裂缝长度、裂缝开度、裂缝倾角的信息,获得大尺度裂缝网络的连通性信息、裂缝切割基质岩块形状信息以及裂缝与基质的接触关系信息;
(三)通过CT扫描、扫描电镜、高压压汞和核磁共振分析小尺度裂缝的连通性,确定基质内部连通性差、或孤立裂缝发育状况、或连通性差和孤立裂缝发育状况;
(四)根据裂缝网络对基质岩块的切割尺寸和接触关系,以及连通性差、或孤立裂缝发育状况、或连通性差和孤立裂缝发育状况,确定裂缝性介质模型。
步骤二:如图2所示,采用非结构化网格对裂缝性介质模型进行离散化处理,获得离散裂缝模型,该离散裂缝模型如图3所示。
步骤三:根据离散裂缝模型确定基质的渗吸方程和裂缝的渗吸方程,并根据含水饱和度和流量关系进行基质有限单元之间、裂缝有限单元之间、基质有限单元与裂缝有限单元之间的耦合,耦合后获得整体的有限元方程。
基质的渗吸方程表达式为:
式中,vm,w为基质内水相流体流速,Pm,c为基质中毛管力,Km为基质的绝对渗透率,km,rw为基质中水相流体相对渗透率,km,ro为基质中油相流体相对渗透率,μw为水相流体年度,μo为油相流体粘度;
进一步表示为:
其中,Sm,w为基岩含水饱和度,基岩内分流量fm,w为:
定义基岩内毛管力扩散系数Dm,e为:
则有:
同时满足质量守恒,则有:
采用混合有限元法求解基质的渗吸方程式(6)和(7)获得基质有限单元上的渗吸方程表达式为:
式中,SmwA,e为基质有限单元e上面含水饱和度,为矩阵Ae的逆矩阵,为基质有限单元e在边li上边含水饱和度,i=1,2,3,φm为基质孔隙度。
裂缝的渗吸方程表达式为:
式中,vf,w为裂缝内水相流体流速,Pf,c为裂缝中毛管力,Kf为裂缝的绝对渗透率,kf,ro为裂缝中油相流体相对渗透率,μo为油相流体粘度,Sf,w为裂缝含水饱和度,裂缝内分流量ff,w为:
定义裂缝内毛管力扩散系数Df,e为:
其数值的大小与裂缝开度有关,当裂缝开度大于0.1mm,因裂缝内部毛管力较小,Df,e数值很小,当裂缝开度小于等于0.1mm时,其内部毛管力不可忽略。则有:
同时满足质量守恒,则有:
采用混合有限元法求解裂缝的渗吸方程式(13)和(14)获得裂缝有限单元上的渗吸方程表达式为:
式中,Sfwl,l为裂缝有限单元l上边含水饱和度,为矩阵Ll的逆矩阵,为裂缝有限单元l在端点di上含水饱和度,i=1,2。
基质有限单元之间的耦合过程为:定义在两个相邻有限单元之间的边界上,通过边界流出第一有限单元的流量等同于通过边界流入第二有限单元的流量,第一有限单元在边界上的含水饱和度等于第二有限单元在边界上的含水饱和度,如此完成基质有限单元之间的耦合。
基质有限单元与裂缝有限单元之间的耦合过程为:定义在裂缝有限单元与相邻的两个基质有限单元之间,流出裂缝有限单元的流量等于通过边界流入到两个基质有限单元的流量之和,如此完成基质有限单元与裂缝有限单元之间的耦合。
裂缝有限单元之间的耦合过程为:定义在裂缝有限单元的节点处,经过该节点流出裂缝有限单元的流量等于经过该节点流入到裂缝有限单元的流量,在裂缝有限单元节点处,各个裂缝有限单元在节点上的含水饱和度相等,根据完成裂缝有限单元之间的耦合。
耦合后获得的整个有限元方程为:
式中,m为基质,f为裂缝,Dm、Bm,m、Bm,f、Bf,f、Mm,m、Mm,f、Mf,f、Df、Mf均为系数矩阵。
步骤四:数值模拟裂缝性介质模型的渗吸过程与规律,预测渗吸作用下饱和度变化和采收率变化。
实施例二:与实施例一不同的是,实施例二中提供的一种裂缝性油藏逆向渗吸采收率预测方法,其步骤一中建立裂缝性介质模型如图4所示,该裂缝性介质模型为基质岩块切割后长方形岩块边缘裂缝首尾相接组成;步骤二中获得的离散模型如图5。其余步骤同实施例一。
本发明渗吸采收率预测方法为数值模拟基质内渗吸过程,预测渗吸作用下饱和度变化和采收率变化。
由图6可知,本发明数值模拟数据的实验结果与现有渗吸采收率方法的数值模拟数据的实验结果相对比,本发明的准确性和有效性更好。
由图8至10可知,当裂缝性介质模型中含有裂缝时,其含水饱和度变化与无裂缝时存在明显不同。
由图12可知,裂缝性介质模型中发育裂缝的长度和位置直接影响其内部渗吸作用下含水率变化及采收率变化。
以上所举实施例仅用为方便举例说明本发明,并非对本发明保护范围的限制,在本发明所述技术方案范畴,所属技术领域的技术人员所作各种简单变形与修饰,均应包含在以上申请专利范围中。
Claims (5)
1.一种裂缝性油藏逆向渗吸采收率预测方法,其特征在于:含有以下步骤:
(一)考虑裂缝对基质岩块切割的不规则性、裂缝与基质之间接触的多样性以及基质中存在的孤立裂缝,建立裂缝性介质模型,该裂缝性介质模型为基质岩块切割后任一岩块边缘裂缝首尾相接组成;
(二)采用非结构化网格对裂缝性介质模型进行离散化处理,获得离散裂缝模型;
(三)根据离散裂缝模型确定基质的渗吸方程和裂缝的渗吸方程,并根据含水饱和度和流量关系进行基质有限单元之间、裂缝有限单元之间、基质有限单元与裂缝有限单元之间的耦合,耦合后获得整体的有限元方程。
(四)数值模拟裂缝性介质模型的渗吸过程与规律,预测渗吸作用下饱和度变化和采收率变化。
2.根据权利要求1所述的裂缝性油藏逆向渗吸采收率预测方法,其特征在于:步骤(一)中,建立裂缝性介质模型的步骤为:
(一)搜集反映裂缝发育状况的岩心、地震、成像测井的资料;
(二)根据岩心统计结果、地震预测结果、成像测井检测结果,统计裂缝发育状况、裂缝发育密度、裂缝长度、裂缝开度、裂缝倾角的信息,获得大尺度裂缝网络的连通性信息、裂缝切割基质岩块形状信息以及裂缝与基质的接触关系信息;
(三)通过CT扫描、扫描电镜、高压压汞和核磁共振分析小尺度裂缝的连通性,确定基质内部连通性差、或孤立裂缝发育状况、或连通性差和孤立裂缝发育状况;
(四)根据裂缝网络对基质岩块的切割尺寸和接触关系,以及连通性差、或孤立裂缝发育状况、或连通性差和孤立裂缝发育状况,确定裂缝性介质模型。
3.根据权利要求1所述的裂缝性油藏逆向渗吸采收率预测方法,其特征在于:步骤(三)中,基质的渗吸方程表达式为:
式中,vm,w为基质内水相流体流速,Pm,c为基质中毛管力,Km为基质的绝对渗透率,km,rw为基质中水相流体相对渗透率,km,ro为基质中油相流体相对渗透率,μw为水相流体年度,μo为油相流体粘度;
进一步表示为:
其中,Sm,w为基岩含水饱和度,基岩内分流量fm,w为:
定义基岩内毛管力扩散系数Dm,e为:
则有:
同时满足质量守恒,则有:
采用混合有限元法求解基质的渗吸方程式(5)和(6)获得基质有限单元上的渗吸方程表达式为:
式中,SmwA,e为基质有限单元e上面含水饱和度,为矩阵Ae的逆矩阵,为基质有限单元e在边li上边含水饱和度,i=1,2,3,φm为基质孔隙度;
裂缝的渗吸方程表达式为:
式中,vf,w为裂缝内水相流体流速,Pf,c为裂缝中毛管力,Kf为裂缝的绝对渗透率,kf,ro为裂缝中油相流体相对渗透率,μo为油相流体粘度,Sf,w为裂缝含水饱和度,裂缝内分流量ff,w为:
定义裂缝内毛管力扩散系数Df,e为:
其数值的大小与裂缝开度有关,当裂缝开度大于0.1mm,因裂缝内部毛管力较小,Df,e数值很小,当裂缝开度小于等于0.1mm时,其内部毛管力不可忽略,则有:
同时满足质量守恒,则有:
采用混合有限元法求解裂缝的渗吸方程式(12)和(13)获得裂缝有限单元上的渗吸方程表达式为:
式中,Sfwl,l为裂缝有限单元l上边含水饱和度,为矩阵Ll的逆矩阵,为裂缝有限单元l在端点di上含水饱和度,i=1,2。
4.根据权利要求1或3任意一项所述的裂缝性油藏逆向渗吸采收率预测方法,其特征在于:步骤(三)中,基质有限单元之间的耦合过程为:定义在两个相邻有限单元之间的边界上,通过边界流出第一个有限单元的流量等同于通过边界流入第二个有限单元的流量,第一个有限单元在边界上的含水饱和度等于第二个有限单元在边界上的含水饱和度,如此完成基质有限单元之间的耦合;基质有限单元与裂缝有限单元之间的耦合过程为:定义在裂缝有限单元与相邻的两个基质有限单元之间,流出裂缝有限单元的流量等于通过边界流入到两个基质有限单元的流量之和,如此完成基质有限单元与裂缝有限单元之间的耦合;裂缝有限单元之间的耦合过程为:定义在裂缝有限单元的节点处,经过该节点流出第一个裂缝有限单元的流量等于经过该节点流入到第二个裂缝有限单元的流量,在裂缝有限单元节点处,各个裂缝有限单元在节点上的含水饱和度相等,如此完成裂缝有限单元之间的耦合。
5.根据权利要求4所述的裂缝性油藏逆向渗吸采收率预测方法,其特征在于:步骤(三)中,耦合后获得的整个有限元方程为:
式中,m为基质,f为裂缝,Dm、Bm,m、Bm,f、Bf,f、Mm,m、Mm,f、Mf,f、Df、Mf均为系数矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610060046.1A CN105676309B (zh) | 2016-01-28 | 2016-01-28 | 裂缝性油藏逆向渗吸采收率预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610060046.1A CN105676309B (zh) | 2016-01-28 | 2016-01-28 | 裂缝性油藏逆向渗吸采收率预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105676309A true CN105676309A (zh) | 2016-06-15 |
CN105676309B CN105676309B (zh) | 2017-10-27 |
Family
ID=56302974
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610060046.1A Expired - Fee Related CN105676309B (zh) | 2016-01-28 | 2016-01-28 | 裂缝性油藏逆向渗吸采收率预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105676309B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108596383A (zh) * | 2018-04-18 | 2018-09-28 | 中国石油大学(北京) | 储层分类的方法及装置 |
CN109884269A (zh) * | 2019-02-02 | 2019-06-14 | 中国石油大学(北京) | 岩心自发渗吸采收率预测方法及系统 |
CN110306960A (zh) * | 2019-05-22 | 2019-10-08 | 西安石油大学 | 双孔双渗介质储层的压裂液渗吸增产模拟方法及装置 |
CN111155974A (zh) * | 2020-02-13 | 2020-05-15 | 西南石油大学 | 缝洞型油藏注水替油的油水置换状态预测方法 |
CN111274528A (zh) * | 2020-03-02 | 2020-06-12 | 中国石油大学(北京) | 储层裂缝渗吸质量预测方法及系统 |
CN112855108A (zh) * | 2021-03-18 | 2021-05-28 | 中国地质大学(北京) | 致密储层滑溜水压裂液渗吸采收率预测方法以及预测装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101424185A (zh) * | 2008-11-24 | 2009-05-06 | 大庆油田有限责任公司 | 注采主流线与分流线油藏剩余油参数特征及关系的地球化学测定方法 |
WO2012112736A2 (en) * | 2011-02-17 | 2012-08-23 | Chevron U.S.A. Inc. | System and method for uncertainty quantification in reservoir simulation |
CN103362500A (zh) * | 2013-08-06 | 2013-10-23 | 中国石油大学(华东) | 碳酸盐岩缝洞型油藏裂缝、溶洞随钻检测系统及方法 |
US20140350902A1 (en) * | 2013-05-21 | 2014-11-27 | IFP Energies Nouvelles | Method for exploiting a fractured medium on the basis of a matched reservoir model for wells chosen by means of an equivalent transmissivity model |
-
2016
- 2016-01-28 CN CN201610060046.1A patent/CN105676309B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101424185A (zh) * | 2008-11-24 | 2009-05-06 | 大庆油田有限责任公司 | 注采主流线与分流线油藏剩余油参数特征及关系的地球化学测定方法 |
WO2012112736A2 (en) * | 2011-02-17 | 2012-08-23 | Chevron U.S.A. Inc. | System and method for uncertainty quantification in reservoir simulation |
US20140350902A1 (en) * | 2013-05-21 | 2014-11-27 | IFP Energies Nouvelles | Method for exploiting a fractured medium on the basis of a matched reservoir model for wells chosen by means of an equivalent transmissivity model |
CN103362500A (zh) * | 2013-08-06 | 2013-10-23 | 中国石油大学(华东) | 碳酸盐岩缝洞型油藏裂缝、溶洞随钻检测系统及方法 |
Non-Patent Citations (3)
Title |
---|
姚军 等: ""天然裂缝性油藏的离散裂缝网络数值模拟方法"", 《石油学报》 * |
张世明 等: ""基于离散裂缝模型的裂缝性介质等效渗透率求解新方法"", 《科学技术与工程》 * |
张娜 等: ""裂缝介质渗流多尺度混合有限元数值模拟研究"", 《力学季刊》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108596383A (zh) * | 2018-04-18 | 2018-09-28 | 中国石油大学(北京) | 储层分类的方法及装置 |
CN108596383B (zh) * | 2018-04-18 | 2021-11-30 | 中国石油大学(北京) | 储层分类的方法及装置 |
CN109884269A (zh) * | 2019-02-02 | 2019-06-14 | 中国石油大学(北京) | 岩心自发渗吸采收率预测方法及系统 |
CN110306960A (zh) * | 2019-05-22 | 2019-10-08 | 西安石油大学 | 双孔双渗介质储层的压裂液渗吸增产模拟方法及装置 |
CN111155974A (zh) * | 2020-02-13 | 2020-05-15 | 西南石油大学 | 缝洞型油藏注水替油的油水置换状态预测方法 |
CN111274528A (zh) * | 2020-03-02 | 2020-06-12 | 中国石油大学(北京) | 储层裂缝渗吸质量预测方法及系统 |
CN111274528B (zh) * | 2020-03-02 | 2021-09-17 | 中国石油大学(北京) | 储层裂缝渗吸质量预测方法及系统 |
CN112855108A (zh) * | 2021-03-18 | 2021-05-28 | 中国地质大学(北京) | 致密储层滑溜水压裂液渗吸采收率预测方法以及预测装置 |
CN112855108B (zh) * | 2021-03-18 | 2022-02-15 | 中国地质大学(北京) | 致密储层滑溜水压裂液渗吸采收率预测方法以及预测装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105676309B (zh) | 2017-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105676309A (zh) | 裂缝性油藏逆向渗吸采收率预测方法 | |
CN104101562B (zh) | 一种确定储层渗透率的方法 | |
CN104634804B (zh) | 一种利用核磁共振t2谱确定储层相对渗透率的方法 | |
CN104453874B (zh) | 一种基于核磁共振的砂砾岩储层含油饱和度的计算方法 | |
CN103198363B (zh) | 一种基于ct孔隙分析的储层产气量预测方法及装置 | |
CN105696986B (zh) | 一种新型复合驱驱油实验/试验模拟方法 | |
CN108763711A (zh) | 一种基于岩心扫描图像分块数值模拟的渗透率预测方法 | |
CN105279790A (zh) | 裂缝网络三维数字岩心建模方法 | |
CN102175832A (zh) | 一种确定典型储层最佳饱和度计算模型的方法 | |
CN109100812B (zh) | 基于核磁共振的岩石孔隙分形维数评价方法及装置 | |
CN103852787A (zh) | 一种砂岩储层成岩地震相表征方法 | |
CN104729972A (zh) | 确定碳酸盐岩气水相对渗透率的方法及装置 | |
CN104751473A (zh) | 碳酸岩多尺度孔洞特征的确定方法及装置 | |
CN112098155B (zh) | 油藏油水岩反应实验装置、方法以及取样位置的确定方法 | |
CN105134185A (zh) | 储层流体性质识别方法 | |
CN106126936B (zh) | 一种致密低渗透储层裂缝有效性的综合评价方法 | |
CN105257284B (zh) | 一种利用元素俘获能谱测井确定凝灰质含量的方法及装置 | |
CN107808049B (zh) | 基于多孔介质三维微观结构模型的dnapl迁移数值模拟方法 | |
CN106897531A (zh) | 一种低渗透石灰岩储层渗透率的定量评价方法 | |
CN104374679B (zh) | 一种确定储层渗透率贡献值的方法 | |
CN106204302A (zh) | 一种计算原始含水饱和度的方法与应用 | |
CN106771071B (zh) | 一种基于油水相渗的密闭取心饱和度校正方法 | |
CN105844011A (zh) | 一种基于毛管模型的渗透率计算方法 | |
Baojun et al. | GIS-based quantitative analysis of orientation anisotropy of contaminant barrier particles using standard deviational ellipse | |
CN105590018A (zh) | 砂泥岩薄互层油藏的油水层判别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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: 20171027 |
|
CF01 | Termination of patent right due to non-payment of annual fee |