孫 慶,武旭斌,李學(xué)德,花日茂
合肥草莓基地土壤重金屬污染風(fēng)險(xiǎn)評(píng)價(jià)
孫 慶,武旭斌,李學(xué)德,花日茂*
(農(nóng)產(chǎn)品質(zhì)量安全安徽省重點(diǎn)實(shí)驗(yàn)室,安徽農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,合肥 230036)
采集合肥市某草莓基地區(qū)域100個(gè)樣點(diǎn)的土壤,測(cè)定表層土壤(深度分別為0~10 cm和10~20 cm)的6種重金屬(Cd、Hg、As、Pb、Cr和Cu),利用地累積指數(shù)方法對(duì)草莓種植土壤重金屬的污染程度和生態(tài)風(fēng)險(xiǎn)進(jìn)行評(píng)價(jià),并結(jié)合多元統(tǒng)計(jì)分析探討重金屬來(lái)源。結(jié)果表明,Pb、Cd、Hg、Cu和Cr雖對(duì)草莓基地部分區(qū)域的土壤有一定程度的累積和污染,但污染風(fēng)險(xiǎn)低;主成分分析表明,Pb、Cd、As、Cu和Hg源于人類(lèi)農(nóng)業(yè)活動(dòng),Cr源于土壤母質(zhì)及風(fēng)化累積。顯示草莓基地整體區(qū)域的土壤重金屬生態(tài)風(fēng)險(xiǎn)低,種植基地土壤重金屬未對(duì)草莓質(zhì)量安全產(chǎn)生影響。
重金屬;土壤;草莓;生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)
草莓為多年生草本植物,因其富含多種營(yíng)養(yǎng)成分和香甜的口感而被譽(yù)為“水果皇后”[1-3]。自其進(jìn)入我國(guó)以來(lái),深受消費(fèi)者的喜愛(ài),因此,草莓種植在我國(guó)快速發(fā)展,目前我國(guó)的草莓種植面積和產(chǎn)量已經(jīng)位居全球第一[4-5]。土壤是草莓種植生產(chǎn)中不可或缺的組成部分,隨著社會(huì)的發(fā)展,土壤各種污染也隨之而來(lái),尤其是土壤重金屬的污染。一些養(yǎng)殖業(yè)生產(chǎn)過(guò)程中的動(dòng)物糞便含有的重金屬被作為有機(jī)肥使用,會(huì)給農(nóng)田土壤帶來(lái)潛在風(fēng)險(xiǎn)。目前,對(duì)土壤的重金屬污染已有眾多學(xué)者進(jìn)行研究[6-7]。本研究以草莓基地的土壤為研究對(duì)象,研究了草莓種植基地土壤中的Cd、Hg、As、Pb、Cr和Cu 6種重金屬的含量分布特征,利用地累積指數(shù)對(duì)土壤重金屬污染程度和生態(tài)風(fēng)險(xiǎn)進(jìn)行評(píng)估,為草莓生產(chǎn)和食用安全,促進(jìn)草莓產(chǎn)業(yè)的健康發(fā)展提供科學(xué)依據(jù)。
某草莓基地位于合肥市區(qū)域內(nèi),該區(qū)域土地平坦,屬于亞熱帶季風(fēng)氣候,氣候條件比較優(yōu)越,雨水充足,光照充足[8]。研究區(qū)域草莓基地的土壤類(lèi)型相對(duì)單一,且具有氣候溫暖濕潤(rùn)的良好條件,因此十分適合種植草莓。目前,草莓產(chǎn)業(yè)已成為當(dāng)?shù)氐奶厣鲗?dǎo)產(chǎn)業(yè),并且已經(jīng)成為當(dāng)?shù)剞r(nóng)業(yè)的主要經(jīng)濟(jì)支柱之一。
1.2.1 樣品采集 調(diào)查研究土壤樣品采集在草莓種植前進(jìn)行,草莓樣品待草莓種植成熟后進(jìn)行采集。土壤采集在5 000畝研究區(qū)域草莓基地內(nèi)設(shè)置100個(gè)取樣點(diǎn),每個(gè)采樣點(diǎn)內(nèi)分別以下0~10 cm和10~20 cm兩種表層土壤進(jìn)行系統(tǒng)采集。采集好的土壤樣品挑出石塊、枯葉等雜物后混勻,再將其風(fēng)干、破碎、過(guò)篩。
1.2.2 樣品測(cè)試 研究測(cè)定的是常見(jiàn)6種重金屬元素:Pb、Cd、Hg、As、Cu和Cr。
土壤重金屬各元素分析測(cè)試方法為:Pb和Cd的測(cè)定采用石墨爐原子吸收分光光度法(GB/T 17141—1997)[9],Hg和As的測(cè)定采用微波消解/原子熒光法(HJ 680-2013)[10],Cu和Cr的測(cè)定采用火焰原子吸收光度法(HJ 491—2019)[11]。
草莓產(chǎn)品重金屬含量按照食品安全國(guó)家標(biāo)準(zhǔn)食品中多元素的測(cè)定。其中Pb、Cd、As、Cu和Cr為食品安全國(guó)家標(biāo)準(zhǔn)(GB 5009.268-2016 第一法)[12]測(cè)定,Hg為食品安全標(biāo)準(zhǔn)(GB 5009.17-2014 第一篇第一法)[13]測(cè)定。
C為樣品土壤中重金屬元素a的濃度;B為當(dāng)?shù)氐耐寥乐亟饘賏的背景濃度,本次研究采用安徽A層土壤元素背景值;1.5為修正指數(shù)。geo≤0(Ⅰ級(jí))未受污染;0<geo≤1(Ⅱ級(jí))輕度污染;1<geo≤2(Ⅲ級(jí))中度污染;2<geo≤3(Ⅳ級(jí))中度污染到重度污染;3<geo≤4(Ⅴ級(jí))重度污染;4<geo≤5(Ⅵ級(jí))重度污染到極度污染;geo>5(Ⅶ級(jí))極度污染。
本次調(diào)查研究選用的元素參比值為安徽省土壤元素算術(shù)平均值[18]。從表1可以看出,地表0~10 cm深度的淺表層土壤中Cd、Hg、As、Pb、Cr和Cu的平均濃度分別為0.118、0.038、8.78、28.02、75.11和28.89 mg·kg-1。除As的平均濃度低于背景值9 mg·kg-1, Cd、Hg、Pb、Cr和Cu的平均濃度均高出了對(duì)應(yīng)的環(huán)境背景值。10~20 cm表層土壤中Cd、Hg、As、Pb、Cr和Cu的平均濃度分別為0.111、0.039、8.86、27.96、75.94和28.48 mg·kg-1,與0~0 cm表層土壤一樣,除了As的平均濃度低于背景值9 mg·kg-1,其余Cd、Hg、Pb、Cr和Cu的平均濃度均高出了對(duì)應(yīng)的安徽省的元素環(huán)境背景值。各個(gè)取樣點(diǎn)0~20 cm土壤中的6種重金屬Cd、Hg、As、Pb、Cr和Cu的平均濃度分別為0.115、0.039、8.82、27.99、75.52和28.69 mg·kg-1。除As的平均濃度低于背景值9 mg·kg-1,其余Cd、Hg、Pb、Cr和Cu的平均濃度均高出了對(duì)應(yīng)的環(huán)境背景值。對(duì)照《土壤環(huán)境質(zhì)量農(nóng)用地土壤污染風(fēng)險(xiǎn)管控標(biāo)準(zhǔn)(試行)(GB15618-2018)》[19],采用較嚴(yán)格的風(fēng)險(xiǎn)篩選值進(jìn)行評(píng)價(jià),所測(cè)6種重金屬100個(gè)取樣點(diǎn)位中的200樣品均未有樣品超出風(fēng)險(xiǎn)篩選值,未超過(guò)土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)。由此可見(jiàn)研究區(qū)域土壤重金屬呈現(xiàn)較低程度的累積。變異系數(shù)可以在一定程度上反映干擾活動(dòng)的影響,根據(jù)變異系數(shù)的分級(jí)標(biāo)準(zhǔn):變異系數(shù)大于30%為強(qiáng)變異,10%至30%為中等變異,小于10%為弱變異。草莓基地土壤中Cd和Hg變異系數(shù)大于30%,為強(qiáng)變異。其余4種重金屬變異系數(shù)在10%和30%之間,為中度變異。由此表明草莓基地土壤重金屬有可能受外界因素的影響。
從表2看出,根據(jù)污染等級(jí)的劃分,草莓基地0~10 cm表層土壤所有采樣點(diǎn)As的geo值為Ⅰ級(jí),這表明土壤未受As的影響。Cd、Hg、Pb、Cr和Cu的geo值范圍均是從Ⅰ級(jí)到Ⅱ級(jí),分別為78%、81%、95%、97%和96%,Ⅱ級(jí)占比分別為22%、19%、5%、3%和4%。草莓基地10~20 cm表層土壤所有采樣點(diǎn)As的geo值為Ⅰ級(jí),表明目前土壤未受As的影響。Cd、Hg、Pb、Cr和Cu的geo值范圍均是從Ⅰ級(jí)到Ⅱ級(jí),分別為86%、85%、97%、97%和97%,Ⅱ級(jí)占比分別為14%、15%、3%、3%和3%。將各個(gè)取樣點(diǎn)0~20 cm表層土壤數(shù)據(jù)進(jìn)行分析評(píng)價(jià),Cd、Hg、As、Pb、Cr和Cu的geo值范圍均是從Ⅰ級(jí)到Ⅱ級(jí),分別為82%、85%、100%、96%、98%和96%,Ⅱ級(jí)的占比分別為18%、15%、0%、4%、2%和4%。3組數(shù)據(jù)表明,重金屬Pb、Cd、Hg、Cu和Cr雖對(duì)草莓基地的土壤有一定的污染,但污染的程度不高,草莓基地大多數(shù)地區(qū)未受到這幾種重金屬的污染。
注:樣本數(shù)=100,GB15618-2018為《土壤環(huán)境質(zhì)量農(nóng)用地土壤污染風(fēng)險(xiǎn)管控標(biāo)準(zhǔn)(執(zhí)行)》[19]。
表 2 土壤重金屬地質(zhì)累積指數(shù)污染分級(jí)統(tǒng)計(jì)
表3 Perason相關(guān)分析結(jié)果
2.3.1 相關(guān)性分析 運(yùn)用SPSS軟件,通過(guò)皮爾遜(Pearson)相關(guān)性分析確定草莓基地6種重金屬濃度之間的相互關(guān)系,初步探究草莓基地的土壤重金屬的可能相似來(lái)源[20-21]。重金屬之間的相關(guān)性越顯著,表明其來(lái)源途徑越接近[22-23]。由表3可以看出,在<0.01水平上,Pb與Cd、As、Cu呈極顯著相關(guān);Cd與Pb、Cu之間呈極顯著相關(guān);As與Pb、Cu呈極顯著相關(guān);Cu與Pb、Cd、As之間存在較強(qiáng)的相關(guān)性。Pb、Cd、As和Cu 4種重金屬元素之間有著相似的來(lái)源途徑,Hg和Cr與其他重金屬之間相關(guān)性不顯著,表明這兩種重金屬元素來(lái)源跟其他重金屬元素有所差別。
表4 土壤重金屬元素主成分分析
表5 土壤重金屬初始因子載荷矩陣
2.3.2 主成分分析 基于相關(guān)性分析結(jié)果,對(duì)6種重金屬進(jìn)行了主成分分析(表4和表5),提取前3個(gè)特征值大于1的主成分,其累計(jì)貢獻(xiàn)率為73.36%,因此結(jié)合草莓基地土壤的利用情況對(duì)這3個(gè)主成分進(jìn)行分析,可以代表本次研究結(jié)果。由表4和表5可以看出,第一主成分的貢獻(xiàn)率為34.46%,在Pb、Cd、As、Cu含量上有較高的正載荷,由相關(guān)性可知這4種重金屬來(lái)源途徑相似,本研究結(jié)果與夏文建等[24]和王飛騰等[25]提出的化肥及有機(jī)肥顯著提高了土壤中Cu、Cd和As等重金屬的含量結(jié)果一致。且在研究區(qū)域草莓基地的各重金屬污染程度較低或無(wú)污染,說(shuō)明第一主成分主要表征為人類(lèi)施用化肥有機(jī)肥的農(nóng)業(yè)活動(dòng)。第二主成分的貢獻(xiàn)率21.51%,在Hg含量上有較高的正載荷,含Hg的農(nóng)藥能給農(nóng)用地土壤帶來(lái)污染[26-27],結(jié)合變異系數(shù),說(shuō)明草莓基地Hg的第二主成分表征為人類(lèi)使用農(nóng)藥影響。第三主成分貢獻(xiàn)率為17.39%,在Cr含量上有較高的正載荷,研究表明Cr主要受土壤母質(zhì)及風(fēng)化累積的作用[28],所以第三主成分主要表征為自然來(lái)源。
由表6可知,草莓樣品檢測(cè)出了6種重金屬中的3種。所有樣品都檢測(cè)出Cu元素,有2個(gè)樣品檢測(cè)出As,還有1個(gè)樣品檢測(cè)出了Cd。Cu是人體的必需元素,表明草莓是含豐富Cu元素的水果,對(duì)照食品中Cu限量衛(wèi)生標(biāo)準(zhǔn)(GB 15199—94)[29],其含量遠(yuǎn)低于水果中Cu含量最高容許限量10 mg·kg-1。根據(jù)食品安全國(guó)家標(biāo)準(zhǔn)食物中污染物限量標(biāo)準(zhǔn)(GB 2762—2017)[30],水果中5種重金屬最高限量(Cd-0.05 mg·kg-1,Hg-0.01 mg·kg-1,As-0.5 mg·kg-1,Pb-0.2 mg·kg-1,Cr-0.5 mg·kg-1)。所有草莓樣品中這5種重金屬的含量低于最高限量。結(jié)合土壤的污染狀態(tài),草莓基地土壤重金屬未對(duì)草莓質(zhì)量安全產(chǎn)生影響。
表6 草莓重金屬含量特征
草莓基地土壤樣品中的所測(cè)重金屬元素Pb、Cd、Hg、Cu和Cr的均值均超出了安徽省的土壤元素背景值,As的均值未超出背景值。在草莓基地100樣品點(diǎn)位中,200份土壤樣品的6種重金屬含量均未超出對(duì)應(yīng)的重金屬風(fēng)險(xiǎn)篩選值,未超過(guò)土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)。表明草莓基地的土壤風(fēng)險(xiǎn)低,一般情況下可忽略其影響。
根據(jù)污染評(píng)價(jià)方法的結(jié)果來(lái)看,6種重金屬Pb、Cd、Hg、Cu、As和Cr對(duì)草莓基地的土壤有一定的污染,但污染的程度不高。草莓基地整體區(qū)域的土壤重金屬的生態(tài)風(fēng)險(xiǎn)低,屬于需要保護(hù)類(lèi)型的土壤。主成分分析表明,Pb、Cd、As、Cu源于人類(lèi)施用化肥、有機(jī)肥的農(nóng)業(yè)活動(dòng),Hg源于農(nóng)業(yè)生產(chǎn)中施用的農(nóng)藥,Cr源于土壤母質(zhì)及風(fēng)化累積。
草莓基地生產(chǎn)的草莓,其含量遠(yuǎn)低于標(biāo)準(zhǔn)食品安全國(guó)家標(biāo)準(zhǔn)食物中污染物限量(GB 2762—2017)[30],表明草莓基地土壤重金屬未對(duì)草莓產(chǎn)品質(zhì)量安全產(chǎn)生影響。
[1] 茍軍. 大棚草莓種植技術(shù)[J]. 農(nóng)民致富之友, 2015(12): 162-163.
[2] KOWALSKA J, KOWALSKA H, MARZEC A, et al. Dried strawberries as a high nutritional value fruit snack[J]. Food Sci Biotechnol, 2018, 27(3): 799-807.
[3] 楊會(huì)苗, 毛清云, 劉蘭英. 北京和保定地區(qū)草莓產(chǎn)業(yè)基本情況調(diào)研報(bào)告[J]. 北方園藝, 2019(6): 170-174.
[4] 張欣馨, 王菲, 李浪, 等. 中國(guó)草莓生產(chǎn)中面臨的主要問(wèn)題及發(fā)展對(duì)策[J]. 中國(guó)林副特產(chǎn), 2016(2): 92-96.
[5] 王鳴謙, 薛莉, 趙珺, 等. 世界草莓生產(chǎn)及貿(mào)易現(xiàn)狀[J]. 中國(guó)果樹(shù), 2021(2): 104-108.
[6] 毛慧, 姚軍, 吳晶. 草莓園土壤中重金屬含量的研究與評(píng)價(jià)[J]. 現(xiàn)代農(nóng)業(yè)科技, 2016(19): 181-182, 184.
[7] 宋偉, 陳百明, 劉琳. 中國(guó)耕地土壤重金屬污染概況[J]. 水土保持研究, 2013, 20(2): 293-298.
[8] 秦倩凝. 地理環(huán)境與城市個(gè)性的交融: 以合肥市為例[J]. 科教文匯(中旬刊), 2013(2): 201-203.
[9] 國(guó)家環(huán)境保護(hù)局規(guī)劃標(biāo)準(zhǔn)處. 土壤質(zhì)量鉛、鎬的測(cè)定石墨爐原子吸收分光光度法:GB/T 17141-1997[S]. 北京:中國(guó)標(biāo)準(zhǔn)出版社,1997.
[10] 環(huán)境保護(hù)部科技標(biāo)準(zhǔn)司. 土壤和沉積物汞、砷、硒、鉍、銻的測(cè)定微波消解/原子熒光法:HJ 680-2013[S]. 北京:中國(guó)環(huán)境科學(xué)出版社,2013.
[11] 生態(tài)環(huán)境部生態(tài)環(huán)境監(jiān)測(cè)司.土壤和沉積物銅、鋅、鉛、鎳、鉻的測(cè)定火焰原子吸收分光光度法HJ 491- 2019[S]. 北京:中國(guó)環(huán)境出版社,2019.
[12] 國(guó)家衛(wèi)生和計(jì)劃生育委員會(huì), 國(guó)家食品藥品監(jiān)督管理總局. 食品安全國(guó)家標(biāo)準(zhǔn)食品中多元素的測(cè)定: GB 5009.268-2016[S].北京:中國(guó)標(biāo)準(zhǔn)出版社, 2017.
[13] 國(guó)家衛(wèi)生和計(jì)劃生育委員會(huì).食品安全國(guó)家標(biāo)準(zhǔn)食品中總汞及有機(jī)汞的測(cè)定:GB 5009.17-2014[S].北京:中國(guó)標(biāo)準(zhǔn)出版社, 2014.
[14] 郭志娟, 周亞龍, 王喬林, 等. 雄安新區(qū)土壤重金屬污染特征及健康風(fēng)險(xiǎn)[J]. 中國(guó)環(huán)境科學(xué), 2021, 41(1): 431-441.
[15] 韋壯綿, 陳華清, 張煜, 等. 湘南柿竹園東河流域農(nóng)田土壤重金屬污染特征及風(fēng)險(xiǎn)評(píng)價(jià)[J]. 環(huán)境化學(xué), 2020, 39(10): 2753-2764.
[16] ZENG Y, LIU J, LI Y Q, et al. Accumulation and health implications of metals in topsoil of an urban riparian zone adjacent to different functional areas in a subtropical city[J]. J Environ Qual, 2021, 50(3): 744-755.
[17] TUME P, ACEVEDO V, ROCA N, et al. Potentially toxic elements concentrations in schoolyard soils in the city of Coronel, Chile[J]. Environ Geochem Health, 2022, 44(5): 1521-1535.
[18] 成杭新, 李括, 李敏, 等. 中國(guó)城市土壤化學(xué)元素的背景值與基準(zhǔn)值[J]. 地學(xué)前緣, 2014, 21(3): 265-306.
[19] 生態(tài)環(huán)境部,國(guó)家市場(chǎng)監(jiān)督管理總局.土壤環(huán)境質(zhì)量農(nóng)用地土壤污染風(fēng)險(xiǎn)管控標(biāo)準(zhǔn)(試行):GB15618-2018[S]. 北京:中國(guó)環(huán)境出版社,2018.
[20] 馮乾偉, 王兵, 馬先杰, 等. 黔西北典型鉛鋅礦區(qū)土壤重金屬污染特征及其來(lái)源分析[J]. 礦物巖石地球化學(xué)通報(bào), 2020, 39(4): 863-870.
[21] 滑小贊, 程濱, 趙瑞芬, 等. 太原市農(nóng)田土壤養(yǎng)分和重金屬累積特征及相關(guān)性分析[J]. 山西農(nóng)業(yè)科學(xué), 2019, 47(6): 1027-1033.
[22] 曾妍妍, 周金龍, 王松濤, 等. 新疆若羌縣農(nóng)田土壤重金屬分布特征及污染評(píng)價(jià)[J]. 干旱區(qū)資源與環(huán)境, 2017, 31(9): 87-91.
[23] 梁雅雅, 易筱筠, 黨志, 等. 某鉛鋅尾礦庫(kù)周邊農(nóng)田土壤重金屬污染狀況及風(fēng)險(xiǎn)評(píng)價(jià)[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2019, 38(1): 103-110.
[24] 夏文建, 張麗芳, 劉增兵, 等. 長(zhǎng)期施用化肥和有機(jī)肥對(duì)稻田土壤重金屬及其有效性的影響[J]. 環(huán)境科學(xué), 2021, 42(5): 2469-2479.
[25] 王騰飛, 譚長(zhǎng)銀, 曹雪瑩, 等. 長(zhǎng)期施肥對(duì)土壤重金屬積累和有效性的影響[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2017, 36(2): 257-263.
[26] 王信凱, 張艷霞, 黃標(biāo), 等. 長(zhǎng)江三角洲典型城市農(nóng)田土壤重金屬累積特征與來(lái)源[J]. 土壤學(xué)報(bào), 2021, 58(1): 82-91.
[27] 孫淑蘭. 汞的來(lái)源、特性、用途及對(duì)環(huán)境的污染和對(duì)人類(lèi)健康的危害[J]. 上海計(jì)量測(cè)試, 2006, 33(5): 6-9.
[28] 馬憲梅, 黃曉飛. 土壤鉻污染現(xiàn)狀及修復(fù)方法研究[J]. 環(huán)境與發(fā)展, 2020, 32(5): 61.
[29] 衛(wèi)生部衛(wèi)生監(jiān)督司.食品中銅限量衛(wèi)生標(biāo)準(zhǔn): GB15199-94[S]. 北京:中國(guó)標(biāo)準(zhǔn)出版社,1994.
[30] 國(guó)家衛(wèi)生和計(jì)劃生育委員會(huì),國(guó)家食品藥品監(jiān)督管理總局.食品安全國(guó)家標(biāo)準(zhǔn):食品中污染物限量:GB 2762- 2017 [S]. 北京:中國(guó)標(biāo)準(zhǔn)出版社,2017.
Risk assessment of heavy metal pollution in soil of Hefei strawberry base
SUN Qing, WU Xubing, LI Xuede, HUA Rimao
(Key Laboratory of Agri-food Safety of Anhui Province, School of Resource and Environment, Anhui Agricultural University, Hefei 230036)
The contents of six heavy metals (Cd, Hg, As, Pb, Cr and Cu) were measured in the surface soil (depth in 0-10 cm and 10-20 cm, respectively) from 100 sampling sites in a strawberry base area in Hefei, and the degree of contamination and ecological risk of heavy metals in strawberry growing soil were evaluated using the ground accumulation index method and combined with multivariate statistical analysis to explore the sources of heavy metals. The results showed that Pb, Cd, Hg, Cu and Cr had accumulated and polluted soil in some areas of strawberry base to a certain extent, but the pollution risk was low. The principal component analysis showed that Pb, Cd, As Cu and Hg were derived from human agricultural activities. Cr originates from soil parent material and weathering accumulation. It showed that the ecological risk of soil heavy metals in the whole area of strawberry base was low, and the heavy metals in the soil have not affected the quality and safety of strawberries.
heavy metals; soil; strawberry; ecological risk assessment
X53; X820.4
A
1672-352X (2022)05-0794-05
10.13610/j.cnki.1672-352x.20221111.015
2022-11-15 9:06:39
[URL] https://kns.cnki.net/kcms/detail/34.1162.s.20221111.1117.030.html
2021-04-19
安徽省農(nóng)業(yè)競(jìng)爭(zhēng)力項(xiàng)目(20803005)資助。
孫 慶,碩士研究生。E-mail:1577588976@qq.com
花日茂,博士,教授。E-mail:rimaohua@ahau.edu.cn