2011年全国大学生数学建模竞赛A题优秀获奖论文—城市表层土壤重金属污染分析

1、城市表层土壤重金属污染分析摘要 本文通过对城市表层土壤重金属污染综进行分析，建立了内梅罗评价模型、内因子分析模型确定重金属对城区的污染情况及抛曲面拟合模型确定污染源位置等一系列数学模型。最后，又提出了对研究城市地质环境的演变模式的可行方法。问题一：建立重金属污染程度处理的内梅罗定量评价模型。导入附件1、附件2中坐标位置，海拔高度及8种金属元素的浓度数据，用Surfer9.0软件依据源文件数据作图，可以分别得到城区区域拟合的区域图，即功能区平面图，及8种金属元素在各城区区域拟合的等高浓度线分布图。从图中清晰的等高线分布，定性分析元素，同时用单因子污染指数，内梅罗综合污染指数法，将该城区内不同区域

2、的重金属的污染程度量化分析，从直观的单因子指数与综合指数中，确定了该城区不同区域重金属污染的程度。问题二：源数据进行内梅罗模型数据处理后，得到评判区域污染的污染指标Pi及由Pi演化的图表，在分块图表中获得出污染元素的区域特性，并进一步利用因子分析模型确定各区域污染的主因子。因子分析模型确定该城区污染原因主要在工业“三废，汽车尾气排放及轮胎磨损。问题三：问题一、二反映出城区污染不仅有点源污染，还有面源污染，结合重金属污染物因雨水冲刷、大气沉降等因素浓度梯度由高到低迁移及海拔由高到低迁移的传播特征，建立抛曲面拟合模型，模型借用物体辐射抛射产生的效果，近似地模拟城区污染程度与地理位置间的关系，通过对

3、最小二乘法修正，得，最终确定有2个可供参考的主要污染源。问题四：对在前三个问题中建立的数学模型做了简要的分析。并对更好地研究城市地质环境的演变模式提出应收集所需研究的元素在土壤中的浓度数据的年度记录，建立扩散模型，扩散方程为：关键词：内梅罗综合污染指数 因子分析法 三维曲面拟合 扩散模型问题一1.1 问题分析问题一是要求分析该城区内不同区域重金属的污染程度，可以根据附件一中列出的采样点的位置、海拔高度以其所属功能区等信息可以用Surfer9.0画出该城区区域图的各分区大致位置，再根据附件二中8种主要重金属元素在采样点处的浓度，可以同样得到这8种重金属分别在该城区中的浓度分布情况。再采用单因子污

4、染指数和综合污染指数法对该城区内不同区域的重金属的污染程度进行评价，从而得出结论。1.2 模型假设（1） 假设取样点的数据能够较好的反映该地区的污染物的浓度；（2） 假设污染源的重金属的浓度在一定的时间内是稳定的。1.3 模型建立和求解通过附件一中列出的采样点位置和海拔高度可以用Surfer9.0拟合出该城区分区的大致位置，如图1.1，图中不同的颜色对应着城区中不同的分区。 图1.1 城区各功能区拟合图通过附件一中采样点的位置坐标和附件二中采样点的各金属的浓度可以用数学软件Surfer拟合出8种金属元素分别对该城区的污染程度图，分别如图1.2图1.9，图中等高线的梯度越大，分布越密集，那么此种

5、金属在此分布范围内就浓度越大。图1.3 Cd浓度分布图图1.2 As浓度分布图 图1.5 Cu浓度分布图图1.4 Cr浓度分布图图1.7 Ni浓度分布图图1.6 Hg浓度分布图图1.9 Zu浓度分布图布图图1.8 Pb浓度分布图布图从图1.2到图1.9只能够直观的大概反映出各种重金属元素对该城区的污染程度，但是还是无法定量的来反映出不同区域重金属的污染程度，因此需要采用单因子污染指数和综合污染指数法对该城区内重金属元素污染状况进行评价，其计算公式为： 式中：为重金属的单因子指数；为重金属的实测浓度。当时表示未受污染，当时表示受到污染，且越大污染越严重。表1.1 城区各种重金属元素平均值重金属元

6、素AsCdCrCuHgNiPbZnPi平均值1.59 1.09 2.70 5.88 0.87 1.44 2.07 从表1.1中可以看出该城区内除了重金属Ni的Pi平均值小于1，其他各金属的Pi平均值都大于1.所以该城区内未受到重金属Ni的污染影响，而其他金属对该城区均有影响，其中重金属Hg对该城区的污染最严重。各重金属元素对城区的污染程度大致为Hg>Cu>Zn>Cd>Pb>Cr>As>Ni。要对该城区内不同区域重金属的污染程度的评价那么需要用到综合污染指数法，其公式为： 式中：为某地区的综合污染指数；为重金属污染物中污染指数最大值,为重金属污染物中污染

7、指数平均值。通过对附录中的重金属浓度和背景值的分析计算，可以得到一区到五区中各种重金属综合作用污染指数，如表1.2所示。 表1.2 城区各分区的各种重金属综合作用污染指数 生活区工业区山区主干道路区公园绿地区P综9.28 0.82 6.46 1.87 表1.3 中国绿色食品开展中心推荐的分级标准进行的污染分级 等级综合污染指数污染程度1平安20.7<P<=1警戒线3轻污染41.5<P<=2中污染52<P<=3重污染6P>3严重污染重金属综合作用污染指数。对照着表1.3的土壤重金属分级标准，可以发现除山区位于警戒线附近，其他的各区均受到不同程度的污染，各

8、区的污染程度依次为工业区>主干道路区>生活区>公园绿地区>山区。工业区和主干道路区的污染程度远远高于其他各区，生活区和公园绿地区的污染程度大致相同，山区的污染程度最小，但是也位于被污染的边缘。 问题二2.1 问题分析利用问题一对所提供污染元素的浓度进行无量纲处理，并取8种元素的调和平均数作为功能区粗评定综合指标。由附录二表建立模型的样本空间，利用样本空间中Pi的调和平均数间的比重关系，确定各功能区占主导因素的污染元素，得到如下百分比图：图表 1图表 2图表 3图表 4图表 5结论与分析从图表中放映出来的信息分析之后，可以得出生活区以Hg，Cu，Pb为主要污染元素，工业区

9、以Hg，Cu为主要污染元素，山区各污染元素相对平等，主干道路区以Hg，Cu为主要污染元素，公园绿地区以Hg，As为主要污染元素。虽然，图表信息表达了区域中起较重要的污染元素，但区域特征及污染原因分析中所得数据集中在Hg，Cu之中，由于调和平均数受极差影响使得所得结果不能让区域性及污染特点量化，故需对该模型进行优化处理。内因子模型建立与求解设有n个样本，n个指标构成样本空间X， ，i = 1,2,n; j = 1,2,m，并对原始数据标准化，标准化公式为，其中为样本中第个样本的第指标值，而和分别为指标的均值和标准差；计算标准化数据的相关系数矩阵，进而消除不同变量的量纲的影响且标准操作数据不会影响

10、变量的相关系数；求出相关系数的特征值和特征向量，进行正交变换，使用方差最大法，使得所得因子载荷两极分化，但不影响旋转后的因子的正交特性。模型通过正交的方差最大旋转法是每个主因子至于最少个数的变量有相关系数，使得足够多的因子负荷均减小，并通过多元统计分析的计算结果，分析及确定主因子。内因子模型结论与分析生活区主因子的求解：由附录一表2-1生活区土壤重金属浓度 得到模型建立的样本空间X，借助统计软件SPSS 11.0计算得出如下表格：表 变量相关联矩阵1As (/g)Cd (ng/g)Cr (g/g)Cu (/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)As (g/g)1

11、.000 Cd (ng/g)0.381 1.000 Cr (g/g)0.238 0.349 1.000 Cu (g/g)0.531 0.499 0.376 1.000 Hg (ng/g)0.293 0.397 0.150 0.198 1.000 Ni (g/g)0.605 0.283 0.527 0.434 0.211 1.000 Pb (g/g)0.450 0.802 0.416 0.502 0.340 0.300 1.000 Zn (g/g)-0.017 0.346 0.412 0.238 0.242 0.334 0.328 1.000 表 特征值和累积奉献率1旋转前旋转后总的特征量占总变

12、量的百分率%累计奉献率%总的特征值占总变量的百分率%累计奉献率%1 23 表 方差极大正交旋转后载荷矩阵1表 旋转前因子载荷矩阵1 因子123As (g/g)Cd (ng/g)Cr (g/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)因子123As (g/g)Cd (ng/g)Cr (g/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g) 由表中所得相关系数观察可得，Ni与As的相关性最好，相关系数为0.605，其次为Cu与As，Ni与Cr相关较好。且在表3-1-2中可以观察得到累积方差为72.797%，分析得到3个主

13、因子，3个主因子提供了源资料的72.797%的信息，以及主因子1与主因子2的方差奉献率均在26.1%左右，从表3-1-3，表3-1-4综合反映出旋转前后载荷的变量结果根本一致，正交因子解说面得出因子1 为Cu与As，因子2为 Ni与Cr，这两组重金属元素为相关性最好的两组元素，可作为评定生活区重金属污染的主因子。工业区主因子的求解：由附录一表3-2 工业区土壤重金属浓度得到模型求解所需的样本空间X，依据生活区主因子的求解方法及如下表格分析求解得，主因子1 Cu与Hg，主因子2 As与Ni相关性最好，可作为工业区重金属污染的主因子。表 变量相关联矩阵2As (g/g)Cd (ng/g)Cr (g

14、/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)As (g/g)1.000 Cd (ng/g)0.329 1.000 Cr (g/g)0.380 0.541 1.000 Cu (g/g)0.153 0.566 0.920 1.000 Hg (ng/g)0.181 0.533 0.902 0.983 1.000 Ni (g/g)0.690 0.489 0.698 0.503 0.479 1.000 Pb (g/g)0.395 0.829 0.675 0.670 0.612 0.578 1.000 Zn (g/g)0.518 0.754 0.695 0.62

15、2 0.590 0.634 0.739 1.000 表 特征值和累积奉献率2旋转前旋转后总的特征量占总变量的百分率%累计奉献率%总的特征值占总变量的百分率%累计奉献率%12表 旋转前因子载荷矩阵2表 方差极大正交旋转后载荷矩阵2 因子12As (g/g)Cd (ng/g)Cr (g/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)因子12As (g/g)Cd (ng/g)Cr (g/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)山区主因子的求解：由附录一表3-3 山区土壤重金属浓度得到模型求解所需的样本空间X，依据山

16、区主因子的求解方法及如下表格分析求解得，主因子1 Ni与Cr、主因子2 Cd与Pb相关性最好，可作山区土壤重金属污染的主因子。表 变量相关联矩阵3As (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 Cd (ng/g)-0.291 1.000 Cr (g/g)0.113 0.066 1.000 Cu (g/g)0.527 0.090 0.364 1.000 Hg (ng/g)0.075 0.246 -0.006 0.505 1.000 Ni (g/g)0.078 0.049 0.945 0

17、.358 -0.045 1.000 Pb (g/g)-0.205 0.766 0.107 0.122 0.226 0.028 1.000 Zn (g/g)-0.176 0.606 0.627 0.252 0.170 0.629 0.590 1.000 表 特征值和累积奉献率3旋转前旋转后总的特征量占总变量的百分率%累计奉献率%总的特征值占总变量的百分率%累计奉献率%123表 方差极大正交旋转后载荷矩阵3表 旋转前因子载荷矩阵3因子123As (g/g)Cd (ng/g)Cr (g/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)因子123As (g/g)C

18、d (ng/g)Cr (g/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)234主干道路区主因子的求解由附录一表3-4 主干道路区土壤重金属浓度得到模型求解所需的样本空间X，依据主干道路区主因子的求解方法及如下表格分析求解得，主因子1 Cr与Cu、主因子2 Pb与Hg相关性最好，可作为主干道路区重金属污染的主因子。表2.2.13 变量相关联矩阵4As (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 Cd (ng/g)0.121 1.000 Cr (g

19、/g)0.139 0.373 1.000 Cu (g/g)0.092 0.424 0.894 1.000 Hg (ng/g)-0.004 0.211 0.012 0.032 1.000 Ni (g/g)0.228 0.351 0.869 0.886 0.040 1.000 Pb (g/g)0.060 0.615 0.428 0.506 0.266 0.396 1.000 Zn (g/g)0.188 0.294 0.395 0.432 0.118 0.503 0.482 1.000 表 特征值和累积奉献率4旋转前旋转后总的特征量占总变量的百分率%累计奉献率%总的特征值占总变量的百分率%累计奉献率

20、%12表6 方差极大正交旋转后载荷矩阵4表 旋转前因子载荷矩阵4因子12As (g/g)Cd (ng/g)Cr (g/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)因子12As (g/g)Cd (ng/g)Cr (g/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)公园绿地区主因子的求解：由附录一表1公园绿地区土壤重金属浓度得到模型求解所需的样本空间X，依据公园绿地区主因子的求解方法及如下表格分析求解得，主因子1 Cr与Cu、主因子2 Pb与Hg相关性最好，可作为公园绿地区重金属污染的主因子。镉Cd镉电镀可以为钢、铁

21、等提供一种抗腐蚀性的保护层, 具有吸附性好且镀层均匀光洁等特点, 因此工业上90%的镉用于电镀、颜料、塑料稳定剂、合金及电池等行业。塑料制造过程中所需的热塑稳定剂的成分中包括镉元素，因此，人类活动的产生的塑料垃圾、燃烧垃圾都会导致Cd元素在进入到大气中，最终沉降到土壤中。同时，在日常生活中电镀生活用品、办公用品等的大量使用至废弃、随意丢放，使得重金属镉进入土壤。在自然界中没有单独的镉矿藏，镉通常与锌共生，镉是铅-锌 矿、铜-铅-锌 矿的伴生元素，锌矿石中含有大量的镉，在锌冶炼的过程中会释放出大量的镉，但农药、医药等生活相关的行业生产过程中产生的工业废水通过在土壤的迁徙、渗透，也导致重金属镉进入

22、土壤。汞Hg汞常温下呈现为液态状，主要来源于照明所用的荧光灯、齿科材料、防“四害药品，并且工业废水中含有大量的Hg元素，虽然经过污水处理，但是其污水中的含量远超过国家污水灌溉的标准，因此，人类生活产生的、人类生产灌溉及工业废水排放、渗透叠加共同造成了土壤的Hg污染。砷AsAs主要来自土壤母质，而生活区和公园绿地区的绿化面积比拟大，一些树木的树龄也比拟大，为了更好的保护树木花草的健康成长，就会常年的施一些含有As的化肥，经过常年的累积就会致使区域的土壤里的As污染。As在自然界分布广，且用途多。例如As的化合物存在陶瓷和玻璃、电子产品，涂料，去污剂，化装品，烟花爆竹等日常生活用品，以及增强金属的

23、防腐性As合金，而当前As化物还正被广泛用做杀虫剂、消毒液、杀菌剂居民在使用这些产品的过程中或其垃圾排放都会使砷的化合物进入城市土壤环境中。此外，各类煤碳中也含有大量的Hg和As，会随着煤碳的燃烧丽进入大气中，最后随大气沉降进入土壤环境中。铅PbPb主要应用与蓄电池、电镀、颜料、橡胶、农药、燃料等制造业，其中铅板制作工艺中排放的酸性废水( pH < 3 )铅浓度最高, 电镀废液产生的废水铅浓度也很高；源于工业“三废排放、污水灌溉和交通运输。其中，燃煤和交通运输可能是其重要的来源交通运输所用各种工具，如汽车、飞机等使用的汽油燃烧后可把含Pb的化合物排入大气，使得机场附近和交通道路两侧的土壤

24、严重污染。铜Cu、锌ZnCu、Zn主要来源于工业“三废排放、污水灌溉，其中交通运输中含铜配件的磨损，轮胎的磨损产生Zn污染，因而造成Cu、Zn土壤污染。铬CrCr主要来自电镀、染料、制药、皮革、颜料等铬化合物制造企业所排放的“三废以及铬矿和金属冶炼、电镀、制革的工业废水、废气和废渣，从而造成Cr的污染。问题三3.1 问题分析 城区重金属污染元素浓度的采用因子分析模型分析处理，结果得出Cr，Ni，Hg，Pb为影响全区污染的主因子，且存在城区影响标高点。问题二反映出污染元素区域性特征具有较强的水溶性特点，如Hg、Cu等，而以Pb、Zn等为主的公路重金属污染物主要分布在主干道路周围，呈现出带状分布特

25、征。故该城区的污染源不只有一个，而有多个；不仅有点源污染，还有面源污染，不可将其单一化看待。再结合重金属污染物因雨水冲刷、大气沉降等自然因素所表现出的按浓度梯度由高到低迁移、按海拔由高到低迁移的传播特征，我们将第一问中得出的各样本点的污染程度综合指数按地理位置坐标进行三维曲面拟合，并作出拟合图形，近似地得到该城区污染程度与地理位置的关系，并确定污染源位置。3.2 抛曲面拟合模型假设1假设问题一中所得出的重金属污染物污染程度综合指数是合理的，能有效地说明样本点的实际污染程度；2假设污染源以面源污染源为主，局部点源污染源集合可用面来表示；3因雨水冲刷，水溶性的重金属污染物易在低洼处富集，其富集程度

26、可能较高，造成二次污染，故假设此类污染物集合为面源污染源，该特性正如物体抛射所成曲线，当物体抛射方向为中心辐射状，将得到曲面，同时在物体抛射过程中，物体在外界情况下的损耗及沉降与该城区污染传播类型相似；4假设该城区范围是可局部延展的，即所拟合出的面源污染源坐标集合可在原附件中所列坐标能够围成的范围之外。3.3 抛曲面拟合模型建立与求解模型确定该城区的重金属污染源，建立抛曲面拟合方程：；其中，x、y分别表示某采样点的横纵坐标，z表示该采样点处的综合污染指数。拟合方程共有10 项，有10 个系数，运算量较大。 Matlab R2021a 中的sftool工具Surface Fitting Tool

27、能够便捷处理、拟合数据，代入相应数据进行拟合运算，解出符合该题的抛曲面拟合方程,该方程拟合所得该城区综合污染指标三维曲面图象及其等高线图取该城区样本中综合污染指标为高，如图3.1和图3.2所示：图3.1 污染源分布三维拟合图象图3.2 污染源分布三维拟合图象等高线图由以上两幅拟合图象中，我们可以看到重金属的污染程度在红色填涂的区域内较高。该城区这样的区域主要有两块，靠近南部为便于观察，图3.1的视角经过调整，与图3.2方位并不对应，这里以图3.2的方位为准，下同的一块较为平缓，说明该处以面污染源为主，而实际上由问题一中图象分析出该处主干道路较多；而位于西南角的另一块重度污染区域曲率较大，说明该

28、处以点污染源为主，而实际上该处有工业区。综上所述，该模型能较好地模拟该城区的污染程度，通过该模型的方程和图象，我们可以确定污染源的位置，并能了解到一定范围的受污染区域内污染源的类型。问题四问题一先通过Surfer软件拟合出城区区域图和城区各金属元素浓度的分布图，可以直观看出各种金属元素对城区的污染情况，然后用单因子污染指数和内梅罗综合污染指数法对问题进行定量的分析，可以更合理,更深刻地来说明问题，得出结论。其中内梅罗综合污染指数法是当前国内外进行综合污染计算的更常用的方法之一，兼顾了单因子污染指数平均值和最高值，突出了较重要的污染物对环境质量的影响。其模型的缺乏之处在于考虑的重金属元素的存在和

29、传播的形式不充分，内梅罗综合指数法主要是考虑土壤中重金属元素的传播，而实际中重金属元素的传播和存在形式是多样的，有可能少量存在于空气中，这些未考虑的因素很可能会对问题的结果产生一定的误差，但误差的影响不大。针对问题二，我们从图表中放映出来的信息分析之后，可以得出生活区以Hg，Cu，Pb为主要污染元素，工业区以Hg，Cu为主要污染元素，山区各污染元素相对平等，主干道路区以Hg，Cu为主要污染元素，公园绿地区以Hg，As为主要污染元素。虽然，图表信息表达了区域中起较重要的污染元素，但区域特征及污染原因分析中所得数据集中在Hg，Cu之中，由于调和平均数受极差影响使得所得结果不能让区域性及污染特点量化

30、，固需对该模型进行优化处理。问题三建立的模型能很好地表达该城区的重金属污染程度分布，且能求解出污染源的大致位置，即点源污染的大致位置。尤其，能通过拟合所得三维图象判断出一定区域内主要包含的污染源类型点源污染源或面源污染源，而且借用Matlab专业数学软件中sftool工具Surface Fitting Tool处理运算数据，极大地提高了运算效率和准确率。但是，在数据分析处理中，当源数据存在污染指数高异常样本点时，拟合结果不能够将该类点剔除，故对全体数据拟合效果较差。不过，这并不影响该模型对所有样本点总体污染情况及污染源大致位置的判断。为更好地研究城市地质环境的演变模式，还应每隔固定的一段时间收

31、集所需研究的元素在土壤中的浓度数据比方说年记录数据。于是，便可利用收集到的数据，建立扩散模型，根据扩散方程；得出某元素的扩散情况，其中c为该元素的浓度值，和为该元素在x、y方向上的扩散系数，S为源汇项。该方程主要适用于人为排放的污染物的扩散迁移的研究，也适用于某特定地点如矿脉的土壤中富集元素的扩散迁移等研究，能一定程度地反映城市地质环境的演变模式。参考文献1 王向东.数学实验(第一版)M.北京.高等教育出版社 2006年2   3 霍晓程，李小平.用最小二乘法拟合曲面方程.赤峰学院学报J.第6期：11-13，2021年6月4 谭和平，陈能武，黄苹，叶善蓉. 川西茶

32、区土壤重金属元素背景值及其评价. 农业环境科学学报J. 2005, 24( 增刊): 196-1995 王树禾.扩散问题的偏微分方程模型. 6 陈玲，夏俊，李宇庆，仇雁翎，赵建夫，李金柱. 上海化学工业区土壤环境背景值调查.生态学杂志.2005，241：65-692附录附录一表3-1生活区土壤重金属浓度功能区As (g/g)Cd (ng/g)Cr (g/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)19.17 287.00 43.94 45.77 168.00 19.70 62.74 223.16 15.72 193.70 80.35 26.57 111

33、.00 19.80 57.64 89.08 111.45 1044.50 94.78 136.97 202.00 22.30 472.48 602.04 17.84 347.90 57.65 97.14 213.00 19.60 70.82 307.24 18.50 614.00 744.46 130.55 156.00 32.80 228.64 1013.47 15.51 257.20 54.64 29.01 104.00 13.20 87.68 223.27 19.39 325.80 172.29 104.89 82.00 31.50 90.90 429.29 14.09 90.50 35

34、.02 11.82 16.00 10.40 29.09 46.84 16.35 532.00 57.51 83.76 191.00 19.50 73.46 297.14 13.50 396.30 138.37 58.97 170.00 24.20 91.76 2893.47 14.29 449.10 67.22 51.64 315.00 15.70 106.97 294.69 17.41 337.30 77.27 248.85 90.00 20.10 99.58 210.00 12.91 265.00 35.66 29.39 24.64 9.23 60.54 122.96 13.30 488.

35、00 51.18 34.55 37.50 10.80 54.62 125.92 16.14 227.00 42.15 67.04 49.29 16.31 34.28 82.96 18.06 113.10 52.40 20.81 65.36 19.69 29.56 62.24 14.69 132.80 36.11 17.28 20.36 15.53 37.03 82.76 12.34 87.60 18.46 9.73 13.93 8.89 24.43 43.37 16.56 245.60 36.73 61.30 55.71 14.18 47.24 114.29 18.23 756.40 42.7

36、3 87.52 63.00 19.26 88.74 184.69 18.90 307.30 54.39 57.21 326.00 25.72 131.93 256.94 15.41 178.90 29.54 23.73 52.00 9.89 49.84 118.88 17.78 315.50 49.76 28.03 550.00 18.95 45.73 109.29 15.62 134.60 25.33 19.10 45.00 11.66 40.50 87.14 16.05 365.00 35.92 30.91 110.00 17.91 44.75 147.35 14.17 310.70 40

37、.70 42.64 58.00 14.99 115.59 177.76 16.26 387.20 38.03 30.06 85.00 21.97 54.98 142.55 15.00 196.50 50.03 18.56 29.00 25.82 33.90 84.80 14.58 129.00 31.09 18.93 38.00 15.10 29.76 69.80 15.41 204.90 40.16 25.86 41.00 16.24 33.28 80.00 17.56 206.60 55.79 24.44 22.00 28.63 29.73 81.73 15.41 250.10 34.32

38、 16.25 15.00 18.53 41.29 90.20 15.83 89.20 54.90 23.28 19.00 21.76 24.95 53.57 15.20 334.30 47.05 40.37 76.00 18.40 61.03 211.47 15.20 113.70 41.85 20.73 12.00 19.00 31.87 57.27 14.38 246.30 29.90 21.84 84.00 14.10 59.00 95.17 15.41 158.00 46.86 24.02 31.00 19.30 36.27 79.24 18.67 245.70 47.79 27.55

39、 35.00 18.40 53.42 98.81 16.47 86.80 41.12 15.46 23.00 15.90 37.53 70.18 17.12 367.80 92.02 49.80 97.00 16.30 41.26 321.12 110.97 248.50 40.61 61.52 81.00 17.20 76.19 168.05 19.81 171.80 75.38 163.20 30.00 26.30 45.27 125.16 12.77 236.20 42.67 16.35 62.00 9.40 41.88 149.52 16.47 282.90 52.68 20.34 2

40、5.00 22.90 32.53 103.50 表3-2 生活区土壤重金属浓度功能区As (g/g)Cd (ng/g)Cr (g/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)26.56 223.90 40.08 25.17 950.00 15.40 32.28 117.35 214.08 1092.90 67.96 308.61 1040.00 28.20 434.80 966.73 29.62 1066.20 285.58 2528.48 13500.00 41.70 381.64 1417.86 221.87 424.50 73.40 59.72

41、1520.00 27.80 83.70 175.71 218.38 630.00 96.68 114.81 645.00 34.80 130.36 1626.02 210.53 635.30 64.03 101.35 190.00 28.30 162.64 615.10 22.34 353.00 24.53 12.70 11.79 9.00 58.80 89.08 29.35 407.50 55.54 61.83 112.00 24.05 66.82 208.27 24.79 289.70 42.12 41.10 57.00 17.49 76.80 237.45 21.61 295.80 15

42、.40 18.35 19.00 4.27 40.42 106.53 24.79 178.60 32.31 29.01 45.00 14.26 47.98 104.90 28.67 423.30 61.41 117.83 216.00 30.30 178.88 293.16 27.12 967.70 36.41 70.71 66.00 16.87 119.35 457.96 24.58 445.70 28.32 132.05 139.00 12.91 87.44 303.06 28.23 420.60 35.81 73.86 561.00 22.49 60.91 241.63 24.58 257

43、.40 42.25 44.84 260.00 19.58 46.90 192.65 28.23 350.30 59.07 36.58 177.00 24.57 52.56 170.00 23.97 306.20 38.30 36.22 115.00 13.54 60.63 145.51 210.74 306.50 53.98 97.05 255.00 27.91 69.40 177.86 23.77 515.60 28.74 77.47 1801.00 15.10 70.18 178.57 210.27 175.90 57.88 37.50 65.00 26.97 45.34 100.41 2

44、8.90 117.50 32.54 20.51 41.00 14.47 41.45 56.33 23.37 379.10 32.13 37.81 49.00 13.74 70.37 208.37 26.69 320.30 87.90 66.99 78.00 21.90 123.60 230.11 25.00 551.10 71.75 71.23 74.00 22.30 111.62 209.50 29.58 201.00 39.11 25.26 125.00 10.70 99.35 90.07 23.17 203.80 26.57 24.29 23.00 11.90 41.69 94.34 2

45、3.77 355.20 26.25 32.87 14.00 9.60 59.22 117.87 23.37 240.50 28.04 22.63 261.00 11.70 35.74 92.46 26.05 193.00 40.60 24.88 27.00 14.40 33.53 84.86 23.77 467.10 49.03 34.44 45.00 15.40 60.83 132.86 25.41 364.70 40.34 40.93 79.00 18.70 83.32 175.34 28.23 409.90 44.67 66.92 80.00 36.00 96.85 197.63 25.41 302.50 34.22 27.60 408.00 14.80 68.70 218.24 27.78 114.50 56.38 26.96 36.00 22.40 31.24 75.91 26.47 165.20 73.40 42.73 40.00 19.70 84.13 95.69 .表3-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)34.09 127.00 27.58 23.