1.一种基于卡方分布的微震事件检测方法,其特征在于,包括:步骤1,输入实测的微震信号序列S;
步骤2,根据所述微震信号序列S的卡方统计量 判断微震事件。具体为:如果所述的卡方统计量 则在所述微震信号序列S的第K点处检测到微震事件,否则未检测到微震事件。其中, 为判断阈值,α为F分布参数。
2.根据权利要求1所述的方法,其特征在于,所述步骤2之前,所述方法还包括:步骤3,求取所述的卡方统计量 和判断阈值
3.根据权利要求2所述的方法,其特征在于,所述步骤3包括:步骤301,对所述微震信号序列S中每一个元素生成相对应的处理序列si,i=1,2,…,N,具体为:si=[S1,S2,…,Si],i=1,2,…,N其中:
si:所述微震信号序列S中第i个元素所对应的处理序列。
Si:所述微震信号序列S中的第i个元素N:所述微震信号序列S中元素的个数步骤302,求取协方差矩阵Ci,i=1,2,…,N,具体为:Ci=[si-mi]T[si-mi],i=1,2,…,N其中
步骤303,对所述协方差矩阵Ci,i=1,2,…,N进行SVD分解,具体为:Ci=UΛUT,i=1,2,…,N其中
U:酋阵。
Λ:对角阵,对角元素为Ci的特征值步骤304,求取观测矢量zi,i=1,2,…,N,具体为:-1/2 T
zi=Λ U[si-mi],i=1,2,…,N。
步骤305,求取所述的卡方统计量 具体为:步骤306,求取所述判断阈值 具体为:其中
F分布的第一自由度, 表示下取整。
χ=ρ+1:参量
F分布参数
Fα[ρ,N-ρ]:第一自由度为ρ,第二自由度为N-ρ的F分布,且F分布参数为α。
4.一种基于卡方分布的微震事件检测系统,其特征在于,包括:获取模块,输入实测的微震信号序列S;
判断模块,根据所述微震信号序列S的卡方统计量 判断微震事件。具体为:如果所述的卡方统计量 则在所述微震信号序列S的第K点处检测到微震事件,否则未检测到微震事件。其中, 为判断阈值,α为F分布参数。
5.根据权利要求4所述的系统,其特征在于,还包括:计算模块,求取所述的卡方统计量 和判断阈值