基于模糊综合评价的土壤重金属污染程度分析毕业论文.pdf

摘要：随着城市经济的快速发展和城市人口的不断增加，人类活动对城市 环境质量的影响日显突出。本文利用2011年全国大学生数学建模竞赛A题 所给数据，借助MATLAB软件绘制了采样点的位置分布图。通过MAP GI S软 件的空间分析功能分别绘制了各重金属的空间分布等值线和表面图。通过 对污染指数、变异系数等的计算，对研究区的重金属污染进行了定量的分 析；建立因子分析模型研究了采样区的重金属污染的主要原因。在以上分 析的基础上，为对该城区不同区域的重金属污染程度进行综合评价，将污 染等级分为清洁、尚清洁、轻污染、中污染以及重污染五级，利用熠值法 确定了八种污染物的权重，然后利用模糊综合评价法对不同区域的的重金 属污染程度进行了评价，通过实证分析表明该方法是合理、可行的。最后，针对以上的分析结果，提出了应对城市土壤重金属污染的对策。关键词：城市土壤；重金属污染；数学软件；因子分析；模糊综合评价Abstrac t:Ac c ording to the development of ec onomy and the c ities inc reased,human ac tivities have a noteworthy influenc e in the environment of urban.Based on the dates of heavy metal elements in soils from CUMCM 2011.The distribution map of the loc ation of sampling points was drawn with MATLAB;the spatial distribution of the c ontour and surfac e maps of various heavy metals were drawn with the spac e analysis func tions of MAP GI S software.By the c alc ulation of the pollution index,c oeffic ient of variation and so on,the 目录1绪论.31.1 研究背景.31.2 研究意义.41.3 国内外研究现状.61.4 本文解决的主要问题.92研究区土壤重金属污染概况.102.1 研究分区与采样点的收集.102.2 采样区重金属元素空间分布.122.3 采样区重金属污染定量分析.182.4 采样区重金属污染传播特征分析.23fac toranalysis;fuzzyc omprehensiveevaluationKeywords:urbansoils;heavymetalpollution;mathematic alsoftware;rational.urbansoil.Theresultsshowthatthemethodmoresc ientificandisusedtoevaluatethec harac teristic sofheavymetalpollutionintheaboveanalysis,themethodofentropyweighingfuzzymathematic sforheavymetalc ontaminationofthesamplingarea.Onthebasisofthefac toranalysismodelwasestablishedtostudythemainreasonstudied;Andthenquantitativeanalysisofheavymetalpollutionoftheareawas.2附录.53参考文献.51谢辞.494控制城市土壤重金属污染的对策.473.4总结.473.3实证分析.413.2模糊综合评价法的步骤.393.1模糊综合评价法的概述.383土壤重金属污染程度的模糊综合评价.382.5采样区重金属污染主要原因分析.25.1绪论1.1 研究背景人类工业文明的发展促进了社会经济得以迅速发展，但工业的发展、城市的扩张、人口的增长已使环境质量日趋恶化，并对农业生产、工业布 局、城市发展、食品安全等带来了一系列问题。土壤作为沉积物的一种存 在形式，是人类生产活动和生存的基础物质资源，是不可再生的资源。土 壤的污染己成为制约区域社会、经济与环境协调发展的瓶颈问题，其中土 壤重金属的污染累积性、生理毒害性、治理难度高、周期性长等问题更为 突出。土壤污染已经成为现代环境污染不可或缺的组成部分。人类在生产 活动中，从自然界取得资源与能量，经过加工、调配、消费最终再以“三 废”的形式直接向土壤或者通过水、大气和生物向土壤中排放和转化。当 输入的污染物质数量超过土壤的容量和自净能力时，必然引起土壤恶化，发生土壤污染。而污染了的土壤又向环境输出污染物质，便引起大气、水 体、生物的进一步污染，从而使环境状况发生变化，环境质量下降，造成 环境污染。由于土壤的组成成分、结构、功能、特性以及土壤在环境生态 系统中的特殊地位和作用，使得土壤污染及不同于大气污染，也不同于水 体污染，而且比他们要复杂的多。土壤是植物，特别是作物的生活环境，作为人类主要食物来源的粮食、蔬菜、家畜、家禽等农副产品，都直接或 间接来自土壤。污染物在土壤中的富集必然引起食物污染进而危害人体健 3.1.2 研究意义重金属污染是人类活动将重金属及其化合物施加到环境中，使环境中 的重金属含量超过一定的标准，导致生态环境日益恶化，最终危及生物生 存和人类健康的现象。重金属对生物具有显著的毒害效应，重金属的浓度 增加可通过减慢光合作用速度，引起禾苗缺水，或由于抑制了有机养分的 矿化而使土壤中氮和磷的供应降低等方式对植物的生长产生抑制作用，重 金属污染对人体和动物的危害有很大的隐蔽性，它可以通过直接的体表接 触，也可以通过生物吸收后以食物链方式进入人体，最终在人体内积累。重金属含量达到一定的量就会对人体造成危害，而且其引起的疾病很难被 4加关注。这些都使污染土壤修复工作迫在眉睫。绿色食品、有机食品、无污染食品的需求愈来愈多，对居住环境问题亦更土壤修复提出更加明确的要求及任务。随着我国居民生活质量的提高，对境的重大威胁。2007年，国务院印发了国家环境保护“十一五”规划，对特别是经济发达地区的重金属污染土壤的问题尤其严重，日益成为当地环大；其中，重金属超标被普遍认为是造成我国土壤污染的主要原因之一，并存、无机有机复合污染的局面；土壤污染途径多，原因复杂，控制难度准体系也未形成；与此同时，我国土壤污染的类型多样，呈现新老污染物导致防治措施缺乏针对性；防治土壤污染的法律还存在空白，土壤环境标点。我国由于土壤污染防治基础薄弱，土壤污染的面积、分布和程度不清,康。土壤污染已成为环境保护研究的核心，污染防治的措施成为研究的重.5一定时间的迁移，对厂区周围的土壤环境也产生了一定的影响。城市土壤的土地受到了严重污染，尤以重金属土壤污染为主。土地中的重金属经过的铅污染。城市建设初期建立的化学、工业企业经过多年发展，企业厂区数开放历史较短、客流量小且相对偏僻的公园而言，表土一般都未见明显史悠久，客流量大且距离市中心较近的公园土壤铅含量明显偏高；对大多程度上影响重金属污染状况。对北京城市公园土壤的铅污染研究发现，历状况、人类活动强度、污染累计时间的长短和距离污染源的远近，在不同人类活动是影响城市土壤污染程度的一个重要因素。不同的土地利用累积，造成植被受害，甚至寸草不生，并会间接引起人畜慢性中毒。硝酸盐和有机氯农药等的污水往往在土壤、植被以至地下水中形成残留和乙醛等有机物的污水极易引起急性中毒；含有无机物如重金属、氟化物、污水所含成分复杂，污水性质不同，对土壤危害程度也不同，如含有三氯为严重，在我国工业化进程较快的城市，土壤的铅含量都非常高。另外，标排污等，使重金属大量沉积于土壤中，其中以铅、锌等金属元素污染最污，使城市土壤化学性质发生重要变化。烟尘、汽车尾气的排放、工业超一系列有关城市土壤污染的研究。城市工业化的发展及与之相伴的工业排1988年成立了城市土壤工作组；美国在上世纪90年代对纽约等城市开始了境质量日益恶化。世界各国对此问题开始予以高度重视，德国土壤学会在质量和人类健康。然而，随着工业的发展和城市化进程加快，城市土壤环城市土壤承载着一定的生态、环境和经济功能，关系到城市生态环境大威胁5,6J检查出来，但是一旦被发现就很难或无法治愈，给人类的生存健康构成巨.1.3 国内外研究现状在我国城市化进程不断加快的形势下，各种各样人类活动引起的城市土 壤重金属积累问题受到关注，我国城市土壤研究在20世纪90年代开始有 零星工作，开展了一些关于郊区土壤重金属污染的初步研究，研究城市土 壤重金属的主要工作在2000年后，现在已有20多个城市进行了土壤重金 属污染状况的调查研究，香港和南京的相关工作取得了一些有意义的结果,2006年国家环保局主持开始全国土壤污染状况调查工作。总体而言，这方 面的研究正处于上升期，今后将更为人们所关注，研究也会更加深化。关 6含量分布、化学形态，进行适当合理的评价具有重要的意义。生态学学者们研究关注的焦点。因此，研究城市土壤的重金属污染的源、响。城市土壤重金属的背景值测定，评价标准以及风险预测，是当前城市生变化，同时对植物的生长发育，体貌特征，生命周期等许多因素产生影至生命划。土壤中的重金属污染不仅会导致植物体内的重金属元素含量产的积累，并通过大气、水体或食物链而直接或间接地威胁着人类的健康甚样的人类活动，将大量的重金属带人城市土壤中，造成这些元素在土壤中断加快，城市生态环境质量日益受到广泛关注叼。在城市环境中，各种各增长的人口对食物和纤维的需要，无暇顾及城市土壤。由于城市化进程不长期以来，土壤研究主要集中在农地和林地土壤上，以确保满足世界不断应者，是城市污染物的汇和源，它关系到城市生态环境质量和人类健康。是城市生态系统的重要组成成分，是城市园林植物生长的介质和养分的供.i i i式中，p为污染物的污染指数；X爪为i污染物的实测值；S为1污染物 i i i的评价标准。尼梅罗综合污染指数的计算公式：P+P 2)/2上综合 L 平均 最大 式中，乙为监测点(评价区)的综合污染指数；P为监测点所有污染 综合 最大物单项污染指数中的最大值；P为监测点所有污染物单项污染指数的平均 平均值。(2)地质累积指数评价地质累积指数(Geoac c umulation index,I geo)是 Muller 于 1969 年提 出定量评价重金属元素污染的参数，是20世纪60年代晚期在欧洲发展起7单项重金属污染指数计算公式：P=x/s综合污染评价，最后评价显示为污染级别。的总体质量。利用该指数评价土壤环境质量分为单因子污染指数评价和作用，能给较严重的单项污染物以较大的权值，能较全面地反映土壤环境时兼顾了单因子污染指数平均值和最高值，可以突出污染较重的污染物的尼梅罗(N emerow)指数评价方法(综合污染指数法)。尼梅罗指数评价同(1)尼梅罗指数评价估。标准与评价项目直接的数学关系，来对土壤的重金属污染状况进行分级评家土壤环境质量标准(GB156181995)中的二级标准值和背景值，建立方法。每种方法的核心思想都是利用一定的标准，通常采用的标准是国主要有尼梅罗指数、地质累积指数、污染负荷指数、潜在生态污染指数等于土壤重金属污染的评价问题，出现了各种各样的环境评价技术方法3，.8r单个重金属的潜在生态危害指数简称屈，其计算公式为：nC表示计算所用的参比值。的实测值；f表式中，c表示某一重金属的污染指数；O表示表层沉积物重金属浓度nfCiCi正Cif单个重金属的污染指数简称C,其计算公式为：的敏感性，对重金属敏感性大的水体R/相对较高。属金属的毒性水平，毒性大的重金属在计算处时所占的比重大;水体对重金重受以下几个因素的影响：表层沉积物中重金属的浓度；重金属的种类；大小及生态危害的方法，即潜在生态危害指数法。潜在生态危害指数（口）值的污染瑞典学者Hakanson于1980年建立了一套应用沉积学原理评价重金属（3）潜在生态污染指数（.）法量分数的自然波动，一般将系数K设定为1.5。nn式中，分别为样品中元素质量分数和该元素的背景值，为消除元素质geo2LnI=logC1（KxB）1中得到了比较广泛的应用口31来的广泛研究沉积物中重金属污染程度的定量指标，在目前地球化学评价.Cii=l i=l n以上评价方法在重金属污染问题中均得到了广泛的应用，但在某些区 域因为评价方法的不一致，存在着污染程度或污染级别的差异性，各自评 价方法都是基于数学模型的结果而对污染程度进行定级，其评价方法与评 价结果各有优劣。1.4 本文解决的主要问题随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环 境质量的影响日显突出。本文利用获得的海量数据资料开展城市环境质量 评价，归纳起来，本文的主要研究内容如下：(1)利用MATLAB软件的绘图功能给出采样点在研究区的位置分布图；(2)利用MAP GI S软件给出研究区重金属污染的空间分布图，分析研 究区土壤重金属的空间分布；(3)分析研究区土壤重金属污染的分布特征和空间变异；(4)通过数据分析，说明重金属污染的主要原因；(5)在研究分析以往评价方法的优缺点的基础上，建立模糊综合评价 模型对研究区土壤重金属污染程度进行综合评价；9RI=Ei，WRI=ZTiJ沉积物中多种重金属潜在生态危害指数简称知，其计算公式为：f式中，人为毒性系数；O某一重金属的污染指数。rrfEi=TiC.2研究区土壤重金属污染概况2.1 研究分区与采样点的收集在研究城市土壤重金属的空间分布时，我们根据国家的标准按照功能 划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区 等，分别记为1类区、2类区、5类区，不同的区域环境受人类活动影 响的程度不同。现对城区土壤地质环境进行调查Ml为此，将所考察的城 区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层 土（010厘米深度）进行取样、编号，并用GP S记录采样点的位置（图1、10（6）针对以上的分析结果，给出控制城市土壤重金属污染的对策。.图1采样点位置分布图图2采样点位置对应的等高线分布图表1研究区主要重金属元素的背景值元素 平均值标准偏 范围差it区取样，其作为该城区表层土壤中元素的背景值（表1）。度数据14另一方面，按照2公里的间距在那些远离人群及工业活动的自然图2）。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓.2.2 采样区重金属元素空间分布为了更加清晰的描绘8种重金属的空间分布，我们应该采用表面绘图和 等值线绘图，从各个角度来表述重金属的分布情况，能更直观的看出其分 布特点。地理信息系统在地球化学领域的应用已经日趋广泛，在制作地球化学异 常图方面已经日渐成熟。地球化学图的制作正是运用了 GI S强大的空间分 析功能，制作地球化学异常图就是运用在MAP GI S的空间分析功能。对于元 素地球化学异常图的制作方法如下：首先,将数据进行整理，整理成为GI S软件的坐标数据，即X,Y。将 测试出来的样品的目标元素含量作为Z值，并将X,Y,Z值之间的空格用 逗号替换，此时应该注意将光标移至文本的最后，以保证成图的准确性。把整理好的带有X,Y,Z值的文本文件输入导入到MAP GI S中，过程是12Zn（Mg/g）69144197P b（Mg/g）3161943N i（Mg/g）12.33.84.719.9Hg（ng/g）3581951Cu（Mg/g）13.23.66.（T20.4Cr（Mg/g）3191349Cd（ng/g）130307 0190As（Mg/g）3.60.91.85.4.13采用以上方法制作的8种主要重金属元素的空间分布如图3-10所示。围修改为“原始坐标范围”。致，不能套合上，因为生成的等值线图的坐标是系统默认的，所以要将范原始数据范围的话，那么就会使生成的等值线图的坐标与边界的坐标不一据中的坐标，这里强调的是要选择原始数据范围的原因在于如果没有选择网格化”，并且使用的范围是“原始数据范围”，即就是按照原来的文本数需要说明的一点是在做等值线时用的网格化方法是“K ring泛克立格法才是关于元素含量的图。程”数据的选择是点属性数据中的元素含量，只有这样做出来的等值线图面的DTM分析做等值线图，需要注意的是在做等值线图的时候，里面的“高生成点图元的目的是用于做等值线图，再在MAP GI S的空间分析模块下保留的小数与刚才文本数据本身的位数相同，即两位。分隔符，这里指定的分隔符是逗号，属性名称所在的行是X,Y,Z,字段所3）选择数据的读取方式按照行读取，从X-Y的顺序，与此同时要指定进行投影变换。件在这里打开，并且按照要求选择“用户投影参数”和“目的投影参数”2）选用主菜单“投影变换一用户文件投影变换”，将整理好的文本文1）打开整理好的文本文件，检查是否数据准确，确保无误。如下:.27-26-25-24:23-22-21 L20_ 19 L 18 _ 17 _ 16 _ 15 _ 14 _ 13 _ 12_ 11 _ 10 _9-8 _ 7-6-5 _4 _ 315001G10001102 0 5000 10000(“)5000 20000 25000图3 As的空间分布(g/g)_1350 _1300 _1250 _1200_1150 _1100 _1050 _1000 _950_900 _850 _800,750 _700,650 _600 _550 _500450 _400 _350 _300250 _200.150.10050图4 Cd的空间分布(ng/g)14.浓度y500(1000)1500)20000 25000750700650600550500450400350300250200150100500图5 Cr的空间分布(p g/g)1900 18001700 16001500 14001300 12001100 10009008007006005004003002001000-1001500)31000)500(0 0 50005000 20000 25000y%图6 Cu的空间分布(p g/g)140001300012000110001000090008000700060005000400030002000100001500)y Q)1000)5000-1000 0 5000 1000*)15000 20000 2500015.图7 Hg的空间分布(ng/g)16.结合该城区五大区域分布等值线图可以得出结论，五大区域分布情况图如图11:150001000050005000 10000 15000 20000 25000图11五大区域分布情况图图中颜色说明：黄色一一生活区 深蓝色一一工业区 浅蓝色一一山区红色一一交通区 绿色一一公园绿地区由图3Tl可以看出，工业区域Cu要素严重超标，此外As、Cd、Pb元素 17.2.3 采样区重金属污染定量分析2.3.1采样区土壤重金属污染指数分析采用单因子指数法和内梅罗综合指数法11分别计算该地区5个区域的重 金属污染指数。计算结果如下表：表2采样土壤重金属污染指数计算结果功 单因子指数 内梅罗能 综合指As Cd Cr Cu Hg N i P b Zn数区18城市附近频繁的人类活动给环境带来较大的影响，污染很严重。主要位于工业区、生活区和交通区。总之，城市周边的工业区生产活动，通区域，这些极高值区主要受特别的人类活动影响。铅的高值区相对较大,N i、P b、Hg和Zn的高值分布基本一致，几个极高值区分散在工业区域和交规律，说明这些元素的空间分布主要受自然因素控制。Cd、As、Cr、Cu、Hg和Zn具有相似的空间分布模式：高值位于工业区，空间模式相对平滑和较轻。从8个变量的空间分布趋势图，可以看出As、Cd、Cr、Cu、N i、P b、超标，几率在60%上（黄色系以上色彩表示），但大部分研究区域污染程度区域的严重污染。其他三种重金属元素，As、Hg及P b,存在特定斑点区域整幅图面上，只有少部分区域超标的概率小于50%,预示了Cd要素在五个Cd、Cr、Cu三幅图上，均出现了概率较高的黄色蓝色区域。特别是Cd,在也存在较明显超标斑块，其余三种重金属超过标准值的斑点相对较少。在.生活 1.7 42 2.230 2.226 3.743 2.658 1.491 2.229 3.435 3.1707区工业 2.014 3.024 1.723 9.662 18.353 1.611 3.001 4.028 13.5330 区山1.122 1.172 1.257 1.312 1.170 1.256 1.179 1.062 1.2531 区主 干道 1.58 6 2.769 1.873 4.713 12.766 1.432 2.049 3.520 9.4262路区公Id绿 1.7 39 2.158 1.408 2.287 3.285 1.243 1.958 2.235 2.7340地区由上表可以看出各个区域均受到了重金属不同程度的污染，其中受污染 程度大小依次为：工业区、主干道路区、生活区、公园区和山区，与我们 问题一所建立的污染综合评价模型吻合。尤其是工业区的Cu、Hg和主干道路区的Hg污染，使得环境质量指数分别达到了 18.353、9.662和12.766,19.名称值值数度景值As7.251.6121.874.2458.482761.66.043.6Cd393.11114.51092.9237.5860.436011.645.41130Cr53.4115.428 5.5844.008 2.381584.2322.8731Cu127.541.272528.48414.94325.34115.6132.9413.2Hg642.3611.79135002244.07349.34775.4431.6835N i19.8 14.2741.78.3742.251390.622.9412.3P b93.0431.24434.88 5.3791.756232.9511.4931Zn27 7.9356.331626.02350.83126.22962.8210.1769该区域土壤中8中重金属元素的平均值均超出了对应元素的背景值，其中元素Hg和Cu的平均值分别为背景值的18.35和9.66倍，说明该地区 受重金属Hg和Cu的污染严重，而As、Cd和P b的平均浓度也超出了背景 值23倍，说明此三种金属的污染也较严重。从整体上说明该区土壤环境 质量已受到相当程度的污染，进而说明该区重金属的污染不容忽视。从820重金属平均最小最大值标准差变异系偏峰度背表3工业区土壤重金属的描述统计结果析。该地区工业区土壤中重金属的统计结果如下表：下面针对污染严重的工业区的统计数据，我们运用统计学的方法进行分2.3.2重点污染区域分析产生的重金属颗粒物的干湿沉降和道路交通。污染已经达到很严重的地步。究其原因，主要重金属污染来源于工业活动.图12工业区As、Cd、Cr Cu元素的浓度与海拔关系图21高，重金属元素的浓度越低。金属元素的浓度比海拔高一些的地区的要高。因此，我们可以得出海拔越元素在山区中的浓度明显比在工业区的浓度高许多，海拔较低的地区的重度以及每个功能区中同一元素在不同海拔的浓度后，不难发现8种重金属的折线图（图12,图13,图14和图15）,对比两种功能区内相同元素的浓以工业区和山区为例，分别作出8种元素在两种功能区内浓度与海拔关系利用MATLAB绘图来分析海拔高度对重金属元素污染浓度的影响。我们2.3.3海拔高度对重金属污染浓度的影响分析正态分布的高峰更加陡峭。元素分布的峰向右倾斜，从其峰度值上看，峰度值均大于0,说明其分布比他区域的环境质量。从偏度值上看，8中元素的偏度值均大于0,说明这些侧面也说明，该地区工业的发展对环境的污染较为严重，甚至影响到了其变化程度的大小。说明该区重金属元素与相邻区域流动较为频繁。从一个异系数均较大，这种变异系数的差异反映了各种元素在当地土壤中含量的中金属的变异系数看，N i的变异系数较小，As和Cd属中等变异，其余变.图图13工业区Hg、N i、P b、Zn元素的浓度与海拔关系ibS a丹 4KM/v4tf*9VMK前as*生X3OCn7VD90谈仪图14山区As、Cd、Cr、Cu元素的浓度与海拔关系22.山1例停费1号3/*|0”电情*界标取图15山区Hg、N i、P b、Zn元素的浓度与海拔关 系图2.4采样区重金属污染传播特征分析2.3.1重金属污染物在各区域间传播的变异性分析利用数据，我们分别计算在不同的区域的重金属污染物的变异系数，如下表（单位：%）表4采样区重金属的变异系数重金属生活区工业区山区主干道路区公园区As变异系数34.2958.4827644.4956.7332.31Cd变异系数63.3560.4360151.4667.618 4.07Cr变异系数156.328 2.3815863.13140.5734.01Cu变异系数95.47325.341161.97193.247 5.1323.24001因素（人为因素）的作用。当C/（c+c）25%时，表明变量的空间变异以结效应.该比值用以反映空间变异影响因素中区域因素（自然因素）和非区域01值c+C之比是反映区域化变量空间异质性程度的重要指标，又称为块金0示系统内的总变异，它是结构性变异和随机性变异之和。块金值C与基台01于最小采样尺度的非连续性变异引起，属于随机变异；基台值C+C通常表0异性15半方差函数在原点处的数值称为块金常数0,它由测定误差和小为分析重金属污染物的传播特征，我们研究重金属污染指数的空间变2.3.2重金属污染指数的空间变异性分析对比较稳定，而其余五种元素变异程度比较大，受人为因素影响较大。流动较大，受人为因素影响较大；公园区内As、Cr、N i变异系数较小，相尤其Cr、Cu、Hg、Zn变异系数都很大，意味着主干道区区内重金属元素的受人为因素影响较小；主干道区内的重金属元素变异程度都在中等以上，异系数均处在中等变异附近，可见山区中的重金属元素变化幅度相对稳定,素的变化程度较大，工业区重金属污染受人为因素影响较大；山区中的变等变异，其余变异系数均较大，说明在工业区Cr、Cu、Hg、P b、Zn五种元活区重金属污染受人为因素影响较大；N i的变异系数较小，As和Cd属中明在生活区土壤内Cd、Cr、Cu、Hg、P b、Zn六种元素的变化程度较大，生可见，生活区中As和N i变异系数较小，其余元素变异系数较大，说Zn变异系数18 7.18126.229642.22158.44149.71P b变异系数104.6691.7562348.5151.207 5.51N i变异系数30.8742.2513967.4866.9132.53Hg变异系数110.6349.347768.0148 7.95195.04.重金属元素c0c+c 0 1c/(c+c)0 0 1As0.0343.9650.008 57 5Cd1.7 8 20.0029614.48 28Cr0.00140.07 980.017 544Cu0.02980.02231.336323Hg0.02130.14510.1467 95N i0.01050.06970.150646P b0.04580.13560.337 7 58Zn0.03540.09120.38 8 158可见该地区重金属元素As、Cr、Hg、N i的空间变异以结构性变异为主,具有很强的空间相关性；P b、Zn为中等程度空间相关；Cd、Cu的空间变异 以随机变异为主，空间相关性很弱。综上分析重金属污染物的传播具有与其他污染物不同的传播途径和传 播规律，并且在传播过程中具有空间变异性。2.5采样区重金属污染主要原因分析为分析采样区土壤重金属污染的主要原因，我们建立因子分析模型1表5采样区土壤重金属污染指数的相关参数染指数的相关参数如下：异为主，变量的空间相关性很弱，利用题目数据，我们可得土壤重金属污001变量为中等程度空间相关；C/(c+o)75%时，变量的空间变异以随机变001构性变异为主，变量具有很强的空间相关性；当25%C/(C+C)X .X为样本相关系数矩阵R的特征值；7!,7，,T1为相应的 1 2 p 1 2 p标准正交化特征向量。设机 j1=-X 4=1/=1-（X 一1）2 ij ii=j=l i=lx和1为分别表示位置1和）处的样本属性值，？为其平均值，W为权 i J ij重矩阵，如果和,在某一阈值范围内W=l;否则，W=0 oij U因子分析法的计算一般经过以下步骤：27/加个公共因子对第个变量方差的贡献称为第，共同度，记为近，.28P b50.0645.8 361.7 447 2.4819.680.8131N i9.9415.9817.26142.54.270.5812.3Hg1629.5450299.7 1160008.575.4435Cu162.9227.9855.022528.482.292.9613.2Cr7 042.0253.51920.8415.321.3131Cd224.99238.7 0302.401619.8400.74130As3.025.315.6830.131.610.533.60元素标准差中值平均值最大值最小值变异系数背景值表6重金属元素富集因子统计表4给出了采样区主要重金属元素的一些基本统计量。及程序运行结果可参见附录。对于因子分析法的计算过程我们通过MATLAB软件进行了编程，程序以（4）确定因子个数，计算因子得分，进行统计分析。而且旋转后的因子仍然正交；（3）进行正交变换，使用方差最大法，其目的是使因子载荷两极分化,征向量；（2）计算标准化数据的相关系数阵，求出相关系数矩阵的特征值和特数；在于消除不同变量的量纲的影响，而且标准化转化不会改变变量的相关系jj本的第j个指标值，而x和6分别为j指标的均值和标准差，标准化的目的jy8ij（1）原始数据的标准化，标准化的公式为z=工:，其中x为第1个样.Zn339.23106.43 1.6969运用MATLAB软件进行因子分析后，可以得到以下结果。首先给出了该城区表层土壤As、Cd、Cr、Cu Hg N i、P b、Zn这8种重金属元素原始含量 数据的相关系数矩阵，如表6所示：表7变量相关矩阵指标AsCdCrCuHgN iP bZnAs1.00000.25470.18 900.15970.06440.31660.28 990.2469Cd0.25471.00000.35240.39670.26470.32940.66030.4312Cr0.18 900.35241.00000.53160.10320.7 1580.38 280.4243Cu0.15970.39670.53161.00000.41670.49460.52000.38 7 3Hg0.06440.26470.10320.41671.00000.10290.298 10.1958N i0.31660.32940.7 1580.49460.10291.00000.30680.4364P b0.28 990.66030.38 280.52000.298 10.30681.00000.4937Zn0.24690.43120.42430.38 7 30.19580.43640.49371.0000利用matlab作聚类分析图如下：（图16）29201.237 60.8232.86.图16 土壤重金属相关元素聚类分析图据表据口图和可见,Cr和N i的相关性最好，相关系数最大，为0.7 158,其次为P b和Cd,相关系数为0.6603,以下依次是Cr和Cu,P d和Cu的相关性 较好，相关系数分别为0.5316和0.52,其它元素之间的相关性一般，从成 因上来分析，相关性较好的元素可能在成因和来源上有一定的关联。因子分析的关键就是利用相关系数矩阵求出相应的因子的特征值和累 计贡献率，我们利用MATLAB统计软件计算得出了表8:表8特征值和累计贡献率旋转前旋转后因总的特 子占总变量的累计贡献总的特占总变量的累计贡献征值百分率/%率/%征值百分率/%率/%F13.560044.499944.49991.826722.830522.8305F21.150214.377458.87731.045413.071135.901630.F30.965112.06337 0.94061.031812.901748.8033F40.76779.59618 0.53671.525119.060067.8633F50.57767.21978 7.75641.003912.55218 0.4154F0.43205.399493.15581.019512.740493.1558在累积方差为93.1558%(90%)的前提下，分析得到6个主因子，可以看到6个主因子提供了源资料的93.1558%的信息，满足因子分析的原 则，而且从上表可以看出旋转前后总的累计贡献率没有发生变化，即总的 信息量没有损失。同时从表还可得出，旋转之后主因子1和主因子4的方差 贡献率均为20%左右，主因子2、主因子3、主因子5和主因子6的方差贡献 率的范围是12.5521%到13.0711%,这可以解释为因子1和因子4可能是该 城区土壤重金属污染的最重要的污染源，因子2、因子3、因子5和因子6对 该城区重金属污染也具有一定的影响。因子分析的主要目的是将具有相似性的因子荷载的各个变量置于一个 公因子之下，正交方差最大旋转使每一个主因子只与最少个数的变量有相 关关系，而使足够多的因子负荷均很小，以便对因子的意义做出更合理的 解释，得到的结果见表9和表10：表9旋转前因子载荷矩阵指标F F1 2F3F4F5F6As0.4257-o.19960.68 100.5508-0.0263-0.0650Cd0.7 107 0.28 140.28 24-0.3221-0.25430.3250Cr0.7 350-0.4440-0.3035-0.0462-0.10960.098 231.Cu0.7 564 -0.1545-0.408 1Hg0.408 4 0.15410.2356N i0.7 228 -0.01390.1997P b0.7 640 -0.157 6-0.2167Zn0.698 8 0.6539-0.0598表10方差极大正交旋转后因子载荷矩阵指标F F1 2F3F4F5F6As0.13510.01600.97 420.12960.08 35-0.0427Cd0.22260.148 30.08 590.917 90.1256-0.0306Cr0.8 590-0.01600.00370.18 7 40.1411-0.2450Cu0.39490.27 310.03110.17 680.1096-0.8 097Hg0.01690.967 10.01630.137 40.0694-0.17 30N i0.8 8 950.05140.19630.10110.1695-0.1031P b0.07 340.06130.16190.7 17 80.27 52-0.4956Zn0.26100.08 190.09660.23480.9165-0.1219由表9和表10可知，旋转前后因子荷载的变量结果基本一致。变量与某 一个因子的联系系数绝对值越大，则该因子与变量关系越密切。正交因子 解说明：因子1为Cr和N i的组合，因子2为Hg,因子3为As,因子4为Cd和P b 的组合，因子5为Zn,因子6为Cu,所以认为Cr和N i、Cd和P b可能是同一个 来源，而且这两组元素正是相关性最好的两组元素。32-0.037 40.1232-0.24100.31430.237 2-0.247 8-0.5146-0.18 980.13670.67 34-0.297 50.448 70.1247-0.36530.137 5.图17第一个主因子图19第