1、二、模型假设1、 假设题目中所给数据可靠无误。2、 假设问题一中各区平均的污染程度可以看做该区的污染程度。3、 假设问题二中只考虑题目中所给的8中重金属，不考虑其它重金属。4、 假设重金属传播特征不受风向等因素影响。三、问题分析（一） 问题1的分析：问题1属于空间分布和综合评价问题，重金属的传播过程是一个扩散的过程，通常物质扩散模型中物质从高浓度向低浓度扩散且其浓度的分布是连续的，据此我们可以用附表中所给的采样点污染数据为基础借助MATLAB软件1进行插值拟合得出8种主要重金属污染物在整个城区的空间分布图。对于该城区内不同区域重金属的污染程度的研究可以借助我国土壤监测技术规范（HJ/T 166
2、-2004）2中推荐的内梅罗综合污染指数法进行评价，求出不同区域重金属的污染等级。（二） 问题2的分析：问题2要求通过数据分析来说明重金属污染的主要原因。首先可以对重金属和海拔进行相关性分析，得出相关矩阵和相关度，再结合问题一求出的结论分析出重金属可能的主要来源和重金属污染的主要原因。（三） 问题3的分析：由问题一的分析我们得知重金属的分布是连续的，同时我们还可以知道物质的扩散是从高浓度向低浓度进行的，在扩散模型中某区域浓度最高的点可能就是扩散源，所以重金属空间分布中的极值点就可能是重金属的传播模型中污染源。因此问题三的求解就转化为在模型一所拟合出的重金属空间分布曲面上搜索极值的问题。搜索极值
3、的现代算法有模拟退火，遗传算法，鱼群算法等多种。考虑的模型中所搜索的域有限，且目标解数目不确定，遍历搜索是较好的方法。得出极值点后再结合国家土壤环境质量标准筛选出污染源。（四） 问题4的分析：首先应对问题一，二，三所建立的模型进行优缺点分析然后根据影响城市演化模型的因素，分析还应搜集的数据以及模型如何建立的问题。四、符号说明符号设定符号说明Pij区域i中第j个重金属的污染分指数Cj第j个重金属的实测浓度Sj第j元素的评价标准PN综合污染指数Pj,ave平均单项污染指数Pj,max最大单项污染指数z浓度分布矩阵注：在此没有设定的符号在下文中会具体说明。五、模型的建立及求解一、问题一的求解： 1.
4、1 用MATLAB软件对所给数据进行插值拟合得出调查区的地形图和8种主要重金属元素在该城区的空间分布图，再用MATLAB软件对所给数据进行分析得出功能区散点图：图1：调查地区的地形图图2：功能区散点图图3：砷和镉在该城区的空间分布图图4：铬和铜在该城区的空间分布图图5：汞和镍在该城区的空间分布图图6：铅和锌在该城区的空间分布图说明：图1的Z轴为海拔高度，X、Y轴为地理坐标值（单位：m）。图2 的X、Y轴为地理坐标值（单位：图3-图6的Z轴为重金属元素的浓度（单位：g/g），X、Y轴为地理坐标值（单位：1.2 模型建立：土壤环境质量单项污染指数主要用来评价某一污染物的污染程度，指数小污染轻，指数
5、大污染则重。但区域内土壤环境质量作为一个整体和外区域进行比较时除用单项污染指数外，还常用综合污染指数。综合污染指数可以综合判断某土壤多种污染物的联合污染效应。目前土壤环境质量评价方法有很多，各有优点和缺点。本文根据我国土壤监测技术规范（HJ/T 166-2004）2中推荐的内梅罗综合污染指数法进行评价。在计算某个区域某种重金属单项污染指数（分指数）的基础上，再计算该区域多种重金属的综合污染指数。单项污染指数和综合污染指数的计算公式如下： （1） （2）当Pij1时，表示土壤未受该因子污染，当Pij 1时，表示土壤受该因子污染。内梅罗综合污染指数反映了各污染物对土壤的作用，同时突出了高浓度污染物
6、对土壤环境质量的影响。根据HJ/T 166-2004,内梅罗综合污染指数的分级标准（见表1），得出各个区域的污染等级。表1：内梅罗综合污染指数的分级标准等级内梅罗污染指数1PN0.7清洁（安全）20.7PN1.0尚清洁（警戒线）31.0PN2.042.0PN3.05PN3.01.3 模型求解：本文以背景值作为评价标准进行求解，用EXCEL对文中所给数据进行分类，把数据分入1类区、2类区、3类区、4类区、5类区。然后得出各个区里面主要重金属含量的平均值，可看作各区中主要重金属含量值。如下表：表2：各区重金属含量的平均值As （g/g）Cd （ng/g）Cr （g/g）Cu （g/g）Hg （ng
7、/g）Ni （g/g）Pb （g/g）Zn （g/g）6.27289.9669.0249.493.0418.3469.11237.017.25393.1153.41127.54642.3619.81277.934.04152.3238.9617.3240.9615.4536.5673.295.71360.0158.0562.21446.8217.6263.53242.856.26280.5443.6430.19114.9915.2960.71154.24然后根据公式（1）、（2）结合MATLAB软件算得各区重金属单项污染指数和综合污染指数，如下表：表3：各区重金属单项污染指数和综合污染指数单
8、项 污 染 指 数AsCdCrCuHgNiPbZn1.74172.23052.22653.74242.65831.49112.22943.43493.17042.01393.02391.72299.662118.35311.61063.00134.02813.53311.12221.17171.25681.31211.17031.25611.17941.06221.25321.58612.76931.87264.712912.76631.43252.04943.51969.42641.73892.1581.40772.28713.28541.24311.95842.23542.7343再由内梅
9、罗综合污染指数的分级标准得出各区的综合污染等级，如下表：表4：各区综合污染等级5 3 轻等污染4 中等污染从表中可以看出，该城区内生活区、工业区、主干道路区属于重污染区，公园绿地区属于中等污染区，山区属于轻度污染区。二、问题二的求解：2.1 模型建立：用SPSS110统计软件对各种重金属元素浓度和海拔做相关性分析，得出各种元素与元素之间和元素与海拔之间的相关系数矩阵及其相关性，结合第一问得出的空间分布图和区域散点图，参照主要重金属含量土壤单项污染的指数，分析得出各重金属污染的主要原因。2.2 模型求解：2.2.1 根据题中所给数据，以As、Cd、Cr、Cu、Hg、Ni、Pb、Zn八种重金属元素
10、浓度和海拔作相关性分析，经SPSS110统计软件进行相关性分析，得出该市表层土壤As、Cd、Cr、Cu、Hg、Ni、Pb、Zn八种重金属原始含量数据和海拔的相关系数矩阵，如图所示图7: 重金属原始含量数据和海拔的相关系数矩阵可见各重金属浓度均和海拔成负相关，即海拔越高，其含各种重金属浓度越低；Cr和Ni的相关性最好，相关系数最大，为0.716，其次为Pb和Cd，相关系数为0660，以下是Cr和Cu的相关性较好，相关系数是0.532，其它元素之间的相关性并不是很好。从成因上来分析，相关性较好的元素可能在成因和来源上有一定的关联。结合第一问中8种主要重金属元素在该城区的空间分布可以看出，Cr和Ni
11、、Pb和Cd可能是来自同一来源。2.2.2 根据空间分布图、区域散点图和主要重金属含量土壤单项污染的指数进行分析：对于Cr和Ni，在来源上关联较密切，该市表层土壤Cr和Ni基本未污染，只有个别点富集程度较高，污染达到中度污染，该富集中心的位置主要分布在生活区周边和主干道区周边，这可能是由于生活废水的排放和交通源汽车尾气的排放等原因造成的。对于Pb和Cd，在来源上关联较密切， Pb和Cd的高含量点主要分布在交通繁忙的主干道路区周边和工业区周边,这可能是因为Pb和Cd来自该市中心交通源汽车尾气的排放、汽车轮胎的磨损和冶炼厂的废水、尘埃和废渣,以及电镀、电池、颜料、塑料稳定剂、涂料工业的废水等。所以
12、可以说Pb和Cd的污染主要是由于主干道污染和工业污染。对于Cu,该市表层土壤Cu基本未污染，只有个别点富集程度较高，污染达到中度污染，该富集中心的位置主要分布在生活区周边，这可能是由城市商业活动、城市居民生活累加到土壤中的Cu。对于Hg, 其高含量点主要分布在交通繁忙的主干道路区周边和工业区周边, Hg污染的一个主要原因是由于燃煤造成的，无论是工业用煤还是居民用煤，而且燃烧方式落后。工业排放也是表层土壤Hg污染的另一个重要来源，主要在大面积污染的几个工业浓集中心。对于Zn,其高含量点也主要分布在交通繁忙的主干道路区周边和工业区周边, 这主要是由于汽车尾气的排放和厂矿企业的三废排放。对于As，该
13、市表层土壤As基本都是轻度或中度污染，只有个别点富集程度较高，该富集中心的位置主要分布在工业区周边，主要来源可能是工厂的废水排放。综上所述，可以认为工业区、主干道路区和生活区的活动是造成该城区表层土壤重金属污染的主要原因。三、问题三的求解：3.1 模型建立：依据问题一得出的各重金属元素在该城区的空间分布，得到浓度分布矩阵Z（Z是100100的矩阵），进而结合MATLAB软件建立搜索模型。Z是100100的矩阵，借鉴元胞的思想建立一个100100规模的二维网格，将元素浓度分布矩阵对应放入，其中每一个元素占据其中一个格子。根据问题分析可知：污染源存在于二维网格中的某些格子中。并且污染源所在格子元素
14、浓度大于周围格子的元素浓度。二维元胞自动机（规则四方网格划分）的邻居通常有以几种形式如图2所示：黑色元胞为中心元胞，灰色元胞为该元胞的邻居。（参考文献4）图8：元胞邻居模型分析三种邻居模型发现第二种模型最适合。第二种邻居模型中污染源存在的格子z（i,j）应满足：对于与外界不相邻的格子 （3） （4） （5） （6） （7） （8） （9） （10）对于边界处的格子理论上应满足以左边界为例 （11） （12） （13） （14） （15）对于顶角处的格子理论上应满足： （16） （17） （18）为了简化模型在此不予考虑，即认为对于边界和顶角处不存在污染源。通过搜索模型可以求出重金属空间分布中的
15、极值点即可能的污染源，再结合国家土壤环境质量标准3（下表5）通过MATLAB软件对极值点进行筛选出，求出重金属的主要污染源。表5：国家土壤环境质量标准级别一级15200903515040100二级2530050050250三级30100040015003.2 模型求解：3.2.1 根据问题一中得出的砷元素在城区的空间分布（见图3）。得到浓度分布矩阵Z（矩阵较大未附出），结合MATLAB软件建立搜索模型进行搜索得出砷元素在空间分布极大值61个。用同样方法得出其他7种重金属在空间分布极大值个数（见表6）。表6：八种金属元素空间分布极大值个数元素As Cd Cr Cu Hg Ni Pb Zn 个数6
16、16057626353583.2.2 运用scatter函数画出各重金属元素空间分布极大值点的散点图使数据可视化。得到各种重金属元素空间分布极大值点的散点图如下：图9：砷、镉、铬、铜的空间分布极大值散点图图10：汞、镍、铅、锌的空间分布极大值散点图3.2.3 结合国家土壤环境三个等级的质量标准通过MATLAB软件对极值点进行分级筛选：首先用国家土壤环境一级质量标准进行筛选，得出筛选结果，再用国家土壤环境二级质量标准对一级指标得出的点进行筛选，依次类推，最终得到筛选结果如表7表7：不同国标等级下的极大值个数一级个数6134539二级个数522017三级个数1118从筛选的结果中选出适当的点作为重
17、金属的主要污染源，所选点个数和点坐标如下列各表：表8：重金属主要污染源个数表9：砷污染源二维坐标及其浓度值15.06123.64116.12123.17530.03218.971X/m1890012900720045001830027600Y/m22003200740078001020012200表10：镉污染源二维坐标及其浓度值1068.81458.61401.91321.91121.41054.924001770051002600340036004000420052001264.410241267.81263.81578.660004800216008600112001140011600表
18、11：铬污染源二维坐标及其浓度值747.81304.81976.76108005600表12：铜污染源二维坐标及其浓度值2759.42609.82622.32565.21391.927003800表13：汞污染源二维坐标及其浓度值1638514487154601542718392333134341341111432169217233000138003300153001560022500870082009200940010600表14：镍污染源二维坐标及其浓度值146.0870.58769.35522200表15：铅污染源二维坐标及其浓度值527.92485.07354.062100表16：锌污
19、染源二维坐标及其浓度值1485.6550.91631.514571749.23092.12801.41965.31961.981009600460058001111.91064.753205859526553793664.52985.5552.15400840098001000011000运用scatter函数画出各重金属元素主要污染源的散点图使数据可视化。各种重金属元素主要污染源的散点图如下：图11:砷、镉主要污染源的散点分布图图12:铬、铜主要污染源的散点分布图图13:汞、镍主要污染源的散点分布图图14:铅、锌主要污染源的散点分布图四、 问题四的求解：优点：解决问题一的第一小问时，我们用M
20、ATLAB对原始数据进行差值拟合。由于所给数据采样点的不规则性，首先使用griddata函数对所给数据进行插值规整得出一个X,Y分别等步长的某种元素的浓度分布矩阵。在规整的浓度分布矩阵基础上分别使用pcolor，contourf，contour，surf等函数绘出了各种重金属元素在城区的空间分布。通过综合比较之后选定三维surf曲面建立重金属元素的空间分布模型，直观明了。解决问题二时，我们用SPSS对各种重金属元素浓度和海拔做因子分析，得出各种元素浓度和海拔相关性，各元素浓度和海拔呈现负相关，正好验证了第一问中求得的山区各重金属浓度最低，污染程度最轻这一结果。模型三中依据模型一中建立的浓度分布矩阵建立了遍历搜索模型。该模型能够有效且快
