城市表层土壤重金属污染程序.docVIP

城市表层土壤重金属污染分析 摘要 土地是人类赖以生存的最基本自然资源之一，随着人类活动的日益频繁，工业、城市污染的加剧和农用化学物质的增加，土壤重金属污染愈发严重。本文主要讨论如何应用对土壤地质环境异常的查证资料来评价城市环境质量及污染源分布情况，并研究人类活动影响下的城市地质环境的演变模式。 针对第一问，利用插值法把所给地形的离散数据平滑化，在所给数据x,y中分别插入100个数据，并求出对应的海拔高度，通过matlab画出该城市3D地形图。依照此法依次处理八种重金属元素浓度采样数据，画出其等高线图。二者综合即可看出重金属元素在该城区的空间分布；通过单因子指数法和内梅罗指数法，可以依次求得各重金属元素在五类区域的污染指数和各区域综合污染指数，对照国家土壤污染标准表，即可得到不同区域重金属污染程度。 针对第二问，为分析重金属污染的主要原因，可根据第一问得出的不同区域内各重金属污染程度表，依次考虑该城区五类区域的污染情况。通过五类区域内部污染原因分析，再考虑该城区地形分布，进行区域间相互污染原因分析，即可得到该城区重金属污染的主要原因。 针对第三问，为求出各重金属的污染源，运用元胞邻居模型思想，我们假设重金属污染物朝各个方向扩散且相对均匀，可以认为污染源的污染程度总是高于周围地区的污染程度。根据插值法，利用matlab算出对应的污染浓度，由此可得各重金属元素在该城区的污染浓度矩阵（100*100）。再利用C++编写搜索程序求得所有污染浓度大于周围采样点污染浓度的点（极大值点），作为污染源的备选点。为了简化模型，我们认为边界点总不是污染源。由于某些地区可能污染浓度普遍偏小，该处极大值点不能作为污染源。因此，需要对照国家土壤环境质量标准，分别对每种重金属元素污染物浓度分级。选取其中适当的点作为各重金属元素的污染源。 关键字：MATLAB 插值 单因子指数法 内梅罗指数法 元胞 一 问题重述 现把某城区分为生活区、工业区、山区、主干道路区及公园绿地区。将城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10 厘米深度）进行取样、编号，并记录采样点的位置。利用化学仪器获得每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，作为表层土壤中元素的背景值。 附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。 现有以下问题： (1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。 (2) 通过数据分析，说明重金属污染的主要原因。 (3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。 (4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？ 二 模型假设 城市处于内陆区，周围没有海洋和大型河流； 重金属污染物朝各个方向扩散且相对均匀； 短期内该城市没有重大的地质变化； 短期内该城市各功能区没有重大变化。 污染源扩散程度与地形无关。 三 符号说明 重金属污染物i的单因子指数 重金属污染物i的实测值的算术平均值 重金属污染物i的背景值的最大值 某区域各重金属单因子指数的平均值 某区域各重金属单因子指数的最大值 重金属污染因子权重 四 问题分析 对于第一问，考虑插值法把所给地形的离散数据平滑化，通过matlab画出该城市3D地形图。依照此法依次处理八种重金属元素浓度采样数据，画出其等高线图。二者综合即可看出重金属元素在该城区的空间分布；通过单因子指数法和内梅罗指数法，可以依次求得各重金属元素在五类区域的污染指数和各区域综合污染指数，对照国家土壤污染标准表，即可得到不同区域重金属污染程度。 对于第二问，为分析重金属污染的主要原因，可依次考虑该城区五类区域的污染情况。通过五类区域内部污染原因分析，再考虑该城区地形分布，进行区域间相互污染原因分析，即可得到该城区重金属污染的主要原因。 对于第三问，为求出各重金属的污染源，运用元胞邻居模型思想，我们假设重金属污染物朝各个方向扩散且相对均匀，可以认为污染源的污染程度总是高于周围地区的污染程度。根据插值法，利用matlab算出对应的污染浓度，由此可得各重金属元素在该城区的污染浓度矩阵（100*100）。再利用C++编写搜索程序求得所有污染浓度大于周围采样点污染