城市表层土壤重金污染问题建模——第二组.docVIP

城市表层土壤重金属污染问题建模 摘 要 本文对某城市城区土壤地质重金属含量进行分析,首先利用数理统计知识和matlab数学软件研究了生活区、工业区、山区、主干道路区及公园绿地区等表层土壤中重金属中As 、Cd、Cr、Cu、Hg、Ni 、Pb、Zn的含量，并应用三维插值得到主要重金属元素在该城区的空间分布图，再采用Mul1er 地积累指数法对其污染状况进行评价，得到结论： （1) 土壤的污染程度为工业区交通区生活区公园绿化区山区，其中工业区部分地区出现强污染。 (2) 不同功能区表层土壤污染的重金属种类存在差异,工业区和交通区以 Cd, Hg为主,山区重金属含量普遍很低，呈无污染。 然后通过八种金属元素主要来源、分布情况及污染特征，得出重金属污染的主要原因： （1）人类活动，包括工业污染、交通尾气排放、城市堆放废弃物、污水灌溉。 （2）与海拔高度有关，海拔越高，其污染共有八个方程，当分别取最大值时，可得到对应的坐标值，即对应重金属的污染源，结果是共有八个污染源。在平面坐标上As与Cd的污染源为（13855，18449），Cr、Cu与Zn为（13855，8260）和Hg、Ni与Pb为（28654，8260）。 本题得出重金属污染情况与特点，能很好应用于城市表层土壤重金属污染分析中，为城市土壤重金属污染防治提供科学依据。 关键词： 重金属污染 Mul1er地积累指数法 多元线性回归分析 污染源matlab数学软件 插值 一 问题重述 城市土壤地质环境查证以及对相关数据处理，是城市环境质量评价的重要依据。现对某城市城区土壤地质环境进行调查, 将城区分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、……、5类区，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10 厘米深度）进行取样、编号，记录采样点的位置，并获得了每个样本所含的多种化学元素的浓度数据。再者，按照2公里的间距在远离人类活动的自然区取样，其值作为该城区表层土壤中元素的背景值。 利用附件1列出了采样点的位置、海拔高度及其所属功能区等信息、附件2列出了8种主要重金属元素在采样点处的浓度、附件3列出了8种主要重金属元素的背景值的条件 要求你们通过数学建模来完成以下任务： (1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。 (2) 通过数据分析，说明重金属污染的主要原因。 (3) 分析重金属污染物的传播特征由此建立模型，确定污染源的位置。 (4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？模型优化有了这些信息，如何建立模型解决问题？ 二 基本假设 8种主要重金属元素在该城区的空间分布暂且定义为每种重金属元素浓度在该城区中的分布； 影响重金属元素在采样点处的浓度的主要因素为采样点的位置（x1与x2）、海拔高度（x3），不考虑其他项目； 设表中的数据是真实可靠，忽略实验测量误差。并具有广泛代表性； 在模型的求解过程当中，考虑到各个变量之间数据差距较大，所以采用了数据归一化处理（MATLAB函数为X=zscore(A)） 各监测指标之间不相关，互相独立； 重金属传播过程中不考虑人为的大量移动； （7） 符号 含义 i=1，2，3；分别代表X轴，Y轴和海拔高度； 8种重金属元素的浓度 对重金属元素浓度起影响的权重，为实数范围内的常量； 变量因素个数（i=3）； 线性方程个数（n=319）； 3.2名词约定 地积累指数：又称Muller 指数，用于研究沉积物及其他物质中重金属污染程度的定量指标，其表达式为：Igeo=log2［Cn/(k ×Bn)］，式中，Cn 是元素n 在沉积物中的含量；Bn 是沉积物中该元素的地球化学背景值；k 为考虑各地岩石差异可能会引起背景值的变动而取的系数(一般取值为1.5)； 四 模型的准备 4.1城区各功能区—8种主要重金属元素平均值 城区各功能区不同重金属平均浓度 　含量 As (μg/g) Cd (μg /g) Cr (μg/g) Cu (μg/g) Hg (μg/g) Ni (μg/g) Pb (μg/g) Zn (μg/g) 生活区 6.27 289.96 69.02 49.40 93.04 18.34 69.11 237.01 工业区 7.25 393.11 53.41 127.54 642.36 19.81 93.04 277.93 山区 4.04 152.32 38.96 17.32 40.96 15.45 36.56 73.29 交通区 5.71 360.01 58.05 62.21 446.82 17.