一种基于随机漫步的地震数据异常检测方法
【技术领域】
[0001 ]本发明涉及一种基于随机漫步的地震数据异常检测方法。
【背景技术】
[0002] 随着科学技术的发展，地震监测产生越来越多的地面监测数据和卫星监测数据， 卫星红外遥感技术作为新兴的对地观测手段已经应用到地震科学的各个领域中，然而对震 前各种监测数据异常的研究一直都缺乏大量的震例及多种遥感数据综合分析.为了更好的 分析和提取这些数据中蕴含的与地震相关的异常信息。本发明专利提出一种基于随机漫步 的地震数据异常检测方法。该发明能够分析时间序列格式的各类地震监测数据，供地震研 究人员研究这些异常跟地震之间的关系。
【发明内容】
[0003] 本发明的目的在于提供一种基于随机漫步的地震数据异常检测方法，能有效的识 别当前数据点前后的异常变化特征，可以用来研究异常特征跟地震发生之前的关系。
[0004] 为实现上述目的，本发明的技术方案是：一种基于随机漫步的地震数据异常检测 方法，首先对时间序列数据根据窗口来提取随机漫步概率分布，然后根据当前数据点前后 的2个窗口的概率变化比值来分析当前数据点前后的数据变化情况，以有效的识别当前数 据点前后的异常变化特征，用来研究异常特征与地震发生之前的关系。
[0005] 在本发明一实施例中，所述当前数据点前后的2个窗口的概率变化比值定义为 CRRW值，CRRW值越大，表示当前数据点前后的变化越大。
[0006] 在本发明一实施例中，该方法适用于时间序列地震数据的异常检测。
[0007] 在本发明一实施例中，该方法具体实现步骤如下： S1:根据等概率分布，选择一个窗口大小ws，判断当前监测数据点的值是否大于或等于 前一个数据点的值，若是，则当前监测数据点的值定义为随机漫步的向右走，否则，当前监 测数据点的值定义为随机漫步的向左走； S2:选择当前监测数据点前ws数据点及后ws数据点，并定义前ws数据点所在的第一窗 口定义为wdl，后ws数据点所在的第二窗口定义为wd2，则通过下述随机游走概率分布公式， 可分别计算获得wdl及wd2的当前数据点根据随机游走的概率分布：
其中，双的范围是
丨是组合运算，参数ja:= |; S3:定义当前数据点前后的2个窗口 wdl、wd2的概率变化比值CRRW:
其中，胃1和_m?2分别为参数碑即窗口大小WS的取值范围的下限和上限，表示对窗 口大小从1到14?2分别计算概率分布情况，然后求和得到CRRW值;CRRW代表当前数据 点前后的概率变化率;CRRW值越大，则表明当前数据点前后的异常变化越大。
[0008] 在本发明一实施例中，所述ml和祕52的取值分别为30和45。
[0009] 相较于现有技术，本发明具有以下有益效果:本发明方法能有效的识别当前数据 点前后的异常变化特征，可以用来研究异常特征跟地震发生之前的关系。
【附图说明】
[0010] 图1为本发明方法的流程图。
[0011] 图2(a)及2(b)为本发明汶川地震分析结果图，其中，图2(a)为表明无法发现异常， 图2(b)表明震前出现大的异常值。
【具体实施方式】
[0012] 下面结合附图，对本发明的技术方案进行具体说明。
[0013] 本发明的一种基于随机漫步的地震数据异常检测方法，首先对时间序列数据根据 窗口来提取随机漫步概率分布，然后根据当前数据点前后的2个窗口的概率变化比值来分 析当前数据点前后的数据变化情况，以有效的识别当前数据点前后的异常变化特征，用来 研究异常特征与地震发生之前的关系。
[0014] 在本发明中，所述当前数据点前后的2个窗口的概率变化比值定义为CRRW值，CRRW 值越大，表示当前数据点前后的变化越大。
[0015] 在本发明中，该方法适用于时间序列地震数据的异常检测。
[0016] 在本发明中，该方法具体实现步骤如下： S1:根据等概率分布，选择一个窗口大小ws，判断当前监测数据点的值是否大于或等于 前一个数据点的值，若是，则当前监测数据点的值定义为随机漫步的向右走，否则，当前监 测数据点的值定义为随机漫步的向左走； S2:选择当前监测数据点前ws数据点及后ws数据点，并定义前ws数据点所在的第一窗 口定义为wdl，后ws数据点所在的第二窗口定义为wd2，则通过下述随机游走概率分布公式， 可分别计算获得wdl及wd2的当前数据点根据随机游走的概率分布： ....
'￡ 其中，Η =柳:，?的范围是到
I组合运算，参数f ; S3:定义当前数据点前后的2个窗口 wdl、wd2的概率变化比值CRRW:
其中，福1和祕P2分别为参数K即窗口大小的取值范围的下限和上限，表示对窗口 大小从篇1到_2分别计算概率分布情况，然后求和得到CRRW值；CRRW代表当前数据 点前后的概率变化率;CRRW值越大，则表明当前数据点前后的异常变化越大。
[0017] 在本发明中，所述鲕1_和彻2的值分别为30和45。
[0018] 以下为本发明的具体实施例。
[0019] 利用汶川地震的长波辐射数据与2.5°X 2.5°分辨率为例，具体的算法步骤如下。
[0020] 根据等概率分布(步行走左到右），我们选择一个窗口大小WS= 30。如果当前的监 测数据值大于或等于以前的数据点，我们认为它是位于右。否则，该值位于左侧。窗口 1 (WD1)为当前数据点之前的WS个数据点和窗口 2(当前数据点之后的WS个数据点）。我们可以 确定rw(WDl)和rw(WD2)为当前数据点根据随机游走的概率分布。然后，我们可以计算每个 数据点的值：
其中，和分别为参数Μ即窗口大小的取值范围的下限和上限，表示对窗口 大小从繼1到分别计算概率分布情况，然后求和得到CRRW值。CRRW代表当前数据 点前后的概率变化率;CRRW值越大，则表明当前数据点前后的异常变化越大。
[0021 ] "rw"的计算方法见下式：
其中，?的范围是一政到〃，是组合运算，参数;公式中的n为窗口的 .2... :2. 大小WS。如何设置窗口大小WS，我们选择WS为从30-45的每一个值，即胃1和的值为 30和45，对他们分别计算CRRW值，然后求他们的均值。具体分析结果见附图2所示。竖线表 示汶川地震事件时间点，图2(a)无法发现异常，图2(b)在震前出现大的异常值。
[0022]以上是本发明的较佳实施例，凡依本发明技术方案所作的改变，所产生的功能作 用未超出本发明技术方案的范围时，均属于本发明的保护范围。
【主权项】
1. 一种基于随机漫步的地震数据异常检测方法，其特征在于：首先对时间序列数据根 据窗口来提取随机漫步概率分布，然后根据当前数据点前后的2个窗口的概率变化比值来 分析当前数据点前后的数据变化情况，W有效的识别当前数据点前后的异常变化特征，用 来研究异常特征与地震发生之前的关系。2. 根据权利要求1所述的一种基于随机漫步的地震数据异常检测方法，其特征在于:所 述当前数据点前后的2个窗口的概率变化比值定义为CRRW值，CRRW值越大，表示当前数据点 前后的变化越大。3. 根据权利要求1或2所述的一种基于随机漫步的地震数据异常检测方法，其特征在 于:该方法适用于时间序列地震数据的异常检测。4. 根据权利要求1所述的一种基于随机漫步的地震数据异常检测方法，其特征在于:该 方法具体实现步骤如下： Sl:根据等概率分布，选择一个窗口大小WS，判断当前监测数据点的值是否大于或等于 前一个数据点的值，若是，则当前监测数据点的值定义为随机漫步的向右走，否则，当前监 测数据点的值定义为随机漫步的向左走； S2:选择当前监测数据点前WS数据点及后WS数据点，并定义前WS数据点所在的第一窗 口定义为Wdl，后WS数据点所在的第二窗口定义为wd2，则通过下述随机游走概率分布公式， 可分别计算获得Wdl及wd2的当前数据点根据随机游走的概率分布：其:串:疑二被咬，I朗雜围是州潮巧，f是组替运算:，:参戳S3:定义当前数据点前后的2个窗口 Wdl、wd2的概率变化比值CRRW:其中，端1巧破垃分别为参数?！即窗口大小WS的取值范围的下限和上限，表示对窗口 大小从HSl到城S2分别计算概率分布情况，然后求和得到CRRW值;CRRW代表当前数据点 前后的概率变化率;CRRW值越大，则表明当前数据点前后的异常变化越大。5. 根据权利要求4所述的一种基于随机漫步的地震数据异常检测方法，其特征在于:所 芽口 14放2的审1{1女^另〇为30芽口45。
【专利摘要】本发明涉及一种基于随机漫步的地震数据异常检测方法。首先对时间序列数据根据窗口来提取随机漫步概率分布，然后根据当前数据点前后的2个窗口的概率变化比值来分析当前数据点前后的数据变化情况，以有效的识别当前数据点前后的异常变化特征，用来研究异常特征与地震发生之前的关系。本发明方法能有效的识别当前数据点前后的异常变化特征，可以用来研究异常特征跟地震发生之前的关系。
【IPC分类】G01V1/30
【公开号】CN105510971
【申请号】CN201610089255
【发明人】孔祥增, 林崧, 陈丽萍, 林新棋, 詹晓珊, 郑之
【申请人】福建师范大学
【公开日】2016年4月20日
【申请日】2016年2月18日