基于粒子成像瞬態(tài)測量技術(shù)的雨滴微物理特性及降雨動能研究

展小云，郭明航，趙 軍，史海靜，稅軍峰

?

基于粒子成像瞬態(tài)測量技術(shù)的雨滴微物理特性及降雨動能研究

展小云，郭明航※，趙 軍，史海靜，稅軍峰

（1. 西北農(nóng)林科技大學(xué)黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室，楊凌 712100； 2. 中國科學(xué)院水利部水土保持研究所，楊凌 712100）

雨滴微物理特性及降雨動能是揭示降雨物理本質(zhì)的重要特征量，亦是開展侵蝕定量分析與建立侵蝕量預(yù)報模型的基礎(chǔ)。采用粒子成像瞬態(tài)測量可視化技術(shù)觀測自然降雨雨滴，結(jié)合計算機視覺識別技術(shù)解算雨滴微物理特性參數(shù)，同時采用虹吸式自記雨量計記錄自然降雨降雨強度。研究表明：該次降雨雨滴以中等粒子為主，雨滴直徑均值為1.52 mm，降落末速度均值為3.47 m/s，其中直徑在1.00～3.00 mm范圍內(nèi)的雨滴占樣本總數(shù)的87.21%。雨滴直徑和降落末速度呈顯著的對數(shù)關(guān)系?；趯崪y的雨滴微物理特性和降雨強度估算降雨動能，該結(jié)果與傳統(tǒng)的經(jīng)驗?zāi)Ｐ凸浪憬Y(jié)果相對誤差均值為7.28%。該方法得到的降雨動能較以往的經(jīng)驗?zāi)Ｐ湍芨鎸嵉姆磻?yīng)雨滴降落過程中的做功大小，為準確計算降雨過程中雨滴所造成的濺蝕量奠定基礎(chǔ)。

侵蝕；圖像處理；物理特性；雨滴；粒子成像；降雨強度；降雨動能

0 引 言

雨滴在降落過程中受到萬有引力、空氣浮力、粘滯阻力、科里奧利力等多種因素的影響，使得雨滴在運動過程中呈現(xiàn)不同的大小、形狀、組成分布、降落速度和動能等微物理特性[1]。在節(jié)水灌溉中，雨滴微物理特性是評價噴灌系統(tǒng)工程質(zhì)量的重要指標之一，在噴灌過程中，直徑小的水滴由于飄逸蒸發(fā)損失，降低水分利用效率；直徑大的水滴降落過程中動能較大，從而傷害作物并導(dǎo)致土壤板結(jié)，降低土壤入滲性能，加劇土壤濺蝕量[2-3]；在水土保持中，雨滴大小和速度是研究土壤水力侵蝕動力學(xué)過程、計算陸面過程中降雨侵蝕力以及制定水土保持工程措施的的基本依據(jù)[4-6]；在物理學(xué)中，雨滴的形狀和折射率是研究雨滴對電磁波散射和衰減影響的重要前提[7-8]；在氣象科學(xué)中，雨滴的尺度、形狀、速度及其譜分布是評估人工模擬降雨效果標準之一，亦是了解自然降雨發(fā)展趨勢，揭示降雨機制的重要參數(shù)[9-10]。如何測量雨滴物理特性及降雨動能成為節(jié)水灌溉、土壤侵蝕、氣象科學(xué)等研究工作的熱點問題[11-12]。

早期國內(nèi)外學(xué)者主要通過濾紙色斑法、面粉團法、浸潤法、動力學(xué)法等監(jiān)測雨滴的微物理特性，但是上述方法普遍存在自動化程度低，實際應(yīng)用局限性大等問題[13-16]。例如，應(yīng)用較多的濾紙色斑法后期數(shù)據(jù)處理冗繁，并且不能對雨滴進行連續(xù)觀測。浸潤法中由于小粒子無法沖破表面張力而浸入液體，使得該方法不能監(jiān)測自然降雨中的小雨滴。19世紀60年代，Clardy等[17]第一次使用雨滴譜儀測量雨滴大小和速度。至此人們開始利用基于光電、聲電等原理研制的雨滴譜儀進行雨滴微物理特性的監(jiān)測，雨滴譜儀得到快速的發(fā)展和應(yīng)用[18-21]。目前應(yīng)用最為廣泛的是Parsivel激光雨滴譜儀和二維雨滴譜儀（2D Vidio Didtromet, 2-DVD）[4,22]。兩者均操作簡單，并可快速地獲取豐富的雨滴微物理特性信息，但是Parsivel激光雨滴譜儀無法區(qū)分同時下落的雨滴粒子，有很大的重疊誤差；2-DVD雨滴譜儀則容易受氣流影響，尤其在強風(fēng)條件下雨滴變形導(dǎo)致雨滴下降速度測量誤差較大。

可見，上述各種方法都存在工作冗繁和誤差較大的問題，并且適用于實驗室的各種雨滴特性測試方法都不能進行實時測量。粒子成像瞬態(tài)測量可視化技術(shù)采樣空間大，測量誤差小，并且能夠快速、準確、實時地觀測雨滴大小、運動速度等雨滴微物理特性。本文利用該技術(shù)和裝置研究自然降雨雨滴微物理特性、降雨動能特征，分析雨滴直徑與降落末速度關(guān)系，在此基礎(chǔ)上基于實測的雨滴微物理特性估算降雨動能，與傳統(tǒng)的統(tǒng)計模型比對，驗證統(tǒng)計模型的準確性和適用性。

1 材料與方法

1.1 試驗環(huán)境

楊凌地處陜西省關(guān)中平原中西部，東經(jīng)107°59′～108°09′，北緯34°14′～34°24′，海拔418.0～540.1 m，地勢北高南低。氣候?qū)倥瘻貛О霛駶櫄夂?，四季分明，年平均氣?2.9 ℃，最高氣溫42.0 ℃，最低氣溫?19.4 ℃，全年無霜期221 d。年均降水量637.6 mm，多集中在7～10月，占多年平均降水量的60%，年平均蒸發(fā)量884 mm。觀測地點位于中國科學(xué)院水利部水土保持研究所科研樓。由于試驗是在無風(fēng)條件下進行的，風(fēng)速對雨滴形狀、速度等微物理特性的影響忽略不計。

1.2 雨滴觀測裝置

采用粒子成像的瞬態(tài)測量可視化技術(shù)測量雨滴微物理特性參數(shù)。該技術(shù)主要是利用CCD工業(yè)像機對降落雨滴進行快拍（曝光時間短）和慢拍（曝光時間長）操作，采集同一時刻的雨滴“靜止”和“拖尾”影像，確定其在特定的時間序列上的空間位置和形狀，進而解算雨滴微物理特性參數(shù)。粒子成像測量系統(tǒng)主要由3個子系統(tǒng)組成，即投影系統(tǒng)、采集系統(tǒng)和采集控制系統(tǒng)，具體功能結(jié)構(gòu)如圖1所示。

1.投影機箱 2.點光源 3.固定架 4.菲涅爾透鏡軌道 5.菲涅爾透鏡 6.菲涅爾透鏡固定框 7.投影機前面板 8.投影機箱前蓋9.采集機箱 10.CCD支架 11.幕布框 12.采集機前面板 13.快拍相機 14.慢拍相機 15.投影幕布 16.控制機箱 17.電路板 18.控制機下蓋 19.交換機 20.出線面板 21.光軸 22.雨滴 23.平行光 24.圓形標靶

該系統(tǒng)主要包括6個工作單元，即點光源、菲涅爾透鏡、投影幕布、CCD工業(yè)相機、相機同步控制器和計算機終端。其中，點光源采用一種色溫為4 300 K左右的飛利浦銀戰(zhàn)士鹵素?zé)簦环颇鶢柧酃忡R將光束擴束后投射到投影幕布，二者構(gòu)成有效采樣空間；CCD工業(yè)相機為德國BASLER acA 1600-20 gm相機，分辨率為1 624 pixels× 1 234 pixels，相機采集頻率設(shè)定為20 幀/s，曝光時間分別設(shè)定為1/4 000和1/500 s，CCD工業(yè)相機將接受到的光強信號轉(zhuǎn)換為數(shù)字圖像信號；采集控制系統(tǒng)由D-Link千兆五端口交換機、同步觸發(fā)板、四芯航空接頭、串口轉(zhuǎn)換裝置等組成，實現(xiàn)圖像數(shù)據(jù)的實時采集、管理和高效率傳輸?shù)?；用戶終端可以實現(xiàn)雨滴形狀、運動軌跡等的記錄和存儲，據(jù)此可以計算和分析雨滴直徑、雨滴降落末速度等。

光學(xué)測量單元采用滿畫幅像素拍攝，像元寬度= 4.1m，用于感應(yīng)點光源產(chǎn)生的平行光。相機使用焦距為8 mm工業(yè)鏡頭，投影幕布距離相機中心的距離為320 mm，菲涅爾透鏡有效面積為230 mm × 210 mm，考慮到相機分辨率則平行光的投影寬度為202.38 mm，同理可以算出平行光的投影高度為266.34 mm。此外，如果投影幕布和菲涅爾透鏡的距離過大或者過小，使得發(fā)射到幕布上的平行光強度較弱或者較強，最終使得幕布背景過于黑暗或者亮度過高，均不能正常顯示雨滴影像，只有保證雨滴在幕布上的影像灰度和幕布本身的灰度差異最大化才能有效獲取雨滴影像。綜合以上2點，將系統(tǒng)的平面采樣空間設(shè)置為200 mm × 200 mm，系統(tǒng)在垂直于投影屏幕方向上的開口距離也為200 mm。該系統(tǒng)相對于目前的2-DVD系統(tǒng)采樣空間增大，提高了采樣的代表性，彌補了2-DVD對大粒子探測能力較差的缺點。

1.3 雨滴影像解譯軟件系統(tǒng)

為了獲取清晰的雨滴影像，目前的研究提出了不同的雨滴微觀特征提取方法[23-24]。本研究針對上述監(jiān)測系統(tǒng)設(shè)計了一套獨有的雨滴微物理特性提取算法，即采用計算機識別技術(shù)對原始雨滴影像進行數(shù)字影像糾正、計算機圖像深度處理等方法剔除噪音，進而勾畫二值化的雨滴影像，在此基礎(chǔ)上匹配快拍和慢拍得到的同名雨滴影像，進而為計算雨滴微物理特性參數(shù)奠定基礎(chǔ)。經(jīng)過大量的對比和計算最終確定獲取的初始雨滴影像中的兩類主要噪音為高斯噪音和尖峰噪音。針對以上兩類噪音設(shè)計了一套噪音處理算法并研發(fā)相應(yīng)的軟件解譯系統(tǒng)，如圖2所示。

圖2 雨滴影像解譯系統(tǒng)

具體算法為：1）統(tǒng)計濾波算法剔除絕大部分的背景噪音；2）滾動濾波算法剔除影像中孤立的噪音，并最大程度地保留影像細節(jié)；3）鄰域平滑濾波算法剔除少量殘留的噪音；4）全局二值化方法勾畫二值化的雨滴影像；5）利用定位標靶信息進行同名雨滴的辨識，噪音剔除的過程影像如圖3所示。雨滴影像包含豐富的細節(jié)信息，如圖4所示，基于處理后的清晰的雨滴影像可以得到雨滴的精細微物理特性。

圖3 噪音剔除過程影像

1.雨滴陰影區(qū) 2. 雨滴外包絡(luò)線 3.雨滴形態(tài)對稱軸（運動方向） 4.雨滴拖尾長度

1.4 試驗設(shè)置

為了檢測上述觀測裝置的準確性，進行鋼珠灑落試驗，觀測鋼珠的直徑和降落末速度，用鋼珠觀測直徑和降落末速度分別與實際直徑和理論降落末速度對比分析以檢測該方法的準確性，結(jié)果表明該方法觀測誤差小，可有效、快速地觀測雨滴微物理特性[25]。本研究利用上述觀測系統(tǒng)在陜西楊凌區(qū)進行外場觀測試驗采集自然降雨的雨滴影像，共采集得到8 448幀圖像，其中2 762幅有效雨滴圖像，因此雨滴對完整捕獲概率為32.7%。同時，試驗采用上海儀器廠生產(chǎn)的SJ型虹吸式自記雨量計測定降雨量，將雨量計布設(shè)在開闊無遮擋物區(qū)域，避免降雨受到建筑物和高大樹木的影響。

2 觀測指標

2.1 雨滴直徑

描述雨滴直徑常用的特征參數(shù)為平均直徑，包括周長變換平均直徑、軌跡平均直徑、體積平均直徑、等圓平均直徑以及幾何平均直徑。本研究中采用幾何平均直徑，相對于其他平均直徑的算法而言，幾何平均直徑計算流程簡單，在平面上充分利用了幾何均值的特點，誤差較小。幾何平均直徑的計算方法是根據(jù)在“靜止”的雨滴輪廓上尋找距離最遠的2個點，其連線作為雨滴最大直徑，經(jīng)該直徑的中點做垂線，交于雨滴輪廓上兩點，該兩點之間的距離即為雨滴第二直徑，最大直徑和第二直徑的幾何均值即為雨滴幾何平均直徑，如圖5a所示，即

式中為雨滴的幾何平均直徑，mm；1為長軸的長度，mm；2為短軸的長度，mm。

注：1和2分別為長軸和短軸的長度，1和2分別為慢拍相機和快拍相機的曝光時間，為雨滴的幾何平均直徑，為雨滴拖尾長度。

Notes:1and2are length of long axis and short axis,1and2are exposure time of snapshot and slowshot camera,is the geometric mean diameter of raindrop,is the tailed length of raindrop.

圖5 雨滴直徑和移動距離算法示意圖

Fig.5 Schematic diagram of diameter and movement distance algorithms of raindrop

2.2 雨滴降落末速度

雨滴在高空的降落過程雖然是一個變加速的過程，但是由于其受力的復(fù)雜性，導(dǎo)致雨滴在到達地面之前必然達到一個收尾速度，加之本研究采用的是瞬態(tài)粒子成像可視化監(jiān)測技術(shù)，因此，雨滴近地面的運動可視為勻速運動。本研究采用外包絡(luò)線提取算法計算雨滴的降落末速度，即利用“拖尾”雨滴影像與同名“靜止”雨滴影像之間的相位差除以快拍和慢拍相機的曝光時間差得到。該算法的計算原理如圖5b所示，即

式中()為雨滴降落末速度，m/s；為雨滴拖尾長度，mm；1為慢拍相機的曝光時間，s；2為快拍相機的曝光時間，s。

2.3 降雨動能

雨滴動能是根據(jù)雨滴大小及其組成計算求得。單個雨滴的動能主要取決于雨滴直徑和降落末速度。本文把雨滴視為一個理想的球體，用下式計算單個雨滴的動能，即

式中為單個雨滴的動能，J；為標準狀態(tài)下水的密度，kg/mm3。

本研究中采用Marshall等[26]提出的M-P指數(shù)描述雨滴分布，即

式中()為單位尺度間隔、單位體積內(nèi)的雨滴數(shù)，個/(m3·mm)；0為約等于8 000的常數(shù)；為降雨強度（I, mm/h）的冪函數(shù)，表達為

在靜風(fēng)的條件下，Hall等[27]的研究表明單位時間單位地表的雨滴數(shù)′()可表達為降落末速度（()）和雨滴譜（()）的乘積形式，即

聯(lián)合（3）式、（4）式、（5）式和（6）式可得降雨總動能（E()）的如下表達式，即

式中E()為單位時間內(nèi)單位地表所接收到的降雨總動能，J/(m2·min)。

3 結(jié)果與分析

3.1 雨滴微物理特性參數(shù)

為了描述雨滴微物理特性的統(tǒng)計特征，首先計算了樣本的偏度值（）和豐度值（），并進行K-S檢驗，結(jié)果表明雨滴直徑和降落末速度均呈偏態(tài)分布（>0.05），對其取常用對數(shù)后，則呈正態(tài)分布，說明雨滴直徑和降落末速度符合對數(shù)正態(tài)分布，所以用幾何均值比算術(shù)均值更能真實地反映樣本的統(tǒng)計特征。從表1可知，雨滴的直徑和降落末速度存在很大的變異性，其變異系數(shù)分別為1.33%和2.30%。其中，雨滴直徑變化范圍為0.42～4.86 mm，幾何均值為1.52 mm；雨滴降落末速度變化范圍為0.24～9.33 m/s，幾何均值為3.47 m/s，并且雨滴降落末速度分布存在明顯的右偏（=1.20），表明其高值較多。圖6將雨滴直徑和降落末速度從最小值到最大值均分成8類，橫坐標為雨滴微物理特性參數(shù)，縱坐標為每類雨滴所占的百分比。由圖6可知，直徑在1.00～3.00 mm范圍內(nèi)的雨滴占樣本總數(shù)的87.21%，而小雨滴（直徑<1.00 mm）和大雨滴（直徑>3.00 mm）僅占樣本總數(shù)的12.79%，可見本次降雨雨滴尺寸以中等粒子為主。雨滴降落末速度集中在1.00～5.00 m/s，占樣本總數(shù)的67.85%，具有極端降落末速度的雨滴（>7.00 m/s）僅占樣本總數(shù)的2.64%。

表1 雨滴直徑和降落末速度的統(tǒng)計特征

圖6 雨滴微物理特性參數(shù)分布

3.2 速度譜分布

雨滴在大氣中受重力、空氣阻力和浮力的綜合作用，其降落過程是一個變加速過程，隨著時間變化加速度逐漸減小，當降落高度足夠大時，雨滴受力達到一個平衡狀態(tài)，此后做勻速運動，此時雨滴的速度即為降落末速度[28]。本研究中，根據(jù)最小二乘法原理對雨滴直徑和降落末速度進行回歸分析，得到雨滴降落末速度與雨滴直徑間呈顯著的對數(shù)關(guān)系（2=0.848，<0.001）（圖7a）?？梢?，雨滴降落末速度隨雨滴直徑的增加先急劇增加，當雨滴直徑大于1.5 mm時，雨滴降落末速度的增加趨勢逐漸平緩（圖7a），這可能是由于較大粒徑的雨滴受到表面張力的影響，改變了其外部形狀，使其所受的空氣阻力增加，從而在一定程度上減弱了雨滴降落速度。圖7b給出了Best模型[29]、Atlas 模型[30]、Uplinger 模型[31]、Gossard[32]模型和Brandes模型[33]及本研究中的實測數(shù)據(jù)的對比情況。由圖7b可知，傳統(tǒng)的經(jīng)驗?zāi)Ｐ团c本研究中的實測數(shù)據(jù)存在較好的一致性，均表現(xiàn)為雨滴降落末速度隨雨滴直徑的增大呈顯著的增加趨勢。相對于傳統(tǒng)的經(jīng)驗?zāi)Ｐ?，本文擬合的速度略小，這可能是由于雨滴飛濺或破碎后經(jīng)過采樣空間造成的，也可能是由于儀器本身采樣空間較大，在采樣系統(tǒng)和投影系統(tǒng)之間出現(xiàn)擾動，從而影響雨滴的運動軌跡，改變了原有的雨滴降落速度。

圖7 雨滴速度譜分布

3.3 降雨動能特征

從表2可以看出，降雨動能隨降雨強度的增大呈顯著的增加趨勢，這與國內(nèi)外基于降雨強度估算降雨動能的經(jīng)驗?zāi)Ｐ褪且恢碌?。基于實測的降雨動能與吳光艷等[34]（2011）和周佩華等[35]估算的降雨動能十分接近，相對誤差分別在3.26%～4.84%和1.23%～8.00%之間，但是明顯高于江忠善等[36]得出的降雨動能，相對誤差最大值高達16.94%，最小相對誤差也高達8.00%。圖8為幾種經(jīng)驗公式曲線與本研究實測數(shù)據(jù)的對比，可見，本研究中降雨強度與降雨動能的關(guān)系曲線位于吳光艷等[34]和周佩華等[35]得出的關(guān)系曲線之間，較江忠善等[36]得出的關(guān)系曲線差異明顯。以上幾種基于降雨強度的經(jīng)驗關(guān)系曲線多是利用色斑法得出的，該方法中雨滴濺落在濾紙上形成一群形狀極不規(guī)則的帶毛刺的圖案，使得依據(jù)色斑形狀估計雨滴粒徑在很大程度上依靠主觀經(jīng)驗，使得該方法的誤差高達6%～14%，并不能真實地反映雨滴微物理特性[3]。可見，傳統(tǒng)的基于降雨強度估算降雨動能的經(jīng)驗?zāi)Ｐ腿狈ζ者m性，而本研究中基于雨滴直徑和降落末速度估算的降雨動能的方法在很大程度上彌補了降雨動能統(tǒng)計模型的不足。

表2 降雨動能統(tǒng)計

圖8 降雨強度與降雨動能關(guān)系

4 結(jié) 論

1）利用粒子成像瞬態(tài)測量可視化技術(shù)進行外場觀測試驗，研究楊凌地區(qū)自然降雨雨滴微物理特性及其降雨動能特征。本次降雨強度較小，87.21%以上的雨滴直徑集中在1.00～3.00 mm范圍內(nèi)，雨滴的平均直徑為1.52 mm，降落末速度均值為3.47 m/s。

2）通過擬合雨滴直徑和降落末速度的關(guān)系分析此次自然降雨的速度譜特征，發(fā)現(xiàn)雨滴降落末速度與雨滴直徑間呈顯著的對數(shù)關(guān)系，相對于以往的經(jīng)驗?zāi)Ｐ停疚臄M合的速度略小，但是本文的擬合關(guān)系與以往的經(jīng)驗?zāi)Ｐ途哂休^好的一致性。究其原因可能是由于較大粒徑的雨滴在降落過程中受表面張力的影響改變了其外部形狀，使其所受的空氣阻力增加，從而減弱了雨滴降落速度的增加程度。

3）通過對降雨強度和雨滴微物理特性的觀測，初步估算了降雨動能，并建立了具有本地化特征的降雨動能和降雨強度之間的統(tǒng)計關(guān)系，即k()=26.820.34?；诮涤陱姸鹊慕涤陝幽芙y(tǒng)計模型具有很大的區(qū)域局限性，通過實測的雨滴微物理特性估算降雨動能，能夠檢驗統(tǒng)計模型的準確性和適用性，為建立侵蝕量預(yù)報數(shù)學(xué)模型提供有關(guān)降雨侵蝕力的依據(jù)。

[1] Michaelides S. Precipitation: Advances in Measurement, Estimation, and Prediction[M]. Berlin: Springer, 2008.

[2] 鞏興暉，朱德蘭，張林，等. 基于2DVD的非旋轉(zhuǎn)折射式噴頭水滴直徑分布規(guī)律[J]. 農(nóng)業(yè)機械學(xué)報，2014，45(8)：128－148.

Gong Xinghui, Zhu Delan, Zhang Lin, et al. Drop size distribution of fixed spray-plate sprinklers with two-dimensional video disdrometer[J]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(8): 128－148. (in Chinese with English abstract)

[3] Mark W S, Timothy R G, James C A. Tillage effects on soil hydraulic properties in space and time: State of the science[J]. Soil and Tillage Research, 2008, 99(1): 4－48.

[4] 朱亞喬，劉元波. 地面雨滴譜觀測技術(shù)及特征研究進展[J]. 地球科學(xué)進展，2013，28(6)：685－694.

Zhu Yaqiao, Liu Yuanbo. Advances in measurement techniques and statistics features of surface raindrop size distribution[J]. Advances in Earth Science, 2013, 28(6): 685－694. (in Chinese with English abstract)

[5] 付玉，李光錄，鄭騰輝，等. 雨滴擊濺對耕作層土壤團聚體粒徑分布的影響[J]. 農(nóng)業(yè)工程學(xué)報，2017，33(3)：155－160.

Fu Yu, Li Guanglu, Zheng Tenghui, et al. Effects of raindrop splash on aggregate particle size distribution of soil plough layer[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(3): 155－160. (in Chinese with English abstract)

[6] 趙龍山，梁心藍，張青峰，等. 裸地雨滴濺蝕對坡面微地形的影響與變化特征[J]. 農(nóng)業(yè)工程學(xué)報，2012，28(19)：71－77.

Zhao Longshan, Liang Xinlan, Zhang Qingfeng, et al. Variation characteristics and effects of splash erosion on slope micro-relief in bare fields[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(19): 71－77. (in Chinese with English abstract)

[7] 劉西川，高太長，劉磊，等. 雨滴微物理特征研究及測量技術(shù)進展[J]. 地球科學(xué)進展，2013，28(11)：1217－1226.

Liu Xichuan, Gao Taichang, Liu Lei, et al. Advances in microphysical features and measurement techniques of raindrops[J]. Advances in Earth Science, 2013, 28(11): 1217－1226. (in Chinese with English abstract)

[8] 劉西川，高太長，劉磊，等. 基于粒子成像測速技術(shù)的雨滴微物理特性研究[J]. 物理學(xué)報，2014，63(2)：469－475.

Liu Xichuan, Gao Taichang, Liu Lei, et al. Research on microphysical property of raindrops based on particle imaging velocimetry technology[J]. Acta Physica Sinica, 2014, 63(2): 469－475. (in Chinese with English abstract)

[9] Szakall M, DiehlK，Mitra S K, et al. A wind tunnel study on the shape, oscillation, and internal circulation of large raindrops with sizes between 2.5 and 7.5 mm[J]. Journal of Atmospheric Sciences, 2009, 66(3): 755－765.

[10] 劉紅燕，雷恒池. 基于地面雨滴譜資料分析層狀云和對流云降水的特征[J]. 大氣科學(xué)，2006(4)：693－702.

Liu Hongyan, Lei Hengchi. Characteristics of rain from stratiform versus convective cloud based on the surface raindrop data [J]. Chinese Journal of Atmospheric Sciences, 2006(4): 693－702. (in Chinese with English abstract)

[11] Martinez-Mena M, Castillo V, Albaladejo J. Relations between interrill erosion processes and sediment particle size distribution in a semiarid Mediterranean area of SE of Spain[J]. Geomorphology, 2002, 45(3): 261－275.

[12] Warrington D N, Mamedov A I, Bhardwag A K, et al. Primary particle size distribution of eroded material affected by degree of aggregate slaking and seal development[J]. European Journal of Soil Science, 2009, 60(1): 84－93.

[13] Bentley W A. Studies of raindrops and raindrop phenomena[J]. Monthly Weather Review, 1904, 32: 450-456.

[14] 李紅，任志遠，袁壽其，等. 高度對色斑法測量雨滴粒徑影響的試驗研究[J]. 中國農(nóng)村水利水電，2006(1)：16－17.

Li Hong, Ren Zhiyuan, Yuan Shouqi, et al. Experimental study on the height effect on raindrop size in splash method [J]. China Rural Water and Hydropower, 2006(1): 16－17. (in Chinese with English abstract)

[15] 舒若杰，高建恩，吳普特，等. 基于CorelDRAW軟件的小流域模型雨滴測量試驗研究[J]. 農(nóng)業(yè)工程學(xué)報，2006, 22(11)：44－46.

Shu Ruojie, Gao Jian’en, Wu Pute, et al. Measurement of raindrop distribution of watershed model using CorelDRAW[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2006, 22(11): 44－46. (in Chinese with English abstract )

[16] 倪際梁，何進，李洪文，等. 便攜式人工模擬降雨裝置的設(shè)計與率定[J]. 農(nóng)業(yè)工程學(xué)報，2012，28(24)：78－84.

Ni Jiliang, He Jin, Li Hongwen, et al. Design and calibration of portable rainfall equipment of artificial simulation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(24): 78－84. (in Chinese with English abstract)

[17] Clardy D E, Tolberrt C W. Electronic disdrometer[J]. Review of Scientific Instruments. 1961, 32(8): 916－920.

[18] Kinnellp I A. Some observations on the Joss-Waldvogel rainfall disdrometer[J]. Journal of Applied Meteorology, 1976, 15(5): 499－502.

[19] Sheppardb E, Joe P I. Comparison of raindrop size distribution measurements by a Joss-Waldvogel Disdrometer, a PMS 2DGSpectrometer，and a POSS Doppler Radar[J]. Journal of Atmospheric and Oceanic Technology, 1994, 11: 874－887.

[20] Wolff D B, Wolff K R. Rain gauge and disdromrter measurements during the Keys Area Microphysics Project (KAMP)[J]. Journal of Atmospheric and Oceanic Technology, 2003, 20: 1460－1477.

[21] 萬和躍，濮江平，劉思瑤. 雨滴譜探測技術(shù)發(fā)展及應(yīng)用概述[J]. 科技創(chuàng)新導(dǎo)報，2015，31：166－167.

[22] 韓文霆，曹培，劉文帥. 人工模擬降雨系統(tǒng)噴頭噴灑雨滴特性測試[J]. 農(nóng)業(yè)機械學(xué)報，2014，45(2)：56－61.

Han Wenting, Cao Pei, Liu Wenshuai. Raindrop characteristics of sprinklers for artificial rainfall system[J]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(12): 56－61. (in Chinese with English abstract)

[23] 朱磊，曹治國，肖陽，等. 基于圖像自動識別技術(shù)的雨滴譜測量方法[J]. 計算機應(yīng)用研究，2015，32(9)：2866－2869，2873.

Zhu Lei, Cao Zhiguo, Xiao Yang, et al. Automatic measurement method of raindrop size distribution based on computer vision[J]. Application Research of Computers，2015, 32(9): 2866－2869, 2873. (in Chinese with English abstract)

[24] 倫向敏，侯一民. 運用迭代最大熵算法選取最佳圖像分割閾值[J]. 計算機工程與設(shè)計，2015, 36(5)：1265－1268，1289.

Lun Xiangmin, Hou Yimin. Optical threshold slection for image segmentation utilizing entropy-max algorithm[J]. Computer Engineering and Design, 2015, 36(5): 1265－1268, 1289. (in Chinese with English abstract)

[25] 郭明航，展小云，趙軍，等. 雨滴物理特性的粒子成像測量技術(shù)研究[J]. 農(nóng)業(yè)機械學(xué)報，2015，46(9)：144－150.

Guo Minghang, Zhan Xiaoyun, Zhao Jun, et al. Measurement of raindrop physical properties with particle imaging measurement technology[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015，46(9): 144－150. (in Chinese with English abstract)

[26] Marshall J S, Palmer W M. The distribution of raindrops with size[J]. Journal of Meteorology, 1948(5): 165－166.

[27] Hall R L, Calder I R. Drop size modification by forest canopies-measurements using a disdrometer[J]. Journal of Geophysical Research-Atmospheres, 1993, 98(D10): 18465－18470.

[28] 劉雅君. 雨滴下落的收尾速度[J]. 大學(xué)物理，2001，19(4)：45－46.

Liu Yajun. Terminal velocity of the falling raindrops[J]. College Physics, 2001, 19(4): 45－46. (in Chinese with English abstract)

[29] Best A C. Empirical formulate for the terminal velocity of water drops falling through the atmosphere[J]. Quarterly Journal of the Royal Meteorologyical Society, 1950, 76: 302－311.

[30] Atlas D, Srivastava R C, Sekhon R S. Doppler radar characteristics of precipitation at vertical incidence[J]. Reviews of Geophysics, 1973, 11(1): 1－35.

[31] Uplinger W G. A new formula for raindrop terminal velocity[C]∥20th Radar Meteorology Conference. Boston: American Meteorological Scoiety, 1977: 389－391.

[32] Gossard E E, Strauch R G, Welsh D C, et al. Cloud layers, particle identification, and rain-rate profiles from ZRVf measurements by clear-air doppler radars[J]. Journal of the Atmospheric and Oceanic Technology, 1992, 9(2): 108－119.

[33] Brandes E A, Zhang G, Vivekanandan J. Experiments in rainfall estimation with a polarimetric radar in a subtropical environment[J]. Journal of Applied Meteorology, 2002,41(6): 674－684.

[34] 吳光艷，吳發(fā)啟，尹武君，等. 陜西楊凌天然降雨雨滴特性研究[J]. 水土保持研究，2011，18(1)：48－51.

Wu Guangyan, Wu Faqi, Yin Wujun, et al. Study on characteristics of natural rainfall in Yangling, Shaanxi Province[J]. Research of Soil and Water Conservation，2011，18(1): 48－51. (in Chinese with English abstract)

[35] 周佩華，竇葆璋，孫清芳. 降雨能量的試驗研究初報[J]. 水土保持通報，1981，1(1)：51－60

[36] 江忠善，宋文經(jīng)，李秀英. 黃土地區(qū)天然降雨雨滴特性研究[J]. 中國水土保持，1983(3)：32－36.

Microphysical features of raindrop and rainfall energy based on particle imaging transient measurement technology

Zhan Xiaoyun, Guo Minghang※, Zhao Jun, Shi Haijing, Shui Junfeng

(1.,,712100,;2.,,712100,)

Microphysical features of raindrop and rainfall energy are the key parameters for study of rainfall physics, which also have great significance in quantitative analysis of soil erosion and in soil erosion prediction model. The existing measurement methods include splash method, immersion method and so on, but most of them have many disadvantages, such as, labor and time consuming, poor real-time response, low precision and so on. Therefore, a new method should be developed urgently. In order to obtain physical properties of raindrops, in this study, the particle imaging transient visual measurement technology, light field and imaging system were designed, image identification, extraction and measurement were investigated, and finally the particle imaging transient visual measurement technology and computer vision identification technology were used to obtain the microphysical features of natural raindrop. The principle of the system could be described as below: a Fresnel lens was installed in the front of the light source. When the lights were passing through the Fresnel lens, the lights from point source would become parallel lights, which would shine on a projecting screen, and then the raindrop would project on the screen during its falling. Specifically, the system consisted of three parts: projection system, image capture system, and image control system. In the image capture system, two cameras were used to capture the raindrop image, one with a fast speed to obtain static image, and the other worked slowly to capture the tailed image. Based on the two images, we calculate the diameter and the fall velocity of raindrops. In order to obtain a clear raindrop image, we must remove the noise in the images. Basically, the image noise removal involved four steps. First, it was statistical filtering; second, the rolling filtering; third, the smooth filtering, and finally, it was the image binarization. Based on the static image and the corresponding tailed image mentioned above, geometric mean value algorithm and outer contour algorithm were used to calculate the diameter and fall velocity of raindrop. In our previous research, we found that the measurement technology had small relative error and it was suitable for the measurement of microphysical features of raindrop. Meanwhile, rainfall intensity was recorded by siphonic pluviograph. The results showed that medium-sized particles were the predominant contributor in the single rainfall. Raindrop diameter and fall velocity in our study were, on average, 1.52 mm and 3.47 m/s, respectively. Specifically, the proportion of raindrops with diameter ranged from 1.00 to 3.00 mm was up to 87.21%. Fall velocity was strong logarithmically related to raindrop diameter, and more precisely, fall velocity grew rapidly with an increase in rainfall diameter when the diameter was below 1.5 mm. As the raindrops fatten, the growth rate of fall velocity was reduced. Moreover, rainfall energy calculated in the present study was compared with the classic statistical model, and the relative error was averaged as 7.28%. In all, microphysical features of raindrop and rainfall energy can be measured precisely by the technique in this study, which sets the basis of estimating rainfall splash erosion.

erosion; image processing; physical properties; raindrop; particle imaging; rainfall intensity; rainfall energy

10.11975/j.issn.1002-6819.2018.02.014

P426; S157

A

1002-6819(2018)-02-0107-07

2017-09-20

2018-01-04

國家自然科學(xué)基金項目（41571269；41503078）；黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室重要方向創(chuàng)新項目（A314021403-C3）

展小云，女，博士，主要從事水土流失過程與水土保持研究。Email：zhanxiaoyun2005@163.com

郭明航，男，陜西乾縣人，研究員，主要從事科研信息化與裝備研究。Email：mhguo@ms.iswc.ac.cn

展小云，郭明航，趙 軍，史海靜，稅軍峰. 基于粒子成像瞬態(tài)測量技術(shù)的雨滴微物理特性及降雨動能研究[J]. 農(nóng)業(yè)工程學(xué)報，2018，34(2)：107－113. doi：10.11975/j.issn.1002-6819.2018.02.014 http://www.tcsae.org

Zhan Xiaoyun, Guo Minghang, Zhao Jun, Shi Haijing, Shui Junfeng. Microphysical features of raindrop and rainfall energy based on particle imaging transient measurement technology[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(2): 107－113. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2018.02.014 http://www.tcsae.org