一种页岩气藏压裂水平井非稳态扩散下产能计算方法技术

本发明专利技术公开了一种页岩气藏压裂水平井非稳态扩散下产能计算方法，包括：（1）收集压裂改造所形成的水力裂缝参数，计算页岩气藏弹性储容比、非稳态扩散窜流系数、无因次渗透率模量、页岩气吸附解吸系数和水力裂缝有限导流能力；（2）在直角坐标系中确定水力裂缝位置，离散水力裂缝并确定各裂缝微元的坐标，根据点源函数理论和压降叠加原理，求取每条水力裂缝缝尖处的压力；（3）根据等效面积原则建立水力裂缝渗流模型，得到每条水力裂缝在水平井筒处的压力；（4）计算在非稳态扩散下的页岩气藏压裂水平井的产能。本发明专利技术原理可靠，操作简便，能够在给定页岩储层参数和水力裂缝参数的前提下，准确预测非稳态扩散方式下页岩气藏压裂水平井的产能。

A Method for Calculating Productivity of Fractured Horizontal Wells in Shale Gas Reservoirs under Unsteady Diffusion

【技术实现步骤摘要】
一种页岩气藏压裂水平井非稳态扩散下产能计算方法
本专利技术涉及油气田开发领域页岩气藏压裂水平井的产能计算方法，特别是涉及一种非稳态扩散方式下页岩气藏压裂水平井的产能计算方法。
技术介绍
随着我国能源需求矛盾日益增强，页岩气作为非常规油气资源的有效补充，备受关注，也成为我国天然气增储上产的关键一环。页岩气在储层中赋存方式多样，运移机理复杂，通常包括吸附解吸、扩散、渗流等多重运移机制，而常规的渗流理论也不适合于页岩气藏，正确且清楚地认识页岩气藏的渗流规律及正确表征页岩气运移机制对指导页岩气藏开发具有重要的意义。基质中的页岩气通过扩散的方式运移到天然裂缝系统，扩散的方式分为拟稳态和非稳态扩散，而实际气藏的开发过程中几乎不会存在拟稳态的扩散方式，气体从基质向天然裂缝运移更符合非稳态扩散方式。由于页岩的低孔、超低渗等特征，水平井分段多簇压裂方法成为了常用的增产措施，而压裂后页岩气藏水平井产能的准确预测是评估压裂改造有效性及页岩气藏是否成功开发的关键，通常建立页岩气藏产能模型的方法有解析法、数值模拟法和半解析法，数值模拟法主要采用计算机程序来近似求解，初期所需气藏参数多，求解迭代时间长，对硬件要求高等特点限制应用；解析法对实际情况进行简化，通过高度抽象、理想化物理模型，利用各种数学方法对模型进行直接求解，产能预测的可靠性主要取决于物理模型的精确程度，且限制条件比较多；半解析法是在解析法的基础之上，从页岩气藏的基本渗流机理出发，建立准确的数学模型，并且通过点源函数理论得到点源解，然后考虑压裂形成的复杂裂缝形态，通过离散水力裂缝建立准确的产能模型，思路简洁直观，操作性强，且能够获得准确的预测结果。综上所述，目前需要的页岩气藏压裂水平井产能预测方法应该具有以下特点：(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来表征应力敏感现象，通过下式计算(LiuM,XiaoC,WangY,etal.Sensitivityanalysisofgeometryformutli-stagefracturedhorizontalwellswithconsiderationoffinite-conductivityfracturesinshalegasreservoirs[J].JournalofNaturalGasScience&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都是表征页岩气吸附解吸特性的参数(LangmuirI.Theadsorptionofgasesonplanesurfacesofglass,micaandplatinum[J].JournalofChemicalPhysics,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条水力裂缝的右翼裂缝本文档来自技高网...

【技术保护点】
1.一种页岩气藏压裂水平井非稳态扩散下产能计算方法，依次包括以下步骤：(1)收集压裂改造所形成的水力裂缝参数，包括水力裂缝间距、缝长、方位角以及裂缝条数，计算页岩气藏弹性储容比ω、非稳态扩散窜流系数λ、无因次渗透率模量γD、页岩气吸附解吸系数σ和水力裂缝有限导流能力CFD；(2)在直角坐标系中确定水力裂缝位置，离散水力裂缝并确定各裂缝微元的坐标，根据点源函数理论和压降叠加原理，求取每条水力裂缝缝尖处的压力；(3)根据等效面积原则建立水力裂缝渗流模型，得到每条水力裂缝在水平井筒处的压力；(4)计算在非稳态扩散下的页岩气藏压裂水平井的产能。

【技术特征摘要】
1.一种页岩气藏压裂水平井非稳态扩散下产能计算方法，依次包括以下步骤：(1)收集压裂改造所形成的水力裂缝参数，包括水力裂缝间距、缝长、方位角以及裂缝条数，计算页岩气藏弹性储容比ω、非稳态扩散窜流系数λ、无因次渗透率模量γD、页岩气吸附解吸系数σ和水力裂缝有限导流能力CFD；(2)在直角坐标系中确定水力裂缝位置，离散水力裂缝并确定各裂缝微元的坐标，根据点源函数理论和压降叠加原理，求取每条水力裂缝缝尖处的压力；(3)根据等效面积原则建立水力裂缝渗流模型，得到每条水力裂缝在水平井筒处的压力；(4)计算在非稳态扩散下的页岩气藏压裂水平井的产能。2.如权利要求1所述的一种页岩气藏压裂水平井非稳态扩散下产能计算方法，其特征在于，所述步骤(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...

【专利技术属性】
技术研发人员：李勇明，吴磊，张烈辉，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：四川,51