[發明專利]虛擬源成像的方法及裝置

在審

• 申請號: 202010670466.8
• 申請日: 2020-07-13
• 公開（公告）號: CN113933894A
• 公開（公告）日: 2022-01-14
• 發明（設計）人: 劉文平;趙楊;周陽田麗;雍銳;左佳卉;吳建發;楊學鋒;石學文;文山師;王星皓
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京三高永信知識產權代理有限責任公司 11138
• 代理人: 董亞軍

1.一種虛擬源成像的方法，其特征在于，所述方法包括：

從勘探區域中的多個檢波器中確定第一檢波器和第二檢波器，將所述第一檢波器作為第一虛擬源；

對于勘探區域中的多個地表震源中的每個地表震源，確定所述第一虛擬源接收的所述地表震源激發的地震波的實際下行波場和所述第二檢波器接收的所述地表震源激發的地震波的第一實際上行波場；

將所述實際下行波場和所述第一實際上行波場進行時間域互相關，得到所述第一虛擬源和所述第二檢波器的互相關道；

將所述實際下行波場與所述互相關道進行時間域褶積，得到第一預測上行波場；

將所述第一預測上行波場與所述第一實際上行波場進行歸一化互相關處理，得到所述地表震源的權重；

根據所述互相關道和所述地表震源的權重，確定所述第一虛擬震源對應的第一虛擬道集；

將所述勘探區域中的多個檢波器中除所述第一檢波器外的至少一個第三檢波器作為第二虛擬源，確定所述第二虛擬源對應的至少一個第二虛擬道集；

根據所述第一虛擬道集和至少一個所述第二虛擬道集，確定所述勘探區域的疊加剖面圖像。

2.根據權利要求1所述的方法，其特征在于，所述根據所述互相關道和所述地表震源的權重，確定所述第一虛擬震源對應的第一虛擬道集，包括：

根據所述互相關道和所述地表震源的權重，通過以下公式一，確定所述第一虛擬震源對應的第一虛擬道集；

公式一：

其中，表示所述地表震源的權重，表示時間域互相關；r_A，r_B和r_S分別表示所述第一檢波器的空間位置、所述第二檢波器的空間位置和所述地表震源的空間位置；t表示時移量；V_w(r_B|r_A；t)表示所述第一虛擬道集；D(r_A|r_S；t)表示所述實際下行波場；U(r_B|r_S；t)表示所述第一實際上行波場。

3.根據權利要求1所述的方法，其特征在于，所述將所述實際下行波場與所述互相關道進行時間域褶積，得到第一預測上行波場，包括：

將所述實際下行波場與所述互相關道進行時間域褶積，通過以下公式二，得到第一預測上行波場；

公式二：

其中，表示所述第一預測上行波場；*表示時間域褶積；t表示時移量；D(r_A|r_S；t)表示所述實際下行波場；表示所述互相關道。

4.根據權利要求1所述的方法，其特征在于，所述將所述第一預測上行波場與所述第一實際上行波場進行歸一化互相關處理，得到所述地表震源的權重，包括：

將所述第一預測上行波場與所述第一實際上行波場分別進行傅里葉變換，得到第二預測上行波場和第二實際上行波場；

根據所述第二預測上行波場和所述第二實際上行波場，通過以下公式三，對所述第二預測上行波場和第二實際上行波場進行歸一化互相關處理，得到所述地表震源的權重；

公式三：

其中，表示所述地表震源的權重；表示第二預測上行波場；U(r_B|r_S；f)表示第二實際上行波場，f表示傅里葉變換的頻率。

5.根據權利要求1所述的方法，其特征在于，所述將所述第一預測上行波場與所述第一實際上行波場進行歸一化互相關處理，得到所述地表震源的權重，包括：

獲取所述地表震源在目標層疊前同相軸的時間、頻率和波數；

根據所述時間，對所述第一預測上行波場進行小波變換得到所述第一預測上行波場對應的第一小波系數，以及，對所述第一實際上行波場進行小波變換得到所述第一實際上行波場對應的第二小波系數；

根據所述頻率和所述波數，對所述第一小波系數和所述第二小波系數進行轉換，得到第一時頻空間數據體和第二時頻空間數據體；

將所述第一時頻空間數據體和所述第二時頻空間數據體進行歸一化互相關處理，得到所述地表震源的權重。