搜索
A题城市表层土壤重金属污染分析随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、……、5类区，不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。现要求你们通过数学建模来完成以下任务：(1)给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。(2)通过数据分析，说明重金属污染的主要原因。(3)分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4)分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？1城市表层土壤重金属污染分析摘要本文主要研究重金属对城市表层土壤污染的问题,我们根据题目所给定的一些数据和信息分析并建立了扩散传播模型、自然沉降模型、对比模型和转换模型解决问题。针对重金属的空间分布问题，先拟出该城区地势图，根据所给数据绘出该地区的三维地势及采样点在其上的空间分布图。再利用MATLAB散乱插值法得到8种重金属元素的空间分布。其次，通过单因子污染指数法和内梅罗综合指数法两种方法，得出城区内不同区域重金属的污染程度：工业区交通区生活区公园绿地区山区。再分析各功能区内8种重金属污染所占比重，结合各重金属来源得出结论：该城区的重金属污染主要原因在于工业释放及交通排放。最后，分别建立大气扩散模型、土壤扩散模型、水体扩散模型三个传播模型，自然沉降模型，结合酸雨对重金属元素的影响，对重金属在土壤中的浓度进行了时间和空间上的拟合。当存在n个污染源时，记000(,,)jjjxyz为第j个污染源坐标，jM是瞬时点源的质量，(,,,)cxyzt是t时刻(,,)xyz处的重金属污染物的浓度，得到与重金属元素污染源位置有关的的偏微分方程：22200031221()()()1(,,,)exp2444(4)()njjjjjxyzxyzMxxyyzzCxyztKttEtEtEtEEE其中,,,xyzEEEk为相关系数。根据取样点重金属浓度变化确定污染源位置,发现城市表层土壤重金属污染主要是由工业区、交通区和生活区的污染造成。为了对重金属浓度和客观因素优化推广，对城市地质环境的演变做更准确的分析和预测，需要搜集不同历史时间段土壤污染程度、土壤浓度的PH值、地下水重金属的浓度、植物中重金属浓度等多重信息。通过指数法评价金属元素在土壤、PH、地下水、植物中的污染程度，采用层次分析法综合评价土壤污染程度，说明地质环境的演变模式。再根据不同历史时间段土壤污染程度，说明城市的地质环境变化模式。关键词：散乱插值法重金属扩散方程单因子污染指数法内梅罗综合指数法2一、问题重述目前，社会经济发展迅速，人口数量不断增加，环境污染现象日显突出，尤其重金属对土壤的污染更受广泛关注。土壤状况直接影响着动植物的生长和安全，甚至通过食物链进入人体，导致一些慢性疾病的发生。对于独立的系统来说，人们的生活和生产将会给环境和土壤造成污染，而且，每一个区域的功能不同，对环境和土壤的污染程度也不同。根据现有的区域——生活区、工业区、山区、交通区、公园绿地区多个取样点的坐标，海拔高度以及8种重金属元素的浓度，参考自然区取样得到的背景值，分析以下问题：1)给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。2)通过数据分析，说明重金属污染的主要原因。3)分析重金属产生后的传播特性，由此建立模型确定污染源的位置。4)考虑重金属浓度和客观因素优化模型并推广，以便对城市地质环境的演变做更准确的分析和预测。二、问题分析1)重金属元素的空间分布及不同区域重金属污染程度根据题中所给采样点的坐标和各重金属的浓度，可用散乱数据插值的方法获得各重金属污染物浓度的空间分布，再用MATLAB软件绘制出重金属在该区的空间分布图及8种重金属的浓度等高线在该地区三维地形图曲面的投影图。可以通过单因子污染指数法和内梅罗综合指数法相互验证，得出城区内不同区域重金属的污染程度。2)重金属污染的主要原因在上题的结论基础上，依据各种重金属在不同区域的浓度大小，通过分析比较，并考虑各区域的特点，可以分析出重金属的主要来源。3)确定污染源的位置由于重金属污染物主要是通过大气、土壤和水体进行传播，同时可能受到酸雨淋溶和冲刷的影响，自然沉降后吸附在土壤中。故受到五方面因素的影响，且影响因素间相互关联，相互影响。可以建立多个传播模型，如大气扩散模型、土壤扩散模型、水体扩散模型，以及自然沉降模型，结合酸雨对重金属元素的影响，对重金属在土壤中的浓度进行时间和空间上的拟合。根据重金属浓度变化得到取样点与污染源的位置关系，得出污染源的位置。34)重金属浓度和客观因素优化模型在研究污染物传播的过程中，本题只给出了一个时间点上的数据,信息量明显不足,应搜集一段时间内的数据。由此建立模型可推出不同时间不同地点污染物的浓度。为确保模型的准确性，还需考虑到土壤浓度的PH值、地下水重金属的浓度、植物中重金属浓度等环境信息，以建立更真实的模型，对城市地质环境的演变做更准确的分析和预测。三、模型假设1)取样点的数据较好的反映了该地区的污染物浓度。2)重金属元素在无穷空间扩散，不计风力影响。3)重金属元素的扩散符合扩散定律，即单位时间通过单法相面积的流量与它的浓度梯度成正比。4)重金属元素在土壤和水中化学反应均匀。5)重金属元素在土壤中处于相对稳定的状态。6)不考虑不同区域内大气、土壤、水体等自然环境差异。四、符号说明符号变量符号说明iP第i个取样点污染指数iC第i个取样点污染浓度测量值jS第j种重金属元素评价标准P内梅罗污染综合指数,,,Cxyzt空间任一点,,xyz在t时刻的污染物浓度K污染物扩散系数1K悬浮泥沙对重金属的吸附速率2K重金属在悬浮泥沙中的解析速率dK重金属离子在土壤中固态和液态的分配系数Q污染物的总污染量4符号变量符号说明D扩散系数R重金属离子在土壤中迁移的阻滞因子N混合沙的吸附量1N悬移质平均吸附量土壤介质的密度泥沙干容重h断面平均水深五、模型的建立及求解1.重金属元素的空间分布及污染程度模型1.1.重金属元素的空间分布：根据已知数据，我们运用MATLAB软件将该地区的三维地势及采样点在其上的综合空间分布图绘制如图1-1所示。图1-1城区地形及各功能区空间分布图图5在此图中，我们能清晰的分出生活区、工业区、山区、交通区、公园绿地区。将生活区、工业区、山区、交通区、公园绿地区用i表示。1,2,3,4,5i根据题目要求，为给出8种主要重金属元素在该城区的空间分布，我们采用等高线的分析方法把各金属的浓度投影到平面坐标，用MATLAB仿真软件分别作出8种重金属元素的空间分布图。图1-2城区地形及各功能区空间分布俯视图图1-3As元素的空间分布图图1-4Cd元素的空间分布图6从图中可以清晰地看出8种主要重金属元素在该城区各区域的空间分布。1.2.重金属元素污染程度：我们利用现有资料以及查找到的信息，利用单因子污染指数法和内梅罗综合指数法来评价城区内不同区域金属的污染程度。1.2.1.单因子污染指数法单因子污染指数法是国内通用的一种重金属污染评价方法，其计算公式如下：图1-9Pb元素的空间分布图图1-10Zn元素的空间分布图图1-5Cr元素的空间分布图图1-6Cu元素的空间分布图图1-7Hg元素的空间分布图图1-8Ni元素的空间分布图7=iijPCS（1.2.1）我们采用《土壤环境质量标准》（GB15618-1995）中国家一级标准作为环境质量的标准（见表1-1）。利用附录二中的数据，我们使用Excel软件做出了城区内不同区域内重金属污染物的单因子污染指数（见表1-2）。Pi1表明未受污染，Pi1表示已经受到污染，Pi数值越大，表明受到的污染越严重。表1-1土壤环境质量标准项目土壤PH值一级自然景值二级PH﹤6.5PH=6.5～7.5PH＞7.5三级PH＞6.5As(mg/kg)1530252030Cd(mg/kg)0.200.300.300.601.00Cr(mg/kg)90150200250300Cu(mg/kg)3550100100400Hg(mg/kg)0.150.30.51.01.5Ni(mg/kg)40405060200Pb(mg/kg)35250300350500Zn(mg/kg)100200250300500表1-2城区不同区域内重金属污染物的单因子污染指数单因子污染指数一类区二类区三类区四类区五类区As0.41803030.4834260.2696060.380.42Cd1.449806821.9655560.761.811.402714Cr0.766871210.590.430.6461580.484844Cu1.413.640.4947791.7862570.86262Hg0.624.2823690.273042.9876160.766611Ni0.458556820.4952920.3863450.4403210.382243Pb1.974467532.658311.0444551.8211351.734531Zn2.370086362.7792750.7329422.4409941.54最大单因子污染指数2.370086364.2823691.0444552.9876161.734531平均单因子污染指数1.183701222.1126930.5494561.5386970.9491968根据表1-2的数据，利用Exce软件做出8种重金属元素在五类区中的单因子污染指数柱状图：由上图可见，在一类区中，Cd、Cu、Pb、Zn的单因子污染指数均高于1,表明一类区已受污染；在二类区中，Cd、Pb、Cu、Hg、Zn的单因子污染指数均高于1，表明二类区已经受到污染；在三类区中，仅Pb的单因子污染指数略高于1，其为受污染区；在四类区中，Cd、Pb、Cu、Hg、Zn的单因子污染指数均高于1，表明其为受污染区；在五类区中，Cd、Pb、Zn的单因子污染指数均高于1，所以五类区均为受污染区。综合以上分析，可以得出，城区的不同区域均受到不同程度的污染。1.2.2.内梅罗污染综合指数法单因子指数可以判断出环境的主要污染因子，但环境是一个复杂的体系，环境污染往往由多个污染因子复合污染导致，故单因子指数法不能全面、综合的反映土壤的污染程度，于是采用内梅罗污染综合指数法进行分析，计算公式如下：22111max2niiiPPPn（1.2.2）根据附件中的数据及单因子污染指数，做出城区不同区域的重金属污染内梅罗综合指数，见表1-3及图1-12。表1-3城区内不同区域的重金属污染内梅罗综合指数内梅罗污染综合指数一类区二类区三类区四类区五类区1.87329363.3765480.8345022.3762831.40图1-118种重金属元素在五类区中的单因子污染指数图9内梅罗综合污染指数全面反映了各污染物对土壤污染的不同程度，同时充分考虑高含量物质对本土壤的影响，结合相关标准，将土壤污染等级进行划分（见表1-4）。表1-4土壤污染等级划分等级P污染程度污染水平1≤