多轴完全匹配层的非分裂实现

完全匹配层吸收边界(PML)已经被证明是非常有效的边界吸收技术,对体波和面波的吸收都具有非常好的效果,已经被广泛应用于弹性波的数值模拟中.但是在某些情况下传统的PML和C-PML(卷积完全匹配层)还是存在一定的问题,比如在一些介质中有固有的不稳定性等等.基于分裂的多轴完全匹配层(M-PML)通过引入多个衰减剖面可以有效地消除不稳定的现象.本文采用C-PML中不分裂的递推卷积方法,对M-PML进行了实现,这样在实现过程中不用分裂变量,与原来的分裂方法相比,直观易懂,方便简单,易于编程,而且可以减少计算量和存储量,节约计算成本.基于一阶速度-应力方程,采用交错网格高阶有限差分法对不同的介质模型进行了数值模拟,结果证明了本文的不分裂M-PML可以消除不稳定性,得到比较好的计算结果.

Lebedev网格改进差分系数TTI介质正演模拟方法研究

本文采用一种新的交错网格-Lebedev网格(LG)进行TTI介质的正演模拟研究,避免了Virieux标准交错网格(SSG)算法在处理TTI、单斜等各向异性介质时波场插值引入的数值误差,提高了模拟精度.在方法实现过程中,本文针对有限差分正演模拟面临的网格频散与边界反射两个关键性问题分别做了优化,并通过模型试算验证了它们的有效性与可行性:(1)结合最小二乘思想推导出新的频散改进差分系数(DIC),该系数比Taylor系数更能有效地压制粗网格引起的数值频散,可以节约内存,提高计算效率;(2)将分裂的多轴完全匹配层(M-PML)吸收边界条件引入到LG算法中,解决了传统PML边界条件在某些各向异性介质中的不稳定现象并且具有较好的边界吸收效果.

孔隙介质与液体界面上波动方程数值模拟的改进方案——MPML边界条件的应用

由于波场数值模拟的计算区域是有限的,为了压制截断边界造成的人工边界反射,通常采用完全匹配层作为吸收边界条件。但是在含液-固界面的弹性介质中进行正演模拟时,完全匹配层边界条件容易造成稳定性方面的问题。因此,该文将多轴完全匹配层应用于该类介质的正演模拟,以改善完全匹配层边界条件引起的不稳定问题。通过在具有弹性海底的海洋环境模型以及充填液体的井孔模型中进行弹性波方程的正演模拟,分别采用传统的完全匹配层和多轴完全匹配层作为吸收边界条件,对比验证了多轴完全匹配层在含液-固界面弹性介质中数值模拟的有效性和稳定性。

基于Lebedev网格的TTI介质二维三分量正演模拟

传统交错网格有限差分法是研究地震波传播规律的一种较为常用的数值模拟方法,但是在交错网格中每个变量的不同分量都是交错定义的,对于没有定义的点需要进行变量插值,从而降低了模拟精度。为此,在前人的研究基础上,推导了TTI介质二维三分量的应力—速度方程,采用Lebedev网格对其进行了高精度的差分离散处理,避免了传统交错网格在处理各向异性介质时波场插值引起的误差,提高了模拟精度,并将多轴完全匹配层吸收边界(M-PML)引入了Lebedev网格。分别对单层TTI介质和含透镜体的复杂TTI介质模型进行了正演模拟,结果表明:(1)由于Lebedev网格对各向异性介质的弹性波方程做离散时无需进行波场插值,与传统交错网格有限差分方法相比,模拟精度更高;对TTI介质进行模拟时可以清晰地观察到纵波、快横波和慢横波,并且快、慢横波的偏振方向相反,在单炮记录中观测到的三种波的速度特征也符合波场传播规律。(2)引入多轴完全匹配层(M-PML)吸收边界条件后,在不影响模拟效果的情况下边界反射现象被有效地压制。

非水平海底情况下海底地震波时域有限差分数值模拟

研究复杂海洋环境中海底地震波激发及其传播特性,对海底物理力学特性研究、资源勘探等具有重要的意义.目前针对时域海底地震波的研究大都局限于水平分层的情况,而实际的海底地质条件比较复杂,基于理想环境假设得出的数值解与实际差别较大.本文在考虑倾斜、隆起等非水平海底模型的情形下,采用时间2阶精度、空间10阶精度的交错网格有限差分方法,同时结合多轴完全匹配层边界条件,对复杂海洋环境下的海底地震波进行时域数值模拟与分析.利用计算得到的声场时域波形,分析了复杂海洋环境下海底地震波的传播特性.结果表明,采用空间高阶精度的交错网格有限差分方法,可改善数值计算中的频散问题;同时采用多轴完全匹配层替代传统的完全匹配层,解决了液-固介质中远距离声场数值模拟不稳定的问题.在含倾斜与隆起构造的复杂海底模型中,海底基岩隆起改变了Scholte波的传播方向,更有利于在较浅深度处接收到Scholte波.

基于M-PML边界的Lebedev网格起伏地表正演模拟方法及稳定性分析

贴体网格作为一种'边界拟合网格',能改善传统有限差分处理起伏地表时由阶梯状网格存在所产生的虚假绕射。此外,完全匹配层(PML)技术被证明是一种能够有效消除虚假反射的重要边界条件,但采用自由地表模拟面波时,传统PML易产生不稳定现象。为实现起伏地表自由边界条件下稳定的地震波场模拟,发展一种基于弹性介质的全交错-Lebedev网格有限差分方法,并将边界条件扩展为多轴完全匹配层(M-PML)。在实现算法的基础上,通过模型试算验证新方法的正确性;并对比几种不同完全匹配层技术的稳定性差异,解释PML边界产生不稳定的机制。正演模拟结果表明,M-PML在处理起伏地表面波模拟中有更好的稳定性。

Lebedev网格黏弹性介质起伏地表正演模拟

在前人研究的基础上,采用一种新的交错网格(Lebedev网格)进行曲坐标系下的黏弹性介质正演模拟,避免了标准交错网格在处理曲坐标系方程时进行波场插值而引入的数值误差,从而提高了模拟精度。在正演模拟的过程中,首先基于广义标准线性固体,推导了曲坐标系下黏弹性介质的波动方程,随后利用在各向异性介质中使用的Lebedev网格有限差分方法对波动方程进行了离散化,在地表附近采用牵引力镜像法来实施自由表面条件,其他三个边界引进多轴卷积完全匹配层技术提高吸收效果,最后通过模型试算分析了黏滞性的引入以及地形起伏对波场的双重影响,并验证了引入的边界条件具有较好的吸收效果。模型试算结果表明,由于黏弹性介质中吸收衰减的影响导致地震波能量降低且主频向低频端移动,同时由于速度频散导致走时差异及波形变化。