基于MIKE21的變電站洪澇災害評估研究

劉 強，孫才華

(國網河南省電力公司，河南 鄭州 450052)

電力行業是對氣象高敏感和高需求性的行業。近年來，隨著氣候變暖加劇，全球水循環加速，頻繁發生的極端天氣氣候事件引發的氣象災害不斷加劇[1-2]，嚴重威脅電網的正常穩定運行。

河南電網居于華中、華北、西北聯網樞紐地位，在“西電東送、南北互供、全國聯網”格局中處于重要位置，河南省共有覆蓋全部地市的500 kV變電站47座，特高壓換流站4座，變電站的安全運行對電網至關重要。然而，隨著城市的發展，早期位于城區遠郊的變電站逐漸已被城市包圍，隨著房屋、道路的修建，變電站周邊下墊面條件發生巨變，其水文、排水、徑流等條件與建站初期已經發生較大變化，在發生設計標準洪水時，變電站存在發生站外洪澇影響及站內積水的風險，抵御洪澇水的能力具有較大不確定因素，存在安全隱患。

在上述背景下，以助力精準防汛為出發點，以國家電網朝歌500 kV變電站為例進行洪澇災害評估研究工作，以期滿足現有防汛管理工作的要求，為運維管理人員的組織、安排、決策提供了支撐。

1 變電站及區域自然地理概況

1.1 變電站概況

朝歌變電站位于鶴壁市東南部浚縣境內，西北距鶴壁淇濱區約9.0 km，東距浚縣城區約14.0 km，西距鉅橋鎮約7.0 km，北距S305省道約250 m；站址西南距離淇河約9.0 km，東南距共渠西滯洪區、共產主義渠分別為2.1、4.9 km。站址區域位于低山丘陵區的崗坡地上，地形平坦開闊，地勢由西北向東南傾斜，自然地面高程+86～+92 m，區域自然條件下排水通暢，自然排水方向為西北至東南。變電站及門口管涵位置如圖1所示。

圖1 變電站及門口管涵位置Fig.1 Schematic diagram of substation location

1.2 區域自然地理及水文氣象概況

鶴壁市地處豫北地區，為海河流域，東、西、北三面與安陽市相連，南面與新鄉市接壤，屬太行山與平原過渡的淺山丘陵地區。地勢西高東低、北高南低，市區以東海拔+100～+200 m；以西為低山丘陵，海拔+800 m以下；以北為丘陵地區，海拔+200 m以上；南部以平原為主。站址區域位于鶴壁市淇濱區與淇縣、浚縣之間。流經該區域內主要河流有淇河、思德河，皆為衛河(共產主義渠)支流，其中對站址防汛工作影響較大的是淇河、共產主義渠及共產主義渠西滯洪區,位置關系如圖2所示。鶴壁市屬暖溫帶半濕潤型季風氣候，由于受地形和季風影響，氣候地區性差異較大。該區的氣候特點是：四季分明、雨熱同季，冬季干冷雨雪少，春季干旱風沙多，夏季炎熱雨充沛，秋季氣爽季節短。多年平均氣溫14.2 ℃，極端最低氣溫-17.5 ℃，極端最高氣溫42.2 ℃，多年平均降水量636.5 mm，歷年最小降水量300.6 mm，最大降水量1 366.5 mm，多年平均無霜期220 d，多年蒸發量1 929 mm。因受季風氣候和地形的影響，區域降水時空分布不均，主要集中在汛期，6—9月降水量占全年降水量的74%；區域分布上西部山區大于東部平原區。

圖2 朝歌變與淇河、共產主義渠、共渠西滯洪區位置關系Fig.2 Location schematic between Zhaoge substation and Qi River，and Gongquxi basin

2 站外防澇評估

2.1 地表水二維數值模型

MIKE系列模型由丹麥水利研究所開發，用于洪水演進、城市內澇[3-4]、水質模擬[5-6]等方面。MIKE21是迄今擁有最多功能的二維水域物理、化學或生物學過程的模擬工具[7-11]。MIKE21是一個專業的工程軟件包，用于模擬河流、湖泊、河口、海灣、海岸及海洋的水流、波浪、泥沙及環境問題[12]。MIKE21水動力學模塊計算原理依據的是描述水流運動的二維非恒定流方程組，共包括3個方程[13]：水流連續性方程、水流沿x方向的動量方程及水流沿y方向的動量方程，形式如下[14-15]：

根據以上方程組，利用迭代法求解即可得到每一時刻在(x，y)處的水位z、水深h以及x、y方向的流速u、v。

2.2 模型概化

為了能夠準確模擬變電站周邊的洪水過程，從整個流域的角度分析區域內的補給、徑流、排泄條件。經過初步分析，模型的北邊界定在S26范輝高速，西邊界為G4京港澳高速和淇河，南邊界為S305省道，東邊界為共產主義渠。模型面積近200 km2，模型范圍如圖3所示。

圖3 初步模型范圍Fig.3 Preliminary model range diagram

將S26范輝高速、G4京港澳高速概化為閉邊界；由于模型范圍內淇河達不到100年一遇，故發生100年一遇洪水時，淇河按漫溢考慮，模型中將淇河概化為固定水頭邊界，水頭值為高出堤頂0.5 m，分段賦值；南邊界及東邊界地勢較低，為模型排水出口，概化為固定水頭邊界，水頭值為地面高程+0.3 m。朝歌變區域高程遠高于淇河堤頂、南邊界以及東邊界高程(>15 m)，因水位邊界影響范圍有限，不足以影響主要研究區域內的模擬成果，故不再對下邊界條件進行率定。

根據原模型的運行結果與流場的分布情況，對模型的邊界進行縮小和優化。模型西邊界根據實際地形及模擬結果取分水嶺位置；北邊界和南邊界沿流場的走向取邊界；東邊界取在共渠西滯洪區邊界位置處。將南、北以及西邊界概化為閉邊界，東邊界概化為固定水頭邊界。模型剖分控制網格邊長35 m以內，面積1 000 m2以內；對變電站周邊300 m范圍重點區域二次加密，網格邊長20 m，面積200 m2以內；對進站道路西側約10 m至東側河溝范圍二次加密，網格邊長0.5 m以內，面積0.2 m2以內。變電站周圍采用航飛地形圖高程點、其余采用河南省1∶5 000高程數據進行網格高程插值，對插值結果進行檢查并修正奇異點高程，得到模型范圍地形云圖。模型范圍內可能阻礙洪水演進的線性地物為S305省道，根據實際高程將其刻畫到模型中，道路上分布的涵洞也根據實際位置尺寸布置。村莊糙率10%，其他區域糙率6%,模型概化示意如圖4所示。

圖4 模型概化示意Fig.4 Schematic diagram of model generalization

2.3 計算方案及結果

模型設置4個計算方案，計算變電站周邊重要位置處的水位、水深、流速等，方案設置見表1，凈雨計算成果見表2。

表1 方案設置Tab.1 Scheme settings

模型共設置9處重要位置水文要素觀測點，分別為朝歌變電站8個拐角的坐標以及變電站進站門口。水位觀測點及增加擋墻如圖5所示。

經模擬計算，各方案最高水位及最大水深成果見表3，表3中水位、水深、流速均為最大值。

表2 計算區域凈雨過程Tab.2 Net rain process of computational domain

圖5 水位觀測點及增加擋墻示意Fig.5 Schematic diagram of water level observation points and added retaining wall

表3 各方案模擬成果Tab.3 Simulation results of each scheme

(1)方案1。變電站門口的排水涵洞過流能力不足，變電站門口淹沒水深達0.81 m。發生100年一遇洪水時，大門口西北側最大來水流量為1.54 m3/s，排水涵管的最大排水流量為0.52 m3/s，排水能力不足。

為解決上述問題，擬采用2種方法進行處理：①增大涵洞直徑；②于進站道路西側高地上修建擋水墻，改變進站道路西北向來流的徑流方向，使雨水沿西圍墻向南排泄。

(2)方案2。在進站道路西側高地上修建擋墻，現狀條件下，暴雨期間，變電站西圍墻和西北圍墻最大淹沒水深不超過0.15 m，考慮到修建圍墻之后水位的上漲以及水流的波動，設置擋墻高度為0.3 m，擋墻的范圍是從進站門口西側轉角一直延伸到省道S305，擋墻長度約235 m。①加入擋墻后，西北側來流沿擋墻及變電站西圍墻、西北圍墻向南流動，最后排入變電站東南側洼地。②增加擋墻后，變電站西圍墻北側水位增高0.08 m，西圍墻南側水位增高0.01 m；水位稍微增加，但增加的幅度很小。③增加擋墻之前，變電站西北角、西南角外流速分別為0.13、0.38 m/s；增加擋墻之后，變電站西北角、西南角外流速分別變為0.33、0.49 m/s。流速稍有增加，但整體流速依然很小，無需考慮對變電站圍墻的沖刷影響。

(3)方案3。當涵洞直徑增大到1.0 m時，涵洞過流能力增加，但仍不能滿足過流能力要求，變電站門口的最大淹沒水深降低為0.58 m，變電站門口積水不再沒過進站道路東側地面匯入河溝。

(4)方案4。當將涵洞的直徑增加到1.6 m，管涵過流能力滿足要求，變電站門口不再受洪水淹沒影響。

暴雨過程流場如圖6所示。

圖6 暴雨過程流場Fig.6 Flow field diagram of rainstorm process

由于增加涵洞直徑投資成本高，工程量大，協調難度大，而擋墻方案具有投資小，施工簡單，效果明顯等優點，故建議采用擋墻方案以解決暴雨期間變電站門口排水能力不足的問題。

若由于外部因素修筑擋墻難以實現，則建議對朝歌變門口現有檢查井改造成截流井，保證洪澇雨季管涵不滿足洪水流量時，洪水通過截流井溢流。截流井后設置排雨水泵池，收集截流井溢流洪水，經水泵提升排入河道，水泵總設計流量為1.5 m3/s。

3 站內洪澇監測預警

3.1 內澇算法模型定制

朝歌站內澇算法模型主要是針對朝歌站的站內排水能力、地形條件等情況進行定制的計算實時站內積水深度的算法，算法輸入為15 min頻率的降雨量數據，通過算法計算之后，輸出當前的站內累計積水深度、排水時間等數值。內澇算法流程如圖7所示。

圖7 算法流程Fig.7 Algorithm flow

3.2 監測預警標準制定

監測預警主要是針對內澇算法模型輸出的當前站內積水深度和降雨強度等結果進行監測，當監測得到的當前站內積水深度或降雨強度超過一定的閾值范圍，系統就會發出預警提示信息，預警閾值范圍的制定主要就是針對上述2個監測量。

(1)站內累計積水深度。站內累計積水深度預警分級見表4。

表4 站內累計積水深度預警分級Tab.4 Warning and classification of accumulated water depth in the station

(2)當前降雨強度。降雨強度預警分級見表5。

表5 降雨強度預警分級Tab.5 Warning and classification of rainfall

4 結論

針對近年來極端天氣和精準防汛的需求，以朝歌站為對象，基于MIKE21水動力模型進行了變電站站外暴雨時整個流域內洪澇水的發生、演進過程，得到變電站周邊區域內的水文要素的時空變化過程，從而準確獲取變電站周邊位置的主要水文數據信息。對不同計算方案進行了對比分析，并給出了站外防洪的具體方案。同時對站內內澇及站內積水進行監測預警研究，其結果可為運維管理人員的組織、安排、決策提供強有力的支撐，對其他變電站的防汛工作也具有一定的借鑒意義。