数学建模论文24对A..docVIP

城市表层土壤重金属污染分析 摘要 由于人类活动的增加，城市表层土壤重金属污染问题日益严重。本文通过分析数据，探讨重金属污染的传播特征及影响因素，据此建立偏微分方程数学模型，寻找城区土壤表层重金属污染源的位置。 对问题一：利用Matlab，地理统计软件ArcGis等绘制重金属浓度的空间分布图（见图4-1至图4-10）；然后采用单因子指数法计算出不同区域内各种金属的污染程度累积贡献率值，并绘制直方图(见图4-12) ；运用内梅罗指数法进行污染级别的综合评价，绘制综合评价结果图（见图4-11）。通过以上数据分析，得到几个土壤表层土壤重金属污染的特征（见正文第10页）。 对问题二：建立城区8种重金属的相关系数矩阵，运用SPSS软件进行主因子分析，寻找污染城市土壤的主要重金属元素以及剔除掉污染权重小（权重低于10%）的元素，再结合重金属元素的地理分布，说明各种重金属污染的主要原因（见正文第13页），可以得到人类活动频繁区域是往往是重污染区，所以人类活动是土壤重金属污染的主要原因。 对问题三：利用数据分析发现的传播特征，利用三元泊松方程建立偏微分方程模型，设是关于样本点空间坐标的重金属浓度函数，则模型为 ， 亦可写为，其中表示向量的散度，表示向量的梯度。在这个模型求解中，为提高操作效率，设海拔为定值，则以上模型退化为二元泊松方程，另由于只需考虑静态的污染状况，所以不考虑散度，因此由模型知只需利用向量的梯度求解。比较平面内样本点的梯度信息，得到几个具有极大梯度值的采样点(比较结果见附录表1)。根据不同功能区的重金属污染元素的浓度随海拔区域的的变化数据表，绘制出每种重金属污染质与海拔高度区域的统计直方分布图（见图6-1至6-5），分析了海拔对污染程度的影响。然后，分别在这几个极值点的局部（以3公里为半径），用二维正态分布函数模拟重金属浓度的变化情况，利用Powell直接法（算法流程图见附录9）求最小二乘法的最小值，以此确定二维正态分布函数的参数（Matlab程序见附录6），进而解得此正态分布的极大值点及其坐标，再确定此点对应的海拔值，此点即为污染源的位置。例如，得到Pb在该城区的部分污染源位置坐标如下： 元素名称 编号 横坐标 纵坐标 海拔 所在功能区 Pb 1 14.1730 11.9410 0.014 5 Pb 2 4.7768 4.6134 0.008 1 Pb 3 22.6240 4.8180 0.027 3 最后对该模型进行正确性与可靠性分析，所建立模型是可靠的(见正文第17页）。 对问题四：分析了模型的优缺点。再增加考虑其他信息对重金属污染的影响，例如风向，土壤介质等，建立带有变系数的偏微分方程模型（见正文第19页）。 关键词：土壤表层重金属污染；泊松方程；正态分布；Powell直接法 一、问题的提出 随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。 按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、……、5类区，不同的区域环境受人类活动影响的程度不同。 (1) 用计算机模拟绘制出等值图像来给出8种主要重金属元素在该城区的空间分布，并用直方分布图分析该城区内不同区域重金属的污染程度。 (2) 通过建立不同重金属污染元素间数据关系建立先关系数矩阵来进行数据分析，说明重金属污染的主要原因。 (3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。 (4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？ 二、问题的分析 调查城市城区土壤地质环境调查。所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（ 厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。 附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。现在我们需要根据这些信息解决以下下的问题： 问题一：拟通过附件给出的城区319个样本点的横、纵坐标，海拔高度及所属功能区，通过软件来绘制出城市地貌轮廓；其次再依据样本点的重金属污染元素的浓度分布，绘制等值二维图得到8种主要重金属元素在该城区的空间分布；最后，根据功能区的划分以及样本点浓度信息绘制柱状图来分析该城区内不同区域