Numerical Simulation of the Brittleness Anisotropy of Laminated Argillaceous Limestone Facies Shale in the Jiyang Depression
-
摘要: 为了掌握济阳坳陷博兴洼陷北部沙四纯上亚段页岩脆性各向异性的规律,选用该区块樊页X井泥灰质纹层页岩岩心,以通过室内试验获得的岩心应力–应变曲线为基准标定岩石力学参数,采用三维数值模拟方法,计算分析了岩心的泊松比、弹性模量和强度参数;并采用4种典型的脆性指数计算方法,计量了所有岩心的脆性指数。研究发现:随着围压升高,各力学参数的各向异性度都呈下降趋势,且弹性参数较强度参数的各向异性对围压变化更为敏感,因此建议采用弹性参数评价岩心力学特性的各向异性;基于能量守恒原理的脆性指数计算方法计算出的脆性指数较为客观;脆性指数随围压升高而显著降低,且随着层理倾角增大,脆性指数总体上呈先降低、后升高的变化趋势,也即在与内摩擦角角度接近的方向,脆性指数最低,0°和90°取心方向上的脆性指数较高。研究结果可为济阳坳陷页岩油储层可压性评价及选井选层提供理论依据。Abstract: Argillaceous limestone facies core samples were drilled from Well Fanye X to reveal the brittleness anisotropy law of shale from the upper Chun submember of the fourth member of Shahejie Formation in the Northern Boxing Subsag of the Jiyang Depression. With the laboratory measured stress–strain curves for rocks as the calibration criterion, a three-dimensional numerical simulation was adopted to calculate and analyze the Poisson’s ratio, Young’s modulus and strength parameters of the cores. Four typical methods were employed to calculate the brittleness index of all cores. Research demonstrates that as the confining pressure increases, the degree of anisotropy of mechanical parameters declines, and the anisotropy of elastic parameters is more sensitive to the variation in confining pressure than that of strength parameters. Therefore, it is recommended to evaluate the anisotropy of core mechanics with elastic parameters. The brittleness index calculation method based on the principle of conservation of energy can objectively evaluate the shale core brittleness. The brittleness index decreases significantly as the confining pressure grows, and it first drops and then rises as the bedding dip increases. In other words, the brittleness index in the direction close to the angle of internal friction is the lowest, while it is relatively higher in the 0° and 90° coring directions. The research results can provide a theoretical basis for evaluating fracturing feasibility in shale oil reservoirs and selecting wells and layers in the Jiyang Depression.
-
Keywords:
- continental facies /
- shale /
- brittleness /
- anisotropy /
- mechanical parameters /
- Jiyang Depression
-
川深1井位于四川省阆中市,是中国石化部署在四川盆地川中隆起北部斜坡带柏垭鼻状构造的一口风险探井,主探震旦系灯影组地层含油气情况,兼探寒武系龙王庙组地层含油气情况,完钻井深8 420.00 m。该井四开采用ϕ241.3 mm钻头钻进,在6 880.00~8 060.00 m井段下入ϕ206.4 mm+ϕ193.7 mm套管,尾管下至井深8 059.50 m,尾管悬挂器位于井深6 527.20 m处。该开次钻遇二叠系、志留系、奥陶系、寒武系地层,累计钻遇14个气层,气层压力系数1.77~1.89,钻井液密度1.97 kg/L,井底静止温度172 ℃。由于地层的储集空间为裂缝和溶洞,井筒内流体与缝洞中气体的置换作用强,液柱压力不能有效压稳气层,钻进过程中油气显示活跃。该井四开井段固井作业属于典型的超深井超高温高压固井,须封住裸眼井段内的多套高压地层,为五开低压地层(压力系数为1.10~1.20)钻进及后期测试提供合格的井筒条件[1]。
高压气井固井作业时,由于水泥浆胶凝失重,环空液柱压力降低,不能压稳气层,气体会窜入水泥环与套管或水泥环与井壁之间的间隙,造成层间互窜甚至窜入井口,导致水泥环密封失效。而且在高温(>110 ℃)环境下,水泥石的强度会衰退,造成水泥石渗透率增大,导致水泥环失去密封能力。为防止气窜,国内外主要采取应用加入胶乳、纳米液硅等防气窜剂的防气窜水泥浆的方法;为防止水泥石强度衰退,主要采用在水泥中添加硅粉的方法。但对于高温高压气井,这些方法均不能有效解决气窜问题。因此,根据川深1井四开固井需求,将高温苯丙胶乳与纳米液硅结合,通过优化设计超高温高密度防气窜水泥浆,并优化固井技术措施,实现了四开超高温高压地层的有效封固,为五开低压地层钻进及后期测试提供了合格的井筒条件。
1. 固井主要技术难点
1.1 气窜风险高、防气窜难度大
1)高压气层防气窜难度大。该井四开累计钻遇14个气层,气层显示活跃。对7 507.00~7 511.00 m、7 661.00~7 666.36 m和7 731.00~7 757.00 m井段钻遇的气层进行了测试,气体上窜速度最大达到458 m/h,最大全烃值99.7%,进出口钻井液密度差0.12 kg/L;裸眼段气层压力最大达到142 MPa,气层压力梯度达到1.89 MPa/100m,潜气窜因子GFR为10.8,属于固井后环空气窜高危井。
2)重叠井段气窜风险高。该井四开固井环空返速设计为1.0 m/s,而重叠段内浆体的环空返速仅为0.49 m/s,顶替效率低,环空会残留钻井液,重叠段的固井质量通常都较差。同时,重叠段为领浆封固,为确保固井施工安全,将领浆稠化时间设计为500 min,顶部强度发展缓慢,气窜风险高。上层套管鞋以下1.50,10.00和70.00 m位置处存在3个气层,重叠段防气窜难度极大。
3)缝洞型地层气液置换效应强。气测显示明显的气层位于茅口组、洗象池群、陡坡寺组、龙王庙组、沧浪铺组、仙女洞组等地层,均为碳酸盐岩地层,裂缝和溶洞较发育,气液置换作用强,依靠液柱压力不能阻止气体进入环空。
1.2 超高温高密度水泥浆性能调整困难
1)超高温条件下高分子水泥添加剂性能会降低。超高温环境下降滤失剂、分散剂、胶乳、缓凝剂等有机高分子水泥添加剂的性能降低,水泥浆的防气窜能力、降滤失性能、高温分散性、稳定性等都得不到有效控制。以降滤失剂为例,国内目前使用的超高温降滤失剂均为AMPS共聚物,当前广泛使用的AMPS降滤失剂高温老化前后(老化条件为在180 ℃条件下养护3 h)的特性黏数测试结果见表1。由表1可知,该类降滤失剂高温老化后的特性黏数损失率均超过25.79%,这说明高分子降滤失剂在高温条件下均会发生高温降解,从而影响水泥浆的性能。
表 1 AMPS降滤失剂高温老化后的特性黏数损失率Table 1. Intrinsic viscosity loss rate of AMPS fluid loss additive after high temperature aging降滤失剂 特性黏数 特性黏数损失率,% 老化前 老化后 BS100 6.276 2.824 73.58 HS-J 5.381 2.978 62.70 SUP102L 2.124 1.776 25.79 2)水泥石需具备长期封固能力。根据API的规定,在温度达到110 ℃时,需在G级油井水泥浆中加入30%~40%的硅粉抑制水泥石强度衰退,但常规加砂水泥浆加入硅粉后并不能满足高温固井需求。室内研究表明:常规加砂水泥石在180 ℃条件下养护14 d后,用扫描电镜可以观察到针状、片状或者粒状产物,胶结物不致密(见图1);高温养护后其强度急剧衰退,降至5 MPa左右;孔隙度和渗透率增大,渗透率可增至10 mD,导致水泥环失去密封能力。同时,常规高温高密度水泥石的弹性模量达到12 GPa以上,表现为硬脆性,难以满足五开钻进及后期射孔、测试等作业对水泥环密封完整性的需求,因此,该井四开固井,要求水泥石具备一定的弹韧性[2]。
1.3 顶替效率难以保证
该井四开所用钻井液的密度为1.97 kg/L,黏度为56 s,动切力为10 Pa,且钻井液中混有3%的原油,固相含量高达34%,难以将井壁上的油膜冲洗干净;钻井液与隔离液的密度差仅为0.05 kg/L,井深7 300.00 m以深井径扩大率仅有1.05%,ϕ206.4 mm小接箍套管无法安装扶正器,套管居中度低,提高顶替效率难度大[3]。
2. 超高温水泥浆优化设计
2.1 设计思路
通过研究不同温度下水泥石强度的发展规律,优化硅粉粒径和加量来抑制水泥石强度在高温下的衰退。针对单一防气窜剂不能有效解决高压气窜的问题,将胶乳和液硅进行复配使用,利用胶乳成膜和液硅堵塞的防气窜特性,增强水泥浆防气窜能力,提高水泥石的弹性,降低其渗透率。
2.2 优化加砂参数控制强度衰退
当井底温度高于150 ℃时,常规加砂水泥石的强度快速衰退,孔隙结构变大,水化产物为以板块状C2SH与CH为主的混合物[3]。解决常规加砂水泥石强度高温快速衰退的主要方法是将硅粉加量增大到50%~60%,并优选粗硅粉和细硅粉的配比,进一步降低水泥石的钙硅(Ca/Si)比,消耗CH与高温下生成的C2SH。选取80目(粗)与200目(细)2种硅粉,将加入不同量硅粉(粗硅粉和细硅粉配比)常规加砂水泥浆形成的水泥石在180 ℃下进行养护,测其在不同养护时间下的抗压强度,结果见图2。由图2可知:加入50%粗硅粉和60%粗硅粉形成的水泥石,在养护到3 d时抗压强度达到最大,随后抗压强度出现衰退,在养护到14 d时抗压强度分别为22.2和26.7 MPa,衰退幅度分别为38.0%和32.4%;80目与200目硅粉复配可以显著抑制水泥石强度的衰退,30%粗硅粉与30%细硅粉复配,养护14 d时抗压强度为34.2 MPa,衰退幅度为4.5%。因此,选用30%粗硅粉与30%细硅粉复配。粗细硅粉复配增大了二氧化硅的比表面积,在二氧化硅加量相同的条件下可以消耗更多的C2SH与CH,其水化产物结构致密(见图3)。
2.3 复合纳米添加剂增强水泥浆的防气窜能力
1)添加胶乳改善水泥浆性能。a. 胶乳能提高水泥浆的防气窜能力。选用玻璃化温度90 ℃的苯丙胶乳提高水泥浆的高温防气窜能力,胶乳粒径为300~400 nm,其具有“成膜”防气窜和“堵塞”防气窜的作用。正压差作用下胶乳颗粒在水泥颗粒间聚集成膜,膜覆盖在滤饼表面可以阻止气窜的发生。同时,胶乳也具有颗粒堵塞作用,可降低水泥石的渗透率。b. 胶乳能改善水泥石的力学性能。胶乳颗粒具有弹性,胶乳水泥浆固化后可以显著降低水泥石的弹性模量。c. 胶乳能进一步降低水泥浆的滤失量、提高水泥浆的高温沉降稳定性[4–8]。
2)纳米液硅改善水泥浆的性能[9]。纳米液硅中含有45%活性纳米二氧化硅微球,粒径10~300 nm,中值粒径160 nm。a. 提高水泥浆的防气窜能力。纳米二氧化硅“硬球体”和胶乳“软球体”填充在水泥颗粒孔隙里,交织在一起,可以增强其堵塞能力,提高水泥浆的高温防气窜能力。b. 提高水泥石高温强度的稳定性。纳米液硅中二氧化硅的比表面积达到了25 m2/g,可以显著提高水泥石高温强度的稳定性。
2.4 超高温高密度水泥浆配方与性能评价
2.4.1 超高温高密度水泥浆配方
依据超高温高密度水泥浆设计方法及性能需求,优选抗高温的降滤失剂(高温黏度损失率小于10%)和缓凝剂(抗温200 ℃,耐温差60 ℃)等,通过优化形成了超高温高密度水泥浆配方:水泥+30.0%硅粉(80目)+30.0%硅粉(200目)+30.0%铁矿粉+2.5%降滤失剂(SCF200L)+10.0%苯丙胶乳(SCJR)+10.0%纳米液硅(SCLS)+x%缓凝剂(SCR–3)+1.0%抑泡剂(SCXP)+48.0%水。其基本性能为:密度2.05 kg/L;稠化时间可调,缓凝剂加量为7.5%和5.0%时,在155 ℃、155 MPa条件下的稠化时间分别为571 和293 min,稠化过渡时间1 min,均为直角稠化;滤失量38 mL;流动度21 cm;流性指数0.75,稠度系数1.28 mPa·sn;自由液0 mL;无沉降;72 h顶部强度15.3 MPa;水泥石在180 ℃条件下养护14 d时的强度可达41 MPa,未见强度衰退。
2.4.2 性能评价
1)防气窜性能。根据稠化时间、稠化过渡时间和API滤失量计算出超高温高密度水泥浆的防气窜系数SPN值为0.20~0.43。使用7150型防气窜模拟分析仪评价了超高温高密度水泥浆的防气窜能力,结果见图4。由图4可知,当水泥浆处于“失重”(胶凝)状态时,未见气窜现象发生。超高温高密度水泥浆形成的水泥石在180 ℃下养护7 d时的气测渗透率为0.008 1 mD,仅为常规水泥石的3%,说明苯丙胶乳和纳米液硅结合防气窜作用显著。
2)水泥石力学性能评价[10–11]。超高温高密度水泥浆形成的水泥石在180 ℃、21 MPa条件下养护7 d,单轴弹性模量为7.54 GPa,较常规水泥石降低了30%。这是由于超高温高密度水泥浆中的胶乳为有机颗粒,其在高温下具有较好的弹韧性;纳米液硅中的纳米二氧化硅颗粒填充在水泥颗粒间,通过参与水化反应,生成CSH凝胶来修补水泥石中的微观缺陷(微裂缝等),改善了水泥石的力学性能,提高了水泥石的弹性形变能力,有利于在井下温度、压力条件下及压裂施工时,保持水泥环的完整性,提高环空密封能力。
3. 现场施工
3.1 固井技术措施
针对该井四开固井防窜难度大、顶替效率低的技术难点,采用了“替净”、“压稳”和“封严”等技术措施来提高固井质量。
1)“替净”技术措施包括:a. 使用性能良好的低黏低切先导浆,占裸眼段长度2 000.00 m;b. 优化浆柱流变性能,控制隔离液屈服值大于钻井液的屈服值、小于水泥浆的屈服值;c. 设计隔离液占裸眼段长度1 500.00 m,冲洗时间25 min;d. 使用洗油型隔离液,隔离液润湿点30%,提高界面水润湿性,增强界面胶结能力。
2)“压稳”技术措施包括:a. 使用分段压稳设计模型,以确保压稳气层,压稳系数设计为1.05;b. 为确保施工过程中全程压稳,采用加重冲洗液,其密度为2.0 kg/L;c. 环空加压5 MPa候凝。3)“封严”技术措施是使用带顶部封隔器的尾管悬挂器,在水泥浆顶替到位后,坐封顶部封隔器,切断气体上窜的通道。
3.2 固井施工
注入30 m3密度1.97 kg/L的低黏低切先导浆;注入20 m3密度2.02 kg/L的加重隔离液;注入密度2.03 kg/L的领浆和尾浆,注入量分别为22和18 m3;下钻杆胶塞,替入1 m3密度2.02 kg/L的压塞液,排量0.50 m3/min;替入26 m3密度1.97 kg/L的井浆,排量1.53 m3/min;替入14 m3密度2.02 kg/L的保护液,排量0.88 m3/min;替入57 m3密度1.97 kg/L的井浆,排量1.50 m3/min;用泵车替入2.8 m3密度1.00 kg/L的清水,排量0.70 m3/min;泄压,放回水断流;下压500 kN,坐封顶部封隔器;憋压5 MPa,验封正常;起钻15柱,正循环洗井一周;关井憋压5 MPa候凝。
3.3 固井质量
候凝72 h后,通井后检测井深7 273.00 m以深固井质量,第一界面优良率94.8%,第二界面优良率为96.4%,整体固井质量达到良好。在五开钻进过程中,未发生气窜等异常现象,满足了超深井超高温高压地层的长效封固需求。
4. 结 论
1)根据胶乳水泥浆防气窜、加入硅粉控制水泥石强度衰退的思路,通过优化形成了适用于超高温高压地层固井的高密度防气窜水泥浆。该水泥浆加入30%的80目硅粉和30%的200目硅粉控制水泥石强度高温衰退,再加入10%纳米液硅改善水泥石高温强度的稳定性,形成的水泥石在180 ℃下养护14 d强度达到了41 MPa,未见强度衰退现象。
2)超高温高压高密度防气窜水泥浆利用高温苯丙胶乳“软球体”和纳米液硅“硬球体”提高水泥浆高温高压下的防气窜能力,防气窜模拟试验未见气窜现象,形成水泥石的气测渗透率仅为0.008 1 mD,单轴弹性模量为7.54 GPa。
3)川深1井四开采用超高温高压高密度防气窜水泥浆,采取“替净”、“压稳”和“封严”等技术措施,实现了对四开井眼内多个气层的有效封隔,为后期钻井、测试等作业提供了安全的井筒条件。
-
表 1 页岩细观岩心力学参数的标定值
Table 1 Calibrated meso-mechanical parameters of shale cores
力学参数 页岩 层理 弹性模量均值/GPa 56 6 单轴抗压强度均值/MPa 600 60 均值度系数 2 2 压拉比 10 15 内摩擦角/(°) 35 20 泊松比 0.21 0.23 残余强度系数 0.1 0.1 -
[1] 邹才能,董大忠,王玉满,等. 中国页岩气特征、挑战及前景 (二)[J]. 石油勘探与开发,2016,43(2):166–178. doi: 10.11698/PED.2016.02.02 ZOU Caineng, DONG Dazhong, WANG Yuman, et al. Shale gas in China: characteristics, challenges and prospects (Ⅱ)[J]. Petroleum Exploration and Development, 2016, 43(2): 166–178. doi: 10.11698/PED.2016.02.02
[2] JOHNSTON J E, CHRISTENSEN N I. Seismic anisotropy of shales[J]. Journal of Geophysical Research: Solid Earth, 1995, 100(B4): 5991–6003. doi: 10.1029/95JB00031
[3] JAEGER J C. Shear failure of anisotropic rocks[J]. Geological Magazine, 1960, 97(1): 65–72. doi: 10.1017/S0016756800061100
[4] NIANDOU H, SHAO J F, HENRY J P, et al. Laboratory investigation of the mechanical behaviour of Tournemire shale[J]. International Journal of Rock Mechanics and Mining Sciences, 1997, 34(1): 3–16. doi: 10.1016/S1365-1609(97)80029-9
[5] KUILA U, DEWHURST D N, SIGGINS A F, et al. Stress anisotropy and velocity anisotropy in low porosity shale[J]. Tectonophysics, 2011, 503(1/2): 34–44.
[6] 陈天宇,冯夏庭,张希巍,等. 黑色页岩力学特性及各向异性特性试验研究[J]. 岩石力学与工程学报,2014,33(9):1772–1779. CHEN Tianyu, FENG Xiating, ZHANG Xiwei, et al. Experimental study on mechanical and anisotropic properties of black shale[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(9): 1772–1779.
[7] 衡帅,杨春和,张保平,等. 页岩各向异性特征的试验研究[J]. 岩土力学,2015,36(3):609–616. HENG Shuai, YANG Chunhe, ZHANG Baoping, et al. Experimental research on anisotropic properties of shale[J]. Rock and Soil Mechanics, 2015, 36(3): 609–616.
[8] 王跃鹏,刘向君,梁利喜. 鄂尔多斯盆地延长组张家滩陆相页岩各向异性及能量演化特征[J]. 岩性油气藏,2019,31(5):149–160. WANG Yuepeng, LIU Xiangjun, LIANG Lixi. Anisotropy and energy evolution characteristics of Zhangjiatan continental shale of Yanchang Formation in Ordos Basin[J]. Lithologic Reservoirs, 2019, 31(5): 149–160.
[9] 洪国斌,陈勉,卢运虎,等. 川南深层页岩各向异性特征及对破裂压力的影响[J]. 石油钻探技术,2018,46(3):78–85. HONG Guobin, CHEN Mian, LU Yunhu, et al. Study on the anisotropy characteristics of deep shale in the Southern Sichuan Basin and their impacts on fracturing pressure[J]. Petroleum Drilling Techniques, 2018, 46(3): 78–85.
[10] 宁凡,邹妞妞,张大权,等. 黔北地区下寒武统牛蹄塘组页岩特征研究[J]. 特种油气藏,2020,27(1):62–67, 135. doi: 10.3969/j.issn.1006-6535.2020.01.009 NING Fan, ZOU Niuniu, ZHANG Daquan, et al. Shale characterization of the Lower Cambrian Niutitang Formation in Northern Guizhou[J]. Special Oil & Gas Reservoirs, 2020, 27(1): 62–67, 135. doi: 10.3969/j.issn.1006-6535.2020.01.009
[11] 范宇,王佳珺,刘厚彬,等. 复杂载荷下深层脆性页岩力学性能及破坏规律实验研究[J]. 特种油气藏,2020,27(4):143–148. doi: 10.3969/j.issn.1006-6535.2020.04.022 FAN Yu, WANG Jiajun, LIU Houbin, et al. Test and analysis of mechanical properties and failure patterns for the deep brittle shale under complex load[J]. Special Oil & Gas Reservoirs, 2020, 27(4): 143–148. doi: 10.3969/j.issn.1006-6535.2020.04.022
[12] HOLT R M, FJÆR E, STENEBRÅTEN J F, et al. Brittleness of shales: relevance to borehole collapse and hydraulic fracturing[J]. Journal of Petroleum Science and Engineering, 2015, 131: 200–209. doi: 10.1016/j.petrol.2015.04.006
[13] RICKMAN R, MULLEN M J, PETRE J E, et al. A practical use of shale petrophysics for stimulation design optimization: all shale plays are not clones of the Barnett Shale[R]. SPE 115258, 2008.
[14] THOMSEN L. On the use of isotropic parameters to understand anisotropic shale behavior[M]//SCHUELKE J. SEG Technical Program Expanded Abstracts 2013. Tulsa: Society of Exploration Geophysicists, 2013: 320-324.
[15] 艾池,仇德智,张军,等. 页岩力学参数测试及脆性各向异性研究[J]. 断块油气田,2017,24(5):647–651. AI Chi, QIU Dezhi, ZHANG Jun, et al. Measurements of mechanical parameters and brittleness anisotropy of shale[J]. Fault-Block Oil & Gas Field, 2017, 24(5): 647–651.
[16] 夏宏泉,梁景瑞,文晓峰. 基于CQ指标的长庆油田长6—长8段致密油储层划分标准研究[J]. 石油钻探技术,2020,48(3):114–119. doi: 10.11911/syztjs.2020064 XIA Hongquan, LIANG Jingrui, WEN Xiaofeng. The standard division of tight oil reservoirs in Chang 6-8 Members of Changqing Oilfield based on CQ index[J]. Petroleum Drilling Techniques, 2020, 48(3): 114–119. doi: 10.11911/syztjs.2020064
[17] 王敏生,光新军,耿黎东. 页岩油高效开发钻井完井关键技术及发展方向[J]. 石油钻探技术,2019,47(5):1–10. WANG Minsheng, GUANG Xinjun, GENG Lidong. Key drilling/completion technologies and development trends in the efficient development of shale oil[J]. Petroleum Drilling Techniques, 2019, 47(5): 1–10.
[18] 张伯虎, 马浩斌, 田小朋, 等.层状页岩力学参数的各向异性研究[J].地下空间与工程学报, 2020, 16(增刊2): 634–638. ZHANG Bohu, MA Haobin, TIAN Xiaopeng, et al. Deep layered shale mechanical parameters anisotropy study[J]. Chinese Journal of Underground Space and Engineering, 2020, 16(supplement 2): 634–638.
[19] 王乔,李虎,刘廷,等. 页岩脆性的表征方法及主控因素[J]. 断块油气田,2020,27(4):458–463. WANG Qiao, LI Hu, LIU Ting, et al. Characterization method and main controlling factors of shale brittleness[J]. Fault-Block Oil & Gas Field, 2020, 27(4): 458–463.
[20] 赖富强,冷寒冰,龚大建,等. 综合矿物组分和弹性力学参数的页岩脆性评价方法[J]. 断块油气田,2019,26(2):168–171, 186. LAI Fuqiang, LENG Hanbing, GONG Dajian, et al. Evaluation of shale brittleness based on mineral compositions and elastic mechanics parameters[J]. Fault-Block Oil & Gas Field, 2019, 26(2): 168–171, 186.
[21] LI L C, TANG C A, Li G, et al. Numerical simulation of 3D hydraulic fracturing based on an improved flow-stress-damage model and a parallel FEM technique[J]. Rock Mechanics and Rock Engineering, 2012, 45: 801–818.
[22] LI Lianchong, ZHAI Mingyang, ZHANG Liaoyuan, et al. Brittleness evaluation of glutenite based on energy balance and damage evolution[J]. Energies, 2019, 12(18): 3421. doi: 10.3390/en12183421
[23] RAHIMZADEH KIVI I, AMERI M, MOLLADAVOODI H. Shale brittleness evaluation based on energy balance analysis of stress-strain curves[J]. Journal of Petroleum Science and Engineering, 2018, 167: 1–19. doi: 10.1016/j.petrol.2018.03.061
[24] 何柏,谢凌志,李凤霞,等. 龙马溪页岩各向异性变形破坏特征及其机理研究[J]. 中国科学(物理学 力学 天文学),2017,47(11):103–114. HE Bo, XIE Lingzhi, LI Fengxia, et al. Anisotropic mechanism and characteristics of deformation and failure of Longmaxi Shale[J]. Scientia Sinica: Physica, Mechanica & Astronomica, 2017, 47(11): 103–114.
-
期刊类型引用(2)
1. 蒋海军,耿黎东,王晓慧,光新军. 国外石油工程碳减排技术与作业管理发展现状及启示. 石油钻探技术. 2022(05): 125-134 . 本站查看
2. 杨建朋. 深层特低渗油藏气驱油藏工程设计. 化工管理. 2021(14): 189-190 . 百度学术
其他类型引用(3)