本公開屬于地震勘探領域，具體涉及三維地震資料多塊拼接處理時差校正方法和裝置。
背景技術：
：在地震勘探中，地震資料的拼接處理進入廣泛應用階段。地震勘探(特別是油氣地震勘探)時單一區(qū)塊的面積多數(shù)比較小，從一百到數(shù)百平方千米不等。隨著地震勘探的不斷深入，這種“郵票式”的資料影響了對地下地質構造的整體宏觀認識與整體構造解釋。目前多數(shù)油田將小勘探區(qū)塊資料進行拼接，以便進行大區(qū)塊的地質研究。由于區(qū)塊資料間存在同一地層波組的地震時差，大面積拼接在一起時產生相位錯動，影響了資料的成像效果，尤其是需要多塊資料拼接時，這一現(xiàn)象更為嚴重。對于地震時差的分析問題，許多學者從不同的角度利用不同的方法作了大量的研究。但從地震資料多區(qū)塊拼接而言，目前采用的方法很多，方法不同其效果存在差異，總結起來主要有如下幾種方法：(1)時移分析法：該方法通過區(qū)塊間的同一地層進行對比，分析兩個地層在剖面上的時差。這種方法在選取地質層時存在人為因素，雖然能夠解決時差問題，但存在較大的誤差。(2)互均衡一致性處理法：該方法通過資料間的地震波組振幅、頻率以及相位等屬性，進行互均衡一致性處理，當波組屬性達到基本一致時分析出資料間的時差，校正后消除時差的影響。(3)匹配濾波分析法：該方法根據(jù)待校正資料的子波特性，分別進行子波分析與提取，將不同的子波進行匹配濾波處理，使得子波特性趨于一致，然后 對資料進行子波褶積處理，得到拼接區(qū)塊的新數(shù)據(jù)體，從而實現(xiàn)區(qū)塊間的時差校正。技術實現(xiàn)要素：本發(fā)明的目的在于提供一種新的三維地震資料多塊拼接處理時差校正的方法和裝置，有效地進行資料拼接時差計算與校正，從而提高地震資料拼接效果。根據(jù)一方面，提出了一種三維地震資料多塊拼接處理時差校正方法，該方法包括以下步驟：從待拼接資料區(qū)塊中確定主區(qū)塊作為起始區(qū)塊，并確定主區(qū)塊為零時差；確定拼接的方向及路線；根據(jù)拼接的方向及路線確定待校正區(qū)塊；獲得使得起始區(qū)塊及待校正區(qū)塊的地震信號的互相關最大的時差，作為起始區(qū)塊與待校正區(qū)塊之間的時差；根據(jù)所獲得的時差，對待校正區(qū)塊的地震信號進行時移，以實現(xiàn)該待校正區(qū)塊與起始區(qū)塊之間的時差校正；以及在所述待校正區(qū)塊不是最后一個待拼接資料區(qū)塊的情況下，根據(jù)所確定的拼接方向和路線確定新的起始區(qū)塊和新的待校正區(qū)塊，重新執(zhí)行上述獲得時差的步驟和時差校正的步驟。另一方面，提出了一種三維地震資料多塊拼接處理時差校正裝置，該裝置包括：用于從待拼接資料區(qū)塊中確定主區(qū)塊作為起始區(qū)塊，并確定主區(qū)塊為零時差的部件；用于確定拼接的方向及路線的部件；用于根據(jù)拼接的方向及路線確定待校正區(qū)塊的部件；用于獲得使得起始區(qū)塊及待校正區(qū)塊的地震信號的互相關最大的時差，作為起始區(qū)塊與待校正區(qū)塊之間的時差的部件；用于根據(jù)所獲得的時差，對待校正區(qū)塊的地震信號進行時移，以實現(xiàn)該待校正區(qū)塊與起始區(qū)塊之間的時差校正的部件；以及用于在所述待校正區(qū)塊不是最后一個待拼接資料區(qū)塊的情況下，根據(jù)所確定的拼接方向和路線確定新的起始區(qū)塊和新的待校正區(qū)塊，重新執(zhí)行上述獲得時差和時差校正的部件。與現(xiàn)有技術相比，本發(fā)明各方面將多塊三維小資料進行拼接處理，區(qū)塊邊界無波組相位錯動痕跡，邊界處同相地質波組達到連續(xù)一致性，地質構造合理， 地層反射強弱關系分明，波組特征明顯，資料信噪比較高，構造成像清楚，斷層斷點清晰，潛山構造完整，不整合強弱及接觸明確，易于后期的整體地質構造研究。附圖說明通過結合附圖對本公開示例性實施方式進行更詳細的描述，本公開的上述以及其它目的、特征和優(yōu)勢將變得更加明顯，其中，在本公開示例性實施方式中，相同的參考標號通常代表相同部件。圖1示出了根據(jù)本發(fā)明的一個實施例的三維地震資料多塊拼接處理時差校正方法的流程圖。圖2示出了根據(jù)一個應用示例的拼接示意圖。圖3示出了根據(jù)一個應用示例的拼接示意圖。圖4示出了根據(jù)一個應用示例的拼接示意圖。圖5示出了根據(jù)一個應用示例的拼接示意圖。圖6示出了未進行時差校正處理的剖面示意圖。圖7示出了針對于圖6相同的待拼接資料區(qū)塊，采用本發(fā)明實施例方法進行時差校正處理后的剖面示意圖。圖8是采用其它方法校正處理后的剖面示意圖圖9示出了針對于圖8相同的待拼接資料區(qū)塊，采用本發(fā)明實施例方法進行時差校正處理后的剖面示意圖。圖10是采用其它方法校正處理后的剖面示意圖圖11示出了針對于圖10相同的待拼接資料區(qū)塊，采用本發(fā)明實施例方法進行時差校正處理后的剖面示意圖。圖12是利用本發(fā)明實施例的方法在某油區(qū)采用21塊地震資料進行拼接處理的拼接方向與路線意圖。圖13是按照圖12所確定的拼接方向與路線的6塊資料拼接線(箭頭為拼 接處)效果示意圖。圖14是按照圖12所確定的拼接方向與路線的4塊資料拼接線(箭頭為拼接處)效果示意圖。具體實施方式下面將參照附圖更詳細地描述本公開的優(yōu)選實施方式。雖然附圖中顯示了本公開的優(yōu)選實施方式，然而應該理解，可以以各種形式實現(xiàn)本公開而不應被這里闡述的實施方式所限制。相反，提供這些實施方式是為了使本公開更加透徹和完整，并且能夠將本公開的范圍完整地傳達給本領域的技術人員。圖1示出了根據(jù)本發(fā)明的一個實施例的三維地震資料多塊拼接處理時差校正方法的流程圖，該方法包括以下步驟：步驟101，從待拼接資料區(qū)塊中確定主區(qū)塊作為起始區(qū)塊，并確定主區(qū)塊為零時差；步驟102，確定拼接的方向及路線；步驟103，根據(jù)拼接的方向及路線確定待校正區(qū)塊；步驟104，獲得使得起始區(qū)塊及待校正區(qū)塊的地震信號的互相關最大的時差，作為起始區(qū)塊與待校正區(qū)塊之間的時差；步驟105根據(jù)所獲得的時差，對待校正區(qū)塊的地震信號進行時移，以實現(xiàn)該待校正區(qū)塊與起始區(qū)塊之間的時差校正；以及步驟106，在所述待校正區(qū)塊不是最后一個待拼接資料區(qū)塊的情況下，根據(jù)所確定的拼接方向和路線確定新的起始區(qū)塊和新的待校正區(qū),重新執(zhí)行上述步驟104和105，即返回步驟104重新執(zhí)行。該實施例利用了最大互相關能量法來進行資料區(qū)塊拼接，使得拼接質量得到了提高。在一個示例中，本領域技術人員可以根據(jù)處理項目要求，來確定所需要拼接的資料區(qū)塊的塊數(shù)，以便滿足地質成像的需求。例如可以是處理區(qū)域所包含 的所有資料區(qū)塊的數(shù)目。確定主區(qū)塊作為起始區(qū)塊在一個示例中，可以在所有拼接的資料內，分析選取地質成像處理效果較好、地質層位明確的區(qū)塊作為主區(qū)塊，同時確定此主區(qū)塊為零時差，其它區(qū)塊以此為參考計算時差。在確定主區(qū)塊時可以主要參考地質解釋人員對資料的客觀評價，盡量選取在這些資料區(qū)塊內相對處理效果較好的作為主區(qū)塊。確定拼接的方向及路線在一個示例中，本領域技術人員可以根據(jù)需要任意制定確定拼接方向及路線的原則。在一個示例中，可以以主區(qū)塊為起始，向四周設計拼接的方向，方向次序可以是相對于主區(qū)塊由近及遠。區(qū)塊較多時可以設計出拼接路線,例如盡可能滿足由近及遠、不跳越區(qū)塊、使路線能夠形成回路，以滿足合理有效拼接的要求。待校正的區(qū)塊以及迭代過程中新的主區(qū)塊和新的待校正的區(qū)塊均可以根據(jù)拼接的方向和路線來確定。獲得時差。發(fā)明人認識到，通常相鄰區(qū)塊的三維地震資料有一部分是重疊的，重疊部分對應的是地下同一區(qū)域的地層，地震反射波的性質基本相同，具有良好的相關性。互相關函數(shù)是對兩個地震信號相似程度的比較，因此，可以利用互相關函數(shù)對同一位置的兩個地震道進行互相關分析，定量描述兩個信號的時間延遲?；ハ嚓P時差定量識別技術的方法原理如下：設x1(t)和x2(t)分別為兩個資料區(qū)塊的地震信號(例如本實施例中的起始區(qū)塊及待校正區(qū)塊的地震信號)，期望地震信號為s(t)，n1(t)和n2(t)分別為相互獨立的均值為0、方差為1的高斯白噪聲，td為延遲時間(即時差)，α為衰減系數(shù)，則地震信號x1(t)和x2(t)可表示為：x1(t)=s(t)+n1(t)x2(t)=αs(t-td)+n2(t)]]>x1(t)和x2(t)的互相關函數(shù)為：Rx1x2(τ)=E[x1(t)x2(t+τ)]=αRss(τ-td)]]>由于|Rss(τ)|≤Rss(0)，因此Rss(τ)在τ-td＝0時為極大值，以此關系可以求得使得最大的td，此時的td就是地震信號x2(t)的校正時移量，也就是地震信號x2(t)相對于地震信號x1(t)的時差。將地震信號x2(t)時移td就實現(xiàn)了時差校正。本實施例中，還可以通過以下示例中的一致性判斷或閉合點分析中的一者或兩者，來進一步提高拼接精度。在一個示例中，在步驟105中進行時差校正之后，還可進一步分析時差校正后的起始資料區(qū)塊及待校正資料區(qū)塊之間的同相波組的頻率、相位與波形的一致性(例如，可以通過專業(yè)分析軟件來分析同一波組的波形對應有無錯位、頻率有無缺失、相位是否一致)。如果一致性高(即波形對應錯位少、頻率缺失少、相位一致性高)，說明時差計算準確合理，可以繼續(xù)執(zhí)行步驟106。反之，如果一致性不高，則可以分析原因，并根據(jù)原因從步驟102-105中的一個開始重復執(zhí)行(或者重復執(zhí)行步驟102-105中的一個或多個)，直至一致性達到要求為止。在一個示例中，如果待拼接的資料區(qū)塊中存在閉合點(就是兩塊資料區(qū)塊位置完全一樣的完全重疊情況)時(如圖5所示)，可以在拼接前確定閉合點，并在步驟105中進行時差校正后判斷所確定的閉合點的閉合程度(即判斷原本完全重疊的資料區(qū)塊在時差校正后的重疊程度)。如閉合程度(即重疊程度)達不到要求(例如不能完全重疊，或重疊程度低于要求)，就需要返回分析前面的拼接點，找出原因，并根據(jù)原因從步驟102-105中的一個開始重復執(zhí)行(或者重復執(zhí)行步驟102-105中的一個或多個)，直至達到閉合(即完全重疊)為止。在一個示例中，可在所有待拼接的資料區(qū)塊拼接完成后，作為一個整體數(shù) 據(jù)輸出所有拼接后的資料區(qū)塊，以便進行后期的處理。本發(fā)明實施例的方法時差計算精度高拼接效果好，實現(xiàn)了利用矢量法進行拼接區(qū)塊間的時差校正，有效的消除了區(qū)塊間的時差，使得資料連片拼接同相疊加，克服了小區(qū)塊間的邊界不一致性問題，提高了拼接資料的地質成像質量。具體應用示例為便于理解，以下給出一些具體應用示例。本領域技術人員應理解，這些應用示例僅為示例性的，其中的任何細節(jié)都非意在限制本發(fā)明。圖2～5是4個不同數(shù)量區(qū)塊拼接方式示意圖，圖2、圖3是3塊拼接示意圖，它們的主區(qū)塊位置不同，拼接方向不同；圖4是5塊資料拼接示意圖，圖中標注了拼接的方向和路線；圖5是6塊資料拼接示意圖，圖中標注了拼接的方向和路線，該拼接方式出現(xiàn)了2個拼接閉合點，此種方式拼接時要考慮閉合點的閉合問題。圖6、圖7是同一塊拼接資料的對比圖。圖6是3塊資料拼接(白色箭頭處為拼接處)，可以看到沒有進行時差校正處理時，區(qū)塊邊界由于時差的存在，使得在邊界處波組雜亂，地層消失，無法進行地層連續(xù)解釋；圖7是采用本發(fā)明實施例方法進行時差處理后，可以看出，邊界處的波組連續(xù)完整，地層連續(xù)，可進行有效解釋。圖8、圖9是一個使用本方法與其它方法的效果對比，是2塊資料拼接(白色箭頭處為拼接處)，圖8是其它時差校正方法處理的剖面，可以看出在拼接處存在誤差，尤其是深部的波組不連續(xù)；地層變得模糊，進行地層連續(xù)解釋存在難度；圖9是采用本發(fā)明實施例的方法進行時差處理后的剖面，可以看出，邊界處的波組連續(xù)完整，地層連續(xù)，可進行追蹤與有效解釋。圖10、圖11是一個使用本發(fā)明實施例的與其它方法的效果對比，是2塊資料拼接(白色矩形框內為拼接處)，圖10是其它方法處理結果，可以看到區(qū)塊拼接處由于時差的存在，使得在邊界處波組不能連續(xù)，地層變得模糊，進行地層連續(xù)解釋存在難度；圖11是采用本發(fā)明實施例的方法進行時差處理后，可以看出，邊界處的波組連續(xù)完整，地層連續(xù)，可進行追蹤與有效解釋。圖12-14是利用本發(fā)明實施例的方法在某油區(qū)采用21塊地震資料進行拼接處理的示意圖。21塊小三維資料拼接的較大資料，拼接區(qū)塊多、難度大，采用本發(fā)明實施例的方法取得了很好的效果。圖12是拼接方向與路線圖，圖13是6塊拼接的一條剖面(白色箭頭為拼接處)，圖14是4塊拼接的一條剖面(白色箭頭為拼接處)。采用本發(fā)明實施例方法進行時差處理后，可以看出：邊界處的波組連續(xù)，地質構造完整，不存在拼接邊界效應，能夠進行地層連續(xù)追蹤與有效解釋。本公開可以是系統(tǒng)、方法和/或計算機程序產品。計算機程序產品可以包括計算機可讀存儲介質，其上載有用于使處理器實現(xiàn)本公開的各個方面的計算機可讀程序指令。計算機可讀存儲介質可以是可以保持和存儲由指令執(zhí)行設備使用的指令的有形設備。計算機可讀存儲介質例如可以是――但不限于――電存儲設備、磁存儲設備、光存儲設備、電磁存儲設備、半導體存儲設備或者上述的任意合適的組合。計算機可讀存儲介質的更具體的例子(非窮舉的列表)包括：便攜式計算機盤、硬盤、隨機存取存儲器(RAM)、只讀存儲器(ROM)、可擦式可編程只讀存儲器(EPROM或閃存)、靜態(tài)隨機存取存儲器(SRAM)、便攜式壓縮盤只讀存儲器(CD-ROM)、數(shù)字多功能盤(DVD)、記憶棒、軟盤、機械編碼設備、例如其上存儲有指令的打孔卡或凹槽內凸起結構、以及上述的任意合適的組合。這里所使用的計算機可讀存儲介質不被解釋為瞬時信號本身，諸如無線電波或者其他自由傳播的電磁波、通過波導或其他傳輸媒介傳播的電磁波(例如，通過光纖電纜的光脈沖)、或者通過電線傳輸?shù)碾娦盘?。這里所描述的計算機可讀程序指令可以從計算機可讀存儲介質下載到各個計算/處理設備，或者通過網絡、例如因特網、局域網、廣域網和/或無線網下載到外部計算機或外部存儲設備。網絡可以包括銅傳輸電纜、光纖傳輸、無線傳輸、路由器、防火墻、交換機、網關計算機和/或邊緣服務器。每個計算/處理設備中的網絡適配卡或者網絡接口從網絡接收計算機可讀程序指令，并轉發(fā) 該計算機可讀程序指令，以供存儲在各個計算/處理設備中的計算機可讀存儲介質中。用于執(zhí)行本公開操作的計算機程序指令可以是匯編指令、指令集架構(ISA)指令、機器指令、機器相關指令、微代碼、固件指令、狀態(tài)設置數(shù)據(jù)、或者以一種或多種編程語言的任意組合編寫的源代碼或目標代碼，所述編程語言包括面向對象的編程語言—諸如Smalltalk、C++等，以及常規(guī)的過程式編程語言—諸如“C”語言或類似的編程語言。計算機可讀程序指令可以完全地在用戶計算機上執(zhí)行、部分地在用戶計算機上執(zhí)行、作為一個獨立的軟件包執(zhí)行、部分在用戶計算機上部分在遠程計算機上執(zhí)行、或者完全在遠程計算機或服務器上執(zhí)行。在涉及遠程計算機的情形中，遠程計算機可以通過任意種類的網絡—包括局域網(LAN)或廣域網(WAN)—連接到用戶計算機，或者，可以連接到外部計算機(例如利用因特網服務提供商來通過因特網連接)。在一些實施例中，通過利用計算機可讀程序指令的狀態(tài)信息來個性化定制電子電路，例如可編程邏輯電路、現(xiàn)場可編程門陣列(FPGA)或可編程邏輯陣列(PLA)，該電子電路可以執(zhí)行計算機可讀程序指令，從而實現(xiàn)本公開的各個方面。這里參照根據(jù)本公開實施例的方法、裝置(系統(tǒng))和計算機程序產品的流程圖和/或框圖描述了本公開的各個方面。應當理解，流程圖和/或框圖的每個方框以及流程圖和/或框圖中各方框的組合，都可以由計算機可讀程序指令實現(xiàn)。這些計算機可讀程序指令可以提供給通用計算機、專用計算機或其它可編程數(shù)據(jù)處理裝置的處理器，從而生產出一種機器，使得這些指令在通過計算機或其它可編程數(shù)據(jù)處理裝置的處理器執(zhí)行時，產生了實現(xiàn)流程圖和/或框圖中的一個或多個方框中規(guī)定的功能/動作的裝置。也可以把這些計算機可讀程序指令存儲在計算機可讀存儲介質中，這些指令使得計算機、可編程數(shù)據(jù)處理裝置和/或其他設備以特定方式工作，從而，存儲有指令的計算機可讀介質則包括一個制造品，其包括實現(xiàn)流程圖和/或框圖中的一個或多個方框中規(guī)定的功能/動作 的各個方面的指令。也可以把計算機可讀程序指令加載到計算機、其它可編程數(shù)據(jù)處理裝置、或其它設備上，使得在計算機、其它可編程數(shù)據(jù)處理裝置或其它設備上執(zhí)行一系列操作步驟，以產生計算機實現(xiàn)的過程，從而使得在計算機、其它可編程數(shù)據(jù)處理裝置、或其它設備上執(zhí)行的指令實現(xiàn)流程圖和/或框圖中的一個或多個方框中規(guī)定的功能/動作。附圖中的流程圖和框圖顯示了根據(jù)本公開的多個實施例的系統(tǒng)、方法和計算機程序產品的可能實現(xiàn)的體系架構、功能和操作。在這點上，流程圖或框圖中的每個方框可以代表一個模塊、程序段或指令的一部分，所述模塊、程序段或指令的一部分包含一個或多個用于實現(xiàn)規(guī)定的邏輯功能的可執(zhí)行指令。在有些作為替換的實現(xiàn)中，方框中所標注的功能也可以以不同于附圖中所標注的順序發(fā)生。例如，兩個連續(xù)的方框實際上可以基本并行地執(zhí)行，它們有時也可以按相反的順序執(zhí)行，這依所涉及的功能而定。也要注意的是，框圖和/或流程圖中的每個方框、以及框圖和/或流程圖中的方框的組合，可以用執(zhí)行規(guī)定的功能或動作的專用的基于硬件的系統(tǒng)來實現(xiàn)，或者可以用專用硬件與計算機指令的組合來實現(xiàn)。以上已經描述了本公開的各實施例，上述說明是示例性的，并非窮盡性的，并且也不限于所披露的各實施例。在不偏離所說明的各實施例的范圍和精神的情況下，對于本
技術領域：
的普通技術人員來說許多修改和變更都是顯而易見的。本文中所用術語的選擇，旨在最好地解釋各實施例的原理、實際應用或對市場中的技術的技術改進，或者使本
技術領域：
的其它普通技術人員能理解本文披露的各實施例。當前第1頁1 2 3