空气污染数学建模.pdfVIP

A.污染气体的传播扩散 摘要 钢铁生产排放的污染气体是造成雾霾的重要原因之一，研究污染气体 的扩散特征，正确模拟污染气体的扩散过程，能够为钢铁生产集团提出更好 的治理管理措施，具有实际意义。 针对问题一：污染气体的排放速度为 300m/s,在不考虑风向风速及高度 影响的情况下，此问题即为二维平面的连续点源扩散问题，由此在二维 xoy 平面上建立连续点源扩散方程模型ut   (uxx  uyy )  f (x , y ,t) ，其中 为气体 扩散系数，本文中取为常数 10，f(x,y,t ) 为污染气体的排放速度，在本文中 恒为 300m/s ；对上述偏微分方程模型，本文采用 ADI 法（Alternating direction implicit，交替方向隐式法）求解出迭代格式，利用 MATLAB 编程， 求出模型一的数值解，并得到任意时刻污染气体的浓度分布情况。通过 SPSS 软件，对附件一所 的原始实际数据与模型一求解得到的模拟值进行 显著性检验，检验结果显示该模型与实际情况吻合。 针对问题二：考虑风向风速对污染气体扩散过程的影响时，在基于对问 题一求解的基础上，在模型一的扩散方程模型中加入风向风速的平流项，由 此得到有风情况下的模型 u   u   u   (u  u )  f (x ,y ,t) ，其中 t 1 x 2 y xx yy  ， 分别为风速在 x, y 方向的分量；对此模型同样采用 ADI 法求出迭代格 1 2 1 式，利用 MATLAB 编程，求出模型二的数值解，并得到任意时刻污染气体 的浓度分布情况。通过 SPSS 软件，对附件二所 的原始实际数据与模型二 求解得到的模拟值进行显著性检验，检验结果显示该模型与实际情况吻合。 针对问题三：考虑有风时增加高度的影响，此 问题即为三维空间的污染 气体扩散问题，考虑到三维模型的编程复杂度，而且污染气体的扩散在 xoy 平面上各向同性，可以将污染气体在y 方向的扩散等价为在 x 方向上的扩散， 此时便只需要建立 xoz 平面上的扩散模型。在基于对问题二求解的基础上， 在模型二的扩散方程中增加高度项，由此得到模型三为 u  2 (z )u  2u  u  f x , z ,t  ，其中 为 z 方向的扩散系数；对该扩散方 t 1 x xx zz 程同样采用 ADI 法求出迭代格式，利用 MATLAB 编程，求出模型二的数值 解，并得到任意时刻污染气体的浓度分布情况。 关键词 ：污染气体 扩散方程 ADI 法 数值解 一、 问题重述 目前，治理雾霾是人们最为关心的热点问题之一。中国社科院发布的《气 候变化绿皮书》中提及，雾霾形成的原因里，重工业、车辆尾气、土方施工都 榜上有名，其中钢铁生产也是造成雾霾的重要原因之一。 2 某钢铁生产集团烟囱污染气体的排放对周边地区大气污染的影响非常 大，为了提出更好的治理管理措施，需要对其污染气体扩散的特征进行分析。 现在，我们需要在三种情况下考虑污染气体的扩散过程： 1. 在不考虑风向和高度影响的情况下，建立模型，模拟某钢铁生产 集团的烟囱排放污染气体的扩散过程，假设烟囱的排放速度为 300m/s。 2. 考虑风向为东北风，平均风速 0.6m/s 的情况下，模拟污染气体的 传播扩散过程。 3. 在考虑风向的基础上增加高度的影响，建立模型，模拟污染气体 的传播扩散过程。 4. 基于上述模型结论，给该钢铁生产集团提供一个污染气体治理建 议报告。 二、 问题