1.一种检测地震电离层tec异常的方法,其特征在于:所述所述检测地震电离层tec异常的方法包括下列步骤:
1)采用savitzky-golay滤波对电离层tec时间序列进行平滑处理,得到tec平滑值;
2)将tec的观测值和步骤1)得到的平滑值做差,得到的tec残差即视为异常。
2.根据权利要求1所述的检测地震电离层tec异常的方法,其特征在于:所述步骤1)的具体实现方式是:
1.1)选取任一各网点tec的时间序列tec(t)(t=1~n);
1.2)采用savitzky-golay滤波对电离层tec时间序列tec(t)(t=1~n)进行平滑处理,得到tec平滑值
3.根据权利要求2所述的检测地震电离层tec异常的方法,其特征在于:所述步骤1.2)中tec平滑值
其中:
t为当前历元;
m和k分别为滤波窗口前后的历元个数;
滤波窗口的长度为(m+k-1),这里选取的长度为9;
ai为滤波系数。
4.根据权利要求3所述的检测地震电离层tec异常的方法,其特征在于:所述滤波系数ai需满足拟合误差的平方和e最小,即:
5.根据权利要求1-4任一所述的检测地震电离层tec异常的方法,其特征在于:所述步骤2)的具体实现方式是:
2.1)选取任一各网点tec的时间序列tec(t)(t=1~n);
2.2)计算tec观测值和步骤1)得到的tec平滑值的差值,所述差值是地震电离层tec扰动的时间序列δtec(t)(t=1~n);
2.3)通过地震电离层tec扰动的时间序列δtec(t)直观地看出tec是否出现异常扰动。
6.根据权利要求5所述的检测地震电离层tec异常的方法,其特征在于:所述步骤2.2)中地震电离层tec扰动的时间序列δtec(t)的表达式是: