陳 鳴，王 鋼，鄭興發(fā)，鄭久興，雷四華

（1.水利部交通部國家能源局南京水利科學(xué)研究院，江蘇 南京 210029；2.水文水資源及水利工程國家重點(diǎn)試驗(yàn)室，江蘇 南京 210029；3.福建寧德市水利局，福建 寧德 352100）

0 引言

我國是水資源十分短缺的國家，人均水資源占有量不到世界人均占有量的 1/4，同時我國又是世界上洪澇災(zāi)害發(fā)生最頻繁的國家之一，有 2/3的國土面積、半數(shù)以上的人口、35% 的耕地、2/3的工農(nóng)業(yè)總產(chǎn)值受到洪水的嚴(yán)重威脅。面對亦旱亦澇的水情，研究利用現(xiàn)有水利工程的優(yōu)化調(diào)度，協(xié)調(diào)缺水與洪災(zāi)之間的矛盾，是投資小、回報高，且不破壞生態(tài)環(huán)境的方法，不失為解決我國水問題的有效途徑之一[1]。

結(jié)合計(jì)算機(jī)技術(shù)、GIS（地理信息系統(tǒng)）、河流水流仿真技術(shù)、流域產(chǎn)匯數(shù)學(xué)模型為一體的流域水庫群水資源調(diào)度系統(tǒng)，融入了全流域綜合考慮，適度承擔(dān)風(fēng)險等現(xiàn)代水庫管理、調(diào)度理念。系統(tǒng)采用圖形界面、自動優(yōu)化計(jì)算、人機(jī)交互等途徑協(xié)調(diào)流域內(nèi)具有防洪、灌溉、發(fā)電、航運(yùn)、養(yǎng)殖、城市及工業(yè)供水等多種功能大中型水庫的不同用水目標(biāo)對水庫蓄水、泄水量和時間要求所發(fā)生的沖突，可避免或減少洪災(zāi)損失，有效提高水資源利用率[2-3]。

1 系統(tǒng)結(jié)構(gòu)及功能

流域水庫群水資源調(diào)度系統(tǒng)采用 C/S（Client/Server）結(jié)構(gòu)，使系統(tǒng)的開放性得以保證?？梢郧度肫渌嚓P(guān)水動力、水文和水庫調(diào)度模型，實(shí)現(xiàn)與已有模型方法的無縫結(jié)合，對整合已有方法、提高工作效率具有明顯的促進(jìn)作用。系統(tǒng)以實(shí)時水雨情庫、歷史水文數(shù)據(jù)庫、空間數(shù)據(jù)庫為基礎(chǔ)；數(shù)據(jù)處理、模型、計(jì)算方法、表格、圖形為分析工具；GIS為交互界面。

系統(tǒng)功能主要包括：系統(tǒng)維護(hù)管理、信息查詢、水情預(yù)報與水庫調(diào)度、全球眼水情視頻。具體功能與操作菜單對應(yīng)如表1所示。

2 系統(tǒng)關(guān)鍵技術(shù)

2.1 框架結(jié)構(gòu)中的分層處理

將系統(tǒng)分為人機(jī)交互控制、計(jì)算應(yīng)用服務(wù)、數(shù)據(jù)基礎(chǔ)服務(wù)3個層次：

1）系統(tǒng)人機(jī)交互控制層（即人機(jī)交互界面），基于 GIS 技術(shù)的信息查詢應(yīng)用模式和可視化多層次的空間圖形操作界面，提供直觀、清晰、方便、靈活的系統(tǒng)操作和控制環(huán)境。

2）系統(tǒng)計(jì)算應(yīng)用服務(wù)層（即系統(tǒng)計(jì)算模型、信息服務(wù)等），依照系統(tǒng)的應(yīng)用種類進(jìn)行功能劃分，由水情信息服務(wù)系統(tǒng)、水情監(jiān)測對比分析、水情預(yù)報分析模型等部分組成。各組成部分采用模塊化開發(fā)，以利于進(jìn)行系統(tǒng)設(shè)計(jì)和開發(fā)的工程化管理，也便于在系統(tǒng)建立過程中進(jìn)行技術(shù)整合、系統(tǒng)總裝、調(diào)試、測試和運(yùn)行管理。

表1 系統(tǒng)功能菜單

3）系統(tǒng)數(shù)據(jù)基礎(chǔ)服務(wù)層（即支持系統(tǒng)運(yùn)行的數(shù)據(jù)環(huán)境），提供支持系統(tǒng)運(yùn)行環(huán)境包括硬件、軟件環(huán)境，數(shù)據(jù)基礎(chǔ)服務(wù)層是軟件環(huán)境組成部分，有水情監(jiān)測、水雨情、模型參數(shù)及成果、系統(tǒng)管理、地理信息空間等數(shù)據(jù)庫，從系統(tǒng)功能層次及其整體結(jié)構(gòu)上看，系統(tǒng)是一個采用多庫結(jié)構(gòu)面向水情監(jiān)測分析的系統(tǒng)。

2.2 GIS 應(yīng)用

水資源調(diào)度系統(tǒng)需要很多基礎(chǔ)數(shù)據(jù)的支持，如水文、河道斷面、防洪工程、區(qū)域數(shù)字地面高程、社會經(jīng)濟(jì)等數(shù)據(jù)。這些數(shù)據(jù)大多與現(xiàn)實(shí)中的某一特定的位置或區(qū)域（點(diǎn)、線或面）相聯(lián)系，而且有些數(shù)據(jù)還具有空間屬性。GIS 是一種綜合分析和處理數(shù)據(jù)空間分布的新技術(shù)，能夠使空間數(shù)據(jù)的處理維護(hù)更為迅速高效。采用嵌入 GIS 控件 mapobject，使 GIS 與模型計(jì)算、數(shù)據(jù)庫等無逢連接。以 GIS 為交互界面輕松實(shí)現(xiàn)圖形導(dǎo)航、放大、縮小、漫游、標(biāo)注、圖層控制等基本功能。流域水資源調(diào)度系統(tǒng)中 GIS 的主要專題圖層如表2所示。

2.3 分布式水情預(yù)報

準(zhǔn)確的水情預(yù)報是水庫合理調(diào)度的前提。水情預(yù)報模型采用基于 GIS 分布式模型。分布式水情預(yù)報模型是指模型結(jié)構(gòu)中各水文要素具有空間的分布特性和隨時間變化的動態(tài)特性，其特點(diǎn)是具有物理基礎(chǔ)，能描述水情的時空變化過程；能及時地模擬出人類活動或農(nóng)作物生長變化對水情的影響[5]。

表2 流域水資源調(diào)度系統(tǒng)中主要專題圖層

分布式水情預(yù)報模型的主要思路是：根據(jù) GIS的行政區(qū)、地形、土壤分布、土地利用等圖層將流域劃分成若干子區(qū)域，對每個子區(qū)域分別輸入不同的降雨、蒸發(fā)，根據(jù)選用的產(chǎn)匯流模型（如新安江模型）計(jì)算徑流過程[6]?；?GIS 的分布式水情預(yù)報模型框架如圖1所示。

2.4 水庫調(diào)度計(jì)算

1）調(diào)洪計(jì)算基本方程。水庫洪水調(diào)度通常根據(jù)預(yù)報入庫洪水過程和水庫運(yùn)行規(guī)則推求下泄流量過程和庫水位過程。水庫調(diào)洪是在水量平衡和動力平衡的支配下進(jìn)行的。

圖1 基于 GIS 的分布式模型框架

在 Δt時段內(nèi)，入庫與出庫水量之差，應(yīng)等于該時段內(nèi)水庫蓄水量的變化。水庫水量平衡方程為：

式中：Q1、Q2分別為時段始末的入庫流量；q1、q2分別為時段始末的出庫流量；V1、V2分別為時段始末的水庫蓄水量；Δt為計(jì)算時段長。式 （1）中Q1、Q2均為已知，根據(jù)起調(diào)條件，調(diào)洪起始水位也為已知，即q1、V1為已知，Δt可選定，未知數(shù)有q2和V2，因此方程不能獨(dú)立求解，還需建立第2個方程。

對于某一水庫，當(dāng)泄洪建筑物的型式、尺寸一定時，泄流能力q取決于泄洪設(shè)施的水頭或水庫蓄水量V，故水庫泄流方程可寫為：

或q-V曲線。從起調(diào)條件開始，逐時段連續(xù)求解上述2方程，即可完成洪水調(diào)度計(jì)算。從不同的關(guān)注點(diǎn)出發(fā)，即有不同的調(diào)度模型。

2）泄量控制調(diào)度模型。泄量控制調(diào)度模型是固定泄量的調(diào)洪計(jì)算，約束條件是泄洪設(shè)備的泄洪能力。計(jì)算規(guī)則是根據(jù)輸入的理想泄量，應(yīng)用水量平衡式計(jì)算相應(yīng)的庫水位過程。在該模型中水庫的最高水位是自由變量，其目的為在保證水庫泄量不超過設(shè)定值的前提下，使水庫的最高水位最低。

將式（1）改寫成：

假設(shè)調(diào)度期的泄流過程q2逐時段連續(xù)計(jì)算V2，通過水位-庫容（H-V）關(guān)系推求庫水位過程。

3）水位控制調(diào)度模型。水位控制調(diào)度模型關(guān)注水量平衡中的Vt因子，將控制水位作為控制約束條件，其目標(biāo)是在保證庫水位控制的前提下，使水庫的下泄量最小。模型計(jì)算的控制條件包括：汛限和調(diào)度期最高控制水位。模型求解準(zhǔn)則為：a）水庫剛漲水時，庫水位未達(dá)到汛限水位，出庫流量等于當(dāng)前出庫流量；b）隨著入庫流量不段增大庫水位也不斷增高，當(dāng)庫水位達(dá)到汛限水位時，下泄流量等于入庫流量；c）當(dāng)入庫流量進(jìn)一步增加達(dá)到并超過調(diào)度期最高控制水位時按泄洪能力泄流；d）當(dāng)庫水位回落到控制水位后，來多少泄多少，即泄洪流量等于入庫流量。

3 系統(tǒng)應(yīng)用

3.1 流域水庫概況

霍童溪位于福建省東北部，上游主流棠口溪，發(fā)源于政和、屏南、周寧3縣交界的鷲峰山脈東麓，自西北折向東南流經(jīng)屏南縣的嶺頭、棠口等地，沿途有金造溪、后垅溪等支流匯入，棠口溪和后垅溪于金鐘渡匯合后稱霍童溪干流，再向東流經(jīng)寧德市區(qū)及、霍童、九都、八都等鄉(xiāng)鎮(zhèn)后注入三沙灣?；敉骱拥篱L126km，流域面積2244km2，平均坡降 6.2‰[6]。

流域內(nèi)主要有后壟一二級、溪尾、洪口等大中型水庫，設(shè)有禮門、嶺下等15個遙測雨量站。水情預(yù)報調(diào)度信息來自于流域內(nèi)遙測雨水情站網(wǎng)?；敉饔虼笾行退畮旒斑b測雨水情站網(wǎng)如表3所示。

表3 霍童溪流域水情遙測站網(wǎng)表

3.2 預(yù)報調(diào)度實(shí)例

2009年8月7～10日流域內(nèi)有1次降雨過程累積降雨191mm，流域時段平均最大降雨11mm/h，出現(xiàn)在 2009-08-09T10:00。洪口水庫預(yù)報作業(yè)時（2009-08-09T10:00）水位 158.40m，依據(jù)預(yù)報結(jié)果 2009-08-10T06:00洪口水庫水位將達(dá) 162.77m 超汛限水位（162.70m），并仍有上漲趨勢，按調(diào)度規(guī)則需棄水泄洪。調(diào)整調(diào)度方案，提前加大發(fā)電流量至253m3/s，進(jìn)騰庫發(fā)電，庫水在計(jì)算時段（24h）內(nèi)漲至162.09m 未達(dá)汛限水位。這樣在計(jì)算的24h 內(nèi)可少棄水1093萬 m3，多發(fā)電2400MW。

4 結(jié)語

目前，我國的水庫調(diào)度多以防洪為重，未考慮洪水的資源性，水庫調(diào)度僅關(guān)注汛期限制水位，在汛期，一旦發(fā)生水庫超過汛期限制水位即要求水庫泄流，降至汛期限制水位。這種粗獷式的調(diào)度往往造成洪水后干旱缺水的局面，嚴(yán)重影響水庫流域內(nèi)生產(chǎn)、生活及生態(tài)[7]。上述計(jì)算實(shí)例表明：應(yīng)用流域水庫群水資源調(diào)度系統(tǒng)可有效提高水資源利用率。

建國以來，我國建成的大中小型水庫總庫容約為450多億 m3。其中許多水庫兼有在防洪、灌溉、發(fā)電、航運(yùn)、養(yǎng)殖、城市及工業(yè)供水等多種功能。對于綜合利用的水庫，不同的用水目標(biāo)對水庫蓄水、泄水的水量和時間要求常常會發(fā)生沖突。例如，發(fā)電希望水庫泄流量越穩(wěn)定越好，在農(nóng)作物生長期內(nèi)的農(nóng)業(yè)灌溉需要較多的水量，防洪則需要在預(yù)測到洪水到來時提前降低水庫的水位。因此，推廣使用流域水庫群水資源調(diào)度系統(tǒng)，可獲得兼顧各種需求的水庫調(diào)度方案，前景廣闊。

[1] 紀(jì)志國，孫慶艷.我國水資源可持續(xù)利用措施研究[J].科技資訊，2010（7）: 137.

[2] 郭生練.水庫調(diào)度綜合自動化系統(tǒng)[M].武漢：武漢水利電力大學(xué)出版社，2000.

[3] 王國利，彭勇，何斌，等.基于 B/S+C/S 模式的防洪調(diào)度決策支持系統(tǒng)研究及其應(yīng)用[J].大連理工大學(xué)學(xué)報，2010,3（2）: 258-263.

[4] 朱新軍，王中根，夏軍，等.基于分布式模擬的流域水平衡分析研究[J].地理科學(xué)進(jìn)展，2008,28（4）: 23-27.

[5] 孫西歡，張柏治，王志璋.汾河上游流域分布式水文模型的構(gòu)建[J].水土保持通報，2008,28（3）: 89-92.

[6] 鄭華鋒，寧德市水電開發(fā)[M].北京：水利水電出版社，2008.

[7] 陳鳴，吳永祥，陸衛(wèi)鮮，等.InfoWorks RS、FloodWorks軟件及應(yīng)用[J].水利水運(yùn)工程學(xué)報，2008,12（4）: 19-24.