2011高教社杯全国大学生数学建模竞赛
承诺书
我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.
我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）：A 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）：_______________ 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)：
日期：2011年9月12日
赛区评阅编号（由赛区组委会评阅前进行编号）：
2011高教社杯全国大学生数学建模竞赛
编号专用页
评 阅 人 评 分 备 注 赛区评阅编号（由赛区组委会评阅前进行编号）：
赛区评阅记录（可供赛区评阅时使用）：
全国统一编号（由赛区组委会送交全国前编号）：
全国评阅编号（由全国组委会评阅前进行编号）：
基于GS技术的城市土壤重金属污染分析
摘要
本文我们主要解决的问题是如何通过污染物在土壤中的传播特征建立模型求解污染源。在模型建立的过程当中，我们主要应用了因子分析法，变异函数模型，kringing插值方法，运用GS软件绘图，SPSS软件处理数据。
问题一，运用GS软件对319个离散数据做出该城市重金属浓度分布图，利用内梅罗综合污染指数法求出各区域的样本重金属污染率，得到综合污染程度。
问题二，采用因子分析法，分别对8种重金属污染物的浓度指标进行了因子分析，运用spss统计软件处理数据，将这8项指标归结为5个公共因子，在此基础上根据不同区域的因子得分对各区域环境污染状况进行了总体比较和评价。
问题三，使用变异函数模型，找到理论变异函数，通过kringing插值，进行交叉检验得到浓度的估计值，此时产生一个标准误差。通过标准误差的范围，确定特异值。通过GS软件画出特异值的空间分布图，从而找到特异值点，此时的特异值点并非污染源。我们再根据污染物的传播特征，以及特异值周围的地势情况，综合确定污染源。 问题四，基于已有的关于土壤污染物的传播特征，我们收集了土壤的PH值，城市的规划现况，地层分布，以及土壤类型等相关信息，在已知条件下，假设未知参量，建立数学模型，定量分析这些因素对污染物在土壤中的浓度变化的影响。
关键字：GS技术 内梅罗综合污染指数法 因子分析法 污染物的传播特征 变异
函数模型 kringing插值方法
1
1、问题的提出
1.1问题背景
随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。 1.2已知条件
某城区319个采样点的位置、海拔高度及其所属功能区等信息，以及8种主要重金属元素在采样点处的浓度，主要重金属元素的背景值。 1.3目标任务
问题一： 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。
问题二：通过数据分析，说明重金属污染的主要原因。
问题三：分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。 问题四：分析所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？
2、问题分析
2.1问题一分析
重金属元素在该区域的空间分布与污染程度紧密相关。空间分布，即要在图上展现出各元素在空间上的浓度分布。可采用GS软件将已知数据拟合画出二维图形。要分析各区域的污染程度，这是一个综合的指标，所以要结合8种污染物共同分析。基于已有的国际通用污染指数测定法——内梅罗综合污染指数法，我们应该分析不同区域重金属的污染程度需处理样本点重金属浓度数据得到污染样本点占该地区总样本点的百分比，结合分布图得出各区域的污染程度结论。 2.2问题二分析
此问题要求我们对已知数据进行数据分析，由于此问题涉及八种元素，如果单独讨论，选择标准复杂，因此要采用因子分析法进行数据统计和分析，能提炼精简出关联度最强、最具代表性的选择标准，由此构建具有一定合理性和可操作性的数字资源选择标准体系。得到标准体系后，我们就需要对数据进一步量化处理，如计算方差，标准差，残差，协方差，公因子方差等。将这些数据经过分析，得到重金属污染的主要原因。 2.3问题三分析
该问题的重点是分析污染物的传播特征，包括传播途径，方式等，再根据他的传播特征建立具体的数学模型，从而定量的来确定污染源。基于已有的土壤污染物传播特征，查阅相关资料我们了解到，污染物的传播途径大概有4类，但是他们分别是基于大气物理传播，植物根叶吸收,以及溶质的水溶程度理论进行的传播。这里我们就不对污染源的确定进行定量的分析，只做定性的说明。又由于我们只知道该城区的海拔高度，而对于城区的城市建筑，河流，气候等情况没有具体的参数，如果自行假设则会增加模型的复杂度，因此，我们就只从地势坡度上来寻找污染源。根据污染物离散点的浓度分布，
2
基于已有的变异函数模型，通过kringing插值法进行交叉检验，GS绘制图形，得到特异值的空间分布图，再定性和定量相结合确定污染源。 2.4问题四分析
首先是根据建立的模型分别其优缺点。其次是查阅相关的资料，基于已有的土壤污染物的传播特征分析，我们在此分析的基础上收集关于土壤PH值，土壤厚度等信息。有了这些信息，我们可以将模型三中定量说明的部分用具体的数据对其定量化，这样就可以得到确定的污染源。
3、模型假设
1）不考虑地质灾害对重金属元素空间分布的影响； 2）不考虑各重金属污染物浓度测量误差；
3）假设不同污染中心向2四周扩散除相互重叠时，污染程度不受重叠影响； 4）污染源的重金属浓度不在增加；
5）取样点的数据较好的反映了该地区的污染物浓度；
4、符号假设
符号 P综 P平均 P最大 ， ， …, p ， ，…, A 意义 综合污染指数(综合反映各污染物对区域土壤的不同作用) 所有单项污染指数平均值 土壤环境中各单项污染指数中的最大值 区域重金属i的单项污染指数 重金属i含量的实测值 重金属i含量的评价起始值 原有变量，是均值为零、标准差为1的标准化变量 原有变量个数 因子变量 因子载荷矩阵 为 在坐标轴 上得投影 特殊因子 显著性水平 单位 mg／kg mg／kg mg／kg mg／kg mg／kg mg／kg / 个 / / / / / 3
2011高教社杯全国大学生数学建模竞赛A题
