錢躍東,王勤耕(南京大學(xué)環(huán)境學(xué)院,污染控制與資源化研究國家重點實驗室,江蘇 南京 210046)

針對大尺度區(qū)域的大氣環(huán)境容量綜合估算方法

錢躍東,王勤耕*(南京大學(xué)環(huán)境學(xué)院,污染控制與資源化研究國家重點實驗室,江蘇 南京 210046)

針對大尺度區(qū)域的特點,提出了將箱模型法、模擬法與線性規(guī)劃法相結(jié)合的大氣環(huán)境容量估算方法.其基本思路是將大尺度區(qū)域劃分為一系列小區(qū)域,將真實污染源簡化為源單元,基于線性規(guī)劃法尋求區(qū)域最大容量及其優(yōu)化配置方案.采用箱模型以保證源單元尺度的環(huán)境質(zhì)量要求,采用模擬法以考慮二次污染對前體物環(huán)境容量的制約,將箱模型容量與模擬法容量作為規(guī)劃法的約束條件.新方法使3種傳統(tǒng)方法相互取長補(bǔ)短,其顯著特征表現(xiàn)為:兼顧了污染源的局地性和區(qū)域性的影響;兼顧了污染物的一次影響和二次影響;大氣環(huán)境容量的規(guī)劃布局具有區(qū)域優(yōu)化特征.應(yīng)用新方法對我國東南沿海某區(qū)域進(jìn)行了大氣環(huán)境容量的評估試驗,結(jié)果表明該方法合理可行.

大尺度區(qū)域；大氣環(huán)境容量；線性規(guī)劃法；箱模型；模擬法；環(huán)境規(guī)劃

近年來,針對大尺度區(qū)域的發(fā)展戰(zhàn)略規(guī)劃及環(huán)境評價日益得到重視.大氣環(huán)境容量(AEC)評估是開展上述工作的重要基礎(chǔ).特別是對于區(qū)域產(chǎn)業(yè)布局、污染物排放總量控制及有關(guān)環(huán)保政策的制定等工作來說,大氣環(huán)境容量往往是不可缺少的決策依據(jù).相對于中小尺度的城市或工業(yè)開發(fā)區(qū),大尺度區(qū)域具有以下特征:空間尺度大,一般在幾百 km以上,下墊面和氣象條件復(fù)雜;排放源多樣,分布不均,且數(shù)量眾多;大氣污染具有區(qū)域性和復(fù)合性,特別是城市群區(qū)域,二次污染過程不可忽略[1-2].上述特征對大氣環(huán)境容量評估提出了特別的要求和挑戰(zhàn).

目前常用的大氣環(huán)境容量估算方法主要有箱模型法(或 A-P值法)、模擬法、線性規(guī)劃法等[3-5].它們各有優(yōu)缺點,均不能很好地滿足大尺度區(qū)域的應(yīng)用需求.如箱模型法在城市、工業(yè)集聚區(qū)等地區(qū)得到了較多應(yīng)用[6-9],主要適用于尺度較小的區(qū)域.模擬法在大氣環(huán)境容量評估工作中得到了廣泛的應(yīng)用[10-11],相比箱模型法,其輸入要求高、計算量大.另外,在容量的區(qū)域配置方面,模擬法一般采用等比例或平方比例削減技術(shù),不具有區(qū)域優(yōu)化特性[4,10].線性規(guī)劃法可以像模擬法一樣較細(xì)致地反映“排放源-受體”的響應(yīng)關(guān)系,同時可以在區(qū)域上對環(huán)境容量進(jìn)行優(yōu)化配置,因此得到了十分廣泛的應(yīng)用[5,9,12,13].但該方法由于受到線性響應(yīng)關(guān)系的制約,一般不能處理非線性過程顯著的二次污染問題.

針對大尺度區(qū)域,本研究提出了一種將線性規(guī)劃法與箱模型法、模擬法有機(jī)地結(jié)合起來的計算方法,使其相互取長補(bǔ)短,以解決上述單一方法應(yīng)用于大尺度區(qū)域時所遇到的困難.

1 研究方法

1.1 基本思路與技術(shù)路線

某種污染物的大氣環(huán)境容量是指所研究的區(qū)域內(nèi),在保證生態(tài)環(huán)境質(zhì)量滿足一定目標(biāo)的前提下,所能容納的污染物最大排放量.本研究的基本思路是將大尺度區(qū)域劃分為一系列小區(qū)域,將真實污染源簡化為源單元;基于線性規(guī)劃法尋求區(qū)域最大容量及優(yōu)化配置方案;采用箱模型以保證源單元尺度的環(huán)境質(zhì)量要求,采用模擬法以考慮二次污染對前體物環(huán)境容量的制約;將箱模型容量與模擬法容量作為規(guī)劃法的約束條件.上述思路的總體技術(shù)路線如圖1所示.

圖1 針對大尺度區(qū)域的大氣環(huán)境容量估算技術(shù)路線Fig.1 Technical route of the AEC estimation method for large-scale region

1.1.1 源單元的劃分 線性規(guī)劃法的一個重要基礎(chǔ)就是建立“源-受體”之間排污量與環(huán)境質(zhì)量的響應(yīng)關(guān)系.對于大尺度區(qū)域來說,由于污染源眾多,不確定性大,源-受體的響應(yīng)關(guān)系難以建立在真實的排放源基礎(chǔ)上.為此,本研究將大尺度區(qū)域劃分為一系列小區(qū)域.為便于資料收集與分析,上述小區(qū)域通常為某級別的行政區(qū)或工業(yè)區(qū).對于任一小區(qū)域,將類型相同、排放高度相近的排放源合并為一個源單元.此外,將區(qū)域內(nèi)一些排放高度高、排放量大、影響范圍廣的工業(yè)點源,如大型火電廠等,作為獨立的源單元.

1.1.2 與箱模型法的結(jié)合 在大尺度區(qū)域下,將一個源單元近似等同于一個點源是合理可行的,但不能保證源單元內(nèi)部的環(huán)境質(zhì)量.為保證局地空氣質(zhì)量滿足一定的要求,本研究在源單元尺度上采用空氣質(zhì)量箱模型,估算該源單元在滿足箱內(nèi)污染物濃度限值時的允許排放量,以下稱之為箱模型容量.將每個源單元的箱模型容量作為規(guī)劃模型的一種約束條件.

1.1.3 與模擬法的結(jié)合 對于大尺度區(qū)域來說,二次污染不可忽略,此時,“源-受體”響應(yīng)關(guān)系往往是非線性的.為了在線性規(guī)劃模型中考慮二次污染對區(qū)域環(huán)境容量的制約,本研究將常規(guī)的模擬法與線性規(guī)劃法相結(jié)合.模擬法的實質(zhì)是基于空氣質(zhì)量模式的多方案模擬,在二次污染因子滿足一定環(huán)境目標(biāo)的前提下,得到各排放源相應(yīng)前體污染物的允許排放量,在此基礎(chǔ)上得到全區(qū)域的某前體物的允許排放量,以下稱之為模擬法容量.將模擬法得到的各源的某(前體)污染物的允許排放量作為線性規(guī)劃法的一種約束條件.

1.2 模型與計算方法

依據(jù) 2.1所述思路與技術(shù)路線,本研究的容量估算模型總體上是基于線性規(guī)劃模型,其目標(biāo)函數(shù)為區(qū)域所能容納的污染物的最大排放量,約束條件包括環(huán)境質(zhì)量約束、基于經(jīng)濟(jì)技術(shù)等條件的排污量約束、基于箱模型容量的排污量約束、基于模擬法容量的排污量約束等.該模型的基本形式如下:

式(1)為目標(biāo)函數(shù).Qi表示源單元的污染物排放量,N為源單元總數(shù).PENALTY表示罰函數(shù),其定義見式(2).x表示罰函數(shù)的系數(shù),其值反映罰函數(shù)在目標(biāo)函數(shù)中的比重.引入罰函數(shù)的意義將在下文中進(jìn)一步介紹.

式(2)為罰函數(shù).Ckj為控制點j的濃度松弛變量(即控制點j的濃度允許超標(biāo)量).之所以要引入松弛變量,主要是因為對環(huán)境質(zhì)量約束來說,常規(guī)的“硬約束”具有明顯的不經(jīng)濟(jì)性,即可能由于局部控制點的制約導(dǎo)致整個地區(qū)環(huán)境容量顯著下降,同時也使削減費(fèi)用大幅上升.王勤耕等[14]提出了基于約束條件影響價格調(diào)整目標(biāo)值的方法,但該方法的計算結(jié)果具有一定的不確定性.本研究通過引入松弛變量,定義了一個反映區(qū)域少數(shù)點超標(biāo)水平的罰函數(shù),并將其加入目標(biāo)函數(shù)中進(jìn)行最小化.wj表示控制點 j的權(quán)重,在不區(qū)分控制點重要性的情形下,一般統(tǒng)一取值為1.

式(3)為區(qū)域一次污染的環(huán)境質(zhì)量約束,其中fij為源i對環(huán)境控制點j的響應(yīng)關(guān)系(即濃度傳輸系數(shù)),Caj為控制點 j的背景濃度,Csj為控制點 j的環(huán)境質(zhì)量目標(biāo).K為控制點的個數(shù).環(huán)境控制點包括人口集中區(qū)、生態(tài)敏感區(qū)和可能的地面濃度高值區(qū).本研究中采用法規(guī)級空氣質(zhì)量模式CALPUFF[15]計算源-受體的響應(yīng)關(guān)系.環(huán)境質(zhì)量控制目標(biāo)主要基于污染物的年均濃度標(biāo)準(zhǔn),對于缺乏法規(guī)標(biāo)準(zhǔn)的污染因子,可參照國內(nèi)外相關(guān)研究成果.

式(4)為依據(jù)源單元i具體的經(jīng)濟(jì)、技術(shù)和管理等條件提出的排放量上下界約束,Qimin和Qimax分別表示源單元i在給定條件下可能的最小和最大排放量.該約束條件可以有效保證容量規(guī)劃方案的現(xiàn)實可行性.

式(5)為基于箱模型法的排污量約束,以此保證源單元尺度的環(huán)境質(zhì)量滿足一定的控制目標(biāo).考慮到重點源的地理位置可能與一般的源單元(區(qū)域)相重合,而箱模型容量的約束應(yīng)考慮箱體所在區(qū)域內(nèi)的所有污染源.為此,用h表示箱體編號(h=1,2,…,H; H≤N,H為箱體總數(shù)),用i∈h表示源單元i位于箱體h內(nèi).Ah表示箱體h的箱模型環(huán)境容量.在忽略干濕沉積和化學(xué)衰變的情形下,箱體h的某污染物的大氣環(huán)境容量為[16]:

式中:Ah的單位為 104t/a;Csh為環(huán)境質(zhì)量目標(biāo),mg/m3;Cah為背景濃度,mg/m3;VE為混合層通風(fēng)系數(shù),m2/s;S為區(qū)域面積,km2.

式(7)為基于模擬法容量的排污量約束,以此保證區(qū)域二次污染滿足一定的控制目標(biāo).Mi表示各源的模擬法容量.目前較嚴(yán)重的二次污染因子主要包括酸沉降、二次氣溶膠粒子、臭氧(O3)等,其前體物包括二氧化硫(SO2)、氮氧化物(NOx)、一氧化碳(CO)、氨(NH3)、揮發(fā)性有機(jī)化合物(VOC)等[17].本研究重點關(guān)注 O3污染對 NOx和VOC排放的制約、細(xì)顆粒物(PM2.5)和酸雨污染對 SO2排放的制約.為此,采用 US EPA 研制的Models-3/CMAQ[18,19]對區(qū)域二次污染物濃度分布進(jìn)行模擬.根據(jù)各源對環(huán)境控制點濃度的貢獻(xiàn)率,按等比例對各源單元的排放量進(jìn)行多方案調(diào)整,直到滿足環(huán)境目標(biāo).這時各源單元的允許排放量為模擬法容量Mi.

式(8)為各控制點的松弛變量的上下界約束,Ckjmin和Ckjmax分別表示控制點j的環(huán)境質(zhì)量允許的最小和最大超標(biāo)量,實用中最小值一般取0.

2 應(yīng)用案例

2.1 案例概況

為了檢驗新方法的合理性和可行性,本研究針對我國東南沿海某區(qū)域,進(jìn)行了大氣環(huán)境容量的評估試驗.該區(qū)域總面積約14.6萬km2,包括13個地級市.該區(qū)域空氣質(zhì)量總體優(yōu)良,首要污染物為可吸入顆粒物,SO2和NO2濃度相對較低.酸雨污染相對嚴(yán)重,個別地區(qū)出現(xiàn)細(xì)顆粒物及O3超標(biāo)現(xiàn)象.本研究重點針對SO2、NO2、PM10和VOC4種污染物進(jìn)行環(huán)境容量評估,除考慮一次污染物的影響外,還考慮酸沉降、臭氧與PM2.5等二次污染對于上述有關(guān)前體污染物排放的制約.

以2007年為基準(zhǔn)年調(diào)研區(qū)域污染源、環(huán)境空氣質(zhì)量以及氣象等資料.根據(jù) 2.1節(jié)的技術(shù)路線對源單元進(jìn)行劃分,區(qū)域共有重點源 20個,一般源單元24個.選取環(huán)境質(zhì)量控制點58個.為從嚴(yán)控制,取各污染因子相應(yīng)標(biāo)準(zhǔn)限值的 75%作為環(huán)境質(zhì)量控制目標(biāo).

2.2 結(jié)果與討論

根據(jù)技術(shù)路線,首先可得到各源單元的模擬法容量及其所屬區(qū)域的箱模型容量.這兩種容量作為線性規(guī)劃法的約束條件進(jìn)入規(guī)劃計算,便可以得到最終的區(qū)域環(huán)境容量估算結(jié)果,同時得到各源單元的允許排放量.為簡化起見,本文僅以SO2和 NO2兩個污染因子為例,對計算結(jié)果進(jìn)行討論,PM10、VOC有類似結(jié)果,不再贅述.

考慮到實際工作中主要基于行政區(qū)開展產(chǎn)業(yè)規(guī)劃和污染物總量控制,對各源單元的規(guī)劃容量按13個市級行政區(qū)進(jìn)行分類匯總,表1和表2分別為SO2和NO2匯總后的結(jié)果.

表1 SO2環(huán)境容量及環(huán)境承載率Table 1 Environmental capacities and environmental bearing- ratios of SO2

表2 NO2環(huán)境容量及環(huán)境承載率Table 2 Environmental capacities and environmental bearing-ratios of NO2

表中“規(guī)劃容量”即為最終的容量估算結(jié)果,“箱模型容量”和“模擬法容量”作為規(guī)劃容量計算的兩種約束條件,也就是說,規(guī)劃法容量將不會大于這兩種容量.“制約因素”反映了容量估算結(jié)果主要受到何種因素的制約:若規(guī)劃容量等于箱模型容量,則說明主要受局地污染的影響;若規(guī)劃容量等于模擬法容量,則主要受某種二次污染的影響(其中包括 O3對于 NOx排放的限制,酸沉降和細(xì)顆粒物污染對于SO2的限制);若規(guī)劃容量顯著小于上述兩種容量,則說明主要受一次污染跨區(qū)域輸送的影響.表1和表2還給出了基準(zhǔn)年(2007)各地區(qū)環(huán)境容量的利用程度(本研究稱之為環(huán)境承載率,定義為污染物的排放量與環(huán)境容量的比值).

從計算結(jié)果可以看出,由于大尺度區(qū)域中空氣污染具有明顯的區(qū)域性和復(fù)合性,大氣環(huán)境容量受到多種因素的制約,具體來說,既有源單元尺度的局地影響,也有大尺度的跨區(qū)域影響;既有一次污染的直接影響,也有二次污染的間接影響.比如,對SO2來說,地區(qū)2的SO2容量主要受到區(qū)域細(xì)顆粒物污染的制約,地區(qū) 3的容量主要受源單元尺度的SO2直接污染的制約,地區(qū)4的容量則主要受SO2跨區(qū)域輸送的制約,地區(qū)10的SO2容量則主要受到區(qū)域酸沉降的制約. 對于 NO2來說,地區(qū) 3的容量主要受源單元尺度直接污染的制約,地區(qū)7的容量主要受跨區(qū)域輸送的制約,其它多數(shù)地區(qū) NO2的排放則主要受區(qū)域性臭氧污染的制約.由此可見,本研究提出的方法可以較全面地考慮大尺度區(qū)域下多種污染因素對于環(huán)境容量的制約.另外,從環(huán)境承載率的評估結(jié)果可以看出,不同地區(qū)環(huán)境容量的利用程度有較大差別,有一些地區(qū)的排放量已經(jīng)接近甚至超過環(huán)境容量,應(yīng)引起有關(guān)方面的關(guān)注,結(jié)合制約環(huán)境容量的主要因素,加強(qiáng)環(huán)境污染的綜合整治.

本研究提出的方法在實際應(yīng)用中應(yīng)關(guān)注以下問題:首先,箱模型法環(huán)境容量的計算結(jié)果往往具有較大的不確定性,除了混合層高度及通風(fēng)系數(shù)的影響以外,污染源分布及排放高度也有很大影響,比如,位于區(qū)域邊緣的高架點源,一般僅對箱內(nèi)污染產(chǎn)生部分影響.對此,必要時可以采用王勤耕等提出的改進(jìn)的箱模型法[20].其次,由于受到計算工作量的限制,模擬法容量估算只能是有限方案的比選,這在一定程度上會增加結(jié)果的不確定性,為此,實際應(yīng)用中應(yīng)開展盡可能多的情景模擬.不過,該問題的根本解決可能需要依賴于非線性規(guī)劃或準(zhǔn)線性規(guī)劃[21]的應(yīng)用研究.另外,環(huán)境控制目標(biāo)是環(huán)境容量估算的重要基礎(chǔ),限于當(dāng)前的認(rèn)識水平和研究能力,缺乏對污染因子累積性影響的足夠關(guān)注,有待于今后進(jìn)一步的加強(qiáng).

3 結(jié)語

本文針對大尺度區(qū)域的大氣環(huán)境容量評估問題,提出了將傳統(tǒng)的線性規(guī)劃法、箱模型法和模擬法有機(jī)結(jié)合的評估方法.該方法充分考慮了大尺度區(qū)域下空氣污染的區(qū)域性和復(fù)合性特征,同時兼顧了多種因素對于大氣環(huán)境容量的制約,有效克服了傳統(tǒng)單一方法所遇到的困難.案例研究結(jié)果表明,該方法合理可行,具有較好的推廣應(yīng)用前景.

[1] 王淑蘭,張遠(yuǎn)航,鐘流舉,等.珠江三角洲城市間空氣污染的相互影響 [J]. 中國環(huán)境科學(xué), 2005,25(2):133-137.

[2] 張 艷,余 琦,伏晴艷,等.長江三角洲區(qū)域輸送對上海市空氣質(zhì)量影響的特征分析 [J]. 中國環(huán)境科學(xué), 2010,30(7):914-923.

[3] HJ/T 131-2003開發(fā)區(qū)區(qū)域環(huán)境影響評價技術(shù)導(dǎo)則 [S].

[4] 李海晶.大氣環(huán)境容量估算及總量控制方法的研究進(jìn)展 [J].四川環(huán)境, 2007,26(1):67-71.

[5] 盧聰景,石曉楓.線性規(guī)劃法計算工業(yè)區(qū)大氣承載力的探討 [J].環(huán)境科學(xué)與技術(shù), 2008,31(10):142-144.

[6] 范紹佳,黃志興,劉嘉玲.大氣污染物排放總量控制 A-P值法及其應(yīng)用 [J]. 中國環(huán)境科學(xué), 1994,14(6):407-410.

[7] 鮑仙華,沈 虹,肖 青,等.區(qū)域環(huán)評中大氣環(huán)境容量及總量控制探討 [J]. 能源研究與信息, 2004,20(2):63-67.

[8] 魯東霞.工業(yè)集聚區(qū)發(fā)展規(guī)劃的大氣環(huán)境影響評價方法研究[J]. 氣象與環(huán)境科學(xué), 2008,31(1):89-92.

[9] Cheng S Y, Li J B, Feng B, et al. A gaussian-box modeling approach for urban air quality management in a Northern Chinese City—II. Pollutant emission abatement [J]. Water Air and Soil Pollution, 2006,178:15-36.

[10] 劉振忠,董 芃,王 麗.城市大氣環(huán)境污染預(yù)測與容量控制方法研究 [J]. 電站系統(tǒng)工程, 2005,21(4):17-19.

[11] 連 軍.多源模式在確定區(qū)域SO2環(huán)境容量中的應(yīng)用 [J]. 環(huán)境科學(xué)與管理, 2007,32(8):26-28.

[12] 王金南,潘向忠.線性規(guī)劃方法在環(huán)境容量資源分配中的應(yīng)用[J]. 環(huán)境科學(xué), 2005,26(6):195-198.

[13] 肖 楊,毛顯強(qiáng),馬根慧,等.基于 ADMS和線性規(guī)劃的區(qū)域大氣環(huán)境容量測算 [J]. 環(huán)境科學(xué)研究, 2008,21(3):13-16.

[14] 王勤耕,李宗愷.城市大氣污染源排污負(fù)荷優(yōu)化分配模型 [J].環(huán)境科學(xué)學(xué)報, 1996,16(4):444-449

[15] Scire J S, Strimaitis D G, Yamatino R J. A user’s guide for the CALPUFF dispersion model [EB/OL]. http://www.src.com/calpuff/calpuff1.htm, 2001-11-03.

[16] 國家環(huán)境保護(hù)局,中國環(huán)境科學(xué)研究院.城市大氣污染總量控制方法手冊 [M]. 北京:中國環(huán)境科學(xué)出版社, 1991.

[17] 曹國良,安心琴,周春紅,等.中國區(qū)域反應(yīng)性氣體排放源清單[J]. 中國環(huán)境科學(xué), 2010,30(7):900-906.

[18] Dennis R L, Byun D W, Novak J H, et al. The next generation of integrated air quality modeling: EPA's Models-3 [J]. Atmospheric Environment, 1996,30(12):1925-1938.

[19] Byun D W, Ching J K S. Science algorithms of the EPA Models-3 community multiscale air quality (CMAQ) modeling system [R].Washington, DC: Office of Research and Development, 1999.

[20] 王勤耕,吳躍明,李宗愷.一種改進(jìn)的 A-P控制法 [J]. 環(huán)境科學(xué)學(xué)報, 1997,17(3):278-283.

[21] Wagner F, Sch?pp W, Heyes C. The RAINS optimization module for the clean air for Europe (CAFE) programme [R]. Laxenburg:International Institute for Applied Systems Analysis, 2006.

An integrated method of atmospheric environmental capacity estimation for large-scale region.

QIAN Yue-dong, WANG Qin-geng*(State Key Laboratory of Pollution Control and Resource Reuse, School of the Environment, Nanjing University, Nanjing 210046, China). China Environmental Science, 2011,31(3)：504～509

According to the features of large-scale regions, a new method was proposed to estimate the atmospheric environmental capacity (AEC), which combines the traditional box-model method, simulation method and linear programming method. The basic idea of the new method is that: the large-scale region is divided into a series of small sub-regions, and pollution sources are classified into different source units. The linear programming method searches the maximum capacity and optimal allocation of the AEC. The box-model method ensures the local air quality in sub-regions meet certain targets. The simulation method limits precursor emissions which cause secondary pollution in the regional scale. Results from the last two methods act as restrictions of the linear programming method. The integrated method makes the three traditional methods to be complementary, and therefore has three significant merits: First, it considers both local and regional influences of pollution sources. Second, it considers both direct and secondary effects of the pollutants. Finally, it can give optimal schemes for the AEC allocation over different sub-regions. A case study conducted for a region in southeast China indicates that the new method is reasonable and feasible.

large-scale region；atmospheric environmental capacity；linear programming method；box-model；simulation method；environmental planning

X823

A

1000-6923(2011)03-0504-06

2010-07-16

國家“863”項目(2006AA06A307);環(huán)境保護(hù)部課題“海西經(jīng)濟(jì)區(qū)重點產(chǎn)業(yè)發(fā)展戰(zhàn)略環(huán)境評價”

* 責(zé)任作者, 教授, Wangqg@nju.edu.cn

錢躍東(1985-),男,浙江上虞人,南京大學(xué)環(huán)境學(xué)院碩士研究生,主要從事大氣環(huán)境容量與承載力的研究.