TI地層聲學(xué)各向異性隨鉆評價方法的數(shù)值研究

王瑞甲,喬文孝,車小花,鞠曉東

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

TI地層聲學(xué)各向異性隨鉆評價方法的數(shù)值研究

王瑞甲1,2,喬文孝1,2,車小花1,2,鞠曉東1,2

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

采用三維有限差分方法模擬偶極子聲源在含鉆鋌的各向異性地層井孔中激發(fā)的聲場,研究地層各向異性對偶極模式波頻散特征的影響,并提出從隨鉆四分量偶極子聲波數(shù)據(jù)中提取地層快、慢橫波慢度及其各向異性的方法。結(jié)果表明,偶極子聲源在各種地層井孔中均可激發(fā)對橫波慢度的方位各向異性敏感的偶極模式波:對于快速地層, F2模式的低頻部分對地層橫波慢度及其方位各向異性最為敏感;對于中速和慢速地層,彎曲波較高頻率部分對地層橫波慢度較為敏感。采用各向同性地層建模的方法可以較為準確地提取地層的快、慢橫波慢度,提取結(jié)果的誤差小于4%。

隨鉆測井;聲波測井;各向異性;偶極子

隨鉆聲波測井是一種發(fā)展迅速,應(yīng)用廣泛的聲學(xué)測量技術(shù)[1-4]。在橫波測量中,人們選用了四極子聲源進行測量[1]。但是,限于四極子聲源的輻射指向性,難以利用其獲得地層橫波的各向異性信息。由于隨鉆四極子聲波測井儀換能器的安裝方式與正交偶極子聲波測井儀相同,且其接收站具備正交偶極子接收功能,通過合理的電路設(shè)計,可以較為方便地實現(xiàn)正交偶極子聲波測井,所以隨鉆正交偶極子測井是進行地層各向異性評價首選方法。關(guān)于偶極子聲源在各向異性地層中激發(fā)的聲場的研究大部分集中在電纜測井條件下的各向異性地層井孔中的聲場特征方面,主要包括3種情況:VTI地層豎直井[5-9]、TI地層斜井和水平井[10-14]和正交各向異性地層井孔[15]。有關(guān)隨鉆條件下偶極子聲源在各向異性地層中激發(fā)的聲場的研究相對較少。Sinha 等[16]研究了多極子聲源在含鉆鋌各向異性地層斜井中激發(fā)的聲場,認為偶極模式波可以用于提取過井軸的兩個正交的平面內(nèi)的剪切模量。王瑞甲等[17-18]采用三維有限差分方法研究了正交偶極子聲源在含鉆鋌TI地層斜井中激發(fā)的聲場,證實了隨鉆條件下偶極子聲源對于地層各向異性的方位敏感性。筆者重點研究地層各向異性的存在對于偶極模式波頻散特征的影響以及從四分量偶極子數(shù)據(jù)中提取快、慢橫波慢度的算法。

1 模擬方法

三維有限差分方法是進行聲場模擬的常用數(shù)值方法。在直角坐標系下,采用應(yīng)力和速度表示的運動方程和本構(gòu)方程[6]分別為

式中,vx、vy、vz分別為x、y、z方向上的速度分量;τxx、τyy、τzz分別為x、y、z方向上的正應(yīng)力;τxy、τyx、τxz為剪切應(yīng)力;ρ為介質(zhì)的密度;gab(a、b=x～z)應(yīng)表示力變化速度的體積源,與體力源fi(i=x～z)組合使用可以模擬各種聲源;C為介質(zhì)的彈性常數(shù)矩陣。網(wǎng)格劃分方式采用中心差分格式的交錯網(wǎng)格實現(xiàn),聲學(xué)模型邊界采用完全匹配層(PML)技術(shù)。由于計算模型較大,采用基于OpenMP和MPI的混合編程技術(shù)實現(xiàn)算法的并行。所采用的有限差分算法的具體實現(xiàn)方法見文獻[17]。

本文中所關(guān)注的地層類型為橫向各向同性(TI)地層,其中地層的TI旋轉(zhuǎn)對稱軸與井軸(z軸)垂直。這種地層可采用5個獨立的彈性參數(shù)描述,

將各向異性地層看作是各向同性地層的微擾,把矩陣C分解為表征參考的各向同性介質(zhì)的矩陣CISO和表征擾動的矩陣C′,即

對于偶極子聲源而言,不同偏振方向的聲源激發(fā)的模式波不同,所以選取的參考各向同性地層的參數(shù)不同。對于HTI地層(井軸與TI旋轉(zhuǎn)對稱軸垂直),選取兩套各向同性地層參數(shù),分別對應(yīng)不同偏振方向的偶極子聲源。對與SH波偏振方向一致的偶極子聲源而言,參考的各向同性地層的縱波速度選擇為,橫波速度選擇為對與SV波偏振方向一致的偶極子聲源,對應(yīng)的各向同性地層的縱、橫波速度分別選擇為和通過對比偶極子聲源在各向異性地層與對應(yīng)的各向同性地層井孔內(nèi)激發(fā)的波形以及模式波的頻散曲線,分析地層各向異性的存在(C′)對于模式波的影響,并進一步討論采用各向同性地層建模從彎曲波中提取快、慢橫波慢度的方法的可行性。

2 數(shù)值模擬結(jié)果及分析

2.1 數(shù)值模型

圖1為含鉆鋌各向異性地層井孔聲學(xué)模型示意圖。該模型為典型的柱狀徑向分層介質(zhì)模型,模型從內(nèi)到外依次為水層、鋼(鉆鋌)、水層和地層,其外徑分別為0.054、0.18、0.24 m和無窮大。模型主計算區(qū)域大小為1 m×1 m×4.8 m,x、y和z方向的空間采樣步長分別為0.005、0.005和0.0125 m。井軸與坐標軸z平行。地層為橫向各向同性介質(zhì),其介質(zhì)對稱軸為S,平行于坐標軸x,與井軸垂直。定義偶極子聲源的偏振方向同y坐標軸的夾角為β。特別地,當(dāng)β=0°時,偶極子源偏振方向同地層快橫波(SH波)偏振方向一致;當(dāng)β=90°時,偶極子源偏振方向同地層慢橫波(SV波)偏振方向一致。聲源加載在距離底界面0.8 m處,采用在鉆鋌外徑處加載兩個振動相位相反的點聲源的方法模擬偶極子聲源。偶極子接收器同樣置于鉆鋌外徑處,源距范圍為2.0～3.5 m,接收器間距為0.15 m。由于本文中重點研究地層模式波,不涉及隔聲及鉆鋌波問題的研究,為壓制鉆鋌波的幅度,在發(fā)射器到源距最小的接收器之間將鉆鋌截斷。

為了分析不同地層情況下偶極子聲源激發(fā)的彎曲波的頻散特性,選取快速、中速、慢速3種地層進行模擬?？焖俚貙硬捎肅otton Vally Shale的彈性參數(shù)[10],慢速地層采用Austin Chalk的彈性參數(shù)[10],中速地層采用實驗室內(nèi)用酚醛樹脂材料制作的模型的TI彈性參數(shù)測量值[17]。數(shù)值模擬中采用的各種地層的彈性參數(shù)見表1。模擬時采用的鉆鋌材質(zhì)的縱、橫波速度分別為5.860和3.130 km/s,密度為7.860×103kg/m3;鉆鋌內(nèi)外流體的縱波速度為1.500 km/s,密度為1.000×103kg/m3。

圖1 含鉆鋌各向異性地層井孔聲學(xué)模型示意圖Fig.1 Schematic diagram of acoustic model for logging while drilling in a HTI formation

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

2.2 快速地層

快速地層情況下,偶極子聲源激發(fā)的聲場的模擬結(jié)果如圖2所示,其中(a)、(b)為β=0°的偶極子聲源激發(fā)的波形及頻散分析結(jié)果,(c)、(d)為β= 90°的偶極子聲源激發(fā)的波形及頻散分析結(jié)果。頻散分析方法采用了文獻[19]中介紹的方法。圖中,實線表示各向異性地層的計算結(jié)果,虛線為各向同性地層結(jié)果,“SH”和“SV”標注的短線分別表示地層SH波和SV波的慢度,“S”標注表示橫波,“F”標注表示依賴地層傳播的偶極子模式波,“C”標注表示依賴鉆鋌傳播的偶極子模式波。在頻率慢度相關(guān)圖中,未標注的相關(guān)條帶均為噪聲信號。

從圖2(a)可見,在各向異性和各向同性地層情況下,偶極子波形的形態(tài)基本一致:波形初至部分幾乎完全重合;對于后續(xù)波部分,各向同性地層情況下的波形到時略早于各向異性地層情況。在圖2(c) 中,兩種地層情況下的偶極子波形幾乎完全重合。上述計算結(jié)果表明,地層各向異性的存在對于快速地層井孔內(nèi)激發(fā)的偶極子聲場的影響較小。

進一步通過頻散分析結(jié)果研究地層各向異性對模式波頻散曲線的影響。從圖2(b)和(d)可觀察到兩種偶極模式波,分別標記為F1和F2。文獻[2, 18,20]的研究結(jié)果表明,這兩種偶極子波為最低階的液浸無限大鉆鋌中的偶極模式波和無鉆鋌充液井孔中最低階偶極模式波耦合而成的兩種新的模式波。文獻[18]通過靈敏度分析的方法對此現(xiàn)象進行了研究,認為兩種模式波在低頻和高頻部分的影響因素不同。對于F1模式,其低頻部分的頻散特性由鉆鋌控制,高頻部分主要受地層參數(shù)影響;對于F2模式,其低頻部分由地層參數(shù)控制,高頻部分由鉆鋌參數(shù)控制。從圖2(b)可見,對于F1的低頻部分(3 kHz以下)和F2的高頻部分(4 kHz以上),各向異性地層和各向同性地層情況下的頻散曲線基本重合,這說明地層各向異性的存在對其影響較小;而對于F1的高頻部分(4 kHz以上)和F2的低頻部分(4 kHz以下),地層各向異性的存在對模式波的頻散規(guī)律則有一定的影響。圖2(d)也有類似的現(xiàn)象。結(jié)合文獻[18]的研究結(jié)果,本文中認為該現(xiàn)象是由于模式波在低頻段和高頻段的影響因素不同造成的:依賴于鉆鋌傳播的部分(F1的低頻部分和F2的高頻部分)受地層參數(shù)影響很小,從而受地層各向異性影響也很小;而依賴地層傳播的部分(F1的高頻部分和F2的低頻部分),則受地層各向異性的影響相對較大。

另外,從圖2(b)和(d)可見,對于不同方向的偶極子聲源和不同的地層,F2低頻部分的慢度均趨近于對應(yīng)的地層體橫波慢度,這說明對于快速地層可以通過F2模式的低頻部分獲得快、慢橫波的慢度。圖2(b)中,在2～4 kHz頻率范圍內(nèi),地層各向異性的存在使得F2模式的頻散曲線向低頻有較大偏移,這使得圖2(a)中兩種地層情況下的偶極子波形存在差異;圖2(d)中,在該頻率范圍內(nèi),地層各向異性的存在使得F2模式的頻散曲線也向低頻偏移,但由于偏移值較小,在圖2(c)中不能夠明顯觀察到。

圖2 偶極子聲源在含鉆鋌快速地層井孔中激發(fā)的波形及頻散分析結(jié)果Fig.2 Waveforms and dispersion analysis results in fluid-filled borehole with a collar surrounded by fast formations excited by dipole sources

2.3 中速地層

圖3是中速地層情況下的數(shù)值模擬結(jié)果,其中(a)、(b)采用的聲源為β=0°的偶極子聲源,(c)、(d)為β=90°的偶極子聲源。由圖3(a)和(c)可見,各向異性地層情況下的彎曲波與對應(yīng)的各向同性地層情況下的彎曲波的形態(tài)基本一致,到時比各向同性地層略晚,且在β=0°時,兩種地層情況下波形的到時差別較大,在β=90°時,差別較小。這說明,對于該地層模型,各向異性的存在使得彎曲波傳播速度略降低,且對于β=0°的偶極子聲源而言更為明顯。從圖3(b)和(d)可見,在低頻段(0～1.2 kHz),基本觀察不到地層各向異性的存在對于彎曲波頻散特性的影響,而在高頻段(1.5～7 kHz),地層各向異性的存在使得彎曲波的慢度增加。這是因為低頻部分的彎曲波主要由鉆鋌參數(shù)控制,而高頻部分則主要受地層參數(shù)的影響。在利用彎曲波反演地層橫波慢度時,應(yīng)選擇對地層參數(shù)敏感的頻率段。對于該地層而言,選擇1.5 kHz以上的頻率是合適的。另外,與快速地層的計算結(jié)果相比,在圖3(b)和(d)中觀察到了地層六極子波(Hexpole mode),但未觀察到其與地層各向異性相關(guān)的信息。

2.4 慢速地層

圖4為慢速地層井孔中偶極子聲源激發(fā)的聲場。從圖4(a)和(c)可見,對于慢速地層的情況,各向異性地層和對應(yīng)的各向同性地層井孔情況下的彎曲波的形態(tài)基本一致,到時略有差別。地層各向異性的存在使得β=0°時彎曲波的到時滯后,而β= 90°時彎曲波的到時提前。從圖4(b)和(d)可見,在低頻段,各向異性地層和各向同性地層兩種情況下的彎曲波頻散曲線基本重合,但隨著頻率增加,二者出現(xiàn)差異,并越來越大。對于β=0°的情況,地層各向異性的存在使得彎曲波慢度增加;對于β=90°的情況,地層各向異性的存在使得彎曲波的慢度降低。如果采用各向同性地層建模,通過頻散校正的方法來計算快、慢橫波的慢度,對于β=0°的情況,反演得到的快橫波慢度會大于真實快橫波慢度;而對于β=90°的情況,反演得到的慢橫波慢度會小于真實的慢橫波慢度。反演結(jié)果的誤差與計算采用的頻率有關(guān),對于該地層而言,頻率越高,誤差越大。選擇恰當(dāng)?shù)念l率可以降低反演結(jié)果的誤差。

圖3 偶極子聲源在含鉆鋌中速地層井孔中激發(fā)的波形及頻散分析結(jié)果Fig.3 Waveforms and dispersion analysis results in fluid-filled borehole with a collar surrounded by intermediate speed formations excited by dipole sources

圖4 偶極子聲源在含鉆鋌慢速地層井孔中激發(fā)的波形及頻散分析結(jié)果Fig.4 Waveforms and dispersion analysis results in fluid-filled borehole with a collar surrounded by slow formations excited by dipole sources

3 各向異性提取方法及數(shù)值驗證

3.1 隨鉆情況下的地層各向異性提取方法

在隨鉆測量條件下,由于井孔中鉆鋌的存在,傳統(tǒng)的基于時域的各向異性提取方法不適用于地層各向異性參數(shù)的提取。本文中提出基于正演模型的隨鉆各向異性提取方法。為簡化反演過程的復(fù)雜度,以各向異性地層快、慢橫波各自對應(yīng)的各向同性地層建立隨鉆測井聲學(xué)模型。在各向同性地層井孔中,彎曲波的慢度取決于地層的密度和縱、橫波慢度、鉆鋌的密度和縱、橫波慢度以及鉆鋌內(nèi)外流體的密度及縱橫波慢度,計算公式為

式中,ω為圓頻率;α為該頻率下彎曲波的慢度;

分別為地層的密度及縱、橫波慢度;ρcol、和SScol分別為鉆鋌的密度及縱、橫波慢度;ρf和分別為鉆鋌內(nèi)外流體的密度及縱橫波慢度;r0、r1和r2分別為鉆鋌內(nèi)、外半徑和井孔半徑。由于實際的鉆鋌是固定的,鉆鋌密度ρcol、鉆鋌縱波慢度鉆鋌橫波慢度鉆鋌內(nèi)外半徑r0和r1及井孔半徑r2均可以看作已知量。另外,鉆鋌內(nèi)外流體的密度及慢度ρf和可以根據(jù)實際情況確定。其他參數(shù),包括地層的密度ρform可以根據(jù)密度測井資料確定,地層的縱波慢度可以根據(jù)單極子縱波的測量結(jié)果確定。據(jù)此,可得到彎曲波的慢度與地層橫波慢度的關(guān)系:

特別地,對于快橫波,選擇為快橫波慢度,對于慢橫波,選擇為慢橫波慢度,其計算公式分別為

對于測得的四分量偶極子波形,首先通過Alford四分量偶極子波形旋轉(zhuǎn)方法獲取快、慢波波形,然后通過波場分離的方法獲得較為純凈的地層偶極子波?？?、慢橫波慢度提取的目標函數(shù)為

式中,(ω)和(ω)分別為快、慢地層偶極子波波形的相位譜;d為接收器源距;m為乘因子;Weight(ω)為加權(quán)系數(shù);ωmin和ωmax分別為計算的最低頻率和最高頻率。有關(guān)加權(quán)系數(shù)的構(gòu)造方法見文獻[21]。使得目標函數(shù)R()和R()取得極大值的和即為所求的快、慢橫波慢度。橫波慢度的各向異性值可以近似地采用下式確定:

綜上所述,從隨鉆正交偶極子資料中提取地層各向異性的算法可以分為以下幾步:

(1)理論計算不同情況下含鉆鋌井孔中偶極子波的頻散曲線,并將數(shù)據(jù)存儲于硬盤中。

(2)通過Alford四分量偶極子波形旋轉(zhuǎn)的方法和波場分離的方法獲得快、慢偶極子波波形。

(3)從硬盤中讀取偶極子波的頻散曲線,并采用插值的方法對頻散曲線進行重構(gòu),然后采用式(11)～(13),獲取地層的快、慢橫波慢度和各向異性值,并確定快橫波方位角。

3.2 模擬數(shù)據(jù)的處理結(jié)果

圖5為快速地層井孔中模擬的偶極子波形數(shù)據(jù)的處理結(jié)果,其中圖5(a)為不同橫波慢度所對應(yīng)的偶極模式波(F2)頻散曲線,圖5(b)為根據(jù)式(11)、(12)計算得到的目標函數(shù)。圖5(a)中的偶極子波的頻散曲線對應(yīng)圖2中的F2模式。從圖5(a)可見, F2模式在低頻段對于橫波慢度變化最敏感,隨頻率增加,敏感性降低。另外,隨著橫波慢度增加,F2模式的頻散曲線向低頻偏移。由此可知,對于快速地層,選擇低頻部分的數(shù)據(jù)可以提高反演結(jié)果的可靠性。對于該地層模型,采用1～2.5 kHz的頻率范圍是合適的。

從圖5(b)可見,目標函數(shù)極大值所對應(yīng)的值,即計算得到的快、慢橫波的慢度值分別為301.1和346.9 μs/m,地層真實的快、慢橫波慢度分別為296.7和346.0 μs/m,反演得到的快、慢橫波慢度誤差分別為1.5%和0.26%。提取的橫波慢度的各向異性值為14.1%,小于真實的各向異性值15.3%。結(jié)合圖2可知,該誤差是由于在反演過程中采用了各向同性地層建模而引入的。

圖5 快速地層模型處理結(jié)果Fig.5 Results of fast formation model

同樣,采用該方法處理了中速地層偶極子聲波波形數(shù)據(jù),如圖6所示。從圖6(a)可見,彎曲波在低頻段對地層橫波慢度變化不敏感,在稍高的頻率段對地層橫波慢度敏感。結(jié)合圖3中的頻散分析結(jié)果,選取理論的彎曲波頻散曲線極小值所對應(yīng)的頻率點作為計算的下限頻率ωmin,計算的頻帶寬度為1 kHz。由于不同橫波慢度下,理論頻散曲線極小值所對應(yīng)的頻率是變化的,所以計算目標函數(shù)所選用的頻率范圍也是變化的。從圖6(b)可見,反演得到的快、慢橫波的慢度分別為688.4和733.1 μs/m。地層真實的快、慢橫波慢度分別為666.4和714.7 μs/m。反演結(jié)果的誤差分別為3.3%和2.6%。提取的橫波慢度的各向異性值為6.3%,略小于真實的各向異性值7.0%。

圖6 中速地層模型處理結(jié)果Fig.6 Results of intermediate speed formation model

圖7為慢速地層井孔偶極子數(shù)據(jù)的處理結(jié)果。由圖4可見,頻率越高,各向異性地層和各向同性地層情況下的彎曲波頻散曲線差別越大。為了盡量減小地層的各向異性對計算結(jié)果的影響,與中速地層一致,選取理論彎曲波頻散曲線極小值所對應(yīng)的頻率點作為計算的下限頻率ωmin,計算的頻帶寬度為1 kHz。從圖7(b)可知,反演得到的快、慢橫波的慢度分別為873.2和954.8 μs/m,地層真實的快、慢橫波的慢度分別為842.4和957.4 μs/m,誤差分別為3.7%和0.27%。該結(jié)果與圖4中的頻散分析結(jié)果一致。另外,提取的橫波慢度的各向異性值為8.9%,真實的各向異性值為12.8%。地層各向異性的存在使得快彎曲波慢度增加,慢彎曲波慢度減小,采用各向同性地層建模的方法提取的各向異性值與真實值有一定差別。

從上述處理結(jié)果可以看出,從模擬數(shù)據(jù)中反演得到的慢橫波的誤差小于快橫波的誤差,這是由于本文中研究的3種地層各向異性的存在對慢彎曲波頻散特性的影響大于其對快彎曲波的影響。這也可以在圖2～4的頻散圖中觀察到。綜上所述,采用各向同性地層建模計算地層快、慢橫波慢度,誤差小于4%,說明該方法是可行的。

圖7 慢速地層模型處理結(jié)果Fig.7 Results of slow formation model

4 結(jié)果討論

正演模型中輸入?yún)?shù)的誤差可能導(dǎo)致反演得到的快、慢橫波慢度和各向異性值出現(xiàn)較大的誤差,結(jié)合隨鉆單極子和四極子的測量結(jié)果進行聯(lián)合反演可以在很大程度上保證計算結(jié)果的可靠性。

在隨鉆測井中,另外一個很重要的問題是實時的數(shù)據(jù)處理和傳輸。四分量偶極子波形數(shù)據(jù)量較大,目前的遙測技術(shù)無法實現(xiàn)波形數(shù)據(jù)的實時傳輸,各向異性提取程序需要在井下進行,從而限制了算法的計算復(fù)雜度。由圖5～7可以看出,在一定的頻率段內(nèi),彎曲波慢度隨地層慢度的變化呈線性。因此,可通過近似公式計算地層各向異性值

式中,ΔSS為快、慢橫波的慢度差;ΔSflex為快、慢彎曲波的慢度差;為主頻附近彎曲波慢度相對于地層橫波慢度的靈敏度。ΔSflex可以通過時域內(nèi)的彎曲波慢度計算或者快、慢波匹配等方法得到。則可以通過理論計算得到,并存儲于井下存儲器中。進一步采用式(13)可近似得到地層橫波慢度的各向異性值。

本文中未涉及鉆井液循環(huán)及鉆進過程中的噪聲、儀器的偏心和鉆鋌波等對測量結(jié)果的影響。這些因素可能對測量結(jié)果有較大的影響。針對這些問題,需要在儀器設(shè)計及信號處理方面做進一步的研究。

5 結(jié) 論

(1)在隨鉆條件下,偶極子聲源在各種地層井孔中均可激發(fā)對橫波慢度的方位各向異性敏感的偶極模式波。對于快速地層,F2模式的低頻部分對地層橫波慢度及其方位各向異性最為敏感;對于中速和慢速地層,彎曲波的高頻部分對對地層橫波慢度較為敏感。

(2)偶極子聲源在各向異性地層和對應(yīng)的各向同性地層井孔中激發(fā)的偶極子波形,形態(tài)相同,慢度略有差別。

(3)采用各向同性地層建模的方法可以較為準確地提取地層的快、慢橫波慢度,提取結(jié)果的誤差小于4%。

[1] TANG X M,CHENG C H.Quantitative borehole acoustic methods[M].San Diego:Elsevier Science Publishing Co Inc,2004.

[2] SINHA B K,SIMSEK E,ASVADUROV S.Influence of a pipe tool on borehole modes[J].Geophysics,2009,74 (3):E111-E123.

[3] 王華,陶果,王兵,等.多極子隨鉆聲波測井波場模擬與采集模式分析[J].地球物理學(xué)報,2009,52(9): 2402-2409.

WANG Hua,TAO Guo,WANG Bing,et al.Wave field simulation and data acquisition scheme analysis for LWD acoustic tool[J].Chinese Journal Geophysics,2009,52 (9):2402-2409.

[4] 陳雪蓮,魏周拓.隨鉆單極子聲反射測井?dāng)?shù)值模擬[J].石油學(xué)報,2012,33(5):835-840.

CHEN Xuelian,WEI Zhoutuo.Numerical simulation of monopole acoustic reflection imaging logging in the logging-while-drilling condition[J].Acta Petrolei Sinica, 2012,33(5):835-840.

[5] SCHMITT D P.Acoustic multipole logging in transversely isotropic poroelastic formations[J].JASA,1989,86 (6):2397-2421.

[6] 張海瀾,王秀明,張碧星.井孔中的聲場與波[M].北京:科學(xué)出版社,2004.

[7] 陳雪蓮,王瑞甲.徑向分層TI孔隙介質(zhì)井孔中激發(fā)的模式波的數(shù)值研究[J].地球物理學(xué)報,2008,51(4): 1270-1277.

CHEN Xuelian,WANG Ruijia.A numerical study on the mode waves excited by multipole sources in the fluidfilled borehole in radially layered transversely isotropic porous medium[J].Chinese Journal Geophysics,2008, 51(4):1270-1277.

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

CHEN Xuelian,WANG Ruijia.Investigating depth of mode waves in the borehole surrounded by transversely isotropic elastic formation[J].Journal of Jilin University (Earth Science Edition),2008,38(3):502-507.

[9] HE X,HU H S.Borehole flexural modes in anisotropic formations:the low-frequency asymptotic velocity[J]. Geophysics,2009,74(4):E149-E158.

[10] SINHA B K.Elastic-wave propagation in deviated wells in anisotropic formations[J].Geophysics,2006,71 (6):D191-D202.

[11] HE X,HU H S,GUAN W,Fast and slow flexural waves in a deviated borehole in a homogeneous or layered anisotropic formation[J].Geophys J Int,2010, 181(1):417-426.

[12] 閆守國,宋若龍,呂偉國,等.橫向各向同性地層斜井中正交偶極子激發(fā)的聲場的數(shù)值模擬[J].地球物理學(xué)報,2011,54(9):2412-2418.

YAN Shouguo,SONG Ruolong,Lü Weiguo,et al.Numerical simulation of acoustic field excited by cross-dipole source in deviated wells in transversely isotropic formation[J].Chinese Journal Geophysics,2011,54 (9):2412-2418.

[13] 王瑞甲,喬文孝,車小花,等.TI地層井孔中正交偶極子聲波測井的物理模擬[J].中國石油大學(xué)學(xué)報:自然科學(xué)版,2012,36(6):62-69.

WANG Ruijia,QIAO Wenxiao,CHE Xiaohua,et al. Physical simulation of cross-dipole acoustic logging in borehole surrounded by TI formation[J].Journal of China University of Petroleum(Edition of Natural Science), 2012,36(6):62-69.

[14] 王瑞甲,喬文孝,鞠曉東,等.不同偏振方向的偶極子聲源在HTI地層井孔中激發(fā)聲場的實驗研究[J].地球物理學(xué)報,2013,56(2):707-717.

WANG Ruijia,QIAO Wenxiao,JU Xiaodong,et al. Experimental study of the acoustic field in the borehole surrounded by HTI formations excited by dipole sources with different orientations[J].Chinese Journal Geophysics,2013,56(2):707-717.

[15] CHENG N Y.Borehole wave propagation in isotropic and anisotropic media:three-dimensional finite difference approach[D].Cambridge:Department of Earth, Atmospheric,and Planetary Sciences,MIT,1994.

[16] SINHA B K,SIMSEK E.Sonic logging in deviated wellbores in the presence of a drill collar:2010 SEG Annual Meeting,Expanded Abstracts,Denver US,October 17-22,2010[C].US Tulsa:SEG,c2010.

[17] 王瑞甲,喬文孝,鞠曉東.利用隨鉆正交偶極子聲波測井評價地層各向異性的數(shù)值研究[J].地球物理學(xué)報,2012,55(11):3870-3882.

WANG Ruijia,QIAO Wenxiao,JU Xiaodong.Numerical study of formation anisotropy evaluation using cross dipole acoustic LWD[J].Chinese Journal Geophysics, 2012,55(11):3870-3882.

[18] 王瑞甲.各向異性地層聲波測井模擬研究[D].北京:中國石油大學(xué)地球物理與信息工程學(xué)院,2012.

WANG Ruijia.Simulation research on acoustic logging in fluid-filled borehole surrounded by anisotropic formations[D].Beijing:College of Geophysics and Information Engineering in China University of Petroleum, 2012.

[19] 王瑞甲,喬文孝,鞠曉東.一種多通道聲波測井信號頻散分析方法[J].測井技術(shù),2012,36(2):135-140.

WANG Ruijia,QIAO Wenxiao,JU Xiaodong.A multichannel acoustic logging signal dispersion analysis method[J].Well Logging Technology,2012,36(2):135-140.

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

CUI Zhiwen.Theoretical and numerical study of modified Biot's models,acousto-electric well logging and acoustic logging while drilling excited by multipole acoustic sources[D].Changchun:College of Physics,Jilin University,2004.

[21] WANG R J,QIAO W X,CHE X H.A new method of extracting shear slowness from borehole mode waves: CISP'10:proceedings of the 2010 3nd International Congress on Image and Signal Processing,Yantai,China, October 16-18,2010[C].US Piscataway:IEEE, c2010.

(編輯 修榮榮)

Numerical study of anisotropy evaluation method by acoustic logging while drilling in TI formations

WANG Ruijia1,2,QIAO Wenxiao1,2,CHE Xiaohua1,2,JU Xiaodong1,2
(1.State Key Laboratory of Petroleum Resources and Prospecting in China University of Petroleum,Beijing 102249,China; 2.Key Laboratory of Earth Prospecting and Information Technology in China University of Petroleum,Beijing 102249,China)

The acoustic field excited by dipole sources in a fluid-filled borehole with a collar surrounded by anisotropic and the corresponding isotropic formations was simulated by a three-dimension finite difference method.The effects of formation anisotropy on dispersions of borehole modes were investigated.The results show that the dipole source can excite dipole modes which are sensitive to the formation anisotropy in all types of formations.For fast formation,the low-frequency part of F2 mode is most sensitive to shear velocities and anisotropy of the formation,while for intermediate speed and slow formations,the higher frequency part is more sensitive to shear velocities.A method of extracting formation fast and slow shear slowness and anisotropy parameters from four-component dipole waveforms is proposed,and is verified using numerical data. The error of extracted result is less than 4%.

logging while drilling;acoustic logging;anisotropy;dipole source

P 631.8

A

1673-5005(2014)04-0065-09

10.3969/j.issn.1673-5005.2014.04.009

2013-06-05

國家自然科學(xué)基金項目(11204380,11374371,11134011,61102102);國家油氣重大科技專項(2011ZX05020-009);中國石油科技創(chuàng)新基金(2013D-5006-0304);中國石油天然氣集團公司項目(2011A-3903,2011B-4001)

王瑞甲(1986-),男,博士,主要從事地球物理測井方法、檢測聲學(xué)等方面的研究。E-mail:wruijia@foxmail.com。

喬文孝(1956-),男,教授,博士,博士生導(dǎo)師,研究方向為聲波測井、檢測聲學(xué)。E-mail:qiaowx@cup.edu.cn。

王瑞甲,喬文孝,車小花,等.TI地層聲學(xué)各向異性隨鉆評價方法的數(shù)值研究[J].中國石油大學(xué)學(xué)報:自然科學(xué)版,2014,38(4):65-73.

WANG Ruijia,QIAO Wenxiao,CHE Xiaohua,et al.Numerical study of anisotropy evaluation method by acoustic logging while drilling in TI formations[J].Journal of China University of Petroleum(Edition of Natural Science),2014,38(4):65-73.