CN103064115A - 一种射线参数域纵波与转换波匹配方法 - Google Patents

Abstract

一种射线参数域纵波与转换波匹配方法，包括：对由地震采集得到并经过偏移成像的纵波与转换波CIP道集进行射线追踪；在射线参数域的纵波与转换波CIP道集上标定和拾取对应同相轴的双程旅行时t和τ计算出初始的纵横波速度比模型γ(t)；对纵波与转换波CIP道集叠加形成纵波和转换波的地震道，在γ(t)的基础上加入预定范围的扰动量形成压缩比函数，将转换波地震道压缩到纵波时间，并将压缩后的地震道振幅进行重构；根据能量谱的极值点位置修正所述的γ(t),然后根据修正后的γ(t)将转换波数据匹配到纵波时间得到匹配数据；利用Gabor变换根据匹配数据产生转换波的时频谱，对每个Gabor变换时窗内的频谱一一进行修正，并利用反Gabor变换生成在时间t的转换波地震数据。

Description

一种射线参数域纵波与转换波匹配方法

技术领域

本发明是关于地球物理勘探技术，特别是关于一种射线参数域纵波与转换波匹配方法。

背景技术

与常规的纵波反射地震勘探技术不同，多分量反射波地震勘探技术利用三分量检波器接收到常规纵波和转换横波，因此能够更有效的解决的地质问题，例如：通过研究多分量地震的各向异性提高地下地质体的成像质量、改善裂缝储藏的描述、通过联合反演和联合及时更有效的对地下介质进行岩性识别和流体预测。纵波数据与转换波数据在时间、空间和波形上都存在差异。因此在多分量地震数据的联合成像、反演和解释中，两种数据的匹配是一项必不可少的关键技术。

目前常规的纵波与转换波的匹配技术的重点是解决两种数据之间存在的时间差异，即将转换波数据从原本的时间域合理的转换到纵波时间域，但是这种做法忽视了两种数据之间的空间差异。而即便在解决两种数据的时间差异的过程中，目前的匹配技术在很大程度上依赖人工作业，这与地震数据处理流程中的动校正速度分析类似。通常通过人工拾取纵波与转换波反射剖面中相对应的主要层位，再以此为依据通过相关分析估算转换波的压缩比，将转换波数据压缩到纵波数据时间。这样虽然能够将两种数据的主要层位进行匹配，但是无法保证两种数据中每个反射层位在时间上的一一对应。另外，当转换波数据被压缩到纵波时间后，原有的波形与频率都发生了非线性畸变，因而导致后续的联合反演、联合解释的可靠性降低。因此纵波与转换波的高精度匹配和后续的保幅处理成为多波地震技术中需要解决的问题。

发明内容

本发明提供一种射线参数域纵波与转换波匹配方法，以在射线参数域进行多波数据高精度的保幅匹配。

为了实现上述目的，本发明提供一种射线参数域纵波与转换波匹配方法，该方法包括：对由地震采集得到并经过偏移成像的纵波与转换波CIP道集进行射线追踪，从偏移距域转换到射线参数域；在所述射线参数域的纵波与转换波CIP道集上标定和拾取对应同相轴的双程旅行时t和τ，并根据t和τ计算出初始的纵横波速度比模型γ(t)；对每一所述纵波与转换波CIP道集叠加形成纵波和转换波的地震道，在所述纵横波速度比模型γ(t)的基础上加入预定范围的扰动量形成一系列新的压缩比函数，利用所述新的压缩比函数将转换波地震道压缩到纵波时间，并将压缩后的地震道振幅进行重构；由压缩及重构的所述转换波及纵波地震道得到能量谱，根据所述能量谱的极值点位置修正所述的纵横波速度比模型γ(t),然后根据修正后的所述的纵横波速度比模型将转换波数据匹配到纵波时间得到匹配数据；利用Gabor变换根据所述匹配数据产生转换波的时频谱，对每个Gabor变换时窗内的频谱一一进行修正，并利用反Gabor变换生成在时间t的转换波地震数据。

进一步地，对于地震反射能量较弱、相关谱结果无法清晰指示的目的层区域，在得到所述匹配数据的基础上对局部速度比函数的梯度进行扫描。

进一步地，进行扫描的局部速度比函数的梯度范围是-1.5到1.5，速度比扰动量从-0.2到0.2。

进一步地，所述的射线追踪为弯曲射线追踪，纵波及转换波的双程旅行时T_pp及T_PS分别为：

$T_{\mathrm{pp}}=\underset{i=0}{\overset{k}{\mathrm{\Σ}}}(u_{p(i-1)}+u_{p\left(i\right)}){\mathrm{dl}}_i$

$T_{\mathrm{PS}}=\frac{1}{2}[\underset{i=0}{\overset{k}{\mathrm{\Σ}}}(u_{p\left(i\right)}+u_{p(i+1)}){\mathrm{dl}}_i+\underset{i=k}{\overset{2k}{\mathrm{\Σ}}}(u_{s\left(i\right)}+u_{s(i+1)}){\mathrm{dl}}_i]$

其中dl_i是在水平坐标x_i和x_i-1之间的射线路径，u_p(i)和u_s(i)分别是第i层介质的纵波和横波慢度：其中u_p(i)＝u_p(2k-i)u_s(i)＝u_s(2k-i)，i=1,Λ,2k-1。

进一步地，根据t和τ计算出初始的纵横波速度比模型γ(t)的依据是：

$\mathrm{\γ}\left(t\right)=\frac{2\mathrm{\τ}-t}{t}.$

进一步地，将压缩后的地震道振幅进行重构的依据为sinc方程：

$u_{\mathrm{PS}}\left(t\right)=\underset{k}{\mathrm{\Σ}}{u}_{\mathrm{PS}}\left(\mathrm{k\Δ\τ}\right)\sin c(\frac{\mathrm{\τ}\left(t\right)}{\mathrm{\Δ\τ}}-k)$

其中u_PS(τ)是在原时间域τ的转换波，u_PS(t)是在纵波时间域的转换波的转换波，Δτ为u_PS(τ)在原时间域的采样间隔，k=-n,-n+1,…,0,…,n-1，通常n>=4（利用sinc方程重构t时间的Ups需要Ups在原来时间tau附近k点的数据）。

进一步地，由压缩及重构的所述转换波及纵波地震道得到能量谱包括：根据下述公式由压缩及重构的所述转换波及纵波地震道得到能量谱：

$c\left(t\right)=\frac{\underset{\mathrm{\ξ}=-T/2}{\overset{T/2}{\mathrm{\Σ}}}{u}_{\mathrm{PS}}(t+\mathrm{\ξ})u_{\mathrm{PP}}(t+\mathrm{\ξ})}{{\left[\underset{\mathrm{\ξ}=-T/2}{\overset{T/2}{\mathrm{\ξ}}}{u}_{\mathrm{PS}}^2(t+\mathrm{\ξ})\underset{\mathrm{\ξ}=-T/2}{\overset{T/2}{\mathrm{\Σ}}}{u}_{\mathrm{PP}}^2(t+\mathrm{\ξ})\right]}^{1/2}}$

其中，u_PP(t+ξ)为纵波振幅，为了增加信号信噪比，在地震到某一时间t附近的小时窗内做相关，时窗长度为T，ξ=-T/2,…,0，…,T/2。

进一步地，所述对每个Gabor变换时窗内的频谱一一进行修正，包括：

在某一以t_c为中心的Gabor变换时窗内根据所述纵横波速度比值γ(t)计算出波形压缩量β(t)，进而根据所述波形压缩量β(t)修正的局部时频谱u_PS(ω(t_c)),将u_PS(ω(t_c))恢复成在对应的原始转换波时间τ_c=t_c/α(t_c)的频谱u_PS(ω(τ_c))；

其中， $\mathrm{\β}\left(t\right)\≡\frac{\mathrm{dt}}{\mathrm{d\τ}}=\frac{2}{1+\mathrm{\γ}\left(t\right)+\mathrm{t\γ}\text{'}\left(t\right)},$

u_PS(ω(t_c))→|β(t_c)|u_PS(β(t_c)ω(t_c))=u_PS(ω(τ_c))。

进一步地，所述预定范围的扰动量为：-0.4到0.4的扰动量。

本发明实施例的有益效果在于，本发明可用于在多波地震勘探数据处理、反演和解释中，纵波与转换波数据在时间和空间上的高精度匹配,并解决了匹配后的转换波的振幅保真的问题。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：

图1为本发明实施例射线参数域纵波与转换波匹配方法流程图；

图2为本发明实施例道集及纵横波速度表示意图；

图3a为本发明实施例与图2中(a)和(c)所示的CIP道集对应的叠后纵波与转换波地震道示意图；

图3b为本发明实施例利用在初始速度比模型上加入一定范围的扰动形成的一系列γ(t)，将转换波地震道压缩到纵波时间的示意图；

图3c为本发明实施例的纵波与转换波的相关谱示意图；

图3d为本发明实施例根据相关谱修正原纵横波速度比模型示意图；

图4a为本发明实施例转换波从原始时间域τ的转换到纵波时间域t波形畸变示意图；

图4b为本发明实施例转换波从原始时间域τ的转换到纵波时间域t频谱畸变示意图

图5a为本发明实施例纵波过井剖面示意图；

图5b为本发明实施例转换波过井剖面示意图；

图5c为本发明实施例匹配后的转换波剖面示意图；

图6a为本发明实施例图5a至图5c的剖面中的井旁地震道示意图；

图6b为本发明实施例利用一系列γ(t)压缩得到的转换波与纵波地震道的相关谱示意图；

图6c为本发明实施例根据相关谱修正原纵横波速度比模型示意图；

图6d为本发明实施例对2.2s附近的局部速度比函数的梯度进行扫描得到的二维谱示意图

图7a至图7d为通过本发明的射线参数域纵波与转换波匹配方法转换波及保幅后的效果示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。

如图1所示，本发明实施例提供一种射线参数域纵波与转换波匹配方法，该射线参数域纵波与转换波匹配方法包括：

步骤S101：对由地震采集得到并经过偏移成像的纵波与转换波CIP道集进行射线追踪，从偏移距域转换到射线参数域。这样保证了沿着某一条射线参数为常数的射线路径的PP波和PS波的地震反射记录来自于地下同一空间点的反射。

为了兼顾效率和精度，本发明的射线追踪法采用在地震层析中普遍使用的方法—弯曲射线追踪法（Wang，2003）。在每个共成像点道集中，速度场仅沿深度方向变化，待计算的路径可以离散为具有固定深度间隔的折线路径。则纵波和转换波的双程旅行时分别为：

$T_{\mathrm{pp}}=\underset{i=0}{\overset{k}{\mathrm{\Σ}}}(u_{p(i-1)}+u_{p\left(i\right)}){\mathrm{dl}}_i$

$T_{\mathrm{PS}}=\frac{1}{2}[\underset{i=0}{\overset{k}{\mathrm{\Σ}}}(u_{p\left(i\right)}+u_{p(i+1)}){\mathrm{dl}}_i+\underset{i=k}{\overset{2k}{\mathrm{\Σ}}}(u_{s\left(i\right)}+u_{s(i+1)}){\mathrm{dl}}_i]---\left(1\right)$

其中dl_i是在水平坐标x_i和x_i-1之间的射线路径，u_p(i)和u_s(i)分别是第i层介质的纵波和横波慢度：其中u_p(i)＝u_p(2k-i)u_s(i)＝u_s(2k-i)，i=1,Λ,2k-1。

根据费马原理，射线路径是使旅行时T最小的路径，因此公式(1)可表示为：

$\frac{\∂T_{\mathrm{PP}}}{\∂x_i}=\frac{(u_{p\left(i\right)}+u_{p(i+1)})}{{\mathrm{dl}}_{i+1}}(x_i-x_{i+1})+\frac{(u_{p\left(i\right)}+u_{p(i-1)})}{{\mathrm{dl}}_i}(x_i-x_{i-1})=0,$

$\frac{\∂T_{\mathrm{PS}}}{\∂x_i}=\left\{\begin{array}{c}\frac{1}{2}[\frac{u_{p\left(i\right)}+u_{p(i+1)}}{{\mathrm{dl}}_{i+1}}(x_i-x_{i+1})+\frac{u_{p\left(i\right)}+u_{p(i-1)}}{{\mathrm{dl}}_i}(x_i-x_{i-1})]=0,i=1,\mathrm{\Λ},k\\ \frac{1}{2}[\frac{u_{s\left(i\right)}+u_{s(i+1)}}{{\mathrm{dl}}_{i+1}}(x_i-x_{i+1})+\frac{u_{s\left(i\right)}+u_{s(i-1)}}{{\mathrm{dl}}_i}(x_i-x_{i-1})]=0,i=k,\mathrm{\Λ},2k-1\end{array}\right.---\left(2\right)$

式中每段射线路径对应的水平坐标x_i和x_i-1可通过解上述线性方程得到。

步骤S102：在所述射线参数域的纵波与转换波CIP道集上标定和拾取对应同相轴（即同一反射层）的双程旅行时t和τ，并根据t和τ计算出初始的纵横波速度比模型γ(t)，

图2为本发明实施例道集及纵横波速度表示意图，如图2所示，(a)表示射线参数域纵波共成像点(CIP)道集，(b)表示在转换波时间域τ的纵横波速度比γ(τ)，(c)表示射线参数域转换波共成像点(CIP)，(d)表示在纵换波时间域t的纵横波速度比γ(t)，(e)表示转换到纵波时间的转换波CIP道集。

步骤S103：对每一所述纵波与转换波CIP道集叠加形成纵波和转换波的地震道，如图3a所示，图3a为与图2中(a)和(c)所示的CIP道集对应的叠后纵波与转换波地震道示意图。在所述纵横波速度比模型γ(t)的基础上加入预定范围的扰动量形成一系列新的压缩比函数，利用所述新的压缩比函数将转换波地震道压缩到纵波时间，并利用sinc方程将压缩后的地震道振幅进行重构，如图3b所示，为利用在初始速度比模型上加入一定范围的扰动形成的一系列γ(t)，将转换波地震道压缩到纵波时间的示意图。sinc方程为：

$u_{\mathrm{PS}}\left(t\right)=\underset{k}{\mathrm{\Σ}}{u}_{\mathrm{PS}}\left(\mathrm{k\Δ\τ}\right)\sin c(\frac{\mathrm{\τ}\left(t\right)}{\mathrm{\Δ\τ}}-k)---\left(3\right)$

其中u_PS(τ)是在原时间域τ的转换波，u_PS(t)是在纵波时间域的转换波的转换波，Δτ为u_PS(τ)在原时间域的采样间隔，k=-n,-n+1，…,0，…,n-1，通常n>=4（利用sinc方程重构t时间的Ups需要Ups在原来时间tau附近k点的数据）。

该步骤中，预定范围的扰动量为：-0.4到0.4的扰动量。

步骤S104：由压缩及重构的所述转换波及纵波地震道得到能量谱，如图3c所示的纵波与转换波的相关谱示意图。然后，根据所述能量谱的极值点位置修正所述的纵横波速度比模型γ(t),然后根据修正后的所述的纵横波速度比模型将转换波数据匹配到纵波时间得到匹配数据。图3d为根据相关谱修正原纵横波速度比模型示意图，图3d中的虚线为原速度比模型，实线为修正后的速度比模型，优化的速度比模型可用来将转换波数据压缩到纵波时间。

能量谱可以由下述公式（4）得到：

$c\left(t\right)=\frac{\underset{\mathrm{\ξ}=-T/2}{\overset{T/2}{\mathrm{\Σ}}}{u}_{\mathrm{PS}}(t+\mathrm{\ξ})u_{\mathrm{PP}}(t+\mathrm{\ξ})}{{\left[\underset{\mathrm{\ξ}=-T/2}{\overset{T/2}{\mathrm{\ξ}}}{u}_{\mathrm{PS}}^2(t+\mathrm{\ξ})\underset{\mathrm{\ξ}=-T/2}{\overset{T/2}{\mathrm{\Σ}}}{u}_{\mathrm{PP}}^2(t+\mathrm{\ξ})\right]}^{1/2}}---\left(4\right)$

其中，u_PP(t+ξ)为纵波振幅，为了增加信号信噪比，在地震到某一时间t附近的小时窗内做相关，时窗长度为T，ξ=-T/2,…,0，…,T/2。

步骤S105：利用Gabor变换根据所述匹配数据产生转换波的时频谱，对每个Gabor变换时窗内的频谱一一进行修正，并利用反Gabor变换生成在时间t的转换波地震数据。

该步骤中，对每个Gabor变换时窗内的频谱一一进行修正，具体包括：在某一以t_c为中心的Gabor变换时窗内根据所述纵横波速度比值γ(t)计算出波形压缩量β(t)，如图4a所示，为转换波从原始时间域τ的转换到纵波时间域t波形畸变示意图。进而根据所述波形压缩量β(t)修正的局部时频谱u_PS(ω(t_c)),将u_PS(ω(t_c))恢复成在对应的原始转换波时间τ_c=t_c/α(t_c)的频谱u_PS(ω(τ_c))，如图4b所示，为转换波从原始时间域τ的转换到纵波时间域t频谱畸变示意图，其中，

$\mathrm{\β}\left(t\right)\≡\frac{\mathrm{dt}}{\mathrm{d\τ}}=\frac{2}{1+\mathrm{\γ}\left(t\right)+\mathrm{t\γ}\text{'}\left(t\right)},$

u_PS(ω(t_c))→|β(t_c)|u_PS(β(t_c)ω(t_c))=u_PS(ω(τ_c))。

图5a至图5c及图6a至图6d图为实际地震数据的匹配示例。图5a为纵波过井剖面示意图。图5b为转换波过井剖面示意图。图5c为匹配后的转换波剖面示意图。图6a为图5剖面中的井旁地震道示意图，从左至右：纵波、转换波、匹配后的转换波、和保幅匹配后的转换波地震道。图6b为利用一系列γ(t)压缩得到的转换波与纵波地震道的相关谱示意图。图6c为根据相关谱修正原纵横波速度比模型示意图。图6d为对2.2s附近的局部速度比函数的梯度进行扫描得到的二维谱示意图。

对于地震反射能量较弱、相关谱结果无法清晰指示的目的层区域，可以在步骤S104之后，在获得的匹配数据的基础上，利用公式（5）对局部速度比函数的梯度进行扫描。较佳地，扫描的局部速度比函数的梯度范围可以是-1.5到1.5，速度比扰动量从-0.2到0.2，从而达到精细匹配的目的。

$\mathrm{\γ}\left(t\right)={\mathrm{\γ}}_c+\mathrm{\Δ\γ}+\frac{\mathrm{d\γ}}{\mathrm{dt}}(t-t_c)---\left(5\right)$

其中，t_c是某一待分析的时窗中点，γ_c≡γ(t_c)为一速度比常值，Δγ是在γ_c基础上的扰动量,dγ/dt是该时窗内的局部速度比函数。

图7a至图7d为通过本发明的射线参数域纵波与转换波匹配方法转换波及保幅后的效果示意图。图7a中从左至右：纵波、转换波、匹配后的转换波、和保幅匹配后的转换波地震道示意图。图7b为转换波匹配前（虚线）后（实线）的振幅谱示意图。图7c为转换波匹配前（虚线）的振幅谱与保幅处理后的振幅谱（实线示意图）。图7d为图7a中第二个子波的放大比较图。

本发明实施例的有益效果在于，本发明提供的匹配方法在射线参数域进行。由于PP波和PS波地震道集是沿着从震源点到反射点的某一条射线参数为常数的射线路径的地震反射记录，保证了这两种波对应的记录来自于地下同一空间点的反射。本发明可用于在多波地震勘探数据处理、反演和解释中，纵波与转换波数据在时间和空间上的高精度匹配,并解决了匹配后的转换波的振幅保真的问题。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种射线参数域纵波与转换波匹配方法，其特征在于，所述的方法包括：

对由地震采集得到并经过偏移成像的纵波与转换波CIP道集进行射线追踪，从偏移距域转换到射线参数域；

在所述射线参数域的纵波与转换波CIP道集上标定和拾取对应同相轴的双程旅行时t和τ，并根据t和τ计算出初始的纵横波速度比模型γ(t)；

对每一所述纵波与转换波CIP道集叠加形成纵波和转换波的地震道，在所述纵横波速度比模型γ(t)的基础上加入预定范围的扰动量形成一系列新的压缩比函数，利用所述新的压缩比函数将转换波地震道压缩到纵波时间，并将压缩后的地震道振幅进行重构；

由压缩及重构的所述转换波及纵波地震道得到能量谱，根据所述能量谱的极值点位置修正所述的纵横波速度比模型γ(t),然后根据修正后的所述的纵横波速度比模型将转换波数据匹配到纵波时间得到匹配数据；

利用Gabor变换根据所述匹配数据产生转换波的时频谱，对每个Gabor变换时窗内的频谱一一进行修正，并利用反Gabor变换生成在时间t的转换波地震数据。

2.根据权利要求1所述的射线参数域纵波与转换波匹配方法，其特征在于，对于地震反射能量较弱、相关谱结果无法清晰指示的目的层区域，在得到所述匹配数据的基础上对局部速度比函数的梯度进行扫描。

3.根据权利要求2所述的射线参数域纵波与转换波匹配方法，其特征在于，进行扫描的局部速度比函数的梯度范围是-1.5到1.5，速度比扰动量从-0.2到0.2。

4.根据权利要求2或3所述的射线参数域纵波与转换波匹配方法，其特征在于，所述的射线追踪为弯曲射线追踪，纵波及转换波的双程旅行时T_pp及T_PS分别为：

$T_{\mathrm{pp}}=\underset{i=0}{\overset{k}{\mathrm{\Σ}}}(u_{p(i-1)}+u_{p\left(i\right)}){\mathrm{dl}}_i$

$T_{\mathrm{PS}}=\frac{1}{2}[\underset{i=0}{\overset{k}{\mathrm{\Σ}}}(u_{p\left(i\right)}+u_{p(i+1)}){\mathrm{dl}}_i+\underset{i=k}{\overset{2k}{\mathrm{\Σ}}}(u_{s\left(i\right)}+u_{s(i+1)}){\mathrm{dl}}_i]$

其中dl_i是在水平坐标x_i和x_i-1之间的射线路径，u_p(i)和u_s(i)分别是第i层介质的纵波和横波慢度：其中u_p(i)＝u_p(2k-i)u_s(i)＝u_s(2k-i)，i=1,Λ,2k-1。

5.根据权利要求4所述的射线参数域纵波与转换波匹配方法，其特征在于，根据t和τ计算出初始的纵横波速度比模型γ(t)的依据是：

$\mathrm{\γ}\left(t\right)=\frac{2\mathrm{\τ}-t}{t}.$

6.根据权利要求5所述的射线参数域纵波与转换波匹配方法，其特征在于，将压缩后的地震道振幅进行重构的依据为sinc方程：

$u_{\mathrm{PS}}\left(t\right)=\underset{k}{\mathrm{\Σ}}{u}_{\mathrm{PS}}\left(\mathrm{k\Δ\τ}\right)\sin c(\frac{\mathrm{\τ}\left(t\right)}{\mathrm{\Δ\τ}}-k)$

其中u_PS(τ)是在原时间域τ的转换波，u_PS(t)是在纵波时间域的转换波的转换波，Δτ为u_PS(τ)在原时间域的采样间隔，k=-n,-n+1,…,0,…,n-1，通常n>=4。

7.根据权利要求6所述的射线参数域纵波与转换波匹配方法，其特征在于，由压缩及重构的所述转换波及纵波地震道得到能量谱包括：根据下述公式由压缩及重构的所述转换波及纵波地震道得到能量谱：

$c\left(t\right)=\frac{\underset{\mathrm{\ξ}=-T/2}{\overset{T/2}{\mathrm{\Σ}}}{u}_{\mathrm{PS}}(t+\mathrm{\ξ})u_{\mathrm{PP}}(t+\mathrm{\ξ})}{{\left[\underset{\mathrm{\ξ}=-T/2}{\overset{T/2}{\mathrm{\ξ}}}{u}_{\mathrm{PS}}^2(t+\mathrm{\ξ})\underset{\mathrm{\ξ}=-T/2}{\overset{T/2}{\mathrm{\Σ}}}{u}_{\mathrm{PP}}^2(t+\mathrm{\ξ})\right]}^{1/2}}$

其中，u_PP(t+ξ)为纵波振幅；ξ=-T/2,…,0，…,T/2，T为时窗长度。

8.根据权利要求7所述的射线参数域纵波与转换波匹配方法，其特征在于，所述对每个Gabor变换时窗内的频谱一一进行修正，包括：

在某一以t_c为中心的Gabor变换时窗内根据所述纵横波速度比值γ(t)计算出波形压缩量β(t)，进而根据所述波形压缩量β(t)修正的局部时频谱u_PS(ω(t_c)),将u_PS(ω(t_c))恢复成在对应的原始转换波时间τ_c=t_c/α(t_c)的频谱u_PS(ω(τ_c))；

其中， $\mathrm{\β}\left(t\right)\≡\frac{\mathrm{dt}}{\mathrm{d\τ}}=\frac{2}{1+\mathrm{\γ}\left(t\right)+\mathrm{t\γ}\text{'}\left(t\right)},$

u_PS(ω(t_c))→|β(t_c)|u_PS(β(t_c)ω(t_c))=u_PS(ω(τ_c))。

9.根据权利要求6所述的射线参数域纵波与转换波匹配方法，其特征在于，所述预定范围的扰动量为：-0.4到0.4的扰动量。

Images