農(nóng)業(yè)退水分層利用裝置設(shè)計(jì)與數(shù)值模擬

孫志林 李眾徹 方詩標(biāo) 許 丹

(浙江大學(xué)港口海岸與近海工程研究所, 杭州 310058)

農(nóng)業(yè)退水分層利用裝置設(shè)計(jì)與數(shù)值模擬

孫志林 李眾徹 方詩標(biāo) 許 丹

(浙江大學(xué)港口海岸與近海工程研究所, 杭州 310058)

從理論可行性出發(fā)，初步設(shè)計(jì)一種分流裝置，再從功能和效率的角度對其進(jìn)行改進(jìn)，最后采用計(jì)算流體力學(xué)軟件Fluent對農(nóng)業(yè)退水分流過程進(jìn)行數(shù)值模擬，對比分析不同流速下的鹽度場和速度場，結(jié)果表明該分流裝置可為最終設(shè)計(jì)提供參考。試驗(yàn)工況分3組進(jìn)行，設(shè)定入口流速作為唯一變量，3組試驗(yàn)均假定相同的入口鹽度分布函數(shù)。截取2個(gè)Y軸剖面的流速場云圖和4個(gè)X軸剖面的鹽度場云圖進(jìn)行對比分析。結(jié)果表明，分流裝置可以對不同鹽度的農(nóng)業(yè)退水進(jìn)行分流。上導(dǎo)流槽得到低鹽度農(nóng)業(yè)退水，下導(dǎo)流槽得到較高鹽度退水。分流的效率取決于裝置尺寸、分層流流速和鹽度分布等因素，可根據(jù)實(shí)際情況進(jìn)行調(diào)整和優(yōu)化。

農(nóng)業(yè)退水； 分流裝置； 排堿渠； 數(shù)值模擬； 鹽堿地

引言

土地鹽堿化嚴(yán)重制約南疆地區(qū)的社會發(fā)展。農(nóng)業(yè)部門一般采取水利、生物、化學(xué)等治理措施改良鹽堿地。水利工程治理措施主要是修建水庫、排水渠道網(wǎng)絡(luò)，建立較為完善的灌溉系統(tǒng)和現(xiàn)代化排水系統(tǒng)[1-10]。生物治理措施主要是培育耐鹽植物，利用生物排水降低地下水位，減少地面蒸發(fā)，防止土地表面積鹽[11-15]。化學(xué)治理措施則是在鹽漬化土地中施加石膏、礦渣等改良劑，降低土地中的鹽堿含量[16-19]。迄今為止，南疆大面積改良土地鹽漬化的方式仍是以“排水洗鹽”為主。然而大面積推廣使用“排水洗鹽”會使未經(jīng)處理的鹽堿水隨意排放，不僅閑置有限的水資源，而且威脅周圍生態(tài)環(huán)境。因此，需要通過更為科學(xué)合理的技術(shù)措施將鹽堿水的礦化度降低，以減輕廢棄鹽堿水對環(huán)境的惡劣影響。

南疆地區(qū)的排堿干渠有著多個(gè)入水口，各個(gè)支排的退水鹽度高低不一，且常年處于太陽高溫輻射量下，具備形成分層流的條件。關(guān)于分層流的數(shù)值模擬，國內(nèi)學(xué)者已經(jīng)進(jìn)行了初步的探索，并取得一定成果。包蕓等[20]運(yùn)用三維斜壓模型對珠江口伶仃洋區(qū)域鹽度分層進(jìn)行了數(shù)值模擬，模擬結(jié)果與實(shí)測資料基本吻合。胡振紅等[21]建立了溫度和鹽度分層流的k-ε模型，進(jìn)行了數(shù)值模擬計(jì)算，計(jì)算結(jié)果正確地反映了分層及紊動特征，與實(shí)測結(jié)果吻合較好。

本文設(shè)計(jì)水體分流裝置[22]，該裝置可將分層水體以一種純物理方式分流。在此基礎(chǔ)上，進(jìn)一步改進(jìn)裝置，提高水體分離的精度并實(shí)現(xiàn)完全自動化。通過Fluent計(jì)算流體力學(xué)軟件對分層流體流經(jīng)分流裝置的過程進(jìn)行數(shù)值模擬，對比分析不同流速下的鹽度場和速度場，為裝置實(shí)際使用提供優(yōu)化建議。

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

1.1 控制方程

控制方程采用雷諾平均方程組(RANS)，連續(xù)方程與動量方程為笛卡爾坐標(biāo)系下的張量形式，湍流模型采用標(biāo)準(zhǔn)k-ε模型，鹽度計(jì)算采用組分輸運(yùn)方程[23]。

連續(xù)性方程為

(1)

式中ρ——流體密度t——時(shí)間ui——速度分量xi——笛卡爾坐標(biāo)系下對應(yīng)的方向

動量方程為

(2)

式中p——壓力u′i、u′j——速度脈動量μ——動力粘性系數(shù)δij——狄克拉函數(shù)uj——速度分量xj、xl——笛卡爾坐標(biāo)系下對應(yīng)的方向

k方程為

Gk+Gb-ρε-YM+Sk

(3)

其中

μt=ρCμk2ε-1Cμ=0.09

式中k——湍動能μt——渦粘性系數(shù)Gk——由層流速度梯度而產(chǎn)生的湍流動能Gb——浮力而產(chǎn)生的湍流動能YM——擴(kuò)散產(chǎn)生的波動Sk——用戶自定義源項(xiàng)σk——k方程的湍流Prandtl數(shù)，取1.0ε——湍流耗散率

ε方程為

(4)

其中C1ε=1.44C2ε=1.92C3ε=tanh|vu-1|

式中C1ε、C2ε、C3ε——常數(shù)項(xiàng)v——沿重力方向速度分量u——平行重力方向速度分量Sε——用戶自定義源項(xiàng)σε——ε方程的湍流Prandtl數(shù)，取1.32

鹽度c控制方程為

(5)

式中Γc——鹽度擴(kuò)散系數(shù)Sc——鹽度源項(xiàng)

1.2 邊界條件

入口邊界設(shè)置為速度入口，速度通過用戶自定義函數(shù)(User defined function，UDF)設(shè)置為充分發(fā)展流動狀態(tài)，湍動強(qiáng)度設(shè)置為5%，同時(shí)鹽度在入口處通過UDF按照設(shè)計(jì)濃度分別加入計(jì)算系統(tǒng)。

出口邊界為無回流自由出口。

在固體壁面附近的粘性底層內(nèi)用標(biāo)準(zhǔn)壁面函數(shù)來處理[24]。

1.3 數(shù)值方法

在笛卡爾坐標(biāo)系下，式(1)～(4)可寫成統(tǒng)一形式

(6)

式中φ——ui、k、εS——源項(xiàng)?！獜V義擴(kuò)散系數(shù)

采用有限控制體積法來離散計(jì)算區(qū)域，然后在每個(gè)控制體積中對微分方程進(jìn)行積分，再把積分方程線性化，得到各未知變量的代數(shù)方程組，最后求解方程組即可求出各未知變量。采用Simple算法對壓力和速度場進(jìn)行耦合計(jì)算。流場迭代穩(wěn)定后計(jì)算鹽度場，模擬中采用單向耦合，鹽度計(jì)算不影響流場計(jì)算。

2 模型設(shè)計(jì)與數(shù)值模擬

2.1 裝置的原理與改進(jìn)

初始分流裝置如圖1所示，該裝置包括左分流裝置和與左分流裝置倒置的右分流裝置。左分流裝置和右分流裝置固定連接且相互平行，固定安裝在河渠的橫截面上。因?yàn)楦×ψ饔茫蟾×η蚝陀腋×η蚍謩e帶動左閘門和右閘門上下平移，由于密度和鹽度呈正相關(guān)關(guān)系，鹽度低的水從左閘門上方流出，鹽度高的水從右閘門下方流出，從而實(shí)現(xiàn)排堿渠不同濃度水體的分流。該裝置用純物理方法分流，具有結(jié)構(gòu)簡單、造價(jià)低廉等特點(diǎn)，目前已獲授權(quán)。

圖1 初始分流裝置Fig.1 Initial shunting device

圖2 帶無線收發(fā)功能的分流裝置Fig.2 Shunting device with wireless transceiver function

初始分流裝置雖優(yōu)點(diǎn)明顯，但也存在分離精度不高、無法遠(yuǎn)程控制等不足，因此在最初的發(fā)明專利基礎(chǔ)上做了進(jìn)一步改進(jìn)，設(shè)計(jì)了一種帶無線收發(fā)功能的鹽水分流裝置[25]。如圖2所示，該裝置包括左分流裝置和右分流裝置。左分流裝置和右分流裝置固定連接且相互平行，固定安裝在排堿渠的橫截面上。因上、下層水體的鹽度不同，導(dǎo)致置于上、下導(dǎo)流槽口邊緣的傳感器產(chǎn)生不同的電信號，進(jìn)而驅(qū)使左右電動機(jī)箱中的電動機(jī)分別帶動左閘門和右閘門進(jìn)行相應(yīng)的上下平移，使鹽度低的水從左閘門上方流出，鹽度高的水從右閘門下方流出，實(shí)現(xiàn)不同濃度鹽水的分流。其后附帶的導(dǎo)流槽可將分離后的液體導(dǎo)流，上層鹽度較低的水可以直接投入農(nóng)業(yè)灌溉，下層較高濃度鹽水則可以曬鹽、提鹽，從而達(dá)到防止土地鹽堿化的目的。裝置左右均附帶無線收發(fā)器，具備實(shí)時(shí)可監(jiān)控以及突發(fā)情況下可人工干預(yù)的功能，能很好地適應(yīng)環(huán)境變化，達(dá)到完全自動化。

2.2 研究對象與試驗(yàn)工況

本文的研究對象為置于排堿渠的分流裝置，實(shí)際模型簡化后的計(jì)算模型如圖3所示，X軸方向是分層流體流動方向，Z軸方向垂直向上，Y軸方向和分層流體流動方向垂直。排堿渠寬4.0 m、深1.5 m、長12.0 m，其橫截面為矩形，在排堿渠X方向2.0 m處，垂直流體流動方向放置分流裝置，分流裝置計(jì)算模型簡化為一塊擋板，在擋板的左上方和右下方各開1個(gè)矩形導(dǎo)流槽口，2個(gè)導(dǎo)流槽口尺寸一致，寬2.0 m、高0.4 m，導(dǎo)流槽口后均連接長10.0 m的導(dǎo)流槽，上下導(dǎo)流槽底部為平坡。各導(dǎo)流槽末端連接水泵，如有需要，可為導(dǎo)流槽中水體的流動提供持續(xù)動力。

圖3 排堿渠放置分流裝置三視圖(單位：m)Fig.3 Three-view graphs of shunting device in alkali draining ditches

對流體計(jì)算域采用六面體網(wǎng)格劃分，考慮在壁面附近和擋板水流分離處流場參數(shù)變化梯度比較大，因此網(wǎng)格劃分時(shí)，此部分采用網(wǎng)格加密處理，在邊界第1層網(wǎng)格設(shè)置為0.005 m，網(wǎng)格擴(kuò)展比為1.2，經(jīng)過驗(yàn)證，設(shè)置符合壁面網(wǎng)格符合標(biāo)準(zhǔn)壁面函數(shù)使用條件。通過網(wǎng)格無關(guān)性驗(yàn)證。

模擬分3組試驗(yàn)工況，第1、2、3組的入口流速分別為u1=0.1 m/s，u2=0.2 m/s，u3=0.5 m/s，這3個(gè)流速為典型試驗(yàn)工況，分別模擬排堿渠中流體在枯水季節(jié)、一般季節(jié)和旺水季節(jié)3種情況下的流動。每組試驗(yàn)除流速條件不一樣外，其他條件均一致。入口鹽度的分布設(shè)置采用有關(guān)文獻(xiàn)的實(shí)測結(jié)果，取底部鹽度是表面鹽度的3倍[20]。分層流體鹽度變化函數(shù)較好地刻畫了鹽度漸進(jìn)變化規(guī)律，函數(shù)曲線如圖4所示。由函數(shù)關(guān)系式可知，水表面(X=0 m)鹽度為1.0 g/L；水底(X=1.5 m)鹽度為3.1 g/L。水底鹽度取3.0 g/L左右，基于絕大多數(shù)植物都不能在此濃度鹽水灌溉下存活的考慮，因此當(dāng)濃度接近3.0 g/L時(shí)，已經(jīng)不適合灌溉，而應(yīng)該分離后作他用。

圖4 入口鹽度分布函數(shù)Fig.4 Distribution function of salinity at inlet

3 結(jié)果與分析

圖5～11是3組不同工況下，分層流體達(dá)到穩(wěn)態(tài)時(shí)，流速和鹽度分布圖。其中，圖5～7為下導(dǎo)流槽(Y=1.0 m剖面)和上導(dǎo)流槽(Y=3.0 m剖面)在X軸方向上達(dá)到穩(wěn)態(tài)時(shí)的流速分布情況。圖8～11分別是X為2.1、4.0、6.0、8.0 m 4個(gè)不同剖面下的鹽度分布。

由圖5～7可看出，無論是u1=0.1 m/s，u2=0.2 m/s，u3=0.5 m/s的何種工況，在上、下導(dǎo)流槽入口處(也就是X=2.0 m處附近)，鹽水的流速快速增加，均達(dá)到各工況下入口流速的4倍左右，但鹽水在導(dǎo)流槽的后續(xù)流動中，速度基本沒有變化。這是由于上、下導(dǎo)流槽橫截面變小，且上、下導(dǎo)流槽都是管流，不存在鹽堿水溢出的可能。排堿渠的寬4.0 m、高1.5 m，橫截面積6.0 m2，而分流裝置上、下導(dǎo)流槽寬2.0 m、高0.4 m，橫截面積之和1.6 m2，因此根據(jù)流體的質(zhì)量守恒定律和連續(xù)性假設(shè)，流量不變，面積減小，必然導(dǎo)致流速增加。

在X=2.0 m處上部，可以明顯看出鹽水的流速快速下降為零，這是由于分流裝置的右截流板發(fā)揮了截流的作用，將分流裝置迎水面右邊濃度較低的稀鹽水截流，因此速度快速降為零。圖5～7的b圖,在X=2.0 m處流速快速降為零，原因和a圖類似，因?yàn)榉至餮b置下部的截流板將分流裝置迎水面左邊濃度較高的鹽水截流。

圖5 第1組 (u1=0.1 m/s)Y軸剖面的流速場Fig.5 Velocity fields of Y-axis transverse section for the first group (u1=0.1 m/s)

圖6 第2組 (u2=0.2 m/s)Y軸剖面的流速場Fig.6 Velocity fields of Y-axis transverse section for the second group (u2=0.2 m/s)

圖7 第3組 (u3=0.5 m/s)Y軸剖面的流速場Fig.7 Velocity fields of Y-axis transverse section for the third group (u3=0.5 m/s)

由圖8～11可以看出，分流裝置能實(shí)現(xiàn)分流的功能：上導(dǎo)流槽的平均鹽度低于下導(dǎo)流槽的平均鹽度，上導(dǎo)流槽分流出低鹽度的水，下導(dǎo)流槽分流出較高鹽度的水。流速越快的鹽水，其X剖面的鹽度分布趨于穩(wěn)定的時(shí)間也越長。

根據(jù)圖4所示鹽度變化函數(shù)，當(dāng)與水面的垂直距離為0.4 m時(shí)，對應(yīng)的鹽度為2.08 g/L。3組工況下，上導(dǎo)流槽X軸橫截面上鹽水的平均鹽度約等于上導(dǎo)流槽槽底對應(yīng)深度的鹽度。該現(xiàn)象說明，上導(dǎo)流槽分流后的鹽水鹽度分布基本取決于入口處的鹽度分布函數(shù)和上導(dǎo)流槽槽底對應(yīng)深度。因此，在排堿渠中使用分流裝置時(shí)，如果要提取低于某一鹽度的分層流體，則應(yīng)該調(diào)高分流裝置上導(dǎo)流槽，使上導(dǎo)流槽槽底對應(yīng)深度的鹽度低于該值，否則提取的分層流體容易鹽度過高，影響分流效果。

圖8 3組不同流速下，X=2.1 m剖面的鹽度場Fig.8 Salinity fields in transverse section of X=2.1 m under three different velocities

圖9 3組不同流速下，X=4.0 m剖面的鹽度場Fig.9 Salinity fields in transverse section of X=4.0 m under three different velocities

相同的X軸截面，不同流速對比可以看出，u1=0.1 m/s工況的分流效果要優(yōu)于u2=0.2 m/s和u3=0.5 m/s工況，進(jìn)一步觀察可以發(fā)現(xiàn)u2=0.2 m/s和u3=0.5 m/s工況之間的分流效果差別不大。該現(xiàn)象說明，為了使上導(dǎo)流槽盡可能分流出濃度較低的鹽水，在實(shí)際的排堿渠中使用分流裝置時(shí)，應(yīng)當(dāng)適當(dāng)調(diào)低上導(dǎo)流槽的流速，因?yàn)榱黧w流速越大越容易從層流向紊流發(fā)展，這會導(dǎo)致流體質(zhì)點(diǎn)的運(yùn)動極不規(guī)則，由于脈動的急劇混摻，流體動量、能量、溫度以及含有物的濃度擴(kuò)散速率都會增大。這樣的流動會破壞原來的鹽度分層現(xiàn)象，使得分流裝置的分流作用減弱。

圖10 3組不同流速下，X=6.0 m剖面的鹽度場Fig.10 Salinity fields in transverse section of X=6.0 m under three different velocities

圖11 3組不同流速下，X=8.0 m剖面的鹽度場Fig.11 Salinity fields in transverse section of X=8.0 m under three different velocities

下導(dǎo)流槽的鹽水流動速度沿程均比較穩(wěn)定，變化不明顯，但是可以看出，隨著鹽水的流動，3個(gè)工況下的下導(dǎo)流槽X截面的鹽度均有所下降。分流裝置下導(dǎo)流槽設(shè)置的目的是為了分流較高濃度的鹽水，若鹽水濃度變小，則會與設(shè)想的結(jié)果背道而馳。進(jìn)一步深入比較各個(gè)工況下導(dǎo)流槽的鹽度分布可以看出，u1=0.1 m/s工況的下導(dǎo)流槽鹽度分布比u2=0.2 m/s和u3=0.5 m/s工況要穩(wěn)定，沿程流動中，鹽水濃度更不易變小。其原因和上導(dǎo)流槽類似，因?yàn)榱黧w流速越大，液體質(zhì)點(diǎn)作不規(guī)則運(yùn)動，互相混摻，軌跡越曲折混亂，下方較高濃度的鹽水和上方低濃度鹽水摻混，致使下方導(dǎo)流槽分流出的鹽水濃度變低。為了提高下導(dǎo)流槽的分流效率，可以減小下導(dǎo)流槽橫截面的寬度。

4 結(jié)論

(1)分流裝置上導(dǎo)流槽旨在分流低鹽度的水體，上導(dǎo)流槽分流后的退水濃度很大部分取決于入口處鹽度分布和上導(dǎo)流槽槽底對應(yīng)的深度。因此，在排堿渠中使用分流裝置時(shí)，若想充分降低上導(dǎo)流槽分流出的鹽水濃度，在實(shí)際工作條件允許的情況下，可以選擇抬高上導(dǎo)流槽的槽底。

(2)降低下導(dǎo)流槽槽底可以提高下導(dǎo)流槽的分流效果，使得其分流的鹽水濃度得到提高，但是這樣會使得下導(dǎo)流槽的橫截面積變小，影響截面的流量，降低分流效率，在實(shí)際使用中需要綜合考慮。

(3)無論對于上導(dǎo)流槽還是下導(dǎo)流槽，降低退水流速都是提高分流效果的有效手段，因此入口處分層鹽水流動速度不應(yīng)過快，否則退水會在流動中快速摻混，分流裝置達(dá)不到分流的目的。但是入口處流速過低，則會導(dǎo)致分流的效率降低。故在實(shí)際使用過程中，需要在分流效率和分流效果方面做一定的權(quán)衡。

1 張珂萌, 牛文全, 汪有科, 等. 微咸水微潤灌溉下土壤水鹽運(yùn)移特性研究[J/OL]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2017, 48(1):175-182.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20170123&flag=1. DOI:10. 6041/j.issn.1000-1298.2017.01.023. ZHANG Kemeng, NIU Wenquan, WANG Youke, et al. Characteristics of water and salt movement in soil under moistube-irrigation with brackish water[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2017, 48(1): 175-182. (in Chinese)

2 BAHCECI I, NACAR A S. Subsurface drainage and salt leaching in irrigated land in southeast Turkey[J]. Irrigation and Drainage, 2009, 58(3): 346-356.

3 姚寶林, 李光永, 葉含春, 等. 干旱綠洲區(qū)膜下滴灌棉田土壤鹽分時(shí)空變化特征研究[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào), 2016, 47(1): 151-161. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20160120&flag=1. DOI：10.6041/j.issn.1000-1298.2016.01.020. YAO Baolin, LI Guangyong, YE Hanchun, et al. Characteristics of spatial and temporal changes in soil salt content in cotton fields under mulched drip irrigation in arid oasis regions[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(1): 151-161. (in Chinese)

4 RITZEMA H P, SATYANARAYANA T V, RAMAN S, et al. Subsurface drainage to combat waterlogging and salinity in irrigated lands in India:lessons learned in farmers’fields[J]. Agricultural Water Management, 2008, 95(3): 179-189.

5 NAVARRO J M, GARRIDO C, FLORES P, et al. The effect of salinity on yield and fruit quality of pepper grown in perlite[J]. Spanish Journal of Agricultural Research, 2010, 8(1): 142-150.

6 郭麗, 鄭春蓮, 曹彩云, 等. 長期咸水灌溉對小麥光合特性與土壤鹽分的影響[J/OL]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2017, 48(1): 183-190. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx? file_no=20170124&flag=1. DOI:10.6041/j.issn.1000-1298.2017.01.024. GUO Li, ZHENG Chunlian, CAO Caiyun, et al. Effect of long-term saline water irrigation on photosynthetic characteristics of winter wheat and soil salt content[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2017, 48(1): 183-190. (in Chinese)

7 王振華, 裴磊, 鄭旭榮, 等. 鹽堿地滴灌春小麥光合特性與耐鹽指標(biāo)研究[J/OL]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2016, 47(4): 65-72, 104. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20160410&flag=1. DOI: 10.6041/j.issn.1000-1298.2016.04.010. WANG Zhenhua, PEI Lei, ZHENG Xurong, et al. Photosynthetic characteristics and salt-tolerance indexes of spring wheat with drip irrigation in saline-alkali soils[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(4): 65-72, 104. (in Chinese)

8 QADIR M, GHAFOOR A, MURTAZA G. Amelioration strategies for saline soils: a review[J]. Land Degradation and Development, 2000, 11(6): 501-521.

9 KELLENERS T J, CHAUDHRY M R. Drainage water salinity of tubewells and pipe drains:a case study from Pakistan[J].Agricultural Water Management, 1998, 37(1): 41-53.

10 RAO K V G K, LEEDS-HARRISON P B. Desalinization with subsurface drainage[J]. Agricultural Water Management, 1991, 19(4): 303-311.

11 AKHTER J, MURRAY R, MAHMOOD K, et al. Improvement of degraded physical of a saline-sodic soil by reclamation with kallar grass(Leptochloafusca)[J]. Plant and Soil, 2004,258(1): 207-216.

12 RAVINDRAN K C, VENKATESAN K, BALAKRISHNAN V, et al. Restoration of saline land by halophytes for Indian soils[J]. Soil Biology & Biochemistry, 2007, 39(10): 2661-2664.

13 REDDY M P, SHAH M T, PATOLIA J S.Salvadorapersica, a potential species for industrial oil production in semiarid saline and alkali soils[J]. Industrial Crops and Products, 2008, 28(3): 273-278.

14 QADIR M, STEFFENS D, YAN F, et al. Sodium removal from a calcareous saline-sodic soil through leaching and plant uptake during phytoremediation[J]. Land Degradation & Development,2003, 14(3): 301-307.

15 MADSEN P A, MULLIGAN D R. Effect of NaCl on emergence and growth of a range of provenances ofEucalyptusciriodora,Eucalyptuspopulnea,EucalyptuscamaldulensisandAcaciasalicina[J]. Forest Ecology Management, 2006, 228(1-3): 152-159.

16 楊軍，孫兆軍，劉吉利，等.脫硫石膏糠醛渣對新墾龜裂堿土的改良洗鹽效果[J].農(nóng)業(yè)工程學(xué)報(bào)，2015， 31(17):128-135. YANG Jun，SUN Zhaojun，LIU Jili，et al. Effects of saline improvement and leaching of desulphurized gypsum combined with furfural residue in newly reclaimed farmland crack alkaline soil[J].Transactions of the CSAE，2015，31(17):128-135. (in Chinese)

17 吳軍虎，陶汪海，王海洋，等.羧甲基纖維素鈉對土壤團(tuán)粒結(jié)構(gòu)及水分運(yùn)動特性的影響[J].農(nóng)業(yè)工程學(xué)報(bào)，2015, 31(2):117-123. WU Junhu, TAO Wanghai, WANG Haiyang，et al. Influence of sodium carboxyl methyl cellulose on soil aggregate structure and soil water movement[J].Transactions of the CSAE，2015，31(2):117-123. (in Chinese)

18 張江輝，白云崗，張勝江，等.兩種化學(xué)改良劑對鹽漬化土壤作用機(jī)制及對棉花生長的影響[J].干旱區(qū)研究，2011, 28(3):384-388. ZHANG Jianghui, BAI Yungang, ZHANG Shengjiang，et al. Effect of two kinds of soil ameliorant on saline soil improvement and cotton growth[J].Arid Zone Research，2011, 28(3):384-388. (in Chinese)

19 王金滿，楊培嶺，張建國，等. 脫硫石膏改良堿化土壤過程中的向日葵苗期鹽響應(yīng)研究[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2005, 21(9):33-37. WANG Jinman, YANG Peiling， ZHANG Jianguo，et al. Salinity effect on sunflower at seedling stage during improving sodic soils reclaimed with by-product from flue gas desulphurization(BFGD)[J].Transactions of the CSAE，2005, 21(9):33-37. (in Chinese)

20 包蕓，任杰. 伶仃洋鹽度高度層化現(xiàn)象及鹽度鋒面的研究[J].水動力學(xué)研究與進(jìn)展, 2005, 20(6):689-693. BAO Yun, REN Jie. Numerical simulation of high resolution on the phenomenon of saline stratification in Lindinyang[J]. Chinese Journal of Hydrodynamics, 2005, 20(6):689-693. (in Chinese)

21 胡振紅, 沈永明, 鄭永紅,等. 溫度和鹽度分層流的數(shù)值模擬[J]. 水科學(xué)進(jìn)展, 2001, 12(4):439-444. HU Zhenhong, SHEN Yongming, ZHENG Yonghong，et al. Numerical simulation of thermal and saline stratified flow[J]. Advances in Water Science, 2001, 12(4):439-444. (in Chinese)

22 孫志林，李眾徹，方詩標(biāo)，等. 排鹽溝咸水分流裝置:中國, 201510392609.2[P].2015-07-02.

23 CHAN L，CHIN C，SORIA J，et al. Large eddy simulation and Reynolds-averaged Navier-Stokes calculations of supersonic impinging jets at varying nozzle-to-wall distances and impinging angles[J]. International Journal of Heat and Fluid Flow, 2014, 47(3):31-41.

24 TIAN L，AHMADI G. Particle deposition in turbulent duct flows—comparisons of different model predictions[J]. Journal of Aerosol Science, 2007, 38(4): 377-397.

25 孫志林, 李眾徹, 方詩標(biāo), 等. 帶無線收發(fā)功能的鹽水分流裝置:中國,201511018244.3[P]. 2015-12-30.

Device Design and Numerical Simulation of Agricultural Wastewater Layered Utilization

SUN Zhilin LI Zhongche FANG Shibiao XU Dan

(InstituteofPort,CoastalandOffshoreEngineering,ZhejiangUniversity,Hangzhou310058,China)

The rainfall in the south region of Xinjiang is rare which caused severe problems of soil salinization. The salt leaching of farmland and wastewater draining from irrigating farmland are rich in sodium which indirectly leads to the ecological pollution. To effectively deal with the excessive salinity of agricultural wastewater from alkali draining ditches, a shunting device was invented. Firstly, the feasibility of the theory to design the preliminary device was introduced. Secondly, improvement of the device in its functionality and efficiency was given. Finally, numerical simulation of the streaming process of agricultural wastewater based on the computational fluid dynamics utilizing the software fluent was obtained. The derived results were used for the ultimate design of the shunting device. The test conditions were categorized into three groups. The inlet flow rate was set as the only variable, and the spatial distribution of inlet salinity among three groups was assumed to be identical. TwoY-axis sectional velocity field cloud diagrams and fourX-axis sectional salinity field cloud diagrams were selected for detailed comparison. It showed that the shunting device can separate water based on salinity distribution, and agricultural wastewater with low salinity was obtained from the lower diversion trench, vice versa. The efficiency of the shunting device depended on its size, stratification flow velocity and salinity distribution. Adjustment and optimization should be subjected to the actual situation.

agricultural wastewater； shunting device； alkali draining ditches； numerical simulation； saline land

10.6041/j.issn.1000-1298.2017.06.026

2016-09-23

2017-01-20

高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金項(xiàng)目(2120101110108)、中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目(2015XZZX002-01)和水體污染控制與治理科技重大專項(xiàng)(2009ZX07424-001)

孫志林(1956—)，男，教授，博士生導(dǎo)師，主要從事水沙動力學(xué)、河口海岸數(shù)值模擬和水資源利用研究，E-mail： oceansun@zju.edu.cn

S24; TV93

A

1000-1298(2017)06-0203-07