搜索
1城市表层土壤重金属污染分析摘要本文基于从某城区划分出的功能区：生活区、工业区、山区、主干道路区、公园绿地区中采集的数据，建立合理的模型，分析了该城市重金属的污染情况，并确定出污染源的位置。针对问题一，用Surfer软件对浓度数据进行克里格插值得出各类重金属元素的空间分布图。在分析城区重金属污染程度时，首先利用Muller指数法结合污染等级的打分，确定出同一功能区内不同元素的污染程度以及不同功能区内同一元素的污染程度。然后基于熵权法确定的各类元素对污染影响的权重，建立多目标模糊综合评价模型，得出各类功能区重金属污染的总体程度：山区无污染，其他地区轻度-中等污染。按模糊评价值得到的污染程度排序为：工业区主干道路区生活区公园绿地区山区。针对问题二，计算不同功能区内8种重金属元素的平均Muller指数值，结合第一问中得出的各类元素对污染影响的权重，确定出各个功能区中的主要污染元素：生活区(Cd、Zn)；工业区(Hg)；山区(Cd)；主干道路区(Hg、Cd)；公园绿地区(Hg、Cd)。然后运用变量聚类法按照重金属元素的污染相似性进行归类。再结合各功能区自身的特点，分析产生重金属污染的主要原因。针对问题三，基于重金属的来源和传播途径进行分析，得出污染的传播特征。综合考虑大气扩散建立高斯扩散方程，土壤迁移建立一维对流弥散方程，并计入海拔对土壤中对流速度的影响，得出大气～土壤传播模型。求解时本文将某个采样点假定为污染源，以其周围距离最近的20个采样点作为传播范围，由模型求解值与实际测量值得出平均相对误差，小于所设阈值则将其作为污染源，依此方法在所有采样点中进行污染源的搜寻。搜寻得到As、Cd、Cr、Cu、Hg、Ni、Pb、Zn污染源的个数分别为：7、6、3、2、3、2、3、4，并得出所有污染源的位置。例如Cu的所有污染源的位置为：(2427,3971)、(3299,6018)。针对问题四，结合实际污染物的传播特征，剖析所建模型的优缺点。基于模型没有考虑水流影响、土壤中纵向深度的传播、土壤吸附作用的主要缺点，收集每月不同地面位置不同土壤深度下各类污染物的浓度、每月平均降水量、土壤相关参数等信息，对模型进行完善改进。本文最后对各功能区污染程度的评价结果和确定污染源时的误差进行检验，剖析了前两问所建模型的优缺点，并就土壤重金属污染防治问题给出具体建议。关键词：克里格插值Muller指数熵权法多目标模糊综合评价变量聚类高斯扩散一维对流弥散2一、问题重述随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、……、5类区，不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值(标准值)。附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。现要求通过数学建模来完成以下任务：(1)给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。(2)通过数据分析，说明重金属污染的主要原因。(3)分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4)分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二、问题分析2.1问题背景分析2.1.1重金属污染物的主要来源及危害随着城市化和工业化进程的加快，城市工业排污量急剧增加，大量重金属污染排向了物环境中，重金属污染日趋严重。在一定条件下，某些重金属(例如汞)还能在某些微生物的作用下转化为毒性更大的有机物质。另外，有毒重金属可以长期停留与积累在环境中，通过食物链逐级富集，最终进入人体，甚至通过遗传或母乳使婴儿受害，主要表现为富集在人体某些器官内形成慢性中毒。例如日本1953—1956年的水俣病(Hg污染)、1955—1972年的骨痛病(Cd污染)以及l961年四日市的哮喘病(SO2和重金属粉尘复合污染)等等事件[1]。因此，重金属污染3问题成为学者们研究的一个焦点。土壤重金属污染是重金属污染的一个分支。土壤重金属污染是指由于人类活动，土壤中的微量有害元在土壤中的含量超过背景值，过量沉积而引起的含量过高，统称为土壤重金属污染。污染土壤的重金属主要包括汞、镉、铅、铬和类金属砷等生物毒性显著的元素，以及有一定毒性的锌、铜、镍等元素。污染主要来自农药、废水、污泥和大气沉降等[2]。2.1.2城市功能区城市功能区是是按照城市各地区功能进行划分的区域。城市功能区是实现经济社会各类职能的重要空间载体。其主要特点为：各类要素高度聚类、辐射带动效应明显、经济社会效益显著、具有明显的“城市名片”效应。按土地使用类型，可将城市划分为：生活区、工业区、山区、主干道路区、公园绿地区。住宅区分布最为广泛，是城市最基本的功能区，主要为人们提供便利的居住和生活条件；工业区一般分布在交通便利的地区的边缘，是城市发展的重要功能区；山区是城市地势较高的地方，人烟相对稀少，可能分布一些果园或者牧场；主干道路区是城市是城市最为繁华的地区，交通便利，人流量大，服务产业分布较多。公园绿地区主要用来进行城市绿化，改善空气质量，为居民提供休息娱乐的场所，一般分布在住宅区密集的地域。2.2问题一的分析问题一要求给出该城市八种污染物在空间的分布状况，之后分析城区内不同区域内重金属的污染程度。其中8种重金属元素为：As(砷)、Cd(镉)、Cr(铬)、Cu(铜)、Hg(汞)、Ni(镍)、Pb(铅)、Zn(锌)。首先是如何给出直观的污染物的空间分布。城市土壤重金属污染分布存在显著空间差异，由于不同的土地利用类型和频繁的人类干扰，同一城市的土壤重金属含量存在显著的空间差异，即使在较短距离内重金属含量也会完全不同。题目中给的是各个采样点的离散的数据，本文考虑画出污染物分布的三维连续图，来表现空间分布。在分析不同区域内重金属的污染程度时，考虑先利用已经研究成熟的评价重金属污染物浓度程度的评价方法，例如单因子指数方法，污染负荷指数法等，对5个功能区8种污染物分别进行分析，判断同一功能区下不同重金属污染物的污染程度，和不同功能区下同一种元素的污染程度。然后建立综合的评价模型，整合8种污染物在某功能区的总体污染情况，通过污染等级对不同功能区的污染程度进行量化。鉴于评价方法种类较多，在选择评价方法时，可以通过对比分析，选出最为合适可行的。2.3问题二的分析问题二要求结合数据的分析，说明重金属污染的主要原因。题目中考虑的重金属有8类，这8类元素对总体重金属污染的影响各不相同。可能某类元素的污染程度较高，但是综合到整个区域污染程度的分析，其所占的比重较小。因此本文考虑计算出各类元素对不同功能区重金属污染影响的贡献值，来确定造成不同区域污染的主要元素。根据理论知识可知，重金属化学性质和物理性质有些较为相似，有些污染来4源可能相同，即污染程度可能相似，可能产生交叉污染，而且8类重金属中按照毒性也可以分为两类：毒性显著（汞、镉、铅、铬、砷），毒性一般（铜、镍、锌）。例如镉和锌是同族元素，往往与锌、铜、铅等共生，在冶炼铜、锌及镀铬工厂中均有相当量的镉排出，造成环境污染。因此本文考虑对每个功能区将8种污染物基于污染浓度进行聚类，将污染物浓度变化趋势相近的元素聚为一类，即认为该类元素具有污染相似性。针对各类功能区的主要污染元素及元素间的污染相似性，结合实际情况分析造成各区域污染的主要原因。2.4问题三的分析题目要求分析重金属污染物的传播特征，基于分析建立传播模型，确定出城市中重金属污染源的位置。结合实际情况可以知道，污染物的传播过程是相当复杂的，传播到土壤中的重金属的来源各有差异，通过查询资料可以得知，污染物主要通过三种方式：大气沉降、土壤扩散、水流迁移进入土壤中。这三种方式都会导致城市产生严重的重金属污染。因为城市的水流分布是一个相当复杂错综的网络，必须结合某一城市的具体情况进行分析，因此本题实现的难度很大。基于这种情况，本文准备只考虑大气传播与土壤迁移的过程对总金属元素浓度的影响，建立以污染源为中心的大气-土壤传播模型。最后结合具体的模型，设计算法，确定出污染源的位置。2.5问题四的分析问题要求分析所建的传播模型的优缺点，就更好的研究地质环境的演变模式，给出需要另外收集的信息，并重新建立改善后的模型。本文针对第三问中所建立的大气-土壤传播模型结合实际生活中污染物的传播特征进行透彻分析，可以得出模型的优缺点。且在重新建立模型方面，本文打算主要针对模型没有考虑或者由于题目信息限制而简化的因素，进行数据收集，然后重新给出污染物传播的方程。三、模型假设1、假设采样点覆盖整个城区，数据能基本反映城市重金属污染情况；2、假设采样点得到的数据是污染物扩散达到稳定后的值；3、假设污染物不是无限扩散的，有一定扩散范围；4、假设各种重金属元素之间的传播不交叉影响；5、假设土壤迁移中分子扩散的作用可以小到忽略不计；6、假设污染源在采样点中进行确定。5四、符号说明符号说明jC采集点j类重金属的浓度jBj类污染物的背景值lgeo地累积指数值ijx第i个采集点点第j类重金属的lgeo值ijr对ijx作标准化处理后的值ijp第j类重金属下第i个观测点占污染程度的比重jej类金属的熵值jgj类金属的差异系数jwj类金属对污染的权重,Cxy土壤中重金属污染物的浓度1,Cxy通过大气传播的重金属的浓度2,Cxy通过土壤扩散的重金属的浓度Q源强u平均风速,yz大气中的扩散系数v土壤中孔隙水流速度D弥散系数x、y、z地面坐标和土壤深度坐标h海拔6pt降雨量函数dR土壤吸附作用的组织系数五、模型的建立与求解5.1问题一模型的建立与求解5.1.1模型准备：数据的处理题目附件中给出的是8种重金属污染物在319个采样点的浓度值，数据量较大，不利于问题的分析，因此对数据进行做如下处理。对五类功能区的对应的数据进行筛选、分类汇总，得到如下8种重金属含量在不同功能区的统计特征：表1：生活区土壤重金属含量的统计特征生活区重金属污染物浓度变异系数最大值最小值平均值标准差As(μg/g)11.452.346.272.150.342903Cd(ng/g)1044.586.8289.96183.680.633467Cr(μg/g)744.4618.4669.02107.891.56317Cu(μg/g)248.859.7349.447.160.954656Hg(ng/g)5501293.04102.91.105976Ni(μg/g)32.88.8918.345.660.308615Pb(μg/g)472.4824.4369.1172.321.046448Zn(μg/g)2893.4743.37237.01443.641.87182表2：工业区土壤重金属含量的统计特征工业区重金属污染物浓度变异系数最大值最小值平均值标准差As(μg/g)21.871.617.254.240.584828Cd(ng/g)1092.9114.5393.11237.580.60436Cr(μg/g)285.5815.453.41440.823816Cu(μg/g)2528.4812.7127.54414.933.253332Hg(ng/g)1350011.71642.362244.083.493493Ni(μg/g)41.74.2719.818.370.