韓　穎，王　博，張飛燕

(1.中原經(jīng)濟(jì)區(qū)煤層(頁巖)氣河南省協(xié)同創(chuàng)新中心，河南 焦作 454000；2.河南理工大學(xué)能源科學(xué)與工程學(xué)院，河南 焦作 454000；3.河南省瓦斯地質(zhì)與瓦斯治理重點實驗室—省部共建國家重點實驗室培育基地，河南 焦作 454000；4.河南理工大學(xué)安全科學(xué)與工程學(xué)院，河南 焦作 454000)

基于啟動壓力梯度的煤層瓦斯流動狀態(tài)快速判識方法

韓穎1，2，3，王博2，張飛燕3，4

(1.中原經(jīng)濟(jì)區(qū)煤層(頁巖)氣河南省協(xié)同創(chuàng)新中心，河南 焦作 454000；2.河南理工大學(xué)能源科學(xué)與工程學(xué)院，河南 焦作 454000；3.河南省瓦斯地質(zhì)與瓦斯治理重點實驗室—省部共建國家重點實驗室培育基地，河南 焦作 454000；4.河南理工大學(xué)安全科學(xué)與工程學(xué)院，河南 焦作 454000)

摘要：為探討低滲煤層瓦斯流動的非達(dá)西特征，基于GSI煤體分類體系與砂巖啟動壓力梯度測試方法，進(jìn)行了不同煤體結(jié)構(gòu)煤樣的啟動壓力梯度測試，建立了啟動壓力梯度與煤體結(jié)構(gòu)的回歸關(guān)系，實現(xiàn)了煤層瓦斯流動狀態(tài)的快速判識。研究表明：在煤礦井下采集煤樣或直接觀測煤壁，獲取其煤體結(jié)構(gòu)，并與GSI煤體分類體系進(jìn)行比對，確定煤樣GSI值后，通過啟動壓力梯度測試并建立其與GSI的關(guān)系，即可實現(xiàn)不同煤體結(jié)構(gòu)煤層啟動壓力梯度的快速獲取與瓦斯流動狀態(tài)的快速判識。

關(guān)鍵詞：啟動壓力梯度；煤層瓦斯；流動狀態(tài)；快速判識；地質(zhì)強(qiáng)度指標(biāo)；低速非線性滲流

長期的低滲透油氣藏開發(fā)實踐表明，流體在低滲多孔介質(zhì)中的流動為帶有啟動壓力梯度的低速非線性滲流[1]。啟動壓力梯度概念由B.A.費勞林(蘇聯(lián)，1951)提出，他認(rèn)為只有當(dāng)實際壓力梯度大于某一臨界值時，流動才能發(fā)生，此臨界值稱為啟動壓力梯度[2]。閆慶來等[3]、吳景春等[4]通過室內(nèi)實驗證實了流體在低滲儲層內(nèi)滲流時存在非線性段及啟動壓力梯度；陳永敏等[5]通過實驗論證了存在滲流啟動壓力和低速滲流時出現(xiàn)非線性的低速非達(dá)西滲流規(guī)律。多年來，眾多研究者就低速非線性滲流的形成機(jī)理[6]、判識標(biāo)準(zhǔn)[7]、數(shù)值解法[8]等問題開展了大量研究，但其對象皆為低滲透油氣藏。

我國煤層滲透率一般在1 md以下，具有“低滲”特點，郭紅玉等[9]通過實驗初步證實了氣體通過低滲煤樣時存在低速非達(dá)西現(xiàn)象和啟動壓力梯度。但是，現(xiàn)行瓦斯抽采技術(shù)大多遵循線性滲流理論—達(dá)西定律，往往忽略了低滲煤層內(nèi)存在的低速非線性滲流及擴(kuò)散現(xiàn)象，導(dǎo)致抽采難易程度評價及抽采工藝選擇出現(xiàn)偏差。針對具有不同煤體結(jié)構(gòu)的煤層，如何快速、準(zhǔn)確地評價其抽采難易程度以及選擇適宜的抽采工藝，以實現(xiàn)瓦斯高效抽采，是煤礦現(xiàn)場亟待解決的問題；而解決該問題的前提及關(guān)鍵，在于對瓦斯在不同煤體結(jié)構(gòu)煤層內(nèi)流動狀態(tài)的準(zhǔn)確把握?；诖?，本文開展了基于啟動壓力梯度的煤層瓦斯流動狀態(tài)快速判識方法研究。

1基于地質(zhì)強(qiáng)度指標(biāo)的煤體結(jié)構(gòu)定量表征方法

1995年，E.Hoek等創(chuàng)立了地質(zhì)強(qiáng)度指標(biāo)[10](Geological Strength Index，GSI)，這是一種巖體分類體系，根據(jù)巖體結(jié)構(gòu)、巖體中巖塊的嵌鎖狀態(tài)與不連續(xù)面質(zhì)量，綜合各種地質(zhì)信息，估算不同地質(zhì)條件下的巖體強(qiáng)度。GSI巖體分類體系[10]如圖1所示。

圖1　GSI巖體分類體系

郭紅玉等[11]分別采用煤體被切割的基質(zhì)塊與裂隙寬度代替圖1中的巖體塊度與不連續(xù)面風(fēng)化狀況，建立了基于GSI的煤體結(jié)構(gòu)定量表征方法。與傳統(tǒng)煤體結(jié)構(gòu)五分法相對應(yīng)的GSI煤體分類體系[11]如圖2所示，GSI取值范圍見表1。

2煤樣啟動壓力梯度測試方法

2.1測試原理

本文測試借鑒吳凡等[12]計算砂巖啟動壓力梯度的方法，其基本思路如下所示。

根據(jù)達(dá)西定律，當(dāng)不考慮啟動壓力梯度時的氣體滲流方程見式(1)。

圖2　GSI煤體分類體系

表1　GSI取值范圍

(1)

式中，v為氣體通過煤芯的流速，m/s；K為煤芯滲透率，m2；p1為入口氣體壓力，Pa；p2為出口氣體壓力，Pa；p0為大氣壓力，101325 Pa；μ為氣體動力粘度，Pa.s；L為氣體流經(jīng)長度，即煤芯高度，m。

當(dāng)存在啟動壓力梯度時，氣體滲流方程見式(2)。

(2)

令v=0，則推導(dǎo)得式(3)。

(3)

因此，啟動壓力梯度計算公式見式(4)。

(4)

式中，λ為啟動壓力梯度，Pa/m。

2.2測試系統(tǒng)與煤樣制備

啟動壓力梯度測試系統(tǒng)主要由高壓N2氣源、RMT-150B巖石力學(xué)試驗機(jī)、樣品室、CY-60型氣體壓力傳感器、氣體質(zhì)量流量計、氣體流量積算儀、數(shù)據(jù)采集及處理儀器構(gòu)成[9]，如圖3所示。

圖3　啟動壓力梯度測試系統(tǒng)

數(shù)據(jù)采集儀采用YJZ-16型靜態(tài)電阻應(yīng)變儀，可與計算機(jī)通訊，實現(xiàn)數(shù)據(jù)實時、自動采集。氣體質(zhì)量流量計為兩種：D07-11CM型流量計，量程為20L/min，與D08-8CM型流量積算儀配套連接；D07-11C型流量計，量程為50mL/min，與D08-8C型流量積算儀配套連接，兩者串聯(lián)工作，實現(xiàn)氣體流量的實時、全程測試。

由于樣品室較小，僅能容納Ф50mm×50mm的煤樣，故將Ⅰ、Ⅱ類煤體結(jié)構(gòu)煤樣直接制備成Ф50mm×50mm的原煤煤芯；對于Ⅲ～Ⅴ類煤體結(jié)構(gòu)煤樣，因其較為破碎、無法直接鉆取煤芯，故依據(jù)密度相等的原則，稱取與Ⅰ、Ⅱ類煤芯相同質(zhì)量的煤樣，加壓制備成Ф50mm×50mm的型煤煤芯。

2.3測試方法與數(shù)據(jù)處理

將制備好的煤芯裝入樣品室，并在其上方放置透氣板；通過RMT-150B試驗機(jī)向煤芯加載4kN軸向力及2MPa圍壓，開啟高壓N2源，使氣體流經(jīng)煤芯后，直接排至大氣；氣體壓力傳感器及氣體質(zhì)量流量計同步、實時測試氣體壓力與流量。為確保系統(tǒng)密封良好，測試過程中氣體壓力不得超過2MPa。

表2　沙曲煤礦1#煤芯測試數(shù)據(jù)

3基于啟動壓力梯度的煤層瓦斯流動狀態(tài)快速判識方法研究

研究表明：當(dāng)瓦斯在煤層中流動時，若瓦斯壓力梯度Δp/L(煤層埋藏深度每增加1m，煤層瓦斯壓力的平均增加值，一般為0.005～0.015MPa/m[13])大于等于啟動壓力梯度，瓦斯將發(fā)生低速非線性滲流；反之則不流動，瓦斯僅通過擴(kuò)散途徑產(chǎn)出；只有當(dāng)啟動壓力梯度為零時，瓦斯才會發(fā)生線性滲流[14]。

根據(jù)前述測試方法，對采自華晉焦煤有限責(zé)任公司沙曲煤礦的一系列煤樣進(jìn)行啟動壓力梯度測試，并參照GSI煤體分類體系(圖2)，對不同煤體結(jié)構(gòu)煤樣的GSI賦值，結(jié)果見表3。

根據(jù)表3中數(shù)據(jù)，對λ與GSI關(guān)系進(jìn)行回歸分析，如圖5所示。可以看出，λ與GSI關(guān)系顯著。當(dāng)GSI=46.58時，λ=0；當(dāng)GSI<46.58時，λ=0.00189(46.58-GSI)1.49665，λ隨GSI的減小而增大，即煤體結(jié)構(gòu)越破碎，瓦斯流動所需的啟動壓力梯度越大；當(dāng)GSI>46.58時，λ=0.00189(GSI-46.58)1.50559，λ隨GSI的增大而增大，即煤體結(jié)構(gòu)越完整，瓦斯流動所需的啟動壓力梯度越大。

圖4　v與關(guān)系回歸分析圖

表3　不同煤體結(jié)構(gòu)煤樣的GSI值與啟動壓力梯度測試數(shù)據(jù)

圖5　λ與GSI關(guān)系回歸分析圖

綜上所述，對于沙曲煤礦而言，只要在井下采集煤樣或直接觀測煤壁，獲取其煤體結(jié)構(gòu)，并與GSI煤體分類體系進(jìn)行比對，確定煤樣GSI值后，即可通過上述關(guān)系得出λ，將λ與Δp/L進(jìn)行比較，即可快速判定煤層瓦斯的流動狀態(tài)：Δp/L≥λ，為低速非線性滲流；Δp/L<λ，為擴(kuò)散；λ=0，為線性滲流。

4結(jié)論

1)基于GSI煤體分類體系與砂巖啟動壓力梯度測試方法，進(jìn)行了不同煤體結(jié)構(gòu)煤樣的啟動壓力梯度測試，并建立了啟動壓力梯度與GSI的回歸關(guān)系。

2)在煤礦井下采集煤樣或直接觀測煤壁，獲取其煤體結(jié)構(gòu)，并與GSI煤體分類體系進(jìn)行比對，確定煤樣GSI值后，通過啟動壓力梯度測試并建立其與GSI的關(guān)系，即可實現(xiàn)不同煤體結(jié)構(gòu)煤層啟動壓力梯度的快速獲取與瓦斯流動狀態(tài)的快速判識。

需要說明的是，本文僅為初步實驗，所得結(jié)論有待進(jìn)一步驗證。下一步需增加實驗樣本，并對啟動壓力梯度測試系統(tǒng)進(jìn)行改進(jìn)，擴(kuò)大充氣壓力范圍；此外，圍壓、氣體黏度及煤樣滲透率、含水飽和度、溫度等因素對啟動壓力梯度的影響需進(jìn)一步深入研究。

參考文獻(xiàn)

[1]黃延章，等.低滲透油層滲流機(jī)理[M].北京：石油工業(yè)出版社，1998.

[2]戈爾布諾夫.A.T.異常油田開發(fā)[M].張樹寶，譯.北京：石油工業(yè)出版社，1987.

[3]閆慶來，何秋軒，尉立崗，等.低滲透油層中單相液體滲流特征的實驗研究[J].西安石油學(xué)院學(xué)報，1990，5(2)：1-6.

[4]吳景春，袁滿，張繼成，等.大慶東部低滲透油藏單相流體低速非達(dá)西滲流特征[J].大慶石油學(xué)院學(xué)報，1999，23(2)：82-84.

[5]陳永敏，周娟，劉文香，等.低速非達(dá)西滲流現(xiàn)象的實驗論證[J].重慶大學(xué)學(xué)報：自然科學(xué)版，2000，23(S1)：59-61.

[6]高樹生，熊偉，劉先貴，等.低滲透砂巖氣藏氣體滲流機(jī)理實驗研究現(xiàn)狀及新認(rèn)識[J].天然氣工業(yè)，2010，30(1)：52-55.

[7]任曉娟，張國輝，繆飛飛.低滲多孔介質(zhì)非達(dá)西滲流啟動壓力梯度存在判識[J].遼寧工程技術(shù)大學(xué)學(xué)報：自然科學(xué)版，2009，28(S1)：273-276.

[8]朱維耀，劉今子，宋洪慶，等.低/特低滲透油藏非達(dá)西滲流有效動用計算方法[J].石油學(xué)報，2010，31(3)：452-457.

[9]郭紅玉，蘇現(xiàn)波.煤儲層啟動壓力梯度的實驗測定及意義[J].天然氣工業(yè)，2010，30(6)：52-54.

[10]Hoek E.，Kaiser P.K.and Bawden W.F.Support of Underground Excavations in Hard Rock[M].Rotterdam，Balkema，1995.

[11]郭紅玉，蘇現(xiàn)波，夏大平，等.煤儲層滲透率與地質(zhì)強(qiáng)度指標(biāo)的關(guān)系研究及意義[J].煤炭學(xué)報，2010，35(8)：1319-1322.

[12]吳凡，孫黎娟，喬國安，等.氣體滲流特征及啟動壓力規(guī)律的研究[J].天然氣工業(yè)，2001，21(1)：82-84.

[13]于不凡，王佑安.煤礦瓦斯災(zāi)害防治及利用技術(shù)手冊(修訂版)[M].北京：煤炭工業(yè)出版社，2005.

[14]馬耕，蘇現(xiàn)波，魏慶喜.基于瓦斯流態(tài)的抽放半徑確定方法[J].煤炭學(xué)報，2009，34(4)：501-504.

Fast recognition method of gas flow state in coal seam based on starting pressure gradient

HAN Ying1，2，3，WANG Bo2，ZHANG Fei-yan3，4

(1.Collaborative Innovation Center of Coalbed Methane and Shale Gas for Central Plains Economic Region(Henan Province)，Jiaozuo 454000，China；2.School of Energy Science and Engineering，Henan Polytechnic University，Jiaozuo 454000，China；3.State Key Laboratory Cultivation Base for Gas Geology and Gas Control(Henan Polytechnic University)，Jiaozuo 454000，China；4.School of Safety Science and Engineering，Henan Polytechnic University，Jiaozuo 454000，China)

Abstract:In order to discuss non-Darcy characteristics of gas flow in low-permeability coal seam，based on coal body classification system and test method of starting pressure gradient in sandstone，starting pressure gradient in coal seam which have different coal structure was measured，the regression relationship between starting pressure gradient and coal structure was acquired，fast recognition of gas flow state in coal seam was realized.Research results show that coal structure can be acquired according to sampling or wall observation in coal mine firstly，GSI value of coal samples will be obtained by use of coal body classification system based on GSI，secondly starting pressure gradient can be tested in laboratory and the relationship between it and GSI will be acquired，then fast acquisition of starting pressure gradient and fast recognition of gas flow state in coal seam which have different coal structure will be realized.

Key words：starting pressure gradient；coal seam gas；flow state；fast recognition；geological strength index；low velocity nonlinear percolation

收稿日期：2015-07-20

基金項目：國家自然科學(xué)基金項目資助(編號：51404093)；河南省瓦斯地質(zhì)與瓦斯治理重點實驗室—省部共建國家重點實驗室培育基地開放基金項目資助(編號：WS2012A09)；國家安全生產(chǎn)監(jiān)督管理總局安全生產(chǎn)重大事故防治關(guān)鍵技術(shù)科技項目資助(編號：henan-0025-2015AQ；henan-0007-2015AQ)；河南省教育廳科學(xué)技術(shù)研究重點項目資助(編號：14A440012)；焦作市科技計劃項目資助(編號：2014400013)；河南理工大學(xué)博士基金項目資助(編號：B2008-15；B2012-093)

作者簡介：韓穎(1980-)，男，山東濟(jì)南人，博士，副教授，主要從事礦井瓦斯防治與瓦斯地質(zhì)領(lǐng)域的教學(xué)與科研工作。

中圖分類號：TD712+.5

文獻(xiàn)標(biāo)識碼：A

文章編號：1004-4051(2016)06-0104-05