孫建孟, 劉海山, 宗成林, 秦瑞寶, 余杰
(1.中國石油大學(xué)地球科學(xué)與技術(shù)學(xué)院, 山東 青島 266580; 2.中海油研究總院, 北京 100027)
煤層氣作為一種新型的天然氣能源,在我國非常規(guī)油氣勘探開發(fā)中具有重要地位[1]。煤的工業(yè)組分變化大、物性差、流體性質(zhì)難以識別,尤其是吸附解析特性認(rèn)識不清和孔隙結(jié)構(gòu)復(fù)雜,嚴(yán)重影響了煤層氣的勘探開發(fā)[2]。煤的吸附能力在一定程度上決定了煤層的含氣量,因此,開展煤的吸附能力研究具有重要的現(xiàn)實意義。
關(guān)于煤層氣吸附能力在國內(nèi)外都有大量研究,早在1977年,KIM[3]根據(jù)等溫吸附實驗結(jié)果與巖心分析提出定量評價煤層吸附能力的KIM方程,KIM方程主要根據(jù)煤的工業(yè)組分建立與含氣量之間的關(guān)系,計算結(jié)果誤差較大,并且需要已知煤的工業(yè)組分。張超謨等根據(jù)壓汞孔喉半徑譜與核磁共振T2譜對比分析,研究結(jié)果表明,核磁共振在煤層氣儲層評價中可以有效地反映孔隙結(jié)構(gòu),能夠?qū)訁?shù)進行有效評價[4-9]。王有智等采用掃描電鏡與低溫氮吸附對煤巖孔隙結(jié)構(gòu)進行研究,結(jié)果顯示,煤巖孔隙類型主要包括屑間孔、殘余組織孔、角礫孔和氣孔,其中氣孔的廣泛發(fā)育表明了煤巖可以大量生烴,又可以為煤層氣儲存富集提供優(yōu)質(zhì)的儲集空間[10-11]。王玥等[12]利用室內(nèi)等溫吸附與現(xiàn)場解吸實驗,研究了煤層氣的儲層參數(shù)對煤層吸附能力的影響。綜上,前人一般利用等溫吸附、低溫液氮吸附實驗和煤的儲層特性評價煤的吸附能力。這些方法的影響因素很多,實驗條件苛刻,具有很大的局限性。目前,國內(nèi)外學(xué)者一致認(rèn)為煤的孔隙特性決定了煤的吸附能力。核磁共振技術(shù)可以反映煤層的微觀孔隙特征,可利用該技術(shù)對煤的吸附能力進行評價,但這方面研究很少。
本文利用核磁共振技術(shù)與低溫液氮吸附實驗的相關(guān)性,得到轉(zhuǎn)化系數(shù)的經(jīng)驗公式[13],將T2譜轉(zhuǎn)換為孔隙結(jié)構(gòu)譜,評價煤層微觀孔隙結(jié)構(gòu),進一步研究煤的吸附能力與微孔隙體積之間的定量關(guān)系,確定煤層的含氣量。
核磁共振技術(shù)具有快速測量巖石孔隙結(jié)構(gòu)且測量的孔隙直徑范圍寬等優(yōu)點,可以識別從毫米級到納米級孔隙。低溫液氮吸附實驗可以直接反映孔隙半徑分布情況,也可以直接評價煤層孔隙結(jié)構(gòu)。對同一煤樣做核磁共振與低溫液氮吸附實驗,既可以檢驗核磁共振對煤層孔隙結(jié)構(gòu)的評價結(jié)果準(zhǔn)確性,也可以將T2譜轉(zhuǎn)換為孔隙結(jié)構(gòu)譜。
實驗煤樣主要取自沁水盆地神府區(qū)塊。儀器使用紐邁MesoMR23-060H-1核磁儀進行T2譜測量,使用美國麥克公司ASAP2020比表面積及孔隙分析儀進行低溫液氮吸附測量。
核磁共振測量原理:外界靜磁場與原子核發(fā)生相互作用,使得氫核1H的核磁矩產(chǎn)生能量分級,再外加一個特定頻率垂直靜磁場方向的射頻磁場,核磁矩就會吸收能量,發(fā)生能級躍遷,產(chǎn)生核磁共振。核磁共振主要測量恢復(fù)到平衡狀態(tài)的弛豫時間,分為縱向弛豫時間T1和橫向弛豫時間T2[14]。橫向弛豫時間T2與孔隙半徑R的關(guān)系成線性
R=3ρT2
(1)
式中,R為孔隙半徑;ρ為巖石橫向表面弛豫率。
低溫液氮實驗原理:氮氣在標(biāo)準(zhǔn)大氣壓p0下冷卻至沸點溫度(-195.8 ℃)變成液體。在相同溫度下,而毛細(xì)管內(nèi)低于p0的壓力下,氮氣就凝結(jié)成液體,這種現(xiàn)象稱為毛細(xì)管冷凝。壓力p/p0與發(fā)生冷凝的孔直徑一一對應(yīng)[15],孔徑越小,產(chǎn)生冷凝所需的壓力也越小,根據(jù)這一原理,孔隙半徑R可以表示為凱爾文方程[16]
(2)
式中,R為孔隙半徑;p為測量壓力;p0為標(biāo)準(zhǔn)大氣壓。
核磁共振實驗主要參數(shù):共振頻率2 MHz回波數(shù)1 024,TE=0.2 ms,TW=6 s,信噪比大于30∶1,T2譜擬合點數(shù)64。
低溫液氮吸附實驗主要參數(shù):孔徑分析范圍0.35~500 nm,微孔區(qū)段的分辨率為0.02 nm,孔體積最小檢測0.000 1 cm3/g。
部分實驗結(jié)果見圖1、圖2。
圖1 SM-14井部分巖心核磁共振T2譜
圖2 SM-14井部分巖心液氮吸附圖
核磁共振T2譜與低溫液氮吸附實驗都是對巖石孔隙結(jié)構(gòu)分布的反映,因此,可以通過對比兩者特征差異,將T2譜轉(zhuǎn)化為孔隙結(jié)構(gòu)譜。從圖1、圖2可以明顯看出核磁共振T2譜的第2個峰與低溫液氮吸附1~100 nm之間的峰相對應(yīng),利用特殊點的相關(guān)性,可以建立兩者之間定量轉(zhuǎn)換關(guān)系(見圖3)。
圖3 利用相關(guān)點法建立T2與孔隙直徑之間的關(guān)系
該地區(qū)煤巖的橫向弛豫時間T2與孔隙直徑之間的關(guān)系大致呈線性關(guān)系
D=0.91212T2
(3)
式中,D為孔隙直徑,nm;T2為橫向弛豫時間,ms。
式(3)本質(zhì)上反映了式(1)的關(guān)系,該地區(qū)煤巖橫向弛豫率大約0.152 nm/ms。核磁共振T2譜縱坐標(biāo)為孔隙度分量,而巖心質(zhì)量體積和巖心總孔隙度可以很容易測得,所以巖樣總的孔隙體積可以計算得到,進而每個T2時間對應(yīng)的孔隙體積可以求得。再利用式(3)的T2與孔隙直徑間的關(guān)系可以將核磁共振T2譜轉(zhuǎn)換為孔隙直徑譜(見圖4)。
圖4 SM-14-17巖心液氮吸附與核磁共振對應(yīng)轉(zhuǎn)換圖
從圖4中可以看出核磁共振與低溫液氮吸附在微孔隙和中孔隙方面有較好的對應(yīng)關(guān)系,但是在反映大孔隙方面,核磁共振比低溫液氮吸附更有優(yōu)勢,這部分壓汞曲線可以與核磁共振譜有較好對應(yīng)關(guān)系。
常規(guī)利用等溫吸附體積法評價煤層吸附能力,主要是測量煤巖的蘭氏壓力與蘭氏體積,這個實驗過程耗時比較長。煤巖的吸附能力主要與煤巖的微孔隙和中孔隙發(fā)育情況有關(guān),因此,可以通過核磁共振技術(shù)測量微、中孔隙發(fā)育大小,研究與吸附能力的關(guān)系。
張錕、李相方等[17-18]研究發(fā)現(xiàn),按照十進制孔隙分類標(biāo)準(zhǔn),可以將孔徑結(jié)構(gòu)按照大小分為大孔(>1 000 nm)、中孔(1 000~100 nm)、過渡孔(100~10 nm)、微孔(<10 nm)。在此,采用孔隙直徑以10 nm為標(biāo)準(zhǔn)[通過式(3)中轉(zhuǎn)換關(guān)系得T2為16 ms],對直徑小于10 nm的微孔體積與等溫吸附實驗測量得到的吸附參數(shù)VL和pL建立關(guān)系,結(jié)果見圖5、圖6。
由圖5可以看出,隨著微孔隙體積的增加,煤樣的最大吸附能力VL逐漸增大,主要是因為微孔的孔徑較小,對甲烷分子有較強束縛能力,且微孔能提供較大的比表面積,對煤層氣吸附能力起到?jīng)Q定性作用。
圖5 VL與小于10 nm的微孔體積關(guān)系圖
圖6 pL與小于10 nm的微孔體積關(guān)系圖
從圖6可以看出隨著微孔隙體積的增加pL逐漸減小,表明微孔隙體積對低壓段的吸附作用起到?jīng)Q定作用。這主要是由于在壓力較低時,孔徑很小,尤其當(dāng)孔隙直徑與甲烷分子大小相當(dāng)時,極易將甲烷分子吸附在煤孔隙表面,使得氣體吸附量在低壓段迅速增長。因此,孔徑小于10 nm的微小孔隙是影響pL的主要因素。
通過圖5、圖6可以看出利用微中孔隙大小與吸附能力有良好的相關(guān)關(guān)系,因此,將小于10 nm的孔隙體積填充甲烷氣的空間,再利用微孔容積法[見式(3)]將其轉(zhuǎn)化為標(biāo)準(zhǔn)狀況下的體積
(4)
式中,V1為小于10 nm的微孔隙體積;ρ為甲烷氣的密度;M為甲烷的摩爾質(zhì)量;V為轉(zhuǎn)換為標(biāo)準(zhǔn)況下的微孔隙空間。
在此,將小于10 nm標(biāo)準(zhǔn)狀況下的微孔隙空間定義為核磁共振計算含氣量,將其與煤巖解析的含氣量進行對比(見圖7)。從圖7中可以看出利用核磁共振譜計算的含氣量與煤心解析出的含氣量有一定差異,大部分都處于45°線附近,只有占少數(shù)的2個點偏離較大。因此,核磁共振實驗對于評價煤層吸附能力在一定精度范圍內(nèi)比較可靠。
在實驗室核磁共振評價吸附能力的基礎(chǔ)上,將計算方法推廣到核磁共振測井,考慮到核磁共振測井存在微孔隙丟失現(xiàn)象,分別嘗試采用10、13 nm和16 nm分別對應(yīng)16、24 ms和33 ms,計算含氣量,再與煤心解析結(jié)果進行對應(yīng)比較(見圖8)。從圖8對比分析可看出,以13 nm計算效果最好,進一步對多塊巖心進行實驗計算對比(見表1)。
圖7 解析測得的含氣量與小于10 nm的孔隙體積的對比關(guān)系
從表1可以看出核磁共振計算的含氣量與煤心測得的含氣量有部分差異,有少部分巖心差異比較大,但大部分誤差都在20%以內(nèi)(平均誤差17.619%,去除差異較大的4塊后平均誤差為12.68%)。這一結(jié)果基本驗證了核磁共振測井評價煤層吸附能力的準(zhǔn)確性。
表1 多塊巖心計算結(jié)果對比結(jié)果數(shù)據(jù)表
*誤差=(計算含氣量-實測含氣量)/實測含氣量
圖8 SM-10以10 nm為標(biāo)準(zhǔn)核磁共振測井計算含氣量與巖心結(jié)果對比圖
(1) 利用核磁共振T2譜測得孔徑小于10 nm的煤層微孔隙體積與等溫吸附參數(shù)之間的良好的線性關(guān)系,應(yīng)用低溫氮吸附與核磁共振孔隙對應(yīng)關(guān)系首先找到納米孔對應(yīng)T2譜的截止值約16 ms。
(2) 評價了實驗室核磁共振譜評價煤吸附能力的應(yīng)用效果和精度,進一步遞推分析應(yīng)用到核磁測井,發(fā)現(xiàn)適當(dāng)增大該截止值到24 ms。
(3) 計算的吸附氣含量與巖心分析結(jié)果具有較好的一致性,說明依據(jù)煤層微、中孔隙體積研究煤的吸附能力是可行的,有推廣應(yīng)用前景。
參考文獻:
[1] 劉之的,楊秀春,陳彩紅,等.鄂東氣田煤層氣儲層測井綜合評價方法研究 [J].測井技術(shù),2013,37(3):289-293.
[2] 楊陸武,孫茂遠(yuǎn).中國煤層氣藏的特殊性及其開發(fā)技術(shù)要求 [J].天然氣工業(yè),2002,22(6):17-19.
[3] KIM A G.Estimating Methane Content of Bituminous Coalbeds from Adsorption Data [Z].US Dept of the Interior,Bureau of Mines,1977: 1-11.
[4] 唐曉梅,何宗斌,張超漠,等.用核磁共振T2分布定量求取孔隙結(jié)構(gòu)參數(shù)的區(qū)域性對比研究 [J].江漢石油學(xué)院學(xué)報,2003,25(4):75-77.
[5] 范英.核磁共振測井在不同油氣藏中的應(yīng)用 [J].國外測井技術(shù),2006,22(1):42-44.
[6] Mullen M J.Log Evaluation in Well Drilled for Coalbed Methane [R].Rocky Mountain Association of Geologists,1989:113-124.
[7] 簡闊,張玉貴,赫少攀,等.構(gòu)造煤甲烷吸附表面能研究 [J].煤田地質(zhì)與勘探,2014(1):31-34.
[8] 陳振宏,鄧澤,李貴中,等.煤層氣等溫吸附/解吸模擬實驗技術(shù)新進展與應(yīng)用 [J].中國石油勘探,2014,19(3):95-100.
[9] 張遵國.煤吸附/解吸變形特征及其影響因素研究 [D].重慶:重慶大學(xué),2015.
[10] 王宇航,王有智.煤儲層特征及其對吸附能力的影響 [J].中國煤炭地質(zhì),2016,28(8):33-38.
[11] 王志戰(zhàn),李新.頁巖油氣層核磁共振評價技術(shù)綜述 [J].波譜學(xué)雜志,2015,32(4):688-699.
[12] 王玥,劉會虎,桑樹勛.沁南東古城區(qū)塊3號煤層吸附解吸特征分析 [J].淮南職業(yè)技術(shù)學(xué)院學(xué)報,2016,16(5):15-20.
[13] 劉天定,趙太平,李高仁,等.利用核磁共振評價致密砂巖儲層孔徑分布的改進方法 [J].測井技術(shù),2012,36(2):119-123.
[14] 張亞蒲.核磁共振技術(shù)在煤層氣儲層評價中的應(yīng)用 [J].石油天然氣學(xué)報,2010,32(4):277-300.
[15] 藺亞兵.基于液氮吸附法對煤的孔隙特征研究與應(yīng)用 [J].煤炭科學(xué)技術(shù),2016,44(3):135-141.
[16] 侯永剛,單民瑜.關(guān)于開爾文公式的推導(dǎo)與應(yīng)用 [J].數(shù)理醫(yī)藥學(xué)雜志,2009,22(1):89-91.
[17] 張錕,侯昌海,趙迪斐,等.煤與頁巖低溫氮吸附孔隙結(jié)構(gòu)特征與分形特征對比——以陽泉地區(qū)山西組15號煤與頁巖為例 [J].科學(xué)技術(shù)與工程,2016(29):68-75.
[18] 李相方,蒲云超,孫長宇,等.煤層氣與頁巖氣吸附/解吸的理論再認(rèn)識 [J].石油學(xué)報,2014,35(6):1113-1129.