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

    齊口裂腹魚上溯洄游過程中的感官圈解譯研究

    2022-08-23 10:11:16鄧月曦唐錫良嚴忠鑾王小明盧晶瑩安瑞冬
    中國環(huán)境科學(xué) 2022年8期
    關(guān)鍵詞:腹魚轉(zhuǎn)角魚類

    鄧月曦,唐錫良,嚴忠鑾,王小明,盧晶瑩,安瑞冬,李 嘉*

    齊口裂腹魚上溯洄游過程中的感官圈解譯研究

    鄧月曦1,唐錫良2,嚴忠鑾2,王小明2,盧晶瑩3,安瑞冬1,李 嘉1*

    (1.四川大學(xué)水利水電學(xué)院,水力學(xué)與山區(qū)河流開發(fā)保護國家重點實驗室,四川 成都 610065;2.中國三峽建工(集團)有限公司,四川 成都 610041;3.中國長江三峽集團有限公司,湖北 武漢 430010)

    選取西南地區(qū)典型魚類齊口裂腹魚(,(32.32±2.73)cm)作為實驗對象,測試魚類在上溯過程中應(yīng)對不同水流環(huán)境時的游泳行為,解譯不同大小網(wǎng)格尺寸下魚類的游泳速度、轉(zhuǎn)向次數(shù)以及轉(zhuǎn)角和轉(zhuǎn)角差與真實游泳軌跡的匹配程度.結(jié)果表明:采用小于或等于0.30m尺寸的網(wǎng)格時,分析得出的魚類上溯運動平均游泳速度、轉(zhuǎn)向次數(shù)以及轉(zhuǎn)角和轉(zhuǎn)角差均與實際情況無顯著性差異,另外魚類更偏向于在0.20~0.40m/s的流速區(qū)間內(nèi)以和水平方向呈-49o~49o的運動搜索范圍進行上溯運動.

    網(wǎng)格尺寸;生態(tài)水流;魚類行為;感官圈

    解析魚類上溯洄游過程中應(yīng)對復(fù)雜流場的行為響應(yīng),如趨流行為及逃逸行為等[1-2],是優(yōu)化過魚設(shè)施設(shè)計[3]、保護洄游性魚類及恢復(fù)河流天然屬性的重要研究內(nèi)容[4-5].基于魚類個體行為學(xué)與水動力學(xué)耦合的數(shù)學(xué)模型為模擬洄游性魚類通過自然水流障礙或水工建筑物提供了理論支撐,同時也為生物粒子模擬構(gòu)建提供了數(shù)據(jù)支持[6].魚的側(cè)線系統(tǒng)能夠獲取魚體周圍的水流信號,評估游泳區(qū)域內(nèi)的水動力條件,進而作出相應(yīng)的行為反應(yīng)[7]. Goodwin等[8]基于此提出了感官圈的概念,即魚類的上溯洄游軌跡模擬定義為以感官圈的形式進行搜索上溯.感官圈特征主要包括網(wǎng)格尺寸劃分、水動力條件趨向性、感官搜索范圍、感官搜索角度4個部分,魚類在合適的網(wǎng)格單元內(nèi)以感官搜索范圍在搜索角度范圍內(nèi)對水動力條件進行分析,選擇偏好的水力條件進行上溯.其中感官搜索范圍(SQD)選擇生物響應(yīng)范圍和CFD模型搜索范圍中的較大值,它與魚體體長呈正相關(guān)關(guān)系,體長越長的魚能夠檢索到更大范圍內(nèi)的水動力學(xué)條件[9].Gao等[10]建立了歐拉-拉格朗日魚類運動模型,但該模型僅考慮湍流動能作為水力刺激因子運用感官圈進行魚類的上溯模擬.Tan等[11]進一步建立考慮湍流動能、應(yīng)變率和速度等多種水力刺激的魚類運動軌跡模型,模擬獲得鰱魚的水力刺激范圍.Zielinski等[12]考慮了感官搜索角度的范圍,假設(shè)魚只被刺激向上游動,即魚不回溯,以魚位置上游180°的半圓弧內(nèi)進行網(wǎng)格節(jié)點的搜索,以減少模擬計算量.Zhu等[13]將感官圈的概念運用于生物顆粒模型來預(yù)測大壩下魚的上溯軌跡,生物顆粒模型將魚視為顆粒,并將流量、速度梯度和湍流動能作為回溯行為的條件,利用控制方程簡化魚的行為,開發(fā)魚的可追溯預(yù)測程序.

    然而上述方法并沒有考慮到網(wǎng)格尺寸的大小對魚類上溯運動感官圈解譯的影響,網(wǎng)格尺寸過大時無法準確解譯魚類有效游泳行為信息,導(dǎo)致模擬結(jié)果與真實行為差異較大,網(wǎng)格尺寸過小時雖能捕捉豐富行為信息,但在大尺度空間中使用時會使得模型計算時間長,效率低下,合理的網(wǎng)格尺寸大小有利于精準、有效的解譯魚類上溯行為.另一方面,目前在模擬魚類上溯運動時常對魚類感官圈定義為360°整圓方向,但魚類的實際運動有可能是呈現(xiàn)的扇形搜索范圍,進一步確定搜索范圍的區(qū)間可以減少預(yù)測魚類軌跡時的計算量.據(jù)此本文研究了不同網(wǎng)格尺寸劃分對齊口裂腹魚上溯軌跡模擬準確程度的影響,選定合適的網(wǎng)格尺寸來反映魚類上溯行為,此外,對齊口裂腹魚在上溯運動過程中的響應(yīng)流速、感官搜索角度的范圍進行進一步探索.本研究為魚類上溯應(yīng)對湍流環(huán)境時的運動軌跡模擬網(wǎng)格尺寸選取提供了一種研究思路,同時也為生物粒子模擬構(gòu)建提供了數(shù)據(jù)支撐.

    1 材料與方法

    1.1 研究對象

    齊口裂腹魚()屬于鯉科、裂腹魚屬的底棲魚類,喜在較低的水溫環(huán)境中活動,一般生活于急緩流交界處,有短距離的河川生殖洄游需求[14-15].實驗所用齊口裂腹魚來自四川水產(chǎn)研究所于2020年10月網(wǎng)捕的成熟魚種,平均體長(32.32±2.73)cm,體重(304.21±77.15)g.實驗過程中隨機挑選生理活性較好的齊口裂腹魚作為實驗對象,為避免魚類環(huán)境適應(yīng)性影響實驗結(jié)果的可靠性,實驗魚不重復(fù)使用.

    1.2 實驗裝置

    本研究在四川大學(xué)水力學(xué)與山區(qū)河流開發(fā)保護國家重點實驗室自制的大型開放式淺水循環(huán)水槽(圖1)內(nèi)進行.水槽由鋼化玻璃和鋼架構(gòu)造,長20.50m、寬4.00m、高1.00m.循環(huán)裝置由4個不同的功能部分組成:上游蓄水池、適應(yīng)區(qū)域、下游蓄水池和循環(huán)供水渠道.水槽的邊墻以及底板均為光滑直線型,以減少流動摩擦和邊界層效應(yīng).

    在上游3個等寬(1.00m)入口處通過進口格柵將游泳測試區(qū)域水流分成低、中和高流速區(qū)3種流速大小不同的區(qū)域,模擬不同水流環(huán)境條件下的魚類上溯運動.

    圖1 魚類上溯行為研究實驗裝置幾何示意

    1.3 實驗方法

    每次試驗開始時,在適應(yīng)區(qū)域放入1條健康的個體魚,以適應(yīng)水流溫度和流速條件.適應(yīng)0.50h后,抬升啟閉攔網(wǎng),讓試驗魚在游泳測試區(qū)域自主游泳,同時使用4臺高清攝像機(WIM SkyStar,30fps)連續(xù)記錄魚類上溯洄游游泳行為.試驗結(jié)束條件為當(dāng)測試魚通過整個游泳測試區(qū)域并最終進入上游進口,或者測試魚在1.00h內(nèi)從未離開過適應(yīng)區(qū)域,測試結(jié)束.

    圖2 實驗方法示意

    魚的運動方向定義為深色和淺色兩種狀態(tài)方向,若運動方向出現(xiàn)顏色的轉(zhuǎn)變,則判定為轉(zhuǎn)向

    本研究選擇50尾齊口裂腹魚進行實驗,實驗期間禁止在水槽周圍活動,避免聲音對實驗魚的干擾.后期通過LoggerPro視頻處理軟件以1/30(s)的幀率記錄魚的位置坐標數(shù)據(jù).將魚的軌跡點分別導(dǎo)入網(wǎng)格大小為0.20m×0.20m、0.25m×0.25m、0.30m× 0.30m、0.40m×0.40m、0.50m×0.50m、1.00m×1.00m的網(wǎng)格中進行分析.如圖2所示,通過對比魚類上溯運動平均游泳速度、轉(zhuǎn)向次數(shù)以及轉(zhuǎn)角和轉(zhuǎn)角差得到合理的齊口裂腹魚上溯運動網(wǎng)格分析尺寸.

    此外,實驗水槽的水流速度通過3D聲學(xué)多普勒測速儀(ADV)進行測量,研究齊口裂腹魚在上溯運動過程中對流速的偏好性,同時探究上溯搜索角度有無特定范圍.

    1.4 數(shù)據(jù)分析

    為研究不同網(wǎng)格尺寸下魚類游泳行為信息的連續(xù)及缺失程度,本研究采用游泳軌跡連續(xù)性、魚類平均游泳速度、上溯軌跡轉(zhuǎn)向、上溯軌跡轉(zhuǎn)角及轉(zhuǎn)角差5類指標對網(wǎng)格有效性進行評估.其中相鄰軌跡點之間的魚類平均速度如式(1)所示:

    圖2c中給出了轉(zhuǎn)角的定義,轉(zhuǎn)角是以橫坐標的正方向作為參考方向的,為利于研究統(tǒng)計,轉(zhuǎn)角范圍以[0o,360o]表示;轉(zhuǎn)角差則定義為相鄰兩次轉(zhuǎn)角之間的差值,其中順時針方向為正(+),稱為右旋,反之為負(-),稱為左旋.

    轉(zhuǎn)角差Δ的計算方法如下式:

    采用SPSS 20.0進行相關(guān)數(shù)據(jù)統(tǒng)計分析,不同網(wǎng)格尺寸下魚的平均游泳速度采用單因素方差(ANOVA)分析.統(tǒng)計值均以平均值±標準誤(Mean±SE)表示,顯著水平為<0.05.

    2 結(jié)果與分析

    2.1 確定適宜網(wǎng)格尺寸

    2.1.1 軌跡連續(xù)性 對魚類運動的研究表明,當(dāng)測試魚進入試驗水體環(huán)境時,它很快就能找到自己喜好的水深(或動水壓力),然后在同一水層中游動.事實上,Rodriguez等[16]記錄的真實魚的軌跡也是在平面上呈現(xiàn)的,即二維上的,因此研究可忽略垂直上的變化.本文采用二維水流速度進行分析(圖3),LoggerPro處理得到的單尾齊口裂腹魚游泳速度區(qū)間為=1.47m/s左右,參考庫朗數(shù)(o)的定義[17],當(dāng)時間步長為1秒時,網(wǎng)格的大小不應(yīng)大于×Δ,即不應(yīng)大于1.47m.因此確定研究網(wǎng)格尺寸為0.20,0.25,0.30,0.40,0.50,1.00m.

    圖3 魚類上溯軌跡圖及游泳速度置信區(qū)間圖

    i>X/表示魚類運動坐標與測試長度范圍(13m)的比值

    將齊口裂腹魚的軌跡點導(dǎo)入不同尺寸大小的網(wǎng)格進行篩選,導(dǎo)入軌跡點圓圈的直徑為魚的平均體長32.32cm,采用計算流體力學(xué)的定義,圈內(nèi)不包含有網(wǎng)格節(jié)點的即視為不能被網(wǎng)格識別,需刪除該軌跡點.最終得到的魚類上溯軌跡點如圖4所示,當(dāng)網(wǎng)格尺寸小于0.30m時,軌跡點連續(xù);當(dāng)尺寸為0.40m時,軌跡點圖出現(xiàn)部分斷裂;而當(dāng)尺寸大于0.40m時,圖像無法完整地描述魚類的上溯運動過程.

    圖4 不同網(wǎng)格尺寸下單尾魚游泳軌跡解譯圖

    D表示網(wǎng)格尺寸

    2.1.2 魚類平均速度對比 通過篩選出的各點以及記錄的間隔時間計算魚的平均游泳速度,結(jié)果如圖5所示,采用小于0.30m的網(wǎng)格尺寸分析齊口裂腹魚的平均游泳速度與不進行網(wǎng)格劃分時沒有顯著性差異(>0.05);采用0.40,0.50,1.00m的網(wǎng)格尺寸計算結(jié)果彼此之間均具有顯著性差異(0.05),且與真實數(shù)據(jù)也具有顯著性差異(0.05).表明當(dāng)網(wǎng)格尺寸大于0.30m時,對齊口裂腹魚的上溯運動描述已經(jīng)出現(xiàn)偏差.

    圖5 50尾魚的游泳速度顯著性差異分析

    2.1.3 上溯軌跡轉(zhuǎn)向?qū)Ρ?以1尾魚為例,將篩選出的魚的運動軌跡點依次連線,當(dāng)網(wǎng)格尺寸為0.25m時(圖6a),齊口裂腹魚在上溯過程中共計轉(zhuǎn)向9次,而網(wǎng)格尺寸增大到0.50m時(圖6b),魚共計轉(zhuǎn)向3次,當(dāng)轉(zhuǎn)向次數(shù)明顯減少時,魚的運動軌跡不連貫,軌跡點不能正確描述上溯運動狀態(tài).統(tǒng)計50尾魚上溯過程的轉(zhuǎn)向次數(shù),如圖6c所示,當(dāng)網(wǎng)格尺寸小于0.30m時,魚的運動轉(zhuǎn)向次數(shù)與實際相似,但網(wǎng)格尺寸大于0.30m時轉(zhuǎn)向次數(shù)明顯減少、異常值增多,與實際情況不符,不能正確反應(yīng)魚的上溯運動軌跡特征.

    2.1.4 上溯軌跡轉(zhuǎn)角及轉(zhuǎn)角差對比 由圖7a、b可見,以1尾魚為例,當(dāng)網(wǎng)格尺寸由0.25m變?yōu)?.50m時,齊口裂腹魚上溯過程中的轉(zhuǎn)角次數(shù)明顯減小,轉(zhuǎn)角大小也發(fā)生較大轉(zhuǎn)變,因此可以通過統(tǒng)計魚類轉(zhuǎn)角的次數(shù)分布和數(shù)目,判別上溯運動軌跡的真實性.

    圖6 魚類上溯轉(zhuǎn)向統(tǒng)計

    如圖8所示,魚類轉(zhuǎn)角的分布規(guī)律與網(wǎng)格尺寸的大小無關(guān),上溯軌跡左旋與右旋大約平分,右旋略大于左旋.網(wǎng)格尺寸為0.20和0.25m時齊口裂腹魚上溯運動的轉(zhuǎn)角、轉(zhuǎn)角差次數(shù)與真實數(shù)據(jù)情況極為相似,網(wǎng)格尺寸為0.30m時轉(zhuǎn)角、轉(zhuǎn)角差次數(shù)略小于真實情況,而網(wǎng)格尺寸為0.40和0.50m時轉(zhuǎn)角、轉(zhuǎn)角差次數(shù)減少了二分之一,1.00m時轉(zhuǎn)角、轉(zhuǎn)角差次數(shù)更是減少了四分之一,不能準確描述到上溯運動軌跡.

    圖7 單尾魚轉(zhuǎn)角示意

    2.2 流速趨向性

    上述研究表明,0.30m是網(wǎng)格尺寸分析的合理界限,因此選取此網(wǎng)格下的軌跡點,以1s為時間步長,分析魚類上溯運動的流速趨向性以及上溯運動角度搜索范圍.

    本文使用聲學(xué)多普勒測速儀提取流場的流速分布,將50尾齊口裂腹魚的軌跡點對應(yīng)水流流速通過雙線性插值法計算得出,統(tǒng)計不同流速對應(yīng)的出現(xiàn)相對頻率如圖9所示.齊口裂腹魚對于流場流速的選擇具有一定的趨向性,它更偏向于0.20~0.40m/s的流速區(qū)間內(nèi)進行上溯運動.

    圖9 50尾魚流速趨向

    2.3 感官上溯搜索角度范圍

    以10°作為統(tǒng)計離散程度,統(tǒng)計魚類上溯運動時不同轉(zhuǎn)角范圍的次數(shù),如圖10所示,50尾魚共計轉(zhuǎn)角4423次,取75%的置信區(qū)間,魚的搜索角度范圍是-49°~49°.因此,認為齊口裂腹魚偏好以與水平方向呈98°的扇形角度搜索范圍進行上溯運動的.

    圖10 50尾魚感官上溯搜索角度

    3 討論及展望

    魚類的運動在時空上具有離散個體的差異性,其游泳軌跡的模擬研究大致可分為兩個過程,首先是感知并評估感官圈內(nèi)的各種水動力學(xué)因子[19],然后根據(jù)反饋的信息作出上溯方向的決策.如Liao等[20]根據(jù)細鱗裂腹魚、圓口銅魚和長薄鰍3種魚類游泳速度與游泳行為之間的關(guān)系,提出了6種對應(yīng)不同游泳行為策略的流速范圍區(qū)間,基于此構(gòu)建了與魚類游泳特征相對應(yīng)的三維游泳軌跡模型;蔣永強等[21]綜合分析流速、湍流動能和湍動耗散率對兩種不同體長鮭魚上溯的影響,并依據(jù)上溯魚類游泳速度與周圍水力因子的相關(guān)性擬定函數(shù)關(guān)系,從而定量計算魚類上溯軌跡,并發(fā)現(xiàn)湍流動能的影響最為顯著,湍動能耗散率影響最小;Tan等[11]利用基于個體的IBM模型來模擬鰱魚運動軌跡,該模型綜合了魚類響應(yīng)的水力刺激因子(湍流動能、速度和應(yīng)變率)以及魚類個體運動的規(guī)則;Arenas等[22]則進一步指出鮭魚可以響應(yīng)除流速外的更多水力刺激因子,如溫度、偏振光、地球磁場、存在的障礙等.本研究齊口裂腹魚上溯應(yīng)對復(fù)雜紊流時網(wǎng)格尺寸劃分對軌跡解譯的影響結(jié)果表明當(dāng)時間間隔為1s,采用小于或等于0.30m尺寸的網(wǎng)格時,分析得出的魚類上溯運動平均游泳速度、轉(zhuǎn)向次數(shù)以及轉(zhuǎn)角和轉(zhuǎn)角差均與實際情況無顯著性差異,而網(wǎng)格尺寸大于0.30m時魚類游泳行為與實際運動軌跡有顯著性差異,因此0.30m的網(wǎng)格尺寸可以較合適地反應(yīng)魚類上溯行為,同時這一網(wǎng)格尺寸大小也最接近魚的體長.本研究以體長范圍為(32.32±2.73)cm,體重范圍為(304.21±77.15)g的齊口裂腹魚為研究對象得出了與體長范圍相關(guān)的適宜感官圈范圍.由于同種魚處于不同生長階段時具有不同游泳行為策略,如轉(zhuǎn)向角度及游泳能力等具有差異性[23-24],故研究結(jié)果可為體長相近范圍內(nèi)的同種魚提供參考,后續(xù)研究可以利用該研究方法結(jié)合不同魚類的不同魚體范圍探索適應(yīng)的感官圈范圍.水力學(xué)數(shù)值模擬計算中,紊流條件越復(fù)雜則網(wǎng)格劃分越精密,而水流平順處網(wǎng)格則相對稀疏,在對魚體周身的水力學(xué)環(huán)境進行紊流數(shù)值模擬時,僅根據(jù)水力學(xué)條件劃分的網(wǎng)格可能不足以滿足魚類捕捉有效水力學(xué)信息的要求,需參考魚類感官范圍對水力學(xué)網(wǎng)格進行局部加密使得魚體能夠捕獲有效的刺激因素完成搜索上溯.此外研究發(fā)現(xiàn),齊口裂腹魚在上溯運動過程中對流速、感官搜索角度的范圍有一定的趨向性,它更偏向于在0.20~0.40m/s的流速區(qū)間內(nèi)以和水平方向呈-49°~49°的扇形感官搜索角度范圍進行上溯洄游.

    魚類上溯行為涉及到魚類的生態(tài)行為學(xué)與水力學(xué)的耦合,是一個比較前沿的學(xué)科交叉問題,魚類在上溯過程中所表現(xiàn)出來的行為也是極其復(fù)雜的,對其上溯行為的研究也是非常具有難度的[25],在跟蹤單尾魚樣本中首先遇到的困難是問題描述中的組合視角.魚道的水力信息被認為是一種刺激物,屬于歐拉流場,而魚對這些刺激物的響應(yīng)屬于拉格朗日參照系.運用Eulerian-Lagrangian-Agent方法(ELAM)對模擬魚類復(fù)雜的生態(tài)過程具有增值效益[26-27].

    生物粒子模型將魚視作粒子來模擬其上溯行為,其原理為首先對魚類所在區(qū)域的環(huán)境流場進行網(wǎng)格劃分,其中每個網(wǎng)格節(jié)點包含一定的水力學(xué)信息,包括上述的流速、流速梯度和湍動能等,然后魚類感知并評估側(cè)線感知范圍內(nèi)的偏好水力學(xué)條件[28],據(jù)此調(diào)整自身游泳策略進行搜索上溯.模型中魚類偏好水力學(xué)條件可根據(jù)AHP層次分析法通過探索魚類游泳行為對水力學(xué)條件的響應(yīng)關(guān)系進行確定[13],此外,利用生物粒子模型預(yù)測的魚類游泳軌跡結(jié)果還受數(shù)值模擬的網(wǎng)格精度、感官圈搜索范圍、水質(zhì)參數(shù)、魚體理化狀態(tài)及各種環(huán)境因子的影響,各種因子共同影響著魚類自主游泳行為.如魚類長時間處于較大紊動和切應(yīng)力環(huán)境下會喪失感官方向,最終導(dǎo)致其游泳能力降低[29];袁佺等[30]指出長薄鰍游泳行為中的回避率隨著水體環(huán)境總?cè)芙鈿怏w過飽和(TDG)的提高而提高;Morrice等[31]得出了河流流量與潮汐對大鱗大馬哈魚幼魚遷移路徑、河口停留時間和生長率具有顯著性的影響;LaBone等[32]采用三維耦合水動力水質(zhì)模型生成的DO場運動算法,提出了基于環(huán)境水體溶解氧水平的3種魚類垂直感觀距離.因此,在后續(xù)模型構(gòu)建的過程中,需綜合考慮魚類感官圈內(nèi)的各種水力學(xué)、水質(zhì)及魚類理化狀態(tài)參數(shù)的影響以優(yōu)化生物粒子模型,真實客觀模擬魚類在不同水流環(huán)境中的運動行為,為建設(shè)高效率過魚設(shè)施及評估過魚效率提供參考.

    4 結(jié)論

    4.1 采用小于或等于0.30m尺寸的網(wǎng)格解譯魚體行為時,魚類上溯運動平均游泳速度、轉(zhuǎn)向次數(shù)以及轉(zhuǎn)角和轉(zhuǎn)角差均與實際情況無顯著性差異.

    4.2 魚類偏向于在0.20~0.40m/s的流速區(qū)間內(nèi)以和水平方向呈-49°~49°的運動搜索范圍進行上溯運動.

    [1] Silva A T,Lucas M C,Castro-Santos T,et al. The future of fish passage science,engineering,and practice [J]. Fish and Fisheries,2018,19(2):340-362.

    [2] 李敏訥,朱海峰,金志軍,等.加速流對鰱、鳙幼魚下行過程中游泳行為的影響 [J]. 水生生物學(xué)報,2018,42(3):571-577.

    Li M N,Zhu H F,Jin Z J,et al. The effect of accelerating flow on swimming behavior of juvenile aristichthys nobilis and hypophthalmichthys molitrix during downstream migration [J]. Acta Hydrobiologica Sinica,2018,42(3):571-577.

    [3] 孫 林,易文敏,張興磊,等.變流場情景下基于魚類游泳特性的上溯通道分析[J]. 水力發(fā)電學(xué)報,2021,40(10):1-9.

    Sun L,Yi W M,Zhang X L,et al. Analysis of upstream migrating channels based on characteristics of fish swimming in variable flow field [J]. Journal of Hydroelectric Engineering,2021,40(10):1-9.

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

    Chen K Q,Chang Z N,Cao X H,et al. Recent progress in studies of overland flow resistance [J]. Journal of Hydraulic Engineering,2012,43(2):182-188.

    [5] 林俊強,彭期冬,黃真理.河流魚類魚卵運動的水力學(xué)研究進展[J]. 水利學(xué)報,2015,46(7):869-876.

    Lin J Q. Peng Q Y,Huang Z L. Review on hydraulics research of fish eggs’ movement in rivers [J]. Journal of Hydraulic Engineering,2015,46(7):869-876.

    [6] Zielinski D P,Voller V R,Sorensen P W. A physiologically inspired agent-based approach to model upstream passage of invasive fish at a lock-and-dam [J]. Ecological Modelling,2018,382:18-32.

    [7] Montgomery J,Coombs S,Halstead M. Biology of the mechanosensory lateral line in fishes [J]. Reviews in Fish Biology and Fisheries,1995,5(4):399-416.

    [8] Goodwin R A,Nestler J M,Anderson J J,et al. Forecasting 3-D fish movement behavior using a Eulerian–Lagrangian–agent method (ELAM) [J]. Ecological Modelling,2006,192(1/2):197-223.

    [9] Coombs S. Signal detection theory,lateral-line excitation patterns and prey capture behaviour of mottled sculpin [J]. Animal Behaviour,1999,58(2):421-430.

    [10] Gao Z,Andersson H I,Dai H,et al. A new Eulerian–Lagrangian agent method to model fish paths in a vertical slot fishway [J]. Ecological Engineering,2016,88:217-225.

    [11] Tan J,Tao L,Gao Z,et al. Modeling fish movement trajectories in relation to hydraulic response relationships in an experimental fishway [J]. Water,2018,10(11):1511.

    [12] Zielinski D P,Voller V R,Sorensen P W. A physiologically inspired agent-based approach to model upstream passage of invasive fish at a lock-and-dam [J]. Ecological Modelling,2018,382:18-32.

    [13] Zhu L,Li J,Deng Y,et al. Based on a biological particle model to predict the trace behavior of fish [J]. Water Supply,2021,21(8):4044- 4057.

    [14] 郭靜瑜,廖伯文,周其航,等.泥沙濃度對齊口裂腹魚生存的影響試驗研究[J]. 泥沙研究,2020,45(4):7-12.

    Guo J Y,Liao B W,Zhou Q H,et al. Experimental study on influence of sediment concentration to survival of schizothorax prenanti [J]. Journal of Sediment Research,2020,45(4):7-12.

    [15] 劉 軍.長江上游特有魚類受威脅及優(yōu)先保護順序的定量分析 [J]. 中國環(huán)境科學(xué),2004,24(4):395-399.

    Liu J. A quantitative analysis on threat and priority of conservation order of the endemic fishes in upper reaches of the Yangtze River [J]. China Environmental Science,2004,24(4):395-399.

    [16] Rodriguez á,Bermúdez M,Rabu?al J R,et al. Optical fish trajectory measurement in fishways through computer vision and artificial neural networks [J]. Journal of Computing in Civil Engineering,2011,25(4): 291-301.

    [17] 田 康,張 堯,李金龍,等.基于OpenFOAM幾何流體體積方法的波浪數(shù)值模擬[J]. 上海交通大學(xué)學(xué)報,2021,55(1):1.

    Tian K,Zhang Y,Li J L,et al. Numerical waves simulation using geometrical VOF method based on OpenFOAM [J]. Journal of Shanghai Jiaotong University,2021,55(1):1.

    [18] 廖琦宇.基于全局軌跡分析的視頻分割與合成算法研究[D]. 合肥:中國科學(xué)技術(shù)大學(xué),2015.

    Liao Q Y. Research on video segmentation and merging based on global trajectory analyze [D]. Hefei: University of Science and Technology of China,2015.

    [19] Bian L. The representation of the environment in the context of individual-based modeling [J]. Ecological Modelling,2003,159(2/3): 279-296.

    [20] Liao L,Chen M,An R,et al. Identifying three-dimensional swimming corridors for fish to match their swimming characteristics under different hydropower plant operations: Optimization of entrance location for fish-passing facilities [J]. Science of the Total Environment,2022,822:153599.

    [21] 蔣永強,楊忠勇,石小濤,等.基于多重水力學(xué)因子的豎縫式魚道中魚類上溯軌跡模擬 [J]. 生態(tài)學(xué)雜志,2018,37(4):1282-1290.

    Jiang Y Q,Yang Z Y,Shi X T,et al. The simulation of fish migratory trajectory in a vertical slot fishway based on multihydraulic indices [J]. Chinese Journal of Ecology,2018,37(4):1282-1290.

    [22] Arenas A,Politano M,Weber L,et al. Analysis of movements and behavior of smolts swimming in hydropower reservoirs [J]. Ecological Modelling,2015,312:292-307.

    [23] Shi X,Ke S,Tu Z,et al. Swimming capability of target fish from eight hydropower stations in China relative to fishway design [J]. Canadian Journal of Fisheries and Aquatic Sciences,2022,79(1):124-132.

    [24] 蔡 露,侯軼群,金 瑤,等.魚游泳能力對體長的響應(yīng)及其在魚道設(shè)計中的應(yīng)用 [J]. 農(nóng)業(yè)工程學(xué)報,2021,37(5):209-215.

    [25] Steel E A,Guttorp P,Anderson J J,et al. Modeling juvenile salmon migration using a simple Markov chain [J]. Journal of Agricultural,Biological,and Environmental Statistics,2001,6(1):80-88.

    [26] Kuli? T,Lon?ar G,Kova?evi? M,et al. Application of agent-based modelling for selecting configuration of vertical slot fishway [J]. Gra?evinar,2021,73(3):235-247.

    [27] Benson T,De Bie J,Gaskell J,et al. Agent-based modelling of juvenile eel migration via selective tidal stream transport [J]. Ecological Modelling,2021,443:109448.

    [28] Kanter M J,Coombs S. Rheotaxis and prey detection in uniform currents by Lake Michigan mottled sculpin (Cottus bairdi) [J]. Journal of Experimental Biology,2003,206(1):59-70.

    [29] 齊 亮,楊 宇,王 悅,等.魚類對水動力環(huán)境變化的行為響應(yīng)特征 [J]. 河海大學(xué)學(xué)報(自然科學(xué)版),2012,40(4):438-445.

    Qi L,Yang Y,Wang Y,et al. Fish behavior characteristics in response to change of hydrodynamic environment [J]. Journal of Hohai University,2012,40(4):438-445.

    [30] 袁 佺,袁 嫄,王遠銘,等.長薄鰍對過飽和總?cè)芙鈿怏w的回避特征研究 [J]. 水生態(tài)學(xué)雜志,2017,38(4):77-81.

    Yuan Q,Yuan Y,Wang Y M,et al. Avoidance characteristics of leptobotia elongate to total dissolved gas supersaturation [J]. Journal of Hydroecology,2017,38(4):77-81.

    [31] Morrice K J,Baptista A M,Burke B J. Environmental and behavioral controls on juvenile Chinook salmon migration pathways in the Columbia River estuary [J]. Ecological Modelling,2020,427:109003.

    [32] LaBone E,Justic D,Rose K,et al. Modeling fish movement in 3-D in the Gulf of Mexico hypoxic zone [J]. Estuaries and Coasts,2019,42(6):1662-1685.

    Interpretation of sensory query distance during upstream migration of.

    DENG Yue-xi1,TANG Xi-liang2,YAN Zhong-luan2,WANG Xiao-ming2,LU Jing-ying3,AN Rui-dong1,LI Jia1*

    (1.State Key Laboratory of Hydraulics and Mountain River Engineering,College of Water Resources and Hydropower,Sichuan University,Chengdu 610065,China;2.China Three Gorges Construction Engineering Corporation,Chengdu 610041,China;3.China Three Gorges Corporation,Wuhan 430010,China).,2022,42(8):3882~3889

    In this paper,an endemic fish located in Southwest China,(32.32±2.73)cm,was selected as the research object to test the volitional swimming behavior during the upstream migration in an open-channel flume with alternative flow regimes,and to interpret the extent to which swimming speed,turning times,angle and angle difference under different mesh sizes can match with the real swimming trajectory. The results showed that when the mesh size was less than or equal to 0.30m,there was no significant difference between the average swimming speed,the number of turns,the angle and the difference of turns with real swimming data. In addition,fish are more inclined to swim in the velocity ranged from 0.20 to 0.40m/s with the search angle ranged from -49o to 49o in horizontal direction.

    mesh size;eco-hydraulics;fish behavior;sensory query distance

    X171

    A

    1000-6923(2022)08-3882-08

    2022-01-10

    國家自然科學(xué)基金資助項目(52122904)

    * 責(zé)任作者,教授,lijia@scu.edu.cn

    鄧月曦(1998-),女,四川綿陽人,四川大學(xué)碩士研究生,主要從事環(huán)境與生態(tài)水力學(xué)研究.

    猜你喜歡
    腹魚轉(zhuǎn)角魚類
    人工養(yǎng)殖2 種裂腹魚生物學(xué)指數(shù)和肌肉營養(yǎng)價值的比較分析
    齊口裂腹魚集群行為對流態(tài)的響應(yīng)
    基于MFCC和ResNet的魚類行為識別
    南疆地區(qū)四種土著裂腹魚類幼魚對氯化鈉和碳酸氫鈉的急性耐受能力
    玩轉(zhuǎn)角的平分線
    三次“轉(zhuǎn)角”遇到愛
    解放軍健康(2017年5期)2017-08-01 06:27:42
    奇妙的古代動物 泥盆紀的魚類
    永春堂贏在轉(zhuǎn)角
    魚類是怎樣保護自己的
    下一個轉(zhuǎn)角:邁出去 開啟“智”造時代
    伦精品一区二区三区| 国产精品综合久久久久久久免费| 欧美丝袜亚洲另类| 国产亚洲91精品色在线| 亚洲乱码一区二区免费版| 在线天堂最新版资源| 亚洲在线自拍视频| 国产 一区 欧美 日韩| 在现免费观看毛片| 丰满乱子伦码专区| 亚洲怡红院男人天堂| 色网站视频免费| 国产 一区 欧美 日韩| 亚洲欧美精品综合久久99| 淫秽高清视频在线观看| 少妇人妻一区二区三区视频| 又粗又爽又猛毛片免费看| 国内少妇人妻偷人精品xxx网站| 天堂中文最新版在线下载 | 日韩国内少妇激情av| 久久久欧美国产精品| 亚洲天堂国产精品一区在线| 国产精品久久电影中文字幕| 国产精品一二三区在线看| 成人av在线播放网站| 一区二区三区免费毛片| 美女内射精品一级片tv| 激情 狠狠 欧美| 日本猛色少妇xxxxx猛交久久| 欧美bdsm另类| 日韩人妻高清精品专区| 国产成人a区在线观看| 国产久久久一区二区三区| a级毛色黄片| 国产淫片久久久久久久久| 一区二区三区免费毛片| av国产久精品久网站免费入址| av在线播放精品| 99热网站在线观看| 在线观看一区二区三区| 看片在线看免费视频| 久久国内精品自在自线图片| 美女xxoo啪啪120秒动态图| 国产成人免费观看mmmm| 人人妻人人澡欧美一区二区| av天堂中文字幕网| av在线亚洲专区| 一级爰片在线观看| 欧美极品一区二区三区四区| 简卡轻食公司| 欧美成人a在线观看| 亚洲精品国产成人久久av| 国产精品一及| 国产 一区精品| 久久精品久久久久久噜噜老黄 | 国产成年人精品一区二区| 日本五十路高清| 国产精品久久久久久久电影| 成年版毛片免费区| av天堂中文字幕网| 青春草视频在线免费观看| 成人av在线播放网站| av国产久精品久网站免费入址| 男女视频在线观看网站免费| 日本免费a在线| 爱豆传媒免费全集在线观看| 久久久久精品久久久久真实原创| 天天躁夜夜躁狠狠久久av| 老司机影院毛片| 熟女电影av网| 伦理电影大哥的女人| h日本视频在线播放| h日本视频在线播放| 国产毛片a区久久久久| 国产午夜精品一二区理论片| 国产探花在线观看一区二区| 亚洲美女搞黄在线观看| 亚洲欧美一区二区三区国产| 亚洲欧美一区二区三区国产| 好男人在线观看高清免费视频| 美女cb高潮喷水在线观看| 久久精品综合一区二区三区| 少妇丰满av| 又粗又爽又猛毛片免费看| 麻豆一二三区av精品| 亚洲精品乱码久久久v下载方式| 午夜福利成人在线免费观看| 久久精品国产亚洲网站| 亚洲av免费高清在线观看| 插阴视频在线观看视频| 国产白丝娇喘喷水9色精品| 欧美潮喷喷水| 日韩一区二区视频免费看| 亚洲av成人精品一区久久| 中文字幕制服av| 又爽又黄无遮挡网站| 人妻少妇偷人精品九色| 一本一本综合久久| 一级毛片我不卡| 一级毛片我不卡| 亚洲五月天丁香| 日本午夜av视频| 国产熟女欧美一区二区| 精品少妇黑人巨大在线播放 | 又粗又爽又猛毛片免费看| 日本五十路高清| 婷婷六月久久综合丁香| 久久人人爽人人爽人人片va| 亚洲av一区综合| 欧美xxxx黑人xx丫x性爽| 亚洲av熟女| 亚洲高清免费不卡视频| 亚洲精品国产av成人精品| 国产亚洲av嫩草精品影院| 性色avwww在线观看| 欧美成人a在线观看| 国语对白做爰xxxⅹ性视频网站| 春色校园在线视频观看| 亚洲综合精品二区| 性色avwww在线观看| or卡值多少钱| 七月丁香在线播放| 国产亚洲91精品色在线| 久久久久久久午夜电影| 亚洲国产欧美人成| 欧美性猛交黑人性爽| 激情 狠狠 欧美| 国产黄片视频在线免费观看| 91精品伊人久久大香线蕉| 国产在视频线精品| 男人的好看免费观看在线视频| 国产精华一区二区三区| 日韩亚洲欧美综合| 乱系列少妇在线播放| 青春草国产在线视频| 秋霞在线观看毛片| 啦啦啦啦在线视频资源| 精品久久久噜噜| 亚洲性久久影院| 国语自产精品视频在线第100页| 一个人看的www免费观看视频| 亚洲无线观看免费| 黄色一级大片看看| 亚洲中文字幕日韩| 日韩av在线免费看完整版不卡| 国产三级在线视频| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品久久久久久精品电影| 国产精品久久久久久精品电影| 欧美一区二区精品小视频在线| 亚洲真实伦在线观看| 色噜噜av男人的天堂激情| 国产成人精品婷婷| 亚洲国产欧洲综合997久久,| 欧美一区二区国产精品久久精品| 欧美一区二区国产精品久久精品| 男女国产视频网站| 最后的刺客免费高清国语| 春色校园在线视频观看| 人妻少妇偷人精品九色| 精品无人区乱码1区二区| 又爽又黄a免费视频| 成人一区二区视频在线观看| 国产大屁股一区二区在线视频| 国产高清三级在线| 日韩成人伦理影院| 日本熟妇午夜| 人妻制服诱惑在线中文字幕| 亚洲欧美日韩东京热| www.av在线官网国产| 色综合色国产| 小说图片视频综合网站| 日本-黄色视频高清免费观看| 亚洲第一区二区三区不卡| 永久网站在线| 成人毛片60女人毛片免费| 永久免费av网站大全| 男的添女的下面高潮视频| 午夜福利高清视频| 亚洲av中文av极速乱| av在线观看视频网站免费| 亚洲成人精品中文字幕电影| 久久国产乱子免费精品| 免费观看a级毛片全部| 日韩精品有码人妻一区| 97超视频在线观看视频| 床上黄色一级片| 欧美一区二区国产精品久久精品| 三级毛片av免费| 精品久久久久久成人av| 国产色婷婷99| 伦理电影大哥的女人| 我的老师免费观看完整版| 国产精品无大码| 日韩欧美三级三区| 久久久精品欧美日韩精品| 欧美性猛交╳xxx乱大交人| 女人十人毛片免费观看3o分钟| 亚洲av免费高清在线观看| 久久精品夜色国产| 非洲黑人性xxxx精品又粗又长| 久久久久久久国产电影| 床上黄色一级片| 国产精品一及| 乱人视频在线观看| 日韩av在线大香蕉| 能在线免费观看的黄片| 国产91av在线免费观看| 国产乱人偷精品视频| 久久久久久久亚洲中文字幕| 国产亚洲精品av在线| 一个人看的www免费观看视频| 人妻夜夜爽99麻豆av| 亚洲精品日韩av片在线观看| 久久亚洲精品不卡| 久久久久久国产a免费观看| 亚洲综合精品二区| 51国产日韩欧美| 晚上一个人看的免费电影| 日本一本二区三区精品| 九草在线视频观看| 蜜臀久久99精品久久宅男| 精品午夜福利在线看| 日韩av在线大香蕉| 成年免费大片在线观看| 大又大粗又爽又黄少妇毛片口| www.av在线官网国产| 日韩欧美 国产精品| 亚洲乱码一区二区免费版| 国产黄片视频在线免费观看| 精品欧美国产一区二区三| 高清毛片免费看| 国产在线男女| 国产探花极品一区二区| 国产精品综合久久久久久久免费| 久久精品久久久久久噜噜老黄 | 国产一区二区亚洲精品在线观看| 毛片一级片免费看久久久久| 国产老妇伦熟女老妇高清| 成年免费大片在线观看| 久久久久久久久久久丰满| 午夜精品一区二区三区免费看| a级一级毛片免费在线观看| 国产成人freesex在线| 欧美最新免费一区二区三区| 国产精品久久久久久久电影| 精品国产三级普通话版| 亚洲欧美日韩无卡精品| 国产精品人妻久久久影院| av视频在线观看入口| av女优亚洲男人天堂| 99热这里只有精品一区| 国产精品不卡视频一区二区| 国产乱人视频| 久久久久久久久久黄片| 亚洲五月天丁香| 久99久视频精品免费| 国产成人aa在线观看| 精品久久久久久久久久久久久| 亚洲最大成人中文| 一区二区三区高清视频在线| 国产私拍福利视频在线观看| 国产黄色小视频在线观看| 精品人妻熟女av久视频| 天堂影院成人在线观看| 91精品一卡2卡3卡4卡| 天堂中文最新版在线下载 | 国产精品av视频在线免费观看| 男女国产视频网站| 色噜噜av男人的天堂激情| 免费av不卡在线播放| 蜜臀久久99精品久久宅男| 国产精品人妻久久久影院| 性插视频无遮挡在线免费观看| 在线免费观看不下载黄p国产| 国产成人精品一,二区| 国产高清不卡午夜福利| 日本午夜av视频| 免费大片18禁| 99在线视频只有这里精品首页| 又黄又爽又刺激的免费视频.| 日韩成人伦理影院| ponron亚洲| 九色成人免费人妻av| 国产亚洲av嫩草精品影院| 如何舔出高潮| 久久久欧美国产精品| 男的添女的下面高潮视频| 国产精品,欧美在线| 最后的刺客免费高清国语| 老师上课跳d突然被开到最大视频| 久久久久国产网址| 在线观看一区二区三区| 嫩草影院新地址| 免费看日本二区| 在线播放无遮挡| 欧美三级亚洲精品| 日韩大片免费观看网站 | 六月丁香七月| 日韩一区二区三区影片| 黄色欧美视频在线观看| 黄色配什么色好看| 91久久精品电影网| 国产男人的电影天堂91| 好男人在线观看高清免费视频| 色视频www国产| 99热这里只有是精品50| 少妇熟女aⅴ在线视频| 国产三级中文精品| 又粗又硬又长又爽又黄的视频| 一区二区三区乱码不卡18| 女人被狂操c到高潮| 国产激情偷乱视频一区二区| 赤兔流量卡办理| 欧美区成人在线视频| 精品一区二区免费观看| 一二三四中文在线观看免费高清| 啦啦啦啦在线视频资源| 人人妻人人看人人澡| 中文欧美无线码| 听说在线观看完整版免费高清| 国产免费视频播放在线视频 | 免费在线观看成人毛片| 亚洲欧美成人精品一区二区| 国语自产精品视频在线第100页| 美女大奶头视频| 国产av不卡久久| 午夜老司机福利剧场| 成人鲁丝片一二三区免费| 国产亚洲av片在线观看秒播厂 | 99久久无色码亚洲精品果冻| 午夜激情欧美在线| 男女那种视频在线观看| 欧美不卡视频在线免费观看| 午夜福利在线在线| 亚洲国产欧洲综合997久久,| 成人无遮挡网站| 国产乱来视频区| 精品久久久久久久人妻蜜臀av| 插逼视频在线观看| 精品一区二区免费观看| 91久久精品国产一区二区成人| 精品欧美国产一区二区三| 中文字幕亚洲精品专区| 看黄色毛片网站| 日本免费一区二区三区高清不卡| 九九热线精品视视频播放| 久久人人爽人人片av| 久久久久久久久中文| 国产女主播在线喷水免费视频网站 | 国产精品一区二区三区四区久久| 久久久久久久久大av| 视频中文字幕在线观看| 韩国高清视频一区二区三区| 亚洲性久久影院| 中文欧美无线码| 少妇裸体淫交视频免费看高清| 色视频www国产| 色5月婷婷丁香| 汤姆久久久久久久影院中文字幕 | 国产成人freesex在线| 黄片无遮挡物在线观看| 乱人视频在线观看| 国产亚洲一区二区精品| 伦精品一区二区三区| 看非洲黑人一级黄片| 国产成人a∨麻豆精品| 亚洲国产精品国产精品| 亚洲在久久综合| 国产高潮美女av| 国产黄a三级三级三级人| 午夜福利视频1000在线观看| 国产在线一区二区三区精 | 亚洲电影在线观看av| 亚洲av.av天堂| 国产精品久久视频播放| 尤物成人国产欧美一区二区三区| 亚洲成色77777| 国产伦精品一区二区三区视频9| 丝袜美腿在线中文| 精品99又大又爽又粗少妇毛片| 黄色一级大片看看| 婷婷色综合大香蕉| 搡老妇女老女人老熟妇| 麻豆av噜噜一区二区三区| 中文资源天堂在线| 自拍偷自拍亚洲精品老妇| 免费看美女性在线毛片视频| 极品教师在线视频| 亚洲成人av在线免费| 欧美一区二区国产精品久久精品| 少妇人妻精品综合一区二区| 亚洲性久久影院| 搡女人真爽免费视频火全软件| 久久久久九九精品影院| 国产激情偷乱视频一区二区| av播播在线观看一区| 亚洲高清免费不卡视频| 国产日韩欧美在线精品| 日本黄大片高清| 舔av片在线| 亚洲国产精品久久男人天堂| 国产熟女欧美一区二区| 男女下面进入的视频免费午夜| 精品久久久久久久末码| 精品国内亚洲2022精品成人| 久久久午夜欧美精品| av在线播放精品| 亚洲天堂国产精品一区在线| 一区二区三区高清视频在线| 少妇裸体淫交视频免费看高清| 最近最新中文字幕大全电影3| 国产三级中文精品| 两个人视频免费观看高清| 精品久久久久久久久久久久久| 精品一区二区三区人妻视频| 舔av片在线| 婷婷色av中文字幕| 美女脱内裤让男人舔精品视频| 色哟哟·www| 18禁动态无遮挡网站| 日本-黄色视频高清免费观看| 精品久久国产蜜桃| 成人鲁丝片一二三区免费| 国内精品一区二区在线观看| 久久99热这里只有精品18| 麻豆成人午夜福利视频| 特级一级黄色大片| 精品久久久久久电影网 | 国产精品不卡视频一区二区| 国产不卡一卡二| 亚洲国产精品成人久久小说| 国产成人一区二区在线| 久久精品久久久久久噜噜老黄 | 特级一级黄色大片| 我的女老师完整版在线观看| 三级国产精品片| 啦啦啦观看免费观看视频高清| 全区人妻精品视频| 久久久精品94久久精品| 女的被弄到高潮叫床怎么办| 精品不卡国产一区二区三区| 少妇的逼好多水| 成人鲁丝片一二三区免费| 亚洲成人av在线免费| 秋霞伦理黄片| 网址你懂的国产日韩在线| 精品久久久久久电影网 | av.在线天堂| 一级二级三级毛片免费看| 国产综合懂色| 国产午夜精品久久久久久一区二区三区| 亚洲综合精品二区| 国产亚洲精品av在线| 日韩制服骚丝袜av| 久久人人爽人人片av| 在线观看66精品国产| 变态另类丝袜制服| 国产精品一二三区在线看| 大香蕉久久网| 少妇的逼好多水| 欧美高清性xxxxhd video| 久久久亚洲精品成人影院| 国产精品久久久久久av不卡| 麻豆一二三区av精品| 九九久久精品国产亚洲av麻豆| 国产精品国产三级国产专区5o | 午夜亚洲福利在线播放| 又黄又爽又刺激的免费视频.| 亚洲欧美精品专区久久| 国产老妇伦熟女老妇高清| 国产精品.久久久| 国国产精品蜜臀av免费| 成人二区视频| 亚洲自拍偷在线| 热99re8久久精品国产| 高清午夜精品一区二区三区| 我要搜黄色片| 十八禁国产超污无遮挡网站| 免费播放大片免费观看视频在线观看 | 淫秽高清视频在线观看| 内地一区二区视频在线| 天天躁夜夜躁狠狠久久av| 日韩高清综合在线| 97超碰精品成人国产| 国产成人免费观看mmmm| 99热这里只有是精品在线观看| 成人亚洲欧美一区二区av| 免费看美女性在线毛片视频| 精品国内亚洲2022精品成人| 国产精品久久久久久精品电影| 日本黄色片子视频| 人妻夜夜爽99麻豆av| 亚洲四区av| h日本视频在线播放| 偷拍熟女少妇极品色| 午夜视频国产福利| 国产伦精品一区二区三区视频9| 中文字幕精品亚洲无线码一区| 国产精品1区2区在线观看.| 噜噜噜噜噜久久久久久91| 尤物成人国产欧美一区二区三区| 日产精品乱码卡一卡2卡三| 日本欧美国产在线视频| 久久久久国产网址| 久久精品综合一区二区三区| 国产91av在线免费观看| 午夜日本视频在线| 国产亚洲av嫩草精品影院| 亚洲色图av天堂| 丰满乱子伦码专区| 国产成人精品一,二区| 久久99热这里只有精品18| 韩国高清视频一区二区三区| 身体一侧抽搐| 18禁在线无遮挡免费观看视频| 国产亚洲精品av在线| 高清av免费在线| 99久久人妻综合| 国产真实伦视频高清在线观看| 简卡轻食公司| 熟女电影av网| 性色avwww在线观看| 天天一区二区日本电影三级| 深爱激情五月婷婷| 精品久久久久久久人妻蜜臀av| 五月玫瑰六月丁香| 国产亚洲av片在线观看秒播厂 | 99久久人妻综合| 欧美日韩在线观看h| 国产成人aa在线观看| 日本-黄色视频高清免费观看| 日日摸夜夜添夜夜爱| 国产精品人妻久久久久久| 亚洲18禁久久av| 久久精品国产亚洲av天美| 最新中文字幕久久久久| 欧美丝袜亚洲另类| 高清午夜精品一区二区三区| 亚洲欧美日韩卡通动漫| 免费av观看视频| 亚洲国产欧美人成| 男女那种视频在线观看| 精品国产三级普通话版| 麻豆精品久久久久久蜜桃| 精品人妻偷拍中文字幕| 2021少妇久久久久久久久久久| 国产精品乱码一区二三区的特点| 18禁在线播放成人免费| 亚洲成av人片在线播放无| 免费观看人在逋| 熟妇人妻久久中文字幕3abv| 校园人妻丝袜中文字幕| 欧美激情国产日韩精品一区| 欧美一区二区亚洲| 看片在线看免费视频| 国产高清视频在线观看网站| 亚洲精品乱码久久久v下载方式| 日本午夜av视频| 99久国产av精品| 成人三级黄色视频| 99热这里只有是精品在线观看| 亚洲精品日韩在线中文字幕| 在线免费十八禁| 身体一侧抽搐| 禁无遮挡网站| 国产免费福利视频在线观看| 性色avwww在线观看| .国产精品久久| 色播亚洲综合网| 岛国毛片在线播放| 亚洲av一区综合| 国产午夜精品论理片| 美女xxoo啪啪120秒动态图| 国产午夜精品一二区理论片| 国产精品精品国产色婷婷| 国产精品人妻久久久久久| 久久久久久久午夜电影| 99热精品在线国产| 天堂√8在线中文| 99热这里只有是精品50| 亚洲av.av天堂| 高清毛片免费看| 日日撸夜夜添| 精品99又大又爽又粗少妇毛片| 免费观看a级毛片全部| 国产乱来视频区| 2022亚洲国产成人精品| 欧美xxxx性猛交bbbb| 日本-黄色视频高清免费观看| av卡一久久| 大香蕉久久网| 天美传媒精品一区二区| 亚洲性久久影院| 欧美一区二区亚洲| 内射极品少妇av片p| 久久综合国产亚洲精品| 伦精品一区二区三区| 日本av手机在线免费观看| 日日摸夜夜添夜夜添av毛片| or卡值多少钱| 国产成人午夜福利电影在线观看| 日本猛色少妇xxxxx猛交久久| 日本一本二区三区精品| 免费电影在线观看免费观看| 在线播放国产精品三级| 最近最新中文字幕免费大全7| 成人美女网站在线观看视频| 五月玫瑰六月丁香| 日本一本二区三区精品| 最近最新中文字幕免费大全7| 国产色爽女视频免费观看| 免费不卡的大黄色大毛片视频在线观看 | 三级国产精品欧美在线观看| 免费观看的影片在线观看|