[發明專利]一種提升海洋深反射地震數據中低頻信號能量的處理方法

有效

• 申請號: 202110791360.8
• 申請日: 2021-07-13
• 公開（公告）號: CN113514889B
• 公開（公告）日: 2022-06-21
• 發明（設計）人: 成谷
• 申請（專利權）人: 中山大學
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 北京集佳知識產權代理有限公司 11227
• 代理人: 沈闖

1.一種提升海洋深反射地震數據中低頻信號能量的處理方法，其特征在于，包括：

獲取觀測地震數據及沉放深度，所述觀測地震數據為包含鬼波和有效波的地震數據，所述沉放深度為地震激發源及接收電纜距離水面的深度；

計算所述觀測地震數據的頻域特性，得到地震數據頻譜；

根據所述沉放深度計算反鬼波濾波器的頻域特性，得到反鬼波濾波器頻譜；

將反鬼波濾波器頻譜的橫軸頻率和縱軸振幅值用兩個一維數組f和A來表達，在振幅數組A中尋找第一個幅值等于1的元素，記錄該位置為j，則數組f中的第j個元素對應的頻率值為高頻端與低頻端的分界點，對所述反鬼波濾波器頻譜的高頻端幅值均賦值為1，保留低頻端頻譜分布，得到只提升低頻信號能量的反鬼波濾波器頻譜，其中，所述高頻端為大于首個幅值為1的頻率點的頻率范圍，所述低頻端為小于首個幅值為1的頻率點的頻率范圍；

將所述地震數據頻譜與所述只提升低頻信號能量的反鬼波濾波器頻譜相乘，得到提升低頻信號能量的地震數據頻譜。

2.根據權利要求1所述的方法，其特征在于，對所述反鬼波濾波器頻譜的高頻端幅值均賦值為1，保留低頻端頻譜分布，得到只提升低頻信號能量的反鬼波濾波器頻譜，包括：

確定首個幅值為1的頻率點，作為高頻端與低頻端的頻率分界點；

將大于所述頻率分界點的頻率范圍作為所述高頻端，將小于所述頻率分界點的頻率范圍作為所述低頻端；

將所述高頻端各個頻率點的幅值修改為1；

將所述低頻端各個頻率點的幅值保持原來數值。

3.根據權利要求1所述的方法，其特征在于，計算所述觀測地震數據的頻域特性，得到地震數據頻譜，包括：

利用傅里葉變換計算觀測地震數據的頻域特性，得到地震數據頻譜。

4.根據權利要求1所述的方法，其特征在于，根據所述沉放深度計算反鬼波濾波器的頻域特性，得到反鬼波濾波器頻譜，包括：

根據所述沉放深度計算所述鬼波與所述有效波到達所述接收電纜的時間差；

根據所述時間差計算鬼波濾波器的頻譜，得到鬼波濾波器頻譜；

根據所述鬼波濾波器頻譜計算反鬼波濾波器的頻譜，得到反鬼波濾波器頻譜。

5.根據權利要求4所述的方法，其特征在于，根據所述沉放深度計算所述鬼波與所述有效波到達所述接收電纜的時間差，包括：

根據所述沉放深度及地震波在水中的傳播速度，計算所述鬼波與所述有效波到達所述接收電纜的時間差；其中，所述時間差由所述地震激發源的沉放深度和所述接收電纜的沉放深度求和后除以所述地震波在水中的傳播速度得到。

6.根據權利要求4所述的方法，其特征在于，根據所述時間差計算鬼波濾波器的頻譜，得到鬼波濾波器頻譜，包括：

計算模擬有效波的頻域特性，得到模擬有效波的頻譜，所述模擬有效波為地震子波和地震激發脈沖信號褶積后的地震波信號；

將所述時間差作為模擬有效波與模擬鬼波之間的時間間隔，所述模擬鬼波為地震子波和鬼波脈沖信號褶積后的地震波信號；

依據所述時間間隔得到模擬有效波受模擬鬼波干擾的模擬受到鬼波影響的地震波，并計算所述模擬受到鬼波影響的地震波的頻譜；

根據所述模擬有效波的頻譜與所述模擬受到鬼波影響的地震波的頻譜，得到鬼波濾波器的頻譜。

7.根據權利要求4所述的方法，其特征在于，根據所述鬼波濾波器頻譜計算反鬼波濾波器的頻譜，得到反鬼波濾波器頻譜，包括：

對所述鬼波濾波器頻譜求倒數，得到反鬼波濾波器頻譜。

8.根據權利要求1所述的方法，其特征在于，還包括：

對所述提升低頻信號能量的地震數據頻譜進行傅里葉反變換至時間域，得到提升低頻信號能量的時間域地震數據。