王 勇，張建豐，吳文勇，趙永安，萬 發(fā) ，廖人寬

(1.西安理工大學(xué)西北旱區(qū)生態(tài)水利工程國家重點實驗室培育基地，西安 710048； 2.中國水利水電科學(xué)研究院水利研究所，北京 100044;3. 運城市尊村引黃灌溉管理局，山西 運城 044000)

0 引 言

在灌溉作物生產(chǎn)系統(tǒng)中廣泛使用表施肥料,目前我國采用地面灌溉方式的面積占98%[1]，畦灌由于簡單易行成為最主要的地面灌水方式之一[2]，且撒施和液施是2種最常用的畦灌施肥方式，其中，撒施是指先將肥料撒在地表，然后灌溉，肥料隨水運動和下滲到土壤中，達到滿足作物對養(yǎng)分的需求方法；而液施是指先將肥料完全溶解，然后注入灌溉水中，達到灌溉施肥目的的方法[3]。黃河流域的氣候條件和水資源狀況決定了農(nóng)業(yè)發(fā)展依賴于灌溉[4]，但由于引黃灌區(qū)田間工程設(shè)施不完善、畦田規(guī)格不合理以及管理粗放等原因?qū)е滤Y源浪費嚴(yán)重[5]，黃河流域灌區(qū)現(xiàn)狀平均灌溉水利用系數(shù)0.49，低于全國平均水平[6]；同時氮肥的利用效率僅為0.35，由此帶來一系列農(nóng)業(yè)管理和環(huán)境問題[7]。在地表水流中水流運動是溶質(zhì)遷移的水力學(xué)基礎(chǔ)[8]，畦長和入畦單寬流量作為畦灌的基本參數(shù)不僅決定著灌水均勻度和水分利用效率[2]，而且對畦灌施肥效果有至關(guān)重要的影響[9]。在我國大部分灌區(qū)的畦長在30～100 m，入畦單寬流量在1～4 L/(s·m)[2]。而引黃灌區(qū)畦灌具有畦長較長、入畦流量大的特點，畦長普遍在150 m以上，甚至有的超過300 m，入畦單寬流量有的達到9.26 L/(s·m)[10,11]。且作為灌溉水源的黃河水中含有大量懸浮顆粒物和營養(yǎng)物質(zhì)，懸浮顆粒物的存在對水中物質(zhì)的遷移轉(zhuǎn)化過程影響巨大[12-15]。

有關(guān)灌溉施肥下地表水流中氮素遷移分布雖已做了大量研究，但大多選擇井水作為灌溉水源且只對其中1～2種氮素的時空變化作了分析，未考慮水質(zhì)、水力等的作用，以及不同氮素之間遷移轉(zhuǎn)化的影響[16-20]。雖然也有部分學(xué)者研究了含氮污染物在高含沙水體中的遷移轉(zhuǎn)化規(guī)律，但主要在室內(nèi)，沒有深入研究自然條件下氮素的遷移轉(zhuǎn)化規(guī)律[21]。因此，開展表施尿素條件下引黃畦灌灌水流中氮素遷移轉(zhuǎn)化規(guī)律的研究具有理論和應(yīng)用價值。

1 材料與方法

1.1 試驗區(qū)概況

田間試驗在山西省尊村引黃灌區(qū)第五提水站北1 500 m玉米地進行。試驗區(qū)位于北緯35°03′、東經(jīng)110°49′,屬暖溫帶大陸性季風(fēng)氣候，多年平均降水量525 mm。地表水資源年均利用量26 177 萬m3，其中提引黃河水16 956 萬m3，約占65%，黃河水是當(dāng)?shù)氐闹饕Y源。液施和撒施化肥是兩種常見的施肥方法，夏玉米(6月初到10月初)和冬小麥(10月初到次年6月初)連作為主的糧食作物種植模式，試驗作物為鄭單958玉米。試驗區(qū)畦田與灌溉施肥基本參數(shù)見表1，液施和撒施尿素的施肥量相同。

1.2 田間試驗設(shè)計

如表1所示，選取施肥方式和地表水流入滲時間作為試驗的兩個因素，設(shè)計4個處理。具體如下：

灌水前，對于撒施，沿畦長方向每隔50 m作為一個施肥單元，化肥分塊均勻撒進畦田；液施則根據(jù)尿素的溶解度、每畝尿素使用量以及灌水時間、泵的流量來計算注入灌溉水中所需水量，其中灌水和取水時間是根據(jù)先前灌水試驗和水流推進過程曲線推算求得的。

表1 田間試驗設(shè)計方案Tab.1 Strategy of field experiment design

注：改水成數(shù)指地面灌溉過程中水流推進的距離占整個田面長的比率；施肥成數(shù)在撒施中指肥料撒在田面中的面積占整個田面面積的比率，而在液施中指注入肥液的時間占灌水時間的比率。

沿程共布設(shè)8個取樣點分別為在距畦首出水口1、50、100、150、200、250、300、345 m處，并提前做好標(biāo)記。

當(dāng)水流推進到各個取樣點時，取開始入滲水樣2 500 mL；并開始計時，等30 min，取各個取樣點水流到達后入滲30 min的水樣2 500 mL。

試驗日期為2016年8月4日，水源為黃河水(總氮3.01 mg/L、氨氮0.4 mg/L、硝氮2.45 mg/L),肥料為尿素(含氮率≥46%)。試驗區(qū)畦田與灌溉施肥基本參數(shù)見表2。

表2 畦田與灌溉施肥基本參數(shù)Tab.2 Basic parameters of irrigation and fertigation

1.3 測試方法

(1)水流測量。使用超聲波明渠流量計(南京寶威儀器儀表有限公司)計量入畦的黃河水流量、秒表計量灌水時間與各個點的水流推進時間。

(2)氮素測量。水樣根據(jù)《水和廢水監(jiān)測分析方法(第四版)》，先使用0.45 μm的濾膜過濾，然后進行酸化使pH≤2、冷藏[22]，以及使用紫外-可見分光光度計HJ-013測量總氮(TN)、流動分析儀設(shè)備(Auto Analyzer Ⅲ,德國Bran+Luebbe)測定氨氮(NH+4-N)與硝氮(NO-3-N)濃度。有機氮和氨氮在水體中都為水相的濃度。

(3)泥沙測量。采用電子天平(上海佑科公司，精度0.01)測定水樣中泥沙質(zhì)量。

1.4 濃度的計算

氮的存在形式包括有機氮(Org-N)和無機氮(Inorg-N)，其在水環(huán)境中的轉(zhuǎn)化作用包括氨化作用、硝化作用和反硝化作用，硝化過程中由亞硝態(tài)氮(NO-2-N)轉(zhuǎn)化為硝氮(NO-3-N)的速度很快，亞硝態(tài)氮很少累積[21]，因此在本文中不考慮亞硝態(tài)氮。各個時刻各個測點氮素的濃度為測試的濃度與黃河水中的濃度之差。則具體濃度計算如下：

CTN=COrg-N+CInorg-N

CInorg-N=CNH+4-N+CNO-3-N

式中：CTN為總氮濃度，mg/L；COrg-N為有機氮濃度，mg/L；CInorg-N為無機氮濃度，mg/L；CNH+4-N為氨氮濃度，mg/L；CNO-3-N為硝氮濃度，mg/L。

1.5 空間分布差異評價方法

采用變異系數(shù)Cv對不同施肥方式和入滲時間下的地表水流中氮素濃度沿畦長的空間分布差異性進行定量描述。其中Cv≤0.1 為弱變異性，0.1

1.6 相關(guān)關(guān)系評價方法

使用相關(guān)系數(shù)r對不同施肥方式和入滲時間下的地表水流中有機氮、氨氮和硝氮濃度之間以及與他們距離、流速、含沙量的相關(guān)關(guān)系的密切程度進行定量描述[24]。

2 結(jié)果與分析

氮素的遷移與轉(zhuǎn)化是一個同時進行的過程，彼此相互聯(lián)系，相互影響，并不是相互獨立的不同方面[18,21,25]。

2.1 氮素的轉(zhuǎn)化

圖1和圖2分別給出了撒施和液施尿素下不同時間地表水流中的總氮、有機氮、氨氮和硝氮的濃度沿畦長的變化的分布圖，從圖中可以看出各個測點氨氮和硝氮的濃度較總氮、有機氮的低，經(jīng)計算無機氮占總氮的比例為14%～30%，有機氮占總氮71%～89%，這說明水中氮素主要以有機氮的形式存在，尿素作為人工合成的有機氮，其水解要經(jīng)歷一個緩慢的過程，同時尿素的水解速率受到有機氮濃度、溫度、pH等影響而不同[26]。并且通過計算氨氮主要在0.2～0.8 mg/L，含量較低，這與有機氮通過氨化作用轉(zhuǎn)化為氨氮的時間較短和氨氮極易轉(zhuǎn)化且易溶于水有關(guān)[21]；硝氮主要在0.5～5 mg/L，較氨氮濃度較高，這可能與適宜的溫度下氨氮迅速轉(zhuǎn)化為硝氮和反硝化速率小于消化速率導(dǎo)致硝氮累積，以及表層土壤中氮素溶解[18]等有關(guān)。

圖1 撒施不同時刻地表水流中氮素濃度沿畦長變化Fig.1 Variation of nitrogen along the border lengthwise in surface water at different time under scattering urea

圖2 液施不同時刻地表水流中氮素濃度沿畦長變化Fig.2 Variation of nitrogen along the border lengthwise in surface water at different time under fertigating urea

表3表明除了液施下有機氮濃度與含沙量沿程正相關(guān)，其余都與含沙量負相關(guān)，且|r|>0.5為中度或高度相關(guān)。由于在撒施條件下各個時間各個測點的氮素的濃度不同，且變化幅度大，難以分析氮素與含沙量之間的關(guān)系；而液施各個時間各個測點的總氮濃度基本相同，變化幅度很小，因此，在此主要分析液施條件下氮素濃度與含沙量的關(guān)系。由圖3(b)可知道含沙量除了局部地區(qū)由于水流流速較大帶起田面表層土壤外都是沿程遞減的[27]，即含沙量越高有機氮含量越高，含沙量對有機氮的轉(zhuǎn)化有抑制作用，這與其他研究者的室內(nèi)研究結(jié)果表明懸浮顆粒物對有機氮的降解有明顯的促進作用的矛盾[13,21],可能是由于水流推進過程中尿素在不斷進行水解，氮素之間在進行不斷轉(zhuǎn)化，以及表層土壤養(yǎng)分隨水流遷移累積[18]等有關(guān),表明引黃灌溉施肥下地表水流中氮素沿畦長的遷移轉(zhuǎn)化不僅僅與含沙量有關(guān)，是眾多因素的作用，如懸浮顆粒物的顆粒組成以及微生物數(shù)量等[21]，氨氮、硝氮與含沙量負相關(guān)，也說明此問題。

圖3 地表水流沿畦長推進速度和含沙量Fig.3 Advancing velocity and sediment content in surface water along the border lengthwise

施肥方式取樣時間有機氮氨氮硝氮撒時開始入滲入滲30min-0．53-0．97-0．78-0．97-0．64-0．91液施開始入滲入滲30min0．690．61-0．58-0．79-0．66-0．66

表4 開始時刻氮素濃度變化與流速變化的相關(guān)系數(shù)Tab.4 Coefficient of correlation between concentrations of nitrogen and velocity in surface water at begin

表5 不同施肥方式下各時刻氮素濃度與距離的相關(guān)系數(shù)Tab.5 Coefficient of correlation between concentrations of nitrogen in surface water and distance

注：*顯著水平0.05，**顯著水平0.01。

白美健等在選用硫酸銨作為試驗肥料時，認(rèn)為表層土壤中氮素的溶解是地表水流中氮素的主要來源，但對比圖1和圖2的可以發(fā)現(xiàn)，在相同時刻液施和撒施在總氮不同情況下，氨氮和硝氮的濃度也不同，一般有機氮濃度越高，氨氮和硝氮的濃度也越高，這表明水流中的氨氮、硝氮不僅來來源于土壤，尿素含量也起著重要作用。這可能是由于硫酸銨分解溫度較高，而尿素在適宜的溫度下完全分解較快[26]，且本次試驗的灌溉水源為夏季黃河水，水溫和懸浮顆粒物較高。

2.2 氮素與流速、距離和時間的關(guān)系

在畦灌地面水流中，溶質(zhì)的遷移受到時間平均流速、紊流速度變化以及分子擴散等影響[28]。在畦灌過程中水流在田面上流動與土壤下滲是同時進行的，從水力學(xué)上看是透水地板上的明渠非恒定流[1]。從表4可以得到在撒施下沿畦長方向氮素濃度變化與流速變化均是正相關(guān)，其中總氮和有機氮的相關(guān)系數(shù)值大于0.5，是中度相關(guān)，這是因為尿素水解轉(zhuǎn)化前不帶電荷，不易被土粒吸附，很容易隨水移動和流失[26]，且有機氮在總氮中所占比率大的原因。然而硝氮為低度相關(guān)，氨氮為極低度相關(guān)，這是可能是由于氨氮帶正電荷，容易被帶負電荷的土壤和懸浮物顆粒物吸附，然而，硝氮帶負電荷，不與帶負電荷的泥沙顆粒結(jié)合，極易隨水流流失[21]；同時從圖1可以看出氮素濃度在畦首的變化比畦尾大，這可能是由于水流在推進過程中開始施推進速度迅速減小，然后逐漸減小。但對于液施，各個氮素的相關(guān)系數(shù)小于撒施的，且總氮和有機氮與流速為正相關(guān)，但相關(guān)系數(shù)較小，這表明在液施情況下總氮和有機氮濃度變化受流速變化作用較小，而氨氮、硝氮濃度變化與流速變化是負相關(guān)，這說明流速變化在其中起的不是主要原因。

從表5可以看出，不同施肥方式下不同時刻氮素濃度除了液施有機氮外其他氮素與畦首距離的相關(guān)系數(shù)都為正數(shù)，且撒施的相關(guān)系數(shù)普遍高于液施的，撒施的相關(guān)系數(shù)都大于0.8，與距離為高度相關(guān)，這表明撒施情況下，氮素在畦首存在明顯的流失，在畦尾存在不斷累積，這可能與畦首流速較大，以及尿素隨水流失，氨氮主要是通過有機氮轉(zhuǎn)化得到，硝氮通過氨氮轉(zhuǎn)化得到，有關(guān)尿素、氨氮和硝氮之間構(gòu)成是一個鏈數(shù)為三的連續(xù)衰減反應(yīng)[29]；并且30 min時的相關(guān)系數(shù)大于開始時的，這可能是由于一方面在撒施開始時尿素為顆粒狀，相對難以沖刷，盡管尿素易溶解于水且速率快， 30 min時全部溶解，很容易流失；另一方面在開始時由于表層土壤含水率較低，土壤非飽和導(dǎo)水率較大，30 min時土壤非飽和導(dǎo)水率較小，水流的流速大于開始時的。對于液施，總氮與距離的相關(guān)系數(shù)小于0.3，為極低度相關(guān)，可能存在畦尾累積，但不明顯，這是由于液施時水泵流量和灌水流量穩(wěn)定，水流中尿素分布均勻，且尿素的水解要經(jīng)歷一個的過程，沿畦長氨氮、硝氮的所占比率比較穩(wěn)定且較小所致；而有機氮與距離為負相關(guān)，氨氮和硝氮為正相關(guān)，相關(guān)程度都為中度以上，這主要是由于在田面水流推進過程中尿素在進行不斷水解，導(dǎo)致有機氮減少，氨氮和硝氮增加。

從圖1可以看出，撒施各測點地表水流中的總氮的濃度隨時間在減小，且畦首氮素濃度降低的幅度大于畦尾。

2.3 氮素的時空分布特征

不同施肥方式下各個時刻氮素遷移轉(zhuǎn)化不同直接導(dǎo)致了氮素分布的不同。圖1、圖2與表5結(jié)果顯示，撒施的氮素時空變異性大于液施，這表明液施較撒施可以提高地表水流中氮素的時空分布的均勻性，彌補撒施在均勻性方面的不足，但對于不同施肥方法不同時間下不同氮素沿畦長分布的變異強度不同，其中同種施肥方式下的總氮、有機氮的變異系數(shù)值相差較小。至于撒施，3個氮素濃度沿畦長在不同時間都是中等變異性，且總氮、有機氮的變異系數(shù)都小于氨氮和硝氮，至于在開始時刻總氮、有機氮和硝氮的變異系數(shù)都小于氨氮，而在30 min時四氮的變異系數(shù)都在0.59～0.69之間，差異不大，但總氮、硝氮變異系數(shù)相對于開始時刻增加了30%多，氨氮有所減小。而對于液施，3個氮素濃度沿畦長在不同時間變異等級不同，總氮與有機氮都為弱變異，氨氮和硝氮都為中等變異，但氨氮遠遠高于硝氮的，在不同時刻總氮、有機氮和硝氮變異系數(shù)變化較小，而氨氮的變化較大且同樣是開始時刻大于30 min的。這表明不同氮素時空分布具有不同步性，各自的影響因素不同。

表6 不同施肥方式下各時刻氮素濃度沿畦長的變異系數(shù)Tab.6 Coefficient of variation of nitrogen in surface water along the border lengthwise

注：a、b表明變量之間的顯著性，如果數(shù)字后面的字母一致，則差異不顯著，否則差異顯著。

3 討 論

對于田間灌溉施肥，相比于有機氮、氨氮和硝氮而言，總氮分布特征直接影響著氮素在土壤中分布的均勻性，作物對氮素的利用率，以及氮素的流失等。

本次試驗由于入畦單寬流量大，而畦埂高度不夠，限制了水流深度、導(dǎo)致水流在田間運動過程中存在少許漏水，對水流的推進產(chǎn)生一定的影響；在畦灌過程中，水流在田面運動其實時間較短，而由于黃河水含沙率較高，在水樣監(jiān)測前要經(jīng)歷靜沉等處理，延長了氮素轉(zhuǎn)化的時間，這可能導(dǎo)致結(jié)果中各個測點中氨氮、硝氮的濃度偏高；且本文中氮素濃度都為水相中的，不涉及顆粒相的，有待進一步研究。

在地面灌溉施肥中，水流的運動是肥料遷移轉(zhuǎn)化的水力基礎(chǔ)[28]，國內(nèi)外學(xué)者主要基于守恒型Saint-Venant方程描述地表水流運動，Richards方程描述田面水流的入滲和用對流彌散方程來描述溶質(zhì)的遷移過程[8,30]。聯(lián)合這些模型能較好地模擬溶質(zhì)濃度沿程變化，彌補地面灌溉施肥上設(shè)計和管理的不足。

4 結(jié) 語

通過分析不同施肥方式下地表水流中的氮素遷移轉(zhuǎn)化和時空分布特征，可以得到如下結(jié)論。

(1)水流中有機氮占總氮71%～89%，總氮主要以有機氮的形式存在，尿素在的水解需要經(jīng)歷一個過程；地表水流中的無機氮不僅來源于表層土壤，尿素的水解也占有重要的部分；懸浮顆粒物不是引起沿程不同氮素變化的主要原因。

(2)撒施下氮素遷移受流速作用較大，氮素在畦首流失，畦尾累積，總氮、有機氮、氨氮和硝氮與距離的相關(guān)系數(shù)都大于0.8，為高度相關(guān)。液施受流速、距離和時間的影響較小。不同氮素在水流中特性不同導(dǎo)致氮素遷移特征不同。

(3)施肥方式對地表水流中氮素的分布特征有顯著影響，撒施的氮素時空變異性大于液施的，不同氮素的時空變異特征不同，液施較撒施可以提高地表水流中氮素的時空分布的均勻性，且撒施30 min時的變異系數(shù)比開始大0.38。

□

[1] 張新燕, 蔡煥杰. 地面灌溉理論研究進展[J]. 中國農(nóng)村水利水電, 2005,(8):6-9.

[2] 闞常慶, 廖梓龍, 龍胤慧. 畦灌技術(shù)研究進展[J]. 水科學(xué)與工程技術(shù), 2012,(3):1-4.

[3] Dillon J, Edinger-Marshall S, Letey J. Farmers adopt new irrigation and fertilizer techniques: Changes could help growers maintain yields, protect water quality[J]. California Agriculture, 1999,53(1):24-31.

[4] 趙銀亮, 宋華力, 毛艷艷. 黃河流域糧食安全及水資源保障對策研究[J]. 人民黃河, 2011,33(11):47-49.

[5] 史學(xué)斌, 馬孝義, 聶衛(wèi)波,等. 地面灌溉的研究現(xiàn)狀與發(fā)展趨勢[J]. 水資源與水工程學(xué)報, 2005,16(1):34-40.

[6] 侯紅雨, 王洪梅, 肖素君. 黃河流域灌溉發(fā)展規(guī)劃分析[J]. 人民黃河, 2013,35(10):96-98.

[7] 張瑜芳，張蔚榛，沈榮開等. 排水農(nóng)田中氮素轉(zhuǎn)化運移和流失[M]. 武漢：中國地質(zhì)大學(xué)出版社，1997:20-28.

[8] Zerihun D, Furman A, Warrick A W, et al. Coupled surface-subsurface solute transport model for irrigation borders and basins. I. Model development[J]. Journal of Irrigation & Drainage Engineering, 2005,131(5):396-406.

[9] 李志新, 許 迪, 李益農(nóng),等. 畦灌施肥地表水流與非飽和土壤水流-溶質(zhì)運移集成模擬Ⅲ:模型應(yīng)用[J]. 水利學(xué)報, 2011, 42(3):271-277.

[10] 張立志, 徐立榮, 徐征和,等. 引黃灌區(qū)畦田灌水技術(shù)試驗研究[J]. 節(jié)水灌溉, 2016,(8):109-112.

[11] 姚 欣,李金山,黃修橋,等.引黃畦灌田間水沙分布規(guī)律[J].農(nóng)業(yè)工程報,2016,32(18):147-152.

[12] Xia X H, Yang Z F, Huang G H, et al. Nitrification in natural waters with high suspended-solid content——A study for the Yellow River[J]. Chemosphere, 2004,57(8):1 017-29.

[13] 張學(xué)青, 夏星輝, 楊志峰. 水體顆粒物對有機氮轉(zhuǎn)化的影響[J]. 環(huán)境科學(xué), 2007, 28(9):1 954-1 959.

[14] Sigleo A C, Macko S A. Carbon and nitrogen isotopes in suspended particles and colloids, Chesapeake and San Francisco Estuaries, U.S.A.[J]. Estuarine Coastal & Shelf Science, 2002,54(4):701-711.

[15] Turner A, Millward G E. Suspended Particles: Their role in estuarine biogeochemical cycles[J]. Estuarine Coastal & Shelf Science, 2002,55(6):857-883.

[16] Abbasi F, Feyen J, Roth R L, et al. Water flow and solute transport in furrow-irrigated fields[J]. Irrigation Science, 2003,22(2):57-65.

[17] 劉 懿, 繳錫云, 張 仙,等. 表施肥料條件下畦灌地表水流中氮素質(zhì)量濃度變化規(guī)律[J]. 灌溉排水學(xué)報, 2013,32(5):34-37.

[18] 白美健, 許 迪, 李益農(nóng),等. 畦灌撒施與液施硫酸銨地表水流和土壤中氮素時空分布特征[J]. 農(nóng)業(yè)工程學(xué)報,2011,27(8):19-24.

[19] 于 飛. 地面施肥灌溉條件下水氮分布試驗研究[D]. 北京：中國水利水電科學(xué)研究院， 2007.

[20] 梁艷萍. 畦灌施肥土壤水氮時空分布試驗研究[D]. 北京：中國水利水電科學(xué)研究院， 2008.

[21] 張學(xué)青. 黃河水體懸浮顆粒物對含氮污染物遷移轉(zhuǎn)化作用的影響[D]. 北京：北京師范大學(xué), 2006.

[22] 魏復(fù)盛. 水和廢水監(jiān)測分析方法[M]. 4版. 北京：環(huán)境科學(xué)出版社, 2002.

[23] 南京農(nóng)科院.田間實驗和統(tǒng)計方法[M]. 北京：農(nóng)業(yè)出版社，1979.

[24] 劉愛利. 地統(tǒng)計學(xué)概論[M]. 北京：科學(xué)出版社, 2012.

[25] 王 康. 非飽和土壤水流運動及溶質(zhì)遷移[M]. 北京：科學(xué)出版社, 2010.

[26] 沈其榮. 土壤肥料學(xué)通論[M]. 北京：高等教育出版社, 2001.

[27] 費良軍, 王文焰. 渾水波涌畦灌田面沉積泥沙沿畦長分布規(guī)律試驗研究[J]. 西安理工大學(xué)學(xué)報, 1999, 15(1):37-41.

[28] Cunge J A, Holly F M, Verwey A. Practical aspects of computational river hydraulics[M]. Pitman Advanced Pub. Program, 1980.

[30] 章少輝, 許 迪, 李益農(nóng),等. 撒施肥料一維畦灌地表水流與溶質(zhì)運移模型構(gòu)建與驗證[J]. 農(nóng)業(yè)工程學(xué)報, 2010,26(12):34-39.