專利名稱：處理地震數(shù)據(jù)的方法和裝置的制作方法
技術(shù)領(lǐng)域：
本發(fā)明涉及處理地震數(shù)據(jù)的方法。其特別涉及包括地震能量傳播的兩個(gè)對應(yīng)模式的地震數(shù)據(jù)的處理方法，例如其中一個(gè)模式由另一模式在地層邊界的部分模式轉(zhuǎn)換而產(chǎn)生。本發(fā)明的方法提供由與地下層具有不同地震屬性的地表面層或接近地表面的層產(chǎn)生的靜態(tài)偏移的估算。本發(fā)明也涉及地震勘測的方法，其包括按照上述方式處理所得到的地震數(shù)據(jù)。本發(fā)明進(jìn)一步涉及處理地震數(shù)據(jù)的設(shè)備。
背景技術(shù)：
圖1是地震勘測裝置的示意圖。在該圖中的勘測裝置是海上勘測裝置，其中地震能量由從拖引船2懸掛于海面下的地震源1發(fā)出。當(dāng)?shù)卣鹪?被開動，就向下發(fā)出地震能量，進(jìn)入地球內(nèi)部并被地貌3反射，該地貌用作地震能量的部分反射體。該反射的地震能量向上通過地球內(nèi)部進(jìn)入海中，并被設(shè)置在海床上的地震接收器4的陣列檢測到。圖1中所示的地震能量路徑通常稱為“主路徑”。例如，通過確定地震能量沿主路徑從地震源1到接收器4的傳播時(shí)間，從而得到關(guān)于地球內(nèi)部的信息。
實(shí)際上在海床和在地球中的不同地質(zhì)組成層之間的邊界處將會出現(xiàn)折射，但為了清楚，在圖1中省去這一點(diǎn)。圖1中只示出了地震能量的一個(gè)反射體，但實(shí)際上地球內(nèi)部中多個(gè)地貌會用作地震能量的部分反射體。而且圖1中只示出了一個(gè)源1，但實(shí)際上海上地震勘測裝置可以具有一個(gè)源陣列。
處理地震數(shù)據(jù)中遇到的一個(gè)問題是地球的地質(zhì)構(gòu)造不統(tǒng)一。實(shí)際上接近地表面(此后稱為“近表面”)的區(qū)域5的屬性通常顯著不同于地下地質(zhì)構(gòu)造6(此后稱為“基層”)的屬性。例如，如果近表面不如基層堅(jiān)固就會出現(xiàn)這種情況。特別的，地震能量在近表面5中的速度會明顯比在基層6中慢，并且因此近表面有時(shí)被稱為“低速層”(LVL)。相比于當(dāng)近表面和基層具有相同的地震屬性時(shí)所記錄的傳播時(shí)間，這種速度上的差別將會導(dǎo)致地震能量從源1傳播到接收器4的時(shí)間偏移，并且該傳播時(shí)間的偏移通常稱為“靜態(tài)偏移”，或簡稱為“靜態(tài)”。
近表面也能夠以其它方式影響入射到接收器上的地震能量。例如，近表面對傳播通過它的地震能量可以具有聚焦或散焦效果，這就可以影響接收器接收到的地震能量的幅度。近表面也可以對傳播通過它的地震能量的幅度進(jìn)行衰減。而且近表面中的塑性衰減和聚焦/散焦效果也可以使地震能量小波的波形變形。
雖然圖1中所示近表面5與海上地震勘測相關(guān)，但是陸地上出現(xiàn)有相同的近表面區(qū)域，從而在基于陸地的地震勘測中也會出現(xiàn)靜態(tài)偏移。
圖1中所示近表面5延伸到地表(海床)，但是并不必須延伸到該表面，并且該近表面5上面可以還有另一個(gè)層面(未示出)。而且，盡管圖1中示出了近表面5與基層6之間定義完好的邊界，但實(shí)際上在基層6與近表面5之間可以是漸進(jìn)的過渡，而不是明確的邊界。
近表面5所產(chǎn)生的靜態(tài)偏移取決于近表面的厚度和地震能量傳播通過該近表面的速度。近表面5的厚度和通過近表面的傳播速度通常出現(xiàn)橫向變化，從而在一個(gè)位置的地震接收器所觀測到的靜態(tài)偏移很可能不同于在另一個(gè)位置的地震接收器所觀測到的靜態(tài)偏移。達(dá)到第一近似值，一個(gè)接收器記錄的全部數(shù)據(jù)集將相對于另一接收器記錄的數(shù)據(jù)要超前或延遲一個(gè)靜態(tài)時(shí)間偏移。
當(dāng)處理地震數(shù)據(jù)的時(shí)候最好要考慮靜態(tài)偏移。除非從地震數(shù)據(jù)中刪除了這些靜態(tài)偏移，否則將存在這樣的不確定性從較深地層的地震事件的到達(dá)時(shí)間變化是由于深度變化，或這些較深地層的橫向位置的原因，或者只是由于近表面5中的傳播效果而產(chǎn)生的。
從圖1中可以看出，地震能量傳播通過近表面5兩次，一次是它從地震源1到反射體3，再一次是它從反射體3到接收器4。兩次橫穿近表面5都會產(chǎn)生靜態(tài)偏移，向下傳播的地震能量通過近表面5時(shí)所產(chǎn)生的偏移就是源側(cè)靜態(tài)偏移，而能量向上傳播通過近表面5所產(chǎn)生的偏移就是接收器側(cè)靜態(tài)偏移。
本發(fā)明的方法使用s模式事件和對應(yīng)的p模式事件來獲得低速層的信息。假定s模式事件和對應(yīng)的p模式事件之間的差別，例如事件的傳播時(shí)間、幅度和/或波形的差別全部或主要在近表面5中產(chǎn)生。
本發(fā)明的一個(gè)實(shí)施例使用s模式事件和對應(yīng)的p模式事件，其中一個(gè)模式是由另一模式的部分模式轉(zhuǎn)換產(chǎn)生的。例如，如果基層6和近表面5之間的物理屬性對比足夠大，當(dāng)?shù)卣鹉芰總鞑ネㄟ^基層6和近表面5之間的邊界時(shí)，p模式的能量傳播和s模式的能量傳播之間發(fā)生顯著的轉(zhuǎn)換，或相反。在下面參照圖2對這種模式轉(zhuǎn)換進(jìn)行更全面的描述。對于近表面和基層之間的邊界上的接收器記錄到的每一p模式事件，這種模式轉(zhuǎn)換的結(jié)果就是得到對應(yīng)的s模式事件。地震源開動后，p模式事件和對應(yīng)的s模式事件將會在不同的時(shí)間出現(xiàn)在接收器4處，因?yàn)樵诮砻嬖趐模式和s模式傳播速度彼此并不相等。然而由于在近表面的低層邊界出現(xiàn)模式轉(zhuǎn)換，p模式事件和對應(yīng)的s模式事件之間的時(shí)間延遲必須同近表面中兩個(gè)模式不同速度的結(jié)果一起產(chǎn)生。p模式事件和對應(yīng)的s模式事件之間的時(shí)間延遲將不會顯著受到基層屬性的影響。于是通過識別地震數(shù)據(jù)中p模式事件和對應(yīng)的s模式事件，和確定p模式事件和對應(yīng)的s模式事件之間的時(shí)間延遲，可以直接確定p和s模式之間穿過近表面的傳播時(shí)間差。
對于任何表示地震數(shù)據(jù)方面的參數(shù)，例如壓力或質(zhì)點(diǎn)運(yùn)動的分量(諸如質(zhì)點(diǎn)位移或質(zhì)點(diǎn)速度分量)，p和s模式之一事件通常會比另一模式事件表現(xiàn)出更強(qiáng)烈。于是通過p模式會主要的出現(xiàn)在一個(gè)參數(shù)，諸如質(zhì)點(diǎn)運(yùn)動的垂直分量，相應(yīng)的被模式轉(zhuǎn)換的s模式將會主要出現(xiàn)在另一參數(shù)中，諸如質(zhì)點(diǎn)運(yùn)動的水平分量。
當(dāng)向下傳播的波經(jīng)過近表面5和基層6之間的界面處的反射時(shí)，也可以出現(xiàn)模式轉(zhuǎn)換。本發(fā)明可以用于這種情況，由于模式轉(zhuǎn)換在反射處出現(xiàn)的地方，反射的信號將包含p分量和對應(yīng)的s分量，這兩個(gè)分量穿過近表面5會具有不同的傳播時(shí)間。
當(dāng)向下傳播的波經(jīng)過近表面5和基層6之間的界面處的臨界折射時(shí)，也可以出現(xiàn)模式轉(zhuǎn)換，以產(chǎn)生沿近表面5和基層6之間的界面?zhèn)鞑サ牡卣鸩āＱ卦摻缑鎮(zhèn)鞑サ牟▽⒃诮砻婕頿模式和s模式，本發(fā)明可以用于這些p和s模式。
共同未決的英國專利申請No.0101936.3和No.0110573.3公開了確定p波形的靜態(tài)偏移和s波形的靜態(tài)偏移之間的差別的方法，其中s波形通過在近表面和基層之間的邊界處的模式轉(zhuǎn)換產(chǎn)生。在地震數(shù)據(jù)中識別成對的p事件和相應(yīng)s事件，就發(fā)現(xiàn)了p事件和相應(yīng)的s事件之間的傳播時(shí)間差。對于特定的接收器，該傳播時(shí)間差是s波形靜態(tài)偏移和p波形靜態(tài)偏移之間的差別。
現(xiàn)在許多地震勘測使用多分量接收器，記錄入射到接收器上的地震能量的兩個(gè)或多個(gè)分量。例如3分量(3-C)地震接收器可以在接收器記錄質(zhì)點(diǎn)運(yùn)動的x-、y-和z-分量，在海上勘測中，4分量(4-C)地震接收器除了可以記錄質(zhì)點(diǎn)運(yùn)動的x-、y-和z-分量之外，還可以記錄水柱的壓力(它是標(biāo)量)。在處理多分量地震數(shù)據(jù)中進(jìn)一步有一個(gè)問題是，接收器可以比其它分量更精確的記錄地震數(shù)據(jù)的某些分量。這個(gè)問題就是矢量失真。產(chǎn)生矢量失真的一個(gè)原因是接收器和地表面之間缺少足夠的耦合，從而在接收器陣列中，陣列的接收器的矢量失真效果可以彼此不同。這種矢量失真將有損于通過處理地震數(shù)據(jù)得到的結(jié)果。

發(fā)明內(nèi)容
本發(fā)明的第一方面是提供處理地震數(shù)據(jù)的方法，該地震數(shù)據(jù)包括地震能量對應(yīng)的第一和第二模式，該方法包括a)處理表示地震數(shù)據(jù)第一參數(shù)的軌跡，該地震數(shù)據(jù)在至少包括對應(yīng)于第一模式的事件的第一接收器處獲得，和處理表示地震數(shù)據(jù)第二參數(shù)的軌跡，該地震數(shù)據(jù)在至少包括對應(yīng)于第二模式的事件的第二接收器處獲得；
b)在處理的數(shù)據(jù)中識別對應(yīng)于部分模式轉(zhuǎn)換的事件；和c)從處理的數(shù)據(jù)中對應(yīng)于部分模式轉(zhuǎn)換的事件的幅度和/或波形獲得信息。
該方法可以用來從對應(yīng)于部分模式轉(zhuǎn)換事件的幅度和/或波形獲得近表面的屬性信息。另外或可替換的，它可以用來從一個(gè)到另一個(gè)接收器獲得耦合強(qiáng)度中的變化信息。
在優(yōu)選實(shí)施例中，本發(fā)明被應(yīng)用到數(shù)據(jù)中，其中通過第一模式在海床層的邊界面的部分模式轉(zhuǎn)換產(chǎn)生第二模式。然而本發(fā)明并不限于此，通?？梢詰?yīng)用到任何事件對，其中一個(gè)處于一個(gè)模式，一個(gè)處于另一模式，由此該兩個(gè)事件的幅度和波形之間的差別主要會在低速層中產(chǎn)生。
優(yōu)選的選擇用于本發(fā)明的地震數(shù)據(jù)參數(shù)，使得p模式將主要的出現(xiàn)在一個(gè)參數(shù)中，而其相對應(yīng)的模式轉(zhuǎn)換的s模式將主要的出現(xiàn)在另一參數(shù)中。該參數(shù)可以是在接收器直接得到的參數(shù)，諸如質(zhì)點(diǎn)運(yùn)動的垂直分量和質(zhì)點(diǎn)運(yùn)動的水平分量?？商鎿Q的，該參考可以是通過處理在接收器得到的地震數(shù)據(jù)而獲得的，例如通過將所得到的地震數(shù)據(jù)分解成它的上行和下行組分或分解成它的p和s組分。
本發(fā)明的方法可以用來處理先前存在的地震數(shù)據(jù)?？商鎿Q的，它可以被包括在用于在得到數(shù)據(jù)的時(shí)候或之后處理數(shù)據(jù)的地震勘測方法中。
本發(fā)明的第二方面提供一種地震勘測的方法，包括步驟將以第一模式傳播的地震能量導(dǎo)向到海床層的邊界面，使得在該邊界面出現(xiàn)地震能量轉(zhuǎn)換成第二模式的部分模式轉(zhuǎn)換；在一個(gè)或多個(gè)接收器得到包括第一和第二模式地震能量的地震數(shù)據(jù)；和根據(jù)如上所定義的方法處理該地震數(shù)據(jù)。
本發(fā)明的第三方面提供一種處理包括第一和第二模式地震能量的地震數(shù)據(jù)的設(shè)備，該設(shè)備包括a)裝置，用于處理表示在至少包括對應(yīng)于第一模式的事件的第一接收器處獲得的地震數(shù)據(jù)第一參數(shù)的軌跡，以及處理表示在至少包括對應(yīng)于第二模式的事件的第二接收器處獲得的地震數(shù)據(jù)第二參數(shù)的軌跡；b)裝置，用于在處理的數(shù)據(jù)中識別對應(yīng)于部分模式轉(zhuǎn)換的事件；和c)裝置，用于從處理的數(shù)據(jù)中對應(yīng)于部分模式轉(zhuǎn)換的事件的幅度和/或波形獲得信息。
該裝置可以包括可編程數(shù)據(jù)處理器。
本發(fā)明的第四方面提供一種存儲介質(zhì)，其包含有用于如上所定義的裝置的數(shù)據(jù)處理器的程序。


現(xiàn)在參照附圖，通過說明性范例的方式將詳細(xì)描述本發(fā)明的優(yōu)選實(shí)施例，其中圖1是海上地震勘測的示意圖；圖2是上行p波在基層和低速層之間的邊界處模式轉(zhuǎn)換的示意說明；圖3所示是在圖2中模式轉(zhuǎn)換出現(xiàn)的地方地震數(shù)據(jù)軌跡的兩個(gè)示意范例；圖4所示是將圖3中的兩個(gè)地震數(shù)據(jù)軌跡彼此相互關(guān)聯(lián)的結(jié)果；圖5所示是英國專利申請No.0101936.3和No.0110573.3的相互關(guān)聯(lián)方法的結(jié)果；圖6(a)是從圖5的結(jié)果中選擇的軌跡的放大視圖；圖6(b)是對圖6(a)的軌跡歸一化的結(jié)果；圖7是從圖6(b)中的軌跡中選擇的事件的幅度變化；圖8是根據(jù)本發(fā)明一個(gè)實(shí)施例的方框圖；和圖9是根據(jù)本發(fā)明一個(gè)實(shí)施例的設(shè)備的示意方框圖。
具體實(shí)施例方式
現(xiàn)在將參照實(shí)施例對本發(fā)明的原理進(jìn)行說明，在該實(shí)施例中通過模式轉(zhuǎn)換產(chǎn)生對應(yīng)的事件對。圖2說明了上行傳播的p波通過基層6和近表面5之間的界面時(shí)出現(xiàn)的模式轉(zhuǎn)換。上行p波部分的以p波7′發(fā)射，但也經(jīng)過部分模式轉(zhuǎn)換成s波7″。假定通過地球內(nèi)部兩層之間深界面的反射、或穿過更深地層的折射產(chǎn)生上行p波。在基層6與近表面5之間的界面也會出現(xiàn)折射，發(fā)射的p波7′的折射角與轉(zhuǎn)換的s波7″的折射角不同。
發(fā)射的p波7′和轉(zhuǎn)換的s波7″都會在接收器陣列產(chǎn)生信號。如圖2所示，發(fā)射的p波7′和轉(zhuǎn)換的s波7″以不同的方向傳播，因此會到達(dá)不同的接收器位置，而不是相同的接收器。然而來自鄰近轉(zhuǎn)換點(diǎn)的類似的被轉(zhuǎn)換的s波會到與達(dá)圖2中所示的被發(fā)射的p波7′相同的接收器。將要這樣假定對于特定接收器所接收的被發(fā)射的p波和被轉(zhuǎn)換的s波，被轉(zhuǎn)換的s波的轉(zhuǎn)換點(diǎn)(近表面基點(diǎn)附近)與p波發(fā)射通過的近表面基點(diǎn)中的點(diǎn)之間的差別并不會顯著影響本發(fā)明的方法所得到的靜態(tài)偏移估算。接收器獲得的數(shù)據(jù)于是將包含由于被發(fā)射的p波7′的到達(dá)而產(chǎn)生的事件，和由于對應(yīng)的被轉(zhuǎn)換的s波7″的到達(dá)而產(chǎn)生的事件。
被發(fā)射的p波7′將主要被記錄在接收器所記錄的質(zhì)點(diǎn)速度的垂直分量上，而被轉(zhuǎn)換的s波7″將主要被記錄在接收器4所記錄的質(zhì)點(diǎn)速度的徑向分量上(徑向方向是投影到海床上的源接收器方向，該方向?qū)⒈欢x為x方向)。這是因?yàn)樵趯?shí)際的地震勘測裝置中，從地球內(nèi)部中向上傳播的波形由于地球內(nèi)部的典型結(jié)構(gòu)和目標(biāo)反射體的典型深度，容易產(chǎn)生30°的角或小于直角。剪切波主要記錄在徑向分量上，因?yàn)橄啾扔趐波沿傳播方向的質(zhì)點(diǎn)運(yùn)動，剪切波的質(zhì)點(diǎn)運(yùn)動垂至于傳播方向。
圖3是表示記錄在接收器4的地震數(shù)據(jù)兩個(gè)參數(shù)軌跡的示意說明，其中該接收器4是如圖2所示的部分模式轉(zhuǎn)換出現(xiàn)的地方。圖3中所示的兩個(gè)參數(shù)是接收器所測得的質(zhì)點(diǎn)運(yùn)動的徑向分量(x分量)，和接收器所測得質(zhì)點(diǎn)運(yùn)動的垂直分量(z分量)。由于s波7″在從上行p波7的轉(zhuǎn)換時(shí)產(chǎn)生，如果基層和近表面之間的界面大致是平坦的，那么該記錄的s波將包含與p到達(dá)信號相似的小波。s波到達(dá)不會與p波到達(dá)同時(shí)出現(xiàn)，而會(相比于p波到達(dá))在一個(gè)時(shí)間延遲dt后到達(dá)，該時(shí)間延遲由近表面5中p波速度和s波速度之間的差產(chǎn)生。該時(shí)間延遲dt就是p波靜態(tài)偏移和s波靜態(tài)偏移之間的差。
在共同未決的英國專利申請No.0101936.3和No.0110573.3的方法中，第一參數(shù)與第二參數(shù)相互關(guān)聯(lián)，或者從第二參數(shù)中卷積展開。通過相互關(guān)聯(lián)記錄圖3中的兩個(gè)軌跡得到的結(jié)果相互關(guān)聯(lián)記錄圖如圖4所示。在該相互關(guān)聯(lián)記錄圖中零時(shí)間偏移有個(gè)峰值，這是因?yàn)閜波7′到達(dá)在測得質(zhì)點(diǎn)運(yùn)動的x分量具有非零的幅度(雖然它主要出現(xiàn)在測得質(zhì)點(diǎn)運(yùn)動的z分量中)。任何在測得質(zhì)點(diǎn)運(yùn)動的垂直分量中出現(xiàn)的s波能量也會產(chǎn)生相互關(guān)聯(lián)記錄圖中零時(shí)間偏移的峰值。該相互關(guān)聯(lián)記錄圖在時(shí)間dt也有個(gè)峰值，該峰值從測得質(zhì)點(diǎn)運(yùn)動的垂直分量中的p波到達(dá)峰值、和測得質(zhì)點(diǎn)運(yùn)動x分量中對應(yīng)的模式轉(zhuǎn)換的s波到達(dá)峰值產(chǎn)生。于是可以從相互關(guān)聯(lián)記錄圖中確定該時(shí)間偏移dt，從而確定p波靜態(tài)和s波靜態(tài)之間的差別。
圖5是共同未決的英國專利申請No.0101936.3和No.0110573.3方法結(jié)果的進(jìn)一步說明。這些結(jié)果通過使用包含大約130個(gè)接收器的接收器線而獲得。對于每一接收器，圖5示出了從記錄在該接收器的質(zhì)點(diǎn)運(yùn)動垂直分量和記錄在接收器的質(zhì)點(diǎn)運(yùn)動x分量(徑向分量)獲得的相互關(guān)聯(lián)記錄圖。對于每一接收器p波靜態(tài)偏移和s波靜態(tài)偏移之間的差別，如由英國專利申請No.0101936.3和No.0110573.3的方法中所確定的，在圖5中用A表示。使用常規(guī)計(jì)數(shù)得到的p波靜態(tài)偏移和s波靜態(tài)偏移之間的差別也示出在圖5中，用B表示。
英國專利申請No.0101936.3和No.0110573.3的方法提供一種估算p波靜態(tài)偏移和s波靜態(tài)偏移之間差別的改進(jìn)方法。然而該方法只是使用了時(shí)間信息，并只能確定p模式到達(dá)和對應(yīng)的s模式到達(dá)之間的時(shí)間延遲。從英國專利申請No.0101936.3和No.0110573.3的方法中不可能確定該時(shí)間延遲是由于近表面的地質(zhì)效果，或者是由于其它因素，諸如接收器和地球之間的耦合、或其它源的矢量失真。
通過考慮相關(guān)數(shù)據(jù)中事件的幅度和/或波形變化，本發(fā)明開始從p模式事件和對應(yīng)的s模式事件之間的時(shí)間延遲得出進(jìn)一步的信息。
現(xiàn)在將參照圖6(a)至7描述本發(fā)明的一個(gè)實(shí)施例。
圖6(a)是從圖5中選擇的軌跡的放大視圖。這些軌跡所示為p模式和s模式的相關(guān)，為了方便稱其為“接收器函數(shù)”。圖6(a)示出了從圖5中的接收器函數(shù)，其為每個(gè)第五接收器，p到達(dá)和s到達(dá)之間的時(shí)間延遲為正值(正時(shí)間延遲表示p模式事件首先到達(dá)接收器，s模式事件隨后，如圖3所示)。
圖6(b)所示為歸一化圖6(a)中的接收器函數(shù)的效果，其分別對每一接收器函數(shù)中峰值的幅度以零時(shí)間延遲進(jìn)行歸一化。也就是，通過歸一化接收器函數(shù)，使得對于圖6(a)中的每一接收器函數(shù)，峰值在零時(shí)間延遲時(shí)具有相同的幅度，從而就從圖6(a)中對應(yīng)的接收器函數(shù)獲得圖6(b)中的每一接收器函數(shù)。如上所示，圖4和5中接收器函數(shù)中零時(shí)間延遲的峰值通過p波到達(dá)泄漏成為測得質(zhì)點(diǎn)運(yùn)動的x分量，通過s波到達(dá)泄漏成為測得質(zhì)點(diǎn)運(yùn)動的y分量。也就是對于特定接收器，接收器函數(shù)中零時(shí)間延遲的峰值幅度表示p波到達(dá)在記錄在該接收器上的質(zhì)點(diǎn)運(yùn)動的x分量上的投影，和s波到達(dá)在質(zhì)點(diǎn)運(yùn)動的y分量上的投影。
通過檢查圖6(a)和6(b)可能會發(fā)現(xiàn)，歸一化接收器函數(shù)對該接收器函數(shù)中事件的幅度并沒有明顯的效果。這表明對于每一接收器，零時(shí)間延遲的峰值幅度大致相等。這又表明對于每一接收器，p波到達(dá)在所得到的質(zhì)點(diǎn)運(yùn)動的x分量上的投影大致恒定。它也表明接收器彼此之間，對于p波耦合強(qiáng)度并沒有明顯的變化。
可以看出圖6(b)中的接收器函數(shù)包含時(shí)間延遲0.3s左右的峰值，雖然接收器之間該峰值出現(xiàn)的準(zhǔn)確時(shí)間延遲稍微有變化。該峰值表示p波到達(dá)和對應(yīng)的模式轉(zhuǎn)換的s波到達(dá)之間的時(shí)間延遲。接收器彼此之間該峰值波形和/或幅度的變化可以提供有關(guān)于s波在近表面中衰減和聚焦效果的信息。
為了方便，由于p波與相應(yīng)的模式轉(zhuǎn)換的s波之間的時(shí)間延遲而引起的接收器函數(shù)中的峰值被稱為“靜態(tài)峰值”。該峰值的幅度和相位也會受到透射系數(shù)變化的影響，和在近表面與基層之間界面出現(xiàn)的模式轉(zhuǎn)換的程度的影響。例如作為材料屬性變化和的沿界面地勢的結(jié)果，不同界面之間透射系數(shù)和轉(zhuǎn)換效率會不同。
對于圖6(b)中所示的每一接收器函數(shù)，圖7示出了該“靜態(tài)峰值”(時(shí)間延遲0.3秒左右的峰值)的幅度，相對于接收器編號設(shè)置。圖7示出了圖5所示每一接收器函數(shù)的結(jié)果，因此包括比圖6(b)中所示更多接收器的數(shù)據(jù)。
可以看出，圖7示出了接收器彼此之間的接收器函數(shù)峰值歸一化的幅度具有相當(dāng)多的變化。圖7中短波長幅度的變化(那是出現(xiàn)在小數(shù)量接收器之間的變化)可能是由于接收器彼此不同的接收器到地球的耦合強(qiáng)度而引起的。圖7中長波長幅度的變化更有可能是由于接收器彼此不同的近表面構(gòu)造而產(chǎn)生的地質(zhì)效果而引起的，諸如由于近表面而產(chǎn)生的衰減或聚焦/散焦的差別。于是，例如使用傅立葉變換分析圖7中的幅度變化可以提供關(guān)于地質(zhì)效果和耦合效果的信息。
本發(fā)明也有可能補(bǔ)償接收器位置彼此不同的接收器與地球耦合強(qiáng)度的變化。由于耦合強(qiáng)度的變化可以使s波到達(dá)相對于對應(yīng)的p波到達(dá)的幅度和相位變形，修正接收器耦合強(qiáng)度的變化是重要的。由于s波到達(dá)主要記錄在接收器得到的質(zhì)點(diǎn)運(yùn)動的x分量中，耦合效率的變化可以使質(zhì)點(diǎn)運(yùn)動記錄的x分量變形，接收器彼此之間的變形量會不同。這在應(yīng)用中特別是問題，例如AVO(幅度-偏移)應(yīng)用，它需要質(zhì)點(diǎn)運(yùn)動x分量的記錄在接收器分布上基本上一致，并且可以與另一個(gè)直接比較。
圖6(b)或圖7的接收器函數(shù)中歸一化的“靜態(tài)峰值”的幅度表示模式轉(zhuǎn)換的s波到達(dá)的強(qiáng)度，如上所示其由近表面產(chǎn)生。“靜態(tài)峰值”的幅度也表示接收器的耦合強(qiáng)度。對于相同的接收器，相對于圖6(b)歸一化的接收器函數(shù)中“靜態(tài)峰值”的幅度，通過歸一化接收器得到的質(zhì)點(diǎn)運(yùn)動x分量原始軌跡的幅度，本發(fā)明有可能補(bǔ)償接收器彼此不同的近表面效果和地球耦合強(qiáng)度的變化。也就是對于每一接收器，歸一化的接收器函數(shù)中“靜態(tài)峰值”的幅度被確定，然后該幅度被用來歸一化接收器的質(zhì)點(diǎn)運(yùn)動x分量的原始軌跡。對于每一接收器重復(fù)該過程——對于該接收器，記錄在接收器的質(zhì)點(diǎn)速度x分量的軌跡被歸一化的接收器函數(shù)中“靜態(tài)峰值”的幅度歸一化。該歸一化處理結(jié)果用于提供歸一化的x分量軌跡，該軌跡在接收器分布上基本上一致，并可以與另一個(gè)直接比較。該處理有效的使得記錄在接收器的質(zhì)點(diǎn)運(yùn)動的x分量相對于記錄在該接收器的質(zhì)點(diǎn)運(yùn)動的垂直分量(其應(yīng)該主要包含p波事件，從而可以相對的不受近表面或耦合變化的影響)被校準(zhǔn)。
另外，如果接收器質(zhì)點(diǎn)運(yùn)動的y分量也像質(zhì)點(diǎn)運(yùn)動的x分量一樣被獲得，將有可能相對于質(zhì)點(diǎn)運(yùn)動的垂直分量校準(zhǔn)質(zhì)點(diǎn)運(yùn)動的y分量。
如果完全3C或4C地震數(shù)據(jù)是有效的，則有可能計(jì)算兩個(gè)相關(guān)圖——質(zhì)點(diǎn)運(yùn)動的x分量相對于z分量的第一相互關(guān)聯(lián)記錄圖，和質(zhì)點(diǎn)運(yùn)動的y分量相對于z分量的第二相互關(guān)聯(lián)記錄圖。每一個(gè)這些相關(guān)圖然后可以被歸一化，使得相關(guān)圖中每一軌跡中峰值在零時(shí)間延遲具有相同的幅度。這兩個(gè)歸一化的相互關(guān)聯(lián)記錄圖然后可以直接相互比較——因?yàn)樵撎幚硐鄬τ谫|(zhì)點(diǎn)運(yùn)動的z分量校準(zhǔn)了質(zhì)點(diǎn)運(yùn)動的x分量，也相對于質(zhì)點(diǎn)運(yùn)動的z分量校準(zhǔn)了質(zhì)點(diǎn)運(yùn)動的y分量，因此等同于相對于質(zhì)點(diǎn)運(yùn)動的y分量校準(zhǔn)質(zhì)點(diǎn)運(yùn)動的x分量。
圖7只包括關(guān)于圖6中軌跡的幅度變化的信息。然而可以看到，圖6(b)的相互關(guān)聯(lián)記錄圖軌跡中“靜態(tài)峰值”的波形隨接收器位置的不同而變化。靜態(tài)峰值波形中的變化也可以提供不同接收器位置的關(guān)于近表面和/或接收器與地球之間耦合強(qiáng)度變化的信息。
從概念上說有可能構(gòu)建顯示出在圖6(b)的歸一化的接收器函數(shù)中“靜態(tài)峰值”波形的某些測量中變化的圖，例如在波形最大高度一半時(shí)的全寬度。這種圖示出這一峰值波形如何隨著接收器位置的不同而變化，并且“靜態(tài)峰值”波形中的這些變化也可以提供關(guān)于近表面地質(zhì)效果和關(guān)于接收器與地球之間耦合強(qiáng)度變化的信息。
可替換的，有可能得出將圖6(b)中一個(gè)接收器函數(shù)中的“靜態(tài)峰值”波形轉(zhuǎn)換為不變形波形的濾波器，諸如軌跡中零時(shí)間延遲的峰值波形。該濾波器然后可以用于原始地震數(shù)據(jù)，例如用于質(zhì)點(diǎn)運(yùn)動x分量的原始軌跡。對于每一接收器可以重復(fù)該處理——每一接收器x分量的軌跡被濾波器轉(zhuǎn)換，該濾波器從該接收器的接收器函數(shù)中的“靜態(tài)峰值”波形得出。該處理應(yīng)該提供在接收器分布上基本上一致的，并可以與另一個(gè)直接比較的歸一化的x分量軌跡。
在上述的實(shí)施例中通過相互關(guān)聯(lián)計(jì)算接收器函數(shù)。然而本發(fā)明并不限于此，任何合適的相關(guān)技術(shù)和算法可以被用來確定該接收器函數(shù)。例如可以通過分量之間去卷積來得到接收器函數(shù)，并且靜態(tài)峰值波形中的幅度和波形信息然后仍然可以按照如上所述相同的方式被使用。實(shí)際上通過去卷積得到的接收器函數(shù)可能會更好，因?yàn)槿ゾ矸e有效的刪除了p波的特定源信息(它存在于入射近表面和/或基層界面)，使得在接收器函數(shù)中只留下由于近表面構(gòu)造和/或矢量失真而產(chǎn)生的效果。
另一個(gè)用于得到接收器函數(shù)的合適技術(shù)的示例是雙相干時(shí)間延遲估算方法，如L.Ikelle在“Geophysics”，Vol 62，p1947(1997)中所述。
在上述本發(fā)明的描述中，接收器所得到的質(zhì)點(diǎn)速度分量在計(jì)算接收器函數(shù)中被用作質(zhì)點(diǎn)運(yùn)動的分量。然而可以替換的使用接收器的質(zhì)點(diǎn)位移分量、或接收器的質(zhì)點(diǎn)加速度分量、或質(zhì)點(diǎn)速度的更高階導(dǎo)數(shù)。
實(shí)際上本發(fā)明并不限于使用質(zhì)點(diǎn)運(yùn)動的分量，所得到的地震數(shù)據(jù)的認(rèn)識適當(dāng)?shù)膮?shù)都可以用于獲得接收器函數(shù)。例如該地震數(shù)據(jù)可以被分解成它的上行和下行分量，或分解成p和s分量，這些可以被用于獲得接收器函數(shù)。
本發(fā)明也可以用于p和s波“回響事件”。在這種事件部分模式轉(zhuǎn)換中，例如根據(jù)傳播通過近表面5與基層6之間的界面，或在其上的反射。然而在回響事件中，p和s波并不直接從接收器出現(xiàn)模式轉(zhuǎn)換的點(diǎn)通過。然而在回響事件中，p和s波在近表面的邊界經(jīng)過一個(gè)或多個(gè)反射，從而在入射到接收器之前多次通過近表面5。
地震數(shù)據(jù)也可以包含由于不是出現(xiàn)在近表面邊界、而是近表面之中的部分模式轉(zhuǎn)換而產(chǎn)生的事件。這可能是由于例如近表面中的成層效果產(chǎn)生的。
在如上所述的實(shí)施例中，p和s波事件之一由部分模式轉(zhuǎn)換產(chǎn)生。然而如上所示本發(fā)明并不限于此，并可以用于任何對應(yīng)的p和s事件對，其中該兩個(gè)事件的幅度和/或波形的差別主要由于近表面5而產(chǎn)生。例如，在使用既發(fā)射p波又發(fā)射s波的地震勘測中，所得到的數(shù)據(jù)將包含由于在近表面5低層上的特定點(diǎn)p波的反射產(chǎn)生的事件，也包含由于在近表面5低層上的該點(diǎn)s波的反射產(chǎn)生的對應(yīng)事件。本發(fā)明可以用于這些事件。應(yīng)該注意到，由于水柱并不支持s波傳播，發(fā)射s和p模式的源需要設(shè)置在海床上，地球表面上或鉆孔中。
圖18是本發(fā)明一個(gè)實(shí)施例的流程圖。首先，在步驟1，在勘測位置得到多分量的震數(shù)據(jù)。該地震數(shù)據(jù)應(yīng)該包括該質(zhì)點(diǎn)運(yùn)動的垂直分量和該質(zhì)點(diǎn)運(yùn)動的至少一個(gè)水平分量。設(shè)置在勘測區(qū)域中不同位置的多個(gè)接收器獲得該數(shù)據(jù)。多分量數(shù)據(jù)可以是全3C或4C地震數(shù)據(jù)。
除了將本發(fā)明用于新得到的地震數(shù)據(jù)，本發(fā)明也可以用于先前存在的勘測數(shù)據(jù)。因此步驟1可以被步驟2替代從存儲器中檢索合適的多分量地震數(shù)據(jù)。
在步驟3，選擇相關(guān)于第i個(gè)接收器的地震數(shù)據(jù)。在步驟4，從第i個(gè)接收器獲得的地震數(shù)據(jù)的兩個(gè)參數(shù)計(jì)算第i個(gè)接收器的接收器函數(shù)。在一個(gè)實(shí)施例中，通過將第i個(gè)接收器得到的質(zhì)點(diǎn)運(yùn)動的垂直分量與第i個(gè)接收器得到的質(zhì)點(diǎn)運(yùn)動的水平分量進(jìn)行相互關(guān)聯(lián)獲得該接收器函數(shù)，其中垂直分量作為該得到的地震數(shù)據(jù)的一個(gè)參數(shù)，水平分量作為該得到的地震數(shù)據(jù)的另一參數(shù)。在該實(shí)施例中，假定水平分量是接收器質(zhì)點(diǎn)運(yùn)動的x分量(徑向分量)，記錄在接收器的質(zhì)點(diǎn)運(yùn)動是質(zhì)點(diǎn)速度，但本發(fā)明并不限于此。
步驟4的結(jié)果給出圖5或圖6(a)中所示類型的第i個(gè)接收器的接收器函數(shù)。
在步驟5，第i個(gè)接收器的接收器函數(shù)相對于接收器函數(shù)零時(shí)間延遲峰值的幅度被歸一化。步驟5的結(jié)果是圖6(b)中所示第i個(gè)接收器的歸一化的接收器函數(shù)。
在步驟6，確定該歸一化的接收器函數(shù)中的“靜態(tài)峰值”。這可以通過例如使用自動選擇方法完成，或者可以通過手動選擇“靜態(tài)峰值”的操作者輸出步驟5的結(jié)果用于顯示來完成。
在步驟7，確定該歸一化的接收器函數(shù)中的“靜態(tài)峰值”幅度。在步驟8，“靜態(tài)峰值”幅度被用來歸一化記錄在第i個(gè)接收器的質(zhì)點(diǎn)運(yùn)動水平分量的原始軌跡。
在步驟9，確定來自接收器陣列中的所有接收器的數(shù)據(jù)是否被處理完。如果確定的結(jié)果是“否”，一個(gè)適當(dāng)?shù)挠?jì)數(shù)器在步驟10增1，然后對于第i+1個(gè)接收器重復(fù)步驟3至9。
當(dāng)在步驟9中得到“是”的確定結(jié)果，這表示來自陣列中的所有接收器的數(shù)據(jù)已經(jīng)被處理完畢。該處理的結(jié)果然后在步驟11可以被輸出。數(shù)據(jù)可以任何適當(dāng)?shù)男问捷敵觥缢梢源蛴〉郊埳稀@示到可視顯示單元上、或可以被存儲到用于順序查找的任何適當(dāng)形式?？商鎿Q的，該數(shù)據(jù)可以提交給附加處理，以確定與勘測位置的地質(zhì)構(gòu)造相關(guān)的信息。
在步驟11可以選擇數(shù)據(jù)輸出的形式，以適合用戶的特定需要。輸出的數(shù)據(jù)可以包括步驟4所產(chǎn)生的接收器函數(shù)、步驟5所產(chǎn)生的歸一化的接收器函數(shù)、步驟7中確定的歸一化的接收器函數(shù)中“靜態(tài)峰值”的幅度、或步驟8所確定的質(zhì)點(diǎn)運(yùn)動水平分量歸一化的軌跡中的一些或全部。
本發(fā)明的方法并不限于按照圖8中的順序執(zhí)行這些步驟。例如選擇“靜態(tài)峰值”的步驟6可以在歸一化接收器函數(shù)的步驟5之前進(jìn)行。而且在圖8中所示的方法中，對于一個(gè)接收器順序執(zhí)行步驟3至8，然后對于另一接收器重復(fù)執(zhí)行等等?？商鎿Q的，有可能對于一個(gè)接收器執(zhí)行步驟3和4，然后對于其它接收器在進(jìn)入步驟5之前重復(fù)步驟3和4。相似的，在進(jìn)入選擇“靜態(tài)峰值”的步驟6之前，可以對于每一接收器執(zhí)行歸一化接收器函數(shù)的步驟5。同樣的，在進(jìn)入確定“靜態(tài)峰值”幅度的步驟7之前，可以對于每一接收器執(zhí)行選擇“靜態(tài)峰值”的步驟6，并且在進(jìn)入歸一化質(zhì)點(diǎn)運(yùn)動的水平分量的步驟8之前，可以對于每一接收器執(zhí)行確定“靜態(tài)峰值”幅度的步驟7。
如果可能識別“靜態(tài)峰值”或與各個(gè)未疊層的接收器函數(shù)中p和s波穿過近表面的時(shí)間差相關(guān)的另一事件，圖8中的方法是適合的。然而在某些情況中，不可能識別各個(gè)接收器函數(shù)中的“靜態(tài)峰值”，在這種情況下，需要包括疊層接收器函數(shù)以增強(qiáng)信噪比的步驟。
如果需要對“靜態(tài)峰值”識別的疊層步驟，圖8中的步驟4將被對第i個(gè)接收器集合計(jì)算所有接收器函數(shù)的步驟代替。第i個(gè)接收器集合包含所有記錄在第i個(gè)接收器的軌跡，并且對于這些軌跡中的每一個(gè)，有可能計(jì)算接收器函數(shù)。該步驟之后是疊層所有第i個(gè)接收器集合所計(jì)算得到的接收器函數(shù)的步驟(圖8中未示出)。然后對第i個(gè)接收器集合的疊層的接收器函數(shù)進(jìn)行圖8中的步驟5至7。
然后對于為其它接收器集合所確定的接收器函數(shù)重復(fù)該疊層處理。
圖9根據(jù)本發(fā)明的可編程設(shè)備5的示意方框圖。該系統(tǒng)包括可編程數(shù)據(jù)處理器8和程序存儲器9，例如是只讀存儲器ROM形式的，用于存儲控制數(shù)據(jù)處理器8的程序以執(zhí)行，例如如圖8中所描述的方法。該系統(tǒng)進(jìn)一步包括非易失性讀/寫存儲器10用于存儲，例如任何在缺少供電電源時(shí)必須保留的數(shù)據(jù)。為數(shù)據(jù)處理器提供隨機(jī)訪問存儲器(RAM)11作為“工作”或“中間結(jié)果暫存器”存儲器。提供輸入接口12，例如用于接收命令或數(shù)據(jù)。提供輸出接口13，例如用于顯示與進(jìn)程相關(guān)的信息和該方法的結(jié)果。用戶定義的數(shù)據(jù)可以通過輸入接口12提供，或可選的通過機(jī)器可讀的存儲器14提供。
用于操作系統(tǒng)和用于執(zhí)行此前所描述的方法的程序存儲在程序存儲器9中，其可以是半導(dǎo)體存儲器，例如熟知的ROM類型的。然而該程序可以在任何其它合適的存儲器媒體中，例如磁數(shù)據(jù)載體9a(諸如“軟盤”)或CD-ROM 9b。
權(quán)利要求
1.一種處理包括對應(yīng)于第一和第二模式的地震能量的地震數(shù)據(jù)的方法，該方法包括a)處理表示在第一接收器處獲得的至少包括對應(yīng)于第一模式的事件的地震數(shù)據(jù)的第一參數(shù)的軌跡，以及處理表示在第一接收器處獲得的至少包括對應(yīng)于第二模式的事件的地震數(shù)據(jù)第二參數(shù)的軌跡；b)在處理的數(shù)據(jù)中識別對應(yīng)于部分模式轉(zhuǎn)換的事件；和c)從處理的數(shù)據(jù)中對應(yīng)于部分模式轉(zhuǎn)換的事件的幅度和/或波形獲得信息。
2.如權(quán)利要求1中的方法，其中通過在地層邊界面的第一模式的部分模式轉(zhuǎn)換產(chǎn)生第二模式。
3.如權(quán)利要求1或2中的方法，其中步驟(c)包括將處理的數(shù)據(jù)中對應(yīng)于部分模式轉(zhuǎn)換的事件的幅度進(jìn)行歸一化。
4.如權(quán)利要求3中的方法，其中步驟(c)包括相對于零時(shí)間延遲時(shí)的處理的數(shù)據(jù)的幅度將處理的數(shù)據(jù)中對應(yīng)于部分模式轉(zhuǎn)換的事件的幅度進(jìn)行歸一化。
5.如權(quán)利要求3或4中的方法，進(jìn)一步包括相對于對應(yīng)部分模式轉(zhuǎn)換事件的歸一化的幅度對表示第二參數(shù)的軌跡進(jìn)行歸一化的步驟。
6.如權(quán)利要求1、2、3、4或5中的方法，進(jìn)一步包括對于第二接收器得到的地震數(shù)據(jù)重復(fù)步驟(a)、(b)和(c)的步驟，第二接收器與第一接收器的位置不同。
7.如前述權(quán)利要求中任一項(xiàng)的方法，其中地震數(shù)據(jù)的第一參數(shù)是各個(gè)接收器得到的質(zhì)點(diǎn)運(yùn)動的垂直分量。
8.如權(quán)利要求7中的方法，其中地震數(shù)據(jù)的第二參數(shù)是各個(gè)接收器得到的質(zhì)點(diǎn)運(yùn)動的水平分量。
9.如權(quán)利要求7中的方法，其中地震數(shù)據(jù)的第二參數(shù)是各個(gè)接收器得到的在地震能量源與各個(gè)接收器之間的方向上的質(zhì)點(diǎn)運(yùn)動的分量。
10.如權(quán)利要求1至6中的任何一項(xiàng)的方法，其中地震數(shù)據(jù)的第一和第二參數(shù)中的至少一個(gè)是通過分解各個(gè)接收器所得到的地震數(shù)據(jù)獲得的。
11.如前述權(quán)利要求中任一項(xiàng)的方法，其中步驟(a)包括將表示第一接收器得到的地震數(shù)據(jù)第一參數(shù)的軌跡與表示第一接收器得到的地震數(shù)據(jù)第二參數(shù)的軌跡進(jìn)行相互關(guān)聯(lián)。
12.如權(quán)利要求1-10中的任何一項(xiàng)的方法，其中步驟(a)包括將表示第一接收器得到的地震數(shù)據(jù)第一參數(shù)的軌跡從表示第一接收器得到的地震數(shù)據(jù)第二參數(shù)的軌跡中去卷積。
13.一種地震勘測的方法，包括將以第一模式傳播的地震能量導(dǎo)向到海床層的邊界面，使得在該邊界面處發(fā)生地震能量轉(zhuǎn)換成第二模式的部分模式轉(zhuǎn)換；在一個(gè)或多個(gè)接收器得到包括第一和第二模式地震能量的地震數(shù)據(jù)；和根據(jù)權(quán)利要求1至12中任一項(xiàng)所定義的方法處理該地震數(shù)據(jù)。
14.一種用于處理包括第一和第二模式的地震能量的地震數(shù)據(jù)的設(shè)備，該設(shè)備包括a)用于處理表示在第一接收器處獲得的至少包括對應(yīng)于第一模式的事件的地震數(shù)據(jù)的第一參數(shù)的軌跡，和處理表示在第一接收器處獲得的至少包括對應(yīng)于第二模式的事件的地震數(shù)據(jù)的第二參數(shù)的軌跡的裝置；b)用于在處理的數(shù)據(jù)中識別對應(yīng)于部分模式轉(zhuǎn)換的事件的裝置；和c)用于從處理的數(shù)據(jù)中對應(yīng)于部分模式轉(zhuǎn)換的事件的幅度和/或波形獲得信息的裝置。
15.如權(quán)利要求14的裝置，包括歸一化裝置，用于將對應(yīng)于部分模式轉(zhuǎn)換的事件的被處理數(shù)據(jù)中的幅度進(jìn)行歸一化。
16.如權(quán)利要求15的裝置，其中所述歸一化裝置適合于相對零時(shí)間延遲時(shí)的相互關(guān)聯(lián)的數(shù)據(jù)的幅度將對應(yīng)于部分模式轉(zhuǎn)換的事件的相互關(guān)聯(lián)的數(shù)據(jù)中的幅度進(jìn)行歸一化。
17.如權(quán)利要求14至16中的任何一項(xiàng)的裝置，包括用于相對于與部分模式轉(zhuǎn)換對應(yīng)的事件的歸一化的幅度將表示第二參數(shù)的軌跡進(jìn)行歸一化的裝置。
18.如權(quán)利要求14至17中的任何一項(xiàng)的裝置，包括可編程數(shù)據(jù)處理器。
19.一種存儲媒體，包含用于如權(quán)利要求18所述裝置的數(shù)據(jù)處理器的程序。
全文摘要
一種處理包括地震能量第一和第二模式的地震數(shù)據(jù)的方法，其中通過在海床層邊界面的第一模式的部分模式轉(zhuǎn)換產(chǎn)生第二模式，該方法包括將在接收器得到的并包括對應(yīng)于第一模式的軌跡與在相同接收器得到的并包括對應(yīng)于第二模式的軌跡進(jìn)行相互關(guān)聯(lián)的步驟。識別在相互關(guān)聯(lián)的數(shù)據(jù)中對應(yīng)于部分模式轉(zhuǎn)換的事件，例如相對于零時(shí)間延遲的相互關(guān)聯(lián)記錄圖中峰值的幅度，歸一化該事件的幅度?？梢詮慕邮掌麝嚵兄薪邮掌鞯臍w一化的相互關(guān)聯(lián)記錄圖中得出關(guān)于由該層所產(chǎn)生的靜態(tài)偏移效果的信息和/或關(guān)于矢量失真的信息。
文檔編號G01V1/36GK1656388SQ03802203
公開日2005年8月17日 申請日期2003年1月20日 優(yōu)先權(quán)日2002年1月25日
發(fā)明者約翰·O·A·羅博特森, 安德魯·柯蒂斯, 德克-簡·范曼恩 申請人:維斯特恩格科地震控股有限公司