專利名稱：一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法
技術(shù)領(lǐng)域：
本發(fā)明涉及地球物理勘探中的地震疊前道集處理方法，尤其涉及一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法。
背景技術(shù)：
地震勘探是利用地下介質(zhì)的彈性差異，通過觀測(cè)和分析大地對(duì)人工激發(fā)地震波的響應(yīng)，推斷地下巖層的性質(zhì)和形態(tài)的地球物理勘探方法。地震勘探是鉆探前勘測(cè)石油與天然氣資源的重要手段，在煤田和工程地質(zhì)勘查、區(qū)域地質(zhì)研究和地殼研究等方面 ，也得到廣泛應(yīng)用。地震勘探主要分為地震資料采集、處理與解釋三個(gè)階段。在地震資料采集階段，多次覆蓋已成為最基本的野外采集方法，所謂的多次覆蓋是指，在地表適當(dāng)?shù)剡x取多個(gè)對(duì)地震波激發(fā)點(diǎn)和接收點(diǎn)，進(jìn)而得到來自地下同一反射點(diǎn)的信號(hào)，將這些信號(hào)組合在一起就構(gòu)成了一個(gè)共反射點(diǎn)道集(一種地震疊前道集)；在地震資料處理階段，以同相軸拉平后的共反射點(diǎn)道集為基礎(chǔ)，對(duì)其進(jìn)行疊加求和(共反射點(diǎn)疊加)可以提高地震資料的信噪比；在地震資料解釋階段，以同相軸拉平后的共反射點(diǎn)道集為基礎(chǔ)，通過分析振幅隨入射角的變化(AVA分析)可以得到地下的各項(xiàng)同性彈性參數(shù)信息，通過分析振幅隨方位角的變化(AVAz分析)可以得到地下的各項(xiàng)異性彈性參數(shù)信息。無論是地震資料處理階段的共反射點(diǎn)疊加，還是地震資料解釋階段的AVA和AVAz分析，都是以同相軸拉平后的共反射點(diǎn)道集為基礎(chǔ)的，同相軸不平會(huì)影響后續(xù)工作的質(zhì)量，甚至得到錯(cuò)誤的結(jié)果。處理上，常根據(jù)共反射點(diǎn)道集的時(shí)距曲線方程進(jìn)行同相軸拉平(動(dòng)校正)，雖然動(dòng)校正可以保證共反射點(diǎn)道集同相軸基本拉平，但受各種因素(如由于地表高程起伏劇烈導(dǎo)致的靜校正不準(zhǔn)確、由于水平層狀各項(xiàng)同性介質(zhì)假設(shè)條件導(dǎo)致的常規(guī)時(shí)距曲線方程誤差)的影響，動(dòng)校正往往無法實(shí)現(xiàn)共反射點(diǎn)道集同相軸精細(xì)拉平，繼而，降低了共反射點(diǎn)疊加的信噪比，并導(dǎo)致由AVA和AVAz分析得到的彈性參數(shù)不準(zhǔn)確甚至錯(cuò)誤。

發(fā)明內(nèi)容
為了解決現(xiàn)有技術(shù)中同相軸不平的問題，從而造成地震勘探結(jié)果不準(zhǔn)確的問題，提供了一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法。本發(fā)明實(shí)施例提供的一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法，包括，步驟101，對(duì)采集的地震數(shù)據(jù)進(jìn)行前期處理；步驟102，抽取想要進(jìn)行同相軸精細(xì)拉平處理的疊前道集，用I行、J列的二維數(shù)組D來表示，其中，I為每個(gè)地震道的采樣點(diǎn)數(shù)，J為疊前道集的地震道數(shù)，二維數(shù)組D，其第i行、第j列的值D(i，j)表示疊前道集的第j個(gè)地震道的第i個(gè)采樣點(diǎn)的振幅值；步驟103，設(shè)定同相軸精細(xì)拉平處理參數(shù)^、^!^、！^，其中^述的心為波形對(duì)比窗口的大小，Nw為整數(shù)；所述的Nm為波形對(duì)比窗口的移動(dòng)量，Nm為整數(shù)；所述的Ns為波形對(duì)比的搜索半徑，凡為整數(shù)；所述的TolgS群容差，是判定是否在某個(gè)波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)的依據(jù)；所述的Toli為個(gè)體容差，是優(yōu)選某個(gè)地震道的相似地震道的依據(jù)；步驟104，以Nw為窗口大小，以Nm為窗口移動(dòng)量，逐窗口判定是否在某個(gè)波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)，并計(jì)算拉平種子點(diǎn)移動(dòng)量；步驟105，計(jì)算拉伸采樣坐標(biāo)；步驟106，根據(jù)拉伸采樣坐標(biāo)實(shí)施同相軸精細(xì)拉平。根據(jù)所述一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法的一個(gè)進(jìn)一步的方面，對(duì)采集的地震數(shù)據(jù)進(jìn)行前期處理中包括，對(duì)采集的地震數(shù)據(jù)進(jìn)行去噪、靜校正、動(dòng)校正的前期處理。根據(jù)所述一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法的再一個(gè)進(jìn)一步的方面，所述疊前道集包括共反射點(diǎn)道集和方位角道集?！じ鶕?jù)所述一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法的另一個(gè)進(jìn)一步的方面，所述逐窗口判定是否在某個(gè)波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)，并計(jì)算拉平種子點(diǎn)移動(dòng)量具體包括，步驟401，對(duì)于第k個(gè)波形對(duì)比窗口，計(jì)算地震道與地震道之間最優(yōu)移動(dòng)量的二維數(shù)組S，其中二維數(shù)組S中的第J1行、第j2列的值S(ji，j2)表示疊前道集的第J1個(gè)地震道和第j2個(gè)地震道的最優(yōu)移動(dòng)量，并計(jì)算地震道與地震道之間最大相關(guān)系數(shù)的二維數(shù)組C，其中二維數(shù)組C中的第J1行、第j2列的值CU1, j2)表示疊前道集的第J1個(gè)地震道和第j2個(gè)地震道的最大相關(guān)系數(shù)；所述的地震道與地震道之間最優(yōu)移動(dòng)量SU1, j2)和地震道與地震道之間最大相關(guān)系數(shù)CU1, j2)的計(jì)算方法為設(shè)^= (k-l)Nm，L = is+Nw，定義集合/ = {l,2,---,i￥w+l{· ^ =,定義長(zhǎng)度為(Nw+1)的向量 f，f 的第i個(gè)元素f(i) =D(is+i-l，定義2NS+1個(gè)長(zhǎng)度為(Nw+1)的向量gl / G Z , gl的第i個(gè)元素玢⑴=D(is+l+i-l，則最大相關(guān)系數(shù)CU1, j2)由以下公式給出
f · gC(j\,j2) =r—r r^-=
ei Vf-fVg,-g;
tools]最優(yōu)移動(dòng)量sa，j2)，即取得ca，j2)時(shí)的I值，由以下公式給出
f . CsSCj1 J2) = argm^x ^——'
leL vf-fVg,-g,步驟402，對(duì)于第k個(gè)波形對(duì)比窗口，判定是否在該波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)，其中，所述的判定是否設(shè)置拉平種子點(diǎn)的方法為求取二維數(shù)組C的所有元素的平均值Cm Cm=^YYjC(JlJ2)
J /,-I Λ=>若cm彡Tolg,則需要在該波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)，記錄該波形對(duì)
f j +i \(N \
比窗口的中心點(diǎn)4^ = round] (k-l)Nm+^-,其中函數(shù)round(x)表示對(duì)x進(jìn)
V2 J
行四舍五入取整，轉(zhuǎn)步驟103，若Cm < Tolg,則不需要在該波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)，記k = k+Ι，轉(zhuǎn)步驟101繼續(xù)；
步驟403，對(duì)于第k個(gè)波形對(duì)比窗口，計(jì)算拉平種子點(diǎn)移動(dòng)量，定義長(zhǎng)度為J的向量m用來存儲(chǔ)每個(gè)地震道的拉平種子點(diǎn)移動(dòng)量，其中，所述的計(jì)算拉平種子點(diǎn)移動(dòng)量的方法為選取基準(zhǔn)地震道，定義集合i = <1,2，···,·/丨，對(duì)于任意一個(gè)定值J1道的相似道集合為之=“ |cc/,j2；^ MJ2 4，則基準(zhǔn)地震道為A =MggfcarCi(Jii)，其中函數(shù)card(人)表示集合Jii的元素的個(gè)數(shù)；優(yōu)選基準(zhǔn)地震道的相似地震道，相似地震道集合為I、，
Jh = {λ |[(λ，λ)乏 To},j2e j]；計(jì)算每個(gè)地震道的拉平種子點(diǎn)移動(dòng)量，令Jmn= mjn ./ ' = i^x J ’貝Ij
^ S(JbJl-SUbJnmi)
m J s(/ft，jm) - s{jb，Jima) , a) > Jmax。
由m(/)，j G Ja線性內(nèi)插求出，當(dāng)/ e其它根據(jù)所述一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法的另一個(gè)進(jìn)一步的方面，在計(jì)算拉伸采樣坐標(biāo)中，定義I行、J列的二維數(shù)組Xm用來存儲(chǔ)拉伸采樣坐標(biāo)，所述的計(jì)算拉伸采樣坐標(biāo)的方法為假設(shè)需要在Q個(gè)波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)，設(shè)其中第q個(gè)中心點(diǎn)以及它所對(duì)應(yīng)的拉平種子點(diǎn)移動(dòng)量分別為i,和q = 1，2，…，Q，定
義集合及υ，···， ·ρ} >Λ+ =RiJ {1,1},則對(duì)于任意一個(gè)定值J j eJ
^/ + Hli(J),''I/ G R
I,當(dāng)f = I
X-(U)=f、卜，。
I, -=\ι = Ι
由Χ_ (/，./),/ e r線性內(nèi)插求出，當(dāng)i e其它根據(jù)所述一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法的另一個(gè)進(jìn)一步的方面，根據(jù)拉伸采樣坐標(biāo)實(shí)施同相軸精細(xì)拉平進(jìn)一步包括，定義I行、J列的二維數(shù)組X用來存儲(chǔ)原始采樣坐標(biāo)，定義I行、J列的二維數(shù)組Dn 用來存儲(chǔ)同相軸精細(xì)拉平處理后的疊前道集；設(shè)置原始采樣坐標(biāo),具體地，X(i, j) = i, i = 1，2，…，I; j = I, 2，...，J ;對(duì)于任意一個(gè)定值j，_/ e J，以X的第j列Xj為自變量、以D的第j列Dj為函數(shù)值，構(gòu)建三次樣條插值函數(shù)7=/x Dj ω，則d_'(u)=/Xj,Dj (XnevVj))。根據(jù)所述一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法的另一個(gè)進(jìn)一步的方面，在根據(jù)拉伸采樣坐標(biāo)實(shí)施同相軸精細(xì)拉平之后還包括步驟107，用Dnew代替步驟102中的D，并在步驟103中設(shè)定新的同相軸精細(xì)拉平處理參數(shù)，再次執(zhí)行步驟104-步驟106。根據(jù)所述一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法的另一個(gè)進(jìn)一步的方面，重復(fù)執(zhí)行所述步驟107兩至三次。根據(jù)所述一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法的另一個(gè)進(jìn)一步的方面，所述在步驟103中設(shè)定新的同相軸精細(xì)拉平處理參數(shù)進(jìn)一步包括，按照參數(shù)設(shè)定的先后次序Nw、Nm、Ns逐漸減小，Tolg, Toli保持不變或增大。通過本發(fā)明實(shí)施例的方法，該方法對(duì)于地震資料的處理和解釋都具有重要的意義，并具有穩(wěn)定、高效的特點(diǎn)。


此處所說明的附圖用來提供對(duì)本發(fā)明的進(jìn)一步理解，構(gòu)成本申請(qǐng)的一部分，并不構(gòu)成對(duì)本發(fā)明的限定。在附圖中 圖I是本發(fā)明實(shí)施例提供的一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法流程圖；圖2是本發(fā)明實(shí)施例提供的一個(gè)方位角道集；圖3是本發(fā)明實(shí)施例提供的移動(dòng)波形對(duì)比窗口的示意圖；圖4是本發(fā)明實(shí)施例提供的用來存儲(chǔ)地震道與地震道之間最大相關(guān)系數(shù)的二維數(shù)組C的示意圖；圖5是本發(fā)明實(shí)施例提供的用來存儲(chǔ)地震道與地震道之間最優(yōu)移動(dòng)量的二維數(shù)組S的示意圖；圖6是本發(fā)明實(shí)施例提供的欲設(shè)置拉平種子點(diǎn)的波形對(duì)比窗口的中心點(diǎn)，及其對(duì)應(yīng)的拉平種子點(diǎn)移動(dòng)量的示意圖；圖7是本發(fā)明實(shí)施例提供的求取拉伸采樣坐標(biāo)的示意圖一；圖8是本發(fā)明實(shí)施例提供的求取拉伸采樣坐標(biāo)的示意圖二 ；圖9a是本發(fā)明實(shí)施例中拉平種子點(diǎn)在處理前的疊前道集上的位置示意圖；圖9b是本發(fā)明實(shí)施例中經(jīng)過一次同相軸精細(xì)拉平處理后的示意圖；圖IOa是本發(fā)明實(shí)施例提供的同相軸精細(xì)拉平處理前的效果圖；圖IOb是本發(fā)明實(shí)施例提供的經(jīng)過一次同相軸精細(xì)拉平處理后的效果圖；圖IOc是本發(fā)明實(shí)施例提供的經(jīng)過二次同相軸精細(xì)拉平處理后的效果圖；圖IOd是本發(fā)明實(shí)施例提供的經(jīng)過三次同相軸精細(xì)拉平處理后的效果圖。
具體實(shí)施例方式為使本發(fā)明實(shí)施例的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚明白，下面結(jié)合實(shí)施例和附圖，對(duì)本發(fā)明實(shí)施例做進(jìn)一步詳細(xì)說明。在此，本發(fā)明的示意性實(shí)施例及其說明用于解釋本發(fā)明，但并不作為對(duì)本發(fā)明的限定。如圖I為本發(fā)明實(shí)施例提供的一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法流程圖。包括步驟101，對(duì)采集的地震數(shù)據(jù)進(jìn)行去噪、靜校正、動(dòng)校正的前期處理操作。通過去噪來消除面波、多次波等干擾波；通過靜校正消除由于地表不平坦、激發(fā)點(diǎn)深度不同等問題對(duì)于時(shí)距曲線形態(tài)的影響；通過動(dòng)校正使時(shí)距曲線與地下界面的形態(tài)一致。步驟102，抽取想要進(jìn)行同相軸精細(xì)拉平處理的疊前道集。這種疊前道集一般是共反射點(diǎn)道集(道與道的入射角不同)，也可以是方位角道集(道與道的方位角不同)，它們都可以用I行、J列的二維數(shù)組D來表示。如圖2所示，其給出了一個(gè)方位角道集。在本步驟中，所述I為每個(gè)地震道的采樣點(diǎn)數(shù)，例如I = 401 ;所述J為疊前道集的地震道數(shù)，例如J = 6。所述的二維數(shù)組D，其第i行、第j列的值D(i，j)表示疊前道集的第j個(gè)地震道的第i個(gè)采樣點(diǎn)的振幅值。步驟103，設(shè)定同相軸精細(xì)拉平處理參數(shù)Nw、Nm、Ns、Tolg、Tolit5其中，所述的Nw為波形對(duì)比窗口的大小，Nw為整數(shù)，例如Nw = 50個(gè)采樣點(diǎn)數(shù)，圖2 縱向單位是采樣點(diǎn)數(shù)。所述的Nm為波形對(duì)比窗口的移動(dòng)量，Nm為整數(shù),例如設(shè)定Nm = 17個(gè)采樣點(diǎn)數(shù)；所述的Ns為波形對(duì)比的搜索半徑，Ns為整數(shù)，例如設(shè)定Ns = 6個(gè)采樣點(diǎn)數(shù)；所述的TolgS群容差，它是判定是否在某個(gè)波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)的依據(jù)，當(dāng)波形對(duì)比窗口內(nèi)地震道相關(guān)系數(shù)的平均值大于Tolg時(shí)，則需要在該波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)，反之亦反，例如設(shè)定Tolg = O. 75 ；所述的Toli為個(gè)體容差，它是優(yōu)選某個(gè)地震道的相似地震道的依據(jù)，當(dāng)波形對(duì)比窗口內(nèi)某個(gè)地震道與基準(zhǔn)地震道的相似系數(shù)大于Toli時(shí)，則該地震道為基準(zhǔn)地震道的相似地震道，例如設(shè)定Toli = O. 85。步驟104，以Nw為窗口大小，以Nm為窗口移動(dòng)量，逐窗口判定是否在某個(gè)波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)，并計(jì)算拉平種子點(diǎn)移動(dòng)量。圖3給出了移動(dòng)波形對(duì)比窗口的示意圖，圖中用梯形框標(biāo)出了最初的5個(gè)波形對(duì)比窗口的位置。具體地，對(duì)于第k個(gè)波形對(duì)比窗口(從道集的第I個(gè)采樣點(diǎn)算起)，可以通過以下步驟4. I到步驟4. 3操作步驟4. 1，對(duì)于第k個(gè)波形對(duì)比窗口，計(jì)算地震道與地震道之間的最優(yōu)移動(dòng)量及最大相關(guān)系數(shù)。定義J行、J列的二維數(shù)組S用來存儲(chǔ)地震道與地震道之間的最優(yōu)移動(dòng)量。定義J行、J列的二維數(shù)組C用來存儲(chǔ)地震道與地震道之間的最大相關(guān)系數(shù)。圖4給出了當(dāng)k = 3時(shí)用來存儲(chǔ)地震道與地震道之間最大相關(guān)系數(shù)的二維數(shù)組C的示意圖，所述的二維數(shù)組C，其第J1行、第j2列的值j2)表示疊前道集的第J1個(gè)地震道和第j2個(gè)地震道的最大相關(guān)系數(shù)。圖5給出了當(dāng)k = 3時(shí)用來存儲(chǔ)地震道與地震道之間最優(yōu)移動(dòng)量的二維數(shù)組S的示意圖，所述的二維數(shù)組S，其第J1行、第j2列的值S U1, J2)表示疊前道集的第J1個(gè)地震道和第j2個(gè)地震道的最優(yōu)移動(dòng)量。所述的地震道與地震道之間最優(yōu)移動(dòng)量和最大相關(guān)系數(shù)，其計(jì)算方法為以計(jì)算 SfJ1, j2)和 CU1, j2)為例，設(shè) is = (k-l)Nm, ie = is+Nw,定義集合 =丨I 2 -,Arw+1}、
i =，定義長(zhǎng)度為(Nw+1)的向量f，f的第i個(gè)元素f(i)=
D(is+i-l, J1) /f / ,定義2NS+1個(gè)長(zhǎng)度為(Nw+1)的向量的第i個(gè)元素玢⑴=D(is+l+i-l，j2)，/ef，則最大相關(guān)系數(shù)Cd j2)由以下公式給出
權(quán)利要求
1.一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法，其特征在于包括， 步驟101，對(duì)采集的地震數(shù)據(jù)進(jìn)行前期處理； 步驟102，抽取想要進(jìn)行同相軸精細(xì)拉平處理的疊前道集，用I行、J列的二維數(shù)組D來表示，其中，I為每個(gè)地震道的采樣點(diǎn)數(shù)，J為疊前道集的地震道數(shù)，二維數(shù)組D，其第i行、第j列的值D(i，j)表示疊前道集的第j個(gè)地震道的第i個(gè)采樣點(diǎn)的振幅值； 步驟103，設(shè)定同相軸精細(xì)拉平處理參數(shù)-^、^、^、！^^、！^^，其中’所述的心為波形對(duì)比窗口的大小，Nw為整數(shù)；所述的Nm為波形對(duì)比窗口的移動(dòng)量，Nm為整數(shù)；所述的Ns為波形對(duì)比的搜索半徑，Ns為整數(shù)；所述的Tolg為群容差，是判定是否在某個(gè)波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)的依據(jù)；所述的Toli為個(gè)體容差，是優(yōu)選某個(gè)地震道的相似地震道的依據(jù)； 步驟104，以Nw為窗口大小，以Nm為窗口移動(dòng)量，逐窗口判定是否在某個(gè)波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)，并計(jì)算拉平種子點(diǎn)移動(dòng)量； 步驟105，計(jì)算拉伸采樣坐標(biāo)； 步驟106，根據(jù)拉伸采樣坐標(biāo)實(shí)施同相軸精細(xì)拉平。
2.根據(jù)權(quán)利要求I所述的方法，其特征在于，對(duì)采集的地震數(shù)據(jù)進(jìn)行前期處理中包括，對(duì)采集的地震數(shù)據(jù)進(jìn)行去噪、靜校正、動(dòng)校正的前期處理。
3.根據(jù)權(quán)利要求I所述的方法，其特征在于，所述疊前道集包括共反射點(diǎn)道集和方位角道集。
4.根據(jù)權(quán)利要求I所述的方法，其特征在于，所述逐窗口判定是否在某個(gè)波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)，并計(jì)算拉平種子點(diǎn)移動(dòng)量具體包括， 步驟401，對(duì)于第k個(gè)波形對(duì)比窗口，計(jì)算地震道與地震道之間最優(yōu)移動(dòng)量的二維數(shù)組S，其中二維數(shù)組S中的第J1行、第j2列的值SU1, J2)表示疊前道集的第J1個(gè)地震道和第j2個(gè)地震道的最優(yōu)移動(dòng)量，并計(jì)算地震道與地震道之間最大相關(guān)系數(shù)的二維數(shù)組C，其中二維數(shù)組c中的第J1行、第j2列的值CU1, J2)表示疊前道集的第J1個(gè)地震道和第j2個(gè)地震道的最大相關(guān)系數(shù)；所述的地震道與地震道之間最優(yōu)移動(dòng)量SU1, j2)和地震道與地震道之間最大相關(guān)系數(shù)CU1, j2)的計(jì)算方法為設(shè)^= (k-l)Nm，L = is+Nw，定義集合/ = {1,2,l = {-Ns,-Ns+l,-,0,-,Ns-l,Ns}，定義長(zhǎng)度為(Nw+1)的向量 f，f 的第i個(gè)元素f(i) = D(is+i-l, jjdef ’定義2NS+1個(gè)長(zhǎng)度為(Nw+1)的向量gp/ei’g!的第i個(gè)元素gji) = D(is+l+i-l, j2)ei，則最大相關(guān)系數(shù)CU1, j2)由以下公式給出 最優(yōu)移動(dòng)量sa，」_2)，即取得ca，j2)時(shí)的I值，由以下公式給出SU J2) = argιηψ , I fvf-fvg；-g； 步驟402，對(duì)于第k個(gè)波形對(duì)比窗口，判定是否在該波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)，其中，所述的判定是否設(shè)置拉平種子點(diǎn)的方法為求取二維數(shù)組C的所有元素的平均值Cm
5.根據(jù)權(quán)利要求4所述的方法，其特征在于，在計(jì)算拉伸采樣坐標(biāo)中，定義I行、J列的二維數(shù)組Xnrat用來存儲(chǔ)拉伸采樣坐標(biāo)，所述的計(jì)算拉伸采樣坐標(biāo)的方法為假設(shè)需要在Q個(gè)波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)，設(shè)其中第q個(gè)中心點(diǎn)以及它所對(duì)應(yīng)的拉平種子點(diǎn)移動(dòng)量分別為i,和'，q= 1，2，一，Q，定義集合及二扎/2,···，^}，及1 =iU{l,/}，則對(duì)于任意一個(gè)定值j，_/ e i
6.根據(jù)權(quán)利要求4所述的方法，其特征在于，根據(jù)拉伸采樣坐標(biāo)實(shí)施同相軸精細(xì)拉平進(jìn)一步包括，定義I行、J列的二維數(shù)組X用來存儲(chǔ)原始采樣坐標(biāo)，定義I行、J列的二維數(shù)組Dn 用來存儲(chǔ)同相軸精細(xì)拉平處理后的疊前道集； 設(shè)置原始采樣坐標(biāo)，具體地，X(i, j) = i, i = 1，2，…，I; j = I, 2，…，J ; 對(duì)于任意一個(gè)定值j，Jei，以X的第j列\(zhòng)為自變量、以D的第j列Dj為函數(shù)值，構(gòu)建三次樣條插值函數(shù)J = /x,Di W ’則D_'(U) = Λ,,Di (XhU)。
7.根據(jù)權(quán)利要求6所述的方法，其特征在于，在根據(jù)拉伸采樣坐標(biāo)實(shí)施同相軸精細(xì)拉平之后還包括步驟107，用Dnew代替步驟102中的D，并在步驟103中設(shè)定新的同相軸精細(xì)拉平處理參數(shù)，再次執(zhí)行步驟104-步驟106。
8.根據(jù)權(quán)利要求7所述的方法，其特征在于，重復(fù)執(zhí)行所述步驟107兩至三次。
9.根據(jù)權(quán)利要求7所述的方法，其特征在于，所述在步驟103中設(shè)定新的同相軸精細(xì)拉平處理參數(shù)進(jìn)一步包括，按照參數(shù)設(shè)定的先后次序Nw、Nffl, Ns逐漸減小，Tolg, Toli保持不變或增大。
全文摘要
本發(fā)明涉及地球物理勘探中的地震疊前道集處理方法，尤其涉及一種針對(duì)地震疊前道集的同相軸精細(xì)拉平處理方法，包括對(duì)采集的地震數(shù)據(jù)進(jìn)行前期處理；抽取想要進(jìn)行同相軸精細(xì)拉平處理的疊前道集，用I行、J列的二維數(shù)組D來表示；設(shè)定同相軸精細(xì)拉平處理參數(shù)；以Nw為窗口大小，以Nm為窗口移動(dòng)量，逐窗口判定是否在某個(gè)波形對(duì)比窗口的中心點(diǎn)設(shè)置拉平種子點(diǎn)，并計(jì)算拉平種子點(diǎn)移動(dòng)量；計(jì)算拉伸采樣坐標(biāo)；根據(jù)拉伸采樣坐標(biāo)實(shí)施同相軸精細(xì)拉平。通過本發(fā)明實(shí)施例的方法，該方法對(duì)于地震資料的處理和解釋都具有重要的意義，并具有穩(wěn)定、高效的特點(diǎn)。
文檔編號(hào)G01V1/28GK102879821SQ201210363749
公開日2013年1月16日 申請(qǐng)日期2012年9月26日 優(yōu)先權(quán)日2012年9月26日
發(fā)明者楊昊, 張研 申請(qǐng)人:中國(guó)石油天然氣股份有限公司