一种基于Kirchhoff积分的快速三维非均匀弹性介质地震波场模拟方法

提出了一种基于Kirchhoff积分的快速地震波场数值模拟算法。在地下弹性介质的速度、反射界面的空间位置及反射率均已知的情况下,可以利用该算法计算三维非均匀弹性介质模型中的反射地震波场。讨论了算法的计算效率和计算精度。针对基于Kirchhoff积分的三维地震波场模拟计算量巨大这一问题,采用了两种尺度网格单元和限制积分区域等手段,不仅降低了计算量,节省了内存,同时还保持了足够的计算精度。利用数值模拟对算法的有效性进行了验证。

非均匀弹性介质地震波场伪谱法正演模拟

基于一阶速度—应力弹性波波动方程,利用傅里叶变换将空间导数变换到频率域求解,对介质波场进行正演模拟。本文结合实际勘探中可能存在的地质问题,设计了多种地质模型,包括背斜模型、向斜模型、透镜体模型和断层--侵入体模型等,进行了伪谱法地震波场正演模拟、合成地震记录的实现,并分析了其波场动力学变化特征。结果表明,应用伪谱法对复杂模型的模拟有较高的精度和可靠性,可直接用于辅助地震野外工作设计、反演解释和资料解释等工作。

Inhomogeneity problem with a sliding interface under remote shearing stress

The problem of an ellipsoidal inhomogeneity embedded in an infinitely extended elastic medium with sliding interfaces is investigated. An exact solution is presented for such an inhomogeneous system that is subject to remote uniform shearing stress. Both the elastic inclusion and matrix are considered isotropic with a separate elastic modulus. Based on Lur’e’s approach to solving ellipsoidal cavity problems through Lamé functions, several harmonic functions are introduced for Papkovich-Neuber displacement potentials. The displacement fields inside and outside the ellipsoidal inclusion are obtained explicitly, and the stress field in the whole domain is consequently determined.