李雪英，王鳳來(lái)，文慧儉，孟青云，宋長(zhǎng)和，馬世忠

（1.東北石油大學(xué) 地球科學(xué)學(xué)院，黑龍江 大慶 163318； 2.大慶油田有限責(zé)任公司 勘探事業(yè)部，黑龍江 大慶

163453； 3.大慶鉆探工程公司測(cè)井公司 吉林事業(yè)部，吉林 松原 138000； 4.吉林油田分公司 新木采油廠(chǎng)，吉林 松原 138000）

0 引言

我國(guó)陸相砂泥巖薄互層沉積的水淹層解釋具有工作量大、專(zhuān)業(yè)化程度高、經(jīng)驗(yàn)性強(qiáng)、不易推廣的特點(diǎn).目前水淹層解釋只是在孤立小層中進(jìn)行，忽略各小層水淹程度在沉積背景中的內(nèi)在規(guī)律性，導(dǎo)致各種水淹層解釋軟件的自動(dòng)化解釋程度不高，符合率偏低.因此，水淹層精細(xì)解釋亟需一套正確的水淹理論和科學(xué)有效的方法加以指導(dǎo)，進(jìn)而建立起一套厚層砂體水淹級(jí)別自動(dòng)判別方法.這對(duì)于減輕現(xiàn)場(chǎng)水淹層解釋強(qiáng)度、提高解釋效率和解釋精度具有重要意義.

我國(guó)主力油田長(zhǎng)期注水開(kāi)發(fā)歷史和特點(diǎn)決定水淹層解釋技術(shù)一直走在世界的前列.大慶油田在1997年提出薄差層泥質(zhì)砂巖的解釋方法和高含水期測(cè)井曲線(xiàn)高分辨率處理方法[1－2]；李全厚等提出基本測(cè)井解釋單元的概念，認(rèn)識(shí)沉積韻律與測(cè)井解釋單元之間的一一對(duì)應(yīng)關(guān)系[3].晁吉俊、鐘蘊(yùn)紫等系統(tǒng)總結(jié)水淹層在各種測(cè)井曲線(xiàn)的測(cè)井響應(yīng)特征[4－5]；劉傳平等提出巖石物理相結(jié)合動(dòng)態(tài)電阻率下降法的水淹層解釋方法[6]；楊景強(qiáng)采用相對(duì)重心、橢圓度、飽滿(mǎn)參數(shù)識(shí)別與劃分不同沉積韻律砂體，采用電阻率曲線(xiàn)和孔隙度曲線(xiàn)重疊準(zhǔn)確判別厚層砂體的高水淹層[7].張?jiān)朴⒌冉o出正、反、復(fù)合、多段多韻律、薄層、薄互層水淹模式，建立薄差層的水淹模式庫(kù)[8].我國(guó)建立厚油層細(xì)分水淹層解釋技術(shù)及薄、差油層水淹層綜合解釋技術(shù)[9－19]，代表目前水淹層研究的國(guó)際水平[20].

在精細(xì)劃分測(cè)井解釋單元和準(zhǔn)確判別沉積韻律類(lèi)型的基礎(chǔ)上，筆者提出以韻律中部為對(duì)比標(biāo)準(zhǔn)層，利用韻律內(nèi)部曲線(xiàn)電性特征相對(duì)可比性，結(jié)合各小層在韻律內(nèi)部的相對(duì)位置，給出各小層水淹級(jí)別的精確解釋方法；引入小波分解系數(shù)表征測(cè)井曲線(xiàn)的光滑程度，增加不同水淹級(jí)別之間可區(qū)分度，為水淹層精細(xì)劃分奠定基礎(chǔ).

1 水淹模式

均質(zhì)韻律厚層砂體巖性、物性較均勻，水淹后表現(xiàn)為均勻或略偏下的水線(xiàn)推進(jìn)型.在均質(zhì)韻律中，低、未水淹比較少見(jiàn)，多以中、高水淹為主，表現(xiàn)為底部為高水淹，中部以中水淹為主.水淹厚度較大，水驅(qū)效率較高（見(jiàn)圖1（a），粉色代表低水淹；淺藍(lán)色代表中水淹；深藍(lán)色代表高水淹）.水淹前后的主要電性特征見(jiàn)表1和圖2.

圖1 不同沉積韻律水淹模式Fig.1 The water flooded mode of different sedimentary rhythm

表1 均質(zhì)韻律地層水淹模式及電性特征Table 1 The water flooded mode and electrical property of homogeneous rhythm

正韻律厚層砂體由下到上巖石粒度逐漸變細(xì)，物性逐漸變差，水淹后表現(xiàn)為偏下的水線(xiàn)推進(jìn)型.在韻律中部，巖性、物性相對(duì)中等，多為中水淹；在韻律頂部多為未、低水淹；底部為高水淹，甚至達(dá)到特高水淹（見(jiàn)圖1（b））.水淹前后的主要電性特征見(jiàn)表2和圖3.

表2 正韻律地層水淹模式及電性特征Table 2 The water flooded mode and electrical property of forward rhythm

圖2 均質(zhì)韻律水淹層典型測(cè)井響應(yīng)特征Fig.2 The typical electrical property of homogeneous water－flooded zone

圖3 正韻律水淹層典型測(cè)井響應(yīng)特征Fig.3 The typical electrical property of forward rhythm water－flooded zone

反韻律厚層砂體巖性、物性由下到上逐漸變好，由于受到水驅(qū)和重力雙重作用，表現(xiàn)為略偏上的水線(xiàn)推進(jìn)型.在韻律中部，多以中水淹為主；底部為低水淹；在韻律頂部高滲透部位易形成高水淹 （見(jiàn)圖1（c））.水淹前后的主要電性特征見(jiàn)表3和圖4.

表3 反韻律地層水淹模式及電性特征Table 3 The water flooded mode and electrical property of inverse rhythm

復(fù)合韻律砂體是正、反2種韻律類(lèi)型的組合，其內(nèi)部的巖性、物性、電性特征和水淹規(guī)律分別受上下部韻律類(lèi)型控制，表現(xiàn)為中部突進(jìn)或較均勻的水線(xiàn)推進(jìn)類(lèi)型（見(jiàn)圖1（d）），韻律中部多為中、高水淹；在韻律頂部和底部為低、未水淹（見(jiàn)圖5）.

圖4 反韻律厚層水淹層模式測(cè)井響應(yīng)特征Fig.4 The typical electrical property of inverse rhythm water－flooded zone

圖5 復(fù)合韻律厚層水淹層模式測(cè)井響應(yīng)特征Fig.5 The typical electrical property of composite rhythm water－flooded zone

多段多韻律多形成于多期河道砂體中，為不同沉積時(shí)期河道的相互切蝕、疊加而成.每期河道之間以巖性、物性?shī)A層相分隔，因此每一個(gè)單韻律段都有自己獨(dú)特的水淹規(guī)律和水線(xiàn)推進(jìn)模式（見(jiàn)圖1（e））.多段多韻律地層在高分辨率聲波、密度、自然伽馬、自然電位曲線(xiàn)上表現(xiàn)為多個(gè)厚度不等的多峰組合，其間被泥質(zhì)、鈣質(zhì)夾層分隔.由于各單期河道砂體內(nèi)部水淹程度不一致，電阻率曲線(xiàn)形態(tài)呈現(xiàn)不同的電性特征，在每一個(gè)單期河道正韻律的高水淹部位電阻率曲線(xiàn)形態(tài)變化較大（見(jiàn)圖6）.

圖6 多段多韻律厚層水淹層模式測(cè)井響應(yīng)特征Fig.6 The typical electrical property of ploy－laminated rhythm water－flooded zone

2 解釋單元?jiǎng)澐峙c原始沉積韻律識(shí)別

在本質(zhì)上，測(cè)井曲線(xiàn)的解釋單元?jiǎng)澐挚赊D(zhuǎn)換為在測(cè)井曲線(xiàn)上尋找穩(wěn)定泥巖隔層.穩(wěn)定泥巖特征明顯，可以根據(jù)厚度和電阻率準(zhǔn)確識(shí)別[21].解釋單元在第一次劃分后，解釋單元分布區(qū)間過(guò)大，其中包含多個(gè)基本沉積單元，不能滿(mǎn)足后續(xù)水淹層精細(xì)劃分的要求.需要在第一次劃分的基礎(chǔ)上，提高電阻率閾值再進(jìn)行下一級(jí)別劃分，直至每一個(gè)測(cè)井解釋單元內(nèi)部只包含一個(gè)基本沉積單元.

在解釋單元?jiǎng)澐值幕A(chǔ)上，采用電阻率曲線(xiàn)聯(lián)合密度測(cè)井曲線(xiàn)，以電阻率曲線(xiàn)平均斜率所對(duì)應(yīng)的傾斜角及密度曲線(xiàn)的相對(duì)重心作為判別參數(shù)，對(duì)解釋單元內(nèi)部的沉積韻律進(jìn)行自動(dòng)劃分與判別[22－23].基本判別準(zhǔn)則：當(dāng)電阻率判別參數(shù)符合均質(zhì)韻律和正韻律特征，即認(rèn)為韻律判別類(lèi)型正確.當(dāng)電阻率判別參數(shù)給出的韻律類(lèi)型是反韻律時(shí)，有兩種可能：一是原始沉積韻律為反韻律沉積，密度曲線(xiàn)的相對(duì)重心一定大于0.5；二是原始沉積韻律為均質(zhì)韻律或者正韻律沉積，由于砂體底部物性較好的地方遭受強(qiáng)水淹，導(dǎo)致電阻率曲線(xiàn)幅值極度降低，從而呈現(xiàn)反韻律沉積，此時(shí)需要觀(guān)察密度曲線(xiàn)的變化，如果此時(shí)密度曲線(xiàn)相對(duì)重心小于0.5，則解釋單元內(nèi)部砂體原始沉積韻律為正韻律沉積；如果近似接近于0.5，則解釋單元內(nèi)部砂體原始沉積韻律為均質(zhì)韻律沉積.

3 對(duì)比標(biāo)準(zhǔn)層

無(wú)論是正韻律、均質(zhì)韻律還是反韻律，中水淹大多位于韻律的中部，與其上部的低水淹和下部的高水淹（都具有相對(duì)低幅值和低幅度差）具有顯著的可比性.選取大慶地區(qū)9口密閉取心井，以巖心分析的水淹級(jí)別為基礎(chǔ)，將層厚大于1m的166個(gè)解釋單元中部?jī)?chǔ)層的水淹級(jí)別分別進(jìn)行統(tǒng)計(jì)：中水淹為116層；高水淹為23層；低水淹為20層；未水淹為7層.解釋單元中部為中水淹的概率達(dá)到70%.中水淹表現(xiàn)較大概率與大慶地區(qū)注水開(kāi)發(fā)歷史密切相關(guān)：首先，經(jīng)過(guò)幾十年的注水開(kāi)發(fā)，大慶地區(qū)的儲(chǔ)層已經(jīng)進(jìn)入開(kāi)發(fā)中后期，厚層已經(jīng)遭受不同程度的水淹，在厚層砂體中，少見(jiàn)未水淹儲(chǔ)層，并且無(wú)論是正韻律、均質(zhì)韻律還是反韻律，其中部巖性、物性與韻律頂、底部位相比都處于中等，介于兩者之間，因此在儲(chǔ)層內(nèi)部出現(xiàn)中水淹的概率較大.其次，中水淹在整個(gè)韻律中的電性特征為相對(duì)高值.這是因?yàn)楦?、特高水淹部位電阻率降低?yán)重，未、低水淹處于韻律巖性較差部位，電性特征受巖性影響較大，呈現(xiàn)低值；由于中水淹巖性、物性介于兩者之間，沒(méi)有遭受到強(qiáng)烈的水淹改造，因此電性特征為相對(duì)高值.

在水淹級(jí)別判別過(guò)程中，選取解釋單元中部作為水淹級(jí)別判定標(biāo)準(zhǔn)，將該小層的水淹級(jí)別定為中水淹，以其特征參數(shù)作為對(duì)比標(biāo)準(zhǔn)，根據(jù)電性、物性、小波分析時(shí)頻特征與標(biāo)準(zhǔn)層之間的相對(duì)差別，確定其他小層的水淹級(jí)別.

4 水淹級(jí)別精細(xì)劃分

首先計(jì)算韻律中部中水淹對(duì)比標(biāo)準(zhǔn)層的特征參數(shù)，包括淺側(cè)向平均電阻率（）、深淺側(cè)向平均幅度差（Δ）、微電位平均電阻率（）、微電極平均電阻率幅度差（）、平均密度）、平均小波系數(shù)［24］，反映測(cè)井曲線(xiàn)的光滑程度，小波系數(shù)越小，曲線(xiàn)越光滑；電測(cè)曲線(xiàn)的小波分解原理另文論述）．

然后計(jì)算韻律內(nèi)部其他各小層的特征對(duì)比參數(shù)，包括淺側(cè)向平均電阻率（）、深淺側(cè)向平均幅度差（）、微電位平均電阻率（、微電極平均電阻率幅度差（?平均密度（?、平均小波系數(shù)（）；

再考慮待判別的各小層與標(biāo)準(zhǔn)層的相對(duì)位置、特征參數(shù)與對(duì)比標(biāo)準(zhǔn)層的大小關(guān)系，以及曲線(xiàn)的光滑程度（由小波系數(shù)反映），從而判斷正韻律、均質(zhì)韻律和反韻律中各小層的水淹級(jí)別（見(jiàn)表4－6）．

對(duì)于復(fù)合韻律，首先根據(jù)平均斜率和相對(duì)重心的變化將復(fù)合韻律分成正韻律和反韻律；然后在每個(gè)單一韻律內(nèi)部分別選擇中部作為對(duì)比標(biāo)準(zhǔn)層，按照正韻律和反韻律的小層水淹級(jí)別評(píng)判方法，給出復(fù)合韻律中各小層的水淹級(jí)別.

在多段多韻律（多期河道切疊砂體）內(nèi)部，尋找?guī)r性、物性?shī)A層，判斷夾層的穩(wěn)定性.如果單韻律（單期河道）之間的夾層厚度大于0.5m，并且微電極曲線(xiàn)的回返小于本地區(qū)的微電極曲線(xiàn)回返的截止值，則可判定該夾層穩(wěn)定，說(shuō)明夾層上下為2個(gè)不同水動(dòng)力單元，可以按照各自韻律類(lèi)型判斷水淹級(jí)別；如果夾層不滿(mǎn)足該標(biāo)準(zhǔn)，則將夾層上下的砂體視為同一個(gè)韻律類(lèi)型判斷水淹級(jí)別.

表4 正韻律水淹級(jí)別判別結(jié)果Table 4 The discrimination parameter of water flooded levels in forward rhythm

表5 均質(zhì)韻律水淹級(jí)別判別結(jié)果Table 5 The discrimination parameter of water flooded levels in homogeneous rhythm

表6 反韻律水淹級(jí)別判別結(jié)果Table 6 The discrimination parameter of water flooded levels in inverse rhythm

5 結(jié)果評(píng)價(jià)

根據(jù)水淹級(jí)別判斷原則與劃分方法，采用C語(yǔ)言編寫(xiě)處理程序，利用大慶地區(qū)11口密閉取心井對(duì)該程序進(jìn)行檢驗(yàn)與測(cè)試（見(jiàn)圖2－6），以巖心分析解釋結(jié)論為標(biāo)準(zhǔn)，將程序自動(dòng)判別的水淹級(jí)別結(jié)果與巖心分析結(jié)果進(jìn)行對(duì)比與校驗(yàn)；對(duì)于存在明顯錯(cuò)誤解釋結(jié)論的層位，采用人機(jī)交互方式對(duì)解釋單元的頂、底界和韻律類(lèi)型加以修正.11口密閉取心井的處理結(jié)果表明，采用該套水淹級(jí)別自動(dòng)判別分析方法，單井嚴(yán)格符合率（與巖心分析結(jié)果一一對(duì)應(yīng)）在75%以上，總符合率達(dá)到78%以上（見(jiàn)表7），為后續(xù)的水淹級(jí)別定量解釋分析奠定基礎(chǔ).

表7 大慶地區(qū)9口密閉取心井水淹解釋結(jié)果Table 7 The water flooded interpretation result of sealed core well in research area

6 結(jié)論

（1）厚層砂體內(nèi)部的水淹規(guī)律主要受儲(chǔ)層原始沉積韻律控制.在同一個(gè)沉積單元內(nèi)部，各有效儲(chǔ)層的空間水淹級(jí)別分布具有規(guī)律性，并且各小層的電性特征具有相對(duì)可比性.沉積單元內(nèi)部水淹規(guī)律高度相關(guān)，外部又彼此獨(dú)立，是儲(chǔ)層水淹層精細(xì)解釋的本質(zhì)機(jī)理.

（2）穩(wěn)定的泥巖夾層對(duì)多期河道砂體內(nèi)部的水淹規(guī)律具有控制與分割作用，夾層的識(shí)別與性質(zhì)的判斷是復(fù)合韻律及多段多韻律的水淹級(jí)別判別的前提和重要基礎(chǔ).

（3）以韻律中部為對(duì)比的標(biāo)準(zhǔn)層，韻律內(nèi)部其他各小層根據(jù)空間位置、測(cè)井曲線(xiàn)的相對(duì)平均值、小波分解系數(shù)精細(xì)判斷相應(yīng)的水淹級(jí)別.現(xiàn)場(chǎng)資料處理結(jié)果表明，該套水淹級(jí)別自動(dòng)判別分析方法可以精確給出沉積韻律內(nèi)部各小層的水淹級(jí)別，總符合率達(dá)到78%以上.

致謝：部分研究?jī)?nèi)容受到東北石油大學(xué)“油氣藏形成機(jī)理與資源評(píng)價(jià)”黑龍江省重點(diǎn)實(shí)驗(yàn)室資助！

[1]荊萬(wàn)學(xué).薄差層泥質(zhì)砂巖水淹層測(cè)井資料綜合解釋方法[J].國(guó)外測(cè)井技術(shù)，1997，12（2）：67－71.Jing Wanxue.Well logging Comprehensive interpretation method in thin and poor shaly sandstone watered－out pay zone[J].Foreign Well Logging Technology，1997，12（2）：67－71.

[2]常明澈，邵豐，謝榮華.非均質(zhì)水淹層測(cè)井與解釋[J].國(guó)外測(cè)井技術(shù)，1997，12（3）：43－55.Chang Mingche，Shao Feng，Xie Ronghua.Well logging and interpretation in anisotropic water－flooded zone[J].Foreign Well Logging Technology，1997，12（3）：43－55.

[3]李全厚，張?jiān)朴?，姜萍，?狀態(tài)空間模型在薄差層水淹層測(cè)井解釋中的應(yīng)用[J].測(cè)井技術(shù)，2001，25（4）：291－293.Li Quanhou，Zhang Yunying，Jiang Ping，et al.Application of state space model to interpretation of thin water flooded zone[J].Well Logging Technology，2001，25（4）：291－293.

[4]晁吉俊.利用測(cè)井資料綜合判斷老君廟油田L(fēng)油層水淹等級(jí)[J].國(guó)外測(cè)井技術(shù)，2002，17（4）：21－28.Chao Jijun.Comprehensive judgment the water－flooded level using by well logging data in L oil reservoir of Laojunmiao oilfield[J].Foreign Well Logging Technology，2002，17（4）：21－28.

[5]鐘蘊(yùn)紫，秦菲莉，朱勝文，等.高含水期水淹層評(píng)價(jià)方法[J].國(guó)外測(cè)井技術(shù)，2002，17（3）：9－12.Zhong Yunzi，Qin Feili，Zhu Shengwen，et al.Evaluation method of water－flooded layer in high water cut stage[J].Foreign Well Logging Technology，2002，17（3）：9－12.

[6]劉傳平，楊青山，楊景強(qiáng)，等.薄差層水淹層測(cè)井解釋技術(shù)研究[J].大慶石油地質(zhì)與開(kāi)發(fā)，2004，23（5）：118－120.Liu Chuanping，Yang Qingshan，Yang Jingqiang，et al.Well logging interpretation technique in thin and poor watered－out pay zone[J].Petroleum Geology ＆ Oilfield Development in Daqing，2004，23（5）：118－120.

[7]楊景強(qiáng)，樊太亮，馬宏宇，等.不同韻律油層水淹層測(cè)井定性識(shí)別方法[J].大慶石油地質(zhì)與開(kāi)發(fā)，2010，29（1）：134－138.Yang Jingqiang，F(xiàn)an Tailiang，Ma Hongyu，et al.Qualitative identifing method of watered－out oil reservoirs with various rhythms by well logging[J].Petroleum Geology ＆ Oilfield Development in Daqing，2010，29（1）：134－138.

[8]張?jiān)朴?，王鐵燕.喇薩杏油田水淹模式庫(kù)的建立及其在水淹層定性解釋中的應(yīng)用[J].測(cè)井技術(shù)，2005，29（6）：548－550.Zhang Yunying，Wang Tieyan.Establishment of watered－out mode based and its application in qualitative interpretation in watered－out reservoir[J].Well Logging Technology，2005，29（6）：548－550.

[9]藺景龍，張慶國(guó)，宋延杰，等.水淹層測(cè)井分析[J].大慶石油學(xué)院學(xué)報(bào)，2001，25（3）：20－23.Lin Jinglong，Zhang Qingguo，Song Yanjie，et al.Water－flooded formation well logging analysis[J].Journal of Daqing Petroleum Institute，2001，25（3）：20－23.

[10]楊景強(qiáng)，樊太亮，馬宏宇，等.利用儲(chǔ)層分類(lèi)進(jìn)行水淹層測(cè)井解釋的方法研究[J].測(cè)井技術(shù)，2010，34（3）：238－241.Yang Jingqiang，F(xiàn)an Tailiang，Ma Hongyu，et al.On log interpretation of water flooded layer with reservoir classification method[J].Well Logging Technology，2010，34（3）：238－241.

[11]宋延杰，王暢，唐曉敏，等.泥質(zhì)砂巖水淹層有效介質(zhì)對(duì)稱(chēng)復(fù)電阻率模型[J].大慶石油學(xué)院學(xué)報(bào)，2009，33（3）：13－17.Song Yanjie，Wang Chang，Tang Xiaomin，et al.Symmetrical effective medium complex resistivity model for water flooded shaly sand formation[J].Journal of Daqing Petroleum Institute，2009，33（3）：13－17.

[12]楊景強(qiáng)，馬宏宇，劉如紅，等.對(duì)水淹層測(cè)井評(píng)價(jià)中幾個(gè)問(wèn)題的認(rèn)識(shí)[J].測(cè)井技術(shù)，2009，23（6）：511－516.Yang Jingqiang，Ma Hongyu，Liu Ruhong，et al.Cognitions on some problems in water－flooded layer log evaluation[J].Well Log－ging Technology，2009，23（6）：511－516.

[13]侯連華，王京紅，劉澤容.水淹層測(cè)井評(píng)價(jià)方法[J].石油學(xué)報(bào)，1999，20（3）：49－55.Hou Lianhua，Wang Jinghong，Liu Zerong.Evaluation of water flooded interval well logging[J].Acta Petrolei Sinica，1999，20（3）：49－55.

[14]李雪英，藺景龍，荊萬(wàn)學(xué)，等.薄差水淹儲(chǔ)層自然電位測(cè)井曲線(xiàn)的校正方法[J].大慶石油學(xué)院學(xué)報(bào)，2007，31（3）：11－14.Li Xueying，Lin Jinglong，Jing Wanxue，et al.The environment factor correction method of spontaneous potential curve in thin and poor water flooded reservoirs[J].Journal of Daqing Petroleum Institute，2007，31（3）：11－14.

[15]王任一，劉曰強(qiáng).基于小波熵水淹層識(shí)別方法探索[J].地球物理學(xué)進(jìn)展，2007，22（5）：1558－1561.Wang Renyi，Liu Yueqiang.The identification method of the water flooded layer based on the wavelet entropy[J].Progress in Geophysics，2007，22（5）：1558－1561.

[16]張慶國(guó)，鮑志東，那未紅.注水開(kāi)發(fā)油田水淹層測(cè)井響應(yīng)特征[J].大慶石油學(xué)院學(xué)報(bào)，2006，30（4）：101－105.Zhang Qingguo，Bao Zhidong，Na Weihong.Characteristics of well－logging response to water－flooded reservoir in water－driven oil field[J].Journal of Daqing Petroleum Institute，2006，30（4）：101－105.

[17]張美玲，孫寶剛，謝磊，等.扶楊油層組含鈣質(zhì)夾層厚油儲(chǔ)層水淹層識(shí)別[J].大慶石油學(xué)院學(xué)報(bào)，2011，35（2）：13－18.Zhang Meilong，Sun Baogang，Xie Lei，et al.Water flooded layer identification of thick oil reservoir of Fuyang oil layer with calcareous intercalation[J].Journal of Daqing Petroleum Institute，2011，35（2）：13－18.

[18]Dabbouk C，Alin L，Williams G，et al.Waterflood in vuggy layer of Middle East reservoir：Displacement physics undersood[R].SPE 78530，2002.

[19]Khisamutdinov A I，Phedorin MA.Numerrical method of evaluationg elemental content of oil－water saturated formations based on pulsed neutron–gamma inelastic log data[R].SPE 104342，2006.

[20]慈建發(fā)，何世明，李振英，等.水淹層測(cè)井現(xiàn)狀與未來(lái)[J].天然氣工業(yè)，2005，25（7）：44－46.Ci Jianfa，He Shiming，Li Zhenying，et al.Present situation of water flooded zone logging and its development trend[J].Natural Gas Industry，2005，25（7）：44－46.

[21]文慧儉，李雪英，孟青云，等.水淹層解釋單元自動(dòng)劃分方法[J].地球物理學(xué)進(jìn)展，2013，28（3）：1516－1521.Wen Huijian，Li Xueying，Meng Qingyun，et al.The automatic division method of the interpretation unit for water－flooded zones[J].Progress in Geophysics，2013，28（3）：1516－1521.

[22]梅紅，張厚福，鐘興水，等.神經(jīng)網(wǎng)絡(luò)技術(shù)在測(cè)井相分析及水淹層判別中的應(yīng)用[J].石油大學(xué)學(xué)報(bào)：自然科學(xué)版，1997，21（3）：24－28.Mei Hong，Zhang Houfu，Zhong Xingshui，et al.Application of neural network to analyzing logging facies and identifying flooding layers[J].Journal of the University of Petroleum，China：Edition of Natural Science，1997，21（3）：24－28.

[23]文慧儉，李雪英，王盛強(qiáng)，等.基于解釋單元?jiǎng)澐值亩嗲€(xiàn)聯(lián)合確定地層原始沉積韻律[J].吉林大學(xué)學(xué)報(bào)：地球科學(xué)版，2013，43（1）：299－305.Wen Huijian，Li Xueying，Wang Shengqiang，et al.Multi－curve jointly determine the original sedimentaryrhythm based on the interpretation unit discrimination[J].Journal of Jilin University：Earth Science Edition，2013，43（1）：299－305.

[24]陳科貴，趙志恒，高偉，等.小波分析在水淹層定性評(píng)價(jià)中的應(yīng)用[J].天然氣工業(yè)，2008，28（8）：49－51.Chen Kegui，Zhao Zhiheng，Gao Wei，et al.Application of wavelet to qualitative evaluation on water flooded layer[J].Natural Gas Industry，2008，28（8）：49－51.