遼河干流石佛寺至盤山閘段河道安全度汛分析

孫成強

(遼寧省水利水電勘測設計研究院有限責任公司，遼寧 沈陽 110006)

遼河石佛寺至盤山閘段河道全長271.6km，河道行洪寬度1060～4720m。兩岸河灘地均為遼河沖積而成，土地肥沃。從二十世紀五六十年代開始，遼河河灘地農業耕種便已形成規模，逐步建成大批民堤和圍堤。二十世紀八十年代遼河整治對全河開展了清障，但實施效果不佳，主行洪區至大堤間套堤依然存在，嚴重影響行洪[1-4]。

目前主行洪通道內所有高于灘面的套堤已全部拆除；遠期清障目標為清除河道范圍內所有套堤，考慮到遼河主行洪通道以外至堤防間套堤基本均為封閉的套堤，套堤內除林地為護堤防浪林權屬為水利以外，其余全部為兩岸居民耕地，耕地多數為集體耕地，部分為基本農田，地方政府和群眾對套堤較為依賴，致使歷次清障均難以達到效果[5-8]。從當下實際情況看，盡管清障從法律和政策層面是合理合法的，但是實際實施可行性很低。考慮到近期完成河道清障難度較大，在河道清障完成前，通過河道安全度汛分析，為下游群眾做好預警預報，指導當地各級政府組織好逃避險工作，具有較大的現實意義，也為遼河干流防汛工作提供可靠的技術支撐。

本次研究首先通過河道洪水演進數學模型分析河段套堤的保證流量，得出確保絕大部分套堤安全的保證流量；再通過石佛寺水庫調洪演算，分析在保證石佛寺水庫下泄流量不大于此套堤保證流量的情況下水庫的最大入庫流量。對小于此流量的洪水進行調節，在近期洪水調度中盡量降低套堤內耕地損失，為遼河干流安全度汛提供參考。

1 分析方法

本次分析計算采用MIKE11數學模型，遼河干流除河口局部段受潮水影響比較大，其余段均為洪水控制區域，遼河河道較寬，發生洪水時流心線比較單一順直，水流橫向運動較少，因此計算模型可采用一維數學模型[9-10]。采用2010年實測洪水位校正后的數學模型以及2013年灘地高程(大斷面位置)成果，根據本次現場調查及測量結果，對于數學模型地形進行率定，推算不同流量級下(恒定流)各河道斷面洪水位。

遼河套堤從二十世紀五六十年代開始建設，逐漸形成大大小小的民堤和圍堤，但由于建設年代久遠，通過現場踏勘大部分套堤存在破堤現象，當地百姓為了便于農業生產，在套堤多處破堤形成道口，另外由于近年對遼河套堤開展了1050線內清障任務，線外則保留原有套堤，因此出現部分套堤斷堤現象，套堤封閉區域被打開，線內保護區域采取了封育措施，將其做為河灘地，納入主行洪區。本次分析過程中，為客觀評價套堤保證水位及相應保證流量，設定將穿越套堤的道口及缺口在防汛期間采取緊急措施，進行臨時封堵，以右岸55#套堤(R131/后皮屯)為例，套堤全長7.3km，實測成果中有1處43.5m斷堤，10處穿堤道口，高差約2m，該種情況濾除10個道口高程點及1處斷堤高程點，選取其余處最低點做為套堤的堤頂高程；對于現狀套堤出現大段堤頂高差1～2m，以右岸58#套堤(R122/季家房)為例，套堤全長1.7km，實測套堤中1+164-1+700處套堤高差約2.0m，長度500m，該種情況下堤頂高程以較低段套堤頂高程為準。

結合歷次洪水套堤破壞情況，綜合分析套堤筑堤土料、尺寸、建設年代等影響套堤安全的因素，選定套堤安全加高為0.3m，推算出各套堤保證水位，利用上述建立完成的河道數據模型，確定各套堤保證流量。

2 成果分析

2.1 現狀套堤保證流量分析

石佛寺水庫至盤山閘段套堤保證流量多集中在1000～2000m3/s，其中1000m3/s以下河段占比8.5%；1000～2000m3/s河段占比51.7%；2000～3000m3/s河段占比23.7%；3000m3/s以上河段占比16.1%。

石佛寺壩下-柳河口段區間保證流量為975～6200m3/s，左、右兩岸均存在5處卡點；柳河口-卡力馬段區間保證流量為1156～3300m3/s，左、右兩岸各存在1處卡點；卡力馬-盤山閘段區間保證流量為673～2600m3/s，左、右兩岸均存在5處卡點。區間保證流量呈遞減趨勢，與區間沿程洪峰折減規律相一致，說明本段區間保證流量較為合理。石佛寺水庫至盤山閘段不同河段保證流量成果見表1。

本次分析過程中，兼顧上下游、左右岸，結合實際地形，綜合考慮在卡點內選取適當的保證流量做為安全泄量值。根據上述22處卡點的保證流量分布，篩選出保證流量小于1000m3/s的卡點，按照保證流量由小到大的順序進行排列，排序結果見表2。結合實測地形圖得出，前6處套堤均處于1050線內，本次分析可將該部分區域視為灘地行洪區，暫不列入保護范圍，而54-1#套堤位于最上游斷面，且不在1050線內，屬于封閉完好套堤，因此將其保證流量975m3/s做為本次分析的安全泄量值。

表1 石佛寺水庫至盤山閘段不同河段保證流量成果表 單位：m3/s

表2 石佛寺水庫至盤山閘段套堤(保證流量小于1000m3/s)成果表

2.2 石佛寺水庫安全度汛調度分析

調洪邊界條件：①空庫迎洪，水庫最大下泄量不超過本次分析的套堤安全泄量值975m3/s。②水庫最高洪水位不超過50.22m。(此水位為水庫100a一遇設計洪水位，決定水庫的壩頂高程。)

通過調洪計算，石佛寺水庫壩址4a一遇洪水(洪峰流量為1310m3/s)剛好滿足上述邊界條件要求；5a一遇以上洪水不滿足上述邊界要求。4a一遇洪水最高庫水位50.08m，小于100a一遇設計洪水位50.22m，滿足壩體自身安全要求，調洪成果見表3，水庫閘門控制運用方式見表4。

表3 石佛寺水庫不同頻率洪水調節成果

表4 2018年度石佛寺水庫汛期調度運行方式

3 結論

(1)本次分析利用2011年實測河道大橫斷及2013年實測套堤測量成果，對2010年洪水進行復核，建立河道數學模型，推算各套堤最大安全泄流量，在兼顧上下游、左右岸套堤的前提下，結合實際地形，綜合考慮提出安全泄量建議值為975m3/s。

(2)根據水文預報，當預報洪水頻率為4a一遇及以下洪水時，按照表4的石佛寺水庫閘門運用方式下泄洪水，在保證壩體自身安全的前提下，可保證下游套堤安全；當預報洪水為5a一遇及以上洪水時，無法保證下游套堤安全。

(3)在水庫泄洪過程中由于洪水具有不確定性，為保證水庫壩體自身安全，當庫水位上升到50.10m，必須加大泄流量，適當多開泄洪閘門，使下泄流量等于或略大于此時入庫流量，保證水庫水位不再升高。

(4)本次采用一維數學模型MIKE11模擬計算套堤影響下的河道過流能力，對局部段地形復雜，套堤分布較多河段考慮不足。今后可進一步細化方案，采用MIKE FLOOD一二維耦合分析計算。