城市表层土壤重金属污染分析.docx

1、城市表层土壤重金属污染分析2011高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的参

2、赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 江西财经大学 参赛队员 （打印并签名）：1. 林礼连 2. 尹宁宁 3. 林顺 指导教师或指导教师组负责人（打印并签名）： 日期： 2011 年 9 月 12 日赛区评阅编号（由赛区组委会评阅前进行编号）：2011高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：城市表层土壤重金属污染分析摘 要目前，人类活动对环境的影响作用越来越大，使得环境质量

3、日益下降。为了更有效的查证城市地质环境的异常现象，本文通过对某城市五个功能区八种重金属的含量进行研究，来确定此城市受污染现状，并建立扩展模型以更全面科学的分析地质环境的演变模式。对于问题一，先对数据进行统计分析，获得重金属在5类区中的超标情况。用surfer8.0软件画等高线的方法分别画出该城市城区土壤8种主要重金属元素的含量浓度分布图，然后用现在普遍使用的综合指数法和内梅罗指数法分析不同区域重金属的污染程度。最后综合得到污染程度：工业区交通区生活区公园绿化区山区。对于问题二，根据污染指数和污染指标两项指标来定量分析各个区域受污染程度，重金属污染主要来源于工业区的工业废水废气废渣、生活区的废弃

4、物、交通要道的汽车尾气，同时由于Hg和Cd两种元素的毒性明显大于其他六种金属而对环境造成了主要危害。对于问题三，由定性的重金属污染物的传播特征，通过matlab软件画出8种重金属密度与地理区域的关系图可得结论，可得结论：污染源就在浓度最高点。用excel进行数据处理，并删除一系列不可能点，获得最后的污染源。重金属污染物的传播特征的定量分析即用三次样条拟合获得金属污染浓度与距离关系。对于问题四，我们综合考虑了机械弥散、对流、扩散和吸附作用对重金属在土壤中运移情况的影响，得出在Z方向上重金属流通量与土壤含水率、弥散度、孔隙的平均流速、土壤溶质浓度、溶质分子在土壤中的扩散系数、土壤最大吸附通量及土壤

5、吸附系数之间的关系（相应符号说明详见后文）：关键字：污染指数 污染指标 综合指数法 内梅罗指数法 传播特征一、问题重述随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、5类区，不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方

6、公里1个采样点对表层土（010 厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。现要求你们通过数学建模来完成以下任务：(1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。(2) 通过数据分析，说明重金属污染的主要原因。(3) 分析重金属污染物的

7、传播特征，由此建立模型，确定污染源的位置。(4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二、数据分析及处理针对问题1，先对数据进行统计分析，利用统计指标衡量城区土壤8种主要重金属元素的浓度特征，获得重金属在5类区中的超标情况。用surfer8.0软件画等高线的方法分别画出该城市城区土壤8种主要重金属元素的含量浓度分布图，并由统计分析表和浓度分布图对重金属污染进行空间分析。污染程度与重金属含量正相关，所以需比较在各区域的含量与标准值之间关系，并把它们的比值污染指数作为标准，污染指数越大，受污染程度越大，反之越小。针对问

8、题2，分析重金属污染的主要原因，首先要考虑的是重金属的含量高于标准值的比率，即污染指数。污染指数越高，说明该地区受这种金属污染越大，反之越小。而不同重金属毒性差别较大，故需同时考虑污染指数和毒性两种因素，我们选择两者之积作为对环境污染的指标，数值越大，危害越大。针对问题3，先对重金属污染物的传播特征进行定性分析，即在不考虑水流、风向、地质等因素的影响时，重金属在土壤中的传播只和距离有关系。分析得到结论：污染源就在浓度最高点。再删除一系列不可能的点，获得最后的污染源。重金属污染物的传播特征的定量分析即用三次样条拟合获得金属污染浓度与距离关系。针对问题4，鉴于问题3中忽略了许多影响重金属在土壤中运

9、移的因素，如土壤紧实程度（土壤密度）、土壤含水率、地质等。而这些因素通过影响机械弥散、对流、扩散和吸附作用来决定重金属的运移，所以我们以流通量为指标，先分别求出机械弥散、对流、扩散和吸附作用对重金属运移的影响模式，再综合考虑得出四种作用共同影响下的重金属运移模式。 三、模型假设和符号设定 3.1模型假设(1)同一种重金属毒性可以按含量叠加；(2)不考虑土壤中其他成分对重金属的毒性的作用；(3)不考虑不同种重金属的相互作用对毒性的影响；(4)附表中所给数据都是有效的，即不考虑偶然因素对浓度的影响；(5)假设山区和公园绿地区不存在污染源。3.2符号设定(元素As、Cd、Cr、Cu、Hg、Ni、Pb

10、、Zn单因子污染指数平均浓度标准值毒性系数污染程度值整体污染程度值变异系数背景值超出背景值比例单因子平均污染指数最大单因子污染指数内梅罗污染指数取样点到污染源的距离取样点位置污染源位置溶质的机械弥散通量土壤含水率溶质的机械弥散系数C土壤溶质浓度Z方向弥散度孔隙的平均流速土壤溶质对流通量土壤水分通量溶质分子在土壤中的扩散通量溶质分子在土壤中的扩散系数土壤吸附作用流通量最大吸附通量K土壤吸附系数四、模型建立与求解4.1问题1的模型建立与求解4.1.1重金属元素空间分布某城市城区土壤8种主要重金属元素的浓度特征，可利用平均浓度、标准差、变异系数、背景值、超出背景值比例五种统计指标来衡量，其中变异系数

11、计算公式： (公式1)超出背景值比例计算公式： (公式2)统计结果见表1。表1 不同功能区城市土壤重金属浓度统计功能区元素As Cd Cr Cu Hg Ni Pb Zn 生活区平均浓度(ng/g)6.27 289.969.02 49.40 93.04 18.34 69.11 237.0标准差0.93093.683.8614变异系数0.140.100.130.070.090.200.090.06背景值3.61303113.23512.33169超出背景值的比例0.741.231.222.741.650.491.222.43工业区平均浓度(ng/g)7.25 393.153.41 127.5642

12、.319.81 93.04 277.9标准差0.93093.683.8614变异系数0.120.080.170.030.010.190.060.05背景值3.61303113.23512.33169超出背景值的比例1.012.020.728.6617.350.612.003.03山区平均浓度(ng/g)4.04 152.338.96 17.32 40.96 15.45 36.56 73.29 标准差0.93093.683.8614变异系数0.220.200.230.210.200.250.160.19背景值3.61303113.23512.33169超出背景值的比例0.120.170.260.

13、310.170.260.180.06交通区平均浓度(ng/g)5.71 360.058.05 62.21 446.817.62 63.53 242.9标准差0.93093.683.8614变异系数0.160.080.160.060.020.210.090.06背景值3.61303113.23512.33169超出背景值的比例0.591.770.873.7111.770.431.052.52公园绿化区平均浓度(ng/g)6.26 280.543.64 30.19 115.015.29 60.71 154.2标准差0.93093.683.8614变异系数0.140.110.210.120.070.

14、250.100.09背景值3.61303113.23512.33169超出背景值的比例0.741.160.411.292.290.240.961.24由上表可知：八种重金属都具有低变异性，变异系数反映数据的离散性，其中As、Cr、Ni的变异系数相对较大，说明土壤中As、Cr、Ni浓度在一定程度上受到人为影响，在空间上有相对的变化。以背景值为标准，生活区除Ni之外，超标率49%，其它元素严重超标，Zn超标最高；工业区Cu、Hg的超标很大，Zn次之，超标最小的是As、Cr、Ni三种重金属；山区超标率极低，最高超标率不超过35%，为Cu，最低为Zn元素；交通区Hg超标严重，Cu、Zn次之；超标率最低

15、为Ni，但也超标43%；公园绿化区Hg超标最大，最低超标率为16%、24%，分别为Cd、Ni。综上所述，生活区、工业区、交通区重金属超标严重，而山区环境最好，公园绿化区次之。用surfer8.0软件画等高线的方法分别画出该城市城区土壤8种主要重金属元素的含量浓度分布图，其中平面（）代表位置，高度代表相应某重金属的密度，越高密度越大。图一 As的浓度分布图 图二 Cd的浓度分布图图三 Cr的浓度分布图 图四 Cu的浓度分布图图五 Hg的浓度分布图 图六 Ni的浓度分布图图七 Pb的浓度分布图 图八 Zn的浓度分布图从空间分布图可以看出：Cd浓度分布分散，As次之，其余6种金属浓度分布有明显的峰值

16、，可观测出最大浓度点的所在位置。由上图结合表1中超标率分析得出：Cr浓度生活区超标最高；Cu浓度工业区超标最高；Hg浓度工业区超标最高，交通区也很高；Ni浓度工业区相对其他区更高；Pb浓度工业区超标最高，生活区次之；Zn浓度工业区超标最高，交通区、生活区次之。4.1.2分析不同区域重金属的污染程度目前分析土壤中重金属的污染程度的评价方法有很多，我们采用较为普遍的两种：综合指数法和内梅罗指数法。（一） 综合指数法：对于同一种重金属而言，超过标准值越高，即污染指数越高，危害就越大。对于不同的金属，由于毒性大小相差较大，即使污染指数相同，危害数也存在较大差异。综合指数法主要是结合单因子污染指数和各金

17、属的毒性系数而得到单个金属的污染程度值，并由每个区域单个金属的综合污染指标累加得到该区域的综合污染指标。其中，单因子污染指数计算方法为： （公式3）为平均值，即实测值的平均值，为标准值，即每种元素的背景值加2倍的标准差。 表2 各区域不同元素的污染指数功能区元素As Cd Cr Cu Hg Ni Pb Zn 生活区平均值(ng/g)6.27 289.96 69.02 49.40 93.04 18.34 69.11 237.01 标准值(ng/g)5.41904920.45119.94397污染指数1.161.531.412.421.820.921.612.44工业区平均值(ng/g)7.25

18、393.11 53.41 127.54 642.36 19.81 93.04 277.93 标准值(ng/g)5.41904920.45119.94397污染指数1.342.071.096.2512.601.002.162.87山区平均值(ng/g)4.04 152.32 38.96 17.32 40.96 15.45 36.56 73.29 标准值(ng/g)5.41904920.45119.94397污染指数0.750.800.800.850.800.780.850.76交通区平均值(ng/g)5.71 360.01 58.05 62.21 446.82 17.62 63.53 242.8

19、5 标准值(ng/g)5.41904920.45119.94397污染指数1.061.891.183.058.760.891.482.50公园绿化区平均值(ng/g)6.26 280.543.64 30.19 114.915.29 60.71 154.2标准值(ng/g)5.41904920.45119.94397污染指数1.161.480.891.482.250.771.411.59为使评价结果更好的反映客观为农业生产服务评价方法应该考虑各种金属的毒性系数，八种重金属的毒性系数如下表3。表3 各种重金属的毒性系数元素AsCdCrCuHgNiPbZn毒性系数103025400510每个区域的整

20、体污染程度值计算公式为 （公式4）可得每个区域的综合污染指标如下表：表4 各区域综合污染指标区域生活区工业区山区交通区公园绿地区综合污染指标177.76652.2281.30467.91178.23根据公式4，各区域的综合污染指标的标准值的计算公式为：所以，分析表三数据并和比较可知，工业区污染最为严重，其次是交通区，生活区和公园绿地区属于轻度污染，山区完全没有污染。（二）内梅罗指数法本题是多种污染物同时污染某一区域土壤的现象，单因子评价难以表示它们的整体污染水平，这时就采用多因子评价方法内梅罗指数法计算土壤各元素在每种类型土地中的污染指数， 并进行污染的分级， 此方法在可反映污染的程度的同时还

21、可以突出污染指数最大的污染物对环境质量的影响和作用。内梅罗指数法，计算公式为： (公式5)其中，为单因子平均污染指数，为最大单因子污染指数。表5 内梅罗污染指数评价标准等级内梅罗污染指数污染等级1清洁（安全）2尚清洁（警戒线）3轻度污染4中度污染5重污染由内梅罗指数法，得到内梅罗综合污染指数表。表6 内梅罗综合污染指数表分级指数As Cd Cr Cu Hg Ni Pb Zn 生活区1.912.042.002.432.161.852.072.44工业区8.969.0268.949.9412.608.939.0379.13山区0.800.830.820.850.830.810.850.80交通区6

22、.246.346.256.568.766.236.286.44公园绿化区1.791.911.711.912.251.681.881.95可得工业区、交通区为重污染地区，生活区中度污染、公园绿化区轻度污染、山区尚清洁。As、Cd、Cr、Cu、Ni、Pb、Zn的综合污染指数大小均依次为工业区交通区生活区公园绿化区山区，Hg的综合污染指数大小依次为工业区交通区公园绿化区生活区山区；由综合指数法和内梅罗指数法综合得到不同区域重金属的污染程度：工业区交通区生活区公园绿化区山区。4.2问题2的模型建立与求解由问题1得出结论进行分析。生活区，Zn和Cu含量超标比较严重，其次是Hg、Cd，这些重金属主要来源于

23、生活垃圾，污染过程主要包括污水处理中产生污泥的堆放、垃圾渗滤液的泄漏、这些污水处理厂产生的污泥中含有大量的重金属，会对土壤造成二次污染，城市垃圾在焚烧过程中产生的飞灰及堆放填埋过程中产生的渗滤液中的重金属也会严重超标。工业区，工业能源大都以煤、石油类为主，它们是环境中Hg、Cu、Zn、Cd等重金属污染的主要来源，并且Hg的含量超标现象尤为严重，高达12.60，其次是Cu、Zn、Cd。在采矿、选矿、冶炼、锻造、加工、运输等工业生产过程中会产生大量的重金属污染。排放的废水、废渣等直接进入水体及土壤中，废气中的重金属经沉降也进入土壤等环境中，从而使得环境中重金属浓度严重超标。山区由于几乎没有人为影响

24、，又不存在垃圾的排放与堆积，在自然环境的自我净化作用下，重金属含量缓慢减少，最终导致山区的重金属少于背景值，环境良好。交通区，Hg含量明显偏高，其次是Cu和Zn，污染源主要是汽车轮胎及尾气。公园绿地区，相对于其他受污染地区，污染程度相对较低，只有Hg含量高于标准值2.25，其他重金属As、Cd、Cu、Pb、Zn略高于标准值，这说明公园绿地人的活动较少，对环境影响不大，环境质量较好。总观五个区域的单因子污染指标和综合污染指数，由于Hg和Cd的毒性系数偏高于其他重金属，故Hg和Cd的污染指标都占较大比重，即对环境形成的危害较大。4.3问题3的模型建立与求解在不考虑水流、风向、地质等因素的影响时，重

25、金属在土壤中的传播只和距离有关系，浓度会随海拔增高而降低、随水平距离增远而降低，即扩散的量随着距离的增加而迅速减少，可以得到浓度最高点为污染源所在地，其他点的重金属是由于污染源的扩散而形成的。结合上面的空间分布图也可以看出距最高点较远的地方浓度加速减少，所以污染源就在浓度最高点。我们在选取过程中，先根据空间分布图选取峰值点，再由多种不合理情况舍弃部分取样点。1、由于山区和公园绿化区污染程度较低，故去掉山区和公园绿化区污染较高点。2、浓度较高点可能集中在某一个小范围之内，处理方法可为选取范围内浓度最高点、舍弃其他的点。用excel进行一系列数据处理，获得污染源的位置。表7 污染源的位置元素x(m

26、)Y(m)海拔(m)功能区As (g/g)As181341004641430.1312696302427423.72474272939221.87Cd21439113834541619.832996018441213.5270822952241123.963951044342967.7Cr3573621351744.461156310298124306.02Cu23833692722528.48Hg27082295224160001369423573341490115248910616413800Ni3299601844142.5Pb4777489781472.483299601844181.48Zn1379796211843760.82932843112412893.4732996018441818.471269630242741635.92已定性分析重金属污染物的传播特征，可使用污染源的位置定量分析重金属污染物的传播特征。在据污染源的范围内选取取样点，距离计算公式： (公式6)其中为取样点位置，为污染源位置。通过三次样条拟合来获得污染浓度与距离关系，三次样条拟合模型： (公式7)应用matlab软件进行编程，得到拟合结果：表8 各种重金属污染浓度与距离的关系AsCdCr