- 相关推荐
重金属传播模型
基于空间统计学模拟方法的城市重金属污染分析
丁宁，孙英君，蔡菲，孟 飞
(山东建筑大学土木学院，济南 250014)
摘 要：本文利用空间统计学的模拟方法，对土壤重金属污染状况进行研究，通过获取多幅模拟结果之间的差异来揭示研究区域土壤重金属污染的整体空间分布形态。文章以上海某区域Cu, Pb,两种重金属为研究对象，利用空间统计学模拟方法对其空间分布予以研究，给出每种重金属空间分布模拟结果。最后，文章以研究区域土壤环境背景上限值为标准，给出研究区域不同土壤重金属空间分布相应级别的不确定性分析结果。
关键词：模拟；地统计学；空间变异；重金属；不确定性
城市的土壤环境质量与人类健康息息相关,随着城市化进程的迅速发展,城市土壤重金属污染日趋引起人们的关注。研究者将汞、镉、铅、铬、砷、锌、铜、镍合称为重金属环境污染元素,尤其是Pb、Cu 及其复合污染最为突出, 是城市和工业发展过程中产生的典型的“污染元素”, 给区域土壤资源的保护、利用与管理带来了严峻挑战。据报道, 中国受重金属污染的耕地面积近2000 万hm , 约占总耕地面积的1/5。土壤重金属污染导致严重的经济损失, 如中国每年就因重金属污染而减产粮食1000多万t, 另外被重金属污染的粮食每年也多达1200万t, 合计经济损失至少200亿元。同时，土壤重金属污染也会导致重金属在植物体中积累, 并通过食物链富集到人体和动物体中, 危害人畜健康, 引发癌症和其他疾病等。其次, 土壤受到重金属污染会导致大气污染、地表水污染, 地下水污染等生态环境问题。
因此开展土壤重金属空间结构特征的研究，特别是对土壤重金属的全局空间分布形态进行模拟分析，从而对土壤重金属污染程度进行定量化的评价，对土壤重金属污染热点予以探测，从而对重金属污染源进行分析研究，及时有效地采取防护、修复措施，对减少城市重金属污染，提高人们的生活质量具有重要的意义。 [2]2[1]
一、土壤重金属空间结构及分布特征研究现状
已有一些研究者对城市重金属污染状况进行分析：Hu等利用地统计学对北京市大兴区土壤重金属含量的空间分布特征进行了研究; 郑袁明, 陈同斌等利用地统计技术对北京市近郊区土壤Ni 的空间结构及分布特征进行了研究; 钟晓兰等以江苏省太仓市为例，运用地统计技术研究了长江三角洲地区土壤重金属污染的空间变异特征;郑海龙等利用地统计学普通克里格插值方法对上海梅山刚体集团附近的土壤重金属空间变异进行研究；李亮亮等通过克里格插
Pilar Burgos值对葫芦岛市连山区、龙港区土壤重金属汞、铬、镍的空间分布进行单因子评价。
等运用Kriging方法对Seville, Spain地区的土壤修复前后重金属、微量元素的变化进行了分析
[8][7][6][5][4][3]; M. Biasioli等利用地统计方法及GIS技术对大型城市发展对土壤重金属的影响加以分析。[9]Li等利用克里格方法，结合GIS技术对香港九龙六种重金属元素（Cd, Cr, Cu, Ni, Pb and Zn）作者一简介：丁宁，男，1961.11，汉，山东济南人，学士，副教授，主要研究方向为测绘科学技术在城市建设中的应用。
[10]的空间分布加以研究，并通过热点探测的方法进行污染源的求解。
总结起来，研究者大多基于GIS技术与空间（地）统计学克里格系列插值方法对城市重金属污染状况进行分析：利用GIS技术进行研究区域及其采样数据、空间分布插值结果的可视化表达；利用空间统计学的变异函数对采样数据的空间异质性予以分析，并通过理论变异函数的不同因子、系数来寻求异质性产生的原因（影响因子），最后给出克里格算法插值得到的某重金属空间分布结果。
二、存在的问题及解决方法
经过克里格插值得到的重金属空间分布结果不一定能够反映研究区域的真实状况，特别是当研究者的目标是发现研究区域的整体分布态势或进行热点探测时。因为克里格算法本身具有一定的平滑作用：即将极大值估计偏小，极小值估计偏大，因此插值后的结果趋于平缓[11]。但在某些问题的解决中，人们需要的恰恰是被克里格方法平滑掉的极值，例如在污染分析中，研究者主要是通过污染金属所表现的极值来发现污染源。此外，克里格方法只能为未采样点提供唯一的“最优”估计值，且仅仅能够提供该点估值的局部不确定性。
土壤重金属污染的全局模拟预测及其空间不确定性分析可通过空间统计学的系列模拟方法加以解决。在此，本文选择序列高斯模拟方法进行研究区域土壤重金属空间分布的模拟预测。该方法通过随机扰动因子，获得不同搜索路径，从而在某一具体点位预测时采用不同的已知采样点及已模拟点估值结果，最后得到相同概率的多个不同模拟结果。这些模拟结果与已知采样点的空间结构（变异函数或直方图分布）一致，但对某个单一的未采样点而言，不同轮次的模拟结果会有所差异，也正是这种差异，体现了模拟结果的空间不确定性问题。
三、应用实例
本文以上海某试验区域为研究对象，采用梅花多点取样方法共获取240个样品。土壤样品经风干、过100目筛后，用HNO3-HClO4-HF消化法制样，土壤中全量Cu、Pb采用X荧光法。 重金属的分析测定中，均加入国家标准土壤标样（GSS-1）进行分析质量控制，分析结果符合质量控制要求。
3.1 土壤重金属空间结构分析
研究区域土壤重金属空间结构分析如下表1所示。
从表1可以看出，研究区域内Cu、Pb两种重金属元素的变程分别为0.81、1.73 km，说明它们受到污水灌溉、施肥等随机因素的影响较大，导致了它们在一个较小的距离范围内存在相关关系。
该比值表1中块金值C0与基台值C0+C1之比是反映区域化变量空间异质性程度的重要指标，
用以反映空间变异影响因素中区域因素（自然因素）和非区域因素（人为因素）的作用作者二简介：孙英君，女，1976.8，汉，山东聊城人，博士，副教授，主要研究方向为空间信息分析技术在城市建设中的应用。 [12]。从
表1可以看出Cu、Pb两种重金属元素的块金常数/基台值均大于0.25小于0.75，属于中等空间相关。
表 1 研究区域土壤重金属的半方差函数模型
Heavy metal Model Nugget(C0) Sill(C0+C1) Range/km C0/(C0+C1)
Cu Spherical 152.450 361.660 0.81 0.422
Pb Exponential 39.829 65.294 1.73 0.609
3.2 土壤重金属空间分布模拟
两种土壤重金属空间分布模拟均在GS+ (v.9)软件环境下完成，下图1所示均是取自1000次模拟结果的平均值。为识别土壤重金属空间分布的整体态势，作者对每一类重金属选择10个以上取值区间成图。同时，本文选取研究区所在城市1987年土壤重金属背景值上限值（按95%置信限）为判断标准，分别是：Cu 40.6mg/kg、Pb 36.1 mg/kg。
对比各重金属标准值，可以看出，研究区域西北区域铜要素严重超标，Pb元素也存在较明显超标斑块。而从研究区域而言，西北、中南偏下部位是受各种重金属污染较为严重的区域。
图1. 研究区重金属空间分布模拟结果
CuPb13]
3.3土壤重金属空间分布模拟不确定性分析
如3.2所述，文章选取研究区域所在城市1987年土壤重金属背景值上限值为阈值，对取得的土壤重金属空间分布模拟结果进行不确定性分析，分别给出每一种金属超出阈值的概率分布图。
由图2可以看出，在Cu要素的概率图上，出现了概率较高的紫色区域，只有少部分区域超标的概率小于50%，预示了铜要素在该研究区域的严重污染。这些区域比较集中在图幅的西北角及其中南偏下部位，与3.2中的分析取得一致结果。在Pb要素的概率图上，存在特定斑点区域超标，几率在60%上（红色系以上色彩表示），但大部分研究区域污染程度较轻。
Prob. of Cu Prob. of Pb
图2. 土壤重金属空间分布不确定性预
四、总结
以上研究，由于资料所限，作者并没有将分析结果与研究区域的自然、社会状况相关联，如污染区域是否与研究区域的厂矿企业分布一致等。但从上文分析可以看出，经过1000次高斯序列模拟，其结果已相对稳定，能够揭示研究区域整体空间分布形态。特别是结合研究区域土壤重金属背景值对模拟结果进行不确定性分析，可以对模拟结果做有效的印证，同时有时还会发现一些在模拟结果中没突现出来的特征污染斑点。将二者结合起来，可以对重金属污染区域进行评价区划，为决策部门对土壤的整治与修复提供科学的依据。
参考文献
[1] 刘勋鑫，王翠红，廖超林等. 土壤重金属特性空间变异性方法研究[J]. 农业网络信息，2007，6：106-108.
[2] 隋红建，吴璇，崔岩山. 土壤重金属迁移模拟研究的现状与展望[J]. 农业工程学报，2006，22(6):197-200.
[3] Hu Ke-lin, Zhang Feng-rong, Li Hong etc.. Spatial Patterns of Soil Heavy Metals in Urban-Rural Transition Zone of Beijing[J]. Pedosphere, 2006, 6:690-698.
[4] 郑袁明，陈同斌，陈煌等. 北京市近郊区土壤镍的空间结构及分布特征[J]. 地理学报，2003，3：470-476.
[5] 钟晓兰，周生路，李江涛等. 长江三角洲地区土壤重金属污染的空间变异特征——以江苏省太仓市为例[J]. 土壤学报，2007，1：33-40.
[6] 郑海龙，陈杰，邓文靖等. 城市边缘带土壤重金属空间变异及其污染评价[J]. 土壤学报，2006，1：39-45.
[7] 李亮亮，张大庚，王延松等. 葫芦岛市连山区、龙港区土壤重金属汞、铬及镍的空间分布及其污染评价[J]. 土壤通报，2006，6：1207-1211.
[8] Pilar Burgos, Engracia Madejo?n, Alfredo Pe?rez-de-Mora etc.. Spatial variability of the chemical characteristics of a trace-element-contaminated soil before and after remedi- ation[J]. Geoderma, 2006, 130:157-175.
[9] M. Biasiolia, R. Barberis, F. Ajmone-Marsan. The influence of a large city on some soil properties and metals content [J]. Science of the Total Environment, 2006, 356: 154–164.
[10] Xiangdong Li, Siu-lan Lee, Sze-chung Wong etc.. The study of metal contamination in urban soils of Hong Kong using a GIS-based approach [J]. Environmental Pollution, 2004, 129:113–124.
[11] Pierre Goovaerts. Geostatistics in Soil Science: State-of-the-Art and Perspectives [J]. Geoderma, 2002, 89:1-45.
[12] 徐建华. 现代地理学中的数学方法[M]. 北京: 高等教育出版社, 2002: 105-115.
[13] 王云，汪雅谷，罗海林等. 上海市土壤环境背景值[M]. 北京:中国环境科学出社, 1992: 37-57.
【重金属传播模型】相关文章：
论传播的分类及其数学模型04-27
基于群体稀疏分布的SIS疾病传播模型04-26
基于信息传播系统的人员疏散模型探讨04-26
HIV传播的人群生态动力学模型05-02
基于二维规则网格的SIRS病毒传播模型04-29
全球气候变化信息传播扩散模型及实证分析04-26
模糊数学模型在土壤重金属污染评价中的应用04-30