馮顯露,劉新穎,孫德智,劉潔琦,霍 達,馬嘉苑,張 真,黨 巖

?

投加顆粒活性炭強化餐廚垃圾的厭氧處理

馮顯露,劉新穎,孫德智,劉潔琦,霍 達,馬嘉苑,張 真,黨 巖*

(北京林業(yè)大學(xué)環(huán)境科學(xué)與工程學(xué)院,水體污染源控制技術(shù)北京市重點實驗室,北京市污染水體源控與生態(tài)修復(fù)技術(shù)高等學(xué)院工程研究中心,北京 100083)

通過投加顆?；钚蕴?GAC)強化直接種間電子傳遞(DIET)進而提升餐廚垃圾的厭氧產(chǎn)甲烷處理效能,并研究了GAC投加導(dǎo)致的微生物群落變化.研究發(fā)現(xiàn),投加了GAC的實驗組反應(yīng)器能夠在更高的有機負荷下(10.4kgCOD/(m3·d))穩(wěn)定運行并維持較高的甲烷產(chǎn)率,不投加GAC的對照組在有機負荷7.8kgCOD/(m3·d)時甲烷產(chǎn)率及pH值均明顯降低,揮發(fā)酸大量積累,反應(yīng)器酸化崩潰.微生物群落結(jié)構(gòu)分析發(fā)現(xiàn),GAC表面富集了大量可以胞外電子傳遞的細菌(占細菌豐度的34%)和可以參與DIET的產(chǎn)甲烷菌(占古菌豐度的88%),表明GAC的加入可以有效富集這兩類微生物的生長,并可能通過GAC強化DIET促進了餐廚垃圾的厭氧消化.

直接種間電子傳遞；厭氧產(chǎn)甲烷；餐廚垃圾；顆?；钚蕴?/p>

餐廚垃圾是城市生活垃圾的主要組成部分,全世界每年產(chǎn)出餐廚垃圾達13億t以上[1].目前,我國處理餐廚垃圾主要通過垃圾填埋和焚燒[2].隨著我國城市生活垃圾分類舉措的大力開展,具有高含水率、高熱值、易腐爛、營養(yǎng)豐富等特點的餐廚垃圾更適合采用厭氧消化手段來單獨處理[3],不但可以實現(xiàn)能源的回收,還可以有效緩解我國垃圾填埋場地緊張、垃圾滲瀝液產(chǎn)量巨大的問題[4].然而,餐廚垃圾厭氧消化的處理效率一直是制約其發(fā)展的一個重要原因[5].

厭氧產(chǎn)甲烷體系中的直接種間電子傳遞(DIET)是近年發(fā)現(xiàn)的不同于傳統(tǒng)種間氫傳遞的全新的產(chǎn)甲烷過程[6].DIET是具有胞外電子傳遞功能的微生物(例如等)將降解有機物產(chǎn)生的電子通過導(dǎo)電菌毛(e-pili)的電子傳遞能力直接傳遞給產(chǎn)甲烷菌,產(chǎn)甲烷菌接受電子并還原CO2的過程.近期的研究發(fā)現(xiàn),DIET可以通過向產(chǎn)甲烷體系中投加導(dǎo)電材料強化DIET過程[7].

針對餐廚垃圾厭氧消化的處理效率問題,本文通過在厭氧消化反應(yīng)器中投加顆?；钚蕴?GAC)來強化餐廚垃圾的處理,考察了反應(yīng)器穩(wěn)定運行所能承受的最高有機負荷(OLR),并對GAC表面、懸浮污泥的微生物群落進行了分析,討論了GAC強化餐廚垃圾厭氧消化的原因.

1 材料與方法

1.1 反應(yīng)器與接種污泥

實驗室反應(yīng)器裝置如圖1所示.

圖1 反應(yīng)器示意

設(shè)置2組每組3個共6個反應(yīng)器,一組額外投加100g顆?；钚蕴?實驗組,顆?；钚蕴砍屎谏?、粒徑為8~20目、比表面積約為1000m3/g、電導(dǎo)率約為26μS/cm),一組不投加材料(對照組).實驗所用反應(yīng)器為有效容積2L的玻璃瓶,其上由橡膠塞密封.每個反應(yīng)器上有一個加(取)樣口和一個出氣口,出氣口連接一個10L氣袋.所有反應(yīng)器在37℃、100rpm轉(zhuǎn)速、避光條件下運行.接種污泥為河北某處理啤酒廢水的厭氧懸浮污泥,TSS為12,000mg/L,VSS/ TSS為0.70.

1.2 餐廚垃圾及反應(yīng)器運行

選用商品狗糧(Vita-one,Nihon,東京,以餐廚垃圾制成)作為處理對象,通過表1的配方配制成餐廚垃圾漿體加入反應(yīng)器內(nèi)進行厭氧消化.

表1 餐廚垃圾漿體成分

注: 以每升含量參考,未標明數(shù)值單位均為mg. DL維生素,DL微量元素具體配方請參考Zhao等[8]的報道.

啟動初期向兩組反應(yīng)器內(nèi)裝入500mL接種污泥,400mL去離子水,100mL配置好的餐廚垃圾漿體.進水OLR逐漸提高,進樣中狗糧混合物濃度為10,20,40,60,80g/L,對應(yīng)OLR分別為1.3, 2.6,5.2,7.8,10.4kgCOD/(m3·d).以24h為一個運行周期,每個周期的前23.5h反應(yīng)器厭氧消化,剩余的0.5h從每個反應(yīng)器內(nèi)排出100mL消化液,再注入100mL新配置的餐廚垃圾漿體,進入下一個運行周期,保證每個OLR下水力停留時間為10d.

1.3 分析項目與方法

pH值采用Denver UB-10pH測定儀;揮發(fā)性脂肪酸(VFAs)的測定采用高效液相色譜儀(Waters e2695,美國);COD利用重鉻酸鉀法測定;氣體體積利用氣體流量計(FMA4000,Omega,美國)測定;氣體用配有氫火焰檢測器(GC-8A)的氣相色譜儀(Agilent 7890A)測定.

1.4 DNA提取與高通量測序

取泥樣1.0mL或顆?；钚蕴?.0g,樣品立即放入液氮中凍結(jié)并儲存于-80℃條件下.樣品提取前利用研缽將樣品破碎,DNA的提取根據(jù)DNA快速提取試劑盒說明書(艾德萊,北京),提取后DNA經(jīng)微量紫外分光光度計(NanoDrop)測定純度,并委托北京奧維森生物技術(shù)公司測序.通過引物27F和394F對其細菌和古菌進行16S rRNA基因擴增[9].高通量測序的原始序列已上傳至NCBI數(shù)據(jù)庫,序列號為SAMN07344859、SAMN07344860.

1.5 數(shù)據(jù)處理

平行實驗均設(shè)置3個平行,利用Micosoft Excel 2013軟件計算平均值(average)和標準偏差(standard deviation),各種指標以平均值±標準偏差在圖中以誤差棒的形式表示.

2 結(jié)果與討論

2.1 反應(yīng)器運行

由圖2可知,兩組厭氧消化反應(yīng)器同步運行, OLR由1.3kgCOD/(m3·d)起逐步提升,投加了GAC的實驗組最終提升至10.4kgCOD/ (m3·d),而對照組則提升至7.8kgCOD/(m3·d),反應(yīng)器崩潰時停止運行.整個運行過程中HRT始終保持10d.

兩組反應(yīng)器在OLR為1.3kgCOD/(m3·d)和2.6kgCOD/(m3·d)時(1~14d)均高效穩(wěn)定運行,pH值維持在7.0~7.5之間,每個運行周期末反應(yīng)體系內(nèi)VFAs含量接近0,剩余COD濃度均低于200mg/L,表明餐廚垃圾厭氧消化產(chǎn)生的VFAs被及時充分的降解.并且兩組反應(yīng)器在這一階段甲烷產(chǎn)率基本一致.然而,當反應(yīng)器OLR進一步提升至5.2kgCOD/(m3·d)時(15~24d),雖然兩組反應(yīng)器仍能穩(wěn)定運行,且甲烷轉(zhuǎn)化率均在86%以上(以實際產(chǎn)甲烷折算成COD計[10]),但是實驗組的甲烷產(chǎn)率呈現(xiàn)出高于對照組的趨勢,對照組的VFAs濃度有上升趨勢,而實驗組仍接近0,實驗組剩余COD濃度均低于200mg/L,而對照組的剩余COD濃度上升至約560mg/L.

進一步提升OLR提升至7.8kgCOD/(m3·d) (25~34d),實驗組反應(yīng)器仍能夠穩(wěn)定運行,運行周期末反應(yīng)體系內(nèi)剩余COD濃度仍低于380mg/L,甲烷轉(zhuǎn)化率維持在80%以上(以實際產(chǎn)甲烷折算成COD計),pH值穩(wěn)定在7.3以上,雖然乙酸、丙酸的濃度有所升高,但明顯低于對照組;相比較而言,在這一負荷下,對照組的甲烷產(chǎn)率直線下降,在短短5d時間內(nèi)由>2000mL/d降低至不足100mL/d,pH值也顯著下降至5.3左右,運行周期末反應(yīng)器內(nèi)剩余COD仍高達5200mg/L以上,乙酸、丙酸、丁酸均大量積累,總量超過了100mmol/L,反應(yīng)器酸化崩潰.此時停止了對照組反應(yīng)器的運行.

圖2 兩組反應(yīng)器的運行情況

(A)甲烷產(chǎn)率;(B)pH值;(C)出水揮發(fā)酸含量

繼續(xù)提升實驗組反應(yīng)器的OLR至10.4kgCOD/(m3·d)(35~44d),此時實驗組pH值維持在7.0以上,運行周期末反應(yīng)器內(nèi)剩余COD濃度上升至3200mg/L,VFAs積累濃度相比于上一負荷提升至約60mmol/L,以丙酸和乙酸為主,丁酸仍維持在極低的水平,且VFAs總量低于對照組在上一負荷酸化崩潰時的積累濃度.此時甲烷產(chǎn)率有所降低,甲烷轉(zhuǎn)化率也降低至65%(以實際產(chǎn)甲烷折算成COD計),但反應(yīng)器仍能在較低處理效率的情況下穩(wěn)定運行.以上實驗結(jié)果表明在厭氧消化反應(yīng)器內(nèi)投加GAC可以顯著提升餐廚垃圾的處理效能,穩(wěn)定運行的最高OLR高達7.8kgCOD/(m3·d),相比于對照組提升了33%以上,且明顯高于其他餐廚垃圾厭氧消化的相關(guān)報道[4,11-12].

2.2 微生物群落分析

通過高通量測序技術(shù)對實驗組的GAC表面、懸浮污泥,對照組的懸浮污泥進行微生物群落結(jié)構(gòu)分析(圖3).對于細菌,3個樣品的群落結(jié)構(gòu)在門的層級上類似,主要分布于厚壁菌門(Firmicutes)、擬桿菌門(Bacteroidetes)、綠彎菌門(Chloroflexi)、互養(yǎng)菌門(Synergistetes)等,其中厚壁菌門的細菌占絕對優(yōu)勢地位,在3個樣品中的豐度均超過68%(圖3A).對細菌群落結(jié)構(gòu)進一步深入到菌科和菌屬的層級進行分析,3個樣品的群落結(jié)構(gòu)差異變得明晰:大量具有胞外電子傳遞能力的細菌在GAC表面和GAC懸浮污泥中富集生長,總豐度分別達到了34.1%和26.7%,而在沒有投加GAC的對照組反應(yīng)器內(nèi),具有胞外電子傳遞能力的細菌的總豐度僅為7.5%(圖3B).具有胞外電子傳遞能力的微生物往往能將降解有機物產(chǎn)生的電子傳遞到細胞外的電子受體,如Fe(III)、電極和產(chǎn)甲烷菌等,使得它們具有潛在的“產(chǎn)電”能力[13].例如,屬于厚壁菌門的、和3個菌屬的菌在GAC顆粒表面(豐度分別為4.4%、3.7%和9.9%)和GAC懸浮污泥中(豐度分別為12.3%、5.5%、5.0%)均富集生長,而在對照組的懸浮污泥中的豐度均低于1%.

圖3 實驗組反應(yīng)器中GAC表面、GAC懸浮污泥和對照組懸浮污泥的微生物群落結(jié)構(gòu)的對比

(A)細菌;(B)具有胞外電子傳遞能力的細菌;(C)古菌

Kim等[13]發(fā)現(xiàn)可以利用多種可發(fā)酵有機物作為電子供體,將降解這些有機物產(chǎn)生的電子通過胞外電子傳遞的形式傳遞給Fe(III);Hernandez-Eugenio等[14]發(fā)現(xiàn)可以在硫單質(zhì)做電子受體時將氨基酸、糖類等有機物降解產(chǎn)生的電子傳遞給硫單質(zhì)并生成H2S,Lovley等[15]也證實微生物能將電子傳遞給硫單質(zhì),通常也能傳遞給Fe(III)或電極;Cord-Ruwisch和Garcia[16]研究發(fā)現(xiàn)同樣具有將有機物降解產(chǎn)生的電子傳遞給硫單質(zhì)的能力.

此外,屬于、-和菌科的菌在GAC表面專性附著生長,其豐度分別為1.2%、6.0%和3.7%,而在GAC懸浮污泥和空白懸浮污泥中均含量極低(<0.3%)或未檢出(圖3B).已報道著名的產(chǎn)電細菌[17]和[18]分別屬于和菌科;而,具有胞外電子傳遞能力[19],則是屬于菌科的代表菌屬.

上述對細菌群落結(jié)構(gòu)的分析表明,厭氧消化反應(yīng)器在投加了GAC后有利于富集這些具有潛在“產(chǎn)電”能力的細菌的生長,而具有導(dǎo)電性的GAC則可能作為電子傳遞載體將這些細菌降解餐廚垃圾過程中產(chǎn)生的電子傳遞給該反應(yīng)器體系中的電子受體——產(chǎn)甲烷菌.截止到目前,通過實驗明確證實能夠通過DIET的形式直接接受電子還原CO2產(chǎn)甲烷的產(chǎn)甲烷菌只有2類:甲烷髦毛菌屬()和甲烷八疊球菌屬()[6,20].從古菌的群落結(jié)構(gòu)可以看出,GAC表面的豐度(74.5%)明顯高于GAC懸浮污泥(38.2%)和對照組懸浮污泥(42.0%),而的含量差異更為明顯,在GAC表面的豐度為14.3%,而在GAC懸浮污泥和對照組懸浮污泥樣品中的豐度僅為0.8%和1.1%(圖3C).

投加GAC的厭氧消化反應(yīng)器處理餐廚垃圾的效能明顯提升,在OLR為10.4kgCOD/(m3·d)的條件下仍能穩(wěn)定運行,而對照組在OLR為7.8kgCOD/(m3·d)時已經(jīng)酸化崩潰.餐廚垃圾易于水解發(fā)酵,中間代謝產(chǎn)物VFAs產(chǎn)生速率較快,其中丙酸等小分子VFAs如不能及時進一步降解去處會在反應(yīng)器中積累,對微生物活性和生長產(chǎn)生抑制作用[21].已有研究表明,在生物反應(yīng)器中投加導(dǎo)電碳材料(如生物炭)或磁鐵礦可以通過強化DIET加速小分子VFAs比如丙酸、丁酸的降解[8,22].本文投加GAC后,厭氧消化反應(yīng)器中這些小分子VFAs的濃度在高OLR條件下仍然維持在較低水平,產(chǎn)甲烷性能良好;而對照組則在OLR7.8kgCOD/(m3·d)時就出現(xiàn)了VFAs的大量積累,短時間內(nèi)產(chǎn)甲烷能力幾乎完全喪失.微生物群落結(jié)構(gòu)的分析結(jié)果可以看出,具有胞外電子傳遞能力的細菌在GAC反應(yīng)器內(nèi)大量富集,這些微生物可以將降解有機物產(chǎn)生的電子傳遞到細胞外,并通過GAC的導(dǎo)電性傳遞給能夠直接接受電子的產(chǎn)甲烷菌,古菌的分析結(jié)構(gòu)進一步驗證了這一猜想,具有直接接受電子能力的產(chǎn)甲烷菌和均在GAC表面明顯富集,表明GAC在反應(yīng)體系中強化了這兩類微生物的富集生長,并且起到了電子傳遞載體的作用,加速了這兩類微生物之間的電子傳遞效率.

截止到目前,通過微生物純培養(yǎng)證實的DIET過程均有產(chǎn)電菌的參與,的體外導(dǎo)電菌毛(e-pili)和細胞色素在DIET過程中起到了決定作用[23-24].然而,Liu等[7]發(fā)現(xiàn)GAC的存在可以解除發(fā)生DIET所依賴的e-pili和細胞色素.因此,GAC的介入為除了以外的其他具有胞外電子傳遞能力的微生物與產(chǎn)甲烷菌之間發(fā)生電子傳遞提供了可能性.對于該猜想,大多數(shù)研究還停留在混合菌群的體系,更直接的證據(jù)需要在今后通過純菌培養(yǎng)實驗進一步驗證.

3 結(jié)論

3.1 投加GAC可以顯著提升厭氧消化處理餐廚垃圾的處理效率,穩(wěn)定運行的最高OLR可達10.4kgCOD/(m3·d),相比于對照組提升了33.3%以上.

3.2 在高有機負荷條件下,投加GAC反應(yīng)器內(nèi)VFAs的積累程度明顯降低,甲烷產(chǎn)率提高,表明投加GAC可以強化VFAs的降解和甲烷的轉(zhuǎn)化.

3.3 實驗組GAC表面和懸浮污泥富集生長了大量具有胞外電子傳遞能力的細菌(豐度分別為34.1%、26.7%)和具有直接接受電子能力的產(chǎn)甲烷菌(豐度分別為18.0%、22.8%),對照組懸浮污泥中上述兩類微生物豐度分別為7.5%和小于1.0%,推測GAC作為電子傳遞載體將這部分細菌降解有機污染物產(chǎn)生的電子傳遞給產(chǎn)甲烷菌,因此強化了餐廚垃圾的處理效率.

[1] Environmental Protection Agency US. Turning Food Waste into Energy at the East Bay Municipal Utility District (EBMUD) [R]. Pacific Southwest, 2007, Region 9U.S.EPA.

[2] 郝曉地,周 鵬,曹達啓.餐廚垃圾處置方式及其碳排放分析 [J]. 環(huán)境工程學(xué)報, 2017,11(2):673-682.

[3] 汪群慧,馬鴻志,王旭明,等.廚余垃圾的資源化技術(shù) [J]. 現(xiàn)代化工, 2004,24(7):56-59.

[4] 王 琦,李亞紅,蔡偉民.提高廚余垃圾厭氧消化甲烷產(chǎn)量的研究進展 [J]. 環(huán)境科學(xué)與技術(shù), 2006,29(B08):130-132.

[5] Zhang C, Su H, Baeyens J, et al. Reviewing the anaerobic digestion of food waste for biogas production [J]. Renewable and Sustainable Energy Reviews, 2014,38:383-392.

[6] Rotaru A E, Shrestha P M, Liu F, et al. A new model for electron flow during anaerobic digestion: direct interspecies electron transfer tofor the reduction of carbon dioxide to methane [J]. Energy & Environmental Science, 2014,7(1):408-415.

[7] Liu F, Rotaru A E, Shrestha P M, et al. Promoting direct interspecies electron transfer with activated carbon [J]. Energy & Environmental Science, 2012,5(10):8982-8989.

[8] Zhao Z, Zhang Y, Holmes D E, et al. Potential enhancement of direct interspecies electron transfer for syntrophic metabolism of propionate and butyrate with biochar in up-flow anaerobic sludge blanket reactors [J]. Bioresource technology, 2016,209:148-156.

[9] Dang Y, Sun D, Woodard T L, et al. Stimulation of the anaerobic digestion of the dry organic fraction of municipal solid waste (OFMSW) with carbon-based conductive materials [J]. Bioresource Technology, 2017,238:30-38.

[10] Timur H, ?zturk I. Anaerobic sequencing batch reactor treatment of landfill leachate [J]. Water Research, 1999,33(15):3225-3230.

[11] Sasaki K, Morita M, Hirano S I, et al. Decreasing ammonia inhibition in thermophilic methanogenic bioreactors using carbon fiber textiles [J]. Applied Microbiology and Biotechnology, 2011,90(4):1555.

[12] Fernández J, Pérez M, Romero L I. Effect of substrate concentration on dry mesophilic anaerobic digestion of organic fraction of municipal solid waste (OFMSW) [J]. Bioresource Technology, 2008,99(14):6075–6080.

[13] Kim G, Hyun M, Chang I, et al. Dissimilatory Fe (III) reduction by an electrochemically active lactic acid bacterium phylogenetically related toisolated from submerged soil [J]. Journal of Applied Microbiology, 2005, 99(4):978-987.

[14] Hernandez-Eugenio G, Fardeau M L, Cayol J L, et al.gen.nov., sp.nov., a novel acetogenic, facultatively sulfur-reducing bacterium [J]. International Journal of Systematic and Evolutionary Microbiology, 2002,52(4):1217-1223.

[15] Lovley D R, Holmes D E, Nevin K P. Dissimilatory Fe (III) and Mn (IV) reduction [M]. Advances in Microbial Physiology, 2004, 49:219-286.

[16] Cord-Ruwisch R, Garcia J L. Isolation and characterization of an anaerobic benzoate-degrading spore-forming sulfate-reducing bacterium,sp. Nov [J]. FEMS Microbiology Letters, 1985,29(3):325-330.

[17] Lovley D R, Giovannoni S, White D, et al.gen. nov. sp. nov., a microorganism capable of coupling the complete oxidation of organic compounds to the reduction of iron and other metals [J]. Archives of Microbiology, 1993,159(4):336-344.

[18] Park D H, Zeikus J G. Improved fuel cell and electrode designs for producing electricity from microbial degradation [J]. Biotechnology and Bioengineering, 2003,81(3):348-355.

[19] Kaksonen A H, Spring S, Schumann P, et al.gen. nov.,sp. nov. a thermophilic sulfate-reducer isolated from a geothermal underground mine in Japan [J]. International Journal of Systematic and Evolutionary Microbiology, 2007,57(1):98-102.

[20] Rotaru A E, Shrestha P M, Liu F,et al. Direct interspecies electron transfer betweenand[J]. Applied and Environmental Microbiology, 2014,80(15):4599-4605.

[21] Barredo M, Evison L. Effect of propionate toxicity on methanogen-enriched sludge,, andat different pH values [J]. Applied and Environmental Microbiology, 1991,57(6):1764-1769.

[22] Cruz Viggi C, Rossetti S, Fazi S, et al.Magnetite particles triggering a faster and more robust syntrophic pathway of methanogenic propionate degradation [J]. Environmental Science & Technology, 2014,48(13):7536-7543.

[23] Lovley D R. Happy together: microbial communities that hook up to swap electrons [J].The ISME Journal, 2017,11(2):327-336.

[24] Lovley D R, Ueki T, Zhang T, et al.: the microbe electric’s physiology, ecology, and practical applications [M]. Advances in Microbial Physiology, 2011,59(1):1-100.

Enhancing anaerobic methanogenic treatment of kitchen wastes by incorporating granular activated carbon.

FENG Xian-lu, LIU Xin-ying, SUN De-zhi, LIU Jie-qi, HUO Da, MA Jia-yuan, ZHANG Zhen, DANG Yan*

(Beijing Key Laboratory for Source Control Technology of Water Pollution, Engineering Research Center for Water Pollution Source Control and Eco-remediation, College of Environmental Science and Engineering, Beijing Forestry University, Beijing 100083, China)., 2018,38(3)：1018~1023

This paper investigated the hypothesis that enhancing anaerobic digestion of kitchen wastes by accelerating DIET through incorporating granular activated carbon (GAC). Besides, changes of microbial communities by the incorporation of GAC were also studied. The experimental results showed that, GAC reactors could operate stably with high methane production rate under the organic loading rate (OLR) as high as 10.4kgCOD/(m3·d). In contrast, for the control reactors without GAC incorporation, the methane production rate and pH declined sharply when the OLR increased to only 7.8kgCOD/(m3·d). In addition, the volatile fatty acids severelyaccumulated and resulted in the treatment efficiency of control reactors deteriorated. Microbial community analysis showed that bacteria capable of extracellular electron transfer (34% of bacterial community) and methanogens known to participate in DIET (88% of archaeal community) were significantly enriched on the GAC surface. It demonstrated that the addition of GAC could enrich these two groups of microbes and enhance the anaerobic digestion of kitchen waste through DIET.

direct interspecies electron transfer；methanogenesis；kitchen waste；granular activated carbon

X705

A

1000-6923(2018)03-1018-06

馮顯露(1996-),女,貴州黔東南人,北京林業(yè)大學(xué)碩士研究生,主要從事廢水的生物處理研究.發(fā)表論文1篇.

2017-08-19

北京林業(yè)大學(xué)青年教師科學(xué)研究中長期項目(2015ZCQ- HJ-01);國家自然科學(xué)基金資助項目(51708031);大學(xué)生創(chuàng)新訓(xùn)練項目(240-131705003)

* 責(zé)任作者, 講師, yandang@bjfu.edu.cn