考虑城市-水库整体风险的汛限水位调整方法研究.pdf

收稿日期：基金项目：国 家 自 然 科 学 基 金 资 助 项 目（，）；河南省高校科技创新人才支持计划项目（）；河南省青年人才托举工程项目（）；国家大坝安全工程技术研究中心开放基金资助项目（）作者简介：史金坤（），男，河南周口人，硕士研究生，研究方向为水工结构 通信作者：王建有（），男，河南洛阳人，正高级工程师，博士，研究方向为水工建筑物的安全性检测与评价：【防洪治河】考虑城市水库整体风险的汛限水位调整方法研究史金坤，段志昌，王建有，葛 巍（郑州大学 黄河实验室，河南 郑州；中国建筑第五工程局有限公司，湖南 长沙）摘 要：针对城市化发展造成水库职能变化的背景下水库汛限水位设置不合理的问题，通过建立通用的优选指标体系，采用基于熵权细化主观权重的指标权重确定方法，构建了综合考虑水库防洪风险和城市供水风险的汛限水位调整方案模糊优选模型。该模型在考虑决策者不同偏好的基础上，根据模糊理论中最大隶属度原则进行方案优属度的计算，从而确定最佳调整方案，并以郑州市尖岗水库为例进行了汛限水位的优化调整。结果表明：在拟定的 个汛限水位方案中，当决策者对防洪目标、供水目标具有不同的偏好时，如两者重要度之比为 、，则最优汛限水位方案分别为、。关键词：汛限水位；模糊优选；风险率；尖岗水库中图分类号：文献标志码：引用格式：史金坤，段志昌，温亚霖，等考虑城市水库整体风险的汛限水位调整方法研究人民黄河，（）：，（，；，）：，：；引言随着城市化的发展，城市区域进一步扩大，使得越来越多的水库进入了城区，水库的防洪安全对于城市来说越来越重要。另外，对于那些缺水型城市而言，其上游的水库一般要承担向城市供水的重任，城市的参与使得水库的职能发生变动。汛限水位是水库在汛期允许兴利蓄水的上限水位，是防洪调度的起调水位，是协调水库运行管理中防洪与兴利矛盾的关键指标。对于此类水库来说，若仍采用其原设计汛限水位，显然是不合适的，需要对其设计汛限水位进行调整。汛限水位的调整会影响水库的防洪风险和城市供水风险，较高的汛限水位可以减小城市供水风险，但会增加水库的防洪风险。相反，较低的汛限水位能减小水库的防洪风险，但会增加城市供水风险。因此，需要在考虑城市和水库整体风险的基础上来调整汛限水位。进行风险效益分析可以更科学直观地体现汛限水位调整所带来的利弊效应，为方案的决策提供依据，已第 卷第 期 人 民 黄 河 ，年 月 ，有的关于汛限水位调整的研究大都集中在风险效益分析和水位控制两方面，其中水位控制聚焦在分期控制和动态控制方面。针对水库职能变化背景下汛限水位的调整也有学者进行了初步尝试。一方面，依据对流域的气候、气象规律的研究，结合降雨预报技术来执行动态变化的汛限水位，但很少考虑水库职能的变化，即使考虑职能变化，也仅仅是对某一特定水库汛限水位调整的必要性进行分析。另一方面，在优选指标的选取上，汛限水位调整所带来的水库供水能力的变化常用供水保证率、可靠度等指标表示，当水库作为备用水源时，缺乏相应的指标表征。鉴于此，本文针对此类水库汛限水位调整方案的决策建立通用的优选指标体系，同时为克服主观权重带来的较大分析负担和客观权重易造成较大误差的不足，采用基于熵权细化主观权重的方法来计算指标权重，建立考虑城市水库整体风险的汛限水位调整方案模糊优选模型，以期为水库汛限水位的优化调整提供一种新思路。考虑城市水库整体风险的汛限水位调整方案优选模型构建 优选指标的确定汛限水位调整方案决策涉及的指标较多，很多指标之间有互相竞争甚至互相矛盾的关系，本文根据研究对象的特性将指标分为两部分，即水库自身风险指标和下游城市风险指标。水库自身风险主要为水库自身的防洪风险，下游城市风险包括下游城市防洪风险和城市供水风险，在防洪风险指标方面，参考已有研究成果，选取较为常用的防洪风险率和易损度指标；在供水风险指标方面，引入适用性更广的供水风险率指标，由此建立汛限水位调整方案优选指标体系，见图。图 汛限水位调整方案优选指标体系水库自身极限防洪风险率（自）常用于水库自身防洪风险的定量分析，即库水位超过水库校核洪水位（校）的概率。通过对各起调水位调洪演算结合频率分析法就可求出不同汛限水位下的水库自身极限防洪风险率 自。自 （校）（）式中：为洪水调节最高水位，为频率。下游城市防洪风险率（下）常用于下游城市防洪风险的定量分析，即水库下泄流量超过下游河道安全流量的概率。一般情况下，当库水位超过防洪高水位时水库下泄流量将超过下游河道安全流量，这样下游城市的防洪风险率就可以用库水位超过防洪高水位的概率来表示。下 （安或 防）（）式中：为水库下泄流量，安为下游河道安全流量；为库水位，防为防洪高水位。易损度表示大坝和下游防护对象遭遇险情的严重程度，对于大坝，当总库容一定时，库水位越高，溃坝造成的后果越严重；对于下游防护对象，水库泄流量越大，可能造成的淹没损失越大。在此采用调洪最高水位超过 校的幅度表示大坝易损度（坝），调洪最高水位超过 防的幅度表示城市易损度（城），其中大坝易损度的最大值为坝顶高程与校核洪水位的高差，若调洪最高水位高于坝顶高程，其大坝易损度为此最大值。抬升汛限水位会增加汛期水库的蓄水量，减少汛期弃水，使得水库在汛期结束后更容易蓄较多的水。因此，可以认为水库在一整年内的存水越多，城市供水风险就越小；相反，存水越少，城市供水风险就越大。不同水库供水形式有所不同（如持续供水、应急供水等），其供水风险的表征就会不同，在此引入一种通用的供水风险表示方法，即假定水库一整年内能蓄到满足城市用水的水量时，此汛限水位下城市就无供水风险。具体为：首先根据来水与用水情况计算水库的最大蓄水量（），见式（），然后可用供水风险率（水）来表示不同汛限水位下的城市供水风险，此指标也可以通过转换来反映不同汛限水位的缺水量，计算公式见式（）。来 用（）式中：来为水库来水量（如非汛期上游来水、调水工程引入水等），用为除供水用途（城市生活用水、紧急备用水等）外的用水量（如生态用水、蒸发渗漏损失水等）。人 民 黄 河 年第 期水满 （满）（满）（）式中：满为满足城市用水需求的水库库容，为满与原设计汛限水位对应库容的库容差。基于熵权细化主观权重的优选指标权重确定本文采用基于熵权法细化主观权重的方法来确定指标权重，首先根据专家的意见确定防洪和城市供水两个目标的相对重要程度，得出两者经验权重，然后在每个目标下通过熵权法分别计算属于该目标的下一级指标权重，进而得出所有指标的权重值。该方法既充分体现了决策者的主观偏好，又不增加决策分析者的负担，同时体现出熵权法依据客观数据赋权的客观性优势，该方法可以更科学地对指标进行赋权。基于模糊优选模型的方案优选模糊优选模型是在模糊数学模糊集理论的基础上发展而来的，其可以充分考虑系统模糊性、不确定性及复杂性，通过隶属函数把目标体系和最优决策方案联系起来，利用模糊优选公式计算相对隶属度，根据定量化的相对隶属度值来评价方案的优劣。由于不同的汛限水位方案对防洪和城市供水的侧重有所不同，最优方案的确定是一个多目标多层次的复杂决策问题，在决策过程中往往存在重要影响因素无法量化的问题，因此模糊优选模型对其有较好的适用性。具体步骤如下。）建立优选指标相对隶属度矩阵。设待优选的汛限水位共有 个方案，每个方案有 个指标值，则有方案集的指标特征值矩阵 ：（）（）式中：为方案 的第 个指标值（，；，）。各指标的量纲不同，需对指标进行归一化处理（本文建立的指标均属越小越优型），得到相对隶属度矩阵 （）。（越小越优型）（）式中：为方案 指标 的隶属度，、分别为方案集中指标 的最大、最小值。）计算优选指标的权重。根据 节可以得到各指标的权重值，从而得出指标权重向量 （）。）计算各方案的优属度：（）（）（）式中：为距离参数，为海明距离，为欧氏距离，本文取 ，即采用欧氏距离进行计算。计算结果表示各方案对于优等方案的隶属度，即各方案的优属度，优属度最大的方案即为最佳方案。实例应用 水库概况尖岗水库位于淮河水系颍河支流贾鲁河干流上游的尖岗村西，距郑州市区 ，水库控制流域面积。由于该水库居于郑州市区上游，下游有多条铁路、公路，还有水厂、电厂等供水及能源设施，还是人口、经济的高度聚集区，防洪位置十分重要，因此按大型水库进行管理。水库现状正常蓄水位为 ，总库容 万；设计洪水标准为 一遇，校核洪水标准为 一遇安全保证值（折合洪水频率为）；汛限水位 ，设计洪水位，校核洪水位 。根据河南省防办编制的河南省大型及重点中型水库汛期调度运用计划，水库主汛期为 月 日 月 日，库水位低于溢洪道底部高程时采用泄洪洞泄洪，库水位高于溢洪道底部高程时由溢洪道自由泄洪。尖岗水库续建时规划该水库功能是以防洪、城市供水、农田灌溉为主，兼顾发电、养鱼等综合利用，由于上游来水不足，失去了发电和灌溉等功能，随着黄河水的引入、南水北调中线工程的通水以及郑州市社会经济的快速发展，水库的主要职能变为防洪和城市紧急供水，而汛限水位多年未变，水库在汛期结束后往往无水可蓄，不利于城市紧急供水。尖岗水库汛限水位调整方案模糊优选模型 防洪风险率与易损度指标考虑到尖岗水库原有特征水位和实际情况，本文拟定、等 个汛限水位方案。根据尖岗水库历年实测洪水资料，将典型洪水按洪峰、洪量同频率进行放大，推求不同频率的设计洪水过程，此处频率 的取值分别为（校核洪水）、。然后以各方案为起调水位，按规划调度方式调节各频率的洪水，可得各方案调洪高水位，见表。人 民 黄 河 年第 期表 不同汛限水位下各频率洪水的调洪高水位方案汛限水位不同频率洪水的调洪高水位 （）（）（）（）（）（）尖岗水库下游贾鲁河经过综合治理后最大防洪标准达到 一遇，因此防洪高水位与设计洪水位相同。根据 中的方法可计算出水库自身极限防洪风险率（自）、下游城市防洪风险率（下）、大坝易损度（坝）和城市易损度（城），结果见表。从表 中可以看出，随着汛限水位的抬升，防洪风险率和大坝、城市易损度都逐渐增大。表 不同汛限水位下的防洪风险率和易损度方案汛限水位 自 下 坝 城（）（）（）（）（）（）供水风险率指标尖岗水库正常蓄水位对应的蓄水量可满足城市紧急用水需求，因此 满为正常蓄水位对应的库容。汛期水位一般都维持在汛限水位，因此尖岗水库最大蓄水量计算中来水主要考虑非汛期的上游来水，此处月增蓄水量的计算中去除 月的数据（见表）。不同水平年上游来水量不同，可选取年径流量与多年平均值接近的平水年进行计算，根据尖岗水库设计平水年的月均流量值（见表），扣除蒸发、渗漏损失和生态用水，损失系数取，可得出水库的最大蓄水量，由式（）可得各汛限水位方案的城市供水风险率（水），结果见表。从表 中可以看出，随着汛限水位的抬升，水库非汛期可蓄到的最大水量逐渐增大，水逐渐减小。表 尖岗水库设计平水年的月均流量及增蓄水量月份平均流量（）增蓄水量 万 表 不同汛限水位下的城市供水风险率方案汛限水位 库容 万 万 水（）（）（）（）（）（）从风险率的角度看，随着汛限水位的抬升，防洪风险逐渐增大，城市供水风险逐渐减小（见图），此处也明显体现出汛限水位调整中水库防洪任务与供水任务的竞争性。图 不同汛限水位下风险率的变化 最优汛限水位方案的确定根据图 中的指标体系和表 调洪演算结果可求出各指标值（见表），则可以建立优选方案的特征值矩阵 ，根据式（）可得特征值的相对隶属度矩阵 。表 尖岗水库汛限水位调整方案优选指标体系汛限水位 自 坝 下 城 水 人 民 黄 河 年第 期 汛限水位调整方案优选指标可以分为两类，一类是防洪目标（自、下、坝和城），另一类是供水目标（水）。当决策者目标重要度偏向防洪目标时，如两者权重之比为 ，在防洪目标下属指标中根据熵权法计算熵权，进而得出最终参与优选的指标权重，结果见表，因此权重向量 （，），模型参数建立完备。表 熵权法细化主观权重计算结果指标熵值熵权计算权重自 坝 下 城 根据式（）可计算出各方案的优属度，结果为 （，），方案（）为最优方案，即最优汛限水位为 。同理，当决策者目标重要度偏向供水目标时，如两者权重之比为 ，此时权重向量 （，），（，），方案（）为最优方案，即最优汛限水位为 。指标特征值和方案优属度计算结果表明，在水库职能变化后综合考虑防洪和供水目标的情况下，相较于原设计汛限水位，在偏好防洪目标或供水目标时，汛限水位分别提高了 和，分别为 和，水库自身极限防洪风险率由 增大到和，下游城市防洪风险率由 增大到和，供水风险率由 减小到 和。在汛期的不同阶段，往往会有不同的目标偏好，决策者越偏好防洪目标最优汛限水位越低，越偏好供水目标最优汛限水位越高，该模型可在决策者不同偏好的基础上进行最优方案的决策。当防洪与供水目标的权重较为合理时，通过模糊优选模型所确定的最优方案最接近理想最优方案，为最符合工程实际状况的方案，因此建议在确定防洪与供水目标相对重要度时，应咨询领域内相关专家，以确定合理的指标权重，从而确定更加合理的汛限水位调整方案。结束语本文综合考虑城市水库整体风险，提出了城市化发展造成水库职能变化背景下水库汛限水位调整的方法，针对水库供水形式的不同，引入通用的城市供水风险表示方法，并应用数学计算方法进行最优汛限水位方案的决策，为该类水库汛限水位的优化调整提供了新的思路。本文只对水库矛盾突出的防洪风险和城市供水风险进行了研究，下一步可以融合发电、航运、灌溉等多个调度目标建立内容和适用性更加广泛的汛限水位控制模型。参考文献：李晓英，张琰，童泽淳水库汛限水位控制多目标协同决策研究水力发电学报，（）：冯平，韩松，李健水库调整汛限水位的风险效益综合分析水利学报，（）：丁伟，周惠成水库汛限水位动态控制研究进展与发展趋势中国防汛抗旱，（）：刘菁稳，玄英姬，陈星汛限水位分期控制方法在石岩水库的应用人民黄河，（）：，（）：秦雪，马长明，陈鹏宇，等岩滩水库汛限水位动态控制域的应用研究人民黄河，（）：张发荣，惠磊双塔水库供水状况及汛限水位调整必要性分析农业科技与信息，（）：，张艳平，张双虎，周惠成，等水库汛限水位实时动态控制方案优选研究中国水利水电科学研究院学报，（）：，吴泽宁，胡彩虹，王宝玉，等黄河中下游水库汛限水位与防洪体系风险分析水利学报，（）：郭军峰，冯平流域下垫面变化对水库防洪风险率的影响水力发电学报，（）：周惠成，张改红，王本德，等基于防洪预报调度调整水库汛限水位的研究水力发电，（）：田林钢，马成功，王绪基于 熵权法的水利工程业主风险模糊评价 人民黄河，（）：，陈守煜系统模糊决策理论与应用大连：大连理工大学出版社，：潘杰义，刘西林科研项目评价指标体系及模糊优选决策模型研究科学学与科学技术管理，（）：于黎雅贾鲁河综合治理工程防洪影响分析河南水利与南水北调，（）：张莉郑州市生态水系水质水量联合调度研究与应用郑州：华北水利水电大学，：【责任编辑 许立新】人 民 黄 河 年第 期