楊懷榮，劉 茂，劉付衍華

南開(kāi)大學(xué)城市公共安全研究中心，天津 300071

利用CALPUFF對(duì)安徽和河南秸稈焚燒的模擬與研究

楊懷榮，劉 茂*，劉付衍華

南開(kāi)大學(xué)城市公共安全研究中心，天津 300071

秸稈焚燒會(huì)產(chǎn)生大量的顆粒物(PM)、氮氧化物、有機(jī)碳、苯以及多環(huán)芳烴等污染物，不僅影響空氣質(zhì)量，危害人體健康，而且大大降低能見(jiàn)度，對(duì)交通運(yùn)輸構(gòu)成威脅.針對(duì)安徽和河南2009年6月嚴(yán)重的秸稈焚燒現(xiàn)象，對(duì)CALPUFF模擬系統(tǒng)和FEPS模型進(jìn)行重新編譯與整合，對(duì)空氣動(dòng)力學(xué)直徑小于10μm的顆粒物(PM10)進(jìn)行擴(kuò)散模擬，得到逐時(shí)ρ(PM10)的煙羽擴(kuò)散，并對(duì)模擬結(jié)果進(jìn)行分析.結(jié)果表明，秸稈焚燒過(guò)程中焚燒點(diǎn)附近的ρ(PM10)較大，研究區(qū)域內(nèi)部分區(qū)域的日均ρ(PM10)大于我國(guó)二級(jí)標(biāo)準(zhǔn)(150μg/m3)甚至三級(jí)標(biāo)準(zhǔn)(350μg/m3).如果考慮二次粒子，其影響程度會(huì)更加嚴(yán)重.

秸稈焚燒;顆粒物;煙羽擴(kuò)散;CALPUFF模擬系統(tǒng);FEPS模型

Abstract:Straw crop burning em its a number of pollutants，including particulate matter(PM)，nitrogen oxides，organic carbons (OC)，benzene，polycyclic aromatic hydrocarbons，etc.These pollutants not only impact air quality and are harm ful to human health，but also greatly reduce visibility and thus pose a threat to transportation safety.In view of the serious straw crop burning phenomenon in Anhui and Henan Provinces in June，2009，the CALPUFF modeling system and FEPS(Fire Emission Production Simulator)model were recompiled and integrated to simulate dispersion of PM with aerodynamic diameter less than 10μm(PM10)from the straw crop burning.Also，hourly PM10mass concentration p lume diffusions were obtained.Simulation results show that high PM10mass concentrations are observed near the fire sources，and that daily average PM10mass concentrations in some areas are higher than the national daily average standard levelⅡ (150μg/m3)and even levelⅢ (350μg/m3).If the formation of secondary particles were also considered，the impacts could be even greater.

Keywords:straw crop burning;particulate matter;plume dispersion;CALPUFF modeling system;FEPSmodel

近年來(lái)我國(guó)對(duì)農(nóng)作物秸稈焚燒采取了多種措施以降低其規(guī)模，但焚燒現(xiàn)象依然難以遏制.農(nóng)作物秸稈燃燒時(shí)會(huì)產(chǎn)生大量的 PM，有機(jī)碳，氮氧化物，苯及多環(huán)芳烴等有害氣體，不僅危害人體健康，造成環(huán)境污染［1］，而且產(chǎn)生大量煙霧，大大降低能見(jiàn)度，威脅交通運(yùn)輸安全.2007年5月華北和江淮流域的糧食主產(chǎn)區(qū)出現(xiàn)大規(guī)模秸稈焚燒現(xiàn)象，使北京等周邊城市的空氣質(zhì)量明顯下降［2］.1997年5月14— 18日，成都雙流機(jī)場(chǎng)由于焚燒秸稈產(chǎn)生大量煙霧造成10個(gè)航班無(wú)法降落，22架進(jìn)出港航班延誤，滯留旅客1 000多人［3］.

目前國(guó)內(nèi)秸稈焚燒的研究主要有遙感對(duì)秸稈焚燒的監(jiān)測(cè)、排放因子的試驗(yàn)研究以及排放清單的估算等，但很少對(duì)污染物的擴(kuò)散進(jìn)行模擬.祝斌等［4］對(duì)農(nóng)作物秸稈燃燒時(shí) PM2.5排放因子進(jìn)行了試驗(yàn)研究，并分別對(duì)水稻、小麥和玉米等農(nóng)作物秸稈在明火和悶火狀態(tài)下的燃燒情況進(jìn)行了試驗(yàn)，測(cè)定了燃燒過(guò)程中PM2.5的排放因子.王書(shū)肖等［5］采用排放因子法，建立了我國(guó)秸稈露天焚燒一次大氣污染物的排放清單，并分析了其時(shí)間和空間的分布特征.厲青等［6］利用衛(wèi)星遙感技術(shù)監(jiān)測(cè)了2007年6月全國(guó)秸稈焚燒狀況，并以某市為例，結(jié)合氣象資料分析秸稈焚燒對(duì)環(huán)境空氣質(zhì)量的影響.方萌等［7］介紹了“3S”技術(shù)在農(nóng)作物秸稈焚燒監(jiān)測(cè)中應(yīng)用的一般原理和方法，同時(shí)還介紹了國(guó)家衛(wèi)星氣象中心近年來(lái)采用“3S”技術(shù)監(jiān)測(cè)農(nóng)作物秸稈焚燒的概況.國(guó)外對(duì)于秸稈焚燒的研究除了有對(duì)排放因子的計(jì)算、遙感的監(jiān)測(cè)外，還有對(duì)焚燒過(guò)程中產(chǎn)生污染物的擴(kuò)散進(jìn)行的計(jì)算機(jī)模擬.JAIN等［8］對(duì)太平洋海岸西北部秸稈焚燒的煙氣擴(kuò)散預(yù)測(cè)系統(tǒng)(ClearSky)進(jìn)行了介紹，ClearSky是數(shù)值煙氣擴(kuò)散預(yù)測(cè)系統(tǒng)，通過(guò)確定秸稈焚燒情景，利用 MM5中尺度氣象預(yù)測(cè)模型和CALPUFF模擬系統(tǒng)，模擬秸稈焚燒產(chǎn)生的 PM2.5濃度擴(kuò)散情況.CHOI等［9］利用 CALPUFF/CALMET/ MM5模擬系統(tǒng)，對(duì)美國(guó)墨西哥邊境秸稈焚燒排放的PM10濃度進(jìn)行了擴(kuò)散模擬，并與PM10的相關(guān)環(huán)境標(biāo)準(zhǔn)進(jìn)行了對(duì)比.

據(jù)秸稈焚燒遙感監(jiān)測(cè)日?qǐng)?bào)發(fā)布的2009年全國(guó)秸稈焚燒分布遙感監(jiān)測(cè)結(jié)果發(fā)現(xiàn)，安徽和河南兩省6月中上旬的秸稈焚燒現(xiàn)象均十分嚴(yán)重.針對(duì)秸稈焚燒嚴(yán)重的現(xiàn)象，筆者借鑒國(guó)內(nèi)外模擬顆粒物擴(kuò)散的經(jīng)驗(yàn)，整合了CALPUFF模型與FEPS模型，并對(duì)CALPUFF模型的部分程序進(jìn)行了重新編譯，模擬了安徽和河南兩省2009年6月3日與5日秸稈焚燒產(chǎn)生的ρ(PM10)擴(kuò)散情況，得到了PM10的煙羽擴(kuò)散濃度，以期為相關(guān)部門(mén)提出相應(yīng)的應(yīng)急對(duì)策以及制訂相關(guān)標(biāo)準(zhǔn)提供依據(jù).

1 方法介紹

安徽和河南兩省秸稈焚燒涉及的區(qū)域較廣，并且氣象和地形條件對(duì) PM10的擴(kuò)散有一定的影響.因此，選擇長(zhǎng)距離中尺度空氣質(zhì)量模型(CALPUFF)比較合適.該模型結(jié)合時(shí)變的氣象資料與復(fù)雜的地形條件，能夠很好地模擬不同尺度區(qū)域污染物的擴(kuò)散情景，使模擬結(jié)果更真實(shí)地反映大氣污染擴(kuò)散的實(shí)際情況［10-11］.CALPUFF模擬系統(tǒng)需要?dú)庀髼l件、地理?xiàng)l件和污染物信息等，其中某些污染物信息由火災(zāi)排放模型 FEPS提供［12］.模擬系統(tǒng)框架見(jiàn)圖1.具體包括:①根據(jù)衛(wèi)星提供的遙感影像或其他信息確定秸稈焚燒位置、面積等參數(shù)，并查找當(dāng)?shù)氐牡乩砗蜌庀髼l件;②將秸稈焚燒當(dāng)成面源處理，利用FEPS模型確定PM10的煙羽抬升，并根據(jù)相關(guān)文獻(xiàn)和經(jīng)驗(yàn)值估算PM10的釋放速率;③利用CALPUFF模型模擬秸稈焚燒排放的ρ(PM10)擴(kuò)散，并對(duì)模擬結(jié)果進(jìn)行分析和可視化.

圖1 CALPUFF模擬系統(tǒng)與FEPS模型整合的模擬系統(tǒng)框架Fig.1 Framework of integrating CALPUFF system and FEPSmodel

2 模擬過(guò)程

2.1 研究區(qū)域與時(shí)段

研究區(qū)域是以116°E，33°N為中心，東西500 km，南北 350 km，總面積 17.5×104km2，轉(zhuǎn)換為L(zhǎng)CC坐標(biāo)范圍為:X(－300～150 km)，Y(－200～150 km)，網(wǎng)格間距為5 km，垂直方向分為10層.研究時(shí)段為2009年6月3日和5日，共2 d.

2.2 秸稈焚燒火點(diǎn)的確定

秸稈焚燒遙感監(jiān)測(cè)日?qǐng)?bào)是環(huán)境保護(hù)部衛(wèi)星環(huán)境應(yīng)用中心基于美國(guó)NASA的Terra和Aqua衛(wèi)星搭載的MODIS傳感器每日觀測(cè)的遙感數(shù)據(jù)，通過(guò)熱異常反演方法并結(jié)合歷史土地利用等數(shù)據(jù)分析制作而成的.Terra衛(wèi)星的過(guò)境時(shí)間大約在每日上午10:30左右，Aqua衛(wèi)星的過(guò)境時(shí)間大約在每日下午13:30左右.MODIS波段多，靈敏度更高，監(jiān)測(cè)頻率高，每天對(duì)同一地方觀測(cè)次數(shù)最多可達(dá)4次，因此能夠探測(cè)到更小的面積(最小面積可達(dá)50 m2)、更多的火點(diǎn).由于衛(wèi)星遙感數(shù)據(jù)空間分辨率僅為1 km，土地利用和行政區(qū)劃數(shù)據(jù)存在精度問(wèn)題，在秸稈焚燒火點(diǎn)的定位會(huì)出現(xiàn)一定的偏差，但研究區(qū)域范圍較大，遙感數(shù)據(jù)可以滿足研究的需要.秸稈焚燒遙感監(jiān)測(cè)日?qǐng)?bào)可提供全國(guó)范圍內(nèi)秸稈焚燒火點(diǎn)的經(jīng)緯度，因此可將各火點(diǎn)的經(jīng)緯度視為秸稈焚燒火點(diǎn)的中心位置.

2.3 地理數(shù)據(jù)

CALMET模塊需要地形數(shù)據(jù)和土地利用數(shù)據(jù)，包括土地類(lèi)型、海拔高度、地表參數(shù)和人為熱導(dǎo)系數(shù).其中，地形數(shù)據(jù)利用子軟件 TERREL處理GTOPO30文件得到TERREL.DAT數(shù)據(jù)文件;土地利用數(shù)據(jù)利用系統(tǒng)子軟件 CTGPROC處理 USGSGlobal-Lambert Azimuthal的 ASIA數(shù)據(jù)，可得 LU. DAT數(shù)據(jù)文件.

2.4 氣象資料

模擬系統(tǒng)中CALMET所需的氣象數(shù)據(jù)包括地面和探空數(shù)據(jù).地面氣象站數(shù)據(jù)采用7個(gè)地面觀測(cè)站2009年6月3日和5日逐時(shí)氣象數(shù)據(jù)，包括風(fēng)向、風(fēng)速、溫度、濕度、壓力、云量和降水等，各地面站的具體信息見(jiàn)表1.探空氣象數(shù)據(jù)采用9個(gè)探空站6月3日和5日每天2次的探空數(shù)據(jù)，探空站具體信息見(jiàn)表2.CALMET模塊可以獲得研究區(qū)域的時(shí)變10層模擬氣象場(chǎng)，用以分析高空以及地面的氣象場(chǎng)對(duì)污染物傳輸和擴(kuò)散的影響，其中風(fēng)速、風(fēng)向受地形條件影響.LCC坐標(biāo)下6月3日08:00的地形風(fēng)場(chǎng)圖見(jiàn)圖2.

表1 地面監(jiān)測(cè)站具體信息Table 1 Basic information about the surface stations

表2 探空站具體信息Table 2 Basic information about the upper air stations

圖2 2009年6月3日08:00地形風(fēng)場(chǎng)圖Fig.2 Simulated surface winds over the CALPUFF domain(08:00 July 3rd，2009)

2.5 秸稈焚燒面源數(shù)據(jù)

秸稈焚燒面源數(shù)據(jù)包括秸稈焚燒火點(diǎn)的空間分布、燃燒持續(xù)時(shí)間、釋放速率和初始垂直擴(kuò)散系數(shù).

2.5.1 秸稈焚燒火點(diǎn)的空間分布

秸稈焚燒火點(diǎn)的空間分布包括火點(diǎn)的位置與焚燒面源的面積.火點(diǎn)位置可根據(jù)秸稈焚燒遙感監(jiān)測(cè)日?qǐng)?bào)提供的數(shù)據(jù)確定，焚燒面積信息需根據(jù)相關(guān)文獻(xiàn)以及實(shí)際情況考慮.何立明等［13］統(tǒng)計(jì)秸稈焚燒像元數(shù)和實(shí)際年焚燒面積的關(guān)系發(fā)現(xiàn)，二者間存在很好的線性關(guān)系，即焚燒面積=0.644 9×火點(diǎn)像元數(shù)+27.365.CHOI等［9］模擬美墨邊境秸稈焚燒時(shí)設(shè)定的面積假定為像元的面積大小(221 m×323 m).SMITH［14］將MODIS監(jiān)測(cè)數(shù)據(jù)和問(wèn)卷調(diào)查結(jié)合起來(lái)，估算2個(gè)月內(nèi)平均秸稈焚燒的面積為(75± 6)hm2，則每 d平均秸稈焚燒為(1.25±0.1)hm2.參考前面3種情況，并為簡(jiǎn)化模擬，假定火點(diǎn)的面積為1 hm2.

2.5.2 秸稈燃燒持續(xù)時(shí)間

JENKINS等［15］的風(fēng)洞試驗(yàn)條件:麥稈點(diǎn)火方式為逆風(fēng)點(diǎn)火，風(fēng)速為 2～3 m/s，空氣溫度為 21～27℃，燃料負(fù)荷為400～500 g/m2.最后試驗(yàn)得到的燃燒蔓延速度為0.2～0.4 m/m in.研究區(qū)域內(nèi)6月3日和5日的平均風(fēng)速、風(fēng)向、溫度和穩(wěn)定度見(jiàn)圖3.比較圖3的氣象條件和JENKINS等［15］的風(fēng)洞試驗(yàn)條件發(fā)現(xiàn)，安徽和河南3日和5日的實(shí)際秸稈焚燒條件與風(fēng)洞的試驗(yàn)條件基本相似，因此可將風(fēng)洞試驗(yàn)得到的結(jié)果作為實(shí)際燃燒速率的初步假設(shè).考慮到6月3日與5日大氣穩(wěn)定度較低，燃燒蔓延速率取風(fēng)洞試驗(yàn)得到的較大蔓延速率(0.4 m/min).假定蔓延速率恒定，1 hm2的秸稈燃燒持續(xù)時(shí)間約為4 h，假定秸稈焚燒時(shí)段為上午08:00—11:00，下午12:00—15:00.

圖3 2009年6月3日研究區(qū)域內(nèi)平均風(fēng)向、風(fēng)速、溫度、穩(wěn)定度隨時(shí)間的變化Fig.3 Average wind direction(a)，average wind speed(b)，average temperature(c)and average atmospheric stability category(d)for the entire burning source area on July 3rd，2009

2.5.3 污染物釋放速率

污染物釋放速率計(jì)算方法:

式中，釋放速率單位為g/(m2·s);排放因子單位為g/kg;可燃物負(fù)荷量單位為kg/m2;燃燒持續(xù)時(shí)間單位為 s.

據(jù)報(bào)道，PM10的排放因子為5.77 g/kg［16］;可燃物負(fù)荷量可以用單位面積秸稈量表示，而秸稈量用產(chǎn)量乘以谷草比表示，其中小麥產(chǎn)量為 0.525 kg/m2，谷草比為 1.366［16］，則可燃物負(fù)荷量為0.717 2 kg/m2;燃燒持續(xù)時(shí)間為4 h.計(jì)算得到釋放速率為2.874×10－4g/(m2·s).

2.5.4 初始垂直擴(kuò)散系數(shù)

秸稈燃燒過(guò)程中釋放的熱量會(huì)產(chǎn)生浮力，浮力會(huì)使秸稈焚燒過(guò)程產(chǎn)生的大部分氣體漂浮并向更廣的時(shí)空尺度擴(kuò)散.CALPUFF可以處理浮力面源、一般面源、點(diǎn)源和體源等，模擬秸稈焚燒最合理的應(yīng)是浮力面源，但因?yàn)槿鄙儆行嬖创怪彼俣群腿紵郎囟鹊葏?shù)，所以考慮由浮力產(chǎn)生的煙羽抬升的一般面源作為替代.其中，煙羽抬升高度的計(jì)算采用FEPS模型計(jì)算煙羽抬升高度的原理，具體公式:

式中，ΔHF-min，i和 ΔHF-max，i分別為第 i小時(shí) FEPS模型中最小與最大經(jīng)驗(yàn)熱氣抬升高度，m;ARi為第i－1小時(shí)與i h內(nèi)的燃燒面積變化率，hm2/h;PBot為最大有效煙羽底部抬升高度，缺省值為4 000 m;PTop為最大有效煙羽抬升高度，缺省值為8 000 m;ARBot與ARTop為基準(zhǔn)變化率，缺省值均為161.876 hm2/h.

1 hm2的秸稈4 h燒完，平均燃燒面積變化率為0.25 m2/h，帶入式(1)，(2)得到 ΔHF-min，i為12.344m，ΔHF-max，i為6.172 m.初始垂直擴(kuò)散系數(shù)(σx)=煙羽抬升高度/3.5［11］，最大初始垂直擴(kuò)散系數(shù)σmax= 3.527 m，最小初始垂直擴(kuò)散系數(shù)σmin=1.764 m.

2.6 CALPUFF模型系統(tǒng)的重新編譯

CALPUFF模型默認(rèn)的可以模擬的最大面源數(shù)是200個(gè)，而6月3日要模擬的火點(diǎn)數(shù)是241個(gè)，6月5日要模擬的火點(diǎn)數(shù)為423個(gè).為滿足模擬需求，對(duì)源程序進(jìn)行重新編譯.首先將 paramsl.puf中最大面源數(shù)限值改成 500個(gè)，然后利用 LF Fortran Professional v7.2重新編譯源程序，最后在CALPUFF主頁(yè)面上執(zhí)行新的參數(shù)文件.

3 結(jié)果與討論

3.1 CALPUFF模擬的逐時(shí)ρ(PM10)

CALPUFF模擬6月3日與5日最大、最小煙羽抬升下逐時(shí)最大ρ(PM10)見(jiàn)表3.由表3可知，3日和5日00:00—07:00因沒(méi)有秸稈焚燒現(xiàn)象，在不考慮背景值的情況下ρ(PM10)為0;3日08:00—23:00最大煙羽抬升下逐時(shí)最大ρ(PM10)范圍為4.996 1～3 040.7μg/m3;最小煙羽抬升下的ρ(PM10)范圍為5.007 8～6 069.2μg/m3.5日08:00—23:00最大煙羽抬升下逐時(shí)的最大ρ(PM10)范圍為 4.285～2 562μg/m3;最小煙羽抬升下的ρ(PM10)范圍為4.284 9～5 122.4μg/m3.對(duì)比3日和5日最大和最小煙羽抬升下的ρ(PM10)，在秸稈燃燒階段最小煙羽抬升下ρ(PM10)是最大煙羽抬升下的2倍，而在其他時(shí)段ρ(PM10)不存在明顯的比例關(guān)系.

表3 最大、最小煙羽抬升下逐時(shí)最大ρ(PM10)Table 3 The maximum hourly PM10mass concentration with different plume rise μg/m3

3.2 CALPUFF模擬ρ(PM10)的時(shí)空分布

3日和5日最大、最小煙羽抬升下ρ(PM10)日均值以及逐時(shí)ρ(PM10)空間分布見(jiàn)圖4～7，其中最大ρ(PM10)標(biāo)在圖4～7中，單位為 μg/m3.由圖4～7可知，ρ(PM10)較高的區(qū)域火點(diǎn)比較密集，并且最大值出現(xiàn)在火點(diǎn)密集的區(qū)域.擴(kuò)散過(guò)程受氣象條件影響明顯，受風(fēng)速及風(fēng)向的影響尤其大.將圖4與圖5，圖6與圖7對(duì)比可知，初始煙羽抬升高度對(duì)擴(kuò)散有明顯的影響，當(dāng)煙羽抬升高度較小時(shí)，擴(kuò)散濃度相對(duì)較大.對(duì)圖4～7各子圖進(jìn)行分析可知，秸稈焚燒時(shí)段逐時(shí)ρ(PM10)要遠(yuǎn)大于日均值和非燃燒時(shí)段23:00的值.

模擬區(qū)域內(nèi)3日火點(diǎn)數(shù)為241個(gè)，5日為423個(gè).根據(jù)模擬結(jié)果，6月3日最大煙羽抬升下 PM10總排放量為3 716.4μg/m3，最小煙羽抬升下 PM10總排放量為6 551.2μg/m3;6月5日最大煙羽抬升下PM10總排放量為9 682.8μg/m3，最小煙羽抬升下PM10總排放量為18 169.1μg/m3.可見(jiàn)，總排放量不僅與煙羽抬升高度有關(guān)，而且還與火點(diǎn)數(shù)有關(guān)，火點(diǎn)越多，總的ρ(PM10)越大.

3日最小、最大煙羽抬升下的最大日均ρ(PM10)分別為323.9和164μg/m3，5日分別為360.1和180.77 μg/m3，可見(jiàn)部分區(qū)域日均ρ(PM10)大于國(guó)家二級(jí)標(biāo)準(zhǔn)的150μg/m3甚至三級(jí)標(biāo)準(zhǔn)的350μg/m3，說(shuō)明秸稈焚燒現(xiàn)象對(duì)空氣質(zhì)量有比較顯著的影響.

3.3 ρ(PM10)與空氣質(zhì)量

6月3日和5日鄭州、合肥和平頂山的首要污染物都為可吸入顆粒物，具體 API數(shù)值和推算得到的ρ(PM10)見(jiàn)表4.由表4可知，各城市5日的空氣質(zhì)量比3日的差，ρ(PM10)也相應(yīng)變大，與秸稈焚燒量有關(guān)，說(shuō)明秸稈焚燒對(duì)空氣質(zhì)量有一定的影響.

圖4 6月3日最小煙羽抬升下日均ρ(PM10)(μg/m3)及11:00，15:00和23:00逐時(shí)平均ρ(PM10)(μg/m3)Fig.4 PM10mass concentrations(μg/m3)in m inimum plume rise on June 3rd，2009

圖5 6月3日最大煙羽抬升下日均ρ(PM10)(μg/m3)及11:00，15:00和23:00逐時(shí)平均ρ(PM10)(μg/m3)Fig.5 PM10mass concentration(μg/m3)in maximum plume rise on June 3rd，2009

圖6 6月5日最小煙羽抬升下日均ρ(PM10)(μg/m3)及12:00，15:00和23:00逐時(shí)平均ρ(PM10)(μg/m3)Fig.6 PM10mass concentration(μg/m3)in minimum plume rise on June 5th，2009

圖7 6月5日最大煙羽抬升下日均ρ(PM10)(μg/m3)及12:00，15:00和23:00逐時(shí)平均ρ(PM10)(μg/m3)Fig.7 PM10mass concentration(μg/m3)in maximum plume rise on June 5th，2009

表4 各城市API指數(shù)與ρ(PM10)Table 4 The air pollutant index and mass concentrations of PM10in three cities

3.4 模型局限性

模型只模擬了安徽和河南6月3日和5日當(dāng)天4 h秸稈燃燒對(duì)周?chē)h(huán)境的影響，如果考慮麥?zhǔn)蘸笃谶B續(xù)幾天或十幾天的燃燒，則監(jiān)測(cè)到的火點(diǎn)要比實(shí)際情況少，以及考慮二次顆粒的形成等因素，那么秸稈焚燒對(duì)周?chē)h(huán)境的影響會(huì)更加嚴(yán)重.

在模擬過(guò)程中一些輸入?yún)?shù)只是簡(jiǎn)單地估算或設(shè)定的缺省值，致使輸入?yún)?shù)不夠精確，所以模型具有一定的不確定性.

4 結(jié)論

a.秸稈焚燒排放的ρ(PM10)受煙羽抬升高度、火點(diǎn)數(shù)以及火點(diǎn)密集程度等因素的影響.煙羽抬升高度越高，ρ(PM10)越低.火點(diǎn)數(shù)越多，總ρ(PM10)越大.火點(diǎn)密集的地方ρ(PM10)較高，最大值出現(xiàn)在火點(diǎn)密集的區(qū)域.

b.在模擬區(qū)域內(nèi)，焚燒階段最小煙羽抬升下逐時(shí)ρ(PM10)是最大煙羽抬升下的2倍，而在其他非焚燒時(shí)段最小和最大煙羽抬升下的逐時(shí)ρ(PM10)不存在明顯的比例關(guān)系.

c.部分區(qū)域日均ρ(PM10)大于國(guó)家二級(jí)標(biāo)準(zhǔn)的150μg/m3甚至三級(jí)標(biāo)準(zhǔn)的350μg/m3，說(shuō)明秸稈焚燒現(xiàn)象對(duì)空氣質(zhì)量有顯著的影響.如果考慮二次粒子和其他條件，影響程度會(huì)更加嚴(yán)重.

［1］MIDDLEBROOK A，TURNER J，SOLOMAN P A.Preface:special issue of atmospheric environment for particulate matter: atmospheric sciences，exposure and the fourth colloquium on PM and human health［J］.Atmos Environ，2004，38(31):5179-5181.

［2］王子峰，陳良富，顧行發(fā).基于 MODIS數(shù)據(jù)的華北地區(qū)秸稈焚燒監(jiān)測(cè)［J］.遙感技術(shù)與應(yīng)用，2008，23(6):611-617.

［3］徐玉宏.我國(guó)秸稈焚燒污染與防治對(duì)策［J］.環(huán)境與可持續(xù)發(fā)展，2009(3):21-24.

［4］祝斌，朱先磊，張?jiān)獎(jiǎng)?，?農(nóng)作物秸稈燃燒 PM2.5排放因子的研究［J］.環(huán)境科學(xué)研究，2005，18(2):29-33.

［5］王書(shū)肖，張楚瑩.中國(guó)秸稈露天焚燒大氣污染物排放時(shí)空分布［J］.中國(guó)科技論文在線，2008，3(5):329-333.

［6］厲青，張麗娟，吳傳慶，等.基于衛(wèi)星遙感的秸稈焚燒監(jiān)測(cè)及對(duì)空氣質(zhì)量影響分析［J］.生態(tài)與農(nóng)業(yè)環(huán)境學(xué)報(bào)，2009，25 (1):32-37.

［7］方萌，張鵬，徐喆.“3S”技術(shù)在農(nóng)作物秸稈焚燒監(jiān)測(cè)中的應(yīng)用［J］.國(guó)土資源遙感，2006，69(3):1-4.

［8］JAIN R，VAUGHAN J，HEITKAMP K，et al.Development of the ClearSky smoke dispersion forecast system for agricultural field burning in the Pacific Northwest［J］.Atmos Environ，2007，41 (32):6745-6761.

［9］CHOIY J，F(xiàn)ERNANDO H JS.Simulation of smoke plumes from agricultural burns:application to the San Luis/Rio Colorado airshed along the U.S./Mexico border［J］.Sci Total Environ，2007，388(1/2/3):270-289.

［10］SCIRE JS，ROBE F R，F(xiàn)ERNAU M E，et al.A user′s guide for the CALMET meteorological model［M］.5 th ed.Concord，MA: Earth Tech Inc，2000.

［11］SCIRE JS，STRIMATIS D G，YAMARTINO R J.A user′s guide for the CALPUFF dispersion model［M］.5 th ed.Concord，MA: Earth Tech Inc，2000.

［12］ANDERSON G K，SANDBERG D V，NORHEIM R A.Fire emission production simulator(FEPS)User′s Guide［M/OL］.［2009-12-26］.http://www.fs.fed.us/pnw/fera/publications/ fulltext/FEPS_User_Guide.pdf.

［13］何立明，王文杰，王橋，等.中國(guó)秸稈焚燒的遙感監(jiān)測(cè)與分析［J］.中國(guó)環(huán)境監(jiān)測(cè)，2007，23(1):42-50.

［14］SMITH R，ADAMS M，MAIER S，et al.Estimating the area of stubble burning from the number of active fires detected by satellite［J］.Remote Sensing of Environment，2007，109(1):95-106.

［15］JENKINS B M，TURN S Q，W ILLIAMS R B.Atmospheric pollutant emission factors from open burning of agricultural and forest biomass by wind tunnel simulations［M］.Sacramento: California State Air Resources Board，1996.

［16］曹?chē)?guó)良，張小曳，王丹，等.秸稈露天焚燒排放的 TSP等污染物清單［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2005，24(4):800-804.

Research and Sim u lation o f Straw Crop Burning in Anhuiand Henan Provinces Using CALPUFF

YANG Huai-rong，LIU Mao，LIUFU Yan-hua
Research Center of Urban Public Safety，Nankai University，Tianjin 300071，China

X513

A

1001－6929(2010)11－1368－08

2010－03－12

2010－06－30

國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目(70833003);國(guó)家“十一五”科技支撐計(jì)劃項(xiàng)目(200603746006)

楊 懷 榮 (1986 －)，女，黑 龍 江 牡 丹 江 人，mapleyb@126.com.

*責(zé)任作者，劉茂(1945－)，男，遼寧撫順人，教授，博導(dǎo)，主要從事城市公共安全研究，lmsafety@nankai.edu.cn