[發明專利]海洋地震數據震源子波提取方法、裝置和存儲介質

有效

• 申請號: 202111083115.8
• 申請日: 2021-09-15
• 公開（公告）號: CN113777652B
• 公開（公告）日: 2022-05-27
• 發明（設計）人: 成谷
• 申請（專利權）人: 中山大學
• 主分類號: G01V1/32
• 分類號: G01V1/32;G01V1/36
• 代理公司: 北京集佳知識產權代理有限公司 11227
• 代理人: 駱英靜

1.一種海洋地震數據震源子波提取方法，其特征在于，包括如下步驟：

從觀測海洋地震數據中截取包含鬼波的直達波窗口，以獲得各地震道的疊加信號W_dg(t)，所述疊加信號W_dg(t)包含有當前地震道的直達波信號W_d(t)和震源鬼波信號W_g(t)；

對任一地震道的疊加信號W_dg(t)進行頻域轉換，獲得當前地震道的頻域疊加信號

計算當前地震道上鬼波與直達波的到時差Δt；

依據當前地震道上鬼波與直達波的到時差Δt，計算不同圓頻率ω時當前地震道對應的1-e^-iωΔt的數值；

依據當前地震道的頻域疊加信號和不同圓頻率ω時當前地震道對應的1-e^-iωΔt的數值，計算當前地震道的頻域直達波信號

對當前地震道的頻域直達波信號進行時域轉換，獲得當前地震道的地震震源子波W_d(t)。

2.如權利要求1所述的海洋地震數據震源子波提取方法，其特征在于，所述從觀測海洋地震數據中截取包含鬼波的直達波窗口，以獲得各地震道的疊加信號W_dg(t)，進一步包括：

以任一地震道上的直達波起跳時間為標準時間值，將各地震道上的起跳時間拉平至所述標準時間值。

3.如權利要求2所述的海洋地震數據震源子波提取方法，其特征在于，以第一個地震道上的直達波起跳時間作為所述標準時間值。

4.如權利要求1所述的海洋地震數據震源子波提取方法，其特征在于，所述對任一地震道的疊加信號W_dg(t)進行頻域轉換，獲得當前地震道的頻域疊加信號具體包括：

利用傅里葉變換計算任一地震道的疊加信號W_dg(t)對應的頻域，獲得當前地震道的頻域疊加信號

5.如權利要求1所述的海洋地震數據震源子波提取方法，其特征在于，所述計算當前地震道上鬼波與直達波的到時差Δt，具體包括：

從觀測系統中獲取震源與水面的垂向距離z_s、拖纜與水面的垂向距離z_g、震源與水聽器的水平距離x_sg，設震源至水聽器的距離為d_sg，震源相對于水面的虛震源至水聽器的距離為d_s′g，地震波在水中的傳播速度為V_w，當前地震道上直達波的到達時間為t_d，當前地震道上鬼波的到達時間為t_g；

依據公式計算當前地震道上直達波的到達時間t_d，及依據公式計算當前地震道上鬼波的到達時間t_g；

依據公式計算當前地震道上鬼波與直達波的到時差Δt。

6.如權利要求1所述的海洋地震數據震源子波提取方法，其特征在于，所述依據當前地震道的頻域疊加信號和不同圓頻率ω時當前地震道對應的1-e^-iωΔt的數值，計算當前地震道的頻域直達波信號具體包括：

依據公式計算當前地震道的頻域直達波信號

7.如權利要求1所述的海洋地震數據震源子波提取方法，其特征在于，所述對當前地震道的頻域直達波信號進行時域轉換，獲得當前地震道的地震震源子波W_d(t)，具體包括：

利用傅里葉反變換計算當前地震道的頻域直達波信號的時域，獲得當前地震道的地震震源子波W_d(t)。