煤層儲存CO2能力的地應(yīng)力響應(yīng)特征

趙中玲

（晉能控股煤業(yè)集團陽泉有限公司，山西 陽泉 045000）

由于化石燃料（煤炭、石油、天然氣等）的高消耗，全球溫室氣體濃度在大幅上升，導(dǎo)致全球變暖等環(huán)境問題[1-4]。為了減緩環(huán)境退化的速度，不影響經(jīng)濟發(fā)展，關(guān)鍵問題是減少溫室氣體，特別是二氧化碳。而碳捕集與封存（CCS）技術(shù)是最有效的方法之一[5-6]，該系統(tǒng)可以捕獲發(fā)電和工業(yè)過程中使用化石燃料產(chǎn)生的90%的二氧化碳排放。近幾十年來，為了提高常規(guī)油氣、煤層氣和頁巖氣的注采采收率，開發(fā)了一系列地質(zhì)封存二氧化碳技術(shù)[7-8]。到目前為止，煤層CO2封存的研究成果主要包括以下5個方面：①不同儲層條件（溫度、壓力、煤變質(zhì)程度、平衡含水率）下純CO2或CH4在煤中的吸附特性[9-10]；②不同儲層條件下煤中混合CO2/CH4吸附特性[11]；③煤中飽和CH4的CO2驅(qū)油行為[12]；④CO2吸附煤基質(zhì)膨脹效應(yīng)[13]；⑤煤層CO2封存潛力評價[14-15]。這些研究將有助于深入了解煤層CO2封存的潛在效益和實際問題[16-17]，但不同應(yīng)力狀態(tài)下煤對CO2吸附的影響進(jìn)行的研究較少[18]。基于此，通過搭建煤體加載吸附實驗系統(tǒng)，對不同應(yīng)力狀態(tài)下煤體對CO2的吸附能力進(jìn)行測試研究。

1 煤樣和測試系統(tǒng)

研究的目的是在機械應(yīng)力作用下，用體積法測定二氧化碳在煤上的吸附等溫線。研究采用從5個硬煤礦采集的煤樣進(jìn)行研究：寧武煤田平朔礦區(qū)（PS），西山煤田古交礦區(qū)（GJ）、沁水煤田晉城礦區(qū)（JC）、霍西煤田霍州礦區(qū)（HZ）和河?xùn)|煤田柳林礦區(qū)（LL）。井下采集的煤樣在實驗室粉碎、篩分出粒度<0.2 mm的顆粒煤。在此，將采用顆粒煤制作成型煤，而后在不同壓力下進(jìn)行煤的CO2吸附分析。

在將粉煤倒入特制的壓模內(nèi)（外徑為40 mm，壁厚為5 mm）內(nèi)，采用液壓萬能壓力機進(jìn)行加載，實驗系統(tǒng)如圖1。試驗系統(tǒng)主要包括3部分：①加載控制系統(tǒng)：主要有伺服千斤頂、壓模罐體及加壓活塞。采用伺服千斤頂對煤樣進(jìn)行雙向加載，為了保證氣密性，在加壓活塞上設(shè)置了密封環(huán)；②變形測試系統(tǒng)：主要有引伸計和數(shù)據(jù)采集裝置，實時采集加載時的變形；③吸附充氣系統(tǒng)：主要有高壓氣瓶（CO2）、真空泵、參考罐及附屬器件，對壓制成型煤樣進(jìn)行CO2吸附性能測試。

圖1 實驗系統(tǒng)Fig.1 Testing system

2 力學(xué)性能測試

在壓制型煤時，為了確定側(cè)向推力系數(shù)η和泊松比υ，需進(jìn)行一些變形測試。為此，將引伸計對稱固定在裝置的側(cè)面（圖1），其目的是測量施加軸向壓應(yīng)力引起的應(yīng)變。應(yīng)力之間的關(guān)系如下：

式中：η為側(cè)向推力系數(shù)，η=υ/（1－υ）；υ為泊松比；σr、σz分別為徑向應(yīng)力和軸向應(yīng)力。

以沁水煤田晉城礦區(qū)（JC）煤樣為例，煤樣壓制過程徑向應(yīng)力和軸向應(yīng)力的關(guān)系如圖2。

圖2 煤樣壓制過程徑向應(yīng)力與軸向應(yīng)力的關(guān)系Fig.2 Relationship between radial stress and axial stress

根據(jù)軸向應(yīng)力和徑向應(yīng)力關(guān)系，得到的煤樣參數(shù)計算結(jié)果見表1。對于JC煤樣，加載/卸載過程應(yīng)變變化如圖3。

表1 煤樣參數(shù)計算結(jié)果Table 1 Calculation results of coal parameters

圖3 加載/卸載過程應(yīng)變變化Fig.3 Strain change during loading/unloading

從圖3可以看出，最大加載力4 kN時應(yīng)變隨加載力的變化，管壁最大應(yīng)變?yōu)?.012 1%。隨后，在卸載時煤樣中存在一定的殘余應(yīng)力，殘余應(yīng)變的大小約為0.000 3%。隨后進(jìn)行在最大加載力10 kN時進(jìn)行加載，管壁最大應(yīng)變?yōu)?.029 6%；然后卸載，殘余應(yīng)變?yōu)?.000 6%。

由此，在已知應(yīng)變值、鋼的楊氏模量和圓柱體的幾何尺寸的情況下，計算煤中的徑向應(yīng)力σr：

式中：E為鋼的彈性模量，MPa；ε為應(yīng)變；D1、D2分別為圓柱體的內(nèi)徑和外徑，m。

對于JC煤樣，加載/卸載過程徑向應(yīng)力變化如圖4。從圖4可以看出，在最大加載力4、10 kN時加載煤樣的殘余徑向應(yīng)力（即卸載力直至為0時對應(yīng)的徑向應(yīng)力）分別為2.51、4.44 MPa。

圖4 加載/卸載過程徑向應(yīng)力變化Fig.4 Change of radial stress during loading/unloading

根據(jù)已知的煤表觀密度ρb和煤基體密度ρs，型煤的孔隙率φ計算如下：

利用Micrometrics公司提供的GeoPyc儀器，采用準(zhǔn)液體法測定煤樣的表觀孔隙率。被壓制煤樣的孔隙率取決于施加的壓應(yīng)力σz，被壓制煤樣孔隙率與軸向應(yīng)力關(guān)系如圖5。

從圖5可以看出，隨著軸向應(yīng)力增大，被壓制煤樣的孔隙率呈冪指數(shù)減小，相關(guān)系數(shù)達(dá)到0.95以上。

圖5 被壓制煤樣孔隙率與軸向應(yīng)力關(guān)系Fig.5 Relationship between porosity and axial stress

PS煤樣：φ＝0.337 1×0.977 5σz（R2=0.987 4）

GL煤樣：φ＝0.322 2×0.971 1σz（R2=0.972 3）

JC煤樣：φ＝0.355 6×0.981 2σz（R2=0.991 5）

HZ煤樣：φ＝0.301 9×0.969 3σz（R2=0.976 6）

LL煤樣：φ＝0.260 3×0.968 6σz（R2=0.953 1）

3 CO2 附結(jié)果與分析

用容量法對每種煤樣在室溫內(nèi)（25℃）下進(jìn)行CO2吸附測試測試過程以下：①在壓制模中裝入松散煤樣，抽真空后進(jìn)行CO2吸附測試；②排出CO2氣體后，在4 kN壓力下壓制煤樣后，根據(jù)給定的壓力載荷σz，計算孔隙率和徑向應(yīng)力σr；③進(jìn)行CO2吸附測量，以確定已知應(yīng)力σr平均值的吸附等溫線；④分別在10、20、40 kN壓力下，重復(fù)步驟②和步驟③；⑤換另一種煤樣，重復(fù)步驟①～步驟④。

在吸附測量之前，采用真空泵對煤樣抽真空（真空度為5×10-2Pa）。而后通過頻率為0.5 Hz的壓力傳感器測量吸附過程中的壓力變化。對于每個型煤的CO2吸附，至少測試5個平衡壓力點，在高達(dá)10 MPa的范圍內(nèi)測定吸附等溫線。以霍西煤田霍州礦區(qū)HZ煤樣的CO2吸附測試結(jié)果為例，HZ煤樣在25℃時的CO2吸附等溫線如圖6。

圖6 HZ煤樣在25℃時的CO2吸附等溫線Fig.6 CO2 adsorption isotherm of HZ coal at 25℃

由圖6可知，隨著軸向應(yīng)力的增加，煤中CO2的吸附量降低。以Langmuir方程對CO2吸附等溫線進(jìn)行擬合，其相關(guān)系數(shù)達(dá)0.99以上，計算如下：

式中：q為平衡壓力p下二氧化碳的吸附量，mL/g；a為Langmuir吸附容量，mL/g；b為煤-氣系統(tǒng)的常數(shù)，MPa-1；p為氣體壓力，MPa。

對不同加載后的煤樣進(jìn)行CO2等溫吸附后，采用Langmuir擬合得到的參數(shù)a、b值及對應(yīng)的軸向應(yīng)力，不同煤樣軸向應(yīng)力對應(yīng)的吸附擬合參數(shù)見表2。

表2 不同煤樣軸向應(yīng)力對應(yīng)的吸附擬合參數(shù)Table 2 Adsorption fitting parameters for axial stress

從表2可以看出，寧武煤田平朔礦區(qū)PS煤樣，其最大吸附容量為24.41 mL/g。西山煤田古交礦區(qū)GJ煤樣對CO2的吸附能力最小。不同煤樣最大吸附量與軸向應(yīng)力的關(guān)系如圖7。

圖7 煤樣最大吸附量與軸向應(yīng)力的關(guān)系Fig.7 Relationship between maximum adsorption capacity of coal samples and axial stress

煤樣的最大吸附量隨軸向應(yīng)力的變化規(guī)律可采用線性方程a=am-kσz來描述，其中am是σz=0時的Langmuir吸附量，k為系數(shù)。朗繆爾系數(shù)在特定煤樣中的可變性見表3。由表3可知，給出了煤樣最大吸附量與軸向應(yīng)力關(guān)系試驗擬合參數(shù)am0和k，相關(guān)系數(shù)R2均達(dá)到0.92以上。其中，k為隨軸向應(yīng)力增大吸附量的遞減率。也就是說，若將軸向應(yīng)力看作地應(yīng)力，煤層的埋深對煤體吸附CO2影響明顯。埋深對煤層CO2最大吸附量影響規(guī)律如圖8。

表3 朗繆爾系數(shù)在特定煤樣中的可變性Table 3 Changeability of Langmuir’s coefficients

圖8 埋深對煤層CO2最大吸附量影響規(guī)律Fig.8 Influence of buried depth on maximum CO2 adsorption capacity in coal seam

如西山煤田古交礦區(qū)GL煤樣，最大吸附量am隨軸向應(yīng)力增大的變化量k值最小（地應(yīng)力每增大1 MPa，最大吸附量減小0.059 3 mL/g）；霍西煤田霍州礦區(qū)HZ煤樣最大吸附量am隨軸向應(yīng)力增大變化量k值最大（地應(yīng)力每增大1 MPa，最大吸附量減小0.114 4 mL/g）。隨著地應(yīng)力增大，不同煤田煤樣最大吸附量變化從0.059 3 mL/g到0.114 4 mL/g不等，平均變化約為0.089 4 mL/g。如果假設(shè)覆蓋巖石的平均密度為ρ=2 500 kg/m3，則1 000 m深度處的地應(yīng)力大約為25 MPa。在此地應(yīng)力下，煤對CO2的吸附量大約減少10.75%（從PS煤樣的7.5%到LL煤樣的13.25%）。

4 結(jié) 語

通過加載吸附試驗裝置，對5種煤樣在不同加載條件下進(jìn)行CO2吸附測試。研究結(jié)果表明：隨著軸向應(yīng)力增大，被壓制煤樣的孔隙率呈冪指數(shù)減小；隨著軸向應(yīng)力的增加，煤中CO2的吸附量呈線性降低。型煤中增大1 MPa的應(yīng)力將Langmuir的最大吸附量減少約0.43%（平均值），即在零應(yīng)力下獲得的煤的最大吸附量。地應(yīng)力對煤體吸附能力的影響非常大，因此在CO2地下儲存方面應(yīng)予以考慮。