基于生活、生产和生态环境用水的水资源配置模型

：从实际需要出发探讨并给出了水资源系统及水资源配置模型的基本概念和定义；研究了水资源系统网络图的绘制方法；给出了用该图概化和表示系统中的需水、水源、工程以及水的运动等复杂关系的方法；从时间结构和空间结构两个方面介绍了全面考虑生活用水、生产用水和生态环境用水要求的、系统反映各种水源及工程供水特点的水资源配置模型的建模思路和技巧，给出了模型的基本任务和主要约束方程．该模型在河南省安阳市水资源可持续利用综合规划中进行了实际应用，推荐的配置方案合理，计算速度很快．该模型可以应用于大型复杂水资源系统．
　　关键词：水资源系统；水资源配置模型；可持续利用；安阳市　　

　　当前国际上水资源系统规划的发展趋势，一是对水资源系统结构、供求关系的描述更深入和具体，更多地采用水文长系列分析方法；二是从区域社会、经济、环境及生态持续发展的动态角度研究水资源工程布局与供需平衡〔1〕．自20世纪80年代初以来，我国开展了一系列大型水资源项目，有力地推动了水资源系统分析理论方法和模型技术的发展和应用．这些项目所建立的复杂水资源配置模型，主要研究的是如何满足生活、生产需水要求，没有同时考虑生活、生产和生态环境的需水要求〔2～6〕．
　　本文模型是在中国水利水电科学研究院研制的水资源配置模型（它反映的因素相对比较全面，在我国应用最广泛）的基础上改进的，全面考虑了生活、生产和生态环境的需水要求．

1　水资源配置模型的基本概念及水资源系统网络图
1．1　水资源系统及水资源配置模型的概念
1．1．1　水资源系统的定义
　　《中国资源百科全书》（2000年）对水资源系统的定义是：在一定区域内由可为人类利用的各种形态的水所构成的统一体．显然这一定义是局限于水的范畴定义的．本文从实际需要出发，给出水资源系统的如下定义：在一定区域内由可为人类及生态环境利用的各种形态的水以及与水存在、运动、演变有关的物体和社会所构成的系统．
1．1．2　水资源配置模型的定义
　　本文对水资源配置模型的定义是：针对以供水为主要目的的水资源系统，以系统分析理论、运筹学方法、知识规则、逻辑推理等为技术基础，对各种工程措施和非工程措施进行适当组合和合理的联合调度运用，以追求系统整体的可持续利用功能最优为目标的计算机模型．
　　该定义有几个要点：首先，水资源配置模型是计算机模型而不是物理模型，其主要用途是合理地安排和调度水资源开发利用的工程措施和非工程措施，配置水资源，使水资源系统的可持续利用功能最优．第二，分析的问题可以是规划阶段的问题，其重点是进行多种水利工程措施和非工程措施组合方案的分析比较，推荐出合理的水利工程布局；也可以是运行管理阶段的问题，其重点是在已有水资源系统条件下合理地调度各种水源和工程使其发挥最大效益．第三，对“水资源系统整体的可持续利用功能最优为目标”要广义地理解，充分体现实际需要和现实可行性．
1．2　水资源系统概化及系统网络图
　　建立水资源配置模型需要先对实际水资源系统进行抽象和概化，并绘制系统网络图．水资源配置模型可概化由五部分构成：①基本物理元素（集合）；②各物理元素量度数据（参数及变量）；③物理元素之间的相互关系（约束条件）和系统协调准则（调度规则和目标函数）；④解决问题的方法；⑤各物理元素在此基础上所处的状态（结果）．
　　水资源系统中各类物理元素（水库工程、拦河闸、引水工程或扬水站、水资源计算单元、河渠道交汇点等）概化为节点，各结点间通过线段（河道或渠道等）连接，形成水资源系统概化网络图．在平面上不方便表示地下水库，则假定每个计算单元都有一个地下水库，统一不在网络图中标出．如果某单元实际上没有地下水，则可令其地下水库的库容和水量为零．水资源配置系统方案的生成以及对每一方案的详细供需平衡分析及流域水循环分析等都以水资源系统网络图为基础．凡是图上标明的水利工程、分区、河流及控制点或断面的有关问题，模型就能够给出详细的分析结果；如果图上没有，一般来说模型就不能给出具体结果．

2　水资源配置模型
2．1　建模的基本思路
　　在水资源系统描述方面，该模型以系统网络图为基础，采用了多水源（地表水、地下水、外调水及污水处理回用水）、多工程（蓄水工程、引水工程、提水工程、污水处理工程以及闸堰湖泊洼淀等）、多传输系统（包括当地地表水传输系统、外调水传输系统、弃水污水传输系统和地下水的侧渗补给与排泄系统）的描述法，使系统中的各种水源、水量在各处的调蓄情况及来去关系都能够得到客观的、清晰的描述，为得到正确结果打基础．
　　从模型的空间结构方面看，要抓住两条主线索：一条是水资源循环运动的线索，即沿流域自上而下，整个研究区域由各个流域及其不同层次的子流域组成．各种水资源的调查评价也是按照流域分区进行的．另一条是社会、经济、生态环境的需水线索，人口、经济、需水、用水、耗水、水利工程等统计资料多是依照各级行政分区统计的．水资源配置模型需要按照流域分区和行政分区获取资料和出成果．按照这两条线索相结合的要求划分计算单元．
　　从模型的时间结构方面看，也是抓住两条主线索：一条是反映水文变化规律的时间线索，也就是历史水文资料系列，即从哪一年到哪一年以及年内的时间层次（如月、旬、日等，根据精度需要和资料获得的可能性确定）；另一条是反映水资源利用规划与管理需要的时间线索，也就是规划水平年以及年内的时间层次．按照这两条线索相结合的要求确定计算时段（最短的时间层次）．
　　在运行方案上考虑了对各个工程组合方案，各水平年的需水量、来水量的预测，以及污水处理与回用能力、节水水平、地下水可开采量、工程运行规则、各种参数等．
　　在结果分析上包括了各系统元素的水平衡分析、系统内各分区的供水能力及供水量分析、供水效益分析、水源利用情况、弃水情况、污水排放情况、工程分水情况、河流与渠道过流情况、河道水质状况、系统发电量等，并对各模拟计算方案进行综合分析比较，寻找出合理可行的规划方案．
　　在软件实现方面，采用GAMS语言编制分析计算核心模块，采用Foxpro语言编制用户界面以及实现数据库管理．
2．2　模型的基本任务
　　水资源配置模型是研究和解决水资源问题的一种重要模型．它的主要任务是：①系统配置方案的生成；②每一系统方案的优化调度及供需平衡分析；③每一方案的水平衡分析；④各个系统方案的综合比较，推荐最佳的配置方案；⑤诊断水资源系统的薄弱环节等等．
2．3　主要约束方程
　　下面给出水资源配置中，任意一年水资源供需平衡设计的主要约束方程．方程中关键符号的含义为：i代表水库；j代表计算单元及与单元对应的地下水库；k代表渠道或河流节点；t代表计算时段；l代表渠道或河流；s代表用水种类，s＝1，2，3，4，5依次为城市用水（包括生活用水和工业用水）、农业用水、河道外生态环境用水、河道内生态环境用水、发电用水；α为降雨入渗补给系数；β为灌溉入渗补给系数；γ为地表供水河流渠道渗漏补给系数；φ为潜水蒸发系数；ls（），ld（），lo（），lg（）依次为地表水供水河流渠道集合、外调水供水河流渠道集合、弃水河流渠道集合、地下水联系集合，括号内有两个参数，分别为起点和终点；u（），d（）分别表示上游对象集合和下游对象集合，括号内有1个参数，为对象，如水利工程、计算单元、节点等．式中没有被运算符号隔开的连续几个字母表示同一个变量，变量名第4个字母A，I，E分别表示农业供水、工业（包括城镇生活）供水和生态供水；V为水库蓄水库容；C为有关系数．
　　a．地表水库当地水平衡方程：
　　　
　

　　c．水库分水量限制．对于有些向多个计算单元供水的地表水库，如果各地区的供水有约定，调度就必须遵守．
　

　　d．计算单元的城市供水平衡方程：
　　

　　f.计算单元的生态供水平衡方程.生态需水分河道内需水与河道外需水，要分别进行平衡．单元河道外生态供水采用方程（6）计算．原则上所有地表水、地下水、外调水、污水回用水及上游弃水、回归水等均可以参加河道外生态需水的平衡．河道外生态需水与生活需水和生产需水要求不一样，生活需水和生产需水要求每一个时段必须平衡，而河道外生态需水弹性要大得多，可以每个时段刚好平衡，也可以不平衡，只要每个时段的供需量在允许的上下限之间，而全年或两三年的供水量与需水量平衡即可．河道内需水的供需平衡约束方程主要有：最小流量约束、年最小过水量约束、河道水质约束等．本模型在安阳市的规划中，对于该市河道内需水选择的是在重点保护河段上建立水量、水质约束方程．方程式比较多、比较复杂，在此省略．
　　
　
2．4　目标函数
　　模型目标函数为