蔣建清，程 超，蔡晶垚，夏 平，王 亮，肖鄔旺

(湖南城市學院 土木工程學院，湖南 益陽 413000)

一種自循環(huán)式人工降雨模擬裝置降雨特性的試驗研究

蔣建清，程 超，蔡晶垚，夏 平，王 亮，肖鄔旺

(湖南城市學院 土木工程學院，湖南 益陽 413000)

結(jié)合天然降雨特性，研制了一套適于巖土工程降雨模型試驗和水土保持的自循環(huán)式人工降雨模擬裝置，并對其降雨特性進行了率定試驗，得到了該裝置降雨強度與控水閥門開度、降雨動能與降雨強度、雨滴平均直徑與濾紙色斑平均直徑等特征量之間的計算關(guān)系式及雨滴直徑累積頻率分布曲線．結(jié)果表明：該裝置實現(xiàn)了雨水自循環(huán)利用，可模擬的降雨強度范圍為4.65 mm/h-406.05 mm/h，降雨均勻度超過85 %，雨滴粒徑和雨滴終點落地速度分別在0.94-2.16 mm、3.69-7.14 m/s之間，模擬降雨與天然降雨的相似度好，為實驗室為在短期內(nèi)開展大批量降雨模型試驗和獲取實驗數(shù)據(jù)提供了保障．

雨水自循環(huán)；降雨模擬裝置；降雨特性；研制；率定

在實驗室使用降雨模擬裝置進行雨水入滲巖土工程穩(wěn)定性、土壤入滲、坡面產(chǎn)流和土壤侵蝕等方面的研究，不僅具有經(jīng)濟性、便捷性、可控性、重復(fù)性等優(yōu)點，還可以控制實驗進程，縮短實驗周期．目前，降雨模擬裝置主要有4種類型，即管網(wǎng)式、針管式、懸線式和噴嘴式[1]．比如，徐向舟等[2]研制的 SX2002管網(wǎng)式降雨模擬裝置，其降雨高度為3.5 m，降雨面積為3.5 m×2.5 m，降雨強度范圍為1.0-4.0 mm/min，降雨均勻度超過80 %，雨滴粒徑為0-3 mm；孫愷等[3]研制針管式人工降雨裝置，其降雨高度為1.6 m，降雨面積為196 250 mm2，降雨強度范圍0-15 mm/h，降雨均勻度超過90 %，雨滴粒徑為0-10 mm；Bowyer-Bower[4]研制的懸線式降雨模擬器，其采用了內(nèi)徑0.7 mm、外徑2-3 mm長15 mm的聚乙烯管，雨滴中徑為2-3 mm；霍云梅等研制的旋轉(zhuǎn)下噴式降雨裝置[5]、蘇溦娜等研制的噴嘴式降雨模擬裝置[6]和張晶玲等研制的搖臂上噴式降雨模擬裝置[7]．各種類型降雨模擬裝置的降雨特性均不盡相同，其中針管式降雨裝置降雨強度控制范圍有限；懸線式降雨裝置在國內(nèi)應(yīng)用不多；噴嘴式降雨裝置主要應(yīng)用于室外，且側(cè)噴式裝置因雨滴直徑較大而受風影響大，下噴式的雨滴落地破壞力不易計算，上噴式的降雨均勻度不夠穩(wěn)定．此外，現(xiàn)有降雨模擬裝置研制尚存在雨水難以自循環(huán)利用、試驗周期較長、設(shè)備安裝調(diào)試復(fù)雜、試樣淋雨平臺不便于上人操作等問題．在裝置的降雨特性率定方面，國內(nèi)外對降雨強度、降雨均勻性、雨滴直徑、雨滴終點落地速度與降雨動能等開展了研究．比如，吳光艷[8]等通過對天然降雨雨滴特征的研究驗證了雨滴直徑特征符合 Best分布函數(shù)；美國學者羅斯等人關(guān)于天然降雨雨滴的研究表明，天然降雨雨滴大小的分布，波動在0-6 mm[9]．但是，不同人工模擬降雨裝置由于其產(chǎn)生降雨的機制不同，亦需有針對性地開展降雨特性試驗研究．

因此，本文借鑒國內(nèi)外降雨模擬裝置的經(jīng)驗，研制了一套自循環(huán)式人工降雨模擬裝置，并開展了該裝置在不同雨強條件下的降雨強度、降雨均勻性、雨滴直徑、雨滴終點落地速度與降雨動能等降雨特性率定試驗，為實驗室在短期內(nèi)開展大批量降雨模型試驗和獲取實驗數(shù)據(jù)提供保障．

1 試驗材料與方法

1.1 人工降雨模擬基本原理

降雨裝置雨滴的動能是模擬天然降雨侵蝕的最好參量[10]．當降雨高度不變且雨滴大小均勻，即雨滴終點落地速度v為常數(shù)時，時間t內(nèi)降雨對淋雨平臺單位面積上試驗土體作用的總動能為

式中：E為時間段t內(nèi)單位面積上的降雨對淋雨平臺作用的總動能(J/(m2·s))；m為時間段t內(nèi)單位面積上的降雨總質(zhì)量(kg/m2)；v為雨滴終點落地速度(m/s)；r為雨滴密度(kg/m3)；I為降雨強度(mm/h)；t為降雨持續(xù)時間(s)；k為常數(shù)．

因此，在人工降雨模擬試驗中，保證降雨高度不變，并使雨滴的大小及其分布基本不變，則只需通過調(diào)整裝置的降雨強度就能控制相應(yīng)的降雨動能，即降雨能量相似可通過相應(yīng)的雨強來表達．另外，理想的降雨模擬裝置所產(chǎn)生的人工降雨的主要降雨特性還包括降雨強度、降雨均勻性、雨滴直徑分布、雨滴降落終點速度等，本文主要針對這些特征進行率定實驗，分析自循環(huán)式人工降雨模擬裝置產(chǎn)生的人工降雨與天然降雨的相似性,為該裝置用于科學實驗提供依據(jù)．

1.2 人工降雨模擬系統(tǒng)

研制的自循環(huán)式人工降雨模擬裝置屬于管網(wǎng)式降雨裝置，由降雨系統(tǒng)、自循環(huán)供水系統(tǒng)、支撐系統(tǒng)和試樣淋雨平臺等四大部分組成，如圖1所示．該裝置由潛水泵供水，水流經(jīng)過供水系統(tǒng)輸送到降雨系統(tǒng)，產(chǎn)生穩(wěn)定降雨．通過調(diào)節(jié)導水管上的控水閥門的開度，調(diào)整降雨強度，從而控制相應(yīng)的降雨動能．

圖1 自循環(huán)式降雨模擬裝置

降雨裝置的降雨系統(tǒng)由降雨管網(wǎng)、導水管和擋雨帷幕組成，降雨面積約 8 m2．降雨管網(wǎng)由PVC管拼接而成，固定在支撐系統(tǒng)頂部，每根PVC管鉆設(shè)降雨噴孔；導水管位于降雨管網(wǎng)兩側(cè)．

降雨裝置的自循環(huán)供水系統(tǒng)由簡易蓄水池、潛水泵、連接蓄水池與降雨系統(tǒng)的進水管及安裝于進水管上的控水閥門和水壓表組成．其自循環(huán)原理為潛水泵將蓄水池內(nèi)的水通過進水管輸送至導水管，再由導水管導入降雨管網(wǎng)，降雨管網(wǎng)中的雨水降到管網(wǎng)正下方的試樣淋雨平臺上，然后通過模板上的透水孔重新回到蓄水池，如此循環(huán)，整個過程既環(huán)保又節(jié)水．蓄水池池壁由固定在支撐系統(tǒng)上的模板組成，其防水層由多層防水帆布與防水薄膜組成；進水管采用PVC材料制成．

降雨裝置的支撐系統(tǒng)由鋼管組裝而成，包括立桿、縱向支撐桿、橫向支撐桿和連接扣件．降雨裝置試樣淋雨平臺由帶透水孔的模板鋪裝而成，通過支撐系統(tǒng)固定在蓄水池上方，不僅能在淋雨平臺同時布置多個雨水入滲侵蝕試驗?zāi)Ｐ?，還能上人安裝試樣模型、進行試驗操作及數(shù)據(jù)采集．

1.3 降雨特性參數(shù)率定取樣點

采用等距取樣法在試樣淋雨平臺上確定共35個取樣點，并在采樣點上擺放雨量計，形成陣列，如圖2所示．待降雨達到穩(wěn)定后，依次進行降雨強度、降雨均勻性、雨滴直徑、雨滴終點落地速度與降雨動能等降雨特性試驗．

圖2 降雨模擬裝置率定實驗取樣點分布

2 試驗結(jié)果與分析

2.1 降雨強度及降雨均勻性

降雨強度是指單位時段內(nèi)的降雨量，采用自制簡易雨量計率定降雨強度，每場降雨結(jié)束后，用量筒測量每個雨量計的降雨量．為防止噴孔被雜質(zhì)堵塞而影響降雨的均勻性，率定前先將供水壓力調(diào)整到較大值，把積存在管網(wǎng)內(nèi)的雜質(zhì)沖出，再在每1個取樣點分別放置雨量計，確保雨滴經(jīng)過漏斗進入相應(yīng)雨量計的量筒中．率定時，在控水閥門不同開度下，分別測出各取樣點上雨量計的降雨量Xi，再求出各取樣點的平均降雨量，控水閥門不同開度下降雨裝置的平均降雨強度為

式中：為試樣淋雨平臺上各取樣點的平均降雨強度(mm/h)；為試樣淋雨平臺上各取樣點的實測平均降雨量(mm)．根據(jù)實測平均降雨強度的計算結(jié)果，繪制平均降雨強度—控水閥門開度的擬合曲線(見圖3)，得到降雨裝置降雨強度的線性擬合計算公式

式中：I為降雨裝置降雨強度的擬合推算值(mm/h)；q為控水閥門開度(°)．式(3)可以指導用于實驗的降雨強度設(shè)計，由式(3)計算可得，裝置可模擬的降雨強度范圍為4.65 mm/h-406.05 mm/h．

圖3 控水閥門開度-降雨強度擬合曲線

裝置的降雨均勻性可根據(jù)同一時段內(nèi)降雨面上多個測點的降雨量來計算，降雨均勻度

式中：K為降雨均勻度(%)；Xi為各取樣點上雨量計的實測降雨量(mm)；n為取樣點個數(shù)，取為35．由式(4)計算可得，在不同降雨強度下，裝置的降雨均勻度均超過85 %．限于篇幅，列舉2例，如下：降雨強度為 30 mm/h時，降雨均勻度為90.91 %，其降雨強度等值線見圖4(a)；降雨強度為150 mm/h時，降雨均勻度達到了96.73 %，其降雨強度等值線見圖4(b)．

2.2 雨滴直徑的率定試驗

雨滴直徑的大小是影響降雨落地速度的重要參數(shù)．理論上，一定高度的天然雨滴受自身重力作用降落時，隨著速度的不斷增大，受到的空氣阻力也不斷增大，從而變形破裂，體積不斷變小，因此雨滴直徑一般不會超過 6 mm．目前，率定雨滴直徑的方法主要有：面粉球法、攝影法以及色斑法．本文采用歷史悠久、適用廣泛的濾紙色斑法，利用不同孔徑的注射針頭模擬雨滴，該方法基于水滴在同種濾紙材料上產(chǎn)生的色斑大小與水滴的真實粒徑大小成正比的假定．預(yù)先率定濾紙色斑直徑與雨滴真實粒徑之間的關(guān)系，然后通過量測降雨裝置的雨滴在相同濾紙材料上形成的色斑大小推算相應(yīng)的雨滴粒徑．

圖4 降雨強度等值線

先將針頭裝在注射器上，并向注射器內(nèi)注入一定量的清水，用電子秤稱量含水注射器的質(zhì)量；然后，將經(jīng)過細致研磨的高錳酸鉀粉末與滑石粉按1︰4(質(zhì)量比)比例攪拌均勻，并盡可能薄地涂于濾紙表面；最后，從濾紙正上方一定高度處將注射器中的清水滴到濾紙上，滴滿一張濾紙并注意使水滴大小均勻且色斑不粘在一起，完成后測定注射器(含剩余水)的質(zhì)量．用不同孔徑的針頭重復(fù)上述步驟若干次，得到若干雨滴色斑(代表性色斑見圖5)．

圖5 雨滴色斑

如果一張濾紙上總共有n個雨滴，可按

計算雨滴實際粒徑的平均直徑式中：為雨滴實際粒徑的平均直徑(mm)；為一張濾紙上所有雨滴的質(zhì)量(g)；m為一張濾紙上雨滴的個數(shù)；r為水的密度；π為圓周率．

根據(jù)試驗中率定的水滴直徑與相應(yīng)的色斑直徑的實測和計算結(jié)果，繪制雨滴直徑—色斑直徑擬合曲線圖(見圖6)．

圖6 雨滴直徑—色斑直徑擬合曲線圖

分析得到雨滴實際粒徑的平均直徑和濾紙色斑直徑平均值之間的推算公式

在降雨裝置模擬降雨試驗中，可通過式(6)由雨滴在濾紙上形成的色斑大小推算雨滴實際粒徑．先將海綿平放在雨滴直徑取樣盒底板上，再將涂好粉末的濾紙平放在海綿上，合上蓋板．調(diào)節(jié)控水閥門，在降雨達到穩(wěn)定后，依次在每個取樣點上放置取樣盒，均勻滑動取樣盒蓋板，當雨滴落到濾紙上的瞬間，立即合上蓋板，并取出取樣盒．經(jīng)過若干次降雨采樣后，可以得到若干降雨裝置雨滴色斑譜．通過分析，可得到裝置的雨滴直徑累積頻率分布曲線(見圖 7)，雨滴粒徑在0.94 mm-2.16 mm之間．

圖7 降雨裝置雨滴直徑累積頻率分布曲線

2.3 雨滴終點落地速度與降雨動能的率定試驗

雨滴終點落地速度是研究降雨動能的重要參數(shù)．采用目前應(yīng)用較廣的公式計算雨滴終點落地速度，即當雨滴直徑時，采用修正的沙玉清公式

式中，為雨滴終點落地速度(m/s)．由式(7)或(8)計算可得，裝置的雨滴終點落地速度的范圍在3.69 m/s-7.14 m/s，與天然降雨終點速度(1.58 m/s-9.44 m/s[11-12])有較好地相似性．

得到所有雨滴樣本的終點落地速度，即可計算單個雨滴的降雨動能．因為裝置降雨穩(wěn)定且都是降雨樣本隨機取樣，可認為一次降雨的能量為多組樣本降雨能量的組合．單位時間單位面積內(nèi)一定降雨強度條件下的降雨動能可用

計算．式中：為單位時間單位面積內(nèi)降雨動能；為一次取樣的雨滴降雨能量和；s為取樣紙的面積；t為在雨中取樣時間．由此計算試驗降雨強度條件下的降雨動能，裝置的降雨動能隨降雨強度的變化如圖8所示，并對降雨動能和降雨強度作回歸分析,可以得到降雨強度與降雨動能的關(guān)系式在圖示的降雨強度范圍內(nèi)，降雨動能與降雨強度有著良好的線性關(guān)系，通過控制降雨強度可以實現(xiàn)降雨動能相似，已知降雨強度就可以計算出其各次降雨降雨動能，用以指導降雨實驗分析．

圖8 降雨動能隨降雨強度的變化趨勢

3 結(jié)論

(1)自循環(huán)式人工降雨模擬裝置模擬降雨與天然降雨的相似程度較高，實現(xiàn)了雨水自循環(huán)利用，裝置可控性好、性能穩(wěn)定，可模擬的降雨強度范圍為4.65 mm/h -406.05 mm/h，降雨均勻度超過85 %，雨滴粒徑在0.94 mm-2.16 mm之間，雨滴終點落地速度在3.69 m/s -7.14 m/s之間，可以滿足人工模擬降雨試驗要求，為實驗室在短期內(nèi)開展大批量降雨模型試驗和獲取實驗數(shù)據(jù)提供了技術(shù)保障．

(2)自循環(huán)式人工降雨模擬裝置的降雨強度—控水閥門開度關(guān)系式為降雨強度與降雨動能的關(guān)系降雨動能與降雨強度有良好線性相關(guān)關(guān)系，可通過調(diào)控雨強實現(xiàn)降雨動能與天然降雨的相似．

[1]夏平, 蔣建清, 蔡晶垚, 等. 人工降雨模擬裝置的研制與工程應(yīng)用進展綜述[J]. 企業(yè)技術(shù)開發(fā), 2015, 34(12): 4-6.

[2]徐向舟, 劉大慶, 張紅武, 等. 室內(nèi)人工模擬降雨試驗研究[J].北京林業(yè)大學學報, 2006, 28(5): 52-58.

[3]孫愷, 張季如. 針管式人工降雨裝置的研究與應(yīng)用[J]. 武漢理工大學學報, 2013, 35(12): 125-129.

[4]Tas B B, Burt T P. Rainfall simulators for investigating soil response to rainfall[J]. Soil Technology, 1989, 2(1): 1-16.

[5]霍云梅, 畢華興, 朱永杰, 等. QYJY-503C人工模擬降雨裝置降雨特性試驗[J]. 中國水土保持科學, 2015, 13(2): 31-36.

[6]蘇溦娜, 田一梅, 高波, 等. 人工模擬降雨裝置的設(shè)計及其參數(shù)率定[J]. 水土保持通報, 2015, 35(6): 120-123.

[7]張晶玲, 周麗麗, 黃毅. 搖臂上噴式降雨模擬器降雨特性分析[J]. 中國水土保持科學, 2016, 14(6): 125-130.

[8]吳光艷, 郝民利, 劉超群, 等. 天然降雨與人工降雨特性的研究[J]. 人民珠江, 2013(2): 5-7.

[9]任樹梅, 劉洪祿, 顧濤. 人工模擬降雨技術(shù)研究綜述[J]. 中國農(nóng)村水利水電, 2003(3): 73-75.

[10]Wischmeier W H, Smith D D. Rainfall energy and its relation to soil loss[J]. Transactions American Geophysical Union, 1958,39(2): 285-291.

[11]吳長文, 徐寧娟. 擺噴式人工降雨機的特性試驗[J]. 南昌大學學報, 1995, 17(1): 58-66.

[12]趙志進, 李桂英. 人工模擬降雨機具和方法的發(fā)展研究與展望 [J]. 中國水土保持學報, 1989(5): 30-33.

(責任編校：徐贊)

Experimental Study on Rainfall Characteristics of Self-circulation Artificial Rainfall Simulator

JIANG Jian-qing,CHENG Chao,CAI Jing-yao,XIA Ping,WANG Liang,Xiao Wu-wang
(School of Civil Engineering, Hunan City University, Yiyang, Hunan 413000, China)

In view of the shortage of the existing rainfall simulation device and the characteristics of the natural rainfall, a self-circulation artificial rainfall simulator for soil and water conservation and geotechnical engineering rainfall model test was developed and its rainfall characteristics were tested. The relationship between the rainfall intensity and the opening degree of water control valve, the rainfall kinetic energy and the rainfall intensity, the average diameter of raindrop and the mean diameter of filter paper were calculated and the accumulated frequency distribution curve of rainfall droplet diameter was obtained through the experimental study of rainfall characteristics. The results show that the similarity between the simulated rainfall and the natural rainfall is quite high, and it can simulate the rainfall intensity of 4.65 mm/h -406.05 mm/h. The rainfall uniformity is more than 85 %, the size of raindrops is between 0.94 mm and 2.16 mm, and the falling speed of the rain drops is between 3.69 m/s and 7.14 m/s. The device has realized the self-recycling of rainwater and provided technical guarantee for the laboratory to carry out large-scale rainfall model tests and obtain experimental data in a short time.

rainwater self-recycling; rainfall simulation device; rainfall characteristics; development;rate

TU 41

A

10.3969/j.issn.1672-7304.2017.02.0001

1672–7304(2017)02–0001–05

2017-02-16

國家自然科學基金項目(51308198)；湖南省教育廳科研項目(14B034)

蔣建清(1979-)，男，湖南寧鄉(xiāng)人，副教授，博士，主要從事巖土工程、結(jié)構(gòu)工程研究，E-mail: lh201314@163.com.