AI革命:地震探测技术如何重塑地球科学研究

1

在2008年1月1日凌晨1点59分,加利福尼亚州卡利帕特里亚发生了一次地震。你可能从未听说过这次地震,即使当时你正住在卡利帕特里亚,也几乎感觉不到任何震动。这次地震的震级为-0.53级,震感大约相当于一辆卡车驶过。然而,这次地震之所以引人注目,并非因为它规模宏大,而是因为它微小——而我们却能够探测到它。

过去七年中,基于计算机视觉的AI工具几乎完全自动化了地震学的一项基本任务:地震探测。这项任务曾经由人工分析师完成,后来由简单的计算机程序接手,而现在则可以由机器学习工具快速自动完成。

这些机器学习工具能够比人类分析师探测到更小的地震,尤其是在城市等嘈杂环境中。地震提供了关于地球构造和未来可能发生灾害的宝贵信息。

"在最佳情况下,当你采用这些新技术时,即使使用同样的旧数据,也像是第一次戴上眼镜,突然能够看清树上的叶子一样,"《地震洞察》newsletter的联合作者Kyle Bradley说道。

我采访了几位地震科学家,他们都一致认为,在这些特定任务上,机器学习方法已经以更好的方式取代了人工工作。

"这确实令人印象深刻,"康奈尔大学教授、Bradley的合著者Judith Hubbard告诉我。

然而,接下来会发生什么则不那么确定。地震探测是地震学的基础部分,但许多其他数据处理任务尚未受到颠覆。最大的潜在影响——一直到地震预测——尚未实现。

"这确实是一场革命,"德克萨斯大学达拉斯分校的教授Joe Byrnes说,"但这场革命仍在继续中。"

地震学家的工作是什么

当某地发生地震时,震动会穿过地面,类似于声波在空气中传播的方式。在这两种情况下,都可以推断出波所穿过的材料特性。

想象一下敲击一面墙以判断它是否是空心的。因为实心墙和空心墙的振动方式不同,你可以通过声音来判断其结构。

地震也是如此。地震波以不同方式穿过不同材料(岩石、石油、岩浆等),科学家利用这些振动来成像地球内部。

科学家传统使用的主要工具是_地震仪_。这些仪器记录地球在三个方向上的运动:上下、南北和东西。如果发生地震,地震仪可以测量特定位置的震动。

传统地震仪

科学家随后处理原始地震仪信息以识别地震。

地震产生多种类型的震动,它们以不同速度传播。其中两种类型——P波(初级波)和S波(次级波)尤为重要,科学家喜欢识别每种波的起始点。

机器学习出现前的地震探测

在良好算法出现之前,地震目录必须手动创建。Byrnes说:"传统上,美国地质调查局这样的实验室会有一群主要由本科生或实习生组成的团队,他们负责查看地震图。"

然而,手动能够发现和分类的地震数量有限。开发有效查找和处理地震的算法长期以来一直是该领域的优先事项——特别是自20世纪50年代初计算机出现以来。

"地震学的历史上,该领域的发展总是随着计算能力的进步而进步,"Bradley告诉我。

然而,传统算法存在一个大问题:它们不容易找到较小的地震,尤其是在嘈杂的环境中。

常见事件的复合地震图

如上图所示,许多不同事件都可能引起地震信号。如果方法过于敏感,可能会错误地将事件检测为地震。这个问题在城市中尤其严重,那里持续不断的交通和建筑噪音可能会淹没小地震。

然而,地震具有特征性的"形状"。上图中的7.7级地震与直升机着陆的波形明显不同。

因此,科学家们的一个想法是从人工标记的数据集中创建模板。如果新波形与现有模板高度相关,则几乎可以确定它是地震。

如果你有足够多的人工标记示例,模板匹配效果非常好。2019年,加州理工学院Zach Ross的实验室使用模板匹配在南加州发现了比先前已知的多10倍的地震,包括本文开头提到的地震。他们发现的所有新地震中,几乎全部都是微小的地震,震级在1级以下。

然而,如果你没有大量预先存在的模板数据集,就难以应用模板匹配。这在南加州不是问题——那里已经基本上完整记录了震级1.7级以上的地震——但在其他地方则是一个挑战。

此外,模板匹配计算成本很高。使用模板匹配创建南加州地震数据集需要200个Nvidia P100 GPU连续运行数天。

必须要有更好的方法。

解析地震Transformer模型

AI探测模型解决了所有这些问题:

  • 它们比模板匹配更快。
  • 因为AI探测模型非常小(约35万个参数,与GPT4.0等大语言模型的数十亿参数相比),它们可以在消费级CPU上运行。
  • AI模型对原始数据集中未代表的区域具有良好的泛化能力。

此外,AI模型还能提供关于不同类型地震震动到达时间的更好信息。对两种最重要的波——P波和S波的到达时间进行计时称为_相位拾取_。它使科学家能够推断地震的结构。AI模型可以在地震探测的同时完成这项工作。

地震探测(和相位拾取)的基本任务如下所示:

地震Transformer模型输出示例

前三行代表不同的振动方向(东西、南北和上下)。给定这三个维度的振动,我们能否确定是否发生了地震,以及何时开始?

我们想要检测初始的P波,它直接从地震源到达。但这可能很棘手,因为P波的回声可能会从其他岩层反射并稍后到达,使波形变得更加复杂。

因此,理想情况下,我们的模型在每个时间步输出三个信息:

  1. 该时刻发生地震的概率。
  2. 第一条P波在该时刻到达的概率。
  3. 第一条S波在该时刻到达的概率。

我们在第四行看到所有三个输出:绿色的检测结果,蓝色的P波到达时间和红色的S波到达时间。(此样本中有两次地震。)

为了训练AI模型,科学家们获取大量标记数据(如上图所示),并进行监督训练。我将描述一个最常用的模型:地震Transformer,它由斯坦福大学团队开发,由S. Mostafa Mousavi领导,他后来成为哈佛大学教授。

与许多地震探测模型一样,地震Transformer借鉴了图像分类的思想。读者可能熟悉AlexNet,这是一个著名的图像识别模型,在2012年掀起了深度学习热潮

AlexNet使用了卷积,这是一种基于物理上接近的像素更可能相关的思想的神经网络架构。AlexNet的第一层卷积将图像分解成小块——每边11个像素——并根据边缘或梯度等简单特征的存在对每个块进行分类。

下一层将第一层的分类作为输入,并检查更高级别的概念,如纹理或简单形状。

每个卷积层分析图像的更大部分,并在更高级别的抽象上操作。在最后的层中,网络查看整个图像并识别"蘑菇"和"集装箱船"等对象。

图像是二维的,因此AlexNet基于二维卷积。相比之下,地震图数据是一维的,因此地震Transformer使用沿时间维度的一维卷积。第一层分析0.1秒时间段的振动数据,而后续层则在 progressively 更长的时间段内识别模式。

很难说地震模型具体提取了什么模式,但我们可以将其类比为使用一维卷积的假设音频转录模型。该模型可能首先识别辅音,然后是音节,然后是单词,然后在逐渐增大的时间尺度上识别句子。

地震Transformer将原始波形数据转换为一系列高级表示,这些表示表明地震和其他地震学重要事件的可能性。随后是一系列反卷积层,它们精确定位地震及其至关重要的P波和S波发生的确切时间。

该模型还在模型中间使用了一个注意力层,用于混合时间序列不同部分的信息。注意力机制在大语言模型中最为著名,它有助于在单词之间传递信息。它在地震图检测中扮演类似角色。地震地震图具有一般结构:P波后跟S波,然后是其他类型的震动。因此,如果某个片段看起来像P波的开始,注意力机制有助于它检查它是否符合更广泛的地震模式。

扩展地震数据规模

地震Transformer的所有组件都来自神经网络文献的标准设计。其他成功的探测模型,如PhaseNet,甚至更简单。PhaseNet仅使用一维卷积来拾取地震波的到达时间。没有注意力层。

根据Byrnes的说法,一般来说,"地震学领域并没有太多需要发明新架构的必要"。源自图像处理的技术已经足够。

那么,是什么使这些通用架构如此有效?数据。大量的数据。

Ars之前曾报道过ImageNet(一个图像识别基准)的引入如何帮助引发深度学习热潮。大型公开可用的地震数据集在地震学中发挥了类似作用。

地震Transformer使用斯坦福地震数据集(STEAD)进行训练,该数据集包含来自全球的120万个人类标记的地震图数据片段。(STEAD的论文明确提到了ImageNet作为灵感来源)。其他模型,如PhaseNet,也在数十万或数百万个标记片段上进行了训练。

斯坦福地震数据集中的所有记录地震

数据和架构的组合效果很好。根据Byrnes的说法,当前的模型在识别和分类地震方面"好得令人难以置信"。通常,机器学习方法发现的地震数量是先前在特定区域识别数量的10倍或更多。你可以在意大利的地震目录中直接看到这一点:

AI与传统方法发现的地震数量对比

AI工具不一定比模板匹配检测到更多地震。但基于AI的技术计算和劳动密集程度低得多,使普通研究项目更容易使用,也更易于在世界各地应用。

总而言之,这些机器学习模型在探测和相位拾取地震方面表现得非常好,几乎完全取代了传统方法,特别是对于较小震级的地震。

AI技术能做什么

地震科学的圣杯是地震预测。例如,科学家们知道西雅图附近将发生大地震,但几乎无法知道它会在明天还是一百年后发生。如果我们能够足够精确地预测地震,以便让受影响地区的人员撤离,那将非常有帮助。

你可能认为AI工具会有助于地震预测,但似乎尚未实现。

康奈尔大学的Judith Hubbard表示,这些应用更加技术化,不那么引人注目。

更好的AI模型为地震学家提供了更全面的地震目录,这解锁了"许多不同的技术",Bradley说。

最酷的应用之一是理解和成像火山。火山活动产生大量小地震,其位置有助于科学家理解岩浆系统的结构。在2022年的论文中,John Wilding及其合作者使用大型AI生成的地震目录创建了夏威夷火山系统的惊人图像。

每个点代表一个单独的地震

他们提供了之前假设的深部Pāhala岩浆复合体与茂纳洛亚浅部火山结构之间岩浆连接的直接证据。你可以在图像中看到这个用箭头标记为Pāhala-Mauna Loa地震活动带的地方。作者还能够将Pāhala岩浆复合体的结构细分为离散的岩浆片。这种级别的细节可能促进更好的地震实时监测和更准确的火山喷发预测。

另一个有前景的领域是降低处理大型数据集的成本。分布式声学传感(DAS)是一种强大的技术,它使用光纤电缆测量电缆全长范围内的地震活动。根据休斯顿大学教授Jiaxuan Li的说法,单个DAS阵列每天可以产生"数百GB的数据"。这些数据可以产生极高分辨率的数据集——足以分辨出单个脚步声。

AI工具使得在DAS数据中非常精确地计时地震成为可能。在DAS数据中引入相位拾取的AI技术之前,Li和他的一些合作者尝试使用传统技术。虽然这些方法"大致可行",但对于他们的下游分析来说不够准确。没有AI,他的许多工作将会"困难得多",他告诉我。

Li也乐观地认为,AI工具未来将能够帮助他在丰富的DAS数据中分离"新型信号"。

并非所有AI技术都取得了成功

与许多其他科学领域一样,地震学家面临着采用AI方法的一些压力,无论这些方法是否与他们的研究相关。

"学校希望你把AI这个词放在一切前面,"Byrnes说,"这有点失控了。"

这可能导致技术上合理但实际上无用的论文。Hubbard和Bradley告诉我,他们看到了很多基于AI技术的论文,这些论文"揭示了对地震工作原理的根本误解"。

他们指出,研究生可能感到压力,要专门学习AI方法,而牺牲了对科学领域基础知识的学习。他们担心,如果这种AI驱动的研究变得根深蒂固,较老的方法将被"无意义地淘汰出局"。

虽然这些问题是真实存在的,并且Understanding AI之前也曾报道过,但我认为它们并不 detract from AI地震探测的成功。在过去的五年中,基于AI的工作流程几乎完全取代了地震学中的一项基本任务,并且取得了更好的效果。

这相当酷。