用并矢Green函數(shù)的矢量本征函數(shù)展開式評價(jià)金屬心軸對多分量感應(yīng)測井響應(yīng)的影響

俞燕明,肖加奇,魏寶君,趙彥偉

(1.中石油長城鉆探工程有限公司測井技術(shù)研究院,北京 102206;2.中國石油大學(xué)理學(xué)院,山東青島 266580)

用并矢Green函數(shù)的矢量本征函數(shù)展開式評價(jià)金屬心軸對多分量感應(yīng)測井響應(yīng)的影響

俞燕明1,肖加奇1,魏寶君2,趙彥偉1

(1.中石油長城鉆探工程有限公司測井技術(shù)研究院,北京 102206;2.中國石油大學(xué)理學(xué)院,山東青島 266580)

基于圓柱坐標(biāo)系下徑向成層介質(zhì)中并矢Green函數(shù)的矢量本征函數(shù)展開式分析金屬心軸對多分量感應(yīng)測井響應(yīng)的影響。為提高計(jì)算精度,考慮各分量線圈系的具體形狀。將金屬心軸作為一層介質(zhì)處理,既可以考慮其電導(dǎo)率有限、也可以考慮其電導(dǎo)率為無窮大的情況。計(jì)算結(jié)果表明,金屬心軸的存在導(dǎo)致各分量感應(yīng)電動勢實(shí)部信號和虛部信號的強(qiáng)度變小。當(dāng)心軸電導(dǎo)率較小時(shí)無論是實(shí)部信號還是虛部信號均不隨心軸電導(dǎo)率的增加而明顯改變。當(dāng)心軸電導(dǎo)率增加到一定臨界值時(shí)感應(yīng)電動勢的實(shí)部信號快速減小并逐漸穩(wěn)定到一個(gè)較低數(shù)值,而虛部信號先逐漸增加再快速減小并逐漸穩(wěn)定到一個(gè)較低數(shù)值。

多分量感應(yīng)測井;金屬心軸;并矢Green函數(shù);矢量本征函數(shù)

多分量感應(yīng)測井儀器[1-5]由多組沿三個(gè)方向彼此垂直的發(fā)射-接收線圈陣列組成,可同時(shí)探測各向異性地層如砂-泥巖薄交互層的水平電阻率和垂向電阻率信息,對復(fù)雜地層尤其是傾斜井或水平井情況下各向異性地層的高精度對比評價(jià)具有重要意義。對于多分量感應(yīng)測井儀器在復(fù)雜地層中響應(yīng)的模擬,已有很多文獻(xiàn)進(jìn)行了分析討論[6-11],但這些數(shù)值模擬對儀器結(jié)構(gòu)考慮的較為簡單,一般將各分量的線圈系作為磁偶極子處理。在多分量感應(yīng)測井儀器的設(shè)計(jì)制造過程中,為增加儀器的機(jī)械強(qiáng)度,在儀器中心一般置有金屬心軸。由于金屬心軸電導(dǎo)率較高,各分量發(fā)射線圈的電流在其表面產(chǎn)生渦流,金屬心軸的渦流與發(fā)射線圈的電流共同在接收線圈處產(chǎn)生感應(yīng)電動勢,從而對接收線圈處的電磁信號產(chǎn)生影響。因而有必要分析金屬心軸對多分量感應(yīng)測井響應(yīng)的影響,并在數(shù)值模擬中考慮線圈系的具體結(jié)構(gòu),從而為實(shí)際儀器的刻度提供準(zhǔn)確的刻度值。文獻(xiàn)[12-13]討論了金屬心軸對僅含軸向線圈系的陣列感應(yīng)測井儀器響應(yīng)的影響。筆者基于圓柱坐標(biāo)系下徑向成層介質(zhì)中并矢Green函數(shù)的矢量本征函數(shù)展開式討論金屬心軸對多分量感應(yīng)測井儀器不同分量線圈系響應(yīng)的影響規(guī)律,計(jì)算考慮金屬心軸和線圈系形狀時(shí)各分量感應(yīng)電動勢隨地層電導(dǎo)率變化的規(guī)律。

1 基本理論

由于多分量感應(yīng)測井儀器心軸的基本形狀為圓柱形,并考慮到地層呈軸對稱分布,故采用圓柱坐標(biāo)系。設(shè)發(fā)射源隨時(shí)間的變化關(guān)系為exp(-iωt),其中ω為角頻率,并假設(shè)場點(diǎn)位置坐標(biāo)為R(r,φ,z)、源點(diǎn)位置坐標(biāo)為R′(r′,φ′,z′)。

1.1 均勻介質(zhì)中并矢Green函數(shù)的矢量本征函數(shù)展開式

均勻介質(zhì)中電流源電型并矢Green函數(shù)PGEJ(R,R′)的矢量本征函數(shù)展開式同PGHM(R,R′),電流源磁型并矢Green函數(shù)PGHJ(R,R′)的矢量本征函數(shù)展開式同-PGEM(R,R′)。

1.2 徑向成2層介質(zhì)中并矢Green函數(shù)的矢量本征函數(shù)展開式

考慮到實(shí)際需要,首先假設(shè)圓柱形徑向成層介質(zhì)由內(nèi)向外在徑向共有2層(圖1),層界面半徑為r1,每層介質(zhì)的參數(shù)為2μlεl,ηl=(l=1,2)?？紤]到多分量感應(yīng)測井儀器的發(fā)射線圈和接收線圈均在心軸外側(cè),只列出源點(diǎn)位置在第2層,場點(diǎn)位置在不同地層時(shí)磁流源和電流源并矢Green函數(shù)的矢量本征函數(shù)展開式。

圖1 線圈系模型簡圖Fig.1 Schematic configuration of coils

當(dāng)場點(diǎn)位置在第1層時(shí):

式中,、、為待定系數(shù)。

當(dāng)場點(diǎn)位置在第2層時(shí):

若將金屬心軸視為理想導(dǎo)體,其電導(dǎo)率為無窮大,則電磁場只存在于介質(zhì)2中,對磁流源并矢Green函數(shù)而言,圓柱界面r1處^er×(R,R′)=0,則

對電流源并矢Green函數(shù)而言,圓柱界面r1處(R,R′)=0,則

無論是金屬心軸的電導(dǎo)率有限還是無窮大,各待定系數(shù)滿足如下關(guān)系:

1.3 多分量感應(yīng)測井響應(yīng)的計(jì)算

多分量感應(yīng)測井儀器的z向線圈可視為圓柱面,x、y方向的線圈可視為長方形(圖1),因而在計(jì)算多分量感應(yīng)測井響應(yīng)時(shí)可分4種情況。

1.3.1z向發(fā)射、z向接收(即zz分量)

設(shè)發(fā)射線圈半徑為rT,中心點(diǎn)軸向坐標(biāo)為zT,匝數(shù)為NT,電流強(qiáng)度為IT,接收線圈半徑為rR,中心點(diǎn)軸向坐標(biāo)為zR,匝數(shù)為NR。發(fā)射線圈沿軸向均勻排列在軸向?qū)挾葹棣T的圓柱面上,接收線圈沿軸向均勻排列在軸向?qū)挾葹棣R的圓柱面上,則發(fā)射線圈上的電流密度可表示為

而地層中的電磁場強(qiáng)度可表示為

式中,D為電流源所在的區(qū)域。將式(3)代入式(4),得

接收線圈處的感應(yīng)電動勢可表示為

將式(6)代入式(7),得

根據(jù)并矢Green函數(shù)的展開特點(diǎn),式(8)的各項(xiàng)積分均可以通過解析方式得到,且并矢Green函數(shù)只含有n=0階項(xiàng)。

1.3.2z向發(fā)射、x或y向接收(即zx分量或zy分量)

設(shè)接收線圈中心點(diǎn)坐標(biāo)為(rR,φR,zR)、高度為ΔzR、寬度為ΔlR,在垂向坐標(biāo)為z、局部橫向位置為l處取一高度為dz、寬度為dl的長方形面元ds,則該面元在圓柱坐標(biāo)系下的坐標(biāo)為

該小面元處的總磁通為

接收線圈的感應(yīng)電動勢為

式(9)中的磁場強(qiáng)度可通過將式(3)代入式(5)得到,即

由于發(fā)射線圈產(chǎn)生的磁場沒有φ分量,將式(10)代入式(9)得

式中除第三重對dl的積分需用數(shù)值計(jì)算外,其他積分均可以通過解析方式得到,且并矢Green函數(shù)只含有n=0階項(xiàng)。

1.3.3x或y向發(fā)射、z向接收(即xz分量或yz分量)

設(shè)發(fā)射線圈中心點(diǎn)坐標(biāo)為(rT,φT,zT)、高度為ΔzT、寬度為ΔlT,在垂向坐標(biāo)為z′、局部橫向位置為l′處取一高度為dz′、寬度為dl′的長方形面元ds′,則該面元在圓柱坐標(biāo)系下的坐標(biāo)為

該面元所繞線圈的磁偶極矩可表示為

而地層中的電磁場強(qiáng)度可表示為

結(jié)合式(12)～(14),得到長方形發(fā)射線圈在地層中產(chǎn)生的電磁場強(qiáng)度為

接收線圈處的感應(yīng)電動勢可表示為

將式(15)代入式(17),得

式中除第一重對dl′的積分需用數(shù)值計(jì)算外,其他積分均可以通過解析方式得到,且并矢Green函數(shù)只含有n=0階項(xiàng)。

1.3.4x或y向發(fā)射、x或y向接收(即xx、xy、yx或yy分量)

將式(16)代入式(9),得

式中,第二重和第四重對dz′、dz的積分可解析得到,第一重和第三重對dl′、dl的積分需用數(shù)值計(jì)算。

2 數(shù)值算例

2.1 多分量感應(yīng)測井響應(yīng)與均質(zhì)地層電導(dǎo)率的關(guān)系

假設(shè)發(fā)射線圈電流IT=1 A、金屬心軸半徑r1= 0.0127 m,分別考察心軸電導(dǎo)率σ1為0和無窮大時(shí)各分量感應(yīng)電動勢隨均質(zhì)地層電導(dǎo)率σ2的變化關(guān)系。

2.1.1zz分量的感應(yīng)電動勢

z向發(fā)射線圈由兩部分組成,匝數(shù)NT=50、半徑rT=0.033 m、軸向?qū)挾圈T=0.052 m、中心點(diǎn)軸向坐標(biāo)分別為zT=±0.0565 m。z向接收線圈由主接收線圈和屏蔽接收線圈兩部分組成,rR=rT。主接收線圈匝數(shù)NR=90、中心點(diǎn)軸向坐標(biāo)zR=1.143 m、軸向?qū)挾圈R=0.09434 m,屏蔽接收線圈匝數(shù)N′R=-40、中心點(diǎn)軸向坐標(biāo)z′R=0.8737 m、軸向?qū)挾圈′R=0.04134 m。取發(fā)射頻率為44 kHz。圖2給出了zz分量感應(yīng)電動勢隨均質(zhì)地層電導(dǎo)率變化關(guān)系。該分量亦采用文獻(xiàn)[18]的方法進(jìn)行了模擬計(jì)算,兩種方法的計(jì)算結(jié)果完全一致,說明了本文方法的正確性。

圖2zz分量感應(yīng)電動勢隨均質(zhì)地層電導(dǎo)率變化關(guān)系Fig.2 Relationship betweenzz-component's electromotive force and homogeneous formation's conductivity

2.1.2zx分量的感應(yīng)電動勢

所采用頻率和z向發(fā)射線圈參數(shù)同上。x向接收線圈系的高度ΔzR=0.036 m、寬度ΔlR=0.04326 m、中心點(diǎn)徑向坐標(biāo)rR=0.02725 m。接收線圈系由方位角分別在φR=0°和φR=180°的兩對完全相同的主接收/屏蔽接收線圈組成。主接收線圈匝數(shù)NR=44、中心點(diǎn)軸向坐標(biāo)zR=1.6 m,屏蔽接收線圈匝數(shù)N′R=-22、中心點(diǎn)軸向坐標(biāo)z′R=1.27 m。圖3給出了zx分量感應(yīng)電動勢隨均質(zhì)地層電導(dǎo)率變化關(guān)系。

圖3zx分量感應(yīng)電動勢隨均質(zhì)地層電導(dǎo)率變化關(guān)系Fig.3 Relationship betweenzx-component's electromotive force and homogeneous formation's conductivity

2.1.3xz分量的感應(yīng)電動勢

取發(fā)射頻率為64 kHz,z向接收線圈參數(shù)同上。x向發(fā)射線圈系由方位角分別在φT=0°和φT=180°處完全相同的兩個(gè)線圈組成,匝數(shù)NT=39、中心點(diǎn)軸向坐標(biāo)zT=0.0 m、中心點(diǎn)徑向坐標(biāo)rT=0.02725 m、高度ΔzT=0.036 m、寬度ΔlT=0.04326 m。圖4給出了xz分量感應(yīng)電動勢隨均質(zhì)地層電導(dǎo)率變化關(guān)系。

圖4xz分量感應(yīng)電動勢隨均質(zhì)地層電導(dǎo)率變化關(guān)系Fig.4 Relationship betweenxz-component's electromotive force and homogeneous formation's conductivity

2.1.4xx分量的感應(yīng)電動勢

取發(fā)射頻率為64 kHz,x向發(fā)射線圈系和x向接收線圈系結(jié)構(gòu)及參數(shù)同上。圖5給出了xx分量感應(yīng)電動勢隨均質(zhì)地層電導(dǎo)率變化關(guān)系。

圖5xx分量感應(yīng)電動勢隨均質(zhì)地層電導(dǎo)率變化關(guān)系Fig.5 Relationship betweenxx-component's electromotive force and homogeneous formation's conductivity

由圖2～5可以看出,由于金屬心軸的存在使發(fā)射線圈和接收線圈的有效面積減小,故各分量感應(yīng)電動勢實(shí)部和虛部的絕對值亦相應(yīng)減小,且不同分量感應(yīng)電動勢減小的規(guī)律相同。當(dāng)然,不同線圈距、不同分量感應(yīng)電動勢的數(shù)值及其隨均質(zhì)地層電導(dǎo)率的變化關(guān)系是有差別的。計(jì)算發(fā)現(xiàn),對同一分量而言,線圈距越大,信號強(qiáng)度越小,并且在電導(dǎo)率越小時(shí)感應(yīng)電動勢開始出現(xiàn)震蕩,說明對地層電導(dǎo)率的探測范圍變小。對比不同分量的感應(yīng)電動勢發(fā)現(xiàn),xx分量對地層電導(dǎo)率的探測范圍最小。由于已考慮到各分量線圈系的具體形狀和金屬心軸的影響,該儀器模型更接近于實(shí)際儀器的結(jié)構(gòu),因而圖2～5的模擬數(shù)據(jù)更接近于儀器測量值,從而可為實(shí)際儀器提供更高精度的刻度值。

2.2 多分量感應(yīng)測井響應(yīng)與鉆井液電導(dǎo)率的關(guān)系

在存在井眼的情況下,采用徑向成3層介質(zhì)模型,第1層介質(zhì)為金屬心軸、第2層介質(zhì)為井眼、第3層介質(zhì)為原狀均質(zhì)地層。設(shè)井眼半徑r2=0.12 m、原狀地層電導(dǎo)率σ3=0.2 S·m-1,線圈系參數(shù)同上。圖6給出了部分分量的感應(yīng)電動勢隨井眼內(nèi)鉆井液電導(dǎo)率σm的變化關(guān)系。由圖6可以得到相同結(jié)論,即在存在井眼的情況下各分量感應(yīng)電動勢實(shí)部和虛部的絕對值均相應(yīng)減小。

圖6 不同分量感應(yīng)電動勢隨鉆井液電導(dǎo)率變化關(guān)系Fig.6 Relationship between different component's electromotive force and drilling mud's conductivity

2.3 多分量感應(yīng)測井響應(yīng)與心軸電導(dǎo)率的關(guān)系

為了考察心軸電導(dǎo)率的不同取值對各分量感應(yīng)電動勢的影響,固定均質(zhì)地層電導(dǎo)率σ2=0.2 S· m-1,取發(fā)射頻率分別為44和68 kHz。圖7給出了zz分量接收線圈系的感應(yīng)電動勢隨心軸電導(dǎo)率的變化關(guān)系。圖8給出了zx分量接收線圈系的感應(yīng)電動勢隨心軸電導(dǎo)率的變化關(guān)系。對其他分量可得到相似規(guī)律。

圖7zz分量感應(yīng)電動勢隨心軸電導(dǎo)率變化關(guān)系Fig.7 Relationship betweenzz-component's electromotive force and mandrel's conductivity

圖8zx分量感應(yīng)電動勢隨心軸電導(dǎo)率變化關(guān)系Fig.8 Relationship betweenzx-component's electromotive force and mandrel's conductivity

由圖7～8可以看出,各分量實(shí)部感應(yīng)電動勢與虛部感應(yīng)電動勢隨心軸電導(dǎo)率的變化關(guān)系是不同的。對各分量實(shí)部感應(yīng)電動勢而言,當(dāng)心軸電導(dǎo)率較小時(shí)其數(shù)值并不隨心軸電導(dǎo)率的增加而明顯改變,只有當(dāng)心軸電導(dǎo)率增加到一定臨界值時(shí)感應(yīng)電動勢數(shù)值快速減小,此后逐漸穩(wěn)定到一個(gè)較低數(shù)值。對各分量虛部感應(yīng)電動勢而言,當(dāng)心軸電導(dǎo)率較小時(shí)其數(shù)值亦不隨心軸電導(dǎo)率的增加而明顯改變,當(dāng)心軸電導(dǎo)率增加到一定臨界值時(shí)感應(yīng)電動勢數(shù)值先逐漸增加再快速減小,此后逐漸穩(wěn)定到一個(gè)較低數(shù)值。心軸電導(dǎo)率的臨界值與頻率有關(guān),頻率越高,心軸電導(dǎo)率的臨界值越小。

3 結(jié) 論

(1)采用基于圓柱坐標(biāo)系下徑向成層介質(zhì)并矢Green函數(shù)的矢量本征函數(shù)展開式可有效模擬含金屬心軸多分量感應(yīng)測井儀器在徑向成層地層中的響應(yīng)。由于在模擬中考慮了金屬心軸對儀器響應(yīng)的影響,不僅能為多分量感應(yīng)測井刻度裝置設(shè)計(jì)提供理論依據(jù),而且能用于高精度的刻度器地層等效讀數(shù)的計(jì)算。

(2)當(dāng)金屬心軸存在時(shí),多分量感應(yīng)測井儀器各分量感應(yīng)電動勢實(shí)部信號和虛部信號的強(qiáng)度變小,且不同分量感應(yīng)電動勢減小的規(guī)律相同。

(3)各分量感應(yīng)電動勢的實(shí)部信號和虛部信號隨心軸電導(dǎo)率的變化關(guān)系不同。當(dāng)心軸電導(dǎo)率較小時(shí)無論是實(shí)部信號還是虛部信號均不隨心軸電導(dǎo)率的增加而明顯改變。當(dāng)心軸電導(dǎo)率增加到一定臨界值時(shí),感應(yīng)電動勢的實(shí)部信號快速減小,并逐漸穩(wěn)定到一個(gè)較低數(shù)值,而虛部信號先逐漸增加再快速減小,并逐漸穩(wěn)定到一個(gè)較低數(shù)值。

[1] KRIEGSHAUSER B,FANINI O,FORGANG S,et al.A new multi-component induction logging tool to resolve anisotropic formations[C/CD].SPWLA 41st Annual Logging Symposium,2000,paper D.

[2] KRIEGSHAUSER B,FANINI O,YU L,et al.Improved shaly sand interpretation in highly deviated and horizontal wells using multicomponent induction log data[C/CD]. SPWLA 42nd Annual Logging Symposium,2001,paper S.

[3] ZHDANOV M,KENNEDY D,PEKSEN E.Foundations of tensor induction well-logging[J].Petrophysics,2001, 42(6):588-610.

[4] ZHANG Z,YU L,KRIEGSHAUSER B,et al.Simultaneous determination of relative angles and anisotropic resistivity using multicomponent induction logging data[C/ CD].SPWLA 42nd Annual Logging Symposium,2001, paper Q.

[5] MALLAN R K,TORRES-VERDIN C.Effects of petrophysical,environmental,and geometrical parameters on multi-component induction measurements acquired in high-angle wells[C/CD].SPWLA 47th Annual Logging Symposium,2006,paper PPP.

[6] 汪功禮,張庚驥,崔鋒修,等.三維感應(yīng)測井響應(yīng)計(jì)算的交錯(cuò)網(wǎng)格有限差分法[J].地球物理學(xué)報(bào),2003, 46(4):561-567.

WANG Gongli,ZHANG Gengji,CUI Fengxiu,et al. Application of staggered grid finite difference method to the computation of 3-D induction logging response[J]. Chinese Journal of Geophysics,2003,46(4):561-567.

[7] 王昌學(xué),楊韋,儲昭坦,等.多分量感應(yīng)測井響應(yīng)的交錯(cuò)網(wǎng)格有限差分法模擬[J].石油大學(xué)學(xué)報(bào):自然科學(xué)版,2005,29(3):35-40.

WANG Changxue,YANG Wei,CHU Zhaotan,et al. Modeling of multicomponent induction logging response by staggered grid finite difference method[J].Journal of the University of Petroleum,China(Edition of Natural Science),2005,29(3):35-40.

[8] WANG H N,SO P,YANG S W,et al.Numerical modeling of multicomponent induction well-logging tools in the crlindrically stratified anisotropic media[J].IEEE Trans on Geosci and Remote Sensing,2008,46(4):1134-1147.

[9] 汪宏年,陶宏根,姚敬金,等.用模式匹配算法研究層狀各向異性傾斜地層中多分量感應(yīng)測井響應(yīng)[J].地球物理學(xué)報(bào),2008,51(5):1591-1599.

WANG Hongnian,TAO Honggen,YAO Jingjin,et al. Study on the response of a multicomponent induction logging tool in deviated and layered anisotropic formations by using numerical mode matching method[J].Chinese Journal of Geophysics,2008,51(5):1591-1599.

[10] 魏寶君,王甜甜,王穎.用磁流源并矢Green函數(shù)的遞推矩陣方法計(jì)算層狀各向異性地層中多分量感應(yīng)測井響應(yīng)[J].地球物理學(xué)報(bào),2009,52(11):2920-2928.

WEI Baojun,WANG Tiantian,WANG Ying.Computing the response of multi-component induction logging in layered anisotropic formation by the recursive matrix method for magnetic-current-source dyadic Green's function[J].Chinese Journal of Geophysics,2009,52 (11):2920-2928.

[11] 張國艷,肖加奇,肖占山,等.三維感應(yīng)測井儀器在三維井眼環(huán)境中的正演模擬[J].中國石油大學(xué)學(xué)報(bào):自然科學(xué)版,2013,37(3):63-67.

ZHANG Guoyan,XIAO Jiaqi,XIAO Zhanshan,et al. Forward modeling of 3D induction logging tool in 3D borehole environment[J].Journal of China University of Petroleum(Edition of Natural Science),2013,37(3): 63-67.

[12] 魏寶君,張庚驥,張中慶,等.金屬心棒對陣列感應(yīng)測井儀器的影響[J].測井技術(shù),1998,22(4):245-249.

WEI Baojun,ZHANG Gengji,ZHANG Zhongqing,et al.Influence of the metal mandrel on the array induction logging tool[J].Well Logging Technology,1998,22 (4):245-249.

[13] 魏寶君,張庚驥.含金屬心棒陣列感應(yīng)測井儀器的井眼校正[J].石油大學(xué)學(xué)報(bào):自然科學(xué)版,1999,23 (3):22-26.

WEI Baojun,ZHANG Gengji.Borehole correction of array induction logging tool with metal mandrel[J].Journal of the University of Petroleum,China(Edition of Natural Science),1999,23(3):22-26.

[14] 戴振鐸,魯述.電磁理論中的并矢格林函數(shù)[M].武漢:武漢大學(xué)出版社,1996.

[15] 魏寶君.一種新型隨鉆電阻率測井儀器的響應(yīng)和刻度[J].地球物理學(xué)報(bào),2007,50(2):632-641.

WEI Baojun.Response and calibration of a new loggingwhile-drilling resistivity tool[J].Chinese Journal of Geophysics,2007,50(2):632-641.

[16] 魏寶君,田坤,張旭,等.用并矢Green函數(shù)的矢量本征函數(shù)展開式評價(jià)偏心對隨鉆電磁波電阻率測井響應(yīng)的影響[J].中國石油大學(xué)學(xué)報(bào):自然科學(xué)版, 2010,34(5):57-62.

WEI Baojun,TIAN Kun,ZHANG Xu,et al.Evaluating influence of eccentricity on response of electromagnetic wave resistivity logging-while-drilling by vector eigenfunction expansion formulae for dyadic Green's functions[J].Journal of China University of Petroleum(E-dition of Natural Science),2010,34(5):57-62.

[17] 魏寶君,王莎莎,歐永峰.用并矢Green函數(shù)的矢量本征函數(shù)展開式模擬隨鉆感應(yīng)測井儀器的響應(yīng)[J].地球物理學(xué)報(bào),2011,54(5):1391-1401.

WEI Baojun,WANG Shasha,OU Yongfeng.Simulating the response of induction logging-while-drilling tools by vector eigenfunction expansion formulae for dyadic Green 's functions[J].Chinese Journal of Geophysics,2011, 54(5):1391-1401.

[18] 魏寶君,王穎,王甜甜.電磁波電阻率儀器的基本理論及其在隨鉆測量中的應(yīng)用[J].地球物理學(xué)進(jìn)展, 2009,24(2):774-781.

WEI Baojun,WANG Ying,WANG Tiantian.The basic theory of electromagnetic wave resistivity instrument and its application to measurement of while-drilling[J].Progress in Geophysics,2009,24(2):774-781.

(編輯 修榮榮)

Evaluating influence of metal mandrel on response of multi-component induction logging by vector eigenfunction expansion formulae for dyadic Green's functions

YU Yanming1,XIAO Jiaqi1,WEI Baojun2,ZHAO Yanwei1
(1.Logging Technique Research Institute of Great Wall Drilling Company in China National Petroleum Corporation, Beijing 102206,China; 2.College of Science in China University of Petroleum,Qingdao 266580,China)

The influence of metal mandrel on the response of multi-component induction logging was analyzed using the vector eigenfunction expansion formulae for dyadic Green's functions of radial-layered medium in cylindrical coordinate system.The shape of each component's coils was taken into account in order to increase the computational precision.The metal mandrel was assumed as a layer of medium whose conductivity can be both finite and infinite.The results show that the intensity of both the real and imaginary parts of each component's electromotive force will decrease because of the introduction of the metal mandrel.There is no distinct change in both the real and imaginary parts of the electromotive force when the amount of increase of the mandrel's conductivity is small.Both the real and imaginary parts of the electromotive force decrease rapidly when the mandrel's conductivity reaches a critical value,and then the real part will gradually achieve a stable lower value. While the imaginary part will first increase and then decrease when the mandrel's conductivity reaches the critical value,it will also gradually achieve a stable lower value eventually.

multi-component induction logging;metal mandrel;dyadic Green's functions;vector eigenfunction

P 631.9

A

1673-5005(2014)04-0057-08

10.3969/j.issn.1673-5005.2014.04.008

2013-10-05

國家“十二五”重大專項(xiàng)(2011ZX05020-004)

俞燕明(1983-),男,工程師,博士,研究方向?yàn)殡姺y井與核磁共振理論方法。E-mail:yuyanming@cnlc.cn。

俞燕明,肖加奇,魏寶君,等.用并矢Green函數(shù)的矢量本征函數(shù)展開式評價(jià)金屬心軸對多分量感應(yīng)測井響應(yīng)的影響[J].中國石油大學(xué)學(xué)報(bào):自然科學(xué)版,2014,38(4):57-64.

YU Yanming,XIAO Jiaqi,WEI Baojun,et al.Evaluating influence of metal mandrel on response of multi-component induction logging by vector eigenfunction expansion formulae for dyadic Green's functions[J].Journal of China University of Petroleum(Edition of Natural Science),2014,38(4):57-64.