城市表層重金屬污染分析

劉學飛 魏 嬌 何 娟 占文雄

（重慶三峽學院數學與統計學院，重慶萬州 404100）

隨著工業發展和城市化進程的加劇，通過交通運輸、工業排放和市政建設等造成城市重金屬污染越來越嚴重.對城市土壤地質環境異常的查證，以及如何應用查證獲得的海量數據資料開展城市環境質量評價，研究城市不同功能區表層土壤重金屬污染特征和污染分布狀況，以便更好地研究人類活動影響下城市地質環境的演變模式，是具有理論價值和應用價值的課題.本文就如何應用查證獲得的有關8種主要重金屬污染的海量數據資料展開城市環境質量評價，研究某城市重金屬污染的分布狀況、污染評價模型以及傳播模型.

1 重金屬污染程度的空間分布及污染程度分析

1.1 重金屬污染程度的空間分布

利用 Matlab軟件編程，[1-3]制作三維網格圖、三維曲面圖、等高線圖和散點圖等可以分析重金屬污染的空間分布.

關于As：從空間三維圖圖2中可以看到，1處的波峰很高說明該處污染情況很嚴重；有2-4處明顯波峰，說明相應位置污染情況比較嚴重；還有1處面積較廣且所處高度稍低：表明該處所受污染情況相對嚴重且污染的范圍較廣；同樣分析二維等高線圖1，圖中有2處等高線之間的間距越來越密集且顏色很深，表明該處受污染情況很嚴重，有2-4處等高線比較密集顏色相對較深表明這兩處的污染情況相對嚴重，還有1處等高線間的距離較密集但是所涉范圍較廣說明該處的污染也較嚴重且污染的面積很廣.結合有關背景數據可知中心污染源集中在（5291，5739），（12696，3024），（17814，10707）處.污染源主要集中在工業區，還有1處污染源位于山林密集區，其污染級別不是特別嚴重.

重金屬Cd、Cr、Cu、Hg、Ni、Pb、Zn的污染程度可作類似分析.

1.2 重金屬污染程度指標的無量綱化與污染分布狀況的統計分析

由于城區8種重金屬元素污染程度分布各不相同，各重金屬元素污染程度的背景值差異也較大，利用地積累指數（Muller指數）[4]是一種很好的選擇. ?= log2[F n/(C × Bn )]，式中Fn表示污染物重金屬元素n的濃度；Bn表示第n種重金屬元素的背景值上限，C為考慮各地地質差異可能會引起背景值的變動而所取的一系列系數（一般取值為1.5），用來表征沉積特征、巖石地質及其他影響.Muller地積累指數評價和分級標準，[4]詳見表1.

圖1 重金屬As分布平面圖

圖2 重金屬As含量分布的空間三維圖

表1 地積累指數分級標準

利用地積累指數，計算各樣本點每種重金屬的地積累指數，判斷他們的污染級別并進行統計分析，從而可得到每種重金屬元素污染的分布狀況.

表 2列出了重金屬 As污染程度統計分析數據.表中數值0表示在該污染級別下不存在觀測樣本點.這是個大樣本事件，可以認為該級別污染很輕微，甚至不存在這種級別的污染，而百分比越大，就說明在該污染級別下涉及的樣本點比較多，污染波及范圍較廣.

表2 As污染程度分布統計表

從表2可以看出各個區域受As的污染程度，其中生活區31.82%的樣本點無污染，63.64%輕度—中度污染，4.55%為中等污染，無強污染和極嚴重污染；工業區3.89%不受該重金屬污染，52.78%受輕度—中度污染，5.56%受中等污染，2.78%受中等—強污染；山區大多數樣本點不受污染，只有15.15%受輕度—中度污染，1.51%受中等污染；主干道路區47.83%樣本點不受污染，50.00%受輕度—中度污染，0.72%受中等污染，1.45%受中等—強污染；公園綠地區大多數受輕度—中度污染，25.71%不受污染，2.86%受中等污染.這種元素對生活區的影響來說相對較輕，工業區受污染最嚴重，污染面積達到了61.11%，其次是生活區、主干道路區，生活區污染面積都達到了50%以上，也就是說這三個區有至少一半的土壤受到該元素的不同程度的污染.其余功能區受污染程度次之.類似的，可以對重金屬Cd、Cr、Cu、Hg、Ni、Pb、Zn污染分布情況進行統計分析.

1.3 內梅羅綜合污染指數與各功能區重金屬污染狀況分析

通過1.1的分析，可粗略地判斷哪幾種元素污染比較大，哪個功能區污染比較嚴重，但是怎樣判斷哪個功能區污染最嚴重，被污染的功能區的土壤哪種重金屬污染最嚴重，所以，不能簡單地進行數據處理，可用兩種方法分別分別分析處理，互相印證互為補充.一是用地累積指標統計分析法分析各功能區污染狀況和污染程度，二是用下面的內梅羅綜合污染指數法.

單因子指數法χn=αn/βn是目前國內土壤重金屬的單項污染指數評價方法之一，χn為重金屬污染物n的環境質量指數；αn為重金屬污染物n的實測值；βn為重金屬污染物的評價標值.χn＞1表示污染；χn=1或χn＜1表示無污染；且χn值越大，則污染越嚴重.為了更全面地反應各重金屬對土壤的污染程度.突出高濃度重金屬對環境質量的影響，采用內梅羅綜合污染指數法.其計算公式為:

式中（αn/βn）max表示重金屬污染指數αn/βn的最大值，（αn/βn）wg表示重金屬污染指數αn/βn的平均值，土壤污染水平分級標準采用國家土壤環境二級標準.土壤污染綜合污染指數分級標準為：綜合污染指數>3為重污染，2~3為中污染，1~2為輕污染，0.7~1為警戒級，≤0.7為安全級.

表6 各重金屬元素綜合污染指標

表6列示了各重金屬的內梅羅綜合污染綜合污染指數，表7列示出各重金屬的在各功能區的內梅羅綜合污染指數，以及各功能區綜合污染加權平均指標，很容易看出哪種重金屬元素污染比較嚴重，或者在哪個功能區的污染比較嚴重.

通過分析，發現污染最嚴重的功能區是主干道路區，其次按照污染程度從大到小依次排序為：工業區、生活區、公園綠地區、山區.之所以主干道路區污染最為嚴重是因為該區土壤表層重金屬元素含量很高，且種類多.

2 重金屬污染的主要原因

根據數據統計分析（并參考文獻[5]-[9]），Hg污染最為嚴重的是主干道路區.公路、鐵路兩側土壤中的重金屬污染，以Pb、Zn、Cd、Cr、Co、Cu等的污染為主，它們來自于含鉛汽油的燃燒，汽車輪胎磨損產生的含鋅粉塵等.隨著時間的推移，公路、鐵路土壤重金屬污染具有很強的疊加性.這也是主干道路區各種重金屬污染都很嚴重的原因.

其次，工業區的重金屬污染也很嚴重，主要是因為經過自然沉降和雨淋沉降進入土壤的重金屬污染，主要以工礦煙囪、廢物堆和公路為中心，向四周及兩側擴散；由城市→郊區→農區，隨著距城市的距離加大而降低，特別是城市的郊區污染較為嚴重.此外，大氣汞的干濕沉降，也可以引起土壤中汞的含量增高.大氣汞通過干濕沉降進入土壤后，被土壤中的粘土礦物和有機物的吸附或固定，富集于土壤表層，或為植物吸收而轉入土壤，造成土壤汞的濃度的升高.還由于部分礦產開發中選礦、冶煉工藝水平落后，個別礦區沒有環保治理設備，廢水、廢氣排放而帶來的大量廢棄物的產生未經處理直接投放環境，而其中的重金屬也隨著自然的沉降、雨水的淋溶等途徑進入土壤，進入正常循環的生態系統，造成重金屬污染，嚴重危害人們的生產生活.

生活區的重金屬污染主要來自生活廢水的排放，化肥農藥的多度使用，化肥中品位較差的過磷酸鈣和磷礦粉中含有微量的As、Cd重金屬元素.含鉛及有機汞的農藥發揮作用的同時也為土壤重金屬污染埋下了禍根，造成土壤的膠質結構改變，營養流失，對農作物的產量及品質都造成極大的不良影響.目前的飼料添加劑中也常含有高含量的Cu和Zn，這使得有機肥料中的Cu、Zn含量也明顯增加，并隨著肥料施入農田.

公園綠地區的重金屬污染主要來自城市綠地區，施用肥料如磷肥等也會帶來Hg的嚴重污染.

山區的重金屬污染主要是來自金屬礦山的開采、冶煉、重金屬尾礦、冶煉廢渣和礦渣堆放等，可以被酸溶出含重金屬離子的礦山酸性廢水，隨著礦山排水和降雨使之帶入水環境（如河流等）或直接進入土壤，都可以間接或直接地造成土壤重金屬污染.

總的來說，隨著城區工業化的進程越來越快，各種重金屬的排放會越來越多，嚴重影響人們的生活健康和環境.

3 重金屬污染傳播模型

3.1 重金屬污染物的傳播特征

根據問題一中建立的重金屬元素分布等高線圖，能夠了解城區重金屬的分布規律.它們是由污染源向四周逐漸擴散的，它們的擴散方式主要有：[5-9]

1）大氣中的重金屬主要來源于工業生產、汽車尾氣排放及汽車輪胎磨損產生的大量含重金屬的有害氣體和粉塵等.它們主要分布在工礦的周圍和公路、鐵路的兩側.大氣中的大多數重金屬是經自然沉降和雨淋沉降進入土壤的.

2）公路、鐵路兩側土壤中的重金屬污染，主要是Pb、Zn、Cd、Cr、Co、Cu的污染為主.它們來自于含鉛汽油的燃燒，汽車輪胎磨損產生的含鋅粉塵等.它們成條帶狀分布，以公路、鐵路為軸向兩側重金屬污染強度逐漸減弱；隨著時間的推移，公路、鐵路土壤重金屬污染具有很強的疊加性.

3）經過自然沉降和雨淋沉降進入土壤的重金屬污染，主要以工礦煙囪、廢物堆和公路為中心，向四周及兩側擴散；由城市→郊區→農區，隨距城市的距離加大而降低，特別是城市的郊區污染較為嚴重.

4）大氣汞的干濕沉降也可以引起土壤中汞的含量增高.大氣汞通過干濕沉降進入土壤后，被土壤中的粘土礦物和有機物的吸附或固定，富集于土壤表層，或為植物吸收而轉入土壤，造成土壤汞的濃度的升高.

5）在河流中不同河段的重金屬污染往往受污染源（礦山）控制，河流同一污染源的下段自上游到下游，由于金屬元素遷移能力減弱和水體自凈化能力的適度恢復，金屬化學污染強度逐漸降低.

3.2 重金屬污染的微分方程傳播模型

3.2.1 考慮煙霧均勻熱傳導，不考慮大氣流動，有微分方程模型：

初始條件為v(x,y,z,0)=Qδ(x,y,z)，v(x,y,z,t)是t時刻在(x,y,z)點處重金屬污染物的濃度，為重金屬元素污染程度的變化規律，Q表示重金含量.

3.2.2考慮吸收-擴散，考慮重金屬污染物從高濃度處向低濃度處擴散，一部分重金屬污染物由于外界因素進入地下被土壤吸收，擴散系數與吸收系數都是常數，設v(x,y,z,t)是t時刻在(x,y,z)點處重金屬污染物的濃度，w(x,y,z,t)是t時刻在(x,y,z)點處被土壤吸收的重金屬污染物的濃度.又設重金屬污染物在各個方向的傳播系數為k，吸收系數h，于是有微分方程模型：

通過變換，上述方程有經典解，但涉及若干待定常數.

4 尋找污染源的向內收索算法

上述微分方程模型，盡管有經典解，但某些參數無法直接確定.為了尋找污染源，可以考慮從污染濃度等高線（污染濃度作立標）污染邊界向內搜尋污染源，并簡化為從等高線節點向內搜索污染源，逐步減少搜索半徑，最終確定污染源.考慮如下算法：

1）k=0，任取初值點（包括二維坐標、濃度）；

2）搜索初值點附近濃度相同的樣本點至少兩個，尋找等高線；

3）擬合前述所得樣本點的等高線，并向內搜索樣本點，如果該樣本點濃度增高，則以該樣本點為初始點，重復第2步，否則按第4）；

4）繼續向內搜索樣本點，重復第3步；

5）如果無法在初始樣本點附近找到新的樣本點，使其滿足濃度增高，則該樣本點即為污染源，輸出該樣本點；

6）循環以上步驟，繼續尋找下一個污染源.

5 模型的分析與評價

對于空間分布幾何模型：該模型能夠直觀顯示各種重金屬的空間分布以及各種重金屬的污染程度，但該模型無法從量化的角度分析其污染情況，而且該幾何模型需熟悉相關軟件.

地累積指標模型：模型較簡單，計算方便，且無量綱化，可以用于各種元素的污染程度比較，但數據計算量比較大，而且污染程度指標的分級是否完全合理有待于進一步實證檢驗.

內梅羅綜合污染指數模型：是綜合反映8種重金屬元素的總的污染程度大小的綜合指標，或者同種元素在某個地區污染程度的大小的綜合指標，（αn/βn）max項對指標值影響較大，從公式看若某地區有重大污染源，那么該指標值較大，有時甚至相當大.事實上，計算出的指標有超過 200的情況.反之，如果該指標值較大和很大，說明該地區有嚴重的污染源；但不能具體得出是哪種重金屬元素產生的污染，以及污染源位置，需要地累積指數模型與該模型結合使用.微分方程模型具有很強的分析特征，在實用方面具有一定的局限.

6 模型的改進

重金屬污染濃度的微分方程模型，初步考慮了重金屬在大氣、土壤中的傳播特征，還應考慮重金屬在液體中的傳播特征，水土流失、植物吸收等傳播特征，對模型進行進一步研究和改進.

[1]姜啟源.數學模型：第三版[M].北京：高等教育出版社，2003.

[2]張德豐.MATLAB數值計算方法[M].北京：機械工業出版社，2010.

[3]汪曉銀，周保平.數學建模與數學實驗[M].北京：科學出版社，2010.

[4]陳莉薇，徐曉春，黃界穎，等.銅陵林沖尾礦庫復墾土壤重金屬含量及污染評價[J].合肥工業大學學報：自然科學版，2011（10）.

[5]鄭喜珅，魯安懷，高翔，等.土壤中重金屬污染現狀與防治方法，中國環境生態網，http://www.eedu.org.cn/Article/academia/papers/s umpapers/200404/115.html.

[6]鄭喜珅，魯安懷，高翔，等.土壤中重金屬污染現狀與防治方法[J].土壤與環境，2002（1）.

[7]王濟.貴陽市表層土壤重金屬污染元素環境地球化學基線研究[D].中國科學院研究生院（地球化學研究所），2004.

[8]劉紹貴.南昌市城郊表層土壤重金屬污染特征研究（EB/OL）.http://wenku.baidu.com/view/e8e 14d10f18583d0496459dd.html，2010-04-02.

[9]伊俊，柳云龍.上海市城郊土壤重金屬空間分布及其污染評價（EB/OL）.http://wenku.baidu.com/view/bd85958b 84868762cbaed505.html，2010.