气液两相状态下的凝析气井产能评价新方法

陆家亮 张 皓 常宝华 曹 雯 孙贺东

1.中国石油勘探开发研究院 2.中国科学院大学 3.中国科学院渗流流体力学研究所

0 引言

凝析气藏既能采出天然气又能采出凝析油，是一种特殊类型的气藏。凝析气藏在开发过程中会发生天然气和凝析油系统复杂的相态变化，在一定温度和压力条件下气液两相共存并不断发生传热、传质，气液两相的体积分数也不断变化。

干气气田和凝析气田的开发最大的区别是地层流体的性质和相态特征不同。当凝析气藏地层压力降到露点压力以下时凝析油会从气相中析出，这些油残留和吸附在岩石颗粒的表面，在超过临界液体饱和度时开始流动，形成油气两相渗流[1]。此时若利用经典产能评价方法[2-5]对凝析气井的产能测试数据进行分析，经常出现二项式产能方程的系数B为负数的异常情况[6]，难以有效评价气井的产能，从而影响对气井生产动态的准确预测。

目前对于凝析气井产能试井的分析有两种思路：①将凝析油产量折算成凝析气产量，即将油气两相流体视为拟单相流体进行产能评价，该方法适用于高气液比情形；②针对地下的气液两相渗流，采用拟压力方法进行产能评价[7-10]，拟压力方法的计算过程相对复杂，需要考虑气液流体性质、相态及油气相渗等因素，但对于低气液比情形，该方法的计算结果更可靠。这两种方法都需进行产能测试，获取产能测试资料后才能进行分析。但对于高压、高产凝析油或产水的气井来说，不宜频繁进行开关井或改变气井的工作制度。在单相稳定点产能二项式方程的基础上[11]，陆家亮等[12]建立了气水两相稳定点的产能评价方法，实现了利用生产数据对产水气井进行产能的实时评价。

在前述研究的基础上，笔者建立了适用于凝析气井的拟单相稳定点产能评价方法（以下简称拟单相法）与气液两相稳定点产能评价方法（以下简称气液两相法）。

1 拟单相法

当井底压力高于露点压力时，地层中为单相气流；当井底压力低于露点压力时，地层近井带由单相气变为油气两相，若凝析油饱和度低于临界流动饱和度，气相流动而油相不流动，此时可认为地层中凝析油流速为零，仍可沿用单相气体渗流理论，只需将井口的凝析液量折算成相应的凝析气量，则有：

式中qt表示折算后的凝析气产量，104m3/d；qg表示井口气产量，104m3/d；qo表示井口凝析油产量，m3/d；γo表示凝析油的相对密度，无量纲。

对于圆形有界干气藏中一口直井，当气井的流动达到拟稳态阶段后，以压力平方表示的产量与压力的关系式为：

将式（1）代入到式（2），拟单相稳定点的二项式产能方程为[13]：

其中

式中pR表示地层压力，MPa；pwf表示井底压力，MPa；T表示储层温度，K； 表示气体平均偏差因子，无量纲；Tsc表示标准条件下温度，K；psc表示标准条件下压力，MPa； 表示天然气平均黏度，mPa·s；K表示储层渗透率，mD；h表示储层有效厚度，m；re表示边界半径，m；rw表示井筒半径，m；S表示表皮系数，无量纲；qt表示折算产气量，104m3/d；D表示非达西流系数，d/104m3。

2 气液两相法

当地层压力低于露点压力时，地层中就会出现油气两相，当凝析油饱和度高于临界流动饱和度时，油相也开始流动，此时所满足的渗流方程是多相渗流方程，引入拟压力函数ψ，即

油气两相渗流方程为：

其中

初始条件为：

定产内边界条件为：

其中

式中ψ表示拟压力，MPa；μgi表示原始条件下气相黏度，mPa·s；ρgi表示原始条件下气相密度，kg/m3；p0表示参考压力，MPa；p表示气藏中任一点压力，MPa；Kro、Krg表示油相、气相的相对渗透率，无量纲；ρo、ρg表示油相密度、气相密度，单位均为kg/m3；μo、μg表示油相、气相的黏度，单位均为mPa·s；r表示储层任一点至井底距离，m；η表示拟导压系数，mD·MPa/（mPa·s）；t表示时间，h； 表示气藏孔隙度，无量纲；So、Sg表示油相、气相的饱和度，无量纲；qtgo表示折算后的气、油两相地层条件下的质量之和，kg/d；ρgsc、ρosc表示气、油标准状态下的密度，单位均为kg/m3。

对于圆形有界封闭气藏中一口直井，当气井的流动达到拟稳定流动阶段后，气藏平均地层压力可以表示为：

计算拟压力函数时应注意：对于凝析气藏而言，当凝析油饱和度高于临界流动饱和度时，凝析油才会流动，此时的拟压力函数应由以下3个部分构成，即

式中pcr、pdew分别表示液相流动临界压力、露点压力，MPa。

相对渗透率和压力的关系可依据稳态理论进行计算，即

其中

式中L表示凝析油摩尔分数，无量纲；V表示凝析气摩尔分数，无量纲。

ρo、ρg、μo、μg、L 及 V 的大小随压力和温度而变化，可由相态闪蒸计算得到。

式（8）变形后可得气液两相稳定点的二项式产能方程为：

其中

3 理论分析与实例

3.1 理论计算与分析

国内较为典型的凝析气田包括牙哈气田、千米桥潜山气田、塔中Ⅰ号气田及迪那2气田等[14-17]，它们在储层类型、物性参数、相态特征等方面存在较大差异（表1）。在同一气田内部，初始的生产气油比变化范围也较大（表1）。

选用上述典型凝析气田的基础参数（表1），结合油气相渗曲线（采用上述凝析气田的油气相渗曲线，进行归一化处理后得到，如图1所示），设定井径为0.1 m，供气半径为1 000 m，表皮系数为10，应用上述产能方程计算不同气油比条件下凝析气井的无阻流量。采用牙哈气田的基础参数，设定井参数，计算凝析气井的IPR曲线。如图2所示，生产气油比R为1 500 m3/m3，当地层压力高于露点压力时，气液两相法与拟单相法所计算的无阻流量相差10.47%；当地层压力低于露点压力后，两种方法所计算的无阻流量相差30.42%。因此，当凝析气藏储层中存在气液两相流动时，无阻流量的评价应以气液两相法计算的结果为准。

结合4个气田的典型参数，采用拟单相法与气液两相法计算凝析气井在不同气油比情形下的无阻流量，统计两种方法计算的无阻流量的差异，即以qAOF1代表采用拟单相法计算的无阻流量，qAOF2代表采用气液两相法计算的无阻流量，计算（qAOF1－qAOF2）与qAOF2的比值。如图3所示，随着生产气油比增大，qAOF1与qAOF2的差异逐渐减小，即在气油比较高的情况下拟单相法与气液两相法计算的无阻流量结果相差较小；当气油比较低时，储层中发生气液两相流动，气相相对渗透率降低，造成无阻流量降低，采用气液两相法计算无阻流量，结果更可靠。

表1 典型凝析气田基础参数表

图1 油气相渗曲线图

图2 露点压力以上、以下的气井IPR曲线图（采用牙哈气田的基础参数）

图3 拟单相法与气液两相法计算无阻流量的差异曲线图

3.2 典型凝析气井的实例分析

牙哈气田的储层为中孔、高渗储层，非均质程度弱、连通性好，凝析油含量高，投产初期气井的生产气油比在300～2 500 m3/m3范围内。牙哈气田A井于2004年投产，初期产气量稳定在（20～30）×104m3/d之间，平均生产气油比为1 800 m3/m3，初始地层压力为57 MPa，投产两个月后进行回压试井，测试数据如表2所示。

应用经典产能评价方法分析上述测试数据，二项式产能方程系数为负数，无法计算无阻流量。采用拟单相法和气液两相法计算该井IPR曲线，如图4所示，无阻流量分别为142.14×104m3/d和102.52×104m3/d，两者相差38.65%，由前述敏感性分析可认为，无阻流量为102.52×104m3/d更可靠。

表2 牙哈气田A井回压试井测试数据表

图4 牙哈气田A井IPR曲线图

千米桥潜山气田主体为中等凝析油含量的饱和型凝析气藏。千米桥潜山气田B井于1999年投产，产量、油压均呈逐渐递减的变化趋势，平均生产气油比为3 000 m3/m3，原始地层压力为43.55 MPa。取8 mm油嘴的测试数据：产油量为56 m3/d、产气量为17×104m3/d、井底流压为37.76 MPa，应用拟单相法和气液两相法计算该井的IPR曲线，如图5所示，无阻流量分别为77.26×104m3/d和60.33×104m3/d，相差28.06%。

图5 千米桥潜山气田B井IPR曲线图

塔中Ⅰ号气田的储层与千米桥潜山气田的储层类似，但洞穴、裂缝更发育，凝析油含量范围广，生产气油比变化较大，经典产能评价方法的适用性差。塔中Ⅰ号气田C井于2010年投产，油、气产量及油压均呈递减的变化趋势，初期平均生产气油比为6 000 m3/m3，初始地层压力为58 MPa。2011年该井的测试数据如下：产油量为10.5 m3/d、产气量为6.3×104m3/d、井底流压为36.24 MPa。采用拟单相法计算的无阻流量为45.62×104m3/d，采用气液两相法计算的无阻流量为41.72×104m3/d，相差9.35%。

迪那2气田的储层为低孔、裂缝性储层，生产气油比相对较高。迪那2气田D井于2009年投产，初始地层压力为106 MPa。2017年该井的测试数据如下：产油量为85.2 m3/d、产气量为72.5×104m3/d、井底流压为83.1 MPa，采用拟单相法计算该井的无阻流量为193.72×104m3/d，采用气液两相法计算该井的无阻流量为181.54×104m3/d，相差6.71%。

如表3所示，4口典型凝析气井无阻流量的计算结果表明：随着生产气油比增大，采用拟单相法与气液两相法计算的结果差异逐渐减小，与理论计算的结果一致。

4 结论

1）采用拟单相法和气液两相法计算凝析气井的无阻流量，避免了采用经典产能评价方法出现产能方程系数为负数而无法计算无阻流量的情况。

表3 典型凝析气井采用拟单相法及气液两相法计算的无阻流量结果表

2）当地层压力高于露点压力时，储层流体以单相流体为主，可以采用拟单相法；当地层压力低于露点压力时，采用拟单相法计算的无阻流量会偏高。

3）当生产气油比较高时，采用拟单相法和气液两相法计算的无阻流量差异较小，且生产气油比越高，差异越小；当生产气油比较低时，采用气液两相法计算的无阻流量比采用拟单相法计算的结果更可靠。

[ 1 ] 李海平, 任东, 郭平, 陈建军. 气藏工程手册[M]. 北京: 石油工业出版社, 2016.Li Haiping, Ren Dong, Guo Ping & Chen Jianjun. Gas reservoir engineering handbook[M]. Beijing: Petroleum Industry Press,2016.

[ 2 ] Pierce HR & Rawlins EL. The study of a fundamental basis for controlling and gauging natural-gas wells, Part 1—computing the pressure at the sand in a gas well: Rept. of Investigations of 2929[R]. Washington DC: US Bureau of Mines, 1929.

[ 3 ] Rawlins EL & Schellhardt MA. Back-pressure data on natural-gas wells and their application to production practices[R]. Washington DC: US Department of the Interior Bureau of Mines, 1935.

[ 4 ] Cullender MH. The isochronal performance methods of determining the flow characteristics of gas wells[J]. Transactions of the AIME, 1955, 204: 137-142.

[ 5 ] Katz DL, Cornell D, Kobayashi R, Poettmann FH, Vary JA, Elenbaas JR, et al. Handbook of natural gas engineering[M]. New York: McGraw-Hill Book Co., Inc., 1959.

[ 6 ] 江同文, 朱忠谦, 熊钰, 阳建平. 凝析气藏循环注气开发动态分析方法与实践[M]. 北京: 石油工业出版社, 2013.Jiang Tongwen, Zhu Zhongqian, Xiong Yu & Yang Jianping. Dynamic analysis method and practice of circulating gas injection in condensate gas reservoir[M]. Beijing: Petroleum Industry Press,2013.

[ 7 ] 何志雄, 孙雷, 李士伦. 凝析气井产能分析新方法[J]. 天然气工业, 1996, 16(5): 32-36.He Zhixiong, Sun Lei & Li Shilun. New method of condensate gas well productivity analysis[J]. Natural Gas Industry, 1996,16(5): 32-36.

[ 8 ] 谢兴礼, 罗凯, 宋文杰. 凝析气新的产能方程研究[J]. 石油学报, 2001, 22(3): 36-42.Xie Xingli, Luo Kai & Song Wenjie. A novel equation for modeling gas condensate well deliverability[J]. Acta Petrolei Sinica,2001, 22(3): 36-42.

[ 9 ] 覃斌, 李相方, 程时清. 考虑气液两相流的凝析气井生产动态研究[J]. 新疆石油地质, 2004, 25(4): 423-426.Qin Bin, Li Xiangfang & Cheng Shiqing. Production performance of condensate gas wells considering two phase fl uid of gas and liquid[J]. Xinjiang Petroleum Geology, 2004, 25(4): 423-426.

[10] 廖发明, 苗继军, 陈文龙, 邓军, 王怒涛. 凝析气井产能和储量计算新方法[J]. 西南石油大学学报(自然科学版), 2012,34(4): 100-104.Liao Faming, Miao Jijun, Chen Wenlong, Deng Jun & Wang Nutao. The new calculation method of condensate gas well production and reservoirs[J]. Journal of Southwest Petroleum University(Science & Technology Edition), 2012, 34(4): 100-104.

[11] 庄惠农. 气藏动态描述和试井[M]. 2版. 北京: 石油工业出版社, 2009.Zhuang Huinong. Dynamic gas reservoir description and well test[M]. 2nded. Beijing: Petroleum Industry Press, 2009.

[12] 陆家亮, 涂梅桂, 刘晓华, 肖前华. 产水气井稳定点二项式产能方程及其应用[J]. 科学技术与工程, 2016, 16(6): 66-69.Lu Jialiang, Tu Meigui, Liu Xiaohua & Xiao Qianhua. Stable point binomial productivity equation of water-bearing gas well and application[J]. Science and Technology and Engineering,2016, 16(6): 66-69.

[13] 张皓, 陆家亮, 曹雯, 孙贺东. 不同计量单位下的气井产能试井常用公式[J]. 油气井测试, 2017, 26(6): 1-7.Zhang Hao, Lu Jialiang, Cao Wen & Sun Hedong. Usual productivity equation of gas well on various units of measurement[J].Well Testing, 2017, 26(6): 1-7.

[14] 李熙喆, 郭振华, 胡勇, 罗瑞兰, 苏云河, 孙贺东, 等. 中国超深层构造型大气田高效开发策略[J]. 石油勘探与开发, 2018,45(1): 111-118.Li Xizhe, Guo Zhenhua, Hu Yong, Luo Ruilan, Su Yunhe, Sun Hedong, et al. Efficient development strategies for large ultra-deep structural gas fi elds in China[J]. Petroleum Exploration and Development, 2018, 45(1): 111-118.

[15] 常志强, 肖香姣, 朱忠谦, 王小培, 王海应. 高温高压气井、凝析气井井筒及近井地带的压力变化规律[J]. 天然气工业,2012, 32(10): 47-51.Chang Zhiqiang, Xiao Xiangjiao, Zhu Zhongqian, Wang Xiaopei & Wang Haiying. Pressure change law at the wellbore and well-adjacent zone of a HTHP or condensate gas well[J]. Natural Gas Industry, 2012, 32(10): 47-51.

[16] 康博, 张烈辉, 王健, 范坤, 王宏. 裂缝—孔洞型碳酸盐岩凝析气井出水特征及预测[J]. 西南石油大学学报(自然科学版),2017, 39(1): 107-113.Kang Bo, Zhang Liehui, Wang Jian, Fan Kun & Wang Hong. Features and forecast of water output in fractured vuggy carbonate condensate reservoir[J]. Journal of Southwest Petroleum University (Science & Technology Edition), 2017, 39(1): 107-113.

[17] 马小明, 徐元洲, 李云, 朱杰, 张玉梅. 千米桥古潜山凝析气藏油气富集规律新认识[J]. 天然气工业, 2009, 29(4): 25-26.Ma Xiaoming, Xu Yuanzhou, Li Yun, Zhu Jie & Zhang Yumei.Research on hydrocarbon enrichment in condensate gas reservoirs in Qianmiqiao buried hill[J]. Natural Gas Industry, 2009, 29(4):25-26.