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

    2013~2015年北京市PM2.5、反應(yīng)性氣體和氣溶膠粒徑的特性分析

    2017-11-07 04:47:51張恒德呂夢(mèng)瑤何佳寶
    中國(guó)環(huán)境科學(xué) 2017年10期
    關(guān)鍵詞:大氣風(fēng)速模態(tài)

    江 琪,王 飛,張恒德,王 飛,呂夢(mèng)瑤,何佳寶

    ?

    2013~2015年北京市PM2.5、反應(yīng)性氣體和氣溶膠粒徑的特性分析

    江 琪1,王 飛2*,張恒德1,王 飛3,呂夢(mèng)瑤1,何佳寶4

    (1.國(guó)家氣象中心,北京100081;2.中國(guó)氣象科學(xué)研究院,北京100081;3.北京市人工影響天氣辦公室,北京100089;4.寧波市環(huán)境監(jiān)測(cè)中心,浙江寧波 315012)

    對(duì)2013~2015北京市PM2.5、氣體組分及氣溶膠譜分布結(jié)合氣象條件進(jìn)行系統(tǒng)的分析,結(jié)果表明:2013~2015年 PM2.5質(zhì)量濃度顯著降低,空氣質(zhì)量呈明顯改善;交通源對(duì)大氣污染的貢獻(xiàn)不能忽略,且呈現(xiàn)出逐漸加重的趨勢(shì);相對(duì)濕度和風(fēng)速對(duì)PM2.5濃度影響的響應(yīng)相對(duì)較快,滯后1h左右對(duì)PM2.5濃度的影響最大,而PM2.5的變化對(duì)溫度的響應(yīng)相對(duì)較慢,一般在改變溫度的7~8h后對(duì)污染物濃度的影響最大.污染日較清潔日氣溶膠的粒徑分布均有不同程度地向大粒子一側(cè)拓寬的現(xiàn)象,清潔日時(shí),峰區(qū)直徑集中在3~5nm,污染天時(shí),峰區(qū)直徑增大至100~200nm;隨降雨量級(jí)逐漸增大,氣溶膠峰值濃度逐漸向積聚模態(tài)方向發(fā)展;2015年2月18至2015年2月21日污染事件中,煙花燃放時(shí)段煙花對(duì)氣溶膠數(shù)濃度的貢獻(xiàn)為26.1%,且伴隨新粒子生成.煙花燃放使積聚模態(tài)粒子和愛(ài)根核模態(tài)粒子出現(xiàn)顯著的峰值,核模態(tài)粒子濃度變化不大.

    PM2.5;反應(yīng)性氣體;氣象要素;氣溶膠粒徑譜;煙花燃放

    北京位于華北平原西北邊緣,三面環(huán)山,地形呈“簸箕型”,不利于污染物擴(kuò)散[1],多種污染物高濃度并存,并發(fā)生復(fù)雜的相互作用,使得大氣呈現(xiàn)出細(xì)顆粒物濃度升高、氧化性增強(qiáng)等特點(diǎn)[2].在污染本質(zhì)上則表現(xiàn)為污染物源和匯之間的相互交錯(cuò),污染轉(zhuǎn)化過(guò)程的耦合作用和環(huán)境影響的協(xié)同效應(yīng)[3-4].在多種污染物中,高濃度的大氣細(xì)粒子(PM2.5)是引起城市霾天氣發(fā)生的主要原因[5-6].工業(yè)、交通、烹飪等一次源排放以及多種氣態(tài)前體物的非均相或均相成核反應(yīng)等均為城市大氣細(xì)粒子的主要來(lái)源[7],其中交通源的貢獻(xiàn)受到更為廣泛的關(guān)注.研究表明,機(jī)動(dòng)車(chē)排放對(duì)北京大氣污染的貢獻(xiàn)可能在逐步增加[8-10].近年來(lái),我國(guó)頻繁發(fā)生的霾污染天氣主要集中在京津冀、長(zhǎng)江三角洲、珠江三角洲等經(jīng)濟(jì)較為發(fā)達(dá)和人口稠密的地區(qū).其中,京津冀地區(qū)[11-12]區(qū)域污染問(wèn)題突出,1980年后霾日數(shù)明顯增加,且多發(fā)生于冬季,12月和1月的霾日數(shù)約達(dá)到全年的30%.史軍等[13]對(duì)長(zhǎng)江三角洲霾氣候數(shù)據(jù)序列的統(tǒng)計(jì)結(jié)果表明,1961~2007整個(gè)長(zhǎng)江三角洲霾日數(shù)基本上都呈增長(zhǎng)趨勢(shì),并以杭州和南京增加最多.珠江三角洲的平均霾日數(shù)大值區(qū)則主要分布在廣東和廣西中東部,且1960~2007年霾日數(shù)呈現(xiàn)顯著的上升趨勢(shì),2008年后有所下降[14].霾的發(fā)生發(fā)展與污染氣體的排放和氣象條件密切相關(guān). Schichtel等[15]分析了美國(guó)1980~ 1995年的霾變化趨勢(shì),認(rèn)為SO2排放量的下降直接影響霾日數(shù).賀克斌等[16]指出,低空O3已逐漸成為大氣污染的核心物質(zhì)之一,O3的含量對(duì)顆粒物的生消也有重要的作用.同時(shí),高濃度的NO2對(duì)重污染的形成起到促進(jìn)作用[17].適宜的氣象條件,特別是靜穩(wěn)天氣和高濕環(huán)境,可以快速形成以PM2.5為典型特征的霾污染過(guò)程,特別是在冬季,重污染過(guò)程頻發(fā),污染物濃度呈現(xiàn)出典型的慢積累、快清除的“鋸齒形”變化特征[18].唐宜西等[19]研究發(fā)現(xiàn),靜小風(fēng)狀況以及逆溫結(jié)構(gòu)是北京地區(qū)霾產(chǎn)生的主要原因,而高濕則會(huì)加劇污染程度.張人禾等[20]指出,2013年我國(guó)中東部大范圍霾天氣中超過(guò)2/3的變化由氣象因子造成,且動(dòng)力和熱力因子貢獻(xiàn)相當(dāng).在霾過(guò)程和對(duì)應(yīng)的清潔過(guò)程中近地層溫度、風(fēng)場(chǎng)、位勢(shì)高度等均顯著不同[21].不同污染背景條件下,氣溶膠的粒徑分布、化學(xué)組分等存在顯著的差異,并決定了其不同的理化特性和環(huán)境效應(yīng)[22].如有機(jī)氣溶膠在清潔天平均貢獻(xiàn)PM1(亞微米細(xì)顆粒物)的約70%,而在重霾污染天,二次無(wú)機(jī)組分貢獻(xiàn)量顯著增加, 超過(guò)50%[23].北京地區(qū)的老化氣溶膠濃度水平以及粒子在積聚模態(tài)大量堆積的分布特點(diǎn)也是造成區(qū)域霾形成的重要原因,且以酸性氣溶膠的貢獻(xiàn)為主[24].近年來(lái),針對(duì)短時(shí)間空氣污染事件及其對(duì)人類健康影響的關(guān)注越來(lái)越多[25].煙花爆竹的集中燃放可產(chǎn)生大量的有機(jī)氣溶膠、重金屬離子、NO、SO2、CO等,可使空氣質(zhì)量在短時(shí)間內(nèi)急劇惡化[26].研究表明[27],煙花燃放時(shí)段氣溶膠數(shù)濃度多集中在100~500nm,并多伴有新粒子生成事件.本文對(duì)2013~2015北京全年及冬季時(shí)段PM2.5、氣體組分及氣溶膠譜分布結(jié)合氣象條件進(jìn)行系統(tǒng)的分析并對(duì)其中一次煙花污染事件進(jìn)行個(gè)例分析,探討大氣污染的生消規(guī)律,以期為防治大氣污染、保護(hù)城市生態(tài)環(huán)境提供科學(xué)依據(jù).

    1 資料來(lái)源

    根據(jù)北京市環(huán)境保護(hù)監(jiān)測(cè)中心網(wǎng)站(http://www.bjmemc.com.cn) 提供的監(jiān)測(cè)點(diǎn)信息,選取了2013至2015年北京城區(qū)國(guó)控環(huán)境監(jiān)測(cè)點(diǎn)(8個(gè)),對(duì)其提供的PM2.5、PM10、SO2、CO、O3、NO2的質(zhì)量濃度進(jìn)行了站點(diǎn)的平均,時(shí)間分辨率為1h.2013~2015年氣象因子數(shù)據(jù)來(lái)自國(guó)家氣象局北京城區(qū)加密觀測(cè)站點(diǎn)(4個(gè))的平均值,包括溫度、壓強(qiáng)、濕度、風(fēng)速、風(fēng)向、能見(jiàn)度,時(shí)間分辨率為1h.SMPS對(duì)氣溶膠數(shù)譜的觀測(cè)資料從2014年11月27日至2015年4月9日,時(shí)間分辨率為3min,測(cè)量粒徑范圍范圍14.6~661.2nm,由北京市氣象局提供.

    2 結(jié)果與討論

    2.1 2013~2015年空氣質(zhì)量整體概況

    2012年3月,國(guó)家環(huán)境保護(hù)部發(fā)布了新的《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)》(GB 3095—2012)[28],將PM2.5等大氣污染物列入了評(píng)價(jià)標(biāo)準(zhǔn).根據(jù)空氣質(zhì)量標(biāo)準(zhǔn),北京市2013~2015年分別有55%、55%和60%的時(shí)段達(dá)到國(guó)家細(xì)顆粒物空氣二級(jí)標(biāo)準(zhǔn)(PM2.5日平均濃度低于75μg/m3),其中空氣質(zhì)量為優(yōu)(PM2.5日平均濃度低于35μg/m3)的比例逐步升高;PM2.5的年平均濃度呈下降趨勢(shì),3年間下降~8.4%,且主要反映在中度以上污染(PM2.5日平均濃度低于115μg/m3)天數(shù)的減少,由2013年的99天下降至2015年的80d.冬季大氣污染的嚴(yán)重程度顯著高于全年平均,2013年冬季PM2.5的平均濃度高達(dá)132.7μg/m3,其中嚴(yán)重污染(PM2.5日平均濃度低于250μg/m3)的比例達(dá)16%. 2013~2015年冬季的整體空氣質(zhì)量也呈轉(zhuǎn)好的趨勢(shì),PM2.5質(zhì)量濃度的平均值呈下降趨勢(shì),達(dá)到國(guó)家細(xì)顆粒物空氣二級(jí)標(biāo)準(zhǔn)的天數(shù)顯著增多.氣象因素對(duì)于污染物的形成、演變、傳輸和消散等過(guò)程都起著至關(guān)重要的作用.氣象因子通過(guò)多種方式改變大氣顆粒物的發(fā)展過(guò)程,如:累積或稀釋,氣粒轉(zhuǎn)化和二次形成等[29].北京市區(qū)污染物濃度的升高通常伴隨弱偏南風(fēng)、高濕的大氣環(huán)境,當(dāng)風(fēng)速增大到3m/s后,對(duì)污染物有明顯的清除作用[29].表1對(duì)2013~2015年全年和冬季時(shí)段PM2.5、PM10、氣體組分和各氣象要素進(jìn)行了統(tǒng)計(jì).除2014年冬季與2013年冬季的平均風(fēng)速持平略有升高外,2013~2015年風(fēng)速的平均值均呈下降趨勢(shì).就全年趨勢(shì)來(lái)看,2013~2015年風(fēng)速的均值略有下降,其中小風(fēng)(<1.5m/s)的比例有所升高,大于3m/s以上風(fēng)速的比例有所降低,同時(shí),偏南風(fēng)在全風(fēng)向中所占比例也呈上升趨勢(shì);相對(duì)濕度、降水、溫度、氣壓等氣象要素均值的年際變化差異不大,綜合上述氣象因素,2013~2015年大氣易出現(xiàn)污染物累積的靜穩(wěn)天氣條件有所升高,但PM2.5質(zhì)量濃度顯著降低,空氣質(zhì)量呈明顯改善,污染減排的貢獻(xiàn)功不可沒(méi).對(duì)于北京地區(qū)冬季,特別是2014~2015年,各類氣象條件均有利于污染物的累積,但輕度以上污染出現(xiàn)的概率顯著降低,也體現(xiàn)出減排的功效.在顆粒物質(zhì)量濃度有所改善的同時(shí),O3的質(zhì)量濃度有所升高,在關(guān)注顆粒物污染的同時(shí),以高濃度O3為典型特征的光化學(xué)污染問(wèn)題也逐漸顯現(xiàn).

    圖1 2013~2015年全年(a)和冬季時(shí)段(b)的空氣質(zhì)量頻數(shù)

    2.2 日變化趨勢(shì)

    硫酸鹽、硝酸鹽等無(wú)機(jī)物是大氣細(xì)粒子的重要組成部分,而這些無(wú)機(jī)鹽的生成主要通過(guò)SO2、NO2等氣態(tài)前體物的光化學(xué)反應(yīng)生成,因而大氣的氧化性對(duì)于PM2.5的濃度有較大影響,將大氣氧化劑 OX(NO2+O3)作為評(píng)價(jià)大氣氧化能力的指標(biāo)[30],其濃度的變化能夠反映區(qū)域及局地大氣的氧化性特征[31].由表1知,冬季O3濃度顯著低于全年均值,大氣氧化能力低于年均值,其中2014年偏低最多,冬季OX較全年平均偏低~18%,大氣的氧化性在16:00左右達(dá)到極大值,在日出前后達(dá)到極小值.由于冬季的大氣污染程度重于其它季節(jié),本文對(duì)2013~2015年各年冬季的污染性氣體、PM2.5質(zhì)量濃度和大氣的氧化性進(jìn)行了日變化對(duì)比(圖2).有研究表明,由機(jī)動(dòng)車(chē)尾氣排放導(dǎo)致城市氮氧化物分擔(dān)率達(dá)39%~79%[32],受交通源影響顯著,NO2、CO的日變化曲線相似,均在早高峰(08:00左右)出現(xiàn)峰值,在16:00左右達(dá)到一日的最低濃度,NO2在晚高峰(20:00左右)也存在顯著的峰值.冬季受采暖影響,SO2濃度高于全年平均近2倍.2013~2015年冬季,SO2均在09:00和13:00前后出現(xiàn)峰值,09:00左右的峰值高于13:00左右,且峰值強(qiáng)度逐年減弱.2013~ 2015年,氣體污染物SO2、CO、NO2都呈逐年遞減的趨勢(shì),特別是2014~2015年,SO2和NO2濃度下降顯著,SO2的日均濃度下降為上一年同期的1/2.由于白天光化學(xué)反應(yīng)速率加大,O3濃度迅速升高,下午16:00左右O3濃度達(dá)到最大值,且2015年冬季的O3濃度值較2013和2014年有所升高.PM2.5濃度受邊界層變化影響表現(xiàn)為夜間高白天低的日變化規(guī)律,均在上午10:00左右和夜間23:00左右出現(xiàn)兩個(gè)峰值,23:00左右出現(xiàn)的峰值濃度為全天最高,且這一峰值濃度自2013至2015年逐年降低.PM2.5峰值濃度的出現(xiàn)時(shí)段并未對(duì)應(yīng)交通高峰時(shí)段,而是延后2h左右,說(shuō)明機(jī)動(dòng)車(chē)對(duì)PM2.5濃度的影響作用主要表現(xiàn)為機(jī)動(dòng)車(chē)排放的氣態(tài)污染物發(fā)生化學(xué)反應(yīng)生成PM2.5中的二次粒子,而非機(jī)動(dòng)車(chē)的一次排放.

    表1 2013~2015年全年及冬季的PM2.5、PM10、氣體組分和各氣象要素的平均值統(tǒng)計(jì)

    機(jī)動(dòng)車(chē)排放可以產(chǎn)生大量的一氧化碳和氮氧化物.本文通過(guò)對(duì)2013~2015年冬季的NO2和CO在早高峰時(shí)段和整個(gè)白天時(shí)段(07:00~18:00)的質(zhì)量濃度進(jìn)行相關(guān)性分析(圖3)探討交通源的作用.早高峰時(shí)段選取07:00~09:00,且排除降水日和周末的影響.早高峰時(shí)段和整個(gè)白天時(shí)段的質(zhì)量濃度相關(guān)性較好,相關(guān)系數(shù)(r)均大于0.6.2013年CO早高峰時(shí)段的濃度較整個(gè)白天時(shí)段增加10%,2014年增大至15%,2015年更是增加到18%;同時(shí),NO2在交通早高峰時(shí)段的貢獻(xiàn)也呈現(xiàn)逐年遞增的趨勢(shì),3年內(nèi)增大了約4%.因而交通源對(duì)污染性氣體的貢獻(xiàn)不能忽略,且呈現(xiàn)出逐漸加重的趨勢(shì),交通源已成為大城市環(huán)境污染的重要來(lái)源之一.

    圖2 2013~2015年各年冬季的污染性氣體、PM2.5質(zhì)量濃度和大氣的氧化性O(shè)X的日變化

    圖3 2013~2015年冬季NO2和CO在早高峰時(shí)段和整個(gè)白天時(shí)段的質(zhì)量濃度的相關(guān)關(guān)系

    2.3 氣象要素對(duì)污染物影響的滯后性分析

    氣象因素作為影響環(huán)境空氣質(zhì)量的一個(gè)重要條件,其與顆粒物濃度的相關(guān)性已經(jīng)達(dá)到了廣泛認(rèn)可[33].氣象要素對(duì)污染物的生消過(guò)程有顯著影響,但影響作用并非同步,通常需要一個(gè)時(shí)間過(guò)程,即氣象要素對(duì)污染物的累積和擴(kuò)散作用有一定的時(shí)間滯后性[34].各類氣象要素中,風(fēng)速通過(guò)影響污染物的水平輸送擴(kuò)散能力對(duì)污染物的質(zhì)量濃度有重要影響,通常大于3m/s的風(fēng)速對(duì)污染物濃度有顯著的清除作用[29];相對(duì)濕度的增高易導(dǎo)致細(xì)粒子吸濕性增長(zhǎng),使得污染物質(zhì)量濃度升高;而溫度則是通過(guò)影響污染物的垂直擴(kuò)散能力對(duì)污染物的質(zhì)量濃度產(chǎn)生影響,當(dāng)大氣出現(xiàn)逆溫時(shí),層結(jié)狀態(tài)較穩(wěn)定,大氣垂直運(yùn)動(dòng)能力差,有利于污染物的積聚.為研究濕度、風(fēng)速和溫度對(duì)PM2.5質(zhì)量濃度影響的時(shí)間效應(yīng),本文給出了2013~2015年全年及冬季的PM2.5濃度與風(fēng)速、溫度滯后1~10h內(nèi)的滯后相關(guān)系數(shù)(表2).由表2可見(jiàn),相對(duì)濕度和溫度與PM2.5濃度呈正相關(guān)性,而風(fēng)速呈負(fù)相關(guān)性.在整個(gè)研究時(shí)段和冬季時(shí)段,相對(duì)濕度和風(fēng)速均在改變1h后和PM2.5濃度的相關(guān)系數(shù)達(dá)到極大值,隨后逐漸降低,其中,冬季PM2.5對(duì)相對(duì)濕度的響應(yīng)顯著的高于全年平均,滯后1h的相對(duì)濕度和PM2.5質(zhì)量濃度的相關(guān)系數(shù)可以達(dá)到0.625.有研究表明,SO2在高濕的狀態(tài)下極易通過(guò)化學(xué)反應(yīng)轉(zhuǎn)化為硫酸鹽并進(jìn)一步影響大氣酸度,使得大氣中酸性粒子增加[35],受冬季燃煤影響,SO2的質(zhì)量濃度為全年平均的2~3倍,因而冬季相對(duì)濕度增大后,硫酸鹽在PM2.5中所占的比例也迅速升高(2.9 μg/m3/10%RH)[36].對(duì)于溫度的變化,PM2.5濃度與溫度的相關(guān)系數(shù)隨滯后時(shí)間的增加逐漸增大,冬季時(shí)段,在滯后7h后相關(guān)系數(shù)達(dá)到最大值,全年的觀測(cè)在滯后8h后相關(guān)系數(shù)最大,而后隨滯后時(shí)效繼續(xù)延長(zhǎng)相關(guān)系數(shù)均逐漸降低.綜上,相對(duì)濕度和風(fēng)速對(duì)PM2.5濃度影響的響應(yīng)相對(duì)較快,1h左右對(duì)PM2.5濃度的影響達(dá)到最大,且對(duì)冬季的影響高于全年平均,而PM2.5對(duì)溫度變化的響應(yīng)相對(duì)較慢,一般在改變溫度的7~8h后影響最大.在進(jìn)行污染物濃度的預(yù)報(bào)預(yù)警時(shí),應(yīng)充分考慮不同氣象要素影響的滯后性.

    表2 PM2.5與濕度(RH)、風(fēng)速(WS)、溫度(T)滯后相關(guān)系數(shù)(r)

    2.3 氣溶膠粒徑譜分布

    為研究不同污染背景條件下氣溶膠數(shù)濃度的分布,依據(jù)PM2.5的日均濃度,將SMPS觀測(cè)期段(2014年11月27日至2015年4月9日)的大氣分為清潔日(PM2.5的日均濃度<75μg/m3)和污染日(PM2.5的日均濃度≥75μg/m3),分別計(jì)算清潔天和污染天條件下氣溶膠數(shù)濃度譜分布(圖4a),計(jì)算過(guò)程去掉了降水日和沙塵日(判別條件:日均能見(jiàn)度低于10km、PM10日均濃度高于150μg/m3、日均風(fēng)速高于3m/s、無(wú)降水).SMPS的觀測(cè)時(shí)段共有34d為清潔天,55d為污染天.本文劃分核模態(tài)粒子、愛(ài)根核模態(tài)粒子及積聚核模態(tài)粒子的粒徑范圍分別為14.6~25.0nm、25.0~ 100.0nm及100.0~661.2nm.由圖4a可見(jiàn),污染時(shí)期較清潔時(shí)期氣溶膠的粒徑分布均有不同程度地向大粒子一側(cè)拓寬的現(xiàn)象,清潔日時(shí),峰區(qū)直徑集中在3~5nm,對(duì)應(yīng)的數(shù)濃度峰值普遍低于20×103個(gè)/cm3,污染天時(shí),峰區(qū)直徑增大至100~ 200nm,積聚核模態(tài)粒子較清潔天時(shí)顯著的增多(平均占總數(shù)濃度的55%),由于積聚核模態(tài)粒子的粒徑更接近可見(jiàn)光波長(zhǎng)范圍而具有強(qiáng)消光效率[23],因而這一模態(tài)粒子濃度的增大可能是導(dǎo)致霾天氣能見(jiàn)度急劇下降的重要原因.SMPS觀測(cè)期間共有9個(gè)降水日,日降水量為0.3~16.2mm,見(jiàn)圖4b.降水對(duì)氣溶膠粒子的清除作用與降水量的大小有直接的關(guān)系,如果大氣處于污染狀態(tài),在降水量較小時(shí),降水的清除作用不顯著,而造成的低層增濕作用有利于氣溶膠粒子吸濕增長(zhǎng),反而使得污染加劇,如2015年1月24日,北京城區(qū)平均出現(xiàn)0.3mm的日降水量,PM2.5的日均濃度高達(dá)219.6μg/m3,較前一日的日均濃度增大了一倍多.只有降水達(dá)到一定的量級(jí),對(duì)氣溶膠粒子才能表現(xiàn)出顯著的清除作用.在降水量級(jí)較小時(shí)(<5mm),通常無(wú)明顯的清除作用,氣溶膠粒徑的峰值與污染天相近,位于100nm附近;而隨降水量級(jí)逐漸增大(35mm),氣溶膠峰區(qū)直徑逐漸向積聚核模態(tài)方向發(fā)展,在本次研究降水量級(jí)最大日(16.2mm),峰值濃度對(duì)應(yīng)的直徑接近200nm.有研究表明,核模態(tài)粒子主要由于布朗擴(kuò)散機(jī)制所主導(dǎo)而被雨滴清除,積聚態(tài)粒子主要受布朗擴(kuò)散和湍流擴(kuò)散的共同作用而被雨滴清除[37].雨天積聚核模態(tài)的粒子減少率在三個(gè)模態(tài)粒子中最低,直徑位于0.2~2μm 之間的區(qū)域通常稱“greenfield gap”[38],很難被雨滴有效清除,因而即便降水量級(jí)較大,但這一區(qū)間粒子的清除率低,使得殘留的氣溶膠粒徑向這一區(qū)間積聚.

    圖4 (a)污染天和清潔天氣溶膠數(shù)濃度譜分布;(b)不同降水天氣條件下氣溶膠數(shù)濃度譜分布

    2.4 2015年污染事件個(gè)例分析

    北京地區(qū)冬季存在多次污染物累積和消散的過(guò)程,在相對(duì)濕度較大,風(fēng)速較低的靜穩(wěn)天氣等不利于污染物擴(kuò)散、稀釋的積累過(guò)程中污染物濃度迅速累積,當(dāng)產(chǎn)生較大規(guī)模的降水和(或)北、西北方向較為潔凈、干燥的空氣團(tuán)對(duì)污染物稀釋作用時(shí),污染物濃度迅速的降低.本文選取了2015年2月18至2015年2月21日的一次污染過(guò)程進(jìn)行個(gè)例分析.2月18日06:00起,靜穩(wěn)天氣形勢(shì)逐漸建立,NO2、CO、SO2等氣體及PM2.5的質(zhì)量濃度均呈現(xiàn)波動(dòng)上升趨勢(shì),2月19日凌晨為除夕夜,受煙花燃放作用影響,污染物累積速率較快,2月19日06:00,PM2.5濃度達(dá)到347μg/m3.19日11時(shí)起,偏南風(fēng)向的風(fēng)速短時(shí)間有所增大,由于較大風(fēng)速的持續(xù)時(shí)間較短,且仍為偏南風(fēng)向,PM2.5濃度有所降低,但降低程度有限,濃度仍維持在100μg/m3左右.19日午后起,相對(duì)濕度逐漸增大,PM2.5開(kāi)始緩慢的累積,上升速率相對(duì)較小,其中PM2.5的增長(zhǎng)率平均~3.2μg/m3/h,21日13時(shí)PM2.5的峰值濃度達(dá)到316μg/m3;21日午后隨風(fēng)速增大和風(fēng)向的轉(zhuǎn)變,PM2.5、NO2、CO、SO2等污染物得到徹底清除,清除過(guò)程僅用了3h就使得PM2.5的質(zhì)量濃度降低到50μg/m3以下,O3的變化趨勢(shì)相反,在其它污染物濃度迅速降低的同時(shí)O3的濃度迅速升高.污染累積較重的時(shí)段(21日15時(shí))500hPa高度場(chǎng)(圖略)受低壓橫槽后西北氣流控制,高空少云,地面夜間有顯著的輻射降溫,在近地層形成逆溫層結(jié),同時(shí),850hPa高度場(chǎng)(圖略)中暖脊的形成使逆溫層厚度向上延伸,在北京等地形成干暖蓋結(jié)構(gòu),水汽和污染物在近地面逐漸積聚使得污染物濃度累積到較高的濃度.

    煙花爆竹的集中燃放可在短時(shí)間導(dǎo)致空氣質(zhì)量的急劇惡化[39-40].2月19日為傳統(tǒng)的除夕佳節(jié),對(duì)應(yīng)一次煙花燃放過(guò)程.有研究表明,煙花燃放時(shí)段積聚模態(tài)100~500nm粒徑段的氣溶膠數(shù)濃度顯著的升高[41-42].本文依據(jù)積聚核模態(tài)粒子濃度變化趨勢(shì)(圖7a)分為煙花燃放集中期(1月19日00:00~04:00),煙花燃放前(1月18日12:00~ 16:00)和煙花燃放后(1月19日13:00~17:00).在煙花集中燃放的時(shí)段,PM2.5、SO2、CO和NO2都出現(xiàn)了顯著的峰值,PM2.5的質(zhì)量濃度在18日23:00~19日01:00的2h內(nèi)由109μg/m3升高到292μg/m3,SMPS測(cè)得的氣溶膠數(shù)濃度也出現(xiàn)了顯著的峰值.煙花對(duì)污染物的貢獻(xiàn)可以通過(guò)選取煙花排放時(shí)間段首末兩點(diǎn)連成直線,與軸形成一個(gè)梯形,對(duì)梯形面積進(jìn)行積分,再對(duì)整個(gè)曲線與軸圍成的面積進(jìn)行積分,二者之間的區(qū)域(見(jiàn)圖5b中灰色與紅色曲線圍城的區(qū)域)近似可以認(rèn)為是煙花燃放的貢獻(xiàn)[43],通過(guò)此方法計(jì)算得到煙花對(duì)氣溶膠數(shù)濃度的貢獻(xiàn)為~26.1%.對(duì)煙花燃放前后和煙花燃放時(shí)段的氣溶膠數(shù)濃度譜進(jìn)行對(duì)比(圖6b),煙花燃放集中時(shí)段較燃放前后氣溶膠峰區(qū)直徑顯著的向增大的方向移動(dòng),積聚模態(tài)粒子和愛(ài)根核模態(tài)粒子出現(xiàn)顯著的峰值,核模態(tài)粒子濃度變化不大.根據(jù)新粒子生成事件主要判定標(biāo)準(zhǔn)[44]:(1)顆粒物的數(shù)譜分布中出現(xiàn)一個(gè)新的模態(tài);(2)這個(gè)新的模態(tài)從核模態(tài)粒徑范圍開(kāi)始出現(xiàn);(3)新的模態(tài)在大氣中存在一段時(shí)間;(4)新模態(tài)呈現(xiàn)增長(zhǎng)趨勢(shì),2月18日06:00~2月19日06:00的污染物累積過(guò)程伴隨一次新粒子的生成過(guò)程.

    為了探討不同方向、流速的氣團(tuán)對(duì)氣溶膠質(zhì)量濃度的影響,本文選取了2月17~2月21日污染過(guò)程中風(fēng)速-風(fēng)向-氣溶膠數(shù)濃度及風(fēng)向-風(fēng)速-PM2.5質(zhì)量濃度進(jìn)行分析,期間,PM2.5的平均質(zhì)量濃度為122.2μg/m3,平均數(shù)濃度為12082個(gè)/cm3.如圖7,圖中不同顏色代表氣溶膠質(zhì)量濃度,不同方向、距離圓心的長(zhǎng)度分別代表風(fēng)向和風(fēng)速.氣溶膠的數(shù)濃度和PM2.5質(zhì)量濃度的高值區(qū)度出現(xiàn)在風(fēng)速較小(<1.5m/s)時(shí)及中等風(fēng)速(1.5~3m/s)的偏南風(fēng)向,此次污染事件是局地污染和來(lái)自偏南方向氣團(tuán)輸送作用的共同結(jié)果.

    圖5 2015年2月18至2015年2月21日污染過(guò)程(a)NO2、O3、CO、SO2;(b)氣溶膠總數(shù)濃度(灰色)及去掉煙花影響后總數(shù)濃度(紅色);(c)相對(duì)濕度(RH)和溫度(Temp);(d)風(fēng)速(WS)、風(fēng)向(WD);(e)氣溶膠數(shù)濃度譜分布;(f)PM2.5和能見(jiàn)度(Vis)的時(shí)間序列

    圖6 (a)2015年2月17~2月21日核模態(tài)、愛(ài)根核模態(tài)和積聚核模態(tài)粒子數(shù)濃度的時(shí)間序列;(b)煙花燃放前(AFW)、后(BFW)及煙花燃放集中時(shí)段(FW)氣溶膠數(shù)濃度譜分布

    3 結(jié)論

    3.1 2013~2015年大氣易出現(xiàn)污染物累積的靜穩(wěn)天氣條件有所升高,但PM2.5質(zhì)量濃度顯著降低,空氣質(zhì)量呈明顯改善,污染減排有一定成效.交通源對(duì)污染性氣體的貢獻(xiàn)不能忽略,且呈現(xiàn)出逐漸加重的趨勢(shì),交通源已成為大城市環(huán)境污染的重要來(lái)源之一.

    3.2 相對(duì)濕度和風(fēng)速對(duì)PM2.5濃度影響的相應(yīng)相對(duì)較快,1h左右對(duì)PM2.5濃度的影響達(dá)到最大,且冬季的影響高于全年均值,而PM2.5的變化對(duì)溫度的響應(yīng)相對(duì)較慢,一般在改變溫度的7~8h后對(duì)污染物濃度的影響最大.

    3.3 污染時(shí)期較清潔時(shí)段氣溶膠的粒徑分布均有不同程度地向大粒子一側(cè)拓寬的現(xiàn)象,清潔日時(shí),峰區(qū)直徑集中在3~5nm,污染天時(shí),峰區(qū)直徑增大至100~200nm;隨降雨量級(jí)逐漸增大,氣溶膠峰值濃度逐漸向積聚核模態(tài)方向發(fā)展.

    3.4 2015年2月18~2015年2月21日污染事件中,煙花燃放時(shí)段煙花對(duì)氣溶膠數(shù)濃度的貢獻(xiàn)為26.1%,且伴隨新粒子生成.煙花燃放時(shí)段氣溶膠峰值濃度向粒徑偏大的方向移動(dòng),積聚模態(tài)粒子和愛(ài)根核模態(tài)粒子出現(xiàn)顯著的峰值,而核模態(tài)粒子濃度變化不大.

    [1] 吳 兌,陳慧忠,吳 蒙,等.三種霾日統(tǒng)計(jì)方法的比較分析[J]. 中國(guó)環(huán)境科學(xué), 2014,34(3):545-554.

    [2] 唐孝炎,張遠(yuǎn)航,邵 敏.大氣環(huán)境化學(xué)[M]. 2版.北京:高等教育出版社, 2006.

    [3] 朱 彤.城市與區(qū)域大氣復(fù)合污染[M]. 北京:化學(xué)工業(yè)出版社, 2005.

    [4] Parrish D D, Zhu T. Clean Air for Megacities [J]. Science, 2009,326(5953):674-675.

    [5] Zhang Q H, Zhang J P, Xue H W. The challenge of improving visibility in Beijing [J]. Atmosphericchemistry and physics, 2010,10(16):7821-7827.

    [6] Sun Y, Zhuang G, Tang A, et al. Chemical characteristics of PM2.5and PM10in haze-fog episodes in Beijing [J]. Environmental Science and Technology, 2006,40:3148-3155.

    [7] Laakso L, Hussein T, Aarni P, et al. Diurnal and annual characteristics of particle mass and number concentrations in urban, rural and Arctic environments in Finland [J]. Atmospheric Environment, 2003,37(19):2629-2641.

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

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

    [10] 楊孝文,周 穎,程水源,等.北京冬季一次重污染過(guò)程的污染特征及成因分析 [J]. 中國(guó)環(huán)境科學(xué), 2016,36(3):679-686.

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

    [12] 許亞宣,李小敏,于華通.等.邯鄲市大氣復(fù)合污染特征的監(jiān)測(cè)研究[J]. 環(huán)境科學(xué)學(xué)報(bào), 2015,35(19):2710-2722.

    [13] 史 軍,崔林麗.長(zhǎng)江三角洲城市群霾的演變特征及影響因素研究[J]. 中國(guó)環(huán)境科學(xué), 2013,33(12):2113-2122.

    [14] 符傳博,丹 麗,唐家翔,等.1960~2013年華南地區(qū)霾污染的時(shí)空變化及其與關(guān)鍵氣候因子的關(guān)系[J]. 中國(guó)環(huán)境科學(xué), 2016, 36(5):1313-1322.

    [15] Schichtel B A, Husar R B, Falke S R, et al. Haze trends over the United States, 1980–1995 [J]. Atmospheric Environment, 2001,35(30):5205-5210.

    [16] 賀克斌,楊復(fù)沫,段鳳魁.大氣顆粒物與區(qū)域復(fù)合污染[J]. 北京:科學(xué)出版社, 2011.

    [17] 程念亮,李云婷,張大偉,等.2013~2014年北京市NO2時(shí)空分布研究[J]. 中國(guó)環(huán)境科學(xué), 2016,36(1):18-26.

    [18] Jia Y, Rahn K A, He K, et al. A novel technique for quantifying the regional component of urban aerosol solely from its sawtooth cycles [J]. Journal of Geophysical Research, 2008,113(D21).

    [19] 唐宜西,張小玲,熊亞軍,等.北京一次持續(xù)霾天氣過(guò)程氣象特征分析[J]. 氣象與環(huán)境學(xué)報(bào), 2013,29(5):12-19.

    [20] 張人禾,李 強(qiáng),張若楠,等.2013年1月中國(guó)東部持續(xù)性強(qiáng)霧霾天氣產(chǎn)生的氣象條件分析[J]. 中國(guó)科學(xué):地球科學(xué), 2014,(1): 27-36.

    [21] 史 軍,崔林麗.長(zhǎng)江三角洲城市群霾的演變特征及影響因素研究 [J]. 中國(guó)環(huán)境科學(xué), 2013,33(12):2113-2122.

    [22] Tang U. Size, composition, and mixing state of individual aerosol particles in a South China coastal city [J]. Journal of Environmental Sciences, 2010,22(4):561-569.

    [23] 江 琪,孫業(yè)樂(lè),王自發(fā),等.應(yīng)用顆粒物化學(xué)組分監(jiān)測(cè)儀(ACSM)實(shí)時(shí)在線測(cè)定致霾細(xì)粒子無(wú)機(jī)和有機(jī)組分[J]. 科學(xué)通報(bào), 2013,58(36):3818-3828.

    [24] Zhang Y M, Sun J Y, Zhang X Y, et al. Seasonal characterization of components and size distributions for submicron aerosols in Beijing [J]. Science China Earth Sciences, 2013,56(5):890-900.

    [25] III C A P, Burnett R T, Thun M J, et al. Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution [J]. Jama the Journal of the American Medical Association, 2002,287(9):1132-1141.

    [26] Vecchi R, Bernardoni V, Cricchio D, et al. The impact of fireworks on airborne particles [J]. Atmospheric Environment, 2008,42(6):1121-1132.

    [27] Zhang M, Wang X, et al. Physical characterizationof aerosol particles during the Chinese New Year’s firework events [J]. Atmospheric Environment, 2010,44(39):5191-5198.

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

    [29] Sun Y L, Wang Z F, Fu P Q, et al. Aerosol composition, sources and processes during wintertime in Beijing, China [J]. Atmospheric Chemistry and Physics Discussions, 2013,13(1): 2077-2123.

    [30] Stephens S, Madronich S, Wu F, et al. Weekly patterns of Mexico City's surface concentrations of CO, NO, PM10and O3during 1986~2007 [J]. Atmospheric Chemistry and Physics, 2008,8(17): 5313-5325.

    [31] 王占山,潘麗波.火電廠大氣污染物排放標(biāo)準(zhǔn)實(shí)施效果的數(shù)值模擬研究[J]. 環(huán)境科學(xué), 2014,35(3):853-863.

    [32] 郝吉明,吳 燁,傅立新,等.北京市機(jī)動(dòng)車(chē)污染分擔(dān)率的研究[J]. 環(huán)境科學(xué), 2001,22(5):1-6.

    [33] 趙晨曦,王云琦,王玉杰,等.北京地區(qū)冬春PM2. 5和PM10污染水平時(shí)空分布及其與氣象條件的關(guān)系[J]. 環(huán)境科學(xué), 2014,35(2): 418-426.

    [34] 肖 嫻,范紹佳,蘇 冉.2011年10月珠江三角洲一次區(qū)域性空氣污染過(guò)程特征分析[J]. 環(huán)境科學(xué)學(xué)報(bào), 2014,34:290-296.

    [35] Kaul D S, Gupta T, Tripathi S N, et al. Secondary organic aerosol: a comparison between foggy and nonfoggy days [J]. Environment Science & Technology, 2011,45(17):7307-7313.

    [36] Sun Y, Wang Z, Fu P, et al. The impact of relative humidity on aerosol composition and evolution processes during wintertime in Beijing, China [J]. Atmospheric Environment, 2013,77:927-934.

    [37] 尚 倩,李子華,楊 軍,等.南京冬季大氣氣溶膠粒子譜分布及其對(duì)能見(jiàn)度的影響[J]. 環(huán)境科學(xué), 2011,32(9):2750-2760.

    [38] Greenfield S M. Rain Scavenging of Radioactive Particulate Matter from the Atmosphere [J]. Journal of the Atmospheric Sciences, 2010,14(14):115-125.

    [39] Shi Y, Zhang N, Gao J, et al. Effect of fireworks display on perchlorate in air aerosols during the Spring Festival [J]. Atmospheric Environment, 2011,45(6):1323-1327.

    [40] Attri A K, Kumar U, Jain V K. Microclimate: Formation of ozone by fireworks [J]. Nature, 2001,411(411):1015-1015.

    [41] Zhang M, Wang X, Chen J, et al. Physical characterization of aerosol particles during the Chinese New Year’s firework events [J]. Atmospheric Environment, 2010,44(39):5191-5198.

    [42] Rissler J, Vestin A, Swietlicki E, et al. Size distribution and hygroscopic properties of aerosol particles from dry-season biomass burning in Amazonia [J]. Atmospheric Chemistry & Physics, 2006,6(2):471-491.

    [43] Jiang Q, Sun Y L, Wang Z, et al. Aerosol composition and sources during the Chinese Spring Festival: fireworks, secondary aerosol, and holiday effects [J]. Atmospheric Chemistry & Physics, 2014,14(14):20617-20646.

    [44] Maso M D, Kulmala M, Riipinen I, et al. Formation and growth of fresh atmospheric aerosols: Eight years of aerosol size distribution data from SMEAR II, Hyyti?l?, Finland [J]. Boreal Environment Research, 2005,10(5):323-336.

    The characteristics of PM2.5、Reactive gas and aerosol Size Distributions of Beijing from 2013 to 2015.

    JIANG Qi1, WANG Fei2*, ZHANG Heng-de1, WANG Fei3, Lü Meng-yao1, HE Jia-bao4

    (1. National Meteorological Centre, Beijing 100081, China;2.Chinese Academy of Meteorological Sciences, Beijing 100081, China;3.Beijing Municipal Weather Modification Office, Beijing 100089, China;4.Environment Monitoring Center of Ningbo 315012, China)., 2017,37(10):3647~3657

    PM2.5, gas composition, aerosol spectral distribution and meteorological conditions were systematically analyzed during 2013~2015 in Beijing. Results showed that the PM2.5mass concentration decreased significantly from 2013 to 2015, and the air quality improved significantly. The contribution of traffic sources to air pollution can not be neglected and tended to increase gradually. The influence of humidity and wind speed on PM2.5concentration was relatively fast and the biggest influence on PM2.5concentration was about 1h time lag. While the PM2.5response to temperature changes was relatively slow, generally after time lag 7~8h before the effect reached to the maximum. The particle size distribution of aerosols moved to the larger size in pollution days than cleaning days. On cleaning days, the diameter of the peak area was concentrated at 3~5nm, and increased to 100~200nm on pollution days. With the increasing of precipitation level, the peak concentration of aerosol gradually developed towards to the accumulation mode. During the pollution process, from February 18 to February 21 in 2015, the contribution of fireworks to aerosol concentration was 26.1% in fireworks time, and accompanied with new particle generation. The discharged of fireworks produced significant peaks in the Accumulation mode and the Aitken mode particles, the Nuclear mode particle concentration changed not much.

    PM2.5;reactive gas;meteorological elements;aerosol particle size spectrum;fireworks

    X51

    A

    1000-6923(2017)10-3647-11

    江 琪(1989-),女,河北唐山人,碩士,主要從事大氣物理及大氣化學(xué)方面的研究.發(fā)表論文10余篇.

    2017-02-27

    科技部大氣污染專項(xiàng)課題(2016YFC0203301);公益性行業(yè)(環(huán)保)專項(xiàng)(201509001);公益性行業(yè)(氣象)科研專項(xiàng)(GYHY201306015);國(guó)家科技支撐計(jì)劃課題(2014BAC16B02)

    * 責(zé)任作者, 工程師, wf2775@sina.com

    猜你喜歡
    大氣風(fēng)速模態(tài)
    大氣的呵護(hù)
    軍事文摘(2023年10期)2023-06-09 09:15:06
    基于Kmeans-VMD-LSTM的短期風(fēng)速預(yù)測(cè)
    基于最優(yōu)TS評(píng)分和頻率匹配的江蘇近海風(fēng)速訂正
    大氣古樸揮灑自如
    大氣、水之后,土十條來(lái)了
    基于GARCH的短時(shí)風(fēng)速預(yù)測(cè)方法
    國(guó)內(nèi)多模態(tài)教學(xué)研究回顧與展望
    考慮風(fēng)速分布與日非平穩(wěn)性的風(fēng)速數(shù)據(jù)預(yù)處理方法研究
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識(shí)別
    由單個(gè)模態(tài)構(gòu)造對(duì)稱簡(jiǎn)支梁的抗彎剛度
    国产色视频综合| 亚洲欧美一区二区三区黑人| 美女午夜性视频免费| 精品福利观看| 精品国产乱子伦一区二区三区| 黄色片一级片一级黄色片| 嫩草影院精品99| 美女午夜性视频免费| 在线观看www视频免费| 精品卡一卡二卡四卡免费| 久久久久国产一级毛片高清牌| 日韩av在线大香蕉| 99精品欧美一区二区三区四区| 神马国产精品三级电影在线观看 | 大陆偷拍与自拍| av欧美777| 亚洲成人久久性| 中文字幕精品免费在线观看视频| 日韩欧美在线二视频| 亚洲一卡2卡3卡4卡5卡精品中文| 叶爱在线成人免费视频播放| 美女免费视频网站| 亚洲国产高清在线一区二区三 | 久99久视频精品免费| 日日干狠狠操夜夜爽| 黄色女人牲交| 悠悠久久av| 超碰成人久久| 国产精品一区二区免费欧美| 亚洲欧美激情综合另类| 免费人成视频x8x8入口观看| 免费高清视频大片| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧美日韩另类电影网站| 国产成人免费无遮挡视频| 亚洲欧美精品综合久久99| 制服诱惑二区| 日本免费a在线| 亚洲熟妇中文字幕五十中出| 成人三级黄色视频| 日韩中文字幕欧美一区二区| 国产精品国产高清国产av| 日韩国内少妇激情av| 性欧美人与动物交配| 午夜福利视频1000在线观看 | 国产亚洲精品综合一区在线观看 | АⅤ资源中文在线天堂| 99re在线观看精品视频| ponron亚洲| 变态另类成人亚洲欧美熟女 | 欧美日韩中文字幕国产精品一区二区三区 | 咕卡用的链子| xxx96com| 精品午夜福利视频在线观看一区| 久久香蕉激情| 日本 欧美在线| 嫩草影院精品99| 国产成年人精品一区二区| 久久国产精品男人的天堂亚洲| 成人av一区二区三区在线看| 亚洲色图av天堂| 久久久久久久久中文| 国产精品乱码一区二三区的特点 | 日本 欧美在线| 久久青草综合色| 精品午夜福利视频在线观看一区| 午夜精品在线福利| 国产高清有码在线观看视频 | 久久精品91无色码中文字幕| 国产精品九九99| av欧美777| 欧美+亚洲+日韩+国产| 日本免费一区二区三区高清不卡 | 国产精品爽爽va在线观看网站 | 欧美激情久久久久久爽电影 | 久久精品国产清高在天天线| 最好的美女福利视频网| 亚洲电影在线观看av| 黄色成人免费大全| 黑人欧美特级aaaaaa片| 欧美色视频一区免费| 国产精品爽爽va在线观看网站 | 99国产精品99久久久久| 夜夜夜夜夜久久久久| 色播亚洲综合网| 国产精品一区二区三区四区久久 | 久久国产精品影院| 亚洲专区字幕在线| 午夜福利高清视频| 免费少妇av软件| 99精品欧美一区二区三区四区| 日本欧美视频一区| 欧美成人午夜精品| 久久久久久久久中文| 亚洲avbb在线观看| 久久九九热精品免费| 看片在线看免费视频| 国产极品粉嫩免费观看在线| 亚洲精品中文字幕一二三四区| 十八禁网站免费在线| 国产一卡二卡三卡精品| 午夜a级毛片| 99国产极品粉嫩在线观看| 看黄色毛片网站| 午夜老司机福利片| 18禁美女被吸乳视频| 99精品欧美一区二区三区四区| 久久人人精品亚洲av| 国产99白浆流出| 女同久久另类99精品国产91| 亚洲午夜精品一区,二区,三区| 国产亚洲精品一区二区www| 欧美大码av| 欧美久久黑人一区二区| 黑人欧美特级aaaaaa片| 级片在线观看| 91大片在线观看| 久久午夜综合久久蜜桃| 色尼玛亚洲综合影院| 色综合欧美亚洲国产小说| 美女 人体艺术 gogo| 一区在线观看完整版| 人人妻人人澡人人看| 国产成人av教育| 69精品国产乱码久久久| 一本久久中文字幕| av超薄肉色丝袜交足视频| 成人亚洲精品一区在线观看| 日本精品一区二区三区蜜桃| 免费av毛片视频| 亚洲成人精品中文字幕电影| 国产高清视频在线播放一区| 91字幕亚洲| 成年版毛片免费区| 曰老女人黄片| 欧美日韩瑟瑟在线播放| 亚洲国产中文字幕在线视频| 三级毛片av免费| 国产成人一区二区三区免费视频网站| 一级片免费观看大全| 欧美日韩亚洲综合一区二区三区_| 久久人妻福利社区极品人妻图片| 黄色视频,在线免费观看| 日韩欧美免费精品| 中文字幕另类日韩欧美亚洲嫩草| 99热只有精品国产| 亚洲成人国产一区在线观看| 亚洲,欧美精品.| www.999成人在线观看| 日本 av在线| 国产av一区二区精品久久| 丰满的人妻完整版| 亚洲激情在线av| 免费av毛片视频| 美女高潮喷水抽搐中文字幕| 黄片小视频在线播放| 亚洲男人天堂网一区| 午夜福利,免费看| av欧美777| 亚洲国产精品sss在线观看| 欧美av亚洲av综合av国产av| 亚洲成人精品中文字幕电影| 亚洲国产看品久久| 女人高潮潮喷娇喘18禁视频| 黄色丝袜av网址大全| 色尼玛亚洲综合影院| 最近最新中文字幕大全电影3 | 无遮挡黄片免费观看| 国产在线观看jvid| 精品国产一区二区三区四区第35| 男女做爰动态图高潮gif福利片 | 亚洲一码二码三码区别大吗| 叶爱在线成人免费视频播放| 中文字幕另类日韩欧美亚洲嫩草| 两个人视频免费观看高清| 人成视频在线观看免费观看| 国产精品久久久久久人妻精品电影| а√天堂www在线а√下载| 国产精品精品国产色婷婷| 精品久久久久久久毛片微露脸| 日韩欧美免费精品| 老鸭窝网址在线观看| 久久久久国内视频| 亚洲人成伊人成综合网2020| 成年女人毛片免费观看观看9| 午夜免费鲁丝| 久久人人爽av亚洲精品天堂| 国产成人精品久久二区二区免费| 老司机在亚洲福利影院| 好男人电影高清在线观看| 一区二区日韩欧美中文字幕| 亚洲第一青青草原| 少妇裸体淫交视频免费看高清 | 亚洲精品国产色婷婷电影| 久久精品国产综合久久久| 久久香蕉精品热| 日本欧美视频一区| 99国产精品免费福利视频| 少妇裸体淫交视频免费看高清 | 香蕉久久夜色| 亚洲色图av天堂| 国产成人欧美| 久久精品国产清高在天天线| 99国产精品免费福利视频| 免费看美女性在线毛片视频| 一边摸一边抽搐一进一小说| 国产av又大| 精品久久蜜臀av无| 日日摸夜夜添夜夜添小说| 手机成人av网站| ponron亚洲| 婷婷六月久久综合丁香| 久久中文字幕人妻熟女| 嫁个100分男人电影在线观看| 日韩一卡2卡3卡4卡2021年| 黄色丝袜av网址大全| 国产伦一二天堂av在线观看| 国产精品98久久久久久宅男小说| 最好的美女福利视频网| 欧美老熟妇乱子伦牲交| av在线天堂中文字幕| 国产亚洲av嫩草精品影院| 极品人妻少妇av视频| 亚洲第一电影网av| 看黄色毛片网站| 一边摸一边做爽爽视频免费| 人人妻,人人澡人人爽秒播| 亚洲精品美女久久av网站| 亚洲av熟女| 国产极品粉嫩免费观看在线| 亚洲国产精品久久男人天堂| 真人做人爱边吃奶动态| 成人特级黄色片久久久久久久| 搡老岳熟女国产| 精品久久久久久久毛片微露脸| a级毛片在线看网站| 最新在线观看一区二区三区| 亚洲中文日韩欧美视频| 免费久久久久久久精品成人欧美视频| 久久国产亚洲av麻豆专区| 国产亚洲av高清不卡| 国产三级在线视频| 在线天堂中文资源库| 麻豆av在线久日| 99久久精品国产亚洲精品| 国内精品久久久久精免费| 手机成人av网站| 亚洲美女黄片视频| 久久香蕉国产精品| 最近最新免费中文字幕在线| 夜夜夜夜夜久久久久| 午夜免费鲁丝| 一级黄色大片毛片| 精品人妻1区二区| 亚洲男人天堂网一区| 日本 欧美在线| 久久精品国产99精品国产亚洲性色 | 国产极品粉嫩免费观看在线| 高清在线国产一区| 久久精品亚洲精品国产色婷小说| 身体一侧抽搐| 午夜福利一区二区在线看| 欧美日韩乱码在线| 国产精品美女特级片免费视频播放器 | 久久亚洲真实| 国产成人av教育| xxx96com| 亚洲最大成人中文| 神马国产精品三级电影在线观看 | 免费一级毛片在线播放高清视频 | 啪啪无遮挡十八禁网站| 女人被狂操c到高潮| 亚洲精品av麻豆狂野| 亚洲 欧美 日韩 在线 免费| 午夜福利18| 亚洲无线在线观看| 亚洲色图 男人天堂 中文字幕| 国产片内射在线| 国产精品1区2区在线观看.| 国产xxxxx性猛交| 淫秽高清视频在线观看| 黄色片一级片一级黄色片| 日日夜夜操网爽| 动漫黄色视频在线观看| 他把我摸到了高潮在线观看| 在线观看午夜福利视频| 女人被狂操c到高潮| 在线国产一区二区在线| 国产高清视频在线播放一区| 国产精品二区激情视频| 男人操女人黄网站| 国产99久久九九免费精品| 久久久久国产精品人妻aⅴ院| av免费在线观看网站| 女人爽到高潮嗷嗷叫在线视频| 级片在线观看| bbb黄色大片| 国产一区二区激情短视频| 免费av毛片视频| 黄色女人牲交| 日韩视频一区二区在线观看| 精品熟女少妇八av免费久了| 久久青草综合色| 精品电影一区二区在线| 国产一区二区三区综合在线观看| 成人手机av| 久久青草综合色| 悠悠久久av| 国产午夜精品久久久久久| 久9热在线精品视频| 中文字幕久久专区| 神马国产精品三级电影在线观看 | 亚洲最大成人中文| 亚洲美女黄片视频| 久久中文字幕人妻熟女| 亚洲成a人片在线一区二区| 色播在线永久视频| av中文乱码字幕在线| 俄罗斯特黄特色一大片| 99久久国产精品久久久| 99久久99久久久精品蜜桃| 精品人妻1区二区| 久久人妻福利社区极品人妻图片| 夜夜夜夜夜久久久久| 亚洲午夜理论影院| 18禁观看日本| 亚洲性夜色夜夜综合| 女人精品久久久久毛片| 久久精品亚洲精品国产色婷小说| 村上凉子中文字幕在线| 久久人妻熟女aⅴ| 亚洲一区高清亚洲精品| 搞女人的毛片| 久久狼人影院| 午夜福利在线观看吧| 制服人妻中文乱码| 悠悠久久av| 午夜福利免费观看在线| 国产精品秋霞免费鲁丝片| 一区二区三区国产精品乱码| 亚洲全国av大片| 欧美日韩亚洲综合一区二区三区_| 桃色一区二区三区在线观看| 亚洲国产欧美网| 天天躁夜夜躁狠狠躁躁| videosex国产| 欧美中文日本在线观看视频| 欧美国产精品va在线观看不卡| 村上凉子中文字幕在线| 狠狠狠狠99中文字幕| 免费久久久久久久精品成人欧美视频| 侵犯人妻中文字幕一二三四区| 黄片小视频在线播放| 麻豆国产av国片精品| 精品久久久久久久人妻蜜臀av | 深夜精品福利| 欧美成狂野欧美在线观看| 看片在线看免费视频| 亚洲少妇的诱惑av| 久久精品国产亚洲av高清一级| 人妻久久中文字幕网| 国产精品亚洲一级av第二区| 岛国视频午夜一区免费看| 可以在线观看的亚洲视频| 免费观看精品视频网站| 天天躁夜夜躁狠狠躁躁| 91av网站免费观看| 午夜久久久在线观看| 亚洲熟女毛片儿| 亚洲精品国产精品久久久不卡| 操美女的视频在线观看| 精品久久久精品久久久| 国产又爽黄色视频| 亚洲精品久久国产高清桃花| 夜夜躁狠狠躁天天躁| 夜夜看夜夜爽夜夜摸| 久久久久久国产a免费观看| 欧美中文日本在线观看视频| 国产精品98久久久久久宅男小说| 国产精品免费一区二区三区在线| 国产熟女午夜一区二区三区| 亚洲国产日韩欧美精品在线观看 | 丁香六月欧美| 国产野战对白在线观看| 免费观看人在逋| 日本三级黄在线观看| 成人av一区二区三区在线看| 成在线人永久免费视频| 成人亚洲精品av一区二区| 熟妇人妻久久中文字幕3abv| 免费少妇av软件| 亚洲色图 男人天堂 中文字幕| 欧美日韩瑟瑟在线播放| 免费高清视频大片| 女人被狂操c到高潮| 久久婷婷成人综合色麻豆| 满18在线观看网站| 久久国产精品人妻蜜桃| 午夜精品国产一区二区电影| 黄色视频不卡| 国产精华一区二区三区| 中文字幕高清在线视频| 日本黄色视频三级网站网址| 精品乱码久久久久久99久播| 久久这里只有精品19| 亚洲色图av天堂| 美女高潮到喷水免费观看| 亚洲第一欧美日韩一区二区三区| 一边摸一边做爽爽视频免费| 日韩视频一区二区在线观看| 久久影院123| 免费看a级黄色片| 亚洲精品久久国产高清桃花| 免费一级毛片在线播放高清视频 | 91麻豆av在线| 悠悠久久av| 一级,二级,三级黄色视频| 高潮久久久久久久久久久不卡| av在线播放免费不卡| 免费在线观看黄色视频的| 激情在线观看视频在线高清| 在线天堂中文资源库| 自拍欧美九色日韩亚洲蝌蚪91| 精品国产亚洲在线| 亚洲av熟女| 一区二区三区国产精品乱码| 两个人免费观看高清视频| 亚洲成a人片在线一区二区| 亚洲国产高清在线一区二区三 | 国产蜜桃级精品一区二区三区| 国产亚洲欧美98| 欧美一区二区精品小视频在线| 亚洲精华国产精华精| 女性生殖器流出的白浆| 麻豆国产av国片精品| 日韩国内少妇激情av| 99精品欧美一区二区三区四区| 久久精品国产99精品国产亚洲性色 | 国产片内射在线| 久久久国产欧美日韩av| 视频在线观看一区二区三区| 黄色毛片三级朝国网站| 久久婷婷人人爽人人干人人爱 | 精品电影一区二区在线| 狠狠狠狠99中文字幕| 看免费av毛片| 国产一区在线观看成人免费| 18禁黄网站禁片午夜丰满| 一二三四在线观看免费中文在| 国产av精品麻豆| 人人澡人人妻人| 久久久久亚洲av毛片大全| 久久久久久久久久久久大奶| 国产成人影院久久av| 中出人妻视频一区二区| 亚洲熟女毛片儿| 俄罗斯特黄特色一大片| 亚洲第一电影网av| 免费看a级黄色片| 给我免费播放毛片高清在线观看| 精品久久久精品久久久| 亚洲av美国av| 中出人妻视频一区二区| 99久久久亚洲精品蜜臀av| 夜夜爽天天搞| 在线国产一区二区在线| 在线永久观看黄色视频| 老司机午夜十八禁免费视频| 亚洲五月色婷婷综合| 日日干狠狠操夜夜爽| 可以在线观看毛片的网站| 两个人免费观看高清视频| 三级毛片av免费| 亚洲成人免费电影在线观看| 欧美日韩乱码在线| 精品国产乱子伦一区二区三区| 中文字幕人成人乱码亚洲影| 日韩有码中文字幕| 亚洲av五月六月丁香网| 国产成人精品在线电影| 精品第一国产精品| 婷婷精品国产亚洲av在线| 青草久久国产| 成人国产综合亚洲| 成在线人永久免费视频| 国产aⅴ精品一区二区三区波| 国产野战对白在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲人成77777在线视频| 国产成人av激情在线播放| 亚洲激情在线av| 成人亚洲精品av一区二区| 99国产极品粉嫩在线观看| 午夜精品久久久久久毛片777| 国产精品免费视频内射| 国产成人精品久久二区二区91| 9191精品国产免费久久| 99国产精品免费福利视频| 成人18禁在线播放| 日本五十路高清| av福利片在线| 久久久久久免费高清国产稀缺| 一级片免费观看大全| 国产伦一二天堂av在线观看| 日韩三级视频一区二区三区| 久久热在线av| 午夜福利一区二区在线看| 精品日产1卡2卡| 久久久久亚洲av毛片大全| 日本一区二区免费在线视频| 亚洲成人免费电影在线观看| 久久亚洲真实| 成人国产综合亚洲| 一夜夜www| 欧美乱码精品一区二区三区| 女性生殖器流出的白浆| 亚洲欧美日韩高清在线视频| 一进一出好大好爽视频| 美女扒开内裤让男人捅视频| 啦啦啦免费观看视频1| 欧美中文日本在线观看视频| 亚洲av成人不卡在线观看播放网| 午夜福利视频1000在线观看 | av在线天堂中文字幕| 国产1区2区3区精品| 亚洲av成人av| 黄片播放在线免费| 久久天堂一区二区三区四区| 多毛熟女@视频| 巨乳人妻的诱惑在线观看| 国产97色在线日韩免费| 欧美+亚洲+日韩+国产| 男人的好看免费观看在线视频 | 国产熟女午夜一区二区三区| 亚洲三区欧美一区| 夜夜爽天天搞| 免费无遮挡裸体视频| 亚洲av熟女| 又黄又爽又免费观看的视频| 满18在线观看网站| 亚洲专区字幕在线| 久久中文字幕一级| 男男h啪啪无遮挡| 国产三级黄色录像| 亚洲色图综合在线观看| 国产精品九九99| 一区二区三区国产精品乱码| 琪琪午夜伦伦电影理论片6080| 波多野结衣巨乳人妻| 午夜福利免费观看在线| 国产亚洲精品av在线| 91成年电影在线观看| 91九色精品人成在线观看| 中文字幕另类日韩欧美亚洲嫩草| 熟妇人妻久久中文字幕3abv| 两个人视频免费观看高清| 欧美激情久久久久久爽电影 | 国产亚洲精品综合一区在线观看 | 亚洲一区二区三区色噜噜| 男人操女人黄网站| 午夜福利高清视频| 国产成人av教育| 亚洲少妇的诱惑av| 狂野欧美激情性xxxx| 在线观看一区二区三区| 国产片内射在线| 精品乱码久久久久久99久播| 50天的宝宝边吃奶边哭怎么回事| 色综合站精品国产| 美国免费a级毛片| 欧美成人性av电影在线观看| 欧美日韩精品网址| 高清毛片免费观看视频网站| 一区二区三区激情视频| 99香蕉大伊视频| 天天躁夜夜躁狠狠躁躁| 亚洲人成伊人成综合网2020| 午夜影院日韩av| 免费久久久久久久精品成人欧美视频| av视频免费观看在线观看| 国产精品一区二区在线不卡| 国产精品一区二区三区四区久久 | 狂野欧美激情性xxxx| 国产成人系列免费观看| 亚洲avbb在线观看| 国产91精品成人一区二区三区| 黄色a级毛片大全视频| 美女高潮到喷水免费观看| 久久 成人 亚洲| 免费观看人在逋| 777久久人妻少妇嫩草av网站| 亚洲人成电影观看| 麻豆久久精品国产亚洲av| 69精品国产乱码久久久| 黑人巨大精品欧美一区二区蜜桃| av片东京热男人的天堂| 巨乳人妻的诱惑在线观看| 日本三级黄在线观看| 深夜精品福利| 亚洲av片天天在线观看| 亚洲性夜色夜夜综合| 色婷婷久久久亚洲欧美| 高清毛片免费观看视频网站| 国产欧美日韩一区二区三区在线| 非洲黑人性xxxx精品又粗又长| 在线视频色国产色| 伦理电影免费视频| 中文字幕人妻丝袜一区二区| 亚洲av片天天在线观看| 国产高清有码在线观看视频 | 国产私拍福利视频在线观看| 两个人看的免费小视频| 人人妻人人爽人人添夜夜欢视频| 日本一区二区免费在线视频| 色哟哟哟哟哟哟|