王露+熊慶+郭彩紅+李夢(mèng)龍+蒲雪梅

摘 要 結(jié)合紫外吸收和偏最小二乘法建立了環(huán)境水樣中2，4，6三硝基甲苯及其分解物的分析方法，在不經(jīng)過(guò)任何預(yù)先分離的情況下，實(shí)現(xiàn)了環(huán)境水樣中2，4，6三硝基甲苯（TNT）及其分解物2，4二硝基甲苯（2，4DNT）和2，6二硝基甲苯（2，6DNT）的快速準(zhǔn)確測(cè)定。采用正交試驗(yàn)設(shè)計(jì)（OAD）配制了25個(gè)訓(xùn)練集樣本和15個(gè)獨(dú)立預(yù)測(cè)集樣本，采用變量種群分析（VCPA）方法挑選特征變量以提高偏最小二乘法（PLS）模型的預(yù)測(cè)效果，預(yù)測(cè)集中每個(gè)物質(zhì)的相關(guān)系數(shù)R2均達(dá)到0.99。將最優(yōu)的VCPA模型用于真實(shí)水樣中TNT以及DNT的檢測(cè)，并采用HPLC方法進(jìn)行了驗(yàn)證。實(shí)驗(yàn)結(jié)果顯示，VCPA模型的回收率與HPLC法相近。紫外光譜法結(jié)合化學(xué)計(jì)量學(xué)方法可作為環(huán)境體系中多成分同時(shí)測(cè)定的一種簡(jiǎn)便、快速、有效的方法。

關(guān)鍵詞 三硝基甲苯； 二硝基甲苯； 紫外光譜； 化學(xué)計(jì)量學(xué)； 環(huán)境水樣

1 引 言

2，4，6三硝基甲苯（TNT）在軍事和商業(yè)中[1]應(yīng)用廣泛。研究表明，TNT可引起生物體中毒[2～4]，而二硝基甲苯（DNT）是TNT的主要分解物，也是涂料、橡膠和顏料等物質(zhì)生產(chǎn)的中間物[5，6]，該類污染物在環(huán)境水中濃度可達(dá)ng/L到 mg/L[7]，長(zhǎng)時(shí)間接觸痕量的DNT會(huì)影響生物體的生長(zhǎng)和存活[8，9]。因此，檢測(cè)環(huán)境水中TNT和DNT具有重要意義。

文獻(xiàn)報(bào)道的用于檢測(cè)空氣、水和土壤中TNT及DNT含量的方法，包括熒光光譜法[10]、表面增強(qiáng)拉曼光譜法（SERS）[11，12]、氣相色譜結(jié)合電子捕獲檢測(cè)法（GCECD）[13]、氣相色譜質(zhì)譜聯(lián)用法（GCMS）[14]和高效液相色譜法（HPLC）[1，15]。這些方法選擇性和靈敏度較高，但是大多需要昂貴的儀器，樣品前處理過(guò)程復(fù)雜，耗時(shí)長(zhǎng)，因此，發(fā)展快速、簡(jiǎn)單和準(zhǔn)確的TNT及DNT定量方法具有重要的實(shí)用價(jià)值。

紫外吸收法是一種簡(jiǎn)單、快速、準(zhǔn)確的定量分析方法[16]，但選擇性不高。如果樣本中物質(zhì)相似或者未知干擾存在的情況下，吸收光譜會(huì)發(fā)生重疊，影響其準(zhǔn)確測(cè)量，因此限制了它在一些復(fù)雜體系中的應(yīng)用。化學(xué)計(jì)量學(xué)可以利用數(shù)學(xué)分離代替化學(xué)分離來(lái)解決光譜重疊的問(wèn)題[17，18]，其中偏最小二乘法是功能強(qiáng)大的多元統(tǒng)計(jì)方法，在樣本數(shù)量較少、光譜重疊的情況下，可以實(shí)現(xiàn)定量定性分析[19，20]。研究表明，通過(guò)消除一些干擾變量或者挑選一些包含重要信息的變量，可以顯著提高模型的效果[21，22]，如競(jìng)爭(zhēng)性自適應(yīng)權(quán)重取樣法（CARS）[23]、迭代偏最小二乘（iPLS）[24]、變量消除法（UVE）[25]、遺傳算法（GA）[26]、迭代保留信息變量法（IRIV）[27]。最近，Yun等[28]提出了一種新的變量挑選方法，即變量結(jié)合種群分析法（VCPA），該方法可充分利用子集之間的統(tǒng)計(jì)信息和變量間的相互關(guān)聯(lián)，在處理近紅外光譜數(shù)據(jù)時(shí)，VCPA方法比CARS、MCUVE、IRIV方法更具有優(yōu)勢(shì)。

TNT分解物中，2，4DNT和2，6DNT受到研究者的廣泛關(guān)注[29，30]，但是目前尚未有結(jié)合計(jì)量學(xué)進(jìn)行分析檢測(cè)的相關(guān)報(bào)道。本研究選擇TNT以及2，4DNT和2，6DNT作為環(huán)境水中炸藥及其分解物的代表物質(zhì)，結(jié)合化學(xué)計(jì)量學(xué)和紫外光譜法對(duì)其進(jìn)行同時(shí)定量測(cè)定。采用正交設(shè)計(jì)方法制備了一系列三組分混合物樣本，其中25個(gè)作為訓(xùn)練集，15個(gè)作為獨(dú)立預(yù)測(cè)集，此外還預(yù)測(cè)了四類真實(shí)水樣（共8個(gè)樣本）。比較了全譜和變量挑選兩種建模方法的預(yù)測(cè)結(jié)果，選取每個(gè)物質(zhì)的最優(yōu)模型用于真實(shí)水樣中目標(biāo)物的預(yù)測(cè)，并采用HPLC方法對(duì)模型的實(shí)用性和準(zhǔn)確性進(jìn)行了驗(yàn)證。研究結(jié)果表明，紫外光譜法結(jié)合變量挑選模型在復(fù)雜體系目標(biāo)物同時(shí)測(cè)定中具有良好的應(yīng)用潛力。

2 實(shí)驗(yàn)部分

2.1 儀器與試劑

Hitachi U2900型紫外分光光度計(jì)、 SHIMADZU LC20AT高效液相色譜儀（日本島津公司），色譜條件：流動(dòng)相為乙腈水（30∶70， V/V），等度洗脫，流速1.0 mL/min，柱溫35℃，檢測(cè)波長(zhǎng)254 nm，整個(gè)分離過(guò)程由CBM20Alite系統(tǒng)控制。

TNT（中國(guó)銀光化學(xué)試劑廠）； 2，4DNT和2，6NDT分析標(biāo)準(zhǔn)品（Aladdin公司上海）；乙腈（色譜純，成都科龍化學(xué)試劑廠）。實(shí)驗(yàn)用水均為超純水。標(biāo)準(zhǔn)工作溶液：準(zhǔn)確稱取適量的TNT，2，4DNT和2，6DNT固體粉末，用乙腈溶解并定容，終濃度為200 μg/mL。

2.2 紫外吸收法測(cè)定單組分硝基苯

采用單組分校正法確定紫外吸收法測(cè)定TNT和DNT的線性范圍。配制系列濃度的工作溶液。以溶劑為空白，測(cè)量各工作溶液在200～400 nm范圍內(nèi)的紫外吸收光譜。根據(jù)各分析物在最大吸收波長(zhǎng)下對(duì)應(yīng)的吸光度與樣品濃度構(gòu)建立校正曲線，得到測(cè)定TNT，2，4DNT和2，6DNT的線性范圍。

2.3 樣本配制

采用三因子五水平的正交試驗(yàn)設(shè)計(jì)（OAD）[31]配制25個(gè)訓(xùn)練集樣本，15個(gè)獨(dú)立測(cè)試樣本。為了驗(yàn)證本方法的實(shí)用性，收集了自來(lái)水（四川大學(xué)化學(xué)館）、池水（四川大學(xué)望江校區(qū)）、河水1（府南河，成都九眼橋）和涂料廠附近的河水2（江安河，雙流蛟龍港）4種水樣，分別用0.45 μm濾膜過(guò)濾。每種水樣配制兩個(gè)樣本。各樣本均加入適量TNT，2，4DNT和2，6DNT （濃度均在線性范圍內(nèi)），模擬真實(shí)環(huán)境水樣，每個(gè)樣本光譜平行測(cè)定3次。表1列出了訓(xùn)練集樣本（1～25號(hào)樣本）、獨(dú)立測(cè)試集（26～40號(hào)樣本）以及真實(shí)水樣（自來(lái)水41～42，池水43～44，河水1 45～46， 河水2 47～48）濃度配比。

2.4 化學(xué)計(jì)量學(xué)方法

偏最小二乘（PLS）算法通過(guò)最大化光譜矩陣X與目標(biāo)物質(zhì)濃度矩陣y之間的協(xié)方差為目標(biāo)，逐一求解模型中的各個(gè)潛變量。將光譜數(shù)據(jù)矩陣與樣本各組分濃度矩陣之間的內(nèi)在聯(lián)系關(guān)聯(lián)起來(lái)，使校正模型更為穩(wěn)健。在面臨樣本數(shù)目少而特征變量較多的問(wèn)題時(shí)，可以有效提取有用化學(xué)信息。

變量模型種群分析（VCPA）[28]是一種結(jié)合模型種群分析和指數(shù)遞減函數(shù)迭代進(jìn)行變量挑選的方法，將挑選出來(lái)的變量用PLS建模。具體步驟如下： （1）隨機(jī)將光譜變量建立K個(gè)子集，并對(duì)這K個(gè)子集建立K個(gè)PLS子模型； （2）統(tǒng)計(jì)K個(gè)子模型的評(píng)價(jià)參數(shù)，選取效果最好的前5%的子模型； （3）將前5%的子模型中的變量進(jìn)行權(quán)重統(tǒng)計(jì)，通過(guò)指數(shù)遞減函數(shù)強(qiáng)制刪除權(quán)重較低的變量； （4）確定剩下變量的最佳組合。

本研究采用以下指標(biāo)對(duì)模型進(jìn)行評(píng)價(jià)：校正集均方根誤差（RMSEC）、預(yù)測(cè)集均方根誤差（RMSEP）、相關(guān)系數(shù)（R2），其計(jì)算公式如下：

所有數(shù)據(jù)處理及計(jì)算在MATLAB（Version 2013a， Math work， Inc）環(huán)境下完成，PLS和VCPA算法代碼來(lái)自于中南大學(xué)中藥現(xiàn)代化研究中心。

3 結(jié)論與討論

3.1 紫外光譜表征

TNT，2，4DNT和2，6NDT的紫外吸收光譜（圖1）最大吸收波長(zhǎng)分別位于228、245和236 nm，3個(gè)物質(zhì)的吸收光譜存在嚴(yán)重的光譜重疊，在不經(jīng)過(guò)預(yù)分離的情況下，使用傳統(tǒng)校正方法難以準(zhǔn)確分析3種物質(zhì)的混合物。 因此，引入化學(xué)計(jì)量學(xué)模型，通過(guò)數(shù)學(xué)分離代替化學(xué)分離來(lái)解決光譜重疊的問(wèn)題。

按照2.2節(jié)中方法，建立了紫外光譜法分別測(cè)定單組分TNT，2，4DNT和2，6NDT的標(biāo)準(zhǔn)曲線方程，表2列出了單組分測(cè)定的一些特征參數(shù)。

表1中的參數(shù)是沒(méi)有經(jīng)過(guò)任何預(yù)處理的樣品得到的結(jié)果，雖然與文獻(xiàn)報(bào)道的一些TNT和DNT測(cè)定方法[3234]相比，檢測(cè)限和線性范圍的濃度偏高，但完全可以滿足國(guó)家標(biāo)準(zhǔn)的測(cè)定需求（國(guó)家《污水綜合排放標(biāo)準(zhǔn)》[35]中規(guī)定第二類污染物最高允許排放濃度硝基苯類為：一級(jí)標(biāo)準(zhǔn)2.0 mg/L、二級(jí)標(biāo)準(zhǔn)3.0 mg/L、三級(jí)標(biāo)準(zhǔn)5.0 mg/L）。

3.2 全譜偏最小二乘模型

為了建立可靠、準(zhǔn)確和穩(wěn)定的模型，本研究用了兩種光譜預(yù)處理方法（SG平滑，一階導(dǎo)數(shù)）處理原始光譜數(shù)據(jù)，并與未經(jīng)處理的數(shù)據(jù)進(jìn)行建模結(jié)果對(duì)比。因此，對(duì)于每個(gè)目標(biāo)物而言，共建立了3個(gè)全譜模型。五倍交叉驗(yàn)證用于評(píng)估模型的質(zhì)量。表3列出了PLS全譜模型的建模結(jié)果。

由表3可知，對(duì)于TNT和2，4DNT， 經(jīng)過(guò)SG平滑預(yù)處理的數(shù)據(jù)建模結(jié)果較好，然而對(duì)于2，6DNT而言，一階導(dǎo)數(shù)處理的數(shù)據(jù)建模效果更好。通過(guò)比較RMSEP和R2發(fā)現(xiàn)，全譜PLS建模預(yù)測(cè)TNT和2，6DNT的R2<0.98，并不理想，這可能是因?yàn)槿V數(shù)據(jù)中包含較多的無(wú)信息變量，因此進(jìn)一步采用變量種群分析（VCPA）方法來(lái)去除無(wú)信息變量，提高模型的預(yù)測(cè)結(jié)果。

3.3 變量挑選模型

同樣采用三種光譜數(shù)據(jù)（原始，SG平滑和一階導(dǎo)數(shù)）對(duì)每個(gè)物質(zhì)建立了三種VCPA模型，預(yù)測(cè)結(jié)果見(jiàn)表4。比較表3和表4的結(jié)果可以發(fā)現(xiàn)，雖然VCPA模型挑選出來(lái)的變量個(gè)數(shù)只有10個(gè)左右，但是每一個(gè)物質(zhì)預(yù)測(cè)結(jié)果的R2>0.98，高于沒(méi)有變量挑選的PLS模型，表明通過(guò)有效變量的挑選可以提高模型的準(zhǔn)確率。最后得到了每個(gè)物質(zhì)的最優(yōu)模型為：TNT，VCPA+SG平滑；2，4DNT，VCPA+SG平滑；2，6DNT，VCPA+一階導(dǎo)數(shù)。

3.4 環(huán)境水樣的檢測(cè)

為了驗(yàn)證上述最優(yōu)模型的應(yīng)用性和可靠性，將構(gòu)建的最優(yōu)模型用于4種環(huán)境水的8個(gè)樣本（每種水樣取兩個(gè)樣本）中TNT、2，4DNT和2，6DNT的濃度預(yù)測(cè)，并同時(shí)采用HPLC方法進(jìn)行驗(yàn)證。表5總結(jié)了環(huán)境水樣的全譜PLS、VCPA和HPLC方法的測(cè)定結(jié)果。

圖2A顯示了真實(shí)樣本的紫外吸收光譜圖，2A中插圖為未加入TNT，2，4DNT和2，6DNT的自來(lái)水、湖水和河水的紫外吸收光譜圖。在200～230 nm波長(zhǎng)范圍內(nèi)有較強(qiáng)的未知背景吸收，而這個(gè)區(qū)域也是3個(gè)目標(biāo)物的主要光譜信息區(qū)。這種背景干擾導(dǎo)致了全譜PLS模型的預(yù)測(cè)效果較差，3個(gè)目標(biāo)物質(zhì)的回收率分別在60.8%～155.2%、73.2%～174.6%和47.2%～141.0%。然而，盡管真實(shí)水樣中存在較大的未知干擾和目標(biāo)物的光譜重疊，VCPA建模卻可以通過(guò)消除無(wú)信息變量和干擾變量仍然取得了滿意的結(jié)果，呈現(xiàn)出強(qiáng)的抗干擾能力，其預(yù)測(cè)值與真實(shí)值和HPLC測(cè)定值非常接近，VCPA模型的平均回收率分別為97.5%、99.1%和98.6%，HPLC方法的平均回收率分別為104.7%、101.0%和102.2%。

圖2B～圖2D顯示了VCPA挑選出來(lái)的光譜點(diǎn)，其中a， b分別代表TNT，2，4DNT和2，6DNT之間的差譜[36]，藍(lán)色和紅色點(diǎn)為VCPA挑選出來(lái)的變量點(diǎn)。圖2B～圖2D的結(jié)果顯示：VCPA模型挑選出來(lái)的光譜點(diǎn)主要分布在3種分析物的強(qiáng)吸收區(qū)域220 nm～300 nm之間，是屬于苯環(huán)的ππ*躍遷的主要吸收帶區(qū)域。 差譜圖顯示， 大多數(shù)VCPA光譜點(diǎn)分布在每種分析物與其它兩種物質(zhì)的差譜峰附近，表明挑選出來(lái)的光譜點(diǎn)大多位于這些分析物差異較大的區(qū)域。雖然VCPA挑選出來(lái)的少量點(diǎn)也位于3種分析物的弱吸收區(qū)域（300～400 nm），并且差譜的差異也比較小，但這些點(diǎn)也同時(shí)位于環(huán)境水背景吸收弱的區(qū)域，并且VCPA測(cè)試結(jié)果顯示這些點(diǎn)對(duì)這3種物質(zhì)的同時(shí)定量有一定貢獻(xiàn)。上述分析證明了復(fù)雜體系中特征變量挑選的重要性，可以通過(guò)變量挑選來(lái)降低干擾，明顯提高預(yù)測(cè)結(jié)果的準(zhǔn)確性。

4 結(jié) 論

結(jié)合紫外吸收與化學(xué)計(jì)量學(xué)中的PLS方法提出了一種簡(jiǎn)單、快速、有效的分析方法，可以在不經(jīng)過(guò)任何預(yù)分離的情況下準(zhǔn)確地測(cè)定環(huán)境水樣中的TNT及其分解物（2，4DNT和2，6DNT）的含量。本研究建立了兩種模型：全譜PLS模型和變量挑選VCPA模型，用于預(yù)測(cè)獨(dú)立測(cè)試集和環(huán)境水樣中目標(biāo)物的含量，并用HPLC方法作為參照。通過(guò)對(duì)比兩種模型和HPLC的測(cè)定結(jié)果發(fā)現(xiàn)，全譜PLS模型在存在干擾的環(huán)境水樣中預(yù)測(cè)準(zhǔn)確度低，而VCPA變量挑選方法通過(guò)有效地選出重要特征變量，剔出無(wú)信息變量，所以在存在未知干擾的情況下依然可以取得準(zhǔn)確的結(jié)果，VCPA模型對(duì)每個(gè)目標(biāo)物的相關(guān)系數(shù)R2>0.99， 真實(shí)環(huán)境水樣目標(biāo)物的回收率為101.0%～104.7%，與HPLC方法的結(jié)果相近，顯示出了良好的實(shí)際應(yīng)用潛力。本研究結(jié)果表明，紫外光譜法與化學(xué)計(jì)量學(xué)方法相結(jié)合的分析策略，可以為環(huán)境中炸藥及相關(guān)物質(zhì)的定量檢測(cè)提供一種簡(jiǎn)便有效的分析方法。

References

1 Gumuscu B， Erdogan Z， Guler M O， Tekinay T. PloS one， 2014， 9（6）： e99230

2 Rieger P G， Knackmuss H J. Basic Knowledge and Perspectives on Biodegradation of 2，4，6Trinitrotoluene and Related Nitroaromatic Compounds in Contaminated Soil. Bdegradation of Nitroaromatic Compounds. Springer US， 1995： 1-18

3 Cenas N， NemeikaiteCeniene A， Sergediene E， Nivinskas H， Anusevicius Z， Sarlauskas J. BBAGen. Subjects， 2001， 1528（1）： 31-38

4 Craig H D， Taylor S. Mar.Technol. Soc. J.， 2011， 45（45）： 35-46

5 Jr Q M， Bazar M A， Mcfarland C A， Perkins， E J， Gust， K A， Gogal， R M. Ecotox. Environ. Safe.， 2007， 26（10）： 2202-2207

6 Mishra M，F(xiàn)eidler A. J. Am. Chem. Soc.， 1998， 120（24）： 6197-6198

7 MonteilRivera F， Beaulieu C， Hawari J. J.Chromatogr. A， 2005， 1066（12）： 177 -187

8 Paden N E， Smith E E， Maul J D， Kendall R J. Ecotox. Environ. Safe.， 2011， 74（4）： 924-928

9 Podlipn R， Pospíilov B， Vaněk T. Ecotox.Environ. Safe.， 2015， 112： 54-59

10 Diaz D G B， Moreno D， Torroba T， Berg A， Gunnars J， Nilsson T. J. Am. Chem. Soc.， 2014， 136（11）： 4125-4128

11 Li J F， Huang Y F， Ding Y， Yang Z L， Li S B， Zhou X S， Fan F R. Nature， 2010， 464： 392-395

12 Gong Z， Du H， Cheng F，Wang C， Fan M. ACS Appl. Mater. Interfaces， 2014， 6（24）： 21931-21937

13 Walsh M E. Talanta， 2001， 54（3）： 427-438

14 Yinon J. J. Chromatogr. A， 1996， 742（12）： 205-209

15 RahimiNasrabadi M， Zahedi M M， Pourmortazavi S M， Heydari R， Rai H， Jazayeri J， Javidan A. Microchim. Acta， 2012， 177（1）： 145-152

16 Perkampus H H. UVVIS Spectroscopy and Its Applications. Springer， Berlin Heidelberg， 1992

17 Gu H W， Wu H L， Yin X L， Li Y， Liu Y J， Xia H. Anal. Chim. Acta， 2014， 848： 10-24

18 Thomas EV，Haaland D M. Anal. Chem.， 1990， 62（10）： 1091-1099

19 Forina M， Lanteri S， Oliveros M C， Millan C P. Anal. Bioanal. Chem.， 2004， 380（3）： 397-418

20 Haddad J E，Miollis F D， Sleiman J B， Canioni L， Mounaix P， Bousquet B. Anal. Chem.， 2014， 86（10）： 4927-4933

21 Mehmood T， Liland K H， Snipen L. Chemom. Intell. Lab. Syst.， 2012， 118（16）： 62-69

22 Shao X， Bian X， Cai W. Anal. Chim. Acta， 2010， 666（12）： 32-37

23 Li H D， Liang Y Z，Xu Q S， Cao D S. Anal. Chim. Acta， 2009， 648（1）： 77-84

24 Lindgren F，Geladi P， Rnnar S， Wold S. J. Chemometr.， 1994， 8（5）： 349-363

25 Centner V， Massart D L， Noord O E D， Jong S d， Vandeginste B M， Sterna C. Anal. Chem.， 1996， 68（21）： 3851-3858

26 Leardi R. J. Chemometr.， 2000， 14（56）： 643-655

27 Yun Y H， Wang W T， Tan M L， Liang Y Z， Li H D， Cao D S， Lu H M， Xu Q S. Anal. Chim. Acta， 2014， 807（1）： 36-43

28 Yun Y H， Wang W T， Deng B C， Lai G B， Liu X B， Ren D B， Liang Y Z， Fan， Xu Q S. Anal. Chim. Acta， 2015， 862： 14-23

29 Wright L K， Zellers ET. Analyst， 2013， 138（2）： 6860-6868

30 Snels M， Venezia T， Belfiore L. Chem. Phys. Lett.， 2010， 489（13）： 134-140

31 Zhu G， Ju H. Anal. Chim. Acta， 2004， 506（2）： 177-181

32 HJ 6002011， Water QualityDetermination of TNT、RDX、DNT Gas Chromatography. The State Environmental Protection Standards of the People′s Republic of China

水質(zhì)梯恩梯、黑索金、地恩梯的測(cè)定氣相色譜法. 中華人民共和國(guó)國(guó)家環(huán)境保護(hù)標(biāo)準(zhǔn). HJ 6002011

33 HJ 5992011， Water QualityDetermination of TNTNCetylpridinium ChlorideSodium Sulfite Spectrophotometric Method. The State Environmental Protection Standards of the People′s Republic of China

水質(zhì)梯恩梯的測(cè)定 N氯代十六烷基吡啶—亞硫酸鈉分光光度法. 中華人民共和國(guó)國(guó)家環(huán)境保護(hù)標(biāo)準(zhǔn). HJ 5992011

34 HJ 5982011， Water QualityDetermination of TNTSodium Sulfite Spectrophotometric Method. The State Environmental Protection Standards of the People′s Republic of China

水質(zhì)梯恩梯的測(cè)定亞硫酸鈉分光光度法. 中華人民共和國(guó)國(guó)家環(huán)境保護(hù)標(biāo)準(zhǔn). HJ 5992011

35 GB 89781996， Integrated Wastewater Discharge Standard. National Standards of the People′s Republic of China.

污水綜合排放標(biāo)準(zhǔn). 中華人民共和國(guó)國(guó)家標(biāo)準(zhǔn). GB 89781996

36 Xu D， Fan W， Lyu H， Liang Y， Shan Y， Li G， Yang Z， Yu L. Spectrochim. Acta A， 2014， 123（7）： 430-435