楊玉峰, 關(guān)威, 崔乃剛, 胡恒山, 鄭曉波

哈爾濱工業(yè)大學(xué)航天學(xué)院， 哈爾濱　150001

?

隨鉆聲波測井FDTD模擬及鉆鋌波傳播特性研究

楊玉峰, 關(guān)威*, 崔乃剛, 胡恒山, 鄭曉波

哈爾濱工業(yè)大學(xué)航天學(xué)院， 哈爾濱150001

摘要幅度大的鉆鋌波掩蓋地層信號是困擾隨鉆聲波測井技術(shù)的難題，認識鉆鋌波特性對于消除或有效降低鉆鋌波至關(guān)重要.記錄巖石動電效應(yīng)引起電磁場的隨鉆動電測井，被認為有望徹底解決鉆鋌波干擾問題.本文采用有限差分法，模擬不同鉆鋌、不同井孔結(jié)構(gòu)的隨鉆聲場，對比分析了單極源鉆鋌波的傳播特性，闡明存在伴隨鉆鋌波動電信號的原因.結(jié)果表明：鉆鋌聲波在沿鉆鋌傳播時向外部介質(zhì)輻射能量，透過井壁進入地層的具有視鉆鋌波速度的聲波，與地層縱橫波一樣可發(fā)生動電轉(zhuǎn)化，因而可導(dǎo)致隨鉆動電測井時產(chǎn)生鉆鋌波速度的電磁信號.計算還表明，高頻情況下，隨鉆聲波測井鉆鋌波呈現(xiàn)兩階模式：具有低頻截止頻率的高階鉆鋌波幅度較小，其速度略低于鉆鋌縱波速度；無截止頻率的低階鉆鋌波幅度較大，其波速在高頻時甚至低于鉆鋌橫波速度.這種頻散特性和多階模式特性，是徑向多分層開波導(dǎo)結(jié)構(gòu)的導(dǎo)波屬性.

關(guān)鍵詞隨鉆聲波測井； 鉆鋌波； 時域有限差分； 聲誘導(dǎo)電磁場

1引言

隨鉆測井是近年來快速發(fā)展的先進測井技術(shù).與傳統(tǒng)的電纜測井方法相比，隨鉆測井是在鉆井過程中實現(xiàn)對地層參數(shù)的探測，并利用這些參數(shù)及時引導(dǎo)和修正鉆井方向，使井的軌道盡量穿行于儲集層，增加油氣開采效率.對大斜度井和水平井而言，很難進行依靠儀器重力牽引的電纜測井，只能隨鉆測井.如今，美國等發(fā)達國家已廣泛采用基于三大測井原理(電、聲和放射性)的隨鉆測井技術(shù).然而，即使目前最先進的隨鉆聲波測井儀仍有待改進.其中最嚴重的問題是：由于厚壁鉆鋌的存在而產(chǎn)生的、被認為從發(fā)射器經(jīng)鉆鋌傳播到鉆鋌表面接收器的聲波(被稱為鉆鋌波)幅度通常大于由地層返回井內(nèi)的聲波信號，導(dǎo)致地層縱、橫波等被掩蓋不易識別，很難準(zhǔn)確地提取地層縱、橫波速度(王華等，2009).

在發(fā)射器和接收器之間的鉆鋌上周期性刻槽，可以使鉆鋌波逐漸衰減，突出地層聲波信號(蘇遠大等，2011).但是，這種方式并不能完全消除鉆鋌波，而且刻槽占用大量鉆鋌空間，降低鉆鋌的強度和剛度，導(dǎo)致應(yīng)力集中，影響鉆鋌的力學(xué)性能.最近，Zhu等(2012)提出利用隨鉆動電測井解決鉆鋌波干擾問題的設(shè)想：鑒于動電效應(yīng)是含流體孔隙巖石所特有的現(xiàn)象，隨鉆動電測井記錄的電磁信號，不會是鉆鋌中聲波引起的，只能是井外孔隙地層中彈性波導(dǎo)致的.如能從這些電磁信號中提取地層縱波和橫波速度，就無需在鉆鋌上刻槽.在Zhu等(2012)實驗記錄的小尺寸模型井隨鉆動電信號中未發(fā)現(xiàn)鉆鋌波波群，并可以提取出縱、橫波速度.然而，Guan等(2013)、鄭曉波等(2014)近期對典型砂巖地層的理論模擬結(jié)果卻顯示，隨鉆動電測井電場全波中依然存在明顯的鉆鋌波波群，盡管與隨鉆聲波測井相比，它相對其他波群的幅度顯著降低.為了闡明隨鉆動電測井的電磁信號中鉆鋌波存在的原因以及隨鉆動電測井是否可用于解決鉆鋌波干擾問題，本文采用時域有限差分(Finite-Difference Time-Domain, FDTD)方法，在鉆鋌不刻槽的前提下模擬隨鉆聲波測井，通過對比不同模型下鉆鋌外表面和井外地層中的波場情況，分析鉆鋌波的傳播特性.

國內(nèi)外很多學(xué)者已對隨鉆聲波測井問題進行了理論模擬研究.Tang等(2002)和崔志文(2004)分別基于彈性固體和孔隙地層模型，理論模擬了多極源隨鉆聲波測井全波.Sinha等(2009)研究了隨鉆聲波測井全波中各分波的頻散和衰減特性.李希強等(2013)和許松等(2014)分別計算了橫觀各向同性地層和裂隙、孔隙并存地層中的多極源隨鉆聲波測井波場.但對于地層存在水平分層或鉆鋌刻槽等復(fù)雜模型，上述解析算法不再適用，必須采用有限差分、有限元等數(shù)值算法.在隨鉆聲波測井的FDTD模擬方面，Wang和Tang(2003)計算了偶極和四極隨鉆聲波測井全波響應(yīng).王華等(2009)考察了聲源頻率對軟地層多極隨鉆聲波測井響應(yīng)的影響.蘇遠大等(2011)分析了刻槽隔聲技術(shù)的應(yīng)用對降低鉆鋌波幅度的效果.

然而，目前關(guān)于隨鉆聲波測井的理論模擬僅限于計算井內(nèi)波場，這無法回答隨鉆動電測井信號中存在伴隨鉆鋌波電磁場的原因.本文通過FDTD計算鉆鋌外表面接收器、井壁處乃至井外地層中的波場情況，認識鉆鋌波的傳播特性，闡明隨鉆動電測井信號中存在鉆鋌波的原因.這不僅對于隨鉆動電測井方法的后續(xù)研究具有重要意義，也有助于更好地降低鉆鋌波的干擾，改進隨鉆聲波測井技術(shù).本文首先簡要給出軸對稱柱坐標(biāo)系下速度-應(yīng)力交錯網(wǎng)格的彈性波差分計算式，并將FDTD模擬的不同地層隨鉆聲波測井全波與半解析的實軸積分法計算結(jié)果進行了對比，以驗證算法的正確性.然后，分別計算鉆鋌置于無限大流體介質(zhì)、硬地層和軟地層隨鉆聲波測井以及鉆鋌截斷的隨鉆聲波測井等模型波場，著重分析鉆鋌波的傳播特性.最后，給出本文結(jié)論.

2FDTD算法實現(xiàn)與驗證

隨鉆聲波測井模型可簡化為圖1所示的柱面分層結(jié)構(gòu)，沿徑向從內(nèi)向外依次為鉆鋌內(nèi)流體、鉆鋌、鉆鋌外流體和井外地層，其中鉆鋌和地層在本文中均看作彈性固體介質(zhì).

建立z軸與井軸重合的柱坐標(biāo)系(r,θ,z)，并采用圖2所示的速度-應(yīng)力交錯FDTD網(wǎng)格(Randall等，1991)離散模型空間.本文僅以置于鉆鋌外表面的單極源情況為例，研究鉆鋌波的傳播特性.這種關(guān)于井軸對稱的固體彈性波場可由4個應(yīng)力分量和2個速度分量描述，其中切應(yīng)力τrz位于網(wǎng)格的四個角點，正應(yīng)力τr r，τz z和τθ θ位于網(wǎng)格中心，徑向速度vr和軸向速度vz分別位于網(wǎng)格軸向與徑向邊線的中點.

圖1　隨鉆聲波測井示意圖

圖2　軸對稱柱坐標(biāo)系下的速度-應(yīng)力交錯有限差分網(wǎng)格

在軸對稱柱坐標(biāo)系下，上述6個分量滿足固體彈性動力學(xué)偏微分方程組，以其中平衡微分方程的徑向分量為例，表示如下：

(1)

其中ρ為固體介質(zhì)的密度.

由于模擬計算僅關(guān)心軸向和徑向幾米范圍內(nèi)的波場(稱為計算區(qū)域)，相比之下，模型沿軸向和徑向可視作無窮大，必須在計算區(qū)域外設(shè)置吸收邊界.本文采用目前公認吸收效果最好的完全匹配層(Perfectly Matched Layer，PML)(Berenger，1994)吸收邊界條件，并采用不分裂場量(Chew et al.，1997；Wang and Tang，2003)的處理方式.相比于傳統(tǒng)的PML，在設(shè)置不分裂場量的PML時，無需將每一個應(yīng)力和速度分量分裂為幾個變量，而是引入擴展坐標(biāo)系，并近似計算出現(xiàn)的卷積積分項.這種方式不僅易于理解，更重要的是，使得PML區(qū)域和計算區(qū)域具有相同的場量FDTD離散化形式，因此可降低算法實現(xiàn)的復(fù)雜性，提高計算效率.以式(1)為例，其顯示差分表達式可寫成

(2)

(3)

(4)

(5)

(6)

上述針對鉆鋌和井外地層的彈性固體差分表達式同樣適用于鉆鋌內(nèi)、外的流體介質(zhì)，但后者中的切應(yīng)力等于零，無需計算，三個正應(yīng)力分量相等(只需計算一個即可)，且等于流體壓力的負值.此外，對于模型中固體/固體、流體/固體等交界面上場量的計算，仍采用均勻介質(zhì)中的差分表達式，只是其中的介質(zhì)參數(shù)取兩邊介質(zhì)參數(shù)平均值(關(guān)威，2009)的方式，其中剪切模量取調(diào)和平均值、其他參數(shù)取算術(shù)平均值.因此，本文算法在包括固體和流體介質(zhì)的計算區(qū)域以及PML區(qū)域采用統(tǒng)一的差分表達式，這便于下文計算多種模型的波場情況.

圖3和圖4分別為采用上述FDTD算法模擬的硬地層(地層橫波速度大于井內(nèi)流體聲速)和軟地層(地層橫波速度小于井內(nèi)流體聲速)中單極源隨鉆聲波測井的全波波形.計算所需的介質(zhì)參數(shù)與Tang等(2002)和崔志文(2004)相同，可見表1.假設(shè)聲源為位于鉆鋌外表面的環(huán)狀聲壓源，源函數(shù)采用余弦包絡(luò)脈沖，其表達式可見文獻(關(guān)威, 2009)，如無特別說明，聲源中心頻率為f0=8.0 kHz.圖3a和圖4a為源距3.0～4.0 m、間距0.2 m六個接收器的歸一化波形，圖3b和圖4b單獨顯示了源距4.0 m處的波形.為驗證本文算法的正確性，圖中還給出了由實軸積分法(Real Axis Integration, RAI)(Tsang and Rader，1979)計算的波形，其中灰色實線和黑色點劃線分別為FDTD和RAI的模擬結(jié)果.可以看到，無論硬地層還是軟地層，兩種方法計算的不同源距全波中，各個波群的幅度和相位均高度吻合.不僅如此，更詳細的對比結(jié)果顯示，在整個聲波測井源頻率(1～20 kHz)范圍內(nèi)，兩種算法模擬的隨鉆聲波測井和電纜聲波測井全波均吻合的很好.這表明，本文FDTD算法和程序的正確性.

表1　隨鉆聲波測井的介質(zhì)參數(shù)

圖3的硬地層隨鉆聲波測井全波中，根據(jù)波群到時以及時間慢度(STC)處理，最早到達波群(a-a)的速度約為4940 m·s-1.由于其波速高于地層縱波速度(3970 m·s-1)，且在鉆鋌縱波和橫波速度之間，可知它不是地層縱波，而是鉆鋌波.后續(xù)的兩個波群依次為地層橫波和偽瑞利波(b-b)以及斯通利波(c-c).而幅度很小的地層縱波被鉆鋌波掩蓋，無法看到.圖4的軟地層隨鉆聲波測井全波中，最早到達的波群為鉆鋌波(a-a)，其波速為4940 m·s-1，與圖3硬地層中的鉆鋌波波速對比，二者相差不大，說明地層特性基本不影響單極源鉆鋌波波速.隨后的三個波群依次為地層縱波(b-b)、內(nèi)斯通利波(c-c)和外斯通利波(d-d)，與崔志文(2004)的計算結(jié)果一致.與硬地層情況相比，軟地層中不存在地層橫波，且由于縱波速度明顯低于鉆鋌波速度，因而可見地層縱波波群.所謂的內(nèi)斯通利波是沿鉆鋌和鉆鋌內(nèi)流體的鉆鋌內(nèi)表面?zhèn)鞑サ乃雇ɡ?，為了區(qū)別于沿井壁傳播的斯通利波，二者分別被稱為內(nèi)、外斯通利波(崔志文，2004).從幅度上看，圖3硬地層中的鉆鋌波幅度略大于地層橫波幅度，而小于斯通利波幅度.而圖4軟地層中的鉆鋌波幅度明顯大于地層縱波和內(nèi)斯通利波，略小于外斯通利波幅度.需要說明的是，本文假設(shè)井內(nèi)為理想流體，地層為彈性固體，且模擬時沒有通過設(shè)置Q值等方式引入衰減，這與井孔內(nèi)充滿泥漿和井外為孔隙地層的實際情況相比，地層縱、橫波以及斯通利波的幅度明顯偏大，而實際隨鉆聲波測井的鉆鋌波對地層波場的干擾更顯著.

圖3　硬地層中FDTD和RAI模擬的隨鉆聲波測井全波對比圖(a) 源距3.0～4.0 m的歸一化波形； (b)源距4.0 m的波形.

圖4　軟地層中FDTD和RAI模擬的隨鉆聲波測井全波對比圖(a) 源距3.0～4.0 m的歸一化波形； (b) 源距4.0 m的波形.

3鉆鋌波傳播特性分析

目前對隨鉆聲波測井鉆鋌波的普遍認識是(崔志文，2004；鄭曉波等，2014)：它是由聲源激發(fā)，沿鉆鋌傳播的幅度較強的鉆鋌儀器波；鉆鋌波是一個無截止頻率的頻散導(dǎo)波，它的波速低于鉆鋌縱波速度，且隨頻率的增大而降低，單極源鉆鋌波波速受地層參數(shù)的影響很小.

既然前人(Guan et al., 2013; 鄭曉波等, 2014)理論模擬的隨鉆動電測井全波中存在顯著的伴隨鉆鋌波電磁場，說明地層中一定存在以鉆鋌波速度傳播的聲場，進而誘導(dǎo)電磁場.然而，基于目前對鉆鋌波的認識還無法解釋鉆鋌波是如何進入井外地層的.因此，與已有研究僅限于分析井內(nèi)接收器記錄的隨鉆聲波測井波場不同，下文將通過綜合考察不同模型下鉆鋌外流體、井壁處乃至井外地層中的鉆鋌波情況，定量分析地層中的鉆鋌波幅度，闡明鉆鋌波的傳播特性.如無特別說明，本節(jié)所有算例的地層均采用表1中的硬地層.

3.1地層中的鉆鋌波幅度

采用與圖3算例相同的介質(zhì)模型(硬地層)和聲源中心頻率(f0=8.0 kHz)，圖5計算了隨鉆聲波測井在井外地層中的歸一化波場.其中，圖5b—5d分別是距離井壁外0.123 m、0.423 m和1.173 m處，不同軸向位置的地層波場.為了對比，圖5a給出了鉆鋌外表面處的波場.通過對上文圖3的分析可知，圖5a鉆鋌外表面全波中的三個波群依次為鉆鋌波(C)、橫波和偽瑞利波(S)以及斯通利波(ST).對比圖5a和圖5b—5d發(fā)現(xiàn)，后者中的三條曲線斜率與前者中的相同.因此，圖5b—5d中的第一個波群與鉆鋌波(C)有關(guān)，而后兩個波群分別是地層橫波(S)和斯通利波(ST).地層中存在與鉆鋌波有關(guān)(以視鉆鋌波速度沿軸向傳播)的波群，說明鉆鋌波在沿鉆鋌傳播的同時，向鉆鋌外流體乃至井外地層中透射能量.事實上，這也不難理解.與之相類似的問題，比如利用套管波幅度進行固井質(zhì)量評價，正是利用水泥膠結(jié)狀況影響套管波能量進入地層的原理.對比圖5b—5d可以發(fā)現(xiàn)，隨著遠離井壁距離的增大，鉆鋌波和地層橫波波群更加清晰，相反斯通利波波群卻逐漸模糊.由于是歸一化波形，這并不代表鉆鋌波和地層橫波的絕對幅度變大，而是它們相對斯通利波的幅度變大.因為斯通利波能量主要集中在井內(nèi)和靠近井壁附近的地層中，相比于鉆鋌波和地層橫波，斯通利波幅度下降的更顯著.此外，圖5b—5d中鉆鋌波的絕對幅度均大于地層橫波，這說明相比于由聲源傳播到井壁、并透射到井外產(chǎn)生的地層橫波，更顯著的波場能量是沿鉆鋌傳播、同時透射到地層的波.因此，在鉆鋌和聲源換能器之間設(shè)置隔震材料有助于降低鉆鋌波的激發(fā).最后，在離井壁最遠的圖5d中，到時位于鉆鋌波和地層橫波之間的地層縱波被凸顯出來，這說明相比于其他三個波群，地層縱波傳播的更遠.

圖5　隨鉆聲波測井環(huán)境下，遠離井壁不同距離的地層波場(a) 鉆鋌外表面; (b) 井壁外0.123 m處; (c) 井壁外0.423 m處; (d) 井壁外1.173 m處.

圖6　比較三種模型(硬地層和軟地層隨鉆聲波測井，鉆鋌置于無限大流體)的鉆鋌波(a) 在鉆鋌外表面處; (b) 徑向距離為0.117 m的井壁位置處.

通過上述分析可以確定，鉆鋌波能量能夠到達井壁，進入地層，并因此誘導(dǎo)電磁場.說明理論上的隨鉆動電測井波形中應(yīng)該存在鉆鋌波波群.然而，既然前人(Zhu et al.，2012)的隨鉆動電測井實驗中未記錄到鉆鋌波波群，這是否意味著到達井外地層中的鉆鋌波能量相對很少，以至于伴隨鉆鋌波電場的幅度足夠小，反而被伴隨地層縱橫波的電場掩蓋呢？因此，有必要考察井壁處波場中的鉆鋌波幅度，并分析井壁界面的存在對鉆鋌波能量輻射的影響.

圖6對比了硬地層和軟地層隨鉆聲波測井以及鉆鋌置于無限大流體的鉆鋌波.其中圖6a和圖6b分別給出在鉆鋌外表面處和徑向距離為0.117 m的井壁所在位置處，軸向源距為4.0 m的時域波形.為了更清楚地顯示鉆鋌波的幅度和相位差異，圖中僅給出了前2.5 ms的鉆鋌波和地層縱波部分.可以看到，無論是在鉆鋌外表面還是在井壁處，鉆鋌置于無限大流體中的鉆鋌波幅度均明顯不同于硬地層或軟地層隨鉆聲波測井模型的情況.這說明，井壁界面的存在和井外介質(zhì)的變化對鉆鋌波幅度的影響很大.仔細觀察后發(fā)現(xiàn)：在井壁處，隨鉆聲波測井與鉆鋌置于無限大流體的鉆鋌波相位基本相同，而前者的鉆鋌波幅度明顯大于后者，另外硬地層和軟地層的鉆鋌波幅度基本相同；然而，在鉆鋌表面處，隨鉆聲波測井與鉆鋌置于無限大流體的鉆鋌波相位差明顯，而前者的鉆鋌波幅度明顯小于后者，另外硬地層和軟地層的鉆鋌波相位差不大，而硬地層鉆鋌波幅度小于軟地層情況.排除從井孔和井外地層返回井內(nèi)的地層縱波與鉆鋌波疊加的因素(由于地層縱波幅度非常小，此疊加效果不顯著)，出現(xiàn)上述現(xiàn)象的主要原因是遇井壁反射的鉆鋌波能量與原有波場的疊加作用.由于反射波場與原有波場在井壁處無相位差，二者疊加后的幅度大于原波場幅度.而鉆鋌表面處，由于二者存在相位差，以至于疊加后的幅度小于原波場幅度.

圖7　不同井徑情況下，鉆鋌波最大幅值沿徑向距離的變化趨勢

經(jīng)過分析，上述疊加作用應(yīng)受井孔半徑影響，也就是說井徑的變化將改變在鉆鋌外表面處接收的鉆鋌波幅度.圖7考察了保持鉆鋌尺寸不變而改變井徑的情況下，軸向源距4.0 m的鉆鋌波幅度隨徑向位置的變化趨勢.圖中橫坐標(biāo)表示接收位置遠離鉆鋌外表面的徑向距離，等于零代表在鉆鋌外表面處接收，縱坐標(biāo)表示鉆鋌波波群的最大幅值.R1、R2、R3、R4和R5分別表示井徑0.105、0.111、0.117、0.129 m以及0.141 m；Fluid表示鉆鋌置于無限大流體的情況.可以看到，隨著井徑的增大，鉆鋌波幅度明顯減小.無論井徑如何變化，井壁處的鉆鋌波幅度均為最大，且明顯大于鉆鋌外表面處的鉆鋌波幅度，而井外地層中鉆鋌波幅度隨徑向距離的增大逐漸減小.對于井徑較小的情況(r=0.105 m和r=0.111 m)，由于鉆鋌外流體層厚度(鉆鋌外表面與井壁之間流體厚度)除以井內(nèi)流體聲速小于鉆鋌波波群相鄰波峰到波谷的時間，鉆鋌外流體層中的鉆鋌波幅度隨徑向距離的減小而單調(diào)遞減.否則，對于井徑較大的情況(r=0.129 m和r=0.141 m)，隨徑向距離的減小，鉆鋌外流體層中的鉆鋌波幅度先減小后增大.那么，實際隨鉆聲波測井儀上置于鉆鋌外表面接收器記錄的鉆鋌波信號不完全是由聲源沿鉆鋌直達接收器的，還包含向外輻射的鉆鋌波能量遇井壁反射回來的部分.由于上述疊加作用，井徑r=0.117 m時鉆鋌外表面處的鉆鋌波幅度(約為0.02 kPa)反而小于井徑r=0.141 m時鉆鋌外表面處的鉆鋌波幅度(約為0.03 kPa)，而且前者鉆鋌波在全波中覆蓋的時間范圍還要小于后者的.也就是說，在鉆鋌尺寸不變的情況下，采用較大的鉆井半徑不僅增加鉆井成本，還有可能增大鉆鋌波幅度，加大鉆鋌波的干擾問題. 因此，實際測井時，應(yīng)根據(jù)井內(nèi)泥漿聲速，選擇合適的鉆井半徑.

由于井壁處的鉆鋌波幅度明顯大于在鉆鋌外表面處，而且從圖6可以看出，對于本文所用介質(zhì)參數(shù)，井壁處鉆鋌波相對于地層縱波的幅度也大于在鉆鋌外表面處的，因此Guan等(2013)和鄭曉波等(2014)模擬的隨鉆動電測井全波中存在明顯的鉆鋌波波群是合理的.不過，由于不同波群動電轉(zhuǎn)換效率的不同，隨鉆動電測井伴隨鉆鋌波電磁場相對于地層縱、橫波和斯通利波的幅度是明顯小于隨鉆聲波測井鉆鋌波的.對于前人(Zhu et al.，2012)隨鉆動電測井實驗中未觀測鉆鋌波的原因，可能是由于該實驗使用了比較細的鉆鋌，鉆鋌外徑與井孔半徑相比較小，而且聲源探出鉆鋌一段距離以及采用了較高的聲源頻率所致.最后需要指出的是，由于由鉆鋌沿徑向輻射的鉆鋌波能量與隨鉆動電測井的伴隨鉆鋌波電磁場幅度有關(guān)，如果在聲源和接收器之間的鉆鋌外表面包裹一層隔聲材料，可以進一步降低隨鉆動電測井的鉆鋌波幅度.

3.2鉆鋌波導(dǎo)結(jié)構(gòu)對鉆鋌波傳播的影響

為了更深入地分析鉆鋌波的傳播及其向外輻射能量的特性，我們設(shè)計了如圖8所示的鉆鋌在末端斷開的假想模型，計算鉆鋌外表面和地層中的波場.通過考察鉆鋌波在鉆鋌斷開前后區(qū)域的變化，分析鉆鋌波傳播和能量輻射特性以及鉆鋌波導(dǎo)結(jié)構(gòu)對它的影響.從圖8a的模型三維示意圖和圖8b的平面FDTD模擬區(qū)域可以看到，相比于圖3算例的硬地層隨鉆聲波測井模型，本模型中鉆鋌末端被截斷.令h表示鉆鋌截斷位置到聲源的軸向距離，本算例中取h=3.6 m.

圖9給出了在鉆鋌外表面處，軸向源距0.6～4.2 m的歸一化波形.為了波場分析的需要，圖10給出了對圖9波形STC處理的時間慢度圖，其中圖10a和圖10b分別是對源距1.8～2.88 m上行波場和3.72～4.2 m波場的STC圖.圖10以及下文STC圖中的白色水平實線Sp、Ss和Sf分別表示鉆鋌縱波、橫波慢度和井內(nèi)流體聲波慢度.圖10a中僅處理上行波場的原因是：為了能更清晰地分辨上行和下行波場的特性，我們使用STC分別給出了上行和下行波，下行波處理結(jié)果未在該圖中給出.

從圖9中可以看到三個由下至上傳播的波群.根據(jù)圖10a的時間慢度圖可知它們依次為速度略低于鉆鋌縱波速度而高于地層縱波速度的鉆鋌波(C)、地層橫波(S)以及速度略低于井內(nèi)流體聲速的斯通利波(ST).事實上，由于本算例在鉆鋌截斷之前的部分與圖3模型相同，通過上文對圖3中波群的分析，也可確定它們依次為鉆鋌波、橫波和斯通利波.當(dāng)這三個波群傳播到位于3.6 m的分界面(鉆鋌截斷的位置)以后，它們的傳播特性各不相同.可以看到，鉆鋌波和斯通利波這類柱面導(dǎo)波到達分界面后發(fā)生反射和透射.依據(jù)圖9中反射波曲線的斜率和圖10b的透射波場時間慢度圖可知，最先到達分界面的鉆鋌波產(chǎn)生的兩個反射波群分別是鉆鋌波(C-rC)和斯通利波(C-rST)，一個透射波群是斯通利波(C-tST).這是因為鉆鋌波能量不僅限于鉆鋌內(nèi)，而是輻射到鉆鋌外流體乃至地層中，在這種由鉆鋌和井孔組成的兩層波導(dǎo)結(jié)構(gòu)中傳播的鉆鋌波遇到界面后將分別產(chǎn)生與這兩個波導(dǎo)結(jié)構(gòu)有關(guān)的反射、透射導(dǎo)波，即鉆鋌波和斯通利波.由于分界面后鉆鋌波導(dǎo)結(jié)構(gòu)的消失，因此透射波場中沒有鉆鋌波.可以看到，圖中鉆鋌波引起的透射斯通利波，其相對幅度明顯強于分界面以前的入射鉆鋌波幅度.這是因為，原本在鉆鋌內(nèi)傳播的、未被鉆鋌外表面接收器記錄的鉆鋌波能量，在鉆鋌截斷后，也將轉(zhuǎn)化為斯通利波，并被后續(xù)的接收器記錄.由于井孔內(nèi)的斯通利波很難進入鉆鋌內(nèi)，因此只看到入射的斯通利波在此界面處產(chǎn)生一個反射(ST-rST)和一個透射斯通利波(ST-tST).相比于鉆鋌波和斯通利波，井孔內(nèi)的界面對沿井壁傳播的臨界折射橫波沒有影響，地層橫波可以自由地通過該界面所在位置并繼續(xù)向前傳播.

圖8　鉆鋌在末端斷開的模型示意圖(a) 三維模型； (b) FDTD平面模擬區(qū)域.

圖9　鉆鋌截斷模型，在鉆鋌外表面處接收的歸一化波形

圖10　STC處理圖9中波場的時間慢度圖(a) 1.8～2.88 m上行波場的時間慢度圖； (b) 3.72～4.2 m波場的時間慢度圖.

鉆鋌截斷的假想模型可以看作是在鉆鋌上刻一個深凹槽，只是在不考慮滿足鉆鋌力學(xué)要求的實際情況下，這個凹槽的深度達到理想極限情況，貫通了鉆鋌.通過上述分析可知，截斷鉆鋌后，鉆鋌波消失.那么，較深的凹槽可以顯著降低鉆鋌波幅度.盡管如此，波場透射后引起的斯通利波仍然掩蓋地層縱波.因此，并不僅僅是刻越深的凹槽就能有效降低鉆鋌波引起的干擾，還需要周期性刻盡可能多的槽，使得透射波(不論是鉆鋌波還是斯通利波)在來回反射中衰減掉.不過，如果沒有達到有效衰減，來回反射和透射的波場甚至?xí)谏w后面到達的橫波和斯通利波.

接下來，我們考察鉆鋌截斷情況下的地層波場情況.圖11是到井壁徑向距離為0.723 m的地層中軸向源距0.6～4.2 m的歸一化波形.圖12是采用STC法處理圖11中波場的時間慢度圖，其中圖12a和12b分別是對1.8～2.88 m處和3.72～4.2 m處上行波場處理結(jié)果.根據(jù)上文對圖5地層中波場的分析以及圖12a可知，圖11中的三個上行波群分別為鉆鋌波(C)、地層橫波(S)和斯通利波(ST)，其中鉆鋌波幅度最強，斯通利波幅度最弱.地層縱波受到鉆鋌波的掩蓋，無法看到.受鉆鋌截斷的影響，鉆鋌截斷位置(圖12中3.6 m處)附近的波場情況很復(fù)雜.考慮到這個界面并不在地層中，所以這些復(fù)雜的波場不是地層波場遇界面的反射和透射，而是與傳播到鉆鋌截斷界面的鉆鋌和井孔導(dǎo)波(鉆鋌波和斯通利波)有關(guān).根據(jù)圖11中波場的到時，在1.0～2.0 ms之間，界面附近的波場是到達界面的鉆鋌波引起的；而在3.0～4.0 ms之間，界面附近的波場是到達界面的斯通利波引起的.通過上述對圖9鉆鋌表面波場的分析，傳播到界面的鉆鋌波和斯通利波將分別產(chǎn)生反射、透射的斯通利波，此外鉆鋌波還將導(dǎo)致反射鉆鋌波(C-rC).然而，由于斯通利波沿徑向衰減很快，圖11中3.0 ms左右、在界面附近明顯可見的波場并不是反射和透射斯通利波向地層中輻射的能量，而是斯通利波在界面處擾動(ST-D)并作為新源向遠處地層的傳播.同樣的，在圖11中1.0 ms左右的波場也含有鉆鋌波在界面擾動(C-D)并以源的形式向地層傳播的結(jié)果.此外，在1.0～2.0 ms左右的下行波場中，還含有反射鉆鋌波向地層中輻射的能量，而上行波場中還可看到不受界面影響的地層臨界折射橫波(S).

3.3地層波場快照

為了顯示鉆鋌波沿鉆桿傳播時向井外地層輻射能量的情況，圖13給出了聲源中心頻率f0=24.0 kHz時，四種模型在1.2 ms時刻的地層波場快照.其中，圖13a—13d分別為鉆鋌在無限大流體中、隨鉆聲波測井、鉆鋌端部截斷的隨鉆聲波測井以及電纜聲波測井的波場快照.采用較高的聲源頻率是為了能夠較清楚地顯示不同波群的傳播情況.

由于隨鉆聲波測井模型的圖13b和圖13c中，波場傳播情況比較復(fù)雜，為了確定波場中的各個波群，并分析它們傳播特性的需要，圖14給出了中心頻率24.0 kHz的隨鉆聲波測井，在鉆鋌外表面處的導(dǎo)波頻散曲線和時域全波波形.當(dāng)聲源中心頻率為8.0 kHz時，結(jié)合圖3的全波波形以及圖14a的頻散曲線可知，全波中存在三種導(dǎo)波，按速度由高到低分別是鉆鋌波、內(nèi)斯通利波和外斯通利波.當(dāng)聲源中心頻率為24.0 kHz時，全波波場中存在兩階鉆鋌波波群：8.0 kHz時已存在的一階鉆鋌波(對應(yīng)圖14b中的b-b波群)沒有截止頻率，在全頻率范圍內(nèi)均被激發(fā)，它的速度隨頻率的增加逐漸降低，24.0 kHz時的速度甚至小于鉆鋌橫波速度；二階鉆鋌波在16.0 kHz后被激發(fā)，它的速度也隨頻率的增加逐漸降低，在24.0 kHz時，它的速度略低于鉆鋌縱波速度，與一階鉆鋌波相比，它的幅度明顯較小(見圖14b中的a-a波群).此外，24.0 kHz時的全波波形中，后續(xù)的三個波群分別為偽瑞利波、內(nèi)斯通利波和外斯通利波(分別見圖14b中的c-c、d-d和e-e波群).

圖11　鉆鋌截斷模型在距井壁徑向距離0.723 m處地層中的歸一化波形

圖12　STC處理圖11中波場的時間慢度圖(a) 1.8～2.88 m上行波場的時間慢度圖； (b) 3.72～4.2 m波場的時間慢度圖.

圖13　聲源中心頻率為24 kHz時，四種模型在1.2 ms時刻的地層波場快照(a) 鉆鋌置于無限大流體介質(zhì); (b) 隨鉆聲波測井; (c) 鉆鋌端部截斷的隨鉆聲波測井; (d) 電纜測井.

圖14　鉆鋌外表面處，中心頻率24.0 kHz單極源隨鉆聲波測井的導(dǎo)波頻散曲線和全波波形(a) 導(dǎo)波頻散曲線； (b) 源距3.0～4.0 m的歸一化全波波形.

通過上述對圖14的分析，再對比圖13a和圖13b可知，鉆鋌置于無限大流體和隨鉆聲波測井這兩種模型中傳播最快的波群為二階鉆鋌波，而且兩種模型的二階鉆鋌波速度接近，說明井外地層參數(shù)對鉆鋌波速度的影響不大.二階鉆鋌波在傳播過程中不斷向地層中透射能量，并分別以地層縱、橫波速度向遠處地層傳播.類似地，后面幅度更大的一階鉆鋌波也不斷地向地層透射能量，并分別以地層縱橫波速度向遠處地層傳播.從圖13c鉆鋌截斷模型的波場可以看到，已到達鉆鋌截斷位置的二階鉆鋌波不再以鉆鋌波速度向前傳播.這是由于，鉆鋌截斷后，不再有維持鉆鋌波存在的鉆鋌波導(dǎo)結(jié)構(gòu).通過對圖9的分析，到達鉆鋌截斷處的鉆鋌波將以斯通利波的形式繼續(xù)沿井孔傳播，但斯通利波向地層中透射的能量有限.圖13c中最上端的圓弧形曲線為二階鉆鋌波在鉆鋌截斷界面激發(fā)的擾動，以源的形式向遠處地層傳播的結(jié)果.

4結(jié)論

本文實現(xiàn)了隨鉆聲波測井的FDTD模擬算法，在不考慮鉆鋌刻槽的前提下，通過計算鉆鋌置于無限大流體、硬地層和軟地層隨鉆聲波測井以及鉆鋌被截斷的隨鉆聲波測井等模型在鉆鋌外表面處和井外介質(zhì)的波場情況，對比分析了單極源鉆鋌波的傳播特性，闡明了隨鉆動電測井波場中存在伴隨鉆鋌波電磁場的原因.得到的主要結(jié)論如下：

(2) 置于鉆鋌外表面處的隨鉆聲波測井接收器記錄的鉆鋌波由沿鉆鋌直達的波場和向外輻射的鉆鋌波被井壁反射的波場疊加而成.井壁處的鉆鋌波幅度大于鉆鋌外表面處.

(3) 占據(jù)井孔大部分空間的完整厚壁鉆鋌波導(dǎo)結(jié)構(gòu)是鉆鋌波產(chǎn)生和維持其存在的原因，鉆鋌截斷后的井內(nèi)和地層中不再有鉆鋌波傳播，但原有鉆鋌波能量轉(zhuǎn)化為斯通利波透射后仍會掩蓋縱橫波等地層波場.

(4) 高頻(24.0 kHz)隨鉆聲波測井波場中存在兩個鉆鋌波波群，低頻時未出現(xiàn)的高階鉆鋌波(具有低頻截止頻率)幅度較小，其速度略低于鉆鋌縱波速度；無截止頻率的低階鉆鋌波幅度較大，其波速隨頻率的增大而減小，甚至可小于鉆鋌橫波速度.

References

Berenger J P. 1994. A perfectly matched layer for the absorption of electromagnetic waves.J.Computat.Phys., 114(2): 185-200.Chew W C, Jin J M, Michielssen E. 1997. Complex coordinate stretching as a generalized absorbing boundary condition.MicrowaveOpt.Technol.Lett., 15(6): 363-368. Cui Z W. 2004. Theoretical and numerical study of modified Biot′s models, acoustoelectric well logging and acoustic logging while drilling excited by multipole acoustic sources[Ph. D. thesis] (in Chinese). Changchun: Jilin University. Guan W. 2009. Simulation of the coupled poroelastic and electromagnetic wavefields in seismoelectric logging [Ph. D. thesis](in Chinese). Harbin: Harbin Institute of Technology.

Guan W, Hu H S, Zheng X B. 2013. Theoretical simulation of the multipole seismoelectric logging while drilling.Geophys.J.Int., 195(2): 1239-1250.

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

Randall C J, Scheibner D J, Wu P T. 1991. Multipole borehole acoustic waveforms: synthetic logs with beds and borehole washouts.Geophysics, 56(11): 1757-1769.

從同一角度進行分析，顯然也可更好地理解國際上關(guān)于教師專業(yè)成長的這樣一項共識，即是對于“在實踐中學(xué)習(xí)”的突出強調(diào)．例如，在論及在職教師的專業(yè)發(fā)展時，由國際數(shù)學(xué)教育委員會組織的專題研究“數(shù)學(xué)教師的專業(yè)教育與發(fā)展”（ICMI Study 15）就采用了“在實踐中和向?qū)嵺`學(xué)習(xí)”（Learning in and from Practice）這樣一個標(biāo)題[3]．

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

Su Y D, Zhang 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.

Tang X M, Wang T, Patterson D. 2002. Multipole acoustic logging-while-drilling. ∥ SEG Technical Program Expanded Abstracts. 364-367.

Tsang L, Rader D. 1979. Numerical evaluation of the transient acoustic waveform due to a point source in a fluid-filled borehole.Geophysics, 44(10): 1706-1720.

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 T, Tang X M. 2003. Finite-difference modeling of elastic wave propagation: a nonsplitting perfectly matched layer approach.Geophysics, 68(5): 1749-1755.

Xu S, Su Y D, Chen X L, et al. 2014. Numerical study on the characteristics of multipole acoustic logging while drilling in cracked porous medium.ChineseJ.Geophys. (in Chinese), 57(6): 1999-2012, doi: 10.6038/cjg20140630.

Zheng X B, Hu H S, Guan W, et al. 2014. Theoretical simulation of the electric field induced by acoustic waves during the seismoelectric logging while drilling.ChineseJ.Geophys. (in Chinese), 57(1): 320-330, doi: 10.6038/cjg20140128.

Zhu Z Y, Wang J, Toks?z M N. 2012. Multipole seismoelectric logging while drilling (LWD) for acoustic velocity measurements. ∥ 74th EAGE Conference & Exhibition Incorporating SPE EUROPEC. Massachusetts Institute of Technology.

附中文參考文獻

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

關(guān)威. 2009. 孔隙介質(zhì)彈性波-電磁場耦合效應(yīng)測井的波場模擬研究[博士論文]. 哈爾濱: 哈爾濱工業(yè)大學(xué).

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

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

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

許松, 蘇遠大, 陳雪蓮等. 2014. 含孔隙、裂隙地層隨鉆多極子聲波測井理論. 地球物理學(xué)報, 57(6): 1999-2012, doi: 10.6038/cjg20140630.

鄭曉波, 胡恒山, 關(guān)威等. 2014. 隨鉆動電測井中聲誘導(dǎo)電場的理論模擬. 地球物理學(xué)報, 57(1): 320-330, doi: 10.6038/cjg20140128.

(本文編輯胡素芳)

基金項目國家自然科學(xué)基金(41204092, 41174110)和中國博士后特別基金(2012T50314)資助.

作者簡介楊玉峰,男,1990年生,哈爾濱工業(yè)大學(xué)固體力學(xué)專業(yè)博士生.主要從事聲波測井的理論和數(shù)值模擬研究.E-mail:xiaoy931@163.com*通訊作者關(guān)威，男，1982年生，哈爾濱工業(yè)大學(xué)固體力學(xué)專業(yè)副教授.主要從事孔隙介質(zhì)中彈性波和震電耦合波傳播及其在地學(xué)中的應(yīng)用研究.E-mail: guanw@hit.edu.cn

doi:10.6038/cjg20160131 中圖分類號P631

收稿日期2014-10-08，2015-09-28收修定稿

FDTD simulation and analysis of the collar wave propagation in acoustic logging while drilling

YANG Yu-Feng, GUAN Wei*, CUI Nai-Gang, HU Heng-Shan, ZHENG Xiao-Bo

SchoolofAstronautics,HarbinInstituteofTechnology,Harbin150001,China

AbstractIt is still an outstanding problem in acoustic logging while drilling (LWD) technique of the large-amplitude collar wave covering the formation waves, thus cognition of the characteristics of the collar wave is essential for eliminating or effectively reducing that wave. The seismoelectric LWD that records the electromagnetic fields induced by acoustic waves in the formation is expected to solve this problem thoroughly. In this study, a finite-difference time-domain algorithm is implemented to simulate the acoustic LWD waves in various collars and borehole models. Based on the propagation characteristics of the collar wave excited by a monopole source, the reason is addressed of the electric field accompanying the collar wave in seismoelectric LWD. It is found that the collar wave radiates energy into the outer media continually when propagating along the collar. Part of the energy with an apparent collar-wave velocity goes into the formation through the borehole wall, which induces electric fields like the formation compressional and shear waves thus can cause a signal of with collar-wave velocity in the seismoelectric LWD. In addition, there are two collar wave groups in the high-frequency acoustic LWD full waveforms. The high-order one having smaller amplitude cuts off in low frequency and it is slightly slower than the compressional wave in the collar. The low-order one having larger amplitude exists in the whole frequency range and it is even slower than the shear wave in the collar. The dispersion and multiple modes are the characteristics of guided waves in cylindrically-layered open wave guide structures.

KeywordsLogging while drilling; Collar wave; Finite-difference time-domain; Acoustic-induced electromagnetic field

楊玉峰, 關(guān)威, 崔乃剛等. 2016. 隨鉆聲波測井FDTD模擬及鉆鋌波傳播特性研究.地球物理學(xué)報，59(1):368-380,doi:10.6038/cjg20160131.

Yang Y F, Guan W, Cui N G, et al. 2016. FDTD simulation and analysis of the collar wave propagation in acoustic logging while drilling.ChineseJ.Geophys. (in Chinese),59(1):368-380,doi:10.6038/cjg20160131.