空气污染数学建模.docx

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