王曉東 ,于建軍 ,劉海東 ,劉應科 ,劉向忠 ,龍昭熹
(1.西山煤電(集團)有限責任公司,山西 太原 030053;2.中國礦業大學 安全工程學院,江蘇 徐州 221116)
煤與瓦斯突出是一種嚴重威脅煤礦安全的動力災害[1]。它是在地應力與瓦斯壓力共同作用下,在很短的時間內將煤與瓦斯向采掘空間大量噴出的現象。煤體是一種典型的大分子有機質,能夠吸附大量瓦斯氣體,并能在突出災害孕育與發展過程中快速解吸與瞬間釋放。因此,煤體的瓦斯吸附、放散特性研究與參數測定對于突出災害危險的預測具有重要意義[2]。
國內外學者對煤體中瓦斯吸附與放散特性進行了較深入的研究[3-4]。BARRER[5]認為吸附和解吸是可逆過程,累積氣體吸附/解吸量與時間平方根成正比;張群等[6]采用平衡水分煤樣進行等溫吸附實驗,實驗結果表明煤的等溫吸附特性主要受煤級影響,隨著煤級增高,表征煤最大吸附能力的郎格繆爾體積逐漸增大;王恩元等[7]從實驗角度得出了Langmuir 吸附/解吸動力學模型;李一波等[8]研究表明瓦斯放散初速度隨粒徑呈對數函數變化規律;李志強等[9]通過擬合解吸實驗數據,提出了動擴散系數的數學表達式,建立了準確描述瓦斯擴散全過程的動擴散系數新模型;賈東旭等[10]研究了不同變質程度煤樣粒徑的改變對其瓦斯放散初速度的影響,煤樣的瓦斯放散初速度隨粒徑的增大呈冪函數關系減少;李成武等[11]研究了承壓煤體瓦斯解吸擴散特性,得出瓦斯解吸量、宏觀微觀有效擴散系數隨軸壓升高呈先下降后波動上升趨勢,且煤樣粒度越大,該系數對應力越敏感。
煤體在長期地質演化作用下,內部發育了豐富的孔隙及裂隙結構[12]。裂隙作為煤內瓦斯流動的主要通道,對瓦斯流動特性有著主要的控制作用,從而影響著煤體內部瓦斯的吸附、解吸及擴散等機制[13]。前人的研究主要集中在煤階、塊度、瓦斯壓力、水分、孔徑等因素對煤體瓦斯吸附、放散特性的影響,而對煤體裂隙結構及發育程度與瓦斯吸附、放散特性關系的研究較少[14]。同時,現場煤樣的參數測試主要包括瓦斯吸附常數與瓦斯放散初速度,通常將現場煤體取樣后轉移至實驗室進行分析,其測試周期長、操作繁瑣且工作量大[15]。而煤體裂隙結構特征則能夠快速識別與及時獲取,且更有利于實現智能化、數字化操作。
因此,現場選取原煤試樣,測試煤體裂隙結構與發育特征,并結合瓦斯吸附、放散實驗結果分析了不同煤樣表面裂隙結構對瓦斯吸附、放散特性的影響,揭示二者之間的對應關系。研究成果能夠為利用煤體裂隙結構特征快速評估瓦斯吸附、放散特性提供新思路,對于煤與瓦斯突出災害的精細預測與防治具有重要意義。
在山西省太原市東曲煤礦選取原煤試樣,試樣共分4 組,分別來源于山西組的2#煤層、4#煤層以及太原組的8#煤層、9#煤層(4 個煤層均為突出煤層),煤樣基本參數見表1。
煤樣表面裂隙結構特征是煤體內部裂隙分布及其發育程度在表面的外延展現。煤體內部裂隙分布難以直接獲取時,可通過對表面裂隙結構特征進行近似描述與表征。采用割機將煤樣表面切割平整,用自來水沖洗表面,得到煤樣表面平整化處理后的實物圖,4 組煤樣表面裂隙分布實物圖如圖1。
圖1 4 組煤樣表面裂隙分布實物圖Fig.1 Physical pictures of surface fracture distribution of four groups of coal samples
為使煤樣表面的裂縫與其他圖像信息進一步分離,對圖像信息進行增益降噪處理。通過提高亮度和對比度,調整灰度曲線,并適當地擦去干擾信息,補充表現不夠充分的裂隙信息,得到比較清晰的裂縫圖像,4 組煤樣圖片增益降噪結果如圖2。煤樣的表面裂隙分布呈現顯著的不規則性、非線性和隨機性,不同組煤樣之間差異較大。總體上,裂隙的分布較為分散,但通過細小孔隙及微小裂隙連接并貫通,這為瓦斯的吸附、解吸與擴散提供了便捷的通道[16]。
圖2 4 組煤樣圖片增益降噪結果Fig.2 Gain and noise reduction results of four groups of coal sample images
通過Matlab 軟件編譯數字圖像程序對圖2 中裂隙的空間分布進行求解,并轉化為矩陣形式[17]。圖像中的像素點以矩陣中的元素值代替,煤樣表面裂隙空間分布二值圖結果示例如圖3,效果十分顯著(“0”表示對應裂隙圖片中該位置含有裂隙,“1”表示該位置不含裂隙)。
圖3 煤樣表面裂隙空間分布二值圖結果示例Fig.3 Example of the results of binary image of spatial distribution of coal surface fractures
研究表明,煤巖體等材料的缺陷、裂隙結構的分布往往具有分形特征[18]。利用分形維數能夠從物理本質上表征煤體表面裂隙的結構特征。因此,采用盒子維數法求解煤體表面裂隙的分形維數[19]。
首先,選擇邊長為ε的正方形盒子,對煤體表面裂隙分布區域進行覆蓋(當煤樣裂隙圖片不為正方形時,進行適當裁剪),統計所需盒子的數量N;然后,改變盒子的邊長為εi,進一步統計該邊長對應的盒子數量為Ni,擬合兩者的關系合數,即可得到分形維度參數D:
基于式(1)對4 組煤樣二值圖對應的矩陣進行求解,擬合得到lgN-lgε關系曲線,對煤樣分形維數曲線求解,斜率即為D的值。裂隙分布的盒子尺寸和盒子總數之間存在良好的線性關系,表明所選煤樣表面裂隙結構具有顯著的分形特征,可以采用分形的方法對裂隙分布特征進行表征。一般地,分形維數越大說明煤體裂隙分布特征越復雜,裂隙發育程度越高。煤樣分形維數求解結果如圖4,圖中4 組煤樣的分形維數D范圍為1.4~1.5 之間,表明原生裂隙已經經過了充分的發育,且發育程度9#煤層>4#煤層>8#煤層>2#煤層。
圖4 煤樣分形維數求解結果Fig.4 Results of solving fractal dimension of coal samples
采用高壓容量法(MT/T 752—1997《煤的甲烷吸附量測定方法》)對4 組煤樣進行瓦斯等溫吸附實驗,并測試瓦斯壓力與瓦斯吸附量。按照AQ 1080—2009《煤的瓦斯放散初速度指標(Δp)測定方法》測定瓦斯放散初速度。實驗裝置為HCA-1型瓦斯吸附裝置與WT-1 型瓦斯擴散速度測定儀。
不同煤樣的等溫(30 ℃)吸附曲線如圖5。
圖5 不同煤樣的等溫(30 ℃)吸附曲線Fig.5 Isothermal (30 ℃) adsorption curves of different coal samples
結果表明,煤體對瓦斯的吸附量隨著壓力的增大而增大,初始階段煤樣瓦斯吸附速率較大;隨著吸附時間的增加,煤體內孔隙表面空間逐漸趨于飽和,瓦斯吸附速率逐漸降低。
基于Langmuir 吸附模型對不同組煤樣的瓦斯等溫吸附曲線進行擬合,得到煤樣吸附常數a、b值以及擬合度R2,各煤樣吸附常數統計結果見表2。各組煤樣擬合度均大于0.99,表明各煤樣瓦斯吸附過程均符合Langmuir 吸附理論,得到的吸附常數可靠。4 組煤樣a取值范圍為34.180~36.920 m3/t,且9#煤層>4#煤層>8#煤層>2#煤層,與分形維數大小順序一致。同一煤體吸附常數a較大時,b相應較小。
表2 各煤樣吸附常數統計結果Table 2 Statistical results of adsorption constants of coal samples
利用瓦斯放散初速度表征煤體的瓦斯放散特性。它與煤的微孔隙結構、孔隙表面性質和大小等因素有關,是鑒定煤層突出危險的關鍵指標。各煤樣瓦斯放散初速度見表3。
表3 各煤樣瓦斯放散初速度Table 3 Initial velocity of gas emission of each coal sample
4 組煤樣瓦斯放散初速度范圍為12~18 mmHg,且9#煤層>4#煤層>8#煤層>2#煤層,與分形維數值、吸附常數a值的大小順序一致。表明煤體的瓦斯放散特性與吸附特性是密切相關的,共同反映了瓦斯在煤體內部孔隙及裂隙中的吸附、運移與釋放的綜合能力。
煤體裂隙分布分形維數與瓦斯吸附和放散特性參數關系擬合結果如圖6。
圖6 煤體裂隙分布分形維數與瓦斯吸附和放散特性參數關系擬合結果Fig.6 Fitting results of relationship between fractal dimension of coal fracture distribution and gas adsorption and emission characteristics parameters
分析可知:分形維數D與瓦斯吸附常數a存在正相關性,擬合函數為a=35D-15,R2=0.672。分形維數反映了煤體裂隙發育程度,而吸附常數a反映了煤體基質表面飽和吸附瓦斯的量,即煤體吸附瓦斯的能力。結果表明,分形維數總體上能夠表征煤體的瓦斯吸附特性與吸附能力。
另外,分形維數與瓦斯放散速度同樣存在正相關性,擬合函數為Δp=71D-87,R2=0.821,擬合程度更高,離散性更小。瓦斯放散初速度衡量的是煤體初始暴露時瓦斯從吸附態轉化為游離態時的初速度,反映了煤體卸壓瞬間釋放瓦斯的能力。統計結果表明,相比吸附特性,分形維數能夠更好地表征瓦斯放散特性,它與突出煤體受采場擾動下瓦斯放散特性與釋放能力的關系更加密切。
在成煤過程中,構造作用會導致原生煤體發生變形與破裂,這是煤體裂隙結構改變與瓦斯運移通道發育的關鍵,進而顯著影響了煤體吸附與釋放瓦斯的能力。瓦斯吸附常數能夠反映原始煤層中煤體吸附瓦斯量的能力,是突出危險存在的“內因”;而瓦斯放散初速度能夠反映煤體在采場擾動下變形破裂后快速釋放瓦斯的能力,是突出危險致災的“外因”。煤體瓦斯的吸附、放散特性是煤與瓦斯突出災害預測的關鍵。
研究表明,煤體表面裂隙分布的分形維數指標反映的是煤體內部裂隙分布特征與發育程度,其值大小能夠表征煤體瓦斯吸附、放散的特性與能力。研究成果為利用煤體裂隙結構特征快速評估煤體瓦斯吸附、放散特性提供了新的研究思路。考慮到采煤現場難以對原煤試樣進行等溫吸附實驗與瓦斯放散實驗,瓦斯吸附、放散特性參數的測定不得不面臨著測試工序復雜、數據測試不及時、參數反饋周期長等諸多困難,而現場對原煤試樣可進行實時拍攝,實物圖片通過井下環網將數字化信息傳輸至地面工作站,可實現實時分析與智能決策,從而指導井下煤與瓦斯突出預警工作。該流程可視為無損測試、準實時反饋,這對于井下煤與瓦斯突出災害的精細預測與防治是十分重要的。
1)對井下原煤試樣的表面裂隙實物照片進行了數字化處理。結果表明:煤體表面裂隙分布信息呈現顯著的不規則性與非線性特征,能夠近似描述與表征煤體內部裂隙結構特征。
2)基于分形理論擬合求解了煤體表面裂隙分布的分形維數。結果表明:分形維數越大,煤體裂隙分布特征越復雜,裂隙發育程度越高;4 組煤樣的分形維數范圍為1.4~1.5,原生裂隙發育程度:9#煤層>4#煤層>8#煤層>2#煤層。
3)實驗測試分析了煤體瓦斯吸附、放散特性參數。結果表明:煤樣吸附常數擬合度均大于0.99,符合Langmuir 吸附理論,吸附常數a取值范圍為34.180~36.920 m3/t,瓦斯放散初速度為12~18 mmHg。上述參數的大小順序均為:9#煤層>4#煤層>8#煤層>2#煤層,與分形維數的大小順序一致。
4)研究分析了分形維數與瓦斯吸附、放散特性參數的相關性。結果表明:東曲煤礦煤體分形維數與瓦斯吸附、放散特性參數呈顯著的正相關性,其能夠表征該礦煤體瓦斯吸附與放散的特性與能力。