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