王 迪，魏 鋒，陳 現(xiàn)，陳 敏，楊 冬

改進(jìn)的西門度模型在飽和度測(cè)井解釋中的應(yīng)用

王 迪1，魏 鋒1，陳 現(xiàn)1，陳 敏1，楊 冬2

（1. 中國石油化工股份有限公司上海海洋油氣分公司勘探開發(fā)研究院，上海 200120；2.中海石油（中國）有限公司湛江分公司，廣東湛江 524057）

由于西門度公式、印度尼西亞公式等具有泥質(zhì)校正的功能，故在含泥質(zhì)砂巖的飽和度測(cè)井解釋中得到了較為廣泛的運(yùn)用。通過研究發(fā)現(xiàn)，對(duì)于西門度公式來說，雖然通過泥質(zhì)校正可以在一定程度上消除泥質(zhì)附加導(dǎo)電作用對(duì)地層電阻率的影響，但在實(shí)際應(yīng)用過程中，由于模型假設(shè)條件和參數(shù)選取等原因，飽和度計(jì)算的精度受到限制：一方面，西門度公式未對(duì)泥質(zhì)中細(xì)粉砂和黏土的附加導(dǎo)電強(qiáng)弱進(jìn)行區(qū)分；另一方面，巖電參數(shù)往往使用的是經(jīng)驗(yàn)值或直接由巖電實(shí)驗(yàn)得到的值，并非為消除了泥質(zhì)導(dǎo)電影響的巖電參數(shù)值。為解決上述問題，從突出黏土在泥質(zhì)附加導(dǎo)電中的重要性出發(fā)，對(duì)原西門度模型進(jìn)行了重新推導(dǎo)，并在此基礎(chǔ)上，總結(jié)出運(yùn)用GR曲線估算黏土含量的方法，同時(shí)，利用巖電實(shí)驗(yàn)數(shù)據(jù)，運(yùn)用多元回歸模型擬合得到消除了泥質(zhì)導(dǎo)電影響的巖電參數(shù)值。以東海陸架盆地X氣田為例對(duì)改進(jìn)的西門度公式進(jìn)行了驗(yàn)證，飽和度計(jì)算精度明顯高于常規(guī)西門度公式的精度，證實(shí)了這一模型的有效性。

西門度公式；黏土含量；巖電參數(shù)；體積物理模型

1942年，Archie等人在實(shí)驗(yàn)的基礎(chǔ)上提出了阿爾奇公式[1-2]，建立了含水飽和度和地層電阻率的相互轉(zhuǎn)換關(guān)系，用于飽和度測(cè)井解釋。由于Archie在實(shí)驗(yàn)中使用的是物性好的純砂巖[3]，故阿爾奇公式在該類儲(chǔ)層中獲得較好的應(yīng)用效果，而在一些低孔低滲含泥質(zhì)砂巖儲(chǔ)層中的應(yīng)用效果并不理想。為實(shí)現(xiàn)泥質(zhì)附加導(dǎo)電校正功能，分散泥質(zhì)砂巖模型[4-5]、層狀泥質(zhì)砂巖模型[6]、混合泥質(zhì)砂巖模型（如西門度公式、印度尼西亞公式）[7]，陽離子交換模型（如W-S模型和雙水模型）等[8-9]飽和度測(cè)井解釋模型應(yīng)運(yùn)而生，其中西門度公式在包括東海陸架盆地在內(nèi)的許多地區(qū)含泥質(zhì)砂巖儲(chǔ)層的飽和度計(jì)算中都得到了廣泛的應(yīng)用。對(duì)于東海陸架盆地X氣田來說，由于地層水礦化度較低（8 000 ~10 000 mg/L），黏土礦物中伊利石—伊蒙混層含量較高，泥質(zhì)附加導(dǎo)電能力對(duì)地層導(dǎo)電性的貢獻(xiàn)不能忽略，因此一直以來，具泥質(zhì)附加導(dǎo)電校正功能的西門度公式在該地區(qū)取得了較好的應(yīng)用效果。近年來，在密閉取心技術(shù)得到廣泛應(yīng)用之后，發(fā)現(xiàn)在運(yùn)用西門度公式對(duì)東海陸架盆地部分油氣田的部分層段進(jìn)行飽和度解釋時(shí)，得到的含水飽和度與密閉取心飽和度差異較為明顯，因此，需要從常規(guī)西門度公式的假設(shè)條件出發(fā)，對(duì)當(dāng)前模型和所取參數(shù)進(jìn)行改進(jìn)和校正。

1 西門度模型的改進(jìn)——并聯(lián)導(dǎo)電模型的優(yōu)化

1.1 常規(guī)西門度公式的推導(dǎo)原理

西門度公式為混合泥質(zhì)砂巖的導(dǎo)電模型，即不考慮儲(chǔ)層中泥質(zhì)的分布形式，只是把泥質(zhì)近似看作是巖性非常細(xì)、可以儲(chǔ)存油氣、符合阿爾奇公式的“粉砂巖”，認(rèn)為泥質(zhì)砂巖的導(dǎo)電性是由這部分“粉砂巖”和純砂巖并聯(lián)導(dǎo)電形成的[7]。由并聯(lián)導(dǎo)電機(jī)理可得巖石的電阻率與泥質(zhì)質(zhì)量分?jǐn)?shù)的關(guān)系為：

式中：Rt為地層電阻率，Ω·m；Vsh為泥質(zhì)質(zhì)量分?jǐn)?shù)，0~1；Rs為含油氣黏土和細(xì)粉砂電阻率，Ω·m；Rsd為純砂巖部分電阻率，Ω·m；

對(duì)于純砂巖部分，其孔隙度：

由阿爾奇公式可得

式中： Sw為含水飽和度，0~1；φe為有效孔隙度，0~1；φsd為純砂巖部分孔隙度，0~1；Rw為地層水電阻率，Ω·m。

對(duì)于泥質(zhì)部分，認(rèn)為其也具有一定的儲(chǔ)存油氣的性能，并近似認(rèn)為這部分細(xì)巖性有效孔隙中含水飽和度等于純巖石孔隙中的含水飽和度[7]，且飽和度指數(shù)n=1，故有：

式中： Rsh為完全含水泥質(zhì)的電阻率，Ω·m。

聯(lián)立以上兩式，得：

為便于計(jì)算機(jī)連續(xù)解釋，上述方程式可寫作：

這就是西門度方程的表達(dá)式。

1.2 基于黏土附加導(dǎo)電的西門度模型優(yōu)化

從上述推導(dǎo)過程中可以看出，西門度方程的得出基于兩個(gè)關(guān)鍵的假設(shè)條件。一是認(rèn)為泥質(zhì)的導(dǎo)電能力是不受其粉砂指數(shù)（泥質(zhì)中細(xì)粉砂所占比例）影響，即對(duì)于相同體積完全含水的泥質(zhì)，不管其中細(xì)粉砂和黏土的比例如何，其導(dǎo)電能力都是相同的。但事實(shí)上，從陽離子交換的機(jī)理可知，泥質(zhì)的附加導(dǎo)電主要由黏土的陽離子交換作用造成，而細(xì)粉砂的陽離子交換能力有限，因此對(duì)于粉砂指數(shù)不同的泥質(zhì)，附加導(dǎo)電能力應(yīng)有較大差別[7]。由于未關(guān)注泥質(zhì)的粉砂指數(shù)，這同時(shí)會(huì)造成解釋時(shí)Rsh值不易準(zhǔn)確讀取。另一方面，認(rèn)為泥質(zhì)和純巖石部分孔隙的含水飽和度相同，但實(shí)際上，由于泥質(zhì)顆粒很細(xì)，其束縛水飽和度會(huì)高于純巖石部分，且對(duì)于親水性較強(qiáng)或本地區(qū)低孔低滲且非均質(zhì)性較強(qiáng)的巖石來說，這種差異會(huì)更加明顯。

為解決上述問題，厘清“泥質(zhì)導(dǎo)電”這一相對(duì)模糊的概念，對(duì)并聯(lián)導(dǎo)電模型進(jìn)行優(yōu)化，把“純巖石”和“泥質(zhì)”并聯(lián)導(dǎo)電模型優(yōu)化為“孔隙流體自由離子”和“黏土陽離子交換”并聯(lián)導(dǎo)電模型，即將“細(xì)粉砂”部分從泥質(zhì)中歸并至巖石骨架中，認(rèn)為計(jì)算得到的Sw為測(cè)井儀器探測(cè)范圍內(nèi)，由粗細(xì)不同砂粒（包括細(xì)粉砂）構(gòu)成的巖石骨架的孔隙中的總含水飽和度，并認(rèn)為極細(xì)的黏土顆粒間對(duì)有效孔隙體積幾乎沒有貢獻(xiàn)的微孔隙不具有儲(chǔ)存和滲流油氣的能力，即Swcl= 1。故西門度方程可改寫為：

式中：Rcl為黏土電阻率，Ω·m；Vcl為黏土質(zhì)量分?jǐn)?shù)，0~1。

需要注意的是，和泥質(zhì)（含細(xì)粉砂）存在多種分布形式（層狀、分散、結(jié)構(gòu)泥質(zhì)）不同，黏土礦物在儲(chǔ)層中基本上均以分散的形式存在于孔隙壁（綠泥石、伊利石）和孔隙空間（伊利石、高嶺石）中，故考慮用分散泥質(zhì)砂巖模型對(duì)西門度公式進(jìn)行改進(jìn)。針對(duì)這一問題，可以參考Hossin的分散泥質(zhì)砂巖模型[4]，認(rèn)為泥質(zhì)砂巖是由“分散的黏土”和“砂巖骨架及骨架中含自由離子的地層水”并聯(lián)導(dǎo)電的。對(duì)于孔隙流體中自由離子導(dǎo)電部分，由阿爾奇公式可得：

將式（8）與式（3）對(duì)比，發(fā)現(xiàn)等號(hào)右邊分子略有差異，其原因在于混合泥質(zhì)砂巖模型認(rèn)為純砂巖部分的孔隙度（含自由離子的地層水所占孔隙度）為sdφ （計(jì)算式如式2），而分散泥質(zhì)砂巖模型中含自由離子的地層水所占孔隙度為eφ。聯(lián)立式（7）和式（8），得到經(jīng)過這一步優(yōu)化的西門度公式為：

和常規(guī)西門度公式讀取Rsh的方式略有不同，為體現(xiàn)泥質(zhì)中黏土的導(dǎo)電能力，Rcl要在全井段自然伽馬值盡可能最高的穩(wěn)定純泥巖處讀取，或通過黏土礦物分析資料得到。

1.3 泥質(zhì)中黏土含量的求取

為實(shí)現(xiàn)含水飽和度的計(jì)算，需要先用測(cè)井資料分別求出泥質(zhì)中的黏土質(zhì)量分?jǐn)?shù)和細(xì)粉砂質(zhì)量分?jǐn)?shù)。一般來說，核磁共振測(cè)井和常規(guī)三孔隙度測(cè)井可以用于黏土質(zhì)量分?jǐn)?shù)的計(jì)算[10-13]，但這類方法往往難以達(dá)到較高的精度：一方面，儲(chǔ)層溫壓條件下不同黏土礦物的骨架值不同且難以確定，另一方面，儲(chǔ)層中質(zhì)量分?jǐn)?shù)本身就不高的黏土對(duì)孔隙度測(cè)井響應(yīng)的貢獻(xiàn)程度很低，且三孔隙度、核磁共振數(shù)據(jù)易受井眼及地層含氣性影響，這使得孔隙度曲線對(duì)黏土質(zhì)量分?jǐn)?shù)敏感程度降低，再者，當(dāng)?shù)貙又叙ね恋V物種類較多且分布呈明顯非均質(zhì)性時(shí)，若將地層組分分得過細(xì)，可能無法構(gòu)成超定方程組進(jìn)行最優(yōu)化求解。包含地球化學(xué)信息的自然伽馬能譜測(cè)井和元素俘獲測(cè)井（ECS）等方法在計(jì)算黏土質(zhì)量分?jǐn)?shù)方面具有一定的優(yōu)勢(shì)[14-17]，但對(duì)于東海陸架盆地X氣田來說，這兩種測(cè)井系列使用較少，用于建模的資料不足。

相對(duì)于三孔隙度曲線，自然伽馬曲線受井眼影響相對(duì)較小，且對(duì)泥質(zhì)（包括黏土）響應(yīng)靈敏，且和自然伽馬能譜、ECS相比，自然伽馬曲線數(shù)據(jù)明顯更容易收集。同時(shí)，粒度分析資料可以提供準(zhǔn)確的細(xì)粉砂和黏土質(zhì)量分?jǐn)?shù)。因此，在這里提出一種運(yùn)用自然伽馬曲線和粒度分析資料估算黏土質(zhì)量分?jǐn)?shù)的方法：先用GR曲線通過常規(guī)方法計(jì)算地層總泥質(zhì)質(zhì)量分?jǐn)?shù)，然后根據(jù)粒度分析資料，分別刻度黏土和細(xì)粉砂對(duì)測(cè)井泥質(zhì)質(zhì)量分?jǐn)?shù)貢獻(xiàn)程度的權(quán)重以及黏土和細(xì)粉砂的比例關(guān)系，建立刻度模型，最后再用自然伽馬曲線和上述刻度系數(shù)估算出儲(chǔ)層中黏土和細(xì)粉砂質(zhì)量分?jǐn)?shù)。

由于細(xì)粉砂和黏土的粒徑不同，因此相同體積的細(xì)粉砂和黏土，對(duì)地層放射性的貢獻(xiàn)程度是不同的。在這里，認(rèn)為測(cè)井解釋得到的泥質(zhì)質(zhì)量分?jǐn)?shù)是地層中黏土和細(xì)粉砂體積的加權(quán)和[18]。由于顆粒越細(xì)，吸附放射性物質(zhì)的能力越強(qiáng)，故顆粒相對(duì)較粗的細(xì)粉砂對(duì)測(cè)井泥質(zhì)質(zhì)量分?jǐn)?shù)貢獻(xiàn)的權(quán)重小于黏土的權(quán)重。為方便計(jì)算，這里不妨設(shè)單位體積黏土對(duì)測(cè)井泥質(zhì)質(zhì)量分?jǐn)?shù)的貢獻(xiàn)權(quán)重為1，單位體積細(xì)粉砂對(duì)測(cè)井泥質(zhì)質(zhì)量分?jǐn)?shù)貢獻(xiàn)權(quán)重為a（a＜1）。則巖石泥質(zhì)質(zhì)量分?jǐn)?shù)的加權(quán)計(jì)算公式為：

式中：Vsh為測(cè)井解釋泥質(zhì)質(zhì)量分?jǐn)?shù)，%；Vfss為細(xì)粉砂質(zhì)量分?jǐn)?shù)，%；Vcl為黏土質(zhì)量分?jǐn)?shù)，%；a為細(xì)粉砂對(duì)泥質(zhì)質(zhì)量分?jǐn)?shù)貢獻(xiàn)的歸一化權(quán)重，無量綱；b為常數(shù)，即巖石中其他粒度成分對(duì)泥質(zhì)質(zhì)量分?jǐn)?shù)測(cè)井響應(yīng)值的貢獻(xiàn)，%。

分別代入粒度分析得到的黏土質(zhì)量分?jǐn)?shù)Vcl、細(xì)粉砂質(zhì)量分?jǐn)?shù)Vfss和測(cè)井解釋得到的泥質(zhì)質(zhì)量分?jǐn)?shù)Vsh，擬合得到系數(shù)a、b分別為0.924 5和1.550 7（擬合圖為圖1a）。

對(duì)于相同層位的砂體，沉積環(huán)境類似，當(dāng)儲(chǔ)層中泥質(zhì)質(zhì)量分?jǐn)?shù)總體變化不大時(shí)，細(xì)粉砂質(zhì)量分?jǐn)?shù)和黏土質(zhì)量分?jǐn)?shù)可能存在穩(wěn)定的函數(shù)關(guān)系[18]。對(duì)X氣田其中一個(gè)目的層段Y3段的粒度分析資料中細(xì)粉砂和黏土的質(zhì)量分?jǐn)?shù)作分析，發(fā)現(xiàn)兩者函數(shù)關(guān)系非常好（圖1b），滿足以下線性關(guān)系式：

圖1 刻度系數(shù)擬合圖

由式（11）得粉砂指數(shù)SI的計(jì)算式為：

聯(lián)立式（10）和（11），可以估算得到Y(jié)3段儲(chǔ)層中的黏土質(zhì)量分?jǐn)?shù)。需要說明的是，細(xì)粉砂和黏土質(zhì)量分?jǐn)?shù)的函數(shù)關(guān)系存在較強(qiáng)的地區(qū)經(jīng)驗(yàn)性，不同地區(qū)甚至不同層位，兩者的比例及相關(guān)系數(shù)可能存在明顯差異。由式（12）可計(jì)算得出，當(dāng)黏土質(zhì)量分?jǐn)?shù)Vcl小于4%時(shí)，該目的層段粉砂指數(shù)始終大于77.5%，即Vcl明顯小于Vsh，因此改進(jìn)的西門度模型可以在一定程度上避免因代入Vsh計(jì)算而往往得到的Sw明顯偏低的結(jié)果[7]。

2 西門度模型的改進(jìn)——m*和n*的定義與求取

2.1 求取消除黏土導(dǎo)電影響的巖電參數(shù)的必要性

巖電實(shí)驗(yàn)的目的，在于借助對(duì)巖心進(jìn)行實(shí)驗(yàn)室測(cè)量來模擬實(shí)際地層，通過對(duì)m、n等參數(shù)的擬合，來得到Rt和Sw的相互轉(zhuǎn)換關(guān)系。擬合公式如下：

巖電實(shí)驗(yàn)所選取的巖心，多多少少具有一定的黏土質(zhì)量分?jǐn)?shù)，黏土的存在會(huì)對(duì)實(shí)驗(yàn)測(cè)量的Rt產(chǎn)生一定的影響，而阿爾奇公式是將含黏土砂巖當(dāng)作純砂巖來處理的，因此通過阿爾奇模型（式13、式14）擬合得到的巖電參數(shù)，實(shí)際上是包含了黏土附加導(dǎo)電影響的，影響程度取決于用于做實(shí)驗(yàn)的巖心的平均黏土質(zhì)量分?jǐn)?shù)水平，在實(shí)際解釋時(shí)，對(duì)儲(chǔ)層中不同深度黏土質(zhì)量分?jǐn)?shù)變化的響應(yīng)不靈敏。

與阿爾奇公式不同，西門度公式將含黏土部分和純砂巖部分的電阻率分開計(jì)算，從式（3）和式（8）中可以看出，巖電參數(shù)m、n值僅僅是用于純砂巖部分電阻率計(jì)算的，僅與純砂巖的物性和潤濕性等性質(zhì)有關(guān)而與黏土無關(guān)，故采用西門度公式計(jì)算含水飽和度時(shí)，所選用的巖電參數(shù)必須為消除了黏土附加導(dǎo)電影響的巖電參數(shù)。若簡(jiǎn)單地代入由阿爾奇模型擬合得到的巖電參數(shù)，則會(huì)導(dǎo)致重復(fù)校正，使計(jì)算得到的Sw偏離真實(shí)值。

2.2 m*和n*求取的基本原理

定義m*和n*分別為消除了黏土附加導(dǎo)電影響的膠結(jié)指數(shù)和飽和度指數(shù)，由式（9）可得：

消除了黏土附加導(dǎo)電影響的膠結(jié)指數(shù)和飽和度指數(shù)，需要運(yùn)用巖電實(shí)驗(yàn)和粒度分析原始數(shù)據(jù)擬合得到。對(duì)式（15）進(jìn)行變形，得擬合公式為：

其中除巖電實(shí)驗(yàn)可直接測(cè)量的參數(shù)外，Rcl可根據(jù)測(cè)井曲線讀值或黏土礦物分析資料得到，Rw可通過實(shí)驗(yàn)水分析得到，Vcl可通過粒度分析資料得到。令a*= 1，將等號(hào)兩邊已知量分別用向量表示，則上式可表示為矩陣向量形式：

式中：A、B、C分別為1列n行向量（n為用于擬合的實(shí)驗(yàn)樣品數(shù)）。

對(duì)于這類向量表達(dá)式，可直接使用多元一次回歸的方法，擬合得到m*、n*分別為1.855、2.017。消除黏土附加導(dǎo)電影響前后的巖電參數(shù)值對(duì)比如表1。從表中可以看出，校正前后巖電參數(shù)差異較為明顯。

表1 消除黏土附加導(dǎo)電影響前后巖電參數(shù)對(duì)比統(tǒng)計(jì)表

3 改進(jìn)模型的應(yīng)用

東海陸架盆地X氣田為一低孔低滲氣田，儲(chǔ)層主要為三角洲前緣水下分流河道砂體，巖性較細(xì)。對(duì)X氣田主要目的層段129個(gè)樣品的黏土礦物X衍射結(jié)果進(jìn)行統(tǒng)計(jì)如表2（表中數(shù)字為某類黏土礦物占黏土礦物總體積分?jǐn)?shù)）。從表中可以看出，陽離子交換能力較強(qiáng)的黏土礦物（伊利石和伊蒙混層）占黏土礦物總量的平均比重較高（51.13%），說明巖石中黏土礦物的附加導(dǎo)電能力不能忽略，故使用具有附加導(dǎo)電性校正功能的飽和度公式進(jìn)行黏土質(zhì)量分?jǐn)?shù)校正很有必要。

表2 黏土類型及所占比重統(tǒng)計(jì)表 %

由于X氣田有關(guān)黏土陽離子交換能力的分析化驗(yàn)資料很少，使得基于陽離子交換的飽和度模型（如W-S模型）的關(guān)鍵參數(shù)不易確定，因此，使用了通過測(cè)井解釋求得的黏土質(zhì)量分?jǐn)?shù)進(jìn)行附加導(dǎo)電性校正的飽和度公式（即西門度公式）。

對(duì)X氣田344個(gè)粒度分析樣品進(jìn)行統(tǒng)計(jì)，發(fā)現(xiàn)主要砂巖儲(chǔ)層段的巖心黏土質(zhì)量分?jǐn)?shù)主要為0~3%（如圖2a），泥質(zhì)質(zhì)量分?jǐn)?shù)為0~12%（如圖2b），Vsh和Vcl差異明顯，故使用了改進(jìn)的西門度公式進(jìn)行逐點(diǎn)處理，突出了“附加導(dǎo)電性是由泥質(zhì)中的黏土礦物造成”這一特征，避免可能出現(xiàn)的因代入Vsh計(jì)算而得到的過于樂觀的結(jié)論。

圖2 黏土質(zhì)量分?jǐn)?shù)及泥質(zhì)質(zhì)量分?jǐn)?shù)所占比重范圍分布圖

分別用常規(guī)西門度公式和改進(jìn)的西門度模型（式15）對(duì)X氣田X-5井目的層段Y3b層進(jìn)行飽和度測(cè)井解釋，并將解釋結(jié)果用密閉取心飽和度分析結(jié)果進(jìn)行標(biāo)定（圖3）。從圖中可以看出，運(yùn)用改進(jìn)后的西門度模型計(jì)算得到的含水飽和度，與密閉取心實(shí)驗(yàn)結(jié)果符合度更高，而常規(guī)西門度公式計(jì)算得到的含水飽和度明顯偏小。經(jīng)對(duì)比，常規(guī)西門度公式計(jì)算結(jié)果的平均絕對(duì)誤差為-6.26%（測(cè)井解釋Sw減巖心分析Sw），改進(jìn)后的西門度公式計(jì)算結(jié)果的平均絕對(duì)誤差為+0.78%，證實(shí)了改進(jìn)后模型的優(yōu)勢(shì)。

圖3 模型改進(jìn)前后飽和度解釋結(jié)果對(duì)比圖

4 結(jié)論

（1）混合泥質(zhì)砂巖模型的一系列假設(shè)條件，諸如泥質(zhì)部分的電阻率增大系數(shù)等于1、泥質(zhì)部分含水飽和度和純砂巖部分相同等，可能與實(shí)際儲(chǔ)層并不相符，對(duì)測(cè)井解釋的精度造成了很大的不確定性，不利于探明儲(chǔ)量的精確確定；改進(jìn)的西門度公式將導(dǎo)電機(jī)制不同的組分拆開來，突出了黏土附加導(dǎo)電機(jī)制的不同，在缺乏陽離子交換能力資料的地區(qū)可當(dāng)做一種陽離子交換模型來使用。

（2）細(xì)粉砂和黏土對(duì)地層放射性的貢獻(xiàn)程度不同，將通過刻度模型求出的黏土質(zhì)量分?jǐn)?shù)Vcl代入飽和度模型中計(jì)算，可以避免代入泥質(zhì)質(zhì)量分?jǐn)?shù)Vsh計(jì)算可能造成的飽和度解釋偏樂觀的結(jié)果。

（3）巖電參數(shù)的擬合和使用需要運(yùn)用相同的模型，否則會(huì)使最后求出的含水飽和度偏離真實(shí)值，經(jīng)比較，消除泥質(zhì)附加導(dǎo)電影響的m、n值會(huì)略大于常規(guī)巖電實(shí)驗(yàn)所得值。

[1]ARCHIE G E. The electrical resistivity log as an aid in determining some reservoir characteristics[J]. Transactions of the AIME， 1942， 146（1）： 54-62.

[2]ARCHIE G E. Electrical resistivity an aid in core-analysis interpretation[J]. AAPG Bulletin， 1947， 31（2）： 350-366.

[3]毛志強(qiáng). 油藏條件下飽和度測(cè)井解釋基礎(chǔ)實(shí)驗(yàn)、理論和應(yīng)用研究[D]. 北京： 中國石油勘探開發(fā)科學(xué)研究院， 1995.

[4]HOSSIN A. Calcul des Saturations en Eau par la Méthode du Ciments Argileux（ formule d’Archie Généralisée）[J]. Bull Assoc Francaise Tech Pet， 1960， 140： 12-14.

[5]DE WHITE L. Relation between Resistivities and fluid Contents of Porous Rocks[J]. Oil and Gas Journal， 1950， 49（16）： 120-132.

[6]POUPON A， LOY M E， TIXIER M P. A Contribution to Electrical Log Interpretation in Shaly Sands[J]. Journal of Petroleum Technology， 1954， 6（6）： 27-34.

[7]雍世和， 張超謨. 測(cè)井?dāng)?shù)據(jù)處理與綜合解釋[M]. 東營： 中國石油大學(xué)出版社， 2007： 174-176.

[8]WAXMAN M H， SMITS L J M. Electrical Conductivities in Oil-Bearing Shaly Sands[J]. Society of Petroleum Engineers Journal，1968， 8（2）： 107-122.

[9]CLAVIER C， COATES G， DUMANOIR J. Theoretical and Experimental Bases for the Dual-Water Model for Interpretation of Shaly Sands[J]. Society of Petroleum Engineers Journal， 1984，24（2）： 153-168.

[10]王建功， 吳剛， 陳燕， 等. 一種復(fù)雜砂巖儲(chǔ)層泥質(zhì)含量評(píng)價(jià)方法[J]. 石油鉆采工藝， 2009， 31（S2）： 48-50.

[11]張德梅， 王桂萍， 婁憲剛， 等. 測(cè)井曲線組合法求取泥質(zhì)含量探討[J]. 測(cè)井技術(shù)， 2011， 35（4）： 358-362.

[12]沈建國. 用密度聲波時(shí)差和自然伽瑪測(cè)井曲線計(jì)算分散泥質(zhì)結(jié)構(gòu)泥質(zhì)含量[J]. 地球物理測(cè)井， 1989， 13（2）： 13-17.

[13]汪愛云， 劉江， 張軍龍. 儲(chǔ)集層泥質(zhì)含量普適性算法探討[J].測(cè)井技術(shù)， 2000， 24（S1）： 502-504.

[14]王功軍， 王冬梅， 陳小波， 等. ECS測(cè)井及其在西江油田的應(yīng)用[J]. 石油天然氣學(xué)報(bào)（江漢石油學(xué)院學(xué)報(bào)）， 2011， 33（7）：101-103.

[15]袁祖貴， 成曉寧， 孫娟. 地層元素測(cè)井（ECS）——一種全面評(píng)價(jià)儲(chǔ)層的測(cè)井新技術(shù)[J]. 原子能科學(xué)技術(shù)， 2004， 38（S1）：208-213.

[16]孫建孟， 李召成. 應(yīng)用自然伽馬能譜測(cè)井確定黏土礦物類型和含量[J]. 石油大學(xué)學(xué)報(bào)（自然科學(xué)版）， 1999， 23（4）： 29-32.

[17]郭影文. 自然伽馬能譜測(cè)井在黏土礦物含量分析中的應(yīng)用[J].石油天然氣學(xué)報(bào)（江漢石油學(xué)院學(xué)報(bào)）， 2008， 30（6）： 268-270.

[18]莫修文， 丁永浩， 王宏建. 葡西地區(qū)油氣儲(chǔ)層泥質(zhì)與黏土含量的計(jì)算方法研究[J]. 吉林大學(xué)學(xué)報(bào)（地球科學(xué)版）， 2006， 36（S）： 139-142.

Application of Modified Simandoux Model in Saturation Logging Interpretation

WANG Di1, WEI Feng1, CHEN Xian1, CHEN Min1, YANG Dong2
（1. Institute of Exploration and Development, SINOPEC Shanghai Offshore Oil & Gas Company, Shanghai 200120, China; 2. Zhanjiang Branch of CNOOC Ltd., Zhanjiang Guangdong 524057, China）

Simandoux Equation and Indonesia Equation, with the function of mud correction, have been widely used in saturation logging interpretation of argillaceous sandstone. Study shows that, as to the Simandoux Equation, the accuracy of saturation calculation is limited due to the limitation of assumed conditions and the parameters selected in the model， although the influence of mud additional conductivity on the formation resistivity can be eliminated by mud correction. On one hand, the conductive ability isn’t distinguished between clay and fine silt sand by Simandoux Equation； On the other hand， the rock-electrical parameters are generally empirical values or obtained by rock-electrical experiment instead of the values corrected after eliminating the influence of mud additional conductivity. In order to solve the shortages mentioned above and emphasize the role of clay in the mud additional conductivity, an improved Simandoux Equation is proposed and, on the basis of this, a method for estimating the clay content by GR curve is summarized. Thus， rock-electrical parameters without the influence of mud additional conductivity are obtained， based on the rockelectrical experiment data and multiple regression model fitting. The verification conducted in X gas field of East China Sea Shelf Basin shows that the saturation calculated by the modified Simandoux Equation is more precise than that of by Simandoux Equation. Keywords： Simandoux Equation; clay content; rock-electrical parameter; volume physical model

P631.8

A DOI:10.3969/j.issn.1008-2336.2016.03.070

1008-2336（2016）03-0070-06

2016-01-28；改回日期：2016-02-19

王迪，男，1989年生，碩士，2015年畢業(yè)于西南石油大學(xué)地球探測(cè)與信息技術(shù)專業(yè)，主要從事測(cè)井資料處理與綜合解釋等工作。E-mail: wangdi.shhy@sinopec.com。