数模论文垃圾焚烧厂的经济补偿问题的数学模型及其分析.docx

垃圾焚烧厂的经济补偿问题的数学模型及其分析姓名 学号 分工潘奇峰 2012300475 建模黄国强 2013302271 编程闫辰毓 2012302229 写作2014年6月8日 目录一.摘要2二.问题的提出3三.问题的分析四.建模过程1）问题一1.模型假设2.定义符号说明3.模型建立4.模型求解2）问题二1.基本假设2.定义符号说明3.模型建立4.模型求解五.模型的评价与改进六.参考文献摘要：垃圾焚烧厂的经济赔偿模型是根据污染物排放对居民产生的潜在影响进行适当的经济赔偿的划分模型，要求根据结合污染物的排放路径，浓度来划分合适的区域对不同区域的居民进行相应的赔偿，同时考虑设备故障情况发生概率下，修正设计的监测方法和补偿方案。对于第一个问题，采用欧拉方程和质量守恒来模拟出不同带状区域的污染物浓度，而由于不同的污染物占到的空气比重是可以推算出的，从而参考空气污染指数（API），根据不同污染物浓度所对应空气比例污染指数（APPI）计算公式中相对应的参数进行计算得到APPI值，最终根据APPPI值的梯度分布对居民进行分梯度的经济赔偿，同时对于流体的成分密度计算还可以用高斯烟羽模型（plume model）对所得结果从另一个角度进行计算从而使模型所得结果更加准确；第二个问题是建立在第一个问题之上的定概率求期望问题，在设备故障的概率情况之下对应了新的赔偿方案，这样的方案中需要确定焚烧炉发生故障的概率模型并且计算出其出故障的概率，在考虑设备出故障的概率情况下因为空气质量的长时间超标或者因为维修而暂停使用造成的空气质量的恶化或者好转对应的赔偿的增加或减少，结合设备正常情况下的赔偿对二者进行加和得到新的赔偿方案。两个问题的模型建立起来较为容易，但是由于是实际问题中的流体扩散，涉及到了很多影响扩散的因素，结果所受的约束条件较多，因而在计算的前期分析占了较大的比重，最终得到了一套可参考性较高的空气检测模式和不同情况下对应的合理经济赔偿方案。关键词：垃圾焚烧 监测方法 经济赔偿 问题的提出： 在垃圾焚烧厂的周围分布着居民区，垃圾焚烧产生的污染物通过烟囱排出，而不同的居民区因为距离焚烧厂的远近以及焚烧厂排出污染物浓度的大小而受到不同程度的影响，加之焚烧厂周边的地形，风向以及污染物扩散过程中减小的梯度等多方面因素的影响，需要建立一个长期监控方案来监测不同居民区区域的污染物浓度，同时对于不同污染影响的居民区进行梯度的经济赔偿，又因为考虑到设备的故障，因此设备故障的概率情况下对于空气质量的监测应该建立新的方案，同时经济赔偿措施也应有所修改，从而不使居民受到损失。这两个问题的出现要求着首先对于空气质量要有一套实时监测方案，并且根据空气质量超标的程度确定不同的经济赔偿方案，另外在第一问的基础上进行概率分析求取平赔偿金额的期望即可。问题的分析： 由题意知，目的是建立一个空气质量监测的模型，根据这个监测模型对于受污染影响的居民进行对应的经济补偿，同时考虑到设备可能出故障的概率情况下新的监测模式和赔偿方案。对于问题一，首先对于最基本的空气质量检测模型的建立，涉及到了污染物的流动传播及其浓度计算，联想到的应该是最基本的扩散运动湍流扩散模型，整个污染物气体扩散过程中满足质量守恒和N-S方程，根据二者可估算不同区域的污染物浓度，从而根据空气质量参数APPI的大小对居民带进行合理划分，不同的区域对应不同的赔偿方案，最后用高斯方程进行辅助验证确定结果的准确性和可参考性；而在问题二中，对于设备概率故障情况下的新方案，建立一个概率模型，在考虑设备出故障的概率情况下因为空气质量的长时间超标或者因为维修而暂停使用造成的空气质量的恶化或者好转对应的赔偿的增加或减少，结合设备正常情况下的赔偿对二者进行加和得到新的赔偿方案。 四：建模过程1） 问题一a.模型假设：1. 污染源不发生移动，即题中所述的三个焚烧炉同时不间断工作，不因客观原因有任何的轮休（除过维修期间）；2. 污染物在扩散路径上呈梯度减小，即满足斐克定律；3. 污染物扩散的空气层面为平流层，不考虑湍流的影响，同时认为地面山面等污染物流动的界面都为无滑移界面；4. 以西南风为主要风向，东北及其周边区域为平原，假定为无风带；5. 污染物排放量不变且每天排放的量相等；6. 流体低俗流动认为是不可压流体；7. 居民区认为是均匀分布；8. 整个过程中，烟囱排放体不发生沉降、分解，不发生任何化学反应等；9. 地面对烟囱排放气体起全反射作用，不发生吸收或吸附作用；b.定义符号说明： 污染气体密度 ,分别是在X,Y,Z方向上的投影 c.模型建立 图1-1如图为Google Earth软件模拟出的垃圾焚烧厂手边的地形环境，由图可以分析得到l 处理厂西南方向有山，而东北方向基本可以看做是平原，风由西南方向吹向东北方向 l 居民区较集中的地方是东部到西北部区域 图 1-2如图1-2所示，为使用ANSYS Fluent简化后的初期模型 ，区域为山区，为模拟的烟囱位置，图中的点都是参照附表4中所给的风向及风俗绘制出的网络图格，在之后的计算中通过表中数字的拟合得到风速不同情况下对应的污染物传播以及每个点出的污染物及其对应的主要成分的浓度大小。 图1-3如图1-3所示，由ANSYS Fluent得到的流场矢径图，图中所对应的空白区域为山所在的位置，所对应的空白区域为烟囱所在的位置，通过这个图模拟出了风从山上吹向地面的过程中不同位置的流速大小，颜色越浅（趋近于绿色）的地方流速越大，例如山两边的绿色区域为流速最大的地方，颜色越深（趋近于蓝色）的地方流速越小，例如整张图的边缘部分都为流速较小的地方。d.模型求解（1）基本理论推导质量守恒定律-连续方程 式1-1流体低速运行可近似认为是不可压流体，满足牛顿第二定律 式1-2 式1-3欧拉公式 式1-4每一种气体组分均满足质量方程，因此得到密度方程 式1-5 其中C为体积浓度，为质量浓度，为扩散系数。 图1-4 考虑风可能来自于八个方向，以及烟囱高高度，我们利用Ansys Geometry软件，将空间抽象化成为一个八棱柱，其高度约为烟囱高度的5倍，宽度约为3Km，以保证在柱体的边缘基本不受流场的影响。 对模拟化之后的空间，利用Ansys ICEM软件，进行进一步处理按照来风方向已经各个面的特性将模型包含的12个面进行进一步的区分：区域名称区域功能各面初始条件设定E东风将从盖面吹入模型、西风将从该面吹出P为一个大气压V=VeEN东北风将从盖面吹入模型、西南风将从该面吹出P为一个大气压V=VenES东南风将从盖面吹入模型、西北风将从该面吹出P为一个大气压V=VesN北风将从盖面吹入模型、南风将从该面吹出P为一个大气压V=VnS南风将从盖面吹入模型、北风将从该面吹出P为一个大气压V=VsW西风将从盖面吹入模型、东风将从该面吹出P为一个大气压V=VwWN西北风将从盖面吹入模型、东南风将从该面吹出P为一个大气压V=VwnWS西南风将从盖面吹入模型、东北风将从该面吹出P为一个大气压V=VwsTOP模拟该区域的天空P为一个大气压,V=0LOW模拟该区域的大地该面为well形（无滑移边界）YC模拟烟囱外壁以及较高的山峦，作为障碍物该面为well形（无滑移边界）YUAN污染源将从该面，以一定的流量垂直向上进入模型V=Vo 表1-1以每个细小网格的中心点的数据来近似代替该网格所包围的空间内的流场数据。如果网格越小，则模型越精细，计算越准确。综合考虑模型要求以及现有计算机的性能之后，我们考虑建立700W个网格进行计算。网格建立完全后 完整的网格模型如下： 图1-5网格建立成功之后，我们利用Ansys Fluent 软件进行求解。结合原有的假设，对于流体我们选择Laminar（层状）模型。经过500次的数据迭代之后，我们得到了最终的效果图： 图1-6（2）所得图像分析经过ANSYS Fluent分析得到 图1-7 由图1-4可以看到，中间红色区域中心是烟囱，外围的曲线为污染物分布的等浓度线，越往外污染气体浓度越低，而且单位长度上的浓度梯度越小，且右边的污染区域大于左边的污染区域。这样的结论为之后按污染物浓度大小划分区域奠定了基础。拟合曲线（月） 图1-8NO拟合曲线（月） 图1-9如图1-5和图 1-6所示，在经过趋近拟合之后发现 和NO在污染气体中所占的比重随着时间的变化曲线与正态分布有较好的相似性，在这样的结果基础上，提出一个参数作为一定污染浓度所对应的赔偿金额即可以得到一个合理的经济赔偿方案。 图1-10 如图1-7所示为由高斯公式确定出的污染物区域分布等值图，之后的经济赔偿方案主要根据不同的污染指数划分的污染区域进行赔偿，每个区域之内对应一个赔偿金额。（3）赔偿方案的确定污染指数污染物浓度（毫克/立方米）APISO2（日均值）NO2（日均值）PM10（日均值）CO（小时均值）O3（小时均值）500.0500.0800.050 50.1201000.1500.1200.150100.2002000.8000.2800.350600.4003001.6000.5650.420900.8004002.1000.7500.5001201.0005002.6200.9400.6001501.200 表1-2如表1-1所示为空气中空气污染指数API 计算的重要指标及其所对应的API值，而API的计算方法如下：设I为某污染物的污染指数，C为该污染物的浓度 式1-6式中： 与 ：在API分级限值表（表1-1）中最贴近C值的两个值，为大于C的限值，为小于C的限值。 与 ：在API分级限值表（表1-1）中最贴近I值的两个值，为大于I的值，为小于I的值。但是API指数作为衡量空气中的污染物浓度给出一个污染指数是科学的，但是因为不同的污染物对于人体的有害程度是不一样的，而产生相同的危害所需要的不同气体的浓度也是不同的，而且在本题中污染物中不同气体的浓度差异比较明显，所以API参数不能完全满足要求，考虑到不同污染物的浓度不同带来的危害不同的因素，我们定义一个新的指标空气比例污染指数APPI（Air Proportion Population Index）。其计算公式如下： 式1-7 其中为污染程度最大的气体成分浓度对应的毒性指数， 为每种气体成分浓度对应的毒性指数。 气体成分毒性指数 0.126 0.3418颗粒物0.5355表1-3如表1-2所示，表格列出了重要的三种污染物的毒性指数，如果需要计算更多气体的污染指数可以进行进一步计算使得结果更加准确同时对于后面的赔偿金额计算更具参考性。（4）赔偿金额的计算公式： 式1-8 其中为第i个区域对应的赔偿金额 K为污染补偿系数，依据实际的污染气体密度确定为第i个区域对应的APPI指数 为可接受的APPI指数（在本次数据模拟中90m高的烟囱对应的为0.5）M为所有赔偿金额总数2） 问题二a. 模型假设1. 表示维修设备所用的天数且为定值；2. 每年的天数按照365计算；3. 维修期间设备不运行，认为空气质量均达标，此时不进行赔偿；（1）焚烧炉出现故障的概率模型设焚烧正常的情况为事件A，正态分布下出现故障的情况为事件B，事件B对应的概率为P（X）P（）=P（B）-P（AB） 0 式2-1又由于P(AB) 0故可认为P（）=P（B） 式2-2即假设中该事件用正态分布趋近是合理的，可认为该小概率事件服从正态分布。 （2）日平均APPI指数的监测模拟 由第一问的图像及结果可以得到主要的污染气体浓度在每个月基本可以认为是服从正态分布的，其对应的APPI指数亦服从正态分布，由此推测每天因为设备故障污染物浓度超出正常标准的概率服从正态分布。 取样 为6:0018:00时间段内的气体APPI指数，取三天的数据对同一指标进行分析，得到的结果可以近似认为 式2-3 即单日的APPI指数超标的概率可以用正态分布来近似趋近，这也辅助证明了第一个模型的建立是合理的。b. 模型计算（1）新的监测方法对每日APPI指数进行监测，如果出现因故障引起的超标，第二天立刻进行维修即可。（2）新的赔偿方案 假设正态分布下焚烧厂设备出现故障的概率为p，对应的设备正常的概率为1-p，则每月可能出现的设备出现故障天数： 式2-4因为维修设备不进行工作的天数设备正常情况下的赔偿金额： 式2-5 因为故障维修产生的赔偿金额变动： 式2-6 其中为空气质量超标之后对应的APPI指数新的赔偿方案下的金额： 式2-7五：模型的改进与评价 本题的模型建立过程中主要考虑到了地理和物理方面的因素对于结果的影响，在建模过程中采用的质量守恒，欧拉公式的推导并且利用高斯烟羽模型对污染气体的浓度的分布进行了较好的模拟，并且在计算过程中拟合出了每个点的动态浓度，对于监测方案的建立有很好的参照作用，最后根据得到的污染物浓度等值图对污染区域进行划分确定出45个赔偿带，而赔偿计算过程中计参照的指标是新建立的APPI指数，由于API指数的主要目的是划分出大气污染等级，在本题中污染气体的浓度都较大，且类型不单一，因此相比较API指数，APPI指数更能反映污染气体的危害大小从而进行赔偿。 而第二问的大部分思想方法是建立在第一问的基础之上的，新的监测方案仍然是取每天的APPI指数进行分析，如果发现超标立刻检查是否有故障，第二天即刻进行修缮，至于新的赔偿方案，需要考虑到设备出故障的概率，维修前APPI超标所带来的赔偿差值，维修过程中因为设备停止运行产生的赔偿差值等多个因素，最后确定出了新的赔偿方案。 图