2.5维声波数值模拟

利用Liner的2.5-D时间域声波方程和高阶有限差分近似,研究了2.5维声波方程在时间-空间域中的衰减特性。简单地质模型(水平二层介质)的数值计算结果证实了下列事实:在Liner的声波方程中的最后一项,即含有1/t2的一项对整个计算结果没有影响,因此可以忽略。这一结论与文献中已经发表的数值结果相吻合。同样的结论对于由Williamson和Pratt提出的2.5-D声波方程也适用。

基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移

时间域常Q黏声波方程,由于含分数阶时间导数项,数值求解需要大量内存,计算效率低,不利于地震偏移的实施.通过一系列近似,可将该方程简化为介质频散效应和衰减效应解耦的分数阶拉普拉斯算子黏声波方程,数值求解内存需求少,计算效率高.本文采用交错网格有限差分逼近时间导数,改进的伪谱法计算空间导数,PML吸收边界去除边界反射,对该方程进行数值离散和地震正演模拟,开展地震数据的黏声介质逆时偏移,实现波场逆时延拓过程中同时完成频散校正和衰减补偿.改善深层构造的成像精度,数值结果表明,基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移可大幅度提高地震模拟计算效率,偏移剖面明显优于常规声波偏移剖面,极大改善深层构造的成像品质.