CN103149586A - 一种倾斜层状粘弹性介质中波场正演模拟方法 - Google Patents

Abstract

本发明公开了一种倾斜层状粘弹性介质中波场正演模拟方法，该方法针对地质模型为含倾斜界面的粘弹性层状介质，将点源激发的球面波分解为平面波，利用矢量波动方程不依赖于坐标系这一性质，在固定坐标系和随倾斜界面变化的动坐标系下交替地讨论波的传播。首先由传统的反射和透射系数推导出单个倾斜界面上平面波的反射和透射系数，得到该界面上的反射及透射波；然后递推求得平面波经过多个倾斜界面时的反射和透射系数。最后由传播后的各平面波合并得到点源激发波场，合并时采用了快速的振荡积分计算方法和一条适合粘弹性介质的积分路径。

Description

一种倾斜层状粘弹性介质中波场正演模拟方法

技术领域

本发明属于地震勘探技术领域，涉及一种模拟方法，尤其是一种倾斜层状粘弹性介质中波场正演模拟方法，其是一种针对含倾斜界面的层状粘弹性介质中地震波场高效精确的正演模拟技术。

背景技术

波场模拟在油气勘探中占有重要地位，尽管已存在多种类型的波场模拟方法，然而，针对要解决的问题的要求，发展波场模拟理论及方法一直是很活跃的研究领域。另外，在油气勘探等领域，为了解释地震图，或进行数据反演，常常需要仅模拟指定类型的波，譬如仅模拟一次反射波，或仅模拟转换波等。常用的波场模拟方法有多种类型，如基于射线理论的方法，数值解方法，广义反射率法以及其推广而得到的反/透射矩阵法等。这些方法各有其优缺点。

射线方法能适用于较复杂的介质模型，可根据事先给定的射线码模拟感兴趣的类型的波，且计算速度快；不足之处是，它只是波动方程的高频近似，用于模拟复杂储层的地震响应具有局限性。采用数值法求解波动方程，是复杂介质中波场模拟的有效途径。常用的数值方法有：有限差分法，有限元法，伪谱法和谱元法等，这类方法可模拟地震波全波场，精度高。但计算量巨大，不能像射线类方法，有选择性的仅模拟我们感性趣的类型的波。

层状介质中波传播理论及波场模拟技术一直受到人们的重视。Aki和Müller等人发展的反射率法，能快速精确地模拟水平层状粘弹性介质的全波场或指定类型的波。Berkhout提出了波场外推法可用于带倾斜界面的弹性介质波场模拟。Chen等人发展的全局广义反/透射矩阵法和边界元法，能模拟带不规则界面的层状弹性介质中全波场。但这些方法不适用于含倾斜界面的层状粘弹性介质。

发明内容

针对上述问题，本发明的目的在于针对含倾斜波阻抗界面的层状粘弹性介质，提出一种倾斜层状粘弹性介质中波场正演模拟方法，该方法可以选择性模拟感兴趣的各种类型的波，如一次反射P波、S波和转换波，且模拟的波场非常精确。

本发明的目的是通过以下技术方案来解决的：

该种倾斜层状粘弹性介质中波场正演模拟方法，包括以下步骤：

1）将点源激励的球面波分解为平面波；

2）将平面波在倾斜层状粘弹性介质中传播，遇到界面时进行如下处理：以该界面为

轴，界面与全局坐标系Z轴的交点为原点建立局部坐标系，垂直于界面向下，将全局坐标系中的入射波转换到局部坐标系中；在局部坐标系下求得反射波和透射波；再将局部坐标系下的反射波和透射波转换到全局坐标系下，得到全局坐标系下的反射系数和透射系数，以及反射波和透射波；在每个界面上根据需要选取波的类型，P波入射时能够选取的波包括反射PP波、反射PS波、透射PP波、透射PS波；S波入射时能够选取的波包括反射SP波、反射SS波、透射SP波、透射SS波；从而模拟所选取的波类型，得到模拟的平面波；

3）由步骤2）模拟得到的平面波根据下式合成点源激发的波场：

$P=\frac{1}{2\mathrm{\π}}{\∫}_{\mathrm{\Γ}}\underset{i=1}{\overset{N-1}{\mathrm{\Σ}}}[R_i\left(p\right)\·\frac{\exp \left[\mathrm{j\ω}(p_{1i}x-q_{1i}z)\right]}{-2\mathrm{jq}}]\mathrm{dp}$

式中，Γ为积分路径；N为平面波的个数；R_i(p)为平面波在多个界面上的透射和反射系数综合；j是虚数单位，ω是角频率；p_1i、q_1i分别为第i个界面反射回第1层的平面波的水平和垂直慢度；p、q分别为点源处平面波的水平和垂直慢度，x和z表示位置坐标。

进一步的，上述步骤2）中，模拟震源为P波时的一次反射P波，设震源和检波器均在第一层，将震源处入射波逐层向下传播的递推步骤为：

①每遇到一个界面i，即以该界面为轴，界面与全局坐标系Z轴的交点为原点建立局部坐标系

垂直于界面向下；将全局坐标系中的入射波φ_i转换到局部坐标系中得

②在局部坐标系下求得反射系数和透射系数

以及反射波

和透射波

③再将局部坐标系下的反射波和透射波转换到全局坐标系下，得到全局坐标系下的反射系数和透射系数

以及反射波φ′_i，i和透射波φ_i+1；

④反射波φ′_i，i逐层向上递推到第一层，遇到界面时执行①-③步，但只取透射波；

⑤透射波φ_i+1继续向下递推，遇到下一个界面时回到步骤①。

进一步，在上述步骤2）中，模拟其它类型的波采用与模拟震源为P波时的一次反射P波的相同方法和步骤进行。

在上述步骤3）中，式中积分方法采用Frazer和Gettrust提出的GFM积分法。

本发明具有以下几点有益效果：

（1）本发明适用于带倾斜界面的层状粘弹性介质，各界面的倾斜方向和倾斜角度可以不同，只需满足界面不相交的条件。

（2）本发明能正确模拟粘弹性介质中波的衰减特征，相对于有限差分法，具有较高的计算效率，更适合层较厚、且层间参数变化剧烈情况。

（3）本发明非常灵活，可以分别模拟各种波，如一次反射P波、SV波、SH波和PS、SP转换波，也可模拟我们感兴趣的多次波，只需将希望模拟的波类型和相应的层次以编码形式给定，入射波入射角超过临界角时会出现首波。

（4）本发明模拟的波场非常接近解析解，可为数值方法解提供一个对比参照。所以本发明在地震资料解释和反演中有重要应用价值和应用前景。

附图说明

图1是本发明针对的倾斜层状粘弹性介质模型；

图2是波入射到界面上时的反射和透射示意图；

图3是本发明中平面波合成点源激励波场时的积分路径示意图。

具体实施方式

本发明针对的地质模型为含倾斜界面的粘弹性层状介质，将点源激发的球面波分解为平面波，利用矢量波动方程不依赖于坐标系这一性质，在固定坐标系和随倾斜界面变化的动坐标系下交替地讨论波的传播。首先由传统的反射和透射系数推导出单个倾斜界面上平面波的反射和透射系数，得到该界面上的反射及透射波；然后递推求得平面波经过多个倾斜界面时的反射和透射系数。最后由传播后的各平面波合并得到点源激发波场，合并时采用了快速的振荡积分计算方法和一条适合粘弹性介质的积分路径。

下面结合附图和实施例对本发明进行详细的描述。

首先本发明的主要符号定义如表1所示：

表1本发明中主要符号约定

图1为倾斜层状粘弹性介质模型示意图，各层介质的密度为ρ_i(i＝1,2,…,N-1)、品质因子为Q_i、P波速度为α_i、S波速度为β_i，各倾斜界面的倾角为θ_i，各倾斜界面与垂直方向的Z轴交点为z_i。点源在第一次层，检波器可以在任意位置。以水平方向为X轴，垂直向下方向为Z轴，建立全局坐标系XOZ。

1）将点源激励的球面波分解为平面波

根据Weyl积分，均匀介质中点源激励的球面波表示为一系列平面波的叠加，二维情况下每个平面波的表达式为：

φ=exp[jω(px+qz)]

上式中省略了时间因子exp(-jωt)，j是虚数单位，ω是角频率，p是水平慢度，q是垂直慢度，x和z表示位置坐标。

2）将平面波在倾斜层状粘弹性介质中传播，遇到界面时进行如下处理：以该界面为

轴，界面与全局坐标系Z轴的交点为原点建立局部坐标系，

垂直于界面向下，将全局坐标系中的入射波转换到局部坐标系中；在局部坐标系下求得反射波和透射波；再将局部坐标系下的反射波和透射波转换到全局坐标系下，得到全局坐标系下的反射系数和透射系数，以及反射波和透射波。在每个界面上可根据需要选取波的类型，P波入射时能够选取的波包括反射PP波、反射PS波、透射PP波、透射PS波；S波入射时能够选取的波包括反射SP波、反射SS波、透射SP波、透射SS波；从而仅模拟我们感兴趣的波类型。例如需模拟震源为P波时的一次反射P波，设震源和检波器均在第一层，将震源处入射波逐层向下传播的递推步骤为：

①每遇到一个界面i，即以该界面为

轴，界面与全局坐标系Z轴的交点为原点建立局部坐标系垂直于界面向下。将全局坐标系中的入射波φ_i转换到局部坐标系中得

②在局部坐标系下求得反射系数

和透射系数

以及反射波

和透射波

③再将局部坐标系下的反射波和透射波转换到全局坐标系下，得到全局坐标系下的反射系数

和透射系数

以及反射波φ′_i，i和透射波φ_i+1；

④反射波φ′_i，i逐层向上递推到第一层，遇到界面时执行①-③步，但只取透射波（上行波）；

⑤透射波φ_i+1继续向下递推，遇到下一个界面时回到步骤①。

3）由平面波合成点源激发的波场；

$P=\frac{1}{2\mathrm{\π}}{\∫}_{\mathrm{\Γ}}\underset{i=1}{\overset{N-1}{\mathrm{\Σ}}}[R_i\left(p\right)\·\frac{\exp \left[\mathrm{j\ω}(p_{1i}x-q_{1i}z)\right]}{-2\mathrm{jq}}]\mathrm{dp}$

其中，R_i(p)为平面波在多个界面上的透射和反射系数综合，如模拟一次P-P反射波时 $R_i\left(p\right)=\underset{n=1}{\overset{i-1}{\mathrm{\Π}}}{T}_{\mathrm{pp},n}^u\·R_{\mathrm{pp},i}^d\·\underset{m=1}{\overset{i-1}{\mathrm{\Π}}}{T}_{\mathrm{pp},m}^d;$ p和q为震源处入射平面波的水平和垂直慢度；p_1i和q_1i为平面波经第i个界面反射后到达检波器所在层时的水平和垂直慢度，它们是入射波水平慢度p的函数。Γ为积分路径，本发明选择如图2所示的积分路径，1/α_i和1/β_i是P波和S波的慢度，路径在第一象限中的拐点B的实部应稍大于max(1/β_i)，对单个接收器，θ=tan^-1(xz_r-z_s|)时积分收敛最快（x、z_r、z_s分别是偏移距和接收器与源的z轴坐标），如是多个接收器同时接收，可取某个折中值。线段OB的倾角应使得OB段非常靠近第一象限内最小极点的下方。积分方法采用Frazer和Gettrust提出GFM(Generalization ofFilon’s Method)积分法。

步骤2）中入射波在界面上产生的反射波和透射波的具体计算方法为：

如图2，设有一水平慢度为p₁的P波φ₁从原点出发在介质I中向下传播，振幅为1，则该平面波可写为：

φ₁=exp[jω(p₁x+q₁z)],

将入射P波由全局坐标系转换到局部坐标系，得到局部坐标系下的入射P波为：

${\widehat{\mathrm{\φ}}}_1=A_1\exp \left[\mathrm{j\ω}({\hat{p}}_1\hat{x}+{\hat{q}}_1\hat{z})\right],$

其中：

A₁=exp[jω(p₁x₀+q₁z₀)], $\left[\begin{array}{c}{\hat{p}}_1\\ {\hat{q}}_1\end{array}\right]=C\left[\begin{array}{c}{p}_1\\ q1\end{array}\right],$ $C=\left[\begin{array}{cc}\cos \mathrm{\θ}& \sin \mathrm{\θ}\\ -\sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right].$

在局部坐标系

下，P波

从介质I向下入射到界面，则在界面处会产生反射回介质I的上行P波

和SV波

以及透射到介质II中的下行P波

和SV波

各反射波和透射波的表达式分别为：

${\widehat{\mathrm{\φ}}}_1^{\′}={\hat{R}}_{\mathrm{pp}}^d\·A_1\exp \left[\mathrm{j\ω}({\hat{p}}_1\hat{x}-{\hat{q}}_{1p}\hat{z})\right],$

${\widehat{\mathrm{\ψ}}}_1^{\′}={\hat{R}}_{\mathrm{ps}}^d\·A_1\exp \left[\mathrm{j\ω}({\hat{p}}_1\hat{x}-{\hat{q}}_{1s}\hat{z})\right],$

${\widehat{\mathrm{\φ}}}_2={\hat{T}}_{\mathrm{pp}}^d\·A_1\exp \left[\mathrm{j\ω}({\hat{p}}_1\hat{x}+{\hat{q}}_{2p}\hat{z})\right],$

${\widehat{\mathrm{\ψ}}}_2={\hat{T}}_{\mathrm{ps}}^d\·A_1\exp \left[\mathrm{j\ω}({\hat{p}}_1\hat{x}+{\hat{q}}_{2s}\hat{z})\right],$

其中

是反射和透射系数。局部坐标系下各波水平慢度相等，垂直慢度由下式给出：

${\hat{q}}_{\mathrm{ip}}=\sqrt{{\mathrm{\α}}_i^{-2}-{\hat{p}}_1^2}.$ ${\hat{q}}_{\mathrm{is}}=\sqrt{{\mathrm{\β}}_i^{-2}-{\hat{p}}_1^2}.(i=\mathrm{1,2})$

相对于局部坐标系，界面是水平的，这种情况的反射和透射系数计算方法可参照相关文献（Muller，1985）。将所有反射波和透射波由局部坐标系转换到全局坐标系，整理变换后各种波的表达式则得到全局坐标系下各反射波和透射波为：

${\mathrm{\φ}}_1^{\′}=R_{\mathrm{pp}}^d\·\exp \left[\mathrm{j\ω}(p_{1p}x-q_{1p}z)\right],$

${\mathrm{\ψ}}_1^{\′}=R_{\mathrm{ps}}^d\·\exp \left[\mathrm{j\ω}(p_{1s}x-q_{1s}z)\right],$

${\mathrm{\φ}}_2=T_{\mathrm{pp}}^d\·\exp \left[\mathrm{j\ω}(p_{2p}x+q_{2p}z)\right],$

${\mathrm{\ψ}}_2=T_{\mathrm{ps}}^d\·\exp \left[\mathrm{j\ω}(p_{2s}x+q_{2s}z)\right],$

式中：

$R_{\mathrm{pp}}^d={\hat{R}}_{\mathrm{pp}}^d\·A_1\exp [-\mathrm{j\ω}(p_{1p}{x}_0-q_{1p}{z}_0)],$

$R_{\mathrm{ps}}^d={\hat{R}}_{\mathrm{ps}}^d\·A_1\exp [-\mathrm{j\ω}(p_{1s}{x}_0-q_{1s}{z}_0)],$

$T_{\mathrm{pp}}^d={\hat{T}}_{\mathrm{pp}}^d\·A_1\exp [-\mathrm{j\ω}(p_{2p}{x}_0+q_{2p}{z}_0)],$

$T_{\mathrm{ps}}^d={\hat{T}}_{\mathrm{ps}}^d\·A_1\exp [-\mathrm{j\ω}(p_{2s}{x}_0+q_{2s}{z}_0)],$

$\left[\begin{array}{c}{p}_{1p}\\ q_{1p}\end{array}\right]=C\left[\begin{array}{c}{\hat{p}}_1\\ {\hat{q}}_{1p}\end{array}\right].$ $\left[\begin{array}{c}{p}_{1s}\\ q_{1s}\end{array}\right]=C\left[\begin{array}{c}{\hat{p}}_1\\ {\hat{q}}_{1s}\end{array}\right].$ $\left[\begin{array}{c}{p}_{2p}\\ q_{2p}\end{array}\right]=C^{-1}\left[\begin{array}{c}{\hat{p}}_2\\ {\hat{q}}_{2p}\end{array}\right].$ $\left[\begin{array}{c}{p}_{2s}\\ q_{2s}\end{array}\right]=C^{-1}\left[\begin{array}{c}{\hat{p}}_2\\ {\hat{q}}_{2s}\end{array}\right].$

其中

即分别为全局坐标系XOZ下粘弹性倾斜界面的P-P反射系数、P-S反射系数、P-P透射系数和P-S透射系数。

对于波由下方向上入射到界面的情况，设全局坐标系下有一P波φ′₂在介质II中向上传播：

φ′₂=exp[jω(p₂x-q₂z)].

则在界面处会产生反射回介质II的下行P波φ₂和SV波Ψ₂，以及透射到介质I中的上行P波φ′₁和SV波Ψ′₁。按照P波从界面上方入射时的处理流程：先将入射波φ′₂转换到局部坐标系下，在局部坐标系里求得各反射波和透射波的传统反射和透射系数；再将局部坐标系下各反射波和透射波转换到全局坐标系，得到全局坐标系下反射波和透射波的表达式分别为：

${\mathrm{\φ}}_2=R_{\mathrm{pp}}^u\·\exp \left[\mathrm{j\ω}(p_{2p}x+q_{2p}z)\right],$

${\mathrm{\ψ}}_2=R_{\mathrm{ps}}^u\·\exp \left[\mathrm{j\ω}(p_{2s}x+q_{2s}z)\right],$

${\mathrm{\φ}}_1^{\′}=T_{\mathrm{pp}}^u\·\exp \left[\mathrm{j\ω}(p_{1p}x-q_{1p}z)\right],$

${\mathrm{\ψ}}_1^{\′}=T_{\mathrm{ps}}^u\·\exp \left[\mathrm{j\ω}(p_{1s}x-q_{1s}z)\right],$

其中：

$R_{\mathrm{pp}}^u={\hat{R}}_{\mathrm{pp}}^u\·A^{\′}\exp [-\mathrm{j\ω}(p_{2p}{x}_0+q_{2p}{z}_0)],$

$R_{\mathrm{ps}}^u={\hat{R}}_{\mathrm{ps}}^u\·A^{\′}\exp [-\mathrm{j\ω}(p_{2s}{x}_0+q_{2s}{z}_0)],$

$T_{\mathrm{pp}}^u={\hat{T}}_{\mathrm{pp}}^u\·A^{\′}\exp [-\mathrm{j\ω}(p_{1p}{x}_0-q_{1p}{z}_0)],$

$T_{\mathrm{ps}}^u={\hat{T}}_{\mathrm{ps}}^u\·A^{\′}\exp [-\mathrm{j\ω}(p_{1s}{x}_0-q_{1s}{z}_0)],$

A′=exp[jω(p₂x₀-q₂z₀)]

$\left[\begin{array}{c}{p}_{2p}\\ q_{2p}\end{array}\right]=C^{-1}\left[\begin{array}{c}{\hat{p}}_2\\ {\hat{q}}_{2p}\end{array}\right],$ $\left[\begin{array}{c}{p}_{2s}\\ q_{2s}\end{array}\right]=C^{-1}\left[\begin{array}{c}{\hat{p}}_2\\ {\hat{q}}_{2s}\end{array}\right],$ $\left[\begin{array}{c}{p}_{1p}\\ q_{1p}\end{array}\right]=C\left[\begin{array}{c}{\hat{p}}_1\\ {\hat{q}}_{1p}\end{array}\right],$ $\left[\begin{array}{c}{p}_{1s}\\ q_{1s}\end{array}\right]=C\left[\begin{array}{c}{\hat{p}}_1\\ {\hat{q}}_{1s}\end{array}\right],$

${\hat{q}}_{\mathrm{ip}}=\sqrt{{\mathrm{\α}}_i^{-2}-{\hat{p}}_1^2},$ ${\hat{q}}_{\mathrm{is}}=\sqrt{{\mathrm{\β}}_i^{-2}-{\hat{p}}_1^2},$ (i=1，2) ${\hat{p}}_1={\hat{p}}_2,$ $\left[\begin{array}{c}{\hat{p}}_2\\ {\hat{q}}_2\end{array}\right]=C^{-1}\left[\begin{array}{c}{p}_2\\ q_2\end{array}\right].$

其中

即分别为全局坐标系下倾斜界面的P-P反射系数、P-S反射系数、P-P透射系数和P-S透射系数。

同理，对于SV波从倾斜界面上方（或下方）入射的情况，其求解与P波入射时相同。

上述步骤3）中GFM积分方法为：

为了利用GFM积分法，先将步骤3）中积分式改写成如下形式(这里只给出求和号里的一项)：

P=∫_Γf(p)exp[sg(p)]dp

其中：

$f\left(p\right)=\frac{{-R}_i\left(p\right)}{2\mathrm{\π}\·2\mathrm{jq}},$ s=jω·max(x,z), $g\left(p\right)=\frac{p_{1i}x-q_{1i}z}{\max (x,z)}.$

GFM梯形法则为：

$\underset{a}{\overset{b}{\∫}}f\left(p\right)\exp \left[\mathrm{sg}\left(p\right)\right]\mathrm{dp}=\left\{\begin{array}{cc}\frac{\mathrm{\δp}}{\mathrm{s\δg}}[\mathrm{\δ}\left({\mathrm{fe}}^{\mathrm{sg}}\right)-\frac{\mathrm{\δ}\left(f\right)\mathrm{\δ}\left(e^{\mathrm{sg}}\right)}{\mathrm{s\δ}\left(g\right)}]& \mathrm{\δg}\≠0\\ \frac{\mathrm{\δp}}{2}[f\left(a\right)e^{\mathrm{sg}\left(a\right)}+f\left(b\right)e^{\mathrm{sg}\left(b\right)}]& \mathrm{\δg}=0\end{array}\right..$

这里b-a等于积分步长，对任意函数h，δh=h(b)-h(a)。

Claims (4)

1.一种倾斜层状粘弹性介质中波场正演模拟方法，其特征在于，包括以下步骤：

1）将点源激励的球面波分解为平面波；

2）将平面波在倾斜层状粘弹性介质中传播，遇到界面时进行如下处理：以该界面为

轴，界面与全局坐标系Z轴的交点为原点建立局部坐标系，

垂直于界面向下，将全局坐标系中的入射波转换到局部坐标系中；在局部坐标系下求得反射波和透射波；再将局部坐标系下的反射波和透射波转换到全局坐标系下，得到全局坐标系下的反射系数和透射系数，以及反射波和透射波；在每个界面上根据需要选取波的类型，P波入射时能够选取的波包括反射PP波、反射PS波、透射PP波、透射PS波；S波入射时能够选取的波包括反射SP波、反射SS波、透射SP波、透射SS波；从而模拟所选取的波类型，得到模拟的平面波；

3）由步骤2）模拟得到的平面波根据下式合成点源激发的波场：

$P=\frac{1}{2\mathrm{\π}}{\∫}_{\mathrm{\Γ}}\underset{i=1}{\overset{N-1}{\mathrm{\Σ}}}[R_i\left(p\right)\·\frac{\exp \left[\mathrm{j\ω}(p_{1i}x-q_{1i}z)\right]}{-2\mathrm{jq}}]\mathrm{dp}$

式中，Γ为积分路径；N为平面波的个数；R_i(p)为平面波在多个界面上的透射和反射系数综合；j是虚数单位，ω是角频率；p_1i、q_1i分别为第i个界面反射回第1层的平面波的水平和垂直慢度；p、q分别为点源处平面波的水平和垂直慢度，x和z表示位置坐标。

2.根据权利要求1所述的倾斜层状粘弹性介质中波场正演模拟方法，其特征在于，步骤2）中，模拟震源为P波时的一次反射P波，设震源和检波器均在第一层，将震源处入射波逐层向下传播的递推步骤为：

①每遇到一个界面i，即以该界面为

轴，界面与全局坐标系Z轴的交点为原点建立局部坐标系垂直于界面向下；将全局坐标系中的入射波φ_i转换到局部坐标系中得

②在局部坐标系下求得反射系数

和透射系数

以及反射波

和透射波

③再将局部坐标系下的反射波和透射波转换到全局坐标系下，得到全局坐标系下的反射系数

和透射系数以及反射波φ′_i，i和透射波φ_i+1；

④反射波φ′_i，i逐层向上递推到第一层，遇到界面时执行①-③步，但只取透射波；

⑤透射波φ_i+1继续向下递推，遇到下一个界面时回到步骤①。

3.根据权利要求1所述的倾斜层状粘弹性介质中波场正演模拟方法，其特征在于，步骤2）中，模拟其它类型的波采用与模拟震源为P波时的一次反射P波的相同方法和步骤进行。

4.根据权利要求1所述的倾斜层状粘弹性介质中波场正演模拟方法，其特征在于，步骤3）中，式中积分方法采用Frazer和Gettrust提出的GFM积分法。

Images