论文部分内容阅读
在地震勘探中,习惯的做法是通过弹性波场分解技术把P波场和S波场分离出来以达到成像的目的。可是这个目标有时很难达到,特别是当速度场较为复杂的时候。有一种有效的多分量分析和模拟方法是直接求解弹性波动方程而得到纯P波场和纯S波场,得到的方程称为分离的弹性波动方程。在本文中,我们首先采用高阶交错网格有限差分法对比了两种这样的波动方程:一阶(速度-应力)和二阶(位移-应力)分离的弹性波动方程,与一阶(速度-应力)和二阶(位移-应力)完全(混合)的弹性波动方程。然后比较了得到的波场快照、共炮点道集的地震剖面和单独的合成地震图。模拟结果表明得到的结果是等价的,不管是一阶或二阶分离的弹性波方程都能获得纯P波场和纯S波场。将纯P波场和纯S波场叠加得到的混合波场与用对应的一阶或二阶完全弹性波动方程计算出的混合波场结果相同。求解这种混合方程的计算代价稍微低于分离方程。分离的方程主要用于实现分离的P和S波场,且可以用来测试多分量波场分离程序的有效性。二阶分离的弹性波方程提供了纯P波场和纯S波场的位移信息,因此是一个很好的选择。其次,我们还使用高阶交错网格有限差分法模拟各向异性介质中的混合波场,可以得到包含准纵波和准横波的混合波场。为了从各向异性介质中分离qP波和qSV波,使用各向同性介质中的散度和旋度算子是没有效的。我们要通过求解Christoffel方程得到的伪求导算子来分离纯qP波和qSV波。虽然伪求导算子可以得到纯qP波和qSV波,但它需要大量内存的计算,且实际模型的模拟比较复杂。