摘 要

摘　要：在异常高压气藏开发过程中，储层岩石常常表现出应力敏感特征。由于压敏储层在有效应力变化范围较大时，渗透率模量将随着有效应力的变化而变化，这样基于常渗透率模量的产能

摘　要：在异常高压气藏开发过程中，储层岩石常常表现出应力敏感特征。由于压敏储层在有效应力变化范围较大时，渗透率模量将随着有效应力的变化而变化，这样基于常渗透率模量的产能方程很难准确评价异常高压气藏的产能特征。为此，针对异常高压气藏的渗流规律，建立了考虑变渗透率模量和高速非达西效应的产能计算模型，并用Newton Raphson方法对模型进行求解。研究及实例分析结果表明：应力敏感作用会使异常高压气藏气井的产量降低，尤其是在井底流压较低时，其产量降低得更为严重；在井底流压较小时，用传统的常渗透率模量计算的产量可能会错误评价应力敏感作用对气井产量的影响；而高速非达西效应将会使异常高压气藏气井的产量降低，随着井底流压的减小，高速非达西效应使气井产量降低的程度更为严重。结论认为：在开发异常高压气藏的过程中应考虑变渗透率模量和高速非达西效应对气井产能的影响。

关键词：异常高压气藏  变渗透率模量  生产能力  应力敏感  高速非达西流动  数学模型

A new method for calculating the productivity of abnormally high-pressure gas reservoirs considering variable permeability modulus

Abstract：In the development of abnormally high pressure gas reservoirs，the reservoir rock generally exhibits the characteristic of stress sensitivity．For such stress sensitive reservoirs，the permeability modulus changes with the stress if the effective stress is of a large Variation range．Therefore，the productivity fornlula based on the constant permeability roodulus is unlikely to accurately evaluate the productivity of such gas reservoirs．In view of this，a productivity calculation roodel was established considering both the variable permeability model and high velocity non Darcy flow effect according to the permeability rule of abnormallv high-pressure gas reservoirs；The Newton-Raphson method was adopted for solving the model．The study resuhs and case analysis showed that the stress sensltlvlty has the effect of reducing the production of abnormally high pressure gas reservoirs．Particularlv，the production decreases heavily when the FBHP(flowing bottomhole pressure)is low．Thus if the traditional calculation method based on the constant permeability modulus is adopted，the evaluated productivity will probably provide unreliable basis for assessing the stress-sensitive effect on the gas well production.Moreover,the high velocity non-Darcy flow effect will reduce the productivity of abnormally high pressure gas reservoirs，and the reduction becomes severer with the decrease of FBHP．Therefore，the effect of such factors as the variable permeability modulus and the high velocity non-Darcy flow effect on well productivity should be considered in thc development of abnormally high-pressure gas reservoirs．

Key words：overpressured gas reservoir，variable permeability modulus，productivity，stress sensitivity，high velocity non Darcy

异常高压气藏作为一种特殊气藏有着广泛的分布，近年来随着常规天然气藏资源的减少以及油气田勘探开发技术的发展，我国相继发现了越来越多的异常高压气藏，如位于塔里木盆地的克拉2气藏和四川盆地的磨溪嘉二段气藏。根据异常高压气藏开发的调研可知^[1-8]，与常规气藏相比，异常高压气藏具有驱动能量大、储层岩石变形对气井产能影响显著等特征。因此对于异常高压气藏的开发需考虑高速非达西效应和应力敏感作用对其产能的影响。目前已有不少学者针对异常高压气藏的特点研究了考虑应力敏感的产能计算模型^[9-14]，然而在描述渗透率随有效应力变化关系时常常假设渗透率模量为常数，这样在气藏开发过程中有效应力变化范围较大时，常渗透率模量就很难准确表征储层的应力敏感特征。虽然已有部分学者认识到用变渗透率模量来描述应力敏感特征更为合理^[15-20]，但是目前却没有见到针对异常高压气藏特征的考虑变渗透率模量影响的产能计算方法。因此笔者针对异常高压气藏的特征，建立了考虑变渗透率模量和高速非达西效应影响的产能模型，并用Newton Raphson方法^[21]对模型进行了求解，分析了变渗透率模量对于异常高压气藏产能的影响，旨在为开发异常高压气藏提供一种新的产能评价手段，以便更加合理高效地开发异常高压气藏。

1　岩石变渗透率模量特征

在异常高压气藏的开采过程中，由于应力敏感作用，储层岩石的渗透率随地层压力的变化关系通常表示为^[22-23]：

K＝K0exp[-a(p0-p)]                 (1)

通常渗透率模量被视为常数，但是目前国内外不少学者通过大量实验发现：当地层压力变化范围较大时，1n(K／K0)与p0-p关系曲线(图1)并非呈现单一的线性关系，而表现为渗透率模量随着有效应力的增加而逐渐变小，具有明显的变渗透率模量特征^[15-20]。为了描述这种变渗透率模量特征，通常采用分段渗透率模量^[15-18]来表征渗透率随有效应力变化规律(图1)。即假定异常高压气藏原始地层压力为p0，原始渗透率为K0则从P0到P的压力变化过程可以根据渗透率模量的不同而划分为m个阶段pi-1—pi(i=1，2，…，m)，每个阶段的渗透率模量为ai(i=1，2，…，m)。

 

2　产能模型的建立及求解

假设有一个圆形外边界的等厚异常高压气藏，一口气井位于其中心并以定产量Qsc生产体在储层中的流动为平面径向流，异常高压气藏原始地层压力为P0，原始渗透率为K0，从气藏供给压力P0到井底流压加(或p0)的压力变化过程可以根据渗透率模量的不同而划分为m个阶段pi-1—pi(i=1，2，…，m)，每个阶段的渗透率模量保持不变，渗透率模量值为ai(pi≤p≤pi-1，i=1，2，…，m)。

根据渗流力学中的稳态平面径向渗流理论可以知道，储层中的压力分布是以井眼为中心的同心圆环。则根据渗透率模量的不同，从井筒半径rm(或rw)到供给半径r0(或re)也可以划分为m段，其分界点依次为rm，rm-1，…，r1，r0，分界点处对应的压力和渗透率分别为pm, pm-1，…，p1，p0，和Km，Km-1，…，K1，K0，其中：

 

由于异常高压气藏驱动能量较大，气体在地层中的流速较高，因此流体的流动不再符合达西定律，而满足Forchheimier二项式方程^[24]：

 

在平面径向渗流过程中气体通过径向距离为r时的渗流速度为：

 

地层条件下的气体密度为：

 

描述孔隙介质的湍流系数为^[25]：

 

在Pi≤P≤i-1内(其对应的径向距离ri≤r≤ri-1)，储层渗透率可以表示为：

 

将式(3)、(4)、(5)和(6)代入式(2)得：

 

对式(7)左右两边在ri≤r≤ri-1内积分，可以得到其在法定单位下的计算公式为：

 

其中：

 

将i＝1，2，…，m分别代入式(8)中，可以得到由m个方程构成考虑变渗透率模量的异常高压气藏产能计算模型。由于该产能计算模型是含有柳个未知数(未知数分别为r1，r2，…，rm-1，Qsc)的非线性方程组。因此采用Newton-Raphson方法^[21]对该产能模型进行求解，从而可以得到给定井底流压下气井的产量，进而分析变渗透率模量和高速非达西作用对异常高压气藏产能的影响。

3　实例分析

以某异常高压气藏为例，其中1口井的资料如下：j＝0.09，K0＝0.708mD，re＝l000m，rw=0.1m，h=15m，m=0.02mPa·s，Z＝1.15，T＝95.5℃，gg=0.57，p0=65MPa。

采用拟三轴应力岩心夹持器测定该井取样岩心在不同有效应力下的渗透率，发现在储层压力范围为0～65MPa时，表现为变渗透率模量特征，研究发现可以用3段渗透率模量来表征渗透率随有效应力变化的关系(表1)。

 

图2为应力敏感作用对气井IPR曲线的影响，从图中可以看到，应力敏感作用会使气井的产量降低，特别是在井底流压较低时，产量降低得越严重；对于渗透率模量计算的IPR曲线，在相同的井底流压下，常渗透率模量越大，气井的产量越低；对于该气井用变渗透率模量计算的IPR曲线与用常渗透率模量计算的IPR曲线在形态上具有相似性，但是当井底流压较小时，用变渗透率模量计算的产量与用常渗透率模量计算的产量差异较大，由于变渗透率模量更能够反映储层的渗流特征，因此用传统的常渗透率模量计算的产量可能会错误估计应力敏感作用对气井产量的影响。

 

图3为r1与r2随井底流压的变化关系图，由于实验发现在有效应力的变化范围内存在着3段渗透率模量，因此从井眼半径r3(或rw)到供给半径r0(或re)之间存在着r2与r1将气藏在径向上分为3个区域，每个区域由于有效应力的不同而呈现出不同的渗透率变化规律，从图中看到，当井底流压小于30MPa时，r1与r2随着井底流压的增大都会减小，当井底流压大于30MPa且小于56MPa时，由于地层中有效应力减小，在径向上则只存在两个区域，且分界处r1随着井底流压的增大会进一步减小，直到当井底流压大于56MPa时，则退化为了常渗透率模量的产能模型。

 

图4为考虑变渗透率模量下高速非达西对气井IPR曲线的影响，从图中可以看到高速非达西效应将会使异常高压气藏气井的产量降低，并且随着井底流压的减小，高速非达两效应使得气井产量降低得更严重。

 

4　结论

1)在异常高压气藏开发过程中，储层岩石存在着应力敏感效应，常渗透率模量在有效应力变化范围较大时难以准确表征储层的应力敏感特征，因此基于常渗透率模量的产能计算方法不能有效评价异常高压气藏的产能特征。

2)建立了考虑变渗透率模量和高速非达西效应的异常高压气藏产能计算模型并用Newton Raphson方法对模型进行求解。

3)实例分析表明：应力敏感作用会使气井的产量降低，特别是在井底流压较低时，产量降低得越严重；在井底流压较小时，用传统的常渗透率模量模型计算的产量可能会错误估计应力敏感作用对气井产量的影响；高速非达西效应将会使异常高压气藏气井的产量降低，随着井底流压的减小，高速非达西效应使气井产量降低的程度越严重。

 

符号说明

P0为原始地层压力，MPa；p为目前地层压力，MPa；K0为原始渗透率，mD；K为目前渗透率，mD；a为渗透率模量，Mpa^-1；r为平面径向距离，m；m为黏度，mPa·s；u为气体在地层中的渗流速度，m/s；b为湍流系数，m^-1；r为地层气体密度，kg/m³；Qsc为地面标准状态下的产量，m³/d；psc为标准大气压，Mpa；Z为气体偏差，无因次；T为地层温度，K；h为储层百度，m；Tsc为标准状态下温度，K；gg为地层气体的相对密度，无因次；R为通用气体常数，Mpa·m³/(kmol·K)；j为孔隙度，无因次。

 

参考文献

[1]卞小强，杜志敏，陈静.二项式物质平衡方程预测异常高压气藏储量[J].西南石油大学学报：自然科学版，2010，32(3)：75-79.

BIAN Xiaoqiang，DU Zhimin，CHEN Jing．A quadratic material balance equation of calculation of initial gas in place for abnormal geo-pressured gas reservoirs[J]．Journal of Southwest Petroleum University：Science＆Technology Edition，2010，32(3)：75-79．

[2]郑琴，陶自强，孟立新．异常高压气藏流固耦合的数值模拟[J]．天然气工业，2013，33(4)：70-74．

ZHENG Qin，TAO Ziqiang，MENG Lixin．Numerical simulation of fluid-solid coupling of abnormally pressured gas reservoirs[J]．Natural Gas Industry，2013，33(4)：70-74．

[3]傅春梅，王怒涛，邹一锋，等．异常高压气井井底压力影响因素分析[J]．西南石油大学学报，2007，29(11)：50-52．

FU Chunmei,WANG Nutao，ZOU Yifeng，et al．Influence on well-bottom pressure for abnormally high pressure well[J]．Journal of Southwest Petroleum University，2007，29(11)：50-52．

[4]刘景东，蒋有录．东濮凹陷中央隆起带北部古近系异常高压与油气成藏的关系[J]．天然气工业，2012，32(12)：30-36．

LIU Jingdong，JIANG Youlu．Relationship between abnormally high pressure and hydrocarbon accumulation of the Paleogene reservoirs in the northern part of central uplift，Dongpu Sag，Bohai Bay Basin[J]．Natural Gas Industry，2012，32(12)：30-36．

[5]李传亮．异常高压气藏开发上的错误认识[J]．西南石油大学学报，2007，29(2)：166-169．

LI Chuanliang．Misunderstanding the performance of abnormal pressure gas reservoir[J]．Journal of Southwest Petroleum University，2007，29(2)：166-l69．

[6]梁冲，尚会昌，卢拥军，等．压前测试评价在异常高压致密气藏应用研究[J]．西南石油大学学报：自然科学版，2011，33(1)：95-101．

LIANG Chong，SHANG Huichang，LU Yongjun，et al．Application and research of pre frae testing technology in abnormal high pressure tight gas reservoir[J]．Journal of Southwest Petroleum University：Science＆ Technology Edition，2011，33(1)：95-101．

[7]李乐忠，李相方．储层应力敏感实验评价方法的误差分析[J]．天然气工业，2013，33(2)：48-51．

LI Lezhong，LI Xiangfang．Error analysis of experimental evaluation methods for stress sensitivity of reservoirs[J]．Natural Gas Industry，2013，33(2)：48-51．

[8]梁斌，张烈辉，李闽，等．用数值模拟方法研究气井产量递减[J]．西南石油大学学报：自然科学版，2008，30(3)：106-109．

LIANG Bin,ZHANG Liehui，LI Min，et al．The gas well productivity decline research by means of numerical sinmlation[J]．Journal of Southwest Petroleum University：Science& Technology Edition，2008，30(3)：106-109．

[9]杨胜来，肖香娇，王小强，等．异常高压气藏岩石应力敏感性及其对产能的影响[J]．天然气工业，2005，25(5)：94-95．

YANG Shenglai，XIAO Xiangjiao，WANG Xiaoqiang，et al．Stress sensitivity of rock and its influence on productivity for gas reservoirs with abnormal high pressure[J]．Natural Gas Industry，2005，25(5)：94-95．

[10]刘启围，江晓敏，程晓旭，等．异常高压气井流入动态特征及产能试井分析方法研究[J]．油气井测试，2008，17(4)：5-7．

LIU Qiguo，JIANG Xiaomin，CHENG Xiaoxu，et al．Study on test analysis method for abnormally pressured gas well about inflow performance characteristic and productivity[J].Well Testing，2008，l7(4)：5-7．

[11]郭晶晶，张烈辉，涂中．异常高压气藏应力敏感性及其对产能的影响[J]．特种油气藏，2010，17(2)：79-81．

GUO Jingjing，ZHANG Liehui，TU Zhong．Stress sensitivity and its influence on productivity in gas reservoirs with abnormally high pressure[J]．Special Oil and Gas Reservoirs，2010，17(2)：79-81．

[12]王建光，廖新维，杨永智．超高压应力敏感性气藏产能评价方法[J1．新疆石油地质，2007，28(2)：216-218．

WANG Jianguang，LIAO Xinwei，YANG Yongzhi．A method for productivity appraisal of supper-high pressure and rock stress-sensitive gas reservoirs[J]．Xinjiang Petroleum Geology，2007，28(2)：216-218．

[13]胥洪俊，范明国，康征，等．考虑渗透率应力敏感的低渗气藏产能预测公式[J]．天然气地球科学，2008，19(1)：145-147．

XU Hongjun，FAN Mingguo，KANG Zheng，et al．A productivity prediction equation considering rock permeability stress-sensitivity in low-permeability gas reservoirs[J]．Natural Gas Geoscience，2008，l9(1)：145-147．

[14]郭春华，周文，冯文光，等．四川河坝异常高压气藏的产能方程及其计算[J]．天然气工业，2009，29(6)：86-88．

GUO Chunhua，ZHOU Wen，FENG Wenguang，et al．Building up and calculation of a deliverability equation for Heba gas reservoir with abnormal high pressure in Sichuan Basin[J]．Natural Gas Industry，2009，29(6)：86-88．

[15]ZHANG M Y，AMBASTHA A K．New insights in Dressure transient analysis for stress-sensitive reservoirs[C]//paper 28420-MS presented at the SPE 69^th Annual Technical Conference and Exhibition，25-28 September l994．New Orleans，Louisiana，USA．New York：SPE，1994．

[16]尹洪军，何应付．变形介质油藏渗流规律和压力特征分析[J]．水动力学研究与进展：A辑，2002，17(5)：538-546．

YIN Hongjun，HE Yingfu．Analysis on flow through porous media and pressure characteristics deformed reservoirs[J]．Journal of Hydrodynamics，Ser．A，2002，17(5)：538-546．

[17]王建忠，姚军，张凯，等．变渗透率模量与双重孔隙介质的压力敏感性[J]．中国石油大学学报：自然科学版，2010，34(3)：80-83．

WANG Jianzhong，YAO Jun，ZHANG Kai，et al．Variable permeability modulus and pressure sellsitivitv of dual porosity medium[J]．Journal of China University of Petroleum：Natural Science Edition，2010，34(3)：80-83．

[18]王建忠，姚军，王子胜．用变渗透率模量分析双重介质油藏压力动态特征[J]．新疆石油地质，2010，31(3)：288-290．

WANG Jianzhong，YAO Jun，WANG Zisheng．Pressure behavior characteristic  analysis  of stress-sensitive reservoirs with variable permeability module[J]．Xinjiang Petroleum Geology，2010，31(3)：288-290．

[19]VAIROGS J，HEARN C L，DAREING D W，et al．Effect of rock stress on gas production from low-perinea bility reservoirs[J]．Journal of Petroleum Technology，1971，23(9)：1161-1167．

[20]VAIROGS J，RHOADES V W．Pressure transient tests in formation having stress sensitive permeability[J]．Journal of Petroleum Technology，l973，25(8)：965-970．

[21]张光澄，张雷．实用数值分析[M]．2版．成都：四川大学出版社，2009：92-97．

ZHANG Guangcheng，ZHANG Lei．Practical numerical analysis[M]．2^nd edition．Chengdu：Sichuan Univcrsily Press，2009：92-97．

[22]蒋海军，鄢捷年，李荣．裂缝性储层应力敏感性实验研究[J]．石油钻探技术，2000，28(6)：32-33．

JIANG Haijun，YAN Jienian，ＬI Rong．Experimental study on stress-sensitivity of fracturing formations[J]．Petroleum Drilling Techniques，2000，28(6)：32-33．

[23]蒋艳芳，张烈辉，刘启国，等．应力敏感影响下低渗透气藏水平井产能分析[J]．天然气工业，20ll，31(10)：54-56．

JIANG Yanfang，ZHANG Liehui，LIU Qiguo，et al.productivity analysis of horizontal wells in low Dermeability gas reservoirs considering stress sensitivity[J]．Natural Gas Industry，2011，31(10)：54-56．

[24]FORCH HEIMER P H．Wasserbewegung durch Boden[J]．Zeitschrift des Vereines Deutscher lngenicure，1901．45(5)：l781-l788．

[25]LI D，SVEC R K，ENGLER T W，et al．Modeling and simulation of the water non-Darcy flow experiments[C]//paper 68822-MS presented at SPE Western Regional Meeting，26-30 March 2001，Bakersfield，California．USA．New York：SPE，2001．

 

本文作者：任俊杰  郭平  王绍平  毕波  汪周华  赵增春

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

  中国石油长庆油田公司第二采油厂

  中国石油新疆油田公司采油一厂

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