孫艷，張逢，胡洪營，2*，牛璋彬

1．清華大學環(huán)境學院，環(huán)境模擬與污染控制國家重點聯(lián)合實驗室，北京 100084

2．清華大學深圳研究生院，國家環(huán)境保護環(huán)境微生物利用與安全控制重點實驗室，廣東深圳 518055

3．住房城鄉(xiāng)建設部城鎮(zhèn)水務管理辦公室，北京 100835

天津市污水處理廠進水水質特征的統(tǒng)計學分析

孫艷1，張逢1，胡洪營1，2*，牛璋彬3

1．清華大學環(huán)境學院，環(huán)境模擬與污染控制國家重點聯(lián)合實驗室，北京 100084

2．清華大學深圳研究生院，國家環(huán)境保護環(huán)境微生物利用與安全控制重點實驗室，廣東深圳 518055

3．住房城鄉(xiāng)建設部城鎮(zhèn)水務管理辦公室，北京 100835

以天津市24座城市污水處理廠實際運行數(shù)據(jù)為基礎，系統(tǒng)分析了進水水質特征及有機物、氮、磷和懸浮物之間的概率分布及相關關系。結果表明，天津市污水中BOD5、CODCr、SS、NH3|N、TN和TP全年濃度均呈正偏態(tài)分布，其月中間值分布范圍分別為95～140、280～370、130～155、15～30、20～35和3～5 mg/L。進水各水質指標間存在較好的一元線性關系，其中BOD5與TP相關性最為顯著，R2為0．983。進水BOD5/CODCr分布在0．4～0．6的概率為39．4%，分布在0．2～0．4的概率為50．2%，說明天津市污水的可生化性較好。BOD5/TN＜4的概率為54．4%，表明多數(shù)情況下進水反硝化碳源不足。BOD5/TP＞20的累積概率為77．9%，表明進水可以滿足生物除磷的需求。TN/TP分布在5～15的概率為73．6%，平均值為9．7，且TN/TP＞5的概率為81．6%，表明進水能夠滿足微生物生長對氮、磷的需求。

污水處理廠;進水水質;水質特征;統(tǒng)計分析

開發(fā)利用非常規(guī)水資源是解決天津水資源短缺的重要手段［4－5］。城市污水作為非常規(guī)水資源的典型代表，其水量穩(wěn)定、供給可靠，對城市污水進行有效處理，達到目標用途要求的水質標準，使之成為多種用途的水資源，是開源節(jié)流、減輕水體污染、改善生態(tài)環(huán)境、解決城市缺水問題的有效途徑之一［6－8］。城市污水水質特征是污水處理廠設計和建設的基本參數(shù)，污水處理廠進水水質的測定分析對污水處理廠的工藝設計和運行管理具有重要意義，掌握污水水質特征能夠為有效提高水處理效果及再生利用提供科學依據(jù)［9－11］。筆者以天津市24座污水處理廠為研究對象，系統(tǒng)分析了進廠原水水質特征，主要包括水質指標BOD5、CODCr、SS、NH3|N、TN和TP的變化規(guī)律和概率統(tǒng)計分析，以及有機物、氮、磷和懸浮物之間的相關關系，以期為天津市污水再生利用，緩解水資源短缺提供可靠數(shù)據(jù)支持。

1 研究方法

以天津市24座污水處理廠為研究對象，以2012年全年實際進廠原水水質數(shù)據(jù)為基礎，主要分析了BOD5、CODCr、SS、NH3|N、TN和TP等的變化規(guī)律，各指標間的相關關系及概率分布。水質分析數(shù)據(jù)來源于住房城鄉(xiāng)建設部城鎮(zhèn)污水處理管理信息系統(tǒng)，統(tǒng)計分析采用SPSS20．0軟件，作圖分析采用Origin9軟件。

2 結果與討論

2.1 進水水質分布特性分析

天津市2012年全年污水處理廠進廠原水主要水質指標BOD5、CODCr、SS、NH3|N、TN和TP的統(tǒng)計分析結果和正態(tài)性檢驗如表1所示。

表1 污水水質統(tǒng)計分析和正態(tài)性檢驗Table 1 Statistical analysis of influent quality characteristics and normality test

根據(jù)表1給出的Kolmogorov－Smirnov統(tǒng)計量、Shapiro－Wilk統(tǒng)計量、偏度和峰度結果，可以對各項水質指標的分布進行正態(tài)性檢驗，當樣本容量N＜ 1 000時，以Shapiro－Wilk檢驗為準。經(jīng)過Shapiro－Wilk檢驗，各項水質指標顯著性水平(Sig．，significance level)均小于0．05，表明各項指標不服從正態(tài)分布。從偏度和峰度系數(shù)可以看出，各污水水質指標的偏度系數(shù)和峰度系數(shù)均大于0，由此可判斷各項水質指標數(shù)據(jù)分布均呈正偏態(tài)分布。因此，在數(shù)據(jù)統(tǒng)計時，對各水質指標分布的平均值、中間值及最頻值分別進行了統(tǒng)計分析。表2列出了各月濃度中間值。

表2 污水水質指標月中間值Table 2 Themonthlymedian of pollution indicators mg/L

2.1.1 進水BOD5分布特性

在三位負責人的帶領下，我們身著工作服，參觀了從上鋁卷到電化學處理、涂布和干燥、在線分切、包裝，直至儲運、物流的完整生產(chǎn)流程。我們看到，在監(jiān)控室里高速生產(chǎn)線的實時數(shù)據(jù)不斷閃現(xiàn)，所有數(shù)據(jù)可以保留2年；生產(chǎn)全過程參數(shù)自動檢測，全球同步；安裝于生產(chǎn)線上的在線檢測系統(tǒng)，可檢測到50微米缺陷，自動切走缺陷部分；采用環(huán)保包裝方式，自動打包。此外，我們了解到，其關鍵材料由總部統(tǒng)一進口；新開發(fā)產(chǎn)品，需經(jīng)過12個月嚴格測試，才能投入正式生產(chǎn)……或許正因為執(zhí)行著最嚴苛的生產(chǎn)要求，無錫工廠才敢于做出最高質量標準的莊重承諾。

進水BOD5濃度逐月變化規(guī)律和概率分布如圖1所示，各月濃度中間值如表2所示。從圖1(a)可以看出，BOD5分布在20～360 mg/L，各月濃度平均值均大于中間值，各月中間值分布在95～140mg/L，其中8月濃度最低，3月濃度最高。由表2可知，中間值最小值為98．0mg/L，最大值為136．3mg/L。從圖1(b)可以看出，進水 BOD5平均值為 125．2 mg/L，中間值為115．0 mg/L。概率分布較高的濃度范圍在50～150 mg/L，其累積概率為75．8%，達到累積頻率95%時，BOD5為253．0 mg/L。

圖1 進水BOD5逐月變化及概率分布(2012年1月—2012年12月)Fig．1 Themonthly change and probability distribution of influent BOD5(January 2012－December 2012)

2.1.2 進水CODCr分布特性

進水CODCr逐月變化規(guī)律和概率分布如圖2所示，各月濃度中間值如表2所示。從圖2(a)可以看出，CODCr分布在30～900 mg/L，1—7月各月濃度平均值大于中間值，8—12月平均值小于中間值。各月中間值分布在280～370 mg/L，其中8月濃度最低，3月濃度最高。由表2可知，中間值最小值為283．0 mg/L;最大值為364．7 mg/L。從圖2(b)可以看出，進水CODCr平均值為334．5 mg/L，中間值為330．5 mg/L。概率分布較高的濃度范圍在150～400 mg/L，其累積概率為77．1%，達到累積頻率95%時，CODCr為575．0 mg/L。

圖2 進水CODCr逐月變化及概率分布(2012年1月—2012年12月)Fig．2 Themonthly change and probability distribution of influent CODCr(January 2012－December 2012)

2.1.3 進水SS分布特性

進水SS濃度逐月變化規(guī)律和概率分布如圖3所示，各月濃度中間值如表2所示。從圖3(a)可以看出，SS濃度分布在20～540 mg/L，各月SS濃度平均值均大于中間值，中間值分布在130～155 mg/L，隨季節(jié)變化不明顯。由表2可知，中間值最小值為134．3 mg/L;最大值為152．0 mg/L。從圖3(b)可以看出，進水SS濃度平均值為162．6 mg/L，中間值為142．0 mg/L。概率分布較高的濃度范圍在40～280 mg/L，其累積概率為89．3%，達到累積頻率95%時，SS濃度為323．0 mg/L。

圖3 進水SS濃度逐月變化及概率分布(2012年1月—2012年12月)Fig．3 Themonthly change and probability distribution of influent SS(January 2012－December 2012)

2.1.4 進水NH3|N分布特性

進水NH3|N濃度逐月變化規(guī)律和概率分布如圖4所示，各月濃度中間值如表2所示。從圖4(a)可以看出，NH3|N濃度分布在5～90 mg/L，各月濃度平均值均大于中間值，中間值分布在15～30 mg/L，其中8月濃度最低，5月濃度最高。由表2可知，中間值最小值為 15．6 mg/L;最大值為 28．3 mg/L。從圖4(b)可以看出，進水NH3|N平均值為26．0 mg/L，中間值為24．0 mg/L。概率分布較高的濃度范圍在5～50 mg/L，其累積概率為93．9%。

圖4 進水NH3|N濃度逐月變化及概率分布(2012年1月—2012年12月)Fig．4 Themonthly change and probability distribution of influent NH3|N(January 2012－December 2012)

2.1.5 進水TN分布特性

進水TN濃度逐月變化規(guī)律和概率分布如圖5所示，各月濃度中間值如表2所示。從圖5(a)可以看出，TN濃度分布在3～85 mg/L，各月中間值分布在20～35 mg/L，除2、4、5月以外，各月濃度中間值均小于平均值。其中8月TN濃度最低，3月濃度最高。由表2可知，中間值最小值為23．5 mg/L;最大值為33．6 mg/L。從圖5(b)可以看出，進水TN平均值為33．6 mg/L，中間值為30．0 mg/L。概率分布較高的濃度范圍在 6～66 mg/L，其累積概率為97．1%。

圖5 進水TN濃度逐月變化及概率分布(2012年1月—2012年12月)Fig．5 Themonthly change and probability distribution of influent TN(January 2012－December 2012)

2.1.6 進水TP分布特性

進水TP濃度逐月變化規(guī)律和概率分布如圖6所示，各月濃度中間值如表2所示。從圖6(a)可以看出，TP濃度分布在0．5～20 mg/L，各月中間值分布在3～5 mg/L，各月濃度中間值均小于平均值。其中8月TP濃度最低，4月和5月濃度最高。由表2可知，中間值最小值為3．4 mg/L;最大值為4．2 mg/L。從圖6(b)可以看出，進水TN平均值為4．2 mg/L，中間值為3．8 mg/L。概率分布較高的濃度范圍0～7 mg/L，其累積概率為95．7%。

2.2 進水水質指標相關性分析

污水中污染物的組成非常復雜，分析各污染物間相關性，了解各污染指標間的相關程度，可以減少監(jiān)測次數(shù)與項目，還可保證監(jiān)測數(shù)據(jù)的真實可靠性［12－13］。針對污水處理廠2012年全年進水水質指標BOD5、CODCr、SS、NH3|N、TN和TP，采用最小二乘法對監(jiān)測數(shù)據(jù)進行回歸計算，確定各指標間的一元線性關系是否存在，建立回歸方程，得到各污染指標之間的相關關系。

圖6 進水TP濃度逐月變化及概率分布(2012年1月—2012年12月)Fig．6 Themonthly change and probability distribution of influent TP(January 2012－December 2012)

各指標間的回歸方程及相關系數(shù)(R2)如表3所示。由表3可見，各項指標間線性相關關系均較好，除 CODCr與 TN的相關關系相對弱些(R2= 0．893)之外，其他指標間 R2均大于0．900，其中BOD5與TP相關性最為顯著，R2為0．983。由此可見，天津市污水中污染物組成相對穩(wěn)定，無明顯劇烈變化，監(jiān)測結果也沒有明顯突變，各水質指標間的相關關系也較為穩(wěn)定和顯著。

表3 污水水質指標間相關和回歸分析Table 3 The correlation and regression analysis of pollution indicators

2.3 進水營養(yǎng)物質的比例關系

2.3.1 進水BOD5/CODCr

污水處理中常用BOD5/CODCr來衡量污水的可生化性。當BOD5/CODCr為0．4～0．6時，污水可生化性較好，適應于生物處理，當BOD5/CODCr為0．2～0．4時，表明污水中存在難生物降解性污染物，當BOD5/CODCr＜0．1時，可生化性差，必須進行預處理后才可進行生化處理［14］。進水BOD5/CODCr的累積分布如圖 7所示。由圖 7可見，BOD5/CODCr為0．05～0．9，平均值和中間值均為 0．4，全年BOD5/CODCr分布在0．4～0．6的概率為39．4%，分布在0．2～0．4的概率為50．2%，BOD5/CODCr＜0．1的概率為1．0%，說明天津城市污水可生化性較好，雖然存在一定的難生物降解性污染物，但仍適應生物處理。

圖7 BOD5/CODCr的概率分布Fig．7 The probability distribution of BOD5/CODCr

2.3.2 進水BOD5/TN

反硝化反應的順利進行需要有充足的碳源，如果碳源不足會嚴重抑制反硝化反應，降低系統(tǒng)總氮去除率;相反如果進水碳源過高會增加系統(tǒng)耗氧量，嚴重抑制硝化反應，降低硝化效率，相應系統(tǒng)總氮去除率也不會太高，因此需維持進水的BOD5/TN［15－16］。一般來說，當污水的BOD5/TKN為4～6時，可認為碳源充足［16－17］。BOD5/TN的概率分布如圖8所示。由圖8可見，BOD5/TN的平均值為4．5，中間值為3．8，BOD5/TN＜4的概率為54．4%，分布在4～6的概率為28．7%，表明多數(shù)情況進水反硝化碳源不足，應注意合理控制外碳源投加量。

圖8 BOD5/TN的概率分布Fig．8 The probability distribution of BOD5/TN

2.3.3 進水BOD5/TP和TN/TP

污水的BOD5/TP是評價采用生物除磷是否可行的主要指標。進水的C/P必須足夠高(特別是聚磷菌厭氧條件下)。若該比值過低，聚磷菌在厭氧池放磷時釋放的能量不能很好地被用來吸收和貯藏溶解性有機物，影響該類細菌在好氧池的吸磷，從而使出水磷濃度升高，一般認為BOD5/TP＞20有較好的磷去除率。比值越大，才能保證聚磷菌有著足夠的基質需求，除磷效果越好［16］。BOD5/TP的概率分布如圖9所示。由圖9可見，BOD5/TP平均值為36．4，中間值為35．1，BOD5/TP＞20的累積概率為77．9%，表明進水在大多數(shù)情況下均可以滿足生物除磷的需求。

圖9 BOD5/TP的概率分布Fig．9 The probability distribution of BOD5/TP

TN/TP的概率分布如圖10所示，一般來說，當污水中BOD5∶N∶P達到100∶5∶1時［16］，即可滿足微生物生長對氮、磷的需求。由圖10可以看出，TN/ TP平均值為9．7，中間值為9．1，分布在5～15的概率為73．6%，且TN/TP＞5的概率為81．6%。由BOD5/TP和TN/TP值可以判斷，氮磷可滿足微生物生長的需求。

圖10 TN/TP的概率分布Fig．10 The probability distribution of TN/TP

3 結論

(1)天津市污水中BOD5、CODCr、SS、NH3|N、TN和TP全年濃度均呈正偏態(tài)分布，月中間值分布在95～140、280～370、130～155、15～30、20～35和3～5 mg/L。各污染指標隨季節(jié)變化呈現(xiàn)出較一致的規(guī)律性。

(2)進水各常規(guī)水質指標間存在較好的一元線性關系，其中BOD5與TP相關性最為顯著，相關系數(shù)(R2)為0．983，說明有機污染物和TP間存在較強的相關關系。CODCr與TN的相關關系相對弱些，R2為0．893，其他指標間R2均大于0．900。

(3)BOD5/CODCr，BOD5/TN，BOD5/TP，TN/TP分析表明，天津城市污水可生化性較好，有機物可以滿足生物除磷的需求，N和P可基本滿足微生物生長對氮、磷的需求，但多數(shù)情況下反硝化碳源不足，應注意合理控制外碳源投加量。

［1］ 賈艷杰，張華兵，王媛．天津可持續(xù)發(fā)展的水資源約束與對策［J］．城市環(huán)境與城市生態(tài)，2003，16(1):76－78．

［2］ 唐福生．天津污水再生利用及水質控制［J］．建設科技，2009 (11):56－57．

［3］ 于佳瀛．天津某中水處理廠污水水質分析及水質預警系統(tǒng)研究［D］．天津:天津大學，2009．

［4］ 孟祥和，張牧，張海生，等．天津工業(yè)節(jié)水與效益分析［J］．資源節(jié)約與環(huán)保，2007，23(3):27－30．

［5］ 高傳昌，劉興．城市非常規(guī)水資源的應用研究進展［J］．灌溉排水學報，2007，26(增刊):68－70．

［6］ QIN X S，GAO F R，CHEN G H，et al．Wastewater quality monitoring system using sensor fusion and machine learning techniques［J］．Water Research，2012，46(4):1133－1144．

［7］ OUYANG Y，NKEDI－KIZZA P，WU T Q，et al．Assessment of seasonal variations in surface water quality［J］．Water Research，2006，40(20):3800－3810．

［8］ WEINRICH L A，JJEMBA P K，GIRALDO E，et al．Implications of organic carbon in the deterioration ofwater quality in reclaimed water distribution systems［J］．Water Research，2010，44(18): 5367－5375．

［9］ 孫迎雪，馬磊，吳光學，等．昆明市合流制排水區(qū)域污水處理廠進水水質特征分析［J］．給水排水，2013，39(3):135－139．

［10］ KRISTENSEN T，ATLAND A，ROSTEN T，et al．Important influent－water quality parameters at freshwater production sites in two salmon producing countries［J］．Aquacultural Engineering，2009，41(2):53－59．

［11］ 胡洪營，吳乾元，黃晶晶，等．再生水水質安全評價與保障原理［M］．北京:科學出版社，2011．

［12］ 張廣強．上海市內陸地表水常規(guī)污染指標相關性研究及應用［D］．上海:東華大學，2009．

［13］ 鄧莉娟，王中琪，張清東，等．氧化溝工藝中污水監(jiān)測指標的相關性分析［J］．工業(yè)用水與廢水，2007，38(6):57－59，88．

［14］ 胡洪營，趙文玉，吳乾元．工業(yè)廢水污染治理途徑與技術研究發(fā)展需求［J］．環(huán)境科學研究，2010，23(7):861－868．

［15］ GANIGUE R，VOLCKE E IP，PUIG S，et al．Impact of influent characteristics on a partial nitritation SBR treating high nitrogen loaded wastewater［J］．Bioresource Technology，2012，111: 62－69．

［16］ 鄭興燦，李亞新．污水除磷脫氮技術［M］．北京:中國建筑工業(yè)出版社，1998．

［17］ 張自杰，林榮忱，金儒霖．排水工程［M］．北京:中國建筑工業(yè)出版社，1999．?

Statistical Analysis of Influent Quality Characteristics of M unicipal W astewater Treatment Plants in Tianjin

SUN Yan1，ZHANG Feng1，HU Hong－ying1，2，NIU Zhang－bin3
1．State Key Joint Laboratory of Environment Simulation and Pollution Control，School of Environment，Tsinghua University，Beijing 100084，China
2．State Environmental Protection Key Laboratory of Microorganism Application and Risk Control，Graduate School at Shenzhen，Tsinghua University，Shenzhen 518055，China
3．Water Administration Office，Ministry of Housing and Urban－rural Development，Beijing 100835，China

The influents of 24 wastewater treatment plants in Tianjin were investigated to analyze the characteristics of water quality，the probability distribution of themain pollution indicators，and the correlations among the organic matter，suspended solids，nitrogen and phosphorus．The results showed that the trend of BOD5，CODCr，SS， NH3－N，TN and TP in the influentswere in accordance with positively skewed distribution，and the concentrations of BOD5，CODCr，SS，NH3－N，TN and TP(monthlymedian)were in the range of95－140 mg/L，280－370 mg/L，130－155 mg/L，15－30 mg/L，20－35 mg/L and 3－5 mg/L，respectively．Simple linear relationships among the six indicators could be observed，where BOD5and TP have a significant correlation，and the correlation coefficients R2was 0．983．The probability distribution of BOD5/CODCrratios between 0．4 and 0．6 was39．4%，and between 0．2 and 0．4 was at50．2%of the cumulative frequency，indicating an adequate biodegradation ability of the influents．The probability distribution of BOD5/TN ratios less than 4 was 54．4%，77．9%of the BOD5/TP ratioswere above 20，which could meet the demand for biological phosphorus removal but not for denitrification．The probability distribution of TN/TP ratios between 5 and 15 was 73．6%with an average of 9．7，and 81．6%of the valueswere above 5．Therefore，the influents could meet the demands of nitrogen and phosphorus formicrobial growth．

wastewater treatment plant;influent quality;characteristics of influent quality;statistical analysis

X703

A

10．3969/j．issn．1674－991X．2014．03．029

1674－991X(2014)03－0173－08

2013－11－21

國家自然科學基金項目(51348011)

孫艷(1977—)，女，助理研究員，博士，主要從事污水再生利用的水質安全評價，sun_yan@tsinghua．edu．cn

*責任作者:胡洪營(1963—)，男，教授，博士，主要從事環(huán)境微生物學與環(huán)境生物技術、污水再生利用的水質安全評價與保障技術等方面的研究，hyhu@tsinghua．edu．cn