冲击荷载对无烟煤微观孔隙影响的分形研究

上一篇

下一篇

王以贤, 梁为民. 冲击荷载对无烟煤微观孔隙影响的分形研究[J]. 高压物理学报, 2020, 34(5): 054203-1. doi: 10.11858/gywlxb.20200528
引用本文: 王以贤, 梁为民. 冲击荷载对无烟煤微观孔隙影响的分形研究[J]. 高压物理学报, 2020, 34(5): 054203-1. doi: 10.11858/gywlxb.20200528
Yixian WANG, Weimin LIANG. Fractal Study on Influence of Impact Load on Microscopic Pore of Anthracite[J]. gaoyawli, 2020, 34(5): 054203-1. doi: 10.11858/gywlxb.20200528
Citation: Yixian WANG, Weimin LIANG. Fractal Study on Influence of Impact Load on Microscopic Pore of Anthracite[J]. gaoyawli, 2020, 34(5): 054203-1. doi: 10.11858/gywlxb.20200528

冲击荷载对无烟煤微观孔隙影响的分形研究

    作者简介: 王以贤(1981–),男,博士研究生,讲师,主要从事爆破工程及岩土工程研究.E-mail:56653648@qq.com .
    通讯作者: 梁为民(1967–),男,博士,教授,主要从事爆破工程及岩土工程研究.E-mail:liangwm@hpu.edu.cn
  • 中图分类号: O347.2

Fractal Study on Influence of Impact Load on Microscopic Pore of Anthracite

    Corresponding author: Weimin LIANG, liangwm@hpu.edu.cn
  • MSC: O347.2

  • 摘要: 为了揭示冲击荷载对无烟煤微观孔隙结构的影响规律,利用分离式霍普金森压杆(SHPB)冲击实验系统模拟了冲击应力在不同衰减过程中的冲击波和应力波,基于冲击前后压汞实验及低温液氮实验测试数据,运用分形理论,研究了赵固二矿不同方向煤体(与层理方向分别呈垂直、平行、45°斜交)冲击前后孔隙结构的分形特征。结果表明:对于渗流孔,冲击荷载提高了瓦斯的渗流与运移速度,对于吸附孔,冲击荷载减小了吸附孔的吸附能力,促进了瓦斯的解吸;分形维数具有明显的冲击方向性,且吸附孔的分形维数明显小于渗流孔;不同方向无烟煤的最佳冲击荷载不同,垂直和斜交层理方向的最佳冲击荷载为51.80 MPa,平行层理方向的最佳冲击荷载为28.46 MPa。研究结果可为冲击荷载促进瓦斯抽采机理的探讨提供参考。
  • 加载中
  • 图 1  煤样图片及尺寸

    Figure 1.  Picture and size of the coal sample

    图 2  实验方案

    Figure 2.  Experimental scheme

    图 3  垂直层理方向渗流孔分形维数的对比

    Figure 3.  Fractal dimension of the seepage hole in the vertical bedding direction

    图 4  平行层理方向渗流孔分形维数的对比

    Figure 4.  Fractal dimension of the seepage hole in the parallel bedding direction

    图 5  斜交层理方向渗流孔分形维数的对比

    Figure 5.  Fractal dimension of the seepage hole in the oblique bedding direction

    图 6  垂直层理方向吸附孔分形维数的对比

    Figure 6.  Fractal dimension of the adsorption hole in the vertical bedding direction

    图 7  平行层理方向吸附孔分形维数的对比

    Figure 7.  Fractal dimension of the adsorption hole in the parallel bedding direction

    图 8  斜交层理方向吸附孔分形维数的对比

    Figure 8.  Fractal dimension of the adsorption hole in the oblique bedding direction

    图 9  不同方向煤体渗流孔与吸附孔分形维数减小率随冲击荷载的变化

    Figure 9.  Fractal dimension reduction rates of the seepage hole and adsorption hole vary with the impact load

    图 10  不同方向煤体的分形维数与冲击荷载之间的关系

    Figure 10.  Relationship between the fractal dimension and the impact load of coal bodies in different directions

    表 1  煤样基本参数

    Table 1.  Basic parameters of coal samples

    No.$\;\rho $/(g·cm−3)fc/MPa CL/(km·s−1)R0,max/%Mass fraction/%
    AshVolatile
    component
    Fixed carbonExiniteVitriniteInertinite
    C01.42317.291.3813.325.626.0583.981981
    P01.46115.131.852
    X01.42211.651.653
    下载: 导出CSV

    表 2  实验数据表

    Table 2.  Impact test data

    DirectionpI/MPapL/MPaNo.Size/mm × mm$ D_{\rm{s} }$$\bar D_{\rm{s} }$$R_{\rm{s} }^2$$ D_{\rm{x} }$$\bar D_{\rm{x} }$$R_{\rm{x}}^2$$d $$\bar d$
    Vertical
    bedding
    00YC12.992.980.972.752.780.95
    YC22.990.992.780.96
    YC32.950.982.810.98
    0.1028.46C1 $ \varnothing $49.75 × 50.332.852.850.952.312.230.950.140.14
    C2$ \varnothing $49.33 × 49.892.920.962.210.930.16
    C3$ \varnothing $49.34 × 50.422.790.952.170.960.12
    0.1532.59C4$ \varnothing $49.42 × 49.572.732.730.992.192.190.950.180.16
    C5$ \varnothing $49.66 × 49.932.780.972.230.960.12
    C6$ \varnothing $49.42 × 50.842.690.942.150.950.18
    0.2041.43C7$ \varnothing $49.39 × 49.932.632.682.052.130.960.180.17
    C8$ \varnothing $49.42 × 49.762.710.982.210.970.17
    C9$ \varnothing $49.45 × 50.282.690.992.130.980.15
    0.3051.80C10$ \varnothing $49.52 × 50.472.562.510.992.052.050.990.170.19
    C11$ \varnothing $49.36 × 49.782.640.962.090.960.19
    C12$ \varnothing $49.27 × 50.272.340.942.010.960.22
    0.5058.70C13$ \varnothing $49.38 × 50.472.742.730.922.112.120.920.330.44
    C14$ \varnothing $49.45 × 50.212.790.942.170.930.57
    C15$ \varnothing $49.49 × 49.442.670.952.070.980.42
    Parallel
    bedding
    00YP12.842.870.992.892.860.95
    YP22.910.962.860.94
    YP32.870.972.830.96
    0.1028.46P1$ \varnothing $49.29 × 50.112.572.570.992.032.030.980.130.13
    P2$ \varnothing $49.36 × 50.372.650.982.050.950.11
    P3$ \varnothing $49.47 × 49.882.500.992.020.960.15
    0.1532.59P4$ \varnothing $49.31 × 49.322.582.580.972.152.151.000.310.25
    P5$ \varnothing $49.33 × 50.282.610.982.190.950.22
    P6$ \varnothing $49.31 × 49.162.560.952.110.970.22
    Parallel
    bedding
    0.2041.43P7$ \varnothing $49.33 × 50.332.642.640.922.202.190.950.180.28
    P8$ \varnothing $49.60 × 49.852.680.932.140.960.22
    P9$ \varnothing $49.31 × 50.112.610.952.240.980.44
    0.3051.80P10$ \varnothing $49.31 × 50.402.812.810.942.212.180.950.370.34
    P11$ \varnothing $49.05 × 49.872.850.952.160.970.30
    P12$ \varnothing $49.42 × 49.872.780.942.180.960.35
    0.5058.70P13$ \varnothing $49.40 × 49.932.892.810.952.822.730.970.540.50
    P14$ \varnothing $49.62 × 49.962.810.992.730.960.49
    P15$ \varnothing $49.32 × 49.792.730.952.630.960.47
    45° oblique bedding00YX12.812.810.942.772.780.95
    YX22.850.952.710.95
    YX32.780.972.850.96
    0.1028.46X1$ \varnothing $49.49 × 50.032.752.750.992.382.380.970.170.16
    X2$ \varnothing $49.27 × 49.632.810.982.280.950.19
    X3$ \varnothing $49.31 × 49.302.680.972.470.960.12
    0.1532.59X4$ \varnothing $49.24 × 49.932.672.650.992.392.310.980.290.24
    X5$ \varnothing $49.27 × 49.332.750.952.310.960.21
    X6$ \varnothing $49.27 × 49.822.540.982.230.950.22
    0.2041.43X7$ \varnothing $49.25 × 49.902.762.700.952.082.120.940.200.22
    X8$ \varnothing $49.60 × 49.382.710.922.170.920.21
    X9$ \varnothing $49.29 × 49.052.640.932.120.930.25
    0.3051.80X10$ \varnothing $49.40 × 50.282.472.550.982.012.060.930.210.18
    X11$ \varnothing $49.38 × 50.882.610.962.070.920.17
    X12$ \varnothing $49.71 × 50.022.580.972.110.950.16
    0.5058.70X13$ \varnothing $49.66 × 50.442.792.730.982.262.280.960.630.71
    X14$ \varnothing $49.73 × 49.912.750.972.310.940.72
    X15$ \varnothing $49.46 × 50.052.640.962.280.930.78
    下载: 导出CSV
  • [1] 高为, 易同生, 金军, 等. 黔西地区煤样孔隙综合分形特征及对孔渗性的影响 [J]. 煤炭学报, 2017, 42(5): 1258–1265. GAO W, YI T S, JIN J, et al. Pore integrated fractal characteristics of coal sample in western Guizhou and its impact to porosity and permeability [J]. Journal of China Coal Society, 2017, 42(5): 1258–1265.
    [2] 秦勇, 邱爱慈, 张永民. 高聚能重复强脉冲波煤储层增渗新技术试验与探索 [J]. 煤炭科学技术, 2014, 42(6): 1–7. QIN Y, QIU A C, ZHANG Y M. Experiment and discovery on permeability improved technology of coal reservoir based on repeated strong pulse waves of high energy accumulation [J]. Coal Science and Technology, 2014, 42(6): 1–7.
    [3] WEI G Y, SHAN Z Y, ZHANG Z M. Research on hydraulic slotting technology controlling coal-gas outbursts [J]. Journal of Coal Science and Engineering, 2008, 14(1): 67–72. doi: 10.1007/s12404-008-0014-4
    [4] 姜永东, 李业, 崔悦震, 等. 声场作用下煤储层渗透性试验研究 [J]. 煤炭学报, 2017, 42(S1): 154–159. JIANG Y D, LI Y, CUI Y Z, et al. Experimental study on characteristics of coal reservoir permeability under acoustic wave [J]. Journal of China Coal Society, 2017, 42(S1): 154–159.
    [5] 赵丽娟, 秦勇. 超声波作用对改善煤储层渗透性的试验分析 [J]. 天然气地球科学, 2014, 25(5): 747–752. ZHAO L J, QIN Y. Experiment on improving the permeability of coal reservoir under ultrasound [J]. Natural Gas Geoscience, 2014, 25(5): 747–752.
    [6] FORRESTAL M J, FREW D J, HANCHAK S J, et al. Penetration of grout and concrete targets with ogive-nose steel projectiles [J]. International Journal of Impacting Engineering, 1996, 18(5): 465–476. doi: 10.1016/0734-743X(95)00048-F
    [7] 穆朝民, 王海露, 黄文尧, 等. 高瓦斯低透气性煤体定向聚能爆破增透机制 [J]. 岩土力学, 2013, 34(9): 2496–2500. MU C M, WANG H L, HUANG W Y, et al. Increasing permeability mechanism using directional cumulative blasting in coal seams with high concentration of gas and low permeability [J]. Rock and Siol Mechanics, 2013, 34(9): 2496–2500.
    [8] 李恒乐, 秦勇, 张永民, 等. 重复脉冲强冲击波对肥煤结构影响的实验研究 [J]. 煤炭学报, 2015, 40(4): 915–921. LI H L, QIN Y, ZHANG Y M, et al. Experimental study on the effect of strong repetitive pulse shockwave on the pore structure of fat coal [J]. Journal of China Coal Society, 2015, 40(4): 915–921.
    [9] CAI Y D, LIU D M, PAN Z J, et al. Pore structure and its impact on CH4 adsorption capacity and flow capability of bituminous and subbituminous coals from Northeast China [J]. Fuel, 2013, 103: 258–268. doi: 10.1016/j.fuel.2012.06.055
    [10] MANDELBROT B B. The fractal geometry of nature [M]. San Francisco: Freeman, 1982: 35.
    [11] BIRD N, DÍAZ M C, SAA A, et al. Fractal and multifractal analysis of pore-scale images of soil [J]. Journal of Hydrology, 2006, 322(1/2/3/4): 211–219.
    [12] DATHE A, EINS S, NIEMEYER J, et al. The surface fractal dimension of the soil-pore interface as measured by image analysis [J]. Geoderma, 2001, 103(1/2): 203–229.
    [13] HUANG S J, YU Y C, LEE T Y, et al. Correlations and characterization of porous solids by fractal dimension and porosity [J]. Physica A: Statistical Mechanics and Its Applications, 1999, 274(3/4): 419–432.
    [14] HILDGEN P, NEKKA F, HILDGEN F, et al. Macroporosity measurement by fractal analysis [J]. Physica A: Statistical Mechanics and Its Applications, 1997, 234(3/4): 593–603.
    [15] YANG F, NING Z F, LIU H Q. Fractal characteristics of shales from a shale gas reservoir in the Sichuan Basin, China [J]. Fuel, 2014, 115: 378–384. doi: 10.1016/j.fuel.2013.07.040
    [16] LIU X J, XIONG J, LIANG L X. Investigation of pore structure and fractal characteristics of organic-rich Yanchang formation shale in central China by nitrogen adsorption/desorption analysis [J]. Journal of Natural Gas Science and Engineering, 2015, 22(1): 62–72.
    [17] 朱汉卿, 贾爱林, 位云生, 等. 蜀南地区五峰-龙马溪组页岩微观孔隙结构及分形特 [J]. 科学技术与工程, 2018, 18(10): 12–19. doi: 10.3969/j.issn.1671-1815.2018.10.003 ZHU H Q, JIA A L, WEI Y S, et al. Microscopic pore structure and fractal characteristics of Wufeng-Longmaxi shale, south Sichuan [J]. Science Technology and Engineering, 2018, 18(10): 12–19. doi: 10.3969/j.issn.1671-1815.2018.10.003
    [18] GARBACZ J K. Fractal description of partially mobile single gas adsorption on energetically homogeneous solid adsorbent [J]. Colloids and Surfaces A: Physicochemical and Engineering Aspects, 1998, 143(1): 95–101.
    [19] LEE G J, PYUN S Ⅰ, RHEE C K. Characterization of geometric and structural properties of pore surfaces of reactivated microporous carbons based upon image analysis and gas adsorption [J]. Microporous and Mesoporous Materials, 2006, 93(1/2/3): 217–225.
    [20] YAO Y B, LIU D M, TANG D Z, et al. Fractal characterization of adsorption-pores of coals from North China: an investigation on CH4 adsorption capacity of coals [J]. International Journal of Coal Geology, 2008, 73(1): 27–42. doi: 10.1016/j.coal.2007.07.003
    [21] XU L J, ZHANG D J, XIAN X F. Fractal dimensions of coals and cokes [J]. Journal of Colloid and Interface Science, 1997, 190(2): 357–359. doi: 10.1006/jcis.1997.4885
    [22] 宋昱, 姜波, 李凤丽, 等. 低-中煤级构造煤纳米孔分形模型适用性及分形特征 [J]. 地球科学, 2018, 43(5): 1611–1622. SONG Y, JIANG B, LI F L, et al. Applicability of fractal models and nanopores’ fractal characteristics for low-middle rank tectonic deformed coals [J]. Earth Science, 2018, 43(5): 1611–1622.
    [23] 王登科, 刘淑敏, 魏建平, 等. 冲击载荷作用下煤的破坏特性试验研究 [J]. 采矿与安全工程学报, 2017, 34(3): 594–600. WANG D K, LIU S M, WEI J P, et al. The failure characteristics of coal under impact load in laboratory [J]. Journal of Mining & Safety Engineering, 2017, 34(3): 594–600.
    [24] 刘晓辉, 张茹, 刘建锋. 不同应变率下煤岩冲击动力试验研究 [J]. 煤炭学报, 2012, 37(9): 1528–1534. LIU X H, ZHANG R, LIU J F. Dynamic test study of coal rock under different strain rates [J]. Journal of China Coal Society, 2012, 37(9): 1528–1534.
    [25] XODOT B B. 煤与瓦斯突出 [M]. 宋世钊, 王佑安, 译. 北京: 中国工业出版社, 1966: 310–318.
    [26] YUAN J, TAKEDA N, WAAS A M. A note on data processing in the split Hopkinson pressure bar tests [J]. Experimental Techniques, 1998, 22(5): 21–24. doi: 10.1111/j.1747-1567.1998.tb02317.x
    [27] SHAN R L, JIANG Y S, LI B Q. Obtaining dynamic complete stress-strain curves for rock using the split Hopkinson pressure bar technique [J]. International Journal of Rock Mechanics and Mining Sciences, 2000, 37(6): 983–992. doi: 10.1016/S1365-1609(00)00031-9
    [28] 陈俊宇, 裴向军, 杜瑞锋, 等. 冲击载荷作用下砂岩的动力学特性及能耗规律 [J]. 科学技术与工程, 2019, 19(31): 304–310. doi: 10.3969/j.issn.1671-1815.2019.31.045 CHEN J Y, PEI X J, DU R F, et al. Dynamic characteristics and energy consumption of sandstone under impact loading [J]. Science Technology and Engineering, 2019, 19(31): 304–310. doi: 10.3969/j.issn.1671-1815.2019.31.045
    [29] KOLSKY H. An investigation of the mechanical properties of materials at very high rates of loading [J]. Proceedings of the Physical Society (Section B), 1949, 62: 676–700. doi: 10.1088/0370-1301/62/11/302
    [30] LUNDBERG B. A split Hopkinson bar study of energy absorption in dynamic rock fragmentation [J]. International Journal of Rock Mechanics and Mining Science & Geomechanics Abstracts, 1976, 13(6): 187–197.
    [31] 赵迪斐, 郭英海, WANG G, 等. 基于分形建模的高煤级煤孔隙结构特征量化表征-以阳泉矿区山西组煤样为例 [J]. 东北石油大学学报, 2019, 43(3): 53–67. doi: 10.3969/j.issn.2095-4107.2019.03.006 ZHAO D F, GUO Y H, WANG G, et al. Quantitative characterization of pore structure characteristics of high quality coal based on fractal modeling: taking coal samples from Shanxi Formation in Yangquan Mining Area as an example [J]. Journal of Northeast Petroleum University, 2019, 43(3): 53–67. doi: 10.3969/j.issn.2095-4107.2019.03.006
    [32] RIGBY S P. Predicting surface diffusivities of molecules from equilibrium adsorption isotherms [J]. Colloids and Surfaces A: Physicochemical and Engineering Aspects, 2005, 262(1/2/3): 139–149.
    [33] 杨峰, 宁正福, 张世栋, 等. 基于氮气吸附实验的页岩孔隙结构表征 [J]. 天然气工业, 2013, 33(4): 135–140. doi: 10.3787/j.issn.1000-0976.2013.04.025 YANG F, NING Z F, ZHANG S D, et al. Characterization of pore structures in shales through nitrogen adsorption experiment [J]. Natural Gas Industry, 2013, 33(4): 135–140. doi: 10.3787/j.issn.1000-0976.2013.04.025
    [34] 吉小峰. 煤中纳米孔隙发育特征及其对气体运移的控制机理研究 [D]. 焦作: 河南理工大学, 2018. JI X F. Development characteristics of nanopores in coal and its controlling mechanism on gas migration [D]. Jiaozuo: Henan Polytechnic University, 2018.
    [35] YAO Y B, LIU D M, TANG D Z, et al. Fractal characterization of seepage-pores of coals from China: an investigation on permeability of coals [J]. Computers & Geosciences, 2009, 35(6): 1159–1166.
    [36] 宋晓夏, 唐跃刚, 李伟, 等. 中梁山南矿构造煤吸附孔分形特征 [J]. 煤炭学报, 2013, 38(1): 134–139. SONG X X, TANG Y G, LI W, et al. Fractal characteristics of adsorption pores of tectonic coal from Zhongliangshan southern coalmine [J]. Journal of China Coal Society, 2013, 38(1): 134–139.
    [37] 林海飞, 刘静波, 严敏, 等. CO2/CH4在煤储层中扩散规律的分子动力学模型 [J]. 中国安全生产科学技术, 2017, 13(1): 84–89. LIN H F, LIU J B, YAN M, et al. Molecular dynamics simulation on diffusion rules of CO2/CH4 in coal reservoir [J]. Journal of Safety Science and Technology, 2017, 13(1): 84–89.
    [38] 董怡静, 韩雨桢, 候泉林, 等. 煤变形产气的力化学机理探讨 [J]. 煤炭学报, 2017, 42(4): 942–949. DONG Y J, HAN Y Z, HOU Q L, et al. Mechanochemistry mechanism of gas generation during coal deformation [J]. Journal of China Coal Society, 2017, 42(4): 942–949.
    [39] 刘运通, 高文学. 爆炸荷载下岩石损伤的数值模拟研究 [J]. 岩石力学与工程学报, 2001, 20(6): 789–792. doi: 10.3321/j.issn:1000-6915.2001.06.007 LIU Y T, GAO W X. Numerical simulation on rock damage under explosion loading [J]. Chinese Journal of Rock Mechanics and Engineering, 2001, 20(6): 789–792. doi: 10.3321/j.issn:1000-6915.2001.06.007
    [40] 唐红梅, 周云涛, 廖云平. 地下工程施工爆破围岩损伤分区研究 [J]. 振动与冲击, 2015, 34(23): 202–206. TANG H M, ZHOU Y T, LIAO Y P. Damage zone of surrounding rock of underground engineering under construction blasting [J]. Journal of Vibration and Shock, 2015, 34(23): 202–206.
  • 加载中
图( 10) 表( 2)
计量
  • 文章访问数:  2661
  • HTML全文浏览数:  2661
  • PDF下载数:  11
  • 施引文献:  0
出版历程
  • 收稿日期:  2020-03-12
  • 刊出日期:  2020-10-05

冲击荷载对无烟煤微观孔隙影响的分形研究

    通讯作者: 梁为民(1967–),男,博士,教授,主要从事爆破工程及岩土工程研究.E-mail:liangwm@hpu.edu.cn
    作者简介: 王以贤(1981–),男,博士研究生,讲师,主要从事爆破工程及岩土工程研究.E-mail:56653648@qq.com
  • 1. 河南理工大学土木工程学院,河南 焦作 454000
  • 2. 河南建筑职业技术学院,河南 郑州 450064

摘要: 为了揭示冲击荷载对无烟煤微观孔隙结构的影响规律,利用分离式霍普金森压杆(SHPB)冲击实验系统模拟了冲击应力在不同衰减过程中的冲击波和应力波,基于冲击前后压汞实验及低温液氮实验测试数据,运用分形理论,研究了赵固二矿不同方向煤体(与层理方向分别呈垂直、平行、45°斜交)冲击前后孔隙结构的分形特征。结果表明:对于渗流孔,冲击荷载提高了瓦斯的渗流与运移速度,对于吸附孔,冲击荷载减小了吸附孔的吸附能力,促进了瓦斯的解吸;分形维数具有明显的冲击方向性,且吸附孔的分形维数明显小于渗流孔;不同方向无烟煤的最佳冲击荷载不同,垂直和斜交层理方向的最佳冲击荷载为51.80 MPa,平行层理方向的最佳冲击荷载为28.46 MPa。研究结果可为冲击荷载促进瓦斯抽采机理的探讨提供参考。

English Abstract

  • 煤体孔隙是煤层气(煤矿瓦斯)赋存、扩散、运移的主要场所[1]。通过外部扰动优化煤体孔径配置、改善煤体孔隙结构是促进煤层气解吸、增强煤层气扩散和运移的主要手段。目前,通过外部扰动改善煤体孔隙结构主要有3种方式:(1)流体冲载,主要包括水力、高能气体压裂等[2-3];(2)物理场激励,主要有静电场、交变电场、电磁场、声场激励等[4-5];(3)冲击荷载,如爆破[6-7]、脉冲冲击波[8]等。

    我国无烟煤煤层气储层孔隙主要发育在纳米级别[9],煤层气含量高,但吸附性强、渗透率低,制约着我国煤层气的产量。因此,要实现煤层气的有效开采,需要在了解煤体微观孔隙形态与分布的基础上,深入探讨外载促进煤体孔隙改变及增强瓦斯渗透性的机理。目前,对于煤体微观孔隙结构的研究大多是针对自然状态下的不同类型煤体,对于冲击荷载作用下煤体微观孔隙结构的研究还鲜见报道。

    分形几何学是由Mandelbrot等[10]于1982年提出和建立的,该理论能够为研究复杂的不规则事物的变化规律提供强有力的工具和途径。大量科研人员应用分形理论研究了土和岩石的孔隙分形[11-13],并得出了不同岩性的孔隙度、比表面积、吸附性及孔隙表面粗糙度、不均匀性等与分形维数的定量关系[14-17],而对于煤体孔隙,特别是外载扰动后煤体孔隙的分形研究相对较少。事实上煤体是具有基质孔隙和裂隙的多孔介质,其复杂的孔隙结构满足自相似性,符合分形几何规律[18-20]。因此,可以将分形参数作为煤体内部复杂孔隙结构的定量表征[21-22],通过对比煤体受载前后的分形参数,探讨煤体外载作用对其孔隙结构的影响规律。

    本研究利用霍普金森冲击实验系统(SHPB)模拟冲击应力在不同衰减过程中的冲击波和应力波,对赵固二矿无烟煤实施不同方向(分别与层理呈垂直、平行、45°斜交方向)、不同速度的冲击加载,利用压汞实验及低温液氮实验测试数据,对冲击前后无烟煤的微观孔隙进行分形研究,力图得出冲击荷载对无烟煤微观孔隙的破坏特征及规律,以期为冲击荷载促进瓦斯抽采规律的研究提供参考。

    • 实验样品取自河南省焦作矿区赵固二矿,属于无烟煤,其显微组分、工业分析结果及基本力学参数如表1所示,其中:C、P、X分别表示垂直、平行、45°斜交层理的煤样,下标0代表冲击前原煤样,$\;\rho$为样品密度,fc为单轴抗压强度,CL为纵波声速,R0, max为最大镜质组反射率。为了尽量减小煤样间的差异,在矿井采煤工作面同一煤层同一地点采集大块完整煤岩。按照GB/T50266—2013 《工程岩体实验方法标准》[23]及GB/T23561.7—2009 《煤和岩石物理力学性质测定方法》[24]钻取直径和高度均为50 mm的圆柱体,如图1所示。利用双端面磨床打磨试件,确保煤样两端面的平整度小于或等于0.05 mm,两端面平行度小于或等于0.02 mm。

    • 本实验方案如图2所示。实验步骤分为4部分。

    • (1)SHPB冲击实验。采用入射杆和透射杆均为变截面杆的SHPB冲击实验系统,杆端直径均为50 mm,入射杆长2 400 mm,透射杆长1 200 mm,杆件材质为钢材,密度为7 800 kg/m3,弹性波速为5 190 m/s,弹性模量为210 GPa。根据冲击波(应力波)随距离衰减规律,结合样品的基本力学参数及试冲情况,设定冲击气压pI分别为0.10、0.15、0.20、0.30和0.50 MPa,利用河南理工大学土木工程学院的SHPB冲击加载系统进行冲击加载实验。首先,调整子弹、入射杆、透射杆呈一条直线,将试件涂抹耦合剂以后装在入射杆与透射杆之间;其次,将尺寸为1 cm × 1 cm的铜整形片粘贴至入射杆靠近子弹一端,按照设定的冲击气压驱动子弹撞击入射杆;最后,采集信号,保存数据。

      (2)取样并干燥处理。为了保证冲击前后压汞实验和低温液氮实验数据的可对比性,实验样品均取自冲击前后圆柱体中轴线附近的煤颗粒,将实验样品放入烘干箱,在110 ℃条件下恒温干燥8 h以后放入干燥器皿中冷却至室温。

      (3)压汞实验。将3 g左右4~6 mm的煤颗粒放入膨胀计中,抽真空以后将膨胀计装入低压站进行低压实验,低压压汞实验结束后将膨胀计移入高压仓进行高压实验,然后结合压汞分析软件记录数据。压汞仪采用美国麦克仪器公司Autopre Ⅳ 9520型全自动压汞仪,工作压力介于0.1~60 000.0 psi(1 psi=0.006 895 MPa),测定孔径范围为3 nm~1 100 μm。

      (4)低温液氮实验。将3 g左右60~80目的煤粉脱气8 h后放入分析站进行低温液氮吸附实验。低温液氮仪为美国麦克仪器公司TriStar Ⅱ 3020型全自动比表面积与孔隙度分析仪,测定孔径范围为0.35~500 nm。

      本实验的孔径划分采用Hodot的十进制划分法[25],即微孔、小孔、中孔、大孔对应的孔径依次为小于10 nm、介于10~100 nm之间、介于100~1 000 nm之间、大于1 000 nm。

    • 根据一维弹性应力波和应力均匀的假设,采用三波法[26-27]计算煤样的动态应力、应变及应变率。子弹以一定速度冲击入射杆撞击试件,首先使试件产生应力、应变并对试件造成一定的损伤,然后形成透射波,最后由吸收系统吸收透射波。该过程中能量传递和转化的表达式为

      式中:${W}$${\sigma}$${\varepsilon}$分别为3种波对应的能量、应力和应变,$V$为试样体积,A0为压杆的横截面积,E为弹性模量,c为波速。

      损伤是外界能量作用于物体的一种客观存在,其量化指标为损伤变量,一般来讲,可以从面积、密度、弹性模量和能量等不同角度定义损伤变量。冲击能量对试件造成的损伤可以通过能量公式计算其损伤变量

      式中:$ d $为损伤变量;$u$为试件的总吸收能密度,即试件的冲击应力-应变曲线所围成的面积,表示为

      鉴于煤体的非均质性及各向异性,同等实验条件下冲击3个试件,取各个参数的平均值进行分析,其中损伤变量为总耗散能密度与总吸收能密度的比值。

      利用式(1)~式(4)可以计算出煤样受冲击后的损伤值[28-30],如表2所示,其中:pI为冲击气压;pL为冲击荷载;煤样尺寸为圆柱体的直径与高度;$D_{\rm{s}} $$\bar D_{\rm{s}}$分别为渗流孔分形维数和渗流孔分形维数平均值;$D_{\rm{x}} $$\bar D_{\rm{x}}$分别为吸附孔分形维数和吸附孔分形维数平均值;$R_{\rm{s}}^2$$R_{\rm{x}}^2$分别为渗流孔和吸附孔分形维数相关系数;$d $$\bar d$分别为损伤变量和平均损伤变量;YC、YP、YX分别代表冲击前垂直、平行和45°斜交层理方向的煤样,下标分别代表不同层理方向的煤样序号。原煤压汞实验和低温液氮实验中的煤样取自冲击前煤块颗粒,不存在煤样尺寸。

    • 高压压汞实验是测量固体孔隙参数的常规方法,孔径测量范围比较大,可以测出从微孔到大孔的所有孔隙。结合Washburn方程,利用Menger海绵分形理论[31]可以定量描述煤体孔隙的粗糙度和不规则程度,从而探析煤体孔隙的吸附、解吸及渗流参数。将Menger模型作为压汞实验数据的分形模型,结合压汞实验数据可以得出煤体的分形维数,计算公式为

      式中:pH为进汞压力,MPa;Vp为压力pH作用下的进汞量,cm3。以ln (dVp/dpH)和ln pH的比值为斜率K,则分形维数D=K+4。

    • 低温液氮实验是将吸附于一定固体表面上的气体量记录为吸附物质的吸附量与相对压力的函数,采用等温线脱附分支上一系列逐步降低的相对压力,或等温线吸附分支上一系列逐步升高的相对压力(也可二者并用)进行计算,从而得到多孔物质材料微观孔隙参数的实验方法。低温液氮数据可以采用FHH模型得出不同冲击荷载、不同方向煤体的分形维数[22]

      FHH模型的分形维数计算基于以下公式

      式中:$V$为平衡压力$p$下吸附的气体体积,cm3${V_0}$为单分子层吸附的气体体积,cm3$p$为平衡压力,MPa;${p_0}$为饱和蒸汽压,MPa;$A$为对数曲线的斜率;$n$为常数;分形维数$D = 3 + A$

    • 根据压汞实验及低温液氮实验的测试数据,结合式(1)、式(2)可以得出不同方向、不同冲击荷载作用下的分形维数(分形维数大小见表2)。通过计算发现:当孔径小于或等于100 nm时,基于压汞实验测试数据的分形维数不在2~3之间;而当孔径大于或等于100 nm时,基于液氮实验测试数据的分形维数也不在2~3之间。根据相关研究[32],以上数据不具有分形意义。这可能与两种方法的测试原理有关,压汞法是根据汞的不浸润性进行孔隙参数测试,进入孔隙直径越小,所需要汞的压力越大,而较大的汞压力会破坏孔隙结构,因此,压汞法对孔径小于或等于100 nm的孔隙参数测量不准确。而低温液氮的测试原理是吸附与凝聚理论,根据孔吸附液氮的体积计算孔隙参数,该方法对孔径大于或等于100 nm的孔隙参数测量不准确[33-34]。所以表2只统计了基于压汞实验数据、孔径大于100 nm孔隙的分形维数和基于液氮实验数据、孔径小于100 nm孔隙的分形维数(为了保证数据的有效统计意义,降低煤体离散性对结果的影响,表2中所有数据均取同等条件下4个煤样数据的平均值)。根据Yao等[20, 35]、宋晓夏等[36]的研究结果,将孔径大于100 nm的孔隙称为渗流孔,孔径小于100 nm的孔隙称为吸附孔。

      不同方向典型分形维数拟合曲线的对比如图3图8所示,相关系数均在0.90以上,拟合效果较好,数据可靠,符合分形规律。

      从分形统计数据可知:垂直层理方向渗流孔的分形维数介于2.51~2.98之间,平均值为2.75,最小值2.51对应的冲击荷载为51.80 MPa,冲击前分形维数最大;平行层理方向渗流孔的分形维数介于2.57~2.87之间,平均值为2.71,最小值对应的冲击荷载为28.46 MPa,冲击前分形维数最大;斜交层理方向渗流孔的分形维数介于2.55~2.81之间,平均值为2.70,最小值对应的冲击荷载为51.80 MPa,冲击前分形维数最大。这说明冲击前无烟煤煤样渗流孔的分形维数大,孔隙结构复杂,孔隙内部表面粗糙,孔壁摩阻力大[37],不利于瓦斯的渗流与运移;冲击后分形维数均小于冲击前,说明冲击荷载改善了孔隙内部结构,降低了孔隙表面粗糙度,减小了孔壁摩阻力,提高了瓦斯的渗流与运移速度。

      吸附孔分形维数存在以下特点:垂直层理方向吸附孔分形维数介于2.05~2.78之间,平均值为2.25,最小值2.05对应的冲击荷载为51.80 MPa,冲击荷载作用前分形维数最大;平行层理方向吸附孔分形维数介于2.03~2.86之间,平均值为2.36,冲击荷载为28.46 MPa时分形维数最小,冲击前分形维数最大;斜交层理方向吸附孔分形维数介于2.06~2.78之间,平均值为2.32,最小值2.06对应的冲击荷载为51.80 MPa,冲击前分形维数最大。这与渗流孔分形规律相似,说明冲击后吸附孔的分形维数均小于冲击前,冲击荷载减小了吸附孔的复杂程度、不规则程度及孔隙表面粗糙度,大大降低了吸附孔的吸附能力,冲击荷载能够将大量的瓦斯由吸附态转化为游离态,从而促进瓦斯的有效解吸。

    • 定义分形维数减小率δ

      式中:${D_1}$为冲击后分形维数,${D_0}$为冲击前分形维数。

      分形维数减小率可以表征冲击荷载对孔隙的改造程度,其值越大说明冲击荷载对该孔隙的改造程度越大,改造效果越好。图9表征的是3个方向渗流孔和吸附孔分形维数减小率与冲击荷载的关系曲线。从图9可以看出,对于渗流孔,在冲击后分形维数最小值相近的前提下(垂直、平行及斜交层理方向的分形维数的最小值分别为2.51、2.57和2.55),垂直和斜交层理方向的分形维数减小率随冲击荷载变化呈波动增加状态,平行层理方向的分形维数减小率随冲击荷载增加呈减小态势,吸附孔的分形维数减小率具有相同的变化规律。这说明冲击荷载改造煤体孔隙具有明显的冲击方向性,对于垂直和斜交层理方向煤体,需要较大冲击荷载才能达到较好的促解增渗效果,而对于平行层理方向煤体,较小冲击荷载就能达到很好的效果。因此,在实际抽采瓦斯过程中,应尽量沿平行层理方向布置炮孔,以达到消耗较少能量就能获得较好促解增渗的目的。

      进一步对比吸附孔与渗流孔的分形维数减小率发现,同一方向、同一冲击荷载下吸附孔的分形维数减小率始终大于渗流孔的分形维数减小率(见图9),说明冲击荷载对吸附孔的改造程度较渗流孔大,即相对于冲击荷载提高瓦斯渗流与运移速度来说,冲击荷载能够更大程度地促进瓦斯的解吸,从而促进大量吸附瓦斯转化为游离瓦斯。这也解释了外载扰动过程中会产生超量煤层气[38]、更容易发生瓦斯爆炸与瓦斯突出的原因。

      图10为不同方向煤体分形维数与冲击荷载关系曲线。从图10可以看出,冲击后渗流孔与吸附孔的分形维数并不随着冲击荷载的增加持续减小,而是存在一个最小值,且不同方向的最小值所对应的冲击荷载不同,但两种孔隙同一方向分形维数最小值所对应的冲击荷载却完全一致,即垂直与斜交层理方向,渗流孔与吸附孔分形维数最小值对应的冲击荷载均为51.80 MPa,平行层理方向分形维数最小值对应的冲击荷载均为28.46 MPa。这说明无论是渗流孔还是吸附孔,同一方向的最佳冲击荷载是相同的,且对应的损伤值均介于0.10~0.20之间(见表2),根据刘运通等[39]、唐红梅等[40]的研究,该损伤值范围正好属于爆破裂隙区,即不管哪个冲击方向,爆破裂隙区是最佳的致裂增渗区。

    • 采用压汞实验及低温液氮实验测试方法,研究了冲击荷载对无烟煤微观孔隙分形参数的影响规律,对比分析了冲击前后无烟煤的分形特征,得出以下主要结论。

      (1)孔径小于或等于100 nm的压汞实验数据及孔径大于或等于100 nm的低温液氮实验数据不具有分形规律,获得了孔径大于100 nm(渗流孔)、基于压汞实验数据的孔隙分形特征和孔径小于100 nm(吸附孔)、基于低温液氮实验数据的孔隙分形特征,且冲击荷载作用后吸附孔的分形维数明显小于渗流孔,说明冲击荷载改善吸附孔的能力较渗流孔强。

      (2)冲击荷载作用后的分形维数明显小于冲击前,说明冲击荷载改善了孔隙内部结构,降低了孔隙的复杂程度、不规则程度及表面粗糙度,减小了孔壁摩阻力。对于渗流孔,冲击荷载提高了瓦斯的渗流与运移速度;对于吸附孔,冲击荷载促进了瓦斯的解吸。

      (3)分形维数减小率能够反映冲击荷载对煤体孔隙的改造程度,不管是渗流孔还是吸附孔,垂直和斜交层理方向的分形维数变化率随冲击荷载变化呈波动增加状态,平行层理方向的分形维数变化率随冲击荷载增加呈减小态势,说明冲击荷载改造煤体孔隙具有明显的冲击方向性,平行层理方向在能量消耗较小的前提下就能达到较好的促解增渗效果。

      (4)不同方向煤体的渗流孔和吸附孔所对应的最佳冲击荷载相同,垂直和斜交层理方向煤体的最佳荷载均为51.80 MPa,平行层理方向煤体的最佳冲击荷载均为28.46 MPa,且对应的损伤值均属于爆破裂隙区。

    参考文献 (40)

目录

/

返回文章
返回