王志強(qiáng)

（中國石油大慶油田有限責(zé)任公司勘探開發(fā)研究院，黑龍江 大慶163712）

0 引言

儲層流體性質(zhì)的準(zhǔn)確識別是新井試油、老井挖潛、提高油氣開發(fā)效果的基礎(chǔ)［1-2］，但低阻油層與高阻水層的出現(xiàn)，給這項工作帶來很大困難。

低阻油層從發(fā)現(xiàn)至今已有40 余年的歷史［3］。隨著人們對低阻油層認(rèn)識的不斷加深，其定義也在不斷改變。最初，人們在生產(chǎn)實踐中發(fā)現(xiàn)，有些產(chǎn)層電阻率低至1～3 Ω·m 時仍能出油，便定義為低阻油層［4-5］；后來，這一概念的絕對電阻率界限值又被擴(kuò)展到0.5 Ω·m［6］；在越來越多的實例分析基礎(chǔ)上，人們認(rèn)識到，對于產(chǎn)油層，是不存在普遍適合的電阻率范圍的。經(jīng)過大量分析，Kulha［7］明確提出，當(dāng)含油氣儲層與泥巖電阻率比值低于1.5 時，為低阻油氣層。一些學(xué)者依此提出以同一油藏體系內(nèi)的電阻增大率定義低阻油氣層［8-9］。對于高阻水層，國內(nèi)外研究報道較少。

低阻油層成因總體上分為工程因素和地質(zhì)因素［10-26］，其機(jī)理研究已頗為成熟； 高阻水層與低阻油層成因相反，但也可劃分為這2 種因素［27-31］。低阻油層識別的常用方法有交會圖、測井新技術(shù)、電阻率校正、常規(guī)測井信息數(shù)理分析等［2，32-37］，識別思路一般是先分析成因，然后針對成因制定相應(yīng)的識別方法。對高阻水層的研究較少，其識別方法多為定性分析或測井新技術(shù)與測井信息的數(shù)理分析［28，38-39］。

DJ 油田儲層為三角洲前緣沉積，砂體厚度小，平面連續(xù)性強(qiáng)，具有孔隙類型多樣、結(jié)構(gòu)復(fù)雜、非均質(zhì)性明顯等特點。Q 段具有較大油氣潛力，但由于儲層地質(zhì)特征復(fù)雜，低阻油層與高阻水層并存，難以識別。本文開展了復(fù)雜地質(zhì)特征儲層流體識別研究，從低阻油層與高阻水層的成因切入，針對其形成機(jī)理，應(yīng)用電阻率校正與分區(qū)評價相結(jié)合的方法，建立的含鈣泥質(zhì)砂巖等效體積模型在類似地質(zhì)條件油藏中可普遍適用，具有一定的探索性。該研究成果有助于降低勘探風(fēng)險，提高開發(fā)效益。

1 地層電阻率影響因素

區(qū)內(nèi)巖心化驗分析資料證實： 低阻油層導(dǎo)電礦物體積分?jǐn)?shù)極小且分布均勻；高阻水層不含吸附油氣；鉆井液性能穩(wěn)定，無工程事故延遲測井情況?？梢?，鉆井工程、導(dǎo)電礦物與吸附油氣等因素不是研究區(qū)儲層電阻率異常的主要成因。

1.1 地層水礦化度

一般情況下，陸相沉積環(huán)境的地層水礦化度相對不穩(wěn)定［14］，即使同一油藏范圍內(nèi)，層內(nèi)礦化度也往往表現(xiàn)出較大差異，這種變化會導(dǎo)致層內(nèi)地層水導(dǎo)電能力的不同。礦化度較高的孔隙水可以形成電導(dǎo)率相對較高的導(dǎo)電網(wǎng)絡(luò)，當(dāng)這種情況存在時，泥質(zhì)體積分?jǐn)?shù)對電阻率的影響將被削弱［40］。

受沉積環(huán)境與后期斷層溝通作用影響，DJ 地區(qū)Q段地層水礦化度平面變化較大，對電阻率測井影響較為明顯（見圖1）。根據(jù)17 口井的地層水化驗資料可知，研究區(qū)內(nèi)地層水礦化度平面差異較大，但其分區(qū)性明顯，以礦化度15，20 g/L 為界，分為低值區(qū)、中值區(qū)與高值區(qū)，在不同區(qū)間內(nèi)，地層水礦化度變化相對穩(wěn)定（見圖2）。這便為流體性質(zhì)的分區(qū)評價提供了可行性。

1.2 泥質(zhì)體積分?jǐn)?shù)

許多黏土礦物（如蒙脫石、伊利石等）具有較高的陽離子交換能力［33］。當(dāng)黏土礦物與孔隙流體產(chǎn)生離子交換時，產(chǎn)生附加導(dǎo)電性，增強(qiáng)地層的導(dǎo)電能力［41-42］。DJ 地區(qū)Q 段形成于三角洲前緣環(huán)境，巖石粒度細(xì)，泥質(zhì)體積分?jǐn)?shù)高，黏土礦物的附加導(dǎo)電作用是導(dǎo)致儲層電阻率降低的一項不容忽略的因素。化驗分析資料顯示，研究區(qū)泥質(zhì)體積分?jǐn)?shù)主要介于5.0%～25.0%，最大為40.7%，最小為1.2%，說明儲層泥質(zhì)體積分?jǐn)?shù)分布不均勻。應(yīng)用120 塊樣品的泥質(zhì)體積分?jǐn)?shù)與深側(cè)向電阻率進(jìn)行擬合，擬合關(guān)系表明，泥質(zhì)體積分?jǐn)?shù)對該區(qū)地層電阻率影響明顯（見圖3）。

圖1 Q 段地層電阻率與地層水礦化度關(guān)系

圖2 Q 段地層水礦化度平面分布

圖3 Q 段地層電阻率與泥質(zhì)體積分?jǐn)?shù)關(guān)系

1.3 碳酸鹽體積分?jǐn)?shù)

碳酸鹽體積分?jǐn)?shù)對儲層電阻率具有一定影響［43-44］。碳酸鹽體積分?jǐn)?shù)的不均勻分布，導(dǎo)致電阻率偏高，無法真實反映地層流體性質(zhì)，容易將水層誤判為油、氣層［45］。巖心分析結(jié)果表明，研究區(qū)碳酸鹽體積分?jǐn)?shù)主要介于5%～25%。

為揭示大情字井青二段地層電阻率與碳酸鹽體積分?jǐn)?shù)的關(guān)系，在礦化度中值區(qū)內(nèi)，對不同流體性質(zhì)的儲層進(jìn)行深感應(yīng)電阻率與碳酸鹽體積分?jǐn)?shù)的關(guān)系擬合（見圖4），但油層由于取樣點過少，不足以擬合兩者的函數(shù)關(guān)系。

圖4 不同流體性質(zhì)儲層電阻率與碳酸鹽體積分?jǐn)?shù)關(guān)系

本次研究利用H45 等13 口試油井的16 個樣品，得到油水同層深感應(yīng)電阻率與碳酸鹽體積分?jǐn)?shù)的函數(shù)關(guān)系為

式中：Rt為地層深感應(yīng)電阻率，Ω·m；VCa為碳酸鹽體積分?jǐn)?shù)；r 為相關(guān)系數(shù)。

利用H82 等16 口試油井的23 個樣品，得到水層深感應(yīng)電阻率與碳酸鹽體積分?jǐn)?shù)的函數(shù)關(guān)系為

上述分析表明，無論地層的流體性質(zhì)如何，隨著碳酸鹽體積分?jǐn)?shù)的增加，地層電阻率均呈指數(shù)增大。如圖5所示，在不區(qū)分流體性質(zhì)的情況下，對深感應(yīng)電阻率與碳酸鹽體積分?jǐn)?shù)進(jìn)行擬合可得：

圖5 儲層電阻率與碳酸鹽體積分?jǐn)?shù)關(guān)系

2 流體識別模型

由于研究區(qū)地質(zhì)背景與實際資料情況決定了雙水模型、S-W 模型等模型的不適用性，因此，本次流體識別的主要思路為： 建立含鈣泥質(zhì)砂巖等效體積模型進(jìn)行純砂巖電阻率校正，采用孔隙度-電阻率交會法判斷儲層流體性質(zhì)。

2.1 泥質(zhì)體積分?jǐn)?shù)計算

研究區(qū)自然伽馬測井對地層泥質(zhì)體積分?jǐn)?shù)的反映較好，本次研究采用自然伽馬曲線計算泥質(zhì)體積分?jǐn)?shù)，其計算公式如下：

式中：GRSH為某一層自然伽馬測井相對值；GR 為某一井段自然伽馬測井值，API；GRmin為某一井段純砂巖自然伽馬值，API；GRmax為某一井段純泥巖地層自然伽馬值，API；Vsh為儲層泥質(zhì)體積分?jǐn)?shù)；GCUR 為與地層年代有關(guān)的經(jīng)驗系數(shù)，取值為2。

2.2 碳酸鹽體積分?jǐn)?shù)計算

碳酸鹽體積分?jǐn)?shù)的測井資料計算方法較多［46-48］。測井參數(shù)與碳酸鹽體積分?jǐn)?shù)的相關(guān)性分析表明，研究區(qū)聲波時差與微球形聚焦測井聯(lián)合參與回歸，所得碳酸鹽體積分?jǐn)?shù)與實測資料符合率最高。

因此，利用H47 等27 口井共86 個巖心分析樣品，通過聲波時差測井、微球形聚焦測井與碳酸鹽體積分?jǐn)?shù)的多元回歸，建立碳酸鹽體積分?jǐn)?shù)的測井資料求取模型：

式中：AC 為聲波時差，μs/m；Rmsf為微球形聚焦測井電阻率，Ω·m。

2.3 含鈣泥質(zhì)砂巖等效體積模型

由地層電阻率影響因素分析可知： 研究區(qū)為含鈣泥質(zhì)砂巖儲層，地層水礦化度分為高值、中值、低值3個區(qū)，各分區(qū)內(nèi)地層水礦化度相對穩(wěn)定。為消除地層水礦化度變化對地層電阻率的影響，在上述3 個分區(qū)內(nèi)分別進(jìn)行流體性質(zhì)識別。下面以中值區(qū)為例，說明本次研究所用的流體識別方法。

針對研究區(qū)儲層泥質(zhì)與碳酸鹽體積分?jǐn)?shù)高的特點，建立含鈣泥質(zhì)砂巖等效體積模型，利用此模型去除泥質(zhì)與碳酸鹽對地層真電阻率的干擾。

如圖6a所示，含鈣泥質(zhì)砂巖儲層中，設(shè)模型總體積為1，泥質(zhì)體積分?jǐn)?shù)為Vsh，碳酸鹽體積分?jǐn)?shù)為VCa，純砂巖地層體積分?jǐn)?shù)為Vs。如圖6b所示，含鈣泥質(zhì)砂巖的電阻可等效為純砂巖地層電阻、泥質(zhì)電阻與碳酸鹽電阻的并聯(lián)電阻，其電導(dǎo)率可表示為

式中：Rt，Rs，Rsh，RCa分別為含鈣泥質(zhì)砂巖電阻率、純砂巖地層電阻率、泥質(zhì)電阻率及碳酸鹽電阻率，Ω·m；Ssi，Sshj，SCak分別為第i 層純砂巖地層截面積、第j 層泥質(zhì)截面積及第k 層碳酸鹽截面積，m2。

圖6 含鈣泥質(zhì)砂巖等效體積模型

推導(dǎo)得

其中，Rsh為純泥巖段深側(cè)向電阻率，取Q 段底部標(biāo)準(zhǔn)泥巖層電阻率中值，為6.3 Ω·m；對于單位體積的純碳酸鹽而言，將VCa=1 代入式（3），可得碳酸鹽電阻率RCa=21.35 Ω·m；Rt可由深側(cè)向電阻率曲線讀取；Vsh可通過式（5）求得；VCa可由式（6）計算得來。Rs即為含鈣泥質(zhì)砂巖地層校正到純砂巖地層所對應(yīng)的電阻率。

阿爾奇公式適用于純砂巖地層，可變形為

式中：a，m，n 為常數(shù)；Rw為地層水電阻率，Ω·m；φ 為孔隙度，小數(shù)；Sw為含水飽和度，小數(shù)。

由式（9）與式（10）聯(lián)合可得

3 模型應(yīng)用

DJ 油田合理選取a=1，m=2，n=2；Rw在地層水礦化度的3 個分區(qū)內(nèi)部穩(wěn)定，可視為常數(shù)。對于地層水礦化度中值區(qū)而言，可以應(yīng)用阿爾奇公式，結(jié)合該區(qū)域內(nèi)試油水層電阻率中值求得，Rw=0.15 Ω·m。

根據(jù)式（11），在地層水礦化度相對穩(wěn)定的區(qū)域內(nèi)，在Rt，Rsh，RCa，Vsh，VCa，φ 等參數(shù)可求或可直接測量的情況下，即可求得含水飽和度Sw。

應(yīng)用該方法對地層水礦化度中值區(qū)試油層段進(jìn)行分析，結(jié)果表明，高阻水層與低阻油層均可通過含水飽和度的求取得到有效識別（見表1）。DJ 油田Q 段地層水礦化度中值區(qū)內(nèi)，油水識別效果顯著，有效解決了流體性質(zhì)識別困難的問題。

表1 Q 段地層水礦化度中值區(qū)流體識別成果

4 結(jié)論

1）DJ 油田儲層地質(zhì)特征復(fù)雜，地層水礦化度、泥質(zhì)體積分?jǐn)?shù)與碳酸鹽體積分?jǐn)?shù)的不均勻分布是形成低阻油層與高阻水層共存的主要原因。

2）基于成因機(jī)理認(rèn)識，單因素校正與分區(qū)評價相結(jié)合的流體識別方法適用DJ 油田，可以有效識別高阻水層與低阻油層。

3）含鈣泥質(zhì)砂巖等效電阻率模型是阿爾奇公式的修正，實質(zhì)是將含鈣與泥質(zhì)對電阻率影響進(jìn)行校正，等效為純砂巖地層，廣泛適用于含鈣泥質(zhì)砂巖儲層。

［1］梁麗梅，喻高明，黎明，等.神經(jīng)網(wǎng)絡(luò)模擬交會圖在H 油田低阻油層流體識別中的應(yīng)用［J］.石油地質(zhì)與工程，2011，25（1）：39-40，43.

［2］申輝林，王敏，柴婧.低孔低滲儲層油水層模糊識別方法［J］.斷塊油氣田，2007，14（3）：79-81.

［3］Granberry R J，Wilshusen R C，Koepf E H.Core analysis and electriclog data gang up on formation-evaluation problems ［J］.Oil and Gas Journal，1964，1（1）：123-128.

［4］Murphy R P，Owens W W.A new approach for low-resistivity sand log analysis［J］.Journal of Petroleum Technology，1972，24（11）：1302-1306.

［5］Tripathi S N，Domangue E J，Murdoch B T.Low-resistivity sand evaluation with the chlorine log［J］.Transactions of the SPWLA：Annual Logging Symposium，1984，25（1）：1-16.

［6］Zemanek J.Low-resistivity hydrocarbon-bearing sand reservoirs［R］.SPE15713，1987.

［7］Kulha J T，Sneider R M.Low resistivity，low contrast pays （Part Ⅱ）：Comparison of examples from Southeast Asia and the Gulf of Mexico［J］.AAPG Bulletin，1994，78（7）：21-24.

［8］歐陽健，王貴文，吳繼余，等.測井地質(zhì)分析與油氣層定量評價［M］.北京：石油工業(yè)出版社，1999：15-20.

［9］Boyd A，Darling H，Tabanou J，et al.The lowdown on low-resistivity pay［J］.Oilfield Review，1995，7（3）：4-18.

［10］Darling H L，Sneider R M.Production of low-resistivity，low-contrast reservoirs，offshore Gulf of Mexico Basin［J］.The Bulletin of the Houston Geological Society，1995，38（1）：12.

［11］黃龍，鮑志東，張文瑞，等.淡水泥漿侵入引起的低阻油氣層實例分析［J］.西南石油大學(xué)學(xué)報：自然科學(xué)版，2008，30（2）：14-18.

［12］程相志，范宜仁，周燦燦，等.基于鉆井液性能優(yōu)化設(shè)計的低對比度油層識別新技術(shù)［J］.中國石油大學(xué)學(xué)報：自然科學(xué)版，2008，32（3）：50-54.

［13］王瑞麗，孫萬華，鄒明生，等.北部灣盆地潿西南凹陷低阻油層成因分析［J］.斷塊油氣田，2010，17（5）：642-645.

［14］李厚義，趙宇芳，蔡林，等.對低電阻率油層的一個誤解及修正［J］.測井技術(shù)，2002，26（2）：148-149.

［15］陳桂菊，高興軍，張為民，等.吉林大情字井油田低阻油層的識別方法研究［J］.斷塊油氣田，2006，13（4）：86-88.

［16］齊陸寧，楊少春，時丕同，等.新北油田館上段儲層電阻率特性及影響因素［J］.斷塊油氣田，2011，18（3）：301-304.

［17］左銀卿，郝以嶺，安霞，等.高束縛水飽和度低阻油層測井解釋技術(shù)［J］.西南石油學(xué)院學(xué)報，2000，22（2）：27-31.

［18］回雪峰，吳錫令，謝慶賓，等.大港油田原始低電阻率油層地質(zhì)成因分析［J］.勘探地球物理進(jìn)展，2003，26（4）：329-332.

［19］谷麗冰，李治平.高束縛水低阻油層識別［J］.內(nèi)蒙古石油化工，2008，34（2）：115.

［20］郭方達(dá)，陳程.大情字井油田流體性質(zhì)識別方法［J］.斷塊油氣田，2008，15（3）：94-96.

［21］王友凈，宋新民，何魯平，等.高尚堡深層低阻油層的地質(zhì)成因［J］.石油學(xué)報，2010，31（3）：426-431.

［22］姚同玉，李繼山，付萍.低阻油層綜合評價參數(shù)研究［J］.斷塊油氣田，2012，19（2）：142-144.

［23］陳東亮，王杰.低阻油層的測井識別技術(shù)及成因分析［J］.內(nèi)蒙古石油化工，2007，33（7）：139-140.

［24］李厚義，黃文革，馬一平，等.低阻油層的識別與認(rèn)識［J］.中國石油勘探，2008，13（1）：68-74.

［25］張沖，毛志強(qiáng)，張超，等.低阻油層的成因機(jī)理及測井識別方法研究［J］.工程地球物理學(xué)報，2008，5（1）：48-53.

［26］Worthington P F.Recongnition and evaluation of low-resistivity pay［J］.Petroleum Geoscience，2000，6（1）：77-92.

［27］張慶國，王碩，馮小東，等.大情字井油田高阻水層成因分析與識別方法［J］.科學(xué)技術(shù)與工程，2011，11（34）：8444-8446，8451.

［28］殷樹軍，閆偉林，楊青山.凝灰質(zhì)儲層復(fù)雜油水層測井識別方法［J］.石油天然氣學(xué)報，2010，32（6）：267-270，535.

［29］嚴(yán)煥德，王天祥，趙為永，等.澀北氣田高阻水層成因及其對儲層參數(shù)的影響［J］.青海石油，2009，27（2）：55-57.

［30］劉紅，李曉華，闞朝暉，等.王集小斷塊油田高阻水層成因分析［J］.內(nèi)蒙古石油化工，2008，34（14）：139-140.

［31］張海燕，王棟梁，周紅濤.塔河油田三疊系低阻油氣層測井評價［J］.斷塊油氣田，2010，17（4）：513-516.

［32］程相志，范宜仁，周燦燦.淡水儲層中低阻油氣層識別技術(shù)［J］.地學(xué)前緣，2008，15（1）：146-153.

［33］尹旭，彭仕宓，陳建文，等.低阻油層判別技術(shù)在吉林油區(qū)的應(yīng)用［J］.大慶石油地質(zhì)與開發(fā)，2008，27（1）：126-129.

［34］連承波，趙永軍，鐘建華，等.基于支持向量機(jī)的低阻油層識別方法及應(yīng)用［J］.石油天然氣學(xué)報，2008，30（1）：80-82.

［35］Lyle W D，Williams D M.Deconvolution of well log data：An innovations approach ［J］.The Log Analyst，1987，28（3）：321-328.

［36］Austin J M，F(xiàn)aulkner T L.Magnetic resonance imaging log evaluates low-resistivity pay［J］.American Oil&Gas Reporter，1993，36（8）：52-57.

［37］田中元，卞德智，陳昊，等.改進(jìn)的PICKETT 法在Y 油田低阻油層識別中的應(yīng)用［J］.石油學(xué)報，2005，26（4）：81-84.

［38］孫宏智，李曉輝.多參數(shù)判別分析技術(shù)在解決油田高阻水層問題中的應(yīng)用［C］//中國地球物理學(xué)會.2000年中國地球物理學(xué)會年刊：中國地球物理學(xué)會第十六屆年會論文集.武漢，2000：93.

［39］李曉輝，張美玲，譚輝紅.高電阻率水層定性識別方法研究［J］.測井技術(shù)，2000，24（增刊1）：491-493，555.

［40］曾文沖.低電阻率油氣層的類型、成因及評價方法（上）［J］.地球物理測井，1991，15（1）：6-12，14.

［41］葉興樹，王偉鋒.閔橋油田阜三段低阻油層成因分析與識別［J］.礦物巖石，2007，27（1）：117-120.

［42］梁春秀.松遼盆地南部低阻油層形成機(jī)理與定量評價［D］.成都：成都理工大學(xué)，2003.

［43］鄧小衛(wèi)，藺志璧，李繼安，等.綜合測井資料在研究新疆吐哈盆地十紅灘鈾礦區(qū)地層巖性及沉積環(huán)境中的應(yīng)用［J］.鈾礦地質(zhì)，2003，19（3）：174-179.

［44］王敏，孫建孟，于華偉，等.G 地區(qū)P 油層砂巖儲層電阻率鈣質(zhì)校正方法研究［J］.石油天然氣學(xué)報，2009，31（1）：61-63.

［45］石油化學(xué)工業(yè)部石油勘探開發(fā)規(guī)劃研究院勘探室編.測井解釋原理與方法［M］.北京：石油化學(xué)工業(yè)出版社，1977：180.

［46］李軍，濮人龍.烏南油田利用沖洗帶電阻率計算碳酸鹽巖含量的新方法［J］.青海石油，2007，25（1）：6-8.

［47］石強(qiáng).鈣質(zhì)砂巖中碳酸鹽含量的計算及校正［J］.測井技術(shù)，1996，20（5）：330-334.

［48］王佩華，彭云彪.2081 砂巖型鈾礦床碳酸鹽含量的計算方法及應(yīng)用［J］.河南理工大學(xué)學(xué)報：自然科學(xué)版，2008，27（3）：282-287.