《2011年全国大学生数学建模大赛A题全国一等奖论文_图文》由会员分享，可在线阅读，更多相关《2011年全国大学生数学建模大赛A题全国一等奖论文_图文（35页珍藏版）》请在金锄头文库上搜索。
1、2011 高教社杯全国大学生数学建模竞赛承 诺 竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从 A/B/C/D 中选择一项填写）： 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印
2、并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 日期： 年 月 日赛区评阅编号（由赛区组委会评阅前进行编号）：2011 高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：1城市表层土壤重金属污染分析摘要随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出，城市土壤表层的重金属污染问题日益严重。本文采用克里格插值法绘制污染物空间分布图，内梅罗评价法对污染状况进行评
3、价，相关分析和主成分分析去探究污染原因，建立二维传播模型确定污染源位置。针对问题 1，为了解决污染物空间分布问题，我们首先采用统计分析的方法大致分析了污染物在各个研究区的分布离散度等情况，然后根据克里格插值法用 件绘制了重金属污染物的空间等值线图并对其进行分析。为了解决污染程度问题，我们采用尼罗德评价法，分析出生活区 染最为严重，整体一般清洁，工业区和交通区污染较为严重，山区和公园绿地大体环境清洁。针对问题 2，为分析重金属污染主要原因，我们首先对重金属污染进行 关分析，相关系数表明大部分重金属污染物之间都有较强的正相关关系，又因信息的重叠度较高，我们采用主成分分析法解析污染来源，发现 r、b
4、、几种重金属污染元素很可能来自农药施用，工业“三废”排放，i 主要来自工业区的污染物排放。针对问题 3，为确定污染源位置，首先我们查阅了污染物的传播特征，根据污染物的扩散和迁移，并由此建立污染物二维传播模型 。22随后用 1行多项式拟合确定传播系数，最后用 程求解得:元素 d u i 7 4 2 5 2 5 7针对问题 4，我们先根据模型分析其优缺点，然后考虑运移方程中的土壤密度，分配系数，土壤含水率，重金属污染的积累等条件，对模型进行优化并建立方程。关键词：空间插值 尼梅罗指数评价 主成分分析 相关系数分析 二维传播微分方程21、问题重述随着城市经济的快速发展和城市人口的不断增加，人类活动对
5、城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为 1 类区、2 类区、5 类区，不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距 1 公里左右的网格子区域，按照每平方公里 1 个采样点对表层土（010 厘米深度）进行取样、编号，并用 录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，
6、按照 2 公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。附件 1 列出了采样点的位置、海拔高度及其所属功能区等信息，附件 2 列出了 8种主要重金属元素在采样点处的浓度，附件 3 列出了 8 种主要重金属元素的背景值。现要求你们通过数学建模来完成以下任务：(1) 给出 8 种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。(2) 通过数据分析，说明重金属污染的主要原因。(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集的信息，以及有了
7、这些信息，建立模型解决问题的方法。2、模型假设1、重金属在土壤中的传播方向仅仅是从高浓度向低浓度传播，不考虑其他的传播方向；32、重金属在土壤中传播忽略高程在其中的作用，重金属污染只看做是迁移、扩散；3、不考虑重金属在土壤中的自我降解能力；4、该城市某个区域土壤是均匀的，其中的水的对流速度和扩散系数都是恒定的。5、题目提供的相关数据皆是真实、可信，在没有特殊情况下测得的；3、符号说明符号 意义 符号 意义i 个监测点重金属浓度的测量值 ),(x，y）下重金属污染浓度重金属污染物的质量尼梅罗污染指数 b土壤密度x，y 方向上的弥散系数 y 方向上的对流速度土壤含水率4、题 1 分析问题 1 题目
8、要求我们给出 8 种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。我们对题目给出的附件进行了数据的预处理和正态分布检验，消除了异常值带来的影响和比例效应。根据附件 1、2 给出的样本位置和样本点污染物浓度我们先对数据进行统计分析，然后用 件绘制空间分布等值线图更加直观的观察重金属元素的空间分布。我们采用内梅罗环境污染指数评价的方法结合土壤环境质量标准，做出了污染评价。题 2 分析问题 2 题目要求我们通过数据分析，说明重金属污染的主要原因。我们查阅资料，查出了我国境内 8 种重金属污染物的主要来源。对 8 种重金属污染物先采用 出的结果表明部分重金属之间表现了较强的
9、正相关关系，信息存在一定的重叠，故我们后又采用主成分分析的方法，求出影响4的综合指标，判断不同的重金属污染原因。题 3 分析问题 3 题目要求我们分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。我们查阅资料，分析了重金属污染物在土壤中的传播特征，据此建立重金属污染物的二维传播模型。在建立模型的过程中我们先根据重金属污染的空间分布图确定预备污染源位置，然后采用 1件确定了参数扩散系数和流动系数，最后用件对方程进行求解得出污染源的位置，并对结果进行分析。题 4 分析问题四，要求我们分析我们建立模型的优缺点并且对模型进行优化，我们发现在上述问题进行求解时，我们并没有考虑土壤密度，分配系数
10、，土壤含水率等条件，所以我们应该对系数进行优化，然后根据二维传播方程进行求解。5、据预处理 （1）对异常数据的处理由于测量的值具有不稳定性，且所产生的异常值会对变异函数有很重要的影响，会影响变异函数模型的精度。所以我们将 所表示的范围认为是正常测量值，而超 的测量值则用 进行替代，将小于 的测量值用 替代。其中 示某一重金属元素所有测得浓度的算术平均值，而 s 则表示标准差。所处理得到的数据如附录一所示。 2（2）正态分布检验本文对八种重金属含量的分布进行了 验,由检验结果发现b、d、r、i 这八种重金属元素全部拒绝正态分布零假设，需要进行对数转换，以满足统计学对数据的要求，减少因异常值带来的
11、影响。型分析 （1）主要重金属元素的空间分布问题 1 要求我们给出 8 种主要重金属元素在该城区的空间分布。根据附件 1、2 给出的污染物浓度和采样点坐标，对其先进行数据处理，然后采用统计分析的方法求出各个重金属污染物的离散程度，再采用克里格空间插值法结合 件绘制金属污染物空间等值线图并对求出的结果进行分析。（2）分析不同区域重金属的污染程度尼梅罗污染指数评价在分析不同区域金属污染物浓度的时候，我们应该先筛选出每个区域所对应的站点，然后先用单因子评价对每个站点的每种污染物做评价，然后用尼梅罗指数法对该区域的所有站点进行分析。 2型的建立与求解（1）主要重金属元素的空间分布1、对 8 种重金属元素样本统计分析处理首先对 8 种重金属元素的样本数据进行统计分析处理，采用 件得出各个元素的结果如下表所示：表 1金属 统计分析结果表重金属 计分析报告研究区 样本数 最大值 最小值 均值 标准差 ) 背景值生活区 统计量特征值得知,生活区、工业区、山区、主干道路区、公园绿地区重金属 高值出现在工业区，工业区的变异系数最大,说明工业区 1金属 统计分析结果表重金属 计分析报告研究区 样本数 最大值 最小值 均值 标准差 ) 背景值生活区 统计量特征值得知,生活区、工业区、山区、主干道路区、公园绿地6区重金属 kg
