CN107563007A - 一种节点流量和管道阻力系数同时调整的供水管网模型快速校正方法 - Google Patents
一种节点流量和管道阻力系数同时调整的供水管网模型快速校正方法 Download PDFInfo
- Publication number
- CN107563007A CN107563007A CN201710667800.2A CN201710667800A CN107563007A CN 107563007 A CN107563007 A CN 107563007A CN 201710667800 A CN201710667800 A CN 201710667800A CN 107563007 A CN107563007 A CN 107563007A
- Authority
- CN
- China
- Prior art keywords
- flow
- node
- resistance coefficient
- representing
- pipeline
- 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
Abstract
本发明提供一种节点流量和管道阻力系数同时调整的供水管网模型快速校正方法,包括:管道阻力系数分组并初始化、节点流量分组并初始化、计算监测点观测值与计算值残差、计算管网监测点灵敏度矩阵、计算参数调整值、更新参数值。本发明具有原创性,是对供水管网模型校正研究领域的一个重要补充,为供水管网的智能化在线监控与管理提供了重要的技术支撑,具有很好的推广和实际工程应用价值。
Description
技术领域
本发明涉及市政工程和城市供水管网领域。
背景技术
近年来,随着人口增长以及城市化发展,我国城市供水管网系统结构日趋复杂,导致供水管网运行管理困难,运行效率较低、漏损严重、爆管频发、供水能耗偏高。为解决这些问题,有必要建立供水管网水力模型,借助水力模型分析,可实现管网运行的智能化管理,并对提高管网运行效率、降低管网漏损、节约泵站能耗具有重要的意义。
管网水力模型在应用之前,需要对模型参数进行校正,即根据监测点实际观察数据调整模型节点流量和管道阻力系数,使得监测点的观测值与模型计算值误差满足模型应用的精度要求。针对模型参数校正问题,国内外广大科研工作者主要基于优化方法(如:遗传算法、粒子群算法)调整模型参数,这些优化方法主要的缺点是计算时间过长,往往需要几小时甚至几天,完全不能满足管网模型在线实时校正的要求。
为解决优化方法的不足,有些科研工作者提出了节点流量反演算法,该算法假定管道阻力系数不变,通过数值计算来推算管网节点流量,以实现管网的校正。这种方法尽管解决了优化方法的计算效率问题,但在校正过程中没有考虑管道阻力系数。实际管网模型的节点流量和管道阻力系数都会影响模型的精度,它们之间是相互关联的,因此,固定一种参数不变,而调节另外一种参数的方法存在不足。而且,目前的反演算法主要是提高了管网主干管流量的模拟精度,而对管网节点流量的校正存在精度较差的缺点,如节点流量校正值存在明显过大或过小的不合理现象。
发明内容
本发明所要解决的技术问题:提供一种节点流量和管道阻力系数同时调整的供水管网模型快速校正方法,以实现管网在线实时模拟,进而实现管网的运行智能化管理,提高了管网运行效率、降低了管网漏损和节约了泵站能耗。
本方法的总体核心技术方案如下:
(1)管道阻力系数分组并初始化:对于一个给定的管道数目(不包含阀门)为np的供水管网,根据管道特征(管材、管龄、管径)对管道阻力系数分组,总数为lC,则管道阻力系数分组矩阵GC大小为np×lC,GC中的元素表示为:
分组之后,每组中所有的管道具有相同的阻力系数,管道阻力系数初始值根据当地管道阻力系数历史资料确定;
(2)节点流量分组并初始化:对于一个给定的节点数目为nn的供水管网,先将单独计量的用水量分配至相关节点,剩余水量根据管线长度按比例分配至对应的节点,然后根据节点用水特性对节点流量进行分组,分组总数为lq,则节点流量分组矩阵Gq大小为nn×lq,Gq中的元素表示为:
上式中,表示节点i的初始分配流量;表示分组j中所有节点的初始分配流量总和,每组中所有节点都具有相同的用水量变化特征;
(3)计算监测点观测值与计算值残差:运行管网水力模拟,计算第k次迭代时(k=1,2,…K),K为设定的最大迭代次数,实际迭代次数根据设定的精度自动确定。压力监测点观测值与计算值残差
流量监测点观测值与计算值残差
式中,NH、NQ为压力、流量监测点的数目;hi、表示第i个压力监测点的观测值、计算值;qj、表示第j个流量监测点的观测值、计算值;T表示向量的转置。
(4)计算管网监测点灵敏度矩阵:根据公式1-1、1-2,计算第k次迭代时(k=1,2,…K),K为设定的最大迭代次数,实际迭代次数根据设定的精度自动确定。监测点偏导数矩阵HC(Ck,qk)、QC(Ck,qk)、Hq(Ck,qk)、Qq(Ck,qk),计算公式为
式中,表示压力监测点压力对管道阻力系数、节点流量的偏导数矩阵;表示流量监测点流量对管道阻力系数、节点流量的偏导数矩阵;GC表示管道阻力系数分组矩阵(公式1-1);Gq表示节点流量分组矩阵(公式1-2)。
(5)计算参数调整值:根据公式1-1~公式1-8,计算第k次迭代时,(k=1,2,…K),K为设定的最大迭代次数,实际迭代次数根据设定的精度自动确定。管道阻力系数调整值△Ck和节点流量调整值△qk,
式中,
表示权重矩阵,对监测点压力和流量数据进行无量纲化处理,矩阵W元素表示压力监测点i的权重,表示流量监测点j的权重,hi表示第i个压力监测点的观测值;qj表示第j个流量监测点的观测值;NH、NQ为压力、流量监测点的数目。
(6)更新参数值:根据公式1-1、公式1-2以及公式1-9计算得到的△Ck、△qk,更新下一次迭代的管道阻力系数及节点流量
然后,判断每根管道阻力系数和每个节点需水量是否在设定的边界范围内
式中,表示管道i(i=1,...,np,np为管道数目)第k+1次迭代的阻力系数;表示节点j(j=1,...,nn,nn为节点数目)第k+1次迭代的流量;和分别表示管道i阻力系数的最小和最大值,即表示管道i阻力系数初始值,c表示管道i阻力系数调整范围,一般取c=15%;和分别表示节点j流量的最小和最大值,即 表示节点j流量初始值,d表示节点j流量调整范围,一般取d=20%。公式1-11和1-12是为了使节点流量和管道阻力系数在合理范围之内,避免参数值过大或过小,进而使校正结果更具有工程意义。
为解决传统模型参数校正方法的瓶颈问题,本发明首次提出了一种节点流量和管道阻力系数同时调整的管网模型快速校正方法。本发明与现有模型参数校正技术相比具有以下优点:①现有模型参数校正技术主要基于优化方法,计算速度很慢,根本无法实现供水管网的在线校正,而本发明方法采用数值计算方法,计算速度快,尤其在解决大型复杂的实际供水管网校正问题时,优势尤为明显。②现有的数值校正方法通常只考虑节点流量,而本发明方法对节点流量和管道阻力系数同时调整,更符合工程实践。③本发明方法能对节点流量和管道阻力系数的校正值进行自动控制,以使其在合理范围内,解决了目前数值校正方法普遍存在校正值过大或过小的不合理问题。比如,目前数值方法校正后的管道阻力系数可能为C=30,显然不符合工程实际,而本发明方法可以将管道阻力系数校正值控制在合理范围内(比如:C=80~110)。本发明具有原创性,是对供水管网模型校正研究领域的一个重要补充,为供水管网的智能化在线监控与管理提供了重要的技术支撑,具有很好的推广和实际工程应用价值。
附图说明
图1是本发明管网模型参数快速校正方法总流程图。
图2是某城市供水管网示意图。
图3是16个临时压力监测点计算值与模拟值对比结果。
具体实施方式
参见图1,本发明的具体实施步骤如下:
(1)输入相关参数:压力监测点观测值h=[h1,h1,...,hNH](NH为压力监测点的数目,hi表示第i个压力监测点的观测值)、流量监测点观测值q=[q1,q2,...,qNQ](NQ为流量监测点的数目,qj表示第j个流量监测点的观测值)、管道阻力系数调整误差ε1、节点流量调整误差ε2、最大迭代次数K、管道阻力系数调整范围c、节点流量调整范围d。
(2)管道阻力系数分组并初始化:首先,根据管道特征(管材、管龄、管径)对管道阻力系数分组,根据公式1-1计算管道阻力系数分组矩阵GC,根据当地历史资料确定管道阻力系数初始值Ck。
(3)节点流量分组并初始化:首先,初始化分配节点流量qk,先将单独计量的大用户和居民用水量按最近节点分配至对应的节点,剩余水量(包括未单独计量居民用水量和漏损量)根据管线长度按比例分配至对应的节点。然后,根据节点用水特性对节点流量进行分组,根据公式1-2计算节点流量分组矩阵Gq。
(4)计算监测点观测值与计算值残差:运行管网水力模拟,根据公式1-3、公式1-4计算第k次迭代时,监测点观测值与计算值残差△Hk、△Qk。
(5)计算管网监测点灵敏度矩阵:根据公式1-5~公式1-8,计算第k次迭代时,监测点灵敏度矩阵HC(Ck,qk)、QC(Ck,qk)、Hq(Ck,qk)、Qq(Ck,qk)。
(6)计算参数调整值:根据公式1-9,计算第k次迭代时,管道阻力系数调整值△Ck和节点流量调整值△qk。然后,判断参数修正值是否在设定的误差范围内,即||△Ck||2<ε1且||△qk||2<ε2,若满足条件或者迭代次数达到设定值K,则停止迭代,输出最终计算结果;若不满足要求,则进入过程(7)。
(7)更新参数值:根据公式1-10,更新下一次迭代的管道阻力系数△Ck+1及节点流量△qk+1。然后,根据公式1-11判断管道阻力系数和节点流量是否超过调整范围,若大于设定的最大值,取其设定的最大值;若小于设定的最小值,取其设定的最小值。
(8)将计算得到的管道阻力系数Ck+1及节点流量qk+1赋予模型中对应的管道及节点,更新管网模型,然后跳转至过程(4)。过程(4)~(8)重复执行,直至满足终止条件。
将本发明的方法应用到某城市供水管网(如图2所示)。为65000人每天提供大约15万吨的饮用水,服务面积超过700平方公里,对应的管网模型包含2189个需水节点和2416根管道,模型现有30个永久压力监测点、16个临时压力监测点和24个永久流量监测点。基于30个永久压力监测点和24个永久流量监测点观测数据,使用本发明提出的校正方法进行计算,经过17次迭代计算后,得出管道阻力系数和节点流量的校正结果,整个计算过程耗时19.94秒。校正之后的模型在30个永久压力监测点和24个永久流量监测点处模拟值与观测值误差极小,压力误差最大仅0.01m,流量误差最大仅0.73%。为了检验模型其余节点压力精度,图3给出了16个临时压力监测点计算值与模拟值对比结果。
应用本发明方法进行模型参数校正之后,16个临时压力监测点现场观测值与模型模拟值之间误差较小,都在模型规定的误差精度范围内(98%压力监测点在±2m范围内,50%压力监测点在±1m范围内)。除此之外,使用本发明方法计算耗时仅十几秒,极大地提高了模型参数校正效率,实现了管网在线实时校正,使模型更好的为供水企业服务,在节约能源、降低产销差的同时,提高了供水的服务水平。
Claims (1)
1.一种节点流量和管道阻力系数同时调整的供水管网模型快速校正方法,其特征在于:
该方法包括步骤:
(1)管道阻力系数分组并初始化:对于一个给定的管道数目(不包含阀门)为np的供水管网,根据管道特征(管材、管龄、管径)对管道阻力系数分组,总数为lC,则管道阻力系数分组矩阵GC大小为np×lC,GC中的元素表示为:
分组之后,每组中所有的管道具有相同的阻力系数,管道阻力系数初始值根据当地管道阻力系数历史资料确定;
(2)节点流量分组并初始化:对于一个给定的节点数目为nn的供水管网,先将单独计量的用水量分配至相关节点,剩余水量根据管线长度按比例分配至对应的节点,然后根据节点用水特性对节点流量进行分组,分组总数为lq,则节点流量分组矩阵Gq大小为nn×lq,Gq中的元素表示为:
上式中,表示节点i的初始分配流量;表示分组j中所有节点的初始分配流量总和,每组中所有节点都具有相同的用水量变化特征;
(3)计算监测点观测值与计算值残差:运行管网水力模拟,计算第k次迭代时(k=1,2,…K),K为设定的最大迭代次数,实际迭代次数根据设定的精度自动确定,压力监测点观测值与计算值残差:
<mrow>
<msub>
<mi>&Delta;H</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>H</mi>
<mn>1</mn>
<mi>k</mi>
</msubsup>
<mo>,</mo>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>H</mi>
<mn>2</mn>
<mi>k</mi>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>h</mi>
<mrow>
<mi>N</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>H</mi>
<mrow>
<mi>N</mi>
<mi>H</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mn>3</mn>
<mo>,</mo>
</mrow>
流量监测点观测值与计算值残差:
<mrow>
<msub>
<mi>&Delta;Q</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>Q</mi>
<mn>1</mn>
<mi>k</mi>
</msubsup>
<mo>,</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>Q</mi>
<mn>2</mn>
<mi>k</mi>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>N</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>N</mi>
<mi>Q</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mn>4</mn>
<mo>,</mo>
</mrow>
式中,NH、NQ为压力、流量监测点的数目;hi、表示第i个压力监测点的观测值、计算值;qj、表示第j个流量监测点的观测值、计算值;T表示向量的转置;
(4)计算管网监测点灵敏度矩阵:根据公式1-1、1-2,计算第k次迭代时(k=1,2,…K),K为设定的最大迭代次数,实际迭代次数根据设定的精度自动确定,监测点偏导数矩阵HC(Ck,qk)、QC(Ck,qk)、Hq(Ck,qk)、Qq(Ck,qk),计算公式为
<mrow>
<msub>
<mi>H</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>H</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>C</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>o</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>G</mi>
<mi>C</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mn>5</mn>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>H</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>H</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>q</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>o</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>G</mi>
<mi>q</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mn>6</mn>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>Q</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>Q</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>C</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>o</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>G</mi>
<mi>C</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mn>7</mn>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>Q</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>Q</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>q</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>o</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>G</mi>
<mi>q</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mn>8</mn>
<mo>,</mo>
</mrow>
式中,表示压力监测点压力对管道阻力系数、节点流量的偏导数矩阵;表示流量监测点流量对管道阻力系数、节点流量的偏导数矩阵;GC表示管道阻力系数分组矩阵(公式1-1);Gq表示节点流量分组矩阵(公式1-2);
(5)计算参数调整值:根据公式1-1~公式1-8,计算第k次迭代时(k=1,2,…K),K为设定的最大迭代次数,实际迭代次数根据设定的精度自动确定,管道阻力系数调整值ΔCk和节点流量调整值Δqk,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;C</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>B</mi>
<mn>1</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>B</mi>
<mn>2</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>B</mi>
<mn>1</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mi>X</mi>
<mo>-</mo>
<msubsup>
<mi>B</mi>
<mn>2</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mi>Y</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;q</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>A</mi>
<mn>1</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msubsup>
<mi>A</mi>
<mn>2</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>A</mi>
<mn>1</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mi>X</mi>
<mo>-</mo>
<msubsup>
<mi>A</mi>
<mn>2</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mi>Y</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mn>9</mn>
<mo>,</mo>
</mrow>
式中,
<mrow>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mi>W</mi>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mi>W</mi>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mi>W</mi>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mi>W</mi>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
<mrow>
<mi>X</mi>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mi>W</mi>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;H</mi>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;Q</mi>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>Y</mi>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mi>W</mi>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;H</mi>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;Q</mi>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
表示权重矩阵,对监测点压力和流量数据进行无量纲化处理,矩阵W元素表示压力监测点i的权重,表示流量监测点j的权重,hi表示第i个压力监测点的观测值;qj表示第j个流量监测点的观测值;NH、NQ为压力、流量监测点的数目,Ck表示管道阻力系数,qk表示节点流量。
(6)更新参数值:根据公式1-1、公式1-2以及公式1-9计算得到的ΔCk、Δqk,更新下一次迭代的管道阻力系数及节点流量,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>G</mi>
<mi>C</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>&Delta;C</mi>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>q</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>G</mi>
<mi>q</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>&Delta;q</mi>
<mi>k</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mn>10</mn>
<mo>,</mo>
</mrow>
然后,判断每根管道阻力系数和每个节点需水量是否在设定的边界范围内,
式中,表示管道i(i=1,...,np,np为管道数目)第k+1次迭代的阻力系数;表示节点j(j=1,...,nn,nn为节点数目)第k+1次迭代的流量;和分别表示管道i阻力系数的最小和最大值,即表示管道i阻力系数初始值,c表示管道i阻力系数调整范围,一般取c=15%;和分别表示节点j流量的最小和最大值,即 表示节点j流量初始值,d表示节点j流量调整范围,一般取d=20%。公式1-11和1-12是为了使节点流量和管道阻力系数在合理范围之内,避免参数值过大或过小,进而使校正结果更具有工程意义。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710667800.2A CN107563007B (zh) | 2017-08-07 | 2017-08-07 | 一种供水管网模型快速校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710667800.2A CN107563007B (zh) | 2017-08-07 | 2017-08-07 | 一种供水管网模型快速校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107563007A true CN107563007A (zh) | 2018-01-09 |
CN107563007B CN107563007B (zh) | 2019-08-27 |
Family
ID=60973989
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710667800.2A Active CN107563007B (zh) | 2017-08-07 | 2017-08-07 | 一种供水管网模型快速校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107563007B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108843977A (zh) * | 2018-06-28 | 2018-11-20 | 武汉新烽光电股份有限公司 | 供水管网实时漏损分析方法、设备、系统及存储介质 |
CN108898512A (zh) * | 2018-07-27 | 2018-11-27 | 苏州市自来水有限公司 | 基于bp神经网络的城市供水管网模型校核方法 |
CN109117985A (zh) * | 2018-07-13 | 2019-01-01 | 上海市水利工程设计研究院有限公司 | 一种基于矩阵与遗传算法的管网压力监测点优化布置方法 |
CN109185708A (zh) * | 2018-11-05 | 2019-01-11 | 泰华智慧产业集团股份有限公司 | 一种供水管网监测方法及系统 |
CN112113146A (zh) * | 2020-09-18 | 2020-12-22 | 同济大学 | 供水管网管道粗糙系数和节点需水量同步自适应校核方法 |
CN112131688A (zh) * | 2019-06-05 | 2020-12-25 | 台湾优化水务股份有限公司 | 供水管网的水力模型校正方法、系统及记录媒体 |
CN112182984A (zh) * | 2020-08-18 | 2021-01-05 | 浙江大学 | 一种基于供水物联网数据同化的污水管网实时模拟方法 |
CN112241767A (zh) * | 2020-11-03 | 2021-01-19 | 广州杰赛科技股份有限公司 | 管网水压监测节点布置方法、装置、设备及存储介质 |
CN115329692A (zh) * | 2022-08-02 | 2022-11-11 | 中国石油大学(北京) | 一种成品油管道水力参数在线校准方法及系统 |
CN115994487A (zh) * | 2023-01-17 | 2023-04-21 | 安徽建筑大学 | 一种基于改进灰狼优化算法的供水管网漏损定位方法 |
CN116150930A (zh) * | 2023-04-20 | 2023-05-23 | 天津智云水务科技有限公司 | 基于3d地图信息的供水管网水力模型节点流量初始化方法 |
CN116187013A (zh) * | 2022-12-29 | 2023-05-30 | 武汉博视曼医疗科技有限公司 | 气腹机自适应负载控制方法、装置、设备及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102867090A (zh) * | 2012-09-13 | 2013-01-09 | 冶金自动化研究设计院 | 基于tbb的并行遗传算法蒸汽管网模型自动校准系统 |
CN102968529A (zh) * | 2012-11-14 | 2013-03-13 | 上海城市水资源开发利用国家工程中心有限公司 | 一种供水管网模型计算结果不确定性区间的量化方法 |
CN103839190A (zh) * | 2014-02-19 | 2014-06-04 | 清华大学深圳研究生院 | 基于压力监测的管网节点流量计量与调度方法 |
CN103983313A (zh) * | 2014-05-30 | 2014-08-13 | 重庆长厦安基建筑设计有限公司 | 一种管网阻力系数确定方法及管网风量测量方法 |
CN106096243A (zh) * | 2016-06-02 | 2016-11-09 | 东北电力大学 | 一种基于伴随矩阵的供水管网泄漏故障反向寻源方法 |
CN106870955A (zh) * | 2017-03-22 | 2017-06-20 | 浙江大学 | 服务于供水管网节点需水量反演的管网监测点优化布置方法 |
-
2017
- 2017-08-07 CN CN201710667800.2A patent/CN107563007B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102867090A (zh) * | 2012-09-13 | 2013-01-09 | 冶金自动化研究设计院 | 基于tbb的并行遗传算法蒸汽管网模型自动校准系统 |
CN102968529A (zh) * | 2012-11-14 | 2013-03-13 | 上海城市水资源开发利用国家工程中心有限公司 | 一种供水管网模型计算结果不确定性区间的量化方法 |
CN103839190A (zh) * | 2014-02-19 | 2014-06-04 | 清华大学深圳研究生院 | 基于压力监测的管网节点流量计量与调度方法 |
CN103983313A (zh) * | 2014-05-30 | 2014-08-13 | 重庆长厦安基建筑设计有限公司 | 一种管网阻力系数确定方法及管网风量测量方法 |
CN106096243A (zh) * | 2016-06-02 | 2016-11-09 | 东北电力大学 | 一种基于伴随矩阵的供水管网泄漏故障反向寻源方法 |
CN106870955A (zh) * | 2017-03-22 | 2017-06-20 | 浙江大学 | 服务于供水管网节点需水量反演的管网监测点优化布置方法 |
Non-Patent Citations (7)
Title |
---|
PILLER O, GILBERT D, VAN ZYL J E.: "Dual calibration for coupled flow and transport", 《WATER DISTRIBUTION SYSTEM ANALYSIS 2010 – WDSA2010》 * |
Q. ZHANG, Z. WU, M. ZHAO, J. QI, Y. HANG AND H, ZHAO.: "Auto-Partitioning of Water", 《JOURNAL OF WATER RESOURCES PLANNING AND MANAGEMENT》 * |
张土乔,许刚等: "给水管网管道摩阻校正方法研究", 《浙江大学学报(工学版)》 * |
张土乔,邵煜: "城镇供水管网漏损监测与控制技术及应用", 《中国环境管理》 * |
张永华, 张土乔, 许刚: "蚁群改进算法用于管道摩阻校正试验点的优化探讨", 《中国农村水利水电》 * |
张清周,黄源等: "给水管网新增压力监测点优化布置方法", 《给水排水》 * |
程伟平、俞亭超、许刚: "广州市大规模供水管网实时水力模型的构建与应用", 《城镇供水》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108843977A (zh) * | 2018-06-28 | 2018-11-20 | 武汉新烽光电股份有限公司 | 供水管网实时漏损分析方法、设备、系统及存储介质 |
CN109117985A (zh) * | 2018-07-13 | 2019-01-01 | 上海市水利工程设计研究院有限公司 | 一种基于矩阵与遗传算法的管网压力监测点优化布置方法 |
CN109117985B (zh) * | 2018-07-13 | 2021-09-17 | 上海市水利工程设计研究院有限公司 | 一种基于矩阵与遗传算法的管网压力监测点优化布置方法 |
CN108898512A (zh) * | 2018-07-27 | 2018-11-27 | 苏州市自来水有限公司 | 基于bp神经网络的城市供水管网模型校核方法 |
CN109185708A (zh) * | 2018-11-05 | 2019-01-11 | 泰华智慧产业集团股份有限公司 | 一种供水管网监测方法及系统 |
CN112131688A (zh) * | 2019-06-05 | 2020-12-25 | 台湾优化水务股份有限公司 | 供水管网的水力模型校正方法、系统及记录媒体 |
WO2022036820A1 (zh) * | 2020-08-18 | 2022-02-24 | 浙江大学 | 一种基于供水物联网数据同化的污水管网实时模拟方法 |
CN112182984A (zh) * | 2020-08-18 | 2021-01-05 | 浙江大学 | 一种基于供水物联网数据同化的污水管网实时模拟方法 |
CN112113146B (zh) * | 2020-09-18 | 2021-09-03 | 同济大学 | 供水管网管道粗糙系数和节点需水量同步自适应校核方法 |
CN112113146A (zh) * | 2020-09-18 | 2020-12-22 | 同济大学 | 供水管网管道粗糙系数和节点需水量同步自适应校核方法 |
CN112241767A (zh) * | 2020-11-03 | 2021-01-19 | 广州杰赛科技股份有限公司 | 管网水压监测节点布置方法、装置、设备及存储介质 |
CN115329692A (zh) * | 2022-08-02 | 2022-11-11 | 中国石油大学(北京) | 一种成品油管道水力参数在线校准方法及系统 |
CN115329692B (zh) * | 2022-08-02 | 2024-03-19 | 中国石油大学(北京) | 一种成品油管道水力参数在线校准方法及系统 |
CN116187013A (zh) * | 2022-12-29 | 2023-05-30 | 武汉博视曼医疗科技有限公司 | 气腹机自适应负载控制方法、装置、设备及存储介质 |
CN115994487A (zh) * | 2023-01-17 | 2023-04-21 | 安徽建筑大学 | 一种基于改进灰狼优化算法的供水管网漏损定位方法 |
CN115994487B (zh) * | 2023-01-17 | 2023-10-03 | 安徽建筑大学 | 一种基于改进灰狼优化算法的供水管网漏损定位方法 |
CN116150930A (zh) * | 2023-04-20 | 2023-05-23 | 天津智云水务科技有限公司 | 基于3d地图信息的供水管网水力模型节点流量初始化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107563007B (zh) | 2019-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107563007A (zh) | 一种节点流量和管道阻力系数同时调整的供水管网模型快速校正方法 | |
CN112128841B (zh) | 一种基于负荷预测和室温反馈修正的全网平衡调节方法 | |
CN105910169B (zh) | 基于机理模型预测控制的城市供热系统热网调节方法及系统 | |
CN107016622B (zh) | 一种含大用户用水信息的城市供水管网节点需水量反演方法 | |
CN106842947B (zh) | 一种城市排水管网的安全运行控制方法 | |
CN112113146B (zh) | 供水管网管道粗糙系数和节点需水量同步自适应校核方法 | |
KR102031714B1 (ko) | 수리해석 기반의 상수도관망 누수탐지 시스템 | |
CN112128842A (zh) | 一种供热管网水力平衡快速调节方法 | |
CN108197725B (zh) | 基于需水量先验信息的供水管网节点需水量校核方法 | |
CN111475913A (zh) | 一种蒸汽动力系统的运行操作优化方法及系统 | |
CN109376925A (zh) | 供水管网节点流量动态自适应优化方法 | |
CN112417662A (zh) | 一种实现集中供热管网系统动态水力优化的方法 | |
CN109558673A (zh) | 一种天然气输配系统适应性分析研究方法 | |
CN103258235A (zh) | 一种基于改进人工鱼群算法的供水管网改扩建优化方法 | |
CN208779553U (zh) | 管网平衡调节系统 | |
CN110244576A (zh) | 一种基于计算机仿真调节的二网平衡方法 | |
CN107559945A (zh) | 一种水力自调节供热系统及自调节方法 | |
CN108800303B (zh) | 一种基于模式识别的精细化节能供热方法 | |
CN113361214A (zh) | 一种基于水位流量数据的明渠控制模型参数辨识方法 | |
CN113065190B (zh) | 一种基于不确定性的居住建筑供热量计算方法 | |
CN110969342A (zh) | 一种平衡小库容水力发电安全与防洪风险管控的方法 | |
CN111859562B (zh) | 一种基于动态阻力特性的压力无关智能调节阀设计方法 | |
CN112381366A (zh) | 一种渐进式可重构的供水管网分区方法 | |
CN117432941A (zh) | 一种水厂供水压力的优化调整方法及系统 | |
CN109782594B (zh) | 一种水务系统安全供水控制器设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |