一个新的考虑启动压力梯度的三项式产能方程

一个新的考虑启动压力梯度的三项式产能方程1

梁斌，张烈辉，刘启国，朱仲义

“油气藏地质及开发工程”国家重点实验室，西南石油大学，成都 (610500 )

E-mail: liangbinswpi@

摘 要：当用常规气藏气井产能方程计算低渗和特低渗气藏的产能时，得到的结果与实际情

况差别很大，甚至不能分析气井产能。其原因是，低渗和特低渗气藏存在启动压力梯度，而常规气藏产能公式没有考虑启动压力梯度，因而产生谬误。文章在Forchheimer 提出的二项式非达西渗流方程中添加了启动压力梯度项，通过数学换元推导出了考虑启动压力梯度的三项式产能方程。与常规二项式产能方程不同的是，三项式产能方程中存在一项启动压力梯度校正项。文章采用试算法求取启动压力梯度校正项，进而求出产能方程。实例分析表明：用常规产能方程分析低渗气藏气井的产能时，得到的产能系数b 小于0；但用新的产能方程得到的(P e 2- P wf 2-c)/q sc 与q sc 关系直线的相关性很好，产能系数都大于0，因而能评价气井产能。另外，通过试算出的启动压力梯度修正系数c ，可以求出启动压力梯度。 关键词：非达西渗流；启动压力梯度；产能方程；低渗气藏

1. 引言

目前我国所发现的很多气藏为低渗和特低渗气藏，气体在此类气藏中流动时，必须克服启动压差后才能流动。当对这类气藏气井进行产能分析时发现，用常规气藏产能公式计算出的产能与实际产能差别很大，甚至得到产能方程系数b<0。究其原因是，常规气藏气井产能公式没有考虑启动压力梯度。因而当评价此类低渗气藏气井的产能时，需要对常规气藏气井产能公式进行修正。

目前，国内很多人研究了低渗透气藏气井产能，并提出了一些比较有价值的新的理论和结论。李治平等人[1]基于Forcheimer 高速渗流方程，建立一个考虑启动压力梯度的形式简洁的三项式产能方程，并用最优化分析方法进行了计算。与常规产能方程相比，其产能方程多了一个修正项D ，以修正启动压力梯度对产能的影响。但通过仔细分析发现，修正项D 的物理意义并不明确，似乎还需要进一步研究。严文德等人[3]通过一个包含启动压力梯度的数学代换，定义了一个考虑启动压力梯度的视表皮系数，得到了一个新的低渗气藏气井的产能方程。但当用其产能方程进行计算时，需要计算一个净产压力平方差以及视表皮系数，这两个参数都与启动压力梯度有关，从而限制了该产能方程在实际中的应用。

2. 新的三项式产能方程

基于Forchheimer 渗流规律的常规稳态产能方程：

2

2

2

sc sc wf e bq aq p p +=− (1) 该产能方程只考虑了粘滞力和惯性力，没有考虑启动压力梯度的影响，不能分析存在启动压力梯度的低渗气藏气井的产能。

考虑启动压力梯度的高速非达西渗流规律：

B v v k

dr dp λβρµ

＋2+= (2)

1

本课题得到教育部博士点基金项目（20040615004）、四川省学术和技术带头人培养基金项目（2200320）和高等学校优秀青年教师科研奖励计划（TRAPOYT ）的资助。

令： B dr

dp

dr dp λ−=' (3)

式(2)可以简化为：

2'v v k

dr dp βρµ

+= (4) 将渗流速度rh

q

v π2=，气体密度和拟压力代入式(4)进行积分，并写为压力平方的形式可以得到：

2

2

'2

'sc sc wf e bq aq P p +=− (5) 假定启动压力梯度为常数，将式(3)从w r 到e r 进行积分，可以得到：

)('

'w e B wf e wf e r r p p p p −−−=−λ (6)

由于边界处的气体还来不及流动，不存在由启动压力梯度造成的附加压降，所以

e e p p ='，由式(6)可以得到：

)('

w e B wf wf r r p p −+=λ (7)

将'wf p 和'

e p 代入式(5)，可以得到：

2

2222])()(2[sc sc w e B w e B wf wf e bq aq r r r r p P p +=−+−−−λλ (8)

令： 2

2)()(2w e B w e B wf r r r r p c −+−=λλ (9) 式(9)表明，c 与启动压力梯度和井底流压有关，它反映了启动压力梯度对气井产能的影响。将式(9)代入式(8)，则得到一个新的考虑启动压力梯度的产能方程：

2

22sc sc wf e bq aq c P p +=−− (10)

3. 新的产能方程的求解

可以采用试算法求解启动压力校正系数c 和产能方程。具体求解步骤是：假定一个c 值，将(P e 2- P wf 2-c)/q sc 与q sc 在直角坐标系中做线性回归，求直线的斜率b 和截距a 。若a>0和b>0，且直线的相关性足够大，则c 值为所求。

根据试算的a 、b 和c 值，由式(10)求出气井无阻流量： b

c p b a a Q e AOF 2)

(422−++−=

(11)

4. 实例分析

为了验证本文提出的三项式产能方程，用常规产能分析方法和本文提出的产能分析方法对我国某特低渗气藏气井产能进行了验算。气井原始地层压力10.64MPa ，地层温度312.5K ，有效厚度8m ，孔隙度14%，含水饱和度60%，井筒半径0.059m ，天然气相对密度0.6028，稳定延时生产240小时，稳定延时产量和井底压力分别为0.79万方/天和4.86MPa ，井距为1000m 。产能试井数据如表1所示。