Study on grey feature of surface strain field during deformation and failure of coal rock mass
-
摘要:
为了表征煤岩变形破坏规律,开展了煤岩试样单轴压缩试验;采用数字散斑观测系统获取表面应变场灰度图像,基于表面应变场灰度图像,计算灰度直方图并提取灰度特征参数,分析灰度演化特征,进而表征煤岩变形破坏特征。结果表明:根据试样表面裂纹出现的时间次序将其分为缓慢破坏型、瞬间破坏型和间隔破坏型;在煤岩破坏过程中,灰度直方图均呈现由高灰度向低灰度区域移动的趋势,不同破坏类型试样的灰度直方图波峰变化存在明显差异;灰度平滑度的变化能有效反映出裂隙的演化规律,灰度能量和灰度熵值反映煤岩变形破坏过程的复杂程度和破坏程度。
Abstract:In order to characterize the deformation and failure law of coal rock, this paper carries out uniaxial compression tests of coal rock specimens, and obtains grey scale images of surface strain field by using digital speckle observation system. Based on the grey scale images of the surface strain field, the grey scale histogram was calculated and the grey scale parameters were extracted to analyze the grey scale evolution characteristics, and then characterize the deformation and damage of the coal rock. The results showed that the specimens were classified into slow failure, instantaneous failure and interval failure according to the time sequence of the appearance of cracks on the surface of the specimens. During the damage process, the grey scale histograms show a tendency to move from high grey level to low grey level, and there are obvious differences in the changes of the grey scale histogram wave peaks of the samples of different damage types; the change of grey scale smoothness can effectively reflect the evolution law of the fissure, while the grey scale energy and grey scale entropy reflect the complexity of the deformation and destruction process of the coal rock and the degree of destruction.
-
CO2-ECBM开采是一种相对高效的煤层气开采技术[1]。煤层中CO2和CH4竞争吸附的微观机制是制约CO2-ECBM的主要机制[2-3]。煤是一种复杂的大分子化合物,因此,从分子水平研究煤层中CO2和CH4的竞争吸附机制至关重要。
与其他大分子有机化合物不同,煤大分子具有分子组成多样化和化学结构复杂的特征[4-5]。因此,合理地构建煤的大分子模型是在分子水平上研究煤的结构和煤中气体的竞争吸附机理的前提条件。目前,常用X射线光电子能谱(XPS)、拉曼光谱(Raman spectroscopy)、高分辨率透射电子显微镜(HRTEM)、13C核磁共振(13C-NMR)和傅里叶变换红外光谱(FTIR)等技术对煤的分子结构进行研究和分析[6-9]。相建华等[10]基于成庄矿煤的工业分析、元素分析、13C-NMR、XPS等实验结果,构建了其大分子结构模型,并采用分子力学和分子动力学方法对所得煤结构模型进行结构优化,结果表明分子间芳香层片的定向排列是高煤级煤结构中短程有序的原因;王绍清等[11]基于HRTEM研究了煤的石墨化过程,指出石墨微晶形成过程还伴随着结构缺陷的演化;NIEKERK等[12]利用HRTEM测量了煤的芳香条纹的大小和分布,并根据13C-NMR在芳香骨架上加入硫、氮、氧、脂肪链侧链和交联键,从而建立了2种南非煤的分子模型;MOKONE等[13]利用岩相分析、元素分析、氦密度、13C-NMR、HRTEM等方法构建了含氧、氮、硫的多环芳烃分子结构,提供了更好、更准确的芳香烃测量方法。随着HRTEM等技术在煤结构分析中的应用,学者们对煤大分子碳骨架结构的认知有了更进一步的提高,这些技术的进步都有利于更全面的表征煤大分子的结构特征,从而为气体在煤中的竞争吸附模拟奠定基础。
学者们针对CO2和CH4在煤中的竞争吸附特性展开了一系列研究。ZHENG等[14]利用核磁共振手段研究了亚烟煤和无烟煤中CO2与CH4的竞争吸附特性,结果表明CO2的注入可使吸附CH4的解吸效率提高14%~26%;WEI等[15]对中煤级煤中CO2和CH4气体的竞争吸附进行了实验研究,结果表明CH4浓度的下降速率随着CO2注入压力的增加表现为先增大后趋于稳定的趋势;DANG等[16]利用分子模拟方法研究了褐煤中CO2和CH4的吸附行为,结果表明CO2对CH4的吸附选择性随着压力的增加表现为先迅速下降后趋于稳定的趋势,这说明气体压力是影响竞争吸附行为的重要因素。此外,也有学者认为气体的分子动力学直径也是CH4和CO2在煤中竞争吸附行为的主导因素[17]。虽然有部分学者采取一系列方法研究了CO2和CH4在煤中的竞争吸附行为,但对于CO2和CH4在分子水平的竞争吸附机制仍需要进一步探讨。
合理的大分子模型是模拟CO2对CH4在煤中竞争吸附性能的基础。为此,基于东曲矿8号煤的大分子结构模型开展了CO2和CH4的单组分与双组分竞争吸附研究,分析了不同吸附摩尔比条件下的吸附量和吸附选择性参数,并从相互作用能和势能分布角度揭示了竞争吸附机理,研究结果对于实现CO2置换CH4和提高煤层气采收率具有重要的意义。
1. 模型结构及模拟方法
东曲矿8号煤的工业分析与元素分析结果为:水分含量2.12%,灰分含量4.06%,挥发分含量13.24%,最大镜质组反射率1.82%,C元素含量91.52%,H元素含量4.04%,O元素含量2.46%,N元素含量1.43%;S元素含量0.56%。
通过HRTEM、13C-NMR、XPS和FTIR构建了东曲矿8号煤的大分子结构(C1258H644O24N16 S2),模型构建过程如图1。
图1中,HRTEM用于确定煤大分子的碳骨架结构,13C-NMR、XPS和FTIR用于确定脂肪结构和官能团结构的分布。通过Materials Studio软件对煤大分子结构进行分子力学与分子动力学优化,力场设置为COMPASS Ⅱ,最终得到东曲矿8号煤的最低能量构型(图1d)。由于煤中存在一些难以清除的小分子[17],其模型的密度(1.40 g/cm3)略低于实测密度(1.42 g/cm3)。在Atom Volume & Surfaces模块中运用半径为1.65 nm的圆球探针对孔隙空间进行探测(图1e),所得煤大分子中的孔隙体积为6493.21 Å3(1 Å=10−10 m)。
采用Materials Studio软件中的吸附模块进行了吸附模拟。模拟温度和最大吸附压力分别设置为298.15 K和30 MPa,根据吸附质气体压力确定气体摩尔比。为了研究CH4和CO2的竞争吸附特性,采用了单组分与不同摩尔比的双组分气体吸附方案。气体吸附模拟参数见表1。
表 1 气体吸附模拟参数Table 1. Gas adsorption simulation parameters编号 摩尔比/% 压力/MPa CH4 CO2 CH4 CO2 1 100 — 30 — 2 — 100 — 30 3 40 60 12 18 4 50 50 15 15 5 60 40 18 12 基于蒙特卡洛方法对煤大分子模型进行CH4和CO2等温吸附模拟,计算固定压力下的等温吸附量来表征气体在特定压力下的吸附性能。模拟过程采用COMPASS Ⅱ力场,平衡步数为100 000 步,生产步数为1 000 000 步,并计算平均吸附量。通过Langmuir公式获得了气体的吸附等温线[18]。朗格缪尔拟合公式为:
$$ V=a b p /(1+b p) $$ (1) 式中:V为气体吸附量,mmol/g;a为压力趋于无穷大时的最大吸附量,mmol/g;b为吸附常数,1/MPa;p为吸附压力,MPa。
2. 模拟结果
2.1 单组分气体吸附
东曲矿8号煤单组分气体的等温吸附曲线如图2,单组分等温吸附曲线的朗格缪尔拟合参数见表2。
表 2 单组分等温吸附曲线的朗格缪尔拟合参数Table 2. Langmuir fitting parameters of single component isothermal adsorption curves气体 a/(mmol·g−1) b/MPa−1 R2/% CH4 4.55 0.43 99.77 CO2 5.79 0.92 99.85 单组分气体吸附条件下CH4和CO2的气体吸附量均随着压力的增大而呈现先快速增加后稳定的趋势,符合I型等温吸附特征。单组分吸附条件下CO2的吸附量显著大于CH4的吸附量,为煤中CH4的置换提供了理论依据。CO2的吸附常数b的值也大于CH4的吸附常数,这说明CO2比CH4更容易达到最大吸附量。主要原因是CO2与煤大分子之间的作用力较强,且CO2分子之间具有较强的相互作用力[19]。
2.2 双组分竞争吸附
2.2.1 气体吸附量
双组分竞争吸附下的等温吸附曲线如图3,双组分竞争吸附等温吸附曲线的朗格缪尔拟合参数见表3。
表 3 双组分竞争吸附等温吸附曲线的朗格缪尔拟合参数Table 3. Langmuir fitting parameters of dual component competitive adsorption isothermal adsorption curves摩尔比(CO2 ∶CH4) a/(mmol·g−1) b/MPa−1 R2/% 60%∶40% 4.21 1.10 99.85 50%∶50% 4.14 1.08 99.91 40%∶60% 4.06 0.95 99.88 双组分竞争吸附下的等温吸附规律同样符合Ⅰ型等温线特征。在CO2/CH4竞争吸附体系中,总吸附量随着CO2的摩尔分数的增大而增大;当CO2的摩尔分数为60%时总吸附量达到最大值(4.21 mmol/g),当CO2的摩尔分数为40%时总吸附量最小(4.06 mmol/g)。对比3种不同摩尔比吸附条件下的吸附常数可以发现,CO2的摩尔分数越大,竞争吸附越容易达到最大吸附量。双组分气体的吸附特性与竞争吸附过程中的吸附选择性有关。
2.2.2 吸附选择性
在竞争吸附体系中,吸附选择性系数是评价气体间竞争吸附的重要指标之一。气体竞争吸附体系中CO2对CH4的吸附选择性系数$S_{{\mathrm{CO}}_2/{\mathrm{CH}}_4} $的定义如下[20]:
$$ S_{\mathrm{CO}_{2} / \mathrm{CH}_{4}}=\frac{x_{\mathrm{CO}_{2}} / x_{\mathrm{CH}_{4}}}{y_{\mathrm{CO}_{2}} / y_{\mathrm{CH}_{4}}} $$ (2) 式中:$ S_{\mathrm{CO}_{2} / \mathrm{CH}_{4}} $为CO2对CH4的吸附选择性系数;$x_{\mathrm{CO}_{2}} 、 x_{\mathrm{CH}_{4}} $为CO2和CH4在吸附相中的摩尔分数;$y_{\mathrm{CO}_{2}} 、 y_{\mathrm{CH}_{4}} $为体积相中CO2和CH4的摩尔分数。
当$ S_{\mathrm{CO}_{2} / \mathrm{CH}_{4}} $>1时,表明CO2容易吸附于煤大分子中,选择性系数越高,吸附效果就越好。CO2和CH4的在煤大分子中的双组分竞争吸附的选择性系数如图4。
对比不同摩尔分数条件下的吸附选择性系数随气体压力的变化规律可以发现:不同条件下CO2对CH4的选择性吸附系数始终大于1;当CO2的摩尔分数为60%和50%时,吸附选择性系数随着气体压力的增大而增大;当CO2的摩尔分数为40%时,吸附选择性系数随着气体压力的增大而呈减小趋势;不同摩尔比的CO2和CH4竞争吸附体系中,吸附选择性系数$S_{{\mathrm{CO}}_2/{\mathrm{CH}}_4} $随着CO2摩尔分数的减小而增大。因此,较低摩尔分数的CO2具有较强的竞争吸附性能,有利于采收率的提高。
2.3 相互作用能
相互作用能指的是吸附前后系统总能量的变化量,CH4和CO2在煤大分子中的吸附主要由范德华能和静电能引起[21]。因此,进一步分析了单组分与双组分竞争吸附条件下的相互作用能的变化。不同吸附体系的相互作用能如图5。
由图5(a)可知,单组分气体吸附条件下,相互作用能随着吸附量的增加而增加,其中,CH4吸附体系的范德华能介于−70.71 kcal/mol和−271.25 kcal/mol之间(1 kcal=4.186 kJ),CO2吸附体系的范德华能介于−68.87 kcal/mol和−254.98 kcal/mol之间,2种单组分气体吸附体系的范德华能并无明显差异;CH4吸附体系的静电能介于−0.44 kcal/mol和−2.01 kcal/mol之间,而CO2吸附体系的静电能介于−34.07 kcal/mol和−140.29 kcal/mol之间。静电能在CO2与煤大分子的吸附过程中发挥了重要的作用。气体在煤中的等量吸附热可用于判断物理吸附与化学吸附[22],与CH4吸附往较低的吸附热不同,CO2吸附往往具有较大的吸附热,这与静电作用引起的化学吸附具有密切的关系。
由图5(b)可知,双组分竞争吸附条件下,双组分竞争吸附体系的范德华能介于−15.34 kcal/mol和−295.03 kcal/mol之间,静电能介于−14.56 kcal/mol和−90.86 kcal/mol之间;不同摩尔比体系的范德华能与静电能均随着CO2摩尔分数的增加而显著增加;如总吸附量为1 mmol/g的条件下,当CO2的摩尔分数由40%增加至60%,吸附体系的范德华能由−15.63 kcal/mol增加至−72.66 kcal/mol,静电能则由−14.56 kcal/mol增加至30.23 kcal/mol;且吸附量越大,由CO2摩尔比增加导致的范德华能和静电能的增加幅度也越大。这是由于体系中总吸附量随着CO2的摩尔分数的增大而增大,从而导致相互作用能的增加。
2.4 势能分布
气体分子的势能分布有利于判别气体在煤大分子中的竞争吸附位点[23]。吸附势的绝对值越大,吸附作用就越强,吸附位点也更具优势。单组分与双组分竞争吸附条件的势能分布如图6。
由图6可知,势能分布形式基本表现为正态分布;单组分CO2的吸附势于−7.05 kcal/mol时达到优势吸附位点,高于CH4的吸附势(−3.95 kcal/mol)。煤中含有一定数量的羧基和羟基的等含氧官能团,而CO2比CH4优先吸附这些极性官能团[10],与煤大分子的相互作用更强,因此,CO2的吸附位点更具优势。
对比3种不同摩尔比条件下的势能分布可以发现:CO2的摩尔比可以改变CH4吸附势的分布,当CO2的摩尔比由60%降低至50%和40%,CH4的吸附势由−4.15 kcal/mol降低至−4.05 kcal/mol和−3.85 kcal/mol。CO2的摩尔比的降低抑制了CH4的吸附。CH4吸附势的降低有利于CO2对CH4的竞争吸附,因此,较低的CO2摩尔比有利于采收率的提高,这与吸附选择性分析中得到的结果是一致的。
3. 结 语
1)单组分吸附模拟中CO2的吸附量和吸附常数显著大于CH4的吸附量。双组分竞争吸附模拟中的总吸附量随着CO2的摩尔分数的增大而增大,吸附常数也表现出相同的规律。
2)不同摩尔比条件下的双组分吸附模拟中CO2对CH4的选择性吸附系数始终大于1。CO2摩尔分数越大,选择性吸附系数越小,较低摩尔分数的CO2具有较强的竞争吸附性能。
3)相互作用能随着吸附量的增大而显著增大,CO2吸附体系中较大的静电能促进了煤大分子对CO2的吸附,因此,不同摩尔比体系的相互作用能随着CO2摩尔分数的增加而显著增加。
4)单组分条件下CO2的吸附势大于CH4的吸附势。双组分竞争吸附条件下CH4的吸附势随着CO2摩尔比的降低而降低,这有利于CO2对CH4的竞争吸附,验证了吸附选择性分析的结果。
-
表 1 试样破坏过程中特征应力点
Table 1 Characteristic stress points in the failure process of the specimens
编号 σa/MPa σb/MPa σc/MPa σd/MPa CS-Slow 2.28 4.41 5.03 4.04 CS-Inst 1.98 4.32 5.74 5.35 CS-Disc 2.27 4.09 5.41 4.56 -
[1] 谢和平. 深部岩体力学与开采理论研究进展[J]. 煤炭学报,2019,44(5):1283−1305. XIE Heping. Research review of the state key research development program of China: Deep rock mechanics and mining theory[J]. Journal of China Coal Society, 2019, 44(5): 1283−1305.
[2] 刘向御,柴肇云,刘绪,等. 循环荷载下粉砂岩孔裂隙扩展及卸载破坏特征[J]. 煤炭学报,2022,47(S1):77−89. LIU Xiangyu, CHAI Zhaoyun, LIU Xu, et al. Pore fracture propagation and unloading failure characteristics of siltite under cyclic loads[J]. Journal of China Coal Society, 2022, 47(S1): 77−89.
[3] 毛灵涛,毕玉洁,刘海洲,等. 基于CT成像和数字体图像相关法的岩石内部变形场量测方法的研究进展[J]. 科学通报,2023,68(4):380−398. doi: 10.1360/TB-2022-0405 MAO Lingtao, BI Yujie, LIU Haizhou, et al. Progress of internal deformation measurement of rock by using CT and digital volume correlation[J]. Chinese Science Bulletin, 2023, 68(4): 380−398. doi: 10.1360/TB-2022-0405
[4] 张海东,陈永峰,赵广臣,等. 单轴压缩下预制端部节理岩桥变形破坏及裂隙扩展机制研究[J]. 煤矿安全,2021,52(9):78−84. ZHANG Haidong, CHEN Yongfeng, ZHAO Guangchen, et al. Study on deformation failure and crack propagation mechanism of prefabricated end joint rock bridge under uniaxial compression[J]. Safety in Coal Mines, 2021, 52(9): 78−84.
[5] 张天军, 景晨, 张磊, 等. 含孔试样孔周破坏的应变局部化特征[J]. 煤炭学报,2020,45(12):4087−4094 ZHANG Tianjun, JING Chen, ZHANG Lei, et al. Strain localization characteristics of perforation failure of perforated specimens[J]. Journal of China Coal Society, 2020, 45(12): 4087−4094.
[6] 刘超,靖洪文,蔚立元,等. 高压流体注入对煤岩变形和破裂特性的影响[J]. 煤炭学报,2022,47(5):2027−2040. LIU Chao, JING Hongwen, YU Liyuan, et al. Effect of high-pressure fluid injection on deformation and fracture characteristics of coal[J]. Journal of China Coal Society, 2022, 47(5): 2027−2040.
[7] 宋义敏,张悦,许海亮,等. 基于非均匀特征的岩石蠕滑与黏滑变形演化研究[J]. 岩土力学,2020,41(2):363−371. SONG Yimin, ZHANG Yue, XU Hailiang, et al. Study on creep-slip and stick-slip deformation evolution of rock based on non-uniform characteristics[J]. Rock and Soil Mechanics, 2020, 41(2): 363−371.
[8] 潘红宇,葛迪,张天军,等. 应变率对岩石裂隙扩展规律的影响[J]. 煤炭学报,2018,43(3):675−683. PAN Hongyu, GE Di, ZHANG Tianjun, et al. Influence of strain rate on the rock fracture propagation law[J]. Journal of China Coal Society, 2018, 43(3): 675−683.
[9] ZHANG T J, ZHANG L, LI S G, et al. Stress inversion of coal with a gas drilling borehole and the law of crack propagation[J]. Energies, 2017, 10(11): 1743. doi: 10.3390/en10111743
[10] 付巍. 基于数字图像技术的相似材料试样破坏试验研究[J]. 煤矿安全,2021,52(8):39−43. FU Wei. Destruction test of similar material sample based on digital image technology[J]. Safety in Coal Mines, 2021, 52(8): 39−43.
[11] 郭明功,双海清,刘思博. 基于散斑监测的抽采钻孔变形破坏特征试验研究[J]. 煤矿安全,2023,54(7):130−136. GUO Minggong, SHUANG Haiqing, LIU Sibo. Experimental study on deformation and failure characteristics of drainage boreholes based on speckle monitoring[J]. Safety in Coal Mines, 2023, 54(7): 130−136.
[12] NAULIA P S, ROY A, WATADA J, et al. A pixel-wise k-immediate neighbour-based image analysis approach for identifying rock pores and fractures from grayscale image samples[J]. Algorithms, 2023, 16(1): 42. doi: 10.3390/a16010042
[13] 邹飞,李海波,周青春,等. 基于数字图像灰度相关性的类岩石材料损伤分形特征研究[J]. 岩土力学,2012,33(3):731−738. doi: 10.3969/j.issn.1000-7598.2012.03.013 ZOU Fei, LI Haibo, ZHOU Qingchun, et al. Fractal features study of rock-like material damage based on gray correlation of digital images[J]. Rock and Soil Mechanics, 2012, 33(3): 731−738. doi: 10.3969/j.issn.1000-7598.2012.03.013
[14] 吕祥锋,朱晨利,宋义敏,等. 煤岩稳定性与变形局部化演化关联性实验研究[J]. 岩石力学与工程学报,2021,40(12):2466−2476. LYU Xiangfeng, ZHU Chenli, SONG Yimin, et al. Experimental study on the correlation between coal rock stability and deformation localization evolution[J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(12): 2466−2476.
[15] 张科,张凯. 裂隙砂岩变形破裂过程中应变场灰度及纹理特征分析[J]. 煤炭学报,2021,46(4):1253−1262. ZHANG Ke, ZHANG Kai. Gray and texture features of strain field for fractured sandstone during failure process[J]. Journal of China Coal Society, 2021, 46(4): 1253−1262.
[16] 马德,赵育云,郭德龙,等. 孔周煤岩体渐进性破坏过程裂隙和应变能演化特征[J]. 煤矿安全,2022,53(3):16−23. MA De, ZHAO Yuyun, GUO Delong, et al. Evolution characteristics of cracks and strain energy during progressive failure of coal and rock masses around the hole[J]. Safety in Coal Mines, 2022, 53(3): 16−23.
[17] 肖杰,刘保国. 水泥石膏相似材料配比试验研究[J]. 岩土工程技术,2015,29(2):65−68,108. doi: 10.3969/j.issn.1007-2993.2015.02.003 XIAO Jie, LIU Baoguo. Research on cement & plaster similar materials proportioning test[J]. Geotechnical Engineering Technique, 2015, 29(2): 65−68,108. doi: 10.3969/j.issn.1007-2993.2015.02.003
[18] MARTIN C D. Seventeenth Canadian Geotechnical Colloquium: The effect of cohesion loss and stress path on brittle rock strength[J]. Canadian Geotechnical Journal, 1997, 34(5): 698−725. doi: 10.1139/t97-030
[19] 李小二,王鹏,靳翔飞,等. 基于CT扫描的煤样图像处理及裂隙形态表征[J]. 煤矿安全,2022,53(2):125−129. LI Xiaoer, WANG Peng, JIN Xiangfei, et al. Coal sample image processing and fracture morphological representation based on CT scanning[J]. Safety in Coal Mines, 2022, 53(2): 125−129.
[20] 王家臣,李良晖,杨胜利. 不同照度下煤矸图像灰度及纹理特征提取的实验研究[J]. 煤炭学报,2018,43(11):3051−3061. WANG Jiachen, LI Lianghui, YANG Shengli. Experimental study on gray and texture features extraction of coal and gangue image under different illuminance[J]. Journal of China Coal Society, 2018, 43(11): 3051−3061.
[21] WANG W X, LI L, HAN Y. Crack detection in shadowed images on gray level deviations in a moving window and distance deviations between connected components[J]. Construction and Building Materials, 2021, 271: 121885. doi: 10.1016/j.conbuildmat.2020.121885
[22] 王越,庞克亮,吴昊天,等. 基于机器学习的煤岩显微组分灰度及纹理特征提取实验[J]. 煤质技术,2023,38(2):56−68,84. doi: 10.3969/j.issn.1007-7677.2023.02.007 WANG Yue, PANG Keliang, WU Haotian, et al. Characteristics of gray and texture features of coal macerals based on machine learning[J]. Coal Quality Technology, 2023, 38(2): 56−68,84. doi: 10.3969/j.issn.1007-7677.2023.02.007