[发明专利]一种基于共振包监测数据的地震预测方法

说明书

技术领域

本发明涉及地震预测预报技术领域，具体地，涉及一种利用地震包体前兆波数据和潮汐波共振包的监测数据，进行地震预测的方法。

背景技术

地震是人类千百年来最严重的自然灾害，人类迫切地希望能够预报地震。尽管有极少数地震曾经被预报，如1976年中国的海城地震，但绝大多数地震仍然无法预报，总是不期而遇。直到钱复业、赵玉林、赵璧如等人发明了HRT波地震预报监测系统(申请号：200820000524.0)，又称PS100地电仪，引入CDMA技术显著提高了抗干扰能力，测试精度提高了两个量阶，通过地电阻率、地电场的测量，利用引潮力产生的谐振共振波(HRT波)，钱复业、赵玉林等人提出了HRT波地震预测法，迎来了同时预测地震三要素(震中、震级和发震时间)的转机。

孙威发明的SW(如应力、谐振和地磁等)监测仪，也提供了地震预测的基本数据，预测了一些地震；刘根深利用地震构造相似性规则(SW指纹法)，几乎每天都在预测地震。因地震机理不明，钱复业、赵璧如等人，视震源为黑箱，借引潮力敲打，形成谐振共振波，发现共振波周期决定震级，快慢波到时差决定震中，月潮规则决定发震时间。PS100地电仪从2004年至2013年，监测了非常丰富的资料，曾经成功预测含反演了20多次地震，预报了如2010年7.1级玉树地震等震例。

但是在申请人在实现本发明的过程中，发现以上这些方法，由于存在理论困难，信号不易识别等诸多问题，仍然没能从根本上解决地震预测问题。

发明内容

为了解决现有技术中存在的信号识别能力差、地震预测准确度低的问题，本发明提出了一种基于共振包监测数据的地震预测方法。

本发明的基于共振包监测数据的地震预测方法包括：

连续采集一个统计周期内的共振包监测数据，所述一个统计周期为自潮汐波共振包第一次出现的前一天至潮汐波共振包最后一次出现的后一天；根据所述共振包监测数据，获取共振波振幅与潮汐波振幅的振幅比序列和共振周期序列；将所述振幅比序列与已知的地震震级及其相应振幅比的序列进行对比，计算未知地震的初步震级；根据所述初步震级和所述共振周期序列，计算未知地震的精确震级。

优选的，还包括：根据所述共振包监测数据，获取共振波振幅序列；根据快慢共振波振幅特征和所述共振波振幅序列，计算快慢波到时差序列；根据所述快慢波到时差序列计算震中距序列。

优选的，还包括：根据所述潮汐波共振包出现的时间和成长时间，计算未知地震的发震时间；根据所述共振包监测数据，获取共振包振幅对时间的导数序列，计算临震时间和预警时间。

优选的，还包括：根据多个前兆台站采集到的共振包快波和慢波到达时间差确定第一震中距，并根据所述震中距计算震中的第一空间坐标经纬度信息；根据多个前兆台站采集到的共振周期序列，计算震中的第二空间坐标经纬度信息；根据所述第一空间坐标经纬度信息和所述第二空间坐标经纬度信息，确定震中。

优选的，所述根据所述快慢波到时差序列计算震中距序列的步骤具体包括：计算快波共振波和慢波共振波到达前兆台站的时间差序列；根据所述时间差序列，计算该前兆台站的震中距。

优选的，所述根据所述潮汐波共振包出现的时间和成长时间，计算未知地震的发震时间的步骤具体包括：从出现潮汐波共振包的当天起算，未知地震的发震时间符合整数规则，将出现第6个潮汐波共振包的当天作为未知地震的发震时间。

优选的，所述根据所述共振包监测数据，获取共振包振幅对时间的导数序列，计算临震时间和预警时间的步骤具体包括：当一个共振包向下一个共振包过度时，计算所述共振包振幅对时间的导数绝对值，将所述导数绝对值趋于极大或无穷时对应的时间作为临震时间，所述临震时间的精度为12h以内，预警时间的精度为临震时间的2-9h以内。

优选的，所述根据多个前兆台站采集到的共振包快波和慢波到达时间差确定第一震中距，并根据所述震中距计算震中的第一空间坐标经纬度信息的步骤具体包括：以所述多个前兆台站为圆心，第一震中距为半径建立多个圆周，将所述多个圆周的交点作为震中的第一空间坐标经纬度信息。

优选的，所述根据多个前兆台站采集到的共振周期序列，计算震中的第二空间坐标经纬度信息的步骤具体包括：根据所述多个前兆台站采集到的共振周期，计算所述多个前兆台站的第二震中距，以所述多个前兆台站为圆心，第二震中距为半径建立多个圆周，将所述多个圆周的交点作为震中的第二空间坐标经纬度信息，所述多个前兆台站采集到的共振周期序列来自同一震源发出的共振信号，所述共振周期相等。