一种基于优化的河流时间水权实时调度管理方法
【专利摘要】本发明涉及一种基于优化的河流时间水权实时调度管理方法，属于环境保护与资源综合利用【技术领域】。根据流域管理部门规定的水量分配目标，确定优化目标和约束条件，通过引水闸门不同时间的开闭定量控制河流时间水权的分配，建立基于河流时间水权的混合整数优化调度模型；确定优化的闸门开闭调度方案，进行闸门关闭的概率计算分析，确定闸门关闭的起始时间和时间长度，得到较优的河流时间水权实时调度管理方法。本方法明确了河流时间水权实时调度管理的数学量化方法，有利于市场经济下水资源的高效利用；根据优化模型能得到更合理的时间水权调度规则，保障最严格水资源管理制度的实施，具有理论意义明确、技术方法先进、便于实际应用的优点。
【专利说明】一种基于优化的河流时间水权实时调度管理方法
【技术领域】
[0001] 本发明涉及一种基于优化的河流时间水权实时调度管理方法，属于环境保护与资 源综合利用【技术领域】。
【背景技术】
[0002] 随着经济社会发展，水资源短缺问题日益严重，人类经济社会与生态环境用水的 矛盾也日益突出，水权制度逐渐成为维护用水秩序、促进市场经济调节水资源高效利用和 环境保护的重要途径。通过水资源使用权界定，不仅可以建立经济社会各部门之间的水资 源权利边界，而且可以较好的界定生态系统与人类社会之间的用水量边界，以保障生态用 水和经济社会可持续发展（郑航，2009)。水权有"水量水权"和"时间水权"两种方式，水量 水权是分配给各用水部门可使用的一定时期内（通常是一年）的总水量，时间水权是分配 给各用水部门在一定时期内（从天然河流中）可取水的时间。水量水权的实时调度管理中， 管理者和用水部门借助水利调蓄工程进行取用总水量的目标控制，手段明确、方法可靠。时 间水权的实时调度管理中，管理者往往是通过制止某一时段非授权用水部门的取水来保证 该时段授权用水部门的取水，实际操作是通过河流上取水闸门的开闭控制各用水部门的水 权使用，称为闭口。然而目前时间水权的管理仍是根据历史因素或经验获得（郑航，2009)， 缺乏科学方法的定量指导，在一些缺乏水利调蓄与控制工程的流域，难以适应当今水资源 三条红线的定量化管理的要求。
[0003] 综上所述，现行时间水权管理存在的问题是缺少明确的数学量化方法，难以协调 不同用水部门之间的用水水量；实时调度操作的计划没有明确合理的计算方法，难以制定 较优的闭口时间，产生了浪费与短缺并存，个别时段的基本经济社会和基本生态环境用水 难以保证等问题。
【发明内容】
[0004] 本发明的目的是提出一种基于优化的河流时间水权实时调度管理方法，结合河流 时间水权管理方式和数学优化模型工具，建立基于河流时间水权的0-1混合整数优化调度 模型，根据引水闸门不同时间的开闭定量控制和优化河流时间水权的调度管理。
[0005] 本发明提出的基于优化的河流时间水权实时调度管理方法，包括以下各步骤：
[0006] (1)根据河流管理部门的要求，确定中游供水区的需水量和下游供水区的 需水量WddCT1，将它们作为河流的供水目标，其中j为河流中游供水区的序号，j = 1，2,…， m，m为河流中游供水区总数，河流下游有一个供水区，每个供水区通过一条渠道从河流引 水，每个供水区通过一条分支回流至河道，中游的每个渠道入口设有闸门，通过河流中游和 下游的控制断面的水量全部供给下游供水区使用，下游供水区直接从河流引水（不受闸门 控制）；
[0007] (2)建立一个基于河流时间水权的0-1混合整数优化调度模型，优化模型的目标 函数为：
[0008]
[0009] 上式中：t为优化模型的时段指标，t = h、……tn，h为第一个时段，tn为模型计 算的最后时段，％为t时段河流中游第j个供水区的缺水量，w dladi(t)为t时段河流 下游供水区的缺水量为权重系数,通过随机或人工生成，满足Α+w；^ = 1 ;
[0010] (3)设定优化模型关于河流、河流中游供水区、河流下游供水区、汇流节点、渠道和 闸门的约束条件，如下：
[0011] 河流：
【权利要求】
1. 一种基于优化的河流时间水权实时调度管理方法，其特征在于该调度管理方法包括 以下各步骤： (1) 根据河流管理部门的要求，确定中游供水区的需水量和下游供水区的需水量 wdde;m,将它们作为河流的供水目标，其中j为河流中游供水区的序号,j = 1,2, -,111,111为河 流中游供水区总数，河流下游有一个供水区，每个供水区通过一条渠道从河流引水，每个供 水区通过一条分支回流至河道，中游的每个渠道入口设有闸门，通过河流中游和下游的控 制断面的水量全部供给下游供水区使用，下游供水区直接从河流引水（不受闸门控制）； (2) 建立一个基于河流时间水权的0-1混合整数优化调度模型，优化模型的目标函数 为：
上式中：t为优化模型的时段指标，t = t(1、……tn，t(1为第一个时段，tn为模型计算的 最后时段，为t时段河流中游第j个供水区的缺水量，wdladt(t)为t时段河流下游供 水区的缺水量为权重系数,通过随机或人工生成，满足Α+w；^ = 1 ; (3) 设定优化模型关于河流、河流中游供水区、河流下游供水区、汇流节点、渠道和闸门 的约束条件，如下： 河流：
河流中游供水区：
河流下游供水区：
汇流节点：
渠道：
闸门：
其中，『υ)、，u)、4、.⑴分别为第t时段向中游第j个供水区供 水的河流入流量、河流出流量，河流供水量和河流损失量，QHv max为河流的过流能力，Qd max 为河流的引水能力，和分别为第t时段河流 中游第j个供水区得到的供水量、耗水量、回归下游河流的水量、需水量和缺水量，wdsup(t)、 Wddem(t)、Wdlack(t)、Wdinc(t)分别为第t时段河流下游供水区的供水量、需水量、缺水量和力口 大量，ν__α)为第t时段汇流节点的出流量，表示第t时段汇流节点的第1个 分支入流量，N为流入同一汇流节点的分支总数，为第t时段向中游第j个供水区 供水的渠道出流量，a (t)为第t时段闸门开启或关闭状态，取值为0或1，β为渠道的损 失系数，取值范围为0-1的实数，上述约束条件中，河流的过流能力、河流的引水能力和渠 道的损失系数由河流管理部门提供； (4) 采用混合整数规划方法求解，以天为时间步长，根据30年水文历史资料中的河流 上游入流量WHv in(t)，求解上述步骤（2)和步骤（3)构成的混合整数优化调度模型，得到 优化的30年内第t时段河流中游第j个供水区的缺水量、河流下游供水区缺水量 Wdlac;k(t)和第t时段闸门开启或关闭状态a (t); (5) 随机或人工生成多组Wl、w2，得到权重集合，，重复步骤（4)，计算得到30年的多组 优化的第t时段河流中游第j个供水区的缺水量、河流下游供水区缺水量Wdladi(t) 和第t时段闸门开启或关闭状态a (t)，形成一个由多组优化的第t时段河流中游供水区的 缺水量、河流下游供水区的缺水量Wdladi(t)和闸门开启或关闭状态α⑴组成的非 劣解集
(6) 根据上述步骤（5)多组优化的第t时段河流中游第j个供水区的缺水量 和第t时段河流下游供水区的缺水量Wdladi(t)，通过下式计算得到河流中游供水区缺水率 Rm_lac；k和河流下游供水区缺水率R<uac；k :
当河流中游供水区缺水率Rm ladt与河流下游供水区缺水率Rd ladt的值相等时（或某一 个特定的协调值），从上述非劣解集W"中检索出与河流中游供水区缺水率Rm ladt和河流下 游供水区缺水率Rd_ladt相等条件下（或某一个特定的协调值）相对应的a (t); (7) 根据上述步骤（6)的a (t)值，利用概率统计方法，得到河流中游一年内每天闸门 关闭的分布概率？、？1 =叭1)/30，叭1)为30年中第1天闸门关闭的次数，1 = 1，2，···， 365 ； (8) 将上述一年内每天闸门关闭的概率pia分为四个调度时期计算每个调度时期内闸 门关闭的天数期望值τχ，作为该调度时期内河流下游供水区授权取水的时间长度：
其中，X = 1，2, 3, 4,分别代表四个调度时期，Mx为每个调度时期的天数； (9) 取每个调度时期内闸门关闭的起始时间Tx(l = tx(l，tx(l为每个调度时期的起始时间， 然后根据上述每个调度时期内闸门关闭的天数期望值T x，依次计算每个调度时期内闸门连 续Τχ天关闭的累积频率Ρχ，得到多个累积频率Ρ χ，
其中，Τχ(ι为满足tx(l彡Τχ(ι彡tx(l+M x-Tx中的正整数，从多个累积频率Ρχ检索出与最大累 积频率Ρχ _相对应的闸门关闭起始时间Τχ(ι ; (10) 根据河流管理部门提供的河流从上游流至下游所需的时间Λ t，以ΤΧ(Γ Λ t作为每 个调度时期闸门关闭的起始时间，得到以Tx(|-At、Tx为优化参数的河流时间水权实时调度 管理方法。
【文档编号】G06Q50/06GK104050537SQ201410284267
【公开日】2014年9月17日 申请日期:2014年6月23日 优先权日:2014年6月23日
【发明者】王忠静, 朱金峰, 郑航, 赵建世, 尚文绣, 黄草 申请人:清华大学