基于calpuff和feps模型的秸秆焚烧烟羽扩散模拟

基于calpuff和feps模型的秸秆焚烧烟羽扩散模拟

近年来，中国采取了大量措施来减少对农业和小麦的燃烧，但这种焚烧现象仍然难以控制。随着大量pm、有机碳、氮、苯和二环芳烃等有害气体燃烧，不仅危害人类健康，而且污染严重。2007年5月，中国北方和江淮流域的主要粮食产区出现了大规模的稻草焚烧，这显著降低了北京和周边城市的相对湿度。从9月14日到8月18日，由于大量燃烧稻草，成都国家执行局无法取消10个订单，22个港口的延误延误，1000多名乘客被安置在港口。目前国内秸秆焚烧的研究主要有遥感对秸秆焚烧的监测、排放因子的试验研究以及排放清单的估算等,但很少对污染物的扩散进行模拟.祝斌等对农作物秸秆燃烧时PM2.5排放因子进行了试验研究,并分别对水稻、小麦和玉米等农作物秸秆在明火和闷火状态下的燃烧情况进行了试验,测定了燃烧过程中PM2.5的排放因子.王书肖等采用排放因子法,建立了我国秸秆露天焚烧一次大气污染物的排放清单,并分析了其时间和空间的分布特征.厉青等利用卫星遥感技术监测了2007年6月全国秸秆焚烧状况,并以某市为例,结合气象资料分析秸秆焚烧对环境空气质量的影响.方萌等介绍了“3S”技术在农作物秸秆焚烧监测中应用的一般原理和方法,同时还介绍了国家卫星气象中心近年来采用“3S”技术监测农作物秸秆焚烧的概况.国外对于秸秆焚烧的研究除了有对排放因子的计算、遥感的监测外,还有对焚烧过程中产生污染物的扩散进行的计算机模拟.JAIN等对太平洋海岸西北部秸秆焚烧的烟气扩散预测系统(ClearSky)进行了介绍,ClearSky是数值烟气扩散预测系统,通过确定秸秆焚烧情景,利用MM5中尺度气象预测模型和CALPUFF模拟系统,模拟秸秆焚烧产生的PM2.5浓度扩散情况.CHOI等利用CALPUFF/CALMET/MM5模拟系统,对美国墨西哥边境秸秆焚烧排放的PM10浓度进行了扩散模拟,并与PM10的相关环境标准进行了对比.据秸秆焚烧遥感监测日报发布的2009年全国秸秆焚烧分布遥感监测结果发现,安徽和河南两省6月中上旬的秸秆焚烧现象均十分严重.针对秸秆焚烧严重的现象,笔者借鉴国内外模拟颗粒物扩散的经验,整合了CALPUFF模型与FEPS模型,并对CALPUFF模型的部分程序进行了重新编译,模拟了安徽和河南两省2009年6月3日与5日秸秆焚烧产生的ρ(PM10)扩散情况,得到了PM10的烟羽扩散浓度,以期为相关部门提出相应的应急对策以及制订相关标准提供依据.1calpuff模拟系统构建安徽和河南两省秸秆焚烧涉及的区域较广,并且气象和地形条件对PM10的扩散有一定的影响.因此,选择长距离中尺度空气质量模型(CALPUFF)比较合适.该模型结合时变的气象资料与复杂的地形条件,能够很好地模拟不同尺度区域污染物的扩散情景,使模拟结果更真实地反映大气污染扩散的实际情况.CALPUFF模拟系统需要气象条件、地理条件和污染物信息等,其中某些污染物信息由火灾排放模型FEPS提供.模拟系统框架见图1.具体包括:(1)根据卫星提供的遥感影像或其他信息确定秸秆焚烧位置、面积等参数,并查找当地的地理和气象条件;(2)将秸秆焚烧当成面源处理,利用FEPS模型确定PM10的烟羽抬升,并根据相关文献和经验值估算PM10的释放速率;(3)利用CALPUFF模型模拟秸秆焚烧排放的ρ(PM10)扩散,并对模拟结果进行分析和可视化.2模拟过程2.1cbd研究时段研究区域是以116°E,33°N为中心,东西500km,南北350km,总面积17.5×104km2,转换为LCC坐标范围为:X(-300~150km),Y(-200~150km),网格间距为5km,垂直方向分为10层.研究时段为2009年6月3日和5日,共2d.2.2数据主要特点秸秆焚烧遥感监测日报是环境保护部卫星环境应用中心基于美国NASA的Terra和Aqua卫星搭载的MODIS传感器每日观测的遥感数据,通过热异常反演方法并结合历史土地利用等数据分析制作而成的.Terra卫星的过境时间大约在每日上午10:30左右,Aqua卫星的过境时间大约在每日下午13:30左右.MODIS波段多,灵敏度更高,监测频率高,每天对同一地方观测次数最多可达4次,因此能够探测到更小的面积(最小面积可达50m2)、更多的火点.由于卫星遥感数据空间分辨率仅为1km,土地利用和行政区划数据存在精度问题,在秸秆焚烧火点的定位会出现一定的偏差,但研究区域范围较大,遥感数据可以满足研究的需要.秸秆焚烧遥感监测日报可提供全国范围内秸秆焚烧火点的经纬度,因此可将各火点的经纬度视为秸秆焚烧火点的中心位置.2.3topo30热导系数的处理CALMET模块需要地形数据和土地利用数据,包括土地类型、海拔高度、地表参数和人为热导系数.其中,地形数据利用子软件TERREL处理GTOPO30文件得到TERREL.DAT数据文件;土地利用数据利用系统子软件CTGPROC处理USGSGlobal-LambertAzimuthal的ASIA数据,可得LU.DAT数据文件.2.4探空气象数据模拟系统中CALMET所需的气象数据包括地面和探空数据.地面气象站数据采用7个地面观测站2009年6月3日和5日逐时气象数据,包括风向、风速、温度、湿度、压力、云量和降水等,各地面站的具体信息见表1.探空气象数据采用9个探空站6月3日和5日每天2次的探空数据,探空站具体信息见表2.CALMET模块可以获得研究区域的时变10层模拟气象场,用以分析高空以及地面的气象场对污染物传输和扩散的影响,其中风速、风向受地形条件影响.LCC坐标下6月3日08:00的地形风场图见图2.2.5秸秆焚烧的面积秸秆焚烧面源数据包括秸秆焚烧火点的空间分布、燃烧持续时间、释放速率和初始垂直扩散系数.秸秆焚烧火点的空间分布包括火点的位置与焚烧面源的面积.火点位置可根据秸秆焚烧遥感监测日报提供的数据确定,焚烧面积信息需根据相关文献以及实际情况考虑.何立明等统计秸秆焚烧像元数和实际年焚烧面积的关系发现,二者间存在很好的线性关系,即焚烧面积=0.6449×火点像元数+27.365.CHOI等模拟美墨边境秸秆焚烧时设定的面积假定为像元的面积大小(221m×323m).SMITH将MODIS监测数据和问卷调查结合起来,估算2个月内平均秸秆焚烧的面积为(75±6)hm2,则每d平均秸秆焚烧为(1.25±0.1)hm2.参考前面3种情况,并为简化模拟,假定火点的面积为1hm2.2.5.2平均风速、风向、温度和稳定度JENKINS等的风洞试验条件:麦秆点火方式为逆风点火,风速为2~3m/s,空气温度为21~27℃,燃料负荷为400~500g/m2.最后试验得到的燃烧蔓延速度为0.2~0.4m/min.研究区域内6月3日和5日的平均风速、风向、温度和稳定度见图3.比较图3的气象条件和JENKINS等的风洞试验条件发现,安徽和河南3日和5日的实际秸秆焚烧条件与风洞的试验条件基本相似,因此可将风洞试验得到的结果作为实际燃烧速率的初步假设.考虑到6月3日与5日大气稳定度较低,燃烧蔓延速率取风洞试验得到的较大蔓延速率(0.4m/min).假定蔓延速率恒定,1hm2的秸秆燃烧持续时间约为4h,假定秸秆焚烧时段为上午08:00—11:00,下午12:00—15:00.2.5.3燃烧持续时间单位的排放因子污染物释放速率计算方法:式中,释放速率单位为g/(m2·s);排放因子单位为g/kg;可燃物负荷量单位为kg/m2;燃烧持续时间单位为s.据报道,PM10的排放因子为5.77g/kg;可燃物负荷量可以用单位面积秸秆量表示,而秸秆量用产量乘以谷草比表示,其中小麦产量为0.525kg/m2,谷草比为1.366,则可燃物负荷量为0.7172kg/m2;燃烧持续时间为4h.计算得到释放速率为2.874×10-4g/(m2·s).2.5.4烟羽升降高度秸秆燃烧过程中释放的热量会产生浮力,浮力会使秸秆焚烧过程产生的大部分气体漂浮并向更广的时空尺度扩散.CALPUFF可以处理浮力面源、一般面源、点源和体源等,模拟秸秆焚烧最合理的应是浮力面源,但因为缺少有效面源垂直速度和燃烧温度等参数,所以考虑由浮力产生的烟羽抬升的一般面源作为替代.其中,烟羽抬升高度的计算采用FEPS模型计算烟羽抬升高度的原理,具体公式:式中,ΔHF-min,i和ΔHF-max,i分别为第i小时FEPS模型中最小与最大经验热气抬升高度,m;ARi为第i-1小时与ih内的燃烧面积变化率,hm2/h;PBot为最大有效烟羽底部抬升高度,缺省值为4000m;PTop为最大有效烟羽抬升高度,缺省值为8000m;ARBot与ARTop为基准变化率,缺省值均为161.876hm2/h.1hm2的秸秆4h烧完,平均燃烧面积变化率为0.25m2/h,带入式(1),(2)得到ΔHF-min,i为12.344m,ΔHF-max,i为6.172m.初始垂直扩散系数(σx)=烟羽抬升高度/3.5,最大初始垂直扩散系数σmax=3.527m,最小初始垂直扩散系数σmin=1.764m.2.6param生存条件重新动力系统的重新动态分析CALPUFF模型默认的可以模拟的最大面源数是200个,而6月3日要模拟的火点数是241个,6月5日要模拟的火点数为423个.为满足模拟需求,对源程序进行重新编译.首先将paramsl.puf中最大面源数限值改成500个,然后利用LFFortranProfessionalv7.2重新编译源程序,最后在CALPUFF主页面上执行新的参数文件.3结果与讨论3.1最大烟羽升降下逐时最CALPUFF模拟6月3日与5日最大、最小烟羽抬升下逐时最大ρ(PM10)见表3.由表3可知,3日和5日00:00—07:00因没有秸秆焚烧现象,在不考虑背景值的情况下ρ(PM10)为0;3日08:00—23:00最大烟羽抬升下逐时最大ρ(PM10)范围为4.9961~3040.7μg/m3;最小烟羽抬升下的ρ(PM10)范围为5.0078~6069.2μg/m3.5日08:00—23:00最大烟羽抬升下逐时的最大ρ(PM10)范围为4.285~2562μg/m3;最小烟羽抬升下的ρ(PM10)范围为4.2849~5122.4μg/m3.对比3日和5日最大和最小烟羽抬升下的ρ(PM10),在秸秆燃烧阶段最小烟羽抬升下ρ(PM10)是最大烟羽抬升下的2倍,而在其他时段ρ(PM10)不存在明显的比例关系.3.2初始烟羽升降高度3日和5日最大、最小烟羽抬升下ρ(PM10)日均值以及逐时ρ(PM10)空间分布见图4~7,其中最大ρ(PM10)标在图4~7中,单位为μg/m3.由图4~7可知,ρ(PM10)较高的区域火点比较密集,并且最大值出现在火点密集的区域.扩散过程受气象条件影响明显,受风速及风向的影响尤其大.将图4与图5,图6与图7对比可知,初始烟羽抬升高度对扩散有明显的影响,当烟羽抬升高度较小时,扩散浓度相对较大.对图4~7各子图进行分析可知,秸秆焚烧时段逐时ρ(PM10)要远大于日均值和非燃烧时段23:00的值.模拟区域内3日火点数为241个,5日为423个.根据模拟结果,6月3日最大烟羽抬升下PM10总排放量为3716.4μg/m3,最小烟羽抬升下PM10总排放量为6551.2μg/m3;6月5日最大烟羽抬升下PM10总排放量为9682.8μg/m3,最小烟羽抬升下PM10总排放量为18169.1μg/m3.可见,总排放量不仅与烟羽抬升高度有关,而且还与火点数有关,火点越多,总的ρ(PM10)越大.3日最小、最大烟羽抬升下的最大日均ρ(PM10)分别为323.9和164μg/m3,5日分别为360.1和180.77μg/m3,可见部分区域日均ρ(PM10)大于国家二级标准的150μg/m3甚至三级标准的350μg/m3,说明秸秆焚烧现象对空气质量有比较显著的影响.3.3秸秆焚烧与空气总有关见表36月3日和5日郑州、合肥和平顶山的首要污染物都为可吸入颗粒物,具体API数值和推算得到的ρ(PM10)见表4.由表4可知,各城市5日的空气质量比3日的差,ρ(PM10)也相应变大,与秸秆焚烧量有关,说明秸秆焚烧对空气质量有一定的影响.3.4模型内部的不确定性模型只模拟了安徽和河南6月3日和5日当天4h秸秆燃烧对周围环境的影响,如果考虑麦收后期连续几天或十几天的燃烧,则监测到的火点要比实际情况少,以及考虑二次颗粒的形成等因素,那么秸秆焚烧对周围环境的影响会更加严重.在模拟过程中一些输入参数只是简单地估算或设定的缺省值,致使输入参数不够精确,所以模型具有一定的不确定性.4秸秆焚烧持续快速发展a.秸秆焚烧排放的ρ(PM10)