姚冠霖,唐明云,2,鄭鵬先,段三壯,張海路
(1.安徽理工大學(xué) 能源與安全學(xué)院,安徽 淮南232001;2.煤礦安全高效開采省部共建教育部重點(diǎn)實(shí)驗(yàn)室,安徽 淮南232001)
滲透率作為反映瓦斯在煤巖中流動(dòng)難易程度的重要指標(biāo),對(duì)礦井瓦斯災(zāi)害的防治有著重要意義。研究表明,煤巖滲透率受溫度、圍巖應(yīng)力、孔隙壓力、流體屬性等因素的影響。溫度主要通過熱效應(yīng),膨脹效應(yīng)等方式影響煤巖滲透率。文獻(xiàn)[1-6]通過不同溫度下實(shí)驗(yàn)研究,得出滲透率隨溫度變化多呈現(xiàn)出指數(shù)函數(shù)關(guān)系以及冪函數(shù)關(guān)系,且包括煤階在內(nèi)的諸多條件均會(huì)影響煤巖滲透特性;而孔隙壓力主要通過改變煤巖吸附特性影響滲透率,文獻(xiàn)[7-10]中均通過實(shí)驗(yàn)均得出了煤巖滲透率隨孔隙壓力變化呈現(xiàn)出冪函數(shù)變化規(guī)律,以上的實(shí)驗(yàn)表明煤巖滲透率對(duì)溫度和孔隙壓力變化都表現(xiàn)出一定的敏感性,文獻(xiàn)[11-14]研究了不同溫度,壓力條件下的滲透率模型,不僅證實(shí)了考慮溫度與孔隙壓力的必要性,同時(shí)也在理論層面驗(yàn)證了前人實(shí)驗(yàn)得出的部分煤巖滲透率變化規(guī)律。以上學(xué)者在對(duì)煤巖滲透率研究方面取得了較好的成果,但其主要是對(duì)溫度或孔隙壓力等單因素作用下的煤巖滲透特性規(guī)律進(jìn)行研究,而且實(shí)驗(yàn)孔隙壓力和溫度普遍偏低,難以準(zhǔn)確反映高溫高壓下煤巖滲透特性,此外,實(shí)驗(yàn)多以甲烷為流體介質(zhì)進(jìn)行實(shí)驗(yàn),但文獻(xiàn)[15-16]都通過實(shí)驗(yàn)證實(shí)吸附氣體對(duì)煤巖滲透特性均存在一定影響。隨著傳統(tǒng)的壓裂技術(shù)越來越難以適應(yīng)生產(chǎn)需求,液氮壓裂作為一種新型的增透技術(shù)正逐步獲得應(yīng)用,可氮?dú)饩烤箤?duì)煤巖滲透率造成怎樣影響卻少有研究[17]。因此研究氮?dú)庾饔孟旅簬r滲透特性對(duì)實(shí)際生產(chǎn)具有一定指導(dǎo)意義。
實(shí)驗(yàn)采集煤巖取自潘集礦區(qū)11-2 煤層,使用巖石鉆孔取樣機(jī)和巖石磨片機(jī)取樣,將煤巖加工為實(shí)驗(yàn)試件,煤巖為直徑50 mm、長度100 mm 的圓柱體。
實(shí)驗(yàn)采用自主研制的RLQZ-2 型煤巖滲透率測試裝置,裝置可進(jìn)行高溫蒸汽或常溫下不同氣體的煤巖滲透率實(shí)驗(yàn)。實(shí)驗(yàn)裝置主要由氣瓶,真空機(jī),夾持器,圍壓泵,恒溫箱,流量計(jì)等組成,RLQZ-2 型煤巖滲透率測試裝置示意圖如圖1。待實(shí)驗(yàn)儀器參數(shù)表盤顯示讀數(shù)穩(wěn)定后,全部過程均在實(shí)驗(yàn)操作桌面完成,由計(jì)算機(jī)采集數(shù)據(jù)及軟件程序處理數(shù)據(jù),確保實(shí)驗(yàn)數(shù)據(jù)準(zhǔn)確可靠性。
圖1 RLQZ-2 型煤巖滲透率測試裝置示意圖Fig.1 Schematic diagram of RLQZ-2 coal rock permeability test device
根據(jù)中華人民共和國石油天然氣行業(yè)標(biāo)準(zhǔn)發(fā)布的SY/T 5336—1996《巖心常規(guī)分析方法》,采用基于達(dá)西定律的穩(wěn)定流法計(jì)算煤巖滲透率,即根據(jù)氣體通過煤巖的穩(wěn)定滲流量和煤巖兩端的滲透壓力差等可測量參數(shù)計(jì)算煤巖的滲透率,具體計(jì)算公式為:
式中:k 為煤巖試件實(shí)測滲透率,m2;Q 為標(biāo)況下的滲流流量,cm3/s;p0為大氣壓力,0.101 MPa;μ為氮?dú)鈿怏w動(dòng)力黏度,Pa·s;L 為試件長度,mm;A為試件橫截面面積,cm2;p1為進(jìn)口端壓力,MPa;p2為出口端壓力,MPa。
不同溫度氮?dú)鈿怏w動(dòng)力黏度不同,滲透率求解時(shí)需要考慮到溫度對(duì)氮?dú)怵ざ鹊挠绊?,不同溫度下氮?dú)鈩?dòng)力黏度見表1。
表1 不同溫度下氮?dú)鈩?dòng)力黏度Table 1 Nitrogen dynamic viscosity at different temperatures
溫度、應(yīng)力狀態(tài)以及天然裂隙等在內(nèi)的各種因素都會(huì)通過影響孔隙率,特別是有效孔隙率來影響滲透率,根據(jù)力學(xué)平衡原理,在孔隙率的研究成果的基礎(chǔ)上有效應(yīng)力方程如下:
式中:σ′為有效應(yīng)力,MPa;α 為孔隙壓縮系數(shù),取1;σi為外應(yīng)力,MPa;p 為孔隙壓力,MPa。
實(shí)驗(yàn)開始時(shí),首先將原煤試件安裝于夾持器中,對(duì)整個(gè)實(shí)驗(yàn)裝置進(jìn)行1 h 抽真空處理,以圍壓泵分別施以5、6、7、8、9 MPa 圍壓,設(shè)定實(shí)驗(yàn)溫度20 ℃,使用恒溫箱對(duì)試件進(jìn)行4 h 的恒溫后進(jìn)行2、3、4、5、6 MPa 孔隙壓力下的滲透率實(shí)驗(yàn)。實(shí)驗(yàn)結(jié)束后更換煤樣,再分別進(jìn)行實(shí)驗(yàn)溫度為40、60、80 ℃下的滲透率實(shí)驗(yàn)。
眾多學(xué)者所進(jìn)行的煤巖滲透特性實(shí)驗(yàn)中,施加的孔隙壓力普遍較低(3 MPa 以內(nèi)),對(duì)高壓狀態(tài)下煤巖滲透特性研究較少,因此進(jìn)行了孔隙壓力最高為6 MPa 的煤巖滲透率實(shí)驗(yàn),根據(jù)恒定有效應(yīng)力下改變孔隙壓力測得的氮?dú)夥€(wěn)定流量,按式(1)計(jì)算出煤巖滲透率,使用數(shù)據(jù)分析軟件對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合,不同溫度下滲透率隨孔隙壓力變化擬合曲線如圖2。
圖2 不同溫度下滲透率隨孔隙壓力變化擬合曲線Fig.2 Fitting curves of permeability with pore pressure at different temperatures
由圖2 可知,在恒定有效應(yīng)力條件下,改變實(shí)驗(yàn)溫度,煤巖滲透率隨孔隙壓力變化呈現(xiàn)出以下規(guī)律:
1)煤巖滲透率在不同溫度下隨孔隙壓力的變化趨勢基本一致,均呈現(xiàn)先快速減小而后逐漸平穩(wěn)的特征。
2)孔隙壓力一定時(shí),煤巖滲透率在不同溫度下表現(xiàn)出較大差異,實(shí)驗(yàn)開始時(shí)的煤巖滲透率隨溫度增加呈下降的趨勢,并且下降量隨溫度的升高而降低。當(dāng)孔隙壓力為2 MPa,與20 ℃低溫區(qū)間相比,當(dāng)溫度分別為40、60、80 ℃時(shí)的煤巖滲透率要低得多,降幅最高為35%,且主要集中在1.7×10-17~2.1×10-17m2這一區(qū)間內(nèi)。
3)孔隙壓力一定時(shí),不同溫度下滲透率在低壓階段表現(xiàn)出較大差異,40 ℃之后的煤巖滲透率曲線出現(xiàn)“斷崖式”下降,且溫度分別為60、80 ℃時(shí)的滲透率曲線比較接近,二者數(shù)據(jù)相差最大僅為6.95%。
在恒定有效應(yīng)力下,改變溫度條件,滲透率隨孔隙壓力的加載呈現(xiàn)出一定的規(guī)律性。在低孔隙壓力階段,煤巖的裂隙-孔隙結(jié)構(gòu)受到的壓縮作用較強(qiáng);同時(shí)因氮?dú)饩哂幸欢ǖ奈叫裕鶕?jù)表面化學(xué)基本理論,煤基質(zhì)吸附氣體后表面化學(xué)能降低,分子間作用力減小引起煤基質(zhì)膨脹,壓縮煤巖滲流通道;附著在基質(zhì)表面的氮?dú)庖矔?huì)摻混其他雜質(zhì)阻塞煤巖內(nèi)部孔隙。以上幾點(diǎn)均會(huì)降低煤巖的孔隙度造成滲透率快速下降。
隨著孔隙壓力的增大,氮?dú)獾奈节呌谄胶?,而煤體內(nèi)部孔隙體積是有限的,它的吸附變形會(huì)越來越小,即孔隙率的減小量會(huì)越來越小,滲透率的變化量也會(huì)越來越?。幻簬r又以微小孔最為發(fā)育,當(dāng)孔隙壓力達(dá)到一定值,煤巖內(nèi)部的原先閉合的微小孔打開,形成新的滲流通道。但孔隙壓力對(duì)滲透率的減小作用更強(qiáng),因此滲透率曲線仍然呈下降態(tài)勢,只不過下降趨勢放緩。
對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行非線性曲線擬合之后,可以得到恒定有效應(yīng)力下,孔隙壓力對(duì)滲透率影響方程:
式中:a1、b1為曲線擬合系數(shù)。
各溫度下曲線的擬合系數(shù)見表2。
根據(jù)不同溫度下測得的氮?dú)夥€(wěn)定流量,按式(1)可計(jì)算出煤巖滲透率,使用origin 對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合,不同孔隙壓力下滲透率隨溫度變化擬合曲線如圖3。
圖3 不同孔隙壓力下滲透率隨溫度變化擬合曲線Fig.3 Fitting curves of permeability with temperature under different pore pressures
分析圖3 可知,煤巖受恒定有效應(yīng)力時(shí),在不同孔隙壓力條件下,滲透率變化呈現(xiàn)如下規(guī)律:
1)煤巖滲透率在不同孔隙壓力下隨溫度變化的趨勢基本一致,均是呈現(xiàn)先快速降低后逐漸平緩的負(fù)指數(shù)關(guān)系。在20~40 ℃的溫度下,滲透率下降趨勢明顯,在40~80 ℃的溫度下,滲透率下降速度降低。
2)當(dāng)溫度一定時(shí),滲透率在隨孔隙壓力增大逐漸降低。在低溫低壓區(qū)間集中在2.4×10-17~2.6×10-17m2之間,而在在高溫高壓區(qū)間,滲透率集中在1.7×10-17~2.2×10-17m2。但總體上滲透率隨孔隙壓力的增加呈下降趨勢。
3)孔隙壓力一定時(shí),煤巖滲透率對(duì)溫度表現(xiàn)出不同的敏感性。即在不同孔隙壓力條件下升高溫度,在相同溫差范圍內(nèi),滲透率減小量越來越??;在溫差相同時(shí),孔隙壓力越大,滲透率減小量也越小。
根據(jù)實(shí)驗(yàn)可知,煤巖的于裂隙-孔隙結(jié)構(gòu)使其具有更強(qiáng)非均質(zhì)性,溫度引發(fā)的熱效應(yīng)通過影響煤巖的物理性質(zhì),使其更容易受應(yīng)力作用產(chǎn)生壓縮變形,很大程度上降低了煤巖孔隙度;由表1 可知,隨著溫度升高,氮?dú)鈩?dòng)力黏度增大,氮?dú)庠诿簬r中滲流變得更加困難,因此在實(shí)驗(yàn)前期滲透率快速下降;但煤巖內(nèi)部的孔隙體積是有限的,受有效應(yīng)力的壓縮作用,孔隙體積的改變量會(huì)越來越小,另一方面溫度升高導(dǎo)致煤巖內(nèi)部的水分等易揮發(fā)物質(zhì)在孔隙壓的作用下排出,之前閉合的孔隙連通,煤巖滲透通道增加,并且溫度升高,煤巖對(duì)氮?dú)獾慕馕芰υ龃?,但由于溫度引起的變形作用較大,但滲透率仍然是呈下降趨勢,但是變化較之前更為平穩(wěn)。
對(duì)滲透率數(shù)據(jù)使用軟件擬合,得到的恒定有效應(yīng)力下,溫度T 對(duì)滲透率k 影響的方程如下:
式中:T 為溫度,℃;a2、b2為曲線擬合系數(shù)。
各孔隙壓力下曲線的擬合系數(shù)見表3。
根據(jù)實(shí)驗(yàn)得出的滲透率數(shù)據(jù),以孔隙壓力為x軸變量,溫度為y 軸變量,滲透率為z 軸變量,使用專業(yè)的數(shù)據(jù)處理軟件,得到煤巖在溫度和孔隙壓力作用下的滲透率曲面圖以及擬合曲線圖,溫度和孔隙壓力共同作用下滲透率變化曲面圖如圖4。
分析可知,影響煤巖滲透率的因素并非等價(jià)存在,在不同階段不同因素影響所占權(quán)重不同導(dǎo)致煤巖滲透率變化出現(xiàn)差異[18]。在恒定有效應(yīng)力的條件下,滲透率在溫度或者孔隙壓力單一因素影響下均是呈現(xiàn)出逐漸降低的趨勢,并且都對(duì)二者表現(xiàn)出一定的敏感性。
圖4 溫度和孔隙壓力共同作用下滲透率變化曲面圖Fig.4 Curved surface diagram of permeability change under the combined action of temperature and pore pressure
圖2 表明,煤巖滲透率隨孔隙壓力增大呈減小趨勢,且減小的趨勢逐漸放緩。煤巖滲透率變化與吸附膨脹密切相關(guān),溫度較低時(shí),熱效應(yīng)并不明顯,該有效應(yīng)力下煤巖并未產(chǎn)生宏觀裂隙,但吸附產(chǎn)生的變形會(huì)造成煤巖內(nèi)部原生孔隙裂隙的縮小,氮?dú)獾倪\(yùn)移通道減小甚至閉合,隨著孔隙壓力進(jìn)一步增大,內(nèi)部壓實(shí)使得相同孔隙壓力差下的壓縮程度降低,滲透率降速放緩。
圖3 表明,煤巖滲透率隨溫度增大呈減小趨勢,且減小的趨勢逐漸放緩。溫度引發(fā)的熱效應(yīng)強(qiáng)化了煤巖的易損狀態(tài)。使其內(nèi)部的微觀結(jié)構(gòu)更容易發(fā)生改變。但隨著實(shí)驗(yàn)溫度的提高,受有效應(yīng)力影響,煤巖的熱膨脹效應(yīng)被作用效果被削弱,導(dǎo)致在高溫條件下,煤巖滲透率對(duì)溫度的敏感性大大降低,這也就解釋了為什么圖2 中不同溫度條件下的煤巖滲透率會(huì)出現(xiàn)“斷崖式”下降。
當(dāng)溫度超過一定閾值后,溫度引發(fā)的效應(yīng)開始占據(jù)主導(dǎo)地位;對(duì)于不同孔隙壓力下,煤巖滲透率隨溫度變化曲線都存在一定降低,并且孔隙壓力為3 MPa 時(shí),煤巖滲透率變化在40 ℃之后明顯更快,孔隙壓力分別為4、5、6 MPa 時(shí)的滲透率曲線的變化呈現(xiàn)出一定的相似性,這一點(diǎn)從擬合系數(shù)上可以看出?;谝陨戏治龊蛯?shí)驗(yàn)結(jié)果,得出3 MPa、40 ℃是實(shí)驗(yàn)孔隙壓力與溫度作用的臨界點(diǎn)。由于溫度和孔隙壓力作用效果存在臨界區(qū)間,所以擬合曲面并不是隨某一因素單調(diào)變化,而是呈現(xiàn)出在擬合曲面的兩翼向中間凹陷區(qū)域過渡的的趨勢。
使用專業(yè)數(shù)據(jù)處理軟件,對(duì)測試得出的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行非線性曲線擬合,方程選取Parabola2D模型,得出溫度和孔隙壓力對(duì)煤巖滲透率共同影響的擬合曲面參數(shù)方程:
式中:a3、b3、c3、d3、k0為曲面擬合系數(shù),分別取-0.318 5、-0.041 5、0.223、2.875×10-4、3.849 6;R2為擬合系數(shù),取0.967 8。
1)與常見的冪函數(shù)曲線不同,在恒定有效應(yīng)力下,煤巖滲透率與孔隙壓力及溫度均呈負(fù)冪指數(shù)函數(shù)關(guān)系,隨著孔隙壓力或溫度的增加,煤巖滲透率先快速降低后逐漸趨于平緩。
2)孔隙壓力和溫度通過不同的作用機(jī)制均對(duì)煤巖滲透率造成一定影響,通過分析曲線的規(guī)律,確定3 MPa、40 ℃是該實(shí)驗(yàn)煤巖孔隙壓力與溫度對(duì)滲透率影響的臨界點(diǎn)。
3)擬合得到滲透率隨孔隙壓力與溫度變化的曲線,并以此得到滲透率與孔隙壓力和溫度之間的函數(shù)關(guān)系。