水平井压裂过程中水泥环界面窜流模拟分析.pdf

油气田开发工程水平井压裂过程中水泥环界面窜流模拟分析陈祝兴1 李 托2 付海峰1 翁定为1 马泽元1 李 军2,3(1.中国石油勘探开发研究院 2.中国石油大学(北京)3.中国石油大学(北京)克拉玛依校区)陈祝兴?,李托,付海峰,等.水平井压裂过程中水泥环界面窜流模拟分析 J.石油机械,2023,51(9):87-93.Chen Zhuxing,Li Tuo,Fu Haifeng,et al.Simulation analysis on channeling at cement sheath interface during fractu-ring of horizontal wells J.China Petroleum Machinery,2023,51(9):87-93.摘要:现有研究不仅未定量分析水平井分段压裂过程中压裂液侵入水泥环界面的窜流长度,而且也未考虑压裂液流速等因素对界面脱黏的影响。为此,基于 Cohesive 单元方法,建立水泥环界面窜流模型,模拟压裂液侵入水泥环界面造成的窜流现象,对比水泥环-地层界面和水泥环-套管界面发生窜流的难易程度,研究水泥环弹性模量、泊松比及压裂液流速等因素对窜流长度的影响,并结合油田实例进行分析。分析结果表明:水泥环-地层界面更容易发生窜流;增大水泥环弹性模量有利于降低界面窜流的风险,水泥环泊松比对界面窜流的影响不大;在压裂施工的过程中若检测到窜流的发生,在保证压裂正常进行的同时,可以适当减小施工排量来降低窜流的风险。所得结论可为现场施工提供借鉴。关键词:水力压裂;水泥环界面;窜流模型;裂缝破裂压力;裂缝扩展压力;敏感性分析中图分类号:TE921 文献标识码:A DOI:10.16082/ki.issn.1001-4578.2023.09.012Simulation Analysis on Channeling at Cement Sheath Interface During Fracturing of Horizontal WellsChen Zhuxing1 Li Tuo2 Fu Haifeng1 Weng Dingwei1 Ma Zeyuan1 Li Jun2,3(1.PetroChina Research Institute of Petroleum Exploration and Development;2.China University of Petroleum(Beijing);3.China University of Petroleum(Beijing)at Karamay)Abstract:The available researches not only do not quantitatively analyze the channeling length of fracturing fluid entering the cement sheath interface during the staged fracturing of horizontal wells,but also do not consider the influence of factors such as fracturing fluid flow rate on debonding of the interface.Therefore,in the paper,the Cohesive unit method was used to build a model of cement sheath interface channeling to simulate the channe-ling phenomenon resulted from fracturing fluid invasion into the cement sheath interface,compare the possibility of channeling at the cement sheath-formation interface and cement sheath-casing interface,and study the influence of factors such as the elastic modulus and Poissons ratio of cement sheath and the flow rate of fracturing fluid on the channeling length;moreover,analysis was conducted combined with field case.The analysis results show that channeling more easily occurs at the cement sheath-formation interface;the increase of elastic modulus of the ce-ment sheath is favorable for reducing the risk of interfacial channeling,and the Poissons ratio of the cement sheath 782023 年 第 51 卷 第 9 期石 油 机 械CHINA PETROLEUM MACHINERY 基金项目:国家自然科学基金企业创新发展联合基金项目“海相深层高温高压钻完井工程基础理论及控制方法”(U19B6003-05);国家自然科学基金青年科学基金项目“基于压裂过程中四维地质力学的套管变形机理及控制方法研究”(52004013);国家自然科学基金青年科学基金项目“页岩气井多级压裂诱发断层滑移量化计算模型与套管变形控制方法研究”(5220040102);中国石油天然气股份有限公司勘探开发研究院横向项目“考虑固井质量的多簇水力压裂数值模拟研究”(RIPED-2022-JS-974)。has little influence on the interfacial channeling;if a channeling is detected during the fracturing,at the same time of ensuring normal operation of the fracturing,the construction displacement can be appropriately reduced to re-duce the risk of channeling.The conclusions provide reference for field construction.Keywords:hydraulic fracturing;cement sheath interface;channeling model;fracture breakdown pressure;fracture propagation pressure;sensitivity analysis0 引 言页岩气储层具有低孔隙度和低渗透率的特点1,而常规开采方法无法满足天然气的商业生产需求。水平井多级压裂作为一种能显著提高低孔、低渗致密砂岩油气藏产量的储层改造措施2,已经广泛应用于页岩气的生产中。在水平井多级压裂的过程中,必须要保证各个压裂段间不发生窜流,从而增强压裂效果3。然而由于压裂施工的影响,套管内流体压力和温度的周期性变化可能导致水泥环密封失效,并在套管-水泥或水泥-地层界面产生微环隙4,压裂液通过微环隙侵入水泥环界面,导致压裂段间产生窜流,从而减弱压裂改造效果,甚至发生压裂失效。K.E.GRAY 等5通过将套管-水泥环界面建模作为接触条件,开发了数值模型来研究水泥环界面的脱黏,该方法允许界面上存在零张力或一定量的张力传递,分别对应无黏结强度和有限黏结强度的情况。B.LECAMPION 等6通过室内试验得到了水泥环界面裂缝扩展的几何形状。WANG W.等7基于内聚力理论,建立了一个嵌入黏结带的耦合三维多孔弹性模型,分析了水力压裂过程中页岩地层和水泥环水力裂缝的扩展情况。FENG Y.及 FENG Y.C.等8-9基于垂直井段,建立了一个三维有限元模型来模拟水泥环界面处的脱黏扩展,从而将界面脱黏的长度、宽度以及周向覆盖范围进行量化。赵效锋等10针对油气井生产过程中固井界面微环隙的产生和发展规律进行仿真试验研究,并对套管-水泥环界面的黏结力进行测量,建立了固井界面微环隙的理论计算方法。彪仿俊等11使用 ABAQUS软件,采用三维流固耦合实体单元和 Cohesive 黏结单元描述岩石和裂缝的行为。从上述研究可以看出,前人的研究主要针对水泥环微环隙的产生、水泥环界面的脱黏情况,并未定量分析水平井分段压裂过程中压裂液侵入水泥环界面的窜流长度,此外,也未考虑压裂液流速等因素对界面脱黏的影响。为此,本文利用 ABAQUS软件,基于 Cohesive 单元方法,建立了水泥环界面窜流三维有限元模型,模拟压裂液侵入水泥环界面造成的窜流现象,对比水泥环-地层界面和水泥环-套管界面发生窜流的难易程度,研究水泥环弹性模量、泊松比及压裂液流速等因素对窜流长度的影响,并结合油田实例进行分析,以期为现场施工提供参考。1 数学物理模型以水泥环-地层界面为例,水泥环-套管界面原理与之相同。为了模拟流体驱动的界面脱黏,需要考虑 2 个关键因素:多孔介质的变形和孔隙流体流动;脱黏裂缝的扩展和裂缝内流体的流动。1.1 水泥环-地层界面间流体的流动方程水泥环和地层之间设置一个 Cohesive 单元,该Cohesive 单元即表示水泥环-地层的界面性质。压裂液则沿着水泥环-地层界面的 Cohesive 单元进行法向流动和切向流动。法向流动的方向垂直于 Co-hesive 单元所在的平面,切向流动的方向平行于Cohesive 单元所在的平面,如图 1 所示。图 1 水泥环-地层界面间流体流动示意图Fig.1 Schematic diagram for fluid flow at cement sheath-formation interface假设流体为不可压缩的牛顿流体,则 Cohesive单元界面中切向上单位长度的体积流量向量 qf可表示为12:qf=w312fpf(1)式中:w 为界面张开厚度,m;pf为界面中流体压力,Pa;f为界面中压裂液黏性系数,Pas。88 石 油 机 械2023 年 第 51 卷 第 9 期界面内流体法向流速为7:vt=ct(pf-pt)vb=cb(pf-pb)(2)式中:pt为界面上表面的孔隙压力,Pa;ct为界面上表面的滤失系数,m/(Pas);pb为界面下表面的孔隙压力,Pa;cb为界面下表面的滤失系数,m/(Pas);vt、vb为流体在界面上、下表面的法向流速,m/s。1.2 流固耦合作用控制方程假设水泥环和地层均为多孔介质,多孔介质中固体骨架的平衡方程为13:v-pwI()dV=ShvdS+VfvdV(3)式中:pw为孔隙压力,Pa;I 为二阶单位张量;h为表面力矩阵,N/m2;为有效力矩阵,Pa;v为虚速度矩阵,m/s;为虚应变率矩阵,s-1;f为体积力矩阵,N/m3;S 为控制体表面积,m2;V为控制体体积,m3。流体在多孔介质中的连续性方程可表示为14:1JtJwnw()+xwnwvw()=0(4)式中:nw为孔隙比,无因次;x 为空间向量,m;w为流体密度,kg/m3;vw为流体渗流速度,m/s;J 为多孔介质体积变化比率;t 为时间,s。多孔介质中的流体服从达西定律,其流速为15-16:vw=-1nwgwkpwx-wg()(5)式中:g 为重力加速度,m/s2;k 为有效渗透流量矩阵,m/s。1.3 Cohesive 单元的损伤模型根据 FENG Y.C.等9的研究,最大名义应力准则能够更加安全地预测界面间裂缝的发展。因此,本文采用最大名义应力准则判断水泥环界面的损伤。当发生损伤时,需满足的条件为:MaxTsT0s,TtT0t,TnT0n=1(6)Tn=Tn,Tn 00,Tn 0(7)式中:为麦考利括号,表示界面不会在纯压缩应力下发生损伤;Tn、Ts、Tt分别为法线方向、第一剪切方向和第二剪切方向的水泥环界面应力,Pa;T0n、T0s、T0t分别为法线方向、第一剪切方向和第二剪切方向的水泥环界面临界强度,Pa。损伤开始后界面刚度降低,当界面刚度下降为1 时,界面完全损伤。本文采用基于能量标准的BK 断裂准则来评估水泥环界面裂缝的损伤演化,具体表达式为:GCn+GCs-GCn()Gs+GtGs+Gt+Gn()=GCs+GCt+GCn(8)式中:为准则指数,无量纲;Gn、Gs、Gt分别为法向、第一剪切方向、第二剪切方向的断裂能,Pam;GCn、GCs、GCt分别为法向、第一剪切方向、第二剪切方向的临界断裂能,J/m2。2 有限元模型2.1 模型建立以水平井某压裂段为研究对象,用 ABAQUS软件建立尺寸为 2 m2 m20 m 的水泥环界面窜流模型,如图 2 所示。图 2 水泥环界面窜流模型Fig.2 Interfacial channeling model of cement sheath 模型中的水泥环、地层均为弹塑性多孔材料,套管为线弹性材料。在套管-水泥环、水泥环-地层界面插入 Cohesive 单元,Cohesive 单元的属性即代表水泥环界面的胶结性能17-18,Cohesive 单元的刚度退化即表示界面发生脱黏。当流体侵入水泥环界面时,Cohesive 单元刚度退化的长度即表示流体侵入的长度。分别在套管-水泥环、水泥环-地层界面的右侧预设 2 个破坏单元,以此表示因射孔或982023 年 第 51 卷 第 9 期陈祝兴,等:水平井压裂过程中水泥环界面窜流模拟分析 固井质量差导致的水泥环界面的初始缺陷。2 个破环单元中间预设一个流体侵入点,压裂液从该点侵入水泥环界面。采用结构化网格和过渡网格的形式对模型进行网格划分,从而提高模型的计算精度及速度。2.2 模型验证B.LECAMPION 等6通过试验得到了水泥环-地层界面裂缝脱黏的几何形状,FENG Y.C.等9的计算结果与试验对比,已经被验证是准确的。本文所建模型采用 FENG Y.C.等9的计算参数,模拟水泥环界面的脱黏过程,以验证模型的准确性。图 3 表示侵入点流体压力随时间的变化情况。从图 3 可见,流体压力迅速增加到峰值,然后下降到一个相对恒定的值。峰值压力对应于裂缝扩展的开始,通常称为裂缝破裂压力(FBP)。相对恒定的压力意味着裂缝已达到稳定扩展阶段,通常将其定义为裂缝扩展压力(FPP)。从图 3 可以看出,FENG Y.C.等9计算出的裂缝破裂压力为 20.7 MPa,裂缝扩展压力为 17.1 MPa;本文模型计算出的裂缝破裂压力为 21.6 MPa,裂缝扩展压力为17.9 MPa。经过比较可知,模型误差在 5%以内,表明本文所建模型能够较为准确地模拟水泥环界面的窜流现象。图 3 数值模型与 FENG Y.C.计算结果对比Fig.3 Comparison of numerical model and FENG Y.C.calculation results2.3 水泥环两界面窜流对比分析FENG Y.C.等9研究了水泥环-地层界面脱黏的情况,并未考虑水泥环-套管界面脱黏的情况,因此,本文在保证参数不变的情况下,在相同材料属性、相同载荷以及相同边界条件下,研究水泥环-套管界面和水泥环-地层界面发生窜流的难易程度,对比结果如图 4 所示。从图 4 可以看出:水泥环-地层界面的裂缝扩展压力为 17.1 MPa,裂缝破裂压力为 21.6 MPa;而水泥环-套管界面的裂缝扩展压力为 22.6 MPa,裂缝破裂压力为 26.8 MPa。因此,水泥环-地层界面比水泥环-套管界面更容易发生窜流。接下来以水泥环-地层界面为主进行研究。图 4 水泥环两界面窜流对比分析Fig.4 Comparative analysis on channeling at 2 interfaces of cement sheath3 实例分析3.1 概况本文以四川盆地泸州区块某井为研究对象,该井有 1 个试油层,分 26 段,段长 43109 m,平均段长 64.23 m;每段射孔 39 簇,共 150 簇,段内簇间距在 8.014.5 m 之间,平均簇间距 10.6 m。压裂资料显示,第 23 段压裂段出现了明显的压窜现象,因此以该压裂段为研究对象建立有限元分析模型。该压裂段的最大水平地应力为 95.8 MPa,最小水平地应力为83.7 MPa,垂向应力为89.9 MPa,初始地层压力为 72.88 MPa,套管内压为 77 MPa,压裂液侵入水泥环界面的流速为 110-5 m3/s。水泥塞-地层界面设置为黏结属性,内聚力刚度为8.5 GPa,黏结强度为 0.42 MPa,界面断裂能为100 J/m2,界面间流体滤失系数为 5.897 10-14 m/(sPa),流体黏度为0.001 mPas。地层、水泥环和套管的材料属性如表 1 所示。表 1 材料属性Table 1 Material properties材料弹性模量/GPa泊松比渗透率/mD孔隙度地层44.30.251.000.2水泥环10.00.170.050.2套管210.00.303.2 结果分析模型采用如下 2 个分析步:地应力平衡分析步,平衡地应力、孔隙压力以及套管内压,时间设置为 1 s;水泥环-地层界面剥离分析步,模拟压09 石 油 机 械2023 年 第 51 卷 第 9 期裂液侵入界面初始缺陷处并发生积聚,当积聚的能量超过界面的临界断裂能时,水泥环-地层界面开始脱黏,压裂液继续向前发生窜流。该分析步时间设置为 600 s,由此可计算出压裂液注入点处随该分析步时间变化的流体压力,以及界面刚度退化指数。当刚度退化指数为 1 时,界面完全脱黏,刚度退化长度即代表压裂液的窜流长度。模拟结果如图5 所示。图 5 窜流长度随时间的变化情况Fig.5 Variation of channeling length with time 图 6 表示压裂液初始侵入点处流体压力随时间的变化情况。横坐标为水泥环-地层界面剥离分析步的时间长度。从图 5 和图 6 中可以看出,水泥环-地层界面的裂缝破裂压力为 41.4 MPa,裂缝扩展压力为 36.1 MPa,压裂液窜流长度为 6.8 m,第23 压裂段的最小簇间距为 8.5 m,因此该压裂段虽然会发生窜流,但不会发生簇间窜流,能够保证压裂的正常进行。图 6 压裂液侵入点处流体压力随时间的变化曲线Fig.6 Variation curve of fluid pressure at fracturing fluid invasion point with time3.3 敏感性分析3.3.1 水泥环弹性模量图 7 和图 8 表示水泥环弹性模量对水泥环-地层界面窜流的影响。图 7 不同水泥环弹性模量下界面窜流长度Fig.7 Interfacial channeling length under different elastic moduli of cement sheath图 8 不同水泥环弹性模量下界面裂缝 扩展压力、破裂压力和窜流长度Fig.8 Interfacial fracture propagation pressure,breakdown pressure and channeling length under different elastic moduli of cement sheath 由图 7 和图 8 可以看出,当弹性模量从 10 GPa 增大到 25 GPa 时,界面的窜流长度从 14.4 m减小到 8.6 m,界面的裂缝破裂压力从 41.9 MPa增大到 55.4 MPa,界面的裂缝扩展压力从 35.1 MPa 增大到 51.1 MPa。结果表明较大的水泥环弹性模量可以增强水泥环界面的抗窜能力,这是因为水泥环的弹性模量越大,其变形能力就越低,水泥环界面张开的位移就越小。3.3.2 水泥环泊松比图 9 和图 10 表示水泥环泊松比对水泥环-地层界面窜流的影响。图 9 不同水泥环泊松比下界面窜流长度Fig.9 Interfacial channeling length under different Poissons ratios of cement sheath192023 年 第 51 卷 第 9 期陈祝兴,等:水平井压裂过程中水泥环界面窜流模拟分析 图 10 不同水泥环泊松比下界面裂缝扩 展压力、破裂压力和窜流长度Fig.10 Interfacial fracture propagation pressure,breakdown pressure and channeling length under different Poissons ratios of cement sheath由图 9 和图 10 可知,当水泥环泊松比从 0.10增加到 0.25,界面的窜流长度从 13.6 m 减小到12.4 m,界面的裂缝破裂压力从 39.9 MPa 增加到45.2 MPa,界面的裂缝扩展压力从 33.1 MPa 增加到 38.4 MPa。结果表明水泥环泊松比越大,其抗窜能力越强,但影响程度较低。3.3.3 压裂液侵入流速图 11 和图 12 表示不同压裂液侵入流速对水泥环-地层界面窜流的影响。图 11 不同压裂液侵入流速下界面窜流长度Fig.11 Interfacial channeling length under different invasion flow rates of fracturing fluid图 12 不同压裂液侵入流速下界面裂缝扩 展压力、破裂压力和窜流长度Fig.12 Interfacial fracture propagation pressure,breakdown pressure and channeling length under different invasion flow rates of fracturing fluid由图 11 和图 12 可以看出,当压裂液侵入流速从 0.5 m3/s 增大到 2.0 m3/s,界面的窜流长度从3.8 m 增加到 14.4 m,界面裂缝的破裂压力和扩展压力无明显变化。结果表明,压裂液侵入流速越大,界面的窜流长度越大,而压裂液侵入流速与排量相关。因此,在压裂施工的过程中若检测到窜流的发生,在保证压裂正常进行的同时,应适当降低排量来降低窜流的风险。4 结 论(1)利用 ABAQUS 软件,基于 Cohesive 单元方法,建立了水泥环界面窜流模型,并对比水泥环-套管界面和水泥环-地层界面的窜流情况。结果表明,水泥环-地层界面更容易发生裂缝扩展,该界面更容易发生窜流。(2)结合现场实例,某井在压裂过程中存在窜流现象,利用本文模型进行窜流分析,模拟出了压裂液的窜流长度,水泥环-地层界面剥离的裂缝扩展压力以及裂缝破裂压力,所得结论可为现场施工提供借鉴。(3)水泥环弹性模量对水泥环界面窜流的影响较大,增大水泥环的弹性模量有利于降低窜流的风险;水泥环泊松比对窜流的影响较小。(4)在压裂施工的过程中若检测到窜流的发生,可以在保证压裂正常进行的同时,适当降低施工排量来降低窜流的风险。参 考 文 献1 GAO Q,CHENG Y F,HAN S C,et al.Numerical modeling of hydraulic fracture propagation behaviors in-fluenced by pre-existing injection and production wellsJ.Journal of Petroleum Science and Engineering,2019,172:976-987.2 LIU K,GAO D L,TALEGHANI A D.Analysis on in-tegrity of cement sheath in the vertical section of wells during hydraulic fracturing J.Journal of Petroleum Science and Engineering,2018,168:370-379.3 PANJWANI S,MCDANIEL J,NIKOLAOU M.Improve-ment of zonal isolation in horizontal shale gas wells:a da-ta-driven model-based approach J.Journal of Natural Gas Science and Engineering,2017,47:101-113.4 DUSSEAULT M B,GRAY M N,NAWROCKI P A.Why oilwells leak:cement behavior and long-term con-sequences CInternational Oil and Gas Conference and Exhibition in China.Beijing,China:SPE,2000:29 石 油 机 械2023 年 第 51 卷 第 9 期SPE 64733-MS.5 GRAY K E,PODNOS E,BECKER E.Finite-element studies of near-wellbore region during cementing opera-tions:part I J.SPE Drilling&Completion,2009,24(1):127-136.6 LECAMPION B,BUNGER A,KEAR J,et al.Inter-face debonding driven by fluid injection in a cased and cemented wellbore:modeling and experiments J.International Journal of Greenhouse Gas Control,2013,18:208-223.7 WANG W,DAHI TALEGHANI A.Impact of hydraulic fracturing on cement sheath integrity;a modelling ap-proach J.Journal of Natural Gas Science and Engi-neering,2017,44:265-277.8 FENG Y,PODNOS E,GRAY K E.Well integrity a-nalysis:3D numerical modeling of cement interface debonding C50th U.S.Rock Mechanics/Geome-chanics Symposium.Houston,Texas:ARMA,2016:ARMA-2016-246.9 FENG Y C,LI X R,GRAY K E.Development of a 3D numerical model for quantifying fluid-driven interface debonding of an injector well J.International Jour-nal of Greenhouse Gas Control,2017,62:76-90.10 赵效锋,管志川,史玉才,等.固井界面微环隙产生机制及计算方法 J.中国石油大学学报(自然科学版),2017,41(5):94-101.ZHAO X F,GUAN Z C,SHI Y C,et al.An assess-ment method for occurrence of micro-annular fractures on cementing interfaces of oil and gas well casing J.Journal of China University of Petroleum(Edition of Natural Sciences),2017,41(5):94-101.11 彪仿俊,刘合,张士诚,等.水力压裂水平裂缝影响参数的数值模拟研究 J.工程力学,2011,28(10):228-235.BIAO F J,LIU H,ZHANG S C,et al.A numerical study of parameter influences on horizontal hydraulic fracture J.Engineering Mechanics,2011,28(10):228-235.12 李宗利,王亚红,任青文.自然营造力作用下岩石单裂纹水力劈裂数值仿真模型 J.岩石力学与工程学报,2007,26(4):727-733.LI Z L,WANG Y H,REN Q W.Numerical simula-tion model of hydraulic fracturing of rock with a single fracture under natural hydraulic power J.Chinese Journal of Rock Mechanics and Engineering,2007,26(4):727-733.13 朱海燕,赵军,邓金根,等.水力压裂过程中微环隙的产生和扩展的数值模拟 J.计算机辅助工程,2013,22(增刊 2):443-447.ZHU H Y,ZHAO J,DENG J G,et al.Numerical simulation on micro-annulus initiation and propagation during hydraulic fracturing J.Computer Aided En-gineering,2013,22(S2):443-447.14 张广明,刘合,张劲,等.油井水力压裂流-固耦合非线性有限元数值模拟 J.石油学报,2009,30(1):113-116.ZHANG G M,LIU H,ZHANG J,et al.Simulation of hydraulic fracturing of oil well based on fluid-solid coupling equation and non-linear finite element J.Acta Petrolei Sinica,2009,30(1):113-116.15 董平川,徐小荷.储层流固耦合的数学模型及其有限元方程 J.石油学报,1998,19(1):64-70.DONG P C,XU X H.The fully coupled mathematical model of the fluid-solid in an oil reservoir and its finite element equations J.Acta Petrolei Sinica,1998,19(1):64-70.16 蒋记伟,李军,柳贡慧,等.裸眼临时弃置井水泥塞-地层界面剥离数值模拟 J.石油机械,2020,48(8):90-96.JIANG J W,LI J,LIU G H,et al.Numerical simu-lation on cement plug/formation interface debonding of open hole temporary abandoned wells J.China Pe-troleum Machinery,2020,48(8):90-96.17 位云生,林铁军,于浩,等.基于嵌入黏聚单元法的页岩储层压裂缝网扩展规律 J.天然气工业,2022,42(10):74-83.WEI Y S,LIN T J,YU H,et al.Propagation law of fracture network in shale reservoirs based on the embe-ded cohesive unit method J.Natural Gas Industry,2022,42(10):74-83.18 欧阳勇,申昭熙,白明娜,等.基于三维模型的开窗侧钻水泥环胶结强度分析 J.石油机械,2022,50(5):16-23,31.OUYANG Y,SHEN Z X,BAI M N,et al.Analysis of cementing strength of sidetracking cement sheath based on 3D model J.China Petroleum Machinery,2022,50(5):16-23,31.第一作者简介:陈祝兴,工程师,生于 1991 年,2017年毕业于美国密苏里科技大学石油工程专业,获硕士学位,现从事非常规油气藏压裂优化设计研究工作。地址:(100083)北京市海淀区。E-mail:1164021189 。通信作者:李 托,现为在读硕士研究生,E-mail:lituo431431 。收稿日期:2023-03-08(本文编辑 杨晓峰)392023 年 第 51 卷 第 9 期陈祝兴,等:水平井压裂过程中水泥环界面窜流模拟分析