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

Lebedev grid finite-difference modeling for irregular free-surface and stability analysis based on M-PML boundary condition

作　　者：黄建平 杨宇 李振春 田坤[2] 李庆洋 

机构地区：[1]中国石油大学地球科学与技术学院,山东青岛266580 [2]中石化胜利油田物探研究院,山东东营257022

出　　处：《中国石油大学学报（自然科学版）》2016年第4期47-56,共10页Journal of China University of Petroleum(Edition of Natural Science)

基　　金：国家'973'课题(2014CB239006;2011CB202402);国家自然科学基金项目(41274124);山东省自然科学基金项目(ZR2011DQ016);中央高校基本科研业务费专项(R1401005A)

摘　　要：贴体网格作为一种'边界拟合网格',能改善传统有限差分处理起伏地表时由阶梯状网格存在所产生的虚假绕射。此外,完全匹配层(PML)技术被证明是一种能够有效消除虚假反射的重要边界条件,但采用自由地表模拟面波时,传统PML易产生不稳定现象。为实现起伏地表自由边界条件下稳定的地震波场模拟,发展一种基于弹性介质的全交错-Lebedev网格有限差分方法,并将边界条件扩展为多轴完全匹配层(M-PML)。在实现算法的基础上,通过模型试算验证新方法的正确性;并对比几种不同完全匹配层技术的稳定性差异,解释PML边界产生不稳定的机制。正演模拟结果表明,M-PML在处理起伏地表面波模拟中有更好的稳定性。Body-fitted grid is usually ter^ned as ＂ boundary conforning grid ＂ , which can eliminate artifacts caused by staircaseapproximation of irregular free surface in classical finite-difference method. In addition, perfectly m atched layer （ PM L）absorbing boundary is approved to be an efficient method to suppress spurious edge reflections. When modeling Rayleighwaves with the existence of free-surface,however, the classical PML algorithm may become unstable. In order to obtain stableseismic wave simulation with an irregular free-surface, this paper proposes a full staggered grid （FSG ）-Lebed evgrid finite-difference method in elastic media, and extends the boundary condition to multiaxial PML （ M -PM L ）. Numerical simulationsdemonstrate the validity of the proposed method. The stability of the tw-o perfectly matched layer techniques is comparedand the instability mechanism of PML is discussed. The num erical results suggest preferred stability of the M-PML techniqueover the classical PML algorithm.

关 键 词：地震波场模拟 起伏地表 Lebedev网格 多轴完全匹配层 稳定性 

分 类 号：P631.4[天文地球—地质矿产勘探]