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

    桂林市霾天氣的潛勢預報模型研究

    2020-07-09 08:03:46龍鳳翔張瑀琳李艷玉龔冬英
    四川環(huán)境 2020年3期
    關鍵詞:潛勢漏報天氣

    龍鳳翔, 張瑀琳, 李艷玉,龔冬英

    (桂林市氣象局,廣西 桂林 541001)

    霾是大量極細微的干塵粒等均勻地浮游在空中,使水平能見度小于10km的空氣普遍混濁現(xiàn)象,使遠處光亮物微帶黃、紅色,使黑暗物微帶藍色[1]。霾屬于天氣現(xiàn)象中的一種,隨著經濟快速發(fā)展、人類污染排放加劇,霾帶來的危害逐漸加重,除影響人們身心健康外,還對交通、旅游、工業(yè)等行業(yè)產生直接影響,加強霾的監(jiān)測,分析霾的特征,發(fā)布霾天氣預報,可以為防御霾、減輕霾的危害發(fā)揮積極作用。

    潛勢預報作為霾預報常見方法,它是選取與預報對象有關的物理因子,運用概率統(tǒng)計的原理,在預報對象與預報因子之間建立某種統(tǒng)計關系,根據(jù)這種統(tǒng)計關系結果預報未來出現(xiàn)的可能性,它的預報對象以二分變量較為常見。目前常見預報方法有逐步回歸法、指標疊加法、支持向量(SVM)方法、決策樹法、Bayes分類法等。如毛程燕[2]等用逐步回歸法建立浙江中西部大霧天氣的潛勢預報模型;梁維亮[3]統(tǒng)計了南寧地區(qū)雷暴天氣和常用對流參數(shù)關系,將對流參數(shù)二值化后,采用邏輯回歸(Logistic)法建立了雷暴天氣預報概率的預報模型;陳豫英[4]用概率回歸方法建立了寧夏春季2-5月沙塵和強沙塵兩個等級的72h沙塵概率預報方程。張慶奎[5]運用指標疊加法建立了阜陽市大霧天氣潛勢預報模型,并取得了較好地預報效果。毛宇清[6]等用SVM分類和回歸方法建立了南京地區(qū)霾日分類預報和有霾日14時能見度預報方程。謝永華[7]等針對霾預報模型少,準確率低缺點,提出將統(tǒng)計預報和數(shù)值預報相結合,用多元逐步與概率回歸法混合的霾預報方法,這種混合的回歸預報算法準確率得到了顯著提高。以上各種潛勢預報方法,都在實際應用中取得了良好的預報效果。

    桂林地處廣西東北部,是國際著名旅游勝地,近年來隨著經濟迅速發(fā)展和城市化推進,霾污染帶來的大氣環(huán)境問題對旅游業(yè)的發(fā)展產生了不利影響,為了降低霾的不利影響,發(fā)布較為準確的霾天氣預報,本文采用基于Logistic和事件概率回歸法、指標疊加法等三種統(tǒng)計方法建立桂林市霾的潛勢預報模型,研究各模型的優(yōu)劣,為桂林霾預警、預報和防治提供重要參考依據(jù)。

    1 資料來源

    預報建模使用的資料包括2012年1月~2015年12月桂林自動氣象站逐時地面常規(guī)氣象觀測資料、桂林高空氣象觀測站逐日08時探空資料、桂林大氣成分站逐時PM2.5小時觀測數(shù)據(jù)。利用上述觀測資料,建立20時起報的桂林市霾未來24小時潛勢預報模型。

    模型檢驗使用的資料有:桂林自動氣象站2016年逐時地面常規(guī)觀測資料、PM2.5觀測資料,逐日08時歐洲中心(ECMWF)細網(wǎng)格(0.25°×0.25°)數(shù)值預報產品。

    2 結果與討論

    2.1 預報因子的選擇

    霾天氣是高濃度大氣顆粒物和特殊氣象條件共同作用的結果[8]。出現(xiàn)霾天氣時常干燥無雨,風速小或靜風,大氣層結穩(wěn)定,近地層伴有逆溫,大氣擴散能力減弱,污染物濃度增大[9-10]。因此霾的預報因子主要從污染物濃度、水汽、風速、大氣層結、降水等幾個方面選取,最終得到預報因子29個。

    2.1.1 污染物濃度因子

    PM2.5質量濃度水平與霾的生成密切相關,其質量濃度上升時會造成大氣能見度惡化,是霾天氣產生的重要因素之一[11~13],故選擇PM2.5質量濃度作為污染物濃度因子。包括前日的PM2.5日均質量濃度(PM2.5Day)、20時PM2.5小時質量濃度(PM2.5Hour)以及PM2.5Day、PM2.5Hour過去24小時變化值(ΔPM2.5Day、ΔPM2.5Hour),通過它們能較好地反映出霾出現(xiàn)前期整個大氣的污染狀況。其中ΔPM2.5Day、ΔPM2.5Hour反映的是過去24小時內污染物濃度變化的趨勢,當PM2.5濃度在過去24小時內有增加時,預示空氣污染加重,對霾的生成有利。

    2.1.2 水汽因子

    包括當日20時地面相對濕度(U)、氣溫露點差(T-Td),霾出現(xiàn)時當日08時高空925hPa氣溫露點差(T-Td925)、850hPa氣溫露點差(T-Td850)以及925hPa 露點溫度(Td925)、850hPa 露點溫度(Td850)、925hPa相對濕度(U925)、850hPa相對濕度(U850)、700hPa相對濕度(U700),當中低層濕度太大時容易出現(xiàn)降水,不利于霾的形成和持續(xù)。

    2.1.3 風速因子

    風速影響著大氣污染物擴散的快慢,尤其是中、低層風速大小直接影響著污染物的濃度,是霾形成的重要指標之一。在風速因子方面,選取當日08時高空925hPa風速(WS925)、850hPa風速(WS850)、700hPa風速(WS700)作為霾發(fā)生時近地層、中低層風狀況。

    2.1.4 層結條件

    霾出現(xiàn)天氣常為靜穩(wěn)天氣,近地層常伴有逆溫,當層結不穩(wěn)定時,大氣湍流加劇,污染物極易得到稀釋和擴散。選取因子時考慮用當日的K指數(shù)(K)、S指數(shù)(S)、地面-1 000hPa溫度差(T-T1 000)、地面-925hPa溫度差(T-T925)、地面-850hPa溫度差(T-T850)、925-850hPa溫度差(T925-850)、850-700hPa溫度差(T850-700)代表了地面、近地層、中層的大氣層結狀況。

    2.1.5 降水因子

    霾出現(xiàn)前期及期間天氣干燥,常持續(xù)無雨天氣或弱降水天氣,當前期、當日降水較大時,則不利于霾生成。霾出現(xiàn)當日降水(R)、前日降水(R24)、前3日累計降水(R72)3個因子能較好反映霾前期、霾期間的降水狀況。

    2.1.6 其他因子

    主要包括前日24小時變溫(ΔT24)、24小時最高氣溫變化(ΔTmax)、24小時最低氣溫變化(ΔTmin)。24小時變溫能反映出冷空氣強弱,當負變溫幅度越大時,代表冷空氣越強,而冷空氣過境則常代來大風、降水等天氣,使得污染物得到擴散。

    2.2 預報因子與霾天氣的相關性研究

    霾預報因子應選取與霾實況顯著相關的因子,對于一些與霾相關不大的因子,應予以剔除,具體因子篩選時,通過兩者相關系數(shù)大小進行取舍。

    在計算相關系數(shù)時,將霾出現(xiàn)與否轉為0,1二值變量,通過點二列相關求出兩者相關系數(shù)。計算得知,PM2.5Day、PM2.5Hour、U、T-Td、R72、WS850與霾出現(xiàn)與否相關性較好,相關系數(shù)(r)均≥0.20,相關系數(shù)最大為PM2.5Hour。在所有初選因子中,ΔTmax、T-T925、T-Td850、T850-700、SI未通過相關性檢查。表1為通過相關性檢查的預報因子及與霾的相關系數(shù),從相關系數(shù)來看,霾出現(xiàn)主要與前日PM2.5濃度、當日降水、前期降水及前期近地層水汽條件密切相關。這些預報因子從不同角度的反映了霾前期的污染程度、霾前期及霾出現(xiàn)時的水汽、層結、風速、降水等條件,均能不同程度的反映了霾天氣潛勢,選取它們建立霾潛勢預報方程,可綜合的判斷出霾天氣出現(xiàn)的可能性。

    2.3 預報因子閾值及二值化

    為獲取上述各預報因子閾值,分別將霾日對應的預報因子累計頻率進行計算,從頻率較大范圍(50%~90%)中的連續(xù)選取閾值。當預報因子值滿足閾值條件時記為1,否則為0,將轉換為二值變量的預報因子與霾出現(xiàn)與否的二值變量計算相關系數(shù),選擇相關系數(shù)最大對應的閾值作為最終確認閾值。各預報因子的最佳閾值、最佳閾值出現(xiàn)概率及二值化后與霾天氣的相關系數(shù)見表1。

    表1 預報因子最佳閾值、最佳閾值出現(xiàn)概率、二值化后及實際值與霾天氣的相關系數(shù)Tab.1 The optimum threshold of forecasting factor,the optimal threshold occurrence probability,correlation coefficient between actual value and haze weather after-binarization

    注:**相關系數(shù)在0.01置信水平下顯著(雙側檢驗); *相關系數(shù)在0.05置信水平下顯著(雙側檢驗)。

    上述結果顯示,預報因子經過二值轉換后與霾天氣均有顯著的正相關,且多數(shù)因子二值化后與霾天氣的相關系數(shù)要比實際相關系數(shù)有明顯增加,僅U700、WS850因子與霾天氣相關系數(shù)有所減小。綜合來看,用最佳閾值對預報參數(shù)進行二值轉換后的因子值比實際預報因子值使用更加合理,能更好地預測霾天氣是否出現(xiàn)。

    在所選預報因子中,有些因子對霾的消空有較好的指示意義,可以作為消空指標使用,如R、R72、T-Td、PM2.5Hour、T-T925。從表2可以看出,當R>5.0mm、R72>25mm、T-T925<1.0℃時霾出現(xiàn)的概率不足5%,T-Td <4℃、PM2.5Hour<20μg/m3時霾出現(xiàn)概率不足10%,因此可以將這些取值作為霾的消空閾值,即達到這些值時,當日霾不易出現(xiàn)。

    表2 霾日各因子消空指標取值及出現(xiàn)天數(shù)、累計百分比Tab.2 The values and occurrence days of blanking index of factors in haze days and cumulative percentage of haze days

    2.4 模型建立

    模型建立使用的資料為2012~2015年霾日前一日20時地面氣象觀測資料、PM2.5資料及當日08時探空資料。

    2.4.1 基于事件概率的多元線性回歸法

    事件概率回歸方法基本的思路是把預報對象和預報因子看成隨機事件并根據(jù)事件出現(xiàn)與否化成0,1的二值變量,把預報事件y與m個預報因子x1, x2, …,xm事件之間的關系看成是在m個因子已經出現(xiàn)的條件下事件出現(xiàn)與否的關系,是用逐步回歸的方法求各預報因子在各種組合狀態(tài)下預報對象出現(xiàn)的條件概率[14]。以X1,X2,...,Xm為因變量,P為自變量,建立線性回歸方程:

    (1)

    式中P為Y=1出現(xiàn)的概率,b0,b1,b2,...,bm為回歸系數(shù),用最小二乘法求出。

    根據(jù)以上方法得到霾概率多元線性回歸方程為:

    Y=0.176+0.174X1+0.125X2+0.102X3

    +0.121X4+0.089X5+0.091X6+0.078X7+0.078X8

    (2)

    (2)式最后納入方程的預報因子為PM2.5Hour(X1)、R(X2)、U(X3)、PM2.5Day(X4)、WS850(X5)、R72(X6)、WS925(X7)、ΔTmin(X8)。

    采用復相關系數(shù)對式(2)進行檢驗,在顯著性水平a=0.05時,臨界值R0.05=0.195,相關系數(shù)臨界值為F> F0.05,說明霾出現(xiàn)概率預報方法在統(tǒng)計意義上是顯著的。

    2.4.2 邏輯回歸法

    邏輯回歸是在線性回歸基礎上套用一個邏輯函數(shù)?;驹硎怯靡唤M觀測數(shù)據(jù)擬合Logistic模型,是研究二分類變量結果Y對一些影響因素(x1,x2,...,xm)之的依存的關系[14]。其中自變量xi稱為危險因素或者暴露因素,可以是分類變量、等級變量或者是連續(xù)變量。

    在m個自變量影響下,發(fā)生的概率記作:

    P=P(Y=1︱x1,x2,…,xm) 0≤P≤1

    (3)

    它與自變量之間的Logistic回歸模型為:

    (4)

    其中,β0為與自變量無關的常數(shù)項,β1,β2,...,βm為回歸系數(shù),表示自變量xi改變一個單位時,對Logistic(P)的貢獻。

    根據(jù)Logistic回歸法原理,利用SPSS 18.0軟件進行邏輯回歸分析,得霾出天氣潛勢預報方程:

    (5)

    其中,P為霾出現(xiàn)的潛勢預報因子a=2.447-0.813 X1+0.017 X2-0.669 X3-0.024 X4-1.027 X5-1.752 X6-0.399 X7-1.441 X8+0.017 X9-0.442 X10-0.213 X11-0.088 X12。其中,X1、X3、X5、X6、X7、X8、X9分別為PM2.5day、 PM2.5hour、T-Td、R、R72、U925、U850原始值轉換后的二值變量,X2、X4、X10、X11、X12分別為PM2.5hour、U、U850、T-Td925、WS850原始值。

    從SPSS 18.0計算結果來看,Cox Snell R方與Nagel Kerke R方值分別為0.311、0.501,兩值擬合效果較好。從另一擬合效果來看,如預測變量分界值設為0.35,當預測沒有霾出現(xiàn)時,預測無霾的準確率為87.4%;而當預測有霾時,預報有霾的準確率達73.6%,總的準確率為84.2%,說明方程總體較穩(wěn)定。

    2.4.3 指標疊加法

    指標疊加法是用來預報災害性天氣的一種方法,該方法已在強對流天氣預報中取得了不錯的預報效果[15-16]。具體思路是選取一定數(shù)量的與霾天氣相關性較好的預報因子,然后確定發(fā)生霾和未發(fā)生霾時的因子閾值,當預報因子超過閾值時,記1,若達不到閾值則記0,當有N個超過閾值時,記Y=N,只要分析Y大于等于霾出現(xiàn)時的臨界值N0,就可以預報霾出現(xiàn)。

    Y=X1+X2+X3+ …XN -1+XN

    (6)

    指標疊加法關鍵是在方程建立前,找到預報因子,方法如下:

    從表1中,選取預報因子標準為二值后與霾相關系數(shù)大于0.3的因子,符合條件因子共9個,建立的潛勢預報方程如下:

    Y=X1+X2+X3+X4+X5+X6+X7+X8+X9

    (7)

    其中X1~X9分別表示PM2.5day、PM2.5hour、U、T-Td、,T-Td925、K、R、R24、R72。當9個因子均達到閾值時,Y取值為9,說明是霾出現(xiàn)的最高閾值,而Y為0時為霾天氣潛勢預報的最小值。

    2.5 預報流程

    利用程序讀取前日地面氣象資料、前日PM2.5資料、當日08時ECMWF細網(wǎng)格數(shù)值預報產品,首先做消空判斷,當達到消空判斷條件時退出計算模型,如未達到則將各預報因子進行二值化處理,后代入潛勢預報模型進行計算,最后得出霾天氣潛勢預報產品,詳細預報流程見下圖。

    2.6 預報檢驗

    預報評分是檢驗預報質量的一個重要手段,它是按時段內發(fā)布的預報與天氣實況評定報對次數(shù)、漏報次數(shù)及空報次數(shù)。根據(jù)2005年中國氣象局頒布的《中短期天氣預報檢驗辦法》,對于災害性天氣落區(qū)預報,主要有TS評分、漏報率(PO)、空報率(FAR)3個指標,計算公式如下:

    圖 桂林霾天氣潛勢預報流程Fig. Haze Weather Potential Forecasting Process in Guilin

    (8)

    (9)

    (10)

    其中,NA為預報正確次數(shù),NB為空報次數(shù),NC為漏報次數(shù)。

    2.6.1 回報檢驗

    在經過對方程的顯著性檢驗后,我們對 2012~2015年的霾天氣實況進行回報,按霾潛勢預報流程,先對消空因子進行條件判斷,當達到消空條件時,則判斷該日無霾,否則代入潛勢預報方程。按此流程得到基于事件概率的多元回歸模型、邏輯回歸模型、指標疊加法三種潛勢預報方程不同概率閾值的霾預報評分,見表3、表4、表5。

    根據(jù)樣本回報的質量評分表可以看出,預報概率不同取值時,預報評分完全不同。其中,基于事件概率的多元回歸模型預測概率P=0.4時,空報率較高,達到47.3%,而在P=0.5時漏報率最高,高達59.2%,同時TS評分也較低。綜合來看,取P=0.45時TS評分最高,因此選取0.45作為預報概率閾值,即:當預報概率大于0.45時,預報未來24h有霾出現(xiàn),否則沒有。而邏輯回歸模型預測概率P從0.30~0.45分別取值的結果來看,差別并不像事件概率回歸法那樣明顯,各概率值對應TS分均能超過50%,空報率、漏報率均在30%左右,其中在概率閾值取0.35時,TS得分最高,大于或者小于該閾值時TS得分稍有下降,因此選取0.35作為事件概率回歸模型預報概率閾值。最后的指標疊加法預報概率閾值分別取4~6,其中閾值取5時TS評分最高,低于該值時空報率高,而高于該值時則是漏報率高, TS評分也低,在閾值為大于7時,漏報率高達73.3%,TS評分只有23.4%。

    表3 基于事件概率的多元回歸模型回報檢驗Tab.3 Multivariate Regression Model Return Test based on Event Probability

    表4 邏輯回歸模型回報檢驗Tab.4 Logistic Regression Model Return Test

    表5 指數(shù)疊加法回報檢驗Tab.5 Exponential superposition Return Test

    比較以上三種預測模型的回報檢驗效果來看,邏輯回歸模型回報檢驗驗效果最好,最穩(wěn)定,TS評分最高達52.5%,漏報率和缺測率對比其他兩種模型在TS評分最高時的漏報率、空報率都要低,而基于事件概率的多元回歸模型回報檢驗效果次之,TS評分最高為50.1%。三種模型中指數(shù)疊加模型回報檢驗效果稍差,TS評分最高只有47.8%,對應空報率、漏報率也較高,預測概率閾值在不同的取值時,TS評分變化幅度大。

    2.6.2 預報檢驗

    為驗證三種預報模式的實際預報效果,取2016年逐日20時地面觀測資料、PM2.5資料、次日08時 ECMWF細網(wǎng)格數(shù)值預報產品預報未來24小時是否出現(xiàn)霾天氣。預報過程同樣先做消空判斷,達到消空條件則預報無霾,否則代入潛勢預報模型。預報試驗時三種模型預報取值如下:事件概率法模型P=0.45時、邏輯回歸法P=0.35、指數(shù)疊加法Y=5則認為有霾出現(xiàn)。通過與實況比較,得出三種模型預報評分(見表6)。由表可以看出,邏輯回歸法同樣在預報檢驗中效果最好,指數(shù)疊加法預報效果較差,其中邏輯回歸法TS評分及準確率分別為50.2%、80.3%,為三種模型中最高,空報率在三種模型中也最低。事件概率法空報率最高,但其漏報率卻是三種模型中最低,只有21.4%。而指數(shù)疊加法的預報效果則介于上述者之間。總體而言三種模型都比較穩(wěn)定,盡管在預報中的準確性較回報檢驗有所下降,但三種霾天氣的潛勢預報模型仍在實際預報有較好地參考價值。

    表6 三種模型試預報效果Tab.6 Prediction effect of three models (%)

    2.6.3 誤差分析

    (1)三種潛勢預報模型對預報的能起到參考作用,但沒有對天氣類型進行劃分,沒有結合霾天氣模型概念,因此對霾預報的準確性有一定影響。

    (2)用于預報所用的次日08時探空資料、日降水是用20時起報的12小時ECMWF細網(wǎng)格數(shù)值預報產品,其數(shù)值與建模時使用的實況探空資料有一定差異,這對模型的準確性有一定影響。此外數(shù)值預報產品,其產品質量一定程度也影響了預報準確性。

    (3)對預報因子進行二分值轉換時,一些因子的閾值界限不像PM2.5、降水等因子那樣劃分明顯,它們出現(xiàn)霾和不出現(xiàn)霾的值范圍重復多,部分因子劃分閾值時所占比例低,影響了預報模式的準確性。

    (4)資料時間序列短,參與建模的年數(shù)只有4年,做出的預報模型代表性還不是很好。

    3 結 論

    本文分析了霾天氣潛勢預報的幾種方法,重點介紹了基于事件概率的逐步回歸、邏輯回歸、指數(shù)疊加等三種方法,并用這三種方法建立了潛勢預報模型,通過回報和試預報檢驗結果來看,三種方法在實際預報中都有較好的參考作用,預報模型穩(wěn)定。對比各自檢驗效果來看,邏輯回歸法在回報檢驗和試預報檢驗中,無論是TS預報評分還是空報率、準確性都優(yōu)于事件概率法和指數(shù)疊加法,但計算較復雜。指數(shù)疊加法預報回報檢驗和試預報檢驗效果差別小,穩(wěn)定性更高,納入預報的因子少,計算時只需將二值化后預報因子簡單相加,運算過程簡單,易操作。綜合來看,三種方法TS預報評分都在50%左右,準確率均超過70%,在預報上都具有可行性、有效性。

    但受限于所用資料時間序列短、資料不全、預報因子篩選欠準確、霾天氣學成因分析缺乏,對霧霾區(qū)分方法手段缺乏,特別是霾維持時間短的現(xiàn)象及霧霾同存現(xiàn)象缺乏足夠分析,導致在預報模型的空報率、漏報率較高,預報的實效性較短。

    猜你喜歡
    潛勢漏報天氣
    賀蘭山沿山地區(qū)強對流潛勢預報方法研究
    臨渙選煤廠浮選效果評價與尾煤再選潛勢分析
    選煤技術(2022年1期)2022-04-19 11:15:00
    天氣冷了,就容易抑郁嗎?
    氨對甲苯二次有機氣溶膠生成潛勢影響的實驗研究
    能源工程(2021年2期)2021-07-21 08:39:54
    誰是天氣之子
    盛暑天氣,覓得書中一味涼
    文苑(2020年7期)2020-08-12 09:36:38
    Weather(天氣)
    “適者生存”?
    雜文月刊(2017年24期)2017-03-03 07:42:55
    各類氣體報警器防誤報漏報管理系統(tǒng)的應用
    傳染病漏報原因分析及對策
    国产欧美日韩一区二区精品| 久久天堂一区二区三区四区| 午夜两性在线视频| 99精品在免费线老司机午夜| 欧美丝袜亚洲另类 | 欧美日韩精品网址| 欧洲精品卡2卡3卡4卡5卡区| 日本a在线网址| 久久精品国产亚洲av香蕉五月| 国产亚洲av高清不卡| 一级毛片高清免费大全| 91大片在线观看| 午夜老司机福利片| 搡老妇女老女人老熟妇| 女性生殖器流出的白浆| 男女下面进入的视频免费午夜 | 午夜老司机福利片| 午夜福利在线观看吧| 亚洲九九香蕉| 欧美日本中文国产一区发布| 欧美黑人精品巨大| 国产一区二区三区视频了| 亚洲情色 制服丝袜| 久久久久久久久中文| 两个人看的免费小视频| 国产精品久久电影中文字幕| 色尼玛亚洲综合影院| 夜夜躁狠狠躁天天躁| aaaaa片日本免费| 日韩精品中文字幕看吧| 国产精品电影一区二区三区| www.熟女人妻精品国产| 中文字幕色久视频| 首页视频小说图片口味搜索| 国产区一区二久久| 一区在线观看完整版| 大香蕉久久成人网| 最近最新免费中文字幕在线| 国产成人影院久久av| 给我免费播放毛片高清在线观看| 精品欧美国产一区二区三| 制服诱惑二区| 欧美人与性动交α欧美精品济南到| 久久午夜亚洲精品久久| 亚洲五月婷婷丁香| 一级,二级,三级黄色视频| 成人手机av| 日韩一卡2卡3卡4卡2021年| 亚洲欧美精品综合久久99| 校园春色视频在线观看| 在线播放国产精品三级| 在线av久久热| 久久人妻av系列| 啦啦啦韩国在线观看视频| 国产亚洲精品一区二区www| 一级a爱视频在线免费观看| 老司机靠b影院| 国产亚洲精品第一综合不卡| avwww免费| 在线免费观看的www视频| 亚洲性夜色夜夜综合| 亚洲aⅴ乱码一区二区在线播放 | 色老头精品视频在线观看| 国产乱人伦免费视频| 悠悠久久av| 久久久久亚洲av毛片大全| 日韩精品免费视频一区二区三区| 欧美激情 高清一区二区三区| 久久精品影院6| 免费女性裸体啪啪无遮挡网站| 91麻豆精品激情在线观看国产| 亚洲国产精品sss在线观看| 久久亚洲精品不卡| 99热只有精品国产| 黄色视频,在线免费观看| 久久精品国产亚洲av高清一级| 桃红色精品国产亚洲av| 欧美色欧美亚洲另类二区 | 搡老岳熟女国产| 一区二区日韩欧美中文字幕| 88av欧美| 亚洲欧美一区二区三区黑人| 午夜福利在线观看吧| 欧美成狂野欧美在线观看| 国产麻豆69| 欧美国产日韩亚洲一区| 99久久99久久久精品蜜桃| 激情在线观看视频在线高清| 亚洲欧美激情在线| 老熟妇仑乱视频hdxx| 午夜福利欧美成人| av有码第一页| 免费搜索国产男女视频| 性少妇av在线| 欧美久久黑人一区二区| 在线十欧美十亚洲十日本专区| 999精品在线视频| 日本欧美视频一区| 19禁男女啪啪无遮挡网站| 国产精品久久久久久人妻精品电影| tocl精华| 亚洲精品国产一区二区精华液| 真人做人爱边吃奶动态| 成人精品一区二区免费| 午夜老司机福利片| 老司机午夜福利在线观看视频| 亚洲一区高清亚洲精品| 91大片在线观看| av电影中文网址| 动漫黄色视频在线观看| 国产真人三级小视频在线观看| 女生性感内裤真人,穿戴方法视频| 亚洲精华国产精华精| 97超级碰碰碰精品色视频在线观看| 国内久久婷婷六月综合欲色啪| 国产成人精品久久二区二区免费| 又紧又爽又黄一区二区| 亚洲成人精品中文字幕电影| 久久精品影院6| 亚洲精华国产精华精| 免费不卡黄色视频| 黄色a级毛片大全视频| 欧美精品啪啪一区二区三区| 日韩精品青青久久久久久| 亚洲性夜色夜夜综合| 久久草成人影院| 色婷婷久久久亚洲欧美| 久热这里只有精品99| 麻豆av在线久日| 9色porny在线观看| 午夜精品国产一区二区电影| 两个人视频免费观看高清| √禁漫天堂资源中文www| 亚洲欧美日韩另类电影网站| 岛国视频午夜一区免费看| 操出白浆在线播放| 日韩国内少妇激情av| 国产成人免费无遮挡视频| 亚洲色图 男人天堂 中文字幕| 老司机午夜十八禁免费视频| 国产精品永久免费网站| 法律面前人人平等表现在哪些方面| 久久久久国产一级毛片高清牌| 天天添夜夜摸| 最近最新免费中文字幕在线| 精品久久蜜臀av无| 级片在线观看| √禁漫天堂资源中文www| 91九色精品人成在线观看| 麻豆久久精品国产亚洲av| avwww免费| 亚洲九九香蕉| bbb黄色大片| 中文字幕最新亚洲高清| 久久精品亚洲精品国产色婷小说| 纯流量卡能插随身wifi吗| 亚洲激情在线av| 满18在线观看网站| 欧美色欧美亚洲另类二区 | 午夜影院日韩av| 俄罗斯特黄特色一大片| 国产精品秋霞免费鲁丝片| 亚洲精品av麻豆狂野| 黄网站色视频无遮挡免费观看| 黑人欧美特级aaaaaa片| 50天的宝宝边吃奶边哭怎么回事| 亚洲av成人一区二区三| 亚洲欧美精品综合一区二区三区| 乱人伦中国视频| 亚洲成av片中文字幕在线观看| 99国产精品一区二区蜜桃av| 性少妇av在线| 麻豆一二三区av精品| 久久精品91蜜桃| 久久精品国产亚洲av香蕉五月| 两个人视频免费观看高清| 国产亚洲欧美98| 亚洲国产欧美网| 少妇粗大呻吟视频| 国产一卡二卡三卡精品| 欧美在线一区亚洲| 亚洲成a人片在线一区二区| 欧美性长视频在线观看| 免费高清视频大片| 在线av久久热| 少妇粗大呻吟视频| 1024视频免费在线观看| 91精品国产国语对白视频| 波多野结衣av一区二区av| 1024视频免费在线观看| 国产伦一二天堂av在线观看| 久久精品影院6| 国产单亲对白刺激| 一区二区三区高清视频在线| cao死你这个sao货| 麻豆一二三区av精品| 亚洲国产看品久久| 欧美人与性动交α欧美精品济南到| 一本大道久久a久久精品| 久久精品国产清高在天天线| 老司机午夜福利在线观看视频| 激情视频va一区二区三区| 大陆偷拍与自拍| 丰满人妻熟妇乱又伦精品不卡| 国产精品综合久久久久久久免费 | 亚洲av电影不卡..在线观看| 国产精品久久久久久亚洲av鲁大| 禁无遮挡网站| 一边摸一边抽搐一进一出视频| 91精品三级在线观看| 欧美色欧美亚洲另类二区 | 女警被强在线播放| 亚洲色图av天堂| 国产午夜精品久久久久久| 老司机午夜福利在线观看视频| 精品一区二区三区视频在线观看免费| 国产成人精品无人区| 人人澡人人妻人| 国产亚洲精品第一综合不卡| 一夜夜www| 亚洲欧美精品综合久久99| 高潮久久久久久久久久久不卡| 国产精品国产高清国产av| 少妇粗大呻吟视频| 成人免费观看视频高清| 午夜老司机福利片| 99久久99久久久精品蜜桃| 嫩草影视91久久| 51午夜福利影视在线观看| 可以免费在线观看a视频的电影网站| 两性午夜刺激爽爽歪歪视频在线观看 | 村上凉子中文字幕在线| 精品国产乱码久久久久久男人| 国产精品亚洲av一区麻豆| 亚洲国产日韩欧美精品在线观看 | 99精品在免费线老司机午夜| 日韩精品青青久久久久久| 久久精品人人爽人人爽视色| 精品一区二区三区四区五区乱码| 女人爽到高潮嗷嗷叫在线视频| 日韩欧美国产一区二区入口| 少妇裸体淫交视频免费看高清 | 天天躁夜夜躁狠狠躁躁| 国产三级黄色录像| 久久久精品国产亚洲av高清涩受| 久久精品亚洲精品国产色婷小说| 国产在线观看jvid| 国产高清有码在线观看视频 | 亚洲精品中文字幕在线视频| 亚洲精品美女久久av网站| 99精品久久久久人妻精品| 久久国产亚洲av麻豆专区| 麻豆久久精品国产亚洲av| 久久久久久大精品| 亚洲成a人片在线一区二区| 亚洲最大成人中文| 日本欧美视频一区| 欧美一级毛片孕妇| 看免费av毛片| 每晚都被弄得嗷嗷叫到高潮| 99国产精品免费福利视频| www.www免费av| 一级毛片女人18水好多| 纯流量卡能插随身wifi吗| 91成年电影在线观看| 首页视频小说图片口味搜索| 老熟妇乱子伦视频在线观看| 午夜免费观看网址| 两性夫妻黄色片| 丁香六月欧美| 无遮挡黄片免费观看| 久久天躁狠狠躁夜夜2o2o| 人妻丰满熟妇av一区二区三区| 波多野结衣高清无吗| 免费看美女性在线毛片视频| 99久久综合精品五月天人人| 少妇熟女aⅴ在线视频| 两个人免费观看高清视频| 午夜福利在线观看吧| 人妻丰满熟妇av一区二区三区| 午夜精品国产一区二区电影| 嫩草影院精品99| e午夜精品久久久久久久| 亚洲av五月六月丁香网| 丰满的人妻完整版| 黄片大片在线免费观看| 一级毛片精品| 最好的美女福利视频网| 亚洲五月天丁香| 9色porny在线观看| 精品电影一区二区在线| 久久久久久久午夜电影| 色综合欧美亚洲国产小说| 少妇 在线观看| www.999成人在线观看| 国产欧美日韩一区二区精品| 午夜日韩欧美国产| 午夜福利18| 母亲3免费完整高清在线观看| 高清在线国产一区| 亚洲专区中文字幕在线| 日韩欧美免费精品| 欧美中文综合在线视频| 亚洲精品国产区一区二| 国产伦人伦偷精品视频| 色尼玛亚洲综合影院| 欧美日韩黄片免| 99在线人妻在线中文字幕| 免费不卡黄色视频| 亚洲一区中文字幕在线| 黄片小视频在线播放| 国产精品秋霞免费鲁丝片| 成年版毛片免费区| 国产精品一区二区三区四区久久 | 男人的好看免费观看在线视频 | cao死你这个sao货| 成人18禁高潮啪啪吃奶动态图| 国产欧美日韩一区二区三区在线| 亚洲av第一区精品v没综合| 国产精品一区二区免费欧美| 国产男靠女视频免费网站| 亚洲av熟女| 丰满人妻熟妇乱又伦精品不卡| 国产亚洲精品久久久久久毛片| 亚洲久久久国产精品| 在线观看舔阴道视频| 操出白浆在线播放| 日日爽夜夜爽网站| 丰满的人妻完整版| 国产熟女午夜一区二区三区| 欧美精品啪啪一区二区三区| 国产一区二区三区综合在线观看| 亚洲av电影不卡..在线观看| 午夜影院日韩av| 给我免费播放毛片高清在线观看| 丝袜美足系列| 欧美老熟妇乱子伦牲交| 美女扒开内裤让男人捅视频| 国产男靠女视频免费网站| 桃红色精品国产亚洲av| 一级黄色大片毛片| 久久亚洲真实| 97人妻天天添夜夜摸| 亚洲九九香蕉| 久久久久久久久中文| 黄片小视频在线播放| 老司机福利观看| 男女午夜视频在线观看| 很黄的视频免费| 999精品在线视频| 满18在线观看网站| 又黄又粗又硬又大视频| 国产一区二区三区在线臀色熟女| 精品国内亚洲2022精品成人| 国产精品久久久久久精品电影 | 女同久久另类99精品国产91| av超薄肉色丝袜交足视频| 欧美 亚洲 国产 日韩一| 成人亚洲精品一区在线观看| 91精品国产国语对白视频| 中文字幕人妻丝袜一区二区| 免费观看精品视频网站| 国产欧美日韩精品亚洲av| 中文字幕最新亚洲高清| 亚洲第一青青草原| 国产aⅴ精品一区二区三区波| 精品一品国产午夜福利视频| 美女扒开内裤让男人捅视频| 国产av一区二区精品久久| 中出人妻视频一区二区| 免费少妇av软件| 天堂影院成人在线观看| 伊人久久大香线蕉亚洲五| 禁无遮挡网站| 激情在线观看视频在线高清| 人成视频在线观看免费观看| 色老头精品视频在线观看| 国产亚洲精品久久久久5区| 怎么达到女性高潮| 99在线人妻在线中文字幕| av天堂久久9| 精品人妻1区二区| 日韩高清综合在线| 少妇熟女aⅴ在线视频| 国内毛片毛片毛片毛片毛片| 一本综合久久免费| 香蕉丝袜av| 搡老熟女国产l中国老女人| 国产成人精品在线电影| 一级毛片高清免费大全| 久久久国产精品麻豆| 可以在线观看毛片的网站| 日韩三级视频一区二区三区| 热re99久久国产66热| av有码第一页| 色播在线永久视频| 99国产综合亚洲精品| av电影中文网址| 国产精品精品国产色婷婷| 色播在线永久视频| 国产亚洲精品久久久久5区| 国产免费男女视频| 亚洲黑人精品在线| 黄色成人免费大全| 18禁裸乳无遮挡免费网站照片 | ponron亚洲| 别揉我奶头~嗯~啊~动态视频| 国产亚洲精品久久久久久毛片| 1024视频免费在线观看| 国产亚洲精品久久久久久毛片| 桃红色精品国产亚洲av| 亚洲国产欧美日韩在线播放| 久久午夜亚洲精品久久| 亚洲成人久久性| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产欧美日韩在线播放| 国产精品秋霞免费鲁丝片| www.自偷自拍.com| 桃色一区二区三区在线观看| 国产成人系列免费观看| 校园春色视频在线观看| 1024香蕉在线观看| 国内精品久久久久精免费| 亚洲精品美女久久久久99蜜臀| 窝窝影院91人妻| 国产一区二区激情短视频| 国产成人免费无遮挡视频| 成人免费观看视频高清| 咕卡用的链子| 亚洲七黄色美女视频| 亚洲一区中文字幕在线| 一级,二级,三级黄色视频| 少妇 在线观看| 18禁裸乳无遮挡免费网站照片 | 国产亚洲精品第一综合不卡| 亚洲av电影不卡..在线观看| 18禁国产床啪视频网站| 午夜免费鲁丝| 女性生殖器流出的白浆| 韩国av一区二区三区四区| 男女之事视频高清在线观看| 免费无遮挡裸体视频| www.www免费av| 老汉色∧v一级毛片| 精品熟女少妇八av免费久了| 午夜免费成人在线视频| 久久久久国产一级毛片高清牌| netflix在线观看网站| 91麻豆精品激情在线观看国产| 午夜精品在线福利| 精品久久久久久久久久免费视频| 女性被躁到高潮视频| 神马国产精品三级电影在线观看 | 在线av久久热| 又黄又爽又免费观看的视频| 国产精品香港三级国产av潘金莲| 好男人在线观看高清免费视频 | 熟女少妇亚洲综合色aaa.| 成年版毛片免费区| 欧美+亚洲+日韩+国产| 久久精品影院6| 国产熟女午夜一区二区三区| 免费一级毛片在线播放高清视频 | 欧美成人一区二区免费高清观看 | 亚洲精品国产一区二区精华液| 黑人巨大精品欧美一区二区mp4| 免费高清视频大片| 97超级碰碰碰精品色视频在线观看| 国产乱人伦免费视频| 国产一卡二卡三卡精品| 久久香蕉激情| 国产男靠女视频免费网站| 欧美成人免费av一区二区三区| 婷婷六月久久综合丁香| 我的亚洲天堂| 嫩草影院精品99| 久久精品亚洲熟妇少妇任你| 在线观看午夜福利视频| 天天躁夜夜躁狠狠躁躁| 婷婷丁香在线五月| 亚洲成人国产一区在线观看| 日日摸夜夜添夜夜添小说| 91av网站免费观看| 一二三四在线观看免费中文在| 青草久久国产| 亚洲av成人不卡在线观看播放网| 久久久久国产一级毛片高清牌| 他把我摸到了高潮在线观看| 麻豆一二三区av精品| 制服丝袜大香蕉在线| 国语自产精品视频在线第100页| 两性夫妻黄色片| 97碰自拍视频| 成人亚洲精品av一区二区| 91麻豆精品激情在线观看国产| 国产1区2区3区精品| 日本a在线网址| 麻豆成人av在线观看| 欧美日韩福利视频一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久精品国产欧美久久久| 老汉色av国产亚洲站长工具| 亚洲一区二区三区不卡视频| 亚洲精华国产精华精| 久久久久久人人人人人| 少妇的丰满在线观看| 1024香蕉在线观看| 精品国产国语对白av| 电影成人av| 欧美亚洲日本最大视频资源| 夜夜夜夜夜久久久久| 久久精品aⅴ一区二区三区四区| 人妻丰满熟妇av一区二区三区| or卡值多少钱| 97人妻天天添夜夜摸| 一级作爱视频免费观看| 亚洲国产精品sss在线观看| 黄色丝袜av网址大全| 视频区欧美日本亚洲| 久久精品亚洲熟妇少妇任你| 久久午夜亚洲精品久久| 国产成人系列免费观看| 50天的宝宝边吃奶边哭怎么回事| 精品国产乱码久久久久久男人| 亚洲欧美精品综合一区二区三区| 99精品久久久久人妻精品| 亚洲欧美日韩无卡精品| av天堂久久9| 国产精品自产拍在线观看55亚洲| 国产乱人伦免费视频| 免费在线观看黄色视频的| 丁香欧美五月| 人人澡人人妻人| 国产av精品麻豆| 亚洲专区国产一区二区| 免费在线观看亚洲国产| 熟女少妇亚洲综合色aaa.| 少妇被粗大的猛进出69影院| 国产精品野战在线观看| 国产精品久久久久久人妻精品电影| 亚洲国产欧美日韩在线播放| 午夜福利成人在线免费观看| 99在线人妻在线中文字幕| 757午夜福利合集在线观看| 男人操女人黄网站| 美女午夜性视频免费| 精品国内亚洲2022精品成人| 亚洲av第一区精品v没综合| 国产精品亚洲一级av第二区| 一区二区三区精品91| 亚洲美女黄片视频| 淫妇啪啪啪对白视频| 叶爱在线成人免费视频播放| 欧美成人午夜精品| 久久久久久亚洲精品国产蜜桃av| 国产91精品成人一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 久久精品成人免费网站| 国产精品影院久久| 久久久久精品国产欧美久久久| 国产一区二区三区在线臀色熟女| 在线观看免费视频日本深夜| 亚洲欧美激情综合另类| 精品国产美女av久久久久小说| 看黄色毛片网站| 亚洲精华国产精华精| 国产精品久久久久久亚洲av鲁大| 国内久久婷婷六月综合欲色啪| videosex国产| 变态另类成人亚洲欧美熟女 | 精品国产一区二区三区四区第35| 91大片在线观看| 丝袜美腿诱惑在线| 久热爱精品视频在线9| 国产极品粉嫩免费观看在线| 两个人免费观看高清视频| 在线观看舔阴道视频| 欧美乱妇无乱码| 啦啦啦免费观看视频1| 麻豆一二三区av精品| 99riav亚洲国产免费| 久久久国产欧美日韩av| 免费在线观看日本一区| cao死你这个sao货| 999久久久国产精品视频| 中文字幕最新亚洲高清| 国产伦人伦偷精品视频| 12—13女人毛片做爰片一| 欧美日韩亚洲国产一区二区在线观看| 国产三级黄色录像| 国产精品免费视频内射| 亚洲成国产人片在线观看| 国产精品二区激情视频| 女性被躁到高潮视频| 精品无人区乱码1区二区| 91麻豆精品激情在线观看国产| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲人成伊人成综合网2020| 免费av毛片视频| 日韩国内少妇激情av| 99久久国产精品久久久| 欧美午夜高清在线| 在线十欧美十亚洲十日本专区| 黄色片一级片一级黄色片| av视频免费观看在线观看| 91精品国产国语对白视频| 母亲3免费完整高清在线观看| 亚洲精品在线观看二区| 国产在线精品亚洲第一网站| 免费搜索国产男女视频| av天堂在线播放| 国产高清videossex| 免费观看人在逋| 午夜视频精品福利| 最近最新中文字幕大全电影3 |