一种供水预测方法及系统的制作方法
【专利摘要】本发明公开了一种供水预测方法及系统，方法包括：分析历史数据并建立供水预测模型，得到输入层参数与输出层输出的供水量之间的函数关系；根据建立的供水预测模型和所需预测日期进行实时供水量预测；分析历史数据并建立供水预测模型的过程包括但不限于采用移动平均法对输入的历史数据进行预处理过程，采用神经网络算法确定供水预测模型的网络拓扑结构过程，根据预处理后的历史数据以及相应的网络拓扑结构确定供水预测模型过程。本发明综合采用了移动平均法和神经网络算法来进行供水量预测，计算的复杂度低，具有明显的周期性和强非线性，实时性好，能准确描述各种因素对供水量的影响，能加快神经网络模型的收敛速度，可广泛应用于供水分析领域。
【专利说明】
-种供水预测方法及系统
技术领域
[0001 ]本发明设及供水分析领域，尤其是一种供水预测方法及系统。
【背景技术】
[0002] 城市供水量预测是指利用过去已有的信息去推求未来某一段时间的城市供水负 荷。城市供水量预测是城市供水系统规划建设和改造优化的重要前期工作，它的准确性与 否将直接影响到工程量大小、工程费用、调度决策的可靠性和实用性。
[0003] 城市供水量预测从时间上分类大致分为两种:短期供水量预测及中长期供水量预 。前者是指根据历史的供水量数据结合影响短期供水量的各种因素对未来较短时间段的 供水量进行预测（如小时预测、日预测或周预测），短期预测的目的是为供水系统的优化调 度管理服务;后者是指根据历史的供水量数据结合影响中长期供水量的各种因数对未来一 年甚至几年的城市供水量进行预测，中长期预测主要为未来长期的规划和城市供水管网改 造服务。城市供水量预测从空间上分类大致分为两种:城市总供水量预测及节点供水量预 测，前者服务于宏观调度;后者服务为城市供水量分配的微观调度。
[0004] 优化调度供水量预测的方法一般分为两大类:解释性预测方法(即回归分析方法） 和时间序列分析方法，前者认为系统的输入参数与输出结果之间存在着某种复杂的因果关 系，故W此关系来构造预测模型进行预测，该类模型对输入参数的精度及可靠性要求较高； 后者是把系统看作一个暗箱，可W不管其影响因素，而只关屯、预测和预测的结果，其预测过 程只依赖于历史观测数据。
[0005] 供水量预测的方法有很多，比较常用的有：指数平滑模型、移动平均（即MA)模型、 自回归（即AR)预测模型、灰色预测模型和神经网络预测等。然而，现有的运些供水量预测方 法存在着较大的缺陷或不足:指数平滑模型需要的数据量大，计算复杂高;移动平均(M)模 型在周期的选择上存在一定的困难性，非线性的特征被一定程度的削弱；自回归(AR)预测 模型具有明显的滞后性，当实际数据发生异常变化时，预测数据因模型的平滑作用而无法 反应出异常变化，使得其在预测一些异常值时造成较大误差甚至失真;灰色预测模型很难 准确描述各种因素对供水量的影响；传统的神经网络模型收敛慢，容易陷入局部最小值的 情形。
【发明内容】
[0006] 为解决上述技术问题，本发明的目的在于:提供一种计算复杂度低、实时、预测精 度高、能准确描述各种因素和收敛快的，对供水量的影响供水预测方法。
[0007] 本发明的另一目的在于:提供一种计算简单、实时、预测精度高、能准确描述各种 因素和收敛快的，供水预测系统。
[000引本发明所采取的技术方案是：
[0009] -种供水预测方法，包括W下步骤：
[0010] S1、分析历史数据并建立供水预测模型，得到输入层参数与输出层输出的供水量 之间的函数关系；
[0011] S2、根据建立的供水预测模型和所需预测日期进行实时供水量预测；
[0012] 所述分析历史数据并建立供水预测模型的过程包括但不限于采用移动平均法对 输入的历史数据进行预处理过程，采用神经网络算法确定供水预测模型的网络拓扑结构过 程，根据预处理后的历史数据W及相应的网络拓扑结构确定供水预测模型过程。
[OOU] 进一步，所述步骤Sl包括；
[0014] S11、根据供水过程的特性确定供水预测模型的输入参数；
[0015] S12、采用移动平均法对输入的历史数据进行预处理；
[0016] S13、采用S层BP神经网络模型建立供水预测模型的网络拓扑结构；
[0017] S14、确定S层BP神经网络模型的计算参数；
[0018] S15、根据预处理后的历史数据W及S层BP神经网络模型相应的网络拓扑结构进 行训练，然后根据训练的结果确定供水预测模型。
[0019] 进一步，所述步骤S12包括：
[0020] S121、对初始样本中输入的历史数据进行平滑处理，得到去除突变异常点后的样 本数据，所述突变异常点满足：
[0021]
[0022] 其中，Qin为第n个样本的历史供水量值，n=l，2,…，N;N为最大样本数，O(Q)
[0023] 为历史供水量样本的均方差，O(Q)的计算公式为：
[0024]
[0025] S122、对去除突变异常点后的样本数据中小范围内的波动点和斜率突然变化的样 本点进行优化，得到优化后的样本数据；
[0026] S123、对优化后的样本数据进行移动平均处理，所述移动平均处理所采用的公式 为：
[0027]
[0028] 其中，耐为优化后的样本数据中第9个样本的样本值，^^为优化后的样本数据中第 q个样本的移动平均值，h为移动平均的周期数；
[0029] S124、对移动平均处理后的所有样本数据进行归一化处理，所述归一化处理所采 用的公式为：
[0030]
[0031] 其中，巧^ >分别为移动平均处理后的所有样本中第j个字段归一化处理前的 值和归一化处理后的值，分别表示所有样本中第j个字段归一化处理前的最大 值和归一化处理后的最大值，。分别表示所有样本中第j个字段归一化处理前的 最小值和归一化处理后的最小值，if >的取范围为[0，1]。
[0032] 进一步，所述供水预测模型的输入参数包括但不限于历史用水量、天气情况、溫 度、湿度、节假日情况、工业情况、城镇居民生活条件、人口和GDP中的任意一种或多种，每个 参数构成输入层的一个节点，输出层输出第二天的供水量，并采用相关性分析的方法对输 入参数进行分析，计算各参数与供水量的之间的相关系数及各参数对供水量的贡献率，然 后选取与供水量相关系数高的参数作为影响日供水量的主要因子，并设定阔值作为神经网 络流出神经元的控制量;所述神经网络算法采用神经元为基本处理单元，所述神经元为多 输入、单输出的非线性器件;所述S层BP神经网络模型的计算参数包括但不限于各层神经 元之间的权值、神经元中的激励函数及数据流出神经元的阔值，所述各层神经元之间的权 值W及数据流出神经元的阔值由random函数随机生成的20种初始的权值矩阵和阔值输入 S层BP神经网络模型计算而来，所述初始的权值和阔值均在(0，1)之间随机取值，所述神经 元中的激励函数采用Log-S i gmo i d函1
[0033] 进一步，所述供水预测模型WS层BP神经网络作为网络拓扑结构，所述输入层与 输出层之间设置有隐含层，所述=层BP神经网络的输入层数目m、输出层数目m和隐含层数 目02满足:n2 = sqd(ni+m)+a，其中，sqd函数为平方根函数，a为常数，且a取1到10的自然 数;所述隐含层的节点根据输入层的节点数确定并通过增减节点数进行训练对比，所述神 经网络算法中每层的输出函数为输入函数关于该层的激发函数的函数;若所述输出层的输 出结果与期望值不符，将输出层的误差转化为各层各个节点的连接权值与阔值的误差，将 输出层的误差逐层向输入层反向传递，把误差分摊到各个节点，计算各个节点的参考误差， 并调整各个节点的连接权值及阔值大小，直到误差最小。
[0034] 进一步，所述步骤Sl 5具体为：
[0035] 根据预处理后的历史数据W及S层BP神经网络模型相应的网络拓扑结构采用训 练精度模式、学习次数模式或者两种模式的结合进行训练，然后根据训练的结果确定供水 预测模型；
[0036] 其中，采用训练精度模式进行训练包括如下步骤:根据S层BP神经网络的网络拓 扑结构计算网络误差，计算训练过程的总误差，判断总误差是否达到设定精度，当判断没有 达到精度时重复上述步骤直至达到设定精度，根据网络误差修改网络权重并完成训练;所 述采用学习次数模式进行训练包括如下步骤:根据S层BP神经网络模型的网络拓扑结构计 算网络误差，根据计算的网络误差修改网络权重，将训练次数加1，判断是否达到最大训练 次数，若判断没有达到最大训练次数则重复上述步骤直至达到最大训练次数完成训练;采 用训练精度模式结合学习次数模式进行训练包括如下步骤:根据S层BP神经网络模型的网 络拓扑结构正向传播计算网络误差，计算训练过程的总误差，判断总误差是否达到设定精 度，判断是否达到最大训练次数，当判断没有达到精度时或没有达到最大训练次数时重复 上述步骤直至达到设定精度和最大训练次数，根据网络误差修改网络权重；
[0037] 所述训练过程的总误差为所有含有预处理后的历史数据样本误差的均方差之和；
[0038] 所述训练过程采用改进的变步长梯度下降法修正各个节点的权值及误差，各个节 点当前时刻的权值修正量A w(t)为：A w(t)=riS+a A w(t-l)，其中，A w(t-l)为各个节点上 一时刻的权值修正量，S为各个节点理论计算得到的权值修正量，a为动量系数，且a G (0， i)，n为自适应调节学习率，n的调整方式为:先设定初始步长，若进行一次迭代后节点的误 差增大，则将步长和n均乘W小于1的调节因子e，然后沿原方向重新进行下一次迭代;若进 行一次迭代后节点的误差减小，则将步长和n均乘W大于1的调节因子0，然后沿原方向重新 进行下一次迭代。
[0039] 进一步，所述步骤S2包括：
[0040] S21、输入预测所需的参数。
[0041] S22、判断所需预测日期是否与上一次预测日期一致，若是，则根据输入的参数和 上一次预测日期的供水预测模型正向计算所需预测日期的供水量，反之，则在执行步骤S15 更新供水预测模型后正向计算所需预测日期的供水量；
[0042] S23、判断所需预测日期是否为节假日，若是，则W历史当日的数据作为参考数据 重新预测所需预测日期的供水量，并将重新预测的供水量进行输出；反之，则直接输出步骤 S22计算出的预测供水量。
[0043] 本发明所采取的另一技术方案是：
[0044] -种供水预测系统，包括：
[0045] 建模模块，用于分析历史数据并建立供水预测模型，得到输入层参数与输出层输 出的供水量之间的函数关系；
[0046] 预测模块，用于根据建立的供水预测模型和所需预测日期进行实时供水量预测；
[0047] 所述分析历史数据并建立供水预测模型的过程包括但不限于采用移动平均法对 输入的历史数据进行预处理过程，采用神经网络算法确定供水预测模型的网络拓扑结构过 程，根据预处理后的历史数据W及相应的网络拓扑结构确定供水预测模型过程；
[0048] 所述建模模块的输出端与预测模块的输入端连接。
[0049] 进一步，所述建模模块包括：
[0050] 输入参数确定单元，用于根据供水过程的特性确定供水预测模型的输入参数；
[0051] 预处理单元，用于采用移动平均法对输入的历史数据进行预处理；
[0052] 拓扑结构建立单元，用于采用S层BP神经网络模型建立供水预测模型的网络拓扑 结构；
[0053] 计算参数确定单元，用于确定S层BP神经网络模型的计算参数；
[0054] 训练单元，用于根据预处理后的历史数据W及S层BP神经网络模型相应的网络拓 扑结构进行训练，然后根据训练的结果确定供水预测模型；
[0055] 所述输入参数确定单元的输出端依次通过预处理单元、拓扑结构建立单元、计算 参数确定单元和训练单元进而与预测模块的输入端连接。
[0化6] 进一步，所述预处理单元包括：
[0057]平滑处理子单元，用于对初始样本中输入的历史数据进行平滑处理，得到去除突 变异常点后的样本数据，所述突变异常点满足：
[0化8;
[0059] 其中，Qin为第n个样本的历史供水量值，n=l，2,…，N;N为最大样本数，O(Q)
[0060] 为历史供水量样本的均方差，O(Q)的计算公式为：
[0061]
[0062] 优化子单元，用于对去除突变异常点后的样本数据中小范围内的波动点和斜率突 然变化的样本点进行优化，得到优化后的样本数据；
[0063] 移动平均处理子单元，用于对优化后的样本数据进行移动平均处理，所述移动平 配巧田的A井电.
[0064]
[0065] 其中，Xq为优化后的样本数据中第q个样本的样本值，^为优化后的样本数据中第 q个样本的移动平均值，h为移动平均的周期数；
[0066] 归一化处理子单元，用于对移动平均处理后的所有样本数据进行归一化处理，所
述归一A巧航要田的公才兩?
[0067]
[0068] 其中，和文^/>分别为移动平均处理后的所有样本中第j个字段归一化处理前 的值和归一化处理后的值，Xjmax和分别表示所有样本中第j个字段归一化处理前的最 大值和归一化处理后的最大值，分别表示所有样本中第j个字段归一化处理前 的最小值和归一化处理后的最小值，的取范围为[0，1];
[0069] 所述输入参数的输出端确定单元依次通过平滑处理子单元、优化子单元、移动平 均处理子单元和归一化处理子单元进而与拓扑结构建立单元的输入端连接。
[0070] 本发明的方法的有益效果是:分析历史数据并建立供水预测模型的过程包括但不 限于采用移动平均法对输入的历史数据进行预处理过程W及采用神经网络算法确定供水 预测模型的网络拓扑结构过程，综合采用了移动平均法和神经网络算法来进行供水量预 ，能够利用较少的历史数据精确预测出未来的供水量，计算的复杂度低;结合神经网络算 法来进行预测，有效克服了单一移动平均算法周期选择困难和非线性特征弱的缺点，具有 明显的周期性和强非线性;采用了移动平均法，去掉了部分噪音，提高了样本的精度;采用 了神经网络算法，能实时反应出实际数据的异常变化，实时性好，且能准确描述各种因素对 供水量的影响;结合了移动平均法，能加快神经网络模型的收敛速度。
[0071] 本发明的系统的有益效果是:建模模块分析历史数据并建立供水预测模型的过程 包括但不限于采用移动平均法对输入的历史数据进行预处理过程W及采用神经网络算法 确定供水预测模型的网络拓扑结构过程，综合采用了移动平均法和神经网络算法来进行供 水量预测，能够利用较少的历史数据精确预测出未来的供水量，计算的复杂度低;结合神经 网络算法来进行预测，有效克服了单一移动平均算法周期选择困难和非线性特征弱的缺 点，具有明显的周期性和强非线性;采用了移动平均法，去掉了部分噪音，提高了样本的精 度;采用了神经网络算法，能实时反应出实际数据的异常变化，实时性好，且能准确描述各 种因素对供水量的影响;结合了移动平均法，能加快神经网络模型的收敛速度。
【附图说明】
[0072] 图1为本发明一种供水预测方法的整体流程图；
[0073] 图2为本发明一种供水预测系统的结构框图；
[0074] 图3为本发明实施例二进行建模的流程图；
[0075] 图4为本发明的实施例二采用训练精度模式结合训练次数进行训练的流程图；
[0076] 图5为本发明实施例二根据训练得到的预测模型进行供水预测的流程图。
【具体实施方式】
[0077] 参照图1，一种供水预测方法，包括W下步骤：
[0078] S1、分析历史数据并建立供水预测模型，得到输入层参数与输出层输出的供水量 之间的函数关系；
[0079] S2、根据建立的供水预测模型和所需预测日期进行实时供水量预测；
[0080] 所述分析历史数据并建立供水预测模型的过程包括但不限于采用移动平均法对 输入的历史数据进行预处理过程，采用神经网络算法确定供水预测模型的网络拓扑结构过 程，根据预处理后的历史数据W及相应的网络拓扑结构确定供水预测模型过程。
[0081] 进一步作为优选的实施方式，所述步骤Sl包括：
[0082] S11、根据供水过程的特性确定供水预测模型的输入参数；
[0083] S12、采用移动平均法对输入的历史数据进行预处理；
[0084] S13、采用S层BP神经网络模型建立供水预测模型的网络拓扑结构；
[0085] S14、确定S层BP神经网络模型的计算参数；
[0086] S15、根据预处理后的历史数据W及S层BP神经网络模型相应的网络拓扑结构进 行训练，然后根据训练的结果确定供水预测模型。
[0087] 进一步作为优选的实施方式，所述步骤S12包括：
[0088] S121、对初始样本中输入的历史数据进行平滑处理，得到去除突变异常点后的样 本数据，所述突变异常点满足：
[0089；
[0090]其中，Qin为第n个样本的历史供水量值，n=l，2,…，N;N为最大样本数，O(Q)为历史 供水量样本的均方差，O(Q)的计算公式为：
[0091]
[0092] S122、对去除突变异常点后的样本数据中小范围内的波动点和斜率突然变化的样 本点进行优化，得到优化后的样本数据；
[0093] S123、对优化后的样本数据进行移动平均处理，所述移动平均处理所采用的公式 为：
[0094]
[00M]其中，Xq为优化后的样本数据中第q个样本的样本值，;^为优化后的样本数据中第 q个样本的移动平均值，h为移动平均的周期数；
[0096] S124、对移动平均处理后的所有样本数据进行归一化处理，所述归一化处理所采 用的公式为：
[0097]
[0098] 其中，XJW和分别为移动平均处理后的所有样本中第j个字段归一化处理前的 值和归一化处理后的值，Xjmax和スym。X分别表示所有样本中第j个字段归一化处理前的最大 值和归一化处理后的最大值，Xjmin和分别表示所有样本中第j个字段归一化处理前的 最小值和归一化处理后的最小值，.的取范围为[0，1]。
[0099] 进一步作为优选的实施方式，所述供水预测模型的输入参数包括但不限于历史用 水量、天气情况、溫度、湿度、节假日情况、工业情况、城镇居民生活条件、人口和GDP中的任 意一种或多种，每个参数构成输入层的一个节点，输出层输出第二天的供水量，并采用相关 性分析的方法对输入参数进行分析，计算各参数与供水量的之间的相关系数及各参数对供 水量的贡献率，然后选取与供水量相关系数高的参数作为影响日供水量的主要因子，并设 定阔值作为神经网络流出神经元的控制量;所述神经网络算法采用神经元为基本处理单 元，所述神经元为多输入、单输出的非线性器件;所述S层BP神经网络模型的计算参数包括 但不限于各层神经元之间的权值、神经元中的激励函数及数据流出神经元的阔值，所述各 层神经元之间的权值W及数据流出神经元的阔值由random函数随机生成的20种初始的权 值矩阵和阔值输入S层BP神经网络模型计算而来，所述初始的权值和阔值均在(0，1)之间 随机取值，所述神经元中的激励函数采用Log-Sigmoid函数:./'W =
[0100] 进一步作为优选的实施方式，所述供水预测模型WS层BP神经网络作为网络拓扑 结构，所述输入层与输出层之间设置有隐含层，所述=层BP神经网络的输入层数目m、输出 层数目m和隐含层数目m满足:n2 = sqd(m+m)+a，其中，sqrt函数为平方根函数，a为常数， 且a取1到10的自然数;所述隐含层的节点根据输入层的节点数确定并通过增减节点数进行 训练对比，所述神经网络算法中每层的输出函数为输入函数关于该层的激发函数的函数； 若所述输出层的输出结果与期望值不符，将输出层的误差转化为各层各个节点的连接权值 与阔值的误差，将输出层的误差逐层向输入层反向传递，把误差分摊到各个节点，计算各个 节点的参考误差，并调整各个节点的连接权值及阔值大小，直到误差最小。
[0101] 进一步作为优选的实施方式，所述步骤S15具体为：
[0102] 根据预处理后的历史数据W及S层BP神经网络模型相应的网络拓扑结构采用训 练精度模式、学习次数模式或者两种模式的结合进行训练，然后根据训练的结果确定供水 预测模型；
[0103] 其中，采用训练精度模式进行训练包括如下步骤:根据S层BP神经网络的网络拓 扑结构计算网络误差，计算训练过程的总误差，判断总误差是否达到设定精度，当判断没有 达到精度时重复上述步骤直至达到设定精度，根据网络误差修改网络权重并完成训练;所 述采用学习次数模式进行训练包括如下步骤:根据S层BP神经网络模型的网络拓扑结构计 算网络误差，根据计算的网络误差修改网络权重，将训练次数加1，判断是否达到最大训练 次数，若判断没有达到最大训练次数则重复上述步骤直至达到最大训练次数完成训练;采 用训练精度模式结合学习次数模式进行训练包括如下步骤:根据S层BP神经网络模型的网 络拓扑结构正向传播计算网络误差，计算训练过程的总误差，判断总误差是否达到设定精 度，判断是否达到最大训练次数，当判断没有达到精度时或没有达到最大训练次数时重复 上述步骤直至达到设定精度和最大训练次数，根据网络误差修改网络权重；
[0104] 所述训练过程的总误差为所有含有预处理后的历史数据样本误差的均方差之和；
[0105] 所述训练过程采用改进的变步长梯度下降法修正各个节点的权值及误差，各个节 点当前时刻的权值修正量A w(t)为：A w(t)=riS+a A w(t-l)，其中，A w(t-l)为各个节点上 一时刻的权值修正量，S为各个节点理论计算得到的权值修正量，a为动量系数，且a G (0， i)，n为自适应调节学习率，n的调整方式为:先设定初始步长，若进行一次迭代后节点的误 差增大，则将步长和n均乘W小于1的调节因子e，然后沿原方向重新进行下一次迭代;若进 行一次迭代后节点的误差减小，则将步长和n均乘W大于1的调节因子0，然后沿原方向重新 进行下一次迭代。
[0106] 进一步作为优选的实施方式，所述步骤S2包括：
[0107] S21、输入预测所需的参数。
[0108] S22、判断所需预测日期是否与上一次预测日期一致，若是，则根据输入的参数和 上一次预测日期的供水预测模型正向计算所需预测日期的供水量，反之，则在执行步骤S15 更新供水预测模型后正向计算所需预测日期的供水量；
[0109] S23、判断所需预测日期是否为节假日，若是，则W历史当日的数据作为参考数据 重新预测所需预测日期的供水量，并将重新预测的供水量进行输出；反之，则直接输出步骤 S22计算出的预测供水量。
[0110] 参照图2，一种供水预测系统，包括：
[0111] 建模模块，用于分析历史数据并建立供水预测模型，得到输入层参数与输出层输 出的供水量之间的函数关系；
[0112] 预测模块，用于根据建立的供水预测模型和所需预测日期进行实时供水量预测；
[0113] 所述分析历史数据并建立供水预测模型的过程包括但不限于采用移动平均法对 输入的历史数据进行预处理过程，采用神经网络算法确定供水预测模型的网络拓扑结构过 程，根据预处理后的历史数据W及相应的网络拓扑结构确定供水预测模型过程；
[0114] 所述建模模块的输出端与预测模块的输入端连接。
[0115] 参照图2,进一步作为优选的实施方式，所述建模模块包括：
[0116] 输入参数确定单元，用于根据供水过程的特性确定供水预测模型的输入参数；
[0117] 预处理单元，用于采用移动平均法对输入的历史数据进行预处理；
[0118] 拓扑结构建立单元，用于采用S层BP神经网络模型建立供水预测模型的网络拓扑 结构；
[0119] 计算参数确定单元，用于确定S层BP神经网络模型的计算参数；
[0120] 训练单元，用于根据预处理后的历史数据W及S层BP神经网络模型相应的网络拓 扑结构进行训练，然后根据训练的结果确定供水预测模型；
[0121] 所述输入参数确定单元的输出端依次通过预处理单元、拓扑结构建立单元、计算 参数确定单元和训练单元进而与预测模块的输入端连接。
[0122] 参照图2,进一步作为优选的实施方式，所述预处理单元包括：
[0123] 平滑处理子单元，用于对初始样本中输入的历史数据进行平滑处理，得到去除突 变异常点后的样本数据，所述突变异常点满足：
[0124]
[0125] 其中，Qin为第n个样本的历史供水量值，n=l，2,…，N;N为最大样本数，O(Q)
[0126] 为历史供水量样本的均方差，O(Q)的计算公式为：
[0127]
[0128] 优化子单元，用于对去除突变异常点后的样本数据中小范围内的波动点和斜率突 然变化的样本点进行优化，得到优化后的样本数据；
[0129] 移动平均处理子单元，用于对优化后的样本数据进行移动平均处理，所述移动平 均处理所采用的公式为：
[0130；
[0131] 其中，Xq为优化后的样本数据中第q个样本的样本值，^为优化后的样本数据中第 q个样本的移动平均值，h为移动平均的周期数；
[0132] 归一化处理子单元，用于对移动平均处理后的所有样本数据进行归一化处理，所 述归一化处理所采用的公式为：
[0133]
[0134] 其中，分别为移动平均处理后的所有样本中第j个字段归一化处理前 的值和归一化处理后的值，?max和?.lna)c分别表示所有样本中第j个字段归一化处理前的最 大值和归一化处理后的最大值，刮。1。和^1。山分别表示所有样本中第j个字段归一化处理前 的最小值和归一化处理后的最小值的取范围为[0，1];
[0135] 所述输入参数的输出端确定单元依次通过平滑处理子单元、优化子单元、移动平 均处理子单元和归一化处理子单元进而与拓扑结构建立单元的输入端连接。
[0136] 下面结合说明书附图和具体实施例对本发明作进一步详细说明。
[0137] 实施例一
[0138] 为了克服现有供水预测所存在的问题，本发明提出了一种由移动平均法与神经网 络预测方法相结合而成的供水预测方法，该方法能够较精确反映出城市供水的规律性;能 够利用较少的历史数据精确预测出未来的供水量，同时在外界条件变化的情况下能保持较 高的预测精度且能提供自动更新校核的智能供水量预测模型。该方法主要包括建模步骤和 实时预测步骤运两大部分。
[0139] 其中，建模步骤分析历史数据并建立供水预测模型，建立输入层参数与输出层输 出的供水量之间的函数关系。实时预测步骤根据建模步骤建立的供水预测模型W及输入参 数实时预测目标预测日期的供水量。
[0140] 建模步骤包括:根据供水过程的特性确定供水预测模型的输入参数、对输入的历 史数据预处理、确定BP神经网络的拓扑结构，确定供水预测模型的计算参数，根据输入历史 数据结合相应的网络拓扑结构进行训练，并根据训练结果确定供水预测模型。
[0141] (1)确定供水预测模型的输入参数。
[0142] 建模步骤的输入参数包括历史用水量、天气情况、溫度、湿度、节假日情况、工业情 况，城镇居民生活条件，人口和GDP中的任意一种或多种。采用相关性分析的方法对上述输 入因子进行分析，计算各因素与供水量间的相关系数及各因素对供水量的贡献率，并选取 与供水量相关系数较高的因素作为影响日供水量的主要因子。每个参数构成输入层的一个 节点，输出层输出预测日期（如第二天）的供水量，神经网络算法采用神经元为基本处理单 元，神经元为多输入、单输出的非线性器件。
[0143] (2)输入的历史数据的预处理。
[0144] 建模步骤中对输入的历史数据预处理尤为关键。初始样本中往往会渗杂一些与研 究问题不相关的因素、相关性不大的因素，运些因素往往会对初始样本造成影响，影响了模 型的精度和稳定性。对于供水预测问题来说，历史流量对外界最为敏感，受异常情况的影响 最为明显，故本发明采用改进的移动平均法对历史流量数据进行预处理，在保证样本真实 性的同时排除一些系统误差和记录误差，提高样本精度。本发明对历史流量数据进行预处 理步骤包括：
[0145] 1)对历史的供水量资料进行平滑处理，W排除偶然因素造成的突变异常点。平滑
[0146] 处理采用W下的3倍均方差方法：
[0147]
[0148] 其中，Qin为历史流量值，将符合公式(2)条件的点从样本中除去，W剩下的样本作 为训练样本输入神经网络模型进行训练。
[0149] 2)对于小范围的波动点和斜率突然变化的样本点，本发明依据厂区的实际运行情 况对其进行优化，W提高输入样本的精度。
[0150] 3)移动平均处理:将去掉异常的数据进行移动平均处理，去掉一些白噪声，W提高 样本的稳定性。移动平均处理所采用的公式为：
[0151]
[0152] ^^为样本的移动平均值，N为最大样本数，h为移动平均的周期数，Xq为去掉异常点 后的样本值。
[0153] 4)将供水量资料作归一化处理。由于输入参数的量级不一致，需要做归一化处理 才能正确反映出输入参数对输出量的影响过程，本发明的归一化过程采用如下式的线性归
一法：
[0154]
[0155] 其中诚min表示所有输入样本第j个字段值的最小值诚max表示所有输入样本第j个 字段值的最大值，将所有样本的第j个字段值归一化到中，范围取[0，1]。
[0156] (3)建立BP神经网络模型，确定BP神经网络的拓扑结构。
[0157] S层BP神经网络可W逼近任意非线性映射关系，而神经网络算法中增加神经网络 的层数可W更进一步降低误差，提高精度，但每增加一层，网络就会变得更加复杂，大大增 加了网络的训练时间。而预测精度的提高实际上也可W通过增加隐含层节点数来实现，而 且其训练的效果也比增加层数的方法更容易观察和调整。因此相比于增加层数的方法，增 加隐含层节点数的方法更值得优先考虑。
[0158] 理论上，隐含层的节点数量可W为任意数值，但在实际的应用中，过多的隐含层节 点会导致一系列的问题。隐含层节点过多会使网络的规模变大，增加网络的复杂性:一方 面，若网络规模相对于训练样本过于庞大，会出现训练结果"太好"的现象，看似训练精度很 高，实际上已经出现了 "过拟合"的现象（即陷入局部极小点而得不到最优点的现象），使得 网络的预测能力大大降低，无法达到预期的效果;另一方面，隐含层节点过多必然会导致更 多的待定参数，网络规模的扩大也给网络的训练带来不小的难度，增加了训练的时间，且降 低了收敛的速度。
[0159] 本发明的供水预测模型的网络拓扑结构采用了 S层BP神经网络，输入层与输出层 之间设置有隐含层;隐含层的节点根据输入层的节点数确定并通过增减节点数进行训练对 比，神经网络算法中每层的输出函数为输入函数关于该层的激发函数的函数;若所述输出 层的输出结果与期望值不符，将输出层的误差转化为各层各个节点的连接权值与阔值的误 差，将输出层的误差逐层向输入层反向传递，把误差分摊到各个节点，计算各个节点的参考 误差，调整各个节点的连接权值及阔值大小，直到误差最小。
[0160] 在BP神经网络模型中，假设中间隐含层的数目为m，输入层数目为m，输出层数目 为m，则运S者的关系满足下式：
[0161] n2 = sqrt(ni+m)+a；
[0162] a = 1-10
[0163] 适当增加和减少输入层、输出层、隐含层的个数组合成多种神经网络结构。如计算 得到输入层、输出层、隐含层的大致数目为7、7、7，则根据试错法对运=层数目进行微调，将 =层数目分别取7、8、9进行组合可生成27种组合;接着，可通过分别尝试法从运27种组合中 找到最为合适的神经网络结构。在本实施例中，当输入层输出层固定时，主要改变的是隐含 层个数。经分析查找后的结果（即分别尝试法的结果)发现，当隐含层数为別寸，模型收敛及 预测能力较好。
[0164] (4)确定供水预测模型的计算参数。
[0165] 神经网络模型中主要的计算参数为各层神经元之间的权值、神经元中的激励函数 及数据流出神经元的阔值。
[0166] 1)各层神经元之间的权值W及数据流出神经元的阔值两个计算参数由random函 数随机生成20种初始的权值矩阵和阔值输入神经网络模型计算，初始的权值和阔值在(0， 1)范围内随机取值(可多次试验找到最为合适的权值和阔值组合）。
[0167] 2)根据供水的实际需要，为了保证每次计算中没有负数的产生影响整个模型计算 的稳定化-擲耐巧掀的佛巧h选择Log-Si卵Oid运一非线性的函数函数：
[016 引
[0169] (5)根据输入历史数据结合相应的网络拓扑结构进行训练，并根据训练结果确定 供水预测模型。
[0170] 模型中训练过程包括采用训练精度模式或学习次数模式的任意一种或多种进行 训练，采用训练精度模式的训练过程包括如下步骤:根据神经网络的网络拓扑结构计算网 络误差，计算训练过程的总误差，判断总误差是否达到设定精度，当判断没有达到精度时重 复上述步骤直至达到精度，根据网络误差修改网络权重并完成训练;采用学习次数模式的 训练过程包括如下步骤:根据网络拓扑结构计算网络误差，根据计算的网络误差修改网络 权重，将训练次数加1，判断是否完成训练次数，若判断没有完成训练次数则重复上述步骤 直至完成训练次数完成训练;训练完成后即可确定供水预测模型。
[0171] 下面对模型中训练过程所设及到的内容进行详细说明：
[0172] 1)正向计算。
[0173] 正向计算所设及的相关理论有：
[0174] a.输入层的输出请(n)等于运个神经网络结构的输入信号x(n)，有：
[0175]
[0176] b.隐含层第i个神经元的输入《;(;!)等于与相应权值《mi(n)的加权和：
[0177]
[017引 C.假设f(.)为Sigmoid函数，则隐含层第i个神经元的输出值)4的)为：
[0179]
[0180] 其中，01为隐含层神经元的阀值。
[0181] d.输出层第j个神经元的输入別对等于诗的加权和：
[0182]
[01831 e .输出房第i个神经元的输出咕(句等于：
[018'
[01化]f.输出层第j个神经元的误差6^11)为：
[0186]
[0187] g.网络的总误差e(n)为：
[018 引
[0189] 2)误差信号反向传播修正。
[0190] 本发明的训练过程采用非线性规划中的最速下降法按误差函数的负梯度方向修 正各个节点的权值及误差，具体包括：
[0191] a.调整隐含层与输出层之间的权值根据最速下降法，需要计算误差对COiJ的 ce{n) 梯度
[0192]
[019；3] 〇ij(n+l)= A 〇ij(n)+?ij(n)[0194] 巧据微分的链式规则有：
[0195]
[0196] 子为一次函数，有：
[0197]
[0199]
[0198] 而输出层传递函数的导数为：
[0200]
[0201]
[0202]
[0203]
[0204]
[0205]
[0206]
[0207]
[020引
[0209]
[0210]
[0211]
[0212]
[0213]
[0214] 3)训练方法和标准
[0215] 本实例中采用训练精度及训练次数共同控制整个训练过程，训练方法使用并行的 训练方式进行训练:神经网络需获得所有的学习样本，计算其所有样本的均方误差和作为 总误差如下式所示：
[0216]
[0217] 为了提高训练精度，本发明采用了 W下两种改进的计算方法计算总误差：
[0218] a.增力喊量项的方法：
[0219] Aw(t)=n5+a Aw(t-l)
[0220] 运时权值修正量A w(t)加上了有关前一时刻权值的修改方向的记忆，使得权值修 正有一定的惯性，具有一定的抗震荡能力，使得迭代过程具有较快的收敛速度。
[0221 ] b.自适应调节学习率的方法：
[0222]
[0223] 若仍采用一般的固定长梯度下降法进行求解，可能会导致两个主要问题:收敛速 度慢和局部极小解。因此，采用改进的变步长算法修正各个节点的权值:先设一个初始步长 no,若一次迭代后误差E增大，则将步长乘W小于1的调节因子0，然后沿原方向重新计算下 一个迭代点；若一次迭代后E减小，则将步长乘W大于1的调节因子0，运样就可W通过误差 的增减来判断、调整学习率，从而提高算法的性能和稳定性。
[0224] 实时预测步骤在预测供水量时，将执行的操作包括：
[0225] (1)特殊节假日处理
[0226] 根据供水量的实际情况，对于春节、国庆和当地的风俗习惯，预测时需要将运些特 殊的节日进行异常处理，W历史当日的数据作为参考来预测供水量。
[0227] (2)自动更新
[0228] 根据建立的供水预测模型在进行预测后，本发明能自动读取数据库中采集到的历 史流量作为新的样本输入到模型中，不断更新预测模型，使模型能随供水条件的变化随时 变化，W适应实际的情况，保证预测的精度。
[0229] 实施例二
[0230] 参照图3-5,本发明的第二实施例：
[0231] 本实施例W贺州市的供水模型为例，贺州城市的供水由担杆岭水厂和河南水厂共 同承担。本发明根据运两个供水厂的实际情况对未来短期的供水量进行预测，来指导自来 水厂进行调度生产。
[0232 ]如图3所示，建模步骤包括：
[0233] 步骤101:分析历史数据并根据供水过程的特性确定供水预测模型的影响因子。
[0234] 本实施例经过相关分析后将高溫度和历史的供水量，节假日情况作为水预测模型 的考虑因子，具体分析过程参考实施例一。
[0235] 步骤102:在确定的影响因子中选择的贡献率高的因子作为参数输入、特殊因子作 特殊处理。
[0236] 本实施例经过分析将高溫度及历史供水量作为输入的神经元输入到BP神经网络 模型中进行计算，节假日作为特殊因子进行特别处理。
[0237] 步骤103:对输入的参数进行预处理；
[0238] 输入参数中，先将流量数据用3倍均方差的方法滤掉异常点，并需要将初次处理过 的流量进行移动平均之后再进行归一化处理，而溫度数据直接进行归一化处理即可
[0239] 步骤104:确定供水预测模型的计算参数，建立输入层参数与输出层输出的供水量 之间的函数关系。
[0240] 该实施例中输入层神经元为2个，输出层神经元为1个，隐含层神经元个数为8个。 由于供水量数据的非负性，本实施例选择Log-Sigmoid函数作为激励函数，能比较准确地反 映出供水的实际过程。
[0241] 步骤105:根据输入的历史数据结合相应的网络拓扑结构进行训练，并根据训练结 果确定供水预测模型。
[0242] 步骤S105中训练的模式可W选用训练精度模式或学习次数模式或二者的结合。
[0243] 如图4所示，其中，采用训练精度模式结合训练次数模式进行训练包括如下步骤：
[0244] 步骤201:根据网络拓扑结构进行正向传播计算；
[0245] 步骤202:计算训练过程的总误差；
[0246] 步骤203:判断总误差是否达到设定精度；
[0247] 步骤204:判断是否达到最大训练次数；
[0248] 步骤205:根据网络误差修改网络权重；
[0249] 当判断没有达到精度时且没有达到最大训练次数时重复上述步骤直至达到精度 并完成训练。
[0250] 步骤206:输入下一组样本进行训练。
[0251 ]步骤S106:根据训练确定的供水预测模型进行供水量预测。
[0252] 如图5所示，步骤106中基于训练模型得到的模型进行预测可进一步分为如下步 骤：
[0253] 步骤301:输入预测所需的参数。
[0254] 步骤302:判断是否与上一次预测日期一致，并根据判断的结果执行步骤。
[0巧5]步骤303:若是，则直接进行正向计算。
[0256] 步骤304:若不是，执行步骤10更新供水预测模型，然后再基于更新后的模型进行 正向计算。
[0257] 步骤305:判断是否为节假日，若是，则执行步骤306,反之，则直接输出预测的结 果。
[0258] 步骤306:根据具体节假日类型进行相应的异常处理，W历史当日的供水数据作为 参考重新进行预测，然后输出重新预测的结果。
[0259] 本实施例利用移动平均法结合改进的BP神经网络的方法对供水量进行短期预测， 通过对过去一年的贺州供水厂的日供水量数据时间序列进行分析和预测，得到了贺州供水 量。通过与传统的预测方法进行对比后发现，采用本发明的方法进行预测有比较明显的周 期性，学习所需数据量少，精度高，能够解决传统方法需要数据量大且计算复杂的问题;相 较于传统单一的移动平均方法而言，本发明在历史流量的数据处理上使用了移动平均的方 法，通过时间序列的分析和试验得到较好的移动平均周期，去掉了部分噪音，提高了样本的 精度;本发明将移动平均的值作为输入因子输入到BP神经网络模型中，由于BP神经网络的 强非线性，城市供水的非线性的特征并未削弱。此外，本发明采用了 BP神经算法，B巧巾经网 络由于其强鲁棒性和容错性，当某些数据发生异常变化时，对整个复杂的神经网络系统的 影响也不大，异常值不会对预测结果产生很大的误差，实时性好;B巧巾经网络由于其强非线 性的特点，整个网络是一个暗箱，能够比较准确反应输入的影响因子和输出的结果之间的 复杂关系，远远优于灰色模型。
[0260] 在传统的神经网络预测的方法上，本发明在输入输出和异常处理的方面进行了改 进，通过平滑处理、优化处理W及移动平均处理等来提高了整个预测模型的精度和稳定性： 考虑到水厂日常运营的可能发生的异常情况(尤其是单一水厂日常供水时，若出现供水设 备故障或管网出现事故，其将会对整个供水量产生较大的影响），故用于分析的样本需要剔 除运些由于异常产生的异常点并且分析其成因。在模型建立和后期更新时，软件程序将自 动除去运些异常的点，因此只需将正确的新样本输入模型，就能保证整个系统的稳定性和 正确性。
[0261] 在理想情况下，历史供水的样本点都是规律性比较强且比较理想的。但实际使用 过程所获取的样本，由于水厂中设备仪表的老化，数据并没有那么理想，故本发明对运些不 符合规律从的点都进行了一一分析，得到了比较准确的样本点，并对于特殊节假日所产生 的影响进行了分析，然后将运些因素作为模型的影响因子进行综合考虑，提高了模型的精 度。
[0262] W上是对本发明的较佳实施进行了具体说明，但本发明创造并不限于所述实施 例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替 换，运些等同的变形或替换均包含在本申请权利要求所限定的范围内。
【主权项】
1. 一种供水预测方法，其特征在于:包括以下步骤： 51、 分析历史数据并建立供水预测模型，得到输入层参数与输出层输出的供水量之间 的函数关系； 52、 根据建立的供水预测模型和所需预测日期进行实时供水量预测； 所述分析历史数据并建立供水预测模型的过程包括但不限于采用移动平均法对输入 的历史数据进行预处理过程，采用神经网络算法确定供水预测模型的网络拓扑结构过程， 根据预处理后的历史数据以及相应的网络拓扑结构确定供水预测模型过程。2. 根据权利要求1所述的一种供水预测方法，其特征在于:所述步骤S1包括： 511、 根据供水过程的特性确定供水预测模型的输入参数； 512、 采用移动平均法对输入的历史数据进行预处理； 513、 采用三层BP神经网络模型建立供水预测模型的网络拓扑结构； 514、 确定三层BP神经网络模型的计算参数； 515、 根据预处理后的历史数据以及三层BP神经网络模型相应的网络拓扑结构进行训 练，然后根据训练的结果确定供水预测模型。3. 根据权利要求2所述的一种供水预测方法，其特征在于:所述步骤S12包括： 5121、 对初始样本中输入的历史数据进行平滑处理，得到去除突变异常点后的样本数 据，所述突变异常点满足：其中，Qin为第η个样本的历史供水量值，n=l，2,…，N; N为最大样本数，σ (Q)为历史供水 量样本的均方差，〇(Q)的计算公式为：5122、 对去除突变异常点后的样本数据中小范围内的波动点和斜率突然变化的样本点 进行优化，得到优化后的样本数据； 5123、 对优化后的样本数据进行移动平均处理，所述移动平均处理所采用的公式为：其中，xq为优化后的样本数据中第q个样本的样本值，g为优化后的样本数据中第q个 样本的移动平均值，h为移动平均的周期数； 5124、 对移动平均处理后的所有样本数据进行归一化处理，所述归一化处理所采用的 公式为：其中，和If分别为移动平均处理后的所有样本中第j个字段归一化处理前的值和 归一化处理后的值，X j max和分别表不所有样本中第j个字段归一化处理前的最大值和 归一化处理后的最大值，Xjmin和i/min:分别表不所有样本中第j个字段归一化处理前的最小 值和归一化处理后的最小值，的取范围为[〇，1]。4. 根据权利要求2所述的一种供水预测方法，其特征在于:所述供水预测模型的输入参 数包括但不限于历史用水量、天气情况、温度、湿度、节假日情况、工业情况、城镇居民生活 条件、人口和GDP中的任意一种或多种，每个参数构成输入层的一个节点，输出层输出第二 天的供水量，并采用相关性分析的方法对输入参数进行分析，计算各参数与供水量的之间 的相关系数及各参数对供水量的贡献率，然后选取与供水量相关系数高的参数作为影响曰 供水量的主要因子，并设定阈值作为神经网络流出神经元的控制量;所述神经网络算法采 用神经元为基本处理单元，所述神经元为多输入、单输出的非线性器件;所述三层BP神经网 络模型的计算参数包括但不限于各层神经元之间的权值、神经元中的激励函数及数据流出 神经元的阈值，所述各层神经元之间的权值以及数据流出神经元的阈值由random函数随机 生成的20种初始的权值矩阵和阈值输入三层BP神经网络模型计算而来，所述初始的权值和 阈值均在（〇，1)之间随机取值，所述神经元中的激励函数采用Log-Sigmoid函数：5. 根据权利要求2所述的一种供水预测方法，其特征在于:所述供水预测模型以三层BP 神经网络作为网络拓扑结构，所述输入层与输出层之间设置有隐含层，所述三层BP神经网 络的输入层数目m、输出层数目m和隐含层数目Π2满足:n2 = sqrt (m+m)+a，其中，sqrt函数 为平方根函数，a为常数，且a取1到10的自然数;所述隐含层的节点根据输入层的节点数确 定并通过增减节点数进行训练对比，所述神经网络算法中每层的输出函数为输入函数关于 该层的激发函数的函数;若所述输出层的输出结果与期望值不符，将输出层的误差转化为 各层各个节点的连接权值与阈值的误差，将输出层的误差逐层向输入层反向传递，把误差 分摊到各个节点，计算各个节点的参考误差，并调整各个节点的连接权值及阈值大小，直到 误差最小。6. 根据权利要求2-5任一项所述的一种供水预测方法，其特征在于:所述步骤S15具体 为： 根据预处理后的历史数据以及三层BP神经网络模型相应的网络拓扑结构采用训练精 度模式、学习次数模式或者两种模式的结合进行训练，然后根据训练的结果确定供水预测 丰旲型； 其中，采用训练精度模式进行训练包括如下步骤:根据三层BP神经网络的网络拓扑结 构计算网络误差，计算训练过程的总误差，判断总误差是否达到设定精度，当判断没有达到 精度时重复上述步骤直至达到设定精度，根据网络误差修改网络权重并完成训练;所述采 用学习次数模式进行训练包括如下步骤:根据三层BP神经网络模型的网络拓扑结构计算网 络误差，根据计算的网络误差修改网络权重，将训练次数加1，判断是否达到最大训练次数， 若判断没有达到最大训练次数则重复上述步骤直至达到最大训练次数完成训练;采用训练 精度模式结合学习次数模式进行训练包括如下步骤:根据三层BP神经网络模型的网络拓扑 结构正向传播计算网络误差，计算训练过程的总误差，判断总误差是否达到设定精度，判断 是否达到最大训练次数，当判断没有达到精度时或没有达到最大训练次数时重复上述步骤 直至达到设定精度和最大训练次数，根据网络误差修改网络权重； 所述训练过程的总误差为所有含有预处理后的历史数据样本误差的均方差之和； 所述训练过程采用改进的变步长梯度下降法修正各个节点的权值及误差，各个节点当 前时刻的权值修正量Δ w(t)为：Δ w(t) =ηδ+α Δ w(t-l)，其中，Δ w(t-l)为各个节点上一时 刻的权值修正量，S为各个节点理论计算得到的权值修正量，α为动量系数，且a e (〇，1)，q为 自适应调节学习率，η的调整方式为:先设定初始步长，若进行一次迭代后节点的误差增大， 则将步长和η均乘以小于1的调节因子β，然后沿原方向重新进行下一次迭代;若进行一次迭 代后节点的误差减小，则将步长和η均乘以大于1的调节因子Θ，然后沿原方向重新进行下一 次迭代。7. 根据权利要求2-5任一项所述的一种供水预测方法，其特征在于:所述步骤S2包括： 521、 输入预测所需的参数。 522、 判断所需预测日期是否与上一次预测日期一致，若是，则根据输入的参数和上一 次预测日期的供水预测模型正向计算所需预测日期的供水量，反之，则在执行步骤S15更新 供水预测模型后正向计算所需预测日期的供水量； 523、 判断所需预测日期是否为节假日，若是，则以历史当日的数据作为参考数据重新 预测所需预测日期的供水量，并将重新预测的供水量进行输出；反之，则直接输出步骤S22 计算出的预测供水量。8. -种供水预测系统，其特征在于:包括： 建模模块，用于分析历史数据并建立供水预测模型，得到输入层参数与输出层输出的 供水量之间的函数关系； 预测模块，用于根据建立的供水预测模型和所需预测日期进行实时供水量预测； 所述分析历史数据并建立供水预测模型的过程包括但不限于采用移动平均法对输入 的历史数据进行预处理过程，采用神经网络算法确定供水预测模型的网络拓扑结构过程， 根据预处理后的历史数据以及相应的网络拓扑结构确定供水预测模型过程； 所述建模模块的输出端与预测模块的输入端连接。9. 根据权利要求8所述的一种供水预测系统，其特征在于:所述建模模块包括： 输入参数确定单元，用于根据供水过程的特性确定供水预测模型的输入参数； 预处理单元，用于采用移动平均法对输入的历史数据进行预处理； 拓扑结构建立单元，用于采用三层ΒΡ神经网络模型建立供水预测模型的网络拓扑结 构； 计算参数确定单元，用于确定三层ΒΡ神经网络模型的计算参数； 训练单元，用于根据预处理后的历史数据以及三层ΒΡ神经网络模型相应的网络拓扑结 构进行训练，然后根据训练的结果确定供水预测模型； 所述输入参数确定单元的输出端依次通过预处理单元、拓扑结构建立单元、计算参数 确定单元和训练单元进而与预测模块的输入端连接。10. 根据权利要求9所述的一种供水预测系统，其特征在于:所述预处理单元包括： 平滑处理子单元，用于对初始样本中输入的历史数据进行平滑处理，得到去除突变异 常点后的样本数据，所述突变异常点满足：其中，Qm为第η个样本的历史供水量值，n=l，2,…，N;N为最大样本数，〇(Q)为历史供水 量样本的均方差，〇(Q)的计算公式为：优化子单元，用于对去除突变异常点后的样本数据中小范围内的波动点和斜率突然变 化的样本点进行优化，得到优化后的样本数据； 移动平均处理子单元，用于对优化后的样本数据进行移动平均处理，所述移动平均处 理所采用的公式为：其中，xq为优化后的样本数据中第q个样本的样本值，€为优化后的样本数据中第q个 Η 样本的移动平均值，h为移动平均的周期数； 归一化处理子单元，用于对移动平均处理后的所有样本数据进行归一化处理，所述归 一化处理所采用的公式为： - -? 其中，和分别为移动平均处理后的所有样本中第j个字段归一化处理前的值和 归一化处理后的值，Xjmax和.1.max_*别表不所有样本中第j个字段归一化处理前的最大值和 归一化处理后的最大值，Xjmin和分别表不所有样本中第.彳个字段归一化处理前的最小 值和归一化处理后的最小值的取范围为[〇，1]; 所述输入参数的输出端确定单元依次通过平滑处理子单元、优化子单元、移动平均处 理子单元和归一化处理子单元进而与拓扑结构建立单元的输入端连接。
【文档编号】G06Q10/04GK105956690SQ201610261368
【公开日】2016年9月21日
【申请日】2016年4月25日
【发明人】代英磊, 篠崎功, 叶景荣
【申请人】广州东芝白云自动化系统有限公司