[发明专利]一种利用中心频率偏移的微震事件检测方法和系统

权利要求书

1.一种利用中心频率偏移的微震事件检测方法，其特征在于，包括：

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

步骤102求取N个K阶信号差分延迟矢量，具体为：第n个K阶信号差分延迟矢量记为ΔSⁿ，其第i个元素为的计算公式为其中，s_i为所述信号序列S的第i个元素；为所述信号序列S的第|i-K|_N个元素；|i-K|_N表示以N为模对i-K取余数运算；K为中频延迟参数，其计算公式为N为所述信号序列S的长度；snr为所述信号序列S的信噪比；表示以N为模下取整运算；n为序列序号，其取值范围为n＝1,2,…,N；

步骤103求取N个离散傅里叶变换系数矢量,具体为：第n个离散傅里叶变换系数矢量记为fⁿ，其计算公式为fⁿ＝DFT[ΔSⁿ]；其中，DFT表示离散傅里叶变换；

步骤104求取N个中心频率偏移，具体为：第n个中心频率偏移记为H_n，其计算公式为其中，k为求和参数，其取值范围为k＝1,2,…,n；f^k为第k个离散傅里叶变换系数矢量；f₀为所述信号序列S的中心频率；为第n个离散傅里叶变换系数矢量fⁿ的均值；

步骤105求取事件判断阈值，具体为：事件判断阈值记为ε₀，其计算公式为

步骤106判断微震事件，具体为：如果所述第n个中心频率偏移H_n大于或者等于所述事件判断阈值ε₀，则在所述信号序列S的第n点处，检测到微震事件，否则未检测到微震事件。

2.一种利用中心频率偏移的微震事件检测系统，其特征在于，包括：

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

模块202求取N个K阶信号差分延迟矢量，具体为：第n个K阶信号差分延迟矢量记为ΔSⁿ，其第i个元素为的计算公式为其中，s_i为所述信号序列S的第i个元素；为所述信号序列S的第|i-K|_N个元素；|i-K|_N表示以N为模对i-K取余数运算；K为中频延迟参数，其计算公式为N为所述信号序列S的长度；snr为所述信号序列S的信噪比；表示以N为模下取整运算；n为序列序号，其取值范围为n＝1，2，…,N；

模块203求取N个离散傅里叶变换系数矢量,具体为：第n个离散傅里叶变换系数矢量记为fⁿ，其计算公式为fⁿ＝DFT[ΔSⁿ]；其中，DFT表示离散傅里叶变换；

模块204求取N个中心频率偏移，具体为：第n个中心频率偏移记为H_n，其计算公式为其中，k为求和参数，其取值范围为k＝1,2,…,n；f^k为第k个离散傅里叶变换系数矢量；f₀为所述信号序列S的中心频率；为第n个离散傅里叶变换系数矢量fⁿ的均值；

模块205求取事件判断阈值，具体为：事件判断阈值记为ε₀，其计算公式为

模块206判断微震事件，具体为：如果所述第n个中心频率偏移H_n大于或者等于所述事件判断阈值ε₀，则在所述信号序列S的第n点处，检测到微震事件，否则未检测到微震事件。