简介： 本发明公开了一种新的球面坐标系下的弹性波正演模拟方法，具体涉及石油勘探领域。本发明通过引入一种新的球面坐标系及球面网格剖分策略，将笛卡尔坐标系下的弹性波方程映射到球面坐标系下，通过波场延拓，应用自由地表边界条件及最佳匹配层边界条件对边界反射进行吸收，再将球面坐标系下的波场快照反变换到笛卡尔坐标系下，输出笛卡尔坐标系下的炮记录和波场快照。本方法克服了均匀网格剖分对地下介质非均质性的影响，既保证了低速层区域满足采样定理，又避免了深部区域的过采样，不仅大幅度地提高了计算效率，还精确地还原了地震波的球面传播特征，为后续复杂地质构造成像、反演工作提供了地震波场计算基础。