徐基祥，McLean B F，宋雪娟

（1.中國石油勘探開發(fā)研究院；2.北京勘源科技有限公司）

近地表地震散射波分離方法

徐基祥1，McLean B F2，宋雪娟1

（1.中國石油勘探開發(fā)研究院；2.北京勘源科技有限公司）

針對常規(guī)去噪處理無法壓制近地表散射波的問題，通過分析近地表地震散射波特征，利用地震干涉測量理論，實(shí)現(xiàn)近地表地震散射波分離。山地地震資料信噪比普遍較低，有效波能量弱，近地表散射波能量強(qiáng)，地震數(shù)據(jù)處理中必須分離近地表散射波，提高地震資料信噪比。相鄰炮集上的近地表地震散射波同相軸相互平行，根據(jù)近地表散射波運(yùn)動學(xué)規(guī)律，通過地震干涉測量，消除地震波傳播過程中相同路徑的影響，近地表散射波得到相干加強(qiáng)，進(jìn)而從原始地震記錄中將其分離。以中國西部山地地震資料為例，通過直達(dá)波場與全波場地震干涉測量，加強(qiáng)散射波能量，進(jìn)行散射波分離，處理后的地震資料基本消除了近地表散射波，驗(yàn)證了近地表地震散射波分離方法的有效性和可行性。圖16參17

山地地震勘探；近地表；地震散射波；散射波分離；地震干涉測量

0 引言

地震勘探中，地震波從激發(fā)到接收一般要經(jīng)過近地表兩次。近地表為不均勻介質(zhì)時(shí)，地震波要經(jīng)歷兩次散射，這種散射為近地表散射。更嚴(yán)重的近地表散射是不均勻近地表對沿地表傳播的地震波（如直達(dá)波、折射波和面波）的散射作用，由于這類波能量強(qiáng)，會使近地表散射影響更加明顯。近地表散射波由近地表不均勻介質(zhì)造成，與勘探目的層無關(guān)（或關(guān)系不大）。山地地震勘探中，不均勻近地表嚴(yán)重影響了地震波傳播和散射，降低了地震資料品質(zhì)，甚至造成復(fù)雜構(gòu)造成像假象。

提高地震勘探成功率的首要任務(wù)是提高地震資料的信噪比。目前地震成像技術(shù)僅將一次反射波作為有效信號，將其他信號都當(dāng)作噪聲，因而提高地震資料信噪比的傳統(tǒng)方法主要是根據(jù)一次反射波與噪聲的差異進(jìn)行噪聲壓制。進(jìn)行噪聲壓制時(shí)，通常把近地表散射波當(dāng)作隨機(jī)噪聲處理，然而近地表散射波有其自身規(guī)律，并非隨機(jī)噪聲。壓制近地表散射噪聲可用逆散射方法來衰減實(shí)際資料中的瑞雷波散射，但該方法反演計(jì)算量大，且只能衰減瑞雷波散射噪聲[1-5]。

吳希光等[6-8]根據(jù)大量地震實(shí)例，總結(jié)了表層散射波的一些規(guī)律，認(rèn)為近地表散射波是造成地震資料信噪比低的根本原因，并提出野外采集中“避強(qiáng)散射波，就弱散射波”的對策，指出增加檢波器垂向組合是提高復(fù)雜地區(qū)地震資料信噪比的基本方法。本文從地震資料室內(nèi)處理角度分析近地表地震散射波分離方法，針對近地表地震散射波同相軸特征，利用地震干涉測量原理，實(shí)現(xiàn)近地表地震散射波分離。

1 近地表地震散射波分離重要性

山地地震勘探中，粗糙起伏地形和不均勻近地表都會產(chǎn)生散射波。散射波嚴(yán)重影響地震數(shù)據(jù)成像質(zhì)量：①降低地震資料信噪比；②近地表散射波在時(shí)間、空間和頻率上都與體波重疊，若其參與常規(guī)偏移成像處理，將會產(chǎn)生假構(gòu)造。研究近地表地震散射波分離方法具有理論和現(xiàn)實(shí)意義，不僅可提高地震資料信噪比，而且能避免復(fù)雜構(gòu)造成像假象。

與常規(guī)地震資料相比，山地地震資料信噪比明顯降低，同相軸連續(xù)性變差，體波信息模糊，散射波更加明顯。圖1所示為實(shí)際山前地震記錄：炮點(diǎn)位于平坦區(qū)，240～410道為平坦區(qū)接收，411～440道為山坡區(qū)接收。地震記錄上直觀可見第410道為分界點(diǎn)：左側(cè)初至波波形穩(wěn)定、線性同相軸光滑連續(xù)；右側(cè)初至波波形較穩(wěn)定但同相軸隨地形變化而跳躍。炮集中初至波信息反映地形起伏或近地表巖性變化情況，光滑地表對應(yīng)高信噪比資料，粗糙地表會產(chǎn)生明顯散射波，因此對應(yīng)于低信噪比資料，必須分離近地表散射波與有效波。

圖1 山前單炮地震記錄

圖2所示為山地近地表速度分布及其單炮地震記錄：圖2a為近地表速度分布圖，表明該區(qū)地形有一定起伏；產(chǎn)生近地表散射波主要原因在于近地表速度突變，如圖2b中Sp和Ss-p，Sp是向左傳播的直達(dá)波遇到速度突變體產(chǎn)生的反向散射波，同相軸呈線性分布，Ss-p是Sp在向右傳播過程中遇到速度突變體再次產(chǎn)生的散射波，同相軸呈雙曲線分布，雙曲線頂點(diǎn)對應(yīng)近地表速度突變點(diǎn)，其旅行時(shí)為Sp到達(dá)該點(diǎn)的時(shí)間。從圖2b中可清楚看出，近地表散射波疊置于地震體波中，且分布范圍廣，頻帶與體波頻帶相當(dāng)。

圖2 近地表速度分布及其單炮地震記錄

近地表線性散射波比較容易識別并壓制，而對于近地表雙曲線散射波，由于其與復(fù)雜構(gòu)造產(chǎn)生的一次反射波非常相似，往往誤將其當(dāng)成有效波，因此難以識別。若近地表散射波與有效波共同成像，必然產(chǎn)生構(gòu)造假象，為此必須在成像處理前，分離近地表散射波。

2 近地表地震散射波運(yùn)動學(xué)規(guī)律

高陡山體、崎嶇地形、戈壁礫石、溶縫溶洞、巖性突變體等不均勻體都是強(qiáng)散射源，其產(chǎn)生的次生地震散射波向四周傳播。以一個簡單模型為例說明近地表散射波的運(yùn)動學(xué)規(guī)律：模型為均勻介質(zhì)，在近地表設(shè)置3個散射源。圖3為近地表散射波傳播路徑示意圖：S為震源，r1—r4為接收點(diǎn)，s1—s3為散射源，s1在測線上，s2在測線外地表上，s3在測線外地下，s2′是s2在測線上的垂直投影，s3′是s3在測線上的垂直投影。紅色箭頭代表沿地表或近地表傳播的直達(dá)波、折射波或面波，藍(lán)色箭頭代表這些波遇到散射源后產(chǎn)生的散射波。

圖3 近地表散射波射線路徑示意圖

對應(yīng)于圖3的地震波時(shí)距曲線如圖4所示，紅線代表沿地表或近地表傳播的直達(dá)波、折射波和面波的時(shí)距曲線，藍(lán)線代表散射波時(shí)距曲線，包括散射直達(dá)波、散射折射波和散射面波時(shí)距曲線。圖4僅顯示一次反向散射波和側(cè)面散射波，沒有考慮多次散射波和轉(zhuǎn)換波，若考慮轉(zhuǎn)換波，就能用不同斜率清晰表示正向散射波。圖4中側(cè)面散射波包括側(cè)面散射折射波和側(cè)面散射直達(dá)波。

圖4 直達(dá)波、折射波和面波散射波同相軸示意圖

以散射源s3為例，分析近地表散射波運(yùn)動學(xué)規(guī)律。震源和接收點(diǎn)在同一水平面上，震源S的坐標(biāo)為（0，0，0），散射源s3的坐標(biāo)為（x0，y0，z0），接收點(diǎn)r3的坐標(biāo)為（x，0，0），背景介質(zhì)速度為v，以散射直達(dá)波為例，由空間幾何關(guān)系可以得到散射波旅行時(shí)t為：

其中t1是從震源S到散射源s3的旅行時(shí)，t2是從散射源s3到接收點(diǎn)r3的旅行時(shí)。若震源和接收點(diǎn)不在同一水平面，接收點(diǎn)r3的坐標(biāo)為（x，0，z），則散射波旅行時(shí)t為：

由此可知近地表地震散射波運(yùn)動學(xué)規(guī)律（以散射直達(dá)波為例）：①地震波從震源到散射源的傳播時(shí)間取決于震源和散射源位置，而從散射源到接收點(diǎn)的傳播時(shí)間取決于散射源和接收點(diǎn)位置，與震源位置無關(guān)，因此改變震源位置，從散射源到接收點(diǎn)的傳播時(shí)間保持不變；②當(dāng)散射源位于測線上，此時(shí)側(cè)面散射波的時(shí)距曲線從雙曲線變成直線，包括正向散射波和反向散射波；③地震波從震源到散射源的傳播時(shí)間決定散射波雙曲線頂點(diǎn)對應(yīng)的時(shí)間，也決定線性正向散射波和反向散射波的起始時(shí)間。

散射直達(dá)波運(yùn)動學(xué)規(guī)律同樣適用于散射折射波和散射面波，因此得到結(jié)論：相鄰炮集上的近地表散射波同相軸平行。該結(jié)論具有普遍性，在復(fù)雜地表區(qū)同樣成立，其為分離近地表地震散射波提供了物理基礎(chǔ)。

3 近地表地震散射波分離理論分析

地震干涉測量法是近地表散射波分離的一種技術(shù)手段。2001年Schuster首次提出“地震干涉測量”概念，隨后系統(tǒng)介紹了地震干涉測量的定義、原理和應(yīng)用[9]，2002年Wapenaar等[10]給出了地震干涉測量格林函數(shù)表達(dá)式，2006年Dong等[11]用地震干涉測量法預(yù)測面波。2008年Halliday和Curtis[12]，2010年Forghani和Snieder[13]闡述了地震干涉測量預(yù)測面波的可行性，同時(shí)Fleury等[14]結(jié)合地震干涉測量和擾動理論，介紹了散射波重建，Snieder等[15-16]和Halliday等[17]用廣義光學(xué)定理說明地震干涉測量可以分離側(cè)面散射面波。

地震干涉測量理論研究已經(jīng)歷十多年發(fā)展，其研究和進(jìn)展為近地表散射波分離提供了新思路。地震干涉測量技術(shù)加強(qiáng)并預(yù)測近地表散射波，采用多道匹配濾波方法從原始地震記錄中自適應(yīng)減去近地表散射波。簡言之，近地表地震散射波分離就是利用近地表散射波運(yùn)動學(xué)規(guī)律，采用地震干涉測量技術(shù)，使近地表散射波得到相干加強(qiáng)，進(jìn)而分離近地表散射波。

3.1 地震干涉測量基本原理

地震干涉測量技術(shù)研究地震成對信號的干涉現(xiàn)象，即兩個接收點(diǎn)的地震記錄互相關(guān)，構(gòu)建一個新的地震響應(yīng)。該響應(yīng)可看作是其中一接收點(diǎn)上激發(fā)，另一接收點(diǎn)接收的響應(yīng)[9]，圖5所示即為地震干涉測量基本原理。

圖5 地震干涉測量基本原理示意圖

圖5a表示地表激發(fā)、地表接收的記錄，S為炮點(diǎn)，g1和g2為接收點(diǎn)，分別接收來自炮點(diǎn)S的信號，旅行時(shí)分別為t1和t2。將g1和g2接收的地震記錄互相關(guān)，得到其互相關(guān)記錄（見圖5b），其旅行時(shí)為t2－t1，此記錄相當(dāng)于在g1處放置震源S*，在g2處放置接收點(diǎn)g*。

互相關(guān)運(yùn)算實(shí)質(zhì)上是相位相減運(yùn)算，因此互相關(guān)型地震干涉測量消除了地震波從S到g1傳播過程的影響，而無需S與g1間介質(zhì)的參數(shù)（假設(shè)信號從S傳播到g2過程中經(jīng)歷了從S傳播到g1的相同路徑），無論其為均勻介質(zhì)還是不均勻介質(zhì)，也無論其為水平地表還是起伏地表，這就是地震干涉測量原理的核心思想。

3.2 散射波干涉測量原理分析

山地地震勘探中，由于山區(qū)地表起伏，因此震源和接收點(diǎn)并不在同一平面內(nèi)。圖6為觀測系統(tǒng)平面示意圖（俯視圖），其中S1—S7為震源在平面上投影位置，A和B為接收點(diǎn)在平面上投影位置，s為散射源在平面上投影位置。

圖6 近地表散射波分離中觀測系統(tǒng)配置示意圖

為便于討論，假設(shè)震源和接收點(diǎn)在同一直線上，且僅考慮地震波傳播旅行時(shí)。如圖7所示，從震源S到接收點(diǎn)B的直達(dá)波旅行時(shí)為tSB、到接收點(diǎn)A的直達(dá)波旅行時(shí)為tSA。當(dāng)A點(diǎn)直達(dá)波與B點(diǎn)直達(dá)波互相關(guān)時(shí)，得到從B到A直達(dá)波旅行時(shí)tBA（tBA= tSA－tSB）。A點(diǎn)直達(dá)波與B點(diǎn)直達(dá)波互相關(guān)，消除了地震波從S傳播到B這段共同路徑的影響。沿測線方向移動震源位置（向左或向右），如圖7所示直達(dá)波時(shí)距曲線也隨震源S移動。只要震源S在接收點(diǎn)B左側(cè)，沿測線方向移動震源，則直達(dá)波時(shí)距曲線會向下或向上移動，但從B點(diǎn)到A點(diǎn)的直達(dá)波旅行時(shí)tBA保持不變。于是將B點(diǎn)左側(cè)的所有炮在B點(diǎn)和A點(diǎn)上的直達(dá)波進(jìn)行互相關(guān)并疊加，則直達(dá)波將會相干加強(qiáng)。

圖7 直達(dá)波時(shí)距曲線示意圖

散射波分為兩種情況：①散射源在測線上，如圖8所示；②散射源不在測線上，如圖9所示（以直達(dá)散射波為例）。

圖8 直達(dá)波及其散射波時(shí)距曲線示意圖

圖9 直達(dá)波及其側(cè)面散射波時(shí)距曲線示意圖

當(dāng)散射源在測線上，如圖8所示，將B點(diǎn)接收的來自散射源s的散射直達(dá)波的旅行時(shí)記為tSsB，將A點(diǎn)接收的來自散射源s的散射直達(dá)波旅行時(shí)記為tSsA。B點(diǎn)直達(dá)波和A點(diǎn)直達(dá)波互相關(guān)得到的旅行時(shí)為tBA，tBA=tSA?tSB，B點(diǎn)直達(dá)波和A點(diǎn)散射波互相關(guān)得到的旅行時(shí)為tBsA，tBsA=tSsA?tSB。由此看出，B點(diǎn)直達(dá)波與A點(diǎn)直達(dá)波及其散射波互相關(guān)，A點(diǎn)直達(dá)波及其散射波旅行時(shí)都減少了tSB，消除了從震源S到接收點(diǎn)B這段共同傳播路徑的影響。

只要S位于接收點(diǎn)B左側(cè)，沿測線方向移動震源S，則直達(dá)波及其散射波時(shí)距曲線會平行向下或向上移動，而tBA和tBsA保持不變。將所有在B點(diǎn)左側(cè)激發(fā)、B點(diǎn)接收的直達(dá)波，分別與A點(diǎn)接收的直達(dá)波和散射直達(dá)波互相關(guān)，每一炮的互相關(guān)結(jié)果對應(yīng)相等，將其分別疊加就會相干加強(qiáng)。由此可知，近地表地震散射波分離正是利用了近地表散射波的運(yùn)動學(xué)規(guī)律，通過地震干涉測量，使近地表散射波相干加強(qiáng)。

當(dāng)散射源不在測線上，此時(shí)散射波為側(cè)面散射波。如圖9所示，s′為散射源在該測線上的垂直投影，黑線為來自震源S的直達(dá)波和來自散射源s的側(cè)面散射直達(dá)波的時(shí)距曲線，紅色虛線為當(dāng)震源S位于B點(diǎn)時(shí)接收的直達(dá)波和側(cè)面散射直達(dá)波的時(shí)距曲線。

依據(jù)近地表散射波運(yùn)動學(xué)規(guī)律，沿測線方向移動震源S，直達(dá)波時(shí)距曲線與散射波時(shí)距曲線分別相互平行。圖9中兩直達(dá)波時(shí)差為tSB，而散射波時(shí)距曲線間的時(shí)差不是tSB，而是tSs?tsB，其中tSs為波從震源S到散射源s的傳播時(shí)間，tsB為散射波從散射源s傳播到B點(diǎn)的時(shí)間。

對于側(cè)面散射源，將所有在B點(diǎn)左側(cè)激發(fā)，B點(diǎn)接收的直達(dá)波與A點(diǎn)接收的側(cè)面散射波（見圖9）互相關(guān)，不能保證每一炮的互相關(guān)得到的波形相位相同。但由波形疊加原理可知，若兩波形間的時(shí)差小于四分之一周期，則仍能形成同相疊加。因此當(dāng)散射源在測線附近時(shí)，只要滿足同相疊加條件，地震干涉測量也能加強(qiáng)側(cè)面散射波。

對多個散射源情況，以上分析仍成立，甚至對多次散射波也可類似分析。同理可以分析側(cè)面散射面波和側(cè)面散射折射波。采用二維或?qū)捑€地震資料，利用地震干涉測量法可分離測線附近的側(cè)面散射波，但不能完全分離所有側(cè)面散射波。

3.3 散射波分離的穩(wěn)相區(qū)分析

Halliday等[17]對震源穩(wěn)相區(qū)分布給出了很好的解釋，提出的地震干涉測量法分區(qū)采用互相關(guān)型和褶積型算法，有助于處理實(shí)際觀測系統(tǒng)地震數(shù)據(jù)。

圖10和圖11分別為一次散射波互相關(guān)型和褶積型干涉測量中震源和散射源穩(wěn)相區(qū)分布。圖10a和圖11a中震源S與接收點(diǎn)A、B及散射源s在同一直線上；圖10b和圖11b中震源S與散射源s在測線外。

圖10 互相關(guān)型干涉測量震源和散射源穩(wěn)相區(qū)分布

圖11 褶積型干涉測量震源和散射源穩(wěn)相區(qū)分布

對于互相關(guān)型干涉測量，震源S在接收點(diǎn)B的左邊，震源S到達(dá)接收點(diǎn)B的直達(dá)波與到達(dá)接收點(diǎn)A的散射波經(jīng)歷了同樣一段路徑（見圖10a），互相關(guān)得到的散射波就像是從B點(diǎn)激發(fā)經(jīng)過散射源s、最后A點(diǎn)接收的散射波，消除了相同傳播路徑的影響。

分析震源分布的穩(wěn)相區(qū)：當(dāng)散射源s在測線上（見圖10a），在震源穩(wěn)相區(qū)內(nèi)激發(fā)地震波，產(chǎn)生的地震記錄經(jīng)互相關(guān)型干涉測量能使散射波相干加強(qiáng)；當(dāng)散射源s在測線外（見圖10b），在震源穩(wěn)相區(qū)內(nèi)激發(fā)地震波，產(chǎn)生的地震記錄經(jīng)互相關(guān)型干涉測量也能使散射波相干加強(qiáng)。分析散射源分布的穩(wěn)相區(qū)：當(dāng)震源S在測線上（見圖10a），散射源穩(wěn)相區(qū)內(nèi)散射源s產(chǎn)生的散射波，經(jīng)互相關(guān)型干涉測量使散射波相干加強(qiáng)；當(dāng)震源S在測線外（見圖10b），互相關(guān)型干涉測量也能使散射源穩(wěn)相區(qū)內(nèi)散射源s產(chǎn)生的散射波相干加強(qiáng)。

同理可分析褶積型干涉測量的穩(wěn)相區(qū)（見圖11）。對于三維地震資料，地震干涉測量可分離散射波（包括部分側(cè)面散射波）。

4 近地表地震散射波分離效果

基于近地表散射波地震干涉測量理論，以中國西部K地區(qū)山地實(shí)際地震資料為例，進(jìn)行近地表散射波分離。

野外地震數(shù)據(jù)采集系統(tǒng)由3條炮線和3條接收線組成：炮線與接收線重合，道間距和線間距為30 m，炮間距為60 m，時(shí)間采樣間隔2 ms，地震記錄時(shí)間長度為8 s。

圖12為K地區(qū)礫石區(qū)單炮地震記錄，震源與接收點(diǎn)不在同一條測線上，炮點(diǎn)位于地形相對平坦的礫石區(qū)。從地震記錄中看出面波及其混響豐富，反向散射波和側(cè)面散射波清晰可見，而體波信息只在其間隱約可見。

圖12 礫石區(qū)單炮地震記錄

圖13為從單炮地震記錄中提取的直達(dá)波場，將其與全波場進(jìn)行地震干涉測量，得到直達(dá)波地震干涉測量炮記錄（見圖14）。直達(dá)波地震干涉測量加強(qiáng)了反向散射波和側(cè)面散射波能量，并預(yù)測了面波及其混響。雖然散射波分布在面波三角區(qū)，但是這些散射波是由于直達(dá)波和折射波在近地表傳播過程中遇到不均勻體而產(chǎn)生的，所以在直達(dá)波干涉測量中得以加強(qiáng)。

圖13 直達(dá)波場

圖14 直達(dá)波地震干涉測量炮記錄

分別采用常規(guī)噪聲壓制技術(shù)和近地表地震散射波壓制與常規(guī)噪聲壓制組合技術(shù)，進(jìn)行散射波分離。常規(guī)噪聲壓制技術(shù)（隨機(jī)噪聲和線性噪聲壓制技術(shù)）處理結(jié)果如圖15所示。近地表地震散射波壓制與常規(guī)噪聲壓制組合技術(shù)處理結(jié)果如圖16所示，其中近地表地震散射波壓制步驟為：①采用本文方法進(jìn)行近地表散射波分離；②利用多道匹配濾波技術(shù)，從原炮集中自適應(yīng)減去分離的近地表散射波。

圖15 常規(guī)噪聲壓制技術(shù)獲得的炮記錄

圖16 近地表地震散射波壓制與常規(guī)噪聲壓制組合技術(shù)獲得的炮記錄

兩種去噪方法都有效突出了體波信息，然而采用常規(guī)噪聲壓制處理，地震記錄中仍然存在明顯的反向散射波和側(cè)向散射波（見圖15），而采用近地表地震散射波壓制與常規(guī)噪聲壓制組合技術(shù)得到的地震剖面中基本消除了近地表散射波（見圖16）。

5 結(jié)論

地震勘探面臨的首要問題就是消除復(fù)雜近地表對地震波的影響，尤其在山地地震勘探中，粗糙起伏地形和不均勻近地表產(chǎn)生的散射波不僅影響地震資料的品質(zhì)，而且會在地震剖面上形成假構(gòu)造形態(tài)，因此有必要分離近地表散射波。

沿不均勻地表傳播的直達(dá)波、折射波和面波都會產(chǎn)生近地表散射波，在地震剖面上表現(xiàn)為線性同相軸和雙曲線同相軸，相鄰炮集上的近地表散射波同相軸相互平行。根據(jù)近地表散射波運(yùn)動學(xué)規(guī)律，采用地震干涉測量技術(shù)，使近地表散射波相干加強(qiáng)，實(shí)現(xiàn)近地表散射波分離。理論分析和實(shí)際資料應(yīng)用表明該方法不僅可分離線性散射波，而且可以分離部分側(cè)面散射波以及多次散射波，提高了地震資料品質(zhì)。該方法可適用于起伏地表和不均勻介質(zhì)，且不需要近地表速度信息。

山地地震資料處理面臨諸多問題，其中復(fù)雜近地表問題是研究難點(diǎn)之一，地震干涉測量技術(shù)為解決該問題提供了一種新思路。

[1]Blonk B,Herman G C,Drijkoningen G G.An elastodynamic inverse scattering method for removing scattered surface waves from field data[J].Geophysics,1995,60(6):1897-1905.

[2]Campman X H,Van Wijk K,Scales J A,et al.Imaging and suppressing near-receiver scattered surface waves[J].Geophysics,2005,70(2):V21-V29.

[3]Herman G C,Perkins C.Predictive removal of scattered noise[J].Geophysics,2006,71(2):V41-V49.

[4]徐基祥,王平,林蓓.地震波逆散射成像技術(shù)潛力展望[J].中國石油勘探,2006,11(4):61-66,79.Xu Jixiang,Wang Ping,Lin Pei.Expectation for technology potential of seismic wave inverse scattering image[J].China Petroleum Exploration,2006,11(4):61-66,79.

[5]徐基祥,楊旭明,胡英,等.一種在地震數(shù)據(jù)處理中衰減瑞雷波散射噪聲的方法:中國,ZL 2008 1 0110902.5[P].2011-01-26.Xu Jixiang,Yang Xuming,Hu Ying,et al.A method for rayleigh scattering noise attenuation in seismic data processing:China,ZL 2008 1 0110902.5[P].2011-01-26.

[6]吳希光,李亞林,張孟,等.復(fù)雜地區(qū)地震資料低信噪比的原因及對策（一）:表層散射波是導(dǎo)致地震資料低信噪比的根本原因[J].天然氣工業(yè),2012,32(1):27-32.Wu Xiguang,Li Yalin,Zhang Meng,et al.Causes of and solutions to a low signal-to-noise ratio of seismic data in complex areas,Part Ⅰ:Scattering of surface wave resulting in low signal-to-noise ratio of seismic data[J].Natural Gas Industry,2012,32(1):27-32.

[7]吳希光,李志榮,李正佳,等.復(fù)雜地區(qū)地震資料低信噪比的原因及對策（二）:表層散射波的基本性質(zhì)及形態(tài)特征[J].天然氣工業(yè),2012,32(2):38-42.Wu Xiguang,Li Zhirong,Li Zhengjia,et al.Causes of and solutions to a low signal-to-noise ratio of seismic data in complex areas,PartⅡ:Basic properties and shape features of surface scattered waves[J].Natural Gas Industry,2012,32(2):38-42.

[8]吳希光,李志榮,張曉斌,等.復(fù)雜地區(qū)地震資料低信噪比的原因及對策（三）:表層散射波壓制的技術(shù)方法[J].天然氣工業(yè),2012,32(3):57-61.Wu Xiguang,Li Zhirong,Zhang Xiaobin,et al.Causes of and solutions to a low signal-to-noise ratio of seismic data in complex areas,Part Ⅲ:Technical methods of surface scattered wave suppression[J].Natural Gas Industry,2012,32(3):57-61.

[9]Schuster G T.Seismic interferometry[M].New York:Cambridge University Press,2009.

[10]Wapenaar K,Ruigrok E,van Der Neut J,et al.Green’s function representation for seismic interferometry by deconvolution[C]// Society of Exploration Geophysicists.SEG technical program expanded abstracts.Denver:Society of Exploration Geophysicists,2010:3972-3978.

[11]Dong Shuqian,He Reiqing,Schuster G T.Interferometric prediction and least squares subtraction of surface waves[C]//Society of Exploration Geophysicists.SEG technical program expanded abstracts.New Orleans:Society of Exploration Geophysicists,2006:2783-2786.

[12]Halliday D,Curtis A.Seismic interferometry,surface waves and source distribution[J].Geophysical Journal International,2008,175(3):1067-1087.

[13]Forghani F,Snieder R.Underestimation of body waves and feasibility of surface-wave reconstruction by seismic interferometry[J].The Leading Edge,2010,29(7):790-794.

[14]Fleury C,Snieder R,Larner K.General representation theorem for perturbed media and application to Green’s function retrieval for scattering problems[J].Geophysical Journal International,2010,183(3):1648-1662.

[15]Snieder R,van Wijk K,Haney M,et al.Cancellation of spurious arrivals in Green’s function extraction and the generalized optical theorem[J].Physical Review E,2008,78(3):0366061-0366068.

[16]Snieder R,Fleury C.Cancellation of spurious arrivals in Green’s function retrieval of multiple scattered waves[J].The Journal of the Acoustical Society of America,2010,128(4):1598-1605.

[17]Halliday D F,Curtis A,Vermeer P,et al.Interferometric ground-roll removal:Attenuation of scattered surface waves in single-sensor data[J].Geophysics,2010,75(2):SA15-SA25.

（編輯 林敏捷 繪圖 劉方方）

A near-surface seismic scattered wave separation method

Xu Jixiang1,McLean B F2,Song Xuejuan1
(1.PetroChina Research Institute of Petroleum Exploration &Development,Beijing 100083,China;2.GeoScience Labs,Ltd,Beijing 100083,China)

Considering that near-surface scattered waves could not be suppressed by conventional denoise method,this paper analyzes the characteristics of scattered waves and achieves the near-surface scattered wave separation by the seismic interferometry theory.Seismic data in mountain area is generally characterized by low signal to noise ratio,weak effective wave energy and strong near-surface scattered wave energy,so the near-surface seismic scattered waves must be firstly separated to improve the signal to noise ratio of seismic data during seismic data processing.The events of the near-surface seismic scattered wave in the adjacent common shot gathers are parallel.According to the near-surface scattered wave kinematics,the seismic interferometry method is used to eliminate the influence of seismic wave propagation in the same path,coherently enhance near-surface scattered waves,and then separate them from original seismic records.Taking the seismic data in mountain area in west China as an example,the energy of near-surface scattered wave is enhanced by direct wave field and full wave field seismic interferometry to separate the waves.As a result,the seismic data after processing generally eliminates near-surface scattered waves,verifying the validity and feasibility of the near-surface seismic scattered wave separation method.

mountain seismic exploration;near-surface;seismic scattered wave;scattered wave separation;seismic interferometry

國家科技重大專項(xiàng)（2011ZX05003-003）；國家重大產(chǎn)業(yè)技術(shù)開發(fā)專項(xiàng)（發(fā)改辦高技[2005]2372號）；中國石油勘探開發(fā)研究院院級項(xiàng)目（2011Y-006）

TE122

：A

1000-0747(2014)06-0705-07

10.11698/PED.2014.06.08

徐基祥（1966-），男，安徽東至人，博士，現(xiàn)為中國石油勘探開發(fā)研究院高級工程師，主要從事地震勘探資料成像技術(shù)研究工作。地址：北京市海淀區(qū)學(xué)院路20號，中國石油勘探開發(fā)研究院物探技術(shù)研究所，郵政編碼：100083。E-mail：xjx6611@petrochina.com.cn

2013-12-09

2014-10-12