張鵬浩 張占松 朱林奇

摘 ? ? ?要： 對(duì)于孔隙結(jié)構(gòu)復(fù)雜的碳酸鹽巖儲(chǔ)層來(lái)說(shuō)，基于中高滲儲(chǔ)層提出的經(jīng)典核磁共振滲透率計(jì)算模型，不能精確的反映低滲儲(chǔ)層的孔隙結(jié)構(gòu)復(fù)雜性，適用性較差。因此，針對(duì)中東H油田低滲碳酸鹽巖儲(chǔ)層，提出了一種能夠反映可動(dòng)流體孔隙結(jié)構(gòu)的滲透率計(jì)算方法。根據(jù)壓汞曲線的分布特征，確立孔徑分類標(biāo)準(zhǔn)，然后通過(guò)建立核磁共振T2譜和壓汞孔喉半徑之間的關(guān)系，將孔徑分類標(biāo)準(zhǔn)轉(zhuǎn)化為核磁共振橫向弛豫時(shí)間標(biāo)準(zhǔn)，并將可動(dòng)流體孔隙精細(xì)的劃分為易流體孔隙和不易流體孔隙?；诓煌紫督M分含量對(duì)滲透率的不同影響，建立了精細(xì)評(píng)價(jià)孔隙結(jié)構(gòu)的核磁共振計(jì)算滲透率新模型。通過(guò)對(duì)比，新方法的計(jì)算精度明顯高于傳統(tǒng)經(jīng)典模型，對(duì)研究區(qū)塊有很好的適用性。

關(guān) ?鍵 ?詞：碳酸鹽巖;低滲儲(chǔ)層;核磁共振;壓汞;流體組分;滲透率

中圖分類號(hào)：TE344 ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼： A ? ? ? 文章編號(hào)： 1671-0460（2019）12-2861-05

Abstract： For carbonate reservoirs with complex pore structure， the classical nuclear magnetic resonance permeability calculation model for medium and high permeability reservoir analysis cannot accurately reflect the pore structure complexity of low permeability reservoirs， so its applicability is poor. Therefore， a permeability calculation method which can reflect the pore structure of movable fluids was proposed for the low permeability carbonate reservoir in the Middle East H oilfield. According to the distribution characteristics of mercury porosimetry curve， the pore size classification criterion was established. Then， the pore size classification criterion was transformed into the transverse relaxation time criterion of NMR by establishing the relationship between the T2 spectrum of NMR and the pore throat radius of mercury porosimetry. The movable fluid pore was divided into easy fluid pore and difficult fluid pore. Based on the different effect of different pore component contents on the permeability， a new permeability calculation model of nuclear magnetic resonance （NMR） for fine evaluation of pore structure was established. The comparison results showed that the accuracy of the new model was much higher than that of the traditional classical model， which had good applicability to the research blocks.

Key words： Carbonate rocks; Low permeability reservoir; Nuclear magnetic resonance; Mercury intrusion; Fluid component; Permeability

滲透率是儲(chǔ)層評(píng)價(jià)和儲(chǔ)量計(jì)算中的重要參數(shù)，滲透率的準(zhǔn)確計(jì)算是油田開發(fā)的關(guān)鍵一步[1-3]。對(duì)于常規(guī)儲(chǔ)層來(lái)說(shuō)，由于孔隙結(jié)構(gòu)簡(jiǎn)單，孔隙度和滲透率相關(guān)性較好，常規(guī)的測(cè)井方法就可以對(duì)滲透率進(jìn)行準(zhǔn)確的計(jì)算。但是對(duì)于巖性復(fù)雜的碳酸鹽巖儲(chǔ)層來(lái)說(shuō)，其孔隙結(jié)構(gòu)復(fù)雜多樣，常規(guī)的方法不能對(duì)滲透率進(jìn)行精確的預(yù)測(cè)。因此，在碳酸鹽巖儲(chǔ)層測(cè)井評(píng)價(jià)中，準(zhǔn)確預(yù)測(cè)滲透率是其中存在的一大難題[4-6]。

核磁共振測(cè)井作為一種能表征儲(chǔ)層孔隙結(jié)構(gòu)的測(cè)井方法，可以將儲(chǔ)集層孔隙中的束縛水體積直觀的表現(xiàn)出來(lái)[7-9]，從而對(duì)滲透率進(jìn)行精細(xì)的評(píng)價(jià)。很多學(xué)者提出了利用核磁共振測(cè)井技術(shù)計(jì)算滲透率的方法，其中最經(jīng)典的有Coates模型和SDR模型[10-12]，以及在此基礎(chǔ)上提出來(lái)的基于核磁共振T2譜分布的滲透率計(jì)算方法[12-14]。但該類模型在表征復(fù)雜孔隙結(jié)構(gòu)儲(chǔ)集層滲流機(jī)理方面，只是通過(guò)巖石宏觀信息來(lái)計(jì)算滲透率，沒(méi)有精細(xì)考慮T2譜的分布，不能充分反映不同孔隙結(jié)構(gòu)對(duì)滲透率的貢獻(xiàn)值不同的情況。尤其是對(duì)于低滲儲(chǔ)層來(lái)說(shuō)，相同的孔隙度，滲透率的大小會(huì)根據(jù)孔隙結(jié)構(gòu)而變化，表現(xiàn)出不同級(jí)別[15，16]。壓汞法繪制的毛管壓力曲線，可以通過(guò)進(jìn)汞過(guò)程表征儲(chǔ)層的孔隙結(jié)構(gòu)，因此，筆者結(jié)合巖心壓汞與核磁共振資料對(duì)滲透率進(jìn)行預(yù)測(cè)。通過(guò)分析壓汞孔喉半徑的分布，確立了小孔和微孔的分類標(biāo)準(zhǔn)，結(jié)合核磁共振T2譜的分布情況，建立了橫向弛豫時(shí)間與孔喉半徑的指數(shù)關(guān)系，將孔徑分類標(biāo)準(zhǔn)轉(zhuǎn)化成橫向弛豫時(shí)間標(biāo)準(zhǔn)，從而確定了不同孔隙組分的核磁共振橫向弛豫時(shí)間范圍，將可動(dòng)流體孔隙組分更精細(xì)的劃分為孔隙半徑較大，連通性較好，對(duì)滲透率起主要貢獻(xiàn)作用的易流體孔隙和孔隙半徑較小，連通性較差，對(duì)滲透率起限制作用的不易流體孔隙。通過(guò)計(jì)算易流體組分和不易流體組分的孔隙含量，并基于兩種孔隙組分和滲透率的關(guān)系，筆者提出了一種精細(xì)劃分孔隙組分的核磁共振計(jì)算滲透率的新方法，計(jì)算精度有了明顯的提高。

1 ?經(jīng)典核磁共振測(cè)井滲透率計(jì)算模型

（1）Coates模型，通過(guò)確定T2截止值，將孔隙分成可動(dòng)流體（VFFI）和束縛流體（VBVI），并將兩種流體孔隙的比值和孔隙度作為影響滲透率的兩大參數(shù)，進(jìn)行線性擬合，得到滲透率的計(jì)算公式：

Coates模型只是將孔隙度直接分為以可動(dòng)流體為主的大孔隙和以束縛流體為主的小孔隙，沒(méi)有考慮不同孔隙對(duì)滲透率的貢獻(xiàn)差異。SDR模型則是不考慮孔隙結(jié)構(gòu)的不同，通過(guò)將T2譜分布平均化，反映孔隙度的整體平均水平，來(lái)計(jì)算滲透率。對(duì)于孔隙結(jié)構(gòu)簡(jiǎn)單的常規(guī)儲(chǔ)層來(lái)說(shuō)，這兩類經(jīng)典模型通?？梢詼?zhǔn)確的求取滲透率。但是對(duì)于孔隙結(jié)構(gòu)復(fù)雜的碳酸巖儲(chǔ)層來(lái)說(shuō)，簡(jiǎn)單的流體劃分和幾何平均，并不能準(zhǔn)確反映出儲(chǔ)集層的儲(chǔ)集特征和滲流機(jī)理，計(jì)算效果不理想。以中東地區(qū)某區(qū)塊A井為例，采用Coates模型和SDR模型計(jì)算滲透率，并分別與核磁實(shí)驗(yàn)測(cè)得的巖心滲透率進(jìn)行對(duì)比，見圖1，可以看到對(duì)于滲透率大于10×10-3μm2的中高滲儲(chǔ)層樣品，兩個(gè)經(jīng)典模型的計(jì)算效果都比較好，但是對(duì)于那些滲透率小于10×10-3μm2的低滲樣品，計(jì)算效果并不理想。因此，對(duì)于滲透率小于10×10-3μm2的低滲樣品，筆者通過(guò)分析對(duì)應(yīng)的壓汞資料和核磁資料，對(duì)這些樣品的孔隙結(jié)構(gòu)進(jìn)行了更加深入的研究。

2 ?巖石孔隙結(jié)構(gòu)對(duì)儲(chǔ)層滲透率的影響

儲(chǔ)層孔隙結(jié)構(gòu)的好壞在孔隙度上會(huì)有對(duì)應(yīng)的體現(xiàn)，但是孔隙度不能作為表征孔隙結(jié)構(gòu)的決定因素。尤其在低滲復(fù)雜巖性儲(chǔ)層中，經(jīng)常出現(xiàn)孔隙度接近的樣品，但是滲透率屬于不同級(jí)別。如表1中，A1、A2是中東H油田兩塊典型的低滲樣品，在巖性和孔隙度接近的情況下，滲透率相差了兩個(gè)數(shù)量級(jí)。對(duì)兩塊樣品的壓汞、薄片等實(shí)驗(yàn)資料進(jìn)行對(duì)比分析，通過(guò)兩塊樣品的孔喉半徑分布特征可以看到，樣品A1的孔喉半徑峰值為0.04 μm，孔型主要是小于0.1 μm的微孔;樣品A2的孔喉半徑峰值為0.15 μm，孔型主要由大于0.1 μm的小孔組成。由樣品對(duì)應(yīng)的薄片資料也可以看出，樣品A1的孔隙類型主要是獨(dú)立存在的，孔隙半徑相對(duì)較小的微孔，不利于巖石孔隙中流體的相互流通，巖石的滲透性差。

樣品A2的孔隙類型主要是連通的，孔隙半徑相對(duì)較大的小孔，有利于巖石孔隙中流體的相互流通，提高了巖石的滲透性。

統(tǒng)計(jì)分析該研究區(qū)塊A井8塊低滲巖心資料的分析結(jié)果，均有上述類似特征。這說(shuō)明了孔隙結(jié)構(gòu)上的差異，是造成這些低滲透率巖心存在滲透率差異的主要原因。因此，有必要對(duì)孔隙結(jié)構(gòu)進(jìn)行精細(xì)研究，得到低滲透儲(chǔ)層的滲流特征，從而對(duì)滲透率進(jìn)行準(zhǔn)確評(píng)價(jià)。

在進(jìn)行核磁共振測(cè)井時(shí)，如果巖石處于水潤(rùn)濕并完全飽和鹽水的狀態(tài)下，那么經(jīng)反演處理得到的核磁T2譜形態(tài)，可以真實(shí)地反映出巖石的孔徑分布狀況[17]。

通過(guò)分析前面不同級(jí)別滲透率的孔喉半徑分布特征的差異，微孔和小孔的臨界孔喉半徑確定為0.1 μm，由公式（3）可得到對(duì)應(yīng)的核磁T2譜為35 ms。微孔和小孔的核磁共振橫向弛豫時(shí)間標(biāo)準(zhǔn)以及壓汞孔喉半徑分布標(biāo)準(zhǔn)見表2。微孔和小孔的準(zhǔn)確劃分，是準(zhǔn)確評(píng)價(jià)低滲透儲(chǔ)層滲透率的前提條件。

3 ?結(jié)合孔隙組分的核磁共振計(jì)算滲透率新模型

在低滲儲(chǔ)層劃分微孔和小孔的基礎(chǔ)上，對(duì)可動(dòng)流體進(jìn)行更精細(xì)的孔隙劃分，明確不同孔隙含量對(duì)滲透率的貢獻(xiàn)。根據(jù)巖心核磁T2譜分布特征的分析，對(duì)滲透率起主要貢獻(xiàn)作用的可動(dòng)流體進(jìn)行對(duì)應(yīng)的弛豫時(shí)間劃分。定義孔隙分量S為T2譜中不同弛豫時(shí)間對(duì)應(yīng)孔隙分量與總孔隙度分量的百分比，S1為可動(dòng)流體孔隙中，橫向弛豫時(shí)間大于35 ms所對(duì)應(yīng)的孔隙分量，稱為易流體孔隙;S2為可動(dòng)流體孔隙中，橫向弛豫時(shí)間小于35 ms所對(duì)應(yīng)的孔隙分量，稱為不易流體孔隙，見圖4。建立S1、S2與滲透率的關(guān)系，見圖5，S1和滲透率呈明顯的正相關(guān)，S1的含量越高，表明孔隙半徑越大，孔隙結(jié)構(gòu)越好，滲透率越大;S2和滲透率的關(guān)系呈負(fù)相關(guān)，S2的含量越高，表明孔隙半徑越小，孔隙結(jié)構(gòu)越差，滲透率越小。

通過(guò)回歸分析，本次研究所選用的參數(shù)分別為：a=0.46，b=-1.21，c=0.07，d=0.04，代入滲透率新模型，低滲巖心滲透率計(jì)算精度達(dá)到了0.943 3，遠(yuǎn)高于經(jīng)典模型Coates模型和SDR模型。

4 ?應(yīng)用效果

B井是該地區(qū)另一口進(jìn)行了巖心核磁實(shí)驗(yàn)的取心井，該井共有10塊低滲樣品。針對(duì)這些樣品的壓汞資料和薄片資料分析發(fā)現(xiàn)，孔隙類型中存在的孔隙半徑小于1 μm，以及連通性較差的孔隙，導(dǎo)致了巖石的滲透率偏小。

對(duì)于經(jīng)典的Coates模型，只是將孔隙類型劃分為可動(dòng)流體和束縛流體，對(duì)滲透率起主要貢獻(xiàn)作用的可動(dòng)流體中包含了孔隙結(jié)構(gòu)較差的孔隙，因此相對(duì)巖心滲透率，該類模型計(jì)算的滲透率往往偏大。而筆者建立的預(yù)測(cè)滲透率新模型，進(jìn)一步對(duì)可動(dòng)流體的孔隙類型進(jìn)行了精細(xì)的劃分，區(qū)分了對(duì)滲透率起主要貢獻(xiàn)作用的易流體孔隙，以及對(duì)滲透率起限制作用的不易流體孔隙。采用新模型對(duì)B井的低滲巖心進(jìn)行預(yù)測(cè)，預(yù)測(cè)結(jié)果見圖7，可以看到新模型的計(jì)算結(jié)果與巖心滲透率基本在45°對(duì)角線上，計(jì)算精度較高，相關(guān)性達(dá)到了0.947 3，說(shuō)明基于孔徑分布和核磁共振T2譜建立的新模型，通過(guò)對(duì)孔隙結(jié)構(gòu)進(jìn)行更加精細(xì)的劃分，對(duì)該地區(qū)低滲透油藏的評(píng)價(jià)有良好的指導(dǎo)作用。

5 ?解釋結(jié)論

（1）對(duì)于孔隙結(jié)構(gòu)和滲流機(jī)理復(fù)雜的碳酸鹽巖儲(chǔ)層來(lái)說(shuō)，基于中高滲儲(chǔ)層提出的經(jīng)典的核磁共振滲透率計(jì)算模型，不能精確表征低滲透滲儲(chǔ)層的孔隙結(jié)構(gòu)，計(jì)算精度較差。

（2）構(gòu)建壓汞孔喉分布半徑和核磁共振橫向弛豫時(shí)間的轉(zhuǎn)換關(guān)系，孔隙結(jié)構(gòu)通過(guò)核磁T2譜刻度，將可動(dòng)流體分為易流體孔隙和不易流體孔隙?；趦烧邔?duì)滲透率的貢獻(xiàn)效果，建立了一種能夠精細(xì)表征孔隙結(jié)構(gòu)，劃分孔隙分量的滲透率計(jì)算新模型，解決了低滲透儲(chǔ)層滲透率的計(jì)算難題。

（3）準(zhǔn)確劃分易流體孔隙和不易流體孔隙的界限，是新模型準(zhǔn)確評(píng)價(jià)滲透率的關(guān)鍵，需要基于一定數(shù)量的壓汞和核磁T2譜實(shí)驗(yàn)資料進(jìn)行統(tǒng)一分析，建立兩者之間的關(guān)系。

參考文獻(xiàn)：

[1]楊釗，王吉晨，萬(wàn)榮暉，高寒.儲(chǔ)層滲透率對(duì)天然氣開采的影響研究[J].當(dāng)代化工，2016，45（01）：34-36.

[2]韓玉嬌，周燦燦，范宜仁，等.基于孔徑組分的核磁共振測(cè)井滲透率計(jì)算新方法——以中東A油田生物碎屑灰?guī)r儲(chǔ)集層為例[J].石油勘探與開發(fā)，2018，45（01）：170-178.

[3]安小平，李相方，程時(shí)清，等.不同方法獲取滲透率的對(duì)比分析[J].油氣井測(cè)試， 2005，14（5）：14-17.

[4]申本科，薛大偉，趙君怡，等.碳酸鹽巖儲(chǔ)層常規(guī)測(cè)井評(píng)價(jià)方法[J].地球物理學(xué)進(jìn)展，2014，29（01）：261-270.

[5]郭振華，李光輝，吳蕾，等.碳酸鹽巖儲(chǔ)層孔隙結(jié)構(gòu)評(píng)價(jià)方法——以土庫(kù)曼斯坦阿姆河右岸氣田為例[J].石油學(xué)報(bào)，2011，32（03）：459-465.

[6]程媛，張沖，陳雨龍，等.基于壓汞資料的碳酸鹽巖儲(chǔ)層滲透率預(yù)測(cè)模型——以扎納若爾油田KT-Ⅰ和KT-Ⅱ含油層系灰?guī)r儲(chǔ)層為例[J].油氣地質(zhì)與采收率，2017，24（03）：10-17.

[7]劉堂宴，王紹民，傅容珊，等.核磁共振譜的巖石孔喉結(jié)構(gòu)分析[J].石油地球物理勘探，2003（03）：328-333.

[8]何雨丹，毛志強(qiáng)，肖立志，等.核磁共振T_2分布評(píng)價(jià)巖石孔徑分布的改進(jìn)方法[J].地球物理學(xué)報(bào)，2005（02）：373-378.

[9]蘇俊磊，孫建孟，王濤，等.應(yīng)用核磁共振測(cè)井資料評(píng)價(jià)儲(chǔ)層孔隙結(jié)構(gòu)的改進(jìn)方法[J].吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），2011，41（S1）：380-386.

[10]肖立志. 核磁共振成像測(cè)井與巖石核磁共振及其應(yīng)用[M]. 北京：科學(xué)出版社，1998：35-62.

[11]石油測(cè)井情報(bào)組編. 測(cè)井技術(shù)新應(yīng)用[M]. 北京：石油工業(yè)出版社，1998：122-150.

[12]周尚文，薛華慶，郭偉，等.川南龍馬溪組頁(yè)巖核磁滲透率新模型研究[J].中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2016，40（01）：56-61.

[13]盧文東，肖立志，李偉，等.核磁共振測(cè)井在低孔低滲儲(chǔ)層滲透率計(jì)算中的應(yīng)用[J].中國(guó)海上油氣，2007（02）：103-106.

[14]徐風(fēng)，白松濤，趙建斌，等.一種基于孔隙分量組合下的滲透率計(jì)算方法[J].石油天然氣學(xué)報(bào)，2013，35（11）：76-80+93+6-7.

[15]李潮流，徐秋貞，張振波.用核磁共振測(cè)井評(píng)價(jià)特低滲透砂巖儲(chǔ)層滲透性新方法[J].測(cè)井技術(shù)，2009，33（05）：436-439.

[16]付金華，石玉江.利用核磁測(cè)井精細(xì)評(píng)價(jià)低滲透砂巖氣層[J].天然氣工業(yè)， 2002 （06）：39-42+9-8.

[17]肖立志，謝然紅，丁娛嬌，等.核磁共振測(cè)井儀器的最新進(jìn)展與未來(lái)發(fā)展方向[J].測(cè)井技術(shù)，2003（04）：265-269.