一种检测地震电离层TEC异常的方法与流程

文档序号:22340399发布日期:2020-09-25 18:12阅读:来源:国知局

技术特征:

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)的表达式是:


技术总结
本发明公开了一种检测地震电离层TEC异常的方法,其步骤为:首先是采用Savitzky‑Golay滤波对原始TEC观测值进行平滑处理,然后将TEC的观测值和平滑值做差,所得残差即为扰动。本发明方法计算简单快速,效果明显,创新性地利用Savitzky‑Golay滤波进行电离层TEC的异常检测,该方法保留了电离层TEC的相对极大值、极小值和宽度等分布特性,减少了日常波动的干扰,大幅提高了检测结果的可靠性。

技术研发人员:祝芙英;林剑;翟笃林;卿芸
受保护的技术使用者:湖北省地震局(中国地震局地震研究所)
技术研发日:2020.06.05
技术公布日:2020.09.25
当前第2页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1