考虑煤基质压缩效应的煤全孔径分布特征研究

2022-03-04刘纪坤任棒王翠霞

工矿自动化 2022年2期

刘纪坤，任棒，王翠霞

(西安科技大学安全科学与工程学院，陕西西安 710054)

0 引言

煤是一种多孔介质，其孔隙结构直接影响煤中瓦斯富集和渗透特性[1-2]。目前，煤孔隙结构表征方法主要有压汞法[3]、低温液氮吸附法[4]、扫描电镜法[5]、显微CT法[6]、小角散射法[7]等。压汞法、低温液氮吸附法可将流体侵入孔隙内表面，但测量的孔径范围存在很大差异，压汞法可准确表征较大孔径的孔隙结构特征，低温液氮吸附法可准确表征较小孔径的孔隙结构特征；扫描电镜法、显微CT法、小角散射法通过图像处理可直观观察孔隙外部，但获取的孔径分布数据不具有代表性。

由于采用单一方法只能对某一孔径范围内的孔隙结构进行表征[8]，联合压汞法和低温液氮吸附法分析煤全孔径分布特征得到广泛应用。江泽标等[9]利用低温液氮吸附试验与压汞试验相结合的方法对CO2致裂前后的煤岩孔隙变化进行定量表征；郝晋伟等[10]基于IUPAC孔隙分类方案，联合使用压汞法、低温液氮吸附法对不同尺度的孔隙结构进行测试，将测试数据进行插值处理，从而获得构造煤孔隙结构的全孔径分布数据；林海飞等[11]采用同孔径处比表面积增量或孔体积增量差值最小的方法，将压汞法和低温液氮吸附法下孔径分布数据进行联合，分析了煤样全孔径段孔隙结构特征；何珊[12]以低温液氮吸附法和压汞法的优势孔径段临界值为连接点，对测量数据进行连接，从而实现煤样全孔径孔隙结构的有效表征。煤基质具有很强的可压缩性，当孔隙内压力增加时，煤基质产生变形甚至破坏，这种现象称为煤基质压缩效应[13-15]。然而上述研究均未考虑煤基质压缩效应对孔径分布测量结果的影响。因此，本文联合压汞法和低温液氮吸附法对不同变质程度煤样进行孔隙结构参数测定，计算煤基质压缩性系数，对孔径分布数据进行校正，进而采用联孔分析原则实现煤全孔径分布特征准确表征。

1 煤样选取及基础参数测定

实验所用4种煤样分别取自阜康矿区的东风福胜矿和鑫龙矿、彬长矿区的小庄矿及澄合矿区的合阳矿，取样煤层均为正在开采煤层。东风福胜矿地处二工河-黄山倒转向斜的东端，鑫龙矿位于阜康市甘河子镇以南7 km处，阜康矿区煤层赋存于侏罗系下统八道湾组地层，组内可采最大瓦斯含量为5.811 m3/t，瓦斯垂直分带位于-170～-450 m。小庄矿地处彬长矿区东部，煤层位于侏罗系中统延安组第一段，主采的4号煤层平均厚度为18.01 m，绝对瓦斯涌出量为4.11 m3/min。合阳矿地处渭北石炭二叠纪煤田，主采的5号煤层埋藏深度为520 m左右，平均厚度为3.5 m，矿井瓦斯含量较低，主要富含贫瘦煤。将采集的煤样密封保存，分别进行工业分析、显微分析和镜质组反射率测定，结果见表1。

表1 煤样基础参数测定结果Table 1 Test results of basic parameters of coal samples

2 压汞实验数据校正

2.1 低温液氮吸附和压汞实验结果对比分析

本文分别采用ASAP 2020比表面积和孔径分析仪器与AutoPoreIV 9510型全自动压汞仪器进行低温液氮吸附实验和压汞实验，2种实验测定的比表面积和孔体积分别见表2和表3(本文采用霍多特(XOJIOT)孔径划分方法，将孔径分为大孔(>1 000 nm)、中孔(100～1 000 nm)、过渡孔(10～100 nm)、微孔(<10 nm))。可看出在微孔段压汞法测得的比表面积是低温液氮吸附法测得结果的37.61～540.53倍，微孔段压汞法测得的孔体积是低温液氮吸附法测得结果的106～568倍，2种实验结果差异较大，主要原因是高压汞对煤基质造成了破坏，使得微孔测量误差增大，因此有必要对压汞实验数据进行校正。

表2 各孔径段比表面积Table 2 Specific surface area of each aperture segment (m2·g-1)

表3 各孔径段孔体积Table 3 Pore volume of each aperture segment (mL·g-1)

2.2 煤基质压缩性系数计算

煤基质压缩性系数是煤基质压缩效应的一个量化参数，对压汞实验数据校正前需要计算煤基质压缩性系数。假设煤基质压缩性系数是常数，根据式(1)计算煤基质压缩性系数，结果见表4。

(1)

式中：km为煤基质压缩性系数，m2/N；Vm为煤基质体积，mL/g；P为汞压，MPa。

2.3 数据校正

压汞法所测得累计进汞体积与汞压关系如图1

表4 煤基质压缩性系数Table 4 Compressibility coefficient of coal matrix

所示。当0.124 MPa≤P≤20 MPa时，累计进汞体积与汞压呈上凸递增关系，表明在该压力段内煤基质被压缩，孔隙结构发生变形，且随着汞压增加，煤基质压缩效应逐渐明显，孔隙结构变形程度逐渐增加；当20 MPa

图1 累计进汞体积与汞压关系曲线Fig.1 Relationship curve of cumulative mercury volume and mercury pressure

将累计进汞体积看作煤基质体积和孔体积2个部分之和，则

(2)

(3)

(4)

ΔVpi=Vobsi-Vobs0-kmVmi(Pi-P0)

(5)

式中：ΔVobs为累计进汞体积变化量，mL/g；ΔVobsi为第i(i=1，2，…，n，n为0.124 MPa≤P≤20 MPa内划分的压力段数)个压力段累计进汞体积变化量，mL/g；ΔVmi为第i个压力段煤基质体积变化量，mL/g；ΔVpi为第i个压力段孔体积变化量，mL/g；Pi为第i个压力段汞压，MPa；P0为初始汞压，MPa；ΔVm为煤基质体积变化量，mL/g；ΔP为汞压变化量，MPa；β为常数；ΔVp0为低温液氮吸附实验所测孔径为6～62.35 nm(根据Wasburn方程[16]计算可知，P=20 MPa对应的孔径为62.35 nm，P=206 MPa对应的孔径为6 nm)时的孔体积变化量，mL/g；Vobsi为第i个压力段累计进汞体积，mL/g；Vobs0为初始汞压下累计进汞体积，mL/g;Vmi为第i个压力段煤基质体积。

假设孔隙均为球状，比表面积变化量为

(6)

式中ΔSpi为第i个压力段比表面积变化量，m2/g。

3 煤样全孔径分布特征

3.1 联孔分析原则

联孔分析是指采用低温液氮吸附法和压汞法均能较为准确地测试出煤样孔隙结构且将测试结果在联孔分界点处进行联合。本文中联孔分界点确定为62.35 nm，采用以下联孔分析原则分析全孔径分布特征：当孔径小于62.35 nm，采用低温液氮吸附法对孔体积和比表面积进行分析；当孔径大于62.35 nm，通过对压汞实验数据校正来分析孔体积和比表面积。

3.2 全孔径比表面积分布特征

全孔径比表面积测定结果如图2所示。可看出煤样微孔比表面积为0.038 3～0.506 5 m2/g，过渡孔比表面积为0.172 8～0.370 5 m2/g，中孔比表面积为0.025 6～0.076 6 m2/g，大孔比表面积为0.002 2～0.009 6 m2/g；中低变质程度煤样(1,2,3号煤样)的比表面积以微孔和过渡孔为主，且过渡孔比表面积最大，微孔和过渡孔是影响中低变质程度煤样总比表面积的主要因素；与1,2,3号煤样相比，4号煤样各孔径段比表面积均有所增加，微孔比表面积增加幅度最大，表明微孔发育最为显著。

图2 全孔径比表面积Fig.2 Specific surface area of full aperture

全孔径比表面积分布特征如图3所示。可看出各煤样比表面积分布趋势大致相同；微孔和过渡孔比表面积增量均较大，表明微孔和过渡孔对比表面积的贡献最大；随着煤样变质程度增加，微孔和过渡孔比表面积增量增幅大于中孔及大孔，表明微孔和过渡孔对比表面积的贡献率增加。

图3 全孔径比表面积分布特征Fig.3 Distribution characteristics of specific surface area of full aperture

3.3 全孔径孔体积分布特征

全孔径孔体积测定结果如图4所示。可看出煤样微孔体积为0.000 044～0.000 554 mL/g，过渡孔体积为0.000 985～0.002 126 mL/g，中孔体积为0.000 87～0.004 51 mL/g，大孔体积为0.003 14～0.009 72 mL/g；大孔体积最大，随着孔径减小，孔体积大幅减小，中孔及大孔对孔体积起主导作用；随着煤样变质程度增加，各孔径体积均呈增大趋势，4号煤样孔体积最大。