劉 潔王 聰魏青松史玉升
(華中科技大學(xué)材料成形及模具技術(shù)國(guó)家重點(diǎn)試驗(yàn)室,湖北武漢 430074)
波動(dòng)水壓參數(shù)對(duì)灌水器水力性能影響試驗(yàn)
劉 潔,王 聰,魏青松,史玉升
(華中科技大學(xué)材料成形及模具技術(shù)國(guó)家重點(diǎn)試驗(yàn)室,湖北武漢 430074)
針對(duì)低壓灌溉系統(tǒng)雖然可以節(jié)省電能但卻容易發(fā)生堵塞的不利情況,將原來(lái)恒定不變的水壓改變?yōu)榘凑找欢l率和振幅波動(dòng)的動(dòng)態(tài)水壓模式,以期增強(qiáng)低壓運(yùn)行時(shí)滴灌灌水器的抗堵塞能力。以基礎(chǔ)水壓、波幅和波動(dòng)周期為變量參數(shù)來(lái)設(shè)計(jì)工作水壓的波動(dòng)模式,通過(guò)正交方式實(shí)施波動(dòng)參數(shù)與流量特性試驗(yàn),分析了3種波動(dòng)參數(shù)對(duì)灌水器流量的影響規(guī)律,結(jié)果表明基礎(chǔ)水壓對(duì)灌水器流量的影響最大。
低壓灌溉;波動(dòng)水壓;滴灌灌水器;灌水器水力性能
灌溉方式選擇將對(duì)滴灌系統(tǒng)的灌溉效率、應(yīng)用成本以及系統(tǒng)的水力性能等產(chǎn)生較大影響,所以選擇合適的工作水壓對(duì)于滴灌系統(tǒng)來(lái)說(shuō)至關(guān)重要,它將直接決定系統(tǒng)的灌溉效果。近年來(lái)低壓灌溉由于對(duì)滴灌系統(tǒng)電力驅(qū)動(dòng)要求較低而得到廣泛關(guān)注,范興科等[1]研究了低壓滴灌條件下提高系統(tǒng)灌水均勻度的途徑;張國(guó)祥等[2]對(duì)滴灌系統(tǒng)滴頭設(shè)計(jì)水頭的取值依據(jù)進(jìn)行了研究;張林等[3]對(duì)低壓滴灌灌水均勻度進(jìn)行試驗(yàn)研究;馬曉鵬等[4]針對(duì)常壓滴灌系統(tǒng)中使用較多的6種國(guó)產(chǎn)單翼迷宮和內(nèi)鑲片式滴灌帶,在低壓條件下進(jìn)行了滴灌帶鋪設(shè)長(zhǎng)度、壓力水頭、地形坡度3個(gè)因素對(duì)滴灌帶灌水均勻系數(shù)影響規(guī)律的試驗(yàn)研究;王威等[5]應(yīng)用流體力學(xué)的一般原理分析了灌水器微小流道中的流體流態(tài)特征,并且利用N-S方程和灌水器的水頭損失描述灌水器微內(nèi)流場(chǎng),建立了微小流道中的流體K-E紊流模型;白丹等[6]對(duì)地下滴灌灌水器的壓力和流量關(guān)系進(jìn)行研究,認(rèn)為壓力和流量是最重要的水力要素之一。
研究者們大多采用Darcy-Weisbach方程和Hazen-Williams經(jīng)驗(yàn)方程計(jì)算管網(wǎng)內(nèi)的水頭損失[7],但是在許多情況下,經(jīng)驗(yàn)方差由于計(jì)算準(zhǔn)確性差和計(jì)算量大的弊端難以滿足工程需求。所以很多學(xué)者對(duì)經(jīng)驗(yàn)方差進(jìn)行優(yōu)化和完善,以提高計(jì)算效率[8-9]。
Duan等[10]研究了低工作水壓和常壓(0~310 kPa)在污水滴灌試驗(yàn)條件下壓力補(bǔ)償式灌水器的補(bǔ)償效果。李永光[11]研究了流態(tài)指數(shù)x對(duì)灌水器壓力補(bǔ)償?shù)囊?guī)律,認(rèn)為當(dāng)x<0.4時(shí)灌水器為壓力補(bǔ)償灌水器,x越接近零則灌水器的補(bǔ)償性能越好。Kang[12]研究了工作水壓對(duì)灌溉均勻度的影響,發(fā)現(xiàn)灌水器的均勻度會(huì)以x=0.5為分界點(diǎn):在0<x<0.5時(shí),工作水壓增加則灌溉均勻度會(huì)出現(xiàn)增加或減少2種情況;當(dāng)x>0.5時(shí),工作水壓的增加將使灌溉均勻度減小。Safi等[13]通過(guò)試驗(yàn)研究了多種類型灌水器在不同工作水壓(最低達(dá)到49 kPa水壓)和鋪設(shè)長(zhǎng)度下的出水均勻度;Li等[14]通過(guò)有限元模擬的方法研究微壓(9.8~49 kPa)條件下迷宮式灌水器內(nèi)的流場(chǎng)狀態(tài),發(fā)現(xiàn)在低壓情況下灌水器內(nèi)流場(chǎng)處于紊流流態(tài)。
低壓灌溉雖然可以在一定程度上節(jié)省電能,但滴灌系統(tǒng)中使用的灌水器因其流道曲折且尺寸微小,在恒定低壓工作條件下容易發(fā)生堵塞現(xiàn)象,從而影響低壓灌溉系統(tǒng)的長(zhǎng)期穩(wěn)定運(yùn)行。在相對(duì)較低工作水壓運(yùn)行基礎(chǔ)上,將原來(lái)恒定不變的水壓改變?yōu)榘凑找欢l率和振幅波動(dòng)的動(dòng)態(tài)水壓模式,利用灌水器內(nèi)流場(chǎng)的連續(xù)波動(dòng)效應(yīng),實(shí)現(xiàn)水流對(duì)固體懸浮物的擾動(dòng)和沖擊,從而增強(qiáng)低壓運(yùn)行時(shí)灌水器的抗堵塞能力。該種思路不同于已有的脈沖式滴灌方式[15-16]。本文研究的波動(dòng)水壓動(dòng)態(tài)模式是通過(guò)自動(dòng)控制技術(shù),實(shí)時(shí)改變和控制加壓水泵的運(yùn)行頻率,使工作水壓值按特定波動(dòng)模式動(dòng)態(tài)變化,時(shí)間頻率可短至秒級(jí),波動(dòng)幅度可低至9.8 kPa水壓左右,屬于一種微觀層面的供壓方式。采用該種波動(dòng)水壓后,宏觀上灌水器的出水形式有異于傳統(tǒng)恒壓狀態(tài),微觀上灌水器內(nèi)流場(chǎng)的流動(dòng)特征發(fā)生了較大的改變。
國(guó)外尚無(wú)學(xué)者涉及波動(dòng)水壓對(duì)灌水器水力性能的研究,只有Assouline等[15]和Elmaloqlou等[16]人提出了脈沖式滴灌方式。而各學(xué)者研究的主要方向是研究具體情況下的灌水器堵塞、低工作水壓滴灌系統(tǒng)和低壓對(duì)灌水器水力性能的影響問(wèn)題。
本文主要針對(duì)波動(dòng)水壓工作模式下滴灌灌水器的流量與壓力特征、灌水均勻度等水力性能的演變及其調(diào)控規(guī)律進(jìn)行深入研究,在微灌綜合水力性能測(cè)試平臺(tái)上對(duì)多種灌水器進(jìn)行波動(dòng)水壓測(cè)試,目的在于揭示清水條件下波動(dòng)水壓對(duì)灌水器水力性能的影響規(guī)律,為該創(chuàng)新構(gòu)思在實(shí)際滴灌工程中獲得應(yīng)用奠定理論與技術(shù)基礎(chǔ)。不同波動(dòng)模式對(duì)灌水器堵塞性能影響的試驗(yàn)研究將在后續(xù)的論文中闡述。
為了研究波動(dòng)水壓對(duì)典型結(jié)構(gòu)灌水器水力性能的影響,選取4種不同灌水器產(chǎn)品:NetafimST16125滴灌帶、NetafimST16150滴灌帶、熱壓折疊式滴灌帶和HWP繞流滴灌帶。各種灌水器的性能參數(shù)見(jiàn)表1。
不同的波動(dòng)模式對(duì)波動(dòng)水壓的影響較大,本文選取4種典型模式(圖1):矩形波形、梯形波形、鋸齒波形和正弦波形。每種波動(dòng)模式的主要影響因素有基礎(chǔ)水壓H、波幅W和波動(dòng)周期T。
在每種波動(dòng)模式下,分別選取H、W和T為變量因素,考慮到波動(dòng)水壓要適合在49 kPa(5 mH2O)低工作水壓下正常運(yùn)行,將每個(gè)變量因素選取如表2所示的3種水平值。依據(jù)L9(33)的正交試驗(yàn)表(表3)來(lái)安排各個(gè)因素,并考慮基礎(chǔ)水壓與波幅的交互作用因素。
試驗(yàn)步驟如下:
a.根據(jù)SL/T 67.1—1994《水利技術(shù)標(biāo)準(zhǔn)匯編》灌溉排水卷節(jié)水設(shè)備與材料-微灌灌水器 滴頭,本文在壓力-流量試驗(yàn)中,選取15個(gè)試驗(yàn)樣安裝在微灌綜合水力性能試驗(yàn)平臺(tái)上,在恒壓條件下測(cè)量15個(gè)試驗(yàn)樣在9.8~98 kPa(1~10 mH2O)每一個(gè)壓力點(diǎn)的出水量,試驗(yàn)時(shí)間不得少于2 min。記錄室溫、水溫、水壓、試驗(yàn)時(shí)間和灌水器出水流量。重復(fù)上述試驗(yàn),保證每個(gè)試驗(yàn)樣2次測(cè)得水量之差不大于2%,取平均值,并計(jì)算流量。按照上述方法,分別測(cè)試上述4種灌水器產(chǎn)品在各工作水壓下的流量,并記錄試驗(yàn)相關(guān)數(shù)據(jù),計(jì)算獲得4種灌水器的壓力-流量特性參數(shù)。
b.在上述壓力-流量試驗(yàn)測(cè)試標(biāo)準(zhǔn)下,按照正交試驗(yàn)方案(表3)進(jìn)行壓力-流量測(cè)試試驗(yàn),記錄試驗(yàn)相關(guān)數(shù)據(jù),計(jì)算獲得4種灌水器的壓力-流量特性參數(shù)。
按照上述步驟a得到9.8~98kPa(1~10mH2O)恒壓條件下4種灌水器的試驗(yàn)數(shù)據(jù),并根據(jù)試驗(yàn)記錄數(shù)據(jù)計(jì)算獲得4種灌水器的壓力與流量試驗(yàn)結(jié)果。
按照上述步驟b得到4種灌水器在4種不同波動(dòng)模式下的試驗(yàn)數(shù)據(jù),并根據(jù)試驗(yàn)記錄數(shù)據(jù)計(jì)算獲得4種灌水器在正交試驗(yàn)條件下的壓力與流量試驗(yàn)結(jié)果。正交表中各因素水平的組合參照表3。
對(duì)上述正交試驗(yàn)數(shù)據(jù)進(jìn)行直觀分析和方差分析,得到鋸齒波動(dòng)模式下4種灌水器的波動(dòng)因素對(duì)流量的影響關(guān)系(表4)。
由表4可以看出,對(duì)于灌水器A,在數(shù)據(jù)分析中的4個(gè)因素——H、W、H∧W和T對(duì)應(yīng)的極差和均方差大小順序?yàn)镠>T>W(wǎng)>(H∧W),該結(jié)果表明:對(duì)灌水器流量影響最大的是基礎(chǔ)水壓,其次為波動(dòng)周期,再次為波幅,影響最弱的是基礎(chǔ)水壓與波幅的交互作用。其最優(yōu)方案為H3W2T1。
對(duì)于灌水器B,各因素對(duì)灌水器流量的極差影響主次順序?yàn)?H>W(wǎng)>(H∧W)>T,即基礎(chǔ)水壓的影響最明顯,波動(dòng)周期影響最小。但是,根據(jù)方差分析可以看出,4個(gè)因素的均方差值在一個(gè)數(shù)量級(jí),說(shuō)明這4個(gè)因素對(duì)灌水器B的流量影響作用比較均衡。這是因?yàn)楣嗨鰾為壓力補(bǔ)償式灌水器,鋸齒波波動(dòng)水壓產(chǎn)生的波動(dòng)效應(yīng)在灌水器內(nèi)部被墊片所吸收。其最優(yōu)方案為H3W2T1。
對(duì)于灌水器C,各因素對(duì)灌水器流量的極差影響主次順序?yàn)?H>W(wǎng)>T>(H∧W),即基礎(chǔ)水壓對(duì)灌水器的流量影響最顯著,而基礎(chǔ)水壓與波幅的互交因素影響最小。方差分析表明,W、(H∧W)、T與H的影響程度相差近3個(gè)數(shù)量級(jí)以上,表明除基礎(chǔ)水壓H以外,其他3個(gè)因素可以歸入誤差。其原因?yàn)?鋸齒形波動(dòng)模式下水壓變化幅度大,且灌水器C的流道尺寸較長(zhǎng),波動(dòng)水壓產(chǎn)生的波動(dòng)效應(yīng)在較長(zhǎng)的流道內(nèi)通過(guò)撞擊和摩擦被弱化。其最優(yōu)方案為H3W1T1。
對(duì)于灌水器D,各因素對(duì)灌水器流量的極差影響主次關(guān)系為:H>W(wǎng)>T>(H∧W),即基礎(chǔ)水壓對(duì)灌水器流量的影響最顯著,由方差分析得到(H∧W)與其他3個(gè)因素相差近3個(gè)數(shù)量級(jí),可以被歸入誤差。其最優(yōu)方案為H3W3T2。
對(duì)正交試驗(yàn)數(shù)據(jù)進(jìn)行直觀分析和方差分析,得到正弦波動(dòng)模式下4種灌水器的波動(dòng)因素對(duì)流量的影響關(guān)系如表5所示。
由表5可以看出,對(duì)于灌水器A,各因素對(duì)灌水器流量的極差影響主次關(guān)系為:H>(H∧W)>W(wǎng)>T,即H對(duì)灌水器流量的影響最顯著,而T對(duì)灌水器流量影響最小。其最優(yōu)方案為H3W3T1。
對(duì)于灌水器B,各因素對(duì)灌水器流量的極差影響主次關(guān)系為:H>(H∧W)>W(wǎng)>T,即H對(duì)灌水器流量的影響最顯著,但是由均方差值看出各因素的影響程度較均衡,即各因素對(duì)灌水器流量都有較大的影響作用。其最優(yōu)方案為H3W3T1。
對(duì)于灌水器C,各因素對(duì)灌水器流量的極差影響主次關(guān)系為:H>(H∧W)>T>W(wǎng),即H對(duì)灌水器流量的影響最顯著,而W的影響最小,且由均方差值看出,W、(H∧W)、T均與H相差2個(gè)數(shù)量級(jí)。其最優(yōu)方案為H3W3T1。
對(duì)于灌水器D,各因素對(duì)灌水器流量的極差影響主次關(guān)系為:H>(H∧W)>W(wǎng)>T,即H對(duì)灌水器流量的影響最顯著,而T的影響最小。其最優(yōu)方案為H3W1T2。
對(duì)正交試驗(yàn)數(shù)據(jù)進(jìn)行直觀分析和方差分析,得到梯形波動(dòng)模式下4種灌水器的波動(dòng)因素對(duì)流量的影響的關(guān)系如表6所示。
由表6可以看出,對(duì)于灌水器A,各因素對(duì)灌水器流量的極差影響主次關(guān)系為:H>(H∧W)>W(wǎng)>T,即H對(duì)灌水器流量的影響最顯著,而T的影響最小,其中T與其他3個(gè)因素相差2個(gè)數(shù)量級(jí),所以可以歸入誤差。其最優(yōu)方案為H3W3T1。
對(duì)于灌水器B,各因素對(duì)灌水器流量的極差影響主次關(guān)系為:H>(H∧W)>W(wǎng)>T,即H對(duì)灌水器流量的影響最顯著,而T的影響最小,從均差值看出各因素對(duì)灌水器流量的影響較均衡。其最優(yōu)方案為H2W3T1。
對(duì)于灌水器C,各因素對(duì)灌水器流量的極差影響主次關(guān)系為:H>T>(H∧W)>W(wǎng),即H對(duì)灌水器流量的影響最顯著,而W的影響最小,這是因?yàn)楣嗨鳛槿切瘟鞯?流道較長(zhǎng),由波幅產(chǎn)生的波動(dòng)效應(yīng)被流道削弱。其最優(yōu)方案為H3W1T2。
對(duì)于灌水器D,各因素對(duì)灌水器流量的極差影響主次關(guān)系為:H>(H∧W)>W(wǎng)>T,即H對(duì)灌水器流量的影響最顯著,T的影響最小。其最優(yōu)方案為H3W2T1。
對(duì)正交試驗(yàn)數(shù)據(jù)進(jìn)行直觀分析和方差分析,得到矩形波動(dòng)模式下4種灌水器的波動(dòng)因素對(duì)流量的影響關(guān)系如表7所示。
由表7可以看出,對(duì)于灌水器A,各因素對(duì)灌水器流量的極差影響主次關(guān)系為:H>(H∧W)>W(wǎng)>T,即H對(duì)灌水器流量的影響最顯著,而T的影響最小。其最優(yōu)方案為H2W2T1。
對(duì)于灌水器B,各因素對(duì)灌水器流量的極差影響主次關(guān)系為:H>(H∧W)>W(wǎng)>T,即H對(duì)灌水器流量的影響最顯著,而T的影響最小。其最優(yōu)方案為H3W2T1。
對(duì)于灌水器C,各因素對(duì)灌水器流量的極差影響主次關(guān)系為:H>(H∧W)>T>W(wǎng),即H對(duì)灌水器流量的影響最顯著,而W的影響最小,且從均方差看出W、H∧W和T與H相差2個(gè)數(shù)量級(jí)。其最優(yōu)方案為H3W2T1。
對(duì)于灌水器D,各因素對(duì)灌水器流量的極差影響主次關(guān)系為:H>(H∧W)>W(wǎng)>T,即H對(duì)灌水器流量的影響最顯著,而T的影響最小,且從均方差看出T與其他3個(gè)因素相差2個(gè)數(shù)量級(jí),所以可以歸入誤差。其最優(yōu)方案為H3W2T1。
清水條件下的波動(dòng)因素——基礎(chǔ)水壓H、波幅W和波動(dòng)周期T對(duì)灌水器流量影響的正交試驗(yàn)結(jié)果表明,波動(dòng)因素影響流量的主次關(guān)系為:基礎(chǔ)水壓H>波幅W>波動(dòng)周期T,基礎(chǔ)水壓對(duì)灌水器流量影響最大,影響率超過(guò)50%。
[1]范興科,吳普特,牛文全,等.低壓滴灌條件下提高系統(tǒng)灌水均勻度的途徑探討[J].灌溉排水學(xué)報(bào),2008,27(1):18-20.(FAN Xingke,WU Pute,NIU Wenquan,et al.The methods of improving system's irrigation uniformity under low-pressure drip irrigation[J].Journal of Irrigation and Drainage,2008,27(1):18-20.(in Chinese))
[2]張國(guó)祥,吳普特.滴灌系統(tǒng)滴頭設(shè)計(jì)水頭的取值依據(jù)[J].農(nóng)業(yè)工程學(xué)報(bào),2005,21(9):20-22.(ZHANG Guoxiang,WU Pute.Determination of the design working head of emitter[J].Transactions of the Chinese Society of Agricultural Engineering, 2005,21(9):20-22.(in Chinese))
[3]張林,吳特普,范興科,等.低壓滴灌灌水均勻度試驗(yàn)研究[J].西北農(nóng)林科技大學(xué)學(xué)報(bào):自然科學(xué)版,2009,37(12):215-220.(ZHANG Lin,WU Pute,FAN Xingke,et al.Experimental research on drip irrigation uniformity at low-pressure[J].Journal of Northwest A&F University:Natural Science Edition,2009,37(12):215-220.(in Chinese))
[4]馬曉鵬,龔時(shí)宏,王建東,等.低壓條件下滴灌帶灌水均勻系數(shù)試驗(yàn)研究[J].灌溉排水學(xué)報(bào),2010,29(4):8-12.(MA Xiaopeng,GONG Shihong,WANG Jiandong,et al.Uniformity coefficient of drip irrigation tapes under low operating pressure [J].Journal of Irrigation and Drainage,2010,29(4):8-12.(in Chinese))
[5]王威,翟治芬.灌水器水力性能參數(shù)的計(jì)算方法[J].安徽農(nóng)業(yè)科學(xué),2008,36(26):11618-11620.(WANG Wei,ZHAI Zhifen.Calculation method for the hydraulic performance parameter of emitter[J].Journal of Anhui Agricultural Sciences,2008, 36(26):11618-11620.(in Chinese))
[6]白丹,宋立勛,王曉愚,等.地下滴灌灌水器流量和壓力關(guān)系的試驗(yàn)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2008,39(8):195-197.(BAI Dan, SONG Lixun,WANG Xiaoyu,et al.Flow rate and pressure relationships tests of subsurface drip irrigation[J].Journal of Agricultural Machinery,2008,39(8):195-197.(in Chinese))
[7]WATTERS G Z,KELLER J.Trickle irrigation tubing hydraulics[M].St Joseph,Michigan:ASAE,1978.
[8]JUANA L,RODRIGUEZ-SINOBAS L,LOSADA A.Determining minor head losses in drip irrigation laterals I:methodology [J].Journal of Irrigation and Drainage Engineering,2002,128(6):376-384.
[9]PROVENZANO G,PUMO D,DIO P D.Simplified procedure to evaluate head losses in drip irrigation laterals[J].Journal of Irrigation and Drainage Engineering-ASCE,2005,131(6):525-532.
[10]DUAN X J,LESIKAR B J,KENIMER A L,et al.Performance of wastewater subsurface drip emitters at low and normal pressure[J].Water Environment Research,2008,80(2):142-148.
[11]李光永.壓力補(bǔ)償灌水器的研究新進(jìn)展與使用中應(yīng)注意的幾個(gè)問(wèn)題[J].節(jié)水灌溉,2000(3):19-21.(LI Guangyong.New advances and several using issues should pay attention to the pressure compensating emitters[J].Water Saving Irrigation,2000 (3):19-21.(in Chinese))
[12]KANG Y.Effect of operating pressures on microirrigation uniformity[J].Irrig Sci,2000,20:23-27.
[13]SAFI B,NEYSHABOURI M R,NAZEMI A H,et al.Water application uniformity of a subsurface drip irrigation system at various operating pressures and tape lengths[J].Turkish Journal of Agriculture and Forestry,2007,31(5):275-285.
[14]LI Yunkai,YANG Peiling,XU Tingwu,et al.Hydraulic property and flow characteristics of three labyrinth flow paths of drip irrigation emitters under micro-pressure[J].Transactions o the ASABE,2009,52(4):1129-1138.
[15]ASSOULINE S,MOLLER M,COHEN S,et al.Soil-plant system response to pulsed drip irrigation and salinity:bell pepper case study[J].Soil Science Society of America Journal,2006,70(5):1556-1568.
[16]ELMALOQLOU S,DIAMANTOPOUOS E.Wetting front advance patterns and water losses by deep percolation under the root zone as influenced by pulsed drip irrigation[J].Agricultural Water Management,2007,90(1/2): 160-163.
Experimental study of effect of fluctuating water pressure factors on hydraulic properties of drip emitters
LIU Jie,WANG Cong,WEI Qingsong,SHI Yusheng
(State Key Laboratory of Materials Processing and Die and Mould Technology, Huazhong University of Science and Technology,Wuhan 430074,China)
Although a low-pressure irrigation system can save power energy,it is prone to clogging.In this study, the constant water pressure was changed to dynamic pressure fluctuating with a certain frequency and amplitude,in order to increase the drip emitter's anti-clogging capability during the runtime at a low pressure.The dynamic model of operating pressure was designed by taking the base pressure,fluctuation amplitude,and fluctuation cycle as variable parameters.With orthogonal experiments,the fluctuation parameters and flow properties were tested.The effects of three fluctuation parameters on the emitter flow were analyzed.The study results show that the water base pressure has the most significant effect on the emitter flow.
low-pressure irrigation;fluctuating water pressure;drip emitter;hydraulic properties of drip emitter
S277
:A
:1000-1980(2014)04-0361-06
10.3876/j.issn.1000-1980.2014.04.015
2013-06 07
國(guó)家自然科學(xué)基金(51109087)
劉潔(1972—),女,廣西南寧人,副教授,博士,主要從事節(jié)水產(chǎn)品快速制造研究。E-mail:hustliuj@mail.hust.edu.cn
河海大學(xué)學(xué)報(bào)(自然科學(xué)版)2014年4期