[发明专利]一种快速预测纵向非均质气藏压裂复杂缝网水平井产能的方法

权利要求书

1.一种快速预测纵向非均质气藏压裂复杂缝网水平井产能的方法，其特征在于，包括以下步骤：

S1：运用测井、地震、室内实验手段获取产能计算参数；

S2：根据步骤S1中获得的测井资料，将储层从纵向上划分为若干个小层；

S3：针对步骤S2中获得的任意一小层，将穿过该小层的裂缝进行离散化处理，利用该小层的产能计算参数，基于实空间源汇函数理论，计算出该层对每一个裂缝微元的压降，再运用压降叠加原理得到该层对裂缝的总压降；

S4：根据达西定律，得到每一个裂缝微元到水平井筒的压降；

S5：将步骤S3、S4得到的计算公式相加，得到该小层的耦合渗流方程组，运用数值迭代法求解得该小层的产量；

S6：对步骤S2中划分的储层内其它小层重复S3～S5步骤，然后将所有小层的产量相加得到水平井总产量。

2.如权利要求1中所述的一种快速预测纵向非均质气藏压裂复杂缝网水平井产能的方法，其特征在于，所述步骤S1～S6存在以下假设：

1)压裂水平井位于储层几何中心；

2)水力裂缝垂直于水平井并且完全贯穿储层；

3)储层内流动为单相气体流动且忽略重力的影响；

4)每小层气体流动互不干扰，并且仅通过水力裂缝流入井筒；

5)不考虑井筒压降。

3.如权利要求1中所述的一种快速预测纵向非均质气藏压裂复杂缝网水平井产能的方法，其特征在于，所述步骤S1中所述产能计算参数包括：

1)气藏参数

气藏长度、气藏宽度、气藏厚度、基质孔隙度、基质渗透率、气藏温度、气藏初始压力、综合压缩系数；

2)裂缝参数

裂缝长度、裂缝宽度、裂缝间距、裂缝渗透率；

3)流体参数

气体黏度、偏差因子、井底流压、标准温度、标准压力。

4.如权利要求1中所述的一种快速预测纵向非均质气藏压裂复杂缝网水平井产能的方法，其特征在于，所述步骤S2中将储层从纵向上划分为若干个小层，主要依据的是各个小层物性参数在测井曲线上反映的差异。

5.如权利要求1中所述的一种快速预测纵向非均质气藏压裂复杂缝网水平井产能的方法，其特征在于，所述步骤S3中任意裂缝微元对该层中任意一点产生的压降可以写成：

式中，p_int为气藏初始压力，MPa；p(x,y,z,t)为该层任意一点的压力，MPa；为第j条裂缝上翼第i个微元的流量，m³/s；μ_g为气体黏度，mPa·s；p_sc为标准压力，MPa；Z为偏差因子；T为气藏温度，K；T_sc为标准温度，K；C_t为综合压缩系数，1/MPa；φ_layer为该层基质孔隙度；上标(a)代表上翼，(b)代表下翼，上下翼方程具有镜像性，只给出上翼方程，下翼方程推导方法一致；

其中，三个方向的基本汇函数为：

式中，x_e、y_e、z_e为该层长度、宽度、厚度，m；x_w、y_w、z_w为点汇坐标，m；x、y、z为该层中任意一点坐标，m；t为时间，s；K_layer为该层基质渗透率，mD。

6.如权利要求1中所述的一种快速预测纵向非均质气藏压裂复杂缝网水平井产能的方法，其特征在于，所述步骤S3中的运用压降叠加原理得到该层对上翼裂缝的总压降为：

式中，为第j条裂缝上翼第i个微元压力，MPa；R(i,j)为气藏压降系数，MPa/m³；该层对下翼裂缝的总压降推导方法一致。