[發明專利]一種基于動態時間規整的地震數據加權疊加方法

有效

• 申請號: 202010529953.2
• 申請日: 2020-06-11
• 公開（公告）號: CN111610566B
• 公開（公告）日: 2022-11-29
• 發明（設計）人: 宋承云;李凌軒;陳韻清;武子騰
• 申請（專利權）人: 重慶理工大學
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 重慶信航知識產權代理有限公司 50218
• 代理人: 穆祥維

1.一種基于動態時間規整的地震數據加權疊加方法，其特征在于，包括以下步驟：

S1、選擇一個參考地震道，該參考地震道設置為疊前地震數據矩陣在空間方向上的算術均值，即：

其中，A_i(t)為第i個疊前地震道A_i在時間t時的采樣值，C為疊前地震道數目，R(t)為參考地震道R在時間t時的采樣值；

S2、通過下式計算每個疊前地震道與參考地震道之間的動態時間規整即DTW匹配點，在CMP道集中尋找與參考地震道的各個采樣點相關聯的采樣點：

其中，γ_k是第k個點對p_k＝(r_n,a_im)的差異，該差異表示連接參考地震道R＝[r₁,r₂,…,r_n,…r_N]與疊前地震道A_i＝[a_i1,a_i2,…,a_im,…,a_iM]兩個道中采樣點的連接線，具體用(r_n-a_im)²來計算，索引n并不一定等于m，且每個點對代表了在以下約束下的一條對齊路徑：

單調性：對齊路徑p_k＝(r_n,a_im)和p_k-1＝(r_n',a_im')的任意兩個相鄰元素，滿足不等式n-n'≥0和m-m'≥0；

連續性：對齊路徑p_k＝(r_n,a_im)和p_k-1＝(r_n',a_im')的任意兩個相鄰元素，滿足不等式n-n'≤1和m-m'≤1；

邊界：第一個對齊路徑是p₁＝(r₁,a_i1)，最后一個對齊路徑是p_k＝(r_N,a_iM)；

S3、根據各相關采樣點的噪聲設置各樣點的權值：

w_p＝1/(d_p)^s

d_p＝|R(t)-X_p|

其中，w_p為X_p的權重，d_p為X_p與參考地震道R(t)之間的距離，t為時間，p＝1，2，...，P，P為關聯點數目，X_p為疊前地震道中與參考地震道中采樣點相關聯的點，s為用于控制權重函數的用戶定義常量；

S4、將參考地震道的每個采樣點更新為與之相關的采樣點的加權平均值：

其中，R(t)即為由加權平均值更新后的參考地震道；

S5、重復步驟S2～S4直到收斂，將參考地震道輸出為最終的疊加地震道。

2.根據權利要求1所述的基于動態時間規整的地震數據加權疊加方法，其特征在于，所述步驟S3中，用于控制權重函數的用戶定義常量s取值為0～2。