一种窄河道储层的产量递减分析方法

摘要

本发明提出了一种窄河道储层的产量递减分析方法。该方法包括：根据储层情况和裂缝参数建立物理模型；裂缝模型和地层模型耦合；裂缝模型和地层模型耦合求解，通过将裂缝离散成N个等长均匀流量段，利用叠加原理，通过裂缝壁面流量和压力连续性条件将储层解析模型与裂缝离散模型进行耦合求解；采用直接调用圆形封闭储层F函数计算裂缝导流能力影响函数，建立归一化裂缝导流能力影响函数，建立获取有限导流裂缝拟稳态常数bDpss计算新模型；基于所述有限导流裂缝拟稳态常数bDpss计算新模型，绘制圆形封闭储层中垂直裂缝井产量递减新图版。将产量递减新图版应用在产量递减分析中。该方法提高了压后气井产量递减分析的精度，适合矿场应用。

说明书

技术领域

本发明属于石油天然气开采领域，尤其涉及油气藏产能评价领域，具体涉及一种窄河道储层的产量递减分析方法。

背景技术

油气藏产量递减分析是一种高效分析油气井实际生产数据的方法，利用无量纲产量递减曲线图版对油气井有效的生产数据进行拟合分析，可以快速计算生产井和储层的某些参数，例如原始地质储量、控制面积、有效渗透率等，进而可以评价油气井和储层。关于油气藏产量递减分析方法，国内外公开发表的期刊文献有很多。文献中目前报道的各种油气藏产量递减分析方法中，主要包括：Arps产量递减方程Fetkovich产量递减方法。Arps产量递减方程是一种全面分析矿场实际生产数据的产量递减方法，Arps递减方程可以简单地分为三种表达形式：双曲形式、指数形式、调和形式。Fetkovich产量递减方法从圆形有界均质地层中一口普通直井入手，将传统试井分析中的不稳定渗流方程引入递减分析之中，使采用Arps产量递减方程建立的递减曲线图版的应用范围更大。但是，Arps产量递减递减方程和Fetkovich产量递减方法只能在定井底流压条件下使用，在变井底流压和变产量生产时并不适用；对于井控半径大，裂缝穿透比(裂缝半长和井控半径的比值)较大的油气藏也不适用。窄河道油气藏垂直井进行水力压裂改造后形成了一条垂直裂缝，大大增加了井与地层的接触面积，同时改善了井底周围储层的渗流条件，油气井在窄河道油气藏中井控半径大，导致裂缝穿透比较大，井底流压和油气井产量也发生改变。目前，还未见专门用于窄河道油气藏的产量递减分析方法。

发明内容

为解决上述技术问题，本发明提供的窄河道储层的产量递减分析方法。包括以下步骤：步骤1，针对窄河道储层具有一条垂直裂缝的直井，根据储层情况和裂缝参数建立物理模型；步骤2，裂缝模型和地层模型耦合；步骤3，裂缝模型和地层模型耦合求解，通过将裂缝离散成N个等长均匀流量段，利用叠加原理，求取每个离散段的无量纲压力，通过裂缝壁面流量和压力连续性条件将油藏解析模型与裂缝离散模型进行耦合求解；步骤4，采用直接调用圆形封闭储层F函数计算裂缝导流能力影响函数，建立归一化裂缝导流能力影响函数，并以此建立获取有限导流裂缝拟稳态常数b_Dpss计算新模型；步骤5，基于所述有限导流裂缝拟稳态常数b_Dpss计算新模型，绘制圆形封闭储层中垂直裂缝井产量递减新图版，结合实际储层数据，将所述圆形封闭储层中垂直裂缝井产量递减新图版应用在产量递减数据拟合分析中。

在所述步骤S1中，所述物理模型包括裂缝模型和地层模型；所述储层情况和所述裂缝参数设置方法如下：

选取地层均质、等厚，上下为不渗透边界、水平方向具有封闭边界的圆形储层；

设置油井产量完全由裂缝产出，一条有限导流裂缝完全贯穿地层，且忽略裂缝两端的流体流入的情况；

设置地层中流体为微可压缩流体，流体在地层和裂缝中流动时符合达西定律；井筒中的流体仅由裂缝中流入，不考虑从地层流入的情况；流体由地层流入裂缝时流量均匀分布。

所述裂缝模型的建立方法如下：

(1)定义如下无量纲参数，所述无量纲参数包括：

无量纲坐标量，无量纲时间，无量纲导流能力，无量纲地层压力，无量纲裂缝压力，无量纲产量，无量纲裂缝截面流量，无量纲线流量

其中x为坐标量、x_f为裂缝半长，m；k为地层渗透率，mD；t为时间，d；φ为孔隙度；μ为气体粘度，mPa·s；C_t为地层压缩系数；k_f为裂缝渗透率，mD；w_f为裂缝宽度，m；h为地层厚度，m；p_i为原始地层压力、p为地层压力、p_f为裂缝压力，MPa；B为流体体积系数；表示流体由地层流入单位长度裂缝的流量、q_c表示t时刻时，任意点x处的裂缝流量，m³/d；

(2)建立裂缝模型的解析式，即建立两翼对称裂缝的拉氏空间压力解析式，所述两翼对称裂缝的拉氏空间压力解析式为：

其中表示拉氏空间井底流压，无因次；表示拉氏空间裂缝压力，无因次；C_fD为无量纲裂缝导流能力，无因次；为拉氏空间无量线流量，无因次；s为拉普拉斯变量；x_D为无量纲坐标量，无因次；v、u为积分系数，无因次。

所述地层模型的建立方法如下：

(1)定义如下无量纲参数，所述无量纲参数包括：

无量纲时间，无量纲柱坐标，无量纲圆形地层半径，无量纲地层压力，

其中k为地层渗透率，mD；t为时间，d；为孔隙度；μ为气体粘度，mPa·s；C_t为地层压缩系数；x_f为裂缝半长，m；r_D为无量纲圆形地层半径，无因次；r为柱坐标，m；r_e为井控半径，m；h为地层厚度，m；p_i为原始地层压力、p为地层压力，MPa；p_D为无量纲地层压力，无因次；q为产量，m³/d；B为流体体积系数；

(2)建立地层模型的解析式，即建立圆形封闭垂直裂缝井地层压力分布在拉普拉斯空间的解析式，所述圆形封闭垂直裂缝井地层压力分布在拉普拉斯空间的解析式为：

其中：为拉普拉斯空间中裂缝方向上的流量分布；s为拉普拉斯变量；I₀(.)为第一类零阶修正贝塞尔函数；I₁(.)为第一类一阶修正贝塞尔函数；K₀(.)为第二类零阶修正贝塞尔函数；K₁(.)——第二类一阶修正贝塞尔函数；r_eD为无量纲井控半径，无因次；x_D为无量纲坐标量，无因次；α为积分系数。

在所述步骤2中，所述裂缝模型和地层模型耦合方法为：

(1)不考虑表皮引起的压降影响，在裂缝和地层相交面上，建立地层压力和裂缝压力相等的关系式：

(2)将所述地层模型的解析式和所述地层压力和裂缝压力相等的关系式代入所述裂缝模型的解析式中，耦合后得到圆形封闭地层中有限导流垂直裂缝井的压力解析式：

其中为拉普拉斯空间中裂缝方向上的流量分布；s为拉普拉斯变量；I₀(.)为第一类零阶修正贝塞尔函数；I₁(.)为第一类一阶修正贝塞尔函数；表示拉氏空间中的井底流压，无因次；为拉氏空间中无量纲线流量，无因次；r_eD为无量纲井控半径，无因次；x_D为无量纲坐标量，无因次；C_fD为无量纲裂缝导流能力，无因次；s为拉普拉斯变量；α、v、u为积分系数，无因次。

在所述步骤3中，所述裂缝模型和地层模型耦合求解是通过将裂缝离散成N个等长均匀流量段，利用叠加原理，求取每个离散段的无量纲压力。通过裂缝壁面流量和压力连续性条件将油藏解析模型与裂缝离散模型进行耦合求解。

对于未知数和采用离散求解方法可以得到一个N+1阶线性方程组：

拉普拉斯空间中，定产压力解和定压产量解的关系式为：

其中A_ij、B_ij表示方程系数；为拉普拉斯空间中裂缝方向上的流量分布，无因次；q_D为无量纲产量，无因次；p_wD为无量纲井底流压，无因次。

利用上式，结合数值反演，得到q_D和t_D的关系，可绘制q_D和t_D的关系曲线图版，q_D和t_D的关系曲线图版前端呈现归一化，后端为发散状态，该曲线图版不利于实际数据拟合。

为建立拟合精度更高的产量递减图版，本发明引入拟稳态常数b_Dpss，下面介绍拟稳态常数b_Dpss的求解方法。

在所述步骤4中，所述圆形封闭储层F函数用于求取无量纲生产指数J_D，所述圆形封闭储层F函数为：

上式中σ(x_D，y_D)和δ(x_D，y_D)分别为：

其中r_eD为无量纲井控半径，无因次；x_D、y_D为无量纲坐标量，无因次。

无量纲生产指数J_D与所述有限导流裂缝拟稳态常数b_Dpss互为倒数关系：

p_wD-p_avgD＝b_Dpss＝1/J_D

其中，P_wD为无量纲井底压力，无因次；p_avgD为无量纲平均地层压力，无因次；b_Dpss为有限导流裂缝拟稳态常数，无因次；J_D为无量纲生产指数，无因次。

在所述步骤4中，所述归一化裂缝导流能力影响函数通过求取的拟稳态常数b_Dpss求得，所述归一化裂缝导流能力影响函数为：

其中，I_x为穿透比，f₀(C_fD)和f₁(C_fD)分别表示穿透比I_x为0和1时的影响函数。

在所述步骤4中，所述有限导流裂缝拟稳态常数b_Dpss计算新模型为：

其中b_Dpss,FC(I_x,C_fD)表示有限导流能力拟稳态常数值；b_Dpss,IC(I_x)表示无限导流能力影响的拟稳态常数值。

在所述步骤5中，还包括重新定义无量纲产量q_Dd和无量纲时间t_Dd，并结合所述有限导流裂缝拟稳态常数b_Dpss计算新模型，定义无量纲产量积分q_Ddi和产量积分导数q_Ddid，所述重新定义无量纲产量为：

q_Dd＝b_Dpss·q_D

所述重新定义无量纲时间为：

所述无量纲产量积分表达式为：

所述无量纲产量积分导数表达式为：

其中，q_Dd为重新定义的无量纲产量，无因次；t_Dd为重新定义的无量纲时间，无因次；q_Ddi为无量纲产量积分，无因次；q_Ddid为产量积分导数，无因次；N_pDd为无量纲累计产量积分，无因次；r_eD为无量纲井控半径，无因次；t_D为无量纲，无因次。

通过上述定义式，并结合无量纲产量q_Dd和t_Dd的关系，进而可以绘制无量纲q_Ddi，q_Ddid和t_Dd曲线的关系曲线，即可绘制圆形封闭储层中垂直裂缝井产量递减新图版

所述绘制圆形封闭储层中垂直裂缝井产量递减新图版应用方法包括：

时间数据处理：对于气井，计算物质平衡拟时间：

产量数据处理：对于气井，计算拟压力规整化产量：

计算拟压力规整化产量积分：

计算拟压力规整化产量积分导数：

在log～log双对数坐标图上分别绘制处理后的产量和处理后的时间的关系曲线，共计三组曲线：(q_g/Δp_p)～t_ca、(q_g/Δp_p)_i～t_ca、(q_g/Δp_p)_id～t_ca；

同时将得到的三组曲线或者将其中两组任意组合，与所述的产量递减新图版进行拟合分析，使得每组曲线能够最大限度地获得较好的拟合效果；

根据拟合结果记录无量纲圆形地层半径r_eD；

在上述拟合曲线上选取任意一个拟合点，记录实际值和理论拟合值，即分别为(t_ca，q_g/Δp_p)、(t_Dd，q_Dd)；

根据记录得到的C_fD和r_eD以及运用有限导流裂缝拟稳态常数b_Dpss计算新模型计算b_Dpss；

综合以上结果，计算油气田参数，原始地质储量：

井控面积、井控半径：

地层有效渗透率：

裂缝半长：

其中t_ca为物质平衡时间，d；μ_g为气体粘度，mPa·s；φ为孔隙度，％；C_g为气体压缩系数，MPa^-1；q_g表示产量，m³/d；G为地质储量，10⁴m³；N_p为累积产油量，10⁴m³；p_p为规整化拟压力，MPa；p_i为原始地层压力，MPa；p_wf为井底压力，MPa；q/△p规整化产量，m³/d/MPa；(q/△p)_i为规整化产量积分，m³/d/Mpa；(q/△p)_id为规整化产量积分导数，m³/d/Mpa；x_f为裂缝半长，m；r_e为圆形地层半径，m；r_eD为无量纲井控半径；S_wi为束缚水饱和度，％；k_g，k_o为油、气有效渗透率，mD；A为井控面积，m²；B_gi，B_oi为原始气、油体积系数；b_Dpss为稳态常数。

本发明提供的窄河道储层的产量递减分析方法克服了常规产量递减分析方法在变井底流压和变产量生产时并不适用；对于井控半径大，裂缝穿透比(裂缝半长和井控半径的比值)较大的油气藏也不适用的缺点；解决了压后水平井瞬时产量高、稳产能力差，传统方法不可靠，合理产量确定无章可循的问题；本发明提供的产量递减新图版提高了压后气井产量递减分析的精度。

本发明的其它特征和优点将在随后的说明书中阐述，并且部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。

附图说明

附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例共同用于解释本发明，并不构成对本发明的限制。在附图中：

图1是本发明窄河道储层的产量递减分析方法的流程示意图；

图2是圆形封闭地层中一口垂直裂缝井示意图；图中：x_f为裂缝半长，m；w_f为裂缝宽度，m；

图3是裂缝中一维流动示意图；图中：q_c(x,t)表示t时刻时，任意点x处的裂缝流量，两端积分表示两端无流量流入裂缝；

图4是裂缝离散示意图；图中：x_Di表示第i的离散单元中心点坐标，q_fDi表示第i离散单元无量纲流量，△x表示每个离散单元长度；

图5是C_fD＝0.1和0.5时的无量纲产量q_D和时间曲线t_D；图中，q_D为无量纲产量，无因次；t_D为无量纲时间，无因次；

图6是某井生产曲线图；图中，q为产量，10⁴m³；p_wf为井底压力，MPa；

图7是某井圆形封闭储层中垂直裂缝井产量递减新图版；图中，q_D为无量纲产量，q_Ddi为无量纲产量积分，q_Ddid为产量积分导数，C_fD为无量纲裂缝导流能力，x_f为裂缝半长，m；r_e为井控半径，m。

具体实施方式

下面结合实施例，对本发明作进一步详细说明，但本发明的实施方式不仅限于此。

图1示出了根据本发明的窄河道储层的产量递减分析方法的流程。如图1所示，本发明提供的一种窄河道储层的产量递减分析方法，该方法包括：

步骤1：针对窄河道储层具有一条垂直裂缝的直井，根据储层情况和裂缝参数建立物理模型；

步骤2：裂缝模型和地层模型耦合；

步骤3：裂缝模型和地层模型耦合求解，通过将裂缝离散成N个等长均匀流量段，利用叠加原理，求取每个离散段的无量纲压力，通过裂缝壁面流量和压力连续性条件将油藏解析模型与裂缝离散模型进行耦合求解；

步骤4：采用直接调用圆形封闭储层F函数计算裂缝导流能力影响函数，建立归一化裂缝导流能力影响函数，并以此建立获取有限导流裂缝拟稳态常数b_Dpss计算新模型；

步骤5：基于所述有限导流裂缝拟稳态常数b_Dpss计算新模型，绘制圆形封闭储层中垂直裂缝井产量递减新图版，结合实际储层数据，将所述圆形封闭储层中垂直裂缝井产量递减新图版应用在产量递减数据拟合分析中。

在所述步骤S1中，所述物理模型包括裂缝模型和地层模型；所述储层情况和所述裂缝参数设置方法如下：

(1)选取地层均质、等厚，上下为不渗透边界、水平方向具有封闭边界的圆形储层；

(2)设置油井产量完全由裂缝产出，一条有限导流裂缝完全贯穿地层，且忽略裂缝两端的流体流入的情况；

(3)设置地层中流体为微可压缩流体，流体在地层和裂缝中流动时符合达西定律；井筒中的流体仅由裂缝中流入，不考虑从地层流入的情况；流体由地层流入裂缝时流量均匀分布。

图2示出了圆形封闭地层中一口垂直裂缝井模型示意图，所述的裂缝模型和地层模型的解析式都是基于此模型建立的。

所述裂缝模型的建立方法如下：

(1)定义如下无量纲参数，所述无量纲参数包括：无量纲坐标量，无量纲时间，无量纲导流能力，无量纲地层压力，无量纲裂缝压力，无量纲产量，无量纲裂缝截面流量，无量纲线流量

其中x为坐标量、x_f为裂缝半长，m；k为地层渗透率，mD；t为时间，s；φ为孔隙度；μ为气体粘度，mPa·s；C_t为地层压缩系数；k_f为裂缝渗透率，mD；w_f为裂缝宽度，m；h为地层厚度，m；p_i为原始地层压力、p为地层压力、p_f为裂缝压力，MPa；B为流体体积系数；表示流体由地层流入单位长度裂缝的流量、q_c表示t时刻时，任意点x处的裂缝流量，m³/d。

(2)建立裂缝模型的解析式，即建立两翼对称裂缝的拉氏空间压力解析式，所述两翼对称裂缝的拉氏空间压力解析式建立方法为：

建立控制方程：

设置初始地层条件：p_fD(x_D，t_D＝0)＝0

设置内边界条件：

设置外边界条件：

同时，裂缝截面流量和裂缝线流量存在以下关系式：

将公式所述裂缝截面流量和裂缝线流量的关系式代入所述的控制方程，并对x_D二次求导后，可以得到中间方程Ⅰ：

将所述的内边界条件和外边界条件及所述初始地层条件公式代入所述中间方程Ⅰ，求解所述中间方程Ⅰ后可以得到中间方程Ⅱ：

当x_D＝0时为井底位置，此时的裂缝压力为井底压力，则有：

p_fD(0)＝p_wD

将公式p_fD(0)＝p_wD代入所述中间方程Ⅱ得到中间方程Ⅲ：

对所述中间方程Ⅲ进行拉普拉斯变换后，可以得到裂缝一翼的拉氏空间压力解析式：

进一步考虑裂缝的对称两翼的情况，则可得到两翼对称裂缝的拉氏空间压力解析式：

其中q_cD为无量纲裂缝截面流量，无因次；表示拉氏空间中的井底流压，无因次；表示拉氏空间中的裂缝压力，无因次；C_fD为无量纲裂缝导流能力，无因次；为拉氏空间中无量线流量，无因次；s为拉普拉斯变量；x_D为无量纲坐标量，无因次；v、u为积分系数。

所述地层模型的建立方法如下：

(1)定义如下无量纲参数，所述无量纲参数包括：无量纲时间，无量纲柱坐标，无量纲圆形地层半径，无量纲地层压力，

其中k为地层渗透率，mD；t为时间，d；为孔隙度；μ为气体粘度，mPa·s；C_t为地层压缩系数；x_f为裂缝半长，m；r_D为无量纲圆形地层半径，无因次；r为柱坐标，m；r_e为井控半径，m；h为地层厚度，m；p_i为原始地层压力、p为地层压力，MPa；p_D为无量纲地层压力，无因次；q为产量，m³/d；B为流体体积系数，无因次；

(2)建立地层模型的解析式，即建立圆形封闭垂直裂缝井地层压力分布在拉普拉斯空间的解析式，所述圆形封闭垂直裂缝井地层压力分布在拉普拉斯空间的解析式建立方法为：

建立地层渗流控制方程：

设置初始地层条件：p_D(r_D，t_D＝0)＝0

设置内边界条件：

设置外边界条件：

将所述的内边界条件和外边界条件及所述初始地层条件公式代入所述层渗流控制方程，并利用Ozkan的点源函数方法可以得到点源解，进一步得到平面源的解析式。在拉普拉斯空间中，圆形封闭储层的点汇压力解析式为：

对所述圆形封闭储层的点汇压力解析式进行叠加积分，得到圆形封闭垂直裂缝井地层压力分布在拉普拉斯空间的解析式为：

其中，为拉普拉斯空间中裂缝方向上的流量分布；s为拉普拉斯变量；I₀(.)为第一类零阶修正贝塞尔函数；I₁(.)为第一类一阶修正贝塞尔函数；r_D为无量纲柱坐标，无因次；r_eD为无量纲井控半径，无因次；t_D为无量纲时间，无因次；x_D、y_D为无量纲坐标量，无因次；p_D为无量纲地层压力，无因次；α为积分系数。

在所述步骤2中，所述裂缝模型和地层模型耦合方法为：

(1)不考虑表皮引起的压降影响，在裂缝和地层相交面上，建立地层压力和裂缝压力相等的关系式：

(2)将所述地层模型的解析式和所述地层压力和裂缝压力相等的关系式代入所述裂缝模型的解析式中，耦合后得到圆形封闭地层中有限导流垂直裂缝井的压力解析式：

其中为拉普拉斯空间中裂缝方向上的流量分布；s为拉普拉斯变量；I₀(.)为第一类零阶修正贝塞尔函数；I₁(.)为第一类一阶修正贝塞尔函数；表示拉氏空间中的井底流压，无因次；为拉氏空间中无量纲线流量，无因次；r_eD为无量纲井控半径，无因次；p_D为无量纲地层压力，无因次；x_D为无量纲坐标量，无因次；C_fD为无量纲裂缝导流能力，无因次；s为拉普拉斯变量；α、v、u为积分系数，无因次。

虽然耦合后得到圆形封闭地层中有限导流垂直裂缝井的压力解析式只有和两个未知数，但直接求解非常困难。因此，需要进行数值离散求解。

在所述步骤3中，所述裂缝模型和地层模型耦合求解是通过将裂缝离散成N个等长均匀流量段，利用叠加原理，求取每个离散段的无量纲压力。通过裂缝壁面流量和压力连续性条件将油藏解析模型与裂缝离散模型进行耦合求解。

图3和图4示出了裂缝中一维流动示意图和裂缝离散示意图，所述裂缝模型和地层模型的耦合求解基于此裂缝中一维流动模型和裂缝离散方法。

对于未知数和采用离散求解方法可以得到一个N+1阶线性方程组：

拉普拉斯空间中，定产压力解和定压产量解的关系式为：

利用上式，结合数值反演，得到q_D和t_D的关系，可绘制q_D和t_D的关系曲线图版，q_D和t_D的关系曲线图版(图5所示)曲线前端呈现归一化，后端为发散状态，不利于实际数据拟合。

为建立拟合精度更高的产量递减图版，本发明引入拟稳态常数b_Dpss，下面介绍拟稳态常数b_Dpss的求解方法。

在所述步骤4中，所述圆形封闭储层F函数用于求取无量纲生产指数J_D，所述圆形封闭储层F函数为：

上式中σ(x_D，y_D)和δ(x_D，y_D)分别为：

其中r_eD为无量纲井控半径，无因次；x_D、y_D为无量纲坐标量，无因次。

无量纲生产指数J_D与所述有限导流裂缝拟稳态常数b_Dpss互为倒数关系：

p_wD-p_avgD＝b_Dpss＝1/J_D

其中，P_wD为无量纲井底压力，无因次；p_avgD为无量纲平均地层压力，无因次；b_Dpss为有限导流裂缝拟稳态常数，无因次；J_D为无量纲生产指数，无因次。

在所述步骤4中，所述归一化裂缝导流能力影响函数通过求取的拟稳态常数b_Dpss求得，所述归一化裂缝导流能力影响函数为：

其中，I_x为穿透比，f₀(C_fD)和f₁(C_fD)分别表示穿透比I_x为0和1时的影响函数。

在所述步骤4中，所述有限导流裂缝拟稳态常数b_Dpss计算新模型为：

其中b_Dpss,FC(I_x,C_fD)表示有限导流能力拟稳态常数值；b_Dpss,IC(I_x)表示无限导流能力影响的拟稳态常数值。

由于实际数据中存在异常值，为了降低产量数据噪声波动的影响，使得生产数据在一定程度上便于拟合分析，且q_D和t_D的关系曲线图版曲线前端呈现归一化，后端为发散状态，也不利于实际数据拟合。

在所述步骤5中，还包括重新定义无量纲产量q_Dd和无量纲时间t_Dd，并结合所述有限导流裂缝拟稳态常数b_Dpss计算新模型，定义无量纲产量积分q_Ddi和产量积分导数q_Ddid，所述重新定义无量纲产量为：

q_Dd＝b_Dpss·q_D

所述重新定义无量纲时间为：

所述无量纲产量积分表达式为：

所述无量纲产量积分导数表达式为：

其中，q_Dd为重新定义的无量纲产量，无因次；t_Dd为重新定义的无量纲时间，无因次；q_Ddi为无量纲产量积分，无因次；q_Ddid为产量积分导数，无因次；N_pDd为无量纲累计产量积分，无因次；r_eD为无量纲井控半径，无因次；t_D为无量纲，无因次。

通过上述定义式，并结合无量纲产量q_Dd和t_Dd的关系，进而可以绘制无量纲q_Ddi，q_Ddid和t_Dd曲线的关系曲线，即可绘制圆形封闭储层中垂直裂缝井产量递减新图版

所述绘制圆形封闭储层中垂直裂缝井产量递减新图版应用方法包括：

(1)时间数据处理：对于气井，计算物质平衡拟时间：

(2)产量数据处理：对于气井，计算拟压力规整化产量：

计算拟压力规整化产量积分：

计算拟压力规整化产量积分导数：

(3)在log～log双对数坐标图上分别绘制处理后的产量和处理后的时间的关系曲线，共计三组曲线：(q_g/Δp_p)～t_ca、(q_g/Δp_p)_i～t_ca、(q_g/Δp_p)_id～t_ca；

(4)同时将得到的三组曲线或者将其中两组任意组合，与所述的产量递减新图版进行拟合分析，使得每组曲线能够最大限度地获得较好的拟合效果；

(5)根据拟合结果记录无量纲圆形地层半径r_eD；

(6)在上述拟合曲线上选取任意一个拟合点，记录实际值和理论拟合值，即分别为(t_ca，q_g/Δp_p)、(t_Dd，q_Dd)；

(7)根据记录得到的C_fD和r_eD以及运用有限导流裂缝拟稳态常数b_Dpss计算新模型计算b_Dpss；

(8)综合以上结果利用绘制圆形封闭储层中垂直裂缝井产量递减新图版，计算油气田参数，原始地质储量：

井控面积、井控半径：

地层有效渗透率：

裂缝半长：

其中t_ca为物质平衡时间，d；μ_g为气体粘度，mPa·s；φ为孔隙度，％；C_g为气体压缩系数，MPa^-1；q_g表示产量，m³/d；G为地质储量，10⁴m³；N_p为累积产油量，10⁴m³；p_p为规整化拟压力，MPa；p_i为原始地层压力，MPa；p_wf为井底压力，MPa；q/△p规整化产量，m³/d/MPa；(q/△p)_i为规整化产量积分，m³/d/Mpa；(q/△p)_id为规整化产量积分导数，m³/d/Mpa；x_f为裂缝半长，m；r_e为圆形地层半径，m；r_eD为无量纲井控半径；S_wi为束缚水饱和度，％；k_g，k_o为油、气有效渗透率，mD；A为井控面积，m²；B_gi，B_oi为原始气、油体积系数；b_Dpss为稳态常数。

下面结合实际生产数据进一步阐述本发明的使用方法。某窄河道低渗透储层某气井储层、流体等相关的基本参数如下表所示：

表1 基本参数表

图6为A井产量数据q(半对数坐标)和历史压力数据p_wf(直角坐标)与时间的生产曲线图。从整体上看，图中数据点有分布规律明显且彼此有较好的相关性。如图7所示，产量递减新图版为三组双对数曲线，即规整化产量(q_g/Δp_p)～t_ca、规整化产量积分(q_g/Δp_p)_i～t_ca、规整化产量积分导数(q_g/Δp_p)_id～t_ca与理论曲线拟合分析图。

利用记录的C_fD和r_eD、时间和产量拟合数据可以计算相关静态或动态参数，例如，天然气原始地质储量G、有效渗透率k、井控面积A及半径r_e、裂缝半长x_f等。

由拟合结果：无量纲导流能力C_fD＝0.5，无量纲圆形地层半径r_eD＝1，时间拟合数据(t_ca/t_Dd)_M＝58.34(d)，产量拟合数据(q_g/Δp_p/q_Dd)_M＝1.3433×10³(m³/d/MPa)。

根据拟合方法计算以下参数：

计算拟稳态常数b_Dpss：

μ＝ln C_fD＝ln(0.5)＝-0.69315，I_x＝1/r_eD＝1/1＝1

计算天然气地质储量G：

计算井控面积A及半径r_e：

计算有效渗透率k，若使用法定单位制，则有：

计算裂缝半长x_f：

通过以上步骤利用该井的产量递减新图版根据拟合方法计算参数，计算参数如表2所示。

表2 某井计算参数汇总表

本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。