全国数学建模竞赛论文参考探索.doc

2013高教社杯全国大学生数学建模竞赛 承 诺 书 我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛 赛区评阅编号（由赛区组委会评阅前进行编号）： 2013高教社杯全国大学生数学建模竞赛 编 号 专 用 页 赛区评阅编号（由赛区组委会评阅前进行编号）： 赛区评阅记录（可供赛区评阅时使用）： 评 阅 人 评 分 备 注 全国统一编号（由赛区组委会送交全国前编号）： 全国评阅编号（由全国组委会评阅前进行编号）：  城市表层土壤重金属污染分析 摘 要 为更好地研究城市地质环境的演变模式，本文通过对城市表层土壤重金属污染情况的分析，得到该研究区域的重金属污染程度、主要污染原因及污染源。 问题一，研究土壤中重金属的空间分布与各功能区的污染程度。首先，建立克里格(Kriging)插值模型，对各区域的地理坐标和土壤重金属浓度进行空间数据插值。然后，利用Suffer10．0专业地理绘图软件中的克里格插值模块，绘出各区域的3维地理分布图和8种主要重金属元素在该城区的空间分布图（见图1-9）。最后，应用Newerow综合污染指数法对各重金属元素的污染情况进行评价，利用Matlab12.0软件，得到了生活区、工业区、山区、主干道路区及公园绿地区的污染比例分布为97.7%、100.0%、86.4%、99.3%、94.3%，重度污染比例分别为40.9%、55.6%、4.5%、56.5%、34.3%。 问题二，研究该城区重金属污染的主要原因。首先，建立基于因子分析的数据分析模型。然后，利用Spss19.0软件中多元统计分析的因子分析法对8种元素进行相关性分析及主成分提取，得到该区域内重金属污染的主要原因是工业污染、交通污染和生活污染。其主要污染元素为Cd、Cr、Cu、Hg、Ni、Pb、Zn。 问题三，研究该城区污染源的位置。首先，利用等浓度线图（见图10）确定重金属污染的传播特征。然后，根据此传播特征及质量守恒定理（即平衡原理），建立一个纯扩散方程模型。最后，采用计算机仿真技术，利用Matlab12.0软件，求解出重金属元素As、Cd、Cr、Cu、Hg、Ni、Pb、Zn的污染源个数分别为4、3、3、2、3、2、2、5个，8种重金属元素的污染源的具体坐标见表。 问题四，更好的研究了该城区地质环境的演变模式。首先，利用已知数据对问题三中的纯扩散方程模型进行检验，发现纯扩散方程模型存在着误差。然后，为更好的研究了该城区地质环境的演变模式，需要收集重金属元素沿各方向的扩散系数，建立偏微分方程模型，提高模型的广泛性。 文章最后，对所建模型进行了优缺点评价，并提出了相应的推广方案。 关键词：重金属污染；Kriging插值；Newerow综合污染指数法；因子分析；扩散方程 1、问题重述 1.1 问题背景 由于近来城市经济快速发展和人口不断增加，人类活动对于城市土壤地质环境的影响日益显著，土壤地质环境已经成为人们越来越关注的焦点和研究的热点[1-2]。近年来，随着工业化与城市化也迅速发展，土壤地质环境受重金属污染已经变成全球性问题[3]。由于含重金属的污染物通过各种途径进入土壤，造成土壤严重污染。土壤重金属污染可影响农作物产量和质量的下降，并可通过食物链危害人类的健康，也可以导致大气和水环境质量的进一步恶化，因此引起世界各国的广泛重视[4]。 重金属，特别是汞、镉、铅、铬等具有显著和生物毒性。它们在水体中不能被微生物降解，而只能发生各种形态相互转化和分散、富集过程(即迁移)[5]。重金属污染的特点是：(1)除被悬浮物带走的外，会因吸附沉淀 作用而富集于排污口附近的底泥中，成为长期的次生污染源；(2)水中各种无机配位体(氯离子、硫酸离子、氢氧离子等)和有机配位体(腐蚀质等)会与其生成络合物或螯合物，导致重金属有更大的水溶解度而使已进入底泥的重金属又可能重新释放出来；(3)重金属的价态不同，其活性与毒性不同。其形态又随pH和氧化还原条件而转化。 1.2 需要解决的问题 按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、、5类区，不同的区域环境受人类活动影响的程度不同（假设）。 现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10 厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值（假设