CN107862457B - 一种水库分期调度规则提取方法 - Google Patents

一种水库分期调度规则提取方法 Download PDF

Info

Publication number
CN107862457B
CN107862457B CN201711068078.7A CN201711068078A CN107862457B CN 107862457 B CN107862457 B CN 107862457B CN 201711068078 A CN201711068078 A CN 201711068078A CN 107862457 B CN107862457 B CN 107862457B
Authority
CN
China
Prior art keywords
reservoir
water level
svm
coefficient
month
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.)
Active
Application number
CN201711068078.7A
Other languages
English (en)
Other versions
CN107862457A (zh
Inventor
方国华
郭玉雪
闻昕
郭枫
丁紫玉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN201711068078.7A priority Critical patent/CN107862457B/zh
Publication of CN107862457A publication Critical patent/CN107862457A/zh
Application granted granted Critical
Publication of CN107862457B publication Critical patent/CN107862457B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Development Economics (AREA)
  • Software Systems (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Operations Research (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Artificial Intelligence (AREA)
  • Quality & Reliability (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Game Theory and Decision Science (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于GRA和PSO‑SVM的水库分期调度规则提取方法,步骤如下:在历史实际调度结果基础上,首先确定决策变量和影响因子属性集,基于GRA筛选分期影响因子;确定训练样本和测试样本,构建SVM模型;基于样本数据,建立PSO‑SVM模拟模型,运用PSO对SVM模型进行参数率定并验证;确定水库分期调度规则。本发明能够较好地继承实际调度过程中的综合效益,实现在来水不确定的前提下对水库实际运行的指导的目的,为水库实际运行提供理论依据和科学支撑。

Description

一种水库分期调度规则提取方法
技术领域
本发明属于水利水电领域的水库调度技术领域,特别是一种基于GRA和PSO-SVM组合的水库分期调度规则提取方法。
背景技术
目前,水库实际调度仍然根据过往经验进行操作,如何从实际调度经验中提取水库调度规则是水库运行调度研究的一个重要方向。调度函数是当下水库调度规则的表现形式之一。调度函数提取方法主要包含统计回归法及智能算法两大类。然而目前大部分调度函数研究并没有进行影响因子筛选操作,迄今为止影响因子的选取方法研究尚在起步。初步选定的影响因子较多,若全部考虑必会增加模型的难度,也会降低模型模拟的准确性,如何从这些因子中找到关键因子,对应用相关性较高的影响因子进行有效地调度规则模拟,具有重要的实际意义。
发明内容
发明目的:针对现有技术中存在的问题,本发明提供一种基于GRA和PSO-SVM组合的水库分期调度规则提取方法,解决目前大部分调度规则研究并没有进行影响因子筛选操作导致模拟效率低问题,同时基于PSO-SVM提供精度较高的水库调度函数,为水电站实际运行提供理论依据和科学支撑。
技术方案:一种基于GRA和PSO-SVM组合的水库分期调度规则提取方法,该方法包括:
步骤1:基于水库历史实际调度结果,确定决策变量和影响因子属性集;
步骤2:基于灰色关联度分析(GRA)筛选分期影响因子;
步骤3:将各个月份下确定的影响因子作为输入向量,决策变量水库逐月末水位Z作为输出向量;确定训练样本和测试样本,分别确定各个月份下的训练样本和测试样本,其中,样本数为M,训练样本数为N,测试样本数为M-N;
步骤4:基于支持向量机建立水库逐月末水位Z和影响因子的逐月SVM模拟模型
Figure GDA0003151511040000011
其中,K(xi,yj)为核函数,xi为由影响因子构成的第i个输入变量,yj为由水库逐月末水位构成的第j个输出向量。本发明选用RBF核参数,b为常数,αi为拉格朗日乘子,i=1,2,…n;
步骤5:基于样本数据,构建PSO-SVM模拟模型,运用PSO对SVM模型进行参数率定并验证,确定水库分期调度规则。
作为本发明进一步的改进,步骤2具体包括:
步骤21,在各个月份下,分别构建由决策变量(逐月末水位Z)构成的参考数列和由影响因子(水库时段初水位Z0、天然来水Q、叠加水位Zα、入库水能Ef、水库蓄能Es、水库入能与蓄能交互项Efs)构成的比较数列,并进行无量纲化处理:
参考数列:X(0)={Z1,Z2,…,Zt…,ZT} (1)
比较数列:
Figure GDA0003151511040000021
式中,t为参考数列个数,t=1,2,...,T;
步骤22,求参考数列与比较数列的灰色关联系数:
Figure GDA0003151511040000022
式中,X(i)为比较数列X的第i行;X(0)为参考数列;X(i)(t),X(0)(t)分别为比较数列X的第i行、参考数列的第t个变量;ρ为分辨系数,一般在0~1中间;
步骤23,计算影响因子相对于决策变量的关联度,其中关联度:
Figure GDA0003151511040000023
为参考数列X(0)(t)和比较数列X(i)(t)在第t点的关联系数。
步骤24,进行关联排序,根据各个月份的关联度排序结果,分别确定逐月份与水库时段出力Nt关联性较强(关联度值/最大关联度值>0.50)的因素作为最终调度规则模拟的影响因子。
作为本发明进一步的改进,步骤5具体包括:
步骤51,设定PSO算法的参数,确定初始种群规模G,搜索空间为D维,全局迭代次数MAXGEN,惯性权重系数w,决策变量离子速度[Vmin,Vmax]、核参数[gmin,gmax]和惩罚参数[Cmin,Cmax],其中,Vmin、Vmax分别为离子速度最小/大值,gmin、gmax分别为松弛变量最小/大值,Cmin、Cmax分别为惩罚参数最小/大值;定义第j个粒子位置表示为向量Xj=(xj1,xj2,…,xjD),速度表示为向量Vj=(vj1,vj2,…,vjD),历史最优位置表示为Qj=(qj1,qj2,…,qjD);
步骤52,令GEN=1,初始化所述离子速度V、惩罚参数C和核参数g,产生G组初始种群,包括初始位置向量与初始速度向量;
步骤53,将初始种群送入所述PSO-SVM模型进行训练,所述SVM模型在训练样本集上的采用纳什系数NSE以及确定性系数R2组合作为适应度函数值f=αNSE+(1-α)R2,NSE、R2分别按式(4)、(5)计算:
Figure GDA0003151511040000031
Figure GDA0003151511040000032
式中,α权重系数,α∈[0,1];
Figure GDA0003151511040000033
为t时刻实际末水位值,m;
Figure GDA0003151511040000034
为t时刻模拟水位值,m;
Figure GDA0003151511040000035
分别为实际末水位平均值和模拟水位平均值,m;T为计算时段。
步骤54,设置初始种群的个体极值
Figure GDA0003151511040000036
及全局极值
Figure GDA0003151511040000037
其中,
Figure GDA0003151511040000038
表示第1次迭代第j个粒子极值;
步骤55,分别利用
Figure GDA0003151511040000039
Figure GDA00031515110400000310
更新粒子当前的速度和位置,并且保证
Figure GDA00031515110400000311
其中,B1和B2为加速因子,B1是用来调节粒子飞向其最优位置的步长,B2是用来调节粒子飞向整个群体最优位置的步长;r1和r2为[0,1]范围之间的随机数;w为权重;
Figure GDA00031515110400000312
表示第t次迭代第j个粒子第n决策变量的飞行速度;
Figure GDA0003151511040000041
表示第t次迭代第j个粒子第n决策变量位置;
Figure GDA0003151511040000042
表示第j个粒子第n决策变量的历史最优位置;qgn表示第n决策变量的历史全局最优位置;
步骤56,根据
Figure GDA0003151511040000043
更新当前个体至最优pbest,同时根据
Figure GDA0003151511040000044
更新群体至最优gbest,其中t表示第t代循环,
Figure GDA0003151511040000045
分别表示第t、t+1次迭代第j个粒子,
Figure GDA0003151511040000046
表示第t、t+1次迭代第j个粒子适应性函数值;
步骤57,令GEN=GEN+1,判断当前迭代次数GEN是否达到全局迭代次数MAXGEN,若GEN=MAXGEN,则输出最优个体,进行步骤58,否则返回步骤53;
步骤58,将最优参数输入所述SVM模型,利用测试样本集对数据进行参数验证,并计算测试样本集的纳什系数NSE以及确定性系数R2,若满足要求参数要求,则输出最优参数,确定最优调度规则,否则转入进入52。
作为本发明进一步的改进,步骤1中的决策变量为水库逐月末水位Z;
作为本发明进一步的改进,步骤1中的影响因子集为:水库时段初水位Z0、天然来水Q、叠加水位Zα、入库水能Ef、水库蓄能Es、水库入能与蓄能交互项Efs
作为本发明进一步的改进,步骤5中确定的水库调度规则需要同时满足训练样本和测试样本的纳什系数NSE以及确定性系数R2:R2>0.5,NSE>0.5。
作为本发明进一步的改进,步骤5中适应度函数值f=αNSE+(1-α)R2中的α依据水库调度决策者对纳什系数和确定性系数的重要程度确定。
本发明采用上述技术方案,具有以下有益效果:
(1)能够较好地继承实际调度过程中的综合效益,实现在来水不确定的前提下对水库实际运行的指导的目的;
(2)将GRA与SVM进行结合实现对水库调度规则的分期提取,能够均衡决策变量与影响因子对于水库调度中的不同时段二者之间的相关性并不一致的问题,并平衡单一模型由模型结构带来的不确定性,提高模型模拟精度。
(3)基于SVM模拟进行水库调度模拟能够提取具体调度函数,实现水库调度规则可操作性;
(4)基于PSO对SVM模型进行参数率定与验证,提高模型模拟准确性。
附图说明
图1是本发明的方法流程图;
图2是步骤3的具体流程示意图;
图3是PSO-SVM模拟模型的长系列实际水位值与实际水位值对比图。
具体实施方式
下面结合具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
附图1为本发明的一种基于GRA和PSO-SVM组合的水库分期调度规则提取方法的流程图,具体步骤为:
步骤1:基于水库历史实际调度结果,确定决策变量和影响因子属性集;
步骤2:基于灰色关联度分析(GRA)筛选分期影响因子,具体包括:
步骤21,在各个月份下,分别构建由决策变量(逐月末水位Z)构成的参考数列和由影响因子(水库时段初水位Z0、天然来水Q、叠加水位Zα、入库水能Ef、水库蓄能Es、水库入能与蓄能交互项Efs)构成的比较数列,并进行无量纲化处理:
参考数列:X(0)={Z1,Z2,…,Zt…,ZT} (1)
比较数列:
Figure GDA0003151511040000051
式中,t为参考数列个数,t=1,2,…,T;
步骤22,求参考数列与比较数列的灰色关联系数:
Figure GDA0003151511040000061
式中,X(i)为比较数列X的第i行;X(0)为参考数列;X(i)(t),X(0)(t)分别为比较数列X的第i行、参考数列的第t个变量;ρ为分辨系数,一般在0~1中间。
步骤23,计算影响因子相对于决策变量的关联度,其中关联度:
Figure GDA0003151511040000062
为参考数列X(0)(t)和比较数列X(i)(t)在第t点的关联系数。
步骤24,进行关联排序,根据各个月份的关联度排序结果,分别确定逐月份与水库时段出力Nt关联性较强(关联度值/最大关联度值>0.50)的因素作为最终调度规则模拟的影响因子。
步骤3:将各个月份下确定的影响因子作为输入向量,决策变量水库逐月末水位Z作为输出向量;确定训练样本和测试样本,分别确定各个月份下的训练样本和测试样本,其中,样本数为M,训练样本数为N,测试样本数为M-N;
步骤4:基于支持向量机建立水库逐月末水位Z和影响因子的逐月SVM模拟模型
Figure GDA0003151511040000063
其中,K(xi,yj)为核函数,xi为由影响因子构成的第i个输入变量,yj为由水库逐月末水位构成的第j个输出向量本发明选用RBF核参数,b为常数,αi为拉格朗日乘子,i=1,2,…n;
附图2为步骤5:基于样本数据,构建PSO-SVM模拟模型,运用PSO对SVM模型进行参数率定并验证,确定水库分期调度规则,具体包括:
步骤51,设定PSO算法的参数,确定初始种群规模G,搜索空间为D维,全局迭代次数MAXGEN,惯性权重系数w,决策变量离子速度[Vmin,Vmax]、核参数[gmin,gmax]和惩罚参数[Cmin,Cmax],其中,Vmin、Vmax分别为离子速度最小/大值,gmin、gmax分别为松弛变量最小/大值,Cmin、Cmax分别为惩罚参数最小/大值;定义第j个粒子位置表示为向量Xj=(xj1,xj2,…,xjD),速度表示为向量Vj=(vj1,vj2,…,vjD),历史最优位置表示为Qj=(qj1,qj2,…,qjD);
步骤52,令GEN=1,初始化所述离子速度V、惩罚参数C和核参数g,产生G组初始种群,包括初始位置向量与初始速度向量;
步骤53,将初始种群送入所述PSO-SVM模型进行训练,所述SVM模型在训练样本集上的采用纳什系数NSE以及确定性系数R2组合作为适应度函数值f=αNSE+(1-α)R2,NSE、R2分别按式(4)、(5)计算:
Figure GDA0003151511040000071
Figure GDA0003151511040000072
式中,α权重系数,α∈[0,1];
Figure GDA0003151511040000073
为t时刻实际末水位值,m;
Figure GDA0003151511040000074
为t时刻模拟水位值,m;
Figure GDA0003151511040000075
分别为实际末水位值和模拟水位平均值,m;T为计算时段。
步骤54,设置初始种群的个体极值
Figure GDA0003151511040000076
及全局极值
Figure GDA0003151511040000077
其中,Pj 1表示第1次迭代第j个粒子极值;
步骤55,分别利用
Figure GDA0003151511040000078
Figure GDA0003151511040000079
更新粒子当前的速度和位置,并且保证
Figure GDA00031515110400000710
其中,B1和B2为加速因子,B1是用来调节粒子飞向其最优位置的步长,B2是用来调节粒子飞向整个群体最优位置的步长;r1和r2为[0,1]范围之间的随机数;w为权重;
Figure GDA00031515110400000711
表示第t次迭代第j个粒子第n决策变量的飞行速度;
Figure GDA00031515110400000712
表示第t次迭代第j个粒子第n决策变量位置;
Figure GDA00031515110400000713
表示第j个粒子第n决策变量的历史最优位置;qgn表示第n决策变量的历史全局最优位置;
步骤56,根据
Figure GDA00031515110400000714
更新当前个体至最优pbest,同时根据
Figure GDA00031515110400000715
更新群体至最优gbest,其中t表示第t代循环,Xj t、Xj t+1分别表示第t、t+1次迭代第j个粒子,f(Xj t)、f(Xj t+1)表示第t、t+1次迭代第j个粒子适应性函数值;
步骤57,令GEN=GEN+1,判断当前迭代次数GEN是否达到全局迭代次数MAXGEN,若GEN=MAXGEN,则输出最优个体,进行步骤58,否则返回步骤53;
步骤58,将最优参数输入所述SVM模型,利用测试样本集对数据进行参数验证,并计算测试样本集的纳什系数NSE以及确定性系数R2,若满足要求参数要求,则输出最优参数,确定最优调度规则(如图3所示),否则转入进入52。
以新安江水电站作为实例,说明发明方法的有效性与合理性。新安江水电站是我国自行设计、自制设备和自己建设的第一座大型水电站,坝址以上河长323km,流域面积10442km2。新安江水库以发电为主,兼有防洪、灌溉、渔业、航运、旅游等综合功能,具有多年调节性能,电站装机容量810.0MW。本发明采用1962年-2008年水文年径流资料,以多年平均运行水位98.4m为起调水位,以月为时段,获取新安江水电站长系列实际调度成果。
基于GRA计算各影响因子对决策变量的关联程度,确定最终影响因子如表1所示。由表1可以看出:水库全年基本上以水库水位及其构成的影响因子所占比重比较大,尤其在供水期(10~1月)优势更为明显,只有因6月份天然来水多,水库来水量及构成的影响要素较优与水库水位及其构成的影响要素。
表2为基于PSO的SVM模型参数率定结果。表3模型模拟值在逐月水位序列的精度评定结果,纳什系数NSE以及确定性系数R2在率定期和验证期均大于0.9,说明模型模拟结果良好。
表1逐月关联度计算以及最终影响因子筛选结果
Figure GDA0003151511040000081
Figure GDA0003151511040000091
表2基于PSO的SVM模型参数率定结果
Figure GDA0003151511040000092
表3模型模拟值在逐月水位序列的精度评定结果
Figure GDA0003151511040000093

Claims (4)

1.一种基于GRA和PSO-SVM组合的水库分期调度规则提取方法,其特征在于,包括以下步骤:
步骤1:基于水库历史实际调度结果,确定决策变量和影响因子属性集;
步骤2:基于灰色关联度分析(GRA)筛选分期影响因子;
步骤21,在各个月份下,分别构建由决策变量构成的参考数列和由影响因子构成的比较数列,并进行无量纲化处理:
参考数列:X(0)={Z1,Z2,…,Zt…,ZT} (1)
比较数列:
Figure FDA0003254522210000011
式中,Z为水库逐月末水位,Z0为水库时段初水位,Q为天然来水,Zα为叠加水位,Ef为入库水能,Es为水库蓄能,Efs为水库入能与蓄能交互项,t为参考数列个数,t=1,2,…,T;
步骤22,求参考数列与比较数列的灰色关联系数:
Figure FDA0003254522210000012
式中,X(i)为比较数列X的第i行;X(0)为参考数列;X(i)(t),X(0)(t)分别为比较数列X的第i行、参考数列的第t个变量;ρ为分辨系数,在0~1中间;
步骤23,计算影响因子相对于决策变量的关联度,其中关联度:
Figure FDA0003254522210000013
为参考数列X(0)(t)和比较数列X(i)(t)在第t点的关联系数;
步骤24,进行关联排序,根据各个月份的关联度排序结果,分别确定逐月份与水库时段出力Nt关联性较强的因素作为最终调度规则模拟的影响因子;
步骤3:将各个月份下确定的影响因子作为输入向量,决策变量水库逐月末水位Z作为输出向量;确定训练样本和测试样本,分别确定各个月份下的训练样本和测试样本,其中,样本数为M,训练样本数为N,测试样本数为M-N;
步骤4:基于支持向量机建立水库逐月末水位Z和影响因子的逐月SVM模拟模型
Figure FDA0003254522210000021
其中,K(xi,yj)为核函数,选用RBF核参数;b为常数,αi为拉格朗日乘子,i=1,2,…n;
步骤5:基于样本数据,构建PSO-SVM模拟模型,运用PSO对SVM模型进行参数率定并验证,确定水库分期调度规则。
2.如权利要求1所述的一种基于GRA和PSO-SVM组合的水库分期调度规则提取方法,其特征在于,步骤5包括以下步骤:
步骤51,设定PSO算法的参数,确定初始种群规模G,搜索空间为D维,全局迭代次数MAXGEN,惯性权重系数w,决策变量离子速度[Vmin,Vmax]、核参数[gmin,gmax]和惩罚参数[Cmin,Cmax],其中,Vmin、Vmax分别为离子速度最小/大值,gmin、gmax分别为核参数最小/大值,Cmin、Cmax分别为惩罚参数最小/大值;定义第j个粒子位置表示为向量Xj=(xj1,xj2,…,xjD),速度表示为向量Vj=(vj1,vj2,…,vjD),历史最优位置表示为Qj=(qj1,qj2,…,qjD);
步骤52,令GEN=1,初始化所述离子速度V、核参数g和惩罚参数C,产生G组初始种群,包括初始位置向量与初始速度向量;
步骤53,将初始种群送入所述PSO-SVM模拟模型进行训练,所述SVM模型在训练样本集上的采用纳什系数NSE以及确定性系数R2组合作为适应度函数值f=αNSE+(1-α)R2,NSE、R2分别按式(4)、(5)计算:
Figure FDA0003254522210000022
Figure FDA0003254522210000023
式中,α为权重系数,α∈[0,1];
Figure FDA0003254522210000024
为t时刻实际末水位值,单位为m;
Figure FDA0003254522210000025
为t时刻模拟水位值,单位为m;
Figure FDA0003254522210000026
分别为实际末水位平均值和模拟水位平均值,单位为m;T为计算时段;
步骤54,设置初始种群的个体极值
Figure FDA0003254522210000031
及全局极值
Figure FDA0003254522210000032
其中,Pj 1表示第1次迭代第j个粒子极值;
步骤55,更新粒子当前的速度和位置,并且保证
Figure FDA0003254522210000033
步骤56,根据
Figure FDA0003254522210000034
更新当前个体至最优pbest,同时根据
Figure FDA0003254522210000035
更新群体至最优gbest,其中t表示第t代循环,Xj t、Xj t+1分别表示第t、t+1次迭代第j个粒子,f(Xj t)、f(Xj t+1)表示第t、t+1次迭代第j个粒子适应性函数值;
步骤57,令GEN=GEN+1,判断当前迭代次数GEN是否达到全局迭代次数MAXGEN,若GEN=MAXGEN,则输出最优个体,进行步骤58,否则返回步骤53;
步骤58,将最优参数输入所述SVM模型,利用测试样本集对数据进行参数验证,并计算测试样本集的纳什系数NSE以及确定性系数R2,若满足要求参数要求,则输出最优参数,确定最优调度规则,否则转入步骤52。
3.如权利要求1所述的一种基于GRA和PSO-SVM组合的水库分期调度规则提取方法,其特征在于,步骤5中确定的水库调度规则需要同时满足训练样本和测试样本的纳什系数NSE以及确定性系数R2:R2>0.5,NSE>0.5。
4.如权利要求1所述的一种基于GRA和PSO-SVM组合的水库分期调度规则提取方法,其特征在于,步骤5中适应度函数值f=αNSE+(1-α)R2中的α依据水库调度决策者对纳什系数和确定性系数的重要程度确定。
CN201711068078.7A 2017-11-03 2017-11-03 一种水库分期调度规则提取方法 Active CN107862457B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711068078.7A CN107862457B (zh) 2017-11-03 2017-11-03 一种水库分期调度规则提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711068078.7A CN107862457B (zh) 2017-11-03 2017-11-03 一种水库分期调度规则提取方法

Publications (2)

Publication Number Publication Date
CN107862457A CN107862457A (zh) 2018-03-30
CN107862457B true CN107862457B (zh) 2021-10-29

Family

ID=61700651

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711068078.7A Active CN107862457B (zh) 2017-11-03 2017-11-03 一种水库分期调度规则提取方法

Country Status (1)

Country Link
CN (1) CN107862457B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109359671B (zh) * 2018-09-18 2021-07-27 华中科技大学 一种水电站水库调度规则的分类智能提取方法
CN109543912B (zh) * 2018-11-29 2021-10-01 中国水利水电科学研究院 基于深度学习的水库最优调度决策模型生成方法
CN111066562A (zh) * 2019-12-24 2020-04-28 北京农业信息技术研究中心 一种葡萄霜霉病预测方法及系统
CN115640918B (zh) * 2022-12-26 2023-04-07 电子科技大学中山学院 电缆温度异常预测方法、装置、介质及设备

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106845712A (zh) * 2017-01-20 2017-06-13 中国水利水电科学研究院 面向大系统水库群的聚合分解调度规则的提取方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106845712A (zh) * 2017-01-20 2017-06-13 中国水利水电科学研究院 面向大系统水库群的聚合分解调度规则的提取方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"基于PSO_SVM 的大坝渗流监测时间序列非线性预报模型";姜谙男等;《水力学报》;20060331;第331-335页 *
"金昌市水库群联合供水优化调度研究";王鹏全;《中国优秀硕士学位论文全文数据库》;20170131;第C037-19页 *

Also Published As

Publication number Publication date
CN107862457A (zh) 2018-03-30

Similar Documents

Publication Publication Date Title
CN107862457B (zh) 一种水库分期调度规则提取方法
CN109214566B (zh) 基于长短期记忆网络的风电功率短期预测方法
Castelletti et al. A multiobjective reinforcement learning approach to water resources systems operation: Pareto frontier approximation in a single run
CN108428017A (zh) 基于核极限学习机分位数回归的风电功率区间预测方法
CN108805434A (zh) 一种基于改进nsga-ⅲ的梯级水电站多目标优化调度方法
CN113468803B (zh) 一种基于改进的woa-gru洪水流量预测方法及系统
CN105243438A (zh) 一种考虑径流不确定性的多年调节水库优化调度方法
CN106022473A (zh) 一种融合粒子群和遗传算法的基因调控网络构建方法
CN113837499A (zh) 一种超短期风电功率预测方法及系统
CN112036632B (zh) 一种基于梯级水库生态发电多目标中长期随机调度模型的优化调度方法
CN110059867A (zh) 一种swlstm结合gpr的风速预测方法
CN105634018B (zh) 一种基于随机响应面法和内点法的含风电随机最优潮流求解方法
CN111553394A (zh) 基于循环神经网络和关注机制的水库水位预测方法
CN114444378A (zh) 一种区域风电集群的短期功率预测方法
CN111008790A (zh) 一种水电站群发电调度规则提取方法
CN106295908A (zh) 一种svm风电功率预测方法
CN105809349A (zh) 一种考虑来水相关性梯级水电站群的调度方法
CN107657349A (zh) 一种水库分期发电调度规则提取方法
CN114498619A (zh) 一种风电功率预测方法及装置
CN116187835A (zh) 一种基于数据驱动的台区理论线损区间估算方法及系统
CN111860974A (zh) 一种基于状态空间和联合分布的干旱多级预测方法
CN117291069A (zh) 一种基于改进de与注意力机制的lstm污水水质预测方法
Wang et al. An improved weighted naive bayesian classification algorithm based on multivariable linear regression model
CN109146196B (zh) 一种居民社区时用水量预测方法
CN113762591A (zh) 一种基于gru和多核svm对抗学习的短期电量预测方法及系统

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