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

    仿自然魚道中卵石墻對池室水力特性改善效果

    2017-09-15 06:17:36李廣寧孫雙科柳海濤趙桂俠鄭鐵剛
    關(guān)鍵詞:魚道不透水隔墻

    李廣寧,孫雙科※,柳海濤,張 超,趙桂俠,鄭鐵剛

    仿自然魚道中卵石墻對池室水力特性改善效果

    李廣寧1,孫雙科1※,柳海濤1,張 超1,趙桂俠2,鄭鐵剛1

    (1. 中國水利水電科學(xué)研究院流域水循環(huán)模擬與調(diào)控國家重點(diǎn)實(shí)驗(yàn)室,北京 100038;2. 天津大學(xué)水利工程仿真與安全國家重點(diǎn)實(shí)驗(yàn)室,天津 300072)

    仿自然魚道是重要的魚道布置方式之一,其特點(diǎn)是利用當(dāng)?shù)睾哟膊牧蠘?gòu)建透水性卵石隔墻,以便更好地模擬出天然河流流態(tài)。為了定量分析透水性卵石墻對魚道水力特性的改善效果,該文針對采用卵石隔墻與不透水隔墻的魚道水力特性進(jìn)行對比研究。結(jié)果表明:采用卵石墻的仿自然魚道,池室內(nèi)流態(tài)豐富,過魚口處流速分布差異明顯,表底流速差達(dá)到0.43~0.61 m/s為多種魚類上溯提供了可能。池室內(nèi)主流區(qū)分布較寬,紊動能較大達(dá)0.06~0.12 m2/s2,同時回流區(qū)減弱,有利于魚類找到主流,實(shí)現(xiàn)上溯。卵石隔墻的不足之處是魚道耗水流量相對較大達(dá)2.86 m3/s,在工程設(shè)計(jì)中需要綜合考慮卵石墻透水特性與耗水流量之間的平衡關(guān)系,找到適宜的布置方案。該研究可為仿自然魚道的水力設(shè)計(jì)提供參考。

    魚道;流速;流量;紊動能;耗水流量;數(shù)值模擬

    李廣寧,孫雙科,柳海濤,張 超,趙桂俠,鄭鐵剛. 仿自然魚道中卵石墻對池室水力特性改善效果[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(15):184-189. doi:10.11975/j.issn.1002-6819.2017.15.024 http://www.tcsae.org

    Li Guangning, Sun Shuangke, Liu Haitao, Zhang Chao, Zhao Guixia, Zheng Tiegang. Improving effect of hydraulic characteristics of nature-like fishway with pools and cobblestone weirs[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(15): 184-189. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2017.15.024 http://www.tcsae.org

    0 引 言

    魚道是幫助魚類克服障礙物(如各種壩、堰、水閘)實(shí)現(xiàn)洄游上溯的設(shè)施[1]。魚道分為技術(shù)型魚道和仿自然魚道,傳統(tǒng)的技術(shù)型魚道(如丹尼爾式魚道、豎縫式魚道等)的結(jié)構(gòu)與流態(tài)相對單一,一般僅適合特定過魚對象[2-6],而在實(shí)際工程中,同一條河流中需要保護(hù)的魚類種類往往是多樣的,且魚類體型相差懸殊、生活習(xí)性和克流能力均不相同[7-9]。與傳統(tǒng)技術(shù)型魚道相比,仿自然魚道構(gòu)建的水流流態(tài)更接近于魚類熟悉的天然狀態(tài),對魚類往往具有更廣的適用性和更高的過魚效率[10-11]。

    仿自然魚道常采用水池-淺灘型式[12],其特點(diǎn)是采用梯級低堰將魚道分隔成一系列的淺灘和深池,淺灘處水深較淺,流速較大,而深池內(nèi)水流緩慢。池室型仿自然魚道[13]是由這種魚道發(fā)展而來,其特點(diǎn)是采用隔墻將魚道分隔成多級水池,墻體通常為卵石墻(數(shù)塊卵石堆疊而成的下寬上窄的墻體),卵石間縫這種分隔方式建造的魚道,不但能夠就地取材,利用漂石與天然河道床沙質(zhì),其水流還兼具水池式魚道與技術(shù)型魚道的雙重特性。

    目前,國外針對仿自然魚道的水力特性已開展了一些研究[14-15],國內(nèi)的相關(guān)研究尚處于起步階段[16]。相比而言,國外學(xué)者比較重視研究各種塊體障礙物對水流結(jié)構(gòu)的影響,如Acharya等[17]針對球狀障礙物、圓柱體和方柱體對水流結(jié)構(gòu)的影響進(jìn)行了模型試驗(yàn),提出為了達(dá)到相似水流條件(弗洛德數(shù)Fr),球體的排列應(yīng)該比方柱體的排列略微緊密一些,但是球體障礙物尾流區(qū)水面連接更加平順。加拿大Alberta大學(xué)的Bretón等[18]研究了蠻石斜坡型魚道中,蠻石間的流速與紊動強(qiáng)度等水力學(xué)參數(shù),認(rèn)為蠻石布置間距應(yīng)該足夠小,以遏制水流的加速,蠻石后尾流區(qū)的紊動能在一定的范圍內(nèi)變化(0.05~0.25 m2/s2)。Baki等[19-20]通過試驗(yàn)研究也得出了類似結(jié)論。Lacey等除了進(jìn)行蠻石周圍的紊流結(jié)構(gòu)的試驗(yàn)研究[21],還實(shí)測了自然河流中典型蠻石周邊的水流流場及紊動能等水力參數(shù),發(fā)現(xiàn)在天然河道中蠻石后的尾流區(qū)內(nèi)往往會形成較強(qiáng)烈的紊動[22]。在數(shù)值模擬方面,Tran等[23]利用淺水方程對蠻石斜坡型魚道的水流特性進(jìn)行了模擬,并進(jìn)行了試驗(yàn)驗(yàn)證。何雨艨等[24]采用三維湍流模型對不同工況的蠻石斜坡型魚道進(jìn)行了數(shù)值模擬與試驗(yàn)驗(yàn)證,認(rèn)為蠻石坎對縫排列比錯縫排列具有更加合理的魚道水力學(xué)特性。

    本文采用Flow-3D軟件,分別針對池室型仿自然魚道內(nèi),采用卵石墻與不透水墻2種型式進(jìn)行三維數(shù)值模擬,對兩者的流態(tài)、流速、紊動能以及過流流量進(jìn)行對比分析,以明晰兩者水力特征的區(qū)別,為工程實(shí)際應(yīng)用提供參考。

    1 數(shù)學(xué)模型及驗(yàn)證

    1.1 基本方程及求解方法

    數(shù)值模擬采用standard k-ε湍流模型[25],該模型基本方程包括水流連續(xù)方程、動量方程、紊動能及其耗散率方程。

    連續(xù)方程:

    動量方程:

    紊動能k方程:

    紊動能耗散率ε方程:

    式中t為時間,s;u、v、w 為x、y、z方向的時均速度分量為坐標(biāo)分量,m/s;ρ為密度,kg/m3;μ、μt分別為黏性系數(shù)和紊動黏性系數(shù),m2/s;p為時均壓強(qiáng),Pa;k為紊動能,m2/s2;ε為紊動能耗散率,m2/s3;Gk為紊動能產(chǎn)生項(xiàng),各項(xiàng)紊流常數(shù)取值為Cμ=0.09,σε=1.3,σk=1.0,C1ε=1.44,C2ε=1.92。

    液體自由表面采用VOF方法進(jìn)行計(jì)算[26],通過計(jì)算水和氣的體積分?jǐn)?shù)來表征物體的形態(tài),aa表示氣體的體積分?jǐn)?shù),aw表示水的體積分?jǐn)?shù),其控制方程

    模型求解采用有限體積法,二階迎風(fēng)格式,壓力-速度耦合采用壓力校正法,離散方程的求解采用GMRES法,時間差分采用全隱格式。

    1.2 水體中障礙物的計(jì)算方法

    計(jì)算模型中對于魚道內(nèi)具有一定透水性的卵石隔墻,采用FAVOR(fractional area/volume ratios)方法[27]進(jìn)行定義。其原理是利用控制體積(網(wǎng)格)的6個面定義出可以有通量通過的局部面積與可自由通過體積的比率,若無可以讓流體通過的面積則交界面不會有對流的通量,即FAVOR技術(shù)提供每個控制體積對于流體與障礙物交界面積的決定方法,網(wǎng)格與物體幾何形狀相互獨(dú)立。盡管模型復(fù)雜,利用FAVOR技術(shù)也可以精確的表示其經(jīng)過網(wǎng)格計(jì)算后的外觀。引入VF=開放體積/單元體積,為網(wǎng)格可通過的通量比率,即通量體積比(fractional volume ratio),AF=開放面積/單元邊面積,為各面可通過的通量比率,即通量面積比(fractional area ratio),網(wǎng)格內(nèi)包含的灰色部分不會有通量通過,定義為障礙物,其余部分通量可通過,如圖1所示,即以網(wǎng)格數(shù)目尺寸來表現(xiàn)原有體型的外觀,并判斷目前網(wǎng)格密度能否精確反應(yīng)模型外觀,避免因網(wǎng)格密度不足造成結(jié)果不精確。

    圖1 FAVOR技術(shù)障礙物定義示意圖Fig.1 Schematic diagram of obstacle in FAVOR technology

    1.3 模型體型及邊界條件

    本文以某擬建魚道工程為基礎(chǔ)開展研究,該魚道工程主要過魚對象為短須裂腹魚,體長0.25~0.36 m,持續(xù)泳速0.59~1.17 m/s,突進(jìn)泳速1.17~2.50 m/s。

    魚道數(shù)值模型的模擬范圍包括水平段和池室段2個部分,水平段長15 m,池室段長40 m,池室段內(nèi)共模擬5級魚池。魚道橫斷面為梯形,底寬5 m,頂寬12.5 m,邊坡系數(shù)1.5,底坡1%。過魚口位于河槽左側(cè),寬度1.0 m,隔墻采用交錯布置型式,開口處出流方向與主槽呈45o夾角。研究中針對池室隔墻為卵石墻和不透水墻2種型式,分別建立數(shù)學(xué)模型。其中,透水卵石墻的卵石粒徑為15~30 cm (當(dāng)?shù)伫Z卵石粒徑范圍),共堆疊3層,底部大,頂端小,呈金字塔狀,依靠自重保持結(jié)構(gòu)穩(wěn)定,總高1.8 m,底厚1.5 m,間距6.5 m;不透水墻體型是將卵石墻替換為的矩形隔墻,其高度、厚度和位置均保持不變。上述布置體型見圖2a與2b。

    計(jì)算模型采用嵌套加密的矩形網(wǎng)格,卵石墻附近網(wǎng)格尺寸0.05 m,其余部位網(wǎng)格尺寸0.1~0.5 m,網(wǎng)格總數(shù)約176萬。計(jì)算網(wǎng)格劃分情況如圖2c所示。

    計(jì)算過程中,模型入口和出口邊界均設(shè)置為固定水深1.4 m;卵石墻與固壁均采用無滑移條件;卵石表面糙率取值,nl=0.015;魚道邊壁及池底糙率參照魚道實(shí)際情況進(jìn)行取值,np=0.03;卵石墻的孔隙率約12%,為縱剖面孔隙所占過流面積的比例。

    圖2 具有不同隔水墻的魚道模型及其網(wǎng)格劃分Fig.2 Fishways with different weir types and computational mesh

    1.4 模型驗(yàn)證

    卵石墻仿自然魚道模型按重力相似準(zhǔn)則設(shè)計(jì),模型比尺λL=5,其結(jié)構(gòu)布置形式與數(shù)值計(jì)算體型一致。仿自然魚道原型底部糙率為0.03,由糙率比尺=λ1/6=1.308可知,模型基礎(chǔ)糙率應(yīng)為n=0.023。模型試驗(yàn)中,魚道基礎(chǔ)底面采用混凝土抹面,糙率為n=0.014,因此需要通過加糙試驗(yàn),達(dá)到目標(biāo)糙率。

    由明渠均勻流公式可知,明渠糙率n可以表達(dá)為式中,A為體型明渠過水面積,m2;Q為流量,m3/s;R為水力半徑,m;J為渠底坡降,%。根據(jù)模型橫斷面尺寸,可知過水面積A=h(b+mh)=h(1+1.5h);對應(yīng)濕周;對應(yīng)的水力半徑為。由此可知,只要知道上游來流量Q和均勻流水深h,即可求解渠道糙率值n。試驗(yàn)中采用了正方形梅花加糙方法,顆粒直徑2.5 cm,顆粒之間間距20 cm,測得模型糙率在0.022~0.023之間,滿足試驗(yàn)精度要求。

    試驗(yàn)中通過調(diào)節(jié)進(jìn)口來流流量及出口尾門的開度,確保魚道上下游水深均為1.4 m。試驗(yàn)測點(diǎn)布置如圖3a所示,其中V1~V6為速度測點(diǎn),位于各級池室過魚口處,A1~E2均為水面線測點(diǎn),位于隔墻上下游滯水區(qū)。圖3b與3c為流速和水面線的試驗(yàn)與計(jì)算結(jié)果對比,流速誤差小于0.1 m/s,水面線誤差小于0.05 m,表明兩者吻合良好。

    圖3 測點(diǎn)布置、流速、水面線計(jì)算值和試驗(yàn)值Fig.3 Arrangement of measuring points and comparison of simulation and test velocity values and water surface elevation

    2 結(jié)果與分析

    2.1 流態(tài)和流速

    根據(jù)魚道設(shè)計(jì)基本原則,魚道池室內(nèi)最大流速一般出現(xiàn)在過魚口附近,其大小決定了魚類是否能夠成功通過該水池,其值以不超過魚類的突進(jìn)流速為宜[28-29]。過魚口處的流速分布在垂向上應(yīng)存在一定差異,可為不同游泳能力的魚類的提供上溯通道。主流區(qū)分布范圍不宜過于集中,兩側(cè)不宜形成大范圍的回流區(qū),以便于休息區(qū)的魚類能夠迅速找到主流,繼續(xù)上溯,池室內(nèi)休息區(qū)的水流具有豐富的流態(tài),以接近自然河流形態(tài)為宜[11],可為魚類提供休息、覓食場所與上溯通道。上述流場特性的影響因素較多,如隔墻類型,底坡變化,運(yùn)行水深等等。

    圖4a和4b為卵石墻池室內(nèi)表層與底層流場計(jì)算結(jié)果。分析表明,其主流流速呈現(xiàn)“表層小,底層大”的分布規(guī)律。卵石墻開口處流速分布具有“表層小,底層大,兩側(cè)小,中間大”的特征,如流速測點(diǎn)V3的表底流速在1.33~1.94 m/s,差值為0.61 m/s,測點(diǎn)V4的表底流速在1.43~1.86 m/s,差值為0.43 m/s。流速量值垂向分布范圍較寬,且介于魚類突進(jìn)泳速1.17~2.5 m/s范圍內(nèi),滿足目標(biāo)魚類上溯需求。主流兩側(cè)未出現(xiàn)大范圍回流區(qū),有利于魚類重新感應(yīng)主流,繼續(xù)上溯。

    圖4c和4d為不透水墻池室內(nèi)表層與底層流場計(jì)算結(jié)果,相比而言,其流態(tài)分布更接近于技術(shù)型魚道。首先,池室開口處流速分布上下基本一致。如流速測點(diǎn)V3的表底流速在1.56~1.64 m/s,差值為0.08 m/s,測點(diǎn)V4的表底流速在1.63~1.69 m/s,差值為0.06 m/s。流速分布較為單一。其次,主流區(qū)較為集中,休息區(qū)內(nèi)分別形成了兩個回流區(qū),主流區(qū)與回流區(qū)之間存在流速接近于0的流速帶,不利于魚類發(fā)現(xiàn)主流,實(shí)現(xiàn)上溯。

    圖4 兩種體型第3級水池(池3)流速分布Fig.4 Velocity distributions of third stage pool (Pool 3) in two schemes

    2.2 紊動能

    水流紊動能(turbulence kinetic energy, TKE)是反映水流紊動強(qiáng)度的指標(biāo),魚道中水流保持適當(dāng)?shù)奈蓜訌?qiáng)度,一方面有利于魚類感應(yīng)主流位置,同時也有利于增加水體含氧量,維持水質(zhì)指標(biāo)[30-31]。圖5為卵石隔墻與不透水條件下,池室內(nèi)表層與底層水體紊動能分布。計(jì)算結(jié)果分析表明,采用不透水墻的仿自然魚道,池室內(nèi)主流區(qū)紊動能僅為0.01~0.06 m2/s2;采用卵石墻的仿自然魚道,受卵石不規(guī)則外形和縫隙間水流的沖擊影響,池室內(nèi)主流區(qū)紊動能可達(dá)0.06~0.12 m2/s2,便于魚類感應(yīng)到主流位置,實(shí)現(xiàn)上溯。卵石墻體型下,水體紊動能高于不透水墻池室體型,其原因主要是水體受到卵石不規(guī)則外形的干擾,以及卵石縫隙間水流的沖擊所致。池室內(nèi)表層水體紊動能高于底層,如卵石墻體型的表層紊動能最大值為0.12 m2/s2,底層最大值為0.06 m2/s2。究其原因在于,魚道池室間存在水面落差,水流跌落對水體表面沖擊所致。總體而言,采用卵石墻的魚道體型,池室內(nèi)紊動能較不透水墻魚道體型有所增加。

    圖5 兩種體型第3級水池紊動能分布Fig.5 Turbulent kinetic energy distributions of third stage pool in two schemes

    2.3 過流流量

    表1為兩種魚道體型各個池室的水位變化與耗水流量。在魚道內(nèi)水位保持相同(對應(yīng)各級池室內(nèi)的水深約1.4 m)的條件下,不透水墻魚道體型的耗水流量為Q=1.39 m3/s,而卵石墻魚道體型的耗水流量為Q=2.86 m3/s,增大了一倍左右。究其原因在于,隔墻內(nèi)的縫隙增加了魚道整體的過流面積,在相同運(yùn)行水深下,卵石墻魚道體型較不透水隔墻魚道需要消耗更多的水量。

    表1 兩種體型各個池室水位及耗水流量Table 1 Pool water level and consumption flow in two schemes

    3 結(jié) 論

    本文采用Flow-3D建立了不同隔墻類型的魚道三維數(shù)學(xué)模型,結(jié)合物理模型試驗(yàn)驗(yàn)證,針對卵石墻仿自然魚道和不透水隔墻魚道中的水力特性進(jìn)行了對比分析,研究表明:

    1)本文建立的三維數(shù)學(xué)模型可以較好地描述透水性卵石墻對魚道內(nèi)流場的影響,數(shù)學(xué)模型計(jì)算結(jié)果與物理模型試驗(yàn)結(jié)果吻合良好,為仿自然魚道的研究提供了可靠的手段。

    2)采用不透水隔墻的魚道,池室開口處流速分布較為單一,表底流速差僅為0.06~0.08 m/s;而采用卵石墻的仿自然魚道,過魚口處流速分布表現(xiàn)出“表層小,底層大,兩側(cè)小,中間大”的特征,且流速分布垂向差異明顯,表底流速差達(dá)到0.43~0.61 m/s,為多種魚類上溯提供了可能。

    3)采用不透水隔墻的魚道,各池室內(nèi)主流較為集中,兩側(cè)形成較大的回流區(qū),兩者之間存在一個低流速帶,不利于魚類回到主流,繼續(xù)上溯;采用卵石墻后,水體受到卵石縫隙間水流的沖擊,流態(tài)豐富,回流區(qū)減弱,有利于魚類重新找到主流,實(shí)現(xiàn)上溯。

    4)采用不透水墻的仿自然魚道,池室內(nèi)主流區(qū)紊動能僅為0.01~0.06 m2/s2;采用卵石墻的仿自然魚道,受卵石不規(guī)則外形和縫隙間水流的沖擊影響,池室內(nèi)主流區(qū)紊動能可達(dá)0.06~0.12 m2/s2,便于魚類感應(yīng)到主流位置,實(shí)現(xiàn)上溯。

    5)采用透水性卵石墻后,魚道耗水流量較之不透水墻體型,增大1倍左右。建議在實(shí)際工程設(shè)計(jì)中,通過改變卵石墻透水特性,尋求水流流態(tài)與耗水流量之間的平衡,制定適宜的布置方案。

    [1] 陳凱麒,常仲農(nóng),曹曉紅,等. 我國魚道的建設(shè)現(xiàn)狀與展望[J]. 水利學(xué)報(bào),2012,43(2):182-188.

    Chen Kaiqi, Chang Zhongnong, Cao Xiaohong, et al. Status and prospection of fish pass construction in China[J].Journal of Hydraulic Engineering, 2012, 43(2): 182-188. (in Chinese with English abstract)

    [2] 董志勇,馮玉平,Ervine A. 同側(cè)豎縫式魚道水力特性及放魚試驗(yàn)研究[J]. 水力發(fā)電學(xué)報(bào),2008,27(6):121-125.

    Dong Zhiyong, Feng Yuping, Ervine A. An experimental study of hydraulic characteristic and fish test in vertical slot fishway to one side[J]. Journal of Hydroelectric Engineering, 2008, 27(6): 121-125. (in Chinese with English abstract)

    [3] Katopodis C, Williams J G. The development of fish passage research in a historical context[J]. Ecological Engineering, 2011, 48(7): 8-18.

    [4] 張國強(qiáng),孫雙科. 豎縫寬度對豎縫式魚道水流結(jié)構(gòu)的影響[J].水力發(fā)電學(xué)報(bào),2012,43(1):151-156.

    Zhang Guoqiang, Sun Shuangke. Effect of slot width on the flow structure of vertical slot fishway[J]. Journal of Hydroelectric Engineering, 2012, 43(1): 151-156. (in Chinese with English abstract)

    [5] Liu M, Rajaratnam N, Zhu D Z. Mean flow and turbulence structure in vertical slot fishways[J]. Journal of Hydraulic Engineering, 2006, 132(8): 765-777.

    [6] 邊永歡,孫雙科. 豎縫式魚道的水力特性研究[J]. 水利學(xué)報(bào),2013,44(12):1462-1467.

    Bian Yonghuan, Sun Shuangke. Study on hydraulic characteristic of flow in the vertical slot fishway[J]. Journal of Hydraulic Engineering, 2013, 44(12): 1462-1467. (in Chinese with English abstract)

    [7] 鄭金秀,韓德舉,胡望斌,等. 與魚道設(shè)計(jì)相關(guān)的魚類游泳行為研究[J]. 水生態(tài)學(xué)雜志,2010,3(5):104-110.

    Zheng Jinxiu, Han Deju, Hu Wangbin, et al. Fish swinmming performance related to fishway design[J]. Journal of Hydroecology, 2010, 3(5): 104-110. (in Chinese with English abstract)

    [8] DVWK. Fish Passes: Design, Dimensions and Monitoring[M]. Food and Agriculture Organization of the United Nations. Rome, 2002.

    [9] 王珂,劉紹平,段辛斌,等. 崔家營航電樞紐工程魚道過魚效果[J]. 農(nóng)業(yè)工程學(xué)報(bào),2013,29(3):184-189.

    Wang Ke, Liu Shaoping, Duan Xinbin, et al. Fishway effect of Cuijiaying navigation-power junction project[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(3): 184-189. (in Chinese with English abstract)

    [10] 吳曉春,史建全. 基于生態(tài)修復(fù)的青海湖沙柳河魚道建設(shè)與維護(hù)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2014,30(22):130-136.

    Wu Xiaochun, Shi Jianquan. Construction and management of fish passage on Shaliu River adjacent to Qinghai Lake based on ecological restoration[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(22): 130-136. (in Chinese with English abstract)

    [11] 孫雙科,張國強(qiáng). 環(huán)境友好的近自然型魚道[J]. 中國水利水電科學(xué)研究院學(xué)報(bào),2012,10(1):41-47.

    Sun Shuangke, Zhang Guoqiang. Environment-friendly fishway in close-to-nature types[J]. Journal of China Institute of Water Resources and Hydropower Research, 2012, 10(1): 41-47. (in Chinese with English abstract)

    [12] Katopodis C, Kells J A, Acharya M. Nature like and conventional fishways: Alternative concepts?[J]. Canadian Water Resources Journal, 2001, 26(2): 211-232.

    [13] 張輝,Kynard Boyd,Junho Ricardo,等. 亞馬遜流域瑪代拉河Santo Antnio 魚道設(shè)計(jì)與建造的啟示[J]. 水生態(tài)學(xué)雜志,2013,34(4):95-99.

    Zhang Hui, Kynard Boyd, Junho Ricardo, et al. Implications from the design and construction of santo antonio fishway in the madeira river in the amazon basin[J]. Journal of Hydroecology, 2013, 34(4): 95-99. (in Chinese with English abstract)

    [14] Marriner B A, Baki A B M, Zhu D Z, et al. The hydraulics of a vertical slot fishway: A case study on the multi-species Vianney-Legendre fishway in Quebec, Canada[J]. Ecological Engineering, 2016, 90: 190-202.

    [15] Jeong Hui Kim, Ju Duk Yoon, Seung Ho Baek, et al. An efficiency analysis of a nature-like fishway for freshwater fish ascending a large Korean river[J]. Water, 2015, 8(1): 1-18.

    [16] 王興勇,郭軍. 國內(nèi)外魚道研究與建設(shè)[J]. 中國水利水電科學(xué)研究院學(xué)報(bào),2005,3(3):222-228.

    Wang Xingyong, Guo Jun. Brief review on research and construction of fishways at home and abroad[J]. Journal of China Institute of Water Resources and Hydropower Research, 2005, 3(3): 222-228. (in Chinese with English abstract)

    [17] Acharya M, Kells J A, Katopodis C. Some hydraulic design aspects of nature-Like fishways[J]. Journal of Water Resources, ASCE, 2004, 5(27): 39-43.

    [18] Bretón F, Bakia B M, Link O, et al. Flow in nature-like fishway and its relation to fish behaviour[J]. Canadian Journal of Civil Engineering, 2013, 40(6): 567-573.

    [19] Baki A B M, Zhu D Z, Rajaratnam N. Mean flow characteristics in a rock-ramp-type fishpass[J]. Journal of Hydraulic Engineering, ASCE, 2014, 140(2): 156-168.

    [20] Baki A B M, Zhu D Z, Rajaratnam N. Turbulence Characteristics In A Rock-Ramp Type Fishpass[J]. Journal of Hydraulic Engineering, 2015, 141(2): 156-168.

    [21] Lacey R W J, Rennie C D. Laboratory investigation of turbulence flow structure around a bed mounted cube at multiple flowstages[J]. Journal of Hydraulic Engineering, ASCE, 2012, 138(1): 71-84.

    [22] Lacey R W J, Roy A G. The spatial characterization of turbulence around large roughness elements in a gravel-bed river[J]. Geomorphology, 2008, 102(2008): 542-553.

    [23] Tran T D, Chorda J, Laurens P, et al. Modelling nature-like fishway flow around unsubmerged obstacles using a 2D shallow water model[J]. Environmental Fluid Mechanics, 2016, 16(2): 413-428.

    [24] 何雨艨,安瑞冬,李嘉,等. 蠻石斜坡型仿自然魚道水力學(xué)特性研究[J]. 水力發(fā)電學(xué)報(bào),2016,35(10):40-47.

    He Yumeng, An Ruidong, Li Jia, et al. Hydraulic characteristics of nature-like fishways of rock-ramp type[J]. Journal of Hydroelectric Engineering, 2016, 35(10): 40-47. (in Chinese with English abstract)

    [25] 鄭鐵剛,孫雙科,柳海濤,等. 基于魚類行為學(xué)與水力學(xué)的水電站魚道進(jìn)口位置選擇[J]. 農(nóng)業(yè)工程學(xué)報(bào),2016,32(24):164-170.

    Zheng Tiegang, Sun Shuangke, Liu Haitao, et al. Location choice of fishway entrance in hydropower project based on fish behavioristics and hydraulics[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(24): 164-170. (in Chinese with English abstract)

    [26] Hirt C W, Nichols B D. Volume of Fluid (VOF) method for the dynamics of free boundaries[J]. Journal of Computational Physics, 1981, 39(1): 201-225.

    [27] FLOW Science. FLOW-3D 1v 9.4 User Manual[M]. FLOW Science, Santa Fe, NM, 2009.

    [28] Wu S, Rajaratnam N, Katopodis C. Structure of flow in vertical slot fishway[J]. Journal of Hydraulic Engineering, 1999, 125(4): 351-360.

    [29] Rajaratnam N, Katopodis C, Solanki S. New design for vertical slot fishways[J]. Canadian Journal of Civil Engineering, 2011, 19(3): 402-414.

    [30] Silva A T, Katopodis C, Santos J M, et al. Cyprinid swimming behaviour in response to turbulent flow[J]. Ecological Engineering, 2012, 44(3): 314-328.

    [31] Silva A T, Santos J M, Ferreira M T, et al. Effects of water velocity and turbulence on the behaviour of Iberian barbel in an experimental pool-type fishway[J]. River Research and Applications, 2011, 27(3): 360-373.

    Improving effect of hydraulic characteristics of nature-like fishway with pools and cobblestone weirs

    Li Guangning1, Sun Shuangke1※, Liu Haitao1, Zhang Chao1, Zhao Guixia2, Zheng Tiegang1
    (1. State Key Laboratory of Simulation of Water Cycle in River Basin, China Institute of Water Resources and Hydropower Research, Beijing 100038, China; 2. State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072, China)

    The nature-like fish way is one of the most promising layouts of fish pass structures, its remarkable characteristics are that the individual cobblestone weirs from local river bed substrates are adopted in the fishway construction to imitate as closely as possible the flow pattern of Natural River. Therefore it has higher running efficiency for fish migration than traditional technical fishways. In order to identify the differences of the hydraulic characteristics between the nature-like fishway and the traditional fishway, we conducted comparative studies. The three-dimensional mathematical models of the fishways with cobblestone weirs and water proof weirs were developed respectively. The cobblestone weirs had porosity about 12%. It was defined as the ratio of pervious area to total area of the vertical section along with the longitudinal direction. Each mathematical cobblestone feature was defined as obstacles with special shapes as the natural one. The gaps between the cobblestones allowed water to pass through existing calculation elements. A hydraulic model of the fishway with cobblestone weirs was designed based on the law of gravity similarity at a scale of 1:5. The fishway structure of the model test was consistent with the mathematical model. The water depth and velocity at typical measurement points were recorded. The results were used to verify the numerical model. In the scenario that the fish way upstream and downstream water depth was 1.4 m, the hydraulic characteristics in the fish ways, such as flow pattern, flow velocity, turbulence kinetic energy, and flow rate were compared. The results showed that the currents coming from the cobblestone gaps improved the flow dynamics conditions in the pools. The velocity distributions of the gaps in the staggered weirs appeared as the velocities of the surface were smaller than the bottom; the velocities of the sides were smaller than that of the middle. Obviously, using permeable weir fishway, the abundance of velocity distribution for fish was quite obvious, and the velocity distribution of impervious weir fishway was close to uniform. Turbulence in the mainstream can attract fish back to migration, the turbulence kinetic energy (TKE) was considered as a potential index to assess the suitability of the fishway structure. The turbulent intensity preferred by fish may depend on their swimming capabilities. In the fishways with cobblestone weirs, the TKE had greatly increased; the TKE on the surface was higher than that near the bottom. As the water surface of each of the pool was kept at the same level, the water consumption flow is very different between the impervious weir and cobblestone weir fishway. The consumption flow of the impervious weir fishway was Q=1.39 m3/s, while the consumption flow of the cobblestone weir fishway was Q=2.86 m3/s, which was increased to about twice. Thus, the nature-like fishways with cobblestone weirs needed more water flow in operating. The flow pattern of the nature-like fishway was rich and varied, and it was closer to the natural river; environment. As such it provided abundant hydrodynamic conditions for different types of migration fish. The effecting factors on the hydraulic characteristics of the nature-like fishway were analyzed, such as partition type, bottom slope, and different operation depth and so on. In general, the simulation technique of nature-like fishway was still at the exploratory stage, and the study should be further developed with a consideration with special life habits of migratory fish. The approach used in this paper can provide reference for the hydraulic design and optimization of nature-like fishway.

    fishways; flow velocity; flow rate; turbulence kinetic energy; consumption flow; numerical simulation

    10.11975/j.issn.1002-6819.2017.15.024

    TV61

    A

    1002-6819(2017)-15-0184-06

    2017-02-09

    2017-07-05

    國家自然科學(xué)基金項(xiàng)目(51679261),國家自然科學(xué)青年基金項(xiàng)目(51709278)

    李廣寧,男,河北邢臺人,工程師,博士,主要從事水力學(xué)研究工作。北京 中國水利水電科學(xué)研究院,100038。Email:94174061@qq.com

    ※通信作者:孫雙科,男,河南南陽人,教授,博士,博士生導(dǎo)師,主要從事水力學(xué)研究工作。北京 中國水利水電科學(xué)研究院,100038。

    Email:sunshuangke@126.com

    猜你喜歡
    魚道不透水隔墻
    基于主筋軸力長期監(jiān)測的連拱隧道淺埋段中隔墻受力特性研究
    四川水泥(2023年1期)2023-01-28 04:19:44
    鍋爐中隔墻中間集箱管座泄漏分析及處理措施
    基于PIT遙測技術(shù)的豎縫式魚道過魚效率及魚類行為分析
    基于無人機(jī)可見光影像與OBIA-RF算法的城市不透水面提取
    蜿蜒式魚道與豎縫式魚道的水力特性對比分析
    紊流特性對魚道結(jié)構(gòu)的影響
    卷宗(2020年16期)2020-08-10 18:21:27
    Landsat8不透水面遙感信息提取方法對比
    兩種主要技術(shù)型魚道的水力特性研究進(jìn)展
    新型高溫隔墻結(jié)構(gòu)的研究與開發(fā)
    工業(yè)爐(2016年3期)2016-10-09 03:35:50
    城市不透水面及地表溫度的遙感估算
    国产精品1区2区在线观看.| av在线亚洲专区| 毛片一级片免费看久久久久| 国产午夜精品论理片| 少妇人妻一区二区三区视频| 久久精品夜色国产| 91在线精品国自产拍蜜月| 日韩av免费高清视频| 成人国产麻豆网| av在线老鸭窝| 国模一区二区三区四区视频| 婷婷六月久久综合丁香| 色5月婷婷丁香| 丝瓜视频免费看黄片| 男女边吃奶边做爰视频| 天堂中文最新版在线下载 | 极品少妇高潮喷水抽搐| 我的女老师完整版在线观看| 国产黄频视频在线观看| 日韩,欧美,国产一区二区三区| 特大巨黑吊av在线直播| 国产精品日韩av在线免费观看| 男女啪啪激烈高潮av片| 看黄色毛片网站| 亚洲欧美日韩卡通动漫| 久久久久久国产a免费观看| 日韩精品有码人妻一区| 大片免费播放器 马上看| 国产成人精品一,二区| 国产成年人精品一区二区| 麻豆久久精品国产亚洲av| 国产白丝娇喘喷水9色精品| 国产伦在线观看视频一区| 午夜免费男女啪啪视频观看| 国产精品伦人一区二区| 日韩av免费高清视频| 国产男人的电影天堂91| 欧美xxⅹ黑人| 能在线免费看毛片的网站| or卡值多少钱| 亚洲18禁久久av| 亚洲第一区二区三区不卡| 亚洲天堂国产精品一区在线| 美女cb高潮喷水在线观看| 日韩国内少妇激情av| 丝瓜视频免费看黄片| 久久草成人影院| 国产亚洲av片在线观看秒播厂 | 亚洲国产欧美在线一区| 国产精品99久久久久久久久| 久久久a久久爽久久v久久| 国产精品嫩草影院av在线观看| 免费大片18禁| 男女下面进入的视频免费午夜| 亚洲av电影不卡..在线观看| 亚洲欧美清纯卡通| 日本av手机在线免费观看| 亚洲乱码一区二区免费版| av福利片在线观看| 久久久久性生活片| 亚洲av成人av| 草草在线视频免费看| 日日啪夜夜撸| 国产精品1区2区在线观看.| 国产欧美另类精品又又久久亚洲欧美| 国产精品1区2区在线观看.| 欧美潮喷喷水| 日本熟妇午夜| 午夜福利网站1000一区二区三区| 午夜福利成人在线免费观看| 天天躁日日操中文字幕| 99九九线精品视频在线观看视频| 色哟哟·www| 色综合站精品国产| 天天躁夜夜躁狠狠久久av| 成人鲁丝片一二三区免费| 丰满少妇做爰视频| 国产三级在线视频| 欧美精品一区二区大全| 日韩av在线免费看完整版不卡| 中文字幕免费在线视频6| 久久久久精品性色| 卡戴珊不雅视频在线播放| 少妇被粗大猛烈的视频| 欧美精品一区二区大全| 身体一侧抽搐| 精品久久久久久久久久久久久| 午夜免费男女啪啪视频观看| 高清欧美精品videossex| 黑人高潮一二区| 偷拍熟女少妇极品色| 亚洲美女视频黄频| 日韩av免费高清视频| 亚洲自偷自拍三级| 日本与韩国留学比较| av免费观看日本| 哪个播放器可以免费观看大片| 亚洲在线观看片| 少妇人妻精品综合一区二区| 亚洲最大成人手机在线| 精品久久久精品久久久| 午夜老司机福利剧场| 六月丁香七月| 国产精品熟女久久久久浪| 久久久精品94久久精品| 国产精品一区二区三区四区久久| 久久久久久久久久久免费av| 精品少妇黑人巨大在线播放| 国产日韩欧美在线精品| 午夜精品国产一区二区电影 | 久久鲁丝午夜福利片| 80岁老熟妇乱子伦牲交| 亚洲av中文字字幕乱码综合| 国产大屁股一区二区在线视频| 国产乱来视频区| 午夜激情福利司机影院| 欧美丝袜亚洲另类| 久99久视频精品免费| 在线免费观看的www视频| 岛国毛片在线播放| 国产一区二区三区综合在线观看 | 日韩欧美精品免费久久| 男人舔女人下体高潮全视频| 国产伦在线观看视频一区| 一区二区三区免费毛片| 欧美bdsm另类| 少妇丰满av| 天堂俺去俺来也www色官网 | 在线观看一区二区三区| 国产熟女欧美一区二区| 蜜桃久久精品国产亚洲av| 18禁动态无遮挡网站| 中文乱码字字幕精品一区二区三区 | 国产探花极品一区二区| 日韩电影二区| 熟妇人妻久久中文字幕3abv| 国产精品一区二区三区四区免费观看| 国产精品一区二区在线观看99 | 免费看a级黄色片| 人妻制服诱惑在线中文字幕| 精品一区二区三卡| 久久精品国产鲁丝片午夜精品| 一级av片app| 午夜福利成人在线免费观看| 国产大屁股一区二区在线视频| 一级爰片在线观看| 十八禁国产超污无遮挡网站| 午夜激情欧美在线| 日韩 亚洲 欧美在线| 女人被狂操c到高潮| 日韩精品有码人妻一区| 精品熟女少妇av免费看| 全区人妻精品视频| 精品久久久久久久久久久久久| av在线老鸭窝| 国产成人一区二区在线| 全区人妻精品视频| 国产精品无大码| 男人舔女人下体高潮全视频| 哪个播放器可以免费观看大片| 亚洲精品视频女| 波野结衣二区三区在线| 80岁老熟妇乱子伦牲交| 成人综合一区亚洲| 午夜免费激情av| 成年人午夜在线观看视频 | 国产三级在线视频| www.色视频.com| av在线蜜桃| av免费观看日本| 亚洲乱码一区二区免费版| 在线天堂最新版资源| 国产一区二区亚洲精品在线观看| 午夜福利在线在线| 亚洲精品视频女| 婷婷色综合大香蕉| av免费在线看不卡| 青青草视频在线视频观看| 一夜夜www| 色播亚洲综合网| 91精品伊人久久大香线蕉| 97热精品久久久久久| 18禁裸乳无遮挡免费网站照片| 亚洲欧洲日产国产| 成人国产麻豆网| 国产一区亚洲一区在线观看| 日本猛色少妇xxxxx猛交久久| 国产av不卡久久| 嫩草影院精品99| 国产免费视频播放在线视频 | 亚洲精品第二区| 免费黄网站久久成人精品| 91狼人影院| 精品少妇黑人巨大在线播放| 黄色日韩在线| 国产爱豆传媒在线观看| 26uuu在线亚洲综合色| 免费播放大片免费观看视频在线观看| 国产麻豆成人av免费视频| 一级毛片我不卡| 非洲黑人性xxxx精品又粗又长| 超碰97精品在线观看| 久久久a久久爽久久v久久| 国产成人精品福利久久| 亚洲经典国产精华液单| 91狼人影院| 一个人看的www免费观看视频| 亚洲精品亚洲一区二区| 永久免费av网站大全| 亚洲激情五月婷婷啪啪| 亚洲精品乱久久久久久| 亚洲aⅴ乱码一区二区在线播放| 午夜精品在线福利| 狂野欧美白嫩少妇大欣赏| 国产精品蜜桃在线观看| 国产日韩欧美在线精品| 中国美白少妇内射xxxbb| 日韩国内少妇激情av| 国产午夜精品论理片| 九草在线视频观看| 日日啪夜夜撸| av免费在线看不卡| 秋霞在线观看毛片| 欧美成人a在线观看| 2021少妇久久久久久久久久久| 麻豆精品久久久久久蜜桃| 国产在线一区二区三区精| 欧美人与善性xxx| 日本免费在线观看一区| 精华霜和精华液先用哪个| 爱豆传媒免费全集在线观看| 又粗又硬又长又爽又黄的视频| 又大又黄又爽视频免费| 三级男女做爰猛烈吃奶摸视频| 亚洲精品视频女| 精品午夜福利在线看| 国产成人a∨麻豆精品| 能在线免费观看的黄片| 国产麻豆成人av免费视频| 中文字幕制服av| 老女人水多毛片| 人妻一区二区av| 日本-黄色视频高清免费观看| 中文欧美无线码| 国产成人精品婷婷| 啦啦啦中文免费视频观看日本| 亚洲欧美中文字幕日韩二区| 国产91av在线免费观看| kizo精华| 久久这里有精品视频免费| 亚洲第一区二区三区不卡| 欧美精品一区二区大全| 热99在线观看视频| 听说在线观看完整版免费高清| 久久久久网色| 91精品一卡2卡3卡4卡| 麻豆久久精品国产亚洲av| 久久久久精品久久久久真实原创| 国产视频首页在线观看| 99久久精品国产国产毛片| 亚洲国产成人一精品久久久| 日韩视频在线欧美| 免费av不卡在线播放| 26uuu在线亚洲综合色| 91精品伊人久久大香线蕉| 中文字幕亚洲精品专区| 免费少妇av软件| 亚洲av日韩在线播放| 亚洲av免费在线观看| 成人无遮挡网站| 亚洲精品成人久久久久久| 精品人妻视频免费看| 丰满少妇做爰视频| 97超视频在线观看视频| 一级a做视频免费观看| 色综合亚洲欧美另类图片| 欧美区成人在线视频| 国产成人freesex在线| 纵有疾风起免费观看全集完整版 | 少妇人妻一区二区三区视频| 欧美性猛交╳xxx乱大交人| 建设人人有责人人尽责人人享有的 | 国产91av在线免费观看| 亚洲国产日韩欧美精品在线观看| 国产av在哪里看| 一个人免费在线观看电影| 啦啦啦韩国在线观看视频| 黄色日韩在线| 久久久久久久久久黄片| 午夜老司机福利剧场| 国产成人a区在线观看| 边亲边吃奶的免费视频| 97在线视频观看| 三级国产精品片| 日韩国内少妇激情av| 国产亚洲一区二区精品| 亚洲av中文字字幕乱码综合| 亚洲av中文av极速乱| 久久6这里有精品| 如何舔出高潮| 中国美白少妇内射xxxbb| 人人妻人人澡欧美一区二区| 亚洲国产精品专区欧美| av国产久精品久网站免费入址| 高清在线视频一区二区三区| 天堂网av新在线| 国产伦在线观看视频一区| 97人妻精品一区二区三区麻豆| 校园人妻丝袜中文字幕| 一级毛片我不卡| 51国产日韩欧美| 晚上一个人看的免费电影| 国产成人精品一,二区| 99re6热这里在线精品视频| 内射极品少妇av片p| 精品久久久久久久久亚洲| 欧美日韩视频高清一区二区三区二| 美女内射精品一级片tv| 水蜜桃什么品种好| 高清毛片免费看| 一级爰片在线观看| 色5月婷婷丁香| 肉色欧美久久久久久久蜜桃 | av在线老鸭窝| 精品酒店卫生间| 亚洲欧美日韩无卡精品| 嫩草影院精品99| 日本黄色片子视频| 色5月婷婷丁香| 天堂中文最新版在线下载 | 久久99热6这里只有精品| av在线老鸭窝| 插逼视频在线观看| 国产成人a∨麻豆精品| 麻豆乱淫一区二区| 久久久久久久久中文| ponron亚洲| 精品一区在线观看国产| 久久人人爽人人片av| 成人高潮视频无遮挡免费网站| 国产成人91sexporn| 毛片女人毛片| 最近最新中文字幕大全电影3| 成人亚洲精品一区在线观看 | 淫秽高清视频在线观看| 欧美性感艳星| 国产成人aa在线观看| 丝袜喷水一区| 欧美日韩一区二区视频在线观看视频在线 | 日韩欧美三级三区| 亚洲国产精品成人久久小说| 成人综合一区亚洲| 国内精品美女久久久久久| 久久人人爽人人爽人人片va| 在线免费观看不下载黄p国产| 女人久久www免费人成看片| 亚洲av成人av| 欧美极品一区二区三区四区| 女人久久www免费人成看片| 亚洲av成人精品一二三区| 国产人妻一区二区三区在| 亚洲丝袜综合中文字幕| 天天躁夜夜躁狠狠久久av| 久久久精品94久久精品| 久久人人爽人人片av| 中文资源天堂在线| av国产免费在线观看| 91在线精品国自产拍蜜月| 精品久久久久久成人av| 超碰av人人做人人爽久久| 三级毛片av免费| 午夜免费男女啪啪视频观看| 国产精品国产三级专区第一集| 美女xxoo啪啪120秒动态图| 日本黄色片子视频| 美女脱内裤让男人舔精品视频| 亚洲欧美成人综合另类久久久| 日日干狠狠操夜夜爽| 免费黄网站久久成人精品| 精品国产一区二区三区久久久樱花 | 街头女战士在线观看网站| 男女视频在线观看网站免费| 乱码一卡2卡4卡精品| 天天躁夜夜躁狠狠久久av| 中国国产av一级| 国产v大片淫在线免费观看| 成年版毛片免费区| 国产精品国产三级国产av玫瑰| 久久国内精品自在自线图片| 99久久精品一区二区三区| 国产精品福利在线免费观看| 69人妻影院| 观看免费一级毛片| 内射极品少妇av片p| 看非洲黑人一级黄片| 春色校园在线视频观看| 小蜜桃在线观看免费完整版高清| 日韩伦理黄色片| 国产精品嫩草影院av在线观看| 国产亚洲精品av在线| 精品一区二区免费观看| 最近的中文字幕免费完整| 秋霞伦理黄片| 久久精品久久久久久噜噜老黄| 91精品伊人久久大香线蕉| 美女被艹到高潮喷水动态| 亚洲精品成人久久久久久| 毛片女人毛片| 国产伦一二天堂av在线观看| 精品久久久久久电影网| 国产亚洲av嫩草精品影院| 国产伦理片在线播放av一区| 日韩大片免费观看网站| 国产精品福利在线免费观看| 久久人人爽人人爽人人片va| 久久人人爽人人片av| 老司机影院毛片| 男女国产视频网站| 99久久人妻综合| 自拍偷自拍亚洲精品老妇| 校园人妻丝袜中文字幕| 精品久久久噜噜| 免费看a级黄色片| 一级毛片我不卡| 高清视频免费观看一区二区 | 最近2019中文字幕mv第一页| 亚洲成色77777| a级毛片免费高清观看在线播放| 午夜精品国产一区二区电影 | 亚洲国产精品国产精品| 精品人妻偷拍中文字幕| 一个人观看的视频www高清免费观看| 黄色欧美视频在线观看| 国产精品福利在线免费观看| 午夜老司机福利剧场| 97超碰精品成人国产| 国产不卡一卡二| 美女cb高潮喷水在线观看| 看黄色毛片网站| 日韩伦理黄色片| 国产成人精品久久久久久| 国产一区二区亚洲精品在线观看| 国产一级毛片在线| 欧美日韩一区二区视频在线观看视频在线 | 欧美高清性xxxxhd video| 午夜福利高清视频| 欧美性感艳星| 亚洲在线观看片| 国产黄色小视频在线观看| 国产高清国产精品国产三级 | 日本黄大片高清| av又黄又爽大尺度在线免费看| 自拍偷自拍亚洲精品老妇| 五月玫瑰六月丁香| 又粗又硬又长又爽又黄的视频| 亚洲一级一片aⅴ在线观看| 国产真实伦视频高清在线观看| 久久精品综合一区二区三区| 国产黄片视频在线免费观看| 啦啦啦中文免费视频观看日本| 视频中文字幕在线观看| 亚洲精品影视一区二区三区av| 国精品久久久久久国模美| 欧美日韩国产mv在线观看视频 | 97超视频在线观看视频| 国产国拍精品亚洲av在线观看| 女人十人毛片免费观看3o分钟| 免费av观看视频| 九九久久精品国产亚洲av麻豆| 一级黄片播放器| 成人国产麻豆网| 色哟哟·www| 欧美成人a在线观看| 国产亚洲精品久久久com| 国产美女午夜福利| 99久久人妻综合| 18+在线观看网站| 美女脱内裤让男人舔精品视频| 亚洲av中文av极速乱| eeuss影院久久| 精品午夜福利在线看| 嫩草影院新地址| 成年人午夜在线观看视频 | 秋霞在线观看毛片| 日韩欧美国产在线观看| 亚洲精品一二三| 国产av码专区亚洲av| 中文欧美无线码| 又粗又硬又长又爽又黄的视频| 免费观看无遮挡的男女| 亚洲人与动物交配视频| 日本三级黄在线观看| 人人妻人人看人人澡| 一级毛片黄色毛片免费观看视频| 国产一区二区三区av在线| 天堂√8在线中文| 丝袜美腿在线中文| 秋霞在线观看毛片| 亚洲三级黄色毛片| 黄色日韩在线| 亚洲国产精品sss在线观看| 色播亚洲综合网| 91狼人影院| 久久国产乱子免费精品| 男女边吃奶边做爰视频| 精品久久久精品久久久| 18禁动态无遮挡网站| 亚洲精品一二三| 成人毛片a级毛片在线播放| 久久99热这里只有精品18| 久久99热这里只频精品6学生| 最近的中文字幕免费完整| 最近最新中文字幕大全电影3| 欧美不卡视频在线免费观看| 蜜臀久久99精品久久宅男| 一边亲一边摸免费视频| 免费观看的影片在线观看| 内地一区二区视频在线| 久久久久久久久久成人| 久久久亚洲精品成人影院| 亚洲18禁久久av| 插阴视频在线观看视频| 国产久久久一区二区三区| 亚洲在久久综合| 只有这里有精品99| 男女国产视频网站| 国产人妻一区二区三区在| 乱码一卡2卡4卡精品| 免费无遮挡裸体视频| 97人妻精品一区二区三区麻豆| 一级二级三级毛片免费看| 三级毛片av免费| 可以在线观看毛片的网站| 日本-黄色视频高清免费观看| 菩萨蛮人人尽说江南好唐韦庄| 一区二区三区四区激情视频| 婷婷色av中文字幕| 国产伦理片在线播放av一区| 亚洲av中文av极速乱| 噜噜噜噜噜久久久久久91| 午夜爱爱视频在线播放| 黄色欧美视频在线观看| 免费观看av网站的网址| av免费在线看不卡| 精品久久久精品久久久| 精品人妻一区二区三区麻豆| 国产精品1区2区在线观看.| 久久精品夜夜夜夜夜久久蜜豆| 搡老妇女老女人老熟妇| 狠狠精品人妻久久久久久综合| 人妻夜夜爽99麻豆av| 嫩草影院新地址| 中文字幕久久专区| 欧美日韩精品成人综合77777| 午夜精品在线福利| 欧美+日韩+精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产熟女欧美一区二区| 综合色av麻豆| 别揉我奶头 嗯啊视频| 中文字幕av成人在线电影| 丝袜美腿在线中文| 赤兔流量卡办理| 亚洲伊人久久精品综合| 天堂俺去俺来也www色官网 | 一区二区三区四区激情视频| 日日摸夜夜添夜夜添av毛片| 久久久精品94久久精品| 国产探花在线观看一区二区| 亚洲av男天堂| 丝瓜视频免费看黄片| 亚洲精品成人久久久久久| 色视频www国产| 街头女战士在线观看网站| 韩国av在线不卡| 国产一级毛片七仙女欲春2| 亚洲国产精品成人综合色| 国产精品不卡视频一区二区| av免费在线看不卡| 18禁动态无遮挡网站| 国产黄色免费在线视频| 国语对白做爰xxxⅹ性视频网站| 搞女人的毛片| 床上黄色一级片| 伊人久久国产一区二区| 亚洲国产精品sss在线观看| 精品人妻熟女av久视频| 美女脱内裤让男人舔精品视频| 啦啦啦中文免费视频观看日本| 午夜免费男女啪啪视频观看| 亚洲av免费在线观看| 国产91av在线免费观看| 国产亚洲5aaaaa淫片| 久久久久网色| 国产精品久久久久久久电影| 亚洲aⅴ乱码一区二区在线播放| 18+在线观看网站| 精品久久久久久久久av| 一级毛片aaaaaa免费看小| 亚洲av电影在线观看一区二区三区 | 亚洲精品成人久久久久久| 欧美性感艳星| 久久久久久久久久成人| 免费看美女性在线毛片视频| av在线亚洲专区| 欧美一区二区亚洲| 亚洲,欧美,日韩| 嫩草影院新地址| 亚洲在久久综合| 中国国产av一级| av专区在线播放| 国产男人的电影天堂91| 亚洲av在线观看美女高潮| 日本色播在线视频| 国产精品一二三区在线看|