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

    2017年8—9月湖州市臭氧污染特征及其生成機(jī)制研究

    2020-05-30 19:22:55常艷春魏莉楊曉紅費(fèi)勇沈?qū)W勇陳江張海燕繆麗娜
    關(guān)鍵詞:湖州市貢獻(xiàn)臭氧

    常艷春 魏莉 楊曉紅 費(fèi)勇 沈?qū)W勇 陳江 張海燕 繆麗娜

    摘要近年來(lái)近地面臭氧問(wèn)題日益凸顯,成為影響空氣質(zhì)量持續(xù)改善的瓶頸.本研究基于2017年8—9月在湖州市城區(qū)開(kāi)展的為期1個(gè)月的臭氧及其前體物揮發(fā)性有機(jī)物(VOCs)和氮氧化物(NOx)在線觀測(cè)數(shù)據(jù),分析了臭氧及其前體物污染特征,利用正矩陣因子分析(PMF)解析了VOCs來(lái)源,并采用基于觀測(cè)的模型(OBM)對(duì)臭氧生成機(jī)制進(jìn)行研究.研究結(jié)果表明:1)觀測(cè)期間湖州市VOCs平均體積分?jǐn)?shù)為(24.78±9.10)×10-9,其中占比最高的組成為烷烴、含氧VOCs(OVOCs)和鹵代烴;2)在臭氧非超標(biāo)時(shí)段,湖州市臭氧生成處于VOCs控制區(qū),而在臭氧重污染期間湖州市處于以VOCs控制為主的過(guò)渡區(qū);3)在臭氧超標(biāo)時(shí)段,對(duì)臭氧生成潛勢(shì)(OFP)貢獻(xiàn)最大的是芳香烴(39.6%),其次是烯烴(21.5%)和OVOCs(19.4%),排名前三的關(guān)鍵組分為甲苯、乙烯和間/對(duì)二甲苯;4)源解析結(jié)果顯示觀測(cè)期間湖州市VOCs的主要來(lái)源是溶劑使用(27.0%)、交通排放(22.7%)、背景+傳輸(19.3%)、工業(yè)排放(16.9%)、汽油揮發(fā)(7.7%)和植物排放(6.4%),重污染過(guò)程期間對(duì)OFP貢獻(xiàn)最大的兩類源是交通排放源和溶劑使用源,貢獻(xiàn)百分比分別為35.1%和30.5%.因此,對(duì)交通排放和溶劑使用方面進(jìn)行控制管理對(duì)湖州市大氣臭氧污染防控有重要意義.關(guān)鍵詞

    臭氧;揮發(fā)性有機(jī)物;敏感性分析;臭氧生成潛勢(shì);源解析

    中圖分類號(hào) X51

    文獻(xiàn)標(biāo)志碼 A

    0 引言

    近年來(lái),我國(guó)大氣中可吸入顆粒物(PM10)、細(xì)顆粒物(PM2.5)、二氧化硫(SO2)和二氧化氮(NO2)濃度已經(jīng)呈現(xiàn)逐年下降趨勢(shì),但是以臭氧(O3)為代表的二次污染卻日益加劇,不僅表現(xiàn)出濃度水平的持續(xù)上升,而且呈現(xiàn)出以城市為中心向周邊地區(qū)蔓延的區(qū)域性特征[1-3],成為影響空氣質(zhì)量持續(xù)改善的瓶頸.近地面O3主要是由氮氧化物(NOx)和揮發(fā)性有機(jī)物(VOCs)在太陽(yáng)光照射下反應(yīng)生成的,對(duì)人體健康和生態(tài)環(huán)境都會(huì)造成不利影響[4-9].

    臭氧與前體物VOCs和NOx之間呈現(xiàn)高度的非線性響應(yīng)關(guān)系[10],因此準(zhǔn)確定量臭氧與前體物之間的關(guān)系是臭氧污染防治對(duì)策制定的基礎(chǔ)[11].不同環(huán)境中的臭氧對(duì)VOCs和NOx變化的響應(yīng)程度都有所不同,其響應(yīng)關(guān)系可用VOCs控制區(qū)、NOx控制區(qū)和過(guò)渡區(qū)來(lái)表征.當(dāng)臭氧的生成處于VOCs控制區(qū)時(shí),降低VOCs排放,能夠有效減少臭氧生成,從而達(dá)到控制臭氧污染的目的,而減少NOx排放效果不顯著,甚至?xí)?dǎo)致臭氧濃度的上升[12],反之,當(dāng)臭氧生成處于NOx控制區(qū)時(shí),減少NOx排放能有效降低臭氧濃度;過(guò)渡區(qū)則表示需要同時(shí)對(duì)VOCs和NOx排放進(jìn)行控制,才能達(dá)到臭氧污染防治的目的.

    研究表明,城市地區(qū)臭氧生成主要受VOCs控制[13],因此對(duì)大氣VOCs開(kāi)展監(jiān)測(cè),解析其來(lái)源,并進(jìn)行控制是當(dāng)前臭氧污染防治的重要手段之一.大氣中的VOCs成分復(fù)雜[14],不同組分之間都具有濃度和反應(yīng)活性的差異,且來(lái)源各不相同.城市環(huán)境中的VOCs主要來(lái)自于與燃料燃燒有關(guān)的機(jī)動(dòng)車尾氣和工業(yè)過(guò)程排放、溶劑的揮發(fā)、天然氣及液化石油氣等氣體泄漏和植被排放等[15-18],其本身各種復(fù)雜的理化性質(zhì)和來(lái)源貢獻(xiàn)給VOCs防控帶來(lái)了很大挑戰(zhàn).

    湖州位于長(zhǎng)江三角洲中心區(qū)域,VOCs排放量在浙江省排第6位[19].近年來(lái)湖州市能源消耗量、機(jī)動(dòng)車保有量和工業(yè)生產(chǎn)量逐步增加.2016—2017年,湖州市以O(shè)3為首要污染物的天數(shù)由70 d增加至75 d,O3已成為影響湖州市環(huán)境空氣質(zhì)量的主要因素.本研究基于2017年8—9月在湖州市城區(qū)開(kāi)展的為期1個(gè)月的臭氧及其前體物揮發(fā)性有機(jī)物(VOCs)和氮氧化物(NOx)在線觀測(cè)數(shù)據(jù),采用基于觀測(cè)的模型和正矩陣因子分析等方法,分析了臭氧及其前體物污染特征,選取臭氧重污染過(guò)程進(jìn)行成因分析,并開(kāi)展了VOCs來(lái)源解析,以期為臭氧污染防治措施的制定提供相關(guān)科學(xué)依據(jù).

    1 資料與方法

    1.1 觀測(cè)資料

    根據(jù)湖州市2016年臭氧日最大8 h第90百分位數(shù)的月變化特征,可以發(fā)現(xiàn)湖州市臭氧污染在7、8、9月最為嚴(yán)重.因此本研究選擇湖州市臭氧污染最嚴(yán)重的時(shí)段,于2017年8月23日—9月23日采用北京大學(xué)和武漢天虹環(huán)保產(chǎn)業(yè)股份有限公司共同研發(fā)的大氣環(huán)境揮發(fā)性有機(jī)物在線監(jiān)測(cè)系統(tǒng)(TH-300B)對(duì)湖州市大氣VOCs開(kāi)展連續(xù)觀測(cè).觀測(cè)站點(diǎn)位于湖州市城區(qū)(120.09°E ,30.89°N),點(diǎn)位周邊主要為行政辦公樓、住宅區(qū)、公園和景區(qū),人口密集,附近無(wú)工業(yè)企業(yè)分布,周邊交通系統(tǒng)較為發(fā)達(dá),屬于典型的中心城區(qū)點(diǎn)位,可代表湖州市城區(qū)大氣環(huán)境質(zhì)量.

    TH-300B采用雙通道采樣,兩路樣品分別在冷凍除水后進(jìn)入兩路捕集柱,在-160 ℃的條件下被冷凍富集.捕集柱被加熱到100 ℃后樣品進(jìn)入色譜柱中分離并分別用氫火焰離子化檢測(cè)器(FID)和質(zhì)譜(MS)進(jìn)行檢測(cè),其中FID檢測(cè)器檢出C2~C5的碳?xì)浠衔铮琈S檢測(cè)器檢出C5~C12的碳?xì)浠衔?、鹵代烴和含氧揮發(fā)性有機(jī)物[20].樣品分析的時(shí)間分辨率為60 min,檢測(cè)分析100種VOCs組分,主要包括28種烷烴、12種烯烴、1種炔烴、16種芳香烴、31種鹵代烴及12種含氧有機(jī)物.

    臭氧和氮氧化物分別采用紫外光度法和化學(xué)發(fā)光法原理進(jìn)行測(cè)量,使用的儀器為美國(guó)Thermo Fisher Scientific公司生產(chǎn)的49i型O3分析儀和42i型NOx分析儀.氣象資料由湖州市氣象局提供,包括風(fēng)向、風(fēng)速、溫度和相對(duì)濕度.

    1.2 模型與方法介紹

    基于觀測(cè)的模型(Observation-Based Model,OBM)是以臭氧及其前體物的實(shí)際觀測(cè)濃度,以及氣象資料作為輸入數(shù)據(jù)來(lái)模擬臭氧生成并分析臭氧與前體物關(guān)系的方法[21].通過(guò)計(jì)算臭氧前體物的相對(duì)增量反應(yīng)活性(RIR),即臭氧生成速率對(duì)時(shí)間的積分(PO3-NO)變化的百分比與特定前體物效應(yīng)變化百分比的比值來(lái)分析臭氧生成機(jī)制及其對(duì)NOx和VOCs的敏感性[22],公式如下:

    RIRS(X)=PO3-NO(X)-PO3-NO(X-ΔX)PO3-NO(X)ΔS(X)S(X),(1)

    式(1)中,X為某特定臭氧前體物,S(X)是前體物X的源效應(yīng),ΔS(X)是源效應(yīng)變化,ΔX代表由于源效應(yīng)變化造成的前體物X質(zhì)量濃度的變化[23].

    臭氧生成潛勢(shì)(Ozone Formation Potentials,OFPs)是基于最大增量反應(yīng)活性(Maximum Incremental Reactivity,MIR)來(lái)量化不同VOCs組分對(duì)臭氧生成貢獻(xiàn)的[24],具體公式如下:

    OFPi=MIRi×[VOC]i,(2)

    其中,[VOC]i表示觀測(cè)的VOCs組分i的體積分?jǐn)?shù)[25-26].

    正矩陣因子分析法(Positive Matrix Factorization,PMF)是一種基于組分觀測(cè)數(shù)據(jù)的多元統(tǒng)計(jì)分析方法,在VOCs源解析研究中得到廣泛應(yīng)用[27-28].該模型通過(guò)輸入觀測(cè)得到的環(huán)境數(shù)據(jù),分析各VOCs組分的變化規(guī)律,從而識(shí)別出主要的VOCs排放源及其化學(xué)組成特征,并計(jì)算各類排放源對(duì)環(huán)境大氣中VOCs濃度的貢獻(xiàn).數(shù)據(jù)矩陣X被分解為矩陣G(源貢獻(xiàn)矩陣)、矩陣F(源成分矩陣)和一個(gè)殘差矩陣E,Q是基于殘差的函數(shù),公式如下:

    eij=xij-∑pk=1gijfkj,(3)

    Q(E)=∑ni=1∑mj=1eijsij,(4)

    式中,xij,gij,fij,eij分別為矩陣X,G,F(xiàn),E的元素,在gij≥0,fij≥0的約束條件下,通過(guò)迭代最小化算法對(duì)Q求解,同時(shí)確定G和F[29].

    2 結(jié)果與討論

    2.1 觀測(cè)期間臭氧及VOCs污染特征分析

    圖 1展示了觀測(cè)期間臭氧及其前體物NOx和VOCs以及其他相關(guān)常規(guī)污染物濃度水平和風(fēng)速風(fēng)向的時(shí)間變化規(guī)律.觀測(cè)期間的主導(dǎo)風(fēng)向?yàn)闁|南偏東風(fēng),平均風(fēng)速為1.4 m/s,風(fēng)速較低,不利于污染物的稀釋擴(kuò)散和遠(yuǎn)距離輸送.湖州市臭氧小時(shí)質(zhì)量濃度最高值達(dá)286 μg/m3,臭氧小時(shí)超標(biāo)次數(shù)高達(dá)20次.臭氧日最大8小時(shí)滑動(dòng)平均值超過(guò)國(guó)家空氣質(zhì)量二級(jí)標(biāo)準(zhǔn)(160 μg/m3)的天數(shù)為9 d,其中9月16—18日為一次典型的臭氧重污染過(guò)程.

    觀測(cè)期間湖州市VOCs總體積分?jǐn)?shù)在7.55×10-9~78.29×10-9之間,整體平均體積分?jǐn)?shù)為(24.78±9.10)×10-9,8月的VOCs體積分?jǐn)?shù)變化較為平穩(wěn),9月的VOCs平均體積分?jǐn)?shù)26.67 ×10-9相比于8月的20.18×10-9有所上升,其中9月22日8時(shí)達(dá)到觀測(cè)期間VOCs體積分?jǐn)?shù)最高值78.29×10-9.VOCs體積分?jǐn)?shù)高值主要出現(xiàn)在風(fēng)速1~2 m/s的氣象條件下,風(fēng)向?yàn)闁|南偏東和西北偏北方向.氮氧化物也呈現(xiàn)出相似的濃度變化特征,9月的NO、NO2平均濃度水平較8月有明顯上升.總體來(lái)看,NOx和VOCs呈現(xiàn)出與臭氧濃度相反的日變化,前者通常在晚上和清晨出現(xiàn)高值,而臭氧在中午出現(xiàn)高值.

    VOCs化學(xué)組成復(fù)雜,會(huì)受到排放源和化學(xué)反應(yīng)的影響.圖2統(tǒng)計(jì)了觀測(cè)期間各類組分在總VOCs體積分?jǐn)?shù)的平均占比.可以看出,烷烴的體積分?jǐn)?shù)值最高為8.17×10-9,占總VOCs體積分?jǐn)?shù)的33.6%;其次是OVOCs,體積分?jǐn)?shù)為5.28×10-9,在總VOCs體積分?jǐn)?shù)中的占比為21.7%;然后是鹵代烴,體積分?jǐn)?shù)為4.05×10-9,占比為16.7%;芳香烴體積分?jǐn)?shù)略低于鹵代烴,為3.53×10-9,占比14.6%;烯烴、乙炔和異戊二烯的體積分?jǐn)?shù)較低,分別占總VOCs體積分?jǐn)?shù)的8.4%、4.4%和0.6%.

    2.2 典型臭氧污染過(guò)程成因分析

    2.2.1 臭氧污染過(guò)程的選取

    在整個(gè)觀測(cè)期間選取一次典型的臭氧重污染過(guò)程(9月16—18日)進(jìn)行更為細(xì)致的分析.圖3展示了此次重污染過(guò)程前后一周的相關(guān)污染物濃度及氣象要素的變化情況,按臭氧質(zhì)量濃度是否超標(biāo)將整個(gè)期間分為超標(biāo)日(09-16—18)和非超標(biāo)日(09-14—15、09-19—21)進(jìn)行分析.在臭氧超標(biāo)日平均風(fēng)速為1.2 m/s,風(fēng)向以西北風(fēng)向和東南風(fēng)向?yàn)橹?,處于靜穩(wěn)天氣,因而本地污染物不易擴(kuò)散.從溫度和相對(duì)濕度角度看,重污染過(guò)程前后一周臭氧小時(shí)濃度的變化趨勢(shì)與溫度的變化趨勢(shì)較為一致,與相對(duì)濕度的變化趨勢(shì)相反,這是因?yàn)楦邷亻g接反映出高輻射強(qiáng)度,輻射越強(qiáng),光化學(xué)反應(yīng)越劇烈,臭氧越易生成[30],而相對(duì)濕度體現(xiàn)的是降水的影響,尤其在9月19日之后湖州出現(xiàn)連續(xù)降雨天氣,此時(shí)氣溫亦下降,導(dǎo)致臭氧濃度均持續(xù)下降.

    在臭氧連續(xù)超標(biāo)期間,臭氧峰值出現(xiàn)在下午的16—17時(shí),這三天的臭氧峰值分別達(dá)到了184、214和278 μg/m3,甚至出現(xiàn)了夜間高值的現(xiàn)象.這一污染過(guò)程風(fēng)速較低,且處在臺(tái)風(fēng)“泰利”的外圍,易造成污染物的積聚,另外高溫低濕的氣象條件為臭氧生成提供了良好條件.從前體物角度來(lái)看,VOCs和NOx于夜間至凌晨的排放濃度也明顯升高,之后在白天的光化學(xué)反應(yīng)消耗作用下濃度降低,可以看到,前體物排放濃度越高,之后的臭氧超標(biāo)情況也越嚴(yán)重.充足的前體物加上適宜的反應(yīng)條件,造成了大氣臭氧濃度升高.

    2.2.2 臭氧生成機(jī)制分析

    為了評(píng)估前體物VOCs和NOx的變化對(duì)臭氧生成的影響,本文利用OBM模型定量分析湖州市重污染期間臭氧與前體物VOCs和NOx之間的非線性關(guān)系.設(shè)置441個(gè)VOCs和NOx減排情景,并計(jì)算不同情景下臭氧日最大8小時(shí)滑動(dòng)平均質(zhì)量濃度,來(lái)繪制臭氧質(zhì)量濃度等值線(EKMA曲線),結(jié)果如圖4所示.

    從圖4中可以看出,超標(biāo)日和非超標(biāo)日臭氧日最大8小時(shí)滑動(dòng)平均值相差較大,非超標(biāo)日臭氧質(zhì)量濃度約為90 μg/m3,超標(biāo)日臭氧質(zhì)量濃度約為188 μg/m3.對(duì)于非超標(biāo)日而言,僅削減VOCs時(shí),隨著削減比例的增大,臭氧濃度逐漸降低;僅削減NOx時(shí),臭氧濃度反而逐漸升高,尤其當(dāng)削減比例升至40%時(shí),臭氧日最大8小時(shí)滑動(dòng)平均值到達(dá)最大值115 μg/m3,隨后需繼續(xù)增大NOx削減比例,臭氧濃度才逐步降低,這說(shuō)明非超標(biāo)日湖州市處于VOCs控制區(qū),與大部分國(guó)內(nèi)研究結(jié)果一致[31-32].而在超標(biāo)日期間,僅削減VOCs時(shí),隨著削減比例的增加臭氧逐漸降低;僅削減NOx時(shí),當(dāng)削減比例達(dá)到10%左右,臭氧質(zhì)量濃度會(huì)略微上升至189 μg/m3,然后NOx削減比例繼續(xù)增大,臭氧濃度就會(huì)明顯下降,說(shuō)明在超標(biāo)日期間,臭氧生成處于以VOCs控制為主的過(guò)渡區(qū),需要VOCs和NOx協(xié)同控制減排,才能有效實(shí)現(xiàn)臭氧污染防治.

    2.3 前體物VOCs關(guān)鍵組分識(shí)別

    VOCs是大氣中的重要污染物,對(duì)臭氧和二次有機(jī)氣溶膠的生成都有著很大影響,且大氣環(huán)境中的VOCs各組分體積分?jǐn)?shù)、化學(xué)反應(yīng)活性差異較大,對(duì)臭氧生成的貢獻(xiàn)程度也參差不齊.如An等在南京的研究發(fā)現(xiàn)[33],VOCs各組分構(gòu)成中對(duì)臭氧生成貢獻(xiàn)最為主要的是烯烴,而Shao等在珠三角地區(qū)的研究結(jié)果表明,芳香烴對(duì)臭氧生成的貢獻(xiàn)更為顯著[34].

    本研究采用臭氧生成潛勢(shì)OFP方法來(lái)量化VOCs中各組分對(duì)臭氧生成的貢獻(xiàn),識(shí)別湖州市大氣中的關(guān)鍵VOCs組分,作為控制近地面臭氧濃度的優(yōu)先考慮物種.根據(jù)式(2)計(jì)算臭氧生成潛勢(shì),得到整個(gè)觀測(cè)期間VOCs各化學(xué)組成的OFP占比情況如圖5所示.可以看出,非超標(biāo)日中,芳香烴對(duì)OFP的貢獻(xiàn)最大,占比為43.6%,其次為烯烴,OFP占比為19.7%,OVOCs的OFP占比略低于烯烴,為18.2%,然后是烷烴,OFP占比為13.1%,天然源、鹵代烴和炔烴的貢獻(xiàn)較小,OFP占比分別為3.0%、1.6%和0.8%.9月16—18日臭氧超標(biāo)期間,湖州市VOCs各化學(xué)組成的OFP占比中,芳香烴依然對(duì)臭氧生成貢獻(xiàn)最大(39.6%),其次為烯烴(21.5%)、OVOCs(19.4%)和烷烴(14.0%),貢獻(xiàn)比例最小的仍然為天然源、鹵代烴和炔烴,分別為2.8%、1.7%和1.0%.VOCs各化學(xué)組成類別的OFP貢獻(xiàn)占比及其排序與非超標(biāo)日相比沒(méi)有明顯變化.

    圖6展示了臭氧污染超標(biāo)期間和非超標(biāo)期間VOCs對(duì)臭氧生成潛勢(shì)貢獻(xiàn)前十的組分.在臭氧重污染過(guò)程中,貢獻(xiàn)最大的10個(gè)關(guān)鍵組分依次是甲苯、乙烯、間/對(duì)二甲苯、丙醛、甲基乙基酮、丙烯、乙苯、鄰二甲苯、異戊二烯和異戊烷.與非超標(biāo)時(shí)段進(jìn)行比較,可以發(fā)現(xiàn)對(duì)臭氧生成起到關(guān)鍵作用的前十種組分基本一致,只排名順序有所差別;貢獻(xiàn)排名前三的組分均為甲苯、間/對(duì)二甲苯和乙烯,其中甲苯一直是超標(biāo)時(shí)段和非超標(biāo)時(shí)段對(duì)臭氧生成貢獻(xiàn)排名第一的物種.

    2.4 VOCs重點(diǎn)排放源識(shí)別

    2.4.1 觀測(cè)期間源解析結(jié)果分析

    利用PMF模型解析VOCs來(lái)源,得出如圖7所示的因子譜圖.因子1中二氯甲烷、三氯甲烷、1,2-二氯乙烷為優(yōu)勢(shì)物種,乙炔和苯也有較高的因子解釋度,這些組分在鋼鐵、煉焦、基礎(chǔ)化學(xué)原料制造等行業(yè)中存在較多排放[35].因此因子1識(shí)別為工業(yè)排放源.因子2中主要組分為C7~C8芳香烴,包括甲苯、乙苯、間/對(duì)二甲苯和鄰二甲苯,此外,支鏈高碳烷烴的解釋度也較高,這可能與當(dāng)?shù)氐男袠I(yè)生產(chǎn)特征有關(guān).因涂料和溶劑使用過(guò)程中會(huì)排放大量的芳香烴[36],所以將因子2識(shí)別為與溶劑使用相關(guān)的排放源.

    因子3中主要組分是乙烷、丙烷、丙酮和2-丁酮等壽命較長(zhǎng)的穩(wěn)定組分,可以代表區(qū)域的背景濃度;此外,丙酮和2-丁酮是二次生成的產(chǎn)物[37],可能來(lái)自其他地區(qū)長(zhǎng)時(shí)間的輸送,因此將因子3識(shí)別為背景+傳輸源.

    因子4中主要物種為異戊二烯.異戊二烯在夏季主要來(lái)自植被排放,冬季可能來(lái)自交通和工業(yè)排放[38],由于在線監(jiān)測(cè)時(shí)間是8—9月(夏季),因此將因子4歸屬為植物排放源.

    因子5中優(yōu)勢(shì)組分為異戊烷、正戊烷,還有部分甲基叔丁基醚(MTBE).MTBE是一種高辛烷值汽油添加劑,主要存在于汽油產(chǎn)品中.汽油揮發(fā)的最主要組分即為異戊烷,其次是正戊烷、正丁烷和異丁烷[28].因此,因子5表征的是油品揮發(fā)源.

    因子6中既有C2~C6烷烴,也有乙烯、丙烯、1-丁二烯等C4烯烴、乙炔、MTBE、苯和甲苯.異戊烷、正戊烷和異丁烷等同時(shí)存在于汽油揮發(fā)和機(jī)動(dòng)車尾氣中[39],而且因子6中還含有汽油燃燒的產(chǎn)物乙烯、丙烯、乙炔、1-丁烯等[40],因而因子6具有明顯的交通排放特征.

    各類源對(duì)總VOCs的貢獻(xiàn)百分比如圖8所示.可以看出,溶劑使用源對(duì)整個(gè)觀測(cè)期間的VOCs體積分?jǐn)?shù)貢獻(xiàn)最高,占比為27.0%,其次是交通排放源(22.7%)和背景+傳輸源(19.3%),以鹵代烴為代表的工業(yè)排放源占比為16.9%,汽油揮發(fā)源和植物排放源的貢獻(xiàn)最低,占比分別為7.7%和6.4%.

    2.4.2 臭氧污染過(guò)程源解析結(jié)果分析

    臭氧污染期間VOCs來(lái)源解析結(jié)果如圖9所示,分別表示6類排放源對(duì)VOCs體積分?jǐn)?shù)和臭氧生成潛勢(shì)OFP的相對(duì)貢獻(xiàn).從圖9中可以看出,重污染過(guò)程期間,工業(yè)排放源和溶劑使用源對(duì)湖州市VOCs體積分?jǐn)?shù)的貢獻(xiàn)最大,分別占比26.0%和25.8%,其次是交通排放,占比為19.2%,汽油揮發(fā)源(12.6%)和背景+傳輸源(11.8%)兩者的貢獻(xiàn)相當(dāng),植被源的貢獻(xiàn)最低,僅為4.6%.

    對(duì)OFP而言,交通排放源和溶劑使用源是最主要的兩大VOCs貢獻(xiàn)源,貢獻(xiàn)百分比分別為35.1%和30.5%.工業(yè)排放源盡管是VOCs體積分?jǐn)?shù)的最大來(lái)源,但其對(duì)OFP的貢獻(xiàn)比例較低,僅有14.2%,這是因?yàn)楣I(yè)排放源的示蹤物主要是鹵代烴,對(duì)臭氧生成的影響不大[41].與工業(yè)排放源相反,植被源對(duì)VOCs體積分?jǐn)?shù)貢獻(xiàn)最小,但其對(duì)OFP的貢獻(xiàn)占比幾乎是VOCs體積分?jǐn)?shù)占比的2倍,這是因?yàn)橹参镌吹氖聚櫸锂愇於┰诖髿庵械姆磻?yīng)活性較高[38],較易生成臭氧.因此,針對(duì)湖州市的臭氧污染過(guò)程,應(yīng)優(yōu)先對(duì)交通排放和溶劑使用源VOCs進(jìn)行控制.

    3 結(jié)論

    1)觀測(cè)期間湖州市VOCs總體積分?jǐn)?shù)在7.55×10-9~78.29 ×10-9之間,整體平均體積分?jǐn)?shù)為(24.78±9.10)×10-9,VOCs體積分?jǐn)?shù)高值主要出現(xiàn)在風(fēng)速1~2 m/s的氣象條件下,烷烴、OVOCs和鹵代烴是VOCs化學(xué)組成中占比最高的類別,所占的比例分別為33.6%、21.7%和16.7%.

    2)OBM模型結(jié)果顯示非超標(biāo)日情況下湖州市為典型的VOCs控制區(qū),而臭氧重污染期間湖州市處于以VOCs控制為主的過(guò)渡區(qū),需要VOCs和NOx協(xié)同控制減排,才能有效實(shí)現(xiàn)臭氧污染防治.

    3)無(wú)論是超標(biāo)時(shí)段還是非超標(biāo)時(shí)段,VOCs化學(xué)組成對(duì)OFP貢獻(xiàn)大小排序均一致,依次為芳香烴>烯烴>OVOCs>烷烴>天然源>鹵代烴>炔烴;排名前三的VOCs關(guān)鍵組分均為甲苯、間/對(duì)二甲苯和乙烯,其中甲苯一直是對(duì)臭氧生成貢獻(xiàn)排名第一的物種.

    4)PMF解析結(jié)果顯示觀測(cè)期間VOCs的主要來(lái)源是工業(yè)排放、溶劑使用、背景+傳輸、植物排放、汽油揮發(fā)和交通排放.臭氧重污染期間,工業(yè)排放源和溶劑使用源對(duì)VOCs體積分?jǐn)?shù)的貢獻(xiàn)最大,分別占比26.0%和25.8%,但計(jì)算各來(lái)源的OFP發(fā)現(xiàn),交通排放源和溶劑使用源對(duì)OFP貢獻(xiàn)最高,貢獻(xiàn)百分比分別為35.1%和30.5%.因此,對(duì)交通排放和溶劑使用進(jìn)行優(yōu)先控制對(duì)湖州市大氣臭氧污染防控有重要意義.

    參考文獻(xiàn)

    References

    [1]

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

    TANG Xiaoyan,ZHANG Yuanhang,SHAO Min.Atmospheric environmental chemistry[M].2nd ed.Beijing:Higher Education Press,2006:272-273

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

    ZHU Tong.Combined pollution of urban and regional atmospheres[M].Beijing:Chemical Industry Press,2005:1-3

    [3] 梁碧玲,張麗,賴鑫,等.深圳市臭氧污染特征及其與氣象條件的關(guān)系[J].氣象與環(huán)境學(xué)報(bào),2017,33(1):66-71

    LIANG Biling,ZHANG Li,LAI Xin,et al.Analysis of the characteristics of ozone pollution and its relationship with meteorological conditions in Shenzhen[J].Journal of Meteorology and Environment,2017,33(1):66-71

    [4] 張紅星,孫旭,姚余輝,等.北京夏季地表臭氧污染分布特征及其對(duì)植物的傷害效應(yīng)[J].生態(tài)學(xué)報(bào),2014,34(16):4756-4765

    ZHANG Hongxing,SUN Xu,YAO Yuhui,et al.Ground-level ozone distribution pattern in summer of Beijing and its foliar injury effect upon plants[J].Acta Ecologica Sinica,2014,34(16):4756-4765

    [5] 趙輝,鄭有飛,魏莉,等.南京大氣臭氧濃度的季節(jié)變化及其對(duì)主要作物影響的評(píng)估[J].環(huán)境科學(xué),2018,39(7):3418-3425

    ZHAO Hui,ZHENG Youfei,WEI Li,et al.Seasonal variation in surface ozone and its effect on the winter wheat and rice in Nanjing,China[J].Environmental Science,2018,39(7):3418-3425

    [6] Francis M,Groves A M,Sun R,et al.CCR2 regulates inflammatory cell accumulation in the lung and tissue injury following ozone exposure[J].Toxicological Sciences,2017,155(2):474-484

    [7] Khaniabadi Y O,Hopke P K,Goudarzi G,et al.Cardiopulmonary mortality and COPD attributed to ambient ozone[J].Environmental Research,2017,152:336-341

    [8] Wang T,Xue L K,Brimblecombe P,et al.Ozone pollution in China:a review of concentrations,meteorological influences,chemical precursors,and effects[J].Science of the Total Environment,2017,575:1582-1596

    [9] Mumawa C L,Levesque S,McGraw C,et al.Microglial priming through the lung-brain axis:the role of air pollution-induced circulating factors[J].FASEB Journal,2016,30(5):1880-1891

    [10] 羅恢泓,袁自冰,鄭君瑜,等.上海夏季臭氧生成機(jī)制時(shí)空變化特征及其影響因素研究[J].環(huán)境科學(xué)學(xué)報(bào),2019,39(1):154-168

    LUO Huihong,YUAN Zibing,ZHENG Junyu,et al.Spatio-temporal variation of summertime ozone formation mechanism in Shanghai and its impact factors[J].Acta Scientiae Circumstantiae,2019,39(1):154-168

    [11] Ou J M,Yuan Z B,Zheng J Y,et al.Ambient ozone control in a photochemically active region:short-term despiking or long-term attainment?[J].Environmental Science & Technology,2016,50(11):5720-5728

    [12] 謝瑞加,陳璋琪.泉州市區(qū)臭氧污染成因初步分析[J].環(huán)境與可持續(xù)發(fā)展,2019,44(1):56-59

    XIE Ruijia,CHEN Zhangqi.Preliminary analysis on the causes of ozone pollution in Quanzhou city[J].Environment and Sustainable Development,2019,44(1):56-59

    [13] 黃美元,徐華英,王庚辰.大氣環(huán)境學(xué)[M].北京:高等教育出版社,2005:214-263

    HUANG Meiyuan,XU Huaying,WANG Gengchen.Atmospheric environment[M].Beijing:Higher Education Press,2005:214-263

    [14] Wang B,Shao M,Lu S H,et al.Variation of ambient non-methane hydrocarbons in Beijing city in summer 2008[J].Atmospheric Chemistry and Physics,2010,10(13):5911-5923

    [15] Kansal A.Sources and reactivity of NMHCs and VOCs in the atmosphere:a review[J].Journal of Hazardous Materials,2009,166(1):17-26

    [16] Huang C,Chen C H,Li L,et al.Emission inventory of anthropogenic air pollutants and VOC species in the Yangtze River Delta region,China[J].Atmospheric Chemistry and Physics,2011,11(9):4105-4120

    [17] An J L,Wang Y S,Wu F K,et al.Characterizations of volatile organic compounds during high ozone episodes in Beijing,China[J].Environmental Monitoring and Assessment,2012,184(4):1879-1889

    [18] Zhang Y J,Mu Y J,Liu J F,et al.Levels,sources and health risks of carbonyls and BTEX in the ambient air of Beijing,China[J].Journal of Environmental Sciences,2012,24(1):124-130

    [19] 吳建,程文,晁娜.浙江省人為源VOCs排放清單[J].上海環(huán)境科學(xué),2015(2):60-66

    WU Jian,CHENG Wen,CHAO Na.An inventory for VOCs emissions from anthropogenic sources in Zhejiang province[J].Shanghai Environmental Sciences,2015(2):60-66

    [20] Wang M,Zeng L M,Lu S H,et al.Development and validation of a cryogen-free automatic gas chromatograph system (GC-MS/FID) for online measurements of volatile organic compounds[J].Analytical Methods,2014,6(23):9424-9434

    [21] 蔣美青,陸克定,蘇榕,等.我國(guó)典型城市群O3污染成因和關(guān)鍵VOCs活性解析[J].科學(xué)通報(bào),2018,63(12):1130-1141

    JIANG Meiqing,LU Keding,SU Rong,et al.Ozone formation and key VOCs in typical Chinese city clusters[J].Chinese Science Bulletin,2018,63(12):1130-1141

    [22] 蘇榕,陸克定,余家燕,等.基于觀測(cè)模型的重慶大氣臭氧污染成因與來(lái)源解析[J].中國(guó)科學(xué):地球科學(xué),2018,48(1):102-112

    SU Rong,LU Keding,YU Jiayan,et al.Exploration of the formation mechanism and source attribution of ambient ozone in Chongqing with an observation-based model[J].Scientia Sinica (Terrae),2018,48(1):102-112

    [23] Cardelino C A,Chameides W L.An observation-based model for analyzing ozone precursor relationships in the urban atmosphere[J].Journal of the Air & Waste Management Association,1995,45(3):161-180

    [24] Zeinali M,McConnell L L,Hapeman C J,et al.Volatile organic compounds in pesticide formulations:methods to estimate ozone formation potential[J].Atmospheric Environment,2011,45(14):2404-2412

    [25] 丁潔然,景長(zhǎng)勇.唐山夏季大氣VOCs污染特征及臭氧生成潛勢(shì)[J].環(huán)境工程,2016,34(6):130-135

    DING Jieran,JING Changyong.Pollution characteristics and ozone formation potential of ambient VOCs in summer in Tangshan[J].Environmental Engineering,2016,34(6):130-135

    [26] 李友平,唐婭,范忠雨,等.成都市大氣環(huán)境VOCs污染特征及其健康風(fēng)險(xiǎn)評(píng)價(jià)[J].環(huán)境科學(xué),2018,39(2):576-584

    LI Youping,TANG Ya,F(xiàn)AN Zhongyu,et al.Pollution characteristics and health risk assessment of atmospheric VOCs in Chengdu[J].Environmental Science,2018,39(2):576-584

    [27] Liu C H,Tian F L,Chen J W,et al.A comparative study on source apportionment of polycyclic aromatic hydrocarbons in sediments of the Daliao River,China:positive matrix factorization and factor analysis with non-negative constraints[J].Chinese Science Bulletin,2010,55(10):915-920

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

    YANG Hui,ZHU Bin,GAO Jinhui,et al.Source apportionment of VOCs in the northern suburb of Nanjing in summer[J].Environmental Science,2013,34(12):4519-4528

    [29] Brown S G,F(xiàn)rankel 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

    [30] 嚴(yán)文蓮,劉端陽(yáng),康志明,等.江蘇臭氧污染特征及其與氣象因子的關(guān)系[J].氣象科學(xué),2019,39(4):477-487

    YAN Wenlian,LIU Duanyang,KANG Zhiming,et al.The characteristics of ozone pollution and its relationship with meteorological factors in Jiangsu[J].Journal of the Meteorological Sciences,2019,39(4):477-487

    [31] Xue L K,Wang T,Gao J,et al.Ground-level ozone in four Chinese cities:precursors,regional transport and heterogeneous processes[J].Atmospheric Chemistry and Physics,2014,14(23):13175-13188

    [32] 張亮,朱彬,高晉徽,等.長(zhǎng)江三角洲夏季一次典型臭氧污染過(guò)程的模擬[J].環(huán)境科學(xué),2015,36(11):3981-3988

    ZHANG Liang,ZHU Bin,GAO Jinhui,et al.Modeling study of a typical summer ozone pollution event over Yangtze River Delta[J].Environmental Science,2015,36(11):3981-3988

    [33] An J L,Zou J N,Wang J X,et al.Differences in ozone photochemical characteristics between the megacity Nanjing and its suburban surroundings,Yangtze River Delta,China[J].Environmental Science and Pollution Research,2015,22(24):19607-19617

    [34] Shao M,Zhang Y H,Zeng L M,et al.Ground-level ozone in the Pearl River Delta and the roles of VOC and NOx in its production[J].Journal of Environmental Management,2009,90(1):512-518

    [35] Srivastava A,Sengupta B,Dutta S A.Source apportionment of ambient VOCs in Delhi city[J].Science of the Total Environment,2005,343(1/2/3):207-220

    [36] Lyu X P,Chen N,Guo H,et al.Ambient volatile organic compounds and their effect on ozone production in Wuhan,Central China[J].Science of the Total Environment,2016,541:200-209

    [37] 王伯光,劉燦,呂萬(wàn)明,等.廣州大氣揮發(fā)性醛酮類化合物的污染特征及來(lái)源研究[J].環(huán)境科學(xué),2009,30(3):631-636

    WANG Boguang,LIU Can,LU Wanming,et al.Character of volatile carbonyl compounds and their source in Guangzhou ambient air[J].Environmental Science,2009,30(3):631-636

    [38] 鄒宇,鄧雪嬌,李菲,等.廣州大氣中異戊二烯濃度變化特征、化學(xué)活性和來(lái)源分析[J].環(huán)境科學(xué)學(xué)報(bào),2015,35(3):647-655

    ZOU Yu,DENG Xuejiao,LI Fei,et al.Variation characteristics,chemical reactivity and sources of isoprene in the atmosphere of Guangzhou[J].Acta Scientiae Circumstantiae,2015,35(3):647-655

    [39] 邵敏,付琳琳,劉瑩,等.北京市大氣揮發(fā)性有機(jī)物的關(guān)鍵活性組分及其來(lái)源[J].中國(guó)科學(xué):地球科學(xué),2005,35(增1):123-130

    SHAO Min,F(xiàn)U Linlin,LIU Ying,et al.Key active components and their sources of volatile organic compounds in Beijing[J].Science in China (Series D),2005,35(sup1):123-130

    [40] Cai C J,Geng F H,Tie X X,et al.Characteristics and source apportionment of VOCs measured in Shanghai,China[J].Atmospheric Environment,2010,44(38):5005-5014

    [41] 周子航,鄧也,譚欽文,等.四川省人為源揮發(fā)性有機(jī)物組分清單及其臭氧生成潛勢(shì)[J].環(huán)境科學(xué),2019,40(4):1613-1626

    ZHOU Zihang,DENG Ye,TAN Qinwen,et al.Speciated VOCs emission inventory and ozone formation potential in Sichuan province[J].Environmental Science,2019,40(4):1613-1626

    Ozone pollution and its formation mechanism in Huzhou during 2017

    CHANG Yanchun1 WEI Li2 YANG Xiaohong1 FEI Yong1 SHEN Xueyong2

    CHEN Jiang1 ZHANG Haiyan1 MIAO Lina1

    1 Huzhou Environmental Protection Monitoring Center Station,Huzhou 313000

    2 Nanjing Intelligent Environmental Sci-Tech Company Limited,Nanjing 211800

    Abstract Surface ozone pollution has become a bottleneck affecting the continuous improvement of air quality in recent years.Based on online observation data of ozone and its precursors including volatile organic compounds (VOCs) and nitrogen oxides (NOx) in Huzhou during August 23rd to September 23rd,2017,the pollution characteristics of ozone and its precursors were analyzed.Furthermore,source apportionment of VOCs and the formation mechanism of ozone were investigated using positive matrix factorization (PMF) and observation based model (OBM),respectively.The results showed that average volume fraction of VOCs in Huzhou was (24.78±9.10)×10-9 during the observation period.Alkanes,oxygenated VOCs (OVOCs),and halogenated hydrocarbons accounted for highest proportion to total VOCs.During the period when ozone concentration did not exceed air quality standard,ozone formation was in VOCs-limited regime,while during ozone pollution days,O3 formation was in transitional regime but still limited by VOCs.During O3 pollution days,aromatic hydrocarbons,olefins,and OVOCs were the three largest contributors to ozone formation potential (OFP),with relative contributions of 39.6%,21.5%,and 19.4%,respectively.The top three key components for OFP were toluene,ethylene,and m,p-xylene.PMF source apportionment results showed that the VOCs in Huzhou during the observation period were from solvent use (27.0%),traffic emission (22.7%),background and transport (19.3%),industrial emission (16.9%),gasoline volatilization (7.7%),and biogenic emission (6.4%).During O3 pollution episodes,traffic emission and solvent use were the two largest contributors to OFP,with relative contributions of 35.1% and 30.5%,respectively.Therefore,reduction in traffic emissions and solvent use is an effective way to control ozone pollution in Huzhou.

    Key words ozone;volatile organic compounds;sensitivity analysis;ozone formation potential;source apportionment

    收稿日期 2020-09-30

    作者簡(jiǎn)介常艷春,女,碩士,高級(jí)工程師,研究方向?yàn)榄h(huán)境監(jiān)測(cè).89996831@qq.com

    猜你喜歡
    湖州市貢獻(xiàn)臭氧
    湖州市庚村陽(yáng)桃規(guī)范化栽培技術(shù)
    湖州市吳興區(qū):構(gòu)建消費(fèi)幫扶機(jī)制 支援工作走深走實(shí)
    人民周刊(2022年24期)2022-02-07 01:30:54
    湖州市吳興區(qū):融入長(zhǎng)三角一體化發(fā)展
    人民周刊(2022年24期)2022-02-07 01:30:52
    文印室內(nèi)臭氧散發(fā)實(shí)測(cè)分析及模擬
    中國(guó)共產(chǎn)黨百年偉大貢獻(xiàn)
    為加快“三個(gè)努力建成”作出人大新貢獻(xiàn)
    My family
    貢獻(xiàn)榜
    海洋貢獻(xiàn)2500億
    商周刊(2017年6期)2017-08-22 03:42:37
    看不見(jiàn)的污染源——臭氧
    两人在一起打扑克的视频| 人人妻,人人澡人人爽秒播| 神马国产精品三级电影在线观看| 日韩欧美一区二区三区在线观看| 丝袜美腿在线中文| 色av中文字幕| 欧美绝顶高潮抽搐喷水| 2021天堂中文幕一二区在线观| 国产一区二区激情短视频| 久久久久国内视频| 男女那种视频在线观看| 麻豆久久精品国产亚洲av| 精品国产三级普通话版| 两个人视频免费观看高清| 国产美女午夜福利| 18+在线观看网站| 日本黄色视频三级网站网址| 男女床上黄色一级片免费看| 亚洲五月婷婷丁香| 午夜精品一区二区三区免费看| 99视频精品全部免费 在线| 国内揄拍国产精品人妻在线| 一区二区三区免费毛片| 午夜久久久久精精品| 婷婷精品国产亚洲av| 亚洲精品成人久久久久久| 香蕉丝袜av| 国内精品久久久久久久电影| 久久久久久久久中文| 日本免费一区二区三区高清不卡| 免费高清视频大片| 国产一区二区三区在线臀色熟女| 成人特级av手机在线观看| 一区二区三区国产精品乱码| 一本精品99久久精品77| 免费观看精品视频网站| 男女那种视频在线观看| 韩国av一区二区三区四区| 国产精品98久久久久久宅男小说| 国产一区二区三区视频了| 亚洲av成人不卡在线观看播放网| 日韩成人在线观看一区二区三区| 97超视频在线观看视频| 亚洲狠狠婷婷综合久久图片| 免费搜索国产男女视频| 在线观看日韩欧美| 久久草成人影院| 亚洲人成网站高清观看| 伊人久久精品亚洲午夜| 一区二区三区高清视频在线| 97人妻精品一区二区三区麻豆| 亚洲熟妇熟女久久| aaaaa片日本免费| 午夜日韩欧美国产| 亚洲国产精品999在线| 超碰av人人做人人爽久久 | 国产在视频线在精品| 亚洲专区国产一区二区| 久久精品影院6| 麻豆国产av国片精品| 天堂网av新在线| 亚洲 欧美 日韩 在线 免费| 亚洲精品国产精品久久久不卡| 一级作爱视频免费观看| 搡老熟女国产l中国老女人| 亚洲精品在线观看二区| 国产黄色小视频在线观看| 欧美中文日本在线观看视频| 三级男女做爰猛烈吃奶摸视频| 99久久精品热视频| 久9热在线精品视频| 啦啦啦观看免费观看视频高清| 不卡一级毛片| 国产一区二区在线av高清观看| 国产欧美日韩一区二区精品| 亚洲午夜理论影院| 精品国产三级普通话版| 国产欧美日韩一区二区精品| 成人欧美大片| 好男人电影高清在线观看| 欧美一区二区国产精品久久精品| 国产单亲对白刺激| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久久久久精品电影| 性色avwww在线观看| 九色国产91popny在线| 亚洲精品在线美女| 综合色av麻豆| 国产视频内射| 婷婷亚洲欧美| 国产伦精品一区二区三区视频9 | 国产精品国产高清国产av| 男人舔女人下体高潮全视频| 日韩高清综合在线| 99热精品在线国产| 最好的美女福利视频网| 欧美激情在线99| 久久伊人香网站| 国产老妇女一区| 深夜精品福利| 狂野欧美白嫩少妇大欣赏| 国产欧美日韩精品一区二区| 丁香欧美五月| 美女cb高潮喷水在线观看| 欧美最黄视频在线播放免费| 麻豆成人av在线观看| 人人妻人人澡欧美一区二区| 亚洲av第一区精品v没综合| 亚洲久久久久久中文字幕| 波野结衣二区三区在线 | 亚洲男人的天堂狠狠| 激情在线观看视频在线高清| 长腿黑丝高跟| 女人被狂操c到高潮| 又黄又爽又免费观看的视频| 97超视频在线观看视频| 亚洲一区高清亚洲精品| 日本成人三级电影网站| 久久久久久大精品| 亚洲男人的天堂狠狠| 中文亚洲av片在线观看爽| 国产高清视频在线观看网站| xxxwww97欧美| 国语自产精品视频在线第100页| 久久久久久久久大av| 男人和女人高潮做爰伦理| 99久久精品一区二区三区| av中文乱码字幕在线| 亚洲欧美一区二区三区黑人| 男女那种视频在线观看| 日韩欧美精品v在线| 欧美成人a在线观看| 真实男女啪啪啪动态图| 国产午夜精品久久久久久一区二区三区 | 琪琪午夜伦伦电影理论片6080| 亚洲人成伊人成综合网2020| 国产在线精品亚洲第一网站| 啦啦啦韩国在线观看视频| 丝袜美腿在线中文| 美女 人体艺术 gogo| 成人av在线播放网站| 叶爱在线成人免费视频播放| 很黄的视频免费| 一区二区三区国产精品乱码| 激情在线观看视频在线高清| 国产三级在线视频| 国产精品一及| a级毛片a级免费在线| 久久久色成人| 中文字幕久久专区| 日韩欧美在线二视频| 久久久久久久久中文| 麻豆成人av在线观看| 国内揄拍国产精品人妻在线| 国产精品久久久人人做人人爽| 狠狠狠狠99中文字幕| 在线观看免费午夜福利视频| 国产精品 欧美亚洲| 亚洲中文日韩欧美视频| 午夜福利高清视频| 变态另类丝袜制服| 中亚洲国语对白在线视频| 国产日本99.免费观看| 午夜日韩欧美国产| 两个人的视频大全免费| 亚洲18禁久久av| 岛国在线免费视频观看| 精品无人区乱码1区二区| 脱女人内裤的视频| 亚洲av中文字字幕乱码综合| 在线观看66精品国产| 精华霜和精华液先用哪个| a级毛片a级免费在线| 狂野欧美白嫩少妇大欣赏| 中文字幕av在线有码专区| 午夜福利在线在线| 有码 亚洲区| 亚洲人成电影免费在线| 夜夜躁狠狠躁天天躁| 国产精品久久久久久久电影 | 国产麻豆成人av免费视频| 99精品在免费线老司机午夜| 日韩中文字幕欧美一区二区| 久久精品国产清高在天天线| 欧美成人一区二区免费高清观看| 国产av一区在线观看免费| 免费无遮挡裸体视频| 男女做爰动态图高潮gif福利片| 精品一区二区三区视频在线 | 精品久久久久久久末码| 成人永久免费在线观看视频| 床上黄色一级片| 少妇的丰满在线观看| 国产视频一区二区在线看| 成人特级av手机在线观看| 亚洲人成网站在线播| 国内精品美女久久久久久| 日本免费a在线| 国产在视频线在精品| 欧美又色又爽又黄视频| 日本一二三区视频观看| 无限看片的www在线观看| 久久中文看片网| 精品无人区乱码1区二区| 在线观看66精品国产| 午夜福利免费观看在线| 真实男女啪啪啪动态图| av欧美777| 一级作爱视频免费观看| 国产视频一区二区在线看| 最近视频中文字幕2019在线8| 色尼玛亚洲综合影院| 日韩精品青青久久久久久| 久久亚洲真实| 亚洲国产日韩欧美精品在线观看 | 亚洲av中文字字幕乱码综合| 色综合亚洲欧美另类图片| 在线免费观看不下载黄p国产 | 伊人久久精品亚洲午夜| 床上黄色一级片| 蜜桃久久精品国产亚洲av| 在线观看av片永久免费下载| 老鸭窝网址在线观看| 麻豆一二三区av精品| 亚洲美女视频黄频| 欧美黄色淫秽网站| 91字幕亚洲| 午夜福利在线观看吧| 少妇的逼水好多| 12—13女人毛片做爰片一| 亚洲精品日韩av片在线观看 | 久久性视频一级片| 色尼玛亚洲综合影院| 网址你懂的国产日韩在线| avwww免费| 首页视频小说图片口味搜索| a级一级毛片免费在线观看| 午夜两性在线视频| 免费高清视频大片| 可以在线观看的亚洲视频| 亚洲午夜理论影院| 久久国产乱子伦精品免费另类| а√天堂www在线а√下载| 99热这里只有精品一区| 一个人看的www免费观看视频| 免费av毛片视频| 一本久久中文字幕| 国产av麻豆久久久久久久| 国产精品影院久久| 听说在线观看完整版免费高清| 色视频www国产| 亚洲国产精品合色在线| 亚洲美女视频黄频| 嫁个100分男人电影在线观看| 国产激情欧美一区二区| 亚洲美女黄片视频| 日本 欧美在线| 日韩大尺度精品在线看网址| 亚洲 国产 在线| 51午夜福利影视在线观看| 中文字幕久久专区| 欧美黄色片欧美黄色片| 亚洲国产精品久久男人天堂| 国产一区二区在线观看日韩 | 久久精品91蜜桃| 国产高清有码在线观看视频| 亚洲国产日韩欧美精品在线观看 | 久久精品亚洲精品国产色婷小说| 亚洲av日韩精品久久久久久密| av黄色大香蕉| 久久精品国产亚洲av香蕉五月| 成年免费大片在线观看| 婷婷精品国产亚洲av在线| 国产精华一区二区三区| 又紧又爽又黄一区二区| 欧美成人a在线观看| 国产伦在线观看视频一区| 国产黄a三级三级三级人| 国产久久久一区二区三区| 亚洲真实伦在线观看| 麻豆国产av国片精品| 亚洲专区国产一区二区| 18禁美女被吸乳视频| 午夜日韩欧美国产| 久久久国产精品麻豆| 国产成人av激情在线播放| 国产蜜桃级精品一区二区三区| 免费看日本二区| 亚洲精品亚洲一区二区| 精品熟女少妇八av免费久了| 999久久久精品免费观看国产| 天天躁日日操中文字幕| 国产毛片a区久久久久| 露出奶头的视频| 一级毛片女人18水好多| 高清日韩中文字幕在线| 天堂影院成人在线观看| 国产免费av片在线观看野外av| 亚洲第一欧美日韩一区二区三区| 12—13女人毛片做爰片一| 国产黄色小视频在线观看| 18禁裸乳无遮挡免费网站照片| 天天添夜夜摸| 国产精品嫩草影院av在线观看 | 亚洲av二区三区四区| 99精品欧美一区二区三区四区| 热99re8久久精品国产| 免费看光身美女| 欧美黑人欧美精品刺激| 在线免费观看不下载黄p国产 | 午夜激情福利司机影院| 日本熟妇午夜| 丰满的人妻完整版| 18禁黄网站禁片午夜丰满| www日本在线高清视频| 国产精品免费一区二区三区在线| 少妇人妻精品综合一区二区 | 国产伦在线观看视频一区| 国产三级黄色录像| 日日摸夜夜添夜夜添小说| 国产免费av片在线观看野外av| 亚洲av不卡在线观看| 日本与韩国留学比较| 狂野欧美白嫩少妇大欣赏| 国产免费一级a男人的天堂| 一本久久中文字幕| 精品一区二区三区视频在线观看免费| 亚洲精品影视一区二区三区av| 久久精品国产亚洲av涩爱 | 老司机福利观看| 在线视频色国产色| 欧美一级毛片孕妇| 久久精品国产自在天天线| 欧美3d第一页| 国产三级在线视频| 成年女人毛片免费观看观看9| 男女午夜视频在线观看| 欧美中文日本在线观看视频| 精华霜和精华液先用哪个| 女人被狂操c到高潮| 男女午夜视频在线观看| 国产日本99.免费观看| 男人舔奶头视频| 国产精品美女特级片免费视频播放器| 欧美一区二区精品小视频在线| 亚洲欧美日韩高清专用| ponron亚洲| 在线观看日韩欧美| 国内少妇人妻偷人精品xxx网站| 久久久国产成人免费| 亚洲一区二区三区色噜噜| 国产亚洲精品久久久com| 一区二区三区激情视频| 午夜老司机福利剧场| 日本成人三级电影网站| 男插女下体视频免费在线播放| 亚洲人成网站在线播放欧美日韩| 亚洲最大成人手机在线| 悠悠久久av| 国产欧美日韩一区二区精品| 国产乱人视频| xxxwww97欧美| 亚洲精品成人久久久久久| 婷婷丁香在线五月| 国产主播在线观看一区二区| 尤物成人国产欧美一区二区三区| 亚洲av美国av| 成人一区二区视频在线观看| 亚洲国产欧美人成| 91在线观看av| 国产黄色小视频在线观看| 免费看光身美女| 69人妻影院| 国产探花极品一区二区| 99久久无色码亚洲精品果冻| 国产精品亚洲一级av第二区| 国产精品嫩草影院av在线观看 | 99久久精品一区二区三区| 国产一区二区激情短视频| 成人精品一区二区免费| 超碰av人人做人人爽久久 | 天堂动漫精品| 欧美一级毛片孕妇| 亚洲欧美日韩东京热| 中文字幕高清在线视频| 日本免费一区二区三区高清不卡| 亚洲av电影在线进入| svipshipincom国产片| 亚洲av成人不卡在线观看播放网| 高清日韩中文字幕在线| 中文字幕精品亚洲无线码一区| 在线观看免费视频日本深夜| 一个人免费在线观看电影| 美女被艹到高潮喷水动态| 国产蜜桃级精品一区二区三区| 国产69精品久久久久777片| 三级国产精品欧美在线观看| 母亲3免费完整高清在线观看| 中文资源天堂在线| 亚洲国产欧美人成| 午夜久久久久精精品| 国产午夜精品久久久久久一区二区三区 | 精品久久久久久久毛片微露脸| 中文字幕av在线有码专区| 噜噜噜噜噜久久久久久91| 91九色精品人成在线观看| 18禁在线播放成人免费| 亚洲国产精品合色在线| 性色avwww在线观看| 不卡一级毛片| x7x7x7水蜜桃| 午夜福利18| 人妻久久中文字幕网| 久久精品国产亚洲av香蕉五月| 国产 一区 欧美 日韩| 日本精品一区二区三区蜜桃| 亚洲av中文字字幕乱码综合| 日韩欧美免费精品| 亚洲av五月六月丁香网| 免费观看的影片在线观看| 免费在线观看影片大全网站| 国内精品久久久久久久电影| 欧美av亚洲av综合av国产av| 在线播放国产精品三级| 哪里可以看免费的av片| 美女被艹到高潮喷水动态| 久久久精品欧美日韩精品| 久久精品国产亚洲av涩爱 | 丰满人妻熟妇乱又伦精品不卡| 黄片小视频在线播放| 国产探花在线观看一区二区| 99国产综合亚洲精品| 搞女人的毛片| 中国美女看黄片| 国产一级毛片七仙女欲春2| 别揉我奶头~嗯~啊~动态视频| 久久精品人妻少妇| 99久国产av精品| 日韩中文字幕欧美一区二区| 成人无遮挡网站| 国产99白浆流出| 无人区码免费观看不卡| 18禁黄网站禁片免费观看直播| 亚洲男人的天堂狠狠| 美女黄网站色视频| 中文字幕高清在线视频| 国产精品免费一区二区三区在线| svipshipincom国产片| 91字幕亚洲| 少妇丰满av| 亚洲黑人精品在线| 韩国av一区二区三区四区| www.999成人在线观看| 尤物成人国产欧美一区二区三区| 亚洲美女黄片视频| 无人区码免费观看不卡| 精品福利观看| 男女之事视频高清在线观看| 亚洲中文日韩欧美视频| 亚洲精品国产精品久久久不卡| 亚洲不卡免费看| 婷婷六月久久综合丁香| 亚洲国产精品久久男人天堂| 日韩 欧美 亚洲 中文字幕| 美女高潮的动态| 伊人久久大香线蕉亚洲五| 高清毛片免费观看视频网站| 欧美日韩瑟瑟在线播放| 国产伦人伦偷精品视频| 亚洲中文日韩欧美视频| 亚洲av成人av| 蜜桃久久精品国产亚洲av| 免费看a级黄色片| 老司机午夜十八禁免费视频| 91字幕亚洲| 无人区码免费观看不卡| 又粗又爽又猛毛片免费看| 男插女下体视频免费在线播放| 欧美日韩综合久久久久久 | 无限看片的www在线观看| 最近最新免费中文字幕在线| 精品人妻偷拍中文字幕| 欧美极品一区二区三区四区| 美女黄网站色视频| 91九色精品人成在线观看| 天天躁日日操中文字幕| 久久久久久九九精品二区国产| 又爽又黄无遮挡网站| 久久精品国产清高在天天线| 亚洲av中文字字幕乱码综合| 亚洲aⅴ乱码一区二区在线播放| 国产极品精品免费视频能看的| 国产日本99.免费观看| 搡女人真爽免费视频火全软件 | 亚洲av美国av| 欧美日韩国产亚洲二区| 成人无遮挡网站| 亚洲 国产 在线| 亚洲人成网站在线播| 老汉色av国产亚洲站长工具| 亚洲最大成人中文| 真人一进一出gif抽搐免费| 色在线成人网| 国产在线精品亚洲第一网站| 超碰av人人做人人爽久久 | 久久久久久久精品吃奶| 12—13女人毛片做爰片一| 老司机午夜福利在线观看视频| 宅男免费午夜| 国产午夜精品论理片| www.www免费av| av天堂在线播放| 亚洲精品乱码久久久v下载方式 | 欧美不卡视频在线免费观看| 久久午夜亚洲精品久久| 亚洲成人中文字幕在线播放| 怎么达到女性高潮| 亚洲av熟女| 国产高潮美女av| 老司机福利观看| 长腿黑丝高跟| 婷婷亚洲欧美| 亚洲成人久久爱视频| 一区福利在线观看| 午夜视频国产福利| 日本 欧美在线| 好男人电影高清在线观看| 亚洲国产中文字幕在线视频| 国产精品日韩av在线免费观看| 少妇高潮的动态图| 一进一出抽搐动态| 国产蜜桃级精品一区二区三区| 中文字幕av在线有码专区| 最近最新免费中文字幕在线| 国产精品影院久久| 精品人妻一区二区三区麻豆 | 国产激情欧美一区二区| 亚洲黑人精品在线| 伊人久久大香线蕉亚洲五| 欧美国产日韩亚洲一区| 九九热线精品视视频播放| 在线观看66精品国产| aaaaa片日本免费| 久久中文看片网| 久久精品亚洲精品国产色婷小说| 午夜精品在线福利| 五月伊人婷婷丁香| 国产激情偷乱视频一区二区| 色精品久久人妻99蜜桃| 免费看a级黄色片| 观看免费一级毛片| 亚洲av第一区精品v没综合| 每晚都被弄得嗷嗷叫到高潮| 欧美极品一区二区三区四区| 欧美日韩中文字幕国产精品一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲人成伊人成综合网2020| 3wmmmm亚洲av在线观看| 色尼玛亚洲综合影院| 国产精品,欧美在线| 母亲3免费完整高清在线观看| 高潮久久久久久久久久久不卡| www.999成人在线观看| 精品一区二区三区av网在线观看| 五月玫瑰六月丁香| 99国产精品一区二区三区| 在线观看日韩欧美| 日韩欧美一区二区三区在线观看| 看黄色毛片网站| 午夜视频国产福利| 男女之事视频高清在线观看| 成人av一区二区三区在线看| 亚洲av不卡在线观看| 国产精品一及| 日日夜夜操网爽| 亚洲欧美一区二区三区黑人| 在线观看av片永久免费下载| 天堂影院成人在线观看| 免费观看人在逋| 三级国产精品欧美在线观看| 欧美中文综合在线视频| av在线天堂中文字幕| 亚洲中文字幕一区二区三区有码在线看| 老司机福利观看| 国产精品爽爽va在线观看网站| 成人特级黄色片久久久久久久| 亚洲欧美精品综合久久99| 欧美一区二区国产精品久久精品| 亚洲人成网站在线播放欧美日韩| 老司机午夜福利在线观看视频| 国产三级黄色录像| 久久香蕉精品热| 99久国产av精品| 久久久国产精品麻豆| 亚洲 欧美 日韩 在线 免费| 欧美另类亚洲清纯唯美| 精品久久久久久久毛片微露脸| 国产99白浆流出| 法律面前人人平等表现在哪些方面| 极品教师在线免费播放| 1000部很黄的大片| 高清在线国产一区| 深爱激情五月婷婷| ponron亚洲| 久久人妻av系列| 欧美激情在线99| 免费看a级黄色片| 窝窝影院91人妻| a在线观看视频网站| 中文字幕人成人乱码亚洲影| 中文字幕高清在线视频| 欧美成狂野欧美在线观看| 久久这里只有精品中国| 国产成+人综合+亚洲专区| 女警被强在线播放|