不同沉積韻律控制的厚層砂體水淹級別精細(xì)劃分

李雪英，王鳳來，文慧儉，孟青云，宋長和，馬世忠

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

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

0 引言

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

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

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

1 水淹模式

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

5 結(jié)果評價(jià)

根據(jù)水淹級別判斷原則與劃分方法，采用C語言編寫處理程序，利用大慶地區(qū)11口密閉取心井對該程序進(jìn)行檢驗(yàn)與測試（見圖2－6），以巖心分析解釋結(jié)論為標(biāo)準(zhǔn)，將程序自動判別的水淹級別結(jié)果與巖心分析結(jié)果進(jìn)行對比與校驗(yàn)；對于存在明顯錯(cuò)誤解釋結(jié)論的層位，采用人機(jī)交互方式對解釋單元的頂、底界和韻律類型加以修正.11口密閉取心井的處理結(jié)果表明，采用該套水淹級別自動判別分析方法，單井嚴(yán)格符合率（與巖心分析結(jié)果一一對應(yīng)）在75%以上，總符合率達(dá)到78%以上（見表7），為后續(xù)的水淹級別定量解釋分析奠定基礎(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ī)律主要受儲層原始沉積韻律控制.在同一個(gè)沉積單元內(nèi)部，各有效儲層的空間水淹級別分布具有規(guī)律性，并且各小層的電性特征具有相對可比性.沉積單元內(nèi)部水淹規(guī)律高度相關(guān)，外部又彼此獨(dú)立，是儲層水淹層精細(xì)解釋的本質(zhì)機(jī)理.

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

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

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

[1]荊萬學(xué).薄差層泥質(zhì)砂巖水淹層測井資料綜合解釋方法[J].國外測井技術(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ì)水淹層測井與解釋[J].國外測井技術(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)空間模型在薄差層水淹層測井解釋中的應(yīng)用[J].測井技術(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]晁吉俊.利用測井資料綜合判斷老君廟油田L(fēng)油層水淹等級[J].國外測井技術(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)紫，秦菲莉，朱勝文，等.高含水期水淹層評價(jià)方法[J].國外測井技術(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)，等.薄差層水淹層測井解釋技術(shù)研究[J].大慶石油地質(zhì)與開發(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)，樊太亮，馬宏宇，等.不同韻律油層水淹層測井定性識別方法[J].大慶石油地質(zhì)與開發(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)朴ⅲ蹊F燕.喇薩杏油田水淹模式庫的建立及其在水淹層定性解釋中的應(yīng)用[J].測井技術(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]藺景龍，張慶國，宋延杰，等.水淹層測井分析[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)，樊太亮，馬宏宇，等.利用儲層分類進(jìn)行水淹層測井解釋的方法研究[J].測井技術(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ì)對稱復(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)，馬宏宇，劉如紅，等.對水淹層測井評價(jià)中幾個(gè)問題的認(rèn)識[J].測井技術(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]侯連華，王京紅，劉澤容.水淹層測井評價(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]李雪英，藺景龍，荊萬學(xué)，等.薄差水淹儲層自然電位測井曲線的校正方法[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).基于小波熵水淹層識別方法探索[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]張慶國，鮑志東，那未紅.注水開發(fā)油田水淹層測井響應(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ì)夾層厚油儲層水淹層識別[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ā)，何世明，李振英，等.水淹層測井現(xiàn)狀與未來[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]文慧儉，李雪英，孟青云，等.水淹層解釋單元自動劃分方法[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ù)在測井相分析及水淹層判別中的應(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)澐值亩嗲€聯(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]陳科貴，趙志恒，高偉，等.小波分析在水淹層定性評價(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.