葉月明，趙昌壘，莊錫進，郭書娟，錢進

1中國石油杭州地質(zhì)研究院，杭州 310023

2中國石化石油物探技術(shù)研究院，南京 211103

3中國科學(xué)院海洋研究所，青島 266071

1 引言

地震波成像的目標(biāo)就是將地表接收到的反射波反向歸位至產(chǎn)生它的反射點上，并由此得到地下構(gòu)造特征.多數(shù)地震偏移成像算法是基于所有散射能量在地下界面僅反射一次的假設(shè)，而實際上，聲波在向地表返回過程中會有部分發(fā)生二次或多次反射，形成多次波.由于海水面存在著強的波阻抗差異，所以多次波在海洋地震資料中尤為嚴(yán)重.多次波的出現(xiàn)經(jīng)常使地震剖面上主要目的層反射難以解釋，特別是沿潛在油氣儲層的微小變化將被多次波嚴(yán)重干擾，所以，長期以來多次波一直被視為噪音并從地震數(shù)據(jù)中去除（郭科等，2010；石穎和邢小林，2011；陳小宏和劉華鋒，2012；石穎等，2012；趙昌壘等，2013），以免在之后的地震資料解釋中造成誤解.

事實上，多次波也是眾多地震信號中的一種，像一次波一樣，多次波也包含了地下界面信息.在復(fù)雜地下構(gòu)造的情況下，有的區(qū)域一次反射波難以照明，而多次波卻可以通過多次反射而照明，包含了一次波所沒有的地下信息，所以，近些年來，國內(nèi)外許多專家學(xué)者都在利用多次波方面做了一系列的研究.

將多次波轉(zhuǎn)變?yōu)橐淮尾ㄔ侔凑諅鹘y(tǒng)偏移方法進行成像是多次波成像方法之一，Shan等 （2003）提出了將含有多次波的數(shù)據(jù)進行互相關(guān)變換為準(zhǔn)一次波，基于雙平方根方程，在炮點-檢波點域進行多次波疊前深度偏移.Shan和Guitton（2004）等通過一次波與多次波的互相關(guān)得到準(zhǔn)一次波，在炮集實現(xiàn)了多次波疊前深度偏移，并應(yīng)用Sigsbee2B模型進行了測試，取得了一定的效果.He等 （2007）提出了地震干涉成像法對自由表面多次波進行成像，并在VSP數(shù)據(jù)中得到了較好的應(yīng)用.Jiang（2006）利用層間多次波對VSP資料進行成像，前提是需要有較強產(chǎn)生層間多次波的界面才能夠較好的對層間多次波進行識別和成像.Vasconcelos等（2008）也利用地震干涉法對墨西哥灣鹽下構(gòu)造進行VSP數(shù)據(jù)的目標(biāo)成像，成像結(jié)果與傳統(tǒng)地面地震成像結(jié)果一致.Berkhout和Verschuur（2006）提出了可以先通過SRME預(yù)測出地表相關(guān)多次波，再將多次波轉(zhuǎn)變?yōu)橐淮尾?，然后再進行多次波成像.（Muijs and Holliger 2005，Muijs etal.，2007）首先將波場分解為上行波場與下行波場，然后再進行基于單程波的下行波場和上行波場延拓，最后通過多維成像得到成像結(jié)果，這種方法的缺點是在前期需要進行波場分離.Brown和Guitton（2005）提出了一種多次波和一次波聯(lián)合成像的方法，稱為多次波和一次波最小二乘法聯(lián)合成像，這種方法利用模型正則化算子將一次波與多次波分離，結(jié)合一次波與多次波的信息來提高信號的保真度.

近年來，隨著計算機硬件的發(fā)展，逆時偏移得到了一定程度的應(yīng)用，Liu等（2011）將逆時偏移震源波場用地震記錄來替代、記錄波場用預(yù)測出的多次波來替代實現(xiàn)了基于逆時偏移的多次波成像，取得了較好的效果，然而由于逆時偏移本身的特點，這種計算方法也是比較耗時的.Zhang和Schuster（2014）提出了基于最小二乘的多次波逆時偏移方法，通過多次迭代運算，改善了逆時多次波偏移方法的成像質(zhì)量，在照明均衡和縱向成像分辨率方面都有了明顯的提高.Wang等（2014）提出了一種不需要預(yù)先預(yù)測多次波的方法，以脈沖震源與原始記錄為下行延拓初始波場，以原始記錄為初始上行波場，避免了多次波預(yù)測所增加的誤差.郭書娟等（2011，2012）提出了基于廣義的炮偏移方法實現(xiàn)地表多次波和一次波的聯(lián)合成像，該方法是將包含有震源的記錄波場做為下行波場，多次波為上行波場來實現(xiàn)一次波與多次波的聯(lián)合成像，由于震源子波包含在記錄波場中，如果子波能量太小，一次波成像比重就較小，如果子波能量太強，多次波成像比重就掩蓋了一次波成像，所以在確定震源子波能量上存在一定的困難.

在前人研究的基礎(chǔ)上，本文系統(tǒng)闡述了多次波在基于單程波疊前深度偏移中的影響，并提出了基于單程波偏移算子的多次波成像方法.該方法在輸入數(shù)據(jù)上采取了多次波逆時偏移類似的做法，將輸入的下行波場用包含多次波的記錄來替代，輸入的上行波場用預(yù)測的表層相關(guān)多次波來替代，采用基于波動方程的單程波偏移算子來實現(xiàn)多次波成像，具有較高的計算效率.為了應(yīng)用多次波成像結(jié)果來彌補一次波成像的不足，通過基于二范式的能量匹配原則將多次波與一次波成像結(jié)果匹配疊加，實現(xiàn)一次波與多次波的聯(lián)合成像，避免了廣義炮偏移聯(lián)合成像中子波能量難以確定的缺點的同時也具有較高的計算效率.簡單模型驗證了多次波成像方法的正確性，Sigsbee2B模型進行的一次波與多次波聯(lián)合成像結(jié)果取得了明顯的改善，尤其是多次波發(fā)育的鹽丘邊界構(gòu)造成像質(zhì)量得到了明顯提高.

2 多次波在單程波疊前深度偏移中的影響

在密度恒定的各向同性完全彈性介質(zhì)中，假設(shè)地震波震源是t＝0時刻激發(fā)的脈沖，則地震波的傳播可以用如下時間-空間域的二維標(biāo)量聲波方程表示為

其中，U（x，z）為空間波場，v（x，z）為介質(zhì)的速度.

基于擬微分算子理論，假設(shè)光滑介質(zhì)中層間多次反射非常弱，或者認(rèn)為多次波被視為噪聲已經(jīng)提前濾除，不考慮地震波的反射和透射得到解耦的單程波方程在頻率域表示為（張關(guān)泉，1993；Zhang etal.，2005）

其中，UU（x，z；w）與UD（x，z；w）分別代表了分裂出的下行波場與上行波場，w為頻率，算子Λ＝稱為擬微分算子，從（2）式就可以推導(dǎo)出分步傅里葉與傅里葉有限差分等經(jīng)典的單程波的波場延拓算子（Stoffa etal.，1990；Ristow and Rūhl，1994）.

圖1是多次波傳播的路徑圖，其中US是震源波場，在地下R處發(fā)生第一次反射，UP是接收到的一次反射波場，一次反射波場在地表被接收后又被反射回地下，在R1處產(chǎn)生二次反射波，返回地表形成一階多次波UM1，類似地產(chǎn)生一系列多次波，最大階為N階多次波UMN.在常規(guī)疊前深度偏移中，輸入下行波場是震源子波，輸入的上行波場是地表接收到的地震記錄，一般認(rèn)為多次波已經(jīng)在前期預(yù)處理中被衰減掉，按照Claerbout教授的入射波的初至與反射波的產(chǎn)生時間相同的成像原則（Claerbout，1971），當(dāng)不存在多次波時，單程波偏移方法成像條件可寫為

圖1 多次反射波傳播路徑Fig.1 Propagation paths of reflective multiples

其中，wmin與wmax代表最小與最大頻率，U＊S（x，z；w）為下行震源波場在頻率域的復(fù)共軛，UP（x，z；w）為一次反射波，為簡便起見，以下公式中波場的表示都省略 （x，z；w），當(dāng)存在多次反射波時，單程波偏移方法實際成像為

其中，第Ⅰ項為一次反射波得到的真實構(gòu)造成像部分，第Ⅱ項為多次波與震源波場互相關(guān)而產(chǎn)生的構(gòu)造假象，由于多次波一般成周期性規(guī)律出現(xiàn)，所以第Ⅱ項低階多次波部分與震源互相關(guān)成像結(jié)果會在有效構(gòu)造區(qū)域產(chǎn)生構(gòu)造假象，影響解釋人員對構(gòu)造的判斷；由于高階多次波產(chǎn)生時間較長，高階多次波與震源互相關(guān)而形成的構(gòu)造假象一般出現(xiàn)在較深部位，容易判斷.由上述可知，在基于單程波的疊前深度偏移之前，應(yīng)該盡可能多的在疊前炮集處理中去除多次波.

3 基于單程波偏移算子的多次波成像

基于波動方程的單程波疊前深度偏移方法相對于射線類的Kirchhoff積分法偏移可以適應(yīng)速度場的強橫向變化，相對于全波逆時偏移具有較高的計算效率和較低的速度依賴性.本文采取類似于劉伊克（2011）提出的逆時偏移多次波成像方法的邊界條件，首先采用地表相關(guān)多次波預(yù)測的方法，預(yù)測出地表相關(guān)多次波；其次用含有多次波的地震記錄作為單程波偏移中的輸入下行波場UD，用預(yù)測出的地表相關(guān)多次波作為單程波偏移中的輸入上行波場UU；最后用式（5）互相關(guān)成像條件成像，公式為

在時間域相關(guān)的兩個信號，在頻率域是一個信號與另一個信號復(fù)共軛的標(biāo)量積.式（5）中的復(fù)共軛項代表了下行波場的正向外推，非復(fù)共軛項代表了上行波場的反向外推.整個多次波成像由三部分構(gòu)成，第一部分Ⅲ項是多次波真實構(gòu)造成像部分，認(rèn)為地表接收到的第N-1階多次波UM（N-1）為震源波場，經(jīng)地層界面RN反射后，在地表接收到的第N階多次波為該源產(chǎn)生的記錄波場UMN，上、下行波場通過互相關(guān)成像得到RN位置的成像值，根據(jù)第Ⅲ項，依次可以得到R1，R2…，RN-1位置處的成像；第二部分Ⅳ項是多次波偏移中產(chǎn)生的噪音，成像條件中的UM（N-2）和UMN分別是下行波場與上行波場，雖然能夠成像，但是所成的像是偏移假象或噪音，因為在UM（N-2）與UMN之間還有一次波場的往返，而單程波疊前深度偏移并沒有考慮到這種往返，還是按照單程波場延拓的方式來延拓波場，在比實際反射層更深的地層處相關(guān)成像，產(chǎn)生偏移假象或噪音；第三部分Ⅴ是不能夠產(chǎn)生成像值的，以UMN·U＊MN為例，因為下行波場UMN的正向延拓意味著波前面時間在變長，而上行波場UMN的反向延拓意味著波前面時間在變小，兩者隨著波場延拓是不可能相遇并相關(guān)成像的，由于沒有偏移噪音產(chǎn)生，這部分對偏移成像影響不大.

4 最小二范式一次波與多次波成像匹配疊加

對于構(gòu)造成像來說，應(yīng)該是最大限度的利用地

下信息來獲取真實構(gòu)造特征.廣義一次波與多次波聯(lián)合成像方法（郭書娟等，2011）將同時含有表層多次波的炮記錄與脈沖震源之和作為下行延拓波場，但是，脈沖震源的能量與炮記錄存在能量差，如果脈沖震源能量太強，聯(lián)合成像結(jié)果主要表現(xiàn)為一次波對成像的貢獻(xiàn)；如果脈沖震源能量太弱，那么聯(lián)合成像的結(jié)果主要表現(xiàn)為多次波對成像的貢獻(xiàn)，需要權(quán)衡脈沖子波的能量是該方法的一個缺陷.本節(jié)采取的思路是先分別得到一次波成像結(jié)果Io（x，z）與多次波成像結(jié)果Im（x，z），然后再基于最小二乘能量匹配原則，求取匹配因子，進行匹配疊加，一次波與多次波成像結(jié)果的優(yōu)勢得到互補，實現(xiàn)聯(lián)合成像.這樣做既利用多次波成像結(jié)果補充了一次波對一些構(gòu)造照明的不足，又避免了廣義一次波多次波聯(lián)合成像方法中震源子波在輸入波場中能量難以確定的缺點.最小誤差能量匹配疊加的關(guān)鍵因素在于匹配因子的求取，本文采取基于最小二范式的匹配濾波方法來求取匹配因子，使得多次波成像結(jié)果與一次波成像結(jié)果數(shù)量級一致，聯(lián)合成像結(jié)果Iom（x，z）可以表示為

其中，a（x，n）是匹配因子，N是信號的長度，匹配因子由式（7）的Toeplitz矩陣決定，可以通過Levinson快速遞推算法求解為

其中，K是匹配因子的長度，K＜N.ψmm（x，z）是多次波成像結(jié)果Im（x，z）的自相關(guān)，ψom（x，z）是一次波成像結(jié)果與多次波成像結(jié)果的互相關(guān)，公式為

5 模型試算

為了驗證本文方法的有效性，本文首先對三層模型進行了基于波動方程的一次波疊前深度偏移試處理和多次波成像，最后對Sigsbee2b模型進行了一次波與多次波聯(lián)合成像試處理，取得了一定的效果.

5.1 三層簡單模型試算

該三層模型的速度場如圖3a所示，第一層、第二層和第三層速度分別是1500m·s-1、2500m·s-1和3500m·s-1.圖2是正演模擬單炮記錄，中間放炮，151道接收，道間距20m，記錄時間3s，采樣率4ms.其中圖2a是含有多次波的單炮記錄，D1是直達(dá)波，R1和R2分別是第一層和第二層的一次反射波，能量較弱的I1反射同相軸是第一層與第二層間的一階層間多次波，M11和M21分別是第一層界面和第二層界面對應(yīng)的地表相關(guān)一階多次波，P1是微屈多次波，M11和M21分別是第一層界面和第二層界面對應(yīng)的地表相關(guān)二階多次波.去除出多次波后得到包含直達(dá)波的一次反射波記錄如圖2b所示，預(yù)測出的地表相關(guān)多次波如圖2c所示.

圖3b是應(yīng)用去除多次波后的記錄（圖2b）做的單程波疊前深度偏移，第一反射層界面L1與第二個反射層界面L2都得到了正確的成像，而且偏移噪音很小.應(yīng)用包含多次波炮集數(shù)據(jù)進行單程波疊前深度偏移的結(jié)果如圖3c所示，雖然可以正確成像L1與L2層界面，但是，由于多次波的存在產(chǎn)生了嚴(yán)重的偏移假象，其中，I是由震源子波延拓的下行波場與層間多次波延拓的上行波場I1產(chǎn)生的偏移假象；M1與M2分別是震源延拓的下行波場與第一層界面和第二層界面對應(yīng)的地表相關(guān)一階多次波M11和M12延拓的上行波場產(chǎn)生的偏移假象.圖3d是應(yīng)用本文方法得到的多次波成像結(jié)果，其中，N1是層間多次波I1與第一個反射層的一階地表相關(guān)多次波M11得到的偏移假象；N2是第二層的反射波R2與第一反射層的一階地表相關(guān)多次波M11得到的偏移假象；L1是實際反射層的成像，它的成像貢獻(xiàn)來自多個上、下行延拓波場的相關(guān)疊加，包括第一層的反射波R1與第一層的一階地表相關(guān)多次波M11的相關(guān)成像、第二層的反射波R2與第二層的一階地表相關(guān)多次波M21的相關(guān)成像、第一層的一階地表相關(guān)多次波與第一層的二階地表相關(guān)多次波的相關(guān)成像、第二層的一階地表相關(guān)多次波與第二層的二階地表相關(guān)多次波的相關(guān)成像等，也就是式（5）中的第Ⅲ部分.類似地，L2層的成像由第二層的一次反射R2與第二層的一階多次波、第二層的一階多次波與第二層的二階多次波等相關(guān)成像得到，這就是為什么多次波相對于一次波具有很廣照明范圍的原因，明顯的表現(xiàn)就是圖3d中L1層成像范圍大廣于一次波對L1層成像（圖3c）的原因.N3、N4和N5都是偏移假象，如式（5）中的第Ⅳ部分所示.

圖2 三層模型正演記錄（a）包含多次波的單炮記錄；（b）一次波記錄；（c）多次波記錄.Fig.2 Forward shot gathers of a three-Layer velocity model（a）Single shot gather with multiples；（b）Primaries；（c）Multiples.

圖3 簡單模型和偏移結(jié)果（a）簡單模型速度場；（b）一次波疊前深度偏移；（c）包含多次波的記錄疊前深度偏移；（d）多次波成像.Fig.3 Simple velocity model and migration（a）Simple velocity model；（b）Migration of primaries；（c）Migration with records including primaries and multiples；（d）Imaging with multiples.

圖4 Sigsbee2B模型速度場與單炮記錄（a）速度-深度模型；（b）包含多次波的單炮記錄；（c）去除地表相關(guān)多次波后的單炮記錄；（d）地表相關(guān)多次波.Fig.4 Sigsbee2Bvelocity model and single shot record（a）Velocity-depth model；（b）Single shot record including multiples；（c）Single shot record after multiples elimination；（d）Surface-related multiples.

圖5 Sigsbee2B模型疊前深度偏移結(jié)果（a）一次波疊前深度偏移；（b）多次波疊前深度偏移；（c）一次波與多次波成像結(jié)果匹配疊加.Fig.5 Pre-stack depth migration results of Sigsbee2Bmodel（a）Pre-stack depth migration with primaries；（b）Pre-stack depth migration with multiples；（c）Matched stack of pre-stack depth migration of primaries and multiples.

圖6 匹配濾波因子Fig.6 Filter operator for matching

圖7 Sigsbee2B模型鹽丘邊界處疊前深度偏移（a）速度模型；（b）一次波疊前深度偏移；（c）多次波疊前深度偏移；（d）一次波與多次波成像結(jié)果匹配疊加.Fig.7 Salt boundary pre-stack depth migration of Sigsbee2Bmodel（a）Velocity model；（b）Pre-stack depth migration of primaries；（c）Pre-stack depth migration of multiples；（d）Matched stack of pre-stack depth migration of primaries and multiples.

5.2 Sigsbee2B模型試算

為了測試本方法對復(fù)雜模型的適用性，下面對國際上通用的多次波數(shù)據(jù)處理測試模型，Sigsbee2B進行試算.該模型的主要特點是存在起伏的海底和高陡鹽丘構(gòu)造（圖4a），圖4b是含有多次波的記錄，圖4c是不含有地表相關(guān)多次波的記錄，圖4d是地表相關(guān)多次波記錄，共3495炮，每炮348道，記錄長度12s.圖5a是對炮集數(shù)據(jù)基于單程波偏移算子傳統(tǒng)疊前深度偏移結(jié)果，整體構(gòu)造得到了較好的成像，但是在鹽下及鹽丘邊界處由于初至反射波照明不足以及鹽丘的遮擋屏蔽作用，使得成像質(zhì)量差，成像能量弱.圖5b是基于本文的方法得到的多次波成像結(jié)果，由于地表相關(guān)多次波主要發(fā)育在具有強速度差的地方，也就是鹽丘構(gòu)造存在的地方，所以多次波成像的主要貢獻(xiàn)就是在鹽丘邊界處，多次波成像使得鹽丘邊界得到了較好的刻畫.對于不存在鹽丘的地方（圖4a左邊部分），不存在強速度差，地表相關(guān)多次波不發(fā)育，所以多次波成像能量非常弱（圖5b左邊部分），這也說明了只有在多次波發(fā)育的地方，多次波對構(gòu)造成像的貢獻(xiàn)才比較大.為了利用多次波對鹽下構(gòu)造成像的優(yōu)勢，通過匹配疊加的方式，進行一次波與多次波聯(lián)合成像.圖6是該匹配濾波因子，匹配因子的長度是49，匹配疊加后的結(jié)果如圖5c所示.圖7b、圖7c和圖7d分別是傳統(tǒng)單程波偏移、多次波成像和聯(lián)合成像結(jié)果在鹽下高陡構(gòu)造處的局部放大，箭頭所指地方結(jié)果多次波成像后得到了明顯改善和提高.

6 結(jié)論

基于單程波偏移算子的疊前深度偏移方法可以適應(yīng)速度場強橫向變化的同時具有較高的計算效率，本文對多次波在疊前深度偏移中的影響做了詳細(xì)的分析，修改了單程波疊前深度偏移的初始條件，地震記錄替代震源子波，預(yù)測得到的地表相關(guān)多次作為上行波場而實現(xiàn)了基于單程波偏移算子的地表相關(guān)多次波成像，發(fā)揮了多次波能夠?qū)Φ叵聵?gòu)造多次照明的優(yōu)勢.使用了基于最小二范式一次波與多次波成像匹配疊加的方法，使得多次波的成像結(jié)果較好的彌補一次波成像的不足，尤其是高陡鹽下構(gòu)造邊界這些一次波照明不足的地方成像效果得到了顯著的改善.該項技術(shù)研究使我們意識到多次波這類“噪音”也是攜帶地下構(gòu)造信息的信號，地球物理工作者可充分發(fā)揮其對地下構(gòu)造多次照明的優(yōu)勢，最大限度的提取地下信息.

Berkhout A J，Verschuur D J.2006.Imaging of multiple reflections.Geophysics，71（4）：SI209-SI220，doi：10.1190／1.2215359.

Brown M P，Guitton A.2005.Least-squares joint imaging of multiples and primaries.Geophysics，70（5）：S79-S89，doi：10.1190／1.2052471.

Chen X H，Liu H F.2012.Comparison between inverse scattering series method and SRME method in free surface related multiple prediction.Progress in Geophysics（in Chinese），27（3）：1040-1050，doi：10.6038／j.issn.1004-2903.2012.03.026.

Claerbout J F.1971.Toward a unified theory of reflector mapping.Geophysics，36（3）：467-481.

Guo K，Guo S，He G Z，etal.2010.A method of multiple blind separation based on ICA.Progress in Geophysics（in Chinese），25（3）：1075-1080，10.3969／j.issn.1004-2903.2010.03.048.

Guo S J，Li Z C，Tong Z Q，etal.2011.Joint imaging of primaries and surface-related multiples based on generalized shot-profile migration.Chinese J.Geophys.（in Chinese），54（4）：1098-1105.doi：10.3969／j.issn.0001-5733.2011.04025.

Guo S J，Li Z C，Tong Z Q，etal.2012.Method and technique for imaging of surface-related multiples.Progress in Geophysics（in Chinese），27（6）：2570-2576，doi：10.6038／j.issn.1004-2903.2012.06.034.

He R Q，Hornby B，Schuster G.2007.3Dwave-equation interferometric migration of VSP free-surface multiples.Geophysics，72（5）：S195-S203，doi：10.1190／1.2743375.

Jiang Z Y.2006.Migration of interbed multiple reflections.∥76th Ann.Internat.Mtg.，Soc.Expl.Geophys.Expanded Abstracts，3501-3505.

Liu Y K，Chang X.2011.Reverse time migration of multiples.∥81st Ann.Internat.Mtg.，Soc.Expl.Geophys.Expanded Abstracts，3326-3331.

Liu Y K，Chang X，Jin D G，etal.2011.Reverse time migration of multiples for subsalt imaging.Geophysics，76（5）：WB209-WB216，doi：10.1190／GEO2010-0312.1.

Muijs R，Holliger K.2005.Prestack depth migration of primary and surface-related multiple reflections.∥75th Ann.Internat.Mtg.，Soc.Expl.Geophys.Expanded Abstracts，2107-2110.

Muijs R，Robertsson J O，Holliger K.2007.Prestack depth migration of primary and surface-related multiple reflections：Part I-Imaging.Geophysics，72（2）：S59-S69，doi：10.1190／1.2422796.

Ristow D，Rūhl T.1994.Fourier finite-difference migration.Geophysics，59（12）：1882-1893.

Shan G J，Guitton A.2004.Migration of surface-related multiples：tests on the Sigsbee2Bdataset.∥74th Ann.Internat.Mtg.，Soc.Expl.Geophys.Expanded Abstracts，1285-1288.

Shan G.2003.Source-receiver migration of multiple reflections.∥73rd Ann.Internat.Mtg.，Soc.Expl.Geophys..Expanded Abstracts，1008-1011.

Shi Y，Jing H L，Li Y.2012.Surface-related multiple suppression effect analysis by feedback iteration approach.Progress in Geophysics（in Chinese），27（4）：1493-1500，doi：10.6038／j.issn.1004-2903.2012.04.024.

Shi Y，Xing X L.2011.Investigation progress on surface-related multiple suppression：review and outlook.Progress in Geophysics（in Chinese），26（6）：2046-2054，doi：10.3969／j.issn.1004-2903.2011.06.020.

Stoffa P L，F(xiàn)okkema J T，De Luna Freire R M，etal.1990.Splitstep Fourier migration.Geophysics，55（4）：410-421.

Vasconcelos I，Snieder R，Hornby B.2008.Imaging internal multiples from subsalt VSP data-examples of target-oriented interferometry.Geophysics，73（4）：S157-S168，doi：10.1190／1.2944168.

Wang Y B，Chang X，Hu H.2014.Simultaneous reverse time migration of primaries and free-surface related multiples without multiple prediction.Geophysics，79（1）：S1-S9，doi：10.1190／GEO2012-0450.1.

Zhang D L，Schuster G T.2014.Least-squares reverse time migration of multiples.Geophysics，79（1）：S11-S21，doi：10.1190／GEO2013-0156.1.

Zhang G Q.1993.System of coupled equations for upcoming and downgoing waves.Acta Mathematicae Applicatae Sinica（in Chinese），16（2）：251-263.

Zhang Y，Zhang G Q，Bleistein N.2005.Theory of true-amplitude one-way wave equations and true-amplitude common-shot migration.Geophysics，70（4）：E1-E10，doi：10.1190／1.1988182.

Zhao C L，Ye Y M，Yao G S，etal.2013.Prediction deconvolution in linear radon domain on the application of ocean multiples attenuation.Progress in Geophysics（in Chinese），28（2）：1026-1032，doi：10.6038／pg20130256.

附中文參考文獻(xiàn)

陳小宏，劉華鋒.2012.預(yù)測多次波的逆散射級數(shù)方法與SRME方法及比較.地球物理學(xué)進展，27（3）：1040-1050.

郭科，郭思，何國柱等.2010.基于獨立分量分析的多次波盲分離方法研究.地球物理學(xué)進展，25（3）：1075-1080.

郭書娟，李振春，仝兆岐等.2011.基于廣義的炮偏移方法實現(xiàn)地表多次波和一次波聯(lián)合成像.地球物理學(xué)報，54（4）：1098-1105.

郭書娟，李振春，仝兆岐等.2012.表層多次波成像方法技術(shù)研究.地球物理學(xué)進展，27（6）：2570-2576.

石穎，井洪亮，李瑩.2012.反饋迭代法壓制表面多次波效果分析.地球物理學(xué)進展，27（4）：1493-1500.

石穎，邢小林.2011.表面多次波壓制的研究進展：回顧與展望.地球物理學(xué)進展，26（6）：2046-2054.

張關(guān)泉.1993.波動方程的上行波和下行波的耦合方程組.應(yīng)用數(shù)學(xué)學(xué)報，16（2）：251-263.

趙昌壘，葉月明，姚根順等.2013.線性拉東域預(yù)測反褶積在海洋多次波去除中的應(yīng)用.地球物理學(xué)進展，28（2）：1026-1032.