一種基于近地表層的折射波靜校正處理方法以及系統(tǒng)的制作方法
【專利摘要】本發(fā)明提供一種基于近地表層的折射波靜校正處理方法及系統(tǒng),所述方法包括:采集地震資料;從所述的地震資料中拾取出折射波初至;確定水平基準(zhǔn)面的靜校正量以及浮動基準(zhǔn)面的靜校正量;根據(jù)所述浮動基準(zhǔn)面的靜校正量確定所述的地震資料對應(yīng)的折射剩余靜校正量;根據(jù)所述的水平基準(zhǔn)面的靜校正量以及所述的折射剩余靜校正量確定所述地震資料的最終靜校正量。本發(fā)明能解決近地表層短波長靜校正問題,與其他靜校正方法(如:模型法或?qū)游龇?相結(jié)合最終解決靜校正問題。
【專利說明】一種基于近地表層的折射波靜校正處理方法以及系統(tǒng)
【技術(shù)領(lǐng)域】
[0001]本發(fā)明關(guān)于地球物理勘探【技術(shù)領(lǐng)域】,特別是關(guān)于地震勘探中的靜校正處理技術(shù),具體的講是一種基于近地表層的折射波靜校正處理方法及系統(tǒng)。
【背景技術(shù)】
[0002]地震勘探中的靜校正處理技術(shù)主要是用來解決地震波在近地表介質(zhì)旅行時間變化影響反射波疊加成像的問題。地震發(fā)生后,地震觀測點(diǎn)最先接收到的波稱初至波。由于各種地震波的傳播速度不同,傳播到觀測點(diǎn)的時間也就有先后。利用初至波的折射靜校正方法或?qū)游鲮o校正方法可以反演出相對精確的近地表模型,解決靜校正問題。地震波波前到達(dá)某個觀測點(diǎn),在觀測點(diǎn)上,檢波器檢測到質(zhì)點(diǎn)振動的時刻稱為波的初至?xí)r間,簡稱初至。
[0003]在地球物理地震勘探中,折射波的靜校正處理方法是解決復(fù)雜區(qū)靜校正問題的主要方法之一。然而,折射波的靜校正在長波長靜校正控制方面存在一定的問題,因此在復(fù)雜區(qū),靜校正主要采用模型法靜校正或?qū)游鲮o校正方法來解決長波長靜校正問題,利用折射波靜校正解決短波長靜校正問題。然而,這兩種方法的結(jié)合存在較大的難度。
[0004]因此,如何提出一種新的折射波靜校正方案,其能夠解決近地表層短波長靜校正問題且能與其他靜校正方法(如:模型法或?qū)游龇?相結(jié)合以致最終解決靜校正問題是本領(lǐng)域亟待解決的技術(shù)難題。
【發(fā)明內(nèi)容】
[0005]為了克服現(xiàn)有技術(shù)中的折射波靜校正處理方案存在的上述缺陷,本發(fā)明提供了一種基于近地表層的折射波靜校正處理方法及系統(tǒng),以一次靜校正(模型法或?qū)游龇o校正)為基礎(chǔ)獲得浮動基準(zhǔn)面靜校正,把浮動基準(zhǔn)面靜校正應(yīng)用到拾取的折射波初至上,確定出折射剩余靜校正,根據(jù)折射剩余靜校正與水平基準(zhǔn)面一次靜校正確定出最終靜校正量,能夠有效解決復(fù)雜近地表地區(qū)長波長靜校正問題和短波長靜校正問題。
[0006]本發(fā)明的目的之一是,提供一種基于近地表層的折射波靜校正處理方法,包括:采集地震資料;從所述的地震資料中拾取出折射波初至;確定水平基準(zhǔn)面的靜校正量以及浮動基準(zhǔn)面的靜校正量;根據(jù)所述浮動基準(zhǔn)面的靜校正量確定所述的地震資料對應(yīng)的折射剩余靜校正量;根據(jù)所述的水平基準(zhǔn)面的靜校正量以及所述的折射剩余靜校正量確定所述地震資料的最終靜校正量,其中根據(jù)所述浮動基準(zhǔn)面的靜校正量確定所述的地震資料對應(yīng)的折射剩余靜校正量包括:將所述浮動基準(zhǔn)面的靜校正量應(yīng)用到所述的折射波初至,得到應(yīng)用浮動基準(zhǔn)面的靜校正量的初至;根據(jù)所述應(yīng)用浮動基準(zhǔn)面的靜校正量的初至確定穩(wěn)定折射層的初至;根據(jù)所述穩(wěn)定折射層的初至確定所述的地震資料對應(yīng)的折射剩余靜校正量,所述的折射剩余靜校正量包括炮點(diǎn)的折射剩余靜校正量以及檢波點(diǎn)的折射剩余靜校正量。
[0007]本發(fā)明的目的之一是,提供了一種基于近地表層的折射波靜校正處理系統(tǒng),包括:地震資料采集裝置,用于采集地震資料;折射波初至拾取裝置,用于從所述的地震資料中拾取出折射波初至;基準(zhǔn)面靜校正量確定裝置,用于確定水平基準(zhǔn)面的靜校正量以及浮動基準(zhǔn)面的靜校正量;折射剩余靜校正量確定裝置,用于根據(jù)所述浮動基準(zhǔn)面的靜校正量確定所述的地震資料對應(yīng)的折射剩余靜校正量;最終靜校正量確定裝置,用于根據(jù)所述的水平基準(zhǔn)面的靜校正量以及所述的折射剩余靜校正量確定所述地震資料的最終靜校正量,其中,所述的折射剩余靜校正量確定裝置包括:靜校正量應(yīng)用模塊,用于將所述浮動基準(zhǔn)面的靜校正量應(yīng)用到所述的折射波初至,得到應(yīng)用浮動基準(zhǔn)面的靜校正量的初至;穩(wěn)定折射層初至確定模塊,用于根據(jù)所述應(yīng)用浮動基準(zhǔn)面的靜校正量的初至確定穩(wěn)定折射層的初至;折射剩余靜校正量確定模塊,用于根據(jù)所述穩(wěn)定折射層的初至確定所述的地震資料對應(yīng)的折射剩余靜校正量,所述的折射剩余靜校正量包括炮點(diǎn)的折射剩余靜校正量以及檢波點(diǎn)的折射剩余靜校正量。
[0008]本發(fā)明的有益效果在于,提供了一種基于近地表層的折射波靜校正處理方法及系統(tǒng),涉及物探地震資料的近地表靜校正方案,以一次靜校正(模型法或?qū)游龇o校正)為基礎(chǔ)獲得浮動基準(zhǔn)面靜校正,把浮動基準(zhǔn)面靜校正應(yīng)用到拾取的折射波初至上,然后選擇穩(wěn)定的折射層并計(jì)算折射波速度,通過計(jì)算的折射波速度進(jìn)行線性動校正,計(jì)算折射剩余靜校正,折射剩余靜校正與水平基準(zhǔn)面一次靜校正的和為最終靜校正量,其能解決近地表層短波長靜校正問題,與其他靜校正方法(如:模型法或?qū)游龇?相結(jié)合最終解決靜校正問題。
[0009]為讓本發(fā)明的上述和其他目的、特征和優(yōu)點(diǎn)能更明顯易懂,下文特舉較佳實(shí)施例,并配合所附圖式,作詳細(xì)說明如下。
【專利附圖】
【附圖說明】
[0010]為了更清楚地說明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案,下面將對實(shí)施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作簡單地介紹,顯而易見地,下面描述中的附圖僅僅是本發(fā)明的一些實(shí)施例,對于本領(lǐng)域普通技術(shù)人員來講,在不付出創(chuàng)造性勞動的前提下,還可以根據(jù)這些附圖獲得其他的附圖。
[0011]圖1為本發(fā)明實(shí)施例提供的一種基于近地表層的折射波靜校正處理方法的流程圖;
[0012]圖2為圖1中的步驟S104的具體流程圖;
[0013]圖3為圖2中的步驟S202的具體流程圖;
[0014]圖4為圖2中的步驟S203的實(shí)施方式一的具體流程圖;
[0015]圖5為圖2中的步驟S203的實(shí)施方式二的具體流程圖;
[0016]圖6為圖1中的步驟S105的具體流程圖;
[0017]圖7為本發(fā)明實(shí)施例提供的一種基于近地表層的折射波靜校正處理系統(tǒng)的結(jié)構(gòu)框圖;
[0018]圖8為本發(fā)明實(shí)施例提供的一種基于近地表層的折射波靜校正處理系統(tǒng)中的折射剩余靜校正量確定裝置400的具體結(jié)構(gòu)框圖;
[0019]圖9為本發(fā)明實(shí)施例提供的一種基于近地表層的折射波靜校正處理系統(tǒng)中的穩(wěn)定折射層初至確定模塊402的具體結(jié)構(gòu)框圖;
[0020]圖10為本發(fā)明實(shí)施例提供的一種基于近地表層的折射波靜校正處理系統(tǒng)中的折射剩余靜校正量確定模塊403的實(shí)施方式一的具體結(jié)構(gòu)框圖;
[0021]圖11為本發(fā)明實(shí)施例提供的一種基于近地表層的折射波靜校正處理系統(tǒng)中的折射剩余靜校正量確定模塊403的實(shí)施方式二的具體結(jié)構(gòu)框圖;
[0022]圖12為本發(fā)明實(shí)施例提供的一種基于近地表層的折射波靜校正處理系統(tǒng)中的最終靜校正量確定裝置500的具體結(jié)構(gòu)框圖;
[0023]圖13為具體實(shí)施例中的地表高程示意圖;
[0024]圖14為具體實(shí)施例中的浮動基準(zhǔn)面(高速頂1500半徑平滑)示意圖;
[0025]圖15為具體實(shí)施例中的浮動基準(zhǔn)面靜校正量示意圖;
[0026]圖16為具體實(shí)施例中的水平基準(zhǔn)面靜校正量示意圖;
[0027]圖17為具體實(shí)施例中的折射剩余靜校正量示意圖;
[0028]圖18為具體實(shí)施例中的最終靜校正量示意圖;
[0029]圖19為具體實(shí)施例中的層析模型法靜校正水平疊加剖面示意圖;
[0030]圖20為具體實(shí)施例中的靜校正水平疊加剖面示意圖。
【具體實(shí)施方式】
[0031]下面將結(jié)合本發(fā)明實(shí)施例中的附圖,對本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述,顯然,所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例,而不是全部的實(shí)施例?;诒景l(fā)明中的實(shí)施例,本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其他實(shí)施例,都屬于本發(fā)明保護(hù)的范圍。
[0032]本發(fā)明的目的是提供一種解決近地表層短波長靜校正問題的方法,與其他靜校正方法(如:模型法或?qū)游龇?相結(jié)合最終解決靜校正問題。
[0033]圖1為本發(fā)明提出的一種基于近地表層的折射波靜校正處理方法的具體流程圖,由圖1可知,在實(shí)施方式一中,所述的方法包括:
[0034]SlOl:采集地震資料。
[0035]在具體的實(shí)施例中,可利用現(xiàn)有儀器采集得到地震資料。
[0036]S102:從所述的地震資料中拾取出折射波初至。
[0037]在具體的實(shí)施方式中,對地震資料進(jìn)行精細(xì)的折射波初至拾取,初至拾取可采用自動拾取或人機(jī)交互拾取。
[0038]S103:確定水平基準(zhǔn)面的靜校正量以及浮動基準(zhǔn)面的靜校正量。
[0039]在具體的實(shí)施方式中,可利用模型法或?qū)游龇?或者其它任何通過近地表速度模型計(jì)算靜校正量的方法)計(jì)算水平基準(zhǔn)面靜校正量和浮動基準(zhǔn)面的靜校正量,浮動基準(zhǔn)面為近地表模型的高速頂界面經(jīng)平滑獲得或平滑地表下移0-200米獲得,平滑半徑為l_5km。
[0040]如圖13所示為具體實(shí)施例中的地表高程示意圖,圖14為具體實(shí)施例中的浮動基準(zhǔn)面(高速頂1.5km半徑平滑)示意圖。在圖13、圖14所示的具體實(shí)施例中,采用層析法計(jì)算靜校正量,浮動基準(zhǔn)面為近地表層析模型的高速頂界面(2000m/s速度界面)經(jīng)平滑獲得,平滑半徑為1.5km,圖15為計(jì)算的層析浮動基準(zhǔn)面靜校正量,圖16為層析水平基準(zhǔn)面靜校正量。
[0041]S104:根據(jù)所述浮動基準(zhǔn)面的靜校正量確定所述的地震資料對應(yīng)的折射剩余靜校正量。圖2為步驟S104的具體流程圖。
[0042]S105:根據(jù)所述的水平基準(zhǔn)面的靜校正量以及所述的折射剩余靜校正量確定所述地震資料的最終靜校正量。圖6為步驟S105的具體流程圖。
[0043]圖2為步驟S104的具體流程圖,由圖2可知,該步驟具體包括:
[0044]S201:將所述浮動基準(zhǔn)面的靜校正量應(yīng)用到所述的折射波初至,得到應(yīng)用浮動基準(zhǔn)面的靜校正量的初至。即把浮動基準(zhǔn)面的靜校正量加上折射波的初至?xí)r間上,即可得到應(yīng)用浮動基準(zhǔn)面靜校正量的初至。
[0045]S202:根據(jù)所述應(yīng)用浮動基準(zhǔn)面的靜校正量的初至確定穩(wěn)定折射層的初至,圖3為步驟S202的具體流程圖,由圖3可知,該步驟具體包括:
[0046]S301:獲取共中心點(diǎn)面元的長以及寬,在具體的實(shí)施方式中,CMP面元的長和寬均為ι-?ο倍的道距。
[0047]S302:獲取控制點(diǎn)的間隔,在具體的實(shí)施方式中,控制點(diǎn)的間隔為l_5km。
[0048]S303:將所述的應(yīng)用浮動基準(zhǔn)面靜校正量的初至通過共中心點(diǎn)道集的方式進(jìn)行顯示,得到初至顯示;
[0049]S304:在所述的初至顯示上進(jìn)行交互選擇,確定出穩(wěn)定折射層;
[0050]S305:根據(jù)所述的穩(wěn)定折射層進(jìn)行插值獲得每個共中心點(diǎn)面元上初至的偏移距范圍,在具體的實(shí)施方式中,可以通過任何常用插值方法。
[0051]S306:根據(jù)所述的偏移距范圍輸出所述穩(wěn)定折射層的初至。
[0052]在圖13所示的具體的實(shí)施例中,利用應(yīng)用浮動基準(zhǔn)面靜校正量的初至,在共中心點(diǎn)道集上交互選擇全區(qū)相對穩(wěn)定折射層,共中心點(diǎn)面元的長和寬都為4倍道距,在控制點(diǎn)上以偏移距范圍的方式選擇,控制點(diǎn)的間隔為2km,通過控制點(diǎn)上的穩(wěn)定折射層線性插值獲得每個面元上初至的偏移距范圍,輸出該穩(wěn)定折射層的初至。
[0053]由圖2可知,步驟S104還包括:
[0054]S203:根據(jù)所述穩(wěn)定折射層的初至確定所述的地震資料對應(yīng)的折射剩余靜校正量,所述的折射剩余靜校正量包括炮點(diǎn)的折射剩余靜校正量以及檢波點(diǎn)的折射剩余靜校正量。
[0055]圖4為步驟S203的實(shí)施方式一的具體流程圖,由圖4可知,該步驟具體包括:
[0056]S401:根據(jù)所述的穩(wěn)定折射層的初至確定所述穩(wěn)定折射層的速度,在具體的實(shí)施方式中,對穩(wěn)定折射層的初至進(jìn)行擬合求取折射層的速度,或者通過其他方法求取折射層速度。
[0057]S402:利用所述穩(wěn)定折射層的速度對所述穩(wěn)定折射層的初至進(jìn)行線性動校正,得到線性動校正后的初至?xí)r間。在具體的實(shí)施方式中,利用穩(wěn)定折射層的速度對穩(wěn)定折射層的初至進(jìn)行線性動校正,得到線性動校正后的初至?xí)r間(LTij)。
[0058]S403:根據(jù)所述線性動校正后的初至?xí)r間構(gòu)建方程組,假設(shè)線性動校正后的初至?xí)r間(LTij)為對應(yīng)炮點(diǎn)的折射剩余靜校正量和檢波點(diǎn)的折射剩余靜校正量的和,所有初至?xí)r間構(gòu)建成大型方程組,如下所示:
[0059]Tsl+Trl = LT11
[0060]Tsl+Tr2 = LT12
[0061]
[0062]Ts2+Trl = LT21
[0063]Ts2+Tr2 = LT22
[0064]...
[0065]Tsi+TrJ = LTij
[0066]其中,Tsi表示第i個激發(fā)點(diǎn)的剩余靜校正量,Trj表示第j個接收點(diǎn)的剩余靜校正量,LTij表示第i個激發(fā)點(diǎn)與第j個接收點(diǎn)初至?xí)r間應(yīng)用浮動基準(zhǔn)面靜校正后再應(yīng)用線性動校正的時間。
[0067]S404:求解所述的方程組,得到折射剩余靜校正量,所述的折射剩余靜校正量包括炮點(diǎn)的折射剩余靜校正量、檢波點(diǎn)的折射剩余靜校正量。在【具體實(shí)施方式】中,求解大型方程組,通過高斯賽德爾迭代方法解方程組得到炮點(diǎn)的折射剩余靜校正量、檢波點(diǎn)的折射剩余靜校正量。
[0068]圖5為步驟S203的實(shí)施方式二的具體流程圖,由圖5可知,該步驟具體包括:
[0069]S501:根據(jù)所述的穩(wěn)定折射層的初至確定所述穩(wěn)定折射層的速度,在具體的實(shí)施方式中,對穩(wěn)定折射層的初至進(jìn)行擬合求取折射層的速度,或者通過其他方法求取折射層速度。
[0070]S502:獲取預(yù)先設(shè)定的第一平滑半徑,在具體的實(shí)施方式中,平滑半徑為l_5km。
[0071]S503:所述穩(wěn)定折射層的速度根據(jù)所述的平滑半徑在平面上平滑,得到平滑后的穩(wěn)定折射層的速度。在圖13所示的具體實(shí)施例中,可將獲得的穩(wěn)定折射層的初至擬合求取折射層的速度,或者通過其他方法求取折射層速度,并在平面上平滑,平滑半徑為2km。
[0072]S504:利用所述穩(wěn)定折射層的速度對所述穩(wěn)定折射層的初至進(jìn)行線性動校正,得到線性動校正后的初至?xí)r間。在具體的實(shí)施方式中,利用穩(wěn)定折射層的速度對穩(wěn)定折射層的初至進(jìn)行線性動校正,得到線性動校正后的初至?xí)r間(LTij)。
[0073]S505:根據(jù)所述線性動校正后的初至?xí)r間構(gòu)建方程組,假設(shè)線性動校正后的初至?xí)r間(LTij)為對應(yīng)炮點(diǎn)的折射剩余靜校正量和檢波點(diǎn)的折射剩余靜校正量的和,所有初至?xí)r間構(gòu)建成大型方程組。
[0074]S506:求解所述的方程組,得到折射剩余靜校正量,所述的折射剩余靜校正量包括炮點(diǎn)的折射剩余靜校正量、檢波點(diǎn)的折射剩余靜校正量。在【具體實(shí)施方式】中,求解大型方程組,通過高斯賽德爾迭代方法解方程組得到炮點(diǎn)的折射剩余靜校正量、檢波點(diǎn)的折射剩余靜校正量。
[0075]圖6為步驟S105的具體流程圖,由圖6可知,該步驟具體包括:
[0076]S601:獲取預(yù)先設(shè)定的第二平滑半徑,在具體的實(shí)施方式中,平滑半徑為l_5km。
[0077]S602:根據(jù)所述的第二平滑半徑對所述的折射剩余靜校正量進(jìn)行平滑處理,得到平滑處理后的折射剩余靜校正量;
[0078]S603:根據(jù)所述的折射剩余靜校正量以及平滑處理的折射剩余靜校正量確定去低頻后的折射剩余靜校正量。在具體的實(shí)施方式中,利用平滑前的折射剩余靜校正量減去平滑后的折射剩余靜校正量即得到去低頻后的折射剩余靜校正量。
[0079]S604:根據(jù)所述的水平基準(zhǔn)面的靜校正量以及去低頻后的折射剩余靜校正量確定所述地震資料的最終靜校正量。在具體的實(shí)施方式中,水平基準(zhǔn)面的靜校正量減去去低頻后的折射剩余靜校正量,即得到最終靜校正量。
[0080]如上所述,即為本發(fā)明提供的一種基于近地表層的折射波靜校正處理方法,對初至應(yīng)用一次靜校正在線性動校正并構(gòu)建方程組,通過高斯賽德爾迭代算法獲得最終折射剩余靜校正量,能夠有效解決復(fù)雜近地表地區(qū)長波長靜校正問題和短波長靜校正問題。
[0081]圖7為本發(fā)明實(shí)施例提供的一種基于近地表層的折射波靜校正處理系統(tǒng)的結(jié)構(gòu)框圖,由圖7可知,所述的系統(tǒng)包括:
[0082]地震資料采集裝置100,用于采集地震資料。
[0083]在具體的實(shí)施例中,可利用現(xiàn)有儀器采集得到地震資料。
[0084]折射波初至拾取裝置200,用于從所述的地震資料中拾取出折射波初至。
[0085]在具體的實(shí)施方式中,對地震資料進(jìn)行精細(xì)的折射波初至拾取,初至拾取可采用自動拾取或人機(jī)交互拾取。
[0086]基準(zhǔn)面靜校正量確定裝置300,用于確定水平基準(zhǔn)面的靜校正量以及浮動基準(zhǔn)面的靜校正量。
[0087]在具體的實(shí)施方式中,可利用模型法或?qū)游龇?或者其它任何通過近地表速度模型計(jì)算靜校正量的方法)計(jì)算水平基準(zhǔn)面靜校正量和浮動基準(zhǔn)面的靜校正量,浮動基準(zhǔn)面為近地表模型的高速頂界面經(jīng)平滑獲得或平滑地表下移0-200m獲得,平滑半徑為l-5km。
[0088]如圖13所示為具體實(shí)施例中的地表高程示意圖,圖14為具體實(shí)施例中的浮動基準(zhǔn)面(高速頂1.5km半徑平滑)示意圖。在圖13、圖14所示的具體實(shí)施例中,采用層析法計(jì)算靜校正量,浮動基準(zhǔn)面為近地表層析模型的高速頂界面(2000m/s速度界面)經(jīng)平滑獲得,平滑半徑為1.5km米,圖15為計(jì)算的層析浮動基準(zhǔn)面靜校正量,圖16為層析水平基準(zhǔn)面靜校正量。
[0089]折射剩余靜校正量確定裝置400,用于根據(jù)所述浮動基準(zhǔn)面的靜校正量確定所述的地震資料對應(yīng)的折射剩余靜校正量。圖8為折射剩余靜校正量確定裝置400的具體結(jié)構(gòu)框圖。
[0090]最終靜校正量確定裝置500,用于根據(jù)所述的水平基準(zhǔn)面的靜校正量以及所述的折射剩余靜校正量確定所述地震資料的最終靜校正量。圖12為最終靜校正量確定裝置500的具體結(jié)構(gòu)框圖。
[0091]圖8為本發(fā)明實(shí)施例提供的一種基于近地表層的折射波靜校正處理系統(tǒng)中的折射剩余靜校正量確定裝置400的具體結(jié)構(gòu)框圖,由圖8可知,折射剩余靜校正量確定裝置400具體包括:
[0092]靜校正量應(yīng)用模塊401,用于將所述浮動基準(zhǔn)面的靜校正量應(yīng)用到所述的折射波初至,得到應(yīng)用浮動基準(zhǔn)面的靜校正量的初至。即把浮動基準(zhǔn)面的靜校正量加上折射波的初至?xí)r間上,即可得到應(yīng)用浮動基準(zhǔn)面靜校正量的初至。
[0093]穩(wěn)定折射層初至確定模塊402,用于根據(jù)所述應(yīng)用浮動基準(zhǔn)面的靜校正量的初至確定穩(wěn)定折射層的初至,圖9為穩(wěn)定折射層初至確定模塊402的具體結(jié)構(gòu)框圖,由圖9可知,穩(wěn)定折射層初至確定模塊402具體包括:
[0094]長寬獲取單元4021,用于獲取共中心點(diǎn)面元的長以及寬,在具體的實(shí)施方式中,CMP面元的長和寬均為1-10倍的道距。
[0095]控制點(diǎn)間隔獲取單元4022,用于獲取控制點(diǎn)的間隔,在具體的實(shí)施方式中,控制點(diǎn)的間隔為I_5km。
[0096]初至顯示單元4023,用于將所述的應(yīng)用浮動基準(zhǔn)面靜校正量的初至通過共中心點(diǎn)道集的方式進(jìn)行顯示,得到初至顯示;
[0097]穩(wěn)定折射層確定單元4024,用于在所述的初至顯示上進(jìn)行交互選擇,確定出穩(wěn)定折射層;
[0098]插值單元4025,用于根據(jù)所述的穩(wěn)定折射層進(jìn)行插值獲得每個共中心點(diǎn)面元上初至的偏移距范圍,在具體的實(shí)施方式中,可以通過任何常用插值方法。
[0099]初至確定單元4026,用于根據(jù)所述的偏移距范圍輸出所述穩(wěn)定折射層的初至。
[0100]在圖13所示的具體的實(shí)施例中,利用應(yīng)用浮動基準(zhǔn)面靜校正量的初至,在共中心點(diǎn)道集上交互選擇全區(qū)相對穩(wěn)定折射層,共中心點(diǎn)面元的長和寬都為4倍道距,在控制點(diǎn)上以偏移距范圍的方式選擇,控制點(diǎn)的間隔為2km,通過控制點(diǎn)上的穩(wěn)定折射層線性插值獲得每個面元上初至的偏移距范圍,輸出該穩(wěn)定折射層的初至。
[0101]由圖8可知,折射剩余靜校正量確定裝置還包括:
[0102]折射剩余靜校正量確定模塊403,用于根據(jù)所述穩(wěn)定折射層的初至確定所述的地震資料對應(yīng)的折射剩余靜校正量,所述的折射剩余靜校正量包括炮點(diǎn)的折射剩余靜校正量以及檢波點(diǎn)的折射剩余靜校正量。
[0103]圖10為折射剩余靜校正量確定模塊403的實(shí)施方式一的具體結(jié)構(gòu)框圖,由圖10可知,折射剩余靜校正量確定模塊403具體包括:
[0104]穩(wěn)定折射層速度確定單元4031,用于根據(jù)所述的穩(wěn)定折射層的初至確定所述穩(wěn)定折射層的速度,在具體的實(shí)施方式中,對穩(wěn)定折射層的初至進(jìn)行擬合求取折射層的速度,或者通過其他方法求取折射層速度。
[0105]線性動校正單元4032,用于利用所述穩(wěn)定折射層的速度對所述穩(wěn)定折射層的初至進(jìn)行線性動校正,得到線性動校正后的初至?xí)r間。在具體的實(shí)施方式中,利用穩(wěn)定折射層的速度對穩(wěn)定折射層的初至進(jìn)行線性動校正,得到線性動校正后的初至?xí)r間(LTij)。
[0106]方程組構(gòu)建單元4033,用于根據(jù)所述線性動校正后的初至?xí)r間構(gòu)建方程組,假設(shè)線性動校正后的初至?xí)r間(LTij)為對應(yīng)炮點(diǎn)的折射剩余靜校正量和檢波點(diǎn)的折射剩余靜校正量的和,所有初至?xí)r間構(gòu)建成大型方程組,如下所示:
[0107]Tsl+Trl = LT11
[0108]Tsl+Tr2 = LT12
[0109]...
[0110]Ts2+Trl = LT21
[0111]Ts2+Tr2 = LT22
[0112]…
[0113]Tsi+TrJ = LTij
[0114]其中,Tsi表示第i個激發(fā)點(diǎn)的剩余靜校正量,IVj表示第j個接收點(diǎn)的剩余靜校正量,LTij表示第i個激發(fā)點(diǎn)與第j個接收點(diǎn)初至?xí)r間應(yīng)用浮動基準(zhǔn)面靜校正后再應(yīng)用線性動校正的時間。
[0115]方程組求解單元4034,用于求解所述的方程組,得到折射剩余靜校正量,所述的折射剩余靜校正量包括炮點(diǎn)的折射剩余靜校正量、檢波點(diǎn)的折射剩余靜校正量。在【具體實(shí)施方式】中,求解大型方程組,通過高斯賽德爾迭代方法解方程組得到炮點(diǎn)的折射剩余靜校正量、檢波點(diǎn)的折射剩余靜校正量。
[0116]圖11為折射剩余靜校正量確定模塊403的實(shí)施方式二的具體結(jié)構(gòu)框圖,由圖11可知,折射剩余靜校正量確定模塊403還包括:
[0117]第一平滑半徑獲取單元4035,用于獲取預(yù)先設(shè)定的第一平滑半徑,在具體的實(shí)施方式中,平滑半徑為l_5km。
[0118]平滑處理單元4036,用于所述穩(wěn)定折射層的速度根據(jù)所述的平滑半徑在平面上平滑,得到平滑后的穩(wěn)定折射層的速度。在圖13所示的具體實(shí)施例中,可將獲得的穩(wěn)定折射層的初至擬合求取折射層的速度,或者通過其他方法求取折射層速度,并在平面上平滑,平滑半徑為2km。
[0119]圖12為最終靜校正量確定裝置500的具體結(jié)構(gòu)框圖,由圖12可知,最終靜校正量確定裝置500具體包括:
[0120]第二平滑半徑獲取模塊501,用于獲取預(yù)先設(shè)定的第二平滑半徑,在具體的實(shí)施方式中,平滑半徑為l_5km。
[0121]第二平滑半徑獲取模塊502,用于根據(jù)所述的第二平滑半徑對所述的折射剩余靜校正量進(jìn)行平滑處理,得到平滑處理后的折射剩余靜校正量;
[0122]去低頻處理模塊503,用于根據(jù)所述的折射剩余靜校正量以及平滑處理的折射剩余靜校正量確定去低頻后的折射剩余靜校正量。在具體的實(shí)施方式中,利用平滑前的折射剩余靜校正量減去平滑后的折射剩余靜校正量即得到去低頻后的折射剩余靜校正量。
[0123]最終靜校正量確定模塊504,用于根據(jù)所述的水平基準(zhǔn)面的靜校正量以及去低頻后的折射剩余靜校正量確定所述地震資料的最終靜校正量。在具體的實(shí)施方式中,水平基準(zhǔn)面的靜校正量減去去低頻后的折射剩余靜校正量,即得到最終靜校正量。
[0124]如上所述,即為本發(fā)明提供的一種基于近地表層的折射波靜校正處理系統(tǒng),對初至應(yīng)用一次靜校正在線性動校正并構(gòu)建方程組,通過高斯賽德爾迭代算法獲得最終折射剩余靜校正量,能夠有效解決復(fù)雜近地表地區(qū)長波長靜校正問題和短波長靜校正問題。
[0125]下面結(jié)合具體的實(shí)施例,詳細(xì)介紹本發(fā)明的技術(shù)方案。
[0126]圖13為該具體實(shí)施例中的地表高程示意圖,圖14為具體實(shí)施例中的浮動基準(zhǔn)面(高速頂1.5km半徑平滑)示意圖。本方案主要包括:
[0127]I)采集地震資料并進(jìn)行初至拾??;
[0128]對地震資料進(jìn)行精細(xì)的折射波初至拾取,初至拾取可采用自動拾取或人機(jī)交互拾取。
[0129]2)計(jì)算一次靜校正量
[0130]利用模型法或?qū)游龇?或者其它任何通過近地表速度模型計(jì)算靜校正量的方法)計(jì)算水平基準(zhǔn)面靜校正量和浮動基準(zhǔn)面的靜校正量,本實(shí)施例中是采用層析法計(jì)算靜校正量,浮動基準(zhǔn)面為近地表層析模型的高速頂界面(2000m/s速度界面)經(jīng)平滑獲得,平滑半徑為1.5km,圖15為計(jì)算的層析浮動基準(zhǔn)面靜校正量,圖16為層析水平基準(zhǔn)面靜校正量。
[0131]3)對初至應(yīng)用浮動基準(zhǔn)面靜校正量
[0132]步驟3)所述的對初至應(yīng)用浮動基準(zhǔn)面靜校正量是把浮動基準(zhǔn)面靜校正量應(yīng)用到初至?xí)r間上,得到應(yīng)用浮動基準(zhǔn)面靜校正量的初至。
[0133]4)初至折射層選擇
[0134]利用應(yīng)用浮動基準(zhǔn)面靜校正量的初至,在CMP道集上交互選擇全區(qū)相對穩(wěn)定的折射層,CMP面元為長和寬都為4倍道距,在控制點(diǎn)上以偏移距范圍的方式選擇,控制點(diǎn)間隔為2km,通過控制點(diǎn)上的偏移距范圍線性插值獲得每個面元上初至的偏移距范圍,輸出該穩(wěn)定折射層的初至。
[0135]5)折射速度分析;
[0136]對獲得的穩(wěn)定折射層的初至擬合求取折射層速度(或者通過其他方法求取折射層速度)并在平面上平滑,平滑半徑為2km。
[0137]6)線性動校正及方程組構(gòu)建
[0138]利用獲得的平滑后的折射速度對穩(wěn)定折射層的初至進(jìn)行線性動校正,得到線性動校正后的初至?xí)r間(LTij),假設(shè)該時間為對應(yīng)炮點(diǎn)和檢波點(diǎn)折射剩余靜校正量的和,所有初至?xí)r間構(gòu)建成大型方程組。
[0139]7)折射剩余靜校正計(jì)算及低頻去除
[0140]求解大型方程組,得到炮、檢點(diǎn)折射剩余靜校正量。對獲得的炮、檢點(diǎn)折射剩余靜校正量進(jìn)行平滑,平滑半徑為3km,然后利用平滑前的量減去平滑后的量得到去低頻后的折射剩余靜校正量,如圖17所示。
[0141]8)最終靜校正計(jì)算
[0142]用步驟2)獲得的水平基準(zhǔn)面靜校正量減去步驟7)獲得的去低頻后的折射剩余靜校正量得到最終的靜校正量,如圖18所示。
[0143]圖19為具體實(shí)施例中的層析模型法靜校正水平疊加剖面示意圖,圖20為具體實(shí)施例中的靜校正水平疊加剖面示意圖。由圖19、20的最終應(yīng)用效果對比可知,應(yīng)用本方法的疊加剖面效果比常規(guī)層析法靜校正效果要好。
[0144]綜上所述,本發(fā)明提出的一種基于近地表層的折射波靜校正處理方法及系統(tǒng),涉及物探地震數(shù)據(jù)的近地表靜校正,以一次靜校正(模型法或?qū)游龇o校正)為基礎(chǔ),對拾取的初至應(yīng)用一次靜校正到浮動基準(zhǔn)面上,然后選擇穩(wěn)定的折射層并計(jì)算折射波速度,通過計(jì)算的折射波速度進(jìn)行線性動校正,計(jì)算折射剩余靜校正,折射剩余靜校正與水平基準(zhǔn)面一次靜校正的和為最終靜校正量,能夠解決短波長靜校正問題。
[0145]本領(lǐng)域普通技術(shù)人員可以理解實(shí)現(xiàn)上述實(shí)施例方法中的全部或部分流程,可以通過計(jì)算機(jī)程序來指令相關(guān)的硬件來完成,所述的程序可存儲于一般計(jì)算機(jī)可讀取存儲介質(zhì)中,該程序在執(zhí)行時,可包括如上述各方法的實(shí)施例的流程。其中,所述的存儲介質(zhì)可為磁碟、光盤、只讀存儲記憶體(Read-Only Memory, ROM)或隨機(jī)存儲記憶體(Random AccessMemory, RAM)等。
[0146]本領(lǐng)域技術(shù)人員還可以了解到本發(fā)明實(shí)施例列出的各種功能是通過硬件還是軟件來實(shí)現(xiàn)取決于特定的應(yīng)用和整個系統(tǒng)的設(shè)計(jì)要求。本領(lǐng)域技術(shù)人員可以對于每種特定的應(yīng)用,可以使用各種方法實(shí)現(xiàn)所述的功能,但這種實(shí)現(xiàn)不應(yīng)被理解為超出本發(fā)明實(shí)施例保護(hù)的范圍。
[0147]本發(fā)明中應(yīng)用了具體實(shí)施例對本發(fā)明的原理及實(shí)施方式進(jìn)行了闡述,以上實(shí)施例的說明只是用于幫助理解本發(fā)明的方法及其核心思想;同時,對于本領(lǐng)域的一般技術(shù)人員,依據(jù)本發(fā)明的思想,在【具體實(shí)施方式】及應(yīng)用范圍上均會有改變之處,綜上所述,本說明書內(nèi)容不應(yīng)理解為對本發(fā)明的限制。
【權(quán)利要求】
1.一種基于近地表層的折射波靜校正處理方法,其特征是,所述的方法包括: 采集地震資料; 從所述的地震資料中拾取出折射波初至; 確定水平基準(zhǔn)面的靜校正量以及浮動基準(zhǔn)面的靜校正量; 根據(jù)所述浮動基準(zhǔn)面的靜校正量確定所述的地震資料對應(yīng)的折射剩余靜校正量;根據(jù)所述的水平基準(zhǔn)面的靜校正量以及所述的折射剩余靜校正量確定所述地震資料的最終靜校正量; 其中根據(jù)所述浮動基準(zhǔn)面的靜校正量確定所述的地震資料對應(yīng)的折射剩余靜校正量包括:將所述浮動基準(zhǔn)面的靜校正量應(yīng)用到所述的折射波初至,得到應(yīng)用浮動基準(zhǔn)面的靜校正量的初至;根據(jù)所述應(yīng)用浮動基準(zhǔn)面的靜校正量的初至確定穩(wěn)定折射層的初至;根據(jù)所述穩(wěn)定折射層的初至確定所述的地震資料對應(yīng)的折射剩余靜校正量,所述的折射剩余靜校正量包括炮點(diǎn)的折射剩余靜校正量以及檢波點(diǎn)的折射剩余靜校正量。
2.根據(jù)權(quán)利要求1所述的方法,其特征是,根據(jù)所述應(yīng)用浮動基準(zhǔn)面的靜校正量的初至確定穩(wěn)定折射層的初至包括: 獲取共中心點(diǎn)面元的長以及寬; 獲取控制點(diǎn)的間隔; 將所述的應(yīng)用浮動基準(zhǔn)面靜校正量的初至通過共中心點(diǎn)道集的方式進(jìn)行顯示,得到初至顯示; 在所述的初至顯示上進(jìn)行交互選擇,確定出穩(wěn)定折射層; 根據(jù)所述的穩(wěn)定折射層進(jìn)行插值獲得每個共中心點(diǎn)面元上初至的偏移距范圍; 根據(jù)所述的偏移距范圍輸出所述穩(wěn)定折射層的初至。
3.根據(jù)權(quán)利要求1所述的方法,其特征是,根據(jù)所述穩(wěn)定折射層的初至確定所述的地震資料對應(yīng)的折射剩余靜校正量包括: 根據(jù)所述的穩(wěn)定折射層的初至確定所述穩(wěn)定折射層的速度; 利用所述穩(wěn)定折射層的速度對所述穩(wěn)定折射層的初至進(jìn)行線性動校正,得到線性動校正后的初至?xí)r間; 根據(jù)所述線性動校正后的初至?xí)r間構(gòu)建方程組; 求解所述的方程組,得到折射剩余靜校正量,所述的折射剩余靜校正量包括炮點(diǎn)的折射剩余靜校正量、檢波點(diǎn)的折射剩余靜校正量。
4.根據(jù)權(quán)利要求3所述的方法,其特征是,在確定出所述穩(wěn)定折射層的速度之后,根據(jù)所述穩(wěn)定折射層的初至確定所述的地震資料對應(yīng)的折射剩余靜校正量還包括: 獲取預(yù)先設(shè)定的第一平滑半徑; 所述穩(wěn)定折射層的速度根據(jù)所述的平滑半徑在平面上平滑,得到平滑后的穩(wěn)定折射層的速度。
5.根據(jù)權(quán)利要求3或4所述的方法,其特征是,根據(jù)所述的水平基準(zhǔn)面的靜校正量以及所述的折射剩余靜校正量確定所述地震資料的最終靜校正量包括: 獲取預(yù)先設(shè)定的第二平滑半徑; 根據(jù)所述的第二平滑半徑對所述的折射剩余靜校正量進(jìn)行平滑處理,得到平滑處理后的折射剩余靜校正量; 根據(jù)所述的折射剩余靜校正量以及平滑處理的折射剩余靜校正量確定去低頻后的折射剩余靜校正量; 根據(jù)所述的水平基準(zhǔn)面的靜校正量以及去低頻后的折射剩余靜校正量確定所述地震資料的最終靜校正量。
6.一種基于近地表層的折射波靜校正處理系統(tǒng),其特征是,所述的系統(tǒng)包括: 地震資料采集裝置,用于采集地震資料; 折射波初至拾取裝置,用于從所述的地震資料中拾取出折射波初至; 基準(zhǔn)面靜校正量確定裝置,用于確定水平基準(zhǔn)面的靜校正量以及浮動基準(zhǔn)面的靜校正量; 折射剩余靜校正量確定裝置,用于根據(jù)所述浮動基準(zhǔn)面的靜校正量確定所述的地震資料對應(yīng)的折射剩余靜校正量; 最終靜校正量確定裝置,用于根據(jù)所述的水平基準(zhǔn)面的靜校正量以及所述的折射剩余靜校正量確定所述地震資料的最終靜校正量; 其中,所述的折射剩余靜校正量確定裝置包括:靜校正量應(yīng)用模塊,用于將所述浮動基準(zhǔn)面的靜校正量應(yīng)用到所述的折射波初至,得到應(yīng)用浮動基準(zhǔn)面的靜校正量的初至;穩(wěn)定折射層初至確定模塊,用于根據(jù)所述應(yīng)用浮動基準(zhǔn)面的靜校正量的初至確定穩(wěn)定折射層的初至;折射剩余靜校正量確定模塊,用于根據(jù)所述穩(wěn)定折射層的初至確定所述的地震資料對應(yīng)的折射剩余靜校正量,所述的折射剩余靜校正量包括炮點(diǎn)的折射剩余靜校正量以及檢波點(diǎn)的折射剩余靜校正量。
7.根據(jù)權(quán)利要求6所述的系統(tǒng),其特征是,所述的穩(wěn)定折射層初至確定模塊包括: 長寬獲取單元,用于獲取共中心點(diǎn)面元的長以及寬; 控制點(diǎn)間隔獲取單元,用于獲取控制點(diǎn)的間隔; 初至顯示單元,用于將所述的應(yīng)用浮動基準(zhǔn)面靜校正量的初至通過共中心點(diǎn)道集的方式進(jìn)行顯示,得到初至顯示; 穩(wěn)定折射層確定單元,用于在所述的初至顯示上進(jìn)行交互選擇,確定出穩(wěn)定折射層; 插值單元,用于根據(jù)所述的穩(wěn)定折射層進(jìn)行插值獲得每個共中心點(diǎn)面元上初至的偏移距范圍; 初至確定單元,用于根據(jù)所述的偏移距范圍輸出所述穩(wěn)定折射層的初至。
8.根據(jù)權(quán)利要求6所述的系統(tǒng),其特征是,所述的折射剩余靜校正量確定模塊包括: 穩(wěn)定折射層速度確定單元,用于根據(jù)所述的穩(wěn)定折射層的初至確定所述穩(wěn)定折射層的速度; 線性動校正單元,用于利用所述穩(wěn)定折射層的速度對所述穩(wěn)定折射層的初至進(jìn)行線性動校正,得到線性動校正后的初至?xí)r間; 方程組構(gòu)建單元,用于根據(jù)所述線性動校正后的初至?xí)r間構(gòu)建方程組; 方程組求解單元,用于求解所述的方程組,得到折射剩余靜校正量,所述的折射剩余靜校正量包括炮點(diǎn)的折射剩余靜校正量、檢波點(diǎn)的折射剩余靜校正量。
9.根據(jù)權(quán)利要求8所述的系統(tǒng),其特征是,所述的折射剩余靜校正量確定模塊還包括: 第一平滑半徑獲取單元,用于獲取預(yù)先設(shè)定的第一平滑半徑; 平滑處理單元,用于所述穩(wěn)定折射層的速度根據(jù)所述的平滑半徑在平面上平滑,得到平滑后的穩(wěn)定折射層的速度。
10.根據(jù)權(quán)利要求8或9所述的系統(tǒng),其特征是,最終靜校正量確定裝置包括: 第二平滑半徑獲取模塊,用于獲取預(yù)先設(shè)定的第二平滑半徑; 平滑處理模塊,用于根據(jù)所述的第二平滑半徑對所述的折射剩余靜校正量進(jìn)行平滑處理,得到平滑處理后的折射剩余靜校正量; 去低頻處理模塊,用于根據(jù)所述的折射剩余靜校正量以及平滑處理的折射剩余靜校正量確定去低頻后的折射剩余靜校正量; 最終靜校正量確定模塊,用于根據(jù)所述的水平基準(zhǔn)面的靜校正量以及去低頻后的折射剩余靜校正量確定所述地震資料的最終靜校正量。
【文檔編號】G01V1/36GK104181599SQ201410408931
【公開日】2014年12月3日 申請日期:2014年8月19日 優(yōu)先權(quán)日:2014年8月19日
【發(fā)明者】祖云飛, 倪宇東, 呂景峰, 馬青坡, 韓淑敏, 張?jiān)虑? 李海翔 申請人:中國石油天然氣集團(tuán)公司, 中國石油集團(tuán)東方地球物理勘探有限責(zé)任公司