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

    2013~2014年北京市NO2時(shí)空分布研究

    2016-06-07 06:21:07程念亮李云婷張大偉劉保獻(xiàn)北京市環(huán)境保護(hù)監(jiān)測中心大氣顆粒物監(jiān)測技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室北京0008北京師范大學(xué)水科學(xué)研究院北京00875中國環(huán)境科學(xué)研究院環(huán)境基準(zhǔn)與風(fēng)險(xiǎn)評估國家重點(diǎn)實(shí)驗(yàn)室北京0002北京市環(huán)境保護(hù)局北京0008
    中國環(huán)境科學(xué) 2016年1期
    關(guān)鍵詞:大氣污染

    程念亮,李云婷,張大偉*,陳 添,董 欣,王 欣,郇 寧,劉保獻(xiàn),閆 賀,孟 凡(.北京市環(huán)境保護(hù)監(jiān)測中心,大氣顆粒物監(jiān)測技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室,北京 0008;2.北京師范大學(xué)水科學(xué)研究院,北京 00875;3.中國環(huán)境科學(xué)研究院,環(huán)境基準(zhǔn)與風(fēng)險(xiǎn)評估國家重點(diǎn)實(shí)驗(yàn)室,北京 0002;.北京市環(huán)境保護(hù)局,北京0008)

    ?

    2013~2014年北京市NO2時(shí)空分布研究

    程念亮1,2,3,李云婷1,張大偉1*,陳 添4,董 欣1,王 欣1,郇 寧1,劉保獻(xiàn)1,閆 賀1,孟 凡2,3**(1.北京市環(huán)境保護(hù)監(jiān)測中心,大氣顆粒物監(jiān)測技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室,北京 100048;2.北京師范大學(xué)水科學(xué)研究院,北京 100875;3.中國環(huán)境科學(xué)研究院,環(huán)境基準(zhǔn)與風(fēng)險(xiǎn)評估國家重點(diǎn)實(shí)驗(yàn)室,北京 100012;4.北京市環(huán)境保護(hù)局,北京100048)

    摘要:根據(jù)2013~2014年北京市NO2監(jiān)測數(shù)據(jù),對比分析了全年及重污染日NO2時(shí)空分布特征.結(jié)果表明:2013年NO2平均濃度為56μg/m3,2014年北京市NO2年均濃度為56.7μg/m3.年均及重污染日NO2月均濃度均呈波浪型分布,日變化呈雙峰型分布;空間分布上北部及西部山區(qū)NO2濃度明顯低于中心城區(qū)及南部地區(qū).NO2濃度與PM2.5、CO、NO呈正相關(guān)關(guān)系,與O3、OX無明顯相關(guān)性;全年NO2光解速率峰值平均在0.105/min左右,重污染日光解速率峰值平均在0.026/min左右;全年及重污染期間,氮氧化速率分別為0.142±0.061、0.190±0.036;高濃度NO2既有利于O3生成,又對重污染的形成起到了促進(jìn)作用;重污染日特定條件下北京市NO2的兩種轉(zhuǎn)化機(jī)制以轉(zhuǎn)化為NO3-過程為主導(dǎo).經(jīng)計(jì)算2000~2014年北京市機(jī)動(dòng)車的保有量與NO2濃度的相關(guān)系數(shù)R為-0.84,機(jī)動(dòng)車N Ox排放量對北京市NO2濃度的變化有顯著的影響.

    關(guān)鍵詞:NO2;時(shí)空分布;北京;NOR;PM2.5;光解速率

    * 責(zé)任作者, 高級工程師, zhangdawei@bjmemc.com.cn

    ** 責(zé)任作者, 研究員, mengfan@craes.org.cn

    NO2既是O3前體物,又可光化學(xué)反應(yīng)生成NO3-,是我國空氣質(zhì)量評價(jià)常規(guī)監(jiān)測項(xiàng)目,大氣中NO2濃度升高影響著生態(tài)環(huán)境及人體健康[1].北京作為我國的首都,僅6000km2面積承載了超過2000×104人、500×104輛機(jī)動(dòng)車以及大量的生產(chǎn)、服務(wù)活動(dòng),氮氧化物排放總量較高,NO2防治形勢較為嚴(yán)峻[2-3].

    近年來,國內(nèi)外對中國中東部特別是北京地區(qū)NO2濃度監(jiān)測及影響因素研究做了較多工作,Verders等[4]利用GOME資料研究NO2垂直柱濃度發(fā)現(xiàn)我國中東部地區(qū)污染較為嚴(yán)重;Richrer 等[5]利用SCIAMACHY資料發(fā)現(xiàn)NO2垂直柱濃度在中國東部以及香港地區(qū)增長顯著;Kang等[6]發(fā)現(xiàn)沙塵天氣的出現(xiàn)會(huì)增加大氣中的 NO2的濃度;Zhang等[7]發(fā)現(xiàn)受中國人為源排放影響的NO2高值區(qū)范圍增加且華北地區(qū)NOx排放量增加了 133%;王郭臣等[8]發(fā)現(xiàn)天津市NO2濃度受氣流影響較大;田賀忠等[9]發(fā)現(xiàn)輻合風(fēng)場的頻繁發(fā)生及不利地形導(dǎo)致了太行山東麓NOx濃度不斷積累;李令軍等[10]分析了1999~2010年北京市大氣中NO2濃度,發(fā)現(xiàn)近年來NO2濃度有所反彈;尉鵬等[11]發(fā)現(xiàn)降水量、邊界層氣壓場是影響NO2濃度季節(jié)分布主要因素;江文華等[12]發(fā)現(xiàn)北京市大氣NO2污染多出現(xiàn)在1月、2月、11月、12 月,與其特有的地形特征、氣候特征密切相關(guān);王躍啟等[13]發(fā)現(xiàn)華北NO2垂直柱濃度分布與經(jīng)濟(jì)發(fā)展有較好的一致性且煤炭產(chǎn)區(qū)和消費(fèi)區(qū)對NO2濃度分布有著重要影響.

    這些研究多基于單個(gè)站點(diǎn)或某個(gè)研究時(shí)段,且所使用的研究數(shù)據(jù)時(shí)空分辨率較粗[14-15].本文利用北京市環(huán)境保護(hù)監(jiān)測中心2013~2014年全市35個(gè)站點(diǎn)NO2逐時(shí)濃度監(jiān)測資料對比分析了年均和重污染狀態(tài)下NO2時(shí)空分布特征及成因,以期為管理部門提供準(zhǔn)確、及時(shí)、全面的信息.

    1 材料與方法

    1.1 站點(diǎn)分布

    北京市位于東經(jīng)115.7°E~117.4°E,北緯39.4°N~41.6°N,地處華北平原西北端,臨近半沙漠化地帶邊緣,地形為簸箕型,三面環(huán)山,平均海拔43.5m,山地一般海拔1000~1500m,地形較不利于污染物擴(kuò)散.國土面積16410.54km2,62%為山區(qū),平原面積僅6000多km2;森林資源總量偏低,全市資源覆蓋率37.6%,大氣自凈功能較弱.位于北緯40°地區(qū),屬溫帶大陸性季風(fēng)氣候,四季分明,夏季高溫多雨,冬季寒冷干燥,近10年年降水量平均不足450mm,又易受沙塵影響,年均降水的80%集中在夏季6、7、8三個(gè)月[16].

    1984年北京市建立了我國第一套環(huán)境空氣質(zhì)量自動(dòng)監(jiān)測系統(tǒng),隨著城市規(guī)模發(fā)展以及城市空間布局的變化,為科學(xué)、合理地反映全市環(huán)境空氣質(zhì)量狀況,站點(diǎn)不斷增加;截至目前,北京市的空氣質(zhì)量地面自動(dòng)監(jiān)測網(wǎng)絡(luò)由35個(gè)站點(diǎn)(圖1)組成,包括1個(gè)城市清潔對照點(diǎn),23個(gè)城市環(huán)境評價(jià)點(diǎn),6個(gè)區(qū)域背景傳輸點(diǎn),5個(gè)交通污染監(jiān)控點(diǎn);35個(gè)站點(diǎn)分布在城六區(qū)(東城、西城、朝陽、海淀、豐臺(tái)、石景山);西北部(昌平、延慶);東北部(懷柔、密云、平谷、順義);東南部(通州、大興、亦莊);西南部(房山、門頭溝),覆蓋所有區(qū)縣

    [17].本研究污染物空間分布使用的是35個(gè)監(jiān)測站點(diǎn),其它濃度分析使用的是東四、官園、天壇、萬壽西宮、奧體中心、農(nóng)展館、萬柳、古城、順義、昌平、懷柔共11個(gè)監(jiān)測點(diǎn)的平均值.

    1.2 監(jiān)測儀器

    北京市污染物監(jiān)測數(shù)據(jù)為北京市環(huán)境保護(hù)監(jiān)測中心自動(dòng)監(jiān)測網(wǎng)絡(luò)逐時(shí)濃度數(shù)據(jù)(http://zx.bjmemc.com.cn/).SO2采用紫外熒光法進(jìn)行監(jiān)測,分析儀為美國T E-43C型分析儀,儀器的標(biāo)定系統(tǒng)由111型零空氣發(fā)生器 (TEI 111)和146型多點(diǎn)校準(zhǔn)儀 (TEI,146C)組成.O3的監(jiān)測均采用美國熱電環(huán)境儀器公司Thermo Fisher生產(chǎn)的49C紫外光度法O3分析儀,O3儀器分析儀最低檢測限1×10-9(體積分?jǐn)?shù)),精度1×10-9(體積分?jǐn)?shù)),零漂0.4%/24h,跨漂±1%/24h、±2%/7d.采用Thermo Fisher 42C化學(xué)發(fā)光NO-NO2-NOx分析儀監(jiān)測NO和NO2濃度,該分析儀最低檢測限0.05×10-9(體積分?jǐn)?shù)),零漂小于0.025×10-9/24h,跨漂:±1%/24h.采用基于氣體濾波紅外吸收法Thermo Fisher 48C分析儀測量CO濃度,儀器最低檢測限為40×10-9(體積分?jǐn)?shù)),零漂小于100×10-9/24h,跨漂為±1%/24h.采用基于微量震蕩天平法的Thermo 1405F系列儀器監(jiān)測PM2.5濃度,24h連續(xù)采樣.采用大氣PM2.5水溶性離子在線測量儀(URG-9000B)對水溶性離子Na+、Mg2+、Ca2+、K+、NH4+、Cl-、SO42-、NO3-、PO43-進(jìn)行在線觀測,其化學(xué)成分分析系統(tǒng)為美國DIONEX公司的Dionex ICS-90 型離子色譜儀,時(shí)間分辨率為60min.各儀器操作流程嚴(yán)格按照《環(huán)境空氣質(zhì)量自動(dòng)監(jiān)測技術(shù)規(guī)范》(HJ/T 193-2005)[18]進(jìn)行,設(shè)備由技術(shù)人員定期檢查并及時(shí)維護(hù)保養(yǎng).氣象資料為北京市觀象臺(tái)地面觀測數(shù)據(jù)(http://cdc.cma.gov.cn/),儀器為荷蘭WAISALA公司的WXT520氣象觀測儀.

    圖1 監(jiān)測站分布和分類Fig.1 Distribution and classification of observation sites

    2 結(jié)果與討論

    2.1 時(shí)間分布

    2.1.1 年際變化 自2000年以來,北京市主要大氣污染物濃度總體呈下降趨勢.SO2濃度下降趨勢較為明顯,累計(jì)降低69%左右.2000~2006年SO2濃度呈波動(dòng)式下降,從80μg/m3下降至53μg/m3;2007年起逐年下降,2014年達(dá)到最低值22μg/m3.

    2000~2014年全市PM10濃度總體呈明顯下降趨勢,累計(jì)下降28%,但近年下降速度有所放緩.其中2000~2006年呈波動(dòng)式變化,2002年濃度最高, 為166μg/m3;2007年以來持續(xù)下降,2013年最低,為108μg/m3;2014年濃度有所回升,為116μg/m3.

    多年來北京市PM2.5整體呈現(xiàn)下降趨勢,2013年北京環(huán)保監(jiān)測中心正式對外發(fā)布PM2.5監(jiān)測數(shù)據(jù),其中歷史數(shù)據(jù)為恢復(fù)計(jì)算的結(jié)果[19],由圖2可知,2008年前PM2.5處于高污染水平,濃度在101~133μg/m3之間,2008年后PM2.5年均濃度下降幅度顯著,其中2013年P(guān)M2.5年均濃度為89.5μg/m3,2014年P(guān)M2.5年均濃度為85.9μg/m3.

    2008~2014年全市CO濃度在1.29~1.67mg/m3之間波動(dòng),其中2010年濃度最高(1.67mg/m3),2014年濃度最低(1.29mg/m3),2014 年CO濃度與2013年(1.47mg/m3)相比,同比降低約12%,其中2008年以來CO濃度年變化率為-34.7μg/m3,累計(jì)降低約10%.

    NO2年均濃度在波動(dòng)中總體呈下降趨勢,其中2000~2002年NO2年均濃度有所上升,2002年達(dá)到最高值76μg/m3;2003~2008年明顯下降,2008年濃度達(dá)到最低值49μg/m3,奧運(yùn)減排措施起到明顯效果;2009~2014年NO2年均濃度趨于平穩(wěn),整體在52~58μg/m3間波動(dòng)但濃度有小幅上升;2000~2014年NO2濃度累計(jì)降幅20%.近兩年,2013年NO2濃度為56μg/m3,2014年NO2濃度為56.7μg/m3;2014年與2013年相比上升約1.3%,但仍超國家二級標(biāo)準(zhǔn)(40μg/m3)42%.

    圖2 2000~2014年北京市主要大氣污染物年際變化Fig.2 Annual variations of major air pollutants in Beijing during 2000~2014

    2.1.2 月變化 根據(jù)環(huán)境保護(hù)部《環(huán)境空氣質(zhì)量指數(shù)(AQI)技術(shù)規(guī)定》(HJ 633-2012)[20]分級方法,0

    圖3為2013~2014年北京市年均及重污染日NO2逐月分布.由圖3可知,北京市NO2的月均濃度整體呈波浪形分布,其中2013年1月份的濃度最高,達(dá)到88.06μg/m3;10月份同樣出現(xiàn)一個(gè)NO2濃度的峰值,峰值濃度為65.77μg/m3;2014年10月份的濃度最高,濃度為72.75μg/m3,2 月份次之,濃度為66.84μg/m3.1~2月份的濃度較高,主要受采暖季燃煤的影響[21];10月份峰值可能與特殊的污染源排放和氣象擴(kuò)散條件有關(guān),2014年10月北京市出現(xiàn)4次典型重污染過程,污染過程中各項(xiàng)污染物濃度明顯上升[22];7月份和8月份的 NO2月均濃度最低,最低濃度為40.21μg/m3.分季節(jié)來看,2013年春季(3~5月)、夏季(6~8月)、秋季(6~8 月)、冬季(12、1、2月)平均濃度分別為55.51、43.61、58.94、68.23μg/m3;2014年春季、夏季、秋季、冬季平均濃度分別為57.53、42.49、63.95、63.05μg/m3,總體呈現(xiàn)出秋、冬季大于春、夏季的濃度分布特征,與李令軍等[10]、王英等[23]對北京NO2濃度監(jiān)測及衛(wèi)星反演的結(jié)果基本一致.

    圖3 北京市NO2月變化曲線Fig.3 Monthly average concentration of NO2in Beijing

    對于重污染日平均濃度月分布而言,2013年NO2月平均濃度在54.00~117.92μg/m3之間;2014 年NO2月平均濃度在39.60~120.50μg/m3之間;2月、11~12月重污染日月平均濃度同比明顯上升,濃度增大13.43~33.00μg/m3.重污染日NO2濃度上升除受區(qū)域源排放(燃煤、機(jī)動(dòng)車)影響外,高濃度污染帶傳輸也起著重要作用[24].

    2.1.3 日變化 圖4為2013~2014年北京市年均及重污染日NO2日變化曲線,全年來看,NO2濃度呈現(xiàn)雙峰型分布,在早上07:00~8:00出現(xiàn)第1個(gè)峰值,在夜間的20:00~22:00之間出現(xiàn)第2個(gè)峰值,第2個(gè)峰值濃度明顯高于第1個(gè)峰值,早上的峰值與交通早高峰比較一致,夜間的峰值可能與夜間NO和O3反應(yīng)生成NO2,而NO2光解反應(yīng)停滯,濃度逐漸積累有關(guān)[25].2014年與2013年NO2日變化相比,02:00~11:00時(shí)濃度明顯偏高,可能與柴油車數(shù)量增大造成的較高排放有關(guān)[26],這些車夜間才允許進(jìn)入北京市五環(huán)內(nèi),造成夜間時(shí)段氮氧化物排放量較大.對于重污染日變化而言,18:00至次日04:00維持較高濃度,下午濃度沒有明顯的降低;與年均濃度相比, 重污染日平均濃度上升約46.62%~112.71%.

    圖4 北京市NO2日變化曲線Fig.4 Diurnal variation of NO2in Beijing

    2.2 空間分布

    圖5為運(yùn)用Kriging插值法[27]得出的2013~2014年均及重污染日北京市NO2空間分布.由圖5可知,NO2年均濃度的空間分布結(jié)構(gòu)具有鮮明的特征,北部及西部山區(qū)NO2濃度水平明顯低于中心城區(qū)及南部地區(qū),懷柔、密云、延慶年均濃度在30μg/m3左右,中心城區(qū)及南部地區(qū)濃度水平在50μg/m3以上.對于重污染日NO2濃度空間分布而言,北部地區(qū)及西部門頭溝等地明顯低于中心城區(qū)及南部地區(qū),前者平均濃度在45μg/m3以下,后者在80~110μg/m3之間.重污染日中心城區(qū)及南部地區(qū)濃度明顯上升與不利氣象條件造成局地積累及污染傳輸有關(guān)[28].2014年與2013年NO2年均濃度相比,在沿八達(dá)嶺高速存在一個(gè)西北方向伸展的NO2增長區(qū),增大約0.60~2.10μg/m3,這可能與西北方向八達(dá)嶺高速、京藏高速的交通污染排放特別是貨運(yùn)交通污染排放有關(guān);2014年與2013年NO2重污染日NO2平均濃度相比, NO2濃度增長區(qū)主要集中在城六區(qū)及南部地區(qū),平均增大2.0μg/m3以上.由于北部及門頭溝等地區(qū)監(jiān)測站點(diǎn)相對較少,插值的濃度分布結(jié)果具有一定的不確定性[29].本文插值后統(tǒng)計(jì)的面積不確定性主要來自:①"簸箕狀"的特殊地形及植被對污染物濃度的影響;②插值方法誤差,插值受采樣點(diǎn)范圍、采樣點(diǎn)密度等參數(shù)影響.

    圖5 北京市NO2空間分布Fig.5 Spatial distribution of NO2in Beijing

    2.3 NO2反周末效應(yīng)

    圖6 北京市城區(qū)工作日和周末NO2、O3小時(shí)濃度變化曲線Fig.6 Diurnal variation of NO2,O3in weekdays and on weekends in urban district in Beijing

    圖6為周末及工作日O3、NO2小時(shí)濃度的日變化曲線,Cleveland等[15]最早提出O3“周末效應(yīng)”概念,即周末O3濃度值高于工作日濃度的現(xiàn)象.由圖6可以看出,北京市存在O3“周末效應(yīng)”,日均值周末比工作日高5%左右,這可能由于周末與工作日機(jī)動(dòng)車行駛數(shù)量造成的 NOx排放強(qiáng)度的明顯變化導(dǎo)致[30].而北京市 NO2濃度存在明顯的“反周末效應(yīng)”,周末NO2濃度比工作日偏低7%左右,與黃志新[31]在上海的觀測結(jié)果相似.北京市 NO2“反周末效應(yīng)”與人類在周末和工作日的作息安排差異密切相關(guān),另外特殊氣象條件也可能有著一定的作用[32].

    2.4 相關(guān)性分析

    圖7為2013~2014年北京市年均及重污染日觀測期間的 NO2、O3、NO、NOx和大氣氧化劑 OX(OX=NO2+O3)小時(shí)濃度的日變化規(guī)律.由圖7可知,對于年均濃度而言,NO呈現(xiàn)出明顯的雙峰變化,NO2在一天之中變化幅度較小;NO在8:00、23:00各有1個(gè)峰值,夜間峰值可能與柴油車較高排放有關(guān),早上的峰值可能與上班時(shí)段機(jī)動(dòng)車排放有關(guān)[33];NO2夜間開始積累并在21:00左右出現(xiàn)1個(gè)峰值,可能因?yàn)橐归gNO和O3反應(yīng)生成NO2,而NO2光解反應(yīng)停滯,濃度逐漸積累; NO2與O3之和OX反映了大氣的氧化程度,整體來看受太陽輻射和光化學(xué)反應(yīng)強(qiáng)度影響,白天O3不斷生成,午后濃度最高,O3日變化呈現(xiàn)規(guī)律的單峰型特征;OX也呈現(xiàn)出單峰型日變化特征,夜間 OX 主要以NO2為主,而正午前后 NO2與O3在OX中的比例大致相當(dāng)[34];PM2.5日變化曲線仍呈現(xiàn)雙峰型分布,在上午10:00或11:00出現(xiàn)第1個(gè)峰值,在夜間的21:00或22:00出現(xiàn)第2個(gè)峰值,第2個(gè)峰值濃度明顯高于第1個(gè)峰值,表明夜間北京市城區(qū)的PM2.5更容易積累或二次生成,可能與夜間的排放源和化學(xué)反應(yīng)有關(guān)

    [35];CO濃度也呈現(xiàn)雙峰型分布,在早上09:00左右出現(xiàn)第1個(gè)峰值,在夜間的23:00左右之間出現(xiàn)第2個(gè)峰值;下午濃度較低,這與下午溫度高、日照強(qiáng)、邊界層高度較高、垂直擴(kuò)散較好有關(guān).對于重污染日平均濃度而言,NO2、CO濃度雙峰型分布不明顯,早上的峰值消失且變化較為平緩;由于重污染日主要集中在冬季,溫度較低,光化學(xué)反應(yīng)較弱,重污染日O3濃度明顯降低, 但PM2.5濃度較高,這些污染物日變化特征及原因還需進(jìn)一步研究.

    空氣中NO3-來源既有一次排放的,也有可能是經(jīng)過大氣光化學(xué)反應(yīng)生成的二次氣溶膠.大氣中的NOx在大氣中經(jīng)大氣光化學(xué)氧化形成NO3-(R1),在NH3存在的情況下生成NH4NO3, 而大氣中的NO3-多以NH4NO3的形式存在,在高溫作用下,NH4NO3可分解為NO3-和NH3(R2).氣相反應(yīng)一般在晴朗、干燥的天氣下發(fā)生;而HNO3主要通過氣相氧化反應(yīng)生成,NO3-是PM2.5的重要化學(xué)組分,因此 NO2濃度對PM2.5濃度有顯著影響.

    圖7 北京市近地層大氣中O3、NO、NO2、NOx、OX、PM2.5濃度小時(shí)變化Fig.7 Diurnal variations of O3, NO, NO2, NOx, OX and PM2.5concentrations near the ground in Beijing

    對2013~2014年NO2小時(shí)濃度與 PM2.5、O3、NO、NOx、CO、OX小時(shí)濃度進(jìn)行相關(guān)性分析發(fā)現(xiàn),NO2與PM2.5、NO、NOx、CO呈正相關(guān),與O3、OX呈負(fù)相關(guān).全年及重污染日NO2濃度與PM2.5濃度相關(guān)系數(shù)分別為0.83、0.86,均達(dá)到 0.70 以上,呈顯著性相關(guān)關(guān)系.全年及重污染日NO2濃度與CO為正相關(guān),相關(guān)系數(shù)分別為0.91、0.46;全年及重污染日NO2濃度與NO為正相關(guān),相關(guān)系數(shù)分別為0.69、0.12;全年NO2濃度與O3濃度呈現(xiàn)顯著的負(fù)相關(guān)性(R=-0.77),重污染日NO2濃度與O3濃度相關(guān)系數(shù)明顯下降(R為-0.04);全年NO2濃度與OX濃度呈負(fù)相關(guān)性(R=-0.60),重污染日NO2濃度與OX濃度相關(guān)系數(shù)明顯下降(R=-0.32).

    重污染日NO2濃度PM2.5濃度相關(guān)性遠(yuǎn)大于與O3濃度的相關(guān)性,原因可能為:重污染日多發(fā)生在冬季,溫度較低,濕度較大,植被VOCs排放量較小,生成O3濃度較低,而重污染日NO2濃度PM2.5濃度相關(guān)系數(shù)有所上升表明重污染日特殊的氣象條件通過R10引發(fā)NO2的二次轉(zhuǎn)化,生成PM2.5濃度中重要的化學(xué)組分NO3-.此外,重污染日NO2與其它污染物相關(guān)性比全年下降表明重污染期間NO2濃度的影響因素增多,特殊的氣象條件可能起著十分重要的作用[36].

    2.5 NO2光解速率與氮氧化速率

    NO2是O3前體物,特定條件下高濃度NO2既有利于O3生成,也引發(fā)NO2向NO3-二次轉(zhuǎn)化,重污染日北京市NO2的兩種轉(zhuǎn)化機(jī)制,是否以轉(zhuǎn)化為NO3-過程為主導(dǎo),研究計(jì)算NO2光解速率與氮氧化速率來說明.

    圖8 北京市NO2光解速率與氮氧化速率日變化曲線Fig.8 Diurnal variation of Photolysis rate of NO2and NOR in Beijing

    一般來說,白天太陽輻射強(qiáng)度變化不大時(shí)NO、NO2和O3這3種物質(zhì)濃度會(huì)達(dá)到光化學(xué)穩(wěn)定態(tài),服從經(jīng)典的光化學(xué)穩(wěn)定態(tài)方程[37],如R4式中K1為方程R12的反應(yīng)速率,是溫度的函數(shù),常用的關(guān)系式為R13.利用實(shí)測的污染物濃度和氣象資料再結(jié)合方程(R3)、(R4)可以求得J2的值(圖8a),由圖8a可知全年及重污染日11:00~15:00 J2最大,全年J2峰值平均在0.105min-1左右,重污染日J(rèn)2峰值平均在0.026min-1左右,一天中重污染日J(rèn)2明顯低于全年平均;本研究計(jì)算的NO2光解速率與安俊琳等[26]、Wang等[38]在北京的計(jì)算結(jié)果相近.

    表1 不同地區(qū)NOR比較Table 1 Comparison of NOR at different districts

    研究表明[39-40]重污染過程中PM2.5主要組分NO3-、SO42-、NH4+質(zhì)量濃度之和占PM2.5濃度的60%以上,3種無機(jī)組分離子主要由氣態(tài)前體物的二次反應(yīng)生成.為了評價(jià)NO3-的生成速率,研究采用氮轉(zhuǎn)化率NOR=[NO3-]/[NO3-]+[NO2]來表征,其中各物質(zhì)均采用摩爾濃度(mol/L).全年及重污染期間,NOR值分別為0.142±0.061、0.190±0.036,全年平均NOR值在一天之中于14:00左右到達(dá)峰值0.166,而重污染日NOR值無明顯日變化,24h變化幅度較小.但重污染日NOR值明顯高于全年平均水平,而重污染日NO2光解速率明顯低于年均值.重污染日特殊的低溫、高濕氣象條件使氮氧化物等一次污染氣體在大氣顆粒物表面通過非均相反應(yīng)轉(zhuǎn)化生硝酸鹽氣溶膠, 白天HNO3主要來源與NOx和CxHy在對流層的光化學(xué)反應(yīng),而夜間主要來自N2O5的反應(yīng),因而硝酸鹽的二次生成加重了重污染的水平,一定程度上表明重污染日北京市NO2的兩種轉(zhuǎn)化機(jī)制中以轉(zhuǎn)化為NO3-過程為主導(dǎo).本研究與其他學(xué)者[41-43]在北京的研究結(jié)果具有一定的可比性.不同城市的NOR結(jié)果[44-47]不盡相同,可能原因?yàn)棰倬值亍^(qū)域污染源排放量的差異造成污染物濃度不同②氣象條件如風(fēng)速、溫度、濕度差異③產(chǎn)業(yè)結(jié)構(gòu)、工業(yè)布局、能源利用的間接影響④地形及下墊面的影響.

    2.6 NO2來源分析

    人類活動(dòng)和自然行為都有可能產(chǎn)生二氧化氮,二氧化氮的自然源主要是閃電和微生物生命活動(dòng),人為污染來源包括汽車尾氣排放、工業(yè)排放、飛機(jī)排放、輪船排放等.近年來隨著汽車保有量量的飛速增長,機(jī)動(dòng)車尾氣排放對NO2的影響日益顯著.張興贏等[48]發(fā)現(xiàn)中國的人口分布與NO2的濃度分布有密切的相關(guān),人口多、人為活動(dòng)頻繁的地區(qū)都是NO2的高值區(qū);陳煥盛等[49]利用嵌套網(wǎng)格空氣質(zhì)量預(yù)報(bào)模式系統(tǒng)NAQPMS研究2006年亞運(yùn)時(shí)段廣州的NO2時(shí)空變化及來源發(fā)現(xiàn)來自本地的NO2的月均貢獻(xiàn)率分別為89.5%;周丹等[50]利用CALPUFF模擬研究了2005年東莞市NO2來源后發(fā)現(xiàn)NO2本地源貢獻(xiàn)占77%,其中機(jī)動(dòng)車貢獻(xiàn)62%;Wang等[51]研究表明2006年11月4~6日在北京中非合作論壇期間機(jī)動(dòng)車的尾氣排放對北京地區(qū)對流層NO2濃度的貢獻(xiàn)量約為70%.國內(nèi)外對PM2.5濃度及組分來源解析研究較多,對NO2來源解析研究較少.本文受數(shù)值模型及排放清單等條件制約,研究利用北京市機(jī)動(dòng)車的保有量、氮氧化物排放量與NO2濃度年變化的相關(guān)關(guān)系初步分析了機(jī)動(dòng)車尾氣排放對NO2濃度的影響(圖9,數(shù)據(jù)來自文獻(xiàn)[16,20]),由圖9可知,2014年底北京市機(jī)動(dòng)車保有量達(dá)到了559.1萬輛,較2006年凈增94%; 2006年以來,北京市機(jī)動(dòng)車保有量增長較快,但此后的2011年,北京市開始實(shí)施汽車總量控制措施和搖號(hào)上牌制度,增長速度放緩.經(jīng)計(jì)算2000~2014年北京市機(jī)動(dòng)車的保有量與NO2濃度的相關(guān)系數(shù)R為-0.84,在機(jī)動(dòng)車保有量增加的情況下,NO2年均濃度總體呈下降趨勢,盡管近兩三年濃度有所回升.

    圖9 2000~2014年北京市NO2年均濃度與機(jī)動(dòng)車保有量Fig.9 Annual variation of NO2concentration and vehicle numbers during 2000~2014 in Beijing

    研究表明[52]北京市NOx主要來源于機(jī)動(dòng)車及工業(yè)排放,機(jī)動(dòng)車已成為北京市氮氧化物排放的首要污染源.本文根據(jù)各車型的年均行駛里程,平均運(yùn)行車速和各車型平均排放因子和保有量,得出了2006~2014年北京市機(jī)動(dòng)車的NOx逐年排放量,計(jì)算公式如下:

    式中:Pj為計(jì)算年份機(jī)動(dòng)車保有量,104輛;j為車型分類;Mj為平均行駛里程,104km;Efj為不同車型的平均排放因子,g/km.

    計(jì)算的2006~2014年北京市機(jī)動(dòng)車的NOx年排放量在7.65×104~15.3×104t,受機(jī)動(dòng)車數(shù)量及油品標(biāo)準(zhǔn)等因素的影響,近兩年機(jī)動(dòng)車排放量在7.5×104t左右,近8年機(jī)動(dòng)車的NOx年排放量與NO2濃度的相關(guān)系數(shù)R為0.77,可見迅速增長的機(jī)動(dòng)車保有量對北京市地區(qū)上空的NO2濃度的變化有顯著的影響.

    3 結(jié)論

    3.1 2013年NO2平均濃度為56μg/m3,2014年北京市NO2年均濃度為56.7μg/m3,2014年與2013年相比上升約1.3%.

    3.2 年均及重污染日NO2月均濃度呈波浪型分布,日變化呈現(xiàn)雙峰型分布;空間分布上北部及西部山區(qū)NO2濃度明顯低于中心城區(qū)及南部地區(qū).3.3 NO2濃度與PM2.5、CO、NO呈正相關(guān),與O3、OX負(fù)相關(guān);全年NO2光解速率峰值平均在0.105/min左右,重污染日光解速率峰值平均在0.026/min左右;全年及重污染期間,氮氧化速率分別為0.142±0.061、0.190±0.036;高濃度NO2既有利于O3生成又對重污染的形成起到了促進(jìn)作用,重污染日特定條件下北京市NO2的兩種轉(zhuǎn)化機(jī)制以轉(zhuǎn)化為NO3-過程為主導(dǎo).

    3.4 2000~2014年北京市機(jī)動(dòng)車的保有量與NO2濃度的相關(guān)系數(shù)R為-0.84,機(jī)動(dòng)車NOx排放量對北京市地區(qū)上空的NO2濃度的變化有顯著的影響.

    參考文獻(xiàn):

    [1] Sachin D G, Vander A, Beig G.Satellite derived trends in N O2over the major global hotspot regions during the past decade and their inter-comparison [J].Environ.Pollut., 2009,157(3):1873-1878.

    [2] Sun Y L, Zhuang G S, Tang A H, et al.Chemical characteristics of PM2.5and PM10in haze-fog episodes in Beijing [J].Environ Sci Technol, 2006,40(10):3148-3155.

    [3] 李令軍,王 英,李金香,等.2000~2010北京大氣重污染研究 [J].中國環(huán)境科學(xué), 2012,32(1):23-30.

    [4] Velders G J M, Granier C, Portmann R W, et al.Global tropospheric NO2column distributions:comparing threedimensional model calculations with GOME measurements [J].Journal of Geophysical Research, 2001,106:12643-12660.

    [5] Richter A, Burrows J P, Granier C, et al.Increase in tropospheric nitrogen dioxide over China observed from space [J].Nature, 2005,473:129-132.

    [6] Kang C H, Kim W H, Ko H J, et al.Asian dust effects on total suspended particulate (TSP) compositions at Gosan in Jeju Island, Korea [J].Atmospheric Research, 2009,94(2):345-355.

    [7] Zhang Q, Geng G N, Wang S W, et al.Satellite remote sensing of changes in N Oxemissions over China:1996-2010 [J].Chinese Science Bull, 2012,57(16):1446-1453.

    [8] 王郭臣,王 玨,信玉潔,等.天津PM10和NO2輸送路徑及潛在源區(qū)研究 [J].中國環(huán)境科學(xué), 2014,34(12):3009-3016.

    [9] 田賀忠,王 艷,趙 丹,等.中國太行山東麓NOx重污染成因分析 [J].科學(xué)通報(bào), 2011,56(18):1464-1469.

    [10] 李令軍,王 英.基于衛(wèi)星遙感與地面監(jiān)測分析北京大氣NO2污染特征 [J].環(huán)境科學(xué)學(xué)報(bào), 2011,31(12):2762-2768.

    [11] 尉 鵬,任陣海,蘇福慶,等.中國NO2的季節(jié)分布及成因分析[J].環(huán)境科學(xué)研究, 2011,24(2):155-161.

    [12] 江文華,馬建中,顏 鵬,等.利用GOME衛(wèi)星資料分析北京大氣NO2污染變化 [J].應(yīng)用氣象學(xué), 2006,17(1):67-72.

    [13] 王躍啟,江 洪,張秀英,等.基于OMI衛(wèi)星遙感數(shù)據(jù)的中國對流層NO2時(shí)空分布 [J].環(huán)境科學(xué)研究, 2009,22(8):992-997.

    [14] Schaap M, Kranenburg R, Curier L, et al.Assessing the sensitivity of the OMI-NO2product to emission changes across Europe [J].Remote Sensing, 2013,5:4187-4208.

    [15] Cleveland W S, Graedel T E, Kleiner B,et al.Sunday and workday variations in photochemical air pollutants in New Jersey and New York [J].Science, 1974,186(4168):1037-1038.

    [16] 中國環(huán)境年鑒編輯委員會(huì).中國環(huán)境年鑒 [M].北京:中國環(huán)境年鑒社, 2000-2014.

    [17] 謝淑艷,王曉彥,吳迓名,等.環(huán)境空氣中PM2.5自動(dòng)監(jiān)測方法比較及應(yīng)用 [J].中國環(huán)境監(jiān)測, 2013,29(2):150-155.

    [18] 環(huán)保部.HJ/T 193-2005 環(huán)境空氣質(zhì)量自動(dòng)監(jiān)測技術(shù)規(guī)范 [S].

    [19] 北京市環(huán)保局.北京市大氣污染現(xiàn)狀及來源解析研究 [R].北京:北京環(huán)保局, 2013.

    [20] 環(huán)保部.HJ 633-2012 環(huán)境空氣質(zhì)量指數(shù)(AQI)技術(shù)規(guī)定(試行) [S].

    [21] 王從梅,楊永勝,李永占,等.2013年1月河北省中南部嚴(yán)重污染的氣象條件及成因分析 [J].環(huán)境科學(xué)研究, 2013,26(7):695-700.

    [22] 程念亮,李云婷,張大偉,等.2014年10月北京市4次典型空氣重污染過程成因分析 [J].環(huán)境科學(xué)研究, 2015,28(2):1-10.

    [23] 王 英,李令軍,劉 陽.京津冀與長三角區(qū)域大氣NO2污染特征 [J].環(huán)境科學(xué), 2012,33(11):3685-3692.

    [24] 任陣海,萬本太,虞 統(tǒng),等.不同尺度大氣系統(tǒng)對污染邊界層的影響及其水平流場輸送 [J].環(huán)境科學(xué)研究, 2004,17(1):2004,7-13.

    [25] Sillman S.The sensitivity of ozone to nitrogen oxides and hydrocarbons in regional ozone episodes [J].Journal of Geophysical Research, 1990,95(D2):1837-1851.

    [26] 安俊琳,王躍思,孫 揚(yáng).氣象因素對北京臭氧的影響 [J].生態(tài)環(huán)境學(xué)報(bào), 2009,18(3):944-951.

    [27] Davis, John C.Statistics and Data Analysis in Geology (3rd Edition) [M].New York: John Wiley & Sons, Inc, 2002:57-61.

    [28] Lal S, Naja M, Subbaraya B H.Seasonal variations in surface ozone and its precursors over and urban site in India [J].Atmospheric Environment, 2000,34:2713-2724.

    [29] He K B, Yang F M, Ma Y L, et al.The characteristics of PM2.5in Beijing, China [J].Atmospheric Environment, 2001,35(29):4959-4970.

    [30] Elkus B, Wilson K R.Photochemical air pollution: weekendweekday differences [J].Atmospheric Environment, 1977,11(6): 509-515.

    [31] 黃志新.上海市郊春季臭氧及其前體物觀測研究 [J].環(huán)境科學(xué)與技術(shù), 2011,34(5):87-99.

    [32] Tang G, Li X, Wang Y, et al.Surface ozone trend details and interpretations in Beijing, 2001-2006 [J].Atmospheric Chemistry and Physics, 2009,9:8813-8823.

    [33] Jacob D J, Logan J A, Murti P P.Effect of rising Asian emissions on surface ozone in the United States [J].Geophysical Research Letters, 1999,26:2175-2178.

    [34] Clapp L J, Jenkin M E.Analysis of the relationship between ambient levels of O3, NO2and NO as a function of NOxin UK [J].Atmospheric Environment, 1995,29:923-946.

    [35] Olcese L E, Palancar G G, Toselli B M.An inexpensive method to estimate CO and NOxemissions from mobile sources [J].Atmospheric Environment, 2001,35:6213-6218.

    [36] 王雪松,李金龍.人為源排放VOC對北京地區(qū)臭氧生成的貢獻(xiàn)[J].中國環(huán)境科學(xué), 2002,22(6):501-505.

    [37] Han S, Bian H, Feng Y, et al.Analysis of the Relationship between O3, NO and NO2in Tianjin, China [J].Aerosol Air Qual.Res, 2011,11:128-139.

    [38] Wang Z S, Zhang D W, Chen T, et al.Spatiotemporal characteristics of NO2in Beijing in 2013 [J].Acta Scientiae Circumstantiae, 2015,35(5):1529-1536.

    [39] Yao X H, Chan C K, Fang M, et al.The water-soluble ionic composition of PM2.5in Shanghai and Beijing, China [J].Atmospheric Environment, 2002,36(26):4223-4234.

    [40] Hu M, He L Y, Zhang Y H, et al.Seasonal variation of ionic species in fine particles at Qingdao, China [J].Atmospheric Environment, 2002,36:5853-5859.

    [41] Wang Y, Zhuang G S, Tang A H, et al.The ion chemistry and the source of PM2.5aerosol in Beijing [J].Atmospheric Environment, 2005,39(21):3771-3784.

    [42] 鄧?yán)?李 紅,柴發(fā)合,等.北京東北部城區(qū)大氣細(xì)粒子與相關(guān)氣體污染特征研究 [J].中國環(huán)境科學(xué), 2011,31(7):1064-1070.

    [43] 金 鑫,程萌田,溫天雪,等.北京冬季一次重污染過程PM2.5中水溶性無機(jī)鹽的變化特征 [J].環(huán)境化學(xué), 2012,31(6):783-790.

    [44] 周 敏,陳長虹,喬利平,等.2013年1月中國中東部大氣重污染期間上海顆粒物的污染特征 [J].環(huán)境科學(xué)學(xué)報(bào), 2013,33(11): 3118-3126.

    [45] 廖碧婷,吳 兌,常 越,等.廣州地區(qū)NO3-、SO42-、NH4+與相關(guān)氣體污染特征研究 [J].環(huán)境科學(xué)學(xué)報(bào), 2014,34(6):1551-1559.

    [46] 李麗珍,沈振興,杜 娜,等.霾和正常天氣下西安大氣顆粒物中水溶性離子特征 [J].中國科學(xué)院研究生院學(xué)報(bào), 2007,24(5):674-670.

    [47] 高曉梅.我國典型地區(qū)大氣PM2.5水溶性離子的理化特征及來源解析 [D].濟(jì)南:山東大學(xué), 2012.

    [48] 張興贏,張 鵬,張 艷,等.近10a中國對流層NO2的變化趨勢、時(shí)空分布特征及其來源解析 [J].中國科學(xué)D輯:地球科學(xué), 2007,37(10):1409-1416.

    [49] 陳煥盛,王自發(fā),吳其重,等.亞運(yùn)時(shí)段廣州大氣污染物來源數(shù)值模擬研究 [J].環(huán)境科學(xué)學(xué)報(bào), 2010,30(11):2145 -2153.

    [50] 周 丹.東莞市大氣環(huán)境污染來源分析與模擬研究 [D].北京:北京化工大學(xué), 2008.

    [51] Wang Y, McElroy M B, Boersma K F, et al.Traffic restrictions associated with the Sino-African summit: reductions of NOxdetected from space [J].Geophysical Reseach Letter, 2007,34, L08814,doi:10.1029/2007GL029326.

    [52] 傅立新,郝吉明,何東全,等.北京市機(jī)動(dòng)車污染物排放特征 [J].環(huán)境科學(xué)研究, 2000,21(3):68-74.

    Spatial and temporal distribution of NO2during 2013~2014 in Beijing.


    CHENG Nian-liang1,2,3, LI Yun-ting1, CHEN Tian4, ZHANG Da-wei1*, DONG Xin1, WANG Xin1, HUAN Ning1, LIU Bao-xian1, YAN He1, MENG Fan2,3**(1.Beijing Key Laboratory of Airborne Particulate Matter Monitoring Technology, Beijing Municipal Environmental Monitoring Center, Beijing 100048, China;2.College of Water Sciences, Beijing Normal University, Beijing 100875, China;3.Chinese Research Academy of Environmental Sciences, Beijing 100012, China;4.Beijing Environmental Protection Bureau, Beijing 100048, China).China Environmental Science, 2016,36(1):18~26

    Abstract:The Spatial and temporal distribution characteristics of NO2during 2013~2014 were analyzed based on the NO2monitoring data released by Beijing Municipal Environmental Monitoring Center.The analytical results showed that the annual average concentration of NO2in Beijing was 56.7μg/m3in 2014and increased 1.3% compared to that of NO2(56μg/m3) in 2013.The monthly average concentration of NO2changed in a wave-shape curve and the diurnal variation of NO2showed bimodal distribution.The annual average concentration of NO2was lower in the northern and western regions while higher in the six city area and southern area; NO2and PM2.5, CO, NO showed a positive correlation while NO2and OX, O3had no significant correlation relationship.Calculated photolysis rate of NO2in the whole year and heavy pollution days was 0.105/min, 0.026/min, respectively.The values of NOR was 0.142±0.061, 0.190±0.036 in the whole year and heavy pollution days, respectively.High concentration of NO2was conducive to the formation of O3and also played a promoting role of heavy air pollution formation.Correlation coefficient between vehicle number and the annuanl concentration of NO2during 2000~2014 in Beijing was -0.84 and the NOxemissions of vehicle had a significant effect on the change of the annual concentration of NO2in the area of Beijing city.

    Key words:NO2;spatial and temporal distribution;Beijing;NOR;PM2.5;photolysis rate

    中圖分類號(hào):X513

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1000-6923(2016)01-0018-09

    收稿日期:2015-06-09

    基金項(xiàng)目:環(huán)保公益性行業(yè)科研專項(xiàng)(201409005);國家科技支撐計(jì)劃課題(2014BAC23B03);北京市科技計(jì)劃課題(Z131100005613046);中心課題(2015-01-02);北京市VOCs基準(zhǔn)監(jiān)測體系評估及環(huán)境水平狀況研究課題

    作者簡介:程念亮(1987-),男,山東泰安人,工程師,博士研究生,主要從事空氣質(zhì)量環(huán)境預(yù)警預(yù)報(bào)及大氣環(huán)境規(guī)劃模擬研究.發(fā)表論文20篇.

    猜你喜歡
    大氣污染
    大氣的呵護(hù)
    軍事文摘(2023年10期)2023-06-09 09:15:06
    太赫茲大氣臨邊探測儀遙感中高層大氣風(fēng)仿真
    什么是污染?
    什么是污染?
    堅(jiān)決打好污染防治攻堅(jiān)戰(zhàn)
    堅(jiān)決打好污染防治攻堅(jiān)戰(zhàn)
    大氣古樸揮灑自如
    大氣、水之后,土十條來了
    對抗塵污染,遠(yuǎn)離“霾”伏
    都市麗人(2015年5期)2015-03-20 13:33:49
    世界知識(shí)畫報(bào)·藝術(shù)視界(2010年9期)2010-12-31 00:00:00
    国产国拍精品亚洲av在线观看| 女生性感内裤真人,穿戴方法视频| 国产私拍福利视频在线观看| 成人午夜高清在线视频| 成年女人毛片免费观看观看9| 欧美最黄视频在线播放免费| 婷婷精品国产亚洲av| 成人三级黄色视频| 嫩草影院精品99| 亚洲av免费高清在线观看| 成年女人毛片免费观看观看9| 久久久久免费精品人妻一区二区| 亚洲 国产 在线| 国产欧美日韩一区二区三| 国产高清有码在线观看视频| 亚洲人成网站在线播| 特级一级黄色大片| 18禁黄网站禁片免费观看直播| 舔av片在线| 精品人妻熟女av久视频| 亚洲中文日韩欧美视频| 国产精品久久久久久亚洲av鲁大| 久久久久久久午夜电影| x7x7x7水蜜桃| 国产精品不卡视频一区二区 | 91久久精品国产一区二区成人| 简卡轻食公司| 18+在线观看网站| 久久香蕉精品热| 一a级毛片在线观看| 久久久久久大精品| xxxwww97欧美| 怎么达到女性高潮| 国产精品不卡视频一区二区 | 国产伦精品一区二区三区视频9| 精品久久久久久久久亚洲 | 少妇人妻一区二区三区视频| 久久热精品热| 又黄又爽又免费观看的视频| 日韩大尺度精品在线看网址| 1000部很黄的大片| 亚洲精品在线美女| 中文字幕av成人在线电影| 亚洲欧美日韩无卡精品| 免费av毛片视频| 国产黄a三级三级三级人| 国产黄a三级三级三级人| 超碰av人人做人人爽久久| 天天躁日日操中文字幕| 久9热在线精品视频| 老司机午夜十八禁免费视频| 久久久久久久久中文| 国产高清视频在线播放一区| 一个人免费在线观看的高清视频| 日本撒尿小便嘘嘘汇集6| 中文字幕av在线有码专区| 精品人妻偷拍中文字幕| 亚洲av日韩精品久久久久久密| 成人性生交大片免费视频hd| 亚洲,欧美,日韩| 亚洲美女黄片视频| 欧美一区二区亚洲| 一卡2卡三卡四卡精品乱码亚洲| 国产午夜精品论理片| 神马国产精品三级电影在线观看| 国产一区二区在线av高清观看| 免费观看的影片在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产精品一区二区三区四区免费观看 | 桃红色精品国产亚洲av| 中文在线观看免费www的网站| 日本成人三级电影网站| 欧美午夜高清在线| 免费观看的影片在线观看| 久久国产精品人妻蜜桃| 国产一级毛片七仙女欲春2| 亚州av有码| 嫁个100分男人电影在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 一个人免费在线观看电影| 久久99热6这里只有精品| 韩国av一区二区三区四区| 亚洲精品一区av在线观看| 成人特级av手机在线观看| 国产淫片久久久久久久久 | 亚洲av电影不卡..在线观看| 国产中年淑女户外野战色| 精品一区二区三区视频在线| 久久久久九九精品影院| 色播亚洲综合网| 免费观看人在逋| 一个人观看的视频www高清免费观看| 国产黄色小视频在线观看| 国产一区二区三区在线臀色熟女| 俄罗斯特黄特色一大片| 99久久无色码亚洲精品果冻| 91在线精品国自产拍蜜月| 亚洲国产欧洲综合997久久,| 色吧在线观看| 欧美日韩中文字幕国产精品一区二区三区| 日本免费a在线| 一本久久中文字幕| 亚洲精华国产精华精| 免费搜索国产男女视频| 欧美黄色淫秽网站| 国产成+人综合+亚洲专区| 97碰自拍视频| 精品一区二区免费观看| 男女床上黄色一级片免费看| 午夜福利免费观看在线| 亚洲国产欧美人成| 成年女人永久免费观看视频| 超碰av人人做人人爽久久| 日韩av在线大香蕉| 日韩有码中文字幕| 国产精品亚洲av一区麻豆| 国内久久婷婷六月综合欲色啪| 国产激情偷乱视频一区二区| 在线十欧美十亚洲十日本专区| 麻豆久久精品国产亚洲av| 亚洲美女视频黄频| a级毛片a级免费在线| 18禁裸乳无遮挡免费网站照片| 中出人妻视频一区二区| 亚洲 欧美 日韩 在线 免费| 美女xxoo啪啪120秒动态图 | 国产高清视频在线观看网站| 久久99热6这里只有精品| 国产精品一区二区三区四区免费观看 | 国产伦在线观看视频一区| 久久国产乱子伦精品免费另类| 女人十人毛片免费观看3o分钟| www.熟女人妻精品国产| 欧美又色又爽又黄视频| 一边摸一边抽搐一进一小说| 欧美日本视频| 久久婷婷人人爽人人干人人爱| 18禁在线播放成人免费| 欧美激情久久久久久爽电影| 人妻丰满熟妇av一区二区三区| 欧美国产日韩亚洲一区| 精品久久久久久久人妻蜜臀av| 欧美日本视频| 欧美潮喷喷水| 搞女人的毛片| 国产精品综合久久久久久久免费| 宅男免费午夜| 国产高清有码在线观看视频| 中文字幕精品亚洲无线码一区| 国产 一区 欧美 日韩| 搞女人的毛片| 午夜免费成人在线视频| 国产成人福利小说| 中亚洲国语对白在线视频| 亚洲久久久久久中文字幕| 国产伦一二天堂av在线观看| 国产一区二区三区视频了| 美女免费视频网站| 少妇人妻一区二区三区视频| 神马国产精品三级电影在线观看| 国内揄拍国产精品人妻在线| 免费高清视频大片| 男女视频在线观看网站免费| 首页视频小说图片口味搜索| 日韩欧美在线乱码| 熟女电影av网| 午夜亚洲福利在线播放| 国产高清三级在线| 18禁黄网站禁片免费观看直播| 一个人免费在线观看电影| 一区二区三区四区激情视频 | 国产免费男女视频| 悠悠久久av| 久久中文看片网| 久久九九热精品免费| 中文在线观看免费www的网站| 一个人看视频在线观看www免费| 亚洲成人免费电影在线观看| 午夜福利欧美成人| 少妇的逼好多水| 久久午夜福利片| 成人性生交大片免费视频hd| 国产精品嫩草影院av在线观看 | 国产精品亚洲av一区麻豆| 性欧美人与动物交配| 国内久久婷婷六月综合欲色啪| 国产探花极品一区二区| 日韩欧美国产一区二区入口| 亚洲熟妇熟女久久| 久久久久九九精品影院| 熟妇人妻久久中文字幕3abv| 久久中文看片网| 中文字幕精品亚洲无线码一区| 成年版毛片免费区| 亚洲国产欧美人成| 99热这里只有精品一区| 97碰自拍视频| 午夜精品在线福利| 欧美日本视频| 亚洲精品粉嫩美女一区| 色av中文字幕| 欧美一区二区精品小视频在线| 亚洲欧美日韩卡通动漫| 夜夜躁狠狠躁天天躁| 久久国产乱子免费精品| 毛片女人毛片| 国产精品1区2区在线观看.| 欧美一区二区国产精品久久精品| 日韩欧美在线二视频| 一进一出抽搐gif免费好疼| 搡老妇女老女人老熟妇| 久久久久久久久中文| 国内揄拍国产精品人妻在线| 毛片女人毛片| 熟妇人妻久久中文字幕3abv| 有码 亚洲区| 中文字幕免费在线视频6| 欧美一区二区精品小视频在线| 狠狠狠狠99中文字幕| 久久久久久久亚洲中文字幕 | 日韩欧美国产在线观看| 人人妻人人澡欧美一区二区| 99久久无色码亚洲精品果冻| 国产精品人妻久久久久久| 国产伦精品一区二区三区视频9| 欧美xxxx性猛交bbbb| 又黄又爽又免费观看的视频| 亚洲av日韩精品久久久久久密| 亚洲成a人片在线一区二区| 精品福利观看| 国产欧美日韩精品亚洲av| 精品一区二区三区视频在线| 国产 一区 欧美 日韩| av专区在线播放| 国内精品久久久久精免费| 白带黄色成豆腐渣| 老熟妇仑乱视频hdxx| 丰满的人妻完整版| 午夜老司机福利剧场| 人妻丰满熟妇av一区二区三区| 人妻夜夜爽99麻豆av| 亚洲熟妇熟女久久| 桃色一区二区三区在线观看| 欧美三级亚洲精品| 老鸭窝网址在线观看| 在线国产一区二区在线| 精品久久久久久久久av| 久久国产精品人妻蜜桃| 九九久久精品国产亚洲av麻豆| 亚洲成a人片在线一区二区| 午夜福利在线观看免费完整高清在 | 国产精品伦人一区二区| 欧美+亚洲+日韩+国产| 欧美xxxx黑人xx丫x性爽| 十八禁国产超污无遮挡网站| 亚洲成人中文字幕在线播放| 桃色一区二区三区在线观看| 亚洲精品亚洲一区二区| 午夜福利视频1000在线观看| 国产乱人伦免费视频| 舔av片在线| 欧美一级a爱片免费观看看| 国产三级在线视频| 51国产日韩欧美| 精品不卡国产一区二区三区| 两个人视频免费观看高清| 美女被艹到高潮喷水动态| 久久久久久久久大av| 在线看三级毛片| 亚洲人成网站在线播放欧美日韩| 97人妻精品一区二区三区麻豆| 国产国拍精品亚洲av在线观看| a在线观看视频网站| 亚洲av电影不卡..在线观看| 欧美日本视频| 精品久久久久久久久亚洲 | 欧美成人一区二区免费高清观看| 直男gayav资源| 欧美日韩黄片免| 亚洲18禁久久av| 可以在线观看的亚洲视频| 少妇的逼水好多| 亚洲成人精品中文字幕电影| 老司机午夜十八禁免费视频| 91狼人影院| 亚洲av成人不卡在线观看播放网| 亚洲av一区综合| 精品久久久久久久久久久久久| 日本五十路高清| 日韩亚洲欧美综合| 日本黄色视频三级网站网址| 一夜夜www| 少妇的逼水好多| 亚洲熟妇中文字幕五十中出| 最近最新中文字幕大全电影3| 久久精品夜夜夜夜夜久久蜜豆| 在线看三级毛片| 老熟妇仑乱视频hdxx| 禁无遮挡网站| 亚洲电影在线观看av| 丝袜美腿在线中文| 欧美精品啪啪一区二区三区| 欧美日韩乱码在线| 看黄色毛片网站| 宅男免费午夜| 精品日产1卡2卡| 人人妻,人人澡人人爽秒播| 午夜精品在线福利| 精品不卡国产一区二区三区| 性色avwww在线观看| 又黄又爽又免费观看的视频| 99国产精品一区二区蜜桃av| 欧美+日韩+精品| 嫩草影院精品99| 免费一级毛片在线播放高清视频| 亚洲综合色惰| 国产av一区在线观看免费| 中文字幕人成人乱码亚洲影| 免费看日本二区| 97碰自拍视频| 老司机午夜十八禁免费视频| 最好的美女福利视频网| 欧美3d第一页| 亚洲美女黄片视频| 免费看日本二区| 人人妻,人人澡人人爽秒播| 简卡轻食公司| 好男人电影高清在线观看| 麻豆成人午夜福利视频| 1000部很黄的大片| 亚洲国产色片| 一进一出抽搐gif免费好疼| 亚洲成av人片免费观看| 在线观看免费视频日本深夜| 在线观看午夜福利视频| 小说图片视频综合网站| 桃红色精品国产亚洲av| 999久久久精品免费观看国产| 一个人观看的视频www高清免费观看| 特级一级黄色大片| 亚洲av熟女| 天天一区二区日本电影三级| 免费观看人在逋| 色5月婷婷丁香| 听说在线观看完整版免费高清| 久久久久久久久中文| 亚洲av日韩精品久久久久久密| 国产欧美日韩精品一区二区| 亚洲国产欧美人成| 9191精品国产免费久久| 制服丝袜大香蕉在线| 国产精品不卡视频一区二区 | 国产乱人伦免费视频| 日日干狠狠操夜夜爽| 日本熟妇午夜| 岛国在线免费视频观看| 精品福利观看| 最后的刺客免费高清国语| 又紧又爽又黄一区二区| 一边摸一边抽搐一进一小说| 精品国产亚洲在线| 精品欧美国产一区二区三| 日本黄色视频三级网站网址| 精品一区二区三区视频在线| 一区二区三区高清视频在线| 精品一区二区三区视频在线| 国内少妇人妻偷人精品xxx网站| 99在线视频只有这里精品首页| 九九在线视频观看精品| 午夜福利视频1000在线观看| 99久久久亚洲精品蜜臀av| 一本久久中文字幕| 每晚都被弄得嗷嗷叫到高潮| 久久久久国内视频| 国产一区二区三区视频了| 亚洲欧美日韩高清专用| 一进一出抽搐gif免费好疼| 成人av在线播放网站| 日本黄大片高清| 麻豆国产av国片精品| 日韩欧美 国产精品| 麻豆国产av国片精品| 精品一区二区三区人妻视频| 亚洲人成电影免费在线| 日本五十路高清| 久99久视频精品免费| 一进一出抽搐gif免费好疼| 亚洲成人久久爱视频| 亚洲av成人不卡在线观看播放网| 老司机午夜十八禁免费视频| 精品久久久久久,| 女人被狂操c到高潮| 一级a爱片免费观看的视频| 国产v大片淫在线免费观看| 国产不卡一卡二| 校园春色视频在线观看| 免费av毛片视频| 亚洲天堂国产精品一区在线| 长腿黑丝高跟| 国产精品不卡视频一区二区 | 麻豆成人av在线观看| 亚洲av熟女| 亚洲精品一卡2卡三卡4卡5卡| 两个人视频免费观看高清| 国产黄片美女视频| 久久九九热精品免费| 国产高清视频在线播放一区| 亚洲成a人片在线一区二区| 人人妻人人澡欧美一区二区| 校园春色视频在线观看| 丰满人妻一区二区三区视频av| 精品久久久久久久久久免费视频| 国产成人福利小说| 91狼人影院| 国产高清有码在线观看视频| 国产老妇女一区| 在现免费观看毛片| 国内少妇人妻偷人精品xxx网站| 午夜激情福利司机影院| 亚洲成人精品中文字幕电影| 亚洲国产精品久久男人天堂| 国产高清视频在线观看网站| 免费av观看视频| 男女之事视频高清在线观看| 一级黄色大片毛片| 国产真实乱freesex| 精品人妻偷拍中文字幕| 欧美中文日本在线观看视频| av国产免费在线观看| 免费黄网站久久成人精品 | 18+在线观看网站| a级毛片a级免费在线| 熟妇人妻久久中文字幕3abv| 18美女黄网站色大片免费观看| 亚洲专区国产一区二区| 国产熟女xx| 免费大片18禁| 成年版毛片免费区| 精品欧美国产一区二区三| 在线观看av片永久免费下载| 欧美黄色淫秽网站| 日韩欧美一区二区三区在线观看| 午夜精品一区二区三区免费看| 色哟哟·www| 在线国产一区二区在线| 女人十人毛片免费观看3o分钟| 在线观看免费视频日本深夜| 精品熟女少妇八av免费久了| 欧美日本亚洲视频在线播放| 日本三级黄在线观看| 免费电影在线观看免费观看| 小蜜桃在线观看免费完整版高清| 国产高清有码在线观看视频| 99视频精品全部免费 在线| 99精品久久久久人妻精品| 最好的美女福利视频网| 国内久久婷婷六月综合欲色啪| 国产精品久久久久久精品电影| 国产高清三级在线| 我要看日韩黄色一级片| 亚洲经典国产精华液单 | 午夜福利18| 一卡2卡三卡四卡精品乱码亚洲| 免费电影在线观看免费观看| 日韩欧美在线乱码| 久久亚洲真实| 每晚都被弄得嗷嗷叫到高潮| 免费看日本二区| 亚洲激情在线av| 看免费av毛片| 精品免费久久久久久久清纯| 亚洲中文字幕日韩| 老司机午夜十八禁免费视频| 婷婷六月久久综合丁香| 国产精品野战在线观看| 搡老妇女老女人老熟妇| 欧美区成人在线视频| 午夜精品久久久久久毛片777| 久久人人精品亚洲av| 91麻豆av在线| 一进一出好大好爽视频| 怎么达到女性高潮| 少妇裸体淫交视频免费看高清| 757午夜福利合集在线观看| 成人无遮挡网站| 精品人妻偷拍中文字幕| 成熟少妇高潮喷水视频| 亚洲精品色激情综合| 国产亚洲精品av在线| 性插视频无遮挡在线免费观看| 999久久久精品免费观看国产| 免费在线观看成人毛片| 99久久精品热视频| 波多野结衣巨乳人妻| 精品免费久久久久久久清纯| 欧美区成人在线视频| 三级毛片av免费| www日本黄色视频网| 观看免费一级毛片| 久久久久久久午夜电影| 网址你懂的国产日韩在线| 麻豆成人午夜福利视频| 亚洲精品成人久久久久久| 老女人水多毛片| a级毛片a级免费在线| 国内精品一区二区在线观看| 最近在线观看免费完整版| 少妇人妻一区二区三区视频| 精品一区二区三区av网在线观看| 欧美日韩乱码在线| 成人特级黄色片久久久久久久| 日韩 亚洲 欧美在线| 五月伊人婷婷丁香| 久久久久久久亚洲中文字幕 | 在线观看66精品国产| 简卡轻食公司| 午夜a级毛片| 日日夜夜操网爽| 女人十人毛片免费观看3o分钟| 身体一侧抽搐| 特大巨黑吊av在线直播| 少妇的逼好多水| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久九九精品影院| 色在线成人网| 成人美女网站在线观看视频| 亚洲av第一区精品v没综合| 亚洲中文字幕一区二区三区有码在线看| 免费在线观看成人毛片| 国产私拍福利视频在线观看| 91狼人影院| 日韩高清综合在线| 午夜影院日韩av| 九色国产91popny在线| 国产一区二区亚洲精品在线观看| 国产黄片美女视频| 国产中年淑女户外野战色| 精品一区二区免费观看| 亚洲国产日韩欧美精品在线观看| 男人的好看免费观看在线视频| 大型黄色视频在线免费观看| 内地一区二区视频在线| 国产精品爽爽va在线观看网站| 久久国产精品影院| 成年免费大片在线观看| 久久中文看片网| 欧美激情国产日韩精品一区| 久久久久九九精品影院| 亚洲美女黄片视频| 真实男女啪啪啪动态图| 黄色视频,在线免费观看| 99视频精品全部免费 在线| 精品欧美国产一区二区三| 日韩人妻高清精品专区| 一级作爱视频免费观看| 免费看a级黄色片| 一卡2卡三卡四卡精品乱码亚洲| 一进一出好大好爽视频| 欧美最黄视频在线播放免费| 两个人视频免费观看高清| 欧美一区二区精品小视频在线| 日韩欧美在线二视频| 国产一区二区三区在线臀色熟女| 在线观看66精品国产| 99久久精品国产亚洲精品| 99热只有精品国产| 在线十欧美十亚洲十日本专区| 成人性生交大片免费视频hd| 久久精品国产亚洲av香蕉五月| netflix在线观看网站| 性色avwww在线观看| 亚洲美女黄片视频| 国产成人a区在线观看| 少妇裸体淫交视频免费看高清| 免费在线观看影片大全网站| 人妻久久中文字幕网| 有码 亚洲区| 久久国产精品影院| 女人被狂操c到高潮| 亚洲熟妇熟女久久| av天堂中文字幕网| 又爽又黄无遮挡网站| 国产精品一区二区三区四区久久| 伦理电影大哥的女人| 在线看三级毛片| 国产成年人精品一区二区| 又黄又爽又免费观看的视频| 少妇人妻精品综合一区二区 | 久久伊人香网站| 国产精品精品国产色婷婷| 国内精品美女久久久久久| 亚洲久久久久久中文字幕| 乱人视频在线观看| 我的老师免费观看完整版| 别揉我奶头 嗯啊视频| 亚洲性夜色夜夜综合| 国产精品一区二区性色av| 精品久久久久久成人av| 亚洲无线在线观看| 久久国产精品影院| 免费看日本二区| 欧美一区二区国产精品久久精品| 全区人妻精品视频| 欧美另类亚洲清纯唯美| 我要看日韩黄色一级片| 人妻夜夜爽99麻豆av| 国产av一区在线观看免费| 国产黄色小视频在线观看| 好男人电影高清在线观看| 又黄又爽又刺激的免费视频.| 一区二区三区四区激情视频 | 精品人妻一区二区三区麻豆 | 少妇丰满av| 色综合婷婷激情| 国产欧美日韩一区二区精品|