基于Remez迭代算法求解差分系数的双相介质自适应有限差分正演模拟

利用传统有限差分方法对基于Biot理论的双相介质波动方程进行数值求解时,由于慢纵波的存在,数值频散效应较为明显,影响模拟精度.相对于声学近似方程及普通弹性波方程,Biot双相介质波动方程在同等数值求解算法和精度要求条件下,其地震波场正演模拟需要更多的计算时间.本文针对Biot一阶速度-应力方程组发展了一种变阶数优化有限差分数值模拟方法,旨在同时提高其正演模拟的精度和效率.首先结合交错网格差分格式推导Biot方程的数值频散关系式.然后基于Remez迭代算法求取一阶空间偏导数的优化差分系数,并用于Biot方程的交错网格有限差分数值模拟.在此基础上把三类波的平均频散误差参数限制在给定的频散误差阈值和频率范围内,此时优化有限差分算子的长度就能自适应非均匀双相介质模型中的不同速度区间.数值频散曲线分析表明:基于Remez迭代算法的优化有限差分方法相较传统泰勒级数展开方法在大波数范围对频散误差的压制效果更明显;可变阶数的优化有限差分方法能取得与固定阶数优化有限差分方法相近的模拟精度.在均匀介质和河道模型的数值模拟实验中将本文变阶数优化有限差分算法与传统泰勒展开算法、最小二乘优化算法进行比较,进一步证明其在复杂地下介质中的有效性和适用性.     

旋转网格和常规网格混合的时空域声波有限差分正演

压制数值频散,提高正演模拟精度,一直是有限差分正演模拟研究的重要内容.基于时空域频散关系的有限差分法,比基于空间域频散关系的传统有限差分法,模拟精度更高.时空域声波方程数值模拟,普遍采用常规十字交叉型高阶有限差分格式.而在频率-空间域,普遍采用旋转网格和常规网格混合的有限差分格式,有效提高了模拟精度和计算效率.本文将频率-空间域混合网格有限差分的思想引入到时空域,提出了时空域混合网格2 M+N型声波方程有限差分方法.首先推导出基于时空域频散关系的混合网格差分系数计算方法,然后进行频散分析、稳定性分析,并和传统高阶、时空域高阶有限差分法对比,结果表明:计算量相同时,新方法能有效压制数值频散,显著提高模拟精度;新方法相比传统2 M阶有限差分法,稳定性增强,与时空域2 M阶有限差分法稳定性基本相当.最后利用新方法进行均匀介质、层状介质、盐丘模型的数值模拟和盐丘模型的逆时偏移,模拟效果和成像质量进一步证实了该方法的有效性和普遍适用性.     

基于纵横波分离FCT弹性波正演频散压制

有限差分方法被广泛应用于地震波数值模拟和传播.传统有限差分法采用Taylor级数展开实现空间偏导数的差分,但该方法会因为网格离散化而产生数值频散,降低地震波模拟的精度.优化差分系数正演方法能在一定程度上压制部分频散,然而纵、横波速度取值差异较大,在弹性波有限差分正演模拟中,在满足纵波最大速度确定的稳定性条件下,浅层低速横波波场往往会产生明显的频散现象.为了削弱弹性波场正演数值频散,提高数值模拟精度,本文首先采用优化差分网格系数降低数值频散,然后再采用通量校正传输(Flux-Correction Transport, FCT)法来进一步压制弹性波场有限差分数值频散.常规的FCT法是对弹性波场直接进行频散压制,但由于弹性波场中纵、横波速度差异明显,横波波场频散明显强于纵波,为了压制横波波场的数值频散,往往需要选取较大的频散压制参数,但这会使频散较弱的纵波产生假象.因此本文提出基于纵横波分离FCT弹性波正演频散压制方法,对分离之后的纵横波场分别选择合适的频散压制参数进行通量校正,可以有效压制数值频散,削弱纵波FCT产生的假象.通过理论分析和数值算例发现,本文方法能有效削弱弹性波场有限差分数值...     

横向各向同性介质地震波场数值模拟研究

地震波场数值模拟是理解地震波在地下介质中的传播特点,帮助解释观测数据的有效手段,而提高计算精度和运算效率是所有波场数值模拟方法研究所追求的目标.有限差分技术是求解波动方程计算效率最高、应用最为广泛的方法之一.但传统的有限差分技术计算过程中的数值频散问题影响了该技术的计算精度与计算效率.本文通过交错网格高阶有限差分技术与通量校正传输方法（Flux|corrected transport method,FCT）相结合, 对横向各向同性介质（Transverse isotropic medium,TI）一阶速度|应力弹性波动方程组进行了数值求解研究.波场快照数值模拟结果表明,本文研究的数值模拟方法与波动方程二阶有限差分方法、交错网格四阶有限差分方法相比,在压制网格数值频散方面有明显的优势,计算精度提高,而且可以利用较大的空间步长,提高计算效率.     

波场模拟有限差分参数优化选取

有限差分方法(Finite-difference Method, FD)广泛用于地震波场数值模拟，但其存在固有的数值频散问题，影响模拟的计算效率和数值精度.本文主要研究了有限差分方法的空间数值频散误差和网格划分精度以及差分算子的关系，基于计算量最小准则，提出了最优化有限差分参数选取流程，为有限差分数值模拟参数选取提供理论指导.本文主要工作包括：(1)提出了空间数值频散正变换过程(Forward Space Dispersion Transform, FSDT)方法，该方法可以高效模拟出不同网格划分精度(波长采样点数)的带有空间数值频散的波场；(2)提出了波场空间数值频散误差衡量准则，可以定量地判断出数值模拟导致的波形频散程度，选取合适的频散误差阈值；(3)研究了给定空间数值频散误差阈值下，差分算子系数、差分算子阶数、网格划分精度与计算量之间的关系.文中基于雷米兹交换方法(Remez Exchange Method, RE)和泰勒级数展开方法(Taylor-series Expansion Method, TE)的差分系数，在空间数值频散误差阈值0.01时，数值模拟了不同差分算子阶数、网...     

用于弹性波方程数值模拟的有限差分系数确定方法

有限差分方法是波场数值模拟的一个重要方法,交错网格差分格式比规则网格差分格式稳定性更好,但方法本身都存在因网格化而形成的数值频散效应,这会降低波场模拟的精度与分辨率.为了缓解有限差分算子的数值频散效应,精确求解空间偏导数,本文把求解波动方程的线性化方法推广到用于求解弹性波方程交错网格有限差分系数;同时应用最大最小准则作为模拟退火(SA)优化算法求解差分系数的数值频散误差判定标准来求解有限差分系数.通过上述两种方法,分别利用均匀各向同性介质和复杂构造模型进行了数值正演模拟和数值频散分析,并与传统泰勒展开算法、最小二乘算法进行比较,验证了线性化方法和模拟退火方法都能有效压制数值频散,并比较了各个算法的特点.     

一种基于全局优化的交错网格有限差分法

在地震波场数值模拟中, 交错网格有限差分技术得到了广泛的应用, 但是在弹性模量变化较大时, 通常会因插值而导致模拟误差增大. 旋转交错网格可以很好地克服这个缺点, 因而适合于各向异性介质正演模拟. 但是对于同样大小的网格单元, 旋转交错网格需要的步长比常规交错网格要大, 这会使梯度和散度算子的误差增大因而更易产生空间数值频散. 针对这些问题, 本文提出了旋转交错网格与紧致有限差分相结合的方法, 并基于模拟退火算法进行全局优化, 压制数值频散, 拓宽波数范围. 数值模拟结果表明, 此方法可以有效地压制数值频散, 且具有较高的模拟精度.      

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

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

基于自适应网格的仿真型有限差分地震波数值模拟

在复杂山地和复杂海底条件下,地表和海底的剧烈起伏对地震波数值模拟提出了更高的要求.常规有限差分法采用矩形网格对模型进行网格剖分,由于矩形网格自身的限制,起伏地表或起伏海底只能由一系列阶梯状折线代替,从而引起人为虚假绕射波.此外,在模拟液-固界面的反射波时,如果界面与网格线不一致,则需要更密的网格才能得到精确的结果.为了解决上述问题,本文将自适应网格生成技术引入到起伏海底速度模型的网格剖分中,采用高阶仿真型有限差分法(MFD)对曲线坐标下的声波方程波进行了数值模拟.利用自适应网格生成技术对速度模型进行网格剖分不仅可以准确地描述模型边界,而且可以有效消除虚假绕射波.高阶仿真型有限差分法可以有效压制频散提高计算精度.模型试算结果表明,本文方法对复杂海底模型具有很好的适应性.     

常Q衰减介质分数阶波动方程优化有限差分模拟

本文基于Kjartansson常Q模型理论,推导了常Q衰减介质中黏声波和黏弹性波的速度-应力方程,并采用基于二项式窗函数的优化交错网格有限差分方法进行了数值模拟,同时引入不分裂的复频移卷积完全匹配层(CPML)吸收边界条件,以消除边界反射.使用基于自适应时间步长记忆方法的中心差分近似时间分数阶导数,与常用的短时记忆方法相比,提高了波动方程的离散化精度和计算效率.通过对比均匀模型下声波的数值解与解析解,验证了算法的精确性,并进一步分析了不同品质因子下地震波的频散及衰减特征.对BP盐丘模型的数值模拟结果可以较好地反映本文数值方法对复杂介质的适应性及频散压制效果.      

弹性波变阶数旋转交错网格数值模拟（英文）

地震波场数值模拟方法对理解和分析地震波的传播规律具有着重要的意义。弹性波动方程能够模拟地下介质的实际情况,为偏移和成像提供有效的依据。在弹性波波场数值模拟中,旋转交错网格数值模拟(RSM)修改了标准交错网格数值模拟(SSM)方法,将同类的参数定义在同样的节点上,拓宽了稳定性条件的约束,但在低速区会出现较严重的频散。变阶数差分方法是自适应空间算子长度方法的一种变化和推广。它以理论频散误差研究为基础,结合实际波场传播的情况进行误差计算,对不同速度匹配不同的差分阶数。本文研究了变阶数旋转交错网格数值模拟(VRSM),即是将变阶数方法应用到RSM中,它可以很好地解决RSM在低速区域的数值频散问题,以及减少不必要的时间损耗;同时讨论了旋转交错网格的理论频散特性,并基于波场分离的方法分析了实际波场传播的频散误差,将原方法的应用范围由声波推广到剪切波,由理论值推广到时变值。在数值模拟试验中,VRSM将被应用于水平层模型和Overthrust模型。通过阶数分配以及相应波场传播效果和计算时间的分析,验证了该方法应用于复杂介质波场模拟中的实用性和有效性。实验的结果表明VRSM能够合理分配不同速度所对应的差分阶数,能保证计算的精确性,并合理控制计算的时间。     

起伏地表QR径向基函数有限差分及其在弹性波逆时偏移中的应用

弹性波逆时偏移不受倾角和偏移孔径的限制,能够实现任意复杂构造的高精度多波成像,是目前最精确的多分量资料偏移成像方法之一.逆时偏移算法的核心是波场延拓,传统波场延拓以水平基准面为边界条件,基于固定采样步长进行规则网格剖分,采用阶梯近似法处理起伏地表和复杂构造界面时会产生台阶散射,严重影响起伏地表复杂构造的成像精度.基于无网格节点模型,定量分析了弹性波模拟中径向基函数有限差分法的频散关系和稳定性条件.基于此,提出一种基于QR径向基函数的高精度有限差分方法,并提出一种优化的起伏地表自适应节点剖分方法,推导了精确的无网格自由边界条件和弹性波无网格混合吸收边界条件,形成了新的基于无网格的起伏地表弹性波数值模拟方法.此外,本文将此无网格径向基函数有限差分方法应用于精确的纵横波场矢量分解公式,实现了起伏地表弹性波逆时偏移成像.通过对高斯山丘模型,起伏凹陷模型和起伏地表Marmousi-2模型进行数值试算,验证了本文方法的有效性和可行性.     

求解二阶解耦弹性波方程的低秩分解法和低秩有限差分法

时间域的波场延拓方法在本质上都可以归结为对一个空间-波数域算子的近似.本文基于一阶波数-空间混合域象征,提出一种新的方法求解解耦的二阶位移弹性波方程.该方法采用交错网格,连续使用两次一阶前向和后向拟微分算子,推导得到了解耦的二阶位移弹性波方程的波场延拓算子.由于该混合域象征在伪谱算子的基础上增加了一个依赖于速度模型的补偿项,可以补偿由于采用二阶中心差分计算时间微分项带来的误差,有效地减少模拟结果的数值频散,提高模拟精度.然而,在非均匀介质中,直接计算该二阶的波场延拓算子,每一个时间步上需要做N次快速傅里叶逆变换,其中N是总的网格点数.为了减少计算量,提出了交错网格低秩分解方法;针对常规有限差分数值频散问题,本文将交错网格低秩方法与有限差分法结合,提出了交错网格低秩有限差分法.数值结果表明,交错网格低秩方法和交错网格低秩有限差分法具有较高的精度,对于复杂介质的地震波数值模拟和偏移成像具有重要的价值.     

三维VTI介质qP波方程频率空间域有限差分高阶加权算子

波动方程有限差分法是波场模拟的一个重要方法,为解决常规有限差分法存在着数值频散的问题,本文从具有垂直对称轴的三维横向各向同性(VTI)介质频率-空间域qP波动方程出发,在常规差分算子的基础上构造了适合三维VTI介质的频率空间域有限差分优化算子,然后利用最优化理论中的Gauss-Newton法求解了优化算子的系数,使差分方程的相速度与波动方程的相速度尽量吻合,从而在理论上使网格数值频散达到极小,精度对比分析及数值测试表明,有限差分优化算子具有较高的波场数值模拟精度,有效地压制了数值频散现象,为三维VTI介质频率一空间域qP波正演模拟研究提供了理论基础.     

求解弹性波动方程的频率域近似解析离散化波场模拟方法

为提高频率域弹性波动方程数值求解的计算效率,本文引入近似解析离散化(NAD)方法将其进行数值离散并得到大型线性代数方程组.在详细分析了相应系数矩阵的稀疏分块结构与数学性质之后,本文提出采用不精确旋转分块三角预处理子加速Krylov子空间迭代方法来快速求解该线性方程组,并利用数值试验证实这种方法在弹性波场模拟方面的数值效率.通过与另外两种经典数值方法(常规有限差分方法和交错网格有限差分方法)对多种介质模型进行波场模拟、数值频散分析以及与解析解的波形对比,NAD方法显示了其在压制数值频散和提高计算效率方面的优势以及对复杂介质模型弹性波场数值模拟的有效性.     

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

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

复杂介质多尺度变网格地震波数值模拟

在应用有限差分法地震波场数值模拟过程中,模拟精度和计算效率是关键问题之一,特别是对于速度变化剧烈的含小尺度地质体的数值模拟尤为重要.为了既能精细刻画介质的局部结构,又能保证模拟的时效性,本文在传统交错网格有限差分算法的基础上,将变网格有限差分方法引入到交错网格高阶差分数值模拟中,对交错网格空间算法进行了改进,并避免了因插值因素影响模拟精度和计算效率.本文采用完全匹配层吸收边界条件实现对可变网格数值模拟的边界处理,根据完全匹配层吸收边界原理和参量的分裂思想,在不同的方向分别加相应的阻尼衰减因子,解决边界的反射问题.此外还分析了数值频散对于模拟结果的影响.文中分别设计了地堑和生物礁两个介质模型对可变交错网格算法的有效性和稳定性进行了验证.数值模拟结果表明,可变交错网格有限差分算法数值模拟能够对地下介质物性的空间变化进行准确刻画,进一步增强数值模拟对复杂介质的适应性,同时为地震数据的波场成像、纵横波联合解释等提供可靠依据.     

基于时空双变网格的起伏地表变坐标系正演模拟方法（英文）

坐标变换法通过将物理空间的曲网格映射为计算空间的矩形网格,将起伏地表转化为水平地表,同时将物理空间的波动方程转化为计算空间的波动方程,在计算空间完成数值模拟,坐标变换的方法对处理起伏自由边界具有较好的适应性和应用效果。本文在传统坐标变换方法的基础上,根据计算区域速度差异采用不同的网格大小和采样时间步长,提出了一种基于时空双变网格的起伏地表坐标变换正演模拟方法。在编程实现算法的基础上,通过典型模型波场模拟试算结果分析可知:(1)变网格方法与常规方法波场模拟误差在0.5%左右;(2)变网格方法计算效率视不同的变网格区域面积及变网格大小可提高几倍量级,在本文模型和计算参数下提高约5倍。(3)在满足模拟精度及频散条件要求下,变网格方法较全局细网格算法能显著节约计算内存。为此,针对起伏地表数值模拟,本文方法具有较高的模拟计算精度和一定的适应性。     

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

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

VTI介质qP波方程高精度有限差分算子

波动方程有限差分法是一种使用广泛的地震波数值模拟方法.但是有限差分法本身固有存在着数值频散问题，会降低地震波场模拟的精度与分辨率.为了克服常规有限差分算子的数值频散，本文针对VTI介质地震波数值模拟问题，构造了频率-空间域qP波波动方程高精度有限差分优化算子，根据最优化理论中高斯-牛顿法确定了高精度有限差分算子的优化系数.利用常规差分算子和高精度优化差分算子对归一化相速度的频散关系精度进行了对比分析，并对均匀各向同性介质和均匀VTI介质中的qP波地震波场进行了有限差分数值模拟，通过频散关系精度分析和波场数值模拟结果表明：有限差分优化算子具有较高的波场数值模拟精度，有效压制了传统有限差分算子数值模拟中的数值频散现象，提高了有限差分算子精度，为VTI介质频率-空间域qP波正演模拟奠定了基础.