基于系统综合评价的城市表层土壤重金属污染分析-数学建模优秀论文

0基于系统综合评价的城市表层土壤重金属污染分析摘要本文针对城市表层土壤重金属污染问题，首先对各重金属元素进行分析，然后对各种重金属元素的基本数据进行统计分析及无量纲化处理，再对各金属元素进行相关性分析，最后针对各个问题建立模型并求解。针对问题一，我们首先利用 EXCEL 和 SPSS 统计软件对各金属元素的数据进行处理，再利用 Matlab 软件绘制出该城区内 8 种重金属元素的空间分布图最后通过内梅罗污染模型： ，其中 为所有单项污染指数的平均值， 为土壤环境中2/1max2P平 均综 平 均 maxP各单项污染指数中的最大值。得到各区内梅罗综合污染指数，进而求得污染程度为：功能区 生活区 工业区 山区 主干道路区 公园绿地区综P094.22769.824.04625.86.1污染程度土壤、作物受中度污染土壤、作物受严重污染 尚清洁土壤、作物受严重污染 超标针对问题二,我们首先利用 EXCELL 软件画出 8 种元素在各个区内相对含量的柱状图，由图可以明显地看出各个区内各种元素的污染情况，然后再根据重金属元素污染来源及传播特征进行分析，可以得出工业区及生活区重金属的堆积和迁移是造成污染的主要原因，Cu 、Hg 、Zn 主要在工业区和交通区如公路、铁路等交通设施的两侧富集，随时间的推移，工业区、交通区的土壤重金属具有很强的叠加性，受人类活动的影响较大。同时城市人口密度，土地利用率，机动车密度也是造成重金属污染的原因。针对问题三，我们从两个方面考虑建模即以点为传染源和以线为传染源。针对以点为传染源我们建立了两个模型：无约束优化模型 ，得到污22yiximD染源的位置坐标 ；有衰减的扩散过程模型得位置坐标（8500,5500），模型6782,5为：， ukzcyubxahu2222针对以线为传染源我们建立了 模型，并通过线性拟合分析线性污染lceY0源的位置。针对问题四，我们在已有信息的基础上，还应收集不同时间内的样点对应的浓度以及各污染源重金属的产生率。根据高斯浓度模型建立高斯修正模型，得到浓度关于时间和空间的表达式 。uteC0在本题求解过程中，我们所建立的模型与实际紧密联系，有很好的通用性和推广性。但在求点污染源时，我们假设只有一个污染源，而实际上可能有多个点污染源，从而使得误差增大，或者使污染源的位置够不准确。1关键词 内梅罗污染模型 无量纲化 相关性 回归模型 高斯浓度模型 一、问题重述俗话说：“一方水土养一方人” 。城市是人类活动最密集的地区，但在废物处理设施仍不发达的绝大多数地区 ,城市及其周边土壤依然发挥着重要的容纳和净化污染物的功能,在强烈的环境负荷冲击下,土壤的服务功能面临极大的威胁,换言之,土壤的缓冲净化功能将接近极限并有被超过的危险,因而将导致严重土壤污染的产生,而其结果将是长远和危险的。随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。我们将城区分为生活区、工业区、山区、主干道路区及公园绿地区五个部分，分别进行土壤地质环境的调查，对城市环境质量做出评价，希望能有效控制重金属污染物的排放及扩散，制定相关措施保护好我们赖以生存的周边环境，根据题意，本文需要解决的问题有：(1) 给出 8 种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。(2) 通过数据分析，说明重金属污染的主要原因。(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4)为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二、模型假设不考虑元素间的相互作用的影响短期内重金属元素的物理、化学变化及迁移对周围环境影响不大假设附录中所给 8 种重金属元素的背景值真实不考虑历史沉积的重金属的影响三、符号说明第 种元素在第 个采样点的浓度（ ） ；ixij 8,.21i第 种元素浓度的平均值（ ） ；,.i第 种元素在第 个采样点无量纲化后的数值（ ） ；ixij 5,.21j第 个功能区重金属 的单项污染指数（ ） ；iPji ,.j第 个功能区重金属 含量的实测值（ ） ；iC2污染距离积；D污染源位置与已知采样点的距离；h给定采样点的坐标；iyx,四、数据处理4.1 对重元素的分析城市工业“三废”排放，金属采矿和冶炼，家庭燃煤，生活垃圾，汽车尾气排放都增加了城市土壤重金属的负荷。重金属污染环境的主要有汞、铅、铬、锌镉、铜等。其中汞的毒性最大，铬、铅、锌等也有相当大毒性。此外还有砷，砷虽不属于金属.但它的毒性与重金属相似，因此归于重金属一类阐述，称为类金属。目前对我国土壤污染比较普遍的重金属有汞、铬、砷。根据该城区重金属污染的情况，下面对重金属在土壤污染中的来源及传播途径作简要介绍。4.1.1 砷元素该元素毒性很低，水体中含砷污染物主要来自砷和含砷金属矿的开采、冶炼，以及和砷化物为原料的玻璃、颜料、药物、纸张的生产都可产生含砷的废水，造成水体的砷污染。砷及砷化物在水中会在水生物体内累积，但累积程度比其他重金属要低。砷和砷化物，一般可通过水、大气和食物进入人体。4.1.2 镉元素当环境受到镉污染后，镉可在生物体内富集，通过食物链进入人体引起慢性中毒。镉的主要污染源是电镀、采矿、冶炼、染料、电池和化学工业等排放的废水。相当数量的镉通过废气、废水、废渣排入环境，造成污染。镉对土壤的污染主要有气型和水型两种。气型污染主要来自工业废气。镉随废气扩散到工厂周围并自然沉降，蓄积于工厂周围的土壤中，可使土壤中的镉浓度达到 40ppm。水型污染主要是铅锌矿的选矿废水和有关工业（电镀、碱性电池等）废水排入地面水或渗入地下水引起。4.1.3 铬元素对水体污染的铬主要来源于电镀、制革、铝盐生产以及铬矿石开采所排放的废水。是我国水体中一种普遍的污染物。水体中铬污染主要是三价铬和六价铬，它们在水体中的迁移转化有一定的规律性。4.1.4 铜元素铜(Cu) 及其化合物在环境中所造成的污染称为铜污染。主要污染来源是铜锌矿的开采和冶炼、金属加工、机械制造、钢铁生产等。冶炼排放的烟尘是大气铜污染的主要来源。世界铜的年迁移量为:岩石风化 20 万吨,河流输送 11 万吨4.1.5 汞元素汞是在常温下唯一呈液态的金属元素。人类活动造成水体汞污染，主要来自氯碱、塑料、电池、电子等工业排放的废水。由于天然本底情况下汞在大气、土壤和水体中均有分布，所以汞的迁移转化也在陆、水、空之间发生。4.1.6 镍元素3镍污染是由镍及其化合物所引起的环境污染。大 部 分 煤 含 有 微 量 镍 ,通 过 燃 烧 过程 被 释 放 出 来 ,这 是 大 气 中 镍 的 主 要 来 源 。 镍 可 以 在 土 壤 中 富 集 。 土 壤 中 的 镍 主 要来 源 于 岩 石 风 化 ， 大 气 降 尘 ， 灌 溉 用 水 （ 包 括 含 镍 废 水 ） ， 农 田 施 肥 ， 植 物 和 动 物残 体 的 腐 烂 等 。 全世界每年镍的迁移状况是：岩石风化量为 320 000 吨，河流输送量为 19 000 吨,开采量为 560 000 吨，矿物燃料燃烧排放 5 600 吨。4.1.7 铅元素铅 对 环 境 的 污 染 ， 一 是 由 冶 炼 、 制 造 和 使 用 铅 制 品 的 工 矿 企 业 ， 尤 其 是 来 自 有色 金 属 冶 炼 过 程 中 所 排 出 的 含 铅 废 水 、 废 气 和 废 渣 造 成 的 。 二 是 由 汽 车 排 出 的 含 铅废 气 造 成 的 ， 汽 油 中 用 四 乙 基 铅 作 为 抗 爆 剂 （ 每 公 斤 汽 油 用 1 3 克 ） ， 在 汽 油 燃烧 过 程 中 ， 铅 便 随 汽 车 排 出 的 废 气 进 入 大 气 ， 成 为 大 气 的 主 要 铅 污 染 源4.1.8 锌元素锌 在 土 壤 中 富 集 ， 会 使 植 物 体 中 也 富 集 而 导 致 食 用 这 种 植 物 的 人 和 动 物 受 害 。金 属 锌 本 身 无 毒 ， 但 在 焙 烧 硫 化 锌 矿 石 、 熔 锌 、 冶 炼 其 他 含 有 锌 杂 质 的 金 属 的 过 程中 ,以 及 在 铸 铜 过 程 中 产 生 的 大 量 氧 化 锌 等 金 属 烟 尘 ， 对 人 有 直 接 的 危 害 。 其 他 如 橡胶 轮 胎 的 磨 损 以 及 煤 的 燃 烧 也 是 大 气 锌 污 染 的 原 因 。 各种工业废水的排放是引起水体锌污染的主要原因。4.2 对基本数据的分析用 EXCELL 软件和 SPSS 统计软件处理数据如表 1 所示：表 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)平均值 6.27 289.96 69.02 49.40 93.04 18.34 69.11 237.01最大值 11.45 1044.50 744.46 248.85 550.00 32.80 472.48 2893.47最小值 2.34 86.80 18.46 9.73 12.00 8.89 24.43 43.37标准差 2.15 183.68 107.89 47.16 102.90 5.66 72.33 443.64生活区变异系数 0.34 0.63 1.56 0.95 1.11 0.31 1.05 1.87平均值 7.25 393.11 53.41 127.54 642.36 19.81 93.04 277.93最大值 21.87 1092.90 285.58 2528.48 13500.00 41.70 434.80 1626.02最小值 1.61 114.50 15.40 12.70 11.79 4.27 31.24 56.33标准差 1.61 114.50 15.40 12.70 11.79 4.27 31.24 56.33工业区变异系数 0.22 0.29 0.29 0.10 0.02 0.22 0.34 0.20平均值 4.04 152.32 38.96 17.32 40.96 15.45 36.56 73.29最大值 10.99 407.60 173.34 69.06 206.79 74.03 113.84 229.80最小值 1.77 40.00 16.20 2.29 9.64 5.51 19.68 32.86标准差 1.80 78.38 24.59 10.73 27.85 10.43 17.73 30.94山区变异系数 0.44 0.51 0.63 0.62 0.68 0.67 0.49 0.42平均值 5.71 360.01 58.05 62.21 446.82 17.62 63.53 242.85最大值 30.13 1619.80 920.84 1364.85 16000.00 142.50 181.48 3760.82最小值 1.61 50.10 15.32 12.34 8.57 6.19 22.01 40.92标准差 3.24 243.39 81.61 120.22 2180.27 11.79 32.53 384.78交通区变异系数 0.57 0.68 1.41 1.93 4.88 0.67 0.51 1.58公 平均值 6.26 280.54 43.64 30.19 114.99 15.29 60.71 154.244最大值 11.68 1024.90 96.28 143.31 1339.29 29.10 227.40 1389.39最小值 2.77 97.20 16.31 9.04 10.00 7.60 26.89 37.14标准差 2.0