重金属污染扩散模型

重金属污染扩散模型
重金属污染扩散模型

扩散模型

摘要:

本文回顾了颗粒物大气扩散模型的应用,概括介绍了现有的不同扩散模型,从简单的箱式模型到复杂的流体动力学模型,以及扩散模型在不同环境中不同方法的适用性,考虑应用尺度、环境复杂性以及浓度参数化的确定。最后,回顾了几个商业的和非商业的粒子扩散软件(packages),详细的介绍了它们的使用过程、在应用方面的优势和局限性。回顾的模型包括:箱式模型(AURORA, CPB和PBM)、高斯模型(CALINE4, HIWAY2, CAR-FMI, OSPM, CALPUFF, AEROPOL, AERMOD, UK-ADMS 和 SCREEN3)、拉格朗日/欧拉模型(GRAL, TAPM, ARIA Regional)、CFD 模型(GRAL, TAPM, ARIA Regional)和气溶胶动力学模型(GATOR, MONO32, UHMA, CIT, AERO, RPM,AEROFOR2, URM-1ATM, MADRID, CALGRID和UNI-AERO)。

1、引言

扩散模型描述大气、扩散和烟羽的化学物理过程,运用数学方程式来计算不同位置的浓度。同时,还有很多关于大气模型、街道峡谷中扩散的研究方法(Vardoulakis et al.,2003),不同模型运用试验气象数据比较(Ellis et al., 2001;Sivacoumar and Thanasekaran, 2001; Hall et al.,2002; Caputo et al., 2003) 的评论文章,这些都是集中在模拟气体的扩散上。

不幸的是,只有很少的研究同时测定了粒子和气体浓度,这两个研究的差别可能部分源于观测的不同。在开放的生境中,几项研究已经揭示气体和粒子浓度变化的相关性。尽管在后来的研究中只有两个采样点,但是Monn et al. (1997)指出了城市环境下PM10和NO2浓度的相关性较差,以及PM2.5和NO2较好的相关性。相反,Clairborn et al. (1995)揭示了SF6 和PM10较好的相关性,尽管只是观测了距离高速公路60m的范围。Roorda-Knape et al.(1998)指出苯、PM2.5和PM10的浓度在公路主干道300m范围内没有明显的降低。这和Hitchins et al. (2000)提出的PM2.5浓度有很小的降低相一致。在那项研究中作者指出高速公路周围颗粒物数量浓度下降比NO2快。Zhu et al.(2002a, b)指出高速路周围颗粒物数量浓度在6-220nm之间,这和CO浓度密切相关。所有这些研究都是在风向垂直的远离道路开放环境中进行的。然而Morawska, (2003; Holmes et al., 2005)已经观测到了局部气体和颗粒物扩散的区别。同时测定的CO和粒子数量浓度显示观测点附近的CO和粒子数量浓度没有很明显的相关性,在观测点周围两个污染物的比较研究指出了不同空间和时间的趋势。Harrison and Jones(2005)在城市观测点的另一个研究显示粒子浓度和NOx相关性较差,但是在人行道监测点两者有很好的相关性,这是由于浓度受扩散影响很小。此外,对(Morawska,2003)很多野外研究验证发现在建筑物附近的粒子数量浓度的垂直廓线不同于气体。这些不同于先前在较复杂环境中进行的研究,在复杂环境中,湍流严重影响风向,污染物的排放也不仅限于简单的线源。一般来说,在开放环境中观测到的气体和粒子浓度相关性很好,可是在较复杂的城市环境中两者有很大的不同。Van Dingenen等在 2004年指出在颗粒物源以交通排放为主的城市环境中所有监测网络中的网格点的PM2.5和PM10的相关系数R2=0.95。可是,

PM10/PM2.5的斜率变化太多而不足以给出一个简单的PM10/PM2.5斜率。在同样

的研究中他们发现年平均粒子浓度和PM10或PM2.5的浓度都没有相关性。这和Harrison et al. (1999)在城市环境中的研究发现相反,他们认为小时粒子数量浓度与PM2.5比与PM10有更好的相关性,尽管3月份两个PM水平和小时粒子浓度都有很好的相关性。

所以,设计成可以模拟像惰性气体的稳定污染物的模型可以模拟一定开放环境中的PM10和PM2.5的浓度,特别是对于较长的平均期以及较大的范围,在较大的范围内由于瞬时粒子形成的短期变化是平缓的。

粒子数量浓度模拟包括气溶胶动力模块嵌接到扩散模型。所以,讨论粒子扩散模拟必须考虑处理粒子扩散各种方法的局限性以及用来评价发生在烟羽中的粒子过程的气溶胶动力模块。更复杂的情况的进一步研究,Lohmeyer (2001)注意到不同模型因为四个因子中的一个不同计算的浓度也不同,即使同一个模型不同组之间的结果也不同。预测浓度和观测浓度的一致性依赖于输入数据的质量。

这些回顾描述了不同模型类型,认识到不同空间尺度从局部到区域模型的特殊需要,以及粒子扩散和气溶胶动力学不同尺度关系的不足。另外,没有全面回顾所有现有模型只是包含了其中的大多数模型和比较重要的模型参数,表1、表2列出了模型输入。

虽然存在几个模型可以模拟粒子扩散,但是对扩散动力学结果没有特殊处理的只能计算一般以PM2.5和PM10的形式存在的粒子质量,不能计算其数量浓度。

此外,粒子有效性研究对很多模型而言是不可用的。考虑到这种情况作者已经试着在气体扩散的有效性研究中突出显示模型的性能。后来很多研究指出非活性气体和大范围中的粒子之间有很好的一致性,有效性研究包括气体在计算粒子质量浓度方面应该是一个好的模型性能指示器,正如前面讨论的。此外,平均气体和粒子浓度之间不同平均次数做比较是困难的,也就是说,不能确定模型预测的气体和粒子浓度之间变化的相关性。

一些包括气溶胶动力学深度处理的局部和区域模型出现了。这些模型中的大部分是非商业化的软件,并且其中的大部分已经和现有的扩散模型耦合,可以提供模拟粒子数量浓度在不同粒径尺寸中变化的软件。这意味着这些模型和扩散模型的性能依赖于包含在动力学模块中的准确和精确的过程。一般可以把气溶胶动力学模块嵌接到不同的扩散模型以适应耦合的扩散软件,使其更好的适合于将来的研究。

2. 模拟方法

2.1 箱式模型

箱式模型的基础是质量守恒。把观测点看作一个箱子,污染物向这里面排放且经过化学物理过程。它需要输入简单的气象状态以及污染物在箱子里面和外面的排放、运动。箱子内部没有明确定义,如果气团混合充分并且浓度均匀就作为箱子处理。箱式模型的一个优点是简化的气象学箱式模型可以包括更多详细的化

学反应模式(e.g. Master Chemical Mechanism)和详细的气溶胶动力学处理,这样可以更好的表达大气中粒子的化学物理性质。然而,输入初始条件后在没有提供当地污染物浓度任何信息的情况下箱式模型可以模拟箱子内污染物的形成。因为这个原因它们不适合模拟局部环境的粒子浓度,在局部环境下浓度和粒子动力学受风场和排放的局部变化影响很大。

2.2 高斯模型

高斯模型广泛的应用于大气扩散模型中,特别是为了常规目的,经常嵌套在拉格朗日和欧拉模型中。高斯模型是基于烟羽在稳定条件下垂直和水平方向的扩散。烟羽的正态分布在较大的尺度上修正,这是因为当混合高度较低时地球表面和边界层湍流反射的影响。烟羽厚度由和决定,和由稳定等级(Pasquill, 1961; Gifford, 1976)或者距离源的传输时间确定。考虑了模拟粒子扩散的烟羽模型很明显的局限性是使用稳定状态近似法,没有考虑污染物传输到受体需要的时间。因此,气溶胶动力学必须通过结果的后处理来计算。而且区域模型一般需要嵌接化学模型,通过二次有机气溶胶(SOA)的形成来准确的预测粒子形成。NOx和SOx化学对确定粒子和臭氧浓度是重要的,一般只用简单的指数衰减计算。更先进的模型用化学后处理可以模拟一些化学输送。尽管大部分的高斯模型只考虑污染物的扩散和平流,近来已经开发了更先进的高斯模型,包括物理过程(例如沉淀)和快速化学反应。而且,高斯烟羽方程假设烟羽之间没有反应,这在城市环境中是很重要的。

已经开发了软件来模拟烟羽内的化学、物理过程和建筑物周围的扩散。建筑物尾流的影响可以通过修正扩散系数和获得。但是,高斯方程不能计算由多重建筑物或在他们的交叉口引起的再循环的影响。

高斯烟羽模型中本身的一些局限可以通过近似当做一系列烟团的排放克服,烟团允许风速变化。这种方法中,每个烟团根据高斯扩散方程扩散,源的全部贡献通过每个烟团关于时间的积分和每个烟团在受体位置的贡献叠加来计算。

为了计算污染物在城市地区上方的浓度,经常使用复合源烟羽。源的性质以及源和受体的高度决定使用不同的方程。

高斯处理另外的一些局限意味着高斯模型不是设计用来模拟低风速条件下

或者监测点靠近源(例如距离小于100m)的扩散。(Benson, 1984; Sokhi et al., 1998)已经证明高斯模型过高预测了低风速下的浓度。混合模型,高斯模型和烟团模型的结合,包括顺风向污染物的扩散以更好的计算低风速条件下的浓度(Sharan et al., 1996; Thomson and Manning, 2001)。另一个局限性是湍流和气象简化处理的结果,所以他们更适合于计算小时污染物浓度。

高斯烟羽方程假设一个均匀风场,如果期待气象学在大距离上变化不推荐其用于远场模拟。Caputo et al. (2003)指出四个高斯模型计算整个下风向区域的非零浓度,所以建议他们应该限制在距离源几千米范围内。

2.3 拉格朗日模型

拉格朗日模型类似于箱式模型,它们定义了一个空气区域作为箱子包括污染物起始浓度。当污染物顺风向输送时拉格朗日模型追踪箱式模型的轨迹。当污染物从x运动到x’时浓度是源项和概率密度函数的产物。

拉格朗日模型合并由于平均流向速度、风组分的湍流和分子扩散导致的浓度变化。

拉格朗日模型适用于平坦地形的均匀和稳定条件 (Oettl et al., 2001; Raza et al., 2001; enkatesan et al.,2002; Tsuang, 2003) 以及复杂地形的非均匀和非稳定条件 (Du,2001; Hurley et al., 2003; Jung et al., 2003)。它可以模拟非线性化学,通过计算每个网格的浓度得到区域上网格浓度的叠加或者把粒子当作扩展的箱子,且模型中光化学模块适用于每个盒子。

气象数据计算风速变化的方差和拉格朗日自相关系数。拉格朗日粒子模型可以通过半随机系数的产生计算扩散系数,这些半随机系数不受稳定等级或者sigma曲线限制,这正是高斯扩散模型的情况。

2.4 计算流体动力学模型

计算流体动力学模型(CFD)以质量和动量守恒为基础,通过在三维模型中用有限差分和有限体积法解Navier–Stokes方程提供了流体运动的复杂分析。湍流的典型计算是用闭合度方法计算出现在动量和污染传输方程中的均质涡流粘度参数,方程假设污染物平均分布在各个方向。这种处理在平坦的边界层变现很好。但是,当存在分层边界层时闭合度方法需要校正,校正包括科里奥利力和上层大气缩减的风切应力,它们会导致涡流粘度估计过高。

Gidhagen et al. (2004)指出不同的CFD模型和所有风向场有很好的一致性,但是论证指出虽然模型输入是一致的,但是在速度和湍流级别上有很大不同。和风洞数据的比较说明这是不同模型封闭机制的结果。

3. 街道环境中扩散模型的综述

Vardoulakis et al. (2003)给出了城市扩散模型的回顾,这里给出模型的简单总结和模拟粒子扩散时适用性的讨论。尽管有一部分的扩散模型用来计算局部环境的城市污染物浓度,但是它们中的很多也包括街道峡谷环境中风流的复杂处理,只有三个模型包括计算粒子扩散的一个模块。

3.1 不包括气溶胶动力学特殊处理的模型

不包括化学形成(成核)或者气溶胶动力学(凝聚、凝结等)模块的扩散模型不能计算粒子数量浓度。模型一般基于质量守恒所以他们可以模拟粒子扩散,在PM2.5和PM10水平。模型一般以和气体相似的方式处理粒子,正如前面的讨论,依赖于平均周期和位置是危险的。但是,既然现在的空气质量规章基于粒子质量浓度,所以简单的粒子扩散模型是必要的,大部分普遍使用的调整模型的性能将会在下面回顾。

3.1.1 箱式模型

AURORA(VITO, Belgium)是一个综合的空气质量模型,它已经用来模拟城市环境惰性和活性气体、粒子的浓度(Mensink et al.,2003)。模型运用稳定状态下的箱式模型来计算街道峡谷中的污染物浓度。模型假设街道浓度一致,但是包括来自于峡谷逆风层顶流动中的湍流间歇现象。在箱子中考虑x和z方向的对流。

CPB(GEOMET)是城市峡谷模型,适用于高度和宽度比率在0.5和2之间的城市峡谷。模型利用街道峡谷内三个不同风向的NO和臭氧反应的简单算法式计算惰性气体和NO2的平均浓度。

3.1.2 高斯模型

用来计算汽车排放物扩散的最普通模型中的两个是CALINE4 (California Department of Transportation)and HIWAY2 (US EPA)。这两个模型以高斯烟羽模型为基础,并且受到短距离和复杂环境的城市扩散模型高斯方程的固有限制。另外,不推荐他们在低风速下模拟。尽管存在这些问题,但是它们已经应用在很多研究中,这是因为他们的使用方便并且不需要太大的计算能力和时间。

这两个模型把交通作为无限线源,并把它们分成垂直于风向的一系列元素。CALINE的垂直扩散参数考虑由车辆引起的热力和机械湍流,但是,HIWAY2只考虑车辆的影响而忽略垂直方向热力湍流的影响。除了上面提到的问题,高斯模型(e.g.CALINE4 AND HIWAY2)缺少在峡谷中模拟需要的精确度,因为建筑物只能通过地表粗糙度的变化来表示。

CAR-FMI(Finnish Meteorological Institute)是基于Luhar and Patil (1989)方程的高斯烟羽模型。它是设计用来计算汽车排放的CO, NO, NO2, NOx 和 PM2.5的小时浓度。大气稳定度用边界层尺度定义。和其他高斯模型相比,CAR-FMI的局限是不能在低风速条件下使用。水平和垂直扩散参数(和)包括周围风速的湍流项,排放速度和车辆。

和CAINE4一样,CAR-FMI使用分段气团法模拟惰性气体和活性气体以及颗粒物的扩散。然而,和CAINE4不同的是它包括三个粒径分组的干沉降处理。Oettl et al. (2001)论证了在芬兰主干道上监测的NOx小时浓度和由CAR-FMI与GRAL 模型预测的浓度一致性相当好。但是,CAR-FMI不能预测低风速条件下的弯曲风流。

OSPM是一个半经验模型,它用高斯烟羽方程求导得出源的直接贡献分布,用箱式模型来计算湍流对浓度的影响(Vignati et al., 1999)。烟羽的逆风扩散是不容忽视的,污染源一般视为无限线源。对线源来说烟羽表达式是沿着由街道风定义的路径嵌接的。

街道风的风向假设是镜反射,街道风可以代表顶层的风。街道水平的风速根据大气风速、风向和表面粗糙度计算。街道峡谷内各种各样湍流过程处理和贡献的变化依赖于大气风速、风向,读者可以关注Vignati etal. (1999) 的完整描

述。模型允许湍流对逆风和峡谷背风面浓度的影响,忽视街道迎风面直接排放的组分。

模型假设交通排放是沿着峡谷均匀分布的,根据经验导出车辆导致的湍流影响。对湍流的简单处理意味着模型不能模拟(wind flow风向)风向的间歇起伏,所以不推荐用于计算时间尺度短于一小时的浓度。此外,模型没有考虑扩散后排出烟羽的冷却,冷却可能对SOA粒子的形成有重要的影响。

OSPM模型用来评价城市街道峡谷中NOx, NO2, O3 和CO的监测数据(Kukkonenet al., 2003)。顶层和街道水平预测的小时平均浓度有很好的一致性。同时,相关性指出可以用原有的背景浓度和预处理气象数据预测小时浓度,没有尝试预测短时间周期的浓度。

在(Ketzel et al., 2000)第二阶段的研究中,OSPM和街道预测值的一致性不是很好,不能准确的预测在两个街道峡谷中不同风向对小时NO2浓度的影响。然而,它可以准确的计算交通车辆对街道年均苯浓度的贡献。这说明了峡谷内的风向简化处理不能反映浓度的短时间变化。

3.1.3 拉格朗日(欧拉)模型

GRAL(Institute for Internal Combustion Enginesand Thermodynamics, Graz, Austria)是一个耦合的欧拉-拉格朗日模型,设计用来模拟不均质风场内惰性组分的扩散。模型的一个局限性是没有考虑粒子的化学形成(e.g. ammonium nitrates,ammonium sulfates)。模型可以计算线源、点源从10min到1h的浓度以及平坦 (Oettl et al., 2005)和复杂的地形(Oettl et al.,2003) 隧道口的浓度。大气稳定度用边界层尺度计算。模型假设在隧道口附近稳定的烟羽抬升是周围空气和隧道流之间温度差的函数。

现在不存在使用GARL的粒子有效性研究,尽管设计者现在已经普遍涉及城市PM10和PM2.5精确的模拟。但是,模型准确模拟来自于变化地形环绕的四个隧道的SF6在风向变化过程中的浓度,尽管三个隧道的计算受低浓度值影响较大。特别的,在Nimomiya隧道附近的监测浓度和预测浓度有明显的不一致(Oettl et al., 2003)。在最近的研究中,模型准确预测了在城市周围的四个监测点的平均小时NOx浓度。在第五个监测点模型没有计算出浓度,因为这个监测点受没有包括在模型中的周围街道的影响很大。尽管模型有较好的性能,但它还没有普遍的应用,因为它有两个经验参数需要使用者根据交通量和隧道位置的细节调整,这可能是Nimomiya研究部分不一致的原因。

3.1.4 计算流体动力学模型

ARIA Local是CFD模型,已经用来计算城市环境中来自于公交车和火车的气体和粒子的真实扩散时间(Moon et al., 1997; Albergel andJasmin, 1998)。模型使用分辨率变化的网格,最小网格尺寸在1m以下,研究区域包括100万节点左右。但是,如果地形是重要的,它推荐使用水平栅格等距单元尺寸。三个不同的湍流模型应用是以或者方法以及变化的气体或者液体特性为基础的,大气稳

定度基于帕斯奎尔分类。污染源包括点,线,面和体的释放,以及从一个连续源或者容器释放产生的扩散。另外,气体中的流动性可以调节根据浮力或者密集的气体。尽管车辆影响导致的湍流没有隐藏的包含在模型中,但是可以通过调整模型参数包含在其中。化学转化可以通过后处理模块模拟。

MISKAM是应用在城市环境中的微尺度扩散模型。典型的,尽管可以模拟较大的功能区,但是模型应用的功能区大约300x300m,在每个方向上有60个非等距网格尺寸。建筑物可以作为块状体处理,模型没有考虑陡峭的地形,也不包括热效应,浮力释放和化学反应。可以用湍流混合因子模拟中性和稳定大气。任意污染物在风洞内产生的年平均浓度和模拟结果的比较显示很好的一致性

(R2=0.97),网格为1.2x1.6x0.14km,但是如果包括城市中心的浓度(R2=0.79)(Balczo et al., 2005) 一致性会降低。模型可以使用几何上改进的网格在一定区域提供更精确的分辨率。

MICRO-CALGRID(Stern and Yamartino,2001)是一个城市尺度光化学模型,运用由MISKAM模型发展的流场和湍流。除了上面讨论的MISKAM的性质,MICRO-CALGRID嵌接了交通感应发射模型MOBILEV(Fige, 1997)、水平和垂直平流、扩散项,一个很大的阻力是基于干沉降、化学转化项的参数化以及来自于CALGRID(在后面4.4节讨论)的SAPRC和CBM-IV。此外,模型通过调整模型细胞的总动能处理车辆导致的湍流,这个动能是车辆在空气中运动时产生的。虽然没有提供统计的和图解的证据,作者指出模型很好的再现了TSP的行为。

3.2 包括详细气溶胶动力学处理的局部气溶胶模型

GATOR(Jacobson, 1997)有一个选项,把移动尺寸或者固定尺寸粒子动力学模块和欧拉扩散方法耦合,计算城市(Jacobson, 1996) 和中尺度(Lu et al., 1997; Jacobson, 2001)环境气体、粒子的扩散。通过均相成核产生的新粒子,均相成核酸雨和来自于海洋环境(Pandis et al., 1994)的水反应参数化,近期的参数显示湿度大于60%(Fitzgerald et al., 1998)。通过半经验化的算法计算凝聚,半经验算法包括布朗扩散(Fuchs, 1964),对流布朗增强(Pruppacher and Klett,1997),重力收集(Pruppacher and Klett, 1997),湍流惯性运动和湍流切应(Saffmann and Turner, 1956)。凝结和蒸发可以从气体表面平衡计算,气体表面平衡是对开尔文效应关于由Zdanovskii, Stokes, Robinson (ZSR) (Robinson and Stokes, 1965)方法预测的气溶胶含水量的校正,包括无机物和有机物的水相解离。只有干沉降用阻力型方法预测(McRae et al.,1982; Russell et al., 1993; Seinfeld and Pandis,1998)。此外,模型会计算太阳辐射量,这对计算散射的光解离和粒子、气体的吸收曲线是很重要的。

Pohjola et al. (2003) 和 Gidhagen et al. (2004)最近公布了一些研究,在这些研究中他们分析了城市环境的粒子扩散。两个研究都是用的粒子动力学模型MONO32,但是分别与不同的扩散模型耦合来计算数量浓度,粒径分布和粒子的化学组成,一个是简单的烟羽模型OPSM(Vignati et al., 1999),一个是CFD 扩散模型STARCD(Gidhagen et al., 2004)。在Gidhagen et al. (2003)公路隧道的研究中,用STARCD很好的模拟了车辆交通排放的粒子数量浓度。当使用基于排放系数的速度时,预测结果和测量数据显示了很好的一致性。当使用恒定的

车辆排放系数时尽管相关性还是很好但是一致性会降低。一个对街道峡谷更深一步的研究(Gidhagen et al., 2004) 利用在顶层10m以上风速和风向,温度和降雨量监测值,使用MONO32准确计算了7到450nm间粒子的小时粒子浓度。这些模拟也说明交通导致的湍流对街道水平的扩散很重要。当存在或者不存在风速大于5m/s时交通导致的湍流,研究没有观察到显著的差异。但是,在2m/s的风速下,因为交通导致的湍流,粒子总的数量浓度会降低10000cm-3。

模型可以计算街道下风向的两边改变NOx浓度的主要因素。当交通导致湍流、凝聚及沉降导致气溶胶浓度变化包含在内时,粒子浓度和测量的小时浓度值显示出很好的一致性。

这些结果证实了MONO32和CFD模型耦合可以准确的预测车辆排放粒子的气溶胶动力学,并且指出当城市模拟使用基于排放系数的速度时计算有改进。近来,几个新的模型因为在操作方面的改进可以预测排气管排放的变化,例如加速或者已经开发的梯度变化,包括MicroFac (US EPA) (Singh et al., 2003), MEASURE (Georgia Tech), CMEM (UC Riverside)and TREFIC (ARIA Technologies).

Pohjola et al. (2003) 运用MONO32检验了排放25s后四个粒径分组的扩散。他们模拟了当存在或者不存在烟羽稀释时各个过程对粒子数量浓度的影响。二元成核、三元成核在模拟时间内都没有影响粒子数量。除去稀释,粒子数量浓度因为凝聚的影响会以数量级降低,因为凝结会有一点增加。当计算包括稀释时,虽然Aitken核型粒子数量降低,但是凝聚对总的粒子数量有负效应,当有机组分的气相浓度超过1010cm-3凝结是唯一重要的因素。

Korhonen et al. (2004) 研发了一个尺寸分离气溶胶动力学模型,UHMA,它的设计包括气溶胶动力学的处理和对新粒子形成、成长的关注。粒子是基于核的体积按尺寸分离的,核是由硫酸、水溶性有机物和不溶组分的变种混合组成的。虽然粒子组成会在不同尺寸分组间变化,但是粒子组成和每个的尺寸等级是一致的。粒子的有机馏分用化合物的集中描述来计算,化合物的性质可以由使用者依据特殊条件调整。

成核包括二元成核(Vehkamakiet al., 2002)和三元成核(Napari et

al.,2002),这取决于大气中氨的混合比,当氨浓度低于0.1ppt时是二元成核。

粒子的生长依赖于粒子的凝聚和凝结。模型中凝结的处理依赖于布朗运动,布朗运动是亚微米粒子凝结的主要因子,会把它作为每个时间段的粒径函数重新计算。

粒子生长也包括低挥发性的有机物凝结到粒子上,NanoKohler理论和氨的吸收,每个时间段的水分是基于粒子硫酸盐和水的可溶性有机质含量间的平衡,修正的ZSR方法是基于在芬兰各地进行的吸水性监测。

粒子干沉降通过尺寸的应用表现,尺寸的应用依赖于Rannik et al. (2003)对沉积速度的处理。

很多研究(Pirjola et al., 1999;Korhonen et al., 2003) 已经证明了UHMA 模型的性能。此外,Korhonen et al.(2004)研究有关于类似在森林上观测到的粒子形成(Makela et al.,2000)的一个新粒子形成的模型。他们指出模型用缩进运动中心的方法很好的预测了总的粒子数量,并且很好的描述了现阶段的生长。然而,模型过高预测了清晨总的数量浓度,并且没有预测出午后粒子数量浓度明显的降低,后面情况的失败是因为监测点气团的性质。

在Korhonen et al. (2003)第二阶段的研究中验证了两个模型化方法的性能,一般在大气模式中用来表示粒径分布,模拟新粒子的形成和生长。

他们发现混合组合方法比单分散方法更适合预测总的粒子浓度,这和在MONO32中应用一样,因为单分散方法不能模拟新形成粒子的凝聚。作者发现,当混合组合方法中的粒子数量降低时,精确解决粒径分布变化的能力会降低。但是在组合方法中越高粒径分辨率的方法需要的计算能力越大,这意味着这个方法不能用在大尺度扩散研究中。较简单的单分散方法可以用在大尺度研究中,正如作者论证的那样:它适合计算不同模式中总的粒子数量浓度和中等粒径的粒子。

2011年数学建模A题 城市表层土壤重金属污染分析

A题城市表层土壤重金属污染分析 摘要 随着城市经济的快速发展和城市人口的不断增加,人类活动对城市环境质量的影响日显突出。通过对城市土壤重金属的调查,应用数学方法对数据进行处理。得到城市环境质量的演变,已是人们日益关注的焦点。 对于问题一,用附件一中给出的数据,用matlab插值法建立三维模型,总共有9个图,一个是取样地点的地形图,另外八个是八种重金属元素的浓度分布图,通过模型图我们可以清楚的看到各种元素不同的空间分布。然后通过均值法,算出不同区域内各种重金属元素的污染程度。 对于问题二,通过对问题一结论的分析得出,生活区和工业区是污染比较厉害的地区。目前我国由于在重金属的开采、冶炼、加工过程中,造成不少重金属如铅、汞、镉、钴等进入大气、水、土壤引起严重的环境污染。人类生活中各种用品都含有不同量的重金属元素,比如说废旧电池,含有较多的汞、铬、锰、铅、镍、锌等重金属。它们通过自然和生物降解,随着雨水进入到土壤和河流当中。 对于问题三,根据前两问的结论分析重金属的传播特征,主要有从高海拔到低海拔,从高浓度区向低浓度区扩散。我们建立扩散模型,求出函数的极值,从而确定污染源的位置。 对于问题四,我们仔细分析了模型的优缺点。为更好地研究城市地质环境的演变模式,还应收集该地区的每年生活、工业等重要污染源的垃圾排放量,以及每年的生物降解量,降雨量对重金属元素扩散的影响,空气污染也应该考虑进去。有了这些数据以后建立因子分析法,回归分析,曲线拟合等模型解决问题。 关键词:插值法、均值法、扩散模型、因子分析、回归分析。 一、问题重述与分析 随着城市经济的快速发展和城市人口的不断增加,人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证,以及如何应用查证获得的海量数据资料开展城市环境质量评价,研究人类活动影响下城市地质环境的演变模

大气污染物扩散模式

第四章 大气扩散浓度估算模式 第一节 湍流扩散的基本理论 一 湍流 1.定义:大气的无规则运动 风速的脉动 风向的摆动 2.类型: 按形成原因 热力湍流:温度垂直分布不均(不稳定)引起,取决于大气稳定度 机械湍流:垂直方向风速分布不均匀及地面粗糙度引起 3.扩散的要素 风:平流输送为主,风大则湍流大 湍流:扩散比分子扩散快105~106倍 二 湍流扩散理论(主要阐述湍流与烟流传播及湍流与物质浓度衰减的关系) 1.梯度输送理论 通过与菲克扩散理论类比建立起来的(菲克定律:单位时间内通过单位断面上的物质的数量与浓度梯 度呈正比) 类比于分子扩散,污染物的扩散速率与负浓度梯度成正比 x C k F ??-= 式中,F — 污染物的输送通量 k — 湍流扩散系数 C — 污染物的浓度 X — 与扩散截面垂直的空间坐标(扩散过程的长度) x C ??— 浓度梯度 要求得各种条件下某污染物的时、空分布,由于边界条件往往很复杂,不能求出严格的分析解,只能是在特定的条件下求出近似解,再根据实际情况进行修正。 2.湍流统计理论 泰勒首先将统计理论应用在湍流扩散上 图4-1显示:从原点O 放出的粒子,在风沿着x 方向吹的湍流大气中扩散。粒子的位置用y 表示,则结论为: ①y 随时间变化,但其变化的平均值为零 ②若从原点放出很多粒子,则在x 轴上粒子的浓度最高,浓席分布以x 轴为对称轴,并符合正态分布。 萨顿实用模式:解决污染物在大气中扩散的实用模式 高斯模式:应用湍流统计理论得出正态分布假设下的扩散模式 3.相似理论 第二节 高斯扩散模式 一 坐标系的建立—右手坐标系

1.原点O :无界点源或地面源,O 为污染物的排放点 高架源,O 为污染物的排放点在地面上的投影点 补充:点源 高架源 连续源 固定源 线源 地面源 间歇源 流动源 面源 2.x 轴:正向为平均风向,烟流中心线与x 轴重合 3.y 轴:垂直于x 轴 4.z 轴:垂直于xoy 平面 二 高斯模式的有关假定 1.污染物浓度在y 、z 轴上的分布为正态分布; )2exp(21 )(22 y y y y f σπ σ-= )2exp(21 )(22 z z z z f σπ σ-= y σ,z σ— 分别为污染物在y 和z 方向上分布的标准差,m 2.全部高度风速均匀稳定,即风速u 为常数; 3.源强是连续均匀稳定的,源强Q 为定值; 4.扩散中污染物是守恒的,不考虑转化,即烟云在扩散过程中没有沉降、化合、分解及地面吸收、吸附作用发生; 0=??t C 5.在x 方向上,输送作用远远大于扩散作用,即 )(x C k x x C u x ????>>??; 6.地面足够平坦。

点污染源空气污染扩散模型

8 点、中午12 点、晚上9 点都没有排放气体,该怎么算,是不是需要找到一个关于时间t的函数,来计算多长时间之后污染还剩下多少 c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.^2)./((sigy+eps).^2)).*(exp(-0.5*(z-H).^2./((sigz+eps).^2))+exp(-0.5*(z+H).^2./((sigz+ eps).^2))); 这个函数对吗?该调用什么函数? 问题: 建立单污染源空气污染扩散模型,描述其对周围空气污染的动态影响规律。 现有河北境内某一工厂废气排放烟囱高50m,主要排放物为氮氧化物。早上9 点至下午 3 点期间的排放浓度为406.92mg/m3,排放速度为1200m3 /h;晚上10 点-凌晨4 点期间 的排放浓度为1160mg/m3,排放速度为5700m3 /h;通过你的扩散模型求解该工厂方圆51 公里分别在早上浓度8 点、中午12 点、晚上9 点空气污染分布和空气质量等级。 源代码 clear all clc [x,y]=meshgrid(0:20:5100,0:20:5100); Q=135.64; z=1.5; H=50; u=1.94; sigy=0.3914238*x.^0.865014; sigz=0.0757182*x.^1.00770; c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.^2)./((sigy+eps).^2)).*(exp(-0.5*(z-H).^2./((sigz+eps).^2))+exp(-0.5*(z+H).^2./((sigz+ eps).^2))); mesh(x,y,c); xlabel('X'),ylabel('Y'),zlabel('C'), clear all clc [x,y]=meshgrid(-5100:20:5100,-5100:20:5100); Q=1836.7; z=1.5; H=50; u=1.7; sigy=0.3914238*x.^0.865014; sigz=0.0757182*x.^1.00770; c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.^2)./((sigy+eps).^2)).*(exp(-0.5*(z-H).^2./((sigz+eps).^2))+exp(-0.5*(z+H).^2./((sigz+ eps).^2))); mesh(x,y,c); xlabel('X'),ylabel('Y'),zlabel('C'), 分享到: 2015-05-29 16:32 提问者采纳 clear all [x,y]=meshgrid(-51000:100:51000,-51000:100:51000); Q=135.64; z=1.5; H=50; u=1.94; sigy=0.3914238*x.^0.865014;

浅谈我国土壤重金属污染现状及修复技术

浅谈我国土壤重金属污染现状及修复技术 土壤是一个开放的缓冲动力学系统,承载着环境中50%~90%的污染负荷[1-2]。随着矿产资源开发、冶炼、加工企业等规模的扩大以及农业生产中农药、化肥、饲料等用量的增加和不合理的使用,致使土壤中重金属含量逐年累积,明显高于其背景值,造成生态破坏和环境质量恶化,对农业环境和人体健康构成严重威胁。重金属在土壤中移动性差、滞留时间长、难降解,可以通过生物富集作用和生物放大作用进入到农牧产品中[3],从而影响产出物的生长、产量和品质,潜在威胁人体健康[4]。本文对我国土壤重金属污染现状进行了简要分析,概述了土壤中重金属的来源,简单介绍了物理修复、化学修复和生物修复技术在土壤重金属污染修复方面的研究进展,以期为土壤重金属污染修复提供参考。 1我国土壤重金属污染现状 随着矿山开采、冶炼、电镀以及制革行业的蓬勃发展,一些企业盲目追逐经济利益,轻视环境保护,再加上农药、化肥、地膜、饲料添加剂等的大量使用,我国土壤中Pb、Cd、Zn等重金属的污染状况日益严重,污染面积逐年扩大,危害人类和动物的生命健康。据报道,2008年以来,全国已发生100余起重大污染事故,其中Pb、Cd、As等重金属污染事故达30多起。据2014年国家环境保护部和国土资源部发布的全国土壤污染状况调查公报显示,全国土壤环境总状况体不容乐观,部分地区土壤污染较重,耕地土壤环境质量堪忧,工矿业废弃地土壤环境问题突出。全国土壤总的点位超标率为16.1%,其中轻微、轻度、中度和重度污染点位比例分别为11.2%、2.3%、1.5%和1.1%。据农业部对我国24个省市、320个重点污染区约548 万hm2土壤调查结果显示,污染超标的大田农作物种植面积为60万hm2,其中重金属含量超标的农产品产量与面积约占污染物超标农产品总量与总面积的80%以上,尤

重金属污染物的传播特征

重金属污染来源、分布、治理方法 摘要:文章阐明了重金属污染物来源与分布,同时对国内外土壤重金属污染治理的研究工作做了系统的综述,提出了土壤中重金属污染物防治的环境矿物学新方法,利用环境矿物材料治理土壤重金属污染物的方法,具有成本低、效果好、无二次污染及有用金属可回收利用等优点,展现出广阔的环境矿物学研究与应用前景。并提醒人们要提高土壤质量意识,保护生态环境。 重金属系指密度4.0以上约60种元素或密度在5.0以上的45种元素。砷、硒是非金属,但是它的毒性及某些性质与重金属相似,所以将砷、硒列入重金属污染物范围内。环境污染方面所指的重金属主要是指生物毒性显著的汞、镉、铅、铬以及类金属砷,还包括具有毒性的重金属锌、铜、钴、镍、锡、钒等污染物。 随着全球经济化的迅速发展,含重金属的污染物通过各种途径进入土壤,造成土壤严重污染。土壤重金属污染可影响农作物产量和质量的下降,并可通过食物链危害人类的健康,也可以导致大气和水环境质量的进一步恶化。因此引起世界各国的广泛重视。目前,世界各国土壤存在不同程度的重金属污染,全世界平均每年排放Hg约1.5万t、Cu为340万t、Pb为500万t、Mn为1500万t、Ni为100万t。中国北方大城市的蔬菜基地和部分商品粮基地也存在着不同程度的重金属污染,如北京、天津、西安、沈阳、济南、长春、郑州等地;。 南方相对较轻,如福州、宁波、上海、武汉、成都等地。土壤重金属污染将会造成生态系统的严重破坏。从中国土壤资源状况看,到2000年底中国人均耕地仅为0.1 hm2,而且随着今后中国经济社会的发展如生态退耕、农业结构调整及自然灾害损毁等,土壤资源将进一步减少。因而如何有效地控制及治理土壤重金属的污染,改良土壤质量,将成为生态环境保护工作中十分重要的一项内容。 重金属污染原理 重金属,特别是汞、镉、铅、铬等具有显著和生物毒性。它们在水体中不能被微生物降解,而只能发生各种形态相互转化和分散、富集过程(即迁移)。重金属污染的特点是:(1)除被悬浮物带走的外,会因吸附沉淀作用而富集于排污口附近的底泥中,成为长期的次生污染源;(2)水中各种无机配位体(氯离子、硫酸离子、氢氧离子等)和有机配位体(腐蚀质等)会与其生成络合物或螯合物,导致重金属有更大的水溶解度而使已进入底泥的重金属又可能重新释放出来;(3)重金属的价态不同,其活性与毒性不同。其形态又随pH和氧化还原条件而转化。(4)在其危害环境方面的特点是:微量浓度即可产生毒性(一般为1~10毫克/升,汞、镉为0.01~0.001毫克/升);在微生物作用会转化为毒性更强的有机金属化合物(如洋-甲基汞);可被生物富集,通过食物链进入人体,造成慢性路线。亲硫重金属元素(汞、镉、铅、锌、硒、铜、砷等)与人体组织某些酶的巯基(-SH)有特别大的亲合力,能抑制酶的活性,亲铁元素(铁、镍)可在人体的肾、脾、肝内累积,抑制精氨酶的活性。六价铬可能是蛋白质和核酸的沉淀剂,可抑制细胞内谷胱甘肽还原酶,导致高铁血红蛋白,可能致癌,过量的钒和锰(亲岩元素)则能损害神经系统的机能。 本文主要从土壤中重金属污染物来源与分布、土壤中重金属污染物的现行治理方法入手,提出土壤中重金属污染物防治的环境矿物学新方法。旨在保护环境,提高土壤的环境质量。 1 土壤中重金属污染物来源与分布

2011年数学建模-重金属污染模型(特等奖作品)

兰州理工大学数学建模面朝大海 2011高教社杯全国大学生数学建模竞赛 承诺书 我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮 件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问 题。 我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他 公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正 文引用处和参考文献中明确列出。 我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反 竞赛规则的行为,我们将受到严肃处理。 我们参赛选择的题号是(从A/B/C/D中选择一项填写): 我们的参赛报名号为(如果赛区设置报名号的话): 所属学校(请填写完整的全名):兰州理工大学 参赛队员 (打印并签名) :1. 2. 3. 指导教师或指导教师组负责人 (打印并签名): 日期: 2011 年 9 月 11 日赛区评阅编号(由赛区组委会评阅前进行编号):

2011高教社杯全国大学生数学建模竞赛 编号专用页 赛区评阅编号(由赛区组委会评阅前进行编号): 全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):

城市表层土壤重金属污染分析 摘要 随着城市经济的快速发展和城市人口的不断增加,人类活动对城市环境质量的影响日显突出。本文根据大气中重金属污染物的传播特征,将高斯扩散模型应用于土壤重金属污染问题,建立数学模型并进行求解。 问题一中,利用题目所提供的重金属元素样本值,借助matlab软件,形象地绘制出了城市地形图以及八种重金属在该城市五个功能区的二维等高线分布图。对所收集的土壤样品分别采用单因子污染指数法与内梅罗(N.L.Nemerow)综合污染指数法对该城市分区进行污染程度评价。 问题二中,通过对所给数据的分析,利用主成分分析法,推断出了重金属污染的主要原因。 问题三中,通过数据分析,发现了金属污染物的传播特征,由此建立重金属高斯扩散模型。利用数形结合的方法最终确定了污染源的位置。 问题四中,通过对建模过程中考虑因素的全面性以及所忽略因素的得当性分析,评估了所建模型的优缺点,确定了更好地研究城市地质环境的演变模式所需要搜集的信息,并给出了更合理建立数学模型的方法。 关键词:重金属污染主成分分析法单因子指数法内梅罗综合污染指数法高斯模型污染源

大气污染扩散模型

第一节大气污染物的扩散 一、湍流与湍流扩散理论 1. 湍流 低层大气中的风向是不断地变化,上下左右出现摆动;同时,风速也是时强时弱,形成迅速的阵风起伏。风的这种强度与方向随时间不规则的变化形成的空气运动称为大气湍流。湍流运动是由无数结构紧密的流体微团——湍涡组成,其特征量的时间与空间分布都具有随机性,但它们的统计平均值仍然遵循一定的规律。大气湍流的流动特征尺度一般取离地面的高度,比流体在管道内流动时要大得多,湍涡的大小及其发展基本不受空间的限制,因此在较小的平均风速下就能有很高的雷诺数,从而达到湍流状态。所以近地层的大气始终处于湍流状态,尤其在大气边界层内,气流受下垫面影响,湍流运动更为剧烈。大气湍流造成流场各部分强烈混合,能使局部的污染气体或微粒迅速扩散。烟团在大气的湍流混合作用下,由湍涡不断把烟气推向周围空气中,同时又将周围的空气卷入烟团,从而形成烟气的快速扩散稀释过程。 烟气在大气中的扩散特征取决于是否存在 湍流以及湍涡的尺度(直径),如图5-7所示。 图5-7(a)为无湍流时,烟团仅仅依靠分子 扩散使烟团长大,烟团的扩散速率非常缓慢, 其扩散速率比湍流扩散小5~6个数量级;图5 -7(b)为烟团在远小于其尺度的湍涡中扩散, 由于烟团边缘受到小湍涡的扰动,逐渐与周边 空气混合而缓慢膨胀,浓度逐渐降低,烟流几乎呈直线向下风运动;图5-7(c)为烟团在与其尺度接近的湍涡中扩散,在湍涡的切入卷出作用下烟团被迅速撕裂,大幅度变形,横截面快速膨胀,因而扩散较快,烟流呈小摆幅曲线向下风运动;图5-7(d)为烟团在远大于其尺度的湍涡中扩散,烟团受大湍涡的卷吸扰动影响较弱,其本身膨胀有限,烟团在大湍涡的夹带下作较大摆幅的蛇形曲线运动。实际上烟云的扩散过程通常不是仅由上述单一情况所完成,因为大气中同时并存的湍涡具有各种不同的尺度。 根据湍流的形成与发展趋势,大气湍流可分为机械湍流和热力湍流两种形式。机械湍流是因地面的摩擦力使风在垂直方向产生速度梯度,或者由于地面障碍物(如山丘、树木与建筑物等)导致风向与风速的突然改变而造成的。热力湍流主要是由于地表受热不均匀,或因大气温度层结不稳定,在垂直方向产生温度梯度而造成的。一般近地面的大气湍流总是机械湍流和热力湍流的共同作用,其发展、结构特征及强弱决定于风速的大小、地面障碍物形成的粗糙度和低层大气的温度层结状况。 2. 湍流扩散与正态分布的基本理论 气体污染物进入大气后,一面随大气整体飘移,同时由于湍流混合,使污染物从高浓度区向低浓度区扩散稀释,其扩散程度取决于大气湍流的强度。大气污染的形成及其危害程度在于有害物质的浓度及其持续时间,大气扩散理论就是用数理方法来模拟各种大气污染源在

土壤重金属污染状况及修复

土壤重金属污染状况及修复 中文摘要:重金属污染因具有毒性、易通过食物链在植物,动物和人体内累积,对生态环境和人体健康构成严重威胁。随着工业快速发展、农药及化肥的广泛使用,农田土壤重金属污染越来越严重,研究农田土壤重金属污染现状及修复技术对农产品安全具有重要意义。综合国内外农田土壤重金属污染状况,农田土壤重金属污染主要来源于固体废弃物堆放及处置、工业废物大气沉降、污水农灌和农用物质的不合理施用。该文综述了国内外有关农田重金属污染土壤修复技术(物理修复、化学修复、生物修复、农业生态和联合修复)的研究进展,并针对各种修复方法,阐述了其原理、修复条件、应用实例及其优缺点,重点论述了植物修复的机理和应用,提出了草本与木本联合修复可有效提高农田土壤重金属复合污染的修复效率,为农田土壤土壤重金属复合污染修复提出了新的途径。最后在对已有研究分析的基础上,提出了联合修复技术(如生物联合技术、物理化学联合技术和物理化学—生物联合技术)可以在一定程度上克服使用单一修复手段存在的缺点,可提高复合污染的修复效率、降低修复成本,未来应深入探索联合修复技术间的相互作用机理,以期为农田土壤重金属综合治理与污染修复提供科学依据。 关键词:农田土壤;重金属;污染;修复技术 Abstract; Heavy metal pollution caused by toxic, easily in the food chain through plants, animals and humans in vivo accumulation of the ecological environment and pose a serious threat to human health. With the rapid development of industry, the widespread use of pesticides and fertilizers, agricultural soil heavy metal pollution is getting worse, research Soil Heavy Metal Pollution and Remediation Technology is important for the safety of agricultural products. Comprehensive Farmland Soil Heavy Metal Contamination at home and abroad, mainly from heavy metals in soils contaminated solid waste deposits and disposal of industrial waste atmospheric deposition, sewage unreasonable application of agricultural irrigation and agricultural materials. This paper reviews the related farmland abroad Heavy Metal Contaminated Soil Research Progress (physical restoration, chemical remediation, bioremediation, ecological agriculture and bioremediation) repair, and for a variety of repair methods, described its principle, to repair the condition, application examples its advantages and disadvantages, Focuses on the mechanism and application of phytoremediation, herbaceous and woody proposed bioremediation can effectively improve the efficiency of heavy metals in soils repair compound contaminated soil farmland soil heavy metals contamination fixes proposed a new way. Finally, the existing research and analysis based on the proposed joint repair techniques (such as bio-technology joint, joint technical and physical chemistry physical chemistry - Biotechnology United Technologies) can overcome the disadvantages of using a single repair means exist to some extent, can improve repair efficiency combined pollution, reduce repair costs, Future should further explore the mechanism of interaction between the United repair techniques, with a view to the comprehensive management of heavy metals in soils and pollution remediation provide a scientific basis. Keywords: Soil; heavy metal; pollution; repair technology 1 土壤中重金属的污染现状 土壤作为开放的缓冲动力学体系,在与周围的环境进行物质和能量的交换过程中,不可避免地会有外源重金属进入这个体系! 重金属对土壤的主要污染途径是工业废渣、废气 中重金属的扩散、沉降、累积,含重金属废水灌溉农田,以及含重金属农药、磷肥的大量施用! 外来重金属多富集在土壤的表层!.工业生产上重金属释放到环境中的主要途径有采矿、冶炼、燃

大气污染物扩散高斯模型模拟

大气污染物扩散的高斯模型模拟:可视化模拟点源大气污染的扩散Gaussian Atmospheric Dispersion Model 突发性大气污染事故时有发生,对大气污染扩散进行模拟和分析,有利于减小事故的危害,减轻人员伤亡和财产损失。高斯扩散模型是国际原子能机构(IAEA)推荐使用于重气云扩散模拟的数学模型,该模型在非重气云扩散的应用日益广泛。高斯扩散模型是描述大气对有害气体的输移、扩散和稀释作用的物理或数学模型,是进行灾害预测和救援指挥的有力手段之一。 高斯扩散模型 高斯模型又分为高斯烟团模型和高斯烟羽模型。大气污染物泄漏分为瞬时泄漏和连续泄漏,瞬时泄漏是指污染物泄放的时间相对于污染物扩散的时间较短如突发泄漏等的情形,连续泄漏则是指污染物泄放的时间较长的情形。瞬时泄漏采用高斯烟团模型模拟,而连续泄漏采用高斯模型烟羽模型模拟。高斯模型适用于非重气云气体,包括轻气云和中性气云气体。要求气体在扩散过程中,风速均匀稳定。 在高斯烟团模型中,选择风向建立坐标系统,即取泄漏源为坐标原点,x轴指向风向,y轴表示在水平面内与风向垂直的方向,z轴则指向与水平面垂直的方向,具体公式见式: (mg/s); x、y、z轴上的扩散系数,需根据大气稳定度选择参数计算得到(m);x、y、z表示x、y、z上的坐标值(m);u 表示平均风速(m/s);t表示扩散时间(s);H 表示泄漏源的高度(m)。 同理,高斯烟羽模型的表达式如: 技术方法 若用高斯模型算出空间每一个点在一个时刻的污染浓度,这个计算量是很大的。因此所设计的系统一般都是采用先进行图层网格化,由高斯模型计算出有限个网格点的上的污染物浓度,在进行空间内插得到面上每一个点的污染物浓度,并由此得到污染物浓度的等值线。整个过程的示意图如图所示

土壤重金属污染的危害及修复教学提纲

土壤重金属污染的危 害及修复

土壤重金属污染的危害及修复 摘要:土壤重金属污染问题越来越引起人们关注,阐明了土壤中重金属污染的来源、污染情况及造成的危 害,主要综述了目前国内专家、学者对土壤污染及生物修复的研究进展,结合我国具体情况,提出一些自己的看法. 关键词:土壤;重金属污染;生物修复 土壤重金属污染是指人类活动将重金属加入到土壤中,致使土壤中重金属含量明显高于原有含量,并造成生态环境恶化 的现象[1].环境污染方面所指的重金属主要指对农作物和人畜生物毒性显著的Hg、Cd、Pb、Cr、以及类金属As,还包括具有毒 性的Zn、Cu、Co,Ni、Sn、V等污染物,后者在常量下对作物和人体是营养元素,过量时则出现危害.加强土壤污染的化学及生态 研究对推动绿色食品和生态农业的发展具有重要意义. 1土壤中重金属元素的来源和污染状况 除了来自于土母质本身的重金属,土壤重金属污染主要来自于人类活动.研究表明:Pb、Cd、Hg、As与大气污染有直接关 系[2].来源于象汽车含铅汽油燃烧排放的尾气、工农业生产、汽车轮胎磨损产生的大量含重金属的气体,它们经过自然沉降和 雨淋进入土壤.公路、铁路两侧土壤中的重金属污染主要是Pb、Cr、Zn,Cu、Co、Cd等,大气汞的干湿沉降也可引起土壤中汞含 量的增高.

城市大量的工业废水流入河道,其中含有的许多重金属离子,随着污水灌溉、污泥施肥而进入土壤.太原、淮阳污灌区土 壤中重金属的含量自污灌以来逐年增高.广州市郊污灌区农田中Pb、Cd、Hg、Cr、As等重金属污染超过临界值,残留超标率分 别达16%、100%、68%、16%和52%[3、4].研究还表明:用城市污水污泥改良土壤,重金属Hg、Pb、Cr的含量明显增加,青菜中 的Pb、Zn、Cu、Cd、Ni也增加[5]. 胡永定[6]通过对徐州荆马河区域土壤重金属污染成因的分析和研究,发现Cd是由垃圾施用和农灌引起的,Pb、Zn、Cu、Cr 是由垃圾施用引起的,As是农田灌溉引起的,Hg是各种途径都有.另外城市生活垃圾、车辆废弃物、垃圾堆放场附近土壤中重 金属的含量都高于当地土壤背景值,如北京郊区某垃圾场周边土壤中Cd含量是对照组的8倍.金属矿山的开采、有色金属的 冶炼排放的废水、重金属冶炼矿渣的堆放,工厂烟囱的排放物等,随着降雨淋溶被带入水环境或直接进入土壤,都会成为土壤 重金属的来源.许多研究表明:随着磷肥、复合肥的大量施用,土壤有效镉的含量在不断增加,作物吸收镉量也相应增加.据马 耀华等对上海地区菜园土研究发现:施肥后,Cd的含量从0.1mgkg- 1上升到 0.32mg kg- 1.魏秀国等人通过对广州市蔬菜地 土壤重金属污染状况调查及评价发现:铅污染最为普遍,其次是砷污染;就污染的程度而言,镉污染最为严重,其次为砷[7].

第五章 大气污染扩散

第五章大气污染扩散 第一节大气结构与气象 有效地防止大气污染的途径,除了采用除尘及废气净化装置等各种工程技术手段外,还需充分利用大气的湍流混合作用对污染物的扩散稀释能力,即大气的自净能力。污染物从污染源排放到大气中的扩散过程及其危害程度,主要决定于气象因素,此外还与污染物的特征和排放特性,以及排放区的地形地貌状况有关。下面简要介绍大气结构以及气象条件的一些基本概念。 一、大气的结构 气象学中的大气是指地球引力作用下包围地球的空气层,其最外层的界限难以确定。通常把自地面至1200 km左右范围内的空气层称做大气圈或大气层,而空气总质量的98.2%集中在距离地球表面30 km以下。超过1200 km的范围,由于空气极其稀薄,一般视为宇宙空间。 自然状态的大气由多种气体的混合物、水蒸气和悬浮微粒组成。其中,纯净干空气中的氧气、氮气和氩气三种主要成分的总和占空气体积的99.97%,它们之间的比例从地面直到90km高空基本不变,为大气的恒定的组分;二氧化碳由于燃料燃烧和动物的呼吸,陆地的含量比海上多,臭氧主要集中在55~60km高空,水蒸气含量在4%以下,在极地或沙漠区的体积分数接近于零,这些为大气的可变的组分;而来源于人类社会生产和火山爆发、森林火灾、海啸、地震等暂时性的灾害排放的煤烟、粉尘、氯化氢、硫化氢、硫氧化物、氮氧化物、碳氧化物为大气的不定的组分。 大气的结构是指垂直(即竖直)方向上大气的密 度、温度及其组成的分布状况。根据大气温度在垂直 方向上的分布规律,可将大气划分为四层:对流层、 平流层、中间层和暖层,如图5-1所示。 1. 对流层 对流层是大气圈最靠近地面的一层,集中了大气 质量的75%和几乎全部的水蒸气、微尘杂质。受太阳 辐射与大气环流的影响,对流层中空气的湍流运动和 垂直方向混合比较强烈,主要的天气现象云雨风雪等 都发生在这一层,有可能形成污染物易于扩散的气象 条件,也可能生成对环境产生有危害的逆温气象条件。 因此,该层对大气污染物的扩散、输送和转化影响最大。 大气对流层的厚度不恒定,随地球纬度增高而降低,且与季节的变化有关,赤道附近约为15km,中纬度地区约为10~12 km,两极地区约为8km;同一地区,夏季比冬季厚。一般情况下,对流层中的气温沿垂直高度自下而上递减,约每升高100m平均降低0.65℃。 从地面向上至1~1.5 km高度范围内的对流层称为大气边界层,该层空气流动受地表影响

重金属污染物的传播特征,以及产生污染的原因

重金属污染物在土壤中的传播特征 重金属系指密度4.0以上约60种元素或密度在5.0以上的45种元素。砷、硒是非金属,但是它的毒性及某些性质与重金属相似,所以将砷、硒列入重金属污染物范围内。环境污染方面所指的重金属主要是指生物毒性显著的汞、镉、铅、铬以及类金属砷,还包括具有毒性的重金属锌、铜、钴、镍、锡、钒等污染物。随着全球经济化的迅速发展,含重金属的污染物通过各种途径进入土壤,造成土壤严重污染。土壤重金属污染可影响农作物产量和质量的下降,并可通过食物链危害人类的健康,也可以导致大气和水环境质量的进一步恶化。因此引起世界各国的广泛重视。目前,世界各国土壤存在不同程度的重金属污染,全世界平均每年排放Hg约1.5万t、Cu为340万t、Pb为500万t、Mn为1500万t、Ni 为100万t[1]。中国北方大城市的蔬菜基地和部分商品粮基地也存在着不同程度的重金属污染,如北京、天津、西安、沈阳、济南、长春、郑州等地;。 南方相对较轻,如福州、宁波、上海、武汉、成都等地。土壤重金属污染将会造成生态系统的严重破坏。从中国土壤资源状况看,到2000年底中国人均耕地仅为0.1 hm2,而且随着今后中国经济社会的发展如生态退耕、农业结构调整及自然灾害损毁等,土壤资源将进一步减少。因而如何有效地控制及治理土壤重金属的污染,改良土壤质量,将成为生态环境保护工作中十分重要的一项内容。 本文主要从土壤中重金属污染物来源与分布、土壤中重金属污染物的现行治理方法入手,提出土壤中重金属污染物防治的环境矿物学新方法。旨在保护环境,提高土壤的环境质量。 1 土壤中重金属污染物来源与分布 土壤中重金属的来源是多途径的,首先是成土母质本身含有重金属,不同的母质、成土过程所形成的土壤含有重金属量差异很大。此外,人类工农业生产活动,也造成重金属对大气、水体和土壤的污染。 1.1 大气中重金属沉降 大气中的重金属主要来源于工业生产、汽车尾气排放及汽车轮胎磨损产生的大量含重金属的有害气体和粉尘等。它们主要分布在工矿的周围和公路、铁路的两侧。大气中的大多数重金属是经自然沉降[2]和雨淋沉降进入土壤的。如瑞典中部Falun市区的铅污染[3],它主要来自于市区铜矿工业厂、硫酸厂、油漆厂、采矿和化学工业产生大量废物,由于风的输送,这些细微颗粒的铅,从工业废物堆扩散至周围地区。南京某生产铬的重工业厂[4]铬污染叠加已超过当地背景值4.4倍,污染以车间烟囱为中心,范围达1.5 km2,污染范围最大延伸下限1.38 km。俄罗斯的一个硫酸生产厂[5]也是由工厂烟囱排放造成S、V、As的污染。 公路、铁路两侧土壤中的重金属污染,主要是Pb、Zn、Cd、Cr、Co、Cu的污染为主。它们来自于含铅汽油的燃烧,汽车轮胎磨损产生的含锌粉尘等。它们成条带状分布,以公路、铁路为轴向两侧重金属污染强度逐渐减弱;随着时间的推移,公路、铁路土壤重金属污染具有很强的叠加性。在宁—杭公路南京段[6]两侧的土壤形成Pb、Cr、Co污染晕带,且沿公路延长方向分布,自公路向两侧污染强度减弱。在宁—连一级公路淮阴段[7]两侧的土壤铅含量增高,向两侧含量逐渐降低,且在地表0~30 cm铅的含量较高。在法国索洛涅地区A71号高速公路[8]沿途严重污染重金属Pb、Zn、Cd,其沉降粒子浓度超过当地土壤背景值

重金属污染源的数学建模

承诺书 我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。 我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃处理。 我们参赛选择的题号是(从A/B/C/D中选择一项填写): A 我们的参赛报名号为(如果赛区设置报名号的话): 所属学校(请填写完整的全名):中国人民解放军国防科学技术大学参赛队员(打印并签名) :1. 李腾骥 2. 袁雪强 3. 赵泉朴 指导教师或指导教师组负责人(打印并签名):杨文强 日期:2011 年9 月9 日赛区评阅编号(由赛区组委会评阅前进行编号):

编号专用页 赛区评阅编号(由赛区组委会评阅前进行编号): 全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):

基于MCMC的城市表层土壤重金属污染分析 摘要:本文针对某城市城区表层土壤重金属污染分析的问题,从各分区重金属元素污染的空间分布、污染程度、污染原因及传播途径等多方面进行了全面的统计与研究,同时采用贝叶斯推理和马尔科夫链蒙特卡洛法(MCMC),利用MATLAB进行编程,对所设计的模型进行求解,得到了各重金属污染源的地理位置坐标,最后对城市地质环境演变模式的问题进行了探究,得出了初步的结论。 首先,根据所给数据信息,利用MATLAB及EXCEL作图,得到了该地区的地势图及各重金属污染程度的分布图。 对于问题1,通过统计各重金属的重污染点在各区域的分布和计算不同区域各中金属污染程度之间的相关系数,得到了各重金属的大致空间分布,并采用改进的内梅罗指数法,得出不同区域重金属的污染程度。 对于问题2,在问题1所得数据及结论的基础上,结合实际情况,对各重金属污染的主要原因进行了较为全面的分析。 对于问题3,通过对污染的传播特征的研究,基于城市的地质特性,得到空气传播为污染主要传播方式的结论,并由此建立关于气体传播的概率模型,采用贝叶斯推理和马尔科夫链蒙特卡洛法,利用MATLAB进行50000次抽样,得到收敛的坐标值,从而计算出各重金属污染源的地理位置坐标分别为:As (2742,7294)、Cd(2382,3693)、Cr(3470,2309)、Cu(2707,2295)、Hg(2708,2294)、Ni(1647,2729)、Pb(2882,3617)、Zn(4152,2299)。 对于问题4,以题中城市土壤重金属污染的演变规律为例进行了一定深度的探究,并对所需信息及模型设计方案进行了初步的讨论。 关键词:内梅罗指数法相关系数贝叶斯推理MCMC法MH算法 一、问题的提出与分析 (一)背景介绍 随着城市经济的快速发展和城市人口的不断增加,人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证,以及如何应用查证获得的海量数据资料开展城市环境质量评价,研究人类活动影响下城市地质环境的演变模式,日益成为人们关注的焦点。 按照功能划分,城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等,分别记为1类区、2类区、……、5类区,不同的区域环境受人类活动影响的程度不同。 现对某城市城区土壤地质环境进行调查。为此,将所考察的城区划分为间距1公里左右的网格子区域,按照每平方公里1个采样点对表层土(0~10 厘米深度)进行取样、编号,并用GPS记录采样点的位置。应用专门仪器测试分析,获得了每个样本所含的多种化学元素的浓度数据。另一方面,按照2公里的间距在那些远离人群及工业活动的自然区取样,将其作为该城区表层土壤中元素的背景值。

重金属传播特征讲解

重金属传播特征 重金属原义是指比重大于 5的金属,包括金、银、铜、铁、铅等,重金属在人体中累积达到一定程度,会造成慢性中毒。对什么是重金属目前尚无严格的定义,化学上跟据金属的密度把金属分成重金属和轻金属,常把密度大于 4.5g/cm3的金属称为重金属。如:金、银、铜、铅、锌、镍、钴、铬、汞、镉等大约 45种。 从环境污染方面所说的重金属是指:汞、镉、铅、铬以及类金属砷等生物毒性显著的重金属。 对人体毒害最大的有 5种:铅、汞、铬、砷、镉。这些重金属在水中不能被分解,人饮用后毒性放大,与水中的其他毒素结合生成毒性更大的有机物或无机物。 在领域中,重金属主要是指对生物有明显毒性的金属元素或类金属元素,如汞、镉、铅、铬、锌、铜、钴、镍、锡、砷等,此类污染物不易被微生物降解。 随着全球经济化的迅速发展, 含重金属的污染物通过各种途径进入土壤, 造成土壤严重污染。土壤重金属污染可影响农作物产量和质量的下降, 并可通过食物链危害人类的健康, 也可以导致大气和水环境质量的进一步恶化。因此引起世界各国的广泛重视。目前, 世界各国土壤存在不同程度的重金属污染,全世界平均每年排放Hg 约 1.5万 t 、 Cu 为 340万 t 、 Pb 为 500万 t 、 Mn 为 1500万 t 、 Ni 为 100万t 。中国北方大城市的蔬菜基地和部分商品粮基地也存在着不同程度的重金属污染, 如北京、天津、西安、沈阳、济南、长春、郑州等地; 。南方相对较轻,如福州、宁波、上海、武汉、成都等地。土壤重金属污染将会造成生态系统的严重破坏。从中国土壤资源状况看, 到 2000年底中国人均耕地仅为 0.1 hm2, 而且随着今后中国经济社会的发展如生态退耕、农业结构调整及自然灾害损毁等, 土壤资源将进一步减少。因而如何有效地控制及治理土壤重金属的污染,改良土壤质量 , 将成为生态环境保护工作中十分重要的一项内容。 重金属污染原理

重金属污染

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。 我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃处理。 我们参赛选择的题号是(从A/B/C/D中选择一项填写): 我们的参赛报名号为(如果赛区设置报名号的话): 所属学校(请填写完整的全名): 参赛队员(打印并签名) :1. 2. 3. 指导教师或指导教师组负责人(打印并签名): 日期: 2012 年 8 月15 日赛区评阅编号(由赛区组委会评阅前进行编号):

2010高教社杯全国大学生数学建模竞赛 编号专用页 赛区评阅编号(由赛区组委会评阅前进行编号): 全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):

城市表层土壤重金属污染分析 摘要 本文运用Matlab、SPSS及Excel软件,综合分析该城区各区重金属污染程度及其重金属污染的主要原因,并根据重金属的传播特征建立数学模型,寻找出污染源,同时更进一步修改模型以更好研究城市地质环境的演变模式。 针对问题一,利用Matlab作出各重金属浓度空间分布图,可看出重金属浓度大多呈高浓度向外逐步扩散现象。同时,利用单因素污染指数公式及污染等级规则,通过模糊综合评价模型,得出各类区域重金属污染程度如下: 针对问题二,利用SPSS软件对各重金属变量进行R型聚类分析,将重金属污染分为四类:Cr、Ni、Cu为一类,Cd、Pb、Zn为一类,而As、Hg各为一类。并计算出各区域各单金属污染指数进行对照分析,说明该城区重金属污染主要来源于工业活动及其生产品和交通主干道污染。 针对问题三,我们根据重金属物质在土壤中扩散的特点,利用有衰减的扩散模型,构造出浓度与坐标的关系,并利用matlab编程,多次试验,合理筛选数据,通过回归分析的方法,求解参数,从而得到不同重金属的不同污染源。 针对问题四,由于问题三求解的局限性,我们在模型三的基础上,增加了时间、水文、人类活动等因素,得出一个更为具体并符合实际的模型。有利于计算土壤中重金属的浓度。 关键词单因素污染指数模糊综合评价 R型聚类分析有衰减的扩散模型扩散规律

相关文档
最新文档