引言

对于低渗透气藏，仅采用单一水平井或压裂垂直井开发往往达不到理想的单井产能和预期的开发效果。水平井多级分段压裂技术能大幅度提高低渗气藏的单井产能，因此，近年来被越来越多地运用于低渗气藏、致密气藏、页岩气藏的开发，并取得了良好的开发效益。

1978年，Cinco-Ley等建立了压裂垂直井的渗流模型，并采用裂缝单元离散的方法对模型进行了求解^[1]，该模型能准确地刻画压裂垂直井的压力和产量动态，但模型仅考虑了单条压裂缝的影响。1993年，Raghavan和Joshi给出了一种多级压裂水平井产能的计算方法，但该模型未严格考虑压裂缝中的流动^[2]。1994年，郎兆新利用位势理论和叠加原理对油藏中的多级压裂水平井(MFHW)产能进行了研究，并推导出了产能计算公式，不过该公式未考虑压裂缝导流能力的影响，只适用于无限导流缝的情形^[3]。2002年，宁正福等在郎兆新研究基础上对MFHW产能计算公式进行了改进，该公式进一步考虑了压裂缝导流能力的影响，能较好地用于MFHW产能的计算，不过，该模型对流体在压裂缝中的流动进行了较大的简化处理，与实际情况偏差较大^[4]。2007年，曾凡辉等利用点汇解和叠加原理为MFHW建立了一种产能计算模型，该模型既考虑了压裂缝引起的势的相互干扰，也考虑了压裂缝导流能力的影响^[5]，不过，采用个数有限的点汇的叠加来代替连续延伸压裂缝的方法不可避免地导致较大的偏差。2009年，王海涛等建立了多级压裂水平井不稳定渗流模型，计算并分析了MFHW井底压力动态及流量分布规律^[6]，但该模型未考虑压裂缝导流能力的影响。同年，Brown将压裂水平井分为3个线状流动区，建立了三线性流模型，分析了MFHW的压力动态，并利用Mukherjee和Economides提出的表皮系数表达式考虑了压裂缝中的流体向水平井筒径向汇聚的影响，但文中未对产量动态进行分析^[7-8]。2010年，曾保全等借助Eclipse软件，利用流线模拟方法，研究了压裂水平井流线分布特征、裂缝产能分布规律、裂缝长度对开发效果的影响^[9]。2014年，张芮菡等基于压裂水平井三线性流模型，结合沃伦-茹特模型，建立了低渗透裂缝性气藏压裂水平井不稳定渗流数学模型，分析了其产量变化动态^[10]，该模型较好地考虑了流体在地层及压裂缝中的流动，但仅适合于早中期的线性流阶段，对晚期流线形态发生变化后则不太适用。2016年，张德良基于非结构网格，采用有限元-有限体积法研究了页岩气藏MFHW水平井产能动态^[11]，模型功能较强，能考虑很多复杂情况，但计算量大、计算耗时。同年，魏明强、段永刚等基于PEBI网格，引入尘气模型建立并推导了综合考虑页多重运移机制下的压裂水平井产量递减数学模型，计算获得了页岩气藏无限导流压裂水平井Blasingame产量递减典型曲线，讨论了相关参数对典型曲线的影响^[12]，但未考虑压裂缝导流能力的影响。此外，还有不少学者对多级压裂水平井产能计算问题作了研究^[13-21]。

由此可见，目前低渗气藏MFHW产能计算方面仍存在一些不足，为此，本文综合利用多种数学方法，建立了一个严格的产能计算数学模型。

1 物理模型

如图 1所示，低渗气藏中有一口多级压裂水平井(MFHW)，图中只画了6条压裂缝，但为了不失一般性，在推导时假设压裂缝条数为$M$，压裂缝半长为$x_{\rm f}$，水平井长度为$L_{\rm h}$。

其他基本假设如下：

(1) 气藏水平等厚，原始地层压力为$p_{\rm i}$；

(2) 考虑水力压裂缝的导流能力；

(3) 忽略重力和毛细管力；

(4) 气藏中为等温渗流。

2 产能计算的数学模型

该产能计算模型由气藏渗流模型和压裂缝渗流模型组成。为了在推导中不产生太多冗长的数值常数，本文模型在推导时采用SI国际单位制。

2.1 气藏渗流模型

由于多级压裂水平井内边界条件的复杂性，直接构建其渗流数学模型时很难准确给出内边界条件，故先从最简单的线汇模型出发，再利用叠加原理获得MFHW在气藏中引起的压力分布表达式。

2.1.1 线汇模型

线汇模型是渗流力学中最基本的渗流模型之一，此处考虑以恒定的井底压力$p_{\rm w}$生产，则无因次线汇模型如下^[22]

渗流微分方程

$ \dfrac{{{\partial ^2}{\psi _{\rm{D}}}}}{{\partial {r_{\rm{D}}}^2}} + \dfrac{1}{{{r_{\rm{D}}}}}\dfrac{{\partial {\psi _{\rm{D}}}}}{{\partial {r_{\rm{D}}}}} = \dfrac{{\partial {\psi _{\rm{D}}}}}{{\partial {t_{\rm{D}}}}} $

(1)

内边界条件(线汇)

$ \mathop {{\rm{lim}}}\limits_{{\xi _{\rm{D}}} \to 0} {\left. {{r_{\rm{D}}}\dfrac{{\partial {\psi _{\rm{D}}}}}{{\partial {r_{\rm{D}}}}}} \right|_{{r_{\rm{D}}} = {\xi _{\rm{D}}}}} = - {\hat q_{\rm{D}}} $

(2)

外边界条件

$ {\left. {{\psi _{\rm{D}}}} \right|_{{r_{\rm{D}}} \to \infty }} = 0 $

(3)

初始条件

$ {\left. {{\psi _{\rm{D}}}} \right|_{{t_{\rm{D}}} = 0}} = 0 $

(4)

式中：$\psi _{\rm{D}}$—地层无因次拟压力，${\psi _{\rm{D}}} = \dfrac{{{\psi _{\rm{i}}} - \psi }}{{{\psi _{\rm{i}}} - {\psi _{\rm{w}}}}}$；

${r_{\rm{D}}}$—无因次径向距离，${r_{\rm{D}}} = \dfrac{r}{{{r_{\rm{w}}}}}$；

${t_{\rm{D}}}$—无因次时间，${t_{\rm{D}}} = \dfrac{{Kt}}{{\phi {\mu _{\rm{i}}}{C_{\rm{i}}}{r_{\rm{w}}}^2}}$；

${\hat q_{\rm{D}}}$—线汇无因次产量，

${\hat q_{\rm{D}}}({t_{\rm{D}}}) = \dfrac{{{p_{{\rm{sc}}}}T}}{{{\mathsf{π}}Kh{T_{{\rm{sc}}}}\left( {{\psi _{\rm{i}}} - {\psi _{\rm{w}}}} \right)}}\hat q(t)$；

${\psi _{\rm{i}}}$—原始地层压力$p_{\rm{i}}$对应的拟压力，Pa/s，

${\psi _{\rm{i}}} = \int\limits_{{p_0}}^{{p_{\rm{i}}}} {\dfrac{{2p}}{{\mu Z}}} {\rm{d}}p$；

${\psi _{\rm{w}}}$—井底压力$p_{\rm{w}}$对应的拟压力，Pa/s，

${\psi _{\rm{w}}} = \int\limits_{{p_0}}^{{p_{\rm{w}}}} {\dfrac{{2p}}{{\mu Z}}} {\rm{d}}p$；

$\psi$—任一压力$p$对应的拟压力，Pa/s，

$\psi = \int\limits_{{p_0}}^p {\dfrac{{2p}}{{\mu Z}}} {\rm{d}}p$；

$p_0$—参考压力，常取为0或一个大气压，Pa；

$r$—径向距离，m；

$r_{\rm w}$—井筒半径，m；

$K$—气藏渗透率，D；

$\phi$—气藏孔隙度，无因次；

$\mu_{\rm i}$—原始条件下的天然气黏度，Pa$\cdot$s；

$C_{\rm i}$—原始条件下天然气的压缩系数，Pa$^{-1}$；

$p_{\rm {sc}}$—地面大气压力，Pa；

$T$—气藏温度，K；

$h$—气藏厚度，m；

$T_{\rm{sc}}$—地面温度，K；

$\hat q(t)$—线汇地面产量，m$^3$/s。

式(1)$\sim$式(4)构成了最基本的定井底压力生产下的线汇模型。

利用Laplace变换可求得上述线汇模型在Laplace空间的解为

$ \overline {{\psi _{\rm{D}}}} = \overline {{{\hat q}_{\rm{D}}}} {{\rm{K}}_0}\left( {\sqrt s {r_{\rm{D}}}} \right) $

(5)

式中：

K$_0$—0阶虚宗量贝塞尔函数；

$s$—Laplace变量。

需要说明的是，本文描述中，变量上面带一杠(如$\overline {\psi _{\rm{D}}}$)表示该变量基于无因次时间${t_{\rm{D}}}$的Laplace变换。

2.1.2 MFHW在气藏中引起的无因次拟压降的表达

由求和形式和积分形式这两种形式的叠加原理，可得$M$条裂缝在气藏中引起的无因次拟压降为

$ \overline {{\psi _{\rm{D}}}} ({x_{\rm{D}}}, {y_{\rm{D}}}, s) = \dfrac{1}{2}\sum\limits_{l = 1}^M {\int\limits_{ - {x_{{\rm{fD}}l}}}^{{x_{{\rm{fD}}l}}} {\overline {{q_{{\rm{fD}}l}}} ({x_{{\rm{wD}}}}, s){{\rm{K}}_0}\left( {\sqrt s \sqrt {{{({x_{\rm{D}}} - {x_{{\rm{wD}}}})}^2} + {{({y_{\rm{D}}} - {y_{{\rm{wD}}l}})}^2}} } \right){\rm{d}}{x_{{\rm{wD}}}}} } $

(6)

式中：${q_{{\rm{fD}}l}}$—第$l$条压裂缝的无因次线密度流量，${q_{{\rm{fD}}l}} = \dfrac{{2{q_{{\rm{f}}l}}{r_{\rm{w}}}{p_{{\rm{sc}}}}T}}{{{\mathsf{π}}Kh{T_{{\rm{sc}}}}({\psi _{\rm{i}}} - {\psi _{\rm{w}}})}}$；

${q_{{\rm{f}}l}}$—第$l$条压裂缝的线密度流量，m$^3$/(s$\cdot$m)；

${x_{\rm{D}}}$—气藏中任一点的无因次$x$坐标，${x_{\rm{D}}} = \dfrac{x}{{{r_{\rm{w}}}}}$；

${y_{\rm{D}}}$—气藏中任一点的无因次$y$坐标，${y_{\rm{D}}} = \dfrac{y}{{{r_{\rm{w}}}}}$；

${x_{{\rm{wD}}}}$—压裂缝的无因次$x$坐标，${x_{{\rm{wD}}}} = \dfrac{{{x_{\rm{w}}}}}{{{r_{\rm{w}}}}}$；

${y_{{\rm{wD}}l}}$—第$l$条压裂缝的无因次$y$坐标，${y_{{\rm{wD}}l}} = \dfrac{{{y_{{\rm{w}}l}}}}{{{r_{\rm{w}}}}}$，$y_{{\rm{w}}l}$—第$l$条压裂缝的$y$坐标；

${x_{{\rm{fD}}l}}$—第$l$条缝的无因次缝半长，${x_{{\rm{fD}}l}} = \dfrac{{{x_{{\rm{f}}l}}}}{{{r_{\rm{w}}}}}$；

${x_{{\rm{f}}l}}$—第$l$条缝的缝半长，m。

2.2 多压裂缝渗流模型

沿用Cinco等建立双翼压裂缝渗流模型时的建模思想^[23-25]，可建立第$j$($j$=1，2，...$M$)条压裂缝的渗流模型如下

$ \dfrac{{{\partial ^2}{\psi _{{\rm{fD}}j}}}}{{\partial {x_{\rm{D}}}^2}} + \dfrac{2}{{{C_{{\rm{FD}}j}}}}{\left. {\dfrac{{\partial {\psi _{\rm{D}}}}}{{\partial {y_{\rm{D}}}}}} \right|_{{y_{\rm{D}}} = {y_{{\rm{wD}}j}} + {W_{{\rm{fD}}j}}/2}} = 0, \\[8pt] {\kern 2em} 0 \leqslant {x_{\rm{D}}} \leqslant {X_{{\rm{fD}}j}}, \;\;{y_{{\rm{wD}}j}} \leqslant {y_{\rm{D}}} \leqslant {y_{{\rm{wD}}j}} + {W_{{\rm{fD}}j}}/2 $

(7)

$ {\left. {\dfrac{{\partial {\psi _{{\rm{fD}}j}}}}{{\partial {x_{\rm{D}}}}}} \right|_{{x_{\rm{D}}} = 0}} = - \dfrac{{\mathsf{π}}}{{{C_{{\rm{FD}}j}}}}{q_{{\rm{D}}j}}({x_{\rm{D}}}, {t_{\rm{D}}}) $

(8)

$ {q_{{\rm{fD}}j}}({x_{\rm{D}}}, {t_{\rm{D}}}) = - \dfrac{2}{{\mathsf{π}}}{\left. {\dfrac{{\partial {\psi _{\rm{D}}}}}{{\partial {y_{\rm{D}}}}}} \right|_{{y_{\rm{D}}} = {y_{{\rm{wD}}j}} + {W_{{\rm{fD}}j}}/2}} $

(9)

式中：

${\psi _{{\rm{fD}}j}}$—第$j$条压裂缝的无因次拟压力；

${C_{{\rm{FD}}j}}$—第$j$条压裂缝的无因次导压系数，

${C_{{\rm{FD}}j}} =\dfrac{{{K_{{\rm{f}}j}}{W_{{\rm{f}}j}}}}{{K{r_{\rm{w}}}}}$；

${K_{{\rm{f}}j}}$—第$j$条压裂缝的渗透率，D；

${W_{{\rm{fD}}j}}$—第$j$条压裂缝的无因次裂缝宽度，

${W_{{\rm{fD}}j}} = {W_{{\rm{f}}j}}/{r_{\rm{w}}}$；

${W_{{\rm{f}}j}}$—第$j$条压裂缝的裂缝宽度，m。

式(7)$\sim$式(9)就构成了压裂缝渗流模型，利用Laplace变换、二重积分可得

$ \overline {{\psi _{{\rm{wD}}}}} - \overline {{\psi _{{\rm{fD}}j}}} ({x_{\rm{D}}}, {y_{{\rm{wD}}j}}, s) = \dfrac{{\mathsf{π}}}{{{C_{{\rm{FD}}j}}}}\left( {{x_{\rm{D}}}\overline {{q_{{\rm{D}}j}}} - \int\limits_0^{{x_{\rm{D}}}} {\int\limits_0^\sigma {\overline {{q_{{\rm{fD}}j}}} ({x_{\rm{D}}}, s){\rm{d}}{x_{\rm{D}}}{\rm{d}}\sigma } } } \right) $

(10)

式中：

${{\psi _{{\rm{wD}}}}}$—无因次井底压力，根据其定义，有

${\psi _{{\rm{wD}}}} = \dfrac{1}{s}$；

${{q_{{\rm{D}}j}}}$—第$j$条压裂缝的无因次产量；

$\overline {{q_{{\rm{D}}j}}} = \dfrac{{{p_{{\rm{sc}}}}T}}{{{\mathsf{π}}Kh{T_{{\rm{sc}}}}\left( {{\psi _{\rm{i}}} - {\psi _{\rm{w}}}} \right)}}{q_j}$；

${q_j}$—第$j$条压裂缝的产量，m$^3$/s。

2.3 气藏渗流模型与多压裂缝渗流模型联立

联立式(6)和式(10)，得

$ \overline {{\psi _{{\rm{wD}}}}} - \dfrac{1}{2}\sum\limits_{l = 1}^M {\int\limits_{ - {x_{{\rm{fD}}l}}}^{{x_{{\rm{fD}}l}}} {\overline {{q_{{\rm{fD}}l}}} ({x_{{\rm{wD}}}}, s){{\rm K}_0}\left( {\sqrt s \sqrt {{{({x_{\rm{D}}} - {x_{{\rm{wD}}}})}^2} + {{({y_{{\rm{wD}}j}} - {y_{{\rm{wD}}l}})}^2}} } \right){\rm{d}}{x_{{\rm{wD}}}}} }= \\ {\kern 6em} \dfrac{{\mathsf{π}}}{{{C_{{\rm{FD}}j}}}}\left( {{x_{\rm{D}}}\overline {{q_{{\rm{D}}j}}} - \int\limits_0^{{x_{\rm{D}}}} {\int\limits_0^\sigma {\overline {{q_{{\rm{fD}}j}}} ({x_{\rm{D}}}, s){\rm{d}}{x_{\rm{D}}}{\rm{d}}\sigma } } } \right), \hspace{3em} j=1, 2, 3, M $

(11)

式(11)为Fredholm积分方程，可通过对内边界进行单元离散来进行求解。

如图 2所示，假设第$j$条压裂缝上第$i$个单元上的中点坐标为$({\hat x_{{\rm{D}}j, i}}, {y_{{\rm{wD}}j}})$，第$j$条压裂缝上第$i$个单元的两个端点坐标分别为$({x_{{\rm{D}}j, i}}, {y_{{\rm{wD}}j}})$和$({x_{{\rm{D}}j, i + 1}}, {y_{{\rm{wD}}j}})$。

离散后的方程为

$ \dfrac{{{\psi _{{\rm{wD}}}}}}{s} - \dfrac{1}{2}\sum\limits_{l = 1}^M {\sum\limits_{v = 1}^N {\int\limits_{{x_{{\rm{D}}l, v}}}^{{x_{{\rm{D}}l, v + 1}}} {\overline {{q_{{\rm{fD}}l, v}}} (s){\rm{G}}({{\hat x}_{{\rm{D}}j, i}}, {x_{{\rm{wD}}}}, {y_{{\rm{wD}}j}}, {y_{{\rm{wD}}l}}, s){\rm{d}}{x_{{\rm{wD}}}}} } } =\\ {\kern 6em} \dfrac{{\mathsf{π}}}{{{C_{{\rm{FD}}j}}}}\left( {{{\hat x}_{{\rm{D}}j, i}}\overline {{q_{{\rm{D}}j}}} - \sum\limits_{v = 1}^{i - 1} {\overline {{q_{{\rm{fD}}j, v}}} } \left( {i - v} \right)\Delta {x_{{\rm{D}}j, i}}^2 - \overline {{q_{{\rm{fD}}j, i}}} \dfrac{{\Delta {x_{{\rm{D}}j, i}}}}{8}} \right) , \\ j=1, 2, ... M; i=1, 2, ..., N $

(12)

其中：

$ {\rm{G}}({x_{\rm{D}}}, {x_{{\rm{wD}}}}, {y_{{\rm{wD}}j}}, {y_{{\rm{wD}}l}}, s) =\\ {\kern 3em} {{\rm{K}}_0}\left( {\sqrt s \sqrt {{{({x_{\rm{D}}} - {x_{{\rm{wD}}}})}^2} + {{({y_{{\rm{wD}}j}} - {y_{{\rm{wD}}l}})}^2}} } \right) +\\ {\kern 4em} {{\rm{K}}_0}\left( {\sqrt s \sqrt {{{({x_{\rm{D}}} +{x_{{\rm{wD}}}})}^2} + {{({y_{{\rm{wD}}j}} - {y_{{\rm{wD}}l}})}^2}} } \right) $

(13)

此外，对于第$j$条压裂缝，还存在如下流量关系

$ \sum\limits_{v = 1}^N {\overline {{q_{{\rm{fD}}j, v}}} \Delta {x_{{\rm{D}}j, v}}} = \overline {{q_{{\rm{D}}j}}} $

(14)

当$i$、$j$变化时，由式(12)和式(14)可得到$N\times M+M$个线性代数方程，未知数(${{q_{{\rm{fD}}l, v}}}$、${{q_{{\rm{D}}j}}}$)的个数也为$N\times M+M$个，方程个数与未知数个数相等，可以封闭求解，此外，在将Laplace空间解反演到真实时间域时可采用Stehfest数值反演^[26]。

求出各压裂缝无因次产量$q_{{\rm D}j}$之后，利用其无因次定义，可求得各压裂缝的有因次产量$q_j$，而气井产量$q$则为各压裂缝产量$q_j$之和

$ q = \sum\limits_{j = 1}^M {{q_j}} $

(15)

3 产能影响因素定量计算及分析

以某低渗透气田T2区块典型储层物性参数为例，定量分析多级压裂水平井产能影响因素。参数情况如下：压裂缝条数$M$=6，压裂缝半长$x_{\rm f}$=40 m，水平井长度$L_{\rm h}$=600 m，$r_{\rm w}$ =0.1 m，地面大气压力$p_{\rm{sc}}$=0.101 MPa，地面温度$T_{\rm{sc}}$=20 ℃，气藏温度$T$=102 ℃，气藏厚度$h$=12 m，原始地层压力$p_{\rm i}$=30.5 MPa，井底压力$p_{\rm w}$=22 MPa，气藏渗透率$K$=0.001 D，裂缝渗透率$K_{\rm f}$ =20 D，气藏孔隙度$\phi$=0.13，天然气黏度$\mu$=0.025 mPa$\cdot$s，天然气的压缩系数$C_{\rm t}$=0.033 MPa$^{-1}$，井筒半径$r_{\rm w}$=0.1 m，裂缝宽度$W_{\rm f}$=0.002 m。

经单位转换后，可利用本文模型进行产能计算分析。讨论某一参数的影响时，只令该参数变化，其他参数除非在图上标注，否则，取上述给定值。

3.1 气层有效厚度$h$

图 3是气层有效厚度对MFHW产量的影响图。从图可以看出，随气层有效厚度的增大，气井产量几乎呈线性增加，故气层有效厚度是影响气井产能的主要因素之一。

3.2 气藏渗透率$K$

图 4是气藏渗透率$K$对气井产量的影响图。可以看出，随渗透率$K$的增大，气井产量也几乎呈线性增加。故气藏渗透率也是影响气井产能的主要因素之一。

3.3 压裂缝条数$M$

图 5是压裂缝条数$M$对气井产量的影响图。从图上可以看出，当裂缝条数$M$增加时，气井产量$q_{\rm{sc}}$相应增大，但它们并不成线性增长关系，在裂缝条数较少时，产量增大的幅度很明显(如图 5所示，压裂缝条数从2条增加至4条时，产能增幅很大)，但随着裂缝条数的增加，产量增大的幅度逐渐减小。因此，对于水平井分段多级压裂来说，虽然可通过增加压裂缝条数来增大气井产量，但也不能一味地增加压裂缝条数。

3.4 压裂缝半长$x_{\rm f}$

图 6是压裂缝半长$x_{\rm f}$对气井产量的影响图。从图上可以看出，压裂缝半长$x_{\rm f}$在较短时，随它的增加，产量增大的幅度较大；压裂缝半长$x_{\rm f}$较长时，随它的增加，产量增大的幅度下降。早期$x_{\rm f}$对流量的影响大于晚期。

3.5 压裂缝导流能力$K_{\rm f}W_{\rm f}$

图 7是压裂缝导流能力$K_{\rm f}W_{\rm f}$对气井产量的影响图。可以看出，当$K_{\rm f}W_{\rm f}$较小时，随$K_{\rm f}W_{\rm f}$的增加，产量增加较大；但当$K_{\rm f}W_{\rm f}$增加到一定值(例如，在T2区块典型参数条件下，当$K_{\rm f}W_{\rm f}$增大到40 D$\cdot$mm)后，产量增大的幅度很小。说明可以通过注入高品质、足够数量支撑剂等方法提高压裂缝导流能力(实际上是提高了压裂缝的$K_{\rm f}$、$W_{\rm f}$)，从而提高气井产能；另外，要综合考虑气井产能的增加幅度与提高$K_{\rm f}W_{\rm f}$时所花费的成本，而不要一味追求过于高的导流能力。

从以上对压裂缝条数$M$、压裂缝半长$x_{\rm f}$、压裂缝导流能力$K_{\rm f}W_{\rm f}$的分析可看出，三者对产能影响的规律比较相似，均是随其增大，产量增大，但前期增速快，后期增速慢。

3.6 各压裂缝流量分布规律

图 8是压裂缝条数$M$=6时，在不同生产时间从地层流入各条压裂缝的流量大小分布图。图中对压裂缝的编号见图 1。

从图 8可以看出，在早期，地层流入各条压裂缝的流量近似相等，但随着生产时间的推移，从地层流入两端压裂缝的流量逐渐大于从地层流入中部压裂缝的流量，即呈现“两端高、中部低”的特点。

4 结论

(1) 为多级压裂水平气井产能计算提供了一个严格的数学模型，并对模型成功地进行了求解，然后以此为基础，对不同因素影响下的多级压裂水平气井产能进行了定量计算和分析。

(2) 随气层有效厚度$h$的增大，气井产量$q$几乎呈线性增大；随气藏渗透率$K$的增大，气井产量$q$也几乎呈线性增大。

(3) 压裂缝条数$M$、压裂缝半长$x_{\rm f}$、压裂缝导流能力$K_{\rm f}W_{\rm f}$对产能影响的规律相似：当$M$、$x_{\rm f}$、$K_{\rm f}W_{\rm f}$值较小时，随它们的增大，产量增大幅度较大，当它们值较大时，产量增大幅度变小。故不能一味地靠增加$M$、$x_{\rm f}$、$K_{\rm f}W_{\rm f}$的大小来增加气井产能，应综合考虑气井产能的增加幅度与增大$M$、$x_{\rm f}$、$K_{\rm f}W_{\rm f}$时所花费的成本。

(4) 在早期，地层流入各条压裂缝的流量近似相等，随生产时间的推移，地层流入端部压裂缝的流量逐渐大于流入中部压裂缝的流量，即呈现“两端高、中部低”的特点。