全国大学生数学建模A题全国优秀论文北京理工大学.docx
- 文档编号:3801924
- 上传时间:2022-11-25
- 格式:DOCX
- 页数:63
- 大小:576.22KB
全国大学生数学建模A题全国优秀论文北京理工大学.docx
《全国大学生数学建模A题全国优秀论文北京理工大学.docx》由会员分享,可在线阅读,更多相关《全国大学生数学建模A题全国优秀论文北京理工大学.docx(63页珍藏版)》请在冰豆网上搜索。
全国大学生数学建模A题全国优秀论文北京理工大学
2011高教社杯全国大学生数学建模竞赛
承诺书
我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.
我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。
我们知道,抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。
如有违反竞赛规则的行为,我们将受到严肃处理。
我们参赛选择的题号是(从A/B/C/D中选择一项填写):
A
我们的参赛报名号为(如果赛区设置报名号的话):
BJ01259
所属学校(请填写完整的全名):
北京理工大学
参赛队员(打印并签名):
1.王山
2.戴琪
3.李永
指导教师或指导教师组负责人(打印并签名):
教练组
日期:
2011年9月12日
赛区评阅编号(由赛区组委会评阅前进行编号):
2011高教社杯全国大学生数学建模竞赛
编号专用页
赛区评阅编号(由赛区组委会评阅前进行编号):
赛区评阅记录(可供赛区评阅时使用):
评
阅
人
评
分
备
注
全国统一编号(由赛区组委会送交全国前编号):
全国评阅编号(由全国组委会评阅前进行编号):
城市表层土壤重金属污染分析
摘要
本文通过对城市表层土壤重金属污染综进行分析,建立了内梅罗评价模型、内因子分析模型确定重金属对城区的污染情况及抛曲面拟合模型确定污染源位置等一系列数学模型。
最后,又提出了对研究城市地质环境的演变模式的可行方法。
问题一:
建立重金属污染程度处理的内梅罗定量评价模型。
导入附件1、附件2中坐标位置,海拔高度及8种金属元素的浓度数据,用Surfer9.0软件依据源文件数据作图,可以分别得到城区区域拟合的区域图,即功能区平面图,及8种金属元素在各城区区域拟合的等高浓度线分布图。
从图中清晰的等高线分布,定性分析元素,同时用单因子污染指数
,内梅罗综合污染指数法
,将该城区内不同区域的重金属的污染程度量化分析,从直观的单因子指数与综合指数中,确定了该城区不同区域重金属污染的程度。
问题二:
源数据进行内梅罗模型数据处理后,得到评判区域污染的污染指标Pi及由Pi演化的图表,在分块图表中获得出污染元素的区域特性,并进一步利用因子分析模型确定各区域污染的主因子。
因子分析模型确定该城区污染原因主要在工业“三废”,汽车尾气排放及轮胎磨损。
问题三:
问题一、二反映出城区污染不仅有点源污染,还有面源污染,结合重金属污染物因雨水冲刷、大气沉降等因素浓度梯度由高到低迁移及海拔由高到低迁移的传播特征,建立抛曲面拟合模型,模型借用物体辐射抛射产生的效果,近似地模拟城区污染程度与地理位置间的关系,通过对最小二乘法修正,得
,最终确定有2个可供参考的主要污染源。
问题四:
对在前三个问题中建立的数学模型做了简要的分析。
并对更好地研究城市地质环境的演变模式提出应收集所需研究的元素在土壤中的浓度数据的年度记录,建立扩散模型,扩散方程为:
关键词:
内梅罗综合污染指数因子分析法三维曲面拟合扩散模型
问题一
1.1问题分析
问题一是要求分析该城区内不同区域重金属的污染程度,可以根据附件一中列出的采样点的位置、海拔高度以其所属功能区等信息可以用Surfer9.0画出该城区区域图的各分区大致位置,再根据附件二中8种主要重金属元素在采样点处的浓度,可以同样得到这8种重金属分别在该城区中的浓度分布情况。
再采用单因子污染指数和综合污染指数法对该城区内不同区域的重金属的污染程度进行评价,从而得出结论。
1.2模型假设
(1)假设取样点的数据能够较好的反映该地区的污染物的浓度;
(2)假设污染源的重金属的浓度在一定的时间内是稳定的。
1.3模型建立和求解
通过附件一中列出的采样点位置和海拔高度可以用Surfer9.0拟合出该城区分区的大致位置,如图1.1,图中不同的颜色对应着城区中不同的分区。
图1.1城区各功能区拟合图
通过附件一中采样点的位置坐标和附件二中采样点的各金属的浓度可以用数学软件Surfer拟合出8种金属元素分别对该城区的污染程度图,分别如图1.2~图1.9,图中等高线的梯度越大,分布越密集,则此种金属在此分布范围内就浓度越大。
图1.3Cd浓度分布图
图1.2As浓度分布图
图1.5Cu浓度分布图
图1.4Cr浓度分布图
图1.7Ni浓度分布图
图1.6Hg浓度分布图
图1.9Zu浓度分布图
布图
图1.8Pb浓度分布图
布图
从图1.2到图1.9只能够直观的大概反映出各种重金属元素对该城区的污染程度,但是还是无法定量的来反映出不同区域重金属的污染程度,因此需要采用单因子污染指数和综合污染指数法对该城区内重金属元素污染状况进行评价,其计算公式为:
式中:
为重金属
的单因子指数;
为重金属
的实测浓度。
当
时表示未受污染,当
时表示受到污染,且
越大污染越严重。
表1.1城区各种重金属元素平均值
重金属元素
As
Cd
Cr
Cu
Hg
Ni
Pb
Zn
Pi平均值
1.05
1.59
1.09
2.70
5.88
0.87
1.44
2.07
从表1.1中可以看出该城区内除了重金属Ni的Pi平均值小于1,其他各金属的Pi平均值都大于1.所以该城区内未受到重金属Ni的污染影响,而其他金属对该城区均有影响,其中重金属Hg对该城区的污染最严重。
各重金属元素对城区的污染程度大致为Hg>Cu>Zn>Cd>Pb>Cr>As>Ni。
要对该城区内不同区域重金属的污染程度的评价则需要用到综合污染指数法,其公式为:
式中:
为某地区的综合污染指数;
为重金属污染物中污染指数最大值,
为重金属污染物中污染指数平均值。
通过对附录中的重金属浓度和背景值的分析计算,可以得到一区到五区中各种重金属综合作用污染指数
,如表1.2所示。
表1.2城区各分区的各种重金属综合作用污染指数
生活区
工业区
山区
主干道路区
公园绿地区
P综
2.09
9.28
0.82
6.46
1.87
表1.3中国绿色食品发展中心推荐的分级标准进行的污染分级
等级
综合污染指数
污染程度
1
P<=0.7
安全
2
0.7
警戒线
3
1
轻污染
4
1.5
中污染
5
2
重污染
6
P>3
严重污染
通过对该城区5个分区的319个取样点,得到如表1.2的各城区的各种重金属综合作用污染指数。
对照着表1.3的土壤重金属分级标准,可以发现除山区位于警戒线附近,其他的各区均受到不同程度的污染,各区的污染程度依次为工业区>主干道路区>生活区>公园绿地区>山区。
工业区和主干道路区的污染程度远远高于其他各区,生活区和公园绿地区的污染程度大致相同,山区的污染程度最小,但是也位于被污染的边缘。
问题二
2.1问题分析
利用问题一对所提供污染元素的浓度进行无量纲处理,并取8种元素的调和平均数作为功能区粗评定综合指标。
由附录二表建立模型的样本空间,利用样本空间中Pi的调和平均数间的比重关系,确定各功能区占主导因素的污染元素,得到如下百分比图:
图表1
图表2
图表3
图表4
图表5
2.1.2结论与分析
从图表中放映出来的信息分析之后,可以得出生活区以Hg,Cu,Pb为主要污染元素,工业区以Hg,Cu为主要污染元素,山区各污染元素相对平等,主干道路区以Hg,Cu为主要污染元素,公园绿地区以Hg,As为主要污染元素。
虽然,图表信息体现了区域中起较重要的污染元素,但区域特征及污染原因分析中所得数据集中在Hg,Cu之中,由于调和平均数受极差影响使得所得结果不能让区域性及污染特点量化,故需对该模型进行优化处理。
2.2.1内因子模型建立与求解
设有n个样本,n个指标构成样本空间X,
,i=1,2,…,n;j=1,2,…m,并对原始数据
标准化,标准化公式为
,其中
为样本中第i个样本的第j指标值,而
和
分别为指标的均值和标准差;计算标准化数据的相关系数矩阵,进而消除不同变量的量纲的影响且标准操作数据不会影响变量的相关系数;求出相关系数的特征值和特征向量,进行正交变换,使用方差最大法,使得所得因子载荷两极分化,但不影响旋转后的因子的正交特性。
模型通过正交的方差最大旋转法是每个主因子至于最少个数的变量有相关系数,使得足够多的因子负荷均减小,并通过多元统计分析的计算结果,分析及确定主因子。
2.2.2内因子模型结论与分析
生活区主因子的求解:
由附录一表2-1生活区土壤重金属浓度得到模型建立的样本空间X,借助统计软件SPSS11.0计算得出如下表格:
As(μ/g)
Cd(ng/g)
Cr(μg/g)
Cu(μ/g)
Hg(ng/g)
Ni(μg/g)
Pb(g/g)
Zn(g/g)
As(μg/g)
1.000
Cd(ng/g)
0.381
1.000
Cr(μg/g)
0.238
0.349
1.000
Cu(μg/g)
0.531
0.499
0.376
1.000
Hg(ng/g)
0.293
0.397
0.150
0.198
1.000
Ni(μg/g)
0.605
0.283
0.527
0.434
0.211
1.000
Pb(μg/g)
0.450
0.802
0.416
0.502
0.340
0.300
1.000
Zn(μg/g)
-0.017
0.346
0.412
0.238
0.242
0.334
0.328
1.000
旋转前
旋转后
总的特征量
占总变量的百分率%
累计贡献率%
总的特征值
占总变量的百分率%
累计贡献率%
1
3.616
45.199
45.199
2.138
26.731
26.731
2
1.133
14.165
59.365
2.09
26.119
52.85
3
1.075
13.432
72.797
1.596
19.947
72.797
因子
1
2
3
As(μg/g)
0.669
-0.646
-0.01
Cd(ng/g)
0.784
0.171
-0.417
Cr(μg/g)
0.643
0.234
0.493
Cu(μg/g)
0.729
-0.246
0.024
Hg(ng/g)
0.492
0.13
-0.437
Ni(μg/g)
0.686
-0.253
0.523
Pb(μg/g)
0.803
0.112
-0.348
Zn(μg/g)
0.501
0.691
0.267
因子
1
2
3
As(μg/g)
0.003
0.513
-0.31
Cd(ng/g)
0.457
-0.105
-0.031
Cr(μg/g)
-0.176
0.123
0.488
Cu(μg/g)
0.067
0.286
-0.045
Hg(ng/g)
0.41
-0.135
-0.101
Ni(μg/g)
-0.282
0.442
0.219
Pb(μg/g)
0.402
-0.042
-0.025
Zn(μg/g)
0.039
-0.258
0.62
由表2.1.1中所得相关系数观察可得,Ni与As的相关性最好,相关系数为0.605,其次为Cu与As,Ni与Cr相关较好。
且在表3-1-2中可以观察得到累积方差为72.797%,分析得到3个主因子,3个主因子提供了源资料的72.797%的信息,以及主因子1与主因子2的方差贡献率均在26.1%左右,从表3-1-3,表3-1-4综合反映出旋转前后载荷的变量结果基本一致,正交因子解说面得出因子1为Cu与As,因子2为Ni与Cr,这两组重金属元素为相关性最好的两组元素,可作为评定生活区重金属污染的主因子。
工业区主因子的求解:
由附录一表3-2工业区土壤重金属浓度得到模型求解所需的样本空间X,依据生活区主因子的求解方法及如下表格分析求解得,主因子1Cu与Hg,主因子2As与Ni相关性最好,可作为工业区重金属污染的主因子。
As(μg/g)
Cd(ng/g)
Cr(μg/g)
Cu(μg/g)
Hg(ng/g)
Ni(μg/g)
Pb(μg/g)
Zn(μg/g)
As(μg/g)
1.000
Cd(ng/g)
0.329
1.000
Cr(μg/g)
0.380
0.541
1.000
Cu(μg/g)
0.153
0.566
0.920
1.000
Hg(ng/g)
0.181
0.533
0.902
0.983
1.000
Ni(μg/g)
0.690
0.489
0.698
0.503
0.479
1.000
Pb(μg/g)
0.395
0.829
0.675
0.670
0.612
0.578
1.000
Zn(μg/g)
0.518
0.754
0.695
0.622
0.590
0.634
0.739
1.000
旋转前
旋转后
总的特征量
占总变量的百分率%
累计贡献率%
总的特征值
占总变量的百分率%
累计贡献率%
1
5.254
65.67
65.67
3.913
48.907
48.907
2
1.263
15.786
81.455
2.604
32.549
81.455
因子
1
2
As(μg/g)
0.518
0.758
Cd(ng/g)
0.786
0.074
Cr(μg/g)
0.916
-0.206
Cu(μg/g)
0.868
-0.463
Hg(ng/g)
0.845
-0.459
Ni(μg/g)
0.767
0.421
Pb(μg/g)
0.858
0.049
Zn(μg/g)
0.859
0.188
因子
1
2
As(μg/g)
-0.268
0.547
Cd(ng/g)
0.088
0.134
Cr(μg/g)
0.236
-0.032
Cu(μg/g)
0.347
-0.203
Hg(ng/g)
0.342
-0.203
Ni(μg/g)
-0.074
0.356
Pb(μg/g)
0.111
0.126
Zn(μg/g)
0.047
0.216
山区主因子的求解:
由附录一表3-3山区土壤重金属浓度得到模型求解所需的样本空间X,依据山区主因子的求解方法及如下表格分析求解得,主因子1Ni与Cr、主因子2Cd与Pb相关性最好,可作山区土壤重金属污染的主因子。
As(μg/g)
Cd(ng/g)
Cr(μg/g)
Cu(μg/g)
Hg(ng/g)
Ni(μg/g)
Pb(μg/g)
Zn(μg/g)
As(μg/g)
1.000
Cd(ng/g)
-0.291
1.000
Cr(μg/g)
0.113
0.066
1.000
Cu(μg/g)
0.527
0.090
0.364
1.000
Hg(ng/g)
0.075
0.246
-0.006
0.505
1.000
Ni(μg/g)
0.078
0.049
0.945
0.358
-0.045
1.000
Pb(μg/g)
-0.205
0.766
0.107
0.122
0.226
0.028
1.000
Zn(μg/g)
-0.176
0.606
0.627
0.252
0.170
0.629
0.590
1.000
旋转前
旋转后
总的特征量
占总变量的百分率%
累计贡献率%
总的特征值
占总变量的百分率%
累计贡献率%
1
3.042
38.022
38.022
2.457
30.715
30.715
2
2.036
25.445
63.467
2.407
30.088
60.802
3
1.549
19.362
82.829
1.762
22.026
82.829
因子
1
2
3
As(μg/g)
-0.009
0.668
0.478
Cd(ng/g)
0.601
-0.678
0.141
Cr(μg/g)
0.761
0.48
-0.365
Cu(μg/g)
0.517
0.469
0.615
Hg(ng/g)
0.324
-0.074
0.747
Ni(μg/g)
0.737
0.499
-0.409
Pb(μg/g)
0.605
-0.63
0.173
Zn(μg/g)
0.905
-0.173
-0.183
因子
1
2
3
As(μg/g)
0.006
-0.202
0.403
Cd(ng/g)
-0.6
0.393
-0.005
Cr(μg/g)
0.41
-0.065
-0.03
Cu(μg/g)
0.043
-0.005
0.488
Hg(ng/g)
-0.174
0.175
0.43
Ni(μg/g)
0.423
-0.082
-0.053
Pb(μg/g)
-0.058
0.379
0.023
Zn(μg/g)
234
0.23
-0.051
主干道路区主因子的求解
由附录一表3-4主干道路区土壤重金属浓度得到模型求解所需的样本空间X,依据主干道路区主因子的求解方法及如下表格分析求解得,主因子1Cr与Cu、主因子2Pb与Hg相关性最好,可作为主干道路区重金属污染的主因子。
As(μg/g)
Cd(ng/g)
Cr(μg/g)
Cu(μg/g)
Hg(ng/g)
Ni(μg/g)
Pb(μg/g)
Zn(μg/g)
As(μg/g)
1.000
Cd(ng/g)
0.121
1.000
Cr(μg/g)
0.139
0.373
1.000
Cu(μg/g)
0.092
0.424
0.894
1.000
Hg(ng/g)
-0.004
0.211
0.012
0.032
1.000
Ni(μg/g)
0.228
0.351
0.869
0.886
0.040
1.000
Pb(μg/g)
0.060
0.615
0.428
0.506
0.266
0.396
1.000
Zn(μg/g)
0.188
0.294
0.395
0.432
0.118
0.503
0.482
1.000
旋转前
旋转后
总的特征量
占总变量的百分率%
累计贡献率%
总的特征值
占总变量的百分率%
累计贡献率%
1
3.751
46.883
46.883
3.149
39.358
39.358
2
1.287
16.084
62.967
1.889
23.609
62.967
因子
1
2
As(μg/g)
0.105
-0.058
Cd(ng/g)
-0.032
0.391
Cr(μg/g)
0.325
-0.1
Cu(μg/g)
0.306
-0.05
Hg(ng/g)
-0.246
-0.525
Ni(μg/g)
0.329
-0.099
Pb(μg/g)
-0.02
0.415
Zn(μg/g)
0.105
0.163
因子
1
2
As(μg/g)
0.235
-0.131
Cd(ng/g)
0.621
0.458
Cr(μg/g)
0.874
-0.319
Cu(μg/g)
0.906
-0.25
Hg(ng/g)
0.17
0.743
Ni(μg/g)
0.888
-0.32
Pb(μg/g)
0.703
0.477
Zn(μg/g)
0.643
0.115
公园绿地区主因子的求解:
由附录一表1公园绿地区土壤重金属浓度得到模型求解所需的样本空间X,依据公园绿地区主因子的求解方法及如下表格分析求解得,主因子1Cr与Cu、主因子2Pb与Hg相关性最好,可作为公园绿地区重金属污
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 全国大学生 数学 建模 全国 优秀论文 北京理工大学
![提示](https://static.bdocx.com/images/bang_tan.gif)