張建利，劉志斌，張?jiān)迄i，李維新，王 赟

（1．中海油研究總院，北京100028；2．中國(guó)科學(xué)院地質(zhì)與地球物理研究所，北京100029；3．中國(guó)地質(zhì)大學(xué)地球物理與信息技術(shù)學(xué)院，北京100083）

兩種海底多分量波場(chǎng)分離方法的模型測(cè)試及改進(jìn)

張建利1，2，劉志斌1，張?jiān)迄i1，李維新1，王 赟3

（1．中海油研究總院，北京100028；2．中國(guó)科學(xué)院地質(zhì)與地球物理研究所，北京100029；3．中國(guó)地質(zhì)大學(xué)地球物理與信息技術(shù)學(xué)院，北京100083）

海上多分量地震勘探能夠提供豐富的海底地層信息，有效降低常規(guī)縱波勘探的多解性，在氣云成像和流體檢測(cè)等方面具有優(yōu)勢(shì)。由于海上多分量資料一般由四分量海底電纜（4C－OBC）采集，可以同時(shí)記錄上行和下行縱、橫波，因此波場(chǎng)分離處理十分關(guān)鍵。首先介紹了SCHALKWIJK和EDME兩種常用的波場(chǎng)分離計(jì)算公式，然后采用海底水平多層介質(zhì)模型對(duì)其進(jìn)行了測(cè)試，并與在波場(chǎng)延拓時(shí)采用亥姆霍茲算子和坡印廷矢量分離得到的上行和下行縱、橫波理論記錄進(jìn)行了對(duì)比；最后，針對(duì)兩種常用波場(chǎng)分離計(jì)算公式無(wú)法剝離上行微屈反射等干擾波的問(wèn)題，根據(jù)海底界面波型轉(zhuǎn)換關(guān)系推導(dǎo)得到上行聲波有效反射波和PP型縱波的計(jì)算公式，可以有效剝離水檢上行分量的微屈反射和上行縱波分量的轉(zhuǎn)換縱波。

波場(chǎng)分離；雙檢合成；4C－OBC；旋轉(zhuǎn)交錯(cuò)網(wǎng)格；坡印廷矢量

多分量地震勘探能夠同時(shí)獲得縱波與轉(zhuǎn)換橫波資料，既包含了地下構(gòu)造信息，也包含了反映巖性、流體特征以及裂縫信息，因而可以降低常規(guī)縱波勘探的多解性［1－3］。20世紀(jì)末，海上轉(zhuǎn)換波地震技術(shù)使多分量地震勘探進(jìn)入發(fā)展高潮，從理論研究階段開(kāi)始轉(zhuǎn)入工業(yè)化應(yīng)用階段［4－6］。國(guó)內(nèi)在鶯歌海盆地進(jìn)行了海上多分量地震勘探的工業(yè)化應(yīng)用，在淺層氣藏構(gòu)造成像方面取得了良好效果［7－8］。近年來(lái)，在陸上多分量地震勘探迅速發(fā)展的同時(shí)［9－12］，海上多分量地震勘探在采集和處理等方面也取得很大進(jìn)展［13－14］。與陸上相比，海底環(huán)境噪聲低，而且降速帶比較簡(jiǎn)單，噪聲和靜校正影響較小，有利于轉(zhuǎn)換波的記錄及后期處理。但另一方面，海上多分量地震勘探一般采用海底電纜四分量檢波器（4C－OBC）記錄數(shù)據(jù)，包括一個(gè)水檢H分量、一個(gè)陸檢垂直分量Z和兩個(gè)相互正交的水平分量X，Y［4］，不僅能夠記錄到上行縱波和橫波，也能記錄到水層多次波等下行波，波場(chǎng)分離和多次波去除遠(yuǎn)比陸上復(fù)雜。因此，波場(chǎng)分離是海上多分量資料處理的關(guān)鍵問(wèn)題。

常用波場(chǎng)分離方法包括平面波分解［15］、Radon變換［16］、偏振濾波［17］、波動(dòng)方程法［18］等。與其它方法相比，平面波分解法基于入射波在海底界面發(fā)生反射與透射時(shí)振幅響應(yīng)與入射角的關(guān)系，不僅可以分離縱波和橫波，也可以分離上行波和下行波，因而成為海底多分量地震波場(chǎng)分離的主要方法。由于水檢記錄海水中的壓力，可以視為放置在海底之上，而陸檢三分量記錄海底介質(zhì)的質(zhì)點(diǎn)偏振速度，可以視為放置在海底之下［15］，因此海底多分量波場(chǎng)分離包括海底之上的聲波分解和海底之下的彈性波分解兩部分。

海底之上的聲波分解也就是雙檢合成，可以得到水檢記錄的上行和下行聲波分量，從而去除檢波點(diǎn)端水層多次波。在地震波垂直入射的假設(shè)下，BARR［19］發(fā)展了x－t域雙檢合成的方法；BALL等［20］、BALE［21］、SOUBARAS［22］、OSEN等［23］在其基礎(chǔ)上進(jìn)行了改進(jìn)。海底之下的彈性波分解即分離上行縱波和上行橫波?；谕干浜头瓷鋾r(shí)的振幅響應(yīng)以及縱橫波的偏振特性，AMUNDSEN等［24］、DONATI等［25］、SOUBARAS［26］和WANG等［27］發(fā)展了多種分離方法，這些方法求取的物理量是縱波勢(shì)和橫波勢(shì)，一般在τ－p域或ω－k域進(jìn)行。近年來(lái)，SCHALKWIJK等［28］和EDME等［29］分別提出了不同的波場(chǎng)分離公式，既包含海底之上的聲波分解，也包含海底之下的彈性波分解，代表了海底波場(chǎng)分離的最新進(jìn)展。為了發(fā)展更完善的波場(chǎng)分離方法，我們通過(guò)正演獲得波場(chǎng)分離的理論解，測(cè)試兩種方法的分離結(jié)果與理論解的差異，并針對(duì)其存在的問(wèn)題進(jìn)行改進(jìn)。

本文首先介紹了SCHALKWIJK和EDME提出的兩種波場(chǎng)分離公式，然后采用旋轉(zhuǎn)交錯(cuò)網(wǎng)格和混合邊界條件模擬得到一個(gè)海底水平多層模型的多分量記錄，同時(shí)在波場(chǎng)延拓的每一步采用亥姆霍茲算子和坡印廷矢量分離的方法得到上行和下行縱、橫波理論記錄，并將兩種方法的分離結(jié)果與理論解進(jìn)行對(duì)比；最后，通過(guò)剝離水檢上行分量的微屈反射和上行縱波分量的轉(zhuǎn)換縱波，對(duì)現(xiàn)有方法進(jìn)行了補(bǔ)充完善。

1 兩種海底多分量波場(chǎng)分離公式

如果不考慮海底界面上的反射、透射和波型轉(zhuǎn)換，海底電纜四分量檢波器能夠記錄的入射波包括水中下行聲波、地下反射界面上行縱波和橫波。如果考慮到水檢可以視為放置在海底之上，陸檢三分量視為放置在海底之下，那么上述3種入射波在海底之上和海底之下存在的所有透射、反射和波型轉(zhuǎn)換關(guān)系如圖1所示。其中，“＋”表示下行波，“－”表示上行波，“R”表示反射，“T”表示透射。由圖1可以看到，位于海底之上的水檢不僅可以記錄到由水層多次波和檢波點(diǎn)端多次波組成的下行聲波Pw＋，也可以記錄到由上行縱波透射形成的縱波、上行橫波透射形成的轉(zhuǎn)換縱波以及下行聲波反射形成的微屈反射Pw－peg組成的上行聲波Pw－，前兩者包含地下反射界面的有效信息，可以記為上行聲波有效分量Pw－pri。位于海底之下的陸檢三分量不僅可以記錄到上行縱波P－和上行橫波S－，也可以記錄到由上行縱波在海底界面反射形成的縱波、上行橫波反射形成的轉(zhuǎn)換縱波與下行聲波透射形成的縱波組成的下行縱波P＋，以及由上行縱波反射形成的轉(zhuǎn)換橫波、上行橫波反射形成的橫波與下行聲波透射形成的轉(zhuǎn)換橫波組成的下行橫波S＋。綜上所述，經(jīng)過(guò)海底界面上的反射、透射和波型轉(zhuǎn)換，海底四分量檢波器可以記錄上行和下行的縱／橫波，上行波和下行波的組成十分復(fù)雜，但其中的有效信號(hào)只有上行縱波和橫波。

根據(jù)上行和下行縱／橫波在海底四分量檢波器上的振幅響應(yīng)關(guān)系，SCHALKWIJK等［28］提出了一種在ω－p域進(jìn)行波場(chǎng)分離的公式，即在海底之上進(jìn)行聲波分解：

圖1 海底之上和海底之下的透射、反射和波型轉(zhuǎn)換關(guān)系

同時(shí)在海底之下進(jìn)行彈性波分解的公式為：

式中：G＝4p2v2s，fqp，fqs，f＋v2s，f（q2s，f－p2）2；Φ是縱波勢(shì)；Ψ是橫波勢(shì)；H為水檢分量；Vx為陸檢水平分量；Vz為陸檢垂直分量；τzz為檢波點(diǎn)附近海底介質(zhì)的正應(yīng)力；τxz為剪應(yīng)力；p為水平慢度；qw為水中聲波的垂直慢度；qp，f，qs，f是海底介質(zhì)中縱波和橫波的垂直慢度；ρw為海水密度；vp，f，vs，f，ρf分別表示海底檢波點(diǎn)處的縱、橫波速度和密度。

根據(jù)上行縱波、上行橫波和下行聲波在海底界面上的透射、反射和波型轉(zhuǎn)換關(guān)系，EDME等［29］提出了在τ－p域進(jìn)行波場(chǎng)分離的公式：

公式（4）上面兩行表示海底之上的聲波分解，下面兩行表示海底之下的彈性波分解。

可以看到，兩種波場(chǎng)分離公式的主要區(qū)別為：SCHALKWIJK公式既可以求取陸檢記錄的上行縱波和橫波，也可以求取下行縱波和橫波；EDME公式考慮到陸檢分量的所有下行縱波和橫波都是由入射波反射或透射形成，公式推導(dǎo)時(shí)將其振幅響應(yīng)歸溯到入射波分量，因而沒(méi)有計(jì)算下行縱波和橫波。

2 模型試算

以海底水平多層各向同性介質(zhì)模型為例，模擬得到水檢H分量、陸檢X分量和Z分量的地震記錄，計(jì)算上行和下行縱／橫波記錄的理論解并作為波場(chǎng)分離測(cè)試的基礎(chǔ)。模型尺寸為1 600m×600m，如圖2所示，第1層為海水，水深120m，其它層為各向同性介質(zhì)，具體物理參數(shù)如表1所示。如此設(shè)置模型參數(shù)，是為了分離直達(dá)波、反射波、轉(zhuǎn)換波以及多次波，以便在炮記錄上識(shí)別。在x＝800m，z＝8m網(wǎng)格點(diǎn)上震源激發(fā)純縱波，并將檢波點(diǎn)等間隔布置在海底界面上，x和z方向的空間采樣間隔為2m，時(shí)間采樣間隔為0．5ms，記錄時(shí)長(zhǎng)為1．0s。采用時(shí)間二階、空間十階差分精度旋轉(zhuǎn)交錯(cuò)網(wǎng)格［30］進(jìn)行波場(chǎng)延拓計(jì)算，在左、右、下邊界構(gòu)造完全匹配層進(jìn)行吸收衰減，同時(shí)采用應(yīng)力鏡像法在上邊界構(gòu)造自由界面［31－32］，最終得到中間放炮多分量地震記錄，如圖3所示。

圖2 海底水平多層介質(zhì)模型

表1 海底水平多層模型參數(shù)

為便于分析，對(duì)圖3中的直達(dá)波、折射波、反射縱波、轉(zhuǎn)換波以及水層多次波進(jìn)行了標(biāo)注，其中，D代表直達(dá)波，R代表折射波，P代表反射縱波，C代表轉(zhuǎn)換橫波，M1和M2分別代表1階和2階水層多次波，M3和M4分別代表炮點(diǎn)端和檢波點(diǎn)端多次波?？梢钥吹剑合滦械乃畬佣啻尾ê蜋z波點(diǎn)端多次波以及上行的炮點(diǎn)端多次波在海底四分量記錄上都存在，表現(xiàn)出明顯的周期性；陸檢水平和垂直分量同時(shí)存在上行和下行的縱、橫波，波組相互混雜。

圖3 海底水平多層介質(zhì)模型多分量記錄

為了計(jì)算上行和下行縱／橫波記錄的理論解，在波場(chǎng)正向延拓的每一步，都將矢量波場(chǎng)分離成上行和下行的縱、橫波波場(chǎng)，并由布置在海底的檢波器記錄，最終得到上行和下行的縱、橫波炮記錄。由于縱波是一種散度場(chǎng)，橫波是一種旋度場(chǎng)，根據(jù)亥姆霍茲分解理論，可以應(yīng)用散度算子和旋度算子來(lái)求取彈性波場(chǎng)中的縱波和橫波成分［33］：

圖4 海底水平多層介質(zhì)模型t＝600ms時(shí)的波場(chǎng)分離快照

其中，F(xiàn)up表示上行波，F(xiàn)down表示下行波。圖4e和圖4f分別為采用坡印廷矢量法分解得到的縱波和橫波的上、下行波分量（其中，紅色表示上行波，藍(lán)色表示下行波）?？梢钥吹?，除上、下行波相交時(shí)以一方為主而造成同相軸截?cái)?，波?chǎng)的上、下行波分量都能得到很好分離。

在波場(chǎng)正向延拓的每一步都進(jìn)行上述分離計(jì)算，最后得到波場(chǎng)分離炮記錄如圖5所示。由圖5可以看到：原來(lái)連續(xù)的直達(dá)波和折射波被準(zhǔn)確分離到上行波和下行波炮記錄；上行縱波炮記錄以反射縱波為主，但仍然存在炮點(diǎn)端多次波，下行縱波炮記錄以直達(dá)波、水層多次波和檢波點(diǎn)端多次波為主；與X分量相比，上行橫波記錄沒(méi)有上行縱波和下行多次波干擾，此時(shí)上行橫波記錄中的折射波為轉(zhuǎn)換折射波，在折射縱波之后到達(dá)，與其存在時(shí)差，同理下行橫波記錄中的轉(zhuǎn)換波主要由上行橫波在海底界面反射形成，與上行橫波也存在時(shí)差。

圖5 海底水平多層介質(zhì)模型波場(chǎng)分離炮記錄理論解

圖6 應(yīng)用SCHALKWIJK公式進(jìn)行波場(chǎng)分離的結(jié)果

我們分別采用SCHALKWIJK公式和EDME公式進(jìn)行波場(chǎng)分離試算，結(jié)果如圖6和圖7所示?？梢钥吹?，兩種方法計(jì)算的上行聲波和下行聲波基本一致，下行聲波沒(méi)有反射縱波和折射波，以下行多次波為主，表明聲波分解準(zhǔn)確；兩種方法得到的上行橫波與理論解基本一致，轉(zhuǎn)換波的到達(dá)時(shí)與上行橫波相同，并保留了轉(zhuǎn)換折射波，殘留的下行直達(dá)波主要由τ－p變換對(duì)波組的補(bǔ)償效應(yīng)引起；如果不考慮“剪刀差”現(xiàn)象，兩種方法得到的上行縱波與理論解也基本一致，但是由于分別將連續(xù)的直達(dá)波和折射波分配到上行波和下行波分量，從而造成了波組截?cái)?，所以容易引起“剪刀差”現(xiàn)象［35］，而且SCHALKWIJK公式的“剪刀差”問(wèn)題更為嚴(yán)重，這可能是由公式具體形式引起。

圖7 應(yīng)用EDME公式進(jìn)行波場(chǎng)分離的結(jié)果

3 波場(chǎng)分離改進(jìn)策略

3．1 微屈反射的剝離

SCHALKWIJK公式和EDME公式的聲波分解雖然準(zhǔn)確，但并不徹底，從圖6和圖7可以看到，上行聲波仍然存在較強(qiáng)的多次波，這是直達(dá)波、水層多次波和檢波點(diǎn)端多次波下行到海底形成的上行微屈反射，即圖1中的，它與地下界面的有效反射共同組成了上行聲波Pw－，由于它不能直接反映海底地層信息，因而也是一種干擾波，需要進(jìn)一步剝離。

考慮到上行聲波包括上行微屈反射和有效反射，下行聲波包括直達(dá)波和下行多次波，而所有的上行微屈反射均由下行波在海底反射形成，若以表示直達(dá)波表示下行多次波，Rpp表示海底反射系數(shù)，那么可以得到以下關(guān)系式：

圖8對(duì)比了采用公式（9）求取的上行聲波有效反射與SCHALKWIJK和EDME公式求取的上行聲波。可以看到，由下行多次波引起的微屈反射M1和M2基本被剝離，表明公式（9）準(zhǔn)確；上行聲波以有效反射為主，還包括地層間多次波P2和炮點(diǎn)端上行多次波M3，此時(shí)在一階水層多次波M1的位置殘留的波組是同時(shí)到達(dá)的地層間多次波P2，由于剝離了一階水層多次波，因而振幅大幅減弱。

3．2 SP型上行縱波的剝離

目前的波場(chǎng)分離方法存在的另一個(gè)問(wèn)題是分離后的上行縱波中仍然存在SP型縱波，如圖9所示，它與全程都以縱波形式傳播的PP型縱波共同組成了上行縱波，但是其射線(xiàn)路徑不對(duì)稱(chēng)，從而影響了速度分析和偏移。

由于SP型縱波與PS型橫波一般同時(shí)到達(dá)，上行橫波包括PS型和SS型，而上行聲波有效反射由PP，SP，PS和SS入射波透射形成，即：

因此，在準(zhǔn)確求取上行橫波的前提下，我們提出一種剝離SP型縱波從而得到純PP型縱波的策略，假設(shè)上行聲波有效反射中由SP，PS以及SS入射波形成的波組與同時(shí)到達(dá)的上行橫波在振幅上存在一定的比值，那么可以由（11）式計(jì)算上行純縱波：

其中，r表示PS入射波所在時(shí)窗內(nèi)上行聲波有效反射與上行橫波的振幅比值。圖9對(duì)比了采用公式（11）計(jì)算得到的上行純縱波與EDME公式求取的上行縱波?？梢钥吹?，部分SP型縱波被剝離，而且不會(huì)出現(xiàn)“剪刀差”問(wèn)題，表明求取上行純縱波的策略有效。

圖8 應(yīng)用SCHALKWIJK和EDME公式求取的上行聲波分量（a）和公式（9）求取的上行有效反射（b）

圖9 應(yīng)用EDME公式求取的上行縱波（a）和公式（11）求取的上行純縱波（b）

4 結(jié)束語(yǔ)

由于海上降速帶比較簡(jiǎn)單，可以采集到比陸上信噪比更高的多分量地震資料，因此，海上多分量地震勘探的發(fā)展和應(yīng)用前景更為廣闊。但是，多次波去除和波場(chǎng)分離是海上多分量地震勘探的關(guān)鍵。本文采用模型正演的方式，得到了上行和下行縱／橫波炮記錄的理論解，并與兩種常用波場(chǎng)分離公式的處理效果進(jìn)行了對(duì)比，具體分析了分離前、后不同波組的變化，并針對(duì)存在的問(wèn)題進(jìn)行了改進(jìn)。獲得以下3點(diǎn)結(jié)論：

1）SCHALKWIJK公式和EDME公式都可以準(zhǔn)確分離上行聲波、下行聲波、上行縱波和橫波，但是由于將連續(xù)的直達(dá)波和折射波分離成下行和上行分量，造成了波組截?cái)啵蚨讦樱璸域計(jì)算得到的上行縱波存在“剪刀差”現(xiàn)象，且SCHALKWIJK公式的“剪刀差”問(wèn)題更為嚴(yán)重。

2）原有波場(chǎng)分離公式計(jì)算的上行聲波并不徹底，在下行多次波的位置仍然存在微屈反射，根據(jù)微屈反射均由下行波在海底反射形成的原理，推導(dǎo)了從上行聲波中剝離微屈反射的公式，從而得到了更純凈的上行有效反射。

3）提出一種由上行聲波有效反射和上行橫波計(jì)算上行純縱波的方法，可以剝離部分SP型縱波，并避免了“剪刀差”問(wèn)題。

下一步研究方向?yàn)閷⒏倪M(jìn)后的波場(chǎng)分離方法用于實(shí)際資料處理。

［1］ 張永剛，王赟，王妙月．目前多分量地震勘探中的幾個(gè)關(guān)鍵問(wèn)題［J］．地球物理學(xué)報(bào)，2004，47（1）：151－155

ZHANG Y G，WANG Y，WANG M Y．Some key problems in the multi－component seismic exploration［J］．Chinese Journal of Geophysics，2004，47（1）：151－155

［2］ 黃中玉，曲壽利，王于靜，等．三維多分量地震資料處理技術(shù)研究［J］．石油物探，2010，49（2）：140－146

HUANG Z Y，QU S L，WANG Y J，et al．The studies on 3－D multi－component seismic data processing techniques［J］．Geophysical Prospecting for Petroleum，2010，49（2）：140－146

［3］ 趙波，王赟，蘆俊．多分量地震勘探技術(shù)新進(jìn)展及關(guān)鍵問(wèn)題探討［J］．石油地球物理勘探，2012，47（3）：506－519

ZHAO B，WANG Y，LU J．Recent advances of multicomponent seismic and some of its key issues［J］．Oil Geophysical Prospecting，2012，47（3）：506－519

［4］ CALDWELL J．Marine multicomponent seismology［J］．The Leading Edge，1999，18（11）：1274－1282

［5］ STEWART R R，GAISER J E，BROWN R J，et al．Converted－wave seismic exploration：methods［J］．Geophysics，2002，67（5）：1348－1363

［6］ BARKVED O，BARTMAN B，GAISER J，et al．The many facets of multicomponent seismic data［J］．Oilfield Review，2004，16（2）：42－56

［7］ 張樹(shù)林，李緒宣，姜立紅．海上多波多分量地震技術(shù)新進(jìn)展與發(fā)展方向［J］．物探化探計(jì)算技術(shù)，2000，22（2）：97－107

ZHANG S L，LI X X，JIANG L H．Improvement and development of China offshore multiwave and multicomponent seismic technique［J］．Computing Techniques for Geophysical and Geochemical Exploration，2000，22（2）：97－107

［8］ 傅旦丹，劉一峰，沈麗麗，等．南海某構(gòu)造二維多分量地震資料處理［J］．石油物探，2004，43（2）：181－184

FU D D，LIU Y F，SHEN L L，et al．Processing the 2D multicomponent seismic data of a structure in South China Sea［J］．Geophysical Prospecting for Petroleum，2004，43（2）：181－184

［9］ LI X Y，ZHANG Y G．Seismic reservoir characterization：how can multicomponent data help？［J］．Journal of Geophysics and Engineering，2011，8（2）：123－141

［10］ 程冰潔，唐建明，徐天吉．轉(zhuǎn)換波三維三分量地震勘探技術(shù)的研究現(xiàn)狀及發(fā)展趨勢(shì)［J］．石油天然氣學(xué)報(bào)，2008，30（2）：235－238

CHENG B J，TANG J M，XU T J．Status and development of converted wave 3D－3Cseismic technologies［J］．Journal of Oil and Gas Technology，2008，30（2）：235－238

［11］ 陳海峰，李向陽(yáng)，趙波，等．三維轉(zhuǎn)換波處理技術(shù)研究及應(yīng)用［J］．石油物探，2013，52（5）：512－518

CHEN H F，LI X Y，ZHAO B，et al．Study on 3Dconverted wave data processing method and application［J］．Geophysical Prospecting for Petroleum，2013，52（5）：512－518

［12］ DONATI M，PIAZZA J L，ROLLET A，et al．3D－3C multicomponent seismic－a successful fracture characterization case study in Algeria［J］．First Break，2016，34（1）：35－47

［13］ PACAL E E．Seismic imaging with ocean－bottom nodes（OBN）：new acquisition designs and the ATLANTIS 4COBN survey［D］．Houston：University of Houston，2012

［14］ AMUNDSEN L，ROBERTSSON J O．Wave equation processing using finite－difference propagators，part 1：wavefield dissection and imaging of marine multicomponent seismic data［J］．Geophysics，2014，79（6）：T287－T300

［15］ MUIJS R，ROBERTSSON J O，HOLLIGER K．Datadriven adaptive decomposition of multi－component seabed seismic recordings［J］．Geophysics，2004，69（5）：1329－1337

［16］ JIANG X X，LIN J，YE F，et al．Separation of P－P and P－SV wavefields by high resolution parabolic Radon transform［J］．Journal of Applied Geophysics，2015，119：192－201

［17］ LU J，WANG Y，YAO C．Separating P－and S－waves in an affine coordinate system［J］．Journal of Geophysics and Engineering，2012，9（1）：12－18

［18］ SUN R，MCMECHAN G A，HSIAO H H，et al．Separating P－and S－waves in prestack 3Delastic seismograms using divergence and curl［J］．Geophysics，2004，69（1）：286－297

［19］ BARR F J．Dual sensor OBC technology［J］．The Leading Edge，1997，16（1）：45－51

［20］ BALL V，CORRIGAN D．Dual sensor summation of noisy ocean－bottom data［J］．Expanded Abstracts of 66thAnnual Internat SEG Mtg，1996：28－31

［21］ BALE R．Plane wave deghosting hydrophone and geophone data［J］．Expanded Abstracts of 68thAnnual Internat SEG Mtg，1998：730－733

［22］ SOUBARAS R．Ocean botton hydrophone and geophone processing［J］．Expanded Abstracts of 66thAnnual Internat SEG Mtg，1996：24－27

［23］ OSEN A，AMUNDSEN L，REITAN A．Removal of water－layer multiples from multi－component sea－bottom data［J］．Geophysics，1999，64（3）：838－851

［24］ AMUNDSEN L，REITAN A．Decomposition of multicomponent sea－floor data into upgoing and downgoing P－and S－waves［J］．Geophysics，1995，60（2）：563－572

［25］ DONATI M S，STEWART R R．P and S wave separation at a liquid－solid interface［J］．Journal of Seismic Exploration，1996，5（2）：113－127

［26］ SOUBARAS R．Multiple attenuation and P－S decomposition of multi－component ocean－bottom data［J］．Expanded Abstracts of 68thAnnual Internat SEG Mtg， 1998：1336－1339

［27］ WANG Y，SINGH S C，BARTON P J．Separation of P and SV wavefields from multi－component seismic data in theτ－pdomain［J］．Applied Geophysics，2006，3（3）：163－168

［28］ SCHALKWIJK K M，WAPENAAR C P A，VERSCHUUR D J．Adaptive decomposition of multicomponent ocean－bottom seismic data into downgoing and upgoing P－and S－waves［J］．Geophysics，2003，68（3）：1091－1102

［29］ EDME P，SINGH S C．Receiver function decomposition of OBC data：theory［J］．Geophysical Journal International，2009，177（3）：966－977

［30］ SAENGER E H，GOLD N，SHAPIRO S A．Modeling the propagation of elastic waves using a modified finitedifference grid［J］．Wave Motion，2000，31（1）：77－92

［31］ 杜啟振，孫瑞艷，張強(qiáng)．組合邊界條件下二維三分量TTI介質(zhì)波場(chǎng)數(shù)值模擬［J］．石油地球物理勘探，2011，46（2）：187－195

DU Q Z，SUN R Y，ZHANG Q．Numerical simulation of three－component elastic wavefield in 2DTTI media in the condition of the combined boundary［J］．Oil Geophysical Prospecting，2011，46（2）：187－195

［32］ JASTRAM C，TESSMER E．Elastic modelling on a grid with vertically varying spacing［J］．Geophysical Prospecting，1994，42（4）：357－370

［33］ DELLINGER J，ETGEN J．Wave－field separation in two－dimensional anisotropic media［J］．Geophysics，1990，55（7）：914－919

［34］ DU Q Z，ZHU Y T，BA J．Polarity reversal correction for elastic reverse time migration［J］．Geophysics，2012，77（2）：S31－S41

［35］ 鞏向博．高精度Radon變換及其應(yīng)用研究［D］．吉林：吉林大學(xué)，2008

GONG X B．High precision Radon transform and its application［D］．Jilin：Jilin University，2008

（編輯：陳 杰）

The model testing and the improvement strategy for two decomposition methods of seabed multi－component seismic recordings

ZHANG Jianli1，2，LIU Zhibin1，ZHANG Yunpeng1，LI Weixin1，WANG Yun3
（1．CNOOC Research Institute，Beijing100028，China；2．Institute of Geology and Geophysics，Chinese Academy of Sciences，Beijing 100029，China；3．School of Geophysics and Information Technology，China University of Geosciences，Beijing100083，China）

Many necessary messages about gas reservoir can be obtained from ocean multi－component seismic prospecting，which effectively reduces the multi－solutions of conventional P－wave prospecting and has the advantages of gas－cloud imaging and fluid detecting．As ocean multi－component data is usually recorded by 4Csensors on the ocean bottom cable，abbreviated as 4C－OBC，the down－going P－wave and S－wave will be

as well as up－going P－wave and S－wave．Therefore，wavefield decomposition is the key problem for ocean multi－component seismic exploration．In this paper，SCHALKWIJK and EDME the two common calculation formulae of wavefield decomposition are introduced firstly．Secondly，both methods are tested on a ocean seabed horizontal layered medium model，and the decomposition results of their testing recordings are compared with the theory solutions，which are calculated by applying Helmholtz operator and Poynting vector at every wavefield extrapolation．Finally，the two common wavefield decomposition calculation formulae cant separate pag－leg reflection of up－going interference wave．The calculation formulae of reflections from up－going PP wave of P－wave is derived based on the waveform，transformation from marine interface not only by separating peg－leg reflection from up－going hydrophone component，but also by separating S－P wave from up－going P－wave component．

wavefield decomposition，dual－sensor seismic data merging，4C－OBC，rotated staggered－grid，Poynting vector

P631

A

1000－1441（2017）03－0373－09

10．3969／j．issn．1000－1441．2017．03．007

張建利，劉志斌，張?jiān)迄i，等．兩種海底多分量波場(chǎng)分離方法的模型測(cè)試及改進(jìn)［J］．石油物探，2017，56（3）：373－381

ZHANG Jianli，LIU Zhibin，ZHANG Yunpeng，et al．The model testing and the improvement strategy for two decomposition methods of seabed multi－component seismic recordings［J］．Geophysical Prospecting for Petroleum，2017，56（3）：373－381

2016－03－28；改回日期：2016－05－09。

張建利（1985—），男，博士，主要從事多分量地震勘探和各向異性研究。

國(guó)家杰出青年科學(xué)基金項(xiàng)目（41425017）資助。

This research is financially supported by the National Science Fund for Distinguished Youth Scholars（Grant No．41425017）．