廖新弟，駱世明，吳銀寶，汪植三
（華南農業(yè)大學 熱帶亞熱帶生態(tài)研究所，廣東廣州 510642）

　　摘要：以風車草人工濕地處理豬場廢水，研究人工濕地處理廢水有機物的動態(tài)模型。人工濕地一年四季連續(xù)運行，秋李和春季濕地以相同進水濃度和不同水力停留時間的方式運行；冬季和夏季以相同水力停留時間和不同進水濃度的方式運行。COD作為有機物降解的建模指標。結果表明，濕地進出水COD服從指數(shù)方程規(guī)律。提出了基于運行溫度和進水濃度的濕地出水COD預測模型，預測結果誤差在10％以內。
　　關鍵詞：人工濕地；風車草；廢水處理；有機物；模型
　　中圖分類號：X703.1　　 文獻標識碼：A　　 文章編號：1009—2455(2004)04—0023—04

A Study of Dynamic Model of Treatment of Organic Matters in Sewage Using Constructed Wetland

LIA0 Xin-di，LUO Shi-ming，WU Yin—bao，WANGZhi-san
(Institute of Tropical and Subtropical Ecology of South China A gricultural University , Guangzhou 510642，China)

　　Abstract：A constructed wetland system (CWS) where cyperas alternifolius had been planted was used for the treatment of sewage from hoggery in order to develop a dynamic model of treatment of organic matters with constructed wetland.The constructed wetland operated continuously at a11 seasons.The wetland operated with the same concentration of influent and different hydraulic retention time (HRT) in autumn and spring，while in winter and summer it operated with the same HRT (3 days) and different concentration of influent.Chemical oxygen demand (COD) was used as an indicator for modeling the degradation of organic matters.The results showed that the COD of the influent and effuent of the wetland followed the rule of exponential equation.A model was proposed for anticipating the COD of the effuent water from the wetland based on the operating temperature and the concentration of influent water，the error of which was within 10％.
　　Key words：constructed wetland；cyperus altrnlifolius；wastewater treatment；organic matter；model

　　人工濕地已廣泛應用于廢水處理，然而，人工濕地的設計和運行常常缺乏模型的指導。本研究以風車草潛流式人工濕地處理養(yǎng)豬場廢水為研究對象，重點討論在同一廢水來源、相同的潛流式人工濕地和統(tǒng)一的水力停留時間下(3 d)，進水有機負荷(以COD表示)和氣溫對COD降解常數(shù)左的影響，以及基于進水有機負荷和氣溫丸值與出水水質預測模型。

1模擬人工濕地的構成

　　 人工濕地的剖面如圖1所示，長度為1 m，寬0.5m，高0.8m，由磚和水泥砌成。濕地內填充粒徑3—5cm的碎石60cm厚作為處理床，處理床上種植風車草，試驗在大棚內進行。風車草植株根系發(fā)達，分蘗數(shù)14-16個／株，高度135-145sm，每組8株，分兩行栽種。濕地進水通過位于濕地前部的進水槽從處理床前端底部分多孔均勻進水，從另一端上部多孔均勻出水。

2COD降解的數(shù)學模型

　　根據反應器中反應物的流動特征和凈化機理，通常把人工濕地看作推流反應器^[1-3]。在理想的推流反應器中，反應物濃度隨空間和時間而變化，并遵循：

　　 Y_t=Y_oexp(-kt)

式中：t——反應器物料的理論停留時間，d；
　　　Y_t——出口COD的質量濃度，mg／L；
　　　Y_o——進口COD的質量濃度，mg／L；
　　　K——COD降解常數(shù)^[4]。

　　在上述方程中，兌反映濕地廢水COD的降解速度，它與廢水性質、進水特征、廢水處理系統(tǒng)的整體特征以及系統(tǒng)運行的環(huán)境條件等有關^[5]。

3實驗方法

　　測試前先運行2個月，以讓處理床掛膜；運行時每隔3 d進1次廢水，進水量90L／室，進水濃度低于當季正式測試時的最低進水濃度。
　　為了研究季節(jié)性氣溫(水溫)對濕地降解COD的影響，將試驗按4個不同季節(jié)進行測試，試驗連續(xù)進行1 a。其中，秋季和春季運行采用同一個進水濃度范圍，廢水在濕地的停留時間(HRT)采用5個方式(5，4，3，2，1 d)，秋季春季試驗平均進水COD的質量濃度為736.29土29.63 mg／L，5種停留時間運行下的水力負荷分別為72，90，120，180和360L／(m²·d)，有機負荷(以COD計)依次為：53.01，66.26，88.35，132.53和265.06 g／(m²·d)。冬季和夏季試驗運行采用同一個停留時間(3d)，5個不同的進水濃度范圍(見表1)。

　　試驗期間各組每天連續(xù)進水，進水是豬場廢水經沉淀和厭氧處理后的出水再按試驗所需用清水調配而成。在同一季節(jié)內，試驗不同停留時間或不同進水水質的運行和測試均重復3次。
　　每一個季節(jié)測試結束，將試驗植物沿床面上30cm剪割掉，各組繼續(xù)在低濃度下繼續(xù)運行。

4 結果與討論

4.1 進水濃度對COD降解系數(shù)A的影響
　　冬、夏季濕地去除COD服從指數(shù)方程規(guī)律。將風車草濕地夏季、冬季停留時間等于3 d時Y_o和Y_t，的5組測定結果分別代人指數(shù)方程Y_t=Y_o·e^(-kθ)，獲得夏季、冬季風車草濕地降解COD的5個k值，見表1。
　　根據表1中Y_o與k的對應關系，可以分別求出冬、夏季COD的降解系數(shù)k(d^-1)與進水COD的質量濃度(Y_o，mg／L)的回歸方程，分別為：
　　夏季：k=0.0003Y_o+0.3727，r=0.9019　 (1)
　　冬季：k=-0.000 02Yo+0.4791，r=0.2083　　 (2)
4.2 溫度對COD降解系數(shù)k的影響
　　將4個進水COD的質量濃度：490.35,867.07，983.38和1440.20mg／L，依次代入冬季k與Y_o的關系式k=-0.000 02 Y_o+0.479 1，可估算出4個k，依次是：0.469 3，0.461 8，0.459 4和0.450 3 d^-1。同理，將上述4個進水COD值代入夏季k與Y_o的關系式k=0.000 3Y_o+0.372 7，求出4個相應的k，依次是：0.519 8，0.632 8，0.721 5和0.725 0d^-1。
　　假設在同一進水濃度下，k與溫度成線性關系，根據上述4種進水濃度下，冬、夏季各4個k值，以及冬、夏季氣溫(分別為17.9℃和27.2℃)，可以分別求出上述4種進水濃度下，k隨溫度θ(℃)變化的直線方程，見表2。

　　根據表2的4種不同進水濃度下、k隨溫度θ變化的估算方程，可以計算出秋季(溫度21.4℃)、上述4種進水濃度下的COD降解常數(shù)&值分別為：0.488 2，0.526 2，0.558 1和0.553 5 d^-1。同理，可以計算出春季(溫度23.8℃)、上述4種進水濃度下的COD降解常數(shù)禿值分別為：0.501 1，0.570 3，0.625 8和0.624 3 d^-1。
　　這樣，4種進水濃度、4個季節(jié)(溫度)下人工濕地的COD降解常數(shù)無值匯總在表3。

4.3　 COD降解常數(shù)k的估算模型

　　根據表3中進水濃度與k的對應數(shù)據，可以作出不同溫度下，k與進水濃度之間的關系圖，如圖2所示。圖2表明：冬季低溫環(huán)境下，進水濃度提高，k則減小(θ1)；其它高溫季節(jié)環(huán)境下，k隨進水濃度提高而增大(θ2，θ3和θ4)。

　　圖2中，不同溫度下，k隨進水濃度變化的回歸方程分別為：
　　方程θ1：k=0.4791-(2×10^-5)y_o，r=1.0000　　 (3)
　　方程θ2：k=0.4656+(7×10^-5)y_o，r=0.8489　　 (4)
　　方程θ3，k=0.4563+0.0001y_o，　 r=0.8735　　(5)
　　方程θ4：k=0.4432+0.0002y_o，　 r=0.8856　　(6)
　　上述θ₁—θ₄的4個方程可用k=a+by_o來表示，方程中a值與溫度關系見圖3，b值與溫度關系見圖4。

　　圖3、圖4可見，a值隨溫度上升而下降，而b值隨溫度上升而增加，其回歸關系為：
　　a=0.5482-0.0039 θ，r=1.000 0　　(7)
　　b=0.0004+(2×10^-5)θ，r=0.992 4　 (8)
　　將式(7)和式(8)代人k=a+by_o中，則k與進水濃度的關系式變?yōu)椋?BR>　　k=(-0.0039θ+0.5482)+[(2×10^-5)θ-0.0004]y_o。　　 (9)
　　若式(7)用a= a₁θ+a₂來表示，則a_l=-0.0039，a₂=0.5482
　　同理，若式(8)用b= b₁θ+b₂來表示，則b₁=2 ×10^-5，b₂=-0.0004
　　將k=a+by_o代入y_t=y_oexp(-kt)，
　　則：y_l=y_oexp[-(a+by_o)t]，再將a= a₁θ+a₂和b= b₁θ+b₂代入，
　　則：y₁=y_oexp{-[( a₁θ+a₂)+ (b₁θ+b₂)y_o]t}　　 (10)
　　式(10)中：y_l為出水COD的質量濃度mg／L；y_o為進水COD的質量濃度mg／L；θ為溫度，℃；t為水力停留時間，d;a₁,a₂，b₁和b₂均為系數(shù)，a₁=-0.0039，a₂=0.5482，b_l=2×10^-5，b₂=-0.0004。
　　這樣，通過式(10)，也可以預測出不同溫度、進水COD濃度和停留時間下的風車草人工濕地豬場廢水處理的出水COD值。

5 實測與模擬結果的比較

　　 式(10)是根據人工濕地實際系統(tǒng)的進水濃度和水力停留時間以及運行溫度環(huán)境，預測出人工濕地的出水水質，式中k是關鍵，其估算模型的精度見表4。

　　表4中，k來自表3，k₁根據式(9)計算而來，k₁誤差指與左比較的結果。
　　表4表明，按式9估計的k比實測偏小，誤差約在10％以內。這說明在溫度和進水濃度已知的條件下，人工濕地k可用k=(-0.0039θ+0.5842)+[(2×10^-5) θ-0.0004]y_o來估計。針對不同水力停留時間的人工濕地出水水質可用y_t=y_oexp{-[(a₁θ+a₂)+(b_lθ+b₂)y_o]t}來預測。盡管本研究只采用COD建立風車草人工濕地去除有機物的動力學模型，但其方法同樣適用于其它植被的人工濕地，因為同樣遵從推流反應器原理；此方法也適用于描述服從指數(shù)模型的其它污染指標的降解過程，例如BOD，只是模型參數(shù)需要重新調整。

6結論

　　人工濕地降解有機物服從指數(shù)方程規(guī)律，模型可以用Y_t=Y_o·e^(-kt)表達。本研究建立了風車草人工濕地COD降解動力學模型dy／dt=-kt(t為時間)，研究闡明了溫度和進水濃度對COD降解常數(shù)兌的影響。在溫度和進水濃度已知的條件下，風車草人工濕地COD降解的兌可用k=(-0.0039θ+0.5482)+[(2×10^-5) θ-0.0004]y_o妁來估計；針對不同水力停留時間(d)的風車草人工濕地出水COD(y_t)可用y_t=y_oexp{-[(a₁θ+a₂)+(b_lθ+b₂)y_o]t}來預測。該模型可用于具推流反應器原理的人工濕地有機物的去除過程，模型參數(shù)需針對現(xiàn)場運行條件進行調整。

參考文獻：

[1]丁廷華.污水蘆葦濕地處理系統(tǒng)示范工程的研究[J).環(huán)境科學，1992，13(2)；8—13.
[2]鄭雅杰.人工濕地系統(tǒng)處理污水新模式的探討[J].環(huán)境科學進展，1995，3(6)：1—8.
[3]諸惠昌，胡紀萃，新型廢水處理工藝一人工濕地的設計方法[J].環(huán)境科學，1993，14(2)：39—43.
[4]唐受印，汪大暈，廢水處理工程[M].北京：化學工業(yè)出版社，1998，8—12.
[5]胡康萍.人工濕地設計中的水力學問題研究[J].環(huán)境科學研究，1991，4(5)；11—16.

作者簡介：廖新弟(1968—)，男，福建連江縣人，副教授，博士，華南農業(yè)大學動物科學學院，電話(020)85280279，xdliao2002@yahoo.com.cn。