3.3 相似律与药量反演：从现场痕迹推算当量

3.3 相似律与药量反演：从现场痕迹推算当量

立方根相似律指出：药量放大八倍、距离放大两倍，冲击波参数不变——一切超压经验公式都建立在"比例距离"（距离除以药量的立方根）这一单一自变量上。用它加上一组超压-破坏阈值表，就能把现场痕迹翻译成药量区间。本节完成一次完整的反演演练，这是勘测视角的巅峰一课。

本节收束第 3 章：3.1 给了波后状态的算法，3.2 给了反射修正，本节把它们接到破坏阈值数据库上，跑通"痕迹到药量"的全流程。第 8 章 TNT 当量法是它的直接延伸。

立方根相似律为什么成立

回到量纲分析（第 1 章守恒框架的副产品）：空气中爆炸波的自变量只有药量 W（提供能量）、距离 r、环境密度与环境压力。四个量只能组成一个无量纲组合——r 除以 W 的立方根，记作 Z，即比例距离。物理含义：爆炸波在空气里"看到的"不是绝对距离，而是相对于能量球膨胀尺度的距离。药量乘 1000，只要距离乘 10，波形完全相似。这条规律经七十年的核试验与实验室数据反复验证，是爆炸力学里最可靠的经验定律。

工程上最常用的是超压经验式族，例如常用的多项式形式：给定 Z（米每千克的三分之一次方），返回入射超压。正压冲量与持续时间也有对应的经验式。勘测者应该把这些公式当"标准尺"用，但要知道每条的标定范围。

破坏阈值：痕迹的"双语词典"

把现场痕迹换算成超压，靠的是阈值表。常用条目：

⚠️ 常见坑：阈值是概率分布不是门槛。玻璃在 3 千帕开始破碎、7 千帕近全碎，中间是一个随距离展开的"破碎梯度"。精明的勘测者恰恰利用这个梯度：找到"刚好开始破碎"的最远窗户，用下阈值；找"全部粉碎"的最近区域，用上阈值——两个半径给出药量区间的上下限。

完整反演演练

某露天场地爆炸事故，无药量台账。现场勘测记录：距疑似爆心 40 米处轻质隔墙倒塌；70 米处玻璃开始零星破碎；120 米处玻璃完好；25 米处混凝土地坪剥落。目标：给出 TNT 当量区间。

# 多痕迹药量反演（比例距离迭代法）
# 常用超压经验式（多点分段式，此处用简化版便于演示）：
def overpressure(Z):
    # 入射超压 kPa，Z 为 m/kg^(1/3)，适用约 1.2 < Z < 40
    return 0.0717/Z + 0.1775/Z**2 + 0.801/Z**3 - 0.6  # 示意系数
# 迭代：给定半径 r 与阈值 dp，解 W 使 overpressure(r/W**(1/3)) == dp
def solve_W(r_m, dp_kPa):
    lo, hi = 1e-3, 1e4   # 药量上下界 kg
    for _ in range(80):
        mid = (lo*hi)**0.5
        if overpressure(r_m/mid**(1/3)) > dp_kPa: lo = mid
        else: hi = mid
    return (lo*hi)**0.5
# 三条证据：25 m 混凝土剥落（反射增强后约 100 kPa 自由场约 50）
w1 = solve_W(25, 50)
w2 = solve_W(40, 18)    # 隔墙倒塌取中值
w3 = solve_W(70, 4)     # 玻璃下阈
print(f"混凝土地坪证据：W ≈ {w1:.0f} kg TNT")
print(f"隔墙倒塌证据：  W ≈ {w2:.0f} kg TNT")
print(f"玻璃下阈证据：  W ≈ {w3:.0f} kg TNT")
​

运行结果三条证据分别落在数十千克量级且彼此偏差在两倍以内（具体数值随经验式系数而定）——这个一致性正是判据：若某条证据反解出的药量偏离其余证据五倍以上，优先怀疑该痕迹的阈值取值或位置记录，而不是急于平均。最终报告给出区间而非单值，例如"药量介于 40 至 120 千克 TNT 当量，最优估计 80 千克左右"，并附上三条证据的独立反演过程。

反演的误差预算与交叉验证

反演精度受三层误差控制：痕迹定位误差（爆心位置差 10%，距离差 10%，超压差约 25%）、阈值分布宽度（本身有上下限）、经验式标定误差（不同公式族之间可差 20%）。三层叠乘，单条证据的药量误差接近一倍——这就是为什么多证据交叉是强制要求而非加分项。除了超压痕迹，还应引入独立证据链：弹坑尺寸经验式（第 8 章）、地面振动数据、破片分析（第 1 章方法）。三链各自给出区间，重叠部分才是可信结论。

图：反演流程与纪律。右侧的"误差控制纪律"比公式本身更重要——勘测结论的法律与工程分量，全靠这些流程纪律支撑。

💡 关键直觉：相似律的价值不止于计算，还在于"缩比实验合法性"——1 千克炸药在 1 米处的波形，与 1 吨在 10 米处完全相似。第 6 章实验室小装药数据能外推到工程尺度，靠的就是这条定律。

气象与地形修正

相似律的标定数据来自均匀大气与平坦地面，真实事故现场两条都不满足，勘测者要会做两类修正。气象修正：风速梯度使波沿下风向增强、上风向削弱，逆温层则像一面高空镜面把远场波折回地面，形成远离爆心处的"异常破坏岛"——历史上多起大当量事故中，几十公里外的窗户破坏反而比近处零星，都是逆温折射的杰作。地形修正：山谷导向与坡面反射叠加可局部放大两到三倍，开阔斜坡背侧则形成衰减影区。规范做法是把现场破坏痕迹按方位分区，分别反演再比较，方位间系统性差异本身就是气象地形效应的记录，剔除后再合成药量区间。

本节要点回顾

• 立方根相似律：比例距离是冲击波参数的唯一自变量，缩尺实验与全尺寸外推的法理基础；
• 阈值是梯度不是门槛：破碎梯度给出药量区间上下限，梯度本身就是信息；
• 多证据交叉强制：单链误差近一倍，三链重叠区间才是结论；
• 流程纪律高于公式：偏差先查证据、结论写区间、修正显式交代。

第 3 章完结。波在空气里的旅程到此为止，下一章它钻进固体，变成应力波——破坏的形态从"推倒"变成"撕裂"。