張憲尚

(1.瓦斯災(zāi)害監(jiān)控與應(yīng)急技術(shù)國家重點實驗室，重慶 400037；2.中煤科工集團重慶研究院有限公司，重慶 400037；3.四川大學，四川 成都 610065)

0 引言

研究煤屑瓦斯解吸擴散特征，對于探索煤礦井下瓦斯涌出、預測煤與瓦斯突出、煤層瓦斯賦存參數(shù)測定及煤層氣產(chǎn)能預測等方面具有重要實際意義[1-4]。國內(nèi)外學者針對不同條件煤屑瓦斯解吸特征進行大量試驗研究并總結(jié)得到不同的經(jīng)驗?zāi)Ｐ?，如Barrer[5]依據(jù)甲烷氣體分子在天然沸石中解吸規(guī)律試驗研究，提出表達氣體累計解吸量和時間關(guān)系的巴雷爾式；Airey[6]研究碎煤瓦斯解吸擴散規(guī)律，得到艾黎經(jīng)驗公式；Bolt等[7]對試驗研究不同變質(zhì)程度煤的瓦斯放散過程并總結(jié)出博特式；孫重旭等[8]認為累計瓦斯解吸量與時間成冪函數(shù)關(guān)系;王佑安等[9]將煤樣解吸環(huán)境快速降至真空，總結(jié)得到出王佑安經(jīng)驗表達公式;宋世釗[10]在譯著中介紹了前蘇聯(lián)的研究煤中瓦斯解吸量與時間關(guān)系的烏斯基諾夫式。

以上經(jīng)驗?zāi)Ｐ褪歉鶕?jù)不同試驗條件和分析角度得到的，存在一定適用范圍并可以在一些特定環(huán)境下相互轉(zhuǎn)化[11-12]。如巴雷爾公式在地勘期間硬煤取芯推算瓦斯損失量方面得到廣泛應(yīng)用，但卻不適于構(gòu)造煤的瓦斯損失量的推算[13]；艾黎式明顯低估了瓦斯解吸開始階段的放散量[11]。根據(jù)瓦斯放散量與時間關(guān)系及公式特點，這些經(jīng)驗?zāi)Ｐ痛笾驴煞譃閮绾瘮?shù)式和指數(shù)式2大類，冪函數(shù)則被認為煤樣暴露后60 s內(nèi)的瓦斯放散速度理想的公式[14]。在實際應(yīng)用過程中，多利用這些經(jīng)驗?zāi)Ｐ蛿M合瓦斯累計解吸量隨時間的變化數(shù)據(jù)并都表現(xiàn)出較高的擬合度，但對各個經(jīng)驗?zāi)Ｐ蛿M合后瓦斯累計解吸量的預測方面關(guān)注較少。

本文在試驗測定煤屑瓦斯累計解吸量的基礎(chǔ)上，取初期不同時間段瓦斯累計放散量數(shù)據(jù)擬合得到經(jīng)驗?zāi)Ｐ蛥?shù)，將模型參數(shù)代入到經(jīng)驗?zāi)Ｐ椭杏嬎愠霾煌瑫r間的瓦斯解吸量，并將其結(jié)果與試驗測定值進行對比。該研究對經(jīng)驗?zāi)Ｐ驮诿旱耐咚褂砍黾懊簩託猱a(chǎn)能預測方面有著重要實際意義。

1 煤屑瓦斯累計解吸量的試驗測定

1.1 試驗樣品

試驗樣品采取寧夏汝箕溝礦3號煤層335工作面煤樣，按標準破碎、篩分成粒度1～3 mm的質(zhì)量100 g樣品備用。樣品的密度、工業(yè)分析及瓦斯等溫吸附參數(shù)見表1。

表1 試驗樣品的工業(yè)分析及等溫吸附參數(shù)

1.2 試驗設(shè)備

煤屑瓦斯的解吸擴散屬于一種傳質(zhì)過程，該過程受煤屑周圍氣體壓力、環(huán)境溫度及煤屑內(nèi)部孔裂結(jié)構(gòu)特征等因素影響，其測定方法一般主要分為重量法和容量法2種，重量法是根據(jù)吸附劑(煤屑)解吸前后質(zhì)量變化測定得到解吸量，而容量法則根據(jù)解吸前后固定體積氣體壓力的變化測定解吸的大小[15]。試驗系統(tǒng)主要是由法國Setaram公司研制的全自動siverts型高壓氣體吸、解吸分析儀PCTPro-evo組成，試驗系統(tǒng)實物如圖1所示。

1-PCTPro-evo儀器；2-試驗樣品池;3-加熱套;4-熱電偶;5-數(shù)據(jù)采集及控制系統(tǒng);6-干式真空泵;7-甲烷氣體氣路;8-氮氣氣體氣路;9-氦氣氣體氣路；10-閥門；11-閥門。圖1 PCTPro-evo吸附解吸試驗系統(tǒng)Fig.1 PCTPro-evo adsorption and desorption test system

1.3 試驗步驟及結(jié)果

1)將待測煤樣取30 g左右，用高精度電子稱詳細稱取煤樣質(zhì)量并記錄，利用錫紙包裹裝入樣品池中并連接儀器，裝好加熱套等。

2)開啟真空泵，開啟吸附解吸測定儀，對儀器進行自檢及氣密性檢測。

3)氣密性檢測合格后，設(shè)定溫度達到60 ℃進行抽真空6 h以上。

4)設(shè)定吸附溫度為30 ℃，待儀器溫度穩(wěn)定后，進行自由空間校準。

5)設(shè)置吸附壓力，并設(shè)置吸附平衡判別標準。

6)樣品吸附平衡后，關(guān)閉閥門10，對儀器進行解吸測定設(shè)置，將參考罐中氣體壓力設(shè)定為大氣壓壓力。

7)開啟閥門11，放空2 min模擬鉆取煤樣的暴露過程，關(guān)閉閥門11，開啟閥門10，進行瓦斯解吸累計量測定，如圖2所示。

圖2 累計瓦斯解吸量Fig.2 Cumulative amount of gas desorption

8)瓦斯解吸測定完成后，關(guān)閉測試軟件及儀器設(shè)備。

2 經(jīng)驗?zāi)Ｐ蛿M合

2.1 經(jīng)驗?zāi)Ｐ突貧w擬合分析

對于含有瓦斯的煤屑解吸放散過程，在煤屑暴露瞬間，只是煤屑外面附近的瓦斯開始涌出。隨著時間的增長，煤屑深部的瓦斯也開始運移，瓦斯運移界面逐漸縮小，瓦斯流動通道長度增加，運動阻力增大，當瓦斯流場長度達到煤屑半徑后，煤屑中心的瓦斯壓力將相應(yīng)地降低。因此，瓦斯煤屑中的放散速度隨著時間的延長而降低，描述這一衰減過程的經(jīng)驗方程有[12]：巴雷爾式、撫順式、重慶式、博特式、艾黎式、文特式、烏斯季諾夫式、指數(shù)式、模擬近似式等。對于特定煤樣，在其吸附解吸環(huán)境一定的情況下，當時間趨于無窮時，其瓦斯解吸量是恒定的，因此在實際的擬合過程中，不同公式的形式可以轉(zhuǎn)化成相同的擬合形式，如文特式與重慶式。由于試驗測定是煤屑暴露后的瓦斯解吸量，因此需要考慮暴露期間的瓦斯損失量，故擬合公式中出現(xiàn)常數(shù)項，并選取暴露時間后解吸5，10，30，60 min的數(shù)據(jù)，利用origin軟件對其進行擬合得到不同經(jīng)驗公式的參數(shù)值，詳見表2。

2.2 回歸擬合過程及效果分析

利用origin軟件采用Levenberg Marquardt算法對瓦斯累計解吸量隨時間變化的數(shù)據(jù)進行回歸擬合分析，其經(jīng)驗公式回歸分析得到的各個參數(shù)受到迭代初始值的限制，因此對同一經(jīng)驗?zāi)Ｐ突貧w分析過程中需要設(shè)置初始值時，保持初始值一致或按照參數(shù)隨時間變化設(shè)置。

表2 各經(jīng)驗?zāi)Ｐ凸綌M合參數(shù)值

表2(續(xù))

從確定系數(shù)R-Square(詳見圖3)值大小可以看出[16]，不同公式擬合過程中確定系數(shù)都接近1，表明模型公式的變量對數(shù)據(jù)解釋能力較強，模型對數(shù)據(jù)擬合效果也較好，其擬合效果好劣程度依次為：巴雷爾式<博特-指數(shù)式<撫順式<模擬近似式<重慶-文特式<烏斯季諾夫式<艾黎式。隨著采取數(shù)據(jù)時間段的增加，重慶-文特式、烏斯季諾夫式及艾黎式的確定系數(shù)R-Square變化平穩(wěn)，擬合效果不會隨著數(shù)據(jù)多少而發(fā)生變化，而其他經(jīng)驗?zāi)Ｐ凸诫S著數(shù)據(jù)的增多而減低，其中巴雷爾式降低尤為明顯。

3 瓦斯解吸量預測對比

為了分析不同經(jīng)驗?zāi)Ｐ凸綄γ盒纪咚菇馕念A測能力，分別取不同時間段模型參數(shù)值代入經(jīng)驗公式得到不同時間的累計瓦斯解吸量，并與試驗測定值進行對比，如圖4所示。

從圖4中可以看出，巴雷爾式明顯高估擬合時間段后的累計瓦斯解吸量，不適合對累計瓦斯解吸量進行預測，而指數(shù)型經(jīng)驗?zāi)Ｐ颓€在擬合時間段后很快趨于平直，而低估擬合后累計瓦斯解吸量；對采用5 min和10min時間段的數(shù)據(jù)擬合參數(shù)的各個經(jīng)驗?zāi)Ｐ椭?，烏斯季諾夫式及重慶式都表現(xiàn)出較高的瓦斯解吸量預測能力，這與文獻[17]結(jié)論相一致，特別是5 min時間段數(shù)據(jù)擬合參數(shù)的烏斯季諾夫式，但從初期的累計瓦斯解吸量對比中可以看出烏斯季諾夫式明顯高估了煤屑瓦斯解吸量，期間預測能力表現(xiàn)不如重慶式；對30 min和60 min時間段的數(shù)據(jù)擬合參數(shù)的各個經(jīng)驗?zāi)Ｐ椭?，重慶式及艾黎式都表現(xiàn)出較好的瓦斯解吸量預測能力，特別是艾黎式在長期的累計瓦斯解吸量中有較好的預測能力。

圖3 不同時間段數(shù)據(jù)回歸分析的確定系數(shù)Fig.3 Determination coefficients for regression analysis of data in different time periods

圖4 經(jīng)驗?zāi)Ｐ凸酵咚菇馕款A測與對比Fig.4 Prediction and comparison of gas desorption amount based on empirical model formulas

為分析不同時間段數(shù)據(jù)量對經(jīng)驗?zāi)Ｐ偷念A測能力影響規(guī)律，特選博特-指數(shù)式、艾黎式、烏斯季諾夫式及重慶-文特式不同時間段擬合計算得到的曲線與實際試驗測定曲線進行對比，如圖5所示。

隨著瓦斯解吸數(shù)據(jù)時間段的增長，博特-指數(shù)式、艾黎式及重慶-文特式計算得到累計瓦斯解吸量曲線越來越接近于實際測定得到瓦斯解吸量曲線，而烏斯季諾夫式隨著瓦斯解吸量數(shù)據(jù)時間段增長而遠離實際測定曲線，其中博特-指數(shù)式和烏斯季諾夫式會隨著時間段延長呈現(xiàn)出規(guī)律性接近或遠離。隨著時間段的增長，對于各個經(jīng)驗?zāi)Ｐ陀嬎愕玫角€而言：博特-指數(shù)式在擬合時間段后很快趨于穩(wěn)定平直，低估煤屑長期瓦斯解吸量，這與指數(shù)衰減的擬合形式有關(guān)；艾黎式接近實際測定曲線速度最快，在長時間段瓦斯解吸數(shù)據(jù)預測累計解吸量表現(xiàn)最好；在利用較短時間煤屑瓦斯解吸量數(shù)據(jù)預測長期瓦斯解吸時，烏斯季諾夫式表現(xiàn)最好。

圖5 不同經(jīng)驗?zāi)Ｐ凸酵咚菇馕A測量隨時間段的變化Fig.5 Variation of gas desorption prediction amount with time periods by different empirical model formulas

4 結(jié)論

1)各個經(jīng)驗?zāi)Ｐ凸皆跀M合不同時間段的煤屑瓦斯累計解吸量時，確定系數(shù)R-Square都接近1，表現(xiàn)出較高擬合效果，其擬合效果好劣程度依次為：巴雷爾式<博特-指數(shù)式<撫順式<模擬近似式<重慶-文特式<烏斯季諾夫式<艾黎式。

2)巴雷爾式不適合用于煤屑的瓦斯解吸量預測，指數(shù)型經(jīng)驗?zāi)Ｐ凸接嬎愕玫酵咚菇馕壳€在擬合時間段后很快趨于平直，但會隨著時間段的增長而接近實際測定曲線，其受制于擬合數(shù)據(jù)時間段長短而低估累計瓦斯解吸量。

3)短時長煤屑瓦斯解吸數(shù)據(jù)推算預測長期瓦斯解吸量時適合采用烏斯季諾夫式，而預測短期時適合采用重慶-文特式；利用較長時間段瓦斯解吸量數(shù)據(jù)推算煤屑瓦斯解吸量宜采用艾黎式。