郭亞萍，李 丹，曹 濱，張 帥，林洪孝，王 剛(.山東農(nóng)業(yè)大學水利土木工程學院，山東 泰安 708；.泰安市大汶河閘壩管理所，山東 泰安 7000)

水系連通對于河流水質(zhì)凈化、防洪排澇以及增加區(qū)域調(diào)蓄和水生態(tài)環(huán)境承載能力等均具有至關(guān)重要的作用。水系連通性的影響因素較多，隨著城市化進程不斷加劇，人類高強度的活動使流域內(nèi)河流出現(xiàn)不同程度的河床淤積、河道堵塞等問題，直接危害了河網(wǎng)水系結(jié)構(gòu)，影響了河流連通性作用的發(fā)揮[1]。孟慧芳的研究中也表明城市化進程加劇導致河網(wǎng)水系中尤其是低等級河道的減少,可降低河流的連通性，此外相關(guān)研究[2]中還表明雨量特征、土地利用狀況以及流域地形等均對水系連通性具有一定的影響。由于在河流的綜合開發(fā)治理過程中，修建水庫、閘壩是普遍采用且較為有效工程措施，人為調(diào)節(jié)了水流時間，改變了水量空間分布，而相關(guān)研究[3-5]也表明對水庫、閘壩水量進行優(yōu)化調(diào)度可以改善流域內(nèi)河流的健康狀況，因此探究水量調(diào)度對水系連通性的影響具有較強的實際意義。

本研究根據(jù)泗河流域水利工程規(guī)劃建設(shè)項目中存在的水庫、閘壩較多，為滿足城市供水、農(nóng)業(yè)灌溉、生態(tài)補水等要求，需對區(qū)域內(nèi)水量進行調(diào)度的實際問題，擬探討以上下游水資源損失最小、最大化攔蓄水量為目標，保證河流水流暢通的前提下，研究對庫區(qū)及閘壩范圍內(nèi)水量進行優(yōu)化調(diào)度時，能否對泗河流域內(nèi)河流水系連通性起到一定的積極作用。

1 研究區(qū)概況

泗河流域發(fā)源于山東省新泰市太平頂，流經(jīng)新泰、泗水、曲阜、兗州、鄒城、高新區(qū)、太白湖新區(qū)、微山，入南陽湖，是南四湖湖東最大入湖河流。流域面積2 357 km2，全長159 km。流域內(nèi)現(xiàn)有水利工程(圖1)中大型水庫1座，中型水庫4座，攔河閘壩25座(現(xiàn)有閘壩11座，規(guī)劃新建閘壩14座)。目前，泗河流域已修建了一系列以蓄、引、防洪等為目的水利工程，雖然發(fā)揮了巨大的社會和經(jīng)濟效益作用，但流域內(nèi)河流遭受人工干擾嚴重，水質(zhì)惡化，河流生態(tài)環(huán)境遭到破壞，降低了流域內(nèi)河流水系的連通性，嚴重阻礙了泗河流域的可持續(xù)發(fā)展。

圖1 泗河流域攔河閘壩分布圖Fig.1 Distribution of Brake and dam on Sihe river basin

2 研究方法

2.1 水量優(yōu)化調(diào)度方案

步驟一：基于泗河流域內(nèi)主干流5庫25閘的分布特點及斷面狀況等因素，設(shè)置目標函數(shù)，以期實現(xiàn)流域內(nèi)梯級閘壩(水庫)攔蓄水量最大化。目標函數(shù)為：

(1)

式中：g(i，t)表示第i個水利樞紐在t時段的攔蓄水量。

步驟二：在滿足泗河流域?qū)嶋H用水量需求的前提下，通過水量平衡約束方程對流域內(nèi)攔蓄水量進行調(diào)度概算。

水量平衡方程：

Vi,t+1=Vi,t+Ii,t+1-(WUi,t+1+WEi,t+1+WAi,t+1+

Wi其他,t+1)-Wi損,t+1-OUTi棄,t+1

(2)

式中：Vi,t+1，Vi,t分別為第i個水利樞紐在t時段初、末的蓄水容積，m3；Ii,t+1為t時段i水利樞紐的來水量，m3；WUi,t+1為t時段i水利樞紐的調(diào)出水量，m3；WEi,t+1為t時段i水利樞紐的生態(tài)補水量，m3；WAi,t+1為t時段i水利樞紐的農(nóng)業(yè)灌溉用水量，m3；Wi其他,t+1為t時段i水利樞紐的其他用水量，m3；Wi損,t+1為t時段i水利樞紐的損失水量，m3；OUTi棄,t+1為t時段i水利樞紐的棄水量，m3。

在進行優(yōu)化調(diào)度計算過程中需同時滿足以下約束條件：變量非負、狀態(tài)約束、水量連續(xù)。其中狀態(tài)約束是指庫壩時段末庫容應(yīng)大于或者等于死庫容，并且小于或者等于興利庫容，以滿足徑流調(diào)節(jié)、水庫供水量的要求；水量連續(xù)主要是保證河流水量流動的連續(xù)性，并將上一級水庫(閘壩)的生態(tài)補水量及棄水量作為下一級閘壩的來水量。

2.2 連通性指標的選取

目前，對于水系連通性的內(nèi)涵尚沒有系統(tǒng)規(guī)范的定義，長江水利委員會編寫的《維護健康長江，促進人水和諧研究報告》中[6]，將水系連通性定義為：河道干支流、湖泊及其他濕地等水系的連通情況，反應(yīng)水流的連通性和水系的連通狀況。在此基礎(chǔ)上，張歐陽[7]等人認為水系連通性就是要滿足兩個基本因素：一是要有能滿足一定需求的保持水流運動的水量；二是要有水流的連接通道。本研究在借鑒相關(guān)研究的基礎(chǔ)上，按照河湖水系的功能和作用選取水資源開發(fā)利用率和本地水資源供需比來表示資源調(diào)配型連通性，選取水質(zhì)指標和水質(zhì)達標率來表示生態(tài)修復(fù)型連通性，選取特大洪水來水量與蓄水能力之比來表示災(zāi)害防御型連通性[8,9]，并根據(jù)山東省水利廳《山東省水資源保護規(guī)劃技術(shù)大綱》選取縱向連通性指標和橫向連通性指標作為評價泗河流域內(nèi)河流結(jié)構(gòu)連通性影響的依據(jù)。

(1)功能連通性指標。

①資源調(diào)配型指標：水資源開發(fā)利用率，反映研究區(qū)水資源開發(fā)利用程度，按照國際標準，若水資源開發(fā)利用率大于60%，則該地區(qū)的水資源開發(fā)利用程度較高，可進一步開發(fā)利用的潛力較小。本地水資源供需比，反映本地缺水程度。

(3)

(4)

②生態(tài)修復(fù)型指標：采用EC(電導率)、五日生化需氧量、化學需氧量和高錳酸鉀指數(shù)來反映研究區(qū)水質(zhì)狀況。采用水質(zhì)達標率來反映水功能區(qū)水質(zhì)的總體狀況。

(5)

③災(zāi)害防御型指標：特大洪水來水量與蓄水能力之比，反映區(qū)域內(nèi)水庫、分滯洪區(qū)等水利工程泄洪能力大小。指標越小，則抵御洪水的能力越強。

(6)

(2)結(jié)構(gòu)連通性指標?？v向連通性是指在河流系統(tǒng)內(nèi)生態(tài)元素在空間結(jié)構(gòu)上的縱向聯(lián)系，橫向連通性指具有連通性的水面面積、長度占統(tǒng)計的水面總數(shù)之比。計算公式如下：

W=N/L

(7)

式中：W是指縱向連通性系數(shù)；N指河流的斷點(或節(jié)點)等障礙物數(shù)量(如閘、壩等，已有過魚設(shè)施的閘壩不在統(tǒng)計范圍之列)；L指連續(xù)河流的有效長度，縱向連通性系數(shù)越大，河流連通性越差。

C=A2/A1

(8)

式中：C指橫向連通性系數(shù)；A1指流域內(nèi)河流水面總面積；A2指流域內(nèi)連通性河流的水面有效面積，橫向連通性指標越大，河流連通性狀況越好。

3 結(jié)果分析與討論

3.1 攔蓄水狀況分析

根據(jù)方案，首先對各閘壩水庫通過截流雨洪水量實現(xiàn)攔蓄水量最大化，然后對水量進行優(yōu)化調(diào)度。結(jié)果表明，通過水量調(diào)度后蓄水情況(多年平均值)可知，豐水年各梯級閘壩年內(nèi)完全蓄滿的月份有7.3個月，占60.6%；部分蓄滿為3.9個月，占32.2%，表明豐水年各級攔河閘壩蓄水狀態(tài)均達到最佳。平水年內(nèi)各梯級閘壩蓄水表現(xiàn)為年內(nèi)完全蓄滿的月份為5.5個月，占45.8%；部分蓄滿為5.6個月，占46.6%，表明平水年各級攔河閘壩能夠達到良好的蓄水狀態(tài)?？菟旮魈菁夐l壩年內(nèi)完全蓄滿的月份為3.0個月，占25%；部分蓄滿為9.0個月，占75%，表明枯水年各級攔河閘壩蓄水狀況差；通過多年平均蓄水狀況可知，總體年內(nèi)完全蓄滿的月份為5.3個月，占44.5%；部分蓄滿為5.3個月，占44.5%，各閘壩蓄水狀況不一，黃陰集閘蓄水狀態(tài)最佳，完全蓄滿的月份占8.3個月，書院壩作為非汛期外調(diào)水功用壩，完全蓄滿的月份為1.8個月；金口壩承擔城市供水任務(wù)，年內(nèi)完全蓄滿月份僅為0.4個月，但部分蓄滿月份為9.9個月，占82.5%，表現(xiàn)為較好的蓄水狀態(tài)?？傮w上，最大化攔蓄水量增加了流域內(nèi)可用水資源量，泗河流域50%總可供水量增加了4 270 萬m3，75%總可供水量增加了3 806 萬m3，95%總可供水量增加了3 109 萬m3，泗河流域內(nèi)可用水資源總量明顯增加(圖2)。

3.2 水量優(yōu)化調(diào)度對水系連通性的影響分析

(1)對功能連通性的影響分析。通過計算分析各功能性指標，得到研究區(qū)水量優(yōu)化調(diào)度前后的功能性特性(表1)。

水資源開發(fā)利用率降低，本地水資源供需比提高。水資源總量由規(guī)劃調(diào)度前的15 310 m3，增加到規(guī)劃調(diào)度后的20 680m3。隨著水資源總量的增加，水資源開發(fā)利用率也由原來的34.7%，減少到25.7%，水資源開發(fā)利用潛力較大?，F(xiàn)狀水平年75%頻率的供需比為0.857，優(yōu)化調(diào)度后，75%頻率的供需比變?yōu)?.903，水資源短缺狀況有所緩解。

圖2 泗河流域水利工程供水量成果Fig.2 Water supply of water conservancy projects on Sihe river basin

表1 水量優(yōu)化調(diào)度前后功能連通性特征變化Tab.1 Change of functional connectivity features before and after the water optimization of scheduling

水質(zhì)狀況有所改善。在我國，目前對于河流保護的主要工作是水質(zhì)的恢復(fù)[10]。基于聯(lián)合調(diào)度改善水質(zhì)狀況在我國大流域內(nèi)也有成功的應(yīng)用[11,12]。在泗河流域內(nèi)進行水量優(yōu)化調(diào)度可增加流域內(nèi)河流的連通性，有助于增強水流的自凈能力與納污能力[13]，達到改善水質(zhì)的效果。同時，通過對流域內(nèi)相關(guān)庫區(qū)閘壩進行水質(zhì)指標監(jiān)測表明(圖4)，在豐水期各水質(zhì)指標明顯低于枯水期，這可能是由于河道內(nèi)河水水量與流速增加[14]，在污染物數(shù)量穩(wěn)定的前提下，由于水體的稀釋作用，增加了水體的環(huán)境容量與自凈能力[15]，降低水質(zhì)的污染狀況，郭巍[16]通過分析COD、NH3-N等與水量的相關(guān)關(guān)系得出，隨著干流流量增加，河段中COD、NH3-N等污染物濃度會有一定程度的下降，與本研究結(jié)果相符。通過對水量調(diào)度前后泗河流域總可供水量進行比對表明，泗河流域50%總可供水量增加了4 270 萬m3，75%總可供水量增加了3 806萬m3，95%總可供水量增加了3 109 萬m3，泗河流域內(nèi)水用水資源總量明顯增加。劉韜[17]在綜述總結(jié)認為，資源具有質(zhì)、量二重性，二者相互影響；徐斌等[15]研究表明，地表水資源量對水質(zhì)具有決定性影響。因此可以認為，水量調(diào)度前后，泗河流域河流水質(zhì)將得到改善，水質(zhì)達標率也由優(yōu)化調(diào)度前的90%預(yù)計達到規(guī)劃調(diào)度后的100%。

圖3 泗河流域水質(zhì)監(jiān)測成果Fig.3 Results of water quality monitoring on Sihe river basin

抵御洪水的能力大大增加。隨著泗河流域水利樞紐的的規(guī)劃與建設(shè)以及水量的優(yōu)化調(diào)度，水利工程最大蓄水量由原來的45 301 萬m3增加到49 447 萬m3，所以特大洪水來水量與蓄水能力之比也由原來的8.39減少到了7.68，說明泗河流域?qū)Φ钟樗哪芰τ兴鰪姟?/p>

(2)對結(jié)構(gòu)連通性的影響分析。本研究在借鑒相關(guān)研究的基礎(chǔ)上，從河流縱向連通性與橫向連通性的角度，對泗河流域內(nèi)通過水量優(yōu)化調(diào)度后的河流連通性進行分析，結(jié)果表明：優(yōu)化調(diào)度前泗河河流的斷點(或節(jié)點)等障礙物數(shù)量為963座(如閘、壩等，已有過魚設(shè)施的閘壩不在統(tǒng)計范圍之列)，考慮到水量不足導致的個別小型支流斷流，經(jīng)計算統(tǒng)計不斷流河流的長度384.6 km，由公式(7)計算得縱向連通性系數(shù)為2.50；優(yōu)化調(diào)度調(diào)整后前泗河河流的斷點或節(jié)點等障礙物數(shù)量為977座(規(guī)劃新增14座橡膠壩)，經(jīng)計算統(tǒng)計連續(xù)河流的長度387.3 km，泗河縱向連通性系數(shù)為2.52，而縱向連通性指標越大，連通性越差，由上述分析可以看出，新增閘壩降低了泗河流域內(nèi)河流連通性，這與羅賢[18]、張歐陽[13]等的研究結(jié)果一致。同時，由濟寧水文局調(diào)查可知，優(yōu)化調(diào)度前泗河具有連通性的水面面積為918 km2，經(jīng)統(tǒng)計的水面面積即泗河流域面積2 357 km2，由公式(8)計算得橫向連通性指標為38.94%，水量優(yōu)化調(diào)度后，保證流域內(nèi)河流連續(xù)，使泗河具有連通性的水面面積增加，經(jīng)濟寧市水文局統(tǒng)計約為1 285 km2，優(yōu)化調(diào)度后橫向連通性指標約為54.52%，河流橫向連通性系數(shù)越大，則河流連通性狀況越好，因此，水量優(yōu)化調(diào)度在一定程度上改善了由于修建水庫、閘壩等對河流連通性造成的不利影響，使河流連通性提高(表2)。

表2 水量優(yōu)化調(diào)度前后結(jié)構(gòu)連通性特征變化Tab.2 Change of structure connectivity characteristics before and after the water optimization of scheduling

優(yōu)化調(diào)度后相對于優(yōu)化調(diào)度前，泗河連通性整體上有所改善，但相對于自然狀況，人為修建的水庫、閘壩等水工建筑物仍降低了河流本身的連通性。作為平原河網(wǎng)，泗河流域內(nèi)河流連通性的降低主要是由于城市化過程中高強度的人類干擾，使廊道的連接率明顯下降[19,20]。河流上修建水庫與閘壩在一定程度上阻斷了水流的連續(xù)性，降低了水體的自凈能力，使同一斷面下的流量減少[13]，然而河網(wǎng)連通狀況還受到閘門與泵站的調(diào)度影響[19]，因此，可以通過對改進已建閘壩的調(diào)度運行方式作為非工程措施改進水系流通性[21]。

4 結(jié) 語

(1)水量優(yōu)化調(diào)度可以改善水系的功能連通性。水量優(yōu)化調(diào)度可以使得研究區(qū)水資源得到合理調(diào)配，改善研究區(qū)河流水質(zhì)，增強研究區(qū)的災(zāi)害防御能力。從而使得研究區(qū)的功能連通性得到改善。

(2)水量優(yōu)化調(diào)度后，可以彌補由于水利工程建設(shè)造成的水系格局連通性變差的狀況，改善河流連通性，為流域水利工程建設(shè)規(guī)劃和流域生態(tài)恢復(fù)提供依據(jù)。

□

[1] 陳云霞, 許有鵬, 付維軍.浙東沿海城鎮(zhèn)化對河網(wǎng)水系的影響[J]. 水科學進展, 2007,18(1):68-73.

[2] Meerkerk A L, Wesemael B V, Bellin N. Application of connectivity theory to model the impact of terrace failure on runoff in semi-arid catchment[J]. Hydrological Processes, 2009,23(19):2 792-2 803.

[3] 劉玉年,施 勇,程緒水,等.淮河中游水量水質(zhì)聯(lián)合調(diào)度模型研究[J]. 水科學進展, 2009,20(2):177-183.

[4] 吳浩云. 大型平原河網(wǎng)地區(qū)水量水質(zhì)耦合模擬及聯(lián)合調(diào)度研究[D]. 南京：河海大學, 2006.

[5] 張松達, 溫進化, 夏夢河. 水庫群以及水庫與河網(wǎng)聯(lián)合調(diào)度實踐與探討[C]∥ 中國原水論壇專輯，2010.

[6] 蔡其華. 維護健康長江 促進人水和諧[J]. 水利建設(shè)與管理, 2005,25(8):55-57.

[7] 張歐陽, 熊 文, 丁洪亮. 長江流域水系連通特征及其影響因素分析[J]. 人民長江, 2010,41(1):1-5.

[8] 崔國韜, 左其亭, 李宗禮,等. 河湖水系連通功能及適應(yīng)性分析[J]. 水電能源科學, 2012,30(2).

[9] 竇 明, 崔國韜, 左其亭,等. 河湖水系連通的特征分析[J]. 中國水利, 2011,(16):17-19.

[10] 董哲仁. 河流健康評估的原則和方法[J]. 中國水利, 2005,(10):17-19.

[11] 李考真, 任淑梅. 地表水水質(zhì)水量聯(lián)合調(diào)度研究——以徒駭河聊城段為例[J]. 聊城師院學報:自然科學版, 1999,(2):72-76.

[12] 程緒水, 沈哲松, 賈 利. 2001-2002年度淮河沙潁河水污染聯(lián)防調(diào)度與經(jīng)驗[J]. 治淮, 2002,(12):3-4.

[13] 張歐陽, 卜惠峰, 王翠平,等. 長江流域水系連通性對河流健康的影響[J]. 人民長江, 2010, 41(2):1-5.

[14] 王錦旗, 王國祥, 鄭建偉. 不同流速下溢流堰對河道水質(zhì)的影響[J]. 人民黃河, 2014,(11):62-64.

[15] 徐 斌, 劉克巖, 時曉飛,等. 地表水資源水質(zhì)水量結(jié)合評價方法研究[J]. 水利規(guī)劃與設(shè)計, 2003,(2):25-30.

[16] 郭 巍. 水質(zhì)水量結(jié)合評價渭河干流(陜西段)水資源變化[J]. 水資源與水工程學報, 2011,22(5):115-120.

[17] 劉 韜. 水資源水質(zhì)水量聯(lián)合評價綜述[J]. 環(huán)境科學導刊, 2012,31(2):73-77.

[18] 羅 賢, 許有鵬, 徐光來,等. 水利工程對河網(wǎng)連通性的影響研究——以太湖西苕溪流域為例[J]. 水利水電技術(shù), 2012,43(9):12-15.

[19] 韓龍飛, 許有鵬, 邵玉龍,等. 城市化對水系結(jié)構(gòu)及其連通性的影響——以秦淮河中、下游為例[J]. 湖泊科學, 2013,(3):335-341.

[20] 孟慧芳, 許有鵬, 徐光來,等. 平原河網(wǎng)區(qū)河流連通性評價研究[J]. 長江流域資源與環(huán)境, 2014,23(5).

[21] 董哲仁, 王宏濤, 趙進勇,等. 恢復(fù)河湖水系連通性生態(tài)調(diào)查與規(guī)劃方法[J]. 水利水電技術(shù), 2013,44(11):8-13.