扩散和沉降耦合作用下霾气溶胶分布的研究_方爱平

扩散和沉降耦合作用下霾气溶胶分布的研究

方爱平　田蓬勃　应　钺　蔡佳音　苟明岳　郭启隆　宋　倩　杨宇驰　喻有理　王小力

（西安交通大学理学院，陕西西安　７１００４９

）摘　要　通过气溶胶的特性，从扩散、沉降、凝聚三个方面对霾微粒进行分析．发现对于霾微粒，

沉降速度采用孤粒子沉降速度是完全可行的；凝聚作用对一般的霾分散系作用很小以至于可以忽略．

在这一前提下，采用对流扩散方程描述扩散和沉降的耦合作用，并给出了积分形式的解析解．参考高斯烟羽模型，给出一种较为典型的初值条件，并在污染物接触地面完全反射的假设下通过数值方法模拟未来的污染物浓度演化．模拟结果显示：（１）在无风的条件下，霾消散的速度是相当缓慢的；（２

）霾的演化中起最主要作用的是沉降．关键词　扩散；沉降；凝聚；扩散系数；Ｓｔｏｋｅｓ速度；

高斯烟羽模型ＳＴＵＤＹ　ＯＮ　ＴＨＥ　ＤＩＳＴＲＩＢＵＴＩＯＮ　ＯＦ　ＡＥＲＯＳＯＬ　ＨＡＺＥ　

ＵＮＤＥＲＴＨＥ　ＣＯＵＰＬＥＤ　ＡＣＴＩＯＮ　ＯＦ　ＤＩＦＦＵＳＩＯＮ　ＡＮＤ　

ＳＵＢＳＩＤＥＮＣＥＦａｎｇ　Ａｉｐｉｎｇ　Ｔｉａｎ　Ｐｅｎｇｂｏ　Ｙｉｎｇ　Ｙｕｅ　Ｃａｉ　Ｊｉａｙｉｎ　Ｇｏｕ　ＭｉｎｇｙｕｅＧｕｏ　Ｑｉｌｏｎｇ　Ｓｏｎｇ　Ｑｉａｎ　Ｙａｎｇ　Ｙｕｃｈｉ　Ｙｕ　Ｙｏｕｌｉ　Ｗａｎｇ　

Ｘｉａｏｌｉ（Ｓｃｈｏｏｌ　ｏｆ　Ｓｃｉｅｎｃｅ，Ｘｉ’ａｎ　Ｊｉａｏｔｏｎｇ　

Ｕｎｉｖｅｒｓｉｔｙ，Ｘｉ’ａｎ，Ｓｈａａｎｘｉ　７１００４９）Ａｂｓｔｒａｃｔ　Ｉｎ　ｔｈｉｓ　ｐａｐｅｒ，ａｃｃｏｒｄｉｎｇ　ｔｏ　ｔｈｅ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏｆ　ｔｈｅ　ａｅｒｏｓｏｌ，ｄｉｆｆｕｓｉｏｎ，ｓｕｂｓｉｄｅｎｃｅａｎｄ　ｃｏａｇｕｌａｔｉｏｎ　ｔｈｅｓｅ　ｔｈｒｅｅ　ａｓｐｅｃｔｓ　ｏｆ　ｈａｚｅ　ｐａｒｔｉｃｌｅｓ　ｗｅｒｅ　ａｎａｌｙｚｅｄ．Ｗｅ　ｆｏｕｎｄ　ｔｈａｔ　ｆｏｒ　ｈａｚｅｐａｒｔｉｃｌｅｓ，ｔｈｅ　ｕｓｅ　ｏｆ　ｓｏｌｉｔａｒｙ　ｐａｒｔｉｃｌｅ　ｓｅｄｉｍｅｎｔａｔｉｏｎ　ｖｅｌｏｃｉｔｙ　ｉｓ　ｆｅａｓｉｂｌｅ；ｃｏｈｅｓｉｖｅ　ｅｆｆｅｃｔ　ｏｎ　ｔｈｅｇｅｎｅｒａｌ　ｈａｚｅ　ｄｉｓｐｅｒｓｉｏｎ　ｅｆｆｅｃｔ　ｉｓ　ｓｏ　ｓｍａｌｌ　ｔｈａｔ　ｃａｎ　ｂｅ　ｉｇ

ｎｏｒｅｄ．Ｉｎ　ｔｈｉｓ　ｃｏｎｔｅｘｔ，ｔｈｅ　ｕｓｅ　ｏｆ　ｃｏｕ－ｐｌｉｎｇ　ｄｉｆｆｕｓｉｏｎ　ａｎｄ　ｃｏｎｖｅｃｔｉｏｎ－ｄｉｆｆｕｓｉｏｎ　ｅｑｕａｔｉｏｎ　ｄｅｓｃｒｉｂｅｓ　ｔｈｅ　ｓｅｔｔｌｅｍｅｎｔ，ａｎｄ　ｔｈｅ　ａｎａｌｙｔｉｃ　ｓｏ－ｌｕｔｉｏｎ　ｏｆ　ｔｈｅ　ｉｎｔｅｇｒａｌ　ｆｏｒｍ．Ｒｅｆｅｒｅｎｃｅ　Ｇａｕｓｓｉａｎ　ｐｌｕｍｅ　ｍｏｄｅｌ　ｔｏ　ｇｉｖｅ　ａ　ｍｏｒｅ　ｔｙｐ

ｉｃａｌ　ｉｎｉｔｉａｌ　ｃｏｎ－ｄｉｔｉｏｎｓ，ａｎｄ　ｕｎｄｅｒ　ｔｈｅ　ａｓｓｕｍｐｔｉｏｎ　ｔｈａｔ　ｔｈｅ　ｇｒｏｕｎｄ　ｃａｎ　ｃｏｍｐｌｅｔｅｌｙ　ｒｅｆｌｅｃｔｅｄ　ｔｈｅ　ｐｏｌｌｕｔａｎｔ，ｔｈｅｆｕｔｕｒｅ　ｅｖｏｌｕｔｉｏｎ　ｏｆ　ｔｈｅ　ｐｏｌｌｕｔａｎｔｓ　ｃａｎ　ｂｅ　ｓｏｌｖｅｄ　ｂｙ　

ｎｕｍｅｒｉｃａｌ　ｍｅｔｈｏｄｓ．Ｔｈｅ　ｓｉｍｕｌａｔｉｏｎ　ｒｅｓｕｌｔｓｓｈｏｗ　ｔｈａｔ：（１）ｉｎ　ｔｈｅ　ｗｉｎｄｌｅｓｓ　ｃｏｎｄｉｔｉｏｎｓ，ｈａｚｅ　ｄｉｓｓｉｐａｔｅｄ　ｓｐｅｅｄ　ｉｓ　ｖｅｒｙ　ｓｌｏｗ；（２）Ｔｈｅ　ｓｅｄｉ－ｍｅｎｔａｔｉｏｎ　ｐｌａｙ

ｓ　ｔｈｅ　ｍｏｓｔ　ｉｍｐｏｒｔａｎｔ　ｒｏｌｅ　ｉｎ　ｈａｚｅ　ｅｖｏｌｕｔｉｏｎ．Ｋｅｙ　ｗｏｒｄｓ　ｄｉｆｆｕｓｉｏｎ；ｓｕｂｓｉｄｅｎｃｅ；ｃｏｈｅｓｉｏｎ；ｄｉｆｆｕｓｉｏｎ　ｃｏｅｆｆｉｃｉｅｎｔ；Ｓｔｏｋｅｓ　ｓｐｅｅｄ；Ｇａｕｓｓｉａｎ　ｐｌｕｍｅ　ｍｏｄｅｌ收稿日期：２０１６－０５－

０１作者简介：方爱平，女，西安交通大学理学院讲师，博士，主要从事量子光学和大学物理的教学与研究工作．通讯作者：田蓬勃，男，副教授，主要从事大学物理和热学教学科研工作，研究方向为光学．ｐ

ｂｔｉａｎ＠１６３．ｃｏｍ引文格式：方爱平，田蓬勃，应钺，等．扩散和沉降耦合作用下霾气溶胶分布的研究［Ｊ］．物理与工程，２０１６，２６（４）：１０４－

１０８． 近年来，

能源消耗量增大，燃烧化石燃料的过程中排放到空气中的化学物质以及颗粒状的污染物（如尘、煤等）不断增加，造成了各种形式的大气污染．其中，霾这一现象尤其引人关注．霾，也称阴霾、灰霾，是指大量烟、尘等微粒悬浮而形成的浑浊现象．霾的核心物质是空气中悬浮的灰尘颗粒，气象学上称为气溶胶颗粒．霾中含有数百种大气化学颗粒物质，可引起多种疾病．此外，霾影响大

气能见度，常引发交通事故．为了治理霾，在实际工作中，需要根据污染源的位置、高度、气象情况等条件正确推算出污染物在大气中的浓度．再以此来制定合理的治理措施．

本文从气溶胶的通性［１

］入手分析影响雾霾分

散系的因素．考虑扩散、沉降、凝聚以及周围环境的影响，

通过分析论证了采用对流扩散方程描述霾分布演化的可行性．以高斯正态烟羽模式为背景，研

究单一污染源排放出的污染物，推算出了周围地区的污染物浓度分布随时间的演化，以供后续使用．１　模型建立

我们基于两个基本假设：（１）

考察的空间足够大以至于可以认为是无界的；（２）环境中无风且温度均匀，通过气溶胶的通性，给出描述扩散，沉降，凝聚的方程．并通过分析合理的简化模型．１．１　控制气溶胶运动的方程

１

）扩散扩散是物质分子从高浓度区域向低浓度区域转移，直到均匀分布的现象．由菲克定律，扩散的速率与物质的浓度梯度成正比．

故在流场瓫内，有对流扩散方程［

２］ ｃ

ｔ＋瓫·!

ｃ＝Ｄ!２ｃ（１）其中，ｃ为质量浓度；Ｄ为扩散系数，

可以由下式确定［

２，３］Ｄ＝

ｋＢ

Ｔ　Ｃｕ３πμ

ｄｐ（２

）其中，ｋＢ为玻尔兹曼常数；Ｔ为绝对温度；μ为动力

黏度系数，在空气中μ＝１．８１×１０－５

Ｐａ·ｓ；ｄｐ为粒

子粒度；Ｃｕ为库宁汉滑移修正系数，

常温常压下，施特劳斯给出的库宁汉修正系数为［

３］

Ｃｕ＝１＋

２λｄｐ１．２５７＋０．４ｅｘｐ－１．１ｄｐ２（）［

］

λ

（３

）其中λ为平均自由程［２］

，有

λ＝λ０

ｐ０

（）ｐＴＴ（）０１＋１１０／Ｔ０１＋１

１０（）

Ｔ（４）其中λ０为标准状态下的空气分子平均自由程，λ０＝

０．０６６５μ

ｍ．２

）沉降在外力作用下，气溶胶粒子与气流之间产生相对运动，产生阻力．假定粒子所受阻力服从斯托克斯定律，则随着粒子和气流的速度差的增大，阻

力也随之增大．

若力场是均匀的，则粒子必有一个最终沉降速度．在重力场中，孤立粒子的最终沉降

速度即Ｓｔｏｋｅｓ速度为［３

］

ｖ＝Ｃｕ

ρｐｄ２

ｐ１８μ

ｇ（５）其中ρ

ｐ为胶粒的质量密度．实际情况中，分散系内的气溶胶粒子会互相影响，实际的沉降速度与

Ｓｔｏｋｅｓ速度会有一定差异．

在单分散体系的假设下，有沉降速度公式［３］

为

珔ｖ＝ｖ（１＋Ｓ ）

（６）其中， 为分散系中粒子的体积浓度；Ｓ为沉降系数，此处Ｓ＝－６．５５．３

）凝聚先考虑简单单分散热凝聚或Ｓｍｏｌｕｃｈｏｗｓｋｉ凝聚：它假设粒子是单分散性的；粒子一旦接触就互相黏附；

这些粒子生长得很慢．在这些假设下推导出气溶胶数量浓度的变化率［

３］

为ｄＮｄｔ

＝－ＫＮ

２

（７）其中，Ｎ为单分子粒子数量密度；Ｋ为凝聚系数，Ｋ＝２　ｋＢＴＣｕ

３μ

．

对变化率积分，有Ｎ（ｔ

）＝Ｎ０

１＋Ｎ０Ｋｔ

（８

） 中位数粒度［３］

为

ｄ（ｔ）＝ｄ０

Ｎ０Ｎｔ（）

１

３

（９） 当粒度为正态分布时，

由Ｌｅｅ和Ｃｈｅｎ提出的计算公式可以得到平均凝聚系数［

３］

珡Ｋ＝２　ｋＢＴ３｛

μ１＋ｅｘｐ　ｌｎ２σ（）ｇ＋２．４９λ（）

ＮＭＤ

ｅｘｐ　０．５ｌｎ２σ（）ｇ＋ｅｘｐ　２．５ｌｎ２

σ（）［］｝

ｇ（１０

）其中，ＮＭＤ为计数中位直径；σｇ为几何标准偏差，

它被定义为

σｇ＝

ｄｐ（

Ｆ＝８４．１％）ｄｐ（

Ｆ＝５０％）（１１

） 当粒度分布符合正态分布时

σｇ＝

ε＋σε

（１２） 由于这些粒子生长得很慢，

在考虑的时间范围内，σｇ的变化不计，

则粒子数为Ｎ（ｔ

）＝Ｎ０

１＋Ｎ０珡Ｋｔ

（１３

） 中位数粒度

ｄ（ｔ）＝ｄ０

Ｎ０Ｎ（ｔ

（）

）１

３

（１４

）１．２　简化模型

对模型中的影响因素进行适度的分析，排除一些影响很小的项，以简化模型，降低求解的难度．１

）忽略凝聚的影响考虑粒度为正态分布时的凝聚情况，考察一个如图１所示的具有代表性的实例，通过估算凝聚速度说明凝聚的效果在我们考虑的时间尺度下可以忽略．

图１　北京西北城区夏季大气ＰＭ１０的粒度分布［４］

仅考虑烟尘集合体，注意到它的频率分布呈

现单峰分布，用正态分布对其进行拟合，通过参数

估计可以得到ＮＭＤ＝０．２９５μｍ；σ

ｇ

＝１．８６．

在标准状态，由Ｌｅｅ和Ｃｈｅｎ给出了公式（１０）

的平均凝聚系数［３］为

珡Ｋ＝６．４１×１０－１６　ｍ－３·ｓ－１

要求得气溶胶的数量密度，胶粒的线度应该

用质量中位径ＭＭＤ，有ＮＭＤ和ＭＮＤ的换算关

系估算［３］

ＭＭＤ＝ＮＭＤ·ｅｘｐ（３ｌｎ２σｇ）＝０．９３７μｍ

考虑ＰＭ１０值为３００μｇ·ｍ－３时的情形，通过

计算此时Ｎ０≈１．５×１０７　ｍ－３，考察粒子数随时间

变化的方程

Ｎ（ｔ）＝

Ｎ０

１＋Ｎ０珡Ｋｔ

（１５）

将粒子数变为初值的１／２所用的时间称为半衰期Ｔ，则由以上讨论可以得到

Ｔ＝１

Ｎ０珡Ｋ≈

１０８ｓ≈３．３年

注意到沉降速度大约为１０－５　ｍ·ｓ－１，如果污染源离地地面的距离为１０１　ｍ量级，则污染源处霾粒子沉降到地面的时间约为１０６ｓ，可见在此情形下凝聚的时间尺度远远大于我们考虑的时间范围，故在以下的讨论中不考虑凝聚的作用．２）忽略多粒子互相作用对沉降的影响

考虑多粒子相互作用时的沉降公式：

珔ｖ＝ｖ（１＋Ｓ ）（１６）注意到ｖ和珔ｖ的差别仅仅在于Ｓ 这一项．仍然考虑文献［２］中的ＰＭ１０样本．ＭＭＤ＝０．９３７μｍ，Ｎ０≈１．５×１０７　ｍ－３，则可以计算出气溶胶体积密度 ≈１０－１１，由于Ｓ＝－６．５５，则Ｓ ≈１０－１０．可以看到，考虑多粒子相互作用的速度珔ｖ和Ｓｔｏｋｅｓ速度ｖ的偏差在１０－１０量级，所有认为斯托克斯速度即沉降速度是合理的．

１．３　模型建立

通过以上的分析，霾分散系的演化可以归结为扩散和Ｓｔｏｋｅｓ沉降的耦合作用．沉降的效果可以等效于对流扩散方程中的对流项．由于考虑的是无风环境，控制霾随时间演化的方程如下表示：

ｃ

ｔ

＋瓫·!ｃ＝Ｄ!２ｃ（１７）其中瓫为沉降速度．

２　模型求解

考虑前面分析过的例子，假设环境处于标准状态，ＰＭ１０值为３００μｇ·ｍ－３．代入一个基于高斯烟羽模型的较具有代表性的初值条件，求解霾分散系的演化情况．

２．１　参数的确定

由中位数粒度和求解假设可以确定扩散系数为

Ｄ＝

ｋＢＴ　Ｃｕ

３πμｄｐ

＝１．１７×１０－１０　ｍ２·ｓ－１

同样可以由公式确定沉降的Ｓｔｏｋｅｓ速度为ｖ＝Ｃｕρｐ

ｄ２ｐ

１８μ

ｇ＝１．２３×１０－５　ｍ·ｓ－１２．２　无源空间对流扩散方程的解析解

为了简化问题，假设考虑上方为无界空间，地面按一定比例吸收运动到地面上的污染物，假设我们已经得到无界边界条件下的浓度场ｕ，那么“反射”的效果即ｚ＜０的ｕ场“翻转”到ｚ＞０的空间中并乘上一个比例系数，再与ｚ＞０处的ｕ场叠加，就可以得到有地面边界条件下的浓度场ｃ控

制ｕ的方程如下为

ｕ

ｔ＝Ｄ ２　ｕ ｚ２＋ｖ ｕ （）

ｚ＋Ｄ１ｒ ｒｒ ｕ （）ｒｕｔ＝０＝Ｃ（ｒ，ｚ）

ｕｒ＝０＝有界，ｕｒ＝∞＝烅

烄烆

０（

１８

） 运用格林函数，

先考虑在ｒ＝０处一小团气体的质量为ｍ０的气体的扩散情况，在球坐标下

ｇ ｔ＝Ｄ!２　

ｇｇ

ｔ＝０＝

ｍ０４πδ（ｒ）

ｒ２

ｇｒ＝０

＝有界，ｇｒ＝∞＝烅烄烆

０

（１９

） 解得

ｇ（ｒ，ｔ）＝ｍ０４π２

ｒ∫

∞０

ｋ·ｓｉｎ（ｋｒ）·ｅｘｐ（－Ｄｋ２

ｔ）ｄｋ（２０

） 取ｍ０

４π

２＝１，

Ｄ＝１，ｔ＝１，２，３浓度随距中心距离ｒ的分布示意图如图２所示

：

图２　时间ｔ取１，２，３ｓ时的浓度分布示意图

加入对流项，相当于坐标系在以ｖ的速度向

－ｚ方向平移，

只需在解中加入一个坐标变换ｒ′＝ｒ２＋（ｖｔ）２

＋２ｒｖ　ｔ·ｃｏｓ槡θ

或

ｒ′＝ｒ－ｒ０－瓫

ｔ则仍有

ｇ（ｒ

′，ｔ）＝ｍ０４π２ｒ

′∫

∞０ｋ·ｓｉｎ（ｋｒ′）·ｅｘｐ（－Ｄｋ２

ｔ）ｄｋ 将由高斯烟云模型得到的初值条件看作ｄｅｌ－

ｔａ函数的叠加，

则整块烟云的浓度场应该是一小团气体扩散产生的浓度场ｇ的叠加．

ｕ（ｒ，ｔ

）＝

（Ｖ）

ｇ（ｒ－ｒ０

－瓫ｔ，ｔ）Ｃ（ｒ０

）ｄＶ其中ｄＶ为ｒ０附近的一个面积微元．

设地面的反射率为ｋ，则由之前提供的计算方法

ｃ（ｒ，ｔ）＝ｕ（ｒ，ｔ）＋ｋ·ｕ（－ｒ，

ｔ）其中ｒ＝（ｘ，ｙ，

ｚ），ｚ＞０．代入初值条件和之前求得的系数，可以得到任意时刻的浓度分布．

２．３　初值条件的确立

对于连续均匀排放的点源，源强为Ｑ（ｍｇ·ｓ－１），离地面的有效排放高度为ｈｅ（

ｍ）．令平均风速ｕ（ｍ·ｓ

－１

）沿ｘ轴方向．当风速较大，稳定之后，可以假设在ｙ、ｚ方向上浓度Ｃ呈正态分布，则扩散公式［

２，３］

为Ｃ（ｘ，ｙ，

ｚ）＝Ｑ２ｕπσｙ

σｚ·

ｅ

－ｙ

２２σ

２

ｙ

·ｅ

－

（ｚ－ｈｅ

）２２σ２ｚ

＋ｅ

－

（ｚ＋ｈｅ

）２２σ２［］

ｚ

（２１

）式中：

Ｃ—下风向任意位置（ｘ，ｙ，ｚ）的污染浓度，ｍｇ

·ｍ－３；σｙ—扩散系数，

ｙ方向的标准差，ｍ；σｚ—扩散系数，

ｚ方向的标准差，ｍ；Ｑ—排放源强，ｍｇ

·ｓ－１；Ｕ—排放口高度处的平均风速，ｍ·ｓ

－１

；ｈｅ—有效排放高度．考虑无风的情形，烟云有向上的初速度，忽略

反射项，

则可以将高斯烟云模型改写为如下形式Ｃ（ｘ，ｙ，

ｚ）＝Ｑ２ｕπσ

２

·

ｅｘｐ－ｒ２

２σ（）

２（ｚ＞ｈｅ）（２１

）其中：ｒ２＝ｘ２＋ｙ２

，

σ为ｒ方向的标准差．这是一个比较具有代表性的情形，若将它作

为初值条件，则它的物理意义为由一点源以较大的初速度向上排放污染物至稳定．之后停止排放，将此处的浓度分布作为初值条件．

为了求解的方便，对初值条件做如下规定：（１）为了避免在场源处出现奇点，在ｚ＞ｈｅ的空间中，对分布

Ｃ（ｘ，ｙ，ｚ）＝Ｑ２ｕπσ

２

·ｅｘｐ－ｒ２

２σ（）

２（ｚ＞ｈｅ） 令σ２＝１＋０．２＊（ｚ－ｈｅ）２

．

（２）为了让浓度场在空间内连续，在ｚ＜ｈｅ的空间中，假设

Ｃ（ｘ，ｙ，ｚ）＝Ｑ２ｕπ

·ｅｘｐ－ｒ２＋（ｚ－ｈｅ）２

（

）

２（ｚ＜ｈｅ）

（

３）为了计算方便，ｕ取１ｍ·ｓ－１

，Ｑ取２πｍｇ·ｓ－１，ｈｅ取１

０ｍ．由此得到初值条件的示意图图３所示．我们将这种形式的的烟云分布作为讨论问题

图３　基于高斯烟羽模型的初值条件

的初值条件．这只是一个值得参考的有一定参考

意义的情形，事实上２．２中的讨论可以处理各种形式的初值条件．

２．４　数值模拟

代入初值条件和之前求得的系数，可以得到任意时刻的浓度分布，代入系数Ｄ＝１．１７×１０－１０　ｍ２·ｓ－１，ｖ＝１．２３×１０－５　ｍ·ｓ－１，当反射率ｋ＝１，ｔ取６，９，１２，２４天，浓度分布如图４所示．可见，在无风的环境下，气溶胶污染物的变化是很缓慢的．此外，我们发现在我们考察的情景下，扩散的效果相较于沉降也非常小．这是因为霾气溶胶的粒度较大，沉降速度大，扩散系数小．霾的演化中起主要作用的是沉降作用

．

图４　当反射率ｋ＝１，第６，９，１２，２４天的浓度分布示意图

３　结语

本文通过分析，略去了凝聚对粒度的影响和

多粒子相互作用对沉降速度的影响，合理地简化

了模型．在考虑沉降和扩散的耦合作用时，将沉降

的效果用对流等效起来，取得了较好的效果．并通

过叠加的方法求出了对流项恒定条件下对流扩散

方程的解析解．找到一种典型的烟云分布，通过数

值方法解出其演化情况．

但是这一模型也有许多局限．它需要温度恒定，

无风或者至少风速恒定这样很强的假设．该模型对污

染物的沉积和化学转化过程只做了十分粗略的处理，

如果在此方面需要考虑时，模式的应用受到限制．

参　考　文　献

［１］　张小曳．中国大气气溶胶及其气候效应的研究［Ｊ］．地球科学

进展，２００７，２２（１）：１２－１６．

［２］　向晓东．气溶胶科学技术基础［Ｍ］．北京：中国环境科学出版

社，２０１２．

［３］　平措．大气污染扩散长期模型的应用研究［Ｄ］．硕士学位论

文．天津：天津大学环境科学与工程学院，２００６．

［４］　邵龙义，时宗波．北京西北城区与清洁对照点夏季大气

ＰＭ１０的微观特征及粒度分布［Ｊ］．环境科学，２００３，２４（５）：

１１－１６．

■