專利名稱：處理地震數(shù)據(jù)的方法和設(shè)備的制作方法
技術(shù)領(lǐng)域：
本發(fā)明涉及處理地震數(shù)據(jù)，尤其是多分量地震數(shù)據(jù)，以從采集到的數(shù)據(jù)中消除不需要的波至(事件，events)的方法。本申請還涉及處理地震數(shù)據(jù)的設(shè)備。本申請還涉及計算去震源特征(de-signature，去小波)和去多次反射(de-multiple)算子，或者計算去震源特征算子的方法和設(shè)備。
背景技術(shù)：
圖1(a)是地震勘探的原理示意圖。地震勘探目的在于提供關(guān)于位于地球內(nèi)部的目標(biāo)地質(zhì)反射體3的信息。
示于圖1的地震勘探包括設(shè)置在地球表面1的地震源4。地震傳感器5(下面稱為接收器)也設(shè)置在地面，離地震源4一定距離。在使用時，使地震源激發(fā)地震能量脈沖，由接收器5檢測發(fā)射出的地震能量。
圖1(a)圖解了陸地地震勘探。但是，這種通用形式的地震勘探不限于陸地，而可以在海洋環(huán)境或者陸海過渡區(qū)進(jìn)行。例如，這樣的海洋地震勘探系統(tǒng)是已知的勘探船拖曳一個或者多個地震源。在這樣的方案中，接收器可以設(shè)置在海床上(所謂的“海底測纜(ocean bottomcable)”勘探或者OBC勘探)，或者接收器也可以由勘探船拖曳著。另外，陸地地震勘探不限于圖1(a)所示的方案，地震源或者地震接收器也可以設(shè)置在地球內(nèi)。例如，在垂直地震剖面(VSP，vertical seismicprofile)地震勘探中，將地震源設(shè)置在地面，接收器設(shè)置在地球內(nèi)部的井孔中。在逆VSP(reverse VSP)勘探中，地震源設(shè)置在井孔中，接收器設(shè)置在地面。
為了方便說明，圖1(a)中僅圖示了一個地震源4和一個地震接收器5。但是，一般而言，實際的地震勘探將包括一個震源陣列和一個接收器陣列。
在地震勘探中，一個問題是地震能量會沿著多條路徑從震源傳播到接收器。這其中的一個原因是除了目標(biāo)反射體之外，在地球中還存在許多其它的反射體。在圖1(a)中，這示意性地以在目標(biāo)反射體3上方的反射體2來表示。這些額外的反射體產(chǎn)生從震源到接收器的地震能量路徑，涉及在目標(biāo)反射體之外的反射體上的反射。存在地震能量的多條路徑的另一個原因是在地球內(nèi)向上傳播的地震能量會在地面1被反射，這是由于地和空氣的不同地震特性。這導(dǎo)致從震源到接收器存在多條地震能量路徑，這涉及在目標(biāo)反射體處的不止一次的反射。這些路徑在接收器所采集到的地震數(shù)據(jù)中導(dǎo)致了不需要的波至。采集到的地震數(shù)據(jù)中與經(jīng)過多次反射的地震能量有關(guān)的波至在下面稱為“多次反射波至”(multiple event)。
圖1(a)圖示了該勘探方案中地震能量的主要路徑，其中，地震能量從震源4到接收器5的路徑只包括在目標(biāo)反射體3的一次反射(為了簡明起見，在圖1(a)中省略了在上覆發(fā)射體2處發(fā)生的折射)。盡管沿著該主要路徑傳播的地震能量在其從震源4到目標(biāo)反射體3的向下的路徑上穿過上覆反射體2，并在從目標(biāo)反射體3到接收器5的向上的路徑上再次穿過該上覆反射體2，該主要路徑并不包括在上覆反射體2上的反射。在理想地震勘探中，接收器5只會檢測沿著主要路徑傳播的地震能量。
在實際的地震勘探中，接收器5會檢測到從震源4沿著除主路徑之外的許多路徑傳播的能量。這些其它地震能量路徑的例子示于圖1(b)到1(d)中。圖1(b)到1(d)圖解了涉及一次以上反射的地震能量路徑，這就是所謂的多次反射波至(multiple events)。在圖1(b)中，從震源4向下傳播的地震能量被上覆反射體2反射，從而使其向上傳播到地球表面1。該地震能量又在地球表面被向下反射，然后入射到目標(biāo)反射體上。這種通常類型的波至就是眾所周知的“震源側(cè)多次反射波至”(“source-leg multiple”)，因為在地震能量從震源到目標(biāo)反射體的路徑上發(fā)生了額外的反射。
圖1(c)圖解了地震能量直接從震源傳播到目標(biāo)反射體3，然后在目標(biāo)反射體3向上反射的地震能量路徑。但是，被反射的地震能量不是直接入射到接收器，而是在地球表面向下反射，然后在到達(dá)接收器之前在上覆反射體2處向上反射。這種類型的地震路徑公知為“接收器側(cè)多次反射波至”(“receiver-leg multiple”)，因為在地震能量從目標(biāo)反射體到接收器的路徑上發(fā)生了額外的反射。
圖1(d)圖解了一種地震路徑，其中，來自震源的地震能量入射到目標(biāo)反射體上，向上反射到地球表面，然后向下反射，在目標(biāo)反射體3上又一次反射，而后到達(dá)接收器。在這樣的地震能量路徑中，在地震能量從震源到目標(biāo)反射體的路徑以及地震能量從目標(biāo)反射體到接收器的路徑之間產(chǎn)生了額外的路徑。
在實際的地震勘探中，在接收器處采集到的地震能量包括與所需的主要路徑1(a)相應(yīng)的波至，但是也包括與不想要的多次反射路徑(multiple paths)比如示于圖1(b)到1(d)的路徑相關(guān)的波至。為了提供關(guān)于目標(biāo)反射體的精確信息，希望能夠從接收器采集到的地震能量中識別和去除多次反射波至(multiple events)。
在陸地勘探中使用的地震源通常是振動或者爆炸震源。如果使用振動器，則可以進(jìn)行多分量勘探(multi-component survey)使用產(chǎn)生三個正交震源運動(兩個在正交的水平方向，一個在垂直方向)的多分量振動器。如果使用的地震接收器能夠記錄三個正交方向的質(zhì)點運動，則可以進(jìn)行3C×3C(或者9C)地震勘探。用于此目的的合適的接收器是能夠測量接收器處的質(zhì)點運動的三個正交分量的接收器，例如包含三個正交地震檢波器的接收器兩個地震檢波器用于測量接收器處質(zhì)點運動的兩個正交水平分量，第三個地震檢波器用于測量接收器處質(zhì)點運動的垂直分量。
地震勘探結(jié)果的分析的另一個問題是多分量振動器(其通過對地面施加推力而工作)發(fā)射三種不同類型的波，也就是P波、Sv波和Sh波(P波是壓力波，Sv和Sh波是剪切波)。多分量振動器發(fā)射的地震能量中，這些不同類型的波的相對幅度隨著地震能量傳播方向而變化。一個多分量接收器記錄三種波類型，其靈敏度取決于接收到的地震能量的入射角。當(dāng)?shù)卣饳z波器測量地球表面的波場的一個分量時，P波和兩種類型的S波都被不加區(qū)別地予以記錄。
這示意性地圖解于圖4(a)中。圖4(a)圖示了一種陸地地震勘探，其中，在地球表面布置一個三分量地震源4(這里是一個多分量振動器)和一個三分量地震接收器5。如圖4(a)所示，地震源4發(fā)射P波和S波(為了簡明起見，僅圖示了一種S波)，接收器5檢測P波和S波。圖4(a)從左至右圖示了四個圖面(i)在振動器4處由水平震源運動產(chǎn)生，并在接收器5處由水平取向的地震檢波器接收的地震能量；(ii)在振動器4處由水平震源運動產(chǎn)生，在接收器5處由垂直取向的地震檢波器檢測到的地震能量；(iii)在振動器4處由垂直震源運動產(chǎn)生，在接收器5處由水平取向的地震檢波器檢測到的地震能量；(iv)在振動器4處由垂直震源運動產(chǎn)生，在接收器5處由垂直取向的地震檢波器檢測到的地震能量。
圖4(a)僅圖解了振動器4處震源運動的一個水平分量和一個垂直分量，以及接收器5處的一個水平地震檢波器分量和一個垂直地震檢波器分量。如前所述，在全多分量勘探中，振動器4還會由在紙面所在平面以外的震源運動產(chǎn)生地震能量，接收器還會包括檢測沿著紙面外的直線的質(zhì)點運動的第三地震檢波器。這樣，在震源處產(chǎn)生的三個正交震源運動和接收器處的三個正交地震檢波器總共產(chǎn)生9個組合。
如前所述，許多地震接收器無區(qū)別地記錄P波和S波，從而，在接收器處采集到的地震記錄跡包括歸因于接收到的P波的波至和歸因于接收到的S波的波至。在許多情況下，希望將地震記錄跡中的P波至與S波至相分離，因為這可以提供關(guān)于地球內(nèi)部的額外的信息。在許多情況下，地質(zhì)結(jié)構(gòu)會對P波和S波有不同的影響。將地震記錄跡中的P波至與S波至相分離的處理通常稱為將地震記錄跡分解為其P分量和S分量。
有許多現(xiàn)有技術(shù)的方法用于將在接收器處采集的地震能量分解為P分量和S分量。同時也有許多現(xiàn)有技術(shù)的方法用于從采集到的地震數(shù)據(jù)中消除多次反射的影響。
對于陸地地震數(shù)據(jù)，C.P.A.Wapenaar et al.在”Decomposition ofmulticomponent seismic data into primary P-and S-wave responses”，Geophys.Prosp.Vol.38，pp633-661(1990)中，以及P.Herrmann在”Decomposition of multicomponent measurements into P-andS-waves”，Ph.D.thesis，Delft University of Technology(1992)中，給出了一種彈性分解方案分解在接收器側(cè)(共炮點道集)和在震源側(cè)(共接收器道集)兩側(cè)的數(shù)據(jù)。應(yīng)用于沿著接收器線(receiver line)的共炮點道集的波場分解用純P波和S波檢測子(detector)取代原來的質(zhì)點速度檢測子。類似地，應(yīng)用于共接收器道集的波場分解用純P波和S波震源取代原來的振動器震源。這樣，總的分解方案提供9個不同的數(shù)據(jù)集來自模擬P、Sv和Sh震源的經(jīng)過P、Sv和Sh處理的道集。
在進(jìn)行波場分解之后，Wapenaar et al.(1990，supra)應(yīng)用一個反演機制(reverse scheme)，從被分解的結(jié)果中消除地球表面的響應(yīng)。在這種技術(shù)中，將地球表面假設(shè)為固體和真空之間的界面。消除地球表面的響應(yīng)就是從被分解的數(shù)據(jù)中消除涉及地球表面反射的所有多次反射路徑，從而給出地球內(nèi)部的一次P波和S波響應(yīng)。
這種現(xiàn)有技術(shù)要求地震源是點震源，并具有已知的子波(wavelet)，或者能夠從采集到的數(shù)據(jù)計算(估計)子波。
E.Holvik和L. Amundsen在”decomposition of multicomponentsea floor data into primary PP，PS，SP and SS waveresponse”，Expanded Abstracts of 68thAnnual Int Mtg of Society ofExploration Geophysicists，pp2040-2043(1998)中將上述現(xiàn)有技術(shù)推廣到了海洋地震勘探的情況其中，在海底布置理想振動器(推力震源，traction sources)和地震檢波器。這種技術(shù)涉及分解所采集到的數(shù)據(jù)，消除由于水層內(nèi)的反射而產(chǎn)生的波至。
K.Matson和A.Weglein在“Removing of elastic interfacemultiples from land and ocean bottom seismic data using inversescattering”in Expanded Abstracts of 66thAnnual Int Mtg of Societyof Acceleration of Geophysicists，pp1526-1529(1996)中，以及K.Matson在“An inverse scattering series method for attenuatingelastic multiples from multicomponent land and ocean bottom seismicdata”，Ph.D.thesis，University of British Columbia(1997)中，發(fā)展了另一種現(xiàn)有技術(shù)。在這種技術(shù)中，使用逆散射級數(shù)(inverse scatteringseries)來發(fā)展從多分量陸地或者洋底地震數(shù)據(jù)衰減多次反射的彈性方案(elastic scheme)。
由Holvik和Amundsen發(fā)展，并由Matson和Weglein再次發(fā)展的所述級數(shù)要求地震源必須是點震源，并且具有已知的子波，或者可以從采集到的數(shù)據(jù)計算(估計)子波。
L.Amundsen在“Geophysics”Vol 66 pp327-341(2001)中提出了一種從布置在水柱中或者海床上的多分量地震接收器采集到的海洋地震數(shù)據(jù)中消除自由表面多次反射的方法。但是，該方法不能消除與海床相關(guān)聯(lián)的多次反射的影響。

發(fā)明內(nèi)容
本發(fā)明的第一方面在于提供一種處理多分量地震數(shù)據(jù)的方法，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，并在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的，該方法包括下列步驟將采集到的地震數(shù)據(jù)分解為上行成分和下行成分；從采集到的地震數(shù)據(jù)的所述下行成分以及環(huán)繞該接收器的介質(zhì)的特性計算一個去震源特征和去多次反射算子。
本發(fā)明的去震源特征和去多次反射算子能夠有效地衰減或者完全消除地震數(shù)據(jù)中蓋層(overburden)的影響。其能夠有效地從地震數(shù)據(jù)中消除與接收器高度以上的任何界面或者接收器高度的任何界面相關(guān)的所有多次反射(術(shù)語“界面”包括任何導(dǎo)致地震能量部分反射的聲學(xué)或者彈性特性不連續(xù))。
其也能有效地從數(shù)據(jù)中衰減或者完全消除震源輻射特性(或者震源“特征(signature)”)的影響。
該方法可以包括用去震源特征和去多次反射算子處理所采集到的地震數(shù)據(jù)，從而衰減或者消除因為多次反射而產(chǎn)生的地震波至。這可以通過利用去震源特征和去多次反射算子處理所采集到的地震數(shù)據(jù)的上行成分而實現(xiàn)?；蛘撸梢詫⑷フ鹪刺卣骱腿ザ啻畏瓷渌阕討?yīng)用于采集到的完整的(也就是未分解的)地震數(shù)據(jù)。
處理所采集到的地震數(shù)據(jù)的步驟可以包括為震源選定所需的地震源特征(seismic signature)。
該方法還可以包括將地震數(shù)據(jù)分解為P波和/或S波數(shù)據(jù)。該分解步驟可以對去震源特征和去多次反射后的數(shù)據(jù)進(jìn)行?；蛘撸梢栽趹?yīng)用去震源特征和去多次反射算子之前將采集到的數(shù)據(jù)分解為P波和/或S波數(shù)據(jù)。
向P波和/或S波數(shù)據(jù)的分解可以是接收器側(cè)分解和/或震源側(cè)分解。
本發(fā)明的第二方面提供一種處理多分量地震數(shù)據(jù)的方法，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的，該方法包括下列步驟將采集到的地震數(shù)據(jù)分解為上行成分和下行成分；從采集到的地震數(shù)據(jù)的初始下行成分和環(huán)繞接收器的介質(zhì)的特性來計算一個去震源特征算子。
在接收器處采集到的地震數(shù)據(jù)的初始下行成分與直接從震源向接收器傳播而未經(jīng)過任何反射的地震能量有關(guān)。因此它們可以用來計算去震源特征算子，以衰減或者消除地震源的輻射特征對采集到的數(shù)據(jù)的影響。
本發(fā)明的上述方法可以應(yīng)用于已經(jīng)存在的地震數(shù)據(jù)?；蛘撸景l(fā)明的方法可以還包括下述步驟在震源位置發(fā)射多分量地震能量；在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)。
本發(fā)明的其它方面體現(xiàn)在獨立權(quán)利要求10、12、13、14、16、23、24、26、27、28和30中。


下面結(jié)合附圖通過舉例描述本發(fā)明的優(yōu)選實施例。附圖中圖1(a)是地震勘探的示意圖，圖示了希望的地震能量路徑；圖1(b)、1(c)和1(d)是地震勘探的示意圖，圖示了地震能量的三種不希望的路徑；圖2是一個流程圖，圖示了根據(jù)本發(fā)明的一個實施例獲得去震源特征和去多次反射算子的方法的可能步驟；圖3(a)是一個真實地震勘探的示意圖；圖3(b)圖示了對應(yīng)于圖3(a)的地震勘探的理想地震勘探；圖3(c)圖示了又一個理想地震勘探，它是示于圖3(b)的理想地震勘探的逆；圖4(a)到4(d)圖解了將本發(fā)明的方法應(yīng)用于陸地地震勘探的效果；圖5(a)到5(b)圖解了將本發(fā)明的方法應(yīng)用于洋底地震勘探的效果；圖6(a)到4(b)圖解了將本發(fā)明的方法應(yīng)用于陸地地震勘探的效果，其中，接收器布置在地球內(nèi)部；圖7(a)到7(d)圖示了原始綜合地震數(shù)據(jù)；圖8(a)到8(d)圖示了圖7(a)到7(d)的綜合地震數(shù)據(jù)的上行分量；圖9(a)到9(d)圖示了圖7(a)到7(d)的綜合地震數(shù)據(jù)的下行分量；圖10(a)到10(d)圖解了利用本發(fā)明的方法從圖7(a)到7(d)所示的綜合地震數(shù)據(jù)獲得的質(zhì)點速度；圖11(a)到11(d)圖示了利用本發(fā)明的方法從圖7(a)到7(d)所示的綜合地震數(shù)據(jù)獲得的PP、SP、PS和SS波場；圖12為本發(fā)明的設(shè)備的示意框圖。
具體實施例方式
本發(fā)明的獲得去震源特征和去多次反射算子的一種方法的原理示于圖3(a)到3(c)中。圖3(a)圖解了一個實際的地震勘探。該地震勘探包括一個位于深度zr的接收器xr14。在其中進(jìn)行地震勘探的介質(zhì)對于小于接收器深度的所有深度(也就是對于所有的z＜zr)歸類為蓋層，對于所有大于接收器深度的深度(也就是對于所有的z＞zr)歸類為“地下層”(subsurface)。蓋層是任意的非均一介質(zhì)。
一個地震源xs(13)設(shè)置在蓋層中，深度小于接收器深度。該地震源發(fā)射上行和下行地震能量。
用蓋層內(nèi)的界面10、11和12的存在示意性地表示蓋層的非均一性質(zhì)。蓋層的地震特性在這些界面處發(fā)生變化。在海洋地震勘探的情況下上界面12可以是海面(即空氣-水體界面)。由于在蓋層中存在這些界面，上行地震能量被在蓋層中向下反射，從而在接收器14所采集到的地震數(shù)據(jù)中出現(xiàn)多次反射波至。
地震源xs(13)既發(fā)射S波又發(fā)射P波。接收器xr(14)不加區(qū)別地記錄接收到的S波和P波。接收器可以是任何測量質(zhì)點運動比如質(zhì)點位移、質(zhì)點速度(質(zhì)點位移的時間導(dǎo)數(shù))或者質(zhì)點加速度(質(zhì)點速度的時間導(dǎo)數(shù))的設(shè)備。通常使用的地震接收器是地震檢波器(geophone)，但是本發(fā)明不限于處理使用一個或者多個地震檢波器作為接收器獲取的地震數(shù)據(jù)。
圖3(b)圖解了對應(yīng)于圖3(a)的理想地震勘探方案。在圖3(b)的理想化方案中，蓋層是均一介質(zhì)，因此來自地震源的上行地震能量繼續(xù)向上傳播而不向下反射，因此在接收器處采集到的地震能量不包含多次反射波至。在圖3(b)中，圖3(a)的蓋層中的界面10、11和12被圖示為虛線，以強調(diào)圖3(b)中的蓋層是均一的。
在圖3(b)的理想地震勘探方案中，地震源13對中于位置xr-，位于深度zr-。地震源的深度zr-是在接收器14所在深度zr上方的無窮小距離E處。該距離E為趨于零的小值。
圖3(c)圖示了一種替代的理想化地震勘探方案。這對應(yīng)于圖3(b)的理想化地震勘探方案，不同之處在于，地震源13和地震接收器14的位置互換了。也就是，在圖3(c)中，地震源13位于深度為zr的位置xr，地震接收器14位于深度為zr-的位置xr-，該深度zr-在源的深度zr上方距離E處。
由于圖3(c)的地震勘探方案除了源和接收器的位置互換外與圖3(b)的地震勘探方案相同，互易定理(theorem of reciprocity)要求示于圖3(b)的地震勘探方案與示于圖3(c)的地震勘探方案產(chǎn)生相同的地震數(shù)據(jù)。
在本發(fā)明的一個實施例中，使用貝蒂(Betti)互易定理將圖3(a)所示類型的實際地震勘探采集到的實際地震數(shù)據(jù)轉(zhuǎn)換為圖3(b)所示類型的理想化地震勘探所應(yīng)獲得的地震數(shù)據(jù)。這能夠消除非均一的蓋層的影響，從而消除多次反射波至。
本發(fā)明還提供分解算子，將所采集的地震數(shù)據(jù)分解為PP、PS、SP和SS分量。
貝蒂互易定理給出了在數(shù)學(xué)或者物理面所包圍的特定空間中所限定的兩個獨立彈性波場之間的積分方程關(guān)系。兩個波場之間的關(guān)系是由介質(zhì)參數(shù)可能存在的差別、源的分布可能存在的差別以及邊界條件可能存在的差別所決定的。在本發(fā)明中，貝蒂定理用于提供能夠用來將具有蓋層響應(yīng)的實際地震勘探(即由于非均一蓋層而產(chǎn)生了多次反射的勘探)中記錄的波場轉(zhuǎn)換為不存在蓋層響應(yīng)的理想化假想地震試驗中應(yīng)當(dāng)記錄到的波場(假想地震勘探方案中的蓋層是均一的，因此不產(chǎn)生響應(yīng))。本發(fā)明使得能夠?qū)⒌卣饠?shù)據(jù)中由于圖1(b)到1(d)所示類型的多次反射導(dǎo)致的不需要的波至消除，從而能夠獲得無多次反射的圖像。本發(fā)明使得能夠從地震數(shù)據(jù)中消除與接收器高度或其上方的任何界面相關(guān)的任何多次反射。
除了源的位置之外，為消除蓋層散射的波，本發(fā)明不需要任何用于計算的源特性。在從實際地震勘探方案到理想化假想地震勘探方案的轉(zhuǎn)換中，消除了實際地震源的輻射特性。假想理想化地震勘探方案中的源被假設(shè)為單點源，其產(chǎn)生相互正交取向的分量，并具有所需的子波。
從獲取的地震數(shù)據(jù)中消除蓋層的響應(yīng)以及多分量震源輻射特性的影響的方法的該實施例在下面稱為“貝蒂去震源特征和彈性去多次反射”(Betti designature and elastic demultiple)。本發(fā)明的貝蒂去震源特征和彈性去多次反射方法具有以下優(yōu)點(a)在消除地震數(shù)據(jù)中由于上行能量在蓋層中的散射而產(chǎn)生的所有波至的同時，保持了一次反射(primary reflection)的幅度；(b)不需要關(guān)于接收器深度以下的介質(zhì)的知識；(c)不需要關(guān)于接收器深度以上的介質(zhì)的知識；(d)僅需要在接收器深度處的局部密度和局部彈性波傳播速度的信息。
在本發(fā)明的一個方法中，記錄的地震數(shù)據(jù)被轉(zhuǎn)換為新的、理想化的、在沒有蓋層響應(yīng)的假想地震勘探(也就是，其中的蓋層是完全均一的)中應(yīng)當(dāng)記錄到的數(shù)據(jù)。該假想地震勘探中的震源是一個點源，具有某種所需的震源輻射特征(子波)。在經(jīng)過了處理的數(shù)據(jù)中，不存在實際震源及其輻射特性的影響。換句話說，新數(shù)據(jù)已經(jīng)被“去震源特征”了。
本發(fā)明的方法可應(yīng)用于在任何使用多分量地震源和多分量地震接收器的地震勘探中所獲得的數(shù)據(jù)。可應(yīng)用本發(fā)明的三種具體的地震勘探方案如下A.3C×3C(或者9C)陸地地震勘探，其中，地震源產(chǎn)生三個正交源運動(通常是兩個水平源運動，一個垂直源運動)。地震波場由一接收器記錄，該接收器包括三個正交的地震檢波器，布置在地球表面或者地球表面以下。這三個正交地震檢波器測量質(zhì)點速度矢量的三個正交分量(通常是x、y和z分量)。
B.3C×4C(或者12C)海洋地震勘探，其中，地震源產(chǎn)生三個正交源運動(通常是兩個水平源運動，一個垂直源運動)。地震波場由一接收器記錄，該接收器包括三個剛好布置在海床以下的正交的地震檢波器和一個剛好布置在海床以上的地震檢波器。所述三個正交地震檢波器測量質(zhì)點速度矢量的三個正交分量，所述一個地震檢波器測量壓力場(是一個標(biāo)量)。
C.3C×6C(或者18C)井孔地震試驗，其中，地震源同樣產(chǎn)生三個獨立的正交源運動。地震源可以布置在陸地上、海底或者水柱中。地震接收器最好測量質(zhì)點速度矢量的三個正交分量，以及垂直推力矢量(vertical traction vector)的三個正交分量或其估計值。
根據(jù)本發(fā)明，從所獲取的地震數(shù)據(jù)的下行成分獲得去震源特征和去多次反射算子。這意味著在接收器處進(jìn)行的測量必須足以允許將所采集到的波場在接收器側(cè)分解為在接收器處的上行和下行波成分。在最為通常的情況中，合適的波場分解要求接收器測量質(zhì)點速度矢量的三個分量以及垂直推力矢量的三個分量。但是，對于陸地震數(shù)據(jù)，垂直推力矢量為零，因此只需要測量質(zhì)點速度矢量的三個分量。
對于洋底地震數(shù)據(jù)，垂直推力矢量的兩個水平分量為零，剛好在海床以下的垂直推力矢量的垂直分量與剛好在海床以上的壓力場的幅度相等、符號相反。這樣，剛好在海床上方設(shè)置一個水聽器就能夠測定垂直推力矢量的垂直分量(推力矢量的唯一非零分量)。
將本發(fā)明應(yīng)用于井孔地震勘探原則上要求必須測量或者知道垂直推力矢量的三個分量。但是，在多數(shù)情況下，垂直推力矢量的三個分量是不知道的。為了克服此問題，通常是在井孔內(nèi)的兩個不同深度記錄三分量質(zhì)點速度。這樣，接收器系統(tǒng)實際上就是一個六分量接收器系統(tǒng)，因此井孔地震勘探可以視為18C勘探。
圖2的流程示了本發(fā)明的方法的主要步驟。
本發(fā)明可以在合適的現(xiàn)有地震數(shù)據(jù)上進(jìn)行。在這種情況下，該方法始于從存儲裝置檢索合適的地震數(shù)據(jù)的步驟(未示出)?；蛘?，該方法可以開始于執(zhí)行多分量震源、多分量接收器地震勘探以獲取適當(dāng)?shù)牡卣饠?shù)據(jù)的步驟(未示出)。
一開始，在步驟S1，確定接收器附近的介質(zhì)的局部彈性特性。對于海洋地震勘探，作為例子，該步驟要求確定接收器附近的海床的局部彈性特性。這些特性可以從測量來確定，或者可以根據(jù)對勘探地點的地質(zhì)結(jié)構(gòu)的知識來估計。
在步驟S2，將已采集到的或者從存儲裝置檢索出的多分量地震數(shù)據(jù)分解為上行和下行波場成分。一般，步驟S2可以與步驟S1同時或者在它之前進(jìn)行。但是，通常最好在步驟S2之前執(zhí)行步驟S1，因為從步驟S1獲得的局部彈性特性然后可以用于將地震數(shù)據(jù)分解為上行和下行波。
在步驟S3，從所述多分量數(shù)據(jù)和介質(zhì)的所述彈性特性計算去震源特征和去多次反射算子。
在步驟S4，選擇地震源的所希望的震源特征。然后將該所希望的震源特征、在步驟S2所獲得的多分量數(shù)據(jù)以及在S3計算的算子用于解一個積分方程，來在步驟S5找到去震源特征并去多次反射后的波場。本質(zhì)上，步驟S5是將圖3(a)所示類型的實際地震勘探所獲得的地震數(shù)據(jù)轉(zhuǎn)換為圖3(b)所示的理想化地震勘探所應(yīng)獲得的數(shù)據(jù)。在一個優(yōu)選實施例中，利用貝蒂定理將實際測量的、圖3(a)的地震勘探方案的面S所圍的波場轉(zhuǎn)換為在圖3(b)所示的理想化地震勘探方案中由同一個面所圍的波場。該面S由示于圖3(a)到3(c)所示的面∑和SR構(gòu)成。所述面∑是一個位于深度zr-，也就是位于接收器(位于深度zr)深度上方無線小距離處的水平平面。所述面SR是半徑為R的半球面。
在步驟S5，可以將所述去震源特征和去多次反射算子應(yīng)用于在步驟S1獲得的上行波場成分?；蛘撸梢詫⑺鋈フ鹪刺卣骱腿ザ啻畏瓷渌阕討?yīng)用于在接收器處原始獲得的地震數(shù)據(jù)。
最后，在步驟S6，將去震源特征和去多次反射后的數(shù)據(jù)分解為如果在地震勘探中使用純P波和S波震源并使用P波和S波接收器會記錄到的一次PP、PS、SP和SS波響應(yīng)。
在一個實施例中，步驟S6具有兩個主要操作。一開始，在水平接收器坐標(biāo)上對多分量數(shù)據(jù)進(jìn)行傅立葉變換。然后在波數(shù)域?qū)ψ儞Q后的數(shù)據(jù)進(jìn)行接收器側(cè)波場分解，得到上行壓力波和剪切波場。這是通過用一個接收器分解矩陣對數(shù)據(jù)進(jìn)行操作而完成的。然后對分解后的數(shù)據(jù)進(jìn)行傅立葉逆變換，得到在接收器處測量的上行壓力波和剪切波場。
然后在水平源坐標(biāo)上對這樣獲得的多分量數(shù)據(jù)進(jìn)行傅立葉變換。然后在波數(shù)域?qū)ψ儞Q后的數(shù)據(jù)進(jìn)行源側(cè)波場分解，得到下行壓力波和剪切波場。這是通過用一個源分解矩陣對數(shù)據(jù)進(jìn)行操作而完成的。然后對分解后的數(shù)據(jù)進(jìn)行傅立葉逆變換，得到在源處發(fā)射的下行壓力波和剪切波場。
在源側(cè)和接收器側(cè)的彈性波場分解的組合給出相當(dāng)于在理想化勘探中獲得的數(shù)據(jù)的地震數(shù)據(jù)，在理想化勘探中，不存在蓋層響應(yīng)，并且使用單分量P波源或者S波源以及單分量P波接收器或者S波接收器。
本發(fā)明的又一個優(yōu)點是消除或者至少部分衰減地滾波噪聲(在陸地地震勘探的情況下)或者斯科特波(Scholte wave)(在海洋地震勘探的情況下)。在陸地地震勘探中，地滾波噪聲通常是最大的噪聲源，因此非常希望消除地滾波噪聲的影響。
如果需要，可以省略步驟S6。或者，可以在應(yīng)用去震源噪聲和去多次反射算子的步驟S5之前進(jìn)行將數(shù)據(jù)分解為PP、PS、SP和SS波場的步驟S6。如果首先進(jìn)行分解為PP、PS、SP和SS波場的操作，該步驟將取代將地震數(shù)據(jù)分解為上行和下行波的步驟S2。
如果在其中進(jìn)行地震勘探的介質(zhì)水平分層，則可以大大簡化貝蒂去震源特征和去彈性多次反射的方案，以及接下來的彈性震源-接收器分解的方案。這可以在τ-P域或者頻率-波數(shù)域方便地進(jìn)行。
本方法最好在共炮點道集上進(jìn)行，或者，如果震源陣列的變化可忽略，最好在共接收器道集上進(jìn)行。在后一種情況下，地球分層的假定可能降低結(jié)果的質(zhì)量，但在大多數(shù)情況下這不顯著。
在利用本發(fā)明的去震源特征以及去多次反射的方法針對蓋層和震源輻射特性的影響校正了地震數(shù)據(jù)之后，以及可選地，在地震數(shù)據(jù)被分解為P波和/或S波成成分之后，然后可以對數(shù)據(jù)進(jìn)行進(jìn)一步處理步驟，例如一個或者多個傳統(tǒng)的處理步驟。
圖4圖解了本發(fā)明在陸地地震勘探中的效果，在該陸地地震勘探中，一個多分量地震源4和一個多分量地震接收器5設(shè)置在地球表面。
圖4(a)圖解了一個真實的地震勘探。當(dāng)激發(fā)震源4時，發(fā)射下行P波和S波。P波用實線表示，S波用虛線表示。表示源的方塊(標(biāo)為“src”)中的箭頭表示源產(chǎn)生的源運動的方向。從源4發(fā)射的地震能量在地球內(nèi)經(jīng)過反射，上行反射波入射到接收器5上。接收器5包括至少兩個正交地震檢波器，表示接收器的方框(標(biāo)為“rec”)中的箭頭表示接收器所測量的質(zhì)點運動的分量。當(dāng)接收器在地球表面記錄波場時，其不加區(qū)別地記錄P波和S波。
圖4(a)圖示了四種情況(i)震源產(chǎn)生水平質(zhì)點運動，接收器記錄在該水平方向的質(zhì)點運動；(ii)震源產(chǎn)生水平質(zhì)點運動，接收器記錄接收到的質(zhì)點運動的垂直分量；(iii)震源產(chǎn)生垂直質(zhì)點運動，接收器記錄接收到的質(zhì)點運動的水平分量；(iv)震源產(chǎn)生垂直質(zhì)點運動，接收器記錄接收到的質(zhì)點運動的垂直分量。
圖4(b)圖解了對圖4(a)所記錄的數(shù)據(jù)應(yīng)用去震源特征和彈性去多次反射過程后的結(jié)果。也就是，圖4(b)圖示了按照圖2的步驟S1到S5對圖4(a)的地震勘探方案采集到的地震數(shù)據(jù)進(jìn)行處理后的結(jié)果。該處理的效果是從采集到的地震數(shù)據(jù)中消除蓋層的響應(yīng)(在此例中，也就是地球表面的響應(yīng))。該處理還消除地震源5的輻射特性的影響。震源在這里是地震能量的點源。一般，能量可以向上傳播而超過震源和接收器，因為用無限的均一介質(zhì)取代了所述蓋層，這如圖4(b)所示。
如上所述，震源4發(fā)射P波和S波。在圖4(a)和4(b)中，接收器5不加區(qū)別地記錄S波和P波。
圖4(c)圖解了將接收器側(cè)的波場分解應(yīng)用于圖4(b)所示的理想化假想地震勘探方案的地震數(shù)據(jù)的結(jié)果。這具有用模擬的純P波檢測器(圖4(c)(i)和圖4(c)(iii))，或者用模擬的純S波檢測器(圖4(c)(ii)和圖4(c)(iv))取代圖4(a)和4(b)中用作接收器的地震檢波器的效果。
最后，圖4(d)圖示了進(jìn)一步應(yīng)用源側(cè)波場分解的結(jié)果。這具有用模擬的純P波震源(圖4(c)(i)和圖4(c)(iii))，或者用模擬的純S波震源(圖4(c)(ii)和圖4(c)(iv))取代圖4(a)到4(d)中同時發(fā)射P波和S波的地震源的效果。
圖4(d)圖解了應(yīng)用圖2所示實施例的步驟S1到S6的結(jié)果。可以看到，該方法消除了蓋層響應(yīng)，消除了震源的輻射特征，將采集到的數(shù)據(jù)分解成了PP、PS、SP和SS波至。
圖5圖解了將本發(fā)明的方法應(yīng)用于洋底地震勘探的結(jié)果。在圖5(a)所示的勘探方案中，分別在海床上設(shè)置一個地震源4和一個地震接收器5。表示震源的方框中的箭頭表示震源產(chǎn)生的源運動的方向。接收器5包含至少兩個正交地震檢波器，表示接收器的方框中的箭頭表示由地震檢波器測得的質(zhì)點運動的分量。(用在實際的洋底地震勘探中的接收器也可以包括一個水聽器，但這在圖5中未圖示)。源運動的水平和垂直方向以及地震檢波器測得的質(zhì)點運動的水平和垂直分量的四個組合示于圖5(a)(i)到圖5(a)(iv)中，對應(yīng)于圖4(a)(i)到圖4(a)(iv)所示的四種情況。
從圖5(a)可以看到，震源發(fā)射下行P波和S波，它們在地球中反射從而入射到接收器5上。圖5(a)還圖解了地震能量的一種可能的水層多次反射路徑，其中，上行P波在水面被反射從而入射到接收器5上。
接收器5中的地震檢波器不加區(qū)別地測量入射的S波和P波。
圖5(b)圖解了對圖5(a)所示的地震勘探方案中獲取的地震數(shù)據(jù)應(yīng)用本發(fā)明的去震源特征和彈性去多次反射方法的效果。可以看到，消除了水層多次反射路徑，因為該發(fā)明從采集到的地震數(shù)據(jù)中消除了蓋層(在此例中，也就是水層)的響應(yīng)。本發(fā)明的方法還消除地震源的輻射特性的影響。圖5(a)的理想化地震勘探方案中的震源是一個點力源。
圖5(b)圖解了將圖1所示的步驟S1到S5應(yīng)用于圖5(a)的地震勘探方案獲得地震數(shù)據(jù)的效果。應(yīng)用本發(fā)明的其它步驟，也就是，應(yīng)用接收器側(cè)分解和源側(cè)分解會提供與圖4(c)(i)到圖4(d)(iv)所示相同的結(jié)果。
圖6(a)和6(b)圖解了將本發(fā)明應(yīng)用于接收器和震源都布置地球內(nèi)的陸地地震勘探的結(jié)果。在圖6(a)所示的地震勘探方案中，地震源4和地震接收器5都布置在地球內(nèi)，接收器的深度大于震源。圖6(a)圖解了地震能量從震源到接收器的主要路徑，其中，下行S波和P波被目標(biāo)反射體3反射而入射到接收器5上。圖6(a)還示出了不希望有的多次反射路徑，其中，上行S波和P波在震源上方的層中反射從而入射到接收器上。
在圖6(a)中，表示震源4的方框中的箭頭表示震源產(chǎn)生的源運動的方向。表示接收器的方框中的箭頭表示包含在接收器中的地震檢波器測到的質(zhì)點運動的分量。跟圖4和圖5中一樣，水平和垂直源運動以及水平和垂直地震檢波器取向的四個組合示于圖6(a)(i)到圖6(a)(iv)中。
圖6(b)圖解了將本發(fā)明應(yīng)用于在圖6(a)的地震勘探方案中獲取的地震數(shù)據(jù)的效果。與在圖4和圖5中一樣，將本發(fā)明的去震源特征和去多次反射方法應(yīng)用于在圖6(a)的地震勘探方案中獲得的數(shù)據(jù)的效果是從獲取的地震數(shù)據(jù)中消除蓋層的影響。同時也消除了地震源5的輻射特征的效應(yīng)。這樣，圖6(b)圖解了一種理想化的地震勘探，其中消除了蓋層響應(yīng)(相當(dāng)于均一蓋層)，并且地震源4是一個點力源。圖6(b)圖解了將圖2的步驟S1到S5應(yīng)用于在圖6(a)的實際地震勘探方案中獲取的地震數(shù)據(jù)的結(jié)果。
應(yīng)用接收器側(cè)分解和源側(cè)分解的效果分別對應(yīng)于圖4(c)(i)到圖4(d)(iv)所示的結(jié)果。
圖7(a)到圖7(d)圖示了原始模擬地震數(shù)據(jù)。該數(shù)據(jù)是針對陸地地震勘探模擬的，其中，震源和接收器均位于地球表面。使用B.L.Kennet在“Seismic Wave Propagation in Stratefied Media”CambridgeUniversity Press(1993)中提出的類型的反射率碼(reflectivity code)模擬地震數(shù)據(jù)。用來模擬地震數(shù)據(jù)的模型由平坦自由表面構(gòu)成，其上布置震源和接收器。在該自由表面之下是一個300m厚的層，具有平坦的下界面。該層是一個各向同性彈性層，P波速率為1800m/s，S波速率為800m/s，密度為200kg/m3。在該層之下是一個各向同性的、無限的半空間，P波速率為2500m/s，S波速率為1500m/s，密度為2200kg/m3。由于該模型是水平分層的，有必要考慮垂直源運動和一個水平源運動，并考慮一個記錄質(zhì)點運動的垂直分量或者水平分量的地震檢波器。這樣，就只需要考慮2C×2C(4C)的情況，而不是完全的3C×3C(9C)的情況。
在圖7(a)到10(d)中，標(biāo)記Vz和Vx分別表示在接收器處的垂直質(zhì)點運動和水平質(zhì)點運動。Fz和Fx分別表示記錄質(zhì)點運動的垂直分量和水平分量的地震檢波器。
圖8(a)到8(d)分別圖解了圖7(a)到7(d)的原始數(shù)據(jù)的上行成分。圖9(a)到9(d)分別圖解了圖7(a)到7(d)的原始地震數(shù)據(jù)的下行成分。圖8(a)到圖9(d)是通過應(yīng)用下面的等式(A-28)、(A-29)、(A-32)和(A-33)獲得的。在圖8(a)到圖9(d)中仍然存在自由表面的效應(yīng)。
圖10(a)到10(d)圖解了消除自由表面和震源特征的效應(yīng)對圖7(a)到7(d)的地震數(shù)據(jù)的效果。這是通過應(yīng)用本發(fā)明的去震源特征和去多次反射方法而實現(xiàn)的。也就是，圖10(a)到10(d)分別圖示了將圖2的步驟S1到S5應(yīng)用于圖7(a)到7(d)所示的原始地震數(shù)據(jù)的效果?？梢钥吹?，圖7(a)到7(d)中與多次反射有關(guān)的地震波至已經(jīng)被從圖10(a)到10(d)所示的地震數(shù)據(jù)中消除了。
圖10(a)到10(d)的地震數(shù)據(jù)中的波至包括所有的PP、PS、SP和SS波至，以及相關(guān)的首波(head-waves)。
圖11(a)到11(d)圖解了將接收器側(cè)分解和源側(cè)分解分別應(yīng)用于圖10(a)到10(d)所示的數(shù)據(jù)的效果。這樣的效果是將地震數(shù)據(jù)分解為PP、SP、PS和SS波至。示于圖11(a)到11(d)的每一個附圖的地震數(shù)據(jù)只包括一個事件以及相關(guān)的首波。
上面針對結(jié)合了去多次反射和去震源特征的方法描述了本發(fā)明。在備選實施例中，本發(fā)明提供僅去震源特征的方法。在此備選方法中，從接收器處地震能量的初波至(initial event)的下行成分計算所述算子。該初波至構(gòu)成從震源不經(jīng)過任何反射，例如沿著圖1(a)中的路徑6，直接傳播到接收器的地震能量。波場的入射部分例如可以由A.Weglein和B.G.Secrest在“Wavelet estimation for amulti-dimensional acoustic or elastic earth”，Geophysics Vol.55，pp902-913(1990)中針對聲學(xué)或彈性地層(acoustic or elastic earth)提出的方法來確定。
另外，如果已知或者可以估計多分量源的震源特征/小波(signature/wavelet)，則本發(fā)明也可用來僅提供去多次反射方法。這樣的方法的主要步驟是(a)根據(jù)本發(fā)明的方法獲得去震源特征和去多次反射組合算子，(b)根據(jù)地震源的已知的或者計算出的輻射特征獲得去震源特征算子，和(c)獲得與在步驟(b)獲得的去震源特征算子相結(jié)合會得到在步驟(a)獲得的去震源特征和去多次反射組合算子的去多次反射算子。
下面對本發(fā)明的一個實施例給出詳細(xì)的數(shù)學(xué)說明。
考慮由面S＝∑+SR圍成的空間V，其具有向外指的法向矢量n，如圖3(a)到3(c)所示?！剖且粋€位于多分量接收器xr14(位于深度zr處)上方無窮小距離處的深度zr-的水平面。為了簡化分析，假設(shè)在深度z以及該深度z下方無窮薄區(qū)域內(nèi)的固體介質(zhì)是均一的、各向同性的。本發(fā)明不需要關(guān)于蓋層(在深度z＜zr處)和地下層(在深度z＞zr)的特性的知識，蓋層和地下層都可以是任意的非均一、各向異性和滯彈性介質(zhì)。當(dāng)接近一個物理表面獲取地震數(shù)據(jù)時，接收器總是位于該物理表面下方無窮小距離處。物理表面一般是地球表面(假設(shè)為固體/真空界面)、洋底(認(rèn)為是固體/流體邊界)或者地球內(nèi)的另一種界面(假定為固體/固體界面)。SR是一個半徑為R的半球面。笛卡爾坐標(biāo)表示為x(ξ，x3)，其中ξ＝(x1，x2)。為了概念表達(dá)的方便，也使用x3＝z。x3軸的正方向向下，平行于n。x1、x2軸在∑平面中。
首先，在實際地震勘探的多分量源和多分量接收器數(shù)據(jù)(這些數(shù)據(jù)包含接收器上方的介質(zhì)的散射響應(yīng))以及在沒有散射響應(yīng)的理想化地震勘探中應(yīng)獲得的希望的多分量源和多分量接收器數(shù)據(jù)之間建立一個積分關(guān)系。
在勘探中使用的地震源假設(shè)單獨地產(chǎn)生三個正交源運動。希望的多分量數(shù)據(jù)是在理想化的假想多分量地震勘探(使用三個單獨作用的點力，當(dāng)接收器上方的介質(zhì)均一、向上延伸到無窮遠(yuǎn)并具有與接收器深度的介質(zhì)相等的參數(shù)時，這三個點力具有希望的、相等的震源特征(子波，signature))中應(yīng)記錄到的數(shù)據(jù)。這樣，在理想化勘探中的蓋層是彈性的各向同性半空間。自然，接收器深度以下的地質(zhì)對于實際的和假想的地震勘探是一樣的。
表1進(jìn)入分別與圖3(a)到圖3(c)所示的狀態(tài)P，H和 相關(guān)的彈性動力學(xué)方程的源和場的符號 在實際的勘探中，地震源是單方向的，或者是振動器或者是力(force)。一個振動器源由一個垂直振動器(張應(yīng)力源)和兩個水平取向的振動器(剪應(yīng)力源)構(gòu)成。位于位置xs的振動器不需要是理想點源。其震源特征不需要是已知的。在假想勘探H和 中，震源是分別在n和m方向的力的點源，震源特征為 考慮按照圖3(a)所示配置的地震勘探，為了方便，下面稱為“實際”地震勘探。由在中心坐標(biāo)xs處以方向n作用、源強度和輻射特征未知的震源引起的，在∑下方無窮小距離處的接收器位置xr處記錄的質(zhì)點速度矢量的第m個分量表示為vmn。類似地，垂直推力矢量的第m個分量表示為smn。實際地震勘探(簡記為“狀態(tài)P”)的源變量和場變量列于表1中。注意，面∑可以也可以不與一個物理面重合。在記錄剛好在地球表面之下進(jìn)行時，或者剛好在洋底之下進(jìn)行時，面∑顯然與物理面重合。
計劃求解的理想波場 和 是當(dāng)接收器高度以上的介質(zhì)是如圖3(b)所示的半空間時，介質(zhì)的來自具有理想震源特征 的三個正交點力的響應(yīng)。在這種情況下，面∑是非物理邊界。對于位于∑上xr-的點力，在位置xr即∑下無窮小距離處記錄所希望的質(zhì)點速度和垂直推力矢量響應(yīng)。在表1中，列出了對于這種假想地震試驗(簡略地標(biāo)為“狀態(tài)H”)的源和場變量。
為了在實際狀態(tài)P和假想狀態(tài)H地震試驗之間建立積分關(guān)系，引入假想“狀態(tài) ”。這是另一種理想化地震勘探，對應(yīng)于圖3(b)所示的理想化地震勘探H，但是源5和接收器6相互調(diào)換了。這樣，在假想“狀態(tài) ”中的波場 和 與狀態(tài)H中的波場互易，遵循互反關(guān)系v^nm(xr-|xr)=v~mn(xr|xr-),---(1)]]>s^nm(xr-|xr)=s~mn(xr|xr-),---(2)]]>這樣， 和 是震源特征為 在圖3(c)所示的面∑下無窮小距離處的位置xr處在m方向作用的點力產(chǎn)生的在面∑上在位置xr-的響應(yīng)。如同在希望的理想化狀態(tài)H中一樣，面∑是一個人為的非物理邊界。下面將利用這樣的事實在假想狀態(tài) 的勘探中獲取的地震數(shù)據(jù)由僅僅從∑以下的地下層散射的上行波至構(gòu)成。另外，由于源在接收器下方無窮小距離處，從源到接收器的直接波模式(direct wavemode)是上行波至。狀態(tài) 的源和場變量列于表1中。
下一步，在面∑+SR→∞圍成的空間V中應(yīng)用貝蒂互易定理(Betti’sreciprocity theorem)，其中SR→∞位于無窮遠(yuǎn)處。例如類似于L.Amundsen et al.在“Elimination of free-sruface related multipleswithout the need of the source wavelet”，Gephpysics，Vol66，pages327-341(2000)中的方案，獲得描述狀態(tài)P和狀態(tài) 之間的關(guān)系的頻域積分方程 面SR當(dāng)R→∞時對于積分沒有貢獻(xiàn)(這是輻射條件，見Y.H.Paoand V.Varatharajulu in“Huygen’s Principle，radiation conditions andintegral formulae for the scattering of elastic waves”，J.Acoust.Soc.Am.Vol.59，pp1361-1371(1976))。通過對∑指定合適的邊界條件，可以簡化等式(3)。在實際狀態(tài)P中，vin和sin是上行波和下行波的和vin=vin(u)+vin(d)---(4)]]>sin=sin(u)+sin(d)---(5)]]>而在假想狀態(tài) 中， 和 是純上行場v^in=v^in(u);v^in(d)=0---(6)]]>s^in=s^in(u);s^in(d)=0---(7)]]>在明確已知上行波和下行波以及它們與垂直推力和質(zhì)點速度矢量的關(guān)系的情況下，通過分析水平波數(shù)域中的問題，這些邊界條件能最為方便地引入等式(3)中。利用帕斯維爾恒等式(Parceval identity)∫-∞∞dξf(ξ)h(ξ)=1(2π)2∫-∞∞dκF(-κ)H(κ),---(8)]]>其中，波數(shù)矢量κ＝(k1，k2)與ξ＝(x1，x2)共軛，從等式(3)得到a~vmn(xr|xs)=1(2π)2∫dκ[S^im(κ,zr-|xr)Vin(-κ,zr|xs)-V^im(κ,zr-|xr)Sin(-κ,zr|xs)].---(9)]]>從求和縮寫法到向量記法，等式(9)寫為a~vmn(xr|xs)=1(2π)2∫dκ[S^mT(κ,zr-|xr)Vn(-κ,zr|xs)-V^mT(κ,zr-|xr)Sn(-κ,zr|xs)],---(10)]]>其中，S^T=(S1,S2,S3)]]>和V^T=(V1,V2,v3)]]>分別是波數(shù)域垂直推力矢量和質(zhì)點速度矢量，上標(biāo)T表示轉(zhuǎn)置。如附件所示，由于假想狀態(tài) 場 和 僅由上行波模式構(gòu)成，它們關(guān)聯(lián)如下S^(κ)=L^SV(κ)V^(κ),---(11)]]>其中， 是一個取決于沿著接收器的分布的局部介質(zhì)參數(shù)ρ、α、β的3×3矩陣。 的要素由附件中的等式(A-17)給出。將等式(11)插入等式(10)中得到a~vmn(xr|xs)=1(2π)2∫dκV^mT(κ,zr-|xr)Rn(D)(-κ,zr|xs),---(12)]]>其中R(D)=L^SVV-S=G-1V(D)---(13)]]>被解釋為由于接收器高度上方的介質(zhì)中的散射導(dǎo)致的蓋層響應(yīng)。另外，矢量V(D)=[V1(D),V2(D),V3(D)]T---(14)]]>包含每一個質(zhì)點速度分量V1、V2和V3上的下行波模式(wavemode)的元素。通常，對于每一個炮點位置，根據(jù)等式(A-29)、(A-31)和(A-33)從質(zhì)點速度矢量和垂直推力矢量在慢度(slowness)和波數(shù)域計算V(D)。所得到的下行反射響應(yīng)的元素為R1(D)=2ρp2+qαqβ{[qα-p22β2(qα-qβ)]V1(D)+p1p2β2(qα-qβ)V2(D)}---(15)]]>R2(D)=2ρp2+qαqβ{p1p2β2(qα-qβ)V1(D)+[qα-p12β2(qα-qβ)]V2(D)}---(16)]]>R3(D)=2ρqβp2+qαqβV3(D).---(17)]]>重新引入求和縮寫法，得到av~mn(xr|xs)=1(2π)2∫dκV^im(κ,zr-|xr)Rin(D)(-κ,zr|xs).---(18)]]>利用帕斯維爾恒等式，在空間域中得到av~mn(xr|xs)=∫ϵdS(ξ)v^im(x|xr)rin(d)(x|xs),---(19)]]>其中rin(d)為Rin(D)的傅立葉變換。
等式(19)給出了假想狀態(tài) 中的場 與狀態(tài)P下記錄到的場vnm之間的積分關(guān)系。注意，除了位置之外，為得到 的積分等式(19)，沒有使用實際源及其輻射特征的任何信息，也沒有使用面∑上方的實際蓋層或者接收器高度下方的實際地下層的特性的任何信息。根據(jù)互易性直接確定假想狀態(tài)H的場 注意，等式(19)中的rin(d)的逆可以解釋為用作下述用途的多維算子(i)消除實際震源特性的決定性去震源特征算子，和(ii)從實際數(shù)據(jù)中消除蓋層響應(yīng)的決定性多次反射衰減算子。
在假想狀態(tài)消除入射波場在假想試驗中所希望的場 可以分解為一個從源到接收器向上傳播的入射波場 以及一個從地下層向上散射的波場 v^nm=v^nm(inc)+v^nm(sc).---(20)]]>所述入射波場在均一介質(zhì)中傳播，因此與震源位置無關(guān)，其可以被表示為V^(inc)(zr-|zr)=a~Gexp(-iκ·ξr),---(21)]]>其中，G是格林張量(Green′s tensor)，相當(dāng)于從下面的等式(A-39)得到的張量。注意到(V^(inc))T(κ,zr-|xr)R(D)(-κ,zr|xs)=a~GT(κ)G-1(-κ)V(D)(-κ)exp(-iκ·ξr),---(22)]]>=a~V(D)(-κ)exp(-iκ·ξr)]]>并且vmn(d)(xr|xs)=1(2π)2∫dκexp(-iκ·ξr)Vmn(D)(-κ,zr|xs),---(23)]]>我們得到a~vmn(u)(xr|xs)=1(2π)2∫dκV^im(sc)(κ,zr-|xr)Rin(D)(-κ,zr|xs).---(24)]]>利用帕斯維爾恒等式，在空間域中得到a~vmn(u)(xr|xs)=∫sdS(ξ)v^im(sc)(x|xr)rin(D)(x|xs).---(25)]]>等式(25)給出了所要找的在假想狀態(tài) 勘探中的散射波場 和狀態(tài)P下的上行和下行場vmn(u)和vmn(d)之間的積分關(guān)系。假想狀態(tài)H下的散射場 由互易關(guān)系(reciprocity relation)得到。
等式(25)是所需要的散射場的第一類弗雷德霍姆積分方程，通過保持固定的接收器坐標(biāo)同時改變源坐標(biāo)，可以得到一組方程，由之解出 波數(shù)域解在源坐標(biāo)ξs和接收器坐標(biāo)ξr上對等式(25)進(jìn)行傅立葉變換，得到a~Vmn(U)(κr,zr|κs,zs)=1(2π)2∫dκV^im(sc)(κ,zr-|κr,zr)Rin(D)(-κ,zr|κs,zs).---(26)]]>通過保持與接收器坐標(biāo)共軛的波數(shù)不變，而使與震源坐標(biāo)共軛的波數(shù)變化，得到一組方程，由之可解出 下行蓋層響應(yīng)場中的正波數(shù)與所希望的場中的負(fù)波數(shù)之間的匹配(coupling)反映了兩個場之間的自相關(guān)過程。
分解為一次PP、PS、SP和SS波響應(yīng)在貝蒂去震源特征和彈性去多次反射步驟之后是彈性波場分解步驟，將等式(25)給出的多分量震源和多分量接收器的去震源特征和彈性去多次反射數(shù)據(jù) 分解為應(yīng)當(dāng)從壓力波和剪切波源和接收器記錄到的一次PP、PS、SP和SS波響應(yīng)。假想數(shù)據(jù)的接收器和源坐標(biāo)分別記為xrh和xsh，其中xrh是面∑上的點，xsh＝xr在其下方無窮小距離處。令矩陣 包含在接收器處來自方向i的力源的上行勢能(up-going potential)，U^=U^P1U^P2U^P3U^SV1U^SV2U^SV3U^SH1U^SH2U^SH3.---(27)]]>令U為包含在接收器處由P、SV和SH源激發(fā)的上行P、SV和SH波勢的矩陣U=UPPUPSVUPSHUSVPUSVSVUSVSHUSHPUSHSVUSHSH.---(28)]]>這樣，該彈性源-接收器波場分解可以分成兩個計算操作。首先，將多分量數(shù)據(jù) 在水平接收器坐標(biāo)ξrh上進(jìn)行傅立葉變換。在波數(shù)域κr中，通過用接收器分解矩陣R左乘 實現(xiàn)到上行壓力波場和剪切波場的接收器側(cè)波場分解U^(κr,zrh|xsh)=R(κr)V^(sc)(κr,zrh|xsh),---(29)]]>其中，
R=1p2+qαqβp1p2-qβ-p1qαp-p2qαp-p-p2p(p2+qαqβ)p1p(p2+qαqβ)0,---(30)]]>其中p＝κr/ω。傅立葉逆變換給出接收器處測量的上行壓力波場和剪切波場 然后，在水平源坐標(biāo)ξsh上對多分量數(shù)據(jù) 進(jìn)行傅立葉變換。在波數(shù)域κs中，通過用源分解算子S對數(shù)據(jù)進(jìn)行右乘而實現(xiàn)到下行壓力波場和剪切波場的源側(cè)波場分解U(xrh|κs,zsh)=U^(xrh|κs,zsh)S(κs),---(31)]]>其中S=-iω-1ρp2+qαqβp1-p1qαp-p2qαpβ2(p2+qαqβ)p2-p2qαpp1qαpβ2(p2+qαqβ)-qβ-p0,---(32)]]>其中p＝κs/ω。傅立葉逆變換給出相當(dāng)于沒有蓋層、具有單分量壓力和剪切波源以及單分量壓力和剪切波接收器的假想勘探的數(shù)據(jù)的數(shù)據(jù)U。
在(κr，κs)空間中，源-接收器分解為U(κr,zrh|κs,zsh)=R(κr)V^(sc)(κr,zrh|κs,zsh)S(κs),---(33)]]>注意，分解算子只是被數(shù)據(jù)乘的矩陣。分解所需的參數(shù)是接收器深度的密度、P波和S波速度。
貝蒂去卷積和源-接收器波場分解水平分層介質(zhì)(“1.5D介質(zhì)”)在水平分層介質(zhì)中，相對于水平源位置ξs，響應(yīng)在橫向上的偏移不變，因此表達(dá)式大大簡化。為此，水平分層介質(zhì)ls有時稱為“1.5D介質(zhì)”。對于實際數(shù)據(jù)和真實地球結(jié)構(gòu)而言，水平分層介質(zhì)的假設(shè)常常得到令人滿意的結(jié)果。
考慮等式(24)，可以設(shè)ξs＝0，這樣V^im(sc)(κ,zr-|xr)=V^im(sc)(κ,zr-|ξ=0,zr)exp(-iκ·ξr).---(34)]]>
相對于ξr對等式(24)進(jìn)行傅立葉變換，并互換積分(interchanging integrals)，得到a~VmnU(κr,zr|ξs=0,zs)=1(2π)2∫-∞∞dκV^im(sc)(κ,zr-|ξ=0,zr)---(35)]]>×Rin(D)(-κ,zr-|ξs=0,zs)∫-∞∞dξrexp(-iξr·(κ+κr))]]>最后一項積分是狄拉克德爾塔函數(shù)。在波數(shù)上積分，使用狄拉克德爾塔函數(shù)特性，并用κ對κr重命名，得到a~VmnU(κ,zr|ξs=0,zs)=V^im(sc)(-κ,zr-|ξ=0,zr)Rin(D)(κ,zr-|ξs=0,zs)---(36)]]>該結(jié)果用矩陣形式表示為V^11(sc)V^21(sc)V^31(sc)V^12(sc)V^22(sc)V^32(sc)V^13(sc)V^23(sc)V^33(sc)(-κ,zr-|zr)=a~V11(U)V12(U)V13(U)V21(U)V22(U)V23(U)V31(U)V32(U)V33(U)(κ,zr|zr)R11(D)R12(D)R13(D)R21(D)R22(D)R23(D)R31(D)R32(D)R33(D)(κ,zr|zr)-1---(37)]]>利用場特性 顯然，通過推廣的頻譜去卷積，獲得所希望的場的散射部分V^11(sc)-V^21(sc)-V^31(sc)-V^12(sc)V^22(sc)-V^32(sc)-V^13(sc)-V^23(sc)V^33(sc)(κ,zr-|zr)=a~V11(U)V12(U)V13(U)V21(U)V22(U)V23(U)V31(U)V32(U)V33(U)(κ.zr|zr)R11(D)R12(D)R13(D)R21(D)R22(D)R23(D)R31(D)R32(D)R33(D)(κ,zr|zr)---(39)]]>從等式(39)可以得出結(jié)論，通過場本身和包含質(zhì)點速度矢量的下行部分的蓋層反射響應(yīng)之間的決定性的頻譜去卷積，獲得了所希望的場的分量。
下一步是將貝蒂去卷積數(shù)據(jù)分解為一次PP、PS、SP和SS波響應(yīng)。這按照“分解為一次PP、PS、SP和SS波響應(yīng)”一節(jié)所述進(jìn)行。
注意，后跟源-接收器分解的1.5D貝蒂去卷積方案可以實現(xiàn)為τ-p或者頻率-波數(shù)域算法。在τ-p域中，對于每一個p記錄跡執(zhí)行去震源特征、多次反射衰減以及源-接收器分解聯(lián)合過程。在頻率-波數(shù)域中，對于每一個頻率和波數(shù)的組合進(jìn)行該聯(lián)合過程。
本發(fā)明的應(yīng)用如上所述，本發(fā)明可應(yīng)用于許多不同的地震勘探方案。下面簡要討論三種尤為重要的勘探類型。
1.陸地表面地震探測本發(fā)明可應(yīng)用于3C×3C(或者說9C)陸地地震勘探，其中，獨立產(chǎn)生三個正交源運動(兩個水平一個垂直)，并用三個正交地震檢波器記錄地震波場，測量質(zhì)點速度矢量的兩個水平和一個垂直分量。
對于在自由表面進(jìn)行的測量，上/下分離大大簡化了所有推力S1、S2和S3變?yōu)榱?，從而下述等?A-28)到(A-33)中的推力項變?yōu)榱恪?br> 2.海床地震探測本發(fā)明可以應(yīng)用于在3C×4C(或者說12C)海床地震試驗中采集的地震數(shù)據(jù)。在所述海床地震試驗中，或者在洋底或者在水柱中分別產(chǎn)生三個正交源運動(兩個水平一個垂直)，在接收器位置用剛好布置在海床以下的三個地震檢波器和剛好布置在海床以上的一個水聽器記錄地震波場。
對于在海床的流體/固體界面上進(jìn)行的測量，上/下行分離大大簡化了推力S1和S2變?yōu)榱?，從而在下面的等?A-28)到(A-33)中，相應(yīng)的推力項變?yōu)榱?。海底上的水聽器測量是S3的負(fù)。
需要注意，通過例如試驗傳統(tǒng)的氣炮(airgun)源，能在水中產(chǎn)生等效于力源的震源。例如，相應(yīng)于垂直點力的響應(yīng)可以通過在同一位置但是在稍微不同的深度獲取同一個源的兩個記錄來獲得。然后，使用有限差近似，可以計算由垂直方向的壓力梯度源導(dǎo)致的響應(yīng)。運動方程(牛頓第二定律)表明，這等效于垂直點力。對于兩個水平方向同樣如此。
3.井孔地震探測在3C×6C(或者說18C)井孔地震試驗中，在陸地、海底或者水柱中單獨地產(chǎn)生三個正交源運動，其中，質(zhì)點速度的三個分量和垂直推力的三個分量是已知的。
在這種勘探方案中，沒有簡化上/下分離步驟的邊界條件。但是，通過在接收器組中增加水聽器，可以獲得附加約束條件(一個水聽器測量質(zhì)點運動的散度)。
質(zhì)點速度和垂直推力矢量以及上行和下行波矢量之間的波數(shù)域基本關(guān)系這部分考慮水平分層的彈性地層。在無源區(qū)，彈性波的傳播使用運動方程和虎克定律(彈性本構(gòu)關(guān)系，elastic constitutive relation)來描述。這可以寫為下列形式的一組一階常微分方程(ordinarydifferential equations)∂3B=-iωAB---(A-1)]]>其中，場矢量B定義為B＝(VT，ST)T(A-2)其中質(zhì)點速度矢量VT＝[V1，V2，V3]和垂直推力矢量ST＝[S1，S2，S3]。彈性系統(tǒng)矩陣為下述形式A=00p11μ0000p201μ0λλ+2μp1λλ+2μp20001λ+2μρ-θp12-μpσpσ-θp1p2000λλ+2μp1-θp1p2ρ-θp22-μpσpσ000λλ+2μp200ρp1p20---(A-3)]]>其中θ＝μ(3λ+2μ)/(λ+2μ)，p2＝p12+p22。
為了表示的方便，省略了各種量與頻率、波數(shù)、深度等的相關(guān)性的明確表達(dá)。例如，由位置xs處的點源產(chǎn)生的、在深度x3記錄的質(zhì)點速度矢量v(ξ，x3，ω，xs)在波數(shù)域中被表示為V或者V(x3)，但應(yīng)理解為V＝V(x3)＝V(ξ，x3，ω，xs)。當(dāng)需要時，會表示出對水平慢度矢量(slownessvector)p＝κ/ω＝(p1，P2)的相關(guān)性。
上行和下行波對于分層地層中彈性場向上行和下行波的分解，需要對給定的波數(shù)和頻率找出系統(tǒng)矩陣A的本征值和本征向量。場向量B可以分解為上行(U)和下行(D)波W＝[UT，DT]T(A-4)其中UT＝[UP，USV，USH]，DT＝[DP，DSV，DSH]，上述分解是通過下述線性變換B＝LW (A-5)其中L是A的局部本征向量矩陣(即，L的每一列是一個本征向量)。等式(A-5)從其上行和下行成分描述了波場B的組成。
給定逆本征向量矩陣L-1，可以通過計算W＝L-1B. (A-6)來計算上行和下行波。
等式(A-6)描述了將波場B分解為上行和下行P波和S波。在進(jìn)行一些簡單然而冗長的計算之后，獲得合成矩陣(compositionmatrix)L=L(p)=LVU(p)-LVU(-p)LSU(p)LSU(-p),---(A-7)]]>和分解矩陣L-1=L-1(p)=LSUT(p)LVUT(p)-LSUT(-p)LVUT(-p),---(A-8)]]>其中
LVU=12-p11ρqαp1pqβρp2p1μqβ-p21ρqαp2pqβρ-p1p1μqβqαρp1ρqβ0,---(A-9)]]>LSU=12-2μp1qαρp1p(ρ-2μp2)1ρqβp2pμqβ-2μp2qαρp2p(ρ-2μp2)1ρqβ-p1pμqβ(ρ-2μp2)1ρqα2μpqβρ0.---(A-10)]]>在無源均一固體中，上行波和下行波滿足以下微分方程∂3UP=-iωqαUP,---(A-11)]]>∂3USV=-iωqβUSV,---(A-12)]]>∂3USH=-iωqβUSH,---(A-13)]]>∂3DP=iωqαDP,---(A-14)]]>∂3DSV=iωqβDSV,---(A-15)]]>∂3DSH=iωqβDSH.---(A-16)]]>當(dāng)下行波變?yōu)榱銜r的推力-質(zhì)點速度矢量關(guān)系利用關(guān)系(A-5)，使D＝0，得到垂直推力矢量和質(zhì)點速度矢量之間的簡單關(guān)系S=L^SVV,---(A-17)]]>其中
L^SV=L^SVL^VU-1.---(A-18)]]>L^SV=ρp2+qαqβqα-p22β2(qα-qβ)p1p2β2(qα-qβ)p1[1-2β2(p2+qαqβ)p1p2β2(qα-qβ)qα-p12β2(qα-qβ)p2[1-2β2(p2+qαqβ)-p1[1-2β2(p2+qαqβ)]-p2[1-2β2([p2+qαqβ)]qβ---(A-19)]]>比例換算為質(zhì)點速度和垂直推力的上行波和下行波不是唯一地限定上行和下行P、SV和SH波?？梢詫⒎至堪幢壤龘Q算，使其具有質(zhì)點速度或者推力的維度，并且其和給出質(zhì)點速度的一個分量或者垂直推力的一個分量。令Vi(Up)表示Vi上的上行波，Vi(Usv)表示Vi上的SV波，等等，則有Vi=Vi(UP)+Vi(USV)+Vi(USH)+Vi(DP)+Vi(DSV)+Vi(DSH)]]>對于垂直推力有類似的等式。另外，可以將上行分量加起來稱為總上行分量，并對下行分量進(jìn)行類似操作。Vi(U)定義為Vi上的上行波之和Vi(U)=Vi(UP)+Vi(USV)+Vi(USH),]]>Vi(D)定義為Vi上的下行波之和Vi(D)=Vi(DP)+Vi(DSV)+Vi(DSH).]]>Si上的總上行波為Si(U)=Si(UP)+Si(USV)+Si(USH),]]>Si上的總下行波為Si(D)=Si(DP)+Si(DSV)+Si(DSH).]]>另外，定義下述矢量V(U)=[V1(U),V2(U),V3(U)]T---(A-20)]]>V(D)=[V1(D),V2(D),V3(D)]T---(A-21)]]>S(U)=[S1(U),S2(U),S3(U)]T---(A-22)]]>S(D)=[S1(D),S2(D),S3(D)]T---(A-23)]]>
可以得出結(jié)論，這些矢量分別與原來限定的上行波和下行波U和D的矢量相關(guān)。
V(U)(p)＝LVU(p)U(p) (A-24)V(D)(p)＝-LVU(-p)D(p) (A-25)S(U)(p)＝LSU(p)U(p) (A-26)S(D)(p)＝LSU(-p)D(p)(A-27)V(D)和V(U)進(jìn)入實際試驗和假想試驗之間的積分關(guān)系，作為已知的場向量。
明確地說，V1上的總上行和下行波由V1、V3和S1給出V1(U)=12V1-p12qα[1-2β2(p2+qαqβ)]V3+12ρqα(p2+qαqβ)S1,---(A-28)]]>V1(D)=12V1-p12qα[1-2β2(p2+qαqβ)]V3-12ρqα(p2+qαqβ)S1.---(A-29)]]>注意，V1上的總上行和下行波與V2、S2以及S3無關(guān)。另外，V2上的總上行和下行波由V2、V3和S1給出V2(U)=12V2-p22qα[1-2β2(p2+qαqβ)]V3+12ρqα(p2+qαqβ)S1,---(A-30)]]>V2(D)=12V2+p22qα[1-2β2(p2+qαqβ)]V3-12ρqα(p2+qαqβ)S1,---(A-31)]]>從而與V1、S2和S3無關(guān)。V3上的總上行和下行波由V1、V2、V3和S3給出V3(U)=12V3+12ρqβ(p2+qαqβ)S3+12qβ[1-2β2(p2+qαqβ)](p1V1+p2V2),---(A-32)]]>V3(D)=12V3-12ρqβ(p2+qαqβ)S3-12qβ[1-2β2(p2+qαqβ)](p1V1+p2V2)---(A-33)]]>導(dǎo)出矩陣關(guān)系等式(A-18)和等式(A-19)的對稱關(guān)系給出[L^SVV-S](p)=LVU-T(-p)LSUT(-p)V(p)-S(p)]]>=LVU-T(-p)[LSUT(-p)V(p)-LVUT(-p)S(p)]---(A-34)]]>方括號中的最后一項利用等式(A-6)和(A-8)看成是根據(jù)下述等式從V和S計算出的下行波矢量D(p)=LSUT(-p)V(p)-LVUT(-p)S(p).---(A-35)]]>在前一小節(jié)中，已經(jīng)描述了通過等式(A-28)和(A-29)，D被關(guān)聯(lián)到下行質(zhì)點速度波矢量V(D)。反演該等式得到D(p)=-LVU-1(-p)V(D)(p),---(A-36)]]>將其代入等式(A-34)得到 =[LVU(-p)LVUT(-p)]1V(D)(p),---(A-37)]]>=G-1(p)V(D)(p)]]>其中G-1=12(LVULVUT)1.---(A-38)]]>顯然，可以得到G-1=2ρp2+qαqβqα-p22β2(qα-qβ)p1p2β2(qα-qβ)0p1p2β2(qα-qβ)qα-p12β2(qα-qβ)000qβ---(A-39)]]>在上面給出的理論中，假設(shè)了在理想化勘探中的震源是點力源。但是，本發(fā)明不限于使用作為點力源的理想化勘探震源，盡管上面給出的詳細(xì)理論需要作少量的修改以適應(yīng)不是點力源的理想化勘探震源。
在上面給出的詳細(xì)說明中，利用貝蒂定理得出去震源特征/去多次反射算子。但是，本發(fā)明不限于用貝蒂定理得出去震源特征/去多次反射算子?？梢詮膹椥詣恿W(xué)波動方程或者彈性動力學(xué)波動方程的表達(dá)得出本發(fā)明的去震源特征/去多次反射算子。貝蒂定理是彈性動力學(xué)波動方程的積分表達(dá)。彈性動力學(xué)波動方程的其它表達(dá)的例子包括彈性基爾霍夫積分或者彈性動力學(xué)表示定理(elastodynimicrepresentation theorem)，以及逆散射級數(shù)方法(inverse scatteringseries method)。
圖12是本發(fā)明的設(shè)備15的示意框圖。該設(shè)備能夠執(zhí)行本發(fā)明的方法。
該設(shè)備15包括一個具有程序存儲器17的可編程數(shù)據(jù)處理器16，該程序存儲器例如是只讀存儲器ROM的形式，存儲用于控制數(shù)據(jù)處理器17執(zhí)行本發(fā)明的方法的程序。該系統(tǒng)還包括例如用于存儲在沒有點源時必須保持的數(shù)據(jù)的非易失性讀寫存儲器18。有一個隨機存取存儲器(RAM)19為數(shù)據(jù)處理器提供一個“工作”或者“便箋式”存儲器。提供一個輸入設(shè)備20，例如用于接收用戶命令和數(shù)據(jù)。提供一個輸出設(shè)備21，例如用于顯示與該方法的進(jìn)程和結(jié)果有關(guān)的信息。輸出設(shè)備可以是例如打印機、可視顯示裝置或者輸出存儲器。
用于根據(jù)本發(fā)明的方法進(jìn)行處理的地震數(shù)據(jù)可以通過輸入設(shè)備20提供，或者，可選地，可以通過可機讀存儲裝置22提供。
用于操作系統(tǒng)并用于執(zhí)行這里所描述的方法的程序存儲在所述程序存儲器17中，后者可以實現(xiàn)為半導(dǎo)體存儲器，例如公知的ROM類型。但是，程序可以存儲到任何合適的存儲介質(zhì)，比如磁性數(shù)據(jù)載體17a(比如軟盤)或者CD-ROM 17b中。
權(quán)利要求
1.一種處理多分量地震數(shù)據(jù)的方法，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，并在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的，該方法包括下列步驟將采集到的地震數(shù)據(jù)分解為上行成分和下行成分；從所述采集到的地震數(shù)據(jù)的所述下行成分以及環(huán)繞該接收器的介質(zhì)的特性計算一個去震源特征和去多次反射算子。
2.如權(quán)利要求1所述的方法，還包括下述步驟利用所述去震源特征和去多次反射算子處理所述采集到的地震數(shù)據(jù)，從而衰減或者完全消除由多次反射導(dǎo)致的地震波至。
3.如權(quán)利要求2所述的方法，其中，處理所述采集到的地震數(shù)據(jù)的步驟包括利用所述去震源特征和去多次反射算子處理所述采集到的地震數(shù)據(jù)的上行成分。
4.如權(quán)利要求2或3所述的方法，其中，處理所述采集到的地震數(shù)據(jù)的步驟還包括為所述震源選定一個所希望的地震源特征。
5.如權(quán)利要求2、3或4所述的方法，還包括下述步驟將處理后的地震數(shù)據(jù)分解為P波和/或S波數(shù)據(jù)。
6.如權(quán)利要求1所述的方法，包括將所述采集到的地震數(shù)據(jù)分解為P波和/或S波數(shù)據(jù)的步驟，將所述采集到的地震數(shù)據(jù)分解為P波和/或S波數(shù)據(jù)的該分解步驟在將所述地震數(shù)據(jù)分解為上行成分和下行成分的步驟之前進(jìn)行。
7.如權(quán)利要求5或6所述的方法，其中，將所述地震數(shù)據(jù)分解為P波和/或S波數(shù)據(jù)的步驟包括在接收器側(cè)分解所述數(shù)據(jù)。
8.如權(quán)利要求5、6或7所述的方法，其中，將所述地震數(shù)據(jù)分解為P波和/或S波數(shù)據(jù)的步驟包括在震源側(cè)分解所述數(shù)據(jù)。
9.一種處理多分量地震數(shù)據(jù)的方法，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的，該方法包括下列步驟將所述采集到的地震數(shù)據(jù)分解為上行成分和下行成分；從所述采集到的地震數(shù)據(jù)的初始下行成分和環(huán)繞接收器的介質(zhì)的特性來計算一個去震源特征算子。
10.一種處理多分量地震數(shù)據(jù)的方法，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的，該方法包括下列步驟將所述采集到的地震數(shù)據(jù)分解為PP成分、PS成分、SP成分和SS成分；從所述采集到的地震數(shù)據(jù)的所述成分和環(huán)繞接收器的介質(zhì)的特性計算一個去震源特征和去多次反射算子。
11.如權(quán)利要求9或10所述的方法，還包括下述步驟利用所述去震源特征和去多次反射算子處理所述采集到的地震數(shù)據(jù)，以衰減或者消除由多次反射導(dǎo)致的地震波至。
12.一種為多分量地震數(shù)據(jù)計算去震源特征和去多次反射算子的方法，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的，該方法包括下列步驟將所述采集到的地震數(shù)據(jù)分解為上行成分和下行成分；從所述采集到的地震數(shù)據(jù)的初始下行成分和環(huán)繞接收器的介質(zhì)的特性來計算一個去震源特征和去多次反射算子。
13.一種為多分量地震數(shù)據(jù)計算去震源特征算子的方法，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的，該方法包括下列步驟將所述采集到的地震數(shù)據(jù)分解為上行成分和下行成分；從所述采集到的地震數(shù)據(jù)的初始下行成分和環(huán)繞接收器的介質(zhì)的特性來計算一個去震源特征算子。
14.一種為多分量地震數(shù)據(jù)計算去震源特征和去多次反射算子的方法，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的，該方法包括下列步驟將所述采集到的地震數(shù)據(jù)分解為PP成分、PS成分、SP成分和SS成分；從所述采集到的地震數(shù)據(jù)的所述成分和環(huán)繞接收器的介質(zhì)的特性計算一個去震源特征和去多次反射算子。
15.如前述任意權(quán)利要求所述的方法，還包括下述步驟在震源位置發(fā)射多分量地震能量；在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)。
16.一種處理多分量地震數(shù)據(jù)的設(shè)備，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，并在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的，該設(shè)備包括將采集到的地震數(shù)據(jù)分解為上行成分和下行成分的裝置；從所述采集到的地震數(shù)據(jù)的所述下行成分以及環(huán)繞該接收器的介質(zhì)的特性計算一個去震源特征和去多次反射算子的裝置。
17.如權(quán)利要求16所述的設(shè)備，還包括利用所述去震源特征和去多次反射算子處理所述采集到的地震數(shù)據(jù)，從而衰減或者完全消除由多次反射導(dǎo)致的地震波至的裝置。
18.如權(quán)利要求17所述的設(shè)備，還包括利用所述去震源特征和去多次反射算子處理所述采集到的地震數(shù)據(jù)的上行成分的裝置。
19.如權(quán)利要求16、17或18所述的設(shè)備，還包括為所述震源選定一個所希望的地震源特征的裝置。
20.如權(quán)利要求17所述的設(shè)備，還包括將處理后的地震數(shù)據(jù)分解為P波和/或S波數(shù)據(jù)的裝置。
21.如權(quán)利要求20所述的設(shè)備，還包括在接收器側(cè)將所述地震數(shù)據(jù)分解為P波和/或S波數(shù)據(jù)的裝置。
22.如權(quán)利要求20或21所述的設(shè)備，還包括在震源側(cè)將所述地震數(shù)據(jù)分解為P波和/或S波數(shù)據(jù)的裝置。
23.一種處理多分量地震數(shù)據(jù)的設(shè)備，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的，該設(shè)備包括將所述采集到的地震數(shù)據(jù)分解為上行成分和下行成分的裝置；從所述采集到的地震數(shù)據(jù)的初始下行成分和環(huán)繞接收器的介質(zhì)的特性來計算一個去震源特征算子的裝置。
24.一種處理多分量地震數(shù)據(jù)的設(shè)備，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的，該設(shè)備包括將所述采集到的地震數(shù)據(jù)分解為PP成分、PS成分、SP成分和SS成分的裝置；從所述采集到的地震數(shù)據(jù)的所述成分和環(huán)繞接收器的介質(zhì)的特性計算一個去震源特征和去多次反射算子的裝置。
25.如權(quán)利要求9或10所述的方法，還包括下述步驟利用所述去震源特征和去多次反射算子處理所述采集到的地震數(shù)據(jù)，以衰減或者消除由多次反射導(dǎo)致的地震波至。
26.一種為多分量地震數(shù)據(jù)計算去震源特征和去多次反射算子的設(shè)備，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的，該設(shè)備包括將所述采集到的地震數(shù)據(jù)分解為上行成分和下行成分的裝置；從所述采集到的地震數(shù)據(jù)的初始下行成分和環(huán)繞接收器的介質(zhì)的特性來計算一個去震源特征和去多次反射算子的裝置。
27.一種為多分量地震數(shù)據(jù)計算去震源特征算子的設(shè)備，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的，該設(shè)備包括將所述采集到的地震數(shù)據(jù)分解為上行成分和下行成分的裝置；從所述采集到的地震數(shù)據(jù)的初始下行成分和環(huán)繞接收器的介質(zhì)的特性來計算一個去震源特征算子的裝置。
28.一種為多分量地震數(shù)據(jù)計算去震源特征和去多次反射算子的設(shè)備，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的，該設(shè)備包括將所述采集到的地震數(shù)據(jù)分解為PP成分、PS成分、SP成分和SS成分的裝置；從所述采集到的地震數(shù)據(jù)的所述成分和環(huán)繞接收器的介質(zhì)的特性計算一個去震源特征和去多次反射算子的裝置。
29.如權(quán)利要求16到28之一所述的設(shè)備，包括一個可編程數(shù)據(jù)處理器。
30.一種包含用于權(quán)利要求29所述設(shè)備的數(shù)據(jù)處理器的程序的存儲介質(zhì)。
全文摘要
本申請公開了一種處理多分量地震數(shù)據(jù)的方法，所述數(shù)據(jù)是通過在震源位置發(fā)射多分量地震能量，并在位于比震源位置更深的位置的多分量地震接收器處采集地震數(shù)據(jù)而得到的。將所述地震數(shù)據(jù)分解為上行成分和下行成分(52)。從所述的地震數(shù)據(jù)的所述下行成分以及環(huán)繞該接收器的介質(zhì)的特性計算一個去震源特征和去多次反射算子(53)。
文檔編號G01W1/00GK1575424SQ02820902
公開日2005年2月2日 申請日期2002年10月28日 優(yōu)先權(quán)日2001年10月26日
發(fā)明者拉斯·艾門德森, 伊基爾·霍爾維克, 喬翰·O·A·羅伯特森 申請人:維斯特恩格科地震控股有限公司