最新消息: USBMI致力于为网友们分享Windows、安卓、IOS等主流手机系统相关的资讯以及评测、同时提供相关教程、应用、软件下载等服务。

基于ARMA模型的自生动态膜生物反应器中试装置出水COD预测

IT圈 admin 31浏览 0评论

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 

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 

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. 

(责任编辑:王颖) 

发布评论

评论列表 (0)

  1. 暂无评论