全面阐述支撑 EarthquakeLive.org 的地球物理目录、Uber H3 空间离散化、地壳应力物理公式,以及经过校准的集成机器学习架构。
分层六边形空间索引,最大限度减少平面与距离畸变
传统的经纬度矩形网格存在严重的极区畸变和不一致的邻域拓扑。EarthquakeLive.org 采用基于正六边形的 Uber H3 离散全球网格系统。
网格直径:约 110 公里
面积:约 11000 平方公里
区域构造带尺度
网格直径:约 40 公里
面积:约 2500 平方公里
局部断层分段尺度
每个六边形聚合空间地震计数、断层距离、b 值斜率和潮汐参数。
刻画应力积累与断层动力学的多维地球物理指标
在正常平衡状态下 b ≈ 1.0(每发生 1 次 M=5 地震,约对应 10 次 M=4 和 100 次 M=3 地震)。当凹凸体闭锁区剪应力不断积累时,微破裂被抑制,b 值降至 0.8 以下,可作为危险性升高的经验代理指标。
由于地震能量随震级变化,弹性回跳位移与 √E 成正比。监测累积应变可以追踪地震平静与前兆性加速矩释放(AMR)。
7 天短期地震频率与 30 天基准背景频率之比。数值升高意味着震群触发或前震级联序列。
震级-频度斜率偏离、弹性应变积累、前震震群加速,以及闭锁断层段上的长期地壳平静。
到 GEM 全球活动断层构造的大地线距离、历史构造滑动速率,以及独热编码的断层运动学类型(俯冲、走滑、逆冲)。
基于天文历建模日月引潮势、朔望周期(大潮与小潮)以及月球轨道近地点相位。
震源深度剖面(浅部地壳与深部贝尼奥夫-和达带)结合 NOAA SWPC 行星际地磁扰动指数(Kp)。
震源较浅(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))计算近岸浅化波高放大效应。
岩浆穿过脆性地壳上升时会使岩石破裂,在火山喷发前产生密集的浅源微震震群。EarthquakeLive 将每个地震事件与史密森尼全球火山计划(GVP)包含 1,214 座全新世火山口的数据库进行匹配(火山口 25 公里半径内,震源深度 ≤ 20 公里),并与 30 天基线地震发生率进行比对,以从构造余震中精准识别真正的岩浆活动。
火山灰云会因硅酸盐玻璃熔融和涡轮叶片磨蚀,对喷气式飞机发动机造成致命灾难性威胁。EarthquakeLive 持续采集全球所有 9 个 ICAO 运营的火山灰咨询中心(VAAC)发布的通告(VAA/VAG),解析涵盖垂直飞行高度层(地表至 FL600+ / 60,000 英尺)在 0h、+6h、+12h 及 +18h 预测周期的空间扩散多边形。
覆盖 12 小时、24 小时、72 小时和 7 天预报时段的经校准集成决策模型
即时的急性前兆与前震危险窗口
短期余震与前震危险性
民防预警的主要业务时间窗
每周构造应变积累的天气图式趋势
| 组件 | 算法 | 预测目标 | 损失函数 / 优化 |
|---|---|---|---|
| 危险性分类器 | 梯度提升决策森林 | 二分类事件:时段内 M ≥ 4.5 | 类别权重平衡,logloss + Brier 损失 |
| 震级回归器 | 稳健的预期震级回归器 | H3 网格内的预期 M_max | Huber 损失(抗离群值) |
| 概率校准 | 单调经验校准器 | P_calibrated ∈ [0, 1] | 最小化 Brier 得分,匹配可靠性曲线 |
在类别高度不平衡的场景下,梯度提升模型的原始输出往往校准不良。保序概率校准保证了经验可靠性:当模型给出 40% 的危险性评分时,历史上每 100 次中恰有 40 次发生地震。
用博弈论方法将风险预测归因到各个具体的物理驱动因子
EarthquakeLive.org 为每个六边形网格计算 TreeSHAP 归因向量。用户与研究人员可点击地图上任意网格,查看每个物理参数的精确百分比贡献。
模型沿新记录的目录持续前推,将预报窗口与实况比对以更新权重,并跟踪特征漂移(PSI)。
0.864
0.078
在时间外推的滚动前推协议下,持续自动地与真实地震目录进行比对评估。
EarthquakeLive 发布的所有预报都带有永久时间戳,并在预报时段结束后与权威目录(USGS ComCat 与 EMSC)比对校验。保序概率校准保证:若某区域被赋予 60% 的概率,则在 100 个此类预报窗口中约有 60 个确实发生地震。
Deck.gl v9 实例化几何渲染、原生 MapLibre 矢量回退、DPR 1.0 显存限制以及穿过反子午线的球面摄像机安全轨迹
对于高性能桌面和移动浏览器,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 度剧烈反向旋转的问题,航线位移通过归一化角位移算法计算:
镜头轨迹具备平滑的三次方缓动、动态俯仰倾角以及专门适配公共大屏展示的自动定格观察周期。
当 18 个实时地震目录中的任何一个录入强震(M ≥ 5.5)时,自动漫游和屏保循环将被实时打断抢占。摄像机镜头平滑直击新震中,弹出关键遥测信息卡,并触发多音频合成警报。
植根于国际地震学文献的严谨理论基础
"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),用于识别地震平静或加速矩释放。
"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."
在系统中的作用: 支撑经过校准的对数正态死亡率曲线与建筑易损性评估曲线。
"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)的衰减规律。
"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 值,确保完全透明。
"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."
在系统中的作用: 计算日月引潮势与朔望对齐周期。
"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.
"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.
"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),用于支持应急准备。
经过验证的科学、运行与应急规程
切换到交互式 WebGL 地图,查看实时断层线、最近地震和经过校准的危险性网格。