基于DMA分区计量的实时产销差率估算方法

目前供水企业的产销差率和漏损率一般都以年为单位进行统计，这种计算方法有明显的滞后性，使得供水企业无法第一时间掌握产销差率实际情况。基于此，提出了实时产销差率的 概念，并给出了基于DMA分区计量的实时产销差率估算方法。实际案例分析结果表明，在不复杂的 DMA 分区背景下，基于 DMA 分区计量的实时产销差率可以直观地体现 DMA 分区管理质量，帮助供水企业了解该区域控漏效果，也提高了对漏水点的检出速度，节约了大量的水资源。

1.  1. 基于 DMA 分区的实时产销差率

      供水企业大多因为供水区域范围广、抄表反馈周期长，亦或是单双月抄表内容不同、季节气候不同，导致每月的产销差率并不能很好地反映真实的 供水管网情况。而且其波动范围较大，易对分析造 成困扰。在局限于现有管网监测设备的投入条件 下，找到一种最低成本的快速评估各 DMA 分区产销差率的方法，有利于供水企业获得更优的成本管 控。若能做到产销差率的实时估算对供水企业是一件很有价值的事情。而难点在于如何得到区域售水量。实时产销差则在于找到一个估算售水量的方法，同时凭借管网和水厂高频的物联采集设备得到 的供水量来实现区域产销差的实时评估。

1.  2. 实时产销差率的概念

       根据国际水协的标准水量平衡计算表，供水企业的产销差是指其供水量与售水量之差，产销差与供水量之比称为产销差率。引申到每个区域的产销差率，亦是如此。通过各 DMA 分区的供水量占比或压力突变可以判断分区内是否发生漏损，然而却无法通过实时评估区域内的产销差而评定该区域的管理质量。

       区域供水量可以通过分析管网结构，采用多组 流量仪计算获得，而区域售水量则不能很准确地得到。因此，研究区域售水量的估算方法很有必要。不考虑计量仪器，区域售水量变化可近似为区域内生活及生产用水量变化。

       首先对企业以往产销差率数据进行移动平均， 以降低各种噪声影响，得到一条相对平稳、贴近真实的产销差率曲线。移动平均法是一种简单的平滑时序数据方法，根据时间序列资料、逐项推移，依次计 算包含一定项数的时序平均值。移动平均计算 产销差率的滑动窗口长度一般采用 12 个月，这减少 了单双月抄收所产生的干扰，包含了一年时间内生 产、生活因素对供水数据产生的特有波动变化。

       假设产销差率用 N 表示，那么 1 － N 就是供水区域内总售水量的占比。则某时间段内供水企业总 供水量与 1 － N 的乘积即是该时间段内供水区域的 总售水量。通过DMA分区内的户表挂载，可以从历史抄 收数据中得到每个DMA分区的售水量，该区域的售水占比就是其售水量与全区域售水量之比。至此，供水区域的总售水量乘以该区域的售水占比，就 能得到目前该区域的售水量估计值，最终可以估计出该区域的产销差。如果出厂水量的采集频次足够密集，且 DMA分区流量仪数据与出厂水量同步采集，即能将分区内的产销差做到实时统计。综上，可以理解为实时产销差率即是为了解决传统产销差率 统计时效滞后、统计不完全而引入的一个新的成本 评估方式。后续假设数据采集频次为 1 次/h、移动平均产 销差率和区域售水占比的计算滑动窗口长度为 12 个月、区域移动平均售水估计量计算周期为 1 周 ( 考虑节假日、生产企业排班等因素，居民生活用水 量和企业生产用水量在 1 周内的变化是相对稳定的) ，进行实时产销差率的估算。

1.  3. 实时产销差率的估算

       首先采用 12 个月的出厂供水量与营收抄收水量计算移动平均产销差。假设供水企业每月抄收水量为 B( 单位为 m3 ) 、每月出厂供水量为 S( 单位为 m3 ) 、移动平均产销差率为N、m'为月份，则N的计算公式如下:

在没有发生特别重大的管网事件影响企业产销差率发生突变的情况下，移动平均产销差率 N 是相对稳定的，变化不大。如果有影响产销差率的突变因素发生，则需对 N 进行修正。

      假设某 DMA 分区的当月售水占比为 P。由于 供水企业无法避免转供水、偷盗追回等特殊用水情况带来的影响，且计算售水占比时没法剔除掉这些 因素，因此 P 也进行移动平均计算以减小噪声。由 于无法明确将追回水量划入现有 DMA分区，故将 其单独作为一个一级分区。不同转供水区域也单独 作为一个一级分区。表1 为供水企业 WW 连续 12 个月各 DMA 分区的移动平均售水占比。

表 1 某供水企业连续 12 个月的 DMA 分区售水占比

区域售水量每月统计 1 次，那么区域售水占比 数据是以月份为采样间隔的时间序列。以表 1 为 例，演示计算移动平均后的售水占比 P。假设一共 有 k 个 DMA 分区( 包括转供水和追回水量)，j 是分 区号( j = 1 ～ k) ; b 是区域售水量( 单位为m3 ) ; m 是 以月计的时间点，以 m 为下标的变量表示在时间点 m 的取值; t 是当前时间点( 当月或当时) ; i 为滑动 窗口的长度，i ＞0，本研究中 i 取 12。那么 P 的计算 公式为:

通过 2． 1 节的描述，将出厂实时供水量设为 x( 单位为 m3 /h) 。区域售水估计量 b ^ ( 单位为 m3 /h) 的计算公式为:

由于各区域用水量受到区域内生活、生产活动 的直接影响，为使区域售水估计量 b ^ 与当前情况的误差尽可能小，x 的采集频率设为每小时采集 1 次， 那么 b ^ 是以小时为采样间隔的时间序列。h 是以小时计的时间点，以 h 为下标的变量表示在时间点 h 的取值。当前时刻区域用水量的移动平均值 b( 单 位为 m3 /h) 计算公式为:

b 建议采用 1 周作为滑动窗口长度，这是因为在 1 周的用水周期内，区域内生活、生产用水量的变 化规律是相对稳定的; 小于 1 周会缺失部分工作、生 活日特有的用水特点，b 值将波动明显; 大于 1 周则 会拉长往期用水特征对当下 b 的影响，使得 b 值过分失准。区域供水量 s 可以直接通过 DMA 区域流量计得到，每小时采集 1 次，时间序列长度为 1 周。移动平均后的区域供水量 s( 单位为 m3 /h) 的计算公式 如下:

最终，区域实时产销差率 n 的估算公式为:

 4. 结论

① 鉴于传统计算产销差率的方法存在滞后 性，不能很好地反映供水企业的实际产销差情况，提出了实时产销差率概念，并给出了基于 DMA 分区计量的实时产销差率估算方法。案例分析表明，基于DMA 分区计量的实时产销差率，可以直观地体现 DMA分区管理质量，帮助供水企业了解该区域 控漏效果，提高供水企业管理绩效。

② 估算实时产销差率时，通过数据移动平均 的方法，可以降低各种因素对数据的影响波动，从而 得到一条相对平稳的产销差率曲线。③ 目前的计算公式是基于没有发生特别重大的管网事件影响企业产销差率发生突变，如果有影响产销差率的突变因素发生，则需对目前的计算公式进行修正，这也是后续的研究方向。