CN106575312B - 用于检测配送网络特别是水配送网络中的异常的方法 - Google Patents
用于检测配送网络特别是水配送网络中的异常的方法 Download PDFInfo
- Publication number
- CN106575312B CN106575312B CN201580040183.2A CN201580040183A CN106575312B CN 106575312 B CN106575312 B CN 106575312B CN 201580040183 A CN201580040183 A CN 201580040183A CN 106575312 B CN106575312 B CN 106575312B
- Authority
- CN
- China
- Prior art keywords
- parameters
- model
- parameter
- value
- primary
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title description 14
- 230000002441 reversible effect Effects 0.000 claims abstract description 16
- 238000005259 measurement Methods 0.000 claims description 65
- 230000035945 sensitivity Effects 0.000 claims description 26
- 230000002159 abnormal effect Effects 0.000 claims description 14
- 238000004088 simulation Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 8
- 239000012530 fluid Substances 0.000 claims description 4
- 238000012876 topography Methods 0.000 claims description 4
- 230000000717 retained effect Effects 0.000 claims 1
- 230000006870 function Effects 0.000 description 12
- 230000000875 corresponding effect Effects 0.000 description 6
- 235000020188 drinking water Nutrition 0.000 description 6
- 239000003651 drinking water Substances 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000001514 detection method Methods 0.000 description 5
- 230000002547 anomalous effect Effects 0.000 description 4
- 238000009530 blood pressure measurement Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 3
- 230000001747 exhibiting effect Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000004913 activation Effects 0.000 description 1
- 230000004931 aggregating effect Effects 0.000 description 1
- 238000004378 air conditioning Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000005293 physical law Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 235000020681 well water Nutrition 0.000 description 1
- 239000002349 well water Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F15/00—Details of, or accessories for, apparatus of groups G01F1/00 - G01F13/00 insofar as such details or appliances are not adapted to particular types of such apparatus
- G01F15/06—Indicating or recording devices
- G01F15/068—Indicating or recording devices with electrical means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0243—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L19/00—Details of, or accessories for, apparatus for measuring steady or quasi-steady pressure of a fluent medium insofar as such details or accessories are not special to particular types of pressure gauges
- G01L19/08—Means for indicating or recording, e.g. for remote indication
- G01L19/083—Means for indicating or recording, e.g. for remote indication electrical
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2117/00—Details relating to the type or aim of the circuit design
- G06F2117/02—Fault tolerance, e.g. for transient fault suppression
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Automation & Control Theory (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Computer Hardware Design (AREA)
- Fluid Mechanics (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Public Health (AREA)
- Testing And Monitoring For Control Systems (AREA)
- Fluid-Pressure Circuits (AREA)
- Control Of Fluid Pressure (AREA)
- Flow Control (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
使用一种网络的水力模型,所述水力模型描述结构数据以及控制配送的流的定律。用与初级输入参数相关的操作数据供给所述模型,并借助模型获得针对初级输出参数的理论值。比较理论值与测量值。在显著偏差的情况下,相对应的初级输出参数在包括新的次级输出参数的反向模型中变为次级输入参数,所述新的次级输出参数相对于正向模型被添加或对应于初级输入参数。优选地,次级输出参数是具有异常值的初级输出参数对其特别灵敏的参数。需要在逐渐地限制次级输出参数的数量的同时,以迭代的方式实施该方法。
Description
技术领域
本发明涉及一种检测流体供应系统,特别是牛顿流体的配送网络,更特别地用于饮用水配送,还例如供热和/或空气调节城市网络中的异常的方法。
背景技术
操作性能处于这种网络的管理的核心地位。借助例如在饮用水运输及配送网络上的水力(hydraulic)异常的检测和定位工具,可以显著改善性能等级。
在这种网络中,可能发生不同类型的事故并且事故可能具有对运营者而言严重的后果。这可以涉及泄漏,对消费者供应的中断,组件的运行故障,在介入期间或者在介入之后的操纵错误,在消费者家中的运行故障、异常表现或恶意动作。
在扩展程度很高并且部分不可达(被掩埋)的网络上,直接的监视是低效率的或者经常是不可能的。通过测量仪器的监视也是不足够的,因为本地异常可以对网络的广泛区域具有影响,或甚至对整个网络具有影响,这使诊断复杂化。
已知基于统计技术的检测方法。这些方法既不允许表示异常的特征也不允许定位异常。
其他已知的检测方法基于使用用于检测和/或定位异常的演化算法(遗传算法)的水力模型。这些演化方法非常消耗计算资源,这使得难以大规模应用。
发明内容
本发明的目的是提出一种用于检测和/或定位和/或表征和/或量化配送网络中的异常的有效方法,该配送网络在必要时可以是巨大的,如饮用水配送网络。
根据本发明,检测配送网络中的异常的方法,包括以下步骤:
-为网络建立包括以下元素的正向模型:
○网络的结构数据,
○操作数据,所述操作数据与被选择作为描述操作场景的初级输入参数的参数相关,
○定律,所述定律将所述结构数据、初级输入参数以及初级输出参数相连接(对模型进行仿真);
-执行正向模型,用于确定初级输出参数的理论值(求解反向问题);
-获得初级输出参数的测量结果;
-比较理论值与测量结果;
-在比较揭示出由于理论值与测量结果之间的大偏差而导致的至少一个异常初级输出参数的情况下,建立反向模型,该反向模型是通过将至少一个异常初级输出参数引入作为反向模型的次级输入参数,从正向模型导出的,除了被保持的初级输出参数之外反向模型还包括至少一个次级输出参数;
-执行反向模型。
优选地,配送网络是流体配送网络。
反向模型的建立可以包括相对于所有可配置的参数使在测量点处的测量结果与仿真之间的偏差函数最小化。
次级输入参数可以将测量值采纳为值。
优选地,由水力仿真引擎实现反向模型的执行。
求解反向问题可以实施使理论输出与测量输出之间的偏差最小化的数字算法,同时使用在正向模型中的相对于输入的输出的灵敏度的计算。
结构数据典型地与网络的地形(topography)、网络组件的性质以及特征等相关。
将“初级输入参数”称为运行参数。典型地,运行参数是时间中的变量。例如,在水配送网络中的流率、压力等,其中,使用该值作为用于正向模型的输入数据。该值可以由测量仪器提供,或以另一种方式已知,例如取决于日期和时刻的消耗统计。通常初级输入参数包括描述在网络的边界(limit)条件下的参数。例如,在水力网络中,网络的供应流率和提供给消费者的流率。
以通常方式,形容词“初级”表示“与正向模型相关”,而形容词“次级”表示“与反向模型相关”。
“场景”是网络的运行情况。例如,在水力网络中,场景以供应储罐中的特定充满水位开始,并且在场景进展期间,在根据在消耗点处取得的流率向网络提供水的同时储罐以特定流率被再次供给。不同的流率根据构成场景的部分的相应的时间图而演变。
特别地,定律(law)是可应用的物理定律。例如,在水力网络中,根据流动速度以及两个点之间的网络的结构数据,网络的所述两个点之间的水头损失(head loss)(压力下降)。
“初级输出参数”是通过执行用于获得其理论值的正向模型的可计算的参数,并且同时初级输出参数是借助提供初级输出参数中的每一个的测量结果的网络仪器可测量的或例如根据之前的统计已知的。
在根据本发明的方法中,比较初级输出参数的理论值与实际值,并且在理论值与测量值之间不一致的情况下,认为在网络中存在异常。
当初级输出参数为“异常”时,这意味着在其理论值与测量值之间存在超过了可接受的错误界限,进入反向模型的求解(执行),对于反向模型,异常初级输出参数变为次级输入参数,对于该次级输入参数将测量的值给定为值。尤其为了使反向模型具有必要的数学自由度,在反向模型中引入由新的次级输出参数构成的至少一个额外值。优选地,根据得自于在执行正向模型后与测量结果的比较的异常,合理地选择这些新的次级输出参数。
反向模型的执行提供直接消息性的结果,例如新的输出参数具有允许用足够的确定性识别异常的理论值,或者如下面将阐述的,允许通过迭代向这种标识转换。
因此,本发明巧妙地利用由“定律”连接的所有变量(运行参数)的一致性,所述“定律”典型地是网络特别是水力网络的表现的方程。
至少一个次级输出参数可以是相对于正向模型被添加到水力模型的参数。
例如,如果怀疑在网络点处的泄漏,可以添加泄漏流率作为借助反向模型寻求确定的次级输出参数。
至少一个次级输出参数还可以对应于相对于正向模型的至少一个删除的数据。例如,如果初级输出参数指示阀门处于打开位置,而通过针对至少一个初级输入参数的理论值与测量值之间的偏差揭示的异常暗示出该阀门可能关闭或部分关闭,则所涉及的初级输入参数被选择作为次级输出参数,并且从用于反向模型执行的输入参数中删除。
在优选的实施方式中,通过应用揭示以下可能性的至少一个准则来选择至少一个次级输出参数:至少一个新的次级输出参数是在测量结果与理论值之间所观察到的偏差中涉及到的可能性。
再次参照之前的示例,在网络的给定点处相对于理论值异常的所测量的流率值,例如很可能表示被认为打开的阀门是关闭的,或者存在严重泄漏,或者在下游的异常逆向流动流率(非限制性地列出)。
根据本发明的非常优选的特性,操作数据、理论值和测量结果中的至少一些包括时间序列,即,值序列中的每一个与时间表相关联。
因此,不仅以统计的方式而且以动态的方式操作正向模型和反向模型。这允许模型进行额外的一致性检查。例如,模型可以具有使储罐水位变化与在网络的不同点测量的流率和相关的定律。此外,在异常情况下,不仅获得在运行时刻的偏差值,还获得偏差函数,所述偏差函数给出取决于全部或部分初级(正向模型的执行)或次级(次级模型的执行)参数(所述参数是输入或输出参数)的偏差。
测量结果,无论被用于获得一些运行参数,还是关系到用于获得要与理论值比较的测量值的输出参数,随后同样地以时间序列的形式被获得。在与理论值比较之前,优选地使测量结果经受平滑而不是在原始状态使用测量结果。因此,偏差函数本身被自动地平滑。
在实施场景初期的初始化步骤中,难以获取一些运行数据。典型地,这是消耗情况。例如,涉及水配送网络,难以逐个小时地知道由网络或由该网络的每个子网络供给的消耗。根据本发明,预先通过历史已知的数据或通过经验代替这种难以实时获取的数据。例如,通过之前的学习,知道在白天和夜晚的每个小时给定子网的估计消耗,连同关于这些值的不确定性信息。至少在初始化阶段,可以使用这些时间序列作为运行数据。
一些结构或运行数据可以是未知的或不确定的。例如,在水配送网络中,压力测量点的高度是重要的,因为高度每额外增加一米对应于压力约减小10kPa,需要将其考虑在内从而使利用控制(govern)流动动态的方程获得的结果合理。管道直径、管道的粗糙度以及影响在流动中的水头损失的其他参数可能是不充分已知的,尤其是在旧设施中。
尤其针对这种情况,可能是明智的,在初始化步骤期间,使用反向模型用于确定不可能或难以直接获得的一些结构或操作数据,根据针对一些初级输出参数的测量值或已知值,在认为已获得这些初级输出参数的值的情况下执行反向模型。
根据本发明的重要方面,在正向模型的执行中整合至少一些初级输出参数对至少一些其他参数的变化的灵敏度的计算,并且对于反向模型,选择呈现异常的初级输出参数对其更加特别灵敏的初级输入参数作为次级输出参数。
灵敏度可以被作为参数相对于其他参数中的每一个的导数(已知将这些参数彼此相连的定律)而理论地获得,或者还通过针对序列的全部值比较输出参数测量结果的变化与其他参数的变化而实验地获得。在估计了灵敏度的参数的值的不同范围中,并且在以下参数的值的不同范围中灵敏度可以是不同的:相对于该参数估计了灵敏度。因此,灵敏度可以被解释为值的表或函数。根据输出参数按以下参数的方向或相应地相反方向演变,灵敏度可以是正的或负的:相对于该参数计算了输出参数的灵敏度。
在优选的实施方式中:
-通过正向模型,根据次级输出参数,获得针对异常初级输出参数的理论值与测量结果之间的偏差值的偏差函数;
-根据灵敏度计算偏差函数相对于参数的梯度;
-通过反向模型的应用来确定新数据的集合,修改初始数据,并且对于新数据而言偏差函数的总体梯度具有零值,指示相对于参数互相抵消的正灵敏度和负灵敏度;
-从这些参数中选择子集合,所述子集合由灵敏度的符号相同的参数形成;
-在将次级输出参数缩小到构成子集合的部分的次级输出参数的同时,以迭代的方式重新开始所述方法;
-并且如此重复直到利用单一参数或小数量的参数定位异常。
为了构成子集合,优选地根据所涉及的最可能的参数的性质、得自于按将异常理论值向相对应的测量值引导的方向的变化的实际值,选择对应于所涉及的最可能的参数的灵敏度的符号。
典型地,当单一初级输出参数呈现出相对于理论值的测量值的异常时,优先考虑测量错误的假设,特别是在呈现出异常的参数仅通过其他测量点与网络的地形边界相连的情况下。在这种情况下,可以在反向模型的应用之前,以检验测量设备开始。然而,可以按迭代的方式应用反向模型从而通过在将异常测量采纳为假设的情况下确认对于其他参数而言没有任何值的组合是可信的,来得出测量错误的结论。如果测量来自于网络边界的测量点,异常有更大的可能来自于超出网络的地形边界的反常元素。例如,消费者将井水注入到网络中,而不是取用由网络供应的水。可能应当重新考虑在消耗点上的操作数据。
在完善的实施方式中,周期性地执行正向模型以及理论值与测量值的比较,并且当至少两个初级输出参数呈现相对于理论值的测量值的异常时,自动激活反向模型的执行。
优选地,模型的元素还包括关于数据及参数的不确定性的指示,并且在比较期间考虑这些不确定性用于将输出参数的测量结果认为正常或相反地是异常的。
具体实施方式
根据下面与关于饮用水配送网络(还称作“水力网络”)的非限制性例子相关的描述,本发明的其他特性及优点将变得更清楚。
初步说明
接下来的描述是关于其包含的任何特性的描述,可以独立于其他特性考虑所述特性,即使这些特性形成相同段落或相同句子的部分,并且接下来的描述是关于这些特性的任何组合的描述,因此,无论与上面介绍的概念结合与否,这种特性或特性的组合是区别于现有技术的并提供了技术效果,并且该特性可以用本说明书中的相同术语或相对而言更一般的术语表达。
词汇表
使用以下词汇表:
时间序列:由时间索引的标量数据的有限序列,该序列通常被恒定的持续时间分隔开。
水力模型:在没有其他明确表示的情况下,表示饮用水运输和/或配送网络的水力模型:数据集合,所述数据集合描述拓扑(网络的图形)、地形(高度)、定义区段的水力表现(长度、内径、材料和/或粗糙度、局部水头损失系数(singular head loss coefficient)、奇点(singularity)等)的区段特征、定义设备的水力表现(直径、水头损失系数、设定值、水泵扬程等)的设备特征、与图形节点相关联的正常消耗、以及描述运行条件的场景。
场景:形势数据或操作数据的集合,其描述与水力模型一起使用从而仿真在特定条件中的水力系统的表现的边界条件,诸如储罐的初始水位,区和/或消费者类型的消耗的调制(modulation)的时间序列,泵、阀门的状态,被控制的活跃设备的设定值,以及/或者伺服规则。
测量结果:来自数据获取系统的时间序列。测量结果与水力模型的特定基本组件(对于流率测量的弧或区段,对于压力测量的节点,对于水位测量的储罐)相关联。
正向模型:水力模型的仿真包括,通过使用在场景中描述的边界条件来数字地求解由水力模型以及其动态演变描述并确定参数的稳定流动方程。结果是在系统任意点处以及任意时刻的流率、速度、水头损失、状态、压力,这些结果通常以时间序列的形式或以图(map)的形式呈现。
反向模型:(求解反向问题)包括数字地求解稳定流动方程以及其在时间上的演变。通过使用一些流率(或速度)和/或压力和/或储罐的水位的测量结果,通过水力模型和场景来相对于模型或场景的一个或多个参数对稳定流动方程以及其在时间上的演变进行描述并确定参数。
远程读取的数据:对于流量计的消耗读数序列,所述数据在1至24小时的周期性被测量,并且一般地每天至少一次被远程传送。
特征:包含对所研究的过程有意义的信息的标量值或时间序列。
综述
根据本发明的方法的实施方式联合地使用来自远程监视系统的操作数据,使得能够构造场景、水力网络系统的水力模型、以及分析模型执行结果的系统。
水力模型化程序整合相对于输入数据(在节点处的消耗和在弧处的水力阻力)的结果(测量)的灵敏度的显式计算,允许计算相对于参数的仿真错误的梯度,然后解决反向问题。
反向模型的使用还允许缓解构造场景所必须的一些测量结果的缺失。
正向模型的使用允许检测系统的表现的异常。
反向模型的使用允许在空间中并且可能地在时间上定位异常的可能起源。
可选的远程读取的数据的使用允许更新在水力模型中描述的消耗的空间分布。
详细描述
被提供作为用于模型的执行、与初级输出参数的测量结果的比较、并且在必要情况下反向模型的执行的输入的元素包括:
-网络的结构数据,其可以包括设备的状态数据和设定值;
-被选择用于构成初级输入参数以及其标量值的时间序列的参数,通过测量结果或者通过对全部或部分订户可用的远程读取的数据,或甚至通过存在的统计获得所述参数;
-被选择用于构成初级输出参数以及这些参数的测量结果的时间序列的参数;
-关于以下内容的信息:测量结果及参数的准确性和不确定性、其相关的权重、被选择作为次级输入参数以及在解决反向问题时的次级输出参数的参数的激活。
基本操作
-根据测量结果和状态序列构造场景
该操作包括将测量结果改编成由水力模型软件可用的形式。例如,水力区的输入/输出流率的时间序列的代数和被转换为该区的额定消耗的调制序列。
-求解正向问题、计算特征
该操作包括利用将系统的水力模型和由之前操作创建的场景作为输入数据,启动水力仿真引擎。所利用的输出是与针对这些参数可用的测量结果可比拟的仿真的时间序列(初级输出参数的理论值)。
-计算特征
针对每个测量点(初级输出参数),基于在测量结果获取步骤计算的残余值序列来计算特征,所述残余值序列被定义为所仿真的值(理论值)与所测量的值(经平滑或未经平滑)之间的差。所述特征指示是否存在异常。
-解决反向问题(一般情况)
解决反向问题是迭代操作,其包括借助优化算法(通过平方和方法或剩余序列的绝对值方法来进行算法运行),相对于可配置的参数集使在测量点处的测量结果与仿真(理论值)之间的偏差函数最小化。在偏差函数中,通过根据先验不确定性以及测量结果的局部经验方差计算的系数对每个序列的每个点进行加权。
在预定的时间窗口(例如,从4小时到24小时)中通过使用加有日期戳的测量结果来实现计算。
最小化问题的每次迭代需要借助计算正向问题来计算偏差函数。所使用的优化方法要求计算偏差函数相对于参数的梯度。根据测量结果相对于水力模型的参数的灵敏度(导数)显式地计算该梯度的值,这是以与求解正向问题相嵌套的方式显式执行的。
本发明的方面之一是,一方面灵敏度的显式计算是与求解正向问题以及使用灵敏度来计算梯度同时进行的;另一方面,是选择与求解出反向问题有关的参数的可能性,例如:
●结构数据(区段以及设备的水力特征,诸如粗糙度或局部水头损失系数(PCS)之类)
●压力和/或水位传感器的刻度(高度)
●消耗的空间分布
●预定的空间分布的一个或更多个消耗类型的时间调制。
本创新的另一个重要元素在于,当识别的参数与模型的多个组件(component)相关联时,使用灵敏度用于定位异常。在这种情况下,当获得了解答时,因为总体梯度在最佳状态时值为零,所以灵敏度为一些组件采取正值,为其他组件采取负值。因此,只需要通过相对于梯度符号相同的分量的两个子集合之一求解来重复所述方法,以便确定异常的位置。由组件的嵌套的子集合如此形成的顺序以逐步准确的方式定位异常。
在大测量密度并且测量的位置对于确定所选参数类型是合理的情况下,定位的准确性更高。
-应用于在初始化步骤期间重构必要测量结果的反向问题
当构成场景(例如区的输入或输出流率)所必要的测量结果不可用时,求解反向模型允许重构该缺失的测量结果,因此聚集了可观察性的一些条件,在这种情况下,存在其值与缺失的测量结果相关的其他测量结果(例如储罐水位的测量结果和/或压力的测量结果)。
因此,通过使用预定的序列来构成初始场景,并且相对于由缺失序列构成的参数求解反向问题。
-应用于检测异常的反向模型的执行
应用于识别诸如粗糙度或PCS(局部水头损失)之类的结构数据,反向问题的求解(反向模型的执行)允许调整水力模型、识别诸如管道的长度或直径的输入错误之类的水力阻力异常、定位在介入之后忘记的关闭的阀门。
应用于识别诸如传感器的高度或设备的设定值之类的结构参数,求解反向问题允许诊断在导致了网络参数之一的修改的介入之后水力模型的代表性的丢失。
应用于识别消耗的空间分布或者一个或更多个消费者类型的调制系数序列,求解反向问题允许定位消耗的异常。每个参数的理论值与测量值之间的偏差构成可用于异常检测的特征或定量标识符。
-特征与异常的解释
一致性标识符:对于每个特征,计算一个或更多个定量标识符:Nash准则、“协定索引”(“Index of Agreement”)、超过错误阈值的率。
异常:对于每个测量结果,相对于操作地设定的阈值,根据测量结果的一致性标识符的位置计算定量标识符(优/良/中/差)。被分类为“中”或“差”的值被识别为异常。
在区上隔离的异常被分类为“测量异常”;在区上观察到多个异常触发“网络异常”。
优点
由上面阐述的组件构成的集合可以与饮用水供应系统的运营者的技术信息系统相连接。以如下方式构造集合:其可以按包括在测量数据的获取频率之间的频率并且一天一次被规律地激活。在这些条件下,异常的检测和特征表示比通常使用的方法快得多,并且显著改善了操作效率。
在针对网络部署水力模型的阶段中,本发明允许极大程度地减少模型调整时间。
第一实施例
本发明被配置用于模型化由两个高位储罐A和B通过重力供给的水力区。
水力模型包括PICCOLO/GANESSA_SIM格式的文件。
所使用的场景生成器为GANESSA_GS。
并入求解正向问题(模型执行)或反向问题的水力仿真引擎是GANESSA_SIM。
以每天一次的频率供给模型的元素是每个储罐的输入流率、每个储罐的水位、每个储罐的输出流率以及在远点处的压力。储罐A的输出流量计具有包括约10%的未被记入的异常。
针对24小时的正向水力模型化产生储罐的输入和输出流率的代数和的合计等。“计算的水位减去测量的水位”特征给出关于流率测量精确性、指示未被记入的输出流率的在时段末的正偏差、或者输入流率的超额记入的信息。根据该特征计算的标识符使得异常凸显出来;针对输出流率测量结果中的每一个的特征的符号指示异常的方向。
利用作为次级输出参数的在每个储罐的输出点固定的额外消耗和系数的时间序列配置的反向模型,针对储罐的每个输出点计算两个强度序列的调整。与储罐A的输出点相关联的序列具有与A的输出流率成比例的轮廓(profile),对应于10%数量级的消耗,而与B相关联的序列具有可忽略的强度。
第二示例
本发明被配置用于模型化由在A和B两个点的主管道供给的水力区。
水力模型包括EPANET格式的文件。
所使用的场景生成器为GANESSA_GS。
并入求解反向问题的水力仿真引擎是GANESSA_SIM。
测量结果是输入(A点和B点)的压力和流率,并且在配送网络上分布的六个附加点处获得压力测量结果。
首先,本发明被离线使用以便调整模型的水力设置(calage):反向模型允许识别压力测量点的准确高度、两个输入阀门的水力阻力、以及构成区的管道的4个主要类型材料的平均粗糙度。
其次,本发明按天的频率被实施。小的附加消耗(例如1l/s)在模型的所有节点上例如按照与节点相连的管道的长度比例分布,并且被分配特定消费者代码(例如“LEAK”),并且将用零值或单一值初始化的时间序列与该“LEAK”代码相关联。随后,相对于“LEAK”的系数的时间序列的系数求解反向问题。随后分析相对应的特征:如果相对应的消耗低于预定阈值,则认为该情况是正常的。在相反情况下,推测出消耗异常。随后,在为对其来说在之前步骤计算的灵敏度为负的个别节点重新分配附加消耗后,重新实施本发明,并且多次重复,直到剩余节点的数量低于预定阈值。随后,相对应节点的最后的集合被报告作为异常定位,并且附加消耗(与“LEAK”相关联)的平均值作为其强度。
当然,本发明不限于所描述的示例。
Claims (13)
1.一种检测流体配送网络中的异常的方法,包括以下步骤:
-为配送网络建立包括以下元素的正向模型:
ο配送网络的结构数据,所述结构数据与所述配送网络的地形相关并包括所述配送网络的设备的状态数据或设定值,
ο操作数据,所述操作数据与被选择作为描述操作场景的初级输入参数的参数相关,
ο定律,所述定律将所述结构数据、初级输入参数以及初级输出参数相连接;
-由水力仿真引擎执行正向模型,用于确定初级输出参数的理论值;
-在以上步骤的全部或部分之前或之后,获得来自初级输出参数的数据获取系统的测量结果;
-比较理论值与测量结果;
-在比较显示出由于理论值与测量结果之间的大偏差而导致的至少一个异常初级输出参数的情况下,相对于一组能配置的参数建立包括使在测量点处的测量结果与仿真之间的偏差函数最小化的反向模型,反向模型是通过将至少一个异常初级输出参数引入作为反向模型的次级输入参数而从正向模型导出的,次级输入参数将测量值采纳为值,并且除了那些被保持的初级输出参数之外反向模型还包括至少一个次级输出参数;
-由水力仿真引擎执行反向模型;
-根据反向模型的执行定位配送网络的异常。
2.根据权利要求1所述的方法,其特征在于,至少一个次级输出参数是相对于正向模型被添加到模型的参数。
3.根据权利要求1或2所述的方法,其特征在于,相对于正向模型,至少一个次级输出参数对应于至少一个删除的数据。
4.根据权利要求1或2所述的方法,其特征在于,通过应用指示以下可能性的准则来选择至少一个次级输出参数:至少一个新的次级输出参数是在测量结果与理论值之间所观察到的偏差中涉及到的可能性。
5.根据权利要求1或2所述的方法,其特征在于,操作数据、理论值和测量结果中的至少一些包括时间序列,即,每个时间序列是一个值序列,值序列中的每个值与时间表相关联。
6.根据权利要求5所述的方法,其特征在于,在能配置的特征性持续时间之后,由根据测量结果获得的值形成的时间序列被平滑。
7.根据权利要求1或2所述的方法,其特征在于,在初始化步骤中,使用反向模型用于确定不可能或难以直接获得的一些结构数据或操作数据,在认为已获得一些初级输出参数的值的情况下执行反向模型。
8.根据权利要求1或2所述的方法,其特征在于,在正向模型的执行中整合至少一些初级输出参数对至少一些其他参数的变化的灵敏度的计算,并且对于反向模型,选择异常初级输出参数对其更加特别灵敏的初级输入参数作为次级输出参数。
9.根据权利要求8所述的方法,其特征在于:
-通过正向模型,根据次级输出参数,获得针对异常初级输出参数的理论值与测量结果之间的偏差值的偏差函数;
-根据灵敏度计算偏差函数相对于参数的梯度;
-通过反向模型的应用来确定新数据的集合,修改初始数据,并且对于新数据而言偏差函数的总体梯度具有零值,指示互相抵消的正灵敏度和负灵敏度;
-从以下参数中选择子集合:相对于所述参数计算了灵敏度,所述子集合由以下参数形成:对所述参数而言灵敏度的符号相同;
-在将次级输出参数限制到构成子集合的部分的次级输出参数的同时,以迭代的方式重新开始所述方法;
-并且如此重复直到利用单一参数或小数量的参数定位异常。
10.根据权利要求1或2所述的方法,其特征在于,当单一初级输出参数呈现相对于理论值的测量值的异常时,优先考虑测量错误的假设。
11.根据权利要求1或2所述的方法,其特征在于,周期性地执行正向模型以及理论值与测量值的比较,并且当至少两个初级输出参数呈现相对于理论值的测量值的异常时,自动激活反向模型的执行。
12.根据权利要求1或2所述的方法,其特征在于,模型的元素还包括关于数据及参数的不确定性的指示,并且在比较期间考虑这些不确定性用于将输出参数的测量结果认为正常或相反地是异常的。
13.根据权利要求1至12中任一项所述的方法的使用,还包括改善配送网络的性能等级。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1457222 | 2014-07-25 | ||
FR1457222A FR3024254B1 (fr) | 2014-07-25 | 2014-07-25 | Procede de detection d'anomalies dans un reseau de distribution, en particulier distribution d'eau |
PCT/IB2015/055581 WO2016012971A1 (fr) | 2014-07-25 | 2015-07-23 | Procede de detection d'anomalies dans un reseau de distribution, en particulier distribution d'eau |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106575312A CN106575312A (zh) | 2017-04-19 |
CN106575312B true CN106575312B (zh) | 2020-06-26 |
Family
ID=52016705
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201580040183.2A Active CN106575312B (zh) | 2014-07-25 | 2015-07-23 | 用于检测配送网络特别是水配送网络中的异常的方法 |
Country Status (8)
Country | Link |
---|---|
US (1) | US10156465B2 (zh) |
EP (1) | EP3172632B1 (zh) |
CN (1) | CN106575312B (zh) |
AU (1) | AU2015293547B2 (zh) |
BR (1) | BR112017000075B8 (zh) |
CL (1) | CL2017000174A1 (zh) |
ES (1) | ES2721272T3 (zh) |
FR (1) | FR3024254B1 (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109325052B (zh) * | 2018-09-29 | 2021-06-04 | 创新先进技术有限公司 | 一种数据处理方法及装置 |
CN112836251A (zh) * | 2021-01-20 | 2021-05-25 | 合兴汽车电子股份有限公司 | 一种产品收缩率准确率的收集方法及系统 |
CN112597263B (zh) * | 2021-03-03 | 2021-05-11 | 浙江和达科技股份有限公司 | 管网检测数据异常判断方法及系统 |
CN113435795B (zh) * | 2021-08-26 | 2021-12-14 | 南方电网科学研究院有限责任公司 | 一种电力设备的状态评估方法、装置、设备及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005233609A (ja) * | 2005-04-28 | 2005-09-02 | Hitachi Ltd | 吸収式冷凍機の異常診断方法及びその装置 |
CN101988608A (zh) * | 2009-07-31 | 2011-03-23 | 法国燃气苏伊士集团 | 用于在流体传输管上分流的分支插接部的组件 |
FR2957097A1 (fr) * | 2010-03-04 | 2011-09-09 | Takadu Ltd | Systeme et procede de surveillance de ressources dans un reseau de distribution d'eau |
CN103155331A (zh) * | 2010-10-14 | 2013-06-12 | 西门子公司 | 用于产生状态信号的方法和装置 |
CN103590444A (zh) * | 2013-11-06 | 2014-02-19 | 中国水利水电科学研究院 | 大型埋地压力输水管道渗漏连续监测方法与爆管预警系统 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2841674C2 (de) * | 1978-09-25 | 1983-09-22 | Heide, Gerhard, Dipl.-Ing., 4006 Erkrath | Verfahren zur Überprüfung auf Leckverluste sowie dabei verwendbarer Meßschacht |
US4712182A (en) * | 1983-03-09 | 1987-12-08 | Hitachi, Ltd. | Method of estimating fracture point of pipe line network |
US6006252A (en) * | 1996-10-08 | 1999-12-21 | Wolfe; Mark A. | System and method for communicating information relating to a network resource |
DE10059567A1 (de) * | 2000-11-30 | 2002-06-13 | Siemens Ag | Verfahren und Vorrichtung zur Berechnung von Prozessgrößen eines industriellen Prozesses |
US10401164B2 (en) * | 2012-10-16 | 2019-09-03 | Exxonmobil Research And Engineering Company | Sensor network design and inverse modeling for reactor condition monitoring |
-
2014
- 2014-07-25 FR FR1457222A patent/FR3024254B1/fr active Active
-
2015
- 2015-07-23 BR BR112017000075A patent/BR112017000075B8/pt active IP Right Grant
- 2015-07-23 ES ES15762728T patent/ES2721272T3/es active Active
- 2015-07-23 EP EP15762728.2A patent/EP3172632B1/fr active Active
- 2015-07-23 AU AU2015293547A patent/AU2015293547B2/en active Active
- 2015-07-23 US US15/328,515 patent/US10156465B2/en active Active
- 2015-07-23 CN CN201580040183.2A patent/CN106575312B/zh active Active
-
2017
- 2017-01-23 CL CL2017000174A patent/CL2017000174A1/es unknown
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005233609A (ja) * | 2005-04-28 | 2005-09-02 | Hitachi Ltd | 吸収式冷凍機の異常診断方法及びその装置 |
CN101988608A (zh) * | 2009-07-31 | 2011-03-23 | 法国燃气苏伊士集团 | 用于在流体传输管上分流的分支插接部的组件 |
FR2957097A1 (fr) * | 2010-03-04 | 2011-09-09 | Takadu Ltd | Systeme et procede de surveillance de ressources dans un reseau de distribution d'eau |
CN103155331A (zh) * | 2010-10-14 | 2013-06-12 | 西门子公司 | 用于产生状态信号的方法和装置 |
CN103590444A (zh) * | 2013-11-06 | 2014-02-19 | 中国水利水电科学研究院 | 大型埋地压力输水管道渗漏连续监测方法与爆管预警系统 |
Also Published As
Publication number | Publication date |
---|---|
EP3172632B1 (fr) | 2019-01-23 |
FR3024254B1 (fr) | 2018-08-03 |
AU2015293547A1 (en) | 2017-02-02 |
AU2015293547B2 (en) | 2019-06-20 |
EP3172632A1 (fr) | 2017-05-31 |
CL2017000174A1 (es) | 2018-01-05 |
US20170205267A1 (en) | 2017-07-20 |
FR3024254A1 (fr) | 2016-01-29 |
CN106575312A (zh) | 2017-04-19 |
ES2721272T3 (es) | 2019-07-30 |
BR112017000075B1 (pt) | 2022-12-27 |
BR112017000075A2 (pt) | 2017-11-14 |
US10156465B2 (en) | 2018-12-18 |
BR112017000075B8 (pt) | 2023-03-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10387590B2 (en) | Techniques for iterative reduction of uncertainty in water distribution networks | |
Sarrate et al. | Sensor placement for leak detection and location in water distribution networks | |
Meseguer et al. | A decision support system for on-line leakage localization | |
Perez et al. | Leak localization in water networks: A model-based methodology using pressure sensors applied to a real network in Barcelona [applications of control] | |
Wu et al. | Water loss detection via genetic algorithm optimization-based model calibration | |
Soldevila et al. | Leak localization in water distribution networks using a mixed model-based/data-driven approach | |
Palau et al. | Burst detection in water networks using principal component analysis | |
EP3568682B1 (en) | Systems and methods for subnetwork hydraulic modeling | |
CN106575312B (zh) | 用于检测配送网络特别是水配送网络中的异常的方法 | |
CN111684238B (zh) | 用于实时发现和解决湿气文丘里流量计问题的系统和方法 | |
Moser et al. | Performance comparison of reduced models for leak detection in water distribution networks | |
CN109716090B (zh) | 一种用于估计配水管网的子网边界处的水流量的改进系统 | |
Marzola et al. | Leakage detection and localization in a water distribution network through comparison of observed and simulated pressure data | |
Vrachimis et al. | Leak detection in water distribution systems using hydraulic interval state estimation | |
Moors et al. | Automated leak localization performance without detailed demand distribution data | |
Takahashi et al. | Water distribution system model calibration under uncertainty environments | |
Okeya et al. | Locating pipe bursts in a district metered area via online hydraulic modelling | |
Cesario et al. | Calibrating water system models | |
JP2014054151A (ja) | 輸送ネットワークのグラフ構造構築装置、グラフ構造構築システム、グラフ構造構築方法及びグラフ構造構築プログラム | |
Sophocleous et al. | Leak detection and localization based on search space reduction and hydraulic modelling | |
상훈전 et al. | Data generation approaches to detect abnormal conditions in water distribution systems | |
JP6018970B2 (ja) | 配水制御装置および方法 | |
Przystałka et al. | Optimal placement of sensors and actuators for leakage detection and localization | |
Rosich et al. | Optimal sensor placement for leakage detection and isolation in water distribution networks | |
Pérez et al. | Leak monitoring |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240123 Address after: Paris France Patentee after: DEGREMONT Country or region after: France Address before: Paris France Patentee before: SUEZ GROUPE Country or region before: France |