沅江托口水電站減水河段魚類棲息生境變化模擬研究

摘要：通過二維水力模擬軟件，分析托口電站減水河段天然流量不同時期（豐、枯水期）與生態(tài)下泄流量下魚類適宜棲息面積變化，以沅河鎮(zhèn)段約3km長的河段為例，根據(jù)代表魚類對水深、流速適應(yīng)度的曲線，在不同工況下，計算出有效棲息地面積（WUA）。計算結(jié)果表明：托口水電站建設(shè)后生態(tài)流量狀況下，相較于建設(shè)前的河道枯水期和豐水期，平均水深減小0.09m和0.59m，河道流速下降70.42%和86.45%，水域面積分別減少6.64%和16.22%;WUA分別減少64.50%和61.72%。

關(guān)鍵詞：生境模擬;魚類棲息地;減水河段;有效棲息地面積

中圖分類號：X174 文獻標識碼：A 文章編號：2095-672X（2020）05-0-04

DOI：10.16647/j.cnki.cn15-1369/X.2020.05.115

Simulation of fish habitat changes in the reach of Tuokou hydropower station of Yuanjiang River

Zhai Lilei

（Datang Yuquhe Hydropower Development Co.，Ltd.，Chengdu Sichuan 610031，China）

Abstract：Based on the two-dimensional hydraulic simulation software， the change of fish suitable habitat area under different periods of natural flow （wet and dry periods） and ecological discharge in the water reducing reach of Tuokou hydropower station is analyzed.Taking the reach of about 3km in Yuanhe town as an example， according to the curve of fish's adaptability to water depth and velocity， the effective habitat area （WUA） is calculated under different conditions.The calculation results show that the average water depth is reduced by 0.09m and 0.59m，the flow velocity is reduced by 70.42% and 86.45%， the water area is reduced by 6.64% and 16.22%， and WUA is reduced by 64.50% and 61.72%， respectively，compared with the low water period and high water period before the construction of Tuokou hydropower station.

Key words：Habitat simulation;Fish habitat;Water reducing reach;Effective habitat area

水電的開發(fā)建設(shè)對河流生態(tài)系統(tǒng)產(chǎn)生的影響主要表現(xiàn)在河流阻隔、局部河段減脫水、庫區(qū)生態(tài)系統(tǒng)的演替及下游水文時空節(jié)律改變、低溫水下泄等方面。河段減水對河流的生態(tài)環(huán)境造成極大的破壞，水生生物棲息地安全性受到脅迫。為緩解減水對河段中魚類的影響，國內(nèi)外普遍采取生態(tài)流量下泄、生境修復(fù)等措施。

魚類棲息地指能夠滿足魚類不同生活階段的不同需要，不僅提供魚類的生存空間，同時還提供滿足魚類生存、生長、繁殖的全部環(huán)境因子。通常包括產(chǎn)卵場、索餌場、越冬場及連接不同生活史階段的洄游通道。

河流棲息地評估為評估下泄流量的有效性發(fā)揮重要作用，棲息地模擬法是常用的河流生物棲息地評估方法，河道內(nèi)流量增量法（Instream Flow Incremental Methodology， 簡稱IFIM）中River2D模型運用較廣泛。River 2D 模型通過水力模型模擬，可以獲得河流局部流速、流向及水深， 根據(jù)魚類適應(yīng)性指數(shù)算出魚類有效棲息地面積。通過將水力學模型與生物信息模型相結(jié)合，定量分析魚類適宜棲息地與流量之間的關(guān)系，對流量變化對魚類棲息地的影響進行評價（黃月群等，2017）。崔玲等（2007）利用River 2D漓江桂林市區(qū)江段進行流速和水深模擬。譚燕平等（2011）對雅礱江錦屏大河灣不同流量下裂腹魚的適宜棲息地面積的變化，得出有效棲息地面積最大的流量。孫嘉寧等（2013）采用River 2D 分析白鶴灘水庫回水支流黑水河在不同蓄水水位下魚類適宜棲息地面積的變化情況。孫瑩等（2015）模擬壩下齊口裂腹魚產(chǎn)卵場，為產(chǎn)卵場地形重塑效果進行評估。李平等（2016）利用River 2D 軟件模擬不同的河漫灘糙率對應(yīng)的魚類有效棲息地面積（WUA），分析了糙率的改變對河流魚類棲息地的影響。

托口水電站壩址位于湖南省沅水干流上游河段，是沅水規(guī)劃梯級的第5級，壩址處多年平均流量為540m3/s。主副壩的修建，形成約9km的減水河段，根據(jù)環(huán)境保護要求，為保證河道內(nèi)的最小生態(tài)需水量（Tennant，1976），采取主壩壩后設(shè)置2組1.5萬KW的生態(tài)機組下泄生態(tài)環(huán)境流量54 m?/s。托口水電站河段中喜靜水緩流生境的魚類增多，喜急流和洄游、半洄游性魚類減少。通過生境模擬分析天然流量不同時期（豐、枯水期）與生態(tài)下泄流量下魚類適宜棲息面積，建立沅河的二維物理棲息地模型，旨在保護沅河減水河段的魚類棲息地，恢復(fù)適宜魚類棲息的生境條件，對保護魚類資源、維護河流生物多樣性具有重要意義，并為河道生境修復(fù)規(guī)劃方案的制定提供基礎(chǔ)數(shù)據(jù)。

1研究河段魚類資源及棲息地現(xiàn)狀

研究河段位于沅江上游，沅河鎮(zhèn)河段河勢蜿蜒，河岸邊灘分布廣泛，魚類資源豐富，隨著上下游梯級電站的興建、過度捕撈、水域污染、采砂活動頻繁，造成魚類資源下降。湘華鯪是湖南特有的大型野生經(jīng)濟魚類，目前僅在沅江和資江發(fā)現(xiàn)分布（卞偉等，2011）。沅江鎮(zhèn)河段是湘華鯪的自然分布江段，托口水電站對其分布產(chǎn)生嚴重的阻隔，根據(jù)環(huán)保要求，拖口水電站每年對湘華鯪進行增殖放流，沅河鎮(zhèn)減水河段是湘華鯪潛在的分布水域，湘華鯪是研究河段中的保護物種，且具有較高的保護價值，本文將湘華鯪作為目標物種。

湘華鯪以攝食硅藻為主，4月下旬性成熟，5月為產(chǎn)卵盛期，具有繁殖期上溯洄游習性，卵微粘性（王賓賢等，1977&1982），喜棲息于有砂、卵石的淺灘或狹流凸石的水域，石上生長著許多的固著藻類，如苔蘚、絲狀藻等。

2 評估方法

2.1 評估方法與模型

IFIM 法是棲息地模擬方法具有代表性一種方法。IFIM 法設(shè)定流量變化對物種數(shù)量和分布造成影響的主要因素為水深、流速、基質(zhì)和覆蓋物，通過調(diào)查分析目標物種對水深、流速、底質(zhì)等因子的適應(yīng)性，繪制水深、流速、底質(zhì)等環(huán)境參數(shù)與適應(yīng)性指數(shù)（介于0～1之間） 之間的適宜性曲線，計算指示物種的有效棲息地面積（ WUA） 。公式如下：

式中CSF（ Vi，Di，Ci） 為每個影響因子的組合適宜性值，其中，i 為劃分的單元數(shù)，V 為流速適宜指數(shù)，D 為水深適宜指數(shù)，C 為河道適宜指數(shù)（ 基質(zhì)和覆蓋物） ;Ai為每個單元的水平面積。本文主要考慮水深和流速的影響。對于棲息地組合適宜性值的確定，采取影響因子適宜指數(shù)相乘的方式，從而體現(xiàn)影響因子的綜合作用。

River2D 模型全稱為河流水動力學和魚類生境的二維平均深度模型是2002年由加拿大淡水研究所、亞伯特大學土木與環(huán)境系、美國地質(zhì)調(diào)查局中部生態(tài)科學中心、亞伯特漁業(yè)部共同研究開發(fā)與測試的棲息地模擬軟件。它將水力模擬和生態(tài)模擬緊密結(jié)合，主要用于微生境（水深、流速、基質(zhì)和河面覆蓋等） 的模擬評價。能反映河流局部的流速和水深分布以及確定流量與魚類潛在適宜棲息地的關(guān)系。

2.2 水深及流速適宜性指數(shù)確定

棲息地適宜性指數(shù) （ Habitat Suitability Indices，HIS）（ Bovee，1986）通常以0～1 之間的數(shù)值來量化物種對某棲息指標的喜好程度，0表示該物種完全不適宜，1 表示該物種完全適宜，值越大表示適宜性越好。

目前湘華鯪相對應(yīng)的水力適宜性研究資料還比較缺乏，本文根據(jù)卞偉等（2011）測定湘華鯪產(chǎn)卵場江段的水深和流速，結(jié)合湘華鯪生態(tài)習性、繁殖習性，確定湘華鯪的最適宜棲息流速為0.86～1.25 m/s，最適宜棲息水深為2～6.8 m。由此確定研究河段湘華鯪在繁殖期對流速、水深的適宜性曲線，如圖1 和圖2 所示。

2.3 模擬工況

沅河鎮(zhèn)河段處于托口水電站減水河段，托口水電站壩址多年平均流量為540 m3/s，枯水期（10-3月）為273.8 m3/s，豐水期（4-7月）為969.3 m3/s，托口水電站運行期主壩下泄生態(tài)流量為54 m3/s。本文以沅河鎮(zhèn)上下游3km河段為例，根據(jù)三種不同河道流量情況，運用River 2D模型對河段生境進行模擬。河段流量工況如表1所示。

通過對3km的模擬河段的河床地形數(shù)據(jù)進行實測了483個點位，并結(jié)合河床遙感數(shù)據(jù)，對部分河段進行修正，在R2D-Bed.exe中根據(jù)三角化不規(guī)則網(wǎng)格法（TIN）對地形特征以及各個節(jié)點參數(shù)進行空間內(nèi)插，形成河流河床文件。R2D-Mesh.exe中由地形文件生成計算離散化網(wǎng)格，網(wǎng)格密度根據(jù)地形數(shù)據(jù)和模型計算結(jié)果調(diào)整，網(wǎng)格間距0.5～10.0m。邊界條件上游采用流量邊界，下游采用水位邊界，河岸邊界采用動邊界方式。

3 生境模擬

通過模擬分析，沅河鎮(zhèn)河段在枯水期、豐水期和生態(tài)流量下泄是的河段的平均水深分別為1.85m，2.35m和1.76m，水深變化不大，變幅在0.6m以內(nèi);平均流速在三種流量工況下分別為0.71m/s，1.55m/s和0.21m/s，河道內(nèi)在生態(tài)流量下泄時，水體流速顯著降低，較豐水期和枯水期分別下降70.42%和86.45%。河道水深和流速分布見圖3。

河道在托口水電站建設(shè)前，沅河鎮(zhèn)河段在枯水期年均273.8 m3/s流量情況下的水域總面積為72.29 hm2，其中適宜性較差的水域面積為19.57 hm2，有效棲息地面積為40.00 hm2，適宜性較好水域面積為36.05 hm2。豐水期年均969.3m3/s流量情況下的水域總面積為80.55 hm2，其中適宜性較差的水域面積為33.52 hm2，有效棲息地面積為37.09 hm2，適宜性較好水域面積為31.00hm2。在豐水期湘華鯪有效棲息地面積較枯水期略有下降。托口水電站建成后，減水河道下泄54 m3/s的生態(tài)流量，水域面積總面積為67.49 hm2，分別較枯水期和豐水期減少6.64%和16.22%;適宜性較差的水域面積為31.73 hm2，較枯水期增加62.14%，較豐水期減少5.35%;適宜性較好水域面積為0.75 hm2，分別較枯水期和豐水期減少97.92%和97.58%;有效棲息地面積為14.20 hm2，分別較枯水期和豐水期減少64.50%和61.72%。

4 結(jié)果與討論

水深是魚類棲息地水力學特征重要的指標之一，水深會影響水體的光線穿透性和水溫等水文條件，同時也會影響魚類、水生植物、底棲生物等生物的種類、密度、分布特征等，從而導(dǎo)致不同水深層魚類食物的種類和數(shù)量存在差異，魚類物種組成和群落數(shù)量也隨之改變（馮廣朋等，2006）。不同生活史階段的魚類對水深也是有選擇特性的（杜浩等，2010），此外，要保持水體連通性、魚類棲息活動，維持一定的水深是必要的（王銳和李嘉，2010）。水流是河流生境的主要決定因素，同時也是生物組成的決定因素，強烈影響了生物結(jié)構(gòu)和行為適應(yīng)性，魚類的習性與流速更是密切相關(guān)，如草魚洄游的適宜流速范圍是0.4～1.0m/s，但超過1.1m/s后，洄游成功率將明顯降低（楊慶，2019）。

模擬分析表明，托口水電站壩下減水河段在流量明顯下降時，水域面積沒有大幅度的減少，較豐水期水位下降約0.5m，水域面積減少13.06hm2，分別減少了25.1%和16.2%。這與大渡河黃金坪水電站（王琴華，2015）、永樂江一級電站（何懷光，2012）等造成減水河段水域面積明顯下降結(jié)果不一致，這些電站減水河段的年平均流量、年平均水深、年平均流速、水域面積等較建壩前明顯下降，河道平均寬度也有所減少，這些因素改變了減水河段水力生態(tài)特征，對河道中魚類及其它生物的生存和繁衍產(chǎn)生了重大影響。沅河鎮(zhèn)減水河段由于河道下游有沙洲，河道深淺相間，當?shù)鼐用褚矠楸ＷC生活需水，對河道分叉處修筑泥土壩，阻隔部分支岔河道水流，保持河段的水面維持在一個較高的水位。但同時，河道內(nèi)流量的大幅度減少，流速顯著下降，使原為山溪性急流的河段變?yōu)殪o緩流水體，河道內(nèi)天然平均流速從0.71～1.55m/s，下降到減水后的平均0.21m/s，減少幅度為70.42%～86.45%。這使得適應(yīng)流水生境的湘華鯪等魚類，在該河段適宜性較好的水域面積幾乎喪失，面積減少在97.5%以上，只在河段上游人工圍筑收窄的河段，保持有流水河段，具有少量適宜性較好的水域;而適宜性一般的水域相對于枯水期和豐水期分別增長了110.00%和118.40%。說明湘華鯪對河段的棲息地適宜度整體顯著下降，經(jīng)統(tǒng)計湘華鯪有效棲息面積從枯水期和豐水期的40.00 hm2和37.09 hm2，下降為14.20 hm2，減少幅度超過60%。

減水河段水深下降、流速減緩，使適應(yīng)流水生境的魚類的適宜生境喪失。在模擬河段中，生態(tài)流量狀況下，僅在人為改造的一小塊河段存在適宜性較好的生境，分布較零散。在黑水河，為緩解減水河段流速過緩、水深偏小、魚類棲息地適宜性較差等問題，孫瑩等（2015）設(shè)計對河道地形重塑，增加目標魚類適宜棲息面積，達到對河道生境進行修復(fù)的目的。改善河道生境條件，國內(nèi)外常用的方法有增加河流彎道、設(shè)置水流擋板、在河床堆置大石塊成堰、布置人工階梯—深潭等（Bockelmann et al， 2004; Roni et al， 2006; Diego &Javier， 2007;余國安等，2008;蒲玉芳等，2013）。托口電站減水河段通過對部分河段進行重塑，改善河道的水流條件，恢復(fù)河流的魚類棲息地適宜性。

5 結(jié)語

本文基于IFIM 原理的棲息地模擬方法，使用River2D 二維水動力學軟件，以托口水電站壩下減水河段的沅河鎮(zhèn)河段為例，進行電站建設(shè)前后棲息地質(zhì)量和空間分布情況模擬，未考慮河床底質(zhì)及餌料分布情況，分析減水河段中湘華鯪的適宜生境顯著減少，適宜性較好的河段基本喪失。本文選取的流速、水深適宜性指數(shù)較粗略，未針對湘華鯪繁殖期生境的適宜性進行分析，湘華鯪產(chǎn)卵需要的生境條件將更加嚴格，這需對目標魚類棲息偏好條件進行更深入的研究。

參考文獻

[1]黃月群，蔡德所，楊培思，等.漓江枯水期流量變化對魚類棲息地的影響模擬[J].生態(tài)科學，2018，37（2）：147-152.

[2]崔玲，陳余道，蔣亞萍.River2D模型及其在漓江桂林市區(qū)段的初步應(yīng)用[J].廣西水利水電，2007（3）：6-9.

[3]譚燕平，王玉蓉，李嘉，等.雅礱江錦屏大河灣減水河段中魚類棲息地模擬研究[J].水電能源科學，2011，29（3）：40–43.

[4]孫嘉寧，張土喬，ZHUD Z，等.白鶴灘水庫回水支流的魚類棲息地模擬評估[J].水利水電技術(shù)，2013，44（10）：17–22.

[5]孫瑩，牛天祥，王玉蓉，等.基于地形重塑的魚類棲息地模擬修復(fù)設(shè)計[J].環(huán)境影響評價，2015，37（3）：29–32.

[6]李平，梁川.洪泛區(qū)林地的開發(fā)對魚類棲息地面積的影響模擬[J].長江流域資源與環(huán)境，2016，25（8）：1286–1291.

[7]王賓賢，劉素文，田習初.湘華鯪的生態(tài)調(diào)查、人工繁殖與苗種培育的初步觀察[J].湖南水產(chǎn)科技，1977（1）：1-9.

[8]王賓賢，劉素文，田習初.湘華鯪的生物學研究[J].1982，6（4）：455-470.

[9]卞偉.湘華鯪的生物學特征及資源動態(tài)研究[J].水生態(tài)學雜志，2011，32（4）：67-73.

[10]馮廣朋，李鐘杰，謝從新，葉少文.湖北牛山湖小型魚類的群落結(jié)構(gòu)及多樣性[J].湖泊科學，2006（03）：299-304.

[11]杜浩，班璇，張輝，危起偉，陳大慶.天然河道中魚類對水深、流速選擇特性的初步觀測——以長江江口至涴市段為例[J].長江科學院院報，2010，27（10）：70-74.

[12]王銳，李嘉.引水式水電站減水河段的水溫、流速及水深變化對魚類產(chǎn)卵的影響分析[J].四川水力發(fā)電，2010，29（02）：76-79.

[13]楊慶，胡鵬，楊澤凡，儲立民，楊靜.草魚洄游的適宜流速條件與適應(yīng)閾值[J].水生態(tài)學雜志，2019，40（04）：93-100.

[14]余國安，王兆印，張康，等.人工階梯深潭改善下切河流水生棲息地及生態(tài)的作用[J].水利學報，2008，39（2）：162-167.

[15]蒲玉芳，吳琳玲，王泰然，等.引水式電站減水河段的第二生態(tài)系統(tǒng)恢復(fù)研究[J].甘肅科技，2013，29（16）：22-24.

[16]Tennant D L. Instream flow regimes for fish， wildlife，recreation and related environmental resources [J]. Fisheries，1976，1（4）：6-10.

[17]Bockelmann BN， Fenrich EK， Lin B， et al. Development of an ecohydraulics model for stream and river restoration[J]. Ecological Engineering，2004，22（4）：227-235.

[18]Roni P， Bennett T， Morley S，et al. Rehabilitation of bedrock stream channels： the effects of boulderweir placement on aquatic habitat and biota[J]. River Research and Applications，2006，22（9）：967-980.

[19]Diego Garcia de Jalon， Javier Gortazar. Evaluation of instream habitat enhancement options using fish habitat simulations： case - studies in the river Pas（ Spain） [J]. Aquatic Ecology，2007，41（3）：461-474.

[20]Bovee K D. Development and evaluation of habitat suitabilitycriteria for use in the instream flow incremental methodology[R]. Washington： United States Fish and Wildlife Service，1986.

收稿日期：2020-04-11

作者簡介：翟立磊（1986-），男，漢族，本科學歷，工程師。