冉雅郡,楊瀟瀟,汪作煒,吳成陽,周 云,夏四清*

?

HRT對CO2-MBfR去除高SO42-廢水中NO3--N的影響

冉雅郡1,楊瀟瀟1,汪作煒2,吳成陽1,周 云1,夏四清1*

(1.同濟大學環(huán)境科學與工程學院,污染控制與資源化研究國家重點實驗室,上海 200092；2.深圳市南山區(qū)環(huán)境保護和水務局,廣東 深圳 518000)

針對含NO3--N與較高濃度SO42-實際工業(yè)廢水處理較難的問題,考察了不同水力停留時間(HRT)下連續(xù)運行的CO2-氫基質膜生物膜反應器(CO2-MBfR)處理模擬廢水和實際工業(yè)廢水的性能,結果表明,2種廢水的出水NO3--N濃度均隨著HRT的減小而增大,模擬廢水中NO3--N的處理效果和電子通量分配比例均優(yōu)于實際廢水,但其電子通量分配的格局基本不變:NO3--N和SO42-的電子通量分別在90.09%~97.49%和2.51%~9.91%左右.要實現實際廢水總氮達到15mg/L的排放標準,需維持HRT不少于10.4h.

二氧化碳-氫基質膜生物膜反應器；實際工業(yè)廢水；水力停留時間；硝酸鹽氮；硫酸鹽

近年來,由NO3--N、SO42-等氧化型污染物所帶來的水污染問題尤為嚴重[1-4].我國對企業(yè)排放的實際工業(yè)廢水中的總氮、硫化物要求日益嚴格,《紡織染整工業(yè)污染物排放標準》(GB4827-2012)[5]中規(guī)定,新建企業(yè)與現有企業(yè)統(tǒng)一執(zhí)行15mg/L的總氮排放標準,硫化物不得在出水中檢出.實際工業(yè)廢水因其COD高、可生化性差等問題,采用傳統(tǒng)處理工藝處理難度大、運行成本高[6-7].

雙膜曝氣膜生物膜反應器(CO2-MBfR)技術是一項清潔的膜處理技術,它是在傳統(tǒng)氫基質膜生物膜反應器(MBfR)的基礎上改進而來的.它的原理是H2和CO2從中空纖維膜內部向外以無泡曝氣的形式傳質擴散,自養(yǎng)菌利用CO2作為無機碳源、H2作為電子供體,以氧化性污染物為電子受體(底物)進行呼吸作用,將污染物還原成無毒或低毒形態(tài)[8-11].通過合理調節(jié)CO2分壓還能中和氧化態(tài)污染物降解過程中產生的OH-,使反應器內pH值控制在適宜氫自養(yǎng)還原菌生長的范圍內,同時避免高硬度實際廢水處理過程中硬度離子的沉淀[12].NO3--N和SO42-在MBfR中被還原成N2和S2-,其中還原NO3--N的菌屬主要為紅環(huán)菌屬、噬氫菌屬等[13],還原SO42-的菌屬主要為脫硫弧菌屬[14],它們以H2作為電子供體還原污染物,具有產生的剩余污泥量少、成本相對較低、清潔無毒、不會產生二次污染等優(yōu)點[1,15]. CO2-MBfR對地下水中的硝酸鹽、高氯酸鹽等污染物均有良好的去除效果[16],但對于實際工業(yè)廢水的處理情況鮮有報道.

本研究利用二氧化碳-氫基質生物膜反應器(CO2-MBfR)處理NO3--N和高濃度SO42-共存下的模擬廢水和實際工業(yè)廢水,探究CO2-MBfR去除高硫酸鹽模擬和實際工業(yè)廢水中的NO3--N實現最大的反硝化效果,并盡可能減少硫化物產生的可行性,考察水力停留時間(HRT)對其處理性能的影響,通過對比模擬和實際廢水NO3--N、SO42-的去除效果,探究實際廢水中兩者的競爭、降解規(guī)律,以期實現最佳HRT運行下的出水達標排放,為CO2-MBfR處理實際工業(yè)廢水提供理論基礎.

1 材料與方法

1.1 實際工業(yè)廢水和模擬廢水

實際工業(yè)廢水來自浙江某公司,其負責收集全縣的生活污水與工業(yè)廢水進行處理,其中工業(yè)廢水以印染廢水為主,以工業(yè)廢水(85%)和生活污水(15%)的比例混合,混合后水體污染物指標主要為COD、氨氮兩種.廢水經過公司內部工藝處理后的出水水質指標如表1所示.出水總氮主要以硝氮形式存在,COD、硫酸根含量略高,水體呈中性,Ca、Mg、K、Na等營養(yǎng)元素豐富,適于微生物生長.

表1 實際工業(yè)廢水成分

表2 模擬廢水培養(yǎng)基成分

模擬廢水以NaHCO3為無機碳源,KH2PO4和Na2HPO4為緩沖介質,并添加微生物生長必需的微量元素,具體參數如表2所示,所有藥品均為分析純.使用培養(yǎng)基前須向水體中曝氮氣3~5min,用以去除進水中的溶解氧.

1.2 氫自養(yǎng)還原菌的馴化

啟動反應器前,需用搖瓶實驗馴化出反應器的種泥.向250mL氣密性塑料血清瓶注入200mL泥水混合液,頂空50mL通入氫氣.泥水混合液來自上海市曲陽污水處理廠的缺氧池污泥,可揮發(fā)性固體(VSS)和總懸浮固體(TSS)分別為(2.13±0.32),(2.68± 0.18)g/L,在打入血清瓶前分別經過生理鹽水清洗1次,磷酸鹽緩沖鹽(KH2PO4+Na2HPO4+KCl)洗滌2次并濃縮,使用無菌注射器分批次打入生理鹽水瓶,維持TSS在3g/L.初始NO3--N濃度為10mg/L,SO42-濃度為30mg/L,待NO3--N濃度至檢測限以下,SO42-還原趨于穩(wěn)定,即可認為馴化完成.

1.3 CO2-MBfR的啟動掛膜與運行

試驗采用圖1所示的CO2-MBfR裝置,反應器的特征參數如表3所示,反應器內膜絲采用PVC(聚氯乙烯)中空纖維膜.利用無菌注射器分別將50mL馴化完成的泥水混合液打入2個完全相同的反應器,用于接種掛膜.當反應器出水穩(wěn)定,且膜絲表面基本覆蓋一層微生物則認為掛膜成功.

圖1 CO2-MBfR反應器及其去除NO3--N過程示意

表3 CO2-MBfR的特征參數

這2個反應器分別進水人工模擬廢水及實際廢水,同時啟動運行,人工模擬廢水中投加 30mg/L NO3--N、1400mg/L SO42-2種氧化態(tài)污染物及必需營養(yǎng)元素,2個完全相同的反應器的運行參數為:pH 7.5、室溫30℃、氫分壓0.06MPa.待反應器掛膜成功并穩(wěn)定運行110d后,保持進水和氫分壓不變,使2個反應器分別按表4的試驗安排改變工況,工況一、二、三的HRT分別為30.12,10.4,5.2h,考察反應器出水水質指標的變化情況.

表4 CO2-MBfR 還原NO3--N研究的運行參數

1.4 取樣與分析

調整反應器水力停留時間,每個工況持續(xù)時間均多于5倍停留時間.每切換到一個新工況后,讓反應器穩(wěn)定24h(多于4倍水力停留時間),可認為其達到擬穩(wěn)態(tài)[15,17],再開始取樣.采用無菌注射器進行定時取樣,經0.45μm水相濾膜過濾后保存在4℃冰箱中.NO3--N、NO2--N、TN通過分光光度法測定,SO42-采用美國戴安公司ICS-1000型離子色譜儀(AS20離子色譜柱)進行測定,pH值由美國哈希公司的 HQ40d型pH計測定.

1.5 基本參數

NO3--N、SO42-的去除效果用去除率、去除通量和當量電子通量表示.其中去除通量(,g/m3×d)指單位膜面積在單位時間內還原氧化態(tài)污染物的質量,計算公式為

式中:為進水流量(m3/d),ie分別表示污染物進水濃度(g/L)、出水濃度(g/L),為膜面積(m2).而當量電子通量(e,mol e-/m2×d)是指單位時間內單位面積的生物膜內某種氧化態(tài)污染物被還原后接受的電子摩爾數,計算公式為

e=/(2)

式中:為轉換因子,NO3--N、SO42-的轉換因子分別為2.8,12g/mol e-.

2 結果和討論

2.1 HRT對CO2-MBfR去除模擬廢水中NO3--N的影響

如圖2所示,隨著水力停留時間的逐漸縮短,出水NO3--N、SO42-濃度呈明顯升高趨勢,NO3--N出水濃度分別穩(wěn)定在0.43,2.12,4.14mg/L,模擬配水中總氮為硝氮與亞硝氮之和,因反硝化較為完全,亞硝氮積累不明顯,總氮出水濃度近似等同于硝氮出水濃度;出水SO42-濃度不斷升高,出水SO42-濃度與進水SO42-濃度之間差距隨HRT的縮短而減小,當HRT為5.2h時,出水SO42-濃度幾乎等于進水SO42-濃度,說明此時硫酸鹽幾乎不被還原.可知,當HRT不斷縮短,NO3--N去除效果下降,可能是因為隨著反應器中處理的污染物負荷不斷增大,反硝化逐漸不完全出現NO2--N的積累[18],此時MBfR無法完成對硝酸鹽的徹底反硝化,硝酸鹽進水負荷超過MBfR可以承受的反硝化負荷.

圖2 不同HRT下模擬廢水中各污染物濃度變化

如圖3所示,在HRT分別為30.12,10.4,5.2h條件下,硝酸鹽的去除率及還原量出現下降,NO3--N去除率分別為98.43%、92.53%和86.66%,NO3--N還原通量分別為0.49,0.47,0.422g/(m2·d).Xia等[18]研究MBfR處理硝酸鹽時發(fā)現,HRT從4h降到0.5h時,總氮去除率明顯下降.Aydin等[19]研究MBfR同時去除硝酸鹽和金霉素,當氫氣壓力為0.01MPa,HRT從15h降到1h時總氮去除效率由100%降到了52%,但當氫氣壓力提高到0.04MPa時,HRT的縮短對硝酸鹽的還原影響并不大,降解率始終大于92%.在本實驗中HRT由30.12h縮短至5.2h,NO3--N的去除率均在85%以上,出水NO3--N濃度均低于5mg/L,說明0.06MPa氫氣壓力足以還原硝酸鹽.

圖3 不同HRT條件下處理模擬廢水各污染物去除率及去除通量

硫酸鹽在本實驗中還原率低,隨著HRT縮短,硫酸鹽的去除率及還原量亦出現下降,SO42-去除率分別為0.8%、0.2%、0.01%,SO42-還原通量分別為0.21,0.1,0.04g/(m2·d).當HRT降到5.2h時,硫酸鹽幾乎不被還原.Chung等[20]在研究MBfR處理鉻酸鹽時也發(fā)現,在鉻酸鹽、硝酸鹽還原效果都較好的情況下,硫酸鹽仍然幾乎沒有去除.究其原因,Ontiveros- Valencia等[2]在研究MBfR中硝酸鹽、硫酸鹽互相影響的機制中的結果表明,當氫氣供應不足時硝酸鹽優(yōu)先被氫氣利用,在氫氣壓力提高到0.2MPa時,硝酸鹽完全被還原,而硫酸鹽只有在氫氣壓力大于0.34MPa時才會開始被還原;當氫氣供應充足時,硫酸鹽的還原還受硝酸鹽負荷的影響,只有硝酸鹽氮的表面負荷低于0.13g N/(m2×d)時硫酸鹽才會被還原.在本實驗中硝酸鹽還原率高而硫酸鹽還原率低,說明0.06MPa氫氣分壓和硝酸鹽的初始濃度利于實現高濃度硫酸鹽廢水中硝酸鹽最大的反硝化而硫酸鹽盡可能不被還原.

2.2 HRT對CO2-MBfR去除實際廢水中NO3--N的影響

同樣維持氫分壓為0.06MPa,改變HRT分別為30.12,10.4,5.2h,考察MBfR中HRT對實際廢水中氧化態(tài)污染物還原的影響及出水總氮達標所需最優(yōu)工況.如圖4所示,隨著水力停留時間的逐漸縮短,出水 NO3--N、SO42-濃度呈明顯升高趨勢,NO3--N出水濃度分別穩(wěn)定在4.53,8.73,15.62mg/L,總氮出水濃度分別穩(wěn)定在13.21,17.59,24.34mg/L,出水SO42-濃度與進水SO42-濃度之間差距隨HRT的縮短而減小,當HRT為5.2h時,出水SO42-濃度幾乎等于進水SO42-濃度.

圖4 不同HRT條件處理實際廢水污染物濃度變化

如圖5所示,在HRT分別為30.12,10.4,5.2h條件時,硝酸鹽的去除率及還原量出現明顯下降, NO3--N去除率分別為90.44%、74.02%和55.29%, NO3--N還原通量分別為0.48,0.41,0.31g/(m2·d).隨著HRT縮短,SO42-去除率及還原量亦出現明顯下降,SO42-去除率分別為1.06%、0.4%、0.1%,SO42-還原通量分別為0.24,0.09,0.02g/(m2·d).與模擬廢水的去除效果相比,同等HRT的條件下,實際廢水硝態(tài)氮的去除效果更差,由前面的分析可知,模擬廢水硝態(tài)氮的去除效果好說明氫氣壓力供應充足,進而說明實際廢水的處理效果受水質本身的影響更大.其中可能的原因較多,Zhai等[21]研究生物膜反應器中硝酸鹽的去除,結果表明硝酸鹽的還原受水的鹽度和COD的影響,鹽度越高、COD值越大,出水硝酸鹽氮的濃度越高,即硝酸鹽處理效果變差.Van Ginkel等[22]研究MBfR處理離子交換鹵水中的硝酸鹽和高氯酸鹽,表明鹽度對硝酸鹽的影響較大,隨著進水鹽度的升高,硝酸鹽的去除率明顯降低.而在本實驗中,實際廢水的COD和鹽度都比模擬廢水高,所以反硝化效果更差可能是由這2個原因導致.想要使實際廢水達標排放并盡可能縮短HRT,在進反應器之前需對實際廢水進行適當的前處理.

圖5 不同HRT條件處理實際廢水各污染物去除率及去除通量

MBfR反應器內的氫氣是通過膜內外的氫氣壓力差為傳質動力,在本實驗中實際廢水硝態(tài)氮去除率明顯低于模擬廢水,說明此時氫氣利用率降低,氫氣的實際通量也會減少[23].

雖然實際廢水的反硝化效果更差,但與模擬廢水的出水水質變化趨勢相同,HRT縮短都會降低兩者的反硝化效果.當HRT為30.12h時,實際廢水的出水總氮低于15mg/L,要利用CO2-MBfR處理工業(yè)廢水實現較好的硝氮去除需控制反應器中水體水力停留時間長于10.4h,且水力停留時間越久,總氮去除效果越佳.當水力停留時間過短時,反應器中污染物負荷會超過其承受負荷,導致出水水質較差.

2.3 不同HRT下CO2-MBfR運行當量電子通量分析

為了更準確地分析同一反應體系中各種不同的污染物(電子受體)對氫氣(電子供體)的競爭關系,計算得出在各HRT下電子在污染物之間的分配情況如表5所示.電子通量分配格局基本不變,NO3--N的當量電子通量分配最大,其次是SO42-.實際廢水中NO3--N的當量電子通量分配(90.09%~97.49%), SO42-的當量電子通量(2.51%~9.91%);模擬廢水中NO3--N當量電子通量分配(90.86%~99.63%),SO42-的當量電子通量(0.37%~9.14%).說明反硝化過程消耗的電子供體數量最多,在同一體系中SO42-對電子供體的競爭性搶奪不占優(yōu)勢,該結論與Ontiveros- Valencia等[2]的研究結果一致.當氫氣供應一定,隨著HRT不斷縮短,反應器中污染物負荷不斷提高,各污染物去除效果逐漸變差,NO3--N大量的爭奪電子供體使得SO42-的還原處于劣勢,對SO42-還原有較大影響,NO3--N仍是反應體系中第一大電子供體消耗者.MBfR中硝酸鹽、硫酸鹽與其他污染物共存時的去除情況也得出類似的結論:硝酸鹽是第一優(yōu)先電子受體,硫酸鹽是第二優(yōu)先電子受體,這兩者是主要的電子供體消耗者[20,24-26].

表5 不同HRT條件下廢水中污染物的電子通量以及分配情況

3 結論

3.1 經80d左右運行的CO2-MBfR處理模擬廢水和實際廢水實驗中,隨著HRT的減小,兩種廢水的出水NO3--N濃度逐步增大,去除率降低,但是模擬廢水中NO3--N的去除效果明顯優(yōu)于實際廢水,實際廢水硝態(tài)氮去除效果差不是因為氫氣分壓供應不足,而可能是受鹽度、COD等水質指標的影響.

3.2 模擬廢水電子通量分配比例優(yōu)于實際廢水,但電子通量分配的格局基本不變:NO3--N和 SO42-的電子通量分別在90.09%~97.49%和2.51%~9.91%左右.H2供應一定時,隨著HRT不斷縮短,反應器中污染物負荷不斷提高,各污染物去除效果逐漸變差,NO3--N大量的爭奪電子供體使得SO42-的還原處于劣勢,對SO42-還原有較大影響,NO3--N仍是反應體系中第一大電子供體消耗者.要實現實際廢水總氮達到15mg/L的排放標準,需維持HRT不少于10.4h.

[1] Nerenberg R, Rittmann B E. Hydrogen-based, hollow-fiber membrane biofilm reactor for reduction of perchlorate and other oxidized contaminants [J]. Water Science and Technology, 2004, 49(11/12):223-230.

[2] Ontiveros-Valencia A, Ziv-El M, Zhao H P, et al. Interactions between nitrate-reducing and sulfate-reducing bacteria coexisting in a hydrogen-fed biofilm [J]. Environmental Science and Technology, 2012,46(20):11289-11298.

[3] Li H, Lin H, Xu X, et al. Simultaneous Bioreduction of Multiple Oxidized Contaminants Using a Membrane Biofilm Reactor [J]. Water Environment Research, 2017,89(2):178-185.

[4] Ginkel S W V, Zhou C, Lien M, et al. Hydrogen-Based Nitrate and Selenate Bioreductions in Flue-Gas Desulfurization Brine [J]. Journal of Environmental Engineering, 2011,137(1):63-68.

[5] GB 4287-2012 紡織染整工業(yè)水污染物排放標準 [S].

[6] 魯 璐,劉漢湖. Fenton試劑預處理實際印染廢水的實驗研究 [J]. 環(huán)境科學與管理, 2008,33(3):89-92.

[7] 彭趙旭,彭永臻,劉旭亮,等.實際污水與模擬污水活性污泥系統(tǒng)的特性差異 [J]. 環(huán)境工程學報, 2012,6(3):749-754.

[8] 楊 昕,夏四清,李海翔,等.氫自養(yǎng)還原菌去除地下水中硫酸鹽的影響因素研究 [J]. 水處理技術, 2010,36(5):38-41.

[9] 夏四清,鐘佛華,張彥浩.氫自養(yǎng)反硝化去除水中硝酸鹽的影響因素研究 [J]. 中國給水排水, 2008,24(21):5-8.

[10] Chung J, Rittmann B E, Wright W F, et al. Simultaneous bio-reduction of nitrate, perchlorate, selenate, chromate, arsenate, and dibromochloropropane using a hydrogen-based membrane biofilm reactor [J]. Biodegradation, 2007,18(2):199-209.

[11] Xia S, Xu X, Zhou C, et al. Direct delivery of CO2into a hydrogen- based membrane biofilm reactor and model development [J]. Chemical Engineering Journal, 2016,290:154-160.

[12] 夏四清,徐曉茵,王晨輝.利用CO2為碳源的新型氫基質生物膜技術去除水中硝酸鹽 [J]. 環(huán)境化學, 2016,35(5):1050-1057.

[13] Zhang Y, Zhong F, Xia S, et al. Autohydrogenotrophic denitrification of drinking water using a polyvinyl chloride hollow fiber membrane biofilm reactor [J]. Journal of Hazardous Materials, 2009,170(1):203- 209.

[14] Zhao H P, Ilhan Z E, Ontiveros-Valencia A, et al. Effects of multiple electron acceptors on microbial interactions in a hydrogen-based biofilm [J]. Environmental Science and Technology, 2013,47(13): 7396-7403.

[15] Rittmann B E, Nerenberg R, Lee K C, et al. The hydrogen-based hollow-fiber membrane biofilm reactor (HFMBfR) for removing oxidized contaminants [J]. Water Science and Technology Water Supply, 2004,49(49):223-230.

[16] 楊瀟瀟,汪作煒,夏四清.雙膜曝氣生物膜反應器除水中硝氮和高氯酸鹽 [J]. 中國環(huán)境科學, 2016,36(10):2972-2980.

[17] Chung J, Li X, Rittmann B E. Bio-reduction of arsenate using a hydrogen-based membrane biofilm reactor [J]. Chemosphere, 2006, 65(1):24-34.

[18] Xia S, Zhong F, Zhang Y, et al. Bio-reduction of nitrate from groundwater using a hydrogen-based membrane biofilm reactor [J]. Journal of Environmental Sciences, 2010,22(2):257-262.

[19] Aydin E, Sahin M, Taskan E, et al. Chlortetracycline removal by using hydrogen based membrane biofilm reactor [J]. Journal of Hazardous Materials, 2016,320:88-95.

[20] Chung J, Nerenberg R, Rittmann B E. Bio-reduction of soluble chromate using a hydrogen-based membrane biofilm reactor [J]. Water Research, 2006,40(8):1634-1642.

[21] Zhai S, Ji M, Zhao Y, et al. Effects of salinity and COD/N on denitrification and bacterial community in dicyclic-type electrode based biofilm reactor [J]. Chemosphere, 2018,192:328.

[22] Van Ginkel S W, Ahn C H, Badruzzaman M, et al. Kinetics of nitrate and perchlorate reduction in ion-exchange brine using the membrane biofilm reactor (MBfR) [J]. Water Research, 2008,42(15):4197-4205.

[23] Tang Y, Zhou C, Ginkel S W V, et al. Hydrogen permeability of the hollow fibers used in H 2 -based membrane biofilm reactors [J]. Journal of Membrane Science, 2012,407–408:176-183.

[24] Chung J, Rittmann B E. Bio-reductive dechlorination of 1,1,1- trichloroethane and chloroform using a hydrogen-based membrane biofilm reactor [J]. Biotechnology Bioengineering, 2007,97(1):52-60.

[25] Xia S, Zhang Z, Zhong F, et al. High efficiency removal of 2-chlorophenol from drinking water by a hydrogen-based polyvinyl chloride membrane biofilm reactor [J]. Journal of Hazardous Materials, 2011,186(2/3):1367-1373.

[26] Chung J, Nerenberg R, Rittmann B E. Bioreduction of selenate using a hydrogen-based membrane biofilm reactor [J]. Environmental Science and Technology, 2006,40(5):1664-1671.

Effects of HRT on the removal of NO3--N from high SO42-wastewater using CO2-MBfR.

RAN Ya-jun1, YANG Xiao-xiao1, WANG Zuo-wei2, WU Cheng-yang1, ZHOU Yun1,XIA Si-qing1*

(1.College of Environmental Science and Engineering, Tongji University, Shanghai 200092, China；2.Environmental Protection and Water Administration Bureau of Shenzhen Nanshan District, Shenzhen 518000, China)., 2018,38(9)：3322~3327

Due to the difficulty of NO3--N removal during treatment of high-sulfate (SO42-) actual industrial wastewater (AIW) treatment, we tested the effects of hydraulic retention time (HRT) on treatment performance of simulated wastewater and AIW using CO2-hydrogen based membrane biofilm reactor (CO2-MBfR). Results show that the effluent concentrations of NO3--N from two kinds of wastewater gradually increased with decreased HRT, and the removal efficiency of NO3--N and electron flux distribution ratio (EFDR) in simulated wastewater treatment were better than that of in AIW treatment, but the electron flux distribution pattern was nearly the same: the EFDR of NO3--N and SO42-were 90.09%~97.49% and 2.51%~9.91%, respectively. Overall, our results reveal that CO2-MBfR could high efficiently remove NO3--N of AIW. From a practical perspective, HRT higher than 10.4h is suggested for achieving the total nitrogen emission standard of 15mg/L in the actual wastewater treatment.

carbon dioxide-hydrogen based membrane biofilm reactor (CO2-MBfR)；actual industrial wastewater；hydraulic retention time (HRT)；NO3--N；SO42-

X703

A

1000-6923(2018)09-3322-06

冉雅郡(1995-),女,重慶墊江人,同濟大學碩士研究生,研究方向為水污染控制及資源化研究.發(fā)表論文1篇.

2018-02-12

國家自然科學基金資助項目(51678422)

*責任作者, 教授, siqingxia@#edu.cn