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

    “7.18”山東暴雨過程分析III:WRF模式邊界層參數(shù)化方案對物理量場的影響

    2016-09-08 06:12:20譚偉才文映方李啟華
    安徽農(nóng)業(yè)科學(xué) 2016年20期
    關(guān)鍵詞:散度邊界層對流

    譚偉才,文映方,李啟華

    (1.75839部隊,廣東廣州 510510;2.貴州省貴陽市氣象局,貴州貴陽 550001;3.解放軍理工大學(xué),江蘇南京 211101)

    ?

    “7.18”山東暴雨過程分析III:WRF模式邊界層參數(shù)化方案對物理量場的影響

    譚偉才1,文映方2,李啟華3

    (1.75839部隊,廣東廣州 510510;2.貴州省貴陽市氣象局,貴州貴陽 550001;3.解放軍理工大學(xué),江蘇南京 211101)

    [目的]分析不同邊界層參數(shù)化方案對此次暴雨過程中物理量場的影響。[方法]利用WRF模式3.4版本,選取不同邊界層參數(shù)化方案對“7.18”山東暴雨過程進(jìn)行敏感性試驗,分析不同邊界層參數(shù)化方案對此次暴雨過程中垂直速度場、水汽通量散度場、相對濕度、對流有效位能和邊界層高度等物理量場的影響。[結(jié)果]不同的邊界層參數(shù)化方案對于暴雨過程中垂直速度場的空間分布及其中心強(qiáng)度具有顯著的影響,從而使得降水中心的分布和強(qiáng)度發(fā)生變化。盡管不同邊界層參數(shù)化方案模擬得到的水汽通量散度在水平分布與實際結(jié)果具有較好的一致性,但其在垂直分布、中心位置及其強(qiáng)度上存在明顯的差異。不同邊界層參數(shù)化方案引起的邊界層高度和對流有效位能的差異與其引起的降水分布差異直接相關(guān)。[結(jié)論]WRF模式中不同邊界層參數(shù)化方案對暴雨過程中不同物理量場模擬效果的影響較大,選擇合適的邊界層參數(shù)化方案能顯著提高對物理量場的模擬效果。

    暴雨;WRF模式;敏感性試驗;邊界層參數(shù)化方案;降水

    在大氣中,邊界層是熱量和水汽源、動量匯,它通過湍流垂直運(yùn)動將熱量和水汽輸送到大氣上層,從而對降水的分布和降水量的大小產(chǎn)生重要的影響。因此研究WRF模式中不同邊界層參數(shù)化方案對于降水模擬的影響是非常必要的。長期以來,國內(nèi)外許多氣象學(xué)者采用不同的數(shù)值模式研究了不同的邊界層參數(shù)化對降水模擬結(jié)果的影響,并取得了一定的研究成果[1-5]。陳靜等[5]利用中尺度非靜力MM5模式,深入分析了對流參數(shù)化方案在中尺度暴雨預(yù)報中的作用,結(jié)果表明,邊界層方案和對流參數(shù)化方案對產(chǎn)生暴雨的3個基本條件(水汽通量散度、垂直速度、不穩(wěn)定層結(jié))的影響很明顯。蔡薌寧等[6]采用非靜力中尺度模式MM5著重研究了不同邊界層參數(shù)化方案對雨量中心強(qiáng)度、雨區(qū)分布的影響,結(jié)果表明,對于不同的邊界層參數(shù)化方案,垂直速度場、水汽通量散度場、渦度場、水平風(fēng)場的散度、θse場均表現(xiàn)出不同的特征。筆者利用WRF模式3.4版本,選取不同邊界層參數(shù)化方案對2007年山東“7.18”暴雨過程進(jìn)行敏感性試驗,分析不同邊界層參數(shù)化方案對此次暴雨過程中物理量場的影響,探討不同邊界層方案造成模擬降水差異的原因和機(jī)理。

    1 模式方案設(shè)計和資料選取

    利用WRF模式3.4版本,選取YSU邊界層參數(shù)化方案、MYJ邊界層參數(shù)化方案、ACM2邊界層參數(shù)化方案、MYNN2.5邊界層參數(shù)化方案、Boulac邊界層參數(shù)化方案以及不加任何邊界層參數(shù)化方案對此次暴雨過程進(jìn)行敏感性試驗。模擬采用單向兩層嵌套,模式水平分辨率分別為36、18 km。垂直方向分為30層,層頂氣壓為100 hPa。初值選用時間間隔為6 h,分辨率為1°×1°的 NCEP 資料,積分時間均為2007年7月18日02:00~19日14:00,共36 h。微物理過程均為Lin方案,積云對流過程均為Kain-Fritsch方案,近地面過程均為Monin-Obukhov方案。

    2 結(jié)果與分析

    2.1不同邊界層參數(shù)化方案對垂直速度場的影響從圖1可以看出,在MYJ方案的模擬結(jié)果中,強(qiáng)上升中心分布在115°~117° E;在117.5°~118.5° E和119°~120° E的對流層高層分別存在2個較弱的上升中心,其中位于西側(cè)的強(qiáng)上升中心向上伸展至100 hPa左右,并在其右側(cè)存在下沉氣流,形成典型的中尺度對流體的垂直運(yùn)動結(jié)構(gòu)。在YSU方案的模擬結(jié)果中,西側(cè)的強(qiáng)上升中心明顯分成了兩部分,在116°~116.5° E 400 hPa出現(xiàn)了一個新的上升中心;與MYJ方案相比,其位于117.5°~119°E的上升中心明顯更加強(qiáng)盛,不僅中心量值更大,垂直伸展和水平范圍均強(qiáng)過MYJ方案;而其東側(cè)的上升運(yùn)動中心較弱,位置更加偏于東側(cè)。與YSU方案類似,ACM2方案也在115°~116.5° E模擬出了2個上升運(yùn)動中心,但其強(qiáng)度明顯較弱,同時其位置更加偏西。在MYNN2.5方案的模擬結(jié)果中,位于116°E附近的上升運(yùn)動中心較前3個方案偏強(qiáng),表現(xiàn)在不僅其中心垂直速度增大并垂直伸展達(dá)100 hPa,水平范圍也擴(kuò)大;與ACM2方案類似,在118.5°E低層700 hPa附近,MYNN2.5方案也模擬出了一個較小的上升中心。在Boulac方案中,其上升運(yùn)動中心的水平范圍是所有方案中最大的,從115.5° E一直延伸至119° E左右,垂直速度大值中心分布也最為密集;西側(cè)主要的上升運(yùn)動伸展達(dá)100 hPa,而在東側(cè)800~200 hPa零星地存在多個上升運(yùn)動中心。在沒有使用邊界層參數(shù)化方案的模擬結(jié)果中,116°~117.5°E的垂直運(yùn)動中心及其強(qiáng)度與使用了邊界層方案的結(jié)果基本一致,但其垂直伸展范圍明顯更?。辉?18.5°~119.5°E模擬出了2個上升運(yùn)動中心,其中高層的位于300 hPa左右,而低層的垂直上升運(yùn)動中心的范圍向下伸展到模式低層,明顯比使用了邊界層方案的結(jié)果偏低。綜上所述,不同的邊界層參數(shù)化方案所得到的垂直運(yùn)動分布情況存在較大區(qū)別,從而直接影響降水以及暴雨中心的分布和強(qiáng)度情況。

    注:a.MYJ方案;b.YSU方案;c.ACM2方案;d.MYNN2.5方案;e.Boulac方案;f.none PBL方案。實等值線表示垂直速度(m/s);虛線表示等壓面。Note:a.MYJ scheme;b.YSU scheme;c.ACM2 scheme;d.MYNN2.5 scheme;e.Boulac scheme;f.none PBL scheme.Solid contour indicates vertical velocity(m/s);dotted line represents the isobaric surface.圖1 2007年7月18日20:00不同邊界層參數(shù)化方案模擬的垂直速度沿36.5°N的經(jīng)向-垂直剖面Fig.1 Meridional-vertical section of vertical velocity obtained by different PBL parameterization schemes along the parallel of 36.5 °N at 20:00 on July 18,2007

    2.2不同邊界層參數(shù)化方案對水汽通量散度場的影響從實際水汽通量散度沿36.5°N緯圈的垂直分布可以看出,主要的水汽輻合中心位于115°E附近950 hPa左右,而在119°E附近的低層和高層分別存在2個較小的水汽輻合中心。在MYJ方案的模擬結(jié)果中(圖2a),位于115°E的水汽輻合中心分成兩部分,且垂直高度比實際結(jié)果較低,位于800 hPa左右;而在119°E附近,水汽輻合中心不明顯,僅存在一條較弱的水汽輻合帶從貼地層伸展至700 hPa左右。在YSU方案的結(jié)果中(圖2b),位于115°E的水汽輻合中心強(qiáng)度較弱,其水平范圍也較窄;而位于119°E附近的水汽輻合中心在垂直方向伸展較低,分布也較零散。在ACM2方案的模擬結(jié)果中(圖2c),位于115°E的水汽輻合帶的垂直伸展范圍和強(qiáng)度與實際結(jié)果基本一致,但在水平方向上的范圍較大,且輻合中心分布也并不明顯;位于119°E的輻合帶比實際結(jié)果偏小,但強(qiáng)度偏強(qiáng)。在MYNN2.5方案的模擬結(jié)果(圖2d)中,位于115°E的水汽輻合帶的位置比實際結(jié)果偏東,其垂直伸展范圍也較低;在119°E附近的輻合帶分布與其他方案相似。在Boulac方案的結(jié)果中(圖2e),位于115°E的水汽輻合帶出現(xiàn)了多個輻合中心,但垂直范圍僅達(dá)800 hPa左右;位于119°E的輻合帶與其他方案大體一致,其中心分布較零散,且垂直伸展范圍較低。在不使用任何邊界層方案的結(jié)果中(圖2f),水汽輻合帶的分布上與實際結(jié)果的差異最大,其位于115°E和119°E附近的輻合中心的強(qiáng)度和位置均與實際結(jié)果存在較大差異,并在低層模擬出了2個較強(qiáng)的虛假水汽輻散帶。從以上分析可知,雖然采用不同邊界層參數(shù)化方案均大致模擬出了水汽通量散度在115°E和119°E附近的輻合中心,但其中心強(qiáng)度和中心垂直伸展范圍存在明顯的差異。

    2.3不同邊界層參數(shù)化方案對不穩(wěn)定層結(jié)的影響從圖3

    注:a.MYJ方案;b.YSU方案;c.ACM2方案;d.MYNN2.5方案;e.Boulac方案;f.none PBL方案。實等值線表示水汽通量散度[10-6g/(cm2·hPa·s)];虛線表示等壓面。Note:a.MYJ scheme;b.YSU scheme;c.ACM2 scheme;d.MYNN2.5 scheme;e.Boulac scheme;f.none PBL scheme.Solid contour indicates moisture flux divergence[10-6g/(cm2·hPa·s)];dotted line represents the isobaric surface.圖2 2007年7月18日20:00不同邊界層參數(shù)化方案模擬的水汽通量散度沿36.5°N的經(jīng)向-垂直剖面Fig.2 Meridional-vertical section of moisture flux divergence obtained by different PBL parameterization schemes along the parallel of 36.5 °N at 20:00 on July 18,2007

    可以看出,所有方案的模擬結(jié)果均在垂直方向上存在上干下濕的對流不穩(wěn)定層結(jié)分布;這種不穩(wěn)定層結(jié)主要分布在114°~115.5° E和116.5°~118° E。在MYJ方案的模擬結(jié)果中,西側(cè)(114°~115.5° E)的對流不穩(wěn)定層結(jié)主要位于對流層低層700~850 hPa,而東側(cè)(116.5°~118° E)的對流不穩(wěn)定層結(jié)主要位于對流層中層,在119° E 700 hPa附近也存在明顯的對流不穩(wěn)定層結(jié)。在YSU方案的模擬結(jié)果中,西側(cè)的相對濕度分布與MYJ方案較為相似,而東側(cè)的不穩(wěn)定層結(jié)在水平方向上具有更大的范圍。在ACM2方案的模擬結(jié)果中,西側(cè)的不穩(wěn)定層結(jié)強(qiáng)度遠(yuǎn)小于前2個方案,而在東側(cè),其水平范圍與YSU方案大體一致。在MYNN2.5方案的模擬結(jié)果中,西側(cè)的不穩(wěn)定層結(jié)強(qiáng)度依然較小,而東側(cè)的高空出現(xiàn)相對濕度低值中心。在Boulac方案的模擬結(jié)果中,西側(cè)的不穩(wěn)定層結(jié)位置相比其他方案明顯偏低,幾乎是從模式低層開始的,且在垂直方向上并不連續(xù)。由此可知,不同的邊界層參數(shù)化方案對暴雨過程中的不穩(wěn)定層結(jié)的垂直結(jié)構(gòu)分布具有重要影響,但從整體來看,不穩(wěn)定層結(jié)均是從對流層中低層開始發(fā)展的,這為暴雨提供了較好的不穩(wěn)定條件。

    注:a.MYU方案;b.YSU方案;c.ACM2方案;d.MYNN2.5方案;e.Boulac方案;f.none PBL方案。實等值線表示相對濕度(%);虛線表示等壓面。Note:a.MYJ scheme;b.YSU scheme;c.ACM2 scheme;d.MYNN2.5 scheme;e.Boulac scheme;f.none PBL scheme.Solid contour indicates relative humidity;dotted line represents the isobaric surface.圖3 2007年7月18日20:00不同邊界層參數(shù)化方案模擬的相對濕度沿36.5°N的經(jīng)向-垂直剖面Fig.3 Meridional-vertical section of relative humidity obtained by different PBL parameterization schemes along the parallel of 36.5°N at 20:00 on July 18,2007

    從圖4可以看出,在所有方案的模擬結(jié)果中,濟(jì)南市地面對流有效位能(CAPE)在18日14:00均達(dá)到了極值,然后隨著暴雨的發(fā)生,對流不穩(wěn)定能量逐漸釋放;與CAPE釋放相對應(yīng),邊界層高度也存在相似的變化趨勢。盡管不同方案中CAPE開始釋放的時間基本一致,但能量極值的大小和釋放的快慢程度存在較大差異。CAPE積累最大的是MYJ方案,最少的是ACM2方案,而不使用任何邊界層方案時其比ACM2方案還低。盡管ACM2方案能量積累最少,但在18日20:00~23:00(即模擬的暴雨峰值出現(xiàn)時段),對流不穩(wěn)定能量釋放最多的也是ACM2方案。邊界層高度的變化也存在較大差異,在18日08:00~14:00不穩(wěn)定能量積累期,所有方案模擬的邊界層高度均增加,當(dāng)對流不穩(wěn)定能量開始釋放后(14:00以后),不同邊界層方案模擬得到的邊界層高度均出現(xiàn)了急劇的下降;而在暴雨結(jié)束后,邊界層高度又開始逐漸恢復(fù)。在暴雨整個過程中,MYJ方案的邊界層高度一直維持著最高的狀態(tài);而當(dāng)CAPE積累到極值時,MYNN2.5方案的邊界層高度是最低的;在暴雨最大降水時段(18日20:00~23:00),ACM2方案和Boulac方案的邊界層高度是最低的。

    圖4 不同邊界層方案模擬的濟(jì)南市上空CAPE(a)和邊界層高度(b)隨時間變化Fig.4 Time-variation of CAPE(a)and boundary layer height(b)obtained by different PBL parameterization schemesover the Jinang City

    3 結(jié)論

    該研究利用WRF模式3.4版本,選取不同邊界層參數(shù)化方案對“7.18”山東暴雨過程進(jìn)行敏感性試驗,分析不同邊界層參數(shù)化方案對此次暴雨過程中垂直速度場、水汽通量散度場、相對濕度、對流有效位能和邊界層高度等物理量場的影響,得到以下主要結(jié)論:

    (1)不同的邊界層參數(shù)化方案對于暴雨過程中垂直速度場的模擬存在很大區(qū)別,主要體現(xiàn)在垂直運(yùn)動在垂直方向和水平方向上的伸展范圍及其中心強(qiáng)度的差異上,從而直接影響降水以及暴雨中心的分布和強(qiáng)度情況。當(dāng)不使用邊界層參數(shù)化方案時,垂直上升運(yùn)動中心的位置明顯較使用參數(shù)化方案得到的垂直運(yùn)動分布偏低。

    (2)盡管不同邊界層參數(shù)化方案模擬的水汽通量散度在水平分布上與實際結(jié)果具有較好的一致性,但其在垂直分布、中心位置及其強(qiáng)度上存在明顯差異。其中,ACM2方案模擬的水汽通量散度分布上與實際結(jié)果最為接近,效果最好。而不使用邊界層方案的模擬結(jié)果與實際水汽通量散度分布之間的差異最大。表明邊界層內(nèi)的物理過程對水汽通量散度分布具有非常重要的影響。

    (3)不同邊界層方案模擬的垂直方向上的不穩(wěn)定層結(jié)分布均大體相似,其中YSU方案和MYJ方案的模擬結(jié)果較為接近,MYNN2.5和ACM2方案的模擬結(jié)果較為接近。不同邊界層方案在對流有效位能的積累和釋放過程上存在明顯的差異。對比暴雨過程中降水的強(qiáng)度變化與對流有效位能、邊界層高度變化可知,不同邊界層方案得到的邊界層高度和對流有效位能變化的差異與其引起的降水分布差異直接相關(guān)。

    [1] ELMOUSTAFA A M,F(xiàn)ARRES H N,ELFAWY M M.Analysis of the role of the planetary boundary layer schemes during a severe convective storm[J].Annales geophysicae,2004,22(6):1861-1874.

    [2] 左志燕.行星邊界層參數(shù)化對我國夏季降水的影響及其成因分析[D].北京:中國氣象科學(xué)研究院,2004:16.

    [3] 徐慧燕,朱業(yè),劉瑞,等.長江下游地區(qū)不同邊界層參數(shù)化方案的試驗研究[J].大氣科學(xué),2013,37(1):149-159.

    [5] 陳靜,薛紀(jì)善,顏宏.物理過程參數(shù)化方案對中尺度暴雨數(shù)值模擬影響的研究[J].氣象學(xué)報,2003,61(2):203-218.

    [6] 蔡薌寧,壽紹文,鐘青.邊界層參數(shù)化方案對暴雨數(shù)值模擬的影響[J].大氣科學(xué)學(xué)報,2006,29(3):364-370.

    Analysis of “7.18” Shandong Heavy Rainfall Process III: the Impact of WRF Model Boundary Layer Parameterization Scheme on Physical Quantity Field

    TAN Wei-cai1, WEN Ying-fang2, LI Qi-hua3

    (1.75839 Troops,Guangzhou,Guangdong 510510;2.Meteorological Bureau of Guiyang City,Guiyang,Guizhou 550001;3. PLA University of Science and Technology, Nanjing, Jiangsu 211101)

    [Objective] The aim was to analyze effects of different boundary layer parameterization schemes on physical quantity field in the heavy rainfall process. [Method] Using WRF mode 3.4 version, sensitivity test was conducted on “7.18” Shandong heavy rainfall process, effects of different boundary layer parameterization schemes on vertical velocity field, water vapor flux divergence field, relative humidity, convective available potential energy and boundary layer height were analyzed. [Result] The results showed that different boundary layer parameterization schemes have a significant impact on spatial distribution of vertical velocity and the strength of rainstorm center, so the location and strength of rainstorm center change. The simulated water vapor flux divergence accords with the actual results on the horizontal distribution, but its center location, vertical distribution and strength have significant differences. The convective boundary layer height and effective potential energy are different in models with boundary layer parameterization schemes, which directly relate to the difference of precipitation distribution. [Conclusion] Different boundary layer parameterization schemes in WRF mode have significant influence on simulation results of various physical quantity field in rainstorm process, using appropriate boundary layer parameterization scheme can improve simulation results of physical quantity field significantly.

    Rainstorm;WRF model; Sensitivity test; Boundary layer parameterization scheme; Precipitation

    譚偉才(1991- ),男,湖南婁底人,碩士,從事中尺度氣象學(xué)研究。

    2016-04-30

    S 165+.2;P 458.2

    A

    0517-6611(2016)20-180-04

    猜你喜歡
    散度邊界層對流
    齊口裂腹魚集群行為對流態(tài)的響應(yīng)
    帶勢加權(quán)散度形式的Grushin型退化橢圓算子的Dirichlet特征值的上下界
    基于HIFiRE-2超燃發(fā)動機(jī)內(nèi)流道的激波邊界層干擾分析
    具有部分BMO系數(shù)的非散度型拋物方程的Lorentz估計
    H型群上一類散度形算子的特征值估計
    H?rmander 向量場上散度型拋物方程弱解的Orlicz估計
    基于ANSYS的自然對流換熱系數(shù)計算方法研究
    二元驅(qū)油水界面Marangoni對流啟動殘余油機(jī)理
    一類具有邊界層性質(zhì)的二次奇攝動邊值問題
    非特征邊界的MHD方程的邊界層
    av在线老鸭窝| 丰满人妻一区二区三区视频av| 美女xxoo啪啪120秒动态图| 男插女下体视频免费在线播放| 久久草成人影院| 赤兔流量卡办理| 国产av码专区亚洲av| 男女下面进入的视频免费午夜| 乱人视频在线观看| 热99re8久久精品国产| 国产成人精品一,二区| 白带黄色成豆腐渣| 久久久久久久久大av| 两个人视频免费观看高清| 色视频www国产| 国产欧美日韩精品一区二区| 久久久亚洲精品成人影院| 又粗又硬又长又爽又黄的视频| 欧美成人一区二区免费高清观看| 国产又黄又爽又无遮挡在线| 国产女主播在线喷水免费视频网站 | 国产成人精品久久久久久| 最新中文字幕久久久久| 女人十人毛片免费观看3o分钟| 色视频www国产| 国产精品人妻久久久久久| 国内精品美女久久久久久| 97超视频在线观看视频| 国产淫片久久久久久久久| 亚洲人成网站高清观看| 欧美一区二区精品小视频在线| 亚洲精品乱久久久久久| 国产精品久久久久久精品电影小说 | 91aial.com中文字幕在线观看| 国产精品日韩av在线免费观看| 精品一区二区三区视频在线| 看黄色毛片网站| 久久热精品热| 国产精品女同一区二区软件| 欧美成人精品欧美一级黄| 特级一级黄色大片| 婷婷六月久久综合丁香| 18禁裸乳无遮挡免费网站照片| 久久久久久久久久久丰满| 亚洲国产精品合色在线| 久久婷婷人人爽人人干人人爱| 欧美zozozo另类| 老司机福利观看| 日韩国内少妇激情av| 国产午夜精品久久久久久一区二区三区| 欧美一区二区国产精品久久精品| 久99久视频精品免费| 欧美潮喷喷水| 久久99精品国语久久久| av视频在线观看入口| 中文字幕亚洲精品专区| 亚洲成人av在线免费| 99久久精品热视频| 老司机影院成人| 成人综合一区亚洲| 国产亚洲5aaaaa淫片| 激情 狠狠 欧美| 秋霞在线观看毛片| 国内精品一区二区在线观看| 国产成人免费观看mmmm| 成人av在线播放网站| 国产精品电影一区二区三区| 国产三级在线视频| 男女视频在线观看网站免费| 亚洲精品一区蜜桃| 成人毛片a级毛片在线播放| 久久人人爽人人爽人人片va| 国产视频内射| 欧美不卡视频在线免费观看| 五月玫瑰六月丁香| 少妇熟女欧美另类| 亚洲欧美成人精品一区二区| 91aial.com中文字幕在线观看| 日本免费a在线| 亚洲精品456在线播放app| 久久鲁丝午夜福利片| 欧美日韩精品成人综合77777| 久久6这里有精品| 亚洲av不卡在线观看| 久久久久久久久久久免费av| 免费不卡的大黄色大毛片视频在线观看 | 久久精品国产亚洲网站| 亚洲精华国产精华液的使用体验| 村上凉子中文字幕在线| 亚洲高清免费不卡视频| 成人欧美大片| 国产一区二区三区av在线| 成人午夜高清在线视频| 国产亚洲av片在线观看秒播厂 | 免费看a级黄色片| 美女xxoo啪啪120秒动态图| 超碰97精品在线观看| 少妇熟女aⅴ在线视频| eeuss影院久久| 69av精品久久久久久| 别揉我奶头 嗯啊视频| 99视频精品全部免费 在线| 日韩大片免费观看网站 | 国产高清不卡午夜福利| 少妇的逼好多水| 少妇熟女欧美另类| .国产精品久久| 免费电影在线观看免费观看| 在现免费观看毛片| 波多野结衣巨乳人妻| 全区人妻精品视频| 寂寞人妻少妇视频99o| 夜夜爽夜夜爽视频| 91精品一卡2卡3卡4卡| 国产精品一区二区性色av| 国产精品av视频在线免费观看| 日韩成人av中文字幕在线观看| 青春草国产在线视频| 国产一区有黄有色的免费视频 | 男人的好看免费观看在线视频| 菩萨蛮人人尽说江南好唐韦庄 | 午夜老司机福利剧场| 国产视频内射| 亚洲三级黄色毛片| 国语自产精品视频在线第100页| 天天躁夜夜躁狠狠久久av| 国产成人91sexporn| 精品免费久久久久久久清纯| 精品一区二区三区视频在线| 日韩欧美国产在线观看| 久久久久国产网址| 日韩欧美 国产精品| 亚洲欧洲国产日韩| 久久久久久久久大av| 亚洲精品国产成人久久av| 国产黄色小视频在线观看| 搡老妇女老女人老熟妇| 在线免费观看不下载黄p国产| 又黄又爽又刺激的免费视频.| 欧美不卡视频在线免费观看| 久久99蜜桃精品久久| 成年女人永久免费观看视频| 中文字幕免费在线视频6| 青春草视频在线免费观看| 免费av不卡在线播放| 国产一区二区三区av在线| 狂野欧美激情性xxxx在线观看| 大香蕉97超碰在线| 高清在线视频一区二区三区 | eeuss影院久久| 国国产精品蜜臀av免费| 精品免费久久久久久久清纯| 嫩草影院入口| 丰满少妇做爰视频| 真实男女啪啪啪动态图| 午夜福利在线在线| 美女cb高潮喷水在线观看| 午夜视频国产福利| 精品国产三级普通话版| 久久人人爽人人片av| 九草在线视频观看| 欧美性猛交黑人性爽| 国产v大片淫在线免费观看| 国产老妇女一区| 久久精品综合一区二区三区| 99久久成人亚洲精品观看| 男的添女的下面高潮视频| 一个人免费在线观看电影| 尤物成人国产欧美一区二区三区| 日日啪夜夜撸| av卡一久久| 国产免费又黄又爽又色| 免费播放大片免费观看视频在线观看 | 欧美成人免费av一区二区三区| 精品不卡国产一区二区三区| 欧美性猛交╳xxx乱大交人| 色尼玛亚洲综合影院| 成人毛片a级毛片在线播放| 一个人免费在线观看电影| videos熟女内射| 久久久午夜欧美精品| 亚洲国产精品成人综合色| 国产色婷婷99| av免费观看日本| 日韩高清综合在线| 日本熟妇午夜| 欧美一区二区国产精品久久精品| 精品99又大又爽又粗少妇毛片| 身体一侧抽搐| 久久6这里有精品| 夜夜爽夜夜爽视频| 精品久久久久久电影网 | 亚洲自拍偷在线| 亚洲五月天丁香| 亚洲av中文字字幕乱码综合| 国产精品久久久久久精品电影| 国产人妻一区二区三区在| 日本熟妇午夜| 国产视频首页在线观看| 亚洲不卡免费看| 日本免费一区二区三区高清不卡| 免费观看性生交大片5| 99热这里只有是精品50| 最新中文字幕久久久久| 桃色一区二区三区在线观看| 1000部很黄的大片| 日韩欧美精品免费久久| 精品少妇黑人巨大在线播放 | 国产伦在线观看视频一区| 最近最新中文字幕免费大全7| 秋霞在线观看毛片| 寂寞人妻少妇视频99o| 日韩成人av中文字幕在线观看| 国产高清三级在线| 草草在线视频免费看| 狂野欧美白嫩少妇大欣赏| 国语自产精品视频在线第100页| 男女边吃奶边做爰视频| 国产久久久一区二区三区| 国产女主播在线喷水免费视频网站 | 日韩高清综合在线| 欧美激情久久久久久爽电影| 国产成人福利小说| 亚洲人成网站高清观看| 九九久久精品国产亚洲av麻豆| 国产中年淑女户外野战色| 国产精品1区2区在线观看.| 淫秽高清视频在线观看| 国产探花在线观看一区二区| 我的女老师完整版在线观看| 女的被弄到高潮叫床怎么办| 2022亚洲国产成人精品| 成人av在线播放网站| 三级经典国产精品| 欧美性感艳星| 看黄色毛片网站| 久久久久久久久大av| 日本黄色片子视频| 91午夜精品亚洲一区二区三区| 亚洲婷婷狠狠爱综合网| 岛国毛片在线播放| 亚洲av一区综合| 天堂网av新在线| 男女下面进入的视频免费午夜| 国产精品美女特级片免费视频播放器| 亚洲电影在线观看av| 免费观看性生交大片5| 欧美潮喷喷水| 亚洲18禁久久av| 精品99又大又爽又粗少妇毛片| 国产三级中文精品| 成人av在线播放网站| 波野结衣二区三区在线| 69人妻影院| 看片在线看免费视频| 九九久久精品国产亚洲av麻豆| 国产视频内射| 亚洲av中文av极速乱| 日本一二三区视频观看| av天堂中文字幕网| 一个人观看的视频www高清免费观看| 亚洲精华国产精华液的使用体验| 男女边吃奶边做爰视频| 青春草视频在线免费观看| 欧美xxxx性猛交bbbb| 免费av不卡在线播放| 我要搜黄色片| 中国美白少妇内射xxxbb| 国产亚洲91精品色在线| 亚洲人成网站高清观看| 亚洲欧美成人综合另类久久久 | 久久这里只有精品中国| 国产精品久久久久久久久免| 成人毛片a级毛片在线播放| 欧美成人一区二区免费高清观看| 成年版毛片免费区| 久久综合国产亚洲精品| kizo精华| 日韩强制内射视频| 秋霞伦理黄片| 久久精品影院6| 高清毛片免费看| 亚洲乱码一区二区免费版| 好男人在线观看高清免费视频| 女人久久www免费人成看片 | 我的老师免费观看完整版| 国产成人a∨麻豆精品| 日日摸夜夜添夜夜添av毛片| 九九热线精品视视频播放| 天堂网av新在线| 精品一区二区免费观看| 亚洲图色成人| 婷婷六月久久综合丁香| 蜜桃久久精品国产亚洲av| 汤姆久久久久久久影院中文字幕 | 国产一级毛片在线| 自拍偷自拍亚洲精品老妇| 国产人妻一区二区三区在| 在线天堂最新版资源| 久久韩国三级中文字幕| 一个人看视频在线观看www免费| 欧美一区二区国产精品久久精品| 国产欧美日韩精品一区二区| 国产成人福利小说| 一级毛片aaaaaa免费看小| 色吧在线观看| 中文字幕人妻熟人妻熟丝袜美| 中文欧美无线码| 久久精品夜色国产| 日韩一本色道免费dvd| 久久久色成人| 国产免费一级a男人的天堂| 国模一区二区三区四区视频| 色网站视频免费| 久久精品久久精品一区二区三区| 久久精品国产亚洲av涩爱| 国产亚洲最大av| 成人综合一区亚洲| 高清av免费在线| 最近最新中文字幕大全电影3| 亚洲av成人精品一二三区| 亚洲欧美精品综合久久99| 婷婷色麻豆天堂久久 | 欧美日韩精品成人综合77777| 三级国产精品片| 久久热精品热| 免费黄色在线免费观看| 男女啪啪激烈高潮av片| 欧美xxxx性猛交bbbb| 日韩欧美 国产精品| 亚洲四区av| 欧美日韩综合久久久久久| 男女视频在线观看网站免费| 波野结衣二区三区在线| 日韩欧美精品免费久久| 全区人妻精品视频| 亚洲最大成人手机在线| 国产免费一级a男人的天堂| 国产亚洲av嫩草精品影院| 老司机影院毛片| 嫩草影院入口| 国产精品.久久久| 亚洲av熟女| 国产成人一区二区在线| 日本免费a在线| 亚洲无线观看免费| 天美传媒精品一区二区| 天天一区二区日本电影三级| av在线播放精品| 精品少妇黑人巨大在线播放 | 在线天堂最新版资源| 欧美成人一区二区免费高清观看| 国产av码专区亚洲av| 丰满少妇做爰视频| 欧美日韩一区二区视频在线观看视频在线 | 国产精品一区二区性色av| av国产久精品久网站免费入址| 亚洲最大成人中文| 免费观看性生交大片5| 国产精品人妻久久久久久| АⅤ资源中文在线天堂| 免费观看精品视频网站| 中文精品一卡2卡3卡4更新| 亚洲天堂国产精品一区在线| 亚洲av.av天堂| 国产91av在线免费观看| 婷婷色综合大香蕉| 晚上一个人看的免费电影| 国产亚洲5aaaaa淫片| 久99久视频精品免费| 天堂中文最新版在线下载 | 欧美日韩综合久久久久久| 91久久精品国产一区二区成人| 国产精品99久久久久久久久| 观看美女的网站| 欧美高清成人免费视频www| 99在线人妻在线中文字幕| 老师上课跳d突然被开到最大视频| 精品国产三级普通话版| 欧美高清性xxxxhd video| 哪个播放器可以免费观看大片| 青春草国产在线视频| 亚洲精品亚洲一区二区| 一级黄片播放器| 久久久亚洲精品成人影院| 国产片特级美女逼逼视频| 欧美三级亚洲精品| 国产精品嫩草影院av在线观看| 亚洲第一区二区三区不卡| 一级毛片我不卡| 插逼视频在线观看| 春色校园在线视频观看| 天堂中文最新版在线下载 | 国产一级毛片七仙女欲春2| 能在线免费看毛片的网站| 黄色日韩在线| 日本一二三区视频观看| ponron亚洲| 秋霞在线观看毛片| 亚洲精品成人久久久久久| 欧美极品一区二区三区四区| 国产精品一及| 啦啦啦韩国在线观看视频| 亚洲国产欧美人成| 久久久久免费精品人妻一区二区| 亚洲av熟女| 国产精品1区2区在线观看.| 最后的刺客免费高清国语| 久久精品91蜜桃| 亚洲国产色片| 人人妻人人澡欧美一区二区| 午夜精品国产一区二区电影 | 人体艺术视频欧美日本| 中文乱码字字幕精品一区二区三区 | 日韩欧美精品v在线| 国产v大片淫在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 久久久久久久久久久丰满| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 禁无遮挡网站| 天天躁夜夜躁狠狠久久av| 国产乱来视频区| 哪个播放器可以免费观看大片| 精品99又大又爽又粗少妇毛片| 精品一区二区免费观看| 大又大粗又爽又黄少妇毛片口| 色哟哟·www| 97超碰精品成人国产| 一级毛片久久久久久久久女| 国产黄片美女视频| 亚洲av免费高清在线观看| 久久久亚洲精品成人影院| 日本黄色片子视频| 久久婷婷人人爽人人干人人爱| 丝袜喷水一区| 日本猛色少妇xxxxx猛交久久| 国产精品人妻久久久久久| 搡老妇女老女人老熟妇| 人体艺术视频欧美日本| 水蜜桃什么品种好| 国产真实伦视频高清在线观看| 99热这里只有精品一区| 免费不卡的大黄色大毛片视频在线观看 | 亚洲成人精品中文字幕电影| 久久精品国产99精品国产亚洲性色| 国产不卡一卡二| 一卡2卡三卡四卡精品乱码亚洲| 草草在线视频免费看| 亚洲在线自拍视频| 日韩强制内射视频| 欧美激情久久久久久爽电影| 禁无遮挡网站| 亚洲精品久久久久久婷婷小说 | 国产精品国产三级专区第一集| 亚洲精品456在线播放app| 久久婷婷人人爽人人干人人爱| 国产成人精品一,二区| 九九热线精品视视频播放| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品影视一区二区三区av| 欧美三级亚洲精品| 麻豆精品久久久久久蜜桃| 欧美bdsm另类| 国产乱人偷精品视频| 日韩三级伦理在线观看| 91精品一卡2卡3卡4卡| 一二三四中文在线观看免费高清| 久久精品91蜜桃| 少妇人妻精品综合一区二区| 少妇被粗大猛烈的视频| 国产日韩欧美在线精品| 听说在线观看完整版免费高清| 成人高潮视频无遮挡免费网站| videossex国产| 中文天堂在线官网| 蜜桃亚洲精品一区二区三区| 亚洲av.av天堂| 免费一级毛片在线播放高清视频| 日日干狠狠操夜夜爽| 在线天堂最新版资源| 亚洲人与动物交配视频| 26uuu在线亚洲综合色| 又黄又爽又刺激的免费视频.| 亚洲精品国产成人久久av| 亚洲国产精品合色在线| 精品欧美国产一区二区三| 日本与韩国留学比较| 亚洲图色成人| 久久久久久久午夜电影| 色综合色国产| 国产精品久久久久久精品电影小说 | 欧美色视频一区免费| 2021天堂中文幕一二区在线观| 亚洲美女搞黄在线观看| 变态另类丝袜制服| 少妇熟女aⅴ在线视频| 毛片女人毛片| 亚洲自拍偷在线| 老师上课跳d突然被开到最大视频| 国产精品日韩av在线免费观看| 看免费成人av毛片| 中文乱码字字幕精品一区二区三区 | 亚洲av成人精品一二三区| 青春草视频在线免费观看| 一个人免费在线观看电影| 高清午夜精品一区二区三区| 人妻系列 视频| 精品一区二区免费观看| 人人妻人人看人人澡| 国产精品久久久久久av不卡| 国产女主播在线喷水免费视频网站 | 在线播放国产精品三级| 久久精品综合一区二区三区| 三级经典国产精品| 最近2019中文字幕mv第一页| 精品国产三级普通话版| 国产伦精品一区二区三区四那| 狠狠狠狠99中文字幕| 午夜福利在线观看免费完整高清在| 人妻系列 视频| 亚洲综合色惰| 中文资源天堂在线| 性色avwww在线观看| 亚洲精品久久久久久婷婷小说 | 91aial.com中文字幕在线观看| 成人漫画全彩无遮挡| 天天一区二区日本电影三级| 老司机影院成人| 日本色播在线视频| 国产探花在线观看一区二区| 久久精品影院6| 成人高潮视频无遮挡免费网站| 一本一本综合久久| 国产 一区精品| 国产免费又黄又爽又色| 69人妻影院| 成人毛片60女人毛片免费| 免费看av在线观看网站| av又黄又爽大尺度在线免费看 | 一边亲一边摸免费视频| 久久鲁丝午夜福利片| av女优亚洲男人天堂| 欧美日韩精品成人综合77777| 啦啦啦观看免费观看视频高清| 免费观看性生交大片5| 国产免费一级a男人的天堂| 国产精品蜜桃在线观看| 国产亚洲91精品色在线| 精品酒店卫生间| 国产成人aa在线观看| 亚洲美女视频黄频| 国产精品人妻久久久久久| 成人无遮挡网站| 国产极品天堂在线| 国产 一区 欧美 日韩| 成人漫画全彩无遮挡| 欧美三级亚洲精品| 久热久热在线精品观看| 麻豆国产97在线/欧美| 亚洲经典国产精华液单| 在线免费观看不下载黄p国产| 国产成人91sexporn| 色5月婷婷丁香| 国产亚洲精品久久久com| 91av网一区二区| 国产精品电影一区二区三区| av卡一久久| 亚洲真实伦在线观看| 乱码一卡2卡4卡精品| 超碰av人人做人人爽久久| 亚洲一区高清亚洲精品| 亚洲成人av在线免费| 三级国产精品片| 亚洲精品乱久久久久久| 看免费成人av毛片| 热99re8久久精品国产| av免费在线看不卡| 久久精品国产亚洲av涩爱| 亚洲国产色片| 免费播放大片免费观看视频在线观看 | 国产探花极品一区二区| 简卡轻食公司| 国产成人91sexporn| 亚洲性久久影院| 久久久色成人| 亚洲欧美成人综合另类久久久 | 欧美成人一区二区免费高清观看| 欧美日本亚洲视频在线播放| 啦啦啦观看免费观看视频高清| 国产精品久久久久久av不卡| 日本五十路高清| 午夜精品一区二区三区免费看| 久久韩国三级中文字幕| 国产亚洲5aaaaa淫片| 亚洲不卡免费看| av在线天堂中文字幕| 麻豆久久精品国产亚洲av| 久久精品国产亚洲av天美| 深爱激情五月婷婷| 亚洲av日韩在线播放| 蜜桃久久精品国产亚洲av| 成人漫画全彩无遮挡| 人妻夜夜爽99麻豆av| 久久久久久久久久久丰满| 观看美女的网站| 久久热精品热| 真实男女啪啪啪动态图| 亚洲aⅴ乱码一区二区在线播放| 久久欧美精品欧美久久欧美| 国产成人a区在线观看| 一级二级三级毛片免费看| 久久久久久久久久久免费av| 中文乱码字字幕精品一区二区三区 | 在线免费观看的www视频| 成人一区二区视频在线观看|