2024年6月11日发(作者:梁丘古韵)
第2O卷第4期
2012年8月
安徽建筑_T-业学院学报(自然科学版)
Journa1 of Anhui Institute of Architecture&Industry
V0I.2O No.4
Aug.2012
基于ARMA模型的自生动态膜生物反应器
中试装置出水COD预测
黄小红
(安徽沃特星水处理运营有限公司,合肥230031)
摘要:根据污水厂大量的在线监测数据对进水COD建立了ARMA预测模型,并结合自生动态膜生物反应
器中试装置,利用莫诺方程求解出适应该装置的进出水COD对应关系,最终达到对装置出水COD进行短期
预测的效果。
关键词:ARMA模型;莫诺方程;COD预测
中图分类号:X328 文献标识码:A 文章编号:1006—4540(2012 J04—075—03
Prediction of effluent COD of Self-Forming Dynamic Membrane Bio-Reactor
pilot-scale equipment Based on the ARMA model
HUANG Xiao—hong
(Anhui Waterstar water Disposal Function Co.,Ltd.,Hefei 230031,China)
Abstract:ARMA model was build to predict influent COD according to a large number of wastewater
treatment plant On-line monitoring data,combined with Self-Forming Dynamic Membrane Bio—Reac—
tor pilot—scale equipment,to solve corresponding relation between influent COD and effluent COD U—
sing Monod equation,finally achieve the short—term prediction of equipment effluent COD.
Key words:ARMA model;Monod equation;COD prediction
l 引 言
ARMA模型对污水厂进水COD进行预测,利用
莫诺方程建立自生动态膜生物反应器进出水
COD对应关系,便于对出水COD的短期预测。
污水经生化处理后,出水有机物排入水体中,
在有溶解氧的条件下,会同时存在微生物呼吸作
用的耗氧与大气向水体复氧的过程。若排人的有
机物量超过水体的环境容量,则好氧速度会超过
复氧速度,水体出现缺氧甚至无氧;在水体缺氧的
2 ARMA模型
ARMA方法是一种精确度较高的短期时间
序列预测方法。ARMA模型[1]的一般形式为:
条件下由于厌氧微生物的作用,有机物被降解为
CH ,CO ,NH。及少量HzS等有害有臭气体,使
Xt=∑ x 一∑ 口 +口 式中, 和 分别
i一1 J=1
水质恶化“黑臭”。为避免水质恶化情况的发生,
对处理后排放水中的化学需氧量进行预测,以便
于及时调整运行参数便显得尤为重要。本文采用
收稿日期:2012—04-28
表示AR部分和MA部分的阶数, 、 分别为各
部分的模型参数,{a )为白噪声序列,可见t时刻
的观测值可以用 时刻之前的观测值x 及 时
作者简介:黄小红(1985一),女,主要研究方向城镇污水处理理论与技术。
76 安徽建筑s-,_lk学院学报(自然科学版) 第2O卷
刻之前进入系统的随机干扰 来逼近,在t时刻
a 是不可观测的随机性干扰,X的随机性正是由
于乜 的随机性所导致的。运用ARMA模型的建
模及预测步骤:
①准备阶段:在得到一组样本数据后,应首先
检验该数据的平稳性及零均值性,若不满足上述
性质,则需作平稳化及零均值化处理;②模式识
别:根据已知的样本数列,判断该序列属于何种模
型,其阶数是多少;③参数估计:在已识别的模型
及其阶数的基础上,对模型参数进行估计,求出初
始模型;④模型检验:用统计检验的方法对初始模
型的合理性进行检验;⑤进行预测:若模型通过检
验,即可用此模型对序列的未来值进行预测。
3 污水厂进水COD预测
现拟用该方法对某污水厂进水COD进行短
期预测,并结合在该污水厂搭建的中试反应器装
置,利用莫诺方程求解出通过反应器后的出水
C0D浓度,最终达到通过对进水C0D建立的预
测模型,进而转化为对出水COD的短期预测的
效果,便于根据对出水值的预测及时调整反应器
的相应操作条件。
考虑到中试现场的取水样时间,故而选取污
水厂2012年4月1日至2011年7月31日每天
上午9时的在线监测进水COD数据,以其中4月
l号至7月24号的数据为原始数据,用于建模。
将收集的数据运用非参数检验法_2]检验其平
稳性,用的统计量为:
一
进猩 =鎏堡 游程标准差 望 一 二 z (]、
式中 一 +1研一LFZN,N 2(2
NIN ̄
-
N)
] ,
(N
—
1)
-J
N—N1+N2,N1为“+”出现的次数,为59;N2为
“
一
”出现的次数,为56;),为游程的总数,为58。经
计算,对于a:==0.05的显著性水平,l Z J一0.086
<1.96a :==10.46,认为该序列为平稳的时间序
列,并对其作零均值化处理,进而计算其自相关
(ACF)及偏相关(PACF)函数,会发现自相关及
偏相关函数均呈现拖尾性,故其适合ARMA
模型。
对平稳系统来说,可以用ARMA(n,n—1)模
型为一般形式建立模型l_3],通过对比最终选择的
预测模型为:
X一一0.354X,一1+0.923X 一2+0.386X 一3
+以 +口 一l (2)
此外,通过计算Ljung--Box统计量Q—N(N+
坚Ⅳ) ^
2)∑(N一是)1pl,式中N为样本数,志为滞后,
一
1
A
L(N)取18, 为{a }的估计值。得Q值为
17.857,小于矮1(14)一21.064及 。5(14)=
23.685,均表明该模型适用。
利用该模型对污水厂进水COD的拟合效果
如图1所示,拟合值的最小相对误差为0.26 ,
平均相对误差为2O.34 ,从图中两条曲线高度
接近可见利用所选模型对原有数据的拟合效果
较好。
犬毅
图1 9.1—9.24号观测值与拟合值对比图
用该模型对污水厂12.25—12.31号的进水
COD进行预测,其效果如图2所示。这里采用平
均相对误差和均等系数评价预测效果的优劣。预
测值的平均相对误差为12.94 ,均等系数为
0.9230,预测结果可视为满意。
l 80
160
140
—
】2o
县100
。80
0
u 60
40
20
0
天数
图2 12.25—12.31号观测值与预测值对比图
4反应器出水COD预测
我们在该污水厂搭建了一座自生动态膜生物
反应器中试装置,此反应器内的流态为完全混合
流[ ,中试系统中的有机底物进行物料平衡,
则有:
』0
S。Q—Qs + 一0 (3)
C
式中S。为进水有机物浓度;Q为处理水量;S 为
第4期 黄小红:基于ARMA模型的自生动态膜生物反应器中试装置出水CIOD预测
一
77
出水有机物浓度;V为反应器有效容积。
莫诺于1942年用纯种的微生物在单一底物
的培养基上进行了微生物增殖速率与底物浓度之
间关系的试验,得到莫诺方程 ]式:
Q 1+Q 3 +Q国 卜3+ +aid
~
——T
(式中以 S0代替(1)式中的x) (7)
即在一定的污泥浓度条件下,水力停留时间
基本恒定时,可根据由ARMA模型对进水COD
一 一x雳一VVma 丽XS丽
一
(4)
的预测,求得经反应器处理后的出水COD的预
式中 为有机底物降解速度;‰ 为有机底物的
最大比降解速度,t一;X为混合液中活性污泥总
量;S为有机底物浓度;K 为饱和常数。
自生动态膜生物反应器中,滤膜的表面会生
成滤饼生物膜,滤饼生物膜对出水COD有一定
的截阻能力,但这种截阻能达到假稳态(或稳定
化),故滤饼生物膜对出水COD的影响忽略不
计。最终可得如下关系式:
Q(So~S ) S 0-Se一 (5)
式中t为水力停留时间。将(4)式取倒数,按直线
方程y:aX+6考虑,利用图解法求解常数值Kz
Vm
==—
ax
—
K 。
利用11月份10天的实测反应装置进出水数
据求解。反应器有效容积为2.5m。。经计算得,
K2—8.13×10_。,则由(4)式得
S 一
1+KzXt
S0
So
1+8.13×10。Xt (6)
式中S。为进水COD浓度,(mg/L);X为污泥浓
度,(mg/L);t为水力停留时间,(d)。
由(5)式并结合进水COD预测模型,可得如
下关系式:
测值。
5结束语
对污水厂进水COD建立了ARMA预测模
型,取得了满意的预测效果,并将中试反应器出水
与污水厂进水COD之间建立了对应关系,可实
现对反应器出水COD的短期预测。有利于根据
预测值是否达标,判断反应器的工作状况,便于对
反应器的操作条件进行及时的、必要的调整,如调
整曝气量的大小、污泥负荷等来保证出水水质。
参考文献
1段智彬,孙恩昌,张延华,董燕基于ARMA模型的网
络流量预测EJ;.中国电子科学研究院学报,2009,4(4):
352—354.
2张树京,齐立心.时间序列分析简明教程EMJ.北京:清
华大学出版社,2003.
3王振龙,胡永宏.应用时间序列分析[M].北京:科学出
版社,2007.
4夏俊林.自生动态膜生物反应器处理城市污水的动态
研究[D].北京:北京科技大学硕士论文,2005.
5张自杰,林荣忱,金儒霖.排水工程下册(第4版)
[M].北京:中国建筑工业出版社,2006.
(责任编辑:王颖)
2024年6月11日发(作者:梁丘古韵)
第2O卷第4期
2012年8月
安徽建筑_T-业学院学报(自然科学版)
Journa1 of Anhui Institute of Architecture&Industry
V0I.2O No.4
Aug.2012
基于ARMA模型的自生动态膜生物反应器
中试装置出水COD预测
黄小红
(安徽沃特星水处理运营有限公司,合肥230031)
摘要:根据污水厂大量的在线监测数据对进水COD建立了ARMA预测模型,并结合自生动态膜生物反应
器中试装置,利用莫诺方程求解出适应该装置的进出水COD对应关系,最终达到对装置出水COD进行短期
预测的效果。
关键词:ARMA模型;莫诺方程;COD预测
中图分类号:X328 文献标识码:A 文章编号:1006—4540(2012 J04—075—03
Prediction of effluent COD of Self-Forming Dynamic Membrane Bio-Reactor
pilot-scale equipment Based on the ARMA model
HUANG Xiao—hong
(Anhui Waterstar water Disposal Function Co.,Ltd.,Hefei 230031,China)
Abstract:ARMA model was build to predict influent COD according to a large number of wastewater
treatment plant On-line monitoring data,combined with Self-Forming Dynamic Membrane Bio—Reac—
tor pilot—scale equipment,to solve corresponding relation between influent COD and effluent COD U—
sing Monod equation,finally achieve the short—term prediction of equipment effluent COD.
Key words:ARMA model;Monod equation;COD prediction
l 引 言
ARMA模型对污水厂进水COD进行预测,利用
莫诺方程建立自生动态膜生物反应器进出水
COD对应关系,便于对出水COD的短期预测。
污水经生化处理后,出水有机物排入水体中,
在有溶解氧的条件下,会同时存在微生物呼吸作
用的耗氧与大气向水体复氧的过程。若排人的有
机物量超过水体的环境容量,则好氧速度会超过
复氧速度,水体出现缺氧甚至无氧;在水体缺氧的
2 ARMA模型
ARMA方法是一种精确度较高的短期时间
序列预测方法。ARMA模型[1]的一般形式为:
条件下由于厌氧微生物的作用,有机物被降解为
CH ,CO ,NH。及少量HzS等有害有臭气体,使
Xt=∑ x 一∑ 口 +口 式中, 和 分别
i一1 J=1
水质恶化“黑臭”。为避免水质恶化情况的发生,
对处理后排放水中的化学需氧量进行预测,以便
于及时调整运行参数便显得尤为重要。本文采用
收稿日期:2012—04-28
表示AR部分和MA部分的阶数, 、 分别为各
部分的模型参数,{a )为白噪声序列,可见t时刻
的观测值可以用 时刻之前的观测值x 及 时
作者简介:黄小红(1985一),女,主要研究方向城镇污水处理理论与技术。
76 安徽建筑s-,_lk学院学报(自然科学版) 第2O卷
刻之前进入系统的随机干扰 来逼近,在t时刻
a 是不可观测的随机性干扰,X的随机性正是由
于乜 的随机性所导致的。运用ARMA模型的建
模及预测步骤:
①准备阶段:在得到一组样本数据后,应首先
检验该数据的平稳性及零均值性,若不满足上述
性质,则需作平稳化及零均值化处理;②模式识
别:根据已知的样本数列,判断该序列属于何种模
型,其阶数是多少;③参数估计:在已识别的模型
及其阶数的基础上,对模型参数进行估计,求出初
始模型;④模型检验:用统计检验的方法对初始模
型的合理性进行检验;⑤进行预测:若模型通过检
验,即可用此模型对序列的未来值进行预测。
3 污水厂进水COD预测
现拟用该方法对某污水厂进水COD进行短
期预测,并结合在该污水厂搭建的中试反应器装
置,利用莫诺方程求解出通过反应器后的出水
C0D浓度,最终达到通过对进水C0D建立的预
测模型,进而转化为对出水COD的短期预测的
效果,便于根据对出水值的预测及时调整反应器
的相应操作条件。
考虑到中试现场的取水样时间,故而选取污
水厂2012年4月1日至2011年7月31日每天
上午9时的在线监测进水COD数据,以其中4月
l号至7月24号的数据为原始数据,用于建模。
将收集的数据运用非参数检验法_2]检验其平
稳性,用的统计量为:
一
进猩 =鎏堡 游程标准差 望 一 二 z (]、
式中 一 +1研一LFZN,N 2(2
NIN ̄
-
N)
] ,
(N
—
1)
-J
N—N1+N2,N1为“+”出现的次数,为59;N2为
“
一
”出现的次数,为56;),为游程的总数,为58。经
计算,对于a:==0.05的显著性水平,l Z J一0.086
<1.96a :==10.46,认为该序列为平稳的时间序
列,并对其作零均值化处理,进而计算其自相关
(ACF)及偏相关(PACF)函数,会发现自相关及
偏相关函数均呈现拖尾性,故其适合ARMA
模型。
对平稳系统来说,可以用ARMA(n,n—1)模
型为一般形式建立模型l_3],通过对比最终选择的
预测模型为:
X一一0.354X,一1+0.923X 一2+0.386X 一3
+以 +口 一l (2)
此外,通过计算Ljung--Box统计量Q—N(N+
坚Ⅳ) ^
2)∑(N一是)1pl,式中N为样本数,志为滞后,
一
1
A
L(N)取18, 为{a }的估计值。得Q值为
17.857,小于矮1(14)一21.064及 。5(14)=
23.685,均表明该模型适用。
利用该模型对污水厂进水COD的拟合效果
如图1所示,拟合值的最小相对误差为0.26 ,
平均相对误差为2O.34 ,从图中两条曲线高度
接近可见利用所选模型对原有数据的拟合效果
较好。
犬毅
图1 9.1—9.24号观测值与拟合值对比图
用该模型对污水厂12.25—12.31号的进水
COD进行预测,其效果如图2所示。这里采用平
均相对误差和均等系数评价预测效果的优劣。预
测值的平均相对误差为12.94 ,均等系数为
0.9230,预测结果可视为满意。
l 80
160
140
—
】2o
县100
。80
0
u 60
40
20
0
天数
图2 12.25—12.31号观测值与预测值对比图
4反应器出水COD预测
我们在该污水厂搭建了一座自生动态膜生物
反应器中试装置,此反应器内的流态为完全混合
流[ ,中试系统中的有机底物进行物料平衡,
则有:
』0
S。Q—Qs + 一0 (3)
C
式中S。为进水有机物浓度;Q为处理水量;S 为
第4期 黄小红:基于ARMA模型的自生动态膜生物反应器中试装置出水CIOD预测
一
77
出水有机物浓度;V为反应器有效容积。
莫诺于1942年用纯种的微生物在单一底物
的培养基上进行了微生物增殖速率与底物浓度之
间关系的试验,得到莫诺方程 ]式:
Q 1+Q 3 +Q国 卜3+ +aid
~
——T
(式中以 S0代替(1)式中的x) (7)
即在一定的污泥浓度条件下,水力停留时间
基本恒定时,可根据由ARMA模型对进水COD
一 一x雳一VVma 丽XS丽
一
(4)
的预测,求得经反应器处理后的出水COD的预
式中 为有机底物降解速度;‰ 为有机底物的
最大比降解速度,t一;X为混合液中活性污泥总
量;S为有机底物浓度;K 为饱和常数。
自生动态膜生物反应器中,滤膜的表面会生
成滤饼生物膜,滤饼生物膜对出水COD有一定
的截阻能力,但这种截阻能达到假稳态(或稳定
化),故滤饼生物膜对出水COD的影响忽略不
计。最终可得如下关系式:
Q(So~S ) S 0-Se一 (5)
式中t为水力停留时间。将(4)式取倒数,按直线
方程y:aX+6考虑,利用图解法求解常数值Kz
Vm
==—
ax
—
K 。
利用11月份10天的实测反应装置进出水数
据求解。反应器有效容积为2.5m。。经计算得,
K2—8.13×10_。,则由(4)式得
S 一
1+KzXt
S0
So
1+8.13×10。Xt (6)
式中S。为进水COD浓度,(mg/L);X为污泥浓
度,(mg/L);t为水力停留时间,(d)。
由(5)式并结合进水COD预测模型,可得如
下关系式:
测值。
5结束语
对污水厂进水COD建立了ARMA预测模
型,取得了满意的预测效果,并将中试反应器出水
与污水厂进水COD之间建立了对应关系,可实
现对反应器出水COD的短期预测。有利于根据
预测值是否达标,判断反应器的工作状况,便于对
反应器的操作条件进行及时的、必要的调整,如调
整曝气量的大小、污泥负荷等来保证出水水质。
参考文献
1段智彬,孙恩昌,张延华,董燕基于ARMA模型的网
络流量预测EJ;.中国电子科学研究院学报,2009,4(4):
352—354.
2张树京,齐立心.时间序列分析简明教程EMJ.北京:清
华大学出版社,2003.
3王振龙,胡永宏.应用时间序列分析[M].北京:科学出
版社,2007.
4夏俊林.自生动态膜生物反应器处理城市污水的动态
研究[D].北京:北京科技大学硕士论文,2005.
5张自杰,林荣忱,金儒霖.排水工程下册(第4版)
[M].北京:中国建筑工业出版社,2006.
(责任编辑:王颖)