李飛鵬，高 雅，張海平**，肖宜華，陳 玲

(1: 同濟(jì)大學(xué)環(huán)境科學(xué)與工程學(xué)院，上海 200092)(2: 同濟(jì)大學(xué)建筑設(shè)計(jì)研究院(集團(tuán))有限公司，上海 200092)

流速對浮游藻類生長和種群變化影響的模擬試驗(yàn)*

李飛鵬1，高 雅2，張海平1**，肖宜華1，陳 玲1

(1: 同濟(jì)大學(xué)環(huán)境科學(xué)與工程學(xué)院，上海 200092)(2: 同濟(jì)大學(xué)建筑設(shè)計(jì)研究院(集團(tuán))有限公司，上海 200092)

采集夏季崇明島中心湖原水，在環(huán)形有機(jī)玻璃水槽開展了不同流速對浮游藻類生長和種群變化影響的研究.相對于靜止水槽，實(shí)驗(yàn)前期不同的流速條件均對浮游藻類的生長產(chǎn)生了一定的抑制作用，實(shí)驗(yàn)后期浮游藻類葉綠素a含量均穩(wěn)定在一定水平，且流動(dòng)水槽中葉綠素a含量略大于靜止水槽，這可能與靜止圍隔中出現(xiàn)大量枝角類浮游動(dòng)物有關(guān)，而流動(dòng)條件則抑制了浮游動(dòng)物的生長，減輕了對浮游藻類的捕食壓力.水體流動(dòng)導(dǎo)致了浮游藻類種群的變化，藍(lán)藻的迅速消失、綠藻和硅藻形成生長優(yōu)勢是主要特征，表明持續(xù)的流動(dòng)條件是導(dǎo)致浮游植物種屬減少和群落結(jié)構(gòu)變化的直接原因.該研究可為調(diào)水引流在河道、湖泊和水庫中的富營養(yǎng)化控制和水華防治提供基礎(chǔ)依據(jù).

流速；浮游藻類；種群變化；模擬試驗(yàn)

流速是水動(dòng)力條件中最基本、最直觀的參數(shù)，對浮游植物的生長和繁殖具有十分重要的影響，不僅影響水體中能量和營養(yǎng)鹽的分配、溶解氧水平、沉積物特性等，而且能夠?qū)λ锏纳砘顒?dòng)產(chǎn)生直接的影響.天然河流、湖泊和水庫中，有關(guān)流速與浮游藻類生物量統(tǒng)計(jì)關(guān)系的研究也常見報(bào)道，流速、流量等水動(dòng)力條件已經(jīng)成為浮游植物生長繁殖、群落更替和誘發(fā)水華現(xiàn)象的重要驅(qū)動(dòng)因子[1-2].對浮游植物個(gè)體而言，每種生命體對其生存的流場環(huán)境都有特定的要求，當(dāng)條件改變時(shí)，生命體都會表現(xiàn)出自由、逃離和失控3種狀態(tài)[3]. 因此，通過研究浮游植物個(gè)體在不同水動(dòng)力條件下的生命規(guī)律，局部改變有害藻類生存的流場，達(dá)到抑制水華暴發(fā)的目的，具有十分重要的意義.

由于自然水體中藻類生長的影響因素較為復(fù)雜，流速效應(yīng)經(jīng)常被其他各種因素所掩蓋而導(dǎo)致研究結(jié)論不明確.室內(nèi)模擬實(shí)驗(yàn)?zāi)軌蛟谂懦渌蛩赜绊懙那闆r下較為精細(xì)地反映流速與浮游植物之間的關(guān)系[4-6]，許多研究者通過研究流速對單一藻種或者混合藻種的影響，認(rèn)為存在一定的臨界流速，當(dāng)流速大于或者小于臨界流速時(shí)，浮游植物的生長速率都會減緩[7].但受空間和時(shí)間尺度限制，其研究結(jié)論缺乏普遍性，結(jié)果差異性較大，仍需大量更大尺度的研究[8-10].朱宜平等[11]、Li等[12]通過野外圍隔實(shí)驗(yàn)，發(fā)現(xiàn)0.03～0.30m/s的流速均對葉綠素a濃度產(chǎn)生抑制作用，并且顯著降低藻類豐度；但受底泥擾動(dòng)和降雨等因素影響，再加上水動(dòng)力產(chǎn)生的剪切力極大地促進(jìn)了水綿的營養(yǎng)生殖，流速效應(yīng)同樣受到一定程度的掩蓋.因此，本研究在野外模擬實(shí)驗(yàn)的基礎(chǔ)上，采用較大尺度的有機(jī)玻璃水槽來模擬夏季富營養(yǎng)水體中流速對浮游藻類生長和種群變化的影響，以期準(zhǔn)確反映流速對浮游植物生長的抑制效果和種群變化規(guī)律.

1 材料與方法

中心湖位于上海市崇明縣豎新鎮(zhèn)前衛(wèi)村，該水體毗鄰?fù)瑵?jì)大學(xué)長江水環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室崇明水環(huán)境研究中心，是由5條相互貫通的河段組成的封閉水域，面積約52000m2，水深1.0～1.3m，全湖總氮(TN)濃度為0.53～4.14mg/L，總磷(TP)濃度為0.056～0.750mg/L. 該水體封閉性較好且具有相對簡單的生態(tài)系統(tǒng)，十分便于野外密集觀測和現(xiàn)場實(shí)驗(yàn)研究.具體信息可參考文獻(xiàn)[13].

實(shí)驗(yàn)裝置采用3個(gè)環(huán)形有機(jī)玻璃水槽，如圖1所示，四周為透明有機(jī)玻璃板，底部為加厚PVC板，高1.3m，置于崇明水環(huán)境研究中心大廳內(nèi)，避開了底泥擾動(dòng)起浮和降雨帶來的影響.從早上至上午10:30左右裝置能接受陽光照射，白天時(shí)段通過日光燈提供光照條件，日光燈于早上6:00開啟，晚上18:00關(guān)閉.實(shí)驗(yàn)用水為中心湖原水，于2010年8月13日取水約480L置于裝置內(nèi)，用2種不同型號潛水泵(堵住增氧口)推流，其中一個(gè)裝置保持靜止，作為對比參照.3個(gè)裝置分別標(biāo)記為S0、S1和S2，對應(yīng)的平均驅(qū)動(dòng)流速分別為0、0.06和0.10m/s，實(shí)驗(yàn)周期約為13d.同時(shí)實(shí)驗(yàn)中還采集中心湖西部的1個(gè)監(jiān)測點(diǎn)(W2)作為對照.實(shí)驗(yàn)期間向3個(gè)裝置中分別加入200ml(15g/L的NaNO3和3g/L的KH2PO4)營養(yǎng)鹽溶液，隔天測1次營養(yǎng)鹽濃度，保證裝置中氮、磷營養(yǎng)鹽濃度維持在高富營養(yǎng)水平(TN濃度大于2mg/L和TP濃度大于0.2mg/L)，自8月15日開始，每天早上(8:00)和晚上(17:00)分別監(jiān)測1次，記錄水溫、溶解氧(DO)、pH、電導(dǎo)率等，并采樣測定葉綠素a(Chl.a)濃度和藍(lán)藻、綠藻、硅藻的生物量比例.

圖1 圍隔實(shí)驗(yàn)裝置及示意圖Fig.1 Enclosure experimental apparatus and its diagram

2 結(jié)果與討論

2.1 物化因素分析

表1 實(shí)驗(yàn)期間不同流速下水溫、pH、電導(dǎo)率和濁度的日變化趨勢

水體流動(dòng)對裝置內(nèi)DO影響較大，DO的日變化如圖2a所示.S1和S2水體中DO濃度在前2d發(fā)生驟降，自第3d開始受水體流動(dòng)影響逐步變大，第5d開始基本保持穩(wěn)定. S2中DO濃度略高于S1，兩者穩(wěn)定后的平均DO濃度分別為7.3和7.7mg/L，顯著高于靜止水槽(DO平均濃度為4.5mg/L).靜止裝置中DO變化波動(dòng)明顯，可能受浮游植物生長影響較大.

2.2 流速對浮游藻類生長的影響

Chl.a濃度變化趨勢如圖2b所示，初始Chl.a濃度為44.22μg/L.浮游藻類從自然水體移到小尺度的水槽，必然要經(jīng)歷一個(gè)適應(yīng)的過程.實(shí)驗(yàn)最初階段，在0.06和0.10m/s流速下，Chl.a濃度在第2d均增加了約1.2倍，靜止水槽中藻類的增長并不明顯，這是由于水流加速了營養(yǎng)鹽向藻細(xì)胞的擴(kuò)散，顯著刺激了浮游藻類生長，此時(shí)水流引起的營養(yǎng)鹽向藻細(xì)胞的擴(kuò)散作用為控制因素.流速對浮游植物的抑制效果在第3～8d表現(xiàn)明顯，靜止圍隔中Chl.a濃度在實(shí)驗(yàn)開始前3d有所降低，之后出現(xiàn)大幅升高，于第6d下午達(dá)到最高值55.24μg/L，隨后逐漸降低至10μg/L左右；兩種流速下浮游藻類的生長規(guī)律相對一致，Chl.a濃度差異不明顯，自第3d開始出現(xiàn)明顯下降，至第6d和第7d降低至10μg/L，隨后小幅升高，穩(wěn)定在20～30μg/L左右.實(shí)驗(yàn)第8d以后，靜止裝置中Chl.a濃度出現(xiàn)明顯下降，這與裝置中出現(xiàn)大量的浮游動(dòng)物關(guān)系密切.鏡檢結(jié)果發(fā)現(xiàn)，此時(shí)裝置中水體清澈透明，藻細(xì)胞數(shù)量很少，枝角類浮游動(dòng)物增多，平均密度超過80ind./L，遠(yuǎn)大于S1和S2中的平均密度20ind./L.朱宜平等[11]在流速為0.002m/s的野外圍隔實(shí)驗(yàn)條件下，發(fā)現(xiàn)水體中有大量的枝角類動(dòng)物生長，主要種類為秀體溞屬(Diaphanosoma)和溞屬(Daphnia),能夠?qū)⒏∮卧孱惒妒炒M.陳偉民等[10]認(rèn)為相對于大水流狀態(tài)，在靜止?fàn)顟B(tài)和微小流動(dòng)狀態(tài)下，枝角類種類和數(shù)量會明顯增多，且以大體型種類為主，常以浮游單細(xì)胞藻類、有機(jī)碎屑等為食，受到這些浮游甲殼動(dòng)物的濾食，浮游藻類的數(shù)量和生物量將大幅降低.在0.06和0.10m/s流速下，Chl.a濃度除在后期出現(xiàn)了小幅度升高外，基本保持穩(wěn)定，這可能是由于在流動(dòng)條件下，流速抑制了浮游動(dòng)物的生長，減小了浮游藻類的被捕食壓力，促進(jìn)了流動(dòng)系統(tǒng)中生態(tài)系統(tǒng)的相對穩(wěn)定.從生物量的角度計(jì)算，在這種尺度的實(shí)驗(yàn)條件下，流速對浮游藻類生長的相對抑制率約可達(dá)到50%.根據(jù)潛水泵自身帶來的剪切效應(yīng)對藻類細(xì)胞的影響，估算由此帶來的抑制率在10%左右[12].

圖2 無底泥圍隔中不同流速下DO(a)和Chl.a(b)的變化趨勢

流速對水體DO濃度的影響是本研究的重要影響因素.實(shí)驗(yàn)過程中堵住潛水泵的增氧口，使得潛水泵的主要功能為制造水流.前2d流動(dòng)水槽中DO濃度有一定的驟降，這與第2d浮游植物的驟增進(jìn)而會消耗大量DO有關(guān).在水流效應(yīng)下，水層之間的攪動(dòng)形成較強(qiáng)的擴(kuò)散梯度，將造成大氣中DO的快速補(bǔ)充，因此自第3d開始DO濃度迅速上升且基本保持穩(wěn)定，這說明水體流動(dòng)引起的DO濃度增加起主導(dǎo)作用，同時(shí)也說明水體中浮游藻類的光合作用能力可能受到了水流的影響；靜止水槽中，DO濃度呈W形波動(dòng)，與該水槽中浮游藻類的Chl.a濃度變化趨勢基本一致，DO的上升與降低略微遲滯于Chl.a濃度的變化，這是由于靜止條件下DO的擴(kuò)散較慢.由于靜止水體中DO濃度主要受浮游藻類光合作用產(chǎn)氧的影響，同時(shí)由于有機(jī)物、浮游動(dòng)物等耗氧，DO濃度經(jīng)常處于較低水平，但靜止水體中Chl.a濃度在第3～8d顯著高于流動(dòng)水體，說明流動(dòng)水槽中流速對浮游藻類的生長產(chǎn)生的抑制作用占據(jù)主導(dǎo)，DO對浮游藻類生長的作用弱于流速的作用.

2.3 流速對浮游藻類種群變化的影響

靜止裝置(S0)和流動(dòng)裝置(S1和S2)中浮游植物種群變化差異明顯，優(yōu)勢種屬變化如表2所示.運(yùn)行前，各裝置優(yōu)勢藻種相同，主要為藍(lán)藻門的微囊藻屬、綠藻門中的纖維藻屬和硅藻中的小環(huán)藻屬.S1和S2形成持續(xù)流動(dòng)后，第5天即觀察到微囊藻消失，且浮游藻種類越來越少，最后以綠藻為主，而靜止圍隔中，微囊藻直到第10d一直存在，這與每日進(jìn)行的Phyto-PAM檢測的種屬變化結(jié)果比較一致(圖3).

在中心湖(W2)監(jiān)測點(diǎn)，前10d主要是藍(lán)藻和硅藻為優(yōu)勢種，平均生物量所占的比例分別為47.7%和52.3%，后3d硅藻少量出現(xiàn)，優(yōu)勢種99%以上為藍(lán)藻，這與微囊藻的生理特性是相符的；S0中的藻種變化在前10d與W2相仿，后3d僅存綠藻中纖維藻屬，S1和S2中的藍(lán)藻分別在第5d和第4d全部消失，之后兩個(gè)裝置中均以綠藻為主、硅藻為輔，此時(shí)S1中綠藻和硅藻的平均生物量比例為77.1%和22.9%，S2中的平均生物量比例為83.7%和16.3%，S2中的環(huán)境條件相對更有利于綠藻生長(圖3).但靜止水槽中，第11d開始微囊藻很快消失.S0、S1和S2后期均以綠藻門中的纖維藻屬為主，這可能與浮游動(dòng)物的選擇性捕食有關(guān).從靜止?fàn)顟B(tài)向持續(xù)流動(dòng)狀態(tài)，浮游植物種屬優(yōu)勢種基本從藍(lán)藻向綠藻轉(zhuǎn)化，這與國外在淡水河流[14]、湖泊[15]中的研究結(jié)論是一致的.從生物學(xué)的意義來看，多數(shù)藍(lán)藻具有氣囊，能夠在水體中垂直移動(dòng)，獲取有利于自身增殖的營養(yǎng)鹽和光照條件，較適宜在靜止或弱紊動(dòng)水體中生長，而多數(shù)綠藻和硅藻自身不能夠運(yùn)動(dòng)，在流動(dòng)或紊動(dòng)的環(huán)境下反而能獲得較大的增殖優(yōu)勢，這也是S0中藍(lán)藻能夠在第10d以后才消失而流動(dòng)水體中藍(lán)藻在第4d或第5d即消失的主要原因.S0中枝角類浮游動(dòng)物多以單細(xì)胞藻類為食，使得圍隔裝置內(nèi)以其難以覓食的纖維藻屬占據(jù)優(yōu)勢地位.

表2 流動(dòng)和靜止圍隔裝置內(nèi)浮游植物的優(yōu)勢種屬變化*

*運(yùn)行前為2010年8月10日中心湖(W2)優(yōu)勢藻種.

水動(dòng)力對浮游藻類群落演替的影響顯著.Leira等[16]指出水動(dòng)力強(qiáng)度變化是影響浮游藻類群落演替和數(shù)量變化的重要原因.Elliott等[17]認(rèn)為水動(dòng)力對藻類群落演替的影響極其顯著，水動(dòng)力條件會對一些浮游藻類生長造成抑制，同時(shí)促進(jìn)其他藻類的生長.本研究中兩組實(shí)驗(yàn)在流動(dòng)的圍隔裝置內(nèi)均發(fā)現(xiàn)了浮游植物優(yōu)勢藻種的顯著變化，藍(lán)藻的迅速消失、綠藻和硅藻形成生長優(yōu)勢是其主要特征.Huisman等[18]通過研究認(rèn)為，隨著水動(dòng)力作用的增強(qiáng)，水體Chl.a濃度增加受抑制，且易下沉的綠藻和硅藻替代能垂直自由運(yùn)動(dòng)的藍(lán)藻而成為優(yōu)勢種，認(rèn)為藻類對光的競爭是造成群落演替的主要原因.本研究中采用透明的有機(jī)玻璃水槽，底部水體也能得到較好的光照條件.與自然水體相比，可能排除了浮游藻類對光的競爭，具有偽空胞的藍(lán)藻無需浮到水面也能獲取足夠光照，因此對光的競爭可能不是引起本研究中藻類種群變化的原因.與所在的自然水體W2相比，W2中藍(lán)藻一直是優(yōu)勢種，流動(dòng)水槽S1和S2分別在第5d和4d以后藍(lán)藻消失，一方面說明流速對藍(lán)藻的抑制作用顯著，另一方面說明空間尺度的改變顯著影響了浮游藻類的生存環(huán)境，由于藻類生理習(xí)性的差異，不同物種對環(huán)境變化的適應(yīng)速度是不同的，因此造成了適應(yīng)相對較快的藻類生長速度加快，而適應(yīng)相對較慢的藻類生長受到抑制，甚至在這樣的環(huán)境中被淘汰.郭蔚華等[19]通過對優(yōu)勢藻種的分析認(rèn)為，在三峽水庫不同類型的水體中，優(yōu)勢硅藻、綠藻、藍(lán)藻對流速表現(xiàn)出多種結(jié)構(gòu)的適應(yīng)性，流速是影響優(yōu)勢硅藻、藍(lán)藻、綠藻組成變化的直接原因.

3 結(jié)論

1) 實(shí)驗(yàn)中0.06和0.10m/s的流速均對浮游藻類的生長具有明顯的抑制作用，總體抑制率可達(dá)50%左右.流速引起的DO濃度快速補(bǔ)充對浮游藻類的作用弱于實(shí)驗(yàn)中0.06和0.10m/s流速的作用，但實(shí)驗(yàn)初期流速加速營養(yǎng)鹽向藻細(xì)胞擴(kuò)散和實(shí)驗(yàn)后期流速對浮游動(dòng)物生長的抑制作用仍會對流速效應(yīng)造成較大影響.

2) 流速對浮游植物優(yōu)勢藻種變化影響顯著.與靜止水槽相比，流速對藍(lán)藻的生長造成明顯的抑制作用，同時(shí)促進(jìn)了綠藻的生長，這與國外眾多研究者在淡水河流和湖泊中的研究結(jié)果相一致.流速對浮游藻類種群變化的影響不僅表現(xiàn)為對光的競爭，而且與藻類對流速多種結(jié)構(gòu)的適應(yīng)性有關(guān).

致謝：本研究的野外工作得到同濟(jì)大學(xué)長江水環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室崇明水環(huán)境研究中心黃德發(fā)工程師和張?jiān)鰟俨┦康拇罅χС?在此表示感謝.

[1] 孔繁翔,馬榮華,高俊峰等.太湖藍(lán)藻水華的預(yù)防、預(yù)測和預(yù)警的理論與實(shí)踐.湖泊科學(xué),2009,21(3):314-328.

[2] 梁培瑜,王 烜,馬芳冰.水動(dòng)力條件對水體富營養(yǎng)化的影響.湖泊科學(xué),2013,25(4):455-462.

[3] 李大美.生態(tài)水力學(xué).北京:科學(xué)出版社,2006.

[4] 顏潤潤,逄 勇,趙 偉等.環(huán)流型水域水動(dòng)力對藻類生長的影響.中國環(huán)境科學(xué),2008,28(9):813-817.

[5] 黃鈺玲,劉德富,陳明曦.不同流速下水華生消的模擬.應(yīng)用生態(tài)學(xué)報(bào),2008,19(10):2293-2298.

[6] 王婷婷,朱 偉,李 林.不同溫度下水流對銅綠微囊藻生長的影響模擬.湖泊科學(xué),2010,22(4):563-568.

[7] Mitrovic SM, Oliver RL, Rees Cetal. Critical flow velocities for the growth and dominance ofAnabaenacircinalisin some turbid freshwater rivers.FreshwaterBiology, 2003,48(1): 164-174.

[8] 古濱河,劉正文,李寬意等譯.湖沼學(xué)——內(nèi)陸水生態(tài)系統(tǒng).北京：高等教育出版社,2011:26-32.

[9] Long TY, Wu L, Meng GHetal. Numerical simulation for impacts of hydrodynamic conditions on algae growth in Chongqing Section of Jialing River, China.EcologicalModelling, 2011,222(1): 112-119.

[10] 陳偉民,陳宇煒,秦伯強(qiáng)等.模擬水動(dòng)力對湖泊生物群落演替的實(shí)驗(yàn).湖泊科學(xué),2000,12(4):343-352.

[11] 朱宜平,張海平,李飛鵬等.水動(dòng)力對浮游生物影響的圍隔研究.環(huán)境科學(xué),2010,31(1):69-75.

[12] Li FP, Zhang HP, Zhu YPetal. Effect of flow velocity on phytoplankton biomass and composition in a freshwater lake.ScienceoftheTotalEnvironment, 2013,447:64-71.

[13] 李飛鵬,張海平,陳 玲.小型封閉水體環(huán)境因子與葉綠素a的時(shí)空分布及相關(guān)性研究.環(huán)境科學(xué),2013,34(10):3854-3861.

[14] Mitrovic SM, Chessman BC, Davie Aetal. Development of blooms ofCyclotellameneghinianaandNitzschiaspp. (Bacillariophyceae) in a shallow river and estimation of effective suppression flows.Hydrobiologia, 2008,596(1):173-185.

[15] Hondzo M, Warnaars TA. Coupled effects of small-scale turbulence and phytoplankton biomass in a small stratified lake.JournalofEnvironmentalEngineering-Asce, 2008, 134(12): 954-960.

[16] Leira M, Cantonati M. Effects of water-level fluctuations on lakes: an annotated bibliography.Hydrobiologia, 2008,613(1):171-184.

[17] Elliott JA, Defew L. Modelling the response of phytoplankton in a shallow lake (Loch Leven, UK) to changes in lake retention time and water temperature.Hydrobiologia, 2012,681(1):105-116.

[18] Huisman J, Sharples J, Stroom JMetal. Changes in turbulent mixing shift competition for light between phytoplankton species.Ecology, 2004,85(11): 2960-2970.

[19] 郭蔚華,李 楠,張 智等.嘉陵江出口段三類水體藍(lán)綠硅藻優(yōu)勢種變化機(jī)理.生態(tài)環(huán)境學(xué)報(bào),2009,18(1):51-56.

Simulation experiment on the effect of flow velocity on phytoplankton growth and composition

LI Feipeng1, GAO Ya2, ZHANG Haiping1, XIAO Yihua1& CHEN Ling1

(1：CollegeofEnvironmentalScienceandEngineering,TongjiUniversity,Shanghai200092,P.R.China)(2:TongjiArchitecturalDesign(Group)Co.,Ltd.,Shanghai200092,P.R.China)

An experimental simulation study was conducted using surface water from Zhongxinhu Lake of Chongming Island in oval plexiglass enclosures. The object of the study was to investigate the effect of flow velocity on phytoplankton growth and composition dynamics. Both flow conditions had inhibition effects on phytoplankton biomass in early experiment compared with the still enclosure and in late experiment the phytoplankton biomass stabilized in a certain level in all enclosures. Moreover, a higher Chl.a concentration was found in both flow enclosures in the late experiment, implying that numerous Cladocera zooplankton occurred in the still enclosure might be the reason that in flow conditions the turbulent flows inhibited the growth of zooplankton and reduced the predation pressure of phytoplankton. The turbulent flow also resulted in the changes of phytoplankton composition with the feature of the rapid disappearance of blue-green algae and the domination of green algae and diatoms. The result indicated that the continuous flow may directly caused the reducing species and composition dynamics of phytoplankton. The study can provide basis of the application of water division in eutrophication management and algae bloom prevention of rivers, lakes and reservoirs.

Flow velocity; phytoplankton; species dynamics; simulation experiment

*國家自然科學(xué)基金項(xiàng)目(51379146，51409190)和中國博士后科學(xué)基金項(xiàng)目(2013M531218)聯(lián)合資助.2014-01-24收稿；2014-04-15收修改稿.李飛鵬(1983～)，男，博士；E-mail:lifeipeng@#edu.cn.

**通信作者；E-mail:hpzhang@#edu.cn.