[發明專利]一種多源地震的混合波場分離方法和裝置

在審

• 申請號: 201911293960.0
• 申請日: 2019-12-16
• 公開（公告）號: CN110967750A
• 公開（公告）日: 2020-04-07
• 發明（設計）人: 邱鐵成;莊祖垠;陳繼宗;焦敘明;王煒;謝濤;徐強;俞小林
• 申請（專利權）人: 中國海洋石油集團有限公司;中海油田服務股份有限公司
• 主分類號: G01V1/30
• 分類號: G01V1/30;G01V1/36
• 代理公司: 北京安信方達知識產權代理有限公司 11262
• 代理人: 解婷婷;龍洪

1.一種多源地震的混合波場分離方法，其特征在于，所述方法包括：

獲取多源地震的混合地震數據，并根據所述混合地震數據創建地質模型；

針對地質模型中的每個層段內的混合地震數據，分別進行多方向掃描，得到每個層段對應的傾角角度；

每個層段采用對應的傾角角度對所述混合地震數據進行中值濾波處理；

當所有層段執行完中值濾波處理后，獲得波場分離的地震數據。

2.根據權利要求1所述多源地震的混合波場分離方法，其特征在于，所述針對地質模型中的每個層段內的混合地震數據，分別進行多方向掃描，得到每個層段對應的傾角角度，包括：

對于每個層段內的混合地震數據分別進行如下操作：

選擇不同的傾角角度分別計算，獲得各傾角角度的方向矢量距離；

確定所述方向矢量距離中最小的值所對應的傾角角度為該層段對應的傾角角度。

3.根據權利要求2所述多源地震的混合波場分離方法，其特征在于，所述每個層段采用對應的傾角角度對所述混合地震數據進行中值濾波處理，包括：

步驟1.將所述混合地震數據進行偽解混處理，得到原始偽分離共炮集地震數據；令i＝0，初始分離結果等于0；

步驟2.i＝i+1，將前一次炮分離結果按照地震多源的激發時間進行混合，再進行偽解混處理，得到第i次偽分離共炮集地震數據；

步驟3.將所述第i次偽分離共炮集地震數據與所述原始偽分離共炮集地震數據相減得到第i次的剩余混疊噪音，按照預先設置的中值濾波窗口長度，進行中值濾波處理得到第i次的有效信號；

步驟4.將中值濾波處理所得到的第i次的有效信號與第i次炮分離結果進行相加，得到第i+1次炮分離結果；

步驟5.將得到多個炮分離結果按照多源地震的激發時間進行組合，并對組合得到的混合炮集記錄進行偽分離，計算該偽分離后的地震數據的信噪比；

步驟6.若所述信噪比達到預置閾值時，則停止計算，輸出炮分離結果；

若所述信噪比未達到預置閾值時，則減小中值濾波窗口大小，跳轉到步驟2。

4.根據權利要求3所述多源地震的混合波場分離方法，其特征在于，所述將所述第i次偽分離共炮集地震數據與所述原始偽分離共炮集地震數據相減得到第i次的剩余混疊噪音，按照預先設置的中值濾波窗口長度，進行中值濾波處理得到第i次的有效信號，包括：

將所述第i次偽分離共炮集地震數據與所述原始偽分離共炮集地震數據相減得到剩余混疊噪音；

將得到的剩余混疊噪音從第i次偽分離共炮集分選到共偏移距道集或共檢波點道集上；

按照預先設置的中值濾波窗口長度，進行多級中值濾波得到有效信號。

5.根據權利要求3所述多源地震的混合波場分離方法，其特征在于，所述計算該偽分離后的地震數據的信噪比包括：采用下式計算信噪比：

其中，S/N表示信噪比，P′表示原始偽分離共炮集地震數據，Γ表示地震多源混合矩陣算子，Γ^H為Γ的共軛轉置矩陣，P_i+1代表第i+1次炮分離結果。

6.一種多源地震的混合波場分離裝置，所述裝置包括：存儲器和處理器；其特征在于：

所述存儲器，用于保存用于多源地震的混合波場分離的程序；

所述處理器，用于讀取執行所述用于多源地震的混合波場分離的程序，執行如下操作：

獲取多源地震的混合地震數據，并根據所述混合地震數據創建地質模型；

針對地質模型中的每個層段內的混合地震數據，分別進行多方向掃描，得到每個層段對應的傾角角度；

每個層段采用對應的傾角角度對所述混合地震數據進行中值濾波處理；

當所有層段執行完中值濾波處理后，獲得波場分離的地震數據。