浙江省数学建模竞赛论文A题 一等奖.docx
《浙江省数学建模竞赛论文A题 一等奖.docx》由会员分享,可在线阅读,更多相关《浙江省数学建模竞赛论文A题 一等奖.docx(62页珍藏版)》请在冰豆网上搜索。
浙江省数学建模竞赛论文A题一等奖
城市表层土壤重金属污染分析
摘要
工业污染物质和生活垃圾的大量输入以及长期使用化肥、农药、污水灌溉等,容易导致重金属元素的蓄积,进而造成土壤重金属的污染。
土壤重金属的污染影响土壤的生态功能,影响植物的生长、农产品的质量,并且通过食物链传递和土壤颗粒物的直接吸入而影响人体健康;土壤重金属污染具有隐蔽性、长期性和不可逆性的特点,越来越受关注。
对于本问题我们主要运用SPSS软件对数据进行处理,从而建立经验模型,研究该城市不同功能区的表层土壤重金属含量特征,计算8种重金属元素的单项污染指数,并成功解决了问题。
针对问题一,我们运用SPSS软件,定制关于x,y,海拔,各元素单项污染指数报表,通过观察和分析报表,得出重金属元素的空间分布;对于污染程度的分析,我们绘制以功能区为x轴、各重金属元素的平均污染指数为y轴的分组条形图,并且对不同功能区的各元素的污染指数制定表格,从图形和表格中准确分析不同区域重金属的污染程度。
针对问题二,我们通过网上大量搜寻8种重金属元素的产生原因、特性、存在状态,及土壤母质的相关属性,绘制了8种重金属元素含量在不同功能区的饼图,说明重金属污染的主要原因。
针对问题三,我们结合问题一和问题二,对各重金属元素的迁移转化进行研究,综合得到重金属污染物的传播特征;建立经验模型,对于x,y和各元素的综合污染指数定制报表,统计出各元素污染情况严重的区域,进而得出各元素的污染源。
本模型优点在于直观和可操作性强。
针对问题四,我们的模型采用经验模型,适用于大范围内长时段重金属转移,对资料依赖性较强,无法阐明重金属迁移与污染源的复杂关系。
因此,需要我们收集该地区土壤的酸碱度PH、降雨量、土壤母质这些相关信息,通过分析,预测模型,从而解决问题。
关键词:
SPSS软件平均污染指数经验模型预测模型
问题重述
随着城市经济的快速发展和城市人口的不断增加,人类活动对城市环境质量的影响日显突出。
对城市土壤地质环境异常的查证,以及如何应用查证获得的海量数据资料开展城市环境质量评价,研究人类活动影响下城市地质环境的演变模式,日益成为人们关注的焦点。
按照功能划分,城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等,分别记为1类区、2类区、……、5类区,不同的区域环境受人类活动影响的程度不同。
现对某城市城区土壤地质环境进行调查。
为此,将所考察的城区划分为间距1公里左右的网格子区域,按照每平方公里1个采样点对表层土(0~10厘米深度)进行取样、编号,并用GPS记录采样点的位置。
应用专门仪器测试分析,获得了每个样本所含的多种化学元素的浓度数据。
另一方面,按照2公里的间距在那些远离人群及工业活动的自然区取样,将其作为该城区表层土壤中元素的背景值。
附件1列出了采样点的位置、海拔高度及其所属功能区等信息,附件2列出了8种主要重金属元素在采样点处的浓度,附件3列出了8种主要重金属元素的背景值。
现要求你们通过数学建模来完成以下任务:
(1)给出8种主要重金属元素在该城区的空间分布,并分析该城区内不同区域重金属的污染程度。
(2)通过数据分析,说明重金属污染的主要原因。
(3)分析重金属污染物的传播特征,由此建立模型,确定污染源的位置。
(4)分析你所建立模型的优缺点,为更好地研究城市地质环境的演变模式,还应收集什么信息?
有了这些信息,如何建立模型解决问题?
模型假设
1、假设只有题中所给的8种元素对该城区造成重金属污染;
2、假设每个样本都能很好地反映该区域的实际情况。
符号说明
i:
i=1,2,3,…,8分别表示As,Cd,Cr,Cu,Hg,Ni,Pb,Zn这8种重金属元素
Pi:
土壤中污染物i的环境质量指数,i=1,2,3,…,8
Ci:
污染物i的实测值(ug/g)
Si:
污染物i的评价标准(ug/g)
Asj:
元素As的污染程度,j=1,2,3,4
Cdj:
元素Cd的污染程度,j=1,2,3,4
Crj:
元素Cr的污染程度,j=1,2,3,4
Cuj:
元素Cu的污染程度,j=1,2,3,4
Hgj:
元素Hg的污染程度,j=1,2,3,4
Nij:
元素Ni的污染程度,j=1,2,3,4
Pbj:
元素Pb的污染程度,j=1,2,3,4
Znj:
元素Zn的污染程度,j=1,2,3,4
问题分析
2-1:
问题1分析
本问主要先研究8种主要重金属元素在该城区的空间分布,再分析该城区不同区域重金属的污染程度。
我们对题目所给的数据进行分析,运用spss软件对x,y,海拔与八种重金属元素中的As元素进行相关性分析,得出x,y,海拔与As元素的相关指数小于1,即x,y,海拔只表示取样点的空间位置,与元素浓度没有一定的关系。
(以As为例,若该元素与x,y,海拔具有相关性则可对其他金属元素进行相关性分析;若无相关性则可说明其他元素也与x,y,海拔之间没有相关性)相关性分析结果见附录1。
因此,为了得到这8种重金属元素在该城区的空间分布,我们分别将x,y,海拔等分成三组,即可得到27个空间区域,再根据单因子指数法[1]计算出各重金属元素的污染指数Pi(其中Pi<=1表示未污染,Pi>1表示污染,然后对污染的部分进行污染程度定义,分为轻度污染,中度污染,重污染)统计出各区域内各元素达到不同污染程度的个数,制成表格,即可得到各元素在该城区的空间分布。
为了得到该城区不同区域重金属的污染程度,我们先用Excel对功能区进行筛选,对不同区域的重金属污染情况进行单独的分析。
在上文我们已求出各元素的污染指数Pi,我们只要将数据进行整理,得到各区各元素的污染指数表,对表中数据进行分析即可得到各区重金属的污染程度。
我们再用Excel绘制出不同重金属元素分布情况的分组条形图,从图中观察不同功能区各重金属元素的分布情况,并对该区的污染程度进行分析,将前后两种分析结果进行对比,便可得到更准确的结论。
2-2:
问题2分析
通过数据分析对重金属污染的主要原因进行说明,数据处理后的图形可以形象地表示出污染情况,然后对污染情况进行分析得到重金属污染的主要原因。
我们通过查找各类资料可以知道造成这8种重金属污染的原因,再从问题1中得到的“各功能区的重金属污染情况图”中了解到各个功能区中哪些重金属的污染较为严重。
我们根据各区的污染情况,绘制出各区污染情况的饼图。
通过分析图形,综合造成这几种污染较为严重的元素的原因,得出该地区重金属污染的主要原因。
2-3:
问题3分析
本问主要分析重金属污染物的传播特征,根据传播特征建立模型,从而确定污染源的位置。
我们从重金属污染物的传播特征入手,通过上网搜索资料等途径找到各元素的传播特征及其化学性质。
对各元素的传播特征进行分析可以发现,大部分重金属污染物为难溶性化合物且迁移能力较弱。
因此,我们可以初步断定重金属浓度较大且较为集中的区域为污染源所在区域。
在问题1的分析中,我们已经得到各元素在该城区的空间分布,且其用计数表示污染程度。
在本问中,我们可以采用相似的方法,由于我们寻找的污染源用x,y即可表示出来,所以只需将x,y,进行等分处理,且本题中与题1相比数据的划分可以再细致些,这样得到的区域范围可以更小更准确。
为了减小计算量,我们将x等分成10组,而将y等分成6组。
分组结束后用SPSS统计各元素在不同区域中超标的个数。
模型建立与求解
3-1:
问题1的模型建立与求解
1)通过分析我们可以先将x,y,海拔这些连续变量离散化,即将x,y,海拔均等分成三组:
x(m)
y(m)
海拔(m)
第一组(x1)
[0,9552]
[0,6150]
[0,103]
第二组(x2)
[9553,19105]
[6151,12301]
[104,207]
第三组(x3)
[19106,28658]
[12302,18452]
[208,311]
根据单因子指数法计算出各重金属元素的污染指数Pi(各元素Pi的计算结果见附录2),Pi<=1表示未污染,Pi>1表示污染,我们可将Pi进行细分:
Pi<=1
13Pi>5
未污染
轻度污染
中度污染
重污染
所有划分结束后,利用SPSS软件,绘制出x,y,海拔与各元素关系的表格,即可从表中了解该城区的空间分布。
As元素与Cd元素的空间分布:
p(As(μg/g))
p(Cd(ng/g))
As1
As2
As3
As4
Cd1
Cd2
Cd3
Cd4
计数
计数
计数
计数
计数
计数
计数
计数
x(m)
x1
y(m)
y1
海拔(m)
z1
11
53
0
0
15
37
6
6
y2
海拔(m)
z1
17
16
2
1
3
23
9
1
y3
海拔(m)
z1
1
6
0
0
2
4
1
0
x2
y(m)
y1
海拔(m)
z1
23
18
1
0
13
23
6
0
y2
海拔(m)
z1
37
21
0
1
19
39
1
0
z2
2
0
0
0
1
1
0
0
y3
海拔(m)
z1
9
21
0
0
11
19
0
0
z2
4
0
0
0
2
2
0
0
z3
1
0
0
0
0
1
0
0
x3
y(m)
y1
海拔(m)
z1
6
3
0
0
6
3
0
0
y2
海拔
z1
25
11
0
0
24
10
1
1
z2
3
1
0
0
3
1
0
0
y3
海拔(m)
z1
8
3
0
0
5
6
0
0
z2
12
0
0
0
9
3
0
0
z3
2
0
0
0
1
1
0
0
Cr元素与Cu元素的空间分布:
p(Cr(μg/g))
p(Cu(μg/g))
Cr1
Cr2
Cr3
Cr4
Cu1
Cu2
Cu3
Cu4
计数
计数
计数
计数
计数
计数
计数
计数
x(m)
x1
y(m)
y1
海拔(m)
z1
27
31
2
4
6
32
15
11
y2
海拔(m)
z1
15
20
1
0
1
13
11
11
y3
海拔(m)
z1
4
3
0
0
1
5
1
0
x2
y(m)
y1
海拔(m)
z1
34
7
0
1
12
25
5
0
y2
海拔(m)
z1
41
18
0
0
11
40
6
2
z2
1
1
0
0
1
1
0
0
y3
海拔(m)
z1
22
8
0
0
9
18
2
1
z2
4
0
0
0
3
1
0
0
z3
1
0
0
0
1
0
0
0
x3
y(m)
y1
海拔(m)
z1
6
3
0
0
3
6
0
0
y2
海拔(m)
z1
32
3
1
0
22
12
1
1
z2
1
3
0
0
4
0
0
0
y3
海拔(m)
z1
10
1
0
0
8
3
0
0
z2
10
2
0
0
12
0
0
0
z3
2
0
0
0
2
0
0
0
Hg元素与Ni元素的空间分布:
p(Hg(ng/g))
p(Ni(μg/g))
Hg1
Hg2
Hg3
Hg4
Ni1
Ni2
Ni3
Ni4
计数
计数
计数
计数
计数
计数
计数
计数
x(m)
x1
y(m)
y1
海拔(m)
z1
9
29
13
13
44
19
0
1
y2
海拔(m)
z1
6
15
6
9
22
14
0
0
y3
海拔(m)
z1
2
3
1
1
4
3
0
0
x2
y(m)
y1
海拔(m)
z1
24
11
3
4
36
6
0
0
y2
海拔(m)
z1
36
18
2
3
47
12
0
0
z2
2
0
0
0
1
1
0
0
y3
海拔(m)
z1
21
8
0
1
24
6
0
0
z2
2
2
0
0
4
0
0
0
z3
0
1
0
0
1
0
0
0
x3
y(m)
y1
海拔(m)
z1
6
3
0
0
7
2
0
0
y2
海拔(m)
z1
30
4
1
1
32
3
1
0
z2
4
0
0
0
1
3
0
0
y3
海拔(m)
z1
10
1
0
0
10
1
0
0
z2
11
1
0
0
12
0
0
0
z3
1
1
0
0
2
0
0
0
Pb元素与Zn元素的空间分布:
p(Pb(μg/g))
p(Zn(μg/g))
Pb1
Pb2
Pb3
Pb4
Zn1
Zn2
Zn3
Zn4
计数
计数
计数
计数
计数
计数
计数
计数
x(m)
x1
y(m)
y1
海拔(m)
z1
16
37
6
5
13
29
10
12
y2
海拔(m)
z1
2
31
3
0
2
20
10
4
y3
海拔(m)
z1
2
5
0
0
2
4
1
0
x2
y(m)
y1
海拔(m)
z1
23
19
0
0
17
21
3
1
y2
海拔(m)
z1
27
30
2
0
29
26
2
2
z2
1
1
0
0
1
1
0
0
y3
海拔(m)
z1
16
14
0
0
13
16
1
0
z2
1
3
0
0
2
2
0
0
z3
1
0
0
0
1
0
0
0
x3
y(m)
y1
海拔(m)
z1
7
2
0
0
5
4
0
0
y2
海拔(m)
z1
26
10
0
0
28
8
0
0
z2
4
0
0
0
4
0
0
0
y3
海拔(m)
z1
6
5
0
0
8
3
0
0
z2
12
0
0
0
12
0
0
0
z3
1
1
0
0
2
0
0
0
2)为了得到该城区不同区域重金属的污染程度,我们先用Excel对功能区进行筛选,对不同区域的重金属情况进行单独的分析。
我们在上文中已求出各元素的污染指数,用同样的方法可求出各区内各元素的污染指数:
As
Cd
Cr
Cu
Hg
Ni
Pb
Zn
1
1.021
1.556
1.077
2.308
4.088
0.860
1.389
2.017
2
1.059
1.642
1.141
2.969
6.647
0.888
1.501
2.246
3
0.987
1.378
0.845
1.666
4.059
0.812
1.214
1.400
4
1.049
1.594
1.092
2.701
5.892
0.866
1.437
2.078
5
1.037
1.553
1.073
2.259
3.977
0.860
1.374
1.981
以下表所示的评价标准对各区各元素的污染程度进行分析。
Pi<=1
13Pi>5
未污染
轻度污染
中度污染
重污染
污染分析结果为:
As
Cd
Cr
Cu
Hg
Ni
Pb
Zn
1
轻度污染
轻度污染
轻度污染
轻度污染
中度污染
未污染
轻度污染
轻度污染
2
轻度污染
轻度污染
轻度污染
轻度污染
重污染
未污染
轻度污染
轻度污染
3
未污染
轻度污染
未污染
轻度污染
中度污染
未污染
轻度污染
轻度污染
4
轻度污染
轻度污染
轻度污染
轻度污染
重污染
未污染
轻度污染
轻度污染
5
轻度污染
轻度污染
轻度污染
轻度污染
中度污染
未污染
轻度污染
轻度污染
用Excel绘制出不同重金属元素分布情况的分组条形图,其中横坐标为功能区,纵坐标为各元素的污染指数。
从图中观察不同功能区各重金属元素的分布情况,并对该区的污染程度进行分析。
图1
由图1可知,城区1中Hg的污染最为严重,为中度污染;城区2中Hg为重污染;城区3中Hg为中度污染;城区4中Hg为重污染;城区5中Hg为中度污染,与表中数据相符,且其他元素的相似度也非常接近,故上述表格中对不同区域重金属的污染程度描述较为准确。
3-2:
问题2的模型建立与求解
由3-1中“各功能区的重金属污染情况图”可以对各功能区的污染情况进行分析,并找出该区重金属污染的主要原因。
分析各区重金属污染原因前先找出能造成各重金属污染的原因,通过上网搜索资料可以了解到造成各金属污染的原因([2][3])有如下方式:
As污染:
砷污染是指由砷或其化合物所引起的环境污染。
砷和含砷金属的开采冶炼,用砷或砷化合物作原料的玻璃、颜料、原药、纸张的生产以及煤的燃烧等过程,都可产生含砷废水、废气和废渣,对环境造成污染。
Cd污染:
镉污染主要来自铅锌矿以及有色金属冶炼、电镀和用镉化合物作原料或触媒的工厂。
镉对土壤的污染主要有气型和水型两种。
气型污染主要由含镉工业废气扩散并自然沉降,蓄集于工厂周围的土壤中,污染范围有的可达数公里。
水型污染主要是铅锌矿的选矿废水和有关工业(电镀、碱性电池等)废水排入地面水或渗入地下水引起的。
水体中镉的污染主要来自地表径流和工业废水。
Cr污染:
主要来源于劣质化妆品原料、皮革制剂、金属部件镀铬部分,工业颜料以及鞣革、橡胶和陶瓷原料等。
Cu污染:
指铜及其化合物在环境中所造成的污染。
主要污染来源是铜锌矿的开采和冶炼、金属加工、机械制造、钢铁生产等。
冶炼排放的烟尘是大气铜污染的主要来源。
Hg污染:
汞及其化合物属于剧毒物质,可在人体内蓄积。
主要来源于仪表厂、食盐电解、贵金属冶炼、化妆品、照明用灯、齿科材料、燃煤、水生生物等。
Ni污染:
是由镍及其化合物所引起的环境污染。
冶炼镍矿石及冶炼钢铁时,部分矿粉会随气流进入大气。
在焙烧过程中也有镍及其化合物排出,主要为不溶于水的硫化镍(NiS),氧化镍(NiO)、金属镍粉尘等,成为大气中的颗粒物。
镍可在土壤中富集,含镍的大气颗粒物沉降、含镍废水灌溉、动植物残体腐烂、岩石风化等都是土壤中镍的来源。
Pb污染:
主要来来源于各种油漆、涂料、蓄电池、冶炼、五金、机械、电镀、化妆品、染发剂、釉彩碗碟、餐具、燃煤、膨化食品、自来水管等。
Zn污染:
是指锌及化合物所引起的环境污染。
主要污染源有锌矿开采、冶炼加工、机械制造以及镀锌、仪器仪表、有机会合成和造纸等工业的排放。
汽车轮胎磨损以及煤燃烧产生的粉尘、烟尘中均含有锌及化合物,工业废水中锌常以锌的羟基络合物存在。
下面对各区受污染情况进行分析:
1区(生活区):
通过生活区受污染情况饼图可以了解到该区主要污染元素为Hg元素,除了Ni元素未污染外,其余各元素均已造成重金属污染。
生活垃圾中有很多含有重金属的危险废物,比如废旧电器、电路板、光管、电池等,这些物品如不集中处理而混在普通垃圾中排放就会带来重金属污染,还有一些化妆品、玻璃、油漆涂料、燃煤等均会造成重金属污染。
图2
2区(工业区):
通过工业区受污染情况饼图可以了解到该区主要污染元素为Cu元素,其他元素如Zn、Hg、Pb等元素的污染情况也比较严重。
由资料可知Cu元素主要污染来源是铜锌矿的开采和冶炼、金属加工、机械制造、钢铁生产等。
冶炼排放的烟尘是大气铜污染的主要来源,工业产业的废弃物,工业废水都含有很大比重的重金属污染物。
图3
3区(山区):
通过工业区受污染情况饼图可以了解到该区主要污染元素为Hg元素(Hg所占比例为56.45%,超过了一半),其他重金属元素所占比例则大致相同。
Hg污染主要来源于仪表厂、食盐电解、重金属冶炼、化妆品、照明用灯、齿科材料、燃煤、水生生物等,而该区只有Hg的比例最大且较为突出,因此排除仪表厂、食盐电解、重金属冶炼等容易造成其他重金属元素污染的原因。
因此,我们可猜想该区有垃圾处理的场地,且该地垃圾中的化妆品、照明用灯、齿科材料等为该地区重金属污染的主要原因。
图4
4区(交通区):
通过交通区受污染情况饼图可以了解到该区主要污染元素有Hg、Cu、Zn等元素。
结合交通区的实际情况与能造成这些元素污染的原因,我们可以初步认为交通区重金属污染的主要原因为机动车量排放的汽车废气、车轮在行驶中脱落的橡胶、汽油燃烧等。
图5
5区(公园绿地区):
通过公园绿地区受污染情况饼图可以了解到该区主要污染元素为Hg元素(Hg所占比例为52.95%,超过了一半),其他重金属元素所占比例则大致相同。
这种情况与山区有点相似,但公园绿地区不可能像山区那样设立垃圾场或其他工业制造厂。
但是,可以假设该绿化区由于人员来往频繁而产生更多的垃圾废弃物,最终导致重金属污染。
图6
3-3:
问题3的模型建立与求解
步骤1.通过搜索资料等方式,我们对重金属污染物的传播特征进行了整理。
各重金属污染物的传播特征[4][5]如下:
重金属一般不易随水淋滤,土壤微生物不能分解,但能吸附于土壤胶体、被土壤微生物和植物所吸收,通过食物链或其它方式转化为毒性更强的物质,对人体健康的危害严重,所以土壤中重金属的污染问题比较突出。
重金属在土壤中的迁移转化受金属的化学特性、土壤的物理特性、生物特性和环境条件等因素影响。
As(砷):
砷是类金属元素,不是重金属。
但从它的环境污染效应来看,常把它作为重金属来研究。
土壤中砷的污染主要来自化工、冶金、炼焦、火力发电、造纸、玻璃、皮革及电子等工业排放的“三废”、冶金与化学工业、含砷农药的使用。
砷主要以正三价和正五价存在于土壤环境中。
其存在形式可分为水溶性砷,吸附态砷和难溶性砷。
由于土壤中砷主要以非水溶性形式