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

    南京北郊VOCs對臭氧和二次有機(jī)氣溶膠潛在貢獻(xiàn)的研究

    2015-11-23 06:37:18安俊琳南京信息工程大學(xué)氣象災(zāi)害預(yù)報(bào)預(yù)警與評估協(xié)同創(chuàng)新中心中國氣象局氣溶膠降水重點(diǎn)開放實(shí)驗(yàn)室江蘇南京210044
    中國環(huán)境科學(xué) 2015年4期
    關(guān)鍵詞:北郊潛勢芳香烴

    林 旭,朱 彬,安俊琳,楊 輝 (南京信息工程大學(xué),氣象災(zāi)害預(yù)報(bào)預(yù)警與評估協(xié)同創(chuàng)新中心,中國氣象局氣溶膠-云-降水重點(diǎn)開放實(shí)驗(yàn)室,江蘇 南京 210044)

    南京北郊VOCs對臭氧和二次有機(jī)氣溶膠潛在貢獻(xiàn)的研究

    林 旭,朱 彬*,安俊琳,楊 輝 (南京信息工程大學(xué),氣象災(zāi)害預(yù)報(bào)預(yù)警與評估協(xié)同創(chuàng)新中心,中國氣象局氣溶膠-云-降水重點(diǎn)開放實(shí)驗(yàn)室,江蘇 南京 210044)

    2013年3月1日~2014年2月28日采用GC5000在線氣相色譜儀對南京北郊大氣環(huán)境中的揮發(fā)性有機(jī)物(VOCs)進(jìn)行了為期一年的連續(xù)監(jiān)測,分析了VOCs的組成特征及季節(jié)、日變化規(guī)律,并結(jié)合PMF受體模型,采用最大增量反應(yīng)活性(MIR)系數(shù)及氣溶膠生成系數(shù)(FAC)分析了VOCs及其各來源的O3和SOA生成潛勢.結(jié)果顯示,南京北郊大氣環(huán)境中TVOCs小時(shí)平均體積分?jǐn)?shù)為45.63×10-9,TVOCs及各組分濃度呈現(xiàn)秋冬季高、夏季低的季節(jié)變化特征和雙峰結(jié)構(gòu)的日變化規(guī)律.SOA總生成量約為2.07μg/m3,芳香烴對其貢獻(xiàn)率最大,占95.93%,其中的苯系物是生成SOA的優(yōu)勢物種;烯烴對OFP的貢獻(xiàn)最大,接近65%;烷烴雖為VOCs的優(yōu)勢組分,卻并非OFP和SOA的主要貢獻(xiàn)者.不同季節(jié)VOCs的受體模型解析結(jié)果顯示,工業(yè)排放和汽車尾氣是南京北郊最主要的VOCs來源.富含苯系物的VOCs來源對SOA的貢獻(xiàn)最大,對OFP貢獻(xiàn)最大的則為富含乙、丙烯及異戊二烯的VOCs來源;春、秋、冬三季汽車尾氣及工業(yè)排放(包括石化工業(yè))二源對大氣中的VOCs濃度、SOA的貢獻(xiàn)及OFP的貢獻(xiàn)影響最為顯著,而夏季溶劑使用及植物源對SOA及OFP的貢獻(xiàn)不容忽視.

    揮發(fā)性有機(jī)物;SOA;OFP;來源;貢獻(xiàn)

    揮發(fā)性有機(jī)物(VOCs)是對流層臭氧(O3)和二次有機(jī)氣溶膠(SOA)的重要前體物[1-2],對城市O3和灰霾等復(fù)合型大氣污染的形成至關(guān)重要[3-5],同時(shí),很多VOCs組分為有毒有害物質(zhì),嚴(yán)重危害人體健康[6].因此,由VOCs引起的大氣環(huán)境污染問題已受到世界上許多國家的廣泛關(guān)注,而其對大氣復(fù)合污染的貢獻(xiàn)亦成為其中的一大研究熱點(diǎn)[7-10].

    國內(nèi)外就VOCs對O3的潛在貢獻(xiàn)方面已展開了大量的研究,Bowman等[11]和Carter等[8]研究表明大氣中VOCs成分復(fù)雜多變且彼此之間臭氧生成潛勢差異極大.吳方等[12]對北京奧運(yùn)期間VOCs濃度變化、臭氧產(chǎn)生潛勢及來源進(jìn)行了研究,結(jié)果顯示北京地區(qū)大氣中O3產(chǎn)生潛力最大的是芳香烴類化合物占47%,其次為烯烴占40%,烷烴最低占13%;臭氧產(chǎn)生潛力最高的物種為1,2,4-三甲苯,其次為間二甲苯,對二甲苯,甲苯,1,3,5-三甲苯和1-丁烯.陳長虹等[13]對上海市城區(qū)VOCs的年變化特征和關(guān)鍵活性組分的研究表明芳香烴對OFP的貢獻(xiàn)率最大為62.75%,其次為烯烴21.70%,VOCs中的關(guān)鍵活性組分是二甲苯、甲苯、乙苯、乙烯、丙烯、反-2-丁烯及異戊二烯等.朱少峰等[14]對深圳大氣VOCs濃度變化特征與化學(xué)反應(yīng)活性進(jìn)行了探討,結(jié)果顯示芳香烴物種是臭氧生成最主要的貢獻(xiàn)者,其次為烯烴和烷烴,且甲苯和間,對-二甲苯是對臭氧生成貢獻(xiàn)最大的2個(gè)物種,占總VOCs的1/3以上.安俊琳等[15]和李用宇等[16]分別就南京北郊VOCs體積分?jǐn)?shù)和光化學(xué)特征進(jìn)行了深入探討,結(jié)果均顯示烯烴對OFP的貢獻(xiàn)最大,芳香烴次之,烷烴對OFP貢獻(xiàn)最小.

    二次有機(jī)氣溶膠(SOA)是細(xì)顆粒物的重要組成部分[17],平均占PM2.5有機(jī)組分質(zhì)量的20%~50%[18],它們是天然源和人為源排放的VOCs(揮發(fā)性有機(jī)物)或SVOC(半揮發(fā)性有機(jī)物)經(jīng)大氣氧化和氣/粒分配等過程而生成的懸浮于大氣中的固體或液體微粒,能影響人體健康、降低能見度及影響氣候變化[19-20].近年來,國內(nèi)外學(xué)者就SOA生成機(jī)制已開展了大量的煙霧箱實(shí)驗(yàn)研究[9-10,21-23],其中使用較為廣泛的是Grosjean等[9-10]在綜合大量煙霧箱實(shí)驗(yàn)數(shù)據(jù)和大氣化學(xué)動(dòng)力學(xué)數(shù)據(jù)的基礎(chǔ)上提出的一套氣溶膠生成系數(shù)(fractional aerosol coefficient,F(xiàn)AC).Barthelmie等[24]、Kourtidis等[25]及Dechapanya等[26]分別在英國哥倫比亞、希臘雅典和美國雷克薩斯州的休斯頓地區(qū)結(jié)合VOCs觀測數(shù)據(jù)及FAC對SOA的生成潛勢做了深入的探討,結(jié)果均顯示芳香烴對SOA的生成貢獻(xiàn)最大.國內(nèi)的呂子峰等[27]和王倩等[28]利用FAC系數(shù)分別對北京和上海的SOA生成潛勢進(jìn)行了探討研究,結(jié)果亦顯示芳香烴是生成SOA的優(yōu)勢組分.王扶潘等[29]利用FAC系數(shù)法對深圳大氣中各芳香烴和異戊二烯的SOA生成潛勢進(jìn)行了研究,結(jié)果顯示甲苯對SOA生成貢獻(xiàn)最大.

    目前,國內(nèi)對SOA生成潛勢的相關(guān)報(bào)道總體比較匱乏,而針對VOCs的氧化對O3和SOA的潛在貢獻(xiàn)的綜合對比分析更是鮮有研究.為此,本文基于在線的GC5000氣相色譜儀對南京北郊大氣中的VOCs進(jìn)行了連續(xù)一年的觀測,并結(jié)合最大增量反應(yīng)活性(MIR)系數(shù)[8]及氣溶膠生成系數(shù)(FAC)[9-10]分析了O3生成潛勢和SOA生成潛勢,并采用PMF受體模型對大氣中的VOCs進(jìn)行了來源解析,在此基礎(chǔ)上進(jìn)一步深入對比分析了各VOCs來源的O3及SOA生成潛勢,以期待為南京市制定相關(guān)政策綜合控制大氣復(fù)合污染提供一定的科學(xué)依據(jù).

    1 材料與方法

    1.1 采樣地點(diǎn)與采樣時(shí)間

    本次研究觀測地點(diǎn)位于南京市浦口區(qū)南京信息工程大學(xué)氣象樓樓頂(32o12'N,118o42'E,海拔高度:62m).采樣點(diǎn)東側(cè)毗鄰主干道寧六路,車流量較大;西南邊為龍王山風(fēng)景區(qū);測點(diǎn)以東3~6km處是包括石化工業(yè)、鋼鐵廠、化工廠和熱電廠等大型企業(yè)的工業(yè)園區(qū).采樣時(shí)間為2013 年03月01日~2014年02月28日.

    1.2 儀器與分析方法

    VOCs觀測采用由德國AMA公司提供的集自動(dòng)采樣、富集和分析功能于一體的GC5000在線氣相色譜儀.該系統(tǒng)主要包括GC5000VOC(單級富集)和GC5000BTX(兩級富集)兩臺自動(dòng)分析色譜儀,可分別測量C2~C6的低沸點(diǎn)物種和C4~C12的高沸點(diǎn)物種,檢測器均為氫火焰離子化檢測器(FID).

    測量原理:環(huán)境空氣樣品通過采樣經(jīng)干燥后直接進(jìn)入分析系統(tǒng),有機(jī)物在富集管的特殊材料中吸附解析后,進(jìn)入毛細(xì)管柱進(jìn)行分離,而后由FID檢測有機(jī)物含量.儀器可檢測C2~C12共計(jì)56種揮發(fā)性有機(jī)物,包括29種烷烴、10種烯烴、16種芳香烴和1種炔烴.系統(tǒng)時(shí)間分辨率為1h.同時(shí)為保證數(shù)據(jù)的有效性和可靠性,觀測期間利用DIM200校準(zhǔn)模塊(稀釋100倍)進(jìn)行4周左右一次的校正,其中校正氣體采用美國環(huán)保署認(rèn)可的Linde Gas North America LLC提供的混合標(biāo)氣. 1.3 VOCs的SOA及O3生成潛勢研究方法

    本研究基于Grosjean的煙霧箱實(shí)驗(yàn)[9],采用氣溶膠生成系數(shù)(FAC)[9-10]來估算環(huán)境大氣中VOCs的SOA生成潛勢:

    式中:SOAp是SOA生成潛勢,μg/m3;VOCs0是排放源排出的初始濃度,μg/m3;FAC是SOA的生成系數(shù).考慮到受體點(diǎn)測得的VOCs往往是經(jīng)過氧化后的濃度VOCst,它與排放源排出的初始濃度VOCs0之間的關(guān)系可通過下式表示:

    式中:FVOCr為VOCs物種參與反應(yīng)的分?jǐn)?shù)(%). Grosjean[9]假設(shè)SOA只在白天生成( 08: 00~17: 00),且VOCs只與OH發(fā)生反應(yīng)生成SOA,公式(1)和(2)中用到的FAC和FVOCr由煙霧箱實(shí)驗(yàn)獲得[9-10].

    采用最大增量反應(yīng)活性(MIR)系數(shù)分析O3生成潛勢(OFP),定義公式如下:

    式中:[VOC]i為實(shí)際觀測的某種VOC大氣環(huán)境濃度;MIR為VOCs物種的最大增量反應(yīng)活性,該系數(shù)采用Carter[8]實(shí)驗(yàn)結(jié)果,計(jì)算結(jié)果見表1.

    1.4 PMF受體模型介紹

    PMF(正矩陣因子分析)由芬蘭赫爾辛基大學(xué)的Dr.Paatero在20世紀(jì)90年代中期開發(fā).隨著PMF2.0,PMF3.0等版本的發(fā)布,該模式已廣泛地應(yīng)用于空氣質(zhì)量中的源解析研究.其原理為:假設(shè)X為n×m矩陣,n為樣品數(shù),m為化學(xué)成分的數(shù)目,那么X可以分解為X=GF+E,其中G為n×p的矩陣,F(xiàn)為p×m的矩陣,p為主要污染源的數(shù)目,E為殘數(shù)矩陣[30].定義

    式中:Xij為第i個(gè)樣本中第j種物質(zhì)的濃度;gik為第k個(gè)因子對第i個(gè)樣本的貢獻(xiàn);fkj為第j種物質(zhì)在第k個(gè)因子中的百分?jǐn)?shù); eij為第i個(gè)樣本中第j種物質(zhì)的殘余;sij為第i個(gè)樣本中第j種物質(zhì)的不確定值估計(jì).

    PMF模式解析過程即要求Q值趨于最小,在gik≥0, fkj≥0的約束條件下,通過迭代最小化算法對Q求解,同時(shí)確定G(污染源貢獻(xiàn))和F(源成分譜)[31].PMF受體模型運(yùn)行前需要為其提供兩類數(shù)據(jù)[32],一類是濃度數(shù)據(jù),一類是不確定值數(shù)據(jù),相應(yīng)的數(shù)據(jù)處理方法見文獻(xiàn)[33-34].

    2 結(jié)果與討論

    2.1 南京北郊VOCs組成及變化特征

    本次實(shí)驗(yàn)共檢測出56種VOCs,包括29種烷烴,10種烯烴,16種芳香烴和1種炔烴(乙炔).表1給出了觀測期間各VOCs平均體積分?jǐn)?shù)的統(tǒng)計(jì)結(jié)果,南京北郊全年TVOCs(總揮發(fā)性有機(jī)物)小時(shí)平均體積分?jǐn)?shù)為45.63×10-9,其中烷烴為21.02×10-9,烯烴為10.01×10-9,芳香烴為10.02× 10-9,乙炔為4.58×10-9.它們分別占TVOCs的46.07%、21.94%、21.96%和10.03%,可見烷烴是南京北郊大氣環(huán)境中含量最為豐富的組分.對VOCs體積分?jǐn)?shù)貢獻(xiàn)率排名前10的物種依次為:乙烷13.17%、乙烯12.29%、乙炔10.04%、丙烷7.54%、苯6.38%、甲苯5.39%、丙烯5.19%、正丁烷4.60%、異丁烷4.12%及乙基苯4.08%,累計(jì)占總VOCs的72.80%,可見南京北郊大氣環(huán)境中VOCs的優(yōu)勢物種為C2~C4的烷、烯烴,乙炔及苯系物,這與李用宇等[16]的研究結(jié)果一致.

    表1 南京北郊VOCs平均體積分?jǐn)?shù)、二次有機(jī)氣溶膠及臭氧生成潛勢Table 1 Average mixing ratio of the VOCs,the potential formation of SOA and O3in the northern suburb of Nanjing

    續(xù)表1

    圖1(a)顯示觀測期間南京北郊大氣中TVOCs體積分?jǐn)?shù)呈現(xiàn)出明顯的季節(jié)變化特征:秋季>冬季>春季>夏季.烷烴和烯烴與TVOCs呈現(xiàn)出一致的季節(jié)變化特征,而芳香烴和乙炔則略有不同,呈現(xiàn)出:冬季>秋季>春季>夏季的變化特征.一般來說,VOCs濃度的季節(jié)性變化主要受3個(gè)因素的影響[35]:(1)VOCs源強(qiáng)的改變;(2)OH自由基的季節(jié)性豐度不同,從而導(dǎo)致光化學(xué)反應(yīng)程度不同;(3)大氣混合狀態(tài)的變化.南京北郊TVOCs及其各組分秋冬季體積分?jǐn)?shù)較大主要是由于相對于夏季來說秋冬季溫度較低,大氣光化學(xué)反應(yīng)不太活躍,不利于VOCs的去除,且大氣層結(jié)較穩(wěn)定,夜間常出現(xiàn)逆溫現(xiàn)象,不利于污染物的擴(kuò)散稀釋.安俊琳等[15]研究顯示南京北郊VOCs的體積分?jǐn)?shù)呈現(xiàn)夏季高,冬季低的季節(jié)變化規(guī)律,與本文結(jié)論存在差異,這主要是由于本文所用數(shù)據(jù)的采樣時(shí)間為2013年3月1日~2014年2月28日,第二屆亞洲青年運(yùn)動(dòng)會于2013年8月16日~8 月24日在南京舉辦,為此南京市政府出臺了一系列環(huán)境質(zhì)量保障臨時(shí)管控措施,并取得不錯(cuò)的成效,因而本次實(shí)驗(yàn)期間夏季VOCs濃度相對偏低.

    污染物濃度的日變化主要受2個(gè)原因的影響,一是污染源強(qiáng)度的變化,二是擴(kuò)散條件(如氣溫、氣壓、風(fēng)速)的改變[36].圖1(b)顯示觀測期間南京北郊TVOCs、烷烴、烯烴、芳香烴及乙炔的日變化曲線呈現(xiàn)出較為一致的雙峰結(jié)構(gòu),峰值主要出現(xiàn)在早上8:00和晚上的20:00左右,與交通高峰時(shí)間基本吻合,極小值出現(xiàn)在14:00.早晨隨著上班高峰的來臨,交通、工業(yè)等人類活動(dòng)逐漸頻繁,使得VOCs體積分?jǐn)?shù)迅速增加達(dá)第1個(gè)峰值,之后,太陽輻射及湍流運(yùn)動(dòng)不斷加強(qiáng),大氣光化學(xué)反應(yīng)程度加劇,VOCs體積分?jǐn)?shù)逐漸下降于14:00出現(xiàn)極小值;傍晚太陽輻射減弱,光化學(xué)反應(yīng)逐步停止,并伴隨交通晚高峰的到來,VOCs濃度水平開始逐步回升并達(dá)到第2個(gè)峰值.此外,VOCs及各組分體積分?jǐn)?shù)均表現(xiàn)出夜大于晝的變化特征,這主要是因?yàn)橐归g沒有光化學(xué)反應(yīng)的消耗及較低的大氣邊界層高度不利于污染物的擴(kuò)散導(dǎo)致的.

    圖1 南京北郊VOCs季節(jié)及日變化特征Fig.1 Seasonal variation and diurnal variation of VOCs in the northern suburb of Nanjing

    2.2 南京北郊VOCs的二次有機(jī)氣溶膠及臭氧生成潛勢

    基于Grosjean的煙霧箱實(shí)驗(yàn)[9],采用氣溶膠生成系數(shù)(FAC)[9-10]估算南京北郊大氣中VOCs的SOA生成潛勢,并結(jié)合最大增量反應(yīng)活性(MIR)系數(shù)[8]分析臭氧生成潛勢(OFP).雖然利用FAC系數(shù)法估算SOA的生成潛勢時(shí)只考慮了VOCs與OH自由基發(fā)生的反應(yīng),而并沒有考慮到NO3自由基和O3,并且本研究也只能測得大氣環(huán)境中的56種VOCs物種,還有其他對SOA有貢獻(xiàn)的VOCs和SVOC無法測得,因而估算的結(jié)果偏低,但是仍能給出SOA生成的大致數(shù)量級,并能指明各SOA前體物的相對貢獻(xiàn),對于識別其關(guān)鍵前體物很有幫助[26-28].表1顯示對SOA有貢獻(xiàn)的VOCs物種共24個(gè),其中烷烴11個(gè),芳香烴13個(gè).SOA的總生成量約為2.07μg/m3,其中烷烴和芳香烴的SOA生成貢獻(xiàn)分別占4.07%和95.93%,與王倩等[28]的研究結(jié)果較為接近.結(jié)合VOCs全年平均體積分?jǐn)?shù)分析可發(fā)現(xiàn)占TVOCs體積分?jǐn)?shù)近一半的烷烴對SOA生成的貢獻(xiàn)率不到5%,而僅占TVOCs體積分?jǐn)?shù)20%左右的芳香烴對SOA的貢獻(xiàn)率則超過95%,且不同季節(jié)(表2)中,芳香烴均為生成SOA的關(guān)鍵組分,其SOA的生成量冬季最高,夏季最低,與其體積分?jǐn)?shù)的季節(jié)變化規(guī)律基本一致.進(jìn)一步分析發(fā)現(xiàn)對SOA生成貢獻(xiàn)率最大的5種VOCs物種分別為甲苯27.28%,乙基苯26.12%,苯9.54%,間、對-二甲苯9.36%和鄰二甲苯5.60%,累積占總貢獻(xiàn)的77.90%.可見芳香烴中苯系物是生成SOA的絕對優(yōu)勢物種.

    圖2 南京北郊環(huán)境大氣中各VOCs組分對總VOCs濃度、SOA及OFP的相對貢獻(xiàn)率Fig.2 Relative contribution of the components of VOCs to the total VOCs concentration,SOA and OFP in the northern suburb of Nanjing

    對VOCs物種進(jìn)行OFP分析后發(fā)現(xiàn),體積分?jǐn)?shù)僅占20%左右的烯烴對OFP的貢獻(xiàn)接近65%,是大氣化學(xué)反應(yīng)活性最強(qiáng)的組分,這與安俊琳等[15]和李用宇等[16]的研究結(jié)果一致.芳香烴對OFP的貢獻(xiàn)僅次于烯烴,其占TVOCs體積分?jǐn)?shù)的百分比和OFP貢獻(xiàn)率均在20%左右;而占體積分?jǐn)?shù)比重最大的烷烴對OFP的貢獻(xiàn)卻最小,僅為13%左右,對SOA生成的貢獻(xiàn)率亦不到5%.不同季節(jié)中(表2)烯烴對OFP的貢獻(xiàn)均最為顯著,芳香烴次之,烷烴最小.做為主要貢獻(xiàn)者,烯烴對OFP的的影響秋季最大,夏季最小,與其體積分?jǐn)?shù)的季節(jié)變化呈現(xiàn)一致的規(guī)律.可見烷烴雖為南京北郊大氣環(huán)境中VOCs的優(yōu)勢組分,卻并非OFP 和SOA的主要貢獻(xiàn)者.

    表2 南京北郊不同季節(jié)各VOCs組分的二次有機(jī)氣溶膠及臭氧生成潛勢Table 2 The potential formation of VOCs components in different seasons to SOA and OFP in the northern suburb of Nanjing

    圖2與表1表明芳香烴對SOA貢獻(xiàn)最大,其貢獻(xiàn)率超過95%;其次為烷烴,貢獻(xiàn)率為4.07%. 對SOA貢獻(xiàn)最大的10種物質(zhì)及其貢獻(xiàn)率分別為:甲苯27.28%、乙基苯26.12%、苯9.54%、間、對-二甲苯9.36%、鄰二甲苯5.60%、鄰乙基甲苯5.1%、1,2,4-三甲基苯4.2%、1,3,5-三甲基苯2.52%、對乙基甲苯2.19%及間乙基甲苯1.45%,累計(jì)貢獻(xiàn)率為93.36%,可見苯系物是生成SOA的絕對優(yōu)勢物種.烯烴是對OFP貢獻(xiàn)最大的組分,貢獻(xiàn)率達(dá)64.46%;其次為芳香烴,貢獻(xiàn)率為20.62%.對OFP貢獻(xiàn)排名前10的物種及各自貢獻(xiàn)率為:乙烯32.65%、丙烯17.52%、異戊二烯8.09%、甲苯5.22%、乙基苯3.95%、間、對-二甲苯3.61%、1-丁烯3.0%、1,2,4-三甲基苯2.98%、鄰二甲苯1.99%及乙炔1.80%,可見C2~C5的烯烴、苯系物及乙炔是生成O3的優(yōu)勢物種.總體而言,雖然在南京北郊大氣環(huán)境中對總VOCs體積分?jǐn)?shù)、SOA及OFP這三者貢獻(xiàn)最大的VOCs組分各不相同,依次為:烷烴、芳香烴及烯烴,但苯系物對三者均有很大貢獻(xiàn),尤其是對SOA貢獻(xiàn)最為顯著,是生成SOA的絕對優(yōu)勢物種,且苯系物多為有毒有害物質(zhì),嚴(yán)重危害人體健康[37],應(yīng)引起高度關(guān)注.

    2.3 南京北郊VOCs來源及其對SOA和O3的潛在貢獻(xiàn)

    為了進(jìn)一步深入探討VOCs對SOA及OFP的貢獻(xiàn)情況,利用PMF受體模型對南京北郊大氣環(huán)境中VOCs分季節(jié)進(jìn)行來源解析,并利用上述的氣溶膠生成系數(shù)法及MIR系數(shù)法就四季VOCs各來源對SOA及OFP的貢獻(xiàn)進(jìn)行進(jìn)一步的分析研究.

    經(jīng)PMF受體模型逐步逼近后,確定了南京北郊春季VOCs的主要來源有5個(gè)(圖3):燃料揮發(fā)源占18.42%、汽車尾氣排放源占27.22%、工業(yè)排放源占32.83%、溶劑使用源占10.71%、植物排放源占10.82%.夏季南京北郊大氣環(huán)境中5種VOCs來源的貢獻(xiàn)比分別為:植物排放源占16.15%、工業(yè)排放源占31.64%、汽車尾氣排放源占21.43%、溶劑使用源占11.74%、化學(xué)工業(yè)源占19.04%.秋季VOCs5種來源的貢獻(xiàn)百分比分別為:燃料揮發(fā)源占20.69%、工業(yè)排放源占25.14%、汽車尾氣排放源占24.14%、溶劑使用源占8.89%、石油化工源占21.15%.冬季南京北郊VOCs共有4個(gè)來源,其貢獻(xiàn)比分布為:燃料揮發(fā)源占19.81%、工業(yè)排放源占41.83%、汽車尾氣排放源占21.35%、溶劑使用源占17.02%.可見工業(yè)排放源和汽車尾氣排放源是南京北郊最主要的VOCs來源,與夏麗[38]的研究結(jié)果較為一致.

    圖4給出了采樣點(diǎn)四季VOCs各來源的SOA貢獻(xiàn)情況,除夏季以外,其余三季對SOA貢獻(xiàn)最大的2種VOCs來源呈現(xiàn)出一致的規(guī)律:汽車尾氣對SOA的貢獻(xiàn)最大,其次為溶劑使用.葉偉紅等[39]認(rèn)為汽車尾氣中苯系物的含量豐富,而溶劑中亦含有大量的BTEX[40],此二源對SOA的貢獻(xiàn)最大表明苯系物是生成SOA的優(yōu)勢物種,這與表1和圖2的結(jié)論一致.而夏季汽車尾氣對SOA的貢獻(xiàn)明顯低于其他季節(jié),這主要是因?yàn)橐环矫?013年夏季亞青會在南京舉辦,南京市政府出臺一系列環(huán)境質(zhì)量保障臨時(shí)管控措施,如延長機(jī)動(dòng)車限行時(shí)間、增加限行車型等,大大降低了汽車尾氣排放源強(qiáng)度;另一方面南京夏季盛行偏南風(fēng),而主要的機(jī)動(dòng)車污染源為采樣點(diǎn)東邊的寧六路,所以夏季汽車尾氣對SOA貢獻(xiàn)明顯降低,同時(shí)溶劑使用源因高溫而具備高揮發(fā)性,對SOA的貢獻(xiàn)較為顯著.

    圖3 南京北郊春、夏、秋、冬四季VOCs來源組成百分比Fig.3 Contribution of various emission sources to VOCs in spring,summer,autumn and winter in the northern suburb of Nanjing

    圖5給出了采樣點(diǎn)四季各VOCs來源的OFP貢獻(xiàn)情況,南京北郊四個(gè)季節(jié)中各VOCs來源對OFP的貢獻(xiàn)沒有較為一致的分布規(guī)律.春、夏、秋、冬四季中對OFP貢獻(xiàn)最大的源分別為:汽車尾氣排放源,植物排放源,石油化工源和工業(yè)源.其中秋季石油化工對OFP的高貢獻(xiàn)率主要與亞青會后工廠全面恢復(fù)生產(chǎn)及秋季盛行東北風(fēng),主要的石化企業(yè)污染源位于采樣點(diǎn)的上風(fēng)向有關(guān).結(jié)合表1和圖2發(fā)現(xiàn)南京北郊VOCs中對OFP貢獻(xiàn)最大的組分為烯烴,其中乙烯、丙烯和異戊二烯是其中最關(guān)鍵的活性物種,其對OFP的累積貢獻(xiàn)率近60%.而陸思華等[41]研究表明乙烯主要存在于石油化工及汽車尾氣中,丙烯主要分布于石油液化氣及汽車尾氣中,而異戊二烯是植物排放的示蹤物.此外乙烯和丙烯對工業(yè)排放源的貢獻(xiàn)也較大.由此表明,雖然南京北郊四季各VOCs來源的OFP生成貢獻(xiàn)沒有較為一致的規(guī)律,但是它們從根本上反應(yīng)出共同的本質(zhì)特征即富含烯烴尤其是乙烯、丙烯及異戊二烯的來源對OFP貢獻(xiàn)最大.對比圖4和圖5可發(fā)現(xiàn),春、夏、秋、冬四季對SOA和OFP的貢獻(xiàn)最大的源分別為:汽車尾氣排放源和汽車尾氣排放源、溶劑使用源和植物排放源、汽車尾氣排放源和石油化工源、汽車尾氣排放源和工業(yè)排放源.此外,圖3亦顯示工業(yè)排放和汽車尾氣排放是南京北郊VOCs最主要的2種排放源.所以,在春、秋、冬三季從源頭方面控制大氣環(huán)境中的VOCs濃度、VOCs對SOA的貢獻(xiàn)及VOCs對OFP的貢獻(xiàn)這三者目標(biāo)是一致的即控制汽車尾氣排放及工業(yè)排放(包括石化工業(yè)).夏季除關(guān)注上述二源外,溶劑使用及植物源因高溫而具備高揮發(fā)性也應(yīng)引起高度關(guān)注.

    圖4 南京北郊四季VOCs來源對SOA的潛在貢獻(xiàn)Fig.4 The potential contribution of the sources of VOCs in different seasons to SOA in the northern suburb of Nanjing

    圖5 南京北郊四季VOCs來源對OFP的潛在貢獻(xiàn)Fig.5 The potential contribution of the sources of VOCs in different seasons to OFP in the northern suburb of Nanjing

    3 結(jié)論

    3.1 觀測期間,南京北郊大氣環(huán)境中TVOCs小時(shí)平均體積分?jǐn)?shù)為45.63×10-9,其中烷烴含量最高,占TVOCs的46.07%;其次為芳香烴、烯烴和乙炔,分別占TVOCs的21.96%、21.94%和10.03 %.TVOCs、烷烴和烯烴的濃度呈現(xiàn)出一致的季節(jié)變化特征:秋季>冬季>春季>夏季;而芳香烴和乙炔則略有不同,呈現(xiàn)出:冬季>秋季>春季>夏季的季節(jié)變化特征.TVOCs及其各組分的日變化規(guī)律呈現(xiàn)出較為一致的雙峰結(jié)構(gòu),峰值出現(xiàn)在早8:00和晚20:00左右,谷值出現(xiàn)在14:00,夜間濃度大于白天.

    3.2 測得的56種VOCs中對SOA有貢獻(xiàn)的共

    24種,其中烷烴11種,芳香烴13種,對生成SOA的相對貢獻(xiàn)率分別為4.07%和95.93%,SOA的總生成量約為2.07μg/m3,對SOA生成貢獻(xiàn)率最大的5種VOCs物種分別為甲苯27.28%,乙基苯26.12 %,苯9.54%,間、對-二甲苯9.36%和鄰二甲苯5.60%,芳香烴中苯系物是生成SOA的絕對優(yōu)勢物種.烯烴對OFP的貢獻(xiàn)最大,接近65%,其次為芳香烴,而烷烴雖為大氣中VOCs的優(yōu)勢組分,卻并非OFP和SOA的主要貢獻(xiàn)者.

    3.3 雖然南京北郊大氣環(huán)境中對總VOCs體積分?jǐn)?shù)、SOA及OFP這三者貢獻(xiàn)最大的VOCs組分各不相同,依次為:烷烴、芳香烴及烯烴,但苯系物對三者均有很大貢獻(xiàn),尤其是對SOA貢獻(xiàn)最為顯著,是生成SOA的絕對優(yōu)勢物種.

    3.4 根據(jù)PMF受體模型的解析結(jié)果,南京北郊春、夏、秋三季的VOCs均有5個(gè)來源,冬季VOCs的來源有4個(gè),不同季節(jié)的源解析結(jié)果一致顯示工業(yè)排放和汽車尾氣是南京北郊最主要的VOCs污染源.

    3.5 富含苯系物的VOCs來源在不同季節(jié)中對SOA的貢獻(xiàn)均較大,如汽車尾氣排放源和溶劑使用源;富含乙烯、丙烯和異戊二烯的VOCs來源則對OFP的貢獻(xiàn)較為顯著,如石油化工源、汽車尾氣排放源和植物排放源.春、秋、冬三季汽車尾氣排放及工業(yè)排放(包括石化工業(yè))此二源對大氣環(huán)境中的VOCs濃度、VOCs對SOA的貢獻(xiàn)及VOCs對OFP的貢獻(xiàn)這三者的影響最為顯著,應(yīng)作為重點(diǎn)源頭加以控制,而夏季溶劑使用源及植物源因高溫而具備高揮發(fā)性,其對SOA及OFP的貢獻(xiàn)亦不容忽視.

    [1]Seinfeld J H, Pandis S N. Atmospheric Chemistry and Physics [M]. New York: Wiley, 1998.

    [2]Blake D R, Rowland F S. Urban leakage of liquefied petroleum gas and its impact on Mexico-City Air-Quality [J]. Science, 1995,269:953-956.

    [3]Derwent R G, Jenkin M E, Saunders S M, et al. Photochemical ozone formation in north west Europe and its contro1 [J]. Atmospheric Environment, 2003,37:1983-1991.

    [4]Xie X, Shao M, Liu Y, et al. Estimate of initial isoprene Contribution to ozone formation potential in Beijing, China [J]. Atmospheric Environment, 2008,42(24):6000-6010.

    [5]鄧雪嬌,王新明,趙春生,等.珠江三角洲典型過程VOCs的平均濃度與化學(xué)反應(yīng)活性 [J]. 中國環(huán)境科學(xué), 2010,30(9):1153-1161.

    [6]USEPA (United States Environmental Protection Agency). Carcinogenic effects of benzene: anupdate [R]. Washington: Prepared by the National Center for Environmental Health, Office of Research and Development, 2002.

    [7]Hansen A B, Palmgren F. VOC air pollutants in Copenhagen [J].The Science of the Total Environment, 1996,189/190:451-457.

    [8]Carter W P L. Development of ozone reactivity scales for volatile Organic compounds [J]. Air and Waste, 1994,44(7):881-899.

    [9]Grosjean D. In situ organic aerosol formation during a smog episode: estimated production and chemical functionality [J]. Atmospheric Environment, Part A, General Topics, 1992,26(6): 953-963.

    [10]Grosjean D, Seinfeld J H. Parameterization of the formation potential of secondary organic aerosols [J]. Atmospheric Environment (1967), 1989,23(8):1733-1747.

    [11]Bowman F M, Seinfeld J H. Fundamental basis of incremental reactivities of organics in ozone formation in VOC/NOx mixtures [J]. Atmospheric Environment, 1994,28:3359-3368.

    [12]吳 方,王躍思,安俊琳,等.北京奧運(yùn)時(shí)段VOCs濃度變化、臭氧產(chǎn)生潛勢及來源分析研究 [J]. 環(huán)境科學(xué), 2010,31(1):10-16.

    [13]陳長虹,蘇雷燕,王紅麗,等.上海市城區(qū)VOCs的年變化特征及其關(guān)鍵活性組分 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2012,32(2):367-376.

    [14]朱少峰,黃曉鋒,何凌燕,等.深圳大氣VOCs濃度的變化特征與化學(xué)反應(yīng)活性 [J]. 中國環(huán)境科學(xué), 2012,32(12):2140-2148.

    [15]安俊琳,朱 彬,李用宇.南京北郊大氣VOCs體積分?jǐn)?shù)變化特征[J]. 環(huán)境科學(xué), 2013,34(12):4504-4512.

    [16]李用宇,朱 彬,安俊琳,等.南京北郊秋季VOCs及其光化學(xué)特征觀測研究 [J]. 環(huán)境科學(xué), 2013,34(8):2933-2942.

    [17]鄭 玫,閆才青,李小瀅,等.二次有機(jī)氣溶膠估算方法研究進(jìn)展[J]. 中國環(huán)境科學(xué), 2014,34(3):555-564.

    [18]Turpin B J, Huntzicker J J. Identification of secondary organic aerosol episodes and qutization of primary and secondary organic aerosol concentrations during SCAQS [J]. Atmospheric Environment, 1995,29(23):3527-3544.

    [19]Hallquist M, Wenger J C, Baltensperger U, et al. The formation,properties and impact secondary organic aerosol:current and emerging issues [J]. Atmospheric Chemistry and Physics, 2009,9(14):5155-5236.

    [20]Carlton A G, Wiedinmyer C, Kroll J H. A review of secondary organic aerosol (SOA) formation from isoprene [J]. Atmospheric Chemistry and Physics, 2009, 9(14):4987-5005.

    [21]Alfarra M R, Paulsen D, Gysel M, et al. A mass spectrometric study of secondary organic aerosols formed from the photooxidation of anthropogenic and biogenic precursors in a reaction chamber [J]. Atmospheric Chemistry and Physics, 2006,6(12):5279-5293.

    [22]Robinson A L, Donahue N M, Shrivastava M K, et al. Rethinking organic aerosols:semivolatile emissions and photochemical Aging [J]. Science, 2007,315(5816):1259-1262.

    [23]George I J, Abbatt J P D. Chemical evolution of secondary organic aerosol from OH-initiated heterogeneous oxidation [J]. Atmospheric Chemistry and Physics, 2010,10(12):5551-5563.

    [24]Barthelmie R J, Pryor S C. Secondary organic aerosols: formation potential and ambient data [J]. The Science of the Total Environment, 1997,205(2/3):167-178.

    [25]Kourtidis K, Ziomas I. Estimation of secondary organic aerosol (SOA) production from traffic emissions in the city of Athens [J]. Global Nest, 1999,1(1):33-39.

    [26]Dechapanya W, Russell M, Allen D T. Estimates of anthropogenic secondary organic aerosol formation in Houston, Texas [J]. Aerosol Science and Technology, 2004,38(S1):156-166.

    [27]呂子峰,郝吉明,段菁春,等.北京市夏季二次有機(jī)氣溶膠生成潛勢的估算 [J]. 環(huán)境科學(xué), 2009,30(4):969-975.

    [28]王 倩,陳長虹,王紅麗,等.上海市秋季大氣VOCs對二次有機(jī)氣溶膠的生成貢獻(xiàn)及來源研究 [J]. 環(huán)境科學(xué), 2013,34(2):424-433.

    [29]王扶潘,朱 喬,馮 凝,等.深圳大氣中VOCs的二次有機(jī)氣溶膠生成潛勢 [J]. 中國環(huán)境科學(xué), 2014,34(10):2449-2457.

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

    [31]Brown S G, Frankel A, Hafner H R. Source apportionment of VOCs in the Los Angeles area using positive matrix factorization [J]. Atmospheric Environment, 2007,41(2):227-237.

    [32]Gary N,Ram V. EPA Positive matrix factorization (PMF) 3.0 Fundamentals&user guide [Z]. USA: EPA, 2008.

    [33]楊 輝,朱 彬,高晉徽,等.南京市北郊夏季揮發(fā)性有機(jī)物的源解析 [J]. 環(huán)境科學(xué), 2013,34(12):4519-4528.

    [34]Reff A, Eberly S I, Bhave P V. Receptor modeling of ambient particulate matter data using positive matrix factorization: review of existing methods [J]. Journal of the Air and Waste Management Association, 2007,57(2):146-154.

    [35]Cheng L, Fu L, Angle R P, et al. Seasonal variations of volatile organic compounds in Edmonton, Alberta [J]. Atmospheric Environment, 1997,31(2):239-246.

    [36]周裕敏,郝鄭平,王海林.北京地區(qū)城鄉(xiāng)結(jié)合部大氣揮發(fā)性有機(jī)物污染及來源分析 [J]. 環(huán)境科學(xué), 2011,32(12):3560-3565.

    [37]馮小瓊,彭 康,凌鎮(zhèn)浩,等.香港地區(qū)2005-2010年VOCs污染來源解析及特征研究 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2013,33(1):173-180.

    [38]夏 麗.南京江北地區(qū)VOCs濃度特征和來源解析 [D]. 南京:南京信息工程大學(xué), 2013.

    [39]葉偉紅,孫曉慧,劉勁松,等.杭州市交通干線空氣中揮發(fā)性有機(jī)物研究 [J]. 中國環(huán)境監(jiān)測, 2009,25(6):87-89.

    [40]Guo H, Wang T, Louie P K K. Source apportionment of ambient non-methane hydrocarbons in Hong Kong: application of a principal component analysis (PCA)/absolute principal component scores (APCS) receptor model [J]. Environmental Pollution, 2004,129(3):489-498.

    [41]陸思華,白郁華,張廣山,等.大氣中揮發(fā)性有機(jī)化合物(VOCs)的人為來源研究 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2006,26(5):757-763.

    Potential contribution of secondary organic aerosols and ozone of VOCs in the Northern Suburb of Nanjing.


    LIN Xu, ZHU Bin*, AN Jun-lin, YANG Hui (Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Nanjing University of Information Science and Technology, Nanjing 210044, China). China Environmental Science, 2015,35(4):976~986

    A continuous observation campaign was carried out with the GC5000 volatile organics online monitoring system from March 1, 2013 to February 28, 2014 in the northern suburb of Nanjing, characteristics of their composition,seasonal variation and diurnal variation were analyzed, PMF, the maximum incremental reactivity coefficient (MIR)and the fractional aerosol coefficients (FAC)were used to estimate the potential formation of secondary organic aerosols (SOA)and O3from VOCs and their sources. The results showed that the hourly average mixing ratio of the TVOCs was 45.63×10-9. There was an obvious seasonal cycle of VOCs, with the maximum in winter and autumn and minimum in summer. Diurnal variation of VOCs mixing ratios showed a very clear bimodal structure. The SOA concentration values obtained by the VOCs were 2.07μg/m3, as the largest contributor, aromatic hydrocarbons accounted for 95.93% and BTEX are the dominant species. Alkenes contributed the largest parts of the ozone formation potential (OFP), closing to 65%. Although alkanes were the most abundant components of VOCs, it is not the main contributor of OFP and SOA. The results from different seasons of receptor model showed that vehicle emissions and industrial emissions were main sources of VOCs in the northern suburb of Nanjing. The sources which contain rich BTEX contributed the largest parts of SOA. Moreover, the sources which contain rich ethylene, propylene and isoprene are the largest contributor of OFP. Vehicle emissions and industrial emissions (including the petrochemical industry) were the main contributor to the concentrations of VOCs, SOA and OFP in spring, autumn and winter. As the influential sources to SOA and OFP, solvent source and plant source should not be overlooked in summer.

    volatile organic compounds;SOA;OFP;sources;contribution

    X511

    A

    1000-6923(2015)04-0976-11

    林 旭(1989-),女,江蘇儀征人,碩士研究生,主要從事大氣環(huán)境與大氣化學(xué)研究.

    2014-09-15

    國家自然科學(xué)基金項(xiàng)目(41275143);江蘇省高校自然科學(xué)研究重大基礎(chǔ)研究項(xiàng)目(12KJA170003);江蘇省“333”高層次人才培養(yǎng)工程項(xiàng)目;江蘇省“六大人才高峰”計(jì)劃項(xiàng)目;國家自然科學(xué)基金項(xiàng)目(41305135)

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

    猜你喜歡
    北郊潛勢芳香烴
    立冬
    漢晉時(shí)期郊祀祭壇的形制演變
    賀蘭山沿山地區(qū)強(qiáng)對流潛勢預(yù)報(bào)方法研究
    臨渙選煤廠浮選效果評價(jià)與尾煤再選潛勢分析
    氨對甲苯二次有機(jī)氣溶膠生成潛勢影響的實(shí)驗(yàn)研究
    能源工程(2021年2期)2021-07-21 08:39:54
    系列芳香烴在土壤中洗脫規(guī)律的分子動(dòng)力學(xué)模擬
    三角形重心和垂心的向量表示的充要條件及其運(yùn)用
    “適者生存”?
    雜文月刊(2017年24期)2017-03-03 07:42:55
    先天性心臟病相關(guān)性肺動(dòng)脈高壓肺組織芳香烴受體的表達(dá)及其與肺血管重構(gòu)的相關(guān)性研究
    道路灰塵中多環(huán)芳香烴污染研究
    熟妇人妻不卡中文字幕| 激情 狠狠 欧美| 成年女人看的毛片在线观看| 久久久亚洲精品成人影院| 亚洲精品,欧美精品| 国产精品成人在线| 日韩中字成人| 黄片wwwwww| 久久久亚洲精品成人影院| 99热国产这里只有精品6| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 在线观看av片永久免费下载| 老司机影院毛片| 婷婷色麻豆天堂久久| 我的老师免费观看完整版| 国产 一区精品| 一区二区三区精品91| 亚洲综合色惰| 亚洲国产欧美人成| 国国产精品蜜臀av免费| 亚洲国产欧美人成| 亚洲精品日本国产第一区| 91精品国产九色| 久久久久久久国产电影| 秋霞伦理黄片| 国产老妇女一区| 99视频精品全部免费 在线| 99久久九九国产精品国产免费| 欧美人与善性xxx| 日韩,欧美,国产一区二区三区| 欧美三级亚洲精品| 亚洲av男天堂| 免费观看性生交大片5| 精品一区二区三区视频在线| 午夜福利视频1000在线观看| 国产v大片淫在线免费观看| 精品一区二区三卡| 亚洲经典国产精华液单| 99精国产麻豆久久婷婷| 深夜a级毛片| 美女主播在线视频| 青春草亚洲视频在线观看| 插阴视频在线观看视频| 久久久久国产精品人妻一区二区| 国产欧美日韩精品一区二区| 国产欧美日韩精品一区二区| 亚洲欧美精品自产自拍| 欧美日韩在线观看h| 少妇人妻 视频| 成人漫画全彩无遮挡| 男女下面进入的视频免费午夜| 精品国产三级普通话版| 久久人人爽av亚洲精品天堂 | 蜜桃亚洲精品一区二区三区| 亚洲不卡免费看| 精品人妻一区二区三区麻豆| 久久久久久久久久久免费av| 五月伊人婷婷丁香| 18禁在线播放成人免费| 一本久久精品| 中文字幕制服av| 久久亚洲国产成人精品v| 97超碰精品成人国产| 欧美日本视频| 亚洲精品乱码久久久久久按摩| 99久久精品热视频| 亚洲高清免费不卡视频| 69av精品久久久久久| 在线亚洲精品国产二区图片欧美 | 麻豆久久精品国产亚洲av| 久久久亚洲精品成人影院| 久久久久久伊人网av| 亚洲精品乱久久久久久| 国产成人a∨麻豆精品| 嫩草影院入口| 如何舔出高潮| 成人综合一区亚洲| 久久久久国产网址| 免费观看在线日韩| 日日啪夜夜撸| 少妇被粗大猛烈的视频| 亚洲内射少妇av| 18禁在线播放成人免费| 禁无遮挡网站| 高清在线视频一区二区三区| 国产老妇女一区| 日韩人妻高清精品专区| 黑人高潮一二区| 爱豆传媒免费全集在线观看| 亚洲va在线va天堂va国产| 一本一本综合久久| 国产精品久久久久久久久免| 国产欧美另类精品又又久久亚洲欧美| 性插视频无遮挡在线免费观看| 男男h啪啪无遮挡| 成人二区视频| 精品99又大又爽又粗少妇毛片| 老女人水多毛片| 久久久久国产网址| 国产爱豆传媒在线观看| a级毛片免费高清观看在线播放| 69人妻影院| 婷婷色综合www| 99热这里只有精品一区| 最近最新中文字幕大全电影3| 成人免费观看视频高清| 26uuu在线亚洲综合色| 久久综合国产亚洲精品| 五月天丁香电影| 熟妇人妻不卡中文字幕| 精品国产露脸久久av麻豆| 永久免费av网站大全| 亚洲国产精品专区欧美| 亚洲精品日本国产第一区| 真实男女啪啪啪动态图| 一级毛片黄色毛片免费观看视频| 国产伦精品一区二区三区视频9| 黄片wwwwww| av免费在线看不卡| 一级毛片久久久久久久久女| 毛片女人毛片| 成人美女网站在线观看视频| 在线精品无人区一区二区三 | 国产欧美日韩精品一区二区| 国产女主播在线喷水免费视频网站| 国产精品人妻久久久影院| 久久人人爽人人爽人人片va| 欧美一区二区亚洲| 国产精品一区www在线观看| 97精品久久久久久久久久精品| 国产成人a区在线观看| 麻豆成人av视频| freevideosex欧美| 亚洲综合色惰| av在线播放精品| 精品久久久久久久末码| 欧美丝袜亚洲另类| 久久精品国产亚洲av天美| 亚洲av中文字字幕乱码综合| 蜜臀久久99精品久久宅男| 国产爱豆传媒在线观看| 97在线视频观看| 女人十人毛片免费观看3o分钟| 国产淫片久久久久久久久| 国产精品久久久久久精品古装| 欧美老熟妇乱子伦牲交| 亚洲欧美精品自产自拍| 国产精品人妻久久久影院| 久久韩国三级中文字幕| 国产精品国产三级国产专区5o| 成人午夜精彩视频在线观看| 午夜福利网站1000一区二区三区| 一级a做视频免费观看| 日韩亚洲欧美综合| 国产视频内射| 亚洲欧美日韩东京热| 亚洲成人精品中文字幕电影| 女人十人毛片免费观看3o分钟| 国产精品99久久99久久久不卡 | 制服丝袜香蕉在线| 精品久久国产蜜桃| 男插女下体视频免费在线播放| av在线app专区| 国产男女超爽视频在线观看| 亚洲精品色激情综合| 一级毛片我不卡| 少妇熟女欧美另类| av专区在线播放| 五月开心婷婷网| 久久久久精品久久久久真实原创| 大香蕉97超碰在线| 国产探花在线观看一区二区| 亚洲天堂av无毛| 国产成人a∨麻豆精品| 日韩免费高清中文字幕av| 国产日韩欧美亚洲二区| 成年女人看的毛片在线观看| 啦啦啦在线观看免费高清www| 日本欧美国产在线视频| 免费看av在线观看网站| 成人午夜精彩视频在线观看| 哪个播放器可以免费观看大片| 日韩成人av中文字幕在线观看| 晚上一个人看的免费电影| 亚洲av在线观看美女高潮| 亚洲成色77777| 高清日韩中文字幕在线| 一级二级三级毛片免费看| 青春草国产在线视频| 在线播放无遮挡| 欧美日韩综合久久久久久| 亚洲一区二区三区欧美精品 | 男人添女人高潮全过程视频| 男女啪啪激烈高潮av片| av免费在线看不卡| 免费观看无遮挡的男女| 最近最新中文字幕大全电影3| 在线a可以看的网站| 亚洲欧洲国产日韩| 最后的刺客免费高清国语| 国产 精品1| 人人妻人人看人人澡| 国产极品天堂在线| 天美传媒精品一区二区| 国产淫片久久久久久久久| 亚洲精品中文字幕在线视频 | 国产精品99久久久久久久久| 两个人的视频大全免费| 亚洲国产欧美在线一区| 日本av手机在线免费观看| 久久久久久久大尺度免费视频| 三级经典国产精品| 99久久九九国产精品国产免费| 熟女电影av网| 成人国产av品久久久| 国产高清国产精品国产三级 | 秋霞在线观看毛片| 中文字幕人妻熟人妻熟丝袜美| 少妇的逼水好多| 午夜福利在线在线| 少妇裸体淫交视频免费看高清| 五月天丁香电影| 国产淫语在线视频| 日本午夜av视频| 日日摸夜夜添夜夜爱| 成年免费大片在线观看| 不卡视频在线观看欧美| 黄片无遮挡物在线观看| 久久99热这里只频精品6学生| 国产精品熟女久久久久浪| 成人综合一区亚洲| 夫妻性生交免费视频一级片| 国产精品人妻久久久影院| 久久久久久久精品精品| 欧美 日韩 精品 国产| 亚洲av二区三区四区| 精品人妻偷拍中文字幕| 九九久久精品国产亚洲av麻豆| 校园人妻丝袜中文字幕| 国产中年淑女户外野战色| 中文字幕免费在线视频6| 亚洲三级黄色毛片| 老司机影院毛片| 亚洲精品乱久久久久久| 久久精品久久久久久久性| 伊人久久精品亚洲午夜| 汤姆久久久久久久影院中文字幕| 国产免费视频播放在线视频| 久久精品熟女亚洲av麻豆精品| 亚洲真实伦在线观看| 亚洲aⅴ乱码一区二区在线播放| av黄色大香蕉| 成人毛片60女人毛片免费| 国产 一区精品| 免费不卡的大黄色大毛片视频在线观看| 日韩成人av中文字幕在线观看| 内地一区二区视频在线| 久久女婷五月综合色啪小说 | 国产久久久一区二区三区| 中文字幕av成人在线电影| 精品少妇久久久久久888优播| 国产精品麻豆人妻色哟哟久久| 欧美日韩在线观看h| 插阴视频在线观看视频| 欧美激情国产日韩精品一区| 大片电影免费在线观看免费| 老女人水多毛片| 国产黄色免费在线视频| 国模一区二区三区四区视频| 久久亚洲国产成人精品v| 精品一区二区免费观看| 少妇熟女欧美另类| 在线看a的网站| 国产精品久久久久久av不卡| 亚洲精品中文字幕在线视频 | 午夜激情久久久久久久| 色视频www国产| 波野结衣二区三区在线| 在线观看av片永久免费下载| 美女xxoo啪啪120秒动态图| 少妇裸体淫交视频免费看高清| av线在线观看网站| 亚洲av中文av极速乱| 91久久精品国产一区二区三区| 免费观看的影片在线观看| 亚洲av中文字字幕乱码综合| 亚洲欧美成人综合另类久久久| av女优亚洲男人天堂| 成人亚洲精品av一区二区| 亚洲欧美一区二区三区黑人 | 亚洲丝袜综合中文字幕| 成人毛片a级毛片在线播放| 男的添女的下面高潮视频| 综合色丁香网| 久久久久久久大尺度免费视频| 国内精品美女久久久久久| 一级毛片久久久久久久久女| 99热6这里只有精品| .国产精品久久| 又大又黄又爽视频免费| 国产精品久久久久久精品电影小说 | 精品午夜福利在线看| 亚洲av免费高清在线观看| 精品亚洲乱码少妇综合久久| 亚洲精品,欧美精品| 国产在视频线精品| 日韩,欧美,国产一区二区三区| xxx大片免费视频| 亚洲精品中文字幕在线视频 | 日韩欧美精品免费久久| 欧美日韩亚洲高清精品| 黄色日韩在线| 久久ye,这里只有精品| 国产在线一区二区三区精| 亚洲精品视频女| 男男h啪啪无遮挡| 人妻 亚洲 视频| 少妇丰满av| 亚洲成人久久爱视频| 国产精品99久久久久久久久| 中文精品一卡2卡3卡4更新| 欧美 日韩 精品 国产| 永久网站在线| 又大又黄又爽视频免费| 麻豆乱淫一区二区| 日韩欧美一区视频在线观看 | 亚洲人成网站在线观看播放| 亚洲精品成人久久久久久| 亚洲色图综合在线观看| 另类亚洲欧美激情| 亚洲av成人精品一二三区| 97精品久久久久久久久久精品| 国产高清不卡午夜福利| 中文精品一卡2卡3卡4更新| 国产精品国产三级国产专区5o| 亚洲,欧美,日韩| 亚洲欧美精品自产自拍| 欧美国产精品一级二级三级 | 啦啦啦中文免费视频观看日本| 精品久久久噜噜| 欧美精品人与动牲交sv欧美| 我要看日韩黄色一级片| 亚洲四区av| 亚洲丝袜综合中文字幕| 婷婷色麻豆天堂久久| 亚洲国产最新在线播放| 在线精品无人区一区二区三 | 日韩三级伦理在线观看| 国产免费一区二区三区四区乱码| 内射极品少妇av片p| 交换朋友夫妻互换小说| 王馨瑶露胸无遮挡在线观看| 欧美高清成人免费视频www| 成人二区视频| 国产午夜精品一二区理论片| 最近最新中文字幕免费大全7| 美女高潮的动态| 边亲边吃奶的免费视频| 欧美日韩视频高清一区二区三区二| 亚洲人成网站高清观看| 视频区图区小说| 欧美zozozo另类| 国产伦理片在线播放av一区| 免费观看a级毛片全部| 国产精品熟女久久久久浪| 成人免费观看视频高清| 国产黄片视频在线免费观看| 亚洲av中文字字幕乱码综合| 国产亚洲91精品色在线| 久久久精品94久久精品| 久久久久九九精品影院| 成人高潮视频无遮挡免费网站| 大片电影免费在线观看免费| 国产乱人偷精品视频| 一个人看的www免费观看视频| 欧美性感艳星| 成人二区视频| 国产成人freesex在线| 亚洲美女视频黄频| 女人被狂操c到高潮| 久久久亚洲精品成人影院| 高清欧美精品videossex| 99久久九九国产精品国产免费| 久久精品久久久久久噜噜老黄| 日韩欧美精品免费久久| 中文字幕人妻熟人妻熟丝袜美| 熟女电影av网| 69av精品久久久久久| 亚洲性久久影院| 精品一区二区免费观看| 18禁裸乳无遮挡免费网站照片| 人人妻人人看人人澡| 又粗又硬又长又爽又黄的视频| 亚洲久久久久久中文字幕| 久久久久国产精品人妻一区二区| 99精国产麻豆久久婷婷| 国产av国产精品国产| 日韩一区二区视频免费看| 深夜a级毛片| 三级国产精品片| 可以在线观看毛片的网站| 亚洲欧美一区二区三区黑人 | 久久99热这里只有精品18| av线在线观看网站| 国产一区亚洲一区在线观看| 老女人水多毛片| 亚洲怡红院男人天堂| 精品亚洲乱码少妇综合久久| 一本色道久久久久久精品综合| 亚洲av成人精品一区久久| 亚洲精品色激情综合| 成人综合一区亚洲| 国产精品秋霞免费鲁丝片| 日本三级黄在线观看| 日日撸夜夜添| 99久久人妻综合| 亚洲国产精品专区欧美| 真实男女啪啪啪动态图| 777米奇影视久久| 精品视频人人做人人爽| 久久久精品欧美日韩精品| 亚洲av成人精品一区久久| 看非洲黑人一级黄片| 国产淫片久久久久久久久| 国产老妇伦熟女老妇高清| 国产精品女同一区二区软件| 高清毛片免费看| 国产午夜精品久久久久久一区二区三区| 五月天丁香电影| 欧美日韩综合久久久久久| 深爱激情五月婷婷| 久久久久久久亚洲中文字幕| 国产黄频视频在线观看| 99热6这里只有精品| 纵有疾风起免费观看全集完整版| 亚洲va在线va天堂va国产| 亚洲av中文字字幕乱码综合| 少妇被粗大猛烈的视频| 99热全是精品| 亚洲第一区二区三区不卡| 嫩草影院精品99| 在线免费观看不下载黄p国产| 亚洲最大成人中文| 少妇人妻久久综合中文| 免费看日本二区| 内射极品少妇av片p| 大陆偷拍与自拍| 舔av片在线| 97超碰精品成人国产| 26uuu在线亚洲综合色| a级毛色黄片| 人人妻人人澡人人爽人人夜夜| 国产精品久久久久久久电影| 777米奇影视久久| 亚洲美女搞黄在线观看| 男女啪啪激烈高潮av片| 欧美成人a在线观看| 美女被艹到高潮喷水动态| 亚洲av国产av综合av卡| 亚洲欧美日韩另类电影网站 | 午夜视频国产福利| 亚洲国产精品专区欧美| 久久久色成人| 国产男女超爽视频在线观看| 国精品久久久久久国模美| 久久久欧美国产精品| 我的女老师完整版在线观看| 一级黄片播放器| 美女脱内裤让男人舔精品视频| 亚洲aⅴ乱码一区二区在线播放| 中文资源天堂在线| 欧美精品国产亚洲| 久热久热在线精品观看| 2021天堂中文幕一二区在线观| 国产美女午夜福利| 亚洲国产高清在线一区二区三| 嫩草影院入口| kizo精华| a级毛片免费高清观看在线播放| 亚洲精品中文字幕在线视频 | 美女高潮的动态| 午夜日本视频在线| 简卡轻食公司| 伦理电影大哥的女人| 综合色av麻豆| 免费黄频网站在线观看国产| 少妇被粗大猛烈的视频| 少妇 在线观看| 久久久久久九九精品二区国产| 国产在视频线精品| 99视频精品全部免费 在线| 成人漫画全彩无遮挡| 国产淫语在线视频| 一级片'在线观看视频| 亚洲精品一二三| 女人十人毛片免费观看3o分钟| 成人免费观看视频高清| 老司机影院成人| 国产乱人视频| 晚上一个人看的免费电影| 日产精品乱码卡一卡2卡三| 国内精品宾馆在线| 波野结衣二区三区在线| 国产男女内射视频| 国产亚洲av片在线观看秒播厂| 在线天堂最新版资源| 男的添女的下面高潮视频| 亚洲最大成人手机在线| 国产黄片美女视频| 一个人观看的视频www高清免费观看| 一级毛片aaaaaa免费看小| 在线观看三级黄色| 国产 一区 欧美 日韩| av又黄又爽大尺度在线免费看| 精品久久久噜噜| av免费在线看不卡| 久久久久久久大尺度免费视频| 蜜臀久久99精品久久宅男| 欧美zozozo另类| 九草在线视频观看| 黄色怎么调成土黄色| 国产精品99久久久久久久久| 亚洲av电影在线观看一区二区三区 | 国产探花在线观看一区二区| 色吧在线观看| 99热这里只有精品一区| av.在线天堂| 在线观看免费高清a一片| 看十八女毛片水多多多| 免费大片18禁| 色5月婷婷丁香| 欧美国产精品一级二级三级 | 最后的刺客免费高清国语| 十八禁网站网址无遮挡 | 汤姆久久久久久久影院中文字幕| 亚洲精品国产av成人精品| 女人久久www免费人成看片| 国产成人免费无遮挡视频| 赤兔流量卡办理| 视频区图区小说| 午夜精品国产一区二区电影 | 青春草亚洲视频在线观看| 伦精品一区二区三区| 国产精品爽爽va在线观看网站| 免费人成在线观看视频色| 中国三级夫妇交换| 久久精品国产亚洲av天美| 国产一区二区在线观看日韩| 久久精品久久久久久久性| 亚洲精品日韩av片在线观看| 麻豆精品久久久久久蜜桃| 日本-黄色视频高清免费观看| 激情 狠狠 欧美| 欧美97在线视频| 精品久久国产蜜桃| 男女无遮挡免费网站观看| 男女国产视频网站| 亚洲性久久影院| 午夜亚洲福利在线播放| 九草在线视频观看| 亚洲va在线va天堂va国产| 国产av码专区亚洲av| 男人爽女人下面视频在线观看| 午夜日本视频在线| 婷婷色av中文字幕| 久久人人爽人人爽人人片va| 伊人久久精品亚洲午夜| 色婷婷久久久亚洲欧美| 又大又黄又爽视频免费| 国产在线男女| 91精品伊人久久大香线蕉| 91午夜精品亚洲一区二区三区| 欧美日韩亚洲高清精品| 久久久精品94久久精品| 亚洲成色77777| 国产视频内射| 欧美精品国产亚洲| 久久久久国产网址| 亚洲国产色片| 搡女人真爽免费视频火全软件| 免费av观看视频| 六月丁香七月| 亚洲精品自拍成人| 真实男女啪啪啪动态图| 国精品久久久久久国模美| 亚洲精品乱码久久久久久按摩| 国产大屁股一区二区在线视频| 国产精品爽爽va在线观看网站| 夜夜看夜夜爽夜夜摸| 久久99热6这里只有精品| 日韩制服骚丝袜av| 午夜免费男女啪啪视频观看| 一个人观看的视频www高清免费观看| 亚洲电影在线观看av| 亚洲av免费高清在线观看| 欧美精品人与动牲交sv欧美| 国产精品久久久久久精品古装| av播播在线观看一区| videos熟女内射| 国产高清有码在线观看视频| 男人舔奶头视频| 中文在线观看免费www的网站| 在线观看国产h片| 国产午夜福利久久久久久| 内地一区二区视频在线| 日韩制服骚丝袜av| 日日摸夜夜添夜夜添av毛片| 午夜福利视频1000在线观看| 国产一区二区三区综合在线观看 | 国产一级毛片在线| 有码 亚洲区| 国产精品蜜桃在线观看| 国国产精品蜜臀av免费| 老司机影院毛片| 少妇 在线观看|