一种高保真声波干扰压制方法
【专利摘要】本发明公开的高保真声波干扰压制方法包括以下步骤：在共炮点道集的每道地震数据中划分出含有声波干扰的时窗；将所述时窗内的数据进行小波分解并等分为多个频带；求取每道地震数据的每个频带的振幅加权中值；将每道地震数据的每个频带的每个样点的噪声分贝值分别与预设阈值进行比较和处理得到每道地震数据的多个频带的衰减后数据，利用样点的振幅值与样点所在频带的振幅加权中值的比值计算样点的噪声分贝值，当样点的噪声分贝值大于或等于预设阈值时则取所述样点之前的3～5个样点的振幅的平均值代替所述样点的振幅值；将所述每道地震数据的多个频带的衰减后数据进行小波重构得到声波干扰压制后的地震数据。
【专利说明】一种高保真声波干扰压制方法
【技术领域】
[0001]本发明属于地震勘探【技术领域】，更具体地讲，涉及一种用于石油地震资料的高保真声波干扰压制方法。
【背景技术】
[0002]在山地地区，由于激发和接收条件的限制，大量地采用浅井组合及坑炮激发，因而导致声波较为发育。声波会在地震记录数据上形成尖锐的强初至，呈窄带出现却比较稳定。声波的特点是能量强、集中在中高频率段且传播速度稳定，因此声波的分布范围极易划定。如果这种强能量的噪声不能得到有效地压制，将直接影响后续的反褶积、速度分析等处理，对提高地震资料的信噪比、分辨率以及可信度是极为不利的。
[0003]目前用于压制声波干扰的手段仍然仅局限于一维滤波、二维滤波、人工编辑等，例如:(1)反褶积方法，其主要通过反褶积技术对声波进行压制，但在声波主频较高时，此压制方法往往不是很理想；(2)切除法:使用内切除法将声波完全剔除，其虽能从根本上消除声波对地震数据的影响并可更好第提高信噪比，但湮没在强噪声干扰中的有效信号也会损失掉，同时还会大大增加工作量，延长处理周期。
[0004]以上压制声波干扰的方法尽管有一些效果，但也极容易伤害和损失有效反射信息。因此，寻找一种有效的衰减声波方法，无论在提高处理质量还是在缩短处理周期，减轻劳动强度方面都是非常具有实际意义的。
【发明内容】
[0005]针对现有技术中存在的不足，本发明的目的在于解决上述技术问题中的一个或多个。
[0006]本发明基于声波具有非常稳定的传播速度以及能量强的特点，提供了一种数据量小、计算效率高且压制声波干扰效果明显的高保真声波干扰压制方法。
[0007]为了实现上述目的，本发明提供了一种高保真声波干扰压制方法，其特征在于，所述方法包括以下步骤:A)在共炮点道集的每道地震数据中划分出含有声波干扰的时窗；B)将每道地震数据中含有声波干扰的时窗内的数据进行小波分解并等分为多个频带；C)统计每道地震数据的每个频带的所有样点的振幅值并求取每道地震数据的每个频带的振幅加权中值；D)将每道地震数据的每个频带的每个样点的噪声分贝值分别与预设阈值进行比较和处理得到每道地震数据的多个频带的衰减后数据，其中，利用所述样点的振幅值与所述样点所在频带的振幅加权中值的比值计算所述样点的噪声分贝值，当样点的噪声分贝值大于或等于预设阈值时则取所述样点之前的3?5个样点的振幅的平均值代替所述样点的振幅值，当样点的噪声分贝值小于预设阈值时则保留所述样点的振幅值；E)将所述每道地震数据的多个频带的衰减后数据进行小波重构得到声波干扰压制后的地震数据。
[0008]根据本发明的高保真声波干扰压制方法的一个实施例，在步骤A中，包括以下子步骤:A1)在所述共炮点道集中拾取声波干扰起始点的起始时间h以及起始偏移距Γ(ι，再利用式I计算第i道地震数据的声波干扰时间ti:
[0009]ti = t0+abs (r1-r0) /v-式I,式中，V为声波的传播速度；
[0010]A2)给定时移量d，划分第i道地震数据的含有声波干扰的时窗为(t1-cbti+d)，其中，i为地震数据的道号且N，N为所述共炮点道集的总道数A+d ( tmax，tmax为所述共炮点道集的最大记录时间。
[0011]根据本发明的高保真声波干扰压制方法的一个实施例，在步骤C中，包括以下子步骤:C1)将第i道地震数据的fk频带的P个样点的振幅值Λ', 按照大小进行升序排序，为所述P个样点的振幅值分别设置加权系数Ws并通过加权设置得到振幅值序列
【权利要求】
1.一种高保真声波干扰压制方法，其特征在于，所述方法包括以下步骤: A)在共炮点道集的每道地震数据中划分出含有声波干扰的时窗； B)将每道地震数据中含有声波干扰的时窗内的数据进行小波分解并等分为多个频带; C)统计每道地震数据的每个频带的所有样点的振幅值并求取每道地震数据的每个频带的振幅加权中值； D)将每道地震数据的每个频带的每个样点的噪声分贝值分别与预设阈值进行比较和处理得到每道地震数据的多个频带的衰减后数据，其中，利用所述样点的振幅值与所述样点所在频带的振幅加权中值的比值计算所述样点的噪声分贝值，当样点的噪声分贝值大于或等于预设阈值时则取所述样点之前的3?5个样点的振幅的平均值代替所述样点的振幅值，当样点的噪声分贝值小于预设阈值时则保留所述样点的振幅值； E)将所述每道地震数据的多个频带的衰减后数据进行小波重构得到声波干扰压制后的地震数据。
2.根据权利要求1所述的高保真声波干扰压制方法，其特征在于，在步骤A中，包括以下子步骤: Al)在所述共炮点道集中拾取声波干扰起始点的起始时间h以及起始偏移距IV再利用式I计算第i道地震数据的声波干扰时间\:
ti = t0+abs (rj-rg) /v-式 I, 式中，V为声波的传播速度； A2)给定时移量d，划分第i道地震数据的含有声波干扰的时窗为(t1-d，h+d)， 其中，i为地震数据的道号且I彡i彡N，N为所述共炮点道集的总道数山+d ( tmax，tmax为所述共炮点道集的最大记录时间。
3.根据权利要求1所述的高保真声波干扰压制方法，其特征在于，在步骤C中，包括以下子步骤: Cl)将第i道地震数据的fk频带的P个样点的振幅值按照大小进行升序排序，为所述P个样点的振幅值分别设置加权系数Ws并通过加权设置得到振幅值序列 WjW2Wpi^TiU)，(“2):(,.，2V.%， 其中，i为地震数据的道号且I < i < N，N为所述共炮点道集的总道数；k为频带的序号且I彡k彡M，M为第i道地震数据的频带总数；p为样点总数且P为正整数，s为样点编号且I彡s彡P ; Xh (i，S)为第i道地震数据的fk频带的第s个样点的振幅值；加权系数Ws为正整数并且表示对应振幅值的存在个数，W1 < W2 <...< wp/2且wp/2 > wp/2+1 > Wp ; C2)对所述振幅值序列的多个振幅值按照大小进行升序排序，当?'为奇数时，
J=I则第i道地震数据的fk频带的振幅加权中值ΜΛ(/)等于升序排序后的振幅值序列的第(Ywt +1) /2个振幅值；当;I；Wi为偶数时，则第i道地震数据的fk频带的振幅加权中值ir~ iS=lμ _等于升序排序后的振幅值序列的第
个振幅值的平均值； C3)重复子步骤Cl至C2，直至得到每道地震数据的每个频带的振幅加权中值。
4.根据权利要求1所述的高保真声波干扰压制方法，其特征在于，在步骤D中，根据式2计算样点s的噪声分贝值dB:
式中，(i，s)为第i道地震数据的fk频带的第S个样点的振幅值，为第i道地震数据的fk频带的振幅加权中值。
【文档编号】G01V1/36GK104133248SQ201410389621
【公开日】2014年11月5日 申请日期:2014年8月8日 优先权日:2014年8月8日
【发明者】曹中林, 张晓斌, 何光明, 陈爱萍, 李振, 陈丹, 杨海涛, 王珑, 张茜 申请人:中国石油集团川庆钻探工程有限公司地球物理勘探公司