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

    氣象因子對(duì)臭氧的影響及其在空氣質(zhì)量預(yù)報(bào)中的應(yīng)用

    2017-09-27 12:35:35蔣昌潭羅財(cái)紅劉姣姣韓世剛
    關(guān)鍵詞:壓差溫差氣壓

    張 燦,蔣昌潭,羅財(cái)紅,劉姣姣,葉 堤,譙 捷,韓世剛

    1.重慶市生態(tài)環(huán)境監(jiān)測(cè)中心,重慶 401147 2.城市大氣環(huán)境綜合觀測(cè)與污染防控重慶市重點(diǎn)實(shí)驗(yàn)室,重慶 401147 3.重慶市環(huán)境監(jiān)察總隊(duì),重慶 401147 4.重慶市氣象服務(wù)中心,重慶 401147

    氣象因子對(duì)臭氧的影響及其在空氣質(zhì)量預(yù)報(bào)中的應(yīng)用

    張 燦1,2,蔣昌潭1,2,羅財(cái)紅1,2,劉姣姣1,2,葉 堤1,2,譙 捷3,韓世剛4

    1.重慶市生態(tài)環(huán)境監(jiān)測(cè)中心,重慶 401147 2.城市大氣環(huán)境綜合觀測(cè)與污染防控重慶市重點(diǎn)實(shí)驗(yàn)室,重慶 401147 3.重慶市環(huán)境監(jiān)察總隊(duì),重慶 401147 4.重慶市氣象服務(wù)中心,重慶 401147

    為提高重慶市臭氧(O3)預(yù)報(bào)準(zhǔn)確率,利用2013—2015年5—10月O3監(jiān)測(cè)數(shù)據(jù)和氣象數(shù)據(jù),通過主成分分析、逐步回歸分析等方法,確定了影響重慶O3濃度的主要?dú)庀笠蛩貫樽罡邷囟?、溫差、太?yáng)輻射、降水量、相對(duì)濕度、水氣壓和壓差;通過基于O3污染水平相似的主要?dú)庀罂刂埔蜃雍Y選和最優(yōu)組合的預(yù)報(bào)結(jié)果優(yōu)化方法,提高了O3預(yù)報(bào)準(zhǔn)確率,使2016年5—8月O3的AQI類別預(yù)報(bào)準(zhǔn)確率由57.7%增至72.4%,O3超標(biāo)的預(yù)報(bào)準(zhǔn)確率由38%增至46%。

    氣象因子;臭氧;預(yù)報(bào)

    環(huán)境空氣質(zhì)量預(yù)報(bào)已經(jīng)成為監(jiān)測(cè)部門一項(xiàng)比較重要的例行工作[1]。我國(guó)目前的環(huán)境空氣污染表現(xiàn)為以大氣細(xì)顆粒物(PM2.5)和臭氧(O3)為代表的復(fù)合型大氣污染,高濃度O3對(duì)眼睛和呼吸道有很強(qiáng)的刺激性,損害人體肺功能和傷害農(nóng)作物,并使大氣能見度降低[2]。O3已經(jīng)成為各大城市除PM2.5之外的造成環(huán)境空氣質(zhì)量超標(biāo)的首要污染物,因此對(duì)O3的準(zhǔn)確預(yù)報(bào)顯得尤為重要。

    目前的預(yù)報(bào)方法主要有數(shù)值模型和統(tǒng)計(jì)模型。統(tǒng)計(jì)模型只需氣象數(shù)據(jù)、污染物監(jiān)測(cè)濃度數(shù)據(jù)即可進(jìn)行預(yù)報(bào)[3],運(yùn)算量小、易于操作,不需要源清單和高性能計(jì)算硬件設(shè)備[4],避免了數(shù)值模型預(yù)報(bào)方法的弊端[3]。與數(shù)值模型相比,統(tǒng)計(jì)模型預(yù)報(bào)結(jié)果準(zhǔn)確率高[5],適合預(yù)報(bào)剛剛起步的城市使用。重慶市從2014年10月開始對(duì)公眾發(fā)布空氣質(zhì)量預(yù)報(bào)結(jié)果,利用的是法蘭德斯技術(shù)研究院 (VITO)與比利時(shí)環(huán)保署合作開發(fā)的基于神經(jīng)網(wǎng)絡(luò)的OVL模型。2014—2015年統(tǒng)計(jì)結(jié)果表明,模型對(duì)O3的AQI類別的預(yù)報(bào)準(zhǔn)確率為67.5%,低于其對(duì)PM2.5的預(yù)報(bào)準(zhǔn)確率73.9%;對(duì)O3超標(biāo)的判斷準(zhǔn)確率僅為36.4%,遠(yuǎn)低于PM2.5的67.9%。重慶市空氣質(zhì)量超標(biāo)天數(shù)均由PM2.5和O3造成,提高O3濃度預(yù)報(bào)準(zhǔn)確率尤其是準(zhǔn)確預(yù)判O3超標(biāo)天數(shù)將在很大程度上提高預(yù)報(bào)準(zhǔn)確率。

    近幾年,北京、天津、上海、廈門、珠三角區(qū)域、沈陽(yáng)、西寧、太原、臨安本底站等都針對(duì)氣象條件對(duì)O3影響進(jìn)行了研究[6-16];但對(duì)O3預(yù)報(bào)方法的研究較少,安俊琳等[6]利用主成分和回歸分析方法對(duì)北京O3進(jìn)行了預(yù)報(bào);劉明花等[9]將多元非線性預(yù)報(bào)模式應(yīng)用在上海地面O3濃度的預(yù)報(bào)中;宋榕榮等[10]利用多元線性回歸法建立了廈門市O3預(yù)報(bào)及評(píng)估系統(tǒng),都得到了相關(guān)性較好的預(yù)報(bào)結(jié)果。重慶也進(jìn)行了O3的相關(guān)研究[17],但目前尚沒有針對(duì)重慶O3預(yù)報(bào)的研究。

    重慶下墊面類型復(fù)雜,高溫高濕,太陽(yáng)輻射強(qiáng)度大,與其他城市不同。統(tǒng)計(jì)模型預(yù)報(bào)不需要源清單,更多的是依賴氣象條件的變化,在不同的天氣類型下,污染物擴(kuò)散方式和濃度分布均有明顯差異,天氣類型對(duì)區(qū)域空氣質(zhì)量有關(guān)鍵作用。由于各個(gè)氣象因子并不是獨(dú)立的而是相互影響、相互制約和協(xié)同作用的,僅單個(gè)分析氣象因素對(duì)O3濃度的影響難免造成信息缺失或誤導(dǎo)。本文在利用主成分分析和逐步回歸分析分析確定了O3的主控氣象因素后,提出了基于O3污染水平相似的主要?dú)庀罂刂埔蜃雍Y選和最優(yōu)組合的預(yù)報(bào)結(jié)果優(yōu)化方法,對(duì)OVL模型進(jìn)行補(bǔ)充,以期提高O3預(yù)報(bào)準(zhǔn)確率,為環(huán)境管理和污染天氣預(yù)警提供及時(shí)準(zhǔn)確的信息。

    1 實(shí)驗(yàn)部分

    1.1數(shù)據(jù)來源和處理

    由于重慶O3污染主要集中在5—10月,因此,O3濃度數(shù)據(jù)利用2013—2015年5—10月主城區(qū)逐日自動(dòng)監(jiān)測(cè)數(shù)據(jù);監(jiān)測(cè)儀器為河北EC9810B O3自動(dòng)監(jiān)測(cè)儀,采用紫外光度法原理。逐日氣象數(shù)據(jù)來自重慶市氣象臺(tái),包括雨量、最高溫度、最低溫度、相對(duì)濕度、最小濕度、水氣壓、氣壓、最低氣壓、最高氣壓、風(fēng)速、風(fēng)向;溫差、壓差由原始數(shù)據(jù)計(jì)算而得,有降水標(biāo)記為“1”,無(wú)降水記為“0”。太陽(yáng)輻射數(shù)據(jù)來自重慶市大氣環(huán)境監(jiān)測(cè)超級(jí)站,監(jiān)測(cè)儀器為荷蘭UVS-AB紫外輻射計(jì),UVS-AB紫外輻射計(jì)通過紫外輻射傳感器和CR1000數(shù)據(jù)采集模塊實(shí)時(shí)測(cè)量透過大氣層的太陽(yáng)紫外輻射UV-A和UV-B,監(jiān)測(cè)到達(dá)地面的太陽(yáng)紫外線輻射強(qiáng)度,共521個(gè)有效樣本。

    1.2分析方法

    1.2.1 主成分分析

    主成分分析(PCA)是對(duì)一組變量進(jìn)行降維的統(tǒng)計(jì)學(xué)方法,即用少量的獨(dú)立變量來代替大量的有相關(guān)性的變量的一種方法?,F(xiàn)在PCA方法中普遍采用最大方差正交旋轉(zhuǎn)法來得到有意義的因子[18-19]。PCA的具體計(jì)算過程大體分3步:奇異值分解;維數(shù)的選擇(一般選擇特征值大于1的提取原則);因子旋轉(zhuǎn)(一般采用方差最大正交旋轉(zhuǎn))[6]。

    1.2.2 逐步回歸分析

    多元線性方程的建立一般采用逐步回歸算法。逐步回歸算法是在所有考慮的參數(shù)中,按其對(duì)因變量Y(大氣污染物預(yù)報(bào)濃度)影響的顯著程度的大小,由大到小逐個(gè)引進(jìn)回歸方程。若已被引進(jìn)回歸方程的參數(shù),在引進(jìn)新參數(shù)后,可能會(huì)由顯著變?yōu)椴伙@著,此時(shí)需將其從回歸方程中剔除,以保證在眾多預(yù)報(bào)參數(shù)中挑選出最佳的關(guān)鍵參數(shù)組合,建立最優(yōu)預(yù)報(bào)方程。污染物預(yù)報(bào)濃度Y與預(yù)報(bào)關(guān)鍵參數(shù)X建立的最優(yōu)回歸方程格式:

    式中:Y為污染物預(yù)報(bào)濃度;BQ為常數(shù)項(xiàng);X1,X2,…,Xn為預(yù)報(bào)關(guān)鍵參數(shù);B1,B2,Bn為關(guān)鍵參數(shù)的系數(shù)[4]。劉閩等[20]采用逐步回歸方法建立了沈陽(yáng)市冬季環(huán)境空氣質(zhì)量統(tǒng)計(jì)預(yù)報(bào)模型。

    2 結(jié)果與討論

    2.1主要影響因子

    2.1.1 主成分分析

    將太陽(yáng)最大輻射、降水量、最高溫度、濕度、最小濕度、水氣壓、最高氣壓、風(fēng)速、風(fēng)向、壓差、溫差、是否下雨等12個(gè)參數(shù)作為變量,進(jìn)行主成分分析,使用方差極大法進(jìn)行旋轉(zhuǎn),選擇KMO(Kaiser-Meyer-Olkin)和Bartlett的球形度檢驗(yàn),得到旋轉(zhuǎn)后的因子載荷矩陣(見表1),Bartlett球檢驗(yàn)顯著性指標(biāo)sig.<0.05為顯著,說明存在因子結(jié)構(gòu)。KMO>0.9,說明非常適合因子分析;0.8

    旋轉(zhuǎn)后的因子載荷矩陣共產(chǎn)生3個(gè)貢獻(xiàn)較大的因子,其特征值均大于1,說明具有統(tǒng)計(jì)顯著性。相關(guān)研究表明,選取因子時(shí)一般選取載荷大于0.50的因子[18]。載荷為負(fù)值時(shí)說明該因子與本項(xiàng)目呈負(fù)相關(guān)關(guān)系。第一個(gè)因子中的高載荷組分主要指示了相對(duì)濕度、溫差、是否下雨,其次為最大太陽(yáng)輻射、雨量、最高溫度,這些因子解釋了36.2%的因變量變化;第二個(gè)因子中的高載荷組分主要指示了水氣壓、最高氣壓,這些因子解釋了21.7%的因變量變化;第三個(gè)因子中的高載荷組分主要指示了壓差、風(fēng)速,這些因子解釋了14.5%的因變量變化。這些因子共解釋了72.4%的因變量,說明最大輻射、降水量、最高溫度,水氣壓、最高氣壓,壓差、風(fēng)速是影響5—10月氣象形勢(shì)的主要因素。

    表1 旋轉(zhuǎn)后的因子載荷矩陣

    用同樣的方法對(duì)O3質(zhì)量濃度小于或等于160 μg/m3的樣本(O3的AQI≤100,O3達(dá)標(biāo)日)以及O3質(zhì)量濃度大于160 μg/m3的樣本(O3的AQI>100,O3超標(biāo)日)分別進(jìn)行主成分分析,結(jié)果見表2、表3。O3達(dá)標(biāo)日產(chǎn)生3個(gè)主要因子,解釋了73.5%的變量,反映了相對(duì)濕度、最高溫度、溫差、降水、水氣壓、最高氣壓、壓差和風(fēng)速等對(duì)O3達(dá)標(biāo)日氣象形勢(shì)的影響,與5—10月相似。

    O3超標(biāo)日產(chǎn)生5個(gè)主要因子:第一個(gè)因子中的高載荷組分主要指示了最高溫度、相對(duì)濕度、最高氣壓;第二個(gè)因子中的高載荷組分主要指示了最小濕度、水氣壓、溫差;第三個(gè)因子中的高載荷組分主要指示了降水量;第四個(gè)因子中的高載荷組分主要指示了最大輻射;第五個(gè)因子中的高載荷組分主要指示了壓差。這5個(gè)因子共解釋了78.7%的變量,反映了最高溫度、溫差、相對(duì)濕度、最高氣壓、壓差、水氣壓、降水量和最大太陽(yáng)輻射對(duì)O3超標(biāo)日氣象形勢(shì)的影響。與O3達(dá)標(biāo)日有所不同,多了“最大太陽(yáng)輻射”的影響,而少了“風(fēng)速”的影響。

    表2 O3達(dá)標(biāo)日旋轉(zhuǎn)后的因子載荷矩陣

    表3 O3超標(biāo)日旋轉(zhuǎn)后的因子載荷矩陣

    2.1.2 逐步回歸分析

    對(duì)521個(gè)樣本經(jīng)過逐步回歸分析,因變量為O3濃度,自變量為雨量、最高溫度、最低溫度、相對(duì)濕度、最小濕度、水氣壓、氣壓、最低氣壓、最高氣壓、風(fēng)速、風(fēng)向、溫差、壓差、是否下雨;利用Durbin-Watson值(DW,取值0~4)檢驗(yàn)回歸模型中殘差獨(dú)立的假設(shè),當(dāng)DW愈接近2時(shí),殘差項(xiàng)間越不相關(guān)。本次計(jì)算得到的DW=1.238,可以基本確定殘差是獨(dú)立的。

    回歸結(jié)果如表4所示,共得到影響O3濃度的4個(gè)主要參數(shù),即溫差、最高溫度、是否下雨和水氣壓。膨脹因子VIF<10,容差大于0.1,說明這4個(gè)變量間不存在共線性,都可用于方程。sig.均小于0.05,具有顯著意義。根據(jù)表4得到O3濃度的逐步回歸方程式:

    式中:Y為O3質(zhì)量濃度,μg/m3;X1為溫差,℃;X2為最高溫度,℃;X3為是否下雨,下雨為“1”,不下雨為“0”;X4為水氣壓,hPa。

    表4 逐步回歸結(jié)果

    注:“/”表示共線性僅針對(duì)變量進(jìn)行統(tǒng)計(jì)。

    由圖1可知,通過該方程得到的O3濃度與實(shí)際值呈顯著正相關(guān)關(guān)系,相關(guān)系數(shù)達(dá)到0.86,通過了α=0.01水平(雙側(cè))的顯著性驗(yàn)證。通過對(duì)521個(gè)數(shù)據(jù)進(jìn)行驗(yàn)證,該方程對(duì)O3的AQI類別的預(yù)報(bào)準(zhǔn)確率為73.9%,雖然對(duì)優(yōu)、良水平的預(yù)報(bào)準(zhǔn)確率較高,但對(duì)O3超標(biāo)的預(yù)報(bào)準(zhǔn)確率僅為24.6%。

    圖1 O3濃度擬合值與實(shí)際值的關(guān)系Fig.1 Fitting value and actual value of O3 concentration

    通過以上方法,已經(jīng)基本確定了重慶主城區(qū)影響O3生成的主要?dú)庀笠蛩貫樽罡邷囟?、溫差、太?yáng)輻射、降水量、相對(duì)濕度、水氣壓、壓差(見表5)。

    表5 不同方法提取出的氣象因子

    注:“√”表示影響較大;“空”表示沒有影響或者影響不顯著。

    安俊琳等[6]分析了北京太陽(yáng)紫外輻射、溫度、相對(duì)濕度和壓力對(duì)夏季O3濃度的影響,表明氣溫和相對(duì)濕度是主要影響因素,其次是紫外輻射。北京超標(biāo)日O3濃度與氣壓、相對(duì)濕度呈負(fù)相關(guān)關(guān)系,與風(fēng)速、溫度呈正相關(guān)關(guān)系[7]。任麗紅等[8]對(duì)北京夏季不同高度O3濃度與氣象條件的研究表明:風(fēng)速在2 m/s以下,近地面O3濃度隨著風(fēng)速增加而增加;風(fēng)速超過2 m/s,O3濃度隨風(fēng)速的變化比較復(fù)雜,其原因有待進(jìn)一步研究。臨安大氣本底站研究也表明,O3與風(fēng)向的關(guān)系較復(fù)雜,與局地生成的O3有關(guān)[16]。劉明花[9]研究發(fā)現(xiàn),上海一年四季對(duì)近地面O3影響最大的氣象因素是太陽(yáng)輻射。廈門市O3污染與溫度、風(fēng)速、相對(duì)濕度等氣象條件相關(guān),且不同季節(jié)氣象因子的影響程度各不相同[10]。沈陽(yáng)風(fēng)速和溫度對(duì)O3濃度影響較大[11]。天津影響O3濃度的主要?dú)庀笠蛩厥菤鉁?、相?duì)濕度和風(fēng)速風(fēng)向[12]。太原市城區(qū)O3與溫度、濕度和風(fēng)速存在極顯著的相關(guān)性[13]。西寧地區(qū)地面O3濃度隨云量狀況、風(fēng)向、風(fēng)速、相對(duì)濕度、氣壓等要素的變化比較明顯[15]。石家莊O3濃度與氣溫呈正相關(guān),與相對(duì)濕度呈負(fù)相關(guān)[21]。從各個(gè)城市的研究結(jié)果來看,溫度、相對(duì)濕度均是影響O3濃度的主要因素;其他因素也因?yàn)榈赜驓庀筇卣?、地形特征的不同而有所不同?/p>

    2.2因子最優(yōu)組合

    圖2給出了氣象因素在不同分段值內(nèi)O3的AQI范圍。

    圖2 O3的AQI在不同氣象因子下的分布Fig.2 Distribution of O3 AQI under meteorological factors

    從圖2可以看出:

    1)最高溫度:溫度在30 ℃及其以上時(shí),出現(xiàn)O3超標(biāo),且主要集中在33~39 ℃;27 ℃及以下,O3的AQI基本為0~50。

    2)溫差:溫差在8 ℃及以下時(shí),O3不會(huì)超標(biāo);溫差在8 ℃以上時(shí),開始出現(xiàn)O3超標(biāo);溫差在6 ℃及以下時(shí),O3的AQI基本為0~50;溫差為10~12 ℃,容易出現(xiàn)高濃度O3,AQI>150。

    3)太陽(yáng)輻射:太陽(yáng)輻射小于或等于40 W/m2時(shí),O3達(dá)標(biāo);太陽(yáng)輻射大于40 W/m2時(shí),開始出現(xiàn)O3超標(biāo),且主要集中在40~70 W/m2;高濃度O3,AQI>150主要出現(xiàn)在40(不含40)~60之間;輻射更高時(shí)(如70 W/m2以上),O3超標(biāo)率反而不高。

    4)降水量:O3超標(biāo)主要發(fā)生在無(wú)降水時(shí);降水量在20 mm以上時(shí),O3的AQI基本為0~50。

    5)相對(duì)濕度:相對(duì)濕度在80%及以下時(shí),都可能出現(xiàn)O3超標(biāo),主要集中在40%~70%;高濃度O3出現(xiàn)在50%~70%;相對(duì)濕度大于80%時(shí),O3的AQI為0~50。

    6)水氣壓:水氣壓在任何范圍內(nèi),都可能出現(xiàn)O3超標(biāo),只是水氣壓在15~20 hPa時(shí),超標(biāo)概率較小。

    7)壓差:壓差在4 hPa以下時(shí),O3基本不會(huì)超標(biāo);壓差在4 hPa以上時(shí),開始出現(xiàn)O3超標(biāo);超標(biāo)主要集中在壓差4~8 hPa。

    單看一種氣象條件下O3的濃度特征,對(duì)預(yù)報(bào)的支撐作用有限。需要對(duì)主要?dú)庀笥绊懸蜃幼霾煌慕M合,找到最大程度(概率)影響O3濃度的最優(yōu)組合。根據(jù)氣象因子預(yù)報(bào)數(shù)據(jù)準(zhǔn)確度高、易獲得的原則,以最高溫度作為首要約束條件,層層篩選其他氣象影響因子,直至得到O3濃度最大相似的組合。通過這種方法,可以得到若干組不同氣象條件組合下對(duì)應(yīng)的O3濃度,這里只列出O3超標(biāo)時(shí)主要?dú)庀笠蜃拥牟糠纸M合(見表6)。

    表6 最優(yōu)化組合示例

    注:“空”表示不對(duì)該氣象參數(shù)進(jìn)行判斷。

    2.3結(jié)果評(píng)估

    利用2016年5—8月的O3預(yù)報(bào)結(jié)果進(jìn)行驗(yàn)證分析,如圖3所示,本研究對(duì)2016年5—8月O3類別的預(yù)報(bào)準(zhǔn)確率為72.4%,高于現(xiàn)有OVL模型的57.7%以及逐步回歸方法的68.3%。對(duì)O3超標(biāo)的預(yù)報(bào)準(zhǔn)確率為46%,高于OVL模型和逐步回歸的38%、31%。

    從預(yù)報(bào)結(jié)果來看,基于O3濃度相似的主要?dú)庀笞顑?yōu)組合的預(yù)報(bào)方法適合重慶本地O3預(yù)報(bào),能得到較高的預(yù)報(bào)準(zhǔn)確率。雖然準(zhǔn)確率有所提高,但該方法也存在各種模型的共同問題,當(dāng)空氣質(zhì)量為O3引發(fā)的輕度污染時(shí),前一日可以較準(zhǔn)確地判斷出首要污染物為O3,但對(duì)空氣質(zhì)量類別容易低估。

    圖3 本研究O3預(yù)報(bào)準(zhǔn)確率與其他方法對(duì)比Fig.3 The accuracy of O3 forecast compared with other methods

    3 結(jié)論

    1)通過主成分分析、逐步回歸分析等方法,確定影響重慶主城區(qū)O3濃度的主要?dú)庀笠蛩貫樽罡邷囟取夭?、太?yáng)輻射、降水量、相對(duì)濕度、水氣壓、壓差。

    2)通過基于O3污染水平相似的主要?dú)庀罂刂埔蜃雍Y選和最優(yōu)組合的預(yù)報(bào)結(jié)果優(yōu)化方法,對(duì)2016年5—8月O3類別的預(yù)報(bào)準(zhǔn)確率為72.4%,高于現(xiàn)有OVL模型的57.7%。對(duì)O3超標(biāo)的預(yù)報(bào)準(zhǔn)確率為46%,高于OVL模型的38%。該方法適合重慶本地O3預(yù)報(bào),能得到較高的預(yù)報(bào)準(zhǔn)確率。

    3)本研究?jī)H考慮了氣象條件的影響,下一步仍需考慮其前體物NO2、VOCs的影響,不斷完善規(guī)律,以提高預(yù)報(bào)準(zhǔn)確性。

    4)樣本量也是重要的影響因素,目前僅利用2013—2015年數(shù)據(jù),畢竟3年的天氣形勢(shì)會(huì)有所不同,需要多年甚至十幾年的監(jiān)測(cè)和氣象數(shù)據(jù),才能更好地反映本地特征。

    5)不同地區(qū)的氣象條件不同,對(duì)O3的影響規(guī)律也不盡相同,需要在分析規(guī)律的基礎(chǔ)上,對(duì)當(dāng)?shù)豋3預(yù)報(bào)進(jìn)行指導(dǎo)。

    [1] 王曉彥,陳佳,朱莉莉,等.城市環(huán)境空氣質(zhì)量指數(shù)范圍預(yù)報(bào)方法初探[J].中國(guó)環(huán)境監(jiān)測(cè),2015,31(6):139-142.

    WANG Xiaoyan,CHEN Jia,ZHU Lili, et al. Discussion on the methods of urban ambient air quality index range forecasting [J]Environmental Monitoring in China, 2015,31(6):139-142.

    [2] 張遠(yuǎn)航,邵可聲,唐孝炎,等.中國(guó)城市光化學(xué)煙霧污染研究[J].北京大學(xué)學(xué)報(bào)(自然科學(xué)版),1998,34(2):392-400.

    ZHANG Yuanhang, SHAO Kesheng, TANG Xiaoyan,et al. The study of urban photochemical smog pollution in China[J]. Acta Scientiarum Naturalium Universitatis Pekinensis (Natural Science Edition), 1998,34(2):392-400.

    [3] 王淑瑩,許榮,尹翠芳,等.基于OPAQ 的城市空氣質(zhì)量預(yù)報(bào)系統(tǒng)研究[J].中國(guó)環(huán)境監(jiān)測(cè),2016,32(3):13-20.

    WANG Shuying, XU Rong,YIN Cuifang, et al. Study on urban air quality forecasting with OPAQ[J]. Environmental Monitoring in China, 2016,32(3):13-20.

    [4] 中國(guó)環(huán)境監(jiān)測(cè)總站.環(huán)境空氣質(zhì)量預(yù)報(bào)預(yù)警方法技術(shù)指南[M].北京:中國(guó)環(huán)境科學(xué)出版社:2014:24-25.

    [5] 汪巍,高愈霄,趙熠琳,等.中小城市空氣質(zhì)量預(yù)報(bào)試驗(yàn)與問題研究[J].中國(guó)環(huán)境監(jiān)測(cè),2015,31(6):134-138.

    WANG Wei,GAO Yuxiao,ZHAO Yilin,et al. Analysis on the experiment and problems of air quality forecast in small and medium-sized cities in China[J].Environmental Monitoring in China,2015,31(6):134-138.

    [6] 安俊琳,王躍思,朱彬.主成分和回歸分析方法在大氣臭氧預(yù)報(bào)的應(yīng)用——以北京夏季為例[J].環(huán)境科學(xué)學(xué)報(bào),2010,30 (6): 1 286-1 294.

    AN Junlin, WANG Yuesi, Zhu Bin.Principal component and multiple regression analysis predicting ozone concentrations: Case study in summer in Beijing[J]. Acta Scientiae Circumstantiae, 2010,30 (6): 1 286-1 294.

    [7] 程念亮,李云婷,張大偉,等.2014年北京市城區(qū)臭氧超標(biāo)日濃度特征及與氣象條件的關(guān)系[J].環(huán)境科學(xué),2016,37(6):2 014-2 051.

    CHENG Nianliang, LI Yunting, ZHANG Dawei, et al. Characteristics of ozone over standard and its relationships with meteorological conditions in Beijing City in 2014[J]. Environmental Science, 2016,37(6):2 014-2 051.

    [8] 任麗紅,胡非,王瑋.北京夏季O3垂直分布與氣象因子的相關(guān)研究[J].氣候與環(huán)境研究,2005,10(2):166-174.

    REN Lihong, HU Fei, WANG Wei. A study of correlations between O3concentration and meteorological factors in summer in Beijing[J]. Climatic and Environmental Research, 2005,10(2):166-174.

    [9] 劉明花.上海市地面臭氧濃度分析及多元非線性預(yù)報(bào)模式研究[D].上海:華東師范大學(xué),2009.

    [10] 宋榕榮,王堅(jiān),張?jiān)鴱|,等.廈門市空氣質(zhì)量臭氧預(yù)報(bào)和評(píng)估系統(tǒng)[J].中國(guó)環(huán)境監(jiān)測(cè),2012,28(1):27-32.

    SONG Rongrong, WANG Jian, ZHANG Zengtao, et al. The air quality forecast and evaluation system of ozone in Xiamen[J]. Environmental Monitoring in China, 2012,28(1):27-32.

    [11] 王闖,王帥,楊碧波,等.氣象條件對(duì)沈陽(yáng)市環(huán)境空氣臭氧濃度影響研究[J].中國(guó)環(huán)境監(jiān)測(cè),2015,31(3):32-37.

    WANG Chuang, WANG Shuai, YANG Bibo, et al. Study of the effect of meteorological conditions on the ambient air ozone concentrations in Shenyang[J]. Environmental Monitoring in China, 2015,31(3):32-37.

    [12] 姚青,孫玫玲,劉愛霞.天津臭氧濃度與氣象因素的相關(guān)性及其預(yù)測(cè)方法[J].生態(tài)環(huán)境學(xué)報(bào),2009,18(6):2 206-2 210.

    YAO Qing, SUN Meiling, LIU Aixia. Analysis and prediction of surface ozone concentration and related meteorological factors in summer in Tianjin[J]. Ecology and Environmental Sciences, 2009, 18(6): 2 206-2 210.

    [13] 溫彥平.太原市近地面臭氧濃度變化規(guī)律及其與氣象要素的關(guān)系[J].環(huán)境工程學(xué)報(bào),2015,9(11):5 545-5 554.

    WEN Yanping. Temporal variations of surface ozone and its relations with meteorological parameters in Taiyuan City[J].Chinese Journal of Environmental Engineering, 2015,9(11):5 545-5 554.

    [14] 張浩月,王雪松,陸克定,等.珠江三角洲秋季典型氣象條件對(duì)O3和PM10污染的影響[J].北京大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,50(3):564-576.

    ZHANG Haoyue, WANG Xuesong, LU Keding, et al. Impact of typical meteorological conditions on the O3and PM10pollution episodes in the Pearl River Delta in Autumn[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2014,50(3):564-576.

    [15] 趙玉成,楊英蓮,蔡永祥,等.西寧地區(qū)近地面臭氧(O3)濃度的變化特征及氣象條件分析[J].青海環(huán)境,2013,23(2):76-89.

    ZHAO Yucheng, YANG Yinglian, CAI Yongxiang, et al. Characteristics of near-surface ozone (O3) concentration and meteorological conditions in Xining[J]. Qinghai Environment, 2013,23(2):76-89.

    [16] 范洋,樊曙先,張紅亮,等.臨安冬夏季SO2、NO2和O3體積分?jǐn)?shù)特征及與氣象條件的關(guān)系[J].大氣科學(xué)學(xué)報(bào),2013,36(1):121-128.

    FAN Yang,FAN Shuxian, ZHANG Hongliang,et al. Characteristics of SO2,NO2,O3volume fractions and their relationship with weather conditions at Linan in summer and winter[J].Transactions of Atmospheric Sciences,36 (1):121-128.

    [17] 劉姣姣,蔣昌潭,宋丹,等.重慶夏季近地面臭氧變化規(guī)律及影響因素分析[J].重慶大學(xué)學(xué)報(bào):自然科學(xué)版,2014,37(8):91-98.

    LIU Jiaojiao, JIANG Changtan, SONG Dan, et al. Analysis of Distribution Characteristics of Surface Ozone and its Influencing Factors in Summer in Chongqing[J]. Journal of Chongqing University: Natural Science Edition, 2014,37(8):91-98.

    [18] GUO H, WANG T, LOUIE PKK. Source apportionment of ambient non-methane hydrocarbons in Hong Kong: application of a principal component analysis (PCA)/absolute principal component scores (APCS) receptor model [J]. Environmental Pollution, 2004,129(3):489-498.

    [19] SINDOSI O A, KATSOU lIS B D, BARTZOKAS A. An objective definition of air mass types affecting Athens, Greece: the corresponding atmospheric pressure patterns and air pollution levels[J].Environmental Technology, 2003,24(8):947-962.

    [20] 劉閩,王帥,林宏,等.沈陽(yáng)市冬季環(huán)境空氣質(zhì)量統(tǒng)計(jì)預(yù)報(bào)模型建立及應(yīng)用[J].中國(guó)環(huán)境監(jiān)測(cè),2014,30(4):10-15.

    LIU Min, WANG Shuai, LIN Hong, et al. The study on establishment and application of winter environment air quality forecasting model of Shenyang [J]. Environmental Monitoring in China,2014,30(4):10-15.

    [21] 潘本鋒,程麟鈞,王建國(guó),等. 京津冀地區(qū)臭氧污染特征與來源分析[J].中國(guó)環(huán)境監(jiān)測(cè),2016,32(5):17-23.

    PAN Benfeng,CHENG Linjun,WANG Jianguo, et al. Characteristics and source attribution of ozone pollution in Beijing-Tianjin-Hebei Region[J]. Environmental Monitoring in China, 2016,32(5):17-23.

    EffectsofMeteorologicalConditionsonOzoneandItsApplicationinAirQualityForecast

    ZHANG Can1,2, JIANG Changtan1,2, LUO Caihong1,2, LIU Jiaojiao1,2, YE Di1,2, QIAO Jie3, HAN Shigang4

    1.Ecological and Environment Monitoring Centre of Chongqing, Chongqing 401147, China 2.Chongqing Key Laboratory of Urban Atmospheric Environment for Integrated Observation and Pollution Prevention and Control, Chongqing 401147, China 3.Chongqing Municipal Environmental Supervise Corps, Chongqing 401147, China 4.Chongqing Meteorological Service Centre, Chongqing 401147, China

    In order to improve the accuracy of O3forecast in Chongqing, O3monitoring concentration and meteorological data from May to October in 2013-2015 were used to analyze the effects of meteorological factors on O3concentration by principal component analysis and stepwise regression analysis. The result showed that O3concentration was influenced by maximum temperature, temperature difference, radiation, precipitation, humidity, vapor pressure and pressure difference. By selection of major meteorological control factors on O3and the optimal combination of control factors based on O3concentration similarity, the accuracy of O3forecast was improved from 57.7% to 72.4% and the accuracy of O3over standard forecast was improved from 38% to 46% during May to August in 2016.

    meteorological factors; ozone; forecast

    X823;X84;X51

    :A

    :1002-6002(2017)04- 0221- 08

    10.19316/j.issn.1002-6002.2017.04.28

    2016-11-11;

    :2017-01-12

    重慶市環(huán)境保護(hù)局環(huán)保科技項(xiàng)目[2015(7號(hào))]

    張 燦(1982-),女,山東金鄉(xiāng)人,碩士,高級(jí)工程師。

    羅財(cái)紅

    猜你喜歡
    壓差溫差氣壓
    手表與溫差
    看不見的氣壓
    幼兒畫刊(2021年5期)2021-12-02 04:24:04
    北方冬季養(yǎng)豬當(dāng)心“四差”
    壓力容器氣壓端蓋注射模設(shè)計(jì)
    模具制造(2019年4期)2019-06-24 03:36:46
    溫差“催甜”等
    在線處理脫丁烷塔壓差高
    低級(jí)熱溫差發(fā)電器的開發(fā)應(yīng)用研究
    電子制作(2018年23期)2018-12-26 01:01:20
    電滲—堆載聯(lián)合氣壓劈烈的室內(nèi)模型試驗(yàn)
    二維氣體動(dòng)力學(xué)中壓差方程的特征分解和簡(jiǎn)單波
    基于改進(jìn)經(jīng)濟(jì)壓差算法的主站AVC的研究
    涩涩av久久男人的天堂| 国产精品国产av在线观看| 亚洲人成77777在线视频| 亚洲熟女精品中文字幕| 黄色怎么调成土黄色| 美女国产高潮福利片在线看| 日韩电影二区| 99精国产麻豆久久婷婷| 七月丁香在线播放| 好男人视频免费观看在线| 女人被躁到高潮嗷嗷叫费观| 国产成人一区二区在线| 亚洲一区二区三区欧美精品| 中国国产av一级| 97人妻天天添夜夜摸| 免费少妇av软件| 大陆偷拍与自拍| 视频区图区小说| 久久ye,这里只有精品| 人妻少妇偷人精品九色| 国产精品国产av在线观看| 男的添女的下面高潮视频| 国产精品人妻久久久影院| 丰满少妇做爰视频| 欧美 日韩 精品 国产| 久久av网站| 亚洲一区二区三区欧美精品| 天堂8中文在线网| 麻豆精品久久久久久蜜桃| 少妇的逼水好多| 边亲边吃奶的免费视频| www.熟女人妻精品国产 | 汤姆久久久久久久影院中文字幕| 亚洲欧洲日产国产| 国产成人aa在线观看| 日韩视频在线欧美| 亚洲精品自拍成人| 久久人人爽av亚洲精品天堂| 欧美日韩亚洲高清精品| 久久女婷五月综合色啪小说| 一本色道久久久久久精品综合| 9191精品国产免费久久| tube8黄色片| 美女脱内裤让男人舔精品视频| 久久这里只有精品19| 一二三四中文在线观看免费高清| 国产爽快片一区二区三区| 久久鲁丝午夜福利片| 亚洲美女黄色视频免费看| 丰满迷人的少妇在线观看| 国产日韩欧美视频二区| 欧美激情国产日韩精品一区| 亚洲欧美成人综合另类久久久| 免费观看无遮挡的男女| 天天躁夜夜躁狠狠躁躁| 亚洲成人av在线免费| 久久这里只有精品19| 卡戴珊不雅视频在线播放| 久久久久久久大尺度免费视频| 在线天堂中文资源库| 日本色播在线视频| 久久99热这里只频精品6学生| 免费看av在线观看网站| 国产白丝娇喘喷水9色精品| 国产一区二区在线观看日韩| 18禁在线无遮挡免费观看视频| 18禁观看日本| 大陆偷拍与自拍| 十分钟在线观看高清视频www| 午夜福利影视在线免费观看| 热re99久久国产66热| 免费黄网站久久成人精品| 少妇高潮的动态图| 日韩大片免费观看网站| 日本欧美国产在线视频| 在线观看免费视频网站a站| 一级片免费观看大全| 国产激情久久老熟女| 人妻少妇偷人精品九色| 成年动漫av网址| 亚洲欧美精品自产自拍| 丰满迷人的少妇在线观看| 国产老妇伦熟女老妇高清| 亚洲三级黄色毛片| 中文欧美无线码| 久久人人爽av亚洲精品天堂| 亚洲av在线观看美女高潮| 人人妻人人澡人人爽人人夜夜| 18在线观看网站| 丰满迷人的少妇在线观看| 久久精品久久久久久久性| 人人妻人人爽人人添夜夜欢视频| 国产欧美日韩综合在线一区二区| 校园人妻丝袜中文字幕| 香蕉精品网在线| 国产xxxxx性猛交| 日韩熟女老妇一区二区性免费视频| 精品人妻偷拍中文字幕| 一本—道久久a久久精品蜜桃钙片| 亚洲综合精品二区| 欧美精品高潮呻吟av久久| 纵有疾风起免费观看全集完整版| 中文字幕最新亚洲高清| 在线观看www视频免费| 久久久久久久久久成人| 国产一区二区在线观看av| 久久精品国产综合久久久 | 国产成人精品婷婷| 麻豆乱淫一区二区| www.av在线官网国产| 日韩精品有码人妻一区| 亚洲av免费高清在线观看| 成年女人在线观看亚洲视频| 成人手机av| 91在线精品国自产拍蜜月| 久久精品国产a三级三级三级| 性高湖久久久久久久久免费观看| 日本av免费视频播放| 日韩不卡一区二区三区视频在线| 人人澡人人妻人| 丝袜脚勾引网站| 午夜福利视频在线观看免费| 久久精品aⅴ一区二区三区四区 | 久久人人爽人人爽人人片va| 大香蕉97超碰在线| 亚洲少妇的诱惑av| 国产精品一区二区在线观看99| 久久久久久久久久久久大奶| 中国三级夫妇交换| av卡一久久| 丰满乱子伦码专区| 男男h啪啪无遮挡| 热re99久久精品国产66热6| av视频免费观看在线观看| 日韩人妻精品一区2区三区| 在线观看三级黄色| 丝袜脚勾引网站| 久热这里只有精品99| 亚洲av欧美aⅴ国产| 18禁裸乳无遮挡动漫免费视频| 久久韩国三级中文字幕| 国产精品久久久久久精品古装| 欧美精品高潮呻吟av久久| 亚洲,欧美精品.| 亚洲成av片中文字幕在线观看 | 亚洲精华国产精华液的使用体验| 国产在线视频一区二区| 午夜福利网站1000一区二区三区| 波多野结衣一区麻豆| 最近最新中文字幕免费大全7| 校园人妻丝袜中文字幕| 在现免费观看毛片| 熟妇人妻不卡中文字幕| 宅男免费午夜| 国产成人精品久久久久久| 五月伊人婷婷丁香| 天美传媒精品一区二区| a级毛片黄视频| 十八禁高潮呻吟视频| 免费在线观看完整版高清| 欧美日韩一区二区视频在线观看视频在线| 国产精品一国产av| 久久久久人妻精品一区果冻| 丰满饥渴人妻一区二区三| 十八禁网站网址无遮挡| 精品国产乱码久久久久久小说| 国产成人免费无遮挡视频| 十八禁高潮呻吟视频| 久久国产精品大桥未久av| 亚洲欧洲国产日韩| 中文字幕最新亚洲高清| 美女国产高潮福利片在线看| 精品人妻一区二区三区麻豆| 男女边摸边吃奶| 久久影院123| 欧美3d第一页| 欧美xxⅹ黑人| 日日摸夜夜添夜夜爱| 老熟女久久久| 一级片免费观看大全| xxxhd国产人妻xxx| 亚洲精品久久成人aⅴ小说| 亚洲精品自拍成人| 咕卡用的链子| 亚洲精华国产精华液的使用体验| 少妇人妻久久综合中文| 麻豆乱淫一区二区| 中文字幕另类日韩欧美亚洲嫩草| 国产片内射在线| 如何舔出高潮| 高清黄色对白视频在线免费看| 国产视频首页在线观看| 亚洲少妇的诱惑av| 日本av免费视频播放| 咕卡用的链子| 国产精品不卡视频一区二区| 爱豆传媒免费全集在线观看| 亚洲精品中文字幕在线视频| 纵有疾风起免费观看全集完整版| 久久 成人 亚洲| 91成人精品电影| 视频在线观看一区二区三区| 久久99热这里只频精品6学生| 免费看av在线观看网站| 天天操日日干夜夜撸| 哪个播放器可以免费观看大片| 国产精品一区二区在线不卡| 国产在线一区二区三区精| 国产探花极品一区二区| 十分钟在线观看高清视频www| 在现免费观看毛片| 久久久久久久久久成人| 黄片无遮挡物在线观看| 国产亚洲精品久久久com| h视频一区二区三区| 只有这里有精品99| 中文字幕精品免费在线观看视频 | 一二三四中文在线观看免费高清| 国产又色又爽无遮挡免| av又黄又爽大尺度在线免费看| 九色亚洲精品在线播放| 欧美日本中文国产一区发布| 久久久久视频综合| 国产精品久久久久久精品古装| 免费高清在线观看日韩| 又粗又硬又长又爽又黄的视频| 秋霞在线观看毛片| 日本免费在线观看一区| 久久 成人 亚洲| 久久国产精品男人的天堂亚洲 | 精品亚洲成国产av| 国产女主播在线喷水免费视频网站| 欧美日韩视频高清一区二区三区二| 亚洲伊人色综图| 这个男人来自地球电影免费观看 | 一级片'在线观看视频| 建设人人有责人人尽责人人享有的| 色94色欧美一区二区| 一级毛片电影观看| 精品亚洲成a人片在线观看| 最新的欧美精品一区二区| 91精品国产国语对白视频| 美女中出高潮动态图| freevideosex欧美| 观看美女的网站| 男女午夜视频在线观看 | 九九爱精品视频在线观看| 日本黄大片高清| 午夜免费鲁丝| 午夜免费男女啪啪视频观看| 视频区图区小说| 免费黄色在线免费观看| 久久精品国产亚洲av涩爱| 成人国产av品久久久| 久久ye,这里只有精品| 亚洲精品一区蜜桃| 91aial.com中文字幕在线观看| 欧美人与性动交α欧美软件 | 99久久综合免费| 日韩精品免费视频一区二区三区 | 亚洲欧美中文字幕日韩二区| 看免费成人av毛片| 黄网站色视频无遮挡免费观看| 大陆偷拍与自拍| 天美传媒精品一区二区| 亚洲人成网站在线观看播放| 日本wwww免费看| 国产日韩欧美视频二区| 伊人亚洲综合成人网| 热re99久久国产66热| 午夜免费男女啪啪视频观看| 免费不卡的大黄色大毛片视频在线观看| 婷婷色av中文字幕| 国产永久视频网站| 日韩免费高清中文字幕av| 99re6热这里在线精品视频| 高清av免费在线| 国产日韩欧美亚洲二区| 中文字幕人妻熟女乱码| 国产亚洲一区二区精品| 日韩一本色道免费dvd| 午夜91福利影院| 久久av网站| 国产国拍精品亚洲av在线观看| 国产深夜福利视频在线观看| 69精品国产乱码久久久| 精品卡一卡二卡四卡免费| 午夜精品国产一区二区电影| 久久久久久久久久成人| 香蕉国产在线看| 中国国产av一级| 国产一级毛片在线| 看免费成人av毛片| 国产麻豆69| 嫩草影院入口| 日韩视频在线欧美| 一区二区三区精品91| 亚洲精品久久午夜乱码| av在线老鸭窝| 男女啪啪激烈高潮av片| 日本vs欧美在线观看视频| 色婷婷久久久亚洲欧美| 久久久精品免费免费高清| 亚洲天堂av无毛| 美国免费a级毛片| 日本欧美视频一区| 国产福利在线免费观看视频| 少妇的逼好多水| 最近2019中文字幕mv第一页| 日韩精品免费视频一区二区三区 | 国产欧美日韩综合在线一区二区| 免费大片18禁| 飞空精品影院首页| 国产白丝娇喘喷水9色精品| 少妇人妻 视频| 毛片一级片免费看久久久久| 欧美成人精品欧美一级黄| www.av在线官网国产| 国产极品粉嫩免费观看在线| 校园人妻丝袜中文字幕| 久久久久久人妻| 一级毛片黄色毛片免费观看视频| 国产亚洲欧美精品永久| 人妻少妇偷人精品九色| 宅男免费午夜| 日韩不卡一区二区三区视频在线| 亚洲av免费高清在线观看| 日韩一区二区视频免费看| freevideosex欧美| 综合色丁香网| 啦啦啦啦在线视频资源| 亚洲在久久综合| 国产精品一国产av| 亚洲精品av麻豆狂野| 精品99又大又爽又粗少妇毛片| 亚洲国产日韩一区二区| 久久精品久久久久久噜噜老黄| 精品少妇内射三级| av在线app专区| 久久久久久久久久人人人人人人| 七月丁香在线播放| 亚洲国产最新在线播放| 久久久国产欧美日韩av| 寂寞人妻少妇视频99o| 七月丁香在线播放| 在线观看国产h片| 性色av一级| 夜夜爽夜夜爽视频| 插逼视频在线观看| 免费久久久久久久精品成人欧美视频 | 亚洲美女视频黄频| 黄色 视频免费看| 国产黄频视频在线观看| av电影中文网址| 一区二区日韩欧美中文字幕 | 91精品伊人久久大香线蕉| 欧美97在线视频| 最新中文字幕久久久久| 狂野欧美激情性bbbbbb| 90打野战视频偷拍视频| 免费av中文字幕在线| 天堂俺去俺来也www色官网| 精品一品国产午夜福利视频| 国产精品麻豆人妻色哟哟久久| 日韩av不卡免费在线播放| 亚洲精品av麻豆狂野| 青春草视频在线免费观看| 国产一区二区在线观看日韩| 五月天丁香电影| 精品人妻熟女毛片av久久网站| 久久这里只有精品19| 国产日韩欧美视频二区| 精品国产一区二区三区久久久樱花| 日本av手机在线免费观看| 免费人成在线观看视频色| 国产麻豆69| 女的被弄到高潮叫床怎么办| 男女啪啪激烈高潮av片| videossex国产| 韩国av在线不卡| 国产亚洲精品第一综合不卡 | 免费大片18禁| 久久久国产欧美日韩av| 在线观看免费日韩欧美大片| 王馨瑶露胸无遮挡在线观看| 欧美日韩亚洲高清精品| 中文字幕人妻丝袜制服| 亚洲av电影在线观看一区二区三区| 熟妇人妻不卡中文字幕| 久久久亚洲精品成人影院| 国产精品久久久久久av不卡| 久久婷婷青草| 久久精品久久久久久噜噜老黄| 下体分泌物呈黄色| 久久久久国产网址| 高清av免费在线| 高清毛片免费看| 91午夜精品亚洲一区二区三区| videossex国产| 深夜精品福利| 欧美精品国产亚洲| 中文字幕免费在线视频6| a级毛片在线看网站| 欧美日韩国产mv在线观看视频| 久久av网站| 欧美精品国产亚洲| 在线观看免费视频网站a站| 欧美另类一区| 男人添女人高潮全过程视频| 男女边吃奶边做爰视频| 国产精品人妻久久久久久| 亚洲情色 制服丝袜| 80岁老熟妇乱子伦牲交| 成人18禁高潮啪啪吃奶动态图| 91精品三级在线观看| 日本爱情动作片www.在线观看| 久久鲁丝午夜福利片| 国产精品久久久久久av不卡| 全区人妻精品视频| 女性被躁到高潮视频| 成人毛片a级毛片在线播放| 国产精品99久久99久久久不卡 | 国产成人精品在线电影| 春色校园在线视频观看| 伊人亚洲综合成人网| 国产日韩欧美视频二区| 国产熟女午夜一区二区三区| 免费黄网站久久成人精品| 亚洲伊人久久精品综合| 国产又爽黄色视频| 欧美另类一区| 欧美日本中文国产一区发布| 精品卡一卡二卡四卡免费| 人体艺术视频欧美日本| 色婷婷久久久亚洲欧美| 亚洲av国产av综合av卡| 999精品在线视频| 99热国产这里只有精品6| 免费黄网站久久成人精品| 美女国产视频在线观看| 久久热在线av| 国产精品不卡视频一区二区| 国精品久久久久久国模美| 久久久久久久久久人人人人人人| 老司机亚洲免费影院| 午夜福利网站1000一区二区三区| 韩国高清视频一区二区三区| 99热国产这里只有精品6| 伊人久久国产一区二区| 一级a做视频免费观看| 午夜视频国产福利| 免费大片18禁| 最近中文字幕2019免费版| 精品久久久精品久久久| 欧美精品人与动牲交sv欧美| 超碰97精品在线观看| 免费在线观看黄色视频的| av黄色大香蕉| 高清av免费在线| 97超碰精品成人国产| 天堂8中文在线网| 亚洲av电影在线进入| 久久国内精品自在自线图片| 国产免费一级a男人的天堂| 在线免费观看不下载黄p国产| 国产精品无大码| 中文字幕免费在线视频6| 成人国产av品久久久| 日韩一区二区三区影片| 男人操女人黄网站| 久久久久久人妻| 一本—道久久a久久精品蜜桃钙片| 久热这里只有精品99| 亚洲情色 制服丝袜| 亚洲国产日韩一区二区| 国产免费一区二区三区四区乱码| 狠狠婷婷综合久久久久久88av| 成年人午夜在线观看视频| 中文天堂在线官网| 亚洲av免费高清在线观看| 亚洲av中文av极速乱| 女的被弄到高潮叫床怎么办| 少妇人妻 视频| 丝袜喷水一区| 黑人猛操日本美女一级片| 在线精品无人区一区二区三| 啦啦啦视频在线资源免费观看| 久久女婷五月综合色啪小说| 亚洲av免费高清在线观看| 日韩人妻精品一区2区三区| 啦啦啦啦在线视频资源| www.色视频.com| 精品人妻在线不人妻| 久久久精品94久久精品| 不卡视频在线观看欧美| 久久午夜福利片| 边亲边吃奶的免费视频| 国产精品久久久久久精品古装| 女人精品久久久久毛片| 午夜免费观看性视频| 国产视频首页在线观看| 午夜视频国产福利| 国产亚洲欧美精品永久| 涩涩av久久男人的天堂| 亚洲熟女精品中文字幕| 啦啦啦啦在线视频资源| 亚洲欧洲国产日韩| 涩涩av久久男人的天堂| 国产视频首页在线观看| 久久精品人人爽人人爽视色| 亚洲欧美一区二区三区国产| 69精品国产乱码久久久| 激情视频va一区二区三区| 欧美精品一区二区免费开放| 亚洲伊人色综图| 一级片'在线观看视频| 国产午夜精品一二区理论片| 亚洲精品久久午夜乱码| 中文字幕另类日韩欧美亚洲嫩草| 丰满饥渴人妻一区二区三| 免费播放大片免费观看视频在线观看| 天天影视国产精品| 亚洲国产精品一区三区| 国产成人精品福利久久| 国产精品熟女久久久久浪| 日韩av在线免费看完整版不卡| 亚洲国产欧美日韩在线播放| 各种免费的搞黄视频| a级片在线免费高清观看视频| 日本爱情动作片www.在线观看| 成人亚洲精品一区在线观看| 老司机亚洲免费影院| 一级毛片黄色毛片免费观看视频| 日韩成人av中文字幕在线观看| 日韩欧美一区视频在线观看| 美女主播在线视频| 国产成人一区二区在线| 亚洲欧洲精品一区二区精品久久久 | 少妇熟女欧美另类| 男女边摸边吃奶| 在线免费观看不下载黄p国产| 日韩人妻精品一区2区三区| 最近2019中文字幕mv第一页| 边亲边吃奶的免费视频| 免费黄色在线免费观看| 美女中出高潮动态图| 国产国拍精品亚洲av在线观看| 考比视频在线观看| 狂野欧美激情性xxxx在线观看| 久久av网站| 久久久久久久精品精品| www.av在线官网国产| 啦啦啦视频在线资源免费观看| 国产白丝娇喘喷水9色精品| av国产久精品久网站免费入址| 中文字幕另类日韩欧美亚洲嫩草| 51国产日韩欧美| 2022亚洲国产成人精品| 亚洲精品中文字幕在线视频| 国产 一区精品| 中文字幕免费在线视频6| 视频区图区小说| 国产乱人偷精品视频| 十分钟在线观看高清视频www| 亚洲第一av免费看| 少妇的逼好多水| 另类亚洲欧美激情| 久久人人爽人人爽人人片va| 久久婷婷青草| av一本久久久久| 午夜福利,免费看| 国产1区2区3区精品| 99视频精品全部免费 在线| 亚洲av免费高清在线观看| 一级,二级,三级黄色视频| 国产日韩欧美视频二区| 一区在线观看完整版| 高清在线视频一区二区三区| 亚洲精品色激情综合| 国产xxxxx性猛交| 在线观看三级黄色| 国产淫语在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产高清三级在线| 亚洲精品,欧美精品| 久久久久久久大尺度免费视频| 美女国产视频在线观看| 亚洲少妇的诱惑av| 各种免费的搞黄视频| 97超碰精品成人国产| 亚洲色图综合在线观看| 9热在线视频观看99| 成人黄色视频免费在线看| 国产亚洲精品第一综合不卡 | 久久国产精品男人的天堂亚洲 | 美女中出高潮动态图| 亚洲av电影在线观看一区二区三区| 在线天堂中文资源库| 欧美日韩视频高清一区二区三区二| 欧美最新免费一区二区三区| 亚洲精品日本国产第一区| 欧美精品一区二区免费开放| 韩国精品一区二区三区 | 青春草视频在线免费观看| 欧美性感艳星| 伦理电影大哥的女人| 国产亚洲精品久久久com| 午夜日本视频在线| 久久久久久人人人人人| 18禁裸乳无遮挡动漫免费视频| 久久久精品免费免费高清| 黑人猛操日本美女一级片| 女性生殖器流出的白浆| a级毛片黄视频| 最新的欧美精品一区二区|