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

    海堤龍口水力計(jì)算方法探討

    2019-09-17 06:21:16
    長江科學(xué)院院報(bào) 2019年9期
    關(guān)鍵詞:海堤龍口潮位

    (浙江省水利水電勘測設(shè)計(jì)院,杭州 310002)

    1 研究背景

    浙江沿海地區(qū)在灘涂保護(hù)與利用開發(fā)中,海堤龍口段的地基大多數(shù)為深厚淤泥軟土地基,受進(jìn)排水要求限制,海堤龍口段拋填石料厚度一般較薄,其預(yù)壓固結(jié)程度不如一般海堤段的預(yù)壓固結(jié)程度,其抗沖刷能力相對較差。龍口寬度太窄則龍口度汛期漲落潮時水流相對較急,容易造成龍口嚴(yán)重沖刷,影響工程安全穩(wěn)定性;龍口寬度太寬則造成堵口難度大,影響工程總體工期,且不易合龍。因此,選擇合適的龍口寬度在深厚淤泥地基海堤工程中顯得尤為重要。

    吳繼偉[1]采用水量平衡法編制了海堤龍口水力計(jì)算軟件,但其未考慮龍口在堵口過程中底坎、口門寬度的動態(tài)變化過程,且對堵口過程中拋石底坎的透水性以及水閘的影響等要素未深入考慮,軟件計(jì)算模型仍有待深入細(xì)化。

    趙庚潤等[2]利用流體有限元軟件Fluent中的k-ε湍流模型對海堤龍口二維流動進(jìn)行了數(shù)值模擬計(jì)算,認(rèn)為落潮時龍口最大流速可達(dá)8.70 m/s,漲潮時龍口最大流速可達(dá)10.60 m/s;其中大流速區(qū)一般位于堤頂或背流側(cè)的區(qū)域;甚至龍口附近會出現(xiàn)水躍或者卷浪型流態(tài)。

    胡志根等[3]采用物理模型試驗(yàn)對大比降、非對稱基床中的合龍進(jìn)行研究,觀測分析堵口合龍進(jìn)占方式、龍口位置、龍口水流流態(tài)及其龍口水力參數(shù)之間的關(guān)系。

    楊磊等[4]通過相關(guān)模型試驗(yàn)對寬戧堤龍口堵口范圍內(nèi)的水流特性進(jìn)行試驗(yàn)研究,認(rèn)為當(dāng)寬戧堤堤頂寬度B與總水頭H0的比值B/H0在311~517之間時,龍口內(nèi)水面形態(tài)將逐漸由單降型流態(tài)漸變?yōu)殡p降型流態(tài);但當(dāng)B/H0>517時,龍口戧堤堤頂?shù)膶挾刃?yīng)逐漸減弱。

    工程實(shí)際海堤龍口計(jì)算中外海潮位等水力邊界條件不斷變化,導(dǎo)致龍口水力計(jì)算較復(fù)雜,一般需采用簡化小程序計(jì)算(如文獻(xiàn)[1])?;诖耍疚囊罁?jù)水力學(xué)理論,結(jié)合MATLAB軟件,采用黃朝煊等[5]研究中給出Runge-Kutta法對海堤龍口水力計(jì)算進(jìn)行分析研究,以便工程應(yīng)用推廣。

    2 龍口水力計(jì)算基本方程

    根據(jù)文獻(xiàn)[6]和文獻(xiàn)[7]可知,由圍區(qū)內(nèi)進(jìn)(出)水量平衡原理得龍口水力計(jì)算基本方程為

    (1)

    q外=q閘+q泄+q滲,

    (2)

    (3)

    (4)

    式中:q內(nèi)為內(nèi)陸流域來水平均流量(m3/s);q外為受外海潮位影響流進(jìn)(出)圍區(qū)的流量(m3/s);q閘為水閘過流流量(m3/s);q泄為口門堰上溢流流量(m3/s);q滲為計(jì)算時段內(nèi)海堤龍口拋石體滲流平均流量(m3/s);V[Z(t)]為t時刻圍區(qū)內(nèi)庫容量(m3);Z(t)為t時刻圍區(qū)內(nèi)水位(m);h(t)為t時刻外海側(cè)潮位(m);F[Z(t)]為圍區(qū)水位為Z(t)時的水面面積;mμ為龍口水力計(jì)算的流量系數(shù);Bω為實(shí)際等效過流計(jì)算寬度;k為參數(shù),堰流時k=1.5。

    由于外海潮位、圍區(qū)庫容等邊界條件為非線性性關(guān)系,因此,以上龍口水力計(jì)算方程為非線性常微分方程,目前還難以給出其解析解。

    3 基于四階Runge-Kutta法的龍口水力計(jì)算

    根據(jù)以上方程式(1)—式(4),代換后可得水力計(jì)算控制方程為

    (5)

    則式(5)可簡化為

    (6)

    可采用經(jīng)典的四階Runge-Kutta法求解控制方程得

    (7)

    其中:

    K1i=fti,Zti;

    K4i=fti+Δt,Zti+ΔtK3i。

    式中: Δt為計(jì)算時間段的步長;Zti為在ti=iΔt時刻的圍區(qū)水位;K1i,K2i,K3i,K4i為i時段的積分參數(shù),由Z(ti),ti等計(jì)算確定。

    根據(jù)以上數(shù)值計(jì)算法,可由外海側(cè)起調(diào)潮位、圍區(qū)側(cè)起調(diào)水位Zt0等邊界條件計(jì)算求解出圍區(qū)水位隨時間的過程線Z(t),其中Runge-Kutta法的求解流程如圖1所示,可利用MATLAB[8-9]軟件數(shù)值編程計(jì)算實(shí)現(xiàn),其中計(jì)算時間步長Δt一般可取0.1~0.5 h。

    圖1 Runge-Kutta法求解過程示意圖Fig.1 Sketch of the solution process of Runge-Kuttamethod

    其中深厚淤泥地基上的海堤龍口典型斷面如圖2所示。

    圖2 深厚淤泥地基上的海堤龍口典型斷面示意圖Fig.2 Schematic diagram of the closure gap section of sea dyke on deep silt foundation

    4 龍口最大水位差的解析分析

    由于外海潮位的近似周期性變化特性,假設(shè)外海側(cè)潮位過程線為

    (8)

    其中,

    式中:h為潮位;T0為潮位半周期,即潮位周期為2T0;H為最大潮差;χ為平均潮位。

    龍口、水閘等總過流能力的計(jì)算公式為

    (9)

    在忽略內(nèi)陸來水情況下,依據(jù)式(1)流量平衡方程得:

    (10)

    (11)

    其中圍區(qū)側(cè)起調(diào)水位取為Z0,參數(shù)Z0為海堤龍口底檻高程;值得說明的是,在落潮時需對式(11)取絕對值。

    將式(11)代入式(10)可得

    (12)

    按照工程計(jì)算經(jīng)驗(yàn),在最不利水力條件即圍區(qū)側(cè)水位與外海側(cè)潮位之間最大水位差情況(如圖3所示)[11]一般出現(xiàn)在高潮位之后,即滿足dH0/dt=0時,對式(11)兩邊求導(dǎo),變化可得

    (13)

    圖3 潮位與圍區(qū)內(nèi)水位關(guān)系曲線示意圖Fig.3 Relation between tide level and water level in the closure area

    在龍口水力計(jì)算中,主要水力參數(shù)是圍區(qū)水位與外海潮位間水位差最大值H0max,其中根據(jù)浙江省水利水電勘測設(shè)計(jì)院龍口水力計(jì)算經(jīng)驗(yàn),最大水位差出現(xiàn)的時間一般為tm=t0+(1.0~1.5)T0,即低潮位之前(0~0.5)T0。記龍口內(nèi)外兩側(cè)水位差最大時刻為tm,此時最大水位差為H0max,相應(yīng)圍區(qū)水位為Z(tm)。

    將式(13)代入式(12),則得

    (14)

    由于式(14)具有2個待求變量tm和Ztm,還需要1個方程才能求解,根據(jù)水量平衡建立方程,即

    (15)

    式中參數(shù)K為圍區(qū)庫容修正系數(shù),一般取值為1.0~1.2。

    由式(14)和式(15),通過數(shù)學(xué)變換,可得以下非線性代數(shù)方程組,即

    根據(jù)方程組(16),求解非線性方程組可解出待求參數(shù)tm和Z(tm)。

    (17)

    依據(jù)浙江省水利水電勘測設(shè)計(jì)院龍口水力計(jì)算經(jīng)驗(yàn),一般μ=1.0~1.50,借助MATLAB軟件計(jì)算求解出μ和tm,并根據(jù)式(8)和式(11)給出最大水位差H0max計(jì)算式。

    設(shè)圍區(qū)側(cè)的庫容曲線滿足二次拋物線,即

    V[Zt]=a2Z2t+a1Zt+a0。(18)

    對式(18)求導(dǎo)可知

    F[Zt]=2a2Zt+a1。

    (19)

    式中常參數(shù)a0,a1,a2根據(jù)實(shí)測庫容關(guān)系曲線擬合得出,將V[Z(t)]以及F[Z(t)]代入以上非線性代數(shù)方程(17),得到最大水位差時龍口圍區(qū)側(cè)水位與龍口寬龍口底坎高程之間的關(guān)系函數(shù)。

    外海潮位與圍區(qū)水位之間的最大水位差為

    (20)

    進(jìn)而依據(jù)《灘涂治理工程規(guī)范》[6]和《海堤工程設(shè)計(jì)規(guī)范》[7]可知最大流速為

    (21)

    根據(jù)式(17)、式(20)及式(21),可對龍口水力特性進(jìn)行計(jì)算分析,并能繪制漲潮、落潮時的圍區(qū)潮位最大水位差與龍口寬度關(guān)系曲線,進(jìn)而能給出龍口段漲潮、落潮時的最大水流流速等水力特性,為工程實(shí)際中深厚淤泥地基海堤龍口抗沖刷設(shè)計(jì)提供計(jì)算依據(jù)。

    5 算例分析

    5.1 項(xiàng)目概況

    某圍墾工程圍墾面積2 296.6 hm2,工程主要由海堤、水閘、施工道路等組成。根據(jù)防護(hù)對象的規(guī)模和重要性,確定工程防潮(洪)標(biāo)準(zhǔn)為50 a一遇。該工程海堤級別為3級,水閘級別同為3級,施工道路、河道等次要建筑物為4級,圍堰等臨時建筑物為5級。海堤地基土層主要為Ⅲ0層淤泥、Ⅲ1層淤泥質(zhì)粉質(zhì)黏土、ⅢSis粉砂、粉土混淤泥、Ⅲ2層淤泥、Ⅳ1淤泥質(zhì)粉質(zhì)黏土、Ⅳ2層粉質(zhì)黏土、黏土組成。

    其中龍口布置原則:①龍口位置宜考慮圍區(qū)涂面較低部位,以利于圍區(qū)潮流進(jìn)出;②龍口位置應(yīng)離水閘有一定距離,以免影響水閘施工和水閘泄流影響堵口;③龍口型式采用“寬高式龍口”龍口。根據(jù)龍口水力條件,為減小龍口流速防止沖刷和簡化保護(hù)措施,目前浙江省大中型圍墾工程和堵港工程的度汛龍口型式均采用“寬高式龍口”,即龍口寬度盡可能寬一些,底檻高程適當(dāng)抬高。

    5.2 龍口水力計(jì)算

    該工程圍堤龍口度汛設(shè)計(jì)潮型,以非汛期5 a一遇設(shè)計(jì)高潮位3.94 m為控制水位,將典型潮位過程作適當(dāng)修正后采用;圍堤度汛期設(shè)計(jì)潮位見表1。

    表1 龍口度汛期設(shè)計(jì)潮位過程線Table 1 Process of design tide level at closure gap in flood season

    其中通過MATLAB軟件對外海設(shè)計(jì)潮位進(jìn)行數(shù)值擬合分析得

    h(t)=0.76+2.67sin0.507t+1.407 。 (22)

    其中潮位過程線與擬合函數(shù)關(guān)系如圖4所示。

    圖4 龍口度汛期設(shè)計(jì)潮位過程線及數(shù)值擬合Fig.4 Process and numerical fitting of design tide level at closure gap in flood season

    該工程圍堤內(nèi)側(cè)庫容與圍區(qū)水位關(guān)系通過實(shí)測地形曲線量算,見表2。

    表2 圍區(qū)庫容-水位關(guān)系

    通過MATLAB軟件數(shù)值擬合,給出庫容-水位關(guān)系曲線(見圖5),擬合函數(shù)為

    V[Z(t)]=39.2Z2(t)+460Z(t)+453 。

    (23)

    圖5 圍區(qū)庫容-水位關(guān)系曲線及數(shù)值擬合Fig.5 Measured and fitted relation between storage capacity of closure area and water level

    根據(jù)本文方法,對該圍墾工程龍口水力條件進(jìn)行計(jì)算研究,其中漲潮、落潮時龍口最大流速Vmax與龍口底高程Z0、龍口寬度B之間關(guān)系曲線族見圖6。

    圖6 龍口最大流速Vmax與龍口底高程、 龍口寬度B之間關(guān)系曲線族Fig.6 Set of curves of maximum flow velocity Vmax at closure gap versus bottom height and width of closure gap

    通過以上數(shù)值計(jì)算分析可知:當(dāng)?shù)讬懜叱淘?.0以上時,龍口處漲落潮流速隨龍口寬度增大而減小,且流速變化較明顯;當(dāng)?shù)讬懜叱倘?.0 m以上時,龍口處漲落潮流速隨龍口寬度變化相對較小。

    根據(jù)龍口布置及控制漲落潮流速等一般要求,同時結(jié)合本工程涂面高程實(shí)際情況,以最大控制流速不大于3 m/s,通過以上數(shù)值計(jì)算,最終設(shè)計(jì)選定龍口底檻高程0.5 m,龍口設(shè)計(jì)寬為300 m。

    6 結(jié)論與建議

    針對深厚淤泥軟土灘涂的海堤龍口抗水流沖刷能力較差的問題,本文對海堤龍口水力學(xué)特性進(jìn)行了研究分析,主要結(jié)論如下:

    (1) 依據(jù)水力學(xué)理論及水量平衡原理建立了海堤龍口水力計(jì)算基本方程,該方程為非線性常微分方程,一般情況下難以給出精確解析解,可依據(jù)數(shù)值分析理論利用經(jīng)典四階Runge-Kutta法求解,結(jié)合MATLAB計(jì)算軟件數(shù)值編程給出數(shù)值解,其中計(jì)算時間步長Δt一般可取0.1~0.5 h。

    (2) 利用正弦函數(shù)擬合外海潮位周期變化特性,利用二次函數(shù)擬合圍區(qū)庫容-水位關(guān)系曲線,結(jié)合相關(guān)數(shù)值等效分析法及海堤龍口計(jì)算經(jīng)驗(yàn),給出了海堤龍口最大水位差與龍口寬度之間的近似解析計(jì)算方法,為工程實(shí)際應(yīng)用提供方便。

    (3)最后通過工程案例,利用MATLAB軟件計(jì)算對比分析,表明該案例中底檻高程在0.0以上時,龍口處漲落潮流速隨龍口寬度增大而減小,且流速變化較明顯;當(dāng)?shù)讬懜叱倘?.0 m以上時,龍口處漲落潮流速隨龍口寬度變化相對較小,與工程實(shí)際經(jīng)驗(yàn)一致。

    值得說明的是,由于漲落潮時口門上水流流態(tài)較復(fù)雜,涉及高淹沒出流情況,因此本文簡化計(jì)算方法難免存在一定不足,筆者將進(jìn)行后續(xù)深入研究。

    猜你喜歡
    海堤龍口潮位
    “龍口”治水
    金橋(2022年9期)2022-09-20 05:51:32
    基于距離倒數(shù)加權(quán)的多站潮位改正方法可行性分析
    2013年龍口歸城遺址調(diào)查、勘探簡報(bào)
    東方考古(2021年0期)2021-07-22 06:28:28
    唐山市警戒潮位標(biāo)志物維護(hù)研究
    『綠色海堤』扮美『飛地漁村』
    紅土地(2019年10期)2019-10-30 03:35:10
    多潮位站海道地形測量潮位控制方法研究
    水位驟降條件下海堤滲流及穩(wěn)定性分析——以寶安區(qū)西海堤為例
    基于改進(jìn)的OLS-RBF模型的感潮河段潮位預(yù)測研究
    基于風(fēng)潮聯(lián)合概率分布的海堤設(shè)防標(biāo)準(zhǔn)分析
    海堤地基加固處理方案比選及設(shè)計(jì)
    欧美在线一区亚洲| 看黄色毛片网站| av国产精品久久久久影院| 亚洲欧洲精品一区二区精品久久久| 日韩大尺度精品在线看网址 | 国产亚洲精品久久久久久毛片| 国产又爽黄色视频| 少妇裸体淫交视频免费看高清 | 午夜两性在线视频| 国产蜜桃级精品一区二区三区| 18禁美女被吸乳视频| 亚洲成a人片在线一区二区| 国产精品永久免费网站| 在线观看舔阴道视频| 女同久久另类99精品国产91| 在线观看一区二区三区| 香蕉丝袜av| 国产精华一区二区三区| 精品久久蜜臀av无| 韩国av一区二区三区四区| 91av网站免费观看| 久久香蕉激情| 欧美久久黑人一区二区| 一二三四在线观看免费中文在| 亚洲伊人色综图| 黄片小视频在线播放| 欧美日韩乱码在线| aaaaa片日本免费| 美女 人体艺术 gogo| 国产欧美日韩精品亚洲av| 99在线人妻在线中文字幕| 久久亚洲真实| 一进一出抽搐动态| 99精品在免费线老司机午夜| 妹子高潮喷水视频| 久久国产乱子伦精品免费另类| 亚洲av成人不卡在线观看播放网| 国产精品亚洲一级av第二区| 丰满饥渴人妻一区二区三| 变态另类成人亚洲欧美熟女 | 国产乱人伦免费视频| 亚洲 国产 在线| 精品国产美女av久久久久小说| 国产av精品麻豆| 香蕉国产在线看| 欧美中文综合在线视频| 免费日韩欧美在线观看| 亚洲成人免费电影在线观看| 国产成人欧美| 热re99久久精品国产66热6| 亚洲色图av天堂| 亚洲欧美一区二区三区久久| 日本免费一区二区三区高清不卡 | 黄片小视频在线播放| 亚洲av熟女| 国产麻豆69| 黄片播放在线免费| 动漫黄色视频在线观看| 搡老岳熟女国产| 国产一区二区三区综合在线观看| 欧美精品一区二区免费开放| 国产区一区二久久| 久久人人爽av亚洲精品天堂| 欧美 亚洲 国产 日韩一| 国产伦一二天堂av在线观看| 少妇粗大呻吟视频| 欧美黄色片欧美黄色片| 大香蕉久久成人网| 国产aⅴ精品一区二区三区波| 无限看片的www在线观看| 99久久久亚洲精品蜜臀av| 国产高清激情床上av| 1024香蕉在线观看| 久久天躁狠狠躁夜夜2o2o| 国产亚洲欧美98| 亚洲欧美日韩高清在线视频| 午夜激情av网站| 无遮挡黄片免费观看| 国产免费现黄频在线看| 欧美精品一区二区免费开放| 精品国产乱子伦一区二区三区| 欧美精品亚洲一区二区| 国产精品久久久av美女十八| 男女之事视频高清在线观看| 精品久久久久久电影网| 丁香六月欧美| 麻豆久久精品国产亚洲av | 在线观看免费午夜福利视频| 亚洲 欧美一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 丁香欧美五月| 久久人妻福利社区极品人妻图片| 国产区一区二久久| 韩国av一区二区三区四区| 成人av一区二区三区在线看| 亚洲第一av免费看| 51午夜福利影视在线观看| 后天国语完整版免费观看| 午夜成年电影在线免费观看| 日韩精品中文字幕看吧| 啦啦啦免费观看视频1| 免费高清在线观看日韩| 日本vs欧美在线观看视频| av网站在线播放免费| 欧美不卡视频在线免费观看 | 成人影院久久| 97碰自拍视频| 国产精品 国内视频| 亚洲精品国产色婷婷电影| 精品国产亚洲在线| 大陆偷拍与自拍| 国产单亲对白刺激| 夜夜看夜夜爽夜夜摸 | 黄色丝袜av网址大全| 亚洲三区欧美一区| 国产av精品麻豆| 一级片'在线观看视频| 在线观看66精品国产| 国产精品二区激情视频| 久久久久久亚洲精品国产蜜桃av| 一级黄色大片毛片| 久久久精品欧美日韩精品| 91麻豆av在线| 老汉色∧v一级毛片| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品美女久久av网站| 国产国语露脸激情在线看| 好看av亚洲va欧美ⅴa在| 99re在线观看精品视频| 久9热在线精品视频| 美女福利国产在线| 亚洲国产精品sss在线观看 | 又黄又爽又免费观看的视频| 两个人免费观看高清视频| 嫁个100分男人电影在线观看| 美女 人体艺术 gogo| 美女福利国产在线| 在线永久观看黄色视频| 国产成人一区二区三区免费视频网站| 欧美精品亚洲一区二区| 99在线人妻在线中文字幕| 国产区一区二久久| 妹子高潮喷水视频| 国产成人精品在线电影| 色综合欧美亚洲国产小说| 一二三四在线观看免费中文在| 88av欧美| 日本vs欧美在线观看视频| 欧美不卡视频在线免费观看 | 波多野结衣高清无吗| 国产精品免费一区二区三区在线| 成年人免费黄色播放视频| 免费在线观看日本一区| 亚洲av片天天在线观看| 亚洲欧美精品综合一区二区三区| 国产欧美日韩一区二区三区在线| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产精品sss在线观看 | 久久狼人影院| 日韩欧美三级三区| 大陆偷拍与自拍| 高清毛片免费观看视频网站 | 国产免费男女视频| 亚洲专区国产一区二区| 久久久久久人人人人人| 一进一出抽搐gif免费好疼 | 精品一区二区三区四区五区乱码| 色婷婷av一区二区三区视频| 性欧美人与动物交配| 丰满迷人的少妇在线观看| av电影中文网址| 国产深夜福利视频在线观看| 宅男免费午夜| 亚洲av第一区精品v没综合| 波多野结衣一区麻豆| 国产日韩一区二区三区精品不卡| 女性被躁到高潮视频| 久久久久久人人人人人| 国产欧美日韩一区二区三| 精品免费久久久久久久清纯| 亚洲欧美激情综合另类| 精品第一国产精品| 777久久人妻少妇嫩草av网站| 色在线成人网| 黄片大片在线免费观看| 中文字幕最新亚洲高清| 国产伦一二天堂av在线观看| 欧美激情 高清一区二区三区| 国产成人av教育| 亚洲av五月六月丁香网| 亚洲av日韩精品久久久久久密| 成年女人毛片免费观看观看9| 9191精品国产免费久久| 久久久久国产精品人妻aⅴ院| 欧美精品亚洲一区二区| 自线自在国产av| 国产精品久久久av美女十八| 男女床上黄色一级片免费看| 欧美日韩亚洲综合一区二区三区_| 亚洲人成电影免费在线| 亚洲情色 制服丝袜| 久久精品亚洲精品国产色婷小说| 欧美+亚洲+日韩+国产| 久久天堂一区二区三区四区| 精品福利永久在线观看| 大码成人一级视频| av电影中文网址| 国产av在哪里看| 女性生殖器流出的白浆| 在线观看免费午夜福利视频| 美女扒开内裤让男人捅视频| 日本黄色日本黄色录像| 久久久久久大精品| 757午夜福利合集在线观看| 亚洲 国产 在线| 国产欧美日韩精品亚洲av| 男人操女人黄网站| 午夜成年电影在线免费观看| 欧美国产精品va在线观看不卡| 女人被躁到高潮嗷嗷叫费观| 日本精品一区二区三区蜜桃| 亚洲人成电影观看| 国产av一区在线观看免费| 不卡av一区二区三区| 天天影视国产精品| 免费观看人在逋| 精品第一国产精品| 亚洲成人精品中文字幕电影 | 亚洲色图综合在线观看| 亚洲成人国产一区在线观看| 国产精品久久电影中文字幕| 97碰自拍视频| 午夜精品久久久久久毛片777| 亚洲一区中文字幕在线| 久久久精品欧美日韩精品| av天堂久久9| 欧美+亚洲+日韩+国产| 国产av精品麻豆| 操美女的视频在线观看| 我的亚洲天堂| 欧美成人午夜精品| 国产不卡一卡二| 满18在线观看网站| 国产成人欧美在线观看| 一本大道久久a久久精品| 国产亚洲精品第一综合不卡| 国产精品一区二区三区四区久久 | 国产激情欧美一区二区| 超碰97精品在线观看| 色综合站精品国产| 校园春色视频在线观看| 国产野战对白在线观看| bbb黄色大片| 免费日韩欧美在线观看| 嫩草影院精品99| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品999在线| 久99久视频精品免费| 变态另类成人亚洲欧美熟女 | aaaaa片日本免费| 亚洲熟女毛片儿| 久久人妻熟女aⅴ| 人妻丰满熟妇av一区二区三区| 日韩精品免费视频一区二区三区| 欧美黄色淫秽网站| 亚洲av成人不卡在线观看播放网| 香蕉国产在线看| 好看av亚洲va欧美ⅴa在| 岛国视频午夜一区免费看| 色综合欧美亚洲国产小说| 嫩草影视91久久| 中文欧美无线码| 国产乱人伦免费视频| 国产av精品麻豆| 亚洲欧美日韩高清在线视频| 国产主播在线观看一区二区| 欧美丝袜亚洲另类 | 亚洲性夜色夜夜综合| svipshipincom国产片| 亚洲成人精品中文字幕电影 | 动漫黄色视频在线观看| 视频在线观看一区二区三区| 国产精品一区二区精品视频观看| 高清黄色对白视频在线免费看| 国产亚洲精品综合一区在线观看 | 91精品国产国语对白视频| 久久中文看片网| 久久久国产成人精品二区 | 成熟少妇高潮喷水视频| 久久久久亚洲av毛片大全| 成人黄色视频免费在线看| 又紧又爽又黄一区二区| 色精品久久人妻99蜜桃| 美女国产高潮福利片在线看| 日本a在线网址| 热re99久久国产66热| 日韩欧美三级三区| 久久热在线av| 久久精品国产亚洲av香蕉五月| www国产在线视频色| 自线自在国产av| 制服诱惑二区| 亚洲精品中文字幕在线视频| 国产精华一区二区三区| 久久精品91蜜桃| а√天堂www在线а√下载| 亚洲午夜理论影院| 天天添夜夜摸| 脱女人内裤的视频| 亚洲激情在线av| 亚洲精品国产一区二区精华液| 夜夜看夜夜爽夜夜摸 | 精品国产乱子伦一区二区三区| 久久狼人影院| av福利片在线| 9色porny在线观看| a在线观看视频网站| 国产高清国产精品国产三级| 欧美精品啪啪一区二区三区| 嫩草影视91久久| 99国产综合亚洲精品| 亚洲成a人片在线一区二区| 热99国产精品久久久久久7| 一区二区日韩欧美中文字幕| 精品福利永久在线观看| 黑人欧美特级aaaaaa片| 老汉色∧v一级毛片| av有码第一页| 亚洲免费av在线视频| 欧美激情久久久久久爽电影 | 亚洲全国av大片| 国产av精品麻豆| 美女 人体艺术 gogo| 亚洲va日本ⅴa欧美va伊人久久| 人人妻人人爽人人添夜夜欢视频| 欧美中文综合在线视频| 午夜福利一区二区在线看| 久久草成人影院| 久久久久亚洲av毛片大全| 一进一出抽搐动态| 波多野结衣高清无吗| 久99久视频精品免费| 12—13女人毛片做爰片一| 国产色视频综合| 亚洲成人免费电影在线观看| 青草久久国产| 女同久久另类99精品国产91| 在线观看免费日韩欧美大片| 超碰97精品在线观看| 久久精品亚洲av国产电影网| 一区二区三区激情视频| 69av精品久久久久久| 久久久久国产精品人妻aⅴ院| 久久久久久久午夜电影 | 国产亚洲精品第一综合不卡| 美女扒开内裤让男人捅视频| 欧美久久黑人一区二区| 一边摸一边抽搐一进一出视频| 在线视频色国产色| a在线观看视频网站| 久久国产亚洲av麻豆专区| 日韩成人在线观看一区二区三区| 国产精品 国内视频| 高清在线国产一区| 97碰自拍视频| 满18在线观看网站| xxx96com| av网站免费在线观看视频| bbb黄色大片| 亚洲国产精品合色在线| 久久中文看片网| 国产精品亚洲av一区麻豆| 午夜视频精品福利| 黄片大片在线免费观看| 亚洲一区中文字幕在线| av在线天堂中文字幕 | 亚洲精品av麻豆狂野| 在线观看免费视频网站a站| 日韩一卡2卡3卡4卡2021年| 一区二区三区激情视频| 国产成人精品无人区| 色播在线永久视频| 长腿黑丝高跟| 最近最新中文字幕大全电影3 | 99热国产这里只有精品6| 欧美成人午夜精品| 极品教师在线免费播放| 国产激情欧美一区二区| 久久中文看片网| 人妻久久中文字幕网| 悠悠久久av| 老鸭窝网址在线观看| 午夜91福利影院| 99精品欧美一区二区三区四区| 精品日产1卡2卡| 涩涩av久久男人的天堂| 午夜精品在线福利| 亚洲午夜精品一区,二区,三区| 大陆偷拍与自拍| 日日夜夜操网爽| 亚洲片人在线观看| 国产麻豆69| 亚洲 欧美 日韩 在线 免费| 国产成人精品久久二区二区91| 日日爽夜夜爽网站| 国产91精品成人一区二区三区| 国产亚洲欧美98| 丰满饥渴人妻一区二区三| 国产伦人伦偷精品视频| xxxhd国产人妻xxx| 日韩欧美三级三区| 欧美日韩亚洲国产一区二区在线观看| 每晚都被弄得嗷嗷叫到高潮| 久久天躁狠狠躁夜夜2o2o| 黑人欧美特级aaaaaa片| 久久国产乱子伦精品免费另类| 一级作爱视频免费观看| 午夜91福利影院| 免费av毛片视频| 久99久视频精品免费| 变态另类成人亚洲欧美熟女 | 欧美激情极品国产一区二区三区| 婷婷丁香在线五月| 12—13女人毛片做爰片一| 亚洲一码二码三码区别大吗| 好看av亚洲va欧美ⅴa在| 日本黄色日本黄色录像| 国产激情久久老熟女| 婷婷精品国产亚洲av在线| 大香蕉久久成人网| 99国产精品一区二区三区| 久久精品人人爽人人爽视色| 精品久久久久久久久久免费视频 | 黄频高清免费视频| 中文字幕高清在线视频| 黄片小视频在线播放| 日本三级黄在线观看| 国产亚洲精品第一综合不卡| 人人妻人人澡人人看| 欧美亚洲日本最大视频资源| 母亲3免费完整高清在线观看| 久久精品亚洲精品国产色婷小说| 国产精品日韩av在线免费观看 | 一区二区三区国产精品乱码| 午夜福利,免费看| 免费av中文字幕在线| 国产不卡一卡二| 午夜老司机福利片| 欧美日韩黄片免| 国产欧美日韩一区二区精品| 欧美亚洲日本最大视频资源| 久久午夜亚洲精品久久| 亚洲欧美精品综合一区二区三区| 性色av乱码一区二区三区2| 亚洲色图av天堂| 亚洲一区二区三区欧美精品| 久久久水蜜桃国产精品网| 黄片播放在线免费| 天天添夜夜摸| 国产精品免费视频内射| 久久香蕉激情| 美女大奶头视频| 日韩欧美在线二视频| 亚洲男人天堂网一区| 99国产精品99久久久久| 国产成人精品久久二区二区91| 黄色成人免费大全| 一级作爱视频免费观看| 美女大奶头视频| 成人影院久久| 久久久久久久午夜电影 | 久久性视频一级片| 99久久人妻综合| 人妻久久中文字幕网| 日韩欧美三级三区| 在线观看舔阴道视频| 久久久久久久久久久久大奶| 欧美日韩亚洲高清精品| 久久久久国产一级毛片高清牌| 精品国产美女av久久久久小说| 日日爽夜夜爽网站| 亚洲av熟女| 精品一区二区三区视频在线观看免费 | 欧美另类亚洲清纯唯美| 一区二区日韩欧美中文字幕| 可以在线观看毛片的网站| 国产伦人伦偷精品视频| 夜夜看夜夜爽夜夜摸 | 黄色 视频免费看| 最新在线观看一区二区三区| 国产91精品成人一区二区三区| 看片在线看免费视频| 中文字幕人妻丝袜制服| 日本黄色视频三级网站网址| 两人在一起打扑克的视频| 国产av一区二区精品久久| 日日夜夜操网爽| 性色av乱码一区二区三区2| 亚洲熟女毛片儿| 精品久久久久久电影网| 一区二区日韩欧美中文字幕| 国产1区2区3区精品| 国产伦人伦偷精品视频| 亚洲成人免费电影在线观看| 精品一区二区三区四区五区乱码| 美女午夜性视频免费| 制服人妻中文乱码| 欧美日韩中文字幕国产精品一区二区三区 | 脱女人内裤的视频| 日韩视频一区二区在线观看| 亚洲va日本ⅴa欧美va伊人久久| 这个男人来自地球电影免费观看| 国产成人欧美| 老司机福利观看| 国产欧美日韩精品亚洲av| 中文字幕精品免费在线观看视频| 午夜免费激情av| 人人妻,人人澡人人爽秒播| 黑人猛操日本美女一级片| 琪琪午夜伦伦电影理论片6080| 日韩精品中文字幕看吧| 麻豆国产av国片精品| 日本欧美视频一区| 精品国内亚洲2022精品成人| 亚洲国产精品一区二区三区在线| 黄色 视频免费看| 精品欧美一区二区三区在线| 国产欧美日韩精品亚洲av| 欧美性长视频在线观看| 男人操女人黄网站| 午夜福利欧美成人| 一级毛片精品| 亚洲黑人精品在线| 自线自在国产av| 国内久久婷婷六月综合欲色啪| 视频区欧美日本亚洲| 国产精品免费视频内射| 久久99热6这里只有精品| 90打野战视频偷拍视频| av福利片在线观看| 亚洲人成伊人成综合网2020| 亚洲激情在线av| 真人做人爱边吃奶动态| 听说在线观看完整版免费高清| 久久久久久久久久成人| 亚洲内射少妇av| 男插女下体视频免费在线播放| .国产精品久久| 亚洲av免费高清在线观看| 久久精品国产亚洲av涩爱 | 俺也久久电影网| 在线a可以看的网站| 午夜精品在线福利| 变态另类成人亚洲欧美熟女| 性插视频无遮挡在线免费观看| 麻豆一二三区av精品| 亚洲av一区综合| 日本成人三级电影网站| 亚洲欧美日韩高清专用| 精品人妻1区二区| 精品日产1卡2卡| 少妇裸体淫交视频免费看高清| 久久久精品大字幕| 亚洲av美国av| 一本精品99久久精品77| 美女高潮喷水抽搐中文字幕| 亚洲av中文字字幕乱码综合| 日韩欧美三级三区| 直男gayav资源| 欧美黑人欧美精品刺激| 中文字幕久久专区| 一级毛片久久久久久久久女| 老司机深夜福利视频在线观看| 欧美成人a在线观看| 久久精品国产亚洲av涩爱 | 18禁裸乳无遮挡免费网站照片| 成人特级黄色片久久久久久久| 高清毛片免费观看视频网站| 国产免费一级a男人的天堂| 久9热在线精品视频| 神马国产精品三级电影在线观看| 99热这里只有是精品在线观看 | 老司机午夜十八禁免费视频| 日韩成人在线观看一区二区三区| 久久久国产成人精品二区| 国产伦精品一区二区三区视频9| 一卡2卡三卡四卡精品乱码亚洲| 色综合亚洲欧美另类图片| 亚洲在线自拍视频| 午夜福利免费观看在线| 熟妇人妻久久中文字幕3abv| 国产高清有码在线观看视频| 在线a可以看的网站| 亚洲国产色片| 久久国产乱子免费精品| 亚洲天堂国产精品一区在线| 亚洲av成人av| 9191精品国产免费久久| 狠狠狠狠99中文字幕| 小说图片视频综合网站| 久9热在线精品视频| 亚洲自拍偷在线| 国内久久婷婷六月综合欲色啪| АⅤ资源中文在线天堂| 久久久久久久久久成人| 18美女黄网站色大片免费观看| 悠悠久久av| 亚洲精品456在线播放app | 精品乱码久久久久久99久播| 国产精品电影一区二区三区| 少妇裸体淫交视频免费看高清| 国产老妇女一区| 亚洲人与动物交配视频| 国产一区二区在线观看日韩| 99久久成人亚洲精品观看|