灵敏度系数:不是数学题,是 Excel 题

【标题】灵敏度系数：不是数学题，是 Excel 题

【摘要】灵敏度系数那列空着或填1，不是不想填，是真的不会算。TÜV SÜD NEL推荐的数值扰动法：一个公式，三列Excel，任意复杂计量公式都适用。

【正文】

本文整理自 TÜV SÜD National Engineering Laboratory 2026年6月测量不确定度网络研讨会，主讲：Colin Likedy，基于 ISO/IEC Guide 98-1:2024（GUM 最新版）。

见过不少不确定度评估报告，那一列标注为"c"或"Sensitivity Coefficient"的栏，有三种填法：

第一种：空着。

第二种：全填1。

第三种：写备注"按线性关系近似处理"，然后空着。

有时候那列直接消失，好像灵敏度系数是可选项。

去年某原油管道项目 FAT 验收，计量撬供应商交了份不确定度评估，结论 ±0.18%（k=2）。业主代表审完直接打回来：温度和压力测量的灵敏度系数未给出，不确定度贡献值无从核实——不满足 API MPMS Ch.13.3，不能作为系统验收的技术文件。

那份报告的计算本身没有错。缺的只是那一列。

TÜV SÜD National Engineering Laboratory今年6月专门做了一场针对灵敏度分析的网络研讨会，主讲人 Colin Likedy 是 NEL 流量测量咨询部门负责人。下面的方法来自那场 webinar。

偏微分在哪里失效

GB/T 27759-2011（等同 ISO 5168:2005）A.4.2 节给了标准答案：ci = ∂y/∂xi。简单公式没问题——质量 = 体积 × 密度，一行偏导就搞定。

但工程公式不是这样的。

API MPMS Ch.12.2的液体计量链：流量计原始脉冲 → 计量系数 MF → GOV → 扣除 BS&W → CTL 温度修正 → CPL 压力修正 → GSV → NSV → 最终质量。每一步都有不确定度分量，每个分量都需要灵敏度系数。

手算偏微分？不现实。

天然气计量更麻烦。AGA-9 本身有解析不确定度公式，但一旦进了流量计算机——超声计速度测量 + AGA-8 压缩因子迭代 + 温度压力修正 + 能量计算——整条链串在一起，偏微分就变成了不可能完成的任务。

这才是"那列空着"的真实原因：不是懒，是真的不知道怎么算。

一个公式，三列 Excel

方法叫数值扰动法（Numerical Method）。原理不复杂：想知道输入量 x 对输出量 y 的影响有多大，不用推公式，直接动它。

具体做法：

将 x 增加一个小量 δx，算出输出 y⁺
将 x 减少 δx，算出输出 y⁻
灵敏度系数：c = (y⁺ − y⁻) / (2δx)

δx 取该输入量的标准不确定度 u。

Excel 三列：x+δx 时的输出、x−δx 时的输出、两者之差除以 2δx。公式往下拉，所有输入量批量完成。

走一遍储罐体积的例子。

储罐直径 4.8 m，标准不确定度 u(d) = 0.024 m。高度 5.3 m。体积公式 V = πd²h/4，当前体积 95.91 m³。

对直径做扰动：

直径取 4.824 m（+0.024 m），计算体积 = 96.868 m³
直径取 4.776 m（−0.024 m），计算体积 = 94.950 m³

灵敏度系数：c(d) = (96.868 − 94.950) / (2 × 0.024) = 1.918 / 0.048 =39.96 m²

解析法偏微分结果：∂V/∂d = πdh/2 = π × 4.8 × 5.3 / 2 =39.96 m²

完全一致。两种方法同时用，可以互相验证——只要数字吻合，计算就没有错。

对高度做同样操作，得到 c(h) = 18.10 m²。把所有分量的 c×u 算出来，用方和根合成，储罐体积最终结果：

95.91 ± 3.13 m³（95% 置信水平）

两个容易踩的坑

数值扰动法本身不难，但两处容易出问题。

坑一：必须在实际测量值处计算。

灵敏度系数不是常数。

同样的储罐，液位高度从 5.3 m 降到 2.3 m，直径的灵敏度系数就从 39.96 m² 变成 17.34 m²——差了将近一倍。原因是 V = πd²h/4 里，体积对直径是非线性关系，不同液位处的局部斜率不同。

解决方法只有一个：用测量时的实际值，不能用量程中点，不能用标称值。

坑二：绝对不确定度配绝对灵敏度，相对不确定度配相对灵敏度。

如果不确定度预算用的是百分比，就必须用相对灵敏度系数 c*：

c* = c × x / y

储罐例子：c(d) = 39.96 × 4.8 / 95.91 =2，c(h) = 18.10 × 5.3 / 95.91 =1。

直觉上理解：体积是直径的平方关系，直径相对不确定度被放大 2 倍；体积是高度的线性关系，倍数为 1。

混用绝对和相对，单位对不上，合成出来的结果再漂亮也是错的。

那列不是装饰

回到被打回的那份 FAT 报告。

补上灵敏度系数之后，不只是合规了。买卖双方发现，温度测量对最终不确定度的贡献，远超双方的预期。后续技术讨论的焦点，从"流量计精度够不够"，转移到了"温度测量方案是不是过于简化"。

这才是灵敏度系数的工程价值：它告诉你，当系统不确定度超标，改哪里有效，改哪里是无用功。

下次做不确定度预算，遇到那列 c：测量模型能写出 Excel 公式就用数值扰动法，δx 取标准不确定度，三列搞定。记得在实际测量值处算，记得对齐绝对和相对。

空着不填，或者全填 1，等于告诉对方——我不清楚这套系统的弱点在哪里。