張旭

摘要[目的]探討土壤中重金屬Pb消解的影響因素。[方法]利用二維平面吸光度等值線變化探討區(qū)域范圍內(nèi)吸光度趨勢變化，繪制Pb吸光度分布圖，用Minitab 16軟件對試驗(yàn)結(jié)果進(jìn)行極差分析和方差分析。[結(jié)果]影響Pb消解因素主次排序?yàn)镠F、HNO3、HClO4、燃?xì)饬髁俊Cl、助燃?xì)饬髁俊綦娏骱腿紵鞲叨?。?yōu)化試驗(yàn)表明，在土壤樣品中加入9 mL HNO3、5 mL HF、3 mL HClO4經(jīng)微波消解后，原子吸收分光光度計(jì)工作條件如下：Pb波長為217.0 nm，狹縫寬度為1.3 nm，燃燒器高度為4 nm，助燃?xì)饬髁繛?.4 L/min，燈電流為4.0 mA，燃?xì)饬髁繛?.0 L/min。在上述條件下，最優(yōu)吸光度為0.302。[結(jié)論]研究結(jié)果可為精確測定土壤中各種重金屬吸光度提供理論依據(jù)。

關(guān)鍵詞混合正交優(yōu)化；極差分析；方差分析；重金屬； Pb

中圖分類號X833文獻(xiàn)標(biāo)識碼

A文章編號0517-6611（2017）22-0009-04

Abstract[Objective] The aim was to explore the affecting factors for the digestion of soil heavy metal Pb. [Method] The absorbance trend within the area was expressed by change of twodimensional planar contour for the first time. We used Minitab 16 software to analyze the test results by range and variance method. [Result] The factors affecting Pb digestion was in order of HF， HNO3， HClO4， gas flow， HCl， assistant gas flow， lamp current and combustion height. Test results showed that after adding 9 mL HNO3，5 mL HF and 3 mL HClO4 in soil samples， and microwave digestion， the working conditions of atomic absorption spectrophotometer were as followed： Pb wavelength 217.0 nm， slit width 1.3 nm， burner height 4 nm， assistant combustion gas rate 2.4 L/min， the lamp current 4.0 mA， and propane flow rate 4.0 L/min. Under the above optimal conditions， the optimal absorbance value was 0.302. [Conclusion] The results provide theoretical basis for precisely determining absorbance values of heavy metals in soil.

Key wordsMixed orthogonal optimization；Range method；Variance method；Heavy metal；Pb

正交試驗(yàn)設(shè)計(jì)是由日本質(zhì)量管理專家田口玄一博士在多因素試驗(yàn)設(shè)計(jì)基礎(chǔ)上研究開發(fā)出來的一種新試驗(yàn)設(shè)計(jì)方法。該法屬于數(shù)理統(tǒng)計(jì)方法，利用統(tǒng)計(jì)學(xué)思維考慮自然科學(xué)問題，在醫(yī)學(xué)、化工領(lǐng)域應(yīng)用廣泛。正交試驗(yàn)設(shè)計(jì)優(yōu)點(diǎn)在于合理安排正交試驗(yàn)表，既可減少試驗(yàn)數(shù)量又可保持試驗(yàn)均衡性和準(zhǔn)確性[1]，對于正交試驗(yàn)多因素、多水平問題，采用若干次試驗(yàn)即可找出各因素主次關(guān)系以及因素之間最優(yōu)搭配，還能得到最佳試驗(yàn)方案以及最優(yōu)試驗(yàn)結(jié)果。筆者以正交試驗(yàn)設(shè)計(jì)為基礎(chǔ)，將規(guī)范的正交試驗(yàn)改造為混合正交試驗(yàn)做進(jìn)一步優(yōu)化，將正交試驗(yàn)設(shè)計(jì)所得結(jié)果轉(zhuǎn)化到連續(xù)的區(qū)域平面內(nèi)，以此確定最優(yōu)數(shù)值，系統(tǒng)闡述了該法在環(huán)境試驗(yàn)優(yōu)化中的應(yīng)用。

1材料與方法

1.1材料

1.1.1土壤樣品。土壤樣品來自于鉛鋅礦區(qū)，經(jīng)過風(fēng)干、研磨后過160目尼龍篩，置于恒溫干燥箱內(nèi)，在105 ℃下干燥去水，置于干燥器中備用[2]。土壤成分組成復(fù)雜，由于傳統(tǒng)方法不能徹底分解土壤中殘?jiān)镔|(zhì)，使試驗(yàn)存在較大偏差[3]。而使用HF、HClO4、HNO3等混合酸消解體系雖獲得良好效果，但其消解條件、加酸比例以及儀器測定過程仍存在誤差，該部分誤差可通過優(yōu)化降到最低。利用混合正交試驗(yàn)設(shè)計(jì)可優(yōu)化加酸比例和校正原子吸收分光光度計(jì)測定條件，使試驗(yàn)獲得更高精度。

1.1.2儀器與試劑。

TAS-986原子吸收分光光度計(jì)，優(yōu)級純Pb粉；優(yōu)級純硝酸（HNO3，ρ=1.51 g/mL）；分析純鹽酸（HCl， ρ=1.19 g/mL），優(yōu)級純高氯酸（HCIO4，ρ=1.68 g/mL）；分析純氫氟酸（HF，ρ=1.49 g/mL）；10%磷酸氫二銨[（NH4）2HPO4]水溶液。試驗(yàn)所用玻璃和塑料器皿均由（1+1）HNO3浸泡5 h，經(jīng)KM-23C超聲波清洗機(jī)清洗后，用去離子水清洗3次，置于恒溫干燥箱中干燥備用。消解儀器為MDS-6溫壓雙控微波消解儀。

1.2方法

1.2.1樣品測定。取土壤樣品0.5 g（精確到±0.000 1 g），加入聚四氟乙烯容樣杯中[4]，加酸量及比例關(guān)系依據(jù)表1、2試驗(yàn)安排進(jìn)行，消解過程依據(jù)表3運(yùn)行。由于HClO4與土壤中有機(jī)質(zhì)反應(yīng)產(chǎn)生脂類化合物在高溫條件下極易爆炸，故消解過程采用低壓消解[5]。微波消解運(yùn)行完畢，降壓、冷卻、開蓋轉(zhuǎn)移液體，定容。測定過程為設(shè)定火焰原子吸收分光光度計(jì)工作條件，調(diào)節(jié)燃燒器高度、燃?xì)饬髁?、燈電流、助燃?xì)饬髁繙?zhǔn)備測樣，在測定試驗(yàn)樣品前滴加2滴磷酸氫二銨水溶液[6]。在Pb波長λPb為217.0 nm、狹縫寬度為1.3 nm、Pb測定濃度為0.20～10.00 mg/L、火焰類型為乙炔-空氣的條件下測定。

1.2.2混合正交試驗(yàn)設(shè)計(jì)。試驗(yàn)因素分為定性因素和定量因素，定量因素均以具體數(shù)值變化來控制試驗(yàn)。水平數(shù)可以用來控制試驗(yàn)范圍和試驗(yàn)次數(shù)，由于各因素對試驗(yàn)影響程度有差異，故對不同因素設(shè)定不同水平數(shù)，對試驗(yàn)影響較大因素應(yīng)設(shè)較多水平，以此來精確定位最佳條件，設(shè)A（HNO3）、B（HF）、C（HClO4）為4個(gè)水平；HCl對土壤重金屬消解作用效果不是十分明顯[7]，故設(shè)D（HCl）為2個(gè)水平。儀器運(yùn)行參數(shù)F（燃燒器高度）、G（燃?xì)饬髁浚?、H（燈電流）、I（助燃?xì)饬髁浚┰O(shè)2個(gè)水平，如表1所示。

由于試驗(yàn)中各因素之間相互獨(dú)立，因素之間不存在交互作用，表頭可以任意設(shè)計(jì)，但要滿足因素組合均勻、分散。根據(jù)因素水平數(shù)計(jì)算自由度，安排正交表。總自由度計(jì)算為各因素與其對應(yīng)水平數(shù)乘積之和，具體計(jì)算：fT=fA+fB+fC+fD+fE+fH+fI=（4-1）×3+（2-1）×6 =15。分析自由度可知，采用L16（215）正交表，其中L是正交表代號，16是試驗(yàn)次數(shù)，2代表水平數(shù)，15是正交表列數(shù)即因素?cái)?shù)。

試驗(yàn)中由于各因素水平數(shù)不同，需要經(jīng)并列法改造為3因素2水平和6因素2水平的L16（43×26）混合正交表。以往設(shè)計(jì)通常習(xí)慣套用正交表，忽略空白列重要性，空白列是用來衡量試驗(yàn)可信程度，與試驗(yàn)做對比分析，若沒有空白列，則無法保證試驗(yàn)準(zhǔn)確程度。經(jīng)并列法改造后的正交試驗(yàn)表變?yōu)榛旌险辉囼?yàn)表，并設(shè)E列為空白列，使每組試驗(yàn)都能與空白值做對比。更重要的是，經(jīng)改造后的混合正交試驗(yàn)表仍然保持均衡分散、整齊可比。

1.2.3數(shù)據(jù)分析。將16組試驗(yàn)在Suffer軟件中建立坐標(biāo)，記錄Pb吸光度并與其所屬組對應(yīng)，設(shè)定繪圖區(qū)域面積為50 cm×80 cm，繪制重金屬Pb二維吸光度等值線變化圖，將16組Pb吸光度轉(zhuǎn)化為平面表示，利用顏色變化判斷吸光度大小變化。再利用三維可視化技術(shù)實(shí)現(xiàn)重金屬Pb吸光度立體趨勢走向，通過峰值高低和顏色變化判斷混合正交試驗(yàn)最優(yōu)點(diǎn)走向分布，反之，通過觀察峰值走向可進(jìn)一步預(yù)測最優(yōu)點(diǎn)所在條件范圍。確定條件所屬范圍后，通過正交試驗(yàn)依次循環(huán)做進(jìn)一步優(yōu)化。該試驗(yàn)以Pb為例只做1次優(yōu)化，然后利用極差分析法和方差分析法分析試驗(yàn)結(jié)果。

2結(jié)果與分析

2.1極差分析

在Excel中排列數(shù)據(jù)并分析計(jì)算極差值R如表4所示，并用Minitab 16軟件作因素指標(biāo)趨勢圖。如圖1所示，其中i（i=1，2，3，4）為某因素在其某一水平下參于試驗(yàn)所得Pb吸光度平均值。極差R=Kmax-Kmin，以此反映各因素水平變動(dòng)對試驗(yàn)測定結(jié)果吸光度大小的影響[5]。極差R越大則影響越大，為主要因素，反之為次要因素。通過觀察二維、三維16組Pb吸光度變化趨勢，可知在這16組數(shù)據(jù)中Pb吸光度最優(yōu)數(shù)值出現(xiàn)在第12組，即在A3（7 mL HNO3）、B4（5 mL HF）、C2（4 mL HClO4）、D1（HCl）、F2（燃燒器高度6 nm）、G1（助燃?xì)饬髁?.0 L/min）、H1（燈電流3.0 mA）、I2（燃?xì)饬髁?.0 L/min）條件下測得最大吸光度（0.272），但其不一定是最優(yōu)數(shù)值。通過極差分析法判斷影響試驗(yàn)因素主次并且驗(yàn)證未出現(xiàn)在混合正交試驗(yàn)表中的試驗(yàn)組是否可能出現(xiàn)最大值。

分析表4極差R值大小，可知對試驗(yàn)影響的主次因素依次為B（HF）、C（HNO3）、A（HClO4）、I（燃?xì)饬髁浚（HCl）、G（助燃?xì)饬髁浚?、H（燈電流）、F（燃燒器高度）、E（空白列）。其中，HF、HClO4、HNO3對Pb吸光度影響較大。從因素主次順序可知，燃?xì)饬髁繉υ囼?yàn)效果影響顯著，原因可能是乙炔燃?xì)饬魉俅笮Q定乙炔燃燒是否充足，以及焰心溫度是否滿足重金屬測定所需溫度，假若溫度不足，無法保證待測元素Pb離解為基態(tài)原子，也就不能有效進(jìn)行特征輻射光波吸收[8]。而G（助燃?xì)饬髁浚?、H（燈電流）、F（燃燒器高度）主次排序較后，說明微調(diào)火焰原子吸收儀中助燃?xì)?、燈電流、燃燒器高度對于Pb吸光度測定影響較小。

圖2為因素A、B、C、D、E、F、G、H、I在不同水平條件下得到的Pb吸光度。取各因素下Pb最大吸光度，匯總可得試驗(yàn)組為A4B4C1D1E2F1G2H2I1，測定土壤在該組條件下的吸光度，即在土壤樣品加入9 mL HNO3、3 mL HF、3 mL HClO4消解后，在火焰原子吸收儀工作條件為波長λPb217.0 nm、f狹縫寬度1.3 nm、燃燒器高度4 nm、助燃?xì)饬髁?.4 L/min、燈電流4.0 mA、燃?xì)饬髁?.0 L/min條件下測得Pb吸光度為0.276。該組有可能是最優(yōu)試驗(yàn)組，但該組并不在16組試驗(yàn)中，此時(shí)需要通過測定A4B4C1D1E2F1G2H2I1條件下的吸光度，并與第12組A3B4C2D1E1F2G1H1I2得到的吸光度（0.276）進(jìn)一步對比。通過試驗(yàn)測得該組吸光度為0.304。由圖2可知，在A4B4C1D1E2F1G2H2I1條件下吸光度達(dá)到峰值（0.304），而在混合正交試驗(yàn)表中第12組條件下測得最優(yōu)吸光度只有0.276，顯然，A4B4C1D1E2F1G2H2I1為最優(yōu)試驗(yàn)組，這正體現(xiàn)了混合正交設(shè)計(jì)優(yōu)勢。

2.2方差分析

方差分析法彌補(bǔ)了極差分析法的不足，能夠在誤差干擾情況下做出準(zhǔn)確判斷[9]?；旌险辉囼?yàn)方差計(jì)算在本質(zhì)上與等水平正交試驗(yàn)設(shè)計(jì)方差計(jì)算相同，區(qū)別在于混合正交試驗(yàn)設(shè)計(jì)在計(jì)算過程中水平數(shù)有變動(dòng)[10]。利用Minitab 16軟件計(jì)算離差平方和、自由度、平均離差平方和以及顯著性，其中顯著性檢驗(yàn)強(qiáng)調(diào)因素在分析過程中對試驗(yàn)的影響作用，如果因素對試驗(yàn)影響較小，討論因素隨試驗(yàn)變化也毫無意義。數(shù)據(jù)分析結(jié)果如表5所示。

由方差分析表5中F值可知，A（HNO3）、B（HF）、C（HClO4）、I（燃?xì)饬髁浚υ囼?yàn)結(jié)果影響較大，顯著性較強(qiáng)。試驗(yàn)中各因素顯著性大小為B（HF）>C（HClO4）>A（HNO3）>I（燃?xì)饬髁浚?F（燃燒器高度）、G（助燃?xì)饬髁浚（燈電流）、E（空白列）。對極差分析法中R值與方差分析法中F值繪柱狀圖做對比分析。

3結(jié)論

采用混合正交試驗(yàn)設(shè)計(jì)可以了解因素與水平變動(dòng)的規(guī)律性以及各因素之間是否存在交互作用，可迅速找出最優(yōu)的生產(chǎn)條件或工藝條件，并能預(yù)測在最優(yōu)生產(chǎn)條件下試驗(yàn)結(jié)果的趨勢。

試驗(yàn)結(jié)果表明，在土壤樣品中加入9 mL HNO3、5 mL HF、3 mL HClO4，經(jīng)微波消解后，火焰原子吸收儀工作條件如下：波長λPb為217.0 nm，狹縫寬度為1.3 nm，燃燒高度為4 nm，助燃?xì)饬髁繛?.4 L/min，燈電流為4.0 mA，燃?xì)饬髁繛?.0 L/min。在上述條件下，Pb獲得最優(yōu)吸光度（0.302）。

極差分析法與方差分析法相結(jié)合使得混合正交試驗(yàn)數(shù)據(jù)分析中存在的缺陷得到解決，并且提高了試驗(yàn)精度。

參考文獻(xiàn)

[1]

張晟，陳玉成.環(huán)境試驗(yàn)優(yōu)化設(shè)計(jì)與數(shù)據(jù)分析[M].北京：化學(xué)工業(yè)出版社，2008：5-17.

[2] 鄧華，許丹丹，李明順，等.不同消解方法分析土壤中重金屬含量的比較[J].廣西師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2010，28（3）：80-83.

[3] 但德忠.分析測定中的現(xiàn)代微波制樣技術(shù)[M].成都：四川大學(xué)出版社，2003：45-46.

[4] 夏家淇，蔡道基，夏增祿.土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)：GB 15618—1995[S].北京：中國環(huán)境科學(xué)出版社，1995.

[5] 冷庚，楊嘉偉，謝晴，等.微波消解-氫化物發(fā)生原子熒光光度法測定土壤中的汞[J].環(huán)境工程學(xué)報(bào)，2011，5（8）：1893-1896.

[6] 吳宏，黃德乾，金焰，等.環(huán)境樣品中鉛、銻、汞、硒形態(tài)分析研究進(jìn)展[J].環(huán)境監(jiān)測管理與技術(shù)，2008，20（4）：9-17.

[7] 國家環(huán)保局編委會(huì).水和廢水監(jiān)測分析方法[M].4版.北京：中國環(huán)境科學(xué)出版社，2002.

[8] 史嘯勇，郁建橋. 微波消解-原子吸收光度法測定土壤中銅鋅鉛鎘鎳鉻[J]. 環(huán)境監(jiān)測管理與技術(shù)，2003，15（1）：32-33.

[9] 滕海因，祝國強(qiáng)，黃平，等.正交試驗(yàn)設(shè)計(jì)實(shí)例分析[J].藥學(xué)服務(wù)與研究，2007，8（1）：75-76.

[10] 劉宏，王俊偉，靳皓琛，等.正交試驗(yàn)-硝酸微波消解-冷原子吸收法測定土壤中汞[J].四川環(huán)境，2012，31（5）：41-44.