Drying Process in Underground Gas Storage Reservoir
-
摘要: 地下储气库在注采循环过程中,井周围区域会发生储层干化,这会对储气库的正常运行产生影响。为此,在分析储层干化机理的基础上,基于大张坨储气库储层参数及地下温度压力条件,建立了单井控制区域的组分模型并模拟储气库周期性注采气过程。模拟结果表明:井周围储层干化的区域会逐年扩大,半径甚至超过百米;储层干化主要发生在注气周期,回采周期对储层的干化作用轻微;储层干化的速率及波及区域的大小取决于注采井的年循环气量及注入和采出天然气中水蒸气含量的差异;不同回采周期采出的天然气中水蒸气的含量呈现规律性变化,从开始接近干气逐渐升高到地下温度压力条件下饱和水蒸气含量,且该过程逐年后延。该研究结果为深入研究储层干化现象对储气库正常运营的影响奠定了基础。Abstract: During the injection/production cycling in underground gas storage,reservoir rock in area around the wellbore will be dried,which might affect its normal operation.Therefor after analyzing the drying mechanism,a two phase two component model for a single well controlled area was established for simulating the gas cycling process with parameters(reservoir properties and underground temperature and pressure)referenced from Dazhangtuo Gas Storage Reservoir.Simulation results show that the drying area around the wellbore would be enlarged year by year,the radius of the dried area may even exceed 100 meters;water evaporation that leads to reservoir rock drying mainly occurs in the period of dry gas injection,the drying effect during the production time could be negligible;the rate of reservoir drying and the size of affected area depends on the annual cycling gas volume for the injection-production well,as well as the water vapor content difference between the injected gas and the produced gas.The water vapor content in the produced gas during various withdrawn periods changes regularly.At the initial stage,it is near to value of the injected dry gas;then it gradually increases to the saturated water vapor content under the downhole temperature and pressure conditions;and this process tends to be postponed year by year.The research result would provide a foundation for studying the effects of reservoir drying on the normal operation of gas storage.
-
Keywords:
- gas storage /
- water vaporization /
- water vapor content /
- water saturation /
- numerical simulation /
- component model
-
塔里木油田2022年完钻超深井(6 000 m以上)占比近90%,固井质量是超深井建井的关键环节,也是保证油气井生产寿命的关键所在,对油气田勘探开发效益和产能建设非常重要。超深井由于裸眼段长,固井段温度、压力差异大,地层流体更为复杂,微间隙普遍存在会造成声幅值变高对固井质量评价产生误导,进而可能影响试油作业方案制定。塔里木油田由于存在多套易漏地层,广泛使用低密度水泥浆,低密度水泥浆会使CBL值升高[1–7],若同时存在微间隙,无法判断是由于微间隙造成的声幅升高还是存在水泥浆窜槽。针对上述难题,塔里木油田对“三高”气井生产套管段实行CBL/VDL和IBC联测的测井方案,以实现固井质量的精确评价。常规声波固井质量评价方法(如水泥胶结测井(CBL)和声波变密度测井(VDL))采用套管波幅度评价水泥的胶结状况,其物理原理是水泥和钻井液的声阻抗差异大,套后胶结水泥时套管波幅度小,若是自由套管,则声幅值明显增大,但声幅测井只能反映井周平均情况;套后胶结低密度水泥时,由于其声阻抗小于常规水泥,有时与井内钻井液声阻抗相当,这会使利用声幅变密度测井评价低密度水泥时的动态响应范围明显变小,灵敏度降低。另外,套管与水泥环之间存在微间隙,但未形成有效连接时,应认为胶结质量好,常规声幅测井易产生误判[8–10]。斯伦贝谢的套后成像测井仪(IBC)是通过斜入射在套管中激发挠曲波与垂直入射激发共振波2种方式联合测井,可较好地识别套管外胶结介质的相态,特别是挠曲波明显提高了套后胶结轻质水泥时固井质量评价的灵敏度[11]。基于此,笔者通过研究CBL套管波、挠曲波及共振波等3种模式波在不同密度水泥、不同微间隙厚度和类型条件下的响应差别,建立了超深井微间隙识别和评价方法,大幅提高了固井质量解释符合率。
1. VDL和IBC测井原理
1.1 VDL测井原理简述
常用的CBL/VDL测井仪采用单发双收声系采集水泥胶结信息,源距分别为914.4和1 524.0 mm,发射频率为20 kHz。其中源距914.4 mm接收换能器采用固定门采集套管波E1波峰幅度(见图1),固定门为在套管尺寸已知后确定的时间窗,取记录波形中此时间窗内的声波峰值,此峰值经过刻度后就是声幅值;源距1 524.0 mm接收换能器采集每一深度点的声波全波列信息,以变密度方式(VDL)记录,用于评价水泥胶结质量。CBL能反映套管−水泥环之间第一界面的胶结情况, 一般认为CBL值越大,第一界面胶结越差;CBL值越小,第一界面胶结越好。VDL不仅能反映第一界面的胶结情况, 而且能反映水泥环−地层之间第二界面的胶结情况(见图2,图中的F和J分别表示发射换能器和接收换能器) 。
1.2 IBC测井原理
IBC套后成像仪(Isolation Scanner)是斯伦贝谢公司用于固井质量及套损评价的仪器[11],它通过结合2种声波技术:经典的脉冲回波技术和最新的挠弯曲波(Lamb波)成像,适用于评价各种类型水泥,其测量原理如图3所示。脉冲回波技术通过井壁反射波幅度对套管内壁的粗糙度进行成像,井壁反射波到时可反演仪器偏心,套管共振波的频率和衰减可分别用于计算套管厚度和反演套后水泥的声阻抗(见图3(a))。脉冲回波测量时换能器辐射的声波垂直入射到套管内壁,套后成像测井仪器的另外一种测量方式是探头辐射的声波斜入射到套管内壁,在套管中激发挠曲波(Lamb波),沿着套管传播的挠曲波(Lamb波)还会向套管两侧辐射声波,利用2个接收换能器记录泄漏的声波信号计算挠曲波(Lamb波)的衰减(见图3(b))。超声回波和挠曲波(Lamb波)的测量声系在仪器中是背靠背放置的。仪器在井下工作时测量声系通过旋转扫描实现对整个环向的套后介质成像,一般旋转一圈采集36个方位的波形数据。脉冲回波测量反演的声阻抗和挠曲波衰减交会识别套后介质相态(见图3(c)),从而实现套后介质气液固成像。
2. CBL/VDL和IBC对微间隙的响应特征
2.1 CBL/VDL数值模拟
采用经典的弹性波传播理论模拟套管与水泥之间不胶结时,是在套管与水泥之间加入一层流体环(见图4)。
流体环的存在阻断了套管与水泥之间的轴向振动(相对于界面为切向振动),即使流体环为厚度0.01 mm的微间隙,也会使套管波幅度明显增强。套管波相对幅度(各胶结模型下首波的第一个正峰与自由套管模型下首波的第一个正峰的比值)随流体环厚度变化的趋势如图5所示,可以看出,微间隙厚度为0.01 mm时,套管波的幅度已接近自由套管状态下(套管与地层间充填水)套管波幅度的50%。
套管后胶结常规水泥和轻质水泥时,套管波的衰减曲线如图6所示。由图6可知,套管后有微间隙存在时,套管波的衰减接近自由套管状态的衰减值,说明CBL/VDL测井对微间隙非常敏感,微间隙的存在会使套管波幅度明显增大,衰减降低,套管后耦合任意类型水泥时的响应特征都是一致的[8–10]。当井下实际存在厚度较小微间隙时,套管与水泥环之间会存在固体颗粒,有时仍可观测到地层波。如图7所示,L1井在图示井段采用密度1.88 kg/L常规水泥固井,第一次测井时钻井液密度为1.52 kg/L,第二次测井时钻井液密度为1.80 kg/L,2次测井采用同一只测井仪器,2次测井时间间隔24 h。该井X 000~X 045 m段第一次测井平均声幅值(第3道黑色曲线)在40以上,地层波也较为明显;但第二次测井平均声幅值(第3道红色曲线)出现明显下降,声幅平均值在10以下,地层波幅度也明显增大,与裸眼井声波时差曲线对应关系更好。这说明该段固井质量第一界面存在较小的微间隙,且存在固体颗粒,所以第一次测井时虽然声幅值较高,但地层波依然可见。
2.2 IBC中的共振波和挠曲波的响应特征
S. Zeroug等人[11]提出了通过控制声束入射角的方式来有效激发不同模式波的方法。图8所示为理论计算的不同角度入射时在浸水钢板可激发的模式波的频散曲线,可见在垂直入射时,套管中主要激发的是高阶对称模式S1、S2及S3等;入射角大于30°时,激发的是套管中的挠曲波[12]。
CBL/VDL测井对微间隙敏感,但是垂直入射下的套管共振回波对微间隙不敏感。图9所示为快速水泥无地层模型下计算的共振波幅度随套管与水泥间流体层厚度的变化趋势,从右侧的放大图可见,微间隙厚0.01~0.10 mm时,对共振波基本无影响,这与斜入射下激发的挠曲波的响应规律不同(见图10),对于快速水泥(高阻抗水泥),一旦出现厚度0.01 mm的微间隙,挠曲波的幅度明显下降;随着微间隙厚度增大,其下降幅度还会逐渐增大;从水泥与地层界面的反射波幅度与直达波的幅度变化趋势正好相反,说明沿着套管传播的挠曲波泄露的能量越大,衰减越快,幅度越低,则从水泥环外界面的反射波幅度就越大。
套管后胶结慢速水泥(轻质水泥或低密度水泥)时计算的共振波幅度随套管与水泥间流体层厚度的变化如图11所示。由图11(b)可知,厚度为0.01~0.10 mm的微间隙对共振波基本无影响,流体环厚度继续增大,共振波的幅度明显增大。
套管后耦合轻质水泥时,若微间隙厚度小于0.30 mm,其对挠曲波幅度变化的影响很小,但从水泥−地层界面反射的SS波幅度减小明显,P/S或S/P波幅度变化不大;随着微间隙厚度增加,挠曲波的幅度增加明显,衰减明显减小(见图12)。
水泥胶结良好和存在微间隙时,挠曲波衰减率随套后水泥阻抗的变化规律如图13所示。从图13可以看出,胶结良好时,随着水泥声阻抗增大,衰减逐渐增大,阻抗值达到4.7 Mrayls时衰减取得最大值;随着水泥声阻抗继续增大,衰减急剧下降,这与较大声阻抗水泥的纵波速度大于套管挠曲波相速度时只向水泥环辐射横波有关。挠曲波的衰减与水泥速度之间的复杂关系,即相同的衰减率对应套管外不同的声阻抗耦合介质,再次说明不能像SBT测井方法仅依赖挠曲波的衰减来评价水泥环第一界面的胶结质量,应建立声衰减−声阻抗图版来评价固井质量。
套管与水泥环之间微间隙厚度小于0.10 mm、水泥阻抗小于4.7 Mrayls时,衰减稍有降低;水泥阻抗大于4.7 Mrayls时,微间隙的出现降低了套后介质的等效声阻抗,使挠曲波的衰减急剧增大。因此,联合CBL/VDL和IBC测井数据可确定井下是否有微间隙存在,CBL高值或VDL变密度图中套管波幅度强,但IBC测量的衰减值相对较高是固体的响应特征,则可判断井下有微间隙存在[13–16]。
3. 实例分析
H1井是塔里木油田H区块一口开发井,图14中井段采用密度1.88 kg/L的常规水泥固井。根据裸眼井测井资料,X 076.50~X 078.50 m井段解释为油层,X 088.50~X 097.00 m井段解释为水层。为精细评价固井质量,目的层固井后同时采集了CBL/VDL和IBC资料。X 078.50~X 088.50 m井段第9道的声幅值为30~45,表明第一界面胶结不好,VDL显示地层波也不明显,综合评价为差,所以射孔后可能导致窜槽。但从IBC测井资料来看,声阻抗值高,挠曲波衰减大,SLG成像图显示全部为固体填充,所以X 078.50~X 088.50 m井段可能存在液态微间隙,导致声幅值高,对油水层来说并不影响水泥环封隔能力。后期对X 076.60~X 078.40 m井段射孔,产油量53.0 m3/d,不含水,未发生窜槽,验证了微间隙识别方法的准确性。
塔里木油田M区块M1井某井段采用密度为1.38 kg/L的低密度水泥固井,为精细评价固井质量,固井后采集了CBL/VDL和IBC资料(见图15)。第9道的声幅值高达50~60,指示第一界面胶结不好,且VDL上接箍明显,地层波不可见,综合评价为胶结差;但从IBC测井资料看,波阻抗平均在4 Mrayls左右,挠曲波衰减平均在120 dB/m左右,指示第一界面胶结良好,套管附近为固体填充。综合CBL/VDL和IBC资料,判断X 545~X 580 m井段存在液体填充的微间隙,不影响水泥环封隔性。X 550~X 575 m井段采用封隔器进行射孔验窜作业,不窜,证明水泥胶结好,CBL高值是由于液体填充微间隙造成的。
4. 结 论
1)CBL/VDL测井时激发的套管波沿套管传播时与套管后物质主要以剪切刚度耦合,因此即便是厚度0.01 mm的微间隙,也会使套管波衰减明显降低,接近套管后为流体时的响应特征。因此,CBL/VDL测井时若声幅值出现大段异常高值,可能存在窜槽或者微间隙,可结合IBC测井进一步验证。
2)IBC测井在套管中产生的挠曲波与套管后物质主要以法向刚度耦合,使套管后胶结轻质水泥时的挠曲波对微间隙不敏感,套管后耦合常规水泥时的微间隙会导致挠曲波的衰减明显增强。因此,CBL/VDL和IBC测井二者结合,当CBL套管波幅度较高时,若挠曲波的衰减较高,则可判定为套管与水泥环之间存在微间隙。
-
[1] KATZ D L,RASIN TEK M.Overview onunderground storage of natural gas[J].Journal of Petroleum Technology,1981,33(6):943-951.
[2] 丁国生,李文阳.国内外地下储气库的现状与发展趋势[J].国际石油经济,2002,10(6):23-26. DING Guosheng,LI Wenyang.Domestic and foreign underground gas storage tanks:current situation and development trends[J].International Petroleum Economics,2002,10(6):23-26. [3] CORNOT-GANDOLPHE S.Underground gas storage in the world-2017 status[R].Rueil-Malmaison:CEDIGAZ,2016.
[4] 丁国生,李春,王皆明,等.中国地下储气库现状及技术发展方向[J].天然气工业,2015,35(11):107-112. DING Guosheng,LI Chun,WANG Jieming,et al.The status quo and technical development direction of underground gas storage in China[J].Natural Gas Industry,2015,35(11):107-112. [5] 魏欢,田静,李建中,等.中国天然气地下储气库现状及发展趋势[J].国际石油经济,2015,23(6):57-62. WEI Huan,TIAN Jing,LI Jianzhong,et al.The status quo and development trend of underground gas storage in China[J].International Petroleum Economics,2015,23(6):57-62. [6] 陆争光.中国地下储气库主要进展、存在问题及对策建议[J].中外能源,2016,21(6):15-19. LU Zhengguang.Major progress in China’s underground gas storage construction,obstacles and countermeasures[J].Sino-Global Energy,2016,21(6):15-19. [7] HAWKINS A B,McCONNELL B J.Sensitivity of sandstone strength and deformability to changes in moisture content[J].Quarterly Journal of Engineering Geology Hydrogeology,1992,25(2):115-130.
[8] ZHOU Zilong,CAI Xin,CAO Wenzhuo,et al.Influence of water content on mechanical properties of rock in both saturation and drying processes[J].Rock Mechanics Rock Engineering,2016,49(8):3009-3025.
[9] McKETTA J J,WEHE A H.Use this chart for water content of natural gases[J].Petroleum Refiner,1958,37(8):153-154.
[10] TERRIGEOL A,COQELET C,CHAPOY A.Water content assessment in acid gas[C].Annual Conference of Gas Processors Association-Europe,Florence,Italy,September 16-18,2015.
[11] 崔立宏,疏壮志,杨树合,等.大张坨地下储气库建设方案[J].西南石油学院学报,2003,25(2):76-79. CUI Lihong,SHU Zhuangzhi,YANG Shuhe,et al.Construction scheme of Dazhangtuo underground gas storage[J].Journal of Southwest Petroleum Institute,2003,25(2):76-79. [12] 何顺利,门成全,周家胜,等.大张坨储气库储层注采渗流特征研究[J].天然气工业,2006,26(5):90-92. HE Shunli,MEN Chengquan,ZHOU Jiasheng,et al.Study on percolation characteristics of reservoirs’ injection-production in Dazhangtuo underground gas storage[J].Natural Gas Industry,2006,26(5):90-92. [13] BAHADORI A.Natural gas processing:technology and engineering design[M].Massachusetts:Gulf Professional Publishing of Elsevier,2014.
[14] FOSS M M.Interstate natural gas:quality specifications interchangeability[R].Texas:Center for Energy Economics,2004.
[15] AS 4564-2005.Specification for general purpose natural gas[S].
[16] DVGW G 260(A). Gasbeschaffenheit.[EB/OL].[2017-05-24].https://www.dvgw-regelwerk.de/plus/#technische-regel/dvgw-g-260/cec1c3.
[17] EIE/06/221/SI2.442663.Intelligent energy-Europe[R].Hungary:University of Miskolc and Hungarian Scientific Society of Energy Economics,2008.
[18] GB 17820-2012.天然气[S]. GB 17820-2012.Natural gas[S]. [19] 肖香姣,姜汉桥,刘曰武,等.高压气井井筒附近地层温度分布计算方法[J].科技导报,2011,29(32):20-24. XIAO Xiangjiao,JIANG Hanqiao,LIU Yuewu,et al.Computing method for the formation temperature distribution of high pressure gas well near wellbore layer[J].Science Technology Review,2011,29(32):20-24. [20] ZULUAGA E,MUNOZ N I,OBANDO G A.An Experimental study to evaluate water vaporization and formation damage caused by dry gas flow through porous media[R].SPE68335,2001.
[21] ZULUAGA E,MONSALVE J C,CASTRILLON J A,et al.Experiments on water vaporization in porous media[J].Journal of Canadian Petroleum Technology,2003,42(7):7-8.
[22] ZULUAGA E,MONSALVE J C.Water vaporization in gas reservoirs[R].SPE84829,2003.
[23] MAHADEVAN J.Flow-through drying of porous media[D].Austin:The University of Texas,2005.
[24] MAHADEVAN J,SHARMA M M,YORTSOS Y C.Flow-through drying of porous media[J].AIChE Journal,2006,52(7):2367-2380.
[25] MAHADEVAN J,SHARMA M M,YORTSOS Y C.Water removal from porous media by gas injection:experiments and simulation[J].Transport in Porous Media,2007,66(3):287-309.
[26] SΦREIDE I,WHITSON C H.Peng-Robinson predictions for hydrocarbons,CO2,N2,and H2S with pure water and NaCl brine[J].Fluid Phase Equilibria,1992,77:217-240.
[28] DUAN Zhenhao,MAO Shide.A thermodynamic model for calculating methane solubility,density and gas phase composition of methane-bearing aqueous fluids from 273 to 523 K and from 1 to 2000 bar[J].Geochimica et Cosmochimica Acta,2006,70(13):3369-3386.
[29] GOLGHANDASHTI H,SAADAT M,ABBASI S,et al.Experimental investigation of salt precipitation during gas injection into a depleted gas reservoir[R].IPTC15179,2011.
[30] GOLGHANDASHTI H,SAADAT M,ABBASI S,et al.Experimental investigation of water vaporization and its induced formation damage associated with underground gas storage[J].Journal of Porous Media,2013,16(2):89-96.
-
期刊类型引用(6)
1. 邓才,王振华,王立垒. 近圆形高含硫气藏辐射状分支水平井产量公式. 石油地质与工程. 2022(06): 72-76 . 百度学术
2. 王苏雯,王寿喜,王力,王勇. 火驱尾气回注井筒流动分布规律. 断块油气田. 2020(01): 99-103 . 百度学术
3. 刘人玮,杨天宇,万宇飞,蒋曙鸿,黄喆. 海上平台热采放喷气扩散数值模拟研究. 石油化工高等学校学报. 2020(03): 74-79 . 百度学术
4. 王永佩. 水平井过热蒸汽循环预热井筒传热模拟. 北京石油化工学院学报. 2019(02): 11-16 . 百度学术
5. 聂帅帅,唐世星,刘可,徐康泰,李江飞,王少征. 数据挖掘诊断X油田低渗透稠油油藏压裂效果的主控因素. 石油地质与工程. 2019(06): 90-94 . 百度学术
6. 王永佩. 超临界水井筒非等温管流数值模拟. 北京石油化工学院学报. 2019(04): 6-8 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 3968
- HTML全文浏览量: 148
- PDF下载量: 5886
- 被引次数: 7