压制地震资料中高频环境噪声的方法技术

提供一种压制地震资料中高频环境噪声的方法，包括：(a)根据叠前单炮地震数据在初至时间前的能量分布确定干扰道和正常道；(b)将干扰道及其左右邻近的正常道组成待处理地震数据，对每道地震数据进行傅里叶正变换，确定频率域的每道地震数据的振幅谱，将振幅谱划分为多个频段，利用自适应门槛值法判断干扰道的各频段是否存在高频环境噪声；(c)当干扰道在特定频段具有高频环境噪声时，线性内插出特定频段的振幅值，并替换特定频段的原振幅值；(d)将经线性内插处理后的干扰道地震数据进行傅里叶反变换，得到压制高频环境噪声后的干扰道地震数据。根据所述方法，能够有效提高地震资料的信噪比，为地震资料的后续处理奠定了有利基础。

【技术实现步骤摘要】

本专利技术总体说来涉及地震资料处理领域，更具体地讲，涉及一种压制地震资料中高频环境噪声的方法。
技术介绍
地震勘探在石油物探中占有重要的地位，但是随着地震勘探环境的复杂性日益增加，地震资料中往往存在许多噪声，噪声对有效波具有较强的干扰，导致地震资料的信噪比低，严重影响地震资料的后续处理以及地震剖面解释的准确性。通常，地震资料中的噪声主要为高频环境噪声，高频环境噪声的特征是频带范围较宽,能量相对集中在较高频段,但其总能量相对于地震波则较弱,而且随时间的变化显得平稳。导致这类高频环境噪声的主要原因是来自接收仪器的热噪声和自然界的各种相对平稳的高频干扰源。目前，压制高频环境噪声方法主要有:一、利用门槛值把含高频环境噪声的整个地震道剔除，以达到去除高频环境噪声的目的,该方法剔除地震道的多少与门槛值有直接关系,对于覆盖次数较低的地震资料来说,由于减少了覆盖次数,会影响浅层有效波同相轴成像，因此，利用该方法压制高频环境噪后得到的地震资料质量并不理想；二、通过比较含高频噪声振幅采样点时窗的频谱与整道频谱的差异,对高频环境噪声振幅进行编辑,从而减少高频环境噪声，但是，由于通过频谱差异很难准确确定高频环境噪声的位置,因此,这种压制高频环境噪声的方法的效果也不理想。综上所述，现有的压制高频环境噪声的方法并不能满足人们对高质量地震资料的需求。
技术实现思路
本专利技术示例性实施例的目的在于提供一种压制地震资料中高频环境噪声的方法，以克服现有技术中在压制高频环境噪声后得到的地震资料质量不理想的问题。根据本专利技术示例性实施例的一方面，提供一种压制地震资料中高频环境噪声的方法，包括：(a)根据地震资料中叠前单炮地震数据在初至时间前的能量分布确定地震资料中包括高频环境噪声的干扰道和不包括高频环境噪声的正常道；(b)将干扰道及其左右邻近的预定数量个正常道组成待处理地震数据，对所述待处理地震数据的每道地震数据进行傅里叶正变换，得到频率域的待处理地震数据，确定频率域的待处理地震数据中每道地震数据的振幅谱，将所述振幅谱划分为多个频段，利用自适应门槛值法判断干扰道的各频段是否存在高频环境噪声；(c)当干扰道在特定频段具有高频环境噪声时，利用干扰道左右邻近的正常道在所述特定频段的振幅值线性内插出干扰道在所述特定频段的振幅值，并将线性内插出的振幅值替换干扰道在所述特定频段的原振幅值；(d)将经线性内插处理后的干扰道地震数据进行傅里叶反变换，得到压制高频环境噪声后的干扰道地震数据。可选地，步骤(a)可包括：(a1)根据地震资料中叠前单炮地震数据中每道地震数据的偏移距和直达波速度确定每道地震数据的初至时间；(a2)基于地震数据采样时间获取每道地震数据从初始时刻到初至时间的初至地震数据；(a3)确定初至地震数据中每道地震数据的平均振幅值；(a4)比较初至地震数据中每道地震数据的平均振幅值与预定能量门槛值的大小，当所述平均振幅值大于所述预定能量门槛值时，确定所述道地震数据所对应的地震道为干扰道。可选地，步骤(b)可包括：(b1)将干扰道及其与其左右邻近的所述预定数量个正常道组成待处理地震数据，并对所述待处理地震数据的每道地震数据进行傅里叶正变换，得到频率域的待处理地震数据；(b2)获取所述频率域的待处理地震数据中每道地震数据中每个频点的振幅值，并基于所述每道地震数据中每个频点的振幅值获得所述道地震数据的振幅谱；(b3)将所述振幅谱按预定频率间隔分为至少一个频段，确定每道地震数据中每个频段的振幅平均值；(b4)利用各道地震数据在相同频段的振幅平均值确定所述频段的振幅中值，并利用所述振幅中值和预定门槛因子确定所述频段的自适应门槛值；(b5)比较干扰道的各频段振幅平均值与其相应的自适应门槛值的大小，当振幅平均值大于自适应门槛值时，确定干扰道在所述频段具有高频环境噪声。可选地，在步骤(c)中，可利用干扰道左右邻近的两个正常道在所述特定频段的振幅值线性内插出干扰道在所述特定频段的振幅值。可选地，在步骤(a1)中，可利用下面的等式确定第i道地震数据的初至时间：tbi=offiv,]]>其中，tbi表示第i道地震数据的初至时间，offi表示第i道地震数据的偏移距，v表示第i道地震数据的直达波速度。可选地，在步骤(a3)中，可利用下面的等式确定初至地震数据中第i道地震数据的平均振幅值：Ei=Σj'=0tbi|yi,j'|tbibs,]]>其中，Ei表示初至地震数据中第i道地震数据的平均振幅值，yi,j'表示初至地震数据中第i道第j'时刻的地震数据，tbi表示第i道地震数据的初至时间，ns表示地震数据的采样时间。可选地，在步骤(b2)中，可利用下面的等式确定频率域的地震数据中第i道第l频点的振幅值：Ai',l＝abs(fft(zi',j))其中，Ai',l表示频率域的待处理地震数据中第i道第l频点的振幅值，其中，l表示频点，1≤l≤PM，PM为最大频点数；zi',j表示待处理地震数据中第i'道第j时刻的地震数据，fft表示傅里叶正变换，abs表示取振幅值。可选地，在步骤(b3)中，可利用下面的等式确定频率域的地震数据中第i'道第k频段的振幅平均值：MeanAi',k=Σl=1+(k-1)×npnp+(k-1)×npAi',lnp,]]>其中，MeanAi',k表示频率域的地震数据中第i'道第k频段的振幅平均值，其中，1≤k≤n，n表示频段个数且为正整数；Ai',l表示频率域的待处理地震数据中第i'道第l频点的振幅值，其中，l表示频点，1≤l≤PM，PM为最大频点数；np表示每个频段中的频点个数，其中，可选地，在步骤(b4)中，可利用下面的等式来确定第k频段的自适应门槛值：Trek＝Meank×F,其中，Trek表示第k频段的自适应门槛值，Meank表示第k频段的振幅中值，F表示所述预定门槛因子。可选地，所述频段可仅包括50赫兹以上的频率。在根据本专利技术示例性实施例的压制地震资料中高频环境噪声的方法中，能够有效压制地震资料中的高频环境噪声，提高了地震资料的信噪比，此外，由于只对干扰道高频段的高频环境噪声进行压制，避开了有效波主频段，从而保护了有效波频率成分的空间分布，使反射波同相轴更加清晰、连续性更好，并且，经压制后的地震资料为地震资料的后续处理奠定了有利基础。附图说明通过下面结合附图进行的详细描述，本本文档来自技高网...

【技术保护点】
一种压制地震资料中高频环境噪声的方法，包括：(a)根据地震资料中叠前单炮地震数据在初至时间前的能量分布确定地震资料中包括高频环境噪声的干扰道和不包括高频环境噪声的正常道；(b)将干扰道及其左右邻近的预定数量个正常道组成待处理地震数据，对所述待处理地震数据的每道地震数据进行傅里叶正变换，得到频率域的待处理地震数据，确定频率域的待处理地震数据中每道地震数据的振幅谱，将所述振幅谱划分为多个频段，利用自适应门槛值法判断干扰道的各频段是否存在高频环境噪声；(c)当干扰道在特定频段具有高频环境噪声时，利用干扰道左右邻近的正常道在所述特定频段的振幅值线性内插出干扰道在所述特定频段的振幅值，并将线性内插出的振幅值替换干扰道在所述特定频段的原振幅值；(d)将经线性内插处理后的干扰道地震数据进行傅里叶反变换，得到压制高频环境噪声后的干扰道地震数据。

【技术特征摘要】
1.一种压制地震资料中高频环境噪声的方法，包括：
(a)根据地震资料中叠前单炮地震数据在初至时间前的能量分布确定地
震资料中包括高频环境噪声的干扰道和不包括高频环境噪声的正常道；
(b)将干扰道及其左右邻近的预定数量个正常道组成待处理地震数据，
对所述待处理地震数据的每道地震数据进行傅里叶正变换，得到频率域的待
处理地震数据，确定频率域的待处理地震数据中每道地震数据的振幅谱，将
所述振幅谱划分为多个频段，利用自适应门槛值法判断干扰道的各频段是否
存在高频环境噪声；
(c)当干扰道在特定频段具有高频环境噪声时，利用干扰道左右邻近的
正常道在所述特定频段的振幅值线性内插出干扰道在所述特定频段的振幅
值，并将线性内插出的振幅值替换干扰道在所述特定频段的原振幅值；
(d)将经线性内插处理后的干扰道地震数据进行傅里叶反变换，得到压
制高频环境噪声后的干扰道地震数据。
2.如权利要求1所述的方法，其中，步骤(a)包括：
(a1)根据地震资料中叠前单炮地震数据中每道地震数据的偏移距和直
达波速度确定每道地震数据的初至时间；
(a2)基于地震数据采样时间获取每道地震数据从初始时刻到初至时间
的初至地震数据；
(a3)确定初至地震数据中每道地震数据的平均振幅值；
(a4)比较初至地震数据中每道地震数据的平均振幅值与预定能量门槛
值的大小，当所述平均振幅值大于所述预定能量门槛值时，确定所述道地震
数据所对应的地震道为干扰道。
3.如权利要求1所述的方法，其中，步骤(b)包括：
(b1)将干扰道及其左右邻近的所述预定数量个正常道组成待处理地震
数据，并对所述待处理地震数据的每道地震数据进行傅里叶正变换，得到频
率域的待处理地震数据；
(b2)获取所述频率域的待处理地震数据中每道地震数据中每个频点的
振幅值，并基于所述每道地震数据中每个频点的振幅值获得所述道地震数据

\t的振幅谱；
(b3)将所述振幅谱按预定频率间隔分为至少一个频段，确定每道地震
数据中每个频段的振幅平均值；
(b4)利用各道地震数据在相同频段的振幅平均值确定所述频段的振幅
中值，并利用所述振幅中值和预定门槛因子确定所述频段的自适应门槛值；
(b5)比较干扰道的各频段振幅平均值与其相应的自适应门槛值的大小，
当振幅平均值大于自适应门槛值时，确定干扰道在所述频段具有高频环境噪
声。
4.如权利要求...

【专利技术属性】
技术研发人员：曹中林，何光明，陈爱萍，罗红明，金德刚，杨海涛，王珑，熊晶璇，李振，
申请(专利权)人：中国石油集团川庆钻探工程有限公司地球物理勘探公司，
类型：发明
国别省市：四川;51