一种基于自回归移动平均模型的震前红外长波辐射异常信息探测方法专利检索-震中自然灾害的预防与保护专利检索查询-专利查询网

一种基于自回归移动平均模型的震前红外长波 辐射异常信息探测方法

阅读：530发布：2020-05-08

专利汇可以提供一种基于自回归移动平均模型的震前红外长波辐射异常信息探测方法专利检索，专利查询，专利分析的服务。并且本发明公开了一种基于自回归移动平均模型时序预测模型 (ARIMA)的震前红外长波辐射 (OLR)异常信息探测方法，该方法预测参考背景值的精度要明显优于经典的时序预测模型，特别是时序预测精度成为预测性能一项重要评判指标，时间序列预测模型拟合效果的好坏直接关系着预测结果精度的高低，该结果将会对后续红外OLR异常信息探测造成重要影响。此外，针对经典时序预测模型方法在确定上下阈值方面存在缺陷的问题，本方法给出了一种较为合理的阈值确定策略，由此利用高精度的参考背景值和合理的上下阈值所探测出异常结果和得出的结论则更加合理。，下面是一种基于自回归移动平均模型的震前红外长波辐射异常信息探测方法专利的具体信息内容。

权利要求

1.一种基于自回归移动平均模型时序预测模型的震前红外长波辐射异常信息探测方法，其特征在于，所述方法包括如下步骤：
1)选取震级大于7.0级且距震中最近且距震前第3至15个月，间隔为1个月的红外OLR建模数据；
2)基于步骤1)中所选震前第3至15个月共13个月的红外OLR建模数据建立时间序列数据，利用时间序列分析法对序列周期进行判别，对该序列进行周期差分运算，直至使其变成平稳序列，并对其进行检验；
3)基于步骤2)中得到的平稳序列，利用ARIMA时序预测模型对其进行拟合和显著性检验，并将通过该检验的建立一个时序红外OLR短期预测模型；
4)利用步骤3)中建立的短期预测模型采用时间序列交叉验证法对其进行建模预测震前1个月的建模背景值，以发震前的时间点为起始点，5天为时间间隔，第30天为终止点。将得到的参考背景值与真实观测值作差，统计得到所占百分比在95％以上的残差值Δ，以此不受扰动正常期参考背景值的残差值Δ作为上下阈值；
5)基于步骤4)中建模预测得到的最佳建模预测区间结果，选取该区间的震前红外OLR样本数据，对新选取的样本数据中的异常数据进行剔除，并利用相邻正常数据内插替换该处异常数据，并利用建立的短期预测模型进行预测得到建模背景值P；
6)计算探测上下阈值L＝P±Δ；
7)将该时期得到的实际观测值与阈值作差得到异常探测结果，超出此范围即为异常扰动。
2.根据权利要求1所述的ARIMA时序预测模型法的震前红外OLR异常信息探测方法，其特征在于，步骤1)中确定大地震的具体方法为：依据地震震级，当震级大于7.0级时确定为大地震；红外OLR数据的位置为震中或者震中附近NOAA格网点。
3.根据权利要求1所述的ARIMA时序预测模型法的震前红外OLR异常信息探测方法，其特征在于，步骤2)中确定建立的红外OLR序列中周期和平稳检验的具体方法为：依据时间序列分析法中自相关图(ACF)和偏自相关图(PACF)来判别其序列周期；依据单位根(P值)检验方法，当P值<0.05，表明该序列为平稳序列。
4.根据权利要求1所述的ARIMA时序预测模型法的震前红外OLR异常信息探测方法，其特征在于，步骤3)中确定利用ARIMA时序预测模型对其进行拟合后通过显著性检验的具体方法为：依据时间序列分析法中白噪声检验法，当白噪声其值大于0.05，表明拟合的模型的残差是白噪声。
5.根据权利要求1所述的ARIMA时序预测模型法的震前红外OLR异常信息探测方法，其特征在于，步骤4)中所占百分比在95％以上的残差值指的是观测值与真值偏差的平方与观测次数比值的平方根RMSE，其用来衡量观测值同真值之间的偏差,能够很好地反映出预测结果的精密程度，计算公式如下：
其中，N是样本量，yt是时间t的时间序列中红外OLR观测值，是时间t的时间序列中红外OLR预测值。
6.根据权利要求1所述的ARIMA时序预测模型法的震前红外OLR异常信息探测方法，其特征在于，步骤5)中对新选取的样本数据中的异常数据进行剔除的具体方法为：依据震中或者震中附近NOAA格网点红外OLR数据为空值，认为数据异常。

说明书全文

一种基于自回归移动平均模型的震前红外长波 辐射异常信息

探测方法

技术领域

[0001] 本发明涉及地震预测领域，具体涉及一种基于自回归移动平均模型(ARIMA)的震前红外红外长波辐射(OLR)异常信息探测方法。

背景技术

[0002] 自前苏联学者分析中亚地区卫星红外遥感图像时意外发现中强地震前存在热异常现象以来，红外热异常变化与地震的关系逐渐得到各国学者的重视，并被期望能够用作短临地震的预报。大量的震例和统计分析表明，大地震在发生的前几天或前一个月的时间里，震区及其附近的红外长波辐射(OLR)均出现不同程度的OLR异常，其特征主要表现在震前一个月内红外OLR有明显的增加趋势。目前常用红外OLR异常信息探测方法有距平法、涡度分析法、指数平滑法等。

[0003] 目前常用红外OLR异常信息探测方法虽具有一定的合理性，但其只考虑了红外OLR序列数据的固有属性，而未顾及其不确定性成分，因而经典的方法预测红外OLR背景值的精度相对比较低。此外，由于经典探测方法在红外OLR异常提取方法上众多，导致同一地震的分析结果各式各样甚至互相矛盾，缺乏合理的统一标准，尤其在确定背景值的阈值时，探测结果的好坏直接关系着预测精度的高低，很容易导致探测错误，而一旦探测错误其将对后续探测结果造成重要影响。

发明内容

[0004] 针对现有技术中存在的问题，本发明的目的在于提供一种基于ARIMA时序预测模型法的震前红外OLR异常信息探测方法，该方法用ARIMA时序预测模型法能够较好地顾及其不确定性成分，以此得到更为准确的红外OLR参考背景值。

[0005] 为实现上述目的，本发明采用以下技术方案：

[0006] 一种基于ARIMA时序预测模型法的震前红外OLR异常信息探测方法，所述方法包括如下步骤：

[0007] 1)选取大地震且距震中最近且距震前第3至15个月，间隔为1个月的红外OLR建模数据；

[0008] 2)基于步骤1)中所选震前第3至15个月共13个月的红外OLR建模数据建立时间序列数据，利用时间序列分析法对序列周期进行判别，对该序列进行周期差分运算，直至使其变成平稳序列，并对其进行检验；

[0009] 3)基于步骤2)中得到的平稳序列，利用ARIMA时序预测模型对其进行拟合和显著性检验，并将通过该检验的建立一个时序红外OLR短期预测模型；

[0010] 4)利用步骤3)中建立的短期预测模型采用时间序列交叉验证法对其进行建模预测震前1个月的建模背景值，以发震前的时间点为起始点，5天为时间间隔，第30天为终止点。将得到的参考背景值与真实观测值作差，统计得到所占百分比在95％以上的残差值Δ，以此不受扰动正常期参考背景值的残差值Δ作为上下阈值；

[0011] 5)基于步骤4)中建模预测得到的最佳建模预测区间结果，选取该区间的震前红外OLR样本数据，对新选取的样本数据中的异常数据进行剔除，并利用相邻正常数据内插替换该处异常数据，并利用建立的短期预测模型进行预测得到建模背景值P；

[0012] 6)计算探测上下阈值L＝P±Δ；

[0013] 7)将该时期得到的实际观测值与阈值作差得到异常探测结果，超出此范围即为异常扰动。

[0014] 进一步，步骤1)中确定大地震的具体方法为：依据地震震级，当震级大于7.0时即可确定为大地震；选取震中或者震中附近的红外长波辐射数据。

[0015] 进一步，步骤2)中确定建立的红外OLR序列中周期和平稳检验的具体方法为：依据时间序列分析法中自相关图(ACF)和偏自相关图(PACF)来判别其序列周期；依据单位根(P值)检验方法，当P值<0.05，表明该序列为平稳序列。

[0016] 进一步，步骤3)中确定利用ARIMA时序预测模型对其进行拟合后通过显著性检验的具体方法为：依据时间序列分析法中白噪声检验法，当白噪声其值大于0.05，表明拟合的模型的残差是白噪声。

[0017] 进一步，步骤4)中所占百分比在95％以上的残差值指的是观测值与真值偏差的平方与观测次数比值的平方根RMSE，其用来衡量观测值同真值之间的偏差,能够很好地反映出预测结果的精密程度，计算公式如下：

[0018]

[0019] 其中，N是样本量，yt是时间t的时间序列中红外OLR观测值，是时间t的时间序列中红外OLR预测值。

[0020] 进一步，步骤5)中对新选取的样本数据中的异常数据进行剔除的具体方法为：依据震中或者震中附近NOAA格网点红外OLR数据为空值，认为数据异常。

[0021] 本发明具有以下有益技术效果：

[0022] 本申请的方法预测背景值的精度要明显优于经典的ETS、TBATS、NNAR、和Holt-Winters等时序预测模型法，特别是经典方法的预测结果存在较大系统偏差，该缺陷将会对后续红外OLR异常信息探测造成重要影响。此外，针对经典方法在确定上下阈值方面存在缺陷的问题，本方法给出了一种较为合理的阈值确定策略，由此利用高精度的背景值和合理的上下阈值所探测出异常结果和得出的结论则更加合理。附图说明

[0023] 图1为本发明实施例的流程图。

具体实施方式

[0024] 下面，参考附图，对本发明进行更全面的说明，附图中示出了本发明的示例性实施例。然而，本发明可以体现为多种不同形式，并不应理解为局限于这里叙述的示例性实施例。而是，提供这些实施例，从而使本发明全面和完整，并将本发明的范围完全地传达给本领域的普通技术人员。

[0025] 如图1所示，本发明提供了一种基于ARIMA时序预测模型法的震前红外OLR异常信息探测方法，所述方法包括如下步骤：

[0026] 1)选取大地震且距震中最近且距震前第3至15个月，间隔为1个月的红外OLR建模数据；

[0027] 2)基于步骤1)中所选震前第3至15个月共13个月的红外OLR建模数据建立时间序列数据，利用时间序列分析法对序列周期进行判别，对该序列进行周期差分运算，直至使其变成平稳序列，并对其进行检验；

[0028] 3)基于步骤2)中得到的平稳序列，利用ARIMA时序预测模型对其进行拟合和显著性检验，并将通过该检验的建立一个时序红外OLR短期预测模型；

[0029] 4)利用步骤3)中建立的短期预测模型采用时间序列交叉验证法对其进行建模预测震前1个月的建模背景值，以发震前的时间点为起始点，5天为时间间隔，第30天为终止点。将得到的参考背景值与真实观测值作差，统计得到所占百分比在95％以上的残差值Δ，以此不受扰动正常期参考背景值的残差值Δ作为上下阈值；

[0030] 5)基于步骤4)中建模预测得到的最佳建模预测区间结果，选取该区间的震前红外OLR样本数据，对新选取的样本数据中的异常数据进行剔除，并利用相邻正常数据内插替换该处异常数据，并利用建立的短期预测模型进行预测得到建模背景值P；

[0031] 6)计算探测上下阈值L＝P±Δ；

[0032] 7)将该时期得到的实际观测值与阈值作差得到异常探测结果，超出此范围即为异常扰动。

[0033] 其中，步骤1)中确定大地震的具体方法为：依据地震震级，当Ms>7.0时即可确定为大地震；红外OLR数据的位置为震中或者震中附近NOAA格网点。

[0034] 步骤2)中确定建立的红外OLR序列中周期和平稳检验的具体方法为：依据时间序列分析法中自相关图(ACF)和偏自相关图(PACF)来判别其序列周期；依据单位根(P值)检验方法，当P值<0.05，表明该序列为平稳序列。

[0035] 步骤3)中确定利用ARIMA时序预测模型对其进行拟合后通过显著性检验的具体方法为：依据时间序列分析法中白噪声检验法，当白噪声其值大于0.05，表明拟合的模型的残差是白噪声。

[0036] 步骤4)中所占百分比在95％以上的残差值指的是观测值与真值偏差的平方与观测次数比值的平方根RMSE，其用来衡量观测值同真值之间的偏差,能够很好地反映出预测结果的精密程度，计算公式如下：

[0037]

[0038] 其中，N是样本量，yt是时间t的时间序列中红外OLR观测值，是时间t的时间序列中红外OLR预测值。

[0039] 步骤5)中对新选取的样本数据中的异常数据进行剔除的具体方法为：依据震中或者震中附近NOAA格网点红外OLR数据为空值，认为数据异常。

[0040] 本发明步骤2)中所涉及的ARIMA时序预测模型方法为本领域的公知常识，在此不再赘述。

[0041] 上面所述只是为了说明本发明，应该理解为本发明并不局限于以上实施示例，符合本发明思想的各种变通形式均在本发明的保护范围之内。

标题	发布/更新时间	阅读量
有发电机的用电产品	2020-05-13	300
地震即时警示方法	2020-05-13	896
一种接受地壳力学信息的元件	2020-05-18	502
一种缓震中底的生产工艺及采用该工艺生产的缓震中底	2020-05-14	128
一种采用Boosting集成学习算法对地震事件进行区分的方法	2020-05-14	216
基于电视机的地震预警信息发布方法、系统及存储介质	2020-05-15	593
基于随机地震动数据库构建的数据交互方法及数据库装置	2020-05-17	654
一种滑坡灾害多源数据一体化综合实时监测方法	2020-05-11	306
一种具有完全自复位功能的双曲面减隔震支座	2020-05-14	213
超微孔吸音板	2020-05-17	232

一种基于自回归移动平均模型的震前红外长波辐射异常信息探测方法

一种基于自回归移动平均模型的震前红外长波辐射异常信息

技术领域

背景技术

发明内容

具体实施方式

该功能需要专业版企业版VIP权限，您可以：