[發(fā)明專利]地震數(shù)據(jù)處理方法和裝置

有效

• 申請?zhí)枺?/td>: 201810478911.3
• 申請日: 2018-05-18
• 公開（公告）號: CN108693558B
• 公開（公告）日: 2020-09-08
• 發(fā)明（設計）人: 杜昕;李皓;李國發(fā);周振曉;張文波;唐博文;王偉
• 申請（專利權）人: 中國石油天然氣集團有限公司;中國石油集團東方地球物理勘探有限責任公司;中國石油大學(北京)
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京三友知識產(chǎn)權代理有限公司 11127
• 代理人: 李輝;徐煥

1.一種地震數(shù)據(jù)處理方法，其特征在于，包括：

獲取目標區(qū)域的地震數(shù)據(jù)；

從所述地震數(shù)據(jù)中提取第一地震子波；

在時間域中對所述第一地震子波進行壓縮處理，以獲取第二地震子波；

根據(jù)所述地震數(shù)據(jù)、所述第一地震子波、所述第二地震子波建立反演目標函數(shù)，其中，所述反演目標函數(shù)為基于反演機理建立的用于進行高分辨率地震數(shù)據(jù)處理的函數(shù)；

根據(jù)所述反演目標函數(shù)，確定處理后的地震數(shù)據(jù)；

其中，在時間域中對所述第一地震子波進行壓縮處理，以獲取第二地震子波，包括：

根據(jù)傅里葉變換伸縮性，按照以下算式對所述第一地震子波進行壓縮處理，以獲取所述第二地震子波：

w_h(t)＝w_t(at)

其中，w_h(t)為第二地震子波，w_t(t)為第一地震子波，t為時間點的編號，a為壓縮因子；

其中，所述第二地震子波為高頻地震子波；

其中，根據(jù)所述地震數(shù)據(jù)、所述第一地震子波、所述第二地震子波建立反演目標函數(shù)，包括：按照以下公式，建立所述反演目標函數(shù)：

其中，minimize表示求取最小值，表示求取L2范數(shù)，y表示反演目標函數(shù)的解析解，W為基于第一地震子波的第一地震子波矩陣，W_h為基于第二地震子波的第二地震子波矩陣，x為地震數(shù)據(jù)的列向量形式，μ為正則化參數(shù)。

2.根據(jù)權利要求1所述的方法，其特征在于，根據(jù)所述第二地震子波，按照以下公式確定所述基于第二地震子波的第二地震子波矩陣：

其中，W_h為基于第二地震子波的第二地震子波矩陣，w_h(i)為第二地震子波中編號為i的時間點的振幅值，n為第二地震子波時間點的總數(shù)。

3.根據(jù)權利要求1所述的方法，其特征在于，根據(jù)所述反演目標函數(shù)，確定處理后的地震數(shù)據(jù)，包括：

根據(jù)所述反演目標函數(shù)，求取所述反演目標函數(shù)的導數(shù)為0時所對應的解析解，并將所述解析解作為所述處理后的地震數(shù)據(jù)。

4.根據(jù)權利要求3所述的方法，其特征在于，按照以下公式求取所述反演目標函數(shù)的導數(shù)為0時所對應的解析解：

y′＝(W^TW+μI)^-1W^Td

其中，(·)^-1表示矩陣求逆，y′為反演目標函數(shù)的導數(shù)為0時所對應的解析解，W為基于第一地震子波的第一地震子波矩陣，W^T為基于第一地震子波的第一地震子波矩陣轉置，I為單位矩陣，d為模型向量，μ為正則化參數(shù)。

5.根據(jù)權利要求4所述的方法，其特征在于，按照以下方式確定所述模型向量：

d＝W_h·x

其中，d為模型向量，W_h為基于第二地震子波的第二地震子波矩陣，x為地震數(shù)據(jù)的列向量。