地震信号动态滤波方法及系统技术方案

公开了一种地震信号动态滤波方法及系统。该方法可以包括：对地震信号进行时频谱分析，获取地震信号的时频谱；根据地震信号的时频谱，获得每个时刻的瞬时谱；分别对每个时刻的瞬时谱进行小波分频滤波，获取对应的小波分频滤波瞬时谱；对小波分频滤波后的瞬时谱进行时频反变换，获得滤波后的时间域地震信号。本发明专利技术充分考虑了地震信号的时变特征，对各时刻的瞬时谱做滤波处理，滤波结果能够充分的体现地震信号的每一时刻的主要成分或者指定的某个特有的成分。

Dynamic filtering method and system of seismic signal

A dynamic filtering method and system for seismic signals are disclosed. The method can include: time-frequency spectrum analysis of seismic signal to obtain the time-frequency spectrum of seismic signal; instantaneous spectrum of each time is obtained according to the time-frequency spectrum of seismic signal; instantaneous spectrum of each time is filtered by wavelet frequency-division filtering to obtain the corresponding instantaneous spectrum of wavelet frequency-division filtering; instantaneous spectrum after wavelet frequency-division filtering is inversely transformed to obtain the filtered instantaneous spectrum. Inter-domain seismic signals. The invention fully considers the time-varying characteristics of seismic signals and filters the instantaneous spectrum of each time. The filtering result can fully reflect the main component of seismic signals at each time or a specified specific component.

【技术实现步骤摘要】
地震信号动态滤波方法及系统
本专利技术涉及地球物理勘探领域，更具体地，涉及一种地震信号动态滤波方法及系统。
技术介绍
在地震信号处理领域中，经常需要对地震信号进行滤波或者分频处理，以突出、获得或者去除某个频率或频段的信号。频率域滤波是最常用的滤波方法，比如带通滤波、低通滤波、高通滤波等。在地震信号处理领域中，褶积也是一种常用的滤波方法，它是在时间域进行的，滤波效果类似于带通滤波，其滤波器的形态由子波振幅谱决定。但是，传统的这种滤波方法在对地震信号进行滤波处理时，一旦滤波器确定之后，对每道地震数据，甚至每个时刻的信号都是采用同一套滤波参数进行处理。而地震数据是一种时变信号，浅层信号的主频较高、频带较宽，而深层信号的主频较低、频带较窄，甚至不同时刻信号的频率分布范围都不同。因此，有必要开发一种地震信号动态滤波方法及系统。
技术实现思路
本专利技术提出了一种地震信号动态滤波方法及系统，充分考虑了地震信号的时变特征，对各时刻的瞬时谱做滤波处理，滤波结果能够充分的体现地震信号的每一时刻的主要成分或者指定的某个特有的成分。根据本专利技术的一方面，提出了一种地震信号动态滤波方法。所述方法可以包括：对地震信号进行时频分析，获取所述地震信号的时频谱；根据所述地震信号的时频谱，获得每个时刻的瞬时谱；分别对每个时刻的瞬时谱进行小波分频滤波，获取对应的小波分频滤波瞬时谱；对所述小波分频滤波后的瞬时谱进行时频反变换，获得滤波后的时间域地震信号。优选地，对每个时刻的瞬时谱进行小波分频滤波，获取小波分频滤波瞬时谱包括：针对所述时频谱每个时刻的瞬时谱，拾取所述瞬时谱能量最大处对应的频率为该时刻的主频，确定小波的主频与时间域形态；通过傅里叶变换计算所述小波的振幅谱，对所述振幅谱进行归一化处理，获得小波分频滤波的带通滤波器；用所述带通滤波器乘以所述瞬时谱，获取所述小波分频滤波瞬时谱。优选地，还包括：以所述瞬时谱能量最大处对应的频率为所述小波的主频。优选地，还包括：设定频率阈值，将距离所述瞬时谱能量最大处对应的频率小于或等于所述频率阈值的范围作为频率范围；将所述频率范围内的任意频率作为所述小波的主频。优选地，所述时频分析为广义S变换，所述时频反变换为广义S变换的逆变换。根据本专利技术的一方面，提出了一种地震信号动态滤波系统，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：对地震信号进行时频分析，获取所述地震信号的时频谱；根据所述地震信号的时频谱，获得每个时刻的瞬时谱；分别对每个时刻的瞬时谱进行小波分频滤波，获取对应的小波分频滤波瞬时谱；对所述小波分频滤波后的瞬时谱进行时频反变换，获得滤波后的时间域地震信号。优选地，对每个时刻的瞬时谱进行小波分频滤波，获取小波分频滤波瞬时谱包括：针对所述时频谱每个时刻的瞬时谱，拾取所述瞬时谱能量最大处对应的频率为该时刻的主频，确定小波的主频与时间域形态；通过傅里叶变换计算所述小波的振幅谱，对所述振幅谱进行归一化处理，获得小波分频滤波的带通滤波器；用所述带通滤波器乘以所述瞬时谱，获取所述小波分频滤波瞬时谱。优选地，还包括：以所述瞬时谱能量最大处对应的频率为所述小波的主频。优选地，还包括：设定频率阈值，将距离所述瞬时谱能量最大处对应的频率小于或等于所述频率阈值的范围作为频率范围；将所述频率范围内的任意频率作为所述小波的主频。优选地，所述时频分析为广义S变换，所述时频反变换为广义S变换的逆变换。本专利技术的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本专利技术的特定原理。附图说明通过结合附图对本专利技术示例性实施例进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施例中，相同的参考标号通常代表相同部件。图1示出了根据本专利技术的地震信号动态滤波方法的步骤的流程图。图2a和图2b分别示出了根据本专利技术的一个实施例的地震信号及其对应的频谱的示意图。图3示出了根据本专利技术的一个实施例的时频谱的示意图。图4a、图4b和图4c分别示出了根据本专利技术的一个实施例的三个不同时刻的瞬时谱、采用传统滤波方法的滤波结果和采用本方法的滤波结果的示意图。图5a、图5b、图5c、图5d和图5e分别示出了根据本专利技术的一个实施例的地震信号、采用24Hz作为滤波参数对地震信号滤波之后的地震信号、滤波参数为每个瞬时谱的主频得到的地震信号、滤波参数为每个时刻瞬时谱的主频-5Hz得到的地震信号、滤波的参数为每个时刻瞬时谱的主频+5Hz得到的地震信号的示意图。具体实施方式下面将参照附图更详细地描述本专利技术。虽然附图中显示了本专利技术的优选实施例，然而应该理解，可以以各种形式实现本专利技术而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本专利技术更加透彻和完整，并且能够将本专利技术的范围完整地传达给本领域的技术人员。图1示出了根据本专利技术的地震信号动态滤波方法的步骤的流程图。在该实施例中，根据本专利技术的地震信号动态滤波方法可以包括：步骤101，对地震信号进行时频分析，获取地震信号的时频谱；步骤102，根据地震信号的时频谱，获得每个时刻的瞬时谱；步骤103，分别对每个时刻的瞬时谱进行小波分频滤波，获取对应的小波分频滤波瞬时谱；以及步骤104，对小波分频滤波后的瞬时谱进行时频反变换，获得滤波后的时间域地震信号。在一个示例中，对每个时刻的瞬时谱进行小波分频滤波，获取小波分频滤波瞬时谱包括：针对时频谱每个时刻的瞬时谱，拾取瞬时谱能量最大处对应的频率为该时刻的主频，确定小波的主频与时间域形态；通过傅里叶变换计算小波的振幅谱，对振幅谱进行归一化处理，获得小波分频滤波的带通滤波器；用带通滤波器乘以瞬时谱，获取小波分频滤波瞬时谱。在一个示例中，还包括：以瞬时谱能量最大处对应的频率为小波的主频。在一个示例中，还包括：设定频率阈值，将距离瞬时谱能量最大处对应的频率小于或等于频率阈值的范围作为频率范围；将频率范围内的任意频率作为小波的主频。在一个示例中，时频分析为广义S变换，时频反变换为广义S变换的逆变换。具体地，对地震信号进行时频分析，获取地震信号的时频谱，时频谱不限于采用任何一种时频分析方法来计算，例如采用传统的短时傅里叶变换、小波变换、广义S变换等均可；确定滤波方式，不需要给定具体的滤波参数，以各时刻的主频为基准，只需要确定滤波参数为小于主频或者大于主频的一个频率间隔，或者能量从主频处的最大能量衰减到最大能量的某个百分比即可；根据地震信号的时频谱，获得每个时刻的瞬时谱，例如，以每个时刻的瞬时谱的主频基准，输出主频处、小于主频5Hz和大于主频5Hz的滤波结果，比如某一时刻瞬时谱的主频为30Hz，那么主频处的输出结果即为采用30Hz作为滤波参数对该时刻的瞬时谱进行滤波；小于主频5Hz的滤波结果即为采用25Hz作为滤波参数对该时刻的瞬时谱进行滤波；同理，大于主频5Hz的滤波结果即为采用35Hz作为滤波参数对该时刻的瞬时谱进行滤波。每一时刻的主频都是不同的，所以确定一个频率间隔就可以获得每个时刻相关的滤波结果。针对时频谱每个时刻的瞬时谱，拾取瞬时谱能量最大处对应本文档来自技高网...

【技术保护点】
1.一种地震信号动态滤波方法，包括：对地震信号进行时频分析，获取所述地震信号的时频谱；根据所述地震信号的时频谱，获得每个时刻的瞬时谱；分别对每个时刻的瞬时谱进行小波分频滤波，获取对应的小波分频滤波瞬时谱；对所述小波分频滤波后的瞬时谱进行时频反变换，获得滤波后的时间域地震信号。

【技术特征摘要】
1.一种地震信号动态滤波方法，包括：对地震信号进行时频分析，获取所述地震信号的时频谱；根据所述地震信号的时频谱，获得每个时刻的瞬时谱；分别对每个时刻的瞬时谱进行小波分频滤波，获取对应的小波分频滤波瞬时谱；对所述小波分频滤波后的瞬时谱进行时频反变换，获得滤波后的时间域地震信号。2.根据权利要求1所述的地震信号动态滤波方法，其中，对每个时刻的瞬时谱进行小波分频滤波，获取小波分频滤波瞬时谱包括：针对所述时频谱每个时刻的瞬时谱，拾取所述瞬时谱能量最大处对应的频率为该时刻的主频，确定小波的主频与时间域形态；通过傅里叶变换计算所述小波的振幅谱，对所述振幅谱进行归一化处理，获得小波分频滤波的带通滤波器；用所述带通滤波器乘以所述瞬时谱，获取所述小波分频滤波瞬时谱。3.根据权利要求2所述的地震信号动态滤波方法，其中，还包括：以所述瞬时谱能量最大处对应的频率为所述小波的主频。4.根据权利要求2所述的地震信号动态滤波方法，其中，还包括：设定频率阈值，将距离所述瞬时谱能量最大处对应的频率小于或等于所述频率阈值的范围作为频率范围；将所述频率范围内的任意频率作为所述小波的主频。5.根据权利要求1所述的地震信号动态滤波方法，其中，所述时频分析为广义S变换，所述时频反变换为广义S变换的逆变换。6.一种地震信号动态滤波系统，其特征在于，该系统包括：...

【专利技术属性】
技术研发人员：陈科，王鹏燕，徐雷鸣，白俊雨，张杨，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11