地震波在地下介质传播过程中由于非弹性衰减的存在将导致能量损失和相位变化, 精确的速度与衰减参数建模对油气识别、提高强衰减介质中地震波成像的质量都起着至关重要的作用.常分数阶拉普拉斯算子黏声方程由于完全分离的速度频散项与振幅衰减项的优势, 以及在强非均质衰减介质中可以高精度求解的特点, 已被应用于速度与衰减参数的建模中.本文将二阶常分数阶拉普拉斯算子黏声方程拆分为等价的一阶方程组, 并在此一阶方程组的基础上推导出新的梯度公式与伴随方程, 建立了一种新的速度与衰减参数同时重建的全波形反演方法.相较于原二阶常分数阶拉普拉斯算子黏声方程建立的全波形反演流程, 数值实验表明, 新建立的反演流程可以有效避免原梯度数值计算中的噪声, 尤其是可以有效提高衰减参数梯度的反演精度, 从而显著提高反演的收敛速度与反演精度.
相似文献伴随状态法初至波走时层析是基于最优化理论的一种层析成像方法,该方法不必进行射线追踪,用两次正演的计算量便可以获得梯度,具有计算效率高、内存占用小等优点.但是其一阶方向在初始模型或观测孔径不理想的情况下往往无法获得正确的反演结果,而二阶方向的实现又比较困难且费时.在伴随状态法的基础上,将走时差替换为定值,再次进行反演,便可以得到类似于射线密度的矩阵,用该矩阵的逆可以方便地进行预条件.基于该方法,本文提出了一种简单易行的预条件伴随状态法初至波走时层析的实现方法.理论模型和实际资料处理结果都表明,该方法既保留了伴随状态法初至波走时层析的优点,又可以克服一阶方向的局限,获得良好的反演效果.
相似文献各向异性介质弹性波方程全波形反演过程中多参数之间的相互耦合,使得弱参数在反演过程中难得到理想的结果.本文以VTI介质为例,在各参数辐射模式分析的基础上,基于改进的散射积分算法实现目标函数梯度的直接求取,进一步构建高斯牛顿方向,实现Hessian矩阵的有效利用,以考虑Hessian矩阵非主对角线元素包含的各参数间的耦合效应,在不使用任何反演策略的情况下实现高精度的VTI介质弹性波方程多参数同步反演.同时,该方法在计算过程中无需存储庞大的核函数矩阵,且无需传统截断牛顿法中额外的正演计算,因此内存占用小,计算效率高.本文数值试验验证了该方法的有效性,为各向异性多参数全波形反演提供了一种新的解决方案.
相似文献