李 超 張 波 許孝凱 張晉言 陳 浩，4 何 曉 王秀明，4

(1 中國(guó)科學(xué)院聲學(xué)研究所 聲場(chǎng)聲信息國(guó)家重點(diǎn)實(shí)驗(yàn)室 北京 100190)

(2 北京市海洋深部鉆探測(cè)量工程技術(shù)研究中心 北京 100190)

(3 中石化勝利石油工程有限公司測(cè)井公司 東營(yíng) 257096)

(4 中國(guó)科學(xué)院大學(xué) 北京 100049)

0 引言

利用陣列聲波測(cè)井中記錄到的反射波來(lái)進(jìn)行井旁構(gòu)造成像的遠(yuǎn)探測(cè)聲波測(cè)井技術(shù)是近年來(lái)聲波測(cè)井領(lǐng)域的重要進(jìn)展之一，它可以探測(cè)裂縫、溶洞等儲(chǔ)層在井外的擴(kuò)展情況，以及水平井的儲(chǔ)層邊界，還可以識(shí)別地震勘探無(wú)法探測(cè)到的小構(gòu)造等[1-2]，目前逐漸成為復(fù)雜油氣藏勘探的有力手段。隨著測(cè)井儀器測(cè)量方式和數(shù)字信號(hào)處理技術(shù)的發(fā)展，出現(xiàn)了利用具有方位接收器的陣列聲波測(cè)井儀器探測(cè)井外構(gòu)造方位信息的技術(shù)[3-5]，以及利用正交偶極聲波測(cè)井儀測(cè)量的反射橫波來(lái)進(jìn)行反射體方位識(shí)別的技術(shù)[2，6-7]，相對(duì)于傳統(tǒng)的單極子探測(cè)技術(shù)，這些新技術(shù)有更好的探測(cè)效果以及方位識(shí)別能力[8-9]，而反射體方位信息對(duì)于儲(chǔ)層評(píng)價(jià)以及工程措施如壓裂改造等至關(guān)重要。

目前，已有的單偶極遠(yuǎn)探測(cè)方位識(shí)別方法均各有優(yōu)勢(shì)。單極子聲源頻率高，不同方位的接收器對(duì)于反射方位敏感，可以根據(jù)波形幅度信息快速判斷方位，但是其探測(cè)深度較淺，方位分辨率受到方位接收器個(gè)數(shù)的影響；而偶極遠(yuǎn)探測(cè)方法探測(cè)深度較大，SH 反射橫波幅度更大，更容易探測(cè)井外深部構(gòu)造，且對(duì)于構(gòu)造走向較為敏感，但其存在180°的方位不確定性[2]。本文通過將單偶極遠(yuǎn)探測(cè)識(shí)別方位的方法相結(jié)合，利用多分量遠(yuǎn)探測(cè)儀器，集合兩類方法的優(yōu)勢(shì)，研究解決反射方位識(shí)別問題的方法。

1 基本原理

1.1 多分量遠(yuǎn)探測(cè)儀器探測(cè)原理

偶極橫波多分量遠(yuǎn)探測(cè)儀器(Victory multicomponent reflected shear wave imager，VMSI)[9]同時(shí)包含了單極子聲源，儀器工作時(shí)單偶極聲源依次激發(fā)聲波，四個(gè)周向上相隔90°的接收器記錄不同方位的四個(gè)分量波形，如圖1所示。單極子記錄MX1、MX2、MY1、MY2 四個(gè)分量，偶極子記錄XX1、XX2、XY1、XY2 以及Y X1、Y X2、Y Y1、Y Y2 共8 個(gè)分量，其中第一個(gè)字母代表聲源，第二個(gè)字母及數(shù)字代表接收器方位。則正交偶極的四分量波形可由式(1)獲得：

圖1 多分量遠(yuǎn)探測(cè)儀器原理圖Fig.1 Schematic diagram of multicomponent remote exploration acoustic logging instrument

1.2 方位識(shí)別方法

當(dāng)反射體到井眼的距離遠(yuǎn)大于井徑時(shí)，反射波可以看作是平面波入射到井孔。以井外不同平面波入射到井內(nèi)被不同方位接收器所接收到的波形特征為基礎(chǔ)，來(lái)說明反射波方位的識(shí)別方法，如圖2所示。圖2(a)顯示了平面波入射井孔以及井內(nèi)四個(gè)不同方位接收器分布的示意圖，采用二維有限差分進(jìn)行模擬，平面波聲源為雷克子波，模擬縱波(單極子)入射時(shí)聲源主頻為8 kHz，橫波(偶極子)時(shí)聲源主頻為4 kHz，模型參數(shù)如表1所示。井孔內(nèi)有一居中的儀器，井孔直徑為0.22 m，儀器直徑為0.1 m，井外地層為硬地層(地層1)。

圖2(b)顯示了平面縱波入射時(shí)四個(gè)接收器接收到的波形，由于頻率較高，且質(zhì)點(diǎn)振動(dòng)方向與波傳播方向相同，儀器的存在使得不同接收器波形幅度不同，靠近入射方向的一側(cè)X1 接收器幅度最強(qiáng)，Y1 和Y2 接收器首波幅度相等，X2 首波最弱，因而可以根據(jù)首波幅度大小快速判斷來(lái)波方位，即反射界面(體)方位，用公式可表示如下：

其中，A表示振幅，AZM表示最大振幅接收器所在方位，若有兩個(gè)接收器振幅都大于其他兩個(gè)且相等，則取兩個(gè)接收器方位的平均。由于方位接收器個(gè)數(shù)的限制，方位分辨率只有45°(360°/2N，N為方位接收器個(gè)數(shù)，本文N= 4)。這里方位分辨率指的是可以分辨的兩個(gè)方位間的最小角度，即兩個(gè)目標(biāo)在這個(gè)角度內(nèi)將不能區(qū)分。

圖2(c)顯示了平面SH 橫波入射時(shí)四個(gè)接收器接收到的波形，可以類比偶極聲源發(fā)射得到的接收波形，由于頻率較低，波形受儀器影響較小，而且SH波入射后與縱波的特征完全不同。質(zhì)點(diǎn)振動(dòng)方向與波傳播方向垂直，當(dāng)入射進(jìn)圓形井孔時(shí)(這里波長(zhǎng)大于井孔尺寸)，橫波轉(zhuǎn)換成井孔內(nèi)縱波(接收波形為井孔內(nèi)聲壓)，Y1和Y2接收器的位置與其傳播方向垂直，且該位置波傳播方向與井孔邊界有夾角，SH 波可以轉(zhuǎn)換成井內(nèi)縱波，又在橫波振動(dòng)方向，所以振幅最大，相位相反。而X1 和X2 接收器平行于波傳播方向，橫波垂直入射井孔，所以振幅為0。這些是由SH 波的偏振特征所致，因而也可以反映方位變化，即反射界面的走向[2]。將式(1)得到的四分量數(shù)據(jù)進(jìn)行旋轉(zhuǎn)，可以獲得SH波和SV波不同方位的波場(chǎng)切片：

其中，角度θ為走向角，對(duì)于遠(yuǎn)場(chǎng)而言，一般SH反射波的幅度大于SV 波[2]，SH 反射能量最大的角度即為反射界面走向角度，即偶極源極化方向和反射界面平行時(shí)，產(chǎn)生的SH 波反射幅度最大。然而由于偶極源的輻射特性，以上方法求得的走向存在180°的不確定性，即無(wú)法區(qū)分在井外的哪一側(cè)。

表1 模型參數(shù)Table1 Model parameters

圖2 不同平面波入射到井孔波形特征Fig.2 Wavefield characteristics in the borehole for different incident plane waves

本文提出應(yīng)用于上述多分量遠(yuǎn)探測(cè)儀器的單極和偶極遠(yuǎn)探測(cè)相結(jié)合的方法來(lái)識(shí)別反射方位，一方面單極反射幅度能較好地區(qū)分方位，但是其方位分辨率低；另一方面，偶極SH 反射能較精確地反演界面走向，但是無(wú)法區(qū)分界面在井哪一側(cè)，將兩者相結(jié)合，可以優(yōu)勢(shì)互補(bǔ)，快速判斷界面方位。本文將以識(shí)別臨井方位為例，通過數(shù)值模擬及實(shí)際資料處理來(lái)說明該方法的效果。

2 數(shù)值模擬實(shí)例

采用三維有限差分方法進(jìn)行數(shù)值模擬來(lái)獲取單偶極遠(yuǎn)探測(cè)波形。模型的xz和xy剖面如圖3所示，兩口相距5 m的井，一口為測(cè)量井，井中有儀器，地層為軟地層(見表1地層2)，另一口井在其15°方位(與x軸夾角)處，儀器結(jié)構(gòu)如圖1所示，單極子聲源主頻為10 kHz，偶極子聲源主頻為3 kHz。兩口井直徑均為0.22 m，儀器直徑為0.1 m。

X方向的接收器在x軸放置，并且X1 偏向有井的一側(cè)，參考方位為X1 接收器方位(0°)。數(shù)值模擬得到的源距相同的不同方位單極子波形見圖4，其中圖4(a)表明單極子反射波有多種模式，在時(shí)間上從前往后依次是P-P 反射縱波、P-S 和S-P 轉(zhuǎn)換波以及S-S 反射橫波，并且不同方位波形幅度有差異，不同方位最大反射波振幅歸一化后的指向性圖見圖4(b)，可以明顯看到方位在X1 接收器一側(cè)，但由于分辨率的限制，無(wú)法獲得準(zhǔn)確的方位角。

模擬得到的偶極子四分量波形見圖5。由于X源方向與臨井的夾角較小，反射縱波、轉(zhuǎn)換波和反射橫波均存在(圖5(a))，這一點(diǎn)和單極子源反射波場(chǎng)類似，而交叉分量反射波幅度很小，Y Y分量主要成分為SH反射橫波，其幅度最大(圖5(d))。利用公式(3)求得不同走向角度的SH 反射波并求出其最大幅度，結(jié)果見圖6。可知最大幅度對(duì)應(yīng)走向角度為79°，這里的走向角度是與x軸方向(X1 接收器方向)的順時(shí)針夾角，轉(zhuǎn)換成模型中的逆時(shí)針方位即為11°與191°，即偶極SH反射波所確定的方位存在180°的不確定性，根據(jù)單極子的結(jié)果(圖4(b))，可以判斷11°為臨井的方位，這與其真實(shí)的方位角15°誤差很小，因而將單偶極子的反射波方位識(shí)別方法聯(lián)合起來(lái)可以消除不確定性，較準(zhǔn)確地識(shí)別出反射方位。

圖3 井旁15°井速度模型Fig.3 Velocity model of a 15°well outside the borehole

圖4 單極子反射及最大幅度指向性Fig.4 Monopole reflection and its maximal amplitude directivity

圖5 偶極子四分量反射波Fig.5 Dipole reflection of four components

圖6 SH 波不同走向角度最大振幅Fig.6 The maximum amplitude of SH reflection for different strike angles

3 實(shí)際資料處理

與數(shù)值模擬模型類似，多分量遠(yuǎn)探測(cè)儀器在勝利油田某一教學(xué)井中測(cè)量，來(lái)探測(cè)井外的另一口井，兩井相距約為5 m，目標(biāo)井方位約在測(cè)量井的西側(cè)，為一套管井，直徑約0.14 m。測(cè)量井段為軟地層，地層參數(shù)與表1中地層2 相近。測(cè)量得到的單極子四分量波形反射波成像結(jié)果見圖7(四個(gè)圖形色標(biāo)值范圍一致)，四個(gè)方位接收器都能探測(cè)到臨井的存在，140～160 m 井段可能由于地層衰減較大沒有測(cè)到反射信號(hào)。由于儀器上提時(shí)旋轉(zhuǎn)，不同深度不同接收器的反射波幅度大小存在不一致，根據(jù)公式(2)，將不同深度的反射波最大幅度以及方位信息進(jìn)行綜合分析，可得圖7(e)的結(jié)果，參考方位為Y1 接收器相對(duì)北極的方位。對(duì)于有較明顯反射的115～125 m 深度段和170～173 m 深度段，可以看出最大振幅方位角在250°～326°這一區(qū)間，位于西側(cè)，但是由于分辨率的限制，無(wú)法確定準(zhǔn)確方位，此外，由于噪聲以及儀器本身旋轉(zhuǎn)的影響，部分深度段方位角出現(xiàn)較大誤差，因此需要將多個(gè)深度點(diǎn)信息結(jié)合起來(lái)進(jìn)行綜合分析。

偶極子四分量反射波成像結(jié)果及走向角分析如圖8所示?？梢姡瓷渲饕性赬X分量上，說明走向在南北向(方位校正后XX指向南北向)，根據(jù)SH 波最大幅度來(lái)進(jìn)行進(jìn)一步的走向角分析(圖8(e))，由圖可知走向角集中在156°附近(方位在66°和246°兩個(gè)方向)，尤其是深度170 m 附近的走向角度比較穩(wěn)定，結(jié)合反射波成像結(jié)果，此處的角度應(yīng)比較可靠。再根據(jù)圖7單極子的結(jié)果，說明246°的方位較為正確，即另一口井在正西稍偏南方向，這也與實(shí)際相符。

圖7 單極子不同方位接收器反射波成像Fig.7 Monopole reflection imaging results for receivers in different directions

圖8 偶極子不同分量反射波成像Fig.8 Dipole reflection imaging results for different components

為了進(jìn)一步說明，圖9和圖10給出了SH 反射波在不同走向角度的剖面，可以看出，反射在150°和165°時(shí)幅度較強(qiáng)，這也與圖8的分析一致。

圖9 偶極子不同角度SH 反射波成像(0°～75°)Fig.9 Dipole imaging results for SH reflection with different strike (0°～75°)

圖10 偶極子不同角度SH 反射波成像(90°～165°)Fig.10 Dipole imaging results for SH reflection with different strike (90°～165°)

4 結(jié)論

多分量遠(yuǎn)探測(cè)聲波測(cè)井方法集成了測(cè)量單偶極不同方位多分量波場(chǎng)的功能。利用單極子不同方位反射波的振幅差異，可以快速判斷反射界面大致方位，但是缺乏準(zhǔn)確性；而偶極SH 反射可以較準(zhǔn)確地識(shí)別走向(和方位垂直)，但是存在180°的不確定性，將兩者結(jié)合，可以快速準(zhǔn)確地識(shí)別出反射界面方位，并且兩種方法的聯(lián)合識(shí)別，提高了結(jié)果的可靠性。本文的研究為一定徑向深度范圍內(nèi)的反射界面方位識(shí)別問題提供了思路，即能識(shí)別遠(yuǎn)探測(cè)方位的距離只受限于單極源的探測(cè)深度，對(duì)于較近反射以及離井眼遠(yuǎn)近都有的反射(與井眼相交傾斜界面)可以識(shí)別，但難以識(shí)別較遠(yuǎn)且與井眼不相交的反射面，同時(shí)也為傳統(tǒng)多極子陣列聲波測(cè)井儀用于遠(yuǎn)探測(cè)時(shí)的信號(hào)記錄方式提供了指導(dǎo)。

致謝感謝中國(guó)科學(xué)院聲學(xué)研究所的賀洪斌高工以及中國(guó)石化勝利石油工程公司測(cè)井公司晁永勝高工為本項(xiàng)研究所提供的優(yōu)質(zhì)數(shù)據(jù)。