2011年数学建模A题__城市表层土壤重金属污染分析.docVIP

1.2 模型建立： 重金属污染评价方法种类较为繁多，从地质环境角度来看，应用于土壤重金属污染评价的方法主要有单因子指数法、地质累计指数法（Muller指数）、内梅罗指数法等。 我们通过查阅资料了解到，地质累积指数[5]不仅反映了元素的自然变化特征, 而且可以判别人为活动对环境影响, 是区分人为活动的重要参数。但是Muller指数只能给出各采样点的某种重金属的污染指数，无法给出各采样点的多种重金属污染综合指数和某种重金属的地区综合指数，即不能全面地反映土壤的污染状况。由此，我们认为不能只用一种评价方法，内梅罗指数[5]不仅考虑到各种影响参数的平均污染状况, 而且特别强调了污染最严重的因子,同时在加权过程中避免了权系数中主观因素的影响, 克服了平均值法各种污染物分担的缺陷。 因此，本小题我们采用综合指数法—地质累积指数与内梅罗指数结合法，来分析该城区内不同区域重金属的污染程度。 目前土壤环境质量评价方法有很多，各有优点和缺点。本文根据我国《土壤监测技术规范》(HJ/T 166-2004)[2]中推荐的内梅罗综合污染指数法进行评价。在计算某个区域某种重金属单项污染指数(分指数)的基础上，再计算该区域多种重金属的综合污染指数。 单项污染指数 （1） （2） 当Pij≤1时，表示土壤未受该因子污染，当Pij 1时，表示土壤受该因子污染。内梅罗综合污染指数反映了各污染物对土壤的作用，同时突出了高浓度污染物对土壤环境质量的影响。根据HJ/T 166-2004,内梅罗综合污染指数的分级标准（见表1），得出各个区域的污染等级。 表1：内梅罗综合污染指数的分级标准 编号 区域名称 综合污染指数 污染程度 1 生活区 重污染 2 工业区 重污染 3 山区 轻度污染 4 交通区 重污染 5 公园绿地区 中度污染 1.3 模型求解： 本文以背景值作为评价标准进行求解，用EXCEL对文中所给数据进行分类，把数据分入1区、2区、3区、4区、5区。然后得出各个区里面主要重金属含量的平均值，可看作各区中主要重金属含量值。如下表： 表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) 1 6.27 289.96 69.02 49.4 93.04 18.34 69.11 237.01 2 7.25 393.11 53.41 127.54 642.36 19.81 93.04 277.93 3 4.04 152.32 38.96 17.32 40.96 15.45 36.56 73.29 4 5.71 360.01 58.05 62.21 446.82 17.62 63.53 242.85 5 6.26 280.54 43.64 30.19 114.99 15.29 60.71 154.24 然后根据公式（1）、（2）结合MATLAB软件算得各区重金属单项污染指数和综合污染指数,如下表(程序见附录)： 表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) 1 1.7417 2.2305 2.2265 3.7424 2.6583 1.4911 2.2294 3.4349 3.1704 2 2.0139 3.0239 1.7229 9.6621 18.3531 1.6106 3.0013 4.028 13.5331 3 1.1222 1.1717 1.2568 1.3121 1.1703 1.2561 1.1794 1.0622 1.2532 4 1.5861 2.7693 1.8726 4.7129 12.7663 1.4325 2.0494 3.5196 9.4264 5 1.7389 2.158 1.4077 2.2871 3.2854 1.2431 1.9584 2.2354 2.7343 再由内梅罗综合污染指数的分级标准得出各区的综合污染等级，如下表： 表4：各区综合污染等级 区域 污染等级 生活区 5 重污染 工业区 5 重污染 山区 3 轻等污染 主干道路区 5 重污染 公园绿地区 4 中等污染 从表中可以看出，该城区内生活区、工业区、主干道路区属于重污染区，公园绿地区属于中等污染区，山区属于轻度污染区。 重金属的污染原因问题 模型建立 用matlab软件对各种