城市表层土壤重金属污染分析 建模论文完整资料doc文档格式.docx

1、现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（010 厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。建立数学建模来完成以下任务：（1） 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重

2、金属的污染程度。（2） 通过数据分析，说明重金属污染的主要原因。（3） 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。（4） 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二、问题分析问题一想要得到重金属元素在该城区的空间分布，并研究该城区内不同区域重金属的污染程度。由于所给采样点不是位于同一个区域，这就需要对所给采样点的数据进行分区域筛选和整合。由于还需要分析不同区域重金属的污染程度，污染程度需要抽象为污染指数，通过比较污染指数与国家背景值，确定污染程度。通过查阅重金属污染程度评价方法的相关资料，主要方法有单

3、因子指数法、内梅罗指数法、地质累积指数法1。问题二中所给元素浓度单位不统一，在进行重金属污染的主要原因分析时会受到量纲和数量级的影响，因此对附件2各区域8种重金属元素浓度进行处理。经过数据的分析，重金属污染的主要原因取决于重金属在各个区域中的含量，由于每个区域的重金属含量不同，我们首先需要确定哪些重金属元素为污染环境的主要因素。主要因素的确定方法有很多，如聚类分析法、因子分析和主成分分析法。问题三为了确定污染源的位置，可以根据污染源的传播特征建立模型来寻找局部最高点。题干研究的是城市表层土壤重金属的污染源，重金属进入土壤的主要方式有空气扩散、重力沉积，水流渗透。通过地质变化原理可假设浓度最大的

4、点是污染源的的位置，因此确定污染源位置即寻找浓度最高点。问题四通过对城市地质环境的相关的资料的查阅对问题三所建立模型进行进一步研究，加入更多可测数据来使模型更加完善精确。三、问题假设（1）假设不考虑元素间的相互作用的影响。（2）假设不考虑历史沉积的重金属的影响。（3）短期内重金属元素的物理、化学变化及迁移对周围环境影响不大。（4）假设用国家土壤环境质量二级标准值与不同功能区各种重金属的实测平均浓度进行对比后得出的评价标准值能准确反映土壤质量情况。（5）各污染源附近区域，其重金属污染源的浓度主要受到该污染源的影响，其他污染源影响较小。四、变量说明符号含义Iij第i个区域元素j的地质累积指数Cij

5、第i个区域元素j的平均值BEj元素j地球化学背景值Sj元素j的国家评价标准Lij第i个区域元素j的分指数Pi第i个区域的综合污染指数Zi第i个原始自变量Fi第i个公因子五、模型的建立、求解和结果分析5.1.问题一5.1.1数据处理及分析由于采样点分布的离散程度大，首先需要应用Excel软件筛选附件1所给采样点数据对其进行筛选，将其整合到五个区域中。利用Matlab软件中三次多项式插值拟合函数绘制出5类区的区域地图（图5.1.1）和三维海拔图（图5.1.2图5.1.1 5类区区域地图然后根据附件2中所给采样点的8种主要重金属元素的浓度应用Matlab软件中三次样条插值拟合函数得到8种主要重金属元

6、素在该城区的区域分布图见图5.1.3-5.1.10。图5.1.2 区域三维海拔图图5.1.3 As元素区域分布图图5.1.4 Cd元素区域分布图图5.1.5 Cr元素空间浓度分布图图5.1.6 Cu元素空间分布图图5.1.7 Hg元素空间分布图图5.1.8 Ni元素空间分布图图5.1.9 Pb元素空间分布图图5.1.10 Zn元素空间分布图5.1.2 地质累积指数模型的建立地质累积指数不仅反映了元素的自然变化特征，而且可以判别人为活动对环境的影响，是区分人为活动的重要参数。地质累积指数法：其中Iij为第i个区域元素j的地质累积指数，Cij为第i个区域元素j的平均值，BEj是元素j地球化学背景值

7、，为各地岩石差异导致的背景值变动系数（一般为1.5）。按受污染程度的强弱，将地质累积指数分为7个级别，0-6级污染程度由弱到强递增，具体分级标准见表5.1.1。表5.1.1 地质累积污染指数分级标准地累积污染指数I值级别污染程度I0无污染0I11轻度污染1I22中度污染2I33中强度污染3I44强度污染4I55强污染到极强污染I56极强污染5.1.3综合污染指数模型的建立单因子指数法可以判断出环境的主要污染因子。土壤环境是一个复杂的体系，土壤污染往往是由多个污染因子复合污染导致。单因子污染指数法在反映本题城区土壤污染程度不够全面综合。因此需将单因子污染指数按一定方法综合起来进行评价，即应用综合

8、污染指数法评价。内梅罗指数法是先求出各因子的分指数，然后求出各分指数的平均值，取最大分指数和平均值计算。将单因子指数法和内梅罗综合指数法结合建立综合污染指数模型2。首先建立单因子污染指数式为：然后应用内梅罗指数法得到各区域新的综合污染指标Pi按受污染程度的强弱，将综合污染指数分为5个级别，1-5级污染程度由弱到强递增，具体分级标准见表5.1.2。表5.1.2 综合污染指数分级标准综合污染指数P值污染水平P0.7安全清洁0.7P1.0警戒线尚清洁1.0P2.0污染物超过起初污染值，作物开始污染2.0主干道路区生活区公园绿地区表5.1.4 生活区8种重金属的分指数3.66.27300.101153

9、5810.0130.290.30.9650221733169.022500.17360004213.249.40351.6606964140.0350.090.21902144112.318.34400.21813258969.110.1740016669237.012001.282508675表5.1.5 工业区8种重金属的分指数7.250.1383101850.391.32442895953.410.102324962127.545.2447629970.642.29190670919.810.27117930293.040.283291476277.931.594866412表5.1.6

10、 山区8种重金属的分指数4.040.0166607480.150.48543448438.960.03634564817.320.1888657210.040.022475715.450.11385515836.560.02536944873.290.032780477表5.1.7 主干道路区8种重金属的分指数5.710.0798501320.361.20910973158.050.1235338562.212.2483911710.401.38254145217.620.19195312163.530.148558004242.851.327136298表5.1.8 公园绿地区8种重金属的分

11、指数6.260.1008982680.280.93220507743.640.0576986330.190.7794364350.110.30185552615.290.10793192460.710.135655577154.240.650704471表5.1.9 各区域综合污染程度评价LmaxLave综合污染指数P0.5992670721.2484057045.2447631.4063843.83960.4854340.1152230.3527912.2483910.8388841.6969070.9322050.3832980.712714警戒级5.2.问题二5.2.1数据处理及分析题中

12、所给元素浓度单位不统一，在进行重金属污染的主要原因分析时会受到量纲和数量级的影响，因此对附件2各区域8种重金属元素浓度进行无量纲化处理。采用均值化方法处理，每一个变量除以该变量的平均值，即标准化以后各变量的平均值都为1，标准差为原始变量的变异系数。该方法在消除量纲和数量级影响的同时，保留了各变量取值差异程度上的信息，差异程度越大的变量对综合分析的影响也越大。然后取其平均值进行对比。采用均值化方法无量纲化处理后，选取8种重金属元素的平均值，用Excel软件画出5类区的8种重金属元素所占百分比饼状图见图5.2.1-图5.2.5。图5.2.1 生活区8种重金属元素所占百分比由图5.2.1可知，生活区

13、Cu元素和Cd元素所占比例较高。图5.2.2 工业区8种重金属元素所占百分比由图5.2.2得，生活区Hg元素所占比例最高，Cu元素次之，是重金属污染的主要原因。图5.2.3 山区8种重金属元素所占百分比观察图5.2.3可得，山区所有重金属元素所占百分比相近。图5.2.4 主干道路区8种重金属元素所占百分比由图5.2.4得，主干道路区Hg元素所占比例最高，Cu元素次之，是重金属污染的主要原因。图5.2.5 公园绿地区8种重金属元素所占百分比由图5.2.5可知，生活区Hg元素所占比例较高。5.2.2 主成分分析模型的建立确定主要因素的方法有很多，如聚类分析、因子分析和主成分分析法。主成分分析法最为

14、客观，是通过对现有数据的分析处理得出主要因子，因此本题选用主成分分析法来确定重金属污染的主要原因。主成分分析的实质是降维，即把多指标转化为少数几个综合指标。在用统计方法研究多变量问题时，变量太多会增加计算量和增加分析问题的复杂性，人们希望在进行定量分析的过程中，涉及的变量较少，得到的信息量较多。设法将原来的变量重新组合成一组新的互相无关的几个综合变量，同时根据实际需要从中可以取出几个较少的综合变量尽可能多地反映原来变量的信息的统计方法叫做主成分分析或称主分量分析，也是数学上用来降维的一种方法3。主成分分析法的计算步骤如下：Step 1:构造样本矩阵，对样本数据标准化；Step 2:对标准化阵Z

15、求相关系数矩阵，公式如下： （5.2.1）（i,j=1,2,p） （5.2.2）Step 3:解样本相关矩阵R的特征方程得到p个特征值，确定主成分; （5.2.3） Step 4:将标准化后的指标变量转换为主成分Fi再对主成分进行加权求和，即得最终评价值。5.2.3模型的求解及结果分析要研究重金属污染的主要原因，首先我们通过将经过无量纲化处理的各区域8种重金属元素浓度数据导入SPSS软件进行主成分分析处理，得到各区域重金属元素的主成分荷载矩阵及方差分析表（见表5.2.1-表5.2.5）。然后根据权重系数值大小降序排列，系数大的对应的重金属元素就是主要污染的主要原因。各区域重金属元素权重系数见表

16、5.2.6。表5.2.1 生活区主成分荷载矩阵及方差分析成分共因子方差初始提取0.669-6.46-0.10.8650.7840.171-0.4170.8180.6430.2340.4930.7120.729-0.2460.0240.5920.4920.13-0.4370.450.686-0.2530.5230.8090.8030.112-0.3480.7780.5010.6910.2670.8表5.2.2 工业区主成分荷载矩阵及方差分析-0.170.9180.8440.5970.5150.6230.3630.9670.9750.1260.9240.9540.1150.7650.3810.7880.6710.5370.7390.5910.6520.774表5.2.3 山区主成分荷载矩阵及方差分析0.006-0.2020.4030.674-0.060.393-0.0050.840.41-0.065-0.030.9430.0430.488-0.1740.1750.430.6680.423-0.0820.0530.959-0.0580.3790.0230.7930.230.0510.883表5.2.4 主干道路区主成分荷载矩阵及方差分析主干道路