污水厂的稳定运行一直是每个污水厂工艺管理人员的终极梦想,但是在实际运行中,运行人员往往面临各种因素导致污水处理的运行不稳定,这些不稳定的起因分析、趋稳控制,水质达标是每次进行工艺调控的重要内容,每个污水厂的运维人员都希望能快速有效的解决污水厂的工艺问题。但是在实际工作中,很多污水厂再面对这些问题的时候,总是无所适从,盲目的做了很多不相关的工作之后,才能摸索出正确的方向,这样的流程过长,出水超标时间也长,现阶段的环保压力的情况下,是不容许污水厂长期进行这种摸索的。为了在污水厂的运行管理上探索更高效精准的控制技术,公众号从这一期开始,用几篇内容来和大家讨论下利用污水厂的数据分析,来指导污水厂的工艺运行的模式,希望能为运行人员从数据分析的角度来改善和调高污水厂的运行管理的效率。
( P. q h6 c7 \5 W! b* Z/ [
) b& q" o4 E: ~: P- ]& p8 D
* w: t0 \5 l0 G1 g
) V0 `% W/ W( ?/ W3 p' P在污水厂的运行中,很多时候会出现工艺的异常的变化,工艺管理人员需要进行准确的判断,并采取相应的措施。这些异常来自于污水厂运行的方方面面,有运行操作的,有设备故障的,有进水水质等等。这些异常在一个运行的污水厂中,是不断地发生并且存在的,运行人员需要掌控这些情况,并作出合理的判断。3 r- q% G" m% T, J: `
$ C3 e9 {3 S4 o& ^一线管理人员都会遇到在污水厂的运行过程中,出水水质从达标的状态忽然就变成超标了,从现场工艺运行操作上又没有特别的操作,或者没有自己认为影响到什么的操作。所以在分析中就会很困惑,这样的情况该如何去解决呢?那么这个时候就需要进行运行数据的分析,从运行数据中找到原因,并采取措施。下面以一个污水厂的运行实例来进行出水异常的情况的数据进行分析探讨。
- J5 ?7 J" R+ E) R
) l7 T' j6 f' v$ v$ Q2 V某污水厂地处西北省份地区,采用的运行工艺为A2O+MBR工艺,设计水量为2万吨,实际进水量平均为1.5~1.6万吨,7月份进水后,利用同地区污水厂污泥进行污泥培养,8月后稳定达标。从运行投产以来,出水水质一直保持稳定达标,进出水的各项指标如下表所示。
1 Z7 T# }$ H, Z+ p/ h3 C |* f3 t, r5 ?: ~: f' S# P
COD(mg/L) | NH3-N(mg/L) | TP ( mg/L) | TN( mg/L) | 水量(m3) | 进水 | 出水 | 进水 | 出水 | 进水 | 出水 | 进水 | 出水 | 进水 | 出水 | 280 | 35.24 | 23.98 | 0.36 | 2.965 | 0.64 | 29.2 | 9.41 | 16721 | 16316 | ( C3 o7 x% w9 |1 S! g" s; s+ C
但是从11月13日开始,出水总氮开始突然升高,从现场的工艺运行操作上上来看,有一些小的操作,但是对整体工艺都影响不大,对厂内的工艺运行,设备运行等等进行了注意排查后,没有发现重点影响总氮超标的问题。总氮开始超标后,厂内采取了增加回流,调整曝气等措施,但是调整的效果不是很明显,总氮仍然持续保持超标,但是其他指标变化并不是很明显。1 r1 S7 i9 I I* |; o
: m. K k1 `" T5 r* f
% u& L# ^9 `3 h5 @4 g+ i, t
8 G. w; m% s' M; n9 m- C/ B对于总氮的去除,在公众号之前的文章中进行了深入的讨论。那么从之前的讨论中我们可以进行一些预判,那就是总氮超标与几个环节有关系,分别是缺氧的反硝化环境,好氧硝化液的回流,反硝化的碳源的补充。根据厂内总氮的超标,围绕着几个方面进行逐一的分析。3 J$ h: h; H: t& s
9 u2 c2 d- U# x) a+ ~
厂内工艺采用的是A2O工艺加MBR,缺氧环境是通过空间上进行的划分实现的,在实际运行中对好氧末端的溶解氧的工艺控制也严格要求,虽然在MBR中采用气洗膜的方式增加了MBR池中的曝气,但是在前期运行中,这种溶解氧经过厌氧区的释放以后,对缺氧区的反硝化反应影响不大。而从硝化液回流的内回流的泵的运行上,近期没有异常的变化,设备由于投用期不长,运行工况良好,没有产生对工艺运行的影响。然后再看反硝化碳源的影响,碳源由于厂内没有每日化验BOD,一般以COD为判断依据,据厂内人员介绍,近期由于上游关停了部分企业,可能会产生一定的影响,为了进一步判断是否由于水质变化造成的影响,对厂内的各项水水质数据进行了统计分析。
8 V1 U; f9 d4 L6 i7 b" c' a
# D" c; O/ E! u0 }2 `, T/ L首先对工艺人员怀疑的进水COD进行分析,通过表格软件对进水的COD进行曲线分析可以看到进水的COD变化情况。
( J6 W( x' ^+ e2 t3 i P2 J* `1 I( S! V
7 Z$ u1 m/ J. f% [# f2 x) I5 t
. ?3 x- I4 Z1 q! B
可以从数据的统计曲线上看到COD的变化有,但是和之前的月份相比,这种变化没有特别明显,为了更好的分析,我们增加了移动平均的趋势线来进一步分析变化的趋势。为了更明显的看到变化趋势,我们把趋势线加粗,数据线变细。9 ^! i# u& @0 o5 f; o% \7 ]
0 s" E; Q- b4 m7 G' x
5 I. [9 S7 h9 F
- \+ J; \; B# z3 ?; @ R5 e' j8 L r
从变化的趋势图来看,进水的COD变化近期内在逐步下降,但是从前期的运行变化的情况上来看,还不能说明反硝化的碳源不足,那么什么原因才是导致总氮超标的原因呢?' R6 a$ g# c2 N4 x0 W& N, |# p
7 n5 Y9 W; f& o8 K$ v* N- ]) K
为了进一步分析,我们把进水的总氮也进行了统计分析。分析图表如下:2 m/ v- y& Q, j7 u# M) O9 l5 X
0 c, l3 _1 |# t( Q
, b; ]1 W: G- R' S2 v- s4 e9 A; Q
* o! P y: F" f" K/ F0 L& i
为了更明显的看到总氮的变化,我们把坐标轴变化一下,坐标起点从10mg/L开始。& s! l8 ?& {1 O3 ]
8 v! h( x* K6 O1 Q; z9 o
8 I% J- N# S8 c, F
, o% S5 e& x" t+ B1 N* N- I
这张图可以更明显的看到进水总氮呈现了一个逐步升高的变化趋势,如果两个曲线叠加之后,是不是能更好的说明问题?利用图表软件里面的一些小功能,我们把进水总氮的变化曲线背景变成透明的,然后进行两张图的叠加,我们能得到这样的图:
9 \% w0 L; N5 _5 i) M/ O+ l6 r3 f; U+ }
' v! B: T2 G4 q) k# q5 H) @1 Y* Y4 k
3 F# p6 _) Z* p7 d4 C: z
注意看两个趋势线的变化情况,总氮是一个逐步增高的趋势,COD是一个逐步下降的趋势,两个曲线在之前的运行中,变化基本上是同步的,增加或者降低都是一致的,但是在曲线上显示从横坐标89开始,也就是从10月30日左右开始出现了两条曲线相互背离的趋势,总氮在持续升高,COD下降。从这张曲线上已经明确显示出相对于之前的运行来说,进水的碳源开始出现不足,导致了出水的总氮开始超标。3 G( {! `0 j0 n: |% Z
( y/ f7 E3 a2 u) ] u$ A* i1 r6 x
2 s8 y8 ~1 r3 I* v! h
3 g( d- W4 g- V% P5 w7 q, }! }3 a为了更明显,可以用颜色区分开他们的变化,当然也可以用更高级的数学工具去分析两者之间的变化,对于污水厂的运行人员来说,简单的表格工具已经可以分析出原因了,也欢迎专业人员进行更深入的数据分析。- u1 p8 j5 J$ x
( k# a! y1 j$ p% X+ w( v7 O2 }
这是对某个厂的总氮超标进行的数据曲线变化的分析,根据这个分析,我们可以明显看到影响超标的主要原因来源是进水水质的变化,但不是进水的碳源下降,而是进水碳源下降,伴随着进水总氮上升,造成的进水碳源的相对不足。在这种分析的结果下,后期的工艺调整或者工艺关注点,工艺运行人员就很能明确了。
, h x }7 y \" H) X9 f
3 Y2 E. K L% M: X+ c& x6 `( Q1 ]数据分析作为现今社会重要的科学手段,对现代社会起到了非常重要的作用,无论从科技发展,工农业发展,商业应用中都有非常多的应用。和我们老百姓常见的就是大数据无时无刻的影响着我们的日常生活,购物,出行,日常生活都会被大数据包围,我们也越来越享受大数据带给我们便捷的生活,而这些大数据背后就是各种各样深入的数据归纳和分析。作为污水厂,其实每天都在不断的产生大量的数据,这些数据对于污水厂的管理,能不能像大数据一样,使我们的管理更加便捷,高效呢?所以我们是需要对污水厂产生的数据进行有效的分析,从而使污水厂的运行管理更加精准和到位。- l, v# ], y0 {' r+ L
0 |5 C& ~/ N3 G R: K! b9 c那么数据分析是如何影响到污水厂的运行管理的呢?我们先来了解数据分析,在没有确定我们需要通过数据分析来了解什么的时候,这些来自于生产一线的数据往往只是一堆数字或者文字,没有任何的贡献。但其实这些数字是来源于每一个污水厂的实际运行的情况,这些数据背后隐藏了每个污水厂的运行管理方方面面,包括活性污泥的各种生存状况。而我们需要做的就是通过数学手段进行数据分析,把这些冷冰冰的数据背后鲜活的污水厂运行的情况挖掘出来,成为我们运行管理的扎实基础。
0 B- F1 N5 {6 m( I% v. O# x6 y. |
下面我们来举一个例子,看看如何用数学手段挖掘数据之间的关联性,对污水厂的运行管理做出帮助。
( Q& d* ~; I3 Q& S4 t2 G# S
# }1 a% s0 g* r# q我们知道在一个污水厂里,进水水质是非常重要的,进水水质中的BOD更为重要,特别是近年来对出水氮磷的要求,从普遍理论上,去除氮磷是对进水中的BOD要有一定的比例关系的。但是对于一般的污水厂来说,BOD的化验难度远远大于COD,很多污水厂都不会去每日进行BOD的化验,但是会长期进行COD的化验。这就出现一个问题,我们能不能用COD来替代BOD,进行每日的进水水质的评估呢?如果进水的COD是400mg/L,我们能不能预测出进水的BOD是多少呢?我们来用一个污水厂的数据进行讨论。+ t! `: a' c6 K8 X7 ] Y. d
) `$ O, Z ^' ?$ T
" M6 h% y* ~7 }! o/ {7 a
' c" F: L7 H3 {% e) ?这是某个污水厂在2004~2007年的COD和BOD的化验数据得月平均数值,一般来说,污水厂由于是专属服务某一个地区的污水处理的,这个地区由于居民生活水平,工业比例,管网建设情况在一定时期内稳定,那么也就是说这个地区的污水的组分是固定的,我们需要通过这些化验数据来做出一个COD和BOD之间的关联关系来,方便我们后期更加快捷方便的通过COD的化验就能预测到BOD的数值,来衡量进水中的碳源。
) X. ~# x2 [. Z5 T. d, I# {8 n Q1 L5 y1 t P. o' S
首先利用EXCEL中的作图功能,这种寻求两个数据之间的关联的方法,一般采用这种散点图,散点图一般用于比较至少两组数值或者两对数据,使用它可以显示值集之间的关系,当我们把某污水厂2004~2007年的COD和BOD月平均值值得数据做成散点图之后,大家可以看到,这个散点图非常明显的显示出两者之间的关系呈线性分布,那么我们现在需要做的就是如何把这个线性关系找到。
`+ T' }, S P! \# M' n- k$ L* \, o# W# r0 h6 ?
2 v7 g; P8 m. S& }1 g, {- s
, E# P' X: z- Q6 p9 Y# u这种线性关系得寻找,在计算机辅助下,可以非常方便快捷的得到,大家可以利用EXCEL里面得数据趋势线的添加方式来进行,选择添加趋势线,在趋势线的选择框里,选择线性,注意下面要把“显示公式”得选项选择上,这样这个公式就能很方便的得出来了。
$ I$ j8 ?' K" l# O' E2 g2 d0 V/ [4 E' V- j9 o3 Z
- _. }/ S9 W ~
7 g, {& K, u5 F- }. N当我们通过几年的数据进行分析后,得出得这个公式8 ?4 x% H8 X1 |& H
* z6 p5 Y3 ~' Z' `0 R' w' s7 y
Y=0.4366X+45.574,换成BOD和COD的关系就是:6 i- _0 v. n1 @% e
, _7 O+ b" I3 r M% O BOD=0.4366COD+45.574
$ U) y Y: l7 s; w# ]8 q- Q! @& K0 C2 j0 X' H* U
能够通过COD来预估进水的BOD,从而研判进水中得碳源是否充足了。
8 |% N. V: v P! T( F
; y; N4 `8 Y; g, T8 O$ \当然这个是快捷得方式,如果是要采用数学推理计算的方式,这种线性关联的方式就要通过线性回归方程来进行计算。4 f5 ]$ U- W. n6 d
4 Q+ L0 ]+ z: W: q8 Z% H% ?, v
第一:用所给样本求出两个相关变量的(算术)平均值:
+ J7 P9 g G5 f' \( ?0 d0 {
2 E2 u3 ^. X: a0 p# rxp=(x1+x2+x3+...+xn)/n) X4 f6 {) p* [. k; r! m
5 }' x1 w$ X7 r5 y9 O& @* E
yp=(y1+y2+y3+...+yn)/n
" y9 _1 a& T* R" v/ `3 n5 }; i! b& L, y. C: x
第二:分别计算分子和分母:(两个公式任选其一)
/ a- S* h8 E: N6 e
$ C( } Y: j/ t7 X) [' _7 Y& X+ H分子=(x1y1+x2y2+x3y3+...+xnyn)-nXpYp r$ W. X7 a' r& s* D
' _4 I& @: J0 J# g3 k8 N7 F/ a( b分母=(x12+x22+x32+...+xn2)-n*x2& \" D, s7 L* [5 U" [
2 q9 K) ` Y' @' O6 e+ b8 B; H
第三:计算b: b=分子/分母
5 g+ O# x" e( V4 `4 T i- v2 ?( N. y: ^- c! {9 J$ t
第四:最后把x,y的平均数xp,yp代入a=Yp-bXp
( m9 B8 w( \7 K$ C7 d/ ]
: F5 `8 K% p* N) G求出a并代入总的公式y=bx+a得到线性回归方程。
+ t7 i. P1 v: t& ?) H$ _' n' k+ ?* V7 [$ D) q; @
为了验证这个数学计算方法,数学方法要求数据样本更多,更多的数据样本可以反应出更准确的数值,我们把某厂的2007年的每日的COD和BOD的值进行统计计算,出于篇幅原因,我们把数据表的部分列出来,然后用EXCEL里面的公式进行上述的数学运算,计算的结果见表后的数据:
0 ]7 b# {6 c0 x3 x! w% U! }. D: `( t# g( u* q' T/ {( l
3 o* |4 C, d0 n9 Z* s0 F
4 A5 z9 [5 \' e& t' x1 U最终计算出2007年的COD和BOD的关系式为:9 x$ Z4 D/ p# g/ U1 l
' z/ Q7 y( q$ x+ w3 Xy=0.3455X+85.56& y. S, z# Y! Y5 M; n A" e
* E# B7 e2 t9 C# w也就是BOD=0.3455COD+85.56 * K" s% C2 T" O
0 @* o3 J s$ B2 X2 ~6 a/ l; b然后我们用EXCEL的散点图和趋势线的方式来进行验证,可以看到这个曲线图的方程是:1 X/ F9 k- I9 \$ g
( z0 `$ u- L+ ]- j) T. l8 b实战:数据分析对污水厂工艺调控的作用 - 环保之家& I, ?9 E% E7 t0 p% z5 T+ A
8 l) f, I0 L u$ w这个方程和我们数学计算的方程基本一致,因此可以知道数学方法和表格方法都是可以完成这个线性方程式的推算的。这个线性方程的得出,可以便捷的通过进水的COD来快捷的判断进水的BOD值,不必要等到5天后出来BOD5再来判断系统的碳源是否充足。8 s& v8 l# L' p/ c' R$ b# d3 ~
0 k) e! W" R7 r# t) l) T |
-
© 声明:本文仅表作者或发布者个人观点,与环保之家[2TECH.CN]无关。其原创性及陈述文字、内容、数据及图片均未经证实,对本文及其全部或部分内容、图片、文字的真实性、完整性、及时性本站不作任何保证或承诺,仅做参考并自行核实。如有侵权,请联系我们处理,在此深表歉意。
|