[發明專利]低頻重構并行Marchenko成像方法

有效

• 申請號: 201710284885.6
• 申請日: 2017-04-27
• 公開（公告）號: CN107102355B
• 公開（公告）日: 2017-12-05
• 發明（設計）人: 靳中原;韓立國;単剛義;胡勇
• 申請（專利權）人: 吉林大學
• 主分類號: G01V1/30
• 分類號: G01V1/30;G01V1/28
• 代理公司: 長春吉大專利代理有限責任公司22201
• 代理人: 王立文

1.一種低頻重構并行Marchenko成像方法，其特征在于，地震干涉法被動源低頻重構地震超越干涉法主動源數據、分炮并行直達波估計和數據驅動面向復雜高速屏蔽地下構造的高精度Marchenko成像，包括以下步驟：

a、在Ubuntu14.04系統平臺上安裝Torque2.4.16，構建PBS批量作業管理系統；

b、對實際采集的被動源地震信號進行地震干涉法處理，基于互相關法，重構與主動源相同觀測系統的被動源地震數據；

c、對被動源地震數據進行多域迭代去噪，得到較高信噪比的被動源數據；補償去噪之后的被動源數據能量，找到被動源低頻端與主動源低頻端在頻譜上的交點，并確定焦點位置的頻率為優勢頻率；

d、結合低通濾波器分別對主動源、被動源在優勢頻率處做低通濾波；

e、設計匹配濾波器，計算匹配濾波算子，把低通濾波后的主動源低頻信息與低通濾波后的被動源低頻信息匹配濾波，得到匹配濾波后的主動源低頻信息；

f、對匹配濾波后的主動源低頻信息和主動源分別做歸一化，然后對二者做加和，然后基于主動源歸一化之前的主動源數據最大振幅進行能量恢復，得到合成后的寬頻帶地震數據；

g、對地震數據預處理，滿足Marchenko迭代求解格林函數的要求；

h、在不精確的模型上進行正演得到地下介質某點到地表布設的檢波器之間估計的直達波響應，經過低頻補償后的寬頻帶地震數據對地下復雜地質體進行成像，此時，PBS主處理器為這個從處理器按順序分配下一個正演計算作業，直到m次正演全部結果；

i、Marchenko成像需要求解成像區域的所有成像點的完全格林函數、上傳格林函數和下傳格林函數，成像區域m個虛擬震源位置對應著m次迭代求解格林函數，當一次Marchenko迭代求解完全格林函數、上傳格林函數和下傳格林函數結束，主處理器將繼續分配后續作業給從處理器處理，直到全部的作業結束；

j、將分解之后得到的上傳格林函數和下傳格林函數用來進行Marchenko成像，PBS主處理器收集記錄該虛擬震源位置的反射系數并記錄該虛擬震源在觀測系統中的坐標，然后為空閑的從處理器繼續分配作業，直到全部的作業計算完成；

k、當系統求得成像區域內全部虛擬震源位置的反射系數后，根據觀測系統，依據成像區域中的網格點坐標的先后順序，串行依次把反射系數排列到網格點的位置，最終得到基于低頻重構后的面向目標區域的數據驅動Marchenko成像結果。

2.按照權利要求1所述的低頻重構并行Marchenko成像方法，其特征在于，g步驟所述的對地震數據預處理，滿足Marchenko迭代求解格林函數的要求，有兩種情況：

A、數值模擬算例中，通過數值計算的方法來計算出寬頻帶地震數據中的直達波、自由表面多次波；

B、在實際地震勘探中，自由表面多次波的消去是通過SRME方法達到，直達波的計算是通過淺層調查進行計算得到。