较复杂介质波场计算中的一种透射边界

传统的人工边界是在声波或弹性波分解为内行和外行分量的基础上设计的.对于较复杂的介质模型（如粘弹介质、弹塑介质、塑性介质模型等）,地震波不能简单地分解为各自满足其运动方程的纵波、横波等成分,且地震波在传播过程中还有衰减现象,因此传统的人工边界已不能有效地用于非弹性的较复杂介质的波场模拟.本文提出了一种适合于更一般的地震模型的波场计算的透射边界条件,并以粘弹介质的一个较简单的二维几何模型为例,通过有限元数值求解这种透射边界和刚性边界下粘弹波的激发和传播,得到地面合成地震记录,证明了这种透射边界用于较复杂介质波场计算是可行的.     

频率-空间域CPML和特征分析法组合边界研究

在地震波场数值模拟过程中,边界反射是影响其模拟结果的一个重要因素。实际地下介质具有各向异性特征,传统的完全匹配层边界（PML）对于小入射角地震波具有良好效果,但该方法并不能有效地吸收低频波和大角度入射波。针对VTI介质边界反射的问题,本文提出在频率-空间域有限差分法数值模拟中采用卷积完全匹配层（CPML）和特征分析法的组合边界条件,并对该组合边界条件进行数值模拟实验和边界反射吸收效果分析,验证所提方法是一种可靠的人工吸收边界条件,能够有效地压制波场模拟过程中产生的边界反射。      

瑞雷波在垂向低速薄层中透射和反射的理论与实验研究

用三维超声模拟研究来自垂向低速薄层反射以及在该薄层中透射的瑞雷波。基于格林函数的近似法计算了瑞雷波的反射和透射系数。非刚性接触的边界条件用于低速层模拟,得到的理论与实验的反射和透射系数的绝对值相当吻合。实验与理论研究结果的某些不一致,特别是相移,可以用理论模型与超声模拟研究的低速层不相适应来解释。     

ABAQUS中动力问题边界条件的选取

数值模拟是研究半无限域中结构地震反应的重要手段,而边界条件的选取则直接影响到数值模拟结果的可靠性.本文应用有限元软件ABAQUS,分别采用固定边界条件、黏性边界条件和远置边界条件对简单弹性半无限介质的波源问题和散射问题进行了数值模拟分析,结果显示:使用固定边界时,波动不能穿过边界面进入半无限介质,这与波动在半无限介质中...     

地震正演数值模拟完全匹配层吸收边界条件研究综述

由于对边界反射具有优秀的吸收效果,完全匹配层(Perfectly Matched Layer,PML)吸收边界条件自被提出就受到了广泛的关注和研究,并发展成为地震正演数值模拟中应用最广泛的边界条件.随着地震正演数值模拟技术的发展,对PML边界条件的研究取得了显著的进展,发展形成了多种PML边界条件,并在声波、弹性波等多种方程的地震正演数值模拟中得到了广泛的应用,取得了良好的效果.但是,对于该领域涌现出的大量的研究成果,缺少系统总结的综述性文献.为此,本文归纳梳理了近年来PML边界条件在地震正演数值模拟领域的研究和应用成果,分阶段地概述了PML边界条件的发展进程.最后讨论总结了各种方法的异同之处,并对地震正演数值模拟领域中PML边界条件的研究方向进行了展望.     

直流电阻率2.5维有限元数值模拟

从稳定电流场满足的基本微分方程出发,推导了地面点源电场电位边值问题对应的等价变分问题.用网格离散积分区域时,以四边形为最小研究单元.针对电场在电源附近衰减快,变化梯度大,数值模拟结果在电源附近位置有较大误差问题,采用边界/奇异校正算法改善模拟结果在电源附近的奇异性问题,同时减少边界条件的影响,提高模拟结果精度.从波数域电位求取空间域电位,利用最优化离散波数进行反付氏变换.模型算例结果验证了该算法的可行性与有效性.     

耦合PML边界条件的三维大地电磁二次场有限差分法

本文将大地电磁场分解为一次场和二次场,应用交错网格有限差分法模拟计算大地电磁二次场,并引入各向异性最佳匹配层(PML)吸收边界条件作为二次场边界条件,实现了耦合PML吸收边界条件的三维大地电磁二次场有限差分正演模拟.为了确保正演的稳定性和效率,QMR求解器和磁感应矢量散度校正技术被用于PML吸收边界条件下系数矩阵的快速求解.三维模型正演响应表明,基于二次场的三维大地电磁有限差分算法具有较高的计算精度和可靠性.通过计算分析不同PML吸收因子条件的大地电磁正演结果,显示在适当的吸收因子下,PML吸收边界条件可较大幅度的减小外边界距离,从而有效的压缩模型求解空间,最终提高三维大地电磁正演模拟的效率.     

弹性波正演模拟中改进的非分裂式PML实现方法(英文)

在弹性波有限差分正演模拟中,吸收边界条件常用来吸收截断边界处引入的不期望边界反射,其中完全匹配层(PML)吸收边界条件被认为是目前最理想的吸收边界条件。但是PML吸收边界条件的传统实现却存在着很大不足:全局分裂式PML吸收边界条件实现简单但是需要占用太多内存;局部分裂式PML吸收边界条件需要考虑多个边界和角点区域,编程实现非常复杂;非分裂式PML吸收边界条件由于涉及卷积运算,计算量很大。本文基于非分裂式PML吸收边界条件,结合复频移伸展函数,提出了一种新的数值实现方法,其计算方程简单、占用内存小、编程实现容易,是对PML介质理论数值实现的改进和完善。     

基于PML边界条件的二阶电磁波动方程GPR时域有限元模拟

完全匹配层(PML)作为一种稳定高效的吸收边界条件,广泛应用于基于一阶电磁波动方程的探地雷达(GPR)数值模拟中.为解决基于二阶电磁波动方程的GPR数值模拟的吸收边界问题,本文借鉴二阶弹性波动方程的PML边界条件构建思想,提出了一种适合二阶电磁波动方程GPR时域有限元模拟的PML边界条件.从二阶电磁波动方程出发,基于复拉伸坐标变换,推导了PML算法的频域表达式;通过合理构造辅助微分方程,得到了PML算法的时域表达式,并以变分形式(弱形式)加载到GPR时域有限元方程中,实现了PML边界条件在二阶电磁波动方程GPR时域有限元模拟中的应用.在此基础上,对比了无边界条件、Sarma边界条件和PML边界条件下均匀模型的波场快照、单道波形、时域反射误差和能量衰减曲线,结果表明:PML边界条件的吸收效果要远优于Sarma边界条件,具有近似零反射系数.一个复杂介质模型的正演模拟验证了PML边界条件在非均匀地电结构中电磁波传播模拟的良好吸收效果.     

有限元模拟中边界条件对计算结果的影响

本文从总电位法和异常电位法两个方面,分别就直流电法有限元中所采用的三类边界条件对计算结果的影响进行了分析,为采用合适的边界条件及选取最佳边界距离来提高正演精度提供了基础.结果表明,混合边界条件精度最高,可大大缩小求解区域而不影响计算精度,其次是Dirichlet边界条件精度较高,但测点越靠边界误差会越大,必须取足够大的边界区域,齐次边界条件的误差比较大,但如果采用非二极装置,通过电位差计算得到的视电阻率,由于无穷远边界对电位差的影响基本消除,视电阻率计算误差与混合边界条件下的接近,在反演中,为了节省计算时间,经常使用齐次边界条件进行有限元正演.     

盾构隧道联络通道地震响应规律研究

联络通道与盾构隧道的连接处是隧道抗震的薄弱环节,但目前有关联络通道抗震性能的研究开展得还很少。对于隧道与联络通道连接部位这类典型的空间结构,三维动力有限元分析是最为有效的方法,而其中的人工边界条件和地震动输入方法是实现这一方法的关键技术之一。本文以显式形式将多次透射边界条件引入三维动力分析控制方程,采用分割算法分别求解边界节点和内部节点的响应,然后以对武汉长江隧道工程为背景,利用这种算法研究了隧道与联络通道连接处的地震响应规律,并对采用柔性搭接代替刚性接头和地基加固两种措施的减振效果进行了计算分析。     

多孔弹性波方程的多尺度波场模拟

本文研究了二维多孔弹性波方程的多尺度波场数值模拟方法.该多尺度方法可采用较粗的网格计算,同时又能反映细尺度上物性参数的变化信息.文中详细阐述了多尺度模拟方法与算法,并推导了相应的计算格式.基本思想是建立粗细两套网格,在粗网格上,基于有限体积方法计算更新波场;在细网格上,计算多尺度基函数,这基于有限元方法通过求解一个局部化问题得到.对含有随机分布散射体的多孔介质模型进行了数值计算,计算中应用了完全匹配层(PML)吸收边界条件,数值结果验证了本文方法和算法的正确性和有效性.     

多源条件下直流电阻率法有限元三维数值模拟中一种近似边界条件

在多源直流电阻率法有限元三维数值模拟中,传统混合边界条件由于刚度矩阵与源点位置相关,无法实现线性方程组的快速回代求解,目前常用齐次边界条件或无限元边界进行替代,虽然实现了快速回代求解,但同时也降低了数值模拟的精度.为了实现快速回代求解,并确保数值模拟的计算精度,本文提出了一种近似边界条件方法.其核心思想是将与场源位置相关的边界系数矩阵从刚度矩阵中分离出来,使得分离后的刚度矩阵与场源位置无关.并将边界系数矩阵与边界处一次电场向量的乘积移到线性方程组右端源项中,当场源位置发生改变时,只需要重新计算右端源项就可以实现快速回代求解.理论模型数值计算表明,在水平地形条件下,本文边界条件数值精度优于混合边界条件;在起伏地形条件下,与齐次边界条件相比,本文边界条件数值结果与混合边界条件吻合度更高.     

高精度双吸收边界条件在地震波场模拟中的应用（英文）

边界反射处理是地震数值模拟中需要研究的重要问题之一。本文将数学科学计算中提出的双吸收边界条件(DABC)(Hagstrom et al.,2014)应用到地震波场数值模拟中。该方法将一种局部高精度吸收边界(ABC)应用于两条平行的人工边界上,从而实现对边界反射的双吸收。本文以二维声波方程为例,设计一个基于DABC的有限差分(FD)正演方案,给出详细推导过程和实现步骤,相对于完全匹配层(PML),其理论分析和实现难度大大降低,稳定性和灵活性亦得到一定程度的增强。文章最后对典型的均质模型和SEG盐丘模型进行有限差分数值模拟实验,实验结果表明,DABC能够得到远远优于Clayton-Engquist边界条件(CEBC)的吸收效果,并与PML吸收效果相当。     

地震波有限差分数值模拟框架下的起伏地表处理方法综述

本文查阅了近30年来国内外关于有限差分方法模拟近地表复杂介质方面的文献,分析了在有限差分法框架下复杂地表地形起伏问题的处理及实现方法.总结了地震波有限差分数值模拟方面的成果.复杂地表类型很多,对地震资料的影响程度也不尽相同,要想完全描述近地表影响因素是不可能的,而关键的问题就是自由边界条件方程的处理.所有处理及实现方法中不外两种:直接法(局部)和间接(全局)法.文中列出的方法,各有特点.我们的目的是寻找合适的方法解决适合的问题.     

2维和2.5维起伏地表直流电法有限差分数值模拟

起伏地表直流电场数值模拟现多采用有限元法,主要是因为其有灵活的处理曲边界的能力,然而有限元法比有限差分法要复杂,如果让有限差分法也同样具有较好的处理曲边界的能力,那数值模拟将变得更为简单.本文通过在非正则内点处采用不等距差分,在起伏地表点处直接实现边界条件,克服转移法人为改变地表形状的弊端,使得基于笛卡尔网格的有限差分...     

三维弹性波数值模拟中改进的NPML研究

在地震波传播数值模拟的过程中,需要使用吸收边界条件达到衰减人为边界反射的目的,其中完全匹配层(PML)吸收边界条件被认为是目前最理想的吸收边界条件.但在三维地震波传播数值模拟中直接应用常规的PML吸收边界条件存在编程实现复杂、计算资源需求高等不足之处.本文基于传统的NPML吸收边界条件,提出了一种改进的方法.该方法通过引入伸展函数和辅助变量,把PML介质中的三维弹性波动方程分解成正常计算项和衰减项,其中衰减项的计算只在PML区域内部进行,其值只与坐标轴的方向相关,可通过离散数值迭代求解.分析和试验结果表明,该方法避免了卷积运算,减少了辅助变量的个数,具有理论方程简洁、编程实现简单、计算资源需求少等优点,能更有效地解决三维地震波传播数值模拟中的吸收边界问题.     

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

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

安新—宽城测深剖面的时域有限差分数值模拟

利用时域有限差分方法进行地震波模拟能够得到丰富的波场信息, 对研究复杂地质构造条件下的地震波波场十分有效。 该文在对时域有限差分方法进行概述的基础上, 以2002年4月施测的安新—宽城地震测深剖面资料为基础, 采用时域有限差分方法对该剖面进行二维全波数值模拟, 同时以理想匹配层作为模拟的吸收边界条件。 将模拟得到的合成地震图, 通过与野外实验得到的记录截面和射线追踪合成记录相比较和讨论, 验证了时域有限差分方法在模拟地震波传播过程的有效性。     

一阶弹性波方程数值模拟中的混合吸收边界条件

Liu和Sen（2010和2012）在地震波场数值模拟中提出一种混合吸收边界条件. 该方法具有计算量小、容易实现及吸收效果好等优点. 但现有的混合吸收边界条件是针对二阶位移-应力方程设计的,存在稳定性问题. 本文首先推导了两种速度-应力单程波方程：二阶Higdon单程波方程和一阶Higdon单程波方程. 进而提出基于一阶弹性波方程的混合吸收边界条件方法. 在内部区域和边界之间引入一个过渡区域,通过单程波与双程波方程平滑过渡来消除人工边界反射. 为了改善混合吸收边界条件的吸收效果和稳定性,我们采用了能同时吸收纵、横波反射的一阶单程波方程和与变量位置有关的加权系数. 为了验证混合吸收边界条件的有效性,将其与常规分裂完全匹配层（PML）方法进行了比较. 数值模拟结果表明,与PML边界条件相比,混合吸收边界条件在耗用更小计算时间和存储量的前提下,可以获得更好的吸收效果. 另外,本文提出的两种混合吸收边界条件中,混合一阶Higdon吸收边界条件具有更好的稳定性.