【技术实现步骤摘要】
本专利技术涉及地球探测与信息，尤其是涉及基于多重同步压缩变换的面波群速度频散能量成像方法。
技术介绍
1、横波速度是近地表岩土工程十分重要的参数，其准确与否对于水文地质与工程地质十分重要。地震面波是无损获取近地表横波速度的重要方法之一。当前地震面波方法主要是利用瑞雷波的相速度反演计算而得，多数研究也集中在相速度频散能量成像方面，而对于群速度的频散能量成像方面开展的研究较少。近年来，国内外一些专家对群速度频散能量开展了研究，主要方法有多重滤波法、小波变换法及基于s变换的时频分析法。但这些方法多存在频散能量图分辨率较低、鲁棒性较差、计算量较大等不足。
2、当前面波群速度频散能量成像方法有多重滤波法、小波变换法及基于s变换的时频分析法。但这些方法多存在频散能量图分辨率较低、鲁棒性较差、计算量较大等不足。如当前较为常用的多重滤波方法对地震记录采用一系列带通滤波器，并计算所得信号的包络来估计群速度。该方法利用的是带通滤波器的中心频率或包络的瞬时频率作为群速度计算时的频率。此方法存在能量聚焦性较差且在低频端群速度计算误差较大的不足。基于s变换的时频分析方法通过调节高斯时窗宽度去获得较高的时频分辨率并提高低频段频散曲线估计的准确度，但在实践中发现该方法在面波多模态时频散成像效果较差、分辨率不高且调节时窗参数十分困难。
技术实现思路
1、本专利技术的目的是提供基于多重同步压缩变换的面波群速度频散能量成像方法，相对于其他方法主要在两个环节进行了改进：（1）从时频分析环节通过采用时频谱大幅
2、为实现上述目的，本专利技术提供了基于多重同步压缩变换的面波群速度频散能量成像方法，包括以下步骤：
3、s1、输入地震面波数据 s( t,x)并对该地震面波数据的每一道进行多重同步压缩变换到时间-频率域，将每道的时间-频率域数据合并到一起产生一个三维时频数据体 ts( t, x, f)；
4、s2、对步骤s1的三维时频数据体 ts( t, x, f)沿频率轴提取频率为 f i时的共频率切片数据 ts( t, x)；
5、s3、设置群速度区间[ v min ,v max]，将在步骤s2得到的共频率切片数据 ts( t, x)采用相似系数法倾斜扫描叠加的方式计算面波群速度区间内不同速度的叠加能量；
6、s4、对频率范围为[ f min ,f max]的共频率切片数据重复步骤s2与步骤s3，将不同频率、不同面波群速度下的叠加能量合并在一起得到最终的面波频散能量图。
7、优选的，在步骤s1中，包括以下步骤：
8、s11、首先对地震面波数据进行振幅增益、滤波及切除处理，增强面波信号，然后对每道地震面波数据进行短时傅里叶变换；
9、s12、设置迭代次数为 n，对短时傅里叶变换结果进行多重同步压缩变换。
10、优选的，在步骤s11中，其中，单道地震信号 s( t)表示为：
11、 （1）
12、式中， k为单组分信号的个数， t为时间， x为偏移距， a k( t)为瞬时振幅， ᵠ k( t)为瞬时相位， i为虚数单位；
13、对应的地震信号短时傅里叶变换为：
14、 （2）
15、式中，为信号 s( t)的短时傅里叶变换， g( u-t)为时窗， u为时移因子，是角频率。
16、优选的，在步骤s12中，其多重同步压缩变换过程如下表述：
17、 （3）
18、式中，为脉冲函数，是多次迭代估算的瞬时频率， η为同步压缩变换的频率。
19、优选的，在步骤s3中，对步骤s2的所有频率范围的共频率切片数据 ts( t, x)在截距时间 t=0处采用相似系数法倾斜扫描叠加的方式计算不同面波群速度时的叠加能量，其中，相似系数法倾斜扫描叠加的公式如下：
20、 （4）
21、式中， sc为不同群速度计算的叠加能量， l为地震道数目 ，m为截距时间 t=0~ m的时窗， r i为偏移距为 x时群速度 v对应的延迟时，即 r i =( x 本文档来自技高网...
【技术保护点】
1.基于多重同步压缩变换的面波群速度频散能量成像方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的基于多重同步压缩变换的面波群速度频散能量成像方法，其特征在于，在步骤S1中，包括以下步骤：
3.根据权利要求2所述的基于多重同步压缩变换的面波群速度频散能量成像方法，其特征在于，在步骤S11中，其中，单道地震信号s(t)表示为：
4.根据权利要求3所述的基于多重同步压缩变换的面波群速度频散能量成像方法，其特征在于，在步骤S12中，其多重同步压缩变换过程如下表述：
5.根据权利要求4所述的基于多重同步压缩变换的面波群速度频散能量成像方法，其特征在于，在步骤S3中，对步骤S2的所有频率范围的共频率切片数据Ts(t,x)在截距时间T=0处采用相似系数法倾斜扫描叠加的方式计算不同面波群速度时的叠加能量，其中，相似系数法倾斜扫描叠加的公式如下：
【技术特征摘要】
1.基于多重同步压缩变换的面波群速度频散能量成像方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的基于多重同步压缩变换的面波群速度频散能量成像方法，其特征在于，在步骤s1中，包括以下步骤：
3.根据权利要求2所述的基于多重同步压缩变换的面波群速度频散能量成像方法，其特征在于，在步骤s11中，其中，单道地震信号s(t)表示为：
4.根据权利要求3所述的基...
【专利技术属性】
技术研发人员：何滔，彭苏萍，崔晓芹，耿恒高，
申请(专利权)人：中国矿业大学北京，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。