李華陽(yáng)，李潮流，周燦燦，胡法龍，李霞，鄒友龍

（中國(guó)石油勘探開(kāi)發(fā)研究院，北京 100083）

0 引 言

據(jù)翟光明、邱中建院士等［1－2］的研究，中國(guó)陸上未來(lái)油氣勘探的主要領(lǐng)域包括前陸盆地、斷陷盆地富油氣凹陷、坳陷盆地大面積地層巖性帶和大型古隆起。從測(cè)井專(zhuān)業(yè)角度看，目前及未來(lái)油氣勘探主體與以往常見(jiàn)油氣藏的區(qū)別主要是儲(chǔ)層的復(fù)雜化，即巖性、孔隙空間類(lèi)型和流體分布特征多樣化，給測(cè)井采集、解釋評(píng)價(jià)與巖石物理研究帶來(lái)很大挑戰(zhàn)。上述4類(lèi)勘探對(duì)象中，復(fù)雜儲(chǔ)層主要包括復(fù)雜碎屑巖（特低孔隙度滲透率、致密砂巖、頁(yè)巖等）、縫洞型碳酸鹽巖、火山巖儲(chǔ)層。其中，以鄂爾多斯盆地、四川盆地為代表的致密砂巖油氣區(qū)勘探開(kāi)發(fā)近2年取得重要突破，已成為目前最為現(xiàn)實(shí)的非常規(guī)油氣資源［3］。該類(lèi)儲(chǔ)層也是本文論述的對(duì)象。

致密砂巖儲(chǔ)層主要特征是孔喉細(xì)小、孔隙度滲透率條件極差、孔隙結(jié)構(gòu)極復(fù)雜。針對(duì)常規(guī)砂巖儲(chǔ)層條件研發(fā)的巖石物理分析技術(shù)和手段已難以滿足致密砂巖儲(chǔ)層實(shí)驗(yàn)精度和時(shí)效性的要求，特別是常規(guī)驅(qū)替實(shí)驗(yàn)不適用于此類(lèi)儲(chǔ)層［3］。以CT掃描和數(shù)值模擬為核心的測(cè)井?dāng)?shù)字巖石物理性質(zhì)研究已逐漸成為解決相關(guān)難題的重要手段，并在滲透率、電性模擬等方面取得重要進(jìn)展，但同時(shí)也遇到了很多挑戰(zhàn)。

1 測(cè)井?dāng)?shù)字巖石物理研究歷程與主要內(nèi)容

利用數(shù)值模擬的方法分析巖石滲流等宏觀物理性質(zhì)的研究早在20世紀(jì)50年代就由Purcell等［4］提出，早期有毛細(xì)管束模型、二維電路網(wǎng)絡(luò)模型等。20世紀(jì)70年代由于計(jì)算機(jī)及圖像分析技術(shù)的出現(xiàn)使得孔隙網(wǎng)絡(luò)模型的研究得到進(jìn)一步發(fā)展，但直到20世紀(jì)90年代以后隨著CT技術(shù)的發(fā)展才真正開(kāi)始了模擬砂巖三維孔隙的數(shù)字巖石物理研究。

美國(guó)科學(xué)院院士Amos Nur教授1977年在斯坦福大學(xué)組建SRB（Stanford Rock Physics ＆ Borehole Geophysics Project）并致力于地層條件下的巖石物理性質(zhì)研究，1995年左右開(kāi)始數(shù)字巖石物理性質(zhì)研究，目前在石油天然氣行業(yè)處于領(lǐng)先地位，引領(lǐng)相關(guān)學(xué)科的發(fā)展方向，并于2007年創(chuàng)建In Grain公司進(jìn)行相關(guān)技術(shù)服務(wù)和咨詢。Schlumberger、Shell和Chevron等服務(wù)公司和油公司都較早地開(kāi)展了數(shù)字巖石物理性質(zhì)方面的研究工作。

中國(guó)石油大學(xué)孫建孟教授、陶果教授和姚軍教授等在2000年以后開(kāi)始利用格子氣自動(dòng)機(jī)、LBM及有限元等方法模擬電學(xué)、聲學(xué)及滲透率等性質(zhì)以及兩相流模擬［5］。畢業(yè)于英國(guó)帝國(guó)理工大學(xué)的董虎博士于2010年在廈門(mén)注冊(cè)成立了數(shù)巖科技公司，專(zhuān)門(mén)從事相關(guān)技術(shù)研發(fā)和服務(wù)。

綜合國(guó)內(nèi)外關(guān)于數(shù)字巖石物理研究的發(fā)展歷程，測(cè)井?dāng)?shù)字巖石物理是依據(jù)物理學(xué)原理和方法，基于反映巖石內(nèi)部孔隙特征的圖像采用數(shù)字模擬手段分析儲(chǔ)層巖石的物理性質(zhì)和輸運(yùn)特性新興學(xué)科，主要內(nèi)容包括2部分。

（1）數(shù)字孔隙格架的構(gòu)建。目前用于構(gòu)建數(shù)字孔隙格架的方法主要有物理實(shí)驗(yàn)法和數(shù)值重建法2種。物理實(shí)驗(yàn)法利用CT或掃描電鏡對(duì)巖心內(nèi)部結(jié)構(gòu)進(jìn)行三維灰度成像，再通過(guò)圖像二值化處理分離孔隙和骨架，以獲取巖石真實(shí)三維孔隙格架。隨著CT掃描技術(shù)和成像技術(shù)的進(jìn)步，該方法的應(yīng)用已經(jīng)越來(lái)越普遍。國(guó)外還有學(xué)者提出采用環(huán)氧樹(shù)脂灌注孔隙并通過(guò)聚焦掃描進(jìn)行孔隙成像的實(shí)驗(yàn)方法，但只能處理巖心薄片［6］。

數(shù)值重建法基于具體的巖心薄片圖像分析數(shù)據(jù)，利用統(tǒng)計(jì)模型或模擬沉積與成巖過(guò)程的方法建立數(shù)字孔隙格架，典型算法有高斯場(chǎng)法、模擬退火法、沉積過(guò)程重建法和馬爾可夫隨機(jī)重建法等［7－11］。需要指出的是，沉積過(guò)程重建方法采用不同尺寸球形、橢球形顆粒堆積并模擬沉積巖的沉積、壓實(shí)、膠結(jié)等成巖過(guò)程，最終得到近似于真實(shí)砂巖的三維孔隙格架，與其他幾種數(shù)值重建方法相比較，其結(jié)果更接近真實(shí)巖心特征［8］。

無(wú)論采用哪種方法，在獲取三維數(shù)字孔隙格架的同時(shí)還必須有效區(qū)分孔隙和喉道，這是對(duì)孔隙網(wǎng)絡(luò)的一種近似簡(jiǎn)化以保證模擬計(jì)算成功的關(guān)鍵。

（2）巖石物理特性的模擬。在獲取三維孔隙格架的基礎(chǔ)上利用多種算法模擬巖石特征參數(shù)（見(jiàn)圖1）。圖1表明，目前的測(cè)井?dāng)?shù)字巖石物理性質(zhì)研究已經(jīng)可以實(shí)現(xiàn)包括電、聲、核磁共振及滲流等主要參數(shù)的模擬，基本涵蓋了常規(guī)實(shí)驗(yàn)測(cè)試項(xiàng)目，可對(duì)致密砂巖的巖石物理特征進(jìn)行較全面的分析研究。

圖1 基于三維孔隙格架的巖石物理性質(zhì)數(shù)值模擬對(duì)象

2 近2年研究進(jìn)展概述

（1）掃描測(cè)試儀器配套能夠獲取各種分辨率圖像。以美國(guó)Xradia公司為代表針對(duì)巖心的CT掃描設(shè)備逐漸配套完善，目前已商業(yè)化應(yīng)用的包括UltraXRM－L200型、Micro XCT－400型、MicroXCT－200型和VersaXRM－500型。美國(guó)GE公司也有類(lèi)似CT設(shè)備，它們的最高成像分辨率分別為50 nm、0.5μm和1μm，所使用的X射線能量一般為120～160 keV，納米CT的X射線能量更低，為8 ke V。此外，美國(guó)FEI公司最新推出了分辨率高達(dá)0.8～3.5 nm級(jí)聚焦離子束掃描電鏡FIBSEM（Focused Ion Beam Scanning Electron Microscopes），通過(guò)對(duì)巖心表面原子級(jí)別的逐層切割輸出幾個(gè)納米級(jí)的高清圖像。

（2）發(fā)展了模擬沉積與成巖過(guò)程的數(shù)字巖心構(gòu)建方法。迄今已發(fā)展了多種較為成熟的數(shù)值重建方法。圖2給出了沉積過(guò)程重建方法的示意圖。該方法綜合考慮了形成沉積巖的顆粒尺寸分布、顆粒形狀、壓實(shí)與膠結(jié)及次生加大等成巖過(guò)程，能夠得到與致密砂巖類(lèi)似的孔隙格架。

（3）孔隙格架的數(shù)學(xué)表征方法不斷完善。巖心CT掃描的圖像是灰度圖像，灰度級(jí)別反映了X射線衰減程度，它主要受礦物密度和X射線傳播距離的影響。致密砂巖一般具有巖屑礦物種類(lèi)雜、含量多的特點(diǎn)［3］，不同類(lèi)型礦物的密度不同，因而其CT圖像一般具有較寬的灰度分布范圍。通過(guò)圖像處理技術(shù)可較好地將孔隙和骨架進(jìn)行分割，得到只包含孔隙（像素值為0）、骨架背景（像素值為1）的二值圖像。

基于分割后的二值圖像，目前已形成多種比較有效識(shí)別孔隙和喉道的算法，如多向掃描法、中軸線法、多面體法和最大球法。應(yīng)用這些處理方法可以建立具有真實(shí)拓?fù)浣Y(jié)構(gòu)的孔隙格架。圖3為最大球法區(qū)分孔隙和喉道的示意圖。在二值圖像中任意孔隙位置放置一個(gè)球體，不斷增加球體半徑直至遇到孔隙壁為止，在每一個(gè)孔隙像素位置均獲得這樣的球體，用局部最大球表示孔隙體［見(jiàn)圖3（a）］，連接相鄰孔隙體的所有球體表示喉道［見(jiàn)圖3（b）］，從而得到三維孔隙格架信息數(shù)據(jù)［見(jiàn)圖3（c）］。

引入歐拉數(shù)表征三維孔隙網(wǎng)絡(luò)的連通性，

式中，N為網(wǎng)絡(luò)中直徑大于某一尺寸的所有孔隙個(gè)數(shù)；C為大于該尺寸的所有連通喉道個(gè)數(shù)；H為溶洞個(gè)數(shù)，砂巖中一般取0；V為樣品體積。

圖2 過(guò)程模擬法構(gòu)建數(shù)字孔隙格架的方法原理示意圖

圖3 最大球法區(qū)分孔隙、喉道的原理示意圖

分析式（1）可知，當(dāng)網(wǎng)絡(luò)中連通喉道個(gè)數(shù)C大于孔隙個(gè)數(shù)N時(shí)，λ為負(fù)值，表示連通程度高；反之，當(dāng)N＞C時(shí)，網(wǎng)絡(luò)中存在孤立孔隙，λ為正值，連通程度低。通過(guò)不斷改變觀測(cè)尺度得到一組λ值，繪制出如圖4所示的歐拉特征圖譜。圖4中歐拉曲線與橫坐標(biāo)的交點(diǎn)對(duì)應(yīng)觀測(cè)尺寸定義為歐拉特征半徑d0，單位為μm，它是表征復(fù)雜砂巖孔喉網(wǎng)絡(luò)連通程度的臨界值，尺寸小于d0的孔喉連通性好，大于d0的孔喉連通性差。歐拉特征半徑是反映復(fù)雜孔隙結(jié)構(gòu)碎屑巖儲(chǔ)層滲流能力和儲(chǔ)層品質(zhì)的關(guān)鍵參數(shù)。

圖4 利用歐拉數(shù)表示孔隙連通性的原理示意圖

除歐拉數(shù)以外，描述孔隙結(jié)構(gòu)的參數(shù)還包括連通孔隙體積比、孔隙尺寸分布函數(shù)、局部孔隙度分布函數(shù)和局部滲流概率函數(shù)等［7］。

（4）形成多種基于數(shù)字孔隙格架的巖石物理屬性模擬方法。基于三維數(shù)字孔隙格架模擬儲(chǔ)層滲透性、電性主要采用基爾霍夫電路節(jié)點(diǎn)法、格子玻爾茲曼方法、隨機(jī)步行者法等，核磁共振性質(zhì)的模擬采用隨機(jī)步行者法，巖石彈性力學(xué)的模擬主要采用有限元法逐像素點(diǎn)分析骨架在應(yīng)力作用下的應(yīng)變分布。基爾霍夫電路節(jié)點(diǎn)法的基本原理如圖5所示。

對(duì)于復(fù)雜電路網(wǎng)絡(luò)，基爾霍夫電流定律（KCL）指出，對(duì)于任意一個(gè)閉合回路或閉合區(qū)域，其中的各段電壓和電流滿足一定的約束關(guān)系［見(jiàn)圖5（a）］，是電路分析的基本理論依據(jù)。根據(jù)該原理，在三維孔隙格架中［見(jiàn)圖5（b）］，假設(shè)兩端施加一定電壓ΔV，應(yīng)用歐姆定律對(duì)于任一孔隙節(jié)點(diǎn)可以計(jì)算出其電流強(qiáng)度I［見(jiàn)圖5（c）］，相鄰孔隙之間滿足KCL定律。對(duì)所有孔隙節(jié)點(diǎn)都得到一個(gè)方程，求解超大方程組就可以計(jì)算電阻率。

圖5 基爾霍夫電路節(jié)點(diǎn)法模擬電阻率的原理示意圖

上述模擬過(guò)程中如果假設(shè)所有孔隙均飽含水（電阻率為Rw），則模擬計(jì)算的是巖石飽含水的電阻率R0。如果進(jìn)一步假設(shè)孔隙由大到小依次飽含油，則依次計(jì)算出不同含水飽和度對(duì)應(yīng)的電阻率Rt，從而模擬得到巖石的I－Sw關(guān)系（見(jiàn)圖6）。圖6中，圓圈代表某一低孔隙度低滲透率砂巖樣品的I－Sw模擬結(jié)果，紅、藍(lán)、綠色3條直線分別代表飽和度指數(shù)n分別取1、2、3的理論結(jié)果，表明該樣品的ISw關(guān)系曲線呈非線性變化。但為了簡(jiǎn)化，其n值也可以近似取2.0。

格子玻爾茲曼方法是采用人工微觀模型，將流體抽象為大量的微觀粒子并按照某種簡(jiǎn)單規(guī)則在離散的格子上碰撞和遷移。通過(guò)對(duì)粒子的運(yùn)動(dòng)進(jìn)行統(tǒng)計(jì)分析可以模擬流體運(yùn)動(dòng)的宏觀特性?；谶@一原理，該方法常被用于模擬巖石的電性和滲透性，圖7為一種格子玻爾茲曼方法模擬計(jì)算常用的D3Q19計(jì)算模型［9］，它假設(shè)粒子在每個(gè)格點(diǎn)上可能存在19個(gè)遷移方向，實(shí)際模擬運(yùn)算時(shí)須對(duì)每一個(gè)方向進(jìn)行分析計(jì)算。

根據(jù)以上介紹可以看出，基爾霍夫方法是通過(guò)簡(jiǎn)化為孔隙、喉道的三維格架建立方程組計(jì)算電阻率，每一節(jié)點(diǎn)上電流的方向是確定的，具體取決于與孔隙連接的喉道個(gè)數(shù)，即使在復(fù)雜孔隙結(jié)構(gòu)的致密砂巖中孔喉比一般都比較小。格子玻爾茲曼方法在每一節(jié)點(diǎn)都要考慮多達(dá)19個(gè)方向，因此，基爾霍夫方法的模擬計(jì)算速度要明顯占優(yōu)勢(shì)，但二者的計(jì)算精度基本相當(dāng)。

圖7 格子玻爾茲曼模擬常用的D3Q19計(jì)算模型示意圖

隨機(jī)步行者方法常用來(lái)模擬巖石的核磁共振特性，基本原理是在孔隙格架中隨機(jī)放置若干個(gè)質(zhì)子，這些質(zhì)子以布朗運(yùn)動(dòng)與孔隙壁碰撞并湮沒(méi)或反彈。通過(guò)設(shè)置一定的時(shí)間步長(zhǎng)計(jì)算質(zhì)子逐漸消亡導(dǎo)致的宏觀磁矩不斷減小的過(guò)程就可以模擬核磁共振信號(hào)。該方法原理簡(jiǎn)單，只要CT分辨率足夠，其模擬的結(jié)果與實(shí)測(cè)巖心的核磁共振信號(hào)基本相當(dāng)。

3 測(cè)井?dāng)?shù)字巖石物理研究面臨的挑戰(zhàn)

（1）分辨率與儲(chǔ)層非均質(zhì)性之間存在矛盾?chē)?yán)重影響模擬精度。致密砂巖儲(chǔ)層主要發(fā)育納米級(jí)孔喉［12］，為提高分辨率以觀測(cè)更細(xì)微孔隙就必須減小樣品尺寸，納米級(jí)掃描要求樣品尺寸小于0.1 cm，而對(duì)具有強(qiáng)烈微觀非均質(zhì)性的砂巖而言，小尺寸樣品的模擬結(jié)果與儲(chǔ)層宏觀巖石物理性質(zhì)相差較大。另一方面，大尺寸樣品的CT掃描只能揭示控制滲流的主要喉道，卻無(wú)法刻畫(huà)對(duì)電性影響顯著的微細(xì)孔喉，該矛盾導(dǎo)致的結(jié)果可以用圖8解釋。圖8是4塊直徑2.54 cm的砂巖樣品采用基爾霍夫電路節(jié)點(diǎn)法模擬的滲透率、地層因素與實(shí)驗(yàn)測(cè)量結(jié)果對(duì)比。可以看出，在14.27μm的CT分辨率下（對(duì)應(yīng)樣品直徑2.54 cm）依據(jù)所能夠觀察到的大孔喉模擬的滲透率與實(shí)驗(yàn)值在數(shù)量級(jí)上基本接近［見(jiàn)圖8（a）］，說(shuō)明大孔喉是控制儲(chǔ)層滲流能力的關(guān)鍵；但地層因素F的模擬結(jié)果與實(shí)驗(yàn)值相差達(dá)2～3個(gè)數(shù)量級(jí)［見(jiàn)圖8（b）］，這主要是由于其中大量納米級(jí)微細(xì)孔喉未被CT探測(cè)到、模擬計(jì)算時(shí)其導(dǎo)電作用被忽略所致。

（2）圖像分割與構(gòu)建數(shù)字孔隙格架過(guò)程的多解性強(qiáng)。圖9（a）、（b）分別是1塊常規(guī)中高孔隙度滲透率砂巖樣品和1塊致密砂巖樣品2 000張CT圖片灰度分布對(duì)比。表明致密砂巖的圖像灰度分布表現(xiàn)為強(qiáng)烈的非均質(zhì)性，圖像處理過(guò)程中很難選取準(zhǔn)確的灰度閾值以實(shí)現(xiàn)二值化分割。圖9（c）顯示典型的致密砂巖中發(fā)育極不規(guī)則的孔隙和喉道，很難用三角形、球形等簡(jiǎn)單的幾何特征描述，因此在致密砂巖中，基于圖像準(zhǔn)確構(gòu)建三維數(shù)字孔隙格架的過(guò)程存在很強(qiáng)的多解性，其結(jié)果也將影響模擬精度。

圖8 4塊樣品基爾霍夫電路節(jié)點(diǎn)法模擬滲透率、地層因素的精度分析

圖9 典型致密砂巖CT圖片灰度分布與復(fù)雜幾何形狀的孔隙特征

（3）缺乏能夠模擬黏土礦物影響的方法模型。黏土礦物一般以分散狀、結(jié)構(gòu)狀或?qū)訝?種形式分布在砂巖儲(chǔ)層中，不同的黏土礦物類(lèi)型和分布形式對(duì)儲(chǔ)層滲透率、電阻率的影響非常顯著［13－14］，特別是類(lèi)似于渤海灣盆地歧口凹陷中深層沙河街組、巖性較細(xì)、泥質(zhì)含量和壓實(shí)程度高、成巖作用強(qiáng)的致密砂巖儲(chǔ)層［15］，這種影響就更加突出。但綜合目前國(guó)內(nèi)外的研究情況，考慮儲(chǔ)層潤(rùn)濕特征影響多相流的模型較多，但尚欠缺考慮黏土礦物影響巖石物理屬性的數(shù)字孔隙格架構(gòu)建及相關(guān)模擬方法的研究。

（4）現(xiàn)有模擬方法的時(shí)效性亟需優(yōu)化完善。受數(shù)據(jù)量大、計(jì)算過(guò)程復(fù)雜等因素制約，現(xiàn)有模擬方法耗時(shí)長(zhǎng)，采用當(dāng)前市場(chǎng)上最好的高性能微機(jī)工作站模擬運(yùn)算1次一般需要近20 h，時(shí)效性差是制約現(xiàn)有方法推廣應(yīng)用的關(guān)鍵，亟需進(jìn)行算法的優(yōu)化完善，開(kāi)發(fā)并行算法程序。

4 研究方向展望

（1）趨勢(shì)云分析技術(shù)。針對(duì)分辨率與儲(chǔ)層非均質(zhì)性之間的矛盾，對(duì)一塊常規(guī)尺寸致密砂巖樣品先進(jìn)行粗掃，根據(jù)粗掃圖像選擇孔隙發(fā)育程度或連通程度不同的多個(gè)部位切割出小尺寸樣品再進(jìn)行精細(xì)掃描，分別基于粗掃圖像和細(xì)掃圖像模擬儲(chǔ)層的巖石物理參數(shù)，最后結(jié)合不同分辨率的模擬結(jié)果進(jìn)行趨勢(shì)分析，判斷滲透率、電阻率等屬性的變化規(guī)律和主要控制因素。

（2）多分辨率三維孔隙網(wǎng)絡(luò)整合。與趨勢(shì)云分析技術(shù)相對(duì)的是多分辨率孔隙網(wǎng)絡(luò)整合，基本思路是對(duì)大尺寸巖心按照相同的思路進(jìn)行逐級(jí)精細(xì)掃描，對(duì)每個(gè)級(jí)別的掃描圖像分別統(tǒng)計(jì)其歐拉數(shù)及孔隙尺寸分布函數(shù)，最后在不改變拓?fù)浣Y(jié)構(gòu)特征的前提下將高分辨率圖像整合到粗掃圖像中，最終獲得能夠反映大尺寸巖心的高分辨率圖像，據(jù)此進(jìn)行模擬分析。

5 結(jié)論與認(rèn)識(shí)

（1）為滿足致密砂巖等非常規(guī)儲(chǔ)層巖石物理研究需要，急需深化和完善測(cè)井?dāng)?shù)字巖石物理模擬與分析技術(shù)。

（2）以致密砂巖為代表的復(fù)雜儲(chǔ)層巖石物理研究不能沿用常規(guī)思路，將高精度實(shí)驗(yàn)?zāi)M與數(shù)值模擬相結(jié)合是一種必然趨勢(shì)。

（3）當(dāng)前的測(cè)井?dāng)?shù)字巖石物理模擬技術(shù)尚不完備，不能完全代替實(shí)驗(yàn)?zāi)M，亟需進(jìn)行理論模型和計(jì)算方法的優(yōu)化完善。

［1］ 翟光明，何文淵.從區(qū)域構(gòu)造背景看我國(guó)油氣勘探方向 ［J］.中國(guó)石油勘探，2005，2（10）：1－8.

［2］ 邱中建，康竹林，何文淵.從近期發(fā)現(xiàn)的油氣新領(lǐng)域展望中國(guó)油氣勘探發(fā)展前景 ［J］.石油學(xué)報(bào)，2002，4（23）：1－6.

［3］ 趙政璋，杜金虎，等.致密油氣 ［M］.北京：石油工業(yè)出版社，2012.

［4］ Purcell W R.Capillary Pressures——Their Measurement Using Mercury and the Calculation of Permeability［C］∥ Trans AIME，1949，186：39－48.

［5］ 陶果，岳文正，謝然紅，等.巖石物理的理論模擬和數(shù)值實(shí)驗(yàn)新方法 ［J］.地球物理學(xué)進(jìn)展，2005，1（20）：4－11.

［6］ Fredrich J T，Menendez B，Wong T F.Imaging the Pore Structure of Geomaterials ［J］.Science，1995，268：276－279.

［7］ 趙秀才.數(shù)字巖心及孔隙網(wǎng)絡(luò)模型重構(gòu)方法研究［D］.青島：中國(guó)石油大學(xué)（華東），2009.

［8］ 姚軍，趙秀才，衣艷靜，等.數(shù)字巖心技術(shù)現(xiàn)狀及展望 ［J］.油氣地質(zhì)與采收率，2005，12（6）：52－54.

［9］ 孫建孟，姜黎明，劉學(xué)鋒，等.數(shù)字巖心技術(shù)測(cè)井應(yīng)用與展望 ［J］.測(cè)井技術(shù)，2012，36（1）：1－7.

［10］王晨晨，姚軍，楊永飛，等.基于CT掃描法構(gòu)建數(shù)字巖心的分辨率選取研究 ［J］.科學(xué)技術(shù)與工程，2013，13（4）：1050－1052.

［11］Jin Guodong，Carlos Torres－Verdin，Chun Lan.Porelevel Study of Grain－shape Effects on Petrophysical Properties of Porous Media［C］∥SPWLA，2009，Paper TTT.

［12］鄒才能.非常規(guī)油氣地質(zhì) ［M］.北京：地質(zhì)出版社，2011.

［13］雍世和，張超謨.測(cè)井?dāng)?shù)據(jù)處理與綜合解釋 ［M］.東營(yíng)：石油大學(xué)出版社，1996.

［14］歐陽(yáng)健，毛志強(qiáng)，修立軍，等.測(cè)井低對(duì)比度油層成因機(jī)理與評(píng)價(jià)方法 ［M］.北京：石油工業(yè)出版社，2009.

［15］中國(guó)石油勘探與生產(chǎn)分公司.低阻油氣藏測(cè)井評(píng)價(jià)技術(shù)及應(yīng)用 ［M］.北京：石油工業(yè)出版社，2009.