本發(fā)明涉及一種地震數(shù)據(jù)處理方法，特別是關(guān)于一種用于時移地震的九面元一致性并行處理方法。

背景技術(shù)：
20世紀80年代末期，時移地震方法被提出并進行了現(xiàn)場應(yīng)用，其主要思想是通過不同時間的多次地震采集，觀測地震響應(yīng)的差異進而確定油藏隨時間的變化，即動態(tài)變化。時移地震是目前進行油氣藏監(jiān)測，合理調(diào)整開發(fā)方案、提高油氣采收率的有效手段，能夠分析研究油氣藏開采過程中造成的儲層流體運動、流體成分、流體飽和度、壓力、孔隙度和溫度等油氣藏特性的變化。面元一致性處理是時移地震處理的重要環(huán)節(jié)之一。由于地震資料重復(fù)采集的時間不同，地面設(shè)施、技術(shù)裝備等因素可能發(fā)生了變化，使得觀測系統(tǒng)、采集參數(shù)很難和原有數(shù)據(jù)完全一樣，另外，目前已有相當(dāng)多的三維地震數(shù)據(jù)，它們不是針對時移地震油藏監(jiān)測而采集的，采集和處理參數(shù)與重復(fù)采集時的參數(shù)很難完全相同，對于三維地震而言主要表現(xiàn)為反射面元的大小和位置不一樣，對于二維數(shù)據(jù)則表現(xiàn)為反射點的位置不同。為了使重復(fù)采集的地震資料具有可比性，需要把來自地下不同反射面元或反射點的地震數(shù)據(jù)校正到同樣的反射面元或反射點，這一處理又稱之為面元重置。面元一致性處理的方法很多，目前主要有兩種方法。第一種是相關(guān)抽道法，這種方法利用最鄰近面元網(wǎng)格上的地震道作為新面元上的地震道，它的缺點是沒有考慮相鄰道振幅和相位的變化；第二種是插值法(線性插值法、F-K域插值法、T-X域動態(tài)求差插值法)，該方法考慮了振幅和相位的變化，但只在線性同相軸和規(guī)則插值時比較實用。每種方法的處理目的都是保證盡量不損害地震振幅信息的前提下，使時移地震數(shù)據(jù)有最佳一致性。地震資料的一致性越高，時移地震解釋的準確度就越高。影響地震資料的一致性的因素包括面元中心點，方位角，偏移距以及震源和檢波器的位置。減小面元中心點位置距離，或者方位角差，或者震源間和檢波器間的距離，都能夠提高地震資料的一致性。對不同時間點采集的兩套地震數(shù)據(jù)進行面元一致性處理時，將采集時間早的數(shù)據(jù)稱為基準數(shù)據(jù)，將采集時間晚的數(shù)據(jù)稱為監(jiān)測數(shù)據(jù)。常用的方法可以在兩套數(shù)據(jù)中依次各取一個面元的地震數(shù)據(jù)進行一致性處理，也可以在基準數(shù)據(jù)中選取一個面元，在監(jiān)測數(shù)據(jù)中選取與該面元相鄰的九個面元進行一致性處理。使用九個面元進行處理，可以更好地提高地震資料的一致性，但對算法、計算資源及時間有著更高的要求，尤其是疊前地震資料數(shù)據(jù)量非常大，實際處理中遇到許多問題。

技術(shù)實現(xiàn)要素：
針對上述問題，本發(fā)明的目的是提供一種用于時移地震的九面元一致性并行處理方法，該方法通過面元重置，提高地震數(shù)據(jù)的一致性，從而為后續(xù)的地震資料處理、解釋和反演提供一個好的數(shù)據(jù)基礎(chǔ)，提高對油氣藏變化監(jiān)測的準確度，從而合理調(diào)整開發(fā)方案、提高油氣采收率。為實現(xiàn)上述目的，本發(fā)明采取以下技術(shù)方案：一種用于時移地震的九面元一致性并行處理方法，其特征在于包括以下步驟：1)根據(jù)已有的同一地區(qū)不同時間采集的兩套地震數(shù)據(jù)，分別作為基準數(shù)據(jù)Db和監(jiān)測數(shù)據(jù)Dm，采用并行節(jié)點k+1個，第1個并行節(jié)點作為主節(jié)點，其余并行節(jié)點為從節(jié)點；從存儲介質(zhì)上向主節(jié)點內(nèi)存中讀入基準數(shù)據(jù)Db中第n個CMP道集，其對應(yīng)的面元為B(i，j)，第n個CMP道集中含有m1道地震數(shù)據(jù)；2)從存儲介質(zhì)讀取監(jiān)測數(shù)據(jù)Dm中的第i-1條、第i條和第i+1條三條主測線數(shù)據(jù)，將數(shù)據(jù)平均分為k段，分別讀入第2～k+1個并行從節(jié)點內(nèi)存中；3)第2～k+1個從節(jié)點并行查找與面元B(i，j)相對應(yīng)的九面元地震數(shù)據(jù)，并將查找到的九面元地震數(shù)據(jù)發(fā)送給主節(jié)點；主節(jié)點得到兩塊數(shù)據(jù)，一塊為基準數(shù)據(jù)Db中第n個CMP道集，含有m道地震數(shù)據(jù)d1；另一塊為監(jiān)測數(shù)據(jù)Dm中的九面元數(shù)據(jù)，含有m9道地震數(shù)據(jù)d9；4)對步驟3)中的兩塊數(shù)據(jù)進行抽取，得到新基準數(shù)據(jù)和新監(jiān)測數(shù)據(jù)完成九面元一致性并行處理。所述步驟4)中，所述兩塊數(shù)據(jù)抽取方法如下：(1)從數(shù)據(jù)d1中選取第l道，在數(shù)據(jù)d9中查找與第l道數(shù)據(jù)的震源距離與檢波器距離之和Dis<Error的道數(shù)據(jù)de；其中，Error為預(yù)先設(shè)定的閾值；(2)如果數(shù)據(jù)de的道數(shù)大于1道，則在數(shù)據(jù)de中查找與第l道數(shù)據(jù)中心點距離最小的道，記為數(shù)據(jù)dm；如果數(shù)據(jù)de的道數(shù)等于1道，則輸出數(shù)據(jù)d1中的第l道數(shù)據(jù)為面元一致性處理后的新基準數(shù)據(jù)輸出de的道為面元一致性處理后的新監(jiān)測數(shù)據(jù)如果數(shù)據(jù)de中沒有數(shù)據(jù)，則轉(zhuǎn)至步驟(4)；(3)如果數(shù)據(jù)dm的道數(shù)大于1道，則在數(shù)據(jù)dm中查找與第l道數(shù)據(jù)方位角相差最小的道，輸出該道為面元一致性處理后的新監(jiān)測數(shù)據(jù)輸出數(shù)據(jù)d1中的第l道數(shù)據(jù)為面元一致性處理后的新基準數(shù)據(jù)(4)如果l<m，則l＝l+1，則返回步驟(1)；如果l＝m，則n＝n+1，返回步驟1)，讀取下一個CMP道集，面元位置也隨之移動，直至基準數(shù)據(jù)Db中所有CMP道集處理完畢。所述存儲介質(zhì)采用磁帶或磁盤。本發(fā)明由于采取以上技術(shù)方案，其具有以下優(yōu)點：1、由于傳統(tǒng)的方法是在兩套數(shù)據(jù)中依次各取一個面元的地震數(shù)據(jù)進行一致性處理，而本發(fā)明所提的方法是在監(jiān)測數(shù)據(jù)中選取九個面元地震數(shù)據(jù)進行一致性處理，提高了處理后地震資料的一致性，為后續(xù)的地震資料處理、解釋和反演提供了良好的數(shù)據(jù)基礎(chǔ)，從而提高對油氣藏變化監(jiān)測的準確度，能夠合理調(diào)整開發(fā)方案，最終提高油氣采收率。2、本發(fā)明采用的九面元并行處理方法對單個節(jié)點的內(nèi)存及計算速度要求不高，但通過并行，能夠處理大量的地震數(shù)據(jù)，快速查找面元數(shù)據(jù)，提高了地震數(shù)據(jù)處理速度，能夠滿足實際生產(chǎn)的需要。本發(fā)明可以廣泛在油藏監(jiān)測中領(lǐng)域中應(yīng)用。附圖說明圖1是本發(fā)明的基準數(shù)據(jù)疊前CMP道集示意圖；圖2是本發(fā)明的面元示意圖，說明工區(qū)內(nèi)面元的分布以及九面元的定義；圖3是本發(fā)明的整體流程示意圖；圖4是本發(fā)明的數(shù)據(jù)抽取方法流程示意圖；圖5是本發(fā)明一致性處理后的新基準數(shù)據(jù)疊前CMP道集；圖6是本發(fā)明一致性處理后的基準數(shù)據(jù)和監(jiān)測數(shù)據(jù)疊后剖面；圖7a是本發(fā)明一致性處理前的原始震源位置示意圖；圖7b是本發(fā)明一致性處理后的震源位置示意圖；圖8a是本發(fā)明一致性處理前的原始偏移距示意圖；圖8b是本發(fā)明一致性處理后的偏移距示意圖；圖9a是本發(fā)明一致性處理前的原始方位角示意圖；圖9b是本發(fā)明一致性處理后的方位角示意圖。具體實施方式本發(fā)明提供一種用于時移地震的九面元一致性并行處理方法，該方法針對同一區(qū)塊不同時間采集的地震數(shù)據(jù)進行處理，地震資料采集后，經(jīng)過前期處理，生成疊前共中心點道集CMP(如圖1所示)，每個道集對應(yīng)一個面元，每個道集中含有m道地震數(shù)據(jù)。如圖2所示，對于某地區(qū)，沿主測線方向有I個面元，沿聯(lián)絡(luò)測線方向有J個面元。對于任一面元B(i，j)，存在8個與其相鄰的面元B(i±1，j±1)，加上自身一共九個面元。當(dāng)面元處于工區(qū)邊界上時，則只存在5個相鄰面元。地震數(shù)據(jù)在存儲介質(zhì)(如磁帶、磁盤等)中是按照面元順序存儲的，每個面元內(nèi)的道數(shù)不同，因而面元數(shù)據(jù)存儲位置沒有特定規(guī)律。面元一致性處理時，需要在數(shù)據(jù)內(nèi)搜索與面元B(i，j)相對應(yīng)的九個面元，而地震數(shù)據(jù)特別巨大，對處理時間也有一定的要求。采用合適的方法讀入九個面元的地震數(shù)據(jù)并及時、高效完成處理，是實際應(yīng)用中所面臨的重要問題。下面結(jié)合附圖和實施例對本發(fā)明進行詳細的描述。如圖3所示，本發(fā)明用于時移地震的九面元一致性并行處理方法具有以下步驟：1)根據(jù)已有的同一地區(qū)不同時間采集的兩套地震數(shù)據(jù)，分別作為基準數(shù)據(jù)Db和監(jiān)測數(shù)據(jù)Dm，采用并行節(jié)點k+1個，第1個并行節(jié)點作為主節(jié)點，其余并行節(jié)點為從節(jié)點。從存儲介質(zhì)(如磁帶、磁盤等)上向主節(jié)點內(nèi)存中讀入基準數(shù)據(jù)Db中第n個CMP道集，其對應(yīng)的面元為B(i，j)，第n個CMP道集中含有m1道地震數(shù)據(jù)；其中，k的取值由地震數(shù)據(jù)的大小及計算機節(jié)點的內(nèi)存大小所決定。2)從存儲介質(zhì)(如磁帶、磁盤等)讀取監(jiān)測數(shù)據(jù)Dm中的第i-1條、第i條和第i+1條三條主測線數(shù)據(jù)，將數(shù)據(jù)平均分為k段，分別讀入第2～k+1個并行從節(jié)點內(nèi)存中。3)第2～k+1個從節(jié)點并行查找與面元B(i，j)相對應(yīng)的九面元地震數(shù)據(jù)，并將查找到的九面元地震數(shù)據(jù)發(fā)送給主節(jié)點；主節(jié)點得到兩塊數(shù)據(jù)，一塊為基準數(shù)據(jù)Db中第n個CMP道集，含有m道地震數(shù)據(jù)d1；另一塊為監(jiān)測數(shù)據(jù)Dm中的九面元數(shù)據(jù)，含有m9道地震數(shù)據(jù)d9。4)對步驟3)中的兩塊數(shù)據(jù)進行抽取，得到新基準數(shù)據(jù)和新監(jiān)測數(shù)據(jù)完成九面元一致性并行處理。上述步驟4)中，如圖4所示，兩塊數(shù)據(jù)抽取方法包括以下步驟：(1)從數(shù)據(jù)d1中選取第l道，在數(shù)據(jù)d9中查找與第l道數(shù)據(jù)的震源距離與檢波器距離之和Dis<Error的道數(shù)據(jù)de；其中，Error為預(yù)先設(shè)定的閾值，表示震源距離與檢波器距離之和的最大值，Error的取值根據(jù)地震數(shù)據(jù)的品質(zhì)所決定。(2)如果數(shù)據(jù)de的道數(shù)大于1道，則在數(shù)據(jù)de中查找與第l道數(shù)據(jù)中心點距離最小的道，記為數(shù)據(jù)dm；如果數(shù)據(jù)de的道數(shù)等于1道，則輸出數(shù)據(jù)d1中的第l道數(shù)據(jù)為面元一致性處理后的新基準數(shù)據(jù)輸出de的道為面元一致性處理后的新監(jiān)測數(shù)據(jù)如果道數(shù)據(jù)de中沒有數(shù)據(jù)，即步驟(1)d1中不存在道數(shù)據(jù)de，則轉(zhuǎn)至步驟(4)；(3)如果數(shù)據(jù)dm的道數(shù)大于1道，則在數(shù)據(jù)dm中查找與第l道數(shù)據(jù)方位角相差最小的道，輸出該道為面元一致性處理后的新監(jiān)測數(shù)據(jù)輸出數(shù)據(jù)d1中的第l道數(shù)據(jù)為面元一致性處理后的新基準數(shù)據(jù)(4)如果l<m，則l＝l+1，則返回步驟(1)；如果l＝m，則n＝n+1，返回步驟1)，讀取下一個CMP道集，面元位置也隨之移動，直至基準數(shù)據(jù)Db中所有CMP道集處理完畢。實施例：本實施例中以四個計算機節(jié)點進行并行處理，計算機節(jié)點之間的數(shù)據(jù)通信采用了MessagePassingInterface，即標(biāo)準消息傳遞界面MPI，但本發(fā)明并不限于使用MPI，也可使用其他方法實現(xiàn)。1)從存儲介質(zhì)(如磁帶、磁盤等)上向主節(jié)點內(nèi)存中讀入基準數(shù)據(jù)Db中第n個CMP道集，其對應(yīng)的面元為B(i，j)，該道集中含有m1道地震數(shù)據(jù)。2)從存儲介質(zhì)(如磁帶、磁盤等)讀取監(jiān)測數(shù)據(jù)Dm中的三條主測線數(shù)據(jù)(第i-1，i和i+1條)，將數(shù)據(jù)平均分為3段，分別讀入第2～4從節(jié)點內(nèi)存中。3)第2～4個從節(jié)點并行查找與面元B(i，j)相對應(yīng)的九面元地震數(shù)據(jù)。并將其發(fā)送給主節(jié)點；主節(jié)點得到兩塊數(shù)據(jù)，一塊為基準數(shù)據(jù)Db中第n個CMP道集，含有m道地震數(shù)據(jù)，稱為d1；另一塊為監(jiān)測數(shù)據(jù)Dm中的九面元數(shù)據(jù)，含有m9道地震數(shù)據(jù)，稱為d9。4)對步驟3)中的兩塊數(shù)據(jù)進行抽取，得到新基準數(shù)據(jù)和新監(jiān)測數(shù)據(jù)完成九面元一致性并行處理：(1)從數(shù)據(jù)d1中選取第l道，在數(shù)據(jù)d9中查找與第l道數(shù)據(jù)的震源距離與檢波器距離之和Dis<Error(Error＝100m)的道數(shù)據(jù)de；(2)如果數(shù)據(jù)de的道數(shù)大于1道，在其中查找與第l道數(shù)據(jù)中心點距離最小的道，記為數(shù)據(jù)dm；如果數(shù)據(jù)de的道數(shù)等于1道，則輸出d1中的第l道數(shù)據(jù)為面元重置后的基準數(shù)據(jù)輸出de的道為面元重置好的監(jiān)測數(shù)據(jù)如果de中沒有數(shù)據(jù)，則轉(zhuǎn)至步驟(4)；(3)如果數(shù)據(jù)dm的道數(shù)大于1道，在其中查找與第l道數(shù)據(jù)方位角相差最小的道，輸出該道為面元重置好的監(jiān)測數(shù)據(jù)輸出d1中的第l道數(shù)據(jù)為面元重置后的基準數(shù)據(jù)(4)如果l<m，則l＝l+1，則返回步驟(1)；如果l＝m，則n＝n+1，返回步驟1)，選取下一個CMP道集，面元位置也隨之移動，直至基準數(shù)據(jù)Db中所有CMP道集處理完畢。綜上所述，如圖5所示，本發(fā)明的用于時移地震的九面元一致性并行處理方法，可以得到新的基準數(shù)據(jù)和監(jiān)測數(shù)據(jù)。如圖6所示，為一致性處理后的基準數(shù)據(jù)和監(jiān)測數(shù)據(jù)疊后剖面。為了驗證本發(fā)明的效果，對一致性處理前后的一個CMP道集內(nèi)地震資料的震源位置，偏移距以及方位角信息進行了分析，如圖7a～圖9b所示。經(jīng)過一致性處理后，震源位置重合度得到了提高，偏移距及方位角之差變小，地震資料的一致性得到了提高。上述各實施例僅用于說明本發(fā)明，各部件的結(jié)構(gòu)、尺寸、設(shè)置位置及形狀都是可以有所變化的，在本發(fā)明技術(shù)方案的基礎(chǔ)上，凡根據(jù)本發(fā)明原理對個別部件進行的改進和等同變換，均不應(yīng)排除在本發(fā)明的保護范圍之外。