基于彈道理論有風(fēng)條件下折射式噴頭噴灌均勻度研究

張以升 朱德蘭, 宋 博 張 林

(1.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院， 陜西楊凌 712100；2.西北農(nóng)林科技大學(xué)中國(guó)旱區(qū)節(jié)水農(nóng)業(yè)研究院， 陜西楊凌 712100)

基于彈道理論有風(fēng)條件下折射式噴頭噴灌均勻度研究

張以升1朱德蘭1,2宋 博1張 林2

(1.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院， 陜西楊凌 712100；2.西北農(nóng)林科技大學(xué)中國(guó)旱區(qū)節(jié)水農(nóng)業(yè)研究院， 陜西楊凌 712100)

為計(jì)算有風(fēng)條件下折射式噴頭水量分布及噴灌均勻度，以彈道軌跡理論為基礎(chǔ)，依據(jù)風(fēng)速分布模型，建立有風(fēng)條件下折射式單噴頭水量分布計(jì)算方法，采用該方法模擬出有風(fēng)條件下Nelson D3000型噴頭倒掛安裝方式下水量分布特性，通過(guò)與實(shí)測(cè)資料進(jìn)行對(duì)比，驗(yàn)證了模擬具有較高的準(zhǔn)確度，可應(yīng)用于有風(fēng)條件下折射式噴頭水量分布計(jì)算。在此基礎(chǔ)上，選用4.76 mm(24號(hào))噴嘴直徑，模擬出不工況下單噴頭水量分布，計(jì)算出組合情況下噴灌均勻度，分析了風(fēng)速、風(fēng)向、噴頭間距、工作壓力和安裝高度5種因素對(duì)噴灌均勻度的影響，并對(duì)蒸發(fā)漂移損失進(jìn)行了分析。結(jié)果表明：95%的置信區(qū)間下，噴頭布置間距對(duì)噴灌均勻度的影響最顯著，其次是安裝高度和噴頭工作壓力，風(fēng)速和風(fēng)向?qū)姽嗑鶆蚨扔绊懖伙@著。風(fēng)速、噴頭工作壓力和安裝高度都會(huì)對(duì)蒸發(fā)漂移損失產(chǎn)生影響，其中工作壓力影響最大。當(dāng)選用Nelson D3000型噴頭在風(fēng)速小于6 m/s的環(huán)境下噴灌時(shí)，應(yīng)將噴頭安裝間距固定在2.13～3.04 m范圍內(nèi)。另外，該安裝間距范圍內(nèi)，噴頭安裝高度和噴灌壓力增大后，噴灌均勻度增大的效果不明顯，因此應(yīng)采用低壓噴灌以降低噴灌系統(tǒng)運(yùn)行成本；考慮到較高的噴頭安裝高度會(huì)產(chǎn)生較大的蒸發(fā)漂移損失，噴灌時(shí)還應(yīng)適當(dāng)降低噴頭安裝高度，以提高噴灌水分利用率。

折射式噴頭； 噴灌均勻度； 水量分布； 有風(fēng)條件； 彈道理論

引言

噴灌水量分布是計(jì)算噴灌均勻系數(shù)及噴灑水利用系數(shù)的基礎(chǔ)[1]。為了設(shè)計(jì)方便，基于彈道理論的單噴頭水量分布模擬方法被廣泛應(yīng)用[2～4]，筆者此前依據(jù)彈道軌跡理論，針對(duì)折射式噴頭特有的噴灑形式，建立了折射式噴頭水量分布計(jì)算模型，可在已知單流道水量分布的前提下，計(jì)算出無(wú)風(fēng)條件下折射式噴頭水量分布特性[5]。然而農(nóng)田噴灌作業(yè)一般都在室外進(jìn)行，大氣溫濕度和風(fēng)速等環(huán)境因素也會(huì)對(duì)水量分布和噴灑水利用系數(shù)產(chǎn)生影響[6-9]，F(xiàn)ACI等[10]通過(guò)室外試驗(yàn)發(fā)現(xiàn)風(fēng)速會(huì)影響Nelson D3000型噴頭平均噴灌強(qiáng)度，并使噴灌區(qū)域產(chǎn)生一定程度的偏移，對(duì)水量分布形狀沒(méi)有明顯影響。TARJUELO等[11]指出，在干燥多風(fēng)的環(huán)境下實(shí)施噴灌作業(yè)，易產(chǎn)生較大的蒸發(fā)漂移損失，因此采用無(wú)風(fēng)條件下水量分布特性進(jìn)行噴灌設(shè)計(jì)會(huì)產(chǎn)生一定的誤差。

為了進(jìn)一步拓寬軟件應(yīng)用范圍，減小室外試驗(yàn)的工作量，本文以彈道軌跡理論為基礎(chǔ)，依據(jù)風(fēng)速分布模型，建立有風(fēng)條件下折射式單噴頭水量分布計(jì)算方法，通過(guò)與實(shí)測(cè)資料對(duì)比，驗(yàn)證模擬的精準(zhǔn)性。在此基礎(chǔ)上，分析風(fēng)速、風(fēng)向、噴頭間距、工作壓力和安裝高度等因素對(duì)4.76 mm噴嘴直徑Nelson D3000型噴頭噴灌均勻度的影響，并對(duì)有風(fēng)條件下噴頭蒸發(fā)漂移損失進(jìn)行計(jì)算，以期為有風(fēng)條件下噴灌提供設(shè)計(jì)參數(shù)。

1 計(jì)算模型及數(shù)據(jù)處理

文中涉及的折射式噴頭出射速度、角度模型、水滴三維運(yùn)動(dòng)模型和水滴運(yùn)動(dòng)蒸發(fā)模型，參考文獻(xiàn)[5]中所述。

1.1 有風(fēng)條件下水量分布計(jì)算方法

VORIES等[12-13]指出，在連續(xù)的空氣環(huán)境中，風(fēng)速隨高度的變化為

(1)

其中

lgd=0.979 3lgz-0.153 6
lgZ0=0.997lgz-0.883

式中W(z)——水滴下降至高度z處的計(jì)算風(fēng)速，m/s

Wa——距地面高度a(本文中a為300 cm)處的設(shè)計(jì)風(fēng)速，m/s

d——粗糙高度，cm

Z0——粗糙系數(shù)，cm

有風(fēng)條件下水量分布計(jì)算如圖1所示[14]，其中n為坐標(biāo)間隔，A、B、C、D、Q為不同粒徑水滴的落點(diǎn)。則水量分布計(jì)算方法如下：

圖1 有風(fēng)條件下單噴頭粒徑(水量)分布Fig.1 Droplet or water distribution under windy condition

(1)確定噴灑區(qū)域范圍。大粒徑水滴受風(fēng)拖拽影響較小，模擬最大粒徑水滴在不同方向上的落地點(diǎn)，以此確定噴灑區(qū)域(折射式噴頭最大水滴直徑分布在最末端)。確定噴灑區(qū)域后，對(duì)噴灑區(qū)域進(jìn)行網(wǎng)格化處理，由于網(wǎng)格化尺寸很大程度上影響到計(jì)算的精度和計(jì)算效率，李永沖等[1]設(shè)虛擬小網(wǎng)格尺寸為1 m，基本能夠滿足計(jì)算精度要求，為了提高軟件的應(yīng)用效果，在程序設(shè)計(jì)過(guò)程中，網(wǎng)格尺寸可根據(jù)所模擬噴頭的噴灑半徑進(jìn)行手動(dòng)輸入，由于所模擬的噴頭噴灑半徑一般不大于10 m，根據(jù)模擬經(jīng)驗(yàn)，0.5 m的網(wǎng)格尺寸基本能夠達(dá)到較高的模擬精度，因此本文網(wǎng)格化尺寸定為0.5 m，以便計(jì)算有風(fēng)條件下不同粒徑水滴落地點(diǎn)坐標(biāo)。

(2)不考慮有風(fēng)條件下水滴碎裂及水滴運(yùn)動(dòng)軌跡相互干擾，采用彈道理論，計(jì)算折射式噴頭有風(fēng)條件下在0～360°范圍內(nèi)水滴粒徑在網(wǎng)格化區(qū)域內(nèi)的分布。模擬過(guò)程中，由于水滴在運(yùn)動(dòng)過(guò)程中產(chǎn)生蒸發(fā)損失，計(jì)算蒸發(fā)時(shí)，用單個(gè)水滴體積減小百分?jǐn)?shù)與該粒徑對(duì)應(yīng)的單位時(shí)間內(nèi)噴灌水量(即噴灌強(qiáng)度)的乘積作為該點(diǎn)單位時(shí)間內(nèi)蒸發(fā)的水量，則噴灌時(shí)的蒸發(fā)損失為所模擬水滴對(duì)應(yīng)的單位時(shí)間內(nèi)總的蒸發(fā)水量與單位時(shí)間內(nèi)測(cè)得的總落地水量的比值。由于水滴在風(fēng)的作用下運(yùn)動(dòng)軌跡發(fā)生了偏移，且水滴直徑越小偏移原落點(diǎn)的距離越遠(yuǎn)，這樣就致使一部分小粒徑水滴不能落入網(wǎng)格化區(qū)域，該部分水滴代表的水量作為飄移損失不進(jìn)行水量分布的計(jì)算。

(3)根據(jù)水滴粒徑分布規(guī)律，結(jié)合不同粒徑水滴對(duì)應(yīng)的水量，對(duì)有風(fēng)條件下單噴頭水量分布進(jìn)行計(jì)算，水量分布計(jì)算原則如圖1所示。假設(shè)某流道射流出的某一粒徑水滴落在Q點(diǎn)，其他流道射流出的不同粒徑水滴分別落在A、B、C、D點(diǎn)，通過(guò)分別對(duì)比Q點(diǎn)與A、B、C、D點(diǎn)的距離，判斷出Q點(diǎn)與A點(diǎn)距離最近，則將Q點(diǎn)處水滴粒徑對(duì)應(yīng)的水量疊加至A點(diǎn)處。

1.2 蒸發(fā)飄移損失計(jì)算方法

蒸發(fā)飄移損失計(jì)算包括兩方面內(nèi)容，首先是水滴運(yùn)動(dòng)過(guò)程中的蒸發(fā)損失率，計(jì)算式為

(2)

圖2 水量分布模擬與測(cè)量值對(duì)比Fig.2 Comparison of measured and simulated water distributions

式中E——蒸發(fā)損失率，%n——模擬的不同粒徑水滴序號(hào)Dci——第i個(gè)水滴D的初始直徑，mDli——第i個(gè)水滴D的落地直徑，mPDi——第i個(gè)水滴D對(duì)應(yīng)的噴灌強(qiáng)度，mm/h

同理，飄移損失率的計(jì)算式為

(3)

式中m——漂移出計(jì)算區(qū)域的不同粒徑水滴序號(hào)PDt——第t個(gè)漂移出計(jì)算區(qū)域水滴對(duì)應(yīng)的噴灌強(qiáng)度，mm/h

由式(2)和式(3)可得蒸發(fā)飄移損失為

L=E+W

(4)

2 有風(fēng)條件下計(jì)算模型驗(yàn)證

單噴頭噴灑面積大，且水量分布試驗(yàn)需要長(zhǎng)時(shí)間測(cè)量，而自然條件下風(fēng)速和風(fēng)向穩(wěn)定性較差，測(cè)量準(zhǔn)確性不能得到保證；創(chuàng)造大區(qū)域風(fēng)洞試驗(yàn)條件，其花費(fèi)又十分昂貴。各種因素限制使得準(zhǔn)確獲取有風(fēng)條件下單噴頭水量分布數(shù)據(jù)較為困難，為驗(yàn)證有風(fēng)條件下單噴頭水量分布模擬的準(zhǔn)確性，根據(jù)FACI的試驗(yàn)條件設(shè)定相應(yīng)的模擬參數(shù)，其中噴嘴直徑為3.8 mm，噴頭工作壓力為140 kPa，模擬風(fēng)速為2.5 m/s，模擬了1 m和2.5 m兩種高度下的水量分布，通過(guò)與其試驗(yàn)值對(duì)比，對(duì)有風(fēng)條件下的單噴頭水量分布和噴灑射程的模擬值與測(cè)量值的差異性進(jìn)行分析。

2.1 單噴頭水量分布

圖2為工作壓力為140 kPa、噴嘴型號(hào)為3.8 mm、風(fēng)速為2.5 m/s水量分布模擬結(jié)果與FACI測(cè)量結(jié)果的對(duì)比，坐標(biāo)點(diǎn)(0，0)為噴頭安裝位置，風(fēng)向?yàn)槠叫杏趛軸方向。水量分布圖顯示，有風(fēng)條件下模擬與實(shí)測(cè)的水量分布與無(wú)風(fēng)條件下類似，均呈現(xiàn)圓形分布，與測(cè)量值相比，模擬值邊界較為規(guī)則，另外，模擬水量分布在周向上呈條帶形分布，這是由于該噴盤(pán)具有36個(gè)流道，噴頭在固定風(fēng)速下工作時(shí)，理想狀況下各射流水束相互無(wú)干擾；實(shí)測(cè)水量分布在周向上水量分布相對(duì)較為均勻，漏噴區(qū)域不明顯，這是因?yàn)闇y(cè)量時(shí)風(fēng)速與風(fēng)向都無(wú)法固定(測(cè)量風(fēng)速變化范圍為0～2.5 m/s)，當(dāng)風(fēng)速和風(fēng)向發(fā)生變化時(shí)，其射流水束落地位置發(fā)生了變化，水束偏移后一定程度上補(bǔ)充了漏噴區(qū)的水量。從圖中還可以看出，水量在迎風(fēng)射程上(y軸正方向)水舌條帶較短，安裝高度越高，效果越明顯，這是由于水束在迎風(fēng)向噴射時(shí)，大粒徑水滴影響較小，迎風(fēng)繼續(xù)飛行，小粒徑水滴受風(fēng)的作用飛向順風(fēng)向區(qū)域，導(dǎo)致迎風(fēng)向水量減少，落地條帶的長(zhǎng)度變短。從圖中還可以看出，噴頭安裝高度增大至2.5 m后，平均噴灌強(qiáng)度減小。通過(guò)對(duì)比分析可知，模擬出的有風(fēng)條件下水量分布規(guī)律基本能夠反映真實(shí)情況。

2.2 單噴頭噴灑射程

有風(fēng)條件下，噴頭安裝高度不同，水量偏離噴頭中心點(diǎn)的距離也有差異，表1給出了有風(fēng)條件下噴頭安裝高度為1 m和2.5 m時(shí)噴灑射程的模擬值與測(cè)量值的對(duì)比結(jié)果，數(shù)據(jù)顯示，安裝高度為1 m時(shí)，噴頭在x軸方向、順風(fēng)向、迎風(fēng)向的射程模擬值分別為5.72、6.48、4.95 m，測(cè)量值分別為5.87、6.36、5.52 m，對(duì)應(yīng)的相對(duì)偏差分別為2.60%、1.90%、-10.30%，安裝高度為2.5 m時(shí)對(duì)應(yīng)的相對(duì)偏差為0.10%、4.40%、0.20%，總體而言，射程模擬與實(shí)測(cè)值的偏差相對(duì)較小，模擬準(zhǔn)確性有一定的保證。

表1 噴灑射程模擬值與FACI計(jì)算值的對(duì)比Tab.1 Comparison of simulated and FACI calculated sprinkler range

3 噴灌均勻度

3.1 噴灌均勻度影響因素分析

以Nelson D3000噴頭為研究對(duì)象，采用該計(jì)算程序，模擬出噴嘴直徑4.76 mm、倒掛安裝方式、不同工作壓力和安裝高度下，風(fēng)速為2、4、6 m/s時(shí)單噴頭水量分布，模擬的因素水平如表2所示，根據(jù)單噴頭水量分布數(shù)據(jù)，采用克里斯琴森均勻度計(jì)算公式[15]，進(jìn)行多噴頭組合疊加計(jì)算組合噴灌均勻度，均勻度計(jì)算式為

(5)

不考慮各因素間的交互作用，采用SPSS軟件分析風(fēng)速Ws、風(fēng)向Wd、工作壓力p、噴頭間距l(xiāng)以及噴頭安裝高度H等因素對(duì)平移式噴灌均勻度的影響，結(jié)果如表3所示。分析結(jié)果顯示，在95%的置信度下，工作壓力、噴頭間距和噴頭安裝高度P值均小于0.05，說(shuō)明這3個(gè)因素對(duì)組合均勻度的影響呈顯著水平，風(fēng)速及風(fēng)向P值分別為0.074和0.098，說(shuō)明風(fēng)速和風(fēng)向?qū)鶆蚨鹊挠绊懖伙@著，與文獻(xiàn)[10]分析結(jié)果一致。另外，從影響噴灌均勻的3個(gè)因素對(duì)應(yīng)的F值來(lái)看，噴頭間距F值最大，說(shuō)明相比噴頭安裝高度和工作壓力，噴頭間距對(duì)噴灌均勻度的影響效果最明顯。

表2 均勻度計(jì)算因素與水平Tab.2 Experimental factors and levels for calculation of CU

表3 各因素對(duì)均勻度影響的方差分析Tab.3 Analysis of variance between each factor and CU

3.2 不同風(fēng)速和風(fēng)向下噴頭間距對(duì)均勻度影響

圖3為噴頭安裝高度為2.0 m，工作壓力分別為100 kPa和200 kPa，不同風(fēng)速和風(fēng)向下噴頭間距對(duì)噴灌均勻度的影響。從圖3中可以看出，不同風(fēng)速和風(fēng)向下，隨著噴頭布置間距的增大，噴灌均勻度均呈現(xiàn)降低的趨勢(shì)。2～5 m噴頭間距下均勻度下降幅度較小，當(dāng)噴頭安裝間距增大至6 m時(shí)，噴灌均勻度降低較為迅速。當(dāng)噴頭壓力為100 kPa時(shí)，噴頭間距為2～5 m時(shí)，2、4、6 m/s風(fēng)速下，均勻度最大降幅分別為7.29%、8.21%和5.30%，增大至6 m后均勻度分別降低了25.89%、17.08%和26.26%(圖3a)，降低幅度較大；噴頭間距由2 m增大至5 m后，30°和60°風(fēng)向下，噴灌均勻度最大降幅分別為8.72%和6.38%，增大至6 m后均勻度分別降低了30.73%和33.23%(圖3b)。當(dāng)壓力增大至200 kPa，噴頭間距在2～5 m范圍內(nèi)，任意風(fēng)速下和風(fēng)向下噴灌均勻度在86.26%以上，噴頭間距增大至6 m后，噴灌均勻度迅速下降至78%以下(圖3c、3d)。從圖中還可以看出，風(fēng)速和風(fēng)向?qū)姽嗑鶆蚨扔绊懢鶡o(wú)顯著性規(guī)律，這與SPSS方差分析得出的結(jié)論一致。

圖3 不同風(fēng)速、風(fēng)向下噴頭間距對(duì)噴灌均勻度的影響Fig.3 Effect of sprinkler spacing on CU under different wind speeds and directions

采用最小二乘法原理對(duì)不同噴頭間距下的均勻度進(jìn)行回歸分析，結(jié)果顯示采用二次多項(xiàng)式回歸系數(shù)較高，關(guān)系式為

C=aL2+bL+c

(6)

式中L——噴頭間距，ma、b、c——待定系數(shù)

采用回歸出的經(jīng)驗(yàn)?zāi)Ｐ?，以均勻度最高為目?biāo)，對(duì)模型求導(dǎo)計(jì)算出最佳噴頭間距，不同安裝高度及工作壓力下的回歸參數(shù)、最佳噴頭間距及最佳噴頭間距下噴灌均勻度如表4所示。從表中可以得出，噴頭安裝高度較低時(shí)，應(yīng)盡量縮小噴頭間距，考慮風(fēng)速和風(fēng)向時(shí)，不同噴頭安裝高度及工作壓力下最佳噴頭安裝間距變化范圍為2.13～3.04 m。

3.3 不同風(fēng)速和風(fēng)向下安裝高度對(duì)噴灌均勻度影響

由以上分析得出噴頭最佳安裝間距為2.13～3.04 m，因此重點(diǎn)分析該間距下噴頭安裝高度對(duì)噴灌均勻度的影響，圖4為噴頭工作壓力為100 kPa，噴頭布置間距為2 m和3 m時(shí)，不同風(fēng)速和風(fēng)向下噴頭安裝高度對(duì)噴灌均勻度影響柱狀圖?？傮w來(lái)說(shuō)，噴頭間距一定時(shí)，噴灌均勻度隨著噴頭安裝高度的增大有增高的趨勢(shì)，然而這種趨勢(shì)并不特別明顯，由此可知，安裝高度的提升對(duì)噴灌均勻度的提高效果不大，因此在有風(fēng)條件下，安裝間距為2～3 m時(shí)，不必為提高噴灌均勻度而刻意增大噴頭安裝高度，增大安裝高度一定程度上還會(huì)造成較大的蒸發(fā)飄移損失，因此在實(shí)際噴灌應(yīng)用中應(yīng)結(jié)合噴灌作物植株高度，盡可能降低噴頭安裝高度。

表4 考慮風(fēng)速風(fēng)向條件下噴灌均勻度隨噴頭間距變化的回歸模型及最佳噴頭間距Tab.4 Optimum spacing calculated by predictive equations for estimation of CU under windy condition

圖4 安裝高度與噴灌均勻度的關(guān)系Fig.4 Relationships between CU and riser height

3.4 不同風(fēng)速和風(fēng)向下工作壓力對(duì)噴灌均勻度影響

同樣選取2 m和3 m噴頭間距，噴頭安裝高度為1 m，分析不同風(fēng)速和風(fēng)向下工作壓力對(duì)噴灌均勻度的影響，結(jié)果如圖5所示。圖中顯示，噴頭間距為2 m和3 m時(shí)，增大噴灌壓力一定程度上可以提高噴灌均勻度，但是提高的尺度有限，由于采用低壓灌溉，很大程度上能夠降低系統(tǒng)運(yùn)行費(fèi)用，因此當(dāng)噴頭安裝間距在2～3 m范圍內(nèi)，盡可能選取較低的工作壓力進(jìn)行噴灌，以降低噴灌系統(tǒng)運(yùn)行成本。

圖5 工作壓力與噴灌均勻度的關(guān)系Fig.5 Relationships between CU and working pressure

4 飄移損失分析

蒸發(fā)漂移損失影響因素較多，其中噴灑水滴直徑、大氣溫濕度和風(fēng)速等環(huán)境因素都會(huì)對(duì)其產(chǎn)生影響。本文在分析蒸發(fā)漂移損失時(shí)，以單噴頭為研究對(duì)象，模擬風(fēng)速、工作壓力以及安裝高度對(duì)蒸發(fā)漂移損失的影響，模擬結(jié)果如圖6所示。由于風(fēng)速增大后小粒徑水滴受風(fēng)的拖拽效果增大，不同工況下蒸發(fā)飄移損失均呈現(xiàn)增高趨勢(shì)；當(dāng)工作壓力為200 kPa、風(fēng)速由4 m/s增大至6 m/s時(shí)，不同高度下飄移損失增加不明顯，這是由于該工作壓力下，4 m/s風(fēng)速已經(jīng)將小粒徑水滴充分吹離目標(biāo)區(qū)域，當(dāng)風(fēng)速增大后對(duì)漂移損失影響較小。當(dāng)風(fēng)速相同時(shí)，飄移損失隨著工作壓力的增大呈現(xiàn)增高的趨勢(shì)，這是因?yàn)閴毫υ龃蠛螅卧诳罩械乃榱演^為明顯，導(dǎo)致在近噴頭處水量增多，在水量分布模擬時(shí)，近噴頭處水量對(duì)應(yīng)的水滴粒徑較小，小粒徑水滴受風(fēng)的作用效果較為明顯，使得小粒徑水滴對(duì)應(yīng)的水量漂移出目標(biāo)區(qū)域，結(jié)果使得蒸發(fā)漂移損失量增加。另外隨著噴頭安裝高度的升高，相同風(fēng)速下蒸發(fā)飄移損失均增大，產(chǎn)生這種現(xiàn)象的主要原因是高度的提升使得噴灑水滴在空氣中的飛行時(shí)間增長(zhǎng)，受風(fēng)的作用持續(xù)增加后，較大粒徑水滴漂移出目標(biāo)區(qū)域，飄移損失量增加。

圖6 不同工況下蒸發(fā)飄移損失率變化曲線Fig.6 Wind drift and evaporation losses under different working conditions

表5給出了風(fēng)速、噴頭工作壓力及安裝高度對(duì)蒸發(fā)飄移損失率影響的分析結(jié)果，數(shù)據(jù)顯示，以上3個(gè)因素對(duì)蒸發(fā)飄移損失率的影響效果均呈顯著水平，從F值大小表明，工作壓力對(duì)蒸發(fā)飄移損失的影響最大，風(fēng)速的影響次之，噴頭安裝高度其影響最小。因此在有風(fēng)條件下噴灌時(shí)，為了降低蒸發(fā)飄移損失，在保證較高噴灌均勻度的前提下，應(yīng)首先考慮降低噴頭工作壓力，以降低水滴碎裂概率，其次還可以降低噴頭安裝高度，以提高噴灌水分利用率。

表5 各因素對(duì)飄移損失影響的方差分析Tab.5 Analysis of variance between each factor and wind drift and evaporation losses

5 結(jié)論

(1)基于彈道理論和風(fēng)速模型，在有風(fēng)條件下模擬出水滴粒徑分布，通過(guò)對(duì)噴灑區(qū)域網(wǎng)格化處理，建立有風(fēng)條件下折射式噴頭水量分布計(jì)算方法。將模擬結(jié)果與FACI試驗(yàn)結(jié)果進(jìn)行了對(duì)比，結(jié)果表明水量分布模擬值與實(shí)測(cè)值基本相符，模擬與實(shí)測(cè)射程的相對(duì)偏差不大于10.03%，模擬的準(zhǔn)確性較高。

(2)計(jì)算了風(fēng)速、風(fēng)向、工作壓力、噴頭間距以及噴頭安裝高度等因素與噴灌均勻度的關(guān)系，采用SPSS軟件分析了以上5種因素對(duì)均勻度的影響，結(jié)果顯示風(fēng)速和風(fēng)向?qū)鶆蚨扔绊懖伙@著，噴頭間距對(duì)噴灌均勻度的影響最顯著，噴頭安裝高度次之，工作壓力影響較小。

(3)回歸了有風(fēng)條件下均勻度與噴頭間距的經(jīng)驗(yàn)?zāi)Ｐ停跃鶆蚨茸罡邽槟繕?biāo)，計(jì)算出最佳噴頭間距2.13～3.04 m，在此間距下，噴頭安裝高度和噴灌壓力的提升對(duì)噴灌均勻度的提高效果不大，因此為降低噴灌系統(tǒng)運(yùn)行成本，應(yīng)優(yōu)先采用低壓噴灌，并適當(dāng)降低噴頭安裝高度。

(4)風(fēng)速、工作壓力及安裝高度對(duì)蒸發(fā)飄移損失影響都較為顯著，其中工作壓力影響最大，風(fēng)速次之，噴頭安裝高度影響最小。因此有風(fēng)噴灌時(shí)，在保證較高噴灌均勻度的前提下，應(yīng)首先考慮降低噴頭工作壓力，以降低水滴碎裂概率，其次還可以降低噴頭安裝高度，以降低蒸發(fā)飄移損失。

1 李永沖, 嚴(yán)海軍, 徐成波, 等. 考慮水滴運(yùn)動(dòng)蒸發(fā)的噴灌水量分布模擬[J/OL]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2013, 44(7): 127-132. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20130722&journal_id=jcsam. DOI：10.6041/j.issn.1000-1298.2013.07.022. LI Yongchong, YAN Haijun, XU Chengbo, et al. Simulation of sprinkler water distribution with droplet dynamics and evaporation[J/OL]. Transactions of the Chinese Society for Agricultural Machinery,2013,44(7):127-132. (in Chinese)2 FUKUI Y, NAKANISHI K, OKAMURA S. Computer evaluation of sprinkler irrigation uniformity[J]. Irrigation Science, 1980, 2(1):23-32.

3 BURILLO G S, DELIRHASANNIA R, PLAYN E, et al. Initial drop velocity in a fixed spray plate sprinkler[J]. Journal of Irrigation & Drainage Engineering, 2013, 139(7):521-531.

4 LI Y, BAI G, YAN H. Development and validation of a modified model to simulate the sprinkler water distribution[J]. Computers & Electronics in Agriculture, 2015, 111:38-47.

5 張以升，朱德蘭，張林，等. 基于彈道軌跡方程的折射式噴頭水量分布計(jì)算模型[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2015,46(12):55-61.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20151208&flag=1. DOI：10.6041/j.issn.1000-1298.2015.12.008. ZHANG Yisheng, ZHU Delan, ZHANG Lin, et al.Water distribution model of fixed spray plate sprinkler based on ballistic trajectory equation[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2015,46(12):55-61. (in Chinese)

6 黃修橋, 廖永誠(chéng). 有風(fēng)條件下噴灌系統(tǒng)組合均勻度的計(jì)算理論與方法研究[J]. 灌溉排水, 1995, 14(1): 12-18.

7 CLARK G A, SRINIVAS K, ROGERS D H, et al. Measured and simulated uniformity of low drift nozzle sprinklers[J]. Transactions of the ASAE, 2003, 46(2): 321-330.

8 O’SHAUGHNESSY S A, URREGO Pereira Y F, EVETT S R, et al. Assessing application uniformity of a variable rate irrigation system in a windy location[J]. Applied Engineering in Agriculture, 2013, 29(4): 497-510.

9 SALVADOR R, BAUTISTA-CAPETILLO C, BURGUETE J, et al. A photographic method for drop characterization in agricultural sprinklers[J]. Irrigation Science, 2009, 27(4): 307-317.

10 FACI J M, SALVADOR R, PLAYAN E, et al. Comparison of fixed and rotating spray plate sprinklers[J]. Journal of Irrigation & Drainage Engineering, 2001, 127(4):224-233.

11 TARJUELO J M, MONTERO J, HONRUBIA F T, et al. Analysis of uniformity of sprinkle irrigation in a semi-arid area[J]. Agricultural Water Management, 1999, 40(2): 315-331.

12 VORIES E D, BERNUTH R D V, MICKELSON R H. Simulating sprinkler performance in wind[J]. Journal of Irrigation & Drainage Engineering, 1987, 113(1):119-130.

13 MONTERO J, TARJUELO J, CARRIN P. SIRIAS: a simulation model for sprinkler irrigation[J]. Irrigation Science, 2001, 20(2):73-84.

14 SEGINER I, NIR D, BERNUTH R D. Simulation of wind-distorted sprinkler patterns[J]. Journal of Irrigation and Drainage Engineering, 1991, 117(2): 285-306.

15 CHRISTIANSEN J E. Irrigation by sprinkling[R]. Berkeley: University of California,1942.

Uniformity of Fixed Spray-plate Sprinkler under Windy Condition Based on Ballistic Simulation

ZHANG Yisheng1ZHU Delan1,2SONG Bo1ZHANG Lin2

(1.CollegeofWaterResourcesandArchitecturalEngineering,NorthwestA&FUniversity,Yangling,Shaanxi712100,China2.InstituteofWaterSavingAgricultureinAridAreas,NorthwestA&FUniversity,Yangling,Shaanxi712100,China)

In order to simulate fixed spray-plate sprinkler water distribution and Christiansen uniformity (CU) under open field conditions, the method for calculating water distribution was proposed by applying droplet size distribution, simulation of droplet distribution utilize droplet trajectory computations and wind-distored model. The accuracy of water distributions of Nelson D3000 spray sprinkler was validated by comparing simulated with observed patterns, it was found that the software appeared to be an effective predictor of sprinkler performance in wind. Under this premise, the water distribution was simulated at different wind speeds with 4.76 mm nozzle diameters (24#), the factors, including wind speed, wind direction, sprinkler spacing, working pressure and riser height, which may affect CU, the wind drift and evaporation losses (WDEL) under different working conditions were analyzed. The results showed that sprinkler spacing had the most important influence on CU at 95% confidence interval, and it was followed by riser height and working pressure, wind speed and wind direction did not significantly affect CU. Increasing wind speed, working pressure or riser height indicated an increase in the OWDEL, and working pressure can dramatically affect the OWDEL on account of smashing large water droplet to pieces. Using this sprinkler, the lateral moving irrigation system can attain a high uniformity coefficient at 2.13～3.14 m sprinkler spacing along the lateral when the wind speed was lower than 6 m/s, and with this range of mounting spacing, increasing sprinkler riser height or working pressure cannot promote a higher CU distinctly, low working pressure was recommended considering little cost in the design of the sprinkler irrigation system. In addition, decreasing the riser height can reduce WDEL effectively, in order to increase water use efficiency, low height cooperated with plant height was reasonable.

fixed spray-plate sprinkler; sprinkler irrigation uniformity; water distribution; windy condition; ballistic theory

10.6041/j.issn.1000-1298.2017.02.012

2016-06-13

2016-07-22

“十二五”國(guó)家科技支撐計(jì)劃項(xiàng)目(2015BAD22B01-02)

張以升(1988—)，男，博士生，主要從事工程水力學(xué)研究，E-mail： hnsqzys@163.com

朱德蘭(1969—)，女，教授，博士生導(dǎo)師，主要從事節(jié)水灌溉技術(shù)研究，E-mail： dlzhu@126.com

S275.5

A

1000-1298(2017)02-0091-07