黃 琳 于增輝 柳 杰

(中海油田服務(wù)股份有限公司油田技術(shù)事業(yè)部 河北三河 065201)

一種基于電成像測(cè)井和雙側(cè)向測(cè)井的聯(lián)合反演方法*

黃 琳 于增輝 柳 杰

(中海油田服務(wù)股份有限公司油田技術(shù)事業(yè)部 河北三河 065201)

利用電阻率測(cè)井理論和反演技術(shù)，以三參數(shù)(侵入帶半徑、侵入帶電阻率、地層電阻率)反演為重點(diǎn)，在常規(guī)雙側(cè)向曲線反演技術(shù)的基礎(chǔ)上，結(jié)合電成像測(cè)井資料，建立了一種基于電成像測(cè)井和雙側(cè)向測(cè)井的聯(lián)合反演方法，即利用電成像測(cè)井探測(cè)深度淺且分辨率高的特性為雙側(cè)向反演提供精細(xì)的地層界面、產(chǎn)狀和侵入帶電阻率信息；運(yùn)用雙側(cè)向測(cè)井資料進(jìn)行迭代反演進(jìn)一步得到侵入帶半徑和地層電阻率參數(shù)。實(shí)際測(cè)井?dāng)?shù)據(jù)的處理結(jié)果表明，本文提出的聯(lián)合反演方法提高了雙側(cè)向反演結(jié)果及測(cè)井解釋的精度，進(jìn)一步提高了常規(guī)測(cè)井資料對(duì)油氣儲(chǔ)層的識(shí)別能力。

電阻率測(cè)井；電成像測(cè)井；雙側(cè)向測(cè)井；聯(lián)合反演；測(cè)井解釋；儲(chǔ)層識(shí)別

在電法測(cè)井領(lǐng)域，地球物理反演是根據(jù)電法儀器在井下觀測(cè)到的物理現(xiàn)象推測(cè)地層內(nèi)部介質(zhì)電性參數(shù)分布狀態(tài)。電法儀器在目的層段的測(cè)井響應(yīng)信息必然包含著井眼、鉆井液侵入以及圍巖層等環(huán)境因素的影響，須進(jìn)一步做校正處理，這是地球物理反演的重要條件。隨著重點(diǎn)勘探井中電測(cè)井曲線的多樣化，利用多種儀器的測(cè)量信息進(jìn)行最優(yōu)化組合成為了可能。不同電法儀器的響應(yīng)機(jī)制雖有差異，但由這些響應(yīng)推斷的地層模型應(yīng)該是相同的;因此，所謂的聯(lián)合反演方法是指利用不同物理機(jī)制的2種或者2種以上的電測(cè)井響應(yīng)曲線進(jìn)行電阻率反演計(jì)算[1-3]。采用聯(lián)合反演方法對(duì)同一地層模型進(jìn)行約束，能夠更好地限制反演過(guò)程中的多解問(wèn)題，增加反演結(jié)果的可靠性。成志剛 等[4]和趙明 等[5]建立了基于高分辨率雙側(cè)向測(cè)井和雙感應(yīng)測(cè)井?dāng)?shù)據(jù)的聯(lián)合反演算法，為提取侵入半徑、侵入帶電阻率和地層電阻率提供了有效手段。劉振華 等[6]建立了適用于感應(yīng)-橫向組合測(cè)井響應(yīng)的聯(lián)合反演算法，研究結(jié)果表明聯(lián)合反演改善了反演結(jié)果的局部收斂性，增加了反演結(jié)果的可信度。國(guó)外研究人員[7-8]同樣給出了相關(guān)的電阻率聯(lián)合反演方法，進(jìn)一步改善了電阻率測(cè)井解釋評(píng)價(jià)。但是，上述參與聯(lián)合反演的電阻率曲線分辨率較低，在利用測(cè)井資料實(shí)現(xiàn)精細(xì)地質(zhì)模型反演方面有待進(jìn)一步提高。

筆者提出了一種新的聯(lián)合反演方法，即基于電成像和雙側(cè)向測(cè)井?dāng)?shù)據(jù)進(jìn)行地層三參數(shù)(侵入帶半徑r、侵入帶電阻率RXO和地層電阻率RT)的聯(lián)合反演方法。常規(guī)的雙側(cè)向迭代反演須使用深、中、淺等3條視電阻率曲線，其中淺電阻率曲線通常使用微球或者微側(cè)向儀器測(cè)量曲線替代，由3條已知曲線迭代求解3個(gè)地層參數(shù)。電成像測(cè)井儀器相比微側(cè)向或微球儀器具有更高的分辨率且徑向探測(cè)深度淺，通過(guò)相關(guān)的數(shù)據(jù)處理方法可以得到侵入帶電阻率信息以及更加精細(xì)的地層層界面、傾角、傾向等信息；因此，基于上述地層信息構(gòu)建井下地層模型，為后續(xù)的雙側(cè)向三維反演提供必要的輸入?yún)?shù)，提高了反演結(jié)果的精度。實(shí)際井的水基電成像和雙側(cè)向測(cè)井資料處理效果驗(yàn)證了本文方法的可靠性和有效性。

1 電成像數(shù)據(jù)處理方法

1.1 侵入帶電阻率曲線合成

雙側(cè)向三維反演須提供侵入帶電阻率作為輸入?yún)?shù)，將這個(gè)待反演參數(shù)作為已知參數(shù)，可減少待反演參數(shù)的數(shù)量，提高反演速度。由于電成像儀器徑向探測(cè)深度淺，其鈕扣電極測(cè)井響應(yīng)主要反映侵入帶電阻率信息，因此侵入帶電阻率RXO可以從電成像數(shù)據(jù)中獲得。本文研究的電成像儀器具有6個(gè)獨(dú)立均勻分布在井筒周向的極板結(jié)構(gòu)，每個(gè)極板上有15個(gè)測(cè)量鈕扣電極。針對(duì)每個(gè)極板得到的15條測(cè)量曲線，采用相關(guān)處理算法在單個(gè)極板內(nèi)合成為1條電測(cè)曲線(圖1)。由圖1可知，合成RXO曲線主要包含深度對(duì)齊、傾角校正、消除局部異常和曲線合成等。經(jīng)過(guò)上述處理，獲得代表不同方位的6條RXO曲線。雙側(cè)向反演中只需要1條能夠反映井周不同方位的綜合侵入帶電阻率響應(yīng)曲線；因此，采用相同的合成方法將6條RXO曲線合成為1條綜合反映井周侵入帶信息的電阻率曲線。

圖1 侵入帶電阻率RXO合成方法Fig .1 Synthetic method of invading zone resistivity RXO

1.2 地層劃分及結(jié)構(gòu)建模

依據(jù)電測(cè)井曲線對(duì)地層進(jìn)行劃分是進(jìn)行油、氣、水層識(shí)別以及測(cè)井解釋評(píng)價(jià)的基礎(chǔ)。雙側(cè)向儀器縱向分辨率較低，只能粗略地對(duì)地層進(jìn)行劃分；這是由儀器自身結(jié)構(gòu)以及探測(cè)特性所決定的。針對(duì)上述不足，提出利用電成像測(cè)井資料進(jìn)行地層劃分及結(jié)構(gòu)建模；這是因?yàn)殡姵上駜x器采用鈕扣電極探頭測(cè)量，能夠提供高分辨率的井壁測(cè)量曲線，利用這些測(cè)量曲線可以實(shí)現(xiàn)地層的精細(xì)劃分及結(jié)構(gòu)建模(圖2)。利用活度分層法[9]和相關(guān)對(duì)比法[10]對(duì)侵入帶曲線及其他輸入曲線進(jìn)行處理，可以得到地層分層數(shù)據(jù)、地層傾角曲線和地層傾向曲線；這些地層分層和產(chǎn)狀信息能夠?yàn)楹罄m(xù)的雙側(cè)向反演提供更為準(zhǔn)確的三維地層正演模型。

圖2 地層劃分及結(jié)構(gòu)建模流程圖Fig .2 Flowchart of stratigraphic division and structural modeling

2 聯(lián)合反演方法

采用馬奎特迭代反演算法[11]構(gòu)造最小二乘目標(biāo)函數(shù)。假定y為實(shí)測(cè)的測(cè)井曲線數(shù)據(jù)，f為地層模型參數(shù)正演響應(yīng)值，其目標(biāo)控制函數(shù)可表述為

(1)

對(duì)式(1)在x(0)附近將f(x)展開(kāi)成Taylor級(jí)數(shù)，并忽略δ的二次項(xiàng)和二次以上的項(xiàng)，則有

(2)

(3)

式(2)～(3)中：P為雅克比矩陣；n為待反演參數(shù)個(gè)數(shù)。

使目標(biāo)控制函數(shù)φ達(dá)到最小修正量δ，其中

(4)

(5)

(6)

式(4)～(6)中：η為阻尼因子;I為單位矩陣。

在給定初始模型參數(shù)值x(0)后，可以計(jì)算得到P、A、g，由此求出δ(0)，進(jìn)一步可得

(7)

再由x(1)計(jì)算得到δ(1)，反復(fù)迭代，直至δ(k)滿足設(shè)定的迭代終止條件為止。

在迭代反演過(guò)程中需要用到二維正演模擬程序和三維正演模擬程序[12]，這樣做的目的是為了提高反演速度，其中地層傾角信息由電成像數(shù)據(jù)處理得到。當(dāng)?shù)貙觾A角較小時(shí)(≤25°)，傾角對(duì)反演結(jié)果影響不大，為了節(jié)省反演時(shí)間，則調(diào)用二維正演模擬程序；當(dāng)?shù)貙觾A角較大時(shí)(>25°)，傾角對(duì)反演結(jié)果影響增大，則調(diào)用三維正演模擬程序。

地層傾角臨界值25°的選取是通過(guò)建立傾斜地層模型，進(jìn)行二維和三維正演模擬并分析2種正演結(jié)果相對(duì)誤差得到的。表1為所建立的傾斜地層模型，圖3、4分別為雙側(cè)向深、淺探測(cè)模式隨地層傾角變化時(shí)對(duì)應(yīng)的響應(yīng)值以及二維和三維正演計(jì)算結(jié)果的相對(duì)誤差對(duì)比，其中上部的3條彩色實(shí)線為深、淺探測(cè)模式下目的層內(nèi)的3個(gè)采樣點(diǎn)在三維正演中的響應(yīng)值，點(diǎn)為對(duì)應(yīng)采樣點(diǎn)在二維正演中的響應(yīng)值，而下部黑色實(shí)線為對(duì)應(yīng)地層傾角時(shí)3個(gè)采樣點(diǎn)在二維與三維正演中的響應(yīng)值相對(duì)誤差的平均值。從圖3、4可以看出，當(dāng)目的層傾角≤25°時(shí)，2種正演方式響應(yīng)值的相對(duì)誤差≤5%；當(dāng)目的層傾角>25°時(shí)，2種正演方式響應(yīng)值的相對(duì)誤差>5%。通過(guò)上述仿真分析可知，在電成像和雙側(cè)向聯(lián)合反演中，可取地層傾角25°作為調(diào)用二維或者三維正演程序的判斷標(biāo)準(zhǔn)。本文提出的電成像和雙側(cè)向聯(lián)合反演流程如圖5所示。

表1 傾斜地層模型Table 1 Inclined stratum model

圖3 傾角變化時(shí)深探測(cè)模式下二維與三維正演 計(jì)算結(jié)果對(duì)比Fig .3 Comparison of 2D and 3D numerical results of deep detection mode with dip angle variation

圖4 傾角變化時(shí)淺探測(cè)模式下二維與三維正演 計(jì)算結(jié)果對(duì)比Fig .4 Comparison of 2D and 3D numerical results of shallow detection mode with dip angle variation

圖5 本文提出的電成像與雙側(cè)向測(cè)井聯(lián)合反演流程Fig .5 Flowchart of joint inversion of electric imaging and dual lateral in the paper

3 實(shí)際資料處理效果

在處理實(shí)際井資料時(shí)，采用逐層反演思路進(jìn)行求解相當(dāng)于在添加約束的前提下把全空間尋找最優(yōu)解的問(wèn)題轉(zhuǎn)化為子空間尋優(yōu)問(wèn)題，可以更好、更快地逼近真實(shí)地層參數(shù)。下面主要通過(guò)2口井的水基電成像和雙側(cè)向測(cè)井資料處理驗(yàn)證本文聯(lián)合反演方法的可靠性和有效性(圖6、7)。

由圖6可知，本文提出的聯(lián)合反演方法得到的地層電阻率RT與雙側(cè)向深、淺電阻率RD、RS的變化趨勢(shì)基本一致，侵入帶電阻率RXO與微側(cè)向電阻率RMLL對(duì)應(yīng)關(guān)系較好，且與侵入深度LXO和GR曲線等的對(duì)應(yīng)關(guān)系也較好，而且反演得到的地層參數(shù)模型再正演結(jié)果與RD、RS的重合度較高，這表明本文聯(lián)合反演方法具有較高的可靠性。此外，本文聯(lián)合反演方法在薄層或薄互層儲(chǔ)層具有更為明顯的優(yōu)勢(shì)。如圖7所示，該井段電成像能夠?qū)崿F(xiàn)精細(xì)的地層劃分并考慮到地層傾斜因素，相對(duì)于只利用雙側(cè)向曲線進(jìn)行分層的常規(guī)二維反演，本文三維聯(lián)合反演的精度有所提高，反演結(jié)果更加接近真實(shí)地層參數(shù)。

圖6 X井聯(lián)合反演和再正演結(jié)果Fig .6 Joint inversion and forward modeling results of Well X

圖7 Y井常規(guī)二維反演和三維聯(lián)合反演結(jié)果對(duì)比Fig .7 Comparison of the results of convention 2D inversion and 3D joint inversion of Well Y

為了進(jìn)一步驗(yàn)證本文聯(lián)合反演方法的有效性，對(duì)Z井中帶有解釋成果(表2)的某井段電成像與雙側(cè)向測(cè)井曲線進(jìn)行了聯(lián)合反演分析(圖8)，可以看出：3 429.6～3 435.5 m及3 442.0～3 448.5 m井段GR值較低，反演結(jié)果表現(xiàn)為侵入較深，滲透性較好；3 456～3 465 m井段GR值較高，反演結(jié)果表現(xiàn)為侵入較淺，滲透性較差。由Z井的解釋結(jié)果(表2)可以看出，反演結(jié)果在儲(chǔ)層位置的滲透性都比較好，與解釋成果一致；以上分析說(shuō)明本文聯(lián)合反演結(jié)果是有效的。

表2 Z井測(cè)井解釋成果表Table 2 Log interpretation results of Well Z

圖8 Z井聯(lián)合反演與解釋成果對(duì)比Fig .8 Comparison of the results of the joint inversion and interpretation of Well Z

4 結(jié)論

提出了一種基于電成像測(cè)井和雙側(cè)向測(cè)井的聯(lián)合反演方法。該方法結(jié)合電成像資料的高分辨率特點(diǎn)，構(gòu)建了更接近地層真實(shí)情況的三維地層模型用于正演計(jì)算，實(shí)現(xiàn)了三維聯(lián)合反演精度的提高；通過(guò)正演仿真考察，提出以25°傾斜角作為二維和三維正演的劃分標(biāo)準(zhǔn)，實(shí)現(xiàn)了聯(lián)合反演速度上的提高。通過(guò)對(duì)實(shí)際井的水基電成像和雙側(cè)向測(cè)井資料的聯(lián)合反演處理，驗(yàn)證了該方法的可靠性和有效性。

[1] 譚茂金，石耀霖，王曉杰.多物理場(chǎng)測(cè)井?dāng)?shù)據(jù)聯(lián)合反演研究進(jìn)展[J].地球物理學(xué)進(jìn)展，2008，23(5)：1520-1525. TAN Maojin，SHI Yaolin，WANG Xiaojie.Progress in joint inversion of well logging data of multi-physics fields[J].Progress in Gophysics,2008,23(5):1520-1525.

[2] 何委徽，王家林，于鵬.地球物理聯(lián)合反演研究的現(xiàn)狀與趨勢(shì)分析[J].地球物理學(xué)進(jìn)展，2009，24(2)：530-540. HE Weihui，WANG Jialin，YU Peng.Overview of the status and prospect of geophysical joint inversion[J].Progress in Geophysics,2009,24(2):530-540.

[3] 敬榮中，鮑光淑，陳紹裘.地球物理聯(lián)合反演研究綜述[J].地球物理學(xué)進(jìn)展，2003，18(3)：535-540. JING Rongzhong，BAO Guangshu，CHEN Shaoqiu.A review of the research for geophysical combinative inversion[J].Progress in Geophysics,2003,18(3):535-540.

[4] 成志剛，孫寶佃，劉振華，等.高分辨率雙側(cè)向測(cè)井和雙感應(yīng)測(cè)井?dāng)?shù)據(jù)聯(lián)合反演研究與應(yīng)用[J].測(cè)井技術(shù)，2010，34(6)：542-547. CHENG Zhigang，SUN Baodian，LIU Zhenhua，et al.Joint inversion of the high resolution dual laterolog data with dual induction logging data and its applications[J].Well Logging Technology,2010,34(6):542-547.

[5] 趙明，高杰，孫友國(guó).常規(guī)電測(cè)井聯(lián)合反演研究與實(shí)際應(yīng)用[J].測(cè)井技術(shù)，2003，27(1)：16-19. ZHAO Ming，GAO Jie，SUN Youguo.Joint inversion study of conventional electrical logging and its applications[J].Well Logging Technology,2003,27(1):16-19.

[6] 劉振華，林豪.感應(yīng)-橫向組合測(cè)井聯(lián)合反演及多解性研究[J].地球物理學(xué)進(jìn)展，2001，16(1)：47-55. LIU Zhenhua，LIN Hao.Joint inversion for combining induction and departure logs[J].Progress in Geophysics,2001,16(1):47-55.

[7] SMITS J W，DUBOURG I,LULING M G,et al.Improved resistivity interpretation utilizing a new array tool and associated inversion processing[J].Resistivity Logging,1998,21(3):831-844.

[8] STRACK K M.Joint inversion of resistivity logs：limit and possibilities[J].Journal of Petroleum Technology,1997,49(11):1234-1236.

[9] 易覺(jué)非.利用活度分層法實(shí)現(xiàn)測(cè)井自動(dòng)地質(zhì)分層[J].石油天然氣學(xué)報(bào)，2007，29(1)：78-80. YI Juefei.Automatic geologic zonation using activity layering method[J].Journal of Oil and Gas Technology,2007,29(1):78-80.

[10] 王曰才，王冠貴.地層傾角測(cè)井[M].北京：石油工業(yè)出版社，1982. WANG Yuecai,WANG Guangui.Formation dip logging[M].Beijing:Petroleum Industry Press,1982.

[11] 馮進(jìn)，張中慶，羅虎.隨鉆電磁波電阻率和電纜電阻率測(cè)井聯(lián)合反演及應(yīng)用[J].測(cè)井技術(shù)，2015，39(1)：21-26. FENG Jin，ZHANG Zhongqing，LUO Hu.Joint inversion of electromagnetic wave resistivity logging while drilling(LWD) and wireline resistivity logging and its applications[J].Well Logging Technology,2015,39(1):21-26.

[12] 張庚驥.電測(cè)井算法[M].北京：石油工業(yè)出版社，2007：1-36． ZHANG Gengji.Electric logging algorithm[M].Beijing：Petroleum Industry Press,2007:1-36.

(編輯：張喜林)

A joint inversion method based on electrical imaging logging and dual laterolog

HUANG Lin YU Zenghui LIU Jie

(COSLWell-TechDivision,Sanhe,Hebei065201,China)

The joint inversion method based on electrical imaging logging and dual laterolog is established with the resistivity logging theory and the inversion technique, especially the inversion of three parameters (invasion radius, invasion resistivity and formation resistivity), in which the electrical imaging logging can provide fine stratigraphic interface, attitude of stratum and invasion resistivity for dual laterolog due to its characteristics of small investigation radius and high resolution. By means of an iterative inversion with dual laterolog data, the invasion radius and formation resistivity are obtained. The processing results of the actual logging data show that the proposed joint inversion improves the accuracy of dual laterolog inversion and well logging interpretation, and further improves oil and gas reservoirs identification ability of conventional well logging data.

resistivity logging; electrical imaging logging; dual laterolog; joint inversion; well logging interpretation; reservoir identification

*“十二五”國(guó)家科技重大專項(xiàng)“大型油氣田及煤層氣開(kāi)發(fā)”課題“三維聲波、油基泥漿電成像、二維核磁成像測(cè)井技術(shù)與裝備(編號(hào)：2011ZX05020-005)”部分研究成果。

黃琳，男，高級(jí)工程師，1993年畢業(yè)于原石油大學(xué)(華東)勘查地球物理專業(yè)，獲學(xué)士學(xué)位，主要從事測(cè)井技術(shù)與儀器應(yīng)用研究工作。地址：河北省三河市燕郊開(kāi)發(fā)區(qū)行宮大街中海油服油田技術(shù)事業(yè)部(郵編：065201)。E-mail：huanglin@cosl.com.cn。

1673-1506(2017)03-0046-06

10.11935/j.issn.1673-1506.2017.03.007

P631.8

A

2016-10-19 改回日期：2017-01-11

黃琳,于增輝,柳杰.一種基于電成像測(cè)井和雙側(cè)向測(cè)井的聯(lián)合反演方法[J].中國(guó)海上油氣，2017,29(3)：46-51.

HUANG Lin,YU Zenghui,LIU Jie.A joint inversion method based on electrical imaging logging and dual laterolog[J].China Offshore Oil and Gas,2017,29(3):46-51.