王瑞甲， 喬文孝*

1 中國石油大學(xué)油氣資源與探測國家重點實驗室, 北京 102249 2 北京市地球探測與信息技術(shù)重點實驗室， 北京 102249

?

VTI地層隨鉆四極子聲波測井?dāng)?shù)值研究

王瑞甲1,2， 喬文孝1,2*

1 中國石油大學(xué)油氣資源與探測國家重點實驗室, 北京 102249 2 北京市地球探測與信息技術(shù)重點實驗室， 北京 102249

數(shù)值計算了四極子聲源在含鉆鋌豎直向橫向各向同性(VTI)地層井孔中產(chǎn)生的各種模式波的頻散曲線、激發(fā)幅度以及對各地層、鉆鋌參數(shù)的靈敏度，合成了四極子陣列波形，研究了快、慢速地層井孔中地層各向異性的存在對四極子波場的影響.數(shù)值計算結(jié)果表明，地層各向異性對鉆鋌波的影響很小，對于地層波的影響較大，且影響因素非常復(fù)雜.僅在一些特殊的頻率點處，地層模式波的控制因素相對較為簡單，如螺旋波的截止頻率處.在慢速地層井孔中，最低階的地層四極子波對距離井壁約2個井孔半徑內(nèi)的地層有明顯的響應(yīng)，可以對該范圍內(nèi)的地層橫波進行層析成像.在快速地層井孔條件下，可以考慮采用與橫波到時相同的波包來評價地層橫波信息.該波包包含了井孔折射橫波、F1和F2模式與井孔折射橫波速度相接近的部分；利用該波包獲得的地層橫波速度基本不受地層各向異性的影響.在慢速地層情況下，螺旋波受地層各向異性的影響較大，建議在實際數(shù)據(jù)處理時，考慮采用基于數(shù)據(jù)的處理方法進行頻散校正.

四極子； 隨鉆測井； 聲波測井； 橫波

1 引言

隨鉆聲波測井可以在鉆進過程中提供實時的地層聲學(xué)信息，在節(jié)約井口占用時間、規(guī)避風(fēng)險等方面有著傳統(tǒng)電纜測井不可比擬的優(yōu)勢(Tang and Cheng, 2004).早期的隨鉆聲波儀器設(shè)計，照搬電纜測井已有的技術(shù)，未對含鉆鋌情況下的井孔波導(dǎo)做細(xì)致的分析，未能開發(fā)出成功的橫波測量儀器.后來，人們逐步認(rèn)識到了含鉆鋌井孔內(nèi)鉆鋌波與地層波的耦合作用，認(rèn)為四極子聲源激發(fā)的螺旋波更適合隨鉆橫波測量，并發(fā)展了采用四極子聲源和接收器的隨鉆橫波測井技術(shù)，實現(xiàn)了任意地層的隨鉆橫波速度評價(Tang et al., 2002;Tang et al., 2003; Kinoshita et al.,2008;Zhang et al., 2010;Yang et al., 2011).在國內(nèi)，目前隨鉆聲波儀器的研發(fā)還處于起步階段，迫切需要充分的理論研究和方法論證(肖紅兵等, 2009; 盧俊強等, 2013).前人已經(jīng)對隨鉆條件下多極子聲源在井孔內(nèi)激發(fā)的聲場做了大量的研究，主要集中在井孔導(dǎo)波的模式分析(Tang et al., 2002;Tang et al., 2003; Wang and Tang, 2003; 王華等, 2009; Sinha et al., 2009；Sinha and Simsek, 2010;Geerits et al，2010; Wang and Tao, 2011; 李希強等,2013)，隔聲體設(shè)計(楊勇等,2009;蘇遠(yuǎn)大等,2011)以及物理模擬研究(Zhu et al., 2008;王兵等, 2012)等諸多方面.實際上，隨鉆四極子測井受到諸多因素的影響，目前的四極子隨鉆聲波測井理論還不完備，在一些基本問題上還存在認(rèn)識不清.這包括地層各向異性對隨鉆四極子波場中各模式的影響，快、慢速地層下隨鉆四極子波場的不同以及螺旋波的徑向探測深度等問題.對這些問題的研究和分析，有助于完善隨鉆四極子橫波測井的理論，為新一代的隨鉆聲波測井儀器的設(shè)計提供必要的理論基礎(chǔ).本文建立了豎直向橫向各向同性(VTI)地層豎直井隨鉆測井聲學(xué)模型，數(shù)值計算了快、慢TI地層情況下四極子聲源激勵的聲場，并嘗試對上述問題進行了研究和解釋.

2 VTI地層豎直井隨鉆測井聲學(xué)模型

隨鉆測井聲學(xué)模型可以用柱狀徑向分層聲學(xué)模型來描述.圖1為隨鉆測井聲學(xué)模型示意圖(Tang and Cheng, 2004;崔志文, 2004)，包括(a)模型的三維示意圖和(b)井孔橫截面示意圖，S表示TI介質(zhì)的對稱軸，井軸與介質(zhì)對稱軸S平行.模型由內(nèi)向外各介質(zhì)依次為水眼(水)、鉆鋌、井孔(水)和地層，各介質(zhì)的外徑分別為r0、r1、r2和無窮大.井孔外地層為無限大TI介質(zhì).鉆鋌位于井孔中央，鉆鋌中間的水眼中充滿水.特別地，在本文的討論中，規(guī)定“井孔”為鉆鋌和地層之間的空間，并假設(shè)鉆鋌在井孔中

圖1 VTI地層豎直井隨鉆測井聲學(xué)模型示意圖(a)及鉆鋌橫截面示意圖(b)Fig.1 Schematic diagram of acoustic model for logging while drilling in a VTI formation (a) and the cross section of a drill collar (b)

完全居中.

3 理論

3.1 井孔內(nèi)外聲場的解

在計算中，我們假設(shè)聲源為安裝在鉆鋌外壁的理想環(huán)狀多極子聲源.另外，假設(shè)水眼內(nèi)代表沿徑向向內(nèi)傳播的縱波的勢函數(shù)的系數(shù)為Awh，鉆鋌內(nèi)代表向內(nèi)和向外傳播的縱波和兩種橫波的勢函數(shù)的系數(shù)分別為Aco、Bco、Cco、Dco、Eco和Fco，井孔中表示向內(nèi)和向外傳播的縱波的勢函數(shù)的系數(shù)分別為Abh、Bbh，地層中表示向內(nèi)和向外傳播的縱波和兩種橫波的勢函數(shù)的系數(shù)分別為Afo、Bfo、Cfo、Dfo、Efo和Ffo，縮寫“co”、“wh”、“bh”、“fo”分別表示鉆鋌、水眼、井孔和地層.在柱坐標(biāo)系(r,θ,z)下，根據(jù)水眼-鉆鋌、鉆鋌-井孔和井孔-地層處的邊界條件(Tang and Cheng, 2004;崔志文, 2004)，得到

(1)

(2)

(3)

式(1)—(3)分別代表水眼-鉆鋌、鉆鋌-井孔、井孔-地層邊界，其中u表示徑向位移，p表示聲壓，τ表示應(yīng)力.聯(lián)立式(1)—(3)得到

M12×12X=S，

(4)

其中，

(5)

其中，f為井孔內(nèi)流體的徑向波數(shù)，Kn和In表示修正的貝塞爾函數(shù)，εn為多極子聲源的權(quán)重因子，當(dāng)n=0時，εn=1，當(dāng)n>0時，εn=2.有關(guān)模型內(nèi)各物理量及矩陣M的詳細(xì)表達式，可參見文獻(王瑞甲, 2012;李希強等, 2013).

通過對式的求解，可得到描述井孔內(nèi)外聲場的勢函數(shù)的系數(shù)：

X=M-1S.

(6)

進一步，可以求得井內(nèi)聲場.與電纜測井的情況不同，在隨鉆測井中，采用兩個勢函數(shù)的系數(shù)來描述井孔聲場.例如儀器外環(huán)上的聲壓場P如下(Tangetal., 2003)：

(7)

其中，ρf表示井孔流體密度，ω表示頻率，S表示聲源的子波函數(shù).

3.2 模式波速度及激發(fā)強度

各模式波的速度、衰減由勢函數(shù)系數(shù)的極點決定(Tang and Cheng, 2002)，即

(8)

式(8)又稱作頻散方程或者周期方程.該方程可以采用牛頓迭代法等數(shù)值方法求解，得到各模式波對應(yīng)的軸向波數(shù)的極點ki，并進一步可得到模式波的傳播速度vi=ω/ki.模式波的激發(fā)函數(shù)由極點對應(yīng)的留數(shù)給出，根據(jù)留數(shù)定理可得

(9)

在TI介質(zhì)井孔聲場研究中，靈敏度分析是研究模式波速度受各地層參數(shù)控制程度的有效手段.對于隨鉆聲波測井而言，靈敏度可以幫助分析模式波的類型及主要影響因素.模式波對于某參數(shù)P的靈敏度定義為歸一化后的相速度對該參數(shù)的偏微分(Tang and Cheng, 2004)：

(10)

模式波速度對于徑向上某位置上的某參數(shù)的靈敏度定義為(陳雪蓮和王瑞甲, 2008)

(11)

(12)

其中c11、c33、c44和c66表示地層剛性系數(shù)，ρ為地層的密度.

4 數(shù)值計算結(jié)果及分析

4.1 數(shù)值模型參數(shù)

計算中采用的TI地層的剛性參數(shù)如表1所示.

表1 地層參數(shù)Table 1 Formation parameters

表2 本文中井孔流體及鉆鋌的參數(shù)Table 2 Parameters of the borehole fluid and collar

本文中的快速、中速和慢速地層分別指其兩種橫波速度(vSH,vSV)均大于、接近或者小于井孔流體聲速的地層.井孔內(nèi)流體及鉆鋌參數(shù)見表2.計算采用的鉆鋌內(nèi)半徑、外半徑及井眼半徑分別為0.027 m、0.090 m和0.120 m.本文同時將TI地層計算結(jié)果與對應(yīng)的各向同性地層計算結(jié)果進行了對比，用于參考的各向同性地層的聲速選擇為TI地層中沿旋轉(zhuǎn)對稱軸方向傳播的縱、橫波的速度，密度與TI地層相同.

4.2 快速地層

圖2顯示了快速地層情況下，四極子波的頻散曲線和激發(fā)強度曲線.在計算的頻率范圍內(nèi)觀察到了四種模式，分別是地層螺旋波(F1、F2)、鉆鋌波(C1、C2).參照Sinha 等 (2009)，采用字符“C”和“F”的分別表示鉆鋌模式和地層模式，采用標(biāo)識符“ISO”表示對應(yīng)的各向同性地層.另外，在下面的頻散曲線圖中，實線表示相速度，虛線表示群速度.

從圖2可見，地層四極子波(F1、F2)在截止頻率處均趨近于地層沿豎直方向傳播的橫波速度vSV.螺旋波(F1)是強頻散的模式，在截止頻率處觀察不到各向異性造成的速度的差異，而在稍高的頻率下，可見各向異性地層螺旋波速度大于對應(yīng)的各向同性地層的螺旋波.F2模式頻散較弱.各向異性地層的F2模式截止頻率小于各向同性地層；相同頻率下，各向異性地層的F2模式的速度大于各向同性地層.未觀察到地層各向異性對鉆鋌模式(C1、C2)的影響.從圖2b可見，地層各向異性對地層波(F1、F2)的激發(fā)強度有輕微影響：F1模式各向異性地層和各向同性地層下的激發(fā)強度基本一致；對于F2模式，在17～27 kHz頻率范圍內(nèi)，各向同性地層下的激發(fā)強度略大于各向異性地層的.地層各向異性對于鉆鋌模式(C1、C2)的激發(fā)強度基本不存在影響.

圖3(a—d)分別是圖2中F1、F2、C1和C2模式的靈敏度分析結(jié)果.從圖3可見，鉆鋌模式波(C1、C2)基本上由鉆鋌縱波(vCP)、鉆鋌橫波速度(vCS)、井孔流體速度(vBF)和鉆鋌內(nèi)流體速度(vCF)控制，地層參數(shù)對其影響較小.地層螺旋波(F1)在截止頻率處由沿豎直方向傳播的橫波速度(vSV)控制，而后隨著頻率增加，F(xiàn)1對于vSV的靈敏度下降很快，對其他地層參數(shù)和井孔流體速度的靈敏度先增加，后減小.高頻處，F(xiàn)1模式主要為井孔流體速度(vBF)控制.F1模式對鉆鋌縱、橫波速度不敏感，靈敏度小于0.1.F2模式在截止頻率處也主要由vSV控制，隨頻率增加，vSV的控制緩慢地減弱，其他參數(shù)對F2的影響增加.一些學(xué)者提出了在快速地層中，采用F2模式而非F1模式來提取沿井軸方向傳播的橫波速度(Tang and Cheng, 2004)，主要基于F2模式頻散較弱，而F1模式是強頻散波的理論.從本文的結(jié)果可見，因為在F2模式的截止頻率，除vSV的其他參數(shù)均對F2模式的速度影響較小(靈敏度均小于0.2)，采用F2模式波可以減小地層各向異性對于提取結(jié)果的影響.崔志文(2004)的理論模擬結(jié)果表明，在他的模型中，F(xiàn)2模式波的群速度和鉆鋌波的群速度相接近，鉆鋌波可能嚴(yán)重干擾F2模式的波形.實際上，鉆鋌波頻散特征與鉆鋌尺寸有關(guān)：當(dāng)鉆鋌外徑不變，內(nèi)徑減小時，頻散曲線向高頻偏移，反之則向低頻偏移(Tang et al., 2003; 崔志文, 2004).在實際儀器設(shè)計中，可選擇合適的儀器結(jié)構(gòu)盡量避開鉆鋌波的干擾.

圖2 隨鉆情況下快速VTI地層四極子聲場模式分析，包括(a)頻散曲線和(b)激發(fā)強度Fig.2 Mode analysis of quadrupole modes in a fast VTI formation, including (a) dispersion curves and (b) excitation intensity

圖3 隨鉆情況下快速VTI地層四極子聲場靈敏度分析，包括(a) F1、(b) F2、(c) C1、(d) C2Fig.3 Sensitivity analysis of quadrupole modes in a fast VTI formation, including (a) F1, (b) F2, (c) C1, (d) C2

4.3 中速地層

圖4顯示，在中速地層計算結(jié)果中觀察到了地層螺旋波(F1)和鉆鋌波(C1)兩種模式，而且鉆鋌波的頻率遠(yuǎn)高于地層螺旋波激發(fā)的頻率，地層螺旋波在截止頻率處的速度為沿豎直方向傳播的橫波速度vSV.從圖4a可見，各向異性地層螺旋波速度高于各向同性地層螺旋波速度，未觀察到地層各向異性對C1模式波速度的影響.相比于各向同性地層(F1-ISO)，各向異性地層螺旋波(F1)激發(fā)強度的主頻向高頻偏移，激發(fā)強度略降低.

圖5顯示了圖4中各模式波的靈敏度分析結(jié)果.螺旋波(F1)在低頻處主要由沿豎直方向傳播的橫波速度vSV控制，而在稍高的頻率內(nèi)vPH、vSH和c13等對螺旋波速度都有一定的影響，說明此時地層各向異性對螺旋波有一定的影響，與圖4一致.另外，地層螺旋波對鉆鋌縱、橫波速度及鉆鋌內(nèi)流體的速度的靈敏度基本為零，說明鉆鋌的彈性參數(shù)對地層螺旋波基本沒有影響.從圖5b可見，鉆鋌波(C1)基本由鉆鋌縱、橫波速度控制，其他參數(shù)對其影響很小.4.4 慢速地層

圖6是慢速地層四極子聲源激發(fā)的各種模式波的頻散曲線和激發(fā)強度曲線.在6a中，可觀察到地層螺旋波(F1)、泄漏模式縱波(F2)和鉆鋌波(C1)三種模式波.地層模式波(F1、F2)受地層各向異性的影響較大：各向異性地層螺旋波速度大于各向同性地層，而泄漏模式縱波的速度小于各向同性地層.與中速地層一致，地層螺旋波在截止頻率處的速度為沿豎直方向傳播的橫波速度vSV.鉆鋌波(C1)存在的頻率高于地層螺旋波，未觀察到地層各向異性對于鉆鋌波的影響.各向異性對于地層螺旋波的激發(fā)幅度有一定的影響，一方面各向異性地層螺旋波的激發(fā)的主頻向高頻偏移，另一方面激發(fā)強度較各向同性地層降低.

圖4 隨鉆情況下中速VTI地層四極子聲場模式分析，包括(a)頻散曲線和(b)激發(fā)強度Fig.4 Mode analysis of quadrupole modes in a middle VTI formation, including (a) dispersion curves and (b) excitation intensity

圖5 隨鉆情況下中速VTI地層四極子聲場靈敏度分析，包括(a)F1、(b)C1Fig.5 Sensitivity analysis of quadrupole modes in a middle VTI formation, including (a) F1, (b) C1

圖6 隨鉆情況下慢速VTI地層四極子聲場模式分析，包括(a)頻散曲線和(b)激發(fā)強度Fig.6 Mode analysis of quadrupole modes in a slow VTI formation, including (a) dispersion curves and (b) excitation intensity

圖7 隨鉆情況下慢速VTI地層四極子聲場靈敏度分析，包括(a) F1、(b) F2、(c) C1Fig.7 Sensitivity analysis of quadrupole modes in a slow VTI formation, including (a) F1, (b) F2, (c) C1

圖7顯示了慢速地層中各模式波的靈敏度分析結(jié)果.由圖7可見，低頻螺旋波(F1)的速度由沿豎直方向傳播的橫波速度(vSV)控制，而在稍高的頻率下，地層各參數(shù)及井孔流體參數(shù)均對地層螺旋波的速度有一定的影響.另外，地層螺旋波對于鉆鋌縱、橫波速度及鉆鋌內(nèi)流體速度不敏感.F2模式主要由沿豎直方向傳播的縱波速度vPV和剛性系數(shù)c13控制，vSV、vSH、vBF以及vPH對其速度也有一定的影響.鉆鋌模式波(C1)主要由鉆鋌縱、橫波速度控制，在9～10 kHz頻率內(nèi)沿水平方向傳播的縱、橫波速度對其有一些影響.

4.5 地層螺旋波的徑向探測特性

圖8(a—c)分別為快速地層、中速地層和慢速地層井孔中，地層螺旋波(F1)對于沿豎直方向傳播的橫波速度vSV的靈敏度.根據(jù)物理學(xué)的相似性原理，對于地層波而言，本征頻率與井孔尺寸基本呈倒數(shù)關(guān)系，即縮小井半徑，地層模式波的頻散曲線、激發(fā)強度及靈敏度曲線向高頻拉伸，反之則向低頻壓縮.因此，本文采用了井孔半徑個數(shù)來表示徑向靈敏度及其探測深度，以便于把結(jié)論推廣到不同尺寸的井孔模型中.從圖8(a—c)可見，距離井壁越遠(yuǎn)，地層螺旋波對于vSV的靈敏度越低，在截止頻率附近，螺旋波對距離井壁大約2個井孔半徑以內(nèi)的地層有響應(yīng).相對于低頻情況，在高頻情況下，靈敏度隨著徑向距離的增加下降速度更快，這說明螺旋波基本呈現(xiàn)低頻探測深，高頻探測淺的特征.在4 kHz左右，快速地層中地層螺旋波對于vSV的靈敏度較高，而隨頻率增加很快降低，這是因為高頻下的快速地層螺旋波的速度主要由井內(nèi)流體速度控制.

可以根據(jù)徑向靈敏度的計算結(jié)果來定量評價模式波的徑向探測點(Investigation Position, IP).本文定義能夠使得該徑向位置之外的地層和該徑向位置之內(nèi)的地層對模式波的貢獻相等的位置為模式波的徑向探測點:

(13)

其中，rIP代表徑向探測點距離井壁的距離，r為任意一點距離井壁的距離，Sen(r,ω)為徑向靈敏度.圖8d提取了各種地層情況下地層螺旋波(F1)的徑向探測點.從圖8d可見，中速地層低頻螺旋波最深的徑向探測點離井壁的距離約為0.8個井孔半徑，而慢速地層為0.45個井孔半徑，在8 kHz下兩種地層的螺旋波的徑向探測點距井壁的距離僅為0.1個井孔半徑左右.

4.6 合成的陣列波形

利用式(7)合成井孔內(nèi)接收的四極子波形，圖9顯示了不同頻率下四極子聲源在井孔內(nèi)激發(fā)的聲場，其中聲源的子波函數(shù)為瑞克子波，接收器源距為2.0 m，S標(biāo)識的直線表示橫波的到時，(a)和(b)分別表示快、慢速各向同性地層.由圖9a可見，觀察到了多個波包，分別以符號W1、W2和W3標(biāo)識.最先到達的為W1波包，到時與地層橫波相同.結(jié)合圖2結(jié)果可知，W1包含多種模式：井孔折射橫波，F(xiàn)1模式截止頻率處的部分和F2模式截止頻率處的部分(10～12 kHz).由于這幾種模式的群速度上幾乎相同，在時域內(nèi)重疊在一起，不能夠分開.W2和W3波包在W1之后到達，主要為F1和F2模式中速度與地層橫波速度差異較大的部分.綜合圖2和圖9可知，W1波包能夠較為準(zhǔn)確地反應(yīng)地層橫波速度，無需進行頻散校正等后續(xù)處理，是快速地層橫波速度測量的首選.對比不同頻率下的波形，可見，聲源主頻在8～12 kHz時，W1波包幅度較大，且受后續(xù)波影響較小，說明對于該地層，儀器合適的工作主頻在8～12 kHz之間.另外，從圖2也可以看到，在該頻率范圍內(nèi)，鉆鋌波的群速度較小，不會影響W1波包.由圖9b可見，在聲源主頻為2 kHz時，螺旋波(F1)波至點與地層橫波相同，隨著頻率增加，螺旋波波至點逐漸后移.對于該模型而言，采用2 kHz的聲源進行測量，可測得與地層橫波速度最接近的螺旋波信號.在實際應(yīng)用中，應(yīng)綜合考慮噪聲、換能器特點等各種因素，選擇合適的頻率.

圖8 隨鉆情況下(a)快速地層、(b)中速地層、(c)慢速地層中螺旋波的徑向靈敏度及(d)徑向探測點Fig.8 Radial sensitivity analysis of screw mode in a (a) fast, (b)middle, (c) slow formation and (d) the investigation depth of Stoneley mode

圖9 不同頻率下，四極子聲源在快(a)、慢(b)地層井孔中產(chǎn)生的聲場Fig.9 Acoustic wave field in fluid-filled borehole produced by a quadrupole source with different frequency

圖10 四極子聲源在各向異性地層與各向同性井孔中產(chǎn)生的聲場的對比，包括(a)快速地層和(b)慢速地層Fig.10 Comparison of the waveforms between TI formations and corresponding isotropic formations produced by a quadrupole source, including (a) fast formation and (b) slow formation

圖10對比了TI地層和各向同性地層情況下四極子聲源在含鉆鋌井孔中產(chǎn)生的聲場，其中實線表示TI地層計算結(jié)果，虛線表示對應(yīng)的各向同性地層的計算結(jié)果.快、慢速地層采用的四極子源的主頻分別為8和2 kHz.圖10a顯示，W1波包到時為地層橫波到時，基本上不受地層各向異性的影響，能夠準(zhǔn)確反映地層的橫波速度；后續(xù)波包受到了地層各向異性的影響，兩種情況下的波形的相位有差別，幅度基本一致.從圖10b可見，在慢速地層情況下，地層螺旋波的幅度和相位均受地層各向異性的影響較大.各向異性的存在使得F1模式到時提前，幅度降低，與圖6中模式分析的結(jié)果相一致.

5 討論

與偶極子聲場不同，從四極子聲場中各模式波的靈敏度分析結(jié)果可見，地層波基本不受鉆鋌彈性參數(shù)影響，鉆鋌波基本不受地層參數(shù)影響，未觀察到地層波和鉆鋌波相互耦合的現(xiàn)象，與前人的結(jié)論一致(崔志文，2004; Sinha et al., ,2009;Geerits et al., 2010).地層各向異性對鉆鋌波的影響很小，對于地層波的影響較大，且影響因素非常復(fù)雜.對于本文研究的三種地層模型，僅在一些特殊的頻率點處，地層模式波的控制因素相對較為簡單，如螺旋波的截止頻率.不過，雖然從這些頻率點可以提取出不受地層各向異性影響的橫波速度信息，但由于該頻率下的模式波的激發(fā)幅度較小，波形信噪比較差，提取結(jié)果可能存在較大的誤差.

螺旋波有著與電纜測井中的彎曲波(陳雪蓮和王瑞甲, 2008)類似的徑向探測特征：低頻探測深、高頻探測淺.與Su 等(2013)的結(jié)果不同，本文計算結(jié)果顯示，對于本文的快速地層，由于稍高頻率下螺旋波(F1)主要為井孔內(nèi)流體的縱波速度控制，對地層橫波速度不敏感，利用F1模式進行橫波徑向?qū)游龀上癫⒉皇且粋€好方法.對于慢速地層，由靈敏度分析結(jié)果可以看到，F(xiàn)1模式對約2個井孔半徑內(nèi)的地層有明顯的響應(yīng)，而且在計算的頻率范圍內(nèi)，地層橫波速度始終是最主要控制參數(shù)，可考慮利用其對地層橫波速度進行層析成像.

對合成波形數(shù)據(jù)的分析表明，在快速地層條件下，W1波包可以給出準(zhǔn)確的地層橫波信息.需要注意的是，該波包包含了多種模式，包括井孔折射橫波、F1模式和F2模式與井孔折射橫波速度相接近的部分.結(jié)合靈敏度分析結(jié)果可知，W1波包基本不受地層各向異性的影響.另外，由于W1模式包含多種模式，其徑向探測特征較為復(fù)雜.對于本文研究的模型，由于該波包的頻率較高，其徑向探測深度較淺.

在慢速地層條件下，低頻地層螺旋波速度與地層橫波速度接近，不受鉆鋌波影響，是軟地層橫波測量的首選模式.本文的計算結(jié)果顯示，該模式受地層各向異性的影響較大.建議在實際處理中，如果考慮地層為各向異性地層，應(yīng)避免采用依賴正演模型的方法(Kimball, 1998; Geerits and Tang, 2003; Wang et al., 2010)，可以考慮采用基于數(shù)據(jù)驅(qū)動的方法(Huang and Yin, 2005; Tang et al., 2010;孫志峰等, 2011)對螺旋波的低頻趨勢線進行估算.另外，上述的模擬結(jié)果顯示，隨鉆四極子測井在快、慢速地層情況下的適用頻率不同，在實際儀器設(shè)計中應(yīng)考慮采用寬頻帶、大功率的發(fā)射器和高靈敏度的接收器，以提高采集數(shù)據(jù)的信噪比，擴大儀器測量的地層類型.

6 結(jié)論

數(shù)值研究了四極子聲源在含鉆鋌TI地層中產(chǎn)生的聲場，主要得到以下結(jié)論.

(1) 與偶極子聲場不同，在四極子聲源激發(fā)的模式波中未觀察到地層波和鉆鋌波相互耦合的現(xiàn)象.地層各向異性對鉆鋌波的影響很小，對于地層波的影響較大，且影響因素非常復(fù)雜.僅在一些特殊的頻率點處，地層模式波的控制因素相對較為簡單，如螺旋波的截止頻率處.

(2) 對于慢速地層，最低階的地層螺旋波對距離井壁約2個井孔半徑內(nèi)的地層有明顯的響應(yīng)，有希望對該范圍內(nèi)的地層橫波速度進行層析成像.

(3) 在快速地層井孔條件下，首至波包含了井孔折射橫波、F1和F2模式中與井孔折射橫波速度相接近的部分多種模式，其到時與橫波的到時相同，可以用于準(zhǔn)確提取地層橫波速度.而且，利用該波包獲得的地層橫波速度基本不受地層各向異性的影響.

(4) 慢度地層螺旋波受地層各向異性的影響較大，建議實際數(shù)據(jù)處理時，如考慮地層為各向異性地層，應(yīng)避免采用依賴正演模型的方法，可以考慮采用基于數(shù)據(jù)的處理方法.

致謝 美國殼牌國際石油公司黃山博士和中石化石油工程研究院吳金平博士對本文的完成提供了一些幫助，并提出了很好的建議和意見，在此表示感謝.

Chen X L, Wang R J. 2008. Investigating depth of mode waves in the borehole surrounded by transversely isotropic elastic formation.JournalofJilinUniversity(EarthScienceEdition) (in Chinese), 38(3): 502-507.Cui Z W. 2004.Theoretical and numerical study of modified Biot′s models, acousto-electric well logging and acoustic logging while drilling excited by multipole acoustic sources (in Chinese) [Ph. D. thesis]. Changchun: Jilin University.

Geerits T W, Tang X M, Hellwig O, et al. 2010. Multipole borehole acoustic theory: Source imbalances and the effects of an elastic logging tool.JournalofAppliedGeophysics, 70(2): 113-143.

Geerits T W, Tang X M. 2003. Centroid phase slowness as a tool for dispersion correction of dipole acoustic logging data.Geophysics, 68(1): 101-107.

Huang X J, Yin H Z. 2005. A data-driven approach to extract shear and compressional slowness from dispersive waveform data. ∥ SEG Technical Program Expanded Abstracts, 384-387.

Kimball C V. 1998. Shear slowness measurement by dispersive processing of the borehole flexural mode.Geophysics, 63(2): 337-344.

Kinoshita T, Takeshi T, Hiroshi N, et al. 2008. Next generation LWD sonic tool. ∥ 14th SPWLA Formation Evaluation Symposium of Japan, Japan.

Li X Q, Chen H, He X, et al. 2013. Analyses on mode waves of acoustic logging while drilling in transversely isotropic formations.ChineseJournalGeophysics(in Chinese), 56(9): 3212-3222, doi: 10.6038/cjg20130933.

Lu J Q, Ju X D, Qiao W X, et al. 2013. Application of digital signal processor to acoustic LWD tool.WellLoggingTechnology(in Chinese), 37(5): 527-530.

Sinha B K, Simsek E, Asvadurov S. 2009. Influence of a pipe tool on borehole modes.Geophysics, 74(3): E111-E123.

Sinha B K, Simsek E. 2010. Sonic logging in deviated wellbores in the presence of a drill collar. ∥ 2010 SEG Annual Meeting, Expanded Abstracts, Denver US, 553-557

Su Y D, Zhuang C X, Tang X M. 2011. LWD acoustic color mode wave attenuation character research and isolator design.ChineseJ.Geophys. (in Chinese), 54(9): 2419-2428, doi: 10.3969/j.issn.0001-5733.2011.09.026.

Su Y D, Tang X M, Zhuang C X, et al. 2013. Mapping formation shear-velocity variation by inverting logging-while-drilling quadrupole-wave dispersion data.Geophysics, 78(6): D491-D498.

Sun Z F, Fan G M, Meng J. 2011. A method of fitting function to extract phase slowness of dispersive wave from sonic logging data.JournalofJilinUniversity(EarthScienceEdition) (in Chinese), 41(4): 1239-1245.

Tang X M, Wang T, Patterson D. 2002. Multipole acoustic logging-while-drilling. ∥ 72nd Ann. Internat. Mtg., Soc. Explor. Geophys., Expanded Abstracts, 364-368.

Tang X M, Dubinsky V, Wang T, et al. 2003. Shear-velocity measurement in the logging-while-drilling environment: modeling and field evaluations.Petrophysics, 44(2): 79-89.Tang X M, Cheng C H. 2004. Quantitative Borehole Acoustic Methods. San Diego: Elsevier Science Publishing Co. Inc.

Tang X M, Li C, Patterson D J. 2010. A curve-fitting technique for determining dispersion characteristics of guided elastic waves.Geophysics, 75(3): E153-E160.

Wang B, Tao G, Chen T R, et al. 2012. Experimental study of acoustic logging while drilling measurements using an asymmetric source.ChineseJ.Geophys. (in Chinese), 55(1): 327-332, doi: 10.6038/j.issn.0001-5733.2012.01.033.

Wang H, Tao G, Wang B, et al. 2009. Wave field simulation and data acquisition scheme analysis for LWD acoustic tool.ChineseJ.Geophys. (in Chinese), 52(9): 2402-2409, doi: 10.3969/j.issn.0001-5733.2009.09.027.

Wang H, Tao G. 2011. Wavefield simulation and data-acquisition-scheme analysis for LWD acoustic tools in very slow formations.Geophysics, 76(3): E59-E68.

Wang T, Tang X M. 2003. LWD quadrupole shear measurement in anisotropic formations. ∥ Technical Program Expanded Abstracts, Dallas, Texas, 309-312.Wang R J, Qiao W X, Che X H. 2010. A new method of extracting

shear slowness from borehole mode waves. ∥ 2010 3rd International Congress on Image and Signal Processing (CISP), vol. 9, 4308-4312.Wang R J. 2012.Simulation research on acoustic logging in fluid-filled borehole surrounded by anisotropic formations (in Chinese) [Ph. D. thesis]. Beijing: China University of Petroleum.Xiao H B, Yang J Z, Ju X D, et al. 2009. Application of V-system in acoustic logging while drilling data denoising.JournalofChinaUniversityofPetroleum(EditionofNaturalScience) (in Chinese), 33(2): 58-62.

Yang Y, Che X H, Li J, et al. 2009. Sound isolation numerical simulation on isolator of logging while drilling sonic tool using finite difference time-domain method.JournalofChinaUniversityofPetroleum(EditionofNaturalScience) (in Chinese), 33(3): 66-70.Yang J Q, Sinha B, Valero H P, et al. 2011. Estimation of the formation shear and borehole fluid slownesses using sonic dispersion data in the presence of a drill collar. ∥ SEG Technical Program Expanded Abstracts, Vol. 30, 464-468.

Zhang Z X, Mochida M, Kubota M, et al. 2010. Shear slowness estimation by inversion of LWD borehole quadrupole mode. ∥SEG Technical Program Expanded Abstracts, 29: 528-532.

Zhu Z Y, Toks?z M N, Rao M, et al. 2008. Experimental studies of monopole, dipole, and quadrupole acoustic logging while drilling (LWD) with scaled borehole models.Geophysics, 73(4): E133-E143.

附中文參考文獻

陳雪蓮, 王瑞甲. 2008. 橫向各向同性彈性地層井孔中模式波的探測深度. 吉林大學(xué)學(xué)報: 地球科學(xué)版, 38(3): 502-507.

崔志文. 2004. 多孔介質(zhì)聲學(xué)模型與多極源聲電效應(yīng)測井和多極隨鉆聲測井的理論與數(shù)值研究[博士論文]. 長春: 吉林大學(xué).

李希強, 陳浩, 何曉等. 2013. 橫向各向同性地層中隨鉆聲波測井模式波分析. 地球物理學(xué)報, 56(9): 3212-3222, doi: 10.6038/cjg20130933.

盧俊強, 鞠曉東, 喬文孝等. 2013. 數(shù)字信號處理器在隨鉆聲波測井儀中的應(yīng)用. 測井技術(shù), 37(5): 527-530.

蘇遠(yuǎn)大, 莊春喜, 唐曉明. 2011. 隨鉆聲波測井鉆鋌模式波衰減規(guī)律研究與隔聲體設(shè)計. 地球物理學(xué)報, 54(9): 2419-2428, doi: 10.3969/j.issn.0001-5733.2011.09.026.

孫志峰, 樊官民, 孟杰. 2011. 一種提取聲波測井頻散波相慢度的適應(yīng)函數(shù)方法. 吉林大學(xué)學(xué)報(地球科學(xué)版), 41(4): 1239-1245.

王兵, 陶果, 陳天潤等. 2012. 非對稱聲源多極子隨鉆聲波測井實驗室測量研究. 地球物理學(xué)報, 55(1): 327-332, doi: 10.6038/j.issn.0001-5733.2012.01.033.

王華, 陶果, 王兵等.2009. 多極子隨鉆聲波測井波場模擬與采集模式分析. 地球物理學(xué)報, 52(9): 2402-2409, doi: 10.3969/j.issn.0001-5733.2009.09.027.

王瑞甲. 各向異性地層聲波測井模擬研究 [博士論文]. 北京: 中國石油大學(xué), 2012.

肖紅兵, 楊錦舟, 鞠曉東等. 2009. V 系統(tǒng)在隨鉆聲波測井?dāng)?shù)據(jù)降噪中的應(yīng)用. 中國石油大學(xué)學(xué)報(自然科學(xué)版), 33(2): 58-62.

楊勇, 車小花, 李俊等. 2009. 基于時域有限差分法的隨鉆聲波測井儀隔聲體隔聲效果的數(shù)值模擬. 中國石油大學(xué)學(xué)報(自然科學(xué)版), 33(3): 66-70.

(本文編輯 胡素芳)

Numerical study on quadrupole acoustic LWD in VTI formations

WANG Rui-Jia1,2, QIAO Wen-Xiao1,2*

1StateKeyLaboratoryofPetroleumResourcesandProspecting,ChinaUniversityofPetroleum,Beijing102249,China2KeyLaboratoryofEarthProspectingandInformationTechnology,Beijing102249,China

Acoustic logging while drilling (LWD) is able to provide real-time acoustic properties of rock formations near bits, and therefore has an incomparable advantage over acoustic wireline logging (WL) in saving wellbore occupation time and avoiding risks. Current quadrupole acoustic LWD is capable of measuring formation shear velocity in arbitrary rock formations. However, the quadrupole acoustic LWD theory is not complete yet, i.e., there is no clear understanding on some fundamentals, which required further study.We derive the dispersion equations for the waveguide consisting of a fluid-filled borehole surrounded by a VTI formation, and numerically calculate the dispersion, amplitude and sensitivity functions of the borehole modes excited by a quadrupole acoustic source. We also synthesize the array waveforms recorded by the receivers in the borehole. The influences of formation anisotropy on borehole quadrupole modes are investigated.Numerical results show that, the formation anisotropy has negligible effects on the collar modes, while it strongly affects the formation modes. Only at a few frequencies, the control factor is simple, such as at the cut-off frequency of screw waves. The formation quadrupole modes of the lowest order show significant response to the formation shear velocity within approximate 2 borehole radii away from the sidewall, and it can be utilized to perform the radial tomography of shear velocity in this radial range. For fast formations, the first arrival includes borehole refracted shear mode, part of F1 and F2 modes, the velocities of which are close to that of shear modes, and it can be utilized to accurately evaluate the formation shear velocity without the influence of formation anisotropy. In slow formations, the formation screw modes are highly influenced by anisotropy.It suggests that, in field data processing, for fast formations, the first arrived wave packet can be utilized to extract shear wave velocity, while for slow formations, the formation screw waves can be processed by data-driven methods to obtain true shear wave velocity of rock formations. Because the operation frequencies for fast and slow formations are different, broadband and high-powered transmitters and high-sensitivity receivers should be adopted in the design of an acoustic LWD tool, to increase the signal-to-noise ratio and to improve application range of the tool.

Quadrupole； LWD； Acoustic logging； Shear waves

國家自然科學(xué)基金項目(11204380、11374371、11134011和61102102)，國家油氣重大科技專項(2011ZX05020-009)，中國石油科技創(chuàng)新基金(2013D-5006-0304、2014D-5006-0307)和中國石油天然氣集團公司項目(2014B-4011、2014D-4105、2014A-3912)聯(lián)合資助.

王瑞甲，男，山東無棣人，漢族，博士，主要從事地球物理測井方法、檢測聲學(xué)等方面的研究工作.E-mail：wruijia@foxmail.com

*通訊作者喬文孝，男，山東禹城人，博士、博士生導(dǎo)師.研究方向：聲波測井、檢測聲學(xué). E-mail:qiaowx@cup.edu.cn

10.6038/cjg20150820.

10.6038/cjg20150820

P631

2014-06-21，2015-03-25收修定稿

王瑞甲，喬文孝. 2015. VTI地層隨鉆四極子聲波測井?dāng)?shù)值研究.地球物理學(xué)報，58(8):2862-2872,

Wang R J, Qiao W X. 2015. Numerical study on quadrupole acoustic LWD in VTI formations.ChineseJ.Geophys. (in Chinese),58(8):2862-2872,doi:10.6038/cjg20150820.