弹性波层析成像波速反演计算方法技术

本发明专利技术涉及一种弹性波层析成像波速反演计算方法。本发明专利技术的目的是提供一种弹性波层析成像波速反演计算方法。本发明专利技术的技术方案是：一种弹性波层析成像波速反演计算方法，其特征在于：将已知波速{V

【技术实现步骤摘要】
弹性波层析成像波速反演计算方法
本专利技术涉及一种弹性波层析成像波速反演计算方法。适用于地球物理、检验检测

技术介绍
弹性波层析成像可分为两大类，一类为基于射线理论的射线层析成像和另一类为基于波动方程的波形层析成像。基于射线层析成像中的反演方法可分为线性反演方法和非线性反演方法两种，线性反演方法是将成像区划分成若干离散单元和波传播节点，根据波线性传播理论，对每条射线建立走时方程，形成线性方程组。此类方程组为一个大型的、稀疏的、常常是病态的线性方程组，一般采用优化迭代法求解，如代数重建法(ART)、联合迭代重建法(SIRT)、奇异值分解法(SVD)、共辄梯度法(CG)和阻尼最小二乘法(DLSQR法)等。迭代法是从某一初始波速开始，建立观测旅行时与理论旅行时之差和波速扰动量满足的线性方程并求解波速的扰动量，修正波速，重新进行正演，逐次迭代，直到反演结果满足某一条件为止，可以看出迭代反演的第一步就是要给定一组初始波速。初始波速的适当与否，关系到迭代的效率好坏，要使一个不适当的初始波速得到与适当初始波速相似的反演结果，可能需要更多次的本文档来自技高网...

【技术保护点】
1.一种弹性波层析成像波速反演计算方法，其特征在于：/n将已知波速{V

【技术特征摘要】
1.一种弹性波层析成像波速反演计算方法，其特征在于：
将已知波速{VR}的离散单元R作为约束条件；
对直射线经过的、单元波速未知的离散单元BPT的波速{VB}，采用反投影技术法进行单元波速的反演计算；对折线射线经过的、单元波速未知的离散单元DT的{VDT}，采用阻尼最小二乘法进行单元波速的反演计算；
再对直射线和折线射线均未经过的离散单元KDB的波速{VKDB}，用克里金插值法进行单元波速的插值计算。


2.根据权利要求1所述的弹性波层析成像波速反演计算方法，其特征在于：采用反投影技术法计算前，先扣除各条直射线经过离散单元R的走时和线段长度，再用反投影技术法进行反演计算，计算各条直射线经过的离散单元BPT的波速{VB}。


3.根据权利要求1所述的弹性波层析成像波速反演计算方法，其特征在于：所述折线射线利用Dijkstra算法在离散单元波速场上搜索各条射线的最短走时路径而形成；
所述离散单元波速场利用离散单元R的波速{VR}、离散单元BPT的波速{VB}和离散单元KB的波速{VKB}建立；
所述离散单元KB为波速未知的离散单元U中直射线未经过的单元；离散单元KB的波速{VKB}利用离散单元R的波速{VR}和离散单元BPT的波速{VB}，通过克里金插值法计算获得。


4.根据权利要求1所述的弹性波层析成像波速反演计算方法，其特征在于：采用阻尼最小二乘法计算前，先扣除各条折线射线经过离散单元R的走时，建立走时线性方程组，再用阻尼最小二乘法反演计算获得离散单元DT的{VDT}。


5.根据权利要求1所述的弹性波层析成像波速反演计算方法，其特征在于：所述离散单元KDB的波速{VKDB}利用离散单元R的波速{VR}、离散单元BPT的波速{VB}、离散单元DT的波速{VDT}和离散单元DB的波速{VDB}，通过克里金插值法计算获得；
所述离散单元DB为波速未知的离散单元U中直射线和折线射线均经过的单元，离...

【专利技术属性】
技术研发人员：陈文华，李涛，王飞，曾超，
申请(专利权)人：杭州国家水电站大坝安全和应急工程技术中心有限公司，中国电建集团华东勘测设计研究院有限公司，
类型：发明
国别省市：浙江;33