一種處理疊前寬線地震數(shù)據(jù)的方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及地質(zhì)勘探領(lǐng)域，具體說涉及一種處理疊前寬線地震數(shù)據(jù)的方法。
【背景技術(shù)】
[0002] 在地震勘探過程中，針對復(fù)雜地況的地區(qū)通常采用疊前偏移成像技術(shù)。由于疊前 偏移成像的成像結(jié)果的好壞直接影響著后續(xù)的解釋和油井的定位，因此疊前偏移成像作為 地震數(shù)據(jù)處理和油氣勘探過程中重要的一環(huán)，在復(fù)雜地區(qū)地震勘探中發(fā)揮了重要作用。
[0003] 在地表復(fù)雜多變的山前帶，由于地質(zhì)數(shù)據(jù)采集激發(fā)和接收條件差，因此導(dǎo)致獲取 到的疊前地震數(shù)據(jù)的信噪比極低。在現(xiàn)有技術(shù)中，通常采用寬線觀測系統(tǒng)采集地震數(shù)據(jù)，通 過增加覆蓋次數(shù)的方式從而獲取較高信噪比的疊前寬線地震數(shù)據(jù)。
[0004] 在現(xiàn)有技術(shù)中，對寬線觀測系統(tǒng)采集到的疊前寬線地震數(shù)據(jù)的常規(guī)處理是在原始 數(shù)據(jù)上進(jìn)行簡單去噪和相應(yīng)預(yù)處理后進(jìn)行橫向疊加，將其轉(zhuǎn)為二維單線數(shù)據(jù)。但是基于上 述方法得到的疊前寬線地震數(shù)據(jù)的處理結(jié)果并不理想，進(jìn)而導(dǎo)致利用上述數(shù)據(jù)處理結(jié)果進(jìn) 行疊前偏移成像的成像結(jié)果也并不理想。
[0005] 因此，為了獲得更為理想的疊前寬線地震數(shù)據(jù)處理結(jié)果，需要一種新的處理疊前 寬線地震數(shù)據(jù)方法。

【發(fā)明內(nèi)容】

[0006] 針對現(xiàn)有疊前寬線地震數(shù)據(jù)處理方法的處理結(jié)果并不理想的問題，本發(fā)明提供了 一種處理疊前寬線地震數(shù)據(jù)的方法，所述方法包含W下步驟：
[0007] 步驟一，獲取疊前寬線地震數(shù)據(jù)；
[0008] 步驟二，基于球面擬合技術(shù)將所述疊前寬線地震數(shù)據(jù)合并為疊前單線地震數(shù)據(jù)。
[0009] 在一實施例中，所述步驟二包含W下步驟：
[0010] 獲取旅行時面步驟，獲取所述疊前寬線地震數(shù)據(jù)中每個樣點對應(yīng)的局部旅行時 面；
[0011] 疊加步驟，根據(jù)所述局部旅行時面對所述疊前寬線地震數(shù)據(jù)進(jìn)行高斯加權(quán)疊加， 從而獲得所述疊前單線地震數(shù)據(jù)。
[0012] 在一實施例中，所述獲取旅行時面步驟包含W下步驟：
[0013] 獲取旅行時面描述步驟，利用球面擬合技術(shù)獲取所述局部旅行時面的描述；
[0014] 獲取運動學(xué)屬性參數(shù)步驟，基于所述局部旅行時面的描述獲取所述疊前寬線地震 數(shù)據(jù)中每個樣點對應(yīng)的運動學(xué)屬性參數(shù)；
[0015] 旅行時面確定步驟，基于所述局部旅行時面的描述W及所述運動學(xué)屬性參數(shù)獲取 所述疊前寬線地震數(shù)據(jù)中每個樣點對應(yīng)的局部旅行時面。
[0016] 在一實施例中，在所述獲取運動學(xué)屬性參數(shù)步驟中，所述運動學(xué)屬性參數(shù)為疊加 能量最大的參數(shù)組合，所述參數(shù)組合包括地表出射地震波的方位角、出射角和局部波前曲 率。
[0017] 在一實施例中，在所述疊加步驟中，針對所述疊前寬線地震數(shù)據(jù)，W目標(biāo)觀測單線 為中必，對于所述目標(biāo)觀測單線上觀測記錄的每個樣點實施所述高斯加權(quán)疊加。
[0018] 在一實施例中，所述方法還包含步驟H，根據(jù)所述疊前單線地震數(shù)據(jù)獲取平面波 數(shù)據(jù)，所述步驟Η包含W下步驟：
[0019] 分解步驟，將所述疊前單線地震數(shù)據(jù)分解為局部平面波數(shù)據(jù)；
[0020] 局部傾斜疊加步驟，對所述局部平面波數(shù)據(jù)進(jìn)行局部傾斜疊加從而獲取所述平面 波數(shù)據(jù)。
[0021] 在一實施例中，所述步驟Η還包含高斯束篩選偏移步驟，利用相似系數(shù)對所述平 面波數(shù)據(jù)進(jìn)行篩選。
[0022] 在一實施例中，所述高斯束篩選偏移步驟包含W下步驟：
[0023] 獲取相似系數(shù)步驟，獲取所述平面波數(shù)據(jù)中平面波分量對應(yīng)的相似系數(shù)；
[0024] 獲取相似系數(shù)闊值步驟，設(shè)定相似系數(shù)的闊值；
[0025] 確定信號平面波分量步驟，基于所述闊值根據(jù)所述平面波分量對應(yīng)的相似系數(shù)確 定信號平面波分量。
[0026] 在一實施例中，在所述確定信號平面波分量步驟中，大于所述闊值的所述相似系 數(shù)所對應(yīng)的平面波分量為所述信號平面波分量。
[0027] 在一實施例中，所述方法還包含步驟四，利用所述信號平面波分量進(jìn)行高斯束延 拓及成像。
[0028] 與現(xiàn)有技術(shù)相比，利用本發(fā)明的處理方法處理疊前寬線地震數(shù)據(jù)在保證數(shù)據(jù)處理 結(jié)果正確性的前提下可W進(jìn)一步提高數(shù)據(jù)的信噪比，從而最終得到更為理想的疊前偏移成 像結(jié)果。
[0029] 本發(fā)明的其它特征或優(yōu)點將在隨后的說明書中闡述。并且，本發(fā)明的部分特征或 優(yōu)點將通過說明書而變得顯而易見，或者通過實施本發(fā)明而被了解。本發(fā)明的目的和部分 優(yōu)點可通過在說明書、權(quán)利要求書W及附圖中所特別指出的步驟來實現(xiàn)或獲得。
【附圖說明】
[0030] 附圖用來提供對本發(fā)明的進(jìn)一步理解，并且構(gòu)成說明書的一部分，與本發(fā)明的實 施例共同用于解釋本發(fā)明，并不構(gòu)成對本發(fā)明的限制。在附圖中：
[0031] 圖1是根據(jù)本發(fā)明一實施例執(zhí)行流程圖；
[0032] 圖2是本發(fā)明一實施例的理論模型圖；
[0033] 圖3是基于本發(fā)明一實施例的理論模型抽取的任意炮道集成像圖；
[0034] 圖4是基于本發(fā)明一實施例的理論模型炮道集加隨機(jī)噪音后的成像圖；
[0035] 圖5是基于本發(fā)明一實施例的理論模型橫向疊加合成的單線炮道集成像圖；
[0036] 圖6是基于本發(fā)明一實施例的理論模型球面擬合合成的單線炮道集成像圖；
[0037] 圖7是基于本發(fā)明一實施例的理論模型橫向疊加后高斯束偏移成像圖；
[0038] 圖8是基于本發(fā)明一實施例的理論模型球面擬合后高斯束偏移成像圖；
[0039] 圖9是本發(fā)明一實施例的任意炮道集實際數(shù)據(jù)成像圖；
[0040] 圖10是基于本發(fā)明一實施例的實際數(shù)據(jù)球面擬合后炮道集成像圖；
[0041] 圖11是基于本發(fā)明一實施例的實際數(shù)據(jù)球面擬合后高斯束偏移剖面結(jié)果成像 圖；
[0042] 圖12是基于本發(fā)明一實施例的實際數(shù)據(jù)球面擬合后克?；舴蚱破拭娼Y(jié)果成像 圖；
[0043] 圖13是基于本發(fā)明一實施例的實際數(shù)據(jù)球面擬合后高斯束偏移成像道集成像 圖；
[0044] 圖14是基于本發(fā)明一實施例的實際數(shù)據(jù)球面擬合后克希霍夫偏移成像道集成像 圖。
【具體實施方式】
[0045] W下將結(jié)合附圖及實施例來詳細(xì)說明本發(fā)明的實施方式，借此本發(fā)明的實施人員 可W充分理解本發(fā)明如何應(yīng)用技術(shù)手段來解決技術(shù)問題，并達(dá)成技術(shù)效果的實現(xiàn)過程并依 據(jù)上述實現(xiàn)過程具體實施本發(fā)明。需要說明的是，只要不構(gòu)成沖突，本發(fā)明中的各個實施 例W及各實施例中的各個特征可W相互結(jié)合，所形成的技術(shù)方案均在本發(fā)明的保護(hù)范圍之 內(nèi)。
[0046] 在現(xiàn)有技術(shù)中，通常采用寬線觀測系統(tǒng)來采集疊前地震數(shù)據(jù)，從而提高數(shù)據(jù)的信 噪比。但是針對寬線觀測系統(tǒng)采集到的疊前寬線地震數(shù)據(jù)，利用現(xiàn)有的數(shù)據(jù)處理方法得到 的數(shù)據(jù)處理結(jié)果并不理想。因此本發(fā)明提出了一種處理疊前寬線地震數(shù)據(jù)的方法W獲取更 為理想的疊前寬線地震數(shù)據(jù)處理結(jié)果。下面基于流程圖來詳細(xì)說明本發(fā)明的一實施例的具 體執(zhí)行過程。
[0047] 附圖的流程圖中示出的步驟可W在包含諸如一組計算機(jī)可執(zhí)行指令的計算機(jī)系 統(tǒng)中執(zhí)行。雖然在流程圖中示出了各步驟的邏輯順序，但是在某些情況下，可不同于此 處的順序執(zhí)行所示出或描述的步驟。
[0048] 如圖1所示，執(zhí)行本發(fā)明的處理方法首先要執(zhí)行步驟S100,獲取疊前寬線地震數(shù) 據(jù)。在送之后就可W對疊前寬線地震數(shù)據(jù)進(jìn)行處理。
[0049] 本發(fā)明的處理方法的中必是基于球面擬合技術(shù)將疊前寬線地震數(shù)據(jù)合并為信噪 比加強(qiáng)的疊前單線地震數(shù)據(jù)。由于在合并過程中首先要確定疊前寬線地震數(shù)據(jù)中每個樣點 對應(yīng)的局部旅行時面。因此首先執(zhí)行步驟S110,獲取旅行時面描述，利用球面擬合技術(shù)獲取 疊前寬線地震數(shù)據(jù)的局部旅行時面的描述。
[0050] 在本實施例中，將疊前寬線地震數(shù)據(jù)視為窄方位的Η維地震資料，利用地震數(shù)據(jù) 的局部同相性，對于疊前寬線地震數(shù)據(jù)中的每個樣點，將其鄰域的數(shù)據(jù)面近似為一個雙曲 面或拋物面。在本實施例中選取式（1)來近似描述疊前寬線地震數(shù)據(jù)的局部旅行時面。
[0051]
[0052] 式（1)中：
[0053] Xm、ym分別為局部鄰域范圍區(qū)域Ω內(nèi)任意一點的橫坐標(biāo)和縱坐標(biāo)；
[0054] X。、y。分別為區(qū)域Ω內(nèi)的中也點橫坐標(biāo)和縱坐標(biāo)；
[00巧]t。為中必點處的地震波旅行時；
[0056] t為區(qū)域Ω內(nèi)任意一點的地震波旅行時；
[0057] W:，稱0是中點矢量；
[005引 為是偏移距矢量；
[0059] V。是近地表速度；
[0060] λ、α、整W分別是地表出射地震波的方位角、出射角和局部波前曲率。
[0061] 由于對于疊前寬線地震數(shù)據(jù)中每個樣點，都對應(yīng)一個由式（1)定義的λ、α、整《 Η參數(shù)構(gòu)成的局部旅行時面，當(dāng)疊加能量為最大時對應(yīng)的局部旅行時面為本發(fā)明處理方法 合并過程中所需的局部旅行時面。接下來執(zhí)行步驟S120,獲取疊加能量為最大時對應(yīng)的 入、α、運W.H參數(shù)。
[0062] 定義疊加能量為最大時λ、α、足、.Ξ個參數(shù)的參數(shù)組合（λ，α，塞W )m。、為運動 學(xué)屬性參數(shù)。步驟S120即為獲取運動學(xué)屬性參數(shù)步驟，基于局部旅行時面的描述獲取疊前 寬線地震數(shù)據(jù)中每個樣點對應(yīng)的運動學(xué)屬性參數(shù)。在本實施例中，基于球面擬合技術(shù)引入 局部球面近似，波前曲率里W則簡化為：
[0065] 首先利用式（3)進(jìn)行關(guān)于參數(shù)λ、α的二維搜索，即掃描不同的（λ, α)。對每 一組（λ, α)沿式（3)描述的旅行時面按用戶自定義的孔徑進(jìn)行疊加，尋找到疊加能量最 大時對應(yīng)的參數(shù)組合（λ，α )m。、，從而得到每個樣點處對應(yīng)的兩個運動學(xué)屬性參數(shù)（λ，α ) m。、；然后再將搜索到的（λ, a)m。、代入式（1)，實施參數(shù)的一維搜索，即掃描不同的基w， 沿式（1)描述的旅行時面按用戶自定義的孔徑進(jìn)行疊加，同樣尋找疊加能量為最大時對應(yīng) 的參數(shù)組合（λ，α，^、）m。、，便可得到每個樣點處最終的運動學(xué)屬性參數(shù)（λ，α，)m。、。
[0066] 接下來執(zhí)行S130,旅行時面確定步驟。基于步驟SllO中得到的局部旅行時面的描 述（式1) W及步驟S120中獲取的運動學(xué)屬性參