弹性波波场P波和S波分解的数值模拟.docVIP

弹性波波场 P 波和 S波分解的数值模拟马德堂Ξ( 长安大学基础课部)朱光明( 长安大学应用地球物理研究所)摘要马德堂, 朱光明.弹性波波 弹性波波场 P 波和 S 波分解的数值模拟 马德堂Ξ ( 长安大学基础课部) 朱光明 ( 长安大学应用地球物理研究所) 摘 要 马德堂, 朱光明. 弹性波波场 P 波和 S 波分解的数值模拟. 石油地球物理勘探, 2003, 38 (5) : 482～ 486 通常采用完全弹性波波动方程进行弹性波波场数值模拟, 只能得到 P 波和 S 波的混合波场。若要进一步获 得 P 波和 S 波的波场, 一般是在得到混合波场之后再进行波场分离, 或用声波方程单独模拟 P 波和 S 波波场, 但是这些方法往往很难确保 P 波和 S 波的振幅不出现畸变。本文给出了完全弹性波波动方程的一种等价方程, 该方程既含混合波场变量, 也含纯 P 波和纯 S 波波场变量。用伪谱法求解该波动方程, 在得到混合波场的同时, 也得到了完全分离的纯 P 波或纯 S 波波场, 并保留了 P 波和 S 波能量相互转换的信息。 模拟波场分析表明, 采 用这种弹性波波场分解的数值模拟方法对认识弹性波的传播规律有重要意义。 关键词 弹性波 波动方程 波场分解 伪谱法 用于 P 波和 S 波分解的 弹性波波动方程 引 言 通常进行弹性波波场数值模拟采用完全弹性波 波动方程, 只能得到 P 波和 S 波的混合波场。 在地 震资料常规处理、偏移成像以及地震反演中, 通常将 波场的垂直分量看成 P 波, 水平分量看成 S 波, 这 种近似只有在表层含低速带时, 才是合理的。从严格 意义上讲, 按 P 波或 S 波处理地震资料时, 必须先 对地震资料做 P 波和 S 波分离。一般作法是把梯度 和旋度算子作用于混合波场, 即可实现 P 波和 S 波 分离的目的, 但这样做却很难确保 P 波和 S 波的振 幅及其之间的能量相互转换信息不出现畸变。 在弹 性波动力学理论中, 通常用标量位的梯度和向量位 的旋度之和表达弹性波波场的 P 波和 S 波分解, 以 此来进行数值模拟是不太方便的。 本文通过直接引入 P 波波场变量和 S 波波场 变量, 给出可实现 P 波和 S 波分解的等价方程; 用 伪谱法求解该波动方程, 同时得到了混合波场以及 完全分离的纯 P 波和纯 S 波波场, 并完全保留了 P 波和 S 波能量的相互转换信息。波场模拟表明, 弹性 波波场分解的数值模拟对认识弹性波的传播规律有 重要意义。 在均匀各向同性弹性介质中, 用位移表示的三 维弹性波波动方程为 2 Θ5 S = Λ) A Η+ ΛA 2S (Κ+ 5t2 或 Θ5 S = 2 2Λ) A Η+ Λ(A 2S - A Η) (Κ+ 5t2 5u i+ 5v j + 5w k; S = (u , v , w ) 为位移场; A 其中: Η= 5x 5y 5z 为梯度算子。 上式还可写成 52S V P A Η+ 2 V S (A S - A Η) 2 2 5t2 = 其中 Κ+ 2Λ V P = Θ Λ V S = Θ V P , V S 分别为纵、横波速度。 用 u , v , w 各位移分量表示的弹性波波动方 程为 Ξ 本文于 2002 年 11 月 19 日收到。 52 u52 u52 v52w52 u52 u52 v52w225t2 =5x 2 +5x 5y ++5y 2 +5z 2 -5x 5y -V PV S5x 5z52w5x 5z52w52 v52 u52 v52 v52 v52 u225t2 =5x 5y +5y 2 ++5x 2 +5z 2 -5x 5y -(1)V PV S5y 5z52w5y 5z52w52 u52 v52w52w52 u52 v22= V P5x 5z +5y 5z ++V S5x 2 + 5y 2 -5x 5z - 5y 5z5t252 z为了表达弹性波场的 P 波和 S 波分解, 现引进新变量A ×S P 关于时间的一阶偏导数为常数, 故A ×S P 关于么为常数, 要么为线性函数。 由波动特性知, A ×S P 关于时间只能为常 52 u 52 u 52 v 52w 52 u 52 u 52 v 52w 2 2 5t2 = 5x 2 + 5x 5y + + 5y 2 + 5z 2 - 5x 5y - V P V S 5x 5z 52w 5x 5z 52w 52 v 52 u 52 v 52 v 52 v 52 u 2 2 5t2 = 5x 5y + 5y 2 + + 5x 2 + 5z 2 - 5x