2011数学建模国赛A题.docx

城市表层土壤重金属污染分析与预测摘要近几十年，由于工业化和城市化的进程加快，土壤重金属污染日益严重。本文针对重金属污染问题，主要包括四个相关问题，分别为确定重金属元素的空间分布，并分析该城区内不同区域重金属的污染程度，说明重金属污染的主要原因，确定污染源的位置，对模型进行评价与扩展并分析该城市地质环境演变模式。我们在合理的假设下，先绘制出重金属元素空间分布图，再根据单因子指数与内梅罗综合污染指数分析方法对各功能区污染程度分析，说明重金属污染的主要原因。接着，在此基础上建立二维污染物空间扩散模型，采用搜索算法在污染较为严重的区域搜索出适量的样本点，再采用多元非线性回归方法对二维扩散模型进行拟合，求解相应的参数。最后通过对第三问模型优缺点的分析，建立时空扩散模型，更加精确反映城市地质环境的演变情况。针对问题一，首先用软件Surfer绘制出该城区的空间分布图，得出该城区大致的地形地貌，再将纵坐标改为八种重金属元素的含量，采用克里格插值，得到各种重金属在城区的分布图，得出各重金属在不同功能区的污染物空间分布特点。再选用单因子指数和内梅罗综合指数评价污染状况，由数据可以反映出工业区和交通区Hg和Cu严重超标，各功能区污染程度分别为：生活区、山区、交通区及公园绿地区受到了轻度污染，工业区受到了重度污染。针对问题二，对各功能区的污染物含量进行分析，分析推测该城市可能存在规模较大或者数量较多的Hg,Cu,Zn，Pb严重超标的工厂。为进一步揭示表层土壤重金属污染来源，我们对三个污染相对严重的功能区的超标较严重的重金属元素Hg,Cu,Zn,Pb进行相关性分析，运用SPSS软件计算出各相关系数，揭示出不同功能区内土壤表层重金属污染的不同来源，分析可得，交通区污染主要来源于汽车尾气排排放，工业区污染主要来源于工厂的超标排放。针对问题三，根据污染物扩散的特征，建立一维污染物扩散模型：，以此推广成污染物的高斯二维空间扩散模型：。再运用搜索算法在局部极值点附近搜索适量样本点，采用多元非线性拟合方法对二维空间扩散模型进行拟合，得到区域的最大值，即为污染源所在地。针对问题四，通过对上述模型优缺点的分析，进一步搜集相关数据，依据重金属传播所具有的富集性和不易分解性等特征，在空间分布模型的基础上引进时间变量，突出时间对于污染物扩散的影响，建立时空扩散模型，以更好地研究地质环境的演化模式。 关键词:空间分布图 克里格插值 单因子指数 内梅罗综合污染指数 相关性分析 高斯二维空间扩散模型 搜索算法 多元非线性拟合问题重述按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、……、5类区，不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10 厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。现要求你们通过数学建模来完成以下任务：(1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。(2) 通过数据分析，说明重金属污染的主要原因。(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？问题分析对于问题一，我们首先用Surfer软件画出该城区简单的空间分布图，得出该地区的大致地形。然后用同样的方法绘出八种重金属元素在该城区的浓度分布图，最后采用最为常用的土地污染评价方法单因子指数法和内梅罗综合污染指数法对五个不同功能区的污染状况进行评价。对于问题二，由问题一得到各功能区不同种类的重金属含量的单因子指数，并由此分析其重金属污染物的主要类型，推测污染物来源的性质。再选择污染较为严重的几种元素进行相关性分析，分析不同功能区污染来源的差异，并推测各功能区重金属污染的主要原因。对于问题三，通过污染物的传播特点，先建立一维污染物空间分布模型，然后推广成二维空间扩散模型。考虑到污染源的位于污染严重的附近区域，仅在污染严重的区域搜索样本点，以简化算法。最后采用多元非线性拟合方法求解二维高斯分布函数，确定局部极值所对应的坐标轴，即污染源所在地。对于问题四，对模型一和模型三存在的优缺点进行评价，通过对模型优缺点进行分析，了解模型存在的问题，收集更多的