各向异性双相介质弹性波场褶积算法数值模拟

从各向同性介质中波场数值模拟的褶积微分算子法出发,推导出了各向异性双相介质中波场传播数值计算的褶积新算法.将常见的二阶微分Biot波动方程用等效的一阶速度—应力双曲方程表示,其中未知的波场向量包括固相和流体的速度分量和应力分量,由此对方程的时间项使用交错网格差分方法计算,而对空间项则采用褶积微分算法进行求解.对各向异性双相介质在单层介质模型和双层介质模型中的波场特征进行了研究.研究的结果显示,在两层介质分界面上当地震波产生反射时能观测到两类纵波和横波,并且在衰减系数大的介质里慢纵波很难见到.     

2.5维地震波场褶积微分算子法数值模拟

早期的褶积微分算子都是基于正反傅立叶变换而实现的,其精度比四阶有限差分的精度稍高,本文将计算数学中的Forsyte广义正交多项式微分算子与褶积算子相结合,构建了一个新的快速、高精度褶积微分算子,其计算结果非常接近实验函数微分的精确值,精度与16阶有限差分的精度相当,远优于错格伪谱法的精确度.另外,2.5维数值模拟比二维模拟可以更真实地模拟三维介质的臬个剖面的波场,并且2.5维地震波模拟的计算量比三维模拟的计算量及计算耗时要大大减少.本文利用基于Forsyte广义正交多项式褶积微分算子法计算2.5维非均匀介质地震波场,模拟结果表明,该算法的计算速度快,计算精度高,能够直观、高效地反映复杂介质中波场的传播规律,并且2.5维波场数值模拟具有更高的计算效率,是一种非常值得深入研究并广泛应用的方法.     

基于辛格式离散奇异褶积微分算子的弹性波场模拟

本文发展了基于辛格式离散奇异褶积微分算子(SDSCD)的保结构方法模拟弹性波场,求解弹性波动方程时,引入辛差分格式进行时间离散,采用离散奇异褶积微分算子进行空间离散.相比于传统的伪谱方法,该方法提高了计算精度和稳定性.数值结果表明SDSCD方法可以有效地抑制数值频散,为解决大尺度、长时程地震波场模拟问题提供了合适的数值方法.     

粘弹性介质地震波传播的褶积微分算子法数值模拟研究

早期的褶积微分算子法都是基于正反傅立叶变换而实现的,其精度比四阶有限差分稍高。本文将计算数学中的Forsyte广义正交多项式微分算子与褶积算子相结合,构建了一个新的快速、高精度褶积微分算子,其计算结果非常接近实验函数微分的精确值,精度与l6阶有限差分相当。粘弹性波动方程更真实地描述了实际地下介质中弹性波的传播规律及其波场特征。本文以二维粘弹性波动方程为例,推导了粘弹性介质波动方程的离散格式,用迭积微分算子法实现了粘弹性介质的地震波场正演模拟,并对其波传播特征进行了分析。计算结果表明该算法能正确模拟粘弹性介质中的地震波,正确地反映粘弹性介质中波场的传播规律。     

地震波数值模拟的褶积微分算子法与伪谱法的对比分析

将基于Forsyte广义正交多项式的褶积微分算子法运用于复杂非均匀介质地震波场模拟中,并将计算结果与伪谱法计算结果进行分析比较。通过二者的计算时间对比发现:在同样的计算条件下,褶积微分算子法的采样时间始终小于伪谱法,这是其进行地震波数值模拟的一个明显优势。通过波场快照的对比,褶积微分算子法的模拟结果与伪谱法数值模拟结果的频散效应相当,可为地震波场的值计算提供一种新的选择。     

各向异性介质地震波场的优化褶积微分算子法数值模拟

在前人工作基础上,通过对窗函数参数进行优化实现了对基于Shannon奇异核理论的交错网格褶积微分算子的优化过程.应用这种优化褶积微分算子方法对各向异性介质进行了数值模拟,讨论了优化褶积微分算子法模拟的PML吸收边界条件以及稳定性条件,分析了弹性波在此类介质中的传播特征,并与高阶交错网格有限差分方法进行了对比.数值实验结果表明,该方法适用于各向异性介质中弹性波场模拟,精度高,稳定性好,是一种研究复杂介质中地震波传播的有效数值方法.     

基于辛格式奇异核褶积微分算子的地震标量波场模拟

本文在对地震波场进行模拟时,采用辛差分格式对波动方程进行时间离散,采用奇异核褶积微分算子对波动方程进行空间离散.该方法尽管增加了一些计算量,但提高了计算精度和稳定性;相对于其他非辛算法,它是全局保结构的,并且具有较强的长时间跟踪能力.该方法为解决大尺度、长时程地震波场的高精度模拟问题提供了一种新的、有效的选择.     

最优化辛格式广义离散奇异核褶积微分算子地震波场模拟

将波动方程变换至Hamilton体系,构造了一种新的保结构算法,即最优化辛格式广义褶积微分算子(OSGCD). 在时间离散上,首先引入了Lie算子设计二级二阶辛格式,基于最小误差原理得到了优化的辛格式. 在空间离散上,引入广义离散奇异核褶积微分算子计算空间微分,提出了一种有效方法优化GCD并得到了稳定的算子系数. 针对本文发展的新方法,给出了OSGCD稳定性条件. 在数值实验中,将OSGCD与多种方法比较,从精度和计算效率两方面分析了OSGCD的计算优势,计算结果也表明OSGCD长时程以及非均匀介质中地震波模拟亦具有较强能力.     

地震波数值模拟中的最优Shannon 奇异核褶积微分算子

如何有效协调地震波数值模拟过程中的计算精度和计算效率,一直是许多地球物理工作者密切关注且努力寻求解决的问题.这一问题求解的好坏,将直接影响到人们正确认识地下构造的地球物理响应.本文从离散Shannon奇异核理论出发,推导了基于Shannon奇异核的交错网格褶积微分算子,分析并求出了影响算子精度窗函数的最优化参数.通过将...     

基于Shannon奇异核理论的褶积微分算子在地震波场模拟中的应用

本文在前人工作的基础上,建立了一种基于Shannon奇异核的交错网格褶积微分算子方法.文中不仅详细讨论了影响算子精度的各种因素,同时也着重分析了其在弹性波模拟中的频散关系和稳定性条件.通过和交错网格有限差分算子比较,发现该算子即使在高波数域也具有较高的精度.均匀介质中的数值试验也表明,该方法9点格式就基本上达到了解析解精度.而分层均匀介质和复杂介质中的地震波数值模拟也同时证实了该方法精度高,稳定性好,是一种研究复杂介质中地震波传播的有效数值方法.     

复杂介质地震波传播的褶积微分算子数值模拟h

将佛尔塞（Forsyte）广义正交多项式微分算子地震波场正演模拟算法运用于复杂非均匀介质模型的波场数值模拟中,并比较了该方法在计算效率和计算精度方面与有限差分方法和伪谱法的差异. 数值结果表明,这种广义正交多项式微分算子法计算速度快、精度高,对计算资源需求低,是一种颇具潜力的数值模拟方法.      

双相介质地震波场数值模拟的迭积微分算子及其PML边界条件

将基于计算数学中Forsyte 广义正交多项式的迭积微分算子引入到地震波动方程的一阶速度--应力方程的空间微分运算中去,并采用时间错格有限差分算子替代传统的差分算子以匹配高精度的空间迭积微分算子,从而发展一种全新的地震波场正演模拟方法,来解决复杂非均匀介质模型中的波场传播问题.为了大幅衰减人工边界引起的反射,本文将完全匹配层(Perfectly Matched Layer,PML)吸收边界条件引入到所构建的方法中,以解决迭积微分算子法的边界问题.以二维波动方程为例,用迭积微分算子法实现了双相介质的地震波场正演模拟,模拟结果表明,双相介质模型较好地解释了含流体孔隙特性.同时也表明迭积微分算子法是一种非常实用、有效的数值模拟方法.     

三维弹性波数值模拟中的吸收边界条件

在地震波传播数值模拟的过程中,需要使用吸收边界条件从而达到衰减人为边界 反射的目的. 本文基于傍轴近似法提出了计算三维弹性波方程的吸收边界条件公式,表示了 各边界面、边棱和角点处波场所满足的单程波方程,并在三维弹性波数值模拟中进行了应用 . 理论模型及三维盐丘地质模型波场切片快照试算结果表明,该吸收边界条件可以有效地吸 收人为边界反射,适用于较大入射角情况,从而消除了边界有效波信息的干扰. 由于采用四 阶近似方程,在保证计算精度的前提下,该方法具有节省计算工作量和易于实现的特点.     

错格傅里叶伪谱微分算子在波场模拟中的应用

本文在传统傅里叶微分矩阵的基础上,对原始微分算子进行改进,引入了错格微分算子的傅里叶伪谱方法．尽管该方法增加了一些计算量,但却极大地提高了计算精度和稳定性．而且,该方法将微分计算过程由传统的傅里叶变换转换为一般的矩阵矢量乘积,大大降低了微分求解过程的复杂程度．在均匀介质中,将错格伪谱微分算子计算的结果和解析解进行比较,结果表明本文算子几乎达到了解析解的精度．而在分层均匀介质中的实验结果同时显示,该方法精度高、稳定性好,是一种研究层状介质中地震波传播的有效数值方法．     

基于广义正交多项式褶积微分算子的地震波场数值模拟方法

地震波场模拟方法研究对于与波动现象有关的地震学问题的重要性是不言而喻的.就目前现有的各种正演算法来说,精度较高的算法（如有限元法、谱元法、高阶有限差分法等）,其计算速度较慢;计算速度较快的算法（如低阶有限差分法、付氏伪谱法等）计算精度却比较低.为了兼顾地震波场模拟的精度与速度,本文推出了一种快速的、高精度地震波场模拟方法（基于Forsyte广义正交多项式的褶积微分算子法）,该方法是以计算数学中的Forsyte广义正交多项式插值函数为基础,构建一个新的褶积微分算子,并将该算子引入到地震波动方程的一阶速度-应力方程的空间微分运算中去,采用时间交错网格有限差分算子替代普通的差分算子以匹配高精度的褶积微分算子,从而构造一种全新的地震波场数值模拟方法.该方法同时具有广义正交多项式方法的高精度和短算子低阶有限差分算法的高速度.通过对算子长度的调节及算子系数的优化,可同时兼顾波场解的全局信息与局部信息.复杂非均匀介质模型中的波场数值模拟实验证实了该方法的可行性及优越性.     

交错网格高阶差分法三维弹性波数值模拟

本文应用交错网格高阶有限差分方法模拟弹性波在三维各向同性介质中的传播。采用时间上二阶、空间上高阶近似的交错网格高阶差分公式求解三维弹性波位移-应力方程,并在计算边界处应用基于傍轴近似法得到的三维弹性波方程吸收边界条件。在此基础上进行了三维盐丘地质模型的地震波传播数值模拟试算。试算结果表明该方法模拟精度高,在很大程度上减小了数值频散,绕射波更加丰富,而且适用于介质速度具有纵向变化和横向变化的情况。     

各向异性弹性波场的有限元数值模拟

研究了各向异性弹性波有限元正演系统的精度和效率问题，提出了一种透射加衰减的组合人工边界方案（吸收边办条件），它对各向异性波具有较好的吸收效果，并且有较好的稳定性能，均匀TI介质中的模拟获得了非常清晰的波场快照，其波场特征与理论分析能够准确吻合，各向异性介质模型的地表地震记录表明，各向异性波炮集记录在波的类型、同相轴形态、能量分布和相位等方面与各向同性波都有很大差别。     

三角网格有限元法声波与弹性波模拟频散分析

本文对声波与弹性波方程进行有限元法离散,构造有限元法频散关系的一般特征值问题,分析了时间离散格式为中心差分的三角网格有限元法声波与弹性波模拟的频散特性. 比较了三种质量矩阵即分布式质量矩阵、集中质量矩阵和混合质量矩阵对有限元法频散的影响;选取四种典型三角网格,分析了混合质量矩阵有限元（MFEM）频散的方向各向异性;数值频散、方向各向异性随插值阶数的增加逐渐减弱,当空间为三阶插值时,频散主要表现为随采样率的变化而几乎无明显方向各向异性, 其频散幅值也较小. 控制其他影响因素不变的情况下,研究了不同波速比介质中弹性波的数值频散. 最后给出了三角网格MFEM的数值耗散性.     

伪谱法弹性波场数值模拟中的边界条件

边界条件问题长期以来一直是困扰地震波数值模拟研究者的一个难题,许多人都提出了有效的方法,但是该问题仍然需要研究.本文针对伪谱法地震波场数值模拟的特殊要求,即所有网格点相互耦合,且傅立叶变换的周期性使得方程在网格边界上很难得到一个吸收的替代方程,研究衰减边界条件中衰减系数的变化对边界反射波衰减效果的影响.数值实验结果说明,随着衰减系数的增大,振幅衰减加快,用很少的过渡带网格就可以使边界上的反射波能量很小.但是如果衰减率过大,就会在传播区和过渡带产生干扰波场.因此,在衰减系数的选择上应采取折中办法,即在取某一较小的衰减系数的前提下,设置较少的过渡带网格数使边界的能量反射最少.