CN106021710B - 基于大气电离层参数的震前卫星轨道异常识别方法 - Google Patents
基于大气电离层参数的震前卫星轨道异常识别方法 Download PDFInfo
- Publication number
- CN106021710B CN106021710B CN201610330097.1A CN201610330097A CN106021710B CN 106021710 B CN106021710 B CN 106021710B CN 201610330097 A CN201610330097 A CN 201610330097A CN 106021710 B CN106021710 B CN 106021710B
- Authority
- CN
- China
- Prior art keywords
- earthquake
- parameter
- parameters
- satellite
- orbit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000005856 abnormality Effects 0.000 title abstract description 7
- 238000001514 detection method Methods 0.000 claims abstract description 10
- 239000005433 ionosphere Substances 0.000 claims description 26
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000011156 evaluation Methods 0.000 claims description 9
- 230000035939 shock Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000012886 linear function Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 230000007774 longterm Effects 0.000 claims 1
- 238000012545 processing Methods 0.000 claims 1
- 239000002243 precursor Substances 0.000 abstract description 7
- 230000000153 supplemental effect Effects 0.000 abstract 1
- 230000002159 abnormal effect Effects 0.000 description 15
- 238000005070 sampling Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000002547 anomalous effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于大气电离层参数的地震前兆卫星轨道异常识别方法，包括：原始电离参数的野值剔除及标准化处理；定义卫星轨道上电离参数偏度和峰度；计算卫星各轨道上参数的偏度和峰度值；使用t location‑scale分布拟合得到分布参数及拟合误差；在不同地震震级下，使用不同阈值划分概率密度函数，建立震前卫星异常轨道识别模型；通过同历史地震目录信息对比，以证明本方法的有效和准确性；使用综合评价方法得到面向多参数的轨道识别模型。本发明的优点是：针对大气层电离参数数据量巨大、信息维度高等特点，给出了一种地震前兆卫星轨道异常的识别方法，可用于卫星轨道异常检测和地震前兆识别。
Description
技术领域
本发明涉及基于大气电离层参数的地震震前卫星轨道异常识别方法，是强震监测大气电离层参数时卫星轨道异常的识别方法，属于航天遥感、地球物理学与计算机科学的交叉领域。
背景技术
我们生活居住的地球每日都有近万次地震发生。作为世界上最严重的自然灾害之一，地震直接导致及间接引发的次生灾害对全球带来了极大的生命和财产损失。我国位于环太平洋地震带和欧亚地震带交汇，是一个地震灾害严重的国家。如何分析、识别地震前兆，提高地震预测精度，对降低地震灾害、减少国家人民损失有着极为重要的意义。
据有关统计，全球每年约发生500多万次地震。近年来全球再次进入地壳板块运动活跃期，屡屡发生高等级强震。仅2015年全球就发生了19次7级以上强震，2016年初至今半年内共发生7级以上强震6次。我国位于亚欧板块、印度洋板块和太平洋板块交界处，欧亚地震带和环太平洋地震带交汇处，地震发生频率高、强度大、分布广、震源浅，是一个地震灾害严重的国家。因此，分析地震前兆中自然环境因素的变化，提高地震预测水平，对降低地震灾害带来的损失具有巨大现实意义。
由于地震突发性强、发生频度高等特点，如何提高地震预报精度和可靠度始终是一个世界性难题。传统方法的空间局限性，无法适用于大面积地域的地震预报。然而随着卫星遥感技术的发展，气体反演、红外遥感等技术被应用于地震的宏观预测，其精度高、重复观察周期短、不受地面条件限制等特点，成为了研究断裂活动性和地震前兆的重要观测手段，极大提高了地震监测能力。
观测表明，在地质学、地球物理学和地球化学等多种地球前兆表现中，电磁异常反应是最为敏感的观测手段之一。在地球100公里至400公里高空，空气极其稀薄，这个高度的空间称为“电离层”。电离层的气体元素都呈游离状态出现，当地磁发生不规则变化时，高层气体元素的核外电子出现丢失或整合现象，空间占用减少，出现了电离层的向下压缩现象。目前研究已证实强震前震区上空的电离层厚度会出现往下压缩的现象，这种压缩会使“太阳风”的高能带电粒子沿着被压缩的电离层向下沉降，从而引起大气电离层参数发生一定程度的突变。从上个世纪70年代至今，已有数百位地震电磁学领域的专家学者研究探索了电离层电磁异常扰动与地震事件间的存在的关联。大量实践研究基础表明部分电磁参数的异常扰动可以作为地震前兆。
对地观测卫星呈周期性围绕地球采集地球表面的参数信息特征，每一条卫星轨道包含了一段完整的信息采集过程。传统的震前异常识别方法大多基于卫星采样点，若采样点附近出现了电离层参数异常(如采样点位于震中附近)，则将该采样点作为异常信息进行处理。但由于地震的偶发性，异常点在全部采样点中所占的比例往往较小，后续震前异常模式识别工作较难进展。而将包含一个或多个异常信息采样点的卫星轨道作为研究单位，即将该轨道标记为震前轨道异常进行研究，在不丢失异常采样点信息的前提下，大大提高了异常数据所占比例，更容易展开震前识别和分析。
发明内容
本发明的目的是基于大气电离层参数，发明一个地震震前卫星轨道异常的识别方法。基于该方法，可有效且准确地识别强震震前记录异常大气电离层参数的卫星轨道，实现了在不丢失异常采样点信息的前提下，扩大异常电离层参数所占比例，提高了强震震前异常识别的能力。
本发明的具体技术方案包括以下几个步骤：
步骤一：数据预处理。针对卫星检测得到的大气电离层参数原始数据，剔除数据测量过程中由于仪器故障等因素产生的野值；卫星通常装载了多个电离层参数的测量仪器，不同参数的量纲往往有所不同，为消除不同量纲属性对数值分析过程造成的影响，利用线性函数对原始数据进行归一化；
步骤二：基于给出的卫星轨道上电离层参数偏度和峰度的新定义，针对预处理后的电离数据计算卫星各轨道上参数的偏度和峰度值；
步骤三：使用t location-scale分布拟合计算得到的各轨道电离层参数的偏度和峰度，得到拟合参数；使用残差平方和SSE(Sum of Squares for Error)作为拟合误差的衡量标准；
步骤四：根据步骤三得到的参数拟合分布，在不同地震震级下划分概率密度函数尾部并记录划分阈值；根据得到的划分阈值序列，建立震前卫星轨道异常识别模型；
步骤五：针对步骤四得到的各电离层参数轨道异常识别模型，结合历史地震目录信息，使用残差平方和SSE和判定系数R2，以证明模型的有效和准确；
步骤六：分别建立不同电离层参数的轨道异常识别模型，根据不同参数模型判定系数R2定义各自权重，使用综合评价方法结合各子模型得到一个综合考虑多参数的识别模型。
本发明的有益效果是：
本发明所提的方法从卫星轨道异常的角度出发，针对监测得到的电离层参数信息进行地震震前的卫星轨道异常变化分析，提出了一种基于大气电离层参数的地震震前卫星轨道异常 识别方法，有效地解决了大气电离异常信息数量少，异常程度低等问题。本发明经多次真实卫星数据测试验证以及相关专家评定，能有效识别强震震前的电离轨道异常。
附图说明
图1是本发明方法的总体流程图。
具体实施方式
下面结合附图和相关算法，对本发明做进一步的说明。
本发明的总体流程如图1所示。
本发明针对大气电离数据，剔除原始参数中野值并做标准化处理；定义了卫星轨道上电离层参数偏度和峰度；计算卫星各轨道上参数的偏度和峰度值；使用t location-scale分布拟合得到分布参数及拟合误差；在不同地震震级下，使用不同阈值划分概率密度函数，建立震前卫星轨道异常识别模型；通过同历史地震目录信息对比，以证明本方法的有效和准确性；使用综合评价方法得到面向多参数的轨道识别模型。具体实施步骤如下：
1.野值剔除及数据标准化
针对本发明中存在的小部分异常数据，使用切比雪夫滤波器进行野值剔除：
式中n表示I型切比雪夫滤波器的阶数，ω0表示滤波器截止频率，ε表示滤波器参数且满足|ε|＜1，Tn(x)＝cos(n·arccos(x))表示n阶切比雪夫多项式。
同样，由于卫星装载的多个测量仪器对不同测量参数的量纲有所不同，为消除不同量纲属性对数值分析过程造成的影响，使用线性函数对原始数据进行归一化。
式中x为原始数据值，xmin和xmax分别表示属性的最大值和最小值，x′为归一化后的数值，其范围在0到1之间。
其过程可用算法1来描述：
算法1：数据野值剔除及其标准化
输入：原始电离层参数D
输出：预处理后的电离层参数
01： D’＝[]；//数据初始化
02： for i＝l：k //共有k个不同的电离层参数
03： X＝D(：，i)； //取出数据D中的第i个属性
04： X’＝Chev(X)； //经过车比雪夫滤波器，实现野值剔除
05： X”＝(X’-min(X’))/(max(X’)-min(X’))； //数据标准化
06： D’＝[D’；X”]； //将预处理后的第i个电离层参数存入D’中
07： endfor
2.卫星轨道上电离层参数的偏度和峰度
本发明完成各电离层参数的预处理后，进一步定义并计算各卫星轨道上参数的偏度和峰度。针对各卫星轨道上电离检测点，每一个检测点对应的地磁纬度为lat，对应电离层参数为f(lat)，同一轨道上所有检测点地磁纬度的集合为LATset，其均值记为下面给出卫星轨道上电离层参数偏度的计算方法：
同样给出卫星轨道上电离层参数峰度的计算方法：
其过程可用算法2来描述：
算法2：计算各卫星轨道上电离层参数的偏度和峰度
输入：原始电离层参数D
输出：计算得到的偏度A3和峰度A4
01： A3＝[]；A4＝[]； //数据初始化
02： for i＝l：k //共有k个不同的电离层参数
03： START＝min(D(：，i).orbitnum)； //标记第i个电离层参数的初始轨道序号
04： STOP＝max(D(：，j).orbitnum)； //标记第i个电离层参数的终止轨道序号
05： X＝[]； //初始化
05： for j＝START：STOP //遍历所有轨道序号
06： for w＝1：length(D(：，i)) //遍历所有轨道
07： if D(w，i).orbitnum＝＝j //属于轨道j
08： X＝[X；D(w，j)]； //存入同一类
09： endif
10： endfor
11： alpha3＝sum(X.para.*(X lat)^3)/(sum(X.para.*X.lat^2))^1.5；//计算偏度值
12： alpha4＝sum(X.para.*(X.lat)^4)/(sum(X.para.*X.lat^2))^2；//计算峰度值
13： A3＝[A3；alpha3]；A4＝[A4；alpha4]；//记录偏度、峰度
14： endfor
15： endfor
3.t location-scale分布拟合
针对计算得到的卫星轨道上电离层参数的偏度和峰度，本发明使用t location-scale分布对其进行拟合。下面给出t location-scale分布概率密度函数的计算表达公式：
式中的μ为位置参数，σ为尺度参数，v为形状参数，Γ(□)为伽马函数。
针对给出的偏度和峰度数据，计算不同数值区间下出现的频数和频率，根据以上给出的t location-scale分布概率密度函数，使用非线性最小二乘法进行参数拟合：
minΣ[y-f(x|μ，σ，v)]2。
式中y表示数据x出现的频率值，f(x|μ，σ，v)表示数据x出对应t location-scale分布的概率密度值。保证拟合值与真实值间的平方和最小，由此估计得到概率分布的参数μ，σ和v。
使用残差平方和SSE(Sum of Squares for Error)作为拟合误差的衡量标准，其计算公式如下：
式中x表示数据的真实值，x′表示对应的拟合值，X表示所有数据x构成的集合。误差平方和SSE越小，表明拟合参数更精确，拟合效果更好。
其具体过程如算法3所述：
算法3：t location-scale分布拟合
输入：计算得到的偏度A3和峰度A4
输出：分布拟合参数MU、SIGMA、NU，拟合残差平方和SSE
01： MU＝[]；SIGMA＝[]；NU＝[]；SSE＝[]； //数据初始化
02： for i＝3：4 //3表示偏度，4表示峰度
03： for j＝l：k //共有k个不同的电离层参数
04： mu＝0；sigma＝0；nu＝0；sse＝0； //数据初始化
05： [mu，sigma，nu，sse]＝dfittool(t location-scale，Ai(：，j))；//使用非线性最小二乘法分布拟合
06： MU＝[MU；mu]；SIGMA＝[SIGMA；sigma]；NU＝[NU；nu]；SSE＝[SSE；sse]；//数据存储
07： endfor
08： endfor
4.建立震前卫星轨道异常识别模型
针对t location-scale分布拟合结果，在不同地震震级下，划分概率密度函数尾部并记录划分阈值。对于一个给定的地震震级M0，查询历史地震目录，以地震发生时刻t为横坐标，地震发生数量为纵坐标，绘制地震发生时间与数量的曲线图划分阈值从概率密度函数尾部两端开始逐步向分布中心靠拢，记录划分阈值为±x时，尾部所包含的所有轨道序号，以轨道时间t为横坐标，轨道总数量N±x为纵坐标，绘制轨道异常时间与数量的曲线图N±x＝g(t)。计算与N±x间的判定系数R2
式中表示N±x的均值。
记录给定震级M0下最小R2对应的±x值作为此时的划分阈值。
根据得到的划分阈值，建立震前卫星轨道异常识别模型：
式中M0表示给定的地震震级，表示截止至时间t共发生震级大于M0的地震数量，a，b和c为分别为模型的尺度参数、形状参数和位置参数，用于控制不同震级限制下的模型的识别精度。
其具体过程如算法4所述：
算法4：建立震前卫星轨道异常识别模型
输入：计算得到的分布拟合参数MU、SIGMA、NU，历史地震目录表EARTHQUAKE
输出：轨道识别模型参数A、B、C
01： A＝[]；B＝[]；C＝[]；THRESHOLD＝[]； //数据初始化
02： for i＝l：k //共有k个不同的电离层参数
03： for M＝5：0.1：8 //震级选择
04： DATA_M＝select(EARTHQUAKE，M)； //筛选震级大于M的历史震例
05： N_X＝[]；N_M＝[]； //数据初始化
06： for w＝5：0.1：M
07： N_X＝[N_X；count(DATA_M＞w)]； //计算并记录DATA_M表中震级大于w的震例数量
08： N_M＝[N_M；model(MU(i，w)，SIGMA(i，w)，NU(i，w))]； //计算并记录模型中震级大于w的震例数量
09： endfor
10： R2＝sum((N_M-mean(N_X))^2)/sum((N_X-mean(N_X))^2)；//计算判定系数R2
11： if R2＞max_r2
12： max_r2＝R2；temp＝w；//记录当前划分阈值
13： endif
14： THRESHOLD＝[THRESHOLD；w]；//记录划分阈值
15： [a，b，c]＝cftool(THRESHOLD，exp，2)；//二阶指数拟合
16： endfor
17： A＝[A；a]；B＝[B；b]；C＝[C；c]；
18： endfor
5.检验模型的有效和准确性
针对得到的震前卫星轨道异常识别模型，查询历史地震目录信息，分别筛选统计截止至时间t所有发生震级大于M0地震的数量。以地震发生时刻t为横坐标，地震发生数量Ny为纵坐标，绘制地震发生时间与数量的曲线图Nt＝f(t)。计算与震前卫星轨道异常识别模型间的残差平方和SSE和判定系数R2，其计算公式如下：
式中表示Nt的均值。残差平方和SSE和判定系数R2越小，表示得到的模型精度越高，轨道异常识别效果越好，更能反映震前异常状况。
其具体过程如算法5所述：
算法5：建立震前卫星轨道异常识别模型
输入：轨道识别模型参数A、B、C，历史地震目录表EARTHQUAKE
输出：残差平方和SSE、判定系数R2
01： SSE＝[]；R2＝[]；//数据初始化
02： for i＝3：4//3表示偏度，4表示峰度
03： for j＝l：k//共有k个不同的电离层参数
04： N_X＝[]；N_M＝[]；//数据初始化
05： for M＝5：01：8//震级选择
06： N_X＝[N_X；count(EARTHQUAKE＞w)]；//计算并记录EARTHQUAKE表中震级大于w的震例数量
07： N_M＝[N_M；model(A(i，M)，B(i，M)，C(i，M))]；//计算并记录模型中震级大于M的震例数量
08： endfor
09： SSE＝[SSE；sum((N_X-N_M)^2)]；//计算并记录SSE
10： R2＝sum((N_M-mean(N_X))^2)/sum((N_X-mean(N_X))^2)；//计算并记录判定系数R2
11： endfor
12： endfor
6.使用综合评价方法得到面向多参数的轨道识别模型
针对不同的电离层参数，已经得到多个含有不同参数值得轨道识别模型。使用一种综合评价的方法结合以上多个子模型，以获得一个面向多参数的轨道识别模型。
现有k个电离层参数分别建立轨道异常识别模型各模型分别计算得到一个判定系数Ri 2。使用综合评价方法设定各子模型的权重
得到最后的综合模型为
其具体过程如算法6所述：
算法6：建立面向多参数的轨道识别模型
输入：轨道识别模型参数A、B、C，判定系数R2
输出：最终模型参数PARA
01： PARA＝[]；W＝[]；N_C＝[]； //数据初始化
02： for i＝l：k //共有k个不同的电离层参数
03： W＝[W；R2(i)/sum(R2)]； //计算并记录子模型权重
04： endfor
05： for M＝5：0.1：8 //震级选择
06： temp＝0； //数据初始化
07： for i＝l；k //共有k个不同的电离层参数
08： temp＝temp+W(i)*model(A(i，M)，B(i，M)，C(i，M)，M)； //子模型加权求和
09： endfor
10： N_C＝[N_C；temp]；
11： endfor
12： PARA＝cftool(N_C，exp，2)； //二阶指数拟合。
Claims (6)
1.一种基于大气电离层参数的地震震前卫星轨道异常识别方法，其特征包括如下步骤：
(1)针对卫星检测得到的大气电离层参数原始数据，剔除数据测量过程中由于仪器故障因素产生的野值，并做标准化处理；
(2)基于给出的卫星轨道上电离层参数偏度和峰度的新定义，针对预处理后的电离数据计算卫星各轨道上参数的偏度和峰度值；
(3)在步骤(2)计算得到的偏度、峰度数值基础上，使用t-location-scale分布拟合得到拟合参数及拟合误差；
(4)根据步骤(3)中得到的参数拟合分布，在不同地震震级下划分曲线尾部并记录划分阈值，建立震前卫星轨道异常识别模型；
(5)针对步骤(4)中得到的模型，结合历史地震目录信息，以证明本方法的有效和准确性；
(6)使用步骤(5)的方法，分别建立不同电离层参数的轨道异常识别模型，使用综合评价方法，结合各子模型得到一个综合考虑多参数的识别模型；
其中，步骤(2)计算卫星各轨道上参数偏度和峰度值，其实现方法如下：
(21)卫星轨道上电离层参数的偏度，针对各卫星轨道上电离层参数检测点，每一个检测点对应的地磁纬度为lat，对应电离层参数为f(lat)，同一轨道上所有检测点地磁纬度的集合为LATset，其均值记为下面给出卫星轨道上电离层参数偏度的计算方法：
其中X为样本集，其均值记为
(22)卫星轨道上电离层参数的峰度，同样给出卫星轨道上电离层参数峰度的计算方法：
2.根据权利要求1所述方法，其特征在于，步骤(1)针对卫星检测得到的大气电离层参数原始数据，剔除数据测量过程中由于仪器故障因素产生的野值，其实现方法如下：
(11)由于卫星测量大气电离层参数是一个长期不间断的过程，期间不可避免产生仪器故障的情况；这导致原始数据中存在小部分不准确的异常数据或野值，使用切比雪夫滤波器能够去除其中大部分噪声数据：
式中n表示I型切比雪夫滤波器的阶数，ω0表示滤波器截止频率，ε表示滤波器参数且满足|ε|＜1，Tn(x)＝cos(n·arccos(x))表示n阶切比雪夫多项式，ω表示滤波器的现有频率，Gn(ω)表示切比雪夫滤波器的幅度与频率关系；
(12)电离层参数测量仪器多样，不同参数的量纲往往有所不同，为消除不同量纲属性对数值分析过程造成的影响，利用线性函数对原始数据进行归一化：
式中x为原始数据值，xmax和xmin分别表示属性的最大值和最小值，x′为归一化后的数值，其范围在0到1之间。
3.根据权利要求1所述的方法，其特征在于，步骤(3)使用t-location-scale分布拟合计算得到拟合参数及拟合误差，其实现方法如下：
(31)t-location-scale分布
下面给出t-location-scale分布概率密度函数的计算表达公式：
式中的μ为位置参数，σ为尺度参数，v为形状参数，Γ(·)为伽马函数；
(32)分布拟合方法
针对给出的偏度和峰度数据，计算不同数值区间下出现的频数和频率，根据步骤(31)给出的t-location-scale分布概率密度函数，使用非线性最小二乘法进行参数拟合：
min∑[y-f(x|μ，σ，v)]2；
式中y表示数据x出现的频率值，f(x|μ，σ，v)表示数据x对应t-location-scale分布的概率密度值；
计算拟合值与真实值间的平方和∑[y-f(x|μ，σ，v)]2，并保证其数值最小；将此刻的参数μ，σ和v作为t-location-scale分布的估计参数；
使用残差平方和，英文全称为Sum of Squares for Error简称SSE，作为拟合误差的衡量标准，其计算公式如下：
式中x表示原始数据值，x′表示对应的拟合值，X表示所有数据x构成的集合；误差平方和SSE越小，表明拟合参数更精确，拟合效果更好。
4.根据权利要求1所述的方法，其特征在于，步骤(4)在不同地震震级下划分概率密度函数尾部并记录划分阈值，建立震前卫星轨道异常识别模型，其实现方法如下：
(41)选取不同震级下的划分阈值
对于一个给定的震级M0，查询历史地震目录，以地震发生时刻t为横坐标，地震发生数量为纵坐标，绘制地震发生时间与数量的曲线图划分阈值从概率密度函数尾部两端开始逐步向分布中心靠拢，记录划分阈值为±x时，尾部所包含的所有轨道序号，以轨道时间t为横坐标，轨道总数量N±x为纵坐标，绘制轨道异常时间与数量的曲线图N±x＝g(t)；计算与N±x间的判定系数R2
式中表示N±x的均值；
记录给定震级M0下最小R2对应的±x值作为此时的划分阈值；
(42)震前卫星轨道异常识别模型
根据步骤(41)得到的划分阈值，建立震前卫星轨道异常识别模型：
式中M0表示给定的地震震级，表示截止至时间t共发生震级大于M0的地震数量，M表示发生的地震震级，a，b和c为分别为模型的尺度参数、形状参数和位置参数，用于控制不同震级限制下的模型的识别精度。
5.根据权利要求1所述的方法，其特征在于，步骤(5)结合历史地震目录信息，以证明本方法的有效和准确性，其实现方法如下：
查询历史地震目录信息，分别筛选统计截止至时间t所有发生震级大于M0地震的数量；以地震发生时刻t为横坐标，地震发生数量Ny为纵坐标，绘制地震发生时间与数量的曲线图Nt＝f(t)；计算与震前卫星轨道异常识别模型间的残差平方和SSE和判定系数R2，其计算公式如下：
式中表示Nt的均值，表示截止至时间t共发生震级大于M0的地震数量，Nt表示的是t时刻地震的发生数量；残差平方和SSE和判定系数R2越小，表示得到的模型精度越高，轨道异常识别效果越好，更能反映震前异常状况。
6.根据权利要求1所述的方法，其特征在于，步骤(6)分别建立不同电离层参数的轨道异常识别模型，使用综合评价方法结合各子模型得到一个综合考虑多参数的识别模型，其实现方法如下：
现有k个电离层参数分别建立轨道异常识别模型各模型分别计算得到一个判定系数式中ai表示第i个电离层参数的尺度参数，bi表示第i个电离层参数的形状参数，ci表示第i个电离层参数的位置参数，t表示t时刻，M表示t时刻的地震级数；使用综合评价方法设定各子模型的权重：
得到综合模型为：
其中，表示位置参数为c时共发生震级大于M0的地震数量，M0表示给定的地震震级，表示电离层参数个数为i时共发生震级大于M0的地震数量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610330097.1A CN106021710B (zh) | 2016-05-13 | 2016-05-13 | 基于大气电离层参数的震前卫星轨道异常识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610330097.1A CN106021710B (zh) | 2016-05-13 | 2016-05-13 | 基于大气电离层参数的震前卫星轨道异常识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106021710A CN106021710A (zh) | 2016-10-12 |
| CN106021710B true CN106021710B (zh) | 2019-08-09 |
Family
ID=57097682
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610330097.1A Active CN106021710B (zh) | 2016-05-13 | 2016-05-13 | 基于大气电离层参数的震前卫星轨道异常识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106021710B (zh) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108627855B (zh) * | 2018-05-15 | 2020-04-28 | 江苏海洋大学 | 台风中心电离层tec异常的滑动三次曲线探测方法 |
| CN108845348B (zh) * | 2018-07-10 | 2024-01-09 | 北京大学深圳研究生院 | 用于地震预报的前兆数据处理方法、地震预报方法和系统 |
| CN109117449B (zh) * | 2018-07-27 | 2022-04-15 | 武汉文网亿联科技有限公司 | 基于非线性最小二乘模型测算网吧安装率的方法 |
| CN109507724B (zh) * | 2018-11-01 | 2019-10-11 | 中国石油大学(华东) | 一种基于非震动态背景场的地震tec异常信息提取方法 |
| CN109740453B (zh) * | 2018-12-19 | 2022-03-29 | 吉林大学 | 一种基于小波变换的卫星磁场数据地震前兆异常提取方法 |
| CN110068857B (zh) * | 2019-04-02 | 2020-02-04 | 吉林大学 | 基于主成分分析的Swarm双星磁场数据地震前兆异常提取方法 |
| CN111830451B (zh) * | 2019-04-19 | 2023-04-18 | 中国科学院光电研究院 | 一种非成像传感器的检验方法 |
| CN110146921B (zh) * | 2019-06-28 | 2021-06-11 | 广东石油化工学院 | 基于狄拉克分布概率的微震事件检测方法和系统 |
| CN110673206B (zh) * | 2019-08-26 | 2020-12-29 | 吉林大学 | 一种基于非负矩阵分解的卫星磁场数据地震异常检测方法 |
| CN113325469B (zh) * | 2020-02-28 | 2023-10-13 | 中国科学院国家空间科学中心 | 一种地震电离层tec关联分析方法 |
| CN111814699B (zh) * | 2020-07-13 | 2023-07-28 | 中国地震局地震预测研究所 | 一种面向swarm电磁卫星数据的深度学习地震预测方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103680069A (zh) * | 2012-09-12 | 2014-03-26 | 李宗洪 | 地震预报机、地震预报系统和地震前兆监测系统 |
| CN104865602A (zh) * | 2014-08-29 | 2015-08-26 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种通过多数据源确定地震震中位置和起震时间的方法 |
-
2016
- 2016-05-13 CN CN201610330097.1A patent/CN106021710B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103680069A (zh) * | 2012-09-12 | 2014-03-26 | 李宗洪 | 地震预报机、地震预报系统和地震前兆监测系统 |
| CN104865602A (zh) * | 2014-08-29 | 2015-08-26 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种通过多数据源确定地震震中位置和起震时间的方法 |
Non-Patent Citations (3)
| Title |
|---|
| Study of the lower hybrid resonance frequency over the regions of gathering earthquakes using DEMETER data;D.I.Vavilov 等;《Journal of Atmospheric and Solar-Terrestrial Physics》;20130830;第100-101卷;第1-12页 * |
| 利用DEMETER卫星数据提取地震异常的方法初探;欧阳新艳;《中国优秀硕士学位论文全文数据库基础科学辑》;20090115(第1期);论文全文 * |
| 基于大气遥感的地震异常特征规律和地震预报模型;臧思聪 等;《2015年（第四届）全国大学生统计建模大赛论文》;20151129;第1-22页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106021710A (zh) | 2016-10-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106021710B (zh) | 基于大气电离层参数的震前卫星轨道异常识别方法 | |
| Hovgaard et al. | Reducing statistical noise in airborne gamma-ray data through spectral component analysis | |
| Pulinets et al. | Correlation analysis technique revealing ionospheric precursors of earthquakes | |
| CN111651941A (zh) | 一种全球电离层电子总含量预测的算法 | |
| CN110673206B (zh) | 一种基于非负矩阵分解的卫星磁场数据地震异常检测方法 | |
| CN114879251B (zh) | 一种基于射线追踪技术的地震定位图形方法及系统 | |
| Piersanti et al. | The Pollino 2012 seismic sequence: clues from continuous radon monitoring | |
| CN103336305A (zh) | 一种基于灰色理论划分致密砂岩储层岩石物理相的方法 | |
| CN111679342A (zh) | 快速寻找盆地隐伏砂岩型铀矿方法 | |
| Zöller et al. | A systematic test on precursory seismic quiescence in Armenia | |
| CN117037421A (zh) | 降雨型滑坡隐患群气象风险预警方法、设备及存储介质 | |
| Zhang et al. | A preliminary evaluation of surface latent heat flux as an earthquake precursor | |
| Wyss et al. | Seismic quiescence precursors to two M7 earthquakes on Sakhalin Island, measured by two methods | |
| Kossobokov et al. | Active zone size versus activity: A study of different seismicity patterns in the context of the prediction algorithm M8 | |
| CN118311645A (zh) | 一种地震前兆信息提取方法、装置、设备及存储介质 | |
| CN118039029A (zh) | 基于机器学习和锆石成分识别花岗岩类型的方法及系统 | |
| CN117310829A (zh) | 一种基于地磁异常数据的优化识别方法 | |
| Sidorov et al. | A combined analysis of geomagnetic data and cosmic ray secondaries for the September 2017 space weather event studies | |
| Lyu et al. | Estimation of PM2. 5 concentrations over Beijing with MODIS AODs using an artificial neural network | |
| Luo et al. | Sequential probability ratio test using scaled time-intervals for environmental radiation monitoring | |
| Jia et al. | The structure and ability of the China seismological gravity monitoring system | |
| CN116381773B (zh) | 一种在地震预报中归一化混杂数据的方法及装置 | |
| Mehmood et al. | Tukey control chart for radon monitoring in relation to the seismic activity | |
| CN116448062B (zh) | 一种桥梁沉降变形检测方法、装置、计算机及存储介质 | |
| CN111830451A (zh) | 一种非成像传感器的检验方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |