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

有效

• 申請號: 201711259621.1
• 申請日: 2017-12-04
• 公開（公告）號: CN108196300B
• 公開（公告）日: 2019-06-11
• 發明（設計）人: 詹仕凡;李磊;陶春峰;萬忠宏;崔京彬;徐晨
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/36
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 李輝;陳剛

1.一種地震數據處理方法，其特征在于，提供有原始地震道數據；所述方法包括：

對所述原始地震道數據中的振幅值進行統計處理，得到概率密度函數；

基于所述原始地震道數據中的振幅值范圍，生成初始隨機數集合，從所述初始隨機數集合中確定符合所述概率密度函數的目標隨機數集合；

基于所述目標隨機數集合，確定所述原始地震道數據對應的多個子波數據，并將所述多個子波數據進行第一疊加，將第一疊加的結果作為擬合地震道數據；

將所述原始地震道數據減去所述擬合地震道數據，得到殘差地震道數據；

當所述殘差地震道數據中的振幅的絕對值的最大值大于或等于預設振幅閾值時，將所述殘差地震道數據乘以指定特征參數得到的結果作為新的原始地震道數據，并生成所述新的原始地震道數據對應的新的擬合地震道數據和新的殘差地震道數據，直至所述新的殘差地震道數據中振幅的絕對值的最大值小于所述預設振幅閾值，并將所述新的擬合地震道數據以及在所述新的擬合地震道數據之前生成的擬合地震道數據進行第二疊加，并將第二疊加的結果作為目標地震道數據；其中，所述指定特征參數與迭代次數相關聯。

2.根據權利要求1所述的一種地震數據處理方法，其特征在于，所述對所述原始地震道數據中的振幅值進行統計處理，得到概率密度函數，包括：

將所述原始地震道數據中數值為負數的振幅值替換為零，得到第一地震道數據；

當所述第一地震道數據中的指定采樣點的振幅值大于預設振幅閾值時，將所述預設振幅閾值作為所述指定采樣點的振幅值，得到第二地震道數據；

對所述第二地震道數據中的振幅值進行統計處理，得到所述概率密度函數。

3.根據權利要求1所述的一種地震數據處理方法，其特征在于，所述基于所述原始地震道數據中的振幅值范圍，生成初始隨機數集合，包括：

采用線性同余法生成多個均勻分布的隨機數；其中，所述多個均勻分布的隨機數的值在所述原始地震道數據中的振幅值范圍內；

將所述多個均勻分布的隨機數構成的集合作為所述初始隨機數集合。

4.根據權利要求3所述的一種地震數據處理方法，其特征在于，所述初始隨機數集合中隨機數的數量的取值范圍為：5000～10000。

5.根據權利要求1所述的一種地震數據處理方法，其特征在于，所述從所述初始隨機數集合中確定符合所述概率密度函數的目標隨機數集合，包括：

采用舍選法從所述初始隨機數集合中確定符合所述概率密度函數的目標隨機數集合。

6.根據權利要求1所述的一種地震數據處理方法，其特征在于，所述基于所述目標隨機數集合，確定所述原始地震道數據對應的多個子波數據，包括：

基于預設子波模型和所述目標隨機數集合，確定目標似然函數模型；

采用期望最大化方法確定所述目標似然函數模型中的子波參數集合中子波參數的值；其中，所述子波參數集合包括多組子波參數，所述子波參數包括子波的振幅系數、中心時間、主頻和相位；

根據所述預設子波模型和所述子波參數集合中子波參數的值，確定所述原始地震道數據對應的多個子波數據；其中，所述子波參數與所述子波數據一一對應，以使得一組子波參數對應一個子波數據。

7.根據權利要求6所述的一種地震數據處理方法，其特征在于，采用下述公式表征所述預設子波模型：

其中，g_γ(x)表示所述預設子波模型，α表示子波的振幅系數，u表示子波的中心時間，ζ表示子波的主頻，表示子波的相位，x表示采樣時間；

以及采用下述公式確定目標似然函數模型：

其中，L(Θ)表示目標似然函數，表示所述多個子波數據中第j個子波數據，x_z表示所述目標隨機數集合中第z個隨機數，P表示所述目標隨機數集合中隨機數的數量；表示所述子波參數集合，α_j表示第j個子波數據中的子波的振幅系數，u_j表示第j個子波數據中的子波的中心時間，ζ_j表示第j個子波數據中的子波的主頻，表示第j個子波數據中的子波的相位，M表示所述多個子波數據中子波數據的數量。