嚴(yán)忠鑾，安瑞冬，李 嘉，李 永，鄧良勝

(1.四川大學(xué)水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室，四川成都 610065;2.四川省紫坪鋪開發(fā)有限責(zé)任公司，四川 成都 610091)

濁度型清渾水交界面識別方法及其在水庫異重流觀測中的應(yīng)用

嚴(yán)忠鑾1，安瑞冬1，李 嘉1，李 永1，鄧良勝2

(1.四川大學(xué)水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室，四川成都 610065;2.四川省紫坪鋪開發(fā)有限責(zé)任公司，四川 成都 610091)

針對傳統(tǒng)異重流觀測方法不利于對異重流進(jìn)行實(shí)時追蹤的問題，基于濁度型清渾水交界面的概念，提出濁度型清渾水交界面識別方法，并將其應(yīng)用于紫坪鋪水庫汛期渾水異重流觀測。該方法利用濁度值突變能判定異重流交界面位置，利用濁度型清渾水交界面描述異重流前鋒厚度和濃度外邊界線，并將異重流交界面位置確定在1.2～3.5 m厚的水層內(nèi)，同時還分析了地形變化對異重流運(yùn)動規(guī)律的影響。異重流原型觀測結(jié)果表明，濁度型清渾水交界面識別方法可以作為水庫汛期監(jiān)測異重流的有效手段。

渾水異重流;清渾水交界面;水庫觀測;水溫

汛期暴雨沖刷河道兩岸，挾帶大量泥沙形成渾水入庫，在適當(dāng)條件下，渾水會潛入庫底形成水庫水沙異重流，并可能運(yùn)動到壩前，通過泄水建筑物排出庫外。但若后續(xù)能量不足，異重流攜帶的泥沙將沿程落淤，減小水庫有效庫容，長此以往勢必影響水庫的使用壽命。因此，研究異重流的運(yùn)動特性能夠?qū)λ畮斓挠俜e情況進(jìn)行預(yù)測，為實(shí)現(xiàn)人工塑造異重流以及水庫異重流多目標(biāo)調(diào)度提供科學(xué)依據(jù)，具有重要的工程意義和學(xué)術(shù)價值。

國內(nèi)外針對異重流的研究開展已久，主要的研究成果集中在模型試驗(yàn)和數(shù)值模擬方面。Alexander等［1］在2002年利用玻璃水槽，通過改變?nèi)肓髁髁亢退燮露冗M(jìn)行異重流試驗(yàn)，得到前鋒厚度與流量、坡度、粒徑及含沙量的關(guān)系;An等［2］在2010年針對渾水異重流潛入點(diǎn)的水力學(xué)特性，在兩相流玻璃水槽中進(jìn)行了16組試驗(yàn)，利用粒子成像測速(PIV)對渾水異重流潛入點(diǎn)附近的流場進(jìn)行了測量，得到了潛入點(diǎn)附近水體的反向速度及回流圖像;Francois等［3］在2004年采用二維數(shù)學(xué)模型模擬了不同底坡傾角條件下的異重流運(yùn)動情況，計(jì)算得到異重流運(yùn)動速度和水流定性特征，與異重流試驗(yàn)現(xiàn)象描述一致;李永等［4］在 2009年采用基于兩相流理論的ASM(algebraic slip mixture)模型，對異重流的潛入、前鋒形態(tài)及推進(jìn)速度進(jìn)行了模擬。目前國內(nèi)已有一些高含沙河流異重流原型觀測方面的成果，如官廳水庫、紅山水庫、小浪底水庫等［5-7］，但鮮有低含沙量河流水庫異重流運(yùn)動觀測成果的報道，究其原因主要是觀測手段的缺乏。對于快速、動態(tài)的水沙異重流運(yùn)動變化過程，傳統(tǒng)方法難以進(jìn)行實(shí)時追蹤，這無疑會延誤異重流的觀測時機(jī)。2004年黃河第3次調(diào)水調(diào)沙試驗(yàn)中采用了光電式清渾水界面探測器對異重流進(jìn)行觀測［8］，該儀器的使用無疑提高了水庫異重流的測量精度和效率。但將清渾水交界面臨界含沙量定義為1kg/m3，含沙量測量精度為±0.2 kg/m3，對于高含沙異重流監(jiān)測追蹤有較強(qiáng)的適用性，而對于低含沙河流異重流的原型觀測還需驗(yàn)證。為此，本文根據(jù)濁度與含沙量的關(guān)系，提出一種適用于低含沙河流異重流原型觀測的方法，該技術(shù)可以快速地確定異重流交界面位置，該方法在紫坪鋪水庫原型觀測中取得了較好的應(yīng)用效果。

1 濁度型清渾水交界面

1.1 定義

在水庫中，當(dāng)渾水初遇清水時，兩種水體的碰撞使得垂向斷面上含沙量分布均勻。隨著后續(xù)渾水的補(bǔ)充及粗顆粒泥沙受重力作用向下運(yùn)動，清渾水體間密度差增大，壓力差在斷面底部達(dá)到最大，使得渾水潛入清水底部形成異重流。至此，交界面由潛入點(diǎn)開始形成，并隨著異重流的運(yùn)動而發(fā)展。而事實(shí)上，交界面并非一個截然分明的面，而是具有一定厚度、密度及流速等物理量梯度分布的過渡層。

在異重流形成后的過渡階段，渾水具有較大動能，含沙量垂向分布均勻，交界面上泥沙擴(kuò)散作用強(qiáng)烈，過渡層較厚且其邊界模糊不清;當(dāng)異重流運(yùn)動達(dá)到穩(wěn)定后，交界面上摻混達(dá)到一定的平衡狀態(tài)，過渡層的邊界逐漸變得清晰。但無論過渡層的厚薄如何，在清水層、過渡層與渾水層之間濁度必然存在沿深度方向的變化。由此本文定義，由過渡層內(nèi)濁度突變位置構(gòu)成的交界面為濁度型清渾水交界面。若在異重流原型觀測中，能通過濁度對異重流進(jìn)行快速定位，就可以迅速確定潛入點(diǎn)位置、前鋒位置及前鋒厚度，對于異重流原型觀測和水庫利用異重流進(jìn)行排沙工作有一定的參考意義。

1.2 識別技術(shù)

以往在水沙異重流原型觀測中，通常以含沙量的變化來分析異重流的運(yùn)動規(guī)律，且異重流渾水水樣的采集和含沙量的測定是異重流原型觀測的重要內(nèi)容。但現(xiàn)場水沙異重流水樣采集工作難度大且耗費(fèi)時間，水樣經(jīng)長時間擱置又容易破壞含沙的原始狀態(tài)，再經(jīng)過實(shí)驗(yàn)室采用如振動管法［9］、超聲波法［10］及 γ射線法［11］等方法分析，對含沙量的測量結(jié)果影響大、效率低且不便于異重流觀測時現(xiàn)場使用。同樣現(xiàn)場測量時交界面位置是通過視覺觀測及人工判斷，這樣做顯然不夠科學(xué)?？梢妭鹘y(tǒng)的異重流觀測方法在觀測精度和效率方面有所欠缺，對于異重流的快速變化很難及時作出反應(yīng)，容易延誤利用異重流清淤排沙的時機(jī)。

本文基于前期研究中進(jìn)行的含沙量與濁度關(guān)系的試驗(yàn)［12］，開展紫坪鋪水庫含沙量與濁度關(guān)系的試驗(yàn)。試驗(yàn)采集岷江紫坪鋪庫區(qū)泥沙，在實(shí)驗(yàn)室內(nèi)烘干，再充分?jǐn)嚢韬笈渲贸?0種含沙量低于2.5 kg/m3的水樣進(jìn)行試驗(yàn)，得到結(jié)果如圖1所示。

圖1 含沙量與濁度的關(guān)系

由圖1可以看出，含沙量與濁度具有很好的線性關(guān)系，根據(jù)試驗(yàn)數(shù)據(jù)可擬合出含沙量與濁度的線性方程如下:

式中:T為濁度，NTU;C為含沙量，mg/L。

濁度型清渾水交界面識別方法正是基于清渾水交界面的研究及含沙量與濁度之間的關(guān)系而建立的觀測方法。應(yīng)用濁度傳感器可自動記錄斷面上垂向濁度分布，通過判斷數(shù)據(jù)突變位置，來實(shí)現(xiàn)清渾水交界面位置的確定，繼而快速反映異重流潛入點(diǎn)位置、前鋒厚度等要素，從而提高異重流的觀測效率。因此，濁度型清渾水交界面識別方法在水庫異重流原型觀測中的應(yīng)用具有可行性。

異重流觀測采用走航方式進(jìn)行，利用全球定位系統(tǒng)GPS進(jìn)行庫區(qū)測量點(diǎn)的定位，采用激光測距儀測量水面寬度。在應(yīng)用濁度傳感器和溫度傳感器記錄斷面上垂向濁度與溫度分布的同時，進(jìn)行水樣采集工作。應(yīng)用濁度傳感器時采用90°散射法，通過觀測水體中散射光線強(qiáng)度來測量水體的渾濁程度，測量范圍為0～1000 NTU，分辨率為0.1 NTU，精度為±5%。

2 水庫異重流觀測

紫坪鋪水利樞紐位于四川省成都市西北60余km的岷江上游都江堰市麻溪鄉(xiāng)，水庫總庫容11.12億m3，工程壩址以上控制流域面積22662km2，多年平均流量469 m3/s，年徑流總量148億m3，水庫水溫呈分層結(jié)構(gòu)，屬不完全年調(diào)節(jié)型水庫。1955—1990年多年平均懸移質(zhì)輸沙量為792萬t，多年平均含沙量為0.572 kg/m3，汛期(7—9月)輸沙量占年輸沙量的98.3%。

由于岷江上游普降暴雨，觀測當(dāng)日紫坪鋪入庫流量由370 m3/s增加到640 m3/s，洪水挾沙入庫往往是異重流運(yùn)動的誘發(fā)因素，本次水庫異重流觀測斷面共有12個，其布置見圖2，斷面距壩里程見表1。

表1 斷面距壩里程

圖2 觀測斷面布置示意圖

2.1 異重流潛入段特性

觀測船行至庫尾1號斷面時，觀測到該斷面上游有大面積的垃圾漂浮帶，且表面水體渾濁，該現(xiàn)象與潛入點(diǎn)的標(biāo)志相符。結(jié)合1號、2號和3號斷面的觀測數(shù)據(jù)分析潛入段的濁度和溫度變化，如圖3和圖4所示。通過分析垂向濁度分布數(shù)據(jù)可以發(fā)現(xiàn)，在1號斷面，水深1.5～2.7m之間的濁度由34NTU增大到131NTU，濁度變化梯度為78.6 NTU/m，根據(jù)清渾水交界面定義，可以判斷該斷面的清渾水交界面位置在水深1.5～2.7 m之間;同樣可以判斷2號與3號斷面的清渾水交界面位置分別在水深2.0～3.6 m和1.3～2.5 m范圍內(nèi)，3個斷面間異重流厚度沿程增大。分析水溫分布數(shù)據(jù)可以發(fā)現(xiàn)，在1號斷面，水深1.5 m處的水溫突然升高，至6.6 m處又開始降低;在2號和3號斷面，水溫的影響主要集中在斷面底部，底部溫度與其斷面上溫度拐點(diǎn)處的溫差均達(dá)到0.2℃，同時在溫度拐點(diǎn)處的濁度也發(fā)生較大變化。由此可以推斷挾沙水流在1號斷面上游已經(jīng)潛入并形成異重流，再經(jīng)過2號斷面，運(yùn)動到3號斷面時已經(jīng)趨于穩(wěn)定，且異重流在影響庫區(qū)垂向水體含沙量分布的同時也對水溫結(jié)構(gòu)造成破壞。

圖3 潛入段各斷面垂向濁度分布

圖4 潛入段各斷面垂向水溫分布

2.2 異重流潛流段及前鋒段特性

7～12號斷面的垂向濁度分布如圖5所示。7、8、9號斷面的垂向濁度突變點(diǎn)位置分別處于水深14.2～16.2 m、22.3～24.6 m 和24.3～27.8 m 的范圍內(nèi)，顯然清渾水交界面位置在潛流段沿程不斷變深。在10號斷面上，濁度突變位置的深度明顯減小，清渾水交界面位于水深18.6～20.3 m之間，而11號和12號斷面的距壩里程分別為5.5 km和1.3 km，兩斷面上垂向濁度最大值僅為32.7NTU，而9號和10號兩斷面上垂向濁度最大值分別達(dá)到1036 NTU和810 NTU，遠(yuǎn)大于11號與12號兩斷面的濁度最大值，因此可以推斷異重流前鋒已經(jīng)運(yùn)動到距壩5.5～7.4 km之間的河段，11號與12號兩斷面同一深度的水體濁度相差不到3 NTU，可以肯定11號斷面的水體未受到異重流的影響且異重流前鋒所在位置更靠近10號斷面。

圖5 潛流段和前鋒段各斷面垂向濁度分布

各段面的垂向溫度分布如圖6所示，分析可以發(fā)現(xiàn)7～10號斷面上垂向水溫拐點(diǎn)位置在水深14.2～25.1 m之間變化，而除10號斷面外，7～9號斷面的溫度拐點(diǎn)位置與其濁度突變位置同步變化，說明10號斷面水溫受異重流影響的時間較短，庫底水溫已受異重流影響而升高，中層水體水溫還未發(fā)生變化或受影響較小。從庫底水溫變化來看，11號斷面在水深35m處的水溫比10號斷面低0.8℃，而與7號斷面之間的溫差可達(dá)1.25℃。由此可見，異重流運(yùn)動對水庫水溫結(jié)構(gòu)的影響由庫底向上發(fā)展，受異重流影響時間越長，水溫結(jié)構(gòu)破壞程度越大。

圖6 潛流段和前鋒段各斷面垂向水溫分布

2.3 地形變化的異重流特性

庫區(qū)濁度、水深與水面寬度的沿程變化如圖7所示。處理數(shù)據(jù)過程中發(fā)現(xiàn)在濁度突變點(diǎn)的水深范圍內(nèi)，相應(yīng)的濁度值在50～150 NTU范圍內(nèi)變動且大部分濁度值接近100 NTU，因此以濁度值為100 NTU的等值線作為此次原型觀測中的清渾水交界面。從1號到4號斷面水庫庫底坡度基本不變，但水面寬度變化較大。在從1號到3號斷面的運(yùn)動過程中，異重流運(yùn)動速度較快，因此其受水面寬度的影響很小，很快趨于穩(wěn)定并保持較好的運(yùn)動形態(tài)，而異重流厚度增大有可能是因較大的流速將沿程淤積的泥沙從庫底卷進(jìn)異重流所致。而從3號到4號斷面沿程為深切的V形河谷，由寬河段進(jìn)入窄河段，受到河槽束窄的作用，異重流的厚度有所增加［13］。

圖7 庫區(qū)濁度、水深與水面寬度的沿程變化

4號到6號斷面之間的河段由壽溪河支流從5號斷面匯入，4號斷面靠近匯口，匯口以下至9號斷面河段水面寬度基本不變，河道相對順直，而4號到6號斷面間異重流厚度的減小，可能是異重流在流經(jīng)4號斷面后于匯口處遇支流匯入，水體碰撞損失部分能量，而水面突然拓寬也降低了流速，水流挾沙能力下降，異重流厚度變薄，而進(jìn)入順直河道后，異重流流速開始增大并形成穩(wěn)定的異重流形態(tài)，但相對4號斷面，異重流厚度減小。6號到9號斷面，水面寬度基本不變，異重流在形態(tài)上趨于穩(wěn)定，異重流厚度隨水深增大而增加。在異重流由9號斷面運(yùn)動至10號斷面的過程中，異重流核心不再沿庫底運(yùn)動，斷面上垂向濁度呈現(xiàn)先增加再減小的趨勢，原因可能是該河段底坡較7號到9號斷面間的河段底坡突然增大，水深增加明顯，而異重流核心部分慣性力較大，依舊能保持原先的運(yùn)動趨勢。

2.4 低含沙量異重流

前文所提及的官廳水庫、紅山水庫及小浪底水庫的多年平均含沙量均大于15 kg/m3［14］，而位于西南山區(qū)岷江上游的紫坪鋪水庫多年平均含沙量僅為0.572 kg/m3。與文獻(xiàn)［15］中小浪底異重流形成的含沙量下限條件(10 kg/m3)比較，位于低含沙河流岷江上游的紫坪鋪水庫的異重流形成時的含沙量要小很多，原因可能是西南山區(qū)天然落差大，岷江上游340 km長河段的天然落差為3009 m，渾水來流在高坡降的河道上運(yùn)動后，在進(jìn)入紫坪鋪水庫時已具有較大的動能，使得渾水潛入形成異重流。可見，該原型觀測在驗(yàn)證濁度型清渾水交界面識別技術(shù)的適用性的同時，填補(bǔ)了西南山區(qū)低含沙量河流水庫異重流觀測資料的空白，進(jìn)一步補(bǔ)充了低含沙量河流異重流觀測資料，對低含沙量河流水庫管理及泥沙運(yùn)動研究有一定的參考意義。

3 結(jié)論

a.基于濁度型清渾交界面概念的交界面識別方法在紫坪埔水庫異重流觀測中應(yīng)用效果良好，成功地捕捉到了低含沙量條件下的異重流現(xiàn)象，并將清渾水交界面的位置確定在1.2～3.5 m厚的水層內(nèi)，同時這也是首次在西南山區(qū)低含沙河流水庫中觀測到異重流現(xiàn)象，是對該區(qū)域異重流資料空白的填補(bǔ)。

b.通過濁度分析發(fā)現(xiàn)，在異重流從潛入到穩(wěn)定的過程中，異重流的厚度受河段水面寬度影響較小，而流速增大使其挾沙能力增強(qiáng)，能將沿程淤積的泥沙卷進(jìn)渾水內(nèi)，并增加自身渾水層厚度;穩(wěn)定后的異重流在支流匯入后，由于水體的碰撞及河寬突擴(kuò)的影響，異重流厚度減小;在異重流能量足夠的情況下，異重流能迅速恢復(fù)先前運(yùn)動形態(tài)。

c.溫度數(shù)據(jù)的分析顯示汛期高溫異重流會破壞水庫的水溫結(jié)構(gòu)，先造成水庫底層水溫升高，水溫影響范圍從庫底向上發(fā)展，且對于受影響較久的水體，溫度拐點(diǎn)靠近濁度突變點(diǎn)位置。

［1］ALEXANDER J，MULDER T.Experimental quasi-steady density currents ［J］.Marine Geol，2002，186(3):195-210.

［2］AN Ruidong，LI Jia.Characteristic analysis of turbidity currents plunging［J］.Journal of Hydrodynamics:Ser B，2010，22(2):274-282.

［3］FRANCOIS B，VINCENT P，ECKART M，et al.Evaluation of simplified approach for simulating gravity currents over slopes of varying angles［J］.Computers ＆ Fluids，2006，35(5):492-500.

［4］李永，李嘉，安瑞冬.水沙兩相流ASM模型在渾水異重流計(jì)算中的應(yīng)用及模型試驗(yàn)研究［J］.四川大學(xué)學(xué)報:工程科學(xué)版，2009，41(4):102-108.(LI Yong，LI Jia，AN Ruidong.The application of water-sediment two phase flow ASM model in the density current of turbid water and related modelteststudy［J］.JournalofSichuan University:Engineering Science Edition，2009，41(4):102-108.(in Chinese))

［5］侯暉昌，焦恩澤，秦芳.官廳水庫1953—1956年異重流資料初步分析［J］.泥沙研究，1958(2):70-94.(HOU Huichang，JIAO Enze，QIN Fang.The preliminary analysis of the turbidity current datas from 1953 to 1956 in Guanting Reservoir［J］.Journal of Sediment Research，1958(2):70-94.(in Chinese))

［6］曹如軒.天然高含沙水流流變與流態(tài)的觀測分析［J］.人民黃河，1983(1):20-24.(CAO Ruxuan.The observation and analysis of rheology and pattern of the natural hyperconcentrated flow［J］.Yellow River，1983(1):20-24.(in Chinese))

［7］侯素珍，焦恩澤，林秀芝，等.小浪底水庫異重流運(yùn)動特

征分析［J］.泥沙研究，2004(5):46-50.(HOU Suzhen，JIAO Enze，LIN Xiuzhi，et al.Analysis on the movement characters of the density current in Xiaolangdi Reservoir［J］.Journal of Sediment Research，2004(5):46-50.(in Chinese))

［8］王慶中，陳富安，王秀清，等.清渾水界面探測器在水庫異重流測驗(yàn)中的應(yīng)用［J］.人民黃河，2006，28(2):67-68.(WANG Qingzhong，CHEN Fu’an，WANG Xiuqing，et al.The application of interface probe of clear-muddy water in the field turbid current test in reservoir［J］.Yellow River，2006，28(2):67-68.(in Chinese))

［9］李廷生.振動管式含沙量測量原理研究［J］.泥沙研究，1985(3):26-35.(LI Tingsheng.A study on the principle of measuring sediment concentration by vibratory pipe［J］.Journal of Sediment Research，1985(3):26-35.(in Chinese))

［10］方彥君，唐懋官.超聲衰減法含沙量測量研究［J］.泥沙研究，1990(2):1-12.(FANG Yanjun，TAN Maoguan.Ultrasonic attenuation method for measuring sediment concentration［J］.Journal of Sediment Research，1990(2):1-12.(in Chinese))

［11］方彥君，張紅梅，程瑛.含沙量測量的新進(jìn)展［J］.武漢水利電力大學(xué)學(xué)報，1999，32(3):55-57.(FANG Yanjun，ZHANG Hongmei，CHENG Ying.Progress of sediment concentration measurement［J］.Engineering Journal of Wuhan University，1999，32(3):55-57.(in Chinese))

［12］彭永勤，李克鋒，曹薇.利用渾濁度推算含沙量方法的限制性與適用性研究［J］.四川環(huán)境，2009，28(1):20-23.(PENG Yongqin，LI Kefeng，CAO Wei.The restriction and applicability for the method of using turbidity to deduce sand content［J］.Sichuan Environment，2009，28(1):20-23.(in Chinese))

［13］范家驊.渾水異重流槽寬突變時的局部摻混［J］.水利學(xué)報，2005，36(1):1-8.(FAN Jiahua.Local entrainment of turbid density current in width abrupt changed channel［J］.Journal of Hydraulic Engineering，2005，36(1):1-8.(in Chinese))

［14］姜乃森.官廳水庫淤積上延問題的初步分析［J］.泥沙研究，1985(1):61-69.(JIANG Naisen.A preliminary analysis of deposition-extension in Guanting Reservoir［J］.Journal of Sediment Research，1985(1):61-69.(in Chinese))

［15］李書霞，張俊華，陳書奎，等.小浪底水庫塑造異重流技術(shù)及調(diào)度方案［J］.水利學(xué)報，2006，37(5):567-572.(LI Shuxia，ZHANG Junhua，CHEN Shukui，et al.Formation of density current by means of reasonable operation of reservoir in Xiaolangdi Project［J］.Journal of Hydraulic Engineering，2006，37(5):567-572.(in Chinese))

Turbid interface identification method of turbidity currents and its application in field observation at Zipingpu

Reservoir

YAN Zhongluan1，AN Ruidong1，LI Jia1，LI Yong1，DENG Liangsheng2(1.State Key Laboratory of Hydraulics and Mountain River Engineering，Sichuan University，Chengdu610065，China;2.Sichuan Province Zipingpu Development Co.，Ltd.，Chengdu610091，China)

The turbidity currents are observed in traditional methods which aren’t suitable for the real-time observation.Therefore，based on the concept of the turbidity-interface between turbid water and clear water，this paper proposes an interface identification method that can be applied in the field observation of turbidity currents at reservoir during floods.By analyzing the sudden change in water turbidity，this method can locate the position of the interface within a water layer thick of about 1.2 m to 3.5 m.Then the front thickness and outer boundary of transition region can be illustrated by using the concept of the interface proposed in the present paper.In addition，the paper also analyzed the impact on the performances of turbidity currents caused by topography variation.On all accounts，the results show that the method is effective for monitoring the turbidity currents at reservoir during the floods.

turbidity currents;interface between turbid water and clear water;reservoir observation;water temperature

TV145+.24

A

1006-7647(2013)06-0071-05

10.3880/j.issn.1006-7647.2013.06.015

國家自然科學(xué)基金(51109148);四川大學(xué)青年教師科研啟動基金(2011SCU11113)

嚴(yán)忠鑾(1988—)，男，福建連江人，碩士研究生，主要從事環(huán)境水力學(xué)研究。E-mail:yanzl0213@163.com

安瑞冬(1983—)，男，天津人，副研究員，博士，主要從事環(huán)境水力學(xué)研究。E-mail:anruidong@scu.edu.cn

2012-01-07 編輯:周紅梅)