《中型垃圾焚烧厂对于周边居民的潜在风险进行分析建模98953664》由会员分享，可在线阅读，更多相关《中型垃圾焚烧厂对于周边居民的潜在风险进行分析建模98953664（36页珍藏版）》请在金锄头文库上搜索。
1、中型垃圾焚烧厂对于周边居民的潜在风险进行分析摘要在国内垃圾焚烧厂出现各种环境污染问题的背景下，本组针对深圳市计划建立的中型垃圾焚烧厂对于周边居民的潜在风险进行分析，确定一种可行的环境动态监控体系，并设计了合理的周围居民风险承担经济补偿方案。对于本问题的研究，我们建立高斯大气扩散模型和概率模型进行总体分析。对于问题1，为了解决垃圾焚烧厂附近污染物浓度随空间位置变化的规律，我们建立了污染物扩散问题中最为经典的高斯大气扩散模型。首先，利用谷歌地球中的经纬度划分，把当地周围划分为多个网格。其次，根据附件4所给数据，用Excel计算11年4月至12年3月的平均风速，风向等指标。随后，使用MATLAB软件
2、编程，根据深圳市气象局所提供的数据，包括风力，降水量，光照时间，得到了08年10年每天白天和夜间的大气稳定度等级。最后，根据大气稳定度，确定污染物浓度函数中的各个参数，进而给出最终的每季度各污染物随空间位置变化的函数，并用MATLAB三维绘图做出函数的图像。建立模型后，为确立合理的经济补偿方案，本组人员查找相关资料，通过计算当地的空气污染指数，按照对居民环境侵害的强度，时间和地点，根据当地的经济赔偿相关标准，设计了合理的经济补偿方案。对于问题2，通过观察附件2中焚烧炉的在线排放监测记录，很容易发现，颗粒物的浓度明显超标，而二氧化硫和氮氧化合物的浓度并不超标。对于颗粒物的监测方法，我们提出了优化
3、措施，建立有关污染物扩散的微分方程模型，即根据透过受污染地区的透过的光的光强来近似判断颗粒物的浓度。为解决故障发生概率下的污染物浓度问题，我们建立概率模型，将排放监测记录中污染物超标的频率近似为除尘装置发生故障的概率。对3月内颗粒物的排放浓度取平均值，在考虑故障发生概率的情况下，再次使用高斯大气扩散模型，得到此时的颗粒物浓度随空间位置的变化规律。最后在发生故障的概率条件下，对问题1中的经济补偿方案进行优化，针对颗粒物所造成的空气污染指数变化实行更为合理的补偿方案。关键词 高斯大气扩散模型，微分方程模型，概率模型目录问题重述1问题分析1模型假设2定义与符号说明2模型的建立与求解3问题1的模型：高
4、斯大气扩散模型3问题2的模型：微分方程模型11模型的优化与评价15参考文献16附录17问题重述垃圾焚烧厂对环境的污染风险与建设投资规模、运行监管力度有直接关系。由于大小的垃圾焚烧炉各有利弊，因此对于不同城市来说，究竟该把大型焚烧厂的建设规模控制在什么水平，这是一个值得研究的课题。深圳市某地点计划建立一个中型的垃圾焚烧厂，计划处理垃圾量1950吨/天（设置三台可处理垃圾650吨/天的焚烧炉，排烟口高度80米，每天24小时运转）。从构建环境动态监控体系、并根据潜在污染风险对周围居民进行合理经济补偿的需求出发，有关部门希望能综合考虑垃圾焚烧厂对周围带来环境污染以及其他危害的多种因素（例如，焚烧炉的污
5、染物排放量、居住点离开垃圾焚烧厂的距离、风力和风向及降雨等气象条件、地形地貌以及建筑物的遮挡程度等等），在进行科学定量分析的基础上，确立一套可行的垃圾焚烧厂环境影响动态监控评估方法，并针对潜在环境风险制定出合理的经济补偿方案。请你在收集相关资料的基础上考虑以下问题：(1) 假定焚烧炉的排放符合国家新的污染物排放标准（参见附件1），根据垃圾焚烧厂周边环境设计一种环境指标监测方法，实现对垃圾焚烧厂烟气排放及相关环境影响状况的动态监控。以你设计的环境动态监控体系实际监控结果为依据，设计合理的周围居民风险承担经济补偿方案。(2) 由于各种因素焚烧炉的除尘装置（如袋式除尘器）损坏或出现其他故障导致污染物
6、的排放增加，致使相关各项指标将严重超标（如：烟尘浓度、二氧化硫、氮氧化物、一氧化碳、二恶英类及重金属等排放超标，附件2给出了一台可处理垃圾350吨/天的焚烧炉正常运作时的在线排放监测记录）。请在考虑故障发生概率的情况下修正你设计的监测方法和补偿方案。问题分析(一) 问题1的分析本问题所研究的环境指标监测方法有着较强的实际意义。通过对周边环境的实时监控，可以让居民了解自己所处的生态环境，并在必要时采取适当的保护措施。本组通过查找相关资料，运用Excel和MATLAB软件，在焚烧炉的排放符合国家新的污染物排放标准的条件下得到了在受影响区域内各个方位的污染物浓度。我们运用高斯大气扩散模型，利用了附件
7、4中的风向和风速，借鉴了深圳市气象局的降水量和日照时间等数据，先得到处理场周围在08年至10年各月的大气稳定度平均值，再用风速，降水，大气稳定度等数据，参考大气扩散模型中的相关计算公式，得到该区域各点的污染物浓度。最后，参考附件3中的生活垃圾焚烧污染控制标准，结合以上对于居民区污染物浓度的分析，研究污染物超标可能导致的对人体的健康威胁，据此设计合理的周围居民风险承担经济补偿方案。(二) 问题2的分析由于问题1中所做的假设“焚烧炉的排放符合国家新的污染物排放标准”过于理想，现实情况下难以做到。由于各种因素焚烧炉的除尘装置（如袋式除尘器）损坏或出现其他故障导致污染物的排放增加，致使相关各项指标将严
8、重超标。附件2给出了3月中烟气的日均排放量，通过计算各污染物严重超标的频率，利用概率模型，可以得到故障发生的概率。在可能发生故障的情况下，再次利用高斯大气扩散模型，在第一问中监测方法可行的条件下，计算出各项指标严重超标时各地的污染物浓度，进而设计出更优化的补偿方案。模型假设1. 假设题目所给的数据真实可靠；2. 在全部高度上，风速是均匀且稳定的。大气污染物与空气没有相对运动，污染物的运动是被动的，且在污染物扩散的过程中，污染物的总质量没有发生变化，符合质量守恒定律，而且污染物本身的化学组成不会与其周围大气中的物质发生化学反应。大气污染物的扩散将在地表被反射出去。3. 大气污染物浓度在 Y 方向
9、和 Z 方向上分布满足正态分布；4. 污染源的污染物浓度分布是连续、均匀、稳定的；定义与符号说明rain_i_j第i年第j个月的降雨量s_i_j第i年第j个月的光照时长R_i第i年的降雨量s_i第i年的光照时长S_i第i年的太阳辐射强度（0代表弱，1代表中，2代表强）w_i第i个月的风速W全年风速day_i第i年白天的大气稳定度night_i第i年晚上的大气稳定度cloud_i第i年云的情况byi大气稳定度为i时的值bzij大气稳定度为i时的值(j=1，2)ave_day按季度的白天大气稳定度平均值ave_night按季度的晚上大气稳定度平均值ave_w按季度的风力的平均值q各项污染源强度c运
10、用高斯模型计算出的浓度表达式IAQIP污染物项目P的空气质量分指数CP污染物项目P的质量浓度值BPHi相应地区的空气质量分指数及对应的污染物项目浓度指数表中与CP相近的污染物浓度限值的高位值BPLo相应地区的空气质量分指数及对应的污染物项目浓度指数表中与CP相近的污染物浓度限值的低位值IAQIHi相应地区的空气质量分指数及对应的污染物项目浓度指数表中与BPHi对应的空气质量分指数IAQILo相应地区的空气质量分指数及对应的污染物项目浓度指数表中与BPLo对应的空气质量分指数模型的建立与求解问题1的模型：高斯大气扩散模型高斯大气污染扩散模型定义如下参数和条件：1、 高斯坐标系原点设定为污染源排放
11、中心点在地表的投射点，X 轴的正方向为该地区的平均风向，Y 轴的正方向为该地区平均风向的切向，在地表的水平面上垂直于 X 轴，Z 轴为垂直于水平地表面 XOY，正方向向上。高斯大气污染物扩散模型的坐标系如图2所示：图12、 理想状态下，高斯污染扩散模型建立的假设a) 大气污染物浓度在 Y 方向和 Z 方向上分布满足正态分布；b) 在全部高度上，风速是均匀且稳定的；c) 污染源的污染物浓度分布是连续、均匀、稳定的；d) 扩散中，污染物的质量是守恒的（不考虑转化等其他因素）；在本题中，需要考虑地形及建筑物等因素对污染物扩散的影响。运用谷歌地球，得到当地周围的地形等高线，近似表示为下图：图2在地形图
12、中可以看出，顺风方向地势较为平坦，海拔较低，居民区较为集中。在本问题中，地势因素对污染物浓度的影响可以近似忽略。定义：我们假定大气污染物的浓度c，而且假定污染物的浓度与其所在的地理位置有关；大气污染源的强度设为 Q；沿 x 轴正方向上的风速设为 u。首先，我们需要利用谷歌地球的经纬度，将污染扩散的预测地区进行网格划分，将该区域划分成a*b的矩阵，矩阵中的每一个点对将对应着一个实际的地理坐标。我们只需要计算出每一个坐标点上的大气污染物浓度，就完成了对该地区大气污染程度的预测。然后，我们还需要创建大气扩散的参数矩阵，以便保存高斯大气污染扩散模型的计算中所需要的参数。第三，我们将读取预测区域当前的气
13、象条件，通过用户输入污染源强度、当地平均风速、所要计算的大气的高度和当时大气稳定程度，来读入当前气象条件。然后根据读入的实际的大气状态来确定高斯系数矩阵中的相应系数计算网格中每一测量点上的污染状况。具体流程如下图所示：图3：高斯大气扩散模型的设立由正态分布假定，可以推定下风向中任意一点的污染物浓度分布，因此大气污染物的浓度 c 将满足下式： （1）其中，方差形式的表达式为： （2） （3）由于单位时间内污染物的质量守恒，由假定可以推出d污染源强的积分公式： （4）4 式中的未知数为：c 代表污染物的浓度，函数 A(x)待定，系数 a，b 待定，通过查表或者将式 4 按照级数展开可得：再代入（2
14、）式，可得： （5） 将 1 式、5 式、6 式代入 4 式中，可得： （6） 将 5 式、6 式代入 1 式中，可得高斯大气扩散模型的最终表达式：（7）由 7 式中的高斯大气污染物浓度表达式可以知道，污染物的浓度除了受到空间位置的影响之外，还会受到大气污染源的强度 q、沿 X 轴正方向的平均风速 u 以及污染扩散系数a,c,f等因素的影响。本系统中，高斯大气污染扩散模型的实现采用的是 matlab 技术，由于 matlab具有高效率的数值计算功能和符号运算功能，可以把编程人员从复杂的数学运算分析过程中解放出来，同时 matlab 还具有比较完备的图形图像处理功能，可以很方便地实现计算结果和编程的可视化，由于 matlab 的语法和思维方式比较接近人类的自然语言,使得初学者容易学习和应用，而且 matlab 提供了功能强大的一些实用的工具，为编程人员提供了大量便于使用的处理工具，所以本项目使用 matlab软件进行编程实现。3、 高斯污染物扩散模型的参数的确定高斯污染物扩散模型的参数是描述大气污染物沿不同方向的扩散程度，其计算公式为：其中参数 a、c、f 是根据不同地区的大气稳定程度而进行确定取值，不同的地区、季节，大气稳定度系数不同。为计算大气稳定度，本组人员对08年至10年的降水
