浙江省数学建模竞赛论文A题 一等奖Word文档格式.docx
《浙江省数学建模竞赛论文A题 一等奖Word文档格式.docx》由会员分享，可在线阅读，更多相关《浙江省数学建模竞赛论文A题 一等奖Word文档格式.docx（62页珍藏版）》请在冰豆网上搜索。
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进行细分：
Pi<
=1
1<
=3
3<
=5
Pi>
5
未污染
轻度污染
中度污染
重污染
所有划分结束后，利用SPSS软件，绘制出x，y，海拔与各元素关系的表格，即可从表中了解该城区的空间分布。
As元素与Cd元素的空间分布：
p（As（μg/g））
p（Cd（ng/g））
As1
As2
As3
As4
Cd1
Cd2
Cd3
Cd4
计数
x1
y1
z1
11
53
15
37
6
y2
17
16
2
1
3
23
9
y3
4
x2
18
13
21
19
39
z2
z3
x3
海拔
25
24
10
8
12
Cr元素与Cu元素的空间分布：
p（Cr（μg/g））
p（Cu（μg/g））
Cr1
Cr2
Cr3
Cr4
Cu1
Cu2
Cu3
Cu4
27
31
32
20
34
7
41
40
22
Hg元素与Ni元素的空间分布：
p（Hg（ng/g））
p（Ni（μg/g））
Hg1
Hg2
Hg3
Hg4
Ni1
Ni2
Ni3
Ni4
29
44
14
36
47
30
Pb元素与Zn元素的空间分布：
p（Pb（μg/g））
p（Zn（μg/g））
Pb1
Pb2
Pb3
Pb4
Zn1
Zn2
Zn3
Zn4
26
28
2）为了得到该城区不同区域重金属的污染程度，我们先用Excel对功能区进行筛选，对不同区域的重金属情况进行单独的分析。
我们在上文中已求出各元素的污染指数，用同样的方法可求出各区内各元素的污染指数：
As
Cd
Cr
Cu
Hg
Ni
Pb
Zn
1.021
1.556
1.077
2.308
4.088
0.860
1.389
2.017
1.059
1.642
1.141
2.969
6.647
0.888
1.501
2.246
0.987
1.378
0.845
1.666
4.059
0.812
1.214
1.400
1.049
1.594
1.092
2.701
5.892
0.866
1.437
2.078
1.037
1.553
1.073
2.259
3.977
1.374
1.981
以下表所示的评价标准对各区各元素的污染程度进行分析。
污染分析结果为：
用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（砷）：
砷是类金属元素，不是重金属。
但从它的环境污染效应来看，常把它作为重金属来研究。
土壤中砷的污染主要来自化工、冶金、炼焦、火力发电、造纸、玻璃、皮革及电子等工业排放的“三废”、冶金与化学工业、含砷农药的使用。
砷主要以正三价和正五价存在于土壤环境中。
其存在形式可分为水溶性砷，吸附态砷和难溶性砷。
由于土壤中砷主要以非水溶性形式
升级会员 