黃石盤水庫地表水環境中水溫影響預測評價

趙得珺

（1.四川省水利水電勘測設計研究院，四川 成都 610072；2.長沙理工大學水利工程學院，湖南 長沙 410114）

0 前言

四川省巴中市黃石盤水庫位于渠江右岸巴河水系南江支流的恩陽河中游河段，壩址位于巴中市恩陽區兩路村與高岸村境內。水庫工程所在恩陽河干流水域規模為中河。水庫具有年調節性能，對河流水文情勢、水溫改變較大。由于恩陽河上沒有水溫監測資料，采用鄰近流域渠江主源州河上的東林水文站資料擬合黃石盤的入庫和壩址水溫。因為水生生物繁殖對環境要求較高，必須滿足一定的水溫、水位、流速、流態、流程等水文條件才能完成繁殖和孵化[1]。產卵時除要求達到一定水溫外，還需要一定的漲水刺激。為了保護黃石盤水庫區域水生生物的繁殖，必須預測評價水溫的影響。

本文主要研究的就是黃石盤水庫工程對水溫的影響，計算出水溫的分布，以及提出對水溫影響減緩的措施。

1 水庫水溫結構

水庫水溫結構類型分為分層型、過渡型、混合型，利用兩種方法研究。1）α-β法，結論是黃石盤水庫為分層型水庫；2）利用密度佛汝得數法，結論是黃石盤水庫是過渡型水庫。最后綜合判斷為季節分層型結構庫水交換次數的方法僅考慮了庫容與來流量，用此類經驗判別法判斷水庫水溫結構類型時考慮的因素比較單一，因此上述判定結果僅可作為是否開展進一步計算的依據，水庫的逐月分層過程、下泄水溫等仍需要采用數學模型作進一步分析。

2 水庫水溫預測

根據《水電水利建設項目河道生態用水、低溫水和過魚設施環境影響評價技術指南（試行）》，采用二維水溫模型[2-5]對黃石盤庫區水溫進行預測。

2.1 寬度平均的立面二維水溫數學模型及其求解方法

對于常態下的水體，可忽略壓力變化對密度的影響，密度與溫度的關系可表示為：

式中：β[1/℃]為等壓膨脹系數；ρ[kg/m3]為密度；T[℃]為溫度；ρs[kg/m3]、Ts[℃]為參考狀態的密度和溫度。

2.2 模型范圍及網格劃分

計算網格單元在主流方向上尺寸為5m～400 m，在水深方向上為0.5 m～1 m，在正常蓄水位下，黃石盤庫區離散為105×59個矩形網格。

2.3 模型驗證

黃石盤水庫尚處設計階段，無法得到該水庫實測水溫資料對所用水溫數值模型進行率定與驗證。為此借用資料豐富的類似工程紫坪鋪水電站的現場水溫觀測資料，對二維水溫模型的應用性能進行檢驗。

模型中需要率定的參數只有熱通量計算中的太陽輻射表面吸收系數β1和太陽輻射在水體中的衰減系數η，分別為0.65和0.5。

2.4 取水方案

取水口設置在大壩左岸的廠房取水壩段內，進水口底板高程349 m，頂高程374.80 m，上游通過交通橋連接與左岸非溢流壩與右岸溢流壩段。取水壩段采用整體式鋼筋混凝土結構，單機單引取水，引水隧洞為有壓隧洞，全長240.61 m，洞徑D=5.6 m，縱坡i=3.07%。引水隧洞長240.61 m，取水閘始端布置2孔攔污柵閘，每孔內設攔污柵一扇，孔口8 m×4.6 m（寬×高）。其后接漸變段，寬度由20 m漸變至5.6 m。進水閘工作門1孔，孔口5.6 m×5.6 m（寬×高），閘頂高程377.40 m。閘室長26.50 m，閘墩厚5.0 m，底板厚5.0 m，采用整體式鋼筋混凝土結構。

2.5 計算條件

水動力條件：由電站運行方式所確定的設計豐、平、枯水年水庫調度過程，給出水庫各設計代表年逐月水位、流量過程，設計豐、平、枯代表年如下：

豐水代表年：1964年5月～1965年4月（P=10%）

平水代表年：1969年5月～1970年4月（P=50%）

枯水代表年：1978年5月～1979年4月（P=90%）

水流邊界條件：上游邊界給來流流量過程、下游出口給水庫下泄流量過程，校核水庫壩前水位；水溫邊界條件：上游來流水溫由鄰近水質監測站觀測的水溫插值給出，水庫水面水溫參考巴中市氣象條件。

2.6 計算工況

對設計豐、平、枯水年的水力條件，多年平均氣象條件，單層取水和疊梁門分層取水[6-8]兩種方案進行組合，形成6組計算工況。

2.7 計算結果與分析

2.7.1 壩前垂向水溫分布

以平水年為例，壩前水溫預測成果見圖1。

圖1 壩前垂向水溫分布（設計平水年）

根據水庫采用二維模型的計算結果，黃石盤水庫一般在3~8月存在一定的垂向溫差，部分月份存在較為明顯的表層溫躍層，而在9月~翌年2月混合較為均勻，據此可判斷黃石盤水庫水溫結構為過渡型。

2.7.2 下泄水溫

圖2 比較了黃石盤電站各典型水文年的下泄水溫與壩址天然水溫的差值。各工況均存在春季低溫水和冬季高溫水問題。春季低溫水方面，平水年和枯水年的下泄低溫水現象相對豐水年較突出。冬季高溫水方面，枯水年相對平水年和豐水年較突出。圖3對比了各典型水文逐日下泄水溫與壩址天然水溫。

圖3 黃石盤電站各典型水文年單層取水時逐日下泄水溫與壩址天然水溫

2.7.3 壩下河道水溫預測

黃石盤下泄水溫作為研究河段的上游邊界，采用增溫率方法計算并計入支流入匯影響可以得到黃石盤壩下的沿程水溫。為方便對比，也采用增溫率方法以壩址天然水溫為基礎計算了壩下各斷面的天然水溫[9]。

在壩下12 km（中華鱉保護區上邊界），單層取水時，與天然水溫相比，春季最大降溫2.6℃（枯水年4月）。疊梁門取水時，枯水年4月降溫0.9℃，最大溫降為枯水年3月的1.2℃，為9.3℃，春季低溫水現象得到有效緩解。

圖4是各斷面在水庫單層取水時的逐日變化過程，各斷面的低溫水降幅呈沿程減小趨勢，在保護區入口處的低溫水降幅最大。圖5是各斷面在水庫疊梁門取水時的逐日變化過程，在保護區入口處的低溫水降幅最大。

圖4 黃石盤壩下12.0 km處逐日水溫變化過程(單層取水)

圖5 黃石盤壩下12.0 km處逐日水溫變化過程(疊梁門取水)

2.8 庫區流場分析

為說明庫區垂向水溫變化與流動的關系，圖6和圖7分別繪出了平水年低溫水降幅最大的5月和高溫水升幅最大的11月的溫度場疊加流場。

圖6 黃石盤水庫平水年5月壩前12 km流場

圖7 黃石盤水庫平水年11月壩前12 km流場

5月來流水溫為20.3℃，庫表水溫為21.4℃，與來流水溫較為接近，所以來流在庫表流動，不能和庫底水體進行摻混。可見，庫尾方向的流動接近表層，而下泄水為水庫中偏下的較冷水體。11月流量為5月的2.5倍，雖水位抬升庫容增大，但來流水溫低于庫區水溫，流動層靠近庫底。受出庫影響，水流在接近壩前時垂向流速逐漸加大，部分水流沿發電孔口出庫，部分形成庫表回流。

3 水溫影響減緩措施

由于黃石盤水庫在豐、平、枯典型水文年運行時均對下游水溫存在明顯的低溫水效應，本工程考慮采用分層取水措施來減緩下泄低溫水的影響。

據《恩陽河中華鱉國家級水產種質資源保護區調整論證報告》，該河段魚類的產卵期主要集中為3月～7月，產卵水溫在16～32℃之間。黃石盤水庫采用單層取水方式時，2月水庫垂向水溫較為均勻，采用分層取水的效果極為有限，且下游無明顯的水生生物水溫需求；6月的水庫下泄水溫已充分滿足產卵水溫需求，無采取分層取水的必要性。因此綜合水庫水溫分布狀況與下游水溫需求，確定在3月～5月采用分層取水措施改善下泄低溫水的影響。

通過在發電進水口前門槽中放置不同數目的疊梁門來滿足分層取水要求，本階段初步擬定的疊梁門分層取水方案為：3.4 m+3 m×3層(共 4層)，門頂高程 372.5 m，底板高程350.6 m。表1對比了黃石盤電站疊梁門方案的下泄水溫。

表1 疊梁門分層取水效果對比（℃）

采用分層取水措施后，黃石盤水庫建成后下泄水溫與天然水溫差異較小，同時豐、平、枯水年下泄水溫在5月分別為20.4℃、20.2℃、19.8℃，滿足作物要求溫度，基本不影響作物正常生長。

4 結論

（1）采用二維水溫模型對黃石盤庫區水溫進行預測，建立模型，網格劃分，模型驗證，取水方案，計算條件和工況，計算壩前垂向水溫分布、下泄水溫、壩下河道水溫預測結果，并提出水溫影響減緩措施，最后進行庫區流場分析。

（2）本工程考慮采用分層取水措施來減緩下泄低溫水的影響。綜合水庫水溫分布狀況與下游水溫需求，確定在3月～5月采用分層取水措施改善下泄低溫水的影響。由于疊梁門可保證在不同運行水位下取得表層溫水，因此為了減緩黃石盤水庫運行對下泄水溫過程的影響，考慮在低溫水影響較為突出的3月～5月采用疊梁門分層取水方案，達到取用上層溫水減緩下泄低溫水影響的目的。采用分層取水措施后，黃石盤水庫建成后下泄水溫與天然水溫差異較小，滿足作物要求溫度，基本不影響作物正常生長。

[1]胡佳純，黃小雪，程香菊.觀音巖水庫水溫結構規律分析[J].中國水運(下半月),2016,6,P201-204.

[2]黃安陽，鄧云，脫友才，等.分層型水庫水溫水質耦合模型研究及其應用[J].人民珠江,2015,1,P8-12.

[3]江春波，周琦，劉健.河道與水庫水溫模擬及預測模型概述[J].環境影響評價,2016,3,P53-56.

[4]邱進生，鄧云，顏劍波，等.水庫水溫數學模型適用性研究[J].環境影響評價,2017,39(2),P57-62.

[5]楊學倩，朱岳明.水庫水溫計算方法綜述[J].人民黃河,2009,31(1),P41-42.

[6]李璐，陳秀銅.水庫分層取水方式減緩下泄低溫水效果研究[J].環境科學與技術,2016,2,P215-218.

[7]傅菁菁，李嘉，芮建良，等.疊梁門分層取水對下泄水溫的改善效果[J].天津大學學報,2014,7.

[8]柳海濤，孫雙科，王曉松，等.大型深水庫分層取水水溫模型試驗研究[J].水力發電學報,2012,1,P129-134.

[9]薛聯芳，顏劍波.水庫水溫結構影響因素及與下泄水溫的變化關系[J].環境影響評價,2016,3,P29-31.