基于Monte Carlo模拟潜在生态危害评价.docVIP

基于Monte Carlo模拟潜在生态危害评价 　　摘要：将Monte Carlo模拟引入到Hakanson潜在生态危害指数评价中，可以有效地解决由于评价区域重金属浓度随机性和不确定性造成的评价风险，将其应用于后湖农场土壤中重金属Mn、Cu、Pb、Zn、Cd的潜在生态危害评估。结果表明，评价区域总体处于中等潜在生态危害等级；Cd处于较高的风险等级，处于较高潜在生态危害及以上的风险概率达到了84.87%。通过敏感性分析，Cd含量的敏感性系数达到了99.1%，对于该区域的潜在生态危害起主导作用。 　　关键词：Monte Carlo模拟；Hakanson潜在生态危害指数；敏感性分析 　　中图分类号： X825文献标志码： A文章编号：1002-1302（2014）06-0320-03 　　收稿日期：2013-09-10 　　基金项目：国家科技支撑计划（编号：2011BAD13B02）。 　　作者简介：曲春风（1981―），男，山东牟平人，硕士，实验师，研究方向为环境生态学。Tel：（0631）5688551；Email：qcf_81@163.com。 　　通信作者：葛长字，博士，副教授。E-mail：changzige@。由于人类的工农业生产活动，使得大量的重金属进入到土壤中，其赋存形态因生物活动而发生改变但不能降解，致使土壤中的重金属得以积聚[1]。土壤中的重金属除污染水质外，还通过食物链对人类健康造成危害[2]。土壤中重金属的环境监测、生态危害评价等成为国内外环境、生态等领域科学工作者关注的热点[3-4]。随着“镉大米事件”的曝光[5]，我国土壤重金属污染问题更加受到重视。目前，普遍采用的土壤重金属污染的评价方法有单因子指数法、地质累计指数法、内梅罗污染指数法、污染负荷指数法和潜在生态危害指数法等[6]。Hakanson提出的潜在生态危害指数法[7]综合考虑了重金属的生态毒性及重金属区域背景值的差异，消除了区域差异影响，在国内外土壤、沉积物评价中得到广泛的应用[8-9]。该评价方法不仅受制于监测样品数量的有限性、区域污染物分布的时空不均匀性以及实测数据的误差性，而且对于重金属生态毒性的权重赋予具有一定的主观性。潜在生态危害指数法不能有效解决土壤中重金属污染评价的不确定性，即评价的风险性难以体现，而评价的风险性是决策的重要依据之一。将风险评价系统的不确定性引入到潜在生态危害指数中，将能更科学地评估土壤重金属污染的现状。本研究将Monte Carlo模拟引入到潜在生态危害指数法中，对土壤中的Mn、Cu、Pb、Zn、Cd等5种重金属的潜在生态危害及其风险性进行评估，并对各种重金属浓度的参数敏感性进行了分析。 　　1材料与方法 　　1.1基于Monte Carlo模拟的潜在生态危害评价方法 　　1.1.1Hakanson潜在生态危害指数法该方法由Hakanson于1980年提出，包含环境化学、生物毒理学、生态学等方面内容，常应用于沉积物及土壤的生态危害评价。基本公式为： 　　Cif=Ci/Cin 　　Cd=∑mi=1Cif 　　Eir=Tir?Cif 　　RI=∑mi=1Eir=∑mi=1Tir?Cif=∑mi=1Tir?Ci/Cin 　　式中：Cif为金属i污染系数，Ci为金属i的实测浓度，Cin为金属i的背景值；Cd为多种金属综合污染程度；Tir为金属i的生物毒性响应因子；Eir为金属i潜在生态风险因子；RI为多种金属潜在生态危害指数。 　　1.1.2Monte Carlo模拟Monte Carlo模拟因著名的摩洛哥赌城而得名，又称为随机模拟，基本思想是根据待求随机问题的变化规律，根据物理现象本身的统计规律，或者人为的构造一个合适的概率模型，依照该模型进行大量的统计试验，使它的某些统计参量正好是待求问题的解[10]。随着计算机技术的快速发展，使得快速进行Monte Carlo模拟成了可能，通过计算机模拟随机过程，进行大量模拟试验，并统计计算结果。 　　1.1.3基于Monte Carlo模拟的潜在生态危害指数评价利用Monte Carlo方法进行土壤中重金属潜在生态危害评估：（1）通过实测重金属浓度确定各种重金属在土壤中分布的概率密度函数；（2）通过Monte Carlo方法对概率密度函数进行随机抽样，获得重金属在土壤中的分布浓度，通过Hakanson潜在生态危害指数的方法，对土壤重金属潜在生态危害进行评估。随机模拟过程及评价结果的计算由Crystal Ball 11.1软件进行。对于不同金属的生物毒性响应因子Tir值，本研究参考经典的Hakanson指数以及国内研究成果[11-13]，设定5种重金属Tir的顺序为：Cd=30Cu=Pb=5Zn=Mn=1。 　　1.2评价区域 　　后湖农场位于