全国大学生数学建模竞赛-城市表层土壤重金属污染分析.doc

PAGE

PAGE1

A题城市表层土壤重金属污染分析（本文2011年全国一第奬）

摘要

本文通过对附件中所给海量数据进行处理分析，找出主要污染原因及污染物传播特性，进而来研究人类活动影响下城市地质环境的变化模式，对地质环境的保护起到积极作用。

对于问题一，先用matlab软件对数据插值拟合，作出各元素在该城区的空间分布图。从图中可直接地观察各元素对不同地点的污染程度，再将附表中所给数据标准化。本文采用目前学术界应用最普遍的单因子浸染指数法和内梅罗综合浸染数法对不同功能区进行逐一分析。得出各个区域重金属污染程度，其中生活区、工业区与主干道路区-重度污染、山区-轻度污染、公园绿地区-中度污染，且它们之间综合污染指数排序为工业区主干道路区生活区公园绿地区山区。

对于问题二，根据问题一中计算出的单因子浸染指数，并作出各功能区不同元素的单因子浸染指数柱状图。通过比较分析发现，某些元素在各功能区单因子指数相近（如As、Ni、Cr），说明同种元素在整个城区的污染原因相同。通过matlab软件分别算出在整个城区和各功能区各元素之间的相关度，可以看出有些元素呈现一定相关性，若相关显著可知其间污染原因可能相同。为了进一步明确不同功能区土壤中重金属的污染来源，再以公园绿地区为例利用主成分分析方法，得出公园绿地区八个变量的全部信息可由四个主成分表示，并分析污染物的来源。根据上述模型进行分析，可知As的主要污染原因为：在某些有色金属的开发和冶炼中，砷化物的广泛利用，煤的燃烧。Ni和Cr相关度较高，主要污染原因有：地球化学成因影响；整个城区的多种人为活动产生的综合污染。除这些元素外，再对其他元素的污染进行分区考虑，根据其它元素之间的相关度可知：生活区Cu、Cd、Pb主要污染原因为生活废弃物；Hg、Zn和其他元素均不相关，说明其主要污染原因是重金属的运移。工业区污染原因主要为工业生产污染。山区轻度污染，主要原因为重金属的运移。主干道路区污染比较复杂，其主要污染原因为交通运输、道路建设、人为活动、金属运移等。公园绿地区处于中度污染，主要污染原因为人类活动和重金属的运移。

对于问题三，综合考虑对流、弥散、扩散吸附和微生物降解等作用,根据质量守恒原理,建立了污染物在饱和及非饱和土壤中运移的对流扩散数学模型:

根据附表数据拟合出扩散系数（=1…10），再利用初始浓度及长时间的稳定态作为初值条件构成微分方程的定解问题，最后采用数值模拟方法预测污染物的迁移过程、迁移范围及浓度分布特征。并由此来分析污染物的传播特征并确定污染源的位置（见表5），结合问题二中的主成分分析的结果，确定最终污染源有7个。

对问题四，问题三模型建立后，发现如果能够收集到城市详细的卫星图、城市历史地质资料、同一位置不同时间的重金属浓度、城市土壤更多的测定值等信息，能够更好的建立污染物迁移速度的演变模式。可以建立通过聚类分析法、人工神经网络的SOFM网络得出城市不同时期的综合评价结果，从而分析出城市地质环境的演变模式。

关键词：环境评价单因子指数内梅罗综合数相关度分析主成分分析对流扩散模型

一、问题重述

随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。

现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其

附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。

现要求你们通过数学建模来完成以下任务：

(1)给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。

(2)通过数据分析，说明重金属污染的主要原因。

(3)分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。

(4)分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？

二、问题分析

本题主要是研究表层土壤重金属污染的问题，主要目的是通过对附表中所给八种主要重金属的采样数据进行分析，从而得到这八种元素的空间分布，并确定各功能区的污染程度；再根据数据分析，先找出各金属间的相关度，确定其是否来自同一污染源，再根据相关分类标准，找出重金属污染的主要原因；然