一种多模式的频散曲线分析方法.docxVIP

一种多模式的频散曲线分析方法

1反演多模瑞利波方法的确定

瑞利波勘探方法是一种废弃的探测方法，用于获得横波速度信息。主要有两种瑞利波勘探方法，一种是静态激励。这种方法每次都会激发单个频率信号，其操作方法相对复杂且需要时间。另一种方法为瞬态激发,这种激发方式在单一冲击下能产生较宽的频带。通过两检波器之间的距离以及在各个频率下互功率谱的相位就能够得到相速度与频率(或波长)之间的关系,该方法被称为表面波谱分析方法(SASW),由于其简便快捷,在工程领域获得了广泛应用。

从实验频散曲线获取介质横波速度剖面有两种方法,一种简单的方法是依据横波速度大约是瑞利波相速度的110%,其有效穿透深度是瑞利波波长的三分之一或二分之一,这种方法在横波速度剖面中会引起一定程度的误差,有效穿透深度通常大于半波长;另一种对横波速度剖面进行解释的方法就是反演方法。不管是哪种方法,首先必须得到实验频散曲线,然后采用某种优化算法拟合理论与实验频散曲线,使两者之间的差异最小。因此,尽可能地构建正确的实验频散曲线是一个非常关键的步骤。

多层介质中的导波可分为两类:表面波和陷模波。表面波主要与自由边界的存在有关,而陷模波与自由边界下的分界面有关,其能量集中在自由边界下某一区域内。层状介质中存在多种导波(表面波和陷模波)模式,在不同频率范围内,不同的模式具有不同的幅度特性,所有这些模式通常都称为瑞利波模式。我们知道,由相关方法得到的频散曲线是所有模式共同作用的结果,这就导致了这样一个问题:实验频散曲线究竟对应于理论上的哪些模式?

在零频附近最先出现的模式通常被称为基阶模式(或第一阶模式),从低频到高频各模式分别被称为第二、第三、第四等高阶模式。通常来说,对于速度随深度递增的剖面,在所有的模式中基阶模式的能量起主导作用,这时可认为实验频散曲线只含有基阶模式,反演时可将基阶模式的频散曲线作为理论频散曲线。然而,对于不规则的速度深度剖面,即横波速度不是随着深度增加而递增的情形,实验频散曲线与基阶模式是不吻合的,这时,必须利用多模式来构建理论上的频散曲线。然而,在这些传统方法中,理论上的频散曲线是通过两个信号的互相关谱得到的,在这个过程中存在相位跳跃现象,即在相位上可以存在2π任意倍的差异,从而导致传播速度的不确定性,当频散曲线发生跳跃时,我们很难判断究竟是真的发生了模式跳跃现象还是由于相位发生了2π的跳跃。这时频散曲线可能会被错误地构建,相应地,得到的横波速度信息等其它结果都会是错误的。

瑞利波多道分析方法可以避免SASW方法中所面临的问题,这种方法可以得到独立的多模频散曲线,使得考虑包含高阶模式的多模瑞利波反演成为可能,这将改善反演的多解性问题,而且在某些情况下,由于基阶模式不能被激发,我们必须考虑高阶模式的瑞利波反演问题。有许多方法用于反演多道的瑞利波信号,提取多模瑞利波的频散曲线,比如频率波数(f-k)分析方法,τ-p变换法等。本文采用f-k方法,在分析层状介质中瑞利波模式的激发和传播特性基础上,针对三种典型的层状介质模型,采用数值模拟的方法研究了各模式的频散特性,详细地分析了有关参数(如接收道数目及源检距等)对各模式频散曲线的影响。

2波场参数选取

考虑均匀各向同性的弹性的水平层状半无界空间,对称点源位于第一层介质中激发P-SV声场。采用传递矩阵方法来计算层状半空间中的弹性波场及瑞利波模式分析。假定所有的接收点都位于水平界面上接收信号,而且接收到的信号都为垂直方向的位移分量。理论上的瑞利波频散曲线可以由传递矩阵方法很方便地得到,对于实验频散曲线,利用多道时域波形采取二维Fourier变换即可得到频率波数域中的表示,并可得到各模式的频散曲线(即实验频散曲线)。在实际激发过程中,声场是三维的,二维Fourier变换的结果具有一定的误差。然而,如果源检距足够大(远大于瑞利波波长),在接收点所得到的声场可以被当作平面波处理,这时二维Fourier变换得到的结果是可靠的。

选取三种类型的层状介质模型,它们的参数取值见表1,其中最后一列表示的是由第一层介质形成的半空间的瑞利波相速度。模型1是一个规则的剖面,横波速度随着深度的增加而增加,模型2中包含一个低速层,模型3是在模型1的基础上加了一个高速层而得到的,模型2和3都是不规则剖面,在大坝裂缝检测以及高速公路检测中经常出现。

所有的接收点在自由界面上成线性排列以接收波形,对每个波形,采样间隔是30μs,总的采样点数是4096。依据采样定律,接收点间的距离d必须小于最小半波长,即:

d≤λmin/2(1)

这意味着能包含在波形中的不发生混叠的最大波数为:

kNyquist=2π/(2d)=π/d(2)

接收点之间的间隔必须小于λmin/2,在f-k域中,导波是与波场最大值相联系的,其相速度可以通过v=2πf/k得到