蒲石河抽水蓄能電站水庫(kù)冰情研究

趙海鏡，劉書寶，張艷紅

(1.中國(guó)電建集團(tuán)北京勘測(cè)設(shè)計(jì)研究院有限公司，北京 100024；2.隆化縣水務(wù)局，河北 承德 068150)

冰凍問題是我國(guó)北方地區(qū)普遍存在的自然現(xiàn)象，目前我國(guó)正在大力發(fā)展新能源產(chǎn)業(yè)，在寒冷地區(qū)規(guī)劃了大量的抽水蓄能電站，這些抽水蓄能電站的規(guī)劃設(shè)計(jì)、運(yùn)行都不同程度地遇到冰凍問題。與常規(guī)水電站不同，抽水蓄能電站具有兩個(gè)顯著特點(diǎn)：①水庫(kù)庫(kù)容一般比較小，水庫(kù)水位漲落幅度較大；②抽水蓄能電站有發(fā)電和抽水兩種主要運(yùn)行方式，在正常運(yùn)行情況下，抽水蓄能電站上下水庫(kù)水位每天一般都要至少經(jīng)歷一個(gè)水位漲落循環(huán)。寒冷地區(qū)抽水蓄能電站水庫(kù)冰情特點(diǎn)具有和常規(guī)電站水庫(kù)冰情不同的特點(diǎn)，為解決寒冷地區(qū)抽水蓄能電站的設(shè)計(jì)和運(yùn)行中遇到的問題，有必要開展抽水蓄能電站水庫(kù)冰情研究[1]。

蒲石河抽水蓄能電站位于北緯40°25′的遼寧省丹東市寬甸滿族自治縣，所處地區(qū)氣候嚴(yán)寒，丹東氣象站最冷月(1月份)平均最低氣溫為-12.8 ℃，極端最低氣溫為-38.5 ℃。據(jù)蒲石河抽水蓄能電站水庫(kù)冰情原型監(jiān)測(cè)成果，該電站水庫(kù)冬季形成的冰蓋會(huì)侵占水庫(kù)庫(kù)容，造成發(fā)電庫(kù)容減小，并在一定程度上影響電站機(jī)組運(yùn)行，甚至對(duì)水工建筑物造成冰凍破壞[2]。本文在研究電站水庫(kù)冰情形成及消長(zhǎng)過程的基礎(chǔ)上，分析電站機(jī)組運(yùn)行和水庫(kù)冰情之間的關(guān)系，提出建議的冬季運(yùn)行方式；力圖得到蒲石河抽水蓄能電站水庫(kù)最大冰厚的計(jì)算方法，進(jìn)而建立能模擬電站水庫(kù)冬季結(jié)冰過程的數(shù)學(xué)模型，為確定水庫(kù)冰凍庫(kù)容及提出防冰害措施提供參考。

1 工程概況

蒲石河抽水蓄能電站(以下簡(jiǎn)稱“蒲石河電站”)總裝機(jī)容量為1 200 MW (4×300 MW)，屬一等大(1)型工程。主要建筑物由上、下水庫(kù)大壩，下水庫(kù)泄洪排沙閘，水道和地下廠房系統(tǒng)等工程組成。

上水庫(kù)位于長(zhǎng)甸鎮(zhèn)東洋河村泉眼溝溝首，上水庫(kù)擋水建筑物為鋼筋混凝土面板堆石壩，正常蓄水位392.0 m，相應(yīng)庫(kù)容1 135萬(wàn)m3，死水位360.0 m，死庫(kù)容95萬(wàn)m3，總庫(kù)容為1 256萬(wàn)m3。壩長(zhǎng)714.0 m，最大壩高78.5 m，壩頂寬10.0 m，壩頂高程395.5 m[3]。

下水庫(kù)位于鴨綠江右岸支流蒲石河干流下游，壩址在長(zhǎng)甸鎮(zhèn)王家街附近。下水庫(kù)正常蓄水位66m，相應(yīng)庫(kù)容2 871萬(wàn)m3，死水位62 m，死庫(kù)容1 616萬(wàn)m3，調(diào)節(jié)庫(kù)容1 255萬(wàn)m3。下水庫(kù)擋水建筑物為混凝土重力壩，擋水壩段分重力壩段、門庫(kù)壩段、導(dǎo)流底孔壩段。壩頂高程為70.1 m，壩頂全長(zhǎng)為336 m(包括泄洪排沙閘和單孔溢流壩段)，最大壩高34.1 m。重力壩段頂寬9.5 m，門庫(kù)壩段壩頂寬度16.5 m。

2 水庫(kù)冰情分析

依據(jù)蒲石河電站2013年～2014年冬季、2014年～2015年冬季和2015年～2016年冬季原型監(jiān)測(cè)資料，對(duì)該電站冰情進(jìn)行分析研究。在冰情觀測(cè)期內(nèi)，電站機(jī)組均正常運(yùn)行。

2.1 水庫(kù)冰情形成及消長(zhǎng)過程分析

2.1.1上水庫(kù)

蒲石河電站上水庫(kù)冬季結(jié)冰形態(tài)為厚冰蓋、薄冰、流冰及碎冰，冰蓋與面板之間有動(dòng)水帶出現(xiàn)，結(jié)冰起始日期在11月中、下旬至12月中旬之間，融冰結(jié)束日期在3月下旬至4月中旬之間，最大冰厚都出現(xiàn)在2月中旬左右。其中，在2013年～2014年度冬季，上水庫(kù)最大冰厚約為25 cm，動(dòng)水帶寬約2～50 cm(其中庫(kù)中心冰蓋與大壩面板之間的動(dòng)水帶寬度為2～10 m，進(jìn)/出水口順?biāo)鞣较騽?dòng)水帶寬約30～50 m)；在2014年～2015年度冬季，上水庫(kù)最大冰厚約為18 cm，動(dòng)水帶寬約3～60 cm(其中庫(kù)中心冰蓋與大壩面板之間的動(dòng)水帶寬度為3～20 m，進(jìn)/出水口順?biāo)鞣较騽?dòng)水帶寬約40～60 m)；在2015年～2016年度冬季，上水庫(kù)最大冰厚約為30 cm，動(dòng)水帶寬約2～40 cm(其中庫(kù)中心冰蓋與大壩面板之間的動(dòng)水帶寬度為2～15 m，進(jìn)/出水口順?biāo)鞣较騽?dòng)水帶寬約20～40 m)。

2.1.2下水庫(kù)和砬子溝水庫(kù)

在冰情觀測(cè)期內(nèi)，蒲石河電站下水庫(kù)冬季結(jié)冰形態(tài)為厚冰蓋，閘門前布置了射流法擾動(dòng)防冰裝置，冰蓋與下水庫(kù)閘門之間有寬約2～8 m的動(dòng)水帶出現(xiàn)。在2013年～2014年冬季、2014年～2015年冬季和2015年～2016年冬季，下水庫(kù)最大冰厚分別為65 cm、40 cm和55 cm，出現(xiàn)在1月上旬～2月上、中旬，動(dòng)水帶寬約1～3 m。下水庫(kù)進(jìn)/出水口順?biāo)鞣较驇?kù)面為無冰的水域，在2013年～2014年冬季水域順?biāo)鞣较蜷L(zhǎng)約450 m，寬為全河寬；在2014年～2015年冬季水域順?biāo)鞣较蜷L(zhǎng)約350 m，寬約為全河寬的3/4；2015年～2016年冬季水域順?biāo)鞣较蜷L(zhǎng)約200 m，寬約為河寬的2/3。結(jié)冰起始日期在11月中、下旬至12月中旬之間，融冰結(jié)束日期在3月中、下旬。

砬子溝水庫(kù)為距離大壩下游約4 km的常規(guī)水庫(kù)，2013年～2014年冬季、2014年～2015年冬季和2015年～2016年冬季，砬子溝水庫(kù)結(jié)冰形態(tài)為厚冰蓋，最大冰厚60 cm、38 cm和55 cm，約出現(xiàn)在2月中旬至3月初之間，無岸冰。結(jié)冰起始日期在11月中旬前后，融冰結(jié)束日期在4月中旬。

2.2 機(jī)組運(yùn)行和水庫(kù)冰情的關(guān)系分析

2013年～2014年冬季和2014年～2015年冬季的運(yùn)行實(shí)踐證明，只要采取包括足夠運(yùn)行次數(shù)和運(yùn)行時(shí)間的冬季運(yùn)行方案，上、下水庫(kù)均不會(huì)形成完整冰蓋，從而不會(huì)對(duì)電站運(yùn)行產(chǎn)生明顯影響。為了保證蒲石河電站冬季運(yùn)行不受冰凍影響，有必要研究電站機(jī)組運(yùn)行(包括運(yùn)行頻次、時(shí)間等)和冰情之間的關(guān)系。

2.2.1研究時(shí)段內(nèi)的電站機(jī)組運(yùn)行情況

蒲石河電站在2013年～2014年冬季和2014年～2015年冬季的最冷月，機(jī)組日均運(yùn)行次數(shù)(抽水次數(shù)和發(fā)電次數(shù)的總和)分別為9.63次和9.98次，機(jī)組日均運(yùn)行時(shí)間分別為36.3臺(tái)·時(shí)和38.6臺(tái)·時(shí)。2013年～2014年冬季和2014年～2015年冬季上、下庫(kù)水位差的變幅分別為1.82、3.23 m。

2.2.2電站機(jī)組運(yùn)行對(duì)上水庫(kù)冰情的影響

2013年～2014年冬季和2014年～2015年冬季，機(jī)組正常運(yùn)行，在未采取庫(kù)水?dāng)_動(dòng)措施的情況下，由于機(jī)組運(yùn)行、水位升降，在庫(kù)岸和冰蓋之間產(chǎn)生大面積的的冰水混合變動(dòng)帶，上水庫(kù)最大冰厚分別為25、18 cm左右。

2.2.3電站機(jī)組運(yùn)行對(duì)下水庫(kù)冰情的影響

2013年～2014年冬季和2014年～2015年冬季，機(jī)組正常運(yùn)行，下水庫(kù)進(jìn)/出水口至攔河壩庫(kù)區(qū)，結(jié)冰形態(tài)為完整的厚冰蓋，最大冰厚分別為65、40 cm；由于閘門前布置了射流法擾動(dòng)防冰裝置，冰蓋與下水庫(kù)閘門之間有寬約2～8 m的動(dòng)水帶出現(xiàn)；在下水庫(kù)進(jìn)/出水口附近出現(xiàn)大面積水域，為無冰水面，有的年份水域可直達(dá)對(duì)岸，將河上、下游冰蓋隔開。下水庫(kù)攔河壩下游約4 km處的砬子溝水庫(kù)最大冰厚分別為60、38 cm。

據(jù)目前已有的冬季運(yùn)行資料和冰情監(jiān)測(cè)資料，在不出現(xiàn)極端寒冷氣溫的條件下，在最冷月，蒲石河電站采用日均運(yùn)行次數(shù)不少于9.61臺(tái)·次、日均運(yùn)行時(shí)間不少于39.6臺(tái)·時(shí)的冬季運(yùn)行方案，上、下水庫(kù)一般不會(huì)形成完整冰蓋，對(duì)電站運(yùn)行不會(huì)產(chǎn)生明顯影響。

2.3 水庫(kù)最大冰厚計(jì)算方法

依據(jù)前述抽水蓄能電站運(yùn)行與水庫(kù)冰情的關(guān)系分析，抽水蓄能電站運(yùn)行情況會(huì)影響到一定范圍內(nèi)的庫(kù)區(qū)最大冰厚。按照最大冰厚是否受電站運(yùn)行影響，抽水蓄能電站水庫(kù)庫(kù)區(qū)可分為不受運(yùn)行影響區(qū)域(指距離進(jìn)/出水口較遠(yuǎn)，冰厚受水位變化影響可忽略不計(jì)的區(qū)域)、受水庫(kù)運(yùn)行影響區(qū)域(指距離進(jìn)/出水口較近，結(jié)冰厚度較小、冰厚受水位變化影響不可忽略的區(qū)域)和水域(指進(jìn)/出水口附近受水位變化影響較大的不結(jié)冰區(qū)域及其他不結(jié)冰區(qū)域)三部分。

經(jīng)過篩選分析，抽水蓄能電站水庫(kù)冰厚主要受到機(jī)組運(yùn)行要素(運(yùn)行次數(shù)、運(yùn)行時(shí)間、上下庫(kù)水位差)、氣溫、水溫、地溫、輻射等因素影響。

2.3.1下水庫(kù)最大冰厚計(jì)算方法

蒲石河電站下水庫(kù)最大冰厚位于不受水位變化影響區(qū)域，最大冰厚采用下式計(jì)算[4]。即

(1)

式中，δip為寒冷地區(qū)抽水蓄能電站水庫(kù)冰厚，m；φi為冰厚系數(shù)，一般可取0.022～0.026(嚴(yán)寒地區(qū)宜取大值)；Im為歷年最大凍結(jié)指數(shù)，℃·d。

2.3.2蒲石河抽水蓄能電站上水庫(kù)最大冰厚計(jì)算方法

據(jù)已有實(shí)測(cè)資料，蒲石河電站上水庫(kù)全庫(kù)面最大冰厚均受到電站運(yùn)行影響。結(jié)合實(shí)際收資情況，歸納蒲石河電站2013年～2014年冬季、2014年～2015年冬季和2015年～2016年冬季水庫(kù)冰情原型監(jiān)測(cè)資料，對(duì)冰厚δip(m)和日均運(yùn)行次數(shù)Nr(次)、日均運(yùn)行時(shí)間Tr(h)、氣溫Ta(℃)、水溫Tw(℃)及上下庫(kù)水位日變幅的絕對(duì)值|ΔH|(m)開展多元回歸分析，對(duì)蒲石河電站上水庫(kù)的19組資料進(jìn)行多元回歸分析，得到與蒲石河電站上水庫(kù)條件類似的受電站運(yùn)行影響區(qū)域的水庫(kù)最大冰厚計(jì)算關(guān)系式

δip=1.84-0.023 1lnNr-0.0193 lnTr-0.427 8ln(Ta+50)-0.027 1lnTw-0.017 0ln|ΔH|

(2)

由式(2)可知：各自變量前的負(fù)號(hào)表示冰厚δip隨日均運(yùn)行次數(shù)Nr、日均運(yùn)行時(shí)間Tr、氣溫Ta、水溫Tw及上下庫(kù)水位日變幅的絕對(duì)值|ΔH|中每一個(gè)參數(shù)的增大而減小；各自變量前的參數(shù)絕對(duì)值大小反映了各自變量對(duì)因變量的影響程度。即，各影響因素對(duì)冰厚δip的影響程度由大到小排序依次是氣溫Ta、水溫Tw、日均運(yùn)行次數(shù)Nr、日均運(yùn)行時(shí)間Tr和上下庫(kù)水位日變幅的絕對(duì)值|ΔH|。該式系根據(jù)蒲石河電站3個(gè)年度冬季的原型監(jiān)測(cè)資料分析得到，所依據(jù)的資料數(shù)量及資料的應(yīng)用范圍有限，其適用性有待于在后期工程應(yīng)用中驗(yàn)證和改進(jìn)。

3 蒲石河電站水庫(kù)冰情數(shù)值模擬

蒲石河電站水庫(kù)冰情數(shù)學(xué)模型由水力計(jì)算模塊、溫度擴(kuò)散模塊、冰蓋熱力消長(zhǎng)模塊3個(gè)模塊組成。模型初步具備人機(jī)互動(dòng)的交互式對(duì)話框界面，并具有圖形顯示庫(kù)區(qū)冰厚模擬結(jié)果的功能，模擬結(jié)果主要包括庫(kù)區(qū)最大冰厚和冰凍庫(kù)容值。

3.1 水力計(jì)算模塊

水力計(jì)算模塊的模型控制方程如下：

(3)

動(dòng)量方程

(4)

(5)

式中，h為單元水深，m；Z為水位，Z=Z0+h，m；Z0為床面高程，m；S為源匯項(xiàng)，m/s；qx、qy分別為x、y方向上的單寬流量，且qx=hu，qy=hv，m3/(s·m)；u，v分別x、y方向上的平均流速，m/s；n為糙率；g為重力加速度。

抽水蓄能電站庫(kù)區(qū)在進(jìn)水和出水過程中庫(kù)區(qū)水位整體上升或下降。因此，在模擬計(jì)算時(shí)必然涉及到動(dòng)邊界處理技術(shù)，亦即干濕床面問題。動(dòng)邊界是指平面計(jì)算區(qū)域中有水和無水區(qū)域的交界線。借鑒洪水流動(dòng)過程中的動(dòng)邊界處理方法，在計(jì)算域包含所有干濕網(wǎng)格的前提下，動(dòng)邊界處理技術(shù)主要考慮將整個(gè)計(jì)算區(qū)域內(nèi)的網(wǎng)格都參與計(jì)算，并用“水位平鋪”方法，將動(dòng)邊界問題轉(zhuǎn)化為網(wǎng)格干濕問題[5]。需要指出的是，由于所有網(wǎng)格都參與了計(jì)算，對(duì)于網(wǎng)格數(shù)目龐大的模擬計(jì)算，應(yīng)盡量將水流始終不能到達(dá)的網(wǎng)格剔除，從而減少計(jì)算上的浪費(fèi)。

3.2 溫度擴(kuò)散模塊

計(jì)算熱傳導(dǎo)及溫度擴(kuò)散問題，需求解能量守恒方程[6]

(6)

式中，Cp為比熱容，水的比熱容一般取值4 200 J/(kg·℃)；T為溫度，℃；k為流體的導(dǎo)熱系數(shù)，一般水體的導(dǎo)熱系數(shù)為0.55 J/(m·℃·s)左右；ST為流體的內(nèi)熱源及由于粘性作用流體機(jī)械能轉(zhuǎn)換為熱能部分，J。

將式(6)轉(zhuǎn)換為平面二維控制方程

(7)

3.3 冰蓋熱力消長(zhǎng)模塊

抽水蓄能電站的冰厚日變化值Δ(hi)d可視作為冰蓋和大氣、冰蓋和水體以及冰蓋和太陽(yáng)輻射的熱交換聯(lián)合影響下的結(jié)果，因此有

(8)

將冰厚日變化值按照冰期天數(shù)求和，即可得到冰期任一天的冰蓋厚度數(shù)值解。式(8)即為本研究提出的抽水蓄能電站冰蓋熱力消長(zhǎng)模塊基本控制方程。

3.4 冰情數(shù)學(xué)模型實(shí)例應(yīng)用計(jì)算結(jié)果分析

3.4.1上水庫(kù)

依據(jù)蒲石河電站上水庫(kù)冰情實(shí)測(cè)資料，取2014年～2015年冰期為實(shí)例應(yīng)用研究時(shí)段。自2014年12月1日起，氣溫轉(zhuǎn)負(fù)，因此將該日期選為冰期計(jì)算的起始日期，此外將2015年1月19日定義為冰期計(jì)算域最末一天，整個(gè)計(jì)算日期共計(jì)59 d。數(shù)學(xué)模型計(jì)算結(jié)果見圖1。

圖1 蒲石河上水庫(kù)冰期50 d模擬計(jì)算庫(kù)區(qū)冰厚結(jié)果

由數(shù)值模擬成果可知，進(jìn)/出水口附近區(qū)域在整個(gè)計(jì)算時(shí)段內(nèi)未有冰蓋形成，庫(kù)區(qū)遠(yuǎn)離進(jìn)出水口的南部區(qū)域的冰蓋厚度則由薄變厚，逐日遞增，且冰厚的增長(zhǎng)速率由塊轉(zhuǎn)慢，至計(jì)算時(shí)段后期已基本趨于穩(wěn)定。上述模擬結(jié)果與原型觀測(cè)結(jié)果基本一致。

3.4.2下水庫(kù)

下水庫(kù)實(shí)測(cè)氣象資料與上水庫(kù)略有差別，本次模擬計(jì)算冰期期間(2016年12月12日～2017年1月19日)，數(shù)學(xué)模型計(jì)算結(jié)果見圖2。

圖2 蒲石河下水庫(kù)冰期39 d模擬計(jì)算庫(kù)區(qū)冰厚結(jié)果

由圖2可見，下水庫(kù)進(jìn)/出水口區(qū)域附近在整個(gè)計(jì)算時(shí)段內(nèi)均無冰蓋形成，這一顯著特征與原型觀測(cè)成果保持一致。

4 結(jié) 語(yǔ)

本文在總結(jié)蒲石河電站水庫(kù)冰情形成及消長(zhǎng)過程、統(tǒng)計(jì)得出電站機(jī)組運(yùn)行和水庫(kù)冰情的關(guān)系的基礎(chǔ)上，得到了如下成果：

(1)提出了蒲石河電站類似工程建議的冬季運(yùn)行方式：冬季日均運(yùn)行次數(shù)不少于9.61臺(tái)·次、日均運(yùn)行時(shí)間不少于39.6臺(tái)·時(shí)，上、下水庫(kù)一般不會(huì)形成完整冰蓋，對(duì)電站運(yùn)行不會(huì)產(chǎn)生明顯影響。

(2)分析了蒲石河電站氣象條件對(duì)冰情的影響因素，提出了與該電站條件類似工程冬季最大冰厚的計(jì)算方法：下水庫(kù)最大冰厚按式(1)計(jì)算，上水庫(kù)最大冰厚按式(2)計(jì)算。

(3)建立了抽水蓄能電站冰情數(shù)學(xué)模型并在蒲石河電站上、下水庫(kù)進(jìn)行了實(shí)例應(yīng)用，冰情模型計(jì)算結(jié)果和原型監(jiān)測(cè)成果基本一致。

本文所得成果可供同類工程借鑒，其適用性有待于在后期工程應(yīng)用中驗(yàn)證和改進(jìn)。