隨鉆方位電磁波儀器探測(cè)電阻率各向異性新方法

楊 震， 文 藝， 肖紅兵

(1.中石化勝利石油工程有限公司鉆井工藝研究院，山東東營(yíng)257017；2.中石化勝利石油工程有限公司測(cè)井公司，山東東營(yíng) 257096)

?

隨鉆方位電磁波儀器探測(cè)電阻率各向異性新方法

楊 震1， 文 藝2， 肖紅兵1

(1.中石化勝利石油工程有限公司鉆井工藝研究院，山東東營(yíng)257017；2.中石化勝利石油工程有限公司測(cè)井公司，山東東營(yíng) 257096)

隨著隨鉆方位電磁波儀器在大斜度井和水平井中的廣泛應(yīng)用，電阻率各向異性已成為影響地質(zhì)導(dǎo)向和地層評(píng)價(jià)準(zhǔn)確性的主要因素之一。以Baker Hughes公司的APR儀器為例，根據(jù)隨鉆方位電磁波儀器多分量的特點(diǎn)，采用數(shù)值模擬方法分析電阻率各向異性對(duì)儀器各分量信號(hào)響應(yīng)的影響，利用對(duì)稱發(fā)射補(bǔ)償測(cè)量來(lái)增強(qiáng)或消除各向異性的影響，并采用正反演方法確定地層電阻率各向異性。模擬結(jié)果表明:磁場(chǎng)信號(hào)軸向分量與電阻率各向異性在相對(duì)井斜角0° ～90°內(nèi)呈單調(diào)遞增關(guān)系；磁場(chǎng)信號(hào)橫向分量與電阻率各向異性呈拋物線規(guī)律，在相對(duì)井斜角為0°和90°時(shí)影響為零；隨鉆方位電磁波電阻率儀器不同分量進(jìn)行組合可以確定地層電阻率的各向異性。對(duì)方位電磁波響應(yīng)曲線數(shù)據(jù)進(jìn)行三參數(shù)反演得到地層水平電阻率、垂直電阻率以及相對(duì)井斜角，利用反演后的測(cè)量資料可以提高地層評(píng)價(jià)和地質(zhì)導(dǎo)向的準(zhǔn)確性。

方位電磁波；電阻率各向異性；補(bǔ)償；反演；地質(zhì)導(dǎo)向；地層評(píng)價(jià)

隨著油田進(jìn)入開發(fā)后期，原來(lái)認(rèn)為沒有工業(yè)開采價(jià)值的小油層、薄油層、斷塊油層和老油田衰竭剩余油藏被重新開發(fā)利用。 據(jù)估計(jì)，有大約30%的油氣存儲(chǔ)在低阻、低對(duì)比度的各向異性頁(yè)巖中[1]，地層微觀顆粒的排列、砂泥薄互層、不同顆?；虿煌紫斗植嫉貙雍土芽p地層都可能引起電性測(cè)量的各向異性[2-4]。傳統(tǒng)的隨鉆電磁波電阻率儀器采用軸向發(fā)射和接收天線，不具備方位探測(cè)特性，因此無(wú)法單獨(dú)確定地層走向方位和電阻率各向異性[5-8]，即使利用多條響應(yīng)曲線也很難消除電阻率各向異性的多解性，給含油飽和度的計(jì)算帶來(lái)了困難；同時(shí)，電阻率各向異性對(duì)地質(zhì)導(dǎo)向中界面距離的預(yù)測(cè)和判斷有比較大的影響，因此探測(cè)和識(shí)別電阻率各向異性變得越來(lái)越重要。

近年來(lái)，國(guó)外公司研發(fā)的隨鉆方位電磁波電阻率儀器提供了不同方向的電磁場(chǎng)分量信息，探測(cè)結(jié)果具有方位特性，可以確定地層的走向方位和各向異性。筆者從儀器測(cè)量信號(hào)出發(fā)，分析和模擬各向異性對(duì)隨鉆方位電磁波儀器信號(hào)響應(yīng)的影響，提出利用對(duì)稱補(bǔ)償測(cè)量方式來(lái)消除或增強(qiáng)電阻率各向異性的影響，并對(duì)隨鉆方位電磁波的電阻率資料進(jìn)行三參數(shù)反演，得到地層電阻率的各向異性，為更好地進(jìn)行資料解釋和地質(zhì)導(dǎo)向提供依據(jù)。

1 電阻率各向異性響應(yīng)特征模擬分析

Schlumberger公司2005年推出的第一代隨鉆方位電磁波電阻率測(cè)井儀PeriScope、Baker Hughes公司2006年推出的隨鉆方位電磁波電阻率測(cè)井儀APR、Halliburton公司2007年推出的隨鉆方位深探測(cè)電磁波電阻率測(cè)井儀ADR以及Weatherford公司2014年推出的GuidewaveTM，都增加了橫向或傾斜天線，具備多分量測(cè)量能力[9-10]。例如，方位電磁波電阻率測(cè)井儀APR[11]在傳統(tǒng)的電磁波電阻率儀器MPR軸向天線(T1—T4、R3和R4)基礎(chǔ)上增加了2個(gè)橫向天線R1和R2，如圖1所示。以儀器軸為z軸，儀器可以同時(shí)測(cè)量同軸電磁場(chǎng)分量Hzz和交叉耦合分量Hzx，儀器工作頻率為2 MHz和400 kHz，可以提供8條對(duì)稱補(bǔ)償電阻率曲線和16條不同方位補(bǔ)償定向電動(dòng)勢(shì)信號(hào)曲線。

圖1 Baker Hughes公司APR儀器線圈排布示意Fig.1 Coil configuration of the Baker Hughes APR tool

1.1 利用Hzz測(cè)量信號(hào)確定地層電阻率的各向異性

在電阻率均勻各向異性地層中，方位電磁波電阻率測(cè)井儀APR軸向接收天線的磁場(chǎng)信號(hào)Hzz可以表示為[12]：

(1)

(2)

(3)

式中：Hzz為zz分量磁場(chǎng)強(qiáng)度，A/m；L為線圈距，m；kh為水平波數(shù)；M為線圈磁矩，A·m2；θ為地層法線與儀器軸向的夾角，rad；β為各向異性因子；Rh為水平電阻率，Ω·m；Rv為垂直電阻率，Ω·m；λ為電阻率各向異性系數(shù)。

2個(gè)軸向接收天線的隨鉆電磁波電阻率幅度比和相位差轉(zhuǎn)換得到地層幅度電阻率和相位電阻率[13-14]，即儀器的視電阻率。由式(1)可知，軸向接收天線信號(hào)同時(shí)受到Rh和β的影響，而不受θ或電阻率各向異性系數(shù)單獨(dú)控制，不同的參數(shù)(水平電阻率、垂直電阻率和井眼地層相對(duì)夾角等)進(jìn)行組合可以得到相同的測(cè)量結(jié)果，因此采用傳統(tǒng)的隨鉆電磁波電阻率儀器測(cè)試數(shù)據(jù)確定地層的各向異性和相對(duì)夾角時(shí)存在很大的不確定性。

真實(shí)地層模型的水平電阻率和垂直電阻率分別為1和4 Ω·m，井眼與地層的相對(duì)夾角為60°時(shí)，儀器響應(yīng)值為Rta；水平電阻率從0.1 Ω·m至10.0 Ω·m，各向異性系數(shù)λ為4，井眼與地層相對(duì)夾角在0°～180°范圍內(nèi)變化時(shí)地層模型響應(yīng)值為Ra，將兩者進(jìn)行對(duì)比，模擬條件為APR儀器T1、R3和R4天線組合、工作頻率為2 MHz時(shí)的視電阻率，對(duì)比結(jié)果如圖2所示(圖中色柱表示不同參數(shù)組合時(shí)儀器響應(yīng)與真實(shí)地層模型響應(yīng)相對(duì)偏差的大小)。相對(duì)偏差的定義為：

(4)

式中：δ為相對(duì)偏差；Ra為井眼與地層相對(duì)夾角變化時(shí)模型響應(yīng)電阻率，Ω·m；Rta為地層真實(shí)模型響應(yīng)電阻率，Ω·m。

圖2 不同地層模型的電阻率響應(yīng)差別示意Fig.2 Resistivity response difference of different formation models

從圖2可以看出：儀器視電阻率在井眼和地層相對(duì)夾角為40°～140°范圍內(nèi)都存在可能的電阻率各向異性和相對(duì)夾角組合，使儀器響應(yīng)與真實(shí)地層模型的響應(yīng)相同，即用幅度電阻率或相位電阻率確定地層電阻率的各向異性存在很強(qiáng)的多解性。

模擬各向異性系數(shù)λ=4情況下，幅度電阻率和相位電阻率響應(yīng)隨井眼和地層相對(duì)夾角的變化情況，結(jié)果如圖3所示。

先前的研究已經(jīng)確定非洲豬瘟病毒不同分離株的最小或中立感染劑量在102和105之間，具體數(shù)量取決于毒株毒力、暴露頻率和入侵方式。目前還沒有中國(guó)流行的非洲豬瘟Georgia分離株在飼料中的感染劑量的報(bào)道。飼料或原科中經(jīng)典豬瘟病毒感染劑量的報(bào)道數(shù)據(jù)也很少。

圖3 各向異性地層模型的電阻率響應(yīng)模擬(λ=4)Fig.3 Resistivity response simulation of anisotropic formation model(λ=4)

從圖3可以看出，在0°～90°范圍內(nèi)，隨著相對(duì)夾角增大，幅度電阻率和相位電阻率單調(diào)遞增；相位電阻率受各向異性的影響比幅度電阻率要大，水平電阻率低時(shí)受到的各向異性影響相對(duì)較小。雖然同軸分量不能確定各向異性地層的走向和相對(duì)夾角，但在相對(duì)夾角大于40°時(shí)卻對(duì)電阻率各向異性有明顯的響應(yīng)，以此在某些確定地層條件下可以確定電阻率的各向異性。

1.2 Hzx測(cè)量信號(hào)確定地層電阻率的各向異性

在均勻各向異性地層中，方位電磁波電阻率測(cè)井儀APR橫向接收天線的磁場(chǎng)信號(hào)Hzx可以表示為[1]：

(5)

式中：Hzx為zx分量磁場(chǎng)強(qiáng)度，A/m。

與圖2采用的地層模型相同，單獨(dú)采用交叉耦合信號(hào)的實(shí)部或虛部很難確定各向異性以及井眼和地層的相對(duì)夾角，實(shí)部和虛部?jī)烧呓M合起來(lái)，對(duì)減少多解性有了較大改善，但仍存在一定的不確定性(見圖4)。當(dāng)電阻率與交叉耦合電動(dòng)勢(shì)信號(hào)組合，即Hzz和Hzx信號(hào)進(jìn)行交會(huì)(見圖5)，則基本能確定地層的各向異性以及井眼與地層的相對(duì)夾角，消除地層模型的不確定性。

圖4 不同地層模型的交叉耦合分量響應(yīng)差別示意Fig.4 Cross coupling component difference of different formation models

圖5 Hzz和Hzx分量組合響應(yīng)差別示意Fig.5 Response difference of Hzz and Hzx combinations

同樣地層條件下，模擬方位電磁波電阻率測(cè)井儀APR橫向接收天線電動(dòng)勢(shì)信號(hào)與相對(duì)夾角的關(guān)系，結(jié)果見圖6。

圖6 各向異性地層的交叉耦合電動(dòng)勢(shì)響應(yīng)(λ=4)Fig.6 Cross coupling electromotive response of anisotropic formation model(λ=4)

由圖6可知，交叉耦合電動(dòng)勢(shì)不再隨著井眼與地層相對(duì)夾角單調(diào)遞增；在夾角為0°(井眼垂直于地層)和90°(井眼平行于地層)時(shí)，儀器響應(yīng)接近于零；在0°～90°范圍內(nèi)存在一個(gè)極值，該極值對(duì)應(yīng)的相對(duì)夾角隨著各向異性系數(shù)的不同而不同。在相對(duì)夾角為90°時(shí)，橫向天線測(cè)量的交叉耦合分量為零，但磁場(chǎng)同軸分量Hzz各向異性反應(yīng)明顯，因此在相對(duì)夾角不太小的情況下可以利用Hzx結(jié)合Hzz確定地層的各向異性。

2 電阻率各向異性響應(yīng)補(bǔ)償及反演方法

2.1 電阻率各向異性響應(yīng)補(bǔ)償

由式(5)可知，均勻各向同性地層中沒有磁力線穿過(guò)橫向接收天線，因此Hzx為零。當(dāng)有地層界面存在時(shí)，由于界面對(duì)電磁場(chǎng)的反射，造成橫向接收天線產(chǎn)生電動(dòng)勢(shì)，幅度隨著儀器與界面之間距離的變化而變化，隨鉆方位電磁波儀器正是利用這一響應(yīng)特點(diǎn)進(jìn)行界面預(yù)測(cè)和判斷[15-16]。但由于地層存在各向異性，改變了這一響應(yīng)特點(diǎn)，即使在均質(zhì)地層中儀器的響應(yīng)也不再為零，給地質(zhì)導(dǎo)向帶來(lái)困難。在探測(cè)地層電阻率各向異性時(shí)，地層界面往往也會(huì)給測(cè)量帶來(lái)不利影響，通過(guò)對(duì)測(cè)量結(jié)果進(jìn)行補(bǔ)償處理，可以將電阻率各向異性和界面的影響進(jìn)行分離。

模擬層狀地層模型中間目的層存在各向異性，水平電阻率Rh=2.5 Ω·m，垂直電阻率Rv=10.0 Ω·m，層厚為6 m。儀器工作頻率為2 MHz，以不同的相對(duì)夾角穿過(guò)界面。忽略地層模型的層狀各向異性，可以應(yīng)用并矢格林函數(shù)得出任意方向磁偶極子的電磁場(chǎng)z分量(z為儀器軸線方向)，進(jìn)而推導(dǎo)出地層的電磁場(chǎng)分布[17]：

(6)

(7)

(8)

(9)

式中：Ez為電場(chǎng)強(qiáng)度z分量，V/m；Hz為磁場(chǎng)強(qiáng)度z分量，A/m；μ為磁導(dǎo)率，H/m；ω為角速度，rad/s；Mh，Mv分別為水平磁矩和垂直磁矩，A·m2；φ為方位角，rad；Jn(·)為n階Bessel函數(shù)；kp為積分變量；kzv,kzh分別為縱向波數(shù)和徑向波數(shù)；σv，σh分別為垂直電導(dǎo)率和水平電導(dǎo)率，S/m。

為了消除井眼偏心的影響，方位電磁波電阻率測(cè)井儀APR采用2個(gè)橫向接收天線信號(hào)的差值來(lái)預(yù)測(cè)和確定地層界面(發(fā)射天線T1、T2相對(duì)于2個(gè)橫向接收天線R1、R2對(duì)稱)，具體表示為[18-19]：

(10)

(11)

式中：HT1，HT2分別為T1和T2發(fā)射時(shí)補(bǔ)償后的磁場(chǎng)強(qiáng)度，A/m；H1，H2分別為橫向接收天線R1和R2接收到的磁場(chǎng)強(qiáng)度，A/m；d1，d2分別為T1與R1的距離以及R1與R2的距離(見圖1)，m。

根據(jù)以上方法得到的模擬結(jié)果如圖7—圖10所示。從模擬結(jié)果可以看出，T1或T2單獨(dú)發(fā)射時(shí)，除了相對(duì)夾角為0°和90°的情況外，在中間目的層中即使儀器遠(yuǎn)離地層界面，由于各向異性的影響使測(cè)量信號(hào)不再為零，給界面距離預(yù)測(cè)帶來(lái)一定困難；但值得注意的是，T1和T2分別發(fā)射時(shí)在地層界面的響應(yīng)類似，但各向異性帶來(lái)的影響正好相反(見圖7、圖8)。因此，將T1和T2分別發(fā)射時(shí)測(cè)得的信號(hào)相加，可以很大程度上消除各向異性的影響(見圖9)；同理，將T1和T2分別發(fā)射時(shí)測(cè)得的信號(hào)相減，可以減小界面的影響，從而增強(qiáng)各向異性的影響(見圖10)。利用儀器的這種響應(yīng)特點(diǎn)，可以對(duì)測(cè)量信號(hào)分別進(jìn)行處理，使測(cè)量資料更方便地分別用于地質(zhì)導(dǎo)向和地層電阻率各向異性的探測(cè)。

圖7 T1發(fā)射時(shí)交叉耦合電動(dòng)勢(shì)響應(yīng)模擬Fig.7 Cross coupling electromotive response simulation at T1 transmission

圖8 T2發(fā)射時(shí)交叉耦合電動(dòng)勢(shì)響應(yīng)模擬Fig.8 Cross coupling electromotive response simulation at T2 transmission

圖9 界面響應(yīng)增強(qiáng)型補(bǔ)償Fig.9 Enhanced compensation for interface response

圖10 各向異性增強(qiáng)型補(bǔ)償Fig.10 Enhanced compensation for anisotropy

2.2 電阻率各向異性三參數(shù)反演

通過(guò)補(bǔ)償測(cè)量，可以得知地層電阻率是各向同性還是各向異性，從而在反演過(guò)程中建立更加準(zhǔn)確的地層模型。在隨鉆測(cè)量時(shí)利用實(shí)時(shí)上傳的2條幅度電阻率、相位電阻率曲線和一條定向電動(dòng)勢(shì)曲線進(jìn)行三參數(shù)反演，常用2層或3層地層模型。將測(cè)井?dāng)?shù)據(jù)表示為向量形式，地層參數(shù)與儀器響應(yīng)可以表示為：

(12)

(13)

式中：Rh為水平電阻率；Rv為垂直電阻率；y為隨鉆方位儀器響應(yīng)值；f為儀器響應(yīng)函數(shù)；α為井斜角；Φ為目標(biāo)函數(shù)。

反演的過(guò)程就是利用最小二乘法逐次修正Rh、Rv和α，最后實(shí)現(xiàn)Φ最小化。

利用圖7—圖10的定向電動(dòng)勢(shì)響應(yīng)曲線以及電阻率響應(yīng)曲線進(jìn)行三參數(shù)反演，得到水平電阻率和垂直電阻率(見圖11)，以及相對(duì)井斜角與真實(shí)井斜角的相對(duì)關(guān)系(見圖12)。

圖11 各向異性反演結(jié)果Fig.11 Inversion results of anisotropy

圖12 相對(duì)井斜角反演結(jié)果Fig.12 Inversion results of relative deviation angle

由反演結(jié)果可知，相對(duì)井斜角為零時(shí)，定向電動(dòng)勢(shì)響應(yīng)和電阻率響應(yīng)均不受電阻率各向異性影響；相對(duì)井斜角不為零時(shí)，電阻率各向異性和相對(duì)井斜角反演結(jié)果與地層模型一致。因此，對(duì)隨鉆方位電磁波電阻率測(cè)量資料進(jìn)行反演，可以確定地層電阻率的各向異性。

3 結(jié) 論

1) 傳統(tǒng)隨鉆電磁波電阻率的同軸耦合電動(dòng)勢(shì)與交叉分量耦合電動(dòng)勢(shì)對(duì)地層電阻率各向異性的響應(yīng)規(guī)律不同，隨鉆方位電磁波電阻率測(cè)井儀APR利用幅度電阻率和相位電阻率，并結(jié)合橫向接收天線測(cè)得的電動(dòng)勢(shì)的實(shí)部和虛部信息進(jìn)行交會(huì)，解決了地層電阻率各向異性多解性的問(wèn)題。

2) 利用橫向接收天線信號(hào)進(jìn)行對(duì)稱發(fā)射補(bǔ)償測(cè)量時(shí)，對(duì)不同發(fā)射天線發(fā)射時(shí)測(cè)量得到的交叉耦合信號(hào)進(jìn)行相加，可以降低各向異性的影響，增強(qiáng)界面響應(yīng)；對(duì)不同發(fā)射天線發(fā)射時(shí)測(cè)量得到的交叉耦合信號(hào)進(jìn)行相減，則會(huì)增強(qiáng)各向異性的影響，降低界面響應(yīng)。

3) 對(duì)方位電磁波電阻率測(cè)量得到的定向電動(dòng)勢(shì)和電阻率曲線進(jìn)行三參數(shù)反演，得到地層水平電阻率、垂直電阻率和相對(duì)井斜角，從而準(zhǔn)確地確定地層電阻率的各向異性。

[1] FANG S，MERCHANT G A，HART E，et al.Determination of structural dip and azimuthal from LWD azimuthal propagation resistivity measurements in anisotropic formations[R].SPE 116123，2008.

[2] YIN Hezhu，KURNIAWAN B.Resistivity anisotropy models and multi-component induction measurements：impact on Swand uncertainty 0 HPV estimation[R].SPWLA-2008-LLLL，2008.

[3] ANDERSON B I，BARBER T D，GIANZERO S C.The effect of crossbedding anisotropy on induction tool response[R].SPWLA-2001-v42n2a6，1998.

[4] 閆建平，蔡進(jìn)功，趙銘海，等.考慮巖石結(jié)構(gòu)信息的砂礫巖體有效儲(chǔ)層測(cè)井判識(shí)方法[J].同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版)，2011，39(9)：1365-1372.

YAN Jianping，CAI Jingong，ZHAO Minghai，et al.Identification method of effective reservoir for glutenite body using well logging based on rock texture[J].Journal of Tongji University(Nature Science)，2011，39(9)：1365-1372.

[5] 高杰，陳木銀，陳雅薇，等.定向井各向異性地層交流電測(cè)井響應(yīng)模擬[J].勘探地球物理進(jìn)展，2003，26(4)：58-61.

GAO Jie，CHEN Muyin，CHEN Yawei，et al.Forward modeling of AC logging response in directional wells of anisotropic formations[J].Progress in Exploration Geophysics，2003，26(4)：58-61.

[6] 楊震，楊錦舟，韓來(lái)聚.隨鉆方位電磁波電阻率成像模擬及應(yīng)用[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版)，2013，43(6)，342-350.

YANG Zhen，YANG Jinzhou，HAN Laiju.Numerical simulation and application of azimuthal propagation resistivity imaging while drilling[J].Journal of Jilin University(Earth Science Edition)，2013，43(6)：342-350.

[7] 王昌學(xué)，周燦燦，儲(chǔ)昭坦，等.電性各向異性地層頻率域電磁響應(yīng)模擬[J].地球物理學(xué)報(bào)，2006，49(6)：1873-1883.

WANG Changxue，ZHOU Cancan，CHU Zhaotan，et al.Modeling of electromagnetic responses in frequency domain to electrical anisotropic formations[J].Chinese Journal of Geophysics，2006，49(6)：1873-1883.

[8] LI Shanjun，CHEN Jiefu，BINFORD T L Jr.Using new LWD measurements to evaluate formation resistivity anisotropy at any dip angle wells[R].SPWLA-2014-EEEE，2014.

[9] LI Qiming，OMERAGIC D，CHOU L，et al.New directional electromagnetic tool for proactive geosteering and accurate formation evaluation while drilling[R].SPWLA-2005-UU,2005.

[10] BITTAR M S，KLEIN J D，RANDY B，et al.A new azimuthal deep reading resistivity tool for geosteering and advanced formation evaluation[R].SPE 109971,2007.

[11] ZHANG Zhiyi，GONGUET C，RAJANI V，et al.Directional LWD resistivity tools and their business impacts[R].SPWLA-2008-FFFF，2008.

[12] HAGIWARA T.A new method to determine horizontal-resistivity in anisotropic formations without prior knowledge of relative dip[R].SPWLA-1996-Q，1996.

[13] 康俊佐，邢光龍，楊善德.多傳播電阻率測(cè)井的探測(cè)能力與資料處理方法[J].石油勘探與開發(fā)，2006，33(3)：97-102.KANG Junzuo，XING Guanglong，YANG Shande.Detective ability and data processing method of the MPR logging[J].Petroleum Exploration and Development，2006，33(3)：97-102.

[14] 孫向陽(yáng)，聶在平，趙延文，等.用矢量有限元方法模擬隨鉆測(cè)井儀在傾斜各向異性地層中的電磁響應(yīng)[J].地球物理學(xué)報(bào)，2008，51(5)：326-333.

SUN Xiangyang，NIE Zaiping，ZHAO Yanwen，et al.The electromagnetic modeling of logging-while-drilling tool in tilted anisotropic formations using vector finite element method[J].Chinese Journal of Geophysics，2008，51(5)：326-333.

[15] OMERAGIC D，LI Qiming，CHOU L，et al.Deep directional electromagnetic measurement for optimal placement[R].SPE 97045，2005.

[16] CALLEJA B，MARKET J，PITCHER J，et al.Multi-sensor geosteering[R].SPWLA-2010-82670，2010.

[17] 肖加奇，張庚驥.水平井和斜度井中的感應(yīng)測(cè)井響應(yīng)計(jì)算[J].地球物理學(xué)報(bào)，1995，38(3)：396-404.

XIAO Jiaqi，ZHANG Gengji.Computation of induction logging response in horizontal and highly-deviated wells[J].Chinese Journal of Geophysics，1995，38(3)：396-404.

[18] CHEMALI R E，CAIRNS P，WANG T，et al.Method for signal enhancement in azimuthal propagation resistivity while drilling：7375530B2[P].2008-05-20.

[19] WANG T，CHEMALI R E.Method of generating a deep resistivity image in LWD measurement：7483793B2[P].2009-01-27.

[編輯 滕春鳴]

A New Method of Detecting while Drilling Resistivity Anisotropy with Azimuthal Electromagnetic Wave Tools

YANG Zhen1, WEN Yi2, XIAO Hongbing1

(1.Drilling Technology Research Institute, Sinopec Shengli Oilfield Service Corporation, Dongying，Shandong，257017，China；2. Well Logging Company, Sinopec Shengli Oilfield Service Corporation, Dongying，Shandong，257096, China)

With the extensive application while drilling of azimuthal electromagnetic wave tools in highly deviated wells or horizontal wells, resistivity anisotropy is currently one of the main factors influencing geosteering and formation evaluation accuracy. Different signals may interface from the azimuthal electromagnetic wave tools. The study used numerical simulation on Baker Hughes’s APR tool to determine the effect of resistivity anisotropy on the response of each component signal. In the simulation, the effect of anisotropy was enhanced or removed by using the symmetrically transmitted compensation measurement. The anisotropy of formation resistivity was calculated from forward and inversion. It showed in the simulation that when relative deviation angle was in the range of 0°-90°, axial component of magnetic field signals increased monotonically with resistivity anisotropy and both followed the parabolic law. When the relative deviation angle was 0° or 90°, the effect would diminish. The anisotropy of formation resistivity could be confirmed effectively by the combination of different signal components of electromagnetic resistivity tools while drilling. Horizontal resistivity, vertical resistivity and relative deviation angle might be obtained from three-parameter inversion on azimuthal electromagnetic wave response data. Formation evaluation and geosteering may be more accurate when based on the inverted measurement data.

azimuthal electromagnetic wave; resistivity anisotropy; compensation; inversion; geosteering; formation evaluation

2015-07-09；改回日期：2016-02-16。

楊震(1982— )，男，山東萊蕪人，2004年畢業(yè)于石油大學(xué)(華東)勘查技術(shù)與工程專業(yè)，2009年獲中國(guó)石油大學(xué)(華東)地質(zhì)資源與地質(zhì)工程專業(yè)博士學(xué)位，高級(jí)工程師，主要從事隨鉆電測(cè)井方法及儀器設(shè)計(jì)方面的研究。E-mail：stingsyzh@126.com。

國(guó)家科技重大專項(xiàng)“低滲透油氣深層高溫高壓隨鉆測(cè)控技術(shù)”(編號(hào)：2016ZX05021-001)資助。

?測(cè)井錄井?

10.11911/syztjs.201603021

P631.8+13

A

1001-0890(2016)03-0115-06