顏劍波，楚凱鋒，張德見(jiàn)，李　璜

(中國(guó)電建集團(tuán)中南勘測(cè)設(shè)計(jì)研究院有限公司，長(zhǎng)沙　410014)

兩種常用經(jīng)驗(yàn)公式法計(jì)算水庫(kù)水溫效果比較研究

顏劍波，楚凱鋒，張德見(jiàn)，李璜

(中國(guó)電建集團(tuán)中南勘測(cè)設(shè)計(jì)研究院有限公司，長(zhǎng)沙410014)

大中型水電水利工程建設(shè)引起的低溫水下泄問(wèn)題是主要的生態(tài)問(wèn)題之一。目前水庫(kù)水溫計(jì)算有經(jīng)驗(yàn)公式法和數(shù)學(xué)模型法兩大類。指數(shù)函數(shù)法和余弦函數(shù)法是目前水庫(kù)水溫計(jì)算兩種常用的經(jīng)驗(yàn)公式方法，也是現(xiàn)有規(guī)范推薦的方法。選取了我國(guó)11個(gè)不同類型的代表性水庫(kù)，利用實(shí)測(cè)資料分析研究了兩種方法在計(jì)算水庫(kù)水溫精度的差異。通過(guò)比較研究，筆者建議在有較好的類比水庫(kù)資料時(shí)采用指數(shù)函數(shù)法，否則建議采用余弦函數(shù)法。

水庫(kù)水溫計(jì)算；經(jīng)驗(yàn)公式法；指數(shù)函數(shù)法；余弦函數(shù)法

1　前　言

水電水利工程特別是大中型水電水利工程建成后，水庫(kù)在沿水深方向上呈現(xiàn)出有規(guī)律的水溫分層，大致表現(xiàn)為：冬季庫(kù)區(qū)趨于等溫分布，下泄水溫較天然情況高；春、夏季庫(kù)表水溫高，庫(kù)底水溫低，下泄水溫較天然情況低。下泄低溫水可對(duì)魚(yú)類造成突出的不利影響。水溫變冷，水體的溶氧量和水化學(xué)成分將發(fā)生變化，影響魚(yú)類和餌料生物的衍生，致使魚(yú)類區(qū)系組成發(fā)生變化。下泄低溫水將使魚(yú)類產(chǎn)卵季節(jié)推遲、影響魚(yú)卵孵化甚至造成不產(chǎn)卵。下泄低溫水還會(huì)降低魚(yú)類新陳代謝的能力，使魚(yú)生長(zhǎng)緩慢。水溫低、餌料生物生長(zhǎng)緩慢，將直接影響魚(yú)類的生長(zhǎng)、育肥和越冬。低溫水會(huì)對(duì)農(nóng)作物產(chǎn)生“冷害”影響，造成減產(chǎn)甚至絕產(chǎn)。原國(guó)家環(huán)境保護(hù)總局2005年12月在北京召開(kāi)了水電水利建設(shè)項(xiàng)目水環(huán)境與水生生態(tài)保護(hù)技術(shù)政策研討會(huì)，會(huì)議形成的“關(guān)于印發(fā)水電水利建設(shè)項(xiàng)目水環(huán)境與水生生態(tài)保護(hù)技術(shù)政策研討會(huì)會(huì)議紀(jì)要的函”(環(huán)辦函〔2006〕11號(hào))提出關(guān)于水庫(kù)水溫的計(jì)算主要有：東勘院計(jì)算方法(本文中稱指數(shù)函數(shù)法)、中國(guó)水科院方法(本文中稱余弦函數(shù)法)等經(jīng)驗(yàn)法，以及垂向一維模型、立面二維模型和三維模型等數(shù)學(xué)模擬法。

總體上，水庫(kù)水溫計(jì)算的方法有兩大類：經(jīng)驗(yàn)公式法和數(shù)值模型法[1，2]。經(jīng)驗(yàn)公式法是基于統(tǒng)計(jì)學(xué)理論建立的方法。數(shù)學(xué)模型法是基于水庫(kù)的熱量平衡機(jī)理建立的方法。理論上，數(shù)值模型法計(jì)算精度較高，但對(duì)基礎(chǔ)資料及計(jì)算參數(shù)率定要求較高[3]。在水電水利工程設(shè)計(jì)實(shí)踐工作中，經(jīng)驗(yàn)公式法仍然是水庫(kù)水溫計(jì)算的主要方法，其中應(yīng)用較多的是指數(shù)函數(shù)法和余弦函數(shù)法?，F(xiàn)行水庫(kù)水溫計(jì)算相關(guān)規(guī)范中，《水電水利工程水文計(jì)算規(guī)范》(DL/T5431-2009)[4]、《水利水電工程水文計(jì)算規(guī)范》(SL278-2002)[5]推薦指數(shù)函數(shù)法，《混凝土拱壩設(shè)計(jì)規(guī)范》(SL282-2003)[6]推薦余弦函數(shù)法。為了比較這兩種經(jīng)驗(yàn)公式法計(jì)算結(jié)果的差異，筆者選取了不同調(diào)節(jié)類型、不同水溫結(jié)構(gòu)的水庫(kù)實(shí)測(cè)資料進(jìn)行研究。

2　指數(shù)函數(shù)法簡(jiǎn)介

指數(shù)函數(shù)法計(jì)算公式如下：

(1)

(2)

(3)

y——水深，m；

m——月份。

3　余弦函數(shù)法簡(jiǎn)介

余弦函數(shù)法計(jì)算公式如下：

(4)

y——水深，m；

τ——時(shí)間，月；

τ0——?dú)鉁啬曛芷谧兓^(guò)程的初始相位，月。緯度高于30°地區(qū)，取τ0=6.5月；緯度低于或等于30°地區(qū)，取τ0=6.7月;

ω——溫度變化的圓頻率，ω=2π/12=0.524。

4　案例分析

4.1代表性工程

考慮指數(shù)函數(shù)法、余弦函數(shù)法建立時(shí)樣本資料的特征，本文選擇我國(guó)北緯30°左右11個(gè)不同類型的代表水庫(kù)進(jìn)行研究，包括不同調(diào)節(jié)性能、不同水溫結(jié)構(gòu)的電站水庫(kù)，見(jiàn)表1。

表1　研究選擇的代表性電站

4.2水溫計(jì)算結(jié)果比較分析

根據(jù)選擇的11個(gè)代表性水庫(kù)水溫實(shí)測(cè)資料，分別采用指數(shù)函數(shù)法、余弦函數(shù)法計(jì)算水庫(kù)壩前逐月水溫，計(jì)算值與與實(shí)測(cè)值誤差見(jiàn)表2。由表2可見(jiàn)，指數(shù)函數(shù)法計(jì)算水庫(kù)垂向水溫誤差基本在3℃以內(nèi)，余弦函數(shù)法計(jì)算水庫(kù)垂向水溫誤差基本在5℃以內(nèi)。不考慮選擇類比水庫(kù)的誤差，理論上指數(shù)函數(shù)法比余弦函數(shù)法有更好的計(jì)算精度。

此外，本文分析了兩種經(jīng)驗(yàn)公式計(jì)算不同類型水庫(kù)水溫分布曲線的差異，見(jiàn)圖1～圖4。其中，圖1新安江水庫(kù)代表多年調(diào)節(jié)的分層型水庫(kù)，圖2鳳灘代表年調(diào)節(jié)的分層型水庫(kù)，柘溪水庫(kù)代表季調(diào)節(jié)的過(guò)渡型水庫(kù)，富春江代表日調(diào)節(jié)的混合型水庫(kù)。從水溫垂向分布曲線形態(tài)來(lái)看，指數(shù)函數(shù)法計(jì)算的水溫分布曲線更加接近實(shí)測(cè)水溫，這種形態(tài)差異是由指數(shù)函數(shù)、余弦函數(shù)不同的插值方法決定的。

表2　指數(shù)函數(shù)法和余弦函數(shù)法計(jì)算水庫(kù)水溫誤差比較表

備注：方案A表示指數(shù)函數(shù)法，方案B代表余弦函數(shù)法?！顢?shù)量越多代表與實(shí)測(cè)水溫分布擬合度越好。

圖1　新安江水庫(kù)Fig.1　Xin’anjiang Reservoir

圖2　鳳灘水庫(kù)Fig.2　Fengtan Reservoir

圖3　柘溪水庫(kù)Fig.3　Zhexi Reservoir

圖4　富春江水庫(kù)Fig.4　Fuchunjiang Reservoir

5　小　結(jié)

水溫是評(píng)價(jià)水環(huán)境和水生生態(tài)環(huán)境變化的一個(gè)重要影響因子。水溫的變化對(duì)水質(zhì)及水生生物、農(nóng)田灌溉和生活用水等將產(chǎn)生重大影響。水庫(kù)水溫計(jì)算方法有經(jīng)驗(yàn)公式法和數(shù)學(xué)模型法兩大類。指數(shù)函數(shù)法、余弦函數(shù)法是現(xiàn)有水庫(kù)水溫經(jīng)驗(yàn)公式計(jì)算兩種主要的計(jì)算方法，也是相關(guān)規(guī)范推薦的水溫計(jì)算方法。指數(shù)函數(shù)法的本質(zhì)是先求得庫(kù)表和庫(kù)底的各月多年平均水溫，然后采用指數(shù)函數(shù)在空間維度上進(jìn)行插值。余弦函數(shù)法的本質(zhì)是先求得庫(kù)表至庫(kù)底任意深度的多年平均水溫和年變幅，然后采用余弦函數(shù)在時(shí)間維度上進(jìn)行插值。通過(guò)選擇11個(gè)不同類型水庫(kù)的實(shí)測(cè)資料分析發(fā)現(xiàn)，不考慮類比水庫(kù)的誤差，理論上指數(shù)函數(shù)法比余弦函數(shù)法具有更好的計(jì)算精度，計(jì)算水溫分布曲線形態(tài)更加貼合實(shí)際。在應(yīng)用這兩種方法進(jìn)行擬建水庫(kù)水溫預(yù)測(cè)時(shí)，如有較好的類比水庫(kù)(庫(kù)容、調(diào)節(jié)性能、地理位置等相似度高的水庫(kù))，本文建議采用指數(shù)函數(shù)法，否則建議采用余弦函數(shù)法。

[1]薛聯(lián)芳.河流水溫?cái)?shù)學(xué)模型的比較研究[J].水電站設(shè)計(jì),1998,14(4)：1.

[2]張仙娥,周孝德,等.水庫(kù)水溫研究方法述評(píng)[J].水資源與水工程學(xué)報(bào),2006,17(3):1-4.

[3]鞠石泉,蘇懷智,等.簡(jiǎn)述水庫(kù)水溫預(yù)測(cè)計(jì)算方法[J].水電能源科學(xué),2004,22(3):74-77.

[4]DL/T5431-2009,水電水利工程水文計(jì)算規(guī)范[S].

[5]SL278-2002,水利水電工程水文計(jì)算規(guī)范[S].

[6]SL282-2003,混凝土拱壩設(shè)計(jì)規(guī)范[S].

Comparative Study of the Water Temperature Calculation of Reservoirs with Two Common Empirical Formulas

YAN Jian-bo, CHU Kai-feng, ZHANG De-jian, LI Huang

(PowerchinaZhongNanEngineeringCorporationLimited，Changsha410014，China)

Low Temperature water discharging problems caused by the construction of large and medium scale hydropower and water conservancy project is one of the major ecological problems. Nowadays, there are two categories method to calculate the water temperature including experienced formula method and mathematical model. The most experienced two methods we had now is Exponential function Method and Cosine Function Method for calculating the water temperature, and they also are the most formal and recommended method. This research project selected 11 representative reservoirs in China, and used actual measuring data of these 11 reservoirs to analysis and research the precision of these two methods and different characteristics. According to our research, we found that we could use Exponential function Method if we have useful and integrated date and information of reservoir, otherwise, we suggest apply Cosine Function Method.

Prediction of water temperature in reservoirs; Empirical formula method; Exponential Function Method Method; Cosine Function Method

2015-06-11

顏劍波(1984-)，男，湖南邵陽(yáng)人，2009年畢業(yè)于河海大學(xué)環(huán)境工程專業(yè)，工程師，研究方向環(huán)境影響評(píng)價(jià)與環(huán)境保護(hù)工程設(shè)計(jì)。

X524

A

1001-3644(2015)06-0117-06