趙鳴雁 孔令仲 鄭艷俠 雷曉輝 權(quán)錦

摘要：針對(duì)采用應(yīng)急同步關(guān)閘措施導(dǎo)致的水位雍高問(wèn)題，分析了多閘門(mén)同步關(guān)閘情況下不同關(guān)閘時(shí)間對(duì)水位雍高的影響。通過(guò)建立一維水力學(xué)模擬模型，進(jìn)行了大量的關(guān)閘工況計(jì)算，研究了上、下游節(jié)制閘分別關(guān)閘動(dòng)作對(duì)閘前水位的影響過(guò)程，結(jié)果表明，下游節(jié)制閘的關(guān)閘過(guò)程中產(chǎn)生的水位漲幅速率遠(yuǎn)大于關(guān)閉完成后的水位漲幅速率。并以此為依據(jù)設(shè)置了不同的下游關(guān)閘時(shí)間，分析上、下游節(jié)制閘同步動(dòng)作情況下關(guān)閘時(shí)間對(duì)閘前最大雍高產(chǎn)生時(shí)間與雍高高度的影響，結(jié)果表明最大雍高出現(xiàn)在上游水波傳遞到下游的時(shí)間或者滯后一小段時(shí)間，且下游關(guān)閘時(shí)間是影響水位雍高的主導(dǎo)因素。通過(guò)對(duì)多渠池聯(lián)調(diào)下其他渠池的閘門(mén)調(diào)控對(duì)研究渠池的水位影響分析，得出了同步關(guān)閘情況下渠池雍水基本只受本渠池上、下游閘門(mén)動(dòng)作的影響的結(jié)論，從而為多閘聯(lián)調(diào)情況下閘前水位優(yōu)化提供了參考。

關(guān)鍵詞：串聯(lián)渠池;同步關(guān)閘;雍水;關(guān)閘時(shí)間

中圖分類(lèi)號(hào)：TV66文獻(xiàn)標(biāo)志碼：A

開(kāi)放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：趙鳴雁

Analysis of the influence of gate closure time on the height of the water level before the control gate in a multi-channel system

ZHAO Mingyan.1，KONG Lingzhong.2，ZHENG Yanxia.3，LEI Xiaohui.4，QUAN Jin.4

（1.Construction and Administration Bureau of South-to-North Water Diversion Middle Route Project，Beijing 100038，China;2.College of Civil Engineering and Architecture，Zhejiang University，Hangzhou 310058;3.Beijing Department of South-to-North Diversion，Beijing 100097，China;4.China Institute of Water Resources and Hydropower Research，Beijing 100038，China）

Abstract：In view of the problem of water level rising caused by emergency synchronous closure of control gates，this paper analyzes the influence of different closure time on the water level under the condition of synchronous closure of multiple gates.We established a one-dimensional hydraulic simulation model，conducted plenty of calculations on closure conditions，and studied the influence of the closure of upstream and downstream control gates on the water level before the gate.The results showed that the water level rise rate during the closure process of the downstream control gate was much greater than that after the completion of the closure.On this basis，we set different downstream gate closure time，and analyzed the influence of the closure time on the height and happening time of the maximum water level during the synchronous closure of upstream and downstream control gates.The results showed that the maximum water level appeared during the transmission of the upstream water wave to the downstream section or a short time after that，and the downstream gate closure time was the dominant factor affecting the water level.After analyzing the influence of gate operation in other pools on the water level of the study pool in a multi-channel system，we concluded that the water level of a pool is basically affected by the operation of the upstream and downstream gates of the pool alone under the condition of synchronous closure.This paper can provide a reference for the optimization of the water level before the gate under synchronous closure.

Key words： multi-channel system;synchronous closure;water level rise;gate closure time

為了解決水資源分布不均的問(wèn)題，國(guó)內(nèi)外已建或在建有許多的大型調(diào)水工程，比如南水北調(diào)工程、美國(guó)的加州調(diào)水工程、法國(guó)的普羅旺斯工程。這些輸水工程大部分輸水段都是采用明渠輸水，并采用明渠兩端的節(jié)制閘進(jìn)行輸水調(diào)控[1-2]。正常運(yùn)行過(guò)程中，明渠的運(yùn)行工況較為穩(wěn)定，一般其分水流量較小，引起的渠道流量或者水位變化很小。因此其節(jié)制閘控制目標(biāo)較為簡(jiǎn)單，達(dá)到輸水目標(biāo)流量并且保證水位盡可能穩(wěn)定在目標(biāo)水位即可[3]。目前針對(duì)渠道運(yùn)行控制方法主要為蓄量補(bǔ)償算法或者預(yù)測(cè)控制算法[4-5]。這些算法都是假設(shè)渠池的流量變化需求提前知道的情況下進(jìn)行的閘門(mén)控制策略計(jì)算[6-7]。但是，當(dāng)工程遇到可能發(fā)生的應(yīng)急事件，比如氣候、洪水、地震、滑坡、建筑物損害、人為失誤和交通事故，就需要進(jìn)行突然的輸水流量調(diào)控[8]。最為嚴(yán)重工況就是在發(fā)生下游段突發(fā)水污染事件的時(shí)候，整個(gè)上游段渠池都需要進(jìn)行大幅度的流量減小[9]。而閘門(mén)的關(guān)閉勢(shì)必會(huì)在渠池下游產(chǎn)生雍水，如果渠池下游雍水過(guò)高則可能導(dǎo)致漫溢。而漫溢會(huì)導(dǎo)致渠池水體淹沒(méi)節(jié)制閘機(jī)房，造成液壓油污染或者其他類(lèi)型的污染。漫溢在渠池的調(diào)控過(guò)程中一般是不允許發(fā)生的[10]。因此，在應(yīng)急調(diào)控過(guò)程中，除了考慮渠池的穩(wěn)定問(wèn)題，還需要考慮關(guān)閘導(dǎo)致的節(jié)制閘水位雍高，節(jié)制閘的調(diào)控較為困難[11]。

目前，針對(duì)渠池下游的突發(fā)的大流量變化下的閘門(mén)調(diào)控方法研究較少。國(guó)外學(xué)者主要是研究了在渠池上游來(lái)水突然減小下的應(yīng)急調(diào)控方法，由于上游調(diào)控產(chǎn)生的漲水波需要一定的時(shí)間才能到達(dá)下游渠池，結(jié)果表明可以通過(guò)閘門(mén)異步調(diào)整來(lái)減小水位雍高[12-13]。部分學(xué)者考慮了利用在線水庫(kù)的庫(kù)容來(lái)進(jìn)行突發(fā)水量變化情況下的閘門(mén)調(diào)控[14]，但是并不是所有的渠池都具有在線調(diào)節(jié)水庫(kù)。中國(guó)學(xué)者則針對(duì)南水北調(diào)中線工程的突發(fā)水污染事件，研究了大型明渠的的應(yīng)急調(diào)控方法[15-16]，并得出應(yīng)急渠段的上游適合采用同步關(guān)閘措施，并且結(jié)合退水閘的使用來(lái)保證工程的安全的結(jié)論。雖然通過(guò)退水閘的使用可以減小渠道雍水，但是退水閘的開(kāi)啟及其不經(jīng)濟(jì)，會(huì)退掉大量的可飲用水。國(guó)內(nèi)已有研究表明在節(jié)制閘的開(kāi)度調(diào)節(jié)過(guò)程中可以通過(guò)延長(zhǎng)關(guān)閘時(shí)間來(lái)減小水位變幅速率[17]。如果在節(jié)制閘的關(guān)閉過(guò)程中，考慮節(jié)制閘的關(guān)閉時(shí)間對(duì)節(jié)制閘的閘前雍水高度的影響，則可以通過(guò)合適的關(guān)閘時(shí)間選擇來(lái)盡可能的降低閘前雍水高度。因此，本文通過(guò)研究節(jié)制閘關(guān)閘時(shí)間對(duì)閘前雍水的影響，分析在應(yīng)急情況下的閘門(mén)關(guān)閘時(shí)間選擇。

1建模方法

（1）渠系建模方程。

對(duì)于輸水明渠而言，其渠道多為長(zhǎng)窄類(lèi)型[18]，因此，可以采用一維水力學(xué)模型來(lái)進(jìn)行較精確模擬。一維方程為圣維南方程組[19]：[JB（{][SX（]?A[]?t[SX）]+[SX（]?Q[]?x[SX）]=q[SX（]?[]

式中：x和t為空間和時(shí)間坐標(biāo);A為過(guò)流面積;Q為流量;h為水深;S0為渠道底坡;g為重力加速度;q為單位長(zhǎng)度渠道上的側(cè)向出流量;Sf為摩擦坡度，定義為

Sf=Q|Q|/K.2（2）

其中K為流量模數(shù)。

（2）節(jié)制閘建模方程。

節(jié)制閘處流量水位關(guān)系過(guò)閘流量方程描述。一般情況下輸水明渠的節(jié)制閘過(guò)閘流量可以假設(shè)可以用與上游水深有關(guān)的方程描述，方程為：

Q=Cdlb[KF（]2gh1[KF）]（3）

式中：Q為流量;l為閘門(mén)開(kāi)度;b為閘門(mén)寬度;h1為上游水深;Cd為過(guò)閘流量系數(shù)，根據(jù)閘門(mén)開(kāi)度、閘門(mén)過(guò)流形式等而定，其取值可參考Henry公式[20]。

在節(jié)制閘開(kāi)度已知的情況下，整個(gè)模型是具有定解的。在模型計(jì)算中，假設(shè)節(jié)制閘的開(kāi)度是根據(jù)關(guān)閘時(shí)間線性減小的。

對(duì)于明渠水流，一般最下游節(jié)制閘水位不至于過(guò)低以至于影響閘門(mén)過(guò)閘能力，基本為明渠緩流，不存在急緩流交替。因此模型的數(shù)值解可采用Preissman隱式差分格式求解[21]。

2研究情景設(shè)置

閘門(mén)關(guān)閉過(guò)程對(duì)水位的影響主要在于流量改變后形成的水波。其中，下游閘門(mén)關(guān)閉過(guò)程中會(huì)形成漲水波，而上游閘門(mén)關(guān)閉過(guò)程會(huì)形成降水波。由于上游的降水波到達(dá)下游節(jié)制閘存在著時(shí)間上的滯后。因此，下游節(jié)制閘閘前水位在上游降水波傳遞過(guò)來(lái)的過(guò)程中會(huì)先上升后下降從而產(chǎn)生雍水現(xiàn)象。為分析上、下游關(guān)閘時(shí)間對(duì)下游水位的影響，首先構(gòu)造了一個(gè)簡(jiǎn)單兩渠池串聯(lián)渠道，通過(guò)設(shè)置兩個(gè)渠池的上、下游閘門(mén)的關(guān)閉時(shí)間，來(lái)分析關(guān)閘時(shí)間對(duì)渠池下游水位的影響。這里，設(shè)置的兩個(gè)渠池的渠池長(zhǎng)度均為15 km，每個(gè)渠池的最下游水深均為3.1 m;初始的工況為最上游進(jìn)口流量為40 m.3/s，中間分水口流量均為0，最下游的出流也為40 m.3/s。此工況下計(jì)算得到的水波的傳遞速度約為4.7 m/s，水波在每個(gè)渠池中的從上游傳播到下游以及從下游傳播到上游的時(shí)間分別約為50 min和70 min，因此可考慮設(shè)置關(guān)閘時(shí)間為10～90 min來(lái)進(jìn)行分析。研究渠段的具體結(jié)構(gòu)和具體參數(shù)見(jiàn)圖1和表1。

2.1單一渠池中節(jié)制閘動(dòng)作對(duì)渠池下游水位的影響分析計(jì)算工況

為了分析單一渠池中，上、下游節(jié)制閘分別關(guān)閉作用下的閘前水位變化過(guò)程，這里選取了渠池1作為研究對(duì)象。為分析上游節(jié)制閘關(guān)閘對(duì)下游水位的影響，假設(shè)節(jié)制閘0的流量從40 m.3/s變?yōu)? m.3/s，而其他節(jié)制閘的開(kāi)度保持不變;為了分析下游節(jié)制閘關(guān)閘對(duì)下游水位的影響，假設(shè)節(jié)制閘1的流量從40 m.3/s變?yōu)? m.3/s，而其他節(jié)制閘的開(kāi)度保持不變。此外，為分析關(guān)閘時(shí)間不同導(dǎo)致的水位變動(dòng)差異，設(shè)置了10 min、30 min兩種關(guān)閘時(shí)間，即節(jié)制閘0或節(jié)制閘1的流量分別在10 min和30 min內(nèi)線性從40 m.3/s變?yōu)? m.3/s。此外，設(shè)置初始的變化開(kāi)始時(shí)間均為10 min。計(jì)算工況見(jiàn)表2。

2.2單一渠池中上、下游節(jié)制閘同時(shí)動(dòng)作對(duì)渠池下游雍水的影響分析計(jì)算工況

而在上、下游節(jié)制閘同時(shí)關(guān)閉的過(guò)程中，由于上游漲水波到達(dá)下游節(jié)制閘前的滯后性，會(huì)導(dǎo)致閘前水位的雍高。為分析在單一渠池中，上、下游節(jié)制閘共同作用下的渠池下游閘前水位雍高過(guò)程，這里選取了渠池1作為研究對(duì)象。假設(shè)節(jié)制閘0以及節(jié)制閘1同時(shí)開(kāi)始關(guān)閘，而節(jié)制閘2的開(kāi)度保持不變，分析節(jié)制閘1的閘前水位雍高過(guò)程。為分析節(jié)制閘關(guān)閘時(shí)間對(duì)閘前水位雍高的影響，同樣需要設(shè)置兩個(gè)不同的關(guān)閘時(shí)間。而由于上游水波傳播時(shí)間到達(dá)下游的時(shí)間具有滯后性，這里計(jì)算滯后時(shí)間大約為50 min。因此，節(jié)制閘1設(shè)置了70 min關(guān)閘，以分析下游節(jié)制閘關(guān)閘時(shí)間大于水波傳播時(shí)間情況下的水位雍高過(guò)程。設(shè)置計(jì)算工況見(jiàn)表3。

2.3串聯(lián)渠池其他渠池節(jié)制閘動(dòng)作對(duì)渠池1下游水位的影響分析計(jì)算工況

而在串聯(lián)渠池中，一般多個(gè)渠池處于聯(lián)動(dòng)狀態(tài)，這個(gè)時(shí)候就需要考慮多渠池節(jié)制閘關(guān)閘對(duì)閘前水位的影響。對(duì)于某一渠池而言，其下游閘前水位除了受渠池上、下游節(jié)制閘影響，也受到上游渠池的節(jié)制閘和下游渠池的節(jié)制閘動(dòng)作影響。以渠池1中的節(jié)制閘1的閘前水位為研究對(duì)象，其即受到節(jié)制閘0上游的節(jié)制閘的關(guān)閘動(dòng)作的影響，也受到節(jié)制閘2以及節(jié)制閘2下游的關(guān)閘動(dòng)作的影響。但是由于節(jié)制閘0上游的節(jié)制閘動(dòng)作產(chǎn)生的水波到達(dá)節(jié)制閘1處要需要傳播兩倍的渠池長(zhǎng)度，其影響滯后時(shí)間很[CM（22]長(zhǎng)，而且在傳播過(guò)去時(shí)節(jié)制閘0的降水波已經(jīng)與節(jié)制閘1產(chǎn)生的漲水波相遇，因此其對(duì)最大水位雍高幾乎沒(méi)有影響。所以這里就只考慮節(jié)制閘2的關(guān)閘對(duì)渠池1中節(jié)制閘1的閘前水位影響。

首先分析閘門(mén)2采取關(guān)閘動(dòng)作，而渠池1中的節(jié)制閘0和節(jié)制閘1開(kāi)度保持不變工況下，閘門(mén)2關(guān)閘對(duì)節(jié)制閘1的閘前水位的影響。設(shè)置工況見(jiàn)表4。

2.4串聯(lián)渠池閘門(mén)都關(guān)閉其他渠池節(jié)制閘動(dòng)作對(duì)渠池1下游雍水的影響分析計(jì)算工況

此外，還需要考慮其他渠池關(guān)閘對(duì)本渠池中的水位雍高過(guò)程的影響。由于節(jié)制閘2的關(guān)閉對(duì)節(jié)制閘1的閘前水位的影響時(shí)間的滯后性較大，如果節(jié)制閘1的關(guān)閘時(shí)間很短，下游水波傳播上來(lái)時(shí)節(jié)制閘1已經(jīng)關(guān)閉，其對(duì)節(jié)制閘1的閘前水位影響幾乎可忽略。因此這里選擇較長(zhǎng)的節(jié)制閘1關(guān)閘時(shí)間。此外還設(shè)置了對(duì)照組：節(jié)制閘0和節(jié)制閘1關(guān)閉，而節(jié)制閘2的開(kāi)度保持不變。以分析在節(jié)制閘同步關(guān)閘過(guò)程中，其他渠池的節(jié)制閘關(guān)閉對(duì)本渠池的閘前水位的影響。具體工況設(shè)置見(jiàn)表5。

采用第2部分的建模思路建立水力學(xué)模擬模型，按照各工況設(shè)定的節(jié)制閘運(yùn)行工況，進(jìn)行非恒定流模擬，分析各工況中節(jié)制閘1的閘前水位變化過(guò)程。假設(shè)初始時(shí)刻渠池處于恒定流運(yùn)行工況，所有的關(guān)閘都是從10 min開(kāi)始。

3結(jié)果分析

3.1單一渠池中節(jié)制閘動(dòng)作對(duì)渠池下游水位的影響

首先分析上、下游關(guān)閘對(duì)渠池下游水位的影響。設(shè)置節(jié)制閘0關(guān)閘時(shí)間分別為10 min和30 min，其節(jié)制閘0的過(guò)閘流量變化過(guò)程與節(jié)制閘1的閘前水位變化過(guò)程見(jiàn)圖2。設(shè)置節(jié)制閘1關(guān)閘時(shí)間分別為10 min和30 min，其節(jié)制閘1的過(guò)閘流量變化過(guò)程與節(jié)制閘1的閘前水位變化過(guò)程見(jiàn)圖3。

在圖2中，節(jié)制閘0在10 min開(kāi)始關(guān)閉，而下游節(jié)制閘閘前水位則在56 min分鐘才開(kāi)始變動(dòng)，說(shuō)明中間間隔著46 min左右的滯后時(shí)間。這說(shuō)明上游關(guān)閘產(chǎn)生的降水波需要一段時(shí)間才能到達(dá)下游閘前。一般渠池上游擾動(dòng)傳遞到下游的時(shí)間為水波的傳遞時(shí)間[22]。對(duì)于上游擾動(dòng)，水波傳播時(shí)間計(jì)算公式為t=[SX（]L[]c+v[WTBX][SX）]，水波波速也用c=[KF（]gh[KF）]近似估計(jì)。在本文設(shè)置工況中，渠池上、下游的水波波速計(jì)算結(jié)果分別為4.1 m/s和4.6 m/s，上、下游水流流速分別為0.8 m/s和0.5 m/s，上游水波傳播到下游時(shí)間計(jì)算為47～50 min。這也說(shuō)明模型的計(jì)算結(jié)果較為準(zhǔn)確。

而且對(duì)比10 min和30 min關(guān)閘時(shí)間下的閘前水位變化過(guò)程可以看出，其產(chǎn)生的水位下降大致是兩階段的線性變化：（1）第1階段和第2階段水位都是線性降低，且第1階段的水位降速略大于第2階段;（2）第2階段不同關(guān)閘時(shí)間帶來(lái)的降速基本接近，水位降幅的差別是在第1階段產(chǎn)生的。

從圖3可以看出，本地節(jié)制閘關(guān)閉后的一段時(shí)間內(nèi)閘前水位大致分為兩階段上升：（1）第1階段和第2階段水位都是線性降低，且第1階段的水位降速略大于第2階段;（2）關(guān)閘結(jié)束后的第2階段的水位漲速也是相同的。其中，關(guān)閘完成之后水位線性增加且增速基本保持一致，符合現(xiàn)階段在渠道集中控制中常用的積分時(shí)滯模型[23-24]的相關(guān)結(jié)論，也說(shuō)明了本文模型模擬結(jié)果的準(zhǔn)確。

而對(duì)比圖2和圖3可以看出，下游節(jié)制閘在關(guān)閘過(guò)程中（關(guān)閘第1階段）引起的水位變動(dòng)更為強(qiáng)烈，水位變化速率遠(yuǎn)大于下游節(jié)制閘關(guān)閘完成后的變化速率以及上游節(jié)制閘關(guān)閘引起的下游渠池水位變化速率。而上游節(jié)制閘關(guān)閘引起的渠池下游水位下降速率略大于或者等于下游節(jié)制閘關(guān)閘完成后引起的水位上升速率。

[BT（3+*6]3.2單一渠池中上、下游節(jié)制閘同時(shí)動(dòng)作對(duì)渠池下游雍水的影響

假設(shè)下游節(jié)制閘的關(guān)閘時(shí)間分別為10 min、30 min和70 min。分析在不同的下游關(guān)閘時(shí)間情況下，上游節(jié)制閘不同的關(guān)閘時(shí)間導(dǎo)致的渠池下游閘前水位變化情況。其結(jié)果見(jiàn)圖4。

從圖4中可以看出，在上、下游節(jié)制閘的共同作用下，渠池下游水位出現(xiàn)波動(dòng)，而閘前的水位雍高大致出現(xiàn)在第一個(gè)波峰處;水位最大雍高出現(xiàn)的時(shí)間大約在60～80 min之間，也就是上游水波傳播到下游節(jié)制閘1滯后一段時(shí)間。即最大雍高主要是由于上、下游水波的傳遞時(shí)間的差值造成的。

同時(shí)，對(duì)比相同上游時(shí)間，下游不同關(guān)閘時(shí)間的情況可以看出，延長(zhǎng)下游關(guān)閘時(shí)間20 min帶來(lái)的水位降幅約為0.06 m;相比之下，在下游相同時(shí)間，縮短上游關(guān)閘使勁20 min帶來(lái)的水位降幅約為0.02 m，這也說(shuō)明下游關(guān)閘時(shí)間是影響下游雍水高度的主要因素。

3.3串聯(lián)渠池其他渠池節(jié)制閘動(dòng)作對(duì)渠池1下游水位的影響

假設(shè)節(jié)制閘0和1開(kāi)度不變，而節(jié)制閘2關(guān)閉，其關(guān)閘時(shí)間分別為10 min、30 min。分析節(jié)制閘2的關(guān)閉對(duì)渠池1的下游水位（節(jié)制閘1的閘前水位）的影響。

從圖5中可以看出，節(jié)制閘2在10 min開(kāi)始關(guān)閉，而節(jié)制閘1閘前水位則在78 min分鐘才開(kāi)始變動(dòng)，說(shuō)明中間間隔著68 min左右的滯后時(shí)間。這是因此下游關(guān)閘產(chǎn)生的漲水波需要一段時(shí)間才能到達(dá)上游節(jié)制閘，其傳播時(shí)間即為渠池除以水波的速率。計(jì)算水波傳播時(shí)間t=[SX（]L[]c-v[SX）]≈70 min。 而在對(duì)比10 min和30 min對(duì)閘1的閘前水深的影響可以看出，節(jié)制閘2的關(guān)閉也會(huì)造成節(jié)制閘1的閘前水位升高，但是關(guān)閘時(shí)間幾乎不會(huì)影響到水位的變化速率。

在節(jié)制閘0、節(jié)制閘1、節(jié)制閘2分別采用10 min關(guān)閘而其他節(jié)制閘不操作工況下的水位過(guò)程，見(jiàn)圖6。從圖6中可以看出節(jié)制閘0和節(jié)制閘1關(guān)閉對(duì)渠池1的閘前水位影響較為強(qiáng)烈，且后期導(dǎo)致的水位變化速率接近，而節(jié)制閘2的關(guān)閉對(duì)渠池1的閘前水位影響極小，且還存在著很大的滯后時(shí)間。因此，需要分析在節(jié)制閘都同步采取關(guān)閘措施情況下的節(jié)制閘2對(duì)渠池1的閘前水位的影響。

3.4串聯(lián)渠池閘門(mén)都關(guān)閉其他渠池節(jié)制閘動(dòng)作對(duì)渠池1下游雍水的影響

從圖7中可以看出，在節(jié)制閘都采取關(guān)閘過(guò)程的時(shí)候，下游節(jié)制閘2的關(guān)閉對(duì)渠池1的下游閘前水位沒(méi)有任何影響。圖7（a）中，節(jié)制閘1在70 min內(nèi)關(guān)閉，由于節(jié)制閘2產(chǎn)生的水波到達(dá)節(jié)制閘1的時(shí)間為76 min，所以節(jié)制閘2關(guān)閘不會(huì)對(duì)閘1的閘前水位產(chǎn)生影響。而在圖7（b）中，盡管節(jié)制閘1的關(guān)閘時(shí)間為90 min，但是由于下游的水波傳遞到節(jié)制閘1時(shí)，節(jié)制閘1的閘前水位雍高較大，而節(jié)制閘1的閘前水位雍高較小，且節(jié)制閘開(kāi)度很小，處于自由閘孔出流狀態(tài)，計(jì)算方程假設(shè)下游水深對(duì)上游流量與水深沒(méi)有影響[25]，所以下游水波不會(huì)對(duì)節(jié)制閘1的閘前水位產(chǎn)生影響。也就說(shuō)明，在同步關(guān)閘時(shí)候，除非是極端情況，下游渠池的關(guān)閘幾乎是不會(huì)對(duì)上游渠池的關(guān)閘產(chǎn)生的雍水有影響的。而由于上游渠池的關(guān)閘產(chǎn)生的水波到達(dá)下游存在很大的時(shí)間滯后，幾乎也不會(huì)對(duì)渠池的關(guān)閘雍水有影響。因此，可以假設(shè)在多渠池串聯(lián)渠系中，同步關(guān)閘產(chǎn)生的最大雍水只與本渠池的上、下游的節(jié)制閘動(dòng)作有關(guān)。

4結(jié)論

（1）渠池上游節(jié)制閘關(guān)閉對(duì)下游水位的影響具有滯后性，滯后時(shí)間為上游水波的傳播時(shí)間。而且其導(dǎo)致的渠池下游水位是呈兩階段線性變化的，關(guān)閘時(shí)間越快第1階段的降速越快，而第2階段的降速在不同關(guān)閘時(shí)間下基本相同。由于第1階段的降速差別導(dǎo)致最終的水位降幅差別。

（2）渠池下游節(jié)制閘關(guān)閉對(duì)渠池下游水位的影響是瞬時(shí)的。而且其導(dǎo)致的渠池下游水位的也是兩階段線性變化的。第一階段和關(guān)閘過(guò)程同步，且第一階段的漲速快于第二階段的。第2階段的漲速在不同關(guān)閘時(shí)間下基本相同。由于第1階段的漲速差別導(dǎo)致最終的水位雍高幅度差別。

（3）在渠池中上、下游閘門(mén)都關(guān)閉的情況，最大水位雍高出現(xiàn)在上游水波傳播到下游的時(shí)間或者之后一小段時(shí)間。且最大水位雍高隨著下游關(guān)閘時(shí)間的延長(zhǎng)以及上游關(guān)閘時(shí)間的縮短而逐漸降低，但是下游關(guān)閘時(shí)間是影響最大雍高的主要因素。

（4）在多串聯(lián)渠池同步關(guān)閘情況，其他渠池的閘門(mén)操作幾乎不會(huì)對(duì)其他渠池的水位雍高產(chǎn)生影響。各個(gè)渠池的最大水位雍高主要受渠池上、下游水位關(guān)閘時(shí)間的影響。在各個(gè)渠池的允許水位雍高不同的情況下，可以利用關(guān)閘時(shí)間優(yōu)化來(lái)盡可能使得多渠池的水位雍高滿足要求。

參考文獻(xiàn)（References）：

[1]皮鈞，熊雁暉.加利福尼亞調(diào)水工程對(duì)我國(guó)調(diào)水工程的啟示[J].南水北調(diào)與水利科技，2004，2（4）：50-52.（PI J，XIONG Y H.Enlightenment from California State Water Project to South-North Water Diversion Project[J].South-To-North Water Transfers And Water Science & Technology，2004，2（4）：50-52.（in Chinese））DOI ：10.3969/j.issn.1672-1683.2004.04.017.

[2]王長(zhǎng)德，阮新建.南水北調(diào)中線總干渠控制運(yùn)行設(shè)計(jì)[J].人民長(zhǎng)江，1999，30 （1）：19-21.（WANG C D，RUAN X J.Design of automatic control operation for main canal of Middle Route S-N Water Transfer Project[J].Yangtze River，1999，30 （1）：19-21.（in Chinese）） ??DOI ：10.3969/j.issn.1001-4179.1999.01.007.

[3]崔巍，陳文學(xué)，姚雄，等.大型輸水明渠運(yùn)行控制模式研究[J].南水北調(diào)與水利科技，2009，7（5）：6-10.（CUI W，CHEN W X，YAO X，et al.Research on canal control of large scale water transfer project[J].South-To-North Water Transfers And Water Science & Technology，2009，7（5）：6-10.（in Chinese）） DOI ：10.3969/j.issn.1672-1683.2009.05.002.

[4]BAUTISTA E，CLEMMENS A J.Volume compensation method for routing irrigation canaldemand changes[J].Journal of Irrigation and Drainage Engineering，2005，131（6）：494-503.DOI ：10.1061/（ASCE） 0733-9437（2005）131：6（494）.

[5]ALVAREZ1 A，RIDAO M A，RAMIREZ D R，et al.Constrained predictive control of an irrigation canal[J].Journal of Irrigation and Drainage Engineering，2013，139（10）：841-854.DOI：10.1061/（ASCE）IR.1943-4774.0000619.

[6]姚雄，王長(zhǎng)德，丁志良，等.渠系流量主動(dòng)補(bǔ)償運(yùn)行控制研究[J].四川大學(xué)學(xué)報(bào)（工程科學(xué)版），2008，40（5）：38-44.（YAO X，WANG C D，DING Z L，et al.Study on canal system operation based on discharge active compensation[J].Journal of Sichuan University（Engineering Science Edition） 2008，40（5）：38-44.（in Chinese）） DOI：10.15961/j .jsuese.2008.05.018.

[7]尚毅梓，吳保生，李鐵鍵，等.渠道分水?dāng)_動(dòng)可預(yù)知算法設(shè)計(jì)與仿真[J].水科學(xué)進(jìn)展，2011，22 （2）：242-248.（SHANG Y Z，WU B S，LI T J，et al.Design and simulation of a foreseeable algorithm for canals[J].Advances in Water Science.2011，22 （2）：242-248.（in Chinese）） ???DOI：10.14042/j.cnki.32.1309.2011.02.004.

[8]SEWELL T，STEPHENS R E，DOMINEYHOWES D，et al.Disaster declarations associated with bushfires，floods and storms in New South Wales，Australia between 2004 and 2014[J].Scientific Reports，2016？，？6.DOI：10.1038/srep36369.

[9]CHENG W P，JIA Y .Identification of contaminant point source in surface waters based on backward location probability density function method[J].Advances in Water Resources，2010，33 （4）：397-410.DOI：10.1016/j.advwatres.2010.01.004

[10]方神光，吳保生，趙剛.南水北調(diào)中線輸水渠道的漫溢現(xiàn)象分析[J].水利水電科技進(jìn)展，2007，27 （4）：11-14.（FANG S G，WU B S，ZHAO G.Analysis of overflowing phenomenon in water transfer channel of middle route of South-to-North Water Transfer Project[J].Advances in Science and Technology of Water Resources.2007，27 （4）：11-14.（in Chinese）） DOI：10.3880/j.issn.1006-7647.2007.04.003.

[11]王浩，雷曉輝，尚毅梓.南水北調(diào)中線工程智能調(diào)控與應(yīng)急調(diào)度關(guān)鍵技術(shù)[J].南水北調(diào)與水利科技，2017，15 （2）：1-8.（WANG H，LEI X H，SHANG Y Z.Key technologies of intelligent control and emergency regulation for The Middle Route of South-To-North Water Diversion Project[J].South-to-North Water Transfers and Water Science & Technology.2017，15 （2）：1-8.（in Chinese）） ??DOI：10.13476/j.cnki.nsbdqk.2017.02.001.

[12]SOLER J，G MEZ M，RODELLAR J，et al.Application of the GoRoSo feedforward algorithm to compute the gate trajectories for a quick canal closing in the case of an emergency[J].Journal of Irrigation & Drainage Engineering，2013，139（12）：1028-1036.DOI：10.1061/（ASCE）IR.1943-4774.0000640.

[13]方神光，吳保生.南水北調(diào)中線干渠閘前變水位運(yùn)行方式探討[J].水動(dòng)力學(xué)研究與進(jìn)展，2007，22 （5）：633-639.（FANG S G，WU B S.Operation of allowing fluctuating of water level in the Middle Route of the South-to-North Water Transfer Channel[J].Journal of Hydrodynamics.2007，22（5）：633-639.（in Chinese）） DOI：10.3969/j.issn.1000-4874.2007.05.014

[14]SOLER J，GAMAZO P，RODELLAR J，et al.Operation of an irrigation canal by means of the passive canal control[J].Irrigation Science，2015，33 （2）：95-106 ?DOI：https：//doi.org/10.1007/s0027.

[15]房彥梅，張大偉，雷曉輝，等.南水北調(diào)中線干渠突發(fā)水污染事故應(yīng)急控制策略[J].南水北調(diào)與水利科技，2014 （2）：133-136.（FANG Y M，ZHANG D W，LEI X H，et al.Emergency control strategy for sudden water pollution accident in the main channel of Middle Route of South-to-North Water Diversion Project[J].South-to-North Water Transfers and Water Science & Technology.2014 （2）：133-136.（in Chinese）） ???DOI：10.13476/j.cnki.nsbdqk.2014.02.031.

[16]王浩，鄭和震，雷曉輝，等.南水北調(diào)中線干線水質(zhì)安全應(yīng)急調(diào)控與處置關(guān)鍵技術(shù)研究[J].四川大學(xué)學(xué)報(bào)（工程科學(xué)版），2016，48 （2）：1-6.（WANG H，ZHENG H Z，LEI X H，et al.Study on key technologies of emergency regulation and treatment to ensure water quality safety of the main canal of Middle Routes of South-to-North Water Diversion Project[J].Journal of Sichuan University（Engineering Science Edition），2016，48 （2）：1-6.（in Chinese）） ???DOI：10.15961/j.jsuese.2016.02.001.

[17]吳輝明，雷曉輝，秦韜，等.南水北調(diào)中線渠系蓄量補(bǔ)償運(yùn)行控制方式[J].南水北調(diào)與水利科技，2015 （4）：788-792.（WU H M，LEI X H，QIN T，et al.Operation control of storage compensation of canal system in the Middle Route of South-To-North Water Diversion Project[J].South-to-North Water Transfers and Water Science & Technology，2015 （4）：788-792（in Chinese）） ??DOI：10.13476/j.cnki.nsbdqk.2015.04.040.

[18]TANG C，YI Y，YANG Z，et al.Water pollution risk simulation and prediction in the main canal of the South-to-North Water Transfer Project[J].Journal of Hydrology，2014，519：2111-2120.DOI：10.1016/j.jhydrol.2014.10.010.

[19]范杰，王長(zhǎng)德，管光華，等.渠道非恒定流水力學(xué)響應(yīng)研究[J].水科學(xué)進(jìn)展，2006，17 （1）：55-60.（FAN J，WANG C D，GUAN G H，et al.Study on the hydraulic reaction of unsteady flows in open channel[J].Advances in Water Science.2006，17 （1）：55-60.（in Chinese）） DOI：10.3321/j.issn：1001-6791.2006.01.009.

[20]HENRY H R.Discussion of diffusion of submerged jets by Albertson M L，Dai Y B，Jensen R A，Rouse H[J].Transaction of ASCE，1950，115：687-694.

[21]田兆偉.一維非恒定流數(shù)學(xué)模型參數(shù)若干問(wèn)題研究[J].人民珠江，2015，36 （3）：16-20.（TIAN Z W.Research on some parameters of one-dimensional unsteady flow mathematical model parameters[J].Pearl River，2015，36（3）：16-20.（in Chinese）） DOI：10.3969/j.issn.1001-9235.2015.03.005.

[22]崔巍，陳文學(xué)，穆祥鵬，等.明渠運(yùn)行前饋控制改進(jìn)蓄量補(bǔ)償算法研究[J].灌溉排水學(xué)報(bào)，2011，30 （3）：12-17.（CUI W，CHEN W X，MU X P，et al.Revised Volume-compensation algorithm for feedforward control of canal[J].Journal of Irrigation and Drainage.2011，30 （3）：12-17.（in Chinese））

[23]CLEMMENS A J，SCHUURMANS J.Modeling of irrigation and drainage canals for controller design[J].Journal of Irrigation & Drainage Engineering，1999，125（6）：338-344.DOI：10.1061/（ASCE）0733-9437（1999）125：6（338）.

[24]CLEMMENS A J，TIAN X，OVERLOOP P J V，et al.Integrator delay zero model for design of upstream water-level controllers[J].Journal of Irrigation & Drainage Engineering，2017，143 （3）：B4015001.DOI：10.1061/（ASCE）IR.1943-4774.0000997.

[25]毛潭，李繼棟，史喆瓊，等.平板閘孔自由出流流量系數(shù)的實(shí)驗(yàn)研究[J].實(shí)驗(yàn)科學(xué)與技術(shù)，2015，13 （3）：6-9.（MAO T，LI J D，SHI Z Q，et al.Experimental research of flow coefficient for free outflow of the flat sluice[J].Experiment Science and Technology.2015，13 （3）：6-9.（in Chinese）） DOI：10.3969/j.issn.1672-4550.2015.03.003.