基于CT掃描的天然氣水合物儲層微觀孔隙結(jié)構(gòu)定量表征及特征分析

董懷民, 孫建孟, 林振洲, 崔利凱, 閆偉超

(1.中國石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院,山東青島 266580; 2.青島海洋科學(xué)與技術(shù)試點國家實驗室海洋礦產(chǎn)資源評價與探測技術(shù)功能實驗室,山東青島 266071;3.中國地質(zhì)大學(xué)地球物理與空間信息學(xué)院,湖北武漢 430074; 4.中國地質(zhì)科學(xué)院地球物理地球化學(xué)勘查研究所,河北廊坊 065000)

儲集層的微觀孔隙結(jié)構(gòu)是指儲集層中巖石的孔隙與喉道的幾何形狀、尺寸分布及聯(lián)通性關(guān)系等[1-2],天然氣水合物儲集層微觀孔隙結(jié)構(gòu)復(fù)雜,其結(jié)構(gòu)特征對水合物的微觀分布、貯存含量及孔隙中流體的滲流特性等具有極其重要的影響。因此,從微觀尺度對水合物儲層的孔隙結(jié)構(gòu)進行系統(tǒng)研究,對水合物數(shù)字巖心構(gòu)建、水合物微觀分布特征及水合物數(shù)字巖心聲電特性數(shù)值模擬等研究具有重要意義。目前,巖石微觀孔隙結(jié)構(gòu)實驗室評價與表征方法呈現(xiàn)多樣化、精細化特征,包括毛管壓力曲線、掃描電鏡、鑄體薄片、聚焦離子束、核磁共振等在內(nèi)的技術(shù)與方法都得到了大量的實際應(yīng)用[3-9]。隨著CT掃描技術(shù)的廣泛應(yīng)用,基于CT成像的高精度儲層表征技術(shù)也得到了專家與學(xué)者的一致認可,大量研究結(jié)果表明,相對于其他方法,CT掃描法是建立巖石微觀孔隙結(jié)構(gòu)并對其進行定量表征的最直接、最準確的方法之一[10-16]。筆者以祁連山木里研究區(qū)水合物儲層為研究對象,運用CT掃描法構(gòu)建三維數(shù)字巖心,提取相應(yīng)的孔隙結(jié)構(gòu)網(wǎng)絡(luò)模型,實現(xiàn)巖心孔隙結(jié)構(gòu)的三維顯示與孔隙結(jié)構(gòu)參數(shù)的定量表征與獲取,此外,梳理巖心不同微觀孔隙結(jié)構(gòu)參數(shù)與水合物儲層物性的相關(guān)關(guān)系。

1 地質(zhì)背景

研究區(qū)位于青海省天峻縣木里煤田內(nèi),該地區(qū)海拔4 000～4 300 m,年平均氣溫約-5.1 ℃,長年凍土層發(fā)育極為廣泛且呈島狀分布,其平均厚度為95 m,部分地區(qū)多年凍土層厚度超過100 m[17]。研究區(qū)大地構(gòu)造處于中祁連與南祁連兩個構(gòu)造單元之間的拗陷區(qū)內(nèi),屬于北祁連深大斷裂體系在燕山期再度復(fù)活而形成的斷陷盆地[18]。研究區(qū)處于木里煤田聚乎更礦區(qū)內(nèi),部分為三疊系組成的背斜,南北兩側(cè)為侏羅系組成的兩個背斜,此背向斜南北兩側(cè)發(fā)育有規(guī)模較大的逆推覆斷裂,而在南北兩個向斜內(nèi)則發(fā)育有規(guī)模較大的北東向剪切斷裂,整個凹陷區(qū)由此被切割成南北分帶、東西分區(qū)構(gòu)造特征的不同塊段[19-21]。出露的礦區(qū)地層主要為第四系、中侏羅統(tǒng)及上三疊統(tǒng)。上三疊統(tǒng)作為基底,其巖性主要包括黑色粉砂巖、泥巖及煤等,與上覆地層侏羅系呈平行不整合接觸特征。中侏羅統(tǒng)則包括木里組及江倉組,木里組下部為辮狀河沖積平原環(huán)境中形成的中粗粒碎屑巖,夾雜碳質(zhì)泥巖或薄煤層,底部發(fā)育有底礫巖,其上部則為湖泊沼澤環(huán)境中形成的灰色粉砂巖、細粒砂巖、灰色細中粒砂巖及粗粒砂巖;江倉組下部為三角洲湖泊環(huán)境中形成的灰色細粒砂巖、中粒砂巖、深灰色泥巖及粉砂巖,其上部為淺湖—半深湖環(huán)境中形成的細碎屑泥巖、粉砂巖,伴隨灰色粉砂巖及菱鐵礦物。第四系分布較為廣泛,主要包括沖積與洪積成因腐植土、砂、礫石,坡積角礫,冰積泥沙、冰層等[21]。

通過前人對祁連山木里地區(qū)氣體組分的實測與分析,并結(jié)合該地區(qū)的凍土層厚度、地溫梯度等地層環(huán)境資料,祝有海等[17]系統(tǒng)地計算了該地區(qū)水合物形成的熱力學(xué)條件。研究結(jié)果表明,該地區(qū)具備較好的水合物形成條件。中國地質(zhì)調(diào)查局在該地區(qū)的實際鉆探試井發(fā)現(xiàn),該地區(qū)水合物主要以“孔隙型”與“裂縫型”兩種狀態(tài)產(chǎn)出,且主要賦存于中侏羅統(tǒng)江倉組,水合物儲層巖性以粉砂巖、泥巖及細砂巖等較為致密的巖石類型為主。由于多期構(gòu)造運動及斷層系統(tǒng)等影響因素的作用,導(dǎo)致儲集層巖石裂縫較為發(fā)育,即該地區(qū)水合物以“裂縫型”賦存類型為主[22]。

2 水合物儲層儲集特征

2.1 儲集空間特征

研究區(qū)巖石較為致密,孔隙連通性較差,儲集空間以溶蝕型孔隙和微裂縫為主,以原生孔隙為輔。溶蝕孔隙可分為粒內(nèi)溶孔(圖1(c))和粒間溶孔(圖1(f)),其形狀不規(guī)則,連通性較差;黏土礦物與巖石粒間孔隙(圖1(d))數(shù)量較多,但大多數(shù)半徑較小,個別半徑較大,連通性相對較差;裂縫(圖1(a))發(fā)育明顯,裂縫開度較大,在二維平面上延伸長度大,伴隨大裂縫常有微裂縫(圖1(e))發(fā)育;而原生孔隙(圖1(b))其邊緣平直,內(nèi)部清潔無雜質(zhì),孔隙半徑相對較大。

2.2 儲集物性特征

研究區(qū)水合物儲層孔隙度主要分布在1.42%～11.26%,平均為5.558%,其中小于6%的達到65%;滲透率主要分布于(0.003 8～0.295 0)×10-3μm2,平均為0.051 4×10-3μm2,其中小于0.1×10-3μm2的占比為91.03%;儲層孔滲相關(guān)性整體較差,不同水合物儲層物性差異較小。研究區(qū)儲集層物性特征表現(xiàn)為巖性致密,物性差,孔滲相關(guān)關(guān)系較差。

圖1 天然氣水合物儲層儲集空間特征Fig.1 Reservoir space characteristics of natural gas hydrate reservoir

3 水合物儲層三維孔隙模型構(gòu)建

數(shù)字巖心三維模型及其孔隙網(wǎng)絡(luò)的準確提取是開展水合物儲層巖心三維孔隙結(jié)構(gòu)定量表征及進行特征分析的基礎(chǔ)。X射線CT掃描由于其具有的直接、準確等特點,被廣泛應(yīng)用于巖石物理數(shù)值模擬研究[23-24]。在構(gòu)建數(shù)字巖心時,對巖心CT掃描直接獲取的二維灰度圖像(圖2(a))選取合理的閾值區(qū)間將灰度圖像二值化,則能夠使其轉(zhuǎn)化為二值圖像(圖2(c)),將二維圖像進行物理空間疊加即可獲取相應(yīng)的三維數(shù)值巖心。

對選取的6塊水合物儲層巖心樣品(巖心編號為A～F,對應(yīng)的巖心實物圖像依次為圖3(d)～(j))的CT掃描圖像進行空間疊加,從而構(gòu)建數(shù)字巖心三維模型,其中巖心樣品A、B的數(shù)字巖心如圖4所示,利用數(shù)值模擬方法實現(xiàn)了孔隙度與滲透率的計算,模擬結(jié)果與巖心巖石物理實驗結(jié)果對比如表1所示。由表1可知，二者具有較好的一致性,其中孔隙度平均絕對誤差為19.8%,平均相對誤差為4.76%,滲透率平均絕對誤差為0.34%,平均相對誤差為14.8%。

圖2 巖心樣品圖像二值化分割示意圖Fig.2 Schematic of image binarization segmentation for samples

圖3 天然氣水合物儲層巖心樣品圖像Fig.3 Images for core samples of natural gas hydrate reservoir

圖4 水合物儲層巖心樣品三維數(shù)字巖心模型Fig.4 Digital core mode of samples in hydrate reservior

基于構(gòu)建的數(shù)字巖心模型,應(yīng)用最大球算法[25-26],提取了與真實巖心對應(yīng)的能夠反映其孔隙空間拓撲結(jié)構(gòu)的孔隙網(wǎng)絡(luò)模型。最大球算法由Blunt與Dong等[25,27]提出,其以數(shù)字巖心孔隙像素中任意一點為球心,不斷增大球體半徑向四周延伸直至球體表面與巖石骨架接觸,包含所有像素的集合即為最大球,不同大小的球體能夠相互重疊,由此,數(shù)字巖心所有孔隙能夠被不同的最大球填滿。該方法目前存在提取的孔隙網(wǎng)絡(luò)模型的孔隙長度偏大而喉道長度偏小的缺陷。閆國亮等[28]對該方法進行了改進,彌補了上述的不足并得到了較好的應(yīng)用。運用上述方法能夠?qū)崿F(xiàn)數(shù)字巖心孔隙結(jié)構(gòu)網(wǎng)絡(luò)的提取,用以區(qū)分孔隙與喉道并確定其連通性關(guān)系。對6塊數(shù)字巖心的孔隙網(wǎng)絡(luò)模型進行提取,圖4(b)與圖4(e)為其中兩塊巖心樣品孔隙結(jié)構(gòu)的三維顯示圖像,圖4(c)與圖4(f)為提取的與巖心對應(yīng)的孔隙網(wǎng)絡(luò)球管模型,其中球狀體代表孔隙,管狀體代表喉道。

由圖4可知,水合物儲層巖心樣品與均質(zhì)性高孔滲砂巖微觀結(jié)構(gòu)特征[29]差異較大,總體而言,水合物儲層巖心孔喉分布不均勻,孔隙連通性差,并且孔隙度低,滲透率低。

表1 巖心樣品物性數(shù)值模擬結(jié)果效果對比Table 1 Comparison of property numerical simulation with rock physics experiment

4 儲層微觀孔隙結(jié)構(gòu)特征

為了實現(xiàn)水合物儲層微觀孔喉結(jié)構(gòu)特征的精細表征,從而綜合分析各個參數(shù)對水合物儲層物性特征的影響規(guī)律,需要定量地提取孔喉半徑、喉道長度、孔喉比、配位數(shù)及形狀因子等孔隙結(jié)構(gòu)參數(shù)?；跇?gòu)建的數(shù)字巖心,采用數(shù)學(xué)統(tǒng)計方法獲取了上述參數(shù),其數(shù)值模擬結(jié)果如圖5所示。

由巖心樣品孔喉參數(shù)分布數(shù)值模擬結(jié)果(圖5)可知,不同物性的水合物儲層,孔隙半徑存在較大差異,巖樣D孔隙半徑分布呈現(xiàn)雙峰狀態(tài),孔隙尺寸分布較廣,剩余巖心樣品則呈單峰狀態(tài),其孔隙半徑主要集中分布于2.5 μm以下,不同巖心樣品孔隙半徑的峰值不同,且概率分布存在差異;同樣,巖心喉道半徑也極為不同,喉道半徑除巖樣F外均呈單峰分布狀態(tài),喉道半徑主要分布于2 μm以下,其峰值約為0.5 μm,喉道尺寸分布范圍隨巖心樣品滲透率的增大而變廣,喉道半徑峰值也隨滲透率的增大而增大,從而在曲線關(guān)系上呈現(xiàn)出明顯的右偏移特征;與之相反,吼道尺寸分布范圍隨巖心樣品滲透率的減小而變窄,喉道半徑峰值隨滲透率的減小而減小;與喉道半徑分布特征相反,喉道長度分布差異相對較小,主要分布于3.7 μm以下,其峰值均約為1 μm,不同樣品的喉道長度分布差異主要在于喉道長度峰值左側(cè)的曲線形態(tài)。

圖5 孔喉參數(shù)分布曲線Fig.5 Pore throat parameters distribution curve

孔喉半徑比,又稱孔喉比,是指在局部范圍內(nèi)孔隙半徑與其相連所有喉道的半徑平均值之比,通過孔喉半徑比分布能夠有效地評價孔隙網(wǎng)絡(luò)模型的均質(zhì)性。局部孔隙空間中孔隙與喉道的差別隨孔喉比的增大而增大,即孔隙空間在微觀尺度上的變化越劇烈;局部孔隙空間的孔隙與喉道差別隨孔喉比的減小而減小,即孔隙空間在微觀尺度發(fā)育越均勻。通過統(tǒng)計巖心孔隙網(wǎng)絡(luò)模型中孔隙數(shù)量及與其相連的喉道數(shù)量,基于此計算巖心樣品的孔喉比分布(圖5)。由圖可知,巖心孔喉比差異不明顯,孔喉比相對較大,局部孔隙空間孔隙與喉道差別較大,孔隙空間在微觀尺度上變化劇烈。

巖心孔隙與喉道的形狀對水合物的微觀分布及孔隙中的流體運移具有非常重要的影響。為了對巖心孔隙與喉道進行定量表征,引入了規(guī)則的幾何形狀并將其用于孔隙與喉道復(fù)雜形狀的表征。形狀因子G定義為

G=A/P2.

(1)

式中,A為孔隙的橫截面積;P為橫截面周長。

6塊巖心樣品的孔喉形狀因子分布如圖6所示。由圖6可知,孔隙與喉道的形狀因子差異較小,表明水合物儲層巖心的孔喉結(jié)構(gòu)基本一致。

圖6 巖樣孔喉形狀因子分布Fig.6 Distribution of throat shape factors in rock samples

巖心孔隙的連通性對孔隙內(nèi)流體的運移、不同物相的微觀分布都有巨大的影響。配位數(shù)作為在微觀尺度上描述巖心孔隙連通性的重要參數(shù),其數(shù)值直接反映孔隙連通性的好壞,配位數(shù)越大,巖心孔隙連通性越好;配位數(shù)越小,巖心孔隙連通性越差。巖心的平均配位數(shù)計算如下式所示:

C=∑nαn.

(2)

式中,C為平均配位數(shù);n為配位數(shù);αn為配位數(shù)為n的孔隙占總孔隙的頻率。

基于數(shù)字巖心技術(shù)模擬巖心統(tǒng)計配位數(shù)分布,6塊巖心樣品的配位數(shù)分布如圖7所示。由圖7可知,水合物儲層巖心配位數(shù)整體較低,表明其孔隙連通性較差,是導(dǎo)致巖心滲透率低的重要因素。

綜上所述,若將1 μm作為劃分納米尺度與微米尺度的標(biāo)準,則研究區(qū)水合物儲層巖石中喉道主要分布于納米尺度,部分分布于微米尺度;與之相對,孔隙則少量處于納米尺度,主要分布于微米尺度,且數(shù)值相對較高。巖心滲透性主要由微米孔喉決定;納米尺度孔喉數(shù)量大、占比高,總體積之和較大。雖然巖心滲流能力較差,但對儲集層的儲集性能仍有重要影響。水合物儲層的物性越差,其納米孔喉發(fā)育越豐富。由此可知,水合物儲層中微米孔喉是影響滲流特性的主控因素,納米孔喉對水合物的儲量具有重要影響。

5 儲層孔隙結(jié)構(gòu)參數(shù)與儲層物性關(guān)系

巖心的孔隙半徑、喉道半徑等參數(shù)直接影響著巖心的孔隙度、滲透率等物性參數(shù)的好壞。通過巖心樣品的孔隙結(jié)構(gòu)參數(shù)與孔隙度、滲透率等儲層物性參數(shù)的關(guān)系曲線,對儲層孔喉結(jié)構(gòu)參數(shù)與儲層物性關(guān)系進行了詳細的分析。

圖8為孔喉結(jié)構(gòu)參數(shù)與儲層孔隙度及滲透率之間的相關(guān)關(guān)系。由圖8可知,水合物儲層平均孔隙半徑與孔隙度及滲透率間均呈現(xiàn)較好的相關(guān)關(guān)系,且相關(guān)性表現(xiàn)為孔隙度大于滲透率;平均喉道半徑與孔隙度及滲透率間同樣呈現(xiàn)較好的相關(guān)性,相關(guān)性也表現(xiàn)為孔隙度大于滲透率,結(jié)果表明,喉道是決定儲層滲透率的關(guān)鍵因素?？缀肀扰c滲透率間存在較好的相關(guān)關(guān)系,與之相反,孔喉比與孔隙度的相關(guān)性相對較差。巖心孔喉比越小,孔喉差異也越小,滲透率卻越高。研究區(qū)水合物儲層孔喉比較大且分布范圍較寬,孔喉比較大的巖心樣品,其孔隙度相對較小,這可能是由于眾多細小喉道圍繞孔隙分布而造成的。此外,由圖可知,基于數(shù)字巖心的孔隙結(jié)構(gòu)參數(shù)模擬結(jié)果與氣測法及數(shù)值模擬擬合結(jié)果較為相近,進一步印證了數(shù)字巖心技術(shù)的準確性,從而表明基于CT掃描對水合物儲層微觀孔隙結(jié)構(gòu)定量表征及特征分析的可靠性。

圖7 巖心樣品配位數(shù)統(tǒng)計分布Fig.7 Statistical distribution of coordination number for core samples

圖8 孔喉結(jié)構(gòu)參數(shù)與儲層孔隙度及滲透率之間的相關(guān)關(guān)系Fig.8 Relationships between pore throat parameters and porosity or permeability

6 結(jié) 論

(1)基于X射線CT掃描技術(shù),獲取了研究區(qū)水合物儲層巖心樣品的二維及三維圖像,構(gòu)建了數(shù)字巖心三維模型并提取了相應(yīng)的孔隙網(wǎng)絡(luò)模型,對孔隙結(jié)構(gòu)參數(shù)進行了數(shù)值模擬與定量提取,實現(xiàn)了水合物儲層巖心微觀孔隙結(jié)構(gòu)的定量表征。

(2)研究區(qū)水合物儲層物性較差,孔喉半徑相對較小,水合物儲集空間以粒內(nèi)溶蝕孔隙、粒間溶蝕孔隙與微裂縫為主,以原生孔隙為輔。不同巖心樣品的孔喉尺寸與形狀略有不同,孔喉主要呈現(xiàn)條帶狀與孤立狀兩種狀態(tài)分布,微孔分布具有非均質(zhì)性,局部較為發(fā)育且主要集中于粒間溶蝕孔隙。喉道尺寸小及連通性差是導(dǎo)致巖心樣品實驗滲透率低的主要原因。

(3)應(yīng)用CT掃描技術(shù)對研究區(qū)水合物儲層巖心樣品的孔隙結(jié)構(gòu)進行了定量表征與相關(guān)參數(shù)的準確提取,為開展水合物數(shù)字巖心構(gòu)建及巖石物理數(shù)值模擬奠定了堅實基礎(chǔ),同時為分析各個因素對水合物巖石物理特性的影響規(guī)律創(chuàng)造了可能,對探索阿爾奇公式計算水合物飽和度的適用性及相關(guān)參數(shù)的確定具有重要意義。