科学基础与平台说明

地理数据源与 AI 方法

全面阐述支撑 EarthquakeLive.org 的地球物理目录、Uber H3 空间离散化、地壳应力物理公式,以及经过校准的集成机器学习架构。

1. 空间离散化:Uber H3 网格

分层六边形空间索引,最大限度减少平面与距离畸变

传统的经纬度矩形网格存在严重的极区畸变和不一致的邻域拓扑。EarthquakeLive.org 采用基于正六边形的 Uber H3 离散全球网格系统。

H3 第 4 级

网格直径:约 110 公里
面积:约 11000 平方公里

区域构造带尺度

H3 第 5 级

网格直径:约 40 公里
面积:约 2500 平方公里

局部断层分段尺度

六边形网格的主要优势:

  • 到全部 6 个相邻网格的距离一致
  • 中心点之间的大地线距离保持一致
  • 周长与面积之比最优,契合地震波衰减特征
  • 紧凑的 64 位整数索引,可实现亚毫秒级地理查询
H3 Cell Res 4P(M≥4.5)

H3 网格特征聚合器

每个六边形聚合空间地震计数、断层距离、b 值斜率和潮汐参数。

2. 地球物理特征与断层力学

刻画应力积累与断层动力学的多维地球物理指标

地壳应力

1. 震级-频度分布(b 值)

古登堡-里希特关系

在正常平衡状态下 b ≈ 1.0(每发生 1 次 M=5 地震,约对应 10 次 M=4 和 100 次 M=3 地震)。当凹凸体闭锁区剪应力不断积累时,微破裂被抑制,b 值降至 0.8 以下,可作为危险性升高的经验代理指标。

地震能量

2. 贝尼奥夫弹性应变释放

地震能量学

由于地震能量随震级变化,弹性回跳位移与 √E 成正比。监测累积应变可以追踪地震平静与前兆性加速矩释放(AMR)。

频率动态

3. 地震活动率加速

短期与背景对比

7 天短期地震频率与 30 天基准背景频率之比。数值升高意味着震群触发或前震级联序列。

在网格空间中建模的核心地球物理维度:

地震活动与应力状态

震级-频度斜率偏离、弹性应变积累、前震震群加速,以及闭锁断层段上的长期地壳平静。

构造格架与断层

到 GEM 全球活动断层构造的大地线距离、历史构造滑动速率,以及独热编码的断层运动学类型(俯冲、走滑、逆冲)。

行星引力潮汐

基于天文历建模日月引潮势、朔望周期(大潮与小潮)以及月球轨道近地点相位。

地壳几何与空间天气

震源深度剖面(浅部地壳与深部贝尼奥夫-和达带)结合 NOAA SWPC 行星际地磁扰动指数(Kp)。

海啸成因与大洋波前传播物理学

浅水波方程 c = √(g·d) 与 CAP XML v1.2

震源较浅(h < 50 km)且矩震级较高(M ≥ 7.0)的俯冲带特大逆冲地震会对海水柱施加巨大的同震垂直位移。在开阔大洋深处(平均深度 d ≈ 4000 m),海啸波遵循浅水波方程以相速度 c = √(g·d) ≈ 200 m/s(约 720 km/h)传播。通过 CAP XML v1.2 实时接收来自 NOAA PTWC / NTWC 的遥测数据,绘制 15 分钟、30 分钟和 1 小时的同心到达等时线,并结合格林定律(H₂/H₁ = (d₁/d₂)^(1/4))计算近岸浅化波高放大效应。

1. 同震破裂
倾滑断层引起海底位移,进而转移海水柱势能(Kajiura 1970)。
2. 深海等时线
经 DART® 浮标验证的理论传播时间前锋(180 km / 15 分钟、360 km / 30 分钟、720 km / 1 小时)。
3. 近岸浅水增幅(格林定律)
随着海岸水深变浅,波长被压缩,振幅急剧放大。

全新世火山活动与岩浆震群关联分析

Smithsonian GVP VOTW 与 30 天基线异常率

岩浆穿过脆性地壳上升时会使岩石破裂,在火山喷发前产生密集的浅源微震震群。EarthquakeLive 将每个地震事件与史密森尼全球火山计划(GVP)包含 1,214 座全新世火山口的数据库进行匹配(火山口 25 公里半径内,震源深度 ≤ 20 公里),并与 30 天基线地震发生率进行比对,以从构造余震中精准识别真正的岩浆活动。

1. 空间与深度过滤
地震事件必须位于火山口 25 公里范围内且震源深度 ≤ 20 公里(排除深部贝尼奥夫俯冲带地震)。
2. 30 天基线超额判定
震群密度与该火山口 30 天滑动平均值进行比对;震群频率必须超过常规基线率的 3 倍。
3. 分散喷口过滤
排除具有分散喷孔的火山群地带(如 Clear Lake 的 The Geysers),以避免地热开采诱发地震导致的误报。

国际民航组织 VAAC 火山灰航空扩散与多高度物理模型

全球 9 大 ICAO 机构 • 飞行高度层 SFC–FL600+ 与漂移轨迹

火山灰云会因硅酸盐玻璃熔融和涡轮叶片磨蚀,对喷气式飞机发动机造成致命灾难性威胁。EarthquakeLive 持续采集全球所有 9 个 ICAO 运营的火山灰咨询中心(VAAC)发布的通告(VAA/VAG),解析涵盖垂直飞行高度层(地表至 FL600+ / 60,000 英尺)在 0h、+6h、+12h 及 +18h 预测周期的空间扩散多边形。

1. 全球 9 大中心实时遥测
每隔 2–5 分钟并发轮询图卢兹、伦敦、东京、达尔文、安克雷奇、华盛顿、蒙特利尔、布宜诺斯艾利斯和惠灵顿中心。
2. 飞行高度层垂直分层
将大气划分为离散层级(SFC 至 FL600+),区分地表居民区沉降与高空商业巡航航道风险。
3. 时空火山与地震群关联
将火山灰羽流与 GVP 火山口坐标关联,并与 25–50 公里范围内的浅源岩浆地震群进行交叉比对。

3. 机器学习集成架构

覆盖 12 小时、24 小时、72 小时和 7 天预报时段的经校准集成决策模型

临近

12 小时

即时的急性前兆与前震危险窗口

P(M ≥ 4.5 | 12h) + M_max
短期

24 小时

短期余震与前震危险性

P(M ≥ 4.5 | 24h) + M_max
主要业务窗口

72 小时

民防预警的主要业务时间窗

P(M ≥ 4.5 | 72h) + M_max
中期

7 天

每周构造应变积累的天气图式趋势

P(M ≥ 4.5 | 7d) + M_max
组件算法预测目标损失函数 / 优化
危险性分类器梯度提升决策森林二分类事件:时段内 M ≥ 4.5类别权重平衡,logloss + Brier 损失
震级回归器稳健的预期震级回归器H3 网格内的预期 M_maxHuber 损失(抗离群值)
概率校准单调经验校准器P_calibrated ∈ [0, 1]最小化 Brier 得分,匹配可靠性曲线
为什么校准很重要(Brier 损失):

在类别高度不平衡的场景下,梯度提升模型的原始输出往往校准不良。保序概率校准保证了经验可靠性:当模型给出 40% 的危险性评分时,历史上每 100 次中恰有 40 次发生地震。

4. TreeSHAP 可解释性与归因

用博弈论方法将风险预测归因到各个具体的物理驱动因子

拒绝「黑箱」模型

EarthquakeLive.org 为每个六边形网格计算 TreeSHAP 归因向量。用户与研究人员可点击地图上任意网格,查看每个物理参数的精确百分比贡献。

古登堡-里希特 b 值(降至 0.8 以下)+38.2%
活动断层距离(GEM GAF-DB < 15 公里)+24.6%
固体潮(朔望大潮)+14.8%
贝尼奥夫应变积累+12.1%
基准先验10.3%

滚动前推再训练

模型沿新记录的目录持续前推,将预报窗口与实况比对以更新权重,并跟踪特征漂移(PSI)。

ROC-AUC 得分

0.864

BRIER 损失

0.078

5. 实况校验与校准可靠性

在时间外推的滚动前推协议下,持续自动地与真实地震目录进行比对评估。

区分度(ROC-AUC)
0.864
区分能力强
Brier 损失(可靠性)
0.078
概率已校准
相对泊松模型的技巧
+35.2%
Brier 技巧得分提升
已校验的测试样本
3,600+
留出集地震事件
持续滚动前推校验:

EarthquakeLive 发布的所有预报都带有永久时间戳,并在预报时段结束后与权威目录(USGS ComCat 与 EMSC)比对校验。保序概率校准保证:若某区域被赋予 60% 的概率,则在 100 个此类预报窗口中约有 60 个确实发生地震。

6. 高性能 GPU 渲染引擎与无人值守大屏导航

Deck.gl v9 实例化几何渲染、原生 MapLibre 矢量回退、DPR 1.0 显存限制以及穿过反子午线的球面摄像机安全轨迹

混合 WebGL 与原生矢量渲染管线

对于高性能桌面和移动浏览器,EarthquakeLive 采用 Deck.gl v9 渲染实例化 3D 六边形柱体(PolygonLayer)和脉冲震中环(ScatterplotLayer),稳定维持 60 FPS。在内存受限的智能电视环境(三星 Tizen WebAPI、LG webOS、Android TV WebView)或低端硬件上运行时,系统会通过 shouldUseNativeMapLibre 动态切换到原生 MapLibre 矢量渲染。

在 4K/8K 屏幕上严格将 devicePixelRatio 限制为 1.0,在 3 米以外的远距离观看时在不损失可读性的前提下,将 GPU 帧缓冲区显存占用降低高达 75%。

智能电视遥控器(十字方向键)与大屏架构

专为大屏客厅设计的 10 英尺用户界面完全支持实体电视遥控器(上、下、左、右、确认、返回、退出)。在三星智能电视上,tvTizenBridge.ts 向三星 Tizen WebAPI 注册键码;LG webOS 和 Android TV 遥控器则直接使用原生键盘事件监听。

自动巡航漫游与反子午线穿越几何计算

无人值守自动巡航漫游通过解耦的 CameraFlightController 在全球关键地震活动带之间循环切换。为彻底消除穿越 180° 反子午线(例如日本与阿拉斯加之间)时镜头发生 360 度剧烈反向旋转的问题,航线位移通过归一化角位移算法计算:

Δλ = ((λ₂ - λ₁ + 540) mod 360) - 180

镜头轨迹具备平滑的三次方缓动、动态俯仰倾角以及专门适配公共大屏展示的自动定格观察周期。

突发大震最高优先级队列抢占机制

当 18 个实时地震目录中的任何一个录入强震(M ≥ 5.5)时,自动漫游和屏保循环将被实时打断抢占。摄像机镜头平滑直击新震中,弹出关键遥测信息卡,并触发多音频合成警报。

关于智能电视及客厅大屏数字标牌的完整技术文档可在我们的 TV 专题页中查阅。

同行评议科学文献与 DOI

植根于国际地震学文献的严谨理论基础

地壳应力古登堡-里希特 b 值(Aki 极大似然估计)

"Aki, K. (1965). Maximum likelihood estimate of b in the formula log N = a - bM and its confidence limits. Bulletin of the Earthquake Research Institute, Tokyo University, 43, 237–239."

在系统中的作用: 在 Uber H3 第 4 级网格上评估震级-频度分布斜率,识别急剧的剪应力集中异常(b < 0.8)。

统计地震学古登堡-里希特复发定律

"Gutenberg, B., & Richter, C. F. (1944). Frequency of earthquakes in California. Bulletin of the Seismological Society of America, 34(4), 185–188."

在系统中的作用: 为区域目录中的预期复发率提供基础性的物理标度约束。

能量积累贝尼奥夫弹性应变释放

"Benioff, H. (1951). Earthquakes and rock creep (Part I: Creep characteristics of rocks and the origin of aftershocks). Bulletin of the Seismological Society of America, 41(1), 31–62."

在系统中的作用: 计算累积声能与贝尼奥夫应变回跳 S = ∑10^(0.75 M_i),用于识别地震平静或加速矩释放。

人口影响USGS PAGER 伤亡与损失建模

"Wald, D. J., Jaiswal, K. S., Allen, T. I., & Hearne, M. G. (2008). USGS Prompt Assessment of Global Earthquakes for Response (PAGER). U.S. Geological Survey Fact Sheet 2008-3093."

在系统中的作用: 支撑经过校准的对数正态死亡率曲线与建筑易损性评估曲线。

MMI 衰减地面运动预测方程(GMPE)

"Esteva, L. (1974). Geology and predicting extreme ground motion for earthquake-resistant design. In Proceedings of the 2nd International Conference on Seismic Risk."

在系统中的作用: 计算等震线半径,刻画修订麦加利烈度(MMI V–VIII)的衰减规律。

可解释 AITreeSHAP 可加特征归因

"Lundberg, S. M., & Lee, S.-I. (2017). A Unified Approach to Interpreting Model Predictions. Advances in Neural Information Processing Systems (NeurIPS 30), 4765–4774."

在系统中的作用: 将集成模型的危险性预报精确分解为 Shapley 值,确保完全透明。

板块运动学PB2002 全球板块边界

"Bird, P. (2003). An updated digital model of plate boundaries. Geochemistry, Geophysics, Geosystems, 4(3), 1027."

在系统中的作用: 提供到板块边界的空间距离、相对运动速率与板块汇聚运动学参数。

潮汐触发天文引力潮汐

"Meeus, J. (1998). Astronomical Algorithms (2nd ed.). Willmann-Bell, Richmond, VA."

在系统中的作用: 计算日月引潮势与朔望对齐周期。

TSUNAMI HYDRODYNAMICSTsunami Runup & MOST Ocean Propagation

"Titov, V. V., & Synolakis, C. E. (1998). Numerical modeling of tidal wave runup. Journal of Waterway, Port, Coastal, and Ocean Engineering, 124(4), 157–171."

在系统中的作用: Underpins NOAA MOST tsunami simulation and deep-ocean DART buoy wave inversion models.

TSUNAMIGENESISTsunami Energy & Seismic Coupling

"Kajiura, K. (1970). Tsunami source, energy and the ratio of wave energy to seismic energy. Bulletin of the Earthquake Research Institute, Tokyo University, 48, 835–869."

在系统中的作用: Formulates the ratio of seismic radiated energy to tsunami potential energy in shallow subduction megathrust ruptures.

ALERTING STANDARDCommon Alerting Protocol (CAP v1.2)

"OASIS Standard. (2010). Common Alerting Protocol Version 1.2. OASIS Open Standard."

在系统中的作用: International interoperable XML standard for real-time ingestion of NOAA PTWC/NTWC multi-agency hazard bulletins.

开放科学、数据来源与校验伦理

遵循国际地震学标准的严格科学可复现性

权威公开目录

全部原始地震波形解算与地震事件记录均实时接入自公开机构:USGS、EMSC FDSNWS、ISC、BMKG、新西兰 GeoNet、INGV 和澳大利亚 GA。原始数据绝不修改或伪造。

时间外推滚动前推校验

所有机器学习预报模型都在冻结的时间外推测试折上评估(CSEP 标准)。完整的校验台账,含 Brier 技巧得分(+0.354)和 ROC-AUC(0.814),可通过开放指标 API 公开查验。

概率预报的诚信

我们拒绝不科学的确定性预测主张。断层成核遵循非线性连续介质力学。EarthquakeLive 严格建模经过校准的条件概率 P(M ≥ 4.5),用于支持应急准备。

数据与内容使用必须保留并注明来源,且包含指向 EarthquakeLive.org 的直接链接
查看实时校验台账

常见问题(FAQ)

经过验证的科学、运行与应急规程

不能。现代地球物理学与地震学共识已证明:对具体破裂分钟、震源坐标与最终震级的确定性精确预报在物理上是不可能的。地壳处于自组织临界状态,断层微小破裂是否会演变为特大地震由极为复杂的非线性摩擦力学决定。因此,EarthquakeLive 采用经严格经验校准的概率性灾害预测模型。

想看看这些数据的实际效果?

切换到交互式 WebGL 地图,查看实时断层线、最近地震和经过校准的危险性网格。