一种提取Rayleigh面波频散曲线的方法
【技术领域】
[0001] 本发明涉及工程地震勘探或地球深部构造探测领域中的Rayleigh面波探测方法 的关键技术，即通过反演地下横波速度结构时Rayleigh面波的频散曲线的提取，我利用地 震记录互相关谱实部曲线和第一类零阶贝塞尔函数的关系，提取出了Rayleigh面波的频 散曲线。
【背景技术】
[0002] Rayleigh面波是一种沿地球自由表面传播的弹性波，其沿地表传播的速度与波 及深度内的介质弹性参数有关，包括介质的密度、纵波和横波速度，而主要的影响参数是 介质的横波速度。Rayleigh波在层状介质中具有频散特征，其传播速度将随着频率的改 变而发生变化，根据Rayleigh面波的这一特征，可以求得地层不同深度的弹性参数。
[0003] 频散曲线的提取Rayleigh面波勘探中的关键步骤，简单来说就是求取对应某频 率下瑞雷面波的相速度问题。传统的提取频散曲线方法的思想为，对不同检波器或地震台 站实际接收到的RayIeigh面波信号做fourier变换得到其相位，计算不同检波器或台站处 的相位差，再根据检波器或台站间的相对距离即可求出某频点的相速度，假设Pl台站处的 Rayleigh面波信号为Rl(t)，p2台站处的Rayleigh面波信号为R2 (t)，分别做傅里叶变换 的到Rl和R2的相位谱Rl(f，(6)和R2(f，(6)，若两台站间的距离为S，则某个频率处的相 速度的可由下式算得：
【主权项】
1. 一种提取频散曲线的方法，其特征在于，所述方法包括： a数据采集： 在研究区域布设台站对或台阵，进行长时背景噪声观测，根据目标探测深度选取合适 的台站间距和采样率； b数据预处理： 对得到的背景噪声信号进行谱白化、滤波等预处理过程，然后将各个台站的微动信号 分段，分段时长根据探测目标深度选取； C互相关运算： 将分段后的信号按相同时间段内两两进行互相关运算，然后将相同台站对的不同时间 段互相关运算的结果进行叠加，叠加后的互相关函数即可视为台站间的经验格林函数； d提取谱实部曲线： 对经验格林函数进行傅里叶变换得到互相关谱，然后提取出谱实部曲线； e计算相速度： 实际数据中计算得到的实部曲线与第一类零阶贝塞尔函数振幅有差异，但是实部曲 线的零点与贝塞尔函数的根仍然一一对应，因此我们将零点处对应频率的相速度值一一求 出； f绘制频散曲线： 将得到的各个零点处对应频率的相速度值绘制在v-f平面坐标内，连接各个频点的相 速度值，即得到Ra^ei曲面波的频散曲线。
2. 根据权利要求1所述的Ra^ei曲面波频散曲线提取方法，其特征在于：所述方法通 过计算互相关函数谱的实部曲线，根据实部曲线零点与零阶贝塞尔函数根的对应关系计算 相速度。
【专利摘要】一种提取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面波的频散曲线。通过贝塞尔函数与互相关谱的关系，有效的提高了频散曲线的计算效率。
【IPC分类】G01V1-28
【公开号】CN104678435
【申请号】CN201410580207
【发明人】李欣欣, 方红萍
【申请人】李欣欣, 方红萍
【公开日】2015年6月3日
【申请日】2014年10月27日