CN107704965A - 一种兼顾生态保护需求的水库兴利调度图优化方法 - Google Patents
一种兼顾生态保护需求的水库兴利调度图优化方法 Download PDFInfo
- Publication number
- CN107704965A CN107704965A CN201710963168.6A CN201710963168A CN107704965A CN 107704965 A CN107704965 A CN 107704965A CN 201710963168 A CN201710963168 A CN 201710963168A CN 107704965 A CN107704965 A CN 107704965A
- Authority
- CN
- China
- Prior art keywords
- mrow
- reservoir
- ecological
- msub
- scheduling
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000005457 optimization Methods 0.000 title claims abstract description 40
- 230000008901 benefit Effects 0.000 claims abstract description 39
- 230000000452 restraining effect Effects 0.000 claims abstract description 32
- 241000251468 Actinopterygii Species 0.000 claims abstract description 10
- 230000005514 two-phase flow Effects 0.000 claims abstract description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 66
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000005065 mining Methods 0.000 claims description 12
- 230000005611 electricity Effects 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 5
- 238000007599 discharging Methods 0.000 claims description 3
- ZZUFCTLCJUWOSV-UHFFFAOYSA-N furosemide Chemical compound C1=C(Cl)C(S(=O)(=O)N)=CC(C(O)=O)=C1NCC1=CC=CO1 ZZUFCTLCJUWOSV-UHFFFAOYSA-N 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 241000196324 Embryophyta Species 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 244000140975 Prunus simonii Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000008400 supply water Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种兼顾生态保护需求的水库兴利调度图优化方法,包括获取水电站水库的基本信息数据;选取代表鱼类,采用物理栖息地模型PHABSIM建立下泄流量与生态保护程度之间的关系,确定生态保护目标对应的生态流量区间;基于发电效益和生态效益建立水库调度图优化模型;以入库径流为水库调度图优化模型的输入,采用DPSA方法求解模型;绘制生态限制线,建立生态限制线和水库调度曲线相结合的水库调度规则,编制多目标优化调度图。本发明根据鱼类不同的生长周期采用不同的流量区间对生态目标进行优化,更符合实际情况;并通过增加生态限制线,来减小生态流量破坏深度,进一步加强对下游生态的保护。
Description
技术领域
本发明涉及水库调度图优化方法,特别是涉及一种兼顾生态保护需求的水库兴利调度图优化方法。
背景技术
水库调度图使用操作简便被广泛的运用于水库调度中,制定合理的水库调度图可以提高水库的社会经济效益并能有效的保护河流生态系统。
目前在调度图优化方面的研究主要考虑发电和供水等经济效益目标,对下游生态环境的保护的考虑较少。
现有的考虑生态的水库调度图通常是选取一个固定的流量作为生态流量进行调度图的优化,不能反映年内生态流量需求的变化,不能有效的保护下游的生态环境。即现有方法难以实现经济效益最大化的同时有效提高河流生态系统保护。
发明内容
发明目的:为解决现有技术的不足,提供一种能保证水库的运行不会造成河流生态系统的严重破坏,使水库运行带来的经济和生态的综合效益最大化,即兼顾生态保护需求的水库兴利调度图优化方法。
技术方案:一种兼顾生态保护需求的水库兴利调度图优化方法,包括以下步骤:
(1)获取水电站水库的基本信息数据;
(2)选取代表鱼类,采用物理栖息地模型PHABSIM建立下泄流量与生态保护程度之间的关系,确定生态保护目标对应的生态流量区间;
(3)基于发电效益和生态效益建立水库调度图优化模型;
(4)以入库径流为水库调度图优化模型的输入,采用DPSA方法求解模型;
(5)绘制生态限制线,建立生态限制线和水库调度曲线相结合的水库调度规则,编制多目标优化调度图。
进一步的,所述步骤(2)中采用物理栖息地模型PHABSIM建立下泄流量与生态保护程度之间的关系,其中,水库下泄流量对应的代表鱼类加权可利用面积的计算公式为:
Ci=f(vi)×f(di)×f(ci) (2)
式中,Ai为研究河段第i个网格的水域面积;Ci为第i个网格的组合适应度因子;f(vi),f(di),f(ci)分别为第i个网格的流速、水深、河床地质适宜度值;根据流量条件计算出各网格不同的流速、水深,由式(1)和式(2)进一步计算出各流量条件下的WUA值,得到流量-WUA关系曲线;
t时段水库下游的生态保护程度Rt的计算公式为:
式中,Rt为t时段水库下游的生态保护程度;WUAt为t时段水库下泄流量对应的代表鱼类加权可利用面积,WUAt,max为t时段对应最大加权可利用面积,建立流量与生态保护程度之间的关系。
进一步的,所述步骤(3)中建立水库调度图优化模型为:
maxE(x)={E1(x),E2(x)}
(4)
s.t x∈G(x)
其中,E(x)表示目标函数集,E1(x)表示生态效益目标函数,E2(x)表示发电效益目标函数,G(x)表示约束条件集。
更进一步的,所述建立水库调度图优化模型的方法包括:
(31)建立生态效益目标:
式中,E1为水库下游年平均生态保护程度,Rt为t时段水库下游的生态保护程度,T为总时段数。
(32)建立发电效益目标:
式中,E2为水库水电站的年发电量,△t为计算时段,k为水电站出力系数,qt为t时段水轮机过机流量,Ht为t时段水电站净水头;
(33)多目标处理方法
对水库调度图进行优化,需同时考虑生态效益目标和经济效益目标,将多目标问题转化成单目标处理,总目标计算公式为:
式中,E为总目标,Ei'为权重系数子目标的均值归一化后的子目标,Ei为权重系数子目标,为权重系数子目标的均值,αi为权重系数,
(4)约束条件集
水库调度图优化模型中优化变量是水库调度图调度线的水位高度,在优化过程中需要满足以下约束条件集:
①调度线形状约束:
式中,为第j条调度线在t时段对应的水库水位,和分别为紧邻第j条调度线之下的第j-1条调度线和紧邻第j条调度线之上的第j+1条调度线在t时段对应的水库水位;和分别为t时段对应的水库允许的最低蓄水位和最高蓄水位;
②水量平衡约束:
Vt+1=Vt+(It-Qt)×△t,t=1,2,…T; (9)
式中,It和Qt分别为t时段水库的入库流量和下泄流量,Vt和Vt+1分别为t时段和t+1时段的水库蓄水量,△t为时段间隔;
③水库库容约束:
Vt min≤Vt≤Vt max,t=1,2,…T; (10)
式中,Vt min和Vt max分别为t时段水库允许的最小蓄水量和最大蓄水量;
④电站出力约束
式中,Nt为t时段的水库出力,和分别为t时段水库允许的最小出力和最大出力;
⑤水库下泄流量约束:
式中,和分别为t时段水库允许的最小下泄流量和最大下泄流量。
进一步的,所述步骤(4)包括:
(41)将常规调度图作为初始调度图C0,调度线条数为S,每条调度线有T个点;
(42)固定除当前调度线L当前阶段的标点以外的其他标点,对当前标点在其可行域范围内进行寻优,其中可行域内任意一个点与其他所有固定标点构成新的调度图,根据长系列历史径流资料,按照出力要求模拟调度过程,计算效益值,取其中最大效益的点作为当前调度线当前时段的标点值;
(43)对当前调度线L下一个阶段的标点进行优化,直到对调度线L所有标点完成优化;
(44)计算本次优化得到的调度线L(m)与上一次优化得到的调度线L(m-1)之间的距离,判断是否满足精度要求,若满足进行步骤(45),否则,重复步骤(42)~步骤(43)直到相邻两轮优化的调度线之间的距离满足精度要求;
(45)从下到上,按步骤(42)~步骤(44)对下一调度线优化,直到调度图中所有调度线完成寻优,并记录优化后的调度图C1及其效益值E1,完成第一轮寻优;
(46)计算本次优化得到的调度图C(n)与上一次优化得到的调度图C(n-1)之间的距离,判断是否满足精度要求,若满足结束迭代,否则,重复步骤(42)~步骤(45)进行迭代寻优,直到前后两轮调度图的距离满足一定精度,整个迭代程序结束;
其中,步骤(44)和步骤(46)中调度线之间的距离和调度图之间的距离定义如下:定义调度线f与调度线g之间的距离为其中fk和gk分别表示在k时段f调度线和g调度线上的水库水位标点值;调度图之间的距离为其中,C1s,k和C2s,k分别对应于调度图C1与C2第s条调度线k时段的标点值。
进一步的,所述步骤(5)包括:
(51)下生态流量限制线绘制方法:选择多个枯水年,从供水期末死水位开始,根据本时段的来水和生态流量区间下限,运用水量平衡公式逆时序推求至蓄水期初,得到水位变化过程,取各枯水年水库水位上包线即为下生态流量限制线,具体推求过程如下式:
Vt=Vt+1+(q1t-It)×△t,t=T+1,T,T-1…2,1 (13)
式中,It为t时段水库的入库流量,Vt和Vt+1分别为t时段和t+1时段的水库蓄水量,q1t为t时段生态流量区间的下限流量,△t为时段间隔;
(52)上生态流量限制线绘制方法:选择多个丰水年,从供水期末死水位开始,根据本时段的来水和生态流量区间上限,运用水量平衡公式逆时序推求至蓄水期初,得到水位变化过程,取各丰水年水库水位下包线即为上生态流量限制线,具体推求过程如下式:
Vt=Vt+1+(q2t-It)×△t,t=T+1,T,T-1…2,1 (14)
式中,It为t时段水库的入库流量,Vt和Vt+1分别为t时段和t+1时段的水库蓄水量,q2t为t时段生态流量区间的上限流量,△t为时段间隔;
(53)调度规则:当水库水位在生态限制区时,电站以该区域的规则进行调度,当生态限制区域与调度图的其他区域重合时,以生态限制区域规则进行调度,其他区域的运行规则同常规调度图。
有益效果:与现有技术相比,本发明的方法具有以下优点:
(1)现有技术在水库调度图生态目标优化方面,通常将生态流量取为一个固定值,不能反映生态要求的年内变化情况,本发明根据鱼类不同的生长周期采用不同的流量区间对生态目标进行优化,更符合实际情况。
(2)在来水较枯(或较丰)的年份,水库水位会处于较低(或较高)的水平,继续按照调度图进行定出力调节,可能会造成生态破坏程度较重的情况,因此,本发明通过增加生态限制线,来减小生态流量破坏深度,进一步加强对下游生态的保护。
附图说明
图1是本发明方法的流程图;
图2是实施例中的常规调度图;
图3是DPSA算法求解模型流程图;
图4是含有生态限制线的水库调度图示意图;
图5是实施例中含有生态限制线的水库优化调度图;
图6是常规调度图模拟结果和优化调度图模拟结果示意图。
具体实施方式
下面以实施例和附图进一步说明本发明的技术方案。
本发明运用物理栖息地模型计算生态流量区间,建立以生态和发电为目标的优化调度模型,增加生态限制线,寻求水库总体生态、发电效益最大,同时单次生态破坏程度最小的多目标调度图。如图1所示,本发明的一种兼顾生态保护需求的水库兴利调度图优化方法,包括以下步骤:
步骤1:获取水电站水库的基本信息数据。
收集水电站水库的资料,某年调节水库的主要功能是发电同时兼顾下游生态保护的要求。正常蓄水位为675m,死水位为640m,可利用库容为8.2亿m3,装机容量为270MW,水轮机最大过流量为270m3/s。
步骤2:选取代表鱼类,采用物理栖息地模型(Physical Habitat SimulationModel,PHABSIM)建立水库下泄流量与生态保护程度的关系,确定生态保护目标对应的生态流量区间。其中,水库下泄流量对应的代表鱼类加权可利用面积的计算公式为:
Ci=f(vi)×f(di)×f(ci) (2)
式中,Ai为研究河段第i个网格的水域面积,Ci为第i个网格的组合适应度因子,f(vi),f(di),f(ci)分别为第i个网格的流速、水深、河床地质适宜度值。根据流量条件计算出各网格不同的流速、水深,由式(1)和式(2)进一步计算出各流量条件下的WUA值,得到流量-WUA关系曲线。
t时段水库下游的生态保护程度Rt的计算公式为:
式中,Rt为t时段水库下游的生态保护程度;WUAt为t时段水库下泄流量对应的代表鱼类加权可利用面积,WUAt,max为t时段对应最大加权可利用面积,建立流量与生态保护程度之间的关系。
实例中选取生态保护程度目标为70%,随着生态流量的增加,WUA值先增加后减少,当生态保护程度选取70%时对应两个WUA值,因此对应两个生态流量,构成70%生态保护程度对应的流量区间,当实际生态流量在流量区间中时生态保护程度满足70%,具体的70%生态保护程度对应的生态流量区间见表1。
表1 70%生态保护程度对应的生态流量区间
步骤3:基于生态效益、发电效益建立水库调度图优化模型
如图2所示,常规兴利调度图一般包含上基本调度线、下基本调度线以及加大出力线等调度线。通过建立以生态、发电效益为目标的调度模型对常规调度图原有的调度线进行优化,目标函数和约束条件的具体形式如下:
(1)生态效益目标:
式中,E1为水库下游年平均生态保护程度,Rt为t时段水库下游的生态保护程度;T为总时段数。
(2)发电效益目标:
式中,E2为水库水电站的年发电量,单位为kW·h;△t为计算时段,单位为h;k为水电站出力系数;qt为t时段水轮机过机流量,单位为m3/s;Ht为t时段水电站净水头。
(3)多目标处理方法
对水库调度图进行优化,需同时考虑生态效益目标和经济效益目标,将多目标问题转化成单目标处理,总目标计算公式为:
式中,E为总目标,Ei'为权重系数子目标的均值归一化后的子目标,Ei为权重系数子目标,为权重系数子目标的均值,αi为权重系数,本实例中α1=α2=0.5。
(4)约束条件
本发明的优化变量是水库调度图调度线的水位高度,在优化过程中需要满足以下G(x)约束条件集:
①调度线形状约束(即:调度线之间不能交叉)
式中,为第j条调度线在t时段对应的水库水位,和分别为紧邻第j条调度线之下的第j-1条调度线和紧邻第j条调度线之上的第j+1条调度线在t时段对应的水库水位;和分别为t时段对应的水库允许的最低蓄水位和最高蓄水位。
②水量平衡约束
Vt+1=Vt+(It-Qt)×△t,t=1,2,…T (8)
式中,It和Qt分别为t时段水库的入库流量和下泄流量;Vt和Vt+1分别为t时段和t+1时段的水库蓄水量;△t为时段间隔。
③水库库容约束
Vt min≤Vt≤Vt max,t=1,2,…T (9)
式中,Vt min和Vt max分别为t时段水库允许的最小蓄水量和最大蓄水量。
④电站出力约束
式中,Nt为t时段的水库出力;和分别为t时段水库允许的最小出力和最大出力。
⑤水库下泄流量约束
式中,和分别为t时段水库允许的最小下泄流量和最大下泄流量。
步骤4:以入库径流为模型的输入,采用逐次逼近动态规划法(DPSA)求解模型。具体算法流程图如图3所示,包括以下步骤:
(1)输入水电站水库基本资料;
(2)将常规调度图作为初始调度图C(0),调度线条数为S,每条调度线有T各点;
(3)固定除当前调度线L当前阶段的标点以外的其他标点,对当前标点在其可行域范围内进行寻优,其中可行域内任意一个点与其他所有固定标点构成新的调度图,根据长系列历史径流资料,按照出力要求模拟调度过程,计算效益值,取其中最大效益的点作为当前调度线当前时段的标点值;
(4)对当前调度线L下一个阶段的标点进行优化,直到对调度线L所有标点完成优化;
(5)计算本次优化得到的调度线L(m)与上一次优化得到的调度线L(m-1)之间的距离,判断是否满足精度要求,若满足进行步骤(6),否则,重复步骤(3)~步骤(4)直到相邻两轮优化的调度线之间的距离满足精度要求;
(6)从下到上,按步骤(3)~步骤(5)对下一调度线优化,直到调度图中所有调度线完成寻优,并记录优化后的调度图C(1)及其效益值E1,完成第一轮寻优;
(7)计算本次优化得到的调度图C(n)与上一次优化得到的调度图C(n-1)之间的距离,判断是否满足精度要求,若满足结束迭代,否则,重复步骤(3)~步骤(6)进行迭代寻优,直到前后两轮调度图的距离满足一定精度,整个迭代程序结束。
其中,步骤(5)和步骤(7)中调度线之间的距离和调度图之间的距离定义如下:定义调度线f与调度线g之间的距离为其中fk和gk分别表示在k时段f调度线和g调度线上的水库水位标点值;调度图之间的距离为其中,C1s,k和C2s,k分别对应于调度图C1与C2第s条调度线k时段的标点值。
实例中以1981~2010年的入库径流实测资料推导优化调度图,其中,DPSA方法参数取值如下:调度线之间距离精度为0.5,调度图之间的距离精度为2。
步骤5:绘制生态限制线,建立生态限制线和水库调度曲线相结合的水库调度规则,编制多目标优化调度图。包括以下步骤:
(a)下生态流量限制线绘制方法:选择多个枯水年,从供水期末死水位开始,根据本时段的来水和生态流量区间下限,运用水量平衡公式逆时序推求至蓄水期初,得到水位变化过程,取各枯水年水库水位上包线即为下生态流量限制线。具体推求过程如式(13):
Vt=Vt+1+(q1t-It)×△t,t=T+1,T,T-1…2,1 (13)
式中,It为t时段水库的入库流量,Vt和Vt+1分别为t时段和t+1时段的水库蓄水量,q1t为t时段生态流量区间的下限流量。
(b)上生态流量限制线绘制方法:选择多个丰水年,从供水期末死水位开始,根据本时段的来水和生态流量区间上限,运用水量平衡公式逆时序推求至蓄水期初,得到水位变化过程,取各丰水年水库水位下包线即为上生态流量限制线。具体推求过程如式(14):
Vt=Vt+1+(q2t-It)×△t,t=T+1,T,T-1…2,1 (14)
式中:It为t时段水库的入库流量,Vt和Vt+1分别为t时段和t+1时段的水库蓄水量,q2t为t时段生态流量区间的上限流量。
(c)调度规则:当水库水位在生态限制区时,电站以该区域的规则进行调度,当生态限制区域与调度图的其他区域重合时,以生态限制区域规则进行调度,其他区域的运行规则同常规调度图。
实例中在步骤4调度图优化的基础上增加生态限制线,以70%生态保护程度的流量区间上界绘制上生态限制线,用70%生态保护程度的流量区间下界绘制下生态限制线,得到图4所示的生态限制线的水库调度图示意图。具体的,上生态限制线其以上的区域为上生态限制区,下生态限制线以下的区域为下生态限制区。当水库水位在生态限制区时,电站以该时段对应的生态流量区间的上界或下界进行下泄,当生态限制区域与调度图的其他区域重合时,以该区域的规则进行调度,其他区域的运行规则同常规调度图。最终编制的含有生态限制线的多目标优化调度图见图5,由下到上依次为下生态限制区,降低出力区、保证出力区、1.5倍加大出力区、2倍加大出力区、满发出力区和上生态限制区。
根据绘制的多目标优化调度图和常规调度图进行水库的模拟,常规调度下,多年平均年发电量为10.99亿kW·h,平均生态保护程度为92.34%,生态破坏频率为6.11%,最低生态保护程度为53.07%。在优化调度图调度下,多年平均年发电量为11.18亿kW·h,比常规增加1.72%,平均生态保护程度为92.61%,比常规增加0.27%;生态破坏频率为1.11%,比常规降低5%;最低生态保护程度为60.84%。结果表明优化调度图可以有效提高水库的发电量及下游整体的生态保护程度,同时,当发生生态破坏情况时有效减小其破坏深度。
水库分别依据常规调度图和优化调度图运行的下泄流量过程如图6所示,位于流量区间之外的流量不能达到70%生态保护程度,用‘o’表示。定义生态保护程度低70%的月份为生态破坏月份,经优化调度生态破坏发生的频率明显降低。另外,在优化调度的下泄流量的年内差异较常规调度减小。
Claims (6)
1.一种兼顾生态保护需求的水库兴利调度图优化方法,其特征在于,包括以下步骤:
(1)获取水电站水库的基本信息数据;
(2)选取代表鱼类,采用物理栖息地模型PHABSIM建立下泄流量与生态保护程度之间的关系,确定生态保护目标对应的生态流量区间;
(3)基于发电效益和生态效益建立水库调度图优化模型;
(4)以入库径流为水库调度图优化模型的输入,采用DPSA方法求解模型;
(5)绘制生态限制线,建立生态限制线和水库调度曲线相结合的水库调度规则,编制多目标优化调度图。
2.根据权利要求1所述的一种兼顾生态保护需求的水库兴利调度图优化方法,其特征在于,所述步骤(2)中采用物理栖息地模型PHABSIM建立下泄流量与生态保护程度之间的关系,其中,水库下泄流量对应的代表鱼类加权可利用面积的计算公式为:
<mrow>
<mi>W</mi>
<mi>U</mi>
<mi>A</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
Ci=f(vi)×f(di)×f(ci) (2)
式中,Ai为研究河段第i个网格的水域面积;Ci为第i个网格的组合适应度因子;f(vi),f(di),f(ci)分别为第i个网格的流速、水深、河床地质适宜度值;根据流量条件计算出各网格不同的流速、水深,由式(1)和式(2)进一步计算出各流量条件下的WUA值,得到流量-WUA关系曲线;
t时段水库下游的生态保护程度Rt的计算公式为:
<mrow>
<msub>
<mi>R</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>WUA</mi>
<mi>t</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>WU</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,Rt为t时段水库下游的生态保护程度;WUAt为t时段水库下泄流量对应的代表鱼类加权可利用面积,WUAt,max为t时段对应最大加权可利用面积,建立流量与生态保护程度之间的关系。
3.根据权利要求1所述的一种兼顾生态保护需求的水库兴利调度图优化方法,其特征在于,所述步骤(3)中建立水库调度图优化模型为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>max</mi>
<mi> </mi>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>{</mo>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mo>&Element;</mo>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,E(x)表示目标函数集,E1(x)表示生态效益目标函数,E2(x)表示发电效益目标函数,G(x)表示约束条件集。
4.根据权利要求3所述的一种兼顾生态保护需求的水库兴利调度图优化方法,其特征在于,所述建立水库调度图优化模型的方法包括:
(31)建立生态效益目标:
<mrow>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<msub>
<mi>R</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,E1为水库下游年平均生态保护程度,Rt为t时段水库下游的生态保护程度,T为总时段数。
(32)建立发电效益目标:
<mrow>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<msub>
<mi>kq</mi>
<mi>t</mi>
</msub>
<msub>
<mi>H</mi>
<mi>t</mi>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,E2为水库水电站的年发电量,△t为计算时段,k为水电站出力系数,qt为t时段水轮机过机流量,Ht为t时段水电站净水头;
(33)多目标处理方法
对水库调度图进行优化,需同时考虑生态效益目标和经济效益目标,将多目标问题转化成单目标处理,总目标计算公式为:
<mrow>
<mi>max</mi>
<mi> </mi>
<mi>E</mi>
<mo>=</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
<mo>,</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</munderover>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<msubsup>
<mi>E</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mover>
<msub>
<mi>E</mi>
<mi>i</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>/</mo>
<msub>
<mi>E</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,E为总目标,E′i为权重系数子目标的均值归一化后的子目标,Ei为权重系数子目标,为权重系数子目标的均值,αi为权重系数,
(4)约束条件集
水库调度图优化模型中优化变量是水库调度图调度线的水位高度,在优化过程中需要满足以下约束条件集:
①调度线形状约束:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>Z</mi>
<mi>t</mi>
<mi>min</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>d</mi>
<mi>t</mi>
<mrow>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>d</mi>
<mi>t</mi>
<mi>j</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>d</mi>
<mi>t</mi>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>Z</mi>
<mi>t</mi>
<mi>max</mi>
</msubsup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mi>T</mi>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mi>k</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,为第j条调度线在t时段对应的水库水位,和分别为紧邻第j条调度线之下的第j-1条调度线和紧邻第j条调度线之上的第j+1条调度线在t时段对应的水库水位;和分别为t时段对应的水库允许的最低蓄水位和最高蓄水位;
②水量平衡约束:
Vt+1=Vt+(It-Qt)×△t,t=1,2,…T; (9)
式中,It和Qt分别为t时段水库的入库流量和下泄流量,Vt和Vt+1分别为t时段和t+1时段的水库蓄水量,△t为时段间隔;
③水库库容约束:
<mrow>
<msubsup>
<mi>V</mi>
<mi>t</mi>
<mi>min</mi>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>V</mi>
<mi>t</mi>
</msub>
<mo>&le;</mo>
<msubsup>
<mi>V</mi>
<mi>t</mi>
<mi>max</mi>
</msubsup>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mi>T</mi>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,和分别为t时段水库允许的最小蓄水量和最大蓄水量;
④电站出力约束
<mrow>
<msubsup>
<mi>N</mi>
<mi>t</mi>
<mi>min</mi>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>N</mi>
<mi>t</mi>
</msub>
<mo>&le;</mo>
<msubsup>
<mi>N</mi>
<mi>t</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msubsup>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mi>T</mi>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,Nt为t时段的水库出力,和分别为t时段水库允许的最小出力和最大出力;
⑤水库下泄流量约束:
<mrow>
<msubsup>
<mi>Q</mi>
<mi>t</mi>
<mi>min</mi>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>Q</mi>
<mi>t</mi>
</msub>
<mo>&le;</mo>
<msubsup>
<mi>Q</mi>
<mi>t</mi>
<mi>max</mi>
</msubsup>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mi>T</mi>
<mo>;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,和分别为t时段水库允许的最小下泄流量和最大下泄流量。
5.根据权利要求1所述的一种兼顾生态保护需求的水库兴利调度图优化方法,其特征在于,所述步骤(4)包括:
(41)将常规调度图作为初始调度图C0,调度线条数为S,每条调度线有T个点;
(42)固定除当前调度线L当前阶段的标点以外的其他标点,对当前标点在其可行域范围内进行寻优,其中可行域内任意一个点与其他所有固定标点构成新的调度图,根据长系列历史径流资料,按照出力要求模拟调度过程,计算效益值,取其中最大效益的点作为当前调度线当前时段的标点值;
(43)对当前调度线L下一个阶段的标点进行优化,直到对调度线L所有标点完成优化;
(44)计算本次优化得到的调度线L(m)与上一次优化得到的调度线L(m-1)之间的距离,判断是否满足精度要求,若满足进行步骤(45),否则,重复步骤(42)~步骤(43)直到相邻两轮优化的调度线之间的距离满足精度要求;
(45)从下到上,按步骤(42)~步骤(44)对下一调度线优化,直到调度图中所有调度线完成寻优,并记录优化后的调度图C1及其效益值E1,完成第一轮寻优;
(46)计算本次优化得到的调度图C(n)与上一次优化得到的调度图C(n-1)之间的距离,判断是否满足精度要求,若满足结束迭代,否则,重复步骤(42)~步骤(45)进行迭代寻优,直到前后两轮调度图的距离满足一定精度,整个迭代程序结束;
其中,步骤(44)和步骤(46)中调度线之间的距离和调度图之间的距离定义如下:定义调度线f与调度线g之间的距离为其中fk和gk分别表示在k时段f调度线和g调度线上的水库水位标点值;调度图之间的距离为其中,C1s,k和C2s,k分别对应于调度图C1与C2第s条调度线k时段的标点值。
6.根据权利要求1所述的一种兼顾生态保护需求的水库兴利调度图优化方法,其特征在于,所述步骤(5)包括:
(51)下生态流量限制线绘制方法:选择多个枯水年,从供水期末死水位开始,根据本时段的来水和生态流量区间下限,运用水量平衡公式逆时序推求至蓄水期初,得到水位变化过程,取各枯水年水库水位上包线即为下生态流量限制线,具体推求过程如下式:
Vt=Vt+1+(q1t-It)×△t,t=T+1,T,T-1…2,1 (13)
式中,It为t时段水库的入库流量,Vt和Vt+1分别为t时段和t+1时段的水库蓄水量,q1t为t时段生态流量区间的下限流量,△t为时段间隔;
(52)上生态流量限制线绘制方法:选择多个丰水年,从供水期末死水位开始,根据本时段的来水和生态流量区间上限,运用水量平衡公式逆时序推求至蓄水期初,得到水位变化过程,取各丰水年水库水位下包线即为上生态流量限制线,具体推求过程如下式:
Vt=Vt+1+(q2t-It)×△t,t=T+1,T,T-1…2,1 (14)
式中,It为t时段水库的入库流量,Vt和Vt+1分别为t时段和t+1时段的水库蓄水量,q2t为t时段生态流量区间的上限流量,△t为时段间隔;
(53)调度规则:当水库水位在生态限制区时,电站以该区域的规则进行调度,当生态限制区域与调度图的其他区域重合时,以生态限制区域规则进行调度,其他区域的运行规则同常规调度图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710963168.6A CN107704965B (zh) | 2017-10-17 | 2017-10-17 | 一种兼顾生态保护需求的水库兴利调度图优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710963168.6A CN107704965B (zh) | 2017-10-17 | 2017-10-17 | 一种兼顾生态保护需求的水库兴利调度图优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107704965A true CN107704965A (zh) | 2018-02-16 |
CN107704965B CN107704965B (zh) | 2019-03-26 |
Family
ID=61184411
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710963168.6A Expired - Fee Related CN107704965B (zh) | 2017-10-17 | 2017-10-17 | 一种兼顾生态保护需求的水库兴利调度图优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107704965B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109118120A (zh) * | 2018-09-10 | 2019-01-01 | 武汉大学 | 考虑水库调度方案可持续性利用的多目标决策方法 |
CN110851977A (zh) * | 2019-11-06 | 2020-02-28 | 武汉大学 | 基于生态流量的供水-发电-生态多目标调度图优化方法 |
CN112132488A (zh) * | 2020-10-10 | 2020-12-25 | 陕西省引汉济渭工程建设有限公司 | 一种耦合建模-优化-优选的水库生态调度方法 |
CN112766593A (zh) * | 2021-01-28 | 2021-05-07 | 中国电建集团北京勘测设计研究院有限公司 | 一种水利水电工程生态调度方案的优化方法 |
CN113052378A (zh) * | 2021-03-24 | 2021-06-29 | 水利部海河水利委员会水资源保护科学研究所 | 一种基于输水效率与生态效益的流域多水库联合调度方法 |
CN117973706A (zh) * | 2024-04-02 | 2024-05-03 | 水利部交通运输部国家能源局南京水利科学研究院 | 融合鱼类生态流量过程水温过程需求的水库生态调度方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106485346A (zh) * | 2016-09-18 | 2017-03-08 | 武汉大学 | 一种混联水库提前蓄水多目标优化调度方法 |
CN106873372A (zh) * | 2017-03-22 | 2017-06-20 | 中国水利水电科学研究院 | 基于防洪调度数据自适应控制的水库防洪调度优化方法 |
-
2017
- 2017-10-17 CN CN201710963168.6A patent/CN107704965B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106485346A (zh) * | 2016-09-18 | 2017-03-08 | 武汉大学 | 一种混联水库提前蓄水多目标优化调度方法 |
CN106873372A (zh) * | 2017-03-22 | 2017-06-20 | 中国水利水电科学研究院 | 基于防洪调度数据自适应控制的水库防洪调度优化方法 |
Non-Patent Citations (4)
Title |
---|
JIANZHONGZHOU等: "Integrated optimization of hydroelectric energy in the upper and middle Yangtze River", 《RENEWABLE AND SUSTAINABLE ENERGY REVIEWS》 * |
王海霞: "考虑生态目标的水库引水与供水联合调度研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
马立亚等: "基于DPSA的梯级水库群优化调度", 《中国水利水电科学研究院学报》 * |
黄显峰等: "基于WUA的戛洒江一级水库生态调度研究", 《中国农村水利水电》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109118120A (zh) * | 2018-09-10 | 2019-01-01 | 武汉大学 | 考虑水库调度方案可持续性利用的多目标决策方法 |
CN109118120B (zh) * | 2018-09-10 | 2020-06-23 | 武汉大学 | 考虑水库调度方案可持续性利用的多目标决策方法 |
CN110851977A (zh) * | 2019-11-06 | 2020-02-28 | 武汉大学 | 基于生态流量的供水-发电-生态多目标调度图优化方法 |
CN110851977B (zh) * | 2019-11-06 | 2023-01-31 | 武汉大学 | 基于生态流量的供水-发电-生态多目标调度图优化方法 |
CN112132488A (zh) * | 2020-10-10 | 2020-12-25 | 陕西省引汉济渭工程建设有限公司 | 一种耦合建模-优化-优选的水库生态调度方法 |
CN112132488B (zh) * | 2020-10-10 | 2024-02-02 | 陕西省引汉济渭工程建设有限公司 | 一种耦合建模-优化-优选的水库生态调度方法 |
CN112766593A (zh) * | 2021-01-28 | 2021-05-07 | 中国电建集团北京勘测设计研究院有限公司 | 一种水利水电工程生态调度方案的优化方法 |
CN112766593B (zh) * | 2021-01-28 | 2023-06-06 | 中国电建集团北京勘测设计研究院有限公司 | 一种水利水电工程生态调度方案的优化方法 |
CN113052378A (zh) * | 2021-03-24 | 2021-06-29 | 水利部海河水利委员会水资源保护科学研究所 | 一种基于输水效率与生态效益的流域多水库联合调度方法 |
CN117973706A (zh) * | 2024-04-02 | 2024-05-03 | 水利部交通运输部国家能源局南京水利科学研究院 | 融合鱼类生态流量过程水温过程需求的水库生态调度方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107704965B (zh) | 2019-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107704965A (zh) | 一种兼顾生态保护需求的水库兴利调度图优化方法 | |
Jiang et al. | Multi-stage progressive optimality algorithm and its application in energy storage operation chart optimization of cascade reservoirs | |
CN110851977B (zh) | 基于生态流量的供水-发电-生态多目标调度图优化方法 | |
CN105243438A (zh) | 一种考虑径流不确定性的多年调节水库优化调度方法 | |
Feng et al. | Optimizing electrical power production of hydropower system by uniform progressive optimality algorithm based on two-stage search mechanism and uniform design | |
CN112184070B (zh) | 一种协同生态流量需求的梯级水电站多目标优化调度方法及系统 | |
CN106502220B (zh) | 梯级泵站输水系统优化运行‑控制耦合协调方法及系统 | |
CN107563538B (zh) | 大电网平台下的面向关键水位控制的多目标水库群调度优化方法 | |
CN106529166A (zh) | 一种基于maepso算法的区域水资源优化配置方法 | |
Zhou et al. | Boosting hydropower output of mega cascade reservoirs using an evolutionary algorithm with successive approximation | |
CN102855591A (zh) | 梯级水库群短期联合发电优化调度方法及系统 | |
CN107609679A (zh) | 一种年调节水库发电调度图的多参数优选绘制方法及系统 | |
CN104965954B (zh) | 一种基于全特性空间曲面建模的梯级负荷经济分配方法 | |
Khan et al. | Optimization and simulation of reservoir operation with sediment evacuation: a case study of the Tarbela Dam, Pakistan | |
CN110400232A (zh) | 考虑电网断面约束的水电站群月度交易计划电量分解方法 | |
CN101714186A (zh) | 兼顾人类和生态需求的供水型水库调度图优化确定方法 | |
CN106873372A (zh) | 基于防洪调度数据自适应控制的水库防洪调度优化方法 | |
CN104504455B (zh) | 一种梯级蓄能控制下水电站群长期优化调度方法 | |
CN106056267A (zh) | 一种并联水库群优化调度方法 | |
Jiang et al. | Research and application of key technologies in drawing energy storage operation chart by discriminant coefficient method | |
CN101714193B (zh) | 面向河流生态系统保护的发电型水库调度函数优化方法 | |
CN109300058A (zh) | 一种特大流域水电站群优化调度两阶段直接搜索降维方法 | |
CN114331033B (zh) | 梯级水库消落期运行水位协同调度方法及装置 | |
Wu et al. | Multi-objective ecological operation model of cascade hydropower reservoirs | |
CN106300441A (zh) | 一种面向电网多级输电断面的水电站机组负荷分配方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190326 |
|
CF01 | Termination of patent right due to non-payment of annual fee |