石娟娟　趙艷玲　何廳廳等

摘要：以某金礦區(qū)周邊土壤重金屬Pb、Cu為研究對(duì)象，采用GIS、地統(tǒng)計(jì)學(xué)相結(jié)合的方法對(duì)該研究區(qū)土壤重金屬的空間結(jié)構(gòu)及變異規(guī)律進(jìn)行研究。結(jié)果表明：Pb、Cu含量均值分別為713.358、108.828 mg/kg，兩者的變異系數(shù)均大于1，屬于強(qiáng)變異。Pb、Cu的空間分布相關(guān)系數(shù)為0.978，達(dá)到極顯著相關(guān)水平，表現(xiàn)出極強(qiáng)的協(xié)同作用。Pb、Cu的含量分布均不符合正態(tài)分布，經(jīng)對(duì)數(shù)變換后，K-S值大于0.05，且P-P圖上2種重金屬的特征數(shù)據(jù)均呈現(xiàn)直線趨勢(shì)。Pb、Cu存在半方差結(jié)構(gòu)，且二者的最優(yōu)擬合模型分別為指數(shù)模型、球狀模型。Pb、Cu的變程分別為450、236 m，Pb的變異程度強(qiáng)于Cu。二者的塊基比分別為0.51×10-3、0.11×10-2，均小于0.25，說明二者有很強(qiáng)的空間自相關(guān)性。Pb含量呈團(tuán)狀分布，Cu含量呈條帶狀分布，均表現(xiàn)出東南高、西北低的特征。

關(guān)鍵詞：空間結(jié)構(gòu)；空間插值；空間自相關(guān)性；重金屬污染

中圖分類號(hào)： X53；S127；X753 文獻(xiàn)標(biāo)志碼： A 文章編號(hào)：1002-1302（2014）07-0373-03

收稿日期：2013-10-24

基金項(xiàng)目：教育部新世紀(jì)優(yōu)秀人才支持計(jì)劃（編號(hào)：NCET-12-0964）。

作者簡(jiǎn)介：石娟娟（1987—），女，河南駐馬店人，碩士研究生，主要從事國土資源遙感研究。E-mail：sjj090312@163.com。土壤重金屬容易通過揚(yáng)塵、植物、人類活動(dòng)等途徑進(jìn)入人體，從而嚴(yán)重影響人類健康。除此之外，重金屬污染還具有環(huán)境危害持久性、地球化學(xué)循環(huán)性、生態(tài)風(fēng)險(xiǎn)性等[1-3]。礦山開采是重金屬污染的重要來源，礦山開發(fā)過程中產(chǎn)生的大量廢棄物如尾礦、廢石隨處堆放，不僅占用大量土地，而且破壞堆置場(chǎng)原有的生態(tài)環(huán)境[4-5]。研究礦區(qū)重金屬的空間結(jié)構(gòu)及分布特征對(duì)于制定重金屬污染修復(fù)方案、提高土壤利用率具有重要意義。重金屬在土壤中的分布是不均一的、連續(xù)變化的，具有高度的空間變異性。傳統(tǒng)的統(tǒng)計(jì)分析方法通常將土壤看作是一個(gè)獨(dú)立的均質(zhì)區(qū)域，通過土壤的統(tǒng)計(jì)數(shù)據(jù)來描述土壤的空間變異，因此存在較多局限性。地統(tǒng)計(jì)學(xué)能定量描述區(qū)域化變量空間的變異特征，克立格插值能可視化土壤特性空間分布格局，并為其空間預(yù)測(cè)提供最優(yōu)無偏估計(jì)[6]。Pb、Cu是土壤重金屬空間變異研究的重要對(duì)象。劉瓊峰等對(duì)城郊農(nóng)田土壤Pb、Cd的空間變異性研究發(fā)現(xiàn)，城市中人類活動(dòng)在一定空間尺度內(nèi)對(duì)近郊區(qū)農(nóng)田土壤Pb、Cd含量有較大的影響，通過半方差函數(shù)分析，發(fā)現(xiàn)空間變異由結(jié)構(gòu)性因素、隨機(jī)因素共同作用[7]。白曉宇等對(duì)銅陵礦區(qū)土壤重金屬元素的空間變異研究發(fā)現(xiàn)，銅陵礦區(qū)土壤重金屬As、Cd、Pb、Zn元素的變異函數(shù)表現(xiàn)為各向異性，其方向性主要受礦床分布控制[8]。本研究以金礦區(qū)周邊土壤重金屬Pb、Cu為材料，利用GIS、地統(tǒng)計(jì)學(xué)分析方法對(duì)該礦區(qū)Pb、Cu的空間結(jié)構(gòu)及變異規(guī)律進(jìn)行分析，以期為該礦區(qū)Pb、Cu污染治理提供依據(jù)。

1材料與方法

1.1研究區(qū)概況

該研究區(qū)屬暖溫帶大陸性雨熱同季的季風(fēng)性干旱氣候，溫差較大，降水較少，氣候比較干燥。受成土母質(zhì)、外界環(huán)境影響，該區(qū)的土壤結(jié)構(gòu)和類型差別較大，由南向北大體可分為棕壤、褐土、沙黏土、壚土4種土壤。

1.2樣品采集與分析

采用網(wǎng)格法布設(shè)采樣點(diǎn)，采集耕作層（0～20 cm）的土壤樣品，除去雜草、草根、礫石等雜物。每采樣點(diǎn)以1點(diǎn)為中心，在周圍50 m范圍內(nèi)采集3處子樣，通過混合四分法，保留 1 kg 樣品，裝入布袋。采樣、樣品保存、樣品處理過程中均采用非金屬容器，避免樣品污染。布設(shè)土壤采樣點(diǎn)時(shí)綜合考慮研究區(qū)成土母質(zhì)、土壤類型、地形條件等因素，采用網(wǎng)格法、分層抽樣法進(jìn)行設(shè)計(jì)。野外樣品采集過程中利用1 ∶50 000地形圖、 GPS進(jìn)行定位，并根據(jù)實(shí)地情況進(jìn)行調(diào)整，共獲得46個(gè)土壤樣點(diǎn)（圖1）。

采集的土壤樣品經(jīng)自然風(fēng)干后，用木棒研碎，過20目篩。室內(nèi)用高鋁缽粉碎均勻，待樣品粒度達(dá)到一定程度后過100目篩。采用等離子體質(zhì)譜法（ICP-MS）測(cè)定土壤樣品中Pb、Cu 2種重金屬的含量。測(cè)定結(jié)果如表1所示。

1.3方法

地統(tǒng)計(jì)學(xué)也被稱為地質(zhì)統(tǒng)計(jì)學(xué)，是以具有空間分布特點(diǎn)的區(qū)域化變量理論為基礎(chǔ)，以半方差函數(shù)為基本工具，運(yùn)用克里格插值法，對(duì)自然現(xiàn)象的空間變異問題進(jìn)行分析研究[9]。本研究采用SPSS19.0軟件對(duì)Pb、Cu的空間含量分布進(jìn)行正態(tài)分布檢驗(yàn)并剔除異常值，采用GS+9.0進(jìn)行模型最優(yōu)擬合，最后在ArcGIS9.3軟件中應(yīng)用Kriging插值模塊，輸入擬合參數(shù)，內(nèi)插生成土壤中Pb、Cu的空間含量分布圖。

3結(jié)論

本研究結(jié)果表明，Pb、Cu含量均值分別為713.358、108.828 mg/kg，均超過了國家土壤環(huán)境二級(jí)標(biāo)準(zhǔn)。兩者的變異系數(shù)均大于1，屬于強(qiáng)變異。Pb、Cu的空間分布相關(guān)系數(shù)為0.978，達(dá)到極顯著相關(guān)水平，表現(xiàn)出極強(qiáng)的協(xié)同作用，而非拮抗。Pb、Cu的含量分布均不符合正態(tài)分布，經(jīng)對(duì)數(shù)變換后，K-S值大于0.05，且P-P圖上2種重金屬的特征數(shù)據(jù)均呈現(xiàn)直線趨勢(shì)，說明對(duì)數(shù)轉(zhuǎn)換后的數(shù)據(jù)符合正態(tài)分布，可用于地統(tǒng)分析。Pb、Cu存在半方差結(jié)構(gòu)，且二者的最優(yōu)擬合模型分別為指數(shù)模型、球狀模型。Pb、Cu的變程分別為450、236 m，Pb的變異程度強(qiáng)于Cu。二者的塊基比分別為 0.51×10-3、0.11×10-2，均小于0.25，說明二者有很強(qiáng)的空間自相關(guān)性，表明Pb、Cu的含量分布主要受人為活動(dòng)的影響。Pb含量呈團(tuán)狀分布，Cu含量呈條帶狀分布，均表現(xiàn)出東南高、西北低的特征，高值區(qū)均分布在東南角的山前沖積斜塬地區(qū)，且空間變異主要發(fā)生在東南到西北的方向。Pb的空間分布圖上圖斑數(shù)目要明顯多于Cu，這與Pb的強(qiáng)變異性密切相關(guān)。

參考文獻(xiàn)：endprint

[1]史貴濤，陳振樓，李海雯，等. 城市土壤重金屬污染研究現(xiàn)狀與趨勢(shì)[J]. 環(huán)境監(jiān)測(cè)管理與技術(shù)，2006，18（6）：9-12，24.

[2]于小俸，唐嶷，胡玉. 酸性土壤重金屬含量的影響因素——以湖南壺瓶山自然保護(hù)區(qū)為例[J]. 廣東農(nóng)業(yè)科學(xué)，2011，38（9）：54-55.

[3]Lu J S，Liu Y，Zhang Z L，et al. Factorial kriging and stepwise regression approach to identify environmental factors influencing spatial multi-scale variability of heavy metals in soils[J]. Journal of Hazardous Materials，2013，261：387-397.

[4]毛竹，張世熔，李婷，等. 鉛鋅礦區(qū)土壤重金屬空間變異及其污染風(fēng)險(xiǎn)評(píng)價(jià)——以四川漢源富泉鉛鋅礦山為例[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2007，26（2）：617-621.

[5]何佳芳，何騰兵. 某鉛鋅冶煉廠廢棄地復(fù)墾整理區(qū)土壤重金屬污染評(píng)價(jià)[J]. 水土保持學(xué)報(bào)，2006，20（2）：97-101.

[6]張博，趙耕毛，劉兆普，等. 江蘇灘涂圍墾區(qū)土壤養(yǎng)分空間變異研究[J]. 江蘇農(nóng)業(yè)科學(xué)，2010（5）：461-464.

[7]劉瓊峰，李明德，吳海勇，等. 城郊農(nóng)田土壤Pb、Cd的空間變異與評(píng)價(jià)研究——以長沙市為例[J]. 長江流域資源與環(huán)境，2012，21（2）：195-203.

[8]白曉宇，袁峰，李湘凌，等. 銅陵礦區(qū)土壤重金屬元素的空間變異及污染分析[J]. 地學(xué)前緣，2008，15（5）：256-263.

[9]黃勇，郭慶榮，任海，等. 地統(tǒng)計(jì)學(xué)在土壤重金屬研究中的應(yīng)用及展望[J]. 生態(tài)環(huán)境，2004，13（4）：681-684.

[10]吳文勇，尹世洋，劉洪祿，等. 污灌區(qū)土壤重金屬空間結(jié)構(gòu)與分布特征[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2013，29（4）：165-173.

[11]姚榮江，楊勁松，楊奇勇. 典型農(nóng)區(qū)耕作土壤重金屬空間變異的穩(wěn)健性分析[J]. 環(huán)境科學(xué)與技術(shù)，2012，35（S2）：1-7.

[12]Wang X J，Zhang Z P. A comparison of conditional simulation，kriging and trend surface analysis for soil heavy metal pollution pattern analysis[J]. Journal of Environmental Science and Health，Part A：Toxic/Hazardous Substances and Environmental Engineering，1999，34（1）：73-89.

[13]孫波，周生路，趙其國. 基于空間變異分析的土壤重金屬復(fù)合污染研究[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2003，22（2）：248-251.

[14]蔡德所，王魁，黃景新，等. 喀斯特峰叢洼地原生林區(qū)土壤礦質(zhì)元素空間異質(zhì)性研究[J]. 中國水土保持，2010（5）：33-36.

[15]葛劍平，郭海燕，仲莉娜. 地統(tǒng)計(jì)學(xué)在生態(tài)學(xué)中的應(yīng)用（Ⅰ）——基本理論和方法[J]. 東北林業(yè)大學(xué)學(xué)報(bào)，1995，23（2）：88-94.

[16]Zhao Y C，Wang Z G，Sun W X，et al. Spatial interrelations and multi-scale sources of soil heavy metal variability in a typical urban-rural transition area in Yangtze River Delta region of China[J]. Geoderma，2010，156（3/4）：216-227.endprint

[1]史貴濤，陳振樓，李海雯，等. 城市土壤重金屬污染研究現(xiàn)狀與趨勢(shì)[J]. 環(huán)境監(jiān)測(cè)管理與技術(shù)，2006，18（6）：9-12，24.

[2]于小俸，唐嶷，胡玉. 酸性土壤重金屬含量的影響因素——以湖南壺瓶山自然保護(hù)區(qū)為例[J]. 廣東農(nóng)業(yè)科學(xué)，2011，38（9）：54-55.

[3]Lu J S，Liu Y，Zhang Z L，et al. Factorial kriging and stepwise regression approach to identify environmental factors influencing spatial multi-scale variability of heavy metals in soils[J]. Journal of Hazardous Materials，2013，261：387-397.

[4]毛竹，張世熔，李婷，等. 鉛鋅礦區(qū)土壤重金屬空間變異及其污染風(fēng)險(xiǎn)評(píng)價(jià)——以四川漢源富泉鉛鋅礦山為例[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2007，26（2）：617-621.

[5]何佳芳，何騰兵. 某鉛鋅冶煉廠廢棄地復(fù)墾整理區(qū)土壤重金屬污染評(píng)價(jià)[J]. 水土保持學(xué)報(bào)，2006，20（2）：97-101.

[6]張博，趙耕毛，劉兆普，等. 江蘇灘涂圍墾區(qū)土壤養(yǎng)分空間變異研究[J]. 江蘇農(nóng)業(yè)科學(xué)，2010（5）：461-464.

[7]劉瓊峰，李明德，吳海勇，等. 城郊農(nóng)田土壤Pb、Cd的空間變異與評(píng)價(jià)研究——以長沙市為例[J]. 長江流域資源與環(huán)境，2012，21（2）：195-203.

[8]白曉宇，袁峰，李湘凌，等. 銅陵礦區(qū)土壤重金屬元素的空間變異及污染分析[J]. 地學(xué)前緣，2008，15（5）：256-263.

[9]黃勇，郭慶榮，任海，等. 地統(tǒng)計(jì)學(xué)在土壤重金屬研究中的應(yīng)用及展望[J]. 生態(tài)環(huán)境，2004，13（4）：681-684.

[10]吳文勇，尹世洋，劉洪祿，等. 污灌區(qū)土壤重金屬空間結(jié)構(gòu)與分布特征[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2013，29（4）：165-173.

[11]姚榮江，楊勁松，楊奇勇. 典型農(nóng)區(qū)耕作土壤重金屬空間變異的穩(wěn)健性分析[J]. 環(huán)境科學(xué)與技術(shù)，2012，35（S2）：1-7.

[12]Wang X J，Zhang Z P. A comparison of conditional simulation，kriging and trend surface analysis for soil heavy metal pollution pattern analysis[J]. Journal of Environmental Science and Health，Part A：Toxic/Hazardous Substances and Environmental Engineering，1999，34（1）：73-89.

[13]孫波，周生路，趙其國. 基于空間變異分析的土壤重金屬復(fù)合污染研究[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2003，22（2）：248-251.

[14]蔡德所，王魁，黃景新，等. 喀斯特峰叢洼地原生林區(qū)土壤礦質(zhì)元素空間異質(zhì)性研究[J]. 中國水土保持，2010（5）：33-36.

[15]葛劍平，郭海燕，仲莉娜. 地統(tǒng)計(jì)學(xué)在生態(tài)學(xué)中的應(yīng)用（Ⅰ）——基本理論和方法[J]. 東北林業(yè)大學(xué)學(xué)報(bào)，1995，23（2）：88-94.

[16]Zhao Y C，Wang Z G，Sun W X，et al. Spatial interrelations and multi-scale sources of soil heavy metal variability in a typical urban-rural transition area in Yangtze River Delta region of China[J]. Geoderma，2010，156（3/4）：216-227.endprint

[1]史貴濤，陳振樓，李海雯，等. 城市土壤重金屬污染研究現(xiàn)狀與趨勢(shì)[J]. 環(huán)境監(jiān)測(cè)管理與技術(shù)，2006，18（6）：9-12，24.

[2]于小俸，唐嶷，胡玉. 酸性土壤重金屬含量的影響因素——以湖南壺瓶山自然保護(hù)區(qū)為例[J]. 廣東農(nóng)業(yè)科學(xué)，2011，38（9）：54-55.

[3]Lu J S，Liu Y，Zhang Z L，et al. Factorial kriging and stepwise regression approach to identify environmental factors influencing spatial multi-scale variability of heavy metals in soils[J]. Journal of Hazardous Materials，2013，261：387-397.

[4]毛竹，張世熔，李婷，等. 鉛鋅礦區(qū)土壤重金屬空間變異及其污染風(fēng)險(xiǎn)評(píng)價(jià)——以四川漢源富泉鉛鋅礦山為例[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2007，26（2）：617-621.

[5]何佳芳，何騰兵. 某鉛鋅冶煉廠廢棄地復(fù)墾整理區(qū)土壤重金屬污染評(píng)價(jià)[J]. 水土保持學(xué)報(bào)，2006，20（2）：97-101.

[6]張博，趙耕毛，劉兆普，等. 江蘇灘涂圍墾區(qū)土壤養(yǎng)分空間變異研究[J]. 江蘇農(nóng)業(yè)科學(xué)，2010（5）：461-464.

[7]劉瓊峰，李明德，吳海勇，等. 城郊農(nóng)田土壤Pb、Cd的空間變異與評(píng)價(jià)研究——以長沙市為例[J]. 長江流域資源與環(huán)境，2012，21（2）：195-203.

[8]白曉宇，袁峰，李湘凌，等. 銅陵礦區(qū)土壤重金屬元素的空間變異及污染分析[J]. 地學(xué)前緣，2008，15（5）：256-263.

[9]黃勇，郭慶榮，任海，等. 地統(tǒng)計(jì)學(xué)在土壤重金屬研究中的應(yīng)用及展望[J]. 生態(tài)環(huán)境，2004，13（4）：681-684.

[10]吳文勇，尹世洋，劉洪祿，等. 污灌區(qū)土壤重金屬空間結(jié)構(gòu)與分布特征[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2013，29（4）：165-173.

[11]姚榮江，楊勁松，楊奇勇. 典型農(nóng)區(qū)耕作土壤重金屬空間變異的穩(wěn)健性分析[J]. 環(huán)境科學(xué)與技術(shù)，2012，35（S2）：1-7.

[12]Wang X J，Zhang Z P. A comparison of conditional simulation，kriging and trend surface analysis for soil heavy metal pollution pattern analysis[J]. Journal of Environmental Science and Health，Part A：Toxic/Hazardous Substances and Environmental Engineering，1999，34（1）：73-89.

[13]孫波，周生路，趙其國. 基于空間變異分析的土壤重金屬復(fù)合污染研究[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2003，22（2）：248-251.

[14]蔡德所，王魁，黃景新，等. 喀斯特峰叢洼地原生林區(qū)土壤礦質(zhì)元素空間異質(zhì)性研究[J]. 中國水土保持，2010（5）：33-36.

[15]葛劍平，郭海燕，仲莉娜. 地統(tǒng)計(jì)學(xué)在生態(tài)學(xué)中的應(yīng)用（Ⅰ）——基本理論和方法[J]. 東北林業(yè)大學(xué)學(xué)報(bào)，1995，23（2）：88-94.

[16]Zhao Y C，Wang Z G，Sun W X，et al. Spatial interrelations and multi-scale sources of soil heavy metal variability in a typical urban-rural transition area in Yangtze River Delta region of China[J]. Geoderma，2010，156（3/4）：216-227.endprint