王興菊,逄曉騰,張衍福,趙華青,劉 凱
(1.山東大學土建與水利學院,濟南250061;2.山東省水利綜合事業(yè)服務中心,濟南250014;3.北京水務咨詢有限公司,北京100038)
河網(wǎng)調(diào)蓄是將河道水系作為天然的雨水調(diào)蓄池,利用調(diào)蓄容積對雨水進行蓄滲和緩釋[1]的一種相對經(jīng)濟而有效的調(diào)蓄措施[2],是影響其水資源承載力以及洪澇災害空間分布的重要因素[3]。汛期河網(wǎng)調(diào)蓄能力的發(fā)揮可以削減洪峰,遲滯洪水過程,緩解城市排水排澇壓力;非汛期,河網(wǎng)蓄水則補給地下水,調(diào)節(jié)河川徑流[2],同時能滿足城市生態(tài)需水及城市景觀的需要[3]。尤其在平原河網(wǎng)地區(qū),河流密集,因而其調(diào)蓄能力對流域水系的生態(tài)補水及防洪排澇具有更加重要的意義。
河網(wǎng)調(diào)蓄能力是一個與河網(wǎng)特征水位相關(guān)的閾值區(qū)間[4],楊柳等[3]采用常水位及警戒水位為河網(wǎng)調(diào)蓄特征水位,對1960s-2010s 蘇州水網(wǎng)區(qū)的河網(wǎng)槽蓄容量、可調(diào)蓄容量、單位面積槽蓄容量及單位面積可調(diào)蓄容量進行分析計算,得到研究區(qū)河網(wǎng)調(diào)蓄能力閾值區(qū)間,為分析研究區(qū)河網(wǎng)對不同暴雨洪水的調(diào)蓄功能承受能力提供依據(jù)。
而目前國內(nèi)外對河網(wǎng)調(diào)蓄能力的研究多局限于計算某一水位下的河網(wǎng)調(diào)蓄能力。王臘春等[5]采用河道非恒定流方法,模擬計算了太湖流域1991年最高水位下的槽蓄容量。王躍峰等[6]以旬平均水位為例,研究了1960-2010年太湖流域武澄錫虞地區(qū)的水系演變與河網(wǎng)調(diào)蓄能力變化,探討了河網(wǎng)調(diào)蓄與河網(wǎng)結(jié)構(gòu)參數(shù)的關(guān)系。袁雯等[7]分析計算了常水位情景下上海市河網(wǎng)的單位面積槽蓄容量及可調(diào)蓄總量,并探究了河流結(jié)構(gòu)對河網(wǎng)調(diào)蓄能力的可能影響。李景保等[8]選取河網(wǎng)槽蓄容量、可調(diào)蓄容量、單位面積槽蓄容量、單位面積可調(diào)蓄容量,對長江荊南三口水系常水位下的調(diào)蓄能力進行分析計算,并剖析了水系結(jié)構(gòu)與調(diào)蓄能力的關(guān)系。
河流水資源結(jié)構(gòu)分解[9]是將宏觀河流水資源規(guī)劃利用和微觀河流水資源結(jié)構(gòu)分析相結(jié)合的一種水量分解方法,是使河流實現(xiàn)定量分解、微觀調(diào)控的一種重要手段。李文義等[9]從生態(tài)、利用、防洪的角度出發(fā),以生態(tài)水位、平灘水位、警戒水位為特征水位,將河流水量劃分為生態(tài)水量、安全水量、風險水量和災害水量四部分。根據(jù)各類水量的功能和特點,確定了不同特征水量的量化方法,并對嫩江下游的江橋和大賚斷面河流水資源結(jié)構(gòu)進行分析計算。許士國等[10]綜合考慮生態(tài)保護、興利運用和防洪減災等因素,將河流水資源劃分為生態(tài)水量、安全水量、風險水量和災害水量,并對松花江哈爾濱水文斷面的水資源結(jié)構(gòu)進行劃分,對其各類水量分配規(guī)律進行分析計算。
桓臺縣地勢平坦,河流縱橫,卻仍存在河道干涸、水資源供需矛盾突出、客水資源利用率低的問題,因此,本文提出了基于河流水資源結(jié)構(gòu)分解的河網(wǎng)調(diào)蓄能力計算方法。通過對桓臺縣河流水資源進行結(jié)構(gòu)分解,并以此為基礎提出桓臺縣河網(wǎng)調(diào)蓄特征水位,對其進行不同特征水位下河網(wǎng)調(diào)蓄能力計算,為桓臺縣河網(wǎng)水資源科學調(diào)控及涵養(yǎng)地下水提供基礎支撐。
桓臺縣位于山東半島中部,地處魯中山區(qū)和魯北平原結(jié)合地帶,全縣總面積509.43 km2。境內(nèi)地形和緩,起伏微小。南部、中部為山前傾斜洪積平原,北部為黃泛平原,地勢南高北低,同時自西南向東北緩傾,比降在1/8 000~1/3 000之間?;概_縣屬暖溫帶季風型大陸性氣候,降水年際、年內(nèi)分配不均。豐枯年降雨量相差懸殊,年內(nèi)降雨量集中在夏秋兩季,春旱秋澇是本地區(qū)的氣候特點。境內(nèi)水系發(fā)達,水網(wǎng)密布,河流均屬于小清河水系。水系包括行洪河道、湖泊(錦秋湖、馬踏湖)、溝渠和排灌兩用渠道。其中,小清河、杏花河、孝婦河、預備河等主要行洪河道自西向東流向,在縣境北部排泄河水,烏河、東豬龍河、西豬龍河、勝利河及澇淄河受地勢影響均為由南向北流向,匯入小清河后流入渤海?;概_縣水系圖如圖1所示。
桓臺縣人均占有可利用水資源量為234 m3/人,低于全省平均水平(334 m3/人),僅占全國人均占有量的10%,屬于嚴重缺水地區(qū)?;概_縣水資源利用主要以地下水資源為主,地下水資源利用占到全縣用水量的60%以上,客水和非常規(guī)水資源沒有得到充分利用。2012年引黃最大量約0.66億m3,利用率僅占引黃指標的82.02%,其余年份引黃利用率更低,究其原因是桓臺縣對水資源沒有足夠的調(diào)蓄能力?;概_縣地下水源包括淺層地下水和深層地下水,淺層地下水主要用于農(nóng)業(yè)灌溉,深層地下水因水質(zhì)較好,長期以來一直作為城鄉(xiāng)生產(chǎn)生活的主要供水水源。近年來,隨著需水量的增加,地下水資源開采量逐年增加,導致水源地及周邊地下水位持續(xù)下降,形成大面積地下水漏斗區(qū)。根據(jù)地下水資源現(xiàn)狀監(jiān)測,淺層地下水年均下降速度為0.4 m/a,深層地下水下降速度最大已達4.23 m/a。地下水位持續(xù)下降,導致地下水環(huán)境不斷惡化,產(chǎn)生地面裂縫、沉降等環(huán)境地質(zhì)問題,并引發(fā)樓房裂縫甚至倒塌等次生災害,同時也對周邊農(nóng)村生活及生產(chǎn)用水產(chǎn)生不利影響。為了有效扼制這種狀況持續(xù)發(fā)展,2015年山東省利用亞行(Asian Development Bank—ADB)貸款地下水漏斗區(qū)域綜合治理示范項目對桓臺縣地下水漏斗區(qū)進行綜合整治。通過項目的實施,形成了由28條主要行洪河道組成的“三橫四縱”的區(qū)域水網(wǎng),水系的連通有利于水文調(diào)蓄、生態(tài)保護等河流功能的發(fā)揮[11],為引用客水資源進行生態(tài)補水和利用雨洪水資源涵養(yǎng)地下水提供了良好的條件。
然而,桓臺縣區(qū)域河網(wǎng)尚缺少生態(tài)補水及合理調(diào)蓄雨洪水資源的依據(jù),非汛期河道干涸,汛末由于懼怕暴雨洪水再次侵襲,出現(xiàn)“想蓄而不敢蓄”最終無水可蓄的尷尬局面。2018年“溫比亞”與2019年“利奇馬”兩次超常規(guī)臺風暴雨,給桓臺縣帶來嚴重的洪澇災害,但暴雨過后河道無水存蓄。
首先,按照行洪等效性原則[12]將研究區(qū)河網(wǎng)概化成28條主要行洪河道;其次,依據(jù)河流水資源結(jié)構(gòu)分解思想,從生態(tài)、經(jīng)濟、防洪的角度出發(fā),分別采用生態(tài)水位、除澇水位、防洪設計水位、堤頂高程為特征水位,將河道分解為生態(tài)水量、安全水量、風險水量、災害水量四種特征水量,建立桓臺縣河流水資源結(jié)構(gòu)分解模型。相應地,上述特征水位即生態(tài)水位、安全水位、風險水位、災害水位。采用Tennant 法[13]計算概化后河道的生態(tài)水量;采用斷面地形法[14],依據(jù)各河道治理工程橫斷面設計圖,利用SOLIDWORKS 建立河道三維立體模型,并計算概化后河道的安全水量、風險水量、災害水量。最后,將生態(tài)水位、安全水位、風險水位作為河網(wǎng)調(diào)蓄特征水位,并選取河網(wǎng)槽蓄容量、可調(diào)蓄容量、單位面積槽蓄容量、單位面積可調(diào)蓄容量為河網(wǎng)調(diào)蓄能力指標,分別將生態(tài)水位、安全水位、風險水位作為桓臺縣生態(tài)調(diào)度、資源調(diào)度、防洪調(diào)度的控制目標水位,計算上述三種調(diào)度情景下4 個河網(wǎng)調(diào)蓄能力指標參數(shù)值,得到桓臺縣河網(wǎng)調(diào)蓄能力閾值區(qū)間。
(1)水系數(shù)據(jù)。基本圖件資料來源于桓臺縣1∶125 000 電子水系圖(350 dpi)、桓臺縣1∶100 000 電子洪水調(diào)度圖(350 dpi)及空間分辨率1.19 m衛(wèi)星遙感影像。
(2)水文、氣象資料。水文、氣象資料來源于《山東省水文圖集》,包括收集的研究區(qū)多年平均徑流量和多年平均蒸發(fā)量數(shù)據(jù)。
(3)河道數(shù)據(jù)。來源于桓臺縣水務局提供的河道治理工程設計報告。
首先,將空間分辨率1.19 m 衛(wèi)星遙感影像圖作為底圖,在ArcGIS10.2的支持下分別對電子水系圖及洪水調(diào)度圖進行地理配準,其中,水系圖及洪水調(diào)度圖精度均高于300 dpi,確保能分辨最細小的河流[15]。隨后,根據(jù)行洪等效性原則[12],結(jié)合桓臺縣自然現(xiàn)狀及資料情況對桓臺縣河網(wǎng)進行概化,增補、移動、刪除部分水系,包括連通性、代表性較差及行洪能力不足的河流,最終將桓臺縣河網(wǎng)概化成28 條主要行洪河道。河網(wǎng)概化并未改變原河道結(jié)構(gòu)特征,因此可采用原有河道特征水位計算各特征蓄水量?;概_縣河網(wǎng)概化圖如圖2所示。
圖2 桓臺縣河網(wǎng)概化圖Fig.2 The figure of the simplified river network in Huantai County
河流水資源結(jié)構(gòu)分解的基本思想是:按照河流水資源承擔的不同功能和作用將其劃分為不同組成部分,為河流水資源合理調(diào)配提供理論依據(jù)。
為方便合理地計算桓臺縣河網(wǎng)調(diào)蓄能力閾值區(qū)間,為其生態(tài)補水、合理存蓄客水及雨洪水資源、防洪減災提供理論依據(jù),本文從河道不同水量的天然屬性出發(fā),以生態(tài)水位、河道設計水位中的除澇水位及防洪設計水位、堤頂高程為特征水位,將河流水量劃分為生態(tài)水量、安全水量、風險水量、災害水量四部分,同時,除澇水位、防洪設計水位、堤頂高程分別也可稱為安全水位、風險水位、災害水位。顯然,安全水量和風險水量即為河網(wǎng)調(diào)蓄水量。河流水資源結(jié)構(gòu)分解模型如圖3所示。
圖3 河流水資源結(jié)構(gòu)分解模型Fig.3 The river water division model
基于河網(wǎng)調(diào)蓄能力的內(nèi)涵,同時考慮桓臺縣河道特征資料的可取性,選取槽蓄容量、可調(diào)蓄容量、單位面積槽蓄容量和單位面積可調(diào)蓄容量4個指標,并進行量化,用其值表征桓臺縣河網(wǎng)調(diào)蓄能力[8]。
(1)槽蓄容量(C)。槽蓄容量是指河流水位在一般情景下,河道所承載的水體的總?cè)萘浚?],其數(shù)字的變化直接反映該研究區(qū)容蓄水資源量的大小。在生態(tài)補水即生態(tài)調(diào)度、合理存蓄客水及雨洪水資源即資源調(diào)度、防洪減災即防洪調(diào)度背景下,槽蓄容量在本文中分別是桓臺縣的河流水位在生態(tài)水位、安全水位、風險水位時河道所承載的水體總?cè)萘?。由于概化后河道為梯形斷面,因此槽蓄容量的計算公式為?/p>
化簡得:
式中:在生態(tài)調(diào)度、資源調(diào)度、防洪調(diào)度中,Ah分別為生態(tài)水位、安全水位、風險水位,m;Dg為河底的相對高程,m;L代表河道長度,m;W代表河流的對應寬度,m;m代表河道斷面邊坡系數(shù),無量綱。
(2)可調(diào)蓄容量(AC)??烧{(diào)蓄容量是指河流在一般情況下,可以連續(xù)地最大限度承載的水體的總?cè)萘浚?],或者說是河道水位由一定水位上升到一定水位時這一段梯形河道所承載的水體的總?cè)萘?,是評價地區(qū)是否易發(fā)生洪澇災害的重要指標[3]。為滿足桓臺縣河網(wǎng)生態(tài)、資源及防洪需求,本文將可調(diào)蓄容量定義為河道水位分別由生態(tài)水位、安全水位到風險水位時這一段梯形河道所承載的水體的總?cè)萘?,其計算公式為?/p>
式中:Cj為河流在風險水位時河流的槽蓄容量;在生態(tài)調(diào)度、資源調(diào)度、防洪調(diào)度中,Ci分別為河流在生態(tài)水位、安全水位、風險水位時河流的槽蓄容量。
(3)單位面積槽蓄容量(SR)。單位面積槽蓄容量是指水位在一般條件下,河網(wǎng)槽蓄容量和區(qū)域水系面積之比。該比值表示研究區(qū)域河網(wǎng)的蓄水能力,比值的大小能直觀地反映出研究區(qū)域河網(wǎng)蓄水能力的相對強弱。其計算公式為:
式中:A為研究區(qū)所對應的水系片區(qū)的面積。
(4)單位面積可調(diào)蓄容量(ASR)。單位面積可調(diào)蓄容量是指河道水位由一定水位上升到一定水位時這一段梯形河道所承載的水體的總?cè)萘颗c所研究區(qū)域水系面積之間的比值。該比值表示研究區(qū)域河網(wǎng)對洪水的調(diào)節(jié)能力。其計算公式為:
基于2.2 所述水文氣象資料、河道數(shù)據(jù),利用Tennant 法及河流水資源結(jié)構(gòu)分解方法得到概化河網(wǎng)的各河道生態(tài)水位、安全水位、風險水位及災害水位4 個特征水位。利用SOLID?WORKS 對相應斷面數(shù)據(jù)進行處理,并采用斷面地形法計算各河道在4個特征水位下的特征水量。
桓臺縣河網(wǎng)四種特征水量計算結(jié)果見表1。
表1 桓臺縣河網(wǎng)特征水量計算結(jié)果 萬m3Tab.1 The calculation results of river network characteristic water quantity in Huantai County
運用式(2)~(5)分別計算該研究區(qū)在生態(tài)調(diào)度、資源調(diào)度、防洪調(diào)度背景下的河網(wǎng)調(diào)蓄能力指標,計算結(jié)果見表2。
表2 生態(tài)調(diào)度、資源調(diào)度、防洪調(diào)度背景下桓臺縣河網(wǎng)調(diào)蓄能力指標計算結(jié)果Tab.2 The calculation results of the storage and flood control capacity indexes of Huantai river network under the background of ecological operation,resource operation and flood control operation
3.2.1 結(jié)果分析
(1)河流水資源結(jié)構(gòu)分解方面:由表1 可知,隨著桓臺縣河網(wǎng)特征水位的抬高,相鄰兩個特征水位之間的特征水量先增加后減少。其中,維持河網(wǎng)基本生態(tài)功能的生態(tài)水量最少,為24.53 萬m3,保證河流安全和經(jīng)濟生活用水的安全水量最多,為1 694.73 萬m3,風險水量及災害水量差值相對不大,約為河道槽蓄量的30%。
(2)河網(wǎng)調(diào)蓄能力方面:由表2 可知,桓臺縣河網(wǎng)單位面積槽蓄容量閾值區(qū)間為[0.05,6.05]萬m3/km2,單位面積可調(diào)蓄容量閾值區(qū)間為[0,6]萬m3/km2,資源調(diào)度背景下河網(wǎng)對洪水資源的調(diào)蓄能力約為生態(tài)調(diào)度背景下的1/2。在防洪調(diào)度背景下河網(wǎng)已無任何調(diào)蓄空間,這是由于此時河網(wǎng)水位已達風險范圍內(nèi)最高限值,若繼續(xù)存蓄水資源將對周邊造成洪澇危害。
(3)河流水資源結(jié)構(gòu)分解與河網(wǎng)調(diào)蓄能力閾值的關(guān)系方面:結(jié)合表1、2可知,不同調(diào)度情景下河網(wǎng)槽蓄容量均為該情景河網(wǎng)控制目標水位以下特征水量之和,河網(wǎng)可調(diào)蓄容量均為該情景河網(wǎng)控制目標水位與風險水位之間特征水量之和。由于安全水量和風險水量差值相對不大、防洪調(diào)度情景下河網(wǎng)控制特征水位即為風險水位,故資源調(diào)度背景下河網(wǎng)對洪水資源的調(diào)蓄能力約為生態(tài)調(diào)度背景下的1/2、防洪調(diào)度背景下河網(wǎng)對洪水資源無調(diào)蓄能力。
3.2.2 討 論
(1)基于河流水資源結(jié)構(gòu)分解的河網(wǎng)調(diào)蓄能力計算方法,相較于現(xiàn)有關(guān)于河網(wǎng)調(diào)蓄能力計算方面的文獻,該方法體現(xiàn)了調(diào)蓄能力應具有的區(qū)間屬性?;概_縣存在地下水超采嚴重、河道干涸、客水資源利用率低及汛期洪澇災害頻發(fā)的問題,水量配置既要考慮通過水量調(diào)配合理調(diào)蓄水資源,實現(xiàn)客水及雨洪水資源高效利用,也要考慮生態(tài)補水和防洪減災的要求。故將其河網(wǎng)按生態(tài)水位、安全水位、風險水位、災害水位進行河流水資源結(jié)構(gòu)劃分,得到滿足相應功能需求的河網(wǎng)調(diào)蓄能力閾值,方法靈活多變,推廣性較強。
(2)基于河流水資源結(jié)構(gòu)分解的河網(wǎng)調(diào)蓄能力與調(diào)蓄量。①結(jié)合河道斷面資料及不同特征水位,桓臺縣河網(wǎng)水資源結(jié)構(gòu)分解后,生態(tài)水量、安全水量、風險水量和災害水量分別為24.53,1 694.73,1 361.48 和1 176.45 萬m3。以上水量分別對應桓臺縣河網(wǎng)的生態(tài)水位、安全水位、風險水位等特征水位,其中生態(tài)水量是從保證河流最低限生態(tài)環(huán)境不被破壞的角度確定的,故量值最低。②考慮河流水資源結(jié)構(gòu)分解與河網(wǎng)調(diào)度的關(guān)系,基于河網(wǎng)調(diào)蓄能力指標體系構(gòu)建,研究不同調(diào)度情景下的河網(wǎng)調(diào)蓄能力?;概_縣河網(wǎng)的槽蓄容量、可調(diào)蓄容量、單位面積槽蓄容量和單位面積可調(diào)蓄容量分別為[24.53,3 080.74]萬m3,[0,3 056.21]萬m3,[0.05,6.05]萬m3/km2和[0,6]萬m3/km2。由于生態(tài)水位、安全水位、風險水位值依次增大,其河網(wǎng)可蓄水空間依次減少,故河網(wǎng)在生態(tài)調(diào)度、資源調(diào)度、防洪調(diào)度下的調(diào)蓄能力依次降低,計算結(jié)果合理、可靠。
(3)基于河流水資源結(jié)構(gòu)分解的不同季節(jié)河網(wǎng)的調(diào)控。當河網(wǎng)水位為安全水位即河網(wǎng)內(nèi)為生態(tài)水量和安全水量時,河網(wǎng)水量能在保證區(qū)域安全的前提下,盡可能地滿足區(qū)域生態(tài)環(huán)境保護及社會經(jīng)濟發(fā)展的需要。當河網(wǎng)水位為風險水位即河網(wǎng)內(nèi)為生態(tài)水量、安全水量、風險水量時,河網(wǎng)能在其風險可控的前提下調(diào)蓄盡可能多的水量,從而提升地下水的涵養(yǎng)效益。因此,非汛期應采取閘壩調(diào)控措施,調(diào)引客水,在安全范圍內(nèi),盡可能地保持較高水位,以實現(xiàn)地下水涵養(yǎng)與生態(tài)環(huán)境保護的效果;而在汛期,防汛部門應在暴雨來臨前及時降水位,增大河道調(diào)蓄容積,提高排澇能力,以避免短歷時強降雨的侵襲。
(4)依據(jù)各河道治理工程橫斷面設計圖,應用SOLID?WORKS 建立的河道三維立體模型還原度較高,相較于其他編程擬合計算方法能更好地反映河道真實的槽蓄容量情況,計算結(jié)果更為準確。
(5)在計算生態(tài)水量時,本文未劃分研究時段,未來的研究還應根據(jù)生態(tài)保護目標提出相應的研究時段劃分方法。如季節(jié)的豐枯、魚類產(chǎn)卵期等。
本文分別采用生態(tài)水位、除澇水位、防洪設計水位、堤頂高程作為生態(tài)水位、安全水位、風險水位、災害水位等四種特征水位,將桓臺縣河流水資源分解成生態(tài)水量、安全水量、風險水量、災害水量等四種特征水量,建立了桓臺縣河流水資源結(jié)構(gòu)分解模型,構(gòu)建了河網(wǎng)調(diào)蓄能力指標體系,研究分析桓臺縣區(qū)域河網(wǎng)調(diào)蓄能力閾值區(qū)間。主要結(jié)論如下:
(1)基于河流水資源結(jié)構(gòu)分解的河網(wǎng)調(diào)蓄能力計算方法,能計算河網(wǎng)滿足不同功能需求的調(diào)蓄能力閾值區(qū)間,充分體現(xiàn)調(diào)蓄能力應具有的區(qū)間屬性,且能根據(jù)各研究對象存在的問題,針對性地對其河流水資源提出不同的劃分方法,得到滿足相應功能需求的河網(wǎng)調(diào)蓄能力閾值,方法靈活多變,推廣性較強。
(2)結(jié)合河道斷面資料及不同特征水位,考慮河流水資源結(jié)構(gòu)分解與河網(wǎng)調(diào)度的關(guān)系,基于河網(wǎng)調(diào)蓄能力指標體系構(gòu)建,研究不同調(diào)度情景下的河網(wǎng)調(diào)蓄能力,為區(qū)域河網(wǎng)的水量調(diào)度提供參考?;概_縣河網(wǎng)水資源結(jié)構(gòu)分解后,生態(tài)水量、安全水量、風險水量和災害水量分別為24.53,1 694.73,1 361.48和1 176.45 萬m3;河網(wǎng)在生態(tài)調(diào)度、資源調(diào)度、防洪調(diào)度下的調(diào)蓄能力分別為[24.53,3 080.74]萬m3,[0,3 056.21]萬m3,[0.05,6.05]萬m3/km2和[0,6]萬m3/km2。
(3)桓臺縣河網(wǎng)在生態(tài)調(diào)度情景下河網(wǎng)可調(diào)蓄容量最大,此時河道內(nèi)的水資源不僅能避免河道干涸,達到生態(tài)補水的效果,且能使客水資源得到較大程度的利用;資源調(diào)度情景下,河道內(nèi)的水資源均為安全水量,在安全水位范圍內(nèi),盡可能地保持較高水位,以實現(xiàn)地下水涵養(yǎng)與生態(tài)環(huán)境保護的效果;防洪調(diào)度情景下,充分發(fā)揮河網(wǎng)水系連通度較高的優(yōu)勢,通過閘壩調(diào)控使各河道盡可能保持在風險水位以下,達到防洪減災的目的。□