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

    異齒裂腹魚(yú)通過(guò)魚(yú)道內(nèi)流速障礙能力及行為

    2018-05-11 07:42:40金志軍馬衛(wèi)忠張?jiān)瑢?/span>陳明曦譚均軍石小濤
    水利學(xué)報(bào) 2018年4期
    關(guān)鍵詞:豎縫水槽魚(yú)類

    金志軍,馬衛(wèi)忠,張?jiān)瑢?,陳明曦,譚均軍,石小濤

    (1.三峽大學(xué) 三峽庫(kù)區(qū)生態(tài)環(huán)境教育部工程研究中心,湖北 宜昌 443002;2.中國(guó)電建集團(tuán) 貴陽(yáng)勘測(cè)設(shè)計(jì)研究院有限公司,貴州 貴陽(yáng) 550081;3.中國(guó)電建集團(tuán) 成都勘測(cè)設(shè)計(jì)研究院有限公司,四川 成都 610072)

    1 研究背景

    魚(yú)道等工程措施有效運(yùn)行是保持河流縱向連通性的重要技術(shù)手段之一[1],魚(yú)類游泳能力科學(xué)定量是魚(yú)道設(shè)計(jì)關(guān)鍵問(wèn)題[2-3]。魚(yú)道設(shè)計(jì)流速根據(jù)魚(yú)類游泳能力而定,設(shè)計(jì)流速過(guò)小,建設(shè)成本增加,魚(yú)道吸引流不夠;設(shè)計(jì)流速過(guò)大,魚(yú)類無(wú)法通過(guò)魚(yú)道流速障礙。魚(yú)道長(zhǎng)度、坡度,休息室大小、擋板、過(guò)魚(yú)口等都應(yīng)基于水流流態(tài)以及魚(yú)類上溯能力進(jìn)行設(shè)計(jì)[4]。因此魚(yú)類通過(guò)流速障礙能力是魚(yú)道設(shè)計(jì)中具有指導(dǎo)意義的主要生態(tài)指標(biāo)[5-6]。

    魚(yú)類通過(guò)流速障礙能力及行為研究有助于理解魚(yú)類能否通過(guò)以及如何通過(guò)魚(yú)道內(nèi)部流速障礙。目前,評(píng)價(jià)魚(yú)類通過(guò)流速障礙能力多定量分析封閉水體(如Loligo Systems游泳能力測(cè)試水槽)中均勻流場(chǎng)下魚(yú)類被迫游泳(Forced swimming)時(shí)的游泳能力、行為和耗氧[7-8]。對(duì)于魚(yú)類游泳能力如何應(yīng)用于過(guò)魚(yú)設(shè)施設(shè)計(jì),也主要集中于游泳能力測(cè)試水槽獲得的魚(yú)類各種游泳速度指標(biāo)值。如涂志英[9]對(duì)細(xì)麟裂腹魚(yú)(體長(zhǎng)BL:9.50~13.20cm)和巨須裂腹魚(yú)(BL:19.50~28.00 cm)進(jìn)行爆發(fā)游泳速度測(cè)試,得到在25℃條件下細(xì)麟裂腹魚(yú)平均爆發(fā)游泳速度為110cm/s,6℃條件下巨須裂腹魚(yú)平均爆發(fā)游泳速度為125cm/s;葉超等[10]對(duì)異齒裂腹魚(yú)進(jìn)行臨界游泳速度和突進(jìn)游泳速度測(cè)試,得到在16℃時(shí)異齒裂腹魚(yú)臨界游泳速度為79.70~144.00cm/s(BL:26.60±11.00 cm),突進(jìn)游泳速度和體長(zhǎng)回歸擬合關(guān)系式為:y=112.106+3.156x(BL:23.60±9.60 cm)。傅菁菁等[11]對(duì)齊口裂腹魚(yú)(BL:29.00±0.10 cm)進(jìn)行突進(jìn)游泳速度測(cè)試,得到在17.2~21.6℃條件下,齊口裂腹魚(yú)突進(jìn)游泳速度為85~153 cm/s,平均值為122 cm/s。張沙龍等[12]對(duì)短須裂腹魚(yú)(BL:23.83±2.47 cm)進(jìn)行臨界游泳速度測(cè)試,得到在12.1~16.1℃條件下,短須裂腹魚(yú)臨界游泳速度為75.04±7.60 cm/s。

    眾多學(xué)者指出,管道均勻流中魚(yú)類被迫游泳測(cè)定方法不僅在生理上缺乏依據(jù),而且缺乏實(shí)踐意義,不能有效地應(yīng)用于魚(yú)類生態(tài)學(xué)定量,如其測(cè)試環(huán)境為自然界幾乎不存在的管道均勻流,魚(yú)類無(wú)法采用各種自主游泳行為,不能反映魚(yú)類在特征流場(chǎng)(如射流和渦流)下通過(guò)流速障礙的能力[13-14]。而在自然界中,魚(yú)類游泳行為是一種較不穩(wěn)定的運(yùn)動(dòng)狀態(tài),階段性的持續(xù)式游泳運(yùn)動(dòng)、靜止及偶發(fā)的爆發(fā)游泳運(yùn)動(dòng)常常相互穿插發(fā)生。以生殖洄游為例,爆發(fā)游泳為魚(yú)類越過(guò)流速障礙到達(dá)產(chǎn)卵場(chǎng)提供保障,持續(xù)式和耐久式游泳狀態(tài)則在魚(yú)類長(zhǎng)距離洄游中發(fā)揮重要作用[15]。少量學(xué)者嘗試了接近自然流態(tài)水體中的魚(yú)類通過(guò)流速障礙能力測(cè)試,但未形成受到普遍認(rèn)可的定量評(píng)價(jià)方法[16]。在應(yīng)用領(lǐng)域,人們對(duì)魚(yú)類上溯通過(guò)流速障礙能力的認(rèn)識(shí)局限,已影響到生境流場(chǎng)設(shè)計(jì),如魚(yú)道流速設(shè)計(jì)和棲息地修復(fù)等工程實(shí)踐[17-18]。盡管在封閉游泳能力測(cè)試水槽進(jìn)行魚(yú)類各種游泳速度指標(biāo)的研究遭質(zhì)疑,但其因具備操作簡(jiǎn)單和技術(shù)成熟的優(yōu)勢(shì)而被廣泛應(yīng)用。封閉游泳水槽條件下測(cè)得的游泳能力在何種程度上能夠真實(shí)反映魚(yú)類通過(guò)復(fù)雜流場(chǎng)下流速障礙能力以及如何用更加精準(zhǔn)的指標(biāo)反映自然界中魚(yú)類游泳能力一直是學(xué)者們探討的問(wèn)題。

    基于游泳能力的單因子(封閉水體中均勻流場(chǎng)下流速值)魚(yú)道設(shè)計(jì)趨勢(shì)導(dǎo)致了我們一直缺乏考慮魚(yú)類對(duì)水力條件的行為反應(yīng),這可能降低過(guò)魚(yú)設(shè)施過(guò)魚(yú)效果。在影響魚(yú)類行為的各種環(huán)境因子中,水力因子被普遍認(rèn)為是影響魚(yú)類游泳行為的主要因素[19-20]。針對(duì)水力因子對(duì)魚(yú)類游泳能力及行為影響的研究已有較多報(bào)道,主要從流速和水流紊動(dòng)等方面來(lái)探討。在不同流速下,各種魚(yú)表現(xiàn)出來(lái)的游泳行為不盡相同,如魚(yú)類在不同流速下通過(guò)改變擺尾頻率與幅度或者采用爆發(fā)-滑行游泳方式來(lái)應(yīng)對(duì)流度障礙[21-22]。而紊動(dòng)水流對(duì)魚(yú)類游泳能力和行為影響主要在于研究紊動(dòng)強(qiáng)度、紊動(dòng)能、雷諾剪切力、渦徑等,如譚均軍等[23]通過(guò)將魚(yú)在魚(yú)道中的上溯軌跡與各水力因子進(jìn)行疊加,研究了鳙魚(yú)和草魚(yú)上溯過(guò)程中的運(yùn)動(dòng)特性、水力偏好以及魚(yú)類運(yùn)動(dòng)軌跡與特定水力因子的相關(guān)程度。

    目前還未形成被廣泛認(rèn)可用于定量評(píng)價(jià)魚(yú)類通過(guò)流速障礙能力的指標(biāo)和方法。結(jié)合魚(yú)類游泳速度指標(biāo)來(lái)探索魚(yú)類自主游泳能力和有關(guān)魚(yú)類游泳行為和水力學(xué)的交叉研究提供了部分研究思路[21]。魚(yú)類通過(guò)魚(yú)道流速障礙能力測(cè)試應(yīng)基于魚(yú)類自主游泳,允許魚(yú)類使用各種游泳行為上溯或休息;其次,魚(yú)類在非均勻流場(chǎng)中如何利用流場(chǎng)將是影響魚(yú)類上溯的核心要素。為此,本文以野生異齒裂腹魚(yú)為研究對(duì)象,首先在游泳能力測(cè)試水槽中測(cè)得試驗(yàn)魚(yú)臨界游泳速度和突進(jìn)游泳速度;再以臨界游泳速度和藏木水電站魚(yú)道豎縫設(shè)計(jì)流速為參考,通過(guò)統(tǒng)計(jì)不同流態(tài)下試驗(yàn)魚(yú)通過(guò)流速障礙成功率、相對(duì)成功率、通過(guò)效率和持續(xù)爆發(fā)游泳時(shí)間,來(lái)定量試驗(yàn)魚(yú)通過(guò)流速障礙能力;同時(shí)通過(guò)將試驗(yàn)魚(yú)上溯軌跡與速度場(chǎng)進(jìn)行疊加,來(lái)探討分析試驗(yàn)魚(yú)如何利用流場(chǎng)達(dá)到上溯的目的。本文研究方法及研究結(jié)論可為魚(yú)道設(shè)計(jì)、改造、評(píng)價(jià)提供依據(jù)。

    2 試驗(yàn)設(shè)計(jì)及試驗(yàn)方法

    2.1 試驗(yàn)水槽

    2.1.1 游泳能力測(cè)試水槽 異齒裂腹魚(yú)臨界游泳速度、突進(jìn)游泳速度均在封閉游泳能力測(cè)試水槽(游泳水槽SW10150,Loligo System;丹麥)中完成[10]。

    2.1.2 魚(yú)類自主游泳能力及行為測(cè)試水槽 魚(yú)類自主游泳能力及行為測(cè)試在長(zhǎng)900 cm、寬40 cm、深30 cm的上端開(kāi)敞式可變坡水槽(圖1(a))中進(jìn)行,該水槽主體由上下回水池、中間實(shí)驗(yàn)水槽、水循環(huán)動(dòng)力系統(tǒng)組成。試驗(yàn)用水由1臺(tái)流量為100 m3/h潛水泵供給,水流經(jīng)消能整流后進(jìn)入槽體實(shí)驗(yàn)段,最后經(jīng)由尾門(mén)排入下游水池循環(huán)。選取水槽水流流態(tài)平順的中間段作為試驗(yàn)測(cè)試段,水槽末端放有百葉柵式尾門(mén)和攔魚(yú)網(wǎng)。試驗(yàn)區(qū)正上方架設(shè)3臺(tái)攝像頭(紅外網(wǎng)絡(luò)攝像頭,焦距8 mm、幀率25 Hz;??低暎?,攝像頭與試驗(yàn)水槽平行架設(shè),記錄魚(yú)類整個(gè)上溯過(guò)程。通過(guò)在水槽放置規(guī)則障礙物束窄水槽過(guò)水?dāng)嗝嬷圃旒绷鲄^(qū)和緩流區(qū),以便研究復(fù)雜流態(tài)下試驗(yàn)魚(yú)通過(guò)流速障礙能力和行為。試驗(yàn)分3種工況,其中工況1水槽坡度為1.10%,工況2和工況3水槽坡度為2.00%。魚(yú)類通過(guò)多級(jí)流速障礙試驗(yàn)(工況1、工況2)中,試驗(yàn)區(qū)固定4個(gè)障礙物(見(jiàn)圖1(b)),其中第1、第2和第3障礙物剖面為上底長(zhǎng)40 cm、下底長(zhǎng)100 cm、高27 cm的等腰梯形,第4個(gè)障礙物(靠近上游回水池)剖面為上底長(zhǎng)40 cm、下底長(zhǎng)125 cm、高27 cm的梯形,形成長(zhǎng)40 cm、寬22 cm的4級(jí)豎縫;魚(yú)類通過(guò)單級(jí)流速障礙試驗(yàn)(工況3)障礙物剖面為上底長(zhǎng)160 cm、下底長(zhǎng)245 cm、高27 cm的梯形,形成了長(zhǎng)160 cm、寬22 cm的單級(jí)豎縫(見(jiàn)圖19(c))。

    圖1 試驗(yàn)水槽正視圖以及各試驗(yàn)工況障礙物布置俯視圖(單位:cm)

    2.2 試驗(yàn)用魚(yú) 試驗(yàn)魚(yú)為異齒裂腹魚(yú),電捕于雅魯藏布江藏木水電站壩下河段(29.8oN,92.3oW),捕獲的試驗(yàn)魚(yú)分批在直徑為2.9 m的鋼化玻璃缸中暫養(yǎng),試驗(yàn)前進(jìn)行饑餓暫養(yǎng)48 h。暫養(yǎng)水取自雅魯藏布江,水溫為14.6±1.1℃,全天不間斷充氧,溶解氧大于6.0 mg/L。從大量漁獲物中挑選出未受傷、體質(zhì)健康的樣本用于試驗(yàn),共147尾。其中18尾(BL=19.41±2.33 cm、濕重Wg=96.87±39.81 g)用于突進(jìn)游泳速度測(cè)試,21尾(BL=20.96±2.66 cm、Wg=137.96±45.16 g)用于臨界游泳速度測(cè)試,39尾(BL=20.02±1.86 cm,Wg=116.45±32.48 g)用于工況1通過(guò)4級(jí)流速障礙能力和行為試驗(yàn),39尾(BL=21.38±2.71cm,Wg=139.05±51.17 g)用于工況2通過(guò)4級(jí)流速障礙能力和行為試驗(yàn),30尾(BL=22.61±2.09 cm,Wg=154.73±43.85 g)用于工況3通過(guò)單級(jí)流速障礙持續(xù)爆發(fā)游泳能力試驗(yàn)。

    2.3 試驗(yàn)方法

    2.3.1 魚(yú)類臨界、突進(jìn)游泳能力測(cè)試 通過(guò)流速遞增法測(cè)試異齒裂腹魚(yú)臨界游泳速度和突進(jìn)游泳速度[10,24]。臨界游泳速度測(cè)試:試驗(yàn)前估計(jì)試驗(yàn)魚(yú)體長(zhǎng),試驗(yàn)魚(yú)在流速為1 BL/s條件下適應(yīng)60 min,適應(yīng)結(jié)束后流速增加1 BL/s,此后,流速每隔20 min增加1 BL/s,直到試驗(yàn)魚(yú)疲勞貼網(wǎng)超過(guò)10 s后結(jié)束試驗(yàn)。突進(jìn)游泳速度測(cè)試:試驗(yàn)方法類似于臨界游泳速度測(cè)試,區(qū)別在于時(shí)間增量t為20 s[8,10]。臨界游泳速度和突進(jìn)游泳速度由下面公式得到:

    式中:Ui為第i尾實(shí)驗(yàn)魚(yú)疲勞時(shí)前一流速,cm/s;▽v為速度增幅,1 BL/s;t為時(shí)間增量,臨界游泳能力測(cè)試為20 min、突進(jìn)游泳能力測(cè)試為20 s;ti為實(shí)驗(yàn)魚(yú)在最高流速下游泳時(shí)間,s;分別表示第i尾實(shí)驗(yàn)魚(yú)臨界游泳速度和突進(jìn)游泳速度。

    2.3.2 魚(yú)類自主游泳能力及行為測(cè)試 每次試驗(yàn)將一尾試驗(yàn)魚(yú)放入水槽下游攔網(wǎng)后適應(yīng)10 min,適應(yīng)結(jié)束后開(kāi)始正式試驗(yàn)。試驗(yàn)魚(yú)通過(guò)第4級(jí)豎縫和通過(guò)長(zhǎng)160 cm的單級(jí)豎縫則試驗(yàn)結(jié)束,且每次試驗(yàn)時(shí)間不超過(guò)60 min。試驗(yàn)水槽底部貼有與魚(yú)體色有較大差異的白色反光膜,以便對(duì)魚(yú)類運(yùn)動(dòng)軌跡進(jìn)行視頻追蹤定位。試驗(yàn)后,截取試驗(yàn)魚(yú)通過(guò)豎縫的游泳視頻,并通過(guò)Logger Pro 3.12軟件提取試驗(yàn)魚(yú)上溯軌跡坐標(biāo)和通過(guò)豎縫所需時(shí)間(精確到0.04 s)。

    通過(guò)多級(jí)流速障礙能力及行為試驗(yàn)過(guò)程中記錄試驗(yàn)魚(yú)通過(guò)每級(jí)豎縫的時(shí)間、成功通過(guò)豎縫次數(shù)、嘗試通過(guò)(魚(yú)頭部進(jìn)入豎縫,但未通過(guò)的情況)次數(shù)以及計(jì)算通過(guò)多級(jí)流速障礙成功率(通過(guò)每級(jí)豎縫試驗(yàn)魚(yú)尾數(shù)占總試驗(yàn)魚(yú)尾數(shù)百分比)、相對(duì)成功率(成功通過(guò)每級(jí)豎縫試驗(yàn)魚(yú)尾數(shù)占成功通過(guò)上一級(jí)豎縫試驗(yàn)魚(yú)尾數(shù)百分比)和通過(guò)效率(每尾試驗(yàn)魚(yú)成功通過(guò)豎縫總次數(shù)占嘗試通過(guò)總次數(shù)和成功通過(guò)總次數(shù)之和的百分比)。為對(duì)比研究不同豎縫流速、不同豎縫長(zhǎng)度下魚(yú)類游泳速度,將長(zhǎng)度為160 cm單級(jí)豎縫劃分為4個(gè)等級(jí),即豎縫長(zhǎng)度為40、80、120、160 cm這4個(gè)等級(jí),同樣記錄試驗(yàn)魚(yú)成功通過(guò)豎縫和嘗試通過(guò)的時(shí)間。

    2.4 水力特性分析 試驗(yàn)水力條件通過(guò)聲學(xué)多普勒點(diǎn)式流速儀(小威龍Vectrino;Nortek)進(jìn)行測(cè)量,試驗(yàn)區(qū)水槽水深在14~22 cm之間,測(cè)點(diǎn)斷面距離槽底6 cm,各測(cè)點(diǎn)測(cè)量頻率為30 Hz,測(cè)量時(shí)間在90~120 s之間。通過(guò)多級(jí)流速障礙能力及行為試驗(yàn)中各樣本點(diǎn)相距3~5 cm,共112個(gè)橫向測(cè)量斷面、9個(gè)縱向測(cè)量斷面,765個(gè)樣本點(diǎn)。通過(guò)單級(jí)流速障礙試驗(yàn)中各樣本點(diǎn)相距3~10 cm,共42個(gè)橫向測(cè)量斷面、9個(gè)縱向測(cè)量斷面,242個(gè)樣本點(diǎn)。通過(guò)Win ADV軟件對(duì)測(cè)點(diǎn)的u、v、w三方向流速數(shù)據(jù)進(jìn)行處理,剔除信噪比小于15、相關(guān)度小于70的數(shù)據(jù)后再進(jìn)行流場(chǎng)分析。

    瞬時(shí)速度與時(shí)間平均流速的差值為脈動(dòng)流速(Win ADV),公式表示為:

    式中:uk為瞬時(shí)速度,cm/s;為時(shí)間平均速度,cm/s;u′k為脈動(dòng)速度,cm/s。

    測(cè)點(diǎn)速度大小表示為:

    豎縫平均流速U為豎縫處各測(cè)點(diǎn)流速均值:

    試驗(yàn)魚(yú)通過(guò)豎縫游泳速度:

    式中:D為豎縫長(zhǎng)度,cm;T為通過(guò)豎縫所需持續(xù)游泳時(shí)間,s。

    曹娜等[25]對(duì)國(guó)內(nèi)2000—2014年經(jīng)過(guò)環(huán)境保護(hù)部環(huán)境工程評(píng)估中心技術(shù)評(píng)估的35個(gè)水利水電建設(shè)項(xiàng)目設(shè)計(jì)或已建魚(yú)道工程進(jìn)行統(tǒng)計(jì)分析,表明35座魚(yú)道除老龍口水利樞紐工程魚(yú)道和西牛航運(yùn)樞紐魚(yú)道的設(shè)計(jì)流速大于1.2 m/s外,其余均不大于1.2 m/s;同時(shí)根據(jù)藏木水電站魚(yú)道豎縫設(shè)計(jì)流速(110 cm/s)[26],試驗(yàn)設(shè)計(jì)了小、中、大3種不同流態(tài)試驗(yàn)工況(見(jiàn)圖2)。工況1豎縫流速為101.55±14.87 cm/s(58.58~ 128.25 cm/s),工況2豎縫速度為114.63± 24.28 cm/s(52.19~ 156.94 cm/s),工況3豎縫流速為137.45±17.63 cm/s(84.72~159.90 cm/s)。試驗(yàn)工況、試驗(yàn)水溫以及各豎縫平均流速如表1。

    3 試驗(yàn)結(jié)果與分析

    3.1 通過(guò)流速障礙能力

    3.1.1 臨界、突進(jìn)游泳速度 魚(yú)類有氧運(yùn)動(dòng)能力評(píng)價(jià)指標(biāo)主要由Ucrit來(lái)表示,而在較短時(shí)間內(nèi)的最大游泳速度即Ubrust是評(píng)價(jià)無(wú)氧運(yùn)動(dòng)主要指標(biāo),這兩個(gè)指標(biāo)是魚(yú)道設(shè)計(jì)的重要參數(shù)。試驗(yàn)結(jié)果表明:在水溫為14.79±1.72℃下試驗(yàn)魚(yú)臨界游泳速度為101.01±20.86 cm/s(69.00~148.00 cm/s)(見(jiàn)圖3);在水溫為15.95±0.62℃下突進(jìn)游泳速度為196.94±21.80 cm/s(145.00~237.00 cm/s)(見(jiàn)圖3)。臨界、突進(jìn)游泳速度測(cè)試結(jié)果與葉超等[10]研究結(jié)果相當(dāng)。

    圖2 各試驗(yàn)工況流速等值線云圖

    表1 試驗(yàn)工況、豎縫平均流速情況

    圖3 異齒裂腹魚(yú)臨界、突進(jìn)游泳速度

    3.1.2 自主上溯游泳能力 對(duì)于成功通過(guò)四級(jí)豎縫的試驗(yàn)魚(yú),統(tǒng)計(jì)其通過(guò)每級(jí)豎縫所需持續(xù)游泳時(shí)間。通過(guò)工況1、工況2各級(jí)豎縫持續(xù)游泳時(shí)間如表2。兩種工況下通過(guò)豎縫持續(xù)游泳時(shí)間具有顯著差異(Mann-WhitneyU=285,P<0.01),其中工況1每尾試驗(yàn)魚(yú)通過(guò)各級(jí)豎縫持續(xù)游泳時(shí)間為0.62±0.28 s,工況2為1.08±0.68 s。試驗(yàn)魚(yú)通過(guò)流速大于其臨界游泳速度(101.01 cm/s)的豎縫所需持續(xù)游泳時(shí)間為0.52±0.34 s。

    表2 工況1、工況2試驗(yàn)魚(yú)通過(guò)每級(jí)豎縫持續(xù)游泳時(shí)間

    通過(guò)單級(jí)流速障礙能力試驗(yàn)中有28尾試驗(yàn)魚(yú)通過(guò)豎縫,1尾嘗試通過(guò)。通過(guò)不同長(zhǎng)度豎縫,對(duì)應(yīng)持續(xù)游泳時(shí)間、平均游泳速度如表3。試驗(yàn)魚(yú)通過(guò)不同長(zhǎng)度豎縫游泳速度為215.18±18.39 cm/s,且無(wú)顯著性差異(one-way ANOVA,P>0.05)。通過(guò)不同豎縫長(zhǎng)度與對(duì)應(yīng)可通過(guò)的流速的關(guān)系可擬合為:y=158.30-0.11x(R2=0.59,P<0.01),通過(guò)不同豎縫長(zhǎng)度與所需持續(xù)爆發(fā)游泳時(shí)間的關(guān)系可擬合為:y=0.06+0.01x(R2=0.61,P<0.01),見(jiàn)圖4。

    表3 通過(guò)不同豎縫長(zhǎng)度對(duì)應(yīng)可通過(guò)流速以及持續(xù)爆發(fā)游泳能力

    圖4 通過(guò)不同障礙長(zhǎng)度所需持續(xù)游泳時(shí)間和可通過(guò)流速

    試驗(yàn)魚(yú)通過(guò)流速為106.05~152.81cm/s的豎縫時(shí),游泳速度無(wú)顯著性差異(one-way ANOVA,P>0.05)(見(jiàn)圖5),值為214.01±30.64cm/s,且與突進(jìn)游泳速度(196.94cm/s)無(wú)顯著性差異(one-way ANO?VA,P>0.05)。圖5中不同字母表示通過(guò)不同流速豎縫時(shí)游泳速度具有顯著性差異(P<0.05);虛線分別表示豎縫流速為臨界游泳速度和游泳速度為突進(jìn)游泳速度??梢?jiàn)在本試驗(yàn)條件下,試驗(yàn)魚(yú)通過(guò)流速大于其臨界游泳速度的豎縫時(shí),以與突進(jìn)游泳速度無(wú)顯著性差異的恒定游泳速度上溯。

    圖5 3種工況下試驗(yàn)魚(yú)通過(guò)不同流速豎縫時(shí)游泳速度

    目前,我國(guó)絕大多數(shù)垂直豎縫式魚(yú)道過(guò)魚(yú)口流速低于1.2m/s[1,25]。本試驗(yàn)結(jié)果表明,93.33%試驗(yàn)魚(yú)以209.43±21.76 cm/s游泳速度成功通過(guò)長(zhǎng)度為160 cm、流速為137.45±17.63 cm/s的豎縫。綜合3種工況,110 cm/s的魚(yú)道豎縫流速對(duì)異齒裂腹魚(yú)上溯不構(gòu)成流速障礙。

    3.2 通過(guò)多級(jí)流速障礙上溯行為

    3.2.1 通過(guò)多級(jí)流速障礙成功率和通過(guò)效率 成功率和相對(duì)成功率是魚(yú)道過(guò)魚(yú)效果整體評(píng)價(jià)和問(wèn)題池室監(jiān)測(cè)的主要指標(biāo)。工況1試驗(yàn)魚(yú)通過(guò)豎縫成功率從第1級(jí)的87.18%降到第4級(jí)82.05%,工況2從第1級(jí)92.31%降到第4級(jí)84.62%(見(jiàn)圖6),兩種工況下成功率無(wú)顯著差異性(one-way ANOVA,P>0.05)。工況1和工況2相對(duì)成功率分別為95.30±5.60%、95.59±3.32%,兩種工況通過(guò)豎縫相對(duì)成功率無(wú)顯著差異性(one-way ANOVA,P>0.05)。兩種工況通過(guò)第1級(jí)豎縫相對(duì)成功率低于通過(guò)第2級(jí)、第3級(jí)、第4級(jí)豎縫相對(duì)成功率;工況1第3級(jí)豎縫和第4級(jí)豎縫流速大于第1級(jí)和第2級(jí)豎縫流速,通過(guò)第3級(jí)豎縫的試驗(yàn)魚(yú)全通過(guò)第4級(jí)豎縫;工況2第2級(jí)豎縫和第3級(jí)豎縫流速大于第1級(jí)和第4級(jí)豎縫流速,通過(guò)第2級(jí)豎縫的試驗(yàn)魚(yú)全通過(guò)第3級(jí)豎縫(見(jiàn)圖6)。

    圖6 工況1、工況2下試驗(yàn)魚(yú)通過(guò)每級(jí)豎縫成功率、相對(duì)成功率

    工況1有31尾試驗(yàn)魚(yú)通過(guò)效率為100%,工況2下有21尾試驗(yàn)魚(yú)通過(guò)效率為100%。工況1、工況2通過(guò)效率分別為97.62±8.23%、84.99±21.38%(見(jiàn)圖7)。兩種工況通過(guò)效率(工況2通過(guò)效率小于工況1)具有顯著差異(Mann-WhitneyU=407.5,P<0.05)的主要原因可能是工況2試驗(yàn)魚(yú)嘗試通過(guò)豎縫次數(shù)高于工況1的嘗試次數(shù)。

    對(duì)試驗(yàn)魚(yú)連續(xù)通過(guò)4級(jí)豎縫所需游泳時(shí)間進(jìn)行統(tǒng)計(jì)。工況1連續(xù)通過(guò)4級(jí)豎縫所需時(shí)間為9.08±4.77 s,工況2為11.73±7.31 s(見(jiàn)圖8)。兩種工況下試驗(yàn)魚(yú)從進(jìn)入第1級(jí)豎縫到通過(guò)第4級(jí)豎縫所需時(shí)間無(wú)顯著性差異(F1,63=2.98;P>0.05)。第1級(jí)豎縫進(jìn)口到第4級(jí)豎縫出口直線距離為460 cm,工況1通過(guò)所需最短時(shí)間為1.96 s,工況2為4.52s。部分試驗(yàn)魚(yú)以極高速度通過(guò)多級(jí)豎縫,可能與上溯過(guò)程中利用流場(chǎng)加快對(duì)地游泳速度,減少上溯時(shí)間有關(guān)。

    圖7 工況1、工況2試驗(yàn)魚(yú)通過(guò)多級(jí)豎縫通過(guò)效率

    圖8 工況1、工況2試驗(yàn)魚(yú)通過(guò)豎縫所需游泳時(shí)間

    若游泳能力是魚(yú)類是否通過(guò)魚(yú)道流速障礙決定性因素,則在低流速工況下成功率應(yīng)更高,而實(shí)際并非如此[27-28]。魚(yú)類必須還愿意且積極地嘗試通過(guò),這受包括生理?xiàng)l件、嗅覺(jué)信號(hào)和魚(yú)對(duì)水流反應(yīng)在內(nèi)的多因素影響[29]。工況1和工況2豎縫流速均接近藏木魚(yú)道豎縫設(shè)計(jì)流速,兩種工況通過(guò)豎縫成功率、相對(duì)成功率以及通過(guò)4級(jí)豎縫所需時(shí)間無(wú)顯著性差異。后期研究中,應(yīng)設(shè)置更廣的豎縫流速范圍,以確定魚(yú)類上溯效果最佳的流速工況。

    3.2.2 不同運(yùn)動(dòng)路徑對(duì)上溯效率的影響 從工況2下33尾成功通過(guò)4級(jí)豎縫的上溯軌跡中提取31條(剔除非連續(xù)上溯軌跡)從第2級(jí)豎縫出口進(jìn)入第3級(jí)豎縫進(jìn)口的上溯軌跡進(jìn)行分析。將上溯軌跡分為3類:通過(guò)低流速區(qū)(回流區(qū))進(jìn)入下一級(jí)豎縫,記為第1類軌跡,共13條(見(jiàn)圖9(a));從主流(高流速區(qū))進(jìn)入下一級(jí)豎縫,記為第2類軌跡,共10條(見(jiàn)圖9(b));其他上溯軌跡記為第3類,共8條(見(jiàn)圖9(c))。統(tǒng)計(jì)各類軌跡通過(guò)池室所需時(shí)間以及各軌跡上溯路徑長(zhǎng)度。結(jié)果表明:第1類軌跡上溯所需時(shí)間小于第2、第3類且具有顯著差異性,第1類上溯路徑長(zhǎng)度小于第2類、第3類軌跡且具有顯著性差異,如表4。

    表4 不同類型上溯軌跡上溯時(shí)間以及上溯路徑長(zhǎng)度

    圖9 工況2下上溯軌跡與第2級(jí)池室水流速度耦合

    相比第2類、第3類上溯軌跡,第1類軌跡處于低流速區(qū)且試驗(yàn)魚(yú)運(yùn)動(dòng)方向與水流方向相同(見(jiàn)圖10)。試驗(yàn)魚(yú)借助自身初始速度以及回流區(qū)水流推動(dòng),無(wú)需擺動(dòng)尾鰭亦可通過(guò)池室,進(jìn)入下一級(jí)豎縫。選擇第1類軌跡上溯的試驗(yàn)魚(yú)不僅上溯效率高于其他上溯路徑,游泳耗能也會(huì)降低。

    4 結(jié)論

    圖10 工況2下第1類上溯軌跡與第2級(jí)池室速度矢量耦合

    本文以臨界游泳速度(101.01 cm/s)和藏木水電站魚(yú)道豎縫流速(設(shè)計(jì)值為110.00 cm/s)為參考,通過(guò)在開(kāi)放游泳水槽內(nèi)加不同束窄梯形體,開(kāi)展兩種底坡條件下4級(jí)短豎縫(工況1和工況2豎縫流度為101.55±14.87 cm/s、114.63±24.28 cm/s,豎縫順?biāo)鏖L(zhǎng)度均為40 cm)和單級(jí)長(zhǎng)豎縫(工況3豎縫流速為137.45±17.63 cm/s、豎縫順?biāo)鏖L(zhǎng)度為160 cm)下異齒裂腹魚(yú)通過(guò)流速障礙能力和行為研究。通過(guò)統(tǒng)計(jì)不同流態(tài)下通過(guò)流速障礙成功率、相對(duì)成功率、通過(guò)效率和持續(xù)爆發(fā)游泳時(shí)間,定量了試驗(yàn)魚(yú)通過(guò)流速障礙能力。同時(shí)通過(guò)將試驗(yàn)魚(yú)上溯軌跡與速度場(chǎng)進(jìn)行耦合,分析了試驗(yàn)魚(yú)利用流場(chǎng)達(dá)到上溯目的的行為。具體結(jié)論如下:

    (1)工況1異齒裂腹魚(yú)通過(guò)豎縫成功率從第1級(jí)的87.18%降到第4級(jí)的82.05%,工況2成功率從第1級(jí)92.31%降到第4級(jí)的84.62%。兩種工況成功率無(wú)顯著差異性(one-way ANOVA,P>0.05)。工況1和工況2相對(duì)成功率分別為95.30±5.60%、95.59±3.32%,無(wú)顯著差異性(one-way ANOVA,P>0.05)。通過(guò)第1級(jí)豎縫相對(duì)成功率低于第2級(jí)、第3級(jí)、第4級(jí)。

    (2)工況1和工況2通過(guò)效率分別為97.62±8.23%、84.99±21.38%,具有顯著性差異(Mann-Whit?neyU=407.5,P<0.01)。工況1試驗(yàn)魚(yú)通過(guò)每級(jí)豎縫時(shí)間為0.62±0.28 s,工況2為1.08±0.68 s,具有顯著差異性(Mann-WhitneyU=285,P<0.01)。工況1和工況2試驗(yàn)魚(yú)從第1級(jí)豎縫進(jìn)口上溯到第4級(jí)豎縫出口所需游泳時(shí)間分別為9.08±4.77 s、11.73±7.31s,且無(wú)顯著性差異(F1,63=2.98;P>0.05);工況1從第1級(jí)豎縫進(jìn)口到第4級(jí)豎縫出口(直線距離為420 cm)所需最短時(shí)間為1.96s,工況2為4.52 s。兩種工況下,試驗(yàn)魚(yú)通過(guò)流速大于其臨界游泳速度(101.01 cm/s)的豎縫所需持續(xù)游泳時(shí)間為0.52±0.34s。

    (3)試驗(yàn)魚(yú)通過(guò)平均流速為106.05~152.81 cm/s的豎縫時(shí),游泳速度無(wú)顯著性差異(one-way ANOVA,P>0.05),值為214.01± 30.64 cm/s,且與突進(jìn)游泳速度(196.94 cm/s)無(wú)顯著性差異(oneway ANOVA,P>0.05)。93.33%試驗(yàn)魚(yú)以209.43±21.76 cm/s游泳速度成功通過(guò)長(zhǎng)度為160 cm、流速為137.45±17.63 cm/s的單級(jí)流速障礙。

    (4)設(shè)計(jì)流速為110.00 cm/s的藏木水電站魚(yú)道豎縫流速對(duì)異齒裂腹魚(yú)上溯不構(gòu)成流速障礙,但綜合上溯效率和魚(yú)道建設(shè)成本,還須設(shè)置覆蓋面更廣的豎縫流速范圍,進(jìn)一步確定魚(yú)類上溯效率最佳的水力條件。

    (5)魚(yú)類游泳軌跡與流場(chǎng)耦合分析表明:不同運(yùn)動(dòng)路徑選擇對(duì)上溯游泳效率具有顯著的影響(one-way ANOVA,P<0.01)。試驗(yàn)魚(yú)通過(guò)借助同向水流的推動(dòng),減少上溯所需時(shí)間和上溯路徑長(zhǎng)度,從而增加上溯效率。在復(fù)雜流態(tài)中,除了流速,紊動(dòng)能、紊動(dòng)強(qiáng)度、雷諾剪切力和渦徑等水力因子也可能影響魚(yú)類游泳行為。在未來(lái)研究中應(yīng)著重于分析復(fù)雜流場(chǎng)下魚(yú)類游泳行為和游泳耗能,以及通過(guò)定量魚(yú)類游泳姿態(tài)失穩(wěn)處或能量消耗突增處的紊流閾值來(lái)修訂魚(yú)道水力設(shè)計(jì)標(biāo)準(zhǔn)值。

    參 考 文 獻(xiàn):

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

    [2]BAINBRIDGE R.Speed and stamina in three fish[J].Journal of Experimental Biology,1960,37:129-153.

    [3]BRETT J R.The respiratory metabolism and swimming performance of young sockeye salmon[J].Journal of the Fisheries Research Board of Canada,1964,21:1183-226.

    [4]HARO A,CASTROSANTOS T,NOREIKA J,et al.Swimming performance of upstream migrant fishes in openchannel flow:a new approach to predicting passage through velocity barriers[J].Canadian Journal of Fisheries&Aquatic Sciences,2004,61(9):1590-1601.

    [5]NATIONS F AOO.Fish Passes:Design,Dimensions and Monitoring[M].Food and Agriculture Organization in Arrangement with DVWK,2002.

    [6]石小濤,陳求穩(wěn),黃應(yīng)平,等.魚(yú)類通過(guò)魚(yú)道內(nèi)水流速度障礙能力的評(píng)估方法分析[J].生態(tài)學(xué)報(bào),2011,31(22):6967-6972.

    [7]KIEFFER J D,ARSENAULT LM,LITVAK M K.Behaviour and performance of juvenile shortnose sturgeon Acipenser brevirostrum at different water velocities[J].Journal of Fish Biology,2009,74:674-682.

    [8]金志軍,陳小龍,王從鋒,等.應(yīng)用于魚(yú)道設(shè)計(jì)的馬口魚(yú)(Opsariichthy sbidens)游泳能力[J].生態(tài)學(xué)雜志,2017,36(9):2678-2684.

    [9]涂志英.雅礱江流域典型魚(yú)類游泳特性研究[D].武漢:武漢大學(xué),2012.

    [10]葉超,王珂,黃福江,等.異齒裂腹魚(yú)游泳能力初探[J].淡水漁業(yè),2013,43(3):33-37.

    [11]傅菁菁,李嘉,安瑞冬,等.基于齊口裂腹魚(yú)游泳能力的豎縫式魚(yú)道流態(tài)塑造研究[J].四川大學(xué)學(xué)報(bào)(工程科學(xué)版),2013,45(3):12-17.

    [12]張沙龍,侯軼群,王龍濤,等.長(zhǎng)絲裂腹魚(yú)的游泳能力和游泳行為研究[J].淡水漁業(yè),2014(5):32-37.

    [13]CASTROSANTOS T,HARO A.Fish Guidance and Passage at Barriers[M]//Fish Locomotion:An Eco-ethologi?cal Perspective.New Hampshire:Science Publishers,2010.

    [14]SANZRONDA F J,CóRDOBA F J B,F(xiàn)UENTES J F,et al.Ascent ability of brown trout,Salmo trutta,and two Iberian cyprinids-Iberian barbel,Luciobarbus bocagei,and northern straight-mouth nase,Pseudochondrostoma duriense-in a vertical slot f i shway[J].Knowledge&Management of Aquatic Ecosystems,2016(417):1-9.

    [15]CASTROSANTOS T.Optimal swim speeds for traversing velocity barriers:an analysis of volitional high-speed swimming behavior of migratory fishes[J].Journal of Experimental Biology,2005,208(3):421-432.

    [16]CASTROSANTOS T,SANZRONDA F J,RUIZLEGAZPI J.Breaking the speed limit—comparative sprinting per?formance of brook trout(Salvelinus fontinalis)and brown trout(Salmo trutta)[J].Canadian Journal of Fisheries and Aquatic Sciences,2012,70(2):280-293.

    [17]KEMP P S.Bridging the gap between fish behaviour,performance and hydrodynamics:an ecohydraulics approach to fish passage research[J].River Research and Applications,2012,28(4):403-406.

    [18]KATOPODIS C,GERVAIS R.Fish swimming performance database and analyses[R].Canadian Science Adviso?ry Secretariat(CSAS)Research Document,2016.

    [19]SILVA A T,SANTOS J M,F(xiàn)ERREIRA M T,et al.Effects of water velocity and turbulence on the behaviour of Iberian barbel(Luciobarbus bocagei,Steindachner 1864)in an experimental pool-type fishway[J].River Re?search&Applications,2011,27(3):360-373.

    [20]SILVA A T,SANTOS J M,F(xiàn)ERREIRA M T,et al.Passage efficiency of offset and straight orifices for upstream movements of Iberian barbel in a pool-type fishway[J].River Research and Applications,2012,28(5):529-542.

    [21]KEMP P S,WILLIAMS J G.Illumination influences the ability of migrating juvenile salmonids to pass a sub?merged experimental weir[J].Ecology of Freshwater Fish,2009,18(2):297-304.

    [22]WINGER P D,HE P,WALSH S.Factors affecting the swimming endurance and catchability of Atlantic cod(Ga?dus morhua)[J].Canadian Journal of Fisheries&Aquatic Sciences,2011,57:1200-1207.

    [23]譚均軍,高柱,戴會(huì)超,等.豎縫式魚(yú)道水力特性與魚(yú)類運(yùn)動(dòng)特性相關(guān)性分析[J].水利學(xué)報(bào),2017,48(8):924-932.

    [24]BRETT J R.Swimming performance of sockeye salmon(Oncorhynchus nerka)in relation to fatigue time and tem?perature[J].Journal of the Fisheries Board of Canada,1967,24(8):1731-1741.

    [25]曹娜,鐘治國(guó),曹曉紅,等.我國(guó)魚(yú)道建設(shè)現(xiàn)狀及典型案例分析[J].水資源保護(hù),2016,32(6):156-162.

    [26]陳靜,郎建,周小波,等.雅魯藏布江藏木水電站魚(yú)道工程設(shè)計(jì)與研究[J].水電站設(shè)計(jì),2017,33(1):52-58.

    [27]CASTROSANTOS T.Quantifying the combined effects of attempt rate and swimming capacity on passage through velocity barriers[J].Canadian Journal of Fisheries and Aquatic Sciences,2004,61:1602-1615.

    [28]CHEONG T S,KAVVAS M L,ANDERSON E K.Evaluation of adult white sturgeon swimming capabilities and applications to fishway design[J].Environmental Biology of Fishes,2006,77:197-208.

    [29]CASTROSANTOS T,SANZRONDA J,RUIZLEGAZPI J.Breaking the speed limit-comparative sprinting perfor?mance of brook trout(Salvelinus fontinalis)and brown trout(Salmo trutta)[J].Canadian Journal of Fisheries and Aquatic Sciences,2013,70:280-293.

    猜你喜歡
    豎縫水槽魚(yú)類
    豎縫位置對(duì)異側(cè)豎縫式魚(yú)道水力特性的數(shù)值模擬與分析
    基于MFCC和ResNet的魚(yú)類行為識(shí)別
    可升降折疊的飲水機(jī)水槽
    裝配式建筑剪力墻豎縫連接形式及其研究現(xiàn)狀
    河南建材(2021年3期)2021-01-07 07:00:41
    異側(cè)豎縫寬度對(duì)魚(yú)道水流結(jié)構(gòu)的影響研究
    帶豎縫鋼管束砼組合剪力墻受力性能模擬分析
    為什么水槽管要做成彎曲狀
    奇妙的古代動(dòng)物 泥盆紀(jì)的魚(yú)類
    要挑好水槽,就看這里了!
    幸福(2016年6期)2016-12-01 03:08:13
    魚(yú)類是怎樣保護(hù)自己的
    亚洲精品国产一区二区精华液| 亚洲在线自拍视频| 欧美中文日本在线观看视频| 啪啪无遮挡十八禁网站| 久久久精品欧美日韩精品| 免费电影在线观看免费观看| 亚洲精品国产精品久久久不卡| 舔av片在线| 久久精品国产亚洲av高清一级| 非洲黑人性xxxx精品又粗又长| 天堂影院成人在线观看| 欧美不卡视频在线免费观看 | 欧美3d第一页| 人人妻人人看人人澡| 人妻夜夜爽99麻豆av| 91麻豆精品激情在线观看国产| 大型av网站在线播放| 最新在线观看一区二区三区| 成人永久免费在线观看视频| 成人三级黄色视频| 久久精品国产亚洲av高清一级| 成人国语在线视频| 人人妻人人看人人澡| 日韩欧美 国产精品| 亚洲av美国av| 1024视频免费在线观看| 日韩欧美国产在线观看| 一区二区三区激情视频| 亚洲,欧美精品.| 夜夜爽天天搞| 国产人伦9x9x在线观看| 日本黄色视频三级网站网址| 长腿黑丝高跟| 国产麻豆成人av免费视频| 午夜精品一区二区三区免费看| 九色国产91popny在线| 国产片内射在线| 国产91精品成人一区二区三区| 黄色视频不卡| 中文在线观看免费www的网站 | 久久久精品国产亚洲av高清涩受| 欧美乱妇无乱码| 又爽又黄无遮挡网站| 国产精品一区二区免费欧美| 精品一区二区三区av网在线观看| 少妇熟女aⅴ在线视频| 老熟妇仑乱视频hdxx| 1024香蕉在线观看| 香蕉丝袜av| 91国产中文字幕| 免费看日本二区| 亚洲人成网站高清观看| 国内精品久久久久精免费| 亚洲精品一卡2卡三卡4卡5卡| 欧美乱色亚洲激情| 成人国产一区最新在线观看| 日本黄色视频三级网站网址| 18禁黄网站禁片免费观看直播| 长腿黑丝高跟| 久久午夜亚洲精品久久| 亚洲国产中文字幕在线视频| 精品高清国产在线一区| 亚洲一区二区三区色噜噜| 精品电影一区二区在线| 亚洲美女视频黄频| 亚洲精品美女久久久久99蜜臀| 床上黄色一级片| 国产精品久久视频播放| 久久草成人影院| 脱女人内裤的视频| 在线a可以看的网站| 国产三级中文精品| 一进一出好大好爽视频| 午夜福利免费观看在线| 色播亚洲综合网| aaaaa片日本免费| 日日夜夜操网爽| 女同久久另类99精品国产91| 搡老岳熟女国产| 97碰自拍视频| 男人舔奶头视频| 老熟妇仑乱视频hdxx| 亚洲激情在线av| 88av欧美| 久久香蕉国产精品| 久久香蕉国产精品| 我要搜黄色片| 国产成人精品久久二区二区91| 亚洲精华国产精华精| 久久精品91无色码中文字幕| 久99久视频精品免费| 美女 人体艺术 gogo| 少妇的丰满在线观看| 国产精品 欧美亚洲| 亚洲精品av麻豆狂野| 老汉色∧v一级毛片| 亚洲,欧美精品.| 亚洲国产精品成人综合色| 久久午夜亚洲精品久久| 亚洲一区高清亚洲精品| 成人av在线播放网站| 色av中文字幕| 国产1区2区3区精品| netflix在线观看网站| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩瑟瑟在线播放| 很黄的视频免费| 嫩草影视91久久| 人妻夜夜爽99麻豆av| 国内少妇人妻偷人精品xxx网站 | 欧美精品一区二区大全| 在线观看av片永久免费下载| 亚洲成人中文字幕在线播放| 99久久中文字幕三级久久日本| 波多野结衣高清作品| 嫩草影院精品99| 亚洲va在线va天堂va国产| 少妇人妻精品综合一区二区 | av在线蜜桃| avwww免费| 观看美女的网站| 亚洲欧美精品综合久久99| 99久久九九国产精品国产免费| 精品一区二区免费观看| 一进一出抽搐动态| 色综合色国产| 亚洲内射少妇av| 在线免费观看不下载黄p国产| 伦精品一区二区三区| 一区福利在线观看| 国产麻豆成人av免费视频| 日韩成人伦理影院| 亚洲av.av天堂| 久久午夜亚洲精品久久| 国内精品一区二区在线观看| 长腿黑丝高跟| 免费看av在线观看网站| 久久精品91蜜桃| 亚洲精品粉嫩美女一区| 欧美3d第一页| 日本黄大片高清| 狂野欧美激情性xxxx在线观看| 99久国产av精品| 国产精品一区二区在线观看99 | 寂寞人妻少妇视频99o| 成年女人永久免费观看视频| 成熟少妇高潮喷水视频| 国产精品综合久久久久久久免费| 日本av手机在线免费观看| 日本在线视频免费播放| 久久精品人妻少妇| 亚洲国产精品成人综合色| 亚洲人成网站高清观看| 赤兔流量卡办理| 99riav亚洲国产免费| 久久久国产成人免费| 成人综合一区亚洲| 亚洲精品亚洲一区二区| 长腿黑丝高跟| 国产精品久久久久久精品电影小说 | 亚洲第一电影网av| 看免费成人av毛片| 22中文网久久字幕| 亚洲熟妇中文字幕五十中出| 国模一区二区三区四区视频| 久久久色成人| 国产一级毛片在线| 免费观看精品视频网站| 大香蕉久久网| 午夜福利在线在线| 精品国产三级普通话版| 黄色视频,在线免费观看| 午夜福利视频1000在线观看| 中出人妻视频一区二区| 久久这里有精品视频免费| 成年版毛片免费区| 国产亚洲5aaaaa淫片| 亚洲无线在线观看| 成人性生交大片免费视频hd| 精品久久久噜噜| 国产免费一级a男人的天堂| 国产国拍精品亚洲av在线观看| 久久韩国三级中文字幕| 男人舔奶头视频| 此物有八面人人有两片| 国产精品人妻久久久影院| 男插女下体视频免费在线播放| 日本三级黄在线观看| 欧美日韩国产亚洲二区| 国产乱人视频| 国产精品人妻久久久影院| 老司机影院成人| 又爽又黄无遮挡网站| 成人欧美大片| 成年版毛片免费区| 九九热线精品视视频播放| 人妻久久中文字幕网| 亚洲欧美精品自产自拍| 岛国在线免费视频观看| 97超碰精品成人国产| 亚洲欧洲国产日韩| a级毛色黄片| 夫妻性生交免费视频一级片| 成年女人永久免费观看视频| 亚洲久久久久久中文字幕| 1000部很黄的大片| 97热精品久久久久久| 97超碰精品成人国产| 亚洲综合色惰| 男人狂女人下面高潮的视频| 永久网站在线| 少妇猛男粗大的猛烈进出视频 | 99热网站在线观看| 成年免费大片在线观看| 可以在线观看毛片的网站| 中国美白少妇内射xxxbb| 嫩草影院入口| 一本一本综合久久| 色综合色国产| 免费大片18禁| 看黄色毛片网站| 久久婷婷人人爽人人干人人爱| 熟女人妻精品中文字幕| 色综合色国产| 天美传媒精品一区二区| 国产激情偷乱视频一区二区| 最新中文字幕久久久久| 桃色一区二区三区在线观看| 日韩人妻高清精品专区| 禁无遮挡网站| 国产精品99久久久久久久久| 久久久久久国产a免费观看| 国产精品永久免费网站| 亚洲精品日韩av片在线观看| 久久精品夜色国产| 午夜福利在线观看免费完整高清在 | 熟女人妻精品中文字幕| 亚洲中文字幕一区二区三区有码在线看| 久久综合国产亚洲精品| 久久久a久久爽久久v久久| 日本与韩国留学比较| 在现免费观看毛片| 在线a可以看的网站| 一级毛片久久久久久久久女| 中文字幕av在线有码专区| 国产精品久久久久久精品电影| 美女 人体艺术 gogo| 九九在线视频观看精品| 日日摸夜夜添夜夜爱| 国产午夜精品一二区理论片| 久久久久久久久大av| 久久久久性生活片| 成人毛片a级毛片在线播放| 亚洲欧美日韩卡通动漫| 国产精品99久久久久久久久| 中国国产av一级| av在线老鸭窝| 草草在线视频免费看| 色视频www国产| 国产极品精品免费视频能看的| 精品免费久久久久久久清纯| 成人午夜精彩视频在线观看| 深夜精品福利| 久久人妻av系列| 免费观看a级毛片全部| 老司机影院成人| 国产精品一区二区在线观看99 | 国产午夜福利久久久久久| 嫩草影院新地址| 国产欧美日韩精品一区二区| 日韩三级伦理在线观看| 尤物成人国产欧美一区二区三区| 少妇被粗大猛烈的视频| 一进一出抽搐动态| 午夜福利高清视频| 日本一二三区视频观看| 人妻制服诱惑在线中文字幕| 久久综合国产亚洲精品| 波多野结衣巨乳人妻| 99久久中文字幕三级久久日本| 老司机福利观看| 神马国产精品三级电影在线观看| 在线免费十八禁| 国产高清不卡午夜福利| 亚洲高清免费不卡视频| 三级经典国产精品| 一个人观看的视频www高清免费观看| 久久久欧美国产精品| 久久精品影院6| 久久精品国产亚洲av天美| 欧美日本亚洲视频在线播放| 午夜激情福利司机影院| 亚洲精品自拍成人| 91在线精品国自产拍蜜月| 国产高清三级在线| 久久久久久久久大av| 欧美另类亚洲清纯唯美| 在线观看一区二区三区| 国产精品久久久久久久久免| 久久人人精品亚洲av| 国产精品久久久久久精品电影小说 | 国产高潮美女av| 六月丁香七月| .国产精品久久| 亚洲av熟女| 欧美+日韩+精品| 97在线视频观看| 日韩欧美国产在线观看| 热99在线观看视频| 色尼玛亚洲综合影院| 97人妻精品一区二区三区麻豆| 国产麻豆成人av免费视频| 人妻夜夜爽99麻豆av| 久久久久九九精品影院| 欧美区成人在线视频| 女人十人毛片免费观看3o分钟| 国产伦在线观看视频一区| 日韩强制内射视频| 午夜免费男女啪啪视频观看| 亚洲最大成人av| 欧美潮喷喷水| 床上黄色一级片| 亚洲在线观看片| 美女黄网站色视频| 国产在线精品亚洲第一网站| 在线免费观看的www视频| 免费黄网站久久成人精品| 18禁黄网站禁片免费观看直播| 色视频www国产| 99国产极品粉嫩在线观看| 亚洲欧美中文字幕日韩二区| 边亲边吃奶的免费视频| 麻豆成人午夜福利视频| 国产高潮美女av| 1024手机看黄色片| 91久久精品电影网| 国产精品久久久久久精品电影小说 | 日产精品乱码卡一卡2卡三| 乱系列少妇在线播放| 两个人的视频大全免费| 精品一区二区三区视频在线| 国内少妇人妻偷人精品xxx网站| 好男人在线观看高清免费视频| 女同久久另类99精品国产91| 亚洲欧美清纯卡通| 欧美日韩一区二区视频在线观看视频在线 | 91在线精品国自产拍蜜月| 热99re8久久精品国产| 一级毛片我不卡| 国产精品人妻久久久久久| 国产伦精品一区二区三区视频9| 性色avwww在线观看| 久久久欧美国产精品| 91麻豆精品激情在线观看国产| 蜜桃亚洲精品一区二区三区| 全区人妻精品视频| 久久综合国产亚洲精品| 国产色爽女视频免费观看| 国产av在哪里看| 亚洲欧美清纯卡通| 日韩一区二区视频免费看| 日日啪夜夜撸| 欧美成人一区二区免费高清观看| 国内精品久久久久精免费| 午夜福利在线在线| 两个人视频免费观看高清| 亚洲在线观看片| 中文欧美无线码| 91久久精品国产一区二区成人| 一级av片app| 波野结衣二区三区在线| 国产男人的电影天堂91| 国产精品一二三区在线看| av天堂中文字幕网| 国产一区二区激情短视频| 中出人妻视频一区二区| 亚洲中文字幕日韩| 久久精品综合一区二区三区| 免费黄网站久久成人精品| 久久热精品热| 69人妻影院| 欧美一级a爱片免费观看看| 精品久久久久久久久久免费视频| 一本精品99久久精品77| 中文亚洲av片在线观看爽| 国产一区二区在线av高清观看| 国产色爽女视频免费观看| 亚洲18禁久久av| 联通29元200g的流量卡| 国产一区二区三区在线臀色熟女| 久久精品国产亚洲av涩爱 | 国产午夜精品一二区理论片| 一级毛片久久久久久久久女| 黄色一级大片看看| 一级毛片我不卡| 免费观看在线日韩| 国产美女午夜福利| 欧美日韩在线观看h| 成人二区视频| 国产伦精品一区二区三区视频9| 高清在线视频一区二区三区 | a级毛片免费高清观看在线播放| 美女被艹到高潮喷水动态| 色哟哟·www| 在线国产一区二区在线| 欧美一区二区亚洲| 久久久久国产网址| 最后的刺客免费高清国语| 中文亚洲av片在线观看爽| 男人和女人高潮做爰伦理| 国产精品一区二区性色av| 可以在线观看的亚洲视频| 日本免费一区二区三区高清不卡| 久久精品国产亚洲网站| 看非洲黑人一级黄片| 蜜桃久久精品国产亚洲av| 久久久久久久久大av| 精品久久久久久久久久免费视频| 久久久精品94久久精品| 国产伦在线观看视频一区| 亚洲欧美日韩卡通动漫| 国产精品久久久久久久电影| 麻豆成人午夜福利视频| 亚洲乱码一区二区免费版| 久久精品国产99精品国产亚洲性色| 精品人妻视频免费看| 精品日产1卡2卡| 亚洲欧美成人综合另类久久久 | 三级国产精品欧美在线观看| 可以在线观看的亚洲视频| 国产亚洲5aaaaa淫片| 晚上一个人看的免费电影| 久久99热这里只有精品18| 久久久久久久亚洲中文字幕| 色哟哟·www| 欧美日韩国产亚洲二区| 亚洲婷婷狠狠爱综合网| 边亲边吃奶的免费视频| 日本撒尿小便嘘嘘汇集6| 婷婷色综合大香蕉| 国产高清三级在线| 国产视频首页在线观看| 免费看光身美女| 一个人观看的视频www高清免费观看| 熟女电影av网| 人妻少妇偷人精品九色| 黄片无遮挡物在线观看| 蜜桃久久精品国产亚洲av| 99热只有精品国产| 欧美高清成人免费视频www| 国产黄片视频在线免费观看| 在线a可以看的网站| 麻豆国产av国片精品| 国产三级中文精品| 欧美另类亚洲清纯唯美| 日韩欧美一区二区三区在线观看| 亚洲欧美日韩高清在线视频| 久久久精品欧美日韩精品| 九草在线视频观看| 午夜免费激情av| 不卡一级毛片| 国产69精品久久久久777片| 看非洲黑人一级黄片| 久久人妻av系列| 丰满乱子伦码专区| 少妇熟女aⅴ在线视频| 青春草亚洲视频在线观看| 成年av动漫网址| 色吧在线观看| 麻豆一二三区av精品| 嫩草影院新地址| 美女脱内裤让男人舔精品视频 | kizo精华| 男女视频在线观看网站免费| 日本黄大片高清| 99热这里只有是精品在线观看| avwww免费| 亚洲欧美日韩高清在线视频| 欧美一区二区亚洲| 国产精品一区二区在线观看99 | 九色成人免费人妻av| 国产精品野战在线观看| 偷拍熟女少妇极品色| 国产免费男女视频| 国产成人影院久久av| 一区二区三区高清视频在线| 干丝袜人妻中文字幕| av国产免费在线观看| 男女那种视频在线观看| 激情 狠狠 欧美| 人人妻人人澡人人爽人人夜夜 | 91久久精品国产一区二区三区| 熟妇人妻久久中文字幕3abv| 能在线免费看毛片的网站| 亚洲一区二区三区色噜噜| 免费在线观看成人毛片| 精品国产三级普通话版| 可以在线观看毛片的网站| 嫩草影院新地址| 免费观看在线日韩| 啦啦啦啦在线视频资源| 亚洲av男天堂| 亚洲在线观看片| 免费大片18禁| 免费观看在线日韩| 久久精品国产99精品国产亚洲性色| 91麻豆精品激情在线观看国产| 中文字幕av在线有码专区| 99久久无色码亚洲精品果冻| 两个人的视频大全免费| 精品日产1卡2卡| 亚洲av.av天堂| 日韩在线高清观看一区二区三区| 午夜福利在线观看吧| 久久久久久九九精品二区国产| 国产精品蜜桃在线观看 | av黄色大香蕉| 91久久精品国产一区二区三区| 九草在线视频观看| 亚洲最大成人手机在线| 91av网一区二区| 国产 一区 欧美 日韩| a级毛片a级免费在线| 国产精品.久久久| 午夜爱爱视频在线播放| 亚洲欧洲国产日韩| 小蜜桃在线观看免费完整版高清| 久久6这里有精品| 成年女人看的毛片在线观看| 99热6这里只有精品| 看非洲黑人一级黄片| 久久久a久久爽久久v久久| 欧美性感艳星| 国产午夜福利久久久久久| 成人二区视频| 长腿黑丝高跟| 国产精品.久久久| 国产美女午夜福利| 久久久成人免费电影| 欧美高清成人免费视频www| 成人无遮挡网站| 此物有八面人人有两片| 乱码一卡2卡4卡精品| 偷拍熟女少妇极品色| 你懂的网址亚洲精品在线观看 | 看免费成人av毛片| 人妻少妇偷人精品九色| 久久婷婷人人爽人人干人人爱| 午夜免费男女啪啪视频观看| 一级毛片我不卡| 男女啪啪激烈高潮av片| 国产69精品久久久久777片| 欧美日韩精品成人综合77777| 亚洲精品乱码久久久v下载方式| 国产高清激情床上av| 色尼玛亚洲综合影院| 91麻豆精品激情在线观看国产| 联通29元200g的流量卡| 久久久精品94久久精品| 麻豆成人av视频| 国产成人a∨麻豆精品| 亚洲成人av在线免费| 春色校园在线视频观看| 天堂中文最新版在线下载 | 成人鲁丝片一二三区免费| 26uuu在线亚洲综合色| 日本三级黄在线观看| 黑人高潮一二区| 国产在线男女| 美女xxoo啪啪120秒动态图| 久久久久国产网址| 亚洲精品久久国产高清桃花| 99视频精品全部免费 在线| 亚洲av免费高清在线观看| 亚洲精品亚洲一区二区| 免费看日本二区| 久久久久久伊人网av| 国产精品麻豆人妻色哟哟久久 | 天堂中文最新版在线下载 | 天堂中文最新版在线下载 | 亚洲成人精品中文字幕电影| 日韩大尺度精品在线看网址| 亚洲欧美日韩卡通动漫| 波多野结衣高清无吗| 级片在线观看| 亚洲av第一区精品v没综合| 国内少妇人妻偷人精品xxx网站| 国产av一区在线观看免费| 亚洲av一区综合| 一进一出抽搐gif免费好疼| 国产av一区在线观看免费| 少妇裸体淫交视频免费看高清| 午夜福利在线观看免费完整高清在 | 日韩欧美 国产精品| 久久精品国产鲁丝片午夜精品| 一级毛片久久久久久久久女| 欧美精品一区二区大全| av黄色大香蕉| а√天堂www在线а√下载| 男女啪啪激烈高潮av片| av在线老鸭窝| 卡戴珊不雅视频在线播放| 人人妻人人澡欧美一区二区| 大又大粗又爽又黄少妇毛片口| 好男人视频免费观看在线| 只有这里有精品99| 啦啦啦韩国在线观看视频| 性插视频无遮挡在线免费观看| 国产精品福利在线免费观看| 国产精品电影一区二区三区| 日韩欧美精品免费久久| 欧美日韩精品成人综合77777| 中文字幕av在线有码专区| 久久九九热精品免费| 婷婷六月久久综合丁香|