一种市政管网水力模型构建方法
【技术领域】
[0001] 本发明涉及市政工程信息技术领域，具体涉及管网水力模型构建方法。
【背景技术】
[0002] 随着城市建设的发展，供水和供热管网每年都会由于扩建而变化，各个管段阻力 情况也会随之变化，这会造成历史数据辨识的模型与现实管网有一定差异，导致模型的精 度难以满足工业要求。因此，建立反映管网实时状态的水力模型十分必要。
[0003] 国外学者早在上世纪80年代应用最小二乘法等对管网参数进行辨识和校核，进而 建立管网水力模型。在这些研究基础上，Savic应用两种流态模型研究管网模拟问题，结果 表明只依据一个水力工况的运行参数得到计算结果的准确程度不如应用非恒定流模型并 在一段实际内观测管网运行参数条件下得到的结果准确程度高。Lansey提出多水力工况条 件下基于梯度法来解决非线性规划技术来估计管网未知参数的方法，以及一种分步的管网 参数估计方法，该方法可以用来估计管网参数测量中包含不确定的变量条件下的管网参数 估计。Datta和Sridharan提出利用加权最小二乘方法解决部分节点压力和管段流量可观测 条件下的管网参数估计问题。进而，Reddy提出了基于高斯一一牛顿法的加权最小二乘方法 来最小化观测点观测值与计算值之间的偏差，该方法应用于小型管网的效果较好。Bush提 出了根据目标函数中各部分的影响大小来确定加权权值的方法对管网参数进行校核。以上 方法均是利用最小二乘原理处理管网参数校核问题，得到管网阻力特性。但这些方法在处 理观测点数目较少及规模较大的管网阻力数辨识问题时往往存在计算效率不高、易出现不 收敛等情况。
[0004] 随着计算技术的不断发展，以及优化算法在工程领域的成功应用，许多学者将遗 传算法应用于管网阻力数辨识和水力模型建模，并取得了相当的研究成果。在这些研究基 础上，Vitkovsky总结前人学者的研究成果，指出准确估计管道内部阻力的重要性，并且对 其中应用的遗传算法进行改进，以此对管网阻力数和水力模型分别进行估计。Ainola提出 了一种供水管网模型校核方法，该方法在引入一个包含管段长度、平均压力和其他系数经 验方程的基础上，建立了观测压力与计算压力差平方和最小为目标函数的数学模型，并通 过线性优化的方法在计算中改变管段阻力数使计算压力和实测压力之间的偏差最小，该方 法在计算之前首先将管段阻力数进行分类处理。Lingireddy提出了一种考虑需求调整因素 的管网模型优化校核方法。该模型以最小化非线性方程转化成的函数为目标函数，包括线 性与非线性的约束条件，通过遗传算法实现给水管网的模型校核。Kapelan等提出将与水力 工况分析软件Epanet 2相结合的全局优化方法应用于管网最小二乘模型校核问题，该方法 最突出的优点是在一步优化计算过程中可以解决具有一定不确定性的管网模型校核问题。 P. K. Sati ja和A. Kumar结合使用了统计技术和压力测试点分析评估方法，同时利用了人工 神经网络技术来预测管网实际模型，该方法有效的缩短了故障定位和修复时间。但该方法 需要大量观测数据样本用以训练。
[0005] 国内学者在这方面也有一定的研究成果，较具代表性的研究如下。信昆仑通过对 管网中管段阻力数进行合理分组，构建了参数校核的最优化数学模型，利用实数型编码的 遗传算法求解该问题并建立管网水力模型，该方法通过对遗传算法的部分改进措施提高了 辨识计算的效率。段焕丰提出了应用改进遗传算法校核管网模型，并通过一个实例验证了 该方法的效率和精度。张海菲提出一种管网仿真系统及水力计算方法，采用牛顿迭代算法 计算管网未知参数。该方法计算效率相对较高，但计算效率和精度会受到算法初始值的影 响。
[0006] 上述方法的核心是最小化观测值与模拟值偏差，一定程度上可实现管网水力模型 的校核。但由于这类方法，如最小二乘方法、优化算法及人工神经网络方法等计算耗时很 长，并不能满足建立管网实时水力模型的要求。
[0007] 以往方法多依据灵敏度等信息布置观测点，并没有考虑阻力数能观测性这一重要 因素。尚需开展最小观测点数目方面的分析方法研究以及满足能观测性条件下的计算和建 模方法。
[0008] 管网监测需要管网参数(压力和流量)观测数据，但仪表观测值一般会存在误差， 这会影响结果的准确度和水力模型的可信度。需要明确观测误差影响模型准确程度的途 径，才能有可能准确估计误差的影响，并对管网水力模型进行相应调整，减小误差影响。
【发明内容】
[0009] 本发明为了解决现有的管网水力模型构建方法存在的计算耗时长、结果误差难于 估计的问题。
[0010] 一种市政管网水力模型构建方法，包括以下步骤：
[0011]步骤1、通过获取管网系统设计施工图纸等设计资料获取管网系统的管网拓扑结 构、节点管段标高、水栗性能曲线及阀门开度信息，建立管网静态数据库；
[0012] 所述的拓扑结构包括节点和管段信息;涉及节点、管段和阻力数:节点是指管网分 支处的连节点，节点所连接管段流量不相同；管段一般是指一条管道，对于两条管道连接处 实际存在分支但分支处没有出流量，也就是说两条相连接且流量始终相同的管道也视为同 一管段;管段阻力数是综合考虑各管道管径、管长、相对粗糙度的特性参数，反映管段阻力 情况，该参数不随管网水力工况的变化而变化；
[0013] 对节点和管段依次进行编号，并录入计算机;利用无向图表达管网中节点和管段 的连接关系及节点信息，完成管网静态数据库录入;管网静态数据库只需一次输入，不会随 水力工况变化而变化；
[0014] 步骤2、根据管网流量方程和管网能量方程，确定最少观测点数目；
[0015]设管段数目为b，独立节点数目为η;水力工况数目为m，流量观测点数目为X，压力 观测点数目为y;
[0016]根据式（1)和(3)所组成的方程组，使方程数目大于或等于未知量数目jpXmin2b-n ; ymin 2 n-(m-l ) · b/m，Xmin和ymin均取整数，最小取1 ;总体数目为Xmin+ymin，在一定条件下流 量观测点可由相同数目的压力观测点替代，但总数目不能减少；比如:流量观测数量Xmin = 4,压力观测数量y_ = 2,满足上述条件即可，此时总数为6;但是如果流量观测数量为3,压 力观测数量为3,总数也为6,在一定条件下也可以满足条件;据此，可以为新建管网提出观 测点布置要求，也可以对既有管网提出观测点增设的要求；
[0017]管网流量方程如公式⑴所示：
[0018]管网流量方程可以表达各管段流量关系，不可观测管段流量可以按照式（1)由可 观测管段流量表达；
[0020] 式中，Aug为与不可观测流量管段相对应的基本关联矩阵子矩阵;Akg为与可观测流 量管段相对应的基本关联矩阵子矩阵;G u为与不可观测流量值列向量;Gk为与可观测流量值 列向量;Q为管网节点出流量列向量;参数上标(k)为水力工况序号，k取不同数值时，表示不 同水力工况；Gf、Gf1、β~分别表示水力工况序号为k时对应的Gu、Gk、Q;
[0021] 联立不同水力工况对应的管网流量方程，构成流量方程组；
[0022] 不可观测管段流量可以由可观测管段流量表达，那么，各工况条件下的所有管段 流量均可以利用可观测管段流量表达;如果式(1)中矩阵A ug为列满秩矩阵，不可观测流量值 可以由可观测流量值唯一确定;否则，说明流量观测点数目不足，在不能进一步提供流量观 测数据的情况下，只能利用方程组广义逆解近似表达流量结果，其准确程度取决于观测点 数目以及位置；
[0023] 根据下一步矩阵运算的需要，将流量列向量转化为对角矩阵形式，具体见式(2);
[0025]式中，为管网各管段流量值对角矩阵;其中，每个元素与流量列向量中的每个 元素一一对应;diag( ·)为将列向量转化为对角矩阵的矩阵函数；
[0026]利用管网能量方程表达管网阻力数、节点压力和管段流量间的关系;针对m个水力 工况，以各管段阻力数和不可观测节点压力为未知量的方程组，构成管网能量方程组;管网 能量方程组如公式(3)所示：
[0027]
[0028] 式中，A为管网基本关联矩阵，Aup为不可观测压力节点所对应的基本关联矩阵子矩 阵，Akp为可观测压力节点所对应的基本关联矩阵子矩阵;S为管网各个管段阻力数列向量