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

    武漢軍運會前后臭氧及其前體物的特征和來源

    2022-07-19 01:11:04燕瑩瑩孔少飛宋藹莉全繼宏祁士華
    中國環(huán)境科學 2022年7期
    關(guān)鍵詞:貢獻率機動車臭氧

    馬 靜,燕瑩瑩*,孔少飛,2,宋藹莉,陳 楠,祝 波,全繼宏,祁士華,2

    武漢軍運會前后臭氧及其前體物的特征和來源

    馬 靜1,燕瑩瑩1*,孔少飛1,2,宋藹莉1,陳 楠2,3,祝 波2,3,全繼宏2,3,祁士華1,2

    (1.中國地質(zhì)大學(武漢)環(huán)境學院,湖北 武漢 430074;2.湖北省大氣復合污染研究中心,湖北 武漢 430074;3.湖北省生態(tài)環(huán)境監(jiān)測中心站,湖北 武漢 430072)

    利用湖北省超級站2019年10~11月的臭氧、NO(= NO + NO2)和102種VOCs物質(zhì)的小時數(shù)據(jù)分析了軍運會期間臭氧污染變化;基于DSMACC箱型模式模擬不同VOCs和NO濃度下臭氧的光化學生成敏感性;采用PMF模型對前體物VOCs進行源解析,并估算不同源類的臭氧生成潛勢.結(jié)果顯示,軍運會保障前臭氧日最大8小時濃度(最大MDA8:219.51mg/m3)超過國家二級標準,保障期臭氧MDA8濃度(135.11mg/m3)明顯下降,保障后濃度回升(140.98mg/m3).軍運會保障前中期臭氧濃度的差異受氣象條件影響更明顯,而保障后臭氧濃度的上升主要是因為前體物濃度的大幅增加.根據(jù)DSMACC模擬的EKMA曲線,武漢市軍運會期間臭氧的光化學生成主要受VOCs濃度變化的影響.進一步對VOCs進行源解析,結(jié)果顯示,保障前VOCs對臭氧生成貢獻較大的源類是燃燒源、石油化工和機動車,分別占23.0%、22.8%和22.5%;保障期間VOCs的主要來源是機動車(38.4%)和燃燒源(25.5%);保障后則主要是石油化工(32.6%)和燃料揮發(fā)(25.7%).三個階段對比發(fā)現(xiàn),軍運會的保障方案對石油化工源減排效果明顯,但對機動車和燃燒源排放的限制效果并不顯著.武漢市應該更注重對燃燒、燃料揮發(fā)和機動車排放的治理.

    臭氧污染;軍運會;源解析;EKMA曲線

    當前我國重點區(qū)域面臨臭氧和細顆粒物復合污染問題,其中臭氧污染日益突出[1].近年來臭氧濃度逐年上升[2-7],以臭氧為首要污染物的污染天氣占比逐漸增加[8],臭氧繼PM2.5之后備受關(guān)注[9].近地層臭氧濃度過高不僅會對農(nóng)作物的生長[10-12]和全球氣候變化[13]造成一定的影響,還會對人類身體健康產(chǎn)生威脅[14-15].我國城市近地層臭氧主要是由氮氧化物(NO)與揮發(fā)性有機化合物(VOCs)以及一氧化碳(CO)等多種前體物在光照條件下發(fā)生一系列光化學氧化反應產(chǎn)生[16],并通過化學消耗和干沉降過程從大氣中去除.

    局地臭氧光化學過程與主要前體物(VOCs和NO)的關(guān)系受到廣泛關(guān)注[17-24].Tang等[25]發(fā)現(xiàn),華北地區(qū)夏季近地面大氣中的臭氧光化學生成由1997年的受NO控制轉(zhuǎn)變成2010年的受VOCs控制.Ding等[26]研究表明,長三角地區(qū)臭氧光化學生成受VOCs控制.對珠江三角洲地區(qū)臭氧污染做敏感性分析,Wang等[27]研究表明珠三角中部內(nèi)陸地區(qū)和周邊沿海地區(qū)的臭氧生成受VOCs控制,而珠三角西南部農(nóng)村地區(qū)的臭氧生成受NO控制.這些研究表明,我國城市地區(qū)臭氧的光化學生成對VOCs高度敏感.因而,VOCs來源解析以及對臭氧生成的貢獻成為近年來研究的熱點問題.已有研究表明,我國東部地區(qū)汽車尾氣排放和生物燃料燃燒占71%,汽油泄漏和溶劑使用占7%,生物質(zhì)燃燒和工業(yè)分別占11%[28];上海地區(qū)汽車尾氣是VOCs最大的排放源[29].

    過去,我國臭氧污染的研究主要集中在四個區(qū)域:京津冀[30]、長三角[31]、珠三角[32]和成渝地區(qū)[33-34].然而,我國中部地區(qū)夏季臭氧污染頻發(fā),形成一個新的區(qū)域性大氣污染核心區(qū)[35-37].作為長江流域經(jīng)濟發(fā)展的中心,武漢及周邊城市群船舶、交通、工業(yè)和農(nóng)業(yè)活動等排放密集使得臭氧前體物污染源復雜,且具有高溫高濕的氣象條件,易于臭氧污染事件形成.近年來武漢地區(qū)臭氧污染呈明顯增加趨勢[38],該地區(qū)臭氧污染及其前體物的特征和來源值得關(guān)注.

    2019年10月18日至27日第七屆世界軍人運動會(簡稱軍運會)在武漢舉行,這次軍運會是世界軍人運動會歷史上規(guī)模最大、參賽人員最多的一次運動會.根據(jù)軍運會環(huán)境空氣質(zhì)量保障工作要求,在生態(tài)環(huán)境部的統(tǒng)一部署下,湖北、河南、安徽等省份確立了保障方案,對武漢市各企業(yè)、各類工地施工和船舶、飛機、燃油非道路移動機械等大氣污染物排放源,嚴格執(zhí)行大氣環(huán)境質(zhì)量管理臨時性措施.保障方案明確會期保障時段為10月8日至10月28日,保障目標為確保會期武漢市環(huán)境空氣質(zhì)量達到優(yōu)良水平.

    類似武漢軍運會這樣大型活動的環(huán)境狀況也有許多學者做過研究.例如,孫天樂等[39]對2008年殘奧會期間黑炭氣溶膠的研究表明,殘奧會期間本地黑炭排放有所降低.王羽琴等[40]研究表明深圳大運會期間PM2.5得到很好的控制,但污染控制措施取消后PM2.5濃度突增.在大型活動期間,采取相應的措施能有效地改善大氣環(huán)境.然而,這些研究大多關(guān)注細顆粒物,少有對臭氧的研究.

    本研究分析對比了軍運會保障期間與保障前后湖北超級站臭氧及其前體物(NO和VOCs)濃度的變化.利用DSMACC箱型模式模擬了不同VOCs和NO濃度下臭氧的生成情況,確定了臭氧光化學生成的敏感性.進一步采用PMF模型對VOCs進行源解析,計算臭氧生成潛勢.本研究的主要創(chuàng)新點為揭示了中部地區(qū)特大城市武漢臭氧污染的主要來源,評估了軍運會保障期間武漢臭氧污染的改善情況及污染減排貢獻.研究結(jié)果可為武漢及周邊地區(qū)臭氧污染的預報預警、臭氧污染成因識別和污染防控政策完善提供科學依據(jù).

    1 數(shù)據(jù)與方法

    1.1 數(shù)據(jù)來源

    本文所使用的觀測數(shù)據(jù)來自武漢市大氣復合污染實驗室(“超級站”),超級站位于湖北省武漢市江漢區(qū)香江路1號,具體地理坐標為(114.2777°E, 30.6042°N).該站點地處漢口中心城區(qū),為二類功能區(qū),觀測點距地面約25m,周圍為居民區(qū),無建筑物遮擋,視線較為開闊,無明顯工業(yè)污染源,監(jiān)測數(shù)據(jù)一定程度上可以代表武漢城市中心城區(qū)大氣污染狀況水平.分析的污染物有臭氧、NO、NO2和102種VOCs物質(zhì),時間段為2019年10月1日~11月7日.本文中,10月8~28日是軍運會保障期,10月1~8日為保障前,10月29日~11月7日為保障后.10月1~7日國慶期間,城市人口、車輛減少,大氣環(huán)境質(zhì)量較好[41].

    1.2 DSMACC模擬

    DSMACC(The Dynamical Simple Model of Atmospheric Chemical Complexity)是一種常用的箱型模式(https://github.com/barronh/DSMACC),可用于研究臭氧及其前體物之間的聯(lián)系.

    采用DSMACC模型模擬武漢地區(qū)不同VOCs和NO濃度下臭氧的光化學生成情況,以表征O3-VOC-NO反應系統(tǒng)的非線性特征.時間設(shè)置為中午(12:00~13:00,臭氧生成效率和臭氧濃度在一天中都是最高的時段[42]),天氣條件設(shè)定為晴空.氣象輸入和污染物初始濃度根據(jù)觀測的平均值設(shè)置.

    1.3 PMF模型

    正矩陣分解(PMF)模型是美國環(huán)保署(USEPA)官方推薦的源解析技術(shù)方法,它被廣泛應用于VOCs源解析[43-44].它可以對數(shù)據(jù)產(chǎn)生更好的擬合,更好的描述污染物的來源[45-46].PMF模型將受體樣本矩陣(×)分解為貢獻矩陣(×)和因子矩陣(×),并使用多次迭代算法使目標函數(shù)Q最小化,進而確定最優(yōu)解析結(jié)果,計算公式如下:

    式中:x表示第個樣本中第個物質(zhì)的質(zhì)量濃度,單位mg/m3;g表示第個來源對第個樣本的貢獻率;f表示第個排放源中第個物質(zhì)的質(zhì)量分數(shù);為排放源的種類數(shù);為需要解析的因子數(shù);e表示第個樣本中第個物質(zhì)的殘差;和分別表示樣本數(shù)和被測物質(zhì)數(shù);u表示第個樣本中第個物質(zhì)的不確定度;Q是基于樣品監(jiān)測誤差的目標函數(shù).

    PMF通過最小二乘法來確定污染源類型及其貢獻,并且對污染源類型及其貢獻具有非負的限制.一般由誤差比例和方法檢出限來確定物種的不確定度.通常物質(zhì)的濃度低于或等于檢出限時,用公式(3)計算;物種的濃度高于檢出限時,用公式(4)計算:

    式中:U為物種的不確定性;MDL為物種的檢出限;為物種的質(zhì)量濃度,mg/m3;EF為誤差分數(shù),本研究中取20%.

    本文中,并非所有被測量的物種都輸入到PMF模型中,而是按照以下方法來選擇適當?shù)腣OCs物種:(1)保留高相關(guān)性物質(zhì);(2)保留能表征揮發(fā)性有機化合物來源的物質(zhì);(3)保留高于檢出限50%的物質(zhì)[47-48].

    1.4 臭氧生成潛勢

    臭氧生成潛勢(OFP)是衡量VOCs的反應活性及對臭氧影響的重要指標,能分析形成臭氧的優(yōu)勢成分.其計算公式為:

    OFP= VOC× MIR(5)

    式中:OFP是某個VOCs物質(zhì)的臭氧生成潛勢; VOC是該VOCs物質(zhì)的濃度;MIR是該VOCs物質(zhì)的反應性,數(shù)值選擇參考前人研究[49].

    本文計算源解析后各因子的臭氧生成潛勢.其計算公式[43]為:

    式中:g為第個源在第個樣本中的貢獻率;f為第個源中第種的濃度.

    2 結(jié)果與討論

    2.1 軍運會期間臭氧及其前體物的變化特征

    2.1.1 武漢市臭氧變化情況武漢地區(qū)臭氧濃度具有明顯的時間變化規(guī)律,年際變化呈現(xiàn)上升的趨勢;在年內(nèi)變化上呈現(xiàn)春夏季高于秋冬的特點[4].由圖1(a)和圖2(a)可知,軍運會保障前臭氧MDA8濃度超標,最大濃度為219.51mg/m3;保障期臭氧MDA8濃度明顯下降,最大濃度為135.11mg/m3;保障后濃度反彈,最大濃度為140.98mg/m3.保障期較之前的臭氧平均濃度下降34.3%,保障后較保障期臭氧平均濃度上升17.0%.

    圖1 (a)軍運會期間臭氧MDA8濃度及其對應時間NOx濃度和TVOC濃度日變化;(b)研究時段溫度、相對濕度、風速和輻射日變化

    2.1.2 武漢市氮氧化物變化情況由圖1(a)和圖2(a)可知,保障期NO和NO2濃度相比于保障前變化不大,NO濃度不超過5.48mg/m3,NO2濃度不超過48.32mg/m3.保障后二者濃度有明顯的升高趨勢,最大濃度分別達到14.14mg/m3和89.45mg/m3,保障后二者總濃度較保障前上升43.2%,保障后NO和NO2出現(xiàn)嚴重反彈,這可能與飛機等非道路排放源和工業(yè)生產(chǎn)有關(guān).軍運會后解除了對非道路排放源和工廠的臨時措施(http://hbj.wuhan.gov.cn/fbjd_19/zc/ gfxwj/202001/t20200106_560629.html).

    圖2 (a)軍運會保障前中后臭氧、NO、NO2和TVOC濃度的平均值,綠色線段表示標準偏差;(b)軍運會保障前中后的溫度、相對濕度、風速和輻射的平均值

    2.1.3 武漢市TVOC變化情況由圖1(a)和圖2(a)可知,保障前和保障期VOCs日均濃度變化波動較小,平均值分別為79.70mg/m3和84.49mg/m3.保障后TVOC濃度明顯升高,最大濃度達到了155.19mg/m3,平均濃度較保障前上升了47.1%,達到117.23mg/m3.

    表1列舉了14種標識性的VOCs物質(zhì),計算了不同保障時期每種物質(zhì)的最大和平均濃度.結(jié)果顯示,除異戊二烯總體濃度變化不大外,其余各VOCs物種(烷烴、烯烴和芳香烴類)在保障前中期濃度差異較小,保障后濃度均有不同程度的升高.

    2.1.4 武漢市天氣條件變化情況由圖1(b)和圖2(b)可知,保障前溫度明顯偏高,最高為34.0℃,平均溫度為24.7℃.保障期除21號溫度較高達到27.6℃外,其余天溫度都較低,平均值為19.4℃.保障后溫度較保障期差異較小,最高為25.6℃,平均值為20.0℃.保障前中后相對濕度略有差異,平均值分別為68.4%、72.2%和58.6%;同時,風速都較小(0.70~2.75m/s).輻射的變化波動較大,1號和19號出現(xiàn)最大值,分別為28.21W/m2和25.81W/m2,保障期輻射平均值為8.55W/m2,較保障前(9.52W/m2)和保障后(9.88W/m2)略小.

    表1 保障前中后部分VOCs的濃度

    2.1.5 臭氧與其前體物和天氣條件的相關(guān)分析由表1可知,保障前臭氧及其前體物相關(guān)性不高,相關(guān)系數(shù)均小于0.10,且都未通過90%的顯著性檢驗,而溫度、相對濕度和輻射與臭氧的相關(guān)性高,相關(guān)系數(shù)都大于0.50,且都通過99%的顯著性檢驗.

    表2 臭氧與其前體物和氣象條件的相關(guān)性及顯著水平

    注:為統(tǒng)計量,99%、95%和90%的T臨界值分別為2.3,1.7,1.3.

    保障期臭氧與前體物的相關(guān)性雖明顯大于保障前,但相關(guān)系數(shù)都較小(<0.4),而臭氧與溫度和相對濕度的相關(guān)性明顯減小,相關(guān)系數(shù)都減小0.20以上.保障后與保障期比,臭氧與其前體物的相關(guān)性更大,相關(guān)系數(shù)都在0.65以上,且都通過99%的顯著性檢驗,而臭氧與溫度、相對濕度和輻射相關(guān)性比保障前小,相關(guān)系數(shù)都減小0.2以上.相較于保障后期,保障前中期臭氧濃度變化受天氣條件影響更強烈,保障前溫度高,輻射強,相對濕度小,有利于臭氧生成.而保障后臭氧增加主要受氮氧化物和VOCs濃度回升的影響,天氣條件與保障前差異不大.

    2.2 武漢市臭氧光化學生成的敏感性

    圖3顯示了DSMACC模擬的臭氧光化學生成敏感性,其中等值線表示不同前體物(TVOC和NO)濃度下臭氧的生成情況,黑色脊線為最大臭氧形成值的連線,TVOC/NO比率為7.1.黑色圓圈為各階段白天(8:00~20:00)觀測的NO和TVOC小時平均濃度.

    圖3 (a)~(c)DSMACC模擬的軍運會保障前中后的臭氧等值線圖;(d)~(f)軍運會保障前中后不同時期的TVOC/NOx比值箱型圖

    由圖3可知,保障前TVOC/NO比率在0.6~4.7范圍內(nèi),保障期TVOC/NO比率在0.2~4.9范圍內(nèi),保障后TVOC/NO比率在0.4~5.8范圍內(nèi).保障前中后臭氧生成都處于VOCs控制區(qū),武漢市的臭氧形成主要受VOCs控制.三個時期臭氧的化學生成區(qū)間差異不大,說明軍運會的管控措施對NO和VOCs的管控力度相差較小,達到了協(xié)同減排的效果.

    2.3 VOCs源解析及對臭氧生成的貢獻

    2.3.1 VOCs來源根據(jù)EPA PMF5.0對VOCs源解析的計算結(jié)果,結(jié)合文獻中典型VOCs源類的標識成分,最終把污染源歸為六類(圖4).

    因子1主要貢獻成分是異戊二烯,異戊二烯是植物光合作用排放的物質(zhì)[50],因此因子1歸為生物源.日變化顯示中午2點出現(xiàn)最大峰值,符合植物源排放的日內(nèi)變化.

    因子2主要貢獻成分是丙烷,其次是鄰二甲苯、間/對二甲苯和乙苯等苯類物質(zhì),還有部分高碳烷烴,其中丙烷是天然氣的工業(yè)加工產(chǎn)物[51],鄰二甲苯、間/對二甲苯來自汽車尾氣排放[43],乙苯主要用于生產(chǎn)其它化學產(chǎn)品,而高碳烷烴主要來自工業(yè)生產(chǎn)[48],這里把因子2歸為石油化工.日變化顯示夜間污染物濃度高于白天,符合工廠夜間排放高于白天.

    因子3主要貢獻成分是乙炔、1,2-二氯乙烷、1,2-二氯丙烷等,1,2-二氯乙烷、1,2-二氯丙烷為機動車尾氣排放[52],乙炔受機動車尾氣、油氣揮發(fā)、溶劑使用和LPG燃料等影響明顯[53],因此因子3歸為機動車尾氣.日變化中白天的值高于夜間值對應日間機動車使用高于夜間.

    因子4主要貢獻成分是1,2,4-三甲基苯、1,3,5-三甲基苯、1,2,3-三甲基苯、1,3-二乙基苯、1,4-二乙基苯,這些物質(zhì)都與溶劑的使用有關(guān)[48][54],因此把因子4歸為溶劑使用.溶劑的排放與工業(yè)密切相關(guān),在夜間出現(xiàn)峰值.

    因子5主要貢獻成分是乙烯、丙烯、1-丁烯,烯烴的主要來源為機動車排放、煤與生物質(zhì)燃燒、石化工業(yè)排放[55],其中丙烯是燃燒過程中排放的主要物質(zhì)[56],因此把因子5歸為燃燒源.日變化對應早上8點和晚上22點有峰值,最低值出現(xiàn)在下午16點.

    因子6主要貢獻成分主要是甲基環(huán)戊烷和大量的C4~C6烷烴,甲基戊烷既是工業(yè)合成的原料,也來自于汽油車尾氣[57-58],C5烷烴是汽油揮發(fā)的重要組成物質(zhì)[58],本文把因子6歸為燃料揮發(fā).日變化顯示夜間高于白天,這符合11:00后由于光化學反應的消耗濃度降低,到18:00光化學反應減弱,濃度升高.

    圖4 PMF源解析的結(jié)果,(a)~(f)為40種VOCs物質(zhì)對不同因子的貢獻率;(g)~(l)為各因子濃度的日內(nèi)變化

    如圖5(a)~(c)所示,保障前對武漢市VOCs貢獻最大的源是機動車排放,貢獻率為29.7% (11.26mg/ m3),其次是石油化工22.1%(8.37mg/m3),燃料揮發(fā)、燃燒源、溶劑和生物源分別為13.7%、11.7%、11.6%和11.1%(5.18、4.44、4.39和4.21mg/m3).保障期較保障前機動車排放(40.0%)和燃燒源(17.1%)貢獻率有所增加,其中機動車貢獻率增加最為顯著,增加10.3%,燃燒源增加5.4%,而石油化工(17.6%)和燃料揮發(fā)(7.3%)分別減少4.5%和6.4%,溶劑使用和生物源減少不明顯,其貢獻率分別為8.4%和9.6%(3.72和4.24mg/m3).保障后石油化工的貢獻率最大(30.4%),其次是燃料揮發(fā)(25.3%)和機動車排放(25.2%),溶劑、燃燒源和生物源貢獻率分別為7.1%、6.5%和5.6%(4.94、4.57和3.89mg/m3),較保障期石油化工和燃料揮發(fā)貢獻率顯著增加,機動車貢獻率雖有明顯下降,但總體占比依然較大.

    圖5 (a)~(c)軍運會保障前中后,各類源對VOCs濃度的貢獻率;(d)~(f)各類源的臭氧生成潛勢

    2.3.2 臭氧生成潛勢如圖5(d)~(f)所示,保障前燃燒源對OFP貢獻最大,貢獻率為23.0%(1.02mg/m3),其次是石油化工22.8%(1.01mg/m3)和機動車排放22.5%(1.00mg/m3),三者的貢獻率相差不大,燃料揮發(fā)、溶劑和生物源貢獻率分別為11.5%、10.7%和9.5%(0.51、0.48和0.42mg/m3).保障期機動車的貢獻率最大,為38.4%(2.16mg/m3),較保障前增加了15.9%,而石油化工為12.7%(0.72mg/m3),較保障前貢獻率降低了10.1%,溶劑和燃料揮發(fā)的貢獻率略有減小,燃燒源和生物源的貢獻率變化不大.保障后石油化工和燃料揮發(fā)的貢獻率顯著增大,分別達32.6%和25.7%(3.81和3.00mg/m3),燃燒源的貢獻率顯著下降,下降到14.1%.總的來說,機動車對臭氧生成的貢獻沒有明顯的減少,保障前中后溶劑使用對OPF的貢獻率也沒有明顯減小的趨勢.武漢軍運會對燃料揮發(fā)和石油化工管控效果明顯,對燃燒源和機動車排放的管控效果較弱,但保障后石油化工和燃料揮發(fā)反彈嚴重.

    根據(jù)前人對2010年上海世博會前后污染物濃度變化和來源解析研究表明[59-61],世博會期間氮氧化物得到了有效的控制,但后期由于天氣和治理力度不足等原因出現(xiàn)反彈,且治理方案對臭氧的減排效果不明顯.對來源的解析表明,汽車尾氣是重點需要關(guān)注的.這與本研究的結(jié)果類似,需要制定長期的管控方案,加強對機動車的管控.

    3 結(jié)論

    3.1 軍運會保障前臭氧日最大8h濃度(最大MDA8: 219.51mg/m3)超過國家二級標準,保障期臭氧MDA8濃度(135.11mg/m3)明顯下降,保障后濃度回升(140.98mg/m3).軍運會保障前和保障期臭氧濃度的差異受氣象條件影響更明顯,而保障后臭氧濃度的上升主要是因為前體物濃度的大幅增加.

    3.2 根據(jù)DSMACC模擬的EKMA曲線,武漢市軍運會期間TVOC/NO比值都在7.1以下,屬于VOCs控制區(qū),臭氧的光化學生成主要受VOCs濃度變化的影響.

    3.3 VOCs源解析結(jié)果顯示,保障前VOCs對臭氧生成貢獻較大的源類是燃燒源、石油化工和機動車,分別占23.0%、22.8%和22.5%;保障期間VOCs的主要來源是機動車(38.4%)和燃燒源(25.5%);保障后則主要是石油化工(32.6%)和燃料揮發(fā)(25.7%).三個階段對比發(fā)現(xiàn),軍運會的保障方案對石油化工源減排效果明顯,但對機動車和燃燒源排放的限制效果并不明顯.

    [1] Kan H, Chen R, Tong S. Ambient air pollution, climate change, and population health in China [J]. Environment International, 2012,42: 10-19.

    [2] Cheng L, Wang S, Gong Z, et al. Regionalization based on spatial and seasonal variation in ground-level ozone concentrations across China [J]. Journal of environmental sciences (China), 2018,67(5):179-190.

    [3] Liu H, Liu J, Liu Y, et al. Spatiotemporal variability and driving factors of ground-level summertime ozone pollution over eastern China [J]. Atmospheric Environment, 2021,265:118686.

    [4] 楊 喆,趙錦慧,劉玉青,等.武漢市臭氧時空分布及環(huán)境變量影響分析[J]. 湖北大學學報(自然科學版), 2021,43(5):522-528.

    Yang Z, Zhao J, Liu Y, et al. Analysis of the temporal and spatial distribution of ozone and the influence of environmental variables in Wuhan [J]. Journal of Hubei University (Natural Science), 2021,43(5): 522-528.

    [5] Li K, Jacob D J, Liao H, et al. Anthropogenic drivers of 2013~2017 trends in summer surface ozone in China [J]. Proceedings of the National Academy of Sciences of the United States of America, 2019,116(2):422-427.

    [6] Wang Y, Gao W, Wang S, et al. Contrasting trends of PM2.5and surface-ozone concentrations in China from 2013 to 2017 [J]. National Science Review, 2020,7(8):1331-1339.

    [7] Chen S, Wang H, Lu K, et al. The trend of surface ozone in Beijing from 2013 to 2019: Indications of the persisting strong atmospheric oxidation capacity [J]. Atmospheric Environment, 2020,242:117801.

    [8] 嚴 剛,薛文博,雷 宇,等.我國臭氧污染形勢分析及防控對策建議[J]. 環(huán)境保護, 2020,48(15):15-19.

    Yan G, Xue W, Lei Y, et al. Situation and control measures of ozone pollution in China [J]. Environmental Protection, 2020,48(15):15-19.

    [9] Anon. Cleaner air for China [J]. Nature Geoscience, 2019,12(7),doi:10. 1038/s41561-019-0406-7.

    [10] Tian H, Ren W, Tao B, et al. Climate extremes and ozone pollution: a growing threat to China's food security [J]. Ecosystem Health and Sustainability, 2016,2(1),e01203.

    [11] Tai A P K, Matin V M, Heald L C. Threat to future global food security from climate change and ozone air pollution [J]. Nature Climate Change, 2014,4(9):817-821.

    [12] Emberson L. Effects of ozone on agriculture, forests and grasslands [J]. Philosophical Transactions of the Royal Society A, 2020,378(2183): 20190327.

    [13] Monks P S, Archibald A T, Colette A, et al. Tropospheric ozone and its precursors from the urban to the global scale from air quality to short-lived climate forcer [J]. Atmospheric Chemistry and Physics, 2015,15(15):8889-8973.

    [14] Anenberg S C, Schwartz J, Shindell D, et al. Global Air Quality and Health Co-benefits of Mitigating Near-Term Climate Change through Methane and Black Carbon Emission Controls [J]. Environmental Health Perspectives, 2012,120(6):831-839.

    [15] Lu X, Hong J, Zhang L, et al. Severe surface ozone pollution in China: a global perspective [J]. Environ. Sci. Technol. Lett., 2018,5(8):487- 494.

    [16] Shen L, Jacob D J, Liu X, et al. An evaluation of the ability of the Ozone Monitoring Instrument (OMI) to observe boundary layer ozone pollution across China: application to 2005~2017 ozone trends [J]. Atmospheric Chemistry and Physics, 2019,19(9):6551-6560.

    [17] Shao M, Lu S, Liu Y, et al. Volatile organic compounds measured in summer in Beijing and their role in ground-level ozone formation [J]. Journal of Geophysical Research, 2009,114(7),D00G06.

    [18] Lu K, Zhang Y, Su H, et al. Oxidant (O3+NO2) production processes and formation regimes in Beijing [J]. Journal of Geophysical Research, 2010,115:1-18.

    [19] Ran L, Zhao C, Xu W, et al. VOC reactivity and its effect on ozone production during the HaChi summer campaign [J]. Atmospheric Chemistry and Physics, 2011,11(10):4657-4667.

    [20] Li L, Chen C, Huang C, et al. Process analysis of regional ozone formation over the Yangtze River Delta, China using the Community Multi-scale Air Quality modeling system [J]. Atmospheric Chemistry and Physics, 2012,12(22):10971-10987.

    [21] Liu Z, Wang Y, Gu D, et al. Summertime photochemistry during CAREBeijing-2007: ROxbudgets and O3formation [J]. Atmospheric Chemistry and Physics, 2012,12(16):7737-7752.

    [22] Ma J, Xu X, Zhao C, et al. A review of atmospheric chemistry research in China: Photochemical smog, haze pollution, and gas-aerosol interactions [J]. Advances in atmospheric sciences, 2012,29(5):1006- 1026.

    [23] Xue L, 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.

    [24] Wang T, Xue L, Brimblecombe P, et al. Ozone pollution in China: A review of concentrations, meteorological influences, chemical precursors, and effects [J]. Sci Total Environ, 2017,575:1582-1596.

    [25] Tang G, Wang Y, Li X, et al. Spatial-temporal variations in surface ozone in Northern China as observed during 2009~2010 and possible implications for future air quality control strategies [J]. Atmospheric Chemistry and Physics, 2012,12(5):2757-2776.

    [26] Ding A, Fu C, Yang X, Sun, et al. Ozone and fine particle in the western Yangtze River Delta: an overview of 1yr data at the SORPES station [J]. Atmospheric Chemistry and Physics, 2013,13(11):5813- 5830.

    [27] Wang X, Zhang Y, Hu Y, et al. Process analysis and sensitivity study of regional ozone formation over the Pearl River Delta, China, during the PRIDE-PRD2004 campaign using the Community Multiscale Air Quality modeling system [J].Atmospheric Chemistry and Physics, 2010,10(9):4423-4437.

    [28] Guo H, Wanga T, Simpsonb I, et al. Source contributions to ambient VOCs and CO at a rural site in eastern China, Atmospheric Environ, 2004,38(27):4551- 4560.

    [29] 蔡長杰,耿福海,俞 瓊,等.上海中心城區(qū)夏季揮發(fā)性有機物(VOCs)的源解析[J]. 環(huán)境科學學報, 2010,30(5):926-934.

    Cai C, Geng F, Yu Q, et al. Source apportionment of VOCs at city centre of Shanghai in summer [J]. Acta Scientiae Circumstantiae, 2010,30(5):926-934.

    [30] Zhang L, Li H, Wu Z, et al. Characteristics of atmospheric volatile organic compounds in urban area of Beijing: Variations, photochemical reactivity and source apportionment [J]. Journal of Environmental Sciences, 2020,95:190-200.

    [31] Qin Y, Li J, Gong K, et al. Double high pollution events in the Yangtze River Delta from 2015 to 2019: Characteristics, trends, and meteorological situations [J]. Science of the Total Environment, 2021,792(12):148349.

    [32] 洪瑩瑩,翁佳烽,譚浩波,等.珠江三角洲秋季典型O3污染的氣象條件及貢獻量化[J]. 中國環(huán)境科學, 2021,41(1):1-10.

    Hong Y, Weng J, Tan H, et al. Meteorological conditions and contribution quantification of typical ozone pollutionduring autumnin Pearl River Delta [J] China Environmental Science, 2021,41(1):1-10.

    [33] 宋夢迪,馮 淼,李 歆,等.成都市臭氧重污染成因與來源解析[J/OL]. 中國環(huán)境科學:1-11[2021-10-26].https://doi.org/10.19674/ j.cnki.issn1000-6923.20210923.004.

    Song M, Feng M, Li X, et al. Causes and sources of heavy ozone pollution in Chengdu [J/OL]. China Environmental Science: - 11 [2021-10-261] ttps://oiorg/10.19674/j.cnki.issn1000-6923.20210923.004.

    [34] Liang B, Yu X, Mi H, et al. Health risk assessment and source apportionment of VOCs inside new vehicle cabins: A case study from Chongqing, China [J]. Atmospheric Pollution Research, 2019,10(5): 1677-1684.

    [35] 白永清,祁海霞,劉 琳,等.武漢大氣能見度與PM2.5濃度及相對濕度關(guān)系的非線性分析及能見度預報[J]. 氣象學報, 2016,74(2): 189-199.

    Bai Y, Qi H, Liu L, et al. Study on the nonlinear relationship among the visibility, PM2.5concentration and relative humidity in Wuhan and the visibility prediction [J]. Acta Meteorologica Sinica, 2016,74(2): 189-199.

    [36] 譚成好,趙天良,崔春光,等.近50年華中地區(qū)霾污染的特征[J]. 中國環(huán)境科學, 2015,35(8):2272-2280.

    Tan C, Zhao T, Cui C, et al. Characteristics of haze pollution over Central China duringthe past 50years [J]. China Environmental Science, 2015,35(8):2272–2280.

    [37] 馬德栗,李 蘭,鞠英芹.湖北省霾日數(shù)氣候特征及夏季典型霾過程氣象因子分析[J]. 環(huán)境科學與技術(shù), 2015,v.38(11):148-153.

    Ma D, Li L, Ju Y, et al. Climatic characteristics of haze days and analysis of summer haze weather event in Hubei province [J]. Environmental Science &Technology, 2015,v.38(1):148-153.

    [38] Yang Y, Liu X, Zheng J, et al. Characteristics of one-year observation of VOCs, NO, and O3at an urban site in Wuhan, China [J]. Journal of Environmental Sciences, 2019,79(5):297-310.

    [39] 孫天樂,何凌燕,曾立武,等.2008北京殘奧會期間大氣黑碳氣溶膠污染特征[J]. 中國環(huán)境科學, 2012,32(12):2123-2127.

    Sun T, He L, Zeng L, et al. Black carbon measurement during Beijing Paralympic Game [J]. China Environmental Science, 2012,32(12): 2123-2127.

    [40] 王羽琴,張元勛,張 陽,等.深圳大運會前后大氣含碳氣溶膠污染特征 [J]. 中國環(huán)境科學, 2014,34(8):1972-1978.

    Wang Y, Zhang Y, Zhang Y, et al. Characterzation of carbonaceous aerosols during and post-Shenzhen UNIVERSIADE period [J]. China Environmental Science, 2014,34(8):1972-1978.

    [41] 周國治,宋冰冰,羅岳平,等.2016年長沙市環(huán)境空氣質(zhì)量的“周末效應”及“假日效應” [J]. 中國環(huán)境監(jiān)測, 2018,34(4):68-76.

    Zhou G, Song B, Luo Y, et al. Holiday effects on ambient air quality in Changsha during 2016 [J]. Environmental Monitoring of China, 2018, 34(4):68-76.

    [42] Yan Y, Zheng H, Kong S, et al.On the local anthropogenic source diversities and transboundary transport for urban agglomeration ozone mitigation [J]. Atmospheric Environment, 2021,245(6):118005.

    [43] Zheng H, Kong S, Chen N, et al. Source apportionment of volatile organic compounds: Implications to reactivity, ozone formation, and secondary organic aerosol potential [J]. Atmospheric Research, 2021, 249:105344.

    [44] Liu Y, Song M, Liu X, et al. Characterization and sources of volatile organic compounds (VOCs) and their related changes during ozone pollution days in 2016 in Beijing, China [J]. Environmental Pollution, 2020,257:113599.

    [45] 鄭 玫,張延君,閆才青,等.中國PM2.5來源解析方法綜述[J]. 北京大學學報(自然科學版), 2014,50(6):1141-1154.

    Zheng M, Zhang Y, Yan C, et al. Review of PM2.5source apportionment methods in China [J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2014,50(6):1141-1154.

    [46] 張延君,鄭 玫,蔡 靖,等.PM2.5源解析方法的比較與評述[J]. 科學通報, 2015,60(2):109-121,1-2.

    Zhang Y, Zheng M, Cai J, et al. Comparison and overview of PM2.5source apportionment methods [J]. Chinese Science Bulletin, 2015,60(2):109-121,1-2.

    [47] 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.

    [48] Shao P, An J, Xin J, et al. Source apportionment of VOCs and the contribution to photochemical ozone formation during summer in the typical industrial area in the Yangtze River Delta, China [J]. Atmospheric Research, 2016,176:64-74.

    [49] 牛真真,孔少飛,嚴 沁,等.薪柴和經(jīng)濟作物秸稈燃燒VOCs排放特征[J]. 環(huán)境科學, 2020,41(3):1107-1115.

    Niu Z, Kong S, Yan Q, et al. Profile characteristics of VOCs from wood and economic crop burning [J]. Environmental Science, 2020, 41(3):1107-1115.

    [50] Zhang L, Li H, Wu Z, et al. Characteristics of atmospheric volatile organic compounds in urban area of Beijing: Variations, photochemical reactivity and source apportionment [J]. Journal of Environmental Sciences, 2020,95(9):190-200.

    [51] McCarthy M C, Aklilu Y, Brown S G, et al. Source apportionment of volatile organic compounds measured in Edmonton, Alberta [J]. Atmospheric Environment, 2013,81:504-516.

    [52] 齊一謹,史來文,陳紹鵬,等.鄭州市金水區(qū)夏季VOCs污染特征、來源及風險評估[J]. 河南科學, 2021,39(9):1504-1512.

    Qi Y, Shi L, Chen S, et al. Pollution characteristic, source apportionment and risk assessment of volatile organic compounds from Jinshui district of Zhengzhou city in summer [J]. Henan Science, 2021,39(9):1504-1512.

    [53] 王成輝,陳軍輝,韓 麗,等.成都市城區(qū)大氣VOCs季節(jié)污染特征及來源解析 [J]. 環(huán)境科學, 2020,41(9):3951-3960.

    Wang C, Chen J, Han L, et al. Seasonal pollution characteristics and analysis of the sources of atmospheric VOCs in Chengdu urban area [J]. Environmental Science, 2020,41(9):3951-3960.

    [54] Ou J, Guo H, Zheng J, et al. Concentrations and sources of non- methane hydrocarbons (NMHCs) from 2005 to 2013 in Hongkong: A multi-year real-time data analysis [J]. Atmospheric Environment, 2015,103:196-206.

    [55] Shen L, Xiang P, Liang S, et al. Sources profiles of volatile organic compounds (VOCs) measured in a typical industrial process in Wuhan, Central China [J]. Atmosphere, 2018,9(8):297.

    [56] Liu Y, Shao M, Fu L, et al. Source profiles of volatile organic compounds (VOCs) measured in China: Part I [J]. Atmospheric Environment, 2008,42(25):6247-6260.

    [57] Sun J, Wu F, Hu B, et al. VOC characteristics, emissions and contributions to SOA formation during hazy episodes [J]. Atmospheric Environment, 2016,141:560-570.

    [58] Liu C, Ma Z, Mu Y, et al. The levels, variation characteristics, and sources of atmospheric non-methane hydrocarbon compounds during wintertime in Beijing, China [J]. Atmospheric Chemistry and Physics, 2017,17(17):10633-10649.

    [59] 黃嫣旻,魏海萍,段玉森,等.上海世博會環(huán)境空氣質(zhì)量狀況和原因分析 [J]. 中國環(huán)境監(jiān)測, 2013,29(5):58-63.

    Huang Y, Wei H, Duan Y, et al. Ambient air quality status and reason analysis of Shanghai World Expo [J] Environmental Monitoring of China, 2013,29(5):58-63.

    [60] Huang K, Zhuang G, Lin Y, et al. How to improve the air quality over megacities in China: pollution characterization and source analysis in Shanghai before, during, and after the 2010 World Expo [J]. Atmospheric Chemistry and Physics, 2013,13(12):5927-5942.

    [61] Li Y, Xu H. Assessment of reductions in emission-driven air pollution during the Beijing Olympic Games, Shanghai World Expo, Guangzhou Asian Games and Wuhan COVID-19 Lockdown [J]. Aerosol and Air Quality Research, 2021,21(9),https://doi.org/10.4209/ aaqr.200644.

    Characteristics and sources of ozone and its precursors aroundthe Wuhan Military Games.

    MA Jing1, YAN Ying-ying1*, KONG Shao-fei1,2, SONG Ai-li1, CHEN Nan2,3, ZHU Bo2,3, QUAN Ji-hong2,3, QI Shi-hua1,2

    (1.School of Environmental Studies, China University ofGeosciences, Wuhan 430074, China;2.Atmospheric Combined Pollution Research Center of Hubei, Wuhan 430074, China;3.Environmental Monitoring Central Station in Hubei Province, Wuhan 430072, China)., 2022,42(7):3023~3032

    The hourly data of ozone, NO(= NO + NO2) and VOCs (102species) during October~November 2019 provided by Hubei Superstationwere used to analyze the variation characteristics of ozone pollutionduring the Wuhan Military Games (WMG). The photochemical regime of ozone was simulated based on DSMACC box model. The source apportionment of VOCs was conducted by PMF model, and the ozone formation potential of different VOCs sourcewas estimated. The results showed that the maximum daily 8-hour average ozone concentration (maximum MDA8:219.51mg/m3) exceeded the national level II standard before the military games. The MDA8 values (135.11mg/m3on average) decreased significantly during the military games, and the concentrations (140.98mg/m3) rise after the WMG. The difference of ozone before and during the military games was mainly affected by meteorological conditions, and the increase of ozone after the WMG was mainly due to the significant increase of precursor emissions. The EKMA curve indicates that the photochemical formation of ozone was controlled by a VOC-limited regime during the study period. The VOCssources that contribute greatly to ozone formation before the military games were combustion, petrochemical industry and motor vehicles, accounting for 23.0%, 22.8% and 22.5%, respectively. The main sources of VOCs during the military games were motor vehicles (38.4%) and combustion (25.5%). While petrochemical industry (32.6%) and fuel volatilization (25.7%)were the main sources after the WMG. The results showed that the prevention strategies of the military games had a significant effect on the emission reduction of petrochemical industry, but the emission mitigation of motor vehicles and combustion was not significant. Wuhan should pay more attention to the control of combustion, fuel volatilization and motor vehicle emissions.

    ozone pollution;Wuhan Military Games;VOCs source apportionment;EKMA curve

    X511

    A

    1000-6923(2022)07-3023-10

    馬 靜(1998-),男,四川南充人,中國地質(zhì)大學(武漢),碩士研究生,主要從事對流層臭氧污染研究.

    2021-12-01

    國家重點研發(fā)計劃(2017YFC0212602);國家自然科學基金資助項目(41905112)

    * 責任作者, 副教授, yanyingying@cug.edu.cn

    猜你喜歡
    貢獻率機動車臭氧
    讓機動車交通安全統(tǒng)籌更
    公民與法治(2022年7期)2022-07-22 07:12:22
    由一起廠內(nèi)機動車事故引發(fā)的思考
    文印室內(nèi)臭氧散發(fā)實測分析及模擬
    一種通用的裝備體系貢獻率評估框架
    鐵路機動車管理信息系統(tǒng)
    電子制作(2019年24期)2019-02-23 13:22:30
    關(guān)于裝備體系貢獻率研究的幾點思考
    看不見的污染源——臭氧
    利用臭氧水防治韭菜遲眼蕈蚊
    中國蔬菜(2016年8期)2017-01-15 14:23:43
    В первой половине 2016 года вклад потребления в рост китайской экономики достиг 73,4 процента
    中亞信息(2016年10期)2016-02-13 02:32:45
    臭氧分子如是說
    在线观看免费午夜福利视频| 最近最新免费中文字幕在线| 欧美在线一区亚洲| 极品人妻少妇av视频| 色尼玛亚洲综合影院| 人人澡人人妻人| 国产精品1区2区在线观看.| 午夜免费成人在线视频| 国产一区二区激情短视频| 1024视频免费在线观看| 成人18禁高潮啪啪吃奶动态图| 少妇熟女aⅴ在线视频| 亚洲色图av天堂| 脱女人内裤的视频| 久久久久久久久免费视频了| 少妇的丰满在线观看| 国产熟女午夜一区二区三区| 69av精品久久久久久| 成人免费观看视频高清| 黄色视频,在线免费观看| 嫩草影院精品99| 免费搜索国产男女视频| 青草久久国产| 日韩一卡2卡3卡4卡2021年| 亚洲男人的天堂狠狠| 免费观看精品视频网站| 午夜福利影视在线免费观看| 国产亚洲av高清不卡| 神马国产精品三级电影在线观看 | 人成视频在线观看免费观看| tocl精华| 精品熟女少妇八av免费久了| 男男h啪啪无遮挡| 黄色 视频免费看| 岛国视频午夜一区免费看| 久久久久九九精品影院| 9色porny在线观看| 天天躁夜夜躁狠狠躁躁| 国产精品精品国产色婷婷| 国产麻豆69| 亚洲精品av麻豆狂野| 久久婷婷成人综合色麻豆| 男女床上黄色一级片免费看| 国产欧美日韩综合在线一区二区| 精品卡一卡二卡四卡免费| 大型黄色视频在线免费观看| 成在线人永久免费视频| 久久久久国产一级毛片高清牌| 成人永久免费在线观看视频| av欧美777| cao死你这个sao货| 多毛熟女@视频| 老司机在亚洲福利影院| 久久中文看片网| 日韩大码丰满熟妇| 亚洲精品av麻豆狂野| 久久久久久人人人人人| 天天添夜夜摸| 久久久久久亚洲精品国产蜜桃av| 老熟妇乱子伦视频在线观看| 国产精品免费一区二区三区在线| 欧美人与性动交α欧美精品济南到| 国产在线精品亚洲第一网站| 国产一区二区三区视频了| 欧美黑人欧美精品刺激| 非洲黑人性xxxx精品又粗又长| 波多野结衣巨乳人妻| www.自偷自拍.com| 大型av网站在线播放| 午夜久久久在线观看| 亚洲一区中文字幕在线| 国产高清激情床上av| 夜夜爽天天搞| www.熟女人妻精品国产| 激情视频va一区二区三区| 天堂动漫精品| 别揉我奶头~嗯~啊~动态视频| 国产精品电影一区二区三区| 午夜精品国产一区二区电影| 俄罗斯特黄特色一大片| 女人被狂操c到高潮| av电影中文网址| 村上凉子中文字幕在线| 一级片免费观看大全| 午夜成年电影在线免费观看| 亚洲成人免费电影在线观看| 黄色视频不卡| 国产精品亚洲av一区麻豆| 天天一区二区日本电影三级 | 久久热在线av| 伊人久久大香线蕉亚洲五| 亚洲av日韩精品久久久久久密| 很黄的视频免费| 一本大道久久a久久精品| 欧美av亚洲av综合av国产av| 高清毛片免费观看视频网站| 国产1区2区3区精品| 亚洲成人久久性| 欧美日韩中文字幕国产精品一区二区三区 | 首页视频小说图片口味搜索| 国产精品99久久99久久久不卡| avwww免费| 色精品久久人妻99蜜桃| 啦啦啦观看免费观看视频高清 | 久久精品国产99精品国产亚洲性色 | 日本免费一区二区三区高清不卡 | 免费在线观看亚洲国产| 国产亚洲av嫩草精品影院| 在线观看日韩欧美| 真人做人爱边吃奶动态| 巨乳人妻的诱惑在线观看| 久久亚洲真实| 亚洲在线自拍视频| 涩涩av久久男人的天堂| 啪啪无遮挡十八禁网站| 脱女人内裤的视频| 亚洲少妇的诱惑av| 18美女黄网站色大片免费观看| 国产日韩一区二区三区精品不卡| 9热在线视频观看99| 国产精品久久久久久精品电影 | 黄色a级毛片大全视频| 黄色 视频免费看| 午夜福利免费观看在线| 免费久久久久久久精品成人欧美视频| 亚洲五月色婷婷综合| 国产91精品成人一区二区三区| 精品欧美国产一区二区三| 国产成人免费无遮挡视频| av天堂在线播放| av视频在线观看入口| 午夜成年电影在线免费观看| 亚洲 国产 在线| av在线天堂中文字幕| 欧美成人免费av一区二区三区| 国语自产精品视频在线第100页| 91麻豆精品激情在线观看国产| 如日韩欧美国产精品一区二区三区| 亚洲精品av麻豆狂野| 国产野战对白在线观看| 一边摸一边抽搐一进一小说| 怎么达到女性高潮| 不卡av一区二区三区| av天堂久久9| 国内精品久久久久久久电影| 亚洲精品久久国产高清桃花| 大型av网站在线播放| x7x7x7水蜜桃| 我的亚洲天堂| 国产精品免费视频内射| 亚洲欧美日韩无卡精品| 亚洲人成77777在线视频| 男人操女人黄网站| 久久久久精品国产欧美久久久| 成人永久免费在线观看视频| 午夜视频精品福利| 国产私拍福利视频在线观看| 母亲3免费完整高清在线观看| 日本vs欧美在线观看视频| 久久久久久久久中文| 嫩草影视91久久| 国产精品一区二区三区四区久久 | 色综合欧美亚洲国产小说| 夜夜夜夜夜久久久久| 正在播放国产对白刺激| 久久精品亚洲熟妇少妇任你| 久久久久亚洲av毛片大全| 免费高清视频大片| 免费看十八禁软件| 精品欧美一区二区三区在线| 国产精品电影一区二区三区| 一区二区三区激情视频| 99国产极品粉嫩在线观看| 午夜福利一区二区在线看| 午夜福利免费观看在线| 电影成人av| 午夜老司机福利片| 久久久久精品国产欧美久久久| 精品卡一卡二卡四卡免费| av免费在线观看网站| 男人舔女人的私密视频| 色在线成人网| 色综合亚洲欧美另类图片| 亚洲性夜色夜夜综合| 国产av又大| 亚洲精品国产一区二区精华液| 黄色毛片三级朝国网站| 麻豆成人av在线观看| 男女午夜视频在线观看| 亚洲成人国产一区在线观看| 男女下面插进去视频免费观看| 国产熟女午夜一区二区三区| 精品久久久久久成人av| 国产精品野战在线观看| 国产精品99久久99久久久不卡| 最好的美女福利视频网| 国产精品,欧美在线| 久久久久久国产a免费观看| 18禁裸乳无遮挡免费网站照片 | 国产精品爽爽va在线观看网站 | 天堂√8在线中文| 狠狠狠狠99中文字幕| 18禁黄网站禁片午夜丰满| 亚洲精品国产色婷婷电影| 国产亚洲精品一区二区www| 亚洲中文字幕日韩| 美女大奶头视频| 久久久久久大精品| 亚洲精品久久国产高清桃花| 黄片播放在线免费| 美女国产高潮福利片在线看| 丝袜在线中文字幕| 两个人免费观看高清视频| 久久久水蜜桃国产精品网| 色精品久久人妻99蜜桃| 亚洲一区高清亚洲精品| 夜夜躁狠狠躁天天躁| 久热爱精品视频在线9| 亚洲精品一区av在线观看| 黄片小视频在线播放| 日本a在线网址| 人人妻人人澡人人看| 亚洲欧美日韩高清在线视频| 久热爱精品视频在线9| 国产又爽黄色视频| 国产精品自产拍在线观看55亚洲| 18禁黄网站禁片午夜丰满| 多毛熟女@视频| 色婷婷久久久亚洲欧美| 99国产极品粉嫩在线观看| 久久精品国产清高在天天线| 欧美日韩亚洲综合一区二区三区_| 香蕉国产在线看| 婷婷丁香在线五月| 精品卡一卡二卡四卡免费| 日本黄色视频三级网站网址| 成人18禁在线播放| 香蕉久久夜色| 午夜精品国产一区二区电影| 国产精品,欧美在线| a级毛片在线看网站| 性欧美人与动物交配| 可以在线观看毛片的网站| 一边摸一边做爽爽视频免费| 日本欧美视频一区| 亚洲av电影在线进入| 国产亚洲av嫩草精品影院| 国产精品永久免费网站| 日韩国内少妇激情av| 自线自在国产av| 人人妻人人爽人人添夜夜欢视频| 免费看十八禁软件| 午夜福利成人在线免费观看| 韩国av一区二区三区四区| 午夜福利高清视频| 亚洲av成人av| 最近最新免费中文字幕在线| 精品久久久久久成人av| 成人国产综合亚洲| 亚洲美女黄片视频| av在线天堂中文字幕| bbb黄色大片| 久久人人爽av亚洲精品天堂| 亚洲av五月六月丁香网| 91成人精品电影| 黄色视频,在线免费观看| e午夜精品久久久久久久| 99国产精品99久久久久| 成人国产一区最新在线观看| 在线天堂中文资源库| 精品国产一区二区久久| 丝袜美足系列| 亚洲成人免费电影在线观看| 男人舔女人下体高潮全视频| 亚洲性夜色夜夜综合| 欧美色欧美亚洲另类二区 | 国产精品二区激情视频| 母亲3免费完整高清在线观看| 精品日产1卡2卡| 激情在线观看视频在线高清| 国产不卡一卡二| videosex国产| 99香蕉大伊视频| 日日爽夜夜爽网站| 女人被躁到高潮嗷嗷叫费观| 99热只有精品国产| 久久九九热精品免费| 9色porny在线观看| 亚洲欧美激情在线| 一进一出好大好爽视频| 欧美成人一区二区免费高清观看 | 久久精品亚洲熟妇少妇任你| 日本欧美视频一区| 一本久久中文字幕| 成人18禁在线播放| 黑人巨大精品欧美一区二区mp4| 97人妻精品一区二区三区麻豆 | 男女下面进入的视频免费午夜 | 黄色女人牲交| 色精品久久人妻99蜜桃| 欧美久久黑人一区二区| 亚洲五月色婷婷综合| 精品乱码久久久久久99久播| 麻豆av在线久日| 久久热在线av| 久久婷婷成人综合色麻豆| 长腿黑丝高跟| 日本vs欧美在线观看视频| 亚洲va日本ⅴa欧美va伊人久久| 日韩精品免费视频一区二区三区| 免费人成视频x8x8入口观看| 国产精品久久电影中文字幕| 咕卡用的链子| 亚洲少妇的诱惑av| 日本a在线网址| 日本三级黄在线观看| ponron亚洲| 久久精品国产亚洲av高清一级| 桃色一区二区三区在线观看| 女性被躁到高潮视频| 在线永久观看黄色视频| 如日韩欧美国产精品一区二区三区| 免费在线观看日本一区| 欧美一区二区精品小视频在线| 琪琪午夜伦伦电影理论片6080| 一进一出抽搐动态| 亚洲午夜理论影院| 人人妻人人澡欧美一区二区 | 国产亚洲精品一区二区www| 美女免费视频网站| 国产成人影院久久av| 两个人看的免费小视频| 999久久久精品免费观看国产| 成人免费观看视频高清| 成年女人毛片免费观看观看9| 咕卡用的链子| 久久国产精品男人的天堂亚洲| cao死你这个sao货| 黄色丝袜av网址大全| 亚洲中文字幕一区二区三区有码在线看 | 一区二区三区激情视频| 国产精品久久久久久精品电影 | 可以免费在线观看a视频的电影网站| 又黄又爽又免费观看的视频| 亚洲激情在线av| 午夜免费激情av| 国产精品野战在线观看| 最好的美女福利视频网| 亚洲黑人精品在线| 久久国产亚洲av麻豆专区| av中文乱码字幕在线| 日韩av在线大香蕉| 夜夜躁狠狠躁天天躁| 999久久久国产精品视频| av在线天堂中文字幕| 亚洲av电影不卡..在线观看| 女同久久另类99精品国产91| 国产欧美日韩精品亚洲av| 操美女的视频在线观看| av中文乱码字幕在线| 两个人视频免费观看高清| 久久热在线av| 身体一侧抽搐| 国产一区二区三区在线臀色熟女| 一卡2卡三卡四卡精品乱码亚洲| 精品少妇一区二区三区视频日本电影| 变态另类丝袜制服| 热99re8久久精品国产| 两性午夜刺激爽爽歪歪视频在线观看 | 日本免费a在线| 亚洲avbb在线观看| 一级a爱视频在线免费观看| 国产欧美日韩精品亚洲av| 国产成人一区二区三区免费视频网站| 亚洲av成人一区二区三| 啦啦啦免费观看视频1| 中亚洲国语对白在线视频| 怎么达到女性高潮| 大型黄色视频在线免费观看| 午夜亚洲福利在线播放| 欧美成人免费av一区二区三区| 老司机在亚洲福利影院| 亚洲男人天堂网一区| 欧美日韩福利视频一区二区| 热re99久久国产66热| 97人妻精品一区二区三区麻豆 | 欧美黑人精品巨大| 香蕉国产在线看| 女人爽到高潮嗷嗷叫在线视频| 精品人妻在线不人妻| 香蕉国产在线看| 91九色精品人成在线观看| 久久精品aⅴ一区二区三区四区| 99在线视频只有这里精品首页| 色播在线永久视频| 巨乳人妻的诱惑在线观看| 夜夜爽天天搞| 日韩有码中文字幕| 免费在线观看视频国产中文字幕亚洲| 久热这里只有精品99| 69精品国产乱码久久久| 两性午夜刺激爽爽歪歪视频在线观看 | av视频在线观看入口| 国产片内射在线| 免费在线观看影片大全网站| 日本免费a在线| 黄片播放在线免费| 亚洲伊人色综图| 免费人成视频x8x8入口观看| 啦啦啦观看免费观看视频高清 | 国产精品秋霞免费鲁丝片| 日韩精品中文字幕看吧| 亚洲男人天堂网一区| aaaaa片日本免费| 黄色成人免费大全| 久久久水蜜桃国产精品网| 欧美成人免费av一区二区三区| 亚洲 欧美 日韩 在线 免费| av有码第一页| 老司机午夜福利在线观看视频| 夜夜看夜夜爽夜夜摸| 人成视频在线观看免费观看| 免费在线观看视频国产中文字幕亚洲| av有码第一页| 国产精品一区二区在线不卡| 亚洲国产毛片av蜜桃av| 日本撒尿小便嘘嘘汇集6| 999久久久国产精品视频| 久久久久九九精品影院| 欧美日韩精品网址| 亚洲一码二码三码区别大吗| 在线国产一区二区在线| 宅男免费午夜| 大码成人一级视频| 妹子高潮喷水视频| 欧美在线一区亚洲| 男男h啪啪无遮挡| 午夜免费观看网址| www.熟女人妻精品国产| 中文字幕人成人乱码亚洲影| 午夜精品久久久久久毛片777| 免费在线观看影片大全网站| 午夜福利欧美成人| 国产区一区二久久| 黑人巨大精品欧美一区二区蜜桃| 午夜精品在线福利| 国产亚洲av高清不卡| 久久精品国产亚洲av高清一级| 人妻丰满熟妇av一区二区三区| 日韩欧美一区视频在线观看| 国产亚洲欧美98| 韩国精品一区二区三区| 99香蕉大伊视频| 久久久久久免费高清国产稀缺| 日本 欧美在线| 欧美精品啪啪一区二区三区| 久久人妻熟女aⅴ| 精品久久久精品久久久| 亚洲熟妇熟女久久| 久久人人爽av亚洲精品天堂| 国产精品久久久人人做人人爽| 午夜福利成人在线免费观看| 午夜福利18| 精品国产一区二区久久| 一区二区三区高清视频在线| a在线观看视频网站| 欧美精品亚洲一区二区| 免费搜索国产男女视频| 又黄又粗又硬又大视频| 少妇裸体淫交视频免费看高清 | 亚洲视频免费观看视频| 久久伊人香网站| 亚洲av成人不卡在线观看播放网| 麻豆国产av国片精品| 亚洲精品在线美女| 1024视频免费在线观看| 九色亚洲精品在线播放| cao死你这个sao货| 超碰成人久久| 免费看美女性在线毛片视频| 国产亚洲精品一区二区www| 亚洲av电影在线进入| 国产精品野战在线观看| 日韩三级视频一区二区三区| 夜夜躁狠狠躁天天躁| 亚洲成av片中文字幕在线观看| 国产片内射在线| 国产高清激情床上av| 久久草成人影院| 天天躁狠狠躁夜夜躁狠狠躁| 在线av久久热| 女人高潮潮喷娇喘18禁视频| 精品无人区乱码1区二区| 亚洲专区中文字幕在线| 99精品久久久久人妻精品| 国产一区二区三区视频了| 欧美成狂野欧美在线观看| 亚洲人成网站在线播放欧美日韩| 夜夜爽天天搞| 在线观看免费日韩欧美大片| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品中文字幕一二三四区| 亚洲午夜理论影院| 国产av精品麻豆| 女人爽到高潮嗷嗷叫在线视频| 国产精品自产拍在线观看55亚洲| 亚洲色图综合在线观看| 欧美日本中文国产一区发布| 午夜免费鲁丝| 亚洲av美国av| 精品第一国产精品| 亚洲av美国av| 国产aⅴ精品一区二区三区波| 一区二区日韩欧美中文字幕| 国产亚洲欧美精品永久| 国产又色又爽无遮挡免费看| 欧美一级毛片孕妇| 精品久久久久久,| 亚洲一区高清亚洲精品| 国产三级在线视频| 久久这里只有精品19| 神马国产精品三级电影在线观看 | 亚洲精品久久成人aⅴ小说| 亚洲精品国产一区二区精华液| 国产不卡一卡二| 成熟少妇高潮喷水视频| 午夜a级毛片| 国产精品久久电影中文字幕| 国产主播在线观看一区二区| 亚洲第一电影网av| 午夜福利,免费看| 欧美久久黑人一区二区| 欧美黄色片欧美黄色片| 欧美久久黑人一区二区| 日韩欧美三级三区| 亚洲情色 制服丝袜| 国产成年人精品一区二区| 9色porny在线观看| 欧美黑人欧美精品刺激| 久久青草综合色| 黄色视频不卡| 男男h啪啪无遮挡| 亚洲人成77777在线视频| 亚洲专区字幕在线| 国产亚洲精品久久久久5区| 天天一区二区日本电影三级 | 亚洲国产精品合色在线| av视频在线观看入口| 女性被躁到高潮视频| 亚洲五月天丁香| 黄网站色视频无遮挡免费观看| 午夜福利在线观看吧| 9热在线视频观看99| 麻豆av在线久日| 成人特级黄色片久久久久久久| 午夜福利视频1000在线观看 | 最近最新免费中文字幕在线| 欧美成人午夜精品| 黄片大片在线免费观看| 亚洲成人久久性| 国产麻豆69| 亚洲成人精品中文字幕电影| 中出人妻视频一区二区| 精品人妻在线不人妻| 久久热在线av| 国产午夜精品久久久久久| 中文字幕av电影在线播放| 99国产综合亚洲精品| 天天一区二区日本电影三级 | 在线观看免费视频日本深夜| tocl精华| 欧美+亚洲+日韩+国产| 免费不卡黄色视频| 亚洲欧美日韩无卡精品| 成人亚洲精品av一区二区| 免费av毛片视频| 亚洲欧美日韩高清在线视频| 亚洲无线在线观看| 亚洲自偷自拍图片 自拍| 男女下面插进去视频免费观看| 国产精品1区2区在线观看.| 一本久久中文字幕| 纯流量卡能插随身wifi吗| 最新在线观看一区二区三区| 男人舔女人下体高潮全视频| 国产日韩一区二区三区精品不卡| 长腿黑丝高跟| 精品久久久精品久久久| 亚洲黑人精品在线| 在线观看66精品国产| 欧美色视频一区免费| 久久人人97超碰香蕉20202| 两人在一起打扑克的视频| 久久狼人影院| 日韩欧美国产在线观看| 亚洲成av片中文字幕在线观看| 色综合婷婷激情| 一夜夜www| 亚洲欧美精品综合一区二区三区| 丝袜人妻中文字幕| 免费在线观看影片大全网站| aaaaa片日本免费| 日韩精品免费视频一区二区三区| 天堂影院成人在线观看| 久久精品人人爽人人爽视色| 亚洲av五月六月丁香网| 亚洲国产精品成人综合色| 在线观看www视频免费| 免费在线观看黄色视频的| 久久人人精品亚洲av| 国产精品 欧美亚洲| 一二三四社区在线视频社区8| 成人欧美大片| 午夜亚洲福利在线播放| 国产乱人伦免费视频|