有限差分正演中PML 边界条件优化算法研究

来源 :2016中国地球科学联合学术年会 | 被引量 : 0次 | 上传用户:lvjieidd
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
有限差分法是进行波动方程正演模拟的重要手段,而边界条件是其关键问题之一,目前应用最广,效果最好的边界条件是完全匹配层(PML)边界条件.然而传统的PML边界条件存在一定缺陷,因而本文研究了PML边界条件的优化算法.本文给出了各向同性介质中波动方程的有限差分格式,并从边界问题入手,介绍了完全匹配层(PML)边界条件的基本思想,对传统的PML及其优化算法进行了详细的推导,进而实现数值模拟,通过对比分析表明,卷积PML和Z变换PML的迭代格式可以有效减小计算量和内存占有量,复频移函数及其高阶形式的引入改善了吸收效果,多轴PML和SMART Layers的提出消除了不稳定性.
其他文献
地震各向异性被证明广泛存在于地下介质中.随着地震勘探技术的发展,深入研究各向介质中的地震波场十分必要,有着重要的理论意义和应用价值.各向异性弹性波控制方程参数多、不
会议
全波形反演理论与技术以其高精度、多参数建模的能力成为当今地球物理勘探领域的研究热点,但在数学上它又是一个病态的非线性问题,其多解性、收敛性以及庞大的计算代价在一定
面波是一种强而广泛存在的规则干扰波,在炮集上呈近似线性分布,能量强、频率低、且振动延续时间长,严重影响中深层有效反射.常规F-K方法虽然能够较好地压制面波,但由于F-K滤
会议
弹性波叠前逆时偏移需要处理多分量数据,而基于射线方法的多分量Kirchhoff偏移在介质情况太复杂的时候并不适用.基于波动方程的逆时偏移方法面临两个难题:P波和S波分离;推导
会议
近年来,TI介质qP波正演模拟受到广泛关注.基于常规耦合型传播方程的正演模拟方法存在伪横波及频散假象干扰.使用Lowrank近似进行qP波一步法波场外推.使用一个解析的源来构建
提高分辨率是反褶积的主要目的之一,它是通过消除大地滤波作用的方式来实现的,处理效果直接影响后续的叠加、偏移以及地质解释,尤其是当前地震勘探工作的不断深入,勘探区域更
在现今的深度偏移方法中,傅里叶法作为一类重要的单程波偏移方法,具有较低数值频散的优势.傅里叶法将空间的偏导数用波数来代替,因此算法稳定且无数值频散.但是为了满足高精
会议
全波形反演方法利用叠前地震波场的运动学和动力学信息重建地层结构,具有揭示复杂地质背景下构造与储层物性的潜力.全波形反演理论和技术以其高精度、多参数建模的能力而备受
会议
盐丘内幕速度的准确性直接影响了盐体及盐下成像质量.全波形反演技术(FWI)利用观测数据和正演数据的差异最小准则,迭代修正初始模型速度,最终可以提供高精度的速度体.本文把
角度道集是公认的没有假象的道集,因此角度道集可以有效的应用于速度分析和逆时偏移成像等.地震波在传播过程中,往往会产生转换波.而随着地震数字处理技术的不断发展,研究的