在问题一中，我们根据附件中给出的数据，运用Matlab数学软件画出了8种元素在城区空间的分布图。查阅相关资料结合分布图，先计算出各区域8种金属的污染程度，再考虑几种金属元素在污染中所占权重，得出各区域重金属的综合污染指数。最终得出该城区内不同区域重金属污染程度，各个区域均为极强的生态危害，这其中工业区最为严重，其后依次为公园绿地区、山区、主干道路区，污染程度最低的是生活区。
在问题二中，通过已知的数据，先画出城区地形图与分布图。然计算全部地区中各元素的平均值、标准差、变异系数、超标率、背景值。结合一问中得到的结果，分析得到的数据，表明重金属污染的主要原因有工业区污水废气排放和生活区的生活垃圾以及装修的涂料器材。公园绿地区、山区、主干道路区等地污染受工业区影响较大。公园土壤和灰尘中重金属的积累主要由人类活动造成。山区Cr,Cd,Pb等元素的污染是由于使用化肥和采用耕种方法不当使得农作物土壤受到了重金属铬和铅的污染。
在第三问中，区域土壤重金属污染具有普遍性、隐蔽性、表聚性、不可逆性等特点。属于分布式参数系统运用非线性优化及遗传算法建立模型，加以数学变换，得到最优解即污染源位置。例（4000,4000）点附近的区域是Cr、Hg、Zn、Ni和Pb的污染源，属重度污染，影响范围大。（22300,10500）点和（13800,9600）点，属轻度污染，影响范围较小。
根据题中所述，第四问为研究城市地质环境的演变模式，我们在原有模型的基础上考虑岩石圈、水圈、大气圈、地势分布对地质环境的影响，收集其数据。用同位素标记法最终建立数学模型，得到重金属污染与土壤的四个重要特性（即土壤紧实度、土壤pH值、土壤有机质以及土壤微生物与的关系）呈相关性的结果。
关键词背景值评价法、沉积法、Matlab、权重比较、遗传反演算法 问题的重述 环境污染与城市经济的快速发展和城市人口的不断增加有着密切的关系，人类活动对城市环境质量的影响日显突出。所以怎样对城市土壤地质环境异常做出有效查证，以及如何应用查证获得的海量数据资料开展对城市环境质量做出评价，以研究人类活动影响下城市地质环境的演变模式等问题，日益成为人们关注的焦点。随着现代经济的高速发展,重金属污染物的排放量与日剧增,环境污染和生态破坏越来越严重,危及到人类的健康。重金属在环境中不会降解、消失,其通过迁移、转化等过程富积到农作物或其它植物,这些富积的重金属可通过食物链等对人类健康造成威胁。区域土壤作物系统重金属特征是当前重金属污染的研究热点。
由于土壤重金属污染具有普遍性、隐蔽性、表聚性、不可逆性等特点，所以现对某城市城区土壤地质环境进行调查。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，不同的区域环境受人类活动影响的程度不同。为此，将所考察的城区划分为间距1公里左右的网格区域，按照每平方公里1个采样点对表层土（010 厘米深度）进行取样、编号。最终得到以下数据每个采样点的位置（见附件1），每个样本所含的多种化学元素的浓度数据（见附件2）。另外在该区的自然区取样，将得到的元素浓度作为该城区表层土壤中元素的背景值（见附件3）。
现要求，根据已给数据查阅相关资料，作出的8种主要重金属元素在该城区采样点的空间分布，并分析出该城区内不同区域重金属的污染程度，说明重金属污染的主要原因，分析重金属污染物的传播特征。建立模型，确定污染源的位置。最后分析所建立模型的优缺点，并指出为了更好地研究城市地质环境的演变模式，还应收集的信息。再根据这些收集的信息优化原有模型解决实际问题。
模型假设 1、 假设土壤的酸碱度不影响污染传播。
2、 假设区地势的总体分布情况不影响污染传播。
3、 假设人为因素引起的污染物传播具有均匀性。
4、 假设初始状态重金属元素分布具有规律性。
符号说明 单项污染指数；
各区域中每项元素浓度平均值；
该地区各项元素背景值；
单项重金属危害系数；
各金属毒性系数（毒性在综合污染中的权重）；
RI不同区域多种重金属的污染程度；
各金属污染源位置；
δ为狄拉克函数；
u污染源扩散介质的速度；
H扩散系数；
K污染物的降解率；
Lx的取值范围（0m-26854m）；
Mi污染物排放强度；
C污染物在x点的排放浓度；
模型的建立与求解 1、问题一的求解过程如下 1.1根据附表1中所给各采样点的坐标值，运用matlab软件画出每种元素在城区的空间分布图（运行程序见附件4）。图形如下，图片顺序依次为As（图1）、Cd（图2）、Cr（图3）、Cu（图4）、Hg（图5）、Ni（图6）、Pb（图7）、Zn（图8）。
图1（As在城区分布图） 图2（Cd在城区分布图） 图3（Cr在城区分布图） 图4（Cu在城区分布图） 图5（Hg在城区分布图） 图6（Ni在城区分布图） 图7（Pb在城区分布图） 图8（Zn在城区分布图） 1.2、由附表2所给各元素在各类区的浓度数值，以及附表3中各元素的背景值，由公式 用excel处理数据，分别计算出每区的各元素单项污染指数，结果（见表1） 生活区 1.74 2.23 2.23 3.74 2.66 1.49 2.23 3.43 工业区 2.01 3.02 1.72 9.66 18.35 1.61 3.00 4.03 山区 1.12 1.17 1.26 1.31 1.17 1.26 1.18 1.06 主干道路区 1.59 2.77 1.87 4.71 12.77 1.43 2.05 3.52 公园绿地区 1.74 2.16 1.41 2.29 3.29 1.24 1.96 2.24 （表1） 1.3、根据各区元素的单项污染指数，由公式 求出单项重金属危害系数。其中为各金属毒性系数（毒性在综合污染中的权重），查阅资料后（参考文献[1][3]），按8 种重金属毒性系数分别为Hg40＞Cd30＞As10＞PbCuNi5＞Cr2Zn1。用excel处理数据后（见表2）。
1.4、得出不同区域多种重金属的污染程度RI，结果见表2。由公式 生活区 工业区 山区 主干道路区 公园绿地区 E（As） 17.42 12 8 6 10 E（Cd） 2416 31663 16632 1523 28055 E（Cr） 38 1024 578 628 704 E（Cu） 68 1750 613 299 522 E（Hg） 1033 16601 7308 5083 14272 E（Ni） 25 101 85 76 75 E（Pb） 95 1786 945 645 1071 E（Zn） 66 18297 5658 4944 10405 RI 3763 71238 31832 26919 55120 （表2） 1.5、从表2中我们可以很容易的横向比较出各个区域8重金属的污染程度，最严重的是工业区，紧随其后的是公园绿地，随后依次是山区、主干道路区，污染程度最低的是生活区。
1.6、再查阅各相关资料（参考资料[4]），找出和RI污染范围，与表2中得出的数据比较，纵向比较单项污染系数和多重金属污染系数RI的值，得出结论。
单项污染系数和多重金属污染系数RI污染程度划分范围（如表3）。
指数类型 所处范围 污染程度 单项金属危害生态系数 公园绿地山区主干道路区 Cd 工业区公园绿地山区生活区主干道路区 Cr 工业区公园绿地主干道路区山区生活区 Cu 工业区山区公园绿地主干道路区生活区 Hg 工业区公园绿地山区主干道路区生活区 Ni 工业区公园绿地山区主干道路区生活区 Pb 工业区公园绿地山区主干道路区生活区 Zn 工业区公园绿地山区主干道路区生活区 （表5） 再根据附表2中的数据，计算城区中各元素的平均值、标准差、变异系数、超标率、背景值和平均值比背景值的倍数。数据见表6。
金属 Asu/g Cdu/g Cru/g Cuu/g Hgu/g Niu/g Pbu/g Znu/g 实际范围 1.61- 30.1 40- 1619.8 15.32-900.84 2.29- 2528.5 8.57- 16000 4.27- 142.5 19.68-472.48 32.86- 3760.8 平均值 5.68 302.4 53.51 55.02 299.7 17.26 61.74 197.07 标准差 3.02 225 70 163 1629．5 9.94 50.06 339.23 变异系数 0.84 1.73 2.26 12.34 46.56 0.81 1.61 4.92 超标率 0.77 0.80 0.81 0.88 0.66 0.75 0.82 0.79 背景值 3.6 130 31 13.2 3