李 軍，路 菁，李 爭(zhēng)，武清釗，南澤宇

(1.中國(guó)石化 石油勘探開發(fā)研究院，北京 100083； 2.中國(guó)石化 江漢油田分公司 勘探開發(fā)研究院，湖北 潛江 433124)

目前，頁(yè)巖氣儲(chǔ)層測(cè)井評(píng)價(jià)取得一定進(jìn)展。在巖心刻度測(cè)井條件下，能夠利用測(cè)井資料定量評(píng)價(jià)巖性組分、有機(jī)質(zhì)含量、總孔隙度、含氣性和巖石力學(xué)性質(zhì)等參數(shù)[1-4]。其中，地層巖性組分、有機(jī)質(zhì)含量、總孔隙度評(píng)價(jià)方法相對(duì)成熟，評(píng)價(jià)精度較高[5-6]。借用煤層氣評(píng)價(jià)方法中朗格繆爾方程來評(píng)價(jià)頁(yè)巖吸附氣含量，而游離氣含量則采用砂泥巖含油氣飽和度測(cè)井模型，如斯倫貝謝公司采用西門杜公式確定含油氣飽和度，進(jìn)而確定頁(yè)巖游離氣含量。由于頁(yè)巖氣儲(chǔ)層孔隙結(jié)構(gòu)和導(dǎo)電機(jī)理復(fù)雜，導(dǎo)致游離氣評(píng)價(jià)模型適應(yīng)性差，評(píng)價(jià)精度較低。基于測(cè)井對(duì)頁(yè)巖總孔隙中各微觀孔隙組分進(jìn)行定量評(píng)價(jià)，目前尚無文獻(xiàn)報(bào)道。通過系統(tǒng)巖心測(cè)試，提出頁(yè)巖氣儲(chǔ)層有機(jī)孔隙、粘土孔隙、碎屑孔隙和微裂縫等“四孔隙”評(píng)價(jià)模型和測(cè)井定量計(jì)算方法，在此基礎(chǔ)上估算游離氣含量，并對(duì)實(shí)際測(cè)井資料進(jìn)行處理和評(píng)價(jià)。

1 “四孔隙”模型建立

頁(yè)巖氣儲(chǔ)層儲(chǔ)集空間類型多，從成因來說，頁(yè)巖孔隙類型分為有機(jī)孔隙和無機(jī)孔隙。有機(jī)孔隙發(fā)育于成熟-高成熟有機(jī)質(zhì)中，是有機(jī)質(zhì)熱演化生烴膨脹作用的結(jié)果[7-10]。無機(jī)孔隙包括粘土孔隙、粉砂質(zhì)等碎屑粒間及粒內(nèi)孔隙、微裂縫孔隙。各微觀孔隙組分形狀和空間尺度不相同。依據(jù)微觀孔隙組分測(cè)井響應(yīng)差異，建立“四孔隙”模型：將總孔隙分為有機(jī)孔隙、粘土孔隙、碎屑孔隙和微裂縫，骨架分為有機(jī)骨架和無機(jī)骨架(圖1a)。圖1b展示頁(yè)巖儲(chǔ)層掃描電鏡(SEM)測(cè)試結(jié)果，顯示有機(jī)孔隙和各種無機(jī)孔隙及空間尺度分布。有機(jī)孔隙呈圓形、橢圓形及不規(guī)則管形，空間尺度以納米級(jí)別為主；粘土孔隙呈片狀、不規(guī)則狀，空間尺度為納米-微米級(jí)；碎屑孔隙(粒間、粒內(nèi))形狀不規(guī)則，以納米-微米級(jí)為主；微裂縫尺度較大，為微米-毫米級(jí)。

圖2為高精度核磁共振實(shí)驗(yàn)，也顯示多重微觀孔隙分布特征。表1展示核磁共振實(shí)驗(yàn)孔隙組分分布情況，基質(zhì)孔隙(中孔、微孔)占80%以上，大孔和微裂縫占20%以下。

基于“四孔隙”模型，利用測(cè)井資料確定總孔隙度、有機(jī)孔隙度、粘土孔隙度、碎屑孔隙度和裂縫孔隙度，為精細(xì)評(píng)價(jià)頁(yè)巖氣儲(chǔ)層奠定基礎(chǔ)。

2 測(cè)井定量計(jì)算方法

2.1 總孔隙度測(cè)井計(jì)算

由于頁(yè)巖氣儲(chǔ)層超低孔滲特征，造成實(shí)驗(yàn)測(cè)定孔隙度、滲透率等物性參數(shù)困難，目前國(guó)內(nèi)外尚無測(cè)定頁(yè)巖孔隙度、滲透率物性參數(shù)的標(biāo)準(zhǔn)方法和步驟，各商業(yè)實(shí)驗(yàn)室由于采用方法和步驟不同，即使對(duì)同一樣品進(jìn)行測(cè)量，測(cè)量結(jié)果也有很大差異[11]。

圖1 頁(yè)巖氣儲(chǔ)層“四孔隙”模型建立Fig.1 ‘Four-pore’ modeling of shale gas reservoira.“四孔隙”體積模型;b.掃描電鏡孔隙類型分析

圖2 頁(yè)巖儲(chǔ)層多重孔隙分布Fig.2 Distribution of multiple micropores of shale gas reservoir

頁(yè)巖儲(chǔ)層總孔隙包括有效孔隙(自由流體孔隙)和束縛孔隙(束縛水孔隙)。影響總孔隙度測(cè)定結(jié)果的關(guān)鍵因素包括：①束縛水和殘余油氣；②測(cè)量氣體介質(zhì)(氦氣、氮?dú)獾?的分子大?。虎勰鄮r吸附作用；④樣品預(yù)處理方法及質(zhì)量損失；⑤圍壓、孔隙壓力[6]。由于這些因素影響，造成測(cè)量結(jié)果差異。氦氣分子具有較小動(dòng)力半徑(表2)，常作為測(cè)量氣體介質(zhì)，依據(jù)氣體平衡方程確定頁(yè)巖總孔隙度。

2.1.1 基于體積模型法確定總孔隙度

測(cè)井地層評(píng)價(jià)中，經(jīng)常采用干粘土、濕粘土和泥質(zhì)三種物質(zhì)的骨架參數(shù)確定總孔隙度和有效孔隙度。干粘土只包括粘土礦物及晶間水，濕粘土則包含干粘土和束縛水(毛管束縛水和薄膜束縛水)，泥質(zhì)則包括濕粘土和少量粉砂質(zhì)。這三種物質(zhì)都有自己的骨架參數(shù)，包括骨架密度、骨架聲波時(shí)差和骨架中子值。利用干粘土骨架參數(shù)確定總孔隙度，而利用濕粘土和泥質(zhì)骨架參數(shù)確定有效孔隙度。依據(jù)體積模型，密度測(cè)井響應(yīng)方程如下：

ρb=Φtρfm+ρorgTOCv+(1-TOCv-Φt)ρm

(1)

式中:ρb,ρm,ρfm和ρorg分別為測(cè)井、干粘土骨架、孔隙中混合流體和有機(jī)質(zhì)密度，g/cm3；Φt為總孔隙度，%；TOCv為有機(jī)碳含量(體積比)，%。TOCv與TOCw(質(zhì)量比)轉(zhuǎn)換關(guān)系為：

(2)

聯(lián)立(1)和(2)式得到：

(3)

表1 頁(yè)巖巖心核磁共振孔隙度Table 1 NMR porosity analysis of shale gas cores

表2 幾種常見物質(zhì)分子直徑Tabel 2 Molecular diameters of several common materials

孔隙中混合流體密度ρfm由下式得到：

ρfm=ρhy+(ρw-ρhy)Swb

(4)

式中：ρhy和ρw分別為油氣和水的密度，g/cm3；Swb為含水飽和度,%。

2.1.2 巖心刻度測(cè)井法確定總孔隙度

依據(jù)巖心刻度測(cè)井方法，建立經(jīng)驗(yàn)公式，確定頁(yè)巖總孔隙度。在頁(yè)巖氣儲(chǔ)層孔隙度測(cè)井系列中，密度測(cè)井受井眼影響大，中子測(cè)井受粘土晶間水影響大。聲波時(shí)差測(cè)井受井眼影響較小，能夠較好地反映頁(yè)巖儲(chǔ)層總孔隙度。涪陵地區(qū)頁(yè)巖巖心氦氣測(cè)量總孔隙度與測(cè)井聲波時(shí)差呈現(xiàn)高度正相關(guān)關(guān)系，相關(guān)系數(shù)達(dá)到0.75以上。其關(guān)系式如下：

Φt=0.062DT-10.959

(5)

式中：DT為聲波時(shí)差，μs/m。

2.2 有機(jī)孔隙度測(cè)井計(jì)算

有機(jī)孔隙發(fā)育于成熟-高成熟有機(jī)質(zhì)中，有機(jī)孔隙度確定首先要確定有機(jī)質(zhì)含量。有機(jī)質(zhì)含量有兩種表示方法，一種是質(zhì)量百分比，常用于實(shí)驗(yàn)室分析；另一種是體積百分比，常用于測(cè)井評(píng)價(jià)中。

2.2.1 有機(jī)質(zhì)含量(TOC)計(jì)算

有機(jī)質(zhì)由于密度低、聲波傳播速度低、含氫指數(shù)高以及不導(dǎo)電等特性，造成鮮明的測(cè)井響應(yīng)特征，同時(shí)富有機(jī)質(zhì)頁(yè)巖還表現(xiàn)為高自然伽馬特征，在自然伽馬能譜測(cè)井響應(yīng)上表現(xiàn)為高鈾含量特征?？梢岳眠@些測(cè)井特征定量計(jì)算TOC。常用測(cè)井計(jì)算TOC方法[12-15]如下(表3)：

涪陵地區(qū)頁(yè)巖TOC與測(cè)井密度之間呈現(xiàn)強(qiáng)負(fù)相關(guān)關(guān)系，相關(guān)系數(shù)達(dá)到0.9以上。

TOC=-14.315 7DEN+39.630 2

(6)

式中：DEN為體積密度，g/cm3。

2.2.2 有機(jī)孔隙度計(jì)算

利用SEM測(cè)試技術(shù)可以直觀地確定有機(jī)孔隙大小及其分布，估算有機(jī)孔隙面孔率。利用平均面孔率對(duì)測(cè)井計(jì)算的TOCv(體積比)進(jìn)行刻度，獲得確定有機(jī)孔隙度(Φorg)關(guān)系式。

(7)

式中：a為刻度系數(shù)，與有機(jī)質(zhì)成熟度有關(guān)，由SEM分析確定。

2.3 粘土孔隙度計(jì)算

粘土孔隙是粘土礦物之間的微細(xì)孔隙，是束縛水的主要賦存空間[16]。從成因上看，粘土束縛水分為粘土礦物表面的薄膜束縛水和微小孔隙中的毛管束縛水。圖3a展示頁(yè)巖氣儲(chǔ)層中實(shí)測(cè)束縛水孔隙度與粘土含量關(guān)系，兩者呈現(xiàn)很強(qiáng)的正相關(guān)關(guān)系，相關(guān)系數(shù)達(dá)到0.89，表明粘土孔隙是頁(yè)巖氣儲(chǔ)層束縛水主要賦存空間。微細(xì)粘土孔隙表面顯現(xiàn)親水性特征，優(yōu)先吸附和儲(chǔ)集水分子，這一點(diǎn)不同于有機(jī)孔隙。圖3b展示束縛水飽和度與有機(jī)質(zhì)之間強(qiáng)烈的負(fù)相關(guān)關(guān)系，相關(guān)系數(shù)達(dá)到0.82。微細(xì)有機(jī)孔隙是有機(jī)質(zhì)成熟、演化和脫水作用的產(chǎn)物，其孔隙表面吸附烴類，其表面潤(rùn)濕性為親油性特征，是烴類的主要儲(chǔ)集空間。

圖3 頁(yè)巖氣儲(chǔ)層束縛水孔隙度與粘土礦物含量及有機(jī)質(zhì)含量關(guān)系Fig.3 Scatter diagram of bound water porosity vs.clay mineral content/TOC for shale gas reservoira.束縛水含量與粘土礦物含量關(guān)系；b.束縛水含量與有機(jī)質(zhì)含量含量關(guān)系

粘土孔隙度關(guān)系式如下：

Φclay=ΦtclayVclay

(8)

式中：Φclay為粘土孔隙度，%，Φtclay為100%粘土含量時(shí)總孔隙度，由鄰近泥巖測(cè)井響應(yīng)確定，%；Vclay為測(cè)井計(jì)算的粘土含量，%。

核磁共振實(shí)驗(yàn)通過測(cè)定孔隙中水分子中氫核橫向弛豫時(shí)間來確定巖石總孔隙度和孔隙結(jié)構(gòu)。由于氦氣分子直徑比水分子小，利用氦氣介質(zhì)能夠測(cè)量到更小的微孔隙，造成氦氣孔隙度比核磁孔隙度大。粘土含量越高，微細(xì)粘土孔隙愈發(fā)育，兩者差異也越大。圖4中2 330～2 350 m井段，粘土含量高，大于50%，微粘土孔隙發(fā)育，氦氣孔隙度與核磁孔隙度差異大，兩者差值與粘土孔隙度相當(dāng)(圖4中第6道2 330～2 350 m)。而在2 390～2 415 m井段，粘土含量低，微細(xì)粘土孔隙不發(fā)育，氦氣孔隙度與核磁孔隙度接近(圖4中第4和第5道)。

圖4 頁(yè)巖總孔隙度及粘土孔隙度計(jì)算Fig.4 Profile of the total porosity and the clay porosity of shale gas reservoir

2.4 裂縫孔隙度計(jì)算

頁(yè)巖中裂縫的發(fā)育大大提高了流體滲濾能力，是形成高產(chǎn)頁(yè)巖氣的重要地質(zhì)因素[17]。井壁成像測(cè)井是研究裂縫最重要的手段，利用成像測(cè)井可以直觀地確定裂縫類型和裂縫產(chǎn)狀，還可定量地確定裂縫密度、寬度和裂縫孔隙度等參數(shù)。此外，利用核磁共振測(cè)井也能定量確定縫洞、裂縫孔隙度[18]。

在實(shí)際生產(chǎn)中，經(jīng)常是利用雙側(cè)向測(cè)井識(shí)別和評(píng)價(jià)裂縫。在頁(yè)巖氣儲(chǔ)層中，由于電阻率較高，通常采用雙側(cè)向電阻率測(cè)井。研究表明，頁(yè)巖儲(chǔ)層中裂縫在雙側(cè)向測(cè)井上造成了顯著的響應(yīng)特征，包括絕對(duì)電阻率降低和深、淺電阻率差異性質(zhì)變化[19]。通過三維有限元法，開展雙側(cè)向測(cè)井正反演研究，確定不同裂縫產(chǎn)狀、不同基巖電阻率等條件下測(cè)井響應(yīng)特性，在此基礎(chǔ)上建立頁(yè)巖氣儲(chǔ)層裂縫孔隙度計(jì)算模型。

(9)

(10)

式中：d1,d2,d3,d4,s1,s2,s3,s4為待定系數(shù)；CLLD為深側(cè)向電導(dǎo)率,S/M;CLLS為淺側(cè)向電導(dǎo)率,S/M，X為裂縫

孔隙度與裂縫充填流體電導(dǎo)率的乘積,S/M;Cb為基巖電導(dǎo)率,S/M。

2.5 碎屑孔隙度計(jì)算

在確定總孔隙度Φt、有機(jī)孔隙度Φorg、粘土孔隙度Φclay和裂縫孔隙度Φfissure后，容易得到碎屑孔隙度Φsd。

Φsd=Φt-Φorg-Φclay-Φfissure

(11)

3 應(yīng)用實(shí)例

對(duì)來自四川盆地一口頁(yè)巖氣儲(chǔ)層測(cè)井資料采用上述方進(jìn)行定量處理和評(píng)價(jià)。地層為志留系龍馬溪組海相地層，測(cè)井資料包括自然伽馬、自然伽馬能譜、自然電位、井徑、雙側(cè)向電阻率、密度、聲波時(shí)差、中子孔隙度等。

基于四孔隙模型，利用測(cè)井?dāng)?shù)據(jù)計(jì)算總孔隙度、有機(jī)孔隙度、粘土孔隙度、碎屑孔隙度和裂縫孔隙度(圖5中第3—第7道)。游離氣主要賦存于有機(jī)孔隙、微裂縫和部分碎屑孔隙中，粘土孔隙基本為束縛水，依據(jù)這種關(guān)系可估算游離氣含量，吸附氣含量則依據(jù)巖心實(shí)驗(yàn)確定的蘭氏方程來確定，游離氣與吸附氣之和為總含氣量(圖5中的第8道)。用這種方法處理得到含氣量與現(xiàn)場(chǎng)測(cè)試結(jié)果趨勢(shì)一致。

圖5 頁(yè)巖氣儲(chǔ)層微觀孔隙組分、含氣量測(cè)井評(píng)價(jià)Fig.5 Profile of micropore components and logging evaluation of gas content of shale gas reservoir

4 結(jié)論

1) 通過關(guān)鍵井巖心巖礦、掃描電鏡、壓汞液氮聯(lián)測(cè)和核磁共振等專項(xiàng)實(shí)驗(yàn)，確定頁(yè)巖氣微觀礦物組分和孔隙組分，建立“四孔隙”評(píng)價(jià)模型。在測(cè)井資料齊備條件下，依據(jù)測(cè)井資料可以計(jì)算總孔隙度、有機(jī)孔隙度、粘土孔隙度、碎屑孔隙度和裂縫孔隙度。

2) 傳統(tǒng)砂泥巖電測(cè)井含油氣飽和度評(píng)價(jià)方法不適用于頁(yè)巖氣儲(chǔ)層。在微觀孔隙組分定量評(píng)價(jià)上，依據(jù)各個(gè)微觀孔隙組分賦存機(jī)制，估算游離氣含量。

頁(yè)巖氣儲(chǔ)層孔隙結(jié)構(gòu)復(fù)雜，導(dǎo)電機(jī)理復(fù)雜，本文旨在基于測(cè)井資料精細(xì)評(píng)價(jià)頁(yè)巖氣儲(chǔ)層，并作了初步探索。本文提出的方法雖取得一定的效果，但方法并不成熟，甚至存在錯(cuò)誤。作者期待拋磚引玉，得到同行專家們幫助和指正！

參 考 文 獻(xiàn)

[1] 郝建飛,周燦燦,李霞,等.頁(yè)巖氣地球物理測(cè)井評(píng)價(jià)綜述[J].地球物理學(xué)進(jìn)展,2012,27(4):1624-1632.

Hao Jianfei,Zhou Cancan,Li Xia,et al.Summary of shale gas evaluation applying geophysical logging [J].Progress in Geophysics,2012,27(4):1624-1632.

[2] 萬金彬,李慶華,白松濤.頁(yè)巖氣儲(chǔ)層測(cè)井評(píng)價(jià)及進(jìn)展[J].測(cè)井技術(shù),2012,36(5):441-446.

Wan Jinbin,Li Qinghua,Bai Songtao.Well-logging evaluation in shale gas reservoir and its advances [J].Well Logging Technology,2012,36(5):441-446.

[3] 莫修文,李舟波,潘保芝.頁(yè)巖氣測(cè)井地層評(píng)價(jià)的方法與進(jìn)展[J].地質(zhì)通報(bào),2011,30:400-405.

Mo Xiuwen,Li Zhoubo,Pan Baozhi.Method and advance of shale gas formation evaluation by means of well logging[J].Geological Bulletin of China,2011,30:400-405.

[4] 劉雙蓮,陸黃生.頁(yè)巖氣測(cè)井評(píng)價(jià)技術(shù)特點(diǎn)及評(píng)價(jià)方法探討[J].測(cè)井技術(shù),2011,35(2):112-116.

Liu Shuanglian,Lu Huangsheng.Evaluation methods and characteristics of log evaluation technology in shale gas[J].Well Logging Technology,2011,35(2):112-116.

[5] 謝灝辰,于炳松,曾秋楠,等.鄂爾多斯盆地延長(zhǎng)組頁(yè)巖有機(jī)碳測(cè)井解釋方法與應(yīng)用[J].石油與天然氣地質(zhì),2013,34(6)：731-736.

Xie Haochen,Yu Bingsong,Zeng Qiunan,et al.TOC logging interpretation method and its application to Yanchang Formation shales,the Ordos Basin[J].Oil & Gas Geology,2013，34(6)：731-736.

[6] Sondergeld C H,Newsham K E,Comisky J T,et al.Petrophysical considerations in evaluation and producing shale gas resources[J].2010,SPE131768:1-29.

[7] 周德華,焦方正,郭旭升,等.川東南涪陵地區(qū)下侏羅統(tǒng)頁(yè)巖油氣地質(zhì)特征[J].石油與天然氣地質(zhì)，2013，34(4)：450-454.

Zhou Dehua,Jiao Fangzheng,Guo Xusheng,et al.Geological features of the Lower Jurassic shale gas play in Fuling area,the southeastern Sichuan Basin[J].Oil & Gas Geology,2013，34(4)：450-454.

[8] 聶海寬,張金川,包書景，等.四川盆地及其周緣上奧陶統(tǒng)-下志留統(tǒng)頁(yè)巖氣聚集條件[J].石油與天然氣地質(zhì)，2012，33(3)：335-345.

Nie Haikuan,Zhang Jinchuan,Bao Shujing,et al.Shale gas accumulation conditions of the Upper Ordovician-Lower Silurian in Sichuan Basin and its periphery[J].Oil & Gas Geology,2012，33(3)：335-345.

[9] 謝忱，張金川，李玉喜，等.渝東南渝科1井下寒武統(tǒng)富有機(jī)質(zhì)頁(yè)巖發(fā)育特征與含氣量[J].石油與天然氣地質(zhì)，2013，34(1)：11-15.

Xie Chen，Zhang Jinchuan，Li Yuxi,et al.Characteristics and gas content of the Lower Cambrian dark shale in Well Yuke 1，Southeast Chongqing[J].Oil & Gas Geology,2013，34(1)：11-15.

[10] 李娟，于炳松，張金川,等.黔北地區(qū)下寒武統(tǒng)黑色頁(yè)巖儲(chǔ)層特征及其影響因素[J].石油與天然氣地質(zhì)，2012，33(3)：364-374.

Li Juan,Yu Bingsong,Zhang Jinchuan，et al.Reservoir characteristics and their influence factors of the Lower Cambrian dark shale in nor-thern Guizhou[J].Oil & Gas Geology,2012，33(3)：364-374.

[11] Rusell W Spears,David Dudus,Andrew Foulds,et al.Shale gas core analysis:strategies for normalization between laboratories and a clear need for standard materials[J].SPWLA 52nd Annual Logging Symposium May 14-18 2011.

[12] Schmoker J W.Determination of organic content of appalachian Devonian shales from formation-density logs [J].AAPG Bulletin,1979,63:1504-1509.

[13] Fertl W H,Chilingar G V.Total organic carbon content determined from well logs [J].SPE Formation Evaluation,1988,3(2).

[14] Rezaee M R,Slatt R M,Sigal R F.Shale gas rock properties prediction using artificial neural network techique and multi regression analysis,an example from a North American shale gas reservoir[C].ASEG Exterded Abstract,2007(1),1-4.

[15] Pemper R R,Han X,Mendez F E,et al.The direct measurement of carbon in wells containing oil and natural gas using a pulsed neutron mineralogy tool[J].SPE,2009,124234.

[16] 陳尚斌,朱炎銘,王紅巖,等.四川盆地南緣下志留統(tǒng)龍馬溪組頁(yè)巖氣儲(chǔ)層礦物成分特征及意義 [J].石油學(xué)報(bào),2011,32(5):775-782.

Chen Shangbin,Zhu Yanming,Wang Hongyan,et al.Characteristics and significance of mineral compositions of Lower Silurian Longmaxi Formation shale gas reservoir in the southern margin of Sichuan Basin[J].Acta Petrolei Sinica,2011,32(5):775-782.

[17] 鞏磊，曾聯(lián)波，李娟，等.南襄盆地安棚淺、中層系特低滲儲(chǔ)層裂縫特征及其與深層系裂縫對(duì)比[J].石油與天然氣地質(zhì)，2012，33(5)：778-784.

Gong Lei,Zeng Lianbo,Li Juan，et al .Features of fractures in shallow-to mid-depth reservoirs with ultra-low permeability and their comparison with those in deep reservoirs in Anpeng oilfield,the Nanxiang Basin[J].Oil & Gas Geology,2012，33(5)：778-784.

[18] 李軍,張超謨,唐小梅,等.核磁共振資料在碳酸鹽巖儲(chǔ)層評(píng)價(jià)中的應(yīng)用[J].江漢石油學(xué)院學(xué)報(bào),2004,26(1):48-50.

Li Jun,Zhang Chaomo,Tang Xiaomei,et al.Application of nuclear magnetic resonance (NMR)data in Carbonate reservoir evaluation[J].Journal of Jianghan Petroleum Institute,2004,26(1):48-50.

[19] Tan Maojin,Wang Peng,Li Jun,et al.Numerical simulation and fracture evaluation method of dual laterolog in organic shale[J].Journal of Applied Geophysics,2014 (100):1-13.