寧波-舟山港船舶排放清單及時空分布特征

尹佩玲,黃爭超,鄭丹楠,王雪松*,田旭東,鄭君瑜,張遠(yuǎn)航(.北京大學(xué)環(huán)境科學(xué)與工程學(xué)院,環(huán)境模擬與污染控制國家重點(diǎn)聯(lián)合實(shí)驗(yàn)室,北京 0087；.浙江省環(huán)境監(jiān)測中心,浙江 杭州 00；.華南理工大學(xué)環(huán)境與能源學(xué)院,廣東 廣州 50006)

寧波-舟山港船舶排放清單及時空分布特征

尹佩玲1,黃爭超1,鄭丹楠1,王雪松1*,田旭東2,鄭君瑜3,張遠(yuǎn)航1(1.北京大學(xué)環(huán)境科學(xué)與工程學(xué)院,環(huán)境模擬與污染控制國家重點(diǎn)聯(lián)合實(shí)驗(yàn)室,北京 100871；2.浙江省環(huán)境監(jiān)測中心,浙江 杭州 310012；3.華南理工大學(xué)環(huán)境與能源學(xué)院,廣東 廣州 510006)

為了掌握寧波-舟山港船舶的污染排放特征,支撐長江三角洲地區(qū)大氣污染控制工作,本研究采用基于船舶活動的排放因子法,收集船舶自動識別系統(tǒng)(AIS)提供的當(dāng)?shù)卮盎顒榆壽E數(shù)據(jù),結(jié)合勞氏船級社(LR)等機(jī)構(gòu)的船舶特征信息,建立了2010年寧波-舟山港船舶排放清單,獲得了船舶污染源排放的空間分布和小時變化特征譜.結(jié)果顯示,2010年寧波-舟山港船舶排放的SO2、NOx、PM10、PM2.5、VOC和CO總量分別為2.16×104、3.50×104、2.59×103、2.35×103、1.50×103和 3.01×103t. 集裝箱船、散干貨船、普通貨船、油輪和客船的排放貢獻(xiàn)占比分別為45.1%～57.7%、8.9%～12.9%、25.3%～40.7%、4.2%～8.1%和0.2%～0.7%.污染物的排放高峰期出現(xiàn)在9:00～14:00,不同類型船舶的逐時排放呈現(xiàn)不同的規(guī)律.排放高值區(qū)為由金塘水道、冊子水道、螺頭水道和蝦峙門水道組成的大型航道以及臨近的核心港區(qū).

寧波-舟山港；船舶排放源；排放清單；時空分布

航海運(yùn)輸在當(dāng)今經(jīng)濟(jì)全球化的發(fā)展中具有不可替代的作用.然而,船舶運(yùn)行消耗柴油和重油等燃料[1],會排放二氧化硫(SO2)、氮氧化物(NOx)、一氧化碳(CO)和碳?xì)?HC)等氣態(tài)污染物以及硫酸鹽、硝酸鹽和重金屬等顆粒態(tài)污染物[2].船舶排放不僅是港口地區(qū)顆粒物的重要來源之一[3-8],還對全球硫酸鹽濃度上升[9]、臭氧濃度增加[10]以及全球輻射強(qiáng)迫均有影響[11].

建立排放清單是評估船舶排放對空氣質(zhì)量的影響程度主要手段之一.1997年 Corbett等[12]率先估算了全球船舶SO2和NOx的排放量,至今歐洲[13]、地中海[14]、長灘島[15]及洛杉磯港[16-17]等多地均有較深入的船舶排放研究報道.

相對而言,國內(nèi)船舶排放清單的研究處于發(fā)展階段.21世紀(jì)初,國內(nèi)采用基于燃料消耗量的排放因子法,估算了全國內(nèi)河[18]、天津港[19]、京津冀內(nèi)河[20]、珠江三角洲[21]及中國水路運(yùn)輸業(yè)

[22]的船舶排放量.后來的研究逐步轉(zhuǎn)變?yōu)椴捎没诖盎顒拥呐欧乓蜃臃?如青島港[23]、廣東省[24]及上海港[25-26]等船舶排放清單,考慮了船舶的基本運(yùn)行模式,以反映船舶的活動特征.國外近年的研究[27-29]顯示,船舶自動識別系統(tǒng)(AIS)能夠提供更精細(xì)和準(zhǔn)確的船舶活動水平數(shù)據(jù),我國香港地區(qū)[30-31]、珠三角[32]、深圳[33]、大連港[34]、上海洋山港[35]、環(huán)渤海[36-37]及長三角和東海[38]等地區(qū)相應(yīng)地在建立船舶排放清單過程中結(jié)合了AIS數(shù)據(jù).

寧波-舟山港是我國東南沿海的重要港口,其貨物吞吐量 2009～2015年一直位居世界首位[39].繁忙船運(yùn)排放的污染物,可能對杭州灣乃至長江三角洲地區(qū)的空氣質(zhì)量造成不利影響.然而目前寧波-舟山港的船舶排放量尚不清晰,對該港口的各類型船舶排放特征及時空分布規(guī)律認(rèn)識不夠充分,因此有必要針對其開展更系統(tǒng)深入的研究.

本文以寧波-舟山港主要港區(qū)為研究區(qū)域,采用基于船舶活動的排放因子法,收集AIS提供的船舶活動數(shù)據(jù),結(jié)合勞氏船級社(Lloyd’s Register of Shipping)等機(jī)構(gòu)的船舶特征信息,建立浙江省寧波-舟山港 2010年主要類型船舶的污染物排放清單,進(jìn)一步研究在船舶類型?行駛模式及引擎功率方面的排放特征及各類船舶排放的時空分布特征,以期為寧波乃至長江三角洲地區(qū)的污染排放管理和大氣污染控制提供科學(xué)支撐.

1 材料和方法

1.1 研究區(qū)域和對象

根據(jù)美國環(huán)保署對船舶排放清單海域邊界的劃分建議[40],本文以寧波-舟山港 6個引航員登船點(diǎn)[41]為圓心、25n mile為半徑構(gòu)成的海域范圍作為研究區(qū)域(圖1).調(diào)研2010年寧波-舟山港客貨運(yùn)輸情況發(fā)現(xiàn),貨物總吞吐量中以散干貨(無包裝貨物,如煤炭、水泥和礦石等)和油品為主,分別占比37%和20%,集裝箱吞吐量占貨物總吞吐量達(dá)15.4%,港口間渡輪的旅客吞吐量為674.3萬人次[42],因此以集裝箱船、散干貨船、普通貨船、油輪和客船5類船舶作為研究對象,估算各類船舶SO2、NOx、PM10、PM2.5、CO和VOC等6種污染物的排放量,數(shù)據(jù)基準(zhǔn)年選為2010年.

圖1 研究范圍示意Fig.1 Sketch map of the study area圓圈代表本研究的寧波舟山港區(qū)域

1.2 船舶污染物排放量估算

1.2.1 估算方法 船舶排放的估算采用基于船舶活動的排放因子法,公式如下:

式中:h,i,j,k分別表示船舶類型、引擎功率、行駛模式和污染物種類;E為船舶的污染物排放總量g/a;VAN為船舶在估算期間的總抵港次數(shù);T為船舶的活動時間,h;P為船舶引擎功率,kW;LF表示船舶引擎的負(fù)載因子;EF為污染物排放因子, g/(kW?h).

船舶在一次抵港過程中行駛速度會發(fā)生變化,其主機(jī)引擎和輔助引擎的利用效率會有所差異,導(dǎo)致污染物排放速率的變化.為了反映船舶活動的各階段排放特征,根據(jù)航行速度將船舶活動過程劃分為4種行駛模式[30]:普通速度行駛(航行速度＞12nmi/h)?慢速行駛(8～12nmi/h)?機(jī)動行駛(1～8nmi/h)和靠泊(＜1nmi/h).對于某種類型船舶,首先分別計算單艘船在4種行駛模式下各類引擎的污染物排放量(主機(jī)引擎在靠泊時關(guān)閉,鍋爐在特定的慢速行駛?機(jī)動行駛和靠泊模式下開啟),4種行駛模式的排放量加和后獲得單艘次該類型船舶的排放量,再乘以該類型船舶2010年的總抵港次數(shù),獲得該類船舶的年排放總量.

1.2.2 排放因子等參數(shù)確定 由于缺乏研究地區(qū)排放因子等相關(guān)參數(shù)的實(shí)測資料,本文主要參考國內(nèi)外相關(guān)研究的結(jié)果[30,40].排放因子的選取由船舶的引擎類型、轉(zhuǎn)速、燃料類型及含硫率共同決定.除客船外,船舶的引擎包括主引擎、輔助引擎和鍋爐3部分,客船不包括鍋爐.引擎的類型?轉(zhuǎn)速和燃料類型通過調(diào)研該港口活動船舶在勞氏船級社中的數(shù)據(jù)獲得,將主機(jī)引擎區(qū)分為低速柴油機(jī)和中速柴油機(jī),輔助引擎和鍋爐類型不細(xì)分.燃料類型區(qū)分為船用柴油?重質(zhì)燃料油和船用汽油(船用汽油僅客船使用),根據(jù)我國船舶燃料的使用現(xiàn)狀[43],將前兩種燃油含硫率分別規(guī)定為1%～2%和3.5%,船用汽油采用文獻(xiàn)默認(rèn)值,采用的排放因子如表1所示.

表1 引擎和鍋爐的排放因子[30,40][g/(kW?h)]Table 1 Emission factors of engines and boilers[30,40][g/(kW?h)]

表2 寧波-舟山港各類型船舶引擎的平均功率*(kW)Table 2 Average power of engines and boilers of marine vessels in Ningbo-Zhoushan Port(kW)

除排放因子外,抵港次數(shù)、引擎功率和負(fù)載系數(shù)等參數(shù)采用如下方法確定:○1船舶抵港次數(shù):通過調(diào)研寧波-舟山港相關(guān)的文獻(xiàn)資料整理獲得;○2引擎功率:主機(jī)引擎功率通過調(diào)研勞氏船級社數(shù)據(jù)獲取,輔助引擎功率利用主輔機(jī)功率的比例進(jìn)行換算,鍋爐缺省值參考國內(nèi)相關(guān)研究成果,各類船舶引擎的平均功率如表2所示;○3負(fù)載系數(shù):表征發(fā)動機(jī)的實(shí)際使用情況,由最大航行速度和實(shí)際航行速度決定,區(qū)分為主引擎負(fù)載系數(shù)和輔助引擎負(fù)載系數(shù).

1.2.3 活動水平獲取 船舶排放的活動水平主要指船舶的航行時間,可通過船舶搭載的AIS獲取其瞬時地理位置、航線及速度等信息后間接計算得到.根據(jù)調(diào)研獲得各港口運(yùn)輸貨物和船舶的構(gòu)成情況,針對不同的船舶類型和噸位,利用船訊網(wǎng)提供的AIS信息查詢平臺,進(jìn)行船舶活動的抽樣調(diào)查,搜索樣本船舶在研究范圍內(nèi)一次進(jìn)出港過程的航行軌跡,記錄間隔4～6min的軌跡數(shù)據(jù),包括速度?經(jīng)緯度及時刻等信息.根據(jù)速度將各軌跡點(diǎn)區(qū)分至 4種行駛模式,表3是各類型船舶4種模式下活動時間的統(tǒng)計結(jié)果.

表3 寧波-舟山港各類型船舶不同行駛模式下的平均活動時間(h)Table 3 Average activity time in different modes of marine vessels in Ningbo-Zhoushan Port(h)

1.3 時空分配處理

根據(jù)上述方法,可以估算出同類型船舶4種行駛模式下的污染物排放量,統(tǒng)計出相應(yīng)行駛模式下的軌跡點(diǎn)數(shù),建立單個軌跡點(diǎn)和排放量的對應(yīng)關(guān)系.利用軌跡點(diǎn)的排放量和經(jīng)緯度信息,進(jìn)行寧波-舟山港船舶排放量的空間分配,應(yīng)用 GIS技術(shù),獲得 2km×2km分辨率的船舶排放清單.結(jié)合軌跡點(diǎn)的排放量和時刻信息,可統(tǒng)計出船舶排放在一天中的小時變化曲線.本文對船舶排放時空分配的處理具體至船舶類型、行駛模式以及污染物類型.

2 結(jié)果與討論

2.1 寧波-舟山港船舶排放清單

2010年寧波-舟山港船舶排放清單如表4所示,船舶污染源 SO2、NOx、PM10、PM2.5、VOC和CO的排放量分別為21574,34594,2593,2349, 1499和3009t,以排放SO2和NOx為主.SO2的排放與燃料中的含硫量直接相關(guān),其排放量大源于船舶使用的燃料含硫率高,本文使用的含硫率為1%～3.5%,國內(nèi)部分船舶所用重油的含硫率甚至高達(dá)4.5%[44].NOx的排放則受到燃料類型、發(fā)動機(jī)類型和負(fù)載等多方面因素影響.

2.1.1 船舶排放特征分析 不同污染物中各類型船舶的排放比重如圖2所示,其中集裝箱船和普通貨船的排放量所占比例最大,兩種船型共占比例為79%～86%,散干貨船、油輪和客船的排放貢獻(xiàn)占比分別為 8.9%～12.9%、4.2%～8.1%和0.2%～0.7%.集裝箱船中遠(yuǎn)洋船舶所占比例高,而遠(yuǎn)洋船舶使用的是重質(zhì)燃料油,含硫率在 1.5%～ 3.5%,所以SO2的排放量高,其次重質(zhì)燃油所含雜質(zhì)較多,燃燒時容易產(chǎn)生較多的顆粒物,故集裝箱船的PM10和PM2.5排放量也很大.另外,集裝箱船主要以中低速柴油機(jī)作為動力,由于大缸徑的低速機(jī)內(nèi)溫度比中小缸徑中高速機(jī)高,NOx生成的濃度更大[45],所以集裝箱船 NOx的排放量很大.集裝箱船單船排放量十分可觀,而普通貨船則是由于抵港次數(shù)多,活動頻繁,產(chǎn)生的污染物排放量也占有較大的比重,為25.3%～40.7%.

表4 2010年寧波-舟山港船舶排放清單(t)Table 4 Marine emission inventory of Ningbo-Zhoushan Port in 2010 (t)

圖2 寧波-舟山港不同類型船舶污染物排放分擔(dān)率Fig.2 Emission contribution rates by marine vessel types in Ningbo-Zhoushan Port

不同動力設(shè)備對各類污染物的排放分擔(dān)率如圖 3所示,總體而言,主機(jī)引擎的排放量最大,其次是輔助引擎,鍋爐最小,貢獻(xiàn)占比分別為49.4%～74.8%,22.6%～33.3%和 2.2%～20.1%.主機(jī)引擎作為船舶活動的主動力,在普通速度行駛、慢速行駛和機(jī)動行駛階段都開啟,運(yùn)行時間長且功率大,所以排放量大.鍋爐排放的主要影響因素是燃料,與運(yùn)行工況關(guān)系不大,主要體現(xiàn)為重質(zhì)燃料油燃燒的排放特征,即對SO2和顆粒物的排放貢獻(xiàn)較大.

圖3 寧波-舟山港船舶不同引擎的污染物排放分擔(dān)率Fig.3 Emission contribution rates by engine types in Ningbo-Zhoushan Port

4種行駛模式對各類污染物的排放分擔(dān)率如圖4所示.普通速度行駛與慢速行駛的排放比重相當(dāng),均約占 30%,其次是靠泊,分擔(dān)率最小的是機(jī)動行駛.普通速度行駛、慢速行駛和機(jī)動行駛的主輔引擎負(fù)載依次降低,燃燒充分程度逐漸下降,因此 NOx和 PM的排放強(qiáng)度減弱,CO和VOC的排放強(qiáng)度增加,可以看出,機(jī)動行駛對VOC和CO的貢獻(xiàn)比重較其他污染物大.SO2的排放與燃料的消耗量直接相關(guān),靠泊時鍋爐開啟以維持船舶的電力供應(yīng),因此SO2在該模式下的排放比重較大.

圖4 寧波-舟山港船舶不同行駛模式的污染物排放分擔(dān)率Fig.4 Emission contribution rates by operation modes in Ningbo-Zhoushan Port

2.1.2 與其他港口船舶排放的比較 對國內(nèi)不同港口船舶排放清單的總量比較,會受到各港口船舶構(gòu)成的嚴(yán)重影響.針對各港口集裝箱船排放的對比情況,如表5所示.

表5 國內(nèi)部分沿海港口集裝箱船污染物排放情況Table 5 The throughputs and emissions of container vessels in domestic costal ports

圖5 2010年國內(nèi)3大集裝箱港口單位吞吐量的排放情況Fig.5 The per-throughput emissions of the domestic top three container vessels ports in 2010

8個港口都采用的是基于船舶活動水平的排放因子法,除上海港外,其余港口都通過AIS獲取活動水平數(shù)據(jù).在沿海和遠(yuǎn)洋的集裝箱船估算上,上海港、上海洋山港、天津港和大連港采用同樣來源的排放因子,僅天津港的顆粒物有所差異,其余4個港口采用同樣來源的排放因子,兩個來源的排放因子差異不大,故排放量的排序與吞吐量排序較對應(yīng).

但香港、珠三角、深圳及本研究通過調(diào)研將集裝箱船區(qū)分了重質(zhì)燃料油和船用柴油等類別,使得平均排放因子降低,而洋山港、天津港和大連港僅選用重質(zhì)燃料油對應(yīng)的排放因子,這是導(dǎo)致香港的排放量較低、天津和大連港的排放量較高的影響原因之一,另外有研究證明大連港的船舶排放可能存在高估情況[37].

2010年我國集裝箱吞吐量前3位港口分別是上海港、深圳港和寧波-舟山港,對應(yīng)的污染物排放量排序相同,說明寧波-舟山港的集裝箱船估算量處于較為合理的水平.進(jìn)一步比較3個港口單位吞吐量的污染物排放量(如圖5所示)可知,寧波-舟山港大于上海港和深圳港,與臨近的上海港更接近.上海港基于調(diào)研結(jié)果在估算沿海船舶時燃油含硫率定為 0.1%～1.5%[26],低于本研究采用的含硫率,使用了燃油校正因子降低了排放因子,而深圳港的集裝箱船未將鍋爐排放納入估算[33].

2.2 寧波-舟山港船舶排放時空分布特征

2.2.1 時間變化特征 各類型船舶排放的小時變化情況如圖6所示(以SO2為例).從所有船型的綜合排放日變化可以看出,污染物的排放集中在日間,高排放時段為 09:00～14:00,夜間及凌晨均有持續(xù)的小幅排放.該日變化特征是5類船舶共同排放的結(jié)果,不同類型船舶排放的小時變化特征則存在一定的差異.

圖6 寧波-舟山港船舶排放的小時變化Fig.6 Hourly profiles of marine emissions by vessel categories in Ningbo-Zhoushan Port

圖7 寧波-舟山港船舶SO2排放空間分布Fig.7 Spatial distribution of total marine SO2emission in Ningbo-Zhoushan Port

集裝箱船和散干貨船的噸位一般較大,并且有固定的船期和掛靠港口,而油輪運(yùn)載的是原油等危險性較高的貨物,因此基于安全的原因,這 3類船舶的活動集中在日間,所以日間排放量較大,而夜間的排放量主要來自于靠泊.普通貨船是寧波-舟山港航行最頻繁的船型,其噸位不大,大部分為中小型貨船,受晝夜的限制較小,因此排放量的小時變化幅度不明顯,維持在3.4%～5.5%.客船主要為來往于舟山-寧波之間的渡輪,從早晨始發(fā)到傍晚停航,且渡輪停航后不需要引擎或鍋爐維持運(yùn)轉(zhuǎn),因此客船只有日間航行期間的排放.

2.2.2 空間分布特征 寧波-舟山港船舶排放的空間分布如圖7所示(以SO2為例,由于SO2的排放量較大,且其他污染物的空間分布與之相似),大型港口如金塘?鎮(zhèn)海?北侖?大榭?穿山島?沈家門?定海和老塘山港口分別位于在寧波市穿山半島和舟山本島上,共用金塘水道?冊子水道以及螺頭水道,形成核心港區(qū),往來船舶頻繁,因此該區(qū)域的污染物排放量十分可觀.5種類型船舶進(jìn)入核心港區(qū)的通道主要有4條,分別是“蝦峙門航道-螺頭水道-港區(qū)”、“佛渡水道-螺頭水道-港區(qū)”、“杭州灣(嘉興方向)-金塘水道-港區(qū)”和“杭州灣(洋山方向)-冊子水道-港區(qū)”(依次如圖8中A～D所示),其中以“蝦峙門航道-螺頭水道-港區(qū)”排放強(qiáng)度最顯著,這是由于此通道通航能力達(dá)到30萬t級,為航行排放量較大的集裝箱船?大型貨船和油輪進(jìn)入港區(qū)的主要路徑.

圖8 進(jìn)入寧波-舟山港核心港區(qū)主要通道Fig.8 The main lanes to enter the main part of Ningbo-Zhoushan Port

各類型船舶排放的空間分布情況如圖9所示(以 SO2為例).集裝箱船和散干貨船的污染物排放集中在寧波市穿山半島側(cè)的港口及附近航道,與寧波-舟山港的總體規(guī)劃[49]相吻合,其規(guī)定在整合運(yùn)輸功能后將集裝箱船運(yùn)輸系統(tǒng)布局在鎮(zhèn)海、北侖、穿山、大榭、金塘、梅山和六橫港區(qū),鐵礦石、煤炭等散雜貨運(yùn)輸布局在北侖、鎮(zhèn)海、大榭島和舟山本島的老塘山等,這些噸位較大的遠(yuǎn)洋船舶途徑蝦峙門外航道、蝦峙門航道、螺頭水道、冊子水道和金塘水道進(jìn)出港區(qū).普通貨船主要航行于寧波-舟山港的港口之間,進(jìn)行沿海的貨物運(yùn)輸,相關(guān)港口和海域排放強(qiáng)度較大,部分排放分布在通往嘉興、洋山乃至上海方向的航線上.油輪的排放分布體現(xiàn)出遠(yuǎn)洋和沿海兼具的特點(diǎn),南面2個遠(yuǎn)洋航道和北向沿海航道均有排放,但排放較多集中在港口附近,這是由于油輪在港口停泊的時間較長,鍋爐和輔助引擎的排放較多,這些港區(qū)除了綜合性的寧波側(cè)港口外還包括以原油、成品油儲存及中轉(zhuǎn)運(yùn)輸為主的老塘山、定海、馬岙和六橫港區(qū).往來于舟山和寧波的渡輪2條主要航線為“鴨白線(舟山鴨蛋山-寧波白峰)”和“沙岙-郭巨線(舟山沙岙-寧波郭巨)”[42]從圖9(d)中能清楚地看出客船在航線上的排放,由于渡輪每日多次往返且停泊的時間較短,所以排放強(qiáng)度主要體現(xiàn)在航行過程中.

2.3 不確定性分析

本文采用評級的方法對各類型船舶的排放因子、活動水平以及排放量進(jìn)行半定量的不確定性評估[21,26,32,34,50-52],為未來進(jìn)一步完善清單提供方向.盡管各類型船舶估算參數(shù)的獲取手段一致,但是受到資料完善程度以及調(diào)研樣本的數(shù)量和質(zhì)量等因素影響,不同類型船舶的排放因子和活動水平的可靠性和本地化程度不同,因此其不確定性也有所差異,具體表現(xiàn)如表6所示.

圖9 寧波-舟山港不同類型船舶SO2排放空間分布Fig.9 Spatial distribution of marine SO2emission by vessel categories in Ningbo-Zhoushan Port

評級從A～E表征數(shù)據(jù)的可靠性或本地化程度逐漸降低,不確定性增加.本研究清單的不確定性來源主要來自抵港次數(shù)、排放因子、活動時間及其他估算參數(shù)4個方面.1)抵港次數(shù):對于集裝箱船、油輪和客船運(yùn)載對象較明確的類型,港口、海事部門及文獻(xiàn)統(tǒng)計資料較為齊全,具有一定的可靠性和本地化程度,相較而言,貨船的數(shù)量繁多且以中小型為主,統(tǒng)計的難度加大,不確定性增加;2)排放因子和其他估算參數(shù):主要參考國內(nèi)外的相關(guān)研究成果,集裝箱船、散干貨船和油輪 3種大型的船舶,與國外同類型船舶的差異較普通貨船和客船(渡輪)小,且遠(yuǎn)洋船舶占有一定的比例.雖然結(jié)合了本地化的引擎和燃料信息調(diào)研,但難以涵蓋船舶的具體特征信息,會帶來較大的不確定性;3)活動時間:針對航線特征較為明顯的集裝箱船、油輪和客船,樣本能夠較大程度地反映船舶實(shí)際的航行狀態(tài),而對于機(jī)動靈活的普通貨船,調(diào)研樣本具有較大的不確定性.基于上述對排放因子和活動水平的不確定性分析,獲得對排放量的不確定性評級(如表 6所示),其中普通貨船排放量的不確定性相對較高.

本研究建立的船舶排放清單具有一定的不確定性,但能夠反映寧波-舟山港船舶污染源排放的基本情況,今后仍需提高活動水平的獲取性和排放因子的本地化,有針對性地在上述方面進(jìn)一步降低清單的不確定性.

表6 2010年寧波-舟山港船舶污染源排放清單的不確定性評級Table 6 Uncertainty rating of the 2010 marine emission inventory in Ningbo-Zhoushan Port

3 結(jié)論

3.1 寧波-舟山港 2010年船舶排放的 SO2、NOx、PM10、PM2.5、VOC和 CO 分別為21574,34594,2593,2349,1499和3009t.

3.2 在各類污染物的排放量中,集裝箱船?普通貨船?散干貨船?油輪和客船的排放貢獻(xiàn)占比分別為 45.1%～57.7%、25.3%～40.7%、8.9%～12.9%、4.2%～8.1%和0.2%～0.7%.不同船舶引擎類型的排放貢獻(xiàn)順序:主機(jī)引擎＞輔助引擎＞鍋爐.不同運(yùn)行模式的排放大小依次為:普通速度行駛＞慢速行駛＞靠泊＞機(jī)動行駛.

3.3 船舶整體排放的高峰期出現(xiàn)在9:00～14:00,這與集裝箱船、散干貨船和油輪的排放規(guī)律較接近,而貨船晝夜變化不大,逐時排放比例維持在3.4%～5.5%,客船則體現(xiàn)出僅日間航行排放的特點(diǎn).

3.4 船舶排放的高值區(qū)出現(xiàn)在由金塘水道、冊子水道、螺頭水道和蝦峙門水道組成的大型航道以及臨近的核心港區(qū),不同類型船舶基于各自的運(yùn)輸功能以及航線體現(xiàn)出不同的空間排放特征.

[1] 翟 君,馮立巖,王 猛,等.氣體燃料發(fā)動機(jī)發(fā)展對中國溫室氣體減排貢獻(xiàn)的生命周期分析 [J]. 中國環(huán)境科學(xué), 2015,35(1): 62-71.

[2] 萬 霖,何凌燕,黃曉鋒.船舶大氣污染排放的研究進(jìn)展 [J]. 環(huán)境科學(xué)與技術(shù), 2013,36(5):57-62.

[3] Pandolfi M, Gonzalez-Castanedo Y, Alastuey A, et al. Source apportionment of PM10and PM2.5at multiple sites in the strait of Gibraltar by PMF: impact of shipping emissions [J]. Environmental Science and Pollution Research, 2011,18(2):260-269.

[4] Minguillón M C, Arhami M, Schauer J J, et al. Seasonal and spatial variations of sources of fine and quasi-ultrafine particulate matter in neighborhoods near the Los Angeles-Long Beach harbor [J]. Atmospheric Environment, 2008,42(32):7317-7328.

[5] Kim E, Hopke P K. Source characterization of ambient fine particles at multiple sites in the Seattle area [J]. Atmospheric Environment, 2008,42(24):6047-6056.

[6] Yuan Z, Lau A K H, Zhang H, et al. Identification and spatiotemporal variations of dominant PM10sources over Hong Kong [J]. Atmospheric Environment, 2006,40(10):1803-1815.

[7] Mazzei F, D'Alessandro A, Lucarelli F, et al. Characterization of particulate matter sources in an urban environment [J]. Science of the Total Environment, 2008,401(1):81-89.

[8] Viana M, Amato F, Alastuey A, et al. Chemical tracers of particulate emissions from commercial shipping [J]. Environmental science & technology, 2009,43(19):7472-7477.

[9] Lauer A, Eyring V, Hendricks J, et al. Global model simulations of the impact of ocean-going ships on aerosols, clouds, and the radiation budget [J]. Atmospheric Chemistry and Physics, 2007, 7(19):5061-5079.

[10] Dals?ren S B, Endresen ?, Isaksen I S, et al. Environmental impacts of the expected increase in sea transportation, with a particular focus on oil and gas scenarios for Norway and northwest Russia [J]. Journal of Geophysical Research: Atmospheres (1984-2012), 2007,112(D2).

[11] Eyring V, Isaksen I S, Berntsen T, et al. Transport impacts on atmosphere and climate: Shipping [J]. Atmospheric Environment, 2010,44(37):4735-4771.

[12] Corbett J J, Fischbeck P. Emissions from Ships [J]. Science, 1997,278(5339):823-824.

[13] Entec. Quantification of Emissions from Ships Associated with Ship Movements between Ports in the European Community.Final Report [R]. United Kingdom: Entec, 2002.

[14] Entec. Ship Emissions Inventory - Mediterranean Sea [R]. United Kingdom: Concawe, 2007.

[15] Starcrest Consulting Group L. Port of long beach air emissions inventory—2009 [R]. The Philippines: The Port of Long Beach, 2010.

[16] Starcrest Consulting Group L. Port of Los Angeles inventory of air emissions-2010 [R]. Los Angeles: The port of Los Angeles, 2011.

[17] Starcrest Consulting Group L. Port of Los Angeles inventory of air emissions-2012 [R]. Los Angeles: The port of Los Angeles, 2013.

[18] 傅立新,程玲琳,粘桂蓮.移動污染源大氣環(huán)境影響研究報告[R]. 北京:清華大學(xué), 2005.

[19] 金陶勝,殷小鴿,許 嘉,等.天津港運(yùn)輸船舶大氣污染物排放清單 [J]. 海洋環(huán)境科學(xué), 2009,28(6):623-625.

[20] 隗 瀟.京津冀非道路移動源排放清單的建立 [C]//2013中國環(huán)境科學(xué)學(xué)會學(xué)術(shù)年會論文集.北京:中國環(huán)境科學(xué)學(xué)會, 2013.

[21] 張禮俊,鄭君瑜,尹沙沙,等.珠江三角洲非道路移動源排放清單開發(fā) [J]. 環(huán)境科學(xué), 2010,31(4):886-891.

[22] 邢 輝,段樹林,黃連忠,等.中國水路運(yùn)輸業(yè)能源消耗與廢氣排放測算 [J]. 中國環(huán)境科學(xué), 2016,36(6):1913-1920.

[23] 劉 靜,王 靜,宋傳真,等.青島市港口船舶大氣污染排放清單的建立及應(yīng)用 [J]. 中國環(huán)境監(jiān)測, 2011,27(3):50-53.

[24] 葉斯琪,鄭君瑜,潘月云,等.廣東省船舶排放源清單及時空分布特征研究 [J]. 環(huán)境科學(xué)學(xué)報, 2014,34(3):537-547.

[25] Yang D, Kwan S H, Lu T, et al. An emission inventory of marine vessels in Shanghai in 2003 [J]. Environmental Science & Technology, 2007,41(15):5183-5190.

[26] 伏晴艷,沈 寅,張 健.上海港船舶大氣污染物排放清單研究[J]. 安全與環(huán)境學(xué)報, 2012,12(5):57-63.

[27] Heather M. Perez R C A R. Automatic identification systems (AIS) data use in marine vessel emission estimation [R]. Baltimore, Maryland: 2009.

[28] Jalkanen J P, Brink A, Kalli J, et al. A modelling system for the exhaust emissions of marine traffic and its application in the Baltic Sea area [J]. Atmospheric Chemistry and Physics, 2009, 9(23):9209-9223.

[29] Johansson L, Jalkanen J, Kalli J, et al. The evolution of shipping emissions and the costs of regulation changes in the northern EU area [J]. Atmospheric Chemistry and Physics, 2013,13(22): 11375-11389.

[30] Ng S K W, Lin C, Chan J W M, et al. Study on marine vessels emission inventory, Final Report [R]. Hong Kong: The Hong Kong Environmental Protection Department, 2012.

[31] Yau P S, Lee S C, Corbett J J, et al. Estimation of exhaust emission from ocean-going vessels in Hong Kong [J]. Science of the total environment, 2012,431:299-306.

[32] 葉斯琪.珠江三角洲地區(qū)船舶排放特征及對區(qū)域空氣質(zhì)量影響的研究 [D]. 廣州:華南理工大學(xué), 2014.

[33] 楊 靜,尹佩玲,葉斯琪,等.深圳市船舶排放清單與時空特征研究 [J]. 環(huán)境科學(xué), 2015,36(4):1217-1226.

[34] 譚建偉,宋亞楠,葛蘊(yùn)珊,等.大連海域遠(yuǎn)洋船舶排放清單 [J]. 環(huán)境科學(xué)研究, 2014,27(12):1426-1431.

[35] Song S. Ship emissions inventory, social cost and eco-efficiency in Shanghai Yangshan port [J]. Atmospheric Environment, 2014, 82:288-297.

[36] 宋亞楠.內(nèi)河和近海船舶排放特性及排放清單研究 [D]. 北京:北京理工大學(xué), 2015.

[37] 邢 輝,段樹林,黃連忠,等.基于 AIS數(shù)據(jù)的渤海灣地區(qū)海運(yùn)排放測算 [J]. 中國環(huán)境科學(xué), 2016,36(3):953-960.

[38] Fan Q, Zhang Y, Ma W, et al. Spatial and seasonal dynamics of ship emissions over the Yangtze River Delta and East China Sea and their potential environmental influence [J]. Environmental Science & Technology, 2016,50(3):1322-1329.

[39] 浙江日報.寧波舟山港實(shí)現(xiàn)新突破-貨物吞吐量連續(xù)七年居世界 首 位 [EB/OL]. http://zjrb.zjol.com.cn/html/2016-01/17/ content_2942826.htm?div=-1.

[40] ICF International. Current methodologies in preparing mobile source port-related emission inventories final report [R]. United States: EPA, 2009.

[41] 寧波海事局. 寧波－舟山港引航員登輪點(diǎn)[EB/OL]. http://www. nbmsa.gov.cn/portal/section/17/2011/11/18/356faa3a17c54c0fa86 174d724a0ca46/.

[42] 中國港口年鑒編輯部.中國港口年鑒2011 [M]. 北京:中國港口雜志社, 2011.

[43] 交通運(yùn)輸部水運(yùn)科學(xué)研究院.中國船舶和港口空氣污染防治概況 [R]. 北京:交通運(yùn)輸部水運(yùn)科學(xué)研究院, 2013.

[44] 航運(yùn)交易公報.關(guān)于船舶污染的困惑 [EB/OL]. http://info. jctrans.com/news/hyxw/20154152129547.shtml.

[45] 孫洲陽,陳景鋒,黃加亮,等.船舶柴油機(jī)運(yùn)行中的排放控制技術(shù)研討 [J]. 航海技術(shù), 2001(4):64-66.

[46] Todd Sax, Andrew A. A critical review of ocean-going vessel particulate matter emission factors [R]. California: California Air Resources Board, 2007.

[47] IVL. Methodology for calculating emissions from ships [R]. Sweden: Swedish Environmental Protection Agency, 2004.

[48] Memo from Chris Lindhjem of ENVIRON, PM Emission Factors [Z]. 2005.

[49] 舟山港航網(wǎng).寧波-舟山港總體規(guī)劃 [EB/OL]. http://port. zhoushan.gov.cn/news/20110726103045697.html.

[50] 張禮俊,鄭君瑜,尹沙沙,等.珠江三角洲非道路移動源排放清單開發(fā) [J]. 環(huán)境科學(xué), 2010,31(4):886-891.

[51] Agency E E. EMEP/EEA air pollutant emission inventory guidebook-2013 [EB/OL]. http://www.eea.europa.eu//publications/ emep-eea-guidebook-2013.

[52] 鄭君瑜,張禮俊,鐘流舉,等.珠江三角洲大氣面源排放清單及空間分布特征 [J]. 中國環(huán)境科學(xué), 2009,29(5):455-460.

Marine vessel emission and its temporal and spatial distribution characteristics in Ningbo-Zhoushan Port.

YIN Pei-ling1, HUANG Zheng-chao1, ZHENG Dan-nan1, WANG Xue-song1*, TIAN Xu-dong2, ZHENG Jun-yu3, ZHANG Yuan-hang1
(1.State Key Joint Laboratory of Environmental Simulation and Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing 100871, China；2.Zhejiang Environmental Monitoring Center, Hangzhou 310012, China；3.College of Environment and Energy, South China University of Technology, Guangzhou 510006, China). China Environmental Science, 2017,37(1)：27～37

A marine vessel emission inventory of 2010 in Ningbo-Zhoushan Port was developed based on an activity-based approach with tracks of local vessels from AIS (automatic identification system) and vessel feature information from LR (Lloyd’s Register of Shipping). The temporal and spatial distribution characteristics of marine vessel emission were acquired. The results showed that the emissions of SO2, NOx, PM10, PM2.5, VOC and CO were 2.16×104, 3.50×104, 2.59×103, 2.35×103, 1.50×103and 3.01×103t, respectively. The contribution percentages of exhausted emissions from container vessel, dry bulk carrier, conventional cargo vessel, oil tanker and cruise/ferry were 45.1%～57.7%, 8.9%～12.9%, 25.3%～40.7%, 4.2%～8.1% and 0.2%～0.7%, respectively. The emissions of pollutants usually reached peaks during 9:00～14:00, and the hourly emission profiles varied with vessel types. The high emissions were found in the large channel composed of Jintang, Cezi, Luotou and Xiazhimen waterways, as well as the surrounding main ports.

Ningbo-Zhoushan Port；marine vessel source；emission inventory；temporal and spatial distribution

X51

A

1000-6923(2017)01-0027-11

尹佩玲(1991-),女,廣東東莞人,北京大學(xué)碩士研究生,主要從事區(qū)域大氣污染控制工作研究.

2016-03-28

國家科技支撐計劃(2014BAC06B02);國家自然科學(xué)基金項(xiàng)目(41175102)

* 責(zé)任作者, 副教授, xswang@pku.edu.cn