基于孔徑組分的核磁共振測(cè)井滲透率計(jì)算新方法
——以中東A油田生物碎屑灰?guī)r儲(chǔ)集層為例

韓玉嬌，周燦燦，范宜仁，李潮流，袁超，叢云海

（1. 中國(guó)石油勘探開發(fā)研究院，北京 100083；2. 中國(guó)石油大學(xué)（華東）地球科學(xué)與技術(shù)學(xué)院，山東青島 266580；3. 中國(guó)石油集團(tuán)長(zhǎng)城鉆探工程有限公司地質(zhì)研究院，遼寧盤錦 124012）

0 引言

滲透率是確定儲(chǔ)集層產(chǎn)能和油藏評(píng)價(jià)的關(guān)鍵指標(biāo)之一，其準(zhǔn)確性關(guān)系到儲(chǔ)集層產(chǎn)能評(píng)價(jià)的準(zhǔn)確程度，但應(yīng)用測(cè)井資料準(zhǔn)確計(jì)算滲透率一直是測(cè)井界的難題[1-3]。尤其對(duì)于復(fù)雜巖性儲(chǔ)集層（巖性、孔隙流體、孔隙結(jié)構(gòu)復(fù)雜），其導(dǎo)電機(jī)理與常規(guī)儲(chǔ)集層差異較大，復(fù)雜的孔隙結(jié)構(gòu)對(duì)儲(chǔ)集層的儲(chǔ)集特征和滲流機(jī)理起主控作用，常規(guī)方法難以準(zhǔn)確計(jì)算其滲透率。核磁共振測(cè)井能夠反映儲(chǔ)集層的孔隙結(jié)構(gòu)和束縛水體積，從而提供精細(xì)的儲(chǔ)集層滲透率計(jì)算模型，如經(jīng)典的SDR模型和Coates模型[4-7]。該類模型的準(zhǔn)確性主要依賴于T2截止值參數(shù)的選取，但截至目前尚無(wú)方法能夠準(zhǔn)確計(jì)算T2截止值，且該模型在表征復(fù)雜孔隙結(jié)構(gòu)儲(chǔ)集層滲流機(jī)理方面有其局限性，導(dǎo)致計(jì)算的滲透率結(jié)果不可靠。近年來(lái)，一些學(xué)者致力于結(jié)合巖心壓汞與核磁共振測(cè)井資料預(yù)測(cè)儲(chǔ)集層滲透率[8-12]，也提出了對(duì)核磁共振T2譜不同組分進(jìn)行劃分的方法，但該方法沒(méi)有深入研究?jī)?chǔ)集層的孔隙結(jié)構(gòu)特征，劃分的依據(jù)不夠充分。

為解決利用核磁共振測(cè)井資料計(jì)算復(fù)雜巖性儲(chǔ)集層滲透率存在的諸多問(wèn)題，本文以生物碎屑灰?guī)r儲(chǔ)集層為例，從孔隙結(jié)構(gòu)評(píng)價(jià)入手，深入挖掘毛管壓力大小與孔喉類型的對(duì)應(yīng)關(guān)系，利用毛管壓力曲線形態(tài)及拐點(diǎn)位置確定了壓汞曲線孔徑分類標(biāo)準(zhǔn)；基于巖心刻度測(cè)井，建立了橫向弛豫時(shí)間與孔喉半徑的“兩段式”關(guān)系，從而確定了不同孔隙組分的核磁共振橫向弛豫時(shí)間范圍；在細(xì)分孔隙組分的基礎(chǔ)上，分析了各個(gè)孔隙分量對(duì)滲透率的貢獻(xiàn)，提出了一種基于核磁共振多組分孔隙分量組合的滲透率計(jì)算新模型。與傳統(tǒng)的SDR模型和Coates模型相比，本文建立的滲透率計(jì)算新模型能更精細(xì)地刻畫儲(chǔ)集層孔隙結(jié)構(gòu)，反映滲流機(jī)理，有效提高滲透率計(jì)算精度，為非均質(zhì)性復(fù)雜儲(chǔ)集層的精細(xì)解釋提供了新的研究思路。

1 復(fù)雜巖性儲(chǔ)集層滲透率主控因素

通常認(rèn)為，巖石的滲透率大小主要取決于其有效孔隙度，即滲透率受巖石骨架顆粒大小、粒徑分布（分選）、排列方向、顆粒充填方式及固結(jié)和膠結(jié)程度等因素影響。對(duì)于常規(guī)儲(chǔ)集層，滲透率與有效孔隙度一般具有較好的相關(guān)性；但在復(fù)雜巖性儲(chǔ)集層中，常出現(xiàn)孔隙度基本一致的巖石樣品，其滲透率有數(shù)量級(jí)上的差別[13]。如在中東 A油田，4塊典型生物碎屑灰?guī)r巖心樣品的薄片、壓汞等實(shí)驗(yàn)分析結(jié)果如表 1所示，在巖性一致、孔隙度相近的情況下，A、B以及C、D兩組樣品孔徑分布分別相似，孔喉半徑中值接近，滲透率也較為接近，而兩組巖心樣品之間孔徑分布及孔喉半徑中值差別較大，滲透率差異也很大。

表1 典型孔隙度相近巖心各分析資料對(duì)比

通過(guò)統(tǒng)計(jì)分析，該油田78塊生物碎屑灰?guī)r巖心分析結(jié)果均有相似特征，這說(shuō)明巖心滲透率的差異主要是由不同巖石樣品孔隙結(jié)構(gòu)上的差異造成的。由此可以得出結(jié)論，對(duì)于孔隙結(jié)構(gòu)多樣的復(fù)雜巖性儲(chǔ)集層，滲透率大小主要受控于孔隙度和孔隙結(jié)構(gòu)，僅用孔隙度參數(shù)不能全面客觀反映巖石的滲流特征，只有加強(qiáng)孔隙結(jié)構(gòu)的研究才能從機(jī)理上認(rèn)識(shí)復(fù)雜巖性儲(chǔ)集層的滲流機(jī)理，從而準(zhǔn)確求取滲透率。

孔隙度和孔隙結(jié)構(gòu)對(duì)滲透率的控制，主要源于孔喉大小及其占比。連通孔喉較粗且占比較高則滲透率較高，反之孔喉較細(xì)且占比較高則滲透率較小。不同大小孔喉對(duì)滲透率的貢獻(xiàn)可由（1）式確定[14]：

孔徑分布與核磁共振T2譜形態(tài)直接相關(guān)，而不同T2組分含量可看作是不同孔徑分量的多少，其累加值就是孔隙度，因此可以看出核磁共振技術(shù)是滲透率計(jì)算的有效方法。

2 經(jīng)典核磁共振測(cè)井滲透率計(jì)算模型的局限性

核磁共振測(cè)井計(jì)算滲透率的經(jīng)典模型主要有 SDR模型和Coates模型。SDR模型通過(guò)總孔隙度和T2幾何平均值來(lái)確定滲透率，Coates模型通過(guò)確定T2截止值，將孔隙分成大孔隙組分的可動(dòng)流體和小孔隙組分的束縛流體兩部分來(lái)建立滲透率計(jì)算模型（見(jiàn)圖 1）。SDR模型認(rèn)為雖然不同大小的孔徑（T2分布不同）對(duì)滲透率的貢獻(xiàn)不同，但將分布平均化可以計(jì)算樣品滲透率。Coates模型也認(rèn)為不同大小的孔徑（T2分布不同）對(duì)滲透率的貢獻(xiàn)不同，但簡(jiǎn)單地認(rèn)為滲透率只是大孔隙的可動(dòng)流體部分和小孔隙的束縛流體部分的兩段式貢獻(xiàn)作用。大量的常規(guī)儲(chǔ)集層為簡(jiǎn)單孔隙結(jié)構(gòu)，孔喉比不大，一般用SDR模型和Coates模型都能得到較為準(zhǔn)確的滲透率。對(duì)于復(fù)雜孔隙結(jié)構(gòu)儲(chǔ)集層，SDR模型和Coates模型計(jì)算滲透率誤差大。以中東A油田為例，圖2是分別利用SDR和Coates模型計(jì)算的滲透率與實(shí)驗(yàn)測(cè)定巖心絕對(duì)滲透率的交會(huì)圖，可以明顯看出計(jì)算誤差很大。壓汞法毛管壓力測(cè)試曲線（見(jiàn)圖3）證明儲(chǔ)集層孔喉分布較寬，不同大小的孔喉連續(xù)分布且對(duì)滲透率的貢獻(xiàn)差異較大，僅用SDR模型的幾何平均值或是Coates模型的“兩段式”貢獻(xiàn)均無(wú)法精確表征其復(fù)雜程度。由此可知，經(jīng)典的核磁共振測(cè)井滲透率模型之所以在復(fù)雜巖性儲(chǔ)集層計(jì)算誤差較大，其局限性在于它們對(duì)復(fù)雜儲(chǔ)集層T2分布所反映的孔隙結(jié)構(gòu)表征精度不夠。

基于以上討論，本文提出多組分滲透率計(jì)算新模型（見(jiàn)圖4），對(duì)孔隙結(jié)構(gòu)進(jìn)行更為精細(xì)的刻畫，以期能夠更好地表征不同孔隙尺寸分布對(duì)滲透率的貢獻(xiàn)，從而提高滲透率計(jì)算精度。

圖1 T2截止值確定的孔隙雙組分

圖2 核磁共振經(jīng)典模型計(jì)算滲透率與巖心滲透率擬合圖

圖3 壓汞法毛管壓力測(cè)試曲線

3 核磁共振多組分模型計(jì)算滲透率的思路和方法

3.1 壓汞孔徑分類標(biāo)準(zhǔn)的確定

利用中東A油田78塊生物碎屑灰?guī)r巖心的高壓壓汞測(cè)試數(shù)據(jù)（測(cè)試條件：溫度25.0 ℃，最高進(jìn)汞壓力235.43 MPa），繪制出潤(rùn)濕相飽和度和對(duì)應(yīng)毛管壓力的雙對(duì)數(shù)關(guān)系圖（見(jiàn)圖5a）。大量的毛管壓力曲線主要呈4類形態(tài)特征，對(duì)應(yīng)的薄片也證實(shí)這4類形態(tài)的孔隙結(jié)構(gòu)具有明顯差異（見(jiàn)圖5b），故將這批樣品的孔隙結(jié)構(gòu)劃分為4類。

毛管壓力曲線形態(tài)主要受孔隙連通程度及孔隙非均質(zhì)性的控制（見(jiàn)圖6），其中橫坐標(biāo)反映孔喉半徑的大小，縱坐標(biāo)為潤(rùn)濕相流體飽和度，反映孔喉的連通性，圖 6左上角區(qū)域?yàn)檫B通性較好的大孔喉，右下角為連通性較差的小孔喉。毛管壓力大小與孔喉半徑的關(guān)系式如下：

圖4 不同T2值確定的多組分孔隙示意圖

圖5 4種孔隙結(jié)構(gòu)巖心的毛管壓力曲線（a）及鑄體薄片（b）特征

圖6 毛管壓力曲線形態(tài)主控因素示意圖

圖7 4類孔隙結(jié)構(gòu)巖心孔徑分布特征

結(jié)合圖7不同孔隙結(jié)構(gòu)孔徑分布特征可以看出，Ⅰ類孔隙結(jié)構(gòu)對(duì)應(yīng)強(qiáng)溶蝕的成巖作用，孔喉分布均勻，排驅(qū)壓力最低，大孔喉最多，鑄體薄片顯示面孔率大，鑄模孔、粒間孔隙為主，連通性較好。Ⅱ類孔隙結(jié)構(gòu)以溶蝕成因?yàn)橹?，膠結(jié)作用亦有明顯影響，孔喉分布為兩段式，拐點(diǎn)左側(cè)粗、中孔喉占比相對(duì)較大且連通性較好，右側(cè)細(xì)、微孔喉相對(duì)較少。Ⅲ類孔隙結(jié)構(gòu)整體以膠結(jié)作用為主，溶蝕作用也有部分影響，孔喉分布亦呈兩段式，相對(duì)于Ⅱ類，拐點(diǎn)明顯靠左，排驅(qū)壓力明顯偏高，說(shuō)明拐點(diǎn)左側(cè)的粗、中孔喉占比較小，拐點(diǎn)右側(cè)的細(xì)、微孔喉分布連續(xù)且占比較大。鑄體薄片顯示面孔率較小，以孤立鑄?？住⒘?nèi)孔和微孔為主，連通性較差。Ⅳ類孔隙結(jié)構(gòu)以細(xì)、微孔喉為主，總體孔隙不發(fā)育，鑄體薄片視域內(nèi)基本為不可見(jiàn)孔隙。

結(jié)合大量巖心毛管壓力曲線的統(tǒng)計(jì)分析結(jié)果，Ⅳ類、Ⅲ類和Ⅱ類孔隙結(jié)構(gòu)的孔喉界限閾值可分別確定為4.90，0.74，0.15 MPa（分別為圖5a中Ⅳ類、Ⅲ類和Ⅱ類孔隙結(jié)構(gòu)的排驅(qū)壓力值），相應(yīng)的孔喉半徑采用毛管壓力公式分別計(jì)算為0.15，1.00，5.00 μm。再結(jié)合前面對(duì)毛管壓力曲線形態(tài)和拐點(diǎn)的分析，可將該地區(qū)儲(chǔ)集層的孔喉大小按照粗孔喉、中孔喉、細(xì)孔喉和微孔喉分成4類，標(biāo)準(zhǔn)如表2所示。

表2 孔喉分布分類標(biāo)準(zhǔn)

孔隙結(jié)構(gòu)從Ⅰ到Ⅳ類，物性逐漸變差，孔喉半徑逐漸變小，鑄體薄片面孔率越來(lái)越小，孔隙連通性越來(lái)越差，儲(chǔ)集空間類型逐漸由連通鑄?？紫蛄?nèi)孔、微孔過(guò)渡。4類孔隙結(jié)構(gòu)的參數(shù)特征如表3所示。

表3 4類孔隙結(jié)構(gòu)參數(shù)表

3.2 基于核磁共振T2譜的孔徑組分劃分

以上孔隙結(jié)構(gòu)的分類僅僅是依據(jù)實(shí)驗(yàn)室數(shù)據(jù)，若要進(jìn)行長(zhǎng)井段連續(xù)孔隙結(jié)構(gòu)類型劃分，必須依賴核磁共振測(cè)井技術(shù)。核磁共振理論認(rèn)為，當(dāng)孔喉中僅飽含水、磁場(chǎng)均勻且回波間隔足夠短時(shí)，T2與孔喉半徑之間具有如下關(guān)系：

可見(jiàn)，核磁共振橫向弛豫時(shí)間與孔喉半徑密切相關(guān)，若能夠確定形狀因子和橫向表面弛豫率，就可以準(zhǔn)確確定孔喉半徑。但橫向表面弛豫率至今還沒(méi)有一個(gè)簡(jiǎn)便易行的確定方法，形狀因子更是因?yàn)槎嗫捉橘|(zhì)形態(tài)的復(fù)雜性而難以簡(jiǎn)單地用球狀或柱狀描述，對(duì)于溶蝕作用強(qiáng)烈的生物碎屑灰?guī)r更是如此。但是本次研究發(fā)現(xiàn)，橫向表面弛豫率與形狀因子的乘積對(duì)于特定巖性的不同孔隙結(jié)構(gòu)類型具有統(tǒng)計(jì)規(guī)律，核磁共振橫向弛豫時(shí)間與孔喉半徑之間的關(guān)系，可以通過(guò)對(duì)實(shí)驗(yàn)樣品進(jìn)行配套的核磁共振和壓汞實(shí)驗(yàn)來(lái)獲得。

對(duì)比每塊樣品的核磁共振橫向弛豫時(shí)間譜和壓汞孔徑分布曲線可知，兩條曲線的總體形態(tài)和趨勢(shì)較為一致（見(jiàn)圖8）。選取對(duì)應(yīng)性較好的19塊樣品，將每塊巖心的T2譜和壓汞孔徑分布曲線的主峰一一對(duì)應(yīng)，其特征峰具有較好的T2-r對(duì)應(yīng)關(guān)系（見(jiàn)圖9），呈正相關(guān)的兩段式分布，說(shuō)明橫向表面弛豫率和形狀因子的乘積在研究區(qū)呈現(xiàn)為兩種類型。橫向弛豫時(shí)間和孔喉半徑r之間的關(guān)系可表示為：

圖8 核磁共振T2譜和壓汞孔徑分布曲線形態(tài)特征對(duì)比圖

圖9 T2-r冪指數(shù)關(guān)系圖

由前述分析知，不同孔隙結(jié)構(gòu)的臨界孔喉半徑值分別為 0.15，1.00和 5.00 μm，由公式（4）、（5）計(jì)算對(duì)應(yīng)的T2值分別為 30，90和 200 ms。不同大小孔隙的核磁共振橫向弛豫時(shí)間分類標(biāo)準(zhǔn)見(jiàn)表 4，據(jù)此標(biāo)準(zhǔn)可利用核磁共振測(cè)井資料將孔隙快速劃分為大、中、小和微孔隙 4個(gè)部分，為滲透率模型的建立奠定基礎(chǔ)。

3.3 多組分滲透率計(jì)算模型的建立

在細(xì)分孔隙組分的基礎(chǔ)上，分析各個(gè)孔隙分量占比對(duì)滲透率的貢獻(xiàn)（見(jiàn)圖10）：圖10a、圖10b的微孔占比和小孔占比均近似與滲透率成負(fù)相關(guān)關(guān)系，只是斜率明顯不同。這說(shuō)明在孔徑分布中，微孔隙和小孔隙所占比例越大，儲(chǔ)集層的滲透性就越差，但兩類孔徑區(qū)間對(duì)滲透率的負(fù)面影響程度不同；圖10c、圖10d的中孔占比、大孔占比均近似與滲透率成正相關(guān)，說(shuō)明中、大尺寸的孔隙組分越多，巖石的孔隙結(jié)構(gòu)越好，滲透性也越好。根據(jù)上述各個(gè)孔隙分量占比對(duì)滲透率貢獻(xiàn)的認(rèn)識(shí)，提出如下基于不同孔隙分量貢獻(xiàn)的滲透率計(jì)算新模型：

表4 不同孔徑核磁共振T2分類標(biāo)準(zhǔn)

實(shí)驗(yàn)室孔滲常規(guī)擬合計(jì)算模型如圖11所示，滲透率計(jì)算新模型和實(shí)驗(yàn)室孔滲常規(guī)擬合模型的計(jì)算結(jié)果與巖心滲透率的對(duì)比分別示于圖 12和圖 13，由這 3張圖可見(jiàn)，基于不同孔隙分量貢獻(xiàn)的滲透率計(jì)算模型的精度遠(yuǎn)高于常規(guī)擬合模型，也高于圖2所示的SDR模型和Coates模型，計(jì)算效果最好。

圖10 核磁共振樣品各類孔隙組分占比與滲透率關(guān)系圖

圖11 實(shí)驗(yàn)室孔滲常規(guī)擬合計(jì)算模型

圖12 滲透率新模型計(jì)算值與巖心滲透率對(duì)比

4 多組分滲透率計(jì)算模型在生物碎屑灰?guī)r儲(chǔ)集層的應(yīng)用

圖13 常規(guī)擬合模型計(jì)算值與巖心滲透率對(duì)比

對(duì)于實(shí)際核磁共振測(cè)井資料，考慮到油氣信號(hào)對(duì)T2譜影響較大且不可忽略，因此首先對(duì)實(shí)際核磁共振測(cè)井資料利用重構(gòu)法進(jìn)行油氣影響校正[15]。具體實(shí)現(xiàn)步驟如下：①結(jié)合工區(qū)實(shí)際情況選取T2截止值，確定核磁共振測(cè)井T2譜的束縛水飽和度Swi；②利用T2截止值將實(shí)測(cè)核磁共振T2譜劃分為束縛水譜和可動(dòng)流體譜；③利用 100%飽和水實(shí)驗(yàn)條件下的巖心核磁共振T2譜和孔隙度測(cè)量值確定束縛水飽和度和T2幾何平均值T2lm，并建立Swi和T2lm的關(guān)系；④根據(jù)步驟①得到的Swi，確定Sw=100%情況下T2的最佳T2lm；⑤采用FFI波形函數(shù)構(gòu)造新的可動(dòng)流體譜，使可動(dòng)流體譜與束縛水譜幅度的累加值等于總孔隙度，調(diào)整其T2譜的位置使整體譜的T2幾何平均值等于T2lm；⑥原始束縛水譜加上步驟⑤構(gòu)造的可動(dòng)流體譜即為最終的油氣校正之后的T2譜。

以 X1井（見(jiàn)圖 14）為例，微孔、小孔、中孔和大孔孔隙組分占比在第 8～11道，基于這4道貢獻(xiàn)的滲透率新模型計(jì)算結(jié)果放在最后一道。為了對(duì)比，最后一道還同時(shí)放入了實(shí)驗(yàn)室滲透率常規(guī)擬合模型和SDR模型、Coates模型計(jì)算的滲透率。對(duì)比可知，常規(guī)滲透率模型（綠線）計(jì)算結(jié)果精度最差，高滲部分計(jì)算結(jié)果明顯偏小，低滲部分又偏大；本文新模型滲透率計(jì)算效果最好（藍(lán)線），其計(jì)算精度也明顯優(yōu)于SDR模型（黃線）和Coates模型（紅線）。

圖14 X1井滲透率計(jì)算效果對(duì)比圖

5 結(jié)論

常規(guī)儲(chǔ)集層孔隙結(jié)構(gòu)簡(jiǎn)單，孔喉比變化不大；但對(duì)于復(fù)雜孔隙結(jié)構(gòu)儲(chǔ)集層，孔隙度和孔隙結(jié)構(gòu)對(duì)儲(chǔ)集特征和滲流規(guī)律有控制性作用，孔徑分布及不同孔徑分量占比決定了最終的滲透率，準(zhǔn)確表征孔隙結(jié)構(gòu)有助于從機(jī)理上認(rèn)識(shí)復(fù)雜巖性儲(chǔ)集層的滲流機(jī)制。

中東A油田生物碎屑灰?guī)r，從孔隙結(jié)構(gòu)角度可分為 4種類型，從Ⅰ到Ⅳ類，其物性逐漸變差，孔喉半徑逐漸變小，面孔率越來(lái)越小，孔隙連通性越來(lái)越差，儲(chǔ)集空間類型也逐漸由連通鑄模孔向粒內(nèi)孔、微孔過(guò)渡?？紫督Y(jié)構(gòu)可由核磁共振測(cè)井刻畫，核磁共振橫向弛豫時(shí)間與孔喉半徑的對(duì)應(yīng)關(guān)系可由巖心樣品的核磁共振與壓汞配套實(shí)驗(yàn)獲得。橫向弛豫時(shí)間與孔喉半徑呈明顯的兩段式關(guān)系，必須采用兩個(gè)方程分段計(jì)算。多組分滲透率計(jì)算新模型基于不同類型孔隙分量對(duì)滲透率的不同貢獻(xiàn)，對(duì)復(fù)雜孔隙結(jié)構(gòu)的微、小、中和大孔喉組分精細(xì)定量刻畫，能更加合理地表征不同孔隙尺寸分布對(duì)滲透率的貢獻(xiàn)差異，可顯著提高滲透率計(jì)算精度。

符號(hào)注釋：

a，b，c，d，e，f——滲透率新模型相關(guān)參數(shù)，可利用樣本點(diǎn)采用優(yōu)化算法迭代獲得，本次研究所用參數(shù)分別為1.46、2.98、0.37、3.27、0.62 和 1.04；dh——井徑，cm；Fs——孔隙形狀因子，無(wú)因次；GR——自然伽馬，API；K——滲透率，10-3μm2；Kc——常規(guī)滲透率模型計(jì)算結(jié)果，10-3μm2；Kn——新模型滲透率計(jì)算結(jié)果，10-3μm2；Kr——實(shí)驗(yàn)室測(cè)量巖心滲透率，10-3μm2；KSDR——核磁共振SDR模型滲透率計(jì)算結(jié)果，10-3μm2；KTIM——核磁共振 Coates模型滲透率計(jì)算結(jié)果，10-3μm2；pc——毛管壓力，MPa；r——孔喉半徑，μm；rj——第j個(gè)區(qū)間的孔喉半徑值，μm；R——相關(guān)系數(shù)；RLLD——深側(cè)向電阻率，?·m；RLLS——淺側(cè)向電阻率，?·m；RMSFL——微球電阻率，?·m；S——潤(rùn)濕相飽和度，%；SP——自然電位，mV；S1，S2，S3，S4——微孔、小孔、中孔、大孔孔隙組分占比，%；SHg——汞飽和度，%；Sw——含水飽和度，%；Swi——束縛水飽和度，%；T2——核磁共振測(cè)井記錄的橫向弛豫時(shí)間，ms；T2lm——T2幾何平均值，ms；αj——第j個(gè)區(qū)間的孔喉頻率，%；ΔKj——第j個(gè)區(qū)間的滲透率貢獻(xiàn)，無(wú)因次；Δt——聲波時(shí)差，μs/m；θ——汞與巖石的潤(rùn)濕角，本文取127°；ρ——核磁共振橫向表面弛豫率，μs/ms；ρr——密度，g/cm3；σ——汞與空氣的界面張力，本文取483×10-3N/m；φ——孔隙度，%；φc——實(shí)驗(yàn)室?guī)r心分析孔隙度，%；φCNL——中子孔隙度，%；φl(shuí)——測(cè)井孔隙度，%。

[1]CHEHRAZI A, REZAEE R. A systematic method for permeability prediction: A petro-facies approach[J]. Journal of Petroleum Science and Engineering, 2012, 82/83(2): 1-16.

[2]王新民, 郭彥如, 付金華, 等. 鄂爾多斯盆地延長(zhǎng)組長(zhǎng)8段相對(duì)高孔滲砂巖儲(chǔ)集層的控制因素分析[J]. 石油勘探與開發(fā), 2005, 32(2):35-38.WANG Xinmin, GUO Yanru, FU Jinhua, et al. Control factors for forming higher porosity and permeability sand stone reservoirs in Chang 8 member of Yanchang Formation, Ordos Basin[J]. Petroleum Exploration and Development, 2005, 32(2): 35-38.

[3]肖佃師, 盧雙舫, 陸正元, 等. 聯(lián)合核磁共振和恒速壓汞方法測(cè)定致密砂巖孔喉結(jié)構(gòu)[J]. 石油勘探與開發(fā), 2016, 43(6): 961-970.XIAO Dianshi, LU Shuangfang, LU Zhengyuan, et al. Combining nuclear magnetic resonance and rate-controlled porosimetry to probe the pore-throat structure of tight sandstones[J]. Petroleum Exploration and Development, 2016, 43(6): 961-970.

[4]WEI D F, LIU X P, HU X X, et al. Estimation of permeability from NMR logs based on formation classification method in tight gas sands[J]. Acta Geophysica, 2015, 63(5): 1316-1338.

[5]DUNN K J, BERGMAN D J, LATORRACA G A. Nuclear magnetic resonance: Petrophysical and logging application[M]. Amsterdam:Elsevier, 2002.

[6]COATES G R, XIAO L Z, PRIMMER M G. NMR logging principles and applications[M]. Houston: Gulf Publishing Company, 2000.

[7]肖立志. 核磁共振成像測(cè)井與巖石核磁共振及其應(yīng)用[M]. 北京:科學(xué)出版社, 1998.XIAO Lizhi. The technology and application of magnetic resonance imaging logging and rock NMR[M]. Beijing: Science Press, 1998.

[8]肖亮, 劉曉鵬, 毛志強(qiáng). 結(jié)合NMR和毛管壓力資料計(jì)算儲(chǔ)層滲透率的方法[J]. 石油學(xué)報(bào), 2009, 30(1): 100-103.XIAO Liang, LIU Xiaopeng, MAO Zhiqiang. A computation method for reservoir permeability by combining NMR log and capillary pressure data[J]. Acta Petrolei Sinica, 2009, 30(1): 100-103.

[9]張麗華, 左義僑, 鞠東陽(yáng), 等. 用核磁 T2譜數(shù)據(jù)計(jì)算復(fù)雜孔隙結(jié)構(gòu)巖石的滲透率[J]. 國(guó)外測(cè)井技術(shù), 2017, 38(1): 10-12.ZHANG Lihua, ZUO Yiqiao, JU Dongyang, et al. Calculation of permeability of rock with complex pore structure by using NMR T2 spectrum data[J]. World Well Logging Technology, 2017, 38(1): 10-12.

[10]徐風(fēng), 白松濤, 趙建斌, 等. 一種基于孔隙分量組合下的滲透率計(jì)算方法[J]. 石油天然氣學(xué)報(bào), 2013, 35(11): 76-80, 93.XU Feng, BAI Songtao, ZHAO Jianbin, et al. A method of calculating permeability based on combination of pore components[J]. Journal of Oil and Gas Technology, 2013, 35(11): 76-80, 93.

[11]肖忠祥, 肖亮. 基于核磁共振測(cè)井和毛管壓力的儲(chǔ)層滲透率計(jì)算方法[J]. 原子能科學(xué)技術(shù), 2008, 42(10): 868-871.XIAO Zhongxiang, XIAO Liang. Method to calculate reservoir permeability using nuclear magnetic resonance logging and capillary pressure data[J]. Atomic Energy Science and Technology, 2008,42(10): 868-871.

[12]盧文東, 肖立志, 李偉, 等. 核磁共振測(cè)井在低孔低滲儲(chǔ)層滲透率計(jì)算中的應(yīng)用[J]. 中國(guó)海上油氣, 2007, 19(2): 103-106.LU Wendong, XIAO Lizhi, LI Wei, et al. An application of NMR logging to calculating permeability in low porosity and low permeability reservoir[J]. China Offshore and Gas, 2007, 19(2):103-106.

[13]TIAB D, DONALDSON E G. Petrophysics: Theory and practice of measuring reservoir rock and transport properties[M]. Houston: Gulf Professional Publishing, 2015.

[14]馮程, 石玉江, 郝建飛, 等. 低滲透復(fù)雜潤(rùn)濕性儲(chǔ)集層核磁共振特征[J]. 石油勘探與開發(fā), 2017, 44(2): 252-257.FENG Cheng, SHI Yujiang, HAO Jianfei, et al. Nuclear magnetic resonance features of low-permeability reservoirs with complex wettability[J]. Petroleum Exploration and Development, 2017, 44(2):252-257.

[15]HOFMAN J, SLIJKERMAIL W, LOOYESTIJN W, et al.Constructing capillary pressure curve from NMR log data in the presence of hydrocarbons[R]. Oslo, Norway: SPWLA 40thAnnual Logging Symposium, 1999.