深度均匀采样梯形网格有限差分地震波场模拟方法

由于重力引起的岩石压实效应,一般来说,地震波传播速度由浅入深整体逐渐增大.梯形坐标系设计可耦合速度由浅入深逐渐增大的变化,该坐标系中均匀网格采样所对应的物理直角坐标系网格由浅入深逐渐增大,也即浅部低速区对应细网格,深部高速区对应粗网格.在梯形坐标系表征波动方程后利用有限差分求解,本文实现一种深度均匀采样、横向采样间隔随深度增加逐渐线性增大的有限差分地震波模拟方法.梯形坐标系波动方程离散后,仍采用常规均匀网格有限差分算法对其求解.由于横向网格大小由浅入深线性增加,本方法可避免不同大小网格区域过渡所产生的虚假反射.梯形坐标系波场模拟浅层精度高,深层横向响应范围广,可有效减少有限差分网格数量.本文提出的方法是在更广义的坐标系下利用有限差分求解波动方程,正交坐标系仅为该梯形坐标系之特例.本文旨在为大速度动态范围深地高效高精度地震波场模拟提供一种思路.

     

基于双变网格算法的地震波正演模拟

为了适应对局部复杂模型的精细模拟,本文实现了可变网格算法,对速度场进行局部加密,从空间上有效地提高模拟精度同时又降低计算机内存需求.但是在数值模拟中,由于稳定性条件的限制,当空间网格变化时,时间稳定性仍然必须满足最短波长的原则,从而增加了时间计算量.为了配合空间可变网格技术,本文对时间层计算也进行了局部变化,提出了双变网格有限差分算法.双变算法从根本上降低了整体的计算量和计算内存.本文还通过一系列的模型试算验证了双变网格算法的稳定性和正确性,并对双变算法进行了误差分析,证明了双变网格算法的精确性.     

二维弹性及粘弹性TTI介质中地震波场数值模拟:四种不同网格高阶有限差分算法研究

本文利用交错网格、辅助网格、旋转交错网格、同位网格有限差分方法分别模拟了二维弹性TTI介质和二维黏弹性TTI介质中的地震波传播.在稳定性条件内,选用不同的网格间距及时间间隔,通过波场快照、合成理论地震图较为系统分析对比了这四种不同网格有限差分数值模拟在计算精度、CPU时间、相移、频散、以及保幅方面的优缺点.数值模拟结果表明:1)这四种不同网格有限差分算法都是很好的波场数值模拟算法;2)就CPU计算时间而言,旋转交错网格有限差分算法的计算效率最高;3)从计算精度来看,同位网格有限差分的计算精度最高;4)从振幅保护方面来看,四种网格的保护振幅的能力相当;5)相移方面,当网格间距增大时,交错网格和旋转交错网格有可能出现相移现象;6)频散方面,同位网格的频散现象不明显.     

旋转时空双变网格有限差分地震波场正演模拟

在交错网格有限差分算法中,模型网格剖分原则与正演计算效率密切相关.当模型存在小型非均质体或者低降速层等情况,为保证精度,满足稳定性条件,需缩小网格步长,导致局部过采样,计算效率低下.为保证模拟精度的同时保持高计算效率,通常采用变化的空间网格与时间步长相结合的高阶有限差分模拟方法对波场进行模拟.然而,时空双变算法存在着交错网格固有缺点,在模拟非均匀性较强的复杂介质波场传播时,需对介质参数进行平均或内插.同时,该算法在空间与时间上的变网格实现极为复杂.为压制变网格引起的虚假反射,提升模拟精度和计算效率,本文在时空双变网格算法的基础上,采用旋转差分角度的方式,提出了旋转时空双变交错网格算法.该方法既保留了旋转网格和双变网格的优势,又简化了时空双变算法流程,更利于推广和应用.

     

基于变换光学有限差分探地雷达数值模拟研究

在采用有限差分方法开展探地雷达复杂目标体精细结构模拟时,为了提高计算精度,常采用非均匀网格对目标区域划分小尺寸的网格,以压制离散网格频散现象和保证有限差分方法的稳定性.常规非均匀网格和自适应亚网格技术在网格剖分数量和粗细网格边界处理上难以达到计算效率和计算精度的均衡.本文根据隐形斗篷(invisible cloak)理论,将基于变换光学(Transformation optics)理论应用于有限差分探地雷达数值计算中.该理论的主要思想是基于目标参数变化而保持电磁场的传播不变性,在坐标变换后,Maxwell方程的形式可以维持不变,而使得相对介电常数与磁导率的表达式变得复杂.通过这种方式可以虚拟地扩大目标体所占的网格节点数,减少背景介质区域的网格数,不增加模型空间的网格总数.另外,这种网格划分方式不但提高了计算效率,同时也可以克服亚网格技术边界反射误差的影响.本文推导实现了基于变换光学的二维有限差分方法,通过典型探地雷达模型测试,对比分析了该方法与常规有限差分、变网格有限差分和自适应亚网格有限差分的优缺点.计算结果验证了基于变换光学的有限差分可用于探地雷达目标精细结构模拟,具有较高的计算精度和计算效率.     

交错网格与旋转交错网格对VTI介质波场分离的影响分析

利用交错网格有限差分和旋转交错网格有限差分进行各向异性介质弹性波场数值模拟时, 质点振动速度分量与应力张量的网格节点定义方式均不相同, 从而对各向异性波场分离效果产生不同的影响. 针对这一问题, 本文以具有垂直对称轴的横向各向同性(VTI)介质的波场分离为例, 首先分析了两种网格的参数定义方式以及VTI介质波场的分离过程; 其次, 详细研究和分析了这两种网格的参数定义方式对各向异性介质波场分离的影响, 并依据波前面连续性以及波场分离效果等方面, 通过数值模拟实验对该影响进行分析验证. 结果表明, 旋转交错网格的参数定义方式更有利于进行各向异性介质波场数值模拟和波场分离.      

用于声波正演模拟的时空域高精度交错网格有限差分方法

地震正演模拟是逆时偏移和全波形反演中的核心问题之一,因为它们都需要高效、高精度地模拟波场正向和反向传播。为了提高数值模拟的精度,人们广泛采用高阶有限差分方法,但是大多数方法仅在空间上具有更高的精度,在时间上只有二阶精度。首先系统介绍时空域高精度交错网格有限差分方法的基本原理,然后利用模型验证方法的有效性,结果表明:时空域高精度交错网格有限差分方法拥有比常规交错网格有限差分方法更低的数值频散。      

基于PML边界条件的高倍可变网格有限差分数值模拟方法

非均匀介质一般为多尺度介质,对非均匀介质进行正演模拟需要多尺度的网格剖分.碳酸盐岩缝洞介质的尺度一般为厘米级甚至毫米级.使用有限差分方法对其进行精细模拟需要差分步长达到缝洞介质的尺度.为了提高有限差分数值模拟方法的精度和效率,使之可以应用于非均匀介质的正演模拟,本文推导了基于PML边界的空间和时间步长同时变化的高倍数可变网格差分格式,步长变化倍数可以达到百倍以上.并且在一般意义的变网格算法的基础上,改进了变网格算法的网格剖分方式,进一步减小了精细尺度模型数值模拟的内存消耗.数值试验表明,该方法可以精细描述毫米尺度的地质体,提高有限差分方法模拟精度,同时也节约了内存,提高了模拟效率.     

高阶交错网格非稳态相移叠后深度偏移

为更好地研究复杂构造和速度分布条件下地震资料的精确偏移.在波场模拟方面,本文采用了高阶交错网格有限差分[1]推导了二维各向同性介质声波方程的数值模拟公式.并给出了声波方程完全匹配层吸收边界[2-5](PML).在偏移方面,本文介绍了一种能适应介质速度横向变化的非稳态相移算子(NSPS)[6]及其叠后深度偏移方法.并对简单的四层模型和Marmousi模型进行试算.偏移结果表明,该方法可以有效地处理复杂条件下的各向同性介质声波偏移问题.     

一种全局优化的隐式交错网格有限差分算法及其在弹性波数值模拟中的应用

在数值模拟中,隐式有限差分具有较高的精度和稳定性.然而,传统隐式有限差分算法大多由于需要求解大型矩阵方程而存在计算效率偏低的局限性.本文针对一阶速度-应力弹性波方程,构建了一种优化隐式交错网格有限差分格式,然后将改进格式由时间-空间域转换为时间-波数域,利用二范数原理建立目标函数,再利用模拟退火法求取优化系数.通过对均匀模型以及复杂介质模型进行一阶速度-应力弹性波方程数值模拟所得单炮记录、波场快照分析表明:这种优化隐式交错网格差分算法与传统的几种显式和隐式交错网格有限差分算法相比不但降低了计算量,而且能有效的压制网格频散,使弹性波数值模拟的精度得到有效的提高.