孟建川 黃渝桂 楊 樂

（中水淮河規(guī)劃設(shè)計(jì)研究有限公司 合肥 230088）

1 引言

在現(xiàn)狀水利工程下，通過水庫與樞紐聯(lián)合調(diào)度，利用河道強(qiáng)迫行洪和最大限度優(yōu)化洪水安排，對(duì)防洪抗洪救災(zāi)決策的快速性和有效性有一定影響，是探討現(xiàn)代洪水災(zāi)害管理的方法。沂河作為沂沭泗流域防洪除澇體系的一部分，是沂沭泗水系中重要骨干河道，沿線城市有臨沂、郯城、邳州、新沂，經(jīng)濟(jì)社會(huì)發(fā)展水平較高，一旦發(fā)生洪水災(zāi)害，將嚴(yán)重危害人民生命財(cái)產(chǎn)安全，社會(huì)影響較大。

2 概述

沂河是沂沭泗水系中最大的山洪河道，河道全長333km，流域面積11820km2。沂河在彭家道口向東辟有分沂入沭水道，分沂河洪水入沭河；在江風(fēng)口辟有邳蒼分洪道，分沂河洪水入中運(yùn)河。

沂河防洪體系包括上游6 座大中型水庫（總控制流域面積5375km2，占臨沂、大官莊以上控制面積的39%）、沂河堤防、控制性樞紐、分洪河道及分沂入沭水道以北應(yīng)急處理區(qū)等。

沂河祊河口以下已按50年一遇防洪標(biāo)準(zhǔn)治理，臨沂至劉家道口、劉家道口至江風(fēng)口、江風(fēng)口至入駱馬湖口段設(shè)計(jì)流量分別為16000m3/s、12000m3/s、8000m3/s。分沂入沭水道已按50年一遇防洪標(biāo)準(zhǔn)治理，設(shè)計(jì)流量4000m3/s。邳蒼分洪道已按50年一遇防洪標(biāo)準(zhǔn)治理，江風(fēng)口閘下至東泇河口設(shè)計(jì)流量4000m3/s、東泇河口以下設(shè)計(jì)流量5500m3/s[1][2]。

通過對(duì)沂河超標(biāo)準(zhǔn)洪水（超50年一遇洪水）調(diào)度運(yùn)用方式的研究，分析其對(duì)沂河洪水安排的影響，從而為沂河100年一遇設(shè)計(jì)洪水調(diào)度預(yù)案制定提供依據(jù)。

3 計(jì)算方法

3.1 水庫調(diào)算

水庫進(jìn)行調(diào)洪調(diào)算，理論方法是基于水力學(xué)的圣維南方程組進(jìn)行簡化[8]，忽略了洪水入庫至泄洪建筑物間的行進(jìn)時(shí)間、沿程流速變化及動(dòng)庫容等的影響，簡化后水庫調(diào)洪計(jì)算的公式即水量平衡方程[4][6]。

式中：

Q1，Q2—分別為計(jì)算時(shí)段初、末的入庫流量（m3/s）；

q1，q2—分別為計(jì)算時(shí)段初、末的下泄流量（m3/s）；

V1，V2—分別為計(jì)算時(shí)段初、末的水庫蓄水量（m3）；

△V—V1與V2之差；

△t—計(jì)算時(shí)段。

3.2 河道演進(jìn)

河道演進(jìn)采用馬斯京根流量演算法的改進(jìn)方法即分段連續(xù)演算法。分段連續(xù)演算法根據(jù)馬斯京根演算方程為線性系統(tǒng)的特點(diǎn)，首先推求上游斷面進(jìn)入一個(gè)單位水量后經(jīng)多個(gè)河段連續(xù)向下游演進(jìn)，在下游斷面形成一個(gè)相應(yīng)的流量過程線。

沂河上游為山區(qū)性河流，干、支流間相互干擾作用不大，可把干、支流各河段視為相互獨(dú)立的無支流河段，求得各自的流量演算參數(shù)，分別把上游站（水庫出流）的入流量演算到下游站（臨沂站），然后疊加即為該斷面的出流過程[3][5][7]。

4 水庫調(diào)度辦法及調(diào)洪成果

4.1 設(shè)計(jì)洪水

沂河臨沂站洪峰及洪量計(jì)算。洪量頻率3d、7d、15d、30d 洪量系列年限與洪峰流量相同。其中1730年洪量用峰量關(guān)系插補(bǔ)，峰量關(guān)系的高水部分無實(shí)測點(diǎn)據(jù)，外延精度差，推得的1730年洪量僅作參考[3]。主要依據(jù)1912年以來的洪水點(diǎn)據(jù)，采用P—Ⅲ型曲線適線，各時(shí)段洪量頻率計(jì)算成果見表1。

表1 沂河（臨沂站）設(shè)計(jì)洪峰、洪量復(fù)核成果表（單位：洪峰m3/s；洪量億m3）

4.2 水庫調(diào)度辦法

沂河臨沂、沭河大官莊來量是按照沂、沭河洪水經(jīng)水庫調(diào)蓄后下泄，匯區(qū)間洪水，并經(jīng)河道演進(jìn)的方法確定的。根據(jù)2019年《淮河流域大型水庫汛期調(diào)度運(yùn)用計(jì)劃》中水庫調(diào)度運(yùn)用辦法進(jìn)行水庫調(diào)洪。

4.3 水庫調(diào)洪水成果

沂河洪水仍采取盡可能東調(diào)原則，洪水安排方案：當(dāng)沂河遇到100年一遇洪水時(shí)，上游水庫按照《淮河流域大型水庫汛期調(diào)度運(yùn)用計(jì)劃》中水庫調(diào)度運(yùn)用辦法進(jìn)行調(diào)洪計(jì)算，匯區(qū)間洪水，并經(jīng)河道演進(jìn)到臨沂、大官莊；沂河臨沂、沭河大官莊洪水按照2012年《沂沭泗河洪水調(diào)度方案》運(yùn)用辦法進(jìn)行安排。調(diào)洪成果中關(guān)鍵控制點(diǎn)（沂河臨沂站、沭河大官莊站）流量成果表見表2。

表2 沂河臨沂站、沭河大官莊站流量成果表

5 調(diào)洪計(jì)算及成果分析

5.1 調(diào)洪原則

根據(jù)《沂沭泗河洪水調(diào)度方案》（以下簡稱《調(diào)度方案》），沂河、沭河洪水盡可能東調(diào)，預(yù)留駱馬湖部分蓄洪容積和新沂河部分行洪能力接納南四湖及邳蒼地區(qū)洪水。遇標(biāo)準(zhǔn)內(nèi)洪水，合理利用水庫、水閘、河道、湖泊等，確保防洪工程安全。遇超標(biāo)準(zhǔn)洪水，除利用水閘、河道強(qiáng)迫行洪外，并相機(jī)利用滯洪區(qū)和采取應(yīng)急措施處理超額洪水，地方政府組織防守，全力搶險(xiǎn)，確保南四湖湖西大堤、新沂河大堤等重要堤防和濟(jì)寧、臨沂、徐州、宿遷、連云港等重要城市城區(qū)的防洪安全，盡量減輕災(zāi)害損失。調(diào)度運(yùn)行辦法表見表3。

表3 劉家道口樞紐、大官莊樞紐的調(diào)度運(yùn)行辦法表

5.2 計(jì)算情景

沂河發(fā)生100年一遇設(shè)計(jì)洪水時(shí)，經(jīng)水庫調(diào)洪后臨沂斷面最大洪峰流量為19389m3/s，通過劉家道口樞紐、大官莊樞紐進(jìn)行洪水安排。根據(jù)現(xiàn)狀河道行洪能力、樞紐調(diào)度情況等因素分四種計(jì)算情景的洪水安排：一是現(xiàn)狀洪水調(diào)度方案；二是沂河江風(fēng)口以下河道強(qiáng)迫行洪（強(qiáng)迫行洪1、2 方案）；三是優(yōu)化水庫調(diào)度運(yùn)行辦法。計(jì)算工況見表4。

表4 計(jì)算工況表

5.3 調(diào)洪計(jì)算成果

沂、沭河上游水庫按照《2019年淮河流域大型水庫汛期調(diào)度運(yùn)用計(jì)劃》中水庫調(diào)度運(yùn)用辦法進(jìn)行調(diào)洪計(jì)算，匯區(qū)間洪水，并經(jīng)河道演進(jìn)到臨沂、大官莊；沂河臨沂、沭河大官莊洪水上面四種工況進(jìn)行洪水安排。

計(jì)算情景一（現(xiàn)狀調(diào)度方案）：沂河為主，沭河相應(yīng)的沂沭河洪水，臨沂最大流量為19389m3/s，當(dāng)臨沂超過16500m3/s（劉家道口閘下泄12000m3/s，彭家道口閘下泄4500m3/s），超額洪水進(jìn)入分沂入沭以北地區(qū)采取應(yīng)急措施處理。分沂入沭最大進(jìn)洪流量2889m3/s，進(jìn)洪為3 個(gè)時(shí)段6 個(gè)小時(shí)，進(jìn)洪量為0.47 億m3。

計(jì)算情景二（強(qiáng)迫行洪1 方案）：沂河為主，沭河相應(yīng)的沂沭河洪水，臨沂最大流量為19389 m3/s，當(dāng)臨沂超過17500m3/s（劉家道口閘下泄13000m3/s，彭家道口閘下泄4500 m3/s），超額洪水進(jìn)入分沂入沭以北地區(qū)采取應(yīng)急措施處理。分沂入沭以北最大進(jìn)洪流量1889 m3/s，進(jìn)洪為3 個(gè)時(shí)段6 個(gè)小時(shí)，進(jìn)洪量為0.26 億m3。

計(jì)算情景三（強(qiáng)迫行洪2 方案）：沂河為主，沭河相應(yīng)的沂沭河洪水，臨沂最大流量為19389 m3/s，當(dāng)臨沂超過18500m3/s（劉家道口閘下泄14000m3/s，彭家道口閘下泄4500 m3/s），超額洪水進(jìn)入分沂入沭以北地區(qū)采取應(yīng)急措施處理。分沂入沭以北最大進(jìn)洪流量889 m3/s，進(jìn)洪為2 個(gè)時(shí)段4 個(gè)小時(shí)，進(jìn)洪量為0.097 億m3。

計(jì)算情景四（水庫優(yōu)化調(diào)度方案結(jié)合河道強(qiáng)迫行洪）對(duì)沂河洪水（臨沂斷面）影響最大的水庫進(jìn)行優(yōu)化調(diào)度，不影響水庫安全的前提下，適當(dāng)調(diào)整水庫運(yùn)行辦法。沂河臨沂最大洪峰流量為19389m3/s，比《淮河流域綜合規(guī)劃》遠(yuǎn)期流量的19000m3/s 略大?？紤]到跋山水庫已完成除險(xiǎn)加固工程，續(xù)建工程對(duì)跋山水庫的控制運(yùn)用辦法進(jìn)行了調(diào)整，把“運(yùn)用計(jì)劃”中庫水位178.54～179.0m水位下泄3120m3/s 調(diào)整為庫水位178.54～180.09m 下泄3120m3/s，100年一遇調(diào)洪水位為180.09m，與水庫設(shè)計(jì)洪水位180.1m 基本相同，但臨沂最大洪峰流量可減小至18331m3/s。跋山水庫調(diào)洪演算成果見表5。

表5 跋山水庫調(diào)洪演算成果表

如采用計(jì)算情景三（沂河江風(fēng)口以下段強(qiáng)迫行洪10000m3/s，邳蒼分洪道4000m3/s，分沂入沭4500m3/s）對(duì)臨沂洪水進(jìn)行安排，通過跋山水庫消峰作用，沂河為主的100年一遇洪水，分沂入沭以北應(yīng)急處理區(qū)理論上可不需要分洪。

根據(jù)四種計(jì)算情景成果分析，采用計(jì)算情景一（現(xiàn)狀調(diào)度方案）、計(jì)算情景二（河道強(qiáng)迫行洪1 方案）、計(jì)算情景三（河道強(qiáng)迫行洪2 方案），分沂入沭以北應(yīng)急處理區(qū)均需滯洪。

6 結(jié)論

（1）當(dāng)沂河發(fā)生100年一遇洪水時(shí)，根據(jù)水庫調(diào)洪計(jì)算和洪水演進(jìn)成果，臨沂斷面設(shè)計(jì)流量為19389m3/s，結(jié)合2012年國務(wù)院批復(fù)的《沂沭泗河洪水調(diào)度方案》中劉家道口樞紐、大官莊樞紐的調(diào)度辦法，分沂入沭以北應(yīng)急處理區(qū)需滯洪；如通過沂河（沂河江風(fēng)口以下8000～10000m3/s）和分沂入沭水道（4000～4500m3/s）強(qiáng)迫行洪，分沂入沭以北應(yīng)急處理區(qū)仍需滯洪；如水庫優(yōu)化調(diào)度結(jié)合沂河、分沂入沭水道強(qiáng)迫行洪（沂河江風(fēng)口以下10000m3/s，分沂入沭水道4500m3/s），分沂入沭以北應(yīng)急處理區(qū)理論上不需要滯洪。

（2）通過對(duì)沂沭河流域河庫聯(lián)動(dòng)調(diào)洪（水庫調(diào)洪、河道演進(jìn)）演算，分析劉家道口樞紐、大官莊樞紐運(yùn)用對(duì)沂河洪水的安排以及分沂入沭以北滯洪區(qū)的分洪口門設(shè)置，同時(shí)分沂入沭以北地區(qū)滯洪對(duì)區(qū)域的影響較大。因此，本次研究成果可為沂河100年一遇洪水調(diào)度預(yù)案制定提供參考