[发明专利]一种提取Rayleigh面波频散曲线的方法

摘要

一种提取Rayleigh面波频散曲线的方法。其包括以下步骤：a、对相距为R的地震检波器对记录的噪声互相关信号进行傅里叶变换，得到频谱W(f)；b、对得到的频谱做Gauss窄带滤波得到WG(f)；c、提取频谱WG(f)的实部曲线并找出其零点，记fn为曲线上的第n个零点处的频率；d、找出第一类零阶贝塞尔函数的根记为Zn；e、根据互相关谱实部曲线与第一类零阶贝塞尔函数的关系求单频点处的Rayleigh面波相速度值V(f)=2πfnR/Zn；f、将算得的单频点相速度连接绘制为一条平滑曲线，做为Rayleigh面波的频散曲线。通过贝塞尔函数与互相关谱的关系，有效的提高了频散曲线的计算效率。

主权项

一种提取频散曲线的方法，其特征在于，所述方法包括：    a数据采集：    在研究区域布设台站对或台阵，进行长时背景噪声观测，根据目标探测深度选取合适的台站间距和采样率；    b数据预处理：    对得到的背景噪声信号进行谱白化、滤波等预处理过程，然后将各个台站的微动信号分段，分段时长根据探测目标深度选取；    c 互相关运算：    将分段后的信号按相同时间段内两两进行互相关运算，然后将相同台站对的不同时间段互相关运算的结果进行叠加，叠加后的互相关函数即可视为台站间的经验格林函数；    d 提取谱实部曲线：    对经验格林函数进行傅里叶变换得到互相关谱，然后提取出谱实部曲线；    e 计算相速度：    实际数据中计算得到的实部曲线与第一类零阶贝塞尔函数振幅有差异，但是实部曲线的零点与贝塞尔函数的根仍然一一对应，因此我们将零点处对应频率的相速度值一一求出；    f 绘制频散曲线：    将得到的各个零点处对应频率的相速度值绘制在V‑f平面坐标内，连接各个频点的相速度值，即得到Rayleigh面波的频散曲线。