地震震中三边定位的物理与几何原理
在地震学中，利用多个地震台站记录到的波形数据来确定震中位置是一项经典方法。由于 P波（初至波/纵波）的传播速度快于 S波（剪切波/横波），因此随着传播距离的增加，这两类地震波到达同一个台站的时间差（即 S−P 走时差）会逐渐增大。
在简化模型中，如果假设地壳介质是均匀的，且地震发生在极浅的地表（浅源地震），那么单一台站观测到的 S−P 走时差 Δ t 就可以直接约束震中到该台站的距离 D。其推导公式如下：
Δ t = D × (1 / Vₛ - 1 / Vₚ)
通过代数变形，可得出距离近似公式：
D = Δ t ÷ (1 / Vₛ - 1 / Vₚ)
在几何上，一个台站提供的 D 仅能确定一个以该台站为圆心、以 D 为半径的距离圆。两个台站的距离圆在球面上通常会相交于两个点。因此，至少需要三个台站才能通过三边交汇确定唯一的震中位置。当引入四个或更多台站时，由于观测误差和模型简化，各圆往往无法完美交汇于一点，此时就需要引入最小二乘法进行数学拟合，以寻找一个使所有台站残差平方和最小的最优解。
传播模型与台站观测输入规范
为了获得准确的拟合结果，用户需要配置合理的传播模型参数并输入 3 至 8 个台站的观测数据。
1. 传播模型参数
- P波速度（
fieldPVelocity）：必须大于 0 且不超过 20 km/s，同时必须大于 S波速度。若输入不合规，系统将提示：Vₚ/Vₛ必须大于 0 且不超过 20 km/s。（其中Vₚ/Vₛ代表具体的输入字段名称）或P波速度必须大于S波速度。 - S波速度（
fieldSVelocity）：必须大于 0 且不超过 20 km/s，同时必须小于 P波速度。 - 走时不确定度 (1σ)（
fieldTimingSigma）：用于评估定位误差，必须大于 0 且不超过 60 秒。若超出范围，将显示：时间不确定性必须大于 0 且不超过 60 秒。
2. 台站观测数据
用户可以手动添加台站（点击“添加台站”），或点击“加载东京示例”快速载入测试数据。若需重置，可点击“清空”。每个台站需提供以下数值：
- 纬度：必须介于 −90° 至 90° 之间。若超出范围，提示：
台站 1：纬度必须介于 −90° 至 90° 之间。（其中1代表具体的台站名称或编号） - 经度：必须介于 −180° 至 180° 之间。若超出范围，提示：
台站 1：经度必须介于 −180° 至 180° 之间。（其中1代表具体的台站名称或编号） - S−P 走时差：必须大于 0 且不超过 2000 秒。若超出范围，提示：
台站 1：请输入大于 0 且不超过 2000 秒的 S−P 走时差。（其中1代表具体的台站名称或编号）
3. 约束规则与异常处理
- 台站数量限制：计算器要求台站数量在 3 到 8 个之间。若不满足，将显示错误：
请使用 3 至 8 个台站。 - 非数值输入：若输入了非数字字符，提示：
input：“abc”不是一个有效的数字。（其中input代表输入出错的界面位置，abc代表用户输入的非数字文本） - 重复坐标：若两个台站坐标完全重合，提示：
台站 1 与 2 的坐标相同。请使用不同的台站位置。（其中1和2代表重合的台站名称） - 几何共线/集群：若所有台站呈一条直线或紧密挤在一点，无法约束二维平面位置，提示：
这些台站的几何分布无法约束一个稳定的二维位置。请在远离当前直线或集群的位置添加一个台站。 - 数值溢出：若计算过程中数值超出计算机处理范围，提示：
数值或中间计算结果超出了支持的数值范围。
最小二乘拟合与误差分析
本工具将地球视为一个半径 R = 6,371 km 的完美球体。
1. 距离与拟合计算
对于每个台站，首先根据 S−P 走时差计算出观测距离（即“S−P 距离”）。接着，通过半正矢公式计算拟合震中到台站的球面距离（即“球面距离”）：
d = R × haversine central angle
震中拟合的数学目标是使所有台站的残差平方和最小化：
震中拟合：使所有台站的 Σ(dᵢ - Dᵢ)² 最小化
2. 残差的定义
残差（columnResidual）定义为拟合球面距离减去由 S−P 走时差推算出的距离。
- 正残差：代表拟合震中到该台站的实际球面距离大于由走时差算出的半径（即拟合点在圆外）。
- 负残差：代表拟合震中到该台站的实际球面距离小于由走时差算出的半径（即拟合点在圆内）。
3. 拟合质量评估
计算器会输出以下指标来评估解算质量：
- 距离拟合质量：根据残差大小分类为“紧密距离拟合”、“交错距离拟合”或“较差距离拟合”。
- 台站几何质量：根据台站相对于震中的空间分布，分类为“强台站几何分布”、“中等台站几何分布”或“弱台站几何分布”。
- 均方根（RMS）距离残差 与 最大距离残差。
教学模型的局限性
本工具是一个旨在展示三边定位和最小二乘法基本原理的教学模型，并非专业的地震目录定位工具。
| 物理因素 | 教学模型（本工具） | 专业地震定位（如 USGS 目录） |
|---|---|---|
| 地球形状 | 完美球体（R = 6,371 km） | 椭球体或局部大地测量坐标系 |
| 速度结构 | 均匀介质、恒定波速 | 随深度变化的分层速度模型或三维速度模型 |
| 震源深度 | 忽略不计（假设为极浅源地震） | 联合求解经度、纬度、深度和发震时刻 |
| 射线路径 | 直线传播（沿球面折算） | 考虑地壳与地幔界面折射、反射的弯曲路径 |
| 误差来源 | 仅考虑用户输入的走时 1σ 不确定度 | 包含台站时钟偏差、震相识别错误、速度模型系统误差等 |
界面输出与公式代入示例
在计算完成后，工具的“公式与数值代入”区域会直观地展示推导步骤：
- 距离因子计算： 距离因子 = Vₚ × Vₛ ÷ (Vₚ - Vₛ) = 每秒 S-P 走时差对应 X km
- 单台距离推导： 台站名称：T 秒 × X 千米/秒 = Y 千米
- 最佳球面拟合：
输出格式如：
最佳球面拟合 = 纬度 A°，经度 B°；均方根（RMS）残差 = C km - 误差外推： 根据输入的走时 1σ 不确定度，结合台站几何分布，外推至水平主轴方向的定位不确定度。
在界面下方，系统会显示方法核对信息。方法与数据源核对日期为 2026年8月6日。状态栏会根据计算状态显示：
修改台站观测数据以重新拟合震中。已基于球面地球模型拟合 stations 个台站。已拟合 stations 个台站，但残差或台站几何分布需要检查。
同时，工具提供“台站几何分布示意图”，并附有说明：“仅显示相对几何关系 —— 无底图、边界或实时地震数据。”
数据隐私与本地处理
本工具充分保护用户的数据隐私。所有输入的台站坐标、S−P 时间以及最终的震中解算结果，均完全在用户的本地浏览器中进行实时计算与处理。任何数据均不会被上传至外部服务器，确保了处理过程的独立性。
常见问题解答 (FAQ)
为什么需要三个台站？
单台的 S−P 测量值只能确定一个距离圆，而无法确定方向。两个圆会在两处相交。通常引入第三个台站即可确定唯一的震中位置，而第四个或更多台站则能通过残差反映出测量偏差。相比于沿直线分布的台站，环绕震源分布的台站能对纬度和经度提供更好的几何约束。
S−P 距离是精确的地表距离吗？
不是。恒定波速公式仅通过到达时间差来估算一条共同的传播路径。将该半径视为球面地表距离，是基于浅源地震和单一速度结构的假设。实际的走时曲线取决于震源深度、震中距离以及每条射线穿过的介质。
正残差或负残差代表什么意思？
残差是指拟合的球面距离减去由 S−P 走时差推算出的距离。正值表示拟合点到该台站的距离大于其 S−P 半径；负值表示距离更近。较大的正负交错残差表明，在所选波速下，各圆无法很好地交汇于同一个清晰的交点。
定位不确定度中包含哪些因素？
它将单一的走时 1σ 不确定度，通过 S−P 距离因子以及拟合点附近的台站几何分布进行外推。该不确定度并未包含速度模型偏差、震源深度、震相识别错误、台站时钟偏差或三维地球结构的影响，因此不能作为完整的精度声明。