TTI介质声波波动方程高阶有限差分正演.pdfVIP

中国地球科学联合学术年会 2014 ·901 · TTI 介质声波波动方程高阶有限差分正演 徐世刚 1,2 刘 洋 1,2 1.中国石油大学（北京）油气资源与探测国家重点实验室 北京 102249 ； 2. 中国石油大学（北京）CNPC 物探重点实验室 北京 102249 各向异性介质地震波传播的研究始终是地震学理论与应用研究的前沿和热点。为了获得弹性介质中 纵波波动方程，可以人为规定横波速度为零（Alkhalifah ，1999），即可得到声波方程。本文在前人对 VTI 介质声波方程研究的基础上（程玖兵等，2013 ），推导出 TTI 介质声波波动方程；对其采用 2N 阶 精度有限差分进行求解，用以压制差分频散，提高模拟精度；将混合吸收边界条件（Liu 和 Sen, 2010） 应用到常规网格有限差分中。通过多个模型的正演模拟，验证了 TTI 介质声波波动方程以及混合吸收 边界的可靠性和实用性。 1.TTI 介质声波波动方程推导 前人对于 VTI 介质波动方程的研究，得到 VTI 介质中的声波波动方程[1]如下： 2 2 2 2 2 ∂ u ∂ ∂ ∂ ∂ x y （1 + 2 ε） v 2 p z ( + )u + (1 + 2 δ)v 2 pz ( + )u ∂t 2 ∂x 2 ∂y 2 xy ∂z 2 ∂y 2 z 2 2 2 ∂ u ∂ u ∂ u （1） z (1 + 2 δ)v 2 pz xy + v 2 pz z ∂t 2 ∂x 2 ∂z 2 δ ε 其中v 代表 qP波垂向传播速度， 和 是 Thomsen 参数。对于 TTI 介质，对称轴方向不再是垂直方向。 pz 但从物理上讲，TTI 介质和 VTI 介质没有本质差异，只不过观测坐标系与描述介质对称性的物理坐标系 不一致而已。通过坐标旋转，可由 VTI介质波动方程推导出 TTI 介质中的对应形式。设对称轴方位角为 ϕ θ ，倾角为 ，则物理坐标X (x , y , z ) 与实际观测坐标 X＝（x,y,z）系满足如下变换关系：X=M 其中： ⎡cos ϕ − sin ϕ 0⎤ ⎡cos θ 0 − sin θ⎤ ⎡cos θ cos ϕ − sin ϕ − sin θ cos ϕ⎤ M ⎢