2011国赛A.doc

城市表层土壤重金属污染建模与分析 摘要 本文通过对样本值的统计分析，画出了土壤重金属污染物的离子等浓度线分布图，利用地质累积指数法对各区域污染等级做出了较全面的评价，分析出各种重金属污染的主要原因，并建立了重金属离子浓度空间分布的数学模型，对城市环境的保护研究有较高的参考价值。 对于问题一，通过对城区土壤采样数据分析，利用SURFER地图软件绘制浓度随空间坐标变化图。其主要步骤是先用克里格插值法划分网格线，再用半方差函数对城市重金属离子做空间等浓度分布曲线。对于污染评价，本文首先想到内梅罗综合评价的方法，虽然它能准确反映各污染物对土壤的不同作用，但无法判断人类活动对重金属元素富集程度的影响。为获得人类活动对城市土壤的影响，文中主要采用地质累积指数法及其污染评级标准对城区土壤中的重金属污染程度进行评价。 对于问题二，按功能区对样本值进行划分，利用统计学相关知识，凭借Excel等软件统计出各种重金属离子在各个功能区的平均值、样本标准差以及各种元素间的相关系数，结合空间离子等浓度分布图得出重金属离子污染的主要原因有，工业废水、废气和废渣排放，煤、石油、天然气燃烧和矿石的焙烧，汽车废气的排放以及Ni-Cd电池废弃物处理不当等。 对于问题三，综合菲克第二定律、对流方程和机械弥散方程建立起一个二元二阶微分方程，在不考虑离子浓度随时间变化的情况下，可将此二元二阶偏微分化简为齐次偏微分方程。由于二元二阶偏微分方程无法求出解析解，于是采用有限差分法求出近似解。算法思想是，先将城区分别沿x，y方向等距离划分，形成矩形网格，然后用网格结点浓度随位置的变化反映离子浓度的空间分布，城区的边界离子浓度可以认为达到背景值。于是对所有内部点建立线性方程组，求出齐次解，最终利用样本值确定常系数，得到各个网格点的数值。最后，利用离子等浓度线图对区域大致分块，要求每块区域中尽量出现一个极值，求出每块区域中的极值，可以认为这个极值为理论上的污染点。另外，还可以利用切面法建立模型，找到污染源位置。其原理是，考虑到重金属离子运移的各向同性，建立离子浓度和到污染源点的距离的函数关系模型，利用拉普拉斯变换和泰勒级数求解出以污染源为原点的近似解析解，最后利用样本值估算污染源点的位置。 对于问题四，收集污染物浓度随时间变化的数据，建立起污染物浓度与时间和空间关系的数学模型。求解模型的过程中，采用了时间步长法和空间步长法将时间和空间区域剖分为若干网格，利用交替方向隐式格式的差分方法，用离散解决定解问题，得到在整个区域上的近似解，具有较高的二阶精度。此模型还可以作为一个预测模型，加以利用可以对污染物的扩散做出有效的预测，进而达到减缓污染物扩散的目的。 关键词：克里格插值，地质累积指数法，菲克定律，扩散方程，拉普拉斯变换，有限差分法 一、问题重述 表示土壤中重金属离子浓度 采样时刻的浓度变化率 、 表示土壤中重金属离子沿x,y方向上的对流速度 、 表示土壤中重金属离子沿x、y方向的弥散系数 、 表示空间沿、方向的距离步长 R 表示土壤的吸附系数 四、问题分析与模型的建立 4.1 八种重金属在该城区的空间分布 图1 SURFER软件处理数据截图 Kriging法是一种最佳空间局部估计方法,作为地统计的主要分析方法,综合考虑变量的随机性和结构性，SURFER将样本值处理，所得半方差函数, 对采样区域土壤中的主成分因子进行普通克里格方法插值, 并作出重金属的等浓度线图。同时将样本值按功能区在同一坐标系上做分类张贴图。 图2 不同区域地形（等高线对应不同的海拔） 图3 As空间分布图 图4 Cd空间分布图 图5 Cr空间分布图 图6 Cu空间分布图 图7 Hg空间分布图 图8 Ni空间分布图 图9 Pb空间分布图 图10 Zn空间分布图 从这八种元素在该城区的空间分布的图中可以看出，每种重金属在城区的分布并不唯一，每种元素离散地分布在城区的不同区域。其中以该城市西南部存在重金属元素，呈现出集中分布的趋势。 4.2 不同区域的污染程度分析 对于污染评价指标，首先想到的是单因子分析综合污染评价即内梅罗综合污染指数评价法，其公式为： （1） 其中表示综合污染指数，表示被评价金属离子污染指数的平均值，表示各种金属离子污染指数中最大值。 但是采用内梅罗综合污染指数评价方法对土壤中重金属污染进行评价，这些方法均能准确反映各污染物对土壤的不同作用，并突出高浓度污染物对环境质量的影响，但无法判断人类活动对重金属元素富集程度的影响。为获得人类活动尤其是工业活动对城市土壤的影响，文中主要采用地质累积指数法[1]对城区土壤中的重金属污染程度进行评价。 地质累积指数法是德国海