2022年全国大学生数学建模竞赛A题——城市表层土壤重金属污染分.docxVIP

城市表层土壤重金属污染分析 摘要 针对问题一，题目要求给出8种重金属元素在该城区的空间分布。本文以采样点所在位置的坐标分别为x、y轴，以该采样点各种元素的浓度为z轴，利用matlab样条工具箱进行“v4”内插，作图得到8种重金属元素的水平面分布图，并将三维图像转化为等值线图，从图4-图11可以直观地看出各种重金属元素在空间的分布情况。结合图像，采用单因子指数法和内梅罗污染指数法对城区内不同区域重金属元素的污染程度进行评估，得到污染程度从高到低依次为：工业区，主干道路区，生活区，公园绿地区，山区。 针对问题三，题目要求建立模型确定污染源的位置。首先从纵向，横向以及考虑自然因素的情况下分析重金属污染物的传播特征，详细情况如模型分析所示。针对这些传播特征，主要考虑对传播路径起主要影响的纵向和横向因素，运用类比的方法，把重金属污染物的传播方式等效为电磁场理论中电场能量的传递方式，即用通量的大小来体现某特定采样点是否为源的可能性。同时给出了线通量的定义及在本问题中的表达式，得到各元素线通量分布及等通量线，并得到各种元素可能的污染源坐标，如表5所示。 针对问题四，题目要求在收集到更多信息的前提下，建立模型以更好地研究城市地质环境的演变模式。本文认为，需要收集的信息还有：气温、降水量、蒸发量、水体污染情况、地面坡度、岩性特征、历史地震活动情况、工程地质。结合在问题一中已获得的土壤污染综合指数，采用主成分分析法建立模型，得到影响城市地质环境的主要因素，从而通过研究这些影响地质环境的主要因素获得城市地质环境的演变模式。 关键词：关键词：重金属元素污染 内梅罗指数 灰色关联分析 线通量模型 主成分分析 一、问题重述 按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、……、5类区，不同的区域环境受人类活动影响的程度不同。 附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。 现要求你们通过数学建模来完成以下任务： (1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。 (2) 通过数据分析，说明重金属污染的主要原因。 (3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。 (4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？ 二、符号说明 三、模型假设 1、假设所采用的采样方法科学合理，所得数据准确。 2、假设污染物的扩散只与浓度有关，且顺浓度梯度方向扩散，无外力参与。 3、假设不同功能区对应一种级别的污染程度，不存在具有相同污染程度的功能区。 4、假设数据是在污染物扩散稳定的情况下测得的。 四、问题分析 4.1 对重金属元素的空间分布和在不同城区的污染程度进行分析 8种重金属元素在空间的分布，用等高线图可以直观地表示出来。本文以采样点所在位置的坐标为x、y轴，用暖色系（灰度高）表示重金属含量高，以冷色系（灰度低）表示重金属含量低。由于题目所给采样点为离散的，故本文在画图前对数据进行了“v4” [1,2]内插处理，利用matlab样条工具箱得到8种重金属元素在该城区的8个空间分布图。 8种重金属元素在该城区各区域的分布各不相同，要得到该城区内不同区域重金属的综合污染程度，还需要一个指标将8种不同的重金属元素的污染程度综合起来。内梅罗污染指数[3]是一种兼顾极值或称为突出最大值的计权型多因子环境质量指数，能够合理有效地评价环境的综合污染程度，内梅罗指数越大说明该地区污染越严重。 4.2 对重金属污染的主要原因进行分析 重金属污染是指由于人类活动，土壤中的微量有害元素在土壤中的含量超过背景值，过量沉积而引起的含量过高，统称为土壤重金属污染[4]。在不同的功能区，重金属的由来有所不同。本文采用灰色关联分析法[5]，通过对各功能区各重金属元素进行关联度分析,得到对每个功能区重金属污染贡献最大的几种元素，并针对功能区的不同特点分析贡献最大的几种重金属元素的由来，从而可以得到各功能区中重金属污染的主要原因。 4.3 对污染源位置确定进行分析 金属污染物的传播路径受到多种因素的影响，如风力，水流，海拔高度……由某一污染源散发出的重金属污染物，主要有以下传播特征： （1）从纵向的角度上看，由于重力的作用，重金属污染物从高处向低处传播，在没有污染源进行补充的情况下，污染物浓度随海拔的减小而下降； （2）从横向的角度上看，重金属污染物从污染源向四周呈放射性传播，在同一海拔高度，污染物浓度随与污染源的距离的增加而下降； （3）考虑自然因素的影响，由于重金属污染物还存在于大气或者水流中，其传播方向与风向和河流等流