[發明專利]一種利用Lipschitz投影的微震信號濾波方法和系統

在審

• 申請號: 202010445175.9
• 申請日: 2020-05-23
• 公開（公告）號: CN111679320A
• 公開（公告）日: 2020-09-18
• 發明（設計）人: 翟明岳
• 申請（專利權）人: 廣東石油化工學院
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.所述一種利用Lipschitz投影的微震信號濾波方法，其特征在于，包括：

步驟101獲取按時間順序采集的信號序列S；

步驟102創建逼近矢量x并賦值為x＝S；創建計數變量k并賦值為k＝0；

步驟103創建中間矢量v并賦值為

步驟104求取所述中間矢量v的梯度值▽v，其計算公式為：▽v＝[0,v₂-v₁,v₃-v₂,···,v_N-v_N-1]。其中，v₁,v₂,v₃,···,v_N-1,v_N分別為所述中間矢量v的第1,2,···,N-1,N個元素；

步驟105求取更新步長t_k并更新所述逼近矢量x。所述更新步長t_k的求取公式為：所述逼近矢量x的更新公式為：x＝v-t_k▽v。其中，ΔT為所述信號序列S的采樣間隔；m₀為所述信號序列S的均值；σ₀為所述信號序列S的均方差；

步驟106求取逼近誤差e，所述逼近誤差e求取公式為：e＝|x-S|；

步驟107判斷所述逼近誤差e是否大于或等于預設閾值ε₀，得到第一判斷結果。如果所述第一判斷結果顯示所述逼近誤差e大于或等于所述預設閾值ε₀，則返回所述步驟105、所述步驟106和所述步驟107，所述計數變量k的值加1；直至所述第一判斷結果顯示所述逼近誤差e小于所述預設閾值ε₀；其中所述預設閾值為ε₀＝0.001；

步驟108記錄濾除了噪聲的信號序列S_NEW，具體為S_NEW＝x。

2.所述一種利用Lipschitz投影的微震信號濾波系統，其特征在于，包括：

模塊201獲取按時間順序采集的信號序列S；

模塊202創建逼近矢量x并賦值為x＝S；創建計數變量k并賦值為k＝0；

模塊203創建中間矢量v并賦值為

模塊204求取所述中間矢量v的梯度值▽v，其計算公式為：▽v＝[0,v₂-v₁,v₃-v₂,···,v_N-v_N-1]。其中，v₁,v₂,v₃,···,v_N-1,v_N分別為所述中間矢量v的第1,2,···,N-1,N個元素；

模塊205求取更新步長t_k并更新所述逼近矢量x。所述更新步長t_k的求取公式為：所述逼近矢量x的更新公式為：x＝v-t_k▽v。其中，ΔT為所述信號序列S的采樣間隔；m₀為所述信號序列S的均值；σ₀為所述信號序列S的均方差；

模塊206求取逼近誤差e，所述逼近誤差e求取公式為：e＝|x-S|；

模塊207判斷所述逼近誤差e是否大于或等于預設閾值ε₀，得到第一判斷結果。如果所述第一判斷結果顯示所述逼近誤差e大于或等于所述預設閾值ε₀，則返回所述模塊205、所述模塊206和所述模塊207，所述計數變量k的值加1；直至所述第一判斷結果顯示所述逼近誤差e小于所述預設閾值ε₀；其中所述預設閾值為ε₀＝0.001；

模塊208記錄濾除了噪聲的信號序列S_NEW，具體為S_NEW＝x。