馬萍,宋鋒惠,史彥江,王健,朱首軍?
(1.西北農(nóng)林科技大學(xué)資環(huán)學(xué)院,712100,陜西楊凌;2.新疆林業(yè)科學(xué)研究院治沙所,830000,烏魯木齊)
伊犁河流域裸露緩坡地放水試驗研究
馬萍1,宋鋒惠2,史彥江2,王健1,朱首軍1?
(1.西北農(nóng)林科技大學(xué)資環(huán)學(xué)院,712100,陜西楊凌;2.新疆林業(yè)科學(xué)研究院治沙所,830000,烏魯木齊)
通過野外徑流小區(qū)實地放水沖刷試驗,研究伊犁河流域裸露緩坡地在不同放水條件下土壤水分變化、徑流泥沙特征,并運用水力學(xué)理論分析坡面徑流水力學(xué)參數(shù)特征。結(jié)果表明:1)土壤入滲率隨著放水流量的增加呈現(xiàn)遞增的趨勢,隨著放水量的增加,土壤表層含水量增大,入滲率隨之減少,坡下部土壤水分變化量高于坡上;2)隨著放水流量的增加,輸沙率呈現(xiàn)遞增的趨勢,含沙量呈先增后減的趨勢,含沙量與放水流量存在較好的冪函數(shù)關(guān)系,隨著沖刷時間的延長,含沙量、輸沙率均呈現(xiàn)遞減的趨勢;3)隨著放水流量的增加,徑流平均流速、雷諾數(shù)、Darcy-Weisbach阻力系數(shù)逐漸增大,弗勞德數(shù)逐漸減小,流速與流量、弗勞德數(shù)與雷諾數(shù)均呈良好的冪函數(shù)關(guān)系,水流為層流、急流狀態(tài)。
裸露緩坡地;放水沖刷;土壤水分;徑流泥沙;水力學(xué)參數(shù)
伊犁河流域的水土資源十分豐富,是全疆最具 有開發(fā)潛力的區(qū)域,在新疆經(jīng)濟發(fā)展和西部大開發(fā)中居于重要的戰(zhàn)略地位[1];但流域內(nèi)由于人為不合理的灌溉方式造成嚴(yán)重的水土流失,尤其是流域內(nèi)的新墾區(qū),土層厚薄不均勻,自然坡降較大,若延續(xù)采用傳統(tǒng)的灌溉方式,勢必會在水土開發(fā)中產(chǎn)生水土流失問題[2]。
土壤侵蝕是水流和土壤相互作用的復(fù)雜物理過程[3],徑流和土壤是水土流失的 2個基本因子[4]。據(jù)多年研究[5-6],坡面徑流是造成水土流失的主導(dǎo)因子,同時坡面徑流的沖刷力是土壤侵蝕的主要動力。對坡面徑流水力學(xué)特性的研究一直受到國內(nèi)外學(xué)者的普遍重視,并取得了一些研究成果[7-13]。由于坡面徑流水力學(xué)特性的復(fù)雜性和試驗技術(shù)的限制,對其理論分析、野外觀測和試驗研究都存在一定困難,在沒有成熟坡面徑流理論的情況下,仍在借鑒水力學(xué)、河流動力學(xué)的原理和方法[14-18]。探索坡面徑流水力學(xué)特性有助于從動力學(xué)角度認識坡面水蝕過程及機制,也為坡面水蝕過程物理模型的建立奠定基礎(chǔ),具有重要的理論和實踐意義。
對伊犁河流域水土流失的現(xiàn)狀調(diào)查顯示,除突發(fā)性暴雨型和融雪型水土流失外,伊犁河流域最普遍、最重要的水土流失來源于該區(qū)不合理的灌溉方式——大水漫灌所造成的水土流失;因此,通過野外徑流小區(qū)實地放水沖刷試驗,深入研究不同放水量、不同放水流量條件下土壤水分變化、徑流泥沙特征及坡面徑流水力學(xué)特性,從而揭示該區(qū)坡面流侵蝕動力機制,為伊犁河流域新墾區(qū)水土資源開發(fā)提供理論支撐和實踐指導(dǎo)。
試驗地位于新疆伊犁哈薩克自治州林業(yè)科學(xué)研究所科技示范基地,地理坐標(biāo)為 E 80°31′~81°43′、N 43°17′~43°57′,海拔約 690m,屬于山前洪積 -沖積傾斜平原。氣候?qū)俅箨懶员睖貛睾透珊禋夂?具有熱量豐富、光照充足、四季分明、夏季炎熱、冬季寒冷,晝夜溫差較大等特點。年降水量 150~250mm,年平均氣溫7.9℃,野生植被以蒿類多年生荒漠植物為主。
伊犁河流域的農(nóng)業(yè)耕作土地大多為近年來大面積新開墾的土地,土壤為灰鈣土,密度為 1.29 g/cm3,有機質(zhì)質(zhì)量分數(shù)為 2.37%,土壤 pH值為8.02,含鹽量 1.18%。按照當(dāng)?shù)剞r(nóng)耕地實際情況,徑流小區(qū)布設(shè)的坡向、坡度與當(dāng)?shù)剞r(nóng)業(yè)耕作土地一致,坡長依實際坡面長度而定。根據(jù)坡面實際面積布設(shè) 9個試驗小區(qū),每個小區(qū)面積均 42m2(21m×2m),南北坡向,坡度 2°。小區(qū)內(nèi)無雜草,為裸地小區(qū)。
試驗裝置由供水系統(tǒng)、試驗沖刷區(qū)、集流系統(tǒng)組成。坡面上方的供水系統(tǒng)由容積為 9m3的水箱、流量調(diào)節(jié)閥門和出水管等幾部分組成。試驗開始時用水泵將水抽到水箱里,水箱設(shè)有溢流孔,使水箱水位水壓恒定,從而保證流量穩(wěn)定。流量按照試驗設(shè)計流量通過閥門進行調(diào)節(jié),水箱里的水通過出水管供給坡面穩(wěn)流槽,穩(wěn)定后的水流流入試驗沖刷區(qū),小區(qū)下方設(shè)有集流桶,用來收集徑流泥沙。試驗裝置如圖 1所示。
圖 1 試驗裝置示意圖Fig.1 Experiment system mode
本試驗?zāi)M當(dāng)?shù)剞r(nóng)業(yè)灌溉制度,設(shè)計 2.25、3.38、4.50m33個放水量,0.45、1.26和 1.80 L/s 3個放水流量,重復(fù) 2次,共進行 18組試驗。放水試驗前用噴壺對小區(qū)進行均勻灑水,使各小區(qū)放水前期土壤含水量(質(zhì)量含水量)基本一致,均為 15%。每次試驗前都對表土進行翻耕,使地表平整,放水后將坡面沖刷出的細溝用耙子耙平,保持各次試驗地表狀況一致。
1.2.1 土壤含水量 放水試驗前后,在小區(qū)坡面上、中、下各設(shè) 3個斷面,在每個斷面上分別取 3個水分觀測點,取其平均值作為該斷面土壤含水量。用土鉆取土,取土深度為 0~20 cm,采用烘干法測定,取土后將取土坑填平。
1.2.2 坡面徑流流速 待坡面產(chǎn)流穩(wěn)定后,用染色法測定坡面徑流流速。將坡面從坡頂至坡底每隔 5 m用紅油漆在小區(qū)左右側(cè)圍梗處做標(biāo)記,待坡面產(chǎn)流穩(wěn)定后,用染色法分別測定 0~5、5~10、10~15、15~20m各斷面的徑流流速,最后取平均值作為小區(qū)全坡面的流速。考慮到用染色法測定的徑流流速為坡面優(yōu)勢流流速,實測流速乘以修正系數(shù) 0.75作為坡面水流的平均流速[19]。
1.2.3 徑流泥沙 放水結(jié)束后將集流桶中收集的渾水充分攪勻,立即采取水樣 1 000m L,在室內(nèi)經(jīng)靜置,過濾之后,采用烘干法測定水樣中的泥沙含量,采用體積法確定泥水總量。
根據(jù)水分平衡公式可知,放水量為徑流量、入滲量與蒸發(fā)量之和,由于水分蒸發(fā)量與徑流量和入滲量相比很小,故忽略水分蒸發(fā)量。土壤入滲率 i為單位時間單位面積上土壤水分下滲量,單位為 mm/min,由上述分析可知,土壤入滲率
式中:W為放水量,m3;R為徑流量,m3;S為小區(qū)面積,m2;t為放水時間,min;1000為單位換算系數(shù)。
徑流深 h是反映水力學(xué)特征的重要因子。由于坡面水流水層極薄,且土壤下墊面條件不斷發(fā)生變化,采用實測法難以準(zhǔn)確測定;因此,假定水流沿坡面均勻分布,用下式[20]計算。
式中:q為單寬流量,m3/(m?min);v為斷面平均流速,m/min;B為過水?dāng)嗝鎸挾?m。
雷諾數(shù)Re是水流型態(tài)的重要判據(jù),其計算公式為
式中:νm為含沙水流的運動黏滯系數(shù),m2/s;由于本試驗區(qū)含沙水流運動黏滯系數(shù)尚無研究,若按清水進行計算勢必會產(chǎn)生大的誤差,故此處借鑒沙玉清[21]公式進行計算。
式中:ν0為同溫清水的運動黏滯系數(shù),ν0=1.0×10-6m2/s;Sv為體積含沙率,%;d50為土壤中值粒徑,mm,取 0.075 mm,數(shù)據(jù)來源于試驗地察布查爾縣土壤調(diào)查。
弗勞德數(shù) Fr作為判別明渠水流流態(tài)的標(biāo)準(zhǔn),其公式為
式中 g為重力加速度,m/s2。
Darcy-Weisbach阻力系數(shù) f主要由顆粒阻力、形態(tài)阻力、波阻力和降雨阻力 4部分組成,最終反映坡面薄層水流受阻情況。
式中:J為水流坡降,可近似取為 sinθ,θ為小區(qū)坡度。
2.1 不同放水處理對坡面土壤水分入滲的影響
試驗設(shè)計通過坡面放水來完成徑流沖刷過程,設(shè)計流量較大,與當(dāng)?shù)卮笏喾绞奖3忠恢?可引起超滲產(chǎn)流。試驗結(jié)果見表 1,可看出,在各放水處理下土壤入滲率在 0.372~1.847 mm/min之間變化,平均值為 0.999mm/min。放水流量為 1.80 L/s時,平均入滲率最大,為 1.591mm/min,比放水流量為 1.26 L/s的高 0.578mm/min。放水流量為 0.45 L/s時,其平均入滲率最低。入滲率隨著放水流量的增加呈現(xiàn)遞增的趨勢。其原因主要是:放水強度超過土壤入滲能力后產(chǎn)生超滲產(chǎn)流,土壤在蓄水未達到飽和之前,放水流量增加,入滲率隨之增加;隨著放水量的增加,土壤表層含水量增大,其入滲率隨之減少。
表 1 不同放水處理坡面土壤水分入滲特征Tab.1 Characteristics of soil infiltration under different treatments
在同一放水量及前期土壤含水量相當(dāng)?shù)臈l件下,同一坡面不同坡位放水試驗前、后土壤水分變化量存在一定差異。圖 2顯示了 3種放水流量處理下各小區(qū)放水試驗前、后表層(0~20 cm)土壤土壤水分變化量,可見相同放水流量下,坡面下部土壤水分變化量高于坡上。這是因為,對于裸地,放水流量達到一定程度,坡面水分來不及向下滲透,產(chǎn)生超滲徑流;但相對坡面上部而言,坡面下部受水體重力作用影響,水分下滲的總量高于坡上。
圖 2 不同坡位試驗前、后土壤水分變化量Fig.2 Changes of soil water content before and after experiment on slope
表 2顯示了不同放水處理坡面徑流、泥沙的特征,可以看出,隨著放水流量的增加,輸沙率呈現(xiàn)遞增的趨勢,放水流量 1.26 L/s是一個轉(zhuǎn)折點。當(dāng)流量小于 1.26 L/s時,輸沙率隨流量增加增長迅速;當(dāng)流量大于 1.26 L/s時,輸沙率隨流量增加增長緩慢。輸沙率增長是因為當(dāng)流量由 0.45 L/s增加到1.80 L/s時,徑流流速增大,因而徑流具有的動能也增大,其攜沙能力增強,沿細溝流走的徑流下蝕和側(cè)蝕作用加強,大量泥沙被攜走,所以輸沙率會呈現(xiàn)增加趨勢。
表2 不同放水處理坡面徑流和泥沙的特征Tab.2 Characteristics of runoff and sediment under different treatments
根據(jù)不同放水量條件下含沙量和放水流量數(shù)據(jù),繪出它們的散點圖(圖 3),發(fā)現(xiàn)含沙量與放水流量之間存在明顯的關(guān)系。通過回歸模擬,得出含沙量與放水流量的冪函數(shù)關(guān)系式為
式中:Ws為含沙量,g/L;Q為放水流量,L/s。
圖3 含沙量和放水流量的關(guān)系Fig.3 Relationship between the sediment concentration and the scouring water discharge
不同放水量沖刷時間不同,由表 2可以看出,隨著沖刷時間的延長,含沙量、輸沙率均呈現(xiàn)遞減的趨勢。主要是因為:在放水過程中會有細溝產(chǎn)生,當(dāng)細溝發(fā)育穩(wěn)定時,水流會沿溝道流走。隨著沖刷時間的延長,沿穩(wěn)定細溝流走的徑流增加,徑流系數(shù)由0.354增大至 0.425。表層土壤在水的浸泡下抗沖性會降低,但隨著沖刷時間的延長,抗沖性降低至一定程度便會相對穩(wěn)定。沿細溝流走的徑流不斷增加而土壤抗沖性相對穩(wěn)定,所以含沙量、輸沙率呈降低趨勢。
在坡面侵蝕過程中,水流水力學(xué)參數(shù)及其流態(tài)變化將對侵蝕過程產(chǎn)生重要影響。徑流流速、徑流深、雷諾數(shù)、弗勞德數(shù)及阻力系數(shù)等水力要素是反映水流動力學(xué)特征的主要參數(shù)[22]。通過了解坡面徑流的水力學(xué)參數(shù)變化特征可以從一定程度上了解坡面土壤侵蝕狀況。坡面徑流水力學(xué)參數(shù)試驗與計算結(jié)果見表 3。
表 3 不同放水處理坡面流水力學(xué)特性試驗結(jié)果Tab. Hydraulic characteristics of surface runoffunder different scouring water discharge
2.3.1 徑流平均流速 根據(jù)坡面水流平均流速和放水流量數(shù)據(jù),繪出它們的關(guān)系曲線(圖 4),可以看出,徑流平均流速隨流量的增加而增大。對曲線進行擬合可知,徑流平均流速與放水流量可表示為
圖4 徑流平均流速與放水流量的關(guān)系Fig.4 Relationship between the average velocity and the scouring water discharge
這一結(jié)果表明對于本試驗而言,坡面徑流的平均流速可以用放水流量的簡單冪函數(shù)關(guān)系式來表示,這與 G.Govers[23](v=3.52Q0.294)、M.Nearing等[24](v=9.802Q0.459)的試驗結(jié)果基本一致。
2.3.2 坡面水流流態(tài) 雷諾數(shù) Re是水流運動狀況(層流和紊流)的重要判別依據(jù),反映了水流的慣性力與粘滯力之比,量綱為 1。根據(jù)明渠均勻流的基本理論,對于矩形斷面的明渠水流層流和紊流的臨界雷諾數(shù)為 500,在雷諾數(shù) 500左右則屬于過渡流。而本試驗各放水處理下的雷諾數(shù)為 102.822~436.102,流態(tài)均屬于層流的范疇。不同處理下雷諾數(shù)隨著放水流量的增大均呈現(xiàn)增大的趨勢,反映出水力侵蝕能力和搬運能力的增大。
在水力學(xué)研究中除了雷諾數(shù)之外,弗勞德數(shù) Fr是表征水流流態(tài)的又一重要參數(shù),它綜合反映了流速和水深的對比關(guān)系。在徑流流量相同的條件下,弗勞德數(shù)越大,說明坡面徑流的流速越大,徑流挾沙能力越強,水深越淺,坡面的徑流侵蝕力越小。根據(jù)明渠水流的判別標(biāo)準(zhǔn),本試驗各放水處理的弗勞德數(shù)均大于 1,均屬于急流狀態(tài),且隨著放水流量的增大而逐漸減小。
現(xiàn)將各放水流量下坡面水流的雷諾數(shù)與弗勞德數(shù)的值點繪于圖 5,可知不同放水流量條件下,坡面水流的弗勞德數(shù)與雷諾數(shù)之間呈較好的冪函數(shù)關(guān)系,對圖中的曲線進行擬和可得
Fr=14.654Re-0.2914R2=0.999
2.3.3 阻力系數(shù)變化特征 Darcy-Weisbach阻力系數(shù)反映了坡面流在流動過程中所受的阻力大小[25]。在流量條件相同的情況下,阻力系數(shù)越大,水流克服坡面阻力所消耗的能量就越大,則用于坡面侵蝕和泥沙輸移的能量就越小,土壤侵蝕就越微弱,反之則土壤侵蝕劇烈。從表 3中可以看出,阻力系數(shù)隨著放水流量的增加而增大,其原因為:流量增加意味著徑流深增大,使相對糙率變小;然而流量的增加促使侵蝕強度增強,對細溝微地貌的塑造更加復(fù)雜化。對于本試驗而言,流量增加對微地貌變化的影響較糙率的影響程度更大,因此,阻力系數(shù)隨流量的增加而增大。
圖5 弗勞德數(shù)與雷諾數(shù)的關(guān)系Fig.5 Relationship between the Froude number and the Reynolds number
1)土壤入滲率隨著放水流量的增加呈現(xiàn)遞增的趨勢,隨著放水量的增加,土壤表層含水量增大,入滲率隨之減少。試驗前、后不同坡位土壤水分變化量存在一定差異,坡下土壤水分變化量高于坡上。
2)隨著放水流量的增加,輸沙率呈現(xiàn)遞增的趨勢,含沙量呈先增后減的趨勢,含沙量與放水流量存在較好的冪函數(shù)關(guān)系:Ws=25.773Q1.1441,決定系數(shù)R2=0.884 7。隨著沖刷時間的延長,含沙量、輸沙率均呈現(xiàn)遞減的趨勢。
3)隨著放水流量的增加,徑流平均流速、雷諾數(shù)、Darcy-Weisbach阻力系數(shù)逐漸增大,弗勞德數(shù)逐漸減小。徑流平均流速與放水流量、弗勞德數(shù)與雷諾數(shù)之間均存在良好的冪函數(shù)關(guān)系,決定系數(shù) R2=0.999。坡面水流在放水流量變化范圍內(nèi)為層流、急流狀態(tài)。
以上結(jié)論是在下墊面為裸地的基礎(chǔ)上得到的,要得到更為一般性的研究結(jié)論,尚需開展一系列不同下墊面條件下的相關(guān)試驗研究。
[1]陳亞寧,賀西安,張小云.流域開發(fā)與土地管理模式:兼談伊犁河流域土地開發(fā)與產(chǎn)業(yè)定位.干旱區(qū)地理,2007,30(4):595-600
[2]楊龍,賀光華.淺談伊犁灌區(qū)建設(shè)對生態(tài)環(huán)境的影響.節(jié)水灌溉,2008(3):34-35
[3]吳普特.動力水蝕實驗研究.西安:陜西科學(xué)技術(shù)出版社,1997:42-44
[4]雷俊山,楊勤科.坡面薄層水流侵蝕試驗研究及土壤抗沖性評價.泥沙研究,2004(6):22-25
[5]劉青泉,李家春,陳力,等.坡面流及土壤侵蝕動力學(xué)(Ⅰ):坡面流.力學(xué)進展,2004,34(3):360-373
[6]劉青泉,李家春,陳力,等.坡面流及土壤侵蝕動力學(xué)(Ⅱ):土壤侵蝕.力學(xué)進展,2004,34(4):493-506
[7]Foster GR,Huggins LF,Meyer L D.A laboratory study of rill hydraulics:velocity relationships.Transactions of ASAE,1984,27(3):790-796
[8]Abrahams AD,LiG,Parsons J.Rill hyd raulics on a semiarid hillslope,southern Arizona.Earth surface processes and land forms,1998,21:35-47
[9]Nearing M,Brad ford M,Parker C.Soil detachment by shallow flow at low slopes.Soil Science Society of America Journal,1991,55(2):339-344
[10]沙際德,蔣允靜.試論生態(tài)侵蝕性坡面薄層水流的基本動力特征.水土保持學(xué)報,1995,9(4):29-35
[11]胡世雄,靳長興.坡面流與坡面侵蝕動力過程研究的最新進展.地理研究,1998,17(3):326-335
[12]張科利,唐克麗.黃土坡面細溝侵蝕能力的水動力學(xué)實驗研究.土壤學(xué)報,2000,37(1):9-15
[13]張光輝.坡面薄層流水動力學(xué)特性的實驗研究.水科學(xué)進展,2002,13(2):159-165
[14]Bunte K,Poesen JW.Effects of rock fragment conversation on erosion and transport of non-cohesive sediment by shallow overland flow.Water resources research,1993,29(5):1415-1424
[15]Poesen JW,Vandaele K,Wesemael B V.Gu lly erosion:importance and model implications,modeling soilerosion by water.NATO ASISeries,1998,I-55:285-311
[16]Larson W E,Lindstrom M J,Schumacher T E.The role of severe storm in soilerosion:A problem needing consideration.Journal of Soil and Water Conservation,1997,52(2):90-95
[17]Olson K P,Laflen JM,Norton L D.Evaluation ofmethods to study soil erosion-p roductivity relationships.Journal of Soil and Water Conservation,1994,49(6):586-590
[18]劉國彬,張光輝.原狀土沖刷法與人工模擬降雨法研究土壤抗沖性對比分析.水土保持通報,1996,16(2):32-37
[19]Li G,Abrahams A D,Atkinson J F.Correction factors in the determ ination ofmean velocity of overland flow.Earth surface processes and land forms,1998,21(6):509-515.
[20]馮浩,吳淑芳,吳普特,等.草地坡面徑流調(diào)控放水試驗研究.水土保持學(xué)報,2005,19(6):23-25
[21]沙玉清.泥沙運動力學(xué)引論.西安:陜西科學(xué)技術(shù)出版社,1996:72-76
[22]潘成忠,上官周平.不同坡度草地含沙水流水力學(xué)特性及其攔沙機理.水科學(xué)進展,2007,18(4):490-495
[23]Govers G.Relationships between discharge,velocity,and flow area for rill eroding loose,non-layered materials.Earth surface processes and landforms,2006,17:515-528
[24]Nearing M,Norton L,Bulgakov A,et al.Hyd raulics and erosion in eroding rills.Water resources research,1997,33(4):865-876
[25]肖培青,鄭粉莉,姚文藝.坡溝系統(tǒng)坡面徑流流態(tài)及水力學(xué)參數(shù)特征研究.水科學(xué)進展,2009,20(2):236-240
Scouring erosion experiment on bare slight slope in Ili River Basin
Ma Ping1,Song Fenghui2,Shi Yanjiang2,Wang Jian1,Zhu Shoujun1
(1.Northwest Sci-Tech University of Agriculture and Forestry,712100,Yangling,Shaanxi;2.Xinjiang Forest-Science Academy,830000,Urumqi:China)
Based on the scouring erosion experiment on bare slight slope in Ili River Basin,water infiltration and the characteristics of runoff and sediment were studied,and the hydraulic parameters characteristics of surface runoff were analyzed.The results showed:1)Soil infiltration rate took on the trend of increasing as the scouring water discharge increased;however,it had a trend of decreasing with the scouringwater amount increasing.A higher variation of soilwater had on the lower slope compared to the upper slope.2)The process of sediment discharge went up as the scouring water discharge increased,and the sediment concentration had a trend from high to low.The relation between the sediment concentration and the scouring water discharge presented a power function.W ith the scouring time continuing the sediment concentration and the sedimentdischarge both went down.3)The average flow velocity,the Reynolds number and the Darcy-Weisbach resistance coefficient increased as the scouring water discharge increased,the Froude number decreased as the scouring water discharge increased.The relation between the average flow velocity and the scouringwater discharge,the Froude number and the Reynolds number both p resented power functions.The water flow p resented a laminar,supercritical flow state.
bare slight slope;water scouring;soilwater;runoff and sediment;hydraulic parameters
2009-06-14
2009-11-20
項目名稱:國家科技支撐計劃項目“伊犁河谷水土流失綜合治理關(guān)鍵技術(shù)開發(fā)與示范”(2007BAC15B07)
馬萍(1985—),女,碩士研究生。主要研究方向:水土保持。E-mail:ping.ma.nwsuaf@gmail.com
?責(zé)任作者簡介:朱首軍(1965—),男,博士,副教授,碩士研究生導(dǎo)師。主要研究方向:水土保持。E-mail:zhushoujun@nwsuaf.edu.cn
(責(zé)任編輯:程 云)