论文定3Word格式.docx

1、CdHg.在生活区、公园绿地区、山区为Zn Hg。关键词： 重金属污染 MATLAB 因子分析法 Excel绘图一、问题的背景与提出随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表

2、层土（010 厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。根据取样得到的各种信息，即采样点的位置、海拔高度及其所属功能区、8种主要重金属元素在采样点处的浓度以及8种主要重金属元素的背景值。我们就是要通过建立合理的模型来完成以下的任务：1. 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。2. 通过数据分析，说明重金属污染的主要原因。3. 分析重金属污染物的传播特征，由此建立模型，确定污

3、染源的位置。4.分析建立模型的优缺点，为更好地研究城市地质环境的演变模式，收集相应的信息。根据相关的信息，建立合理的模型来解决问题。二、模型的假设1、假设在每一个采样点采集的土壤样品都能充分的代表该点周围土壤所含重金属的浓度，且测量的结果比较准确。 2、假设在远离人群及工业活动2公里的自然区所取的土壤样品中，测得的重金属元素含量与标准土壤中重金属元素的含量相同。 3、假设在采集样本点的土壤所含重金属浓度不受其它区域的影响。 4、假设研究的各地区重金属污染的浓度相对稳定。 5、假设在查证获得的海量数据中能充分的反应该区域的重金属污染程度。6、假设采样点土壤表面重金属污染物分布均匀。7、 假设各个

4、区域的重金属污染物不发生大幅度的变化。三、符号说明符号意义综合污染指数因子分析法的相关系数土壤重金属元素质量分数综合污染程度土壤重金属元素质量分数背景值特征根重金属元素的单项污染系数特征向量重金属元素的实测质量分数计算土壤参数所需比值重金属元素的毒性系数四、模型的建立与求解1.模型一的建立与求解：（1）模型一的建立与求解：为了得到8种主要重金属元素在该城区的空间分布，我们将数据进行统计分析，得到如表一所示的结果。表一:该城市表层土壤重金属浓度元素最小值最大值平均值背景值标准差变异系数累积强度As（ug/g）1.6130.135.683.60.90.1581.578Cd（ng/g）40 1619

5、.80302.4130300.0992.326Cr（ug/g）15.32920.8453.513190.1681.726Cu（ug/g）2.292528.4855.0213.20.0654.168Hg（ng/g）8.5716000299.718350.0278.563Ni（ug/g）4.27142.5017.263.812.30.2201.403Pb（ug/g）19.68472.4861.7460.0971.992Zn（ug/g）32.863760.82201.2014690.072.916注：累积强度=平均值/背景值 变异系数=标准差/平均值由表一可知，与该城区的土壤背景值相比，As,Cd,

6、Cr,Cu,Hg,Ni,Pb,Zn各元素在该城区发生了明显的积累，已经构成不同程度的污染。其累积强度依次为HgZnNi。W（Ni）平均为17.26ug/g，与土壤背景值12.3ug/g基本持平，由此可以看出Ni的人为累积量较小，但其最大值明显高于土壤背景值，说明部分区域出现富积现象。As,Cr,Pb的累积强度为1.578,1.726,1.992，说明他们在很大一部分区域存在一定程度的富积。其余四种重金属元素的累积强度都很大，说明在该城区的富积程度很高。而其在不同的区域，富积程度有一定的差异。由于各重金属元素间的变异系数差别很大，他们对土壤的污染有着较大的差异。我们将附件中的数据分在各个区域进行

7、处理，求得各种重金属元素在各个区域的平均浓度以及八种重金属元素在这个城区的平均浓度。通过这些平均值，我们在Excel中分析数据，通过插入图表的方式，将八种重金属元素在这五个区域的分布图描绘出来。以下几个图形即为八种主要重金属元素在该区域的空间分布图。 从图中可以看出Hg,Cd,Zn在各个区域的质量分数相对较大，他们的变化较为明显。而且从图中可以知道在工业区，交通区比较明显。这可能与“三废”的排放，人流量，车流量有着较大的关系。无论在那个区域Ni和As的含量都相对较小而且变化都不明显。根据空间分布图，我们可以知道各区域内重金属的污染程度。表现为工业区最为严重，其次是交通区，生活区，工业绿地区，最

8、后是山区。其中在工业区和交通区八种重金属的污染程度表现为znCd Hg.在生活区、工业绿地区、山区为Zn 根据上述分析可以知道在不同的区域，他们的污染源有所不同。（2）模型一的改进： 为了进一步的说明不同区域内重金属的污染程度，我们通过用内梅罗综合污染指数全面的反映各重金属元素对土壤的影响。这样兼顾了在不同区域内单元素污染指数平均值和最大值，突出了污染指数最大的污染物以及最严重的污染区域。其计算公式为： （1）其中土壤污染评价标准为：P0.7为安全；0.71为警戒线；12为轻度污染：2表三是根据（1）式而得到的污染指数，而图二则是根据表二画出的图形。AsCdCrCuHgNiPbZn生活区2.1

9、034.11811.4897.3827.071.877.38816.182工业区3.3454.0325.11443.107652.2066.1058.778山区2.1872.2093.6373.383.953.92.5322.405交通区4.8315.64515.40633.84390.9196.923.23420.694公园绿地区2.134.092.10211.3862.7581.7393.9579.654表二：各区污染指数表图二：各重金属污染指数根据表二和图二我们得知只有少部分重金属在部分区域出现轻度污染。其中山区，公园绿地区为中度污染，生活区污染表较严重，交通区，工业区为重度污染。2.模

10、型二的建立与求解：（1）模型二的建立与求解：通过数据的分析以及模型一的相关信息，我们知道工业区和交通区的污染最严重；公园绿地区和山区的污染较小。根据公式 （2） （3）我们可以得到如表四所示的五个区域的单项污染系数以及综合污染程度。城市分布pb1.7422.232.2263.7422.6581.4912.2293.43519.7532.0143.0241.7239.66218.3531.6113.0014.02843.4161.1221.1721.2571.3121.1701.2561.1791.0629.531.5862.7691.8734.71312.7661.4332.0493.5230

11、.7092.1581.4082.2873.2851.2431.9582.23516.313表三：土壤重金属污染单项系数和综合污染程度从表三可以明确的知道他们的综合污染程度由重到轻依次是：工业区，交通区，生活区，公园绿地区，山区。这与模型一中的金属污染程度的结论完全吻合，由此不仅说明模型一的准确型，而且我们可以通过这样的数据分析，进一步说明重金属污染的主要原因。由于工业区的污染程度最高，我们可以知道工厂的“三废”排放是这个区域内重金属污染的最主要原因。工厂排放的大量含有Hg，Cd，zn，Cu等重金属的灰尘降落地表，导致该区域污染严重。交通区的污染严重，是因为交通区除了较大的车流量外人流量也比较大

12、。Tamd等研究发现在交通区内的土壤中重金属元素的含量与交通流量有显著的相关性。因此交通流量是交通区了重金属污染的主要原因。刘廷良等认为，汽车轮胎的添加剂中含有Zn这种重金属，因此轮胎磨损产生的粉尘是土壤中Zn污染的主要来源。在生活区中，人们的活动比较频繁。各种形式的活动都会在一定程度上对该区域的土壤产生影响。生活区中市民垃圾的堆放，各种废弃物的排放会导致Zn，Cr,Pb,Cu等重金属的出现。在公园绿地区和山区，虽然污染程度不算太高，但还是受到一定程度的重金属污染。在这两个区域中无论是累积强度还是单项污染系数都只有一两种重金属的较高。从这样的数据分析我们可以知道，这两个区域内重金属污染的主要原

13、因是人流量。以上都是我们从不同的区域来看的，为了进一步的说明重金属污染的主要原因，我们在模型二的改进中从各种元素的分布情况来分析重金属污染的主要原因。（2）模型二的改进:为了进一步的说明重金属污染的主要原因，我们从原始数据中分析各种重金属的分布情况，可以进一步的说明污染的原因。As是浓度最低的重金属，然而在生活区和公园绿地区的平均含量是最高的，由此我们可以得到As主要来自土壤母质。而且生活区和公园绿地区绿化面积大，一些树木的年龄也比较大，为了促使树木健康成长，会人为的施加含As的化肥。经过常量的累积，As的含量就会逐渐的增加，导致As在这两个区域平均含量最高。这就是As这种重金属污染的主要原因

14、。Cd的浓度是最高的，而且在生活区，交通区，公园绿地区污染较严重。它主要来自污水的灌溉，垃圾的堆放，大量的使用农药，化肥等带来的污染以及工业废水的排放。在生活区，居民主要以务农为业，经过长期的耕种，施加农药化肥以及灌溉污水就会导致Cd的含量增加。主要干路区以及交通区的Cd含量也较高。而在公园绿地区中，公园与花园中大量的绿化过程中污水污泥堆肥的广泛使用也会导致Cd含量的增加。这就是Cd这种重金属污染的主要原因。从数据中我们可以的到Cr不严重，平均最低浓度为22.9，最高为75.37.它的污染在生活区较明显，由此可以说明Cr的污染主要是来自土母质。Hg的污染也很严重，工业区或者是大型的污染企业是H

15、g污染单一的来源，“三废”的排放是该重金属污染的主要来源。在工业区和交通区，都会导致Hg含量的增加。比如火力发电，冶炼金属，水泥的制造，垃圾的处理等都是Hg污染的主要原因。Ni的污染也不是太严重，主要就是来自土壤的母质。Pb的来源广泛，“三废”的排放，污水的灌溉和交通运输都是其最主要的来源。在工厂中，大量的烧煤，会增加Pb的含量。在各种交通运输中，汽油，柴油等得燃烧会导致Pb这种化合物的大量排放，最终导致土壤严重污染。有机肥，化肥和农药大量的使用，是土壤中Cu和Zn污染的主要途径。在矿工企业，农业活动以及交通运输的共同作用下，使得这两种重金属含量增加。因此，在工业区，交通区表现得比较明显。在这

16、个模型中，我们首先通过对五个区域进行分析，得出污染的主要原因。在模型的改进中，我们对八种元素进行了分析，也得出污染的主要原因。通过这两方面的分析综合我们最终可以得到重金属污染的主要原因。“三废”的排放，土壤母质，污水灌溉，大量的使用农药、化肥、有机肥，各种物质的燃烧，人流量，车流量，汽车尾气的排放等就是这八种重金属污染的主要原因。3.模型三的建立与求解：利用因子分析法定量解析重金属元素对土壤污染的特征来确定污染源的位置：因子分析法的因子变量具有明显的解释性，可以最大限度地发挥专业分析的作用，它把众多的指标综合成几个为数较少的指标，这些指标即因子指标，对变量间相关关系探测，将原始变量进行分类，将

17、相关性高的变量分为一组，用共性因子代替该组变量。通过对多个变量的相关系数矩阵的研究找出同时支配所有变量的主因子，将城区作为主要影响变量的因子，从而通过重金属的传播特征，根据这样的方法，找出给地区的污染源，进一步确定污染源的位置。如下表所示，我们将各个区域作为主因子，将重金属污染的来源作为非公性因子：表四：各因子高载荷指标因子高载荷指标（非共性因子）因子命名第一主因子生活废水生活区F1生活垃圾第二主因子工业废水工业区F2工业废气工业废渣第三主因子农业农药产生山区F3第四主因子交通工具废气排放交通区F4第五主因子人们娱乐产生垃圾公园绿地区F5宠物粪便根据选取的五个主因子，构建重金属污染的综合指标指

18、数F （4）Fi， i=1.5为各主要因子的表达式，表示第i个因子其方差贡献率的比重 、经计算得各重金属元素在城区的平均值如下表所示：As （g/g）Cd （ng/g）Cr （g/g）Cu （g/g）Hg （ng/g）Ni （g/g）Pb （g/g）Zn （g/g）6.2700289.960069.020049.400093.040018.340069.1100273.01007.2500393.110053.4100127.5400642.360019.8100277.93004.0400152.320038.960017.320040.960015.420036.560073.2900公路

19、区5.7100360.010058.050062.2100446.820017.620063.5300242.85006.2400280.540043.640030.1900144.990015.290060.7100154.2400表五：各重金属元素在城区的平均值相关系数：用MATLAB软件（R=corrcoef（A）将数据标准化后的相关系数矩阵为用MATLAB软件（V,D=eig（C1）将数据标准化后的相关系数矩阵求得特征根和特征向量分别如下：将这一分析推广至八个变量存在五个潜在公因子可建立一下模型：用MATELB软件（A=COV（R）求出方差贡献率：土壤重金属污染元素指标的载荷、特征根及

20、方差贡献率矩阵如下表：0.00010.00950.00070.00390.00860.86110.06070.33370.16890.01420.01860.01270.18380.07960.01972.08510.06511.02420.41280.00020.01400.00160.00750.00340.00230.04080.72000.09080.28530.1558特征值0.00270.00290.10270.79114.1006方差贡献率0.000540.000580.0220540.0158220.82012累计方差贡献率0.001120.023170.034750.8548

21、7方差贡献率=（为特征值）由主因子分析法知污染物累计方差贡献率达85.487可知研究的主因子具有很高的依赖性，影响污染物传播途径的有主因子，共性因子，非共性因。由此可知污染物的主要来源，生活区主要带来污染的重金属是Hg，在工业区主要带来Cd，Cu，Hg,Pb，Zn等重金属元素的污染，在山区主要带来的是Pb，Zn等重金属污染（主要原因是农作物使用农药造成的）在交通区主要是CU，Hg,Zn等重金属的污染（主要原因是轮胎中含有大量的Zn，由于摩擦产生污染），在公园绿地区主要带来污染的重金属是Cd，Hg.Zn .由上述表格所得数据分析，通过对城区单项重金属测定就可以初步知道重金属污染物的主要来源，由此

22、可以得出工业对城区土壤重金属污染有非常大的关系。所以各重金属污染源的位置为工业区，交通区附近。4.模型四的建立与求解：在问题一中，我们采用图表和数据分析，通过计算得到各重金属元素的变异系数，累积强度等进一步的说明了不同区域重金属的污染程度。这样将图形数据结合在一起解决问题，大大提高了结论的准确度。在问题二中，我们通过计算单项污染系数和综合污染指数，更加的说明了问题结论的准确性。而其我们首先通过对五个区域重金属相关数据的分析，得到了在这些区域内重金属污染的主要原因；我们又分析了八种重金属污染的主要原因以及分布的区域， 问题一二我们都是通过分析数据，结合图表而得到结论的。这样的模型其优点在于通过数

23、据与图形结合，增加了结论的准确率，使得它更具有说服力。其缺点就是有一定的局限性。问题三中，我们采用了主因子分析法，这种方法的优点在于各综合因子的权重不是主观赋值而是根据各自的方差贡献率大小来确定的，方差越大，变量越重要，从而具有较大的权重；相反，方差较小的权重所对应的权重也就越小。这样就避免了人为确定权重的随意性，使得结论唯一，而其较为客观合理。此外，因子分析的整个过程运用了计算机软件方便快捷的进行，可操作性强。因此，这种方法是一种科学、实用、简便的方法。主因子分析法的缺点在于将变量减少，而其本身的工作量大，在分析时，是以整体进行分析，对基础的数据的准确度要求较高。而其在分析中数据出现错误时不易被发现，这就需要分析时前后进行核对以及对数据的推理复核。这会在一定程度上增加工作量。为了更好的研究城市地质环境的演变模式，我们应该从许多客观因素考察及其分析，并总结得出相关合理的考察依据。人类活动导致自然生态系统向人为生态系统转变已经有了很多年的历史，这相对于城市的地质环境演变