资源描述
Click to edit Master title style,Click to edit Master text styles,Second level,Third level,Fourth level,Fifth level,*,城市表层土壤重金属污染分析,组员：孙英峰,高丽,彭小伟,指导老师：任红萍,目录,1.,采用地下水模拟系统（,GMS,）软件,对样点分布、地表高程等高线及各种重金属质量浓度等值线进行拟合，描述重金属污染的空间分布。并且用单因子污染指数和内梅罗综合指数不同角度分析各区域重金属污染程度。,2.,用,spss,对各区重金属进行,ward,聚类分析，并分析各重金属污染物的相关性，据此推断各区重金属污染的主要原因。,3.,通过合理假设，,我们建立重金属污染物在土壤中的传播模型为一阶线性齐次偏微分方程,，确定污染源位置。,4.,模型的评价及完善。,由地表高程图可知，该城市总体为东北方地势最高，沿西南方向，地势逐渐降低，西面地势较为平坦，中部存在小盆地，东部及东北为高山地。,重金属元素在城区的空间分布,采样点分布图,As,浓度等值线图,Cd,浓度等值线,Cr,浓度等值线,Cu,的浓度等值线,Hg,的浓度等值线,Ni,的浓度等值线,Pb,的浓度等值线,Zn,的浓度等值线,1.2,各,区域重金属的污染程度,1.2.1,单因子污染指数和综合污染指数的定义,单因子污染指数：,式中，,为土壤中污染物,i,的环境污染指数；,为污染物,i,的实测值；,为污染物,i,的背景值。,综合污染指数：,式中，,为土壤污染中单因子污染指数最大的值；,为土壤污染物中单因子污染指数的平均值。,生活区,As,、,Ni,为轻度污染，,Cd,、,Cr,、,Hg,、,Pb,为中度污染，,Cu,、,Zn,为重度污染。,工业区,单因子污染指数五个功能区内最高，只有,Cr,、,Ni,、,As,污染较轻，其他的都很严重。,山区,只有,Zn,是严重污染，其他都属于轻度污染。,交通区,主要是,Hg,污染，其单因子污染指数高达，,As,、,Cr,、,Ni,为轻度污染，,Cd,、,Pb,为中度污染，,Cu,、,Hg,、,Zn,为重度污染。,公园绿地,As,、,Cr,、,Ni,、,Pb,属于轻度污染，,Cd,、,Cu,、,Zn,为中度污染，,Hg,为严重污染。,1.2.2,单因子污染指数结合污染程度等级,单因子污染指数,生活区,工业区,山区,公园绿地区,As,1.7418,2.0143,1.1234,1.5856,Cd,2.2305,3.0239,1.1717,2.7693,Cr,2.2264,1.7229,1.2568,1.8727,Cu,3.7427,9.6618,1.3119,4.7133,Hg,2.6583,18.353,1.1702,12.7664,Ni,1.4912,1.6107,1.2564,1.4323,Pb,2.2292,3.0013,1.1792,2.0495,Zn,3.4349,4.0279,4.0279,3.5196,等级,污染指数,污染程度,污染水平,1,P0.7,安全,清洁,2,0.7P1.0,警戒级,尚清洁,3,1.0P2.0,轻污染,土壤污染超过其背景值，,轻污染,4,2.0P3.0,重污染,土壤污染相当严重,单因子污染指数和综合污染指数,总体上，各区的污染物排序为：,生活区：,CuZn HgPbCdCrAsNi,工业区：,HgCuZnPbCrCdAsNi,山区：,Zn CuNiHgCrPb CdAs,主干道路,HgCuZnCrCdNiAsPb,公园绿地区：,HgZnCu CdPbAsCrNi,由综合污染指数指标中可以看出，总体来说，五个区的污染综合程度依次为：工业区,交通区,生活区,公园绿化区,山区。,2.,重金属污染的主要原因,公园绿地区及生活区,ward,聚类分析,公园绿地,区,ward,聚类法树状图,生活区,ward,聚类法树状图,由,上两图,可知，公园绿地区,Cd,，,Zn,与,Pb,有很大的相关性，,Cu,，,Hg,与,Cr,有很大的相关性，,Ni,与其它重金属的相关性不是很明显。与生活区相比较，这两个区域污染有很大的相似性,。由于公园绿地区与生活区相邻，生活区排放的各种生活污水以及生活垃圾中含有大量重金属。,交通区,ward,聚类分析,由,上图,可知交通区重金属污染中,Cr,，,Ni,与,Cu,有很大的相似性，,Pb,与,Cd,关联性很大，,Hg,与其它金属的污染关系不是很大。研究证实,Pb,和,Cd,在来源上关系较密切，,Pb,主要来源于汽车尾气排放，而汽车轮胎磨损可以释放,Cd,，,Cu,可以来源于汽油和车体的磨损，从而可以推断出其可能主要受汽车交通运输活动的影响，其次,Hg,污染不仅与燃煤污染有关，在很大的程度上还与汽车尾气的排放有关,。,工业区,ward,聚类分析,通过,上图,知,Zn,，,Ni,，,Cd,和,Cr,污染可能是由同一种污染源造成的，有研究表明，厂矿企业的三废排放中含有大量的,Zn,和,Cr,，由此我们可以认为工业区的工业排放是重金属,Zn,和,Cr,的主要原因。,山区,ward,聚类分析,由上图知，,Cu,，,Hg,和,Zn,的相关性比较大，而,Cd,和,Pb,有较大的相关性。,3,.,重金属污染物的传播模型,地表水的流动作用是重金属在土壤表层中传播的主要原因，故假设重金属传播以对流形式为主，海拔因素在建立方程时,先不,考虑，而是在确定传播方向时候考虑，将问题简化为二维，二维对流方程为：,其中，,为某重金属的的质量浓度，,为渗流速度沿两个坐标轴的分量。设地下水渗流速率是常数，变化的只是方向，输移方向沿着该点的海拔梯度方向的反方向。经查阅资料一般可设渗流速率为,设,.,把采样时刻设为,，将问题一所得的浓度空间分布函数,作为定解条件：,这时，传播模型为一阶线性齐次偏微分方程的初值问题：,易知方程的,通,解为,令 ，得 ，,即为初始时刻的重金属浓度的分布，此函数的极大值点即为污染源的位置。,由于本问题所给数据中没有涉及采样时间，我们设采样时刻为污染物排放后的半年，即,。,由问题一所得的每种重金属的分布图可知采样时刻浓度的最大值点，结合地表等高线图，查找每个极大值点附近等高线的分布，估算出该点的梯度方向，污染物在地表土壤的传播方向是从高海拔向低海拔，和梯度反向，我们可以由此确定,a,b,，比如由图,3,知，,As,的浓度有一个极大值点为（,4500,7500,），由图,1,估算其梯度的方向余弦为,则可设，,故这个污染源为,(4507.5,7495.7),，同理可以求出,As,的另外两个污染源的位置为,(12506,2994), (18500,10491.4),。,类似可得其它七种重金属的污染源,.,4.,模型的评价及完善,模型的优缺点,我们建立传播模型时，认为污染物传播是借助于地表水的渗透，所以引用对流方程，简化渗流速率为常数，渗流方向和该点的地表高程,梯度,有关，所以模型合理且精确，以第一题求出的浓度拟合函数为定解条件，解偏微分方程的初值问题即可确定污染源的位置。但由于建模时，,假设污染物是在一个时刻以点源的形式排放的；,未考虑重金属在土壤中的吸附、转化等化学过程；,未区分城市功能区、河流等诸多因素；,未考虑重金属在大气中的传播。,由计算结果来看，计算出的污染源的位置和采样得到的浓度最大值点相距,较近,，,其,原因是没有考虑河流、生活污水等承载重金属污染物的流体的流动所引起的迁移，且未考虑在交通区汽车尾气排放及汽车轮胎磨损产生大量,含重金属的,有害气体和粉尘经自然沉降和雨淋进入土壤，所以模型还有待完善。,4.2,如何更好地研究城市地质环境的演变模式,要准确地分析该城区的地质演变模式，首先数据需要在一定时间间隔做一次采样，也即得到动态的数据，观察同一采样点在不同时刻重金属浓度的变化情况；其次需要考虑降雨量数据、水文地质条件、地下水埋藏和分布规律、土壤中的含水量等数据。比如：,各,区土壤中含水量应该是有区别的，而含水量的多少影响渗流速率，我们可以根据疑似污染源所处的区域改变渗流速率，如,Cd,的两个污染源,(2006,，,2994), (21500,11491.5),，前者属于工业区，后者属于山区，若根据区域调整渗流速率，可以认为后者的位置为,(21500,11495),更合理。,可以根据不同时刻的采样数据，分别求解污染源的位置，观察污染源的位置移动轨迹，推测其精确位置。,待完善,谢谢各位评委！,
展开阅读全文