中等變質煤孔隙特征隨揮發分變化規律

柴天興，靳樂軍

（山西陽城陽泰集團伏巖煤業有限公司，山西 陽城 048100）

·試驗研究·

中等變質煤孔隙特征隨揮發分變化規律

柴天興，靳樂軍

（山西陽城陽泰集團伏巖煤業有限公司，山西 陽城 048100）

在中等變質煤的低溫液氮吸附實驗和工業分析實驗的基礎上，分析了煤的總孔體積、總比表面積和孔隙率隨揮發分變化的關系。研究發現：對于中等變質程度的肥煤和氣煤，隨著揮發分變化范圍為28．02%～40．39%，煤的總孔體積變化范圍為0．002 44～0．010 356 mL／g，煤的孔隙總比表面積的變化范圍為0．534～3．965 m2／g，煤的孔隙率的變化范圍為1．24%～7．33%；總孔體積隨揮發分的變化呈冪函數關系，而總比表面積隨揮發分呈二次函數關系；總孔體積、總比表面積和孔隙率隨揮發分增加均會呈現增加的趨勢，并均在揮發分為28%左右時達到最小值。

中等變質煤；孔隙結構；揮發分；變化規律

煤是一種具有孔隙裂隙雙重結構的多孔介質，煤中的瓦斯80%～90%以上呈吸附狀態存在于孔隙結構表面，因此，對于煤孔隙結構的研究可以更好地了解瓦斯在煤層中的賦存情況，同時可以了解煤層中瓦斯運移和聚集的規律，對于煤與瓦斯的突出機理也有十分重要的意義。對煤的孔隙結構的研究當前主要是通過對孔隙率（真、假密度）、表面特性以及孔徑分布研究，用于研究物質孔隙結構表征的主要方法有：等溫線法（N2或CO2）、壓汞法、氦比重法、NMR旋轉-松弛測量法、Wood＇s合金侵入法、小角度X射線散射法、小角度中子散射法、可視區域電磁波譜發射法等［1］。但是每一種研究法由于其原理及分析儀器等技術條件的局限性，都很難完全描述出煤的真實結構特征，因此在今后的研究過程中應該綜合利用各種手段，盡量得出真實的煤的孔隙結構情況。

各國學者對煤的孔隙結構進行了很多研究：White等［2］研究結果表明，隨著煤階的提高，孔隙度略有下降趨勢；煤的孔體積與煤種、煤的粒度以及所采用的測量有關。Rivera-utrilla［3］研究結果表明，CO2測得的煤的吸附面積大大高于N2的，并且煤的變質程度越高，N2進入該孔隙內的擴散速率越低，煤的表面積、孔徑分布與煤階、顯微組分有關。幾部等［4］發現瓦斯涌出大的砂川煤礦、大島煤礦煤的比表面積要比瓦斯涌出少的幌內煤礦、太平洋煤礦幾乎要大兩倍。邱介山等［5］研究結果表明，隨著變質程度的從低到高變化煤的CO2表面積曲線和N2表面積曲線先減少后增加，在Cdaf=85%和Vdaf=30%左右出現最小值。劉常洪［6］研究結果表明，相同變質程度的煤，不論何種煤巖組分，微孔變化不大；孔容和總比表面積隨著變質程度的變化有著“U”型的規律。陳萍等［7］研究得出，煤中孔隙是從小至分子級（孔徑約0．86 nm）大至無上限的較連續的孔系統，其中孔形結構可分為：開放型透氣性孔；一端封閉的不透氣孔；細頸瓶形孔。煤中孔徑約小于3 nm的孔大多為一端封閉的孔；細頸瓶形孔的瓶頸直徑也在3 nm左右。

現在國內用于測定孔隙結構比較成熟的技術是壓汞法和等溫吸附法，其它方法還不夠成熟。而壓汞法主要研究7．5 nm以上的孔隙結構，對于比7．5 nm更小的孔壓汞法還不能進行描述，而煤中微孔是瓦斯吸附的主要場所之一，因此，壓汞法對于微孔的研究比較欠缺。而對于等溫吸附法，可以采用的吸附氣體有N2和CO2。雖然用CO2能測出較小的孔徑的孔隙，但是用CO2，有較高的沸點和臨界溫度，在測定比表面積時實驗條件及煤本身的復雜性質使得實驗確定中孔及微孔的表面積、孔體積及孔徑分布特性的氣體吸附法存在許多困難。而用N2做吸附劑時，不僅可以滿足對煤微孔研究的需要，而且實驗條件也較容易實現。同時鑒于其它方法還不是很成熟，所以，本次實驗采用N2做吸附劑對中等變質程度煤的孔隙結構進行研究。

1 實驗原理

在液氮溫度（77K），N2分子進入待測樣品中產生多層吸附。在樣品內部多個點上的力能夠達到平衡，而在樣品的表面則不同，有剩余的表面自由能。因此，當分子與樣品的表面接觸時，便為其表面所吸附。當向經過預處理的樣品管注入氮氣后，氮氣便會在煤樣表面吸附，隨著相對壓力的增高，吸附量逐漸增加。煤在吸附低溫氮氣的過程中，氮氣分子（分子直徑約為4．3×10-10m）首先在煤的微孔里發生毛細填充以及在較大孔壁上進行單分子層吸附，隨著相對壓力的升高，單分子層排滿，吸附層加厚，當相對壓力達到了壓力表設定的相對壓力的最高值時，毛細凝聚便達到了最大孔半徑，吸附與凝聚的過程便告結束。吸附的機理為微孔填充，填充的過程為在孔壁上一層又一層地筑膜［8］。低溫氮吸附過程的曲線圖見圖1。

圖1 低溫氮等溫吸附曲線圖

低溫氮吸附實驗主要是在測定N2（77K）吸附等溫線的基礎上，利用BET吸附理論測定比表面積；當BET線性較差時，可用αs來求中孔的比表面積，同時αs還可求微孔體積；由BJH法測定中孔孔徑分布及中孔孔容［9］。

本次低溫液氮實驗采用的是SSA-4200孔隙比表面積測定儀。它執行的標準是GB1T 19587-2004，主要輸出報告有吸附及脫附等溫線、BJH孔體積分布圖、BJH孔比表面積分布圖、總孔體積、BET比表面積、BJH累積比表面積；可以對孔隙微觀結構的比表面積、孔容、孔半徑進行全面的分析。儀器采用軟件操作可以對兩個煤樣同時進行分析，測試的孔徑為3．5×10-10～2 000×10-10m；測試的比表面積下限為0．000 1m2／g～無上限；壓力測定范圍為0～120 kPa，精度為：±0．1%FS。

2 煤樣的制備及實驗測定

實驗煤樣的8個煤樣主要來自兩淮礦區，并對取來的煤樣進行了編號分別為1＃、2＃、3＃、4＃、5＃、6＃、7＃和8＃。工業分析實驗采用中華人民共和國國家標準煤的工業分析方法（GB／T 212-2008），實驗煤樣粒徑為0．074～0．2 mm；為了排除煤樣粒徑對孔隙測定的影響，煤的低溫液氮吸附實驗也采用煤樣粒徑為0．074～0．2 mm進行實驗。工業分析實驗結果見表1。

表1 實驗煤樣工業分析表

由表1可知，本次實驗的煤樣主要是中等變質程度的煤，煤的揮發分變化范圍28．02%～40．39%。根據俞啟香［9］利用揮發分劃分煤種的方法，1＃、2＃、3＃、4＃煤樣的揮發分變化范圍為28．02%～30．57%，變質程度為肥煤；4＃、5＃、6＃、7＃、8＃煤樣的變化范圍為34．22%～40．39%，變質程度為氣煤。

3 結果與討論

由于煤中孔隙結構受到煤變質程度的影響，而揮發分可以作為煤變質程度的指標之一，因此，通過研究煤的孔隙結構與揮發分的變化關系，從而可以更好地認識煤的孔隙結構特征隨變質程度變化規律。中等變質程度肥煤和氣煤孔隙特征分布表見表2。

表2 中等變質程度肥煤和氣煤孔隙特征分布表

由表2可知，對于中等變質程度的肥煤和氣煤，隨著揮發分變化范圍為28．02%～40．39%，煤的總孔體積變化范圍為0．002 44～0．010 356 mL／g，煤的孔隙總比表面積的變化范圍為0．534～3．965 m2／g，煤的孔隙率的變化范圍為1．24%～7．33%。

為了更好地研究中等變質程度煤孔隙的總孔體積和總孔比表面積隨揮發分的變化情況，因此，采用不同的曲線進行數值擬合分析并進行比較，本次擬合采用的曲線主要有線性函數、二次函數、乘冪函數和指數函數。擬合的具體效果見表3。

表3 不同函數下總孔體積、總比表面積隨揮發分變化擬合效果表

由表3可知，對于中等變質程度的肥煤和氣煤，根據相關指數的大小可知，總孔體積隨揮發分的變化關系擬合效果為：乘冪函數＞指數函數＞線性函數＞二次函數；而總比表面積隨揮發分的變化關系擬合效果為：二次函數＞線性函數＞指數函數＞乘冪函數。因此對總孔體積隨揮發分的變化乘冪函數y= 0．000 000 7x2．5509535可以實現較好的擬合效果，而對于總比表面積隨揮發分的變化關系二次函數y= 0．020 1x2-1．121 5x+16．57可以實現較好的擬合效果。為了更好地進行分析，做出了中等變質煤總孔體積、總比表面積與揮發分變化關系最佳效果圖，見圖2，圖3。

圖2 中等變質煤總孔體積與揮發分變化關系圖

圖3 中等變質煤總比表面積與揮發分變化關系圖

從圖2和圖3可以得出，在中等變質程度的肥煤和氣煤部分（揮發分在28%～40%），隨著揮發分的從低到高，總孔體積會呈現增加的趨勢，并且當揮發分在28%左右時，總孔體積值最小值為0．003 442 mL／g；而對于中的變質程度煤的總比表面積，隨著揮發分的增加（28%～40%）總比表面積也呈現上升的趨勢，在揮發分為28%左右時值最小值為0．926 4 m2／g。總孔體積和總比表面積隨揮發分的變化關系的大體趨勢是一樣，但是曲線變化規律與擬合曲線卻不一樣，這主要是比表面積不僅受到總孔體積的影響，而且還受到孔徑分布的影響。

在研究完中等變質程度煤總孔體積、總比表面積與揮發分的關系后，再分析孔隙率與比揮發分的關系，根據表2中的實驗數據并采用線性方程進行擬合，具體見圖4。

圖4 中等變質煤孔隙率與揮發分變化關系圖

從圖4可以看出，中等變質程度的肥煤和氣煤揮發分從28%～40%的增加過程中，孔隙率呈現出增加的趨勢，并且揮發分大概在28%左右時孔隙率的值最小為1．774 4%，這和俞啟香［10］的研究成果基本相同。

4 結 論

基于中等變質程度煤的低溫液氮和工業分析的實驗結果，主要得到了以下實驗結論：

1）對于中等變質程度的肥煤和氣煤，隨著揮發分變化的范圍為28．02%～40．39%，煤的總孔體積變化范圍為0．002 44～0．010 356 mL／g，煤的孔隙總比表面積的變化范圍為0．534～3．965 m2／g，煤的孔隙率的變化范圍為1．24%～7．33%。

2）對于中等變質程度的肥煤和氣煤，總孔體積隨揮發分的變化乘冪函數y=0．000 000 7x2．5509535可以實現較好的擬合效果，而對于總比表面積隨揮發分的變化關系二次函數y=0．020 1x2-1．121 5x+ 16．57可以實現較好的擬合效果。

3）對于中等變質程度的肥煤和氣煤，隨著揮發分的從低到高，總孔體積、總比表面積和孔隙率均會呈現增加的趨勢，并均在揮發分為28%左右時達到最小值；總孔體積、總比表面積和孔隙率最小分別為0．003 442 mL／g、0．926 4 m2／g、1．774 4%。

［1］ 孟憲明．煤孔隙結構和煤對氣體吸附特性研究［D］．青島：山東科技大學，2004．

［2］ White，Clarkson CR．Adsorption potential theories to coal methane adsorption isotherms at elevated temperature and pressure［J］．Carbon，1997（35）：1689-1705．

［3］ Rivera-utrilla，J．E．Fitzgerald．Adsorption of methane，nitrogen，carbon dioxide and their mixtures on wet Tiffany coal［J］．Fule，2004（84）：2351 -2363．

［4］ 幾部等．關于煤的表面反應—吸附熱及比表面積測定［J］．日本北海道礦山學會會志，1967（8）：29-30．

［5］ 邱介山，郭樹才．中國一些煤的比表面積變化規律［J］．燃料化學學報，1991（10）：250-259．

［6］ 劉常洪．煤孔隙結構的實驗研究［J］．煤礦安全，1993（8）：1-4．

［7］ 陳 萍，唐修義．低溫氮吸附法與煤中微孔隙特征的研究［J］．煤炭學報，2001（10）：552-556．

［8］ 田英姿，陳克復．用壓汞法和氮吸附法測定孔徑分布及比表面積［J］．中國造紙，2004：（4）：21-23．

［9］ 近藤精一，石川達雄，安部郁夫．吸附科學（李國希譯）［M］．北京：化學工業出版社，2005：41-48．

［10］ 俞啟香．礦井瓦斯防治［M］．徐州：中國礦業大學出版社，1992：1-19．

Research on Change Law of Pore Structure Characteristic of Middling Metamorphism Coal along with the Volatile

Chai Tian-xing，Jin Le-jun

Based on the cryogenic liquid nitrogen adsorption experiment and industrial analysis experiment of the middling metamorphism coal，the relationship among total pore volume，the specific surface area and porosity ratio of coal along with the volatile law is further studied．And it is found that for the middling metamorphism fat coal and gas coal，with the volatile variation range from 28．02%to 40．39%，the total pore volume of the coal changes from 0．002 44 to 0．010 356 mL／g，the total surface area changes from 0．534 to 3．965 m2／g，and range of porosity ratio varies from 1．24%to 7．33%．The total pore volume changes along with the volatiles presented as exponent function，while the total surface area changes with the volatile as quadratic function．The total pore volume，total specific surface area and porosity ratio will show an increasing trend with volatile，and the minimum value appears when the volatile is 28%．

Middling metamorphism coal；Pore structure；Volatile；Change law

TD984

B

1672-0652（2013）09-0007-04

2013-06-07

柴天興（1966—），男，山西陽城人，2011年畢業于太原理工大學，工程師，主要從事煤礦安全生產技術管理工作（E-mail）jlejun＠163．com