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

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

中圖分類號(hào)：X174 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：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è)對(duì)河流生態(tài)系統(tǒng)產(chǎn)生的影響主要表現(xiàn)在河流阻隔、局部河段減脫水、庫區(qū)生態(tài)系統(tǒng)的演替及下游水文時(shí)空節(jié)律改變、低溫水下泄等方面。河段減水對(duì)河流的生態(tài)環(huán)境造成極大的破壞，水生生物棲息地安全性受到脅迫。為緩解減水對(duì)河段中魚類的影響，國內(nèi)外普遍采取生態(tài)流量下泄、生境修復(fù)等措施。

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

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

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

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

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

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

2 評(píng)估方法

2.1 評(píng)估方法與模型

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

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

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

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

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

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

2.3 模擬工況

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

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

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)流量下泄時(shí)，水體流速顯著降低，較豐水期和枯水期分別下降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é)果與討論

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

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

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

5 結(jié)語

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

參考文獻(xiàn)

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

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

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

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

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

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

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

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

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

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

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

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

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

[14]余國安，王兆印，張康，等.人工階梯深潭改善下切河流水生棲息地及生態(tài)的作用[J].水利學(xué)報(bào)，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-），男，漢族，本科學(xué)歷，工程師。