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

    基于WRF-波浪風(fēng)暴潮耦合模型的深圳市極端風(fēng)暴潮安全分析

    2022-07-01 06:24:44林澤輝
    黑龍江水利科技 2022年5期
    關(guān)鍵詞:風(fēng)暴潮風(fēng)場(chǎng)波浪

    林澤輝

    (惠州市惠陽(yáng)區(qū)永良堤圍管理所,廣東 惠州 516007)

    我國(guó)東南沿海地區(qū)氣象水文條件復(fù)雜,水文方面,海水水位受天體引潮力作用而周期性漲落。氣象方面,西太平洋熱帶地區(qū)氣壓受到擾動(dòng)作用成渦,并不斷發(fā)展移動(dòng),這些氣旋登陸后會(huì)對(duì)人類社會(huì)造成巨大的破壞。另一方面,臺(tái)風(fēng)所造成的大氣壓強(qiáng)場(chǎng)變化也會(huì)導(dǎo)致海面異常起伏,產(chǎn)生風(fēng)暴潮,當(dāng)風(fēng)暴潮恰好疊加上天文大潮,海平面會(huì)巨幅上升,這對(duì)海岸防護(hù)、河口擋潮等水利工程帶來(lái)了巨大挑戰(zhàn)。因此研究臺(tái)風(fēng)過(guò)程前后的波浪以及增水具有重要意義。文章利用…模型對(duì)深圳地區(qū)水域在臺(tái)風(fēng)作用下的水動(dòng)力特征進(jìn)行了數(shù)值模擬研究。

    1 研究方法

    1.1 區(qū)域概況

    深圳市海岸線長(zhǎng)約250km,東起大亞灣,西至珠江口,增水>50cm的風(fēng)暴潮過(guò)程年平均出現(xiàn)2.1次。沿岸分布有機(jī)場(chǎng)、重點(diǎn)港口、大型住宅區(qū)、核電站等大型企業(yè),人口密度極高。

    深圳市年平均風(fēng)速2.6m/s,冬季各月風(fēng)速較大為3.0m/s左右,夏季各月較小約2.0m/s。大風(fēng)日數(shù)(風(fēng)速≥17m/s,即風(fēng)力相當(dāng)于8級(jí)或8級(jí)以上)年平均為7.3d。全年各月均出現(xiàn)大風(fēng),其中以7-9月居多,共占全年大風(fēng)日數(shù)的44%,其次為6月和10月,由于受臺(tái)風(fēng)或熱帶低壓影響,常出現(xiàn)大風(fēng),因此,大風(fēng)出現(xiàn)多的月份,也是臺(tái)風(fēng)季節(jié)。

    1.2 數(shù)據(jù)資料

    影響風(fēng)暴增水量級(jí)的因素有很多,包括臺(tái)風(fēng)氣壓場(chǎng)分布、臺(tái)風(fēng)風(fēng)場(chǎng)的方向與大小以及近岸地形變化因素。在從我國(guó)東南沿海登陸的臺(tái)風(fēng)里,影響到陽(yáng)江至汕尾海域的臺(tái)風(fēng)均可引起深圳海域的風(fēng)暴潮增水。目前在深圳海域,國(guó)家海洋局共設(shè)有兩個(gè)長(zhǎng)期驗(yàn)潮站,分別是赤灣和鹽田站,緊鄰深圳大亞灣岸線的有位于惠州市的惠州海洋站,各站位置見(jiàn)圖1。收集上述站點(diǎn)自建站以來(lái)臺(tái)風(fēng)期間的潮位資料統(tǒng)計(jì)歷史風(fēng)暴潮增水情況。廣東省水文局1964年在深圳市赤灣港內(nèi)建設(shè)有長(zhǎng)期水文觀測(cè)站,1986年國(guó)家海洋局南海分局在相距不到1km處建設(shè)了赤灣海洋站,并于1986年開(kāi)始投入業(yè)務(wù)化觀測(cè)。鹽田海洋站位于大鵬灣鹽田港海域,由國(guó)家海洋局南海分局建設(shè)和管理,2003年投入使用?;葜莺Q笳疚挥诖髞啚碁稠斘鱾?cè)荃灣港區(qū),由國(guó)家海洋局南海分局建設(shè)和管理,2006年建成投入使用,他們的分布位置如圖1所示。

    圖1 資料站位位置示意圖

    1.3 研究方法

    早期人們主要通過(guò)經(jīng)驗(yàn)總結(jié)歸納的方法研究風(fēng)暴潮,觀察臺(tái)風(fēng)天時(shí)海塘外潮位變化研究風(fēng)暴潮的成因。20世紀(jì)中期,隨著計(jì)算機(jī)技術(shù)的發(fā)展,人們開(kāi)始利用數(shù)值模擬研究風(fēng)暴潮。1956年W·Hansen首次利用計(jì)算機(jī)對(duì)歐洲北海的最大風(fēng)暴潮增水進(jìn)行了研究。隨著數(shù)值模擬技術(shù)的進(jìn)步,Jelesnianski建立了研究風(fēng)暴潮的SPLASH模式以及基于二維水動(dòng)力的SLOSH模式,在不斷地摸索研究后,該模式的預(yù)報(bào)精度不斷提高。近年來(lái),不斷發(fā)展的模型計(jì)算實(shí)現(xiàn)了對(duì)復(fù)雜地形影響的考慮,提高了計(jì)算的準(zhǔn)確性[1]。其次,有關(guān)風(fēng)暴增水與天文潮的非線性擬合機(jī)制也變得更加清晰,這也提高了風(fēng)暴潮數(shù)學(xué)模型的準(zhǔn)確性。我國(guó)自20世紀(jì)60年代也開(kāi)展了關(guān)于風(fēng)暴潮的相關(guān)研究。馮士筰所著的《風(fēng)暴潮導(dǎo)論》系統(tǒng)闡述了風(fēng)暴潮的概念,形成機(jī)制,并介紹了早期的風(fēng)暴潮數(shù)值預(yù)報(bào)方法。趙長(zhǎng)進(jìn)等利用ADCIRC潮流與SWAN波浪耦合模式對(duì)長(zhǎng)江口附近海域的風(fēng)暴潮過(guò)程進(jìn)行了模擬。黃本勝等用SELFE模型構(gòu)建了南?!榻涌陔p重嵌套模型,研究了“山竹”臺(tái)風(fēng)過(guò)程中珠江河口的風(fēng)暴增水與波浪增水分布。

    在波浪模擬方面,目前工程界使用較為廣泛的數(shù)學(xué)模型有兩類,一類是波浪流體動(dòng)力學(xué)模型,另一類是動(dòng)譜密度方程。緩坡方程模型與Boussinesq方程模型是主要的波浪流體力學(xué)模型。經(jīng)典的Boussinesq僅僅適用于淺水,并不適用于水深變化大的水域和不規(guī)則波的情況。因此在工程中適用范圍不廣。動(dòng)譜密度方程的基礎(chǔ)是能量平衡理論,在模擬過(guò)程中可提高海岸波浪場(chǎng)的動(dòng)力作用過(guò)程的計(jì)算效率。該方程以動(dòng)譜密度來(lái)描述波浪變形,考慮了對(duì)波-波非線性相互作用、波流相互作用、波浪破碎和底部摩擦引起的能量耗散等過(guò)程,通過(guò)將源項(xiàng)添加到方程中,可以全面考慮動(dòng)力作用的復(fù)雜過(guò)程,使其能完成大尺度區(qū)域的波浪產(chǎn)生、成長(zhǎng)、傳播和變形的計(jì)算過(guò)程。其中對(duì)物理源項(xiàng)的處理和研究是波作用量守恒的關(guān)鍵,波作用量守恒方程發(fā)展到現(xiàn)在已經(jīng)到了第三代,逐漸開(kāi)發(fā)出WAVEWATCH、SWAN、MIKE21SW、TOMAWAC等模型。

    1.4 模型設(shè)置

    1.4.1 模型介紹

    文章采用的是DELFT3D-FLOW與SWAN的耦合模型,DELFT3D是一款由荷蘭DELTRA開(kāi)發(fā)的開(kāi)源水動(dòng)力模擬軟件,在風(fēng)暴潮模擬方面已經(jīng)廣泛應(yīng)用。

    Delft3D-Flow是基于淺水方程和Boussinesq假設(shè),采用有限差分法求解斜壓Navier-Stokes方程的模型,模型控制方程組分別為:

    質(zhì)量守恒方程:

    (1)

    式中:ζ為水位(總水深);d為到參考平面的當(dāng)前水深(凈水深);U、V-x、y為方向上的速度分量;Q為單位面積上質(zhì)量強(qiáng)度。

    動(dòng)量守恒方程:

    (2)

    (3)

    式中:f為科氏力參數(shù);g為重力加速度;νh為動(dòng)態(tài)水平渦流黏度;τsx,τsy為x、y方向上的海平面風(fēng)壓分量;τbx,τby為x、y方向上的底部的剪應(yīng)力分量;ρ0為參考密度;ρ'為不規(guī)則密度。

    SWAN模型的控制方程為:

    (4)

    式中:Cx,Cy,Cσ,Cθ為二維動(dòng)譜密度在x、y、σ、θ為空間上的傳播速度;σ為相對(duì)頻率;θ為波向;S為能量源項(xiàng)。

    1.4.2 模型設(shè)置

    1)模型范圍:

    文章所采用的模型囊括整個(gè)珠江三角洲區(qū)域,見(jiàn)圖2。

    圖2 模型區(qū)域示意圖

    2)地形水深:

    模擬范圍內(nèi)的地形水深采用的是中華人民共和國(guó)海事局2016年繪制的海圖插值數(shù)據(jù)。將地形水深數(shù)據(jù)轉(zhuǎn)化到統(tǒng)一基面(平均海平面)后,用IDW插值方法插值到網(wǎng)格點(diǎn)上,得到模型的水深地形。

    3)開(kāi)邊界設(shè)置:

    開(kāi)邊界采用天文潮水位,主要考慮M2、N2、S2、K2、K1、O1、P1、Q1等八個(gè)天文分潮以及MS4、M4、M6等淺水分潮共11個(gè)分潮的作用。分潮的調(diào)和分析常數(shù)由2016-2019年的潮位實(shí)測(cè)數(shù)據(jù)經(jīng)MATLAB中的T_TIDE工具包調(diào)和分析所得。模型上部采用流量邊界,輸入數(shù)據(jù)為梧州、石角、博羅三個(gè)水文站的2016-2019洪季平均流量。

    4)底摩擦設(shè)置:

    模型根據(jù)水深范圍設(shè)置不同大小的糙率,水深越大糙率越小。

    1.4.3 風(fēng)場(chǎng)輸入

    文章的風(fēng)場(chǎng)采用美國(guó)國(guó)家大氣研究中心(NCAR)等機(jī)構(gòu)共同開(kāi)發(fā)的中尺度大氣預(yù)報(bào)系統(tǒng)。WRF模式適用于臺(tái)風(fēng)模擬,可以通過(guò)嵌套網(wǎng)格運(yùn)算來(lái)達(dá)到高精度模擬的目的。本次采用了WSM6微物理方案、KainFritsch卷積云方案、Monin-Obukhov表面層方案、Noah陸面層方案、MYJ大氣邊界層方案、RRTMG輻射方案模擬臺(tái)風(fēng)。

    2 模型驗(yàn)證

    2.1 風(fēng)場(chǎng)驗(yàn)證

    首先利用實(shí)測(cè)資料對(duì)所建立的風(fēng)暴潮-波浪耦合模型進(jìn)行率定驗(yàn)證。選取1713號(hào)臺(tái)風(fēng)(2017年第13號(hào)臺(tái)風(fēng)—天鴿,其移行路徑如圖3所示)作為模型邊界條件,為了證明WRF模式的準(zhǔn)確性其在風(fēng)暴潮模擬方面的優(yōu)勢(shì),文章也用其他論文常用的經(jīng)驗(yàn)風(fēng)場(chǎng)模型做了驗(yàn)證對(duì)比。經(jīng)驗(yàn)風(fēng)場(chǎng)采用的是宮崎正衛(wèi)模型作為移動(dòng)風(fēng)場(chǎng),HOLLAND模型作為梯度風(fēng)場(chǎng),他們的公式分別為:

    圖3 “天鴿”臺(tái)風(fēng)移行路徑

    (5)

    (5)

    式中:B為HOLLAND模型氣壓參數(shù);f為科氏力參數(shù),與模擬區(qū)域緯度有關(guān);R為最大風(fēng)速半徑,與中心氣壓有關(guān)。

    圖4為WRF模型模擬風(fēng)場(chǎng)值以及經(jīng)驗(yàn)風(fēng)場(chǎng)模擬值與實(shí)測(cè)風(fēng)場(chǎng)的對(duì)比,可以看出,WRF模型模擬的風(fēng)場(chǎng)在峰值以及過(guò)程上都比經(jīng)驗(yàn)風(fēng)場(chǎng)模型更貼近實(shí)際,因此利用WRF模型生成的臺(tái)風(fēng)風(fēng)場(chǎng)作為風(fēng)暴潮-波浪耦合模型的邊界條件更為合理。

    圖4 WRF風(fēng)場(chǎng)及模型風(fēng)場(chǎng)與實(shí)測(cè)風(fēng)場(chǎng)對(duì)比

    2.2 風(fēng)暴潮及波浪率定驗(yàn)證

    以1713號(hào)臺(tái)風(fēng)(2017年第13號(hào)臺(tái)風(fēng))為邊界條件分別驅(qū)動(dòng)珠江口風(fēng)暴潮模型以及珠江口波浪-風(fēng)暴潮耦合模型,得到如圖5所示的赤灣站與鹽田站水位驗(yàn)證曲線,可以看出,考慮波浪后的水位過(guò)程曲線更符合實(shí)際,最大誤差絕對(duì)值為13cm,平均相對(duì)誤差為11.3%,標(biāo)準(zhǔn)差為2.94,符合模擬精度要求。

    圖5(a) 赤灣站不同模型下水位曲線

    圖5(b) 鹽田站不同模型下水位曲線

    3 深圳市沿岸漫堤風(fēng)險(xiǎn)評(píng)估

    3.1 可能最大風(fēng)暴潮模型計(jì)算

    為了尋求可能影響深圳市的最大風(fēng)暴潮對(duì)深圳市沿岸的影響,采用概率論法進(jìn)行計(jì)算可能最大臺(tái)風(fēng)中的近中心最大風(fēng)速Vmax,臺(tái)風(fēng)中心氣壓P0。

    利用中國(guó)氣象局(CMA)1949-2015年熱帶氣旋最佳路徑集(Best-track)數(shù)據(jù),選取出研究區(qū)域400km直徑范圍內(nèi)對(duì)本研究區(qū)域有重要影響的臺(tái)風(fēng),以所選取臺(tái)風(fēng)的最小臺(tái)風(fēng)中心氣壓值作樣本。采用極值I型分布計(jì)算1000年一遇的臺(tái)風(fēng)中心氣壓值作為可能最大臺(tái)風(fēng)的中心氣壓,計(jì)算結(jié)果如圖6(a)所示,即深圳區(qū)域1000a遇的臺(tái)風(fēng)中心氣壓P0=880hPa。

    圖6(a) 影響深圳的臺(tái)風(fēng)中心最低氣壓極值I型分布

    圖6(b) 影響深圳的臺(tái)風(fēng)中心最大風(fēng)速極值I型分布

    利用中國(guó)氣象局(CMA)1949-2015年熱帶氣旋最佳路徑集(Best-track)數(shù)據(jù),取深圳市海灣區(qū)域400km直徑范圍內(nèi)歷年路經(jīng)本海區(qū)的臺(tái)風(fēng)最大風(fēng)速值作樣本,對(duì)于沒(méi)有臺(tái)風(fēng)經(jīng)過(guò)研究區(qū)域的年份,該年Vmax取為臺(tái)風(fēng)Vmax系列中的平均值。采用極值I型分布計(jì)算1000a一遇的最大風(fēng)速值作為可能最大臺(tái)風(fēng)的最大風(fēng)速值,計(jì)算結(jié)果如圖6(b)所示,深圳區(qū)域1000a一遇的臺(tái)風(fēng)中心最大風(fēng)速Vmax=85m/s。另外,臺(tái)風(fēng)路徑選取為東南向西北由大鵬灣登陸的路徑。

    輸出堤前代表點(diǎn)(具體見(jiàn)圖7)的最高潮位與堤頂高程進(jìn)行比對(duì),了解沿岸淹沒(méi)的情況,由于模式計(jì)算結(jié)果是相對(duì)平均海平面的,為便于比對(duì),需要根據(jù)當(dāng)?shù)仄骄F矫媾c85高程的關(guān)系,將計(jì)算結(jié)果轉(zhuǎn)成以85黃海基面為起算面,其中西段提防采用赤灣站高程關(guān)系,即計(jì)算結(jié)果加上0.54m;東段提防采用鹽田站高程關(guān)系進(jìn)行轉(zhuǎn)換,即計(jì)算結(jié)果加上0.66m。表1列出堤前統(tǒng)計(jì)的最高潮位及堤頂高程。從代表點(diǎn)的最高潮位和海堤高程比較來(lái)看,在可能最大熱帶氣旋影響下,深圳沿岸均可能發(fā)生漫堤。

    圖7 堤前代表點(diǎn)位置示意圖

    表1 可能最大風(fēng)暴潮漫堤分析 (單位:m,基面:85高程)

    續(xù)表1 可能最大風(fēng)暴潮漫堤分析 (單位:m,基面:85高程)

    3.2 海堤防御標(biāo)準(zhǔn)

    對(duì)于大部分沿海地區(qū)來(lái)說(shuō),1000a一遇的標(biāo)準(zhǔn)是過(guò)高的。但是深圳作為特大型城市,其風(fēng)暴潮安全具有重要意義,因此考慮1000a一遇情形下的風(fēng)暴潮安全是有必要的,在此情形下得出的淹沒(méi)情況有利于應(yīng)急管理部門采取相應(yīng)的防范措施。

    4 結(jié) 論

    1)利用WRF中長(zhǎng)尺度氣象預(yù)報(bào)模式作為波浪-風(fēng)暴潮耦合模型的風(fēng)場(chǎng)邊界條件輸入,提高了風(fēng)暴潮模型的精確度

    2)在珠江口波浪-風(fēng)暴潮模型驗(yàn)證實(shí)測(cè)資料準(zhǔn)確的基礎(chǔ)上,利用極值I型分布確定1000a一遇臺(tái)風(fēng)參數(shù)對(duì)影響深圳市海岸堤防未來(lái)發(fā)生的最大可能風(fēng)暴潮進(jìn)行計(jì)算分析。文章研究成果可為深圳市風(fēng)暴潮安全防范提供重要技術(shù)支撐,同時(shí)也為海洋水利防災(zāi)應(yīng)急部門提供決策參考。

    猜你喜歡
    風(fēng)暴潮風(fēng)場(chǎng)波浪
    波浪谷和波浪巖
    基于FLUENT的下?lián)舯┝魅S風(fēng)場(chǎng)建模
    2012年“蘇拉”和“達(dá)維”雙臺(tái)風(fēng)影響的近海風(fēng)暴潮過(guò)程
    防范未來(lái)風(fēng)暴潮災(zāi)害的綠色海堤藍(lán)圖
    科學(xué)(2020年4期)2020-11-26 08:27:00
    基于多變量LSTM神經(jīng)網(wǎng)絡(luò)模型的風(fēng)暴潮臨近預(yù)報(bào)
    波浪谷隨想
    去看神奇波浪谷
    “最美風(fēng)場(chǎng)”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    側(cè)向風(fēng)場(chǎng)中無(wú)人機(jī)的飛行研究
    廣東省風(fēng)暴潮時(shí)空分布特征及重點(diǎn)城市風(fēng)暴潮風(fēng)險(xiǎn)研究
    日本色播在线视频| 亚洲四区av| 建设人人有责人人尽责人人享有的| 国产免费又黄又爽又色| 国产av精品麻豆| 免费在线观看黄色视频的| 亚洲欧洲精品一区二区精品久久久 | a级毛色黄片| 亚洲国产看品久久| 永久网站在线| 美女 人体艺术 gogo| 成人影院久久| 建设人人有责人人尽责人人享有的| 国产蜜桃级精品一区二区三区 | 中文亚洲av片在线观看爽 | 看免费av毛片| 亚洲精品在线美女| 久久久久精品人妻al黑| 欧美日韩一级在线毛片| av国产精品久久久久影院| 国产精品99久久99久久久不卡| 亚洲欧美色中文字幕在线| 一个人免费在线观看的高清视频| 在线观看日韩欧美| 一级毛片高清免费大全| 亚洲伊人色综图| 日韩熟女老妇一区二区性免费视频| 亚洲熟女精品中文字幕| 国产激情欧美一区二区| 伊人久久大香线蕉亚洲五| 成年人免费黄色播放视频| 亚洲av欧美aⅴ国产| 中出人妻视频一区二区| 岛国在线观看网站| 精品国产一区二区三区四区第35| 久久99一区二区三区| 亚洲一区高清亚洲精品| 欧美精品啪啪一区二区三区| 免费少妇av软件| 男女之事视频高清在线观看| 精品午夜福利视频在线观看一区| 老司机在亚洲福利影院| 宅男免费午夜| 亚洲一卡2卡3卡4卡5卡精品中文| 国产激情欧美一区二区| 极品少妇高潮喷水抽搐| 国产男女超爽视频在线观看| 首页视频小说图片口味搜索| 老司机亚洲免费影院| tube8黄色片| 99国产综合亚洲精品| 日本一区二区免费在线视频| 午夜91福利影院| 日本五十路高清| 成人免费观看视频高清| 国产成人一区二区三区免费视频网站| av电影中文网址| 极品少妇高潮喷水抽搐| 交换朋友夫妻互换小说| 搡老乐熟女国产| 黑人巨大精品欧美一区二区mp4| 80岁老熟妇乱子伦牲交| 国产精品1区2区在线观看. | 欧美黄色淫秽网站| 欧美激情久久久久久爽电影 | 国产亚洲精品第一综合不卡| 18禁美女被吸乳视频| 久久久久视频综合| 无限看片的www在线观看| videos熟女内射| 俄罗斯特黄特色一大片| 亚洲第一欧美日韩一区二区三区| 99国产精品一区二区三区| 国产xxxxx性猛交| 大型av网站在线播放| 91麻豆av在线| 亚洲欧美日韩另类电影网站| 亚洲精品久久成人aⅴ小说| 交换朋友夫妻互换小说| 免费久久久久久久精品成人欧美视频| 国产视频一区二区在线看| 18禁国产床啪视频网站| 无人区码免费观看不卡| 99国产精品一区二区三区| 啦啦啦免费观看视频1| 国产精品久久久久久人妻精品电影| 亚洲精品国产一区二区精华液| 曰老女人黄片| 天堂动漫精品| 99精品久久久久人妻精品| 国产激情欧美一区二区| 亚洲欧美色中文字幕在线| 国产成人精品无人区| 美女扒开内裤让男人捅视频| 热re99久久精品国产66热6| 男女床上黄色一级片免费看| 国产成人免费观看mmmm| 亚洲欧美精品综合一区二区三区| 在线天堂中文资源库| 国产极品粉嫩免费观看在线| 一本大道久久a久久精品| 亚洲国产精品sss在线观看 | 69精品国产乱码久久久| 亚洲欧美一区二区三区久久| 一区二区日韩欧美中文字幕| 黄片大片在线免费观看| 亚洲成人手机| 丝瓜视频免费看黄片| 中文字幕人妻丝袜一区二区| 黄色视频,在线免费观看| 天堂中文最新版在线下载| 国产极品粉嫩免费观看在线| 免费在线观看视频国产中文字幕亚洲| 国产又色又爽无遮挡免费看| 如日韩欧美国产精品一区二区三区| 亚洲精品久久午夜乱码| 涩涩av久久男人的天堂| 国产成人精品久久二区二区91| av一本久久久久| 丰满饥渴人妻一区二区三| 亚洲熟妇熟女久久| 成年女人毛片免费观看观看9 | 国产亚洲精品久久久久5区| 欧美日韩亚洲国产一区二区在线观看 | 天天躁狠狠躁夜夜躁狠狠躁| 国产黄色免费在线视频| 一级毛片高清免费大全| 日韩欧美一区视频在线观看| 99精品久久久久人妻精品| 久久精品国产综合久久久| 亚洲欧美色中文字幕在线| av超薄肉色丝袜交足视频| xxx96com| 99国产极品粉嫩在线观看| 亚洲欧美一区二区三区黑人| av天堂久久9| 纯流量卡能插随身wifi吗| 天天添夜夜摸| 在线播放国产精品三级| 老熟女久久久| 在线观看午夜福利视频| 欧美乱码精品一区二区三区| 欧美午夜高清在线| 很黄的视频免费| a级毛片黄视频| 亚洲,欧美精品.| 日本一区二区免费在线视频| 97人妻天天添夜夜摸| 高清视频免费观看一区二区| 女人精品久久久久毛片| 精品国产亚洲在线| 亚洲欧美一区二区三区黑人| 亚洲精品久久午夜乱码| av在线播放免费不卡| 国产精华一区二区三区| 一本综合久久免费| 日韩欧美三级三区| 日本vs欧美在线观看视频| 欧美精品啪啪一区二区三区| 久久香蕉国产精品| 国产精品欧美亚洲77777| 999精品在线视频| 久久国产精品大桥未久av| 一级片'在线观看视频| 亚洲成人免费电影在线观看| 性色av乱码一区二区三区2| av一本久久久久| 国产欧美日韩一区二区三区在线| 亚洲av欧美aⅴ国产| 高潮久久久久久久久久久不卡| 免费看a级黄色片| 国产精品久久久久久人妻精品电影| 国产成人av激情在线播放| av天堂在线播放| 亚洲五月天丁香| 精品少妇久久久久久888优播| 精品亚洲成a人片在线观看| 夜夜爽天天搞| 在线国产一区二区在线| 久久精品国产清高在天天线| 一进一出抽搐动态| 又紧又爽又黄一区二区| 国产有黄有色有爽视频| 制服人妻中文乱码| 亚洲国产毛片av蜜桃av| 久99久视频精品免费| 久久精品亚洲熟妇少妇任你| 欧美激情极品国产一区二区三区| 久久久久精品国产欧美久久久| 好男人电影高清在线观看| 国产成人精品久久二区二区91| 我的亚洲天堂| 久久精品国产亚洲av香蕉五月 | av欧美777| 69精品国产乱码久久久| 亚洲欧美日韩另类电影网站| 午夜亚洲福利在线播放| 日日爽夜夜爽网站| 18禁国产床啪视频网站| 久久久久久亚洲精品国产蜜桃av| 久久香蕉精品热| 黑丝袜美女国产一区| 日韩欧美在线二视频 | 精品高清国产在线一区| 国产亚洲av高清不卡| 黄片小视频在线播放| 操美女的视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲欧美一区二区三区黑人| 在线国产一区二区在线| 1024香蕉在线观看| 欧美日韩亚洲国产一区二区在线观看 | 男女床上黄色一级片免费看| 午夜成年电影在线免费观看| 精品亚洲成国产av| 亚洲一卡2卡3卡4卡5卡精品中文| 国产亚洲av高清不卡| 大陆偷拍与自拍| 欧美精品av麻豆av| 精品一区二区三卡| 黄色片一级片一级黄色片| 一边摸一边抽搐一进一出视频| 天堂√8在线中文| 欧美日韩亚洲国产一区二区在线观看 | 国产欧美亚洲国产| 最新美女视频免费是黄的| 色在线成人网| 女警被强在线播放| av免费在线观看网站| videosex国产| 高清在线国产一区| 国产亚洲精品第一综合不卡| av福利片在线| 久久久久久免费高清国产稀缺| 免费观看人在逋| 黄片播放在线免费| 国产精品国产av在线观看| 两人在一起打扑克的视频| 在线观看www视频免费| 黄频高清免费视频| netflix在线观看网站| 亚洲成人国产一区在线观看| 熟女少妇亚洲综合色aaa.| 嫩草影视91久久| 亚洲人成电影观看| 曰老女人黄片| 亚洲av第一区精品v没综合| 精品少妇久久久久久888优播| avwww免费| 后天国语完整版免费观看| 色婷婷久久久亚洲欧美| 午夜免费观看网址| 国产黄色免费在线视频| 18禁裸乳无遮挡免费网站照片 | 日韩熟女老妇一区二区性免费视频| 在线av久久热| 亚洲av片天天在线观看| 免费在线观看影片大全网站| 国产精品九九99| 亚洲第一欧美日韩一区二区三区| 亚洲中文字幕日韩| 老司机影院毛片| √禁漫天堂资源中文www| 老司机靠b影院| 日本a在线网址| 制服人妻中文乱码| 老熟女久久久| 久久人妻熟女aⅴ| 国产97色在线日韩免费| 免费日韩欧美在线观看| 精品久久久久久久毛片微露脸| xxx96com| 在线免费观看的www视频| 成人三级做爰电影| 日本欧美视频一区| 亚洲成国产人片在线观看| 一进一出抽搐gif免费好疼 | 久久久久久久久免费视频了| 日韩欧美一区二区三区在线观看 | 亚洲专区中文字幕在线| 欧美日韩一级在线毛片| 每晚都被弄得嗷嗷叫到高潮| 后天国语完整版免费观看| 变态另类成人亚洲欧美熟女 | 亚洲av电影在线进入| 我的亚洲天堂| 另类亚洲欧美激情| svipshipincom国产片| 国产精品久久久av美女十八| 青草久久国产| 国产日韩一区二区三区精品不卡| 国产欧美日韩一区二区三区在线| 亚洲精华国产精华精| 女人爽到高潮嗷嗷叫在线视频| 黄色视频不卡| 国产高清国产精品国产三级| 精品国产美女av久久久久小说| 精品午夜福利视频在线观看一区| 国产在线观看jvid| 十分钟在线观看高清视频www| 国产欧美日韩一区二区三区在线| 中文欧美无线码| 欧美日韩精品网址| 69av精品久久久久久| 黑人操中国人逼视频| 欧美激情极品国产一区二区三区| 亚洲 欧美一区二区三区| 热99久久久久精品小说推荐| 国产精品欧美亚洲77777| 天天影视国产精品| 一区二区三区精品91| 欧美日韩国产mv在线观看视频| 宅男免费午夜| 最近最新中文字幕大全电影3 | 久久久久国产一级毛片高清牌| 婷婷精品国产亚洲av在线 | 国产一区在线观看成人免费| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品国产a三级三级三级| 制服诱惑二区| 老司机靠b影院| 国产精品成人在线| 日韩中文字幕欧美一区二区| 国产成人影院久久av| a级毛片黄视频| 天堂中文最新版在线下载| 91在线观看av| 国产有黄有色有爽视频| 又黄又粗又硬又大视频| а√天堂www在线а√下载 | 波多野结衣一区麻豆| 日本黄色视频三级网站网址 | 欧美+亚洲+日韩+国产| 国产欧美日韩一区二区精品| 女人被躁到高潮嗷嗷叫费观| 久久久国产成人免费| 999久久久精品免费观看国产| 日韩大码丰满熟妇| 亚洲精品在线美女| 欧美色视频一区免费| 999久久久精品免费观看国产| a在线观看视频网站| 亚洲欧美日韩另类电影网站| 超色免费av| 女性被躁到高潮视频| 欧美黑人欧美精品刺激| 亚洲人成电影观看| 国产精华一区二区三区| 午夜91福利影院| 国产欧美日韩精品亚洲av| 大码成人一级视频| 欧美日韩黄片免| а√天堂www在线а√下载 | 久久精品aⅴ一区二区三区四区| 欧美成狂野欧美在线观看| 亚洲av成人一区二区三| 成人免费观看视频高清| 国产1区2区3区精品| 女人爽到高潮嗷嗷叫在线视频| 国产精品.久久久| 免费一级毛片在线播放高清视频 | 国产又色又爽无遮挡免费看| 法律面前人人平等表现在哪些方面| 久久精品熟女亚洲av麻豆精品| 高清av免费在线| 巨乳人妻的诱惑在线观看| 久久精品国产99精品国产亚洲性色 | 大码成人一级视频| 久久精品熟女亚洲av麻豆精品| 欧美另类亚洲清纯唯美| 丰满的人妻完整版| 精品久久久久久久久久免费视频 | 精品国产乱码久久久久久男人| 天天躁日日躁夜夜躁夜夜| 国产精品自产拍在线观看55亚洲 | 日韩人妻精品一区2区三区| 色老头精品视频在线观看| 亚洲av熟女| 色婷婷久久久亚洲欧美| 中亚洲国语对白在线视频| 久久国产精品影院| 桃红色精品国产亚洲av| 美女视频免费永久观看网站| 在线观看免费午夜福利视频| 十八禁高潮呻吟视频| 精品久久久久久电影网| 亚洲少妇的诱惑av| 777米奇影视久久| 久久国产精品男人的天堂亚洲| 精品午夜福利视频在线观看一区| 国产淫语在线视频| 国产aⅴ精品一区二区三区波| 久久天堂一区二区三区四区| 黄色怎么调成土黄色| 精品人妻熟女毛片av久久网站| 亚洲 欧美一区二区三区| 国产亚洲欧美在线一区二区| 欧美乱码精品一区二区三区| 男女下面插进去视频免费观看| 国产精品一区二区精品视频观看| 久久精品亚洲av国产电影网| 精品电影一区二区在线| 欧美黄色淫秽网站| 久久亚洲精品不卡| 精品一区二区三区av网在线观看| 五月开心婷婷网| 在线观看舔阴道视频| 狠狠婷婷综合久久久久久88av| www日本在线高清视频| 交换朋友夫妻互换小说| 成人特级黄色片久久久久久久| tocl精华| 嫩草影视91久久| e午夜精品久久久久久久| 欧美av亚洲av综合av国产av| 欧美国产精品一级二级三级| 在线免费观看的www视频| 50天的宝宝边吃奶边哭怎么回事| 91成年电影在线观看| netflix在线观看网站| 国产精品av久久久久免费| 高清欧美精品videossex| 69av精品久久久久久| 亚洲欧美日韩另类电影网站| 欧美成狂野欧美在线观看| 国产欧美日韩一区二区三区在线| 十分钟在线观看高清视频www| 99re在线观看精品视频| 久久影院123| 男女床上黄色一级片免费看| 国产成人精品久久二区二区免费| 亚洲精品中文字幕一二三四区| av免费在线观看网站| 久久性视频一级片| 18禁观看日本| 大码成人一级视频| 18禁国产床啪视频网站| 国产1区2区3区精品| 叶爱在线成人免费视频播放| 777久久人妻少妇嫩草av网站| 国产成人影院久久av| 免费在线观看影片大全网站| 最新在线观看一区二区三区| 18禁观看日本| tocl精华| 午夜影院日韩av| 国产精品免费大片| 天堂动漫精品| 一本综合久久免费| 亚洲精品自拍成人| 三上悠亚av全集在线观看| 亚洲欧美激情综合另类| 一区福利在线观看| 亚洲在线自拍视频| 精品久久久久久久毛片微露脸| 久久精品亚洲熟妇少妇任你| 超色免费av| 999久久久国产精品视频| 成人国语在线视频| 两个人免费观看高清视频| 丰满迷人的少妇在线观看| 大香蕉久久网| 欧美国产精品va在线观看不卡| 一边摸一边抽搐一进一小说 | 精品久久蜜臀av无| 少妇被粗大的猛进出69影院| 免费在线观看日本一区| 免费观看人在逋| 侵犯人妻中文字幕一二三四区| 欧美人与性动交α欧美精品济南到| 亚洲少妇的诱惑av| 国内久久婷婷六月综合欲色啪| 亚洲久久久国产精品| 91在线观看av| 久久精品91无色码中文字幕| 欧美 日韩 精品 国产| 色婷婷av一区二区三区视频| 国产蜜桃级精品一区二区三区 | 久久精品国产亚洲av高清一级| 一二三四在线观看免费中文在| 欧美人与性动交α欧美精品济南到| 国产亚洲精品第一综合不卡| 欧美精品av麻豆av| 露出奶头的视频| 女性被躁到高潮视频| 91国产中文字幕| 精品人妻1区二区| 夫妻午夜视频| 女人久久www免费人成看片| 欧美精品啪啪一区二区三区| 如日韩欧美国产精品一区二区三区| 日本五十路高清| 色尼玛亚洲综合影院| 国产一区二区激情短视频| 欧美 亚洲 国产 日韩一| 亚洲人成电影观看| av中文乱码字幕在线| 国内久久婷婷六月综合欲色啪| 岛国毛片在线播放| 亚洲一区中文字幕在线| 亚洲一卡2卡3卡4卡5卡精品中文| 天天影视国产精品| 午夜免费鲁丝| 亚洲精品中文字幕一二三四区| 亚洲精品在线美女| 久久精品国产99精品国产亚洲性色 | 久久 成人 亚洲| 亚洲av熟女| 免费在线观看影片大全网站| 婷婷丁香在线五月| 亚洲精品美女久久av网站| 99国产精品一区二区三区| 久久性视频一级片| 国产一区有黄有色的免费视频| 欧美久久黑人一区二区| 中文字幕另类日韩欧美亚洲嫩草| 精品国产超薄肉色丝袜足j| av片东京热男人的天堂| 国产精品免费一区二区三区在线 | 亚洲中文av在线| 曰老女人黄片| av福利片在线| 亚洲精品久久成人aⅴ小说| 日日爽夜夜爽网站| 一级片免费观看大全| 自线自在国产av| 国产成人啪精品午夜网站| 自线自在国产av| 国产成人啪精品午夜网站| 女性生殖器流出的白浆| 老熟妇乱子伦视频在线观看| 国产国语露脸激情在线看| 两个人免费观看高清视频| 亚洲中文日韩欧美视频| 日韩人妻精品一区2区三区| 欧美+亚洲+日韩+国产| 成人18禁在线播放| 免费高清在线观看日韩| 亚洲熟女毛片儿| 十分钟在线观看高清视频www| 久久中文看片网| 午夜福利影视在线免费观看| 嫁个100分男人电影在线观看| 十八禁高潮呻吟视频| 国产日韩一区二区三区精品不卡| 免费观看人在逋| 亚洲精品国产区一区二| 久久久国产一区二区| 国产一区二区激情短视频| 欧美精品亚洲一区二区| 热99re8久久精品国产| 人妻 亚洲 视频| 少妇猛男粗大的猛烈进出视频| 一级,二级,三级黄色视频| 国产深夜福利视频在线观看| 欧美大码av| 欧美日韩乱码在线| 飞空精品影院首页| 夫妻午夜视频| 极品教师在线免费播放| 不卡一级毛片| 午夜久久久在线观看| 国产av精品麻豆| 一区二区三区激情视频| 中出人妻视频一区二区| 国产一区二区激情短视频| 一二三四在线观看免费中文在| 性少妇av在线| 中文字幕人妻丝袜一区二区| 国产精品国产高清国产av | 午夜福利乱码中文字幕| 国产1区2区3区精品| 久久久久国内视频| 午夜影院日韩av| 怎么达到女性高潮| av超薄肉色丝袜交足视频| 女人高潮潮喷娇喘18禁视频| 我的亚洲天堂| 久久亚洲真实| videos熟女内射| 亚洲精品久久午夜乱码| 日韩人妻精品一区2区三区| 亚洲少妇的诱惑av| 黄色怎么调成土黄色| 久久亚洲真实| videos熟女内射| 免费观看精品视频网站| 老鸭窝网址在线观看| 涩涩av久久男人的天堂| 精品国内亚洲2022精品成人 | 免费在线观看视频国产中文字幕亚洲| 99香蕉大伊视频| 另类亚洲欧美激情| 老司机深夜福利视频在线观看| 操出白浆在线播放| 宅男免费午夜| 亚洲一区中文字幕在线| 性少妇av在线| 久久久久精品人妻al黑| 变态另类成人亚洲欧美熟女 | 日本撒尿小便嘘嘘汇集6| 中文字幕人妻丝袜一区二区| 午夜激情av网站| 亚洲一区高清亚洲精品| 大香蕉久久网| 中文亚洲av片在线观看爽 | 捣出白浆h1v1| 又紧又爽又黄一区二区| 大码成人一级视频| 国产精品一区二区在线观看99| 一区二区三区激情视频| 国产男女内射视频| 天天添夜夜摸| 午夜视频精品福利| 国产片内射在线|