通过预测泵站流量实现污水网管投药量有效控制的方法

　　申请日2013.10.23

　　公开(公告)日2014.02.19

　　IPC分类号G05B13/04; G05D7/06

　　摘要

　　本发明公开了一种通过预测泵站流量实现污水网管投药量有效控制的方法，首先判断泵站是否装有输入流量检测装置，若有则进行下一步，若没有则需先将泵站的事件式离散的输出流量转换为连续的泵站输入流量;之后利用ARMA模型对泵站的输入流量进行建模和预测，同时利用递推最小二乘法对ARMA模型进行迭代辨识，最终获得泵站输入流量的预测结果，最后根据泵站输入流量的预测结果，及时调整投药，从而实现污水网管内细菌的有效抑制，为污水管网的投药优化提供更多的便利。

　　权利要求书

　　1.一种通过预测泵站流量实现污水网管投药量有效控制的方法，其特征在于：首先判断泵站是否装有输入流量检测装置，若有则进行下一步，若没有则需先将泵站的事件式离散的输出流量转换为连续的泵站输入流量;之后利用 ARMA模型对泵站的输入流量进行建模和预测，同时利用递推最小二乘法对ARMA 模型进行迭代辨识，最终获得泵站输入流量的预测结果，最后根据泵站输入流量的预测结果，及时调整投药，从而实现污水网管内细菌的有效抑制。

　　2.根据权利要求1所述的一种通过预测泵站流量实现污水网管投药量有效控制的方法，其特征在于，包括以下步骤：

　　1)判断泵站是否装有输入流量检测装置，若有则直接进行步骤2)，若没有则需进行如下步骤后，再进行步骤2)：

　　①泵站输入流量Qin的计算

　　 $Q_{in} = Q_{out} + \frac{dV}{dt}$

　　其中，Qin是输入流量，Qout是输出流量，是泵站内污水体积的变化;

　　②利用了3σ的方法对离群点进行检测

　　mean|Qin|≤3σ

　　此外，受天气影响泵站的输入流量主要有晴天、中雨和暴雨三种工作模式，离群点的检测只需将各个模式下的均值按上述公式进行分类计算;

　　③泵站的输入流量是一个随机变量所代表的自回归过程,在同一样本区间内的一个变量可以作为它们过去值的线性函数

　　A(Z-1)y(t)=C(z-1)v(t)

　　其中，Z是滞后因子,z-1y(t)=y(t-1),v(t)是0均值的高斯白噪声，y(t)是观测数据,

　　 $A (Z^{- 1}) = 1 + a_{1} z^{- 1} + . . . + a_{n_{a}} z^{- n_{a}}$

　　 $C (z^{- 1}) = 1 + c_{1} z^{- 1} + . . . + c_{n_{c}} z^{- n_{c}}$

　　其中，na，nc是A(z-1)，C(z-1)的阶次;

　　2)差分阶次辨识

　　ARMA模型存在着明显的随机序列的特性，由自相关函数ACF和偏相关函数 PACF决定，ACF和PACF若能收敛，说明这是一个稳定序列，在ACF和PACF的分析过程中，其收敛行为非常重要，ACF的计算如下：

　　 $ρ_{j} = \frac{Σ_{t = 1}^{n - j} (y (t) - \overset{&OverBar;}{y}) (y (t + j) - \overset{&OverBar;}{y})}{Σ_{t = 1}^{n} {(y (t) - \overset{&OverBar;}{y})}^{2}}$

　　其中,t=1,2,…,n;j=1,2,…，n,-1;n是数据序列的个数，是y(t)的均值，ACF的分析主要关于ARMA模型AR部分的模型阶次和差分阶次，而另一个可以描述数据序列相关性能的是偏相关系数PACF，PACF计算如下：

　　Φ11=ρ1

　　 $Φ_{j + 1, j + 1} = (ρ_{j} + 1 - Σ_{i = 1}^{j} ρ_{j} + 1 - i Φ_{j, i}) {(1 - Σ_{i = 1}^{j} ρ_{j} Φ_{ji})}^{- 1}$

　　Φj+1，i=Φji-Φj+1，j+Φj，j+1-i

　　其中，i=1,2,…,j;Φjj是关于j的函数即是PACF，通常正确的差分阶次能使得随机变量围绕均值附近波动，同时ACF能迅速衰退到0且在95%的区间内，但如果ACF经过长期的衰减仍然无法达到稳定，则说明模型需要更高的差分阶次，其中0阶的差分如下：

　　 $X (t) = a_{1} X (t - 1) + . . . + a_{n_{a}} X (t - n_{a}) + c_{1} v (t - 1) + . . . + c_{n_{c}} v (t - n_{c}) + v (t)$

　　其中，X(t)=y(t)-μ是0阶差分，即y(t)=X(t)+μ，而对于1和2阶来说，分别等于y(t-1)-y(t-2)和y(t-1)-2y(t-2)+y(t-3);

　　3)差分模型参数和阶次辨识

　　①模型的阶次可通过统计学F分布测试和贝叶斯信息原则BIC参数离线方式获得，模型的阶次过高会增加模型的计算量，反之会降低模型的预测精度，因此，需要选择合适的模型参数，利用BIC参数进行模型阶次辨识：

　　 $BIS (n) = N \ln σ_{a}^{2} + n \ln N$

　　其中，N是数据点的个数，n=na+nc是ARMA模型的模型阶次，而 $σ_{a}^{2} = \frac{1}{N - n} Σ_{t = n + 1}^{N} {(y_{t} - Σ_{i = 1}^{na} a_{i} y_{t - i} - Σ_{i = 1}^{nc} b_{i} v_{t - i})}^{2},$ N是总的样本个数;

　　②在每次BIC的确认过程中，参数ai，ci由递推最小二乘法RELS辨识获得，而参数矩阵θ和则表示如下：

　　 $θ = {[a_{1}, a_{2}, . . ., a_{n_{a}}, c_{1}, c_{2,, . . .,} c_{n_{c}}]}^{T}$

　　其中，T是转置，而v(t)则由下式估计：

　　同时由上述公式估计得到，

　　 $\hat{θ} (t) = \hat{θ} (t - 1) + L (t)$

　　其中，L(t)和P(t)分别是估计参数、增益矩阵和协方差矩阵，L(t)和P(t) 可视为的过度参数;

　　4)ARMA模型的多步预测

　　在辨识阶次和参数的基础上，多步预测变量具体如下式所示：

　　当1

　　 $\hat{y} (t + k | t) = Σ_{i = 1}^{n_{a}} a_{i} \hat{y} (t + k - i | t) + Σ_{i = k}^{n_{c}} c_{i} \hat{v} (t + k - i)$

　　其中， $\hat{v} (t), \hat{v} (t - 1), . . ., \hat{v} (t + k - n_{c})$ 由噪声发生函数产生得到;

　　当k>nc时，

　　 $\hat{y} (t + k | t) = Σ_{i = 1}^{n_{a}} a_{i} \hat{y} (t + k - i | t)$

　　当t+k-i≤t时，假定L>nc，多步最优预测 k-1，2，...L可以在所辨识的参数ai，ci的基础上获得;

　　5)基于上述ARMA模型泵站入水流量的预测值，确定本网段的投药量：

　　M=KP·Ypre

　　其中，M为单位时间内投药量，Kp为比例系数矩阵，Ypre为泵站入水流量预测序列，而比例系数矩阵由历史数据回归得到。

　　说明书

　　通过预测泵站流量实现污水网管投药量有效控制的方法

　　技术领域

　　本发明涉及城市污水处理的技术领域，尤其是指一种通过预测泵站流量实现污水网管投药量有效控制的方法。

　　背景技术

　　在城市污水处理系统中，污水管网起到污水收集、运输的作用。由于污水管网覆盖面积广，污水管道需要在地下埋深，为了减少后续输送管道的埋深，在污水输送的中途设置提升泵站，污水经提升后再通过污水管道送到更远的地方(污水处理厂或下一个污水提升泵站)，所以污水提升泵站是污水收集中转站。泵站泵的启动和停止基于本泵站进水池水位上下限的限制值，即液位达到上限值开泵而在下限时停泵，而且对于大多数泵站来说，为了减少泵的故障率，对单机运行有时间限制，所以，各个泵站一般会在抽水达到最低水位时候进行水泵相互间的切换。由于泵的开停使得污水在管网长期滞留，从而导致厌氧条件的长期存在，使得存在于生物膜或管网沉淀物中的细菌代谢，特别是SRB细菌代谢会产生了大量的H2S。与此同时，H2S在从液态到气态的转化过程中会产生管网腐蚀、恶臭，这些都会导致安全生产事故。解决该问题的一种主要方式是在污水管网中投加药剂，如氢氧化镁，铁盐，硝酸盐，压硝酸盐等等。而污水在管网的停留时间(HRT，水力停留时间)与泵站出水流量有很大的关系，因此，通过预测流量从而控制投药量的在线方法给投药控制系统的改善提供了可能性。

　　发明内容

　　本发明的目的在于克服现有技术的不足，提供一种安全可靠的通过预测泵站流量实现污水网管投药量有效控制的方法。

　　为实现上述目的，本发明所提供的技术方案为：一种通过预测泵站流量实现污水网管投药量有效控制的方法，首先判断泵站是否装有输入流量检测装置，若有则进行下一步，若没有则需先将泵站的事件式离散的输出流量转换为连续的泵站输入流量;之后利用ARMA模型对泵站的输入流量进行建模和预测，同时利用递推最小二乘法对ARMA模型进行迭代辨识，最终获得泵站输入流量的预测结果，最后根据泵站输入流量的预测结果，及时调整投药，从而实现污水网管内细菌的有效抑制。

　　所述通过预测泵站流量实现污水网管投药量有效控制的方法，包括以下步骤：