涪陵页岩气井射孔及滑移段套管变形特征分析.pdf

石油管工程涪陵页岩气井射孔及滑移段套管变形特征分析易先中1,2 孙养清1 万继方3 陈 辉4 姚秀田4 吴霁薇1 张仕帆5(1.长江大学机械工程学院 2.湖北省智能油气钻采装备企校联合创新中心3.中国石油集团工程技术研究院有限公司非常规油气工程研究所 4.中国石油化工股份有限公司胜利油田分公司孤岛采油厂 5.中国地质大学未来技术学院)易先中,孙养清,万继方,等.涪陵页岩气井射孔及滑移段套管变形特征分析 J.石油机械,2023,51(9):114-123.Yi Xianzhong,Sun Yangqing,Wan Jifang,et al.Analysis on casing deformation characteristics in perforation and slip section of wells in Fuling shale gas Block J.China Petroleum Machinery,2023,51(9):114-123.摘要:涪陵区块页岩气老井在分段多级压裂后滑移段套管变形严重,将会影响后期优化改造效果,且目前对该区块套管变形特征分析较少。为此,通过结合现场测井数据,建立不同射孔参数、水泥环材料属性、页岩滑移量的套管-水泥环-地层耦合模型,并采用经典 Lame 公式验证模型的正确性。数值分析结果表明:139.70 mm10.54 mm 的 P110 无损套管受周向挤压与上覆压力时的应力变化率相近,但射孔后套管的周向挤压应力变化加快,此时套管易产生较大变形,45、90相位角射孔比 30、60、120和 180相位角射孔对套管强度的影响小;45方向页岩滑移下,套管最大等效应力位于滑移面上接触点,微小的滑移量将导致套管严重变形,页岩滑移量为4.243 mm 时,套管的总变形量为 3.957 mm;增大水泥环在页岩滑移处的弹性模量能减少套管的最大等效应力,但套管的总形变量也相对增加。水泥环弹性模量为 7 GPa 时,减小或增大泊松比至 0.26 以上均能降低套管的最大等效应力,但对预防套管变形效果不明显。设计的套管受内压地面试验装置及搭建室内试验平台,可为套管变形分析提供试验指导。关键词:涪陵页岩气;套管变形;水泥环;射孔;滑移;数值分析中图分类号:TE931 文献标识码:A DOI:10.16082/ki.issn.1001-4578.2023.09.016Analysis on Casing Deformation Characteristics in Perforation and Slip Section of Wells in Fuling Shale Gas BlockYi Xianzhong1,2 Sun Yangqing1 Wan Jifang3 Chen Hui4 Yao Xiutian4 Wu Jiwei1 Zhang Shifan5(1.School of Mechanical Engineering,Yangtze University;2.Intelligent Oil and Gas Drilling and Production Equipment Enter-prise-School Joint Innovation Center of Hubei Province;3.Unconventional Oil and Gas Engineering Research Institute,CNPC Engineer-ing Technology R&D Company Limited;4.Gudao Oil Production Plant,Sinopec Shengli Oilfield Company;5.School of Future Tech-nology,China University of Geosciences)Abstract:After sectional multistage fracturing of shale gas old wells in Fuling block,the casing seriously de-forms in the slip section,which affects the later optimization and stimulation effect;furthermore,there is little a-nalysis on the deformation characteristics of casing in the block.Therefore,combined with field well log data,a casing-cement sheath-formation coupling model with different perforating parameters,cement sheath material prop-411 石 油 机 械CHINA PETROLEUM MACHINERY2023 年 第 51 卷 第 9 期基金项目:国家自然科学基金项目“地面旋转导向钻井新方法及其控制机理”(51974035、U1262108、U1762214);湖北省重点研发计划项目“智能钻井多维动力集成钻具研发及应用”(2020BAB055);长江大学非常规油气省部共建协同创新中心开放基金(UOG2022-26、UOGBX2022-04、UOGBX2022-05)。erties and shale slippage was built,and the classical Lame formula was used to verify its correctness.The numeri-cal analysis results show that the stress variation rates of 139.70 mm10.54 mm P110 nondestructive casing under circumferential compression and overlying pressure are similar,but after perforation,the circumferential compres-sion stress of the casing changes fast,the casing is prone to significant deformation,and the impact of 45 and 90 phase angle perforation on the casing strength is smaller than that of 30,60,120 and 180 phase angle perfo-ration;under shale slip in the 45 direction,the maximum equivalent stress of the casing is located at the contact point on the slip plane,a small slippage causes severe deformation of the casing;when the shale slips for 4.243 mm,the total deformation of the casing is 3.957 mm;the increase of the elastic modulus of cement sheath can re-duce the maximum equivalent stress of the casing at the shale slip point,but the total deformation of the casing also increases correspondingly;when the elastic modulus of cement sheath is 7 GPa,the decrease or increase of the Poissons ratio to more than 0.26 can reduce the maximum equivalent stress of the casing,but it has little effect on preventing deformation of the casing.In addition,a ground experiment device for casing under internal pressure was designed,and an indoor experiment platform was built to provide experimental guidance for casing deformation analysis.Keywords:Fuling shale gas;casing deformation;cement sheath;perforation;slip;numerical analysis0 引 言2021 年我国累计开采页岩气 230 亿 m3,其中涪陵页岩气田年生产能力达 100 亿 m3,是国内第一个实现商业开发的页岩气田1。涪陵页岩气井压裂过程中套管变形问题显著,带来了页岩气开采成本增加、产量大幅降低等不良后果,严重制约着我国页岩气的开发进程2。近几年,众多学者针对套管变形问题进行了研究。刘鹏林等3认为断层滑移是套管剪切变形的主要因素,并因此建立了一种考虑断层上、下部地层的断层滑移模型及滑移量计算方法。陈朝伟等4对威远、长宁和昭通区块 9 口井的 37 处 24 臂井径测井数据进行统计处理,分析该区块的套管变形特征,结果表明剪切导致的套管变形占 60%。李皋等5通过研究页岩膨胀对套管变形的影响,认为密集射孔将释放因页岩膨胀而导致的膨胀应力,这是套管变形的主要原因。姜广彬等6建立了考虑上覆地层位移与出砂空洞的射孔套管模型,分析认为,出砂井射孔套管在轴向载荷作用下的变形不同于在径向挤压作用下的变形,微小的上覆地层就能导致射孔套管发生大变形。综上所述,大部分学者认为套管变形主要是压裂过程中产生的裂缝以及地层对套管所产生的应力等造成的,并且地质条件及材料属性影响着套管的形变特征。为了缓解涪陵区块套管变形问题,结合前人对套管变形的研究,本文以焦页 4HF 井为例,针对套管变形最为严重的射孔段及滑移段建立套管-水泥环-地层耦合模型,开展数值分析,探究套管的形变特征及其规律。1 套管变形情况图 1a 为油气开采示意图。由图 1a 可看出,油气开采阶段套管受力情况复杂,需结合现场情况进行分析。通过涪陵区块出现不同程度套管变形 8 口井的 40 臂井径、磁测厚、磁定位显示,射孔井段所在套管受压裂影响较大,套管变形严重,其中 3口井钻塞遇阻。结合测井地质参数分析认为,压裂作用导致该区块部分断层地应力改变,引起断层被激活而对套管产生剪应力7,如图 1b 所示。由图 1 可知,压裂导致套管变形的原因主要体现在以下 2 方面:(1)射孔时发出的巨大冲击力使套管承受巨大载荷,射孔周围的套管因此而产生变形和裂纹。若持续输出较大载荷,将会导致套管疲劳裂纹扩展,最终损坏套管8。套管外侧的酸液、压裂液等各种液体,会对套管产生腐蚀,也会极大地削弱材料抗挤压强度,同时这些液体还对套管施以挤压力,使套管变形损坏9。(2)套管在射孔完井后出现多处孔洞,使圆柱表面的均匀性遭到破坏,导致套管的强度降低。同时,射孔位置及其附近的套管也是地层滑移、压裂、挤压力等外力作用最强的井段,加之油井长期开采,容易造成套管损坏10。5112023 年 第 51 卷 第 9 期易先中,等:涪陵页岩气井射孔及滑移段套管变形特征分析 图 1 套管井下及滑移示意图Fig.1 Schematic diagram for downhole casing and slip 由此可见,引起套管变形的主要原因是射孔后套管强度降低,地层滑移及多次压裂等。因此探究以上情况下套管变形的特征对涪陵页岩气的高效开采具有指导意义。2 井下套管力学分析井下管柱力学行为十分复杂,必须结合实际的工程约束条件进行研究11。本文与焦页 4-4HF、10-SHF、21-5HF 井作业报告相结合展开分析计算。由 Ariys 应力函数可设非均匀应力场中套管、水泥环、地层系统的应力函数为12:=A+Blnr+Cr2+Dr2lnr+Kr2+Fr4+M/r2+H()cos2(1)式中:为系统的应力函数;A、B、C、D、F、K、M 和 H 为无因次系数;r 为径向距离,mm;为周向角,()。应力分量函数为:r=1rr+1r222=2/r2r=r1r()(2)将式(1)和式(2)联立即可得到套管、水泥环、地层的各应力分量表达式:i=Bir2+2Ci-2Ki+6Mir4+4Hir2()cos2i=-Bir2+2Ci+2Ki+6Mir4+12Fir2()cos2ri=2Ki+6Fir2-6Mir4-2Hir2()sin2(3)式中:i 取 1、2、3,分别代表套管、水泥环和地层;i为系统径向应力,MPa;i为系统周向应力,MPa;ri为系统剪切应力,MPa;Bi、Ci、Ki、Fi、Mi、Hi分别代表套管、水泥环和地层与上述力学模型有关的常数。3 数值分析模型3.1 模型建立以 139.70 mm10.54 mm 的 P110 管材为分析对象,建立套管模型进行分析,射孔参数如表1 所示。表 1 焦页 4HF 井射孔参数Table 1 Perforating parameters of Well Jiaoye 4HF分段序号射孔簇数簇长/m段长/m孔密/(孔m-1)相位角/()孔径/mm1-221.53.020609.53-1531.03.020609.5 为了便于分析计算,假设条件如下:射孔过程中射孔方向垂直套管壁面、孔眼的位置与设计位置一致;射孔之后形成的孔眼标准,未产生椭圆等其他形状,同时套管的孔眼没有裂纹和毛刺产生;射孔段及周围套管为壁厚均匀的圆柱形。套管的弹性模量 E1取 210 GPa,泊松比为0.3,屈服强度为 758 MPa;水泥环的弹性模量 E2为 7 GPa,泊松比为 0.22;围岩的弹性模量 E3取35.5 GPa,泊松比为 0.18。为了便于分析套管受不同类型地应力作用的情况,采用如图 2a 所示的矩形地层模型,分析由于地层坍塌及滑移情况下上覆地层施加于模型的垂直应力 v与水平主应力h,采用如图 2b 所示的柱形地层模型分析周向挤压 p1和内压 p2下套管的变形特征13。模型几何611 石 油 机 械2023 年 第 51 卷 第 9 期参数如下:地层为 1 m1m,水泥环外径 199.70 mm,套管外径 139.70 mm,套管内径 118.62 mm。图 2 套管-水泥环-地层组合模型Fig.2 Casing-cement sheath-formation composite model3.2 网格划分用 SolidWorks 软件建模与装配后,将所建模型导入 Workbench 中进行多区域网格划分,网格大小为 35 mm。为了提高网格质量,对水泥环网格加密为 10 mm,对套管网格加密为 5 mm。根据网格的正交质量显示,最小为 0.846 39,最大为 1,平均值为 0.973 7,满足数值分析对准确性的要求。3.3 模型验证套管模型如图 3 所示。采用弹性力学理论计算14与数值分析相结合的方式进行模型验证。图 3 套管力学分析与边界条件Fig.3 Mechanical analysis and boundary conditions of casing 圆筒内外受均布压力属于轴对称应力问题,可以采用轴对称问题的通解。结合应力边界条件和位移单值条件,可得方程的 Lame 解15:r=-1-r21/r21-r21/r22p1=-1+r21/r21-r21/r22p1(4)式中:r1为套管内径,mm;r2为套管外径,mm;p1为套管外压,MPa。为了计算套管的理论应力值,取 r 等于 r3带入公式(4),可采用薄壁圆筒方法简化计算,当 r1等于 r3时,r1趋于 0,则式(4)可化为:r1=01=-21-r21/r22p1(5)式中:r3为套管中径,mm;r1为套管径向距离为r1的正应力,MPa;1为套管径向距离为 r1的周向应力,MPa。若将 r1=59.31 mm、r2=69.85 mm、p1=20 MPa 带入式(5),可得 1=143.36 MPa。图 4 所示最大压力取 143.39 MPa,考虑到 r1/r3小于 1,故实际计算值略小于仿真计算值,验证了模型的正确性。7112023 年 第 51 卷 第 9 期易先中,等:涪陵页岩气井射孔及滑移段套管变形特征分析 图 4 套管数值分析与理论分析Fig.4 Numerical and theoretical analysis of casing4 套管变形特征分析4.1 地应力对套管变形的影响基于柱形、矩形 2 种地层受压情况,以焦页4HF 井射孔井段为例,分析套管与水泥环的最大等效应力。载荷约束如图 5 所示:A、a 分别为上覆压力与挤压力;B、b、c 为固定约束面。套管与水泥环、水泥环与地层间的摩擦因数取 1=2=0.2。在压裂车的高压作业下,地质条件将发生较大变化,取围岩压力范围为 20100 MPa,探究套管与水泥环的应力变化规律。图 6 为 40 MPa 外压下套管射孔前后等效应力云图。图 5 载荷约束Fig.5 Load constraints图 6 40 MPa 外压下套管射孔前后等效应力云图Fig.6 Cloud chart for equivalent stress before and after casing perforation under 40 MPa external pressure 如图 6a 所示,40 MPa 外压下,套管在上覆压力下呈椭圆形,孔眼处出现应力集中现象,其应力较无损套管最大应力值增大 55%;如图 6b 所示,套管在周向挤压下与未挤压交界面出现颈缩现象,射孔段最大应力较无损套管增大 82%。其应力变化情况如图 7 所示。由图 7a 可知,水泥环在射孔前后其等效应力变化较小。由图 7b 可知,无损套管受周向挤压与上覆压力的等效应力曲线相近,线性拟合斜率分别为 6.65 与 6.79,但射孔套管在周向挤压下的线性斜率为 13,此时套管易产生较大变形。4.2 射孔段套管变形特征保持 20 孔/m 的射孔密度不变,建立 30、45、60、90、120和 180共 6 种不同相位角,620 mm 8 种不同孔径的套管模型。结合现场数据,取地层挤压力为 20 MPa,套管内压 70 MPa,探究射孔参数对套管变形的影响规律。为了减小边界约束对射孔分析的影响,取套管长 2 m,并且使 1 m 射孔段位于模型中段,如图8 所示。811 石 油 机 械2023 年 第 51 卷 第 9 期图 7 射孔前后水泥环与套管等效应力曲线Fig.7 Equivalent stress curves of cement sheath and casing before and after perforation图 8 不同相位角射孔模型Fig.8 Perforation models with different phase angles 为了直观地显示套管在不同射孔角度与孔眼直径下的变形情况,隐藏水泥环与地层后,放大 250倍变形系数,取 18 mm 孔径下的套管等效应力云图,如图 9 所示。为了探究其变化规律,将所有的射孔参数通过数值仿真进行计算,在不改变信号趋势和宽度的情况下提高数据的精度,并通过 Savitzky-Golay 卷积平衡算法来平滑数据,结果如图 10 所示。在不同相位角的情况下,随着孔径的增大套管的最大等效应力随之不断增大。45、90相位角射孔下套管等效应力值最小,且当孔径为 16 mm 时,水泥环最大等效应力最小。图 9 不同相位角射孔套管等效应力云图Fig.9 Cloud chart for equivalent stress of casing perforated at different phase angles图 10 不同相位角射孔和孔径对套管与水泥环的最大等效应力曲线Fig.10 Maximum equivalent stress curves of casing and cement sheath under perforation with different phase angles and perforation sizes9112023 年 第 51 卷 第 9 期易先中,等:涪陵页岩气井射孔及滑移段套管变形特征分析 4.3 滑移段套管变形特性根据对套管变形情况的统计分析,在出现断层、天然裂缝以及层理发育区,套管变形的概率较高,占比约为 61.7%16。采用 2 m1 m1 m 的地层模型进行分析,设套管内压为 75 MPa,页岩滑移面为斜向下 45方向,如图 11 所示。为了便于分析,设置裂缝面的滑移距离 dy=dx,滑移距离分别为 0.5、1.0、1.5、2.0、2.5 和3.0 mm,共 6 组。为了便于分析滑移条件下套管的变形情况,将其 30 倍比例放大后如图 12 所示。套管在滑移面处出现变形,将严重影响后续的施工效果。由以上数值分析可知,在 45页岩滑移状况下,套管的最大等效应力位于断面上接触点,如图 13 所示。由图 13 可知,套管变形量与最大等效应力随滑移量增大而增大,沿 45 方向滑移4.243 mm,套管的最大等效应力值为 634.91 MPa,总变形量为 3.957 mm。因此预防和减小页岩滑移是防止套管形变的重要措施。图 11 页岩滑移模型Fig.11 Shale slip model图 12 45方向不同滑移量套管等效应力云图Fig.12 Cloud chart for equivalent stress of casing under different slippages in the 45 direction图 13 页岩滑移量对套管与水泥环的形变特征影响Fig.13 Deformation characteristics of shale slippage on casing and cement sheath4.4 水泥环对套管变形的影响在滑移段套管变形特征分析的基础上,通过改变水泥环性能来探究其对套管变形的影响。为了直观地展示变形特征,分析不同水泥环泊松比在 X、Y方向滑移 1 mm 时套管的最大等效应力,放大 40 倍后取不同泊松比的套管应力云图,如图 14 所示。021 石 油 机 械2023 年 第 51 卷 第 9 期图 14 不同水泥环泊松比的滑移套管等效应力云图Fig.14 Cloud chart for equivalent stress of slip casing with different Poissons ratios of cement sheath 为了减小天然裂缝、水力压裂等对套管变形的影响,仅采用提高套管钢级和加厚套管壁厚等措施效果非常有限,水泥环的材料属性很大程度也影响套管的使用寿命17。因此,设水泥环弹性模量为7 GPa,取 0.20.4 共 10 种条件下的泊松比进行数值分析,结果如图 15 所示。图 15 水泥环泊松比对套管变形特征影响曲线Fig.15 Influence curve of Poissons ratio of cement sheath on casing deformation characteristics 由图 15a 可知,随着泊松比减小套管形变量不断减小,但效果并不显著。因此,通过增大水泥环泊松比来预防套管变形不可取。由图 15b 可知,当水泥环弹性模量为 7 GPa,泊松比为 0.26 时,套管的应力值最大,减小泊松比至 0.26 以下或撤去水泥环,套管等效应力有一定的降低。水泥环弹性模量与泊松比对套管变形特征的共同作用效果如图 16 所示。图 16 水泥环弹性模量与泊松比对套管变形特征影响曲线Fig.16 Influence curve of elastic modulus and Poissons ratio of cement sheath on casing deformation characteristics 由图 16 可见,增大水泥环的弹性模量在页岩滑移处能在一定程度减小套管的最大等效应力,但因此也改变了地层力的传递系数,使套管的总变形量相对增大。弹性模量减小时,出现应力峰值时的泊松比值也相应降低。因此,需结合现场地质条件及施工要求合理选择水泥环的材料属性,以减轻套管形变与损坏的状况。5 套管变形地面试验装置为了测试受损套管承受内压时的应力应变情况,搭建了地面液压试验台。试验台主要包括CB200-5.5 超高压试压泵、井口装置、套管堵头、连接器、压力表、高压管线、封隔器和应变片等。1212023 年 第 51 卷 第 9 期易先中,等:涪陵页岩气井射孔及滑移段套管变形特征分析 试验步骤如下:将试验套管与井口装置套管接口端连接,井口装置另一端通过特定的接口与堵头连接,并通过内置封隔器将套管另外一侧堵住。井口装置两侧接口可与高压软管连接以施加液压。由于套管变形尺寸较小,可在安装完套管后在其表面贴上应变片,后续可通过示波器实时观察套管各位置的应力变化,从而推断出套管变形特征。该试验装置与连接情况如图 17 所示。图 17 套管压力试验现场连接情况Fig.17 Field connection status of casing pressure test 通过向装置施加 10、20、30、40 及 50 MPa 压力,并且经过保压后判断该装置是否满足密封性要求,试验合格后,可将该装置应用于后续的套管变形研究。6 结 论结合涪陵区块套管变形的数据及地质参数,以焦页 4-HF 井分段压裂为例,通过建立套管-水泥环-地层耦合模型,分析射孔参数、45方向页岩滑移、水泥环材料属性对套管变形的影响规律及变形特征,结论如下:(1)139.7 mm10.54 mm 的 P110 无损套管受周向挤压与上覆压力时的应力变化率相近,但射孔后,套管在周向挤压下应力变化较快,此时套管易产生较大变形。射孔套管孔眼内边缘为应力集中点,45、90相位角射孔相较于 30、60、120和180相位角射孔对套管强度的影响较小,且当孔径为 16 mm 时,水泥环所受最大等效应力最小。(2)45方向页岩滑移情况下,套管最大等效应力位于滑移面上接触点,微小的滑移量将导致套管严重变形。当页岩滑移量 4.243 mm 时,套管的总变形量为 3.957 mm。(3)增大水泥环的弹性模量时,在页岩滑移处能减小套管的最大等效应力,但套管的总变形量也相对增大。当水泥环弹性模量为 7 GPa 时,减小或增大泊松比至 0.26 以上均能有降低套管的最大等效应力,但对预防套管变形效果并不明显。(4)设计了套管受内压地面试验装置,搭建了室内试验平台,经过测试其保压效果合格后,可为套管变形分析提供试验平台保障。参 考 文 献1 HUANG Q,LIU Y,PENG J Z,et al.Utilizing shale gas drilling cuttings as admixture in cement mortars:a case study in Fuling,Chongqing,China J.Envi-ronmental Science and Pollution Research,2022,29(16):24362-24369.2 席岩,李军,梅博文,等.体积压裂过程中岩性界面对套管应力影响研究 J.石油机械,2018,46(5):115-121.XI Y,LI J,MEI B W,et al.Effect of lithology inter-face on casing stress during stimulated reservoir volume fracturing J.China Petroleum Machinery,2018,46(5):115-121.3 刘鹏林,李军,席岩,等.页岩断层滑移量计算模型及 影 响 因 素 研 究 J.石 油 机 械,2022,50(8):74-80.LIU P L,LI J,XI Y,et al.Study on calculation model of shale fault slip and its influencing factors J.China Petroleum Machinery,2022,50(8):74-80.4 陈朝伟,房超,朱勇,等.四川页岩气井套管变形特征及受力模式 J.石油机械,2020,48(2):126-134.CHEN C W,FANG C,ZHU Y,et al.Deformation characteristics and stress modes of casings for shale gas wells in Sichuan J.China Petroleum Machinery,2020,48(2):126-134.5 李皋,李泽,简旭,等.页岩膨胀应变及固井质量对套管变形的影响研究 J.特种油气藏,2021,28(3):139-143.LI G,LI Z,JIAN X,et al.Study on the effect of shale expansion strain and cementing quality on casing deformation J.Special Oil&Gas Reservoirs,2021,28(3):139-143.221 石 油 机 械2023 年 第 51 卷 第 9 期6 姜广彬,郭晓强,张瑞霞,等.出砂井射孔套管变形分析及射孔参数优化 J.济南大学学报(自然科学版),2022,36(1):86-91,99.JIANG G B,GUO X Q,ZHANG R X,et al.Deform-ation analysis of perforation casings and optimization of perforation parameters in sand production wells J.Journal of University of Jinan(Science and Technolo-gy),2022,36(1):86-91,99.7 李军,赵超杰,柳贡慧,等.页岩气压裂条件下断层滑移及其影响因素 J.中国石油大学学报(自然科学版),2021,45(2):63-70.LI J,ZHAO C J,LIU G H,et al.Assessment of fault slip in shale formation during hydraulic fracturing and its influence factors J.Journal of China University of Petroleum(Edition of Natural Sciences),2021,45(2):63-70.8 付盼.威远地区套管变形物理模拟试验研究 D.北京:中国石油大学(北京),2019.FU P.Physical simulation test study on casing deforma-tion in Weiyuan area D.Beijing:China University of Petroleum(Beijing),2019.9 李进,龚宁,李早元,等.射孔完井工况下固井水泥环破坏研究进展 J.钻井液与完井液,2016,33(6):10-16.LI J,GONG N,LI Z Y,et al.Progress in studying cement sheath failure in perforated wells J.Drilling Fluid&Completion Fluid,2016,33(6):10-16.10 邓永刚,李进,林家昱,等.套管井射孔优化设计新方法及应用 J.化工管理,2021(10):111-112,159.DENG Y G,LI J,LIN J Y,et al.Research and ap-plication of perforation optimization design method for cased wells J.Chemical Enterprise Management,2021(10):111-112,159.11 高德利,黄文君.井下管柱力学与控制方法若干研究进展 J.力学进展,2021,51(3):620-647.GAO D L,HUANG W J.Some research advances in downhole tubular mechanics and control methods J.Advances in Mechanics,2021,51(3):620-647.12 王耀锋,李军强,杨小辉.套管-水泥环-地层系统应力分布规律研究 J.石油钻探技术,2008,36(5):7-11.WANG Y F,LI J Q,YANG X H.Stress distribution in casing-cement-stratum J.Petroleum Drilling Techniques,2008,36(5):7-11.13 吴烁,房军,王宴滨,等.一种非均匀地应力下椭圆井眼套管-水泥环-地层力学特性半解析解J.科 学 技 术 与 工 程,2022,22(3):997-1003.WU S,FANG J,WANG Y B,et al.Semi analytical solution for mechanical properties of elliptical borehole casing-cement-stratum system in non-uniform in-situ stress field J.Science Technology and Engineer-ing,2022,22(3):997-1003.14 TIMOSHENKO S P,GOODIER J N.Theory of elas-ticity M.3rd ed.London:McGraw-Hill Book Company,1970.15 张皓斌.厚壁圆筒强度及开孔补强的分析与讨论J.化工设计,2020,30(6):19-24.ZHANG H B.Analysis and discussion on strength of thick wall cylinder and opening reinforcement J.Chemical Engineering Design,2020,30(6):19-24.16 廖仕孟,桑宇,宋毅,等.页岩气水平井套管变形影响段分段压裂工艺研究及现场试验 J.天然气工业,2017,37(7):40-45.LIAO S M,SANG Y,SONG Y,et al.Research and field tests of staged fracturing technology for casing de-formation section in horizontal shale gas wells J.Natural Gas Industry,2017,37(7):40-45.17 陈朝伟,石林,项德贵.长宁-威远页岩气示范区套管变形机理及对策 J.天然气工业,2016,36(11):70-75.CHEN C