《城市表层土壤重金属污染优秀论文》由会员分享，可在线阅读，更多相关《城市表层土壤重金属污染优秀论文（32页珍藏版）》请在金锄头文库上搜索。
1、城市表层土壤重金属污染分析摘要本文主要是由采样数据对该城市八种重金属元素污染的程度进行分析.首先我们对采样数据进行预处理,主要是修改异常数据和利用 图做正态性检验.然后对第一问我们对八种重金属含量分成五段分别来表示不同的浓度,再利用Matlab软件根据取样点的坐标,做出各种重金属元素的空间分布图,并用不同的颜色表示重金属元素的浓度图中重金属浓度由高到低分别用黑、蓝、红、黄、绿表示.第二问为确定污染原因,我们从两方面来坐分析,一方面利用Excel对五种功能区作单因子方差分析,得出不同功能区的污染程度有非常显著性的差异,并得到绝大数都是由于工业污染造成的；另一方面,求出它们的尼梅罗综合评价指数也说
2、明绝大数都是由于工业污染造成的.第三问为找出传播特征和确定污染源,我们主要是利用地统计学的半方差函数或称变异函数模型,先通过GS+V9.0软件拟合估计出半方差函数的三个参数,并由此确定各重金属污染的传播特征,再进行克里格插值来推出它们的污染源,另外我们还求出各种重金属元素以与海拔间的相关系数,从相关性的角度说明它们的传播特征和污染源的确定.最后说明该模型的缺点以与需要改进的地方.关键词：正态性检验；单因子方差分析；半方差函数；克里格插值；一、 问题提出随着城市经济的快速发展和城市人口的不断增加,人类活动对城市环境质量的影响日显突出.对城市土壤地质环境异常的查证,以与如何应用查证获得的海量数据资
3、料开展城市环境质量评价,研究人类活动影响下城市地质环境的演变模式,日益成为人们关注的焦点.按照功能划分,城区一般可分为生活区、工业区、山区、交通区与公园绿地区等,分别记为1类区、2类区、5类区,不同的区域环境受人类活动影响的程度不同.现对某城市城区土壤地质环境进行调查.为此,将所考察的城区划分为间距1公里左右的网格子区域,按照每平方公里1个采样点对表层土0-10 厘米深度进行取样、编号,并用GPS记录采样点的位置.应用专门仪器测试分析,获得了每个样本所含的多种化学元素的浓度数据.另一方面,按照2公里的间距在那些远离人群与工业活动的自然区取样,将其作为该城区表层土壤中元素的背景值.附件1列出了采
4、样点的位置、海拔高度与其所属功能区等信息,附件2列出了8种主要重金属元素在采样点处的浓度,附件3列出了8种主要重金属元素的背景值.现要求你们通过数学建模来完成以下任务： 给出8种主要重金属元素在该城区的空间分布,并分析该城区内不同区域重金属的污染程度. 通过数据分析,说明重金属污染的主要原因. 分析重金属污染物的传播特征,由此建立模型,确定污染源的位置. 分析你所建立模型的优缺点,为更好地研究城市地质环境的演变模式,还应收集什么信息？有了这些信息,如何建立模型解决问题？二、 问题分析本文对该城市八种金属元素污染程度的分析,首先是主要将异常的数据进行处理,我们采用统计中常用的原则,即以采用数据超
5、出范围的为异常数据,其中 、 分别为采用数据中各重金属的样本均值和样本标准差.当采用数据大于时,以代替；当采用数据小于时,以代替.然后利用SPSS作 图做正态性分析.处理完数据之后,再利用Matlab软件根据取样点的坐标,做出各种重金属元素的空间分布图,并用不同的颜色表示重金属元素的浓度,从而解决第一问.对于第二问为确定重金属污染的主要原因,我们对城市五种不同的功能区进行八种重金属污染指数的单因子方差分析,分析这五种功能区的污染程度是否有显著性差异,若差异很显著,则主要看是那个因子即功能区影响最大,从而确定污染的原因,另一方面比较这八种元素在不同功能区的尼梅罗综合指数,找出最低最高的尼梅罗综合
6、指数,从而根据该数据也可以得到第二问的答案.对于第三问,为分析出重金属污染物的传播特征,并由此确定污染源的位置,我们主要根据地统计学半方差函数模型和克里格插值法来做的.其总的思想就是根据已知采样点提供的信息拟合半方差函数,然后根据该函数对未知点进行估计和模拟.因此,半方差函数是我们利用的主要工具,一方面,利用半方差函数对重金属的空间分布进行结构分析和变异性分析,并由此得出重金属污染物的传播特征；另一方面应用前面分析的结果,利用克里格法进行估值,并由此对其它未知点污染程度进行预测得出图形,并由此确定污染源污染最严重的区域.当然为了更好确定传播特征和污染源的位置,我们对各重金属元素和海拔的相关性分
7、析,通过数据不难看出他们的传播特征在一定程度上也相似,但是也有些相关性相对较弱,如As和Cd,说明他们的传播特征有一定的区别.但总体上,重金属的含量之间都是正相关的,说明所有重金属的传播有一定的共同点.而所有重金属的含量与海拔之间都是负相关的,说明,随着海拔越低,他们的重金属含量越高,表明重金属的传播向下沉. 最后对所建模型进行评价和改进.三、 模型的假设1) 假设所有数据都是真实可靠的；2) 假设污染源为重金属含量最高的地方；3) 假设问题三中重金属含量是平稳的,即假设第三问中半方差函数的值只与距离h有关,而与其他无关.四、 符号说明：在图中表示工业区；：在图中表示生活区；：在图中表示交通区
8、；：在图中表示山区；：在图中表示公园绿地区；：变异函数；：是分隔距离为时的所有观测样本对总数；：为块金常数；：为拱高；：为变程.五、 模型建立与求解一采样数据的预处理1对异常数据的处理异常数据值在第三问中对重金属污染的空间变异特性有很重要的影响,会影响变异函数理论模型的精度.我们采用统计中常用的原则,即以采用数据超出范围的为异常数据,其中、分别为采用数据中各重金属的样本均值和样本标准差.当采用数据大于时,以代替；当采用数据小于时,以代替.2对采用数据的正态性检验本文的第二问的各功能区重金属污染的单因子方差分析和第三问的重金属传播的空间变异特征的研究模型都是基于正态分布的,因此有必要对八种重金属
9、元素污染的样本数据进行正态分布检验.本文采用的是利用SPSS软件的 图来检验他们是否服从正态分布若样本数据基本近似于在一条直线上就服从正态分布,若不服从正态分布,则对试验数据进行对数转换,看是否服从正态分布.经过反复比较,本文八种重金属中,As直接符合正态分布,其它七种重金属Cd、Cr、Cu、Ni、Hg、Pb、Zn经对数转化后服从正态分布.它们的 PLOT图如下：As的正态检验图1Cd的正态检验图2Cr的正态检验 图3 Cu的正态检验 图4Hg的正态检验 图5Ni的正态检验 图6 Pb的正态检验 图7Zn的正态检验 图8二问题1模型的建立与求解首先,根据附件一的样本,我们Matlab编程可以由
10、此作出该城市城区的样本取样点图,并标出这些样本点所属的不能功能区如图9.为了标出重金属元素在该城区的空间分布,并在不同区域重金属的污染程度,现对这八种重金属的样本等分成5组,分别用不同的颜色标注它们的污染程度,具体如下.样本点所属不能功能区 图9As的含量范围1.6114.75ug/g,由As的空间分布结构图可以看出,As污染最严重的地方主要在地图的左下方的地区,特别是工业区.土壤重金属元素As含量空间分布图10Cd的含量范围为40977.359ng/g,由Cd的空间分布结构图可以看出如图11,Cd污染最严重的地方主要分布在地图的左下图,特别是工业区和生活区,说明Cd的含量高主要由人类活动造成
11、的和工厂产生的废弃物.土壤重金属元素Cd含量空间分布图11 Cr含量的范围在15.32263.515ug/g,由Cr的空间分布结构图可以看出如图12,Cr污染最严重的地方主要分布在左下图,特别是生活区,说明城市市民的生活、大气降尘对重金属Cr的含量影响比较大.土壤重金属元素Cr含量空间分布图12重金属Cu 的分布范围为2.29543.762ug/g,由Cu的空间分布结构图可以看出如图13,Cu污染最严重的地方主要分布在地图的左下图,特别是工业区.说明Cu主要来源是工厂的含有重金属Cu的废弃物排放,而在山区和公园绿地区Cu的含量均在自然背景值之中.土壤重金属元素Cu含量空间分布图13Hg的含量范
12、围为8.575188.33ng/g,由Hg的空间分布结构图可以看出如图14,Hg污染最严重的地方主要在地图的左方,特别是工业区,说明Hg含量较高值区域主要是工厂排放污染物造成的. 土壤重金属元素Hg含量空间分布图14重金属Ni的分布范围为4.2747.08ug/g,由Ni的空间分布结构图可以看出如图15.Ni污染最严重的地方是主要在地图中左下方,特别是工业区.土壤重金属元素Ni含量空间分布图15Pb的含量范围为19.68211.915ug/g,由Pb的空间分布结构图可以看出如图16.Pb污染最严重的地方主要分布在地图的左下方,特别是工业区和生活区.Pb污染主要是人类活动与工厂排放污染物造成的.
13、土壤重金属元素Pb含量空间分布图16Zn的含量范围为 32.861218.9ug/g,由Zn的空间分布结构图可以看出如图17,Pb污染最严重的地方主要分布在地图的左下方,特别是工业区和生活区.这些地方Zn的含量高主要是因人类活动与工厂排放废水废渣引起的. 土壤重金属元素Zn含量空间分布图17三问题2模型的建立与求解样本数据的统计分析为确定重金属污染的主要原因,我们对城市五种不同的功能区进行八种重金属污染指数的单因子方差分析,分析这五种功能区的污染程度是否有显著性差异,若差异很显著结果表明它们的差异都很显著,则主要看是那个因子即功能区影响最大,从而确定污染的原因.由表18-1,F=12.9115
14、远大于临界值2.400402,这表明不同功能区域污染的程度有显著性差异.再由表18-2可以看出,工业区的As平均污染最高,然后是生活区和公园绿地.这说明As的污染主要是因为工业区工厂废水废渣的排放和居民日常生活所引起的污染.差异源SSdfMSFP-valueF crit组间261.2779465.3194912.91159.6E-102.400402组内1588.5313145.059016总计1849.809318As元素单因素方差分析表18-1组观测数求和平均方差生活区44275.96.2704554.623349工业区36250.29876.95274311.49077山区66266.9
15、14.0440913.237455交通区138761.52815.518324.65617公园绿地35219.236.2637144.094718As元素单因素方差分析表18-2由表19-1,F=22.94427远大于临界值2.400402,说明不同功能区域污染的程度很显著.再由表19-2可以看出,工业区的Cd平均污染最高,然后是生活区和交通区.这说明Cd的污染主要是因为工业区工厂废水废渣的排放、居民日常生活与车辆尾气排放所引起的污染.差异源SSdfMSFP-valueF crit组间31.5398447.88496122.944271.2E-162.400402组内107.90833140.343657总计139.4482318Cd元素单因素方差分析表19-1组观测数求和平均
