图1：调查地区的地形图
图2：功能区散点图
4
图3：砷和镉在该城区的空间分布图
图4：铬和铜在该城区的空间分布图
图5：汞和镍在该城区的空间分布图
5
图6：铅和锌在该城区的空间分布图
说明：图1的Z轴为海拔高度，X、Y轴为地理坐标值（单位：m）。
图2 的X、Y轴为地理坐标值（单位：m）。
图3-图6的Z轴为重金属元素的浓度（单位：μg/g），X、Y轴为地理坐标值（单位：m）。
1.2 模型建立：
土壤环境质量单项污染指数主要用来评价某一污染物的污染程度，指数小污染轻，指数大污染则重。但区域内土壤环境质量作为一个整体和外区域进行比较时除用单项污染指数外，还常用综合污染指数。综合污染指数可以综合判断某土壤多种污染物的联合污染效应。
目前土壤环境质量评价方法有很多，各有优点和缺点。本文根据我国《土壤监测技术规范》(HJ/T 166-2004)[2]中推荐的内梅罗综合污染指数法进行评价。在计算某个区域某种重金属单项污染指数(分指数)的基础上，再计算该区域多种重金属的综合污染指数。单项污染指数和综合污染指数的计算公式如下：
Pij=Cj/Sj （1）
PN=(P2j,ave+P2j,max)/2 （2）
当Pij≤1时，表示土壤未受该因子污染，当Pij >1时，表示土壤受该因子污染。内梅罗综合污染指数反映了各污染物对土壤的作用，同时突出了高浓度污染物对土壤环境质量的影响。根据HJ/T 166-2004,内梅罗综合污染指数的分级标准（见表1），得出各个区域的污染等级。
表1：内梅罗综合污染指数的分级标准 等级 1 2 3 4 5 内梅罗污染指数 PN≤0.7 0.7＜PN≤1.0 1.0＜PN≤2.0 2.0＜PN≤3.0 PN＞3.0 6
污染等级 清洁（安全） 尚清洁（警戒线） 轻度污染 中度污染 重污染
1.3 模型求解：
本文以背景值作为评价标准进行求解，用EXCEL对文中所给数据进行分类，把数据分入1类区、2类区、3类区、4类区、5类区。然后得出各个区里面主要重金属含量的平均值，可看作各区中主要重金属含量值。如下表：
表2：各区重金属含量的平均值 As Cd Cr Cu Hg Ni Pb Zn 区域 (μg/g) (ng/g) (μg/g) (μg/g) (ng/g) (μg/g) (μg/g) (μg/g) 1 2 3 4 5 6.27 289.96 69.02 49.4 93.04 18.34 69.11 237.01 93.04 277.93 36.56 73.29 7.25 393.11 53.41 127.54 642.36 19.81 4.04 152.32 38.96 5.71 360.01 58.05 6.26 280.54 43.64 17.32 40.96 15.45 62.21 446.82 17.62 30.19 114.99 15.29 63.53 242.85 60.71 154.24 然后根据公式（1）、（2）结合MATLAB软件算得各区重金属单项污染指数和综合污
染指数，如下表：
表3：各区重金属单项污染指数和综合污染指数 区域 1 2 3 4 5 As Cd Cr 单 项 污 染 指 数 Cu Hg 2.6583 1.1703 3.2854 Ni Pb Zn 综合污染指数 3.1704 13.5331 1.2532 9.4264 2.7343 1.4911 2.2294 3.4349 1.2561 1.1794 1.0622 1.2431 1.9584 2.2354 1.7417 2.2305 2.2265 3.742