川東南龍馬溪組頁(yè)巖孔隙分形特征

徐　勇，呂成福，陳國(guó)俊，仲佳愛(ài)，楊　巍，薛蓮花（.甘肅省油氣資源研究重點(diǎn)實(shí)驗(yàn)室/中國(guó)科學(xué)院油氣資源研究重點(diǎn)實(shí)驗(yàn)室，蘭州730000；2.中國(guó)科學(xué)院大學(xué)，北京00049）

川東南龍馬溪組頁(yè)巖孔隙分形特征

徐勇1，2，呂成福1，陳國(guó)俊1，仲佳愛(ài)1，2，楊巍1，2，薛蓮花1
（1.甘肅省油氣資源研究重點(diǎn)實(shí)驗(yàn)室/中國(guó)科學(xué)院油氣資源研究重點(diǎn)實(shí)驗(yàn)室，蘭州730000；2.中國(guó)科學(xué)院大學(xué)，北京100049）

基于低溫氮?dú)馕椒中螏缀螌W(xué)研究方法，對(duì)川東南龍馬溪組頁(yè)巖儲(chǔ)層進(jìn)行了孔隙分形特征研究，并運(yùn)用分形FHH模型計(jì)算了孔隙分形維數(shù)，討論了分形維數(shù)與孔隙結(jié)構(gòu)、有機(jī)碳含量、頁(yè)巖礦物成分之間的關(guān)系。研究表明：川東南龍馬溪組頁(yè)巖以中孔為主，孔隙內(nèi)部特征以墨水瓶狀孔和狹縫狀孔為主；頁(yè)巖樣品分形維數(shù)為2.629 2～2.898 0，反映了頁(yè)巖孔隙結(jié)構(gòu)復(fù)雜和非均質(zhì)性強(qiáng)的特征；頁(yè)巖平均孔徑越小、微小孔隙越多，孔隙結(jié)構(gòu)越復(fù)雜，孔表面越不規(guī)則，比表面積和分形維數(shù)則越大；有機(jī)碳含量和微孔發(fā)育程度對(duì)分形維數(shù)影響較大。通過(guò)分形維數(shù)可定量描述孔隙結(jié)構(gòu)的復(fù)雜程度和不均一性，為研究頁(yè)巖孔隙結(jié)構(gòu)的分布特征提供了思路。

頁(yè)巖儲(chǔ)層；低溫氮?dú)馕?；分形特征；龍馬溪組；川東南地區(qū)

0　引言

頁(yè)巖氣已成為非常規(guī)油氣的重要組成部分，頁(yè)巖儲(chǔ)層評(píng)價(jià)受到廣泛關(guān)注［1-3］?？紫督Y(jié)構(gòu)是儲(chǔ)層評(píng)價(jià)的核心，也是影響頁(yè)巖氣儲(chǔ)集和滲流能力等的重要因素［4-6］。因此，研究頁(yè)巖孔隙結(jié)構(gòu)對(duì)頁(yè)巖含氣性、富集規(guī)律以及勘探開(kāi)發(fā)均具有重要意義［7-8］。頁(yè)巖作為一種復(fù)雜的多孔介質(zhì)，孔隙結(jié)構(gòu)復(fù)雜，孔徑分布廣泛，其孔徑分類方法因劃分依據(jù)的不同而各有差異。Loucks等［9］將頁(yè)巖孔隙劃分為微孔（孔徑≥0.75 μm）和納米孔（孔徑＜0.75 μm）；筆者沿用國(guó)際理論和應(yīng)用化學(xué)協(xié)會(huì)（IUPAC）的孔隙分類方法［10］，將頁(yè)巖孔隙劃分為微孔（孔徑＜2nm）、中孔（2nm＜孔徑＜50nm）和大孔（孔徑＞50 nm）。

目前，研究頁(yè)巖孔隙結(jié)構(gòu)的儀器和方法較多［11］，主要為氬離子拋光-高分辨率場(chǎng)發(fā)射掃描電鏡、納米CT、雙粒子束系統(tǒng)、壓汞法及氣體吸附法等。氬離子拋光-高分辨率場(chǎng)發(fā)射掃描電鏡統(tǒng)計(jì)性差，范圍局限，不能有效反映頁(yè)巖儲(chǔ)層的孔隙結(jié)構(gòu)［12］；納米CT受三維圖像算法的影響，分析的結(jié)果具有局限性［13］；雙粒子束盡管可以表征納米級(jí)儲(chǔ)集空間，但該方法破壞了樣品的性質(zhì)［14］；壓汞法一般用于分析連通的大孔，但由于頁(yè)巖樣品在進(jìn)汞過(guò)程中的高壓使得樣品形成人工裂縫，這將對(duì)測(cè)定結(jié)果產(chǎn)生影響［15］；氣體吸附法已經(jīng)被廣泛用于表征頁(yè)巖儲(chǔ)層孔隙結(jié)構(gòu)，可獲得比表面積、孔徑分布和孔隙度等孔隙結(jié)構(gòu)參數(shù)［16］。頁(yè)巖孔隙結(jié)構(gòu)的不均一性強(qiáng)，用傳統(tǒng)的歐氏幾何原理根本無(wú)法對(duì)其進(jìn)行定量研究，而分形理論為頁(yè)巖孔隙結(jié)構(gòu)的定量研究提供了一種全新方法，它可以定量求得分形體的空間展布特征參數(shù)并定量表征孔隙結(jié)構(gòu)參數(shù)。國(guó)外學(xué)者Pfeifer等［17］最早利用分子吸附法證明儲(chǔ)層巖石孔隙具有分形特征；Krohn［18］利用掃描電鏡觀察巖石樣品，發(fā)現(xiàn)砂巖、頁(yè)巖及碳酸鹽巖中0.2～5.0 μm的孔隙具有較好的分形性質(zhì)，分形維數(shù)為2.27～2.89；國(guó)內(nèi)學(xué)者傅雪海等［19］也用分形理論研究了煤的孔隙結(jié)構(gòu)。目前對(duì)頁(yè)巖孔隙分形特征及分形維數(shù)與孔隙結(jié)構(gòu)參數(shù)的研究較少?；诖?，筆者根據(jù)低溫氮?dú)馕綄?shí)驗(yàn)結(jié)果，采用分形FHH模型，研究龍馬溪組頁(yè)巖孔隙結(jié)構(gòu)分形特征，并探討分形維數(shù)與孔隙結(jié)構(gòu)參數(shù)、有機(jī)碳含量及頁(yè)巖礦物成分的關(guān)系，這對(duì)研究復(fù)雜頁(yè)巖孔隙結(jié)構(gòu)的分布特征具有一定意義。

1　實(shí)驗(yàn)樣品與分析方法

1.1實(shí)驗(yàn)樣品

本文的實(shí)驗(yàn)樣品采自川東南地區(qū)下志留統(tǒng)龍馬溪組長(zhǎng)寧雙河—荷葉村采石場(chǎng)的新鮮泥頁(yè)巖樣品（CN-SH-01—CN-SH-12為由深到淺每隔0.5 m采樣），并采用CS-34碳硫測(cè)定儀測(cè)定了12個(gè)實(shí)驗(yàn)樣品，其中有機(jī)碳質(zhì)量分?jǐn)?shù)為0.9%～6.3%，平均為3.7%。全巖X射線衍射分析測(cè)試采用Rigaku Dmax-2500型X射線衍射儀。分析表明，樣品黏土含量較高，其體積分?jǐn)?shù)為13.3%～58.0%，平均為39.9%；石英體積分?jǐn)?shù)為15.7%～65.3%，平均為32.5%；此外，還含少量的長(zhǎng)石、碳酸鹽等礦物。樣品礦物成分分析結(jié)果如表1所列。

表1　頁(yè)巖有機(jī)碳含量與礦物成分Table 1Mineral composition and organic carbon content of shale samples

1.2分析儀器與方法

低溫氮?dú)馕綄?shí)驗(yàn)由中國(guó)科學(xué)院蘭州化學(xué)物理研究所測(cè)試完成，采用的是ASAP2020吸附儀，比表面積測(cè)試范圍下限為0.000 5 m2/g，孔徑分析范圍為0.35～500.00 nm，檢測(cè)依據(jù)為GB/T19587—2004。為了除去樣品中殘余的水分及雜質(zhì)，實(shí)驗(yàn)前樣品經(jīng)70℃的熱處理，并進(jìn)行12 h的抽真空。吸附介質(zhì)為純度99.9%的氮?dú)?，?195.8℃下測(cè)定不同相對(duì)壓力下氮?dú)獾奈搅?。以相?duì)壓力為橫坐標(biāo)，單位樣品質(zhì)量吸附量為縱坐標(biāo)，繪制低溫氮?dú)馕?脫附等溫曲線圖。根據(jù)樣品的吸附-脫附曲線，采用BJH法分析孔徑分布和孔體積，比表面積則利用BET多點(diǎn)法計(jì)算。

圖1　頁(yè)巖樣品低溫氮?dú)馕?脫附曲線Fig.1The adsorption-desorption isothermal chart of shale samples

2　頁(yè)巖氮?dú)馕?脫附曲線及孔徑分布

頁(yè)巖低溫氮?dú)馕?脫附曲線如圖1所示。頁(yè)巖的吸附（等溫）曲線與IUPAC分類的IV型吸附（等溫）曲線相近［20］。從圖1可以看出：吸附曲線與脫附曲線并不重合，產(chǎn)生滯后環(huán)。根據(jù)IUPAC分類標(biāo)準(zhǔn)［20］，將12個(gè)頁(yè)巖樣品的滯后環(huán)分為兩大類：一類為H2型，吸附-脫附曲線在相對(duì)壓力小于0.45時(shí)接近重合，而在相對(duì)壓力較高（p/p0＞0.45）時(shí)，吸附-脫附曲線分支出現(xiàn)明顯的滯后環(huán)，CN-SH-03，CNSH-04，CN-SH-06，CN-SH-07及CN-SH-12的滯后環(huán)均出現(xiàn)在相對(duì)壓力為0.45左右；另一類與H4型相近，在整個(gè)相對(duì)壓力段，吸附-脫附曲線分支基本平行，滯后環(huán)不明顯，如CN-SH-01，CN-SH-02，CN-SH-05，CN-SH-08，CN-SH-09，CN-SH-10及CN-SH-11。吸附-脫附曲線分支存在滯后環(huán)主要是由于頁(yè)巖在吸附-脫附過(guò)程中，微孔、中孔或大孔形態(tài)以及連通性存在明顯的差異（圖版Ⅰ）。H2型滯后環(huán)的吸附曲線分支由于發(fā)生毛管凝聚而逐漸上升，脫附曲線分支在相對(duì)壓力較低時(shí)突然下降，反映的孔隙以墨水瓶狀孔為主；H4型滯后環(huán)吸附-脫附曲線均上升緩慢，在相對(duì)壓力靠近1時(shí)吸附量突然增加，此類曲線對(duì)應(yīng)的孔隙一部分是均勻的狹縫狀孔，而另一部分是錐狀孔。頁(yè)巖樣品的Ⅳ型吸附等溫線和滯后環(huán)都反映出主體孔隙為中孔，也含有一定的微孔和大孔，大孔造成了在接近飽和蒸汽壓力時(shí)吸附不飽和的現(xiàn)象。

從表2可看出：分形BET模型中的比表面積為4.930～30.632 cm2/g，平均為18.650 cm2/g；孔體積為0.015 95～0.092 50 m3/g，平均為0.041 55 m3/g；根據(jù)BJH模型計(jì)算的平均孔徑為2.450～5.888 nm，平均為4.329 nm。

表2　頁(yè)巖樣品孔隙結(jié)構(gòu)參數(shù)Table 2Pore structure parameters of shale samples

根據(jù)BJH模型得出的樣品孔徑分布特征如圖2所示。從圖2可看出：頁(yè)巖孔徑分布曲線具有多峰分布的特征，CN-SH-07峰值孔徑主要為2～5 nm和40～60 nm，CN-SH-08峰值孔徑主要為35～60 nm。根據(jù)IUPAC的分類標(biāo)準(zhǔn)，確定頁(yè)巖主體孔隙為中孔，同時(shí)也含有一定的微孔和大孔。

圖2　頁(yè)巖樣品孔徑分布曲線Fig.2Pore diameter distribution of shale samples

3　孔隙分形維數(shù)

據(jù)文獻(xiàn)［21］報(bào)道，法國(guó)數(shù)學(xué)家曼德勃羅特首次提出分形幾何，它研究的是自然界中極其零碎復(fù)雜，但又具有自相似性、自反演性的體系，可以用分形維數(shù)來(lái)表示。分形維數(shù)的確定方法有多種，如氣體吸附法［22］、壓汞法［23］、小角度X射線衍射法［24］及掃描電鏡法等［25］。其中氣體吸附法中的分形Langmuir模型、分形BET模型和分形FHH模型等［26］已被廣泛運(yùn)用，特別是分形FHH模型現(xiàn)已應(yīng)用到多孔介質(zhì)的孔隙結(jié)構(gòu)特征的分形維數(shù)計(jì)算中，其計(jì)算公式如下：

式中：v為平衡壓力p下吸附的氣體分子體積，cm3/g；v0為單分子層吸附氣體的體積，cm3/g；K為待定常量；A為ln（v）和ln［ln（p0/p）］的雙對(duì)數(shù)曲線的斜率，它取決于頁(yè)巖的孔隙分形維數(shù)（D）和吸附機(jī)制；p0為氣體吸附的飽和蒸汽壓力，MPa。頁(yè)巖吸附-脫附曲線分支在相對(duì)壓力為0.45～1.00時(shí)出現(xiàn)滯后環(huán)，因此采用0.45～1.00的數(shù)據(jù)進(jìn)行孔隙分形維數(shù)計(jì)算。

不同類型多孔材料的吸附特征存在差別。國(guó)內(nèi)外學(xué)者對(duì)分形維數(shù)的計(jì)算方法存在爭(zhēng)議，并提出了以下2種計(jì)算方法［27-28］：

根據(jù)前人的研究成果［29］，確定巖石的分形維數(shù)為2～3時(shí)，符合孔隙系統(tǒng)的分形意義，可反映巖石孔表面或孔結(jié)構(gòu)的非均質(zhì)性。當(dāng)分形維數(shù)接近于2時(shí)，孔隙表面越規(guī)則；當(dāng)分形維數(shù)接近于3時(shí)，孔隙表面越不規(guī)則，表明儲(chǔ)層孔隙結(jié)構(gòu)越復(fù)雜，非均質(zhì)性越強(qiáng)；當(dāng)分形維數(shù)小于2時(shí)，已經(jīng)脫離了孔表面或孔結(jié)構(gòu)分形的意義［30］。筆者通過(guò)上述2種方法的計(jì)算（表3）發(fā)現(xiàn)式（2）計(jì)算的結(jié)果更符合孔隙系統(tǒng)的分形意義。

表3　基于分形FHH模型理論分形維數(shù)Table 3Pore fractal dimension derived from fractal FHH model

從頁(yè)巖樣品的孔隙分形計(jì)算圖（圖3）可以看出：低溫氮?dú)馕襟w積與相對(duì)壓力（p/p0＞0.45）的雙對(duì)數(shù)曲線的擬合性較好，擬合的相關(guān)系數(shù)（R2）為0.825～0.998，說(shuō)明頁(yè)巖孔隙具有很好的分形性質(zhì)。分形維數(shù)為2.629 2～2.898 0（表3），這與頁(yè)巖的孔表面不規(guī)則性或粗糙性有關(guān)［31］，代表了不同孔隙具有不同的分形特征?？紫督Y(jié)構(gòu)分形維數(shù)越大，孔隙結(jié)構(gòu)的復(fù)雜程度則越高，非均質(zhì)性就越強(qiáng)。

圖3　低溫氮?dú)馕襟w積與相對(duì)壓力（p/p0＞0.45）的雙對(duì)數(shù)曲線Fig.3Plots of lnV and ln［ln（p0/p）］（p/p0＞0.45）reconstructed from the N2gas adsorption isotherms

4　結(jié)果與討論

4.1頁(yè)巖分形維數(shù)與孔隙結(jié)構(gòu)的關(guān)系

頁(yè)巖分形維數(shù)不但反映了頁(yè)巖孔隙結(jié)構(gòu)的復(fù)雜性，而且也反映了頁(yè)巖儲(chǔ)層不同的非均質(zhì)結(jié)構(gòu)所遵循的自相似性規(guī)律。從圖4可看出：分形維數(shù)與比表面積具有正相關(guān)性（R2=0.81）［圖4（a）］，即比表面越大，分形維數(shù)也越大；分形維數(shù)與孔體積相關(guān)性較好（R2=0.70）［圖4（b）］。然而，分形維數(shù)與平均孔徑具有負(fù)相關(guān)性（R2=0.77）［圖4（c）］，即平均孔徑越小，分形維數(shù)越大。分形維數(shù)與微孔比表面積具有較好的正相關(guān)性（R2=0.76）［圖4（d）］，說(shuō)明頁(yè)巖分形維數(shù)與微孔發(fā)育程度有關(guān)。

圖4　分形維數(shù)與孔隙結(jié)構(gòu)參數(shù)的關(guān)系Fig.4Relationship between fractal dimension and pore structure parameters

4.2頁(yè)巖分形維數(shù)與有機(jī)碳含量的關(guān)系

從圖5（a）可以看出：頁(yè)巖分形維數(shù)與有機(jī)碳含量具有正相關(guān)性（R2=0.76）。隨著有機(jī)碳含量的增加和頁(yè)巖生烴能力的增強(qiáng)，生成的微孔數(shù)量逐漸增多，即有機(jī)碳含量與微孔發(fā)育程度具有正相關(guān)性（R2=0.73）［圖5（b）］，使頁(yè)巖微觀孔隙結(jié)構(gòu)復(fù)雜化，分形維數(shù)增大。Loucks等［32］研究美國(guó)Barnett頁(yè)巖發(fā)現(xiàn)，微孔主要發(fā)育在有機(jī)質(zhì)中。

圖5　TOC與分形維數(shù)、微孔比表面積關(guān)系Fig.5Relations of TOC with fractal dimension and micropore specific surface parameters

4.3頁(yè)巖分形維數(shù)與礦物成分的關(guān)系

從圖6（a）和圖6（b）可以看出：分形維數(shù)與石英、黏土含量并沒(méi)有明顯的相關(guān)性，而與方解石含量具有負(fù)相關(guān)性（R2=0.75）［圖6（c）］，與黃鐵礦含量具有正相關(guān)性（R2=0.72）［圖6（d）］。這是因?yàn)樘妓猁}巖中主要發(fā)育原生殘留孔和溶蝕孔，孔徑比較大，直徑一般為1～2 μm（圖版Ⅰ-1），而黃鐵礦集合體中草莓狀晶間孔、有機(jī)質(zhì)孔發(fā)育，孔徑比較小，直徑為10～200 nm，因此納米級(jí)孔隙的發(fā)育使孔隙結(jié)構(gòu)趨于復(fù)雜（圖版Ⅰ-2）。

圖6　分形維數(shù)與礦物成分的關(guān)系Fig.6Relationship between fractal dimension and mineral composition

6　結(jié)論

（1）川東南地區(qū)龍馬溪組頁(yè)巖主體孔隙為中孔，同時(shí)含有一定的微孔和大孔。孔隙內(nèi)部特征以墨水瓶狀孔和狹縫狀孔為主。

（2）頁(yè)巖孔隙具有很好的分形特征，根據(jù)孔隙的分形維數(shù)可以定量表征孔隙復(fù)雜程度和非均質(zhì)性。川東南地區(qū)龍馬溪組頁(yè)巖孔隙的分形維數(shù)為2.629 2～2.898 0，反映了復(fù)雜的頁(yè)巖孔隙結(jié)構(gòu)，以及孔表面的強(qiáng)非均質(zhì)性。

（3）頁(yè)巖分形維數(shù)與比表面積、孔體積具有正相關(guān)性，而與平均孔徑具有負(fù)相關(guān)性；分形維數(shù)與有機(jī)碳含量具有正相關(guān)性；方解石、黃鐵礦較石英和黏土礦物更易影響頁(yè)巖孔隙結(jié)構(gòu)的分布特征。

（References）：

［1］Curtis M E，Cardott B J，Songdergeld C H，et al.Development of organicporosityintheWoodfordshalewith increasing thermal maturity［J］.International Journal of Coal Geology，2012，103：26-31.

［2］胡海燕.富有機(jī)質(zhì)Woodford頁(yè)巖孔隙演化的熱模擬實(shí)驗(yàn)［J］.石油學(xué)報(bào)，2013，34（5）：820-825. Hu Haiyan.Porosity evolution of the organic-rich shale with thermal maturity increasing［J］.Acta Petrolei Sinica，2013，34（5）：820-825.

［3］聶海寬，邊瑞康，張培，等.川東南地區(qū)下古生界頁(yè)巖儲(chǔ)層微觀類型與特征及其對(duì)含氣量的影響［J］.地學(xué)前緣，2014，21（4）：331-343. Nie Haikuan，Bian Ruikang，Zhang Peixian，et al.Study micro-types characteristics of shale reservoir of the lower Paleozoic in southeast Sichuan basin，and their effects on the gas content［J］.Earth Science Frontiers，2014，21（4）：331-343.

［4］楊峰，寧正福，胡昌蓬，等.頁(yè)巖儲(chǔ)層微觀孔隙結(jié)構(gòu)特征［J］.石油學(xué)報(bào)，2013，34（2）：301-311. Yang Feng，Ning Zhengfu，Hu Changpeng，et al.Characterization of microscopic pore structures in shale reservoirs［J］.Acta Petrolei Sinica，2013，34（2）：301-311.

［5］徐祖新.基于CT掃描圖像的頁(yè)巖儲(chǔ)層非均質(zhì)性研究［J］.巖性油氣藏，2014，26（6）：46-49. Xu Zuxin.Heterogeneity of shale reservoirs based on CT images［J］. Lithologic Reservoirs，2014，26（6）：46-49.

［6］Mastalerz M，Schimmelmann A，Drobniak A，et al.Porosity of Devo-nian and Mississippian New Albany Shale across a maturation gradient：Insights from organic petrology，gas adsorption，and mercury intrusion［J］.AAPG Bulletin，2013，97：1621-1643.

［7］Ross D J K，Bustin R M.The importance of shale composition and pore structure upon gas storage potential of shale gas reservoirs［J］. Marine and Petroleum Geology，2009，26（6）：916-927.

［8］王玉滿，董大忠，楊樺，等.川南下志留統(tǒng)龍馬溪組頁(yè)巖儲(chǔ)集空間定量表征［J］.中國(guó)科學(xué)：地球科學(xué)，2014，44（6）：1348-1356. Wang Yuman，Dong Dazhong，Yang Hua，et al.Quantitative characterization of reservoir space in the Lower Silurian Longmaxi Shale，southern Sichuan［J］.Science China：Earth Sciences，2014，44（6）：1348-1356.

［9］Loucks R G，Reed RM，Ruppel S C，et al.Spectrum of pore types and networks in mudrocks and adescriptives classification for matrixrelated mudrocks pores［J］.AAPG Bulletin，2012，96（6）：1071-1098.

［10］Gregg S J，Sing K S W.Adsorption，surface area and porosity［M］. 2nd Edition.London：Academic Press，1982.

［11］何建華，丁文龍，付景龍，等.頁(yè)巖微觀孔隙成因類型研究［J］.巖性油氣藏，2014，26（5）：40-45. He Jianhua，Ding Wenlong，F(xiàn)u Jinglong，et al.Study on genetic type of micropore in shale reservoir［J］.Lithologic Reservoirs，2014，26（5）：40-45.

［12］田華，張水昌，柳少波，等.壓汞法和氣體吸附法研究富有機(jī)質(zhì)頁(yè)巖孔隙特征［J］.石油學(xué)報(bào)，2012，33（3）：419-427. Tian Hua，Zhang Shuichang，Liu Shaobo，et al.Determination of organic-richshaleporefeatures by mercury injection and gas adsorption methods［J］.Acta Petrolei Sinica，2012，33（3）：419-427.

［13］郝樂(lè)偉，王琪，唐俊.儲(chǔ)層巖石微觀孔隙結(jié)構(gòu)研究方法與理論綜述［J］.巖性油氣藏，2013，25（5）：123-128. Hao Lewei，Wang Qi，Tang Jun.Research progress of reservoir microscopic pore structure［J］.Lithologic Reservoirs，2013，25（5）：123-128.

［14］崔景偉，鄒才能，朱如凱，等.頁(yè)巖孔隙研究新進(jìn)展［J］.地球科學(xué)進(jìn)展，2012，27（12）：1319-1325. Cui Jingwei，Zou Caineng，Zhu Rukai，et al.New advances in shale porosity research［J］.Advances in Earth Science，2012，27（12）：1319-1325.

［15］焦堃，姚素平，吳浩，等.頁(yè)巖氣儲(chǔ)層孔隙系統(tǒng)表征方法研究進(jìn)展［J］.高校地質(zhì)學(xué)報(bào)，2014，20（1）：151-161. Jiao Kun，Yao Suping，Wu Hao，et al.Advances in characterization of poresystemofgasshales［J］.GeologicalJournalofChinaUniversities，2014，20（1）：151-161.

［16］楊正紅.氣體吸附法進(jìn)行孔徑分析進(jìn)展——密度函數(shù)理論（DFT）及蒙特卡洛法（MC）的應(yīng)用［J］.中國(guó)粉體工業(yè)，2005（6）：36-42. Yang Zhenghong.Density functional theory（DFT）and the Monte Carlo method（MC）applications-pore distribution analyze progress in gas adsorption［J］.China Powder Industry，2005（6）：36-42.

［17］Pfeifer P，Avnir D.Chemistry in nointegral dimensions between two and three［J］.ChemPhys，1983，79（7）：3369-3558.

［18］Krohn C E.Sandstone fractal and euclidean pore volume distributions［J］.GeophysRes，1988，93（B4）：3286-3296.

［19］傅雪海，秦勇，薛秀謙，等.煤儲(chǔ)層孔、裂隙系統(tǒng)分形研究［J］.中國(guó)礦業(yè)大學(xué)學(xué)報(bào)：自然科學(xué)版，2001，30（3）：225-228. Fu Xuehai，Qin Yong，Xue Xiuqian，et al.Research on fractals of pore and fracture-structure of coal reservoirs［J］.Journal of China University of Mining&Technology，2001，30（3）：225-228.

［20］徐如人，龐文琴，于吉紅，等.分子篩與多孔材料化學(xué)［M］.北京：科學(xué)出版社，2004. Xu Ruren，Pang Wenqin，Yu Jihong，et al.chemistry-zeolites andporous materials［M］.Beijing：Science Press，2004.

［21］楊宇，孫晗森，彭小東，等.煤層氣儲(chǔ)層孔隙結(jié)構(gòu)分形特征定量研究［J］.特種油氣藏，2013，20（1）：31-33. Yang Yu，Sun Hansen，Peng Xiaodong，et al.Quantitative study on fractal characteristics of the structure of CBM reservoir［J］.Special Oil&Gas Reservoir，2013，20（1）：31-33.

［22］Yin Y B.Adsorption isotherm on fractally porous materials［J］.Langmuir，1991，7（2）：216-217.

［23］劉義坤，王永平，唐慧敏，等.毛管壓力曲線和分形理論在儲(chǔ)層分類中的應(yīng)用［J］.巖性油氣藏，2014，26（3）：89-93. LiuYikun，WangYongping，TangHuimin，etal.Applicationofcapillary pressure curves and fractal theory to reservoir classification［J］. Lithologic Reservoirs，2014，26（3）：89-93.

［24］Patricia T，Nora Y K C，Chan H K，et al.Limitation of determination of surface fractal dimension using N2adsorption isotherms and modified frenkel-halsey-hill theory［J］.Langmuir，2003，19（7）：2632-2638.

［25］賀承祖，華明琪.儲(chǔ)層孔隙結(jié)構(gòu)的分形幾何描述［J］.石油與天然氣地質(zhì)，1998，19（3）：15-23. He Chengzu，Hua Mingqi.Fractal geometry description of reservoir pore structure［J］.Oil&Gas Geology，1998，19（1）：15-23.

［26］Rigby S P.Predicting surface diffusivities of molecules from equilibrium adsorption isotherms［J］.Colloids and Surfaces A：Physicochemical and Engineering Aspects，2005，262（1）：139-149.

［27］Qi H，Ma J，Wong P.Adsorption isotherms of fractal surfaces［J］. Colloids and Surfaces A：Physicochemical and Engineering Aspects，2002，206（1）：401-407.

［28］Pyun S I，Rhee C K.An investigation of fractal characteristics of mesoporous carbon electrodes with various pore structures［J］. Electrochimica Acta，2004，49（24）：4171-4180.

［29］張憲國(guó)，張濤，林承焰.基于孔隙分形特征的低滲透儲(chǔ)層孔隙結(jié)構(gòu)評(píng)價(jià)［J］.巖性油氣藏，2013，25（6）：40-45. Zhang Xianguo，Zhang Tao，Lin Chengyan.Pore structure evaluation of lowpermeabilityreservoirbasedonporefractalfeatures［J］.Lithologic Reservoirs，2013，25（6）：40-45.

［30］宋曉夏，唐躍剛，李偉，等.中梁山南礦構(gòu)造煤吸附孔分形特征［J］.煤炭學(xué)報(bào)，2013，38（1）：134-139. Song Xiaoxia，Tang Yuegang，LI Wei，et al.Fractal characteristics of adsorption pores of tectonic coal from Zhongliangshan southern coalmine［J］.Journal of China Coal Society，2013，38（1）：134-139.

［31］Yao Y，Liu D，Tang D，et al.Fractal characterization of adsorptionpores of coals from North China：An investigation on CH4adsorption capacity of coals［J］.International Journal of Coal Geology，2008，73（1）：27-42.

［32］Loucks R G，Reed R M，Ruppel S C，et al.Morphology，genesis，and distribution of nanometer-scale pores in siliceous mudstones of the Mississippian Barnett Shale［J］.Journal of Sedimentary Research，2009，79（12）：848-861.

圖版Ⅰ

（本文編輯：楊琦）

Fractal characteristics of shale pores of Longmaxi Formation in southeast Sichuan Basin

Xu Yong1，2，Lu Chengfu1，Chen Guojun1，Zhong Jiaai1，2，Yang Wei1，2，Xue Lianhua1
（1.Key Laboratory of Petroleum Resources，Gansu Province/Key Laboratory of Petroleum Resources Research，Institute of Geology and Geophysics，Chinese Academy of Sciences，Lanzhou 730000，China；2.University of Chinese Academy of Sciences，Beijing 100049，China）

：

Based on fractal theory and methods and the nitrogen adsorption results of shale samples from the Lower Silurian Longmaxi Formation in southeast Sichuan Basin，the fractal characteristics of pores were studied，the fractal dimension was calculated by using the fractal Frenkel-Halsey-Hill（FHH）model，and the relations of fractal dimension with pore structure，organic carbon content and shale composition were discussed.The results show that mesopores are the main pores in shale，and mainly consist of ink-bottle-like and slit-shaped pore.The fractal dimension of shale samples is between 2.629 2 and 2.898 0，which reflects the complexity and heterogeneity of the shale porosity.The smaller the averageporediameter is，themorethemicroporesare，themorecomplextheporestructureis，themoreirregular thespecific shale face is，the larger the specific surface face is，and the higher the fractal dimension is.Organic carbon content and thedegreeof microporedevelopment impact alot onthefractal dimension.

shale reservoir；nitrogen adsorption；fractal characteristics；Longmaxi Formation；southeast Sichuan Basin

TE121.2

A

1673-8926（2015）04-0032-08

2015-02-15；

2015-04-15

國(guó)家自然科學(xué)基金項(xiàng)目“鄂爾多斯盆地三疊系陸相頁(yè)巖微孔隙特征與天然氣賦存方式研究”（編號(hào)：41272144）和中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)（B類）“頁(yè)巖氣勘探開(kāi)發(fā)基礎(chǔ)理論與關(guān)鍵技術(shù)”（編號(hào)：XDB10010300）聯(lián)合資助

徐勇（1986-），男，中國(guó)科學(xué)院大學(xué)在讀博士研究生，研究方向?yàn)閮?chǔ)層地質(zhì)學(xué)。地址：（730000）甘肅省蘭州市城關(guān)區(qū)東崗西路382號(hào)。E-mail：xymjjx@163.com

陳國(guó)?。?967-），男，博士，研究員，主要從事儲(chǔ)層沉積學(xué)方面的研究工作。E-mail：gjchen@lzb.ac.cn。