Study on crack propagation law of hydraulic fracturing of coal rock complex based on cohesive element
-
摘要:
为分析煤岩复合体中水力裂缝扩展路径及演化规律,从应力场、煤岩交界面强度以及煤层层理倾角出发,基于黏聚力单元法建立了煤岩复合体水力压裂模型,分析了不同工况下煤岩复合体水力裂缝扩展规律。研究结果表明:煤岩复合体水压裂缝在上覆岩层中起裂后,随着垂向与水平应力差的逐渐增大,水力压裂裂缝受应力场的控制作用越来越显著;煤岩交界面强度较小时,水力裂缝在交界面中难以形成水压聚集达到煤层的开裂条件,煤岩交界面强度较大时能够在煤岩交界面中扩展一定距离并完成憋压,当压力达到煤层开裂条件时水压裂缝穿越煤岩交界面并诱导煤层起裂扩展;在煤层层理倾角30°条件下水力裂缝沿主应力方向形成主裂缝并沿层理方向形成次生分支裂缝;在不同地质条件下,应力场、煤岩交界面和煤体层理对水力裂纹的控制作用不同,可结合三者对水力裂纹的控制作用使其形成复杂的裂纹网络结构。
Abstract:In order to further analyze the propagation path and evolution law of hydraulic fractures in coal rock complex, from the perspective of stress field, coal rock interface strength and coal seam bedding dip angle, a hydraulic fracture model of coal rock complex is established based on cohesion unit method, and the propagation law of hydraulic fractures in coal rock complex under different working conditions is analyzed. It is found that after the hydraulic fracture of coal rock complex cracks in the overlying strata, with the gradual increase of the stress difference between σy and σx, the hydraulic fracture is more and more significantly controlled by the stress field. When the strength of coal rock interface is small, it is difficult for hydraulic cracks to form water pressure accumulation in the interface to reach the cracking conditions of coal seams. When the strength of coal rock interface is large, it can expand a certain distance in the coal rock interface and complete pressure holding. When the pressure reaches the cracking conditions of coal seams, water pressure cracks pass through the coal rock interface and induce the crack initiation and expansion of coal seams. Under the condition of 30° dip angle of coal seam bedding, hydraulic fractures form main fractures along the direction of main stress and secondary branch fractures along the direction of bedding. Under different geological conditions, stress field, coal rock interface and coal bedding have different control effects on hydraulic cracks, so they can be effectively combined to form a complex crack network structure.
-
Keywords:
- coal rock complex /
- hydraulic fracturing /
- interface strength /
- cohesive element /
- crack propagation
-
煤层气是指储存在煤层中以甲烷为主的一种烃类气体,其在煤储层中主要以吸附态的形式存在,故煤层对甲烷吸附性能的优劣将直接影响煤层气的抽采效率[1-2]。酸化技术是强化煤层气抽采的重要措施之一,目前国内外已开展了一系列煤层酸化增透的室内实验和现场应用。苏现波等[3]发现酸液能有效溶蚀堵塞煤层的矿物、改善煤层气的导流通道。WANG等[4]研究了酸化前后煤样孔隙结构的变化,结果表明:酸化后煤样中微孔的孔体积下降,中孔和大孔的数量及孔体积明显增加。马永元等[5]测定了酸化前后煤样孔隙结构和吸附能力的变化,发现酸处理能改变煤样的孔隙组成,使煤样的吸附孔数量减少、吸附能力降低。原文杰[6]测定了不同酸化条件下煤样孔隙结构和吸附特性的变化,发现在适当的酸化条件下煤样的微孔及过渡孔具有向中、大孔转化的趋势,且酸化后煤样对甲烷吸附能力降低。
随着研究的深入,一些学者逐渐把目光聚焦在了煤的微观分子结构上,探索酸液对煤体分子结构的影响。王芳芳等[7]使用盐酸对煤样进行了酸化处理,发现酸化后煤样的富氢指数明显减小而芳构化指数和富氧指数明显增大。贺成杰等[8]使用盐酸与氢氟酸依次对煤进行处理,发现酸化后煤样的苯环取代方式发生改变、含氧官能团的含量增加。张锐等[9]使用氢氟酸对煤样进行处理,发现酸化后煤样中含氧官能团的含量上升,煤样的润湿性增强。
上述研究从微细观结构的变化剖析了酸液对煤体理化结构的影响,并通过孔隙结构的变化探究了酸液对煤体吸附特性的影响机理,但尚未进一步揭示酸化作用下微细观结构相互转化的内在机制及其对煤体吸附能力的影响。因此,本文以贵州矿区无烟煤为研究对象,使用氢氟酸处理煤样,拟通过低温氮气吸附实验、红外光谱实验(FTIR)及等温吸附实验探究酸化前后煤样微观结构和吸附能力的变化,并借助吸附势理论分析酸化作用下煤样吸附特性的演化规律,从而深入理解酸液对煤体吸附特性的影响机制。
1. 实 验
实验所用无烟煤取自贵州省金沙县林华煤矿M9煤层的1091掘进面,取样点埋深约296 m,标高+787 m。在开采面的新鲜暴露处取得所需煤块后,即用保鲜膜包裹送至实验室加工。经过破碎筛分后,选取粒径在0.18~0.25 mm范围的煤样,进行低温氮气吸附实验和等温吸附实验,并挑选粒径小于0.075 mm的煤样用于FTIR测试。
1)酸化实验。对实验煤样开展了X射线衍射实验并对谱图进行物相检索,结果表明:煤样中硅酸盐矿物含量较高,约占总矿物组分的91.60%,而碳酸盐矿物含量较低,仅占4.80%。张小东等[10]研究发现,HF可溶蚀煤中大部分矿物,酸化溶蚀效果显著,故考虑使用HF对煤样进行酸化处理。使用去离子水将质量分数为40%的HF稀释至6%,随后将煤样浸泡于6%的HF溶液中12 h,反应完毕后过滤并转入干燥箱中,在85℃的条件下干燥至恒重,处理完毕后放入试样袋中密封备用。
2)低温氮气吸附实验。实验采用3H-2000PS型静态容量法比表面积及孔结构分析仪进行低温氮气吸附实验。将煤样在90℃的条件下抽真空预处理4 h,去除孔内吸附杂质,然后用纯度大于99.99%的氮气作为吸附质进行实验,得到吸附/脱附曲线。根据实验结果,分别利用BET、BJH理论模型计算煤样的比表面积、平均孔径等孔隙结构参数。
3)FTIR实验。FTIR实验采用Nicolet 6700型傅里叶变换红外光谱仪测定,扫描范围为400~4000 cm−1,分辨率为4 cm−1,每次测定共扫描32次。取少量煤粉与溴化钾按1:100的比例充分混合研磨,将混合物在10 MPa的压力下压制30 s,成片后放入仪器进行扫描。
4)等温吸附实验。实验采用HCA-1型高压容量法瓦斯吸附装置进行等温吸附实验,实验温度为30℃,吸附质为纯度大于99.99%的甲烷气体。在吸附实验中设定6个平衡压力点,当监测到甲烷压力值在0.5 h内变化不超过0.01 MPa时,即认为达到了吸附平衡。
2. 实验结果
2.1 低温氮气吸附
酸化前后煤样吸附/脱附曲线如图1,图中p为氮气平衡压力;p0为氮气饱和蒸汽压力。由图1可知,酸化后煤样的吸附/脱附曲线的回滞现象明显小于酸化前煤样,说明经过酸处理后煤样的孔隙连通性得到改善。同时,相较于酸化前,酸化后煤样对液氮的吸附量显著降低,降幅约为36.35%。
低温氮气吸附实验可以表征煤样的孔隙发育情况,运用BET和BJH理论模型分析得到实验煤样的孔隙结构参数,酸化前后煤样孔隙结构特征见表1。
表 1 酸化前后煤样孔隙结构特征Table 1. Pore structure characteristics of coal samples before and after acidfication煤样 平均孔径/
nm比表面积/
(m2·g−1)孔体积/
(cm3·g−1)最可几孔径/
nm酸化前 2.64 45.0614 0.0026 1.48 酸化后 8.97 8.4283 0.0096 28.09 由表1可知,酸化后煤样的平均孔径、孔体积、最可几孔直径均有不同幅度的增大,分别为酸化前煤样对应值的3.40、3.69、18.98倍,而比表面积下降了81.96%。宋申等[11]指出最可几孔径代表煤中出现概率最大的孔径,可以反映孔径的分布情况。同时,戚灵灵等[12]研究发现煤样的比表面积主要由微、小孔贡献,其数量越多则比表面积越大、甲烷吸附空间越大。结合表1可以发现,酸化后煤样的平均孔径、最可几孔径、总孔体积增大而比表面积降低,说明酸化后煤样中的矿物被溶蚀,微、小孔数量大幅减少,煤样出现了显著的扩孔增容现象。
2.2 FTIR分析
研究表明,煤的微观分子结构对其甲烷吸附、解吸能力有显著影响,烷基侧链能促进煤对甲烷的吸附而含氧官能团会抑制其对甲烷的吸附[13-14]。酸化前后煤样红外光谱谱图如图2。
由图2可知,酸化后煤样在含氧官能团和羟基官能团及芳香结构波段上的吸收峰强度较酸化前大幅增加,说明在酸性条件下煤中官能团会积极参加氧化反应,从而改变煤中官能团的含量及分布。考虑某一波段的吸收峰可能是由多个吸收峰叠加而成,故需要对谱图进行分峰拟合,以深入分析酸化处理对煤分子结构的影响。根据官能团归属情况,将实验煤样的红外光谱谱图分为4部分:芳香结构、含氧官能团、脂肪结构及羟基结构,并依次对各波段进行分峰拟合[15]。
2.2.1 芳香结构分峰拟合
煤中芳香结构的特征吸收峰主要集中在700~900cm−1波段,酸化前后煤样芳香结构分峰拟合结果图3。
由图3可知,酸化后苯环单取代对应的吸收峰消失,苯环二取代所占比例下降,而苯环三取代、四取代及五取代所占比例上升,说明在酸液作用下煤中芳香结构发生了取代反应,苯环单取代和二取代逐渐向多取代转变,使得苯环上的氢原子数量减少、侧链数量增多。
2.2.2 含氧官能团分峰拟合
煤中含氧官能团的特征吸收峰主要集中在1000~1800 cm−1波段,酸化前后煤样含氧官能团分峰拟合结果如图4。研究表明,含氧官能团对煤样吸附特性影响较大:一方面,含氧官能团会占据吸附位点、减少甲烷的吸附空间[16];另一方面,含氧官能团的亲水性较强,会加剧煤中的气水竞争吸附[17]。由图4可知,酸化后煤样中的含氧官能团的种类和占比大幅增加,说明酸化过程中伴随着含氧官能团的生成与转化且酸化后煤样对甲烷的吸附能力降低。同时,XING等[18]通过分子模拟发现含氧官能团脱除后,煤晶胞中的总孔体积增加,故推测酸化作用下煤中官能团发生破坏、转化会导致微、小孔的孔径增大,从而导致实验煤样的比表面积降低、甲烷的吸附空间减少。
2.2.3 脂肪结构的谱峰拟合
煤中脂肪结构的特征吸收峰主要集中在2800~3000 cm−1波段,酸化前后煤样脂肪结构分峰拟合结果见表2。
表 2 酸化前后煤样脂肪结构分峰拟合结果Table 2. Peak fitting results of coal sample fat structure before and after acidification煤样 波数/cm−1 归属 面积占比/% 酸化前 2826.46 -CH2对称振动 1.26 2853.83 -CH2对称振动 20.19 2878.85 -CH3对称振动 8.86 2900.12 -CH不对称振动 12.85 2921.84 -CH2不对称振动 34.53 2937.38 -CH3不对称振动 10.24 2959.58 -CH3不对称振动 12.06 酸化后 2827.13 -CH2对称振动 3.08 2853.20 -CH2对称振动 20.27 2877.47 -CH3对称振动 11.24 2899.00 -CH不对称振动 15.47 2921.01 -CH2不对称振动 30.05 2935.85 -CH3不对称振动 9.68 2958.18 -CH3不对称振动 10.21 表2中各脂肪结构含量变化幅度各异,为进一步探究酸化作用下脂肪结构的变化,引入脂肪结构参数$ {\mathrm{I}}_{4} $来表征脂肪链长度,其值越大则脂肪链长度越长,计算公式如下:
$$ {\mathrm{I}}_{4}=\frac{{\mathrm{A}}_{\left({\mathrm{C}\mathrm{H}}_{2}\right)}}{{\mathrm{A}}_{\left({\mathrm{C}\mathrm{H}}_{3}\right)}} $$ (1) 式中:$ {\mathrm{A}}_{\left({\mathrm{C}\mathrm{H}}_{2}\right)} $为2920~3000 cm−1波段的-CH2含量,%;$ {\mathrm{A}}_{\left({\mathrm{C}\mathrm{H}}_{3}\right)} $为2920~3000 cm−1波段的-CH3含量,%。
酸化前后,煤样脂肪结构参数$ {\mathrm{I}}_{4} $的值分别为1.55、1.51,酸化后下降了2.43%。表明酸化后煤样的脂肪链长度缩短、支链化程度增加。结合图3、图4进行分析,认为酸化作用下连接在苯环上较长的脂肪链发生了破坏重组,生成了数量更多的含氧官能团及较短的侧链。
2.2.4 羟基结构的谱峰拟合
煤中羟基结构的特征吸收峰主要集中在3000~3600 cm−1波段,酸化前后煤样羟基结构分峰拟合结果见表3。
表 3 酸化前后煤样羟基结构分峰拟合结果Table 3. Fitting results of hydroxyl structure peaks in coal samples before and after acidification煤样 波数/cm−1 归属 面积占比/% 酸化前 3197.32 环氢键 49.68 3320.74 羟基醚氢键 21.27 3439.68 自缔合羟基氢键 19.27 3540.19 羟基$ \mathrm{\pi } $氢键 9.78 酸化后 3179.24 环氢键 36.40 3414.54 自缔合羟基氢键 53.55 3512.85 羟基π氢键 10.05 由表3可知,酸化后煤中环氢键、羟基π氢键的含量均有较大幅度的减小、羟基醚氢键甚至完全消失,而自缔合羟基氢键的含量增加至酸化前对应值的2.52倍。上述3种氢键在酸性环境下不稳定、易被破坏,且酸化处理会促进部分羟基结构的生成。
2.3 等温吸附实验
将等温吸附实验所得实验数据进行langmuir拟合,可得酸化前后实验煤样的等温吸附曲线及吸附常数a、b,酸化前后煤样的等温吸附曲线如图5。
由图5可知,实验煤样的甲烷吸附量随压力增大而增多。但随着压力不断增大,煤样的吸附速率在气压约1 MPa处骤降,并不断放缓直至吸附平衡。相较于酸化前,酸化后煤样的a、b值分别下降了14.02%和23.58%,说明酸化后煤样对甲烷的吸附量及吸附速率下降,酸液抑制了煤样的吸附能力。
2.3 吸附势特性曲线
吸附势理论是Polanyi于1914年提出的适用于物理吸附的热力学理论,可以反映物体吸附1mol吸附质所需要做的功[19]。以吸附势理论中的吸附势为纵轴、吸附相体积为横轴,即可得到吸附势特性曲线,且同一试样在不同温度下所得曲线唯一。运用吸附势理论来研究酸化前后煤样吸附特性的变化规律,可以预测某一压力条件下煤样的吸附量,进而为煤层酸化改造及煤层气抽采提供一定参考。
由吸附势理论可知,吸附势与压力和温度的表达式如下[20]:
$$ \varepsilon ={\int }_{{P}_{i}}^{{P}_{0}}\frac{RT}{P}dP=RT\mathrm{ln}\frac{{P}_{0}}{{P}_{i}} $$ (2) 式中:$ \varepsilon $为吸附势,J/mol;$ {p}_{0} $为气体饱和蒸汽压力,MPa;$ {p}_{i} $为气体吸附时的平衡压力,MPa;R为常数,取8.314 J/(mol·k);T为吸附平衡温度,取303.15 K。
由于甲烷在煤体表面吸附时的温度处于临界温度之上,故采用Dubinin建立的超临界条件下的虚拟饱和蒸汽压力的经验计算公式代替[21],即:
$$ {p}_{0}={p}_{\mathrm{c}}{\left(\frac{T}{{T}_{c}}\right)}^{2.7} $$ (3) 式中:$ {p}_{\mathrm{c}} $为甲烷的临界压力,取4.62 MPa;$ {T}_{\mathrm{c}} $为甲烷的临界温度,取190.6 K。
吸附相体积代表了气体在煤中占据的空间,计算表达式如下[22]:
$$ \omega =\frac{{M}_{g}{V}_{g}}{22\;400{\rho }_{ad}} $$ (4) 式中:$ \omega $为吸附相体积,cm3/g;$ {M}_{g} $为甲烷的摩尔质量,g/mol;$ {V}_{g} $为平衡条件下的气体吸附量,cm3/g;$ {\rho }_{ad} $为气体吸附相密度,g/cm3。
其中,气体吸附相密度$ {\rho }_{ad} $可由经验公式得[23]:
$$ {\rho }_{ad}={\rho }_{b}\mathrm{e}\mathrm{x}\mathrm{p}[-0.002\;5(T-{T}_{c}\left)\right] $$ (5) 式中:$ {\rho }_{b} $为沸点下甲烷的密度,取0.424 g/cm3。
将等温吸附实验的实验数据代入式(2)~式(5),可得到吸附势和吸附相体积。酸化前后煤样的吸附势特性曲线如图6。
由图6可知,吸附等体积甲烷时,酸化后煤样的吸附势小于酸化前煤样的对应值,说明酸化后煤样对甲烷的吸附能力降低。同时,取拟合公式的截距来表征最大吸附势和最大吸附相体积,可得酸化前最大吸附势和最大吸附相体积为1.53×104 J/mol和11.26×10−2 cm3/g,而酸化后煤样的对应值分别为1.43×104 J/mol和9.47×10−2 cm3/g,下降幅度分别为6.40%和15.86%,说明酸化后煤样对甲烷的吸附能力降低、甲烷吸附空间减少。煤中甲烷的主要吸附场所为微、小孔,酸处理后煤样的平均孔径及最可几孔径增大,说明酸化后部分微、小孔的孔径增大,印证了FTIR实验中官能团反应会导致扩孔的结论。
经酸处理后,煤样对甲烷的吸附能力降低,这主要是由酸化作用对煤样微观结构的改造所致。酸化过程中,煤中活跃的有机组分及矿物质积极参与反应,煤样的表面性质及孔隙结构发生改变,煤样的亲甲烷能力下降、甲烷的吸附空间减少,进而表现为酸化后煤样对甲烷吸附能力降低。
3. 结 论
1)酸化处理后,煤中的矿物被溶蚀,平均孔径及总孔体积增大、比表面积降低、甲烷的吸附空间减小,酸液对煤体孔隙具有明显的扩孔增容作用。
2)酸处理后,苯环结构逐渐向多取代苯环发展、部分氢键结构及脂肪侧链被破坏、含氧官能团的含量上升。
3)酸处理后,煤样的吸附常数a、b值分别下降了14.02%和23.58%,说明酸化后煤样的对甲烷的极限吸附量和吸附速率下降,酸处理能够抑制煤对甲烷的吸附作用。
4)吸附势理论表明:吸附等体积甲烷时,酸化后煤样的吸附势小于酸化前煤样的对应值,且酸化作用下实验煤样的极限吸附势和吸附相体积分别下降6.40%和15.86%,说明酸化后煤样亲甲烷能力降低、甲烷的吸附空间减少,印证了酸化作用下煤中官能团反应会导致微、小孔体积增大的结论。
-
表 1 数值计算模型主要物理力学参数表
Table 1 Table of physical and mechanical parameters of numerical calculation model
岩体 弹性模
量/GPa泊松比 孔隙比 渗透率/
10−3μm2内摩擦
角/(°)抗拉强
度/MPa煤层 3.50 0.29 0.055 4 0.039 3 27 1.13 岩层 17.50 0.25 0.055 4 0.019 3 30 4.13 交界面 2.00 0.25 0.065 4 0.058 9 25 0.73 -
[1] 柳占立,庄茁,孟庆国,等. 页岩气高效开采的力学问题与挑战[J]. 力学学报,2017,49(3):507−516. LIU Zhanli, ZHUANG Zhuo, MENG Qingguo, et al. Problems and challenges of mechanics in shale gas efficient exploitation[J]. Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(3): 507−516.
[2] 王耀锋,何学秋,王恩元,等. 水力化煤层增透技术研究进展及发展趋势[J]. 煤炭学报,2014,39(10):1945−1955. WANG Yaofeng, HE Xueqiu, WANG Enyuan, et al. Research progress and development tendency of the hydraulic technology for increasing the permeability of coal seams[J]. Journal of China Coal Society, 2014, 39(10): 1945−1955.
[3] 丁秀丽,张雨霆,张传健,等. 隧洞穿越活动断层应对措施及其适应性研究综述[J]. 隧道与地下工程灾害防治,2019,1(1):20−35. DING Xiuli, ZHANG Yuting, ZHANG Chuanjian, et al. Review on countermeasures and their adaptability evaluation to tunnels crossing active faults[J]. Hazard Control in Tunnelling and Underground Engineering, 2019, 1(1): 20−35.
[4] 种照辉,闫仑,高喜才,等. 天然裂隙储层水力裂缝扩展及渗流形态研究[J]. 采矿与安全工程学报,2020,37(6):1263−1273. CHONG Zhaohui, YAN Lun, GAO Xicai, et al. Mechanism of hydraulic fracture propagation and seepage pattern in naturally fractured reservoirs[J]. Journal of Mining & Safety Engineering, 2020, 37(6): 1263−1273.
[5] 李五忠,孙斌,孙钦平,等. 以煤系天然气开发促进中国煤层气发展的对策分析[J]. 煤炭学报,2016,41(1):67−71. LI Wuzhong, SUN Bin, SUN Qinping, et al. Analysis on coal-bed methane development based on coal measure gas in China and its countermeasure[J]. Journal of China Coal Society, 2016, 41(1): 67−71.
[6] 谢英刚,孟尚志,万欢,等. 临兴地区煤系地层多类型天然气储层地质条件分析[J]. 煤炭科学技术,2015,43(9):71−75. XIE Yinggang, MENG Shangzhi, WAN Huan, et al. Analysis on geological conditions of multi type natural gas reservoir in coal measure strata of Linxing Area[J]. Coal Science and Technology, 2015, 43(9): 71−75.
[7] 高杰,侯冰,陈勉,等. 岩性差异及界面性质对裂缝起裂扩展的影响[J]. 岩石力学与工程学报,2018,37(S2):4108−4114. GAO Jie, HOU Bing, CHEN Mian, et al. Effects of rock strength and interfacial property on fracture initiation and propagation[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(S2): 4108−4114.
[8] TAN P, Jin Y, HOU B, et al. Experiments and analysis on hydraulic sand fracturing by an improved true triaxial cell[J]. Journal of Petroleum Science and Engineering, 2017, 158: 766−774. doi: 10.1016/j.petrol.2017.09.004
[9] TAN P, Jin Y, HAN K, et al. Analysis of hydraulic fracture initiation and vertical propagation behavior in laminated shale formation[J]. Fuel, 2017, 206: 482−493. doi: 10.1016/j.fuel.2017.05.033
[10] LIU Z, JIN Y, CHEN M, et al. Analysis of non-planar multi-fracture propagation from layered-formation inclined-well hydraulic fracturing[J]. Rock Mechanics and Rock Engineering, 2016, 49: 1747−1758. doi: 10.1007/s00603-015-0872-1
[11] 马衍坤,刘泽功,周健,等. 基于孔壁应变发展规律的压裂孔三阶段起裂特征试验研究[J]. 岩土力学,2016,36(8):2151−2158. MA Yankun, LIU Zegong, ZHOU Jian, et al. Study of tri-stage fracturing characteristic of borehole based on strain of hole-wall in hydraulic fracturing process[J]. Rock and Soil Mechanics, 2016, 36(8): 2151−2158.
[12] 马衍坤,刘泽功,田厚. 煤岩水力压裂过程中钻孔应变发展特征的试验研究[J]. 安全与环境学报,2015,15(4):47−51. MA Yankun, LIU Zegong, TIAN Houqiang. Experimental study of the features of the holeboring strain development in the hydraulic fracture process[J]. Journal of Safety and Environment, 2015, 15(4): 47−51.
[13] XING P, YOSHIOKA K, ADACHI J. Lattice simulation of laboratory hydraulic fracture containment in layered reservoirs[J]. Computer and Geotechnics, 2018, 100: 62−75. doi: 10.1016/j.compgeo.2018.03.010
[14] 武鹏飞,梁卫国,廉浩杰,等. 大尺寸煤岩组合体水力裂缝越界形成缝网机理及试验研究[J]. 煤炭学报,2018,43(5):1381−1389. WU Pengfei, LIANG Weiguo, LIAN Haojie, et al. Mechanism and experimental investigation of the formation of hydro-fracture system by fracturing through the interface of large-size coal-rock[J]. Journal of China Coal Society, 2018, 43(5): 1381−1389.
[15] 付世豪,侯冰,夏阳,等. 多岩性组合层状储层一体化压裂裂缝扩展试验研究[J]. 煤炭学报,2021,46(S1):377−384. FU Shihao, HOU Bing, XIA Yang, et al. Experimental research on hydraulic fracture propagation in integrated fracturing for layered formation with multi-lithology combination[J]. Journal of China Coal Society, 2021, 46(S1): 377−384.
[16] GUO J C, LUO B, LU C, et al. Numerical investigation of hydraulic fracture propagation in a layered reservoir using the cohesive zone method[J]. Engineering Fracture Mechanics, 2017, 186: 195−207. doi: 10.1016/j.engfracmech.2017.10.013
[17] ZOU Y S, MA X F, ZHOU T, et al. Hydraulic fracture growth in a layered formation based on fracturing experiments and discrete element modeling[J]. Rock Mechanics and Rock Engineering, 2017, 50(9): 2381−2395. doi: 10.1007/s00603-017-1241-z
[18] LI L C, XIA Y J, HUANG B, et al. The behaviour of fracture growth in sedimentary rocks: a numerical study based on hydraulic fracturing processes[J]. Energies, 2016, 169(9): 1−28.
[19] 牛钦环. 含水煤样剪切破坏损伤规律试验研究[D]. 徐州: 中国矿业大学, 2017. [20] 王伸. 裂隙煤岩多级水力裂缝形成过程及机理研究[D]. 焦作: 河南理工大学, 2019. [21] WANG S, LI D, MITRI H, et al. Numerical simulation of hydraulic fracture deflection influenced by slotted directional boreholes using XFEM with a modified rock fracture energy model[J]. Journal of Petroleum Science and Engineering, 2020, 193: 107375. doi: 10.1016/j.petrol.2020.107375