污染物扩散模型深圳数学建模

1、赛区评阅编号（由赛区组委会填写）：2015高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了全国大学生数学建模竞赛章程和全国大学生数学建模竞赛参赛规则（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站 下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、 网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果 或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正 文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛

2、的公正、公平性。如 有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号（从 A/B/C/D中选择一项填写）： C我们的报名参赛队号（12 位数字全国统一编号）： 参赛学校（完整的学校全称，不含院系名）：参赛队员（打印并签名）：1.章成俊2.杨超3谢锦指导教师或指导教师组负责人（打印并签名）：日期：年月日2015高教社杯全国大学生数学建模竞赛编号专用页赛区评阅记录（可供赛区评阅时使用）：评阅人nn备注送全国评阅统一编号（由赛区组委会填写

3、）：全国评阅随机编号（由全国组委会填写）：对垃圾处理厂污染的动态监控及居民补偿摘要城市垃圾处理问题是一个世界性难题。目前垃圾焚烧正逐步成为中国垃圾处理的 主要手段之一。 本论文构根据题目设置的垃圾处理厂规模， 建立了环境动态监控体系， 并根据潜在污染风险对周围居民进行了合理经济补偿的设计。对于问题（1），为了实现对垃圾焚烧厂烟气排放及相关环境影响状况的动态监控， 本论文在高斯烟羽模型的基础上进行改进，引入温度、降雨对污染物扩散的影响，建 立了新的污染物扩散模型。本论文创新性的提出了风雨影响指数M ,用来衡量风向、降雨对颗粒物扩散的影响。 本论文将抽象的污染物含量形象化， 利用空气污染指数 AP

4、I 描述具体的污染程度及其给周围居民带来的影响。并且从不同角度给出了模型检验， 验证了所建模型的准确性。对于问题（ 1）具体赔偿方案的制定，在综合考虑了不同方位风向频率、受污染时 间、受污染程度的基础上，本论文使用了层次分析法，并且进行了一致性检验，使得 赔偿方案具有说服力。通过 MATLAB 编程，计算出当政府和垃圾处理厂共支付风险 赔偿金为 N 时，得出居住地的每位居民应得的赔偿金额计算公式。 对于监测点的设置， 经计算共需 21 个，具体布置情况见后文。对于问题（ 2），在题目所述的发生事故的情况下，对污染物的具体含量进行了合 理的预测与假设。模拟出酸性物质与颗粒物的影响范围，并根据具体

5、的污染程度设置 不同的污染区。对每个污染区的不同情况设置更改监测点的设置，并且在问题（1）的基础上对居民的经济补偿进行合理修改。关键词： 高斯烟羽模型，层次分析法，空气污染指数，烟气抬升公式一、 问题重述“垃圾围城”是世界性难题，在今天的中国显得尤为突出。数据显示，目前全国 三分之二以上的城市面临“垃圾围城”问题，垃圾堆放累计侵占土地75 万亩。因此，垃圾焚烧正逐步成为中国垃圾处理的主要手段之一。然而，由于政府监管不力、投资 者目光短浅等多方面的原因，致使前些年各地建设的垃圾焚烧电厂在运营中出现了环 境污染问题，给垃圾焚烧技术在我国的推广造成了很大阻力，许多城市的新建垃圾焚 烧厂选址都出现因居

6、民反对而难以落地的局面。在垃圾焚烧厂运行监管方面，目前主 要是在垃圾焚烧厂内进行测量监控，缺少从周边环境视角出发的外围动态监控，因而 难以形成为民众所信服的全方位垃圾焚烧厂环境监控体系。深圳市某地点计划建立一个中型的垃圾焚烧厂，计划处理垃圾量1950 吨/天（设置三台可处理垃圾 650 吨/天的焚烧炉，排烟口高度 80 米，每天 24 小时运转）。从构 建环境动态监控体系、并根据潜在污染风险对周围居民进行合理经济补偿的需求出发， 有关部门希望能综合考虑垃圾焚烧厂对周围带来环境污染以及其他危害的多种因素 （例如，焚烧炉的污染物排放量、居住点离开垃圾焚烧厂的距离、风力和风向及降雨 等气象条件、地形

7、地貌以及建筑物的遮挡程度等等） ，在进行科学定量分析的基础上，确立一套可行的垃圾焚烧厂环境影响动态监控评估方法，并针对潜在环境风险制定出 合理的经济补偿方案。请你在收集相关资料的基础上考虑以下冋题：(1) 假定焚烧炉的排放符合国家新的污染物排放标准(参见附件1),根据垃圾焚烧厂周边环境设计一种环境指标监测方法，实现对垃圾焚烧厂烟气排放及相关环境影 响状况的动态监控。以你设计的环境动态监控体系实际监控结果为依据，设计合理的 周围居民风险承担经济补偿方案。(2) 由于各种因素焚烧炉的除尘装置(如袋式除尘器)损坏或出现其他故障导致污染物的排放增加，致使相关各项指标将严重超标(如：烟尘浓度、二氧化硫、

8、氮氧 化物、一氧化碳、二恶英类及重金属等排放超标，附件2给出了一台可处理垃圾350吨/天的焚烧炉正常运作时的在线排放监测记录)。请在考虑故障发生概率的情况下修 正你设计的监测方法和补偿方案。模型假设3#(1) .假设污染物在所建的坐标轴的y、z风向上分布为正态分布;(2) .假设全部高度风速均匀稳定；(3) .假设污染物在扩散中不会相互转化；(4) .假设源强是连续均匀稳定的；(5) .假设附件4所给风速距地面10米处所测得风速；三、符号约定C:任意点的污染物浓度，单位：mg/m3;Q:源强，单位时间内污染物排放量，单位：mg/s ；二y :侧向扩散系数；-z：竖向扩散系数；U :排放口的平均

9、风速，单位：m/s；H :烟囱的有效高度；x:污染源排放点至下风向上任一点的距离，单位：m；m；y :污染物的中心轴在直角水平方向上到任一点的距离，单位：z:从地表到任一点的高度，单位：m；Hs：烟囱的几何高度，单位：m；=H :污染物的抬升高度，单位：m;P,:大气压力，单位： KPa ；Qh：烟气热释放率，单位：KJ/sQv:实际排烟率，单位：m3/s；Vs：排气筒出口处污染物排出速度，单位：m/s ；Ts：烟气出口温度，单位：K;T :烟气出口温度与环境温度差，单位：K;。Pj(i =4,5.2,3; j =1,.,8):表示第i月的风向为j( j代表风向分别从正北方向沿顺时 针依次到到

10、西北方向)M :风雨影响指数；四、模型建立与求解4.1问题一: 4.1.1问题一分析根据垃圾焚烧厂周边环境设计一种环境指标监测方法，以实现对垃圾焚烧厂烟气 排放及相关环境影响状况的动态监控，这首先需要调查出各类污染源对这片地区的污 染影响。而在高斯烟羽扩散模型上进行改进，即可模拟出污染源对周边空气影响的大 致范围及污染程度。而后，根据具体的污染程度及污染物类型，选取合适的位置设置 监控点，并制定相应的周围居民风险承担经济补偿方案。4图1.对于问题一的分析流程图4.12模型建立（1）.高斯烟羽模型基本形式的数学表达式为:C(x,y, z,H)4z_H)2e出4z h)2erH 二 Hs H公式（

11、1）公式（2）其中C为任意点的污染物浓度，Q为单位时间内污染物排放量,U为排放口的平均风速，H为烟囱的有效高度，Hs为烟囱的几何高度， H为污染物的抬升高度，匚y与二z参数污染物在排放出来时，具有初始动量，且受空气浮力，所以污染物会有抬升高度。 根据烟气抬升公式，烟气抬升高度 厶H为：H -noQh1Hn2U公式（3）Q0.35PaQ-T公式（4）Qh T700公式（5）H12 1.5VsD 0.01QhU0.048 Qh -1700U公式（6）公式（7）2(1.5VsD +0.01Qh )lH =U以上公式中，Qh为烟气热释放率，Qv为实际排烟率，Vs为排气筒出口处污染物排出速度，Ts为烟气

12、出口温度，汀为烟气出口温度与环境温度差，no、n1、分别 为参数，具体值见下表：心、n、n2的选取Qh（单位：KJ/s）地表状况（平 原）nnn2Qh 2100农村或城市远 郊区1.4271/32/3城市及近郊区1.3031/32/32100 兰 Qh 35K农革或城市远 郊区0.3323/52/5城市及近郊区0.2923/52/5表1.n。、n1、rb取值参考一般已知地面10m高处的风速U10,而烟囱出口处的风速 U应取该处实测值。无 实测值时，可用以下幕指数法求出：公式（8）式中U、U10分别为烟囱高度H处和地面10m高度处10mir平均风速，m/s;当 H240m时，按240m计算。P为

13、风速高度指数，依赖于大气稳定度和地面粗糙度。应 根据观测结果，利用统计学方法求出。根据具体的观测数据，也可采用风速随高度变化的对数律或其它半经验分式。无实测数据时，P值可按下表选取。由附件5推测,此垃圾焚烧厂在深圳市的乡村。地区各稳定等级下的P值ABCDE或F城市0.10.150.20.250.3乡村0.070.070.10.150.25表2.各稳定度等级下的P值在本文中，根据广东沿海地区大气稳定度的分类方法探讨1 丄估算出深圳市大气稳定度等级频率，具体数值如下表:稳定度等 级ABCDEF频率1.833.8624.2748.6315.245.55表3.深圳市大气稳定等级频率对于参数二y与二z的

14、取值参考下表得出:大气稳定度66A0.22x心+0.0001 x)0.50.2 xB0 50.16x心+0.0001 x)0.12xC0.11x/(1+0.0001 x)0.50.08x/(1+0.0002x)0.5D0 50.08x心+0.0001 x)0 50.06x心+0.0015x)E0.06x/(1+0.0001 x)0.50.03x/(1+0.0003x)F0 50.04x/(1+0.0001x)0.016x/(1+0.0003x)表4.参数二y与二z的取值当污染热释放率Qh大于或等于2100KJ/S且污染物温度与环境温度的差值T大于或等于35K时，污染物抬升高度-H采用公式（3）

15、、（ 4）计算；当1700KJ/sQh 2100KJ/s时，污染物抬升高度 的 采用公式（5）、（6）计算；当Qh 700KJ/S或心TC35K时, 污染物抬升高度丄H采用公式（7）计算。此模型具体的坐标系建立如下图：图2.高斯烟羽模型坐标系的建立示意图利用高斯烟羽扩散模型可以表示出下风向任意一点处污染物浓度。因此，在题目 中给定的八个方向分别建立坐标轴，然后分片区计算污染程度。即八个方向因焚烧炉 的排污在324天内分别所受到的污染。因为空气污染物对居民生活造成的影响并不只是简单的线性关系，仅仅依据污染 物浓度去进行周围居民风险承担经济补偿方案设计较片面。所以本文还引入了另一评 价指标一一空气

16、污染指数（API）。API是由国家规定的定量描述空气质量状况的无量纲 指数，用以表征空气污染程度。API越高，则相应区域空气污染越严重。本文根据具体的S、Nx、颗粒物浓度计算出空气污染指数（API）0然后根据 高斯烟羽模型所计算出的污染物浓度及八个方位具体的 API值确定合理的周围居民风 险承担经济补偿方案。根据不同的空气污染指数等级利用层次分析法得出相应权重，再结合不同风向的 频率确定最终的周围居民风险承担经济补偿方案。4.1.3模型求解（1） 在MATLAB中建立高斯烟羽模型，得出八个方向的污染物浓度分布图如下:气体扩散尿度平面團50：0400030002000WOOWOO2D0D3000

17、4000enoo西冻（m-r 69#图3.空气中SQ的浓度#5000400030002000WOO-1000-2000-3000-400050004000z；o:o2000WOO-1000-2000-3000-4000气体扩散浓度平面團2016141210o西冻(m)图4.enoo空气中NOx的浓度气体扩散尿度平面團0西冻(m)enoo211图5.空气中颗粒物的浓度12#由以上图片观察，正东、东北和正北方向污染物浓度较高 八个方位的空气污染指数（API）分布如下图：5000 L1m11150000-4000-200002000400060Q0西-东图6.各个方位的API数值参照国家的对空气污染

18、指数（API）划分的等级，本文对八个方向也划分污染等 级。如下图：#13SOTojj-3000-2TO-1000100UOTO30TO300重污染健康人运动耐受力 降低，有明显强烈症 状，提前出现某些疾 病老年人和病人应当 留在室内，避免体力 消耗，一般人群应避 免户外活动表5.API指数及其对应的空气质量和人类生活的影响利用层次分析法，计算不同空气污染等级对用的权重。设定API数值越高及其对1234应的赔偿重要性越重要。得出正互反矩阵:1/21231/31/212141/31/21一通过MATLAB计算得出相应的赔偿权重:如下表:API3002003001002000100赔偿权重a0.46

19、730.27720.16010.0954表6.不同API等级对应的赔偿权重经计算，CI=0.0103 CR=0.0115 RI=0.9Q CIvRI且 CR3002003001002000100级别 一级区二级区三级区四级区赔偿权重口0.46730.27720.16010.0954表6区域内烟尘浓度严重超标，根据严重程度进行赔偿危险程度65544331级别丁一级区二级区三级区四级区赔偿权重0.46730.27720.16010.0954表7烟尘泄漏X轴浓度分布200015001000500亘TT0-500-10001500-2000 _ 2 050010001600200025003000西-东图13