[发明专利]一种基于卡方分布的微震事件检测方法和系统

权利要求书

1.一种基于卡方分布的微震事件检测方法，其特征在于，包括：

步骤1,输入实测的微震信号序列S；

步骤2，对所述微震信号序列S中每一个元素生成相对应的处理序列

s_i,i＝1,2,…,N，具体为：

s_i＝[S₁,S₂,…,S_i],i＝1,2,…,N

其中：

s_i：所述微震信号序列S中第i个元素所对应的处理序列

S_i：所述微震信号序列S中的第i个元素

N：所述微震信号序列S中元素的个数；

步骤3,求取协方差矩阵C_i,i＝1,2,…,N，具体为：

C_i＝[s_i-m_i]^T[s_i-m_i],i＝1,2,…,N

其中

步骤4,对所述协方差矩阵C_i,i＝1,2,…,N进行SVD分解，具体为：

C_i＝UΛU^T,i＝1,2,…,N

其中

U：酋阵

Λ：对角阵，对角元素为C_i的特征值；

步骤5,求取观测矢量z_i,i＝1,2,…,N，具体为：

z_i＝Λ^-1/2U^T[s_i-m_i],i＝1,2,…,N；

步骤6，求取卡方统计量具体为：

步骤7，求取判断阈值具体为：

其中

F分布的第一自由度，表示下取整

χ＝ρ+1：参量

F分布参数

F_α[ρ,N-ρ]：第一自由度为ρ、 第二自由度为N-ρ的F分布，且F分布参数为α；

步骤8,根据所述微震信号序列S的卡方统计量判断微震事件；具体为：如果所述的卡方统计量则在所述微震信号序列S的第K点处检测到微震事件，否则未检测到微震事件；其中，为判断阈值，α为F分布参数。