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

权利要求书

1.一种提取被动源面波频散曲线的方法，其特征在于，所述方法包括：

获取目标侧线上的面波数据组；所述目标侧线上设有一个震源和N个节点；所述面波数据组包括：N组面波数据；一组面波数据对应一个所述节点；所述面波数据均包括：震源的位置、节点的位置、震源处的地震波的能量和节点处的地震波的能量；所述面波数据组是根据炮集组数据确定的；

对所述面波数据组进行频率域的转换，得到多组功率谱数据；所述功率谱数据，包括：震源的功率谱和节点的功率谱；

计算各组所述功率谱数据中所述震源的功率谱和所述节点的功率谱的互相关，得到多组地震记录数据；所述地震记录数据包括震源和节点间的相关地震波能量；

对各组所述地震记录数据进行时域的转换，得到多组地震数据；

将各组所述地震数据进行拉东变换以及频率域转换，得到多个面波频谱数据；每个所述面波频谱数据均包含一个频率和对应的多个所述相关地震波能量；

对于任一所述面波频谱数据，将所述面波频谱数据中多个所述相关地震波能量中的最大值，确定为所述面波频谱数据中所述频率对应的面波相速度的值；

根据多个所述面波相速度的值绘制面波频散曲线。

2.根据权利要求1所述的提取被动源面波频散曲线的方法，其特征在于，所述对所述面波数据组进行频率域的转换，得到多组功率谱数据，具体包括：

对所述面波数据组进行滤波，得到滤波后的面波数据组；

对所述滤波后的面波数据组进行傅里叶变换，得到多组功率谱数据。

3.根据权利要求1所述的提取被动源面波频散曲线的方法，其特征在于，所述互相关的计算公式为：

其中，表示地震记录数据的实部；为时间域格林函数G(x_A,x_B,ω)的傅里叶变换；ρ为质量密度；c为地震波的传播速度；为震源的功率谱；为节点的功率谱；·为空间的整体平均，为功率谱数据；为震源的功率谱的转置；ω为震动角速度；x_A为震源的位置；x_B为节点的位置。

4.根据权利要求1所述的提取被动源面波频散曲线的方法，其特征在于，所述拉东变换的计算公式为：

其中，ρ为斜率，τ为截距，u(t,x)为地震数据，u(τ,p)为对地震数据进行拉东变换后得到的数据。