可变网格与局部时间步长的高阶差分地震波数值模拟

提高计算精度与效率是所有地震波正演方法所追求的目标.本文通过将变化的空间网格与变化的时间步长技术相结合,提出一种空间网格大小与时间步长均可任意变化的高阶有限差分模拟方法.一系列数值试验表明,该方法在保证模拟精度的同时,显著提高了模拟的效率.这种可变空间网格与局部时间步长的模拟方法,能够精细刻画含孔缝洞介质以及横向变化剧烈介质的微小结构,减小地震波模拟误差,提高介质细微结构情况下的地震波传播模拟精度与效率.     

各向异性介质弹性波多参数全波形反演

各向异性介质弹性波方程全波形反演过程中多参数之间的相互耦合,使得弱参数在反演过程中难得到理想的结果.本文以VTI介质为例,在各参数辐射模式分析的基础上,基于改进的散射积分算法实现目标函数梯度的直接求取,进一步构建高斯牛顿方向,实现Hessian矩阵的有效利用,以考虑Hessian矩阵非主对角线元素包含的各参数间的耦合效应,在不使用任何反演策略的情况下实现高精度的VTI介质弹性波方程多参数同步反演.同时,该方法在计算过程中无需存储庞大的核函数矩阵,且无需传统截断牛顿法中额外的正演计算,因此内存占用小,计算效率高.本文数值试验验证了该方法的有效性,为各向异性多参数全波形反演提供了一种新的解决方案.

     

一种新的预条件伴随状态法初至波走时层析

伴随状态法初至波走时层析是基于最优化理论的一种层析成像方法,该方法不必进行射线追踪,用两次正演的计算量便可以获得梯度,具有计算效率高、内存占用小等优点.但是其一阶方向在初始模型或观测孔径不理想的情况下往往无法获得正确的反演结果,而二阶方向的实现又比较困难且费时.在伴随状态法的基础上,将走时差替换为定值,再次进行反演,便可以得到类似于射线密度的矩阵,用该矩阵的逆可以方便地进行预条件.基于该方法,本文提出了一种简单易行的预条件伴随状态法初至波走时层析的实现方法.理论模型和实际资料处理结果都表明,该方法既保留了伴随状态法初至波走时层析的优点,又可以克服一阶方向的局限,获得良好的反演效果.

     

菲涅耳体地震层析成像分辨率研究

层析成像分辨率的研究,不仅可以帮助分析层析方法的反演能力,评价层析反演的效果,还可以帮助指导层析参数设置,优化观测系统设计等.本文对比研究了前人提出的两种菲涅耳体层析成像分辨率的计算方法,并针对其存在的问题进行了优化.文中通过对二维理论模型的定量计算,总结了菲涅耳体地震层析成像分辨率的一些规律,并将其与射线层析的分辨率进行了对比.文中最后给出分辨率研究对观测系统与层析成像的指导作用,尤其是对层析中平滑窗设置、变网格层析中模型剖分原则给出了具体的建议.     

初至波菲涅尔体地震层析成像

根据地震波传播的有限频理论,对于某个特定震相的观测信息,不仅射线路径上的点对该信息具有影响,射线领域上的其他点对接收信息也具有影响,这种影响可以用核函数来表达.本文基于波动方程的Born近似与Rytov近似,给出了非均匀介质情况下初至地震波振幅与走时菲涅尔体层析成像单频、带限层析核函数的计算方法.通过对均匀介质情况下初至地震波菲涅尔体层析成像核函数解析表达式的理论模型实验与分析,给出了不同维度振幅、走时单频菲涅尔体的空间分布范围,进而给出了带限菲涅尔体边界的确定方法.将本文的走时菲涅尔体层析成像理论应用于表层速度结构反演中,理论模型试验与实际资料处理结果表明,初至波菲涅尔体地震层析成像方法比传统的初至波射线层析成像理论具有更高的反演精度.