一种小区供水漏损量检测方法与流程

本发明属于供水漏损检测领域，具体是一种小区供水漏损量检测方法。

背景技术：

在大表远传技术的支持下，水司管理人员常常采用mnf法进行漏损自动监测，利用这些流量数据进行小区漏损情况评估，对小区供水漏损监管起到积极作用^[1][2]。实际应用中，需要对真实漏损进行较为精准的检测。由于受到夜间用水干扰、大表分钟级采样的局限，现有技术手段往往存在困难。

文献[3]利用μ-2σ表征剔除异常值后的dma最小流量，尝试消除实测数据的误差和干扰影响而得出小区的真实漏损水量。事实上，漏损流量并非恒定值，每天夜间流量的统计量μ和σ不同，每天的μ-2σ波动还不小，在实际应用中给水司管理人员带来困惑。文献[4]采用混合截尾正态分布的样本统计量及其导数轨迹偏离标准图方法来分离漏失量，算法复杂，特别是非稳定泄漏必须用到导数轨迹，偏差较大。

在实际应用中，针对小区供水非稳定泄漏问题，需要对正常用水流量和漏失流量进行快速、准确分离的检测方法。

[1]李岚,吴珊,寇晓霞,等.基于独立计量区的夜间最小流量的研究进展[j].给水排水,2018,44(6):135-141；

[2]徐哲,陈晖,何必仕.一种改进spc的小区供水漏损监测预警方法:cn201910699653.6[p].2019-12-13；

[3]郝志萍,侯煜堃,崔昌.计量小区(dma)夜间最小流量解析方法探讨与案例研究[j].南水北调与水利科技,2013(04):128-131+199；

[4]steveng.buchbergerandgayatrinadimpallil.leakestimationinwaterdistributionsystembystatisticalanalysisofflowreading[j].joumalofwaterresourcesplanningandmanagement,2004,130(4):321-329。

技术实现要素：

本发明针对小区供水中存在的非稳定泄漏问题，对高频率、高精度的流量计夜间实测流量数据进行用水流量和漏失流量分离，提出了一种小区供水漏损量检测方法，包括以下步骤：

步骤1、小区低用水前提下，连续快速采集入水口供水流量数据

在小区低用水时段，采用高频率、高精度流量计连续快速检测供水流量，采样精度1级以上，采样频率freq在0.1hz～10hz区间，一般选1hz，采样时长不小于1小时，选择其中供水量最小的1小时数据序列q(t)，0<t≤3600/freq。其中低用水时段，一般为夜间2:00-5:00。如小区有多个入口，检测时可以将次要入口关闭。

步骤2、计算小区供水流量序列的统计量和流量累积概率

计算供水流量序列q(t)的最大/最小值qmax/qmin，均值μ及标准差σ，在流量[qmin,qmax]区间计算分组流量频数f(bandi)、累积概率p(q)。这里，bandi为分组序列，q为流量变量(qmin≤q≤qmax)。

步骤3、检验小区供水流量数据q(t)的可分离性

针对分组流量频数f(bandi)，低通滤波平滑得到频数分布曲线，找到第一个波峰f(bandmax)和波谷f(bandmin)。若累积概率p(bandmin)>阈值pot(一般pot≥0.1)，则进入步骤4；否则，提高采样频率，或优化调整低用水检测时段，返回步骤1。

步骤4、优化小区分离正常用水流量和漏失流量

正常用水流量qu(t)与漏失流量ql(t)，是2个相互独立的随机变量，故假设：qu～n(μu,σu²)，ql～n(μl,σl²)，则供水流量q(t)：

q(t)＝qu(t)+ql(t)(1)

根据两个正态分布随机变量的线性组合仍服从正态分布性质，则：

μ＝μu+μl(2)

σ²＝σu²+σl²(3)

鉴于上述考虑，分离正常用水流量和漏失流量的优化求解如下：

步骤4.1以满足式(2)、(3)为优化目标，设立目标函数

则，优化目标为minj。

现设截断流量qt初始值为bandmax流量，迭代次数k＝1，迭代次数限值kt(100<kt<1000，一般取200)，迭代精度ε＝σ/100，j初始值j(0)＝u+σ；

步骤4.2利用截断流量qt将供水流量序列数据q(t)拆分成2部分，小于截断流量qt的为漏失流量为主和少量用水流量的混合流量q1’，大于截断流量qt的是用水流量和漏失流量的混合流量q2’。

步骤4.3以q1’作为漏失估计流量ql’，计算得到其统计量ul、σl²，；针对混合流量q2’，随机剔除漏失估计流量n(μl,σl²)，得到正常用水估计流量qu’，计算其统计量μu、σu²。

步骤4.4根据公式(4)计算目标函数值j(k)，比较并记录当前迭代次数下，jmin及μl‘、σl‘值。判断jmin<ε或迭代次数k≥kt？若满足，则步骤5，否则，qt增加δ(δ＝qmax/1000)，迭代次数k增1，返回步骤4.2。

步骤5、计算小区漏失流量

根据步骤4得到的μl‘、σl‘，按95％置信区间，计算得到漏失流量上下边界值ql＝μl‘±2σl‘。

本发明通过对高频率、高精度的流量计夜间实测流量数据进行用水流量和漏失流量分离，实现小区供水非稳定泄漏的漏损流量精准估计，完善了小区供水漏损量检测手段。

附图说明

图1：本发明方法流程示意图；

图2：实施例供水流量累积概率及频数分布图。

具体实施方式

下面通过实施例，并结合附图，对本发明的技术方案作进一步具体的说明。

实施例：本实施例的一种小区供水漏损量检测方法，如图1所示，包括以下步骤：

步骤1、小区低用水前提下，连续快速采集入水口供水流量数据

针对某小区，在夜间2:00-5:00低用水时段，采用1hz、j1等级(1级精度)的超声流量计，连续快速检测管径dn100的单入口供水流量，采样时长3小时，选择其中供水量最小的1小时数据序列q(t)，0<t≤3600。该小区供水流量累积概率及频数分布，如图2所示。

步骤2、计算小区供水流量序列的统计量和流量累积概率

计算供水流量序列q(t)的最大qmax＝67.74l/min，最小值qmin＝5.33l/min，均值μ＝36.16l/min，及标准差σ＝10.98l/min，在流量[5.33,67.74]区间计算分组流量频数f(bandi)、累积概率p(q)。这里，流量分组以0.5l/min为间隔，从5.0l/min至68l/min共分126组band，q为流量变量(5.33≤q≤67.74)。

步骤3、检验小区供水流量数据q(t)的可分离性

针对分组流量频数f(bandi)，低通滤波平滑得到频数分布曲线，找到第一个波峰f(band9)和波谷f(band17)。此时，此时，band9＝9.5l/min，band17＝13.5l/min，累积概率p(13.5)＝0.11>阈值pot(pot取0.1)，则进入步骤4；否则，提高采样频率，或优化调整低用水检测时段，返回步骤1。

步骤4、优化小区分离正常用水流量和漏失流量

分离正常用水流量和漏失流量的优化求解方法如下：

步骤4.1设立目标函数优化目标为minj。

现设截断流量qt初始值为band8流量9l/min，迭代次数k＝1，迭代次数限值kt＝200，迭代精度ε＝σ/100＝0.1098，j初始值j(0)＝μ+σ＝47.14。

步骤4.2利用截断流量qt将供水流量序列数据q(t)拆分成2部分，小于截断流量qt的为漏失流量为主和少量用水流量的混合流量q1’，大于截断流量qt的是用水流量和漏失流量的混合流量q2’。

步骤4.3以q1’作为漏失估计流量ql’，计算得到其统计量μl、σl²，；针对混合流量q2’，随机剔除漏失估计流量n(ul,σl²)，得到正常用水估计流量qu’，计算其统计量μu、σu²。

为了体现随机剔除，混合流量q2’可以选择多次随机产生漏失估计流量n(μl,σl²)进行减除，取σu最小的那一组。

步骤4.4根据公式(4)计算目标函数值j(k)，比较并记录当前迭代次数下，jmin及μl‘、σl‘值。判断jmin<ε或迭代次数k≥kt？若满足，则步骤5，否则，qt增加δ＝qmax/1000＝0.06774，迭代次数增1，返回步骤4.2。

经过上述步骤4.1-4.4循环计算，得到优化结果：jmin＝0.1009，此时，μl‘＝8.85l/min、σl‘＝1.08l/min。

步骤5、计算小区漏失流量

根据步骤4得到的μl‘、σl‘，按95％置信区间，计算漏失流量上下边界值ql＝μl‘±2σl，得到漏失流量区间[6.69，11.02]。

经过利用该小区用户智能分表同期实测，得到正常用水量qu(t)，用供水流量q(t)减去正常用水量qu(t)，得到漏失流量ql(t)，再计算得到ql～n(8.90,1.14)，其95％置信区间为[6.62,11.18]。可见运用本发明方法的实例结果占实际漏失流量区间的94.96％，具有较高检测精度。

本文中所描述的具体实施例仅仅是对本发明方法作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。