[发明专利]一种利用近端梯度的微震信号滤波方法和系统

权利要求书

1.所述一种利用近端梯度的微震信号滤波方法，其特征在于，包括：

步骤101获取按时间顺序采集的信号序列S；

步骤102求取固定迭代步长t，求取公式为其中，SNR为所述信号序列S的信噪比，ΔT为所述信号序列S的采用间隔；

步骤103求取中间矢量初始化值x₀。所述中间矢量初始化值x₀的第j个元素表示为其计算公式为j为第一元素序号，其取值范围为j＝1,2,···,N，为所述信号序列S的均方差，m₀为所述信号序列S的均值，表示均值为m₀，方差为的高斯随机函数，N为所述信号序列S的长度；

步骤104求取近端梯度初始化值g₀，创建迭代控制参数k并赋值为1。所述近端梯度初始化值g₀的第1个元素表示为所述近端梯度初始化值g₀的第i个元素表示为其中，i为第二元素序号，其取值范围为i＝2,3,···,N，s_i为所述信号序列S的第i个元素，s_i-1为所述信号序列S的第i-1个元素，ΔT为所述信号序列S的采样间隔；

步骤105求取中间矢量第k步值x_k，求取近端梯度第k步值g_k。所述中间矢量第k步值x_k求取公式为x_k＝x_k-1+tg_k-1；所述近端梯度第k步值g_k求取公式为

步骤106求取逼近误差e，所述逼近误差e求取公式为：

步骤107判断所述逼近误差e是否大于或等于预设阈值ε₀，得到第一判断结果。如果所述第一判断结果显示所述逼近误差e大于或等于所述预设阈值ε₀，则所述迭代控制参数k的值加1并返回所述步骤105、所述步骤106和所述步骤107，所述迭代控制参数k的值加1；直至所述第一判断结果显示所述逼近误差e小于所述预设阈值ε₀；其中所述预设阈值为ε₀＝0.001；

步骤108求取滤除了噪声的信号序列S_NEW，求取公式为

2.所述一种利用近端梯度的微震信号滤波系统，其特征在于，包括：

模块201获取按时间顺序采集的信号序列S；

模块202求取固定迭代步长t，求取公式为其中，SNR为所述信号序列S的信噪比，ΔT为所述信号序列S的采用间隔；

模块203求取中间矢量初始化值x₀。所述中间矢量初始化值x₀的第j个元素表示为其计算公式为j为第一元素序号，其取值范围为j＝1,2,···,N，为所述信号序列S的均方差，m₀为所述信号序列S的均值，表示均值为m₀，方差为的高斯随机函数，N为所述信号序列S的长度；

模块204求取近端梯度初始化值g₀，创建迭代控制参数k并赋值为1。所述近端梯度初始化值g₀的第1个元素表示为所述近端梯度初始化值g₀的第i个元素表示为其中，i为第二元素序号，其取值范围为i＝2,3,···,N，s_i为所述信号序列S的第i个元素，s_i-1为所述信号序列S的第i-1个元素，ΔT为所述信号序列S的采样间隔；

模块205求取中间矢量第k步值x_k，求取近端梯度第k步值g_k。所述中间矢量第k步值x_k求取公式为x_k＝x_k-1+tg_k-1；所述近端梯度第k步值g_k求取公式为

模块206求取逼近误差e，所述逼近误差e求取公式为：

模块207判断所述逼近误差e是否大于或等于预设阈值ε₀，得到第一判断结果。如果所述第一判断结果显示所述逼近误差e大于或等于所述预设阈值ε₀，则所述迭代控制参数k的值加1并返回所述模块205、所述模块206和所述模块207，所述迭代控制参数k的值加1；直至所述第一判断结果显示所述逼近误差e小于所述预设阈值ε₀；其中所述预设阈值为ε₀＝0.001；

模块208求取滤除了噪声的信号序列S_NEW，求取公式为