CN103592962A - 通过预测泵站流量实现污水网管投药量有效控制的方法 - Google Patents
通过预测泵站流量实现污水网管投药量有效控制的方法 Download PDFInfo
- Publication number
- CN103592962A CN103592962A CN201310503431.5A CN201310503431A CN103592962A CN 103592962 A CN103592962 A CN 103592962A CN 201310503431 A CN201310503431 A CN 201310503431A CN 103592962 A CN103592962 A CN 103592962A
- Authority
- CN
- China
- Prior art keywords
- pumping plant
- sigma
- input flow
- parameter
- model
- 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.)
- Granted
Links
Images
Abstract
本发明公开了一种通过预测泵站流量实现污水网管投药量有效控制的方法,首先判断泵站是否装有输入流量检测装置,若有则进行下一步,若没有则需先将泵站的事件式离散的输出流量转换为连续的泵站输入流量;之后利用ARMA模型对泵站的输入流量进行建模和预测,同时利用递推最小二乘法对ARMA模型进行迭代辨识,最终获得泵站输入流量的预测结果,最后根据泵站输入流量的预测结果,及时调整投药,从而实现污水网管内细菌的有效抑制,为污水管网的投药优化提供更多的便利。
Description
技术领域
本发明涉及城市污水处理的技术领域,尤其是指一种通过预测泵站流量实现污水网管投药量有效控制的方法。
背景技术
在城市污水处理系统中,污水管网起到污水收集、运输的作用。由于污水管网覆盖面积广,污水管道需要在地下埋深,为了减少后续输送管道的埋深,在污水输送的中途设置提升泵站,污水经提升后再通过污水管道送到更远的地方(污水处理厂或下一个污水提升泵站),所以污水提升泵站是污水收集中转站。泵站泵的启动和停止基于本泵站进水池水位上下限的限制值,即液位达到上限值开泵而在下限时停泵,而且对于大多数泵站来说,为了减少泵的故障率,对单机运行有时间限制,所以,各个泵站一般会在抽水达到最低水位时候进行水泵相互间的切换。由于泵的开停使得污水在管网长期滞留,从而导致厌氧条件的长期存在,使得存在于生物膜或管网沉淀物中的细菌代谢,特别是SRB细菌代谢会产生了大量的H2S。与此同时,H2S在从液态到气态的转化过程中会产生管网腐蚀、恶臭,这些都会导致安全生产事故。解决该问题的一种主要方式是在污水管网中投加药剂,如氢氧化镁,铁盐,硝酸盐,压硝酸盐等等。而污水在管网的停留时间(HRT,水力停留时间)与泵站出水流量有很大的关系,因此,通过预测流量从而控制投药量的在线方法给投药控制系统的改善提供了可能性。
发明内容
本发明的目的在于克服现有技术的不足,提供一种安全可靠的通过预测泵站流量实现污水网管投药量有效控制的方法。
为实现上述目的,本发明所提供的技术方案为:一种通过预测泵站流量实现污水网管投药量有效控制的方法,首先判断泵站是否装有输入流量检测装置,若有则进行下一步,若没有则需先将泵站的事件式离散的输出流量转换为连续的泵站输入流量;之后利用ARMA模型对泵站的输入流量进行建模和预测,同时利用递推最小二乘法对ARMA模型进行迭代辨识,最终获得泵站输入流量的预测结果,最后根据泵站输入流量的预测结果,及时调整投药,从而实现污水网管内细菌的有效抑制。
所述通过预测泵站流量实现污水网管投药量有效控制的方法,包括以下步骤:
1)判断泵站是否装有输入流量检测装置,若有则直接进行步骤2),若没有则需进行如下步骤后,再进行步骤2):
①泵站输入流量Qin的计算
③利用了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)是观测数据,
其中,na,nc是A(z-1),C(z-1)的阶次;
2)差分阶次辨识
ARMA模型存在着明显的随机序列的特性,由自相关函数ACF和偏相关函数PACF决定,ACF和PACF若能收敛,说明这是一个稳定序列,在ACF和PACF的分析过程中,其收敛行为非常重要,ACF的计算如下:
其中,t=1,2,…,n;j=1,2,…,n,-1;n是数据序列的个数,是y(t)的均值,ACF的分析主要关于ARMA模型AR部分的模型阶次和差分阶次,而另一个可以描述数据序列相关性能的是偏相关系数PACF,PACF计算如下:
Φ11=ρ1
Φj+1,i=Φji-Φj+1,j+1Φj,j+1-i
其中,i=1,2,…,j;Φjj是关于j的函数即是PACF,通常正确的差分阶次能使得随机变量围绕均值附近波动,同时ACF能迅速衰退到0且在95%的区间内,但如果ACF经过长期的衰减仍然无法达到稳定,则说明模型需要更高的差分阶次,其中0阶的差分如下:
其中,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参数进行模型阶次辨识:
其中,N是数据点的个数,n=na+nc是ARMA模型的模型阶次,而 N是总的样本个数;
其中,T是转置,而v(t)则由下式估计:
同时由上述公式估计得到,
其中,L(t)和P(t)分别是估计参数、增益矩阵和协方差矩阵,L(t)和P(t)可视为的过度参数;
4)ARMA模型的多步预测
在辨识阶次和参数的基础上,多步预测变量具体如下式所示:
当1<k≤nc时,
其中, 由噪声发生函数产生得到;
当k>nc时,
5)基于上述ARMA模型泵站入水流量的预测值,确定本网段的投药量:
M=KP·Ypre
其中,M为单位时间内投药量,Kp为比例系数矩阵,Ypre为泵站入水流量预测序列,而比例系数矩阵由历史数据回归得到。
本发明与现有技术相比,具有如下优点与有益效果:
1、利用泵站储液池的工作原理,可以将事件式离散的输出流量转换为连续的输入流量,以便于连续建模方法的处理;
2、利用ARMA模型对连续的输入流量进行建模和预测,在此基础上再次利用泵站储液池的工作原理预测泵站的输入流量,进而为污水管网的投药优化提供更多的便利。
附图说明
图1为基于ARMA模型流量预测管网投药系统图。
图2为两个小型泵站输出流量的预测结果对比图。
具体实施方式
下面结合具体实施例对本发明作进一步说明。
本实施例所述的通过预测泵站流量实现污水网管投药量有效控制的方法,以小型泵站为例,实验数据来自某污水管网中的两个小型泵站,如图1所示,给出了单管投药系统控制图,由于小型泵站一般只在管道输出监控流量,而输入流量没有仪表监控,因此,在实施时需将事件式离散的输出流量转化为连续的输入流量,在此基础上用ARMA模型对连续的流量数据进行建模,最后对不超过6小时的输出流量进行预测。流量的单位均为立方米/天。训练和预测样本各1天数据共135天和96天(包括输出流量和当时的泵站液位,其中液位所得数据是按体积的百分比获得),储水池的大小分别为100立方米和130立方米,具体的实现步骤如下:
1)判断泵站是否装有输入流量检测装置,若有则直接进行步骤2),若没有则需进行如下步骤后,再进行步骤2):
①泵站输入流量Qin的计算
其中,Qout已知,采样时间为0.0104天,V则对应相应的泵站体积;
②利用了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)是观测数据,
其中,na,nc是A(z-1),C(z-1)的阶次。
2)差分阶次辨识
ARMA模型存在着明显的随机序列的特性,由自相关函数ACF和偏相关函数PACF决定,ACF和PACF若能收敛,说明这是一个稳定序列,在ACF和PACF的分析过程中,其收敛行为非常重要,ACF的计算如下:
其中,t=1,2,…,n;j=1,2,…,n,-1;n是数据序列的个数(分别为135×96和96×96个数据点),是y(t)的均值,ACF的分析主要关于ARMA模型AR部分的模型阶次和差分阶次,而另一个可以描述数据序列相关性能的是偏相关系数PACF,PACF计算如下:
Φ11=ρ1
Φj+1,i=Φji-Φj+1,j+1Φj,j+1-i
其中,i=1,2,…,j;Φjj是关于j的函数即是PACF,通常正确的差分阶次能使得随机变量围绕均值附近波动,同时ACF能迅速衰退到0且在95%的区间内,但如果ACF经过长期的衰减仍然无法达到稳定,则说明模型需要更高的差分阶次,其中0阶的差分如下:
其中,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),由于上述i=19和j=18的时候ACF或PACF基本能处于收敛状态,因此具体的0阶差分如表1所示:
表1:泵站1和泵站2的ACF和PACF参数
从上表可知对于泵站1,PACF有明显衰减而ACF无,说明模型倾向于AR模型,而且PACF衰减到7的时候已经达到95%的限制,故na=7.nc=0。而对于泵站2,类似的PACF有明显衰减而ACF无,说明模型倾向于AR模型,而且PACF衰减到7的时候已经达到95%的限制,故na=3.nc=0。通过类似的方法计算1和2阶的有明显过差分的ACF和PACF,即ACF和PACF从1或2开始基本都是负值。
3)差分模型参数和阶次辨识
①模型的阶次可通过统计学F分布测试和贝叶斯信息原则BIC参数离线方式获得,模型的阶次过高会增加模型的计算量,反之会降低模型的预测精度,因此,需要选择合适的模型参数,利用BIC参数进行模型阶次辨识:
其中,N是数据点的个数,n=na+nc是ARMA模型的模型阶次,而 N是总的样本个数;
其中,T是转置,而v(t)则由下式估计:
由前述辨识可知对于泵站1和泵站2,其na分别为7和3,而nc均为0,再根据上述步骤①和②对历史数据对模型进行辨识,最后得出参数θ,经过辨识,两个泵站的参数分别如表2所示:
表2:泵站1和泵站2的辨识参数
4)ARMA模型的多步预测
在辨识阶次和参数的基础上,多步预测变量具体如下式所示:
当1<k≤nc时,
其中, 由噪声发生函数产生得到;
当k>nc时,
而在本实施例中,nc=0,k>nc
5)基于上述ARMA模型泵站入水流量的预测值,确定本网段的投药量:
M=KP·Ypre
其中,M为单位时间内投药量,Kp为比例系数矩阵,这里对于泵站1和2分别为0,07和0.08,Ypre为泵站入水流量预测序列,而具体的比例系数矩阵由历史数据回归得到。
在采用以上方案后,本发明能有效确保污水管网投药量的准确控制,充分考虑了泵站工作方式不同天气状况下的水流模式,可以对泵站的输出流量做长达4~6小时的预测,从而为污水管网的投药优化提供更多的便利,值得推广。
以上所述之实施例子只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。
Claims (2)
1.一种通过预测泵站流量实现污水网管投药量有效控制的方法,其特征在于:首先判断泵站是否装有输入流量检测装置,若有则进行下一步,若没有则需先将泵站的事件式离散的输出流量转换为连续的泵站输入流量;之后利用ARMA模型对泵站的输入流量进行建模和预测,同时利用递推最小二乘法对ARMA模型进行迭代辨识,最终获得泵站输入流量的预测结果,最后根据泵站输入流量的预测结果,及时调整投药,从而实现污水网管内细菌的有效抑制。
2.根据权利要求1所述的一种通过预测泵站流量实现污水网管投药量有效控制的方法,其特征在于,包括以下步骤:
1)判断泵站是否装有输入流量检测装置,若有则直接进行步骤2),若没有则需进行如下步骤后,再进行步骤2):
①泵站输入流量Qin的计算
其中,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)是观测数据,
其中,na,nc是A(z-1),C(z-1)的阶次;
2)差分阶次辨识
ARMA模型存在着明显的随机序列的特性,由自相关函数ACF和偏相关函数PACF决定,ACF和PACF若能收敛,说明这是一个稳定序列,在ACF和PACF的分析过程中,其收敛行为非常重要,ACF的计算如下:
其中,t=1,2,…,n;j=1,2,…,n,-1;n是数据序列的个数,是y(t)的均值,ACF的分析主要关于ARMA模型AR部分的模型阶次和差分阶次,而另一个可以描述数据序列相关性能的是偏相关系数PACF,PACF计算如下:
Φ11=ρ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)=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参数进行模型阶次辨识:
其中,N是数据点的个数,n=na+nc是ARMA模型的模型阶次,而 N是总的样本个数;
其中,T是转置,而v(t)则由下式估计:
4)ARMA模型的多步预测
在辨识阶次和参数的基础上,多步预测变量具体如下式所示:
当1<k≤nc时,
其中, 由噪声发生函数产生得到;
当k>nc时,
5)基于上述ARMA模型泵站入水流量的预测值,确定本网段的投药量:
M=KP·Ypre
其中,M为单位时间内投药量,Kp为比例系数矩阵,Ypre为泵站入水流量预测序列,而比例系数矩阵由历史数据回归得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310503431.5A CN103592962B (zh) | 2013-10-23 | 2013-10-23 | 通过预测泵站流量实现污水网管投药量有效控制的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310503431.5A CN103592962B (zh) | 2013-10-23 | 2013-10-23 | 通过预测泵站流量实现污水网管投药量有效控制的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103592962A true CN103592962A (zh) | 2014-02-19 |
CN103592962B CN103592962B (zh) | 2016-05-04 |
Family
ID=50083153
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310503431.5A Expired - Fee Related CN103592962B (zh) | 2013-10-23 | 2013-10-23 | 通过预测泵站流量实现污水网管投药量有效控制的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103592962B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106595788A (zh) * | 2016-12-20 | 2017-04-26 | 河海大学 | 基于多参数耦合修正的大型泵站流量监测方法 |
CN109313416A (zh) * | 2016-06-14 | 2019-02-05 | 苏伊士集团 | 对流体网络进行动态控制的方法、计算机程序产品和系统 |
CN110110890A (zh) * | 2019-03-28 | 2019-08-09 | 杭州电子科技大学 | 基于elman神经网络的日污水量预测方法 |
CN113568440A (zh) * | 2021-07-21 | 2021-10-29 | 上海市城市建设设计研究总院(集团)有限公司 | 基于多数据源分析的大型污水系统辅助调度方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004250939A (ja) * | 2003-02-19 | 2004-09-09 | Toshiba Corp | 雨水排水ポンプ運転支援装置 |
CN101424943A (zh) * | 2008-12-18 | 2009-05-06 | 浙江大学 | 基于互联网的污水流量跟踪协调优化控制方法及其系统 |
CN102013050A (zh) * | 2010-10-22 | 2011-04-13 | 上海海事大学 | 一种污水流量变化模式辨识的mcmc方法 |
CN102032935A (zh) * | 2010-12-07 | 2011-04-27 | 杭州电子科技大学 | 城市排水汇流管网污水泵站流量的软测量方法 |
-
2013
- 2013-10-23 CN CN201310503431.5A patent/CN103592962B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004250939A (ja) * | 2003-02-19 | 2004-09-09 | Toshiba Corp | 雨水排水ポンプ運転支援装置 |
CN101424943A (zh) * | 2008-12-18 | 2009-05-06 | 浙江大学 | 基于互联网的污水流量跟踪协调优化控制方法及其系统 |
CN102013050A (zh) * | 2010-10-22 | 2011-04-13 | 上海海事大学 | 一种污水流量变化模式辨识的mcmc方法 |
CN102032935A (zh) * | 2010-12-07 | 2011-04-27 | 杭州电子科技大学 | 城市排水汇流管网污水泵站流量的软测量方法 |
Non-Patent Citations (3)
Title |
---|
RAMON GANIGUE等: "On-line control of magnesium hydroxide dosing for sulfide mitigation in sewers", 《AWA OZWATER 2012 CONF.》, 31 December 2012 (2012-12-31) * |
王凤仙等: "城市污水量预测模型及方法综述", 《河南科学》, vol. 27, no. 4, 30 April 2009 (2009-04-30), pages 483 - 487 * |
黄大荣等: "基于ARMA和小波变换的交通流预测模型研究", 《计算机工程与应用》, vol. 42, no. 36, 15 January 2007 (2007-01-15), pages 191 - 194 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109313416A (zh) * | 2016-06-14 | 2019-02-05 | 苏伊士集团 | 对流体网络进行动态控制的方法、计算机程序产品和系统 |
CN106595788A (zh) * | 2016-12-20 | 2017-04-26 | 河海大学 | 基于多参数耦合修正的大型泵站流量监测方法 |
CN106595788B (zh) * | 2016-12-20 | 2018-12-11 | 河海大学 | 基于多参数耦合修正的大型泵站流量监测方法 |
CN110110890A (zh) * | 2019-03-28 | 2019-08-09 | 杭州电子科技大学 | 基于elman神经网络的日污水量预测方法 |
CN113568440A (zh) * | 2021-07-21 | 2021-10-29 | 上海市城市建设设计研究总院(集团)有限公司 | 基于多数据源分析的大型污水系统辅助调度方法 |
CN113568440B (zh) * | 2021-07-21 | 2023-09-01 | 上海市城市建设设计研究总院(集团)有限公司 | 基于多数据源分析的大型污水系统辅助调度方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103592962B (zh) | 2016-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103197544B (zh) | 基于非线性模型预测的污水处理过程多目标控制方法 | |
Foley et al. | Dissolved methane in rising main sewer systems: field measurements and simple model development for estimating greenhouse gas emissions | |
CN103886369B (zh) | 一种基于模糊神经网络的出水总磷tp预测方法 | |
CN107368623B (zh) | 一种基于径流污染控制的调蓄池容积确定方法 | |
CN103592962B (zh) | 通过预测泵站流量实现污水网管投药量有效控制的方法 | |
Ganigué et al. | Improved sulfide mitigation in sewers through on-line control of ferrous salt dosing | |
CN101786721B (zh) | 城市污水处理厂出水有机物浓度的随机过程预测方法 | |
CN110054274A (zh) | 一种净水絮凝沉淀投药控制技术 | |
CN105574326A (zh) | 一种基于自组织模糊神经网络的出水氨氮浓度软测量方法 | |
Chen et al. | Real-time multistep prediction of sewer flow for online chemical dosing control | |
Olsson et al. | Operation and control of wastewater treatment—a Scandinavian perspective over 20 years | |
CN102831324A (zh) | 一种基于区间数理论河流环境容量的计算方法 | |
CN104090974A (zh) | 展延水库后续来水的动态数据挖掘方法及系统 | |
CN106200381A (zh) | 一种根据处理水量分阶段控制水厂运行的方法 | |
CN113608444A (zh) | 基于自适应预测控制的污水处理控制方法 | |
Sun et al. | Advanced integrated real-time control of combined Urban drainage systems using MPC | |
CN113077847A (zh) | 在线模型水质转换方法、系统、电子设备及介质 | |
CA3035131A1 (en) | Treatment of wastewater | |
CN116681174A (zh) | 污水处理生化池曝气量预测方法、系统、设备及介质 | |
CN103942600B (zh) | 基于尖峰自组织径向基神经网络的污泥膨胀预测方法 | |
CN113642256B (zh) | 一种基于主题数据库的实用闸门启闭优化方法 | |
Zaher et al. | Application of the IWA ADM1 model to simulate anaerobic digester dynamics using a concise set of practical measurements | |
Meng et al. | Integrated modelling and control of urban wastewater systems | |
De Haas | The energy versus nitrous oxide emissions nexus in wastewater treatment systems | |
Stentoft et al. | Smart, Data-Driven Aeration in Water Resource Recovery Facilities for Balancing Wastewater Treatment and Energy Consumption |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160504 Termination date: 20211023 |