素,是以何种规律相互影响、相互作用而导致地震的产生,在目前的科学界并无定论。
欧洲气象中心地震科研团队提供过来的这个地震预测模型也只是运用当前主流的各种学术理论,并根据大数据分析进行了各个影响因子的赋权。
这个数学模型本身没什么错误,可以看得出确实是一群精英科学家花了无数心血的结晶,尤其是根据当初秦克的建议,在数学模型里加入了星球内部摩擦系数与N-S方程的应用,使得它的准确率应该是领先于当今世界的。
但秦克仔细看罢所有的资料与数据分析结果后,总有种“缺了点什么、也不够完善”的直觉。
一个数学模型,哪怕本身没问题,但如果在建立之初就缺失了几个关键的参数,就会导致数据模型的失效,尤其是像地震这样的自然灾害预测模型,在理论方面并不完善,完全是基本过往数据进行分析后建立起来的,更别说这个模型当中还较多地运用了“古登堡-里克特复发关系式”(M=a-bN)的回归算法,即对历史数据进行拟合,形成拟合方程,接下来使用该方程对新数据进行预测。
作为数学建模之神的秦克对此有极深的体会,在数学建模过程中“回归算法”是常用的方法,但正因为太过常用,秦克本人反倒极少用到,因为这很容易导致“经验主义”和“主观拟合”,反倒影响到推演预测的准确率。
秦克对自己的神级直觉很有信心,既然直觉“缺了点什么,也不够完善”,那这个地震预测模型多半就“缺了点什么,也不够完善”,这个推演与预测结果也多半是有问题的。
秦克随手写下震级的公式M=lg(A/T)max 1.66lgΔ 3.5,这里A表示地震面波最大地动位移,取两个水平分向地动位移的矢量和,单位为微米(μm);T是相应周期,单位为秒(s);Δ为震中距,单位为度(°)。
他沉思了好一会,随即示意微光分析一下欧洲气象中心传来的历史地震数据。
这些数据极为详细,除了每次地震的震级(M)、发震时刻、纬度(°)、经度(°)、深度(千米)外,还有主震强震动记录,包括震中距、最大峰值加速度,强震台震中距、水平方向与垂直方向加速度峰值、速度峰值、仪器地震烈度、电离层斜测数据、电离层观测网络,以及余震数据,甚至还有新构造块体划分边界、构造地貌单元划分边界、中新生代盆地、前第四纪大地构造单元划分与主要构造界线、新生代岩浆岩分布、数字高程模型地貌阴影