数模论文垃圾焚烧发电.docx
- 文档编号:24549536
- 上传时间:2023-05-28
- 格式:DOCX
- 页数:21
- 大小:380.81KB
数模论文垃圾焚烧发电.docx
《数模论文垃圾焚烧发电.docx》由会员分享,可在线阅读,更多相关《数模论文垃圾焚烧发电.docx(21页珍藏版)》请在冰豆网上搜索。
数模论文垃圾焚烧发电
2014深圳杯全国大学生数学建模夏令营
承诺书
我们仔细阅读了中国大学生数学建模夏令营的竞赛规则.
我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。
我们知道,抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。
如有违反竞赛规则的行为,我们将受到严肃处理。
我们参赛选择的题号是(从A/B/C/D中选择一项填写):
C
我们的参赛报名号为(如果赛区设置报名号的话):
所属学校(请填写完整的全名):
参赛队员(打印并签名):
1.
2.
3.
指导教师或指导教师组负责人(打印并签名):
日期:
2014年5月28日
赛区评阅编号(由赛区组委会评阅前进行编号):
2014深圳杯全国大学生数学建模夏令营
编号专用页
赛区评阅编号(由赛区组委会评阅前进行编号):
赛区评阅记录(可供赛区评阅时使用):
评
阅
人
评
分
备
注
全国统一编号(由赛区组委会送交全国前编号):
全国评阅编号(由全国组委会评阅前进行编号):
垃圾焚烧厂的经济补偿问题
摘要
在垃圾焚烧厂运行监管方面,目前主要是在垃圾焚烧厂内进行测量监控,缺少从周边环境视角出发的外围动态监控。
本文采用一种大气污染物的扩散模型-----高斯模式,对垃圾焚烧厂周围居民点的污染物浓度进行合理的预测,进而形成各居民点污染程度的评价,对不同等级的污染情况制定不同的赔偿标准。
本文主要采用微分方程方法,研究2011年4月到2012年4月垃圾焚烧厂所在地每天的风向、风速、降雨量等记录数据,以焚烧炉的污染物排放量、居住点离开垃圾焚烧厂的距离以及风速、风向、降雨量等气象因素为主要影响因素,研究这些因素对各居民点污染物浓度的影响,建立扩散方程,形成每个居民点每天每种污染物浓度的记录,进而计算空气污染指数API,然后以API为依据,形成污染程度的等级划分,进而形成不同污染等级的赔偿标准。
本文中用高斯模型计算得到的各地污染物浓度与实地监控得到的数据基本吻合,用计算得到的数据和我们检验用到的实测数据分别代入后续的API计算公式,所得到的污染分级大致形同,说明本文建立的模型与实际情况较为相符,具有可行性。
对于后续模型的改进,我们提出加入地形起伏等因素的影响(考虑到建筑物高度对污染物扩散的影响相对有限,所以在高度对模型影响时只考虑地形因素),这需要在原有扩散方程中加入kh因子,其中k为影响系数(假设h对浓度的影响是线性的,后续根据模型与实际浓度吻合情况考虑h的非线性影响),但由于时间限制,我们只提出了初步的想法,没能最终形成修正后的扩散方程。
关键词:
高架连续点源扩散模式,空气污染指数API,动态监控
1、问题重述
1.1问题背景
“垃圾围城”是世界性难题,在今天的中国显得尤为突出。
数据显示,目前全国三分之二以上的城市面临“垃圾围城”问题,垃圾堆放累计侵占土地75万亩。
因此,垃圾焚烧正逐步成为中国垃圾处理的主要手段之一。
然而,由于政府监管不力、投资者目光短浅等多方面的原因,致使前些年各地建设的垃圾焚烧电厂在运营中出现了环境污染问题,给垃圾焚烧技术在我国的推广造成了很大阻力。
这就需要客观上对垃圾焚烧厂所排出的废气对本地区空气污染情况进行评价。
但是在垃圾焚烧厂运行监管方面,目前主要是在垃圾焚烧厂内进行测量监控,缺少从周边环境视角出发的外围动态监控,因而难以形成为民众所信服的全方位垃圾焚烧厂环境监控体系。
基于上述情况,根据现有数据,运用数学建模的方法,对城市垃圾焚烧厂空气污染做出准确的分析,在进行科学定量分析的基础上,确立一套可行的垃圾焚烧厂环境影响动态监控评估方法和具体的赔偿标准
1.2问题提出
深圳市某地点计划建立一个中型的垃圾焚烧厂,计划处理垃圾量1950吨/天(设置三台可处理垃圾650吨/天的焚烧炉,排烟口高度80米,每天24小时运转)。
从构建环境动态监控体系、并根据潜在污染风险对周围居民进行合理经济补偿的需求出发,有关部门希望能综合考虑垃圾焚烧厂对周围带来环境污染以及其他危害的多种因素(例如,焚烧炉的污染物排放量、居住点离开垃圾焚烧厂的距离、风力和风向及降雨等气象条件、地形地貌以及建筑物的遮挡程度等等),在进行科学定量分析的基础上,确立一套可行的垃圾焚烧厂环境影响动态监控评估方法,并针对潜在环境风险制定出合理的经济补偿方案。
要求在收集相关资料的基础上考虑以下问题:
(1)假定焚烧炉的排放符合国家新的污染物排放标准,根据垃圾焚烧厂周边环境设计一种环境指标监测方法,实现对垃圾焚烧厂烟气排放及相关环境影响状况的动态监控。
以设计的环境动态监控体系实际监控结果为依据,设计合理的周围居民风险承担经济补偿方案。
(2)由于各种因素焚烧炉的除尘装置(如袋式除尘器)损坏或出现其他故障导致污染物的排放增加,致使相关各项指标将严重超标(如:
烟尘浓度、二氧化硫、氮氧化物、一氧化碳、二恶英类及重金属等排放超标)。
在考虑故障发生概率的情况下修正监测方法和补偿方案。
2、模型假设
1)以右手坐标系(食指—x轴;中指—y轴;拇指—z轴)建立高斯模型的坐标系,原点:
为无界点源或地面源的排放点,或者高架源排放点在地面上的投影点;x为主风向;y为横风向;z为垂直向
2)全部高度风速均匀稳定;
3)源强是连续均匀稳定的;
4)扩散中污染物是守恒的(不考虑转化);
5)污染物浓度在y、z风向上分布为正态分布;
6)不考虑建筑物高度对污染物浓度的影响,即z=0,
3、符号说明
:
平均风速,m/s;
:
源强,g/s;
:
侧向扩散参数,污染物在y方向分布的标准偏差,m;
:
竖向扩散参数,污染物在z方向分布的标准偏差,m;
:
垃圾焚烧厂烟囱的高度;
:
横风向坐标;
:
主风向坐标;
、
:
在API分级限值表(表1)中最贴近C值的两个值,C大为大于C的限值,C小为小于C的限值;
、
:
在API分级限值表(表1)中最贴近I值的两个值,I大为大于I的值,I小为小于I的值;
:
测点逐时污染物浓度,n为测点的日测试次数;
:
某一污染物的日均值;
:
第i种污染物的污染指数;
4.问题分析
4.1问题一的分析
空气中的污染物通常都会受到风力降雨等自然因素以及扩散距离的影响,风力和降雨越大,距离污染源越远,污染物的浓度就会越小。
利用高架连续点源扩散模式,计算在风力和降雨影响下不同距离处居民区空气污染物的浓度,并把污染浓度换算成API指数。
以API为依据,结合当地的经济发展水平和国家相关规定算出垃圾焚烧厂附近居民应得的经济补偿。
下图为问题一建模的流程图.
4.2问题二的分析
由于各种因素焚烧炉的除尘装置(如袋式除尘器)损坏或出现其他故障导致污染物的排放增加,致使相关各项指标将严重超标(如:
烟尘浓度、二氧化硫、氮氧化物、一氧化碳、二恶英类及重金属等排放超标。
所以在焚烧炉因为故障而导致污染物排放量超出国家标准的时候,应该另外评估监测指标和赔偿标准。
要考虑到故障机器个数对污染物浓度的不同影响,就需要重新考虑污染物初始浓度,计算出此时的API值,与正常工作条件下的API值进行比较,最后将故障的概率考虑成权重从而计算出最后的补偿金额。
5问题的建立与求解
5.1对问题一的建模与求解
5.1.1扩散方程的建立
下图是高斯扩散模型的坐标系,
(5.1.1)
式(5.1.1)是高架连续点源扩散模式,为简单起见,我们忽略高度对居民点污染物浓度的影响,即令z=0,得到地面浓度模式:
(5.1.2)
由附件4及深圳市气象局网站提供的2011年4月至2012年4月的数据记录表,得到每天的风向、风速、降雨量等数据,由附件一得到垃圾焚烧厂烟囱每天排出的污染物浓度。
我们要根据烟囱排出的污染物浓度进一步得到各污染物的源强q.
对于源强的计算,假设烟囱是一个规则的圆柱体,底面直径为5m,高为80m.
烟囱形状如左图所示,则
其中,q为源强,单位是g/s;c为污染物浓度,单位是mg/m3;v为烟囱体积,单位是m3,r=5m为底面半径,h=80m为烟囱高度;t为烟雾在排出之前在烟囱内平均停留的时间,由附件五知t大约为3s.
通过谷歌地图我们找到垃圾焚烧厂附近居民点的分布情况,我们选定6个居民区作为监控对象,这6个居民区分别是罗田社区、东方社区、光明新区西田社区、碧光社区、辅城坳社区、和禾花社区,为了叙述方便,将这6个社区分别记为社区A、B、C、D、E、F,他们距离焚烧厂的距离分别为3km、4.3km、5.5km、4.1km、1.6km、2.5km.
5.1.1.1大气扩散系数(
,
)的确定
扩散参数是表征湍流扩散剧烈程度的物理量,是影响污染物浓度的重要参数。
首先我们考虑使用P-G曲线法确定大气扩散系数,该法的要点是首先根据帕斯奎尔划分大气稳定度的方法来确定大气稳定度级别;然后从图4-4和图4-5中查得对应的扩散参数σy和σz;最后将σy、σz代入前面介绍的高斯扩散模式中,就可估计出各种情况下的浓度值。
帕斯奎尔在1961年推荐一种仅需要常规气象观测资料就能估算σy,σz的方法,吉福德(Gifford)进一步将它制成应用更方便的图表。
应用观测到的风速、云量、云状和日照等天气资料,将大气扩散稀释能力分为6个等级:
A—极不稳定,B—不稳定,C—弱不稳定,D—中性,
E—弱稳定,F—稳定。
若稳定级别为A~B,则表示按A、B级的数据内插。
可见,稳定度需要根据各居民点每天的风速和降雨量(是否阴天)确定。
下一步,要根据P-G曲线得到
,
。
但在实际的建模过程中我们发现利用P-G图虽然比较直观,但读数精度差,且不便于使用计算机计算,因此我们考虑使用一组经验公式作为代替,通过查阅资料,我们找到了一种较为精确的公式【1】,
(5.1.3)
这个公式是在Martin【2】提出的
(5.1.4)
基础上进一步拟合得到的,该经验公式和拟合结果可以在0.01~100.00km范围内较好地表示水平和垂直方向上的大气扩散系数。
接下来,我们需要确定(5.1.3)中系数a、b、c、g、e。
公式(5.1.3)在水平方向上扩散系数的拟合结果如下表:
在垂直方向上扩散系数的拟合结果为:
5.1.1.2计算出6个居民点在高斯模型中的坐标(x,y)
以垃圾焚烧发电厂为坐标圆心,将6个居民区的方位用角度表示,即用0°表示正东方向,45°表示东北方向,90°表示正北方向,以此类推,结合各居民点距离中心点的距离建立居民区的空间分布图。
依据上文假设知道,高斯模型里的x轴方向为正风向,在不同风向下,有不同的x轴,也就有不同的六组位置坐标。
注:
0°(180°)表示x轴正向为0°,对应的,此时风向为180°,即西风。
上表即不同风向(x轴)下居民点的地理坐标。
5.1.1.3某种污染物浓度的最终确定
由附件四的风向数据得到垃圾焚烧发电场附近的风向频率图,
显然,225°,即西南方向的风频最高,相应的,位于45°,即东北方向的居民点受污染情况应该最严重,理论上,我们最终由高斯方程得到的东北方向居民点E的污染物浓度应
该最高(当然还应考虑每个风向上平均风速的大小,故最终得到的污染最严重居民区未必一定是E)。
我们将垃圾焚烧发电厂释放的各种污染物浓度、各居民点每天在各自大气稳定度级别下的
,
值,以及各居民点每天确定风向下地理坐标值(x,y)带入式(3)和式
(1),最终将得到每个居民点每天的各种污染物浓度。
5.1.2空气污染指数API的计算
API(AirPollutionIndex的英文缩写)是空气污染指数,我国城市空气量日报API分级标准及不同空气污染指数下空气质量如下表:
表1空气污染指数对应的污染物浓度限值
污染指数
污染物浓度(毫克/立方米)
API
SO2
(日均值)
NO2
(日均值)
PM10
(日均值)
CO
(小时均值)
O3
(小时均值)
50
0.050
0.080
0.050
5
0.120
100
0.150
0.120
0.150
10
0.200
200
0.800
0.280
0.350
60
0.400
300
1.600
0.565
0.420
90
0.800
400
2.100
0.750
0.500
120
1.000
500
2.620
0.940
0.600
150
1.200
表2空气污染指数范围及相应的空气质量类别
空气污染指数
API
空气质量
状况
对健康的影响
建议采取的措施
0~50
优
可正常活动
51~100
良
101~150
轻微污染
易感人群症状有轻度加剧,健康人群出现刺激症状
心脏病和呼吸系统疾病患者应减少体力消耗和户外活动
151~200
轻度污染
201~250
中度污染
心脏病和肺病患者症状显著加剧,运动耐受力降低,健康人群中普遍出现症状
老年人和心脏病、肺病患者应在停留在室内,并减少体力活动
251~300
中度重污染
>300
重污染
健康人运动耐受力降低,有明显强烈症状,提前出现某些疾病
老年人和病人应当留在室内,避免体力消耗,一般人群应避免户外活动
设I为某污染物的污染指数,C为该污染物的浓度。
则:
(5.1.5)
式(5)中:
C大与C小:
在API分级限值表(表1)中最贴近C值的两个值,C大为大于C的限值,C小为小于C的限值。
I大与I小:
在API分级限值表(表1)中最贴近I值的两个值,I大为大于I的值,I小为小于I的值。
API的计算过程如下:
求某污染物在某个居民点每一测点的日均值
(5.1.6)
式中:
Ci为测点逐时污染物浓度,n为测点的日测试次数
2.求某一污染物整个居民点的日均值
(5.1.7)
式中:
l为整个居民点的监测点数。
2.将各污染物的日均值分别代入API基本计算式所得值,便是每项污染物的API分指数。
3.选取API分指数最大值为某个居民点的API
对于上述过程,我们简化为:
用高斯模型算出的各污染物浓度替代
,最终也是选取最大的API值作为各居民点用于评价污染程度的API.
5.1.3制定赔偿标准
环境补偿是指对环境产生不良影响或破坏的生产者、开发者、经营者应对环境污染、生态破坏进行补偿,对环境资源由于现在的使用而放弃未来价值进行补偿。
建立环境补偿机制应遵循以下原则:
一是谁利用谁补偿,谁受益或者谁损害谁付费;二是有利于受损地区与受益地区的共同发展;三是满足需要与现实可行相结合;四是有利于促进环卫设施的顺利建设正常运行。
具体的赔偿分为两部分,一部分为固定赔偿,实施补偿时按垃圾处理量补偿方法操作,参考国内某垃圾焚烧发电厂【3】每吨垃圾处理的补偿费(补偿对象为所有受影响居民)为40元/吨。
第二部分为动态赔偿,即根据每个居民点每天的污染程度进行不同金额的赔偿,赔偿标准如下
(5.1.8)
其中I为每个居民点每天的API值,Imin和Imax为每天各居民点中API的最小值和最大值。
即每天赔偿的上限是10元【3】,API越大,赔偿金额越多。
处理2011年4月至2012年4月的数据,我们得到每个居民点每人总赔偿金额如下:
动态赔偿
固定赔偿元/吨*
A
B
C
D
E
F
赔偿数额/元
1515
1414
1211
1412
2759
2403
2
日均赔偿金额/元
11.1
10.8
10.2
10.8
14.9
13.8
*固定赔偿金额为垃圾焚烧发电厂每处理一吨垃圾需要支付的赔偿,假设每个居民区有100名居民,已知垃圾处理厂每天处理1950吨垃圾,故每户居民每天平均得到的固定赔偿为2*1950/600=6.5元。
日均赔偿金额为动态赔偿每天的平均值加上换算后的每天平均固定赔偿。
最终从赔偿金额来看,风频最高方向的E居民区和距离垃圾焚烧发电厂最近的F居民区受污染情况最为严重,年赔偿金额最高;距离焚烧厂最远的C居民区,相对来说受污染相对较少,得到的赔偿金额也最少。
这与我们先前的预测和实际情况相符,说明我们的模型具有一定的可行性。
接下来,我们要考虑由于各种因素导致焚烧炉的除尘装置(如袋式除尘器)损坏或出现其他故障导致污染物的排放增加,致使相关各项指标将严重超标(如:
烟尘浓度、二氧化硫、氮氧化物、一氧化碳、二恶英类及重金属等排放超标)。
在这种故障发生概率的情况下修正我们的监测方法和补偿方案。
5.2、对问题二的建模和求解
污染物的传播模型和问题一的传播模型相同,仍然采用高斯模型中的高架点源传播模式,只是在故障状态下的污染物的初始浓度不同。
在考虑机器故障的条件下,参照附件2中,350吨的焚烧炉的污染物浓度,重新计算初始污染物浓度。
(5.2.1)
其中k(0,1,2,3)为出现故障的机器数目,
为一台350吨焚烧炉的污染物的排放某种污染物的浓度,
为国家标准的某种污染物的最大浓度值。
风向,降雨量等限制条件和问题一相同。
3台焚烧炉故障的概率应该服从二项分布,计算公式为
(5.2.2)
我们将概率p(k)和1-p(k)的数值分别作为故障状态和正常状态下补偿金额的权重,从而计算出最后的补偿金额。
(5.2.3)
式(5.2.3)中
为最终的补偿金额,
为焚烧炉正常工作时的补偿金额,即问题一的模型得到的补偿金,
为故障状态下补偿金额。
查阅相关资料得知垃圾焚烧炉一年内正常工作的时间大约为8000小时,有760小时时间处于故障停工状态。
所以一台焚烧炉的故障率近似等于0.0868,即p=0.0868。
表5.2.1年赔偿金额
机器故障数目
动态赔偿
固定赔偿元/吨*
A
B
C
D
E
F
1
1515
1414
1211
1418
2759
2403
3
2
1515
1414
1211
1418
2759
2403
3.5
3
1515
1414
1211
1418
2759
2403
4
居民因为机器故障而应该得到的固定补偿应该为
K=1:
3*1950/600=9.8元。
K=23.5*1950/600=11.4元
K=34*1950/600=13元
下表为故障状态下补偿金额
表5.2.2故障状态下日均补偿金额表.
机器故障数目
日均赔偿金额
人均固定赔偿/元
A
B
C
D
E
F
1
14.4
14.1
13.5
14.1
18.2
17.1
9.8
2
16
15.7
15.1
15.7
19.8
18.7
11.4
3
17.6
17.3
16.3
17.3
21.4
20.3
13
由公式5.2.2可以计算出P
(1)=0.22p
(2)=0.02p(3)=6.5*10^-4
所以在考虑机器故障的条件下,最后的补偿金额应该是:
表5.2.3最终补偿金额
机器故障数目
日均补偿金额
A
B
C
D
E
F
1
11.8
11.5
10.9
11.5
15.6
14.5
2
11.2
10.9
10.3
10.9
15.0
13.9
3
11.1
10.8
10.2
10.8
14.9
13.8
6、模型的评价
6.1对问题一的评价
在考虑污染物浓度扩散的时候,由于建筑物等物体的干扰很微弱且无法量化计算,所以我们对问题作了简化。
通过把天气因素和居民区的位置综合考虑,计算出居民区处的API指数,并根据居民区所处位置的API指数情况进行相应的经济补偿。
我们的模型优点是进行动态监测时,只需要将实时的污染物排放情况录入该模型,则污染物浓度和相应赔偿情况就会自动生成。
2011年4月至2012年4月的数据表明E、F居民点的补偿额度最高,C点补偿额度最小,这是因为E风频最高和F距离垃圾焚烧发电厂最近,两者受污染情况最为严重,所以年赔偿金额最高,C点距离垃圾焚烧发电厂最远,受污染程度最小,所以我们的模型具有一定的可靠性,这也是模型的优点之一。
在问题一的假设下,我们可以计算出焚烧厂一年需要支付的补偿金额为2350000元左右。
根据附件5可以知道,采用“倾斜往复式逆推机械式炉”时,一吨垃圾的利润为85.5元,一吨垃圾的政府补贴为最少为45元,所以在不考虑工资,机器损耗等经营成本的情况下,焚烧发电一年的毛利润为8千万左右,所以赔偿金额只占工厂产值的3%,工厂是可以接受的。
6.2对问题二的评价
这个问题本身及其建模过程都非常简单,我们着重考虑的是机器正常运行过程中出现故障后污染物的增量和相应经济补偿标准计算方法的改变。
这种分析对于一个模型,特别是在实际生产工作中是很重要的。
由表6.1、表6.2、表6.3可知在考虑到焚烧炉可能发生故障的情况下个居民点的补偿额度都有相应增加。
本文的两个问题的分析基本上都是建立在直接赔偿的方式进行,即结合项目规模和当地经济发展水平确定每吨补偿若干元和受影响的人口数,从运营单位收入(主要是垃圾处理费收入)或政府财政中补偿给当地受影响的社区。
考虑到这种补偿方式过于单一,我们想增加间接补偿的方式,间接补偿方法包括以下几种方式:
垃圾处理副产品回馈(包括灰渣、废铁、电、汽、热);公共设施共享;发包工程、采购、用工的优先;实行部分税收减免;经营收入下降的补偿以及其他间接补偿方式(包括免费体检、扶贫、助学等)。
基于以上分析,我们的模型进一步改进的方向是采用直接补偿结合间接补偿的方式进行补偿。
直接补偿是按人口和垃圾焚烧厂处理垃圾当量进行补偿,即对补偿范围内的居民和土地一定标准进行长期补偿。
该方式一方面可使受影响的居民个体直接受益,另一方面也弥补了项目周边补偿范围内土地利用价值降低而减少的经济收入,居民感受度良好。
间接补偿推荐按照公用设施共享的方法进行补偿,该补偿方式不仅使当地居民受益,也使焚烧厂的职工受益。
为了尽量减少污染物的排放,在更换仪器的时间上,工厂应该将不同设备的检修时间错开,避免过多的焚烧炉处于不可工作状态。
其中袋式除尘器的使用寿命经查询在3年左右,一般在20个月左右进行更换。
通过我们的计算,也可以看到其实污染物对居民区的影响是很小的,即使在仪器出故障的情况下,污染物对人的影响也是很小的,居民没有必要过度惊慌。
在补偿方案的选择上,发电厂也可以选择为周边居民提供每月的免费用电数额,或者现金加免费用电额度的补偿方案。
7、参考文献
[1]李玉平计算大气扩散系数的一组经验公式[M].北京:
北京理工大学学报第29卷第10期,2009.
[2]Martin
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数模 论文 垃圾 焚烧 发电
![提示](https://static.bdocx.com/images/bang_tan.gif)