让每一次判断，保留它成立的条件。
余震之间是面向历史地震目录的独立分析工具。它把数据审计、参数实验、空间观察和证据复验接成一条可回查的工作流程。
研究问题：描述何时改变方向？
“震后第三周的记录更少”是一项依赖条件的目录描述。震级门槛、空间半径、窗口位置与震级类型会影响所选事件。工具关注这些分析选择如何影响描述，以及第三方能否根据同一快照重建完整结果。
目前研究定位是可复现的地震目录分析与科学软件工程。后续可检验的问题包括：单一默认视图是否掩盖方向改变；记录稀少时，倍率是否产生过度解释；同时展示反例和条件是否改善使用者对证据边界的判断。前两项可通过目录实验研究，最后一项需要真实用户研究，不能由功能测试代替。
实验方法：预先说明参数，完整报告所有格。
| 要素 | 实现与口径 |
|---|---|
| 分析对象 | 六例 USGS ComCat 冻结目录，共 4,145 条原始记录。展示震后 [0, 28) 天，排除主震；未逐条判定区域事件的余震归属。 |
| 时间定义 | 一天为主震时刻起连续 24 小时，采用 UTC，窗口左闭右开。A、B 等长；共同平移保留两窗间隔，不自动裁切越界窗口。 |
| 默认协议 | A [7, 14)、B [14, 21) 天；门槛 3 / 3.5 / 4 / 4.5 / 5；半径 100 / 200 / 300 km；两窗共同平移 −1 / 0 / +1 天。每例 45 格，六例共 270 格。 |
| 逐格结果 | A、B 计数与条/天、B−A、B/A、共有事件、去重并集及事件 ID。A=0 时倍率不定义；双零独立标识，结果中不删格。 |
| 阅读门槛 | 默认每窗至少 5 条，可调整。分别报告所有格与达到门槛的格的分母。它是阅读规则，不是显著性、统计效能或完备性标准。 |
| 空间定义 | 筛选使用半径 6371.0088 km 的球面大圆距离；质量审计根据经纬度重算。显示地形与统计筛选分开，地图不作为烈度或实测断层图。 |
不同震级类型保持原值，没有同尺度化。Mw 家族筛选合并名称以 mw 开头的报告类型，仍不意味着这些数值已校准。不同地区收录能力不同，不能依据计数或“方向一致格数”对地区危险性排序。
每个参数格复用同一目录中的事件，不是独立重复。这里不输出置信度、发生概率或显著性 p 值；空记录也不等于没有地震发生。
证据：从原始字节到逐格事件 ID。
- 读取随包原始 CSV，核对下载清单中的 SHA-256。命令行进一步逐字段核对 CSV 与派生目录。
- 检查重复 ID、非法筛选字段、主震参数、UTC 相对日、球面距离与查询时段。审计失败时阻止实验。
- 完整运行所选参数空间，记录每格的条件、计数和两个窗口的事件 ID；所有摘要从这些格重算。
- 对科学数据快照和实验包分别计算摘要。导入后重新运行协议，逐项核对来源、审计、全部结果和事件 ID。
SHA-256 可以检测内容变化。它不是第三方签名，也不证明作者观察文字正确。重新计算依赖本地冻结目录，无法据此证明外部机构始终保留相同目录版本。
node scripts/audit-data.mjs node scripts/reproduce.mjs --output work/reproduction node scripts/reproduce.mjs --verify work/reproduction/experiment.json
验证状态与后续研究
| 层面 | 已有能力 | 下一步需要的证据 |
|---|---|---|
| 数据与计算 | 固定快照、CSV 字段核对、边界条件测试、实验包重算与篡改拒绝。 | 独立领域人员复核数据适用性；用独立实现与更多目录做复现。 |
| 软件交付 | 本机项目管理、批量参数实验、后台计算、命令行复现与独立静态部署。 | 更广浏览器与设备覆盖、长期部署运行监测、外部使用反馈。 |
| 科学表达 | 计数、倍率、空记录、混合震级、重叠窗口及审计状态明确显示。 | 预先约定使用者任务、评价指标和分析方法，再开展真实用户研究。 |
| 地震研究 | 历史区域目录的描述性探索与证据追溯。 | 依具体科学问题开展完备性分析、震级统一、余震识别和模型检验；当前没有这些研究成果。 |
工程验证不等于地震学结论成立，也不等于基金评审或资助。项目没有声明国家自然科学基金立项、机构背书、专家认可或已验证的教学提升。
2.1 的连续研究流程
研究预设只选择明确的参数空间，不改变计数算法。逐格查验从本地冻结目录解析事件 ID，并保留 A、B 归属。界面筛选仅影响查阅，逐格 CSV 始终导出该格完整并集、条件与来源摘要。实验恢复副本先按原快照复验，2.0 实验包继续兼容。取消或失败不会用迟到的结果覆盖已有实验。
资料与实现贡献
改进集中在数据口径一致的参数实验、失败时阻止继续的审计、逐格事件绑定及跨浏览器和命令行复验。目录、球面距离、SHA-256、三维渲染与常规计数均使用已有技术，不宣称为原创地震预测算法。AI 辅助参与代码与文档实现，保留披露。
USGS 目录接口与参数 · ComCat 文档 · 数据、字体、三维与 AI 资源说明
保留科学演示的来源
原有门槛回放由脚本读取同一冻结目录逐帧计数；真实高程、报告深度、均匀介质 P/S 波模型与 AI 概念影像各自标识。它们承担空间理解或教学解释，不代替观测与统计证据。许可待确认的区域构造原图未随 2.0 部署包提供。
视觉、地形与台站工具
首页三幅地形及井下观测图为 AI 概念视觉，不表示实测地质结构、站点或反演结果。PNG 原件保持 1672 × 941；画面采用裁切视口，原始字节未缩小或重编码。艺术图与三维真实高程分别使用独立入口。
真实高程来自 Mapzen / USGS / NOAA / Copernicus；目录点计数、窗口与来源均在工作台或导出文件中保留。西藏巴青已有真实地形，尚未接入批量目录。
台站工作台支持本地台站与事件 CSV、SAC v6/v7 等间隔时序、人工标记、窗内测量及交接。参考到时采用直线路径与显式 Vp/Vs，不是走时表定位、预警或自动震相识别。未去仪器响应，幅值按原始单位解释。
图像原件与 SHA-256 · 数据和资源许可 · SAC 格式定义
3.0 台站工作台和区域对照新增西藏定日 2025 年目录：USGS 80 条原始区域记录。包含既有六例时共 7 例、4,225 条；旧目录研究入口继续保留六例默认 270 格。中国地震台网报告 M6.8 与 USGS 报告 M7.1 分别保留，不作震级同尺度替换。
波形科研与 P/S 分析
科研工作台使用本机 ObsPy、SciPy、PyWavelets 与随包 PhaseNet ONNX 模型，支持 SAC、MiniSEED、REFTEK130 和 StationXML。输入单位、分量映射、校准、频带、窗口和来源摘要写入每次研究记录。
噪声、P、S 均使用半开分窗；RMS 信噪比与峰值能量信噪比分别报告。当前 PDF 的样本是完整重叠噪声块，不是长期小时 PPSD。小波色阶为系数功率，不是 PSD。PhaseNet 返回待复核候选；STA/LTA 是触发，TauP 是一维模型理论走时，均不替代人工观测。
导出 SVG、高清 PNG、测量/候选 CSV 和带摘要 JSON。参数变更后必须重新计算。三张宣传概念图不参与科研计算。