供水管网动态模型分类及应用分析

2024-01-13 02:44

本文主要是介绍供水管网动态模型分类及应用分析,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

       当供水管网中发生启停泵、快速关阀等事件时, 延时时段模拟 (即准稳态模型) 不能准确预测系统的瞬时动态变化, 而需要采用更为准确复杂的瞬变流动态模型。为明确多种动态模型之间的差异, 探讨和分析了供水管网动态模型的分类、模型理论以及在管网运行管理中的应用。结果表明, 准稳态模型适用于水力状态缓慢变化的情况;弹性水锤模型可准确模拟系统的瞬时动态特性, 包括流体惯性和系统弹性 (即流体可压缩性和管壁弹性) ;刚性水柱模型只考虑了流体的惯性作用, 可作为前两者的过渡。在实际应用中, 3种动态模型应结合水力状态变化情况和模拟准确性需求合理选用, 为供水管网系统提供更精细的运行管理方式。

   目前, 计算机模拟分析已经成为供水管网系统设计、运行和管理的基本工具, 用于监测或预测供水管网在不同运行工况下的压力、流量、水池水位等水力状态。在管网建模中, 一般认为供水管网的运行状态是缓慢变化的, 可以通过连续模拟整时段的平均状态来表示系统的动态变化, 如间隔时间5 min或15min的延时时段模拟 (Extended Period Simulation, EPS) 。但是, 实际的管网运行状态并非是完全缓慢变化的过程, 存在着很多种水力瞬变过程, 比如水泵的启闭、阀门的快速调节、节点用水量的突然波动等。当水力瞬变过程 (即非恒定流) 发生时, 有可能引起系统压力和流量的不正常波动, 给系统的安全运行带来威胁。而常规的延时时段模拟方法不能准确预测到瞬变波动, 需要进一步采用更为准确的瞬变流计算模型。目前国内大多数管理人员容易忽视这点。因此, 本文主要探讨供水管网动态水力模型的分类、模型理论及在不同场景下的应用, 为供水管网系统的安全运行管理提供更加全面的技术指导和理论补充。

1 供水管网动态模型分类

   有压管网系统中的水体流动实质上是一种非恒定流动, 其水力状态变化可根据强度划分为缓慢变化和快速变化过程。一般来说, 流体状态的缓慢变化可认为系统的水力状态是轻微非恒定, 流体惯性和系统弹性 (包括流体可压缩性和管壁弹性) 的影响可忽略不计, 而流体状态的快速变化需要考虑到流体惯性和系统弹性的影响。因此, 对应于不同的状态变化强度假设, 供水管网动态模型可分为三大类:准稳态模型 (即延时时段模拟) ———假设流体动态作用可以忽略不计, 适用于模拟缓慢变化过程中的水力状态;弹性水锤模型———考虑流体的惯性和系统弹性的影响, 适用于模拟快速变化过程中的水力瞬变过程;刚性水柱模型———介于上述两种模型之间, 只考虑流体的惯性作用, 适用于模拟系统的弹性作用影响较小时的水力瞬变状况[1]。3种动态模型的具体分类方法见图1。

图片

 图1 管网动态模型分类示意   

   如图1所示, 随着模型准确性要求的提高, 模型的复杂度随之提高。与此同时, 模型求解的时间尺度逐步由小时级或分钟级的平均状态“精细”到以秒计的水力瞬态过程。准稳态模型可以快速得到系统在长时段的连续状态变化;而弹性水锤模型可以准确地预测系统在短时段的瞬时状态变化情况;刚性水柱模型介于两者之间, 可以作为两者的过渡模型。另外, 弹性水锤理论尚处于逐步发展完善的状态, 非恒定摩阻、管道粘弹性等项的加入会使模型的准确度和复杂度得到进一步的提升。由此可见, 供水管网动态模拟需要在模型复杂度和准确性之间作出合理的均衡以选择最合适的模型。而如何在两者之间作出最合理的均衡是学者和专业人员们所必须面临的问题, 尤其是针对复杂的管网系统尚缺乏明确的理论和技术指导。一般来说, 我们根据事件所引起的水力状态变化强度特性来选择合理的动态模型, 如图1中的水力状态变化时间尺度所示。

2 动态模型理论

2.1 基本数学模型

   描述有压管流的一维非恒定流动过程的连续方程和运动方程表达如式 (1) 所示[2]:

图片

   式中H———瞬时水头, m;

   Q———瞬时流量, m3/s;

   X———距离, m;

   T———时间, s;

   f———管道非恒定流水力摩阻系数, 假定只考虑稳态摩阻, 则为达西-魏斯巴赫系数;

   D———管道的直径, m;

   A———管道截面积, m2;

   a———水锤传播速度, m/s;

   g———重力加速度, m/s2。

   式 (1) 为准线性双曲型偏微分方程组, 一般采用特征线法 (Method of Characteristic, MOC) 求解一系列边界条件下的压力和流量波动情况, H=H (T, X) 和Q=Q (T, X) , 即弹性水锤模型。

   如果流体为非恒定不可压缩流体, 假设波速a→∞ (接近于一个真实的刚性管道系统) , 则式 (1) 中的连续方程简化为Q/X=0。由此, 瞬态流量函数简化为式 (2) :

图片

   则式 (1) 中的运动方程可重新整理为:

图片

   式中L———管道长度, m;

   Hf———管道阻力损失, m;

   u、d———分别表示管道的上游节点和下游节点。

   显然, 式 (3) 表示管道两端的瞬态压差等于管道的阻力损失和流体加速水头 (即惯性作用) 之和, 即刚性水柱模型。

   进一步, 如果忽略流体的惯性作用 (dQ/dT→0) , 式 (3) 可简化为式 (4) :

图片

   即是常见的达西-魏斯巴赫水头损失公式。利用式 (4) 表示管道上下游之间水头的缓慢变化, 即是准稳态模型。

2.2 模型求解方法

   供水管网动态模型求解是对相应的控制方程组进行求解:弹性水锤模型求解流体的连续方程和运动方程, 即式 (1) ;刚性水柱模型和准稳态模型求解管段能量方程和节点流量连续性方程, 即式 (3) 或式 (4) 与节点连续性方程。弹性水锤模型的求解一般是将双曲型偏微分方程组转化常微分方程组进行求解, 如常用的特征线法, 具体可参考经典文献[3, 4]。准稳态模型是稳态单时段模型的动态连续, 可对稳态模型的求解方法进行适当改进, 比如广泛应用的全局梯度算法 (Global Gradient Algorithm, GGA) 。考虑到管网状态的动态变化, Todini[5]和Giustolisi等[6]将可变水位水池引入到GGA中, 提出了更为一般化的全局梯度算法 (GeneralGGA, G-GGA) 用于准稳态动态模拟。进一步, Nault等[7]将刚性水柱模型与GGA整合, 形成了可模拟刚性水柱的全局梯度算法 (Rigid Water Column-GGA, RWC-GGA) 。本文选用G-GGA用于准稳态模型, RWC-GGA用于刚性水柱模型, MOC用于弹性水锤模型。

3 应用实例分析

3.1 模型描述

   为表达3种不同类型动态模型在供水管网中的应用, 选择一个实例管网进行分析。实例管网选取自经典水力计算软件EPANET 2.0中的示例管网Net1, 如图2所示。为方便进行计算分析, 对管网属性进行适当修改, 未提及的部分与原模型相同。各管段属性如表1所示。水池初始水位和水池直径分别为295m和15m。水泵出口阀门直径为450mm, 初始为全开状态, 阻力系数为1.0。该水力模型模拟了管网在24h的水力状态变化, 其中高位水池作为水力调蓄设施决定着水泵的启停动作。当水池水位接近最高水位时 (大于300m) , 水泵停止工作;当水池水位接近最低水位时 (小于290 m) , 水泵重新启动。由此, 水池水位和水泵出口压力在1天中的变化过程可用延时时段模拟得到, 如图3所示。

图片

 图2 管网拓扑结构示意   

 表1 管网管段信息    

图片

   从图3可知, 水池水位的变化引起水泵在10:33关闭和在19:19开启的动作。在准稳态动态模型 (即EPS) 中, 水泵的启停动作引起水泵出口处的水头发生“跳跃”变化, 水力波动过渡过程被简化忽略。然而在实际运行中, 如果水泵启停动作不当, 可能会引起管网产生明显的水力瞬变波动, 甚至造成严重事故。这时仍然采用准稳态动态模型进行预测显然不合适。下面将通过不同动态模型的模拟计算, 分析动态水力模型之间的差异性以及传统的准稳态动态模型的局限性。

图片

图3 水池水位和水泵出口处水头在24h内的变化曲线  

3.2 动态模型分析

   假设水泵在启停动作中遵循闭阀操作的原则, 即在水泵出口阀门完全闭合的情况下执行停泵或开泵操作。因此, 阀门动作过程关乎管网系统的动态变化情况。以10:33的停泵关阀动作为例, 对比分析不同的关阀动作过程下的管网水力状态的动态变化过程。采取3种不同的两阶段关阀方案, 分别是: (1) 30s关阀总历时, 其中6s关闭80%, 24s关闭20%; (2) 90s关阀总历时, 其中18s关闭80%, 72s关闭20%; (3) 300s关阀总历时, 其中60s关闭80%, 240s关闭20%。

   图4表示了水泵出口处在3种不同关阀方案下动态模型的模拟计算结果 (弹性水锤模型中, 管道波速为1 000m/s, 时间步长为0.1s;准稳态模型和刚性水柱模型的时间步长均为1s) 。根据弹性水锤计算结果, 随着关阀时间的延长, 停泵关阀事件所产生的动态压力波动逐渐减小, 由30s关阀时的55 m左右降低至300s关阀时的5 m左右。同时, 刚性水柱模型的计算结果 (初始压力波相) 逐渐接近弹性水锤模型, 说明在长时段关阀过程中, 系统弹性可忽略不计, 即刚性水柱模型可替代弹性水柱模型以减少计算量。当停泵关阀动作足够缓慢 (如300s关阀) , 流体的真实状态虽仍存在动态波动, 但波动很小 (见图4c) , 实际应用中采用准稳态模型近似代替仍可保证足够精度。由此可见, 停泵关阀过程中的系统动态变化应根据水力状态变化强度特性选用相应的动态模型进行准确模拟。

   进一步, 将刚性水柱模型和弹性水锤模型与准稳态模型相结合, 在启停泵动作发生时采用, 以准确连续地表达管网系统在24h的动态变化情况。停泵关阀采用90s历时的两阶段关阀方案 (18s关闭80%, 72s关闭20%) , 启泵开阀采用与之相反的开阀方案 (即72s开启80%, 18s开启20%) 。最终, 水泵出口处的水头连续变化状态如图5所示。

图片

   图4 水泵出口处水头波动情况   

图片

   图5 动态模型修正的水泵出口处水头动态变化过程   

   由图5可知, 经瞬态动态模型修正后的总水头变化曲线更能准确地反应系统状态在全尺度时间范围内的动态变化情况, 包括系统状态的平稳变化过程和快速变化过程。由此, 3种动态模型的结合使用可以为管网的运行管理提供更加精细的管理方案, 保障系统的合理安全运行。

4 结论

   (1) 供水管网动态水力模型可分为三大类:准稳态模型、刚性水柱模型和弹性水锤模型。准稳态模型 (即延时时段模拟) 通过将单时段静态模型连接起来形成连续时段模型, 可近似模拟缓慢变化的水力过程;刚性水柱模型和弹性水锤模型分别部分和全部地考虑了流体的动态特性, 即前者考虑流体的惯性作用, 后者同时考虑流体惯性和系统弹性的作用, 因而模型准确度和复杂度都得到了提升。

   (2) 通过动态模型在实例管网中的应用发现:在快速水力状态变化事件中, 应采用弹性水锤模型准确模拟系统的瞬时状态变化;随着水力状态变化尺度增大, 刚性水柱模型可近似替代弹性水锤模型以减少计算量;在水力状态变化足够缓慢的情况下可以采用准稳态模型进行水力状态的模拟与预测。

   (3) 实例应用结果表明了3种动态模型在不同应用场景下的模拟效果, 也展示了一种将3种模型相互结合用于系统全尺度连续动态模拟的形式, 可以更加全面的预测系统的动态变化。据此提出, 管理人员应根据管网实际运行工况 (即水力状态变化强度) 和模拟准确度需求结合使用3种动态模型, 为供水管网系统提供更精细的运行管理方式。

文章来源icon-default.png?t=N7T8https://mp.weixin.qq.com/s?__biz=Mzg3NzkxNTI1MA==&mid=2247485556&idx=1&sn=32b08f7df59948bc77a40bb1a2288579&chksm=cf1af569f86d7c7f549613ff4be835e70a022ec722db7ea140454cddad35df135ac63da8c95a&token=2109378225&lang=zh_CN#rd

这篇关于供水管网动态模型分类及应用分析的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/599979

相关文章

大模型研发全揭秘:客服工单数据标注的完整攻略

在人工智能(AI)领域,数据标注是模型训练过程中至关重要的一步。无论你是新手还是有经验的从业者,掌握数据标注的技术细节和常见问题的解决方案都能为你的AI项目增添不少价值。在电信运营商的客服系统中,工单数据是客户问题和解决方案的重要记录。通过对这些工单数据进行有效标注,不仅能够帮助提升客服自动化系统的智能化水平,还能优化客户服务流程,提高客户满意度。本文将详细介绍如何在电信运营商客服工单的背景下进行

中文分词jieba库的使用与实景应用(一)

知识星球:https://articles.zsxq.com/id_fxvgc803qmr2.html 目录 一.定义: 精确模式(默认模式): 全模式: 搜索引擎模式: paddle 模式(基于深度学习的分词模式): 二 自定义词典 三.文本解析   调整词出现的频率 四. 关键词提取 A. 基于TF-IDF算法的关键词提取 B. 基于TextRank算法的关键词提取

基于人工智能的图像分类系统

目录 引言项目背景环境准备 硬件要求软件安装与配置系统设计 系统架构关键技术代码示例 数据预处理模型训练模型预测应用场景结论 1. 引言 图像分类是计算机视觉中的一个重要任务,目标是自动识别图像中的对象类别。通过卷积神经网络(CNN)等深度学习技术,我们可以构建高效的图像分类系统,广泛应用于自动驾驶、医疗影像诊断、监控分析等领域。本文将介绍如何构建一个基于人工智能的图像分类系统,包括环境

水位雨量在线监测系统概述及应用介绍

在当今社会,随着科技的飞速发展,各种智能监测系统已成为保障公共安全、促进资源管理和环境保护的重要工具。其中,水位雨量在线监测系统作为自然灾害预警、水资源管理及水利工程运行的关键技术,其重要性不言而喻。 一、水位雨量在线监测系统的基本原理 水位雨量在线监测系统主要由数据采集单元、数据传输网络、数据处理中心及用户终端四大部分构成,形成了一个完整的闭环系统。 数据采集单元:这是系统的“眼睛”,

认识、理解、分类——acm之搜索

普通搜索方法有两种:1、广度优先搜索;2、深度优先搜索; 更多搜索方法: 3、双向广度优先搜索; 4、启发式搜索(包括A*算法等); 搜索通常会用到的知识点:状态压缩(位压缩,利用hash思想压缩)。

第10章 中断和动态时钟显示

第10章 中断和动态时钟显示 从本章开始,按照书籍的划分,第10章开始就进入保护模式(Protected Mode)部分了,感觉从这里开始难度突然就增加了。 书中介绍了为什么有中断(Interrupt)的设计,中断的几种方式:外部硬件中断、内部中断和软中断。通过中断做了一个会走的时钟和屏幕上输入字符的程序。 我自己理解中断的一些作用: 为了更好的利用处理器的性能。协同快速和慢速设备一起工作

性能分析之MySQL索引实战案例

文章目录 一、前言二、准备三、MySQL索引优化四、MySQL 索引知识回顾五、总结 一、前言 在上一讲性能工具之 JProfiler 简单登录案例分析实战中已经发现SQL没有建立索引问题,本文将一起从代码层去分析为什么没有建立索引? 开源ERP项目地址:https://gitee.com/jishenghua/JSH_ERP 二、准备 打开IDEA找到登录请求资源路径位置

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

Andrej Karpathy最新采访:认知核心模型10亿参数就够了,AI会打破教育不公的僵局

夕小瑶科技说 原创  作者 | 海野 AI圈子的红人,AI大神Andrej Karpathy,曾是OpenAI联合创始人之一,特斯拉AI总监。上一次的动态是官宣创办一家名为 Eureka Labs 的人工智能+教育公司 ,宣布将长期致力于AI原生教育。 近日,Andrej Karpathy接受了No Priors(投资博客)的采访,与硅谷知名投资人 Sara Guo 和 Elad G

hdu1394(线段树点更新的应用)

题意:求一个序列经过一定的操作得到的序列的最小逆序数 这题会用到逆序数的一个性质,在0到n-1这些数字组成的乱序排列,将第一个数字A移到最后一位,得到的逆序数为res-a+(n-a-1) 知道上面的知识点后,可以用暴力来解 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#in