基于压力表接口的高压天然气管道颗粒物采样性能研究：一种安全实用的替代采样方法

《Fuels》：Particulate Matter Sampling Performance of High-Pressure Natural Gas Pipelines Based on Pressure Gauge Ports: A Study of a Safe and Practical Alternative Sampling Approach

　　

针对高压天然气管道中在线颗粒物采样存在的安全限制和操作挑战，研究人员提出并系统研究了一种利用现有压力表接口的非侵入式安全采样替代方案。通过实验和计算流体动力学（CFD）数值分析，研究了速度比R（采样速度与管道气体速度之比）、管道流速V₀和接口方向对采样效率η（采样浓度与真实浓度之比）的影响。结果表明，通过压力表接口采样表现出与传统等速采样相反的趋势：效率随R显著增加。在R=2时，低于10μm的颗粒被有效收集，较大颗粒的捕获显著改善。流场分析表明，效率提升源于上游的再循环区、下游高压区的扩展以及高R条件下的“壁面撞击-反弹-再捕获”机制。基于实验数据建立了经验效率公式，预测误差低于15%。在输气站现场测试中，应用该修正将评估的过滤效率从9.4%提高到50.2%，有效恢复了真实运行条件。这项研究为天然气管道中安全、方便和准确的颗粒物监测提供了可靠的理论基础和实用技术方案。

**研究背景与问题**
高压天然气管道是能源系统的核心载体，管道内颗粒物（主要源于腐蚀产物、地层碎屑、焊接残留及旋转机械磨损）不仅导致仪器损坏、流量计精度下降、过滤分离设备堵塞，还会引发压缩机和燃气轮机内部部件的侵蚀磨损，造成巨大经济损失和安全风险。准确测量颗粒物浓度与粒径分布是评估管道完整性、优化净化工艺及实施预测性维护的关键前提。然而，现有标准等速采样方法（如ISO 9096、EPA Method 5）在高压、易燃易爆的天然气工业环境中面临严重安全挑战：安装标准等速采样喷嘴需对管道进行直接开孔，操作成本高且存在泄漏爆炸风险，难以用于常规测量和大规模调查，导致管道内颗粒物数据长期匮乏。为避免直接开孔风险，研究人员探索了基于声学、光学或静电感应的在线监测技术，但这些设备昂贵、需校准，且在高压复杂组分天然气环境中适应性有限。利用现有工艺接口（如压力表接口）进行采样具有自然安全优势，但该接口结构（采样方向垂直于管道气流，引压管内径6–15 mm）导致严重非等速采样状态，颗粒轨迹因惯性效应发生显著偏离，采样效率大幅降低且难以预测。经典壁面采样理论（如Hangal & Willeke模型）假设薄壁孔口，无法解释真实压力表引压管的厚壁腔体结构及“撞击-反弹-再捕获”机制。因此，针对压力表接口这一特定工程结构，定量研究非等速采样效率规律、多因素耦合机制及可靠修正方法，成为推动该便捷采样方法工程应用的关键知识空白。为此，研究人员开展了本研究，旨在提供一种安全、实用、非侵入式的颗粒物采样解决方案，论文发表在《Fuels》。

**主要技术方法**
研究人员采用实验与计算流体动力学（CFD）数值模拟相结合的方法。实验装置包括固体粉尘气溶胶发生部分（使用ISO 12103 A2标准试验粉尘，密度2650 kg/m3）、气溶胶采样部分（含水平管和垂直管两个测试段，压力表接口与标准等速采样喷嘴间距500 mm，引压管和等速喷嘴内径均为8 mm，主管内径68 mm）、在线颗粒测量部分（基于米氏散射原理的Palas Welas 3000光学颗粒计数器）及条件控制调节部分。CFD采用Realizable k-ε湍流模型闭合方程，基于欧拉-拉格朗日方法单向耦合模拟颗粒运动，仅考虑流体曳力和重力，忽略虚拟质量力、压力梯度力、Basset力和Saffman升力（量级分析表明其影响可忽略），采用Haider和Levenspiel非球形颗粒曳力系数关联式，颗粒-壁面相互作用假设完全非弹性捕获（基线模型）。网格为全六面体结构化网格（中网格约72万单元），通过网格独立性验证。现场验证在输气站进行，利用垂直向上的压力表接口，管道直径900 mm、压力2.5 MPa、温度23–25°C、气体流速2 m/s，采样速度4 m/s（R=2），使用高压光学颗粒计数器（OPC）在线检测，颗粒主要成分为磁铁矿、针铁矿和钠长石（平均真密度3563 kg/m3）。

**研究结果**
**3.1 速度比的影响**
通过实验和CFD模拟，在固定管道流速V₀=2 m/s条件下，研究了R=1、2、3时的采样效率η随粒径d_p的变化。结果表明，采样效率η受速度比R显著控制，且与粒径d_p密切相关。等速采样（R=1）时，η随粒径增大单调递减；R=2时，采样速度增加增强了惯性收集，显著提高了各粒径颗粒的η，大颗粒提升更明显；R=3时进一步增大。对于≥3μm颗粒，增大R有助于提升η；对于≤3μm颗粒，η在不同R下均≥90%。流场分析揭示，高R下上游形成再循环区、下游高压区向外扩展，并出现“壁面撞击-反弹-再捕获”机制，使大颗粒有效捕获概率增加。颗粒入口截面分布分析表明，增大R使有效采样范围向管道中心扩展，提高了采样代表性。

**3.2 管道流速的影响**
通过CFD数值模拟，在固定R=2条件下，研究了管道平均流速V₀分别为2、3、5、7、10 m/s时的采样效率。结果表明，对给定粒径，采样效率随管道流速增大而降低。例如，8μm颗粒在V=2 m/s时η=86.5%，当V增大至10 m/s时降至49.0%。流场分析显示，不同流速下采样口附近流场具有自相似性（归一化速度与压力分布高度一致），但流速增大使颗粒Stokes数Stk增大，惯性效应增强，颗粒轨迹偏离流线更严重，导致捕获概率降低。入口截面分布分析表明，流速增大使有效采样范围向管壁移动，对大颗粒的空间约束效应更明显，需进一步增大R以补偿。

**3.3 压力表接口方向的影响**
通过实验，在V₀=2 m/s、R=2条件下，比较了水平与垂直向上两种接口方向的采样效率。结果显示，由于重力影响，垂直向上方向的η低于水平方向，且随粒径增大差异更显著（5–8μm颗粒垂直方向低5–15%）。大颗粒惯性大，向上运动需克服重力做功，更难进入采样管。现场验证中因条件限制使用垂直向上接口，经验公式（基于水平方向实验）可能引入额外不确定性（1–5μm颗粒估计±10%）。

**3.4 压力表接口采样效率经验公式**
基于实验数据（78个数据点，粒径1–8μm，R=1、2、3，V₀=2 m/s），采用非线性最小二乘回归对Hangal & Willeke 90°采样效率公式进行修正，得到压力表接口几何下的经验公式：η = 1 – (1 – 1/(1 + 2.4·Stk^0.3·R))·(1 – 1/R)。该公式适用于0.002 < Stk < 0.1和1 < R < 3，决定系数R2=0.94，预测误差一般在±15%以内。对Stk>0.1（对应粒径约>4.5μm）外推可能低估效率，需进一步验证。

**讨论与结论**
讨论部分指出，本研究在机理理解和工程应用上取得明确进展，但存在局限性：经验公式主要基于实验室标准粉尘，对真实管道中复杂成分颗粒的适用性及垂直方向接口的修正方法需进一步发展。未来工作应扩展公式至高Stokes数范围，并开展多站点长期现场验证以增强模型普适性。研究结论翻译如下：本研究结合实验测量、CFD模拟和现场验证，系统研究了高压天然气管道压力表接口的颗粒物采样。压力表接口提供了一种安全、非侵入式、实用的采样替代方案，无需管道改造，特别适用于传统等速采样不可行的在线颗粒检测。与经典等速行为相反，压力表接口的采样效率η随速度比R增大而增加；在R=2时，低于10μm颗粒被有效收集，5–8μm颗粒的捕获随R增大显著改善，这归因于再循环区的形成和“壁面撞击-反弹-再捕获”机制。基于实验数据（1–8μm颗粒，R=1–3，V₀=2 m/s）建立了经验效率公式η = 1 – (1 – 1/(1 + 2.4·Stk^0.3·R))·(1 – 1/R)，适用于0.002 < Stk < 0.1和1 < R < 3，预测误差一般在±15%以内（R2=0.94）；在Stk=0.1外推时不推荐使用。在输气站现场应用（2.5 MPa，垂直方向接口）表明，原始未修正数据因非等速采样偏差严重低估颗粒浓度，应用经验修正后，过滤器入口和出口浓度分别从53.0和48.0 mg/m3增至551.1和274.4 mg/m3，过滤效率从9.4%提高到50.2%，反映了实际过滤器退化状态。对于垂直方向接口，经验公式（基于水平方向实验）可能引入额外不确定性（1–5μm颗粒估计±10%），未来工作应开发方向特定修正因子。经验公式主要基于实验室标准粉尘（密度2650 kg/m3，粒径1–8μm），对复杂成分和更高Stokes数颗粒的适用性需进一步研究；此外，CFD模拟假设完全非弹性颗粒-壁面捕获，仅用于定性趋势分析。未来研究应纳入恢复系数和临界捕获速度以更好模拟反弹-再捕获过程，并在多站点进行长期现场验证以增强模型普适性。总体而言，本研究为利用现有基础设施进行高压天然气管道颗粒物监测提供了具有理论依据、操作安全且实际可靠的技术方案。