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

    基于豎縫隔板的大底坡魚道水力優(yōu)化數(shù)值模擬

    2023-11-28 02:14:04黃姣鳳尹龍斌張家強(qiáng)蔣昌波
    水利水電科技進(jìn)展 2023年6期
    關(guān)鍵詞:豎縫魚道隔板

    鄧 斌,黃姣鳳,徐 拓,尹龍斌,張家強(qiáng),蔣昌波

    (1.長沙理工大學(xué)水利與環(huán)境工程學(xué)院,湖南 長沙 410114; 2.水沙科學(xué)與水災(zāi)害防治湖南省重點(diǎn)實(shí)驗(yàn)室,湖南 長沙 410114; 3.洞庭湖水環(huán)境治理與生態(tài)修復(fù)湖南省重點(diǎn)實(shí)驗(yàn)室,湖南 長沙 410114)

    閘、壩等水利設(shè)施的興建在給防洪、灌溉、緩解缺水、開發(fā)清潔能源等方面帶來巨大效益[1]的同時會截?cái)嗪恿?導(dǎo)致完整的河流碎片化[2],阻斷魚類洄游通道并影響其種群繁殖[3]。魚道作為一種生態(tài)型水利工程,能夠重建河流連通性,幫助魚類克服障礙恢復(fù)自然洄游[4],對保護(hù)魚類資源及水生生物多樣性具有一定的積極作用。魚道主要有仿自然型魚道(仿自然魚梯、魚坡和過魚旁路等)和技術(shù)型魚道(豎縫式魚道、池堰式魚道和丹尼爾魚道等)兩種類型[5],其中豎縫式魚道因結(jié)構(gòu)簡單、能適應(yīng)多變的上下游水位、消能效果顯著等特點(diǎn),成為應(yīng)用最廣泛和最有效的魚道型式之一[6]。

    研究發(fā)現(xiàn),我國已建的豎縫式魚道大多底坡較緩,為1/100~1/50[7]。而底坡坡度決定了魚道的總長度,是影響魚類上溯時間、工程建設(shè)成本及規(guī)劃布置的主要因素。一般來說,魚道底坡越陡,總長度越短,占地面積就越小,建設(shè)成本也越低。但是大底坡結(jié)構(gòu)所形成的高速度、高湍動水流流場不利于游泳能力較弱的魚類洄游上溯[8]。因此,考慮河流生態(tài)修復(fù)及工程經(jīng)濟(jì)性,國內(nèi)外學(xué)者針對大底坡下豎縫式魚道的水力優(yōu)化開展了相關(guān)研究,如部分學(xué)者基于物理模型試驗(yàn)[9-10]和數(shù)值模擬[11-14]方法,分別研究了在各種大底坡(1/20~1/7)結(jié)構(gòu)下,豎縫式魚道豎縫下游處設(shè)置不同直徑、個數(shù)及排列方式的圓柱體障礙物對池室水力特性的影響,并通過比較池室流速和池室湍動特性等相關(guān)參數(shù)來評估不同圓柱體配置方式的阻流效果。此外,也有部分學(xué)者研究了在豎縫處增設(shè)隔板對池室水流分布的影響,如Pena等[15]認(rèn)為在底坡為1/10的豎縫式魚道豎縫底部設(shè)置不同長度的隔板可以減小池室最大流速,但同時會產(chǎn)生比較復(fù)雜的流場結(jié)構(gòu)并增加池室湍動能;Mallen-Cooper等[16]基于原型觀測研究了底坡為1/20時在豎縫式魚道豎縫中間位置設(shè)置隔板對魚道過魚效率的影響,發(fā)現(xiàn)與未經(jīng)改造的豎縫式魚道相比,其小魚通過效率提升了6~13倍,這種通過在豎縫式魚道結(jié)構(gòu)水流最集中的部位設(shè)置消能措施以減少池室流速和湍流的方法比較有效、簡單和經(jīng)濟(jì),但是鮮見其相關(guān)的池室水力特性研究報(bào)道,而水力特征分布是探討魚類上溯行為機(jī)理的依據(jù),也亟須對其進(jìn)行相關(guān)研究。

    綜上所述,通過在豎縫式魚道結(jié)構(gòu)內(nèi)增設(shè)局部消能措施(如圓柱體、豎縫隔板等)可改變池室流場分布,并且關(guān)于不同配置方式的圓柱體對池室水力特性的影響已進(jìn)行了一定的研究,然而對于豎縫隔板所產(chǎn)生的影響,目前相關(guān)研究尚不夠全面、系統(tǒng),還有待進(jìn)一步探索與分析。因此,本文從水力、工程經(jīng)濟(jì)等角度出發(fā),在Mallen-Cooper等[16]的研究基礎(chǔ)上,通過建立魚道池室三維水動力數(shù)值模型,深入分析了5種大底坡J(1/40、1/30、1/20、1/15、1/10)條件下,豎縫隔板及其不同安裝高度h(0、0.04、0.08m,為豎縫隔板底部到池室底板的距離)對池室水力特性的影響,并闡明了豎縫隔板影響和優(yōu)化大底坡魚道池室流場分布的內(nèi)在原因。通過對池室水力特性的數(shù)值模擬研究,以期獲取關(guān)于這類魚道結(jié)構(gòu)運(yùn)行的基本知識,可為改進(jìn)更有效、經(jīng)濟(jì)的大底坡魚道結(jié)構(gòu)提供相關(guān)技術(shù)依據(jù)。

    1 數(shù)學(xué)模型

    1.1 控制方程

    本研究采用計(jì)算流體力學(xué)商業(yè)軟件FlOW-3D對豎縫式魚道內(nèi)的水動力進(jìn)行模擬,采用有限體積法求解三維雷諾平均納維-斯托克斯(RANS)方程。湍流模型選用重整化群(RNG)k-ε模型,該模型能夠較準(zhǔn)確地模擬復(fù)雜水工結(jié)構(gòu)內(nèi)的流場分布[17]。

    采用流體體積(VOF)方法模擬自由液面,即氣相和液相之間的界面。VOF方法通過求解整個計(jì)算域內(nèi)各網(wǎng)格單元的體積分?jǐn)?shù),并進(jìn)行持續(xù)記錄,構(gòu)建出自由液面。其表達(dá)式為

    (1)

    式中:aw為水的體積分?jǐn)?shù);t為時間;ui為流體速度;xi為坐標(biāo)。當(dāng)aw=0時,計(jì)算域內(nèi)全部充滿空氣;當(dāng)aw=1時,計(jì)算域內(nèi)全部充滿水;當(dāng)0

    1.2 計(jì)算域及網(wǎng)格劃分

    數(shù)學(xué)模型計(jì)算域包括魚道入流段、出流段及5級池室,如圖1所示。豎縫寬度b為0.05m,單個池室凈長L=10b=0.5m,凈寬B=8b=0.4m,池室長寬比為10/8,底坡J分別取1/40、1/30、1/20、1/15、1/10。豎縫隔板尺寸為0.05m×0.01m×0.05m(長×寬×高),其中高度與豎縫寬度相等,布置方式如圖2所示。為減少進(jìn)出口邊界對池室流場的影響,取第三級池室作為主要研究對象,采用六面體結(jié)構(gòu)化網(wǎng)格對模型進(jìn)行網(wǎng)格劃分,網(wǎng)格單元邊長設(shè)為0.006m。

    圖1 魚道數(shù)學(xué)模型平面示意圖

    圖2 豎縫隔板布置三維示意圖

    1.3 邊界條件及初始條件

    考慮到水流流動的物理特性及工程模擬的相似性,魚道上游進(jìn)口及下游出口邊界條件均采用壓力邊界,并設(shè)置相應(yīng)的水深,魚道進(jìn)口水深為0.195m,出口水深為0.245m;模型頂部設(shè)置為壓力進(jìn)口邊界,其相對壓強(qiáng)大小為0,水體體積分?jǐn)?shù)設(shè)置為0,以便空氣可自由進(jìn)入;魚道左右邊墻及底部設(shè)置為固壁無滑移邊界。模型初始水深H設(shè)置為0.2m,以使計(jì)算更快收斂并更加符合實(shí)際情況。模型采用廣義最小殘量(GMRES)隱式求解器計(jì)算,初始時間步長設(shè)置為0.001s,最小時間步長設(shè)置為10-7s。

    1.4 數(shù)值模擬與物理模型試驗(yàn)結(jié)果對比分析

    為了評估數(shù)學(xué)模型模擬魚道池室水動力特性的能力,通過調(diào)整數(shù)學(xué)模型中魚道邊墻及底部的表面粗糙度rs(0.010、0.012、0.015、0.018),采用與數(shù)學(xué)模型尺寸相一致的物理模型對各表面粗糙度下的流速進(jìn)行了對比驗(yàn)證。物理模型試驗(yàn)在長沙理工大學(xué)水利實(shí)驗(yàn)中心的變坡水槽中進(jìn)行,水槽尺寸為16m×0.4m×0.5m(長×寬×高),水槽底面及兩側(cè)均為透明的鋼化玻璃。為得到充分發(fā)展的水流,將模型布置在水槽中后段部分,并通過多普勒聲學(xué)流速儀測量了第三級池室z=0.5H水深平面處的速度場分布,池室流速測點(diǎn)布置如圖3所示。對比分析了底坡為1/20且無豎縫隔板情況下,距上游長擋板背水面距離0.22m斷面線處的物理模型與數(shù)學(xué)模型流速分布(圖4(a))及不同水深豎縫斷面中心點(diǎn)連線處的流速分布(圖4(b)),并通過決定系數(shù)R2對數(shù)學(xué)模型計(jì)算精度進(jìn)行評價。驗(yàn)證結(jié)果表明,當(dāng)rs=0.015時,數(shù)學(xué)模型結(jié)果與物理模型結(jié)果吻合較好,R2都在0.9以上,數(shù)學(xué)模型能夠較準(zhǔn)確地模擬魚道池室內(nèi)的水力特性。

    圖3 第三級池室測點(diǎn)布置(單位:m)

    2 結(jié)果與分析

    根據(jù)上述數(shù)學(xué)模型,分別對各工況進(jìn)行了模擬計(jì)算,為減少魚道進(jìn)出口邊界條件對計(jì)算結(jié)果的影響,以中間第三級池室作為主要研究對象進(jìn)行水力分析與討論。

    2.1 不同底坡池室水力參數(shù)對比

    在豎縫式魚道的水力設(shè)計(jì)中,池室最大流速vmax和池室湍動能k是影響過魚效果最主要的水力參數(shù)[18]。由于豎縫式魚道在相鄰池室間會產(chǎn)生水位落差,vmax一般出現(xiàn)在豎縫位置附近,而豎縫是魚類進(jìn)行洄游的必經(jīng)之地,因此,為使魚類能夠順利通行,vmax必須小于魚的爆發(fā)速度[19]。k表示脈動流速的動能,是衡量水流湍動程度的重要指標(biāo)[20];若魚類上溯路徑中的水流k過大,會導(dǎo)致魚類特別是游泳能力較弱的魚類在上溯過程中難以保持方向,同時增加魚類游泳耗能甚至對魚類造成不同程度的傷害[21],這表明池室中適當(dāng)?shù)膋分布也是魚類成功上溯的關(guān)鍵。因此,優(yōu)化大底坡魚道池室水力分布的主要目標(biāo)是減小池室最大流速和池室湍流。

    表1 不同底坡最大流速和平均湍動能下降率

    2.2 豎縫隔板對池室水力特性的影響

    為了進(jìn)一步闡明豎縫隔板對魚道池室流場分布的影響機(jī)理,選取第三級池室距底板0.03m(豎縫隔板下方水深處,除h=0m時)及0.15m(豎縫隔板上方水深處)的水平面作為主要研究對象,具體分析了底坡為1/20時,不同豎縫隔板配置方式在Q=19.36m3/h和Q=10.23m3/h(Q為魚道入口流量)下的池室水力特性分布。

    2.2.1流態(tài)

    圖5為不同流量條件下豎縫隔板不同配置方式的三維流線分布:當(dāng)豎縫處無隔板時(圖5(a)(e)),水流通過上游豎縫后,以射流的形式進(jìn)入池室,并沿池室隔板導(dǎo)角方向向前流動,射流到達(dá)池室中心后受下游擋板墩頭影響,其斷面逐漸收縮并同時折向下游豎縫進(jìn)入下一級池室,呈現(xiàn)出反S形的主流,并在主流兩側(cè)產(chǎn)生了兩個尺度相當(dāng)?shù)幕亓鲄^(qū);而對于大底坡和流量較大的情況,這種流型的主流流速(0.5~0.6m/s)過大,會導(dǎo)致魚類上溯耗能的增加并不利于游泳能力較弱的魚類通行。豎縫處設(shè)置隔板后(圖6(b)(c)(d)(f)(g)(h)),水流也以射流形式流入池室,但整體流動模式發(fā)生了顯著變化:豎縫隔板的存在改變了射流流向,形成了利用池室整體容積進(jìn)行低湍流能量轉(zhuǎn)換、自身回轉(zhuǎn)的主流,使得主流路徑更加彎曲并呈Ω形,該流型通過射流混合消能的方式,可有效延長大底坡、大流量條件下的水流能量耗散距離,在池室內(nèi)形成面積較大的低流速區(qū)(0.25~0.35m/s),使魚類能以持續(xù)游泳速度和耐久游泳速度相交替的游泳方式[22]游過各級魚道池室,有利于減小魚類上溯耗能;同時,在池室的上游角和主流的凸出部分產(chǎn)生了兩個尺度不一的回流區(qū),給魚類洄游提供休息區(qū)域。而改變豎縫隔板安裝高度,流態(tài)呈現(xiàn)相似的分布規(guī)律。

    圖6 不同水深處池室流速分布

    圖6為Q=19.36m3/h時豎縫隔板不同安裝情況下不同水深處的池室流速分布。與無隔板相比,3種隔板安裝高度下的池室流速都顯著下降,且具有相似的分布情況;當(dāng)h為0時,在0.03m水深處(圖6(b)),由于豎縫隔板的阻擋,其豎縫下游處的流速分布較緩;而豎縫隔板上下水深處的流速分布(圖6(c)(d)(f)(g)(h))大體保持一致,這表明改變豎縫隔板安裝高度,對其上下水深處的流速分布無顯著影響。

    2.2.2主流區(qū)水力特性

    圖7 最大流速軌跡線

    圖8 最大流速沿程分布曲線

    2.2.3豎縫斷面流速

    分析發(fā)現(xiàn),不同流量條件下豎縫斷面流速分布規(guī)律具有較好的相似性。因此,圖9僅給出了Q=19.36m3/h時H=0.03m和H=0.15m處豎縫斷面流速沿豎縫寬度b方向的分布規(guī)律。在豎縫處設(shè)置隔板使豎縫中間位置(b′/b=0.3~0.9,其中b′為豎縫上的點(diǎn)到隔板的距離)的流速顯著下降,而對靠近兩側(cè)擋板處的流速影響較小;同時,當(dāng)h=0m時,魚道底部由于豎縫隔板的阻擋,使得魚道底部(H=0.03m)豎縫流速為0,因此圖9(a)未給出其流速分布,當(dāng)h=0.04m時,對豎縫隔板上方水深處(H=0.15m)的流速影響較小,當(dāng)h=0.08m時,豎縫隔板上下水深處(H=0.03m、0.15m)的流速都有明顯降低,這表明豎縫隔板能顯著減小其安裝位置附近水深處的豎縫斷面流速,且h越大,這種影響范圍也越大?;诖?初步認(rèn)為可將隔板安裝在豎縫的中上方水深處,以保障魚類能夠更加順利地通過豎縫完成洄游。

    圖9 不同水深處豎縫斷面流速分布

    2.2.4湍動能

    圖10為Q=19.36m3/h時豎縫隔板不同安裝情況在不同水深處的池室湍動能分布:不同工況池室湍動能的分布與其所對應(yīng)的池室流型分布具有較好的一致性;在不同水深處(除豎縫隔板阻擋水深處),湍動能最大值均出現(xiàn)在豎縫附近,同時水流通過豎縫后由于主流區(qū)流速較大,兩側(cè)回流區(qū)流速較小,造成主流邊界流速梯度較大,因此在主流及其左右邊界的剪切區(qū)域亦表現(xiàn)出較大的湍動能,并隨著水深的增加,池室中湍動能的分布規(guī)律大致相同,但相同位置的湍動能有所增大;與無隔板相比,豎縫處設(shè)置隔板可降低池室湍動能的大小,并使其在豎縫下游位置區(qū)域內(nèi)迅速衰減,產(chǎn)生連續(xù)的低湍動區(qū)域(0.007~0.01J/kg),這給魚類特別是游泳能力較弱的魚類提供了一個較好的上溯條件;同流速分布一樣,改變豎縫隔板安裝高度,對豎縫隔板上下水深處池室湍動能的分布規(guī)律無顯著影響。

    圖10 不同水深處池室湍動能分布

    2.2.5應(yīng)變率

    應(yīng)變率(sr)為水流各方向速度梯度的絕對值之和。相關(guān)研究[24-25]認(rèn)為水流速度梯度的大小是影響魚類側(cè)線系統(tǒng)感知水流刺激并產(chǎn)生運(yùn)動的重要因素,因此,應(yīng)變率也是研究魚道水力特性的重要參數(shù)之一。圖11為Q=19.36m3/h時不同隔板安裝情況下不同水深處的池室應(yīng)變率分布。與池室湍動能分布一樣,應(yīng)變率分布也取決于池室流型的分布,并沿主流兩側(cè)邊界有逐漸減小的趨勢;同時在豎縫位置內(nèi)側(cè)以及靠近豎縫處下游位置的兩側(cè)呈現(xiàn)相對較大的數(shù)值,即在這些區(qū)域水流速度變化顯著、流速梯度大(豎縫隔板阻擋水深處除外);設(shè)置豎縫隔板會使應(yīng)變率的分布在豎縫下游處發(fā)生偏轉(zhuǎn),但對應(yīng)變率的大小沒有顯著影響,這表明魚類能在小流速和低湍動的水流條件下感知到相同的水流刺激向前運(yùn)動,可促進(jìn)魚類的洄游上溯;而改變豎縫隔板的安裝高度,對豎縫隔板上下水深處應(yīng)變率的分布規(guī)律無明顯影響。

    圖11 不同水深處池室水力應(yīng)變率分布

    3 結(jié) 論

    a.對于不同的大底坡坡度,豎縫處設(shè)置隔板可有效減小池室最大流速及池室平均湍動能:池室最大流速下降率分布范圍為2.01%~21.49%,池室平均湍動能下降率分布范圍為3.77%~42.19%。

    b.底坡為1/20時,池室主流從反S形轉(zhuǎn)變?yōu)槔贸厥艺w容積進(jìn)行低湍流能量轉(zhuǎn)換并自身回轉(zhuǎn)的Ω形,并且其主流最大流速衰減率達(dá)到61.24%~61.82%。

    c.改變豎縫隔板安裝高度對隔板上下水深處的豎縫流速分布存在一定的影響,而對流速、湍動能、水力應(yīng)變率的分布規(guī)律無顯著影響,因此為了使魚類能夠更加順利地通過豎縫,可將豎縫隔板安裝在豎縫的中上方水深位置。

    d.在豎縫處設(shè)置隔板以優(yōu)化大底坡魚道形成的高流速、高湍動的水流流場具有一定的可行性,但是還需進(jìn)一步結(jié)合過魚試驗(yàn)研究來評估魚類對這些水動力條件的上溯行為反應(yīng),以論證研究結(jié)果的工程應(yīng)用價值。

    猜你喜歡
    豎縫魚道隔板
    豎縫位置對異側(cè)豎縫式魚道水力特性的數(shù)值模擬與分析
    基于PIT遙測技術(shù)的豎縫式魚道過魚效率及魚類行為分析
    基于有限元分析的隔板工藝支撐技術(shù)研究
    壓力容器隔板的一種設(shè)計(jì)方法
    蜿蜒式魚道與豎縫式魚道的水力特性對比分析
    裝配式建筑剪力墻豎縫連接形式及其研究現(xiàn)狀
    河南建材(2021年3期)2021-01-07 07:00:41
    異側(cè)豎縫寬度對魚道水流結(jié)構(gòu)的影響研究
    紊流特性對魚道結(jié)構(gòu)的影響
    卷宗(2020年16期)2020-08-10 18:21:27
    帶豎縫鋼管束砼組合剪力墻受力性能模擬分析
    橫隔板參數(shù)對裝配式箱梁寬橋受力的影響
    久久久久人妻精品一区果冻| 男女边摸边吃奶| 久热这里只有精品99| 三上悠亚av全集在线观看| 国产成人aa在线观看| 91精品国产国语对白视频| 国产片内射在线| 啦啦啦在线观看免费高清www| 欧美bdsm另类| 亚洲精品aⅴ在线观看| 亚洲第一区二区三区不卡| 18禁国产床啪视频网站| 久久狼人影院| 黑人巨大精品欧美一区二区蜜桃 | 97超碰精品成人国产| 乱码一卡2卡4卡精品| 韩国精品一区二区三区 | 一级毛片黄色毛片免费观看视频| 热99久久久久精品小说推荐| 欧美精品高潮呻吟av久久| 免费黄网站久久成人精品| 99热全是精品| 国内精品宾馆在线| 亚洲 欧美一区二区三区| 久久久国产一区二区| 色吧在线观看| 欧美日韩精品成人综合77777| 亚洲精品自拍成人| 日本-黄色视频高清免费观看| 九草在线视频观看| 欧美97在线视频| 国产一区有黄有色的免费视频| 丝瓜视频免费看黄片| 国产精品久久久久久久久免| 久久精品夜色国产| 在线 av 中文字幕| 国产成人a∨麻豆精品| 亚洲 欧美一区二区三区| 伦理电影免费视频| 日本猛色少妇xxxxx猛交久久| 有码 亚洲区| 少妇的丰满在线观看| 久久久久久伊人网av| 久久精品久久久久久噜噜老黄| 日本午夜av视频| 69精品国产乱码久久久| 成人手机av| 91国产中文字幕| 欧美精品一区二区免费开放| 欧美亚洲 丝袜 人妻 在线| 美女大奶头黄色视频| 欧美精品一区二区大全| 国产免费一级a男人的天堂| 亚洲国产精品一区二区三区在线| 国产免费一级a男人的天堂| 国产女主播在线喷水免费视频网站| 欧美成人午夜免费资源| 中文字幕免费在线视频6| 蜜桃国产av成人99| 妹子高潮喷水视频| 在现免费观看毛片| 久久午夜福利片| 好男人视频免费观看在线| 丰满少妇做爰视频| 丝袜美足系列| 精品少妇内射三级| 91精品国产国语对白视频| 丝袜美足系列| 久久国内精品自在自线图片| 亚洲av免费高清在线观看| 春色校园在线视频观看| 日本爱情动作片www.在线观看| 精品第一国产精品| 校园人妻丝袜中文字幕| 久久 成人 亚洲| 啦啦啦视频在线资源免费观看| 9191精品国产免费久久| 日韩免费高清中文字幕av| 国产免费一区二区三区四区乱码| 热re99久久国产66热| 久久 成人 亚洲| videossex国产| 热99国产精品久久久久久7| 色婷婷久久久亚洲欧美| 国产片特级美女逼逼视频| 成年人午夜在线观看视频| 99久国产av精品国产电影| 91精品国产国语对白视频| 久久午夜福利片| 国内精品宾馆在线| 国产精品99久久99久久久不卡 | 欧美日韩一区二区视频在线观看视频在线| 99视频精品全部免费 在线| 日韩熟女老妇一区二区性免费视频| 90打野战视频偷拍视频| 国产高清不卡午夜福利| 一级毛片黄色毛片免费观看视频| 精品福利永久在线观看| 91午夜精品亚洲一区二区三区| 男女啪啪激烈高潮av片| 最新的欧美精品一区二区| 成人黄色视频免费在线看| 中文字幕人妻熟女乱码| 亚洲熟女精品中文字幕| 日本午夜av视频| 国产一区二区在线观看av| 亚洲精品色激情综合| 精品午夜福利在线看| 国产精品三级大全| 黄片播放在线免费| 亚洲精品久久午夜乱码| 少妇猛男粗大的猛烈进出视频| 看免费av毛片| 美女福利国产在线| 99久久精品国产国产毛片| 亚洲av免费高清在线观看| 一二三四在线观看免费中文在 | 国产色爽女视频免费观看| 成年女人在线观看亚洲视频| 午夜免费观看性视频| 蜜桃在线观看..| 亚洲国产日韩一区二区| 亚洲成色77777| 亚洲av日韩在线播放| a级毛片在线看网站| a级毛片在线看网站| 免费人妻精品一区二区三区视频| 国产一区二区三区av在线| 久久精品国产自在天天线| www.色视频.com| 一本—道久久a久久精品蜜桃钙片| 色哟哟·www| 天美传媒精品一区二区| 亚洲婷婷狠狠爱综合网| 日韩欧美精品免费久久| 国产精品久久久久成人av| 国产精品一区www在线观看| 欧美 日韩 精品 国产| 欧美精品高潮呻吟av久久| 午夜免费观看性视频| 99re6热这里在线精品视频| 最近中文字幕2019免费版| 国产精品久久久久久久电影| 亚洲av欧美aⅴ国产| 欧美精品一区二区免费开放| 99热网站在线观看| 精品一区在线观看国产| 国产综合精华液| 黄色毛片三级朝国网站| 亚洲av免费高清在线观看| 亚洲av综合色区一区| tube8黄色片| 国产成人午夜福利电影在线观看| 女人久久www免费人成看片| 久久国产精品大桥未久av| 22中文网久久字幕| 亚洲av欧美aⅴ国产| 制服丝袜香蕉在线| 国产精品无大码| 桃花免费在线播放| www.熟女人妻精品国产 | 精品亚洲成国产av| 久久久久久久亚洲中文字幕| 亚洲精品一二三| 久久女婷五月综合色啪小说| 性高湖久久久久久久久免费观看| 国产熟女欧美一区二区| 香蕉丝袜av| 亚洲欧美成人精品一区二区| 国产一区二区激情短视频 | 97人妻天天添夜夜摸| 蜜桃国产av成人99| 久久国产精品男人的天堂亚洲 | 大片电影免费在线观看免费| 日韩成人av中文字幕在线观看| 又黄又粗又硬又大视频| www日本在线高清视频| a级毛色黄片| 国产成人欧美| 晚上一个人看的免费电影| 一个人免费看片子| 国产在线视频一区二区| 午夜福利网站1000一区二区三区| 久久99精品国语久久久| 一区二区av电影网| 一级片'在线观看视频| 2018国产大陆天天弄谢| 免费大片18禁| 老司机影院成人| 久久久久久人妻| 青春草视频在线免费观看| 伦精品一区二区三区| 黄色一级大片看看| 天天操日日干夜夜撸| 男女边摸边吃奶| 日韩大片免费观看网站| av天堂久久9| 毛片一级片免费看久久久久| 男人爽女人下面视频在线观看| 超碰97精品在线观看| 另类精品久久| 丝瓜视频免费看黄片| 人妻一区二区av| 免费观看av网站的网址| 国产一区二区三区av在线| 亚洲综合色网址| 啦啦啦中文免费视频观看日本| 欧美人与性动交α欧美软件 | 天堂中文最新版在线下载| 黄片无遮挡物在线观看| 色婷婷av一区二区三区视频| 男女高潮啪啪啪动态图| 一本大道久久a久久精品| 99热全是精品| 欧美97在线视频| 国产成人精品在线电影| 只有这里有精品99| 亚洲一级一片aⅴ在线观看| 久久精品久久精品一区二区三区| 一本大道久久a久久精品| 亚洲国产成人一精品久久久| 亚洲人成网站在线观看播放| 免费av不卡在线播放| 18在线观看网站| 日本vs欧美在线观看视频| 天堂8中文在线网| 国产片特级美女逼逼视频| 飞空精品影院首页| 精品一品国产午夜福利视频| 久久久久精品人妻al黑| 最新的欧美精品一区二区| 久久久久视频综合| 99久久综合免费| 啦啦啦啦在线视频资源| 午夜福利乱码中文字幕| 欧美日韩综合久久久久久| 最近最新中文字幕大全免费视频 | 日本色播在线视频| 久久久国产一区二区| 亚洲精品国产av成人精品| 晚上一个人看的免费电影| 久久精品aⅴ一区二区三区四区 | 久久久久久久亚洲中文字幕| 精品国产一区二区三区久久久樱花| 91aial.com中文字幕在线观看| 大码成人一级视频| 国产成人免费无遮挡视频| 日韩一区二区三区影片| 日韩电影二区| 一级毛片黄色毛片免费观看视频| 亚洲天堂av无毛| 成人影院久久| 国产免费一区二区三区四区乱码| 天天躁夜夜躁狠狠躁躁| 蜜桃在线观看..| 久久久亚洲精品成人影院| 18禁裸乳无遮挡动漫免费视频| 毛片一级片免费看久久久久| 亚洲成人一二三区av| 秋霞伦理黄片| 亚洲av日韩在线播放| av片东京热男人的天堂| 亚洲色图综合在线观看| 日韩 亚洲 欧美在线| 日本欧美视频一区| 最近的中文字幕免费完整| 亚洲少妇的诱惑av| 成人亚洲精品一区在线观看| 18在线观看网站| www.熟女人妻精品国产 | 国产免费现黄频在线看| 一区二区三区精品91| 丰满乱子伦码专区| 大香蕉97超碰在线| 日韩制服丝袜自拍偷拍| 深夜精品福利| 日韩一本色道免费dvd| 少妇的逼好多水| av女优亚洲男人天堂| 免费黄频网站在线观看国产| 90打野战视频偷拍视频| 亚洲av成人精品一二三区| 久久97久久精品| 麻豆精品久久久久久蜜桃| 国产午夜精品一二区理论片| 在线天堂最新版资源| 国产成人91sexporn| 99久久综合免费| 久久久久久久国产电影| 亚洲欧美成人精品一区二区| 男女边摸边吃奶| 一级毛片电影观看| 久久精品夜色国产| 午夜免费观看性视频| 精品福利永久在线观看| 日本猛色少妇xxxxx猛交久久| 日韩 亚洲 欧美在线| 少妇被粗大的猛进出69影院 | 菩萨蛮人人尽说江南好唐韦庄| 精品卡一卡二卡四卡免费| 99久久人妻综合| 国产欧美亚洲国产| 香蕉国产在线看| 久久99热这里只频精品6学生| 欧美 日韩 精品 国产| 国产在线免费精品| 99热国产这里只有精品6| 国产成人av激情在线播放| 黄网站色视频无遮挡免费观看| 日韩一区二区视频免费看| 22中文网久久字幕| 精品一区在线观看国产| 韩国av在线不卡| 国产毛片在线视频| 91aial.com中文字幕在线观看| 人人妻人人添人人爽欧美一区卜| 亚洲五月色婷婷综合| av视频免费观看在线观看| 亚洲精品乱码久久久久久按摩| 中国国产av一级| 51国产日韩欧美| 日韩精品免费视频一区二区三区 | 色吧在线观看| 成人无遮挡网站| 亚洲伊人色综图| 免费观看av网站的网址| 亚洲综合色网址| 欧美精品亚洲一区二区| 蜜桃在线观看..| 香蕉国产在线看| 中文字幕av电影在线播放| 晚上一个人看的免费电影| 免费av不卡在线播放| 女人精品久久久久毛片| 国产在视频线精品| 国产亚洲av片在线观看秒播厂| 国产成人精品在线电影| 日韩大片免费观看网站| 女性被躁到高潮视频| 精品一品国产午夜福利视频| 亚洲精品视频女| 国产视频首页在线观看| 国产无遮挡羞羞视频在线观看| 免费看光身美女| 日本av免费视频播放| 亚洲精品av麻豆狂野| 熟妇人妻不卡中文字幕| 亚洲人成77777在线视频| 亚洲成国产人片在线观看| 曰老女人黄片| 日韩在线高清观看一区二区三区| 香蕉国产在线看| 午夜福利视频精品| 老女人水多毛片| 交换朋友夫妻互换小说| 亚洲欧洲日产国产| 另类精品久久| videosex国产| 日韩三级伦理在线观看| 熟女人妻精品中文字幕| 欧美+日韩+精品| 国产国语露脸激情在线看| 国产亚洲精品久久久com| 免费少妇av软件| 午夜日本视频在线| 欧美日韩视频高清一区二区三区二| 久久久久久久国产电影| 九九在线视频观看精品| 精品久久蜜臀av无| 亚洲人与动物交配视频| 亚洲精品美女久久久久99蜜臀 | 国产国语露脸激情在线看| 最后的刺客免费高清国语| 青春草视频在线免费观看| freevideosex欧美| 人人澡人人妻人| 女的被弄到高潮叫床怎么办| 国产男女超爽视频在线观看| 亚洲精品乱久久久久久| 亚洲欧美一区二区三区国产| 国产成人一区二区在线| 观看av在线不卡| 性色avwww在线观看| 欧美97在线视频| 2022亚洲国产成人精品| 一区二区三区四区激情视频| 久久久国产欧美日韩av| 国产女主播在线喷水免费视频网站| 青春草亚洲视频在线观看| 内地一区二区视频在线| 高清黄色对白视频在线免费看| 亚洲欧洲日产国产| 免费看av在线观看网站| 黑丝袜美女国产一区| 啦啦啦在线观看免费高清www| 99热这里只有是精品在线观看| av有码第一页| 欧美精品国产亚洲| 男女午夜视频在线观看 | 热re99久久国产66热| 日韩电影二区| 香蕉国产在线看| 好男人视频免费观看在线| 成年女人在线观看亚洲视频| 成人手机av| 精品久久久精品久久久| 热99久久久久精品小说推荐| 99香蕉大伊视频| 人妻人人澡人人爽人人| 亚洲三级黄色毛片| 超色免费av| 最近手机中文字幕大全| 热99久久久久精品小说推荐| 国产熟女午夜一区二区三区| 大片免费播放器 马上看| 亚洲丝袜综合中文字幕| 日本与韩国留学比较| 毛片一级片免费看久久久久| 亚洲国产看品久久| 日本av免费视频播放| 欧美人与性动交α欧美软件 | 有码 亚洲区| 国产欧美日韩一区二区三区在线| 高清黄色对白视频在线免费看| 精品第一国产精品| 精品亚洲乱码少妇综合久久| 男的添女的下面高潮视频| 亚洲精品av麻豆狂野| 天天影视国产精品| 日本色播在线视频| 精品一区二区三区视频在线| 精品少妇黑人巨大在线播放| av免费在线看不卡| 国产精品一二三区在线看| 国产av精品麻豆| www.熟女人妻精品国产 | 久久99蜜桃精品久久| 欧美日韩国产mv在线观看视频| 国产av一区二区精品久久| 日韩制服骚丝袜av| 欧美精品亚洲一区二区| 成人亚洲精品一区在线观看| 亚洲精品国产色婷婷电影| 久久久久久人妻| 欧美日韩成人在线一区二区| 午夜91福利影院| 欧美97在线视频| 最黄视频免费看| 最近手机中文字幕大全| 久久精品国产亚洲av涩爱| 自拍欧美九色日韩亚洲蝌蚪91| 99热这里只有是精品在线观看| 99视频精品全部免费 在线| 制服人妻中文乱码| av线在线观看网站| 国产精品偷伦视频观看了| 大片免费播放器 马上看| 亚洲人成网站在线观看播放| 日日撸夜夜添| 欧美丝袜亚洲另类| 永久网站在线| 免费观看无遮挡的男女| 国产av一区二区精品久久| 国产成人免费观看mmmm| 在线观看www视频免费| 久久99一区二区三区| 色视频在线一区二区三区| 男人添女人高潮全过程视频| 下体分泌物呈黄色| 欧美97在线视频| 日本wwww免费看| 激情视频va一区二区三区| 国产免费一级a男人的天堂| 亚洲欧美日韩卡通动漫| 考比视频在线观看| 一级毛片 在线播放| 美女视频免费永久观看网站| 插逼视频在线观看| 搡老乐熟女国产| 日本欧美视频一区| 波多野结衣一区麻豆| 欧美日韩av久久| 久久久久久久久久久免费av| 天堂8中文在线网| 久久精品熟女亚洲av麻豆精品| 少妇人妻 视频| 亚洲欧美日韩另类电影网站| 日本免费在线观看一区| 久久青草综合色| 久久这里有精品视频免费| www.色视频.com| 纯流量卡能插随身wifi吗| 男女下面插进去视频免费观看 | 国产毛片在线视频| 国产男人的电影天堂91| 飞空精品影院首页| 日韩一本色道免费dvd| videossex国产| 亚洲精品美女久久av网站| 国产男女内射视频| 亚洲精品一区蜜桃| 中文欧美无线码| 99热6这里只有精品| 日本-黄色视频高清免费观看| 久久韩国三级中文字幕| 人妻一区二区av| 亚洲内射少妇av| 少妇精品久久久久久久| 国产精品无大码| 99热6这里只有精品| 亚洲一码二码三码区别大吗| 精品卡一卡二卡四卡免费| 亚洲久久久国产精品| 九色成人免费人妻av| 男女午夜视频在线观看 | 大香蕉97超碰在线| 久热这里只有精品99| 亚洲国产精品专区欧美| 国产精品久久久久久精品电影小说| 宅男免费午夜| 一级片'在线观看视频| 精品亚洲成a人片在线观看| 国产国语露脸激情在线看| 成人二区视频| 最新的欧美精品一区二区| 国产精品一区www在线观看| 男女免费视频国产| 亚洲精品国产av成人精品| 亚洲中文av在线| 黄色怎么调成土黄色| 人人妻人人爽人人添夜夜欢视频| 国产亚洲午夜精品一区二区久久| 国产一区有黄有色的免费视频| 少妇人妻久久综合中文| 欧美成人精品欧美一级黄| 国产极品天堂在线| 18禁观看日本| 国产成人91sexporn| 一个人免费看片子| 午夜激情久久久久久久| 国产av一区二区精品久久| 亚洲国产精品一区三区| 亚洲在久久综合| 国产1区2区3区精品| 色网站视频免费| 热99国产精品久久久久久7| 日韩人妻精品一区2区三区| 久久久精品区二区三区| 精品国产乱码久久久久久小说| 伊人久久国产一区二区| 菩萨蛮人人尽说江南好唐韦庄| tube8黄色片| 亚洲欧洲日产国产| 亚洲,欧美精品.| 一区二区三区四区激情视频| 成人亚洲欧美一区二区av| www日本在线高清视频| 日韩av在线免费看完整版不卡| 丰满饥渴人妻一区二区三| 一区在线观看完整版| 男人舔女人的私密视频| 成人国产av品久久久| 亚洲激情五月婷婷啪啪| 欧美丝袜亚洲另类| 最新中文字幕久久久久| 丁香六月天网| 一级a做视频免费观看| 国产毛片在线视频| 男男h啪啪无遮挡| 成年美女黄网站色视频大全免费| 最近2019中文字幕mv第一页| 两性夫妻黄色片 | 在线天堂最新版资源| 国产精品一区www在线观看| 2022亚洲国产成人精品| 久久人人97超碰香蕉20202| 久久久久视频综合| 日韩电影二区| 久久人人爽人人爽人人片va| 国产成人精品一,二区| 国产麻豆69| 十八禁高潮呻吟视频| 七月丁香在线播放| 国产亚洲精品第一综合不卡 | 韩国高清视频一区二区三区| 两性夫妻黄色片 | 国产又色又爽无遮挡免| 国产精品欧美亚洲77777| 欧美变态另类bdsm刘玥| 亚洲综合色网址| 日韩精品免费视频一区二区三区 | 性色avwww在线观看| 香蕉丝袜av| 久久国内精品自在自线图片| 欧美xxⅹ黑人| 久久久久网色| 亚洲精品aⅴ在线观看| 欧美xxxx性猛交bbbb| 汤姆久久久久久久影院中文字幕| 最近的中文字幕免费完整| 这个男人来自地球电影免费观看 | 久久免费观看电影| 亚洲欧美日韩另类电影网站| 在线观看免费日韩欧美大片| 婷婷成人精品国产| 免费播放大片免费观看视频在线观看| 乱人伦中国视频| 一本—道久久a久久精品蜜桃钙片| 国产日韩欧美在线精品| 国产男女内射视频| 国产成人午夜福利电影在线观看| 黑人欧美特级aaaaaa片| 99国产精品免费福利视频| 黄片播放在线免费| 黄色一级大片看看|