2022年多水源,多用户大型水资源系统优化模型研究2.pdf

期刊文献多水源、多用户大型水资源系统优化模型研究黄昉1许文斌2郑建青3 （1.浙江省水利水电勘测设计院， 杭州 310002 ；2.浙江省水利厅， 杭州 310009 ； 3. 温州市水利水电局，温州 325000 ）摘要：本文提出一种将水资源系统顺序决策问题转换成有约束非线性优化问题的实用模型模拟权重系数模型。以具有多水源、多用户、多级串并联性质的宁波城市供水区水资源系统为研究对象，通过与采用模拟法和多维增量动态规划法（MDDDP）相结合的 “混合模型”之成果进行对比， 表明模拟权重系数模型具有方法简便、成果明瞭之优点。为此类课题的研究开辟了一条新途径。关键词：水资源系统；模拟权重系数模型；水库调度线1 水资源系统的特点对于多水源、多用户、多级串并联的大型水资源系统的优化调度，以往常采用增量动态规划模型1、分层模型2、系列模型3以及模拟技术与多维增量动态规划相结合的混合模型4。为了避免“维数灾”，先将多个水源工程概化合并，计算完后再进行分解。但对于多级串并联之水资源系统，亦只能采用反复迭代，力求逼近最优解。这种水资源系统的特点有：（1）大系统由几个具有相对独立的供、 排水网络的子系统组成。 各子系统的供水工程由河网和水库组成，且存在“多维”问题。有的子系统维数高达十几维。各子系统之间联通成网的水厂管道可相互调配各子系统的蓄水工程供水量。（2）各子系统均有多种水源： 水库蓄水、 非水库控制山区迳流、 平原地表迳流、 可供开采的地下水以及生活、工业和农业回归水。（3）用水部门分：“确保”、一般工业、畜牧业、农业灌溉、环境、航运等六个用水户。其中：“确保用水户”是指对国计民生极为重要，必须予以最大限度保证的用水部门，包括城乡人民生活、绿化、消防和重要工矿用水等。根据重要程度和单方供水效益差异，各用水户有不同的供水保证程度要求，优化成果需提供各水库的各用户供水调度线，以指导实际调度。（4）各子系统的供水工程均为多级串并联。 水库的弃水是河网的水源之一， 而有些水库却是从所在子系统内河网抽水蓄库的“调峰”水库。（5）水库有管道直供水厂、 重要工矿以及放水至河网再供给各用水户两种供水方式。水厂有部分从河网取水、部分由水库直供两种取水方式。2 数学模型与模型求解2.1 建模思路对于复杂的水资源系统，常规分析方法是采用模拟法，该法的优点在于不管系统多么错综复杂，只要事先确定调度原则和选择好有较佳代表性的确定型迳流系列就能顺利得出分析结果。工程应用中，常由有经验的设计人员定出有限套可能的调度原则方案，通过方案比较确定推荐方案，成果清晰、明瞭。但此法存在成果随调度方案确定者的经验多寡而优劣。对供水系统来讲，优化调度的目标是使其供水效益最大或各用水户的缺水损失总值最小（有时简化采用缺水量总和最小）。在实际调度中，一般通过使调节性能好的供水工程多蓄水，调节性能差的多供水，充分利用各供水工程的库容和提高其库容的重复利用率来增加供水量。由于河网和江道型水库汇水面积大， 调蓄性能较差， 应尽量先用该部分蓄水。 若当前时段该部分蓄水能满足各用户需要，除部分直接由水库供水的用户外，水库群不供水；反之，不足部分即为需要水库群提供的水量。由此精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 1 页，共 8 页 - - - - - - - - - - 期刊文献第一作者简介：黄昉（ 1952） ，男，浙江杭州人，高级工程师，从事水利规划工作。可见，河网和江道型水库的调度方式经定性分析即可得以确定，并可得出各时段要求水库群供水的总需供水量。假定水库 r 时段 i 的需供水量 gri等于水库群总需供水量Gi乘上供水权重系数 ri，据此将水库群各时段的总需供水量分配到各水库， 用模拟方法计算出各用水户在整个计算系列中的总缺水量，再根据使缺水损失总值最小的原则对供水权重系数寻优，从而求得最优解。但ri与时段 i 有关, 仍是顺序决策问题。 若能进一步把对 ri寻优的问题转换成对另一个与时段i 无关、但能确定ri值的参数进行寻优的话，则复杂的系统顺序决策问题将被转换成简单的非线性优化问题。2.2 供水权重系数和优化供水幂数增加供水、减少弃水量的目标是使得各库容系数不同、兴利库容不同的水库能够同时蓄满和同时将蓄水供完。水库时段初蓄水量以及兴利库容与多年平均迳流量的比值库容系数是决定水库当前时段供水量大小的两个重要因素。供水权重系数ri是水库 r 时段 i的需供水量 gri在该时段水库群的总需供水量Gi中所占的比例， 并且 gri应与库容系数 r成反比，与时段初库容 Vri成正比。本文采用ri=gri/lr1gri=（1/ r）*Xr*Vri*Gi/lr1 （1/r）*Xr*Vri*Gi =（1/r）*Xr*Vri/lr1 （1/r）*Xr*Vri （1）式中：Xr为水库 r 的优化供水幂数 ,l 为水库群的水库数。从式（ 1）可见： r越大, ri越小，反之 ri越大，使得库容系数大的水库充分利用其库容多蓄水，库容系数小的多供水来提高其库容的重复利用率；Vri越大, ri越大，反之 ri越小，使得当前蓄水量较大的水库多供水，腾出库容容纳来水，以减少水库弃水。 r与 Vri是两个相互制约的因素， 存在着它们在 ri中各占多少比重可使目标函数最优之优化问题。幂数 Xr起着调整 r与 Vri所占比重的作用。通过对Xr寻优，可得出最佳供水权重系数，从而确定水库群的最优调度方式。由于Xr与时段 i 无关, 通过上述转换，已将系统顺序决策问题转变成了非线性优化问题。2.3 供水调度线为满足各用水户不同供水保证程度的要求，需按n 个用水户设置 n档供水调度线将水库容积自上而下划分为n 个区。根据保证程度要求的高低确定用水户优先级的高低，各用水户供水调度线按优先级从低到高自调度图的上方往下排列。1 区为公共供水区（亦为优先级最低的用水户供水区），当库水位位于此区时按各用水户需求供水；2 区n 区为限制供水区。设用水户K为一般工业用水户， k 区库容是用水户 k 的预留库容，则当水库水位位于k 区时，仅对一般工业及比一般工业用水户优先级高的用水户 （如“确保” 用水户）供水，而中断比一般工业优先级低的各用水户 （如畜牧、农业灌溉、环境、航运等用水户）之水库供水。用水户k 可代表一个或几个具有相同供水保证程度要求的用水户。确定各水库的各用水户预留库容有二个主要问题：一是在哪些水库中设置哪几档预留库容；二是如何确定各用水户供水调度线。2.3.1选择设置预留库容的水库按是否需要管道直供将各用水户分为两类。第一类是有管道直供要求的“确保”用水户，第二类是其余用水户。对于第一类用水户，选择设置“确保”预留库容的水库与“确保”用水户的相对地理位置、 水库性质和库容大小、 “确保”用水量要求以及专供管道的长度、规模、投资等因素有关。可先通过对以上因素进行综合分析，列出可能设置“确保”预留库容的水库组合方案，然后进行技术经济综合比较，选出较优方案作为设置“确保”预留库容的水库推荐方案。对于第二类用水户，由于不涉及供水管道投资比较，可令所有水库均设置相应的预留库容，通过对各档预留库容最大值进行寻优，确定各用水户的供水调度线。2.3.2 供水调度线各时段各用水户的预留库容值与水库的入库迳流、供水量和库容等因素有关。而且，可供给优先级较低用水户之入库迳流量是比其优先级高的用水户用过后剩下的部分。用水户k 各精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 2 页，共 8 页 - - - - - - - - - - 期刊文献时段初总预留库容初始值Vik，按下式逆序倒推确定：精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 3 页，共 8 页 - - - - - - - - - - 期刊文献 Vi+1,k+Gik-Qi （Vi+1,k+GikQi）Vik= (2) O（Vi+1，k+GikQi）式中：Vik、Vi+1,k为时段 i 初和末用水户 k 要求在水库群设置的总预留库容， Gik为时段 i 内用水户 k 要求水库群供水的总供水量， Qi为时段 i 内比用水户 k 优先级高的用水户用过后剩下的水库群总入库迳流，且 Qi0。先由（2）式逐年逐日逆序倒推确定各用水户各年的预留库容线簇的外包线作为水库群各用水户的初始总预留库容线，再按各水库的兴利库容与水库群总兴利库容的比例将总预留库容线逐时段分配到各水库作为各用水户的初始供水调度线。令水库r 用水户 k 的最大预留库容Mrk=yrkUk (3) 式中：yrk为水库 r 用水户 k 的最大预留库容优化系数， Uk为用水户 k 的多年平均日用水量。通过对yrk寻优，寻优时用水户 k的供水调度线的各时段值随jrkjrkyy/1(j 为寻优迭代序号 )的比值同倍比放大或缩小，即可确定各水库、各用水户的最优供水调度线。2.4 目标函数和约束条件2.4.1 目标函数对于上述供水系统，其优化调度应以供水效益最大或以因缺水造成的损失最小作为目标函数。因供水效益涉及到产业结构等诸多因素而较难确定，一般采用以系统缺水损失最小作为目标函数。模拟权重系数模型的目标函数数学表达式可写为：min. E=nkmi11CkWik（4）式中：i 为以天为单位的计算时段序号， m 为计算系列的总天数， k 为用水户序号， n为用水户个数，Ck为用水户 k 的单位缺水量所造成的平均经济损失值，Wik为时段 i 用水户 k 的缺水量。且有：Wik=Gik-Dik=lr1(gikr-dikr)= lr 1(ri*Gik-dikr) （5）式中：Gik和 Dik分别为用水户 k 在时段 i 要求水库群供水的总量和水库群的实际供水总量，gikr和 dikr分别为用水户 k 在时段 i 要求水库 r 的供水量和水库 r 的实际供水量。因 dikr的值与时段 i 水库 r 的蓄水量 Vir、入库迳流 qir、水库的渗漏蒸发损失 eir以及时段 i 用水户 k 在水库 r 的预留库容 Pikr有关，而 Pikr又与 yrk有关，故 dikr=f(Vir,qir,eir,yrk) 。又因ri是r、xr、 Vir的函数，所以：Wik=F(Gik,Vir,qir,eir, r,Xr,yrk)。采用确定型迳流系列进行模拟顺序决策时，且eir显式处理为仅与 Vir有关，则 Gik,Vir,qir，eir, r均为已知，通过对 Xr、yrk寻优即可使目标函数达到最小。由于 Xr、yrk与时段 i 无关，模型属有约束非线性规划范畴。目标函数无解析表达式，宜采用直接法进行求解。若系统的目标函数以缺水量最小为目标，只需令各Ck=1即可。2.4.2 约束条件决策约束：决策变量为 xr和 yrk，两者均要求为非负数，即要求xr0、yrk0。边界条件：以丰水年汛期末为计算系列的起始点，以库满作为各供水工程（包括河网、湖泊、江道型水库及各水库）的初始条件。3 成果对比宁波市城市供水区属典型的多水源、 多用户、多级串并联的大型水资源系统， 系统概化图见图 1。该系统年际最大与最小年降水量比值达2.79，年内雨季三个月降水量占全年的67% ，农业用水占系统总用水量的 60% 以上。选用 1960-1980 年 21年系列，逐日进行供需平衡分析。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 4 页，共 8 页 - - - - - - - - - - 期刊文献用“权重模型”进行优化调配的2000年水平年推荐方案的xr、yrk的决策成果和最大预留库容精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 5 页，共 8 页 - - - - - - - - - - 期刊文献专用供水管道系统山区产水宁东片平原产水宁东片农业灌溉环境、航运群库片宁东水它山堰弃水抽水入库弃 水供 水河网东片宁山区产水宁西片平原产水宁西片群库水片西宁宁网河片 西厂宁片西水来自宁片东水厂自来工东业片一般宁东户用片确保宁户用保确片宁北确宁一般工业宁片北户用保片西般宁业工片西一厂宁片北来水自群丈库水片陆丈网河片陆州湾杭山区产水丈陆片群宁库水片北山区产水宁北片平原产水宁北片抽水入库弃 水供 水丈陆片平原产水供 水弃 水弃 水供 水宁网河片北水道江江姚抽 水供 水供 水图1 宁波市水资源系统概化网络图工业、畜牧业宁东片工业、畜牧业环境、航运农业灌溉片西宁农业灌溉环境、航运工业、畜牧业宁北片农业灌溉环境、航运工业、畜牧业丈陆片供水弃水抽 水弃 水库萧镇堰弃水供水供水供水供 水供水抽 水弃 水弃 水弃 水精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 6 页，共 8 页 - - - - - - - - - - 期刊文献值见表 1，计算系列的各年度缺水量成果见表2。用模拟法和多维增量动态规划法（MDDDP）相结合的“混合模型”进行优化调配的2000年水平年推荐方案的各年度缺水量成果4见表 2。应该指出，“权重模型”以不随来水和需水变化过程而变的Xr、yrk指导调度，影响了调度成果的优性。但由于水资源系统的多维性和多级串并联性，“混合模型”也存在离散精度问题影响了成果的优性。从表 2 可见，21 年系列的总缺水量“权重模型”比“混合模型”少4484万 m3；另外“权重模型”的缺水量较为分散，而“混合模型”的缺水量较为集中，缺水量集中将增加治理方案的难度。1967年至 1969年三年中，“权重模型”的缺水量为56904万 m3，“混合模型”的缺水量为64246万m3，“混合模型”的缺水量比“权重模型”多7342万 m3。4 结语本文提出的“权重模型”有较强的实用性。它特别适用于多级反馈、多级串并联的复杂供水系统，具有方法简便、建模迅速、计算工作量小等特点。能通过已发生的水文系列分析得出Xr*、yrk*，由 y*rk定量地确定各供水工程能适应各用水户有不同保证率要求的优化供水调度图，并结合 Xr*可定量地确定各供水工程当前时段的供水量，指导实时调度。表 1 2000年水平年推荐方案“权重模型”决策成果水库名称库容系数Xr* Yrn*(日)最大预留库容(万 m3) “确保”供水日最大供水量（万 m3/日）千亩岙8.53 10.284 14.25 870 25 亭下0.68 2.082 46.12 2750 20 横山0.41 1.500 横溪0.84 2.249 三溪浦0.54 1.936 19.29 1150 8.3 东钱湖0.59 1.988 新路岙0.50 1.894 宁东小(一)型0.75 2.155 周公宅0.84 2.249 164.35 9800 45 皎口0.67 2.072 宁西小(一)型0.55 1.946 十字路2.89 4.391 岗山2.06 3.524 16.77 1000 19 宁北小(一)型1.08 2.500 注：1、设置“确保”预留库容的水库由考虑专供管道投资等因素的方案比较确定。2、表中仅列了“确保”用水户的最大预留库容优化系数yrn*。参考文献1 M Heidart ,VT Chow et al ,Discrete Differential Dynamic Programming Approach to Water Resources Systems Optimization ,Water Resources Research , vol7，No2，April 1971. 2 YYHaimes , Hierarchical Analyses of Water Resources System ,1977. 3 张超等，系列模型方法在城市水资源开发利用规划中的应用，水利学报，1993年第 2 期。4Shen Peijun ，et al,Application of a Mixed Model to the Optimal Operation of Water Resources in the 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 7 页，共 8 页 - - - - - - - - - - 期刊文献Coastel River Network Area,Journal of Hydraulic Engineering,Vol.1,No.4,1992. 表 2 “权重模型”与“混合模型”优化成果对比权重模型混合模型日历年缺水量（万 m3）水文年缺水量（万 m3）1960 0 60.11-61.10 0 1961 0 61.11-62.10 0 1962 0 62.11-63.10 111 1963 1982 63.11-64.10 2298 1964 1409 64.11-65.10 0 1965 0 65.11-66.10 0 1966 0 66.11-67.10 20421 1967 27906 67.11-68.10 43745 1968 28986 68.11-69.10 80 1969 14 69.11-70.10 202 1970 0 70.11-71.10 1063 1971 3139 71.11-72.10 0 1972 0 72.11-73.10 0 1973 0 73.11-74.10 0 1974 0 74.11-75.10 0 1975 0 75.11-76.10 0 1976 0 76.11-77.10 0 1977 0 77.11-78.10 0 1978 0 78.11-79.10 0 1979 0 79.11-80.10 0 1980 0 80.11-60.10 0 合计63436 67920 A STUDY ON THE OPTIMAL OPERTION MODEL FOR A MUL TI WATER SOURCES，MUL TIUSERS，LARGE SCALE WATER RESOURCES SYSTEM Huang Fang 1Xu Wenbin2Zheng Jianqing3(1、Zhejiang Design Institute of Water Conservancy & Hydroelectric Power ，Hangzhou 310002；2、Zhejiang Provincial Bureau of Water Resources ，Hangzhou 310009；3、Wenzhou Municipal Bureau of Water Conservancy & Hydroelectric Power ，Wenzhou 325000 ）Abstract In this paper a simulation and distributive coefficient model is established for water resources system which is characterized by multi-water sources, multi-users, multi-step series-parallel connection. The simulation and distributive coefficient model puts forward a practical one by converting system sequence policy-making topic into nonlinear programming topic. Good advantages are clearly shown from simulation distributive coefficient model both in simple method and clear-cut results, which is contrasted with the results of mixed model including a simulation model and a multi-dimensional discrete differential dynamic programming model. Key Words water resources system, simulation and distributive coefficient model, reservoir regulation lines. 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 8 页，共 8 页 - - - - - - - - - -