闸门调度运行模式探讨I.pdfVIP

南水北调中线干渠闸门调度运行方式探讨 方神光，吴保生，傅旭东 (清华大学 水沙科学与水利水电工程国家重点实验室，水利水电工程系，北京 100084) 摘 要：本文选取天津分水口为例，采用建立的南水北调中线电子渠道模型，模拟计算了该 分水口分水流量变化时分别采用时序控制和同步控制进行闸门调度时干渠中的水流过渡过程。 结果显示，时序控制较同步控制可以大幅度缩短干渠水力过渡时间并减小闸前水位变化幅度， 但闸下流量和闸门开度变化幅度较同步控制要大。若采用分水口之间的分水流量进行区间调度 以维持渠首流量不变，应用时序控制时可以更进一步缩短干渠水力过渡时间。闸前常水位运行 模式下，不管采用哪种调度方式，都会影响到进入惠南庄前池的流量，因此这要求前池具有一 定的调蓄容积以避免水位大幅度的变化。 关键词：南水北调中线工程；非恒定流；数学模型；时序控制；同步控制； 中图分类号：TV 133.2 文献标识码：A 南水北调中线工程从丹江口水库陶岔引水，最终至北京团城湖，输水总干线全长 1 267 km， 渠首设计引水流量 350 m3/s 。由于中线主干渠流量大、线路长，在保证安全稳定输水的前提下， 如何调度各节制闸开度进行适时适量供水是即将建成的中线渠道面临的主要问题之一。王长德 和阮新建[1]针对中线总干渠的控制运行指出，中线工程有必要考虑以下游控制运行为主的自动 控制方式，由于实现该目标涉及到各学科的理论和方法，因此需要组织各学科专家协调攻关。 但不管是采用何种自动控制方式，其最根本依托点是对主干渠中水流波动传播过程的了解和掌 握，以此通过联合调度各闸门来对其进行响应。进行闸门联合调度的首要条件是确定中线输水 渠道的运行方式，考虑到渠道安全和稳定供水等因素，长委勘测规划设计研究院推荐中线干渠 采用下游常水位运行方式[2]，各节制闸前常水位值的大小主要根据水头分配以及恒定流水力计 算来确定；另外，多个闸门的调度运行模式主要有三种方法：顺序控制法、同步控制法和时序 控制法，由于顺序控制法可能会造成水量传递缓慢和水位的大幅度波动，一般仅适用于需水变 化简单的中小型输水系统[3]，对中线输水干渠而言，应用中央控制系统对干渠所有闸门实现同 步控制或时序控制更有利于渠道安全运行以及对渠道中的流量变化进行快速响应。因此，本文 在长委勘测规划设计研究院推荐的常水位运行方式基础上，设想为响应天津分水口分水流量增 大或减小时，对中线输水干渠各节制闸调度分别采用同步控制法和时序控制法的优缺点进行了 ________________________ 收稿日期： 基金项目：中国博士后科学基金资助项目（20060390078 ）； 国家自然科学创新研究群体科学基金资助项目） 作者简介：方神光（1978-），男（汉），湖北，博士后。 E-mail: fangsg@mail.tsinghua.edu.cn 模拟计算、分析和比较。所得研究成果可供即将建成的南水北调中线干渠的闸门调度运行和管 理进行借鉴和参考。 1 数学模型 中线输水渠道平均底坡为 1/25 000，属于缓坡渠道，并认为受节制闸启闭或取水口流量变 化影响的渠道水流流态为非恒定渐变流，其流动特性可采用水位z 、流量Q和时间t表示的 Saint-Venant方程组表述，主要由连续方程和运动方程构成[4]，方程具体形式可参考各文献和资 料。此处采用上下游流量严格守恒的Preismann 四点差分隐格式对圣维南方程组进行离散化处 理[3]，将圣维南方程组转化为一组差分方程组，采用追赶法进行迭代求解。所用非恒定流求解 计算软件为作者采用FORTRAN语言编写。 2 边界条件 2.1 渠道两端的边界条件 模拟计算的中线输水渠段起始位置为陶岔渠首，终止位置是河北渠段终点（进入惠南庄前 池），渠道中通过 70% 的设计流量，各分水口分水流量也为设计分水流量的 70%，模拟计算各 工况时认为渠道末端水位保持不变。 2.2 内边界条件的处理方法 中线输水渠道中存在大量的节制闸、分水口和退水口，在采用追赶法求解非恒定圣维南方 程组的离散代数方程组时，为使追赶法计算能自始至终的进行下去，必须根据具体水力特性将 其作为内部边界条件进行特殊处理。内边界条件通常都包含两个相容条件：流量的连续性条件 和能量守恒条件（或动量守恒条件）[5]，应用这两个条件即可将节制闸、分水口和退水口的影 响恰当的反应在计算程序中。对闸门的出流流态和