孟令儀,徐夢潔,李小曼,,劉 勤
(1.南京農(nóng)業(yè)大學公共管理學院,江蘇南京210095;2.中國科學院南京土壤研究所,江蘇南京 210008)
基于多種方法的蘇南鎮(zhèn)域土壤重金屬污染評價
孟令儀1,徐夢潔1,李小曼1,2,劉 勤2
(1.南京農(nóng)業(yè)大學公共管理學院,江蘇南京210095;2.中國科學院南京土壤研究所,江蘇南京 210008)
以蘇南小鎮(zhèn)為依托,采用地累積指數(shù)法、污染負荷指數(shù)法和內(nèi)梅羅指數(shù)結(jié)合土壤環(huán)境質(zhì)量標準的改進法對研究區(qū)域土壤重金屬污染進行評價,比較各種方法的適用性。結(jié)果表明,內(nèi)梅羅指數(shù)結(jié)合土壤環(huán)境質(zhì)量標準的改進法綜合了人為因素和自然因素的作用情況,分析結(jié)果更為客觀,且污染分布的區(qū)分度較高,兼有地累積指數(shù)法和污染負荷指數(shù)法的優(yōu)點,更適用于小尺度地區(qū)的土壤重金屬污染評價。
重金屬污染評價;改進方法;小尺度;蘇南
近年來,重金屬污染已經(jīng)引起研究人員的廣泛關(guān)注,其主要來源包括工業(yè)污染、金屬礦山開采、污水灌溉、固體廢棄物處置等方面[1]。我國在 “十二五”規(guī)劃通過了《重金屬污染綜合防治規(guī)劃》,說明對重金屬污染問題的重視已經(jīng)上升到國家層面。污染物在土壤中的聚集往往因人為活動的增強呈現(xiàn)出高度的空間變異性[2]。此外,土壤重金屬元素含量也會隨著時間的推移而發(fā)生改變[3],使得土壤重金屬污染問題更為復(fù)雜。由于土壤重金屬不易移動降解,進入生物體后會積累下來,并對生物體產(chǎn)生不利影響,甚至可以通過食物鏈影響人類健康[4]。因此,了解土壤中重金屬的含量與空間分布,并進行污染評價,對于土壤重金屬污染防治乃至維護人體健康都具有十分重要的意義。研究表明,人為來源是導(dǎo)致土壤重金屬污染問題的首要因素[5]。
在城鎮(zhèn)化發(fā)展飛速的蘇南地區(qū),因鄉(xiāng)鎮(zhèn)企業(yè)大力發(fā)展和監(jiān)管力度缺乏而導(dǎo)致的村鎮(zhèn)地區(qū)土壤重金屬污染問題十分嚴重[6]。鄉(xiāng)鎮(zhèn)是連接農(nóng)村和城市的橋梁與紐帶,在接受大中型城市輻射的同時,帶動著周圍廣大農(nóng)村經(jīng)濟和社會的發(fā)展,并促進農(nóng)民生活水平的提高[7]。隨著蘇南鄉(xiāng)鎮(zhèn)企業(yè)的快速發(fā)展,環(huán)境問題也隨之而來。由于鄉(xiāng)鎮(zhèn)監(jiān)管力度薄弱,鄉(xiāng)鎮(zhèn)企業(yè)的廠房設(shè)施一般較為簡陋,重金屬物質(zhì)排放不合理,土壤受到人為因素的嚴重改變[8],使得土壤重金屬污染成為蘇南鄉(xiāng)鎮(zhèn)環(huán)境的嚴峻問題[9]。由于重金屬在土壤中遷移緩慢,較小的尺度更能反映污染的真實性,所以開展以村鎮(zhèn)為單位的小尺度土壤重金屬污染研究顯得更為重要。
土壤重金屬污染研究中,評價方法也是十分重要的研究內(nèi)容。目前,國內(nèi)外應(yīng)用較多的土壤重金屬評價方法有內(nèi)梅羅指數(shù)法[10]、環(huán)境風險指數(shù)法[11]、灰色聚類法[12]等,其中內(nèi)梅羅指數(shù)法應(yīng)用最為廣泛,并對該方法進行了改進。在土壤重金屬污染評價時,所采用的基準值往往簡單劃一,未能因地制宜,反映研究區(qū)域的實際情況;研究人員也往往只使用一種評價方法,評價結(jié)果唯一,缺乏可比性,難以從中擇優(yōu)。同時,重金屬污染評價的范圍一般趨向大中尺度,小尺度的研究較為少見。筆者以蘇南小鎮(zhèn)為依托,采用地累積指數(shù)法、污染負荷指數(shù)法和內(nèi)梅羅結(jié)合土壤環(huán)境質(zhì)量標準的改進法對研究區(qū)域的重金屬污染進行評價,探討小尺度區(qū)域重金屬污染評價的適宜方法,旨在為政府部門以及環(huán)保部門的工作提供參考。
1.1 研究區(qū)概況 研究區(qū)域Y鎮(zhèn)位于江蘇省C市,地處我國經(jīng)濟最活躍的區(qū)域——上海經(jīng)濟圈中心,總面積約為12.67 km2。境內(nèi)四季分明,隸屬亞熱帶季風氣候,年平均氣溫為16 ℃;雨量充沛,水網(wǎng)稠密,年平均降水量為1 200 mm,土壤類型主要為水耕人為土[13]。Y鎮(zhèn)高速發(fā)展工業(yè)、苗木種植、水產(chǎn)養(yǎng)殖、農(nóng)業(yè)觀光旅游,其中工業(yè)產(chǎn)值占該鎮(zhèn)總產(chǎn)值的90%以上。工廠主要集中在東部,尤其是東南部地區(qū),主要有紡織廠、印染廠、機械制造廠、五金廠、溫室制造廠等。經(jīng)調(diào)查,Y鎮(zhèn)有大型銅管廠、金屬材料廠和電鍍廠,很多工廠的廢水、廢渣中都含有重金屬。
1.2 樣品采集與前期處理 采用格網(wǎng)布點法與半隨機布點法相結(jié)合的采樣方法。在鄉(xiāng)鎮(zhèn)工業(yè)區(qū)范圍以外的研究區(qū)內(nèi)以500 m格網(wǎng)布點;鄉(xiāng)鎮(zhèn)工業(yè)區(qū)內(nèi)由于廠房眾多,并且排列不規(guī)則,同時要避開道路等設(shè)施不便于采用格網(wǎng)的形式,因此采取加密布點的方式,加密時一部分考慮了主要污染源,一部分隨機布點,盡量覆蓋景觀草坪、居民地、工廠、水稻田、養(yǎng)殖水面等不同的土地利用類型。采樣時間為2014年10月21— 22日,用土壤樣品采樣器(型號為WN0306)采取用于重金屬分析的土壤樣品。采樣時用手持式GPS精確定位,采集表層0 ~ 20 cm土壤。其中,每個樣點分別向4個方向輻射約2 m(5 m×5 m方格內(nèi)),共采集5次。將5份樣品均勻混合制成該樣點的樣品,并標號,共采集110個樣。工廠附近200 m以內(nèi)的點占43.64%,農(nóng)作物及附近50 m以內(nèi)的點(包含居民區(qū)內(nèi)零碎菜地、農(nóng)業(yè)科技園等)占22.72%。
采集的土壤樣品帶回實驗室,置于陰涼處,室溫下使其自然風干。將風干的土壤樣品用木棒壓揉捻碎過20目的尼龍篩,并分成2份,1份樣品進行pH的測定,另1份樣品研磨后過100目的尼龍篩并充分混合均勻,裝瓶并貯存于干燥器內(nèi)備用。
1.3 樣品分析 土樣pH的測定使用便攜式pH計(SG8-FK)采用玻璃電極法進行測定,水土比為2.5∶1。測定土壤中Cr、Cu、Zn、Ni、Pb、Cd、As、Hg等重金屬含量,其中Cr、Cu、Zn、Ni、Pb、Cd前處理采用了HF-HCLO4-HNO3三酸高溫消解,Cr、Cu、Zn、Ni采用電感耦合等離子光譜法(ICP-AES)測定[14],Pb和Cd含量使用液相色譜電感耦合等離子體質(zhì)譜儀(HPLC-ICP-MS)測定,As和Hg含量使用GB/T 22105.1-3—2008王水溶液沸水浴消化原子熒光光度計法(AFS-230E、AFS-9700)測定。重金屬分析測試過程中均用GSS- 4、GSS-5標準參考土樣進行全程質(zhì)量控制。
1.4 主要評價方法
1.4.1 土壤重金屬污染評價標準。土壤背景值是土壤環(huán)境評價的基準值,同時也是土壤污染態(tài)勢預(yù)報、土壤環(huán)境容量計算、土壤環(huán)境質(zhì)量標準確立以及制定國民經(jīng)濟發(fā)展規(guī)劃的重要基礎(chǔ)數(shù)據(jù)[15]。為真實反映Y鎮(zhèn)的實際情況,筆者采用當?shù)赝寥辣尘爸底鳛樵u價基準(表1)。
該研究采用地累積指數(shù)法、污染負荷指數(shù)法和內(nèi)梅羅指數(shù)法結(jié)合土壤環(huán)境質(zhì)量標準的修正方法等多種方法評價Y鎮(zhèn)重金屬污染狀況。
表1 土壤質(zhì)量標準及背景值
1.4.2 地累積指數(shù)法。1969年,德國科學家Muller[18]提出地累積指數(shù)法,用于量化評價沉積物重金屬污染程度的指標[19],并根據(jù)計算得到的地累積指數(shù)I將評價結(jié)果的等級分為:無—中、中—強、強—極強。這種劃分使污染程度更為詳細,提醒人們及時控制污染,防止其惡化到更高級別。
1.4.3 污染負荷指數(shù)法。污染負荷指數(shù)法是Tomlinson等提出的基于土壤背景值和采樣點多種元素的一種求積的污染評價方法。該方法的評價公式為:
CFi=Ci/C0i
(1)
(2)
(3)
其中,Ci為元素i的實測值(mg/kg);C0i為元素i的評價標準(mg/kg);n為評價元素的個數(shù);m為評價點的個數(shù)(即采樣點的個數(shù));CFi為某單一金屬最高污染系數(shù);PLI為某點污染負荷指數(shù);PLIzone為評價區(qū)域污染負荷指數(shù)。若PLI<1,則污染等級0,無污染;若1 ≤PLI< 2,則污染等級 Ⅰ,中等污染;若2≤PLI<3,則污染等級Ⅱ,強污染;若PLI≥3,則污染等級Ⅲ,極強污染。
1.4.4 基于內(nèi)梅羅指數(shù)法結(jié)合土壤環(huán)境質(zhì)量標準的修正方法。內(nèi)梅羅指數(shù)法(綜合污染指數(shù)法)是一種兼顧極值或突出最大值的計權(quán)型多因子環(huán)境質(zhì)量指數(shù)法[20]。該方法包括單項污染指數(shù)法和綜合指數(shù)法,其中單項污染指數(shù)法能夠準確反映某種重金屬元素的污染程度,計算公式及分級標準參見相關(guān)文獻[1]。
直接采用內(nèi)梅羅指數(shù)法評價,僅在公式中依據(jù)當?shù)赝寥辣尘爸祷蛲寥拉h(huán)境質(zhì)量評價標準來計算,安全系數(shù)較小,且規(guī)范性和可比性欠佳。筆者在研究區(qū)域進行了評價試點,并最終確定采用內(nèi)梅羅指數(shù)結(jié)合土壤環(huán)境質(zhì)量標準修正評價法。此方法是對評價體系及公式都作出改進,而不是僅對公式進行改進,使得評價結(jié)果更加合理,不會過寬或過嚴,使得評價值比較充分合理地分布在相應(yīng)區(qū)間,可比性增強,更符合評價的實際需要。將國家土壤環(huán)境質(zhì)量標準(GB 15618—1995)的三級標準分別作為評價的各級污染的初始值,將土壤重金屬污染分為4級,并在劃分級別上也進行相應(yīng)調(diào)整。
采用上述方法分別對研究區(qū)土壤中的8種重金屬元素進行污染評價,對采樣點進行分級統(tǒng)計。為了更直觀地反映研究區(qū)的整體污染狀況,選用插值效果較好的徑向基函數(shù)法的高次曲面函數(shù)對樣點土壤環(huán)境質(zhì)量分類評價結(jié)果進行空間插值,并運用ArcGIS軟件進行等值線圖的繪制。
2.1 基于地累積指數(shù)的重金屬污染評價 地累積指數(shù)是通過扣除重金屬背景含量從而獲取由人為原因所引起的重金屬積累。該研究選取蘇南地區(qū)土壤背景值作為基準。由表2可知,污染強度最高的是Cd元素。除Ni、Cr和As以外,其他元素都達到中—強、強度污染,其中,Hg、Cd強度污染所占比例較高,而Hg、Pb、Cd在中—強及中度污染所占比例非常高,尤其是Pb中度污染的樣點達到68個,而Hg在中—強度污染樣點達到35個。
由表3可知,各元素的累積指數(shù)變化范圍最大的是Cd,說明Cd局部污染非常嚴重,區(qū)域污染程度參差不齊。I值均值從大到小排序依次為Pb、Cd、Hg、Cu、Ni、Zn、Cr、As。Pb元素極值污染程度不高,但整體污染嚴重,都大于1級,地累積指數(shù)法將其歸為中度污染。
表3 地累積指數(shù)均值評價結(jié)果
Table 3 Evaluation results of average values of geoaccumulation index
重金屬元素ElementsofheavymetalsI值等級范圍DegreeofIvalue評價結(jié)果(均值)Evaluationresults(mean)累積污染程度AccumulativepollutiondegreeCu0~30.45無—中Ni0~20.21無—中Zn0~40.02無—中Cr0~2-0.38無Cd0~51.49中Pb1~41.69中As0~1-0.70無Hg0~41.35中
2.2 基于污染負荷指數(shù)的重金屬污染評價
2.2.1 最高污染系數(shù)(CF)。以蘇南地區(qū)背景值作為評價標準,對采樣點8種重金屬的最高污染系數(shù)進行統(tǒng)計。從表4可以看出,8種元素的污染強度從強到弱依次為Hg、Pb、Cd、Cu、Zn、Ni、Cr、As,這與地累積指數(shù)的污染程度均值排序較為一致。這說明Hg和Pb的污染最為嚴重,對區(qū)域污染負荷貢獻最大。此外,各采樣點元素的污染系數(shù)極差較大,變化范圍大,元素之間相差幾倍甚至幾十倍,而標準偏差偏小,所以元素含量分布可能存在一定規(guī)律性,值得深究。
2.2.2 點污染評價。根據(jù)公式和分級標準可得到8種重金屬點污染負荷指數(shù)PLI、4種毒性強的重金屬點污染負荷指數(shù)PLIHg-As(Hg、Pb、Cd、As)及4種毒性相對較弱的重金屬點污染負荷指數(shù)PLICu-Cr(Cu、Ni、Zn、Cr)的統(tǒng)計結(jié)果。從8種重金屬污染負荷指數(shù)來看,沒有無污染級別的樣品,說明Y鎮(zhèn)土壤重金屬污染已經(jīng)十分嚴峻。從4種毒性強的重金屬點污染負荷指數(shù)來看,強污染以上樣品占到73.64%,極強污染占33.64%,說明毒性強的重金屬污染很嚴重。
表4 采樣點8種重金屬的最高污染系數(shù)分析
Table 4 The analysis of the highest pollution coefficient for 8 kinds of heavy metals in sampling points
重金屬元素Elementsofheavymetals平均值Mean最大值Maximum最小值Minimum極差Range標準差StandarddeviationCu2.196.620.915.710.95Ni1.875.780.085.690.77Zn1.8717.920.5817.341.94Cr1.265.070.055.020.70Cd5.1537.720.9536.774.25Pb5.1816.492.3414.152.10As0.931.780.521.260.15Hg5.4619.550.4119.154.18
從圖1可以看出,Hg、Pb、Cd、As綜合污染程度不僅從覆蓋范圍上,而且從污染級別上都要遠大于Cu、Ni、Zn、Cr。毒性強的元素污染總體分布在東北、西南區(qū)域,絕大部分污染程度在強污染以上,只有西北部地區(qū)污染在中污染水平,其余污染都在強污染以上。
2.2.3 區(qū)域污染評價。4種毒性強的重金屬區(qū)域污染負荷指數(shù)為2.91;4種毒性相對較弱的重金屬區(qū)域污染負荷指數(shù)為1.58;8種重金屬區(qū)域污染負荷指數(shù)為2.15。研究區(qū)重金屬污染整體為強污染,盡管毒性較弱的重金屬污染負荷指數(shù)只有中等污染,但是毒性強的重金屬污染指數(shù)接近極強污染,對生物健康危害極大。因此,綜合評價結(jié)果為強污染。
2.3 基于內(nèi)梅羅指數(shù)結(jié)合土壤環(huán)境質(zhì)量標準的重金屬污染評價 由表5和表6可知,Hg和Cd的污染程度較大。在中度污染部分,Ni、Hg、Cd的比例較修正前大幅提高,在輕度污染方面,除了Cd下降外,Hg和Ni基本持平,其余元素較修正前均有所上升,其中上升最為明顯的是Pb、Cu、Zn。改進后的多因子綜合污染評價,與修正前相比,P綜值從1.98上升為2.20,即從原來的輕度污染上升為中度污染。
圖1 點污染負荷指數(shù)評價結(jié)果Fig.1 Evaluation results of point pollution load index
表6 改進后單因子綜合污染指數(shù)評價
Table 6 The evaluation results of modified comprehensive pollution index of single factor
重金屬元素ElementsofheavymetalsPjPjmaxP綜2污染程度PollutiondegreeCu0.501.511.12輕度Ni1.093.362.50中度Zn0.545.233.72重度Cr0.281.110.81輕度Cd1.9914.5910.41重度Pb0.351.120.83清潔As0.330.630.50清潔Hg1.786.374.68重度
從評價結(jié)果來看,地累積指數(shù)法相對客觀,分級更加準確,沒有夸大局部值的效果,且相對污染負荷指數(shù)法更能反映人為因素和自然因素的綜合作用情況,分析結(jié)果更為客觀;污染負荷指數(shù)法比地累積指數(shù)法更能反映實際情況,且污染分布的區(qū)分度較高,更加形象清晰。
結(jié)合土壤環(huán)境質(zhì)量標準的內(nèi)梅羅指數(shù)法,改進后的單項污染指數(shù)污染程度普遍上升一個級別,并且增強了區(qū)分度,污染評價結(jié)果更為準確,對于Pb、Cd污染指數(shù)普遍較高的能夠合理地區(qū)分污染程度,而不是全部歸為無污染的級別,改進方法比改進前的規(guī)范性和可比性有所提高。改進后的指數(shù)計算方式更為科學、規(guī)范,不會出現(xiàn)過大或過小的值,不會使最后的污染指數(shù)偏差較大。改進后的單因子綜合指數(shù)排序發(fā)生了較大的變化,Zn位居第2,說明它保留了內(nèi)梅羅指數(shù)重視極值的特點;Pb和Cu排名有所提前,說明它也改進了內(nèi)梅羅指數(shù)不關(guān)注中間值及普遍污染程度的缺點。整體來看,污染程度較改進前有所上升,如Cu、Cr和Pb的污染級別升高,其他元素基本不變,這種評價方法更符合實際情況,優(yōu)勢十分明顯。內(nèi)梅羅結(jié)合土壤環(huán)境質(zhì)量標準的改進法兼有前2種方法的優(yōu)點。
該評價結(jié)果表明,毒性強的元素污染總體分布在東北和西南區(qū)域。人類活動密集的區(qū)域往往是污染來源最大的區(qū)域,同時也是污染風險較大的地區(qū),例如紡織廠、印染廠、機械制造廠、五金廠、溫室制造廠等。其中農(nóng)作物種植區(qū)的污染危害程度最大,往往需要保護,其余地區(qū)則需要進行防治。因此,在Y鎮(zhèn)這樣的小尺度區(qū)域,重金屬污染評價適宜采用內(nèi)梅羅指數(shù)的改進方法,因為該方法更符合實際,對于制訂針對性調(diào)控措施具有更大的參考價值。
[1] 范拴喜,甘卓亭,李梅娟,等.土壤重金屬污染評價方法進展[J].中國農(nóng)學通報,2010,26(17):310-315.
[2] LIN Y P.Multivariate geostatistical methods to identify and map spatial variations of soil heavy metals[J].Environmental geology,2002,42(1):1-10.
[3] 吳塹虹,戴塔根,方建武,等.長沙、株洲、湘潭三市土壤中重金屬元素的來源[J].地質(zhì)通報,2007,26(11):1453-1458.
[4] 陳懷滿.土壤-植物系統(tǒng)中的重金屬污染[M].北京:科學出版社,1996.
[5] HU Y N,CHENG H F.Application of stochastic models in identification and apportionment of heavy metal pollution sources in the surface soils of a large-scale region[J].Environmental science & technology,2013,47(8):3752-3760.
[6] 肖智,劉志偉,畢華.土壤重金屬污染研究述評[J].安徽農(nóng)業(yè)科學,2010 38(33):18812-18815.
[7] 方如康.環(huán)境學詞典[M].北京:科學出版社,2003:530.
[8] 盧瑛,龔子同,張甘霖,等.南京城市土壤重金屬含量及其影響因素[J].應(yīng)用生態(tài)學報,2004,15(1):123-126.
[9] 馬成玲,周健民,王火焰,等.農(nóng)田土壤重金屬污染評價方法研究:以長江三角洲典型縣級市常熟市為例[J].生態(tài)與農(nóng)村環(huán)境學報,2006,22(1):48-53.
[10] 潘勇軍,陳步峰,肖以華,等.廣州市城市森林土壤重金屬污染狀況及其評價[J].生態(tài)環(huán)境,2008,17(1):210-215.
[11] 沈體忠,肖杰.綠色食品產(chǎn)地土壤重金屬生態(tài)風險預(yù)警評估:以湖北省天門市為例[J].農(nóng)產(chǎn)品質(zhì)量與安全,2013(6):65-68.
[12] 何廳廳,趙艷玲,李建華,等.基于GCM_CB模型的土壤重金屬污染評價[J].環(huán)境工程學報,2012,6(10):3781-3786.
[13] RAN J,WANG D,WANG C,et al.Using portable X-ray fluorescence spectrometry and GIS to assess environmental risk and identify sources of trace metals in soils of per-urban areas in the Yangtze Delta region,China [J].Environmental science:Processes & impacts,2014,16(8):1870-1877.
[14] LEIVUORI M.Heavy metal contamination in surface sediments in the Gulf of Finland and comparison with the Gulf of Bothnia[J].Chemosphere,1998,36(1):43-59.
[15] 夏立江,王宏康.土壤污染及其防治[M].上海:華東理工大學出版社,2001.
[16] 錢新鋒,沈國清.蘇南河網(wǎng)地區(qū)河岸帶土壤重金屬污染生態(tài)風險評價[J].南京師大學報(自然科學版),2012,35(4):78-83.
[17] 國家環(huán)境保護局南京環(huán)境科學研究所.土壤環(huán)境質(zhì)量標準:GB 15618—1995[J].北京:中國標準出版社,1996.
[18] MULLER G.Index of geoaccumulation in sediments of the Rhine River[J].Geojournal,1969,2(3):108-118.
[19] MULLER G.Die schwermetallelstung der sedimente des necars und seiner nebenflusse:Eine bestandsaufnahme [J].Chemiker-Zeitung,1981,105:156-164.
[20] NEMEROW N L.Industrial water pollution:Origins,characteristics and treatment[M].London:Addison-Wesley Pub Co,1978.
Evaluation of Soil Heavy Metal Pollution in Southern Towns of Jiangsu Province by Multiple Methods
MENG Ling-yi1, XU Meng-jie1,LI Xiao-man1,2et al
(1. College of Public and Administration, Nanjing Agricultural University, Nanjing, Jiangsu 210095;2. Institute of Soil Science,Chinese Academy of Sciences,Nanjing,Jiangsu 210008)
Taking a small town in southern areas of Jiangsu Province as an example, the soil heavy metal pollution in research area was evaluated by using geo-accumulation index method, pollution load index method and the modified method of Nemero index combined with soil environmental quality standards. The feasibility of these methods was discussed. The results showed that the modified method of Nemerow index combined with soil environmental quality standards combined with human factors and natural factors. This method had advantages of more objective evaluation results,higher discriminabiltity of pollution distribution. And the modified method of Nemerow index had the advantages of both geo-accumulation index method and pollution load index method,so this method was most suitable for the evaluation of soil heavy metal pollution in a small scale region.
Evaluation of heavy metal pollution; Improvement method; Small scale; Southern area of Jiangsu Province
國家科技支撐計劃項目(2012BAJ24B06)。
孟令儀(1992- ),女,江蘇南京人,碩士研究生,研究方向:土地利用與土地規(guī)劃。
2016-09-19
S 158.4
A
0517-6611(2016)34-0114-05