Analysis on the elastoplastic of cement sheath in underground gas storage wellbore
-
摘要: 为了分析连续注采工况下地下储气库井筒水泥环的受力与变形规律,基于弹塑性力学的基本原理,以套管-水泥环-围岩组合体为研究对象,结合温度、孔隙压力及地层远场应力等影响因素,推导出注采工况下受内外压作用的水泥环的弹塑性应力和位移的解析解。实例计算表明:注气中水泥环会发生塑性变形,内外壁面径向接触应力随内压非线性增加;注采中的井筒温度和孔隙压力变化是影响水泥环受力的重要因素;水泥环弹性极限内压随温度的降低和孔隙压力的增大而增大;水泥环与套管、围岩的径向接触应力随温度的升高和孔隙压力的增大而增大。研究成果为优化储气库运行参数、控制“微环空”的产生提供了技术支持。Abstract: In this paper, the casing-cement sheath-surrounding rock assembly was taken as the study object to reveal the stress and deformation of cement sheath in underground gas storage wellbore in the case of continuous injection and production. The analytical solutions to the elastoplastic stress and displacement of cement sheath subjected to internal and external pressure in the working conditions of injection and production were derived based on the basic principles of elastoplastic theory combined with the influential factors, e.g. temperature, pore pressure and far-field in-situ stress. The calculation examples show that in the process of gas injection, the cement sheath may deform plastically and the radial contact stress on the inner and outer wall increases nonlinearly with the increase of internal pressure. The variation of temperature and pore pressure during gas injection and production has significant effect on the stress of cement sheath. The elastic limit internal pressure of cement sheath increases with the increasing of pore pressure and the decreasing of temperature. As temperature and pore pressure increase, the radial contact stress of cement sheath on the casing and surrounding rock increases. The research results provide the technological support for optimizing underground gas storage operating parameters and preventing the formation of micro-annulus.
-
Keywords:
- underground gas storage /
- cement sheath /
- elastoplastic theory /
- analytical solution /
- pore pressure
-
地下储气库井筒水泥环服役条件复杂, 在天然气的循环注入和采出作业过程中, 井筒温度、气体内压、运行压力及地应力场会交替变化。在多期注采后, 受交变载荷作用水泥环的完整性可能失效[1-3], 导致气体泄漏, 引起环空持续带压。其失效形式包括水泥与地层、水泥与套管胶结面密封失效形成的“微环空”、水泥环发生径向开裂和塑性变形等。此外, 作业载荷、热应力、钻井液清除不彻底、地层表面粗糙、水泥的失水收缩均会引起水泥环胶结界面密封失效[4-5]。
关于井筒水泥环失效机理与预测模型的研究, 众多学者基于弹性和弹塑性力学原理将其转化为平面应变问题, 建立解析模型并求解。Thiercelin等[6-10]基于线弹性假设, 假设地层、水泥及套管均为各项同性均质线弹性材料, 建立了水泥环受力的解析模型, 对水泥环与井筒的界面压力及影响因素进行分析。李军等[11-12]考虑水泥的塑性变形, 基于弹塑性力学基本原理分析了套管-水泥环-地层组合的变形与失效。相关试验研究为验证与完善理论分析提供了依据。唐汝众等[13]通过水泥环的高内压物理模拟试验, 测定了双向应力作用下的应变值和破坏载荷。有限元分析是研究水泥环受力与完整性失效的有效方法。Andrade等[14-15]建立了水泥环受力的有限元模型, 分析了材料力学性能、几何形状及井筒载荷特别是热应力等因素对水泥环受力与失效的影响, 并与解析结果进行对比。此外, 水泥环的受力与失效会受原地应力、温度及压力变化的影响, 需要考虑非均匀地应力场、温度场的影响[16-18]。上述研究均未结合地下岩石的孔隙弹性特征分析水泥的弹塑性受力, 但是在实际地下储气库注采中孔隙压力、温度发生变化, 会导致原地应力场变化[19]。
基于弹塑性力学理论, 考虑注采中的孔隙压力、温度及远场地应力等因素, 分别采用热弹性力学、弹塑性力学及孔隙弹性理论分析套管、水泥及地层的受力。结合储气库实际工况条件分析了均匀地应力场中水泥环的弹塑性受力与变形, 对于优化地下储气库运行参数, 保持井筒完整性具有一定指导意义。
1. 模型建立
1.1 基本假设
套管内压和温度变化会使套管径向膨胀, 地层围岩在远场应力作用下对水泥环产生挤压并限制其径向膨胀, 在水泥环第一、第二胶结界面产生径向接触应力, 可能导致水泥环剪切或拉伸失效(图 1, a、b、c、d分别为套管内半径、套管外半径、水泥环外半径、井周围岩外半径, pc1为水泥环与套管胶结面接触压力, pc2为水泥环与围岩胶结面的径向接触压力, pf为近井远场径向应力, pi为套管内压)。基于热-孔隙弹性和弹塑性理论, 假设: ①套管为完全弹性体, 水泥环为理想弹塑性体, 地层为孔隙弹性介质; ②固井质量优良, 水泥环与套管、地层完全胶结, 无窜槽; ③从近井围岩到井筒沿径向的温度没有变化, 当工况变化以后, 同一深度各组件的温度变化值相同, 即∂(ΔT)/∂r=0;④地层为各向同性介质, 水平地应力为均匀应力场, 随孔隙压力同步均匀变化; ⑤研究对象为直井, 套管-水泥环-地层组合体的变形属于平面应变问题, εz=0。
1.2 套管和地层受力
对于套管-水泥环-地层组合体, 套管为各向同性完全弹性材料, 考虑温度引起的热膨胀效应, 应力-应变本构方程为:
(1) 式中: εr、εθ、εz分别为径向、周向、轴向应变; Es为套管的杨氏模量, MPa; σr、σθ、σz分别为径向、周向、轴向应力, MPa; μs为套管的泊松比; αs为套管的热膨胀系数, ℃ -1; ΔT为温度变化值, ℃。
对于平面应变问题, 轴向应力σz可以简化为:
(2) 代入式(1), 得到套管周向应变方程:
(3) 对于孔隙弹性材料(地层围岩), 考虑孔隙压力的影响, 应力-应变本构关系为:
(4) 式中: Ef为地层围岩的杨氏模量, MPa; μf为地层围岩的泊松比; β为近井围岩的Biot系数; pp为地层孔隙压力, MPa; αf为地层围岩的热膨胀系数, ℃-1。
同理, 根据平面应变原理对轴向应力进行简化, 回代至式(4), 得到孔隙弹性材料的周向应变方程:
(5) 式(3)和式(5)分别为套管(弹性材料)和地层岩石(孔隙弹性材料)的本构方程。根据厚壁圆筒弹性理论, 套管和地层的径向、周向应力可以通过拉梅公式计算, 设定应力应变符号以拉伸为正。对于套管
(a < r < b), 其径向和周向应力分布为:
(6) 式中: r为套管横截面上某一点的半径, mm。几何方程为:
(7) 将式(6)代入式(3), 并结合式(7)可以得到套管外边界(r=b)的位移ur-casing:
(8) 同理, 对于井周围岩(c < r < d), 根据拉梅公式能够确定周向和径向应力, 并结合孔隙弹性介质本构方程式(5)和几何方程式(7), 得到地层围岩内壁面(r=c)的位移uf:
(9) 在实际应用中, 考虑到地层应力扰动区的尺度远大于井眼本身, 即c << d, 则式(9)可以简化为:
(10) 1.3 水泥环弹塑性变形
将水泥环视为理想弹塑性体, 采用岩石力学中广泛应用的摩尔-库伦屈服准则, 屈服条件可以表述为:
(11) 式中: C为水泥石的内聚力, MPa; φ为水泥石的内摩擦角, (°)。
随内压增加, 水泥环将首先从内侧开始屈服并逐渐向外扩展。设弹塑性分界面半径为rp, 水泥环壁厚为t, 塑性区厚度为tp。根据理想弹塑性厚壁圆筒理论分析水泥环弹塑性应力与变形(图 2)。
将屈服条件与平衡方程联立并结合边界条件(r=b, σr=-pc1), 可得塑性区(b < r < rp, rp=b+tp)的应力分布[11, 20]:
(12) 其中, 当pc1 < pc2 时,
; pc1≤pc2 时, 。当r=rp时, 水泥环弹塑性分界面的应力pm满足σr =-pm, 可得:
(13) 弹性区应力分布可以通过拉梅公式计算得到。
当r=rp时, 水泥环弹塑性分界面上的应力同时满足弹性区应力分布和屈服准则式(11), 可得:
(14) 其中, 当pc1 > pc2 时, D=1;当pc1≤p c2 时, D=-1。
同前述套管位移推导过程, 水泥环弹性区的径向位移ur-e 可由弹性厚壁筒平面应变问题理论求解:
(15) 式中: E c为水泥环的杨氏弹性模量, MPa; μc为水泥环的泊松比; αc为水泥环的热膨胀系数, ℃-1。
根据塑性理论, 塑性体应变为0, 结合非关联流动法则, 由体积弹性定律[21]可得水泥环塑性区径向位移同前述套管位移推导过程, 水泥环弹性区的径向ur-p的微分方程:
(16) 边界条件: 当r=rp时, ur-e(rp)=ur-p(rp)。
将式(12)代入以上微分方程, 结合边界条件和式(13)进行积分并确定积分常数, 可得塑性区的位移ur-p :
(17) 据此可计算得到水泥环塑性区内边界(r=b)的径向位移ur-cement1:
(18) 根据式(15), 水泥环外边界(r=c)的径向位移ur-cement2为:
(19) 综上, 套管-水泥环-地层紧密相连, 不同介质交界的内外壁面位移均需参考式(15)进行求解, 然后根据面处的位移连续, 由式(8)、式(10)、式(18)、式(19)位移连续条件即可求得水泥环的受力与变形。可得水泥环-套管和水泥环-地层胶结面的位移连续条件, 并结合式(13)、式(14)得到套管-水泥环-井周围岩组合的弹塑性变形方程组:
(20) 求解该方程组可得水泥环的受力和位移分布情况。水泥环内壁刚开始发生塑性变形时, r
p=b, 求解方程组可得套管临界内压值, 称为弹性极限内压pie; 当rp=c时, pm=pc2, 水泥环全部发生塑性变形, 相应的临界内压称为塑性极限内压pis, 可通过逐渐增大pi, 由该方程组试算rp直至满足|rp-c < δ|求得, 其中δ为精度要求, 取10-9。 当内压小于pie 时, 水泥只发生弹性变形, 水泥环的内外壁面位移均需参考式(15)进行求解, 然后根据位移连续条件即可求得水泥环的受力与变形。
2. 实例计算
2.1 具体实例
某地下储气库注采井深2 500 m, 对应深度处远场地层温度为90 ℃。套管外径为177.8 mm, 内径为159.4 mm; Es=207 GPa, 泊松比为0.3, αs=1.21×10-5℃-1 。水泥环采用G级水泥, 内摩擦角为30°, 内聚力为6.0 MPa; 外径为222.2 mm; Ec=25 GPa, μc=0.25, αc=0.7×10-5 ℃-1。井周围岩为石灰岩, Ef=40 GPa, μf=0.35, β=0.9;αf=0.5×10-5 ℃-1。相国寺石炭系气藏地下储气库运行压力(即注采中孔隙压力pp)范围为11.7~28.0 MPa。注采中, 近井底位置的套管内压最大, 水泥环最易失效, 因此, 对井底位置的水泥环进行受力分析。
注气中井底套管-水泥环-围岩组合体温度变化∆T=-20 ℃, pp=11.7 MPa, 对于非蠕变地层, 围岩远场应力对井周的挤压较弱, pf=10 MPa。计算得到储气库水泥环弹性极限内压pie=33.42 MPa, 水泥环塑性极限内压pis=47.65 MPa(图 3)。在该工况下, 应该控制最高套管内压不超过33.42 MPa, 避免水泥环发生塑性变形, 在内压卸载中可能造成胶结面分离, 形成“微环空”。
分别采用弹性和弹塑性理论计算水泥环内外壁面径向接触应力随内压的变化(图 4)。水泥环第一、第二胶结面的径向接触应力随套管内压的增加而增加, 且套管-水泥接触力大于水泥-地层接触力。根据弹塑性解, 水泥环径向接触应力随内压非线性增加。当内压超过弹性极限内压时, 接触压力的弹塑性解大于弹性解, 显然在地下储气库的运行设计中采用弹塑性理论更安全, 且弹塑性解可以得到水泥环的弹性极限内压, 为地下储气库井筒控制水泥环塑性变形提供了依据。
取套管内压pi=45.0 MPa, 由水泥环的应力和位移沿径向分布(图 5、图 6)可知, 水泥环应力和位移呈非线性分布, 其中径向应力沿半径减小, 周向应力沿半径先减小后增大, 在弹塑性分界面上发生突变, 该工况下, 二者均为压应力, 径向位移沿半径逐渐减小。
2.2 温度影响
保持其他参数不变, 计算不同温度变化下水泥环的弹性极限套管内压和内外壁面径向接触应力随内压变化(图 7、图 8)。随温度的升高, 水泥环弹性极限内压降低, 水泥环更容易发生塑性变形。套管内压相同时, 随温度升高水泥环内外壁面径向接触应力增加, 这是由于热膨胀在水泥胶结面上产生了较大的热应力。因此在注气过程中, 水泥环弹性极限内压将升高, 注入低温度气体有利于提高许用注气压力。
2.3 孔隙压力影响
保持温度变化ΔT=-20℃, 计算不同孔隙压力下水泥环塑性区扩展和内外壁面径向接触应力随内压变化规律(图 9、图 10)。随地层孔隙压力的增加, 水泥环弹性极限套管内压增加, 水泥环承压能力提高, 水泥环与套管和围岩的径向接触应力增加。这是由于孔隙压力增加, 近井围岩承受的有效应力减小, 地层位移减小, 对水泥环的限制增强。因此在储气库运行过程中, 随着气体的注入, 水泥环内外壁面径向接触应力增加。孔隙压力对水泥环的受力影响较大, 在储气库水泥环完整性分析中应该考虑孔隙压力变化。
3. 结论
(1) 考虑地下储气库运行中的温度、孔隙压力变化, 基于弹塑性力学原理, 推导了具有Mohr-Coulomb屈服条件的受内外压水泥环的弹塑性应力和位移的解析解。
(2) 由弹塑性理论得到水泥环的应力、位移呈非线性分布, 水泥环内外壁面径向接触应力随内压非线性增加。由弹塑性解可得水泥环的弹性极限内压, 有助于控制地下储气库运行中的“微环空”。
(3) 注采作业引起套管-水泥环-地层的温度、孔隙压力变化是影响水泥环受力和失效的重要因素。水泥环与套管和围岩的径向接触应力随温度的升高和孔隙压力的增大而增大。
(4) 远场应力对水泥环受力和变形的影响较大, 在此仅考虑了均匀地应力场的情况, 建议继续对非均匀地应力场下水泥环的弹塑性应力与变形作深入研究。
-
[1] 魏东吼, 董绍华, 梁伟. 地下储气库完整性管理体系及相关技术应用研究[J]. 油气储运, 2015, 34(2): 115-121. http://yqcy.paperonce.org/oa/darticle.aspx?type=view&id=20150201 WEI D H, DONG S H, LIANG W. Applied research of integritymanagement system and related technologies of underground gasstorage[J]. Oil & Gas Storage and Transportation, 2015, 34(2): 115-121. http://yqcy.paperonce.org/oa/darticle.aspx?type=view&id=20150201
[2] 李丽锋, 赵新伟, 罗金恒, 等. 盐穴地下储气库失效分析与预防措施[J]. 油气储运, 2010, 29(6): 407-410. http://yqcy.paperonce.org/oa/darticle.aspx?type=view&id=20100602 LI L F, ZHAO X W, LUO J H, et al. Failure analysis ofunderground salt cavern gas storage and countermeasures[J]. Oil & Gas Storage and Transportation, 2010, 29(6): 407-410. http://yqcy.paperonce.org/oa/darticle.aspx?type=view&id=20100602
[3] 刘洋, 严海兵, 余鑫, 等. 井内压力变化对水泥环密封完整性的影响及对策[J]. 天然气工业, 2014, 34(4): 95-98. https://www.cnki.com.cn/Article/CJFDTOTAL-TRQG201404020.htm LIU Y, YAN H B, YU X, et al. Negative impacts of boreholepressure change on cement sheath sealing integrity andcountermeasures[J]. Natural Gas Industry, 2014, 34(4): 95-98. https://www.cnki.com.cn/Article/CJFDTOTAL-TRQG201404020.htm
[4] UGWU I. Cement fatigue and well integrity with applicationto life of well prediction[D]. College Station: Texas A & MUniversity, 2008: 17-18.
[5] LADVA H K J, CRASTER B, JONES T G J, et al. The cementto-formation interface in zonal isolation[J]. SPE Drilling & Completion, 2005, 20(3): 186-197.
[6] THIERCELIN M J, BEMARD D, BARET J F, et al. Cementdesign based on cement mechanical response[J]. SPE Drilling & Completion, 1998, 13(4): 266-273.
[7] 段志祥, 石坤. 水泥环和混凝土对储气井的加强与固定作用[J]. 油气储运, 2013, 32(3): 287-290. http://yqcy.paperonce.org/oa/darticle.aspx?type=view&id=201303000013 DUAN Z X, SHI K. The reinforcement and cementation ofcement sheath to gas storage well[J]. Oil & Gas Storage and Transportation, 2013, 32(3): 287-290. http://yqcy.paperonce.org/oa/darticle.aspx?type=view&id=201303000013
[8] 赵效锋, 管志川, 吴彦先, 等. 均匀地应力下水泥环应力计算及影响规律分析[J]. 石油机械, 2013, 41(9): 1-6. https://www.cnki.com.cn/Article/CJFDTOTAL-SYJI201309001.htm ZHAO X F, GUAN Z C, WU Y X, et al. Cement sheath stresscalculation and analysis of influence law under even groundstress[J]. China Petroleum Machinery, 2013, 41(9): 1-6. https://www.cnki.com.cn/Article/CJFDTOTAL-SYJI201309001.htm
[9] 艾池, 李静, 李自平, 等. CO2埋存井注入过程中水泥环应力完整性研究[J]. 科学技术与工程, 2013, 13(8): 2057-2061. https://www.cnki.com.cn/Article/CJFDTOTAL-KXJS201308008.htm AI C, LI J, LI Z P, et al. Research on cement sheath stressintegrity of CO2 buried well in the process of injection[J]. ScienceTechnology and Engineering, 2013, 13(8): 2057-2061. https://www.cnki.com.cn/Article/CJFDTOTAL-KXJS201308008.htm
[10] 张智, 许红林, 刘志伟, 等. 气井环空带压对水泥环力学完整性的影响[J]. 西南石油大学学报(自然科学版), 2016, 38(2): 155-161. https://www.cnki.com.cn/Article/CJFDTOTAL-XNSY201602020.htm ZHANG Z, XU H L, LIU Z W, et al. The effect of sustainedcasing pressure on the mechanical integrity of cement sheath ingas wells[J]. Journal of Southwest Petroleum University(Science & Technology Edition), 2016, 38(2): 155-161. https://www.cnki.com.cn/Article/CJFDTOTAL-XNSY201602020.htm
[11] 李军, 陈勉, 柳贡慧, 等. 套管、水泥环及井壁围岩组合体的弹塑性分析[J]. 石油学报, 2005, 26(6): 99-103. https://www.cnki.com.cn/Article/CJFDTOTAL-SYXB200506022.htm LI J, CHEN M, LIU G H, et al. Elastic-plastic analysis ofcasing-concrete sheath-rock combination[J]. Acta PetroleiSinica, 2005, 26(6): 99-103. https://www.cnki.com.cn/Article/CJFDTOTAL-SYXB200506022.htm
[12] 陈朝伟, 蔡永恩. 套管-地层系统套管载荷的弹塑性理论分析[J]. 石油勘探与开发, 2009, 36(2): 242-246. https://www.cnki.com.cn/Article/CJFDTOTAL-SKYK200902017.htm CHEN Z W, CAI Y E. Study on casing load in a casing-stratumsystem by elastoplastic theory[J]. Petroleum Exploration & Development, 2009, 36(2): 242-246. https://www.cnki.com.cn/Article/CJFDTOTAL-SKYK200902017.htm
[13] 唐汝众, 张士诚, 俞然刚. 高内压下固井水泥环损坏机理研究[J]. 钻采工艺, 2012, 35(1): 14-16. https://www.cnki.com.cn/Article/CJFDTOTAL-ZCGY201201007.htm TANG R Z, ZHANG S C, YU R G. Research on mechanism ofcement sheath damage under high inner pressure[J]. Drilling & Production Technology, 2012, 35(1): 14-16. https://www.cnki.com.cn/Article/CJFDTOTAL-ZCGY201201007.htm
[14] ANDRADE J D, SANGESLND S. Cement sheath failuremechanisms: Numerical estimates to design for long-term wellintegrity[J]. Journal of Petroleum Science & Engineering, 2016, 147(11): 682-698.
[15] TEODORIU C, UGWU I, SCHUBERT J. Estimation ofcasing-cement-formation interaction using a new analyticalmodel[C]. Barcelona: SPE EUROPEC/EAGE AnnualConference and Exhibition, 2010: 131335.
[16] 殷有泉, 蔡永恩, 陈朝伟, 等. 非均匀地应力场中套管载荷的理论解[J]. 石油学报, 2006, 27(4): 133-138. https://www.cnki.com.cn/Article/CJFDTOTAL-SYXB200604029.htm YIN Y Q, CAI Y E, CHEN Z W, et al. Theoretical solutionof casing loading in non-uniform ground stress field[J]. ActaPetrolei Sinica, 2006, 27(4): 133-138. https://www.cnki.com.cn/Article/CJFDTOTAL-SYXB200604029.htm
[17] 王宴滨, 高德利, 房军. 套管-水泥环-地层多层组合系统受力特性分析[J]. 应用力学学报, 2014, 31(3): 387-392, 490. https://www.cnki.com.cn/Article/CJFDTOTAL-YYLX201403016.htm WANG Y B, GAO D L, FANG J. Mechanical characteristics analysis of casing-cement ring-formation multilayer compositesystem[J]. Chinese Journal of Applied Mechanics, 2014, 31(3): 387-392, 490. https://www.cnki.com.cn/Article/CJFDTOTAL-YYLX201403016.htm
[18] LI Y, YUAN J P, QI F Z. Analysis of cemented casingmechanical failure under arbitrary in-situ stress field couplingeffects of downhole pressure and temperature[C]. Tianjin: IADC/SPE Asia Pacific Drilling Technology Conference andExhibition, 2012: 155895.
[19] 梁何生, 闻国峰, 王桂华, 等. 孔隙压力变化对地应力的影响研究[J]. 石油钻探技术, 2004, 32(2): 18-20. https://www.cnki.com.cn/Article/CJFDTOTAL-SYZT200402006.htm LIANG H S, WEN G F, WANG G H, et al. The effects of porepressure alterations on the in-situ stress[J]. Petroleum DrillingTechniques, 2004, 32(2): 18-20. https://www.cnki.com.cn/Article/CJFDTOTAL-SYZT200402006.htm
[20] 李铀, 李锶, 白世伟. 理想弹塑性厚壁圆筒弹塑性应力场再研究[J]. 岩石力学与工程学报, 2002, 21(6): 897-899. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200206031.htm LI Y, LI S, BAI S W. Re-research on the stress field of perfectelasto-plastic thick-wall tube[J]. Chinese Journal of RockMechanics and Engineering, 2002, 21(6): 897-899. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200206031.htm
[21] 王仁,熊祝华,黄文彬. 塑性力学基础[M]. 北京: 科学出版社, 1982: 208-211. WANG R, XIONG Z H, HUANG W B. Fundamentals of plasticity [M]. Beijing: Science Press, 1982: 208-211.