[發明專利]變窗口加權地震數據尖峰噪聲壓制方法

有效

• 申請號: 201910824963.6
• 申請日: 2019-09-02
• 公開（公告）號: CN110542927B
• 公開（公告）日: 2020-07-03
• 發明（設計）人: 姜弢;王京椰;王鑫;林君;晁云峰;鄭晶涵;胡秋月;覃鍇;周琪;楊大鵬;鄭凡;孫淑琴
• 申請（專利權）人: 吉林大學
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 長春吉大專利代理有限責任公司 22201
• 代理人: 王立文

1.一種變窗口加權地震數據尖峰噪聲壓制方法，其特征在于，包括以下步驟：

a、讀取一道含多尖峰噪聲干擾的地震信號x(n)(n＝1,2,...,N)，N為地震信號總采樣點數，通過觀察分析含尖峰噪聲的時域波形，確定尖峰噪聲的最短持續時間t_p，信號采樣率為Fs；

b、計算地震信號x(n)的特征函數

利用遞歸法計算地震信號振幅的短時變化量

長時變化量

式中L_l和L_s為整數，且滿足0.01Fs≤L_s≤9Fs·t_p且L_l/L_s＝2；定義

c、定義尖峰噪聲的檢測閾值

d、滿足ψ(n)＞η所有連續采樣點區間E＝{E₁,E₂,...,E_K}，K為超過閾值的連續采樣點區間的個數，第k個連續采樣點區間為和為該連續采樣點區間起始采樣點序號和終止采樣點序號,k∈[1,K]；

e、對每一個E_k判斷其長度是否滿足條件

若滿足條件，該區間則被認為存在尖峰噪聲，否則視為不存在尖峰噪聲；

f、經過步驟e判斷得出尖峰噪聲的連續采樣點區間E′＝{E′₁,E′₂,...,E′_J}，J為尖峰噪聲個數，為第j個尖峰噪聲段，和為該尖峰噪聲所在連續采樣點區間的起始采樣點序號和終止采樣點序號,j∈[1,J]；定義第j個尖峰噪聲的持續采樣點數的一半為

為向上取整運算；

g、對于每一個E′_j的定義兩個權值

式中

h、對于每一個E′_j的采樣點，在其兩側各取L_j窗長的信號

令

式(10)中Υ()為計算中值的函數；x′(n)為壓制尖峰噪聲后的信號。