CN105550405A - 一种市政管网水力模型构建方法 - Google Patents
一种市政管网水力模型构建方法 Download PDFInfo
- Publication number
- CN105550405A CN105550405A CN201510891150.0A CN201510891150A CN105550405A CN 105550405 A CN105550405 A CN 105550405A CN 201510891150 A CN201510891150 A CN 201510891150A CN 105550405 A CN105550405 A CN 105550405A
- Authority
- CN
- China
- Prior art keywords
- pipeline
- flow
- pipe network
- hydraulic
- centerdot
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种市政管网水力模型构建方法,涉及市政工程信息技术领域,具体涉及管网水力模型构建方法。本发明为了解决现有的管网水力模型构建方法存在的计算耗时长、结果误差难于估计的问题。本发明通过获取管网系统设计资料建立管网静态数据库;并根据管网流量方程和管网能量方程,确定最少观测点数目;采集多工况下的节点压力和管段流量观测数据;确定各个管段中的流量方向表达流量方向有向图;然后提取出各个管段阻力数,将管段阻力数数值赋予前述有向图,完成管网水力模型构建并验证管网水力模型可信度,直至得到满足管网分析需要的管网水力模型;本发明适用于市政工程信息技术领域中的管网水力的模型。
Description
技术领域
本发明涉及市政工程信息技术领域,具体涉及管网水力模型构建方法。
背景技术
随着城市建设的发展,供水和供热管网每年都会由于扩建而变化,各个管段阻力情况也会随之变化,这会造成历史数据辨识的模型与现实管网有一定差异,导致模型的精度难以满足工业要求。因此,建立反映管网实时状态的水力模型十分必要。
国外学者早在上世纪80年代应用最小二乘法等对管网参数进行辨识和校核,进而建立管网水力模型。在这些研究基础上,Savic应用两种流态模型研究管网模拟问题,结果表明只依据一个水力工况的运行参数得到计算结果的准确程度不如应用非恒定流模型并在一段实际内观测管网运行参数条件下得到的结果准确程度高。Lansey提出多水力工况条件下基于梯度法来解决非线性规划技术来估计管网未知参数的方法,以及一种分步的管网参数估计方法,该方法可以用来估计管网参数测量中包含不确定的变量条件下的管网参数估计。Datta和Sridharan提出利用加权最小二乘方法解决部分节点压力和管段流量可观测条件下的管网参数估计问题。进而,Reddy提出了基于高斯——牛顿法的加权最小二乘方法来最小化观测点观测值与计算值之间的偏差,该方法应用于小型管网的效果较好。Bush提出了根据目标函数中各部分的影响大小来确定加权权值的方法对管网参数进行校核。以上方法均是利用最小二乘原理处理管网参数校核问题,得到管网阻力特性。但这些方法在处理观测点数目较少及规模较大的管网阻力数辨识问题时往往存在计算效率不高、易出现不收敛等情况。
随着计算技术的不断发展,以及优化算法在工程领域的成功应用,许多学者将遗传算法应用于管网阻力数辨识和水力模型建模,并取得了相当的研究成果。在这些研究基础上,Vitkovsky总结前人学者的研究成果,指出准确估计管道内部阻力的重要性,并且对其中应用的遗传算法进行改进,以此对管网阻力数和水力模型分别进行估计。Ainola提出了一种供水管网模型校核方法,该方法在引入一个包含管段长度、平均压力和其他系数经验方程的基础上,建立了观测压力与计算压力差平方和最小为目标函数的数学模型,并通过线性优化的方法在计算中改变管段阻力数使计算压力和实测压力之间的偏差最小,该方法在计算之前首先将管段阻力数进行分类处理。Lingireddy提出了一种考虑需求调整因素的管网模型优化校核方法。该模型以最小化非线性方程转化成的函数为目标函数,包括线性与非线性的约束条件,通过遗传算法实现给水管网的模型校核。Kapelan等提出将与水力工况分析软件Epanet2相结合的全局优化方法应用于管网最小二乘模型校核问题,该方法最突出的优点是在一步优化计算过程中可以解决具有一定不确定性的管网模型校核问题。P.K.Satija和A.Kumar结合使用了统计技术和压力测试点分析评估方法,同时利用了人工神经网络技术来预测管网实际模型,该方法有效的缩短了故障定位和修复时间。但该方法需要大量观测数据样本用以训练。
国内学者在这方面也有一定的研究成果,较具代表性的研究如下。信昆仑通过对管网中管段阻力数进行合理分组,构建了参数校核的最优化数学模型,利用实数型编码的遗传算法求解该问题并建立管网水力模型,该方法通过对遗传算法的部分改进措施提高了辨识计算的效率。段焕丰提出了应用改进遗传算法校核管网模型,并通过一个实例验证了该方法的效率和精度。张海菲提出一种管网仿真系统及水力计算方法,采用牛顿迭代算法计算管网未知参数。该方法计算效率相对较高,但计算效率和精度会受到算法初始值的影响。
上述方法的核心是最小化观测值与模拟值偏差,一定程度上可实现管网水力模型的校核。但由于这类方法,如最小二乘方法、优化算法及人工神经网络方法等计算耗时很长,并不能满足建立管网实时水力模型的要求。
以往方法多依据灵敏度等信息布置观测点,并没有考虑阻力数能观测性这一重要因素。尚需开展最小观测点数目方面的分析方法研究以及满足能观测性条件下的计算和建模方法。
管网监测需要管网参数(压力和流量)观测数据,但仪表观测值一般会存在误差,这会影响结果的准确度和水力模型的可信度。需要明确观测误差影响模型准确程度的途径,才能有可能准确估计误差的影响,并对管网水力模型进行相应调整,减小误差影响。
发明内容
本发明为了解决现有的管网水力模型构建方法存在的计算耗时长、结果误差难于估计的问题。
一种市政管网水力模型构建方法,包括以下步骤:
步骤1、通过获取管网系统设计施工图纸等设计资料获取管网系统的管网拓扑结构、节点管段标高、水泵性能曲线及阀门开度信息,建立管网静态数据库;
所述的拓扑结构包括节点和管段信息;涉及节点、管段和阻力数:节点是指管网分支处的连节点,节点所连接管段流量不相同;管段一般是指一条管道,对于两条管道连接处实际存在分支但分支处没有出流量,也就是说两条相连接且流量始终相同的管道也视为同一管段;管段阻力数是综合考虑各管道管径、管长、相对粗糙度的特性参数,反映管段阻力情况,该参数不随管网水力工况的变化而变化;
对节点和管段依次进行编号,并录入计算机;利用无向图表达管网中节点和管段的连接关系及节点信息,完成管网静态数据库录入;管网静态数据库只需一次输入,不会随水力工况变化而变化;
步骤2、根据管网流量方程和管网能量方程,确定最少观测点数目;
设管段数目为b,独立节点数目为n;水力工况数目为m,流量观测点数目为x,压力观测点数目为y;
根据式(1)和(3)所组成的方程组,使方程数目大于或等于未知量数目,即xmin≥b-n;ymin≥n-(m-1)·b/m,xmin和ymin均取整数,最小取1;总体数目为xmin+ymin,在一定条件下流量观测点可由相同数目的压力观测点替代,但总数目不能减少;比如:流量观测数量xmin=4,压力观测数量ymin=2,满足上述条件即可,此时总数为6;但是如果流量观测数量为3,压力观测数量为3,总数也为6,在一定条件下也可以满足条件;据此,可以为新建管网提出观测点布置要求,也可以对既有管网提出观测点增设的要求;
管网流量方程如公式(1)所示:
管网流量方程可以表达各管段流量关系,不可观测管段流量可以按照式(1)由可观测管段流量表达;
式中,Aug为与不可观测流量管段相对应的基本关联矩阵子矩阵;Akg为与可观测流量管段相对应的基本关联矩阵子矩阵;Gu为与不可观测流量值列向量;Gk为与可观测流量值列向量;Q为管网节点出流量列向量;参数上标(k)为水力工况序号,k取不同数值时,表示不同水力工况;分别表示水力工况序号为k时对应的Gu、Gk、Q;
联立不同水力工况对应的管网流量方程,构成流量方程组;
不可观测管段流量可以由可观测管段流量表达,那么,各工况条件下的所有管段流量均可以利用可观测管段流量表达;如果式(1)中矩阵Aug为列满秩矩阵,不可观测流量值可以由可观测流量值唯一确定;否则,说明流量观测点数目不足,在不能进一步提供流量观测数据的情况下,只能利用方程组广义逆解近似表达流量结果,其准确程度取决于观测点数目以及位置;
根据下一步矩阵运算的需要,将流量列向量转化为对角矩阵形式,具体见式(2);
式中,为管网各管段流量值对角矩阵;其中,每个元素与流量列向量中的每个元素一一对应;diag(·)为将列向量转化为对角矩阵的矩阵函数;
利用管网能量方程表达管网阻力数、节点压力和管段流量间的关系;针对m个水力工况,以各管段阻力数和不可观测节点压力为未知量的方程组,构成管网能量方程组;管网能量方程组如公式(3)所示:
式中,A为管网基本关联矩阵,Aup为不可观测压力节点所对应的基本关联矩阵子矩阵,Akp为可观测压力节点所对应的基本关联矩阵子矩阵;S为管网各个管段阻力数列向量;Pu为不可观测压力节点压力值列向量,Pk为可观测压力节点压力值列向量;Z为节点标高列向量;DH为管段包含水泵扬程列向量;abs(·)为变量绝对值函数;参数上标(1),(2),…,(m)表示该参数为对应水力工况下的参数;
如果式(3)中的系数矩阵为列满秩矩阵,各管段阻力数可以唯一确定;否则,说明节点压力观测点数目不足,在不能进一步提供节点压力观测数据的情况下,只能利用方程组广义逆解近似表达阻力数,其准确程度同样取决于观测点数目以及位置;
步骤3、采集多工况下的节点压力和管段流量观测数据;
步骤4、管段流量表达:
将一定数目的水力工况条件下采集的管段流量观测值代入式(1),求解其对应方程组的广义逆解,以该解作为各对应水力工况不可观测管段的流量值,依据所得到的流量值确定各个管段中的流量方向,将前述表达管网拓扑结构的无向图补充为可表达流量方向有向图;本方法提出可以唯一确定各个管段流量,即本步骤中的计算结果可以表达不可观测管段流量的准确唯一解;
步骤5、阻力数表达及建立水力模型:
将不同水力工况条件下可采集的节点压力观测值以及步骤4的得到的管段流量值及流量方向代入式(3),求解其对应方程组的广义逆解,从该解中提取出各个管段阻力数;将管段阻力数数值赋予前述有向图,完成管网水力模型构建;该模型可以表达管网拓扑结构、阻力情况和流量信息等;应用本方法可以确定不同水力工况数目条件下的最小压力观测点数目;如果实际可观测压力节点数目大于等于上述最小节点数目,各个管段阻力数可以唯一确定,本步骤中的计算结果理论上可以表达阻力数的唯一解,即准确值;如果可观测压力节点数目少于上述最小节点数目,本步骤中的计算结果实际上是各个管段阻力数的最小二乘解或极小范数解,一般比较接近实际准确值;该模型可用于热力工况分析及管网用水量预测等后续分析中;
步骤6、管网水力模型可信度验证:
管网水力模型构建完成后,利用管网运行参数观测数据中其他水力工况条件下的管网运行参数观测数据,将观测值与依据水力模型计算值进行比对;如果其间总体偏差小于设定值ε,认为步骤5中所构建水力模型可以满足水力工况分析需要;否则,利用其他水力工况观测数据重复步骤3-5,直至得到满足管网分析需要的管网水力模型。
本发明具有以下有益效果:
1、本发明通过分步表达,将阻力数计算这一非线性问题转化为线性方程组求解问题,利用线性化的流量和能量方程分步求解计算。直接求解方程只需一次计算,用以取代以往校核方法,可以降低计算耗时90%以上,可用于大型管网水力模型构建等研究。
2、本发明将多水力工况观测数据作为已知条件,利用方程组广义逆解表达各个管段阻力数有利于分析观测误差对所建立水力模型可信度的影响,结果误差易于估计。
3、本发明能够指出满足管网阻力数能观测性,即可以唯一确定阻力数准确值的最小压力和流量观测点数目,保证建立准确的管网水力模型的情况下削减不必要的测量仪表设置成本。
4、本发明在管网观测点数目不足时,通过计算线性化的流量和能量方程组广义逆解,获得较准确的阻力数估计值,构建较接近实际情况的管网水力模型。
附图说明
图1为本发明的技术流程图。
具体实施方式
具体实施方式一:结合图1说明本实施方式,
一种市政管网水力模型构建方法,包括以下步骤:
步骤1、通过获取管网系统设计施工图纸等设计资料获取管网系统的管网拓扑结构、节点管段标高、水泵性能曲线及阀门开度信息,建立管网静态数据库;
所述的拓扑结构包括节点和管段信息;涉及节点、管段和阻力数:节点是指管网分支处的连节点,节点所连接管段流量不相同;管段一般是指一条管道,对于两条管道连接处实际存在分支但分支处没有出流量,也就是说两条相连接且流量始终相同的管道也视为同一管段;管段阻力数是综合考虑各管道管径、管长、相对粗糙度的特性参数,反映管段阻力情况,该参数不随管网水力工况的变化而变化;
对节点和管段依次进行编号,并录入计算机;利用无向图表达管网中节点和管段的连接关系及节点信息,完成管网静态数据库录入;管网静态数据库只需一次输入,不会随水力工况变化而变化;
步骤2、根据管网流量方程和管网能量方程,确定最少观测点数目;
设管段数目为b,独立节点数目为n;水力工况数目为m,流量观测点数目为x,压力观测点数目为y;
根据式(1)和(3)所组成的方程组,使方程数目大于或等于未知量数目,即xmin≥b-n;ymin≥n-(m-1)·b/m,xmin和ymin均取整数,最小取1;总体数目为xmin+ymin,在一定条件下流量观测点可由相同数目的压力观测点替代,但总数目不能减少;比如:流量观测数量xmin=4,压力观测数量ymin=2,满足上述条件即可,此时总数为6;但是如果流量观测数量为3,压力观测数量为3,总数也为6,在一定条件下也可以满足条件;据此,可以为新建管网提出观测点布置要求,也可以对既有管网提出观测点增设的要求;
管网流量方程如公式(1)所示:
管网流量方程可以表达各管段流量关系,不可观测管段流量可以按照式(1)由可观测管段流量表达;
式中,Aug为与不可观测流量管段相对应的基本关联矩阵子矩阵;Akg为与可观测流量管段相对应的基本关联矩阵子矩阵;Gu为与不可观测流量值列向量;Gk为与可观测流量值列向量;Q为管网节点出流量列向量;参数上标(k)为水力工况序号,k取不同数值时,表示不同水力工况;分别表示水力工况序号为k时对应的Gu、Gk、Q;
联立不同水力工况对应的管网流量方程,构成流量方程组;
不可观测管段流量可以由可观测管段流量表达,那么,各工况条件下的所有管段流量均可以利用可观测管段流量表达;如果式(1)中矩阵Aug为列满秩矩阵,不可观测流量值可以由可观测流量值唯一确定;否则,说明流量观测点数目不足,在不能进一步提供流量观测数据的情况下,只能利用方程组广义逆解近似表达流量结果,其准确程度取决于观测点数目以及位置;
根据下一步矩阵运算的需要,将流量列向量转化为对角矩阵形式,具体见式(2);
式中,为管网各管段流量值对角矩阵;其中,每个元素与流量列向量中的每个元素一一对应;diag(·)为将列向量转化为对角矩阵的矩阵函数;
利用管网能量方程表达管网阻力数、节点压力和管段流量间的关系;针对m个水力工况,以各管段阻力数和不可观测节点压力为未知量的方程组,构成管网能量方程组;管网能量方程组如公式(3)所示:
式中,A为管网基本关联矩阵,Aup为不可观测压力节点所对应的基本关联矩阵子矩阵,Akp为可观测压力节点所对应的基本关联矩阵子矩阵;S为管网各个管段阻力数列向量;Pu为不可观测压力节点压力值列向量,Pk为可观测压力节点压力值列向量;Z为节点标高列向量;DH为管段包含水泵扬程列向量;abs(·)为变量绝对值函数;参数上标(1),(2),…,(m)表示该参数为对应水力工况下的参数;
如果式(3)中的系数矩阵为列满秩矩阵,各管段阻力数可以唯一确定;否则,说明节点压力观测点数目不足,在不能进一步提供节点压力观测数据的情况下,只能利用方程组广义逆解近似表达阻力数,其准确程度同样取决于观测点数目以及位置;
步骤3、采集多工况下的节点压力和管段流量观测数据;
步骤4、管段流量表达:
将一定数目的水力工况条件下采集的管段流量观测值代入式(1),求解其对应方程组的广义逆解,以该解作为各对应水力工况不可观测管段的流量值,依据所得到的流量值确定各个管段中的流量方向,将前述表达管网拓扑结构的无向图补充为可表达流量方向有向图;本方法提出可以唯一确定各个管段流量,即本步骤中的计算结果可以表达不可观测管段流量的准确唯一解;
步骤5、阻力数表达及建立水力模型:
将不同水力工况条件下可采集的节点压力观测值以及步骤4的得到的管段流量值及流量方向代入式(3),求解其对应方程组的广义逆解,从该解中提取出各个管段阻力数;将管段阻力数数值赋予前述有向图,完成管网水力模型构建;该模型可以表达管网拓扑结构、阻力情况和流量信息等;应用本方法可以确定不同水力工况数目条件下的最小压力观测点数目;如果实际可观测压力节点数目大于等于上述最小节点数目,各个管段阻力数可以唯一确定,本步骤中的计算结果理论上可以表达阻力数的唯一解,即准确值;如果可观测压力节点数目少于上述最小节点数目,本步骤中的计算结果实际上是各个管段阻力数的最小二乘解或极小范数解,一般比较接近实际准确值;该模型可用于热力工况分析及管网用水量预测等后续分析中;
步骤6、管网水力模型可信度验证:
管网水力模型构建完成后,利用管网运行参数观测数据中其他水力工况条件下的管网运行参数观测数据,将观测值与依据水力模型计算值进行比对;如果其间总体偏差小于设定值ε,认为步骤5中所构建水力模型可以满足水力工况分析需要;否则,利用其他水力工况观测数据重复步骤3-5,直至得到满足管网分析需要的管网水力模型。
具体实施方式二:
本实施方式的步骤3所述的采集多工况下的节点压力和管段流量观测数据的具体过程如下:
利用消防、夜间及正常等水力工况下管网中节点和管段处设置的压力和流量测试装置获取压力和流量观测数据。
其他步骤和参数与具体实施方式一相同。
具体实施方式三:
本实施方式所述利用消防、夜间及正常水力工况下管网中节点和管段处设置的压力和流量测试装置获取压力和流量观测数据中的夜间为2:00-3:00,即夜间水力工况为2:00-3:00之间的水力工况。
其他步骤和参数与具体实施方式二相同。
具体实施方式四:
本实施方式步骤6中所述的其他水力工况条件下的管网运行参数观测数据包括节点压力和管段流量。
其他步骤和参数与具体实施方式一至三之一相同。
具体实施方式五:
本实施方式步骤4所述的求解方程组的广义逆解的过程是利用Matlab软件工具箱实现的。
其他步骤和参数与具体实施方式一至四之一相同。
具体实施方式六:
本实施方式步骤5所述的求解方程组的广义逆解的过程是利用Matlab软件工具箱实现的。
其他步骤和参数与具体实施方式一至五之一相同。
Claims (6)
1.一种市政管网水力模型构建方法,其特征在于包括以下步骤:
步骤1、通过获取管网系统设计资料获取管网系统的管网拓扑结构、节点管段标高、水泵性能曲线及阀门开度信息,建立管网静态数据库;
所述的拓扑结构包括节点和管段信息;
对节点和管段依次进行编号,并录入计算机;利用无向图表达管网中节点和管段的连接关系及节点信息,完成管网静态数据库录入;
步骤2、根据管网流量方程和管网能量方程,确定最少观测点数目;
设管段数目为b,独立节点数目为n;水力工况数目为m,流量观测点数目为x,压力观测点数目为y;
根据式(1)和(3)所组成的方程组,使方程数目大于或等于未知量数目,即xmin≥b-n;ymin≥n-(m-1)·b/m,xmin和ymin均取整数,最小取1;
管网流量方程如公式(1)所示:
式中,Aug为与不可观测流量管段相对应的基本关联矩阵子矩阵;Akg为与可观测流量管段相对应的基本关联矩阵子矩阵;Gu为与不可观测流量值列向量;Gk为与可观测流量值列向量;Q为管网节点出流量列向量;参数上标(k)为水力工况序号,k取不同数值时,表示不同水力工况;Q(k)分别表示水力工况序号为k时对应的Gu、Gk、Q;
联立不同水力工况对应的管网流量方程,构成流量方程组;
将流量列向量转化为对角矩阵形式,具体见式(2);
式中,为管网各管段流量值对角矩阵;其中,每个元素与流量列向量中的每个元素一一对应;diag(·)为将列向量转化为对角矩阵的矩阵函数;
利用管网能量方程表达管网阻力数、节点压力和管段流量间的关系;针对m个水力工况,以各管段阻力数和不可观测节点压力为未知量的方程组,构成管网能量方程组;管网能量方程组如公式(3)所示:
式中,A为管网基本关联矩阵,Aup为不可观测压力节点所对应的基本关联矩阵子矩阵,Akp为可观测压力节点所对应的基本关联矩阵子矩阵;S为管网各个管段阻力数列向量;Pu为不可观测压力节点压力值列向量,Pk为可观测压力节点压力值列向量;Z为节点标高列向量;DH为管段包含水泵扬程列向量;abs(·)为变量绝对值函数;参数上标(1),(2),…,(m)表示该参数为对应水力工况下的参数;
步骤3、采集多工况下的节点压力和管段流量观测数据;
步骤4、将一定数目的水力工况条件下采集的管段流量观测值代入式(1),求解其对应方程组的广义逆解,以该解作为各对应水力工况不可观测管段的流量值,依据所得到的流量值确定各个管段中的流量方向,将前述表达管网拓扑结构的无向图补充为可表达流量方向有向图;
步骤5、将不同水力工况条件下可采集的节点压力观测值以及步骤4的得到的管段流量值及流量方向代入式(3),求解其对应方程组的广义逆解,从该解中提取出各个管段阻力数;将管段阻力数数值赋予前述有向图,完成管网水力模型构建;
步骤6、管网水力模型构建完成后,利用管网运行参数观测数据中其他水力工况条件下的管网运行参数观测数据,将观测值与依据水力模型计算值进行比对;如果其间总体偏差小于设定值ε,认为步骤5中所构建水力模型可以满足水力工况分析需要;否则,利用其他水力工况观测数据重复步骤3-5,直至得到满足管网分析需要的管网水力模型。
2.根据权利要求1所述的一种市政管网水力模型构建方法,其特征在于步骤3所述的采集多工况下的节点压力和管段流量观测数据的具体过程如下:
利用消防、夜间及正常水力工况下管网中节点和管段处设置的压力和流量测试装置获取压力和流量观测数据。
3.根据权利要求2所述的一种市政管网水力模型构建方法,其特征在于所述利用消防、夜间及正常水力工况下管网中节点和管段处设置的压力和流量测试装置获取压力和流量观测数据中的夜间为2:00-3:00,即夜间水力工况为2:00-3:00之间的水力工况。
4.根据权利要求1、2或3所述的一种市政管网水力模型构建方法,其特征在于步骤6中所述的其他水力工况条件下的管网运行参数观测数据包括节点压力和管段流量。
5.根据权利要求4所述的一种市政管网水力模型构建方法,其特征在于步骤4所述的求解方程组的广义逆解的过程是利用Matlab软件工具箱实现的。
6.根据权利要求5所述的一种市政管网水力模型构建方法,其特征在于步骤5所述的求解方程组的广义逆解的过程是利用Matlab软件工具箱实现的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510891150.0A CN105550405B (zh) | 2015-12-07 | 2015-12-07 | 一种市政管网水力模型构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510891150.0A CN105550405B (zh) | 2015-12-07 | 2015-12-07 | 一种市政管网水力模型构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105550405A true CN105550405A (zh) | 2016-05-04 |
CN105550405B CN105550405B (zh) | 2018-03-20 |
Family
ID=55829594
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510891150.0A Expired - Fee Related CN105550405B (zh) | 2015-12-07 | 2015-12-07 | 一种市政管网水力模型构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105550405B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106647329A (zh) * | 2016-12-16 | 2017-05-10 | 中国电力科学研究院 | 一种供热管网等效建模方法 |
CN106777731A (zh) * | 2016-12-26 | 2017-05-31 | 广东技术师范学院 | 一种燃气管网水力计算方法 |
CN107016622A (zh) * | 2017-03-22 | 2017-08-04 | 浙江大学 | 一种含大用户用水信息的城市供水管网节点需水量反演方法 |
CN108679458A (zh) * | 2018-07-03 | 2018-10-19 | 安徽建筑大学 | 一种供水管网压力相关漏损定位方法 |
CN109214549A (zh) * | 2018-08-01 | 2019-01-15 | 武汉众智鸿图科技有限公司 | 一种基于图论的供水管网辅助dma分区方法和系统 |
CN109325206A (zh) * | 2018-09-10 | 2019-02-12 | 柳创新 | 一种降雨径流模型参数优化方法 |
CN110795810A (zh) * | 2019-11-08 | 2020-02-14 | 苏州智品信息科技有限公司 | 一种在线水力模型生成方法 |
CN110953685A (zh) * | 2019-12-23 | 2020-04-03 | 珠海格力电器股份有限公司 | 通用的空调系统的水力仿真计算方法、存储装置、处理器 |
CN115456432A (zh) * | 2022-09-20 | 2022-12-09 | 厦门市政水务集团有限公司 | 供水管道重要性评估方法及计算机可读存储介质 |
CN116579272A (zh) * | 2023-07-13 | 2023-08-11 | 天津市津安热电有限公司 | 一种基于稳态水力平差计算的管网渗漏定位方法和系统 |
CN117195778A (zh) * | 2023-11-08 | 2023-12-08 | 天津市津安热电有限公司 | 供热管网水力仿真模型参数辨识修正方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101692273A (zh) * | 2009-09-28 | 2010-04-07 | 北京工业大学 | 城市排水管网在线水力模型建模方法 |
CN102063566A (zh) * | 2010-12-15 | 2011-05-18 | 冶金自动化研究设计院 | 一种水力热力耦合仿真模型的多气源蒸汽管网计算系统 |
CN102867090A (zh) * | 2012-09-13 | 2013-01-09 | 冶金自动化研究设计院 | 基于tbb的并行遗传算法蒸汽管网模型自动校准系统 |
CN103294847A (zh) * | 2013-04-12 | 2013-09-11 | 杭州电子科技大学 | 基于水力平差的供水管网模型模糊辨识方法 |
CN103970610A (zh) * | 2014-04-30 | 2014-08-06 | 浙江传媒学院 | 一种供水管网中节点流量的监控方法 |
CN104573143A (zh) * | 2013-10-12 | 2015-04-29 | 苏州热工研究院有限公司 | 一种管网仿真系统及水力计算方法 |
-
2015
- 2015-12-07 CN CN201510891150.0A patent/CN105550405B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101692273A (zh) * | 2009-09-28 | 2010-04-07 | 北京工业大学 | 城市排水管网在线水力模型建模方法 |
CN102063566A (zh) * | 2010-12-15 | 2011-05-18 | 冶金自动化研究设计院 | 一种水力热力耦合仿真模型的多气源蒸汽管网计算系统 |
CN102867090A (zh) * | 2012-09-13 | 2013-01-09 | 冶金自动化研究设计院 | 基于tbb的并行遗传算法蒸汽管网模型自动校准系统 |
CN103294847A (zh) * | 2013-04-12 | 2013-09-11 | 杭州电子科技大学 | 基于水力平差的供水管网模型模糊辨识方法 |
CN104573143A (zh) * | 2013-10-12 | 2015-04-29 | 苏州热工研究院有限公司 | 一种管网仿真系统及水力计算方法 |
CN103970610A (zh) * | 2014-04-30 | 2014-08-06 | 浙江传媒学院 | 一种供水管网中节点流量的监控方法 |
Non-Patent Citations (6)
Title |
---|
LIU Y X, ZOU P H.: "Pipe Resistance Identification of an Experimental Heat-supply Network Based on Minimal Norm and GA Methods", 《INFORMATION AND COMPUTING (ICIC), 2010 THIRD INTERNATIONAL CONFERENCE ON. IEEE》 * |
LIU Y, ZOU P, KANG Y.: "Study on Impedance Identification of Pipes in Heat-Supply Networks", 《INFORMATION ENGINEERING AND COMPUTER SCIENCE, 2009. ICIECS 2009. INTERNATIONAL CONFERENCE ON. IEEE》 * |
刘永鑫,邹平华: "供热管网阻力系数辨识方法比较研究", 《全国暖通空调制冷2010年学术年会论文集》 * |
刘永鑫: "基于矩阵论的供热管网阻力系数辨识研究", 《中国博士学位论文全文数据库工程科技II辑》 * |
张颖: "银川市供水管网水力数学模型的建立及应用", 《宁夏工程技术》 * |
王海,王海鹰等: "供热管网中管段阻力系数的辨识方法", 《计算物理》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106647329A (zh) * | 2016-12-16 | 2017-05-10 | 中国电力科学研究院 | 一种供热管网等效建模方法 |
CN106777731A (zh) * | 2016-12-26 | 2017-05-31 | 广东技术师范学院 | 一种燃气管网水力计算方法 |
CN107016622B (zh) * | 2017-03-22 | 2020-08-04 | 浙江大学 | 一种含大用户用水信息的城市供水管网节点需水量反演方法 |
CN107016622A (zh) * | 2017-03-22 | 2017-08-04 | 浙江大学 | 一种含大用户用水信息的城市供水管网节点需水量反演方法 |
CN108679458A (zh) * | 2018-07-03 | 2018-10-19 | 安徽建筑大学 | 一种供水管网压力相关漏损定位方法 |
CN109214549A (zh) * | 2018-08-01 | 2019-01-15 | 武汉众智鸿图科技有限公司 | 一种基于图论的供水管网辅助dma分区方法和系统 |
CN109214549B (zh) * | 2018-08-01 | 2020-10-20 | 武汉众智鸿图科技有限公司 | 一种基于图论的供水管网辅助dma分区方法和系统 |
CN109325206A (zh) * | 2018-09-10 | 2019-02-12 | 柳创新 | 一种降雨径流模型参数优化方法 |
CN109325206B (zh) * | 2018-09-10 | 2023-03-24 | 柳创新 | 一种降雨径流模型参数优化方法 |
CN110795810A (zh) * | 2019-11-08 | 2020-02-14 | 苏州智品信息科技有限公司 | 一种在线水力模型生成方法 |
CN110953685A (zh) * | 2019-12-23 | 2020-04-03 | 珠海格力电器股份有限公司 | 通用的空调系统的水力仿真计算方法、存储装置、处理器 |
CN110953685B (zh) * | 2019-12-23 | 2020-10-27 | 珠海格力电器股份有限公司 | 通用的空调系统的水力仿真计算方法、存储装置、处理器 |
CN115456432A (zh) * | 2022-09-20 | 2022-12-09 | 厦门市政水务集团有限公司 | 供水管道重要性评估方法及计算机可读存储介质 |
CN116579272A (zh) * | 2023-07-13 | 2023-08-11 | 天津市津安热电有限公司 | 一种基于稳态水力平差计算的管网渗漏定位方法和系统 |
CN116579272B (zh) * | 2023-07-13 | 2023-09-08 | 天津市津安热电有限公司 | 一种基于稳态水力平差计算的管网渗漏定位方法和系统 |
CN117195778A (zh) * | 2023-11-08 | 2023-12-08 | 天津市津安热电有限公司 | 供热管网水力仿真模型参数辨识修正方法 |
CN117195778B (zh) * | 2023-11-08 | 2024-02-20 | 天津市津安热电有限公司 | 供热管网水力仿真模型参数辨识修正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105550405B (zh) | 2018-03-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105550405A (zh) | 一种市政管网水力模型构建方法 | |
Xu et al. | Short‐term inflow forecasting using an artificial neural network model | |
CN109791637B (zh) | 用于布水系统的基于模型的控制的方法和装置 | |
CN104200005A (zh) | 基于神经网络的桥梁损伤识别方法 | |
Sanikhani et al. | Estimation of daily pan evaporation using two different adaptive neuro-fuzzy computing techniques | |
CN104200265A (zh) | 一种改进的基于神经网络的桥梁损伤识别方法 | |
CN111259607B (zh) | 一种河湖过渡区水文边界界定方法 | |
CN107869653A (zh) | 一种管道流量敏感性矩阵漏损检测方法 | |
CN109615860A (zh) | 一种基于非参数贝叶斯框架的信号交叉口状态估计方法 | |
Nasir et al. | Measurement error sensitivity analysis for detecting and locating leak in pipeline using ANN and SVM | |
Dayev et al. | Approximation of the discharge coefficient of differential pressure flowmeters using different soft computing strategies | |
Machekposhti et al. | Modeling climate variables of rivers basin using time series analysis (case study: Karkheh River basin at Iran) | |
Arsene et al. | Mixed simulation-state estimation of water distribution systems based on a least squares loop flows state estimator | |
CN103617563B (zh) | 一种基于地统计空间分析理论的供水管网无监测节点压力确定方法 | |
CN103353295A (zh) | 一种精确预测大坝坝体垂直变形量的方法 | |
CN103279030B (zh) | 基于贝叶斯框架的动态软测量建模方法及装置 | |
Wang et al. | Two-dimension monthly river flow simulation using hierarchical network-copula conditional models | |
Rohaimi et al. | 3 Hours ahead of time flood water level prediction using NNARX structure: Case study pahang | |
CN101609411B (zh) | 一种基于复杂度的嵌入式软件功耗bp神经网络建模方法 | |
Torres et al. | Nonlinear estimation of a power law for the friction in a pipeline | |
Arsene et al. | Modelling and simulation of water systems based on loop equations | |
Preis et al. | Online hydraulic state prediction for water distribution systems | |
CN105302948A (zh) | 单相流体网络建模方法 | |
Sebayang et al. | Dependable flow modeling in upper basin Citarum using multilayer perceptron backpropagation | |
CN110705187B (zh) | 通过最小二乘算法即时在线仪表校验和诊断方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180320 Termination date: 20181207 |