用于三分量地震背景噪声的时频双域P波提取及映射成像方法与流程

本发明涉及地震勘探，尤其涉及一种用于三分量地震背景噪声的时频双域p波提取及映射成像方法。

背景技术：

1、地震背景噪声成像技术是二十一世纪以来地震学领域最重要的突破性进展，该技术利用连续地震记录中普遍存在的“背景噪声”，通过单台站内或不同台站间连续信号的自(互)相关提取经验格林函数，并用以提取地质构造和属性信息。相比于传统地震探测技术，背景噪声成像技术从随机场源微扰信号中提取地下介质信息，具备成像分辨率高、观测时效好、重复性好、能够实现无损勘察等优势。近年来，随着该方法的不断深入发展，地震背景噪声成像逐渐从地壳、地幔与岩石圈等地球深部结构探测延伸到固体矿集区和城市地下空间全要素成像等浅地表勘察领域，并取得了日渐广泛的应用。

2、在众多地震背景噪声成像技术当中，以利用面波频散效应进行横波速度反演为主，原因在于地表附近震源(人类活动)所激发的能量多为面波能量；相对而言，地震背景噪声中的体波能量(主要指p波)属于弱信号，往往不予充分利用。将体波信号从地震背景噪声中分离出来有益于后续开展背景噪声反射波地震勘探研究，并提供更加丰富的地下参数信息及成像结果。现有分离技术中，通常基于面波和体波信号的频谱差异进行波型分离，然而由于频谱混叠效应，这类技术的分离效果不佳。

技术实现思路

1、为了解决现有技术存在的问题，本发明提供一种分离效果好、能优化地质解译效果、保证设计和施工安全的用于三分量地震背景噪声的时频双域p波提取及映射成像方法。

2、为此，本发明采用以下技术方案：

3、一种用于三分量地震背景噪声的时频双域p波提取及映射成像方法，包括以下步骤：

4、s1，对原始三分量地震背景噪声数据进行预处理，并设置时间窗口长度；

5、s2，对三分量地震背景噪声数据进行时间域极化滤波，分离纵波和横波，包括以下分步骤；

6、s21，计算时间窗口内三分量x、y、z地震背景噪声数据的协方差矩阵c：

7、

8、其中：

9、xi、yi、zi分别为x、y、z三分量数据的第i个数；

10、分别为x、y、z三分量数据的均值；

11、k为时间窗口内采样点数；

12、s22，对步骤s21得到的协方差矩阵进行奇异值分解，得到协方差矩阵的奇异值和奇异值向量(v1,v2,v3)；

13、s23，根据所述奇异值向量计算最大极化方向的倾角angle，

14、

15、其中，v1,x、v1,y、v1,z分别为奇异值向量v1中的三个分量；

16、s24，根据最大极化方向倾角进行p波、横波判别：

17、如果某一个台站接收到的倾角angle接近0度，则该时间窗口所对应的波形为p波；如果某一个台站接收到的倾角接近90度，则该时间窗口所对应的波形为横波；

18、s3，对上述时间域极化滤波后的数据进行频域数据筛选，进一步消除背景噪声数据中的低频面波信号，突显p波能量，包括：

19、s31，对数据进行带通滤波，保存有效频带范围；

20、s32，对噪声数据进行自相关计算，并自动增益，进行能量均衡；

21、s33，将增益后的自相关记录转换到频率域，对其功率谱进行统计分析，计算方差s2，得到对应的功率谱，方差s2的计算公式如下：

22、

23、其中，f(f)为自相关记录的功率谱；n为功率谱中的频率点数；为功率谱均值；

24、s34，根据功率谱的统计性规律对频域数据进行筛选：根据体波功率谱比较平滑、面波功率谱比较震荡的性质，认定方差s2大于一定阈值的信号为面波，其余信号为p波；

25、s4，输出处理后的背景噪声数据；

26、s5，对步骤s4输出的处理后的背景噪声数据进行映射成像计算及映射剖面输出，得到勘探区域的地震映射结果。

27、上述的步骤s5中，映射成像的计算公式为：

28、

29、其中index为x,y或z时分别表示地震背景噪声的x、y或z分量，iindex(t)为映射成像结果，m为地震背景噪声记录分割的时间窗口个数，rindexrindex(t)表示相应分量某一时间窗口内信号的自相关函数。

30、优选的是，步骤s31中所述有效频带范围为3hz-40hz；步骤s34中所述一定阈值为0.7。

31、步骤s1所述的预处理包括：

32、s11，进行趋势项去除操作；

33、s12，进行第一次归一化：根据设置的阈值对绝对值大于阈值的数据进行置零，消除信号中的野值和均衡不同地震道之间的能量；

34、s13，进行第二次归一化处理；

35、s14，设置时间窗口长度。优选的时间窗口长度为500ms。

36、本发明提供了一种新的应用于三分量地震背景噪声的时-频双域体波提取技术，基于此，进一步提供了地震映射成像方法。通过获取勘探区域中设置的每个地震仪采集到的三分量背景噪声信号，并在时-频双域分别进行极化滤波和数据筛选，弱化了地震背景噪声中的面波能量,增强了体波能量。与现有技术相比，本发明具有以下有益效果：

37、(1)本发明的用于地震背景噪声的时-频双域p波提取方法充分利用了面波、体波在时间域的偏振极化差异以及在频率分布上的不同，最大限度地弱化了地震背景噪声中的面波、突显了体波能量，分离效果极佳；

38、(2)本发明的基于地震背景噪声的映射成像方法，额外为被动源勘察提供了一种成像结果，进一步优化了地质解译效果，保证了设计、施工安全。

39、(3)相比于传统人工源映射成像技术，本发明方法具有绿色无损、易实施、成本低等优势。

技术特征：

1.一种用于三分量地震背景噪声的时频双域p波提取及映射成像方法，包括以下步骤：

2.根据权利要求1所述的用于三分量地震背景噪声的时频双域p波提取及映射成像方法，其特征在于，步骤s5中映射成像的计算公式为：

3.根据权利要求1所述的用于三分量地震背景噪声的时频双域p波提取及映射成像方法，其特征在于：步骤s31中所述有效频带范围为3hz-40hz。

4.根据权利要求1所述的用于三分量地震背景噪声的时频双域p波提取及映射成像方法，其特征在于：步骤s34中所述一定阈值为0.7。

5.根据权利要求1所述的用于三分量地震背景噪声的时频双域p波提取及映射成像方法，其特征在于，步骤s1所述的预处理包括：

6.根据权利要求5所述的用于三分量地震背景噪声的时频双域p波提取及映射成像方法，其特征在于：步骤s14所述时间窗口长度为500ms。

技术总结
本发明公开了一种用于三分量地震背景噪声的时频双域P波提取及映射成像方法，包括：S1，对原始三分量地震背景噪声数据进行预处理，并设置时间窗口长度；S2，对三分量地震背景噪声数据进行时间域极化滤波，分离纵波和横波；S3，对上述时间域极化滤波后的数据进行频域数据筛选，进一步消除背景噪声数据中的低频面波信号，突显P波能量；S4，输出处理后的背景噪声数据；S5，对步骤S4输出的处理后的背景噪声数据进行映射成像计算及映射剖面输出，得到勘探区域的地震映射结果。该方法最大限度地弱化了地震背景噪声中的面波、突显了体波能量，其分离效果极佳，进一步优化了地质解译效果，具有绿色无损、易实施、成本低等优势。

技术研发人员：王官超,牛永效,祁晓雨,秦海旭,李国和,齐春雨,李志华,刘占峰,赵广茂,许广春,陈承申
受保护的技术使用者：中国铁路设计集团有限公司
技术研发日：
技术公布日：2024/1/13