3. 压力瞬变分析基础
在长输管道运行中，只要流体速度发生非定常变化，压力就会随之波动，这种压力、流量随时间剧烈变化的现象统称为压力瞬变，工程上也常叫水击或水力过渡过程。它与稳态工况下"压力沿程缓慢下降"完全不同：瞬变是压力波在管道内往返传播、叠加、衰减的动态过程，时间尺度常在秒到分钟级，空间上表现为沿管线的压力波列。理解瞬变，是判断"一次压力异常究竟是泄漏、调阀还是泵机组动作"的前提。
一、基本概念与产生机理
瞬变的本质是动量守恒被打破。稳态时，管段两端压差与摩阻平衡；一旦流速在短时间内改变，流体惯性使压力必须重新分布，形成压力波。以阀门快速关小为例：下游流速骤降，紧邻阀门的流体被"憋住"，动能转化为压能，产生一个向上游传播的升压波；同时阀门下游侧则出现降压波。波到达上游端（如泵出口、储罐）后发生反射，再向下游传播，如此往复，直到摩阻把能量耗散掉。
工程上常用儒可夫斯基（Joukowsky）公式估算直接水击的压强增量：
Δp = ρ · a · Δv
ΔH = a · Δv / g （以水头表示，ΔH 单位 m）
其中：ρ 为介质密度(kg/m³)，a 为压力波速(m/s)，
Δv 为流速变化量(m/s)，g 取 9.81 m/s²举例（示例参数，非实测）：某原油管道 ρ≈850 kg/m³，波速 a≈1050 m/s，若阀门动作使流速在极短时间内下降 0.5 m/s，则 Δp≈850×1050×0.5≈0.45 MPa，折合水头约 53 m。可见即使流速变化不大，瞬变压力也可能相当可观，这正是泄漏检测中必须把瞬变与泄漏区分开的原因。
是否属于"直接水击"，取决于阀门关闭时间 Tc 与相长 Tr=2L/a 的关系（L 为管段长度）：Tc ≤ Tr 为直接水击，压力增量接近儒可夫斯基值；Tc > Tr 为间接水击，压力增量被削弱。这一判据直接决定瞬变分析的建模精度要求。
二、传播速度与衰减特性
压力波速不是常数，它由流体可压缩性与管壁弹性共同决定：
a = sqrt( K/ρ / (1 + (K·D)/(E·e)·C) )
K：流体体积弹性模量；D：管内径；e：壁厚；
E：管材弹性模量；C：管道约束系数（薄壁锚固管约取 1−μ²）对成品油管道，a 通常在 900～1200 m/s；对天然气，因气体可压缩性大，波速显著偏低，常为 300～450 m/s。波速直接决定相长：一条 50 km 的原油管道，a=1050 m/s 时 Tr=2×50000/1050≈95 s，意味着一次阀门动作引发的压力波要在管内往返近 1.5 分钟才完成一个周期。
衰减来自两方面：一是摩阻耗散，压力波在传播中不断被管壁剪切消耗；二是波形弥散与反射叠加，在分输点、变径、弯头、阀室处产生部分反射与透射，使陡峭波头逐渐变缓。工程上表现为：距离扰动源越远，压力波峰值越低、上升沿越平缓。这一"幅值随距离衰减、波形随传播展宽"的特征，是判断扰动源位置的重要线索。
三、瞬变分析在泄漏检测中的角色
泄漏检测系统（如负压波法、压力梯度法、瞬态模型法）本质上都在"读"压力瞬变信号。泄漏发生时，泄漏点处产生一个向两端传播的负压波，其幅值与泄漏速率、波速、介质密度相关；通过两端传感器捕获负压波到达的时间差 Δt，结合波速即可估算泄漏位置：x = (L + a·Δt)/2（x 为距上游端距离，Δt 为下游到达时刻减上游到达时刻，需先统一时间基准与符号约定）。
但现场压力信号里同时混有调阀、启停泵、分输切换等正常操作引起的瞬变。因此瞬变分析的核心任务不是"看到压力变化就报警"，而是建立正常瞬变的特征基线：包括典型操作下的幅值范围、相长、传播方向与衰减规律，再据此把"符合基线"的操作瞬变剔除，把"不符合基线、且具有负压波特征"的事件识别为疑似泄漏。实际做法通常包括：
- 建立波速与相长基线：按管段材质、壁厚、介质物性计算 a，并用已知操作（如计划性调阀）反算校核，形成各管段 Tr 参考值。
- 提取瞬变特征量：压力下降幅值、上升/下降沿时间、到达时刻差、是否出现往返反射、衰减速率。
- 设定判别阈值：以正常操作瞬变的统计分布（如幅值 P95、到达时差范围）为界，超出且方向一致者进入疑似泄漏队列。
- 多源交叉验证：结合流量平衡、上下游压力梯度变化、阀门状态信号，排除单一传感器误报。
综上，压力瞬变既是管道运行中必须防范的水力风险，也是泄漏检测最灵敏的"信号源"。掌握产生机理、波速计算与衰减规律，才能把压力曲线上的每一次波动，翻译成可判断、可定位、可复核的工程结论。