CN110210157B - 一种页岩气藏压裂水平井非稳态扩散下产能计算方法 - Google Patents
一种页岩气藏压裂水平井非稳态扩散下产能计算方法 Download PDFInfo
- Publication number
- CN110210157B CN110210157B CN201910498213.4A CN201910498213A CN110210157B CN 110210157 B CN110210157 B CN 110210157B CN 201910498213 A CN201910498213 A CN 201910498213A CN 110210157 B CN110210157 B CN 110210157B
- Authority
- CN
- China
- Prior art keywords
- hydraulic fracture
- fracture
- shale gas
- mth
- pressure
- 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
- 238000009792 diffusion process Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000001179 sorption measurement Methods 0.000 claims abstract description 12
- 230000035699 permeability Effects 0.000 claims abstract description 11
- 238000003795 desorption Methods 0.000 claims abstract description 9
- 238000003860 storage Methods 0.000 claims abstract description 9
- 230000004048 modification Effects 0.000 claims abstract description 5
- 238000012986 modification Methods 0.000 claims abstract description 5
- 230000005465 channeling Effects 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000004519 manufacturing process Methods 0.000 claims description 18
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 239000007789 gas Substances 0.000 description 93
- 239000011159 matrix material Substances 0.000 description 13
- 238000004458 analytical method Methods 0.000 description 7
- 230000007246 mechanism Effects 0.000 description 6
- 230000005012 migration Effects 0.000 description 6
- 238000013508 migration Methods 0.000 description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000003337 fertilizer Substances 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000000007 visual effect Effects 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
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- 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
- E21B47/00—Survey of boreholes or wells
- E21B47/10—Locating fluid leaks, intrusions or movements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种页岩气藏压裂水平井非稳态扩散下产能计算方法，包括：（1）收集压裂改造所形成的水力裂缝参数，计算页岩气藏弹性储容比、非稳态扩散窜流系数、无因次渗透率模量、页岩气吸附解吸系数和水力裂缝有限导流能力；（2）在直角坐标系中确定水力裂缝位置，离散水力裂缝并确定各裂缝微元的坐标，根据点源函数理论和压降叠加原理，求取每条水力裂缝缝尖处的压力；（3）根据等效面积原则建立水力裂缝渗流模型，得到每条水力裂缝在水平井筒处的压力；（4）计算在非稳态扩散下的页岩气藏压裂水平井的产能。本发明原理可靠，操作简便，能够在给定页岩储层参数和水力裂缝参数的前提下，准确预测非稳态扩散方式下页岩气藏压裂水平井的产能。
Description
技术领域
本发明涉及油气田开发领域页岩气藏压裂水平井的产能计算方法，特别是涉及一种非稳态扩散方式下页岩气藏压裂水平井的产能计算方法。
背景技术
随着我国能源需求矛盾日益增强，页岩气作为非常规油气资源的有效补充，备受关注，也成为我国天然气增储上产的关键一环。页岩气在储层中赋存方式多样，运移机理复杂，通常包括吸附解吸、扩散、渗流等多重运移机制，而常规的渗流理论也不适合于页岩气藏，正确且清楚地认识页岩气藏的渗流规律及正确表征页岩气运移机制对指导页岩气藏开发具有重要的意义。基质中的页岩气通过扩散的方式运移到天然裂缝系统，扩散的方式分为拟稳态和非稳态扩散，而实际气藏的开发过程中几乎不会存在拟稳态的扩散方式，气体从基质向天然裂缝运移更符合非稳态扩散方式。由于页岩的低孔、超低渗等特征，水平井分段多簇压裂方法成为了常用的增产措施，而压裂后页岩气藏水平井产能的准确预测是评估压裂改造有效性及页岩气藏是否成功开发的关键，通常建立页岩气藏产能模型的方法有解析法、数值模拟法和半解析法，数值模拟法主要采用计算机程序来近似求解，初期所需气藏参数多，求解迭代时间长，对硬件要求高等特点限制应用；解析法对实际情况进行简化，通过高度抽象、理想化物理模型，利用各种数学方法对模型进行直接求解，产能预测的可靠性主要取决于物理模型的精确程度，且限制条件比较多；半解析法是在解析法的基础之上，从页岩气藏的基本渗流机理出发，建立准确的数学模型，并且通过点源函数理论得到点源解，然后考虑压裂形成的复杂裂缝形态，通过离散水力裂缝建立准确的产能模型，思路简洁直观，操作性强，且能够获得准确的预测结果。
综上所述，目前需要的页岩气藏压裂水平井产能预测方法应该具有以下特点：(1)对产量进行动态预测前期所需参数少，计算结果准确；(2)考虑页岩气藏的复杂渗流机理，在非稳态扩散方式下，对早期产量进行准确预测。
发明内容
本发明的目的在于提供一种页岩气藏压裂水平井非稳态扩散下产能计算方法，该方法原理可靠，操作简便，能够在给定页岩储层参数和水力裂缝参数的前提下，准确预测非稳态扩散方式下页岩气藏压裂水平井的产能。
为达到以上技术目的，本发明采用以下技术方案。
将页岩气藏看作是双重介质系统，主要包括页岩基质系统和天然裂缝系统，而页岩气在储层中的渗流又分为两个阶段，第一阶段为气体在页岩储层中解吸后通过非稳态扩散的方式运移到天然裂缝系统，第二阶段为天然裂缝中气体进入水力裂缝，通过水力裂缝流入井筒直至产出地面的过程。
一种页岩气藏压裂水平井非稳态扩散下产能计算方法，依次包括以下步骤：
(1)收集压裂改造所形成的水力裂缝参数，包括水力裂缝间距、缝长、方位角以及裂缝条数，计算页岩气藏弹性储容比ω、非稳态扩散窜流系数λ、无因次渗透率模量γD、页岩气吸附解吸系数σ和水力裂缝有限导流能力CFD；
(2)结合收集的水力裂缝参数，在直角坐标系中确定水力裂缝位置，离散水力裂缝并确定各裂缝微元的坐标，然后根据点源函数理论和压降叠加原理，求取每条水力裂缝缝尖处的压力；
(3)根据等效面积原则建立水力裂缝渗流模型，得到每条水力裂缝在水平井筒处的压力；
(4)计算在非稳态扩散下的页岩气藏压裂水平井的产能。
所述步骤(1)中计算的参数包括页岩气藏弹性储容比ω、非稳态扩散窜流系数λ、无因次渗透率模量γD、页岩气吸附解吸系数σ和水力裂缝有限导流能力CFD，过程如下：
1)页岩储层发育基质孔隙和天然裂缝，是页岩气在储层中的主要储存空间，为了描述基质和天然裂缝弹性储容能力的大小，定义页岩气藏的弹性储容比ω，通过下式计算(李晓平.地下油气渗流力学[M].北京:石油工业出版社,2007:153)：
式中
ki——天然裂缝系统的初始渗透率，m2；
μi——气体初始黏度，Pa·s；
φ——天然裂缝系统孔隙度，无因次；
Cgi——气体初始压缩系数，Pa-1；
h——页岩气藏的厚度，m；
qsc——页岩气井初始产量，m3/s。
2)页岩基质系统的气体通过拟稳态扩散的方式向天然裂缝系统运移，定义扩散窜流系数λ表征此阶段运移特征，通过下式计算(李晓平.地下油气渗流力学[M].北京:石油工业出版社,2007:154)：
式中
Lref——参考长度，本发明中选取水平井筒半径，m；
D——扩散系数，m2/s；
R——页岩基质球块的半径，m。
3)页岩气藏发育大量的天然裂缝，在储层压力降低时，会发生应力敏感效应，因此定义无因次渗透率模量γD来表征应力敏感现象，通过下式计算(Liu M,Xiao C,Wang Y,etal.Sensitivity analysis of geometry for mutli-stage fractured horizontalwells with consideration of finite-conductivity fractures in shale gasreservoirs[J].Journal of Natural Gas Science&Engineering,2015,22:182-195)：
式中
Z——气体偏差因子，无因次；
μ——气体黏度，Pa·s；
T——页岩储层的温度，K；
psc——标准状况下压力(我国石油行业规定取值为0.101325MPa)，Pa；
Tsc——标准状况下温度(我国石油行业规定标准状况温度为293K)，K；
γ——应力敏感系数，Pa-1；
p——天然裂缝系统的压力，Pa。
4)页岩气藏基质孔隙表面吸附有大量的吸附气，吸附气解吸后对页岩气藏的运移有很大影响，因此引入吸附解吸系数σ表征页岩气藏吸附能力的大小，计算如下：
式中——Langmuir拟压力，Pa/s；
——页岩气藏初始拟压力，Pa/s；
pL——页岩气的Langmuir压力，Pa；
VL——页岩气的Langmuir体积，sm3/m3；
pi——页岩气藏初始压力，Pa。
公式(4)～(6)中的Langmuir压力pL和Langmuir体积VL都是表征页岩气吸附解吸特性的参数(Langmuir I.The adsorption of gases on plane surfaces of glass,micaand platinum[J].Journal of Chemical Physics,2015,40(12):1361-1403)，而Langmuir拟压力是利用拟压力定义把pL转换为拟压力形式。
5)经过水力压裂形成的水力裂缝具有有限的导流能力，水力裂缝导流能力CFD通过下式计算：
式中kF——水力裂缝渗透率，m2；
wF——水力裂缝宽度，m；
CFD——水力裂缝导流能力，无因次。
所述步骤(2)过程如下：根据水力裂缝参数，在平面直角坐标系中绘制出水力裂缝的位置分布情况，把每条水力裂缝分割成若干微元，确定离散后每个水力裂缝微元及每条水力裂缝缝尖的坐标位置：
水力裂缝微元中心坐标
水力裂缝微元节点坐标
水力裂缝缝尖的坐标为水力裂缝两翼末端微元的节点坐标，即j＝1或2N。
式中——第i条水力裂缝的第j个微元的中心横坐标，m；
——第i条水力裂缝的第j个微元的中心纵坐标，m；
xi,j——第i条水力裂缝的第j个微元的节点横坐标，m；
yi,j——第i条水力裂缝的第j个微元的节点纵坐标，m；
yi——第i条水力裂缝与水平井筒相交时的纵坐标，m；
xfli——第i条水力裂缝的左翼裂缝长度，m；
xfri——第i条水力裂缝的右翼裂缝长度，m；
N——水力裂缝单翼离散微元数；
αi——第i条水力裂缝与y轴的夹角，度。
然后根据点源函数理论和压降叠加原理，通过拉氏变换得到第m条水力裂缝缝尖处的压力计算公式(孔祥言.高等渗流力学[M].合肥:中国科学技术大学出版社,1999:171-172；张建国,杜殿发,侯健,等.油气层渗流力学[M].东营:中国石油大学出版社,2009:79-84)：
式中
x,y——地层中某一点的横、纵坐标，m；
xi,j,xi,j+1——第i条水力裂缝上第j个裂缝微元节点的横坐标，m；
yi——第i条水力裂缝与y轴的交点坐标值，m；
ym——第m条水力裂缝与y轴的交点坐标值，m；
xw,yw——水力裂缝微元的中心横、纵坐标，m；
xwDi——积分变量，无因次；
xflm——第m条水力裂缝左翼裂缝长度，m；
RD——某一裂缝微元中心到地层某一点的距离，无因次；
——第i条裂缝的无因次流量；无因次；
s——拉普拉斯变量，无因次；
αm——第m条水力裂缝与y轴的夹角，°；
M——页岩气藏水平井压裂产生水力裂缝的条数；
K0——零阶第二类修正Bessel函数。
所述步骤(3)过程如下：
首先根据等效面积原则建立水力裂缝渗流模型，得到气体在水力裂缝中作平面径向流时的渗流公式(Mukherjee H,Economides M J.A parametric comparison ofhorizontal and vertical well performance[J].SPE Formation Evaluation,1991,6(2):209-216)：
πre 2＝(xfrm+xflm)h (15)
根据式(15)和(16)可得：
式中re——等效为平面径向流后的供给半径，m；
rw——水平井筒半径，m；
xflm——第m条水力裂缝的左翼缝长度，m；
xfrm——第m条水力裂缝的右翼缝长度，m；
pfm——第m条水力裂缝缝尖处的压力，Pa；
pw——水平井筒处压力，Pa；
qfm——第m条水力裂缝的流量；m3/s。
考虑水力裂缝的有限导流能力，通过拉氏变换和扰动变换的方法可得到拉氏空间下的水力裂缝渗流模型(孔祥言.高等渗流力学[M].合肥:中国科学技术大学出版社,1999:202-206；Pedrosa O A.Pressure Transient Response in Stress-SensitiveFormations[J].SPE California Regional Meeting,1986)，得到第m条水力裂缝在水平井筒处的压力计算公式：
式中
xfrm——第m条水力裂缝右翼裂缝长度，m；
——无因次形式的第m条水力裂缝尖产生的压力；
——无因次形式的第m条水力裂缝在水平井筒处的压力；
——第m条水力裂缝的无因次流量。
将第m条水力裂缝缝尖处的压力计算式代入，得到第m条水力裂缝在水平井筒处的压力表达式：
所述步骤(4)包括以下内容：
1)假设水平井筒中的渗流不存在阻力，即第m条水力裂缝在水平井筒处的压力等于井底压力步骤(3)水力裂缝在水平井筒处的压力计算式中含有第m条水力裂缝的无因次流量和井底压力然后结合流量约束条件，并令m＝1,2,3,...,M，可得到包含 的M+1个方程组合而成的方程组，将步骤(1)计算的参数代入，求解得到页岩气井定流量生产时的井底压力
2)根据拉氏空间下页岩气井定流量生产时的井底压力与定井底压力生产时的气井流量之间的转换关系，可得到拉氏空间下定井底压力生产时的页岩气井无因次流量然后通过数值反演和无因次流量定义，得到真实空间下的页岩气井流量q，即页岩气藏压裂水平井在非稳态扩散下的产能。
具体过程如下：
在式(19)中同时包含第m条水力裂缝的无因次流量和井底压力令m＝1,2,3,…,M，则可以得到包含和共M+1个未知量的M个方程；然后结合流量约束条件：
可以得到M+1个方程组成的方程组，用矩阵的形式表示为：
其中m＝1,2,3,...,M；i＝1,2,3,...,M。
通过高斯法(李庆扬,王能超,易大义.数值分析(第5版)[M].清华大学出版社,北京:2008)求解矩阵方程组(21)，可以得到页岩气井定流量生产时的井底压力
根据拉氏空间下页岩气井定流量生产时的井底压力与定井底压力生产时的气井流量之间的转换关系(Van Everdingen AF,Hurst W.The application of theLaplace transformation to flow problems in reservoirs[J].Journal of PetroleumTechnology,1949,1(12):305-324)，得到拉氏空间下定井底压力生产时的页岩气井无因次流量
通过数值反演(同登科,陈钦雷.关于Laplace数值反演Stehfest方法的一点注记[J].石油学报,2001,22(6):91-92)得到真实空间下的页岩气井无因次流量qD，最后利用无因次流量定义计算页岩气井流量q，即页岩气藏压裂水平井在非稳态扩散下的产能：
q＝qD×qsc (26)
本发明所涉及的公式较多，若一一列出会显得累赘，为了保证本发明的简洁直观性，在此只列出了计算产量时必要的计算公式和参数，对于未给出的计算公式和参数则写出了相应的参考文献。
与现有技术相比，本发明的有益效果为：
(1)该方法考虑页岩气藏的实际运移机理，通过非稳态扩散的方式流入天然裂缝，并且考虑天然裂缝的应力敏感效应，利用面积等效原则建立考虑有限导流能力的水力裂缝渗流模型，采用点源函数方法和压降叠加原理求出水力裂缝在水平井筒处的压力表达式，然后通过求解矩阵方程组得到拉氏空间下的井底压力，最后通过转换关系得到页岩气井的产量。该方法求解产量思路清晰，操作简便。
(2)该方法建立的模型考虑因素比较全面，求解速度快，前期不需要大量的储层参数，并且计算准确度高。
附图说明
图1为页岩气藏多级压裂水平井裂缝离散模型示意图。
图2为某页岩气井生产600天的实际产量递减曲线和计算得到产量曲线的对比图。
具体实施方式
以下结合附图和现场实际数据对本发明作进一步说明。
以某口页岩气井为例，该井所在的页岩储层厚度为30m，页岩储层平均孔隙度为0.08，页岩基质渗透率为2.5×10-9mD，天然裂缝系统的渗透率为0.015mD，可见该页岩气藏为超低渗、低孔储层，需要进行压裂改造才能产生工业气流，页岩气藏的具体参数如表1所示。
表1某口页岩气井基本参数表
一种页岩气藏压裂水平井非稳态扩散下产能计算方法，依次包括以下步骤：
步骤1，计算中间参数弹性储容比ω、非稳态扩散窜流系数λ、无因次渗透率模量γD、页岩气吸附解吸系数σ、裂缝导流能力CFD，方便后面方程求解；
步骤2，根据收集的水力裂缝参数，在平面直角坐标系中绘制出水力裂缝的位置分布情况，把每条水力裂缝分割成若干微元(图1所示)，然后根据式(8)～(11)确定离散后每个水力裂缝微元及每条水力裂缝缝尖的坐标位置，再得到所有水力裂缝生产时对各条水力裂缝尖端所产生的压力表达式；
步骤3，根据每条水力裂缝缝尖处的压力，得到水力裂缝在水平井筒处的压力表达式；
步骤4，将式(12)和式(18)联立可得到页岩气藏所有水力裂缝在水平井筒处的压力，结合流量约束条件，通过高斯法求解矩阵方程组式(21)，可得到页岩气井定流量生产时的井底压力然后根据式(25)计算得到定井底压力条件下页岩气井拉氏空间下的无因次流量，然后采用数值反演方法得到真实空间下页岩气井无因次流量，最后利用无因次流量定义式(26)，可得到页岩气井流量。
从图2可以看出，利用本方法计算得到的页岩气井产量能够很好地吻合实际生产曲线，误差比较小，说明本发明提供的页岩气藏压裂水平井非稳态扩散下的产能计算方法较为合理，而且还可以直观地预测页岩气井后期产量递减趋势，指导压裂水平井参数优化设计以及评估压裂改造有效性。
Claims (2)
1.一种页岩气藏压裂水平井非稳态扩散下产能计算方法，依次包括以下步骤：
(1)收集压裂改造所形成的水力裂缝参数，包括水力裂缝间距、缝长、方位角以及裂缝条数，计算页岩气藏弹性储容比ω、非稳态扩散窜流系数λ、无因次渗透率模量γD、页岩气吸附解吸系数σ和水力裂缝有限导流能力CFD；
(2)在直角坐标系中确定水力裂缝位置，离散水力裂缝并确定各裂缝微元的坐标，根据点源函数理论和压降叠加原理，求取每条水力裂缝缝尖处的压力，过程如下：根据水力裂缝参数，在平面直角坐标系中绘制出水力裂缝的位置分布情况，把每条水力裂缝分割成若干微元，确定离散后每个水力裂缝微元及每条水力裂缝缝尖的坐标位置：
水力裂缝微元中心坐标
水力裂缝微元节点坐标
水力裂缝缝尖的坐标为水力裂缝两翼末端微元的节点坐标，即j＝1或2N；
式中——第i条水力裂缝的第j个微元的中心横坐标，m；
——第i条水力裂缝的第j个微元的中心纵坐标，m；
xi,j——第i条水力裂缝的第j个微元的节点横坐标，m；
yi,j——第i条水力裂缝的第j个微元的节点纵坐标，m；
yi——第i条水力裂缝与水平井筒相交时的纵坐标，m；
xfli——第i条水力裂缝的左翼裂缝长度，m；
xfri——第i条水力裂缝的右翼裂缝长度，m；
N——水力裂缝单翼离散微元数；
αi——第i条水力裂缝与y轴的夹角，度；
然后根据点源函数理论和压降叠加原理，通过拉氏变换得到第m条水力裂缝缝尖处的压力计算公式：
式中
x,y——地层中某一点的横、纵坐标，m；
xi,j,xi,j+1——第i条水力裂缝上第j个裂缝微元节点的横坐标，m；
yi——第i条水力裂缝与y轴的交点坐标值，m；
ym——第m条水力裂缝与y轴的交点坐标值，m；
xw,yw——水力裂缝微元的中心横、纵坐标，m；
xwDi——积分变量，无因次；
xflm——第m条水力裂缝左翼裂缝长度，m；
RD——某一裂缝微元中心到地层某一点的距离，无因次；
——第i条裂缝的无因次流量；无因次；
s——拉普拉斯变量，无因次；
αm——第m条水力裂缝与y轴的夹角，°；
M——页岩气藏水平井压裂产生水力裂缝的条数；
K0——零阶第二类修正Bessel函数；
Lref——参考长度，选取水平井筒半径，m；
(3)根据等效面积原则建立水力裂缝渗流模型，得到每条水力裂缝在水平井筒处的压力，过程如下：
首先根据等效面积原则建立水力裂缝渗流模型，得到气体在水力裂缝中作平面径向流时的渗流公式：
式中re——等效为平面径向流后的供给半径，m；
rw——水平井筒半径，m；
xflm——第m条水力裂缝的左翼缝长度，m；
xfrm——第m条水力裂缝的右翼缝长度，m；
pfm——第m条水力裂缝缝尖处的压力，Pa；
pw——水平井筒处压力，Pa；
qfm——第m条水力裂缝的流量；m3/s；
h——页岩气藏的厚度，m；
μ——气体黏度，Pa·s；
kF——水力裂缝渗透率，m2；
wF——水力裂缝宽度，m；
考虑水力裂缝的有限导流能力，得到拉氏空间下的水力裂缝渗流模型，得到第m条水力裂缝在水平井筒处的压力计算公式：
式中
xfrm——第m条水力裂缝右翼裂缝长度，m；
——无因次形式的第m条水力裂缝尖产生的压力；
——无因次形式的第m条水力裂缝在水平井筒处的压力；
——第m条水力裂缝的无因次流量；
再将第m条水力裂缝缝尖处的压力计算式代入，得到第m条水力裂缝在水平井筒处的压力计算式：
(4)计算在非稳态扩散下的页岩气藏压裂水平井的产能，包括以下内容：
1)假设水平井筒中的渗流不存在阻力，第m条水力裂缝在水平井筒处的压力等于井底压力水力裂缝在水平井筒处的压力计算式中含有第m条水力裂缝的无因次流量和井底压力结合流量约束条件，令m＝1,2,3,...,M，得到包含 的M+1个方程组合而成的方程组，求得页岩气井定流量生产时的井底压力
2)根据拉氏空间下页岩气井定流量生产时的井底压力与气井流量之间的转换关系，得到拉氏空间下页岩气井无因次流量从而得到真实空间下的页岩气井流量q，即页岩气藏压裂水平井在非稳态扩散下的产能。
2.如权利要求1所述的一种页岩气藏压裂水平井非稳态扩散下产能计算方法，其特征在于，所述结合流量约束条件为：
根据页岩气井定流量生产时的井底压力通过下式得到拉氏空间下页岩气井无因次流量
通过数值反演得到真实空间下的页岩气井无因次流量qD，再通过下式计算页岩气井流量q，即页岩气藏压裂水平井在非稳态扩散下的产能：
q＝qD×qsc
式中qsc——页岩气井初始产量，m3/s。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910498213.4A CN110210157B (zh) | 2019-06-10 | 2019-06-10 | 一种页岩气藏压裂水平井非稳态扩散下产能计算方法 |
| RU2019131430A RU2732868C1 (ru) | 2019-06-10 | 2019-10-04 | Способ расчета продуктивности горизонтальных скважин в залежах сланцевого газа при гидроразрыве в условиях нестационарной диффузии |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910498213.4A CN110210157B (zh) | 2019-06-10 | 2019-06-10 | 一种页岩气藏压裂水平井非稳态扩散下产能计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110210157A CN110210157A (zh) | 2019-09-06 |
| CN110210157B true CN110210157B (zh) | 2019-12-20 |
Family
ID=67791642
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910498213.4A Active CN110210157B (zh) | 2019-06-10 | 2019-06-10 | 一种页岩气藏压裂水平井非稳态扩散下产能计算方法 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN110210157B (zh) |
| RU (1) | RU2732868C1 (zh) |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110805436B (zh) * | 2019-10-09 | 2021-12-03 | 中国石油大学(北京) | 基于压降数据的单段裂缝产液贡献率评价方法及设备 |
| CN111307684A (zh) * | 2019-12-04 | 2020-06-19 | 西南石油大学 | 一种用于计算微纳米孔隙中气体渗透率的分子模拟方法 |
| CN112051182B (zh) * | 2020-09-07 | 2022-01-28 | 西南石油大学 | 一种快速预测不同深度页岩储层吸附甲烷能力的方法 |
| CN111911148B (zh) * | 2020-09-16 | 2021-06-29 | 西南石油大学 | 一种模拟页岩天然裂缝中气液两相渗流特征的方法 |
| CN112084718B (zh) * | 2020-09-16 | 2021-05-04 | 西南石油大学 | 基于渗流差异的页岩气藏单相气体三孔三渗模型构建方法 |
| CN111980695B (zh) * | 2020-09-16 | 2021-03-12 | 西南石油大学 | 基于不同占比的页岩气藏单相气体双孔双渗模型构建方法 |
| CN112082922B (zh) * | 2020-09-18 | 2021-03-16 | 西南石油大学 | 一种矩形平板大模型岩样平面渗流渗透率的确定方法 |
| CN112360422B (zh) * | 2020-12-08 | 2022-04-08 | 西南石油大学 | 一种页岩气藏压裂水平井产量预测方法及系统 |
| CN112878982B (zh) * | 2020-12-31 | 2022-03-01 | 西南石油大学 | 一种考虑裂缝长期导流能力的深层页岩气产能预测方法 |
| CN112924354B (zh) * | 2021-01-22 | 2022-09-20 | 西南石油大学 | 一种基于气体扩散系数的页岩储层应力敏感实验评价方法 |
| CN113029898B (zh) * | 2021-02-22 | 2022-04-15 | 西南石油大学 | 一种裂缝动态导流能力、基岩供气能力的测试装置及方法 |
| CN114065113B (zh) * | 2021-11-19 | 2024-02-27 | 西北大学 | 一种自然裂缝储层多段压裂水平井干扰试井分析方法 |
| CN114166698B (zh) * | 2022-02-08 | 2022-04-26 | 中国矿业大学(北京) | 一种竞争吸附下煤微孔道气体扩散能力评估方法和系统 |
| CN115270663B (zh) * | 2022-09-27 | 2022-12-09 | 中国石油大学(华东) | 一种描述页岩基质与天然裂缝间非稳态窜流的方法及系统 |
| CN115906428B (zh) * | 2022-11-02 | 2023-07-28 | 北京科技大学 | 生产动态数据的页岩凝析气藏凝析油产量预测方法及装置 |
| CN115994500B (zh) * | 2023-03-22 | 2023-06-02 | 北京科技大学 | 评估页岩气井压裂裂缝导流能力动态变化的方法及系统 |
| CN117634347A (zh) * | 2023-11-29 | 2024-03-01 | 北京科技大学 | 一种判断多段压裂水平井井间干扰强度的方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106351651A (zh) * | 2016-08-26 | 2017-01-25 | 中国石油天然气股份有限公司 | 页岩气井产能的预测方法及装置 |
| CN107066674A (zh) * | 2017-01-18 | 2017-08-18 | 西南石油大学 | 计算页岩气藏体积压裂水平井非稳态产量的方法 |
| CN108240214A (zh) * | 2018-01-18 | 2018-07-03 | 西南石油大学 | 一种页岩气藏水平井分段压裂压后产能预测方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7272973B2 (en) * | 2005-10-07 | 2007-09-25 | Halliburton Energy Services, Inc. | Methods and systems for determining reservoir properties of subterranean formations |
-
2019
- 2019-06-10 CN CN201910498213.4A patent/CN110210157B/zh active Active
- 2019-10-04 RU RU2019131430A patent/RU2732868C1/ru active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106351651A (zh) * | 2016-08-26 | 2017-01-25 | 中国石油天然气股份有限公司 | 页岩气井产能的预测方法及装置 |
| CN107066674A (zh) * | 2017-01-18 | 2017-08-18 | 西南石油大学 | 计算页岩气藏体积压裂水平井非稳态产量的方法 |
| CN108240214A (zh) * | 2018-01-18 | 2018-07-03 | 西南石油大学 | 一种页岩气藏水平井分段压裂压后产能预测方法 |
Non-Patent Citations (4)
| Title |
|---|
| Production forecast of a multistage fractured horizontal well by an analytical method in shale gas reservoir;Yongming Li等;《Environmental Earth Sciences》;20190423;第78卷;271-290 * |
| 页岩气水平井增产改造体积评价模型及其应用;任岚等;《开发工程》;20180831;第38卷(第8期);47-56 * |
| 页岩气水平井多段压裂缝—井筒流动耦合模型研究;张月娟;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20171115;B019-238 * |
| 页岩气藏压裂水平井压力及产量动态特征研究;程子洋;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20150915;B019-103 * |
Also Published As
| Publication number | Publication date |
|---|---|
| RU2732868C1 (ru) | 2020-09-24 |
| CN110210157A (zh) | 2019-09-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110210157B (zh) | 一种页岩气藏压裂水平井非稳态扩散下产能计算方法 | |
| CN109441422B (zh) | 一种页岩气井间距优化开采方法 | |
| Zhao et al. | Latest research progresses in network fracturing theories and technologies for shale gas reservoirs | |
| Wu et al. | A model for multiple transport mechanisms through nanopores of shale gas reservoirs with real gas effect–adsorption-mechanic coupling | |
| Cao et al. | Combined impact of flow regimes and effective stress on the evolution of shale apparent permeability | |
| Feng et al. | Apparent permeability model for shale oil with multiple mechanisms | |
| Zhang et al. | Employing a quad-porosity numerical model to analyze the productivity of shale gas reservoir | |
| Li et al. | Composition-transient analysis in shale-gas reservoirs with consideration of multicomponent adsorption | |
| CN104612635A (zh) | 煤层群瓦斯联合抽采达标预评判方法 | |
| Liu et al. | Numerical modeling of gas flow in coal using a modified dual-porosity model: a multi-mechanistic approach and finite difference method | |
| Ding et al. | The analytical model of hydraulic fracture initiation for perforated borehole in fractured formation | |
| Li et al. | A new method for production data analysis in shale gas reservoirs | |
| Gao et al. | Effect of shale matrix heterogeneity on gas transport during production: A microscopic investigation | |
| Xu et al. | Production performance of horizontal wells with dendritic-like hydraulic fractures in tight gas reservoirs | |
| Liu et al. | Well type and pattern optimization method based on fine numerical simulation in coal-bed methane reservoir | |
| Wang et al. | Experimental investigation of reservoir fluid interlayer crossflow through fracture during the drainage stage of coal measure gas well | |
| CN114201932A (zh) | 一种复杂情况下的致密油藏压裂井试井模拟方法 | |
| Zhang et al. | Gas adsorptions of geological carbon storage with enhanced gas recovery | |
| Jiang et al. | Characterization of wormhole growth and propagation dynamics during cold heavy oil production with sand (CHOPS) processes by integrating rate transient analysis and a pressure-gradient-based sand failure criterion | |
| Wang et al. | Optimization of rock failure criteria under different fault mechanisms and borehole trajectories | |
| Lei et al. | Characterization of multiphase flow in shale oil reservoirs considering multiscale porous media by high-resolution numerical simulation | |
| Yan et al. | A coupled matrix-fracture productivity calculation model considering low-velocity non-Darcy flow in shale reservoirs | |
| Liu et al. | Study on effective extraction radius of directional long borehole and analysis of the influence mechanism | |
| Shi et al. | A coupled nonlinear flow model for particle migration and seepage properties of water inrush through broken rock mass | |
| Rongze et al. | A review of empirical production decline analysis methods for shale gas reservoir |
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 |