回归系数对a a和β a进行计笪并ο
[0102]在使用记录于该K-NET的地震的地震动数据进行回归分析的情况下，该结果为Y=3.561og10X-0.07，因此上述预测公式的a a为3.56，β a为-0.07。
[0103]此外，由于以丽I的指标值来表示各个地震的晃动的大小，因此，对于该指标值(以下称为“计算值(MMIa) ”)的计算，使用Wald等在非专利文献I所提议的计算公式。
[0104]计算公式:MMIa= a 1g10 (Amax) + β
[0105]在此，Amax为地震动最大速度的绝对值。
[0106]此外，α为 3.66，β 为-1.66。
[0107]接下来，使用记录于K-NET的地震波形数据和上述预测公式以及计算公式，分别对预测值(MMIap)和计算值(MMIa)相对于时间的时间历程变化进行模拟，并对该模拟的结果进行了比较，因此对该比较结果进行说明。
[0108]如图3所示，此次作为探讨对象的2323例地震波形数据中，预测值(MMIap)以及计算值(MMIa)均显示为丽I的指标值5.5等级以上的例子有450例。
[0109]其中，根据上述模拟，预测值(MMIap)比计算值(MMIa)先达到丽I的指标值中5.5等级的例子有384例，相反地，计算值(MMIa)先达到的例子有66例。
[0110]并且，进一步对上述384例的地震波形数据进行探讨的话，则如图4所示可知:预测值(MMIap)比计算值(MMIa)在O秒以上不足2秒的范围内先达到丽I的指标值中5.5等级的地震波形数据有126例。
[0111]此外可知:平均来说，预测值(MMIap)比计算值(MMIa)大约早5.4秒达到修正的麦加利地震烈度的5.5等级。
[0112]S卩，本实施方式的地震预测装置I使用丽I作为地震动指标，能够在地震动的初期震动阶段，提前预测是否会成为具有需要进行警报大小的晃动。
[0113]3.地震警报处理
[0114]接下来，参照图5，对由本实施方式的地震预测装置I执行的地震警报处理A进行说明。
[0115]若将地震预测装置I的未图示的电源开关接通，本实施方式的地震警报处理A便会开始，之后，在每个采样周期反复执行直到电源切断。
[0116]在该地震警报处理A中，最初将执行SlO的加速度获取处理。
[0117]在该SlO中，执行以下处理，S卩，作为通过加速度获取部10而执行的处理，输入表示从加速度传感器装置3输入的地震动3个方向(东西、南北、上下)的加速度分量的模拟信号，并进行采样。
[0118]然后，在该SlO中执行以下处理，S卩，向加速度记录部14和地震发生判断部20输出采样后而生成的表示地震动上下方向加速度分量的数字信号；向地震发生判断部20输出表不东西方向加速度分量的数字信号以及表不南北方向加速度分量的数字信号。加速度记录部14对关于该数字信号所表示的地震动上下方向加速度分量的、上下加速度信息进行记录。
[0119]接下来，执行S12的丽lap计算处理。
[0120]在该S12中，执行以下处理，即，作为通过预测值计算部16而执行的处理，从记录于加速度记录部14中的上下加速度信息中，使用最大加速度值(Aumax)，来计算预测值(MMIap)，其中该最大加速度值(Aumax)为上下方向加速度分量的绝对值中最大的绝对值。
[0121]接下来，在S14中，执行以下处理，S卩，作为通过地震发生判断部20而执行的处理，根据由加速度获取部10转换成数字信号的地震动3个方向的加速度分量，来计算观测点地震动的加速度。
[0122]接下来，执行S16的处理。
[0123]在该S16中，执行对是否发生有地震进行判断的处理。
[0124]该S16执行作为通过第I警报部18而进行的处理，具体地，将执行以下处理，即，判断储存在标志存储区域20a中的标志，是表示正在发生地震的“I”还是表示未发生地震的正常状态的“O”。
[0125]如果在该S16中，判断标志是“0”，即判断为“正常状态”(S16 ??是)，则接下来执行S18的处理。如果在S16中，判断标志是“1”，即判断为“正在发生地震”(S16:否)，则接下来执行S22的处理。
[0126]在S18中，执行以下处理，S卩，对于观测点地震动的加速度，判断该加速度的绝对值是否大于上述地震发生基准值。
[0127]该S18将通过地震发生判断部20来执行。
[0128]在该S18中，当地震动加速度的绝对值大于地震发生基准值时，即发生有地震的情形下(S18:是)，执行将储存在标志存储区域20a中的标志从“O”变更为“I”的处理(S20)，并结束该地震警报处理A，再次执行SlO以下的处理。
[0129]另一方面，当地震动加速度的绝对值小于地震发生基准值时，即未发生地震的情形下(S18:否)，则立即结束该地震警报处理A，再次执行SlO以下的处理。
[0130]接下来，对在S16中判断标志是“1”，即判断为“正在发生地震”的情形下(S16:否)要执行的S22的处理进行说明。
[0131]该S22执行以下处理，即，作为通过第I警报部18而执行的处理，判断由S12计算出的预测值(MMIap)是否为作为警报基准的警报基准值以上，即，判断是否为MMI中5.5等级以上。
[0132]在该S22中，若判断预测值(MMIap)大于警报基准值，如上所述可知，平均来说，在设置有地震预测装置I的观测点，为实际发生MMI中5.5等级以上的晃动的5.4秒之前。
[0133]因此，在S22中，当判断预测值(MMIap)为丽I中5.5等级以上时(S22:是)，则接下来执行S24的处理，执行从第I警报部18向外部警报装置5发送警报信号的处理。并且，在该S24后执行S27的处理。
[0134]另一方面，在S22中，当判断预测值(MMIap)不足丽I中5.5等级时(S22:否)，接下来执行S27的处理。
[0135]在S27中，与S18相反，执行对观测点的地震动加速度的大小是否小于预先规定的地震发生基准值进行判断的处理。
[0136]该S27通过地震发生判断部20来执行。在该S27中，与S18相同，对于观测点的地震动的加速度，则执行对该加速度的绝对值是否为上述地震发生基准值以下进行判断的处理。
[0137]通过该S27，当地震动加速度的绝对值为基准值以下时(S27:是)，则执行将储存在标志存储区域20a中的标志从“I”变更为“O”的处理(S28)，并结束该地震警报处理A，再次执行SlO以下的处理。
[0138]另一方面，当地震动加速度的大小大于基准值时(S27:否)，则立即结束该地震警报处理A，再次执行SlO以下的处理。
[0139]4.本实施方式的地震预测装置的特色性作用效果
[0140]如上所述可知:若对过去发生的地震的地震动的预测值(MMIap)和计算值(MMIa)进行比较，则在地震的初期震动阶段，预测值(MMIap)比计算值(MMIa)早5.4秒达到警报基准值。
[0141]因此，如果利用本实施方式的地震预测装置1，使用丽I作为地震动指标，在地震动的初期震动阶段，便能够对需要进行警报的地震的发生提前进行预测。
[0142]像这样如果能够对需要进行警报的地震的发生提前进行预测，例如在以270km行驶的高速铁路的情形下，由于在大的晃动到来前5.4秒便能够进行制动，因此，如果利用本实施方式的地震预测装置1，与在大的晃动到来后进行制动的情形相比，能够比其少行驶375m而停车。
[0143]此外，当通过本实施方式的地震预测装置I预测到需要进行警报的地震时，也可以使升降机停在就近层，使不耐晃动的机器停止，或通过扬声器或电视向人们进行讯息通知等。
[0144]并且，本实施方式的地震预测装置I由于使用MMI对需要进行警报的地震的发生提前进行预测，因此能够进行国际上都易于理解的地震预测。
[0145]此外，本实施方式的地震预测装置I仅在预测值(MMIap)超过预先规定的地震发生基准值的情形下才会进行警报(S22 — S24)，因此能够抑制在不需要进行警报的地震发生时无谓地进行警报。
[0146](第2实施方式)
[0147]接下来，对本发明的第2实施方式进行说明。
[0148]在本实施方式中，仅就与第I实施方式的不同点进行说明。其中，在以下对第2实施方式进行说明的叙述中，将第2实施方式称为本实施方式而进行说明。
[0149]1.地震预测装置I
[0150]如图6所示，本实施方式的地震预测装置I具备调整系数设定部22，这一点不同于第I实