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