4. 孔板差压变送器阻尼时间对脉动流计量的校核

脉动流下阻尼时间常数的响应衰减机理

孔板流量计的输出差压与流量呈平方关系，当管道中存在压缩机往复、调压阀振荡或用户负荷快速切换时，差压信号会叠加周期性脉动分量。差压变送器内部的阻尼电路本质上是一个一阶低通滤波器，其时间常数 τ 决定了变送器对脉动分量的衰减能力，同时也决定了它对真实阶跃变化的跟随速度。

一阶阻尼环节的幅频特性为：

A(f) = 1 / √(1 + (2πfτ)²)

其中 f 为脉动频率，τ 为阻尼时间常数。当 2πfτ ≫ 1 时，脉动幅值被显著压缩。问题在于：孔板流量方程中差压经开方后才得到流量，而开方是凹函数，脉动差压的均值经开方后并不等于脉动流量的均值——阻尼把差压脉动"压平"后，流量计算机读到的差压更接近脉动均值，但真实流量应按瞬时差压开方后再平均。两者之差就是阻尼引起的计量偏低量。

校核计算：用示例参数量化偏低量

设孔板设计差压 ΔP₀ = 25 kPa，脉动为近似正弦分量，幅值 Δp = 4 kPa，频率 f = 1.2 Hz，流量系数按 q = k·√ΔP 处理（k 取 1，只比较相对值）。对每个阻尼时间常数，先算衰减系数 A，再算阻尼后差压的均方根开方值，与"先开方后平均"的真实值比较。

import math

dP0 = 25.0      # 设计差压 kPa
dp  = 4.0       # 脉动幅值 kPa
f   = 1.2       # 脉动频率 Hz
N   = 2000      # 一个周期采样点数

def true_mean_flow():
    # 先开方后平均：真实平均流量
    s = 0.0
    for i in range(N):
        dP = dP0 + dp * math.sin(2*math.pi*i/N)
        s += math.sqrt(dP)
    return s / N

def damped_flow(tau):
    # 一阶低通离散仿真：先阻尼后开方
    dt = 1.0 / (f * N)          # 采样步长
    alpha = dt / (tau + dt)     # 一阶滤波系数
    y = dP0                     # 初始值取均值
    s = 0.0
    for i in range(N):
        x = dP0 + dp * math.sin(2*math.pi*i/N)
        y = y + alpha * (x - y) # 阻尼响应
        s += math.sqrt(y)
    return s / N

q_true = true_mean_flow()
for tau in (0.1, 0.5, 2.0):
    q_d = damped_flow(tau)
    eps = (q_true - q_d) / q_true * 100
    print(f"tau={tau:>4}s  q_true={q_true:.4f}  q_d={q_d:.4f}  偏低率={eps:.2f}%")

示例运行结果（由上述给定输入直接推出，非实测）：

现场校核操作与结果判读

1. 记录当前阻尼设置：从变送器手操器或组态软件读取阻尼时间常数 τ，同时记录流量计算机中的差压量程与开方配置。
2. 采集脉动差压：用高采样率记录仪（采样间隔不大于脉动周期的 1/20）在变送器输出端连续记录 2~3 分钟差压，统计均值 ΔP₀、峰峰值与主频 f。
3. 计算衰减系数：按 A(f)=1/√(1+(2πfτ)²) 估算当前 τ 对主频脉动的衰减，若 A 明显小于 1 且脉动幅值占比大，则存在偏低风险。
4. 对比验证：将 τ 临时调至最小（如 0.05 s），稳定后读取流量计算机输出，与调整前比较。若小 τ 下流量读数上升且更接近按瞬时差压开方平均的结果，说明原阻尼设置造成偏低。
5. 确定最终设置：在抑制噪声与保证响应之间折中，通常取 τ 使主频脉动衰减后残余幅值不影响开方平均，同时阶跃响应时间满足计量响应要求。