[发明专利]一种基于小波变换的卫星磁场数据地震前兆异常提取方法

说明书

本发明地震监测领域，为一种基于小波变换的卫星磁场数据地震前兆异常提取方法，包括：读取卫星的磁场数据，并根据标志位选取有效数据；根据地磁指数选取地磁活动较平静的静磁轨道；通过差分和离散小波变换对静磁轨道的磁场三分量数据进行预处理，去除幅值大和变化缓慢的静磁场部分，得到磁场数据的变化情况；利用主成分分析对预处理后的磁场三分量数据进行降维处理，去除冗余信息的同时保留信号中主要的特征；对降维后的磁场数据进行连续小波变换，通过小波系数定义轨道能量强度，并利用其进行异常轨道提取。本发明弥补了单独使用其中某一分量不能充分使用数据有用信息和其异常提取结果不能全面反映地震前兆对磁场数据影响的不足。

技术领域

本发明地震监测领域，具体地来讲为一种基于小波变换的卫星磁场数据地震前兆异常提取方法。

背景技术

地震前兆指地震发生前出现的异常现象，其中包括震源及附近物质发生的地磁、地电、重力等地球物理异常。其观测方法可以包括地面观测和卫星观测两种类型，从全球范围来看地面观测台站分布较为稀疏，而卫星观测弥补了地面观测观测范围有限，且需要铺设大量台站的不足。同时岩石圈-大气层-电离层耦合机制的提出，表明地震前兆产生的异常可以通过岩石圈、大气层传播至电离层，并对其产生影响。因此，利用卫星数据提取地震前兆异常是可行的。磁场在地震孕育的过程中会发生异常变化，卫星测量的磁场数据包括标量总场和三分量矢量数据，矢量数据反映了不同方向磁场的变化情况，包含了更多可以反映地震前的异常情况的信息。

CN106021710A公开了一种基于大气电离层参数的震前卫星轨道异常识别方法，包括原始电离参数的野值剔除及标准化处理；在不同地震震级下，使用不同阈值划分概率密度函数，建立震前卫星异常轨道识别模型；使用综合评价方法得到面向多参数的轨道识别模型等内容。针对大气层电离参数数据量巨大、信息维度高等特点，给出了一种地震前兆卫星轨道异常的识别方法,可用于卫星轨道异常检测和地震前兆识别。

CN107356969A公开了一种基于卫星热红外数据及GIS的地震前兆分析方法，采用热红外数据加GIS的方式进行地震预报分析，能提高准确率，同时用于分析的数据来源较广，具有较好的研究效果。拓展了地震前兆的研究方向，利用GIS方式可将所有对地震发生可能产生影响的量化因素全都纳入空间运算当中，不仅仅局限于地质系数以及地质构造的分析。

Mehdi Akhoondzadeh等利用中值与四分位法分析了2016年4月16日的厄瓜多尔地震，对2015年11月1日至2016年4月30日的Swarm卫星群Alpha星的磁场三分量数据分别进行处理，提取出了相应的前兆异常。Angelo De Santis等分析了2015年4月25日的尼泊尔地震，对地震前后一个月内的东向磁场分量异常按天进行累计，并通过S曲线进行拟合发现存在异常个数快速增加的现象，并将加速度最大处作为前兆异常天。上述技术分别对磁场三分量数据单独进行处理或仅使用其中一个分量进行处理研究，使用的数据中存在较多冗余信息或者使用的数据不全面，未能充分高效的利用三分量数据中含有的有用信息进行地震前兆异常的提取，其异常结果仅由某一分量的数据得到，未能充分全面的反映地震前兆对磁场的影响。

发明内容

本发明所要解决的技术问题在于提供一种基于小波变换的卫星磁场数据地震前兆异常提取方法，通过主成分分析方法对磁场三分量数据进行降维处理，减少冗余数据并保留三分量数据的主要特征，将更多有用信息转换到降维数据中，并通过连续小波变换得到其时频域能量分布，进一步定义并求取轨道能量强度，有效提取降维数据中的异常轨道，排除非震影响因素后，得到地震前兆异常。本发明可以充分利用磁场三分量数据中的有用信息，全面的反映地震前兆对磁场数据的影响，弥补了上述技术单独使用其中某一分量不能充分使用数据有用信息和其异常提取结果不能全面反映地震前兆对磁场数据影响的不足。

本发明是这样实现的，

一种基于小波变换的卫星磁场数据地震前兆异常提取方法，该方法包括：

a、读取卫星的磁场数据，并根据标志位选取有效数据；