• <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又大| 亚洲成a人片在线一区二区| xxx96com| 麻豆成人av在线观看| 亚洲精品在线观看二区| videosex国产| 无人区码免费观看不卡| 看片在线看免费视频| 妹子高潮喷水视频| 久久天堂一区二区三区四区| av网站免费在线观看视频| 国产亚洲欧美98| 一级毛片高清免费大全| 黑人巨大精品欧美一区二区蜜桃| 一二三四社区在线视频社区8| 麻豆av在线久日| 欧美成人免费av一区二区三区 | 久久天堂一区二区三区四区| 大香蕉久久成人网| 亚洲精品国产色婷婷电影| 咕卡用的链子| 免费黄频网站在线观看国产| 99精国产麻豆久久婷婷| 在线av久久热| 欧美国产精品va在线观看不卡| 日韩免费av在线播放| 中文字幕人妻丝袜一区二区| 久久九九热精品免费| 在线播放国产精品三级| 精品少妇一区二区三区视频日本电影| 在线永久观看黄色视频| 高清毛片免费观看视频网站 | 欧美久久黑人一区二区| 老汉色∧v一级毛片| 18禁国产床啪视频网站| 国产av又大| 精品一区二区三区av网在线观看| 久久婷婷成人综合色麻豆| 精品人妻在线不人妻| 成在线人永久免费视频| 高清黄色对白视频在线免费看| 欧美乱妇无乱码| 久久香蕉激情| 亚洲精品乱久久久久久| www.精华液| 人妻丰满熟妇av一区二区三区 | 女人被躁到高潮嗷嗷叫费观| 免费黄频网站在线观看国产| 久久精品国产清高在天天线| 捣出白浆h1v1| 老司机影院毛片| 男女之事视频高清在线观看| 国产精品久久视频播放| 99在线人妻在线中文字幕 | 别揉我奶头~嗯~啊~动态视频| 亚洲一卡2卡3卡4卡5卡精品中文| 老司机午夜十八禁免费视频| 成人三级做爰电影| 久久国产精品人妻蜜桃| 午夜免费鲁丝| 搡老熟女国产l中国老女人| 激情在线观看视频在线高清 | 中国美女看黄片| av福利片在线| 99久久综合精品五月天人人| 久久ye,这里只有精品| 看黄色毛片网站| 亚洲欧美一区二区三区久久| 亚洲精品成人av观看孕妇| 久久 成人 亚洲| 一级毛片女人18水好多| 一级毛片女人18水好多| 国产在视频线精品| 国产亚洲av高清不卡| ponron亚洲| 国产不卡一卡二| 亚洲熟妇熟女久久| 亚洲欧美一区二区三区黑人| 久久亚洲精品不卡| 久久国产精品男人的天堂亚洲| 欧美老熟妇乱子伦牲交| 亚洲中文av在线| 亚洲精品美女久久久久99蜜臀| a级片在线免费高清观看视频| 午夜福利乱码中文字幕| 国产一区在线观看成人免费| 亚洲熟妇熟女久久| 久久久久久久精品吃奶| 午夜精品久久久久久毛片777| 大片电影免费在线观看免费| 黄色怎么调成土黄色| 久久精品亚洲熟妇少妇任你| 国产成人欧美| 精品人妻在线不人妻| 亚洲精华国产精华精| www.熟女人妻精品国产| 每晚都被弄得嗷嗷叫到高潮| 性色av乱码一区二区三区2| 久久久久国产一级毛片高清牌| 老司机亚洲免费影院| 久久久精品免费免费高清| 国产精品久久久久久人妻精品电影| 国产精品成人在线| 国产成人一区二区三区免费视频网站| 无限看片的www在线观看| 日本欧美视频一区| 大香蕉久久成人网| 成人永久免费在线观看视频| 国产精品久久久人人做人人爽| 国产黄色免费在线视频| 淫妇啪啪啪对白视频| 丝袜在线中文字幕| 国产有黄有色有爽视频| 亚洲欧洲精品一区二区精品久久久| 麻豆乱淫一区二区| 国产成人精品在线电影| www.精华液| 久久久精品区二区三区| 亚洲性夜色夜夜综合| 久久久久视频综合| 激情视频va一区二区三区| 国产一区二区激情短视频| 国产精品秋霞免费鲁丝片| 高清欧美精品videossex| 国产主播在线观看一区二区| 午夜精品在线福利| 亚洲一区二区三区不卡视频| 亚洲视频免费观看视频| 精品久久久久久电影网| 超碰成人久久| 久久天堂一区二区三区四区| 久久香蕉激情| 超碰成人久久| 日本黄色日本黄色录像| 久久久国产成人精品二区 | 999久久久国产精品视频| 极品教师在线免费播放| 精品第一国产精品| 在线观看一区二区三区激情| 久久精品亚洲av国产电影网| 美女福利国产在线| 日韩欧美一区视频在线观看| 欧美日韩视频精品一区| 超碰成人久久| 国产午夜精品久久久久久| aaaaa片日本免费| 国产精品乱码一区二三区的特点 | 中文字幕制服av| 天堂动漫精品| 99国产综合亚洲精品| 国产精品偷伦视频观看了| 多毛熟女@视频| 动漫黄色视频在线观看| 少妇裸体淫交视频免费看高清 | 午夜成年电影在线免费观看| 国产深夜福利视频在线观看| 亚洲国产精品一区二区三区在线| 黑人欧美特级aaaaaa片| 91字幕亚洲| 中文字幕另类日韩欧美亚洲嫩草| 国产在线精品亚洲第一网站| tube8黄色片| 午夜福利在线观看吧| 午夜免费鲁丝| 国产主播在线观看一区二区| 亚洲欧美激情综合另类| 中文亚洲av片在线观看爽 | 老司机亚洲免费影院| 亚洲第一青青草原| 久久久精品国产亚洲av高清涩受| 亚洲欧美精品综合一区二区三区| 露出奶头的视频| 亚洲综合色网址| 午夜两性在线视频| 中文字幕人妻丝袜制服| 国内久久婷婷六月综合欲色啪| 亚洲精品中文字幕在线视频| 满18在线观看网站| 一级毛片精品| 午夜福利一区二区在线看| 免费观看a级毛片全部| 亚洲一区二区三区不卡视频| 久久久久精品人妻al黑| 欧美另类亚洲清纯唯美| 亚洲少妇的诱惑av| 久久九九热精品免费| av电影中文网址| 国产免费av片在线观看野外av| 久久久国产成人免费| 午夜福利影视在线免费观看| 亚洲精品美女久久久久99蜜臀| 亚洲精品乱久久久久久| 999久久久精品免费观看国产| 国产一区二区三区在线臀色熟女 | 一边摸一边做爽爽视频免费| 国产精品一区二区在线观看99| 国产精品欧美亚洲77777| 老汉色∧v一级毛片| 性色av乱码一区二区三区2| 首页视频小说图片口味搜索| 亚洲午夜理论影院| 激情视频va一区二区三区| 国产又色又爽无遮挡免费看| 大型av网站在线播放| 精品少妇久久久久久888优播| 黑人巨大精品欧美一区二区蜜桃| videosex国产| 在线观看日韩欧美| 成人影院久久| 新久久久久国产一级毛片| 一级,二级,三级黄色视频| 日日爽夜夜爽网站| 丰满的人妻完整版| 两人在一起打扑克的视频| 亚洲欧美精品综合一区二区三区| 欧美日韩乱码在线| 高潮久久久久久久久久久不卡| 精品亚洲成a人片在线观看| 一区二区三区精品91| 无限看片的www在线观看| 日本精品一区二区三区蜜桃| 国产精品久久电影中文字幕 | 90打野战视频偷拍视频| 久久久久精品国产欧美久久久| av天堂在线播放| 1024视频免费在线观看| 中文字幕精品免费在线观看视频| 国产欧美日韩精品亚洲av| 久久草成人影院| 国产精品久久久久久精品古装| 国产成人精品在线电影| 两人在一起打扑克的视频| 99精品在免费线老司机午夜| 欧美成人午夜精品| 1024香蕉在线观看| 欧美日韩成人在线一区二区| 欧美国产精品一级二级三级| 成年人午夜在线观看视频| 99re在线观看精品视频| 18禁美女被吸乳视频| 99久久综合精品五月天人人| 老司机亚洲免费影院| 无遮挡黄片免费观看| 久久人妻福利社区极品人妻图片| 日本撒尿小便嘘嘘汇集6| 色播在线永久视频| 国产深夜福利视频在线观看| 亚洲性夜色夜夜综合| 人妻久久中文字幕网| 亚洲熟女毛片儿| 黄色视频不卡| 精品第一国产精品| 在线永久观看黄色视频| 国产男女内射视频| 韩国精品一区二区三区| 嫁个100分男人电影在线观看| 99精品在免费线老司机午夜| 人人妻,人人澡人人爽秒播| 亚洲av电影在线进入| 多毛熟女@视频| 国产一区二区三区视频了| 99re6热这里在线精品视频| 国产欧美日韩综合在线一区二区| 一级a爱视频在线免费观看| 精品国产一区二区三区久久久樱花| 精品第一国产精品| 一区福利在线观看| 91麻豆av在线| 日韩成人在线观看一区二区三区| 日本wwww免费看| 欧美日韩黄片免| 男女下面插进去视频免费观看| 美女福利国产在线| 亚洲久久久国产精品| 亚洲国产中文字幕在线视频| 精品久久久久久,| 女警被强在线播放| 欧美日韩黄片免| 亚洲五月天丁香| 搡老熟女国产l中国老女人| 捣出白浆h1v1| 高潮久久久久久久久久久不卡| 两个人免费观看高清视频| 99热网站在线观看| 午夜免费成人在线视频| 日本一区二区免费在线视频| 999精品在线视频| 母亲3免费完整高清在线观看| 91在线观看av| 国产免费现黄频在线看| 脱女人内裤的视频| 深夜精品福利| 91老司机精品| 欧美激情 高清一区二区三区| 国产男女内射视频| tube8黄色片| x7x7x7水蜜桃| 国产在线一区二区三区精| 麻豆国产av国片精品| 亚洲精品一卡2卡三卡4卡5卡| 午夜亚洲福利在线播放| 欧美日韩国产mv在线观看视频| 99re6热这里在线精品视频| 伦理电影免费视频| 日韩欧美一区视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产不卡一卡二| 国产精品亚洲av一区麻豆| 纯流量卡能插随身wifi吗| 18禁裸乳无遮挡免费网站照片 | 黄色视频不卡| 精品国内亚洲2022精品成人 | 欧美日韩亚洲国产一区二区在线观看 | 欧美最黄视频在线播放免费 | 不卡av一区二区三区| 亚洲精品久久成人aⅴ小说| 国产精品成人在线| 黄色毛片三级朝国网站| 亚洲伊人色综图| 亚洲色图 男人天堂 中文字幕| 国产精品亚洲一级av第二区| 一级片免费观看大全| www日本在线高清视频| 久久久国产精品麻豆| 一级毛片女人18水好多| 日韩三级视频一区二区三区| 精品少妇久久久久久888优播| 欧美老熟妇乱子伦牲交| 最近最新中文字幕大全电影3 | aaaaa片日本免费| 18禁观看日本| 99精品久久久久人妻精品| 不卡一级毛片| 久久午夜亚洲精品久久| 飞空精品影院首页| 国产精品美女特级片免费视频播放器 | 欧美精品av麻豆av| 在线观看免费视频网站a站| 国产单亲对白刺激| 国产xxxxx性猛交| 夜夜夜夜夜久久久久| 亚洲 国产 在线| 久99久视频精品免费| 欧美 日韩 精品 国产| 成年人黄色毛片网站| 国产亚洲一区二区精品| 国产精品成人在线| 悠悠久久av| 如日韩欧美国产精品一区二区三区| 午夜激情av网站| 精品乱码久久久久久99久播| 欧美成狂野欧美在线观看| 最近最新中文字幕大全电影3 | 日韩熟女老妇一区二区性免费视频| 精品久久久久久久久久免费视频 | 久久中文字幕人妻熟女| 国产成人精品久久二区二区免费| 黄色视频,在线免费观看| 黑人操中国人逼视频| 日韩视频一区二区在线观看| 精品久久久精品久久久| 亚洲精品中文字幕一二三四区| 免费在线观看视频国产中文字幕亚洲| 欧美黑人欧美精品刺激| 99久久精品国产亚洲精品| 欧美精品高潮呻吟av久久| 国产一区二区三区在线臀色熟女 | 欧美亚洲日本最大视频资源| 一夜夜www| 亚洲欧美日韩高清在线视频| 大型黄色视频在线免费观看| 日本a在线网址| xxxhd国产人妻xxx| 人妻 亚洲 视频| 男人操女人黄网站| 亚洲国产欧美日韩在线播放| 国产精品美女特级片免费视频播放器 | 国产在线精品亚洲第一网站| 在线天堂中文资源库| 一级毛片精品| 大陆偷拍与自拍| 亚洲专区字幕在线| av网站免费在线观看视频| 欧美性长视频在线观看| 岛国毛片在线播放| 欧美日韩亚洲国产一区二区在线观看 | 黄色女人牲交| 久久热在线av| 在线免费观看的www视频| 大码成人一级视频| 亚洲精品在线美女| 黄频高清免费视频| 亚洲色图综合在线观看| 亚洲av日韩精品久久久久久密| 日日摸夜夜添夜夜添小说| 久久精品国产亚洲av香蕉五月 | 亚洲伊人色综图| 亚洲色图 男人天堂 中文字幕| 老司机靠b影院| 女警被强在线播放| 久久国产乱子伦精品免费另类| 可以免费在线观看a视频的电影网站| 在线观看一区二区三区激情| 国产精品 欧美亚洲| 国产av一区二区精品久久| 最新美女视频免费是黄的| 亚洲成av片中文字幕在线观看| av国产精品久久久久影院| 亚洲精品粉嫩美女一区| 搡老熟女国产l中国老女人| 性色av乱码一区二区三区2| 一边摸一边抽搐一进一出视频| 亚洲黑人精品在线| 丰满饥渴人妻一区二区三| 制服诱惑二区| 大型黄色视频在线免费观看| 国产成人欧美在线观看 | 两个人看的免费小视频| 国产三级黄色录像| 一本大道久久a久久精品| 99国产精品一区二区蜜桃av | 女警被强在线播放| 久久国产乱子伦精品免费另类| 成熟少妇高潮喷水视频| 亚洲第一青青草原| 一级毛片高清免费大全| 亚洲欧美一区二区三区黑人| 亚洲精品成人av观看孕妇| 制服诱惑二区| 亚洲少妇的诱惑av| 亚洲国产中文字幕在线视频| 一本一本久久a久久精品综合妖精| 国产精品电影一区二区三区 | 侵犯人妻中文字幕一二三四区| 又大又爽又粗| 成人av一区二区三区在线看| 欧美丝袜亚洲另类 | 亚洲伊人色综图| 老司机午夜福利在线观看视频| 一区二区日韩欧美中文字幕| 一级毛片女人18水好多| 99国产综合亚洲精品| 成年人午夜在线观看视频| 18在线观看网站| 悠悠久久av| 精品熟女少妇八av免费久了| 不卡一级毛片| 18禁裸乳无遮挡免费网站照片 | 国产午夜精品久久久久久| 午夜91福利影院| 99国产精品一区二区蜜桃av | 精品无人区乱码1区二区| 国产成人系列免费观看| 国产精品一区二区在线不卡| 五月开心婷婷网| 国产精品久久久久成人av| 啪啪无遮挡十八禁网站| 精品国产超薄肉色丝袜足j| 国产成人av教育| 免费少妇av软件| www.熟女人妻精品国产| av中文乱码字幕在线| 日本wwww免费看| 免费少妇av软件| 亚洲欧洲精品一区二区精品久久久| 一二三四社区在线视频社区8| 9191精品国产免费久久| 国产成人精品久久二区二区91| 老司机深夜福利视频在线观看| 露出奶头的视频| 中文字幕高清在线视频| 国产91精品成人一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 极品教师在线免费播放| 国产精品免费视频内射| 国产精品美女特级片免费视频播放器 | 精品国产国语对白av| 满18在线观看网站| 久久久久久免费高清国产稀缺| 日本一区二区免费在线视频| 精品一区二区三卡| 午夜免费鲁丝| 俄罗斯特黄特色一大片| 亚洲精品成人av观看孕妇| 免费看十八禁软件| 亚洲专区国产一区二区| 夜夜爽天天搞| 国产亚洲精品一区二区www | 久久人人爽av亚洲精品天堂| 久久天堂一区二区三区四区| 老司机福利观看| 国产又爽黄色视频| 很黄的视频免费| 国产精品综合久久久久久久免费 | 777久久人妻少妇嫩草av网站| 天堂中文最新版在线下载| 99久久99久久久精品蜜桃| 午夜亚洲福利在线播放| 在线观看www视频免费| 国产精品成人在线| 欧美日韩黄片免| 制服诱惑二区| 丝袜在线中文字幕| 性色av乱码一区二区三区2| 色播在线永久视频| 日韩熟女老妇一区二区性免费视频| 日日爽夜夜爽网站| 99精品欧美一区二区三区四区| 国产精品久久电影中文字幕 | 久久久精品免费免费高清| 91麻豆av在线| 国产高清视频在线播放一区| 中亚洲国语对白在线视频| 国产视频一区二区在线看| 美女 人体艺术 gogo| 精品国产一区二区三区四区第35| 9色porny在线观看| 成在线人永久免费视频| 一级a爱视频在线免费观看| 国产国语露脸激情在线看| 巨乳人妻的诱惑在线观看| 高潮久久久久久久久久久不卡| 免费观看人在逋| 国产亚洲一区二区精品| 国产成人欧美在线观看 | 天天躁日日躁夜夜躁夜夜| 最近最新免费中文字幕在线| 国产深夜福利视频在线观看| 高清在线国产一区| 免费高清在线观看日韩| 欧美成人午夜精品| 最新在线观看一区二区三区| 精品国产美女av久久久久小说| 真人做人爱边吃奶动态| 一级片免费观看大全| 伊人久久大香线蕉亚洲五| 国产精品免费大片| 久久久国产精品麻豆| 一级a爱视频在线免费观看| 美女 人体艺术 gogo| 日韩欧美一区二区三区在线观看 | 天天躁夜夜躁狠狠躁躁| 女警被强在线播放| 精品亚洲成a人片在线观看| 婷婷丁香在线五月| 少妇被粗大的猛进出69影院| 国产片内射在线| 黑丝袜美女国产一区| 少妇 在线观看| 极品人妻少妇av视频| 亚洲精品一二三| 国产在视频线精品| 成人黄色视频免费在线看| 国产一区有黄有色的免费视频| 人妻丰满熟妇av一区二区三区 | 黑人操中国人逼视频| 男人的好看免费观看在线视频 | 精品国产一区二区三区久久久樱花| 亚洲欧美激情在线| 国产97色在线日韩免费| www.999成人在线观看| av不卡在线播放| tocl精华| 国产精品香港三级国产av潘金莲| 最近最新中文字幕大全电影3 | 真人做人爱边吃奶动态| 亚洲av日韩精品久久久久久密| 欧美性长视频在线观看| 一夜夜www| 亚洲熟女毛片儿| 久久这里只有精品19| 国产不卡一卡二| 久久久久久久久免费视频了| 亚洲精品国产精品久久久不卡| 看黄色毛片网站| 亚洲欧美一区二区三区黑人| 999久久久国产精品视频| 国产片内射在线| 中文字幕人妻丝袜制服| 欧美 日韩 精品 国产| 熟女少妇亚洲综合色aaa.| 亚洲av电影在线进入| 精品国产超薄肉色丝袜足j| 大片电影免费在线观看免费| 国产成人精品在线电影| 久久久精品区二区三区| 校园春色视频在线观看| 精品久久久久久久毛片微露脸| 日本黄色视频三级网站网址 | 美女高潮喷水抽搐中文字幕| 超碰成人久久| 人人妻人人澡人人看| 久久天堂一区二区三区四区| av不卡在线播放| 视频区图区小说| 欧美黄色片欧美黄色片| 欧美不卡视频在线免费观看 | 亚洲综合色网址| 18禁国产床啪视频网站| 91精品国产国语对白视频| 男女之事视频高清在线观看| 午夜两性在线视频| 99久久国产精品久久久| 久久久国产欧美日韩av| 在线播放国产精品三级| av有码第一页| 色老头精品视频在线观看| 久久九九热精品免费| 精品人妻1区二区| 熟女少妇亚洲综合色aaa.| 美女高潮喷水抽搐中文字幕| 午夜91福利影院| 中国美女看黄片| 男人舔女人的私密视频|