基于數(shù)據(jù)分析的自來水廠參數(shù)自適應(yīng)次氯酸鈉投加算法*

彭一恒 章紅 劉衛(wèi)斌

開發(fā)設(shè)計

基于數(shù)據(jù)分析的自來水廠參數(shù)自適應(yīng)次氯酸鈉投加算法*

彭一恒1章紅1劉衛(wèi)斌2

（1.江漢大學(xué)物理與信息工程學(xué)院，湖北 武漢 430056 2.武漢市白鶴嘴水廠，湖北 武漢 430040）

針對某自來水廠的實際情況，通過數(shù)據(jù)分析獲得次氯酸鈉投加的特征和規(guī)律，設(shè)計參數(shù)自適應(yīng)次氯酸鈉投加算法?；谀乘畯S制水工藝，利用歷史數(shù)據(jù)庫中進(jìn)出水流量、清水庫液位、源水濁度和反沖洗狀態(tài)等參數(shù)，采用趨勢擬合等方法確定各參數(shù)與次氯酸鈉投加流量的關(guān)系，并將次氯酸鈉的投加控制轉(zhuǎn)換成清水庫的液位預(yù)測控制。經(jīng)仿真，該算法理論庫氯值與實際庫氯值的平均絕對誤差為0.060 mg/L，可實現(xiàn)對庫氯的穩(wěn)定控制。

自動控制；次氯酸鈉投加；數(shù)據(jù)分析

0　引言

加氯環(huán)節(jié)作為水廠制水工藝最重要的環(huán)節(jié)之一，起到消毒滅菌作用。目前，國內(nèi)大多數(shù)水廠采用手動加氯方法，其效率和成本難以保證。而現(xiàn)階段的自動加氯方法大都采用PLC設(shè)備及PID算法[1-5]，大滯后的加氯控制系統(tǒng)穩(wěn)定性和可靠性無法保證。隨著生活水平的提高，人們對水質(zhì)的要求也不斷提升，高效率、低成本、可靠和穩(wěn)定的自動加氯控制系統(tǒng)成為當(dāng)前趨勢。自來水通常采用次氯酸鈉消毒滅菌。次氯酸鈉相較于氯氣的消毒設(shè)備更簡單，安全性更高[6]。

本文以某水廠的升級改造為契機，針對該水廠次氯酸鈉投加控制要求，設(shè)計基于數(shù)據(jù)分析的自來水廠參數(shù)自適應(yīng)次氯酸鈉投加算法。首先，介紹該水廠的制水工藝流程，通過其歷史數(shù)據(jù)，分析次氯酸鈉投加特征和規(guī)律，發(fā)現(xiàn)改造前次氯酸鈉投加控制缺陷；然后，提出參數(shù)自適應(yīng)次氯酸鈉投加算法，通過進(jìn)出水流量、清水庫液位、源水濁度和反沖洗狀態(tài)等參數(shù)預(yù)測清水庫液位，從而實現(xiàn)次氯酸鈉投加流量自動控制；最后，經(jīng)仿真和實際運行證明該算法可行。

1　水廠消毒滅菌過程

改造前，某水廠采用傳統(tǒng)的手動加氯氣方式實現(xiàn)自來水的消毒滅菌。改造過程中，利用10%的次氯酸鈉溶液替代氯氣，與傳統(tǒng)加氯氣相比，消毒能力基本相同，但提高了安全性，降低了運輸成本，增加了消毒反應(yīng)時間。次氯酸鈉水解后分解成次氯酸根和鈉離子，次氯酸根與水中氫離子結(jié)合形成次氯酸，起到消毒作用；與細(xì)菌反應(yīng)主要生成NH2Cl，NHCl2及NHCl3，稱為化合氯。未反應(yīng)的氯主要存在形式為OCl+，HOCl，Cl2，稱為游離氯，游離氯在清水庫和出廠水中分別稱為庫氯和余氯。

出廠自來水要求余氯在0.9 mg/L~1.2 mg/L內(nèi)，以達(dá)到自來水流入城市管網(wǎng)中還有消毒滅菌作用。某水廠制水工藝流程如圖1所示，大致可分為3個階段：

1）源水消毒沉淀階段：在源水池中加入礬和臭氧，臭氧可以氧化源水中多種有機物和無機物，滅毒殺菌；礬可以與源水中的雜質(zhì)結(jié)合，形成沉淀，此階段主要目的是對源水進(jìn)行沉淀和消毒；

2）過濾階段：在石英砂、碳濾池中，用石英砂濾除礬無法結(jié)合沉淀的雜質(zhì)，并通過活性炭吸附雜質(zhì)，此階段的主要作用是去除水中剩余雜質(zhì)；

3）清水消毒階段：通過變頻泵向清水庫和調(diào)水庫投加次氯酸鈉溶液，完成對出廠自來水的消毒滅菌。

該廠改造前，次氯酸鈉投加以余氯控制為主，庫氯控制為輔，自來水中的氯含量控制不夠穩(wěn)定，控制過程存在滯后，可能出現(xiàn)余氯過高或過低情況，影響出廠自來水水質(zhì)。

本文設(shè)計為以庫氯控制為主，采用參數(shù)自適應(yīng)次氯酸鈉投加算法，既可提高消毒滅菌效率又能保證出廠自來水的安全。

圖1　某水廠制水工藝流程圖

2　次氯酸鈉投加特征分析

改造前，該水廠次氯酸鈉投加通過調(diào)整各次氯酸鈉投加變頻泵的開度，控制出廠自來水余氯在要求范圍內(nèi)，控制效率低，穩(wěn)定性差。改造后，通過上位機監(jiān)控界面可實時監(jiān)控次氯酸鈉投加參數(shù)的變化；同時，可以通過上位機內(nèi)部算法調(diào)整各次氯酸鈉投加變頻泵工作，以提高控制效率和控制穩(wěn)定性。

本文算法設(shè)計過程中，收集該水廠2019年6月~2019年8月近十萬條數(shù)據(jù)，包括進(jìn)出水流量、清水庫和調(diào)水庫補加氯的次氯酸鈉投加流量、變頻泵開度、庫氯、余氯、清水庫液位、源水濁度和反沖洗狀態(tài)等。通過對比分析發(fā)現(xiàn)：進(jìn)出水流量、變頻泵開度、源水濁度和次氯酸鈉投加控制關(guān)聯(lián)度較強。具體分析如下：

1）次氯酸鈉投加對庫氯的影響存在延時

以2019年8月15日的清水庫2為例，次氯酸鈉投加流量與庫氯變化關(guān)系圖如圖2所示，進(jìn)出水流量折線圖如圖3所示。為更直觀分析，庫氯做擴大100倍處理。由圖3可知：在7:00~23:00之間進(jìn)出水流量無明顯變化。分析圖2中11:00~17:30和18:15~20:40兩個時間段次氯酸鈉投加流量變化與庫氯變化之間的關(guān)系可知：當(dāng)設(shè)定次氯酸鈉投加流量變化時，啟動變頻泵工作，經(jīng)過25 min左右?guī)炻戎挡砰_始變化，再經(jīng)過20 min庫氯值才基本穩(wěn)定。因此，次氯酸鈉投加與庫氯值變化存在延時，延時響應(yīng)時間約25 min，庫氯值達(dá)到穩(wěn)定時間約45 min。

圖2　2019年8月15日次氯酸鈉投加流量與庫氯變化關(guān)系圖

圖3　2019年8月15日進(jìn)出水流量折線圖

2）進(jìn)出水流量可預(yù)測次氯酸鈉投加流量

圖4和圖5是2019年8月27日和28日兩天的進(jìn)出水流量與庫氯變化曲線。為更直觀分析，庫氯做擴大15000倍處理。進(jìn)出水流量的變化與城市生活用水量密切相關(guān)，當(dāng)用水量增加時，需加大供水，出水流量增大，要求制水量增加，從而進(jìn)水量增加（加車）；當(dāng)用水量減小時，需減少供水，出水量降低，可減少制水量，從而減少進(jìn)水量（減車）。

圖4　2019年8月27日進(jìn)出水流量與庫氯變化曲線

圖5　2019年8月28日進(jìn)出水流量與庫氯變化曲線

對比圖4和圖5可知：加車時間一般在清晨，減車時間一般在深夜，但不同日期的加減車時間不盡相同。再進(jìn)一步分析圖4和圖5的庫氯變化曲線，發(fā)現(xiàn)加減車時，庫氯延時一段時間后出現(xiàn)明顯變化，這主要是因為加減車會使清水庫液位發(fā)生明顯變化，從而影響庫氯值的變化。若要保證庫氯值穩(wěn)定，需改變次氯酸鈉投加流量。由此可見，根據(jù)進(jìn)出水流量可預(yù)知清水庫液位變化，從而預(yù)測庫氯穩(wěn)定情況下次氯酸鈉投加流量。綜上分析，該水廠的庫氯變化與進(jìn)出水流量有關(guān)，而次氯酸鈉投加與庫氯變化之間存在滯后。如果能夠進(jìn)一步明確進(jìn)出水流量與庫氯變化之間的關(guān)系，即可通過控制次氯酸鈉投加實現(xiàn)對庫氯的控制。

3　算法設(shè)計

該水廠清水庫次氯酸鈉投加反應(yīng)具有大滯后、延時特性，通過進(jìn)出水流量可預(yù)測清水庫液位的變化，進(jìn)而可預(yù)測次氯酸鈉投加流量。因此，次氯酸鈉投加算法將次氯酸鈉投加流量問題轉(zhuǎn)化為清水庫液位預(yù)測控制問題，解決次氯酸鈉投加大滯后和延時問題。具體思路為：根據(jù)進(jìn)出水流量、源水濁度、投礬流量、反沖洗時間、加次氯酸鈉流量、實際液位等參數(shù)，預(yù)測清水庫液位；再根據(jù)庫氯期望值（水廠工程師根據(jù)經(jīng)驗，將庫氯期望值設(shè)定在0.6 mg/L~0.8 mg/L，本文設(shè)定為0.7 mg/L），確定次氯酸鈉投加流量。次氯酸鈉投加控制系統(tǒng)結(jié)構(gòu)圖如圖6所示。次氯酸鈉投加控制器由參數(shù)融合模塊和投加流量計算模塊組成，其中參數(shù)融合模塊主要完成傳感器采集參數(shù)的融合計算；投加流量計算模塊根據(jù)融合參數(shù)計算當(dāng)前次氯酸鈉投加流量。

圖6　次氯酸鈉投加控制系統(tǒng)結(jié)構(gòu)圖

次氯酸鈉投加算法實現(xiàn)步驟如下：

1）通過采集進(jìn)水流量、出水流量、源水濁度、投礬流量、反沖洗時間、加次氯酸鈉流量、當(dāng)前穩(wěn)定液位等數(shù)據(jù)，進(jìn)行數(shù)據(jù)融合來預(yù)測20 min后的清水庫液位。

當(dāng)設(shè)定次氯酸鈉投加流量變化時，庫氯值發(fā)生變化約滯后25 min，故預(yù)測20 min后的液位可以作為庫氯值滯后變化的超前補償。因進(jìn)水流量、出水流量、源水濁度、投礬流量、反沖洗時間、加次氯酸鈉流量、當(dāng)前穩(wěn)定液位等參數(shù)均對清水庫20 min后的液位有影響，通過對數(shù)據(jù)庫中的數(shù)據(jù)分析比對，采用趨勢擬合法，即將上述參數(shù)的序列值作為自變量，20 min后清水庫液位的序列值作為因變量，建立時間序列的回歸模型，并選擇耦合度較好的擬合函數(shù)。擬合函數(shù)為

0.5+Clflow]/100+dis_now (1)

2）在液位預(yù)測準(zhǔn)確的情況下，通過清水庫液位變化計算次氯酸鈉投加流量，計算公式為

3）分別根據(jù)進(jìn)出水流量變化對液位進(jìn)行預(yù)測，計算次氯酸鈉的投加流量。

進(jìn)水流量對應(yīng)的液位預(yù)測公式為

出水流量對應(yīng)的液位預(yù)測公式為

次氯酸鈉投加流量控制公式為

次氯酸鈉投加算法設(shè)計（式(1)~式(5)）用到的相關(guān)參數(shù)符號定義如表1所示。

表1 算法參數(shù)符號定義

4　算法仿真與調(diào)試

為保證次氯酸鈉投加控制系統(tǒng)的安全性和穩(wěn)定性，將次氯酸鈉投加算法調(diào)試過程分為3個步驟。

步驟2：算法參數(shù)整定合理后，轉(zhuǎn)入全自動控制模式。在全自動模式下，上位機對次氯酸鈉投加流量自動轉(zhuǎn)為變頻泵開度，并控制變頻泵動作，實現(xiàn)次氯酸鈉的自動投加。

步驟3：在全自動控制模式下，采集運行數(shù)據(jù)與日志進(jìn)行分析，對異常數(shù)據(jù)所對應(yīng)各項參數(shù)重新梳理和推導(dǎo)，找出影響結(jié)果的因素，并根據(jù)實際情況，對算法進(jìn)行優(yōu)化。

采用2019年10月15日~10月31日共17天數(shù)據(jù)作為仿真數(shù)據(jù)，進(jìn)行次氯酸鈉投加算法仿真調(diào)試，具體步驟如下：

1） 20 min后的液位預(yù)測

以2019年10月19日的數(shù)據(jù)為例，在式(1)中，取值1 = 0.012，在Clean_W[19]={1,2,3,4,5,6,7,8,9,10,9,8,7,6,5,4,3,2,1}的情況下，預(yù)估液位與實際液位曲線如圖7所示。此時兩者的平均絕對誤差為2.32 cm，由于清水庫液位一般在4 m ~5 m，因此液位預(yù)測結(jié)果在數(shù)據(jù)驗證中較為準(zhǔn)確。

圖7　2019年10月19日預(yù)估液位與實際液位對比曲線

2）通過液位變化計算次氯酸鈉投加流量

式(2)中，要確定2值，需通過現(xiàn)有的次氯酸鈉投加流量、液位預(yù)測值和庫氯來推斷。液位預(yù)測值誤差為0時，次氯酸鈉預(yù)測量最準(zhǔn)確。

圖8　清水庫液位預(yù)測庫氯值與實際庫氯值對比曲線圖

兩者的平均絕對誤差為0.59 mg/L，誤差值較大。由圖8可知：在18:30~22:00時間段內(nèi)，庫氯理論值與實際值較為吻合，但其他時間兩者相差較大。如圖9所示，通過分析進(jìn)水流量、出水流量、源水濁度等對液位預(yù)測有影響的自變量后發(fā)現(xiàn)，進(jìn)出水流量在18:30~22:00較為穩(wěn)定，沒有較大的階梯變化。

圖9　進(jìn)水流量與出水流量對比曲線圖

3）通過進(jìn)出水流量對應(yīng)的液位變化計算次氯酸鈉投加流量

由于通過液位預(yù)測計算次氯酸鈉投加流量與進(jìn)水量、出水量有一定的關(guān)系，故液位預(yù)測值由進(jìn)水流量對應(yīng)的液位預(yù)測值、出水流量對應(yīng)的液位預(yù)測值和穩(wěn)定液位3部分構(gòu)成。

在式(5)中，取3 = 1.5，4 = 0.01，此時進(jìn)水流量對應(yīng)的液位預(yù)測庫氯值與實際庫氯值對比曲線圖如圖10所示，兩者的平均絕對誤差為0.060 mg/L，達(dá)到較為理想的狀態(tài)。

圖10　根據(jù)進(jìn)水流量加權(quán)預(yù)測庫氯值與實際庫氯值對比曲線圖

次氯酸鈉投加算法實現(xiàn)了該水廠次氯酸鈉自動投加控制，通過進(jìn)出水流量、源水濁度等參數(shù)對清水庫液位進(jìn)行趨勢擬合，分別計算進(jìn)水流量、出水流量、穩(wěn)定液位三者所對應(yīng)的預(yù)測液位，并進(jìn)行加權(quán)運算，得出給定庫氯要求下次氯酸鈉的投加流量。通過仿真可知，該算法能夠初步達(dá)到預(yù)想效果，但是還存在一些問題：

1）由于沒有對清水庫的進(jìn)水量進(jìn)行控制，2個清水庫隨機進(jìn)水，導(dǎo)致每個庫的進(jìn)水量隨機；

2）有時出現(xiàn)加次氯酸鈉變頻泵堵塞的情況，導(dǎo)致變頻泵開度最大時流量偏低，無法達(dá)到設(shè)置值；

3）清水庫的次氯酸鈉投加點與庫氯值監(jiān)測點相隔較近，導(dǎo)致在次氯酸鈉與水混合不均勻時就直接檢測；

4）當(dāng)庫氯值穩(wěn)定在0.5 mg/L時，由于進(jìn)水流量增加，算法運算后將加次氯酸鈉流量從30 L/h改變至40 L/h，經(jīng)過一段時間的反應(yīng)，理論上庫氯將穩(wěn)定在0.7 mg/L～0.8 mg/L左右，但穩(wěn)定后庫氯值為1.3 mg/L，與理論推斷不符，推斷變頻泵采集的次氯酸鈉投加流量與實際的次氯酸鈉投加流量可能有一定的出入。

5　自適應(yīng)改進(jìn)算法

為達(dá)到更好的控制效果，本文設(shè)計了自適應(yīng)次氯酸鈉投加算法，對應(yīng)的參數(shù)自適應(yīng)控制系統(tǒng)結(jié)構(gòu)圖如圖11所示。其中，庫氯值穩(wěn)定性判斷主要是利用前60 min的庫氯值，求平均值與給定庫氯值比較得到偏差；權(quán)值產(chǎn)生器根據(jù)庫氯偏差進(jìn)行積分加權(quán)求取次氯酸鈉投加流量的權(quán)重值。經(jīng)測試，權(quán)值產(chǎn)生器中的積分加權(quán)值在0.5~1.0之間取值，選取0.7為最佳。權(quán)值產(chǎn)生器輸出權(quán)值的范圍在0.5~2.0之間。該算法在實際運行中達(dá)到了次氯酸鈉自動投加下庫氯穩(wěn)定可控。

圖11　參數(shù)自適應(yīng)次氯酸鈉投加控制系統(tǒng)結(jié)構(gòu)圖

6　結(jié)語

本文通過數(shù)據(jù)分析次氯酸鈉投加的特征和規(guī)律，設(shè)計了參數(shù)自適應(yīng)次氯酸鈉投加算法。該算法在仿真中達(dá)到較為理想狀態(tài)，在實際運行中可實現(xiàn)庫氯穩(wěn)定可控，從而降低控制成本，減輕員工工作量，同時保證自來水的出廠品質(zhì)。針對不同水廠的典型特征，只要修改一定的參數(shù)，該算法就可以推廣到其他水廠，使加氯環(huán)節(jié)變得簡單、可靠和穩(wěn)定。

[1] SINGER P C. Disinfection by-products in US drinking waters: past, present and future[J]. Water Science & Technology: Water Supply, 2004, 4(1):151-157.

[2] 劉躍輝.某凈水廠藥劑投加系統(tǒng)自動化控制改造研究[D].長春:吉林大學(xué),2017.

[3] 謝崇錦.水廠加藥自動控制系統(tǒng)的設(shè)計與應(yīng)用探討[J].科技風(fēng),2018(18):170.

[4] 王建平,陸翔敏.水廠加藥自動控制系統(tǒng)的設(shè)計與應(yīng)用[J].凈水技術(shù),2016,35(S1):140-143.

[5] 趙士琦.水廠加氯系統(tǒng)自動控制的設(shè)計與實現(xiàn)[J].工業(yè)控制計算機,2012,25(7):116,118.

[6] 封雷.氯氣消毒改造為次氯酸鈉消毒在大型水廠中的應(yīng)用[J].中國科技投資,2017(11):325-326.

Parameter Adaptive NaClO Dosing Algorithm for Waterworks Based on Data Analysis

Peng Yiheng1Zhang Hong1Liu Weibin2

(1.Physics and Electronic Information Engineering of Jianghan University, Wuhan 430056, China 2.Baihezui Waterworks, Wuhan 430040, China)

Based on the actual situation of a waterworks, this paper uses data analysis technology to analyze the characteristics and laws of NaClO dosing in waterworks, and designs an parameter adaptive NaClO dosing algorithm. Based on the waterworks making process, utilizing the historical databases which contains in and out of the water flow, reservoir's water level, source water turbidity and backwash state parameters, Using the methods of trend fitting analysis of the various parameters and the relationship between NaClO addition quantity, to make NaClO dosing controlling convert to reservoir prediction of liquid level controlling. In the process of debugging, according to the actual circumstances of the waterworks, design a set of parameter adaptive algorithm of automatic dosing NaClO. In the simulation, the average absolute error between the theoretical chlorine value and the actual chlorine value of the algorithm is 0.060 mg/L, which reaches the ideal state. In the actual operation, the stable control of the reservoir's free chlorinein can be realized.

automatic control; NaClO dosing; data analysis

江漢大學(xué)產(chǎn)學(xué)研用橫向課題（H2019120308175693）

彭一恒，男，1996年生，碩士研究生，主要研究方向：智能機器人。

章紅（通信作者），女，1969年生，教授，主要研究方向：計算機控制系統(tǒng)，網(wǎng)絡(luò)化控制系統(tǒng)及其故障診斷與容錯控制等。E-mail: ujjxa@163.com

TP273.2

A

1674-2605(2020)05-0005-07

10.3969/j.issn.1674-2605.2020.05.005