范思雨， 龍?zhí)煊澹?**

(1.重慶大學(xué) 三峽庫(kù)區(qū)生態(tài)環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室，重慶 400045；2.重慶大學(xué) 低碳綠色建筑國(guó)際聯(lián)合研究中心，重慶 400045)

不同出流方式下垂直流人工濕地?cái)?shù)值模擬*

范思雨1， 龍?zhí)煊?，2**

(1.重慶大學(xué) 三峽庫(kù)區(qū)生態(tài)環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室，重慶 400045；2.重慶大學(xué) 低碳綠色建筑國(guó)際聯(lián)合研究中心，重慶 400045)

建立了二級(jí)垂直流人工濕地的二維非恒定微生物動(dòng)力學(xué)模型，針對(duì)淹沒(méi)出流方式、不同的等間隔淹沒(méi)出流方式和定水頭出流方式，采用數(shù)值模擬的方法研究二級(jí)垂直流人工濕地處理重慶市徑流雨水的效果;結(jié)果表明：在相同工況下，間隔較大的淹沒(méi)出流方式能達(dá)到最大的有效體積比和水力效率，水流在三層填料中與各種類微生物接觸更加充分、持久，出流方式促進(jìn)了溶解氧的傳質(zhì)，復(fù)氧后溶解氧濃度回升更快、濃度峰值更大、高濃度維持時(shí)間更長(zhǎng)，因此促進(jìn)了由微生物主導(dǎo)的生物化學(xué)降解反應(yīng)，二級(jí)垂直流濕地對(duì)COD、氮和磷的去除效果也優(yōu)于其他出流方式。

人工濕地；出流方式；CW2D；數(shù)值模擬

0 前 言

人工濕地污水處理技術(shù)因其顯著的低能耗、低維護(hù)成本、處理范圍廣和良好的景觀效應(yīng)而被廣泛應(yīng)用[1]，又因其較強(qiáng)的抗沖擊負(fù)荷能力而在蓄納、處理城市徑流雨水上有很好的應(yīng)用前景[2-6]。目前大多數(shù)人工濕地的設(shè)計(jì)與運(yùn)行依然建立在經(jīng)驗(yàn)的基礎(chǔ)上，對(duì)其全運(yùn)行過(guò)程中復(fù)雜的物理-化學(xué)-生物協(xié)同作用認(rèn)知的“黑箱”現(xiàn)象依然存在，而這一作用與人工濕地內(nèi)的水流規(guī)律密切相關(guān)[5]。在模擬研究中將濕地的出流方式抽象成邊界條件，邊界條件的設(shè)定在流體力學(xué)的計(jì)算中極大地影響了水流規(guī)律，因此對(duì)人工濕地出流方式的探討和研究對(duì)提高其處理效率、指導(dǎo)工程設(shè)計(jì)具有重要意義。

目前對(duì)于人工濕地進(jìn)/出水方式的研究主要是兩個(gè)方面：一是討論穿孔管在進(jìn)水處和出水處的配/集水的均勻性；二是討論進(jìn)/出水口在垂直方向上的位置(如上、中、下)及其組合搭配對(duì)配水均勻性和處理效果的影響。張濤等[7]將潛流濕地實(shí)驗(yàn)裝置的進(jìn)水側(cè)從上到下安裝3個(gè)進(jìn)水口，出水側(cè)從上到下安裝4個(gè)出水口，將這些進(jìn)/出水口兩兩組合以示蹤劑實(shí)驗(yàn)為手段考察潛流濕地的水力效率，實(shí)驗(yàn)結(jié)果表明進(jìn)水口和出水口都在上層是潛流濕地的水力效率和有效體積比均大于其他組合。孔德川等[8]采用計(jì)算流體力學(xué)軟件Fluent中的多孔介質(zhì)模型對(duì)分層式潛流人工濕地內(nèi)部的流場(chǎng)進(jìn)行模擬，并分析了不同進(jìn)出口位置對(duì)濕地水力學(xué)行為的影響，結(jié)果表明當(dāng)進(jìn)水口置于濕地基質(zhì)表層、出水口置于基質(zhì)層的中上部時(shí)系統(tǒng)的水力效率較高。牛瑞華等[9]通過(guò)對(duì)水平潛流人工濕地9種不同進(jìn)出水方式進(jìn)行NaCl脈沖示蹤實(shí)驗(yàn)，研究進(jìn)出水位置對(duì)濕地水力效率影響，結(jié)果表明：將出口位置設(shè)置在上部，將轉(zhuǎn)角處設(shè)計(jì)為圓角有助于減小“死區(qū)”面積，提高污染物去除效率。Gunter Langergraber[10]對(duì)人工濕地的數(shù)值模擬展開了廣泛深入的研究，2002年提出人工濕地?cái)?shù)值模擬模型(CW2D，Constructed Wetland 2D)，2005年與HYDRUS開發(fā)者Jirka Simunek合作將CW2D模型以拓展程序包的形式應(yīng)用于HYDRUS-2D中。

在查閱資料的范圍內(nèi)，還未存在利用CW2D人工濕地模型系統(tǒng)研究濕地水平面出流方式對(duì)垂直流潛流人工濕地處理效果和水流規(guī)律影響的先例，而在流體數(shù)值模擬研究中出流方式(水流邊界條件)深刻地影響了水流規(guī)律。利用CW2D人工濕地模型模擬二級(jí)垂直流潛流人工濕地在不同出流方式下的運(yùn)行特征及處理效果，研究探討不同出流方式對(duì)其處理效果的影響，以進(jìn)一步揭示二級(jí)垂直流潛流人工濕地運(yùn)行規(guī)律并為其工程設(shè)計(jì)提供理論指導(dǎo)。

1 材料與方法

1.1 數(shù)學(xué)模型

1.1.1 水分運(yùn)移模型

土壤水分運(yùn)動(dòng)模型，用修正的Richards方程描述：

其中，θ為土壤體積含水量；t為時(shí)間，單位視具體問(wèn)題而定；x、y、z為空間坐標(biāo)，z向上為正，單位視具體問(wèn)題而定；K(θ)為非飽和土壤導(dǎo)水率；D(θ)為非飽和土壤水?dāng)U散率[11]。

采用Van Genuchten-Mualem方程來(lái)模擬土壤水分特性和滲透系數(shù)(導(dǎo)水率)：

其中，θ(h)為土壤體積含水量；K(h)為土壤非飽和滲透系數(shù)(導(dǎo)水率)；θr為土壤剩余體積含水率；θs為土壤飽和體積含水率；α和n是經(jīng)驗(yàn)擬合參數(shù)(或曲線性狀參數(shù))(m= 1-1/n)；Ks為飽和導(dǎo)水率；Se為有效體積含水量；l為經(jīng)驗(yàn)擬合參數(shù)，通常取平均值0.5[12]。

1.1.2 溶質(zhì)運(yùn)移模型

溶質(zhì)在土壤中受對(duì)流和彌散兩種運(yùn)動(dòng)的影響，在模型中用對(duì)流-彌散方程表示溶質(zhì)的運(yùn)移：

其中，ci為溶質(zhì)濃度；q為水流通量；Di為擴(kuò)散度[12]；ri代表微生物的反應(yīng)。

1.2 參數(shù)確定

1.2.1 水分運(yùn)移參數(shù)

水分運(yùn)移模型采用單孔隙模型中的Van Genuchten-Mualem模型，不考慮水分滯后效應(yīng)。Van Genuchten-Mualem公式中的各項(xiàng)參數(shù)首先由LANGERGRABER G[10]開發(fā)的并經(jīng)過(guò)各個(gè)尺度實(shí)驗(yàn)驗(yàn)證的濕地模型中獲取，然后經(jīng)過(guò)多次反算優(yōu)化，濕地濾層從上到下3種材料的水力特征參數(shù)值見表1。

表1 水力特征參數(shù)

1.2.2 CW2D參數(shù)

在CW2D中包含微生物對(duì)有機(jī)污染物、氮、磷的好氧、缺氧運(yùn)移轉(zhuǎn)化和降解過(guò)程。HYDRUS-CW2D中針對(duì)濕地的多種物質(zhì)反應(yīng)運(yùn)移模型包括了12種物質(zhì)和9個(gè)反應(yīng)過(guò)程。

其中12種物質(zhì)分別是：溶解氧、易生物降解COD(CR)、難生物降解COD(CS)、生物惰性COD(CI)、非自養(yǎng)微生物、亞硝化單胞菌種、硝化菌種、氨氮、亞硝態(tài)氮、硝態(tài)氮、氮?dú)?、無(wú)機(jī)磷。

9個(gè)生物化學(xué)過(guò)程分別是：水解、非自養(yǎng)菌的好氧生長(zhǎng)、基于亞硝態(tài)氮的反硝化、基于硝態(tài)氮反硝化、非自養(yǎng)菌的細(xì)胞溶解、基于氨氮的亞硝化單胞菌的好氧生長(zhǎng)(硝化第一步)、亞硝化單胞菌的細(xì)胞溶解、基于亞硝態(tài)氮的硝化細(xì)菌的好氧生長(zhǎng)(硝化第二步)、硝化細(xì)菌的細(xì)胞溶解。

CW2D模型參數(shù)首先采用軟件系統(tǒng)自帶數(shù)據(jù)，根據(jù)模擬結(jié)果再加以微調(diào)，如果沒(méi)有豐富的生物化學(xué)動(dòng)力學(xué)基礎(chǔ)或大量實(shí)驗(yàn)結(jié)果支撐，不建議大范圍或大幅度修改這些參數(shù)。

1.3 模型驗(yàn)證

表2 模型驗(yàn)證所使用的進(jìn)水各物質(zhì)濃度

2 模擬研究

2.1 幾何模型的建立

二維的二級(jí)垂直流潛流人工濕地模型尺寸見圖2。濾料分為3層，填料尺寸從上至下分為中、細(xì)、粗，上層(4～8 mm)和下層(16～32 mm)均為碎石，中層為砂(0.06～4 mm)。進(jìn)水方式為左側(cè)進(jìn)，右側(cè)出，人為控制進(jìn)水流量，模擬時(shí)間4 h，在最初和運(yùn)行2 h后兩次進(jìn)水。

圖2 濕地幾何模型示意圖Fig.2 Schematic diagram of the 2-stage vertical flow constructed wetlands

利用有限單元法完成以上數(shù)學(xué)模型的計(jì)算，時(shí)間采用隱式差分法離散[13]。在計(jì)算區(qū)域內(nèi)采用三角形網(wǎng)格劃分，為了實(shí)現(xiàn)準(zhǔn)確的網(wǎng)格劃分，在繪制幾何模型時(shí)，除了給出濕地模型外輪廓線之外，還在深度方向上按照一定梯度附加一些橫線以適當(dāng)加密網(wǎng)格并使網(wǎng)格劃分得更有秩序，在一、二級(jí)濕地交界處也適當(dāng)加密網(wǎng)格，以提高計(jì)算精確度和準(zhǔn)確度。

在模擬區(qū)域不同深度處設(shè)置6處觀察點(diǎn)，根據(jù)之前多次模擬結(jié)果的規(guī)律，觀察點(diǎn)設(shè)定基本按照此濕地中的水流方向及污染物的遷移方向，因而能大致反映此點(diǎn)處同一水平面上濃度的均值，這點(diǎn)處濃度的變化情況也能反映全過(guò)程污染物濃度隨時(shí)間的變化情況。觀察點(diǎn)情況：觀察點(diǎn)1,一級(jí)濕地水面0 cm處(靠近進(jìn)水處)；觀察點(diǎn)2,一級(jí)濕地距水面以下5 cm處；觀察點(diǎn)3,一級(jí)濕地距水面以下15 cm處；觀察點(diǎn)4,一級(jí)濕地距水面以下40 cm處；觀察點(diǎn)5,二級(jí)濕地距水面以下10 cm處；觀察點(diǎn)6,二級(jí)濕地水面0 cm處(靠近濕地最遠(yuǎn)端出口處)。

2.2 初始條件

需要設(shè)定的初始條件為上述12種物質(zhì)的濃度分布，模擬采用人工濕地處理重慶市徑流雨水。進(jìn)水中各物質(zhì)濃度查閱資料獲得[14-16]，進(jìn)水濃度采用雨水的EMC值，初始條件設(shè)定按劃分的網(wǎng)格點(diǎn)賦值，形成二維濃度分布圖。模擬中進(jìn)水的各物質(zhì)濃度見表3。

表3 模擬中的進(jìn)水濃度

2.3 邊界條件

2.3.1 水流邊界條件

將不同的出流方式在模擬中抽象成不同的邊界條件，主要討論水流邊界改變對(duì)濕地水流規(guī)律及處理效果的影響，按照人工濕地實(shí)際設(shè)計(jì)經(jīng)驗(yàn)并結(jié)合工程實(shí)踐，在HYDRUS中采用4種邊界條件組合：進(jìn)、出水兩側(cè)為大氣-滲出面邊界,沒(méi)有進(jìn)行任何改進(jìn)，沿寬度均勻布水的進(jìn)水方式和自由淹沒(méi)出流的出水方式，進(jìn)/出水兩側(cè)均暴露在大氣中；大氣-不完全滲出面邊界一,出水側(cè)間隔出水；大氣-不完全滲出面邊界二,出水側(cè)間隔增大為前者的兩倍；大氣-定水頭邊界,出水側(cè)保持水面高度不變，即一直維持最初的水位。其中大氣-滲出面邊界簡(jiǎn)稱為A-S邊界，大氣-定水頭邊界簡(jiǎn)稱為A-C邊界。

控制其他變量，以上各種邊界條件均對(duì)應(yīng)同一進(jìn)水工況：600 L/d，每2 h進(jìn)水50 L；在軟件中對(duì)應(yīng)的時(shí)變邊界條件(Time Variable Boundary Conditions)處設(shè)置，模擬其中4 h內(nèi)的2次進(jìn)水過(guò)程。

2.3.2 溶質(zhì)運(yùn)移邊界條件

溶質(zhì)運(yùn)移出水進(jìn)水側(cè)的邊界條件均為第三類邊界條件(third-type)，沿邊界的溶質(zhì)通量是確定時(shí)使用第三類條件。Genuchten V等[16]建議盡量調(diào)用第三類邊界條件，因?yàn)榈谌愡吔鐥l件更接近現(xiàn)實(shí)情況，并且在模擬系統(tǒng)中保留了溶質(zhì)質(zhì)量。

3 結(jié)果與討論

3.1 不同邊界條件對(duì)處理效果的影響

3.1.1 A-S邊界

在此邊界條件下的溶解氧、易生物降解COD(CR)、氨氮、硝態(tài)氮和無(wú)機(jī)磷隨時(shí)間變化曲線如圖3。

從圖3中所示可以明顯看出兩次進(jìn)水這一過(guò)程，首次進(jìn)水是在0.05 h處，可以看到觀察點(diǎn)1(第一級(jí)濕地水面0 cm處)和觀察點(diǎn)2(非飽和帶5 cm處)溶解氧濃度在進(jìn)水后迅速升高；而觀察點(diǎn)3(水面下15 cm)則由于氧氣的溶解和傳遞過(guò)程，溶解氧濃度緩慢升高。

易生物降解COD也因?yàn)閮纱芜M(jìn)水而有兩次升高和下降的過(guò)程，可以看到水下15 cm處因?yàn)槲⑸餄舛容^高并且溶解氧相對(duì)充足CR濃度在進(jìn)水后迅速下降。

IP濃度與溶解氧濃度變化相符，當(dāng)還未進(jìn)水時(shí)溶解氧濃度下降，聚磷菌釋磷IP升高，氧濃度升高后磷被聚磷菌吸收，同樣的,隨著深度增加溶解氧成為限制因素，聚磷菌吸磷作用放緩(15 cm處)。

圖3 A-S邊界條件下各組分濃度隨時(shí)間變化圖Fig.3 Simulated time series of 5 components on A-S boundary condition

通過(guò)圖4可知，在A-S邊界條件下，在兩次進(jìn)水后，CR濃度顯著提高，通過(guò)濕地作用，CR濃度下降，出水面CR濃度降至50 mg/L及以下。

圖4 A-S邊界條件時(shí)CR濃度變化圖Fig.4 The concentration changing of CR on A-S boundary condition

3.1.2 A-S1和A-S2邊界

這兩個(gè)邊界是將A-S邊界的出水邊界設(shè)置成等間隔出水形式，A-S2邊界的間隔是A-S1的兩倍，這樣設(shè)計(jì)的目的是減少出水側(cè)的“死水區(qū)”，因?yàn)椤八浪畢^(qū)”主要集中在出水口一側(cè)，并且是在遠(yuǎn)離出水口的位置[8]，這樣間隔的出水設(shè)置，是為了最大程度地調(diào)動(dòng)整個(gè)出水面。圖5是A-S、A-S1、A-S23種不同邊界條件下觀察點(diǎn)3(第一級(jí)濕地水下15 cm處)溶解氧、易生物降解COD(CR)、氨氮、硝態(tài)氮和無(wú)機(jī)磷隨時(shí)間變化曲線的比較。

圖5 在3種邊界條件下觀察點(diǎn)3處各組分濃度變化對(duì)比圖Fig.5 The comparison of simulated time series of 5 components on A-S/A-S1/A-S2 boundary conditions

有效體積比(e)表示濕地中水流能夠有效通過(guò)的體積占總體積的比例，一般e越大，表示人工濕地中用于去除污染物的有效體積比例越大。其定義為

其中Veffective為有效體積，Vtotal為總體積[19]；兩種間隔設(shè)置的出水方式減少了出水一側(cè)的死水區(qū)，提升了濕地的有效面積，因此濕地的去除效率隨之提升。

如圖6所示，3種邊界條件下，顯然A-S2最大限度地利用了出水面，而A-S和A-S1邊界在此時(shí)都有“死水區(qū)”的存在，這也說(shuō)明在A-S2邊界下，濕地的有效體積比最大，其用于處理污染物的體積也最大。

通過(guò)比較A-S1和A-S2可以看出，A-S2邊界在濕地處理過(guò)程中雖不明顯但還是相對(duì)優(yōu)于A-S1邊界的，這是因?yàn)槌鏊g隔較小的A-S1邊界，仍存在“死水區(qū)”，A-S2邊界的出水口和間隔都增大了一倍，出流阻礙減少，水頭損失亦減少，進(jìn)而減少了不良干擾、提高了氧氣的傳質(zhì)和污染物的去除效率。此外，在模擬計(jì)算時(shí)，A-S2邊界收斂速度最快，也說(shuō)明在這個(gè)邊界條件下水流模型更加準(zhǔn)確、穩(wěn)定。

圖6 3種邊界條件在同一時(shí)刻出水面情況對(duì)比圖Fig.6 The comparison of effluent surface of 3 boundary conditions on the same time

3.1.3 A-C邊界

滲流邊界是假設(shè)滲流面不飽和處出水流量一直保持為零，因此除了考慮出流側(cè)設(shè)定為滲流邊界外，模擬研究還采用了定水頭邊界，定水頭邊界是指設(shè)定為定水頭邊界點(diǎn)的壓力水頭始終不變等于初始值，即水位不變，這與滲流面邊界的假設(shè)相同。在這種條件下將A-C與A-S2邊界在同等工況下各組分濃度變化作比較，得到了一組完全重合的曲線，說(shuō)明保持出流面飽和狀態(tài)，保持水位不變?cè)贑W2D模型中模擬二級(jí)垂直流潛流濕地能得到同樣的結(jié)果，但是考慮到滲流面邊界最大程度地還原了實(shí)際濕地運(yùn)行情況，因此推薦使用滲流面邊界。

4 結(jié) 論

利用CW2D/HYDRUS-2D模型對(duì)不同邊界條件下二級(jí)垂直流潛流人工濕地的水分、溶質(zhì)時(shí)空運(yùn)移和人工濕地中多組分生化反應(yīng)進(jìn)行模擬，以實(shí)驗(yàn)室實(shí)測(cè)數(shù)據(jù)驗(yàn)證二級(jí)垂直流模型，結(jié)果表明：模型有較好的仿真度，能夠基本正確模擬濕地中的水鹽變化和污染物降解過(guò)程，也符合污染物去除的一般規(guī)律，可以用于模擬二級(jí)垂直流潛流人工濕地流處理過(guò)程。

在不同邊界條件下，對(duì)濕地運(yùn)行過(guò)程進(jìn)行模擬研究，結(jié)果表明：在大氣-滲流面邊界二(A-S2)條件下，濕地中生化反應(yīng)效果最佳，此邊界條件提供了最好的水流狀態(tài)，有效體積比最大，進(jìn)而有相對(duì)良好的溶解氧傳遞和分布，促進(jìn)好氧硝化細(xì)菌的作用，能將氨氮快速大量地轉(zhuǎn)化成硝態(tài)氮；此外相對(duì)優(yōu)良的水流環(huán)境和氧分布也促進(jìn)了COD和磷的去除效率。

因此，在人工濕地出水面設(shè)置合理的等距間隔可以顯著提升進(jìn)濕地的處理效果。等距間隔可以在實(shí)際工程中搭配景觀設(shè)計(jì)實(shí)現(xiàn)。重慶市初期雨水中COD含量較高，夏季雨水形成徑流流量較大，依靠大容量景觀化的人工濕地接納、處理雨水徑流能達(dá)到良好的處理效果，既節(jié)約了處理成本又美化了城市環(huán)境。

[1] 王曉雪, 李杰, 鐘成華, 等.人工濕地用于三峽庫(kù)區(qū)小城鎮(zhèn)生活污水處理的研究[J]. 重慶工商大學(xué)學(xué)報(bào)(自然科學(xué)版), 2011, 28(2):195-198

WANG X X, LI J, ZHONG C H. Study on Sewage Treatment of Small Towns in Three Gorges Reservoir Area by Constructed Wetlands[J]. Journal of Chongqing Technology and Business University (Natural Science Edition), 2011,28(2):195-198

[2] PALFY T G, LANGERGRABER G. The Verifycation of the Constructed Wetland Model No. 1 Implementation in HYDR-US using Column Experiment Data[J]. Ecological Engineering, 2014,68(7): 105-115

[3] RIZZO A, LANGERGRABER G, GALVAO A, et al. Modelling The Response of Laboratory Horizontal Flow Con-structed Wetlands to Unsteady Organic Loads with HYDRUS-CWM1[J]. Ecological Engineering,2014, 68(7): 209-213

[4] HABERL R, GREGO S, LANGERGRA-BER G, et al. Constructed Wetlands for the Treatment of Organic Pollutants[J]. Journal of Soils and Sediments, 2003,3(2):109

[5] KADLEC R H,WALLACE. Treatment Wetlands[M]. Boca Raton：CRC Press,2004

[6] CAMPBELL C S, OGDEN M.Constructed Wetlands in the Sustainable Landscape[M]. New York：John Wiley& Sons，1999

[7] 張濤, 宋新山. 不同運(yùn)行條件下潛流人工濕地的水力效率分析[J]. 環(huán)境污染與防治, 2010,32(10):37-40

ZHANG T, SONG X S. Hydraulic Efficiency of Subsurface Flow Constructed Wetlands under Different Operating Conditions[J]. Environmental Pollution & Control,2010, 32(10):37-45

[8] 孔德川, 丁愛(ài)中, 鄭蕾, 等. 分層式潛流人工濕地水力學(xué)特性數(shù)值模擬與分析[J]. 環(huán)境工程學(xué)報(bào), 2011,5(4):741-744

KONG D C, DING A Z, ZHENG L, et al. Numerical Simulation and Analysis of Hydraulic Characteristics of Layered Subsurface Flow Constructed Wetland[J]. Chinese Journal of Environmental Engineering,2011, 5(4):741-744

[9] 牛瑞華, 宋新山, 王宇暉. 進(jìn)出水位對(duì)水平潛流人工濕地水力效率影響可視化分析研究[J]. 環(huán)境污染與防治, 2014, 36(4):11-16

NIU R H, SONG X S, WANG Y H. Hydraulic Efficiency of Horizontal Subsurface Flow Constructed Wetlands with Different Inlet-outlet Location Configuration and Its Visual Analysis[J]. Environmental Pollution & Control, 2014, 36(4):11-16

[10] LANGERGRABER G, IMUNEK J. Modeling Variably Saturated Water Flow and Multicomponent Reactive Transport in Constructed Wetlands[J]. Vadose Zone Journal, 2005,4(4):924-938

[11] 馬海燕, 王昕, 張展羽, 等. 基于HYDRUS-3D的微咸水膜孔溝灌水鹽分布數(shù)值模擬[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2015(2):137-145

MA H Y, WANG X, ZHANG Z Y, et al. Numerical Simulation of Water-salt Distribution under Brackish Water Film Hole Furrow Irrigation Based on HYDRUS-3D Model[J]. Transactions of the Chinese Society of Agricultural Machinery, 2015, 46(2):137-145

[12] FOURNEL J, FORQUET N, MOLLE P, et al. Modeling Constructed Wetlands with Variably Saturated Vertical Subsurface-flow for Urban Stormwater Treatment[J]. Ecological Engineering, 2013,55(6):1-8

[13] 李耀剛, 王文娥, 胡笑濤. 基于HYDRUS-3D的涌泉根灌土壤入滲數(shù)值模擬[J]. 排灌機(jī)械工程學(xué)報(bào), 2013,31(6):546-552

LI Y G, WANG W E, HU X T. Numerical Simulation of Soil Water Infiltration Under Bubbled Root Irrigation Based on HYDRUS-3D[J]. Journal of Drainage and Irrigation Machinery Engineering, 2013, 31(6)：546-552

[14] 肖海文. 城市徑流特征與人工濕地處理技術(shù)研究[D]. 重慶：重慶大學(xué), 2010

XIAO H W. Study on Characteristic of Urban Runoff and the Technology of Constructed Wetland for Urban Runoff Treatment[D].Chongqing： Chongqing University, 2010

[15] 謝瑤.重慶市海綿城市建設(shè)技術(shù)模式研究[D].重慶：重慶大學(xué)，2016

XIE Y. Technical Mode of Chongqing Sponge City Construction[D]. Chongqing：Chongqing University, 2016

[16] 李彥庭.生態(tài)濾池示范工程對(duì)趙家溪水質(zhì)凈化效果研究[D].重慶：重慶大學(xué)，2016

LI Y T. Research on Purification Effects of Ecological Filter Demonstration Project in Zhaojia Creek[D]. Chongqing: Chongqing University, 2016

[17] GENUCHTEN M T V, PARKER J C. Boundary Conditions for Displacement Experiments through Short Laboratory Soil Columns[J]. Soil Science Society of America Journal, 1985,48(4):703-708

[18] THACKSTON E L, SHIELDS F D, SCHROEDER P R. Residence Time Distributions of Shallow Basins[J]. Journal of Environmental Engineering, 1987,113(6):1319-1332

Modeling Treatment Efficiency of Two-stage Subsurface Vertical Flow Constructed Wetlands with Different Ways of Effluent

FANSi-yu1,LONGTian-yu1,2

(1. Key Laboratory of Eco-Environment of Three Gorges Region of Ministry of Education, Chongqing University, Chongqing 400045, China; 2.National Centre for International Research of Low-carbon and Green Buildings, Chongqing University, Chongqing 400045, China)

A two-dimension convection-dispersion-biodynamics model was given to simulate the treatment efficiency of two-stage subsurface vertical flow constructed wetlands for urban stormwater runoff of Chongqing. Focus on the following 4 types of boundary conditions: Atmospheric Boundary-Seepage Face, Atmospheric Boundary-incomplete Seepage Face 1 (uniformly-spaced effluent), Atmospheric Boundary- incomplete Seepage Face 2 (A-S2’s space is 2 times larger than A-S1’s), Atmospheric Boundary- Constant Head. The simulated results show that under the same working condition, the A-S2 boundary condition has an obviously better performance on the transportation of dissolved oxygen and the treatment efficiency of COD,nitrogen and phosphorus removal , since the specific effluent mode offers a better hydraulic characteristic of the constructed wetland.

constructed wetlands; effluent mode; CW2D; numerical simulation

X703

A

2017-03-10；

2017-04-18.

“十二五”國(guó)家科技支撐計(jì)劃重點(diǎn)項(xiàng)目(2011BAD31B03).

范思雨(1992-) ，女，四川成都市人，碩士研究生，從事人工濕地?cái)?shù)值模擬研究.

**通訊作者：龍?zhí)煊?1960-)，女，重慶人，教授，博士生導(dǎo)師，從事水污染控制與水環(huán)境模擬研究.

責(zé)任編輯：田靜