煤层水锁效应的消除及其对甲烷解吸特性的影响

刘 军，杨 通，王立国,2，陈向军

(1.河南理工大学 安全科学与工程学院，河南 焦作 454003；2. 河南省瓦斯地质与瓦斯治理重点试验室，河南 焦作 454003；3. 煤炭安全生产与清洁高效利用省部共建协同创新中心，河南 焦作 454003)

我国作为世界上已知的第三大煤层气储量国，总量约3.7×1013m3，拥有庞大的煤层气资源[1-2]。然而，我国煤层“三低一强”——低压力、低渗透率、低饱和度及非均质性强的特征导致瓦斯采前抽采困难[3-4]。目前，多数高瓦斯突出矿区采用水力压裂[5-6]、水力割缝[7-9]、水力挤出[10]以及水力冲孔[11-12]等水力措施作为一种增加煤层透气性、提高抽采效果的技术措施。

由于毛细管力引起水锁效应的产生，在水力化措施研究与应用过程中，在一定程度上抑制瓦斯抽采。从水锁机理上进行分析，含瓦斯煤为水锁效应提供了自然条件、物质条件和压力条件[13]。胡友林[14-16]、吴家浩[17]等通过试验对不同含水率煤样的解吸时间与渗透率的变化情况进行评价，从中分析煤层气储层水锁损害机理；宋金星等[18-19]提出压裂液中加入亲水性表面活性剂，减小压裂液与煤表面的接触角，降低毛管压力，增强压裂液的可排性，进而减缓水锁效应。然而，亲水型表面活性剂仅能在一定程度上减小毛细管力，并未消除毛细管力产生的阻碍作用，甚至实现将毛细管力由阻力向助力的转变，进而促进瓦斯的抽采；同时，对于孔径较小的微孔及纳米级的孔径而言，其产生的毛细管力仍较大，不利于瓦斯抽采。

在凝析气藏领域，LI等[20],FENG等[21],FAHES等[22]将气润湿反转剂用于解决凝析气藏出现的水锁损害，岩心润湿性实现了反转，由优先液湿转变为疏油疏水的优先气湿，且渗透率大幅度提高，为提高气井产量提供了重要的理论依据。通过利用气润湿反转技术将固体表面润湿性进行改变，消除毛细管力产生的水锁损害，进而增加气井产量。相比于凝析气藏的储层特征，煤层是孔隙结构更加复杂、渗透性更低的多孔介质；同时，煤体又具有亲水性的特征。在水力压裂过程中外在水进入煤体后，产生的水锁伤害更加严重。因此，采用气润湿反转技术提高煤层瓦斯抽采效果具有一定的借鉴性和可行性。

鉴于此，笔者引入气润湿反转技术，优选具有疏水性、低表面张力的气润湿反转剂，通过接触角和表面张力的测定，直观地验证表面活性剂处理前后的效果转变。然后测定1 h内不同甲烷平衡压力下干燥/含水煤样的甲烷解吸量和解吸速率，定量地分析气润湿反转消除水锁效应的影响效果。利用解吸公式进行拟合，预测极限解吸量，且拟合效果良好。揭示气润湿反转剂在防水锁增产的作用机理，为我国煤层瓦斯“增透促抽”技术提供一定的理论依据和指导。

气润湿反转剂具有疏水、低表面能的特征，通过物理吸附的形式附着在煤体表面，降低了煤体表面张力、增加接触角至90°以上，改变煤体表面特征，实现煤体表面润湿性的反转，进而消除水锁损害。毛细管力是评价水锁损害严重程度的关键指标，外来流体进入煤体等多孔介质后，毛细管力的大小和方向影响着水锁的严重程度以及返排过程中瓦斯的抽采难易。其中，毛细管压力表达式：

Pc=2σcosθ/r

(1)

式中，Pc为毛细管压力，Pa/cm；σ为表面张力，mN/cm；θ为接触角，(°)；r为毛细管半径，cm。

式(1)可知，毛细管半径作为煤体的固有特征，难以改变；利用表面活性剂可以实现煤体张力和接触角的改变。接触角和表面张力是影响毛细管力大小和方向的关键因素。润湿性的变化情况可以通过对接触角和表面张力的测定直观地呈现出来。

气润湿反转技术消除煤层水锁效应，增加瓦斯抽采效果的原理如下：由表面物理化学理论可知，固体的表面能越高，越容易被液体所润湿[23-24]；同时，表面自由能的大小也体现了固体表面吸附能力的强弱，表面能越高，其表面的吸附能力越强[25-27]。具有疏水、低表面能特征的气润湿反转剂分子中包含亲水性的羟基(—OH)和憎水性的高分子化合物氟碳烷基(—CF)。其中，亲水性的羟基(—OH)通过物理吸附的形式附着在煤体表面，而憎水性的氟碳烷基(—CF)朝向远离煤体表面的方向。气润湿反转剂分子在煤体表面的吸附，降低了煤体表面自由能，减弱了煤体表面对瓦斯的吸附能力，增加了游离瓦斯量；同时，由于煤体表面润湿性由液润湿转变为气润湿，毛细管力产生的方向发生变化，实现了“反毛细管效应”，从而由阻力变成动力，促进瓦斯的抽采。

选用的煤样取自晋城玉溪煤矿，其最大原始瓦斯含量为18.53～25.59 m3/t，最大原始压力为1.76～2.90 MPa，煤层透气性系数为0.134～0.26 m2/(MPa2·d)，百米钻孔自然初始瓦斯涌出强度为0.055 9～0.091 0 m3/(min·hm)，钻孔自然瓦斯流量衰减系数为0.042～0.046 d-1。利用不同孔径大小的标准筛进行筛分，筛选出粒度大小为200目以下的煤粉用于煤样片的制备；筛选出粒度大小为60～80目的煤样，用于煤样的工业分析试验和解吸试验。煤样工业分析参数见表1。

表1 煤样工业分析

笔者在前期进行不同类型表面活性剂的筛选，优选出具有疏水、低表面能的气润湿反转剂[28]。为进一步研究气润湿反转剂处理前后对煤体解吸特性的影响，在本文中采用的表面活性剂包括全氟烷基季胺碘化物(FC134)和氟碳表面活性剂(FC117)。

试验装置包括JMBJS20002型电子天平、压力机(自制)、电热鼓风干燥箱、JC2 000C1型接触角测量仪、含瓦斯多元气体置换系统。JMBJS20002型电子天平用于称量出质量为0.5 g的煤粉，确保压制出的煤样片形状、大小一致。压力机用于煤样片的压制工作，压力载荷设置为120 kN，稳压30 min，制成直径为25 mm，厚度为5 mm 的煤样片[29]。电热鼓风干燥箱用于煤样片和颗粒煤(60～80目)的干燥，其温度设置范围为0～400 ℃。JC2 000C1型接触角测量仪包括样品台、成像系统以及计算机与相应的检测软件系统，可实现接触角和表面张力的精准测量。含瓦斯多元气体置换系统包括3套独立的煤样室以及对应的传感器、参考罐、气瓶、真空泵、计算机与相应的软件系统。

通过JC2000C1型接触角测量仪及其分析系统，如图1所示，利用躺滴法测定2种气润湿反转剂在不同浓度梯度下处理后煤样的接触角以及利用悬滴法对表面活性剂溶液的表面张力进行测定[30]。其中，表面活性剂的浓度依次为0、0.2%、0.5%、0.8%、1.0%和1.5%。

图1 JC2000C1型接触角测量仪Fig.1 JC2000C1 contact angle meter

采用躺滴法进行接触角的测定时，为了避免水对气润湿反转剂反转结果的影响，首先利用喷散法将不同浓度的表面活性剂溶液喷洒在煤样片表面，放置在105 ℃鼓风干燥箱中干燥24 h以上，直至煤样质量不再变化，取出冷却备用[31]。同时，将剩余的表面活性剂溶液按照上述的操作，将不同浓度表面活性剂处理后的煤样按照浓度大小的不同进行排序，并贴上表面活性剂种类和浓度的标签，以便表面张力的测定和记录。

为了进一步研究气润湿反转剂对煤体解吸效果的影响，利用含瓦斯多元气体置换系统对气润湿反转剂消除煤层水锁效应、提高瓦斯抽采效果进行定量研究，具体试验步骤：

1)首先将取自玉溪煤矿煤样进行破碎，筛选出60～80目的颗粒煤6份。

2)选择3份煤样，分别浸泡在蒸馏水、FC117和FC134表面活性剂溶液中，浸泡48 h；然后放置在105 ℃鼓风干燥箱中干燥48 h以上，直至煤样质量不再变化，取出冷却备用。

3)检查含瓦斯多元气体置换系统中煤样室及管路的气密性、传感器与软件系统的运行情况，确保整个系统能够正常工作。

4)将煤样装进煤样罐，联通煤样罐与真空泵系统，进行抽真空工作24 h，抽至20 Pa以下。

5)抽真空结束之后，首先进行甲烷平衡压力为0.5 MPa的充气工作，吸附时间12 h 以上，且2 h内煤样罐中压力不发生变化。

7)依此类推进行甲烷平衡压力为1.0,1.5,2.0和2.5 MPa的充气工作，重复步骤5)和6)。

8)选择步骤1)的另外煤样3份，称量完全干燥条件下的煤样质量，计算其含水率为1%条件下的质量。分别浸泡在蒸馏水、FC117和FC134表面活性剂溶液中，浸泡48 h；然后放置在60 ℃鼓风干燥箱中干燥48 h，称量煤样质量，观察是否达到预定含水质量。若离预定质量较近时，每隔30 min，称量一次，直至达到预定的含水质量，取出冷却备用。

9)参照干燥煤样解吸步骤，重复步骤3)～7)。

3.1.1 接触角测试结果与分析

将经过不同浓度表面活性剂处理后的煤样片依次放置在接触角测量仪的样品台，采用蒸馏水作为滴定溶液，分别测定其接触角，如图2所示。

图2 接触角的测定Fig.2 Measurement of contactangle

如图2所示，FC117和FC134两种表面活性剂均实现了煤体润湿性的反转。当表面活性剂浓度从0增加至0.5%时，接触角逐渐增加，并在溶液浓度为0.5%时，润湿性由液润湿向气润湿的转变(接触角>90 °)；表面活性剂溶液浓度继续增加，接触角不再增加，维持在一个相对稳定的范围，甚至出现小幅度地下降。出现煤样接触角变化规律的主要原因为

1)当表面活性剂浓度从0增加至0.5%时，FC117和FC134两种表面活性剂处理后的煤体表面接触角增加速率较快。在该阶段，表面活性剂分子较少，能够有序的吸附在煤体表面，表面活性剂分子中的亲水基团朝着在煤体表面的方向吸附；疏水基团朝着远离煤体表面的方向，并将抑制周围水分子在煤体表面的吸附。随着表面活性剂浓度的增加，煤体表面吸附的表面活性剂分子逐渐增加，疏水基团在煤体表面抑制水分子的能力逐渐提高，煤体表面接触角也逐渐增大。

2)当表面活性剂浓度等于0.5%时，煤体表面吸附的表面活性剂分子趋于饱和，疏水基团在煤体表面抑制水分子的能力达到顶峰，煤体表面接触角也逐渐达到最大值；若继续增加表面活性剂浓度，多余的表面活性剂分子与正常吸附在煤体表面活性剂分子发生紊乱，疏水基团之间产生吸附、缠绕，使得多余表面活性剂分子中的亲水基团朝向背离煤体表面的方向，对水分子产生吸附作用，进而降低了煤体表面是疏水能力，煤体表面接触角会产生小幅度地降低。

对于多组测量数据,mij表示第 i 组的mj,由于每一组测量数据都有测量误差,故有相应的测量残差方程组如式(3)所示:

3.1.2 表面张力测试结果与分析

将经过不同浓度表面活性剂溶液分别装至用于滴定装置的注射器中，利用悬滴法对表面活性剂溶液进行表面张力的测定。每次滴定工作结束之后使用蒸馏水进行多次冲洗，避免测定溶液对待测溶液的结果产生影响，见表2。

表2 表面张力的测定

由表2可知，随着FC117和FC134两种表面活性剂溶液浓度的增加，溶液的表面张力先迅速降低，再逐渐趋于稳定。当表面活性剂浓度从0增加至0.2%时，表面张力从初始的73 mN/m分别降至23.124 mN/m和19.07 mN/m，表面张力大幅度地降低；当表面活性剂浓度从0.2%增加至1.5% 时，表面张力变化幅度较小，保持着较小程度的降低，并趋于稳定。

综上分析，对FC117和FC134两种表面活性剂从接触角和表面张力两个方面进行性能测定，两种表面活性剂均具有疏水性、低表面张力的性能，符合预期的效果。其中，当两种表面活性剂溶液浓度为0.5%时，可以实现煤样表面润湿性的转变，并具有较低的表面张力。因此，在后续的煤样解吸试验中，可将溶液浓度为0.5%的表面活性剂作为对解吸试验所用的煤样进行处理的最佳使用浓度。

根据2.4节试验步骤中的2)～5)，其中，2种表面活性剂溶液的使用浓度设置为0.5%。依次对气润湿反转剂处理前后的煤样进行不同吸附平衡压力条件下的吸附试验。根据《煤的甲烷吸附量测定方法》，计算充入煤样罐中的甲烷量Qci和吸附平衡后剩余空间的游离甲烷量Qdi，其表达式分别为

(2)

(3)

式中，Qci为充入煤样罐的甲烷标准体积，cm3；P1i,P2i为分别为充气前后参考罐内绝对压力，MPa；Z1i,Z2i分别为P1,P2压力下及t1时甲烷的压缩系数，1/MPa；t1为室内温度，℃；V0为参考罐及连通管标准体积，cm3；Vd为剩余体积，cm3；t3为试验温度，℃；

充入煤样罐的甲烷量扣除煤样罐中剩余体积的游离甲烷量即为压力段内煤样吸附甲烷量ΔQi：

ΔQi=Qci-Qdi

(4)

每克煤压力段内的吸附量为

(5)

式中，Gr为煤样可燃物质量，g。

通过计算得到气润湿反转剂处理前后干燥煤样在不同吸附平衡压力条件下的吸附量，如图3所示。

图3 干燥煤样吸附量Fig.3 Adsorption quantity of dry coal samples

由图3可知，随着甲烷吸附平衡压力的增加，气润湿反转剂处理前后煤样的吸附量均逐渐增加，且增加量均呈现逐渐减小趋势。经过气润湿反转剂FC117和FC134处理后的煤样吸附量始终小于未处理的煤样。其中，产生该现象的原因主要为：一方面，气润湿反转剂分子自身占有一定的体积，随水溶液一同进入煤体内部。在干燥过程中，水分通过蒸发作用排出煤体，而气润湿反转剂分子则滞留在煤体内部，进而减小甲烷分子的吸附空间；另一方面，气润湿反转剂具有低表面张力的特性，在吸附在煤体内部孔隙裂隙的表面，降低其表面自由能，减小对甲烷分子的吸附能力，进而使得吸附量下降。

由数学中的不定式性质可知，若在相同的吸附量的条件下，未处理煤样的解吸量和解吸速率小于处理后的煤样，则表明处理后的煤样解吸效果较好；而实际的吸附试验中，气润湿反转剂处理前煤样的吸附量大于处理后煤样的吸附量，则表明气润湿反转剂处理后的煤样解吸率更高、解吸效果更加显著。为了定量的分析气润湿反转剂对煤样解吸效果的影响，在下述进行的解吸量和解吸速率的比较中，均假定在相同的吸附量条件下进行定量比较。

为了得到气润湿反转剂处理前后干燥煤样的解吸规律，根据2.4节试验步骤中的5)～7)，依次对不同吸附平衡压力条件下的煤样进行解吸试验，如图4所示。

由图4a可知，随着甲烷吸附平衡压力的增加，3种处理方式的煤样解吸量均逐渐增加，且增加量均呈现逐渐减小趋势。经过气润湿反转剂FC117和FC134处理后的煤样解吸量始终大于未处理的煤样。在甲烷吸附平衡压力为0.5 MPa条件下，表面活性剂FC134的表现更加优异；在其余的甲烷吸附平衡压力，表面活性剂FC117的效果均更加显著。

图4 在1 h内干燥煤样的解吸规律Fig.4 Desorption law of dry coal samples within 1 hour

对干燥煤样进行解吸试验的测定，其目的主要是从定量的角度分析气润湿反转剂通过降低煤体表面自由能，增加游离瓦斯量，进而促进瓦斯解吸。其中，在不同吸附平衡条件下，气润湿反转剂处理后的煤样解吸量均有所增加，增长率基本上呈先增加再降低的趋势。在吸附平衡压力达到1.5 MPa 时，增长率达到最大值，分别达到近48%和34%，见表3。

表3 干燥煤样解吸量的增长率

低表面能是表面活性剂促进甲烷解吸的重要特征，其作用原理为：具有低表面自由能的表面活性剂分子通过物理吸附的方式有序的吸附在煤体内部的孔隙、裂隙表面，降低煤体表面对甲烷的吸附能力，增加游离瓦斯量。同时，煤体表面吸附能力的下降，减弱瓦斯流动过程中的阻力，有利于甲烷的解吸和运移。

甲烷气体在煤体内部孔隙、裂隙中解吸和运移过程中，主要受到煤体表面自身产生的摩擦阻力和气体本身解吸产生的动力共同作用。当甲烷吸附平衡压力较小时，煤体表面的摩擦阻力占主导作用。表面活性剂通过物理吸附的方式吸附在煤体表面，降低了煤体表面的吸附能力，增加了游离瓦斯量和气体在煤体表面运移的流平性。因此，当吸附平衡压力处于0.5～1.5 MPa时，随着平衡压力的增加，表面活性剂处理后的煤样解吸量增长率呈现增加的趋势。当吸附平衡压力处于1.5～2.5 MPa时，气体本身产生的动力占主导作用。随着吸附平衡压力的增加，甲烷解吸产生的动力逐渐增强，未处理的煤样受摩擦阻力产生的影响逐渐减小，因此，气润湿反转剂处理后的煤样解吸量增长率呈现逐渐减小的趋势。

通过对干燥煤样解吸量和解吸速率的分析，经过气润湿反转剂FC117和FC134处理后的煤样解吸量和解吸速率均有所提升。定量地分析了气润湿反转剂通过物理吸附的方式吸附在煤体表面，降低煤体表面自由能，增加煤样解吸量，为气润湿反转剂提高瓦斯抽采效果提供理论依据。

根据2.4节试验步骤中的8)～9)，依次对气润湿反转剂处理前后的含水煤样进行不同吸附平衡压力条件下的吸附试验。利用式(2)—式(5)，通过计算得到气润湿反转剂处理前后含水煤样在不同吸附平衡压力条件下的吸附量，如图5所示。

图5 含水煤样吸附量Fig.5 Adsorption quantity of water content coal samples

由图5可知，随着甲烷吸附平衡压力的增加，气润湿反转剂处理前后煤样的吸附量均逐渐增加，且增加量均呈现逐渐减小趋势。经过气润湿反转剂FC117和FC134处理后的煤样吸附量始终小于未处理的煤样，其变化规律与干燥煤样的吸附量一致。其中，含水煤样的吸附量略小于干燥煤样的吸附量，其主要原因是在含水率为1%的煤样中，少量的水滞留在煤体内部。一方面，由于水的体积侵占了原本甲烷分子存在的空间；另一方面，由于煤样自身的亲水性，水在煤体表面的吸附能力强于甲烷气体，使得部分吸附状态的甲烷分子被置换和驱替。因此，含水煤样的甲烷吸附量略小于干燥煤样的吸附量。

为了得到气润湿反转剂处理前后含水(1%含水率)煤样的解吸规律，依次对不同吸附平衡压力条件下的煤样进行解吸试验，通过对比干燥与含水煤样的解吸量和解吸速率，进而定量分析气润湿反转剂消除水锁效应，如图6所示。

由图6a可知，随着甲烷吸附平衡压力的增加，3种处理方式的含水煤样解吸量均逐渐增加，且增加量均呈现逐渐减小趋势。经过气润湿反转剂FC117和FC134处理后的含水煤样解吸量始终大于未处理的煤样，与干燥煤样解吸量的规律基本一致。同时，在3种不同处理方式条件下，干燥煤样的解吸量均大于含水条件下的解吸量，其结果表明在含水煤样解吸过程中存在一定程度的水锁效应。

对含水煤样进行解吸试验的测定，其目的主要是从定量的角度分析气润湿反转剂通过将煤体润湿性由亲水性转变为疏水性，消除因毛细管作用所引起的水锁效应，进而提高甲烷解吸效果。

其中，在不同吸附平衡条件下，气润湿反转剂处理后的含水煤样解吸量均有所增加，增长率基本上呈先增加再降低的趋势。在吸附平衡压力达到1.0 MPa 时，水锁损害程度最大，表面活性剂消除水锁效果最好，表面活性剂增加煤样解吸量最佳，增长率达到最大值，分别达到近38%和25%，见表4。煤样解吸量的变化规律与上述3.2节干燥煤样解吸量的变化规律基本上是一致的，作用机理亦与之相同。

表4 含水煤样解吸量的增长率

由图6b可知，随着甲烷吸附平衡压力的增加，3种处理方式的含水煤样解吸速率均有所增加。经过气润湿反转剂FC117和FC134处理后的煤样解吸速率也始终大于未处理的煤样。其中，表面活性剂FC117处理后的煤样解吸速率始终保持着显著的效果。同时，由于水锁效应产生的影响，不同处理方式的含水煤样解吸速率均小于干燥煤样的解吸速率。

图6 在1 h内含水与干燥煤样解吸规律的对比Fig.6 Comparison of desorption law between water content and dry coal samples within 1 hour

通过对干燥与含水煤样解吸量和解吸速率的分析，经过气润湿反转剂——FC117和FC134处理后的煤样解吸量和解吸速率均有所提升，其中FC117的效果更加显著。同时，经过表面活性剂处理后，煤样的润湿性发生转变，水锁效应得到缓解和消除。含水煤样的解吸量和解吸速率也均大于未处理的干燥煤样。定量地分析了表面活性剂改变煤体表面润湿性，由亲水性向疏水性转变，缓解水锁损害，进而增加煤样解吸量。

通过上述试验数据的对比，经过气润湿反转剂处理后，在1 h内煤样的解吸量和解吸速率均显著提高。为预测煤样的极限解吸量，可以根据煤样解吸量随时间的变化规律，利用拟合曲线得到最终的解吸量。

通过利用解吸公式——王佑安式对前1个小时的累计煤样解吸量与时间的关系曲线进行拟合，得到极限解吸量和拟合度，见表5和表6。

表5 气润湿反转前后煤样解吸量

表6 拟合度(R2)

由表5可知，利用王佑安式拟合不同平衡压力条件下的解吸量与时间关系曲线，拟合得到累计解吸量。其中，干燥煤样的解吸量大于对应含水煤样的解吸量，气润湿反转剂处理后的煤样解吸量大于处理前的煤样解吸量，拟合结果与1 h内解吸量的变化规律是一致的。同时，根据拟合结果可以预测气润湿反转剂处理前后的极限解吸量。

由表6可知，利用王佑安式拟合不同平衡压力条件下的解吸量与时间关系曲线，得到曲线的拟合度。其中，气润湿反转剂处理前后含水煤样的拟合效果最好，拟合度均大于0.99。在干燥煤样中，FC134处理后的煤样拟合度最好，不同吸附平衡压力条件下的拟合度均大于0.99；其次是未处理和FC117，拟合度也比较好，除去未处理煤样中吸附平衡压力为1.0 MPa,其余均在0.98以上。

1)分别对FC117和FC134处理后的煤样进行接触角和表面张力进行测定。其中，接触角从未处理的30°增加至90°以上，实现了煤体润湿性由亲水性向疏水性转变；表面活性剂溶液的表面张力由初始的72.8 mN/m下降至25 mN/m以下，较大幅度地降低了煤体表面张力。

2)通过接触角和表面张力的测定，当表面活性剂浓度为0.5%时，两种表面活性剂处理后的煤样可实现润湿性的转变和表面张力大幅度地降低。将吸附解吸试验中所用的煤样在该浓度下的表面活性剂溶液中进行处理，并将其作为最佳使用浓度。

3)通过对气润湿反转剂处理前后煤样吸附量的计算，对比两种表面活性剂处理前后的煤样吸附量变化情况。其中，经过FC117和FC134处理后的干燥/含水煤样均小于未处理煤样的吸附量。

4)对经过FC117和FC134处理后的干燥/含水煤样进行解吸量和解吸速率的测定，两种表面活性剂处理后煤样的解吸效果均有显著的效果。在气润湿反转剂处理后的煤样吸附量小于未处理的，而解吸量和解吸速率大于未处理的，则表明气润湿反转剂处理后的煤样可以消除水锁效应，增加解吸率，提高解吸效果。

5)利用解吸公式对解吸量进行拟合，得到极限解吸量和拟合度。其中，采用王佑安式拟合的效果较好，拟合度基本上都在0.98以上。

煤炭科学技术2022年9期