发明名称 可控震源地震数据的强脉冲噪声检测及压制方法
摘要 本发明涉及一种可控震源地震数据的强脉冲噪声检测及压制方法。依据数理统计的随机误差理论,计算可控震源地震记录的各采集道基准振幅,输入一道原始非相关数据,对所有峰值及谷值数据分别排序,去除峰值和谷值中的最大值和最小值,定义强脉冲噪声检测的正阈值,在地震数据中检测强脉冲信号,对强脉冲噪声位置的地震数据清零或噪声振幅压缩。经试验,本发明公开的可控震源地震数据的强脉冲噪声检测及压制方法广泛适用于具有较强背景噪声,难以得到高信噪比数据的可控震源地震数据预处理过程,图1和图2展示了采用强脉冲噪声压制技术前后的地震数据,可以看出本发明对改善地震数据质量改善效果明显,能有效恢复被强脉冲噪声淹没的地震波信号。
申请公布号 CN102879825A 申请公布日期 2013.01.16
申请号 CN201210346894.0 申请日期 2012.09.18
申请人 吉林大学 发明人 姜弢;贾海青;徐学纯;林君;葛利华;杨志超;宋健;宋健
分类号 G01V1/36(2006.01)I 主分类号 G01V1/36(2006.01)I
代理机构 长春吉大专利代理有限责任公司 22201 代理人 王立文
主权项 1.一种可控震源地震数据的强脉冲噪声检测及压制方法,其特征在于,包括以下步骤:a、输入一道原始非相关数据,自动拾取该道信号的全部峰值和谷值数据;b、对所有峰值及谷值数据分别排序,去除这些峰值和谷值中的若干个最大值和最小值,求得峰值平均值和谷值平均值,作为该道数据的峰值基准振幅<img file="FDA00002156981500011.GIF" wi="58" he="70" />和谷值基准振幅<img file="FDA00002156981500012.GIF" wi="77" he="63" />具体计算公式如下:<maths num="0001"><![CDATA[<math><mrow><mover><msub><mi>A</mi><mi>p</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><msub><mi>n</mi><mn>1</mn></msub></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mn>1</mn></msub></munderover><msub><mi>p</mi><mi>i</mi></msub><mo>.</mo><mo>.</mo><mo>.</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mover><msub><mi>A</mi><mi>v</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><msub><mi>n</mi><mn>2</mn></msub></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mn>2</mn></msub></munderover><msub><mi>v</mi><mi>i</mi></msub><mo>.</mo><mo>.</mo><mo>.</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中:n<sub>1</sub>为检测到的峰值个数,n<sub>2</sub>为检测到的谷值个数,i为峰值点序号,j为谷值点序号,存在i∈[1,n<sub>1</sub>],j∈[1,n<sub>2</sub>],p<sub>i</sub>为峰值数据,v<sub>j</sub>为谷值数据;c、定义强脉冲噪声检测的正阈值<img file="FDA00002156981500015.GIF" wi="238" he="70" />和负阈值<img file="FDA00002156981500016.GIF" wi="254" he="63" />其中N为权值,原则上N值应大于等于3,一般建议N值取为3至9之间;d、以阈值A<sub>+</sub>和A<sub>-</sub>为门限,采用比较算法,在地震数据中检测强脉冲信号,具体地,令u为采样点序列,每道地震记录用s(u)表示,则所有满足判定条件s(u)≥A<sub>+</sub>…………………………………(3)或s(u)≤A<sub>-</sub>…………………………………(4)的采样点处均存在强噪声;e、对强脉冲噪声位置的地震数据清零或噪声振幅压缩,然后用互相关检测获得脉冲形式的可控震源地震数据。
地址 130012 吉林省长春市前进大街2699号