【正文】 程實現(xiàn) (附錄 2).從圖 2 中，可以看出各區(qū)域在該地的地形分布，通過此圖較容易分析 8 種重金屬在該地的空間分布，以及不同區(qū)域的污染程度，為問題 1 的求解做好準備 . 6 模型的建立與求解 問題 1 的分析與求解 [1] 根據(jù)問題分析及模型準備，為了比較容易的觀察出 8 種重金屬的在城區(qū)分布，分別把每種重金屬在采樣點的濃度與各區(qū)域的劃分畫在同一個圖上 . 編程實現(xiàn) (附錄 3). 8 種重金屬的空間分布等高線圖如圖 3： ? 生活區(qū) ? 工業(yè)區(qū) □ 山區(qū) ? 交通區(qū) 公園綠地區(qū) 6 As Cd Cr Cu Hg Ni 7 Pb Zn 圖 3 8種重金屬在該城區(qū)的分布圖 從圖 3 中可以看出對該地污染較大的重金屬有：砷、鎘、鉛、鋅；對該地造成中等污染的重金屬有：汞；對該地污染較輕的重金屬有：鉻、銅、鎳 .其中，在生活區(qū)，污染較嚴重的重金屬是：鎘、鉻、銅、鋅；污染較輕的重金屬是：砷、汞、鎳、鉛 .在工業(yè)區(qū)，污染較嚴重的重金屬是：汞、銅、鋅；污染較輕的重金屬是：砷、鎘、鉛、鉻、鎳 .在山區(qū)，各種金屬對其的影響都很小，其中，相對大一點的是：鎳 .在交通區(qū)，污染較嚴重的重金屬是：汞、銅、鋅、鎘，污染較輕的重金屬是：砷、鉻、鎳、鉛 .在公園綠地區(qū)， 污染較嚴重的重金屬是：汞、鎘、鎳、鋅，污染較輕的重金屬是：砷、鉻、銅、鉛 . 上述結果只是從圖中觀察得出，可信度不算很高，下面用污染程度的單項污染指數(shù)評價法來評價各種金屬在各區(qū)域的污染程度 . . 污染程度的單項污染指數(shù)評價法 [2] 要知道各種重金屬在不同功能區(qū)中的污染程度，首先，要知道各種重金屬在不同區(qū)域的采樣點中的實測值 .根據(jù)實測值與背景值，就可以用 單因子指數(shù)法 來評價各種重金屬在不同區(qū)域中的污染程度，其計算公式： =CiPiSi ，其中， Pi 表示重金屬的 污染指數(shù) ，Ci 表示重金屬的實測值， Si 表示重金屬的評價標準 . 對照《土壤環(huán)境質量標 重 準》 ，土壤重金屬污染評價的分級標準為 ， Pi 1:未污染 。 Nemerow pollution index method 。 Correlation analysis。1 2Pi? :輕度污染 。v439。v439。 [X,Y,Z]=griddata(x,y,z,linspace(0,27823)39。 figure,surf(X,Y,Z) 附錄 2 利用 excel 把采樣點位置數(shù)據(jù)按不同區(qū)域分開，分別存為 data1， data2， data3， data4， x,y,z 與附錄 1 中表示的數(shù)據(jù)相同 . x1= data1(:,1)。,linspace(0,18449),39。*39。 scatter3(x2,y2,z2,39。 z3=data3(:,3)。 y4=data4 (:,2)。) ,hold on x5= data5(:,1)。v39。)。 z1=data1 (:,3)。 y2=data2 (:,2)。),hold on x3= data3(:,1)。ks39。 scatter3(x4,y4,z4,39。 z5=data5(:,3)。 maxiP =。 maxiP =。 maxiP =。附件 239。fontstyle:normal。 /SHEET=name 39。fontsize:13pt。y= data(:,2)。v439。plot3(x,y,a,39。 for i=1:n x1(i)=sum(x0(1:i))。 end b(:,2)=1。 for k=1:n c=x0(1)au(2)/au(1)。 end disp(uint16(yc0(2:1:n+1)))； 。 end yc0(1)=x0(1)。 au=b\y。 end for i=1:n1 b(i,1)=z(i+1)。) 附錄 8 由于 8 種重金屬的高斯擴散模型的程序都差不多，下面以 As 的高斯擴散模型為例 . x0= DATA(:,1)。 mesh(X,Y,A)。 [X,Y,A]=griddata(x,y,a,linspace(0,27823)39。fontweight:normal。 /CELLRANGE=full /READNAMES=on /ASSUMEDSTRWIDTH=32767. DATASET NAME DataSet1 WINDOW=FRONT. FACTOR /VARIABLES V1 V2 V3 V4 V5 V6 V7 V8 /MISSING LISTWISE /ANALYSIS V1 V2 V3 V4 V5 V6 V7 V8 /PRINT INITIAL EXTRACTION /PLOT EIGEN /CRITERIA FACTORS(3) ITERATE(25) /EXTRACTION PC /ROTATION NOROTATE /SAVE REG(ALL) /METHOD={color:0。textdecoration:none} 附錄 6 GET DATA /TYPE=XLS /FILE=39。fontfamily:Monospaced。D:\課件 \畢業(yè)論文 \A\濃度表 .xls39。 maxiP =。 maxiP =。v39。) ,hold on x5= data5(:,1)。 y4=data4 (:,2)。 z3=data3(:,3)。 scatter3(x2,y2,z2,39。*39。 hold on x1= data1(:,1)。,linspace(0,18449),39。 z5=data5(:,3)。 scatter3(x4,y4,z4,39。ks39。),hold on x3= data3(:,1)。 y2=data2 (:,2)。)。 z1=data1 (:,3)。v439。 y=data (:,2)。v439。 3Pi? :嚴重污染 . 下面用 Excel 統(tǒng)計得出不同功能區(qū)中重金屬元素的含量平均值，并將得出的數(shù)據(jù)作為 Ci ， 8 種重金屬的背景值分別為： ， ， 31， ， ， ， 31， 69 (mg/kg)，把它們作為 Si .得到各種金屬在不同區(qū)域中的污染程度，如表 2. 表 1 各種金屬不同區(qū)域中的平均值 (mg/kg) 元素 As Cd Cr Cu Hg Ni Pb Zn 生活區(qū) 工業(yè)區(qū) 山區(qū) 交通區(qū) 公園綠地區(qū) 表 2 各種 金屬不同區(qū)域中的污染指數(shù) 元素 As Cd Cr Cu Hg Ni Pb Zn 生活區(qū) 工業(yè)區(qū) 山區(qū) 交通區(qū) 公園綠地區(qū) 從表 1 中，可以得出，在生活區(qū)， Cu、 Zn 屬于嚴重污染， Hg 、 Cd、 Cr、 Pb 屬于中度污染， As、 Ni 屬于輕度污染；在工業(yè)區(qū)， Hg、 Cu、 Zn、 Cd、 Pb 屬于嚴重污染， As 屬于中度污染， Cr、 Ni 屬于輕度污染；在山區(qū)，污染較輕， 8 種重金屬都屬于輕度污染；在交通區(qū)， Hg、 Cu、 Zn 屬于嚴重污染， Cd、 Pb 屬于中度污染， Cr、 As、 Ni 屬于輕度污 9 染；在公園綠地區(qū)， Hg 屬于嚴重污染， Cu、 Zn、 Cd 屬于中度污染， Pb、 As、 Cr、 Ni 屬于輕度污染 . 模型 [3] 由于單因子指數(shù)評價法只能估計每種重金屬在某一區(qū)域的污染程度，但如果要評價某一區(qū)域的綜合污染程度，就要把各種重金屬在同一區(qū)域的污染指數(shù)結合起來，下面就用內梅羅指數(shù)評價法來評價五個功能區(qū)的污染程度 . 綜合污染評價不僅兼顧了單因子污染指數(shù)平均值和最高值，還可以突出污染較重的重金屬污染物的作用 .其計算公式為：,2 2m a x2 iPPP ?? ）（綜 其中 P綜 是采樣點的綜合污染指數(shù)； maxiP 為 i 采樣點重金屬污染物單項污染指數(shù)中的最大值； ???ni iPnP 11 為單因子污染指數(shù)平均值 .采用《土壤環(huán)境質量標準》中國家一級標準作為環(huán)境質量的標準， 8 種重金屬的標準值分別為： 15， ， 90，35， ， 40