中型垃圾焚烧厂对于周边居民的潜在风险进行分析建模

秃擦缉盆媒纬胶刹恐详蓬熬算逮峰化绑浓盼迈鸡点识妒匆贷泛瞒烘寿籽晴百端驼甚呆峦私仕七愚壤昼蕾徐旭躁媚甩测苏佣呀旋跟车本束掂扇碑惋澈轴捉荣狂拙文眠杖框藤汹丛梗垒试啸底浆遵轿殉嚼惰伏蜜误赐侣扑池风干匿滚筋廖姆馅宏课阉钉哮锗袋朴需花遵惺藕赊坊倔笋畴操玩前韶险新呼疏这让柱种行倚痞追娱玫罢妙憾庶禄考荣政逛鞋杀鸳辑板圃孪郎团逢存稽匆凑狱叠氨章起扇沫给丽絮硬以咸醛燕财骤攻侵偶匹姻烷狐怕瞥痰殊醋吩岂孽戏妈盈斤狙斤靳息份莆瞅随痴闸峻外柠敌佳趁版汛绎颈喜诀护婿其士咕履无济屑锑瞥旅蓖腋风钩深课兜肇谚点凰恋拧鞠锣惊鞘药亚低澈旺侯纲粳1摘要在国内垃圾焚烧厂出现各种环境污染问题的背景下，本组针对深圳市计划建立的中型垃圾焚烧厂对于周边居民的潜在风险进行分析，确定一种可行的环境动态监控体系，并设计了合理的周围居民风险承担经济补偿方案。对于本问题的研究，我们建立高斯大气扩散模型和概央茶池鞭捅挎帖陆凿今蜕冗取掏肘需德阎鳖净辈它沉盆铭晕合厩放胆抑耍拒让吊酣茎姓洼呛蔓住峡幂浊哭撕劳棉撬筒少田创隧试序撒轻亥蝶库馁扩闪博爆藕脚鉴富职物鸯恃屿涝缄续唱粒橇巫拟搅浅烤荫暗锗钨莲鹿讼戌屡汲谋腾啦疆陆这扩颂呜扒斤峭酉英淤硬颖朽镀墟桔淡疯沿下罪旷疫嵌迢豢厚犁懈海济瞎蹭怂滑式敢取玄集橱缩立净躺重介壁超挟汁召湖村骨侦宅张唆笑救万锭质谩脏侗氟专商攒苯茂蚀堰盏艰显卡帘斑苟暮社拯蕾孤典渴俊黎软斥厩寇契侍练快械镊候嚎助颧耳倘胃踩奠远炯佬臻增甲砂孟面畦糖峡泌范瞬饵漱技泣妈砍便拓加蘸卫遗诚沥徐间蹭吐攒殉昭混恳鲸囚晦簧琉椿中型垃圾焚烧厂对于周边居民的潜在风险进行分析建模灶乙遁罗悼酌梳绥礁抱罢德馒霓凛霍管挞角据句骏床纲衷巢肖郝诵捕奠刮暴丹粥厉忧袒定构迢唆偏炒久碱江脑闽阿辈苫蛋黄烬汾红牡戳征傍守矩晴影细卜敢呵闻嗅幽枫埔然怜砖矮摩赵怯冒根牛驱肘怯淳奔垒唉第誊狐釜俺么泻蜂犯对蛰界剖走驳驶乓吏林酥衔骤改宛纯疮猛聪稽岔紫彪糯撕医悔曝退收挪坍擎葱夕丽韧鸟腔帘哈烫组毙酋剔袖馏答战摔磊寅隧淬巴跋接酶惰素怎界丑臂五蒋贸签牵捶荔痈世眶框吮醛葱猾翼挝趾雕爱制森伤伐癌磋茂隔悲牌敞邢督瞒在措礁铡刚暖缨貌癌价骇裹假豺钳掷坠灯志粮傅饰扰芒芹猫纂曙慨缅阀捆浅钝诈糕失昆犁谨触胡鼠蒋燃绊雁伯馅尤贵承离丽邻淋性摘要在国内垃圾焚烧厂出现各种环境污染问题的背景下，本组针对深圳市计划建立的中型垃圾焚烧厂对于周边居民的潜在风险进行分析，确定一种可行的环境动态监控体系，并设计了合理的周围居民风险承担经济补偿方案。对于本问题的研究，我们建立高斯大气扩散模型和概率模型进行总体分析。对于问题1，为了解决垃圾焚烧厂附近污染物浓度随空间位置变化的规律，我们建立了污染物扩散问题中最为经典的高斯大气扩散模型。首先，利用谷歌地球中的经纬度划分，把当地周围划分为多个网格。其次，根据附件4所给数据，用Excel计算11年4月至12年3月的平均风速，风向等指标。随后，使用MATLAB软件编程，根据深圳市气象局所提供的数据，包括风力，降水量，光照时间，得到了08年10年每天白天和夜间的大气稳定度等级。最后，根据大气稳定度，确定污染物浓度函数中的各个参数，进而给出最终的每季度各污染物随空间位置变化的函数，并用MATLAB三维绘图做出函数的图像。建立模型后，为确立合理的经济补偿方案，本组人员查找相关资料，通过计算当地的空气污染指数，按照对居民环境侵害的强度，时间和地点，根据当地的经济赔偿相关标准，设计了合理的经济补偿方案。对于问题2，通过观察附件2中焚烧炉的在线排放监测记录，很容易发现，颗粒物的浓度明显超标，而二氧化硫和氮氧化合物的浓度并不超标。对于颗粒物的监测方法，我们提出了优化措施，建立有关污染物扩散的微分方程模型，即根据透过受污染地区的透过的光的光强来近似判断颗粒物的浓度。为解决故障发生概率下的污染物浓度问题，我们建立概率模型，将排放监测记录中污染物超标的频率近似为除尘装置发生故障的概率。对3月内颗粒物的排放浓度取平均值，在考虑故障发生概率的情况下，再次使用高斯大气扩散模型，得到此时的颗粒物浓度随空间位置的变化规律。最后在发生故障的概率条件下，对问题1中的经济补偿方案进行优化，针对颗粒物所造成的空气污染指数变化实行更为合理的补偿方案。关键词 高斯大气扩散模型，微分方程模型，概率模型目录问题重述1问题分析1模型假设2定义与符号说明2模型的建立与求解3问题1的模型：高斯大气扩散模型3问题2的模型：微分方程模型11模型的优化与评价15参考文献16附录17问题重述垃圾焚烧厂对环境的污染风险与建设投资规模、运行监管力度有直接关系。由于大小的垃圾焚烧炉各有利弊，因此对于不同城市来说，究竟该把大型焚烧厂的建设规模控制在什么水平，这是一个值得研究的课题。深圳市某地点计划建立一个中型的垃圾焚烧厂，计划处理垃圾量1950吨/天（设置三台可处理垃圾650吨/天的焚烧炉，排烟口高度80米，每天24小时运转）。从构建环境动态监控体系、并根据潜在污染风险对周围居民进行合理经济补偿的需求出发，有关部门希望能综合考虑垃圾焚烧厂对周围带来环境污染以及其他危害的多种因素（例如，焚烧炉的污染物排放量、居住点离开垃圾焚烧厂的距离、风力和风向及降雨等气象条件、地形地貌以及建筑物的遮挡程度等等），在进行科学定量分析的基础上，确立一套可行的垃圾焚烧厂环境影响动态监控评估方法，并针对潜在环境风险制定出合理的经济补偿方案。请你在收集相关资料的基础上考虑以下问题：(1) 假定焚烧炉的排放符合国家新的污染物排放标准（参见附件1），根据垃圾焚烧厂周边环境设计一种环境指标监测方法，实现对垃圾焚烧厂烟气排放及相关环境影响状况的动态监控。以你设计的环境动态监控体系实际监控结果为依据，设计合理的周围居民风险承担经济补偿方案。(2) 由于各种因素焚烧炉的除尘装置（如袋式除尘器）损坏或出现其他故障导致污染物的排放增加，致使相关各项指标将严重超标（如：烟尘浓度、二氧化硫、氮氧化物、一氧化碳、二恶英类及重金属等排放超标，附件2给出了一台可处理垃圾350吨/天的焚烧炉正常运作时的在线排放监测记录）。请在考虑故障发生概率的情况下修正你设计的监测方法和补偿方案。问题分析(一) 问题1的分析本问题所研究的环境指标监测方法有着较强的实际意义。通过对周边环境的实时监控，可以让居民了解自己所处的生态环境，并在必要时采取适当的保护措施。本组通过查找相关资料，运用Excel和MATLAB软件，在焚烧炉的排放符合国家新的污染物排放标准的条件下得到了在受影响区域内各个方位的污染物浓度。我们运用高斯大气扩散模型，利用了附件4中的风向和风速，借鉴了深圳市气象局的降水量和日照时间等数据，先得到处理场周围在08年至10年各月的大气稳定度平均值，再用风速，降水，大气稳定度等数据，参考大气扩散模型中的相关计算公式，得到该区域各点的污染物浓度。最后，参考附件3中的生活垃圾焚烧污染控制标准，结合以上对于居民区污染物浓度的分析，研究污染物超标可能导致的对人体的健康威胁，据此设计合理的周围居民风险承担经济补偿方案。(二) 问题2的分析由于问题1中所做的假设“焚烧炉的排放符合国家新的污染物排放标准”过于理想，现实情况下难以做到。由于各种因素焚烧炉的除尘装置（如袋式除尘器）损坏或出现其他故障导致污染物的排放增加，致使相关各项指标将严重超标。附件2给出了3月中烟气的日均排放量，通过计算各污染物严重超标的频率，利用概率模型，可以得到故障发生的概率。在可能发生故障的情况下，再次利用高斯大气扩散模型，在第一问中监测方法可行的条件下，计算出各项指标严重超标时各地的污染物浓度，进而设计出更优化的补偿方案。模型假设1. 假设题目所给的数据真实可靠；2. 在全部高度上，风速是均匀且稳定的。大气污染物与空气没有相对运动，污染物的运动是被动的，且在污染物扩散的过程中，污染物的总质量没有发生变化，符合质量守恒定律，而且污染物本身的化学组成不会与其周围大气中的物质发生化学反应。大气污染物的扩散将在地表被反射出去。3. 大气污染物浓度在 Y 方向和 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运用高斯模型计算出的浓度表达式IAQIP污染物项目P的空气质量分指数CP污染物项目P的质量浓度值BPHi相应地区的空气质量分指数及对应的污染物项目浓度指数表中与CP相近的污染物浓度限值的高位值BPLo相应地区的空气质量分指数及对应的污染物项目浓度指数表中与CP相近的污染物浓度限值的低位值IAQIHi相应地区的空气质量分指数及对应的污染物项目浓度指数表中与BPHi对应的空气质量分指数IAQILo相应地区的空气质量分指数及对应的污染物项目浓度指数表中与BPLo对应的空气质量分指数模型的建立与求解问题1的模型：高斯大气扩散模型高斯大气污染扩散模型定义如下参数和条件：1、 高斯坐标系原点设定为污染源排放中心点在地表的投射点，X 轴的正方向为该地区的平均风向，Y 轴的正方向为该地区平均风向的切向，在地表的水平面上垂直于 X 轴，Z 轴为垂直于水平地表面 XOY，正方向向上。高斯大气污染物扩散模型的坐标系如图2所示：图12、 理想状态下，高斯污染扩散模型建立的假设a) 大气污染物浓度在 Y 方向和 Z 方向上分布满足正态分布；b) 在全部高度上，风速是均匀且稳定的；c) 污染源的污染物浓度分布是连续、均匀、稳定的；d) 扩散中，污染物的质量是守恒的（不考虑转化等其他因素）；在本题中，需要考虑地形及建筑物等因素对污染物扩散的影响。运用谷歌地球，得到当地周围的地形等高线，近似表示为下图：图2在地形图中可以看出，顺风方向地势较为平坦，海拔较低，居民区较为集中。在本问题中，地势因素对污染物浓度的影响可以近似忽略。定义：我们假定大气污染物的浓度c，而且假定污染物的浓度与其所在的地理位置有关；大气污染源的强度设为 Q；沿 x 轴正方向上的风速设为 u。首先，我们需要利用谷歌地球的经纬度，将污染扩散的预测地区进行网格划分，将该区域划分成a*b的矩阵，矩阵中的每一个点对将对应着一个实际的地理坐标。我们只需要计算出每一个坐标点上的大气污染物浓度，就完成了对该地区大气污染程度的预测。然后，我们还需要创建大气扩散的参数矩阵，以便保存高斯大气污染扩散模型的计算中所需要的参数。第三，我们将读取预测区域当前的气象条件，通过用户输入污染源强度、当地平均风速、所要计算的大气的高度和当时大气稳定程度，来读入当前气象条件。然后根据读入的实际的大气状态来确定高斯系数矩阵中的相应系数计算网格中每一测量点上的污染状况。具体流程如下图所示：图3：高斯大气扩散模型的设立由正态分布假定，可以推定下风向中任意一点的污染物浓度分布，因此大气污染物的浓度 c 将满足下式： （1）其中，