管道材料与液体物性对水锤波传播及瞬态压力响应的影响：基于OFAT敏感性分析与决策树模型的定量评估

《Materials》：Influence of Pipe Material and Liquid Properties on Water-Hammer Wave Propagation and Transient Pressure Response

　　

摘要专业翻译 pressurized pipelines 中的水锤瞬变过程由液体可压缩性与管壁弹性响应的耦合作用共同决定。本研究基于一维特征线法（Method of Characteristics, MOC）模型，并引入基于卷积形式的非定常摩擦公式，定量评估了管道材料、系统几何参数及液体物性变化对瞬态压力响应的影响。该模型应用于层流和湍流条件下的实验铜管系统。研究采用单因素逐一变化（one-factor-at-a-time, OFAT）敏感性分析方法，对杨氏模量（Young's modulus）、泊松比（Poisson's ratio）、管壁厚度、内半径、液体密度、体积弹性模量（bulk modulus）及运动黏度等参数分别独立施加±5%和±10%的扰动。结果表明，液体密度、液体体积弹性模量和管道杨氏模量对压力波速（pressure-wave celerity）及瞬态响应的时间演化具有显著影响；管道内半径对压力响应具有显著作用，而黏度和体积弹性模量的影响在湍流条件下更为突出。材料、管道几何及液体物性的综合变化使压力波速范围约为1159～1479 m/s，表明输入参数的不确定性对预测瞬态行为具有实质性影响。进一步采用多值逻辑（multiple-valued logic, MVL）决策树和分类回归树（classification and regression tree, CART）模型对OFAT结果进行分析，以识别主导参数-响应关系及阈值。研究结果强调，准确表征管道材料与液体物性对于可靠的水锤预测及 pressurized pipeline 系统评估具有重要意义。

论文解读

一、研究背景与问题

水锤（water hammer）现象是 pressurized pipeline 系统中因流速突变而产生的压力瞬变过程，其传播特性由液体可压缩性与管壁弹性变形的耦合作用共同决定。准确预测水锤压力波速及瞬态压力响应，对于高圧液压系统、区域供热装置、核反应堆冷却系统及长距离输水管线的安全设计至关重要。然而，在实际工程中，管道材料的力学参数（如杨氏模量、泊松比）、几何尺寸（内半径、壁厚）以及液体物性（密度、体积弹性模量、运动黏度）往往难以精确获取，且这些参数随温度、压力及制造公差而变化，导致预测结果存在显著不确定性。现有研究虽已认识到管道材料对水锤行为的影响，但对液体物性与管道参数的系统性敏感性比较仍较为缺乏。此外，传统敏感性分析方法多局限于参数逐一比较，缺乏将定量敏感性评估与可解释性决策规则相结合的系统框架。因此，本研究旨在通过整合OFAT敏感性分析、MVL决策树和CART回归模型，系统评估管道材料、几何参数及液体物性对水锤瞬态响应的影响，并建立可解释的参数-响应关系模型。

二、研究内容与主要结论

研究人员基于一维特征线法（MOC）建立水锤瞬变模型，并引入基于卷积形式的非定常摩擦公式，以实验铜管系统（管长L = 98.11 m）为对象，分别在层流（初始流速v₀ = 0.066 m/s）和湍流（v₀ = 0.94 m/s）条件下开展数值模拟。采用OFAT方法对7个关键参数（液体密度ρ、管道内半径R、壁厚e、泊松比ν_P、运动黏度ν、液体体积弹性模量K、杨氏模量E₀）分别施加±5%和±10%的扰动，以第15个压力峰值的时间T₁₅和幅值p₁₅作为响应指标。研究得出以下主要结论：

第一，对于时间响应T₁₅，参数影响排序为ρ > K > E₀ > e > R > ν_P > ν，且该排序在层流和湍流条件下基本一致，表明压力波传播的时间结构主要由液体压缩性和管壁弹性决定，黏性效应的影响相对有限。

第二，对于压力响应p₁₅，管道内半径R在两种流态下均为最主导参数；但在湍流条件下，液体体积弹性模量K和运动黏度ν的影响显著增强，其归一化影响系数分别从层流条件下的0.233%和0.311%增至湍流条件下的2.207%和1.806%，表明流体物性对压力响应的贡献具有流态依赖性。

第三，当管道参数与液体物性综合变化时，压力波速范围为1159～1479 m/s，差异达320 m/s，说明输入参数不确定性可显著改变水锤事件的预测时间结构。

三、关键技术方法

本研究采用的主要技术方法包括：（1）一维特征线法（MOC）数值模型，用于求解水锤瞬变控制方程；（2）基于卷积积分的非定常摩擦公式，用于表征瞬变流动中的历史相关黏性耗散效应；（3）单因素逐一变化（OFAT）敏感性分析方法，独立扰动各参数并量化其对响应的影响；（4）多值逻辑（MVL）决策树，将连续响应转化为定性分类并提取可解释的IF-THEN规则；（5）分类回归树（CART）模型，通过递归方差最小化准则识别主导分裂变量和响应阈值。模型基准参数取自Adamkowski和Lewandowski在格但斯克IMP PAN开展的实验研究，管道离散为N = 32个计算段。

四、研究结果

4.1 OFAT敏感性分析

OFAT分析结果表明，时间响应T₁₅对液体密度ρ最为敏感，其次为体积弹性模量K和杨氏模量E，壁厚e和内半径R的影响较低，泊松比ν_P和运动黏度ν的影响可忽略。压力响应p₁₅的敏感性则呈现不同的模式：层流条件下R为主导参数，其次为ν和K；湍流条件下R仍为主导，但K和ν的影响显著增大。值得注意的是，层流条件下ρ和ν_P的归一化影响系数为零，这归因于压力值以MPa为单位保留三位小数时掩盖了微小变化，而非物理上无影响。

4.2 规则基决策结构

基于OFAT排序结果，研究分别针对T₁₅和p₁₅构建了四套独立的规则基决策树（层流/湍流×时间/压力）。对于T₁₅，ρ和K提供了Short与Long类别的清晰分离：ρ减小导致Short响应，ρ增大产生Long响应；K的变化方向相反。杨氏模量E的扰动始终落在Medium类别内。对于p₁₅，R提供最清晰的类别分离：R减小对应Low响应，R增大对应High响应；ν的影响方向相反，K呈渐进正向影响。层流条件下p₁₅值集中在约1.283×10⁶～1.290×10⁶ Pa的窄范围内，而湍流条件则扩展至约1.468×10⁶～1.525×10⁶ Pa。

4.3 CART回归模型

CART回归模型独立构建于四个响应案例。对于T₁₅，两种流态下ρ均被识别为第一分裂变量，K提供进一步细化，与OFAT排序一致。对于p₁₅，R为第一分裂变量，K在参数空间中间区域提供细化；但层流条件下CART模型选择K作为第二重要变量，与OFAT排序中ν为第二重要变量存在细微差异，这反映了两类模型数学机制的不同。模型性能方面，决定系数R²介于0.657～0.717之间，表明模型可解释大部分响应变异但仍有约28%～34%未解释。留一交叉验证（LOOCV）的R²值仅为0.028～0.115，表明模型的预测泛化能力有限，应将其定位为探索性和解释性工具而非独立验证的预测模型。

五、讨论与结论

研究讨论指出，时间响应的参数排序在层流与湍流间基本一致，说明压力波传播时间结构主要由液体压缩性和管壁弹性主导；而压力响应的敏感性则显著依赖流态，湍流条件下流体可压缩性和黏性效应的贡献大幅增强。研究强调，杨氏模量和泊松比不应仅被视为次要输入参数，它们通过管壁柔度贡献直接耦合管道材料的力学行为与输送液体的可压缩性。

研究结论可概括为：时间响应T₁₅主要由液体密度、液体体积弹性模量和管道杨氏模量控制；压力响应p₁₅以管道内半径为最主导参数，且湍流条件下体积弹性模量和运动黏度的影响显著增强；参数综合变化导致压力波速在1159～1479 m/s范围内变动，差异达320 m/s，表明材料、几何及液体物性不确定性可显著改变水锤事件的预测时间结构。本研究的核心贡献在于整合OFAT定量敏感性评估、定性规则基决策结构和CART回归树建模三种互补方法，为层流和湍流条件下水锤参数影响比较提供了系统框架。但需指出，受限于OFAT单参数扰动设计及有限的数据集规模（每流态29个配置），该框架应视为在研究所涉参数范围内的探索性和解释性方法，尚不足以作为独立验证的工程预测工具。未来工作应考虑代表性工程管道材料、温度依赖的材料与液体物性、更宽参数范围以及因子设计或响应面实验设计，以更全面地解析多参数交互效应。