0.   引　　言

原生结构煤受构造应力破坏后，形成不同破坏程度的煤体结构，不同煤体结构煤的孔隙结构有明显差异性。孔隙内表面是吸附气存在的主要空间，研究不同煤体结构煤的孔隙结构、形态和连通性特征，有助于了解煤破坏变形后储层气体的吸附/解吸、扩散等运移规律的差异机理。目前对于多孔介质结构的研究方法主要有射线法和流体注入数据分析法。其中射线法包含光学显微镜、扫描电镜、透射电镜、核磁共振、CT扫描和小角X射线法；流体注入分析法目前常用的有压汞法、低温液氮吸附分析、低压二氧化碳吸附分析等^[1-2]。流体注入法中常用的是低温液氮和低压二氧化碳分析，利用高压和低压下气体在不同孔径内充填和吸附方式的差异而表现出来的不同形态的吸附−脱附曲线，从而推测煤中孔隙分布特征，并根据压力和吸附数据结合数值模型计算出多孔介质孔径、比表面积和孔隙容积等特征^[3-4]。低温液氮吸附分析测试孔径变化范围为1.7～300 nm^[5]，而由于CO₂分子直径相对较小，常用来表征小于2 nm孔隙的孔径分布、比表面积和孔隙容积特性^[6]。图像观测法和流体注入分析法常被用来研究煤中孔隙结构特征，与原生结构煤相比，经过破坏变形后的煤孔隙结构各向异性增强，构造变形会促进煤中大孔隙向小孔隙的尺度转化。例如，王佑安等^[7]通过压汞法对比了不同井田不同演化程度的构造变形煤样，发现煤破坏变形过程中微孔并未受到改变，中孔和过渡孔显著增加；但是JU等^[8]认为构造应力会破坏直径小于2 nm的微孔。随着煤破坏变形的增强，2～10 nm和10～20 nm孔径所占的比例随着破坏强度的增加而增加^[9]。王涛等^[10]研究液氮吸附分析对构造煤孔结构影响时发现5 nm以下的孔受影响较弱，5～10 nm孔隙显著增加。孟召平等^[11]研究发现随着煤体结构破坏程度的增高，孔容和比表面积也相应增大，孔容主要由中孔贡献，比表面积主要由微孔贡献。分形维数可精细化定量表征煤孔隙结构的复杂性行为^[12]，部分学者通过建立煤孔隙分形维数与变质程度、煤岩煤质和吸附能力等之间的定量关系以解释不同煤储层的物性特征及其演化机理^[13-15]。降文萍等^[16]研究发现孔隙结构分形维数随着煤破碎程度的增强而增大，解释了煤中孔隙系统非均质性强、连通性差的原因。金毅等^[17]基于压汞法对比了煤样受压前后基质孔隙的非均质性变化及分形特征。

目前对孔径分类方法有2种：一种是基于IUPAC提出的分类方法，它将孔隙分为微孔（<2 nm）、中孔（2～50 nm）和大孔（>50 nm）^[18]。HODOT在1966年提出了另一种广泛使用的方法，将孔径分为微孔（<10 nm）、过渡孔（10～100 nm）、中孔（100～1000 nm）和大孔（>1000 nm）^[19]。对于HODOT分类，一些学者研究发现煤中小于10 nm的微孔的孔径分布差异性较大，研究表明，小于2 nm的孔隙对于比表面积和孔容的贡献是2～50 nm的数百倍^[20]，说明煤中大量微孔分布于小于2 nm范围，将小于2 nm的孔隙划分到小于10 nm的微孔范围是不合理的，需要单独分类。小于2 nm的孔隙是煤层气主要的吸附孔，决定着甲烷的吸附量。IUPAC虽然将小于2 nm的孔隙单独划分为微孔，但是与原生结构煤相比，对构造煤的中孔（IUPAC）孔容贡献主要是来自于2～10 nm的孔隙，构造变形会将更多的10～50 nm的孔隙变形破坏到10 nm以下^[9]。因此，研究煤岩变形过程中孔隙在2 nm以下和2～10 nm的变化就显得非常重要。由于构造煤小于2 nm孔隙结构测量困难，现有可靠的测量数据极少。对于不同煤体结构煤孔隙结构的演化规律及分形特征的研究也相对较少。因此，基于低压二氧化碳和低温液氮吸附数据，分析不同煤体结构煤孔隙结构的演化及其分形特征，揭示了煤破坏变形后孔隙结构的分形规律，进一步探讨了分形维数与煤中甲烷吸附-扩散的关系，为构造煤区瓦斯防治和煤层气高效开发提供理论支撑。

1.   试验方法和计算模型

1.1   样品及预处理

不同煤体结构煤样品来自于沁水盆地赵庄井田，目标煤层为二叠系下统山西组3号煤层。煤岩类型以半亮煤、光亮煤为主，不同煤体结构样品断口明显，断裂面较为平整且无拉伸或压扭现象，样品仍保持了相当的硬度，为脆性变形煤。煤体结构包含原生结构、碎裂结构、碎粒结构和糜棱结构。随着煤破坏程度增强，煤中层理结构消失，破碎煤岩颗粒和变形孔−裂隙增多（图1）。煤样的基础数据见表1，不同煤体结构之间煤岩煤质及变质程度的差异性较小。在对不同煤体结构样品进行宏观和微观观察后，将样品粉碎并研磨至40～60目（0.250～0.425 mm），测试样品质量在2～3 g。在液氮沸点温度（101.3 kPa，77.35 K）下，以低温液氮为吸附质，用ASAP 2020孔隙率分析仪测定大于2 nm的孔隙结构特征，如比表面积、孔容和孔径分布。进行吸附分析前，将样品在110 ℃真空下脱气24 h，以除去样品中吸附的挥发性物质。随后将样品放入分析仪并设定测试压力，仪器软件将根据多个压力点下的吸附−解吸量，通过经典的计算模型，如BET和BJH计算比表面积和孔隙容积。由于N₂的动力学直径较大，以及与孔壁的四极矩作用，N₂很难进入小于2 nm的微孔^[21]，氮气的测试孔径通常分布在1.7～300 nm。而CO₂具有较小的分子动力学直径（0.33 nm）和较高的饱和蒸气压（237 K时为3.5 MPa），可用于计算小于2 nm的超微孔孔径。低压二氧化碳测试试验样品和制备过程与氮气吸附试验相似，但孔径分布和孔隙体积的计算参数模型是通过DFT理论模型计算，该模型考虑了表面粗糙度和各向异性对超微孔的影响。另取样品研磨至60～80目（0.178～0.250 mm），质量100 g左右，通过容积法进行甲烷的等温吸附试验，试验温度为35 ℃，压力范围为1～8 MPa。

图  1  不同煤体结构宏观和微观煤岩

Figure  1.  Macroscopic and microscopic coal petrography of different coal structures

下载: 全尺寸图片 幻灯片

表  1  样品工业分析及显微组分

Table  1.  Proximate analysis of samples and macerals of coal structures

样品	Ro/%	工业分析/%			镜质组体积分数/%	惰质组体积分数/%
		Mad	Ad	Vdaf
原生结构	2.23	0.98	14.22	11.71	65.80	33.77
碎裂结构	2.50	1.07	13.14	10.57	65.93	34.07
碎粒结构	2.48	0.90	12.14	9.17	72.07	27.93
糜棱结构	2.43	0.96	12.47	11.87	63.84	35.91

下载: 导出CSV 

| 显示表格

1.2   分形模型

分形行为可定量表征多孔介质孔隙结构的复杂程度，它的实质是吸附剂通过吸附不同分子直径的气体，再结合数值模型测量不同尺度下多孔介质吸附量，然后依据吸附量与压力之间的拟合关系求取分形参数以定量表征孔隙结构的不规则度。CO₂和N₂由于分子直径的差异，可以进入不同孔隙空间，通常依据CO₂吸附行为表征小于2 nm的孔隙结构分形特征，N₂则表征大于2 nm的孔隙分形特征。WANG等^[22]通过分析CO₂吸附量和相对压力（μ=P/P₀）之间的关系，提出超微孔的表面分形规律可通过下式表征：

式中：D_m为超微孔表面分形维数；C为常数项；$A\left(\mu \right)=\dfrac{-\displaystyle\int _{N\left(\mu \right)}^{{N}_{{\rm{max}}}}\mathrm{ln}\left(\mu \right){\rm{d}}N\left(\mu \right)}{r{\left(\mu \right)}^{2}}$，$B(\mu)=\dfrac{{\left({N}_{{\rm{max}}}-{N}_{\mu }\right)}^{\frac{1}{3}}}{r\left(\mu \right)}$，其中N_max为最大相对压力下吸附量，cm³/g；N(μ)为相对压力μ时的吸附量，cm³/g；不同压力阶段CO₂覆盖的孔隙空间不同，测得的平均孔径r与相对压力μ之间函数关系可用开尔文公式表征：

式中：σ为表面张力；v为摩尔体积，22.4 L/mol；R为理想气体常数，8.314 J/(mol·K)；T为绝对温度，K。

依据低温液氮吸附数据计算的分形维数被证明是表征多孔介质中孔隙几何结构的有效方法，通过定量定义分形维数D₁和D₂来表征孔隙表面粗糙度和孔隙结构各向异性特征^[14]。根据Frenkel-Halsey-Hill（FHH）模型，基于液氮吸附数据，可通过以下方法计算分形维数D：

式中：V/V₀为平衡压力P下的相对吸附量；C为lnV/V₀随ln（ln（P₀/P））函数变化关系投影图中y的截距；A为用于计算分形维数D的幂律指数；P₀为饱和压力，MPa。通过拟合数据的斜率，通常有2种方法用于测量煤中孔隙结构分形维数D：

由于式（5）计算的分形维数D值一般小于2，这偏离了煤孔隙结构的定义^[23]。本次采用方程（4）来表征大于2 nm的孔隙的分形维数，该值通常介于2～3。

1.3   扩散模型

通过不同煤体结构煤样等温吸附试验数据中甲烷解吸动力学数据随时间和空间变化关系，结合非稳态菲克第二扩散定律（式（6）），从而定性和定量表征不同煤体结构的扩散性能：

式中：r为扩散半径，mm；C为吸附质浓度，g/cm³；D为扩散系数，m²/s；t为时间，s。理想条件下数值模型式（6）的解析解可表示为

式中：M_t为在时间t内扩散气体的总量，cm³；M_∞为无限时间内的总解吸气体量，cm³；r为扩散路径长度，mm；D为扩散系数，m²/s。其中，扩散路径（r）通常表征煤粒半径，但是根据等温吸附试验规范，煤样经过破碎后，其半径难以估算，通常取其等效半径r=0.1075 mm。但是不同煤体结构的煤粒半径分布不均，用等效半径计算存在误差。因此，此次论文中用有效扩散系数表征煤体结构扩散性能，即D_e=D/r², s⁻¹，忽略等效半径的计算影响，提高表征准确性。

2.   试验结果和分析

2.1   液氮和二氧化碳测试结果

通过N₂和CO₂吸附分析可测量煤中不同尺度孔隙结构分布特征。不同煤体结构的液氮吸附-脱附曲线显示出相似的形态（图2a），对应于IUPAC多孔固体分类的IV型等温线和H₃型滞后环^[24]。图中明显的滞后回线说明在低压区（P/P₀<0.2）、中压区（P/P₀=0.2～0.9）和高压区（P/P₀>0.9），孔隙中对于N₂的吸附与脱附过程是不同的。图2显示N₂的吸附量随煤破坏程度的增强而增大，糜棱结构煤对N₂的吸附量最大（1.93 cm³/g），其次是碎粒结构（1.76 cm³/g）、碎裂结构（1.53 cm³/g）和原生结构（0.85 cm³/g）。煤中气体主要吸附于孔隙内表面，经过破坏变形后的煤体结构比表面积增大（表2），吸附能力增强。

图  2  不同煤体结构N₂和CO₂吸附−脱附曲线

Figure  2.  Adsorption-desorption curves of N₂ and CO₂ analysis

下载: 全尺寸图片 幻灯片

表  2  液氮、二氧化碳及等温吸附实验测试结果

Table  2.  Results of low temperature nitrogen, pressure carbon dioxide and isothermal adsorption experiment

样品	平均孔径/ nm	BET比表面积 / (m2·g−1)	BJH孔容/ (cm3·g−1)	CO2测试 孔径/nm	DFT比表面积/ (m2·g−1)	DFT孔容/ (cm3·g−1)	Langmuir吸附量/ (cm3·g−1)	Langmuir压力/ MPa
原生结构	24.671	0.2141	0.0011	0.524～1.083	58.086	0.0196	32.26	1.64
碎裂结构	14.972	0.3251	0.0013	0.524～1.083	61.838	0.0208	34.48	1.62
碎粒结构	21.049	0.5683	0.0024	0.507～1.101	78.433	0.0267	34.72	1.56
糜棱结构	10.791	1.5716	0.0030	0.489～1.083	105.808	0.0340	36.50	1.65

下载: 导出CSV 

| 显示表格

由于N₂和CO₂吸附试验在测试温度、分子动能、分子直径和测试机理等方面存在差异性，因此CO₂可以比N₂进入更小的孔隙（<2 nm）空间。图2b显示4种煤体结构CO₂吸附等温线属于IUPAC分类的Ⅰ型^[24]。吸附等温线均随相对压力的增大而增大，且对于CO₂的吸附量整体上随着煤体结构破坏程度的增强而增大。不同煤体结构的超微孔比表面积和孔容均随着煤破坏程度的增强而逐渐增大，变化范围58.086～105.808 m²/g和为0.0196～0.0340 cm³/g（表2）。超微孔在构造煤中的分布差异性较大，特别是糜棱结构煤超微孔贡献的容积和比表面积几乎是原生结构的2倍。表2中糜棱结构煤的测试超微孔孔径范围在0.489～1.083 nm，大于原生结构煤0.524～1.083 nm的范围，破坏变形越强烈的煤超微孔孔径分布要相对较小。

2.2   不同煤体结构孔隙结构分布

2.2.1   微孔、中孔和大孔分布

微孔（2～10 nm）、中孔（10～100 nm）和大孔（>100 nm）对于孔容和比表面积的贡献率在不同煤体结构之间存在着显著差异，如图3a和图3b所示。随着煤体结构破坏程度增强，大孔贡献的孔容和比表面积比例降低，微孔和中孔（10～50 nm）比例增加，说明大孔和部分中孔（50～100 nm）被破坏成中孔（10～50 nm）和微孔。不同煤体结构的总BET比表面积主要由微孔（2～10 nm）和中孔（10～100 nm）贡献，平均百分比为88.96%，尤其糜棱结构煤在2～10 nm范围内增长率最大，说明构造变形严重影响了2～10 nm的孔隙分布。图4为阶段比表面积随孔径变化，在10～300 nm范围内，四种煤体结构阶段比表面积与煤的破碎程度呈正相关性，在2～10 nm范围内，碎粒结构和糜棱结构孔隙复杂性增强，对比表面积贡献率提升。综合分析随着煤体结构破坏程度增加，孔容和比表面积均增加，以50 nm为分界点，对于孔容的贡献主要是分布于50～300 nm的孔，而对于比表面积则是分布于2～50 nm的孔。构造应力破坏了煤中大孔，形成了相对较多的小孔，导致大孔贡献率逐渐降低，而微孔和中孔（10～50 nm）的贡献比例逐渐增加。

图  3  不同煤体结构阶段孔径贡献孔容和比表面积百分比

Figure  3.  Differential pore volume and specific surface area contribution versus pore sizes

下载: 全尺寸图片 幻灯片

图  4  不同煤体结构阶段比表面积随孔径分布

Figure  4.  Specific surface area of coal structures versus pore sizes

下载: 全尺寸图片 幻灯片

2.2.2   超微孔分布

图5为不同煤体结构的超微孔孔径分布与阶段比表面积的关系，可将超微孔孔径分布划分为R₁（0.45～0.65 nm）、R₂（0.65～0.8 nm）和R₃（0.8～1.0 nm）3个孔径分布区间（表3），其中对不同煤体结构超微孔比表面积贡献较大的孔径峰值（P1，P2和P3）分布均位于这3个区间内，峰值1和2分布在R₁，峰值3位于R₃，不同煤体结构中的超微孔主要分布在0.45～0.65 nm和0.80～1.0 nm。但是CO₂测试结果显示不同煤体结构煤的超微孔孔径分布区间的范围以及各分布区间对应的累积比表面积值有显著差异。首先，随着煤破坏强度的增加，超微孔孔径变化范围越来越大，如碎裂结构煤的测试超微孔孔径介于0.54～0.93 nm之间，而在糜棱结构煤中，探测孔径甚至小于0.5 nm，变化范围为0.49～0.93 nm。此外，R₁和R₃区间对应的累计比表面积值随煤破碎程度的增加而增加，这些都表明原生结构煤在破坏过程中新形成的超微孔增多，特别是分布在0.65 nm以下的超微孔所占比例较高，贡献的比表面积值较大，例如，R₁区间内糜棱结构煤的累积比表面积为72.12 m²/g，是其它煤体结构的2倍左右。

图  5  不同煤体结构超微孔阶段比表面积随孔径分布变化

Figure  5.  Differential ultra-micropore specific surface area versus pore size of coal structures

下载: 全尺寸图片 幻灯片

表  3  不同煤体结构超微孔分布区间及对应的累积比表面积值

Table  3.  Distribution of ultra-micropore and cumulative SSAs in different intervals

样品	孔径R1			孔径R2			孔径R3
	区间/nm	比表面积/(m2·g−1)		区间/nm	比表面积/(m2·g−1)		区间/nm	比表面积/(m2·g−1)
原生结构	0.54～0.65	34.36		0.66～0.80	9.60		0.80～0.91	14.12
碎裂结构	0.54～0.63	33.88		0.63～0.80	13.99		0.80～0.93	13.96
碎粒结构	0.52～0.63	43.77		0.63～0.86	20.38		0.87～0.99	14.28
糜棱结构	0.49～0.65	72.12		0.66～0.82	15.11		0.82～0.93	18.57

下载: 导出CSV 

| 显示表格

2.3   不同煤体结构孔隙的分形规律

图6为不同煤体结构基于液氮和FHH模型的分形曲线，图中显示在ln（ln（P₀/P））=−0.5处有2个明显的线性分界段。通过曲线拟合得到了在相对压力区间0～0.5和0.5～1内的分形维数D₁和D₂，由于吸附机理和测试方法的不同，它们分别代表了煤中大于2 nm的孔表面粗糙度和孔隙结构的分形特征。研究发现，在较低压力下，D₁表征了微孔中气体分子与孔表面之间范德华力的作用，但随着压力增大，孔隙表面吸附多层气体分子，由于大量吸附质多层覆盖，表面光滑性增大，此外由于中孔和大孔的毛细凝聚作用，D₁不能反映孔结构特性，而与表面粗糙度有关，而分形维数D₂与孔结构均质性相关^[25]。D₁值越大，煤的表面越粗糙，D₂值越小，煤的孔隙结构分布越均匀，复杂程度越低。不同煤体结构计算的D₁和D₂见表4，D₁和D₂的值都介于2～3。D₁的变化范围为2.7291～2.921，随着煤体结构破坏强度的增大而增大，说明构造煤孔隙表面在破坏过程中变的比原生结构煤更不规则、更粗糙。D₂由2.6336下降到2.4077，整体上D₂值随煤体结构破坏程度的增加呈降低趋势，但D₂最低的是碎粒结构煤，这可能与数据拟合度较低有关（R²=0.250）。D₂值越大表明原生结构煤的孔隙结构越复杂，而经历构造变形后的煤孔隙结构将被简单化，孔径分布变窄。上述分析表明，构造破坏严重的煤表面粗糙度加剧，但孔隙结构非均质性逐渐减小。

图  6  不同煤体结构大于2 nm孔隙分形曲线

Figure  6.  Fractal analysis of coal structures based on N₂ isotherms

下载: 全尺寸图片 幻灯片

依据CO₂吸附量和相对压力的数据拟合关系，图7得到了不同煤体结构超微孔表面分形维数D_m，从原生结构煤到糜棱结构煤，超微孔分形维数D_m变化从2.4033～2.4249（表4），呈现上升的趋势，随着煤体结构破坏程度的增加，超微孔界面几何结构复杂性增强，这一定程度上影响了不同煤体结构煤的吸附性。

表  4  基于液氮和二氧化碳吸附数据计算的分形维数

Table  4.  Fractal dimension (D_1, D₂ and D_m) based adsorption analysis calculation

样品	P/P0<0.5	P/P0>0.5	D1	D2	Dm
原生结构	y = −0.2709x + 0.364	y = −0.3554x + 0.2532	2.7291	2.6336	2.4033
碎裂结构	y = −0.2399x + 0.3691	y = −0.4413x + 0.173	2.7601	2.5587	2.4071
碎粒结构	y = −0.079x + 0.2946	y = −0.5923x − 0.2843	2.921	2.4077	2.4096
糜棱结构	y = −0.0832x + 0.1773	y = −0.5487x − 0.1847	2.9168	2.4513	2.4249

下载: 导出CSV 

| 显示表格

图  7  不同煤体结构小于2 nm孔隙分形曲线

Figure  7.  Fractal curves of pores with different coal structures less than 2 nm

下载: 全尺寸图片 幻灯片

2.4   孔隙分形特征对煤中甲烷吸附、扩散影响

甲烷主要吸附在煤基质孔隙的内表面，尤其是超微孔、微孔和中孔为煤中甲烷提供了主要吸附空间。不同的煤体结构之间，不同孔径贡献的吸附比表面积分布会有明显的变化，一定程度上会影响不同煤体结构储层吸附能力。超微孔贡献的比表面积值要远大于微孔、中孔和大孔，为煤中的主要的吸附孔。超微孔分形维数D_m与Langmuir V_L呈正相关性，相关系数为0.8225（图8a），即D_m越大，吸附性越强。超微孔分形维数D_m与孔隙表面结构有关，随着煤体结构破坏程度增强，超微孔分形维数D_m增大，构造煤孔隙表面粗糙性增强，煤中可提供的吸附点位增加，吸附能力提升。超微孔分形维数D_m也与小于2 nm的孔隙贡献的比表面积值呈正相关性（图8a），相关系数为0.9563，进一步说明超微孔分形维数D_m可定量表征煤中内表面积的几何特征。

图  8  分形维数与Langmuir V_L、扩散系数、比表面积和孔容关系

Figure  8.  Relationships between fractal dimensions and Langmuir V_L, diffusion coefficient, specific surface area and pore volume

下载: 全尺寸图片 幻灯片

煤储层压力降低至临界解吸压力之下，煤中气体发生解吸，并立即进入扩散阶段。微孔和中孔作为扩散过程的主要空间，影响着煤中气体的扩散行为^[26]。不同煤体结构孔隙结构分形维数与有效扩散系数关系如图8b、图8c，图中显示分形维数D₁与有效扩散系数呈正相关性，相关系数分别为0.6233，而分形维数D₂与有效扩散系数呈负相关性，相关系数为0.7634。分形维数D₁增大，表征煤中微孔和超微孔贡献的吸附比表面积增加，破坏较严重的煤体结构中可吸附的气体量增加，一旦发生解吸，煤中气体解吸量较高，扩散空间内气体浓度差相对较高，扩散速率较快。此外，分形维数D₁与孔容呈正相关性，相关系数为0.9166。由于煤中气体总是选择从阻力最小的空间发生运移，构造煤破坏变形后，不同尺度孔隙空间的孔容均增大，扩散阻力降低，扩散系数增大。分形维数D₂与孔容呈负相关性，即分形维数D₂增大，煤中孔隙结构复杂性增强，连通性变差，降低了气体进出孔隙的效率。随着煤体结构破坏强度的增大，分形维数D₂降低，煤中孔隙将被简单化且连互相通，孔隙容积增大，气体扩散效率增高，有效扩散系数增大。

3.   结　　论

1）N₂和CO₂的吸附试验分析可以相对较完整地表征不同煤体结构不同尺度下的孔隙结构特征。煤的孔隙容积主要由中孔（50～100 nm）和大孔（>100 nm）贡献，对比表面积的贡献则主要是超微孔（<2 nm）、微孔（2～10 nm）和中孔（10～50 nm）。

2）随着煤体结构破坏强度的增加，50～300 nm的孔隙所占比例逐渐减小，2～50 nm的微孔、中孔和小于2 nm的超微孔增加，超微孔主要分布在0.45～0.65 nm和0.80～1.0 nm。

3）N₂和CO₂的吸附量随煤体结构破坏程度的增大而增加，吸附由大到小顺序为：糜棱结构煤>碎粒结构煤>碎裂结构煤>原生结构煤。

4）D₁和D₂分别表征大于2 nm孔隙表面和结构分形维数，D_m表征小于2 nm孔隙界面分形维数；随着煤体结构破坏程度的增加，D₁增大，D₂降低，D_m增大，原生结构煤经历构造变形后的孔隙结构将被简单化，孔隙表面粗糙度增强。D_m与Langmuir V_L和对应的表面积程正相关性，煤中吸附孔主要由超微孔贡献；有效扩散系数、孔容与分形维数D₁呈正相关性，与D₂成负相关性，表明经历构造破坏后煤孔隙结构连通性好，孔隙容积大，气体进出效率高，扩散效率增大。