pm2.5成因、传播及控制

题 目 PM2.5 成因、传播及控制 摘 要 本文对 成因、传播及控制进行研究，找出形成因素，并通过建立高斯模型2.5PM 分析扩散情况，结合实际问题得出最佳控制方案。 问题一要求对 的形成因素进行定量描述。查找资料了解 的主要形成. 2.5PM 因素有 、 、 、 ，利用 软件结合相关性分析法对这几个因素进2SO10C2NOSP 行分析，得出各自之间的相关密切程度，由相关性公式： 22()xyr 可得出结论： 与 和 两变量之间是中度相关，其中与 相关度更2.5PM10 CO 高，相关值达到 ，而与 和 这两变量是呈弱相关。0.82S 问题二要求在仅考虑大气运动影响下，建立 扩散传播的数学模型。采用大.5PM 气污染的经典高斯扩散模型通过查阅相关资料，发现连续点源的平均气体流，其浓度 分布符合正态分布规律，可知污染物浓度在 、 轴上的分布符合正态分布。取烟云轴yz 线为 方向（平均风向） ，得出无界情况下空间某一点的浓度函数： 。x (,)xyzH 问题三要求在问题一和问题二的基础上，综合最低成本和最佳效果两因素，提出 最佳控制方案。考虑到这几个因素之间会存在某些联系，选择对几个因素建立多目标 系统动力学模型，利用 软件进行求解，求得减少量和投入费用，并通过比对得Matlb 知将综合治理与专项治理相结合能有效地解决实际问题。 关键词 相关性分析 高斯扩散模型 多目标系统动力学模型2.5PM 1 一、问题背景与重述 1.1问题背景 （细颗粒物）这个大家原本很陌生的专有名词，因为 2011 年 10 月以来我2.5PM 国多地灰霾天气造成严重大气污染，而迅速成为社会热词。之后， 被纳入新的2.5PM 环境空气质量标准，并在一些地方开始监测。2012 年全国“两会”上， 首次写 入政府工作报告。全社会对 的关注和重视，折射出当前我国环境污染的严峻性。2.5PM 环境是人的生存之本、发展之基。让百姓喝上干净的水，呼吸清洁的空气，在良好的 环境中生产生活，是关乎社会发展和人民福祉的大事。现阶段，环境问题不仅是经济 问题、发展问题，而且是政治问题、民生问题等。 是指大气中直径小于或等于 微米的颗粒物，也称为可入肺颗粒物。它的2.5PM. 直径还不到人的头发丝粗细的 。虽然 只是地球大气成分中含量很少的组分，1/202.5P 但它对空气质量和能见度等有重要的影响。与较粗的大气颗粒物相比， 粒径小，2.5PM 富含大量的有毒、有害物质且在大气中的停留时间长、输送距离远，因而对人体健康 和大气环境质量的影响更大。 年 月，国务院同意发布新修订的环境空气质量 标准增加了 监测指标。2.5 1.2问题重述 为更加有效的治理 ，试建立数学模型，分析我国的 环境问题的现状：.PM2.5P 1、通过查阅资料，分析出 的形成因素，并建立数学模型进行定量描述；25 2、在仅考虑大气运动影响下，建立 扩散传播的数学模型；. 3、通过 的形因模型及传播模型，综合防治成本和治理效果两方面因素，提出最.5 优的控制方案。 二、问题分析 2.1问题一的分析 要对 的形成因素建立数学模型进行定量描述，通过查阅网上资料，了解到2.5PM 的形成因素主要有 、 、 、 。考虑到要分析我国现状，数据太. 2SO10PMC2NO 大，所以可以以徐州为例。 采用 简单相关关系方法对选取的数据进行分析，分析 与这几个因earson 2.5PM 素之间的相关程度。可以利用 对各变量进行相关性分析和显著性检验，得出相关 性矩阵，从而得出 与其形成因素之间的相关密切程度。然而考虑到简单相关系.5 数法可能不能真实反映出两变量之间的相关性，通常可能受到不止一个变量的影响， 选择再对其进行偏相关性分析，有效地分析出 、 、 、 、 几个2.5P2SO10C2NO 变量间净相关的强弱程度。 2.2问题二的分析 要在仅考虑大气运动条件下，分析 扩散的情况。为此选取发电厂作为研究.M 点，引入风速变量对 浓度的影响，选取高斯扩散模型进行分析。在不考虑垂直2.5P 风速，假设空间 浓度服从高斯分布的前提下，运用高斯模型合理的计算出发电M 厂周边地区的污染物浓度。 2.3问题三的分析 出现和形成与国民经济的发展和国家工业化进程有着密不可分的关系，基2.5PM 于问题一和问题二的研究，要提出最优控制方案。考虑到要成本最低，即对 治2.5PM 理费用最低，可以建立多目标优化模型，再利用软件求出结果。除此之外还需考虑治 理效果最优，则可以对 分别进行综合治理和专项治理。2.5P 对于综合治理与专项治理既可以分别研究，也可以结合两者进行研究，比较两者 方式的治理时间与收益。最优控制方案则可以结合综合治理和专项治理，专项治理加 速治理，综合治理巩固治理结果。 三、模型假设 结合本题的实际，为了确保模型求解的准确性与合理性，我们排除一些干扰因素， 提出以下几点假设： 1、假设所使用的数据能客观地反映现实情况，且数据准确真实； 2、假设污染物在输送扩散中质量守恒，且污染源的源强均匀，连续； 3、假设环境为无边界空间，污染物在扩散过程中无遮挡物； 4、假设风的平均流场稳定，风速均匀，风向平直； 5、假设污染物的浓度在 轴上的变化分布符合正态分布；yz、 6、假设污染物扩散服从扩散定律，即单位距离通过单位法向面积的流量与它的浓度梯 度成正比； 7、在探究 的形因过程中忽略工地扬尘的作用；2.5PM 8、扩散过程中不考虑云团内部温度的变化，忽略热传递、热对流及热辐射。 四、符号说明 为了便于问题的求解，我们给出以下符号说明： H 离地面的有效排放高度i 扩散参数u 平均风速 坐标轴grad 梯度 空间域 相关系数Q 的量2.5PMq 单位时间排放物div 散度记号Cpm 浓度.,xyz 空间任意一点的污染物浓度0 污染源排放的污染物总量 3 五、模型的建立与求解 5.1问题一的模型建立与求解 5.1.1模型建立 工业作业、工地扬尘和燃煤发电等活动产生大量的污染物，致使 成为徐州2.5PM 污染物的重要成员之一，灰霾现象频繁发生，时刻危害着人体健康。通过上网查阅资 料，得知主要因素有： 、 、 、 。考虑到要分析我国情况，数据太2SO10PMC2NO 大，那么现在就徐州为例，选取 2015 年 1 月至 3 月的数据进行分析，由于各变量之间 没有顺序之分，那么就可以采用二元变量相关分析法中的 简单相关系数方法，earson 分析其与这几个因素之间的相关密度程度。 相关系数的计算公式如下： Pearson22xyr 5.1.2模型求解 利用 对各变量进行相关性分析和显著性检验，得到如下表：SP 表 1 简单相关分析表2.5PM2SO10PMCO2N 相关性earson1.00 0.51 0.86 0.87 0.712.5M 显著性（双侧） 0.00 0.00 0.00 0.00 0.00 相关性P0.51 1.00 0.56 0.64 0.63 2SO显著性（双侧） 0.00 0.00 0.00 0.00 0.00 相关性earson0.86 0.56 1.00 0.82 0.7410 显著性（双侧） 0.00 0.00 0.00 0.00 0.00 相关性P0.87 0.64 0.82 1.00 0.80 C显著性（双侧） 0.00 0.00 0.00 0.00 0.00 相关性earson0.71 0.63 0.74 0.80 1002NO 显著性（双侧） 0.00 0.00 0.00 0.00 0.00 由上表得到各指标的相关性矩阵和 与各因素之间的相关性矩阵，分别如下2.5PM 所示： 2.510.54.8620.7318. 4.6=014,.86273,0.18 862.17301.07.81PMr r 根据两矩阵中数据显示，这几个指标两两之间相关密度程度都较高，且相关值都 为正数，由此可知各指标之间都是成正相关关系。通过分析 与各因素之间的相2.5PM 关性矩阵，可以发现它与 和 高度相关，与 和 成中度相关。10PMCO2SNO 然而简单相关系数法可能不能真实反映出两变量之间的相关性，通常可能受到不 止一个变量的影响，那么现在就需要利用相关性分析的另一方法即偏相关系数分析法， 有效地分析出 、 、 、 、 几个变量间净相关的强弱程度。2.52S2 利用样本数据进行偏相关性分析，在分析变量 和变量 之间的净相关时，当控1x2 制了变量 的线性作用后， 和 之间的一阶偏相关系数为：3x1x212332r 同样利用 软件对这几变量进行偏相关性分析，得到下表：SP 表 2 偏相关分析表2O10PMCO2N 偏相关性值 0.533 0.562 -0.031 -0.189 部分相关性值 0.256 0.277 -0.013 -0.079 置信区间 (0.210,0.432) (42.024,81.424) (-0.624,0.466) (-0.469,0.028) 通过分析上表， 和 偏相关性值在区间 之间，且 的置信区间10PC0.5,8（ ） 最大，而和偏相关性值为负值，得到下面结论： 与 和 两变量之间是中2P1C 度相关，其中与 相关度更高，相关值达到 ，而与 和 这两变量是呈弱O.2NOS 相关。 5.2问题二的模型建立与求解 5.2.1模型建立 实际生活中汽车尾气、工业废气、燃煤发电和秸秆燃烧等过程中都会燃烧排放大 量 ，再由气象、地形等原因对外扩散，导致出现令人惊愕的城市灰霾天气现象，2.5PM 对于它的扩散问题展开研究。由于排放源不一，故现在就以工厂为例，以发电厂烟囱 正下方的地面为坐标原点 ，平均风向为 轴，指向下风向，铅直方向为 轴，0,（ ） xz 轴则水平垂直于风向轴，建立空间坐标系如图 1，则 排放点 距有效地面的y 2.5PMO 高度为 ，则排放点位置坐标为 。H,OH 图 1 空间系示意图 记时刻为 时，空间任意一点 的浓度为 ，根据假设可知单位时间t 2.5PM,Cxyzt 通过单位法向面积的流量与浓度梯度呈正比，则有： 5 iqgradC 是扩散系数， 表示梯度，负数表示由浓度高向浓度低的地方扩散。,ixyzgrad 考察空间域 ,其体积为 ，包围的曲面为 ， 为一规则的球面，则在 内通V1s2 ,t 过 的流量为： 11tQqndts 空间域 内 的增量为：2.5PM,VCxyztCxyztV 从污染源释放 的总量为：.00tQdtp 5.2.2模型求解 根据“质量守恒定律”和“气体泄漏连续性原理” ，单位时间内通过所选曲面 的s 向外扩散的 与 曲面内 质量的增量之和，等于排放源在单位时间内向外2.5PMs2.5P 排放的颗粒物。有： 012+Q 即： 10,t tVdtCxyztCxyztdVqndssp 又根据曲面积分的 公式：Gaus sVqndivq 则有： 0 0,limlimtt tkdivgraCtCxyztCxyzt 由上式得： 0VVVdivqt P 由以上公式并利用积分中值定理得： 222,i xyzCCCdvgar txyzt 上式是无界区域的抛物线型偏微分方程，根据假设，初始条件为作用在坐标原点 的点源函数，记作： ,0,zQxz 考虑到在做题时忽略了地面的作用，然而在点源的实际扩散中，污染物可能受到 地面污染物的阻挡，因此现在考虑地面对扩散的影响，假定污染物在扩散过程中，质 量保持不变，到达地面时不发生沉降或化学反应而全部反射。这可认为地面就像镜子 一样，对放射性气体起全放射作用，可用“像源法”处理，如下图建立三个坐标系， 一是以泄漏源（实源）为坐标原点，二是以泄漏源在地面的投影点为原点， 点是空p 间的任意一点，坐标为 ，三是以泄漏源。关于地面的像对称源（像源）为原点,xyz 。把点放射性气体浓度看成两部分（实源和像源）作用之和，由此可推导出： 图2 空间坐标系 其中，实源的贡献率为： 22,2yzyzQCxyzHexpexkH 虚源的贡献率为： 22,2yzyzxyzexpexk 由此可推导出： 222,2yzyz zQCxyzHexpexexpkH 式中 是污染物在 方向分布的标准偏差，是距离 的函数， 是污染物在 方y yz 向分布的标准偏差，是距离 的函数，单位为 。m 已知高斯扩散模型的特定条件是平均风速一定，那么在现实生活中不可能做到这 样的情况，那么现在讨论风速的问题。 （1）无风点源的解 此时 ，根据前面的分析结果可知：0u222,0,i xyzCCCdvgar txyzt 对于瞬时点源，应满足如下初始条件和边界条件： 初始条件： 0Hx时 ， 处 ，处 ， 边界条件： C时 ， 得到解为： 221,2yzxyzqxyzepu， 为了简化模型，令 、 与 的关系为：x12yzx 7 式中 分别为常系数。1212， ， ， （2）有风瞬时点源的解 假定有一定常的平均风速，在坐标原点时释放污染物，该污染物会随风飘动，并 因扩散而不断胀大。这个污染物的中心相对于固定坐标系的位置为 ，如再取,xuto 一个移动坐标系，原点设在移动污染物的中心，则利用坐标移动得到有风瞬时的扩散 公式： 21!, +22xyzzyzutq nCxyzHexpur， (3)有风连续点源的解 由于污染源的持续作用，乐意认定扩散过程为定常态，即，浓度仅是空间坐标的 函数。而且在风速大于高斯扩散模型所假定的均速时认为方向上风的平流输送作用引 起的物质质量大于扩散作用，这时忽略 。2icx 根据题意和假设可知 源强，即单位时间内排放的污染物， ，由此可得：q /gs21,yzyzqCxyzepu 5.3问题三的模型建立与求解 5.3.1 的综合治理2.5PM 治理问题是一个全方位的、多元化的问题，涉及到很多方面，它的出现和 形成与国民经济的发展和国家工业化进程有着密不可分的关系。工业燃煤、工地扬尘、 火力发电等作业都会造成 量和浓度的增加。在问题一和二的基础上，可以了解2.5P 到它的几个影响因素，然而这几个因素之间还存在着紧密联系，因此这里可以建立多 目标模型探究合理的治理方案。 建立多目标模型： minCpgd2 105%.60.8ijijiCpmMstQPp 其中， 表示投入治理资金，单位为百万元， 表示第 个项目经过 年后总投Mijij 资数， 表示 年后 浓度的降低值。Cpmi2.5P 利用 软件计算得到五年治理投资的费用如下表所示：atlb 表3 五年治理投资费用 年份 01212013204 投入费用 8757568综合治理 减少量2.5P11 8 36 24 19 投入费用 38 33 24 18 16专项治理 减少量.M32 27 25 17 14 根据求解结果和过程， 综合治理主要集中在环境保护和基础设施建设上。2.5PM 绿化环境是环境改建的重要部分，通过对外绿化以增加对扬尘的吸收，同时对绿化带 进行定期清理，以做到最大化吸收，这部分占投资费用的 。除此之外还需扩大交30% 通网络，发展公共交通，推行清洁能源的使用，这部分所占的比例是 。对于基础45 设施建设这方面就必须从技术上入手，增加科研投入与创新，改进技术，推广煤炭清 洁化利用技术，优化能源结构，减少化工燃料的使用。同时本着优胜劣汰的原则，淘 汰一些在城市上风向且对大气污染较重的企业，努力控制源头这部分约占 。2 5.3.2 的专项治理2.5PM 在治理单位费用已知的情况下，对徐州市 进行专项治理。通常在专项治理2.5PM 过程中要求成本最低，效果最显著，所以在各个变量之间的相互关系的前提下，同样 地找出约束条件，建立多目标系统动力学模型，找出最核心的影响因素，当地的实际 情况来给管理部门相应的规划安排。 建立多目标系统动力学模型 minaxCp2105%.60.8ijijiCpmMstQPp 在专项治理下，徐州市治理费用和 浓度减少量如下表：2.5 表 4 专项治理各年投入产出表 年份 2010110132014 投入费用(百万) 68 63 60 55 59 减少浓度2.5PM48 35 28 34 18 由表4中数据分析可知专项治理投入费用一定的情况下， 减少量逐渐减少，2.5PM 这就导致了对该问题有监管，有管理，却不尽人意的局面，因此徐州市在进行专项治 理的同时也要配合着进行综合治理。 因此接着分别考虑专项治理、综合治理以及综合治理与专项治理相结合情况下 治理时间与收益的对比图，如下图所示：2.5P 图 3 不同治理模式收益对比图 9 由上图所示，仅做到综合治理不易于较危机的情况，仅做到专项治理，容易出现 “有监管有治理”的局面，然而同时进行综合治理和专项治理，更有效地减少污染物 的排放，两者相互约束，大大提高了工作效率，获得预期的效果。 六、模型检验 问题一采用的是简单相关和偏相关分析法，利用 软件对结果进行显著性检验。SP 假设总体相关系数为 ，则有： ，0:H1: 对简单相关系数的统计检验是 统计量：t2-rn 如果 ,则拒绝 。显著差异；否则接受 ，无显著差异。所以得到各变量的nt0 0 显著性检验值均为 ，此时为极显著相关。 同样对偏相关分析法进行显著性检验，一般假设显著性水平为 ，通过带值可0.5 以知道 值小于 ，显著性较好。p.5 七、模型的评价与改进 7.1模型的评价 7.1.1模型的优点 1、相关性分析较为准确地分析出几个因素之间的相关密切程度，比较简明直接； 2、利用相关性分析，既减少了需要处理的数据量，又保证了 监测点的代表性与2.5PM 科学性，从而大大简化了问题； 3、利用高斯扩散模型将实际问题转化为抛物线偏微分方程的问题，较为理性地分析出 扩散时的情况；.5PM 4、分析过程中排出了一些不必要的麻烦，使得研究更加简单； 、利用多目标优化模型，有条理地约束分析。 7.1.2模型的缺点 1、题目要求分析全国情况，而本文只选取徐州为观察点，对 1 至 3 月的数据进行分析， 只能说明一部分情况，不能准确说明实际情况； 2、在研究扩散问题时没有考虑工地扬尘的作用和下垫面、地形、湿度等因素的影响， 不能做到准确分析； 3、对于距离较远的扩散检测不到，容易被其他污染源干扰。 7.2模型的改进 1、相关性分析法较为准确地分析出各自之间的相关密切程度，然而在这里结果存在模 糊性，可以利用主成分分析找出最主要因素。 2、题目要求分析我国 PM2.5 的情况，然而本文是以徐州为例，也只选择了 1 至 3 月的 数据，存在一定的局限性，所以在继续探究时应选取更多的数据，更多的城市，更 贴切地解决问题。 3、问题二中选用的高斯扩散模型是在风速一定的基础上进行探究，忽略了热对流、热 传递和热辐射等情况的影响，所以为了进一步探究扩散情况，需要把这些因素考虑 在内。 八、模型推广 相关性分