屈信忠, 唐文生, 姜明忠, 丁曉軍, 廖春, 宋祖勇, 劉青文

(1.中國石油青海油田公司勘探開發(fā)研究院, 甘肅 敦煌 736202;2.北京泰姆泰克石油科技開發(fā)有限公司, 北京 昌平 100220)

0 引 言

高鹽油藏由于地層水礦化度高,導(dǎo)電能力較強(qiáng),油砂體普遍呈現(xiàn)較低的電阻率[1]。油層水淹后,電阻率進(jìn)一步下降,油層、水層、水淹層的區(qū)分識(shí)別難度進(jìn)一步加大。多數(shù)油藏開發(fā)過程中采用污水回注方法,注入水礦化度無規(guī)律性,油藏混合液礦化度的確定難度大,地層水導(dǎo)電特征難以量化,成為水淹層飽和度計(jì)算的主要難題[2]。

油砂體水淹后,如果發(fā)生礦化度的明顯改變,在測(cè)井資料上具有一定的響應(yīng)特征。例如利用自然電位測(cè)井的幅度變化、基線偏移程度等特征定量評(píng)價(jià)混合液電阻率[3-4]。理論上,這種方法可行[5],但是實(shí)際效果并不令人滿意。對(duì)于巖性均勻的厚砂體,自然電位可在一定程度上反映儲(chǔ)集層礦化度變化[6],但多數(shù)油藏的砂體厚度和巖性變化很大,自然電位變化幅度主要是厚度、滲透性、黏土組分等的體現(xiàn),資料本身還需要大量的影響因素校正[7],視自然電位很難體現(xiàn)礦化度變化。由于混合液礦化度難以量化,在剩余油飽和度的解釋中,電測(cè)井資料應(yīng)用受到嚴(yán)重制約,出現(xiàn)了諸多借助于數(shù)學(xué)統(tǒng)計(jì)等手段進(jìn)行飽和度計(jì)算的方法探討[8-12],或借助于非電法資料認(rèn)識(shí)飽和度的方法研究[13-14]。這也是經(jīng)典的Archie基礎(chǔ)模型自提出與發(fā)展70余年來,在油藏開發(fā)解釋應(yīng)用中遇到的主要技術(shù)障礙之一[15-16]。

本文結(jié)合柴達(dá)木盆地多個(gè)大中型油藏水淹層的研究,首次提出用于水淹層飽和度解釋的礦化度驅(qū)替-交換理論模型,在此基礎(chǔ)上建立開發(fā)期水淹層分級(jí)評(píng)價(jià)體系,成為解決開發(fā)期飽和度計(jì)算和水淹定量評(píng)價(jià)的重要技術(shù)手段。研究認(rèn)為,注水開發(fā)過程中油藏混合液礦化度的變化主要受注入水驅(qū)替和離子交換2種因素控制。水淹初期,礦化度變化主要受孔隙體積中的注入水比例影響;隨著注水時(shí)間和注水倍數(shù)增加,注入水和原生水之間的離子交換程度不斷增加。通過建立礦化度驅(qū)替-交換物理模型,確定出混合液礦化度與含水飽和度、原生水礦化度、注入水礦化度的關(guān)系,聯(lián)合電測(cè)井資料,采用迭代計(jì)算的方法,消除飽和度解釋中的雙解問題,大幅度提高飽和度和礦化度的計(jì)算精度。以該理論模型為核心,建立涵蓋水淹層飽和度計(jì)算、含水率預(yù)測(cè)、水淹層分級(jí)解釋為主的評(píng)價(jià)體系,并在多個(gè)大中型油藏投入生產(chǎn)應(yīng)用,取得了良好的應(yīng)用效果。

1 注入開發(fā)油藏礦化度的驅(qū)替-交換理論模型

1.1 注水驅(qū)替過程中礦化度變化的體積模型

油藏水驅(qū)過程中礦化度的變化基于巖石物理體積模型建立(見圖1)。假定單位體巖石積內(nèi)孔隙度為φ,原始含水飽和度為Swi,原始地層水礦化度為穩(wěn)定值(礦化度記為P0),在注入水(礦化度記為Pin)的驅(qū)替過程中,一部分烴類體積被注入水替換,含水飽和度升高到Sw時(shí),替換進(jìn)入的注入水體積為φ×(Sw-Swi),地層內(nèi)混合液礦化度(記為Pmix)是原生部分和替換部分的體積加權(quán)平均值。其中原生水體積包含的離子量為φSwiP0,被替換體積包含的離子量為φ(Sw-Swi)Pin,單位巖石體積內(nèi)混合液礦化度Pmix表述為

(1)

式中,φ為巖石孔隙度,小數(shù);Sw為含水飽和度,小數(shù);Swi為束縛水飽和度,小數(shù);P0為原始地層水礦化度,mg/L;Pin為注入水礦化度,mg/L;Pmix為水淹層混合液礦化度,mg/L。

可見,注水驅(qū)替過程中,混合液礦化度與原生水礦化度、注入水礦化度、飽和度變化情況密切相關(guān)。對(duì)于清污水回注的油藏而言,通過及時(shí)監(jiān)測(cè)注入水礦化度,可以為混合液礦化度的量化計(jì)算奠定基礎(chǔ)。

1.2 注水驅(qū)替過程中的離子交換

不同濃度和類型的水體接觸,將逐步發(fā)生離子交換作用。在水淹初期,注入水占比例較小,與原生水的接觸體積和接觸時(shí)間有限,這種交換作用相對(duì)較弱。隨著開發(fā)程度深入,注入水驅(qū)替倍數(shù)的增加或注入水在巖石中滯留的時(shí)間加長(zhǎng),注入水和巖石內(nèi)部的原生水發(fā)生離子交換作用的強(qiáng)度增加[17]。

圖1 注水驅(qū)替開發(fā)過程中水體礦化度變化的體積模型

無論是巖石內(nèi)部部分水淹還是相對(duì)均勻的水淹,注入水體積和原生地層水的交換作用應(yīng)當(dāng)是發(fā)生在2類流體接觸的條件下。油層條件下,原生水的離子存在于黏土結(jié)晶水和毛細(xì)管滯留水為主的賦存狀態(tài)下。當(dāng)巖石局部被水洗后,部分烴體積被替換,這部分體積內(nèi)的束縛水雖可以視為無流動(dòng)能力,但與注入水形成接觸關(guān)系,存在發(fā)生離子交換的基本條件。

假定巖石水洗后含水飽和度為Sw,束縛水飽和度為Swi,殘余油飽和度為Sor,則巖石被水洗的體積比例(記為Vb)可以表述為

(2)

理論上,我們可以假定在沒有水洗的巖石體積中,可流動(dòng)的巖石喉道由烴類占據(jù),其中的束縛水和注入水沒有相接觸的條件,離子交換可以忽略。那么原生水體積中可能與注入水接觸的比例(即有條件發(fā)生離子交換的比例),即可由式(2)表示(Swi、Sor參數(shù)通過相關(guān)巖石物理研究模型進(jìn)行預(yù)測(cè))。

水洗初期, 原生水礦化度為P0,原生含水體積為φSwi;隨著水洗程度或注入水倍數(shù)增加,這部分礦化度逐步通過離子交換發(fā)生改變,最終趨向變化的極值是接近注入水礦化度Pin。因此,在注入驅(qū)替過程中,離子交換的程度是一個(gè)變化的動(dòng)態(tài)過程,即其中一部分原生水礦化度逐步趨向接近注入水礦化度的過程,而未水洗部分的巖石中原生水礦化度可以視為保持不變。

在這個(gè)動(dòng)態(tài)變化的過程中,由于可交換的水體比例是有限的,可以先考慮計(jì)算終極趨向條件下的礦化度變化,即認(rèn)為被水洗的部分巖石,由于注入水倍數(shù)的增加,存在接觸關(guān)系的那部分束縛水體積,其礦化度最終趨向等于注入水礦化度,從而簡(jiǎn)化其模型計(jì)算過程。

因此,在巖石終極水洗條件下,接觸注入水的那部分原生水體,礦化度趨向等于注入水礦化度,結(jié)合公式(2),這時(shí)公式(1)混合液礦化度表述變化為

(3)

針對(duì)尕斯庫勒油藏,選擇了一個(gè)原始礦化度為160 000 mg/L、束縛水飽和度30%、殘余油飽和度為28%的油層進(jìn)行分析(見圖2)。圖2中模擬計(jì)算出了注入水礦化度60 000～180 000 mg/L條件下的含水飽和度—混合液礦化度關(guān)系。在注入水礦化度低于原始水礦化度的條件下,隨含水飽和度升高,混合液礦化度成持續(xù)下降的曲線。注入水礦化度越低,下降速度越快。

統(tǒng)計(jì)表明,由于多數(shù)油田采用清污水回注,油田開發(fā)過程中采出水礦化度都有所下降,回注污水礦化度高于原始地層水礦化度的情況基本不會(huì)出現(xiàn)。以尕斯庫勒油田為例,原始地層水礦化度160 000 mg/L左右,回注水礦化度在70 000～130 000 mg/L之間。

1.3 離子交換程度的影響分析

注水驅(qū)替過程中,離子替換量、離子交換量是控制混合液礦化度的主要因素。其中離子交換過程難以判別和量化分析,圖2采用的簡(jiǎn)化算法(即認(rèn)為束縛水中接觸注入水的部分逐步趨近注入水礦化度)更符合水淹中后期的情況。但是,模擬計(jì)算表明,即使在水淹初期,由于砂體內(nèi)部水洗巖石比例有限,無論接觸部分離子交換程度如何,對(duì)混合液礦化度影響都比較小。以圖2條件為例,表1模擬計(jì)算了注入水礦化度為100 000 mg/L、注入水飽和度5%～15%時(shí)的混合液礦化度變化。假如水淹初期注入水飽和度達(dá)到10%時(shí), 在不考慮交換現(xiàn)象時(shí),計(jì)算地層混合液礦化度為145 000 mg/L;如果離子交換非常充分,模擬計(jì)算的礦化度為141 000 mg/L?？梢?種情況下差異很小,使用公式(3)計(jì)算的混合液礦化度誤差并不明顯;開發(fā)后期,由于隨著水淹程度加強(qiáng),離子交換的程度必然會(huì)越來越強(qiáng),更趨近于交換充分條件下的計(jì)算結(jié)果,利用公式(3)計(jì)算的混合液礦化度誤差也應(yīng)該更小?？梢?礦化度離子驅(qū)替-交換模型是一個(gè)非常符合油藏實(shí)際開發(fā)規(guī)律的巖石物理模型。另一方面,水體離子交換是一個(gè)雙向過程,實(shí)際上是束縛水中部分水體與注入水中的部分水體接觸形成的局部動(dòng)態(tài)離子平衡。由于注入水處在一個(gè)動(dòng)態(tài)的更替過程,經(jīng)模擬計(jì)算不同注入倍數(shù)下的影響,認(rèn)為離子交換對(duì)注入水水體形成的影響較小,可以忽略不計(jì)。

圖2 不同礦化度驅(qū)替過程中混合液礦化度模擬計(jì)算圖版

表1 水淹初期離子交換作用對(duì)礦化度影響分析

2 利用驅(qū)替-交換模型認(rèn)識(shí)水淹層電性特征

2.1 水驅(qū)過程巖石電阻率理論計(jì)算

通過礦化度驅(qū)替-交換模型,可以確定不同巖石條件下的礦化度—飽和度關(guān)系。進(jìn)一步將混合液礦化度轉(zhuǎn)換為某一溫度(油藏深度)下混合液電阻率,計(jì)算出特定巖石條件下的電阻率—飽和度關(guān)系,為油藏開發(fā)工作提出指導(dǎo)和支持。電阻率計(jì)算依據(jù)Archie公式

(4)

式中,Rw,z為混合液電阻率,Ω·m,由混合液礦化度Pmix換算而來;Rt為模擬計(jì)算的儲(chǔ)層電阻率,Ω·m;a、b、m、n等巖電參數(shù)通過巖心實(shí)驗(yàn)取得;f(Pmix)表示Rw,z可以借助理論圖版或公式從混合液礦化度Pmix轉(zhuǎn)換而來。

圖3 利用離子驅(qū)替-交換模型進(jìn)行電阻率變化的理論計(jì)算(某特定油層條件)

2.2 水驅(qū)過程中電阻率變化規(guī)律

圖3是針對(duì)某一典型油層的計(jì)算結(jié)果。該油層孔隙度16%,油層溫度65 ℃,原始含油飽和度70%,原始礦化度160 000 mg/L。模擬了注入水礦化度20 000～160 000 mg/L幾種條件,結(jié)合驅(qū)替-交換理論模型,計(jì)算出了油層在水淹過程中電阻率變化曲線。注入水礦化度在100 000 mg/L以上時(shí),隨著水淹程度的深入,油層電阻率基本呈現(xiàn)單調(diào)下降形態(tài),電阻率從5.9 Ω·m最終下降到1.0～2.0 Ω·m之間,說明這些注水礦化度條件下電阻率受混合液礦化度影響較小,礦化度的計(jì)算誤差很難影響到對(duì)水淹級(jí)別的判斷。

當(dāng)注入水礦化度低于80 000 mg/L時(shí),電阻率曲線隨飽和度變化呈現(xiàn)“S”形態(tài),與國內(nèi)外多年來相關(guān)實(shí)驗(yàn)室測(cè)量和研究結(jié)果吻合[18-20]。采用較低礦化度驅(qū)替開發(fā),在水淹初期電阻率下降,水淹中后期出現(xiàn)平臺(tái)或升高。尤其使用礦化度小于50 000 mg/L的注入水開發(fā)時(shí),強(qiáng)水淹層的電阻率會(huì)與油層初期持平甚至更高,造成流體性質(zhì)的電性特征混淆。例如,使用20 000 mg/L注入水長(zhǎng)期驅(qū)替,理論上上述儲(chǔ)集層電阻率由油層的5.9 Ω·m下降到3.7 Ω·m,隨后在中強(qiáng)水淹條件下最高升至9.5 Ω·m,在殘余油狀態(tài)下,電阻率仍達(dá)到6.8 Ω·m,高于原始油層電阻率。

模擬計(jì)算結(jié)果說明,高鹽油藏注入水礦化度應(yīng)盡量接近原生水礦化度,尤其是避免使用鹽度低于原始礦化度50%以下的注入水,以免導(dǎo)致不同流體性質(zhì)儲(chǔ)層的電性特征混淆。

3 水淹層解釋中的應(yīng)用

3.1 通過模型迭代計(jì)算確定水淹層飽和度

目前,在柴達(dá)木盆地多個(gè)注水開發(fā)20年以上的主力油藏中,建立了以礦化度的驅(qū)替-交換模型為紐帶的水淹層解釋評(píng)價(jià)方法,在生產(chǎn)中發(fā)揮了重要作用。同時(shí)加強(qiáng)了模擬水驅(qū)過程的巖石電性實(shí)驗(yàn)研究、利用測(cè)井資料在單砂體間進(jìn)行巖石潤(rùn)濕性差異研究[21],以提高基礎(chǔ)參數(shù)精度。

在水淹層解釋中,由于混合液礦化度和含水飽和度都是未知變量,且都是主要的導(dǎo)電控制因素,通過驅(qū)替-交換模型,可量化表征各類儲(chǔ)層的飽和度與礦化度的關(guān)系,結(jié)合電測(cè)井解釋進(jìn)行迭代計(jì)算或聯(lián)合求解,求取水淹層混合液礦化度和含水飽和度:①首先選擇注水礦化度或原始地層礦化度為Pmix初始值(或采用后面計(jì)算出的Pmix),轉(zhuǎn)換為混合液電阻率Rw,z,進(jìn)而通過Archie模型求取含水飽和度

(5)

②基于計(jì)算出的含水飽和度,代入驅(qū)替-交換模型計(jì)算地層混合液礦化度,即通過式(3)進(jìn)行混合液礦化度計(jì)算;③將計(jì)算出的混合液礦化度Pmix值與前值進(jìn)行對(duì)比,按一定步長(zhǎng)增加或減小礦化度選值(亦即調(diào)整Rw,z值),循環(huán)利用公式(5)迭代計(jì)算求取含水飽和度。

當(dāng)?shù)?jì)算過程中,礦化度變化導(dǎo)致的飽和度變化小于設(shè)定誤差時(shí),表明解釋的混合液礦化度、含水飽和度滿足驅(qū)替-交換模型的要求,同時(shí)又與電測(cè)井響應(yīng)計(jì)算的飽和度一致,迭代計(jì)算結(jié)束(或者某個(gè)飽和度值時(shí),計(jì)算的出混合液礦化度與參與計(jì)算的輸入礦化度間誤差小于設(shè)定條件)。

也可以直接以聯(lián)合方程組(6)求解,得出混合液礦化度和飽和度

(6)

3.2 解釋應(yīng)用

表2對(duì)比了該區(qū)一典型油層選用不同混合液礦化度時(shí)飽和度的計(jì)算偏差。該層孔隙度為16%、原始礦化度170 000 mg/L、地層溫度65 ℃、測(cè)井電阻率為3.5 Ω·m。計(jì)算結(jié)果顯示,若地層混合液礦化度選值70 000～170 000 mg/L之間變化,則計(jì)算出對(duì)應(yīng)的含水飽和度在37.1%～54.5%之間大幅度變化,礦化度不確定性嚴(yán)重影響了解釋結(jié)果和對(duì)水淹狀況的認(rèn)識(shí)。

礦化度的驅(qū)替-交換模型為礦化度的確定提供了理論依據(jù)和實(shí)現(xiàn)手段。以表2提供的油層為例,對(duì)應(yīng)表中已經(jīng)計(jì)算出的含水飽和度,按照模型公式(3)求取對(duì)應(yīng)的混合液礦化度,結(jié)果在136 500～91 000 mg/L之間變化(見表2)。只有在礦化度假定值110 000 mg/L左右時(shí),與模型推算值113 000 mg/L最為接近。也就是說實(shí)際地層混合液礦化度應(yīng)當(dāng)為11 000 mg/L左右,含水飽和度約44.4%,才能同時(shí)滿足驅(qū)替-交換模型和Archie解釋模型(實(shí)際計(jì)算中采用計(jì)算機(jī)迭代或方程求解,結(jié)果為混合液礦化度116 000 mg/L,含水飽和度43.7%)。

表2 某典型油層條件下不同混合液礦化度選值導(dǎo)致飽和度偏差對(duì)比

該模型為注水開發(fā)過程中飽和度量化評(píng)價(jià)奠定了基礎(chǔ)。進(jìn)一步結(jié)合相滲等巖石物理實(shí)驗(yàn)成果,可以建立以飽和度和含水率為主的水淹層解釋標(biāo)準(zhǔn)和分級(jí)方法[22-24]。目前在多個(gè)高鹽油藏水淹層解釋評(píng)價(jià)標(biāo)準(zhǔn)中,將水淹層解釋劃分為油層(含水率<10%)、1級(jí)水淹(含水率=10%～40%)、2級(jí)水淹(含水率=40%～60%)、3級(jí)水淹(含水率=60%～80%)、4級(jí)水淹(含水率>80%)、水層等6個(gè)級(jí)別,并依據(jù)研究成果分別在各油藏細(xì)化了解釋標(biāo)準(zhǔn)。解釋成果中提供了含水飽和度、可動(dòng)水飽和度、混合液礦化度、含水率等重要參數(shù)的預(yù)測(cè)結(jié)果。

圖4 YX××井水淹層處理成果圖

圖4是某井解釋成果圖示例。展示的層組原始地層礦化度140 000 mg/L左右,原始油層電阻率2.0～20.0 Ω·m,注入水礦化度85 000 mg/L左右。處理解釋顯示,V1-10小層水淹較強(qiáng)(2到4級(jí)),模型計(jì)算混合液礦化度上部100 000 mg/L左右,含水率85%左右;下部混合液礦化度110 000 mg/L,含水率52%,解釋成果直觀地展示了礦化度變化、可動(dòng)水狀況、各砂體量化的水淹級(jí)別等情況。

基于驅(qū)替-交換模型建立起的高鹽油藏水淹層解釋體系,實(shí)現(xiàn)了油藏水淹層6個(gè)級(jí)別分級(jí)解釋和砂體內(nèi)部細(xì)分解釋。該模型在生產(chǎn)實(shí)踐中得到了檢驗(yàn)和驗(yàn)證,在柴達(dá)木盆地多個(gè)大中型油氣藏開發(fā)期解釋中發(fā)揮了重要作用,已經(jīng)成為剩余油挖潛和水淹層飽和度定量評(píng)價(jià)的重要技術(shù)手段之一。

4 結(jié) 論

(1) 飽和度和礦化度2個(gè)變量的量化,是注水開發(fā)油藏剩余油飽和度解釋中面對(duì)的雙解問題。首次提出的礦化度驅(qū)替-交換巖石物理模型,創(chuàng)新和完善了開發(fā)期飽和度定量計(jì)算的理論方法,有效破解了注水開發(fā)油藏剩余油評(píng)價(jià)中長(zhǎng)期存在的主要技術(shù)難題。

(2) 以礦化度的驅(qū)替-交換模型為核心,結(jié)合電測(cè)井資料和巖石物理實(shí)驗(yàn)、雙相滲流實(shí)驗(yàn)等研究,建立起儲(chǔ)集層含水飽和度、可動(dòng)水飽和度、殘余油飽和度等重要參數(shù)的定量解釋。通過飽和度評(píng)價(jià)和含水率預(yù)測(cè),達(dá)到水淹層分級(jí)解釋,可以形成注水開發(fā)油藏量化、系統(tǒng)化的水淹層評(píng)價(jià)體系。

(3) 模型計(jì)算表明,注入水礦化度低于原始礦化度50%后,水淹層電阻率將打破隨著水淹程度加強(qiáng)而單調(diào)下降的趨勢(shì)。若注水礦化度進(jìn)一步下降,水淹層電阻率將出現(xiàn)升高趨勢(shì),導(dǎo)致油水性質(zhì)識(shí)別難度加大。生產(chǎn)中應(yīng)控制好注入水礦化度的配比,避免水淹層“S”型電性形態(tài)的出現(xiàn)。

(4) 注入水礦化度是用于驅(qū)替-交換模型計(jì)算的重要參數(shù),而采出水礦化度可作為模型可靠性監(jiān)測(cè)和驗(yàn)證的重要資料。生產(chǎn)中應(yīng)做好基礎(chǔ)資料采集工作,加強(qiáng)回注水和采出水礦化度監(jiān)測(cè)。

參考文獻(xiàn):

[1] 孫建孟, 王克文, 朱家俊. 濟(jì)陽坳陷低電阻率儲(chǔ)層電性微觀影響因素研究 [J]. 石油學(xué)報(bào), 2006, 27(5): 61-65.

[2] 劉中奇, 崔琳, 董婷, 等. 孤島油田水淹層地層水電阻率計(jì)算方法研究 [J]. 測(cè)井技術(shù), 2012, 36(1): 37-40.

[3] 李漢川, 荊萬學(xué), 張海寧, 等. 基于近井壁自然電位的大慶油田水淹儲(chǔ)層地層水電阻率計(jì)算方法 [J]. 測(cè)井技術(shù), 2017, 41(1): 41-45.

[4] 楊瑞明, 王良書, 徐正順, 等. 大慶油田高含水期自然電位的變異和校正 [J]. 石油勘探與開發(fā), 2003, 30(1): 60-61.

[5] 牛超群, 張守謙, 郭余峰, 等. 自然電流測(cè)井機(jī)理的研究 [J]. 石油學(xué)報(bào), 1987, 8(4): 56-64.

[6] 唐文生, 歐陽健, 劉興禮, 等. 天然水淹油藏的測(cè)井分析 [J]. 測(cè)井技術(shù), 1995, 19(6): 421-427.

[7] 蔡志杰, 陳娓, 李大潛, 等. 三維自然電位測(cè)井的數(shù)學(xué)模型與求解方法 [J]. 測(cè)井技術(shù), 2017, 41(1): 20-25.

[8] 高興軍, 于興河, 李勝利, 等. 利用神經(jīng)網(wǎng)絡(luò)技術(shù)預(yù)測(cè)剩余油分布 [J]. 石油學(xué)報(bào), 2005, 26(3): 60-63.

[9] 豐全, 羅娜. 大慶長(zhǎng)垣水淹層測(cè)井評(píng)價(jià)技術(shù)的研究 [J]. 石油學(xué)報(bào), 1998, 19(3): 68-72.

[10] 宋考平, 吳玉樹, 計(jì)秉玉. 水驅(qū)油藏剩余油飽和度分布預(yù)測(cè)的函數(shù)法 [J]. 石油學(xué)報(bào), 2006, 27(3): 91-95.

[11] 唐文忠, 陳彬, 翟雨陽, 等. 深度延遲人工神經(jīng)網(wǎng)絡(luò)判別水淹程度 [J]. 石油與天然氣地質(zhì), 2001, 22(3): 249-252.

[12] CHEN S L. The Application of Comprehensive Fuzzy Judgement in the Interpretation of Water-flooded Reservoirs [J]. Fuzzy Mathematics, 2001, 9(3): 739-743.

[13] ZHANG X A, WANG Z W. Distinguishing Oil and Water Layers by Interpreting Acoustic Logging Data with Changing Well Diameters [J]. Geophysical Prospecting, 2015, 63(3): 669-679.

[14] 李來林, 陳小宏, 李延峰. 利用振幅譜差識(shí)別水淹層的變化 [J]. 地球物理學(xué)報(bào), 2007, 50(6): 1899-1904.

[15] ARCHIE G E. The Electrical Resistivity Log as an Aid in Determining Some Reservoir Characteristics [J]. AIME, 1942, 146: 54-61.

[16] 孫建孟, 王克文, 李偉. 測(cè)井飽和度解釋模型發(fā)展及分析 [J]. 石油勘探與開發(fā), 2008, 35(1): 101-107.

[17] 許長(zhǎng)福, 劉紅現(xiàn), 錢根寶, 等. 克拉瑪依礫巖儲(chǔ)集層微觀水驅(qū)油機(jī)理 [J]. 石油勘探與開發(fā), 2011, 38(6): 725-732.

[18] FAN Y R, DENG H G, ZHOU C C, et al. Experimental Study and Theoretical Analyses of Formation Resistivity under Fresh Water Flooding [C]∥SPWLA 42nd Annual Logging Symposium, 2011.

[19] 唐文生, 符曉穎, 張高久, 等. 南陽油田高阻水淹層成因與解釋方法淺探 [J]. 測(cè)井技術(shù), 1998, 22(3): 53-57.

[20] DANG C, LONG N, NGUYEN N, et al. Mechanistic Modeling of Low Salinity Water Flooding [J]. Petroleum Science & Engineering, 2016, 146: 209.

[21] 李春霞, 唐文生, 張超謨, 等. 用測(cè)井資料確定儲(chǔ)層潤(rùn)濕性及對(duì)應(yīng)的水淹層評(píng)價(jià)模型 [J]. 測(cè)井技術(shù), 2010, 34(3): 284-288.

[22] SHADDEL S, TABATABAE-NEJAD S A, FATHI S J. Low-salinity Water Flooding: Evaluating The Effect of Salinity on Oil and Water Relative Permeability, Wettability, and Oil Recovery [J]. Special Topics & Reviews in Porous Media, 2014, 5(2): 133-143.

[23] XU F, MU L X, WU X H, et al. New Expression of Oil/Water Relative Permeability Ratio vs Water Saturation and Its Application in Water Flooding Curve [J]. Energy Exploration & Exploitation, 2014, 32(5): 817-830.

[24] SHENG C W, DONG W. The Research for Method to Determine Interpretation Standard of Water Flooded Layer at the High Water-cut Stage [J]. Advances in Petroleum Exploration and Development, 2015, 9(2): 103-106.