储气库注采参数与管柱临界屈曲载荷分析

练章华 徐帅 丁建东 张强 潘众

练章华, 徐帅, 丁建东, 张强, 潘众. 储气库注采参数与管柱临界屈曲载荷分析[J]. 石油钻采工艺, 2018, 40(1): 131-136. doi: 10.13639/j.odpt.2018.01.021
引用本文: 练章华, 徐帅, 丁建东, 张强, 潘众. 储气库注采参数与管柱临界屈曲载荷分析[J]. 石油钻采工艺, 2018, 40(1): 131-136. doi: 10.13639/j.odpt.2018.01.021
LIAN Zhanghua, XU Shuai, DING Jiandong, ZHANG Qiang, PAN Zhong. Analysis on the injection and production parameters of underground gas storage and the critical buckling load of pipe string[J]. Oil Drilling & Production Technology, 2018, 40(1): 131-136. doi: 10.13639/j.odpt.2018.01.021
Citation: LIAN Zhanghua, XU Shuai, DING Jiandong, ZHANG Qiang, PAN Zhong. Analysis on the injection and production parameters of underground gas storage and the critical buckling load of pipe string[J]. Oil Drilling & Production Technology, 2018, 40(1): 131-136. doi: 10.13639/j.odpt.2018.01.021

储气库注采参数与管柱临界屈曲载荷分析

doi: 10.13639/j.odpt.2018.01.021
基金项目: 国家自然科学基金项目“极端条件下气井油套管端力学行为及其螺纹密封机理研究”(编号:51574198);中国石油天然气股份有限公司重大项目“华北油田储气库建设配套技术研究”(编号:2011-HB-Z220)
详细信息
    作者简介:

    练章华(1964-),1994年毕业于西南石油学院机械工程专业,获博士学位,现从事CAD/CAE/CFD、套管损坏机理、管柱力学及射孔完井等教学与科研工作,教授、博士、博士生导师。通讯地址:(610500)四川省成都市新都区西南石油大学。电话:028-83032210。E-mail:cwctlzh@swpu.edu.cn

    通讯作者:

    徐帅,西南石油大学在读研究生,主要从事管柱力学研究。通讯地址:(610500)四川省成都市新都区西南石油大学。电话:028-83032210。E-mail:527386780@qq.com

  • 中图分类号: TE931.2

Analysis on the injection and production parameters of underground gas storage and the critical buckling load of pipe string

图(7) / 表 (1)
计量
  • 文章访问数:  492
  • HTML全文浏览量:  102
  • PDF下载量:  362
  • 被引次数: 0
出版历程
  • 修回日期:  2017-12-15
  • 刊出日期:  2018-01-20

储气库注采参数与管柱临界屈曲载荷分析

doi: 10.13639/j.odpt.2018.01.021
    基金项目:  国家自然科学基金项目“极端条件下气井油套管端力学行为及其螺纹密封机理研究”(编号:51574198);中国石油天然气股份有限公司重大项目“华北油田储气库建设配套技术研究”(编号:2011-HB-Z220)
    作者简介:

    练章华(1964-),1994年毕业于西南石油学院机械工程专业,获博士学位,现从事CAD/CAE/CFD、套管损坏机理、管柱力学及射孔完井等教学与科研工作,教授、博士、博士生导师。通讯地址:(610500)四川省成都市新都区西南石油大学。电话:028-83032210。E-mail:cwctlzh@swpu.edu.cn

    通讯作者: 徐帅,西南石油大学在读研究生,主要从事管柱力学研究。通讯地址:(610500)四川省成都市新都区西南石油大学。电话:028-83032210。E-mail:527386780@qq.com
  • 中图分类号: TE931.2

摘要: 针对华北油田储气井井身结构及其采气过程,开展了储气井井口油压、产量与油管柱屈曲变形的深入研究。在前人研究的基础上,讨论了管柱在井筒中发生正弦屈曲或螺旋屈曲变形的3个临界载荷值以及这3个临界载荷之间的屈曲状态,在此基础上,得到了划分不稳定区的新临界载荷,推导出了计算油管底部轴向力的数学模型。基于带有封隔器的管柱力学分析和理论研究、封隔器管柱中的各种效应变化计算数学模型以及建立的管柱临界屈曲载荷的数学模型,用VB2013开发出了一套储气井注采管柱力学安全性评价软件,应用该软件对华北油田S-x储气井注采管柱进行了研究,为储气井井口油压、产量等参数的优选及注采管柱安全性评价提供依据。

English Abstract

练章华, 徐帅, 丁建东, 张强, 潘众. 储气库注采参数与管柱临界屈曲载荷分析[J]. 石油钻采工艺, 2018, 40(1): 131-136. doi: 10.13639/j.odpt.2018.01.021
引用本文: 练章华, 徐帅, 丁建东, 张强, 潘众. 储气库注采参数与管柱临界屈曲载荷分析[J]. 石油钻采工艺, 2018, 40(1): 131-136. doi: 10.13639/j.odpt.2018.01.021
LIAN Zhanghua, XU Shuai, DING Jiandong, ZHANG Qiang, PAN Zhong. Analysis on the injection and production parameters of underground gas storage and the critical buckling load of pipe string[J]. Oil Drilling & Production Technology, 2018, 40(1): 131-136. doi: 10.13639/j.odpt.2018.01.021
Citation: LIAN Zhanghua, XU Shuai, DING Jiandong, ZHANG Qiang, PAN Zhong. Analysis on the injection and production parameters of underground gas storage and the critical buckling load of pipe string[J]. Oil Drilling & Production Technology, 2018, 40(1): 131-136. doi: 10.13639/j.odpt.2018.01.021
  • 储气库运营中, 受储气库周期性注气、采气和由此带来的温度、压力周期性变化的影响[1], 再加上管柱腐蚀引起的强度下降, 注采管柱将长期处于周期性变化的复杂应力状态。研究管柱注采中的应力响应特征及其屈曲问题, 对于预防交替注采引起的注采管柱的安全性问题、疲劳失效等具有重要意义。

    从1950年Lubinski[2-3]首次关于垂直井管柱正弦屈曲理论方法的研究以来, 有很多学者对油管柱、钻杆柱、套管柱以及带封隔器管柱的屈曲问题进行了大量的研究工作。在Lubinski之后, Mitchell[4-5]、Cunha [5-7]、Miska[8]和Hammerlindl[9-11]等对油管柱在井筒内的屈曲形式进行了更深入的研究, 对油管柱无屈曲、正弦屈曲或螺旋屈曲等形式, 推导出了其判别准则公式。国内学者[12-20]在Lubinski[2-11]等研究成果的基础上, 对封隔器管柱的受力、应力及屈曲变形进行了系统分析和讨论。2001年, 高德利[12]考虑了封隔器作为固支端对管柱屈曲的影响, 给出了受井眼约束管柱在弯矩载荷作用下的正弦屈曲和螺旋屈曲的临界载荷。高国华[13]、冷继先[14]、董蓬勃[15]等在前人的研究基础上, 讨论了油管柱屈曲的“构形”。近几年, 宋周成、练章华、李敬元等[16-19]对高产气井管柱动力学损伤、屈曲以及冲蚀等问题也开展了大量的理论研究和相关的应用软件开发工作。2016年李子丰[20]对管柱复杂力学问题的进展进行了深入、系统地研究和论述, 取得了不少可借鉴的成果。

    在前人研究的基础上, 笔者针对储气库注采管柱的静动力学问题、屈曲问题开展深入研究, 探讨储气井中带封隔器管柱屈曲形态问题, 基于管柱力学基础理论、管柱屈曲理论建立计算数学模型, 在VB2013平台上开发储气库注采管柱安全评价软件, 为井口油压、产量等参数的优选与注采管柱安全性评价提供理论依据。

    • 管柱屈曲形态或构形的判别是研究管柱屈曲损伤的一个重要依据, Lubinski[2-3]之后, Mitchell [4-5]、Cunha J.C.[6-7]和Miska[8]等对油管柱在井筒内的屈曲形式进行了更深入地研究, 基于前人的研究成果, 管柱的屈曲形态可以划分为4种形态[17]:直线状态、正弦屈曲状态、正弦和螺旋屈曲状态、螺旋屈曲状态, 其轴向压缩临界载荷判别公式见表 1

      表 1  油管柱屈曲形态判别公式

      Table 1.  Discrimination formula of tubing string buckling form

      基于表 1的轴向压缩临界载荷判别公式, 本文将判断管柱的屈曲形态划分为Fcr1Fcr2Fcr33个临界值, Fcr1 < Fcr2 < Fcr3, 如图 1。当管柱所受的压缩力F < Fcr1时, 管柱呈直线形态, 即没有屈曲; 当管柱所受的压缩力Fcr1 < F < Fcr2时, 管柱为正弦稳定的屈曲形态; 当管柱所受的压缩力Fcr2 < F < Fcr3时, 管柱为正弦和螺旋不稳定的屈曲形态; 当管柱所受的压缩力F > Fcr3时, 管柱为螺旋稳定的屈曲形态。

      图  1  管柱屈曲形态与临界载荷示意图

      Figure 1.  Sketch of string buckling form and critical load

      设 $F_{0}=\sqrt{\frac{E I w \sin a}{r}}$, 则

      $$ \begin{array}{c} F_{\mathrm{crl}}=2 F_{0} \end{array} $$ (1)
      $$ \begin{array}{c} F_{\mathrm{cr} 2}=(2 \sqrt{2}+3.75) F_{0} / 2 \end{array} $$ (2)
      $$ \begin{array}{c} F_{\mathrm{cr} 3}=4 \sqrt{2} F_{0} \end{array} $$ (3)

      由于在实际管柱屈曲问题的分析中, 对于不稳定区管柱的屈曲形态不能确定, 给理论分析带来不便, 所以提出如果压缩载荷落在不稳定区的60%的载荷以内, 为正弦屈曲, 否则为螺旋弯曲屈曲, 其示意图见图 2。大量的研究证明[5, 12-13, 21], 该假设的百分比符合实际管柱的屈曲形态, 因此提出了一个新的临界值载荷Fcrm计算公式

      图  2  管柱屈曲形态与新临界载荷Fcrm示意图

      Figure 2.  Sketch of string buckling form and new critical load Fcrm

      $$ F_{\mathrm{crm}}=F_{\mathrm{cr} 2}+\left(F_{\mathrm{cr} 3}-F_{\mathrm{cr} 2}\right) \times 0.6=4.71 F_{0} $$ (4)

      图 2把管柱屈曲划分为了直线形态、正弦屈曲和螺旋屈曲, 主要是把不稳定区进行了划分, 为文中储气井注采管柱力学安全评价软件的开发提供了理论依据, 在现场实际工况中, 将会尽量避免不稳定区的压缩载荷, 同时也将尽量避免正弦屈曲和螺旋弯曲屈曲形态的发生, 主要目的是为了更好地避免管柱的屈曲发生, 为注采管柱现场的注采参数优化或优选、管柱安全运行及其安全性评价提供理论依据。

    • 带有封隔器油管柱受力如图 3所示, 油管柱与封隔器的关系可分为3类:自由移动、有限移动、不能移动。

      图  3  有封隔器油管柱的受力示意图

      Figure 3.  Schematic force of tubing string with packer

      本文重点研究油管柱受封隔器约束不能移动的受力情况, 此时油管底部轴向力Fa*由下式表示

      $$ F_{\mathrm{a}}^{*}=F_{\mathrm{a}}+F_{\mathrm{p}} $$ (5)

      油管底部轴向力Fa*是判断管柱屈曲形态的重要参数, 将Fa*与封隔器管柱临界载荷进行对比, 即可对管柱屈曲形态进行分类, 式中Fa由活塞效应产生, Fp是封隔器对管柱产生的力, 其计算比较复杂, 将重点分析探讨。

      在本文研究中, 封隔器对管柱的活塞效应ΔL1、螺旋弯曲效应ΔL2、鼓胀效应ΔL3和温度效应ΔL4称为4种基本效应, 其计算模型见文献[2]。由于封隔器不能移动, 在坐封封隔器时, 有时将一部分重量上提, 降低工作过程中封隔器处油管柱的压力, 尽量使整个油管柱处于拉伸状态, 即称此力为提拉力Ft, 此力一般发生在井中的压力和温度变化之前, 油管在提拉力的作用下将伸长ΔL5; 管柱中, 由于流体的流动引起的摩阻效应使管柱轴向变形为ΔL6, 因此共有6种效应, 且变形增量小于2%, 属于线弹性问题, 满足叠加原理, 在这6种效应共同作用下, 管柱的总长度变化ΔL

      $$ \Delta L=\Delta L_{1}+\Delta L_{2}+\Delta L_{3}+\Delta L_{4}+\Delta L_{5}+\Delta L_{6} $$ (6)
    • 活塞效应发生在管柱变径处, 由于受力面积变化导致管柱受力发生变化, 进而影响管柱轴向变形, 活塞效应的力学模型见图 4

      图  4  活塞效应力学模型

      Figure 4.  Mechanical model of piston effect

      模型中由下向上作用的力为

      $$ F_{\mathrm{a}}^{\prime}=\left(A_{\mathrm{p}}-A_{\mathrm{i}}\right) p_{\mathrm{i}} $$ (7)

      由上向下的力为

      $$ F_{\mathrm{a}}^{\prime \prime}=\left(A_{\mathrm{p}}-A_{\mathrm{o}}\right) p_{\mathrm{o}} $$ (8)

      其合力为

      $$ F_{\mathrm{a}}=F_{\mathrm{a}}^{\prime}+F_{\mathrm{a}}^{\prime \prime}=\left(A_{\mathrm{p}}-A_{\mathrm{i}}\right) p_{\mathrm{i}}-\left(A_{\mathrm{p}}-A_{\mathrm{o}}\right) p_{\mathrm{o}} $$ (9)
    • 关于Fp的求解, 前人主要用叠加法或者图解法, 这些方法都比较繁琐。本文采用解析法, 并用VB2013编程, 直接得到Fp。采用图 3b有限移动的封隔器, 假设去除封隔器处的限制, 用式(6)计算出在各种效应作用下油管的总长度变化ΔL。要求取Fp, 就要将油管恢复到之前在封隔器中的位置, 即ΔL产生轴向变形ΔLp, 此时需要在油管底部施加轴向力Fp, 如果Fp对油管下端的作用是压缩(Fp > 0), 那么油管长度变化为虎克定律长度变化ΔLp1与螺旋弯曲长度变化[2]ΔLp2的和; 如果Fp对油管底部是拉伸力(Fp < 0), 那么油管不会产生螺旋弯曲, 只有虎克定律长度变化, 见式(10)。

      $$ \left\{\begin{array}{l} \Delta L_{\mathrm{p}}=\Delta L_{\mathrm{p} 1}+\Delta L_{\mathrm{p} 2}=-\frac{L}{E A_{\mathrm{s}}} F_{\mathrm{p}}-\frac{r^{2}}{8 E I W} F_{\mathrm{p}}^{2} \quad\left(F_{\mathrm{p}} \gt 0\right) \\ \Delta L_{\mathrm{p}}=\Delta L_{\mathrm{p} 1}=-\frac{L}{E A_{\mathrm{s}}} F_{\mathrm{p}} \quad\left(F_{\mathrm{p}} \lt 0\right) \end{array}\right. $$ (10)

      根据式(10)可求出Fp

      $$ \left\{\begin{array}{l} F_{\mathrm{p}}=4 E I W\left[\sqrt{\left(\frac{L}{E A_{\mathrm{s}}}\right)^{2}-\frac{r^{2}}{2 E I W} \Delta L_{\mathrm{p}}}-\frac{L}{E A_{\mathrm{s}}}\right) / r^{2} \quad\left(\Delta L_{\mathrm{p}} \lt 0\right) \\ F_{\mathrm{p}}=-\frac{E A_{\mathrm{s}}}{L} \Delta L_{\mathrm{p}} \quad\left(\Delta L_{\mathrm{p}} \gt 0\right) \end{array}\right. $$ (11)

      式(11)中的ΔLp等于式(6)中各种效应作用下的总变形ΔL, 将求出的Fp代入式(5)即可得到油管柱给定工况下的油管底部轴向力Fa*, 即为表 1中的F

    • 基于上文研究结果, 用VB2013开发出了一套具有自主知识产权的储气井注采管柱力学安全性评价软件, 应用该软件对华北油田S-x储气井带封隔器管柱屈曲形态进行了分析评价。

    • 华北油田S-x井为直井, 井深4 258.2 m, 生产套管Ø177.8 mm, 注采管柱为Ø88.9 mm P110油管, 内径76 mm, 可取式封隔器位置为4 185 m, 为单一尺寸的完井管柱, 地层压力为48 MPa, 采气初期, 井口压力为29.96 MPa, 实际采气量为67×104 m3/d, 井底温度为156 ℃。天然气相对密度0.671, 黏度0.02523 mPa·s, 井底流压41.928 MPa, 井口实际温度为101.88 ℃, 环境温度–29~41 ℃。

    • 基于S-x储气井实际注采参数, 应用自主开发的软件, 计算得到该井采气时, 油管底部轴向力Fa*为105.12 kN, 临界值Fcr1为53 kN, 临界值Fcr2为87.16 kN, 过渡临界值Fcrm为124.81 kN, 临界值Fcr3为149.91 kN, 油管底部轴向力Fa*是个“+”的数值, 则油管底部受压, 且大于临界值Fcr2, 小于过渡临界值Fcrm, 即Fcr2 < Fa* < Fcrm, 管柱屈曲形态为不稳定正弦。

      图 5为油管内摩阻压降、井口温度随产量的预测变化曲线, 可以看出, 在给定井口压力29.96 MPa、油管长度4 185 m的工况下, 实际采气量67×104 m3/d时对应的井口温度为101.88 ℃, 摩阻压降为4.235 MPa。

      图  5  油管内摩阻压降、井口温度随产量的变化关系

      Figure 5.  Relationship of frictional pressure drop inside the tubing and wellhead temperature vs.production rate

      应用本文开发的软件, 预测华北油田S-x储气井不同井口油压下油管底部轴向力Fa*随产量的变化, 结果见图 6。可以看出, 在井口压力一定的工况下, 随着产量的增加, 油管底部轴向力先逐渐增加, 达到某一最大值后开始逐渐减小; 随着井口油压的增加, 油管底部轴向力逐渐减小, 曲线整体向下移动。图中红色虚线对应的纵坐标为正弦屈曲临界载荷Fcr1, 当油管底部轴向力的最大值等于正弦屈曲临界载荷53 kN时, 管柱屈曲的临界井口油压为35.68 MPa, 即在华北油田S-x储气井油管柱生产工况下, 当井口油压大于35.68 MPa时, 采气管柱不会发生屈曲变形。图中井口油压为35.68 MPa的绿色点划线为管柱不发生屈曲变形的临界曲线, 当产量小于30×104 m3/d时, 油管底部轴向力为负值, 即管柱受轴向拉力; 产量控制在30×104 m3/d左右时, 油管底部轴向力接近于0, 管柱不会发生压缩或者拉伸变形; 当产量大于30×104 m3/d时, 油管底部轴向力为正值, 即油管受到轴向压缩力。根据该曲线分析结果控制该井产量在30×104 m3/d左右, 可以使油管柱受到的轴向力最小, 当然最终要结合油田的产量和储量进行综合考虑。另外, 从图 6临界曲线还可看出, 当采气量为(75~90)×104 m3/d时, 油管柱底部的轴向力最大, 即管柱受力最恶劣, 实际优选采气量时建议尽量避开该采气量范围。

      图  6  不同井口油压下产量对油管底部轴向力的影响

      Figure 6.  Effect of production rate on the axial force at the bottom of tubing under different tubing pressures

      根据图 6结果, 拟合出华北油田S-x储气井井口油压为35.68 MPa时, 油管底部轴向力随产量变化的临界曲线关系式(12), 其相关系数R2=0.99999。

      $$ F_{\mathrm{b}}(Q)=-88.118+3.7589 Q-0.0271 Q^{2}+2.708 \times 10^{-5} Q^{3} $$ (12)

      在对管柱注采参数的分析中, 将图 6中不同井口油压下, 随产量变化的油管底部轴向力的最大值提取出来进行分析。从图 7中可以直观地分析, 当井口油压大于35.68 MPa时, 油管底部轴向力的最大值将小于正弦屈曲临界载荷, 管柱不会发生屈曲变形; 当井口油压为32.04~35.68 MPa时, 管柱处于正弦屈曲状态; 当井口油压为27.94~32.04 MPa时, 管柱处于不稳定正弦屈曲状态; 当井口油压为20.36~27.94 MPa时, 管柱处于不稳定螺旋屈曲状态; 当井口油压小于20.36 MPa时, 管柱处于螺旋屈曲状态。

      图  7  油管底部最大压力随井口油压的变化关系

      Figure 7.  Relationship of maximum pressure at the bottom of tubing vs.wellhead tubing pressure

    • (1) 油管柱在井筒中的屈曲形态分为正弦屈曲和螺旋屈曲, 基于正弦屈曲或螺旋屈曲变形的3个临界载荷值以及这3个临界载荷之间的屈曲状态, 提出了判别管柱正弦屈曲和螺旋屈曲混合的不稳定状态的临界载荷计算新公式, 基于鲁宾斯基管柱力学理论, 推导出了油管底部轴向力计算数学模型, 能够对不稳定区的管柱屈曲形态进行判断。

      (2) 基于带有封隔器的管柱力学分析和理论研究、封隔器管柱中的各种效应变化计算数学模型以及本文建立的管柱临界屈曲载荷的数学模型, 开发出了自主知识产权的储气井注采管柱力学安全性评价软件, 对S-x井注采管柱进行分析, 拟合出了油压35.68 MPa时油管底部轴向力随产量变化的临界曲线。在临界油压35.68 MPa时, 建议实际采气中尽量避开采气量75~90万m3/d, 即采气量应该低于或高于此范围为优选采气量, 对现场生产有指导意义。

      (3) 储气库注采过程是一个循环进行的过程, 注气过程中, 注入的气体是冷流体, 井口温度只有35℃, 井口处附加的温度应力很小, 对注采管柱的屈曲行为影响不大, 而采气过程中, 采出气体是在地下储集的热流体, 井口温度达到107.3 ℃, 温度过高会使管柱产生线应变并且受到井筒的约束在管柱内产生应力, 对管柱的屈曲行为有一定影响, 故本文重点强调采气过程中的采气参数对管柱屈曲行为的影响。

      符号说明:

      F为轴向压缩力, N;r为油管柱环空间隙, m;I为油管柱惯性矩, m4; E为油管柱弹性模量, MPa; w为管柱单位长度的重量, N/m; α为井斜角, °;Fa*为油管底部轴向力, N;Fa为管柱底部活塞力, N;Fp为封隔器对管柱产生的力, N;po为油管外压力, MPa; pi为油管内压, MPa; Ai为油管内截面积, mm2; Ao为油管外截面积, mm2; Ap为封隔器密封腔的横截面积, mm2; As为油管横截面积, mm2; Δpi为封隔器处油管内的压力变化, MPa; ;Δpo为封隔器处环形空间的压力变化, MPa; L为管柱长度, m。

参考文献 (21)

目录

    /

    返回文章
    返回