[发明专利]一种利用泛化梯度的微震信号重构方法和系统

权利要求书

1.所述一种利用泛化梯度的微震信号重构方法，其特征在于，包括：

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

步骤102求取泛化矢量初始化值x⁰，迭代控制参数k赋值为1。所述泛化矢量初始化值x⁰求取公式为x⁰＝S-m₀。其中，m₀表示所述信号序列S的均值；

步骤103求取迭代步长t₀，具体为其中，ΔT为所述信号序列S的采样间隔；SNR为所述信号序列S的信噪比；

步骤104求取梯度矢量第k-1步值g^k-1，具体为：所述梯度矢量第k-1步值g^k-1的第1个元素的值为0；所述梯度矢量第k-1步值g^k-1的第i个元素的值为其中，为泛化矢量第k-1步值x^k-1的第i个元素；为所述泛化矢量第k-1步值x^k-1的第i-1个元素；i为元素序号，其取值范围为i＝2,3,···,N；N为所述信号序列S的长度；

步骤105求取泛化梯度矢量第k-1步值G^k-1，具体为：其中，Proh[x^k-1-t_kg^k-1]为泛化映射，其值为其中，z为中间参数矢量；

步骤106求取泛化矢量第k步值x^k，具体为：x^k＝x^k-1-t₀G^k-1；

步骤107求取逼近误差e，所述逼近误差e求取公式为：e＝|x_k-x_k-1|；

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

步骤109求取滤除了噪声的信号序列S_NEW，具体为S_NEW＝x^k+[S-m₀]^T[S-m₀][x^k-S]。

2.所述一种利用泛化梯度的微震信号重构系统，其特征在于，包括：

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

模块202求取泛化矢量初始化值x⁰，迭代控制参数k赋值为1。所述泛化矢量初始化值x⁰求取公式为x⁰＝S-m₀。其中，m₀表示所述信号序列S的均值；

模块203求取迭代步长t₀，具体为其中，ΔT为所述信号序列S的采样间隔；SNR为所述信号序列S的信噪比；

模块204求取梯度矢量第k-1步值g^k-1，具体为：所述梯度矢量第k-1步值g^k-1的第1个元素的值为0；所述梯度矢量第k-1步值g^k-1的第i个元素的值为其中，为泛化矢量第k-1步值x^k-1的第i个元素；为所述泛化矢量第k-1步值x^k-1的第i-1个元素；i为元素序号，其取值范围为i＝2,3,···,N；N为所述信号序列S的长度；

模块205求取泛化梯度矢量第k-1步值G^k-1，具体为：其中，Proh[x^k-1-t_kg^k-1]为泛化映射，其值为其中，z为中间参数矢量；

模块206求取泛化矢量第k步值x^k，具体为：x^k＝x^k-1-t₀G^k-1；

模块207求取逼近误差e，所述逼近误差e求取公式为：e＝|x_k-x_k-1|；

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

模块209求取滤除了噪声的信号序列S_NEW，具体为S_NEW＝x^k+[S-m₀]^T[S-m₀][x^k-S]。