李世陽,王 立,Prasanta K Kalita,馬 放

(1.城市水資源與水環(huán)境國家重點(diǎn)實(shí)驗(yàn)室(哈爾濱工業(yè)大學(xué))，哈爾濱 150090; 2.伊利諾伊大學(xué)香檳分校 農(nóng)業(yè)與生物技術(shù)工程學(xué)院，美國 Illinois 61802)

淺層排水系統(tǒng)降雨過程氮磷垂直遷移數(shù)值仿真

李世陽1,2,王 立1,Prasanta K Kalita2,馬 放1

(1.城市水資源與水環(huán)境國家重點(diǎn)實(shí)驗(yàn)室(哈爾濱工業(yè)大學(xué))，哈爾濱 150090; 2.伊利諾伊大學(xué)香檳分校 農(nóng)業(yè)與生物技術(shù)工程學(xué)院，美國 Illinois 61802)

淺層排水系統(tǒng)是一種快速排出農(nóng)田積水的有效設(shè)施，但其同時(shí)也改變了水分及面源污染物如氮磷在農(nóng)田系統(tǒng)遷移的速度和途徑.為了對淺層排水農(nóng)田系統(tǒng)中的氮磷垂直遷移過程有全面了解,利用有限元分析技術(shù)對已監(jiān)測的淺層排水系統(tǒng)單元土壤柱進(jìn)行數(shù)值仿真.結(jié)果表明，在5 h 4 mm/h降雨結(jié)束后，降雨量沒有達(dá)到使地表到排水管道之間土壤層完全飽和，從而使隨之溶解的氮磷污染物不足以滲透到排水管，導(dǎo)致沒有氮磷的流失.降雨時(shí)間延長至10 h后，降雨量使土壤層達(dá)到飽和，從而使溶解的污染物產(chǎn)生流失，其中硝酸鹽的流失量明顯大于磷酸鹽，硝酸鹽流失量為9～10 mg/L,磷酸鹽流失量為0.05～0.1 mg/L.在24 h降雨事件后，排水管中硝酸鹽質(zhì)量分?jǐn)?shù)達(dá)最大值.在降雨時(shí)間不斷增加的過程中，硝酸鹽流失呈指數(shù)增長，而磷酸鹽只是線性增長，形成該現(xiàn)象的原因與污染物的本身性質(zhì)有關(guān).

淺層排水系統(tǒng)；面源污染；垂直遷移；有限元；數(shù)值仿真

淺層排水系統(tǒng)廣泛存在于平原地區(qū)的集約化農(nóng)田地域，尤其在美國中部地區(qū)，約70%的農(nóng)田具有淺層排水系統(tǒng)[1]，在中國的西北和江南地區(qū)也有不成規(guī)模的農(nóng)田排水系統(tǒng)，該系統(tǒng)能夠提高產(chǎn)量，但同時(shí)加速了污染物進(jìn)入水體的速度[2].淺層排水系統(tǒng)存在的農(nóng)田氮磷污染物的遷移過程一直是專家學(xué)者們關(guān)注的問題，美國學(xué)者開發(fā)了淺層排水水分垂直遷移過程模型DRAINMOD分析排水管網(wǎng)中水量的大小，并計(jì)算污染物流失量[3-6].但對于污染物在土壤層中遷移的微觀過程并沒有相應(yīng)的研究和報(bào)道，國內(nèi)學(xué)者對該問題的研究也處于初級(jí)階段.本文利用有限元分析的方法[7-10]，采用ANSYS有限元分析軟件對包含淺層排水系統(tǒng)的土壤柱進(jìn)行簡化并進(jìn)行數(shù)值仿真[11-12]，進(jìn)而分析土壤層中溶解性氮磷污染物的遷移過程及規(guī)律.

1 研究區(qū)域基本情況

Little Vermilion River流域位于美國伊利諾伊州東部中心位置，流域中心坐標(biāo)為(40°06′21.45″N, 87°41′34.12″W)，是本研究的對象區(qū)域.該地區(qū)的土壤類型包括粉砂壤土(Drummer)和粉壤土(Flangan)[2]，該地區(qū)特征實(shí)驗(yàn)點(diǎn)實(shí)測得到的初始條件如表1所示.讀取當(dāng)?shù)亟涤曛噩F(xiàn)期曲線，選取5 h降雨為基準(zhǔn)，降雨強(qiáng)度為4 mm/h，同時(shí)保持該強(qiáng)度，將降雨歷時(shí)提高到10～48 h，觀察樣方中氮磷遷移變化.

表1 遷移過程初始條件與邊界條件Tab.1 Initial and boundary condition of the transport process

為簡化計(jì)算，對樣方進(jìn)行相應(yīng)假設(shè)，流化床的設(shè)計(jì)為重復(fù)對稱結(jié)構(gòu)，且在縱向延伸至無限遠(yuǎn)，所以，樣方可以簡化成切面二維模型.本文主要關(guān)注流化床內(nèi)氮磷遷移的質(zhì)量分?jǐn)?shù)變化及遷移速度，因此，將降雨過程簡化成初始流速，省略布水裝置，通過計(jì)算等效水力傳導(dǎo)率，簡化了土壤分層過程計(jì)算量.大大降低計(jì)算難度，節(jié)省計(jì)算時(shí)間.出口處為壓力出口，流體(進(jìn)料、進(jìn)水)進(jìn)口處為速度入口，其余默認(rèn)設(shè)為穿孔墻邊界條件.

2 簡化淺層排水土壤柱結(jié)構(gòu)

為了有效計(jì)算垂直方向遷移過程，利用有限元分析軟件ANSYS?對選取的點(diǎn)遷移樣方進(jìn)行網(wǎng)格化，并通過求解邊界條件的數(shù)值分析方法對網(wǎng)格內(nèi)“三傳”過程進(jìn)行求解從而得到精確的遷移過程描述.首先，將土壤樣方簡化成土壤顆粒流化床，其中土壤粒徑為平均粒徑，流化床的設(shè)計(jì)如圖1所示.

圖1 土壤流化床設(shè)計(jì)示意Fig.1 Sketch of fluidized bed

3 結(jié)果與分析

通過實(shí)測數(shù)據(jù)知土壤表層的硝酸鹽、磷酸鹽質(zhì)量分?jǐn)?shù)分別為20和0.3 mg/kg，進(jìn)一步對樣方進(jìn)行近似，硝酸鹽與磷酸鹽均分布在表層20 cm土壤層內(nèi)，池容量為無限大.首先研究5 h降雨過程中土壤層中硝酸鹽質(zhì)量分?jǐn)?shù)的變化，由圖2可以看出，降雨5 h后，排水管線中并沒有滲流水分流入及硝酸鹽質(zhì)量分?jǐn)?shù)變化發(fā)生.

但在地表到淺層排水管處，質(zhì)量分?jǐn)?shù)已經(jīng)發(fā)生變化，質(zhì)量分?jǐn)?shù)分層過程也符合水位降曲線.由于土壤等效水力傳導(dǎo)率低于降雨速率，對遷移過程造成額外壓力，水流趨向于在淺層排水管上方形成快速滲透區(qū).這也是淺層排水能夠快速排出累積在農(nóng)田的多余水分的原因，同時(shí)，由于該驅(qū)動(dòng)力的存在，土壤中硝酸鹽也產(chǎn)生了濃差，如圖2所示.在同一水平線5 h降雨后最大濃差2 mg/kg，但5 h降雨及4 mm/h降雨強(qiáng)度所產(chǎn)生的驅(qū)動(dòng)力依然不足以使得水分到達(dá)排水管，水分只向下運(yùn)動(dòng)了33.5 cm.由于土壤也含有水分，界面質(zhì)量分?jǐn)?shù)差則是傳質(zhì)過程驅(qū)動(dòng)力，傳質(zhì)速率由濃差及傳質(zhì)系數(shù)決定，但此時(shí)降雨停止，減少了壓力，水分的運(yùn)動(dòng)只能通過重力作用，遷移過程也將減慢.

圖2 降雨5 h硝酸鹽 5 h后點(diǎn)垂直方向質(zhì)量分?jǐn)?shù)變化

Fig.2 Vertical nitrate concentration change in soil after 5 h rainfall event

5 h降雨停止后，隨著重力作用，水分繼續(xù)向下遷移，24 h后硝酸鹽開始出現(xiàn)在排水管中，其質(zhì)量分?jǐn)?shù)為3.99～4.99 mg/kg.由于有坡度存在，遠(yuǎn)處依然有排水管沒有水分流入，但只是理論上的情況，實(shí)際情況中并沒有無限遠(yuǎn)的管道存在，因此，管體將全部接收排水.由圖3可以看出，24 h后，由于降雨歷時(shí)只有5 h，質(zhì)量分?jǐn)?shù)差有明顯的分層，因?yàn)樵?4 h內(nèi)隨著水分的重力下滲并伴隨著水分被不飽和土壤吸收，水分向下遷移的速度不斷減慢，這為傳質(zhì)過程提供了足夠的時(shí)間和空間，使得質(zhì)量分?jǐn)?shù)分層現(xiàn)象明顯，而且質(zhì)量分?jǐn)?shù)分布也與地下水位的趨勢相同，在具有淺層排水管的位置水位出現(xiàn)下降，而在無排水管位置層呈現(xiàn)均勻分布趨勢，因此，在兩條管線直接地下水位呈現(xiàn)“拱形”.而質(zhì)量分?jǐn)?shù)的分布也與水分分布規(guī)律相同.

當(dāng)整個(gè)過程持續(xù)48 h后，硝酸鹽質(zhì)量分?jǐn)?shù)分布如圖4所示.隨著時(shí)間的推移，由于無持續(xù)水分的輸入，上層土壤質(zhì)量分?jǐn)?shù)無進(jìn)一步的擴(kuò)散變化，下層土壤隨水分重力下滲，使得質(zhì)量分?jǐn)?shù)分層被壓縮，但由于水量小，壓縮效果并不明顯.質(zhì)量分?jǐn)?shù)分布基本與24 h后接近，并無顯著變化，而管內(nèi)質(zhì)量分?jǐn)?shù)為4～5 mg/kg.

圖3 降雨5 h硝酸鹽24 h后點(diǎn)垂直方向質(zhì)量分?jǐn)?shù)變化

Fig.3 Vertical nitrate concentration change in soil after 24 h rainfall event

圖4 降雨5 h硝酸鹽48 h后點(diǎn)垂直方向質(zhì)量分?jǐn)?shù)變化

Fig.4 Vertical nitrate concentration change in soil after 48 h rainfall event

在5 h降雨結(jié)束后，磷的質(zhì)量分?jǐn)?shù)分布趨勢與硝酸鹽類似，因?yàn)槿芙鈶B(tài)磷酸鹽在遷移過程上與溶解態(tài)硝酸鹽并無差異.在初始質(zhì)量分?jǐn)?shù)為非飽和溶液的條件下，磷酸鹽在降雨發(fā)生后主要發(fā)生的遷移為傳質(zhì)過程，其中溶解過程是限制步驟，由圖5可以看出，磷酸鹽在降雨5 h后的質(zhì)量分?jǐn)?shù)分布，相對于硝酸鹽分布較密集，這是由溶解速度慢造成的，而質(zhì)量分?jǐn)?shù)降低的位置依然與硝酸鹽類似，出現(xiàn)在淺層排水管設(shè)置位置，但由于磷酸鹽初始質(zhì)量分?jǐn)?shù)較低，在圖內(nèi)顏色差異較大，但實(shí)際最大濃差僅為0.01 mg/kg.但該種分析過程并不能體現(xiàn)溶解過程，對于磷酸鹽這種微溶的物質(zhì)，傳質(zhì)現(xiàn)象在短時(shí)間內(nèi)并不明顯.水分也同樣沒有到管，水分到管時(shí)間理論上與硝酸鹽相同.當(dāng)水分到管時(shí)，磷酸鹽的濃差將比硝酸鹽更大，也是由溶解度低造成的.

如圖6，經(jīng)過24 h的傳質(zhì)過程后，濃差的分層加大，最大濃差為0.2 mg/kg，管內(nèi)質(zhì)量分?jǐn)?shù)為0.03 mg/kg.經(jīng)過48 h的滲透過程后，磷酸鹽的質(zhì)量分?jǐn)?shù)分布基本與24 h后類似，只是濃差進(jìn)一步擴(kuò)大，這與硝酸鹽遷移過程類似.在無進(jìn)一步的水分輸入后，水分通過重力作用向下運(yùn)動(dòng)，但速度會(huì)不斷減慢，這是由于下層土壤的不飽和度高于上層，土壤的張力更大，吸水性更強(qiáng)，使得水分更不易于向下遷移，從而濃差開始被壓縮.但由于初始質(zhì)量分?jǐn)?shù)不高，造成傳質(zhì)速率不高，從而導(dǎo)致濃差壓縮過程不明顯，依然保持與24 h后的狀態(tài)接近，從而隨著所有能夠到管并排走的水分都排走后，土壤中磷酸鹽的遷移將停滯在無水分運(yùn)動(dòng)之后.同時(shí)伴隨轉(zhuǎn)化過程的發(fā)生，在遷移過程中雖然也伴隨著轉(zhuǎn)化，但由于時(shí)間較短，轉(zhuǎn)化所造成的影響相對較小，但在遷移過程停滯后，轉(zhuǎn)化起到主要作用，其中包括土壤中有效磷與穩(wěn)態(tài)磷的動(dòng)態(tài)平衡過程，有效磷和溶解態(tài)磷之間的動(dòng)態(tài)平衡過程，磷的各個(gè)形態(tài)的相互轉(zhuǎn)化，土壤中微生物吸收作用等，以上作用導(dǎo)致了磷酸鹽遷移的停滯，但轉(zhuǎn)化不包括在本研究的討論中.

圖5 降雨5 h磷酸鹽5 h后點(diǎn)垂直方向質(zhì)量分?jǐn)?shù)變化

Fig.5 Vertical phosphate concentration change in soil after 5 h rainfall event

圖6 降雨5 h磷酸鹽24 h后點(diǎn)垂直方向質(zhì)量分?jǐn)?shù)變化

Fig.6 Vertical phosphate concentration change in soil after 24 h rainfall event

圖7 降雨10 h硝酸鹽10、24、48 h后點(diǎn)垂直方向質(zhì)量分?jǐn)?shù)變化

Fig.7 Vertical nitrate concentration change in soil after 10, 24, 48 h for 10 h rainfall event

進(jìn)一步分析圖7(a)、(b)、(c)，在10 h降雨結(jié)束后，管內(nèi)已經(jīng)出現(xiàn)質(zhì)量分?jǐn)?shù)變化，管內(nèi)的質(zhì)量分?jǐn)?shù)為5～6 mg/kg，此時(shí)最大濃差達(dá)8～9 mg/kg.在沒有達(dá)到最大質(zhì)量分?jǐn)?shù)的土壤層中，濃差也明顯高于5 h降雨過程，這是由降雨量不同直接造成的結(jié)果.流量直接影響了硝酸鹽的遷移，而且流量越大加速了硝酸鹽的遷移.流量的增加與硝酸鹽的遷移并不是線性關(guān)系，因?yàn)橄跛猁}的易溶性，在流失量的層面上限制條件為流失通量，在降雨量增大后，飽和土壤的透水性要明顯大于不飽和土壤，因?yàn)樵谂趴胀寥烂?xì)孔內(nèi)空氣后，同時(shí)由于淺層排水管的作用，飽和土壤的水滲透只是取決于土壤粒徑，土壤粒徑越大透水速度越快，因此，在土壤飽和后硝酸鹽的流失通量將達(dá)最大值.

24 h后管內(nèi)質(zhì)量分?jǐn)?shù)持續(xù)增大，達(dá)8～9 mg/kg，最大濃差為13～14 mg/kg，因?yàn)榻涤贻^之前增加1倍，增加了硝酸鹽的遷移.圖7(c)為48 h后質(zhì)量分?jǐn)?shù)分布，可以看出，管內(nèi)質(zhì)量分?jǐn)?shù)已達(dá)14～15 mg/kg，最大濃差也達(dá)15～16 mg/kg，通過觀察動(dòng)態(tài)變化過程，在30 h后管內(nèi)質(zhì)量分?jǐn)?shù)分布基本穩(wěn)定在該狀態(tài)，并持續(xù)至48 h.

對于10 h降雨過程磷酸鹽的遷移過程如圖8(a)、(b)、(c)所示，其質(zhì)量分?jǐn)?shù)變化過程與硝酸鹽并無顯著差異，變化趨勢類似，主要差別在于磷酸鹽遷移過程的限制步驟為溶解及傳質(zhì)過程，溶解度低且濃差低是造成磷酸鹽遷移慢的主要原因，但如果水量充分，時(shí)間足夠，在最終分布形式上磷酸鹽的樣式應(yīng)與硝酸鹽相同.

圖8 降雨10 h磷酸鹽10、24、48 h后點(diǎn)垂直方向質(zhì)量分?jǐn)?shù)變化

Fig.8 Vertical phosphate concentration change in soil after 10, 24, 48 h for 10 h rainfall event

進(jìn)一步觀察48 h連續(xù)降雨過程后氮磷遷移過程，圖9為硝酸鹽在48 h、4 mm/h降雨過程中，0、5、24、48 h后質(zhì)量分?jǐn)?shù)的分布狀況.可以看出，在連續(xù)降雨過程中，硝酸鹽的遷移隨時(shí)間不斷增大，直至達(dá)到最大質(zhì)量分?jǐn)?shù)，在24 h內(nèi)已經(jīng)有最大質(zhì)量分?jǐn)?shù)在管內(nèi)出現(xiàn)，說明硝酸鹽的遷移過程直接與流量相關(guān)，流量越大遷移量越大，遷移過程直接符合水分遷移過程.

為了體現(xiàn)硝酸鹽與磷酸鹽在垂直方向上遷移過程的不同,將質(zhì)量分?jǐn)?shù)分布標(biāo)尺設(shè)置與硝酸鹽一致(以大質(zhì)量分?jǐn)?shù)范圍為準(zhǔn)).由圖10可以看出，在0、5、24、48 h的遷移過程中，磷酸鹽在質(zhì)量分?jǐn)?shù)層面的變化并不顯著，但溶解的磷酸鹽的遷移也隨同水分遷移規(guī)律，磷酸鹽的垂直遷移過程不只取決于流量大小.其遷移也由降雨時(shí)間、初始池溶解態(tài)磷質(zhì)量分?jǐn)?shù)、有效磷含量、施肥時(shí)間等多種因素決定，因此，磷隨水分流失過程是磷遷移途徑中不顯著的一部分.正是因?yàn)榱椎牟灰兹苄?，多以穩(wěn)定態(tài)存在，大量磷的遷移則是由吸附在土壤顆粒的多種磷形式隨泥沙遷移過程中產(chǎn)生，并伴隨吸附解析、溶解溶出、水解化合等作用，最終進(jìn)入水體，因此，本研究只關(guān)注溶解態(tài)的磷遷移過程.進(jìn)一步討論在淺層排水管中，兩種物質(zhì)的質(zhì)量分?jǐn)?shù)隨時(shí)間變化過程，結(jié)果如圖11所示.在持續(xù)降雨過程中，硝酸鹽最快在6 h后到達(dá)管內(nèi)，6～18 h在持續(xù)的水分輸入下，質(zhì)量分?jǐn)?shù)呈指數(shù)增長趨勢，18 h后質(zhì)量分?jǐn)?shù)達(dá)到穩(wěn)定增長狀態(tài)，并超過20 mg/kg，48 h計(jì)算結(jié)束.如果繼續(xù)計(jì)算，沒有排出的多余水分在重力作用下將持續(xù)滲出，質(zhì)量分?jǐn)?shù)會(huì)在某一時(shí)刻達(dá)到峰值，然后再不斷隨水量減少，傳質(zhì)過程的驅(qū)動(dòng)力減弱，入管質(zhì)量分?jǐn)?shù)不斷下降，最終停止.而磷酸鹽的遷移能力相對硝酸鹽低得多，首先，活性磷在土壤中的比例相對較低，再加上溶解過程緩慢等原因使得磷酸鹽在垂直遷移過程不顯著.

圖9 48 h連續(xù)降雨中硝酸鹽遷移過程Fig.9 Nitrate transport pattern after 48 h continuous rainfall event

圖10 48 h連續(xù)降雨過程中磷酸鹽遷移過程Fig.10 Nitrate transport pattern after 48 h continuous rainfall event

圖11 管內(nèi)氮磷質(zhì)量分?jǐn)?shù)隨時(shí)間變化

Fig.11 Nitrate and phosphate concentration variation along with the time

4 結(jié) 論

1)在空間遷移過程點(diǎn)層面垂直遷移途徑中，在短時(shí)間降雨事件，氮磷隨水分滲透過程的形式一致，但表現(xiàn)出不同的特征.氮遷移的主要形式為硝酸鹽形態(tài)，硝酸鹽溶解度大，在雨量不斷增大過程中，硝酸鹽遷移量呈指數(shù)增長趨勢，直至水量輸入輸出平穩(wěn)后，硝酸鹽遷移量也達(dá)到平穩(wěn).

2)磷的遷移形態(tài)為溶解態(tài)磷酸鹽，磷酸鹽的溶解度低，磷各種形態(tài)的平均吸附系數(shù)高，因此，需要一定的水量才能打破磷酸鹽的水解平衡態(tài)，導(dǎo)致磷酸鹽遷移量及速度較硝酸鹽慢.

3)淺層排水管網(wǎng)中，氮磷質(zhì)量分?jǐn)?shù)的變化與土壤層中的變化過程相近，但在水分沒有滲透到排水管網(wǎng)土層時(shí)，管網(wǎng)內(nèi)氮磷質(zhì)量分?jǐn)?shù)并沒有變化，在降雨量適中時(shí)，氮磷的質(zhì)量分?jǐn)?shù)在管網(wǎng)中都呈線性增加的趨勢，而當(dāng)水量進(jìn)一步增大時(shí)，氮的質(zhì)量分?jǐn)?shù)呈現(xiàn)指數(shù)增加，而磷的流失質(zhì)量分?jǐn)?shù)依然為線性增加.

[1] ALEXANDER R B, SMITH R A, SCHWARZ G E, et al. Differences in phosphorus and nitrogen delivery to the gulf of Mexico from the Mississippi river basin[J]. Environ Sci Technol, 2008,42(3): 822-830.

[2] KALITA P K, COOKE R A C, ANDERSON S M, et al. Subsurface drainage and water quality: The Illinois experience[J]. Trans ASABE, 2007,50(5): 1651-1656.

[3] BROWN R A, SKAGGS R W, HUNT W F. Calibration and validation of DRAINMOD to model bioretention hydrology[J]. J Hydrol, 2013,486: 430-442.

[4] ALE S, GOWDA P H, MULLA D J, et al. Comparison of the performances of DRAINMOD-NII and ADAPT models in simulating nitrate losses from subsurface drainage systems[J]. Agric Water Manage, 2013,129: 21-30.

[5] TIAN S Y, YOUSSEF M A, SKAGGS R W, et al. DRAINMOD-FOREST: Integrated modeling of hydrology, soil carbon and nitrogen dynamics, and plant growth for drained forests[J]. J Environ Qual, 2012,41(3): 764-782.

[6] SKAGGS R W, YOUSSEF M A, CHESCHEIR G M. A drainmod-based method to estimate effects of drainage water management on annual nitrogen loss to surface water[J]. Trans ASABE, 2012,55(3): 799-808.

[7] MOHSIN R, MAJID Z A, YUSOF M Z. Safety distance between underground natural gas and water pipeline facilities[J]. Reliab Eng Syst Safe, 2014,131: 53-60.

[8] DAI Z L, HUANG Y, CHENG H L, et al. 3D numerical modeling using smoothed particle hydrodynamics of flow-like landslide propagation triggered by the 2008 Wenchuan earthquake[J]. Eng Geol, 2014,180: 21-33.

[9] NGUYEN T P O, HELBLING D E, BERS K, et al. Genetic and metabolic analysis of the carbofuran catabolic pathway in Novosphingobium sp KN65.2[J]. Appl Microbiol Biot, 2014,98(19): 8235-8252.

[10]YANG X F, RICHMOND M C, SCHEIBE T D, et al. Flow partitioning in fully saturated soil aggregates[J]. Transport Porous Med, 2014,103(2): 295-314.

[11]HU D W, LI L, LIU H, et al. Design and control of rotating soil-like substrate plant-growing facility based on plant water requirement and computational fluid dynamics simulation[J]. Ecol Eng, 2014,64: 269-275.

[12]HONG S W, LEE I B, SEO I H, et al. Measurement and prediction of soil erosion in dry field using portable wind erosion tunnel[J]. Biosyst Eng, 2014,118: 68-82.

Numerical simulation of nitrogen and phosphorus vertical transportation in subsurface drainage area

LI Shiyang1, WANG Li1, Prasanta K Kalita2, MA Fang1

(1.State Key Laboratory of Urban Water Resource and Environment (Harbin Institute of Technology), Harbin 150090, China;2.Department of Agriculture and Biological Engineering，University of Illinois at Urbana-Champaign, Illinois 61802，USA)

Subsurface drainage system is an efficient system to drain extra flood out from the farmland. But this system also changed the transport pathway and speed of water and the dissolved pollution in farmland. To get a full understanding of vertical transportation of nitrate and phosphate in subsurface drainage agricultural watershed, finite element analysis technology was used to make the numerical simulation of mass transfer of soil column. Results show that there is no water and dissolved pollutant in pipe after 5 h rainfall with intensity of 4 mm/h, because the rainfall in 5 h is not enough to saturate the soil layer from surface to pipe. But if extending the rainfall duration to 10 h, the soil became saturated and then the nitrate and phosphate loss observed. The nitrate loss was extremely higher than phosphate (nitrate concentration is 9-10 mg/L，phosphate is 0.05-0.1 mg/L). After 24 h rainfall event, the nitrate loss reached to the maximum. The nitrate loss in tile water intensity was at an exponential growth but phosphate just increased at a liner growth. This phenomenon may be caused by the differences of material properties.

subsurface drainage system; non-point pollution; vertical transportation; finite element analysis; numerical simulation

10.11918/j.issn.0367-6234.2017.02.008

2015-01-22

國家自然科學(xué)基金(51179041);水體污染控制與治理科技重大專項(xiàng)(2013ZX07201007);城市水資源與水環(huán)境國家重點(diǎn)實(shí)驗(yàn)室自主課題(2014TS05)

李世陽(1987—)，男，博士研究生； 王 立(1978—)，女，副教授，博士生導(dǎo)師； 馬 放(1963—)，男，教授，博士生導(dǎo)師

馬 放，mafang@hit.edu.cn

X254

A

0367-6234(2017)02-0044-05