偶極橫波遠(yuǎn)探測在裸眼井及套管井中接收波場的方位響應(yīng)特征

李楊虎， 唐曉明,3*， 蘇遠(yuǎn)大,3

1 中國石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院， 青島 266580 2 中國石油大學(xué)(華東)深層油氣重點(diǎn)實(shí)驗(yàn)室， 青島 266580 3 青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室海洋礦產(chǎn)資源評(píng)價(jià)與探測技術(shù)功能實(shí)驗(yàn)室， 青島 266071

0 引言

近年來，偶極橫波遠(yuǎn)探測成像測井技術(shù)取得長足發(fā)展，其在裂縫、斷層、溶蝕孔洞等地質(zhì)構(gòu)造的探測中取得了良好的應(yīng)用效果.Tang(2004)于2004年首次提出了偶極橫波遠(yuǎn)探測成像方法，該方法通過充液井中的偶極聲源向地層輻射橫波，并在井中接收經(jīng)反射體反射回來的橫波進(jìn)行井周成像.與常規(guī)的單極縱波法相比(Hornby, 1989; 薛梅, 2002; 喬文孝等, 2004; 陶果等, 2008)，偶極橫波法具有以下優(yōu)勢：一、偶極聲源具有方位靈敏度，采用多分量發(fā)射和接收可以確定反射體的走向；二、偶極聲源的工作頻率為3 kHz左右的低頻，因而其徑向探測深度更深.此后，國內(nèi)外學(xué)者針對(duì)偶極橫波遠(yuǎn)探測技術(shù)開展了大量的研究工作，其現(xiàn)已成為非常規(guī)油氣勘探中尋找油氣甜點(diǎn)的利器(Tang and Patterson, 2009; 唐曉明和魏周拓, 2012; 蘇遠(yuǎn)大等, 2014; Tang et al., 2014, 2016; Hirabayashi et al., 2017; 宮昊等, 2018; Bennett, 2019; Lee et al., 2019; Xu et al., 2019; 董經(jīng)利等, 2020; 古希浩等, 2020; 劉匯鑫等, 2021).

研究井孔對(duì)聲場的調(diào)制作用對(duì)偶極橫波遠(yuǎn)探測具有重要意義，包括井內(nèi)偶極聲源在地層中的輻射特性以及井孔對(duì)入射波的方位接收響應(yīng).對(duì)于偶極聲源在井外的輻射聲場，國內(nèi)學(xué)者已經(jīng)做了充分研究.Tang等(2014)給出了偶極聲源在井外無限大地層中輻射橫波的位移遠(yuǎn)場漸近解及其輻射指向性表達(dá)式.在此基礎(chǔ)上，曹景記等(2014)和唐曉明等(2019)研究了偶極聲源在裸眼和套管井外的橫波輻射特征，并證實(shí)了在套管井中進(jìn)行偶極橫波遠(yuǎn)探測的可行性.井孔對(duì)入射波的接收響應(yīng)可分為平面波和球面波入射兩種情況，對(duì)于平面波入射引起的井內(nèi)聲場的方位響應(yīng)，人們已經(jīng)進(jìn)行了深入研究，見White(1953)、Schoenberg(1986)、Peng等(1993)以及Hirabayashi等(2017)發(fā)表的研究成果.然而，聲波遠(yuǎn)探測中入射波場是以球面波的形式向井孔傳播，對(duì)于球面波入射引起的井內(nèi)聲場的方位響應(yīng)尚未做深入研究.雖然前人能夠?qū)ε紭O橫波遠(yuǎn)探測的接收聲場進(jìn)行快速模擬(Tang et al., 2014; 許家旗和胡恒山, 2019)，但這些研究中都使用了互易定理來求解井孔的球面波入射響應(yīng)，該方法主要用于分析井中聲源和接收器位置同軸線的情況，不能對(duì)實(shí)際測井時(shí)儀器接收到的方位波場進(jìn)行有效模擬.針對(duì)此，Li等(2021)在虛源法的基礎(chǔ)上，結(jié)合柱面波展開法(Tang and Cheng，2004)，求解得到了充液裸眼井中偶極橫波方位遠(yuǎn)探測接收聲場的解析解，為方位聲波遠(yuǎn)探測的模擬奠定了理論基礎(chǔ).

目前，利用偶極橫波遠(yuǎn)探測進(jìn)行井周成像或反射體方位識(shí)別的相關(guān)方法主要是根據(jù)井中接收波形的幅度變化特征(Tang and Patterson, 2009; Hirabayashi et al., 2017; 李盛清等, 2020; 許家旗和胡恒山, 2020; Li et al., 2021)，但充液井中遠(yuǎn)探測波場的方位接收響應(yīng)還缺乏詳盡的理論分析.因此，對(duì)偶極橫波遠(yuǎn)探測聲場在井內(nèi)流體中產(chǎn)生的徑向位移和聲壓的方位響應(yīng)特征進(jìn)行深入研究具有重要意義.另一方面，基于全球各地油田中套管井的大量存在以及被證實(shí)的過套管偶極橫波遠(yuǎn)探測的可行性(唐曉明等, 2019)，本文將裸眼井中遠(yuǎn)探測聲場的方位解析解延展到套管井中，研究了不同膠結(jié)情況套管井中接收波場的方位響應(yīng)特征，并與裸眼井情況進(jìn)行了對(duì)比分析.最后，利用三維有限差分法證實(shí)了聲場幅值方位響應(yīng)中的180°反轉(zhuǎn)現(xiàn)象.本文的研究對(duì)偶極橫波遠(yuǎn)探測測井時(shí)裸眼和套管井中接收聲場的方位響應(yīng)特性有了更為全面的認(rèn)識(shí)，為偶極橫波方位遠(yuǎn)探測的技術(shù)發(fā)展、儀器研制以及現(xiàn)場解釋提供了理論支持.

1 充液井中偶極橫波方位遠(yuǎn)探測接收聲場的解析解

圖1a給出了當(dāng)井外存在反射體時(shí)，充液裸眼或套管井中偶極聲波遠(yuǎn)探測模擬所用的坐標(biāo)系.在測井儀器的偶極聲源上方環(huán)向布置了多組偏心距為e的接收器，其接收橫截面如圖1b中模型所示.按照不同的聲源—接收器組合方式，可模擬不同的測量工作模式，如正交偶極發(fā)射—正交偶極接收模式，正交偶極發(fā)射—多方位接收模式，本文采用多方位接收模式對(duì)井中接收波場的方位響應(yīng)進(jìn)行研究.偶極聲波遠(yuǎn)探測包含三個(gè)過程，即充液井孔內(nèi)偶極聲源在井外的輻射、井外地質(zhì)體的反射以及充液井孔對(duì)反射波的接收調(diào)制.在輻射距離遠(yuǎn)大于波長的遠(yuǎn)場輻射條件下，圖1a中偶極聲源輻射到地層中的SH和SV橫波位移勢函數(shù)和Γ的遠(yuǎn)場漸近解為(唐曉明和魏周拓, 2012)：

(1)

當(dāng)輻射橫波與地層中的反射體作用后，可將反射回井孔的入射橫波波場表示為(Tang et al., 2014, 2016)：

(2)

式中，RD為聲源在井外遠(yuǎn)處的輻射函數(shù)，由式(1)的中括號(hào)里的部分組成；RF為聲波在反射體處的反射系數(shù)，可由佐普利茲方程求得(Aki and Richards,

圖1 基于虛源法的偶極聲波遠(yuǎn)探測示意圖(a)及方位接收井孔截面示意圖(b)Fig.1 The schematic model of dipole single-well reflection imaging based on the virtual source analogy (a) and the cross section of borehole receiving configuration using azimuthally spaced off-axis receivers (b)

1980)；S、RD、RF都是關(guān)于圓頻率ω的函數(shù)；D為聲波在地層中總的傳播距離，其傳播路徑為一條折線.Li等(2021)通過將遠(yuǎn)探測聲場看作是聲源在反射體向外一側(cè)的鏡像點(diǎn)的輻射聲場，此時(shí)聲波的傳播路徑D轉(zhuǎn)化為一條直線，如圖1a中所示，從而式(2)中的球面波傳播因子可以展開為多極子柱面波疊加的形式(唐曉明和鄭傳漢, 2004)：

×eik(z-z0)dk，

(3)

將式(3)代入式(2)，對(duì)于圖1a中所示的偶極聲源指向，入射SH和SV橫波的位移勢函數(shù)表達(dá)式為：

×cos[n(φ-φ0)],

(4)

其中，上標(biāo)i表示入射波場；下標(biāo)SH和SV表示計(jì)算聲源的輻射函數(shù)RD和反射系數(shù)RF時(shí)應(yīng)分別對(duì)應(yīng)SH和SV橫波.

考慮徑向包含多層的井眼模型，假設(shè)井壁與地層之間共有N-1圓柱層，其中第m層為固體層，第j層為液體層.當(dāng)入射橫波與井眼相互作用時(shí)，根據(jù)赫姆霍茲方程，在柱坐標(biāo)系(r,φ,z)下井內(nèi)流體中的位移場uf、第j層流體中的位移場uf(j)和第m層固體中散射波的位移場us(m)分別為：

(5)

(6)

在裸眼井情況下，式(5)中的模型由井內(nèi)流體和地層兩部分組成，此時(shí)井壁處的邊界條件為：徑向位移和正應(yīng)力連續(xù)，環(huán)向和軸向的剪切應(yīng)力為零，表達(dá)式為：

(7)

式中，a為井孔半徑；上標(biāo)s和f分別表示地層中的散射波場和流體中的透射波場.

由式(4)、(5)、(6)計(jì)算出入射波、散射波和透射波的位移和應(yīng)力，代入式(7)的邊界條件得到兩個(gè)矩陣方程：

(8)

為了便于描述，將式(8)中方程左側(cè)與井中透射波和地層中散射波有關(guān)的矩陣記為矩陣M，右側(cè)與地層中入射波有關(guān)的向量分別記為向量b和b′.矩陣M，向量b和b′中的元素表達(dá)式詳見附錄A.

(9)

其中，矩陣元素N11～N41及向量c和c′可由方程(8)中與之相對(duì)應(yīng)的部分將井孔半徑置換為套管內(nèi)徑后得到；G矩陣為6×6階復(fù)數(shù)矩陣，計(jì)算G矩陣所用到的T矩陣的元素表達(dá)式詳見附錄B.

在套管和地層之間添加一流體層來模擬自由套管的情況，此時(shí)由井中流體-套管-流體環(huán)-地層組成的模型中所包含的3個(gè)界面均為流-固界面.與固-固界面不同，對(duì)流-固界面的每一個(gè)邊界條件都要單獨(dú)進(jìn)行處理，通過連接邊界條件，最后整合得到矩陣方程：

(10)

式中，H為12×12矩陣，d和d′為12×1向量；系數(shù)向量中的上標(biāo)1、2和3表示該振幅系數(shù)分別對(duì)應(yīng)套管、流體環(huán)和地層介質(zhì).矩陣H，向量d和d′中的元素表達(dá)式詳見附錄C.

(11)

其中，ρf為流體密度.對(duì)于遠(yuǎn)場輻射，虛源的徑向距離r0遠(yuǎn)大于彈性波的波長，滿足|sr0|?1的條件，由此可對(duì)式(11)中的波數(shù)積分運(yùn)用最速下降方法求解(Tang and Patterson, 2009; Li et al., 2021)，整理得到：

(12)

(13)

式(13)指出，井軸上接收到的流體聲壓只有單極項(xiàng)(n=0)貢獻(xiàn)，而徑向位移只有偶極項(xiàng)(n=1)貢獻(xiàn).

綜上所述，利用式(12)和式(13)便可計(jì)算偶極橫波遠(yuǎn)探測聲場在裸眼和套管井情況下所引起的井內(nèi)流體的徑向位移和聲壓的方位響應(yīng).

2 算法對(duì)比驗(yàn)證及分析

我們首先利用三維有限差分法驗(yàn)證上述解析解的正確性，采用膠結(jié)良好的套管井模型.計(jì)算所用的井孔及地層參數(shù)如表1所示，其中反射體設(shè)為地層1和地層2之間傾角為45°的界面，位于0°方位，虛源的徑向距離為10 m.聲源上方1 m處共有8個(gè)接收器間隔45°環(huán)向布置，接收器的偏心距離均為0.035 m.聲源采用中心頻率為3 kHz的Ricker子波，差分模擬時(shí)的空間網(wǎng)格和迭代時(shí)間步長分別為0.005 m和0.48 μs.

表1 模型計(jì)算參數(shù)Table 1 Model parameters for calculation

通過將偶極聲源指向y(x)來模擬SH(SV)橫波入射，圖2和圖3給出了對(duì)比結(jié)果.其中，虛線為有限差分計(jì)算的結(jié)果，實(shí)線為解析解計(jì)算的結(jié)果.由于有限差分法模擬的是井內(nèi)和井外的全波波場，而解析解式(12)只模擬了井外的反射波部分，因此圖2和圖3中只顯示了反射橫波的到時(shí)區(qū)間.當(dāng)聲源指向y向時(shí)，產(chǎn)生純SH入射波.從圖2中可以看出，無論是井內(nèi)流體的徑向位移(圖2a)還是聲壓(圖2b)，采用解析解得到的漸近結(jié)果與有限差分模擬的數(shù)值結(jié)果幾乎完全重合，這說明解析解是準(zhǔn)確可靠的.圖2中也顯示出了井中接收波形的方位特征：當(dāng)接收器指向與反射體走向平行(即90°和270°)時(shí)，接收到的徑向位移和聲壓的幅度最大；而當(dāng)接收器指向垂直于反射面(即0°和180°)時(shí)，接收不到徑向位移和聲壓；此外，x軸上下兩側(cè)接收到的徑向位移和聲壓的極性均反生反轉(zhuǎn).上述分析表明，SH橫波入射引起的井內(nèi)流體的徑向位移和聲壓響

圖2 SH橫波入射引起的井內(nèi)流體響應(yīng)的漸近解(實(shí)線)和三維有限差分(虛線)的計(jì)算結(jié)果對(duì)比(a) 徑向位移； (b) 聲壓.Fig.2 Comparison of the borehole-fluid response calculated using the analytical solution (solid curves) and the 3D finite-difference method (dashed curves) for SH-wave incidence(a) Radial displacement; (b) Pressure.

應(yīng)對(duì)井外反射體方位的靈敏度高，且兩者特征具有較好的一致性，說明在進(jìn)行遠(yuǎn)探測時(shí)采用反射波中的SH分量更具優(yōu)勢；但該分析也表明了其存在的180°方位不確定性.

當(dāng)聲源指向x向時(shí)，利用解析解式(12)計(jì)算的是純SV橫波入射，但實(shí)際上聲源在產(chǎn)生SV橫波的同時(shí)還會(huì)產(chǎn)生P波入射，因此有限差分結(jié)果中還包含一些小振幅的P-SV和SV-P轉(zhuǎn)換波(如圖3中4.5～5.5 ms所示).盡管如此，圖中顯示兩種方法計(jì)算的SV橫波反射結(jié)果吻合較好.與圖2中SH橫波的入射結(jié)果相比，SV橫波入射引起的井內(nèi)流體的徑向位移和聲壓的方位特征存在較大差異.對(duì)于徑向位移而言(圖3a)，當(dāng)接收器靠近反射體時(shí)接收到的幅度最大，當(dāng)接收器指向平行于反射面(即90°和270°)時(shí)接收到的幅度最小但不為零；y軸左右兩側(cè)接收到的位移的極性反生反轉(zhuǎn).而聲壓(圖3b)則表現(xiàn)出較強(qiáng)的方位性，靠近反射體的接收器接收到的聲壓幅度最大，并且不同方位接收的聲壓極性相同.結(jié)果表明，在進(jìn)行遠(yuǎn)探測時(shí)，考慮反射波中的SV分量有助于識(shí)別井外反射體的方位.

圖3 SV橫波入射引起的井內(nèi)流體響應(yīng)的漸近解(實(shí)線)和三維有限差分(虛線)的計(jì)算結(jié)果對(duì)比(a) 徑向位移； (b) 聲壓.Fig.3 Comparison of the borehole-fluid response calculated using the analytical solution (solid curves) and the 3D finite-difference method (dashed curves) for SV-wave incidence(a) Radial displacement; (b) Pressure.

在套管井復(fù)雜模型下所表現(xiàn)出的良好的吻合性，驗(yàn)證了解析解的正確性和普適性.然而，兩種方法的計(jì)算效率相差很大，為了更準(zhǔn)確地模擬接收器的偏心效應(yīng)以及套管的影響，在三維有限差分建模中必須使用毫米級(jí)的網(wǎng)格尺寸，這需要花費(fèi)數(shù)10 h才能完成模擬.相比之下，式(12)的解析解只需要數(shù)秒便可以得到準(zhǔn)確的結(jié)果，為偶極橫波方位遠(yuǎn)探測提供了一種快速高效的模擬方法.

為了進(jìn)一步詳細(xì)考察圖2和圖3中所示波形的方位響應(yīng)特征，我們?cè)黾恿私邮掌鞯膫€(gè)數(shù)并縮小其環(huán)向布置間隔.考慮到偶極橫波反射成像中存在的180°方位不確定性問題，圖4和圖5的極坐標(biāo)圖分別給出了SH和SV橫波入射時(shí)兩個(gè)相反反射體方位(φ0=0°，180°)的計(jì)算結(jié)果(圖中實(shí)線)，其他參數(shù)與圖2中相同，結(jié)果按不同方位接收波形的最大幅值進(jìn)行歸一化處理.如極坐標(biāo)圖中所示，其結(jié)果在0°～180°和180°～360°之間具有對(duì)稱性，因此直角坐標(biāo)系中只顯示了0°～180°的結(jié)果.為了對(duì)比，同時(shí)模擬了當(dāng)接收器居中時(shí)(即位于井軸上)接收的結(jié)果，如圖中虛線所示.從圖4中可以看出，對(duì)于SH橫波入射，接收器居中時(shí)接收到的流體徑向位移的方位幅值響應(yīng)關(guān)于90°和270°對(duì)稱，并且接收不到聲壓；而當(dāng)接收器偏心時(shí)，接收到的徑向位移和聲壓的幅度最大值所在方位均由90°和270°向靠近反射體一側(cè)偏移.對(duì)于SV橫波入射(圖5)，接收器居中時(shí)接收到的流體徑向位移和聲壓的方位幅值響應(yīng)均關(guān)于90°和270°對(duì)稱，并且聲壓為一常數(shù)值；而接收器偏心時(shí)接收到的徑向位移的幅度最小值所在方位由90°和270°向靠近反射體一側(cè)偏移，并且靠近反射體的接收器將接收到最大幅值的徑向位移和聲壓.接收器居中時(shí)聲場響應(yīng)具有的對(duì)稱現(xiàn)象可從式(13)中得到，這種對(duì)稱性正是我們所熟知的180°不確定性.相比之下，接收器偏心接收數(shù)據(jù)所展示出的非對(duì)稱性對(duì)井外反射體的方位更加敏感，其方位靈敏度有助于解決180°不確定問題.

圖4 對(duì)于兩個(gè)相反方位的反射體，SH橫波入射引起的井內(nèi)波場的方位響應(yīng)(a)、(b) 徑向位移； (c)、(d) 聲壓.Fig.4 The azimuthal response of waveforms inside borehole for SH-wave incidence for the 0° (red) and 180° (black) reflector positions. For comparison, the on-axis result (dashed curves) is also plotted(a) and (b) Radial displacement; (c) and (d) Pressure.

圖5 對(duì)于兩個(gè)相反方位的反射體，SV橫波入射引起的井內(nèi)波場的方位響應(yīng)(a)、(b) 徑向位移； (c)、(d) 聲壓.Fig.5 The azimuthal response of waveforms inside borehole for SV-wave incidence for the 0° (red) and 180° (black) reflector positions. For comparison, the on-axis result (dashed curves) is also plotted(a) and (b) Radial displacement; (c) and (d) Pressure.

接下來我們將重點(diǎn)分析不同接收器偏心距、聲源頻率和地層類型下裸眼及套管井中接收波場的方位響應(yīng)特征.

3 裸眼和套管井中接收波場方位響應(yīng)的影響因素研究

3.1 接收器偏心距對(duì)井中接收波場的影響

圖6和圖7分別給出了裸眼井、自由套管和套管膠結(jié)良好情況下，SH和SV橫波入射時(shí)，不同偏心距接收器接收到的井內(nèi)流體徑向位移和聲壓的方位幅值響應(yīng)曲線.計(jì)算所用聲源的中心頻率為3000 Hz，源距為3 m，反射體傾角為60°，位于0°方位，虛源的徑向距離為30 m，井孔和地層參數(shù)見表1(采用地層1).對(duì)于0.001 m的微小偏心距，三種情況下SH和SV橫波入射引起的井內(nèi)流體的徑向位移響應(yīng)相對(duì)于90°方位幾乎對(duì)稱，與圖4和圖5中接收器居中接收的結(jié)果相似，此時(shí)徑向位移的測量結(jié)果主要以偶極成分為主，這也說明較小的偏心距對(duì)測量結(jié)果的影響可以忽略.隨著偏心距離從0.001 m增加到0.05 m，當(dāng)SH橫波入射時(shí)，裸眼井和自由套管情況下井中接收徑向位移和聲壓的振幅最大值所在方位均逐漸偏離反射體位置，而膠結(jié)良好套管井中呈現(xiàn)出偏離方位相反的現(xiàn)象；當(dāng)SV橫波入射時(shí)，三種情況下井中接收徑向位移的振幅最小值所在方位均逐漸靠近反射體位置，并且其振幅最小值隨著偏心距的增加而逐漸增大.此外，SV橫波入射時(shí)井內(nèi)聲壓的振幅最大值所在方位不受接收器偏心距的影響，但受井類型的影響，在裸眼井和膠結(jié)良好套管井中，靠近反射體一側(cè)的接收器接收到的聲壓幅值最大，而在自由套管中則是遠(yuǎn)離反射體一側(cè)的接收器將接收到最大幅值的聲壓.值得注意的是，雖然加套管之后井中方位接收響應(yīng)的變化規(guī)律不盡相同，但與裸眼井情況相比，自由套管和套管膠結(jié)良好時(shí)井中接收波場的方位響應(yīng)受接收器偏心距的影響更大，說明套管井中的波場響應(yīng)對(duì)井外反射體方位的靈敏度更高.

圖6 快速地層條件下，SH橫波入射時(shí)，不同偏心距接收器接收到的井內(nèi)波場的方位響應(yīng)(a) 裸眼井； (b) 自由套管； (c) 膠結(jié)良好.Fig.6 The azimuthal response of waveforms inside borehole received by receivers with different eccentric distances for SH-wave incidence in the fast formation(a) Open hole; (b) Cased hole with free pipe; (c) Cased hole with good bonding.

圖7 快速地層條件下，SV橫波入射時(shí)，不同偏心距接收器接收到的井內(nèi)波場的方位響應(yīng)(a) 裸眼井； (b) 自由套管； (c) 膠結(jié)良好.Fig.7 The azimuthal response of waveforms inside borehole received by receivers with different eccentric distances for SV-wave incidence in the fast formation(a) Open hole; (b) Cased hole with free pipe; (c) Cased hole with good bonding.

3.2 聲源頻率對(duì)井中接收波場的影響

圖8和圖9是聲源頻率變化時(shí)，SH和SV橫波入射到裸眼井、自由套管和膠結(jié)良好套管井中引起的井內(nèi)流體的徑向位移和聲壓的方位幅值響應(yīng)曲線.計(jì)算時(shí)接收器的偏心距離為0.035 m，其他參數(shù)與圖6中相同.從圖中可看出，對(duì)于100 Hz的低頻情況，SH和SV橫波入射引起的井內(nèi)流體的徑向位移響應(yīng)相對(duì)于90°方位幾乎對(duì)稱.這是因?yàn)轭l率非常低時(shí)，聲波波長遠(yuǎn)大于井徑，使井孔的調(diào)制作用明顯減弱，說明低頻時(shí)可忽略接收器的偏心效應(yīng).隨著頻率從100 Hz增加到5000 Hz，當(dāng)SH橫波入射時(shí)，自由套管和膠結(jié)良好情況下井中接收徑向位移和聲壓的振幅最大值所在方位均先靠近后偏離反射體位置，但在裸眼井情況下則表現(xiàn)出先靠近后偏離再靠近反射體位置的現(xiàn)象.當(dāng)SV橫波入射時(shí)，隨著頻率增加，三種情況下井中接收徑向位移和聲壓的方位響應(yīng)較為復(fù)雜，其幅度最大值所在方位發(fā)生了180°反轉(zhuǎn).并且某些頻率下，聲壓最大值所在方位既不位于靠近反射體一側(cè)，也不位于遠(yuǎn)離反射體一側(cè)(比如裸眼井在500 Hz和1000 Hz時(shí)，以及膠結(jié)良好套管井在500 Hz、1000 Hz和2000 Hz時(shí)).上述結(jié)果表明，雖然聲源頻率能夠顯著改變井內(nèi)接收波場的方位響應(yīng)特征，但方位響應(yīng)隨頻率呈復(fù)雜變化.

圖8 快速地層條件下，聲源頻率不同時(shí)，SH橫波入射引起的井內(nèi)波場的方位響應(yīng)(a) 裸眼井； (b) 自由套管； (c) 膠結(jié)良好.Fig.8 The azimuthal response of waveforms inside borehole at different frequencies for SH-wave incidence in the fast formation(a) Open hole; (b) Cased hole with free pipe; (c) Cased hole with good bonding.

圖9 快速地層條件下，聲源頻率不同時(shí)，SV橫波入射引起的井內(nèi)波場的方位響應(yīng)(a) 裸眼井； (b) 自由套管； (c) 膠結(jié)良好.Fig.9 The azimuthal response of waveforms inside borehole at different frequencies for SV-wave incidence in the fast formation(a) Open hole; (b) Cased hole with free pipe; (c) Cased hole with good bonding.

3.3 地層類型對(duì)井中接收波場的影響

除了接收器偏心距和聲源頻率外，地層類型也會(huì)影響井中接收波場的方位響應(yīng).圖10和圖11對(duì)比了快速和慢速地層條件下，SH和SV橫波入射到裸眼井、自由套管和膠結(jié)良好套管井中引起的井內(nèi)流體的徑向位移和聲壓的方位幅值響應(yīng).計(jì)算時(shí)接收器的偏心距為0.035 m，聲源主頻為3000 Hz，其他參數(shù)與圖6中相同.從圖10中可以看出，當(dāng)SH橫波入射時(shí)，裸眼井和自由套管井中流體徑向位移和聲壓的振幅最大值所在方位在快速地層條件下均由90°向遠(yuǎn)離反射體一側(cè)偏移，而在慢速地層中均向靠近反射體一側(cè)偏移，呈現(xiàn)出相反的現(xiàn)象；在快速和慢速地層條件下，膠結(jié)良好時(shí)的結(jié)果均由90°向靠近反射體一側(cè)偏移.當(dāng)SV橫波入射時(shí)(圖11)，無論在快速還是慢速地層中，三種情況下井中流體徑向位移的振幅最小值所在方位的變化趨勢相同，均由90°向靠近反射體一側(cè)偏移，并且靠近反射體的接收器將接收到最大幅值的徑向位移；相對(duì)于快速地層，在慢速地層中三種情況下的聲壓方位響應(yīng)均反生180°反轉(zhuǎn).上述分析表明，當(dāng)根據(jù)井內(nèi)接收波場的幅值信息對(duì)井外反射體進(jìn)行方位識(shí)別時(shí)，需要考慮地層性質(zhì)的影響.

圖10 快速和慢速地層條件下，SH橫波入射引起的井內(nèi)波場的方位響應(yīng)(a) 裸眼井； (b) 自由套管； (c) 膠結(jié)良好.Fig.10 The azimuthal response of waveforms inside borehole for SH-wave incidence in the fast and slow formations(a) Open hole; (b) Cased hole with free pipe; (c) Cased hole with good bonding.

圖11 快速和慢速地層條件下，SV橫波入射引起的井內(nèi)波場的方位響應(yīng)(a) 裸眼井； (b) 自由套管； (c) 膠結(jié)良好.Fig.11 The azimuthal response of waveforms inside borehole for SV-wave incidence in the fast and slow formations(a) Open hole; (b) Cased hole with free pipe; (c) Cased hole with good bonding.

4 井中聲場方位幅值響應(yīng)反轉(zhuǎn)現(xiàn)象的正確性驗(yàn)證

上一節(jié)研究了偶極橫波遠(yuǎn)探測在裸眼和套管井中接收波場的方位響應(yīng)特征，我們注意到對(duì)于不同的模型參數(shù)，井中接收波場的方位幅值特征存在較大差異，在某些情況下甚至?xí)霈F(xiàn)相反的現(xiàn)象，過去文獻(xiàn)中鮮有類似報(bào)道.為了證實(shí)這一現(xiàn)象，以SV橫波入射到裸眼井中為例，進(jìn)一步通過三維有限差分法模擬不同頻率下井中接收波場的方位響應(yīng)特征.差分模擬時(shí)聲源上方共有32個(gè)接收器等間隔環(huán)向布置，其他參數(shù)與圖3中相同.考慮聲源頻率為1500 Hz和4000 Hz兩種情況，三維有限差分的計(jì)算結(jié)果分別如圖12中圓圈和三角形標(biāo)記所示，其中圖12a、b分別為流體徑向位移和聲壓的方位響應(yīng).為了對(duì)比，圖中同時(shí)給出了利用式(12)得到的解析結(jié)果(實(shí)線和虛線).從圖中可以看出，無論頻率為1500 Hz還是4000 Hz時(shí)，兩種方法計(jì)算得到的流體徑向位移和聲壓的方位響應(yīng)結(jié)果均吻合較好.此外，在不同頻率下，兩種方法均顯示出了相反的結(jié)果：當(dāng)頻率為1500 Hz時(shí)，靠近(遠(yuǎn)離)反射體一側(cè)的接收器接收到的徑向位移(聲壓)幅值最大；而當(dāng)頻率為4000 Hz時(shí)，則是遠(yuǎn)離(靠近)反射體一側(cè)的接收器將接收到最大幅值的徑向位移(聲壓).三維有限差分計(jì)算結(jié)果證實(shí)了方位幅值響應(yīng)存在180°反轉(zhuǎn)現(xiàn)象，指示利用井中接收波場幅值特征進(jìn)行井外反射體方位識(shí)別時(shí)需要考慮接收器偏心距、聲源頻率和地層性質(zhì)等因素的影響，這對(duì)于儀器研制和現(xiàn)場解釋具有重要的指導(dǎo)意義.

圖12 快速地層條件下，聲源頻率為1500 Hz和4000 Hz時(shí)，SV橫波入射引起的井中接收波場的方位響應(yīng)(a) 徑向位移； (b) 聲壓.為了對(duì)比，同時(shí)給出了漸近解(曲線)和三維有限差分(標(biāo)記)的計(jì)算結(jié)果.Fig.12 The azimuthal response of waveforms inside borehole at frequencies 1500 Hz and 4000 Hz for SV-wave incidence in the fast formation(a) Radial displacement; (b) Pressure. The curves and circles are calculated using analytical (curves) and 3D finite-difference (markers) modeling, respectively.

5 結(jié)論

本文采用解析法和三維有限差分法對(duì)偶極橫波遠(yuǎn)探測在裸眼和套管井中接收波場的方位響應(yīng)特征進(jìn)行了研究，得到以下認(rèn)識(shí)和結(jié)論：

(1)結(jié)合虛源和球面波的柱面波展開，采用最速下降法得到了充液裸眼和套管井中偶極橫波方位遠(yuǎn)探測接收波場的漸近解析解.通過解析法和三維有限差分法對(duì)套管井中方位接收波場的模擬對(duì)比，驗(yàn)證了漸近解的正確性和普適性.

(2)SH橫波入射引起的井內(nèi)流體的徑向位移和聲壓的方位響應(yīng)特征具有較好的一致性，而SV橫波入射時(shí)二者特征卻存在較大差異.值得注意的是，SV橫波入射時(shí)，井內(nèi)波場的方位幅值響應(yīng)對(duì)井外反射體位置具有唯一性，即在常用偶極測井頻率下，靠近或遠(yuǎn)離反射體的接收器接收到的波形幅值最大.這表明在進(jìn)行聲波遠(yuǎn)探測成像解釋時(shí)，考慮井內(nèi)接收信號(hào)中的SV分量有助于識(shí)別井外反射體的方位，但需要注意接收換能器記錄的信號(hào)類型(徑向位移/聲壓)的不同響應(yīng).

(3)偶極橫波遠(yuǎn)探測井內(nèi)接收波場的方位響應(yīng)特征受井眼環(huán)境、接收器偏心距、聲源頻率和地層性質(zhì)等多個(gè)因素的影響，在不同參數(shù)下方位幅值響應(yīng)甚至?xí)霈F(xiàn)反轉(zhuǎn)現(xiàn)象，漸近解析解和三維有限差分?jǐn)?shù)值解的結(jié)果對(duì)比證實(shí)了這一現(xiàn)象，實(shí)際數(shù)據(jù)解釋時(shí)需要考慮這些因素的影響.此外，在影響井內(nèi)聲場方位響應(yīng)的眾多因素當(dāng)中，最重要的主控因素是接收器偏心距.正是由于偏心接收才有文中不同于偶極180°對(duì)稱性的變化現(xiàn)象(即圖4—圖12所示結(jié)果)，其物理實(shí)質(zhì)是聲場在偏離井心位置上的復(fù)雜三維變化；而居中時(shí)三維變化就降到了二維，復(fù)雜性減少但(180°)對(duì)稱現(xiàn)象出現(xiàn).

(4)相對(duì)于裸眼井，套管井中接收波場的方位響應(yīng)對(duì)井外反射體位置的靈敏度更高，受接收器偏心距、聲源頻率和地層性質(zhì)的影響也更大，說明在套管井中進(jìn)行偶極橫波方位遠(yuǎn)探測更有利于井外反射體的方位識(shí)別.

(5)現(xiàn)有的偶極橫波測井儀將儀器兩側(cè)的接收換能器記錄的信號(hào)之差作為測量信號(hào)，此時(shí)無論換能器記錄的是流體徑向位移還是聲壓，相減后的測量信號(hào)中都主要包含偶極成分，可以得到反射體的走向信息；接收器的偏心效應(yīng)包含了反射體的傾向信息，有助于解決偶極橫波遠(yuǎn)探測方位的180°不確定性問題.

附錄A 裸眼井情況下矩陣方程中的元素表達(dá)式

矩陣M中各元素的具體表達(dá)式為：

M11=-nIn(fa)/a-fIn+1(fa),

M12=nKn(pa)/a-pKn+1(pa),

M13=nKn(sa)/a,

M14=ik[nKn(sa)/a-sKn+1(sa)],

M21=ρfω2In(fa),

M22=2ρβ2[(n2-n)Kn(pa)/a2+pKn+1(pa)/a]

+ρ(2k2β2-ω2)Kn(pa),

M23=2nρβ2[(n-1)Kn(sa)/a2-sKn+1(sa)/a],

M24=2ikρβ2{[(n2-n)/a2+s2]Kn(sa)

+sKn+1(sa)/a},

M31=0,

M32=2nρβ2[(1-n)Kn(pa)/a2+pKn+1(pa)/a],

M33=-ρβ2{[2(n2-n)/a2+s2]Kn(sa)

+2sKn+1(sa)/a},

M34=2iknρβ2[(1-n)Kn(sa)/a2+sKn+1(sa)/a],

M41=0,

M42=2ikρβ2[nKn(pa)/a-pKn+1(pa)],

M43=iknρβ2Kn(sa)/a.

M44=-(s2+k2)ρβ2[nKn(sa)/a-sKn+1(sa)].

(A1)

向量b和b′中的元素由于入射波的類型不同而有所區(qū)別.當(dāng)SH橫波入射時(shí)，向量b和b′中各元素的具體表達(dá)式為(省略公共因子S(ω)·RDSH(ω)·RFSH(ω))：

b1=-nεnIn(sa)sin(nφ0)/a,

b2=-2nεnρβ2[(n-1)In(sa)/a2

+sIn+1(sa)/a]sin(nφ0),

b3=εnρβ2{[2(n2-n)/a2+s2]In(sa)

-2sIn+1(sa)/a}sin(nφ0),

b4=-iknεnρβ2In(sa)sin(nφ0)/a,

b′1=nεnIn(sa)cos(nφ0)/a,

b′2=2nεnρβ2[(n-1)In(sa)/a2+sIn+1(sa)/a]

×cos(nφ0),

b′3=-εnρβ2{[2(n2-n)/a2+s2]In(sa)

-2sIn+1(sa)/a}cos(nφ0),

b′4=iknεnρβ2In(sa)cos(nφ0)/a.

(A2)

當(dāng)SV橫波入射時(shí)，向量b和b′中各元素的具體表達(dá)式為(省略公共因子S(ω)·RDSV(ω)·RFSV(ω))：

b1=-ikεn[nIn(sa)/a+sIn+1(sa)]cos(nφ0),

b2=-2ikεnρβ2{[(n2-n)/a2+s2]In(sa)

-sIn+1(sa)/a}cos(nφ0),

b3=-2iknεnρβ2[(1-n)In(sa)/a2-sIn+1(sa)/a]

×cos(nφ0),

b4=εn(s2+k2)ρβ2[nIn(sa)/a+sIn+1(sa)]cos(nφ0),

b′1=-ikεn[nIn(sa)/a+sIn+1(sa)]sin(nφ0),

b′2=-2ikεnρβ2{[(n2-n)/a2+s2]In(sa)

-sIn+1(sa)/a}sin(nφ0),

b′3=-2iknεnρβ2[(1-n)In(sa)/a2

-sIn+1(sa)/a]sin(nφ0),

b′4=εn(s2+k2)ρβ2[nIn(sa)/a+sIn+1(sa)]sin(nφ0)，

(A3)

附錄B 套管井(膠結(jié)良好)情況下T矩陣的元素表達(dá)式

T矩陣的元素表達(dá)式為：

T11=nIn(pR)/R+pIn+1(pR),

T12=nKn(pR)/R-pKn+1(pR),

T13=nIn(sR)/R,

T14=nKn(sR)/R,

T15=ik[nIn(sR)/R+sIn+1(sR)],

T16=ik[nKn(sR)/R-sKn+1(sR)],

T21=-nIn(pR)/R,

T22=-nKn(pR)/R,

T23=-nIn(sR)/R-sIn+1(sR),

T24=-nKn(sR)/R+sKn+1(sR),

T25=-iknIn(sR)/R,

T26=-iknKn(sR)/R,

T31=ikIn(pR),

T32=ikKn(pR),

T33=0,

T34=0,

T35=-s2In(sR),

T36=-s2Kn(sR),

T41=2ρβ2[(n2-n)In(pR)/R2-pIn+1(pR)/R]

+ρ(2k2β2-ω2)In(pR),

T42=2ρβ2[(n2-n)Kn(pR)/R2+pKn+1(pR)/R]

+ρ(2k2β2-ω2)Kn(pR),

T43=2nρβ2[(n-1)In(sR)/R2+sIn+1(sR)/R],

T44=2nρβ2[(n-1)Kn(sR)/R2-sKn+1(sR)/R],

T45=2ikρβ2{[(n2-n)/R2+s2]In(sR)

-sIn+1(sR)/R},

T46=2ikρβ2{[(n2-n)/R2+s2]Kn(sR)

+sKn+1(sR)/R},

T51=2nρβ2[(1-n)In(pR)/R2-pIn+1(pR)/R],

T52=2nρβ2[(1-n)Kn(pR)/R2+pKn+1(pR)/R],

T53=-ρβ2{[2(n2-n)/R2+s2]In(sR)

-2sIn+1(sR)/R},

T54=-ρβ2{[2(n2-n)/R2+s2]Kn(sR)

+2sKn+1(sR)/R},

T55=2iknρβ2[(1-n)In(sR)/R2-sIn+1(sR)/R],

T56=2iknρβ2[(1-n)Kn(sR)/R2+sKn+1(sR)/R],

T61=2ikρβ2[nIn(pR)/R+pIn+1(pR)],

T62=2ikρβ2[nKn(pR)/R-pKn+1(pR)],

T63=iknρβ2In(sR)/R,

T64=iknρβ2Kn(sR)/R,

T65=-(s2+k2)ρβ2[nIn(sR)/R+sIn+1(sR)],

T66=-(s2+k2)ρβ2[nKn(sR)/R-sKn+1(sR)].

(B1)

附錄C 套管井(自由套管)情況下矩陣方程中的元素表達(dá)式

矩陣H中各元素的具體表達(dá)式為(未列出的元素取值為0)：

H11=-nIn(fr1)/r1-fIn+1(fr1),

H58=-nIn(fr2)/r2-fIn+1(fr2),

H59=-nKn(fr2)/r2+fKn+1(fr2),

H69=ρfω2Kn(fr2),

H98=-nIn(fa)/a-fIn+1(fa),

H99=-nKn(fa)/a+fKn+1(fa),

H10,8=ρfω2In(fa),H10，9=ρfω2Kn(fa),

(C1)

其中，r1和r2分別為套管的內(nèi)半徑和外半徑；T的右上標(biāo)表示計(jì)算時(shí)半徑R的取值.

根據(jù)入射波類型的不同(即SH或SV橫波入射)，向量d和d′中的元素也有所區(qū)別，其各元素的具體表達(dá)式為(未列出的元素取值為0)：

d9=b1,d10=b2,d11=b3,d12=b4,

d′9=b′1,d′10=b′2,d′11=b′3,d′12=b′4.

(C2)

附錄D 井軸處的聲壓和徑向位移

首先求解井軸處的聲壓，第一類變型貝塞爾函數(shù)具有如下性質(zhì)：

(D1)

其中，n為階數(shù).接收器位于井軸時(shí)，式(12)中的r=0，由式(D1)可知，此時(shí)聲壓的多極響應(yīng)中只有n=0這一項(xiàng)有值，其余全為0，表達(dá)式為：

(D2)

求解井軸處的徑向位移時(shí)，需要用到第一類變型貝塞爾函數(shù)的遞推公式：

(D3)

其中，x為變量.將式(D3)代入式(12)的徑向位移表達(dá)式中可得：

×[In-1(f0r)+In+1(f0r)].

(D4)

n=0時(shí)，式(D4)中將出現(xiàn)負(fù)數(shù)階貝塞爾函數(shù)，可根據(jù)整階貝塞爾函數(shù)的奇偶性進(jìn)行計(jì)算，公式為：

I-n(x)=In(x).

(D5)

令式(D4)中的r=0，將式(D1)和(D5)代入式(D4)中，得到井軸處的徑向位移為：

(D6)

式(D6)表明，井軸處的徑向位移只有偶極項(xiàng)(n=1)貢獻(xiàn).