[发明专利]一种瑞雷面波速度-波数分析方法

权利要求书

1.一种瑞雷面波速度-波数分析方法，其特征在于，包括如下步骤：

第一步，将含有m道，每道有n个采样点，时间采样率为                                                的单炮瑞雷面波地震记录读取到二维数组g中；

第二步，对g进行二维傅里叶变换，得到f-k谱，二维离散傅里叶正变换为：

       

  

其中，k为波数，为波数采样率，j₀=1，2，…，m，为波数采样序号；f为频率，为频率采样率，l₀=1，2，…，n，为频率采样序号；j=1，2，…，m，为道序号；l=1，2，…，n，为采样点序号；为道间距，为最小偏移距，x为偏移距；为时间采样率，t为采样时间； 为虚数单位；

第三步，建立传播速度（v）与频率（f）、波数（k）之间的关系，具体如下：

                                                

第四步，依据瑞雷面波的传播速度（v）与频率（f）、波数（k）之间的关系，将f-k谱中对应（k，f）坐标的能量映射到v-k谱中对应（k，v）坐标的位置，同时选定v-k谱输出的坐标范围并成图；

第五步，基于v-k谱，从左至右依次拾取v-k谱中的强能量团，并提取对应的（k，v）坐标值，形成速度-波数（v-k）曲线，便实现了瑞雷面波的速度-波数（v-k）分析。

2.根据权利要求1所述的一种瑞雷面波速度-波数分析方法，其特征在于，包括如下步骤：

第一步，将含有100道，每道有1000个采样点，时间采样率为4ms的单炮瑞雷面波地震记录读取到二维数组g中；

第二步，对g进行二维傅里叶变换，得到f-k谱，二维离散傅里叶正变换为：

       

  

其中，k为波数，为波数采样率，j₀=1，2，…，100，为波数采样序号；f为频率，为频率采样率，l₀=1，2，…，1000，为频率采样序号；j=1，2，…，100，为道序号；l=1，2，…，1000，为采样点序号；为道间距，为最小偏移距，x为偏移距；为时间采样率，t为采样时间； 为虚数单位；

第三步，建立传播速度（v）与频率（f）、波数（k）之间的关系，具体如下：

                                                

第四步，依据瑞雷面波的传播速度（v）与频率（f）、波数（k）之间的关系，将f-k谱中对应（k，f）坐标的能量映射到v-k谱中对应（k，v）坐标的位置，同时选定v-k谱输出的坐标范围并成图，其中传播速度（v）的范围为：200～800m/s，波数（k）范围为0～0.05/m^-1；

第五步，基于v-k谱，从左至右依次拾取v-k谱中的强能量团，并提取对应的（k，v）坐标值，形成速度-波数（v-k）曲线，便实现了瑞雷面波的速度-波数（v-k）分析。