CN110532592B - 一种缝洞油气藏压裂井大溶洞试井解释方法 - Google Patents
一种缝洞油气藏压裂井大溶洞试井解释方法 Download PDFInfo
- Publication number
- CN110532592B CN110532592B CN201910635926.0A CN201910635926A CN110532592B CN 110532592 B CN110532592 B CN 110532592B CN 201910635926 A CN201910635926 A CN 201910635926A CN 110532592 B CN110532592 B CN 110532592B
- Authority
- CN
- China
- Prior art keywords
- karst cave
- pressure
- fracture
- dimensionless
- well
- 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.)
- Active
Links
- 238000012360 testing method Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 29
- 239000004215 Carbon black (E152) Substances 0.000 title claims description 10
- 229930195733 hydrocarbon Natural products 0.000 title claims description 10
- 150000002430 hydrocarbons Chemical class 0.000 title claims description 10
- 238000013178 mathematical model Methods 0.000 claims abstract description 32
- 230000008878 coupling Effects 0.000 claims abstract description 12
- 238000010168 coupling process Methods 0.000 claims abstract description 12
- 238000005859 coupling reaction Methods 0.000 claims abstract description 12
- 230000009466 transformation Effects 0.000 claims abstract description 11
- 239000003921 oil Substances 0.000 claims description 72
- 238000004519 manufacturing process Methods 0.000 claims description 27
- 239000010779 crude oil Substances 0.000 claims description 18
- 230000006870 function Effects 0.000 claims description 18
- 238000003860 storage Methods 0.000 claims description 13
- 239000011148 porous material Substances 0.000 claims description 12
- 230000000694 effects Effects 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 230000007797 corrosion Effects 0.000 claims description 9
- 238000005260 corrosion Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 239000012530 fluid Substances 0.000 claims description 8
- 230000002500 effect on skin Effects 0.000 claims description 7
- 230000014509 gene expression Effects 0.000 claims description 7
- 230000006835 compression Effects 0.000 claims description 6
- 238000007906 compression Methods 0.000 claims description 6
- 238000009795 derivation Methods 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 6
- 238000004090 dissolution Methods 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 241000232299 Ralstonia Species 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000011065 in-situ storage Methods 0.000 claims description 3
- 230000000704 physical effect Effects 0.000 claims description 2
- 238000010206 sensitivity analysis Methods 0.000 abstract description 2
- 241001282153 Scopelogadus mizolepis Species 0.000 abstract 1
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 239000000243 solution Substances 0.000 description 13
- 230000005465 channeling Effects 0.000 description 5
- 238000005553 drilling Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 235000019994 cava Nutrition 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 150000002500 ions Chemical class 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000003129 oil well Substances 0.000 description 2
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 description 1
- 230000007704 transition Effects 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
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (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)
Abstract
本发明公开了一种缝洞油气藏压裂井大溶洞试井解释方法，包括以下步骤：S1、建立压裂井缝洞油藏试井物理模型；S2、建立溶洞及溶洞泄油区径向渗流无因次数学模型；S3、建立压裂缝线性渗流无因次数学模型及其与径向渗流模型的耦合模型；S4、数学模型拉氏变换并求解；S5、绘制典型试井双对数曲线；S6、实测试井数据的拟合解释；S7、研究溶洞半径的影响；S8、研究主裂缝长度的影响。本发明的有益效果是：建立了缝洞油气藏压裂井大溶洞试井解释方法，考虑井筒压裂形成大尺度压裂缝渗流通道，压裂缝沟通大溶洞，将大尺度溶洞视为势体，以溶洞为中心形成双重介质溶洞泄油区，通过溶洞泄油区的径向渗流与压裂缝的线性流动耦合建立数学模型，对模型进行拉普拉斯变换求解，绘制典型试井双对数曲线，并对模型进行参数敏感性分析。
Description
技术领域
本发明涉及地质勘探技术领域，具体涉及一种缝洞油气藏压裂井大溶洞试井解释方法。
背景技术
缝洞型油气藏中发育有大量的天然裂缝与溶蚀溶孔与大尺度溶洞，这三类介质都是良好的油气储集空间。目前的地震识别与钻井技术还无法完全保证钻遇缝洞储集体。存在生产井未钻遇溶洞的情况，常采用压裂等措施产生人工裂缝以沟通有效缝洞储集体。对于此类压裂井缝洞油藏，识别溶洞体积与压裂缝长度对指导后期生产具有重要的指导意义。
从60年代以来，许多学者提出并发展了多重介质模型碳酸盐岩油气藏试井理论，认为油气藏可视多几种介质系统在空间上的叠加，各系统间存在窜流。然而，多重介质试井模型的基础是连续性假设，此类试井模型的连续性假设对于缝洞发育尺度大且离散分布的缝洞型油气藏并不完全适用，连续介质模型试井图版对缝洞油气藏矿场实测曲线往往难以拟合，且无法解释出大尺度裂缝与大溶洞参数。20世纪以来，一些学者考虑缝洞离散分布，建立了井钻遇溶洞试井模型、缝洞串并联组合等试井模型成功解释出了缝洞参数。有学者建立了井钻遇溶洞，溶洞外为双重连续介质的试井模型并求得了解析解。但这些模型的缝洞组合方式简单，适用范围有限，不能适用于压裂井缝洞油气藏。在缝洞型油气藏的生产实践中，存在生产井未钻遇有效缝洞区，导致油井产能低。综上所述，大尺度缝洞发育的缝洞型油气藏不完全多重介质连续性假设，现有的钻遇溶洞模型和缝洞组合试井模型又无法解决压裂油气藏试井问题。
本发明在研究压裂油气藏缝洞发育特征基础上，建立了一种新的试井模型与方法，模型考虑溶洞外区为双重介质，引入溶洞泄油区概念，对压裂缝中的线性渗流与溶洞径向渗流区中的径向渗流进行耦合建立数学模型，对模型进行了耦合求解，绘制了典型试井曲线，分析了流动特征并进行了参数敏感性分析。
发明内容
本发明的目的在于克服现有技术的缺点，提供一种建立了缝洞油气藏压裂井大溶洞试井模型与方法，考虑溶洞外区为双重介质，引入溶洞泄油区概念，对压裂缝中的线性渗流与溶洞径向渗流区中的径向渗流进行耦合建立数学模型，对模型进行了耦合求解，绘制了典型试井曲线，分析了流动特征并进行了参数敏感性分析。
本发明的目的通过以下技术方案来实现：一种缝洞油气藏压裂井大溶洞试井解释模型与方法，它包括以下步骤：
包括以下步骤：
S1、建立压裂井缝洞油藏试井物理模型；
S2、建立溶洞及溶洞泄油区径向渗流无因次数学模型；
S3、建立压裂缝线性渗流无因次数学模型及其与径向渗流模型的耦合模型；
S4、数学模型拉氏变换并求解；
S5、绘制典型试井双对数曲线；
S6、实测试井数据的拟合方法，编写解释程序；
S7、研究溶洞半径的影响；
S8、研究主裂缝长度的影响。
所述步骤S1中物理模型的基本假设条件如下：
(1)以溶洞为中心形成泄油区，溶洞外发育有连续分布的天然裂缝和溶蚀溶孔，视为连续介质；(2)地层中有一口直井以定产量生产，井筒通过压裂形成的主裂缝与溶洞相连，储层流体为单相原油，溶洞被原油充填；(3)溶洞泄油区中，原油通过溶蚀溶孔向裂缝窜流，再向大溶洞发生窜流，最后通过压裂缝流入井筒；(4)原油在储层裂缝中的流动符合达西温流动；忽略毛管力与重力的作用；(5)初始时刻地层压力处处相，于原始地层压力；(6)溶洞、裂缝、溶孔和流体均为微可压缩，压缩系数为常数；(7)考虑井筒储集效应与表皮效应的影响。
在所述步骤S2中建立溶洞及溶洞泄油区径向渗流无因次数学模型的步骤如下：
模型以油藏模型为例，若为气藏，则将模型中的所有压力参数P替换为气藏拟压力,下标与下标含义均不变，流体物性参数采用对应气体参数即可。所述拟压力计算方法为：
式中，ψ为真实气体拟压力(对应于压力p),Mpa2/(mPa.s)；p为压力，MPa；μ为气体粘度，mPa.s；Z为气体偏差因子；p0为任意选定某一参考压力。
定义如下无因次量：
无因次压力pfD、pfrD、pvrD、pvD：
式中：pfD、pfrD、pvrD分别为压裂缝、径向泄油区天然裂缝、径向泄油区中的溶孔无因次压力；pf、pfr、pvr分别为压裂缝、径向泄油区天然裂缝、径向泄油区中的溶孔压力，MPa；
pv为大溶洞等势体压力，MPa；pvD为大溶洞等势体无因次压力。pi为地层原始压力，MPa；h为储层厚度，m；qo为生产井产量，m3/d；Bo为原油体积系数，m3/m3；μ为原油黏度，mPa.s；Kf为压裂缝渗透率，μm2。
无因次生产时间tD：
式中：t为井生产时间，h；为井无因次生产时间；φf为压裂缝孔隙度，无因次；Cft为压裂缝综合压缩系数，MPa-1，rw为生产井半径，m。
无因次地层厚度hD：
无因次压裂缝距离xD：
式中：x为主裂缝上其点距井口坐标，m；为主裂缝上点距井口无因次坐标；
无因次溶洞泄油区径向距离rD：
式中：r为径向溶洞泄油区点距溶洞中心坐标，m；
无因次泄油区半径reD：
式中：re为溶洞泄油区半径。
无因次溶洞半径rvD：
式中：rv为溶洞半径，m；
弹性储容比ω:
式中：φfr、φvr分别为径向天然裂缝、溶孔孔隙度，无因次；Cfrt、Cvrt为径向天然裂缝、溶孔综合压缩系数，MPa-1；
无因次窜流系数λ
式中：α为形状因子，m-2；Kfr、Kvr分别为径向天然裂缝、溶孔渗透率，μm2。
无因次井筒储集系数CD
式中：C为井筒储集系数，MPa-1。
无因次压裂缝横截面积AD：
Af为压裂缝横截面积，m2。
径向流区与线性流区流动系数之比MD：
模型自定无因次参数η1、η2：
溶洞以外泄油区裂缝与溶孔微分控制方程：
外边界条件包括无限大条件、封闭边界条件和定压边界条件：
1)无限大：
pfrD(∞,tD)＝0 (3a)
2)封闭边界：
3)定压边界：
pfrD(reD,tD)＝0 (3c)
内边界条件：
初始条件：
在所述步骤S3中建立压裂缝线性渗流无因次数学模型及其与径向渗流模型的耦合模型的步骤如下：
压裂缝微分控制方程：
内边界条件：
外边界条件(耦合条件)：
初始条件：
在所述步骤S4中数学模型拉氏变换并求解的步骤如下：
将径向泄油区无因次数学模型拉氏变换可得：
外边界条件：
1)无限大：
2)封闭边界：
3)定压边界：
内边界条件：
将主裂缝数学模型进行拉氏变换可得：
式中，s为拉氏变量。
将式(11)代入式(10)可得：
其中：
式中，f(s)为特征函数。
为简化方程，令：
式中，σ为自定中间变量。
溶洞径向渗流控制方程(17)的通解为：
对其求导可得：
式中，A0、B0为待定参数，I0为零阶第一类修正贝塞尔函数，K0为零阶第二类修正贝塞尔函数，I1为一阶第一类修正贝塞尔函数，K1为一阶第二类修正贝塞尔函数。
将式(20)代入内边界条件式(13)可得：
以封闭外边界为例，将式(21)代入外边界条件(25b)可得：
A0σI1(reDσ)-B0σK1(reDσ)＝0 (23)
裂缝的渗流微分控制方程(14)求通解可得：
式中：A1、B1为待定参数。
对式(24)求导可得：
将式(25)代入压裂缝内边界条件(15)可得：
再将式(25)代入外边界条件(16)可得：
溶洞压力可表示为：
式中，xfD为无因次压裂缝长度。
共有5个方程式(22)、(23)、(26)、(27)、(28)，和5个未知数(A0,B0,A1,B1,pvD)，联立求解出这五个参数的表达式；
当xfD＝rwD＝1时，井底流压可表示为：
存在井筒存储效应和表皮效应时，利用Duhamel原理结合井储和表皮系数的定义计算出生产井以定产量生产时的井底流压表达式：
在所述步骤S5中绘制典型试井双对数曲线的步骤如下：采用stehfest数值反演方法，将计算所得的式(30)进行反演计算，可得到真实空间中压力与压力导数典型试井曲线；
在所述步骤S6中理论曲线与实测数据拟合方法的步骤如下：将现场实测井底压力和导数曲线和建立的试井解释模型计算的曲线对比，其误差作为目标函数：
式中：pc为模型计算得到的井底压力，MPa；pi *为现场实测井底压力，MPa；n为实验数据个数。
在所述步骤S7中研究溶洞半径的影响步骤如下：改变模型的溶洞半径并保持其他参数不变，得到在不同溶洞半径下的无因次压力与压力导数的对比曲线；
在所述步骤S8中研究主裂缝长度的影响步骤如下：改变模型的主裂缝长度并保持其他参数不变，得到在不同主裂缝长度下的无因次压力与压力导数的对比曲线。
本发明具有以下优点：
1、本方法针对缝洞型油藏建立了缝洞油藏压裂井大溶洞试井模，考虑在井筒附近形成有压裂缝渗流优势通道，通过压裂缝沟通大溶洞，在溶洞周围，发育有天然裂缝和溶蚀溶孔，采用双重连续介质油藏模型，将大尺度溶洞视为势体，以溶洞为中心形成溶洞泄油区，通过溶洞泄油区的径向注与主裂缝的线性流动耦合进行模型求解，对数学模型无因次化并进行拉普拉斯变换求解；
2、通过数值反演得到无因次压力与压力导数试井典型图版，编制相应试井解释软件对现场实测数据进行拟合解释。
3、本方法能够解释方法能适用于更多缝洞组合类型的油藏，从而进一步提高油井产能。
附图说明
图1为压裂井油藏示意图；
图2为压裂井试井物理模型图；
图3为真实空间中压力与压力导数典型试井曲线图；
图4为程序设计流程图；
图5为在不同溶洞半径rvD下的无因次压力与压力导数的对比曲线图；
图6为在不同主裂缝长度x1D下的无因次压力与压力导数的对比曲线。
具体实施方式
下面结合附图对本发明做进一步的描述，本发明的保护范围不局限于以下所述：
一种缝洞油气藏压裂井大溶洞试井解释方法，它包括以下步骤：
S1、建立压裂井试井物理模型：
在径向油藏中，生产井通过压裂缝连接大溶洞，以溶洞为中心存在双重介质泄油区，如图1所示为压裂井油藏示意图；物理模型的基本假设条件如下：
(1)以溶洞为中心形成泄油区，溶洞外发育有连续分布的天然裂缝和溶蚀溶孔，视为连续介质；(2)地层中有一口直井以定产量生产，井筒通过压裂形成的主裂缝与溶洞相连，储层流体为单相原油，溶洞被原油充填；(3)溶洞泄油区中，原油通过溶蚀溶孔向裂缝窜流，再向大溶洞发生窜流，最后通过压裂缝流入井筒；(4)原油在储层裂缝中的流动符合达西温流动；忽略毛管力与重力的作用；(5)初始时刻地层压力处处相，于原始地层压力；(6)溶洞、裂缝、溶孔和流体均为微可压缩，压缩系数为常数；(7)考虑井筒储集效应与表皮效应的影响。
S2、建立溶洞及溶洞泄油区径向渗流无因次数学模型：
模型与方法以油藏模型为例，若为气藏，则将模型中的所有压力参数P替换为气藏拟压力ψ即可，下标与下标含义均不变，所述拟压力计算方法为：
式中，ψ为真实气体拟压力(对应于压力p),Mpa2/(mPa.s)；p为压力，MPa；μ为气体粘度，mPa.s；Z为气体偏差因子；p0为任意选定某一参考压力。
定义如下无因次量：
无因次压力pfD、pfrD、pvrD、pvD：
式中：pfD、pfrD、pvrD分别为压裂缝、径向泄油区天然裂缝、径向泄油区中的溶孔无因次压力；pf、pfr、pvr分别为压裂缝、径向泄油区天然裂缝、径向泄油区中的溶孔压力，MPa；
pv为大溶洞等势体压力，MPa；pvD为大溶洞等势体无因次压力。pi为地层原始压力，MPa；h为储层厚度，m；qo为生产井产量，m3/d；Bo为原油体积系数，m3/m3；μ为原油黏度，mPa.s；Kf为压裂缝渗透率，μm2。
无因次生产时间tD：
式中：t为井生产时间，h；为井无因次生产时间；φf为压裂缝孔隙度，无因次；Cft为压裂缝综合压缩系数，MPa-1，rw为生产井半径，m。
无因次地层厚度hD：
无因次压裂缝距离xD：
式中：x为主裂缝上其点距井口坐标，m；为主裂缝上点距井口无因次坐标；
无因次溶洞泄油区径向距离rD：
式中：r为径向溶洞泄油区点距溶洞中心坐标，m；
无因次泄油区半径reD：
式中：re为溶洞泄油区半径。
无因次溶洞半径rvD：
式中：rv为溶洞半径，m；
弹性储容比ω:
式中：φfr、φvr分别为径向天然裂缝、溶孔孔隙度，无因次；Cfrt、Cvrt为径向天然裂缝、溶孔综合压缩系数，MPa-1；
无因次窜流系数λ
式中：α为形状因子，m-2；Kfr、Kvr分别为径向天然裂缝、溶孔渗透率，μm2。
无因次井筒储集系数CD
式中：C为井筒储集系数，MPa-1。
无因次压裂缝横截面积AD：
Af为压裂缝横截面积，m2。
径向流区与线性流区流动系数之比MD：
模型自定无因次参数η1、η2：
溶洞以外泄油区裂缝与溶孔微分控制方程：
外边界条件包括无限大条件、封闭边界条件和定压边界条件：
1)无限大：
pfrD(∞,tD)＝0 (3a)
2)封闭边界：
3)定压边界：
pfrD(reD,tD)＝0 (3c)
内边界条件：
初始条件：
S3、建立压裂缝线性渗流无因次数学模型及其与径向渗流模型的耦合模型：
压裂缝微分控制方程：
内边界条件：
外边界条件(耦合条件)：
初始条件：
S4、数学模型拉氏变换并求解
将径向泄油区无因次数学模型拉氏变换可得：
外边界条件：
1)无限大：
2)封闭边界：
3)定压边界：
内边界条件：
将主裂缝数学模型进行拉氏变换可得：
式中，s为拉氏变量。
将式(11)代入式(10)可得：
其中：
式中，f(s)为特征函数。
为简化方程，令：
式中，σ为自定中间变量。
溶洞径向渗流控制方程(17)的通解为：
对其求导可得：
式中，A0、B0为待定参数，I0为零阶第一类修正贝塞尔函数，K0为零阶第二类修正贝塞尔函数，I1为一阶第一类修正贝塞尔函数，K1为一阶第二类修正贝塞尔函数。
将式(20)代入内边界条件式(13)可得：
以封闭外边界为例，将式(21)代入外边界条件(25b)可得：
A0σI1(reDσ)-B0σK1(reDσ)＝0 (23)
裂缝的渗流微分控制方程(14)求通解可得：
式中：A1、B1为待定参数。
对式(24)求导可得：
将式(25)代入压裂缝内边界条件(15)可得：
再将式(25)代入外边界条件(16)可得：
溶洞压力可表示为：
式中，xfD为无因次压裂缝长度。
共有5个方程式(22)、(23)、(26)、(27)、(28)，和5个未知数(A0,B0,A1,B1,pvD)，联立求解出这五个参数的表达式；
当xfD＝rwD＝1时，井底流压可表示为：
存在井筒存储效应和表皮效应时，利用Duhamel原理结合井储和表皮系数的定义计算出生产井以定产量生产时的井底流压表达式：
S5、绘制典型图版曲线：
采用stehfest数值反演方法，将上述无因次井底流压表达式(30)进行反演计算，可得到真实空间中压力与压力导数典型试井曲线如图3所示，从附图3中可以看出：主裂缝沟通溶洞模型试井曲线可以划分为九个流动阶段。
第Ⅰ阶段：井筒储集效应阶段，压力与压力导数曲线都呈斜率为1的直线；
第Ⅱ阶段：表皮效应阶段，压力与压力导数曲线向下弯曲；
第Ⅲ阶段：主裂缝线性流动段，压力与压力导数曲线呈斜率为0.5的平等线；
第Ⅳ阶段：大溶洞向主裂缝窜流段，压力导数曲线下凹，形成一个凹子；
第V阶段：溶洞窜流向溶洞泄油区径向流动过渡段，特征为压力导数下降，逐渐过渡到0.5线；
第VI阶段：溶洞外泄油区裂缝径向流动段，压力导数为0.5水平线；
第VII阶段：溶洞外泄油区溶孔向裂缝窜流段，压力导数曲线下凹，形成第二个凹子；
第VIII阶段：系统总径向流动段，压力导数曲线为0.5线；
第IX阶段：边界响应段。无限大地层压力导数呈0.5线，封闭边界压力导数线快速上升，与压力线重合，呈斜率为1的直线，定压边界压力导数线快速下降到零；
S6、试井解释程序的设计：
采用优化算法对现场实测采出示踪剂浓度曲线进行拟合，以解释地层参数。通过遗传算法模拟自然界生物的“优胜劣汰”，进化过程包括遗传、交叉、变异、自然选择等，最终得到满足要求的最优解。将现场实测井底压力和导数曲线和建立的试井解释模型计算的曲线对比，其误差作为目标函数：
式中：pc为模型计算得到的井底压力，MPa；pi *为现场实测井底压力，MPa；n为实验数据个数。
S7、研究溶洞半径的影响：
改变模型的溶洞半径并保持其他参数不变，得到在不同溶洞半径下的无因次压力与压力导数的对比曲线如图5所示，从附图5中可以得出结论：溶洞半径增大，对应溶洞向主裂缝的窜流更容易，压力导数线形成的凹子下凹更深，溶洞半径越大，流动阶段过渡到0.5线越滞后。
S8、研究主裂缝长度的影响：
改变模型的主裂缝长度并保持其他参数不变，得到在不同主裂缝长度下的无因次压力与压力导数的对比曲线如图6所示，从附图6可以得出结论：主裂缝长度增大，主裂缝线性流动时间越长，大溶洞窜流形成的凹子出现的时间越滞后，且在相同溶洞体积情况下，溶洞距离越大，形成的凹子下凹幅度越小。
综上所述，本方法针对缝洞型油藏建立了缝洞油藏压裂井大溶洞试井模，考虑在井筒附近形成有压裂缝渗流优势通道，通过压裂缝沟通大溶洞，在溶洞周围，发育有天然裂缝和溶蚀溶孔，采用双重连续介质油藏模型，将大尺度溶洞视为势体，以溶洞为中心形成溶洞泄油区，通过溶洞泄油区的径向注与主裂缝的线性流动耦合进行模型求解，对数学模型无因次化并进行拉普拉斯变换求解；通过数值反演得到无因次压力与压力导数试井典型图版，编制相应试井解释软件对现场实测数据进行拟合解释。因此能够解释方法能适用于更多缝洞组合类型的油藏，更精确解释此类油气藏的储层参数，从而更好的指导后续的生产。
Claims (7)
1.一种缝洞油气藏压裂井大溶洞试井解释方法，其特征在于，包括以下步骤：
S1、建立压裂井缝洞油藏试井物理模型，物理模型的基本假设条件如下：
(1)以溶洞为中心形成泄油区，溶洞外发育有连续分布的天然裂缝和溶蚀溶孔，视为连续介质；(2)地层中有一口直井以定产量生产，井筒通过压裂形成的主裂缝与溶洞相连，储层流体为单相原油，溶洞被原油充填；(3)溶洞泄油区中，原油通过溶蚀溶孔向裂缝窜流，再向大溶洞发生窜流，最后通过压裂缝流入井筒；(4)原油在储层裂缝中的流动符合达西温流动；忽略毛管力与重力的作用；(5)初始时刻地层压力处处相等于原始地层压力；(6)溶洞、裂缝、溶孔和流体均为微可压缩，压缩系数为常数；(7)考虑井筒储集效应与表皮效应的影响；
S2、建立溶洞及溶洞泄油区径向渗流无因次数学模型，建立溶洞及溶洞泄油区径向渗流无因次数学模型的步骤如下：
模型以油藏模型为例，若为气藏，则将模型中的所有压力参数P替换为气藏拟压力,下标与下标含义均不变，流体物性参数采用对应气体参数即可，所述拟压力计算方法为：
式中，ψ为真实气体拟压力，Mpa2/(mPa.s)；P为压力，MPa；μ为气体粘度，mPa.s；Z为气体偏差因子；P0为任意选定某一参考压力；
定义如下无因次量：
无因次压力pfD、pfrD、pvrD、pvD：
式中：pfD、pfrD、pvrD分别为压裂缝、径向泄油区天然裂缝、径向泄油区中的溶孔无因次压力；pf、pfr、pvr分别为压裂缝、径向泄油区天然裂缝、径向泄油区中的溶孔压力，MPa；
pv为大溶洞等势体压力，MPa；pvD为大溶洞等势体无因次压力； pi为地层原始压力，MPa；h为储层厚度，m；qo为生产井产量，m3/d；Bo为原油体积系数，m3/m3；μ为原油黏度，mPa.s；Kf为压裂缝渗透率，μm2；
无因次生产时间tD：
式中：t为井生产时间，h；为井无因次生产时间；φf为压裂缝孔隙度，无因次；Cft为压裂缝综合压缩系数，MPa-1，rw为生产井半径，m；
无因次地层厚度hD：
无因次压裂缝距离xD：
式中：x为主裂缝上其点距井口坐标，m；为主裂缝上点距井口无因次坐标；
无因次溶洞泄油区径向距离rD：
式中：r为径向溶洞泄油区点距溶洞中心坐标，m；
无因次泄油区半径reD：
式中：re为溶洞泄油区半径；
无因次溶洞半径rvD：
式中：rv为溶洞半径，m；
弹性储容比ω:
式中：φfr、φvr分别为径向天然裂缝、溶孔孔隙度，无因次；Cfrt、Cvrt为径向天然裂缝、溶孔综合压缩系数，MPa-1；
无因次窜流系数λ
式中：α为形状因子，m-2；Kfr、Kvr分别为径向天然裂缝、溶孔渗透率，μm2；
无因次井筒储集系数CD
式中：C为井筒储集系数，MPa-1；
无因次压裂缝横截面积AD：
Af为压裂缝横截面积，m2；
径向流区与线性流区流动系数之比MD：
模型自定无因次参数η1、η2：
溶洞以外泄油区裂缝与溶孔微分控制方程：
外边界条件包括无限大条件、封闭边界条件和定压边界条件：
1)无限大：
pfrD(∞,tD)＝0 (3a)
2)封闭边界：
3)定压边界：
pfrD(reD,tD)＝0 (3c)
内边界条件：
初始条件：
S3、建立压裂缝线性渗流无因次数学模型及其与径向渗流模型的耦合模型；
S4、数学模型拉氏变换并求解；
S5、绘制典型试井双对数曲线；
S6、实测试井数据的拟合方法，编写解释程序；
S7、研究溶洞半径的影响；
S8、研究主裂缝长度的影响。
2.根据权利要求1所述的一种缝洞油气藏压裂井大溶洞试井解释方法，其特征在于：在所述步骤S3中建立压裂缝线性渗流无因次数学模型及其与径向渗流模型的耦合模型的步骤如下：
压裂缝微分控制方程：
内边界条件：
外边界条件：
初始条件：
3.根据权利要求1所述的一种缝洞油气藏压裂井大溶洞试井解释方法，其特征在于：在所述步骤S4中数学模型拉氏变换并求解的步骤如下：
将径向泄油区无因次数学模型拉氏变换可得：
外边界条件：
1)无限大：
2)封闭边界：
3)定压边界：
内边界条件：
将主裂缝数学模型进行拉氏变换可得：
式中，s为拉氏变量；
将式(11)代入式(10)可得：
其中：
式中，f(s)为特征函数；
为简化方程，令：
式中，σ为自定中间变量；
溶洞径向渗流控制方程(17)的通解为：
对其求导可得：
式中，A0、B0为待定参数，I0为零阶第一类修正贝塞尔函数，K0为零阶第二类修正贝塞尔函数，I1为一阶第一类修正贝塞尔函数，K1为一阶第二类修正贝塞尔函数；
将式(20)代入内边界条件式(13)可得：
以封闭外边界为例，将式(21)代入外边界条件(25b)可得：
A0σI1(reDσ)-B0σK1(reDσ)＝0 (23)
裂缝的渗流微分控制方程(14)求通解可得：
式中：A1、B1为待定参数；
对式(24)求导可得：
将式(25)代入压裂缝内边界条件(15)可得：
再将式(25)代入外边界条件(16)可得：
溶洞压力可表示为：
式中，xfD为无因次压裂缝长度；
共有5个方程式(22)、(23)、(26)、(27)、(28)，和5个未知数A0,B0,A1,B1,pvD，联立求解出这五个参数的表达式；
当xfD＝rwD＝1时，井底流压可表示为：
存在井筒存储效应和表皮效应时，利用Duhamel原理结合井储和表皮系数的定义计算出生产井以定产量生产时的井底流压表达式：
4.根据权利要求1所述的一种缝洞油气藏压裂井大溶洞试井解释方法，其特征在于：在所述步骤S5中绘制典型试井双对数曲线的步骤如下：
采用stehfest数值反演方法，将计算所得的式(30)进行反演计算，可得到真实空间中压力与压力导数典型试井曲线。
5.根据权利要求1所述的一种缝洞油气藏压裂井大溶洞试井解释方法，其特征在于：在所述步骤S6中理论曲线与实测数据拟合方法的步骤如下：
将现场实测井底压力和导数曲线和建立的试井解释模型计算的曲线对比，其误差作为目标函数：
式中：pc为模型计算得到的井底压力，MPa；pi *为现场实测井底压力，MPa；n为实验数据个数。
6.根据权利要求1所述的一种缝洞油气藏压裂井大溶洞试井解释方法，其特征在于：在所述步骤S7中研究溶洞半径的影响步骤如下：
改变模型的溶洞半径并保持其他参数不变，得到在不同溶洞半径下的无因次压力与压力导数的对比曲线。
7.根据权利要求1所述的一种缝洞油气藏压裂井大溶洞试井解释方法，其特征在于：在所述步骤S8中研究主裂缝长度的影响步骤如下：
改变模型的主裂缝长度并保持其他参数不变，得到在不同主裂缝长度下的无因次压力与压力导数的对比曲线。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910635926.0A CN110532592B (zh) | 2019-07-15 | 2019-07-15 | 一种缝洞油气藏压裂井大溶洞试井解释方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910635926.0A CN110532592B (zh) | 2019-07-15 | 2019-07-15 | 一种缝洞油气藏压裂井大溶洞试井解释方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110532592A CN110532592A (zh) | 2019-12-03 |
| CN110532592B true CN110532592B (zh) | 2022-04-22 |
Family
ID=68660290
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910635926.0A Active CN110532592B (zh) | 2019-07-15 | 2019-07-15 | 一种缝洞油气藏压裂井大溶洞试井解释方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110532592B (zh) |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113294147A (zh) * | 2020-02-24 | 2021-08-24 | 中国石油化工股份有限公司 | 一种考虑重力因素影响的单洞型断溶体储层试井解释方法 |
| CN111428425B (zh) * | 2020-03-19 | 2020-12-29 | 西南石油大学 | 一种页岩油藏变裂缝渗透率分段压裂水平井产能计算方法 |
| CN111984918A (zh) * | 2020-04-07 | 2020-11-24 | 合肥辰工科技有限公司 | 一种利用压裂施工停泵压力数据获得裂缝高度方法 |
| CN111581786B (zh) * | 2020-04-19 | 2021-02-09 | 东北石油大学 | 用于分析缝洞串联模式双孔复合储层的试井解释模型的试井解释方法 |
| CN111553067B (zh) * | 2020-04-22 | 2022-08-23 | 中国石油大学(北京) | 页岩油藏压裂水平井不停井试井测试和解释方法、装置 |
| CN113626967B (zh) * | 2020-05-07 | 2022-08-26 | 中国石油化工股份有限公司 | 一种考虑应力敏感的缝洞型储层产能确定方法及系统 |
| CN114086934B (zh) * | 2020-08-03 | 2024-03-19 | 中国石油化工股份有限公司 | 一种试井解释方法、装置、存储介质以及计算机设备 |
| CN114086933B (zh) * | 2020-08-03 | 2024-03-19 | 中国石油化工股份有限公司 | 一种试井解释方法、装置、存储介质以及计算机设备 |
| CN112100890B (zh) * | 2020-09-10 | 2022-08-23 | 中国石油天然气股份有限公司 | 一种压裂水平缝起裂扩展全三维数学模型的计算方法 |
| CN114427427B (zh) * | 2020-09-24 | 2024-06-18 | 中国石油化工股份有限公司 | 一种提高缝洞型油藏波及范围的缩缝分流方法 |
| CN112257147B (zh) * | 2020-10-13 | 2022-08-23 | 中国石油天然气股份有限公司 | 一种压驱井合理焖井时间评价方法 |
| CN112377182B (zh) * | 2020-11-30 | 2023-08-22 | 中国石油天然气股份有限公司 | 大型孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 |
| CN112377185B (zh) * | 2020-11-30 | 2023-08-22 | 中国石油天然气股份有限公司 | 单洞与多重介质复合型储层参数敏感性分析方法及装置 |
| CN113111582B (zh) * | 2021-04-12 | 2022-08-23 | 中国石油大学(北京) | 基于深度学习的页岩储层试井智能解释分析方法及装置 |
| CN114065113B (zh) * | 2021-11-19 | 2024-02-27 | 西北大学 | 一种自然裂缝储层多段压裂水平井干扰试井分析方法 |
| CN114925632B (zh) * | 2022-05-26 | 2023-09-01 | 西南石油大学 | 一种缝洞型气藏产能测试的动态模拟方法 |
| CN115828787B (zh) * | 2022-12-13 | 2023-09-08 | 西南石油大学 | 一种改进的mdh试井分析方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016192077A1 (zh) * | 2015-06-04 | 2016-12-08 | 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 | 一种致密气压裂水平井数值试井模型建立求解方法 |
| CN107563027A (zh) * | 2017-08-21 | 2018-01-09 | 西南石油大学 | 用于分析多缝洞单元串联结构的试井解释模型与方法 |
| CN108798634A (zh) * | 2018-07-09 | 2018-11-13 | 合肥辰工科技有限公司 | 一种缝洞型储层单洞-单层均值试井解释方法 |
| CN108920849A (zh) * | 2018-07-09 | 2018-11-30 | 合肥辰工科技有限公司 | 一种缝洞型储层试井确定地层中洞流量方法 |
| CN109138974A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏离散数值试井分析方法及系统 |
-
2019
- 2019-07-15 CN CN201910635926.0A patent/CN110532592B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016192077A1 (zh) * | 2015-06-04 | 2016-12-08 | 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 | 一种致密气压裂水平井数值试井模型建立求解方法 |
| CN109138974A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏离散数值试井分析方法及系统 |
| CN107563027A (zh) * | 2017-08-21 | 2018-01-09 | 西南石油大学 | 用于分析多缝洞单元串联结构的试井解释模型与方法 |
| CN108798634A (zh) * | 2018-07-09 | 2018-11-13 | 合肥辰工科技有限公司 | 一种缝洞型储层单洞-单层均值试井解释方法 |
| CN108920849A (zh) * | 2018-07-09 | 2018-11-30 | 合肥辰工科技有限公司 | 一种缝洞型储层试井确定地层中洞流量方法 |
Non-Patent Citations (2)
| Title |
|---|
| Dual Porosity and Dual Permeability Modeling of Horizontal Well in Naturally Fractured Reservoir;Ren-Shi Nie 等;《Transport in Porous Media》;20111102;第92卷;213-235 * |
| 大尺度溶洞发育的缝洞型油藏试井解释模型研究;尹洪军 等;《特种油气藏》;20180802;第25卷(第5期);84-88 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110532592A (zh) | 2019-12-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110532592B (zh) | 一种缝洞油气藏压裂井大溶洞试井解释方法 | |
| CN104533370B (zh) | 压裂水平井油藏、裂缝、井筒全耦合模拟方法 | |
| CN107701180B (zh) | 一种基于密闭取心的原始油藏含水饱和度计算方法 | |
| WO2016192077A1 (zh) | 一种致密气压裂水平井数值试井模型建立求解方法 | |
| US10495782B2 (en) | System, method and computer program product for determining placement of perforation intervals using facies, fluid boundaries, geobodies and dynamic fluid properties | |
| CN110593865B (zh) | 一种油藏缝洞特征参数试井解释方法 | |
| CN105089582A (zh) | 基于井下流量控制设备的油藏数值模拟方法及装置 | |
| CN112682016B (zh) | 油气田薄互层储层穿层压裂参数的确定方法及压裂方法 | |
| CN105822298A (zh) | 基于产气指数的页岩气层绝对无阻流量的获取方法 | |
| CN110656915B (zh) | 一种页岩气多段压裂水平井多工作制度产能预测方法 | |
| CN110107277B (zh) | 一种获取碳酸盐岩油藏中井钻遇溶洞的体积的方法 | |
| CN112541287A (zh) | 疏松砂岩压裂充填防砂增产调剖一体化设计方法 | |
| CN111339481B (zh) | 一种煤层气井的试井分析方法 | |
| CN115587674A (zh) | 油藏改建储气库扩容达产过程气井动态产能预测方法 | |
| CN114580100A (zh) | 压裂水平井全井筒压力计算方法、设备和计算机可读储存介质 | |
| CN112304838B (zh) | 基于岩心三轴压缩的油页岩储层可改造性室内评价方法 | |
| CN111950111B (zh) | 一种适用于底部开放的碳酸盐岩储层动态分析方法 | |
| CN112836442A (zh) | 一种页岩油水平井网老井水力裂缝注液量的确定方法 | |
| CN115017841B (zh) | 一种确定断溶体合采储层缝洞空间结构的方法及系统 | |
| CN107704646B (zh) | 一种致密储层体积改造后的建模方法 | |
| CN116738672A (zh) | 一种复杂生物礁底水气藏数值模拟模型的建立方法 | |
| CN115559714A (zh) | 一种多分支断溶体油藏缝洞参数的确定方法及系统 | |
| CN112035993A (zh) | 一种底部定压的碳酸盐岩储层测试评价方法 | |
| CN113449353B (zh) | 考虑地层流体渗流作用的水平井井壁稳定分析方法 | |
| CN113919111B (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 |