[發明專利]彈性波層析成像波速反演計算方法

有效

• 申請號: 202010941185.1
• 申請日: 2020-09-09
• 公開（公告）號: CN112099084B
• 公開（公告）日: 2022-05-10
• 發明（設計）人: 陳文華;李濤;王飛;曾超
• 申請（專利權）人: 杭州國家水電站大壩安全和應急工程技術中心有限公司;中國電建集團華東勘測設計研究院有限公司
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 杭州九洲專利事務所有限公司 33101
• 代理人: 韓小燕

1.一種彈性波層析成像波速反演計算方法，其特征在于：

將已知波速{V_R}的離散單元R作為約束條件；

對直射線經過的、單元波速未知的離散單元BPT的波速{V_B}，采用反投影技術法進行單元波速的反演計算；對折線射線經過的、單元波速未知的離散單元DT的{V_DT}，采用阻尼最小二乘法進行單元波速的反演計算；

再對直射線和折線射線均未經過的離散單元KDB的波速{V_KDB}，用克里金插值法進行單元波速的插值計算。

2.根據權利要求1所述的彈性波層析成像波速反演計算方法，其特征在于：采用反投影技術法計算前，先扣除各條直射線經過離散單元R的走時和線段長度，再用反投影技術法進行反演計算，計算各條直射線經過的離散單元BPT的波速{V_B}。

3.根據權利要求1所述的彈性波層析成像波速反演計算方法，其特征在于：所述折線射線利用Dijkstra算法在離散單元波速場上搜索各條射線的最短走時路徑而形成；

所述離散單元波速場利用離散單元R的波速{V_R}、離散單元BPT的波速{V_B}和離散單元KB的波速{V_KB}建立；

所述離散單元KB為波速未知的離散單元U中直射線未經過的單元；離散單元KB的波速{V_KB}利用離散單元R的波速{V_R}和離散單元BPT的波速{V_B}，通過克里金插值法計算獲得。

4.根據權利要求1所述的彈性波層析成像波速反演計算方法，其特征在于：采用阻尼最小二乘法計算前，先扣除各條折線射線經過離散單元R的走時，建立走時線性方程組，再用阻尼最小二乘法反演計算獲得離散單元DT的{V_DT}。

5.根據權利要求1所述的彈性波層析成像波速反演計算方法，其特征在于：所述離散單元KDB的波速{V_KDB}利用離散單元R的波速{V_R}、離散單元BPT的波速{V_B}、離散單元DT的波速{V_DT}和離散單元DB的波速{V_DB}，通過克里金插值法計算獲得；

所述離散單元DB為波速未知的離散單元U中直射線和折線射線均經過的單元，離散單元DB的波速{V_DB}取離散單元DT的波速{V_DT}。