基于GIS的惠州市土壤多環芳烴空間分布特征研究

王珂，楊國義，王玉瑯，張天彬，萬洪富*

1. 廣東省生態環境技術研究所，廣東 廣州 510650；2. 珠江流域水土保持監測中心站，廣東 廣州 510611 3. 華南農業大學資源環境學院，廣東 廣州 510642

基于GIS的惠州市土壤多環芳烴空間分布特征研究

王珂1，3，楊國義1，王玉瑯2，張天彬1，萬洪富1*

1. 廣東省生態環境技術研究所，廣東 廣州 510650；2. 珠江流域水土保持監測中心站，廣東 廣州 510611 3. 華南農業大學資源環境學院，廣東 廣州 510642

為了解惠州市土壤中多環芳烴的空間分布特征，采集了該市不同區縣42個不同土地利用類型的表層土壤樣品，并結合地理信息系統（GIS）技術分析惠州市土壤多環芳烴的空間分布特征，探清惠州市土壤多環芳烴的殘留現狀、來源和擴散趨勢等，為惠州市土壤多環芳烴的污染治理奠定理論基礎。結果表明，土壤中多環芳烴總量介于35.40~534.5 μg·kg-1，平均值為123.09 μg·kg-1?？臻g分布分析結果顯示，Chr和Bkf具有強烈的空間相關性；Nap、Ace、Acy、Flu、Phe、Ant、Fla、Baa、Daa和Bbf具有中等空間相關性；Pyr、Bap、I1p和Bgp空間相關性很弱。趨勢分析顯示，16種多環芳烴在東西和南北方向上呈現出明顯的增減規律性。Nap和Bkf由東向西、由南向北逐漸升高，而Flu、Pyr、Baa、Bbf、Ipy、Daa則恰好相反；Acy、Ace、Ant均由東向西、由南向北急劇升高，與Phe、Fla的變化趨勢恰好相反；Chr和Bap在東西、南北方向沒有明顯趨勢變化；Bgp在東西方向上，由東向西逐漸下降，而在南北方向上，呈明顯的倒U形。利用GIS軟件地統計模塊進行了普通克里格（Ordinary Kriging）插值，結果顯示16種多環芳烴組分分布特征各不相同，Ace、Acy和Ant等3種多環芳烴組分含量呈現出自西北向東南方向遞增的趨勢；Phe、Fla、Pyr、Bkf、Chr、I1p、Daa和Bgp等則恰好相反；Nap、Bbf和Bap 等沒有明顯的分布趨勢。Nap最高含量區面積很小，Bap高含量區面積較大，Baa在整個惠州市含量均較低，但在惠城區出現明顯的島狀高含量分布區。以上結果反映了不同區域多環芳烴來源的不同，而多環芳烴的來源與當地經濟社會發展情況及產業結構關系密切，主要包括工業“三廢”、交通廢氣、石油燃燒和、生活垃圾焚燒等。

土壤；多環芳烴；空間分布；GIS；惠州市

多環芳烴（PAHs）是廣泛分布于環境中的一類典型持久性有機污染物，因其具有強烈的致癌、致畸、致突變作用及光致毒效應（孫紅文等，1998），而受到國內外環境科學界的普遍關注（鄭太輝等，2014）。土壤作為一種極其重要的環境介質，是包括多環芳烴在內的持久性有機污染物的遷移、轉化、累積中轉站，因此，土壤多環芳烴污染引起了世界各國的廣泛關注（Trapido，1999；Mielke et al.，2001）。以往針對土壤多環芳烴的研究主要集中在多環芳烴含量及污染評價方面（Honda et al.，2007；Fu et al.，2003；Maliszewska-Kordybach，1996）。國外借助GIS技術進行土壤多環芳烴空間分布的研究已有很多成果（Carlon et al，2001；Colin et al.，1996），國內鄭一等（2003）、段永紅等（2005）等曾運用GIS技術對天津地區表層土壤多環芳烴進行了空間分析，結果表明，天津地區土壤多環芳烴普遍受燃煤釋放、燃油和其他工業活動排放等的影響；彭馳（2009）對北京市土壤中多環芳烴進行了空間分析與生態評價，通過比值法確定北京土壤中多環芳烴主要來源于煤燃燒和化石燃料燃燒，等等。

珠江三角洲地區乃至廣東省范圍內，有關土壤污染物的研究以重金屬為主，其次為鄰苯二甲酸酯和有機氯農藥等，而有關土壤中多環芳烴的研究報道卻不是很多，已有的研究也主要圍繞湛江、東莞、深圳等城市展開。例如，馬瑾等（2012）2624基于多元地統計和GIS對湛江市土壤多環芳烴的污染狀況和空間分布進行分析，結果表明，湛江市土壤存在不同程度的多環芳烴污染，主要集中在經濟發達地區，受人類活動影響大，水田更容易吸附和累積多環芳烴。黃超勝（2012）對貴嶼及周邊地區進行研究，結果表明，貴嶼及周邊農業土壤多環芳烴污染是多種污染源綜合作用的結果，其中燃燒源（包括電子垃圾燃燒）是其主要來源。張天彬等（2005）、馬瑾等（2011）曾對東莞市土壤多環芳烴進行了研究，前者結果表明，東莞市土壤中多環芳烴主要來自化石燃料高溫燃燒產物的大氣沉降，來源于石油的多環芳烴輸入只占小部分；而后者來源則主要是燃燒源、石油源和煤燃燒源，說明隨著時間與產業的變化，多環芳烴的污染也會發生很大的變化。張天彬等（2008）對深圳表層土壤中多環芳烴的研究結果顯示，深圳土壤中大部分多環芳烴的檢出率在80%以上，致癌性多環芳烴占總量的51.9%，主要來源于燃燒來源，部分為石油來源。陳來國等（2004）對廣州周邊菜地中多環芳烴污染現狀進行了分析，結果表明廣州周邊菜地多環芳烴屬中度污染。然而，目前對于惠州市土壤多環芳烴的污染研究卻未見任何報道。

土壤多環芳烴空間分布特征研究通常采用土壤要素空間插值中最常用的克里格插值法（高俊峰等，2002）。該方法是借助已有樣點的分布、密度、大小及其與待估樣點相互間的空間分布位置等幾何特征和空間結構，賦于已有樣點一定的系數，再進行加權平均以確定未采樣區域的值，由此形成空間分布圖?？臻g趨勢能夠很好地反映空間物體在空間區域上變化的主體特征，忽略局部的變異特征，揭示空間物體的總體規律。趨勢面分析是根據空間抽樣數據，擬合成一個數學曲面，從而反映空間分布的變化情況。趨勢分析圖中的每一條豎線均代表一個數據點的值和位置，這些點被投影到一個東西向和一個南北向的正交平面上，通過投影點可以得到一條最佳擬合線，借助這條擬合線可以得知空間物體在特定方向上的變化趨勢。如果該線是直的，表示空間物體不存在變化趨勢。

塊金值（C0）與基臺值（C0+C）之比可以反映系統變量的空間相關性程度。C0/(C0+C)<25%，說明系統變量具有強烈的空間相關性；25%75%，說明系統變量空間相關性很弱（王政權，1999）。空間相關性強，說明該污染物可能存在區域性污染，即污染的面積比較大?？臻g相關性弱，則說明該污染物可能僅存在局部性污染，或者存在較明顯的點狀污染源。

選取位于珠江三角洲的惠州市作為研究區域，借助GIS技術和克里格插值法等分析該市土壤多環芳烴含量及空間分布特征，并對影響土壤多環芳烴空間分布的因素進行探討，旨在確定惠州市表層土壤中多環芳烴的含量及殘留現狀，以期為惠州市土壤多環芳烴的治理提供理論基礎。

1 研究區概況

惠州市地處北緯22°24′~23°57′、東經113°51′~ 115°28′之間，位于廣東省東南部，珠江三角洲東北端，南臨南海大亞灣，毗鄰廣州、東莞、深圳、香港，設有大亞灣經濟技術開發區和仲愷高新技術產業開發區兩個國家級開發區，屬珠三角經濟區。該市屬亞熱帶海洋性季風氣候區，冬半年（非汛期，10—次年3月）主要受大陸季風的影響，盛行東北風，而夏半年（汛期，4—9月）主要受西南和東南季風的影響，盛行東南風，臺風較多。氣候溫和，年平均氣溫在21 ℃左右。南北多丘陵，中部多平原和臺地。東江、西枝江橫貫境內。

2015年，惠州市實現地區生產總值（GDP）3140.03億元，增長9.0%。人均GDP為66231元，按平均匯率折算為10634美元。惠州是一座快速崛起的現代工業新城。在中海殼牌100萬噸乙烯、中海油1200萬噸煉油、TCL液晶模組、德賽、華陽集團、科銳LED等龍頭項目帶動下，石化與電子兩大支柱產業成為跨越發展、轉型升級的強大引擎，形成了“2+4”的工業支柱產業格局，即數碼、石化兩大支柱和服裝、制鞋、水泥和汽車及零部件的4個比較具有區域優勢的產業，已成為首批國家級電子信息產業基地、國家級火炬計劃“數碼視聽產業基地”和國家重要石油化工基地。隨著工業的快速發展，該市不少區域出現了嚴重的環境污染問題。

2 材料與方法

2.1 布點原則

根據惠州市工農業布局、“三廢”排放狀況、灌溉水類型、土壤類型確定采樣點，考慮技術力量與財力條件，采樣點重點布設在典型之地塊。

2.2 采樣方法

采取0~20 cm表層土壤，采用“X”法隨機采取5個以上采樣點，充分混合成1個代表樣，采用四分法留取1 kg土壤樣品，共采集表層土壤樣品42個。土壤樣品置于布袋中，并置于冰箱中低溫保存。每個樣品在采集過程中，都利用全球定位系統（GPS）進行定位，采樣點分布見圖1。

2.3 分析方法

圖1 惠州市土壤采樣點分布圖Fig. 1 Location of sampling sites in Huizhou City

多環芳烴的測定參考馬瑾等的文獻（2012）。分析前，將土壤樣品從冰箱中取出，于低溫下風干，用瑪瑙研缽研磨后過60目篩。準確稱取土樣20 g于索式抽濾筒中，在250 mL平底燒瓶中加入200 mL二氯甲烷、2 g經活化的銅片以及回收率指示物（萘-d8、苊-d10、?-d12和苝-d12），在46 ℃恒溫水浴鍋上進行索式抽提48 h，冷卻循環水溫調節為22 ℃，回流速度控制在5~6 time·h-1。將提取液置于旋轉蒸發儀上濃縮至1 mL，加入10 mL正己烷轉換劑，再濃縮至1~2 mL。濃縮樣過層析柱進行分離純化，凈化柱為進口層析柱（內徑1 cm），采用正己烷濕法裝柱，由下至上依次為硅膠12 cm，氧化鋁6 cm，無水硫酸鈉2 cm。先用15 mL正己烷淋洗出烷烴類化合物，再用70 mL二氯甲烷-正己烷（3∶7，V/V）淋洗出多環芳烴。淋洗液經正己烷轉換溶劑后，用高純氮氣吹至0.2 mL，加入內標物進行定量分析。采用HP6890Ⅱ GC/HP5973 MSD進行分析測定，采用惠普化學工作站（HP Chemstation）進行數據處理。主要分析了被美國環保總署（USEPA）列為優先控制的16種多環芳烴（PAHs）的含量（Berset et al.，1998），分別為萘（Nap， 2環）、苊（Acy，3環）、二氫苊（Ace，3環）、芴（Flu，3環）、菲（Phe，3環）、蒽（Ant，3環）、熒蒽（Fla，4環）、芘（Pyr，4環）、苯并（a）蒽（Baa，4環）、?（Chr，4環）、苯并（b）熒蒽（Bbf，5環）、苯并（k）熒蒽（Bkf，5環）、苯并（a）芘（Bap，5環）、茚并（1, 2.3-cd）芘（Ipy，6環）、二苯并（a, h）蒽（Daa，5環）、苯并（ghi）苝（Bgp，6環）。多環芳烴按其物化性質可以分為兩類（彭馳，2009）：一類為帶有2~3個苯環的低分子量芳烴，這些化合物易揮發，對水生生物有一定的毒性；第二類為帶有4~6個苯環的高分子量芳烴，不易揮發，具有致癌、致突變作用。

按照USEPA要求，在每個樣品中都加入了回收率指示物（4種回收率指示物萘-d8、苊-d10、?-d12和苝-d12的平均回收率分別為50.86%、69.89%、88.95%和83.68%），對整個分析流程的回收率進行監控。同時每15個樣品做1個空白對照，加標空白，基質加標，以保證分析過程中的準確性。

2.4 數據處理

數據描述性統計分析在SPSS 13.0軟件中完成。應用ArcGIS 9.0軟件地統計模塊擬合多環芳烴各組分的理論模型（擬合主要采用有塊金效應的球狀模型），并根據理論模型提供的參數，對未采樣點的區域化變量的取值進行線性無偏最優估計，生成惠州市土壤16種優控多環芳烴的克里格插值空間分布圖。

3 結果與討論

3.1 土壤多環芳烴的統計性描述

惠州市土壤中多環芳烴基本統計結果見表1。

由表1可知，惠州市土壤中16種優控多環芳烴均有不同程度的檢出，其含量也存在很大的差異，其中含量最高的組分Bbf均值達到22.87 μg·kg-1，而含量最低的Ace均值僅為0.73 μg·kg-1。多環芳烴總量介于35.40~534.5 μg·kg-1之間，平均值為123.09 μg·kg-1，大大超過了Edwards（1983）提出的土壤中內源性PAHs的含量（1~10 μg·kg-1），即使是最低∑PAHs含量也接近了土壤內源性多環芳烴含量的最高值，因此惠州市土壤多環芳烴主要為外源性輸入。

表1 惠州市土壤多環芳烴基本統計結果Table 1 Statistical results of soil PAHs in Huizhou City μg·kg-1

根據荷蘭Maliszewska-Kordybach（1996）對土壤中多環芳烴污染的分類標準（<200 μg·kg-1為未污染，200~600 μg·kg-1為輕度污染，600~1000 為μg·kg-1中度污染，>1000 μg·kg-1為嚴重污染），惠州市表層土壤中多環芳烴污染屬于輕度污染?！锻寥牢廴痉乐涡袆佑媱潯罚ê喎Q“土十條”）中對農用地進行了劃分，規定未污染和輕微污染土壤劃分為優先保護類，輕度和中度污染的劃分為安全利用類，重度污染的劃為嚴格管控類。所以，針對惠州市表層土壤多環芳烴污染，應將其列為優先保護并安全利用。

變異系數可以反映變量的離散程度，也可以粗略地反映變異強度，本研究中Ace、Ant、Baa、Bbf、Bkf、Ipy、Bgp等7種多環芳烴組分存在中等變異（10%100%）。

3.2 土壤多環芳烴空間趨勢分析

借助ArcGIS軟件對惠州市土壤16種優控多環芳烴進行趨勢分析，分析結果見圖2。

圖2顯示16種優控多環芳烴的趨勢分布如下：Nap和Bkf由東向西、由南向北逐漸升高，而Flu、Pyr、Baa、Bbf、Ipy、Daa則恰好相反；Acy、Ace、Ant的趨勢均為由東向西、由南向北急劇升高，與Phe、Fla的趨勢恰好相反。Chr和Bap在東西、南北方向近似于一條直線，趨勢變化不明顯。Bgp的變化趨勢比較特殊，在東西方向上，由東向西逐漸下降，而在南北方向上呈明顯的倒U形。

3.3 土壤多環芳烴空間分布特征

圖2 惠州市土壤多環芳烴趨勢面分析透視面圖Fig. 2 Trend analysis of soil PAHs in Huizhou City

利用半變異函數，分別采用不同的模型進行擬合，得到模型的相關參數值，在對不同參數模型的比較可以參考預測誤差中的幾個指標，符合如下4個標準的模型是最優的：標準平均值最接近于0；均方根預測誤差最??；平均標準誤差最接近均方根預測誤差；標準均方根預測誤差最接近于1。擬合理論模型及相關參數見表2。

從表2可知，C0/(C0+C)<25%的有Chr和Bkf，即這兩種多環芳烴組分具有強烈的空間相關性，污染面積較大，可能存在區域性污染；Nap、Ace、Acy、Flu、Phe、Ant、Fla、Baa、Daa和Bbf等10種多環芳烴組分C0/(C0+C)介于25%~75%之間，具有中等空間相關性，表明空間異質性是隨機因素和空間結構共同作用的結果，即土壤性質（如土壤母質、成土作用、地形等）和人類活動（如工業發展、耕作制度、施肥等）共同作用的結果。Pyr、Bap、I1p和Bgp等4種多環芳烴組分C0/(C0+C)>75%，空間相關性很弱，可能存在局部性污染，或者存在較明顯的點狀污染。

表2 惠州市土壤多環芳烴變異函數理論模型及其相關參數Table 2 Theoretic model and parameters of semevariance of soil PAHs in Huizhou City

根據所得半方差理論模型和參數進行克里格插值，得出惠州市土壤多環芳烴的空間分布圖（圖3）。

根據克里格插值圖可將惠州市土壤16種優控多環芳烴空間分布特征歸為4類：

（1）Ace、Acy和Ant等3種多環芳烴組分含量呈現自西北向東南方向遞增的趨勢。

（2）Phe、Fla、Pyr、Bkf、Chr、I1p、Daa和Bgp等8種多環芳烴組分含量總體上呈現自西北向東南遞減的趨勢，其中，Chr和I1p分別在博羅縣和龍門縣出現明顯的島狀高含量分布區，Daa在龍門縣和博羅縣分別出現1個島狀高含量區。

（3）Nap、Bbf和Bap等3種多環芳烴組分沒有明顯的分布趨勢，呈現出零星的島狀分布的高含量區，Nap最高含量區很小，出現在惠陽與惠東縣交界區域；而Bap高含量區面積較大，主要分布在博羅縣境內。

（4）Baa在整個惠州市含量均較低，但在惠城區呈現出明顯的島狀高含量分布區，而且以惠城區為中心，隨著距離增大，Baa濃度逐漸降低；在惠東縣也出現1條狹長的高含量分布帶。

一般而言，低環的多環芳烴主要起源于自然界和沒有經過燃燒的石油類產品，而高環的多環芳烴主要來自于各類燃燒源（彭馳，2009）。Ace、Acy、Ant、Phe、Fla、Pyr、Bkf、Chr、I1p、Daa、Bgp等11種PAHs在惠州市所呈現出的截然不同的空間分布特征，反映了惠州市不同區域多環芳烴來源的不同。惠州市東南部惠東縣、惠陽縣和惠城區土壤多環芳烴以低環（2環、3環）為主，目前，由于地理優勢，惠州市將石化產業列為了重點發展產業，石化基地位于大亞灣石化區，該石化區有兩大龍頭項目，即中國海油惠州煉油項目和殼牌南海石化項目，McCready et al.（2000）研究認為，低環PAHs主要來自石油源，因此，大亞灣石化產業與惠州市東南部土壤低環多環芳烴污染存在密切關系。

西北部的博羅縣和龍門縣土壤多環芳烴則以高環為主，高環多環芳烴中的Phe、Fla和Pyr是煤燃燒的主要產物（Harrison et al.，1996），博羅、龍門兩縣存在較多的熱電廠、冶金廠等，這些企業的生產均以煤為燃料，故其是該地區高環多環芳烴的主要來源。另外，Bkf是柴油動力燃燒的指示物，這可能是龍門縣和博羅縣土壤中Bkf含量普遍偏高的原因之一。

惠城區是惠州市的市中心區域，經濟發達，汽車保有量多，人類經濟社會活動對環境影響深刻。根據McCready et al.（2000）的研究，燃燒源產生的多環芳烴以高環為主，Baa為4環多環芳烴，因此，惠城區土壤Baa含量偏高與人類經濟社會活動存在密切關系。另外，Simcik et al.（1999）的研究認為Baa是天然氣燃燒的標志物?；葜菔杏?002年率先在大亞灣建設了全國首批天然氣電廠，天然氣發電與過去的燃煤發電相比，所排放的廢棄物較少，更環保，效益也更高，但本次研究發現在大亞灣地區土壤中Baa的含量較高，結合Simcik的研究結果，推測天然氣燃燒對大亞灣地區土壤Baa含量偏高有一定貢獻。

圖3 惠州市土壤多環芳烴克里格插值圖Fig. 3 Kriging interpolation map of soil PAHs in Huizhou City

Bap是一種典型的熱解成因多環芳烴，其來源主要包括工業生產、垃圾焚燒及失火、機動車尾氣排放等，各排放源的貢獻比例分別為62.3%、26.8%、10.9%（岳敏等，2003；Thomas et al.，2004），可見工業生產對Bap排放的貢獻最大，這可能也是博羅縣土壤Bap普遍含量偏高的一個重要原因。

另外，對多環芳烴的克里格插值分析很好地印證了空間趨勢分析的結果。

4 結論

（1）惠州市土壤中16種優控多環芳烴均有不同程度的檢出，其含量也存在很大的差異，Ace、Ant、Baa、Bbf、Bkf、Ipy、Bgp等7種多環芳烴組分存在中等變異（10%100%），說明惠州市多環芳烴污染受人類活動影響很大。

（2）惠州市東南部惠東縣、惠陽縣和惠城區土壤多環芳烴以低環（2環、3環）為主，其污染源主要來自大亞灣石化產業、天然氣燃燒、汽車尾氣等；西北部的博羅縣和龍門縣土壤多環芳烴則主要以高環為主，其污染源主要來自工業生產、機動車尾氣排放、熱電廠、冶金廠等。

（3）16種多環芳烴的空間分布規律表現為：Nap和Bkf由東向西、由南向北逐漸升高，而Flu、Pyr、Baa、Bbf、I1p、Daa則恰好相反；Acy、Ace、Ant由東向西、由南向北急劇升高，與Phe、Fla的變化趨勢恰好相反；Chr和Bap在東西、南北方向沒有明顯趨勢變化；Bgp變化趨勢比較特殊，在東西方向上，由東向西逐漸下降，而在南北方向上，呈明顯的倒U形。

（4）惠州市土壤中16種多環芳烴的組分分布特征各不相同，在東西和南北方向上呈現出明顯的增減規律性。

BERSET J D, HOLZER R, HANI H. 1998. Solving coelution problems on a semctic——liquid crystalline polysiloxane capillary column for the determination of priory polycyclic zromatic hydrocarbons in environmental samples [J]. Journal of Chromatography A, 823(1-2): 179-187.

CARLON C, CRITTO A, MARCOMINI A, et al. 2001. Risk based characterization of contaminated industrial site using multivariate and geostatisti-cal tools [J]. Environmental Pollution, 111(3): 417-427.

CHEN L G, RAN Y, XING B S, et al. 2005. Contents and sources of polycyclic hydrocarbons and organochlorine pesticides in vegetable soils of Guangzhou, China [J]. Chemosphere, 60(7): 879-890.

COLIN P, FROIDEVAUX R, GARCIA M, et al. 1996. Integrating geophysical data for mapping the contamination of industrial sites by polycyclic aromatic hydrocarbons: a geostatistical approach [J]. ASTM Special Technical Publication, 1270: 69-87.

EDWARDS N T J. 1983. Polycyclic aromatic hydrocarbons (PAH’s) in the terrestrial environment: a review [J]. Environ Qual, 12(4): 427-441.

FU J M, MAI B X, SHENG G Y, et al. 2003. Persistent organic pollutants in environment of Pearl River Delta, China: an overview [J]. Chemosphere, 52(9): 1411-1422.

HARRISON R M, SMITH D J T, LUHANA L. 1996. Source appointment of atmosphere polycyclic aromatic hydrocarbons collected from an urban location in Birmingham, UK [J]. Environmental science & Technology, 30(3): 825-832.

HONDA K, MIZUKAMI M, UEDA Y, et al. 2007. Residue level of polycyclic aromatic hydrocarbons in Japanese paddy soils from 1959 to 2002 [J]. Chemosphere, 68(9): 1763-1771.

MALISZEWSKA-KORDYBACH B. 1996. Po1ycyclic aromatic hydrocarbons in agricultural soils in Poland: preliminary proposals for criteria to evaluate the level of soil contamination [J]. Applied Geochemistry, 11(1-2): 121-127.

MALIZEWSKA-KORDYBACH B. 1996. Polycyclic aromatic hydrocarbons in agricultural soils in Poland: preliminary proposals for criteria to evaluate the level of soil contamination [J]. Applied Geochemistry, 11(1-2): 121-127.

MCCREADY S, SLEE D J, BIRCH G F, et al. 2000. The distribution of polycyclic aromatic hydrocarbons in surficial sediments of Sydney harbor, Australia [J]. Marine pollution bulletin, 40(11): 999-1006.

MIELKE H W, WANG G, GONZALES G R, et al. 2001. PAH and metal mixtures in New Orleans soils and sediments [J]. The Science of Total Environment, 281(1-3): 217-227.

SIMCIK M F, EISENREICH S J, LIOY P J. 1999. Source apportionment and source/sink relationships of PAHs in the coastal atmosphere of Chicago and Lake Michigan [J]. Atmospheric Environment, 33(30): 5071-5079.

THOMAS D. BUCHELI, FRANZISKA BLUM, ANDRE DESAULES, et al. 2004. Polycyclic aromatic hydrocarbons, black carbon,and molecular markers in soils of Switzerland [J]. Chemosphere, 56(11): 1061-1076.

TRAPIDO M. 1999. Polycyclic Aromatic Hydrocarbons in Estonian soil: contamination and profiles [J]. Environmental pollution, 105(1): 65-74.

陳來國, 冉勇, 麥碧嫻, 等. 2004. 廣州周邊菜地中多環芳烴的污染現狀[J]. 環境化學, 23(3): 341-344.

段永紅, 陶澎, 王學軍, 等. 2005. 天津表土中多環芳烴含量空間分布特征與來源[J]. 土壤學報, 42(6): 942-947.

高俊峰, 曹慧. 2002. GIS在土壤空間分析中的應用[J]. 土壤, 34(4): 206-209.

黃超勝. 2012. 貴嶼及周邊地區農業土壤中多環芳烴的空間分布研究和生態風險評價[D]. 廣州: 暨南大學.

馬瑾, 邱興華, 周永章, 等. 2011. 東莞市農業土壤多環芳烴污染及其空間分布特征研究[J]. 北京大學學報(自然科學版), 47(1): 149-158.

馬瑾, 周永章, 萬洪富. 2012. 基于多元地統計和GIS的湛江市土壤多環芳烴污染狀況及空間分布特征研究[C]//2012中國環境科學學會學術年會論文集(第三卷): 2613-2625.

彭馳, 王美娥, 廖曉蘭. 2010. 城市土壤中多環芳烴分布和風險評價研究進展[J]. 應用生態學報, 21(2): 514-522.

彭馳. 2009. 北京市土壤多環芳烴分布特征與風險評價[D]. 長沙: 湖南農業大學.

孫紅文, 李書霞. 1998. 多環芳烴的光致毒效應[J]. 環境科學進展, 6(6): 1-11.

王政權. 1999. 地統計學及其在生態學中的應用[M]. 北京: 科學出版社: 35-149.

岳敏, 谷學新, 鄒洪. 2003. 多環芳烴的危害與防治[J]. 首都師范大學學報(自然科學版), 24(3): 40-44.

張天彬, 萬洪富, 周健民, 等. 2008. 深圳表層土壤中多環芳烴的污染特征及來源[J]. 生態環境, 17(3): 1032-1036.

張天彬, 楊國義, 萬洪富, 等. 2005. 東莞市土壤中多環芳烴的含量、代表物及其來源[J]. 土壤, 37(3): 265-271.

鄭太輝, 冉勇, 陳來國. 2014. 東江流域農村土壤中多環芳烴的分布特征及其健康風險評估[J]. 生態環境學報, 23(4): 657-661.

鄭一, 王學軍, 李本綱, 等. 2003. 天津地區表層土壤多環芳烴含量的中尺度空間結構特征[J]. 環境科學學報, 23(3): 311-316.

A study on soil PAHs spatial distribution characters of Huizhou City based on GIS [J]. Ecology and Environmental Sciences, 26(4): 700-707.

WANG Ke, YANG Guoyi, WANG Yulang, ZHANG Tianbin, WAN Hongfu. 2017.

A study on Soil PAHs Spatial Distribution Characters of Huizhou City Based on GIS

WANG Ke1,3, YANG Guoyi1, WANG Yulang2, ZHANG Tianbin1, WAN Hongfu1*
1. Guangdong Institute of Eco-environment Sciences & Technology, Guangzhou 510650, China;
2. Central Station of Soil and Water Conservation Monitoring of Pearl River Basin, Guangzhou 510611, China;
3. College of Natural Resources and Environment, South China Agricultural University, Guangzhou 510642, China

In order to understand the spatial distribution characteristics of polycyclic aromatic hydrocarbons (PAHs) in soils of Huizhou, 42 soil samples of different land use types in different counties were collected. The spatial distribution characteristic of polycyclic aromatic hydrocarbons (PAHs) in soils was analyzed by using geographic information system (GIS) to analyze the residual status, source and diffusion trend of polycyclic aromatic hydrocarbons (PAHs) and provide theoretical support for pollution control of PAHs in Huizhou City. The results showed that the content of PAHs in soils ranged from 35.40 to 534.5 μg·kg-1with the average value of 123.09 μg·kg-1. The results of spatial distribution analysis showed strong spatial correlation between Chr and Bkf; moderate spatial correlation was occurred among Nap, Ace, Acy, Flu, Phe, Ant, Fla, Baa, Daa and Bbf. In contrast, the spatial correlation among Pyr, Bap, I1p and Bgp was very weak. The trend analysis showed obvious regularity for concentrations of 16 polycyclic aromatic hydrocarbons (PAHs). For instance, concentrations of Nap and Bkf in Huizhou soils gradually increased from east to west and from south to north, while Flu, Pyr, Baa, Bbf, Ipy and Daa were opposite; Acy, Ace, Ant increased acutely from east to west and from south to north, but Phe and Fla were opposite to them. No change was observed for Chr and Bap in both of the two directions. Concentrations of Bgp in Huizhou soils decreased gradually from east to west, and showed a clear inverted U-shape in the north-south direction. Ordinary Kriging interpolation was carried out using GIS software. The results showed that the distribution characteristics of 16 polycyclic aromatic hydrocarbons (PAHs) were different. The concentrations of Ace, Acy and Ant increased from northwest to southeast, but the reverse trend was observed for Phe, Fla, Pyr, Bkf, Chr, I1p, Daa and Bgp. No obvious trend occurred for Nap, Bbf and Bap. Highest content area of Nap was very small, but high content area of Bap was relatively larger. The concentrations of Baa in all soil samples of Huizhou were low, but a ‘hot-pot’ area was observed in Huicheng District. The above results reflected different sources of PAHs in different regions. The sources of PAHs are usually related to the economic and social development, and industrial structure in local area, mainly including the emission of “three wastes”, traffic exhaust, oil burning, waste incineration and so on.

soil; polycyclic aromatic hydrocarbon; spacial distribution; GIS; Huizhou City

10.16258/j.cnki.1674-5906.2017.04.023

X53; X144

A

1674-5906（2017）04-0700-08

王珂, 楊國義, 王玉瑯, 張天彬, 萬洪富. 2017. 基于GIS的惠州市土壤多環芳烴空間分布特征研究[J]. 生態環境學報, 26(4): 700-707.

國土資源部農業地質與生態地球化學調查項目（基[2005]011-16）；國家重點研究計劃973子項目（2002CB410810）；廣東省科技計劃項目（2010B031800015）

王珂（1986年生），女，碩士研究生，研究方向為土壤污染控制及修復。E-mail: wangke@soil.gd.cn

*通信作者：萬洪富（1944年生），男，研究員，博士生導師，研究方向為土壤污染控制及修復。E-mail: hfwan@soil.gd.cn

2016-09-20