不同煤级煤液相侵入效应低场核磁共振实验研究

摘要:我国煤层普遍具有低透气性、高瓦斯含量的特性，在低透气性煤层增透方面，煤层注水、水力压裂、水力割缝等水力化技术得到了广泛应用，并取得了良好的瓦斯治理效果。然而水作为外侵液进入煤体,堵塞了瓦斯流动通道，降低了瓦斯解吸量，产生了水锁效应。为分析水力化技术造成水锁效应的内在机理，利用压汞实验分析了煤样孔容分布规律，以及利用扫描电镜分析了原始、饱水、饱CMC溶液煤样的微观结构，基于低场核磁共振技术研究了煤样在饱水状态以及饱CMC溶液状态下的液相滞留效应,并根据曲线相似度法分析了孔径与束缚流体饱和度的相似度。研究结果表明：CMC溶液可以溶解煤中的矿物质增加煤孔隙裂隙以及降低水在煤体表面的表面张力，从而达到解除水锁效应的目的；随着煤变质程度的增大,丁2截止值在逐渐减小,4截止值的数值与煤样孔径大小呈负相关;煤样的束缚流体饱和度远大于自由流体饱和度,煤样在饱水状态下的束缚流体饱和度比饱CMC溶液状态下高；高变质程度的煤大孔孔容少、微孔孔容多，使得水在煤孔隙中的毛细管力大,最终造成高阶煤的水锁效应严重;大孔孔容是影响束缚流体饱和度的主控因素,微孔起到正向促进作用，得到束缚流体饱和度S与大孔孔容*、微孔孔容*的耦合关系式：S=94.861 078.96VA+261.24VDO滞留在煤体内的束缚水阻塞了瓦斯流动通道，是造成水锁效应的根本原因，增加煤层的孔隙裂隙以及选用合适的表面活性是减缓煤层水锁效应的有效措施。

关键词：水力化技术；水领效应；核磁共振;束缚流体;液相滞留效应

 

我国大部分煤层具有低渗透性、高瓦斯的特点⑴，随着煤矿开采深度的增加，煤层的透气性将进一步减小⑵，这为瓦斯灾害的防治带来了新的困难。在低透气性煤层增透方面，水力化技术得到广泛应用,如煤层注水⑶、水力压裂⑷、水力割缝⑸、水力掏槽⑹、水力冲孔⑺、水力挤岀⑻等。袁亮院士⑼指出水力化技术是深部高瓦斯突出煤层卸压增透及区域瓦斯治理的关键技术。然而，水力化技术引起的煤层水锁效应会降低瓦斯解吸量、阻碍瓦斯流动,这也引起了相关学者的关注[10-12]。水锁效应的研究最早出现在石油领域[13]，它是指外来流体沿着孔隙通道进入储层深部并在油气/水界面形成附加的毛细管阻力，如果储层能量不足以克服毛细管阻力，外来流体就会堵塞通道,导致储层较无外液侵入条件下渗透率降低,从而导致油气回采率降低的现象[14]。

国内外学者针对水锁效应进行了深入分析,特别是在油气藏水锁效应方面。CHAKRABORTY等发现水锁效应对地层渗透率产生重要影响，严重时将会导致页岩的渗透率下降99.5%。BIKASH等〔⑹发现油藏储层岩石表面的润湿性影响地层的渗透率。贺成祖等'切指出毛细管力是影响水锁效应的最主要的因素。吕渐江等〔⑻利用相渗曲线来研究低渗储层的水锁效应。在煤层水锁效应研究方面，张国华等研究渗透剂溶液侵入条件下对瓦斯自然解吸的影响，发现外液的侵入能够降低瓦斯的解吸量'明，并提出了利用水锁效应治理工作面超限的思路I測。朱借等⑵］利用不同表面活性剂水基溶液阻止煤体内部孔隙内的吸附瓦斯转变为游离瓦斯，以降低落煤瓦斯涌出速度。胡友林等⑵［从煤的多孔性、外来液侵入、毛细管力3个方面分析了煤层水锁机理，利用表面活性剂通过降低侵入液表面张力、增大接触角两个方面减少对煤层的水锁伤害。不论油气藏还是煤层，水锁效应皆降低了其渗透率，影响了油气资源的开釆效率,广大学者对其影响因素进行了深入分析和研究。针对煤层，相关学者利用水锁效应封堵瓦斯运移通道,将不利效应转变为瓦斯治理措施，是一种技术创新。然而，水力化措施的根本目的在于抽采瓦斯,降低煤层瓦斯含量,最大化利用瓦斯资源，因此解除煤层水锁效应才能充分发挥水力化措施的作用。

国内外学者研究水锁伤害时，通常使用外来流体滞留前后地层渗透率下降的相对值、储层初始含水饱和度、界面张力等参数分析水锁效应。然而,煤层水锁伤害的本质在于外液侵入煤体内部孔隙后阻塞瓦斯运移通道,造成瓦斯难以被抽出，这些参数难以体现煤层水锁效应的本质。目前，关于煤矿水锁效应基础研究方面正处于起步阶段，深入研究煤层水锁效应机理对矿井瓦斯防治具有十分重要的意义。笔者利用低场核磁共振技术，研究煤样在饱水状态下的液相滞留效应,分析束缚流体饱和度随煤变质程度的变化规律,揭示煤层水锁效应产生的内在机理，为解除水锁效应提供理论依据。

1实验部分

1.1材料与试剂

本实验选用蒸憶水对煤样进行饱和处理,模拟煤样在水入侵条件下产生的水锁效应，选择具有缓解水锁效应的梭甲基纤维素作为对比试剂⑵］，实验所用试剂见表1。

 

1.2实验样品

本实验采用6种煤样，分别来自淮南顾桥矿(GQ)、义马千秋矿(QQ)、枣庄付村矿(FC)、平顶山首山一矿(SS)、登封新丰矿(XF)、晋城大宁矿(DN)。利用长沙开元仪器股份有限公司生产的5E-MAC IV全自动工业分析仪进行煤样的工业分析,利用德国蔡司(ZEISS)显微镜光度计分析煤样的镜质组反射率,测试结果见表2。

 

1.3煤样饱水(液)处理

将煤矿现场取回来的新鲜煤样制成巾25 mmx 50 mm的圆柱，同种煤样制备10个，为保证试样的均一性,钻取的10个煤样均来自同1块煤炭，保证实验试件无裂隙、缝隙，在视觉上无差别,并做好标记。实验利用NM-VSD型真空饱和装置对煤样进行饱水(液)处理，装置原理图如图1所示。饱水(液)步骤如下：

(1)对煤样进行编号,在温度60龙的干燥箱中干燥24 h；

(2)分别从6种煤样中随机选择1块，放入洁净的烧杯，并将烧杯放入饱和缸内，利用真空密封胶将饱和缸与真空饱和装置密封；

(3)关闭进液开关，打开抽气开关，利用真空饱和装置对煤样抽真空8 h；

(4)在储液罐中加入适量蒸储水，进而关闭抽气开(4)关,打开进水开关，等水不再进入烧杯内时,关闭进液开关，再次打开抽气开关，再抽真空4 h,使煤样完全饱水；

(5)更换新煤样，将蒸储水换为质量分数为0.4%的CMC溶液，重复步骤(1)-(4),将新煤样再次饱和。

 

1.4实验方法

煤样微观孔隙结构

釆用压汞仪和扫描电镜研究煤样的微观孔隙结构，压汞实验分析6种原始煤样的孔隙结构,煤样为粒径3~6 mm；利用扫描电镜分析GQ原始煤样、饱水煤样、饱CMC煤样的微观形貌。压汞仪为美国Quantachrome Instruments公司的PoreMaster 33全自动压汞仪，扫描电镜为美国FEI公司生产的Quanta 250环境扫描电子显微镜。

(2)测试设备及方法

核磁共振测试采用纽迈公司生产的MMR-60型核磁共振分析仪，试样采用1.3节中饱水(液)处理过的煤样，核磁共振仪参数设置见表3。

2结果及分析

2.1煤样微观孔隙特征

2.1.1压汞实验

通过压汞实验，得到了6种煤样的孔容分布见表4,从表4可以看出，小孔所占比例最大，高达56%~62%,大孔所占比例在4.36%-14.29%,中孔所占比例在6.55%~11.98%，微孔所占比例在11.52%~30.55%。DN煤样的总孔容最大，QQ煤样的总孔容最小。随着变质程度的提高，煤的大孔孔容减少，微孔孔容以及总孔容增大，小孔孔容先减小后增大，而中孔孔容变化规律不明显。

 

2.1.2扫描电镜

通过扫描电镜得到了GQ煤样的表面微观特征，结果如图2所示。从图2中可以看出，原始煤样的表面矿物质较多，且填充到煤的裂隙孔隙如图2(a)所示；经过饱水处理后，煤样表面的部发育，分矿物质被溶解，裂隙边缘钝化，如图2(b)所示；经过饱CMC溶液处理后，煤样表面的矿物质不仅中，可以被溶解，煤基质内部的部分矿物质也被溶解,使更多的矿物质暴露出来，煤的棱角也更加平滑,如图2(c)所示。CMC可以溶解矿物晶体使孔隙增加了煤的孔隙尺寸和体积，这对降低水锁效应具有重要意义。

 

在煤层卸压增透方面，水力化技术得到了广泛应用,水作为媒介对煤体产生压裂、割缝、冲刷等作用,并进入煤体孔隙,从而产生水锁效应。为解除水锁效应，可以添加CMC改变水的性质,CMC一方面可以溶解煤中的矿物质，增加煤孔隙裂隙；另一方面改变水在煤体表面的润湿性，降低表面张力。这两方面的作用都可以降低水在煤孔隙中的毛细管力，从根本上达到解除水锁效应的目的。

2.2液相滞留效应

为了分析水力化技术中煤层的液相滞留效应,利用核磁共振设备分别对6种煤样在饱水、饱CMC溶液状态下进行测试，结果如图3所示。

 

从图3可以看出,6种煤样的弛豫时间匕曲线呈现出3个峰，前2个峰高而尖锐，第3个峰矮而平缓。煤样在饱水状态下形成的弛豫时间4曲线比饱CMC溶液的曲线要高，与横坐标所围的面积要大。图3中纵坐标的幅度表示信号的强弱,这种信号的强弱表示煤样中氢原子核'H含量的多少，亦即表明煤中水的多少。饱水状态下，测得4曲线高，说明煤样中水含量高，亦即说明在饱水状态下,水进入煤体的数量多，液相滞留效应明显，水锁伤害越严重。匕截止值将匕分布曲线分为束缚流体和自由流体两部分°】，前者所占的比例称为束缚流体饱和度,后者所占的比例称为自由流体饱和度。为了更好的分析这种现象，利用核磁共振软件得到T2截止值以及流体饱和度，所得结果见表5。

由表5可知,束缚流体饱和度远大于自由流体饱和度，束缚流体饱和度在93%以上，而自由流体饱和度则小于7%。煤样在饱水状态下所测的束缚流体饱和度高于饱CMC溶液,这说明在饱水条件下，滞留在煤体中的水要多，这也意味着水锁效应也越严重。随着煤变质程度的增大,匕截止值逐渐减小,GQ煤样的％截止值最大，为21.54 ms；DN煤样的T2截止值最小，为14.17 ms0结合煤样孔隙结构的变化规律，从中可以推断出,煤的孔径越小,4截止值就越小。为了更加直观的表达变质程度对匕截止值和束缚流体饱和度的影响，现将其随镜质组反射率的变化关系作图，结果如图4所示。

不同煤阶煤在饱水下&截止值和束缚流体饱和度Fig.4 T2 cutoff value and irreducible fluid saturation of different samples of coals under the saturation of water由图4可知,随着变质程度的增大，煤样的束缚流体饱和度呈现先上升后趋于平缓的趋势，这说明高阶煤比低阶煤更容易产生水锁效应。根据前人研究结果可知[，7'25],水锁效应产生的根本原因是毛细管力。其中,毛细管力的大小可由任意界面拉普拉斯方程网表示，即

 

式中，Pc为毛细管力，Pa；b为溶液表面张力，mN/m；0为接触角，（。）；r为孔隙半径,nm。

结合压汞实验分析结果可以看出,高变质程度的煤大孔孔容少、微孔孔容多，使得水在煤孔隙中的毛细管力大,最终造成高阶煤的水锁效应严重。

2.3孔径与束缚流体饱和度关联系分析

曲线相似度法根据两条曲线的相似程度判断对应变量之间的关联性⑵】，本文采用此法分析孔径对束缚流体饱和度的影响程度。在已知因变量变化趋势前提下,可以利用曲线相似度法分析哪个自变量是控制因变量的主控因素。这其中的主要思想是通过对已知数据进行数学变换，得到相应的向量式，根据因变量与自变量向量式的相似度的大小，研究控制因变量变化的主控因素。

曲线相似度法确定主控因素，主要包括以下4个步骤：

（1） 确定因变量:为分析束缚流体饱和度与孔容特性的关联性，将CMC溶液状态下测定的束缚流体饱和度作为因变量进行分析。

（2） 数据向量化处理：束缚流体饱和度S=[92.12 93.38 93.92 94.57 95.43 95.66]'「，大孔（A）孔容Va=[0.003 1 0.002 3 0.001 8 0.001 5 0.001 4 0.001 2]t,中孔（B）孔容VB=[0.002 6 0.001 8 0.001 8 0.001 7 0.002 9 0.001 8卩,小孔（（：）孔容Vc=[0.013 5 0.011 6 0.011 2 0.012 3 0.015 1 0.016 1卩,微孔（D）孔容V°=[0.002 5 0.003 2 0.004 6 0.005 1 0.007 2 0.008 4]T,总孔容（E）孔容VE=[0.002 5 0.003 2 0.004 6 0.005 1 0.007 2 0.008 4卩。

（3） 求列向量的相似程度。

（4） 主控因素确定:通过表6可以看出，从各级孔径孔容对束缚流体饱和度的影响程度来看,大孔孔容的变化趋势与束缚流体饱和度变化趋势相似度最高,高达-0.972 4,说明两者呈现较好的负相关性,两者的相似程度要高于其他孔隙孔容和总孔容。由此推断，大孔孔隙是影响束缚流体饱和度的主控因素。同时可以看出，微孔、总孔容与束缚流体饱和度的正相关性很高，皆为0.958 50为分析大孔、微孔、总孔容单一因素对束缚流体的影响，利用Origin进行线性拟合分析，结果如图5所示。

 

通过图5可以看出，束缚流体饱和度随着大孔孔容的增加呈线性降低趋势,随着微孔孔容的增加呈线性增加趋势。这说明束缚流体饱和度跟煤样的孔隙结构相关,微孔数量越多，束缚流体饱和度越大，进入煤体的水的数量越多，水锁效应也就越严重。虽然束缚流体饱和度与总孔容的相关性较高,但线性关系较差，因此在后续关联分析中不再考虑总孔容对束缚流体饱和度的影响。束缚流体饱和度的主控因素是大孔数量，但微孔数量同样起到重要作用，为分析两者对束缚流体饱和度的影响，利用SPSS对不同煤阶煤的大孔孔容、微孔孔容与束缚流体饱和度进行回归拟合，得到束缚流体饱和度与大孔孔容、微孔孔容的耦合关系式：S=94.86-1 078.96/+261.24昨，砂=0.9870通过上述耦合关系式可以看出，大孔孔容的权值要远远大于微孔孔容，这也反映出大孔孔容是决定束缚流体饱和度的决定性因素。

2.4基于液相滞留效应的煤层水锁机理分析

煤是一种复杂的多孔介质,具有双重孔隙结构,水在外界压力作用下沿着煤的孔隙通道进入煤体内部，水通过竞争吸附使吸附态的瓦斯变为游离态，游离态的瓦斯遇到水形成弯向瓦斯的水/瓦斯界面，在表面张力作用下,弯液面上会形成一个指向瓦斯的毛细管力。当煤层能量以及抽采负压不足以克服毛细管力时，形成水锁效应'闾，瓦斯就不会被抽出。从中也可以看出，水最终进入煤体内部的动力是毛细管力，毛细管力越大,水进入煤体的数量就会越多，滞留在煤体内部的水分就越多。滞留在煤体内部水分的多少可以通过束缚流体饱和度的大小来反映,进而可以通过束缚流体饱和度来衡量煤层水锁效应的大小。

3结论

(1)随着煤变质程度的增大，煤的微孔数量呈增加趋势截止值在逐渐减小;CMC溶液通过溶解煤中矿物质增加孔隙，以及改变煤表面润来降低水在煤孔隙中的毛细管力。

(2)煤样的束缚流体饱和度远大于自由流体饱和度，束缚流体滞留在煤体内阻塞了瓦斯流动通道,这是引起水锁效应的根本原因,可以用束缚流体饱和度来衡量煤层水锁效应的大小。

(3)大孔孔容是影响束缚流体饱和度的主控因素,微孔对束缚流体饱和度起到正向促进作用，得到束缚流体饱和度与大孔孔容、微孔孔容的耦合关系式：S=94.86-1 078.96VA+261.24VD。

 

参考文献(References):

[1]YUAN Liang.Theory and practice of integrated coal production and gas extractionf J].International Journal of Coal Science&Technolo-gy,2015,2(1):3-11.

[2]蔡永博，王凯，袁亮，等.深部煤岩体卸荷损伤变形演化特征数值模拟及验证[JL煤炭学报,2019,44(5)：1527-1535.CAI Yongbo,WANG Kai,YUAN Liang,et al.Numerical simulation and verification of unloading damage evolution characteristics of coal and rock mass during deep mining]J].Journal of China Coal Socie-ty,2019,44(5):1527-1535.

[3]赵振保.变频脉冲式煤层注水技术研究[J].采矿与安全工程学报,2008,25(4):486-489.

ZHAO Zhenbao.Study o£technology of variable-frequency pulse water iniiision into coal seam[J].Journal of Mining&Safety Engineering,2008,25(4)；486-489.

[4]ZHAI Cheng,LI Ming,SUN Chen,et al.Guiding-controlling technology of coal seam hydraulic fracturing fractures extension]J].International Journal of Mining Science and Technology,2012,22(6)：831-836.

[5]SHEN Chunming,LIN Baiquan,ZHANG Qizhi,et al.Induced drillspray during hydraulic slotting of a coal seam and its influence on gas extraction[J].International Journal of Mining Science and Technology,2012,22(6)：785-791.

[6]魏国营，郭中海,谢伦荣，等•煤巷掘进水力掏槽防治煤与瓦斯突出技术[J].煤炭学报,2007,32(2)：172-176.

WEI Guoying,GUO Zhonghai,XIE Lunrong,et al.Hydraulic slotting technology to prevent coal and gas outburst during heading excavation[J].Journal of China Coal Society,2007,32(2):172-176.

[7]郝富昌，孙丽娟,左伟芹.考虑流变特性的水力冲孔孔径变化规律及防堵孔技术[J],煤炭学报,2016,41(6)：1434-1440.

HAO Fuchang,SUN Lijuan,ZUO Weiqin.Hydraulic flushing aperture variation and anti-blocking technology considering rheological property[J].Journal of China Coal Society,2016,41(6):14341440.

[8]陈向军，王兆丰,程远平，等.水力挤出消突技术在水井头煤矿掘巷中的应用〔JL煤炭科学技术,2012,40(3):49-52.

CHEN Xiangjun,WANG Zhaofeng,CHENG Yuanpingtet al.Application of outburst elimination technology with hydraulic pressing to seam gateway driving in shuijingtou mine[J].Coal Science and Technology,2012,40(3)：49-52.

[9]袁亮,林柏泉，杨威.我国煤矿水力化技术瓦斯治理研究进展及发展方向[J].煤炭科学技术,2015,43(1)=45-49.

YUAN Liang,LIN Baiquan,YANG Wei.Research progress and development direction of gas control with mine hydraulic technology in China coal mine[J].Coal Science and Technology,2015,43(1)：45-49.

[10]倪冠华，李钊，解宏超.基于核磁共振测试的煤层水锁效应解除方法[J].煤炭学报,2018,43(8):2280-2287.

NI Guanhua,LI Zhao,XIE Hongchao.Experimental research on the methods to remove water blocking effect based on nuclear magnetic resonance testmgf J].Journal of China Coal Society,2018,43(8):2280-2287.

[11]张国华，梁冰，侯凤才，等.渗透剂溶液水锁效应对瓦斯解吸影响实验研究[JL科技导报,2011,29(22)=45-50.

ZHANG Guohua,LIANG Bing,HOU Fengcai,et al.Impact of penetrant solution water lock efiection on gas desorption]J].Science&Technology Review,2011,29(22)：45-50.

[12]吴家浩，王兆丰,苏伟伟，等.自吸水分对煤中瓦斯解吸的综合影响[J].煤田地质与勘探,2017,45⑴：35-40.

WU JiahaotWANG Zhaofeng,SU Weiwei,et al.Comprehensive influence of self-adsorbed water on gas desorption in coal[J].Coal Geology&Exploration,2017,45(1)：35-40.

[13]HOLDITCH S A.Factors affecting water blocking and gas flow from hydraulically fractured gas wells[J].Journal of Petroleum Technology,1979,31(12):1515-1524.

[14]赖南君，叶仲斌,刘向君，等.低渗透致密砂岩气藏水锁损害室内研究[J].天然气工业,2005,25(4),125-127.

LAI Nanjun,YE Zhongbin,LIU Xiangjun,et al.In-house study on water locking damage of tight sand gas reservoirs with low permeability[J].Natural Gas Industry,2005,25(4)•125-127.

[15]CHAKRABORTY N,KARPYN Z T,LIU S,et al.Penneability evolution of shale during spontaneous imbibition[J].Journal of Natural Gas Science and Engineering,2017,38：590-596.

[16]BIKASH D Saikia,JAGANNATHAN Mahadevan,DANDINA N Rao.Exploring mechanisms for wettability alteraticm in low-salinity waterfloods in carbonate rocks[J].Journal of Petroleum Science and Engineering,2018,164：595-602.

[17]贺承祖，胡文才.浅谈水锁效应与储层伤害[J].天然气工业，1994,14(6):36-38.

HE Chengzu,HU Wencai.Elementary lecture on water lock effect and reservoir damage]J].Natural Gas Industry,1994,14(6)：36-38.

[18]吕渐江，唐海，吕栋梁，等.利用相渗曲线研究低渗气藏水锁效应的新方法[J].天然气勘探与开发,2008,31(3)：7,58-61.Lt)Jianjiang,TANG Hai,LC Dongliang,et al.New research method of water blocking effect in low penneability gas reservoir by relativepermeability curve[J].Natural Gas Exploration and Development,2008,31(3):7,58-61.

[19]张国华，梁冰，毕业武.水锁对含瓦斯煤体的瓦斯解吸的影响[J].煤炭学报,2012,37(2)：253-258.

ZHANG Guohua,LIANG Bing,BI Yewu.Impact of water lock on gas desorption of coal with gas[J].Journal of China Coal Society,2012,37(2)：253-258.

〔20]张国华，鲁婷，梁冰，等.基于水锁机理的瓦斯超限防治理论〔J〕.黑龙江科技学院学报,2010,20(2):103-106.

ZHANG Guohua,LU Ting,LIANG Bing,et al.Theory on preventing and controlling gas exceeding limit based on water lock mechanism[J].Journal of Heilongjiang University of Science and Technology,2010,20(2):103-106.

[21]朱借，张标，邱泉玉.表面活性剂降低瓦斯涌出强度的性能指标研究[J].华北科技学院学报,2010,7(1)：32-35.

ZHU Kai,ZHANG Biao,QIU Quanyu.Researchers on performance of surfactant about decreasing gas emission intensity]J].Journal of North China Institute of Science and Technology,2010,7(1)：32-35.

[22]胡友林，乌效鸣.煤层气储层水锁损害机理及防水锁剂的研究[J].煤炭学报,2014,39(6):1107-1111.

HU Youlin,WU Xiaoming.Research on coalbed methane reservoir water blocking damage mechanism and anti-water blocking[J].Journal of China Coal Society,2014,39(6)：1107-lll 1.

[23]LIU Qian,GUO Yusen,AN Fenghua,et al.Water blocking effect caused by the use of hydraulic methods for permeability enhancement in coal seams and methods for its removal[J].International Journal of Mining Science and Technology,2016,26(4)♦615-621.

[24]张世民，孙银锁,章丽莎，等.核磁共振技术在非冻结土体孔隙分析中的应用[J].人民长江,2019,50(11)：183-188.

ZHANG Shimin,SUN Yinsuo,ZHANG Lisha,el al.Application of nuclear magnetic resonance technology in pore analysis of unfrozen soil[J].Yangtze River,2019,50(11)：183-188.

[25]崔晓飞.水锁伤害影响因素分析及解除方法研究现状[J].内蒙古石油化工,2011,37(10)：16T8.

CUI Xiaolei.Research of affective fact and cleaning up of waterblock damage in low penneability oil and gas field[J].Inner Mongolia Petrochemical Industry,2011,37(10):16-18.

[26]胡英.物理化学(第5版)[M].北京:高等教育出版社,2007：515-533.

[27]陈向军，刘军，王林，等.不同变质程度煤的孔径分布及其对吸附常数的影响[J].煤炭学报,2013,38(2):294-300.

CHEN Xiangjun,LIU Jun,WANG Lin,et al.Influence of pore size distribution of different metamorphic grade o£coalon adsorption constant[J].Journal of China Coal Society,2013,38(2)：294-300.