 鲜花( 0)  鸡蛋( 0)
|
摘 要 报导了锰砂净水 石英砂净水过程中使用计算机实时采集和定量分析絮凝池中絮体图像,算出絮体沉淀“等效直径”作为控制聚合氯化铝加注率的目标值,并用进水流量和沉淀水浊度信号自动修正目标设定值。加注量为进水流量与加注率之积。实验数据表明,该方法能适应水质和其他因素的变化自动调整加注量,达到了稳定沉淀水浊度的目的。
! d* Z2 ]2 T1 P+ C 关键词 絮体图像 等效直径 聚合氯化铝加注量 计算机控制 5 B! o7 t% }- c) ~) R7 ~
在锰砂净水 石英砂净水自动化中,聚合氯化铝加注量的自动控制是一个既关键又难解决的问题。因为影响聚合氯化铝加注量的因素很多,如原水浊度、水温、流量、碱度、氨氮、耗氧量和锰砂净水 石英砂净水设备的负荷、状态等,目前国内外使用的方法尚不完善。因此,需寻求效果好、适用范围广、性能先进、运行方便、造价低、维护简单的自动控制方式。
7 ?* ?! ]5 d: S" `9 W# |4 [0 m5 t2 U
1 定量分析絮体形状
2 i) i! X1 i e
8 U' ~4 s( @# A" _& Y9 _) Z 确定适当的聚合氯化铝加注量,关键是要找出一个滞后时间较短,而与沉淀水浊度相关性又较好的参量作为目标值来控制。" ?/ T; |/ L3 f
1.1 絮体沉淀特性和沉淀水浊度的关系
+ o S, U; `+ T& p. g 从锰砂净水 石英砂净水过程可知,沉淀水浊度与原水加聚合氯化铝后形成的絮体特征和沉淀有关,絮体形成得越好,沉淀越充分,沉淀水浊度越小。在一定沉淀条件下,沉淀水浊度和絮体的沉淀特性密切相关。 ]: M: B9 }- w+ _' w8 `3 V3 f
絮体的沉降规律是比较复杂的,常简化用颗粒沉降的Stokes公式来描述:8 R8 N7 r) a0 Y% R$ w2 p4 t
1 R) v& {( `% }# C3 {5 X4 @ v=(ρs-ρ)gds2/18μ (1)* d, W7 g- ]" c7 D- U7 w y5 d
式中 v——絮体沉降速度,cm/s: d2 q/ |8 G) X
ρs——絮体体积质量,g/cm32 _# z1 K5 g3 k( g
ρ——水的体积质量,g/cm3: B4 V2 R5 [5 _9 e1 f; ]
dS——絮体直径,cm
$ z9 ] U) {2 x$ e- l' c) B6 d; D μ——水的粘滞系数,g/(cm·s)6 t0 H. Y! C7 J& E7 M0 y; S7 O$ ]6 J
g——重力加速度,980 cm/s27 _6 c) I! G) R, C# `. P
进一步的研究[1~6]表明,絮体粒径增加时,体积质量相应减小,其关系式为:) A; H* d4 Z% I& G( X; \
ρS-ρ=dS-kp (2)
5 i+ x8 [ |5 R& d 式中kP——系数,1.2~1.5,决定于聚合氯化铝加注率与原水水质6 i2 q' G8 X ~0 e: a% D
5 s: D2 h) R/ X" ^! D$ @. N; P4 u d
v=gds(2-kp)/18μ (3)0 E$ Z' b) \6 j9 {6 d
. O9 @9 n: L; f2 U7 c 上述分析均假设絮体为球状颗粒,而实际絮体基本上是不规则状态,其沉降速度显然应比同体积的球状絮体慢一些。絮体的大小、形状可反映在絮体图像上,因此通过分析絮体的图像,可以得到一个与沉淀水浊度相关性很好的参量。用它作为目标值来控制聚合氯化铝加注量可使滞后时间大大缩短。' W! y) @. X m* O$ s2 H
1.2 定量分析絮体图形,计算等效直径; X( v6 e* [; K4 \& x9 G: _6 m
为了从絮体的二维图形综合出与絮体沉淀速度有关的参量,给絮体图形定出了以下四个特征:①表示絮体大小的絮体面积s;②与絮体形状有关的絮体周长l;③与絮体松散程度有关的絮体中间空出面积s0;④絮体的长宽比m。这些特征基本表示了絮体的特性,且易于计算。最后按下式折算成称之为“等效直径”Φ的参量: }" P$ F, C e. C1 Z( y
2 s8 h2 d! e$ ? Φ=2(s/π)0.5×[1-k1(1-2(sπ)/l)]×[1-k2(1-1/m)]×(1-k3s0/s) (4)9 p/ N4 W4 R3 e8 h, J) \2 S
0 K$ m# e5 S2 V, w& u o
式中 k1、k2、k3——周长、长宽比、中空面积的折扣系数,是0~1的小数,为0时不打折扣,为1时折扣最大,可根据实际情况选取Φ越大,沉降速度越快。当絮体为一标准圆形时,式(4)的等效直径即为实际直径。4 P) m0 h$ ~, i# g
絮体的图像是通过传感头在絮凝池中直接采集的,水不停地流过传感头的取样窗,取样窗的水域面积为26mm×20mm,厚度仅为3.5mm,以减少二维图像中两个絮体重叠的可能性。在该截面中一般包含了数十到上百个絮体。系统每5s采集一幅图像,按式(4)计算出每个絮体的等效直径,每5min得到60幅图像中所有絮体的等效直径及其数值分布情况,取其分布中某一部分加权平均算出平均等效直径。试验表明,在沉淀条件不变的情况下,平均等效直径与沉淀水浊度有很好的相关性。5 Y$ O8 i) U% ?1 t, i
1.3 等效直径控制加注率4 ^4 g' A3 U. \; ~ |/ K$ J. r& u
加注率是等效直径的实测值与设置值之差,通过数字PID(比例、积分、微分)运算后得到的,其递推式为:; X9 q" f+ w2 @7 ~5 G: E( [
8 J" N* A5 q+ K# L# [9 k
Δp(k)=p(k)-p(k-1)=kp[e(k)-e(k-1)]+kie(k)+kd[e(k)-2e(k-1)+e(k-2)] (5)
1 V7 t5 b j& s, e- ^" H 式中 Δp(k)——第k次采样时加注率的修正值
, y0 i# y# ]) w1 d; L p(k-1)——第(k-1)次采样时的加注率
+ s& T; r1 C0 j# i& ^( O! c e(k-1)——第(k-1)次采样时的偏差
9 Y3 ?/ x6 e- D( T6 M p(k)——第k次采样时的加注率
3 H8 }- G- ^* j4 _ e(k)——第k次采样时的偏差
8 B# p! V0 b4 r4 H, p$ u) n# R kp——比例系数3 E# ~( b9 P; t! l! M/ I- `
ki——积分系数 g% y) ?1 F& Y, w! m# T+ J
kd——微分系数% \+ u: h8 `9 P1 R) z
由于聚合氯化铝的加注到絮体的形成有一定的滞后,仅用一般的PID算法难以得到较好的动态性能,因此还增加了Smith预估控制,在假设系统为一阶惯性纯滞后系统的前提下,该方法增加一纯滞后补偿环节,该补偿环节和原系统一起构成的广义对象不再具有纯滞后,利用广义对象的输出信号作为反馈信号进行控制,可解决纯滞后问题。具体做法是在偏差计算时增加一修正项e1S(k):5 M1 p2 ~5 a# O1 t
0 Z) m1 n! T1 a9 e' K5 ?. u9 S! C e(k)=e1(k)+e1S(k) (6)
/ S3 r& f3 q2 G3 H 式中 e1(k)——实际偏差
9 e. v- s4 g+ ^5 r4 _* V M' W# q( r- y: o) ]" S, E2 N" y0 \
e1S(k)的算式为:, d: N Q6 l) j, t" s
, Z& N( e E( y: ]& O) }' B3 B e1S(k)=K[ν(k)-ν(k-1)] (7)
2 e" D* d6 u+ b* E( I/ ~9 \ v(k)=(1-T/Tm)v(k-1)+p(k-1)T/Tm (8)& u) V2 r/ ~ Q2 {" H
式中 T——采样周期
& |! |4 ~( [. H6 K9 s K——系统增益
* z: E# X" h( j+ _% ? Tm——系统惯性时间常数+ U' @ A0 g/ f/ C; ^
l与系统滞后时间τ的关系为:
4 L5 w0 S) o+ m1 S( r! D1 t2 _1 [ τ=lT (9)
: {5 q% l1 `5 Z5 l( c 第k次采样时的加注量Q(k)为:: J. h* S+ |& v; Y5 }' |2 Y& b: H
Q(k)=p(k)×L(k) (10)
8 c* M4 L( ~8 I3 l! P. s0 `# ] 式中L(k)——第k次采样时的进水流量 z ~+ A2 V' ]* ]8 C+ c6 w
1.4 条件变化的反馈自动调整
$ F, g: V, d) y3 j7 ^8 [ 沉淀条件不变时,等效直径与沉淀水浊度有很好的对应关系,但当沉淀条件变化时,等效直径与沉淀水浊度的对应关系会有变化。这时就需要调整等效直径的设定值,也就是说该方法和前述多数方法一样,不能模拟沉淀池。为此,碱式氯化铝,当流量有变化时,每次采样该系统能自动调整等效直径的设定值,以消除流量变化引起的沉淀条件变化对设定值的影响。一般除流量外,沉淀条件的变化都较缓慢,为此该系统还使用沉淀水浊度反馈来自动调整设定值。做法是每30min将沉淀水浊度的测量值与目标值之差加上Smith预估修正项,经PID运算后改变等效直径的设定值。' j- i- m. ]1 y4 h
; d/ {; X( Z+ O8 a! q2 Z2 控制系统的硬件和软件
6 O/ F r* `. ` w* y" a) z! j3 f) V
( }# j. k) h6 h 系统硬件如图1所示,586主机通过图像接口将絮体图像信号数字化后送入内存;通过模拟接口采集4~20mA进水流量和沉淀水浊度信号,并输出4~20mA电流信号控制聚合氯化铝加注泵。6 ?4 t- p4 ?% G% l3 u. N) b
' A& u8 t2 L* X
软件的主要功能是:
/ R0 K; {' \( h) ~( F ① 将采集的絮体活动图像实时显示在计算机屏幕上;9 B5 `4 o. G0 J
② 对絮体图像进行边缘增强、数字滤波、二值化处理、连通性判别,算出每个絮体的s,l,s0,m,最后按式(4)算出絮体的平均等效直径;
! D% Y/ W+ k( [" g ③ 采集进水流量、沉淀水浊度信号;
: N0 o0 ~! A& ?/ X8 I ④ 按式(5)~(9)算出聚合氯化铝加注量,并通过模拟接口输出;3 ~9 `5 }4 B5 C0 U5 y0 e
⑤ 在屏幕上显示采集和计算出的各种数据并实时更新;
9 n8 p9 g! v; c ⑥ 各种参数如P、I、D参数,系统延迟时间,惯性时间,等效直径,沉淀水浊度的设定值,以及絮体图像的对比度和亮度等,都可通过下拉式菜单自行设定,以适应不同生产设备和工艺的需要;2 C! m# J2 K3 E+ E! K
⑦ 所测得和计算出的结果及时间、日期等数据,每5min一次自动存入硬盘,可存10年。) ^+ P( ]& ?3 P8 n, y
) [( g. B: C+ ?
3 试验结果
4 k K# c9 g) \( g" V* _8 v- Y
& Y% n( R& m& J: n 本系统在上海南汇县航头水厂进行了现场试验,试验池为15×104t/d的隔板絮凝水平沉淀池,第一阶段只记录絮体数据、沉淀水浊度、进水流量、聚合氯化铝加注量等参数,不控制聚合氯化铝加注量。然后分析这些数据,在使等效直径与沉淀水浊度相关性最好的前提下,得到计算平均等效直径的参数和统计方法。第二阶段以所得参数和方法计算等效直径,按式(5)~(9)进行加注量控制。按该厂生产要求,沉淀水浊度控制设定值为5NTU,4~6 NTU为符合要求。图(2)为1998年4月25日—5月1日一周内进行加注量控制后的沉淀水浊度、流量、絮体等效直径与加注量的对应关系曲线图。经每5min一次数据统计,一周内沉淀水浊度高于6 NTU的有72次,占3.57%;低于4NTU的有20次,占0.99%;沉淀水浊度总合格率为95.44%,取得了较好的控制效果。从图(2)可看出,4月27日、28日和5月1日进水流量无大变化,但加注量却从加注泵冲程的70%左右降到了50%左右,而沉淀水浊度基本相同,说明本系统能适应水质和其他因素的变化,有效控制加注量,稳定沉淀水浊度。
8 x( S9 j! E5 x
0 w' N0 c5 m+ A7 q" c4 结论% O" ]0 i7 ^8 h, L: z
" U1 k& \0 m, v* m% N4 m ① 使用计算机实时采集和定量分析絮体图像并算出加注率的方法,能有效地控制聚合氯化铝的加注量,将沉淀水浊度稳定在一定范围内,达到保证水质、降低生产成本的目的。虽然实时图形分析运算量非常大,但目前计算机技术的发展已完全能满足运算速度要求,且成本也较低。# A' ^" D9 o5 r# @/ R/ [
9 `" A; Q0 w, Z+ x5 @0 N
② 从聚合氯化铝加注到絮体形成一般需10~20min,这就是滞后时间,比流动电流测定(SCD)法要长。但因流量因素已单独处理,滞后时间已能适应水质变化的要求,且使用了Smith预估控制,系统动态性能有一定改善,满足了使用要求。一般检测的参量越靠前,滞后时间越短,对该参量后系统的模拟性就越差。SCD的测量在絮凝前,因此不能模拟絮凝和沉淀条件,当絮凝或沉淀条件变化时,设定值都需要改变[1]。因此选择检测参量时,滞后时间应综合考虑,不能绝对地说越短越好。. x! J' |% m2 X. c8 M$ `
③ 本文所述的“等效直径”,只是用较简单的运算尽可能较准确地来描述絮体沉淀特性的一个参数,并不表示与某直径物体的沉淀特性“等效”。本系统使用计算机定量分析绒体图形,可得到各种絮体参数供统计研究使用,随着生产数据的大量积累,等效直径的算法也可不断改进和完善。 |
|