一種調(diào)整用于疊前偏移成像的地震數(shù)據(jù)能量一致性的方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及地震勘探資料處理技術(shù)，具體地，涉及一種調(diào)整用于疊前偏移成像的 地震數(shù)據(jù)能量一致性的方法。
【背景技術(shù)】
[0002] 目前，隨著油氣勘探的深入和數(shù)字計算技術(shù)的飛速發(fā)展，疊前偏移已取代疊后偏 移，成為解決復雜地質(zhì)條件精細成像的常規(guī)處理技術(shù)。特別是基于射線理論的Kichhoff?積 分法疊前偏移方法，在目前工業(yè)生產(chǎn)中應用最為廣泛。
[0003] Kichhoff積分法偏移的原理是沿橢球面對地層振幅求和，相鄰面元能量不均通常 會產(chǎn)生畫弧現(xiàn)象，影響地震資料的成像精度。疊前偏移對道集數(shù)據(jù)的能量一致性要求相對 較高，為了消除地震波在傳播過程中球面擴散、吸收衰減、震源能量差異、檢波器耦合差異 以及不同時間采集的區(qū)塊間的能量差異對反射波振幅的影響，現(xiàn)有的處理流程中通常采用 幾何擴散補償、地表一致性振幅補償和剩余振幅補償來消除能量差異，盡量提高振幅保真 度，保持地震波振幅能量的相對強弱關(guān)系。
[0004] 為更好的改善能量不均對疊前時間偏移成像的影響，目前在工業(yè)化生產(chǎn)處理中， 正逐步得到重視。吳清嶺等提出基于覆蓋次數(shù)的疊前振幅歸一化處理，改善振幅一致性處 理對連片處理的影響；張學濤等為了解決在地震資料連片處理中，能量不均嚴重影響疊前 時間偏移成像質(zhì)量的問題，采用了 CMP域逐點自動能量調(diào)整法；張偉等在假設(shè)共偏移距道 集覆蓋次數(shù)均勻情況下，分析了能量一致性對偏移成像影響，研究認為如果處理過程中振 幅補償不均就會造成分辨率橫向變化。綜合采用幾何擴散補償、地表一致性振幅補償和剩 余振幅補償來消除能量差異，利用監(jiān)控手段保證補償效果，可以有效提高偏移成像質(zhì)量。

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

[0005] 目前，地震資料的處理中通常仍是采用補償前后的道集、純波疊加剖面以及均方 根振幅屬性的監(jiān)控，判斷能量補償?shù)男Ч?br>[0006] 但是，在實現(xiàn)本發(fā)明的過程中，本發(fā)明人發(fā)現(xiàn)：為有效提高偏移成像質(zhì)量，目前的 這些監(jiān)控方法對于疊后偏移是能滿足的，但對于疊前偏移仍是不足的。原因在于疊前偏移 的工業(yè)化算法設(shè)計，其在實現(xiàn)過程中沒有考慮觀測系統(tǒng)設(shè)計和切除等原因造成的成像點間 的覆蓋次數(shù)差異以及炮檢距的分布的均勻性帶來的能量差異。為利用疊前偏移方法更好的 提高成像質(zhì)量，一方面要充分利用質(zhì)量監(jiān)控做好偏移前振幅補償；另一方面還要從疊前偏 移的工業(yè)化算法著手，彌補算法的不足，在偏移前做好道集能量調(diào)整，以解決考慮觀測系統(tǒng) 設(shè)計和切除等原因造成的成像點間的覆蓋次數(shù)差異以及炮檢距的分布的均勻性帶來的能 量差異。
[0007] 為此，非常需要一種調(diào)整用于疊前偏移成像的地震數(shù)據(jù)能量一致性的方法。
[0008] 為了實現(xiàn)上述目的，本發(fā)明實施例提供一種調(diào)整用于疊前偏移成像的地震數(shù)據(jù)能 量一致性的方法，包括：
[0009] 步驟I，采集地震勘探數(shù)據(jù)；
[0010] 步驟2,對所述地震勘探數(shù)據(jù)進行預處理，并分選成共中心點道集；
[0011] 步驟3,將所述共中心點道集分組；
[0012] 步驟4,針對每一組都進行如下處理：
[0013] 步驟41，按照如下公式計算考慮覆蓋次數(shù)的有效疊加地震數(shù)據(jù)：
[0014]
[0015] 步驟42,按照如下公式計算不考慮覆蓋次數(shù)的有效疊加地震數(shù)據(jù)：
[0016]
[0017] 步驟43,按照如下公式計算能量調(diào)整因子：
[0018] opz (t) = szl(t)/sz2(t)
[0019] 步驟44,按照如下公式將所述能量調(diào)整因子應用于該組的所有地震數(shù)據(jù)：
[0020] A-；(〇 = .Yy (/)〇/;(/)
[0021] 其中，z表示該組在所述共中心點道集中的編號；
[0022] szl (t)表示考慮覆蓋次數(shù)的有效疊加地震數(shù)據(jù)；
[0023] t表示時間；
[0024] N表示覆蓋次數(shù)；
[0025] j表示該組的地震數(shù)據(jù)的編號；
[0026] X]⑴表示該組的地震數(shù)據(jù)；
[0027] sz2 (t)表示不考慮覆蓋次數(shù)的有效疊加地震數(shù)據(jù)；
[0028] opz (t)表示能量調(diào)整因子；
[0029] 表示應用能量因子調(diào)整后的地震數(shù)據(jù)；
[0030] 步驟5,得到能量調(diào)整后的、用于疊前偏移成像的地震數(shù)據(jù)。
[0031] 借助于上述技術(shù)方案，本發(fā)明利用考慮覆蓋次數(shù)的有效疊加除以不考慮覆蓋次數(shù) 的有效疊加計算疊前偏移前能量調(diào)整因子，不僅考慮了覆蓋次數(shù)的影響，更考慮到了實際 疊加數(shù)據(jù)的有效值，是純數(shù)據(jù)驅(qū)動，能量調(diào)整更穩(wěn)定，容易理解和操作。應用本發(fā)明進行疊 前偏移前能量調(diào)整，更加符合疊前偏移的工業(yè)化算法，可獲得更好的地震勘探成像質(zhì)量。
【附圖說明】
[0032] 為了更清楚地說明本發(fā)明實施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對實施例描述 中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本發(fā)明的一些 實施例，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動的前提下，還可以根據(jù)這些附 圖獲得其他的附圖。
[0033] 圖1是本發(fā)明示例性方法的具體流程圖；
[0034] 圖2是本發(fā)明實施例一的具體流程圖；
[0035] 圖3是未應用本實施例進行的疊前時間偏移成像結(jié)果；
[0036] 圖4是應用本實施例進行的疊前時間偏移成像結(jié)果。
【具體實施方式】
[0037] 下面將結(jié)合本發(fā)明實施例中的附圖，對本發(fā)明實施例中的技術(shù)方案進行清楚、完 整地描述，顯然，所描述的實施例僅僅是本發(fā)明一部分實施例，而不是全部的實施例?；?本發(fā)明中的實施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其他 實施例，都屬于本發(fā)明保護的范圍。
[0038] 本領(lǐng)域技術(shù)技術(shù)人員知道，本發(fā)明的實施方式可以實現(xiàn)為一種系統(tǒng)、裝置、設(shè)備、 方法或計算機程序產(chǎn)品。因此，本公開可以具體實現(xiàn)為以下形式，即：完全的硬件、完全的軟 件（包括固件、駐留軟件、微代碼等），或者硬件和軟件結(jié)合的形式。
[0039] 根據(jù)本發(fā)明的實施方式，提出了一種調(diào)整用于疊前偏移成像的地震數(shù)據(jù)能量一致 性的方法。
[0040] 在本文中，需要理解的是，附圖中的任何元素數(shù)量均用于示例而非限制，以及任何 命名都僅用于區(qū)分，而不具有任何限制含義。
[0041] 示例件方法
[0042] 本示例性方法提供一種調(diào)整用于疊前偏移成像的地震數(shù)據(jù)能量一致性的方法，如 圖1所示，該方法包括：
[0043] 步驟S11，采集地震勘探數(shù)據(jù)。
[0044] 步驟S12,對地震勘探數(shù)據(jù)進行預處理，并分選成共中心點道集。
[0045] 可選地，該步驟可以對地震勘探數(shù)據(jù)進行包括靜校正，去噪，球面擴散補償、地層 吸收衰減、地表一致性等能量補償，以及子波一致性校正等預處理。
[0046] 步驟S13,將共中心點道集分組。
[0047] 可選地，該步驟可以選擇按炮檢距分組，也可以選擇炮檢距和方位角分組，也可以 按炮檢距向量片技術(shù)（Offset Vector Tile，簡稱0VT)進行分組。
[0048] 分組后，每組在該共中心點道集中都有唯一的編號。
[0049] 步驟S14,針對每一組都進行如下處理：
[0050] 步驟S141，按照如下公式計算考慮覆蓋次數(shù)的有效疊加地震數(shù)據(jù)：
[0051]
(公式 1)
[0052] 步驟S142,按照如下公式計算不考慮覆蓋次數(shù)的有效疊加地震數(shù)據(jù)：
[0053]
(公式 2)
[0054] 步驟S143,按照如下公式計算能量調(diào)整因子：
[0055] opz(t) = szl(t)/sz2(t) (公式 3)
[0056] 步驟S144,按照如下公式將能量調(diào)整因子應用于該組的所有地震數(shù)據(jù)：
[0057] X11 (t) = X1(I