2011数学建模A题全国一等奖论文(城市重金属污染及地质演变).pdf

2011 高教社杯全国大学生数学建模竞赛 承 诺 书 我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网 上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。 我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的 资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参 考文献中明确列出。 我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规 则的行为，我们将受到严肃处理。 我们参赛选择的题号是（从A/B/C/D 中选择一项填写）： 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）：中国石油大学（北京） 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名) ： 日期： 年 月 日 赛区评阅编号（由赛区组委会评阅前进行编号）： 2010 高教社杯全国大学生数学建模竞赛 编 号 专 用 页 赛区评阅编号（由赛区组委会评阅前进行编号）： 赛区评阅记录（可供赛区评阅时使用）： 评 阅 人 评 分 备 注 全国统一编号（由赛区组委会送交全国前编号）： 全国评阅编号（由全国组委会评阅前进行编号）： 城市重金属污染及地质环境演变分析 摘要 在问题1 中，我们首先根据离散的300 余个样本的坐标及其相应的重金属浓度，用 matlab 分别拟合绘制8 种重金属元素浓度随坐标位置变化的浓度分布渐变图，用不同颜 色代表不同的污染程度，定性直观地看出 8 种重金属元素在该城区的空间分布。同时， 我们引入污染指数 ，利用地累积指数法对8 种元素在5 块区域的污染程度进行定量计 I 算，并且对污染程度进行分级，不仅定量分析了不同区域的重金属污染情况，而且对第 一部分的拟合图做了检验，两部分吻合很好。 在问题2 中，我们建立主成分分析模型，运用主成分分析法，用较少的新变量 （主 成分）代替原来较多的旧变量 （8 种重金属元素），既减少了变量个数，又保证了新变量 能够覆盖大部分信息。这样便可以比较容易地通过主成分挖掘出各污染元素与区域功能 之间的关系，进而确定重金属污染物的来源。 在问题3 中，我们分析了重金属污染物在地层中的扩散属于“分子扩散”而无“对 流扩散”，符合菲克扩散定律的传播特征。我们根据空间中任一体积元内金属污染物的 物料守恒，结合反映金属元素传播特征的菲克扩散定律，建立金属元素在地层中的扩散 模型，得到浓度关于坐标的拉普拉斯方程。我们采用有限差分法进行数值求解，在问题 一浓度分布图的基础上，利用直观浓度最大点附近的若干已知样本数据作为边界条件， 借助Excel 进行迭代运算，使所得到的数值解与真实情况最大误差为0.001 。再用matlab 筛选出浓度最大点的坐标即为要求的污染源位置。由此分别求得8 种重金属的污染源。 在问题4 中，我们首先分析了研究地质环境转换的具体内容，也就是既要研究不同 区域地质环境状况的“横向对比”，又要研究不同时间地质环境状况的“纵向对比”，对 前三问建立的模型进行了优缺点分析。并且提出要从物质组成、地质结构、动力作用三 个方面补充收集 14 类信息，建立层次分析模型，得到了一套地质环境状况评价指标， 也就是各因素所占权重（权重的变化可以看出演化的过程）及综合评价指标。再将不同 时期的评价体系进行对比，即可得出城市地质环境在的演变模式。 关键词：地污染指数；主成分分析；菲克定律；层次分析