2011全国数学建模大赛二等奖A1城市表层土壤重金属污染分析.docVIP

城市表层土壤重金属污染分析 摘要 本文主要对城市表层土壤重金属污染程度及原因进行了分析，并讨论了确定污染源位置的模型。 对于问题（1），首先根据重金属元素浓度的变异系数，将8种金属元素按空间分布分为较均匀的一类（包括As、Cd、Ni、Pb），和不均匀的一类（包括Cr、Cu、Hg、Zn）。并综合海拔等高线和四点样条插值的三维曲面图，反映8种金属浓度在立体空间中的分布情况。接着，建立了污染评价标准，并用加权的思想对现有的内梅罗综合指数法进行了改进，利用改进后的综合指数对不同区域中重金属的污染程度进行量化，得到结果为：生活区、工业区、交通区、生活区的污染程度为重污染；公园绿地区的污染程度为中污染；山区的污染程度为轻污染。 对于问题（2），通过主成分分析确定各功能区土壤的主要污染金属元素，并结合问题（1）中的结论和功能区的分布图，分析出不同功能区重金属污染的主要原因：生活区为交通污染；工业区为工业废物；山区为经大气传播的污染物；交通区为汽车尾气和轮胎摩擦；公园绿地区为绿色植物对金属元素的富集作用。 对于问题（3），根据重金属污染物在土壤中的传播特征，利用了第一扩散定律和高斯定理，建立偏微分方程模型，根据污染物的传播规律对同一重金属浓度进行分析，找出污染源所处的位置，并确定了污染源是点污染源、面污染源还是线污染源。其中以As为例，存在两个位置分别为（18134，10046，41），（12696，3024，27）的点污染源，坐标（4742，7293，9）附近的半径为203m的圆面污染源。 对于问题（4），问题（3）模型中的最大缺点和不足之处是，只考虑了重金属污染物在土壤中的扩散作用。因此，在能够收集对流，弥散，吸附过程的信息的前提下，建立了一个重金属在土壤中传播的二次偏微分方程模型。 本文的最后，对模型的优缺点进行了评价，并提出了改进方向。 关键词： 土壤重金属污染；内梅罗综合指数法；偏微分方程模型 一、问题提出与重述 1.1 问题提出 土壤是一种重要资源，也是人类赖以生存和发展的最根本的物质基础。但是近年来, 随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出，所引起的土壤污染问题也越来越严重。重金属在土壤中的积累，不仅直接影响土壤理化性状、降低土壤生物活性，而且通过食物链数十倍富集，通过多种途径直接或间接地威胁人类健康。 持久性有机污染物( POPs)、重金属、石油类污染物、农药、化肥等,是比较典型的污染。其中,重金属由于其具有隐蔽性和潜伏性,不易降解,具有不可逆性和长期性,造成的后果尤其严重。所以,目前研究人类活动影响下城市地质环境的演变模式，成为环境科学与技术研究的前沿与重点。 1.2 问题重述 城区按照功能划分，一般可分为生活区、工业区、山区、主干道路区及公园绿地区等。不同的区域环境受人类活动影响程度不同。现对某城市城区土壤地质环境进行调查，将所考察的城区划分为间距1公里左右的网格子区域，按每个网格1个采样点对表层土进行取样、编号，并用GPS记录采样点的位置，应用专门仪器获得每个样本所含的多种化学元素的浓度数据。另外，在远离人群及工业活动的自然区取样，作为该城区表层土壤中元素的背景值。 根据所给数据，通过数学建模来完成以下任务： (1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。 (2) 通过数据分析，说明重金属污染的主要原因。 (3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。 (4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？ 二、问题分析 本题要解决的问题是城市表层土壤重金属污染程度及人为因素对重金属含量的影响。 针对问题（1），首先需要对附录中的8种重金属元素浓度离散的样本值进行插值，得到连续的重金属元素浓度值。对于“空间分布”，不能只考虑x、y方向上的分布情况，同时要结合各样本点的海拔，分析在三维空间中的分布。对于不同区域重金属的污染程度，通过资料，我们得知，重金属污染物的主要来源为工业区、交通区及生活区，因此，这三个区域的污染程度应该大于山区和公园绿地区。另外，可以通过比较各区域中重金属的绝对浓度以及相对于背景值的相对浓度来分析该地区的污染程度。 针对问题（2），考虑根据金属元素之间的相关性来对各功能区土壤重金属污染来源作初步判断，成显著相关的元素很有可能来自同一污染源。进而可以用主成分分析进一步确定各功能区土壤的主要污染金属元素，并需要结合问题（1）中的结论和功能区的分布图，来对重金属污染的主要原因进行分析说明。 针对问题（3），由于是对土壤进行了采样，所以我们认为重金属污染物的传播特性是对重金属在土壤中的传播特性的研究。污染物在土壤中的传播主要