首页 > 专利 > 杭州电子科技大学 > 基于SCADA数据和实时模型的城市供水管网爆管在线定位方法专利详情

基于SCADA数据和实时模型的城市供水管网爆管在线定位方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2016-07-15
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2016-12-21
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2019-07-12
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2036-07-15
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN201610557790.2 申请日 2016-07-15
公开/公告号 CN106157173B 公开/公告日 2019-07-12
授权日 2019-07-12 预估到期日 2036-07-15
申请年 2016年 公开/公告年 2019年
缴费截止日
分类号 F17D5/06 主分类号 F17D5/06
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 0
权利要求数量 1 非专利引证数量 0
引用专利数量 0 被引证专利数量 0
非专利引证
引用专利 被引证专利
专利权维持 6 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 徐哲、李玉全、蔡华强、熊晓锋、何必仕 第一发明人 徐哲
地址 浙江省杭州市下沙高教园区2号大街 邮编 310018
申请人数量 1 发明人数量 5
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
浙江杭州金通专利事务所有限公司 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
王佳健
摘要
本发明公开了一种基于SCADA数据和实时模型的城市供水管网爆管在线定位方法。本发明首先监测压力波谷,确定疑似爆管时间边界和校核测点,然后根据每次爆管发生后测点的流量异常率变化,选择测点当中异常率变化明显的测点,按照它们所在位置,划定矩形区域,作为实时优化校核模型搜索爆管点的初始区域。最后目标函数以连续时间内压力测点的压力模拟值和实际值、流量测点的模拟值与实际值之间的差值绝对值之和作为评价标准,迭代求解最小值,获得最可能的爆管节点和爆管流量。本发明能实现对爆管等异常事件的高精度快速侦测定位,降低爆管对供水管网系统造成的危害。
  • 摘要附图
    基于SCADA数据和实时模型的城市供水管网爆管在线定位方法
  • 说明书附图:图1
    基于SCADA数据和实时模型的城市供水管网爆管在线定位方法
  • 说明书附图:图2
    基于SCADA数据和实时模型的城市供水管网爆管在线定位方法
  • 说明书附图:图3
    基于SCADA数据和实时模型的城市供水管网爆管在线定位方法
  • 说明书附图:图4
    基于SCADA数据和实时模型的城市供水管网爆管在线定位方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2019-07-12 授权
2 2016-12-21 实质审查的生效 IPC(主分类): G06Q 50/06 专利申请号: 201610557790.2 申请日: 2016.07.15
3 2016-11-23 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.基于SCADA数据和实时模型的城市供水管网爆管在线定位方法,其特征在于该方法包括以下步骤:
步骤1、监测压力波谷,确定疑似爆管时间边界和校核测点
根据各测压点在爆管发生前后的压力变化的先后顺序和压力波谷,得出测点压力因爆管发生变化的起始时间、截止时间,从而获得用于实时模型优化校核的起始时间t0,持续时间T以及有压力值反馈的时刻k,选择负压波变化前几位的测点作为校核测点,设其个数为m;
步骤2、初步定位
爆管发生前后,测点与水源压差的比值反映了爆管发生后测点附近的流量异常情况,用公式表示为:
其中, 表示管网中各测点i在爆管发生后k时刻时附近区域的流量异常率,t0表示爆管发生起始时刻, 表示爆管发生后k时刻测点i与水源水头的压差, 表示t0时刻测点i与水源水头的压差;
根据每次爆管发生后测点的流量异常率变化,选择测点当中异常率变化前几位的测点,按照它们所在位置,划定矩形区域,作为实时优化校核模型搜索爆管点的初始区域;
步骤3、精确定位
目标函数以连续时间内压力测点的压力模拟值和实际值、流量测点的模拟值与实际值之间的差值绝对值之和作为评价标准,迭代求解最小值,获得最可能的爆管节点和爆管流量;具体为:
其中, 为压力测点i在t时刻的优化模拟值, 为压力测点SCADA监测值;
为测流管线j的优化模拟值, 为测流管线j的SCADA监测值,T为持续时间,m为由步骤1选定的压力测点数目,n为供水管网中管线流量测点数目,αi、βj是分配的权重系数;
通过模型模拟某节点的流量大幅度变化来模拟爆管现象,以节点位置以及节点漏失量作为优化过程的自变量,通过遗传算法优化求解;
设定遗传算法参数,包括种群规模、进化代数、选择、交叉、变异方法、交叉概率、变异概率;确定染色体编码方式,以空间自然搜索的方式编码染色体,即将节点的二维坐标,连续时间区间的爆管流量作为染色体基因位,具体如下表所示:
Cordx Cordy Demand1 Demand2 Demand3 … Demandn
每个节点对应一个二维坐标,反之,每个二维坐标可映射到一个距离其最近的节点;在搜索时,改变二维坐标大小,即可改变目标节点;Demand1代表优化校核第一个时刻的节点爆管流量大小,Demand2则代表下一时刻,后面同理;
确定基于压力和流量拟合误差最小化的目标函数,利用获取的管网边界数据建立实时模型,迭代进化模拟优化实时模型,求得精确的爆管点和爆管漏损量。
说明书

技术领域

[0001] 本发明属于城市供水管网领域,涉及供水管网的爆管异常定位,具体是一种基于SCADA数据和实时模型的城市供水管网爆管在线定位方法。

背景技术

[0002] 爆管漏水在供水管网中普遍存在,具有突发性,会导致供水系统运行异常,影响居民生活和企业生产,甚至会危害人身安全。然而,爆管侦测与定位较为复杂,涉及爆管发生时管网系统动态状况,受到测点少、测量精度低、采集传输实时性差、水力模型不可靠等因素影响。目前管网爆管定位方法主要有基于在线数据分析(如瞬变流逆分析等)方法和基于微观模型(如离线模型异常校核)方法,前一类方法基于在线数据,可一定程度上实现爆管流量初步估计和爆管点区域划定,但易受管网噪声影响,易漏诊,且定位误差大;后一类方法虽可在理论上实现精确定位,但离线模型往往精度较低,存在固有误差源,同时,由于爆管会对管网系统产生破环性影响,直接影响管网各个区域压力态势,同时会对入水口水压造成波动,传统离线模型忽略这部分变化,易产生较大定位误差。

发明内容

[0003] 有鉴于此,本发明针对供水管网在线定位的需要,提出了一种基于SCADA数据和实时模型的城市供水管网爆管在线定位方法。
[0004] 为实现上述目的,本发明采取以下步骤:
[0005] 步骤1 、监测压力波谷,确定疑似爆管时间边界和校核测点
[0006] 根据各测压点在爆管发生前后的压力变化的先后顺序和压力波谷,得出测点压力因爆管发生变化的起始时间、截止时间,从而获得用于实时模型优化校核的起始时间 ,持续时间T以及有压力值反馈的时刻k,选择负压波变化比较明显的测点作为校核测点,设其个数为m。
[0007] 步骤2、初步定位
[0008] 爆管发生前后,测点与水源压差的比值(压降比)反映了爆管发生后测点附近的流量异常情况,用公式表示为:
[0009]
[0010] 表示管网中各测点 在爆管发生后 时刻时附近区域的流量异常率, 表示爆管发生起始时刻, 表示爆管发生后 时刻测点 与水源水头的压差,表示 时刻测点 与水源水头的压差。
[0011] 根据每次爆管发生后测点的流量异常率变化,选择测点当中异常率变化明显的测点,按照它们所在位置,划定矩形区域,作为实时优化校核模型搜索爆管点的初始区域。
[0012] 步骤3、精确定位
[0013] 目标函数以连续时间内压力测点的压力模拟值和实际值、流量测点的模拟值与实际值之间的差值绝对值之和作为评价标准,迭代求解最小值,获得最可能的爆管节点和爆管流量。具体为:
[0014]
[0015] 其中, 为压力测点 在 时刻的优化模拟值, 为压力测点SCADA监测值;为测流管线 的优化模拟值, 为测流管线 的SCADA监测值,T为优化校核时间段,为由步骤1选定的压力测点数目,n为供水管网中管线流量测点数目, 、 是分配的权重系数,具体由管网实际情况经验值给定。
[0016] 通过模型模拟某节点的流量大幅度变化来模拟爆管现象,以节点位置以及节点漏失量作为优化过程的自变量。通过遗传算法优化求解。
[0017] 设定遗传算法参数,包括种群规模、进化代数、选择、交叉、变异方法、交叉概率、变异概率;确定染色体编码方式,以空间自然搜索的方式编码染色体,即将节点的二维坐标,连续时间区间的爆管流量作为染色体基因位,具体如下表1所示:
[0018] 表1
[0019]…
[0020] 每个节点对应一个二维坐标,反之,每个二维坐标可映射到一个距离其最近的节点。在搜索时,改变二维坐标大小,即可改变目标节点。 代表优化校核第一个时刻的节点爆管流量大小, 则代表下一时刻,后面同理。
[0021] 确定基于压力和流量拟合误差最小化的目标函数,利用获取的管网边界数据建立实时模型,迭代进化模拟优化实时模型,求得精确的爆管点和爆管漏损量。
[0022] 本发明的有益效果:本发明能实现对爆管等异常事件的高精度快速侦测定位,降低爆管对供水管网系统造成的危害。

实施方案

[0027] 为使本发明的技术易于理解,下面结合附图给出一个实例,对本发明的具体实施方式作进一步详细描述。
[0028] S市供水公司于2015年4月3日上午9点到12点,做了 5次模拟爆管实验。本实例以该S市做爆管实验的小型DMA管网为研究对象,如图1所示,该管网包含5293个普通用水节点,5415条管道,14个内部压力监测点,18个大表流量测点,2条内部测流管线,图中标注的41922、26124和26125为管网中的3个入水口,压力和流量均有测量。已知左侧矩形方框区域曾发生过一次爆管事件,左侧粗原点标记的节点为内部压力测点,中间圆圈圈起来的两条管线是内部测流管线。
[0029] 步骤1、监测压力波谷,确定疑似爆管时间边界和校核测点
[0030] 根据SCADA系统获取爆管发生的大致时间段和压力数据,绘制爆管发生时间段内14个测点的压力曲线图,如图2所示,确定14个测点压力普遍开始下降的起始时间、降幅较大的主要节点、下降到压力波波谷的平均时间。
[0031] 分析图2可知,14个测点压力普遍开始下降的时间是第4分钟开始,第8分钟普遍下降到压力波波谷。其中,测点25727、23770、1193、32048下降幅度较大,下降速度也较快,其余测点之间压力变化大致相同,且降幅较小。
[0032] 根据14个测压点的时间边界,确定实时模型起始时间 =4,连续校核时间区间T=4,选择25727、23770、1193、32048作为校核目标测点,测点数目m=4,忽略其他变化不明显的压力测点。爆管发生后,测点采集到压力值的时刻k取值为:1,2,3,4,……,10。
[0033] 步骤2 初步定位
[0034] 依据以下公式,求测点的流量异常率
[0035]
[0036] 表示管网中各测点 在爆管发生后 时刻时附近区域的流量异常率, 表示爆管发生起始时刻, 表示 时刻测点 与水源水头的压差, 表示爆管发生后 时刻测点 与水源水头的压差。
[0037] 求得爆管发生后10分钟内流量异常率如下表2所示
[0038] 表2
[0039]测点/异常率/时间 1 2 3 4 5 6 7 8 9 10
605 1.00 1.13 1.27 1.47 1.40 1.49 1.47 1.49 1.38 1.34
850 1.00 1.04 1.04 1.04 1.00 1.04 1.00 0.99 0.99 0.99
1193 1.00 1.39 1.34 1.40 1.34 1.40 1.40 1.35 1.37 1.30
1244 1.00 1.06 1.03 1.04 1.03 1.04 1.00 1.03 1.01 1.01
3120 1.00 1.02 1.32 1.45 1.40 1.41 1.41 1.56 1.38 1.32
3487 1.00 1.05 1.12 1.37 1.39 1.45 1.45 1.47 1.47 1.29
5604 1.00 1.02 1.29 1.34 1.32 1.29 1.44 1.44 1.24 1.20
23770 1.00 1.05 1.02 1.03 1.00 1.02 0.94 0.99 0.97 0.97
25555 1.00 1.06 1.01 1.02 1.00 1.01 0.94 0.97 0.95 0.98
25727 1.00 1.03 1.02 1.00 1.02 1.00 0.95 0.95 0.95 0.91
27115 1.00 1.01 1.00 1.01 0.97 0.92 0.93 0.97 0.99 0.97
28211 1.00 1.05 1.02 1.01 1.00 1.02 0.96 0.98 0.99 0.98
32048 1.00 1.31 1.27 1.30 1.24 1.33 1.33 1.28 1.25 1.22
40084 1.00 1.06 0.99 1.04 0.99 0.96 0.97 1.00 1.00 0.99
[0040] 绘制14个测点异常率变化曲线图,如图3所示。找出异常率变化较大的几个节点,作为初步搜索区域。
[0041] 统计得出异常率值较大,变化较明显的几个测点分别是:1193、32048、605、3487和3120。异常率值得分布区间是(1.2,1.6)。
[0042] 根据选取节点,在EPANET供水管网模型软件中获得其具体坐标。将异常率大于1的这几个主要测点,作为优化校核模型的初步搜索边界,也就是自变量中坐标参数的取值边界;具体是将异常率大于1的测点的节点坐标进行排序,选择X坐标中最小,最大值,Y坐标中最小、最大值作为优化校核模型的左右、下上四个搜索边界。爆管流量的初步取值范围依据历史压力异常率对应的爆管流量幅度区间经验确定。如下表3所示
[0043] 表3
[0044]爆管序号/值/参数 Xmin Xmax Ymin Ymax 爆管流量范围
1 57,811.39 59,674.42 18,223.87 20,246.66 0 1,000.00
~
[0045] 步骤3.精确定位
[0046] 依据步骤1和步骤2确定优化校核测点、优化校核时间区间、爆管搜索范围和爆管漏损量取值范围进行精确定位。
[0047] 采集第4分到第8分的SCADA数据,具体是出、入口压力数据、流量数据、优化目标测点的压力数据、管道状态数据、大表用户流量数据等。其中,入口压力数据在线更新为水库水头,大表流量数据直接更新为对应节点的需水量数据,依据管道的开度,直接更新对应管道为开或者关。
[0048] 以爆管发生时间段的测点实时模拟值与实测值之差,作为优化校核目标函数,具体为:
[0049]
[0050] , 依据压力拟合差和流量拟合差的比值区间经验给定,该管网是1/10 的关系,则 设为1, 为0.1,以平衡两者在误差累计中的权重,管网中有两条流量测点,故n=2。
[0051] 编码、选择、交叉和变异方法分别采用实数编码、轮盘赌选择、实数交叉和单点变异;进化代数 ;种群规模 ;交叉和变异概率分别为 , 。染色体的具体形式如表4所示
[0052] 表4
[0053]X坐标 Y坐标 t时刻需水量 t+1时刻需水量 t+2时刻需水量 …
57,822.39 59,624.42 345.21 452.12 321.78 …
[0054] 利用EPANET水力求解引擎将选取的SCADA在线数据按优化时间对应更新到管网模型上去,实现实时模拟。
[0055] 通过遗传算法迭代进化,改变染色体基因,评价适应度函数,实现实时校核,最终获得最好的染色体,从而得到最可能的爆管坐标和爆管流量,方案流程图如图4所示。利用EPANET获取所有DMA管网的节点坐标,求得与爆管坐标欧氏距离最近的节点坐标,确定该节点即为爆管节点。
[0056] 根据以上方法,求得5次爆管定位偏差见表5所示,求解得到的爆管漏损量见表6。
[0057] 表5求解爆管点与实际爆管点距离偏差
[0058]实验序号/值/参数 爆管求解点 实际点 欧式距离差值
1 1671 26551 249.10966
2 618 29675 73.932308
3 75 25575 631.7071
4 17137 31733 478.93182
5 54199 25133 997.32585
[0059] 表6最优解对应的爆管流量
[0060]实验序号/值/参数 最优解 第1分钟 第2分钟 第3分钟
1 1671 166.1869903 58.48658 10.98238
2 618 172.1317902 185.5464 160.7959
3 75 11.65619662 12.06562 11.26899
4 17137 10.86489503 13.45882 12.09842
5 54199 10.83961735 125.496 10.37903

附图说明

[0023] 图1为S市小型DMA管网;
[0024] 图2为14个测点的压力曲线图;
[0025] 图3为14个测点异常率变化曲线图;
[0026] 图4为爆管定位具体方案流程图。
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号