Nonsmooth Hydraulics, Smooth Control: System Theory Framework for Analyzing Water Networks

技术摘要：非光滑水力学与平滑控制：用于分析水网系统的系统理论框架

问题陈述
配水网络（WDNs）是受控液压系统，其中泵速和阀门模式引入了控制输入、代数约束以及切换事件。尽管经典的水力分析算法和运行控制算法已相当成熟，但关于配水网络为何会表现出特定行为的严格控制理论理解仍然有限。现有工具（例如 EPANET）依赖于准稳态模拟或扩展周期模拟，这些方法在离散时间步长上求解代数系统，往往忽略管道惯性并将设备切换视为瞬时跳变。这导致在将配水网络作为连续时间系统进行分析时，存在局部正则性、可微性、稳定性和可控性分析方面的空白，特别是在泵和阀门切换模式引入非光滑性从而使得标准控制分析复杂化的情况下。

方法论
作者基于连续时间非线性微分代数方程（DAE）模型开发了一个综合框架。该方法论分为四个主要阶段：

1. DAE 建模：本文推导了一个完整的非线性 DAE 模型，该模型以未约化的状态空间形式保留了管道、泵和阀门的流量，节点和储罐的水头，储罐的蓄水量以及节点需水量。该模型介于稳态代数求解器和完整水锤瞬态模型之间，利用“刚性水柱”近似，既包含了管道惯性和储罐蓄水量，又忽略了声波传播。
2. 正则性与平滑化：作者在固定控制下分析了 DAE 的局部适定性，证明了若代数雅可比矩阵非奇异，则系统是局部指数为 1 的，并存在唯一轨迹。他们表明，泵和阀门的硬切换会破坏状态轨迹的可微性。为解决这一问题，作者引入了一种五次平滑正则化技术。该技术将 abrupt 的控制更新替换为窄时间窗口内的平滑过渡，确保约化后的向量场保持可微，从而适用于连续时间分析。
3. 图论线性化：平滑后的 DAE 在运行点附近进行线性化。作者表明，该线性化诱导了一种加权图拉普拉斯结构。该拉普拉斯的权重源自运行点处的增量水力导纳（水头损失曲线的斜率），从而明确了网络互连和能量耗散特性。
4. 鲁棒性与裕度分析：该框架量化了参数不确定性（例如管道粗糙度、泵系数）如何影响线性化模型。利用隐函数定理的论证，作者建立了平衡点偏移和矩阵摄动的二次误差界。此外，他们还推导了显式界限，说明这些摄动如何影响局部稳定性裕度（指数衰减率）和有限时域可控性裕度。

主要贡献
本文对文献做出了四项主要贡献：

1. 统一非线性 DAE 模型：一个连续时间 DAE 模型，集成了管道惯性、储罐蓄水量和设备切换，而无需将系统简化为纯代数形式。作者证明了在固定控制下系统的局部指数为 1 的正则性以及唯一的水力演化。
2. 平滑控制正则化：一种处理设备切换非光滑性的方法，通过用可微的五次插值替换硬性的开/关过渡。这使得可以使用标准的连续时间控制分析工具（线性化、稳定性、可控性），而这些工具原本会在切换瞬间失效。
3. 图论稳定性分析：证明了线性化后的平滑 DAE 由加权拉普拉斯矩阵支配。只要每个活跃子网都锚定在水库上，这种结构就允许基于李雅普诺夫函数的局部渐近稳定性证明。该分析直接将网络拓扑与能量耗散率联系起来。
4. 参数灵敏度与裕度界限：建立了参数摄动下线性化的显式误差界，并推导了稳定性和可控性裕度得以保持的条件。这包括一种实用程序，用于根据参数对系统裕度的关键性对网络参数进行排序。

结果
该理论框架通过在六个基准网络（包括 EPANET 示例和 KIOS 基准）上的广泛案例研究进行了验证，使用了与 EPyT 接口的 Python 实现。

• 模型保真度：平滑 DAE 轨迹与 EPANET 模拟结果高度吻合。对于固定控制场景，水头和流量偏差极小（水头偏差在厘米级）。对于 24 小时平滑控制模拟，误差保持在可接受范围内，偏差主要集中在切换事件周围，因为平滑窗口引入了有限斜坡而非瞬时跳变。
• 稳定性验证：数值测试证实，线性化 DAE 矩阵的特征值位于开左半平面，这与理论上的李雅普诺夫稳定性证明一致。加权拉普拉斯结构成功捕捉了能量耗散特性，导纳权重的分布与耗散率的分布相关联。
• 参数灵敏度：误差界分析（过程 1 和 2）成功识别出，对于适度的参数摄动（例如，粗糙度变化高达约 10%），线性化仍然有效。至关重要的是，灵敏度分析揭示了对稳定性和可控性最关键的管道，并不一定是流量最高的那些，而是增量水头损失斜率对参数变化最敏感的那些。
• 需水量影响：在剧烈的需水量变化下，系统保持了局部稳定性和正的可控性裕度，尽管随着需水量增加，泵的权威性（增益）略有下降，这与高流量下更高的耗散效应一致。

意义与主张
本文声称提供了一个系统理论基石，用于分析配水网络水力学，从而弥合了运行模拟与控制理论之间的鸿沟。

• 它超越了“黑盒”模拟，提供了可解释的、图结构的线性化，揭示了拓扑和水力参数如何共同决定稳定性和可控性。
• 它解决了配水网络中非光滑性的具体挑战，提供了一条数学上严谨的路径（通过平滑化），将连续时间控制分析应用于具有切换设备的系统。
• 作者强调，他们的贡献并非孤立地发明新的控制理论要素，而是将这些基础（DAE 理论、图拉普拉斯、摄动界限）定制并组装成一个专门针对配水动力学的统一框架。
• 这项工作通过量化不确定性和运行点变化如何影响系统裕度，使网络分析、缓解策略和系统设计成为可能，提供了一种在不依赖纯启发式方法的情况下识别关键管道以进行校准和设计的工具。

局限性与未来工作
作者承认，该理论是局部的，仅针对特定的运行点和活跃设备模式。该模型目前未包含水锤效应、压力依赖性泄漏或压力依赖性需水量。此外，虽然该框架诊断了稳定性和可控性裕度，但它并未规定具体的系统级修改措施（例如管道升级或控制节点布置）以改善这些裕度。未来的工作建议将分析扩展到完整的非线性 DAE，纳入更复杂的水力现象，并基于这些理论见解开发网络修改方法。