[發明專利]基于空間信息索引的多道匹配追蹤方法、存儲介質及終端

有效

• 申請號: 201910252884.2
• 申請日: 2019-03-29
• 公開（公告）號: CN110082826B
• 公開（公告）日: 2020-10-27
• 發明（設計）人: 張向濤;陳維濤;劉道理;宗兆云;印興耀;許璐;李坤
• 申請（專利權）人: 中國海洋石油集團有限公司;中海石油（中國）有限公司深圳分公司;中國石油大學（華東）
• 主分類號: G01V1/48
• 分類號: G01V1/48;G01V1/50;G06F16/901;G06F16/903
• 代理公司: 深圳市瑞方達知識產權事務所(普通合伙) 44314
• 代理人: 張約宗;張秋紅

1.一種基于空間信息索引的多道匹配追蹤方法，其特征在于，包括：

處理獲取的多道信號得到最佳匹配原子以及最佳匹配原子的振幅參數；

構建空間時變原子庫用于對鄰道進行分解，所述空間時變原子庫中時間參數為變化量；所述空間時變原子庫中瞬時頻率和瞬時相位為固定量，且與所述最佳匹配原子相同；

計算所述空間時變原子庫確定預設道信號的最佳原子，由所述最佳原子的中心時間替換所述空間時變原子庫中時間參數，替換后的空間時變原子庫繼續對下一道信號進行多道搜索；

所述處理獲取的多道信號得到最佳匹配原子以及最佳匹配原子的振幅參數包括：

對獲取的多道信號的第i道信號進行希爾伯特變換，得到瞬時頻率、瞬時相位以及包絡振幅信息；

對信號的最大包絡振幅位置進行分解進而構建動態原子庫，搜索所述動態原子庫得到最佳匹配原子；

計算所述最佳匹配原子的振幅參數；

所述對信號的最大包絡振幅位置進行分解進而構建動態原子庫包括：

所述動態原子庫D＝{g_γ(t)}_γ∈Γ用下式表示：

其中，式中γ為原子g(t)的控制參數，t為時間，i為當前被分解的地震信號的道號，u_i為原子的中心時間即信號的最大包絡振幅所在的時間位置，n＝1:N，N為每次迭代得到的最優時頻原子的個數；ω和表示地震信號瞬時頻率及瞬時相位的先驗信息；U[ω(u_i),δ_ω]、為頻率和相位的搜索鄰域；δ_ω、為相應參數搜索半徑；

所述搜索所述動態原子庫得到最佳匹配原子包括：

所述最佳匹配原子滿足：

其中，f_i(t)代表第i道地震信號，為內積運算符；

所述最佳匹配原子的振幅參數為：

所述空間時變原子庫

表示為：

其中，W為搜索時窗，決定所述空間時變原子庫的大小；△u為采樣間隔，n＝1:N,N為原子的個數，頻率ω和相位與式(2)中最佳匹配原子的參數大小相同；

所述計算所述空間時變原子庫確定預設道信號的最佳原子包括：

所述計算所述空間時變原子庫中所有原子與第i+1道信號f_i+1(t)的振幅索引A_n＝1:1:M，最大振幅sup A_n＝1:1:M所指示的原子為第i+1道信號的最佳原子；

其中，為第i+1道信號第一次多道搜索確定的最佳原子，a_i+1為所述最佳原子的振幅參數。

2.根據權利要求1所述的基于空間信息索引的多道匹配追蹤方法，其特征在于，在所述替換后的空間時變原子庫繼續對下一道信號進行多道搜索之后，所述方法還包括：

對所有道進行第一次多道搜索的振幅索引滿足式(6)時，停止第一次多道搜索，開始對第i道信號進行第二次多道分解；

sup A_k+1≤α sup A_k (6)

式中，0α≤1為多道搜索閾值，A_k和A_k+1分別為第一次多道搜索中第k道和第k+1道地震信號最佳原子的振幅索引值。

3.一種計算機可讀存儲介質，其上存儲有計算機程序，其特征在于，所述計算機程序被處理器執行時實現如權利要求1-2中任意一項所述的基于空間信息索引的多道匹配追蹤方法。

4.一種終端，其特征在于，所述終端包括處理器，所述處理器用于執行存儲器中存儲的計算機程序時實現如權利要求1-2中任意一項所述基于空間信息索引的多道匹配追蹤方法的步驟。