共查询到10条相似文献，搜索用时 957 毫秒
1.
《地球物理学进展》2016,(6)
在地震波传播数值模拟的过程中,需要使用吸收边界条件达到衰减人为边界反射的目的,其中完全匹配层(PML)吸收边界条件被认为是目前最理想的吸收边界条件.但在三维地震波传播数值模拟中直接应用常规的PML吸收边界条件存在编程实现复杂、计算资源需求高等不足之处.本文基于传统的NPML吸收边界条件,提出了一种改进的方法.该方法通过引入伸展函数和辅助变量,把PML介质中的三维弹性波动方程分解成正常计算项和衰减项,其中衰减项的计算只在PML区域内部进行,其值只与坐标轴的方向相关,可通过离散数值迭代求解.分析和试验结果表明,该方法避免了卷积运算,减少了辅助变量的个数,具有理论方程简洁、编程实现简单、计算资源需求少等优点,能更有效地解决三维地震波传播数值模拟中的吸收边界问题. 相似文献
2.
Liu和Sen（2010和2012）在地震波场数值模拟中提出一种混合吸收边界条件. 该方法具有计算量小、容易实现及吸收效果好等优点. 但现有的混合吸收边界条件是针对二阶位移-应力方程设计的,存在稳定性问题. 本文首先推导了两种速度-应力单程波方程：二阶Higdon单程波方程和一阶Higdon单程波方程. 进而提出基于一阶弹性波方程的混合吸收边界条件方法. 在内部区域和边界之间引入一个过渡区域,通过单程波与双程波方程平滑过渡来消除人工边界反射. 为了改善混合吸收边界条件的吸收效果和稳定性,我们采用了能同时吸收纵、横波反射的一阶单程波方程和与变量位置有关的加权系数. 为了验证混合吸收边界条件的有效性,将其与常规分裂完全匹配层（PML）方法进行了比较. 数值模拟结果表明,与PML边界条件相比,混合吸收边界条件在耗用更小计算时间和存储量的前提下,可以获得更好的吸收效果. 另外,本文提出的两种混合吸收边界条件中,混合一阶Higdon吸收边界条件具有更好的稳定性. 相似文献
3.
弹性波正演模拟中改进的非分裂式PML实现方法(英文) 总被引：6，自引：3，他引：3
在弹性波有限差分正演模拟中,吸收边界条件常用来吸收截断边界处引入的不期望边界反射,其中完全匹配层(PML)吸收边界条件被认为是目前最理想的吸收边界条件。但是PML吸收边界条件的传统实现却存在着很大不足:全局分裂式PML吸收边界条件实现简单但是需要占用太多内存;局部分裂式PML吸收边界条件需要考虑多个边界和角点区域,编程实现非常复杂;非分裂式PML吸收边界条件由于涉及卷积运算,计算量很大。本文基于非分裂式PML吸收边界条件,结合复频移伸展函数,提出了一种新的数值实现方法,其计算方程简单、占用内存小、编程实现容易,是对PML介质理论数值实现的改进和完善。 相似文献
4.
本文将大地电磁场分解为一次场和二次场,应用交错网格有限差分法模拟计算大地电磁二次场,并引入各向异性最佳匹配层(PML)吸收边界条件作为二次场边界条件,实现了耦合PML吸收边界条件的三维大地电磁二次场有限差分正演模拟.为了确保正演的稳定性和效率,QMR求解器和磁感应矢量散度校正技术被用于PML吸收边界条件下系数矩阵的快速求解.三维模型正演响应表明,基于二次场的三维大地电磁有限差分算法具有较高的计算精度和可靠性.通过计算分析不同PML吸收因子条件的大地电磁正演结果,显示在适当的吸收因子下,PML吸收边界条件可较大幅度的减小外边界距离,从而有效的压缩模型求解空间,最终提高三维大地电磁正演模拟的效率. 相似文献
5.
《应用地球物理》2017,(2)
地震波场数值模拟中不可避免地会出现边界反射,一般采用吸收边界条件以压制人工边界反射。目前常用的分裂式完全匹配层(PML)边界条件需要在边界处进行特殊处理,尤其是在三维情况下需要将变量分裂为三个分量,增加了数值模拟的计算时间和内存占用量。与分裂式PML吸收边界条件相比,混合吸收边界条件(HABC)具有易于实现、计算量小和吸收效果好等优点,可以提高三维波动方程数值模拟的计算效率。本文将基于一阶Higdon单程波方程的混合吸收边界条件从二维计算域发展到三维,提出了适用于三维弹性波数值模拟的混合吸收边界条件。均匀模型以及复杂模型的三维数值模拟结果表明,混合吸收边界条件与传统的完全匹配层边界条件相比,具有效率高、吸收效果好的优势。 相似文献
6.
完全匹配层(PML)作为一种稳定高效的吸收边界条件,广泛应用于基于一阶电磁波动方程的探地雷达(GPR)数值模拟中.为解决基于二阶电磁波动方程的GPR数值模拟的吸收边界问题,本文借鉴二阶弹性波动方程的PML边界条件构建思想,提出了一种适合二阶电磁波动方程GPR时域有限元模拟的PML边界条件.从二阶电磁波动方程出发,基于复拉伸坐标变换,推导了PML算法的频域表达式;通过合理构造辅助微分方程,得到了PML算法的时域表达式,并以变分形式(弱形式)加载到GPR时域有限元方程中,实现了PML边界条件在二阶电磁波动方程GPR时域有限元模拟中的应用.在此基础上,对比了无边界条件、Sarma边界条件和PML边界条件下均匀模型的波场快照、单道波形、时域反射误差和能量衰减曲线,结果表明:PML边界条件的吸收效果要远优于Sarma边界条件,具有近似零反射系数.一个复杂介质模型的正演模拟验证了PML边界条件在非均匀地电结构中电磁波传播模拟的良好吸收效果. 相似文献
7.
《中国科学:地球科学》2021,(8)
完全匹配层(PML)是处理波动方程数值模拟中模型边界反射问题的常用方法,在时间域有限差分(FD)波场数值模拟中得到了广泛的应用.早期的PML技术存在对于近切入射和低频反射吸收效果较差的问题,由此产生了复频移PML(CFS-PML)模型边界反射压制法,其中基于递归卷积技术的复频移PML(CPML),由于其算法在时间域的高效性,被广泛应用于时间域一阶方程的数值模拟.一阶波动方程CPML边界条件是非常成熟的方法,但是二阶方程, CPML方法还需要进一步发展研究.目前,二阶方程的CPML处理方法需要引入辅助方程或者辅助变量,导致计算量增加,从而对数值模拟效率产生一定的影响.本文提出了一种分区域两步法CPML实施策略,简称为TS-CPML. TS-CPML通过对二阶CPML控制方程进行两步有限差分进行数值模拟,从而避免计算引入辅助变量的新二阶偏微分方程或其他辅助方程,具有计算效率高、编程易于实现等特点.我们将该方法应用于二阶声波方程和TTI介质二阶拟声波方程组,数值模拟显示,在相同参数下,由模型边界产生的反射波振幅约为传统二阶方程CPML方法的50%以上,验证了TSCMPL技术的优越性. 相似文献
8.
数值频散和边界反射是频率域模拟时需要解决的两个重要问题.然而,受计算效率和分解阻抗矩阵时的内存占用量的制约,提高有限差分算子长度或增加有限差分网格数目均不是提高频率域模拟精度的最优解决方案.本文首先分析了数值频散产生的理论机制,在此基础上,推导了一种“波数补偿”的声波方程表达式来压制数值频散,并给出其物理意义,有效地改善了数值频散问题,提高了模拟精度;在边界问题上,本文采用多轴卷积完全匹配层(MCPML)边界条件代替传统的完全匹配层(PML)边界条件,快速吸收边界内的残余能量,压制边界反射.结合改进声波方程和MCPML边界条件,给出了一种高精度的频率域声波方程有限差分格式.数值模拟结果表明,在不增加计算量和内存占用量的前提下,本文研究的方法、正演精度高、波场模拟清晰、无干扰反射,是一种可靠高效的频率域模拟方法. 相似文献
9.
全局反射误差分析是深入研究探地雷达(GPR)吸收边界条件吸收效率的有力工具.基于常规完全匹配层(PML)的标准交错网格有限差分算法必须满足严格的CFL条件限制,即在单位时间步长内,不容许电磁波传播的距离超过单元网格尺寸.为了提高主区域所有网格节点的计算效率,并有效地吸收传播后期出现的低频隐失波,提出基于非分裂递归卷积完全匹配层(UCPML)的旋转交错网格(RSG)GPR正演算法.该算法采用不同的网格交错策略,并在边界条件中引入了吸收低频隐失波的自由可变因子,使得该算法允许选取较大的时间步长,提高了计算效率,并且实现了对低频隐失波的高效吸收.本文首先给出了RSG差分格式,推导了基于UCPML的RSG差分更新方程,探讨了数值色散的稳定性条件,然后以绕射现象严重的衰减夹层数值模拟为例,分别从波场快照、单道波记录、时间域反射误差(TDRE)、频率域反射误差(FDRE)四个方面分析了UCPML与常规PML的全局反射误差,说明了基于UCPML和RSG的GPR正演算法能更有效地吸收低频隐失波. 相似文献
10.
《地震地质》2020,(3)
完全匹配层吸收边界条件通常可以很好地吸收模型边界的地震反射波,但对于横向各向同性介质的模拟效果欠佳,且尚在发展阶段。为此,文中推导了横向各向同性介质中弹性动力学波动方程,给出了施加完全匹配层(PML)吸收边界条件下时间域二阶、空间域十阶精度的高阶交错网格的有限差分形式,并分别建立了均匀的垂直向对称轴的横向各向同性介质(VTI介质)和倾斜向对称轴的横向各向同性介质(TTI介质)模型。计算结果表明,对于对称轴为任意角度的横向各向同性介质,当PML边界层厚度达到一定的数值时,可以很好地抑制人工边界所产生的地震波反射效应,且PML的吸收效果不会被入射角与入射波频率影响。 相似文献