• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    關中城市群道路移動源氣態(tài)污染物排放特征

    2019-02-27 09:12:32陶雙成鄧順熙郝艷召高碩晗熊新竹孔亞平
    中國環(huán)境科學 2019年2期
    關鍵詞:小客車氣態(tài)貨車

    陶雙成,鄧順熙,郝艷召,高碩晗,熊新竹,孔亞平

    關中城市群道路移動源氣態(tài)污染物排放特征

    陶雙成1*,鄧順熙2,郝艷召3,高碩晗1,熊新竹1,孔亞平1

    (1.交通運輸部科學研究院,北京 100029;2.長安大學環(huán)境科學與工程學院,陜西 西安 710064;3.長安大學汽車學院,陜西 西安 710064)

    基于MOVES模型對參數(shù)進行本地化修正,計算機動車排放的氣態(tài)污染物排放因子,以2012年為基準年建立關中城市群的道路移動源常規(guī)和非常規(guī)氣態(tài)污染物的排放總量清單,并得到不同車型、不同城市區(qū)域的機動車污染物排放分擔率.結(jié)果顯示,關中城市群的道路移動源常規(guī)氣態(tài)污染物中一氧化碳(CO)的排放量為45.40萬t,氮氧化物(NO)為8.190萬t、二氧化硫(SO2)為0.420萬t、氨(NH3)為0.10萬t;非常規(guī)氣態(tài)污染物中非甲烷碳氫化合物(NMHC)的排放量為4.168萬t,甲醛(HCHO)為0.057萬t、乙醛(CH3CHO)為0.027萬t、丙烯醛(C3H4O)為0.004萬t、1,3-丁二烯(C4H6)為0.012萬t、苯為0.090萬t、甲烷(CH4)為0.123萬t、氧化亞氮(N2O)為0.004萬t.此外,各城市按照排放分擔率從高至低依次為西安(50%)、渭南(23%)、咸陽(含楊凌)(12%)、寶雞(10%)和銅川(5%).本研究還發(fā)現(xiàn)污染物排放分擔率在不同車型中差異顯著,其中NO排放以重型貨車(33.85%)和中型貨車(21.21%)為主;SO2、醛類物質(zhì)在重型貨車中排放分擔率分別為31.31%和30%;而CO、NMHC、C4H6、苯和CH4的排放主要來自小客車(分別為32.86%、17.55%、26.64%、26.45%和38.85%)和摩托車(分別為32.64%、55.21%、43.29%、49.04%和30.97%);NH3的小客車和重型貨車排放分擔率分別為49.5%和31.31%.

    道路移動源排放;氣態(tài)污染物;排放特征;MOVES模型;關中城市群

    隨著我國城鎮(zhèn)化進程的推進,城市連片區(qū)域環(huán)境空氣污染問題日益嚴重,主要表現(xiàn)為區(qū)域重污染天氣事件頻發(fā),顆粒物總體濃度高、化學組成復雜,臭氧(O3)季節(jié)性、區(qū)域性污染問題日趨嚴重,嚴重影響人民群眾的身心健康和日常生產(chǎn)生活[1-2].由于道路移動源是城市群內(nèi)環(huán)境空氣污染時空變化的主要影響因素之一[3],摸清區(qū)域道路移動源排放特征是識別和解析污染來源、支撐空氣污染預警預報[4]、分析解釋觀測結(jié)果和制定減排控制方案的重要基礎.國內(nèi)已經(jīng)在珠江三角洲地區(qū)[5]、長株潭地區(qū)[6]、京津冀地區(qū)[7-8]、海峽西岸經(jīng)濟區(qū)[9]等大氣污染主要區(qū)域,以及廊坊[10]、杭州[11-12]、成都[13]、南昌[14]、青島[15]等典型城市開展了機動車污染排放特征研究,初步摸清了機動車排放常規(guī)氣態(tài)污染物一氧化碳(CO)、氮氧化物(NO)與二氧化硫(SO2)以及顆粒物(PM2.5和PM10)排放特征,部分研究還計算了氧化亞氮(N2O)、二氧化碳(CO2)等溫室氣體排放清單[16].這些研究為厘清典型區(qū)域和城市的機動車污染物排放狀況,合理制定機動車污染防治策略提供了理論支持和決策依據(jù).已有研究多集中在機動車常規(guī)氣態(tài)污染物和顆粒物排放特征方面,對非甲烷碳氫化合物(NMHC)、醛類物質(zhì)等非常規(guī)氣態(tài)污染物研究較少,研究表明機動車排放的非甲烷碳氫化合物(NMHC)、甲醛(HCHO)、乙醛(CH3CHO)、丙烯醛(C3H4O)、1,3-丁二烯(C4H6)、苯(C6H6)、氧化亞氮(N2O)和氨(NH3)等氣態(tài)污染物具有較強的光反應活性[17],在光照或與羥基(OH·)、O3結(jié)合條件下會迅速發(fā)生反應,是城市光化學煙霧污染O3和過氧乙酰硝酸酯(PAN)以及二次氣溶膠形成的重要前體物,在城市光化學煙霧污染[18-20]和二次氣溶膠[21]形成中起著關鍵作用.另據(jù)報道,醛類物質(zhì)中的甲醛、乙醛、丙烯醛和苯甲醛等都是重要的致癌物質(zhì)[22-23],隨著區(qū)域性機動車保有量的爆發(fā)式增長,城市大氣中醛類物質(zhì)等氣態(tài)污染物對環(huán)境和人體健康的危害會越來越嚴重,必須引起足夠重視.

    關中城市群是陜西省的工業(yè)、經(jīng)濟、文化和教育中心,但天然的“U”字型河谷地形極不利于本地大氣污染物的擴散,易形成污染物累積并發(fā)生二次污染[24].近年來,關中地區(qū)光化學煙霧污染頻發(fā)、夏季O3濃度超標明顯,O3已經(jīng)超過PM2.5和PM10成為該地區(qū)當前占比最高的污染物[25],隨著機動車保有量的迅速增加,機動車排放對以西安為中心的城市群內(nèi)大氣光化學煙霧污染的貢獻愈發(fā)顯著并將嚴重影響該區(qū)域的經(jīng)濟可持續(xù)發(fā)展和居民健康.建立關中城市群道路移動源氣態(tài)污染物排放清單,分析機動車排放的常規(guī)和非常規(guī)氣態(tài)污染物的基本特征和分擔率,是做好區(qū)域大氣污染防治特別是光化學煙霧等二次污染防控工作的基礎,目前針對關中城市群機動車排放氣態(tài)污染物的研究比較缺乏.本研究在對MOVES模型輸入?yún)?shù)本土化修正的基礎上,以2012年為基準年建立了關中城市群道路移動源常規(guī)和非常規(guī)氣態(tài)污染物排放清單,研究了分車型、城市區(qū)域的機動車氣態(tài)污染物排放特征與分擔率,以期為決策者制定機動車污染防治措施提供依據(jù),并為城市群光化學煙霧污染防控提供支持.

    1 材料與方法

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

    關中城市群是以西安為中心,包括咸陽、寶雞、渭南、銅川等地級城市以及楊凌農(nóng)業(yè)示范區(qū),本研究以這5市(其中將楊凌區(qū)并入咸陽市一起考慮)為研究對象,關中城市群地理位置如圖1所示.

    圖1 關中城市群地理位置示意

    研究建立關中城市群道路移動源常規(guī)和非常規(guī)氣態(tài)污染物排放清單,主要包括機動車排放的CO、NO、NH3與SO2等4種常規(guī)氣態(tài)污染物和非甲烷碳氫化合物(NMHC)、甲醛(HCHO)、乙醛(CH3CHO)、丙烯醛(C3H4O)、1,3-丁二烯(C4H6)、苯(C6H6)、甲烷(CH4)和氧化亞氮(N2O)等8種非常規(guī)氣態(tài)污染物.研究對象為汽車、摩托車和低速汽車三大類道路移動排放源,不包括非道路移動源(如農(nóng)用機械、工程機械等)的排放.

    1.2 清單建立方法

    機動車污染物排放量的計算基于排放因子法,主要涉及三個參數(shù):機動車保有量、年均行駛里程和機動車排放因子.估算采用以下式(1)[26]:

    式中:EQ為污染物年排放總量,g;EF為型車的排放因子,g/(km·輛);P為型車的保有量,輛;VKT為型車的單車年均行駛里程,km.

    1.3 排放因子

    污染物排放因子反映了一個地區(qū)某種車型的污染物排放水平,獲得地區(qū)不同類型車輛單車污染物排放因子EF是進行道路移動源污染物排放清單計算的基礎.近年來美國環(huán)境保護署基于VSP分布開發(fā)了最新一代排放模型—MOVES模型,該模型能夠運用于國家、城市和路段等尺度的排放模擬研究[27],模擬精度較MOBILE模型等更高,其中增加了NMHC、NO、NO2、甲醛、乙醛、丙烯醛、苯等光化學煙霧前體物的計算,提高了計算結(jié)果與光化學煙霧反應化學機制的匹配性.郝艷召等[28-29]利用該模型研究了西安市機動車排放清單和輕型車顆粒物排放,邱兆文等[30]利用該模型模擬評估了平面交叉口區(qū)域的機動車排放污染,基于前期研究,本文利用關中城市群內(nèi)的實地調(diào)研和測試數(shù)據(jù)對MOVES模型進行本土化修正后作為道路移動源氣態(tài)污染物排放因子計算模型.論文對車齡分布、燃油信息、平均速度分布、道路類型分布、行駛里程時間分布、氣象信息、燃油和發(fā)動機分類、機動車保有量等8項具有區(qū)域性特征的參數(shù)進行了關中地區(qū)本土化修正,占模型可修正參數(shù)總數(shù)的80%,I/M制度和匝道比例2項參數(shù)采用模型默認值.

    1.3.1 車齡分布 模型中車齡分布參數(shù)是指研究區(qū)域內(nèi)在用機動車不同登記年代的車輛占同類型在用機動車保有量的百分比.論文通過統(tǒng)計年鑒、現(xiàn)場抽樣調(diào)查(調(diào)查有效樣本1332輛,其中小客車240輛,大客車450輛,載貨汽車642輛)獲得關中城市群道路機動車車齡分布參數(shù).其中城市群內(nèi)不同車齡載貨汽車分布比例見圖2.

    圖2 關中城市群載貨汽車車齡分布

    1.3.2 燃油信息 根據(jù)調(diào)查,關中地區(qū)機動車使用的燃油類型主要有汽油、柴油和壓縮天然氣(CNG)三種,車用燃油大部分為國III標準.本文參照國家標準《車用汽油》(GB17930-2011)[31]、《車用柴油》(GB19147-2009)[32]和《車用壓縮天然氣》(GB18047-2000)[33]各類燃油的組分指標,選取與MOVES模型給定的燃油組分中最接近的燃油配方,用我國燃油標準中的參數(shù)將模型參數(shù)進行修正,主要采用常規(guī)汽油(ID為10)和常規(guī)柴油(ID為20).

    1.3.3 平均速度和行駛里程時間分布 論文通過在測試車輛(客車20輛,貨車10輛)上安裝GPS設備(Multifunction GPS Data Logger),實時收集不同類型監(jiān)測車輛的瞬時速度、經(jīng)緯度、時間和海拔高度等數(shù)據(jù)信息,觀測時間覆蓋了工作日和周末,24h連續(xù)監(jiān)測,包含全部4種道路類型,試驗期間獲得有效樣本量為小型客車7161645s,大型客車1154592s,重型貨車198103s.通過跟蹤監(jiān)測獲得的GPS數(shù)據(jù)與GIS(Map info)解譯的道路地圖匹配,將CSV格式數(shù)據(jù)通過Access導入到GIS地圖,進而計算出平均速度分布比例、行駛里程(VMT)時間分布比例.其中行駛里程的周和小時分布比例見圖3和圖4.

    圖3 行駛里程周分布

    圖4 行駛里程日分布

    1.3.4 道路類型分布 模型中道路類型分布是指機動車在不同類型道路上行駛里程分布比例.論文將關中地區(qū)典型道路類型與MOVES模型中的道路類型進行匹配,其中快速路對應城市限制入口道路(如西安市二環(huán)、三環(huán)等)、非快速路對應城市非限制入口道路(如城市內(nèi)主干路、次干路、支路等)、高速路對應鄉(xiāng)村限制入口道路(如連霍高速等)、非高速路對應鄉(xiāng)村非限制入口道路(如一級公路、二級公路等).通過測試車輛GPS數(shù)據(jù)與GIS地圖匹配,按照車型-道路類型-車輛行駛里程的格式分類,統(tǒng)計分析出不同類型車輛在不同類型道路上的實際行駛距離,并按車型進行具體比例計算.關中城市群內(nèi)典型機動車在不同道路上的行駛里程分布情況,具體如圖5所示.

    1.3.5 氣象信息 根據(jù)陜西省氣象局西安基準站(E108°56′,N34°18′)的常規(guī)地面氣象觀測數(shù)據(jù),得到關中地區(qū)每月的平均溫度和濕度變化,具體見表1,相關數(shù)據(jù)輸入模型.

    圖5 不同類型道路上的行駛里程分布

    表1 關中地區(qū)月平均溫度和相對濕度

    1.3.6 燃油和發(fā)動機類型 關中地區(qū)不同機動車燃油類型占比情況如圖6所示.由圖6可以看出,小型客車絕大部分為汽油車,占比達到96.36%;中型客車和輕型貨車中汽油車比例也較高,分別為61.5%和65.6%;重型貨車絕大部分為柴油車,占比達到87.88%;中型貨車中柴油車略占主導,占比達到60.77%;而大型客車中汽油、柴油和其他燃油車基本各占1/3,相關數(shù)據(jù)輸入模型.

    1.4 機動車保有量

    根據(jù)西安、寶雞、銅川、咸陽(含楊凌)、渭南等城市歷年的統(tǒng)計年鑒[34-38]和陜西省環(huán)境監(jiān)測站的統(tǒng)計數(shù)據(jù),按照用途、載重量和軸數(shù)等參數(shù),將我國典型的機動車類型與MOVES模型中設定的車型進行匹配,得到2012年關中城市群各市的機動車保有量如表2所示.其中,載客汽車、載貨汽車詳細車型構(gòu)成情況如下,載客汽車中以小型客車為主,占比達97.05%,中型客車和大型客車的比例很低,分別為1.47%和1.48%.載貨汽車中輕型貨車最多,占52.94%,重型貨車比例次之為28.24%,中型貨車的比例最低占18.82%.

    圖6 不同類型汽車燃油類型比例

    表2 關中城市群機動車保有量(2012年)(輛)

    1.5 年均行駛里程

    圖7 各車型實調(diào)年均行駛里程

    通過對機動車檢測站、客運站和貨運站實地抽樣調(diào)查(調(diào)查有效樣本1332輛,其中小客車240輛,大客車450輛,載貨汽車642輛)獲得車輛行駛里程和車齡信息數(shù)據(jù),得到各車型年均行駛里程見圖7.對于中型客車、摩托車、低速汽車年均行駛里程,以中國環(huán)境科學研究院統(tǒng)計數(shù)據(jù)[39]和國家環(huán)境保護部《道路機動車排放清單編制技術指南(試行)》和《非道路移動污染源排放清單編制技術指南(試行)》中推薦數(shù)據(jù)作為參考值.

    2 結(jié)果與討論

    2.1 排放因子特征

    利用輸入?yún)?shù)本土化修正后的MOVES模型,對關中城市群內(nèi)小型客車、中型客車、大型客車、小型貨車、中型貨車、大型貨車、摩托車等7種車型氣態(tài)污染物排放因子進行計算.低速汽車(主要包括農(nóng)用三輪和四輪車)的NO和CO排放因子采用《非道路移動源排放清單編制技術指南(試行)》推薦值,NMHC、NH3和SO2等排放因子參考計算得到的輕型柴油貨車的值.分車型污染物排放因子見表3.

    由表3可以看出,中型貨車NO排放因子最高,是小客車的18倍之多,大型客車和重型貨車次之.摩托車和中型貨車CO排放因子較其他車型明顯高,分別為小客車CO排放因子的5倍和4倍,分析認為這主要與不同車型的燃料類型、車齡和車況等因素的差異有關.從機動車直接排放的NMHC來看,摩托車排放因子最高,其次為中型貨車,分別為小客車排放因子的15倍和4.5倍,摩托車直接排放的甲烷、1,3-丁二烯和苯的單車排放因子也明顯高于其他車輛,分別為相應小客車排放因子的4倍、9倍和9.2倍之多,這與摩托車其發(fā)動機燃燒不完全、車輛劣化嚴重等因素有關,提高摩托車單車污染控制技術水平、加強摩托車生產(chǎn)和在用車輛污染物排放限值達標控制、推行城市限摩措施[40]等,對關中城市中心區(qū)的有機污染物排放控制不容忽視.據(jù)研究[41-42],NH3已成為影響我國大城市空氣質(zhì)量的重要污染物,其能夠與城市大氣中高濃度的NO2、SO2和揮發(fā)性有機酸等物質(zhì)反應生成二次氣溶膠,會導致細粒子污染問題發(fā)生.機動車尾氣排放的NH3是大氣環(huán)境中氨氣和銨離子的重要來源之一[43].關中地區(qū)各種車型都會產(chǎn)生一定的NH3,其排放因子在0.014~0.024g/km之間,中型客車和中型貨車的NH3排放因子相對較高,研究表明汽油車排放的NH4+約有56.32%未被強酸根中和,是重要的活性NH4+排放源[44],機動車NH3排放將對城市環(huán)境空氣二次氣溶膠形成產(chǎn)生顯著影響.

    表3 機動車污染物排放因子(g/km)

    將本文研究得到的關中地區(qū)小客車(汽油) CO、NO和HC(NMHC+CH4)污染物排放因子與清華大學2008~2010年在北京、廣州和深圳實測的輕型汽油轎車排放因子進行比較,將甲醛、乙醛、苯等非常規(guī)氣態(tài)污染物排放因子與北京理工大學輕型汽油車臺架試驗排放因子進行比較,見表4.從表4中可以看出,關中地區(qū)機動車CO、NO和HC(NMHC+CH4)污染物的排放因子水平整體偏高,處于研究實測值歐Ⅱ和歐Ⅲ之間,非常規(guī)氣態(tài)污染物甲醛和苯的排放因子與北京理工大學實測汽油車排放因子比較一致,模擬得到的乙醛排放因子值較實測值偏低.

    表4 輕型汽油轎車排放因子結(jié)果對比(g/km)

    2.2 總量排放特征

    根據(jù)各城市不同類型機動車活動水平數(shù)據(jù)和計算得到的排放因子,按照式(1)計算得到2012年為基準年的關中城市群道路移動源污染物年排放總量清單,計算結(jié)果見表5.結(jié)果表明,關中城市群道路移動源排放的常規(guī)污染物中CO的年排放量最高,達到45.4萬t;其次是NO,排放量為8.19萬t,是2012年陜西省NO非電力排放量(54.68萬t)[47]的16.28%;SO2排放量相對較低,為0.42萬t,是2012年陜西省SO2非電力排放量(57.62萬t)[47]的0.73%; NH3的年排放量為0.1萬t,其中西安市機動車NH3排放量為0.054′104t是其2013年人為源大氣氨排放量[48]的1.14%,對西安市二次氣溶膠形成具有重要作用,隨著機動車保有量的增加其對城市NH3排放的貢獻量將不斷增大,城市機動車的NH3排放越發(fā)受到關注,將是城市機動車污染減排的重要對象之一.非常規(guī)污染物中NMHC為4.168萬t,典型有機污染物中甲烷年排放量最高,達到0.123萬t,其次是苯為0.09萬t,甲醛為0.057萬t,乙醛為0.027萬t,丙烯醛為0.004萬t,1,3-丁二烯為0.012萬t.隨著機動車排放非常規(guī)污染物的增加,城市大氣中NO和NHMC、醛類、烯烴類物質(zhì)等發(fā)生反應形成二次污染的耦合過程更加復雜,會使城市大氣污染防控策略制定更加艱難.

    表5 2012年氣態(tài)污染物年排放總量清單(′104t)

    2.3 城市分擔率

    由圖8可以看出,2012年各城市的污染物排放分擔率受機動車保有量、機動車活動水平影響明顯,不同城市排放分擔率差異較大.NO、CO、NH3和SO2等常規(guī)氣態(tài)污染物的排放量分擔率占比西安最高,分別達到49.88%、47.82%、55.56%和50%,這與西安市機動車保有量占關中地區(qū)總量的48.98%基本一致,其次是渭南,咸陽(含楊凌)、寶雞,銅川分擔率最低.非常規(guī)氣態(tài)污染物如NMHC、甲醛、乙醛、丙烯醛等的排放城市分擔率也呈現(xiàn)相同特點.分析認為,西安市、渭南市機動車排放的NO和有機污染物區(qū)域分擔率總體較高,從城市光化學煙霧形成潛勢來看,在相同氣象條件下,這兩個城市更易發(fā)生區(qū)域性光化學煙霧污染天氣.研究表明城市機動車排放污染物主要集中在水平高度1m左右的人體呼吸區(qū)域,中心城區(qū)因高密度路網(wǎng)結(jié)構(gòu)而成為典型的高強度排放區(qū)域[49],關中城市群建成區(qū)路網(wǎng)密度高、人口聚集,機動車污染物排放控制形勢更嚴峻.研究認為要減少關中城市群道路移動源污染物排放總量,著重加強西安市和渭南市機動車污染物排放控制(占總排放量的61.99%~75.47%)是必要的選擇.

    圖8 關中機動車氣態(tài)污染物排放的城市分擔率

    2.4 車型分擔率

    根據(jù)關中城市群道路移動源污染物排放總量,結(jié)合不同城市、不同車型污染物排放量,計算得到關中地區(qū)不同車型氣態(tài)污染物排放分擔率.由圖9可以看出,不同車型對不同污染物排放分擔率差別較大,其中重型貨車對NO的排放分擔率最高,達到33.85%,其次為中型貨車,小型客車和大型客車的分擔率分別占10%以上;小型客車和摩托車是CO的主要排放源,兩種車型排放CO占總排量的60%以上;重型貨車對SO2的排放分擔率最高,達到了31.31%,這與其以柴油發(fā)動機為主有直接關系,其次為小型客車、中型貨車和大型客車,排放分擔率分別為21.15%、15.0%和13.89%;小型客車對NH3的排放分擔率達到50.25%,相關研究中也指出小型客車是機動車氨排放的主要貢獻者[50],這主要與燃料的不完全燃燒及小客車安裝尾氣催化裝置等有關,其中尾氣催化裝置產(chǎn)生的NH3更多,可能存在氨逃逸等問題.研究結(jié)果表明,機動車排放常規(guī)氣態(tài)污染物減排控制中,應當以不同車型的總量貢獻為出發(fā)點,在NO減排中應加強重型貨車、中型貨車和大型客車的污染物控制,CO減排則需要重點考慮小客車和摩托車,NH3減排控制以小客車、重型貨車為主,SO2減排控制則重點以重型貨車、小型客車、中型貨車為主.

    圖10 關中機動車非常規(guī)氣態(tài)污染物排放分車型分擔率

    圖10是機動車排放有機污染物NMHC、甲烷、甲醛、1,3-丁二烯、丙烯醛、苯等分車型排放分擔率情況.由圖10可以看出,不同車型對不同有機污染物的貢獻差異較大,總體來看摩托車和小客車是有機污染物排放的主要車型,低速汽車、重型貨車和中型貨車在醛類物質(zhì)排放中貢獻顯著.摩托車對NMHC和甲烷的排放分擔率分別達到55.21%和30.97%,摩托車NMHC排放貢獻量是小客車的3.15倍,是重型貨車的7.72倍.小客車和摩托車對1,3-丁二烯和苯的排放分擔率高達68.93%和75.59%,這與小客車保有量高、摩托車燃燒不充分導致單車排放高等因素有關.重型貨車、中型貨車、低速汽車和摩托車是甲醛、丙烯醛等醛類物質(zhì)的重要來源,四類車型的總排放量占甲醛、丙烯醛年總排量的76.71%和78.14%.研究結(jié)果表明,機動車排放的NMHC和甲烷總量減排中應加強摩托車的有機污染物排放控制,甲醛、丙烯醛等醛類污染物減排中需重點對重型貨車、中型貨車、低速汽車和摩托車進行控制,1,3-丁二烯和苯污染物減排以控制小客車和摩托車為主.

    2.5 不確定性分析與排放清單對比

    建立道路移動源排放清單過程中需要采集的數(shù)據(jù)較多,各類數(shù)據(jù)的獲取方法也不同,不確定性主要來源于機動車活動水平數(shù)據(jù)缺失以及排放因子的本地化差異.本研究機動車保有量數(shù)據(jù)來自各城市統(tǒng)計年鑒,不確定性較小,但由于我國機動車統(tǒng)計口徑與MOVES模型的車型設定存在差異,車型匹配過程中可能會存在一定的誤差.另外,受限于區(qū)域外流動機動車數(shù)據(jù)的可得性,本研究未考慮關中城市群區(qū)域外注冊車輛的跨域流動帶來的排放,增大了排放清單的不確定性.本研究采用的年均行駛里程數(shù)據(jù)基于現(xiàn)場抽樣調(diào)查獲得,同時部分車型數(shù)據(jù)參考了相關文獻,可能無法精確反映關中城市群內(nèi)機動車年均行駛里程狀況.本研究基于MOVES模型模擬的排放因子綜合考慮了車齡分布、燃油信息、道路類型、平均速度、行駛里程、氣象信息、燃油和發(fā)動機類型等本地化因素,排放因子較為可靠.

    表6 關中地區(qū)機動車排放氣態(tài)污染物清單對比(′104t/a)

    目前,關于關中城市群道路移動源氣態(tài)污染物排放的研究較少.通過文獻調(diào)研收集總結(jié)了相關學者對關中典型城市機動車常規(guī)氣態(tài)污染物排放清單的研究結(jié)果,如表6所示.從表6可以看出,本研究結(jié)果與王靜晞等[51]、閆東杰等[52]、蘇航等[48]的研究結(jié)果在數(shù)量級上一致,但略有差異,CO和NH3的數(shù)據(jù)一致性較好,NO差異較大,這種差別來自于活動水平設定的不同以及排放因子選取的差異.從活動水平設定上看,本文重點考慮到了NO主要貢獻車型大客車和載貨汽車年均行駛里程對排放造成的影響,通過現(xiàn)場調(diào)查等方式獲取了相對準確的行駛里程數(shù)據(jù).從排放因子選取上看,本文基于輸入?yún)?shù)本地化修正后的MOVES模型模擬的排放因子與王靜晞等通過匯總國內(nèi)的相關研究獲得的排放因子存在著固有差異.

    3 結(jié)論

    3.1 在對MOVES模型參數(shù)進行關中本地化修正的基礎上,研究了常規(guī)氣態(tài)污染物、非常規(guī)氣態(tài)污染物的排放因子,其中CO排放因子最高,NO和有機污染物排放因子次之,SO2排放因子較低.各車型中中型貨車NO排放因子最高,其次是重型貨車和低速汽車.摩托車有機污染物排放因子最高,其次為中型貨車.將本研究得到的關中地區(qū)機動車CO、HC和NO常規(guī)污染物的排放因子水平與歐洲標準相比整體偏高,非常規(guī)污染物甲醛和苯的排放因子與國內(nèi)同類實驗研究結(jié)果比較一致,乙醛排放因子略低.

    3.2 2012年關中城市群道路移動源常規(guī)氣態(tài)污染物排放量CO為45.4萬t,NO為8.19萬t,SO2為0.42萬t,NH3為0.1萬t,主要集中在西安市(占關中城市群的50%).機動車排放非常規(guī)氣態(tài)污染物中有機污染物年排放量NMHC為4.168萬t,典型有機污染物年排放量中甲烷最高,其次是苯,分別為0.123萬t和0.09萬t.

    3.3 關中地區(qū)對NO和SO2的年排放量分擔率最高的車型為重型貨車,分別達到33.85%和31.31%;小型客車和摩托車是CO的主要排放源,占總排量的60%以上,小型客車對NH3的排放分擔率達到50.25%.關中地區(qū)有機污染物排放分擔率較高的車型是摩托車和小客車,其對NMHC的排放分擔率分別達到55.21%和17.55%.對甲醛、丙烯醛等醛類物質(zhì)貢獻較大的車型是重型貨車、中型貨車、低速汽車和摩托車,四者貢獻率之和超過75%.

    3.4 關中各城市的污染物排放分擔率與城市機動車保有量、機動車活動水平有關,總體呈現(xiàn)出西安最高、銅川最低.

    [1] Wang Y, Liu H W, Mao G Z, et al. Inter-regional and sectoral linkage analysis of air pollution in Beijing–Tianjin–Hebei (Jing-Jin-Ji) urban agglomeration of China [J]. Journal of Cleaner Production, 2017,165:1436-1444.

    [2] Du Y Y, Sun T S, Peng J, et al. Direct and spillover effects of urbanization on PM2.5concentrations in China's top three urban agglomerations [J]. Journal of Cleaner Production, 2018,190:72-83.

    [3] Cariolet J M, Colombert M, VuilletM, et al. Assessing the resilience of urban areas to traffic-related air pollution: Application in Greater Paris [J]. Science of The Total Environment, 2018,615:588-596.

    [4] Itahashi S, Uno I, Kim S. Seasonal source contributions of tropospheric ozone over East Asia based on CMAQ-HDDM [J]. Atmospheric Environment, 2013,70:204-217.

    [5] 程曉梅,劉永紅,陳泳釗,等.珠江三角洲機動車排放控制措施協(xié)同效應分析[J]. 中國環(huán)境科學, 2014,34(6):1599-1606. Cheng X M, Liu Y H, Cheng YZ, et al. A comparative co-benefit analysis of the implements of vehicle emissions control policy in Pearl River Delta [J]. China Environmental Science, 2014,34(6):1599-1606.

    [6] 李貝睿,劉 湛,尤翔宇,等.長株潭區(qū)域機動車尾氣排放清單及特征分析[J]. 環(huán)境科學與技術, 2016,39(11):167-173. Li B R, Liu Z, YouX Y, et al. Emission and characteristics of vehicle exhausts in Changsha-Zhuzhou-Xiangtan area of Hunan Province [J]. Environmental Science & Technology, 2016,39(11):167-173.

    [7] 楊 雯,王學軍,張倩茹.基于高精度清單的京津冀機動車控污政策研究[J]. 中國環(huán)境科學, 2018,38(8):2803-2810. Yang W, Wang X J, Zhang Q R. Vehicle pollution control policy study on the Beijing-Tianjin-Hebei region, based on a high-precision emission inventory [J]. China Environmental Science, 2018,38(8): 2803-2810.

    [8] 郎建壘,程水源,韓力慧,等.京津冀地區(qū)機動車大氣污染物排放特征[J]. 北京工業(yè)大學學報, 2012,38(11):1716-1723. Lang J L, Cheng S Y, Han L H, et al. Vehicular Emission Characteristics in the Beijing-Tianjin-Hebei (BTH) Region [J]. Journal of Beijing University of Technology, 2012,38(11):1716-1723.

    [9] 魯斯唯,胡清華,吳水平,等.海峽西岸經(jīng)濟區(qū)大氣污染物排放清單的初步估算[J]. 環(huán)境科學學報, 2014,34(10):2624-2634.Lu S W, Hu Q H, Wu S P, et al. Establishment of air pollutant emission inventory in the West Coast of Taiwan Strait [J]. Acta Scientiae Circumstantiae, 2014,34(10):2624-2634.

    [10] 李瑞芃,吳 琳,毛洪鈞,等.廊坊市區(qū)主要大氣污染源排放清單的建立[J]. 環(huán)境科學學報, 2016,36(10):3527-3534.Li R P, Wu L, Mao H J, et al. Establishment of an air pollutant emission inventory for Langfang urban area [J]. Acta Scientiae Circumstantiae, 2016,36(10):3527-3534.

    [11] 李新興,孫國金,王孝文,等.杭州市區(qū)機動車污染物排放特征及分擔率[J]. 中國環(huán)境科學, 2013,33(9):1684-1689. Li X X, Sun G J, Wang X W, et al. Vehicle exhausts emission characteristics and contributions in Hangzhou district [J]. China Environmental Science, 2013,33(9):1684-1689.

    [12] 唐 偉,楊 強,黃 成,等.基于大數(shù)據(jù)分析和IVE模型的杭州市機動車污染物排放變化特征研究[J]. 環(huán)境科學學報, 2018,38(1): 71-78. Tang W, Yang Q, Huang C, et al. Study on characteristics of pollutant emission from motor vehicles in Hangzhou based on large data analysis and IVE model [J]. Acta Scientiae Circumstantiae, 2018,38(1): 71-78.

    [13] 周子航,鄧 也,吳柯穎,等.成都市道路移動源排放清單與空間分布特征[J]. 環(huán)境科學學報, 2018,38(1):79-91.Zhou Z H, Deng Y, Wu K Y, et al. On-road mobile source emission inventory and spatial distribution characteristics in Chengdu [J]. Acta Scientiae Circumstantiae, 2018,38(1):79-91.

    [14] 張景文,咸 月,陳報章.南昌市移動源排放清單研究[J]. 環(huán)境科學學報, 2017,37(7):2449-2458. Zhang J W, Xian Y, Chen B Z. Study on the emission inventory from mobile source in Nanchang [J]. Acta Scientiae Circumstantiae, 2017, 37(7):2449-2458.

    [15] 孫世達,姜 巍,高衛(wèi)東.青島市機動車排放清單與空間分布特征[J]. 中國環(huán)境科學, 2017,37(1):49-59. Sun S D, Jiang W, Gao W D. Vehicle emission inventory and spatial distribution in Qingdao [J]. China Environmental Science, 2017,37(1): 49-59.

    [16] 唐 偉,何 平,楊 強,等.基于IVE模型和大數(shù)據(jù)分析的杭州市道路移動源主要溫室氣體排放清單研究[J]. 環(huán)境科學學報, 2018, 38(4):1368-1376. Tang W, He P, Yang Q, et al. Study on greenhouse gas emission inventory of road source in Hangzhou based on IVE model and large data analysis [J]. Acta Scientiae Circumstantiae, 2018,38(4):1368- 1376.

    [17] Seinfeld J H, Pandis S N. Atmospheric chemistry and physics: from air pollution to climate change [M]. 3rded. Hoboken, New Jersey: John Wiley & Sons, Inc., 2016.

    [18] 陶雙成,鄧順熙,劉 寧,等.機動車排放對城市光化學煙霧形成的影響[J]. 世界科技研究與發(fā)展, 2015,37(1):21-25. Tao S C, Deng S X, Liu N, et al. Influence of vehicle exhaust on production of cities photochemical smog [J]. World Sci-tech R&D, 2015,37(1):21-25.

    [19] Yang W D, Zhang Q Y, Wang J L, et al. Emission characteristics and ozone formation potentials of VOCs from gasoline passenger cars at different driving modes [J]. Atmospheric Pollution Research, 2018,9: 804-813.

    [20] Ling Z H, Zhao J, FanS J, et al. Sources of formaldehyde and their contributions to photochemical O3formation at an urban site in the Pearl River Delta, southern China [J]. Chemosphere, 2017,168:1293- 1301.

    [21] Yin S S, Huang Z J, Zheng J Y,et alCharacteristics of inorganic aerosol formation over ammonia-poor and ammonia-rich areas in the Pearl River Delta region, China [J]. Atmospheric Environment, 2018, 177:120-131.

    [22] Sellakumar A R, Snyder C A, Solomon J J, et alCarcinogenicity of formaldehyde and hydrogen chloride in rats [J]. Toxicology and Applied Pharmacology, 1985,81:401-406.

    [23] McGregor D B, Baan R A, Partensky C, et al. Evaluation of the carcinogenic risks to humans associated with surgical implants and other foreign bodies – a report of an IARC Monographs Programme Meeting [J]. European Journal Cancer, 2000,36:307-313.

    [24] 黃少妮,袁 媛,井 宇,等.陜西關中地區(qū)冬季一次重霾污染過程及氣象條件影響分析[J]. 干旱氣象, 2016,34(6):1037-1046.Huang S N, Yuan Y, Jing Y, et al. Characteristics of a severe haze pollution process in winter of 2013 and meteorological cause in Guanzhong of Shaanxi [J]. Journal of Arid Meteorology, 2016,34(6): 1037-1046.

    [25] 陜西日報.關中平原地區(qū)整體空氣污染較重臭氧污染需警惕[EB/OL]. http://news.xiancity.cn/system/2018/05/30/030576902.shtml, 2018-5-30. Shaanxi Daily. The air pollution in Guanzhong Plain is serious, Ozone pollution needs to be vigilant [EB/OL]. http://news.xiancity.cn/system/ 2018/05/30/030576902.shtml, 2018-5-30.

    [26] 環(huán)境保護部.道路機動車大氣污染物排放清單編制技術指南(試行) [EB/OL]. http://www.zhb.gov.cn/gkml/hbb/bgg/201501/t20150107_ 293955.htm, 2014-12-31. Ministry of Environmental Protection. Technical guidelines for compiling air pollutant emission inventory of vehicles (trial implementation) [EB/OL]. http://www.zhb.gov.cn/gkml/hbb/bgg/ 201501/t20150107_293955.htm, 2014-12-31.

    [27] USEPA. User’s guide for MOVES (motor vehicle emission simulator) 2014 [R]. Washington, D. C.: USEPA, 2014.

    [28] 郝艷召,鄧順熙,邱兆文,等.基于MOVES模型的西安市機動車排放清單研究[J]. 環(huán)境污染與防治, 2017,39(3):227-231.Hao Y Z, Deng S X, Qiu Z W, et al. Vehicle emission inventory for Xi'an based on MOVES model [J]. Environmental Pollution and Control, 2017,39(3):227-231.

    [29] 郝艷召,鄧順熙,邱兆文,等.基于MOVES的輕型車顆粒物排放來源和特征分析[J]. 環(huán)境工程學報, 2015,9(8):3915-3922.Hao Y Z, Deng S X, Qiu Z W, et al. Particle sources and characteristics of light duty vehicle based on MOVES model [J]. Chinese Journal of Environmental Engineering, 2015,9(8):3915-3922.

    [30] 邱兆文,張廣昕,郝艷召.基于運行模式分布的平面交叉口汽車排放評估 [J]. 鄭州大學學報(工學版), 2014,35(6):86-90. Qiu Z W, Zhang G X, Hao Y Z.Estimation for vehicles emissions of at-grade intersections based on operating mode [J]. Journal of Zhengzhou University (Engineering Science), 2014,35(6):86-90.

    [31] GB17930-2011 《車用汽油》[S]. GB17930-2011 《Vehicle Gasoline》[S].

    [32] GB19147-2013 《車用柴油(Ⅳ)》[S]. GB19147-2013 《Vehicle Diesel (IV)》[S].

    [33] GB18047-2000 《車用壓縮天然氣》[S]. GB18047-2000 《Compressed Natural Gas for Vehicles》[S].

    [34] 西安市統(tǒng)計局.西安市2000~2012統(tǒng)計年鑒[R]. 西安:西安市統(tǒng)計局, 2001-2013. Xi'an Statistical Bureau. Xi'an Statistical Yearbook 2000~2012 [R]. Xi'an:Xi'an Statistical Bureau, 2001-2013.

    [35] 寶雞市統(tǒng)計局.寶雞市2000~2012統(tǒng)計年鑒[R]. 寶雞:寶雞市統(tǒng)計局, 2001-2013. Baoji Statistical Bureau. Baoji Statistical Yearbook 2000~2012 [R]. Baoji: Baoji Statistical Bureau, 2001-2013.

    [36] 銅川市統(tǒng)計局.銅川市2000~2012統(tǒng)計年鑒[R]. 銅川:銅川市統(tǒng)計局, 2001-2013. Tongchuan Statistical Bureau. Tongchuan Statistical Yearbook 2000-2012 [R]. Tongchuan: Baoji Statistical Bureau, 2001-2013.

    [37] 渭南市統(tǒng)計局.渭南市2000~2012統(tǒng)計年鑒[R]. 渭南:渭南市統(tǒng)計局, 2001-2013. Weinan Statistical Bureau. Weinan Statistical Yearbook 2000~2012 [R]. Weinan: Weinan Statistical Bureau, 2001-2013.

    [38] 咸陽市統(tǒng)計局.咸陽市2000~2012統(tǒng)計年鑒[R]. 咸陽:咸陽市統(tǒng)計局, 2001-2013. Xianyang Statistical Bureau. Xianyang Statistical Yearbook 2000~2012 [R]. Xianyang: Xianyang Statistical Bureau, 2001-2013.

    [39] 何立強,胡京南,謝淑霞,等.2010年中國機動車CH4和N2O排放清單[J]. 環(huán)境科學研究, 2014,27(1):28-35. He L Q, Hu J N, Xie S X, et al. CH4and N2O emission inventory for motor vehicles in China in 2010 [J]. Research of Environmental Sciences, 2014,27(1):28-35.

    [40] 程曉梅,劉永紅,徐偉嘉,等.佛山限摩對摩托車出行特征及減排的影響分析[J]. 環(huán)境科學與技術, 2013,36(6L):400-404. Cheng X M, Liu Y H, Xu W J, et al. Analysis of urban motorcycle travel characteristics and emissions reduction of motorcycle restriction in Foshan [J]. Environmental Science & Technology, 2013,36(6L): 400-404.

    [41] 李麗平,石金輝,李非非,等.青島大氣中HNO3、HNO2、NH3及PM2.5中氮組分的濃度特征和氣-粒平衡關系[J]. 環(huán)境科學學報, 2014, 34(11):2869-2877. Li L P, Shi J H, Li F F, et al. Concentrations of nitric acid, nitrous acid, ammonia gases and partner nitrogen ions in PM2.5and gas-particle partitioning analysis in Qingdao [J]. Acta Scientiae Circumstantiae, 2014,34(11):2869-2877.

    [42] Behera S N, Betha R, Balasubramanian R. Insights into chemical coupling among acidic gases, ammonia and secondary inorganic aerosols [J]. Aerosol and Air Quality Research, 2013,13(4):1282- 1296.

    [43] 楊 俊.廣州市機動車尾氣中氨氣的排放因子研究 [D]. 廣州:暨南大學, 2011. Yang J. Studyonammoniaemissionfactorsof vehiclesexhaustin Guangzhoucity [D]. Guangzhou: Jinan University, 2011.

    [44] 石愛軍,馬俊文,耿春梅,等.北京市機動車尾氣排放PM10組分特征研究[J]. 中國環(huán)境監(jiān)測, 2014,30(4):44-49.Shi A J, Ma J W, Geng C M, et al. Characteristics of chemical composition of particulate matter (PM10) from Beijing vehicle [J]. Environmental Monitoring in China, 2014,30(4):44-49.

    [45] Huo H, Yao Z L, Zhang Y Z, et alOn-board measurements of emissions from light-duty gasoline vehicles in three mega-cities of China [J]. Atmospheric Environment, 2012,49:371-377.

    [46] 游秋雯,葛蘊珊,尤可為,等.汽油車非常規(guī)污染物排放特性研究[J]. 環(huán)境科學, 2009,30(2):335-341.You Q W, Ge Y S, You K W, et al. Unregulated Emissions from the Gasoline Vehicle [J]. Environmental Science, 2009,30(2):335-341.

    [47] 綠色和平組織.江蘇、湖北、陜西、四川四省大氣污染防治政策評估報告[EB/OL]. http://www.greenpeace.org.cn/4-provinces- against-air-pollution/, 2016-6-15. Greenpeace. Assessment report on air pollution control policies in Jiangsu, Hubei, Shaanxi and Sichuan provinces [EB/OL]. http://www. greenpeace.org.cn/4-provinces-against-air-pollution/, 2016-6-15.

    [48] 蘇 航,閆東杰,黃學敏,等.西安市人為源大氣氨排放清單及特征[J]. 環(huán)境科學, 2016,37(11):4117-4123. Su H, Yan D J, Huang X M, et al. Inventory and characteristics of anthropogenic ammonia emission in Xi'an [J]. Environmental Science, 2016,37(11):4117-4123.

    [49] 王孝文,田偉利,張清宇.杭州市機動車污染物排放清單的建立[J]. 中國環(huán)境科學, 2012,32(8):1368-1374. Wang X W, Tian W L, Zhang Q Y. Development of motor vehicles emission inventory in Hangzhou [J]. China Environmental Science, 2012,32(8):1368-1374.

    [50] 張 雙,王阿婧,張增杰,等.北京市人為源氨排放清單的初步建立 [J]. 安全與環(huán)境學報, 2016,16(2):242-245. Zhang S, Wang A J, Zhang Z J, et al. On creating an anthropogenic ammonia emission inventory in capital Beijing [J]. Journal of Safety and Environment, 2016,16(2):242-245.

    [51] 王靜晞,曹國良,韓 蕾,等.關中地區(qū)人為源大氣污染物排放清單研究[J]. 安全與環(huán)境學報, 2015,15(5):282-286. Wang J X, Cao G L, Han L, et al. Probe to the anthropogenic air pollutant emission inventory in Guanzhong Region, Shaanxi [J]. Journal of Safety and Environment, 2015,15(5):282-286.

    [52] 閆東杰,蘇 航,黃學敏,等.西安市人為源揮發(fā)性有機物排放清單及研究[J]. 環(huán)境科學學報, 2017,37(2):446-452. Yan D J, Su H, Huang X M, et al. VOCs emission inventory of anthropogenic sources in Xi'an [J]. Acta Scientiae Circumstantiae, 2017,37(2):446-452.

    Vehicle emission characteristics of gaseous pollutants in Guanzhong urban agglomeration.

    TAO Shuang-cheng1*, DENG Shun-xi2, HAO Yan-zhao3, GAO Shuo-han1, XIONG Xin-zhu1, KONG Ya-ping1

    (1.China Academy of Transportation Sciences, Beijing 100029, China;2.School of Environment Science and Engineering, Chang’an University, Xi’an 710064, China;3.School of Automobile,Chang’an University, Xi’an 710064,China)., 2019,39(2):542~553

    This research calculated the emission factors of conventional and unconventional vehicle-emitted gas pollutants, based on the MOVES model with all parameters corrected and simulated to the situation of Guanzhong urban agglomeration. And the motor vehicle pollutant emission inventoryin 2012 of the research area was established, including the contributions from various motor vehicle types in each part of the area.The results showed that the amounts of conventional vehicle-emitted gas pollutants were 81.9kt, 454k, 1k and 4.2k tons of NOx, CO, NH3and SO2, respectively. While the amounts of unconventional gas pollutants were 41.68kt, 0.57kt, 0.27kt, 0.04kt, 0.12kt, 0.95kt, 1.23kt and 0.04kt of NMHC, formaldehyde, acetaldehyde, propylene aldehyde, 1, 3-butadiene and benzene, CH4and N2O, respectively. In addition, the shares of major vehicle-emitted gas pollutants in this area ranking from highest to lowestwere Xi’an(50%), Weinan (23%), Xianyang (including Yanglin) (12%), Baoji(10%), Tongchuan(5%), respectively. Furthermore, the shares of different types of vehicle varied significantly, among whichthe main source of NOxemission wereheavy trucks and medium trucks withcontributionaround 33.85% and 21.21%, respectively. While heavy trucks contributed the most to the emission of sulfur dioxide (~31.31%) and aldehydes(~30%). Specifically, small-size buses and motorcycles covered the most to the emission of CO, NMHC, 1,3-butadiene, benzene and methane, with 32.86%, 17.55%, 26.64%, 26.45%, 38.85% for the former and 32.64%, 55.21%, 43.29%, 49.04%, 30.97% for the later, respectively. And small-size buses and heavy trucks were the maincontributors to NH3emission with 49.5% and 31.31%, respectively.

    vehicleexhaust;gaseous pollutants;emission characteristics;MOVES;Guanzhong urban agglomeration

    X511

    A

    1000-6923(2019)02-0542-12

    陶雙成(1981-),男,甘肅白銀人,副研究員,博士,主要研究方向為交通環(huán)境與安全技術.發(fā)表論文10余篇.

    2018-07-17

    中央級公益性科研院所基本科研業(yè)務費項目(20160609, 20180614);國家自然科學基金項目(51478045);陜西省科技統(tǒng)籌創(chuàng)新工程項目(2012KTZB03-01-04)

    * 責任作者, 副研究員, tsc504@126.com

    猜你喜歡
    小客車氣態(tài)貨車
    貨車
    幼兒畫刊(2023年12期)2024-01-15 07:06:14
    北京小客車搖號又創(chuàng)新高3076人搶一個指標
    ISO/TS 19880-1:2016氣態(tài)氫加注站第1部分一般要求標準解讀
    氣態(tài)燃料發(fā)動機相關發(fā)明專利(三)
    氣態(tài)燃料發(fā)動機相關發(fā)明專利(二)
    全省實行小客車總量調(diào)控管理
    氣態(tài)燃料發(fā)動機相關發(fā)明專利
    貨車也便捷之ETC新時代!——看高速公路貨車ETC如何實現(xiàn)
    北京今年小客車指標配額15萬
    汽車縱橫(2017年3期)2017-03-18 23:17:14
    推貨車里的愛
    學與玩(2017年6期)2017-02-16 07:07:24
    中国国产av一级| 久久综合国产亚洲精品| 乱码一卡2卡4卡精品| 成人影院久久| 精品人妻偷拍中文字幕| 久热久热在线精品观看| 免费在线观看完整版高清| 卡戴珊不雅视频在线播放| 国产国语露脸激情在线看| 欧美人与性动交α欧美精品济南到 | 久久这里有精品视频免费| 大香蕉久久网| 免费大片黄手机在线观看| 18禁观看日本| 精品视频人人做人人爽| 亚洲,欧美精品.| 啦啦啦视频在线资源免费观看| 国产成人精品无人区| 丁香六月天网| 亚洲婷婷狠狠爱综合网| 国产黄色免费在线视频| 丝瓜视频免费看黄片| 亚洲欧洲日产国产| 亚洲中文av在线| 欧美xxxx性猛交bbbb| 免费观看a级毛片全部| 巨乳人妻的诱惑在线观看| 咕卡用的链子| 波多野结衣一区麻豆| 日韩在线高清观看一区二区三区| 老司机影院成人| 极品人妻少妇av视频| www.av在线官网国产| 国产亚洲欧美精品永久| 国产高清三级在线| 久久久久久伊人网av| 亚洲在久久综合| 天堂俺去俺来也www色官网| 精品久久蜜臀av无| 国产探花极品一区二区| 亚洲五月色婷婷综合| 精品国产乱码久久久久久小说| 欧美变态另类bdsm刘玥| 欧美人与性动交α欧美软件 | 国产男女内射视频| 久久人人爽人人片av| 女人精品久久久久毛片| www.av在线官网国产| 天天躁夜夜躁狠狠躁躁| 亚洲精品久久午夜乱码| 日韩人妻精品一区2区三区| 日韩制服丝袜自拍偷拍| 国产精品免费大片| 午夜福利在线观看免费完整高清在| 精品99又大又爽又粗少妇毛片| 国产日韩一区二区三区精品不卡| 日本免费在线观看一区| 久久人人97超碰香蕉20202| 亚洲av电影在线观看一区二区三区| 国产成人精品久久久久久| 狠狠精品人妻久久久久久综合| 成人18禁高潮啪啪吃奶动态图| 国产精品国产三级专区第一集| 少妇高潮的动态图| 国产精品国产三级国产专区5o| 欧美国产精品一级二级三级| 成年美女黄网站色视频大全免费| av在线app专区| 国产亚洲精品久久久com| 精品人妻偷拍中文字幕| 99热全是精品| 黑人猛操日本美女一级片| 男女下面插进去视频免费观看 | 人妻 亚洲 视频| 日韩免费高清中文字幕av| 亚洲欧洲国产日韩| 国产亚洲av片在线观看秒播厂| 色视频在线一区二区三区| 国产一区二区三区av在线| 精品亚洲乱码少妇综合久久| 日本猛色少妇xxxxx猛交久久| 99视频精品全部免费 在线| 色婷婷久久久亚洲欧美| 国产成人av激情在线播放| av免费观看日本| 国产一区有黄有色的免费视频| 黄色毛片三级朝国网站| 精品熟女少妇av免费看| 亚洲av免费高清在线观看| 在线观看三级黄色| 亚洲一码二码三码区别大吗| 久久久a久久爽久久v久久| 亚洲,一卡二卡三卡| 久久99精品国语久久久| 欧美激情国产日韩精品一区| 亚洲美女搞黄在线观看| 成人无遮挡网站| 大码成人一级视频| 日韩中文字幕视频在线看片| 1024视频免费在线观看| 国产精品久久久av美女十八| 侵犯人妻中文字幕一二三四区| xxxhd国产人妻xxx| 亚洲成国产人片在线观看| 国产极品天堂在线| 国产欧美亚洲国产| 国产高清三级在线| 人妻人人澡人人爽人人| 色哟哟·www| 中文天堂在线官网| 深夜精品福利| 高清在线视频一区二区三区| 色婷婷久久久亚洲欧美| 午夜免费鲁丝| 十分钟在线观看高清视频www| 97在线人人人人妻| 男人爽女人下面视频在线观看| 又粗又硬又长又爽又黄的视频| 亚洲欧美成人精品一区二区| av在线app专区| 麻豆乱淫一区二区| 91精品伊人久久大香线蕉| av女优亚洲男人天堂| www日本在线高清视频| 五月天丁香电影| 日本av手机在线免费观看| 桃花免费在线播放| 成人毛片a级毛片在线播放| 欧美日本中文国产一区发布| 天天躁夜夜躁狠狠躁躁| 免费日韩欧美在线观看| 亚洲精品国产色婷婷电影| 日韩中字成人| 国产成人精品久久久久久| 在线天堂中文资源库| 欧美+日韩+精品| 国产欧美亚洲国产| 免费高清在线观看日韩| 99国产综合亚洲精品| 亚洲一区二区三区欧美精品| 纵有疾风起免费观看全集完整版| 成人毛片a级毛片在线播放| 欧美精品人与动牲交sv欧美| 在线看a的网站| 中文字幕另类日韩欧美亚洲嫩草| 精品熟女少妇av免费看| 国产成人精品无人区| 大香蕉久久成人网| 国产一区二区在线观看日韩| 成人手机av| 免费人妻精品一区二区三区视频| www.av在线官网国产| 久久人人97超碰香蕉20202| 欧美成人午夜精品| 另类精品久久| 久久久久久人人人人人| 香蕉精品网在线| 久久热在线av| 18禁观看日本| 国产乱人偷精品视频| 国产亚洲一区二区精品| 国产无遮挡羞羞视频在线观看| 欧美日韩亚洲高清精品| 超碰97精品在线观看| 久热久热在线精品观看| 十八禁高潮呻吟视频| 免费黄频网站在线观看国产| 亚洲丝袜综合中文字幕| kizo精华| 91精品三级在线观看| 中文字幕另类日韩欧美亚洲嫩草| 丝袜喷水一区| 哪个播放器可以免费观看大片| 青春草国产在线视频| 又黄又粗又硬又大视频| av天堂久久9| 日本猛色少妇xxxxx猛交久久| 国产欧美另类精品又又久久亚洲欧美| 街头女战士在线观看网站| 你懂的网址亚洲精品在线观看| av有码第一页| 国产探花极品一区二区| 99国产精品免费福利视频| 日本91视频免费播放| 丝袜脚勾引网站| 黄色毛片三级朝国网站| 制服人妻中文乱码| 亚洲欧洲精品一区二区精品久久久 | 亚洲欧美成人综合另类久久久| 香蕉国产在线看| 亚洲精品日韩在线中文字幕| 伊人亚洲综合成人网| 男女免费视频国产| 久久99一区二区三区| 午夜福利网站1000一区二区三区| 男的添女的下面高潮视频| 一级片免费观看大全| 国产黄色免费在线视频| 如日韩欧美国产精品一区二区三区| 亚洲美女视频黄频| 精品久久蜜臀av无| 人人澡人人妻人| 国产在视频线精品| 久久精品国产综合久久久 | 久久精品国产亚洲av涩爱| 国产免费现黄频在线看| 免费观看av网站的网址| 满18在线观看网站| 亚洲av综合色区一区| 黄色 视频免费看| 男人舔女人的私密视频| 999精品在线视频| 久久韩国三级中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 18在线观看网站| 国产精品一区二区在线观看99| 亚洲欧美精品自产自拍| 十八禁网站网址无遮挡| 秋霞伦理黄片| 青春草国产在线视频| 亚洲精品中文字幕在线视频| www日本在线高清视频| 亚洲 欧美一区二区三区| 一级毛片 在线播放| 国产片特级美女逼逼视频| xxxhd国产人妻xxx| √禁漫天堂资源中文www| 七月丁香在线播放| 久久久久精品久久久久真实原创| 午夜福利网站1000一区二区三区| 免费看光身美女| 精品人妻在线不人妻| 精品亚洲乱码少妇综合久久| 亚洲精品日本国产第一区| 亚洲国产精品成人久久小说| 99久国产av精品国产电影| 99久国产av精品国产电影| 美女福利国产在线| 欧美另类一区| 亚洲精品国产色婷婷电影| 国产激情久久老熟女| 大码成人一级视频| 亚洲一级一片aⅴ在线观看| 国产成人免费观看mmmm| 亚洲美女视频黄频| 十分钟在线观看高清视频www| 免费大片18禁| 亚洲精品久久午夜乱码| 一二三四中文在线观看免费高清| 欧美另类一区| 国产精品国产三级国产av玫瑰| 这个男人来自地球电影免费观看 | 成年女人在线观看亚洲视频| 女人被躁到高潮嗷嗷叫费观| 51国产日韩欧美| 曰老女人黄片| 亚洲成人手机| xxxhd国产人妻xxx| 美女福利国产在线| 久久 成人 亚洲| 97人妻天天添夜夜摸| 久久久久国产精品人妻一区二区| 国产老妇伦熟女老妇高清| 51国产日韩欧美| 韩国av在线不卡| 欧美成人精品欧美一级黄| 久久久精品免费免费高清| 久久av网站| 亚洲,欧美,日韩| 久久女婷五月综合色啪小说| 久久久久久久久久久免费av| 婷婷色麻豆天堂久久| 亚洲国产av新网站| 国产永久视频网站| 欧美激情 高清一区二区三区| 亚洲精品av麻豆狂野| 2018国产大陆天天弄谢| 日韩成人伦理影院| av有码第一页| 国产熟女欧美一区二区| 亚洲一码二码三码区别大吗| 在线天堂最新版资源| 交换朋友夫妻互换小说| 国产极品粉嫩免费观看在线| 777米奇影视久久| 久久久国产一区二区| 亚洲av电影在线观看一区二区三区| 看非洲黑人一级黄片| 久久久久久伊人网av| 激情五月婷婷亚洲| 男女高潮啪啪啪动态图| 欧美精品av麻豆av| 狂野欧美激情性bbbbbb| 欧美成人精品欧美一级黄| 久久国产精品男人的天堂亚洲 | a级毛片在线看网站| www.色视频.com| 日本-黄色视频高清免费观看| 夫妻性生交免费视频一级片| 欧美97在线视频| 日韩制服丝袜自拍偷拍| 女人被躁到高潮嗷嗷叫费观| 建设人人有责人人尽责人人享有的| 国产精品99久久99久久久不卡 | 两性夫妻黄色片 | 五月开心婷婷网| 国产69精品久久久久777片| 18禁观看日本| 亚洲av在线观看美女高潮| 久久精品国产亚洲av天美| 久久久久久人人人人人| 人妻少妇偷人精品九色| 久久精品夜色国产| av不卡在线播放| 男人爽女人下面视频在线观看| 久久国内精品自在自线图片| 欧美国产精品va在线观看不卡| 日本欧美国产在线视频| 在线观看国产h片| 久久精品人人爽人人爽视色| 国产精品 国内视频| 亚洲综合色惰| 欧美少妇被猛烈插入视频| av卡一久久| 亚洲五月色婷婷综合| 99re6热这里在线精品视频| 日韩熟女老妇一区二区性免费视频| 九九在线视频观看精品| 国产黄色免费在线视频| 精品卡一卡二卡四卡免费| 亚洲图色成人| 精品少妇黑人巨大在线播放| 日本黄色日本黄色录像| 街头女战士在线观看网站| 亚洲av国产av综合av卡| 免费高清在线观看视频在线观看| 国产免费福利视频在线观看| 飞空精品影院首页| 亚洲国产精品成人久久小说| 免费黄色在线免费观看| 亚洲精品第二区| 国产69精品久久久久777片| 亚洲国产成人一精品久久久| videos熟女内射| 大香蕉久久成人网| 日韩中字成人| 久久久久久伊人网av| 国产成人精品婷婷| 一级毛片电影观看| 久久综合国产亚洲精品| 亚洲情色 制服丝袜| 高清不卡的av网站| 欧美日韩国产mv在线观看视频| 国产成人免费无遮挡视频| 国产免费一区二区三区四区乱码| 久久精品国产综合久久久 | 高清不卡的av网站| 夜夜骑夜夜射夜夜干| 在线观看免费高清a一片| 大香蕉97超碰在线| 亚洲综合精品二区| 国产 精品1| 天堂俺去俺来也www色官网| 少妇熟女欧美另类| 狠狠精品人妻久久久久久综合| 纯流量卡能插随身wifi吗| 伊人久久国产一区二区| 又黄又爽又刺激的免费视频.| tube8黄色片| 丝瓜视频免费看黄片| 午夜影院在线不卡| 国产午夜精品一二区理论片| 欧美日本中文国产一区发布| 免费av不卡在线播放| 成人毛片a级毛片在线播放| 亚洲精品aⅴ在线观看| 在线观看免费日韩欧美大片| 制服诱惑二区| 色吧在线观看| 免费av中文字幕在线| 婷婷色综合大香蕉| 久久这里只有精品19| 日韩制服骚丝袜av| 视频区图区小说| www.色视频.com| 制服人妻中文乱码| 丝袜美足系列| 久久韩国三级中文字幕| 99久国产av精品国产电影| a 毛片基地| 亚洲国产成人一精品久久久| 婷婷色av中文字幕| tube8黄色片| 久久久精品免费免费高清| 亚洲精品美女久久久久99蜜臀 | 亚洲经典国产精华液单| 日本黄色日本黄色录像| 精品国产一区二区三区久久久樱花| 国产有黄有色有爽视频| 人体艺术视频欧美日本| 久久99蜜桃精品久久| 亚洲国产av影院在线观看| 午夜免费鲁丝| 婷婷色综合大香蕉| 久久久国产一区二区| 在线免费观看不下载黄p国产| tube8黄色片| 亚洲,欧美,日韩| 99精国产麻豆久久婷婷| 少妇人妻 视频| 美女福利国产在线| 97精品久久久久久久久久精品| 国产日韩欧美亚洲二区| 伊人久久国产一区二区| 少妇的逼水好多| 全区人妻精品视频| 国产精品无大码| 国产精品欧美亚洲77777| 国产亚洲午夜精品一区二区久久| 亚洲欧美日韩卡通动漫| 自线自在国产av| 精品少妇黑人巨大在线播放| 香蕉丝袜av| 日本色播在线视频| 激情五月婷婷亚洲| 亚洲高清免费不卡视频| 少妇人妻精品综合一区二区| 中文字幕制服av| 亚洲美女搞黄在线观看| kizo精华| 亚洲欧洲国产日韩| 婷婷色麻豆天堂久久| 国产成人精品婷婷| 国产老妇伦熟女老妇高清| 狠狠精品人妻久久久久久综合| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 9191精品国产免费久久| 亚洲精品一区蜜桃| 亚洲av日韩在线播放| 成人国语在线视频| 亚洲精品日本国产第一区| 久久精品国产亚洲av天美| 赤兔流量卡办理| 岛国毛片在线播放| 伊人亚洲综合成人网| 亚洲国产欧美在线一区| 啦啦啦视频在线资源免费观看| 激情视频va一区二区三区| 亚洲欧美色中文字幕在线| 久久av网站| 这个男人来自地球电影免费观看 | 母亲3免费完整高清在线观看 | 精品午夜福利在线看| 成人无遮挡网站| 深夜精品福利| 最黄视频免费看| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看免费高清a一片| 久热久热在线精品观看| 国产精品久久久久久久电影| 久久久久视频综合| 久久午夜综合久久蜜桃| 国内精品宾馆在线| 极品少妇高潮喷水抽搐| 成人免费观看视频高清| 国产成人91sexporn| av在线播放精品| 秋霞在线观看毛片| 日韩一区二区三区影片| 新久久久久国产一级毛片| 亚洲精品日韩在线中文字幕| av天堂久久9| 久久久久精品性色| 在线观看一区二区三区激情| 精品熟女少妇av免费看| 七月丁香在线播放| 在线观看国产h片| 人妻少妇偷人精品九色| 18禁动态无遮挡网站| 涩涩av久久男人的天堂| 18+在线观看网站| 精品国产一区二区久久| 国产免费视频播放在线视频| 最近的中文字幕免费完整| 少妇的逼水好多| 交换朋友夫妻互换小说| 午夜福利视频精品| 九九爱精品视频在线观看| 久久av网站| 欧美精品一区二区免费开放| av天堂久久9| 久久99热6这里只有精品| 成人综合一区亚洲| 夫妻性生交免费视频一级片| 久久久久久伊人网av| 亚洲婷婷狠狠爱综合网| 新久久久久国产一级毛片| 宅男免费午夜| 成人国产av品久久久| 黄片无遮挡物在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲人与动物交配视频| 80岁老熟妇乱子伦牲交| 边亲边吃奶的免费视频| 国产福利在线免费观看视频| 国产亚洲最大av| 在线观看人妻少妇| av免费观看日本| 成年人午夜在线观看视频| 亚洲中文av在线| 色婷婷av一区二区三区视频| 亚洲国产色片| 90打野战视频偷拍视频| 黄色怎么调成土黄色| 国产色爽女视频免费观看| 国产精品 国内视频| 99久久精品国产国产毛片| 国产在线视频一区二区| 美女国产视频在线观看| 欧美日韩av久久| 亚洲精品粉嫩美女一区| 大片电影免费在线观看免费| 亚洲成人免费电影在线观看| 国产一区二区三区在线臀色熟女 | 午夜影院日韩av| 这个男人来自地球电影免费观看| 脱女人内裤的视频| 韩国精品一区二区三区| 黄色视频,在线免费观看| 亚洲欧洲精品一区二区精品久久久| 在线观看www视频免费| 黄色视频,在线免费观看| 操美女的视频在线观看| 在线观看舔阴道视频| 午夜福利一区二区在线看| 精品熟女少妇八av免费久了| 韩国av一区二区三区四区| 日本欧美视频一区| 免费看十八禁软件| 亚洲午夜理论影院| 国产高清激情床上av| 久久午夜亚洲精品久久| 美女 人体艺术 gogo| 成年人午夜在线观看视频| 亚洲情色 制服丝袜| www.自偷自拍.com| 操出白浆在线播放| 久久国产精品影院| 岛国毛片在线播放| 操美女的视频在线观看| 欧美激情久久久久久爽电影 | 久久精品aⅴ一区二区三区四区| 最新的欧美精品一区二区| 丰满人妻熟妇乱又伦精品不卡| 国产亚洲一区二区精品| 久久国产精品大桥未久av| 欧美人与性动交α欧美精品济南到| 搡老熟女国产l中国老女人| 国产精品国产av在线观看| 久久这里只有精品19| 欧美日韩亚洲高清精品| 青草久久国产| 亚洲成人免费电影在线观看| 最新的欧美精品一区二区| 麻豆成人av在线观看| 日韩人妻精品一区2区三区| 免费在线观看视频国产中文字幕亚洲| 欧美日韩亚洲综合一区二区三区_| 久久久久精品国产欧美久久久| 日本五十路高清| 午夜精品在线福利| 欧美不卡视频在线免费观看 | 国产高清视频在线播放一区| 男人操女人黄网站| 日本黄色视频三级网站网址 | 欧美日韩视频精品一区| 岛国在线观看网站| 国内毛片毛片毛片毛片毛片| 久久久久久久久免费视频了| 国产不卡一卡二| 亚洲精品国产精品久久久不卡| 亚洲情色 制服丝袜| 国产av一区二区精品久久| 国产高清国产精品国产三级| 午夜福利,免费看| 国产精品电影一区二区三区 | 久久久久久久久免费视频了| 亚洲熟妇熟女久久| 在线国产一区二区在线| 99国产综合亚洲精品| av不卡在线播放| 国产成人精品久久二区二区91| 久久久久国产精品人妻aⅴ院 | 韩国av一区二区三区四区| 国产精品av久久久久免费| 亚洲专区字幕在线| 久久中文字幕一级| 高清欧美精品videossex| 久久人妻熟女aⅴ| 欧美精品啪啪一区二区三区| 国产成人精品久久二区二区91| 欧美激情久久久久久爽电影 | 热99国产精品久久久久久7| 亚洲av第一区精品v没综合| 好看av亚洲va欧美ⅴa在| 国产伦人伦偷精品视频| av电影中文网址| 91成年电影在线观看| 亚洲欧美一区二区三区久久| 日韩欧美免费精品| 成人永久免费在线观看视频| 日韩视频一区二区在线观看| netflix在线观看网站| 90打野战视频偷拍视频| 黑人操中国人逼视频|