数学建模全国奖论文.docx

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

2、： 所属学校（请填写完整的全名）： 参赛队员（打印并签名）：1. 2. 3. 指导教师或指导教师组负责人 （打印并签名）： 日期：2011 年9月12 日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分OOOOOO备注OOOO全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：基于系统综合评价的城市表层土壤重金属污染分析摘要本文针对城市表层土壤重金属污染问题，首先对各重金属元素进行分析，然后对各 种重金属元素的基本数据进行统计分析及无量纲化处理， 再对各金属元素进

3、行相关性分析，最后针对各个问题建立模型并求解。针对问题一，我们首先利用EXCEL和SPSS统计软件对各金属元素的数据进行处理， 再利用Matlab软件绘制出该城区内8种重金属元素的空间分布图最后通过内梅罗污染/ 2 2 卍 2模型：卩综=P平均+Pmax )，其中P平均为所有单项污染指数的平均值， Pmax为土壤环境中2各单项污染指数中的最大值。得到各区内梅罗综合污染指数，进而求得污染程度为:功能 区生活区工业区山区主干道路区公园绿地区P综2.09049.27690.82426.46251.8688污染土壤、作物受土壤、作物受尚清洁土壤、作物受超标程度中度污染严重污染严重污染针对问题二,我们首

4、先利用EXCELS件画出8种元素在各个区内相对含量的柱状图，由图可以明显地看出各个区内各种元素的污染情况，然后再根据重金属元素污染来 源及传播特征进行分析，可以得出工业区及生活区重金属的堆积和迁移是造成污染的主 要原因，Cu、Hg、Zn主要在工业区和交通区如公路、铁路等交通设施的两侧富集，随 时间的推移，工业区、交通区的土壤重金属具有很强的叠加性， 受人类活动的影响较大 同时城市人口密度，土地利用率，机动车密度也是造成重金属污染的原因。针对问题三，我们从两个方面考虑建模即以点为传染源和以线为传染源。针对以点为传染源我们建立了两个模型：无约束优化模型 D = m lx i - x f ly i

5、- y f，得到污染源的位置坐标5567 ,6782 ；有衰减的扩散过程模型得位置坐标(8500,5500 )，模型为:2 2,2 ：一 u 2 ：一 u , 2b 2 - c 2 k u：x:：y :z针对以线为传染源我们建立了丫 =u be 模型，并通过线性拟合分析线性污染源的位置。针对问题四，我们在已有信息的基础上，还应收集不同时间内的样点对应的浓度以 及各污染源重金属的产生率。根据高斯浓度模型建立高斯修正模型，得到浓度关于时间 和空间的表达式C二C。上。在本题求解过程中，我们所建立的模型与实际紧密联系， 有很好的通用性和推广性 但在求点污染源时，我们假设只有一个污染源，而实际上可能有多

6、个点污染源，从而使 得误差增大，或者使污染源的位置够不准确。关键词 内梅罗污染模型 无量纲化 相关性 回归模型 高斯浓度模型一 问题重述随着城市经济的快速发展和城市人口的不断增加， 人类活动对城市环境质量的影响 日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开 展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关 注的焦点。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区 等，分别记为1类区、2类区、5类区，不同的区域环境受人类活动影响的程度 不同。现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距 1

7、 公里 左右的网格子区域，按照每平方公里 1 个采样点对表层土（ 010 厘米深度）进行取样、 编号，并用 GPS 记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的 多种化学元素的浓度数据。另一方面，按照 2公里的间距在那些远离人群及工业活动的 自然区取样，将其作为该城区表层土壤中元素的背景值。现要求你们通过数学建模来完成以下任务：1） 给出 8 种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属 的污染程度。2） 通过数据分析，说明重金属污染的主要原因。3） 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。4） 分析你所建立模型的优缺点，为更好地研究城市

8、地质环境的演变模式，还应收集 什么信息？有了这些信息，如何建立模型解决问题？二 基本假设1） 假设土壤中重金属元素以离子态存在，不存在别的形态；2） 假设扩散过程已经达到恒稳状态，离子浓度不在变化；3） 取样点的数据较好的反应了该地区的污染物浓度；4） 假设在确定污染源的位置时， 只考虑元素在二维平面上的扩散， 不考虑在垂直方 向（海拔）的扩散；5） 假设该区域土壤都为粘土；6） 假设取样获得数据是处于标准状态下；三符号说明符号符号的意义土壤中污染物i的污染指数土壤中污染物才的背景值利用单因子质量评价模型所得到的百分比构造评价矩阵任意两种元素i,j的相关系数因素重要程度系数内梅罗综合指数Dt重

9、金属离子在粘土介质中得扩散系数重金属元素的来源方向污染源距被选取点的距离Qyt 污染源的位置Jc表示的是土壤溶质对流通量表示土壤水分通量土壤的平均孔隙流速土壤的体积含水率Jd溶质在自由水中的分子扩散通量Do表示的是溶质分子在自由水中的扩散系数Ds表示的是溶质分子在土壤中的扩散系数J ds表示的是溶质分子在土壤中的扩散通量Jh表示的是溶质的机械弥散通量Dh溶质的机械弥散系数Z弥散度J dh表示由于水动力弥散引起的溶质运移通量D表示的是有效水动力弥冃攵系数四、数据处理4.1对重元素的分析城市工业“三废”排放，金属采矿和冶炼，家庭燃煤，生活垃圾，汽车尾气排放都 增加了城市土壤重金属的负荷。重金属污染

10、环境的主要有汞、铅、铬、锌镉、铜等。其 中汞的毒性最大，铬、铅、锌等也有相当大毒性。此外还有砷，砷虽不属于金属 但它的毒性与重金属相似，因此归于重金属一类阐述，称为类金属。目前对我国土壤污染比 较普遍的重金属有汞、铬、砷。根据该城区重金属污染的情况，下面对重金属在土壤污 染中的来源及传播途径作简要介绍。4.1.1砷元素该元素毒性很低，水体中含砷污染物主要来自砷和含砷金属矿的开采、冶炼，以及 和砷化物为原料的玻璃、颜料、药物、纸张的生产都可产生含砷的废水，造成水体的砷 污染。砷及砷化物在水中会在水生物体内累积，但累积程度比其他重金属要低。砷和砷 化物，一般可通过水、大气和食物进入人体。4.1.2

11、镉元素当环境受到镉污染后，镉可在生物体内富集，通过食物链进入人体引起慢性中毒。 镉的主要污染源是电镀、采矿、冶炼、染料、电池和化学工业等排放的废水。相当数量 的镉通过废气、废水、废渣排入环境，造成污染。镉对土壤的污染主要有气型和水型两 种。气型污染主要来自工业废气。镉随废气扩散到工厂周围并自然沉降，蓄积于工厂周 围的土壤中，可使土壤中的镉浓度达到 40ppm。水型污染主要是铅锌矿的选矿废水和有 关工业（电镀、碱性电池等）废水排入地面水或渗入地下水引起。4.1.3铬元素对水体污染的铬主要来源于电镀、制革、铝盐生产以及铬矿石开采所排放的废水。 是我国水体中一种普遍的污染物。水体中铬污染主要是三价铬

12、和六价铬，它们在水体中 的迁移转化有一定的规律性。4.1.4铜元素铜（Cu）及其化合物在环境中所造成的污染称为铜污染。主要污染来源是铜锌矿的开 采和冶炼、金属加工、机械制造、钢铁生产等。冶炼排放的烟尘是大气铜污染的主要来 源。世界铜的年迁移量为：岩石风化20万吨,河流输送11万吨4.1.5汞元素汞是在常温下唯一呈液态的金属元素。人类活动造成水体汞污染，主要来自氯碱、 塑料、电池、电子等工业排放的废水。由于天然本底情况下汞在大气、土壤和水体中均 有分布，所以汞的迁移转化也在陆、水、空之间发生。4.1.6镍元素镍污染是由镍及其化合物所引起的环境污染。大部分煤含有微量镍 ,通过燃烧过程被释放出来,这

13、是大气中镍的主要来源。镍可以在土壤中富集。土壤中的镍主要 来源于岩石风化，大气降尘，灌溉用水（包括含镍废水） ，农田施肥，植物和动物残体的腐烂等。全世界每年镍的迁移状况是：岩石风化量为 320 000吨，河流输送量为19 000吨，开采量为560 000吨，矿物燃料燃烧排放5 600吨。4.1.7铅元素铅对环境的污染，一是由冶炼、制造和使用铅制品的工矿企业，尤其是来自 有色金属冶炼过程中所排出的含铅废水、废气和废渣造成的。二是由汽车排出的 含铅废气造成的，汽油中用四乙基铅作为抗爆剂（每公斤汽油用 13克），在汽油燃烧过程中，铅便随汽车排出的废气进入大气，成为大气的主要铅污染源4.1.8锌元素锌

14、在土壤中富集，会使植物体中也富集而导致食用这种植物的人和动物受害。 金属锌本身无毒，但在焙烧硫化锌矿石、熔锌、冶炼其他含有锌杂质的金属的过 程中，以及在铸铜过程中产生的大量氧化锌等金属烟尘，对人有直接的危害。其他 如橡胶轮胎的磨损以及煤的燃烧也是大气锌污染的原因。各种工业废水的排放是引 起水体锌污染的主要原因。4.2对基本数据的分析用EXCELL软件和SPSS统计软件处理数据如表1所示：表1功能区丿元糸As（卩 g/g）Cd（n g/g）Cr（卩 g/g）Cu（卩 g/g）Hg（ng/g）Ni（卩 g/g）Pb（卩 g/g）Zn（卩 g/g）生 活 区平均值6.27289.9669.0249.

15、4093.0418.3469.11237.01取大值11.451044.50744.46248.85550.0032.80472.482893.47最小值2.3486.8018.469.7312.008.8924.4343.37标准差2.15183.68107.8947.16102.905.6672.33443.64变异系数0.340.631.560.951.110.311.051.87工业区平均值7.25393.1153.41127.54642.3619.8193.04277.93最大值:21.871092.90285.582528.4813500.0041.70434.801626.02最

16、小值1.61114.5015.4012.7011.794.2731.2456.33标准差1.61114.5015.4012.7011.794.2731.2456.33变异系数:0.220.290.290.100.020.220.340.20山 区平均值4.04152.3238.9617.3240.9615.4536.5673.29取大值10.99407.60173.3469.06206.7974.03113.84229.80最小值1.7740.0016.202.299.645.5119.6832.86标准差1.8078.3824.5910.7327.8510.4317.7330.94变异系数0

17、.440.510.630.620.680.670.490.42交 通 区平均值5.71360.0158.0562.21446.8217.6263.53242.85最大值30.131619.80920.841364.8516000.00142.50181.483760.82最小值1.6150.1015.3212.348.576.1922.0140.92标准差3.24243.3981.61120.222180.2711.7932.53384.78变异系数0.570.681.411.934.880.670.511.58公 园 绿 地 区平均值16.26280.5443.6430.19114.9915

18、.2960.71154.24最大值11.681024.9096.28143.311339.2929.10227.401389.39最小值2.7797.2016.319.0410.007.6026.8937.14标准差2.02235.8414.8422.68224.284.9745.84230.92变异系数0.320.840.340.751.950.330.761.504.3元素浓度的无量纲化处理在利用SPSS统计软件数据进行聚类分析的时候，因为单位不统一需要进行无量纲 化处理，我们采用均值化方法，即每一个变量除以该变量的平均值，即xixi标准化以后各变量的平均值都为1,标准差为原始变量的变异系

19、数。该方法在消除量纲 和数量级影响的同时，保留了各变量取值差异程度上的信息，差异程度越大的变量对综 合分析的影响也越大。4.4重金属元素间的相关性分析研究土壤中重金属的相关性可以推测重金属的来源是否相同 ，若重金属含量有显著的相关性，说明有相同来源的可能性较大，否则来源可能不止一个.我们用积差法来计算各重金属之间的相关系数，所谓积差法就是用两个变量的协方 差与两个变量的标准差的乘积之比： x _ x y _ y表2重金属元素间的相关系数重金属AsCdCrCuHgNiPbZnAs1Cd*0.2551Cr*0.189*0.3521Cu0.160 *0.397 *0.532 *1Hg0.0640.2

20、65 *0.1030.417 *1Ni0.317 *0.329 *0.716 *0.495 *0.1031Pb*0.290*0.660*0.383*0.520*0.298*0.3071Zn0.247 *0.431 *0.424 *0.387 *0.196*0.436 *0.494 *1由元素间的相关系数（见表2）我们对这些元素进行粗略的分组，大致分为以下两组:A组：Cr，Ni， CuB 组：Pb, Cd, Zn而对于As、Hg由相关系数表可见，其相关系数较小，我们认为相关参数小的元素间没 有关系，所以将其各自单独一组。下面我们建立回归模型图像验证它们之间的函数关系：显然Cr，Ni和Cu及Pb,

21、 Cd和Zn显示属于适度空间相关性，反映区域因素（土壤母质） 对其含量的影响较大，而As、Hg元素则属于低空间相关性，说明其受到人为因素（工 业布局施肥灌溉和土地利用方式等）作用较强。五、模型的分析、建立与求解5.1问题一5.1.1重金属的空间分布由附件中所给的数据，我们考虑将各采样点的坐标和重金属的浓度建立对应关系, 利用Matble软件画出等高线来体现该城区8种重金属的空间分布。图3城市土壤Cr的空间分x布。特征图4城市土壤的空间分X%特征从图中可以看出：该城市土壤中 As元素的分布没有出现明显的富集，整体有从西 向东递减的趋势（见图1）。说明人类活动对As元素的分布影响不大。所以可以推断

22、城市 土壤中这种元素主要是自然来源，另外它的浓度在中国土壤背景值范围内，这说明它的 含量可能主要受成土母质影响。该城市土壤中Cd元素的分布没有出现明显的富集，整体浓度偏差不大（见图 2）。对比数据可以看出，整个城市除边缘部分外 Cd的浓度都明显高于背景值的范围。可知该城市Cd污染很严重。该城市土壤中Cr和Cu两种元素含量的空间分布规律比较相似（见图3、图4）,表现 为,在西南部形成一个明显峰值，并且西部Cr和Cu的浓度远远超出背景值的范围。从整 体上看，两种元素的浓度在东部和中部都为零，显然西部高于中东部。该城市土壤中Hg和Zn两种元素含量的空间分布规律比较相似（见图5、图8）,表现 为,一个

23、峰值区出现在西南部，一个峰值区出现在中南部，另外一个峰值区出现在中部。 当然，两图也存在着不同之处。Hg元素除峰值区外，其他部分的浓度大都为零，而 Zn元素除峰值区外，还有整个的西部浓度远远超出背景值的范围，其余部分的浓度为零。该城市土壤中Ni和Pb两种元素含量的空间分布规律比较相似（见图6、图7），虽 然Pb比Ni多出一个明显峰值区，但可以看到两种元素的峰值区都分布在西南部。并且 两元素在西部的浓度明显高出背景值范围，而在东部的浓度都非常接近背景值的范围， 整体有从西向东递减的趋势。同时，土壤中重金属元素的空间分布还显示出各元素的异常分布区具有地理趋势的 相似性，指示其可能受共同的污染源影响

24、5.1.2不同功能区的污染程度为了求得各功能区的污染程度，我们建立了内梅罗多因子污染综合评价模型，我们 首先求得单项污染指数式为：PiCiSiSi为重金属i式中：Pi为区域重金属i的单项污染指数；Ci为重金属i含量的实测值;含量的起始评价值，其中起始评价值为所给重金属元素的背景值加上两倍的标准差，即Si =ai ;若Pi 1，则表示该区域受到污染。多项污染综合指数式为2 2P平均* P max式中：P综为综合污染指数-综合反映各重金属对区域土壤的不同作用； P平均为所有单项污染指数的平均值；Pmax为土壤环境中各单项污染指数中的最大值。根据模型分别计 算出该城区五个功能区的污染指数值， 然后与

25、内梅罗综合污染指数的分类标准 （见表3） 相比较得出五个功能区的污染程度。表3内梅罗综合污染指数的分类标准P综污染等级污染程度1P综 0.7安全清洁20.7 P综 1 .0警戒线尚清洁3XOCP综兰2.0轻污染超标42.0 P综兰3.0中污染土壤、作物受中度污染5P综 A3.0重污染土壤、作物受严重污染由Matlab软件计算得1800016CD012000-10000e呦JOOO2000表4各功能区的污染程度功能区生活区工业区山区主干道路区公园绿地区P综2.09049.27690.82426.46251.8688污染程度土壤、作物受 中度污染土壤、作物受 严重污染尚清洁土壤、作物受 严重污染超

26、标由表4可得污染程度排序为：工业区 主干道路区生活区公园绿地区山区5.2问题二5.2.1无量纲化处理在利用SPSS统计软件数据进行聚类分析的时候，因为单位不统一需要进行无量纲 化处理，我们采用均值化方法，即每一个变量除以该变量的平均值，即xiXi -，X标准化以后各变量的平均值都为1,标准差为原始变量的变异系数。该方法在消除量纲 和数量级影响的同时，保留了各变量取值差异程度上的信息，差异程度越大的变量对综 合分析的影响也越大。然后取其平均值进行对比。5.2.2重金属污染的原因分析对附件中所给数据进行分析，采用均值化方法无量纲化处理后，我们选取8种重金 属元素的平均值，用EXCELLO件画出其相

27、对含量的柱状图：图5各功能区重金属元素无量纲化后数据比较As Cd(2) Cr(3) Cu Hg(5) Ni(6) Pb(7) Zn(8)8种重金属元素经无量纲化处理之后可以相互比较，由柱状图可以看出：生活区中 Cr含量最大,Zn,Pb含量次之；工业区中Cu, Hg含量较高，且其他重金属也相对较多； 山区中Ni含量最大，Cr，As次之，且8种重金属含量都相对较少；交通区中各重金属 含量都比较集中，Hg最多；公园绿地区中As含量较多，Pb,Cd含量次之。我们还可以横向来比较8个元素在五个功能区中的相对含量：As、Cr、Ni在各功能 区的含量相对均匀，而 Cu Hg在工业区中的含量与其它区差别明显

28、； Cd Zn、Pb、Cu Hg在工业区中的含量和山区相比差异较大，说明它们受人类活动的影响相对明显。重金属具有富集性，很难在环境中降解且 Hg Cd Cr、Pb具有很强的生物毒性。 由于工业生产产生大量废物，汽车尾气排放，汽车轮胎磨损产生大量含重金属的有害气 体和粉尘，Cu、Hg、Zn主要在工业区和交通区如公路、铁路等交通设施的两侧富集， 随时间的推移，工业区、交通区的土壤重金属具有很强的叠加性。 Pb Cd As、Cr主要在生活区富集，生活区内由于含重金属的废弃物堆积，施用含 Pb Cd As等的农药，不合理的施用化肥以及大量使用农用塑料薄膜都会造成重金属污染， 综合风向等多种因素的影响可能会发生迁移且富集，从而造成公园绿地区的重金属污染。Ni元素在各个功 能区中的分布相对均匀，山区分析的 8种重金属中Ni含量最多，所以说工业区及生活 区污染的富集和迁移，都有可能