2011数学建模国赛A.doc

城市表层土壤重金属污染分析 摘 要 城市表层土壤重金属污染严重影响人的生活和生产，在已知319个样本的坐标以及各重金属污染物浓度情况下，按照各问题的要求，我们建立相应的模型来对具体问题进行求解。 首先为了了解该城市的地形，我们使用了MATLAB软件绘制各分区分布地形图,根据区域化变量理论，运用克里金插值法对未采样点的区域化变量的值进行线性无偏最优估计，绘制出各污染元素的拟合分布图。在分析不同区域受污染程度时，以题中所给的元素背景值以及国家土壤环境质量准则为基础制定污染等级评判标准，据此进行单因子指数评价以及模糊数学评价分析，得出每个区域内主要污染元素以及污染程度。 接下来，为了分析污染的原因，我们将单因子指数的结果与聚类分析算法的结果进行比较，得出生活区内铬、铜、镉、锌含量较高；工业区内的主要污染元素为镉、铬、汞、镍、铅、锌；山区各重金属含量均较低；交通区内镉、铜、汞污染较重；在公园绿地区内，镉、汞、铅、锌的污染较为严重。 更进一步，为了确定每一种重金属元素的来源区域，我们结合重金属污染物在土壤内易沉积，不易传播的特性以及当地海拔高度，通过主成分分析法和相关分析法计算出污染物的所在区域，并与surfer的绘图点进行比较，发现两者有很好的吻合度，找到了Cd的污染主要来自工业区和交通区；Cr污染主要来自交通区和生活区；Cu污染主要来自交通区和工业区；Hg污染主要来自交通区；Ni污染主要来自交通区；Pb污染主要来自生活区和工业区；Zn污染主要来自交通区，生活区，工业区。对于污染物传播特性的分析，我们首先用SPSS对各区域内各污染因子进行正态检验，然后对于正态分布和非正态分布的因子我们各给出一种函数描述其传播特性。 对于问题四，我们提出建立城市污染物数据库，记录城市污染物随时间和空间的变化情况，有利于对未知点进行时间和空间的预测。 关键字：重金属污染、克里金插值法、多因素分析法、聚类分析、主成分分析法 问题重述 城市是人类社会不断发展的产物,是人类重要的生活环境,也是人类生产生活的主要场所. 随着经济和社会发展速度的不断加快,城市环境问题,尤其是重金属作为一个威胁人类生存和城市发展的重要的城市环境污染物,引起了学者、专家们的关注.随着工农业发展和城市化进程的加剧, 通过农业灌溉、工业排放、交通运输、市政建设和大气沉降等造成城市土壤重金属的污染日趋严重。由于重金属有污染物在土壤中移动性差、滞留时间长、不能被微生物降解的特点以及具有极强的生物富集性,能通过生物链转移,直接或间接地危害人体健康，因此进行针对土壤重金属污染程度的研究和原因的探析对于保护人类健康和促进经济的可持续发展都具有十分重要的意义。 本题在5个不同的城市功能区中随机选取319个采样样本，并应用专门仪器测试分析，获得了每个样本所含的8种主要重金属污染元素的浓度数据。下面我们结合数据利用数学建模，展开对城市表层土壤重金属污染的分析，以期解决下列问题： 问题一，确定8种主要重金属元素在该城区的空间分布，并分析不同区域重金属的污染程度。 问题二，通过数据分析，说明重金属污染的主要原因。 问题三，分析重金属污染物的传播特征，对此建模，确定污染源的位置。 问题四，分析建立模型的优缺点，为更好地研究城市地质环境的演变模式，还需要那些信息，并且如何建立模型解决问题。 二、问题分析 从总体上看, 城市土壤中重金属含量要明显高于郊区及远离城市的农田土壤的含量, 但是近年来, 随着城市环境的综合治理, 部分企业向城郊搬迁, 加之乡镇企业蓬勃发展, 小型乡镇企业数量剧增, 使得重金属污染出现了明显的郊区化趋势。另外,城市重金属污染还表现出个别区域的某些重金属元素的严重富集( 如加油站、冶炼厂附近等) ，以及交通干线两侧污染较为严重等特征。城市土壤受人为扰动较大, 挖掘、搬运、堆积、混合和大量的废弃物填充, 使土壤结构和剖面发育层次十分混乱, 所以重金属在土壤剖面特征不同于自然土壤或表层耕作的农业土壤, 分布无规律性。 由于所给采样点较多，而且各采样点在三维空间随机分布，很难确切了解到生活区，工业区等的分布情况，我们首先用MATLAB画出各采样点的平面和空间坐标，达到对各区分布的初步了解，有利于对下面各问题的分析。图形如下： 图1 各功能区分布的平面图 图2 各区分布的立体图 由图1图2可以判断出山区主要分布在远离生活区和工业区的高纬度地区；生活区与工业区相邻，与交通区交叉分布；公园绿地区较少，分布在生活区和交通区之间。由于各区相近或者交叉，所以它们的污染物极有可能也存在相关关系。 由于实际情况下，重金属的分布规律还受时间的影响，但在这里我们不考虑重金属的分布随时间的变化，只研究在空间的分布规律。 城市表层土壤重金属污染分析的问题，要对给定的有限的数据进行分析与处理，来为后续各个问题的