2013年数学建模分析思路.docVIP

一 问题重述 对某城市城区土壤地质环境进行调查。按照每平方公里1个采样点对表层土（0~10 厘米深度）进行取样、编号，并用GPS记录采样点的位置。分析获得每个样本所含的多种化学元素的浓度数据。按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。 通过数学建模来完成以下任务： (1) 分析8种主要重金属元素在该城区的空间分布及其对该城区内不同区域重金属的污染程度。 (2) 分析数据，说明重金属污染的主要原因。 (3) 分析重金属污染物的传播特征，建立模型分析污染源的位置。 (4) 分析建立模型的优缺点，并说明如何建立模型可以使其更好地研究城市演变。 二 模型假设 1 .外界环境对重金属的传播影响是均匀的， 三 模型的建立及分析求解 问题一 显然本题要求我们对该城区进行横向和纵向的分析， 先进行纵向分析，研究各种元素随海拔的上升其浓度的变化，用matlab软件绘图分析，图如下： 图一 As浓度-海拔 图二 Cd浓度-海拔 图三 Cr浓度-海拔 图四 Cu浓度-海拔 图五 Hg浓度-海拔 图六 Ni浓度-海拔 从图表可以看出，各种元素随海拔的升高，其变化并不是太明显， As Cu Ｈg等八种元素随海拔的升高，其浓度不断降低，不断趋于正常水平，高海拔地区的重金属元素少，而低海拔地区的重金属比较 图七 Pb浓度-海拔 图八 Zn浓度-海拔 多。特别是在海拔100米以下，重金属元素占有比例明显高于正常水平，用数据拟合方法对这些离散量进行拟合，可以发现，随海拔的上升，这些元素的量与海拔对应有二元函数的关系。这些关系式分别如下： 由函数关系式可以推断，随着海拔的上升，这些元素的量又会慢慢的有所增加，这些量在一定程度上不会减少，会不断地累积起来，而且这些元素随海拔的变化大致是一样的，这说明各种元素的变化是均匀的，它们以固定的比例在变动。 纵向分析，利用matlab软件绘制各元素浓度对应坐标的三维图，如下： As浓度分布图 1-1 Cd浓度分布图1-2 Cr浓度分布图 1-3 Cu浓度分布图1-4 假设以横轴正向为东，纵轴正向为南为方向坐标，由图1-1分析，As元素在整个城区的含量分布波动较大，五个地方出现较大的峰值，中西部的元素含量相对较大，由图1-2可以看出，该地区的Cd元素变化比As元素的波动更大，很多地方该元素都比较多，但主要集中在西部，南部，西南部，在东部的小部分地区也有较明显的变化。在图1-3和1-4中，Cr和Cu元素的分布规律大致相同，这两种元素都主要集中在西南部地区，且分布较多，峰值明显，由西向东，元素的分布趋于平缓，没有较大波动，元素分布较为均匀，Hg 元素有三个地方出现峰值，表明其主要分布在该城区的中部，南部和西南部，在其 Hg浓度分布图 1-5 Ni浓度分布图 1-6 Pb浓度分布图 1-7 Zn浓度分布图 1-8 它地区分布较少，而其相互间没有太大差异，Ni元素主要分布在西边和东边，其它地区有较小的波动。由图1-7分析，很明显可以看出，Pb元素至西向东含量不断减少，主要集中于城区西部，其量相对较多。Zn元素在东部和北部分布较少，主要集中在中部，中西部。从整体上看，各种元素在西南部都比较集中，含量比较大。 污染程度分析 对给出的数据在各区域求其平均值，可以在一定程度上表示该地区的重金属含量，表如下； 编号 As (μg/g) Cd (ng/g) Cr (μg/g) Cu (μg/g) Hg (ng/g) Ni (μg/g) Pb (μg/g) Zn (μg/g) 1 6.27 289.96 69.02 49.40 93.04 18.34 69.11 237.01 2 7.25 393.11 53.41 127.54 642.36 19.81 93.04 277.93 3 4.04 152.32 38.96 17.32 40.96 15.45 36.56 73.29 4 5.71 360.01 58.05