[發明專利]地震數據處理方法和裝置

有效

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

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

獲取目標區域的地震數據；

從所述地震數據中提取第一地震子波；

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

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

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

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

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

w_h(t)＝w_t(at)

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

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

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

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

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

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

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

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

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

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

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

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

d＝W_h·x

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