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

    西安第十四屆全國運動會及殘?zhí)貖W會前后PM2.5溯源研究

    2024-05-23 18:00:04李香凝王羽琴高燕王競錚臧匯宇謝林花
    陜西科技大學(xué)學(xué)報 2024年3期

    李香凝 王羽琴 高燕 王競錚 臧匯宇 謝林花

    文章編號:2096-398X2024)03-0027-11

    (1.陜西科技大學(xué) 環(huán)境科學(xué)與工程學(xué)院, 陜西 西安 710021; 2.隆基綠能光伏科技(西咸新區(qū))有限公司, 陜西 西安 710021)

    摘 要:基于2021年9月~11月西安市環(huán)境空氣質(zhì)量監(jiān)測數(shù)據(jù),分析十四運會及殘?zhí)貖W會前中后期西安市區(qū)域PM2.5的污染特征,利用HYSPLIT模型、潛在源貢獻(xiàn)因子法(PSCF)和權(quán)重濃度軌跡法(CWT)確定氣團(tuán)輸送路徑、潛在源區(qū)及其對PM2.5質(zhì)量濃度的貢獻(xiàn).結(jié)果表明:研究期間西安市以東南風(fēng)為主導(dǎo),PM2.5濃度在十四運會期間最低,殘?zhí)貖W會結(jié)束后最高.PM2.5的氣團(tuán)輸送路徑以陜南附近的本地輸送、偏東和偏南方向的近距離輸送以及西北方向的遠(yuǎn)距離輸送為主.PM2.5的潛在源區(qū)主要分布在河南湖北、貴州四川重慶、新疆內(nèi)蒙地區(qū)以及陜南各城市及周邊省份交界處.

    關(guān)鍵詞:十四運; 細(xì)顆粒物; 后向軌跡; 潛在源分析; 濃度權(quán)重分析

    中圖分類號:X513??? 文獻(xiàn)標(biāo)志碼: A

    Traceability study of PM2.5 before,during and after the 14th National Games of the People′s Republic of China & the National Games for Persons with Disabilities in Xi′an City

    LI Xiang-ning1, WANG Yu-qin1*, GAO Yan1,2, WANG Jing-zheng1, ANG Hui-yu1, XIE Lin-hua1

    1.School of Environmental Science and Engineering, Shaanxi University of Science & Technology, Xi′an 710021, China; 2.Longi Green Energy Photovoltaic Technology Xixian New Area) Co., Ltd., Xi′an 710021, China)

    Abstract:Based on the ambient air quality monitoring data of Xi′an from September to November in 2021,the pollution characteristics of regional PM2.5 in Xi′an before,during and after the 14th National Games of the People′s Republic of China & the National Games for Persons with Disabilities were analysed.In addition,the air mass transport pathways,potential source areas and their contributions to PM2.5 mass concentrations were also determined by the Hybrid Single-Particle Lagrangian Integrated Trajectory(HYSPLIT) model,potential source contributing factor PSCF) method and weighted concentration trajectory CWT) method.The results showed that during the observation,Southeasterly winds dominated in Xi′an,and the PM2.5 concentration was lowest during the 14th National Games of the People′s Republic of China NG) and highest after the end of the National Games for Persons with Disabilities PG).The air mass transport pathways of PM2.5 were dominated by the local transport near Southern Shaanxi,the proximity transport in the Easterly and Southerly directions,and the long-range transport in the Northwesterly direction.Potential source regions of PM2.5 were mainly located in the cities and areas of Henan,Hubei,Guizhou,Sichuan,Chongqing,Xinjiang Uygur Autonomous Regions,Inner Mongolia Autonomous Region,Southern Shaanxi,and the junction of cities in Southern Shaanxi and neighboring provinces.

    Key words:the 14th National Games of the People′s Republic of China; fine particulate matter; backward trajectory; potential source analysis; concentration weight analysis

    0 引言

    隨著我國城市化和工業(yè)化進(jìn)程飛速發(fā)展,大氣環(huán)境問題日趨嚴(yán)峻,根據(jù)國家環(huán)保部發(fā)布的2019-2021年生態(tài)環(huán)境公報,西安市已連續(xù)3年成為全國空氣質(zhì)量相對較差的城市[1],主要污染物為PM2.5,已被認(rèn)為是區(qū)域性污染物[2].區(qū)域大氣污染已經(jīng)成為當(dāng)前社會發(fā)展面臨的重大環(huán)境問題之一[3].

    西安市大氣污染兼具本地污染累積、周邊跨區(qū)域傳輸和不利氣象條件等共同作用的特征[4].目前,混合單粒子拉格朗日綜合軌跡(Hybrid Single-Particle Lagrangian Integrated Trajectory,HYSPLIT)模型已被廣泛用于區(qū)域污染物分布、擴散和傳輸特征研究[5-10].污染物潛在源分析可利用后向軌跡確定污染源的地理位置及空間分布,研究主要集中在我國長三角[8,11,12]、珠三角[13-15]和京津冀[16-18]等重點大氣污染區(qū)域,對于中部內(nèi)陸地區(qū)顆粒物潛在源分析的研究卻較為鮮見.對西安市溯源的研究,主要集中在污染物時空分布特征[19,20]、成分分析[21-26]、氣象影響因素分析[27-31]等方面,研究區(qū)域氣團(tuán)傳輸對大氣污染特征的影響以及顆粒物潛在源分析的研究還存在明顯不足.

    當(dāng)前,我國舉辦大型峰會等活動的數(shù)量與日俱增.活動期間,主辦方通常會采取聯(lián)防聯(lián)控等臨時性管控措施保障環(huán)境空氣質(zhì)量[32].杭州市G20峰會管控期間空氣質(zhì)量總體優(yōu)良,PM2.5污染物濃度同比下降40%[33];APEC會議期間北京地區(qū)的平均風(fēng)速和相對濕度優(yōu)于會前會后[34].

    為保障十四運會及殘?zhí)貖W會期間的空氣質(zhì)量,西安市各區(qū)貫徹執(zhí)行《西安市大氣污染防治條例》,以污染物濃度、氣象條件以及模型模擬為主的大氣環(huán)境質(zhì)量分析是評估重大活動舉辦期間管控有效性的重要方法.基于此,本文通過收集十四運及殘?zhí)貖W會前中后期PM2.5環(huán)境監(jiān)測數(shù)據(jù)和氣象要素等常規(guī)觀測數(shù)據(jù),利用HYSPLIT模型、潛在源貢獻(xiàn)因子法(Potential Source Contribution Function,PSCF)和權(quán)重濃度軌跡法(Concentration Weighted Trajectory,CWT)等方法,探討了西安市PM2.5的傳輸途徑和潛在來源,定量識別不同區(qū)域之間的傳輸貢獻(xiàn),以期為揭示西部區(qū)域大氣污染分布特征、來源解析提供數(shù)據(jù)支撐.

    1 實驗部分

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

    研究時期為十四運會和殘?zhí)貖W會觀測期間(2021年9月2日~11月6日),分為6個階段,分別是十四運之前(9月2日~14日,NGB)、十四運期間(9月15日~27日,NGD)、十四運之后(9月28日~10月10日,NGA)、殘?zhí)貖W會之前(10月14日~21日,PGB)、殘?zhí)貖W會期間(10月20日~29日,PGD)和殘?zhí)貖W會之后(2021年10月30日~11月6日,PGA).觀測期間6種主要空氣污染物PM2.5、PM10、SO2、NO2、CO和O3的逐小時質(zhì)量濃度數(shù)據(jù)來自西安市氣象數(shù)據(jù)網(wǎng)以及中國環(huán)境監(jiān)測總站的全國城市空氣質(zhì)量數(shù)據(jù)發(fā)布平臺(https://www.aqistudy.cn/historydata/).風(fēng)速、氣溫和降水量的日均數(shù)據(jù)來自國家氣象科學(xué)數(shù)據(jù)中心(http://data.cma.cn).后向軌跡模型數(shù)據(jù)來自美國國家預(yù)報中心的格式為“grib2”的“FNL全球分析資料”和全球資料同化系統(tǒng)(ftp://arlftp.arlhq.noaa.gov/pub/archives/gdas1),分辨率為0.25°×0.25°[35-37].

    1.2 HYSPLIT模型

    本研究利用HYSPLIT模型模擬西安市顆粒物的輸送軌跡、擴散路徑和來源[38-40].模擬高度選取500 m[41,42],西安市市區(qū)(34.27°N,108.93°E)作為目標(biāo)模擬區(qū)域,時間分辨率為1 h,進(jìn)行每日24個時次UTC的72 h后向軌跡計算,以反映出污染物的輸送特征.本研究使用MeteoInfo軟件及TrajStat插件中的聚類方法對觀測期間的氣流軌跡進(jìn)行聚類分析,整個觀測期間共產(chǎn)生1 512條軌跡.

    1.3 潛在源貢獻(xiàn)因子分析法

    PSCF利用氣流軌跡確定污染區(qū)域及空間分布,將氣流軌跡覆蓋的污染源區(qū)全部劃為0.25°×0.25°水平的網(wǎng)格(i,j)[43-47],定量描述研究期間每個網(wǎng)格概率場的軌跡來源強度,識別對污染物濃度影響大的區(qū)域.PSCF即是經(jīng)過網(wǎng)格的污染軌跡(Xij)與所有區(qū)域氣流總軌跡(Yij)之比,計算公式如式(1)所示:

    PSCFij=XijYij(1)

    PSCF是一種條件概率,值大代表網(wǎng)格污染軌跡數(shù)多,該網(wǎng)格對應(yīng)區(qū)域的污染源多.為了減小分母Yij引起的不確定性,引入權(quán)重函數(shù)Wij來修正PSCF值,用WPSCF表示,詳見公式(2)所示:

    WPSCF=Wij×PSCF(2)

    Wij函數(shù)定義如公式(3)[48,49]所示:

    Wnij)=[JB{][HL2:1,;2,Y]1.0080

    1.4 權(quán)重濃度軌跡分析法

    PSCF只能反映污染軌跡經(jīng)過每個網(wǎng)格的概率,即潛在源區(qū)貢獻(xiàn)率大小,不能確定研究區(qū)域的污染程度.而CWT則可以定量估算每個網(wǎng)格上污染物的濃度.基于代表經(jīng)過某網(wǎng)格單元(i,j)的空氣團(tuán)對西安市觀測點空氣質(zhì)量的貢獻(xiàn)權(quán)重的CWTij值,可以揭示出西安市PM2.5濃度網(wǎng)格化地域分布規(guī)律.某網(wǎng)格單元(i,j)內(nèi)的CWTij值可用以下公式進(jìn)行計算[38,50]:

    CWTij=1∑Ml=1τijl∑Ml=1Clτijl(4)

    式 4) 中:l和M分別表示網(wǎng)格單元編號和軌跡數(shù)量;C表示軌跡l到達(dá)觀測點時大氣空氣質(zhì)量的濃度水平,μg·m-3;τijl表示軌跡l在網(wǎng)格單元(i,j)內(nèi)停留的時間.本研究設(shè)置CWT網(wǎng)格單元的大小為0.25°×0.25°.

    CWT法不會因總軌跡個數(shù)少而存在不確定性,在評估遠(yuǎn)距離傳輸貢獻(xiàn)時,網(wǎng)格點內(nèi)軌跡的距離節(jié)點減少,產(chǎn)生誤差,也要通過權(quán)重函數(shù)Wij消除少量端點網(wǎng)格的不確定性,降低Cij誤差,修正后的CWT用WCWT表示,詳見公式(5)所示:

    WCWT=Wij×CWT(5)

    2 結(jié)果與討論

    2.1 污染特征分析

    十四運及殘?zhí)貖W會觀測期間西安市主要空氣污染物質(zhì)量濃度及氣象因素日均變化見圖1所示.PM2.5和PM10的平均質(zhì)量濃度呈現(xiàn)先降低后持續(xù)升高的趨勢,十四運會期間PM2.5和PM10的平均質(zhì)量濃度顯著低于其他時期,分別為12.79±9.44 μg·m-3和21.53±15.89 μg·m-3,表明十四運會期間空氣污染控制措施效果顯著.殘?zhí)貖W會結(jié)束后PM2.5和PM10的平均質(zhì)量濃度達(dá)到最高值,分別為61.75±33.57 μg·m-3和112.51±33.91 μg·m-3,11月4日PM2.5和PM10的日均濃度最高,分別達(dá)107 μg·m-3和156 μg·m-3,污染物濃度的升高可能與短期污染控制措施的停止有關(guān).

    PM2.5/PM10代表細(xì)顆粒物占可吸入顆粒物之比,幾乎均>0.5,殘?zhí)貖W會之后為0.54±0.19,與往年同期PM2.5在PM10中占比一致[51].PM2.5和PM10濃度分別有3天和2天超過《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)(GB 3095-2012)》二級標(biāo)準(zhǔn).研究期間PM2.5平均濃度為53.84±44.38 μg·m-3,低于袁夢娜[52]2018年9~11月西安市PM2.5濃度66.7 μg·m-3,表明十四運會期間各項空氣污染最高級別保障措施效果顯著,之后污染物濃度的升高可能與短期污染控制措施的停止有關(guān).

    基于SPSS 22.0對氣象參數(shù)與污染物濃度進(jìn)行Pearson相關(guān)性分析,如圖2所示.SO2、NO2及CO是典型的污染源識別物.整個觀測期間西安市SO2、NO2、CO濃度變化與顆粒物呈顯著正相關(guān),而SO2、NO2、CO主要來自燃煤燃燒、機動車尾氣排放[53]以及碳基燃料的不完全燃燒[54],說明PM2.5、PM10的濃度受燃燒源和機動車尾氣排放源的影響較大.

    觀測期間不同時期西安市大氣污染物濃度變化如圖3所示.SO2和NO2平均濃度在十四運會期間最低,分別為5.38±1.11 μg·m-3和17.99±5.11 μg·m-3,在殘?zhí)貖W會后最高,是最低值的1.38和2.81倍.CO濃度整體呈上升趨勢,在殘?zhí)貖W會期間最高,濃度為0.80±0.11 mg·m-3.研究期間NO2濃度整體上升,O3濃度基本下降,十四運會之前O3濃度最高,為66.63±30.14 μg·m-3,可能是O3的前體物質(zhì)增多而促進(jìn)O3生成[55].

    2.2 氣象特征分析

    整個觀測期間西安市主要氣象因素均值如表1所示.期間主要風(fēng)向為東南方向,以東南偏東方向為主,平均溫度17.96±5.29 ℃,平均相對濕度68.3±15.1%,風(fēng)速2~21 m·s-1,存在明顯波動.整體平均風(fēng)速較低,為4~5級,十四運和殘?zhí)貖W會前期的風(fēng)速相對較高,平均風(fēng)速分別為8.63±4.20和8.36±4.04 m·s-1,大氣水平運動較強,空氣流動較快,有利于污染物擴散.

    溫度影響污染物垂直擴散,觀測期間溫度整體呈下降趨勢,十四運會之前平均氣溫最高,為24.86±2.99 ℃.相關(guān)性分析表明溫度與PM10在0.05水平上呈負(fù)相關(guān),溫度與O3濃度顯著正相關(guān),Pearson系數(shù)為0.665,溫度與O3的回歸分析如圖4所示.殘?zhí)貖W會之后進(jìn)入11月,平均氣溫低,為14.08±4.24 ℃,大氣較穩(wěn)定,空氣在地表聚集不宜垂直擴散,不利于污染物稀釋,導(dǎo)致污染物濃度較高;溫度低,光照少,光化學(xué)反應(yīng)變?nèi)酰琌3含量相對變少.

    降雨影響顆粒物的大氣停留時間,圖1表明研究期間西安市相對濕度偏高,整體先增大后降低,在9月17日和23日達(dá)到最大值96%,表1說明十四運會期間平均相對濕度達(dá)82.7±15.9%.十四運會期間降雨充足,9月15~27日,10月3~9日發(fā)生連續(xù)性降雨,水汽條件好,雨水蒸騰促進(jìn)氣體交換,降低污染物濃度,降雨時PM2.5與PM10濃度分別降低2.88倍和6.55倍,起到顯著疏散稀釋與凈化潔凈作用.相關(guān)性分析表明相對濕度與PM2.5(R=-0.250)、PM10(R=-0.430)、NO2(R=-0.475)呈負(fù)相關(guān)關(guān)系,均通過0.05水平顯著性檢驗,表明降雨及較大的相對濕度可有效沖刷大氣中顆粒物及NO2,利于污染物稀釋擴散.PM2.5中20%~80% 的有機物是水溶性的[56],相對濕度高更易于NH3、NOx等氣體污染物吸附于顆粒物表面,降雨時隨雨滴沉降至近地面,濃度降低.而殘?zhí)貖W會之后進(jìn)入11月,天氣干燥降雨量少,水汽條件差,缺少雨水對大氣的清潔作用,污染物難聚集沉降,濃度較高.

    研究期間風(fēng)向風(fēng)速如圖5所示,結(jié)合圖1,觀測期間的主要風(fēng)向較穩(wěn)定,從東北、東方和西南偏南方向往西安市匯入污染物,以東北向為主.平均風(fēng)速4~5級,11月6日風(fēng)速最大,為21 m·s-1,10月6日風(fēng)速最小,為2 m·s-1,十四運會和殘?zhí)貖W會期間平均風(fēng)速低,分別為7.85±3.0 m·s-1和6.77±3.16 m·s-1,氣流穩(wěn)定,污染傳輸少;十四運會之前和殘?zhí)貖W會之前及之后平均風(fēng)速高,分[JP3]別為8.63±4.20、8.36±4.04和8.59±5.63 m·s-1,風(fēng)力等級高,空氣流動快大氣水平運動強,利于污染物擴散及匯入,外來氣團(tuán)攜帶污染物多,導(dǎo)致西安市大氣組分更復(fù)雜.相關(guān)性分析表明風(fēng)速與O3呈正相關(guān),風(fēng)將顆粒物吹散,光照和太陽輻射加強,有利于O3生成.同時,平流層底與對流層頂進(jìn)行物質(zhì)交換,高濃度的O3向下傳輸累積,冷鋒過境使CO、NO2減少.風(fēng)速與CO、NO2呈負(fù)相關(guān),對氣態(tài)污染物的濃度有稀釋作用.

    2.3 污染物傳輸特征分析

    2.3.1 氣團(tuán)后向軌跡及聚類分析

    依據(jù)十四運及殘?zhí)貖W會觀測期間各時期PM2.5濃度軌跡的空間分布進(jìn)行聚類,得到3~6個主要代表路徑及其垂直傳輸高度.觀測期間各時期西安市后向軌跡聚類分析結(jié)果及其對應(yīng)的污染物濃度如圖6和表2所示.

    十四運會之前各氣團(tuán)對應(yīng)的PM2.5濃度均值為軌跡1>3>6>4>2>5,PM2.5/PM10低.本地PM2.5排放較少,受傳輸影響較大.軌跡1、3、6的PM2.5濃度較高,為十四運會之前PM2.5的主要輸送路徑,占總軌跡的51.3%,主要集中于陜南,來自貴州、重慶北部,軌跡1和軌跡3路徑短,傳輸速[HJ1.7mm]度慢.軌跡1沿陜西南部經(jīng)漢中安康抵達(dá)西安市,占比28.5%,軌跡3沿江蘇、安徽北部、河南抵達(dá)西安市,占比16.7%.軌跡5沿著新疆北部途經(jīng)內(nèi)蒙古西部、寧夏銀川市經(jīng)銅川抵達(dá)西安市,占比16.4%,屬于長距離傳輸,對應(yīng)氣團(tuán)移動速度快,擴散條件好,較清潔.十四運會期間氣團(tuán)軌跡對應(yīng)的PM2.5均值3>4>1>2.來自貴州的軌跡3和來自新疆的軌跡4氣團(tuán)攜帶的污染物濃度高,屬于遠(yuǎn)距離輸送.來自東部地區(qū)河南北部、山西南部的軌跡1和經(jīng)過甘肅中部、內(nèi)蒙南部、寧夏后輸送到西安市的軌跡2都屬于短距離運輸,分別占總軌跡的48.4%和23.1%,污染物濃度低,相對清潔.十四運會期間管控措施最為嚴(yán)格且多雨,導(dǎo)致PM2.5濃度在所有時期中最低.十四運會之后氣團(tuán)對應(yīng)的PM2.5濃度均值為軌跡4>1>3>2.來自陜西商洛附近的軌跡1占比高達(dá)60.6%,屬于本地運輸,是十四運會之后的主要輸送通道,由于路徑較短,對應(yīng)的PM2.5濃度相對較高,該氣團(tuán)長時間停留在西安市,加劇了顆粒物的本地累積;此外,位于陜南地區(qū)的秦嶺山脈阻擋了部分氣團(tuán)的移動和限制了污染物的擴散,導(dǎo)致本地污染較為嚴(yán)重.軌跡3攜帶的顆粒物濃度較低,傳輸路徑長.軌跡4為短距離運輸,由貴州、經(jīng)重慶四川到達(dá)西安市,說明大氣環(huán)境較為穩(wěn)定,有利于顆粒物累積.

    殘?zhí)貖W會觀測期間西安市后向軌跡聚類分析、占比及污染物平均質(zhì)量濃度如圖7和表3所示.殘?zhí)貖W會之前西安市大部分氣團(tuán)來自西北方向的遠(yuǎn)距離傳輸和東南偏東方向的短距離傳輸,分別占總軌跡的53.7%和46.4%,其中來自西北方向的軌跡又可分為軌跡3和軌跡1,途徑空氣較為清潔的新疆、甘肅、內(nèi)蒙古、寧夏和銅川地區(qū),攜帶的污染物少,加上風(fēng)速較大,利于污染物的擴散和稀釋.來自河南和陜西東部的軌跡2氣團(tuán)對應(yīng)的PM2.5濃度最高,表明河南和陜西東部是殘?zhí)貖W會之前西安市高濃度顆粒物的重要污染源區(qū).

    殘?zhí)貖W會期間,氣團(tuán)軌跡主要分6類,均以短距離傳輸為主,詳見圖7(b).其中,軌跡2氣團(tuán)源于陜西延安、經(jīng)河南和陜西商洛,后進(jìn)入西安市,占比28.7%,軌跡4來源于陜西漢中,占比22.4%,軌跡1源于內(nèi)蒙古,經(jīng)陜西南部一路輸送到西安市,占比20.3%.軌跡1、2和4氣團(tuán)所對應(yīng)的PM2.5濃度相對較高,可能是的原因是較低的風(fēng)速和穩(wěn)定的氣象條件抑制了污染物擴散和稀釋,進(jìn)而造成污染物積聚.來自西北方向的軌跡3占比少,僅為7.3%,攜帶的PM2.5濃度相對較低.殘?zhí)貖W會之后,PM2.5濃度顯著高于殘?zhí)貖W會開始前后,達(dá)61.75±33.57 μg·m-3.根據(jù)圖7(e),殘?zhí)貖W會結(jié)束后氣團(tuán)主要來自東南方向的河南和陜西渭南商洛.其中軌跡1氣團(tuán)在總軌跡中占比最多,為32.8%,從河南西部出發(fā),經(jīng)湖北北部、商洛地區(qū),從東南方向到達(dá)西安市,傳輸距離短,風(fēng)速較小,不有利于顆粒物污染的擴散,其對應(yīng)的PM2.5平均濃度最高,為97.19 μg·m-3,其次是軌跡6.軌跡2占比23.4%,對應(yīng)PM2.5平均濃度為52.27 μg·m-3,其輸送距離同樣較短,從西安北部到達(dá)銅川,又經(jīng)渭南轉(zhuǎn)回西安市,該氣團(tuán)在關(guān)中地區(qū)有所停滯,存在明顯的繞流,是關(guān)中地區(qū)二次污染發(fā)生的重要因素[57].

    2.3.2 污染氣團(tuán)潛在源區(qū)分析

    依據(jù)PM2.5二級標(biāo)準(zhǔn)限值75 μg·m-3,利用PSCF分析計算十四運及殘?zhí)貖W會觀測期間不同時期PM2.5的WPSCF值.不同時期的WPSCF分布存在顯著差異,WPSCF值越大表明該區(qū)域越有可能成為西安市污染物的潛在源區(qū).圖8表明十四運會前WPSCF在江蘇安徽河南中部湖北西北部、陜南、山西南部和甘肅寧夏的交界處最高,其WPSCF>0.6,結(jié)果表明,除了當(dāng)?shù)匚廴驹磁欧磐?,這些地區(qū)的氣團(tuán)輸送和污染物傳播對西安市大氣污染也有很大影響.十四運會期間及之后污染格網(wǎng)數(shù)量相當(dāng),分散著零星PM2.5潛在源高值區(qū),整體上WPSCF值較低.十四運會期間陜西省南部、四川重慶交界處及貴州南部的WPSCF>0.5.十四運會之后WPSCF>0.6的區(qū)域主要在貴州.

    殘?zhí)貖W會前WPSCF均在0.5左右,其高值區(qū)分布較少,在陜南和湖北北部團(tuán)狀集中分布.陜南安康商洛及周邊城市交接處如湖北十堰WPSCF>0.5,進(jìn)一步驗證了周邊污染物傳輸影響西安市大氣污染.殘?zhí)貖W會期間WPSCF高值區(qū)顯著增加,顆粒物污染的主要源區(qū)是西安市的鄰近區(qū)域,包括陜西南部、山西南部、河南西部、湖北北部以及四川重慶交界區(qū)域,WPSCF>0.6.殘?zhí)貖W會后污染格網(wǎng)數(shù)量少,位于東南偏東方向的河南和商洛線是PM2.5高值潛在源區(qū),主要集中在陜西南部、山西南部、河南西部、湖北北部,WPSCF>0.6.陜西南部、湖北西北部、四川省與重慶市交匯處、貴州南部等地區(qū)是西安市十四運會觀測期間顆粒物主要潛在源區(qū).殘?zhí)貖W會觀測期間顆粒物的主要源區(qū)是陜西省南部、山西南部、河南西部和湖北北部等偏東區(qū)域.

    2.3.3 顆粒物權(quán)重濃度軌跡分析

    利用CWT法定量不同污染源區(qū)對西安市PM2.5的貢獻(xiàn)情況,WCWT高代表貢獻(xiàn)大.圖9表明觀測期間西安市PM2.5的WCWT高值區(qū)呈點狀或團(tuán)狀分布,十四運會之前主要在陜南安康附近以及重慶北部,WCWT值在30 μg·m-3左右.西北和東部區(qū)域WCWT貢獻(xiàn)較低,小于9 μg·m-3.十四運會期間WCWT最低,大于20 μg·m-3的區(qū)域集中在貴州、四川重慶交界地區(qū)、陜西南部、河南西部和湖北北部.十四運會后WCWT高值區(qū)域主要分布在貴州地區(qū),WCWT>30 μg·m-3.

    殘?zhí)貖W會之前WCWT高值區(qū)主要在陜西南部、重慶市北部、四川東北部、河南西部和湖北北部,其中,漢中、巴中、洛陽和十堰等地區(qū)是WCWT>35 μg·m-3的主要貢獻(xiàn)區(qū)域.殘?zhí)貖W會期間河南西部、湖北北部、陜南和四川重慶交界地區(qū)均為西安市的潛在貢獻(xiàn)源區(qū),PM2.5濃度貢獻(xiàn)達(dá)到30~60 μg·m-3.而殘?zhí)貖W會之后,WCWT高值區(qū)域主要集中在陜西省附近,呈東南方向帶狀分布,其中WCWT>100 μg·m-3的地區(qū)主要集中在陜西南部、河南西部和湖北北部.

    3 結(jié)論

    本文以十四運及殘?zhí)貖W會為例,探討重大會議、賽事等活動期間的管控減排措施和氣象條件與污染物濃度變化的內(nèi)在聯(lián)系,并基于HYSPLIT、PSCF和CWT等方法精確識別污染物的潛在貢獻(xiàn)源區(qū),以期為相關(guān)部門后續(xù)研究區(qū)域聯(lián)防聯(lián)控下的污染物排放規(guī)律和碳減排策略的制定提供理論和實際數(shù)據(jù)支撐.得到的主要結(jié)論如下:

    (1)十四運會期間PM2.5的平均質(zhì)量濃度最低,為12.79±9.44 μg·m-3,說明十四運會期間各項空氣污染控制措施效果顯著;殘?zhí)貖W會結(jié)束后PM2.5濃度顯著高于其他時期,可能與短期污染控制措施的停止和不利的氣象條件有關(guān).

    (2)基于后向軌跡聚類分析,觀測期間的氣團(tuán)主要來自遠(yuǎn)距離的新疆內(nèi)蒙地區(qū)、近距離的河南湖北貴州四川重慶地區(qū)和陜西南部.

    (3)基于PSCF和CWT的分析,觀測期間PM2.5的潛在貢獻(xiàn)源區(qū)主要集中在陜南各地市周邊省份交界地帶,且濃度水平受偏南和偏東區(qū)域傳輸?shù)挠绊戄^大.

    參考文獻(xiàn)

    [1] 中華人民共和國生態(tài)環(huán)境部.2019-2021中國生態(tài)環(huán)境狀況公報[DB/OL].https://www.mee.gov.cn/hjzl/sthjzk/,2022-12-20.

    [2] 李莎莎.利用HYSPLIT模式探析西安市大氣顆粒物的輸送路徑.智能城市,2017,35):62-65.

    [3] San José R,Pérez J L,Pérez L,et al.Effects of climate change on the health of citizens modelling urban weather and air pollution.Energy,2018,165:53-62.

    [4] 王 芳,陳東升,程水源,等.基于氣流軌跡聚類的大氣污染輸送影響.環(huán)境科學(xué)研究,2009,226):637-642.

    [5] Liao T T,Wang S,Ai J,et al.Heavy pollution episodes,transport pathways and potential sources of PM2.5 during the winter of 2013 in Chengdu China).Science of the Total Environment,2017,584-585:1 056-1 065.

    [6] Wang Y Q,Wang M M,Li S P,et al.Study on the oxidation potential of the water-soluble components of ambient PM2.5 over Xi′an,China:Pollution levels,source apportionment and transport pathways.Environment International,2020,136:105 515.

    [7] Draxler R R.NOAA Technical memorandum ERL ARL-224 description of the hysplit 4 modeling system[R].Silver Spring Maryland:Air Resources Laboratory,1997:1-22.

    [8] 王 茜.利用軌跡模式研究上海大氣污染的輸送來源.環(huán)境科學(xué)研究,2013,264):357-363.

    [9] hao M F,Huang? S,Qiao T,et al.Chemical characterization,the transport pathways and potential sources of PM2.5 in Shanghai:Seasonal variations.Atmospheric Research,2015,158:66-78.

    [10] 高 陽,韓永貴,黃曉宇,等.基于后向軌跡模式的豫南地區(qū)冬季PM2.5來源分布及傳輸分析.環(huán)境科學(xué)研究,2021,343):538-548.

    [11] Wang Y F,Wang H Y,hang S H.A weighted higher order network analysis of fine particulate matter PM2.5) transport in Yangtze River Delta.Physica A:Statistical Mechanics and its Applications,2018,49615):654-662.

    [12] 李 敏,何紅弟,郝楊楊,等.上海市大氣環(huán)境中PM2.5/PM10時空分布特征.云南大學(xué)學(xué)報:自然科學(xué)版,2019,412):323-332.

    [13] 潘月云,陳多宏,葉斯琪,等.廣東省大氣污染典型案例特征及其影響因素分析.安全與環(huán)境工程,2017,242):58-66.

    [14] 黃 俊,廖碧婷,王春林,等.廣州市PM2.5污染特征及潛在貢獻(xiàn)源區(qū)分析.環(huán)境科學(xué)與技術(shù),2019,424):115-124.

    [15] 李婷苑,鄧雪嬌,范紹佳,等.2010年廣州亞運期間空氣質(zhì)量與污染氣象條件分析.環(huán)境科學(xué),2012,339):2 932-2 938.

    [16] 沈洪艷,呂總璞,師華定,等.基于HYSPLIT模型的京津冀地區(qū)大氣污染物輸送的路徑分析.環(huán)境工程技術(shù)學(xué)報,2018,84):359-366.

    [17] 鈐偉妙,張艷品,陳 靜,等.石家莊大氣污染物輸送通道及污染源區(qū)研究.環(huán)境科學(xué)學(xué)報,2018,389):3 438-3 448.

    [18] hang? Y,Wong M S,Lee? H.Estimation of potential source regions of PM2.5 in Beijing using backward trajectories.Atmospheric Pollution Research,2015,61):173-177.

    [19] 于英翠,李利敏,左亞杰,等.西安市大氣污染物時空分布特征及其趨勢評估.環(huán)境工程,2018,3610):165-169.

    [20] 李杏茹,白 羽,陳 曦,等.北京冬季重污染過程大氣細(xì)顆粒物化學(xué)組成特征及來源分析.環(huán)境化學(xué),2018,3711):2 397-2 409.

    [21] 曹 寧,黃學(xué)敏,祝 穎,等.西安冬季重污染過程PM2.5理化特征及來源解析.中國環(huán)境科學(xué),2019,391):32-39.

    [22] Thamban N M,Lalchandani V,umar V,et al.Evolution of size and composition of fine particulate matter in the Delhi megacity during later winter.Atmospheric Environment,2021,267:118 752.

    [23] 高庚申,遲 峰,毛金群,等.貴陽市大氣細(xì)顆粒物中多環(huán)芳烴時空分布特征.環(huán)境監(jiān)控與預(yù)警,2019,114):43-47.

    [24] 謝運興,唐 曉,郭宇宏,等.新疆大氣顆粒物的時空分布特征.中國環(huán)境監(jiān)測,2019,351):26-36.

    [25] 王 成,閆雨龍,謝 凱,等.陽泉市秋冬季PM2.5化學(xué)組分及來源分析.環(huán)境科學(xué),2020,413):1 036-1 044.

    [26] 吳 丹,張立平,夏廣鋒,等.沈陽市大氣中PM2.5來源解析研究.環(huán)境監(jiān)測管理與技術(shù),2020,321):18-22.

    [27] Wang J,Xing J,Mathur R,et al.Historical trends in PM2.5-related premature mortality during 1990-2010 across the Northern Hemisphere.Environment Health Perspective,2017,1253):400-408.

    [28] 李星諭,毛 瑤,陳展樂,等.華中地區(qū)冬季灰霾天氣下PM2.5中重金屬污染特征及健康風(fēng)險評價:以湖北黃岡為例.環(huán)境科學(xué),2021,4210):4 593-4 601.

    [29] Mai X W,hou H F,Li Y Y,et al.Associations between ambient fine particulate PM2.5) exposure and cardiovascular disease:Findings from the China Health and Retirement Longitudinal Study CHARLS).Environmental Science & Pollution Research International,2022,29:1-8.

    [30] Bianca S M,Virginia C,Hwanmi L,et al.Similar polycyclic aromatic hydrocarbon and genotoxicity profiles of atmospheric particulate matter from cities on three different continents.Environmental and Molecular Mutagenesis,2020,615):560-573.

    [31] Sánchez Piero J,Moreda Pieiro J,Concha Graa E,et al.Inhalation bioaccessibility estimation of polycyclic aromatic hydrocarbons from atmospheric particulate matter PM10):Influence of PM10 composition and health risk assessment.Chemosphere,2021,263:127 847.

    [32] 伏晴艷,李健軍,田旭東,等.我國重大活動環(huán)境空氣質(zhì)量監(jiān)測預(yù)報評估技術(shù)體系的發(fā)展與思考.中國環(huán)境監(jiān)測,2020,362):1-9.

    [33] 周德榮,田旭東,蔡 哲,等.G20峰會期間長三角區(qū)域臭氧變化及管控效果評估.中國環(huán)境監(jiān)測,2020,362):41-49.

    [34] 張礁石,陸亦懷,桂華僑,等.APEC會議前后北京地區(qū)PM2.5污染特征及氣象影響因素.環(huán)境科學(xué)研究,2016,295):646-653.

    [35] 楊燕燕.甘肅西北關(guān)鍵城市空氣污染特征及顆粒物潛在源區(qū)研究[D].蘭州:蘭州大學(xué),2020.

    [36] hao N,Wang G,Li G H,et al.Air pollution episodes during the COVID-19 outbreak in the Beijing-Tianjin-Hebei region of China:An insight into the transport pathways and source distribution.Environmental Pollution,2020,267:115 617

    [37] hao ,Liu R,hang .Characteristics of winter haze pollution in the Fenwei Plain and the possible influence of EU during 1984-2017.Earth and Space Science,2020,76):1-13.

    [38] 黃光球,雷 哲.西安市大氣顆粒物PM2.5的輸送路徑和潛在源分析.云南大學(xué)學(xué)報自然科學(xué)版),2019,416):1 191-1 200.

    [39] Draxler R R,Hess G D.An overview of the HYSPLIT-4 modelling system for trajectories,dispersion and deposition.Australian Meteorological Magazine,1998,474):295-308.

    [40] Stein A F,Draxler R R,Rolph G D,et al.NOAA′s HYSPLIT atmospheric transport and dispersion modeling system.Bulletin of the American Meteorological Society,2015,96:2 059-2 077.

    [41] 鄭 悅,程 方,張 凱,等.保定市大氣污染特征和潛在輸送源分析.環(huán)境科學(xué)研究,2020,332):260-270.

    [42] 段時光,姜 楠,楊留明,等.鄭州市冬季大氣PM2.5傳輸路徑和潛在源分析.環(huán)境科學(xué),2019,401):86-93.

    [43] Cheng M D,Hopke P ,Barrie L,et al.Qualitative determination of source regions of aerosol in Canadian high arctic.Environmental Science & Technology,1993,2720):63-71.

    [44] Begum B A,im E,Jeong C H,et al.Evaluation of the potential source contribution function using the 2002 Quebec forest fire episode.Atmospheric Environment,2005,39:3 719 -3 724.

    [45] Michael L A,Linb C J,Peter M,et al.Atmospheric mercury near Salmon Falls Creek Reservoir in southern Idaho.Applied Geochemistry,2008,233):438-453.

    [46] Ferhat ,Ismail A,Omar A.Long-range potential source contributions of episodic aerosol events to PM10 profile of a megacity.Atmospheric Environment,2009,4336):5 713-5 722.

    [47] Salvador P,Artinano B,Alonso D G,et al.Identification and characterization of sources of PM10 in Madrid Spain) by statistical methods.Atmospheric Environment,2004,38:435-447.

    [48] Xu X,Akhtar U S.Identification of potential regional sources of atmospheric total gaseous mercury in Windsor,Ontario,Canada using hybrid receptor modeling.Atmospheric Chemistry and Physics,2010,1015):7 073-7 083.

    [49] Polissar A V,Hopke P ,Harris J M.Source regions for atmospheric aerosol measured at Barrow,Alaska.Environmental Science & Technology,2001,3521):4 214-4 226.

    [50]Bhuyan P ,Bharali C,Pathak B,et al.The role of precursor gases and meteorology on temporal evolution of O3 at a tropical location in northeast India.Environmental Science & Pollution Research International,2014,2110):6 696-6 713.

    [51] 李婷苑,鄧雪嬌,范紹佳,等.2010年廣州亞運期間空氣質(zhì)量與污染氣象條件分析.環(huán)境科學(xué),2012,339):2 932-2 938.

    [52] 袁夢娜.西安大氣碳?xì)馊苣z光學(xué)性質(zhì)及其來源解析[D].西安:中國科學(xué)院大學(xué)中國科學(xué)院地球環(huán)境研究所),2022.

    [53] Jion M M M F,Jannat J N,Mia M Y,et al.A critical review and prospect of NO2 and SO2 pollution over Asia:Hotspots,trends,and sources.Science of the Total Environment,2023,876:162 851.

    [54] Ling J X,hang G,Liu L X,et al.Observation of atmospheric CO2 and CO in a low-carbon pilot city:Insight into CO2 sources and regional transport.Atmospheric Pollution Research,2022,137):101 585.

    [55] 楊 帆.降雨對大氣顆粒物和氣態(tài)污染物的清除效率及機制[D].南昌:南昌大學(xué),2015.

    [56] 顧澤平.大氣細(xì)顆粒物有機質(zhì)組成的變化規(guī)律及其在源解析中的應(yīng)用[D].上海:上海大學(xué),2010.

    [57] 李培榮,姚 靜,張煦庭,等.西安冬季不同環(huán)流形勢下大氣細(xì)顆粒物的輸送特征研究.陜西氣象,20224):7-14.

    【責(zé)任編輯:蔣亞儒】

    基金項目:陜西省科技廳自然科學(xué)基礎(chǔ)研究計劃項目(2022JM-076)

    作者簡介:李香凝(1999—),女,河北辛集人,在讀碩士研究生,研究方向:大氣環(huán)境化學(xué)

    通訊作者:王羽琴(1988—),女,陜西榆林人,副教授,博士,研究方向:大氣環(huán)境化學(xué),wangyuqin@sust.edu.cn

    av福利片在线观看| 级片在线观看| 久久99蜜桃精品久久| 我的老师免费观看完整版| av福利片在线观看| 国产精品久久久久久久电影| 插逼视频在线观看| 老熟妇乱子伦视频在线观看| 欧美日本视频| 啦啦啦韩国在线观看视频| 中文在线观看免费www的网站| 亚洲在线自拍视频| 一级av片app| 国产亚洲欧美98| 欧美zozozo另类| 亚洲精品国产成人久久av| 成人漫画全彩无遮挡| 一级毛片电影观看 | 亚洲一区高清亚洲精品| 日本欧美国产在线视频| 两个人视频免费观看高清| 女人十人毛片免费观看3o分钟| 国产高清有码在线观看视频| 日韩三级伦理在线观看| 亚洲精品成人久久久久久| 99久久久亚洲精品蜜臀av| 在线观看av片永久免费下载| 麻豆久久精品国产亚洲av| 内射极品少妇av片p| 内射极品少妇av片p| 国产白丝娇喘喷水9色精品| 伦精品一区二区三区| 久久鲁丝午夜福利片| www.av在线官网国产| 亚洲成人av在线免费| 看片在线看免费视频| 欧美激情久久久久久爽电影| 日韩欧美在线乱码| 亚洲成av人片在线播放无| 国产色爽女视频免费观看| 精品久久久久久成人av| 九九爱精品视频在线观看| 精品久久久久久久久久久久久| 欧美高清性xxxxhd video| 一本精品99久久精品77| 麻豆精品久久久久久蜜桃| 最近2019中文字幕mv第一页| 日本一二三区视频观看| 国产精品美女特级片免费视频播放器| 能在线免费观看的黄片| av在线播放精品| 在线观看美女被高潮喷水网站| 超碰av人人做人人爽久久| 少妇裸体淫交视频免费看高清| 亚洲国产欧洲综合997久久,| 久久精品综合一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 免费看av在线观看网站| 春色校园在线视频观看| 亚洲无线在线观看| 综合色丁香网| 99久久无色码亚洲精品果冻| 禁无遮挡网站| 中国美女看黄片| 亚洲欧美日韩无卡精品| 99久久精品一区二区三区| 日韩视频在线欧美| 日韩一区二区视频免费看| 一个人看视频在线观看www免费| 白带黄色成豆腐渣| 国产亚洲精品久久久久久毛片| 身体一侧抽搐| 黄片无遮挡物在线观看| 亚洲自拍偷在线| 亚洲精品日韩av片在线观看| 三级国产精品欧美在线观看| 免费观看a级毛片全部| 五月伊人婷婷丁香| 久久亚洲精品不卡| 黑人高潮一二区| 床上黄色一级片| 国产一区二区三区在线臀色熟女| av视频在线观看入口| 欧美性猛交╳xxx乱大交人| 中文字幕熟女人妻在线| 欧美zozozo另类| 国模一区二区三区四区视频| 午夜爱爱视频在线播放| 免费在线观看成人毛片| 亚洲国产精品国产精品| 一卡2卡三卡四卡精品乱码亚洲| 亚洲久久久久久中文字幕| 成人综合一区亚洲| 亚洲天堂国产精品一区在线| 国产探花极品一区二区| 国产爱豆传媒在线观看| 少妇丰满av| 亚洲国产欧洲综合997久久,| 日韩三级伦理在线观看| 日本爱情动作片www.在线观看| 18禁裸乳无遮挡免费网站照片| 嫩草影院入口| 日韩av在线大香蕉| 精品人妻熟女av久视频| 偷拍熟女少妇极品色| 日本一二三区视频观看| 人妻系列 视频| 中文在线观看免费www的网站| 欧美日韩国产亚洲二区| 人人妻人人澡人人爽人人夜夜 | 九九热线精品视视频播放| 99riav亚洲国产免费| 久久精品国产亚洲网站| 亚洲精品乱码久久久v下载方式| 看黄色毛片网站| 国产老妇女一区| 2022亚洲国产成人精品| 国产精品久久久久久亚洲av鲁大| 欧美日韩在线观看h| 国产成人精品婷婷| 免费看美女性在线毛片视频| 亚洲高清免费不卡视频| 国产人妻一区二区三区在| 99热这里只有是精品在线观看| 国产av一区在线观看免费| 国产伦精品一区二区三区四那| av在线亚洲专区| 欧美最新免费一区二区三区| 波多野结衣高清无吗| 在线观看一区二区三区| 日韩欧美三级三区| 听说在线观看完整版免费高清| 一区二区三区四区激情视频 | 国产伦精品一区二区三区四那| 特级一级黄色大片| 亚洲第一区二区三区不卡| 国产精品无大码| 长腿黑丝高跟| 成人漫画全彩无遮挡| 午夜精品在线福利| 麻豆精品久久久久久蜜桃| 一本精品99久久精品77| 啦啦啦韩国在线观看视频| 亚洲最大成人中文| 国产精品久久久久久久电影| 久久99蜜桃精品久久| 五月玫瑰六月丁香| 国产黄片视频在线免费观看| 欧美国产精品一级二级三级| 毛片一级片免费看久久久久| av不卡在线播放| 久久综合国产亚洲精品| 久久久久久久精品精品| 婷婷色综合www| 精品人妻熟女av久视频| 伊人久久国产一区二区| 久久久欧美国产精品| 美女脱内裤让男人舔精品视频| 男女边摸边吃奶| 搡女人真爽免费视频火全软件| 99久久中文字幕三级久久日本| 人人澡人人妻人| 精品国产国语对白av| 国产亚洲一区二区精品| 极品少妇高潮喷水抽搐| 国产成人午夜福利电影在线观看| 亚洲成色77777| 男女啪啪激烈高潮av片| 草草在线视频免费看| 日产精品乱码卡一卡2卡三| 久久久午夜欧美精品| 久久久久久伊人网av| 黑人猛操日本美女一级片| 丰满少妇做爰视频| 亚洲精品成人av观看孕妇| 又黄又爽又刺激的免费视频.| 亚洲精品美女久久av网站| 亚洲精品日韩在线中文字幕| 桃花免费在线播放| 中文字幕最新亚洲高清| 老司机亚洲免费影院| 97在线人人人人妻| 成年人免费黄色播放视频| 一级,二级,三级黄色视频| 婷婷色综合大香蕉| 各种免费的搞黄视频| 少妇丰满av| 一本色道久久久久久精品综合| 亚洲欧洲国产日韩| 午夜久久久在线观看| 久久97久久精品| 精品久久久久久久久亚洲| 制服丝袜香蕉在线| 能在线免费看毛片的网站| 两个人的视频大全免费| 国产成人精品久久久久久| 永久网站在线| 亚洲伊人久久精品综合| 日韩强制内射视频| 在线 av 中文字幕| 老熟女久久久| 国产精品.久久久| 美女主播在线视频| 黑人猛操日本美女一级片| 国产亚洲精品久久久com| 三级国产精品片| 国产亚洲精品第一综合不卡 | 亚洲精品久久成人aⅴ小说 | 在线看a的网站| 桃花免费在线播放| 久久韩国三级中文字幕| 老司机影院成人| 蜜桃久久精品国产亚洲av| 亚洲性久久影院| 日韩人妻高清精品专区| 啦啦啦在线观看免费高清www| 国产精品免费大片| 嫩草影院入口| 国产一级毛片在线| 久久亚洲国产成人精品v| 亚洲欧美一区二区三区黑人 | 欧美老熟妇乱子伦牲交| 国产黄色免费在线视频| 亚洲精品乱久久久久久| 国产成人免费无遮挡视频| 国产精品蜜桃在线观看| 一个人看视频在线观看www免费| 国产成人精品无人区| 大码成人一级视频| 亚洲欧美清纯卡通| 丝袜在线中文字幕| 熟女av电影| 国产黄色免费在线视频| 少妇人妻 视频| 97在线视频观看| 亚洲欧美清纯卡通| 亚洲av在线观看美女高潮| av有码第一页| 又大又黄又爽视频免费| 国产国拍精品亚洲av在线观看| 2021少妇久久久久久久久久久| 九九爱精品视频在线观看| a级片在线免费高清观看视频| 美女大奶头黄色视频| 日本vs欧美在线观看视频| 啦啦啦中文免费视频观看日本| 少妇 在线观看| 成年女人在线观看亚洲视频| 制服丝袜香蕉在线| 少妇人妻 视频| av在线播放精品| 国产一级毛片在线| 久久鲁丝午夜福利片| 久久久午夜欧美精品| 99九九线精品视频在线观看视频| 啦啦啦啦在线视频资源| 热99久久久久精品小说推荐| 亚洲精品自拍成人| 又粗又硬又长又爽又黄的视频| 国产精品一区www在线观看| 新久久久久国产一级毛片| 美女脱内裤让男人舔精品视频| 桃花免费在线播放| 亚洲综合精品二区| 97在线人人人人妻| 国产熟女欧美一区二区| 日韩中文字幕视频在线看片| 人妻人人澡人人爽人人| 国产精品麻豆人妻色哟哟久久| 国产熟女欧美一区二区| 亚洲av男天堂| 中文字幕精品免费在线观看视频 | 日韩视频在线欧美| 人成视频在线观看免费观看| 亚洲国产欧美在线一区| 中文乱码字字幕精品一区二区三区| 黄色欧美视频在线观看| 寂寞人妻少妇视频99o| 久久久精品94久久精品| 国产69精品久久久久777片| 蜜桃在线观看..| 边亲边吃奶的免费视频| 亚州av有码| 简卡轻食公司| 少妇被粗大的猛进出69影院 | 国产成人一区二区在线| av在线观看视频网站免费| 午夜免费鲁丝| 亚洲国产成人一精品久久久| 亚洲中文av在线| 精品熟女少妇av免费看| 男女边吃奶边做爰视频| 九色亚洲精品在线播放| 熟女av电影| 婷婷色综合大香蕉| 久久99热这里只频精品6学生| 久久久久久久国产电影| 久久av网站| 中文天堂在线官网| 国产高清国产精品国产三级| 中文字幕亚洲精品专区| 国产高清不卡午夜福利| 亚洲内射少妇av| 久久99热6这里只有精品| 亚洲精华国产精华液的使用体验| 这个男人来自地球电影免费观看 | 插阴视频在线观看视频| videos熟女内射| 成人毛片60女人毛片免费| 伦理电影免费视频| 亚洲精品中文字幕在线视频| 亚洲性久久影院| 婷婷色综合大香蕉| 最近手机中文字幕大全| 日本欧美视频一区| 男女边摸边吃奶| 久久99热这里只频精品6学生| 亚洲精品一区蜜桃| 国产av码专区亚洲av| 国产片内射在线| 51国产日韩欧美| 日韩欧美一区视频在线观看| 人成视频在线观看免费观看| 少妇的逼好多水| 韩国av在线不卡| 国产日韩欧美在线精品| 亚洲精品自拍成人| 国产成人a∨麻豆精品| 久久久久久久久久久丰满| 国产白丝娇喘喷水9色精品| 日韩伦理黄色片| 亚洲婷婷狠狠爱综合网| 高清毛片免费看| 欧美三级亚洲精品| 久久影院123| 国产精品蜜桃在线观看| 99视频精品全部免费 在线| 丰满乱子伦码专区| 少妇熟女欧美另类| 91精品三级在线观看| 色吧在线观看| 国产在视频线精品| 两个人免费观看高清视频| 国产精品一国产av| 18禁在线无遮挡免费观看视频| 精品亚洲成国产av| 看十八女毛片水多多多| 成人国语在线视频| 亚洲欧美中文字幕日韩二区| 少妇的逼水好多| 亚洲精品日韩av片在线观看| 精品卡一卡二卡四卡免费| 成人二区视频| 国产探花极品一区二区| 九九久久精品国产亚洲av麻豆| 热99国产精品久久久久久7| 免费黄频网站在线观看国产| 久久久久久久久久久久大奶| 18禁在线播放成人免费| 麻豆乱淫一区二区| 欧美+日韩+精品| av福利片在线| 精品少妇内射三级| 99久久人妻综合| 欧美日韩一区二区视频在线观看视频在线| 国产精品人妻久久久影院| 亚洲av综合色区一区| 五月开心婷婷网| 日韩三级伦理在线观看| 人妻人人澡人人爽人人| 大又大粗又爽又黄少妇毛片口| 亚洲国产精品一区二区三区在线| 精品一区在线观看国产| 嘟嘟电影网在线观看| 全区人妻精品视频| 一级二级三级毛片免费看| 一级爰片在线观看| 肉色欧美久久久久久久蜜桃| 2018国产大陆天天弄谢| 国产乱来视频区| 欧美日韩综合久久久久久| 国产一区亚洲一区在线观看| 亚洲五月色婷婷综合| 大香蕉久久网| 久久精品久久久久久噜噜老黄| 春色校园在线视频观看| 国产av一区二区精品久久| 午夜av观看不卡| 婷婷色综合www| 久久久国产欧美日韩av| 亚洲精品亚洲一区二区| 少妇猛男粗大的猛烈进出视频| 日本欧美视频一区| 欧美bdsm另类| 亚洲精品乱码久久久久久按摩| 2018国产大陆天天弄谢| 亚洲五月色婷婷综合| 精品视频人人做人人爽| 黄片无遮挡物在线观看| 热re99久久国产66热| 国产成人午夜福利电影在线观看| 美女主播在线视频| 熟女人妻精品中文字幕| 亚洲性久久影院| 桃花免费在线播放| 日韩一区二区三区影片| 国产有黄有色有爽视频| 啦啦啦中文免费视频观看日本| 国产色婷婷99| 又大又黄又爽视频免费| av不卡在线播放| 亚洲无线观看免费| 色吧在线观看| 午夜激情福利司机影院| 岛国毛片在线播放| 女人精品久久久久毛片| 老司机亚洲免费影院| 免费大片黄手机在线观看| 亚洲内射少妇av| 亚洲五月色婷婷综合| 少妇高潮的动态图| 一本色道久久久久久精品综合| 蜜桃国产av成人99| 大话2 男鬼变身卡| 久久97久久精品| 热99久久久久精品小说推荐| 日韩制服骚丝袜av| 欧美xxxx性猛交bbbb| 波野结衣二区三区在线| 欧美少妇被猛烈插入视频| 纯流量卡能插随身wifi吗| 午夜福利网站1000一区二区三区| 中文字幕av电影在线播放| 精品人妻在线不人妻| av免费观看日本| 国产黄色视频一区二区在线观看| videosex国产| 国产综合精华液| 亚洲三级黄色毛片| 少妇丰满av| a级毛色黄片| 国产精品久久久久久久电影| 五月玫瑰六月丁香| 国产69精品久久久久777片| 精品酒店卫生间| 亚洲国产av新网站| 少妇 在线观看| videossex国产| 国产伦理片在线播放av一区| 色婷婷av一区二区三区视频| 麻豆精品久久久久久蜜桃| 满18在线观看网站| 91久久精品国产一区二区三区| 99久久人妻综合| 国产精品一国产av| 人妻少妇偷人精品九色| 女性生殖器流出的白浆| 丰满迷人的少妇在线观看| 亚洲国产毛片av蜜桃av| 亚洲av成人精品一二三区| 亚洲欧美日韩另类电影网站| 精品熟女少妇av免费看| 中文字幕av电影在线播放| 中国美白少妇内射xxxbb| 久久久久久久久大av| 欧美丝袜亚洲另类| 天堂中文最新版在线下载| 人人妻人人添人人爽欧美一区卜| 国产高清三级在线| 精品熟女少妇av免费看| 亚洲av电影在线观看一区二区三区| 少妇精品久久久久久久| 如日韩欧美国产精品一区二区三区 | 日韩一本色道免费dvd| 高清欧美精品videossex| 日韩精品有码人妻一区| 丁香六月天网| 国产精品偷伦视频观看了| 99热全是精品| 99久久人妻综合| 亚洲成人av在线免费| 亚洲高清免费不卡视频| 久久精品国产亚洲av涩爱| 精品午夜福利在线看| 超碰97精品在线观看| 91国产中文字幕| 亚洲人与动物交配视频| 国产成人精品久久久久久| 午夜老司机福利剧场| 亚洲性久久影院| 青春草亚洲视频在线观看| 久久毛片免费看一区二区三区| 国产亚洲午夜精品一区二区久久| 免费少妇av软件| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产成人一精品久久久| 国产黄频视频在线观看| 少妇 在线观看| 国产精品 国内视频| 欧美日韩av久久| av.在线天堂| 丰满少妇做爰视频| 国产精品一区二区在线观看99| 91在线精品国自产拍蜜月| 国产伦理片在线播放av一区| 黑人欧美特级aaaaaa片| 日韩在线高清观看一区二区三区| 97在线视频观看| 精品人妻在线不人妻| 国产成人精品婷婷| 国产亚洲精品第一综合不卡 | 亚州av有码| 一个人免费看片子| 一本—道久久a久久精品蜜桃钙片| 免费黄网站久久成人精品| 婷婷色综合大香蕉| 成年av动漫网址| 久久久久精品久久久久真实原创| 亚洲欧美一区二区三区国产| 交换朋友夫妻互换小说| 欧美xxxx性猛交bbbb| 22中文网久久字幕| a 毛片基地| 91午夜精品亚洲一区二区三区| 久久久久网色| 最近中文字幕2019免费版| 久久国产亚洲av麻豆专区| 欧美激情极品国产一区二区三区 | 精品少妇内射三级| 成年人免费黄色播放视频| 天美传媒精品一区二区| 欧美激情国产日韩精品一区| 黄色怎么调成土黄色| 国产成人freesex在线| 国产日韩一区二区三区精品不卡 | 街头女战士在线观看网站| 嘟嘟电影网在线观看| 久久精品国产自在天天线| 内地一区二区视频在线| av播播在线观看一区| 成人二区视频| 欧美人与性动交α欧美精品济南到 | av国产精品久久久久影院| 免费黄频网站在线观看国产| 99热网站在线观看| 欧美 日韩 精品 国产| 成人手机av| 在现免费观看毛片| 人妻夜夜爽99麻豆av| 中文乱码字字幕精品一区二区三区| 午夜福利网站1000一区二区三区| 亚洲精品日本国产第一区| 纵有疾风起免费观看全集完整版| 国产黄色视频一区二区在线观看| 熟女电影av网| 成人漫画全彩无遮挡| 一个人免费看片子| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av在线观看美女高潮| 十分钟在线观看高清视频www| 国产高清三级在线| 国产精品人妻久久久影院| 亚洲精品久久成人aⅴ小说 | 天堂俺去俺来也www色官网| 亚洲av中文av极速乱| 啦啦啦啦在线视频资源| 国产精品无大码| 午夜精品国产一区二区电影| 日日撸夜夜添| 亚洲图色成人| 女人精品久久久久毛片| 亚洲三级黄色毛片| 亚洲国产欧美在线一区| 丝袜喷水一区| 国产亚洲欧美精品永久| 亚洲av成人精品一区久久| av又黄又爽大尺度在线免费看| 在线观看免费视频网站a站| 久久久国产一区二区| 欧美一级a爱片免费观看看| 欧美少妇被猛烈插入视频| 在线播放无遮挡| 国产亚洲最大av| 日韩成人伦理影院| 欧美日韩视频高清一区二区三区二| 国产av一区二区精品久久| 久久午夜综合久久蜜桃| 欧美日本中文国产一区发布| 嫩草影院入口| 纵有疾风起免费观看全集完整版| 韩国av在线不卡| 国产黄色免费在线视频| 日本黄色片子视频| 亚洲av成人精品一二三区| 久久久久久久亚洲中文字幕| 亚洲国产av新网站| 亚洲第一av免费看| 2018国产大陆天天弄谢| 日日摸夜夜添夜夜添av毛片| 久久99热这里只频精品6学生| 久久女婷五月综合色啪小说| 婷婷色综合www| 亚洲精品日韩在线中文字幕| 女的被弄到高潮叫床怎么办| av免费在线看不卡| 国产av一区二区精品久久| 久久99热这里只频精品6学生| 伊人久久精品亚洲午夜| 王馨瑶露胸无遮挡在线观看| 岛国毛片在线播放| 最近手机中文字幕大全| 亚洲激情五月婷婷啪啪| 国产无遮挡羞羞视频在线观看| 国产欧美日韩综合在线一区二区| 亚洲av在线观看美女高潮|