2011年数学建模A题城市表层土壤重金属污染分析.docVIP

城市表层土壤重金属污染分析 摘要：本文根据单因子指数分析法、内梅罗指数分析法对城市表层土壤重金属污染程度进行评估，利用因子分析法、主成分分析法构建数学模型，并使用Matlab和Spss软件进行数值分析。问题三中研究的KMO值为0.801，较适合做因子分析。卡方检验结果表明，Bartlett球形检验的卡方统计值为1215.773（P0.000），拒绝原假设，相关矩阵不是单位阵。通过以上两项统计指标的检验表明本研究适合进行因子分析。所以此法对开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式提供了科学依据。具体步骤如下： 问题一：结合附件中给出的数据，用matlab软件绘制出三维等高线图，总共有9个图，一个是取样区域的地形图，另外八个是八种重金属元素在该城区的浓度分布图。通过图形我们可以清楚地看到各种重金属元素在不同空间的污染程度，然后利用单因子指数分析法和内梅罗指数分析法，得出不同功能区域的重金属污染程度从高到低依次为：工业区，交通区，生活区，公园绿地区，山区。 问题二：利用Spss软件对所给数据采用因子分析法得出相关系数矩阵、成分矩阵，通过比较各元素之间的相关性，结合主成分分析，判断得出土壤重金属污染的主要原因来自于工业污染，交通污染次之。 问题三：根据重金属的相关性和传播途径，分析重金属的传播特征，同样利用Spss软件对所给数据采用因子分析法、主成分分析法，建立数学模型，并结合各主因子的方差贡献率以及各主因子内部主要指标的载荷系数，构造的综合评价模型 问题四：仔细分析了模型的优缺点，为更好地研究城市地质环境的演变模式，我们又考虑了土壤PH值、降雨量、土壤有机质含量三个因素，并随机产生10组数据，与问题三中综合得分前十名的相关数据结合，同样采用问题二、三中的数据分析方法建立数学模型得出各元素之间的相关系数、得分系数矩阵。结果表明工业区X=2383m，Y=3692m，海拔为7m的样本点总得分最高，这个结果比第三题的估计值更准确，说明了本模型的可靠性。 关键词：单因子指数、内梅罗指数、因子分析法、主成分分析法、土壤重金属污染 1、问题提出 城市土壤作为构成城市环境的重要组成部分，是人类日常生活和生产的物质基础。土壤中含有多种矿物质，就As、Cd、Cr、Cu、Hg、Ni、Pb、Zn等重金属来讲，适当的浓度对动植物体生长有很大的帮助。但是由于人类活动等原因，造成土壤重金属在不同的功能区发生了不同程度的污染。该文通过对城市不同区域的土壤重金属各方面的分析，探寻城市化过程中城市土壤重金属的分布特征和影响因素，对开展城市环境质量评测等都会有很大的帮助。 2、问题重述 随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。 按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、……、5类区，不同的区域环境受人类活动影响的程度不同。 现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10 厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。 附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。 现通过数学建模来完成以下任务： (1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。 (2) 通过数据分析，说明重金属污染的主要原因。 (3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。 (4) 分析所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？ 3、问题分析 3.1 问题1的分析： 问题1要求给出8种重金属在该城市区的空间分布，并分析该城区内不同区域重金属的污染程度。 附件1已给出了调查取样点的地区坐标位置，附件2已给出了8种重金属元素在各个样点处的浓度，附件3是8种重金属在该区域的背景值。 对于以上已知条件，我们用Matlab软件吸收各个样点的位置值，并结合每种重金属在各个样点处的浓度画出8种重金属在该城区内的空间分布图即可。然后采用国内普遍的污染指数评测方法：单因子污染指数法和内梅罗污染指数法 [1]，我们通过计算式可以评价出该城区8种重金属单因子污染程度和不同