一种水光出力联合预测方法与流程

1.本发明属于水光联合发电技术领域，更具体地说，涉及一种水光出力联合预测方法。


背景技术：

2.化石能源枯竭和环境污染加剧了能源产业结构的转变和开发可再生能源发电的进程，截止2016年底，全球累计可再生能源(不含水电)装机达921gw，我国水资源丰富，水电装机容量大，基于水电调节能力大，启停灵活等优点，将水电和其他形式的新能源结合形成互补系统，能够减轻新能源大规模并网对电网产生的冲击，促进可再生能源的消纳。
3.光伏系统和梯级电站结合形成梯级水光联合发电系统，能够提高电力系统的稳定性和安全性，利用具有调节能力的水电站补偿光伏系统出力能够减轻光伏并网对电网产生的冲击，促进光伏的消纳，降低弃光率；由于光伏电站发电收到调度周期内天气变化的影响，其随机性给电力系统的日前发电调度和负荷分配带来了困难，同时，针对光伏电站出力的频繁变化，还需要对梯级水电站的出力进行及时的调控。
4.因此，目前亟需一种能有效预测水光联合发电系统出力的方法，为电网的日前调度计划提供可靠的参考，由此来提高系统并网的稳定性和安全性。


技术实现要素：

5.为了解决上述技术问题至少之一，根据本发明的一方面，提供了一种水光出力联合预测方法，包括如下步骤：
6.s10、采集光伏电站数据，采集目标光伏电站历史n1年中与目标预测日相同日期所处区间n2天，共n1×
n2＝n天的日光伏功率曲线数据及日辐照量数据，并获取目标光伏电站的出力上限p
pv,max
；
7.s20、光伏电站出力预测，利用步骤s10采集的数据进行目标预测日光伏电站出力预测，得到预测结果p
pv
(t)；
8.s30、采集梯级水电站数据，采集各目标梯级水电站历史年中平水年与目标预测日相同日期的区间径流量数据，作为对应目标梯级水电站目标预测日的区间径流量，并获取各目标梯级水电站的参数信息；
9.s40、梯级水电站出力优化，利用步骤s30采集的数据及步骤s20的预测结果p
pv
(t)，通过粒子群优化算法，优化各目标梯级水电站的出力，得到优化结果pi″
(t)，i
″
＝1,2,
…
,d，d为目标梯级水电站个数；
10.s50、构建水光出力联合模型，
11.p
total
(t)＝p1(t) p2(t)
…
pi″
(t)
…
pd(t) p
pv
(t)，
12.其中，p
total
(t)为t时刻各目标梯级水电站与目标光伏电站的总出力，p1(t)、p2(t)、...、pi″
(t)、...、pd(t)分别为各目标梯级水电站t时刻的出力，p
pv
(t)为目标光伏电站t时刻的出力。
13.根据本发明实施例的水光出力联合预测方法，可选地，所述步骤s30要采集的数据如下：
14.第i
″
个梯级水电站的水头hi；
15.第i
″
个梯级水电站的发电效率ηi；
16.梯级水电站群中目标梯级电站的个数d；
17.第i
″
个梯级水电站的出力上限pi″
，max
；
18.第i
″
个梯级水电站的允许发电流量上限
19.第i
″
个梯级水电站的允许发电流量下限
20.第i
″
个梯级水电站第t时刻的区间径流量
21.第i
″
个梯级水电站的初始库容vi″
，1
；
22.第i
″
个梯级水电站的最小库容vi″
，min
；
23.第i
″
个梯级水电站的最大库容vi″
，max
。
24.根据本发明实施例的水光出力联合预测方法，可选地，所述步骤s40梯级水电站出力优化的具体步骤为：
25.s41、构建目标函数min(f)，
[0026][0027]
其中，
[0028][0029][0030][0031]
式中，p
total
(t)为第t时刻光伏电站与各梯级水电站联合的总出力，p
aver
为整个调度周期内光伏电站与各梯级水电站联合的每小时平均出力，为第i
″
个梯级水电站第t时刻的发电流量，pi″
(t)为第i
″
个梯级水电站第t时刻的出力；
[0032]
s42、建立约束，
[0033]
0≤p
pv，t
≤p
pv，max
；
[0034]
0≤pi″
(t)≤pi″
，max
；
[0035][0036]vi
″
，min
≤vi″
，t
≤vi″
，max
；
[0037]
式中，vi″
，t
为第i
″
个梯级水电站第t时刻的库容；
[0038]
s43、设置粒子群参数，将设为未知变量，粒子群算法维数d＝d，随机产生粒子总数m＝100，最大迭代次数t
′
max
＝500，学习因子c1＝c2＝1.2，初始优化权重ω＝0.9；
[0039]
s44、开始优化，
[0040]
(1)在t时刻第z个粒子的位置为其中分别对
应各未知变量；
[0041]
在t时刻第z个粒子的速度为
[0042]
在t时刻的第z个粒子的最优位置为
[0043]
在t时刻种群搜索的全局最优位置为
[0044]
使用rand函数对粒子进行初始化，确定初始粒子位置和速度；
[0045]
(2)将粒子位置带入目标函数min(f)，求得的值为粒子适应度，并以此作为粒子个体最优值，将粒子个体最优值在全局范围内做比较，数值最小的粒子个体最优值作为全局最优值，并记录其对应的粒子位置；
[0046]
(3)按照如下模型对粒子位置和速度进行更新，更新过程中相关参数要符合s42的约束，模型为，
[0047][0048][0049]
其中，c1＝c2＝1.2，r1和r2是两个相互独立的随机数，取值范围(0，1)，ω＝0.9；是t时刻第个z粒子的个体最优值；是t时刻z粒子的全局最优值；为更新后的粒子速度，为更新后的粒子位置；
[0050]
(4)将更新后的粒子位置带入目标函数min(f)，重复步骤(2)～(3)，迭代t
′
max
次后，结束迭代，以此时全局最优值对应的粒子位置的值为各未知变量的最优解，由此得到后，结束迭代，以此时全局最优值对应的粒子位置的值为各未知变量的最优解，由此得到的最优值；
[0051]
s45、获取梯级水电站最优出力，
[0052]
将步骤s44获得的最优的值代入式得到的pi″
(t)值即为第i
″
个梯级水电站第t时刻的最优出力。
[0053]
根据本发明实施例的水光出力联合预测方法，可选地，所述步骤s10中，与目标预测日相同日期所处区间n2天具体为：与目标预测日相同日期及其前后各一周，共15天。
[0054]
根据本发明实施例的水光出力联合预测方法，可选地，所述步骤s20光伏电站出力预测包括如下步骤：
[0055]
s21、k-means聚类分析预测，利用步骤s10中n天的日光伏功率曲线，通过k-means聚类分析法预测目标预测日的光伏功率曲线，得到预测结果p1；
[0056]
s22、bp神经网络模型预测，构建bp神经网络模型，用步骤s10中采集的n天的日辐照量数据训练bp神经网络模型，用训练好的模型预测目标预测日的光伏功率曲线，得到预测结果p2；
[0057]
s23、预测结果整合，构建整合模型，将步骤s20和步骤s30的预测结果整合，整合模型如下，
[0058]
p
pv
＝p1ω1 p2ω2；
[0059]
其中，p
pv
为最终预测的日光伏功率曲线，ω1为预测结果p1所占权重，ω2预测结果p2所占权重；
[0060]
根据最终预测的日光伏功率曲线即可获得目标预测日光伏电站的出力p
pv
(t)。
[0061]
根据本发明实施例的水光出力联合预测方法，可选地，所述步骤s21的k-means聚
类分析预测具体步骤如下：
[0062]
s210、用步骤s10中采集的n天的日光伏功率曲线作为样本数据集{y1，y2，
…
，yn}，yi表示第i天的日光伏功率曲线，从中随机选取四个样本数据作为四个初始聚类中心，分别代表晴、多云、阴及雨4种典型天气类型下的日光伏功率曲线模式；
[0063]
s211、分别计算每个样本数据至各初始聚类中心的欧氏距离d
ic
，与对应初始聚类中心的欧氏距离最小的样本即归为该初始聚类中心所属簇，将样本数据更新成四簇；
[0064]
所述欧氏距离计算模型如下，
[0065]dic
＝|y
i-xc|，
[0066]
其中，yi为第i个样本，xc表示第c簇的聚类中心；
[0067]
s212、对产生的四簇样本重新计算对应簇的聚类中心，通过如下模型计算，
[0068][0069]
其中，nc为第c簇的样本个数；
[0070]
s213、重复步骤s22～步骤s23，直到达到最大迭代次数100时终止重复；
[0071]
s214、获取目标预测日的天气预报信息，确定天气预报中四种典型天气类型对应的时段，根据不同时段对应的天气类型，选取对应簇聚类中心的日光伏功率曲线中对应时段的曲线，拼接形成目标预测日的光伏功率曲线，即为预测结果p1。
[0072]
根据本发明实施例的水光出力联合预测方法，可选地，所述步骤s22的bp神经网络模型预测具体步骤如下：
[0073]
s220、用步骤s10采集的n天的日辐照量数据作为训练样本数据集一{x1，x2，
…
，xn}，xi表示第i天的辐照量数据，用步骤s10采集的m天的日光伏功率曲线作为样本数据集二{y1，y2，
…
，yn}，yi表示第i天的日光伏功率曲线；
[0074]
s221、采用如下计算模型对训练样本数据归一化处理，
[0075][0076][0077]
其中，x
min
为训练样本数据集一中的最小值，x
max
为训练样本数据集一中的最大值，y
min
为训练样本数据集二中的最小值，y
max
为训练样本数据集二中的最大值；
[0078]
s222、创建bp神经网络模型；
[0079]
s223、用归一化处理后的训练样本数据训练bp神经网络模型，设置最大迭代次数t＝1000次，允许误差ε＝0.00001，学习速率η＝0.01，当达到最大迭代次数或满足精度范围时，结束训练；
[0080]
s224、用训练好的bp神经网络模型预测目标预测日的光伏功率曲线。
[0081]
根据本发明实施例的水光出力联合预测方法，可选地，所述步骤s222的bp神经网络模型包括输入层、隐含层和输出层，输入层节点n
′
＝24，隐含层节点l＝1，输出层节点m＝24，隐含层阈值aj＝[a1，a2，
…
，a
l
]，输出层阈值bk＝[b1，b2，
…
，bm]，输入层与隐藏层间的连接权值为ω
ij
，隐含层与输出层之间的连接权值为ω
′
jk
，其中：
[0082]
隐含层输出hj，
[0083]
[0084]
输出层输出ok，
[0085][0086]
连接权值更新，
[0087]
ω
ij
(t 1)＝ω
ij
(t) η[(1-β)d(t) βd(t-1)]，
[0088]
ω
′
jk
(t 1)＝ω
′
jk
(t) η[(1-β)d(t) βd(t-1)]，
[0089]
阈值更新，
[0090][0091]bk
(t 1)＝bk(t) (y
k-ok)，
[0092]
式中，j为bp神经网络模型实际输出值和期望输出值的均方误差，t指当前时刻，t 1指下一时刻，t-1指前一时刻，β为动量因子，0≤β＜1，yk为节点k对应的的值。
[0093]
根据本发明实施例的水光出力联合预测方法，可选地，所述步骤s23中，ω1和ω2的取值为通过粒子群优化算法进行优化后的最优值。
[0094]
根据本发明实施例的水光出力联合预测方法，可选地，所述步骤s23中，ω1＝ω2＝0.5。
[0095]
本发明的水光出力联合预测方法，通过预测目标预测日光伏电站的出力，然后采用粒子群优化算法，相应优化各目标梯级水电站的出力，最终构建的水光出力联合模型对梯级水电站、光伏电站联合系统的发电计划有着指导价值，能有效减轻光伏并网对电网产生的冲击，促进光伏的消纳，降低弃光率。
附图说明
[0096]
为了更清楚地说明本发明实施例的技术方案，下面将对实施例的附图作简单地介绍，显而易见地，下面描述中的附图仅仅涉及本发明的一些实施例，而非对本发明的限制。
[0097]
图1示出了本发明的水光出力联合预测方法流程图；
[0098]
图2示出了本发明的bp神经网络模型结构示意图；
[0099]
图3示出了实施例2中三种预测结果与实际日光伏功率曲线比较图；
[0100]
图4示出了实施例2中水光联合系统出力图。
具体实施方式
[0101]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例的附图，对本发明实施例的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。
[0102]
目前，水光联合发电系统的调度存在诸多困难，光伏发电站的出力受到天气频繁变化、辐照强度的影响，使得光伏发电具有随机性，并网时易对电网产生冲击，因此需对梯级水电站的出力进行及时的调控调度，在对水光联合发电系统目标日制定调度计划时，如果能对水光联合发电系统的水光联合出力进行有效的预测，提前构建目标预测日的水光联合出力模型，将对标预测日电网的调度计划提供可靠的参考，由此能有效提高系统并网的稳定性和安全性，促进光伏的消纳，降低弃光率。
[0103]
实施例1
[0104]
本实施例的水光出力联合预测方法，如图1所示，包括如下步骤：
[0105]
s10、采集光伏电站数据，
[0106]
采集目标光伏电站历史n1年中与目标预测日相同日期所处区间n2天，共n1×
n2＝n天的日光伏功率曲线数据及日辐照量数据，更具体地说，是采集对应日期0:00～23:00的小时光伏功率曲线数据及小时辐照量数据，并获取目标光伏电站的出力上限p
pv，max
；
[0107]
进一步地，本步骤中与目标预测日相同日期所处区间n2天具体为：与目标预测日相同日期及其前后各一周，共15天，做此选择的目的是，这n2天与目标预测日日期相近，数据的相关性较高，以此作为训练数据或初始数据有更高的参考价值，能有效提高预测精度。
[0108]
s20、光伏电站出力预测，
[0109]
利用步骤s10采集的数据进行目标预测日光伏电站出力预测，得到预测结果p
pv
(t)，具体步骤如下，
[0110]
s21、k-means聚类分析预测，利用步骤s10中m天的日光伏功率曲线，通过k-means聚类分析法预测目标预测日的光伏功率曲线，得到预测结果p1；
[0111]
s22、bp神经网络模型预测，构建bp神经网络模型，用步骤s10中采集的n天的日辐照量数据训练bp神经网络模型，用训练好的模型预测目标预测日的光伏功率曲线，得到预测结果p2；
[0112]
s23、预测结果整合，构建整合模型，将步骤s20和步骤s30的预测结果整合，整合模型如下，
[0113]
p
pv
＝p1ω1 p2ω2；
[0114]
其中，p
pv
为最终预测的日光伏功率曲线，ω1为预测结果p1所占权重，ω2预测结果p2所占权重；
[0115]
根据最终预测的日光伏功率曲线即可获得目标预测日光伏电站的出力p
pv
(t)。
[0116]
s30、采集梯级水电站数据，
[0117]
采集各目标梯级水电站历史年中平水年与目标预测日相同日期的区间径流量数据，作为对应目标梯级水电站目标预测日的区间径流量，并获取各目标梯级水电站的参数信息；
[0118]
本步骤要采集的数据如下：
[0119]
第i
″
个梯级水电站的水头hi；
[0120]
第i
″
个梯级水电站的发电效率ηi；
[0121]
梯级水电站群中目标梯级电站的个数d；
[0122]
第i
″
个梯级水电站的出力上限pi″
，max
；
[0123]
第i
″
个梯级水电站的允许发电流量上限
[0124]
第i
″
个梯级水电站的允许发电流量下限
[0125]
第i
″
个梯级水电站第t时刻的区间径流量
[0126]
第i
″
个梯级水电站的初始库容vi″
，1
；
[0127]
第i
″
个梯级水电站的最小库容vi″
，min
；
[0128]
第i
″
个梯级水电站的最大库容vi″
，max
。
[0129]
s40、梯级水电站出力优化，
[0130]
利用步骤s30采集的数据及步骤s20的预测结果p
pv
(t)，通过粒子群优化算法，优化各目标梯级水电站的出力，得到优化结果pi″
(t)，i
″
＝1，2，...，d，d为目标梯级水电站个数；
[0131]
粒子群优化算法具体方法为，
[0132]
s41、构建目标函数min(f)，为了能够保证水光联合系统并网后不会对电网造成较大的冲击，对日前调度进行优化，因此本实施例以优化目标为水光联合系统总出力波动最小化构建了目标函数，
[0133][0134]
其中，
[0135][0136][0137][0138]
式中，p
total
(t)为第t时刻光伏电站与各梯级水电站联合的总出力，p
aver
为整个调度周期内光伏电站与各梯级水电站联合的每小时平均出力，为第i
″
个梯级水电站第t时刻的发电流量，pi″
(t)为第i
″
个梯级水电站第t时刻的出力，一整个调度周期为24小时；
[0139]
s42、建立约束，
[0140]
0≤p
pv，t
≤p
pv，max
；
[0141]
0≤pi″
(t)≤pi″
，max
；
[0142][0143]vi
″
，min
≤vi″
，t
≤vi″
，max
；
[0144]
式中，vi″
，t
为第i
″
个梯级水电站第t时刻的库容通过下式获得，
[0145][0146][0147][0148]
式中，为第i
″
个梯级水电站第t时刻的入库流量，为第i
″
个梯级水电站第t时刻的出库流量，为第i
″
个梯级水电站第t时刻的出库流量，为第i
″
个梯级水电站第t时刻的弃水量，vi″
，t 1
为第i
″
个梯级水电站第t 1时刻的库容，由于第i
″
个梯级水电站的初始库容已知，可通过关系式来获取各时刻下的库容；
[0149]
s43、设置粒子群参数，将设为未知变量，粒子群算法维数d＝d，随机产生粒子总数m＝100，最大迭代次数t
′
max
＝500，学习因子c1＝c2＝1.2，初始优化权重ω＝0.9；
[0150]
s44、开始优化，
[0151]
(1)在t时刻第z个粒子的位置为其中分别对应各未知变量；在t时刻第z个粒子的速度为在t时刻的第z个粒子的
最优位置为在t时刻种群搜索的全局最优位置为使用rand函数对粒子进行初始化，确定初始粒子位置和速度；
[0152]
(2)将粒子位置带入目标函数min(f)，求得的值为粒子适应度，并以此作为粒子个体最优值，将粒子个体最优值在全局范围内做比较，数值最小的粒子个体最优值作为全局最优值，并记录其对应的粒子位置；
[0153]
(3)按照如下模型对粒子位置和速度进行更新，更新过程中相关参数要符合s42的约束，模型为，
[0154][0155][0156]
其中，c1＝c2＝1.2，r1和r2是两个相互独立的随机数，取值范围(0，1)，ω＝0.9；是t时刻第个z粒子的个体最优值；是t时刻z粒子的全局最优值；为更新后的粒子速度，为更新后的粒子位置；
[0157]
(4)将更新后的粒子位置带入目标函数min(f)，重复步骤(2)～(3)，迭代t
′
max
次后，结束迭代，以此时全局最优值对应的粒子位置的值为各未知变量的最优解，由此得到后，结束迭代，以此时全局最优值对应的粒子位置的值为各未知变量的最优解，由此得到的最优值；
[0158]
s45、获取梯级水电站最优出力，
[0159]
将步骤s44获得的最优的值代入式得到的pi″
(t)值即为第i
″
个梯级水电站第t时刻的最优出力；
[0160]
本实施例的步骤s40，以水光联合系统总出力最小化为目标，且考虑到水光联合系统运行过程中功率平衡、发电单元出力限制、水库库容限制等约束，建立起粒子群优化模型，最终得到优化后的由此得出梯级水电站的最优出力。
[0161]
s50、构建水光出力联合模型，
[0162]
p
total
(t)＝p1(t) p2(t)
…
pi″
(t)
…
pd(t) p
pv
(t)，
[0163]
其中，p
total
(t)为t时刻各目标梯级水电站与目标光伏电站的总出力，p1(t)、p2(t)、...、pi″
(t)、...、pd(t)分别为各目标梯级水电站t时刻的出力，p
pv
(t)为目标光伏电站t时刻的出力；
[0164]
通过步骤s50构建的水光出力联合模型，可根据各目标梯级水电站t时刻的出力对目标预测日当天的各梯级水电站进行调度调控，由此能有效降低出力波动，从而减轻光伏并网对电网产生的冲击。
[0165]
进一步地，由于光伏电站在目标预测日的出力，受天气等因素影响随机性较大，而本实施例的水光出力联合模型及各梯级水电站的最优出力均需基于光伏电站在目标预测日的出力来确定，因此，有效预测光伏电站在目标预测日的出力，对本实施例构建的水光出力联合模型的可靠性起着重要作用，基于此，本实施例还对光伏电站出力预测步骤进行了改进。
[0166]
本实施例的光伏电站出力预测中，s21的k-means聚类分析预测，具体步骤如下：
[0167]
s210、用步骤s10中采集的n天的日光伏功率曲线作为样本数据集{y1，y2，
…
，yn}，yi表示第i天的日光伏功率曲线，从中随机选取四个样本数据作为四个初始聚类中心，分别
代表晴、多云、阴及雨4种典型天气类型下的日光伏功率曲线模式；
[0168]
s211、分别计算每个样本数据至各初始聚类中心的欧氏距离d
ic
，与对应初始聚类中心的欧氏距离最小的样本即归为该初始聚类中心所属簇，将样本数据更新成四簇；
[0169]
所述欧氏距离计算模型如下，
[0170]dic
＝|y
i-xc|，
[0171]
其中，yi为第i个样本，xc表示第c簇的聚类中心；
[0172]
s212、对产生的四簇样本重新计算对应簇的聚类中心，通过如下模型计算，
[0173][0174]
其中，nc为第c簇的样本个数；
[0175]
s213、重复步骤s22～步骤s23，直到达到最大迭代次数100时终止重复；
[0176]
s214、获取目标预测日的天气预报信息，确定天气预报中四种典型天气类型对应的时段，根据不同时段对应的天气类型，选取对应簇聚类中心的日光伏功率曲线中对应时段的曲线，拼接形成目标预测日的光伏功率曲线，即为预测结果p1；
[0177]
本实施例步骤s21的k-means聚类分析预测，通过k-means聚类分析的方法，将大量具有参考价值的光伏出力的样本数据分成典型的几类，然后利用天气预报，建立起聚类中心与各时段天气状况的映射关系，充分考虑到日内天气频繁变化造成的日光伏功率曲线的多样化和随机性的误差因素，通过步骤s210～s214，以k-means聚类分析 天气预报的模式极大的降低了误差因素的影响，更加准确的追踪日内气象变换频繁的光伏功率曲线的变化，提高了预测结果p1的精度。
[0178]
本实施例的光伏电站出力预测中，s22的bp神经网络模型预测，具体步骤如下：
[0179]
s220、用步骤s10采集的n天的日辐照量数据作为训练样本数据集一{x1，x2，
…
，xn}，xi表示第i天的辐照量数据，用步骤s10采集的m天的日光伏功率曲线作为样本数据集二{y1，y2，
…
，yn}，yi表示第i天的日光伏功率曲线；
[0180]
s221、采用如下计算模型对训练样本数据归一化处理，
[0181][0182][0183]
其中，x
min
为训练样本数据集一中的最小值，x
max
为训练样本数据集一中的最大值，y
min
为训练样本数据集二中的最小值，y
max
为训练样本数据集二中的最大值；此步骤的目的是由于辐照数据和光伏功率数据的数量级差别较大，为避免产生病态问题，加快训练过程，采用归一化处理将数据统一映射到[0，1]之间；
[0184]
s222、创建bp神经网络模型，如图2所示，本实施例中创建的bp神经网络模型包括输入层、隐含层和输出层，输入层节点n
′
＝24，隐含层节点l＝1，输出层节点m＝24，隐含层阈值aj＝[a1，a2，
…
，a
l
]，输出层阈值bk＝[b1，b2，
…
，bm]，输入层与隐藏层间的连接权值为ω
ij
，隐含层与输出层之间的连接权值为ω
′
jk
，
[0185]
其中：
[0186]
隐含层输出hj，
[0187]
[0188]
f为隐含层激励函数
[0189]
输出层输出ok，
[0190][0191]
连接权值更新，
[0192]
ω
ij
(t 1)＝ω
ij
(t) η[(1-β)d(t) βd(t-1)]，
[0193]
ω
′
jk
(t 1)＝ω
′
jk
(t) η[(1-β)d
′
(t) βd
′
(t-1)]，
[0194]
阈值更新，
[0195][0196]bk
(t 1)＝bk(t) (y
k-ok)，
[0197]
式中，j为bp神经网络模型实际输出值和期望输出值的均方误差，具体为n
′
为训练样本总数，n
′
＝n，为样本p在节点k的目标输出值，为样本p在节点k的实际输出值，t指当前时刻，t 1指下一时刻，t-1指前一时刻，β为动量因子，0≤β＜1，yk为节点k对应的的值；通过训练bp神经网络模型，不断的更新阈值和权值，使得模型的实际输出值和目标输出值的均方误差不断减小，最终完成训练的模型能有较高的预测精度；
[0198]
s223、用归一化处理后的训练样本数据训练bp神经网络模型，设置最大迭代次数t＝1000次，允许误差ε＝0.00001，学习速率η＝0.01，当达到最大迭代次数或满足精度范围时，结束训练；
[0199]
s224、用训练好的bp神经网络模型预测目标预测日的光伏功率曲线；
[0200]
本实施例步骤s22的bp神经网络模型预测，选取历史年度与目标预测日相近日期的历史数据做参考，能有效提高预测精度，bp神经网络模型预测在预测时，充分考虑到辐照强度对光伏出力的影响，利用辐照强度数据为bp神经网络模型预测的输入值，利用光伏出力数据为bp神经网络模型预测的目标输出值，不断训练优化bp神经网络模型预测，提高了预测结果p2的精度。
[0201]
进一步地，本实施例步骤s21和步骤s22互不干扰，可以同时进行，提高本方法的预测效率。
[0202]
本实施例的光伏电站出力预测中，在预测结果整合时，对p1和p2两种预测结果所占权重取值通过粒子群优化算法进行了优化，优化目的为了提高整合模型的预测精度，而平均绝对误差mae、均方误差mse和均方根误差rmse通常作为评判预测模型精度的三个指标，因此本实施例以这三个精度指标最小化作为目标，构建了目标函数f，
[0203]
f＝(f1 f2 f3)/3，
[0204]
其中，
[0205]
f1＝min(mae)，
[0206]
f2＝min(mse)，
[0207]
f3＝min(rmse)，
[0208][0209][0210][0211]
式中，mae为平均绝对误差精度指标，mse为均方误差精度指标，rmse为均方根误差精度指标，t
′
为整个运行周期，为第i
′
时刻的预测光伏功率，为第i
′
时刻的实际光伏功率；
[0212]
然后确定约束条件，
[0213]
ω1 ω2＝1，
[0214]
0≤ω1≤1，
[0215]
0≤ω2≤1；
[0216]
然后设置粒子群参数，将ω1和ω2设为未知变量，粒子群算法位数d＝2，随机产生粒子总数为100，最大迭代次数设为100，学习因子、初始优化权重与步骤s43相同；然后采用与步骤s44相同的方法开始优化，经过100次迭代，最终得到优化后的ω1和ω2的值均为0.5，由此整合模型为p
pv
＝0.5
·
p1 0.5
·
p2ω2。
[0217]
本实施例的光伏电站出力预测中，基于k-means聚类分析和天气预报相结合的预测方法能够考虑到日内天气频繁变化造成的日光伏功率曲线的多样化和随机性，通过结合预测日当天的实时天气预报信息提高光伏出力预测精度；通过使用历史辐照数据和历史光伏出力数据建立的bp神经网络预测模型充分考虑辐照强度对光伏出力的影响；再利用粒子群优化算法优化后的权重构建整合模型，将两种预测结果整合，得到的最终预测结果将天气变化及辐照强度对光伏出力预测结果的误差影响降至最低，使p
pv
具有很高的精度，从而能进一步确保本实施例构建的水光出力联合模型的可靠性。
[0218]
实施例2
[0219]
本实施例利用水光出力联合预测方法对某地水光联合系统出力进行预测，选取2019年5月15日为预测日，其中，光伏电站为100mw，梯级水电站有三个，装机容量分别是10mw，200mw和200mw，水库最大容量分别是1.9*104m3，12.16*104m3和2.93*104m3，水头分别是70m，90m和77m。
[0220]
首先对光伏电站的出力进行预测，得到预测结果如图3所示，从图3中预测结果与预测日实际值的对比表明，本发明的方法对光伏电站出力的预测精度很高，然后对三个梯级水电站的出力进行优化，优化后各梯级水电站出力如图4所示，根据水光出力联合模型，得出水光联合系统各时刻的总出力如图4所示；各时刻水光联合系统的出力结果如表1。
[0221]
表1.水光联合系统出力结果
[0222][0223]
结合表1和图4可以看出，水光联合系统的总出力曲线平稳，该调度方案下的系统出力波动仅有1.015，联合系统并入电网时不会产生较大的波动；在0：00-5：00由于没有辐照光伏系统出力为0，由三个梯级水电站承担系统出力；在6：00时刻光伏系统开始发电，联合系统总出力略有上升，随后梯级水电及时降低出力来容纳光伏出力，联合系统于7：00再次趋于平稳；18：00时刻光伏系统几乎不再发电，联合系统再次出现较小波动，梯级水电增大出力使得联合系统出力再次趋于平稳；在本发明水光出力联合模型的指导下，仅在光伏系统发电的两个临界点产生较小的波动，总体出力波动的标准差仅为1.015，由此表面本发明方法的有效性、可靠性。
[0224]
本发明所述实例仅仅是对本发明的优选实施方式进行描述，并非对本发明构思和范围进行限定，在不脱离本发明设计思想的前提下，本领域工程技术人员对本发明的技术方案作出的各种变形和改进，均应落入本发明的保护范围。