本博客下载链接包含修改的word版本, 可免费下载阅览学习, 也可作为数学建模相关课程作业修改上交:
链接：/s/1pFe64VPDomPDtIQ7gTLV-A?pwd=qtnc 提取码：qtnc
A题 城市表层土壤重金属污染分析
摘要
本文通过对附件中所给海量数据进行处理分析，找出主要污染原因及污染物传播特性，进而来研究人类活动影响下城市地质环境的变化模式，对地质环境的保护起到积极作用。
对于问题一，先用matlab软件对数据插值拟合，作出各元素在该城区的空间分布图。从图中可直接地观察各元素对不同地点的污染程度，再将附表中所给数据标准化。本文采用目前学术界应用最普遍的单因子浸染指数法和内梅罗综合浸染数法对不同功能区进行逐一分析。得出各个区域重金属污染程度，其中生活区、工业区与主干道路区-重度污染、山区-轻度污染、公园绿地区-中度污染，且它们之间综合污染指数排序为工业区>主干道路区>生活区>公园绿地区>山区。
对于问题二，根据问题一中计算出的单因子浸染指数，并作出各功能区不同元素的单因子浸染指数柱状图。通过比较分析发现，某些元素在各功能区单因子指数相近（如As、Ni、Cr），说明同种元素在整个城区的污染原因相同。通过matlab软件分别算出在整个城区和各功能区各元素之间的相关度，可以看出有些元素呈现一定相关性，若相关显著可知其间污染原因可能相同。为了进一步明确不同功能区土壤中重金属的污染来源，再以公园绿地区为例利用主成分分析方法，得出公园绿地区八个变量的全部信息可由四个主成分表示，并分析污染物的来源。根据上述模型进行分析，可知As的主要污染原因为：在某些有色金属的开发和冶炼中，砷化物的广泛利用，煤的燃烧。Ni和Cr相关度较高，主要污染原因有：地球化学成因影响；整个城区的多种人为活动产生的综合污染。除这些元素外，再对其他元素的污染进行分区考虑，根据其它元素之间的相关度可知：生活区Cu、Cd、Pb主要污染原因为生活废弃物；Hg、Zn和其他元素均不相关，说明其主要污染原因是重金属的运移。工业区污染原因主要为工业生产污染。山区轻度污染，主要原因为重金属的运移。主干道路区污染比较复杂，其主要污染原因为交通运输、道路建设、人为活动、金属运移等。公园绿地区处于中度污染，主要污染原因为人类活动和重金属的运移。
对于问题三，综合考虑对流、弥散、扩散吸附和微生物降解等作用, 根据质量守恒原理, 建立了污染物在饱和及非饱和土壤中运移的对流扩散数学模型:
根据附表数据拟合出扩散系，再利用初始浓度及长时间的稳定态作为初值条件构成微分方程的定解问题，最后采用数值模拟方法预测污染物的迁移过程、迁移范围及浓度分布特征。并由此来分析污染物的传播特征并确定污染源的位置（见表5），结合问题二中的主成分分析的结果，确定最终污染源有7个。
对问题四，问题三模型建立后，发现如果能够收集到城市详细的卫星图、城市历史地质资料、同一位置不同时间的重金属浓度、城市土壤更多的测定值等信息，能够更好的建立污染物迁移速度的演变模式。可以建立通过聚类分析法、人工神经网络的SOFM网络得出城市不同时期的综合评价结果，从而分析出城市地质环境的演变模式。
关键词：环境评价 单因子指数 内梅罗综合数 相关度分析 主成分分析 对流扩散模型
