以概率分析方法確定北方地區(qū)某水廠的優(yōu)化運(yùn)行方案

徐勇鵬¹，崔福義¹，劉國(guó)平²
（1.哈爾濱工業(yè)大學(xué)市政環(huán)境工程學(xué)院，黑龍江 哈爾濱，150090；2. 哈爾濱紹和供水有限公司，黑龍江 哈爾濱，150001）

摘 要：我國(guó)東北地區(qū)水質(zhì)受氣候影響，處理難度較大，研究北方水廠的優(yōu)化運(yùn)行規(guī)律能夠保證水質(zhì)與水量，并且降低運(yùn)行費(fèi)用。采用概率分析方法分析整理了東北地區(qū)S水廠的連續(xù)三年的生產(chǎn)運(yùn)行數(shù)據(jù)，探討了水處理系統(tǒng)優(yōu)化運(yùn)行規(guī)律，并繪制出該水廠混凝工藝運(yùn)行工況圖和最經(jīng)濟(jì)沉后水濁度指示圖。根據(jù)兩圖可以確定使該水廠優(yōu)化運(yùn)行的沉后水生產(chǎn)控制目標(biāo)值，進(jìn)而確定經(jīng)濟(jì)的混凝劑投加量。按相應(yīng)的最經(jīng)濟(jì)沉后水濁度優(yōu)化水廠運(yùn)行工況，體現(xiàn)了明顯的經(jīng)濟(jì)效益。
關(guān)鍵詞：優(yōu)化運(yùn)行，運(yùn)行費(fèi)用，最經(jīng)濟(jì)沉后水濁度，最佳沉后水濁度 設(shè)想

　　我國(guó)東北地區(qū)受氣候影響，水質(zhì)條件隨季節(jié)變化很大。每年冬季有4～5個(gè)月的冰封期，原水在這一時(shí)期呈現(xiàn)低溫低濁特性；夏季汛期時(shí)，原水濁度較高。近些年松花江水受到了嚴(yán)重的污染，水中的有機(jī)物含量高。這種北方河流水體，水質(zhì)條件隨季節(jié)變化很大，不同季節(jié)的水質(zhì)具有不同的混凝和絮凝體分離特性，給凈水處理系統(tǒng)增加較大難度^[[xv]]。隨著飲用水水質(zhì)標(biāo)準(zhǔn)的提高，該種水質(zhì)的處理難度又有增加。為了使濾后出水達(dá)到要求，多數(shù)供水企業(yè)強(qiáng)化混凝工藝盡力將沉后水濁度控制在2NTU以下，有些水廠甚至要求低于1NTU。較低的沉后水濁度勢(shì)必導(dǎo)致混凝投入過大，由此可能導(dǎo)致給水處理系統(tǒng)的運(yùn)行費(fèi)用增加^[[xvi]]。
　　 由于保證水質(zhì)與水量、降低運(yùn)行成本始終是貫穿水廠運(yùn)行的目標(biāo)，因此在這種情況下，用優(yōu)化的思想指導(dǎo)水廠運(yùn)行是具有重要意義的^[[xvii]]。本文結(jié)合哈爾濱S水廠多年的生產(chǎn)數(shù)據(jù)，采用概率分析的方法，分析水處理系統(tǒng)優(yōu)化運(yùn)行規(guī)律，尋找最經(jīng)濟(jì)沉后水濁度^[[xviii]]，并用以指導(dǎo)生產(chǎn)。

1 水廠概況和分析的基本方案

　　哈爾濱S水凈水能力為225000m³/d。工藝流程見圖1。S水廠以地表水為水源，水質(zhì)運(yùn)行目標(biāo)要求出廠水濁度不高于1NTU，COD_Mn不超過3mg/L。
　　 收集S水廠在三年的生產(chǎn)運(yùn)行數(shù)據(jù)2萬多組，并根據(jù)原水水質(zhì)狀況，將濁度和水溫分別劃分成不同范圍，分成幾類水質(zhì)，分析處理S水廠的生產(chǎn)記錄數(shù)據(jù)，探討S水廠生產(chǎn)工藝各單元環(huán)節(jié)的運(yùn)行規(guī)律，并尋找適合S水廠處理系統(tǒng)優(yōu)化優(yōu)化運(yùn)行的有效途徑。

圖1 S水廠工藝流程圖

　　圖2是水廠原水水質(zhì)隨著季節(jié)的變化曲線。從圖中可以看出從每年的4月中旬到10月中旬水溫在6℃以上，濁度在20～300NTU之間，pH值在7左右，COD_Mn為5～10mg/L，平均為8mg/L，將此期間水質(zhì)稱為常溫水質(zhì)。從每年的10月下旬到第二年4月上旬松花江處于冰封期，水溫在0～5℃，濁度在5～30NTU，原水的平均COD_Mn為5mg/L，是典型的低溫低濁水質(zhì)。

圖2 1997年原水水質(zhì)

　　從圖中發(fā)現(xiàn)，原水水質(zhì)的諸多參數(shù)中，原水溫度和濁度的變化范圍較大，許多文獻(xiàn)中也提到過，這兩個(gè)參數(shù)對(duì)給水處理各工藝的處理效率有較大影響^[[xix],[xx]]，為了方便研究雜亂的生產(chǎn)數(shù)據(jù)，將溫度和濁度各分成幾個(gè)階段。溫度分成（0～5）、（5～13）、（13～21）、（21～28）4個(gè)階段，分別用字母A、B、C、D表示；濁度分成（0～10）、（10～50）、（50～100）、（100以上）4個(gè)階段，分別用數(shù)字1、2、3、4表示。交叉組合，一共分成16組情況，三年內(nèi)每組原水水質(zhì)出現(xiàn)的頻率列在表1中。在16種情況中出現(xiàn)過12種，其中出現(xiàn)概率最大的是A1和A2組，三年內(nèi)總共出現(xiàn)的概率超過37%。取出現(xiàn)概率大于2.0%的情況作為分析研究?jī)?nèi)容，按照在表中的的位置，分別定義為A1、A2、B2、B3、B4、C2、C3、C4、D3、D4等10組，A系列表示低溫水質(zhì)，B、C、D代表常溫水質(zhì)。

不同水質(zhì)出現(xiàn)的概率　　　　　 表1

2 單元優(yōu)化運(yùn)行分析

2.1 混凝工藝
　　S水廠投加混凝劑硫酸鋁（AS），并用助凝劑活化硅酸（Si）強(qiáng)化混凝。由于混凝劑和助凝劑的質(zhì)量投加量的配合比例（以下簡(jiǎn)稱配比）對(duì)混凝效果有影響，適當(dāng)?shù)呐浔扔锌赡芙档?0%左右的混凝劑消耗，并降低沉后水濁度1/3～1/2^[[xxi]]，因此以下著重討論配比和沉后水濁度的關(guān)系。
　　 由于生產(chǎn)數(shù)據(jù)比較分散，平均結(jié)果可能掩蓋某些因素。而在實(shí)際生產(chǎn)中經(jīng)常出現(xiàn)配比相同，沉后水濁度不同的現(xiàn)象。在分析隨機(jī)性較大的生產(chǎn)數(shù)據(jù)時(shí)，應(yīng)考慮在相同配比下不同沉后水濁度出現(xiàn)的概率。將配比分成5段：（5～15）、（15～25）、（25～35）、（35～45）、（45～55）；把所有出現(xiàn)過的沉后水濁度分成7個(gè)區(qū)間：（0，1）、（1，1.5） 、（1.5，2）、（2，2.5） 、（2.5，3） 、（3，3.5）、（3.5，4），分別計(jì)算在相同配比下，沉后水濁度在各個(gè)區(qū)間內(nèi)出現(xiàn)的概率，結(jié)果見圖3。

a （C3組）　　　　　　　　　　　　　　 b （A1與A2）

圖3 各種同配比下不同沉后水濁度出現(xiàn)的概率

　　圖3a是以C3為例的原水水質(zhì)條件下，沉后水濁度在不同區(qū)間內(nèi)出現(xiàn)概率的曲線。。從圖中可以看出，曲線的最大值和配比成反比，配比增加，峰值向數(shù)值軸方向移動(dòng)，即混凝劑和助凝劑的配合比例在（5～15）內(nèi)，沉后水濁度在（2.5～3）NTU的區(qū)間內(nèi)出現(xiàn)的概率最大，當(dāng)配比增加到（35～45）時(shí)，出現(xiàn)概率最大的沉后水濁度降低至（1.5～2.0）NTU，隨著配比值的增加，出現(xiàn)概率最大點(diǎn)的沉后水濁度繼續(xù)降低。同樣B、C、D系列等統(tǒng)計(jì)結(jié)果與圖4a相似。圖3b為低溫原水條件下的情況。同圖3a一樣，每條曲線都有一個(gè)最大值，該點(diǎn)對(duì)應(yīng)著該投加配合比例下，出現(xiàn)概率最大的沉后水濁度。與圖3a不同的是：1）最高概率的沉后水濁度出現(xiàn)在（2.5～3.5）NTU內(nèi)，而且隨著配比增加，最大值的位置在此區(qū)間內(nèi)有向遠(yuǎn)離數(shù)值軸方向移動(dòng)的微弱趨勢(shì)，即隨著配比的變化，出現(xiàn)概率最大的濁度有升高的趨勢(shì)；2）對(duì)于常溫水質(zhì)：B、C、D系列水質(zhì)，混凝藥劑配比在25以上，低于2NTU的沉后水濁度出現(xiàn)概率較大；而對(duì)于低溫低濁水質(zhì)，混凝藥劑配比在15～25之間，較低的沉后水濁度出現(xiàn)概率大。
　　同樣，可以得出其它各組水質(zhì)的配比與沉后水濁度的出現(xiàn)概率之間相互關(guān)系的曲線，并將各組數(shù)據(jù)匯總在一起，同時(shí)計(jì)算各沉后水濁度區(qū)間混凝投加量與配比的平均值，結(jié)果見圖4和圖5。在圖4中分類軸（x）和數(shù)值軸（y）的網(wǎng)格線圍成的每一個(gè)網(wǎng)格代表一個(gè)工況，圖中的數(shù)據(jù)點(diǎn)所處的網(wǎng)格是生產(chǎn)上出現(xiàn)概率最大的工況區(qū)域。圖中曲線表明，為了降低沉后水濁度，宜采用高配比（25～45）進(jìn)行混凝，也可以根據(jù)圖4中的指示，依據(jù)沉后水濁度目標(biāo)值，確定獲得該沉后水濁度可能性最大的混凝劑和助凝劑投加量的配比值，然后根據(jù)圖5中（以C3水質(zhì)為例）混凝劑投加量指示圖，確定該配比下的能夠獲得沉后水濁度目標(biāo)值可能性最大的混凝劑投加量。例如，沉淀出水目標(biāo)值為2NTU，在每年的6月份水溫一般為15～20℃，濁度一般為60～100NTU，是C3類原水水質(zhì)，根據(jù)圖4中的C3曲線可以確定沉后水濁度2NTU時(shí)對(duì)應(yīng)的配比應(yīng)為25～35。并結(jié)合圖5中2NTU的曲線確定混凝劑量為40～45mg/L，也可知助凝劑量為1.4～1.5mg/L。

圖4 配比和出現(xiàn)概率最高的沉后水濁度情況圖

圖5 混凝劑投加量指示圖（C3）

2.2沉淀和過濾工藝
　　綜合可變費(fèi)用中包含沉淀池排泥費(fèi)用和反沖洗水費(fèi)用^[2]，需要合理確定排泥周期和濾池周期。
　　文獻(xiàn)^[[xxii]]提到沉淀池排泥水費(fèi)在給水處理系統(tǒng)綜合可變費(fèi)用中所占的比例較小，可以采用簡(jiǎn)單的時(shí)間模式控制沉淀池排泥周期和排泥歷時(shí)。通過在不同水質(zhì)期間進(jìn)行排泥試驗(yàn)，結(jié)果表明，即低溫低濁水質(zhì)期間，將排泥周期定為3h，排泥歷時(shí)定為1.5min，其它情況下，將排泥周期定為2h，排泥歷時(shí)定為2min，即能保證沉淀出水水質(zhì)，又比較經(jīng)濟(jì)節(jié)水。
　　由于取自松花江的原水水質(zhì)隨季節(jié)性的變化較大，濾池的運(yùn)行周期不能固定不變，需要根據(jù)不同水質(zhì)條件確定與其相應(yīng)的周期。圖6是三年內(nèi)10組不同原水水質(zhì)條件下，每個(gè)沉后水濁度區(qū)間內(nèi)濾池周期平均值的統(tǒng)計(jì)結(jié)果。

圖6 過濾周期和沉后水濁度的關(guān)系

　　從圖中可知，過濾周期隨著沉后水濁度的升高而降低；每種水質(zhì)的每個(gè)沉后水區(qū)間內(nèi)對(duì)應(yīng)著不同的過濾周期，選取該值作為下一步計(jì)算經(jīng)濟(jì)費(fèi)用的參考值。反沖洗時(shí)間根據(jù)水廠數(shù)據(jù)，低溫低濁期間和其它的水質(zhì)期間處理工藝的運(yùn)行參數(shù)不同，因此需要分別考慮。根據(jù)兩種情況下水廠生產(chǎn)記錄，水廠反沖洗排水濁度低于10NTU作為目標(biāo)值，確定常溫水質(zhì)期間反沖洗參數(shù)為，單獨(dú)氣洗2min，氣水混合洗4 min，水洗4 min；低溫低濁期間的沖洗時(shí)間為氣沖2min，氣水混合沖4min和水沖5min。

3 系統(tǒng)優(yōu)化運(yùn)行分析

3.1 綜合可變費(fèi)用和最經(jīng)濟(jì)沉后水濁度
　　那么，如何確定沉后水濁度的生產(chǎn)控制目標(biāo)值呢？下面將結(jié)合給水處理系統(tǒng)運(yùn)行過程中的可變運(yùn)行費(fèi)用，討論使綜合可變運(yùn)行費(fèi)用最低的最經(jīng)濟(jì)沉后水濁度的變化確定方法。
　　綜合可變費(fèi)用是混凝劑費(fèi)用、助凝劑費(fèi)用、排泥水費(fèi)用和反沖洗水費(fèi)用的綜合，S水廠的混凝劑是硫酸鋁，單價(jià)為630元/噸，助凝劑活化硅酸的單價(jià)為290元/噸，排泥水單價(jià)＝水資源費(fèi)（0.3元/m³）+藥劑單耗，反沖洗水費(fèi)的單價(jià)為1.5元/m³。根據(jù)原水分類方法，將10組水質(zhì)的各項(xiàng)可變費(fèi)用計(jì)算并繪圖（如圖7，以C3水質(zhì)為例），發(fā)現(xiàn)：在綜合可變費(fèi)用中，混凝劑的費(fèi)用所占的比例最大；助凝劑的費(fèi)用和反沖洗水費(fèi)相當(dāng)，各占10%左右；排泥水費(fèi)用比例最小，只有4%左右。因此，降低運(yùn)行費(fèi)用的關(guān)鍵是降低混凝劑費(fèi)用，排泥水費(fèi)所占的比重最小，可以采用簡(jiǎn)單的方式確定排泥周期和排泥歷時(shí)^[8]。

圖7 可變費(fèi)用與沉后水濁度的關(guān)系曲線（C3）

　　曲線上使綜合可變費(fèi)用最小的沉后水濁度稱為最經(jīng)濟(jì)沉后水濁度，10組水質(zhì)的最經(jīng)濟(jì)沉后水濁度列于表，統(tǒng)計(jì)結(jié)果列于表2中。由于水溫是影響混凝沉淀效果的重要因素，對(duì)比表2中相同水溫下不同原水濁度對(duì)應(yīng)的最經(jīng)濟(jì)沉后水濁度，發(fā)現(xiàn)1）在低溫條件下（A系列），隨著原水濁度升高，最經(jīng)濟(jì)沉后水濁度降低，這是由于A1水質(zhì)是較難處理的水質(zhì)；2）B和C系列水質(zhì)隨著原水濁度的升高最經(jīng)濟(jì)沉后水濁度也升高，其原因是，原水濁度高的藥劑量大，因此最經(jīng)濟(jì)沉后水濁度升高；3）D系列水質(zhì)經(jīng)常出現(xiàn)在夏季洪水期，運(yùn)行工況很不穩(wěn)定，應(yīng)該視具體情況而定，因此，表中的數(shù)據(jù)不能說明其規(guī)律。
　　 經(jīng)過上述分析，可知，在實(shí)際生產(chǎn)中，可以確定最經(jīng)濟(jì)沉后水濁度為沉后水濁度目標(biāo)值，并通過圖4和5就可以制定混凝劑和助凝劑的投加量。

3.2 有機(jī)物指標(biāo)對(duì)最經(jīng)濟(jì)沉后水濁度的約束
　　最經(jīng)濟(jì)沉后水濁度可以使現(xiàn)行水處理系統(tǒng)的可變費(fèi)用最低，但從圖2看，生活污水和工業(yè)廢水的大量排放使水源的污染非常嚴(yán)重，COD_Mn平均在8mg/L左右，最高可至14mg/L，對(duì)給水處理系統(tǒng)的優(yōu)化運(yùn)行條件產(chǎn)生影響，不僅要使濾后水濁度低于1NTU，而且有機(jī)物指標(biāo)耗氧量也要低于3mg/L。嚴(yán)重污染在一定程度上會(huì)導(dǎo)致濾后水濁度（取過濾周期內(nèi)的平均值）的增加和隨之帶來的有機(jī)物泄漏^[[xxiii]] ，因此，給水處理系統(tǒng)優(yōu)化運(yùn)行受有機(jī)物約束，滿足約束條件的沉后水濁度為最佳沉后水濁度。圖8是10中水質(zhì)條件下的濾后水COD_Mn和沉后水濁度的關(guān)系圖。為了保證濾后水COD_Mn低于3mg/L，要求沉后水濁度不超過一個(gè)限值，將有機(jī)物約束條件計(jì)算結(jié)果列于表2中。從結(jié)果中發(fā)現(xiàn)，最經(jīng)濟(jì)沉后水濁度皆小于有機(jī)物指標(biāo)約束的限值，上述最佳沉后水濁度在數(shù)值上等于最經(jīng)濟(jì)沉后水濁度。

圖8 平均濾后水CODMn和沉后水濁度的關(guān)系

3.3 沉后水濁度的生產(chǎn)控制目標(biāo)
　　前面的分析結(jié)果表明，在不同原水水質(zhì)期間，使運(yùn)行費(fèi)用最低的最經(jīng)濟(jì)沉后水濁度不同。將表2中的最經(jīng)濟(jì)沉后水濁度繪制圖9中，在生產(chǎn)實(shí)踐中可以根據(jù)圖中的指示確定沉后水濁度生產(chǎn)控制目標(biāo)。
　　圖9a中的曲面是S水廠的最經(jīng)濟(jì)沉后水濁度指示圖，圖9b是滿足有機(jī)物約束條件的沉后水濁度范圍。從兩圖中可知圖9b中的濁度值均較圖9a中的濁度值偏高，圖9a中的最經(jīng)濟(jì)沉后水濁度在數(shù)值上等于滿足有機(jī)物約束的最佳沉后水濁度。

最經(jīng)濟(jì)（最佳）沉后水濁度　　　　　　　　　　　　 表2

　　同時(shí)，從兩圖中可知，隨著原水水溫和濁度的變化，最經(jīng)濟(jì)沉后水的變化曲面是一個(gè)中間低、四周高的下凹形的曲面，當(dāng)原水濁度和溫度較低或者較高時(shí)所對(duì)應(yīng)的最經(jīng)濟(jì)沉后水濁度處于曲面的高處，此范圍的原水水質(zhì)的混凝效果一般，最經(jīng)濟(jì)的沉后水濁度較高。當(dāng)原水溫度較高和濁度適中的情況下，最經(jīng)濟(jì)沉后水濁度處在曲面的低凹處，此時(shí)的原水條件混凝的效果較好，最經(jīng)濟(jì)沉后水濁度相對(duì)偏低。對(duì)于地處北方地區(qū)的S水廠，原水水質(zhì)變化較大，不同時(shí)期的沉后水濁度目標(biāo)值不同，為使水處理系統(tǒng)運(yùn)行最接近優(yōu)化狀態(tài)，可以參考圖9指示，制定生產(chǎn)運(yùn)行中的沉后水濁度目標(biāo)值。然后根據(jù)圖4和圖5中的指示，確定混凝劑和助凝劑的投加量。

　　圖9是在前面介紹的原水分類的基礎(chǔ)上，將有機(jī)物指標(biāo)作為約束條件最終獲得的，給水處理系統(tǒng)優(yōu)化運(yùn)行還受到比較多的條件約束，可以參考有機(jī)物約束的討論方法確定最經(jīng)濟(jì)沉后水濁度。
　　以上原水的分類方法比較適合北方地區(qū)季節(jié)變化較大的河流水體，如果對(duì)于湖泊水體或者南方地區(qū)溫度變化不明顯的水體，則需要在原水分類方法上找到影響因素較大的水質(zhì)參數(shù)進(jìn)行分類，繪制出最經(jīng)濟(jì)沉后水濁度圖。
　　在對(duì)三年的生產(chǎn)運(yùn)行數(shù)據(jù)的分析之后，在S水廠調(diào)整運(yùn)行條件以最經(jīng)濟(jì)沉后水濁度運(yùn)行1年。圖11是2000年優(yōu)化運(yùn)行結(jié)果與前三年平均運(yùn)行結(jié)果的對(duì)比圖。從圖中可知，除在6、7、8、9和12月份，2000年的綜合可變費(fèi)用稍高一些，在一年的大部分時(shí)間都較前三年的平均水平每月每噸水大約可以節(jié)省0.01元；一年可以較前三年節(jié)省約5%的綜合可變費(fèi)用。

圖10 經(jīng)濟(jì)比較

4 結(jié)論

　　1.用概率分析方法討論了混凝劑和助凝劑投加量配合比例和沉后水濁度的對(duì)應(yīng)關(guān)系，得出常溫水質(zhì)混凝工藝的優(yōu)化運(yùn)行工況圖，可以依據(jù)沉后水濁度目標(biāo)值確定混凝劑和助凝劑的投加量，以及低溫低濁水質(zhì)混凝工藝優(yōu)化運(yùn)行的混凝劑和助凝劑的投加量確定方法。
　　2.水廠的綜合可變費(fèi)用60%以上是混凝劑費(fèi)用，沉淀池排泥費(fèi)用最少，只占4%左右，控制混凝是優(yōu)化運(yùn)行的關(guān)鍵。在不同水質(zhì)條件下，綜合可變費(fèi)用隨著沉后水濁度變化的趨勢(shì)不同，根據(jù)二者相互關(guān)系分別歸納出使綜合可變費(fèi)用最低的最經(jīng)濟(jì)沉后水濁度的范圍。
　　3.符合有機(jī)物約束的最佳沉后水濁度在數(shù)值上等于是最經(jīng)濟(jì)沉后水濁度。根據(jù)不同的原水水質(zhì)分組，按相應(yīng)的最經(jīng)濟(jì)沉后水濁度優(yōu)化水廠運(yùn)行工況，體現(xiàn)了明顯的經(jīng)濟(jì)效益。

參考文獻(xiàn)

1.霍明昕,劉馨遠(yuǎn). 低溫低濁水質(zhì)特性的分析. 中國(guó)給水排水, 1998, 14（6）: 33~34.
2 石明巖,崔福義,趙天慧. 給水處理系統(tǒng)高效經(jīng)濟(jì)運(yùn)行研究之二——相關(guān)關(guān)系的試驗(yàn)研究. 哈爾濱建筑大學(xué)學(xué)報(bào). 1999,32（6）:9~11.
3. 田一梅等.水處理系統(tǒng)優(yōu)化運(yùn)行. 中國(guó)給水排水. 1999,15（5）:5~9.
4. 趙天慧.給水處理系統(tǒng)優(yōu)化運(yùn)行的初步研究. 哈爾濱建筑大學(xué)碩士論文. 1996,（6）:14~16
5. 孟凡良,崔福義. 低溫低濁地表水處理技術(shù)的探討. 哈爾濱商業(yè)大學(xué)學(xué)報(bào). 2003, 19（4）:188~190
6. 王云波,劉金河,譚萬春等.PDA2000監(jiān)測(cè)低溫低濁水混凝試驗(yàn)研究. 長(zhǎng)沙交通學(xué)院學(xué)報(bào). 2003,19（1）:54～57
7. David J. Pernitsky. Drinking water coagulation with polyaluminum coagulants mechanisms and selection guidelines. Subintted to the Graduate School of the University of Massachusetts Amherst in partial fulfillment of the requriments for the degree of Doctor of Philosophy, 2001,9: 44~55.
8. 崔福義，徐勇鵬，趙志偉等. 東北地區(qū)水廠沉淀池排泥控制模式探討. 給水排水.2003, 29（5）:10~13.
9. David Hall, Caroline S.B. & Fitzpatrick. Suspension on centration profiles during rapid gravity filter backwashing. Chemical Engineering Journal. 2000, 80:197~201.