完全匹配层吸收边界条件的研究

边界反射是数值模拟中的一个关键问题。构建完全匹配层(PML)吸收边界条件是消除边界反射的理想方法之一。此方法是在研究区域的边界上加入吸收层,使边界上传入吸收层的波随传播距离按指数规律衰减,不产生任何反射,达到消除边界反射的目的。首先给出了二维声波方程完全匹配层的公式,然后通过地面地震模型和井间地震模型数值算例与常规吸收边界条件进行了比较。结果表明,PML边界条件具有吸收更干净、相位对应关系好且能够吸收各种角度的边界反射等优点,从而验证了完全匹配层吸收效果的优越性。     

TI介质地震波场数值模拟边界条件处理

地震波场数值模拟研究中,边界条件处理非常重要,它直接影响模拟地震数据的精度.针对TI(Transversely Isotropic)介质一阶速度—应力弹性波动方程组,利用时间域分裂法给出了PML(Perfectly matched layer,完全匹配层)吸收边界条件,其基本思想是在所研究的区域外引入吸收层,地震波从有效计算区域传播到吸收层时不发生任何反射,且在吸收层内按传播距离呈指数衰减,不发生反射.利用交错网格高阶差分技术对PML吸收边界条件求解,对比不同边界条件数值模拟结果表明,PML吸收边界条件较应用广泛的Clayton吸收边界有更好的吸收效果,能更好地满足地震数值模拟的需要.     

声学FDTD法完全匹配层数和衰减系数对衰减性能的影响

分析了声学时域有限差分(FDTD)法中完全匹配层(PML)层数、衰减系数对吸收效果的影响.通过建立研究PML边界吸收性能的模型发现:PML边界层层数和衰减系数的选择相互关联;边界层数越多,PML的吸收效果越好,但计算量增大;层数一定时,需要为各衰减层设置合适的衰减系数,以减少衰减系数差异带来的数值反射;在实际计算中,应...     

起伏地表有限元叠前逆时偏移完全匹配层

在地震波传播的数值模拟中,用有限区域代替求解无限区域地震波传播规律是一种有效手段。目前完全匹配层吸收边界条件是边界吸收最好的选择。通常正演模拟不考虑自由表面吸收条件,但在逆时偏移中,地表反射的存在会引起假的同相轴。推导了位移形式的有限元二阶声波方程的完全匹配层吸收边界条件公式,并应用于起伏地表声波方程有限元数值模拟与逆时偏移成像;通过数值分析并与传统的Clayton-Engquist吸收边界条件做了比较。计算结果表明,方法取得了理想的边界吸收效果,可以消除起伏地表反射引起的假同相轴。     

波动方程正演模拟边界条件的比较分析

通过数值模拟研究了透明边界、Clayton-Engquist边界和完全匹配层边界的吸收效果,得出如下结论:在反射角和频率相同的情况下,完全匹配层边界条件效果最好,ClaytonEngquist边界效果次之,而透明边界条件的效果最差。以边界条件对100 Hz模型边界垂直反射的吸收效果来衡量,完全匹配层边界条件与Clayton-Engquist边界条件的效果分别是透明边界条件的16.5倍和3.5倍。在40 Hz的低频范围内,或者在反射角65°的情况下,Clayton-Engquist边界相对透明边界的吸收效果相对优势显著变弱。而完全匹配层边界的吸收效果则在150Hz频率范围内和75°反射角范围内始终保持稳定的相对优势。     

CFS-PML在探地雷达FDTD法模拟中应用

探地雷达(GPR)数值模拟中,为了使截断边界处不产生人为的反射波,需要加载边界条件.目前常用的完全匹配层(PML)边界对隐失波、低频波、掠角波等不能很好地吸收.为此,在时域有限差分(FDTD)法中引入复频移卷积完全匹配层(CFS-PML).从二维TM波的Maxwell方程组出发,阐述了CFS-PML边界条件的原理,推导了CFS-PML边界加载的离散差分公式,给出了CFS-PML边界条件的参数选取原则.应用均匀介质中GPR波场快照,对比了截断处没有加载边界条件、加载3种不同边界条件的边界处反射波的强弱,实例说明单轴各向异性完全匹配层(UPML)、CFS-PML边界条件具有较好的吸收效果.为了进一步对比UPML与CFS-PML的吸收效果,设置了双层介质狭长模型,通过波场快照与全局反射误差证明了CFS-PML对隐失波、低频波、掠角波的吸收较UPML效果更佳.     

基于M-PML边界的Lebedev网格起伏地表正演模拟方法及稳定性分析

贴体网格作为一种"边界拟合网格",能改善传统有限差分处理起伏地表时由阶梯状网格存在所产生的虚假绕射。此外,完全匹配层(PML)技术被证明是一种能够有效消除虚假反射的重要边界条件,但采用自由地表模拟面波时,传统PML易产生不稳定现象。为实现起伏地表自由边界条件下稳定的地震波场模拟,发展一种基于弹性介质的全交错-Lebedev网格有限差分方法,并将边界条件扩展为多轴完全匹配层(M-PML)。在实现算法的基础上,通过模型试算验证新方法的正确性;并对比几种不同完全匹配层技术的稳定性差异,解释PML边界产生不稳定的机制。正演模拟结果表明,M-PML在处理起伏地表面波模拟中有更好的稳定性。     

基于复频移完全匹配层的辛龙格库塔算法探地雷达正演

为了消除探地雷达(ground penetrating radar,GPR)正演模拟的截断边界处的非物理反射,需要对模拟区域边界进行特殊处理.目前应用效果较好的边界处理方法为加载完全匹配层吸收边界条件.而复频移完全匹配层(complex frequency shift perfectly matched layer,C...     

瑞雷波数值模拟中的边界条件及模拟实例分析

应用2×12阶高精度交错网格有限差分法,建立震源位于地表时瑞雷波数值模拟的自由边界条件:同时,采用PML吸收边界条件以消除模拟时的边界反射.通过对均匀各向同性介质及2层介质模型的计算,验证所采用的自由边界条件和吸收边界条件的有效性及波场模拟的正确性.研究结果表明:模拟结果与解析结果相对误差在5%以内,并指出在同一精度下通过提高差分阶数比减少网格间距所需计算时间少.在2层介质模拟的基础上,对含有软弱夹层的3层介质进行模拟,得到更加接近实际情况的地震记录,频散曲线产生多个模式,并出现"之"字型回折.     

起伏地表复杂地质模型的格子法声波方程正演模拟

利用声波方程格子法,结合非规则网格完全匹配层(PML)吸收边界,对具有起伏地表、强速度对比的模型进行了正演模拟.模拟结果表明,采用依据波速改变尺度的非则网格对模型进行离散化,不仅能够精确刻画地表、反射界面,避免阶梯近似带来的数值散射现象,而且能有效克服规则网格离散化带来的空间、时间过采样问题,节省了存储空间,提高了模拟效率.非规则网格PML方法以及相应的坐标转换方法的应用,可对所关心的区域周围设置曲线吸收边界,避免了常规PML方法中处理拐角时的麻烦.