确定产水气井流入动态关系的新方法

摘 要

摘　要：当地层为单相气体渗流时，气井二项式产能方程系数A、B中的渗透率为地层渗透率；当地层为气水两相流动时，则系数A、B中的渗透率为气相渗透率。为确定产水气井的产能方程，基于

摘　要：当地层为单相气体渗流时，气井二项式产能方程系数A、B中的渗透率为地层渗透率；当地层为气水两相流动时，则系数A、B中的渗透率为气相渗透率。为确定产水气井的产能方程，基于气体地下稳态渗流理论建立的产能方程，在没有产能试井资料的情况下，提出了利用产水气井生产动态参数确定产水气井流入动态关系的新方法。该新方法主要包括：首先根据垂直管气水两相流压力计算方法，结合产水气井动态参数，确定对应的井底流压；进而利用初期不产水时的不稳定试井解释成果，结合目前的生产动态资料及压力资料，确定对应的气相渗透率；最后，根据确定的气相渗透率，结合产水气井气相二项式产能方程，确定产水气井流入动态关系。利用该新方法对产水气井进行了实例计算，分别获得了不产水及产水情况下的二项式产能方程。计算结果表明，气井产水使得二项式产能方程系数A、B变大，气井无阻流量大大降低。

关键词：气井  产水  流入动态  气水两相流  相渗透率  产能方程  无阻流量  数学模型

A new method of determining the inflow performance relationship in water-producing gas wells

Abstract：When it comes to a single phase flow in the strata，the permeability in binomial coefficients A and B of binomial deliver ability equatlon for a gas well stands for formation permeability，but when it comes to a gas water two phase flow，it turns to be gas phase permeability,Based on the deliverability equation established upon the steady seepage theory of gas underground，and in the absence ot dellVerabllltY test data for a gas well，this paper presented a new method of determining inflow performance relationshid of water-producing gas wells by using dynamic production parameters of water-producing wells．This new method mainly includes the following contents：First,the corresponding bottomhole flowing pressure is determined according to the gas&water flowing pressure calculatlon method ot a vertlcal pipe and dynamic parameters of a water producing gas well．Second，the corresponding gas permeability ls determlned according to transient well test results while non-water-producing at the preliminary stage in combination with dynamlc Productlon parameters and pressure data at the current stage．Finally，the inflow performance relationshiD of a water producing gas well ls determlned according to the calculated gas permeability and binomial deliverability equation for a water-producing gas well.This method was put into practice for the calculation of a water-producing well，and the binomial deliverability eauations were obtained for a gas well with and without water production respectively．The calculation results show that when there is water Produced in a gas well，the coefficients A and B in the binomial deliverability equation for the well are made greater and its open ilow capacity is remarkably reduced．This study provides a new idea for determining the deliverability equation of a water-producing gas well．

Key words：gas well，water producing，inflow performance，gas-water two-phase flow，phase permeability，deliverability eauation．open flow capacity，mathematical model

气井产能方程的确定方法[1-2]一通常分为由产能试井分析获得产能方程与根据渗流理论推导出产能方程表达式并确定产能方程两大类[3]。迄今为止，气水两相渗流的研究相对较少，郝玉鸿等[4]提出了利用复合模型对地层产水时气井的产能方程进行研究；李晓平等[5-6]从渗流力学的基本原理出发，推导出气水两相流井的二项式产能方程，做出了不同气水比影响下的稳态及瞬态流入动态关系曲线；贺遵义[8]建立了气液两相渗流数学模型，提出了用一点法进行动态计算的无因次方程和无因次曲线图版；姜必武等[8]考虑表皮系数和非达西效应的影响，建立了气水两相渗流的产能方程；常彦荣等[9]。基于气水井的流入动态研究成果，将其与流出动态相耦合，探讨了气水井的流入流出动态关系；朱光亚等[10]在考虑低渗透气藏的气、水渗流规律不遵循达西定律的情况下，推导得到了低渗透砂岩气藏气水两相稳态径向渗流问题的半解析解。另一些学者进行了气水两相渗流的实验研究，胡勇等[11]通过CT成像技术和核磁共振技术，研究了火山岩储层中气水两相渗流的微观机理，研究发现，即使在储层含水饱和度较高的情况下，气相依然具有一定的渗流能力，但是气相渗透率随含水饱和度的增加而迅速下降；鄢友军等[12]制作了激光刻蚀微观鲕粒模型，进行了气水两相驱替实验，实验研究表明：对于鲕粒模型，水驱气时卡断、绕流是形成封闭气的主要原因。

笔者基于气体地下稳态渗流理论建立的产能方程，提出了利用产水气井生产动态参数确定产水气井产能方程的理论与方法，并用l口实例气井进行了计算，分别获得了不产水及产水情况下的流入动态关系曲线，这为产水气井产能方程的确定提供了一种新思路。

1　产水气井产能方程

气井单相渗流的稳定二项式产能方程为：

 

其中

 

式中PR为地层压力，MPa；Pwf为井底流压，MPa；A、B为二项式系数；qsc为日产气量，m3／d；T为地层温度，K；`m为平均黏度，mPa·S；`Z为平均偏差因子；K为绝对渗透率，mD；h为地层有效厚度，m；re为供给半径，m；rw为井底半径，m；S为表皮系数；b为湍流引起的惯性阻力系数；ge为天然气的相对密度。

当气井产水，地层之中为气水两相流之时，对气相而言，式(1)中的参数渗透率应该变成气相渗透率。即对于气水两相流动而言，气井稳定状态条件下以压力平方形式表示的气相二项式产能方程为：

 

 

其中

 

式中A1、B1为气水两相流时气相二项式系数；Kg为气相渗透率，mD；b1为气水两相流时湍流引起的惯性阻力系数。

在没有产能试井资料的情况下，利用曰前地层压力、井口油压、相应的日产气量、日产水量、初期不稳定试井资料解释所获得的参数，确定气水两相流动时气相二项式产能方程式(2)的步骤为：①根据井口油压、相应的口产气量、日产水量及相关参数，由气水两相流井筒模型，确定对应的井底流压；②根据目前的地层压力，计算所得井底流压、日产气量及相关参数，结合式(2)可求得对应的气相渗透率；③将求得的气相渗透率代入即呵求得二项式产能系数A1、B1，从而可确定目前地层压力下气水两相流动时的气相二项式产能方程。

2　井底流压的计算

当井筒之中存在气水两相时，根据井VI油压、对应压力下的日产气量、日产水量及相关参数，在确定对应井底流压的过程之中，用Hagedorn Brown模型来计算气液两相流的井底压力优于其他模型[13-16]。

与单相流动相同，根据垂直管流中流体的质量守恒与动量守恒原理可得：

 

从式(3)出发，考虑到现场实际单位的换算，可以得到：

 

其中

 

 

对于不产油的产水气井而言，M1、r1可简化为：

 

式中P为井筒中任意点的压力，MPa；D为井深，m；rm为气液混合物的密度，kg／m3；um为气液混合物的流速，m／s；g为重力加速度，m／s2；fm为两相摩阻系数；d为油管内径，m；qt为地面产液量，m3／d；Mt为地面标况下，每生产1m3液体伴生油、气、水的总质量，kg／m3；H1为持液率；r1为液体密度，kg／m3；go、gg、gw分别为油、气、水的相对密度；Rs为天然气的溶解气油比；Bo、Bw分别为油、水的体积系数；WOR为生产水油比；rw、rg分别为水、气的密度，kf／m3；¦w、¦o分别为所产液体中的含水率、含油率；GLR为气液比。

在选用此模型进行计算的时候主要采用了迭代的方法。首先选定井筒长度增量(AD)，然后根据公式：

 

首先在相应变化的长度区间内假设一个DP，计算出在此长度和压力区间内的物性参数，利用式(7)计算出DP1。将DP1与假设的DP进行比较，如果满足精度要求则进行下一个长度区间的计算，否则令DP＝DP1，再次进行计算，如此反复，直到满足要求为止。

3　气相渗透率计算

由气水两相流动时，气井稳定状态条件下以压力平方形式表示的气相二项式产能方程式：

 

令

 

将其代换后可得：

 

由式(9)解得气相渗透率的表达式为：

 

由式(10)即可确定目前地层压力、井底流压、相应日产气量、日产水量条件下的气相渗透率。

将计算所得的气相渗透率代入相应式中，即可求得二项式产能系数A1、B1，从而可确定目前地层压力下气水两相流动时气相的二项式产能方程。

4　实例计算

某水驱气藏中1口气井的实际生产资料、测压资料、初期不产水时不稳定试井解释成果、井身资料和实验室数据如下：井半径0.1m，油管直径88.9mm，产层厚度35.5m，地层中部温度405.49K，气体相对密度0.574，供给边界半径720m，表皮系数1.98，气井深度4887.7m，渗透率l.356mD，目前地层压力55.12MPa，产气量5.6×104m3／d，产水量l8.9m3／d，井口油压31.9MPa，平均偏差因子1.24，气体平均黏度0.0365mPa·S。根据上述理论计算可得，目前生产条件下的井底流压为48.4MPa，气相渗透率为0.59mD。考虑水影响时，二项式产能方程系数A1为0.012365、Bl为1.10×10-9、无阻流量为24.06×104m3／d；不考虑水影响时，二项式产能方程系数A为0.005353、B为3.14×10-10、无阻流量为54.99×104m3／d。计算结果表明，地层中的气水两相流动使得气井的无阻流量下降了56.25％。考虑及不考虑水影响的气井流入动态关系曲线如图1所示。

 

由图l可以看出，不考虑水影响时计算的无阻流量要远远大于考虑水影响时计算所得的无阻流量。根据单相气流的渗流理论与气水两相的渗流理论分析可知，此计算结果是可信的。这也表明气井的产水将大大地降低气井的产能。

5　结论

1)笔者提出了利用产水气井生产动态参数确定产水气井产能方程的理论与方法，利用该方法对l口实例产水气井进行了计算，分别获得了不产水及产水情况下的二项式产能方程，本研究为产水气井产能方程的确定提供了一种新思路，为气井合理工作制度的确定提供了依据。

2)研究表明，气井产水使得地层中的流动由单相气流过渡到气水两相流，渗流阻力增大。二项式产能方程系数A、B均变大，气井无阻流量大大降低。

 

参考文献

[1]A1-ATTAR H，Al-ZUHAIR S．A general approach for deliverability calculations of gas wells[J]．Journal of Petroleum Science and Engineering，2009，67(3)：97-104．

[2]CHASE R W，ANTHONY T M．A simplified method for determining gas-well deliverability[J]．SPE Reservoir EngineerinG，1988，3(3)：1090-1096．

[3]庄惠农．气藏动态描述和试井[M]．北京：石油工业出版社，2004：52-69．

ZHUANG Huinong．Dynamic description and well test of gas reservoir[M]．Beijing：Petroleum Industry Press，2004：52-69．

[4]郝玉鸿，李治平．地层产水影响气井产能的定量分析[J]．油气井测试，l998，7(4)：ll-13．

HAO Yuhong，LI Zhiping．A quantitative analysis of for mation water affecting gas well productivity[J]．Well Tesring，1998，7(4)：11-13．

[5]李晓平，赵必荣．气水两相流井产能分析方法研究[J]．油气井测试，2001，10(4)：8-10．

LI Xiaoping，ZHAO Birong．Study of productivity analysis method for gas water well[J]．Well Testing，2001，10(4)：8-10．

[6]李晓平，胡勇．气水同产井瞬态流入动态关系曲线探讨[J]．天然气工业，2001，21(3)：65-67．

LI Xiaoping，HU Yong．A discussion on the transient in flow performance relationship curve of gas water well[J]．Natural Gas Industry，2001，2l(3)：65-67．

[7]贺遵义．气水两相流入井的一种数理模型及一点法求解[J]．钻采工艺，2002，25(5)：28-31．

HE Zunyi．A mathematical model of gas-water percolation and its calculation method[J]．Drilling&Production Technology，2002，25(5)：28-31．

[8]姜必武，丘陵．含水气藏合理产能新方法研究[J]．天然气工业，2005，25(12)：80-82．

JIANG Biwu，QIU Ling．New way to determine reasonable productivity of water bearing gas reservoir[J]．Natural Gas Industry，2005，25(12)：80-82．

[9]常彦荣，李允，李晓平，等．气水同产井流入流出动态关系曲线研究[J]．西南石油学院学报，2006，28(3)：11-12．

CHANG Yanrong，LI Yun，LI Xiaoping，et al．Research of production nodal analysis method for gas-water wells[J]．Journal of Southwest Petroleum Institute，2006，28(3)：11-l2．

[10]朱光亚，刘先贵，高树生，等．低渗透气藏气水两相渗流模型及其产能分析[J]．天然气工业，2009，29(9)：67-70．

ZHU Guangya，LIU Xiangui，GAO Shusheng，et al．A gas-water two-phase flow model for low permeability gas reservoirs and deliverability analysis[J]．Natural Gas Industry，2009，29(9)：67-70．

[11]胡勇，朱华银，万玉金，等．大庆火山岩孔隙结构及气水渗流特征[J]．西南石油学院学报，2007，29(5)：63-65．

HU Yong，ZHU Huayin，WAN Yujin，et ai．The porestructure and the properties of water-gas flowing through porous media in a volcanic rock in Daqing[J]．Journal of Southwest Petroleum Institute，2007，29(6)：63-65．

[12]鄢友军，陈俊宇，郭静姝，等．龙岗地区储层微观鲕粒模型气水两相渗流町视化实验及分析[J]．天然气工业，2012，32(1)：64-66．

YAN Youjun，CHEN Junyu，GUO Jingshu，et al．A visualized experiment on gas-water two-phase seepage through oolitic reservoirs in the Longgang gas field，Sichuan Basin[J]．Natural Gas Industry，2012，32(1)：64-66．

[13]HAGEDORN A，BROWN K．Experimental study of pressure gradients occurring during continuous two phase flow in small diameter vertical conduits[J]．Journal of Petroleum Technology，1965，17(4)：475-484．

[14]ORKISZEWSKI J．Predicting two phase pressure drop in vertical pipe[J]．Journal of Petroleum Technology，1967，19(6)：829-838．

[15]BEGGS H D．Production optimization using nodal analysis[M]．1st ed．Tulsa：Oil＆Gas Consultants International Publication，l991.

[16]HASAN A R，KABIR C S．Aspects of wellbore heat transfer during two phase flow[J]．Oil Production&Facilities，1994，9(3)：211-216．

 

本文作者：李晓平  胡俊坤  王玉忠  张烈辉

作者单位：“油气藏地质及开发工程”国家重点实验室·西南石油大学

  西南石油大学研究生学院

  中国石油华北油田公司第二采油厂