共查询到20条相似文献，搜索用时 62 毫秒
1.
弹性波正演模拟中改进的非分裂式PML实现方法(英文) 总被引：6，自引：3，他引：3
在弹性波有限差分正演模拟中,吸收边界条件常用来吸收截断边界处引入的不期望边界反射,其中完全匹配层(PML)吸收边界条件被认为是目前最理想的吸收边界条件。但是PML吸收边界条件的传统实现却存在着很大不足:全局分裂式PML吸收边界条件实现简单但是需要占用太多内存;局部分裂式PML吸收边界条件需要考虑多个边界和角点区域,编程实现非常复杂;非分裂式PML吸收边界条件由于涉及卷积运算,计算量很大。本文基于非分裂式PML吸收边界条件,结合复频移伸展函数,提出了一种新的数值实现方法,其计算方程简单、占用内存小、编程实现容易,是对PML介质理论数值实现的改进和完善。 相似文献
2.
《地震工程与工程振动》2016,(2)
在波动方程有限差分波场数值模拟中,为了使计算得到的边界波场值更接近于真实的边界波场值,需要使用边界条件以减少来自计算区域边界的人为反射能量。传统的分裂式完全匹配层(SPML)吸收边界条件不能有效吸收掠射波,后来发展的不分裂卷积完全匹配层(CPML)能够较好地处理这个问题,并且CPML在处理边界问题时无需对波场进行非物理分裂。本文针对Kelvin-Voigt模型的黏弹性介质,采用高阶交错网格有限差分方法进行数值模拟,同时利用CPML吸收边界条件进行边界处理。数值模拟结果表明:与弹性波相比,大角度入射产生的低频掠射波对黏弹性波的影响更大;CPML吸收边界条件对低频掠射波的吸收效果比传统的SPML吸收边界条件更好。 相似文献
3.
Liu和Sen（2010和2012）在地震波场数值模拟中提出一种混合吸收边界条件. 该方法具有计算量小、容易实现及吸收效果好等优点. 但现有的混合吸收边界条件是针对二阶位移-应力方程设计的,存在稳定性问题. 本文首先推导了两种速度-应力单程波方程：二阶Higdon单程波方程和一阶Higdon单程波方程. 进而提出基于一阶弹性波方程的混合吸收边界条件方法. 在内部区域和边界之间引入一个过渡区域,通过单程波与双程波方程平滑过渡来消除人工边界反射. 为了改善混合吸收边界条件的吸收效果和稳定性,我们采用了能同时吸收纵、横波反射的一阶单程波方程和与变量位置有关的加权系数. 为了验证混合吸收边界条件的有效性,将其与常规分裂完全匹配层（PML）方法进行了比较. 数值模拟结果表明,与PML边界条件相比,混合吸收边界条件在耗用更小计算时间和存储量的前提下,可以获得更好的吸收效果. 另外,本文提出的两种混合吸收边界条件中,混合一阶Higdon吸收边界条件具有更好的稳定性. 相似文献
4.
《应用地球物理》2019,(4)
无网格有限差分法能有效提高数值模拟的几何灵活性,且无需网格映射或复杂的网格生成过程。RBF-FD (基于径向基函数的有限差分)是最常用的无网格有限差分法之一,可以准确模拟地震波在非矩形计算域中的传播。本文提出适于弹性波方程无网格有限差分数值解的PML (完全匹配层)吸收边界条件,可以应用于非矩形速度模型的边界。但是PML吸收边界对近掠射波、低频波的吸收效果不好。为此,我们继续提出适于弹性波方程无网格有限差分数值解的CFS-PML(复频移完全匹配层)吸收边界条件。本文所提两种边界条件均是通过构造辅助微分方程,得到不分裂时域表达式,具有存储量小、便于编程实现的特点。模拟结果表明,两种PML吸收边界条件都能有效地消除无网格有限差分数值模拟的人工边界反射。此外,本文所提CFS-PML相较PML对近掠射波和损耗波的吸收效果更好。 相似文献
5.
常规完全匹配吸收边界(PML)对以近掠射角入射到界面上的波以及低频波、损耗波都会产生虚假边界反射.基于递归积分的不分裂复频移PML算法,利用复频移拉伸函数,极大地改善了PML边界条件的性能,我们进一步推导出基于递归积分的不分裂复频移PML弹性波方程交错网格高阶差分法,对长条形介质模型进行数值模拟,与常规PML算法进行比较说明该算法对以掠射角入射到PML界面的波以及PML层内损耗波的吸收效果. 相似文献
6.
《地球物理学进展》2016,(6)
在地震波传播数值模拟的过程中,需要使用吸收边界条件达到衰减人为边界反射的目的,其中完全匹配层(PML)吸收边界条件被认为是目前最理想的吸收边界条件.但在三维地震波传播数值模拟中直接应用常规的PML吸收边界条件存在编程实现复杂、计算资源需求高等不足之处.本文基于传统的NPML吸收边界条件,提出了一种改进的方法.该方法通过引入伸展函数和辅助变量,把PML介质中的三维弹性波动方程分解成正常计算项和衰减项,其中衰减项的计算只在PML区域内部进行,其值只与坐标轴的方向相关,可通过离散数值迭代求解.分析和试验结果表明,该方法避免了卷积运算,减少了辅助变量的个数,具有理论方程简洁、编程实现简单、计算资源需求少等优点,能更有效地解决三维地震波传播数值模拟中的吸收边界问题. 相似文献
7.
在地震波场数值模拟过程中,边界反射是影响其模拟结果的一个重要因素。实际地下介质具有各向异性特征,传统的完全匹配层边界（PML）对于小入射角地震波具有良好效果,但该方法并不能有效地吸收低频波和大角度入射波。针对VTI介质边界反射的问题,本文提出在频率-空间域有限差分法数值模拟中采用卷积完全匹配层（CPML）和特征分析法的组合边界条件,并对该组合边界条件进行数值模拟实验和边界反射吸收效果分析,验证所提方法是一种可靠的人工吸收边界条件,能够有效地压制波场模拟过程中产生的边界反射。 相似文献
8.
完全匹配层吸收边界在地震波模拟中已广泛使用,但常用的场分裂格式完全匹配层吸收边界(SPML)和传统的不分裂完全匹配层吸收边界(NPML)对极低频入射波或大角度入射波的边界吸收效果不好.一种无需分裂和显式卷积计算的完全匹配层吸收边界(CPML)不仅能够解决常规PML吸收边界的不足,而且具有存储量小、计算效率高、易于编程实现的特点.本文将这种完全匹配层(CPML)吸收边界引入到孔隙弹性介质速度-应力格式的旋转交错网格有限差分算法中,对完全匹配层吸收边界参数进行数值分析,得到一组优化的参数.孔隙弹性介质数值模拟结果表明这种不分裂卷积完全匹配层的吸收效果优于常规完全匹配层. 相似文献
9.
完全匹配层(PML)作为一种稳定高效的吸收边界条件,广泛应用于基于一阶电磁波动方程的探地雷达(GPR)数值模拟中.为解决基于二阶电磁波动方程的GPR数值模拟的吸收边界问题,本文借鉴二阶弹性波动方程的PML边界条件构建思想,提出了一种适合二阶电磁波动方程GPR时域有限元模拟的PML边界条件.从二阶电磁波动方程出发,基于复拉伸坐标变换,推导了PML算法的频域表达式;通过合理构造辅助微分方程,得到了PML算法的时域表达式,并以变分形式(弱形式)加载到GPR时域有限元方程中,实现了PML边界条件在二阶电磁波动方程GPR时域有限元模拟中的应用.在此基础上,对比了无边界条件、Sarma边界条件和PML边界条件下均匀模型的波场快照、单道波形、时域反射误差和能量衰减曲线,结果表明:PML边界条件的吸收效果要远优于Sarma边界条件,具有近似零反射系数.一个复杂介质模型的正演模拟验证了PML边界条件在非均匀地电结构中电磁波传播模拟的良好吸收效果. 相似文献
10.
在进行地震波模拟计算的过程中用有限的计算区域模拟地下无限空间,需要进行边界截断.为了在边界处不产生虚假反射影响模拟结果,需要引入吸收边界条件.本文采用的近似完全匹配层是一种新型非分裂完全匹配层,计算效率较高.同时相比于其他非分裂完全匹配层,其还具有不改变方程的形式、易于实现等优势.但是当入射波角度较大,边界吸收效果变弱,且残留在边界中的能量使近似完全匹配层变得极其不稳定.多轴复频移近似完全匹配层的提出就是为了改善对大角度入射波的吸收并且提高边界的稳定性.通过实验模拟和矩阵特征值灵敏度来研究多轴复频移近似完全匹配层的吸收效果及稳定性.结果表明该方法不仅能够吸收掠入波,而且对常规入射波的吸收也得到提升,同时拥有更好的稳定性. 相似文献
11.
SH波场中完全匹配层吸收边界研究(英文) 总被引：4，自引：1，他引：3
无限域中的波动方程数值模拟往往需要稳定有效的吸收边界来消除人为边界截断所引起的虚假反射.本文首先写出了全匹配层(PML)内SH波的波动方程推导结果,并给出了方程的Crank-Nicolson计算格式与其中窄间导数2阶,6阶,10阶精度的有限差分算法以及伪谱法算法.然后设计了均匀各向同性介质模型和分层溶洞模型并引入图像处理中的信噪比(SNR)概念来定量研究边界吸收效果同PML宽度、不同精度算法的关系.数值结果表明当匹配层宽度比较薄时可以用低精度的有限差格式来获得比较好的吸收效果,当匹配层比较宽时,采用高精度的算法可以获得很好的吸收效果.最后对"反射系数"进行了讨论,指出"反射系数"的不足和文中用SNR来定量衡量吸收边界 相似文献
12.
《地震地质》2020,(3)
完全匹配层吸收边界条件通常可以很好地吸收模型边界的地震反射波,但对于横向各向同性介质的模拟效果欠佳,且尚在发展阶段。为此,文中推导了横向各向同性介质中弹性动力学波动方程,给出了施加完全匹配层(PML)吸收边界条件下时间域二阶、空间域十阶精度的高阶交错网格的有限差分形式,并分别建立了均匀的垂直向对称轴的横向各向同性介质(VTI介质)和倾斜向对称轴的横向各向同性介质(TTI介质)模型。计算结果表明,对于对称轴为任意角度的横向各向同性介质,当PML边界层厚度达到一定的数值时,可以很好地抑制人工边界所产生的地震波反射效应,且PML的吸收效果不会被入射角与入射波频率影响。 相似文献
13.
《中国科学:地球科学》2021,(8)
完全匹配层(PML)是处理波动方程数值模拟中模型边界反射问题的常用方法,在时间域有限差分(FD)波场数值模拟中得到了广泛的应用.早期的PML技术存在对于近切入射和低频反射吸收效果较差的问题,由此产生了复频移PML(CFS-PML)模型边界反射压制法,其中基于递归卷积技术的复频移PML(CPML),由于其算法在时间域的高效性,被广泛应用于时间域一阶方程的数值模拟.一阶波动方程CPML边界条件是非常成熟的方法,但是二阶方程, CPML方法还需要进一步发展研究.目前,二阶方程的CPML处理方法需要引入辅助方程或者辅助变量,导致计算量增加,从而对数值模拟效率产生一定的影响.本文提出了一种分区域两步法CPML实施策略,简称为TS-CPML. TS-CPML通过对二阶CPML控制方程进行两步有限差分进行数值模拟,从而避免计算引入辅助变量的新二阶偏微分方程或其他辅助方程,具有计算效率高、编程易于实现等特点.我们将该方法应用于二阶声波方程和TTI介质二阶拟声波方程组,数值模拟显示,在相同参数下,由模型边界产生的反射波振幅约为传统二阶方程CPML方法的50%以上,验证了TSCMPL技术的优越性. 相似文献
14.
边界反射处理是地震数值模拟中需要研究的重要问题之一。本文将数学科学计算中提出的双吸收边界条件(DABC)(Hagstrom et al.,2014)应用到地震波场数值模拟中。该方法将一种局部高精度吸收边界(ABC)应用于两条平行的人工边界上,从而实现对边界反射的双吸收。本文以二维声波方程为例,设计一个基于DABC的有限差分(FD)正演方案,给出详细推导过程和实现步骤,相对于完全匹配层(PML),其理论分析和实现难度大大降低,稳定性和灵活性亦得到一定程度的增强。文章最后对典型的均质模型和SEG盐丘模型进行有限差分数值模拟实验,实验结果表明,DABC能够得到远远优于Clayton-Engquist边界条件(CEBC)的吸收效果,并与PML吸收效果相当。 相似文献
15.
复频移完全匹配层(Complex Frequency-Shifted PML,CFS-PML)在长时间时域计算中对凋落波、倏失波具有好的吸收效果,并被广泛应用于时域有限差分模拟中.而本文采用卷积方法将CFS-PML应用于时域有限元求解GPR波动方程的数值模拟中.论文以TM波为例,推导了基于CPML(Convolutional PML)边界的时域有限元GPR波动方程求解公式,采用Newmark-β方法对时间导数进行离散,有效改善了时域有限元GPR数值计算程序的稳定性.并以狭长模型为例,开展了CPML边界中关键参数m、R和κ的选取实验,通过对比反射误差大小确定了综合最优参数组合.相同时刻UPML与CPML波场快照、3个检测点的反射误差比较,说明CPML较UPML具有更好的吸收效果.最后,采用非规则四边形网格对1个复杂GPR模型进行剖分,应用加载CPML边界条件的FETD程序对该模型进行了正演,得到了二维剖面法、宽角法正演GPR剖面图,说明非规则四边形对复杂模型的良好适应性,基于CPML边界条件的FETD可有效减少边界反射误差,能实现对任意复杂不规则模型的正演模拟. 相似文献
16.
《地球物理学进展》2015,(4)
由于对边界反射具有优秀的吸收效果,完全匹配层(Perfectly Matched Layer,PML)吸收边界条件自被提出就受到了广泛的关注和研究,并发展成为地震正演数值模拟中应用最广泛的边界条件.随着地震正演数值模拟技术的发展,对PML边界条件的研究取得了显著的进展,发展形成了多种PML边界条件,并在声波、弹性波等多种方程的地震正演数值模拟中得到了广泛的应用,取得了良好的效果.但是,对于该领域涌现出的大量的研究成果,缺少系统总结的综述性文献.为此,本文归纳梳理了近年来PML边界条件在地震正演数值模拟领域的研究和应用成果,分阶段地概述了PML边界条件的发展进程.最后讨论总结了各种方法的异同之处,并对地震正演数值模拟领域中PML边界条件的研究方向进行了展望. 相似文献
17.
数值频散和边界反射是频率域模拟时需要解决的两个重要问题.然而,受计算效率和分解阻抗矩阵时的内存占用量的制约,提高有限差分算子长度或增加有限差分网格数目均不是提高频率域模拟精度的最优解决方案.本文首先分析了数值频散产生的理论机制,在此基础上,推导了一种“波数补偿”的声波方程表达式来压制数值频散,并给出其物理意义,有效地改善了数值频散问题,提高了模拟精度;在边界问题上,本文采用多轴卷积完全匹配层(MCPML)边界条件代替传统的完全匹配层(PML)边界条件,快速吸收边界内的残余能量,压制边界反射.结合改进声波方程和MCPML边界条件,给出了一种高精度的频率域声波方程有限差分格式.数值模拟结果表明,在不增加计算量和内存占用量的前提下,本文研究的方法、正演精度高、波场模拟清晰、无干扰反射,是一种可靠高效的频率域模拟方法. 相似文献
18.
多分量联合逆时偏移最佳匹配层吸收边界 总被引：3，自引：2，他引：1
有限空间内的波动方程逆时偏移需要利用有效的边界处理技术用以消除人工截断对偏移结果产生的影响。本文以横向各向同性介质弹性波速度-应力方程为基础,依据传统分裂式最佳匹配层（Perfect Matched Layer,PML）吸收边界技术的思想,推导了应用于逆时偏移的完全匹配层波动方程,并给出了其高阶交错网格有限差分格式。针对由边界处向计算区域内传播的＂反射波＂,以及地震记录排列两端地震同相轴突变对计算区域的影响这两方面问题,本文给出了逆时偏移中吸收层的布设方式。模型和实际资料的弹性波叠前多分量联合逆时深度偏移结果表明本文的边界处理方法取得了较好的吸收效果,获得了好的联合偏移成像结果。 相似文献
19.
在地震波传播数值模拟的过程中,需要使用吸收边界条件从而达到衰减人为边界 反射的目的. 本文基于傍轴近似法提出了计算三维弹性波方程的吸收边界条件公式,表示了 各边界面、边棱和角点处波场所满足的单程波方程,并在三维弹性波数值模拟中进行了应用 . 理论模型及三维盐丘地质模型波场切片快照试算结果表明,该吸收边界条件可以有效地吸 收人为边界反射,适用于较大入射角情况,从而消除了边界有效波信息的干扰. 由于采用四 阶近似方程,在保证计算精度的前提下,该方法具有节省计算工作量和易于实现的特点. 相似文献