棲息地適宜度模型在魚類生境中的應(yīng)用

董謀

摘要：基于DHIMIKE軟件中的MIKE21FM模塊建立二維水動力學(xué)模型，模擬在18000 m3/s和25000 m3/s流量下，研究江段魚類棲息地的流速、水深環(huán)境參數(shù)適宜度的時空分布特征，并以2019年5月勘測船實測的水文數(shù)據(jù)對模型的參數(shù)進(jìn)行驗證。結(jié)果表明：流速和水深是影響三峽水庫變動回水區(qū)魚類棲息生境的關(guān)鍵因子，且魚類棲息地適宜度受流速和水深共同的影響顯著。

關(guān)鍵詞：魚類;生境特性;適宜度指標(biāo);三峽水庫;變動回水區(qū)

中圖分類號：X174? ? ? ? 文獻(xiàn)標(biāo)識碼：A

河流棲息地，又稱為河流水生生物棲息地，是各種魚類群落賴以生存、繁衍的空間與環(huán)境，與水生生物多樣性密切相關(guān)，棲息地是河流生態(tài)系統(tǒng)的重要組成部分[1]。魚類棲息地包括產(chǎn)卵場、索餌場、越冬場以及洄游通道等，魚類的“三場一通道”是魚類完成其生活史過程所必需的水域范圍[2]，需要河流形貌、水環(huán)境、水文特性和水力學(xué)特征共同構(gòu)造復(fù)雜的物理環(huán)境，從而為魚類提供各種棲息、覓食、繁殖和掩蔽場所[3]?，F(xiàn)有研究認(rèn)為魚類的適宜度模型用于棲息地適宜度評估，指標(biāo)從幾個到幾十個不等[4]，但魚類對各指標(biāo)的適應(yīng)度隨保護(hù)目標(biāo)種類、生理階段和研究區(qū)域不同而存在差異，這造成了棲息地的評價模型十分復(fù)雜且難以通用[5]，如本研究河段的魚類聚集水深遠(yuǎn)大于其他流域。為減輕航道整治工程對魚類的影響，有必要獲取魚類棲居的位置。根據(jù)近兩年對的魚群分布調(diào)查結(jié)果，文章對研究江段魚類生境的流速和水深因子加以討論。

1? 平面二維水流數(shù)學(xué)模型的建立

1.1? 控制方程

天然河道由于地形復(fù)雜，蜿蜒曲折，在網(wǎng)格劃分過程中往往存在較大難度，為有效克服計算域邊界起伏變化幅度較大而引起的問題，本文選用計算網(wǎng)格與河道邊界貼合的方法進(jìn)行模擬計算（見公式1）：

其中，α=x2φ+y2φ;γ=x2μ+y2μ;J=xμyφ-xφyμ;P、Q為調(diào)節(jié)因子。若水域中的水體維持有勢運動，其流線簇與勢線簇必定會保持正交，由此可以導(dǎo)出以網(wǎng)格間距的變化作為調(diào)節(jié)因子的貼體正交曲線坐標(biāo)方程：

（1）水流連續(xù)方程（見公式2）

（2）μ方向動量方程（見公式3）

其中，μ、φ分別是正交曲線坐標(biāo)系中的兩個正交曲線坐標(biāo);u、ν分別是沿μ、φ方向上的流速;h表示水深;H表示水位;Cμ、Cφ是正交曲線坐標(biāo)系中的拉梅系數(shù);、;σμμ、σμφ、σφμ、σφφ表示紊動切應(yīng)力：

其中，Vt表示紊動黏性系數(shù)，即Vt=CμK2/ε。

式中，K及ε分別代表紊動動能及紊動耗散系數(shù)，由K以及ε的運輸方程確定。

1.2? 工程河段網(wǎng)格剖分

工程河段模擬范圍上起和尚灘水道（上游航道里程536 km），下止豐都土地盤（下游航道里程488 km），全長48 km。初始網(wǎng)格與擴散光順法使用的網(wǎng)格劃分基本上相同，網(wǎng)格總量為5000×51，其中縱向X沿河流流動方向網(wǎng)格劃分5000，橫向Y為河流寬度方向網(wǎng)格劃分51，地形數(shù)據(jù)為2019年實測數(shù)據(jù)。總體上在靠近河流流場的區(qū)域網(wǎng)格較為密集，在河流的網(wǎng)格壁面處網(wǎng)格較為稀疏;在寬度較小的順直型河流斷面處網(wǎng)格較密，在寬度較大的彎曲型或分汊型河流斷面處網(wǎng)格較為稀疏，見圖1。

2? 模型驗證

2.1? 驗證所采用的資料

模型驗證采用長江重慶航運工程勘察設(shè)計院2015年6月實測的涪陵至豐都水道地形圖，測圖長河段比例1∶2000，1954北京坐標(biāo)系，吳淞高程。采用2019年5月清溪水文站實測水面線進(jìn)行水位驗證，Q=14200 m3/s，H=149.23 m。根據(jù)水深實測值與二維模型計算的模擬值之間的水位相對誤差，調(diào)整研究江段模型中的曼寧系數(shù)（Manning Number），渦黏系數(shù)（Eddy viscosity）以及進(jìn)出口和邊界條件，使水位相對誤差最小。渦黏系數(shù)設(shè)定為全局分布，率定結(jié)果為0.28。

2.2? 水位驗證

利用實測水位和流速對河床糙率進(jìn)行率定，以清溪場附近水位模型驗證結(jié)果顯示，二維模型計算的模擬值和實測水面線走勢相同，水位偏差在±0.4 m以內(nèi)。

3? 模型結(jié)果

從其適宜度指標(biāo)的水深分布、流速分布，來進(jìn)一步分析適宜度指標(biāo)在研究河段中的分布規(guī)律。

4? 結(jié)論

魚類群落整體上以礁石后深沱江段和分汊區(qū)域為高適宜區(qū)。當(dāng)流量為18000 m3/s時，高適宜區(qū)域基本分布在深沱和平綏壩附近的區(qū)域，棲息的深度可作為生境營造的重要指標(biāo);當(dāng)流量為25000 m3/s時，平綏壩逐漸被水流淹沒，隨著流量的繼續(xù)增大，研究河段區(qū)域的水深越來越大，主航道的表層流速也越來越快，此時，魚類在水深適中處及河漫灘附近流速適宜度更高，可以有效減少水流紊亂的沖擊，但是隨著流量的繼續(xù)增加，兩岸河漫灘和平綏壩區(qū)域被淹沒深度越來越大，魚類群落的棲息地適宜范圍又逐漸減小。因此流速和水深是影響變動回水區(qū)魚類棲居的關(guān)鍵因子，且都有一定的適宜范圍值。

參考文獻(xiàn)

[1]熊飛，劉紅艷，段辛斌，等.長江上游江津江段吻鮈的種群參數(shù)[J].淡水漁業(yè)，2018，48（06）：27-32.

[2]王震，Arunjith Thundiparambil Sathrajith，謝松光，等.金沙江梯級大壩運行和三峽水庫運行水位增高對長江上游干流寡鱗飄魚仔魚豐度和分布的影響[J].水生生物學(xué)報，2019，43（03）：606-611.

[3]楊志，唐會元，朱迪，等.三峽水庫175 m試驗性蓄水期庫區(qū)及其上游江段魚類群落結(jié)構(gòu)時空分布格局[J].生態(tài)學(xué)報，2015，35（15）：5064-5075.

[4]徐薇，楊志，陳小娟，等.三峽水庫生態(tài)調(diào)度試驗對四大家魚產(chǎn)卵的影響分析[J].環(huán)境科學(xué)研究，2020，33（05）：1129-1139.

[5]付翔，付成，付世建.五種淡水魚類幼魚游泳能力的比較[J].生態(tài)學(xué)雜志，2020，39（05）：1629-1635.