基于ArcGIS的污染場地土壤鉻空間格局分析

方敦 吳稱意 韋曉宇 龍小蘭

摘要 以某鉻鹽生產廢棄場地為例，運用空間數(shù)據(jù)管理平臺ArcGIS以普通克里格插值法、趨勢分析法評估污染物鉻（Cr）在場地土壤中的空間分布格局并實現(xiàn)其可視化表達。普通克里格插值分析表明場地污染物Cr主要集中在中上層土壤且Cr污染在各層中部污染較其他區(qū)域嚴重，Cr在場地中的空間分布與原廠工藝布局關聯(lián)度較高；趨勢分析與克里格插值空間分析結果具有較好的一致性，場地空間數(shù)據(jù)分層統(tǒng)計，合理地利用了空間插值分析中區(qū)域變量的自相關假設，提高了評估成果的可靠性。

關鍵詞 土壤；污染場地；空間分析；克里格插值； ArcGIS

中圖分類號 S181.3；X53；X825 文獻標識碼 A 文章編號 0517-6611（2014）16-05181-05

當前，城市化是中國經(jīng)濟社會發(fā)展的重要戰(zhàn)略部署，城市化過程中由于產業(yè)結構和土地利用空間調整已不可避免地涉及到騰退工業(yè)土地再利用問題[1]。土壤重金屬污染的隱蔽性和空間分異性[2]使得工業(yè)場地再利用具有潛在的土壤環(huán)境健康風險，因此，如何有效利用科學可靠的技術手段查清污染物在場地土壤空間的污染程度和立體空間分布格局是確保污染場地再利用環(huán)境安全性和經(jīng)濟合理性的關鍵。提高騰退工業(yè)場地再開發(fā)過程的環(huán)境安全性以準確評價場地污染物的空間分布格局為前提，然而，污染物與場地土壤環(huán)境的相互作用隨區(qū)域分布的不確定性給準確認識污染物的空間分布帶來了巨大困難。

地理信息系統(tǒng)（Geographic Information System，GIS）是環(huán)境污染物空間分析的有效工具[3-5]，地統(tǒng)計模塊是連接GIS與地統(tǒng)計學的橋梁。地統(tǒng)計學是以區(qū)域變化量為基礎，借助變異函數(shù)，研究自然現(xiàn)象的空間變異與空間結構的理論與方法，是研究空間不確定性的有效方法[2，6-7]?；诖耍撗芯恳糟t鹽污染場地為例，運用地統(tǒng)計學普通克里格（Ordinary Kriging）插值法和趨勢分析法以分層統(tǒng)計方式開展基于ArcGIS的污染物鉻的空間格局分析研究，以期在污染物空間格局分析的可靠性和直觀性方面為場地環(huán)境評估及管理提供思路和技術方法參考。

1 材料與方法

1.1 場地概況 研究場地為北方某鉻鹽生產企業(yè)退役場地，位于37°42′～38°13′ N，113°～114°18′ E，下覆土壤pH在6.5～10.9之間。原企業(yè)歷史上均采用鈣焙燒法生產紅礬鈉，生產車間主要有原料庫房、混料車間、煅燒車間、浸取車間、成品車間、硫酸儲存區(qū)、鉻渣堆場及鍋爐房等。研究區(qū)域地勢由西南向東北傾斜，出露地層為古生界寒武系張夏組一段（∈Z1）、二段（∈Z2）巖層及新生界第四系（Q）沉積物，第四系亞粘土層厚度3.0～25.9 m不等，其下為致密巖石層，節(jié)理裂隙不發(fā)育，完整性好，地下水埋深300 m以上。研究場地土壤基本理化性質如表1所示。

1.2 樣品采集 布點以原廠歷史沿革、工藝布局和地層條件為參考，采樣布點分階段采用專業(yè)判斷和系統(tǒng)布點相結合的方法，第一階段在具有最大污染可能性的車間中心區(qū)域設置采樣點，第二階段對整個研究區(qū)域劃分30 m×30 m網(wǎng)格進行系統(tǒng)網(wǎng)格布點并根據(jù)第一階段樣品測試結果對重點區(qū)域適當加密，場地共設置水平土壤采樣點27個，采樣點編號為S1～S27。同時，根據(jù)水文地質條件和現(xiàn)場X射線熒光光譜（Niton XL2 XRF）檢測結果設置垂向采樣深度為5.0 m的規(guī)定深度采集剖面樣且各階段采樣點數(shù)據(jù)均按上層（0～1.0 m）、中層（1.0～3.0 m）和下層（3.0～5.0 m）內采樣點的算術平均值分層統(tǒng)計，各樣點垂向采樣遵循“等量”原則。采樣過程中每個定深度采樣后對鉆頭徹底清洗，樣品于帶聚四氟乙烯密封墊的玻璃瓶中4 ℃下保存。

1.3 分析測試方法 土壤樣品經(jīng)風干、除雜、磨細、過篩（10目、60目、100目）等前處理后備用，但用于六價鉻（Cr6+）測定的樣品除剔除砂礫石和植物根系外不作處理。研究場地土壤樣品分析測試指標按常規(guī)指標和特征污染指標兩類，依據(jù)原企業(yè)歷史沿革深入開展污染識別結果進行設計。該研究僅以其特征污染物鉻（Cr）和六價鉻（Cr6+）為研究對象進行探討，Cr和Cr6+分別按JY/T 015.1996[8]、USEPA 3060A[9]及USEPA 7196A[10]操作規(guī)程采用電感耦合等離子體原子發(fā)射光譜（ICP.AES）和堿消解-紫外/可見分光光度法分析測定。分析測試過程質量控制手段包括隨機抽樣平行測試、全程序空白測試及樣品加標回收率測定。

1.4 空間分析

1.4.1 樣本正態(tài)性檢驗。樣本正態(tài)性是地統(tǒng)計插值法的基本假設之一，正態(tài)QQPlot圖是檢驗樣本數(shù)據(jù)是否符合正態(tài)分布的有效方法[11]，在ArcGIS地統(tǒng)計分析模塊（Geostatistical Analyst）中易于實現(xiàn)，通過模塊下的探索性數(shù)據(jù)分析（Explore Data）即可實現(xiàn)。

1.4.2 Kriging插值空間分析?？臻g分析是GIS的重要功能，可提取和傳輸空間信息。污染場地在理論上可認為是具有漸變特征的連續(xù)空間。由于原企業(yè)工藝布局和污染物在土壤中的遷移轉化等因素對場地污染物含量分布的影響，使得場地污染物含量的空間分布同時具有隨機性和內部相似性，即場地污染物含量的空間分布具有區(qū)域變化量的特征。

隨機性因素可由系統(tǒng)性采樣點對區(qū)域總體特征的代表性而最大限度地在樣本數(shù)據(jù)中得以體現(xiàn)，因此，空間分析主要應解決污染場地區(qū)域變量數(shù)據(jù)的內部相似性問題。地統(tǒng)計學Kriging插值是一種以半方差函數(shù)理論和結構分析為基礎的對區(qū)域變化量進行線性、無偏、最優(yōu)估計的一種局部空間插值法[12]。Kriging插值在插值分析中不僅考慮未知樣點與已知樣點的距離遠近而且考慮了兩者的空間方位關系，其未知樣點與已知樣點的關系如（1）式[12]所示?；谝陨戏治觯瑢龅匚廴疚锏你t含量屬性數(shù)據(jù)可采用地統(tǒng)計學空間內插法進行空間分析并采用連續(xù)的空間漸變模型來實現(xiàn)，該研究應用ArcGIS 9.3空間分析模塊中的Geostatistical Analyst工具對研究場地污染物進行地統(tǒng)計學Ordinary Kriging空間插值分析。2 結果與分析

2.1 樣本數(shù)據(jù)正態(tài)性 Ordinary Kriging插值法屬于地統(tǒng)計學插值方法，樣本數(shù)據(jù)正態(tài)分布是地統(tǒng)計學插值分析的基本假設，因此，用于空間插值的樣本數(shù)據(jù)需符合正態(tài)分布假設[11]，若不符合應盡可能采用對數(shù)（Log）、冪變換（BOX.COX）等可逆變換使其符合正態(tài)分布的基本假設[13]。在正態(tài)QQ圖中，若數(shù)據(jù)的分布近似成一條直線，則數(shù)據(jù)服從正態(tài)分布[14]。由圖1～3可知，經(jīng)Log變換后的研究場地各層數(shù)據(jù)均近似呈直線分布，表明研究場地調查數(shù)據(jù)經(jīng)Log變換后符合Kriging插值法的基本假設；圖2顯示在低值區(qū)和高值區(qū)出現(xiàn)了個別異常點，S4、S5、S7及S25號樣點位于場地鉻渣堆存區(qū)（防滲），其檢測值相對樣本整體而言具有跳躍性，結合場地實際，研究判定其為局部真實離群值[14]，與場地原工藝布局等因素有關，需要指出的是，在插值分析過程中，對S4、S5、S7及S25號具有跳躍性的局部真實離群值樣點數(shù)據(jù)采用隔離-疊加措施處理使其盡可能符合真實空間分布且不影響整體的結構性，提高空間分析的準確性，即在GIS中將修剪（Clip）掉異常點后的插值結果與原插值柵格平面經(jīng)異常點邊界面數(shù)據(jù)提?。‥xtract）的結果疊加。

2.2 鉻污染空間分布 場地污染物Cr的Ordinary Kriging插值空間分布如圖4～6所示。從空間分布圖可知：研究場地污染物Cr集中分布在中上土層中，上層Cr污染呈現(xiàn)出污染嚴重、污染覆蓋范圍大的特點，垂向具有由上至下污染范圍逐漸縮小的整體空間分布特征；場地整體污染Cr空間分布情況與原廠工藝布局關聯(lián)度較高，中上層污染程度大的區(qū)域與場地原工藝布局基本一致；場地主要污染深度范圍在0～3.0 m之間，下層土壤Cr含量未超出標準值，場地污染土壤修復范圍分層疊加空間分布（圖7）更直觀地體現(xiàn)了Cr在場地中的這種層間空間分布特征。下層Cr空間分布（圖6）與中上層（圖4、圖5）相比，插值結果具有顯著差異，表明下層土壤Cr含量在土壤污染的意義上來說在垂向上已脫離了研究場地垂向污染范圍自相關的結構性特征。

場地地勢由東南向西北傾斜，在降雨沖刷作用下將Cr運移至場地西北面致使整個場地第一層在長期沖刷和淋溶作用下被嚴重污染；圖5、圖6顯示，場地中、下層污染范圍和程度有較大幅度的縮減，下層土壤Cr含量狀況良好，這與場地土壤性質、pH及Cr的特性有關[15-17]，是土壤pH、土壤膠體電負性對六價鉻負離子的向下遷移推動作用、降雨淋溶作用、土壤有機質對三價鉻和六價鉻的吸附、還原吸附共同作用的結果，同時六價鉻在土壤介質中通過形成CaCrO4、BaCrO4、PbCrO4等難溶鹽[18]也是阻礙場地Cr垂向遷移的原因之一[16]。一方面，研究場地土壤pH呈較強堿性為土壤中Cr形成沉淀和殘渣態(tài)提供了有利條件[19]，各層次土壤中的Cr均具有被沉淀和吸附進而遷移能力降低的趨勢；另一方面，Cr6+是研究場地主要污染物之一，Cr6+在環(huán)境中除與金屬形成難溶鹽外主要以Cr2O2-7、HCrO-4和CrO2-4等陰離子形式存在[18]，由于土壤膠體電負性而使其具有向下垂向遷移的能力，但同時，土壤有機質對Cr6+的還原作用在促進其垂向遷移方面又與土壤膠體形成拮抗作用。

2.3 趨勢分析 全局趨勢分析認為研究表面由全局趨勢和隨機短程變異兩部分組成，可通過空間抽樣數(shù)據(jù)擬合一個數(shù)學曲面來反映空間分布的變化情況，在ArcGIS Geostatistical Analyst中通過Trend Analysis工具即可實現(xiàn)。趨勢分析工具將樣點值投影到一個由東西向和南北向平面構成的正交平面上后通過多項式在投影平面上擬合曲線形成趨勢分析圖[14]，若擬合曲線平直則表示該數(shù)據(jù)內部不存在明顯趨勢[11，20]。由趨勢分析圖（圖8～圖10）可知：研究場地上、中層在東南-西北向和西南-東北向上污染物Cr含量分布呈現(xiàn)“U”形趨勢，即呈中部高四周相對較低的整體趨勢，趨勢分析結果與研究場地主要污染工藝在中部集中的原廠工藝布局具有一致性，與Kriging插值空間分析結果相符。而趨勢分析圖（圖10）呈現(xiàn)平滑過渡的直線，表明下層土壤Cr含量在場地空間內不具有明顯趨勢性，結合Kriging插值分析結果可判定其已脫離了研究場地Cr污染物分布范圍的區(qū)域變量范疇。

3 結論

當前，環(huán)境問題已成為阻礙我國可持續(xù)發(fā)展的重要因子，我國也越來越重視污染土壤的治理修復問題，然而，要使治理修復工作有條不紊行之有效地開展，準確而直觀的場地評估結論具有不可替代的作用。該研究在空間數(shù)據(jù)分層統(tǒng)計和樣本數(shù)據(jù)檢驗的基礎上通過Kriging插值分析和趨勢分析發(fā)現(xiàn)研究場地土壤鉻污染主要集中在第一、二土層中，第三層土壤幾乎不受污染，從插值分析和趨勢分析來看均表現(xiàn)為第三層與前兩層土壤的屬性數(shù)據(jù)脫離了自相關的結構性特征，場地土壤污染調查中分層數(shù)據(jù)統(tǒng)計方式可確保研究場地污染物空間分布分析的可靠性，是污染場地評價空間分析有益的實例借鑒。

參考文獻

[1] 謝劍，李發(fā)生.中國污染場地的修復與再開發(fā)的現(xiàn)狀分析[R].華盛頓： 世界銀行，2010：1-4.

[2] 朱文轉，夏北成，江學頂.惠州市土壤重金屬污染景觀空間格局特征研究[J].中國環(huán)境科學，2009，29（8）：885-891.

[3] HOEK G，F(xiàn)ICHER P，BRANDT P V D，et al.Estimation of long-term average exposure to outdoor air pollution for a cohort study on mortality[J].Journal of Exposure Analysis and Environmental Epidemiology，2001，11（6）：459-469.

[4] BEELEN R，HOEK G，F(xiàn)ISCHER P，et al.Estimated long-term outdoor air pollution concentrations in a cohort study[J].Atmospheric Environment，2007，41：1343-1358.

[5] BANJA M，COMO E，MURTAJ B，et al.Mapping air pollution in urban Tirana area using GIS[C]//International Conference SDI.Skopje，Macedonia，2010.

[6] GRAMATICA P，BATTAINI F，GIANI E，et al.Analysis of mosses and soils for quantifying heavy metal concentrations in sicily： A multivariate and spatial analytical approach[J].Environ Sci & Pollut Res，2006，13 （1）：28-36.

[7] 謝正苗，李靜，徐建明，等.杭嘉湖平原土壤中氟元素的空間分布特征[J].中國環(huán)境科學，2005，25（6）：719-723.

[8] 國家教育委員會.JY/T 015-1996 感耦等離子體原子發(fā)射光譜方法通則[S].1996.

[9] US EPA Method 3060A.Alkaline digestion for hexavalent chromium[S].1996.

[10] US EPA Method 7196A.Chromium，hexavalent（colorimetric）[S].1992.

[11] DASHTPAGERDI M M，SADATINEJAD S J，BIDAKI R Z，et al.Evaluation of air pollution trend using GIS and RS applications in south west of Iran [J].Journal of the Indian Society of Remote Sensing，2014，42（1）：179-186.

[12] MARKO K，AL.AMRI N S，ELFEKI A M M.Geostatistical analysis using GIS for mapping groundwater quality： case study in the recharge area of Wadi Usfan，western Saudi Arabia[J].Arabian Journal of Geosciences，2013，13：1156-1157.

[13] WANG J，CHEN Z L，SUN X J，et al.Quantitative spatial characteristics and environmental risk of toxic heavy metals in urban dusts of Shanghai，China [J].Environ Earth Sci，2009，59：645-654.

[14] JOHNSTON K，HOEF J M V，KRIVORUCHKO K，et al.Using ArcGIS Geostatistical Analyst[M].New York Street： ESRI Inc，2001.

[15] HON.WAH L M，YUK.WAI T A，CHUN.CHIU C，et al.Speciation study of chromium，copper and nickel in coastal estuarine sediments polluted by domestic and industrial effluents[J].Marine Pollution Bulletin，1997，34（11）：949-959.

[16] 榮偉英，周啟星.鉻渣堆放場地土壤的污染過程、影響因素及植物修復[J].生態(tài)學雜志，2010，29（3）：598-604.

[17] 陳英旭，何增耀，吳建平.土壤中鉻的形態(tài)及其轉化[J].環(huán)境科學，1994，15（3）：53-56.

[18] 尹晉，馬小東，孫紅文，等.電動修復不同形態(tài)鉻污染土壤的研究[J].環(huán)境工程學報，2008，2（5）：684-689.

[19] 桂新安，楊海真，王少平，等.鉻在土壤中的吸附解吸研究進展[J].土壤通報，2007，38（5）：1007-1012.

[20] PRABAHARAN S，MANONMANI R，RAMALINGAM M，et al.Geostatistical modeling for groundwater pollution in Salem，tamilnadu.A GIS based approach [J].International Journal of Engineering Science and Technology，2011，3（2）：1273-1278.