专利名称：一种非汛期水量协同优化调度方法
技术领域：
本发明属于水资源保护和水资源管理领域，涉及一种非汛期水量协同优化调度方法，具体地说是涉及一种供水型水库动态整体联合优化调度的方法，尤其是兼顾生态、生活、生产以及灌溉区的大中型供水型水库的动态整体优化调度方法。
背景技术：
水资源短缺、水环境日益恶化已成为制约社会经济发展的重要因素。以水库群为主的优化调度研究，是实现水资源调控的主要工程措施，意在解决非汛期各用水部门之间的矛盾，满足其基本要求。本发明通过群智能技术，分析各个环节水库群之间的相互作用， 协调统一河流流域各个水库，经济合理地利用各种水资源，包括地表水和地下水、农业用水和城市(工业、生活)用水、经济用水和生态用水等，解决供水与需水、用水和防污之间的水资源管理的矛盾，以获得综合利用的最大效益，解决实时动态统一水利调度问题。这项发明不同于单一地只考虑生态或发电量对水库水量分配调度的问题，比如以下发明专利“面向河流生态系统保护的发电型水库调度函数优化方法”(申请号 200910265477. 1)，考虑的是水库调度与河流生态、水库发电量之间最优化的关系。发明专利“水库预报调度的灰色动态规划方法”(申请号200710133917. 9)，主要研究的是水库调度分级与最优发电量调度路径问题。发明专利“兼顾人类和生态需求的供水型水库调度图优化确定方法”(申请号20091(^6M78. 6)，建立了最小生态需水和水库调度曲线相结合的水库调度规则，确定水库调度的生态效益目标和社会经济效益目标度量方法。发明专利“一种供水型水库生态调度技术”(申请号201010100909. 6)，以适宜生态需水贴近度表征水库调度的生态效益，以供水短缺率表征水库调度的社会经济效益，将生态效益和社会经济效益作为水库调度的目标进行优化调度，应用动态规划法对这些优化模型进行求解。以上这些专利都涉及到了水库优化调度问题，但是都没有具体研究非汛期流域水量统一调度问题；更加没有将“三生”用水(生产、生态、生活)与各级水库群作为一个统一调度整体进行研究。
发明内容
本发明的目的是提供一种非汛期水量协同优化调度方法，具体地说是提供一种供水型水库动态整体联合优化调度的方法，尤其是兼顾生态、生活、生产以及灌溉区的大中型供水型水库的动态整体优化调度方法，利用群智能最优化方法，解决多级流域水量协同优化调度，寻求理想的多水库群间调度方案，使生产、生态、生活、灌溉等各部门在整个分析期内的保证供水缺水量最小，满足最优水量分配。本发明根据流域水系、区域分布、水文水质监测站的数据记录，划定研究河段、设置控制断面，确定研究河段两岸的生活、生态、生产各方用水区域。通过研究控制断面间各大中型水库分布，依据水利规划、水资源规划、工农业发展布局，结合水系、地形、行政辖区及典型调查等情况综合分析，划分多级供水区域，结合 “支流水库并联化、干流水库串联化，，思想，建立水库群混联模型，应用群智能优化方法统一
6协同调控，均衡水量分配。本发明的一种非汛期水量协同优化调度方法，通过分析非汛期水利分布，构建网络概化图，然后以河流干流为主研究对象，获取各河流干流的最小需水量，在此基础上，将河道的水流量按照河流走向，由控制断面划分成各个河段，建立各河段最小需水流量模型， 采集各河段的区间入流量，计算出各控制断面最小需水流量，各断面最小需水流量等于各河段最小需水流量减去同河段的区间入流量，监测每个控制断面的径流量，当某控制断面的实时径流量小于该控制断面最小需水流量时，则启动预警，并采用协同粒子群智能优化算法，动态协同调控水库群，进行整体水量分配规划，并计算出各控制断面径流量，输出最终水量分配计划，获得非汛期流域水量调度方案，按该调度方案实施非汛期水量协同优化调度；其中，所述的网络概化图是指待研究流域的干流、支流、控制断面和水库群分布的网络概化图；具体步骤如下(1)根据水文站、水库坝址、水系分布及行政区界，标记待研究河段干流、支流，划分河段，确定控制断面；分析待研究河道的水库分布和管辖区域，抽象出待研究流域的干流、支流、控制断面和水库群分布的网络概化图；河段划分主要根据水文站、水库坝址、水资源分区及行政区界划定；河段水量调度监测点则以河流上的水文站或大型水库坝址标注，又称为控制断面，其中水文站是水量调度的断面，可进行流量监测；大型水库坝址断面主要控制水库下泄流量；(2)以河流干流为主研究对象，获取河道的产汇流、河道外用水量、河道的工、农、 林、生活回水量以及河道内最小需水量，具体如下(I)控制断面间河道的径流量，又称产汇流降雨到径流过程受到许多因素的影响与制约，人们一般把降雨形成流域出口断面的径流过程概化为两个阶段产流和汇流； 降雨转化为净雨的过程为产流过程；净雨沿地面和地下汇入河网，并经河网汇集形成区间流域出口断面的径流过程，称为区间流域汇流过程，两者相应的计算合称流域产汇流计算。 非汛期，表土的下渗能力强，流域内壤中流和地下径流丰富，适宜用产汇流模型进行产汇流计算。产汇流模型涉及蒸散发计算、产流计算、分水源计算和汇流计算等，模拟得到单元面积到达流域出口的出流过程；将各个单位的出流过程线性叠加，得到区间流域总的出流过程。(II)河道外最小需水量生活用水主要包含居民用水、牲畜用水两部分，同时需考虑城镇居民、农村居民用水额定需求不同，以及农村大、小牲畜用水额定需求差异。生活用水量=城镇居民数X 城镇居民用水定额+农村居民数X农村居民用水定额。工业用水定额用万元总产值取水量和工业用水重复利用率进行衡量，其值的大小与水资源条件、工业产业结构、工艺水平等多种因素有关。工业用水量=万元总产值取水量χ总产值。农业用水定额即各种作物用水定额之和。各种作物用水量等于播种面积X灌溉定额。其中作物播种面积=种植结构系数χ播种面积。河道外生态环境需水量包括每一个河段采取人工措施的水土保持治理生态环境需水；每一个河段提高城市公共绿化覆盖率生态环境需水，可根据以往经验假定每年整体用水量保持某个小的系数增长。林牧渔用水林牧渔用水所占比例小，用水定额假定不变，林牧渔用水为林牧渔规模X林牧渔用水定额(每亩每年的用水定额预测)。(III)工业、农业、林牧渔、生活回水量各个控制断面区域的工业、农业、林牧渔、生活回水量与其用水量成一定的比例。 于是，回水量=各用水量χ各回水系数。(IV)河道内最小需水量包括河道基本生态需水量和河道内环境需水(河道冲淤输沙需水量和河道自净、稀释需水量)，以及航运需水量。对通航有要求的河道，水库设计死水位时已考虑了满足河道通航需求，故计算河道内需水量时可不再单独考虑航运需水量。 在这里，取河道基本生态需水流量和河道内环境需水流量的最大值，得到河道内最小需水流量。a.河道内基本生态需水模型河道内生态需水通常是指河流为了维持某一特定生态系统的基本生态功能，河道应保持的流量。设河道内基本生态需水为Wb，综合90%保证率最枯月平均流量和河口生态需水(7Q10法)确定。《国家制定地方水污染物排放标准的技术原则和方法》(GB3838-83) 中规定一般河流采用近10年最枯月平均流量或90%保证率最枯月平均流量作为河流最小流量。鉴于我国水环境现状，对于河口生态需水，90%保证率最枯月平均流量计算要求过高，一般会采用7Q10法完成河口生态需水量计算。b.河道内环境需水模型河道内环境需水包括河道冲淤输沙需水量和河道自净、 稀释需水量(i)设河道冲於输沙需水为Ws，采用汛期最小输沙量法得到 Ws=^
权利要求
1. 一种非汛期水量协同优化调度方法，其特征是由两部分构成第一部分获取研究流域的河流、水库整体流量方法；主要包括通过分析非汛期水利分布，构建网络概化图；然后以网络概化图中河流干流为主研究对象，考虑生态、生活、工业、农业各方的用水和回水信息，获取各河流干流的最小需水量；在此基础上，将河道的水流量按照河流走向，由控制断面划分成各个河段，建立各河段最小需水流量模型，采集各河段的区间入流量，计算出各控制断面最小需水流量；第二部分在获取各控制断面最小需水流量之后，实时监控控制断面的径流量，当某控制断面的实时径流量小于该控制断面最小需水流量时，则启动预警，并采用协同粒子群智能优化算法，动态协同调控流域水库群，进行整体水量分配规划，并计算出水量分配规划后各控制断面径流量，输出最终水库水量分配计划，获得非汛期流域水量调度方案；其中，所述的网络概化图是指待研究流域的干流、支流、控制断面和水库群分布的网络概化图；具体步骤如下(1)根据水文站、水库坝址、水系分布及行政区界，标记待研究河段干流、支流，划分河段，确定控制断面；分析待研究河道的水库分布和管辖区域，抽象出待研究流域的干流、支流、控制断面和水库群分布的网络概化图；(2)以河流干流为主研究对象，获取河道的产汇流、河道外用水量、河道的工、农、林、生活回水量以及河道内最小需水量；(3)建立河段最小需水流量模型，将河道的水流量按照河流走向，由控制断面划分成各个河段；河道最小需水流量等于河道各河段的产汇流和+各种回水量_河道外需水量_河道内最小需水量；将河道最小需水流量根据以往各个河段分配系数来确定各河段最小需水流量；(4)获取河段的区间入流量区间入流量有别于工、农、林、生活回水量，它主要指一些小支流的流入干流的流量；计算河段区间入流量时，考虑到缺水期小支流基本断流，故只考虑流域河段上集水面积大于 500km2的支流，对于集水面积500km2以下的支流，一般不计小支流的入流量；1)对于集水面积大于2000km2的支流，由于已计算了河道外及河道内需水流量，故支流流入干流河段的区间入流量，取其支流最下游一个河段的最小流量减去同一河段的河段外耗水量；2)对于集雨面积小于2000km2的支流，由于河段相对于其他河段来说，社会经济发展速度相对较慢，河道外耗水也较小，因此计算区间入流量时将忽略河道外耗水量；此时，河道内用水主要是生态需水量；因此，对于集水面积小于2000km2的支流河段，流入干流的区间入流量将直接取支流河段的生态流量；(5)计算控制断面最小需水流量为确定水量调度启动条件与方法，需将前述计算结果的河段最小需水流量换算成该河段控制断面的最小需水流量，以便进行监测；由于以河段上断面作为本河段的控制断面，因此各河段控制断面最小需水流量等于各河段最小需水流量减去同河段的区间入流量；(6)监测每个控制断面的径流量，分析干流、支流各个控制断面区间区域的来水量和回水量，如果控制断面的实时径流量小于控制断面最小需水流量，则启动预警；(7)预警启动后， 采用协同粒子群智能优化算法，动态协同调控水库群，进行整体水量分配，其详细过程如下(I)在一个D维空间中，设定种群规模为m、迭代次数为iter、加速因子C1= C2 = 2.0， 权重因子w;针对水量调度问题，设置为粒子位置矢量，其各个位置矢量为水库对应的出库流量为< (粒子i情况下j库t时段出库流量)，0广表示如下οτα,……，……，d..4...4，《，…《…，‘，…，‘乂，…，⑦(1)其粒子位置是由速度向量if"来改变的，if"定义如下Viter =Cv1 V2 ...... VD V1 V2 ---V ··· VD ··· V1 V2 ... VD)⑵yiV^l, vi\->, vH , vij, vij,vIJ, vIJ , , vik, vik, , vik)其中<分别为粒子i情况下j库t时段水库的库容。上述及以下叙述中，i代表粒子，i = 1，2···，πι，； j代表水库，j = 1,2,…，k，k为水库总数量；t代表时段，t=l，2，…，D，D为总调度期。随机初始化m个粒子，所有的粒子组成一个特征空间，每个粒子的位置向量都是解空间中一个潜在的可行解；在每一次迭代中，粒子通过两个“极值”来更新自己的位置，一个是粒子本身的最好解， 称为个体极值点Pbest，另一个是整个种群目前的最好解，称为全局最优解，初始化pbest =gbest = 0 ；(II)计算各粒子的适应度值； 特征选择算法的适应度函数设计为Fit{g;er ) = max(^+^2-D1)(3)其中=R1表示每个控制断面径流量；R2表示水库下泄流量转换到控制断面的调节水流量；若以月作为调度，D1表示每个控制断面月均最小需水流量； 其约束条件如下水量平衡约束=vT1+^ ——乓(4)水库库容约束出库流量约束^Qj^(6)非负约束<，ν；1 , ^>0(7)式(4)中V^1、<分别为j库t时段初、末库容《为j库t时段的天然入库流量；h为水量与流量之间的转换系数，h = 2. 592 X IOV1 ·Α〗为j库t时段的损失水量4为j库t时段出库流量；式(5)中。》、Vt分别为j库t时段允许的最小、最大库容；式(6)中f .，j max^iZmm0;max分别为j库t时段最小、最大下泄流量能力。0丨_一般取枯水期调节流量，0丨max分别为 t时刻j库的库容量对应的最大下泄流量；达max求解方法利用水库的库容量与水位，水位与下泄流量对应关系表(水库设计参数)，确定水库的库容量与下泄流量(下泄能力)的函数关系；根据水库的t时刻的库容量，来确定水库最大的下泄能力；(III)更新全局极值gbest；比较各个粒子适应度大小，依据下式(8)挑选出适应度值F最大的粒子，并记录下其序号，这个新粒子适应度是否优于全局极值gbest，如果优于则更新全局极值gbest，i = i+1 ；F(Qi)=HiaxFit (8)(IV)判断是否满足终止条件终止条件为迭代次数i已达到最大迭代次数iter，或者 I Fit(Qr)-FitiQr)) \<eps，其中印s为非常小的正数。如果满足终止条件，输出，gbest ； 否则转向下一步；(V)i是否恒等于m余3取整的值当条件不成立时，判断新粒子适应度是否优于Pbest适应度，如果优于则将pbest设为新位置，并微调w，更新粒子的速度与位置v;x=wvr+《r(P:t-an+c^^z-on⑷Q-" =Q-er +V；^(10)当条件成立时，将全局极值对应的粒子作为父代，按照变异操作速度向量^r中各个分量的 <<巧·，^:]，有VU=aQ)mm+{l-a)Q)^(11)Vtjmwi= PQtjmmHl-(12) α = (l+rand())/2 (13) β = (l-rand())/2 (14)randO为
的一个随机数，另夕卜，若《>V;-时，K =^max，当《<"mm时，VhVU .依据式(9)，(10)重新更新粒子i的位置和速度； ，(VI)转向第(II)步循环，直到跳出循环输出gbeSt，Q^，获得年调度期各个水库水量分配的最优调度方案；(8)依据年调度期内各个水库水量分配的调度方案，计算出各控制断面径流量；依据动态水库群调控，将每个水库调节水量依据分配系数，将水量分配到各个控制断面，其中分配系数由各个水库的控制范围以及控制的流域面积来决定由于水库只能对位于控制河段下游的控制断面起到调节作用，因此分别考虑每个水库的调控作用时，对于无法调节的断面，其分配系数设置为0 ；对于能够调节的断面，其分配系数设置为该断面控制的流域面积与所有能够调节的断面流域面积和之比；(9)输出最终水库水量分配计划和各控制断面径流量，获得非汛期流域水量调度方案， 按该调度方案实施非汛期水量协同优化调度。
2.根据权利要求1所述的一种非汛期水量协同优化调度方法，其特征在于，所述的河段划分主要根据水文站、水库坝址、水资源分区及行政区界划定。
3.根据权利要求1所述的一种非汛期水量协同优化调度方法，其特征在于，所述的非汛期水量预警依据控制断面的缺水情况为三个等级1级预警，即黄色预警，表示轻度干旱，；II级预警，即橙色预警，表示中度干旱；III级预警，即红色预警，表示严重干旱；当控制断面来水量等于或小于1. 1倍断面最小流量时，发布黄色预警；当控制断面来水量等于或小于1. 1倍断面最小流量与河道外生产性需水流量之差时，发布橙色预警；当控制断面来水量等于1.1倍断面基本 生态环境流量，则发布红色预警。
全文摘要
本发明涉及一种非汛期水量协同优化调度方法，先从河流流域的生态、生活、工业、农业等整体上考虑各方的用水和回水信息、并通过计算河道产汇流、河道外需水量以及河道内最小需水量，确定河段和控制断面最小需水流量。通过对流域河道进行实时水量监控，当控制断面流量低于其最小需水流量时，自动输出缺水预警等级；最后，在不同预警等级下，通过协同粒子群智能优化算法，动态整体调控水库群，达到水量统一调度功能。此方法充分考虑了生态、生活、工业、农业的用水，依据非汛期“先确保生活，后满足农工”的原则，确保流域各个用水部门用水的均衡。
文档编号G06Q10/00GK102156914SQ20111007825
公开日2011年8月17日 申请日期2011年3月30日 优先权日2011年3月30日
发明者丁永生, 李荣昉, 杨永生, 王伟, 程丽俊 申请人:东华大学