我只想当一个安静的学霸 第423节

次日清晨沈奇醒来，去到办公室继续他的计算，为此他查询了一些基础性的地质学文献。

“除了数学，其他科目都是坑。”

沈奇为了学好物理，已将化学升到了7级，现在将地理升到5级。

地震学属于物理学旗下的品牌，是固体地球物理学的一个分支。

沈奇将地理升到5级的职业级，多了解一些地质学方面的常识，对于地震学乃至固体地球物理学的深入研究会有帮助。

“地球岩流圈内部的对流取决于地球内核裂变反应产生的能量，这个过程导致板块边缘压力的产生，由于边缘的摩擦，压力能量间歇性释放而非持续性，岩板块相互挤压而非滑动，这种突发的挤压造成了我们所说的地震。”

“我需要一些地震指标参数，以建立一种全新而有效的数学模型。”沈奇有目的性的开始查阅资料。

地震强度m定义为m=lgs，这个指标非常关键，mark。

在美国南加州地震网的官网上，可以查到1983年至今的地震时序。

对于尺寸为lgs＞1.5的地震，符合一个极其简单明了的幂律衰减规律：N（S＞s）∝s^-B，s≥10^1.5，B＞0.95.

对于菲奖得主沈奇而言，这个定律太过简单，这是著名的关于地震尺寸频率的古登堡-里克特定律，里氏XX级地震中的里氏，就是古登堡-里克特定律中的美国地震学家里克特。

地震强度lgs＜2的地震通常被称为微震，人类难以感知，只能通过地震仪记录。

美国南加州地震网上的绝大部分数据采集来自微震，微震也会产生能量，能量是另一个重要指标，沈奇mark能量指标，以及美国南加州地震网所记录的，自1983年以来每一次地震发生的具体时间和坐标经纬度。

孕震过程的流体具有不可忽视的重要性，只不过流体运动在其孕震、发震、余震的过程中，难以采集精确数据，这是困扰地震界及物理界的一大难题。

沈奇忽然间清醒了，以微积分方式描述流体运动的N-S方程不仅适用于天空和海洋，它同样适用于地表之下！

理论上来说是这样的。

沈奇未加证实的理论。

一直到现在，沈奇几乎是依靠直觉在做一些有可能毫无意义的事情。

复杂性的表象需要严密的数学语言给出定量化的内在支撑。

沈奇重新梳理昨天的计算结果，全新而有效的数学模型需要一些核心工具作为灵魂。

这具灵魂应该是一个或一组数学公式、表达式、等等。

N-S方程尚未找到通用的求解公式，以得到通解。

特解倒是有一百多个，基于其中一个特解，沈奇花费数周时间，瞒着欧叶和所有人，一个人推导、计算、验证，最终得到一个复杂的等式：

1/2d/dt∫Ω∣A^s/2ω∣^2dx+v∫Ω∣A^s+1/2ω∣^2dx+α∫Ω∣A^s/2ω∣^2dx+……（B（u，u），A^sω）=＜χ1，χ2＞

根据这个核心等式，建立一个新的数学模型，沈奇编写了一个计算机程序。

系统中没有单独开通计算机科目，数学中的一些分支包含部分计算机知识与技能。

沈奇用基础的Fortran90语言编写程序，分为设定参数个数及其误差范围的主程序和利用模型计算＜χ1，χ2＞的子程序，来拟合从美国南加州地震网上下载的参数数据。

“激动人心的时刻，到了。”

沈奇深深呼吸一口，重重敲击回车键，启动程序。

程序开始执行，沈奇自编自导的构想是否具备实际意义，答案即将揭晓。

程序瞬间生成大量的矩阵、常微分方程和偏微分方程！

沈奇倒吸一口冷气！

靠他自己手动计算这些矩阵和方程，不吃不喝不睡，至少需要连续计算十年以上。

加上欧叶，两人一起干，也得五年。

主机U风扇急速转动，U快要被玩坏了。

“坚持，再坚持一会儿啊！”

沈奇大吼到，他这台高性能电脑刚换不久，好几万美元的机器不可能这么弱鸡。

好几万美元的机器汗流浃背、泪流满面。

大量的原始数据洪流还在汹涌的倾泻进来，它们是美国南加州地震网自1983年以来采集到的所有真实地震数据，以及沈奇独家发明的基于N-S方程的流体虚拟数据。

呜呜呜！

风扇狂转，它已超越五档，直奔八档而去。

U在哭泣、哀嚎，它正处于一生中最难熬的时刻。

艰难时光持续了半个小时，U终于解脱了。

机箱猛然抖动了一下，遂恢复平静。

微观、宏观、分子、地壳，流体、能量、间歇性的摩擦，看似无关的事情却又遵循某种神秘规律，沈奇自编自导的这出数理戏剧落下帷幕。

程序启动之前，沈奇甚至怀疑这是荒诞剧。

程序结束之后，沈奇鼓起了掌，这是写实题材，真实感人。

大数据计算最终得出一组极其简单的结果：＜55，-167＞、＜55，167＞、＜-55，-167＞、＜-55，167＞

沈奇早已购置了一个地球仪，就摆在他的办公桌上。

“目前我只能得到这种程度的计算结果，毕竟我用的是N-S方程的特解。”