城市表层土壤重金属污染分析(定稿).doc

城市表层土壤重金属污染分析 摘要：通过采集该区以采样点的位置作为自变量, 采样点处每种重金属的浓度为因变量利用MatlabR2009a进行插值拟合, 给出了8重金属在该城区的空间分布, 5个功能区, 用Excel 2010绘制出了每个功能区中8重金属含量柱状图, 并对结果进行了分析. 针对该城区不同重金属的污染程度, 污染因等问题, 对5个功能区的数据分别采用对应因子分析法进行处理, 结果表明每个功能区均可分为四类, 每类都包含有限采样点和相应的污染重金属, 给出了每个功能区中重金属的种类, 用单因子污染指数评价法和内梅罗污染综合指数评价法定量地了5个功能区中8种重金属的污染程度, 结果表明, 不同功能区单因子污染指数都大于1, 该地区不同功能区土壤重金属污染表现为富集现象, 其中Hg和Cu污染程度最高, 平均含量分别为背景值的7.65和4.35倍. 功能区污染程度的趋势大致为HgCuZnCdPbCrAsNi, 生活区中Cu, Zn, 工业区中的Cd, Cu, Hg, Zn, 主干道路区中的Cu, Hg, Zn公园绿地区中的Hg污染均达到重度污染. 内梅罗污染综合指数评价中以工业区, 主干道区综合指数为高, 工业区的内梅罗污染综合指数最高, 主要原因是重金属Hg污染严重所致. 根据土壤中重金属污染物的传播特征, 对采样点进行, 分为无污染、轻度污染、中度污染和重度污染四类；然后对于重度污染区的样本点按照坐标进行了二次模糊C-均值聚类, 得到空间相近聚类点；利用二次聚类点的标准化后的重金属含量作为权重, 建立了重金属污染物的传播特征模型, 确定出了污染源的位置. 最后, 对所建模型的优缺点作了客观的评价, 关键字：城市表层土壤；重金属；污染分析；因子分析法；模糊C均值聚类分析；Matlab一、问题重述 随着城市经济的快速发展和城市人口的不断增加, 人类活动对城市环境质量的影响日显突出. 对城市土壤地质环境异常的查证, 以及如何应用查证获得的海量数据资料开展城市环境质量评价, 研究人类活动影响下城市地质环境的演变模式, 日益成为人们关注的焦点. 按照功能划分, 城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等, 分别记为1类区、2类区、……、5类区, 不同的区域环境受人类活动影响的程度不同. 现对某城市城区土壤地质环境进行调查. 为此, 将所考察的城区划分为间距1公里左右的网格子区域, 按照每平方公里1个采样点对表层土（0~10 厘米深度）进行取样、编号, 并用GPS记录采样点的位置. 应用专门仪器测试分析, 获得了每个样本所含的多种化学元素的浓度数据. 另一方面, 按照2公里的间距在那些远离人群及工业活动的自然区取样, 将其作为该城区表层土壤中元素的背景值. 附件1列出了采样点的位置、海拔高度及其所属功能区等信息, 附件2列出了8种主要重金属元素在采样点处的浓度, 附件3列出了8种主要重金属元素的背景值. 现通过数学建模来完成以下任务： (1) 给出8种主要重金属元素在该城区的空间分布, 并分析该城区内不同区域重金属的污染程度. (2) 通过数据分析, 说明重金属污染的主要原因. (3) 分析重金属污染物的传播特征, 由此建立模型, 确定污染源的位置. (4) 分析你所建立模型的优缺点, 为更好地研究城市地质环境的演变模式, 还应收集什么信息？有了这些信息, 如何建立模型解决问题？ 二、问题假设与符号说明 根据题意, 可以进行如下假设： 1．不考虑时间对重金属的因素2．3．． 根据题意需要, 我们对用到的符号做以下说明： 是土壤中重金属的污染指数； 是土壤中重金属的实测值； 是土壤中重金属的评价标准是是个变量的平均值是个变量的标准差是标准是相似系数是隶属度矩阵； 是聚类中心. 三、问题分析及方法 城市不同功能区具有不同的土地利用特征, 人类活动所施加的影响也不同. 选取该城区的生活区、工业区、山区、主干道路区及公园绿地5个功能区为研究对象, 共布设319个采样点, 其中生活区布设44个、工业区36个、山区66 个、主干道路区138个、公园绿地35个. 由于该地区各个功能区土壤重金属元素含量的数据特征完全符合因子分析的要求, 因此我们对5个功能区土壤重金属的污染状况均采用因子分析法进行分析. 在这里以Hg、Cd、Pb、As、Cu、Cr、Ni、Zn 八种重金属作因子分析, 这样在解释各指标变化异常时可以着重讨论综合指标因子, 同时为该地区重金属污染成因的解释提供一定的理论依据. 因子分析是多元统计数学中的一种方法, 它根据多个实测变量之间的相互关系, 运用数学变换, 将多个变量转变为少数几个线性不相关的综合指标, 从而简化数据处理, 其目的在于对大量观测数据, 用较少的有代表性的因子来说明众多变