数学建模.docx

1、数学建模城市表层土壤重金属污染分析模型与研究方案小组成员：10级金融学院8班 10级金融学院8班 10级会计学院12班 城市表层土壤重金属污染分析的模型与研究方案摘要本文针对城市表层土壤重金属污染问题求出重金属元素在该城区的空间分布，建立重金属污染分析评价模型，并作出了相应的研究方案。问题一，利用Matlab作图做出重金属污染物在该城区不同功能区含量的空间分布，及重金属元素在各个功能区的总体分布图，接着采用单因子污染指数法和综合污染指数法来评价城区土壤重金属污染程度分析出不同重金属元素的污染程度，得出以下结论：生活区是中度污染，工业区是中度污染，山区是警戒线，主干道路区是重度污染，公园绿地区是

2、轻度污染；重金属元素主要分布在地势较低的主干道路区、生活区、工业区。问题二，利用主成分分析法将8种主要重金属元素在问题一中受污染严重的工业区、主干道路区、生活区这三个区分别进行提取。分析得主干道路区提取出Pb、 Cd，生活区提取出Cr、 Zn 、Ni，工业区提取出来Ni,同时结合各主要重金属元素的相关系数大小，得出重金属污染的主要原因是由工业、生活、交通这三个方面所引起的。问题三，通过观察每种超标的重金属与工业区、主干道路区、生活区之间的关系的分布图，我们可以判断重金属浓度与所在地理位置的大致关系，然后我们利用Matlab对他们的浓度与地理位置进行拟合，得出重金属浓度与地理位置的具体关系，根据

3、他们之间的数值关系，我们可以判断出污染源的大致位置及污染源的数目问题四，通过增加地质构造、岩性、海拔高度、地形坡度、植被覆盖类型及覆盖率、土壤类型、水资源、地质灾害、人类活动影响及环境质量这10个评价因子，并根据此10个评价因子构建了生态地质环境脆弱性评价指标体系，然后运用熵权模糊数学综合评判的方法，对该城市生态地质环境脆弱性及脆弱性分区进行评价。关键词：重金属元素 因子分析 单因子污染指数 内梅罗综合指数法 数值拟合 熵权模糊数学一、 问题重述随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展

4、城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、5类区，不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（010 厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。附件1列出了采

5、样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。现要求你们通过数学建模来完成以下任务：(1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。(2) 通过数据分析，说明重金属污染的主要原因。(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二、 引言土壤重金属污染是指由于人类活动将重金属带入土壤中，致使土壤中重金属含量明显高于其自然背景含量

6、，并造成生态破坏和环境质量恶化的现象 。由于土壤重金属污染具有多源性、隐蔽性和一定程度上的长距离传输性以及污染后果的严重性，故有人形象地将土壤污染称之为 “化学定时炸弹 ” 。随着工业化、城市化的发展，土壤环境污染问题越来越突出。土壤中重金属不仅对植物的生长造成影响，还通过食物链危害人类健康。受重金属污染的土壤进入到大气和水体中，还会造成大气污染和水污染等环境问题。对土壤中的重金属污染进行分析和评价，目前已有很多报道，因此，土壤重金属污染以及与其相关的评价研究越来越引起人们的关注。同时，近几年，国内外学者针对不同地区已进行了很多土壤环境调查及质量评价方面的研究工作，以土壤重金属元素为例，其主要

7、研究方法为：有针对性地对某城市采取土壤样品，并分析采样点土壤的污染物含量，在此基础上，从多个方面进行土壤重金属污染程度评价，并分析可能的污染源。三、 模型的假设（一）假设所给的采样点的数据的可靠精确的（二）污染源的重金属浓度不再增加四、 符号说明X表示平均值；S表示标准偏差；X+2S表示污染起始值值；表示土壤中污染物i的污染指数；表示金属污染物i的实测值平均值；表示污染物i的评价标准，即是污染起始值；表示污染物i的单因子污染指数；表示最严重污染物的单因子指数；P表示综合污染指数；A表示各个脆弱性等级的区间边界五、 问题一、二的分析与模型的建立5.1 问题一的分析与求解5.1.1 重金属元素在不

8、同区域的空间分布利用Matlab软件，在原有该城区的地形数据上做出以下该城区的地形图：根据原有的附录数据，利用Matlab软件做出砷、镉、铬、铜、汞、镍、铅及锌八种重金属元素在该城区的不同区域的分布图，分别如下图所示：砷（As）元素在不同功能区的含量显示图（1）注：图（1）中，右边的竖直彩条从下到上表示砷（As）的浓度由高到低，以下镉（Cd）、铬（Cr）、铜（Cu）、汞（Hg）、镍（Ni）、铅（Pb）和锌（Zn）在不同功能区的含量显示图类似，见附录一利用Matlab软件做出该城区内生活区、工业区、山区、主干道路区及公园绿地区的分布结构如下图（9）：图（9）注释：区域生活区工业区山区主干道路区公

9、园绿地区符号*.o+p综上可知，八种主要重金属元素在地势较低的主干道路区、工业区及生活区含量都很高，分布较集中。5.1.2 某城市不同功能区重金属含量统计值某城市土壤重金属含量统计表，如下表（1）所示：功能区项目砷(As)镉(Cd)铬(Cr)铜(Cu)汞(Hg)镍(Ni)铅(Pb)锌(Zn)生活区全距2.34-11.4586.80-1044.5018.46-744.469.73-248.8512.00-550.008.89-32.8024.43-472.4843.37-2893.47平均值6.27289.9669.0249.493.0418.3469.11237.01标准差2.15183.68

10、107.8947.16102.95.6672.32443.64工业区全距1.61-21.87114.50-1092.9015.40-285.5812.70-2528.4811.79-13500.004.27-41.7031.24-434.8056.33-1626.02平均值7.24393.1153.41127.54642.3619.8193.04277.93标准差4.24237.5744414.942244.078.3785.37350.83山区全距1.77-10.9940.00-407.6016.20-173.342.29-69.069.64-206.795.51-74.0319.68-11

11、3.8432.86-229.80平均值4.04152.3238.9617.3240.9615.4536.5673.29标准差1.878.3824.5910.7327.8510.4317.7330.94主干道路区全距1.61-30.1350.10-1619.8015.32-920.8412.34-1364.858.57-16000.006.19-142.5022.01-181.4840.92-3760.82平均值5.71360.0158.0562.21446.8217.6263.53242.85标准差3.04239.2475.43111.522010.5611.7932.53384.78公园绿地

12、区全距2.77-11.6897.20-1024.9016.31-96.289.04-143.3110.00-1339.297.60-29.1026.89-227.4037.14-1389.39平均值6.26280.5443.6430.19114.9915.2960.71154.24标准差2.02235.8414.8422.68224.284.9745.84230.92主要重金属元素的背景值，如下表（2）所示：元素平均值（X）标准偏差（S）污染起始值（X+2S）As (g/g)3.60.95.4Cd (ng/g)13030190Cr (g/g)31949Cu (g/g)13.23.620.4Hg

13、 (ng/g)35851Ni (g/g)12.33.819.9Pb (g/g)31643Zn (g/g)691497不同功能区土壤中重金属元素的密度平均值，如下表（3）所示：功能区项目砷(As)镉(Cd)铬(Cr)铜(Cu)汞(Hg)镍(Ni)铅(Pb)锌(Zn)生活区平均值6.27289.9669.0249.493.0418.3469.11237.01工业区平均值7.24393.1153.41127.54642.3619.8193.04277.93山区平均值4.04152.3238.9617.3240.9615.4536.5673.29主干道路区平均值5.71360.0158.0562.21

14、446.8217.6263.53242.85公园绿地区平均值6.26280.5443.6430.19114.9915.2960.71154.24由表（3）得出以下直方图（）：图（）从图（）可知：砷（As）含量顺序是工业区生活区主干道路区公园绿地区山区；镉（Cd）含量顺序是工业区生活区主干道路区公园绿地区山区；铬（Cr）含量顺序是生活区主干道路区工业区公园绿地区山区；铜（Cu）含量顺序是工业区主干道路区生活区公园绿地区山区；汞（Hg）含量顺序是工业区主干道路区公园绿地区生活区山区；镍（Ni）含量顺序是工业区生活区主干道路区山区公园绿地区；铅（Pb）含量顺序是工业区生活区主干道路区公园绿地区山区；

15、锌（Zn）含量顺序是工业区主干道路区生活区公园绿地区山区5.1.3 该城区内不同区域重金属的污染程度（1）土壤重金属污染程度评价方法本文采用单因子污染指数法和综合污染指数法来评价该城区土壤重金属污染程度，污染程度等级用污染指数来表示。单因子污染指数法公式如下: 综合污染指数采用内梅罗综合指数法：P=根据以上单因子污染指数和综合污染指数的计算方法得出以下表格，分析出各个功能区的污染程度及重金属含量的高低。土壤污染评价等级如下表： 污染指数P0.70.7P11P22P3P3污染等级安全警戒线轻度污染中度污染重度污染经过计算可的五个区的八种重金属元素的单因子污染指数如下表：元素生活区工业区山区主干道

16、路公园绿地区As (g/g)1.161.340.751.061.16Cd (ng/g)1.532.070.801.891.48Cr (g/g)1.411.090.801.180.89Cu (g/g)2.426.250.853.051.48Hg (ng/g)1.8212.600.808.762.25Ni (g/g)0.921.000.780.890.77Pb (g/g)1.612.160.851.481.41Zn (g/g)2.442.870.762.501.59五个区的重金属综合污染指数（P）如下表：功能区生活区工业区山区主干道路区公园绿地区P2.099.280.826.461.87（2）综合

17、评价将各区单因子污染指数与土壤污染分级标准对比可得到如下结果：生活区：Ni属于警戒级As，Cd，Cr，Hg，Pb都属于轻度污染，Cu，Zn属于中度污染。工业区：Ni属于警戒级；As，Cr属于轻度污染；Cd，Pb，Zn属于中度污染，Cu，Hg属于严重超标，属于重度污染。山区：As，Cd，Cr，Cu，Hg，Ni，Pb，Zn均处于警戒级。主干道路区：Ni属于警戒级，As，Cd，Cr ，Pb，属于轻度污染，Zn属于中度污染，Cu，Hg严重超标，属于重度污染。公园绿地区：Cr，Ni属于警戒级，As，Cd，Cu，Pb，Zn，属于轻度污染，Hg属于中度污染。从综合分析指数分析可得到不同区的污染程度如下表：功

18、能区生活区工业区山区主干道路区公园绿地区污染程度中度污染重度污染警戒线重度污染轻度污染5.2 问题二的分析与求解土壤中重金属元素主要有自然来源和人为干扰输入两种途径。在自然因素中，成土母质和成土过程对土壤重金属含量的影响很大，在各种人为因素中，则主要包括工业、生活和交通等来源引起的土壤重金属污染。5 .2.1 利用因子分析各种金属元素表A个金属元素的相关矩阵相关矩阵AsCdCrCuHgNiPbZn相关As10.2550.1890.160.0640.3170.290.247Cd0.25510.3520.3970.2650.3290.660.431Cr0.1890.35210.5320.1030.

19、7160.3830.424Cu0.160.3970.53210.4170.4950.520.387Hg0.0640.2650.1030.41710.1030.2980.196Ni0.3170.3290.7160.4950.10310.3070.436Pb0.290.660.3830.520.2980.30710.494Zn0.2470.4310.4240.3870.1960.4360.4941Sig.（单侧）As000.0020.126000Cd0000000Cr0000.033000Cu0.002000000Hg0.12600.03300.03300Ni00000.03300Pb000000

20、0Zn0000000由相关矩阵可看出Cr和Ni的相关性最好，是0.716，其次是Cd和Pb，相关系数是0.660，Cu和Cr，Pb和Cu相关系数较好分别为0.532，0.520;Ni和Cu，Zn和Pb的的相关系数分别为0.495；0.494.，其他的相关性不太好，从成因上分析，相关性好的元素可能在成因和来源上有一定的关联。解释的总方差成份初始特征值提取平方和载入旋转平方和载入合计方差的 %累积 %合计方差的 %累积 %合计方差的 %累积 %13.5644.544.53.5644.544.51.82722.83422.83421.1514.37758.8771.1514.37758.8771.5

21、2519.06341.89730.96512.06370.9410.96512.06370.9411.04513.06854.96540.7689.59680.5370.7689.59680.5371.03212.89867.86350.5787.2287.7560.5787.2287.7561.0212.74480.60760.4325.39993.1560.4325.39993.1561.00412.54993.15670.3013.76996.92480.2463.076100旋转成份矩阵a成份12345As.131.154.023.970.076Cd.170.877.082.102.1

22、32Cr.882.209.002.014.146Cu.614.362.505-.022.029Hg.013.134.953.030.083Ni.864.089.019.222.196Pb.195.832.191.121.208Zn.271.285.101.091.904在累积方差90%的前提下分析得到五个主因子.从表B 可知旋转后主因子1和主因子2的方差贡献率均在20%左右，主因子3到主因子方差贡献率在12.54913.068之间，这可解释为因子1和因子2可能为该城市土壤重金属污染的最重要污染源。其他四个因子对该城市重金属污染有重要影响，由旋转成分矩阵可知因子1为Cr、Ni、Cu，因子2为Cd

23、和Pb，因子3为Hg，因子4为As，因子5为Zn。Cr、Ni和Cu，Cd和Pb，可能是同一来源，并且这两组元素的相关性最好。5.2.2 土壤重金属污染的主要原因分析通过刚才对主要重金属元素的数据分析可知某些重金属类污染元素异常与主干道路区、工业和生活区等在空间分布上往往具有很好的对应性，并且它们地处的位置比较低容易使重金属离子在降雨、降雪和自身所载水分的作用下被淋融下来，进入土壤后造成土壤的重金属污染，因此可以初步判定人为干扰是这类异常的主要成因。下面再借助第一问中各功能区的各重金属元素含量顺序的前提下就对造成土壤重金属污染的主要原因进行分析：（1）由于Cr的主要来源是电镀、燃料、制药、皮革、

24、颜料等铬化合物制造企业中排放的废物，有报道南京某合金厂周围土壤中的Cr大大超过土壤背景值，Cr污染以工厂烟囱为中心，范围达到1.5 km2。Ni的主要来源是工业污染和矿山开采，因此可以推测重金属Cr的主要来源是工业污染。（2）由于化工生产、金属制造等多种工业原料中都含有Pb、Cd，释放后进入土壤，煤与原油中也含有Cd,在燃烧过程中释放到大气中，而汽车尾气是城市中Pb的主要来源之一，FAKAYODE和OWOLABI研究了尼日利亚不同交通密度公路边表层土壤中Pb、Cd的分布，结果表明，重金属含量在车流密度大的公路两侧土壤中要高于车流密度小的公路两侧土壤，且随着距公路距离的增大，重金属含量快速降低，

25、到距公路50 m左右的地方，重金属含量基本降低到背景值水平。因此可以推测重金属Pb、Cd的主要来源是工业污染、交通污染和生活污染。（3）由于化工产业等排放的未经处理或处理不彻底的废水中含有Hg元素，而各种化妆品厂排放的工业废水中也含有大量的Hg，因此可以推测重金属Hg的主要来源是工业污染。（4）由于土壤中As的来源可以分为自然源和人为源，前者主要是地质岩石中的As；后者主要来自农药、化肥和矿山冶炼，MESHALIKNA等研究了俄罗斯一硫酸生产厂周围土壤中元素的污染及其空间变异后发现，在距烟囱12 km外的土壤中仍能监测到高含量的S、V和As。因此可以推测重金属As的主要来源是工业污染，（5）由

26、于汽车尾气与生活垃圾中中含有大量的铜，等研究发现在公路两侧50m的距离有被污染的痕迹，每月每平方米累积的易溶性污染物在440 g。进入环境的强度顺序为：Cu、Pb、Co、Fe和Zn。对武汉市垃圾堆放场，杭州铬渣堆放区附近土壤中重金属含量的研究发现，这些区域土壤中Cd、Hg、Cr、 Cu、Zn、Pb、As等重金属含量均高于当地土壤背景值。因此可以推断重金属Cu的主要污染来源是交通污染与生活污染。（6）由于Zn是汽油添加剂的成分之一，汽车部件的电镀材料，润滑剂以及轮胎等均不同程度都含有Zn，而机动车辆排放的废气以及车辆构造磨损物都会释放出Zn, 另外工业制造也会产生大量的Zn,例如沈阳冶炼厂冶炼

27、Zn 的过程中产生的矿渣主要含Zn和Cd，1971年开始堆放在一个洼地场所，其浸入液中Zn、Cd含量分别达6.6103 mgL-1和7.5103 mgL-1，目前已扩散到离堆放场700 m以外的范围，重金属污染物质量浓度是以同心圆状分布因此可以推测重金属Zn地主要来源是交通污染和工业污染。综合上述结论可得：该城区重金属污染的主要原因是工业区中废物的排放、原料的使用、生产的产品等多种方面都含有重金属元素，汽车尾气的排放、车辆构造的磨损物等也会释放重金属元素；生活区中密集的人口对煤炭的大量燃烧、废水的不当处理也会产生重金属元素。即工业，交通，密集的人口活动是重金属污染的主要原因。六、 重金属污染物污染源分析模型6.1 问题分析与模型的建立6.1.1问题分析若要了解污染源的传播特征并确定污染源的位置，必须要知道该区域受什么污染以及污染情况由前两问的