[發明專利]提取下行子波和衰減參數的方法、裝置

有效

• 申請號: 202010063248.8
• 申請日: 2020-01-20
• 公開（公告）號: CN113138419B
• 公開（公告）日: 2022-05-10
• 發明（設計）人: 李彥鵬;王靜
• 申請（專利權）人: 中國石油天然氣集團有限公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/40
• 分類號: G01V1/40;G01V1/50;G01V1/30
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 王天堯;湯在彥

1.一種提取下行子波和衰減參數的方法，其特征在于，包括：

對下行縱波數據中初至時窗長度為設定值的地震道數據進行傅里葉變換，獲得傅里葉變換后的下行縱波數據和多道下行縱波對數頻譜；

在多道下行縱波對數頻譜中，去除下行子波對數頻譜，獲得滿足設定要求的多道子波校正后的下行縱波對數頻譜，所述設定要求為多道下行縱波對數頻譜的參數數據與多道子波校正后的下行縱波對數頻譜的參數數據的差值在設定范圍內，所述參數數據包括截距和斜率；

根據多道子波校正后的下行縱波對數頻譜的參數數據，對傅里葉變換后的下行縱波數據進行傅里葉反變換，獲得下行子波；

根據縱波初至時間和多道子波校正后的下行縱波對數頻譜的參數數據，獲得衰減參數，所述衰減參數包括真振幅恢復因子和地層品質因子；

在多道下行縱波對數頻譜中，去除下行子波對數頻譜，獲得多道子波校正后的下行縱波對數頻譜，包括：對當前每道下行縱波對數頻譜進行拉平處理；對拉平處理后的當前多道下行縱波對數頻譜進行疊加，獲得下行子波對數頻譜；將當前每道下行縱波對數頻譜減去下行子波對數頻譜，獲得新的每道下行縱波對數頻譜；在新的每道下行縱波對數頻譜的參數數據與當前每道下行縱波對數頻譜的參數數據的差值不在設定范圍內時，將新的每道下行縱波對數頻譜替換當前每道下行縱波對數頻譜，迭代執行以上步驟，直至新的每道下行縱波對數頻譜的參數數據與當前每道下行縱波對數頻譜的參數數據的差值在設定范圍內，確定最新的多道下行縱波對數頻譜為多道子波校正后的下行縱波對數頻譜。

2.如權利要求1所述的提取下行子波和衰減參數的方法，其特征在于，還包括：

根據VSP地震波場數據，獲得下行縱波數據和縱波初至時間。

3.如權利要求2所述的提取下行子波和衰減參數的方法，其特征在于，還包括：

根據縱波初至時間，對下行縱波數據進行拉平處理；

對下行縱波數據中初至時窗長度為設定值的地震道數據進行傅里葉變換，包括：

對拉平處理后的下行縱波數據中初至時窗長度為設定值的地震道數據進行傅里葉變換。

4.如權利要求1所述的提取下行子波和衰減參數的方法，其特征在于，根據多道子波校正后的下行縱波對數頻譜的參數數據，對傅里葉變換后的下行縱波數據進行傅里葉反變換，獲得下行子波，包括：

根據多道子波校正后的下行縱波對數頻譜的參數數據，對傅里葉變換后的下行縱波數據的實部和虛部進行校正，獲得校正后的下行縱波數據；

確定校正后的下行縱波數據的有效頻段；

對所述有效頻段之外的預設范圍內的數據點進行斜坡處理；

對所述有效頻段和預設范圍之外的下行縱波數據進行清零處理，獲得目標下行縱波數據；

對目標下行縱波數據進行傅里葉反變換，獲得下行子波。

5.如權利要求1所述的提取下行子波和衰減參數的方法，其特征在于，根據縱波初至時間和多道子波校正后的下行縱波對數頻譜的參數數據，獲得衰減參數，包括：

根據縱波初至時間和多道子波校正后的下行縱波對數頻譜的截距，確定真振幅恢復因子；

根據縱波初至時間和多道子波校正后的下行縱波對數頻譜的斜率，獲得地層品質因子。

6.如權利要求5所述的提取下行子波和衰減參數的方法，其特征在于，根據縱波初至時間和多道子波校正后的下行縱波對數頻譜的截距，確定真振幅恢復因子，包括：

對縱波初至時間的對數和多道子波校正后的下行縱波對數頻譜的截距進行線性擬合，確定線性擬合的斜率為真振幅恢復因子。

7.如權利要求5所述的提取下行子波和衰減參數的方法，其特征在于，所述地層品質因子包括等效品質因子和層品質因子；

根據縱波初至時間和多道子波校正后的下行縱波對數頻譜的斜率，獲得地層品質因子，包括：

根據縱波初至時間和每道子波校正后的下行縱波對數頻譜的斜率，獲得每道對應的等效品質因子；

根據縱波初至時間和每道對應的等效品質因子，獲得每道對應的層品質因子。