王漢闖，陳生昌＊，陳國新，梁東輝，佘德平

1 浙江大學地球科學系，杭州 310027

2 中國石油化工股份有限公司石油物探技術研究院，南京 210014

1 引言

采集效率與處理效率的提高是地震采集技術一直追求的目標.傳統(tǒng)觀測系統(tǒng)的設計受一定的放炮時間間隔門檻值的制約，施工周期都較長；為了縮短采樣周期，在觀測系統(tǒng)設計的時候總是盡可能地減少激發(fā)炮數，從而導致激發(fā)震源點不足、資料品質差.近年來國外興起的多震源激發(fā)技術（Bagaini，2006；Beasley et al.，1998；Moerig et al.，2002）則可以克服這個缺陷，大大地提高觀測效率，縮短采樣周期，從而降低勘探成本，并已經在近幾年地震數據采集中逐步得到應用（Aaron et al.，2009；Beasley，2008；Berkhout et al.，2009；Berkhout et al.，2008；Blacquiere et al.，2009；Hampson et al.，2008），在地震勘探中表現出了高效、低成本等巨大的優(yōu)越性.多震源地震方法在海洋與陸地高效地震數據采集中都開始得到應用，所使用的震源主要為可控震源（Bagaini，2006；Moerig et al，2002）.陸上多震源采集如獨立同步掃描方法（Howe et al.，2009；Howe et al.，2008）（ISS）、遠距離同步掃描方法（Bouska，2010）（DS3）等一般用于三維大面積地震數據采集，多臺震源車間隔一定的距離獨立工作；為了減少彼此之間的干擾，各個震源之間一般間隔較遠，極大地提高了生產率并一定程度上壓制了震源之間的干擾.相比陸地勘探，海洋多震源地震方法（Beasley，2008；Hampson et al，2008；Ting and Zhao，2009）應用更加廣泛，震源主要使用氣槍震源，放置震源的船按照一定的排列進行隨機激發(fā)，放置拖纜檢波器束的船同時接收來自各個震源的波場數據.這種工作方式極大地提高了工作效率，并可以進行寬方位角地震數據采集.所謂多震源就是在多個不同位置上激發(fā)震源，可以是同步激發(fā)也可以是非同步激發(fā)（如線性時間延遲、隨機激發(fā)等），因此得到的地震記錄是一種由多個震源產生的混合地震記錄，如何對這種多震源地震數據進行處理是亟待研究的問題.

對于多震源地震數據，若可以準確地分離為常規(guī)單震源地震數據，就可以按照現有的方法進行成像處理.多個震源激發(fā)的情況下，對某一個震源的波場來講，來自其他震源的波場可以認為是噪音.因此可以把多震源地震數據的分離問題看作一個去噪問題來研究，通過幾何濾波（Beasley，2008；Beasley et al，1998）、變換域濾波（Akerberg et al.，2008）、多方向矢量中值濾波（Huo et al.，2009）等方法可以達到分離多震源地震數據的目的；也可以把多震源數據的采集過程看成一個數學過程，把多震源地震數據的分離問題轉化為一個反演問題來研究，采取反演（Ikelle，2007）與迭代反演的方法（Mahdad et al.，2011）求取單震源地震數據；多次掃描（含一次掃描）多震源地震數據的分離問題（Wang et al.，2013）也得到了深入討論.多震源地震數據分離后就得到了各個震源的地震數據，然后就可以按照常規(guī)單震源地震偏移成像方法進行處理.

在目前的單震源地震勘探方法技術中，波動方程疊前深度偏移在實際數據處理中獲得了越來越多的認同，共炮點道集數據疊前深度偏移成像可以獲得高質量的偏移成像結果，但因邊界處理而增加計算量等原因致使它的計算效率很低.該方法的計算效率還與炮道集數有關，炮道集數越多，計算效率越低（Berkhout，2008）.針對波動方程共炮點道集偏移成像的上述特點，人們提出了一系列改進方法，如平面波偏移成像方法（Berkhout，1992；Chen et al.，2012；Rietveld et al.，1992；Zhang et al.，1999；陳樹民等，2001）、相位編碼（Romero，et al.，2000；Sun，et al.，2002）方法與小波束方法（Wu and Chen，2001）等.陳生昌根據地震數據的線性疊加原理提出了地震數據廣義合成方法（陳生昌和馬在田，2006），又把合成函數選取為隨機函數，發(fā)展為廣義隨機合成偏移成像方法（陳生昌等，2012）.按照這種方法，人們可以根據不同的地質情況和要求通過人工合成，得到各種不同的人工合成震源和地震數據道集，如平面波震源和數據道集、局部平面波震源和道集以及面向目標的人工合成震源和數據道集，對這些人工合成的震源與地震數據可以直接進行偏移成像處理，大大提高了地震數據偏移成像的計算效率.這些都為發(fā)展更高效、更準確、更具針對性的波動方程疊前深度偏移成像方法技術奠定了基礎.與廣義合成震源類似，多震源就是在多個不同位置上激發(fā)震源，是一種按照一定格式編碼的多個震源，所形成的波場是多個震源波場疊加的混合波場；不同的是廣義合成震源與廣義合成地震數據是通過人工手段對常規(guī)單炮地震數據所做的處理，目的是為了提高室內數據處理效率，而多震源是按照一定的采集系統(tǒng)要求設計的，多震源地震數據是在新的采集系統(tǒng)下直接采集得到的.由此可見，多震源數據的偏移成像可以采取與廣義合成地震數據相似的處理方法.

本文一方面把多震源地震數據分離問題作為一個數學反演問題來研究，通過多震源地震數據分離方法把多震源地震數據分離為常規(guī)單震源地震數據，然后按照常規(guī)單程波偏移成像方法進行偏移成像處理；另一方面在陳生昌提出的廣義隨機合成道集偏移成像方法的基礎上，結合Berkhout等提出的Blended Source和Blended Data及其偏移成像的概念（Berkhout，2008；Berkhout et al.，2009），把多震源看作經過編碼的廣義震源，多震源地震數據也就是這種廣義震源激發(fā)條件下的地震數據，提出了多震源地震數據直接偏移成像方法.把本文提出的兩種處理多震源數據的成像方法應用于數值模擬的多震源地震數據中，并與常規(guī)單震源地震方法成像結果進行對比，結果表明本文方法是正確的、有效的.

2 多震源地震波場的傳播與模擬

2.1 多震源地震波場的傳播

二維密度均一的頻率域聲波方程可以表示為

這里，v（x，z）為位置r（x，z）處的聲波速度，u（x，z，ω）為位置r處的位移場或壓力場，rs（xs，zs）為震源位置，S（rs，ω）為震源項.

多震源是用多個震源在不同的位置上隨機延遲激發(fā)的組合（即經過編碼的廣義震源），表達式可以寫為

這里，Γ（rs，ω）為地震波場混合矩陣算子：

其中t（rs）是位置rs上各個震源的隨機延遲激發(fā)時間，α（rs）為與炮點坐標有關的權系數.為了保持多震源中各個震源函數的振幅一致性，一般令權系數α（rs）＝1，有

由式（1）和（2）可知，多震源激發(fā)條件下的波動方程可以表示為

設rd（xd，zd）為檢波器位置，由式（1）和（5）可知，多震源Sbl（rs，ω）激發(fā)條件下檢波器接收到地震記錄Pbl＝ubl（xd，zd，ω）和常規(guī)單震源激發(fā)條件下的地震記錄P＝u（xd，zd，ω）之間的關系就可以表示為可以看出，多震源地震數據是常規(guī)單震源地震數據的線性疊加，疊加算子就是多震源的混合矩陣算子.

同時，由式（1）和（5）也可以看出，描述常規(guī)單震源和多震源激發(fā)的地震波場的波動方程在形式上是一致的，區(qū)別主要是兩者右端的震源項.把震源項改為多震源，方程就可以描述多震源激發(fā)條件下的地震波場的傳播規(guī)律.

2.2 多震源地震波場的模擬

地震波場數值模擬是研究復雜介質中地震波傳播規(guī)律的有效途徑之一，這里在前一小節(jié)多震源地震波場傳播理論的基礎上進行多震源波場的數值模擬.多震源數值模擬的觀測系統(tǒng)在常規(guī)觀測系統(tǒng)基礎上重新設計，如圖1所示，設觀測排列上總的震源個數為nshot，把整個震源排列分為N部分，每部分的震源個數為M（則nshot＝N×M）；每個部分的中第i（i＝1，2，…，M）個震源一起組成一組一起激發(fā)的多震源（即震源1、M＋1、2M＋1、…、（N－1）M＋1等N個震源組成一個多震源，震源2、M＋2、2M＋2、…、（N－1）M＋2等N個震源組成一個多震源…，共有M對組合）.檢波器的分布范圍固定，即在測區(qū)范圍固定好檢波器，然后移動炮點進行激發(fā)而不再移動檢波點，第一個檢波點與第一震源點的距離為xl.這樣，在所有M個多震源激發(fā)觀測之后就可以得到M個多震源地震數據，也就完成了整個測區(qū)的觀測.與常規(guī)單震源地震采集系統(tǒng)相比，多震源地震采集系統(tǒng)一次激發(fā)多個炮點，明顯加快了采集效率.

圖1 多震源地震觀測系統(tǒng).整個觀測區(qū)域分成N部分，每個部分有M個震源，每個部分第i個一起組成一個多震源（i＝1，2，…，M），共M個多震源.xl為第一個震源與第一個檢波點之間的距離，檢波點位置固定不變Fig.1 Multi－source seismic acquisition system.Partition the whole acquisition region to Nsubintervals，each of which contains Msources.The i－th source of each subinterval（i＝1，2，…，M）form a blended sources，and Mblended sources have been formed in total by this way.xlis the distance from the 1st source to the 1st receiver，and all the receivers are fixed in the acquisition system

由式（2）—（3）及圖1所示觀測系統(tǒng)可以看出，影響多震源的因素主要有3個：（1）多震源中激發(fā)震源的個數；（2）多震源中各個激發(fā)震源位置；（3）多震源中各個震源的隨機延遲激發(fā)時間，設置合適的多震源參數對多震源模擬波場的特征有著重要影響.多震源地震波場數值模擬主要是改變震源函數，即使用多個震源同步激發(fā)，然后按照波動方程進行傳播，最后得到多震源地震數據.

如圖2a所示層狀速度模型，從上到下層速度依次為3000、3500、4000m／s，在圖1所示的觀測系統(tǒng)中，整個測區(qū)（0～4000m）設置400個震源、400個檢波器，檢波器排列的范圍和震源排列的范圍相同（xl＝0），炮間距和道間距都為10m，令N＝2，M＝200，即每次使用兩個震源一起激發(fā)，共可獲得200個多震源地震記錄.以第101個多震源（1000m處的s1震源點和3000m處的s2震源點）為例，兩個震源的延遲激發(fā)時間分別為368ms和789ms，記錄時間點為4001，長度為4s，時間采樣率為1ms，得到的模擬地震記錄如圖2b所示.可以看出，多震源地震記錄是來自兩個單震源激發(fā)的混合地震記錄，同相軸相互交織在一起，形成一種復雜的地震記錄.

圖2 層狀速度模型及其所對應的含兩個震源激發(fā)的地震記錄（a）層狀速度模型，從上到下層速度依次為3000m／s、3500m／s、4000m／s；（b）含兩個震源激發(fā)的混合地震記錄，兩個震源的延遲激發(fā)時間分別為368ms和789ms.Fig.2 The layered velocity model and one of the blended shot gathers with 2simultaneous sources（a）The layered velocity model with the velocity of the three layers are 3000m／s，3500m／s and 4000m／s from the top to the bottom；（b）One of the blended shot gathers with 2simultaneous sources whose random time delays are 368ms and 789ms respectively.

3 多震源數據分離方法

多震源地震波場中，對于欲得到的那個震源波場而言，其他震源產生的波場可以看作是噪音，若去除了這些噪音也就實現了多震源地震數據的分離.下面結合2.2節(jié)模擬得到的200個含兩個震源的多震源地震記錄（其中一炮如圖2b所示）為例說明多震源地震數據的分離方法.

把式（6）寫為離散形式，則多震源地震數據可以表示為

式中，k＝1，2，…，N，每個單震源波場Pk隨機延遲激發(fā)得到的多震源地震波場Pbl是通過混合矩陣Γk對各個單震源波場Pk做時移然后再疊加.

由于地震波場滿足線性疊加原理，多震源波場中的交叉噪音在其本來的觀測道集中是顯現不出來的.考慮到每個震源波場都具有一定的隨機延遲時間，可以反時移消除Pbl中施加在各個單震源波場上的隨機延遲時間，這一步可通過混合矩陣Γk預分離處理，即有

其中，是Γk的共軛轉置，從上式可知：預分離結果不僅包含了要求的Pk也包含了其他震源波場Pi（i≠k）（交叉噪音波場），但相對于Pk消除了隨機延遲激發(fā)時間.如圖2所示模型試驗中，200個含兩個震源激發(fā)的地震記錄數據經過預分離得到了400個單震源地震記錄，每個記錄都消除了延遲時間，處于同一反射層的同相軸頂點（圖3a中的點A和圖3b中的點B）處于同一時間線上.

由于多個震源激發(fā)的波場包含在采集到的多震源地震數據中，且隨機延遲激發(fā)時震源編碼不是正交的，通過預分離處理的交叉噪音波場Pi（i≠k）在其他的人工分選道集（如共檢波點道集、共中心點道集和共偏移距道集等）中并不滿足線性疊加原理，就明顯地表現出來了（Mahdad et al，2011），隨機延遲激發(fā)方式產生的交叉噪音變?yōu)殡S機噪音.這里把分選到共檢波點域，即預分離得到的400個地震記錄分選到共檢波點域，對某一欲分離數據來說則來自其他震源的波場都作為隨機噪音的形式出現在道集中，如圖3c、圖3d所示；

在人工分選的道集中對預分離波場進行去噪處理，任何能用來消除隨機噪音的濾波方法都可以應用.在人工分選的道集中進行去噪處理后，再進行抽道集就可得到消除了交叉噪音波場Pi（i≠k）后的預分離波場P′k＝－Pi（i≠k）.如果中交叉噪音波場Pi（i≠k）能得到完全消除，則P′k就是待求的第k號震源的單震源波場Pk，即得到了分離波場P.這里用中值濾波進行隨機噪音去除，對于圖2b所示的多震源記錄，分離之后得到如圖3e、圖3f所示的兩個單震源地震記錄.

4 多震源地震數據偏移成像

在多震源地震波場傳播理論、數值模擬及數據分離的基礎上，針對多震源地震數據的特點，本文提出了兩種多震源地震數據的偏移成像方法：分離后再偏移成像及直接偏移成像處理.

圖3 圖2b所示的含兩個震源激發(fā)多震源地震數據分離圖示（a）和（b）是預分離結果；（c）和（d）是分選到共檢波點域中的結果；（e）和（f）是最終分離結果.圖中各圖都去除了隨機延遲激發(fā)時間，時間變?yōu)?s.Fig.3 The separation of the blended shot gather with 2simultaneous sources that shown in Fig.2b（a）and（b）are the pseudo－debelending results，and that in common receiver domain are shown in（c）and（d），（e）and（f）are the final separation results.The random time delays of all the seismic shot gathers are eliminated，and the record length become 3s.

4.1 分離后再偏移成像

完成了多震源地震數據的分離處理之后，就可以按照常規(guī)偏移成像方法進行處理.設震源深度為z0，檢波器深度為zd，深度方向上用下行波方程正向外推多震源波場，即利用波場算子（陳生昌等，2001）W（zn，z0）將震源波場從z0傳播到zn；同時用上行波方程反向外推多震源記錄波場，即利用W＊（zn，zd）將檢波器波場從zd傳播到zn，也就獲得了深度zn上的震源波場S（x，zn，ω）和檢波器波場P（x，zn，ω），這個過程可以表示為

其中＊表示復共軛，應用時間一致性成像條件提取偏移成像結果，則深度zn處的像可以由如下公式給出：

式中，x表示橫向位置，zn表示深度，ω表示頻率，Re表示取實部運算，I（x，zn）表示 （x，zn）處的像.

4.2 直接偏移成像

多震源激發(fā)條件下，多震源地震波場滿足波動方程（5），因此可采用與常規(guī)單震源激發(fā)地震偏移成像相似的技術方案，差別在于震源波場進行傳播時，震源項是多個單震源的組合，記錄波場反向傳播時，使用的是混合地震記錄，即來自各個震源的波場記錄.

設多個震源的深度為z0，并令式（2）的深度為z0，即rs＝ （x，z0），據此就可以得到多震源地面激發(fā)的波場Sbl，記為Sbl（x，z0，ω）.這個過程可以寫為

其中，Γ（x，z0，ω）＝α（x，z0，ω）ejωt（x，z0），表示地面處多震源波場的混合算子.把式（11）表示的混合震源取代方程（9）的第一個式子中的震源S（x，z0，ω）進行震源波場傳播，即

由此就可以得到各個深度層zn上的多震源波場Sbl（x，zn，ω）.

設檢波器深度為zd，則記錄波場Pbl可以寫為Pbl（x，zd，ω）.用Pbl（x，zd，ω）取代方程（9）中第二個式子的記錄波場P（x，zd，ω）進行檢波器波場的反向傳播，即使用單程波外推算子的共軛W＊（zn，zd）將檢波器波場從zd傳播到zn，可得深度zn上的檢波器波場Pbl（x，zn，ω）：

得到了深度zn上的震源波場Sbl（x，zn，ω）和檢波器波場Pbl（x，zn，ω）之后，根據時間一致性成像條件，深度zn處的像可以用如下公式給出：

式中，x表示橫向位置，ω表示頻率，＊表示復共軛，Re表示取實部運算，Ibl（x，zn）表示 （x，zn）處的像.經過深度外推就可以得到各個深度的成像結果，最后輸出整個成像剖面.

由以上可知，多震源成像和常規(guī)單震源成像類似，只要分別在多個震源加上隨機激發(fā)時間的影響，即通過混合矩陣算子把多個位置上激發(fā)的單震源組合為多震源，然后通過波場傳播得到地下各個區(qū)域的多震源波場，結合多震源地震記錄波場，應用時間一致性成像原理可以得到最終的成像結果.

下面從理論上分析多震源激發(fā)對成像結果的影響.在式（7）中，Pk可以寫為

這里Lk是第k個震源對應的傳播模擬算子.把上式代入式（7）可得

其中，m表示地下真實模型.偏移成像算子ΦT就可以定義如下的形式：

則偏移成像結果就是

從式（18）就可以看出，第一項是來自各個單震源的波場的成像結果，而第二項就是多個震源之間震源波場與記錄波場間不匹配而產生的交叉干擾像（Dai and Schuster，2009）.若Γi和Γj存在比較大的差異，即震源i和j距離較遠或隨機延遲時間差別大，則可以通過疊加和偏移很好地壓制交叉干擾像.

由于偏移結果中包含了各個震源波場的貢獻，因此從理論上說和常規(guī)單震源偏移結果是可以對比的，若不考慮交叉干擾像則直接成像結果和常規(guī)單震源成像結果的幅值是一致的.野外實際生產過程中設計觀測系統(tǒng)時盡可能把隨機延遲時間加大一些、同步激發(fā)的震源距離遠一些，這樣就可以明顯降低交叉干擾像，通過偏移成像和疊加又壓制了其中部分隨機干擾，直接偏移成像可以取得和常規(guī)偏移成像方法相差不多的效果.

5 數值試驗

為了驗證上節(jié)提出的兩種多震源地震數據偏移成像方法的正確性和有效性，本文設計圖4a所示的速度模型進行試驗.該模型是中國石化石油物探技術研究院根據實際地質資料建模得到的，模型的淺層中部有一個類似地層尖滅的斷層；模型深層是一個背斜構造，利于油氣富集與存儲.該模型是對實際地質模型的一個簡化，可以一定程度上反映斷層和背斜構造在油氣勘探中的影響.模型整體呈層狀構造，淺層速度?。?000m／s），深層速度大（4700m／s）.

模型網格點數：nx＝3200、nz＝220，網格間距：dx＝dz＝10m.試驗中的多震源觀測系統(tǒng)如圖1所示，單炮數nshot＝250，炮間距為80m，共1600道，道間距為20m，第一個檢波點與第一個震源點的距離xl＝6000m，每一個混合炮的觀測范圍都等于常規(guī)單炮激發(fā)總的觀測范圍.這里取一起激發(fā)震源個數N＝2，得到125個多震源地震記錄，記錄長度為8s，時間采樣率為1ms；震源使用主頻為30Hz的雷克子波.按照上述觀測系統(tǒng)進行數值模擬，得到125個含兩個震源的多震源地震記錄，其中一個如圖4b所示.

5.1 分離后再偏移成像試驗

本文對上述125個含兩個震源的多震源地震數據按照文中所述方法進行分離處理，得到了250個單震源地震記錄.分離后，各個記錄波場都消除了隨機延遲激發(fā)時間的影響，記錄長度減少為6s.如圖4b所示的多震源地震記錄分離結果見圖5，可以看到本文方法很好地實現了多震源地震記錄的分離，得到了相對應的單震源地震數據，隨機延遲激發(fā)引起的交叉噪音得到了很好的壓制.對分離之后的地震數據進行偏移成像處理，成像結果如圖6a所示，可以看出分離后再偏移成像由于疊前偏移成像過程中的疊加，震源隨機延遲激發(fā)引起的分離波場中的隨機噪音在偏移成像處理中得到了很好的壓制，去除了部分隨機噪音，整體剖面連續(xù)性好、構造清晰.

為了便于比較，這里對常規(guī)單震源地震方法也進行了試驗，觀測范圍和多震源觀測范圍相同，單炮數nshot＝250，震源同樣使用主頻為30Hz的雷克子波，記錄長度為6s，時間采樣率為1ms，圖6b為常規(guī)單震源偏移成像結果，可以看出，多震源分離后再成像結果相比常規(guī)單震源偏移成像結果，模型吻合度、構造連續(xù)性等方面基本一致，證明本文提出的多震源地震數據分離后再偏移成像方法是正確可行的.

圖5 圖4b所示多震源地震記錄經過分離得到的兩個單震源地震記錄Fig.5 The two separated seismic shot gathers from the blended shot gather that shown in Fig.4b

圖6 二維斷層－背斜構造模型偏移成像結果（a）125個含兩個震源的多震源地震記錄分離得到250個地震記錄后再偏移成像結果；（b）常規(guī)單震源地震記錄偏移成像結果.Fig.6 The migration results of 2Dfault－anticline model（a）The migration result of the 250separated shot gathers from 125blended shot gathers with 2simultaneous sources；（b）The conventional migration result by single－source method.

另外，為了驗證激發(fā)震源個數對偏移成像的影響，本文在圖1觀測系統(tǒng)的基礎上，使用5個震源一起激發(fā)，得到了50個多震源地震記錄，同樣完成了對整個測區(qū)的觀測，其中一個地震記錄如圖7a所示（從左至右震源的延遲時間分別為：1407、304、1114、1982和1778，單位ms），可以看出多個單震源記錄經過延遲疊加在一起，形成了更為復雜的地震記錄.圖7b是按照本文方法分離出來的對應第2個震源的地震記錄，與兩個震源激發(fā)地震記錄的分離結果相比，所需要的主要信息已經呈現出來，但有部分殘余干擾存在于分離后的記錄中；圖7c是相應觀測條件下只在第二個震源位置激發(fā)的常規(guī)單炮地震記錄，可以看出分離結果還有部分殘留，但主要同相軸等信息已經分離出來了.分離后再偏移成像結果如圖7d所示，可以看到成像過程的疊加作用壓制了分離結果中的部分隨機殘留干擾，剖面連續(xù)性比較好、位置信息比較準確；但也出現了一些大的假象和干擾，主要集中在淺層（圖中圈出部分），這主要是由于分離結果中的殘留波場所致；但因為疊前偏移的疊加過程，淺層由于地震波傳播距離較小，在成像點上的疊加次數少，而在深層地震波傳播距離較大，在成像點上的疊加次數多，這種疊加作用壓制了大部分的殘留波場.

在主頻1.6GHz的Intel（R）Xeon（R）CPU、內存8G的計算機上進行多震源地震記錄分離并進行偏移成像處理，兩種觀測系統(tǒng)下多震源記錄的分離時間大體一樣為75s，偏移成像時間為28773s，總的計算時間為28848s.

5.2 直接偏移成像試驗

為了驗證本文所述直接偏移成像方法對多震源地震數據的成像效果，這里對數值模擬得到的125個含兩個震源的多震源地震數據直接進行偏移成像處理，結果如圖8a所示，可以看到這種情況下的直接偏移成像結果比圖6a所示的分離后再偏移成像結果信噪比更高，與圖6b所示的常規(guī)單震源偏移成像精度相當.這說明當同步激發(fā)震源個數較少（這里為2）、激發(fā)震源相距較遠時，多震源地震直接成像效果可以取得與單震源地震方法成像相當的精度.

圖7 二維斷層－背斜構造模型多震源分離后再偏移成像試驗圖示（a）一個含5個震源激發(fā)的多震源地震記錄；（b）圖（a）的分離結果中的第二個震源對應的地震記錄；（c）和圖b相對應的常規(guī)單震源地震記錄；（d）分離得到250個單震源記錄的成像結果.Fig.7 The migration result after the blended seismic data separation of 2Dfault－anticline model（a）One of the blended shot gathers with 5simultaneous sources；（b）The 2nd shot gather of the separation results from the data that shown in Fig.a；（c）The conventional single－source shot gather on the location of the 2nd source corresponding to Fig.b；（d）The migration result of the 250separated shot gathers.

類似地，這里也對如圖7a所示的50個含5震源激發(fā)的多震源地震記錄進行了直接偏移成像試驗，結果如圖8b所示.從成像結果中可以看出直接偏移成像方法對模型的主要構造取得了很好的成像結果，但由于同步激發(fā)震源個數的增多以及隨機延遲激發(fā)時間的影響，整體剖面有一些隨機干擾；同樣由于淺層地震波在成像點上的疊加次數少，深層地震波在成像點上的疊加次數多的原因，這種干擾在淺層尤為明顯，深層干擾較小、效果較好；整體成像剖面也會出現一些大的假象，在成像剖面中分布比較分散（圖中圈出部分）.因此，多震源勘探中必須合理地設置觀測系統(tǒng)，盡可能使同步激發(fā)的震源相隔地遠一些，降低相互之間的干擾，提高成像精度.為了更清楚地比較直接偏移成像的交叉噪音干擾，這里取成像結果中2070～2130網格上的成像結果進行對比分析，道集顯示如圖9所示，（a）為兩個震源激發(fā)成像結果的道集顯示，（b）為5個震源激發(fā)成像結果的道集顯示，可以清楚地看到隨著同步激發(fā)震源個數的增多，成像結果中淺層的干擾也隨著增多，深層也有一定的變化，但干擾像明顯弱于淺層的.

同樣在主頻1.6GHz的Intel（R）Xeon（R）CPU、內存8G的計算機上進行多震源地震記錄分離并進行偏移成像處理，對125個含兩個震源的多震源地震數據直接偏移成像處理計算時間為19057s，對50個含5個震源的多震源地震數據直接偏移成像處理計算時間為7344s，相比分離后再偏移成像的計算時間28848s都提高了計算效率，也說明隨著震源數目的增加，在提高野外采集效率、降低數據采集成本的同時，室內處理時間明顯地減少，降低了處理成本.

5.3 復雜模型試驗及振幅保真性分析

Marmousi模型是法國石油研究院推出的模型數據，是檢驗波動方程疊前深度偏移成像算法標準模型數據.該模型構造復雜，可以用來驗證本文方法對復雜模型的適應性.

圖8 二維斷層－背斜構造模型直接偏移成像結果（250個單震源，不同的激發(fā)方式）（a）125個含兩個震源的多震源地震記錄成像結果；（b）50個含5個震源的多震源地震記錄成像結果.Fig.8 The direct migration result of 2Dfault－anticline model（different blended shooting methods with a total number of 250sources）（a）The direct migration result from 125shot gathers with 2simultaneous sources；（b）The direct migration result from 50shot gathers with 5simultaneous sources.

圖9 圖8所示的直接偏移成像結果中第2040—2140網格位置的道顯示（a）兩個震源同步激發(fā)的情況；（b）5個震源同步激發(fā)的情況.Fig.9 The traces from 2040to 2140of direct migration result that shown in Fig.8（a）The condition of blended shooting with 2sources；（b）The condition of blended shooting with 5sources.

Marmousi速度模型網格大小為nx＝243，nz＝750，網格間距dx＝25，dz＝4（單位：m）.本試驗仍利用如圖1所示的多震源觀測系統(tǒng)，取單震源個數nshot＝240，炮間距為25m，接收道數為241道，道間距為25m，第一個檢波點與第一個震源點的距離xl為25m，每一個混合炮的觀測范圍都一樣，且等于常規(guī)單炮激發(fā)總的觀測范圍.震源使用主頻為20Hz的雷克子波，一起激發(fā)的震源個數N＝2，激發(fā)次數M＝120，記錄長度為4s，時間采樣率為4 ms.按照上述觀測系統(tǒng)進行數值模擬，得到120個含兩個震源的多震源地震記錄，其中一個（第38個）混合地震記錄如圖10b所示，按照本文所述方法進行分離處理，得到了240個單震源所對應的地震記錄，圖10b所示的混合記錄對應的分離結果如圖10c和圖10d所示.可以看到，分離結果中消除了相應的延遲時間，得到了單震源記錄的主要信息；由于波場記錄的復雜程度較高，分離結果中還有部分殘留.

對通過分離120個混合地震記錄得到的240個地震記錄進行偏移成像處理，結果如圖11a所示.為了便于對比分析，本文還對常規(guī)單震源模擬得到的240炮地震數據進行了偏移成像處理，結果如圖11b所示.可以看到，分離后再偏移成像的結果構造清晰，和常規(guī)地震方法相差不大，本文提出的分離后再進行偏移成像的方法對于復雜模型也是適應的.

應用本文所述直接偏移成像方法，對如圖10的120個含兩個震源的混合地震記錄進行了直接偏移成像處理，結果如圖12所示.對比圖11可以看出，直接偏移成像方法對于復雜模型的成像效果也是不錯的，對陡傾等復雜構造都可以很好地成像，本文提出的直接偏移成像方法對復雜模型也是適應的.不足之處是偏移結果中出現了一些隨機噪音，信噪比有待進一步提高.

圖11 Marmousi模型多震源數據分離后偏移與常規(guī)單震源偏移對比試驗（a）分離得到的240炮地震數據偏移結果；（b）常規(guī)單震源激發(fā)得到的240炮地震數據偏移結果.Fig.11 The comparison of migration after separation and conventional single－source method（a）The migration result from 240separated shot gathers；（b）The migration result from 240shot gathers of conventional single－source method.

為了驗證多震源激發(fā)成像方法的振幅保真性，這里以圖11b和圖12所示的偏移結果為例，從以下三點對常規(guī)單震源成像方法和直接偏移成像方法的效果進行對比分析：（1）取偏移成像結果的2450m位置上的一列數據進行分析.兩種方法的幅值變化曲線如圖13a所示，其中藍線表示常規(guī)成像方法成像結果的各個深度的幅值，紅線表示直接成像方法成像結果的各個深度的幅值（下同）.可以看出兩種方法的成像結果數值上有很好的擬合，在一些小的變化劇烈的地方（如圖中A、B、C和D點等處）的擬合會有稍微差異.（2）取480m深度層上x方向各個位置的數據進行分析.兩種方法的幅值變化曲線如圖13b所示，可以看出兩種方法的成像結果擬合的很好，在振幅突變點（如圖中A、B點）也有比較好的效果；直接偏移成像方法在某些地方（如圖中C點）會有一定的震蕩出現.（3）取1640m深度層上的各個x方向的數據進行分析.兩種方法的幅值變化曲線如圖13c所示，與淺層（上一種情況）相比，兩種成像方法的結果整體上擬合的更好，也一定程度上從側面說明了深層的成像效果要優(yōu)于淺層.

圖12 Marmousi模型直接偏移成像結果（120個含兩個震源的地震記錄）Fig.12 The direct migration result of Marmousi model（120blended shot gathers with 2simultaneous sources）

通過以上對比分析可知，本文提出的兩種對多震源混合地震數據的處理方法對復雜構造模型也是適應的，直接成像方法在一定程度上相比常規(guī)單震源地震成像方法是保幅的.

圖13 Marmousi模型常規(guī)成像方法和直接成像方法結果幅值對比曲線（a）2450m橫向位置上的數據幅值對比曲線；（b）480m深度位置上的數據幅值對比曲線；（c）1640m深度位置上的數據幅值對比曲線.其中藍線表示常規(guī)成像方法，紅線表示直接成像方法.Fig.13 The amplitude curves of the conventional migration result and the direct migration result of Marmousi model（a）The curves on the horizontal position of 2450m；（b）The curves on the depth of 480m；（c）The curves on the depth of 1640m.The blue curve represents the conventional migration result and the red one represents the direct migration result.

6 結論

多震源地震勘探方法技術是當今地震勘探的一個熱點方向，主要是為了提高勘探效率、降低勘探成本，是一種有著良好發(fā)展前景的方法技術，特別是對于三維大面積探區(qū)，同步激發(fā)震源個數可以設置的更多，道集數也就成倍減少，后續(xù)的處理效率也就會有大幅度提升.本文從多震源的本質出發(fā)，基于多震源波場的傳播規(guī)律以及波場的線性疊加原理，提出了多震源數據分離后再偏移成像和多震源地震數據的直接偏移成像方法，兩種方法在數值模擬數據的模型試驗中均得到了滿意的成像效果，表明本文方法是正確的、有效的，可以應用于多震源地震數據的偏移成像處理流程中.主要結論有以下幾點：

（1）當多震源中激發(fā)震源較少（或相距較遠）時，來自各個震源的波場相互干擾相對較弱，分離后再偏移成像和直接偏移成像方法都可以取得和常規(guī)單震源地震方法相當的效果，總體精度比較高，這說明本文提出的多震源地震數據的分離后再偏移成像方法與直接偏移方法可以滿足疊前深度偏移要求.

（2）當多震源中激發(fā)震源較多（或相距較近）時，深層由于疊前深度偏移中疊加次數增多的原因，兩種方法偏移成像效果都呈現深層較淺層好的特點；多震源地震記錄的分離剖面上殘留的一些隨機干擾在偏移成像后已經得到了一定程度的壓制，但仍有部分假象出現在淺層；直接偏移成像剖面中假象分布比較分散.

（3）多震源地震數據的直接偏移成像方法相比常規(guī)單震源地震成像方法是保幅的，特別是在同步震源距離較遠時，交叉干擾像會比較弱，成像的保幅性更好.

總的來說，直接偏移成像省去了分離多震源地震數據的步驟，計算效率相比常規(guī)偏移成像和分離后再偏移成像有了很大程度的提高；特別是隨著多震源中激發(fā)震源個數增多，偏移成像的計算效率會更高，但偏移成像效果會有所下降.實際工作中可以把兩種成像方法結合起來，直接進行偏移成像處理仍是多震源地震勘探數據處理的一個重要的發(fā)展方向；本文是對多震源激發(fā)條件的混合波場成像問題進行一些探索研究，多震源地震方法還處于發(fā)展階段，速度建模方法還有待研究，這也是我們下一步的工作方向.

Aaron P，van Borselen R，Fromyr E.2009.Simultaneous Sources：A Controlled Experiment On Different Source Configurations.79th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，56－60.

Akerberg P，Hampson G，Rickett J，et al.2008.Simultaneous source separation by sparse Radon transform.78th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，2801－2805.

Bagaini C.2006.Overview of simultaneous vibroseis acquisition methods.76th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，70－74.

Beasley C J.2008.A new look at marine simultaneous sources.The LeadingEdge，27（7）：914－917.

Beasley C J，Chambers R E，Jiang Z J.1998.A new look at simultaneous sources.68th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，133－135.

Berkhout A J.1992.Areal shot record technology.Journalof SeismicExploration，1（3）：251－264.

Berkhout A J.2008.Changing the mindset in seismic data acquisition.TheLeadingEdge，27（7）：924－938.

Berkhout A J，Blacquière G，Verschuur D J.2009a.The concept of double blending：Combining incoherent shooting with incoherent sensing.Geophysics，74（4）：59－62.

Berkhout A J，Blacquière G G，Verschuur E G.2008.From simultaneous shooting to blended acquisition.78th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，2831－2838.

Berkhout A J，Verschuur D J，Blacquière G.2009b.Seismic imaging with incoherent wavefields.79th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，2894－2898.

Blacquiere G，Berkhout G，Verschuur E.2009.Survey design for blended acquisition.79th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，56－60.

Bouska J.2010.Distance separated simultaneous sweeping，for fast，clean，vibroseis acquisition.GeophysicalProspecting，58（1）：123－153.

Chen S C，Cao J Z，Ma Z T.2001.Prestack depth migration method based on quasi－linear born approximation.Chinese JournalofGeophysics－ChineseEdition（in Chinese），44（5）：704－710.

Chen S C，Ma Z T.2006.Generalized synthesis of seismic data and its migration.ChineseJournalofGeophysics－ChineseEdition（in Chinese），49（4）：1144－1149.

Chen S C，Wang H C，She D P.2012.Migration of generalized random synthesis of seismic data.OilGeophysicalProspecting（in Chinese），47（6）：868－872.

Chen S M，Liu H，Li Y M.2001.The area shot method for land 3D seismic acquisition.ProgressinGeophysics（in Chinese），16（3）：58－67.

Dai W，Schuster J.2009.Least－squares migration of simultaneous sources data with a deblurring filter.79th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，2090－2994.

Hampson G，Stefani J，Herkenhoff F.2008.Acquisition using simultaneous sources.TheLeadingEdge，27（7）：918－923.

Howe D，Foster M，Allen T，et al.2009.Independent simultaneous sweeping in Libya－full scale implementation and new developments.79th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，109－111.

Huo S D，Luo Y，Kelamis P.2009.Simultaneous sources separation via multi－directional vector－median filter.79th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，31－35.

Ikelle L.2007.Coding And Decoding：Seismic data modeling acquisition and processing.77th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，66－70.

Howe D，Foser M，Allen T，et al.2008.Independent simultaneous sweeping－a method to increase the productivity of land seismic crews.78th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，2826－2830.

Mahdad A，Doulgeris P，Blacquiere G.2011.Separation of blended data by iterative estimation and subtraction of blending interference noise.Geophysics，76（3）：9－17.

Moerig R，Barr F J，Nyland D L.2002.Simultaneous shooting using cascaded sweeps.72th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，74－76.

Rietveld W E A，Berkhout A J，Wapenaar C P A.1992.Optimum seismic illumination of hydrocarbon reservoirs.Geophysics，57（10）：1334－1345.

Romero L A，Ghiglia D C，Ober C C，et al.2000.Phase encoding of shot records in prestack migration.Geophysics，65（2）：426－436.

Sun P Y，Zhang S L，Liu F.2002.Prestack migration of areal shot records with phase encoding.72nd Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，1172－1175.

Ting C O，Zhao W.2009.A simulated wide azimuth simultaneous shooting experiment.79th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，76－80.

Wang H C，Chen S C，Zhang B，et al.2013.Separation method for multi－source blended seismic data.AppliedGeophysics，10（3）：251－264.

Wu R S，Chen L.2001.Beamlet migration using Gabor－Daubechies frame propagator.63rd EAGE Annual Meeting，Extended Abstracts，1356－1359.

Zhang G，Zhang W S，Hao X J.1999.Prestack depth migration with common－shot and synthesis－shot records.69th Ann.Internat Mtg.，Soc.Expi.Geophys.，Expanded Abstracts，1469－1472.

附中文參考文獻

陳樹民，劉洪，李幼銘.2001.適于窄線三維地震資料的面炮方法.地球物理學進展，16（3）：58－67.

陳生昌，王漢闖，佘德平.2012.地震數據廣義隨機合成的偏移成像.石油地球物理勘探，47（6）：868－872.

陳生昌，馬在田.2006.廣義地震數據合成及其偏移成像.地球物理學報，49（4）：1144－1149.

陳生昌，曹景忠，馬在田.2001.基于擬線性Born近似的疊前深度偏移方法.地球物理學報，44（5）：704－710.