H11第四次论文研读报告嫌犯围捕.docVIP

城市表层土壤重金属污染分析（特等奖） 主要方法：第一问：地累积指数；克里金插值；ARCGIS软件。 图1各种区的分布图 代码：clear clear,clc A=xlsread(地点及功能区.xls); s1=A([1:44],:) s2=A([45:80],:) s3=A([81:146],:) s4=A([147:284],:) s5=A([285:319],:) plot(s1(:,1),s1(:,2),o,s2(:,1),s2(:,2),r*,s3(:,1),s3(:,2),gs,s4(:,1),s4(:,2),k+,s5(:,1),s5(:,2),m.) hold on legend(生活区,工业区,山区,交通区,公园区) hold on xi=linspace(min(A(:,1)),max(A(:,1)),500); yi=linspace(min(A(:,2)),max(A(:,2)),500); x=(A(:,1)); y=(A(:,2)); z=(A(:,3)); [X,Y,Z]=griddata(x,y,z,xi,yi,cubic) contour(X,Y,Z) 第一问： 4.3 用污染负荷指数的分析各功能区的污染程度 地积累指数分析了每种重金属的污染分布情况，为了评定不同功能区的重金属的综合污染程度，定义最高污染物指数： （4.2） 式中，为第i种金属浓度的实测值，第i种金属的背景值。 如果某点由多种重金属共同污染，该点的污染负荷指数为： （4.3） 为单个点的污染负荷指数，n为金属污染物种数。 由式（4.3），对有单点组成的区域来说，污染负荷指数为： （4.4） 为区域的污染负荷指数，m为采样点个数。 度如表4-3。 表4-3 各功能区污染程度 生活区 工业区 山区 主干道区 公园绿地区 污染负荷指数 1.83 2.35 1.06 1.92 1.58 中等污染 强污染 中等污染 中等污染 中等污染 ，观测值分别为和，那么相关系数可以表示为： 值向量间的相关性 As (μg/g) Cd (ng/g) Cr (μg/g) Cu (μg/g) Hg (ng/g) Ni (μg/g) Pb (μg/g) Zn (μg/g) As (μg/g) 1.000 .199 .440 .416 .359 .538 .281 .306 Cd (ng/g) .199 1.000 .420 .682 .503 .306 .774 .754 Cr (μg/g) .440 .420 1.000 .632 .374 .759 .491 .587 Cu (μg/g) .416 .682 .632 1.000 .567 .538 .754 .792 Hg (ng/g) .359 .503 .374 .567 1.000 .316 .587 .580 Ni (μg/g) .538 .306 .759 .538 .316 1.000 .329 .482 Pb (μg/g) .281 .774 .491 .754 .587 .329 1.000 .818 Zn (μg/g) .306 .754 .587 .792 .580 .482 .818 1.000 重现之后数据没有差别。 可以看到相关性比较大的元素有：Zn和Pb（0.818），Zn和Cd（0.754），Pb和Cd（0.774），Cu和Pb（0.754），Cu和Cd（0.682），Cu和Cr（0.792），也就是说，Zn、Pb、Cd、Cu之间的两两相关非常密切，相关系数均在0.6以上。Cr和Ni的相关性比较大，达到了0.759。As和 Hg是两种比较独立的元素，与其他几种因素的指数的空间分布相关性不是非常显著。 3. 对标准化浓度数据的相关分析 本题给出的数据中，不同的金属浓度变化范围相差很大，为了使不同的数据能在同一标准上比较，需要对原始的浓度数据进行标准化。标准化是将变量中的观察值减去该变量的平均数，然后除以该变量的标准差。即 标准化浓度数据的相关性分析 值向量间的相关性 As (μg/g) Cd (ng/g) Cr (μg/g) Cu (μg/g) Hg (ng/g) Ni (μg/g) Pb (μg/g) Zn (μg/g) As (μg/g) 1.000 .255 .189 .160 .064 .317 .290 .247 Cd