CN103500365A - 光伏发电功率预测方法和*** - Google Patents
光伏发电功率预测方法和*** Download PDFInfo
- Publication number
- CN103500365A CN103500365A CN201310430694.8A CN201310430694A CN103500365A CN 103500365 A CN103500365 A CN 103500365A CN 201310430694 A CN201310430694 A CN 201310430694A CN 103500365 A CN103500365 A CN 103500365A
- Authority
- CN
- China
- Prior art keywords
- intensity
- prediction
- solar radiation
- constantly
- photovoltaic generation
- 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
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Photovoltaic Devices (AREA)
Abstract
一种光伏发电功率预测方法和***,采用相关性分析方法对历史数据进行分析,确定辐射强度预测相关时刻和发电功率预测相关时刻。采用BP神经网络对太阳辐射强度预测样本和光伏发电功率预测样本进行训练,得到太阳辐射强度预测模型和光伏发电功率预测模型。利用太阳辐射强度预测模型计算预测日预测时刻的太阳辐射强度,利用光伏发电功率预测模型计算预测日预测时刻的光伏发电功率。利用灰色关联度分析方法去除历史数据中关联度较低的辐射强度预测相关时刻的太阳辐射强度,提高了太阳辐射强度预测精度。采用BP神经网络良好的非线性函数逼近能力,对太阳辐射强度预测样本和光伏发电功率预测样本进行训练建立预测模型,提高了预测模型的预测精度。
Description
技术领域
本发明涉及电力***技术领域,特别是涉及一种光伏发电功率预测方法和***。
背景技术
光伏发电是利用太阳能的一种有效方式,但光伏发电***的输出功率具有不连续性和不确定性的特点。光伏发电***并网运行以后会对电网产生周期性的冲击,光伏***输出功率的扰动将有可能影响电网的稳定。因此,对光伏电站的输出功率进行预测有助于统筹安排常规电源和光伏发电的协调配置,适时及时的调整调度计划,合理安排电网运行方式。
光伏发电功率预测一般采用人工智能方法,主要包括神经网络、模糊预测、数据挖掘、支持向量机等。无论采用哪一种方法,都存在因数据采集准确度低,而导致光伏发电功率预测精度低的缺点。
发明内容
基于此,有必要提供一种预测精度高的光伏发电功率预测方法和***。
一种光伏发电功率预测方法,包括以下步骤:根据历史数据,采用相关性分析方法获取辐射强度预测相关时刻和发电功率预测相关时刻,所述历史数据包括太阳辐射强度、温度和光伏发电功率,所述辐射强度预测相关时刻指太阳辐射强度与预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻,所述发电功率预测相关时刻指太阳辐射强度和光伏发电功率与所述预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻;提取所述历史数据中所述辐射强度预测相关时刻的太阳辐射强度,并进行灰色关联度分析,得到太阳辐射强度预测样本;采用BP神经网络对所述太阳辐射强度预测样本进行训练,得到太阳辐射强度预测模型;将预测日所述辐射强度预测相关时刻的太阳辐射强度作为所述太阳辐射强度预测模型的输入,计算预测日所述预测时刻的太阳辐射强度;提取所述历史数据中预测日之前所述发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及所述预测时刻的温度和太阳辐射强度,得到光伏发电功率预测样本;采用BP神经网络对所述光伏发电功率预测样本进行训练,得到光伏发电功率预测模型;将预测日所述发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及所述预测时刻的太阳辐射强度和温度作为所述光伏发电功率预测模型的输入,计算预测日所述预测时刻的光伏发电功率。
在其中一个实施例中,所述根据历史数据,采用相关性分析方法获取辐射强度预测相关时刻和发电功率预测相关时刻的步骤,包括以下步骤:
提取所述历史数据中所述预测日之前,所述预测时刻及预测时刻前设定时间段内的太阳辐射强度;
根据分别计算所述设定时间段内各时刻点的太阳辐射强度,与所述预测时刻的太阳辐射强度的相关系数,其中,cov(X,Y)=E[X-E(X)]E[Y-E(Y)],X为预测时刻的太阳辐射强度,Y为设定时间段内一时刻点的太阳辐射强度,ρXY表示X和Y的相关系数,cov(X,Y)表示X和Y的协方差,E(X)和E(Y)分别表示X和Y的期望,D(X)和D(Y)分别表示X和Y的方差;
提取所述设定时间段内,太阳辐射强度与所述预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻,得到所述辐射强度预测相关时刻;
提取所述历史数据中所述预测日之前,所述预测时刻的光伏发电功率及预测时刻前设定时间段内的太阳辐射强度和光伏发电功率;
根据分别计算所述设定时间段内各时刻点的太阳辐射强度和光伏发电功率,与所述预测时刻的光伏发电功率的相关系数,其中cov(X',Y')=E[X'-E(X')]E[Y'-E(Y')],X'为预测时刻的光伏发电功率,Y'为设定时间段内一时刻点的太阳辐射强度或光伏发电功率,ρX'Y'表示X'和Y'的相关系数,cov(X',Y')表示X'和Y'的协方差,E(X')和E(Y')分别表示X'和Y'的期望,D(X')和D(Y')分别表示X'和Y'的方差;
提取设定时间范围内,太阳辐射强度和光伏发电功率与所述预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻,得到所述发电功率预测相关时刻。
在其中一个实施例中,所述根据历史数据,采用相关性分析方法获取辐射强度预测相关时刻和发电功率预测相关时刻的步骤,包括以下步骤:
提取所述历史数据中所述预测日之前,所述预测时刻及预测时刻前若干个等间隔时刻点的太阳辐射强度;
根据分别计算所述若干个等间隔时刻点的太阳辐射强度,与所述预测时刻的太阳辐射强度的相关系数,其中,cov(X,Y)=E[X-E(X)]E[Y-E(Y)],X为预测时刻的太阳辐射强度,Y为若干个等间隔时刻点中一时刻点的太阳辐射强度,ρXY表示X和Y的相关系数,cov(X,Y)表示X和Y的协方差,E(X)和E(Y)分别表示X和Y的期望,D(X)和D(Y)分别表示X和Y的方差;
提取所述若干个等间隔时刻点中,太阳辐射强度与所述预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻点,得到所述辐射强度预测相关时刻;
提取所述历史数据中所述预测日之前,所述预测时刻的光伏发电功率及预测时刻前若干个等间隔时刻点的太阳辐射强度和光伏发电功率;
根据分别计算所述若干个等间隔时刻点的太阳辐射强度和光伏发电功率,与所述预测时刻的光伏发电功率的相关系数,其中cov(X',Y')=E[X'-E(X')]E[Y'-E(Y')],X'为预测时刻的光伏发电功率,Y'为若干个等间隔时刻点中一时刻点的太阳辐射强度或光伏发电功率,ρX'Y'表示X'和Y'的相关系数,cov(X',Y')表示X'和Y'的协方差,E(X')和E(Y')分别表示X'和Y'的期望,D(X')和D(Y')分别表示X'和Y'的方差;
提取所述若干个等间隔时刻点中,太阳辐射强度和光伏发电功率与所述预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻,得到所述发电功率预测相关时刻。
在其中一个实施例中,所述提取所述历史数据中所述辐射强度预测相关时刻的太阳辐射强度,并进行灰色关联度分析,得到太阳辐射强度预测样本的步骤,包括以下步骤:
获取预测日所述辐射强度预测相关时刻的太阳辐射强度,得到参考序列;
获取预测日之前所述辐射强度预测相关时刻的太阳辐射强度,并对同一日的太阳辐射强度建立序列,得到多个比较序列;
根据对所述参考序列和比较序列进行无量纲化处理,得到无量纲化参考序列和无量纲化比较序列,其中表示序列Xi(k)中的最大值,且i=0时X0(k)为参考序列,i不为0时Xi(k)为比较序列,xi(k)表示序列Xi(k)对应的无量纲化序列;
根据 计算所述无量纲化参考序列和无量纲化比较序列的灰色关联系数,其中x0(k)为无量纲化参考序列,xi(k)为无量纲化比较序列,ρ为分辨系数,ξ0i(k)表示无量纲化参考序列x0(k)中第k个值,与无量纲化比较序列xi(k)中第k个值的灰色关联系数;
根据计算所述比较序列与参考序列的关联度,其中ξ0i(k)表示无量纲化参考序列x0(k)中第k个值,与无量纲化比较序列xi(k)中第k个值的灰色关联系数,r0i表示参考序列X0(k)与比较序列Xi(k)的关联度;
根据比较序列与参考序列关联度的大小,对比较序列进行排序,提取前预设个数的比较序列,得到所述太阳辐射强度预测样本。
在其中一个实施例中,所述采用BP神经网络对所述太阳辐射强度预测样本进行训练,得到太阳辐射强度预测模型的步骤,包括以下步骤:
初始化输入层与隐含层的连接权值、隐含层与输出层的连接权值、隐含层各节点的阈值以及输出层各节点的阈值;
获取所述太阳辐射强度预测样本中多日的辐射强度预测相关时刻的太阳辐射强度,并将同一天的所述辐射强度预测相关时刻的太阳辐射强度作为一个学习向量,得到学习向量集;
将所述学习向量集中的一个学习向量作为所述输入层的输入,计算所述隐含层的输出向量;具体为
根据所述隐含层的输出向量计算所述输出层的输出向量;具体为
f2(x)=purelin(x)=x
根据所述输出层各节点的校正误差,对隐含层与输出层的连接权值,及输出层各节点的阈值进行修正;具体为
其中α为学习速率,Δvjt表示对所述隐含层与输出层的连接权值vjt的修正值,Δγt表示对所述输出层各节点的阈值γt的修正值;
根据所述隐含层各节点的校正误差,对输入层与隐含层的连接权值,及隐含层各节点的阈值进行修正;具体为
其中β为学习速率,Δwij表示对所述输入层与隐含层的连接权值wij的修正值,Δθj表示对所述隐含层各节点的阈值θj的修正值;
判断所述学习向量集中的学习向量是否训练完毕;
若否,则返回所述将所述学习向量集中的一个学习向量作为所述输入层的输入,计算所述隐含层的输出向量的步骤;
若是,则判断全局误差是否小于误差阈值;若否,则返回所述获取所述太阳辐射强度预测样本中多日的辐射强度预测相关时刻的太阳辐射强度,并将同一天的所述辐射强度预测相关时刻的太阳辐射强度作为一个学习向量,得到学习向量集的步骤,若是,则得到所述太阳辐射强度预测模型。
一种光伏发电功率预测***,包括:
相关时刻计算模块,用于根据历史数据,采用相关性分析方法获取辐射强度预测相关时刻和发电功率预测相关时刻,所述历史数据包括太阳辐射强度、温度和光伏发电功率,所述辐射强度预测相关时刻指太阳辐射强度与预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻,所述发电功率预测相关时刻指太阳辐射强度和光伏发电功率与所述预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻;
辐射强度预测样本建立模块,用于提取所述历史数据中所述辐射强度预测相关时刻的太阳辐射强度,并进行灰色关联度分析,得到太阳辐射强度预测样本;
辐射强度预测模型建立模块,用于采用BP神经网络对所述太阳辐射强度预测样本进行训练,得到太阳辐射强度预测模型;
预测时刻辐射强度计算模块,用于将预测日所述辐射强度预测相关时刻的太阳辐射强度作为所述太阳辐射强度预测模型的输入,计算预测日所述预测时刻的太阳辐射强度;
发电功率预测样本建立模块,用于提取所述历史数据中预测日之前所述发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及所述预测时刻的温度和太阳辐射强度,得到光伏发电功率预测样本;
发电功率预测模型建立模块,用于采用BP神经网络对所述光伏发电功率预测样本进行训练,得到光伏发电功率预测模型;
预测时刻发电功率计算模块,用于将预测日所述发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及所述预测时刻的太阳辐射强度和温度作为所述光伏发电功率预测模型的输入,计算预测日所述预测时刻的光伏发电功率。
在其中一个实施例中,所述相关时刻计算模块包括:
第一提取单元,用于提取所述历史数据中所述预测日之前,所述预测时刻及预测时刻前设定时间段内的太阳辐射强度;
第一计算单元,用于根据分别计算所述设定时间段内各时刻点的太阳辐射强度,与所述预测时刻的太阳辐射强度的相关系数,其中,cov(X,Y)=E[X-E(X)]E[Y-E(Y)],X为预测时刻的太阳辐射强度,Y为设定时间段内一时刻点的太阳辐射强度,ρXY表示X和Y的相关系数,cov(X,Y)表示X和Y的协方差,E(X)和E(Y)分别表示X和Y的期望,D(X)和D(Y)分别表示X和Y的方差;
第一处理单元,用于提取所述设定时间段内,太阳辐射强度与所述预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻,得到所述辐射强度预测相关时刻;
第二提取单元,用于提取所述历史数据中所述预测日之前,所述预测时刻的光伏发电功率及预测时刻前设定时间段内的太阳辐射强度和光伏发电功率;
第二计算单元,用于根据分别计算所述设定时间段内各时刻点的太阳辐射强度和光伏发电功率,与所述预测时刻的光伏发电功率的相关系数,其中cov(X',Y')=E[X'-E(X')]E[Y'-E(Y')],X'为预测时刻的光伏发电功率,Y'为设定时间段内一时刻点的太阳辐射强度或光伏发电功率,ρX'Y'表示X'和Y'的相关系数,cov(X',Y')表示X'和Y'的协方差,E(X')和E(Y')分别表示X'和Y'的期望,D(X')和D(Y')分别表示X'和Y'的方差;
第二处理单元,用于提取设定时间范围内,太阳辐射强度和光伏发电功率与所述预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻,得到所述发电功率预测相关时刻。
在其中一个实施例中,所述相关时刻计算模块包括:
第一提取单元,用于提取所述历史数据中所述预测日之前,所述预测时刻及预测时刻前若干个等间隔时刻点的太阳辐射强度;
第一计算单元,用于根据分别计算所述若干个等间隔时刻点的太阳辐射强度,与所述预测时刻的太阳辐射强度的相关系数,其中,cov(X,Y)=E[X-E(X)]E[Y-E(Y)],X为预测时刻的太阳辐射强度,Y为若干个等间隔时刻点中一时刻点的太阳辐射强度,ρXY表示X和Y的相关系数,cov(X,Y)表示X和Y的协方差,E(X)和E(Y)分别表示X和Y的期望,D(X)和D(Y)分别表示X和Y的方差;
第一处理单元,用于提取所述若干个等间隔时刻点中,太阳辐射强度与所述预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻点,得到所述辐射强度预测相关时刻;
第二提取单元,用于提取所述历史数据中所述预测日之前,所述预测时刻的光伏发电功率及预测时刻前若干个等间隔时刻点的太阳辐射强度和光伏发电功率;
第二计算单元,用于根据分别计算所述若干个等间隔时刻点的太阳辐射强度和光伏发电功率,与所述预测时刻的光伏发电功率的相关系数,其中cov(X',Y')=E[X'-E(X')]E[Y'-E(Y')],X'为预测时刻的光伏发电功率,Y'为若干个等间隔时刻点中一时刻点的太阳辐射强度或光伏发电功率,ρX'Y'表示X'和Y'的相关系数,cov(X',Y')表示X'和Y'的协方差,E(X')和E(Y')分别表示X'和Y'的期望,D(X')和D(Y')分别表示X'和Y'的方差;
第二处理单元,用于提取所述若干个等间隔时刻点中,太阳辐射强度和光伏发电功率与所述预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻,得到所述发电功率预测相关时刻。
在其中一个实施例中,所述辐射强度预测样本建立模块包括:
参考序列建立单元,用于获取预测日所述辐射强度预测相关时刻的太阳辐射强度,得到参考序列;
比较序列建立单元,用于获取预测日之前所述辐射强度预测相关时刻的太阳辐射强度,并对同一日的太阳辐射强度建立序列,得到多个比较序列;
无量纲化处理单元,用于根据对所述参考序列和比较序列进行无量纲化处理,得到无量纲化参考序列和无量纲化比较序列,其中表示序列Xi(k)中的最大值,且i=0时X0(k)为参考序列,i不为0时Xi(k)为比较序列,xi(k)表示序列Xi(k)对应的无量纲化序列;
灰色关联系数计算单元,用于根据 计算所述无量纲化参考序列和无量纲化比较序列的灰色关联系数,其中x0(k)为无量纲化参考序列,xi(k)为无量纲化比较序列,ρ为分辨系数,ξ0i(k)表示无量纲化参考序列x0(k)中第k个值,与无量纲化比较序列xi(k)中第k个值的灰色关联系数;
关联度计算单元,用于根据计算所述比较序列与参考序列的关联度,其中ξ0i(k)表示无量纲化参考序列x0(k)中第k个值,与无量纲化比较序列xi(k)中第k个值的灰色关联系数,r0i表示参考序列X0(k)与比较序列Xi(k)的关联度;
关联度排序单元,用于根据比较序列与参考序列关联度的大小,对比较序列进行排序,提取前预设个数的比较序列,得到所述太阳辐射强度预测样本。
在其中一个实施例中,所述辐射强度预测模型建立模块包括:
初始化单元,用于初始化输入层与隐含层的连接权值、隐含层与输出层的连接权值、隐含层各节点的阈值以及输出层各节点的阈值;
学习向量集建立单元,用于获取所述太阳辐射强度预测样本中多日的辐射强度预测相关时刻的太阳辐射强度,并将同一天的所述辐射强度预测相关时刻的太阳辐射强度作为一个学习向量,得到学习向量集;
隐含层输出向量计算单元,用于将所述学习向量集中的一个学习向量作为所述输入层的输入,计算所述隐含层的输出向量;具体为
输出层输出向量计算单元,用于根据所述隐含层的输出向量计算所述输出层的输出向量;具体为
f2(x)=purelin(x)=x
第一修正单元,用于根据所述输出层各节点的校正误差,对隐含层与输出层的连接权值,及输出层各节点的阈值进行修正;具体为
其中α为学习速率,Δvjt表示对所述隐含层与输出层的连接权值vjt的修正值,Δγt表示对所述输出层各节点的阈值γt的修正值;
第二修正单元,用于根据所述隐含层各节点的校正误差,对输入层与隐含层的连接权值,及隐含层各节点的阈值进行修正;具体为
其中β为学习速率,Δwij表示对所述输入层与隐含层的连接权值wij的修正值,Δθj表示对所述隐含层各节点的阈值θj的修正值;
判断单元,用于判断所述学习向量集中的学习向量是否训练完毕,并在所述学习向量集中的学习向量未训练完毕时,控制所述隐含层输出向量计算单元将所述学习向量集中的一个学习向量作为所述输入层的输入;
控制单元,用于在所述学习向量集中的学习向量训练完毕时,判断全局误差是否小于误差阈值,并在全局误差不小于误差阈值时,控制所述学习向量集建立单元获取所述太阳辐射强度预测样本中多日的辐射强度预测相关时刻的太阳辐射强度,并将同一天的所述辐射强度预测相关时刻的太阳辐射强度作为一个学习向量,得到学习向量集;若全局误差小于误差阈值,则得到所述太阳辐射强度预测模型。
上述光伏发电功率预测方法和***,采用相关性分析方法对历史数据进行分析,确定辐射强度预测相关时刻和发电功率预测相关时刻。采用BP神经网络对太阳辐射强度预测样本和光伏发电功率预测样本进行训练,得到太阳辐射强度预测模型和光伏发电功率预测模型。利用太阳辐射强度预测模型计算预测日预测时刻的太阳辐射强度,利用光伏发电功率预测模型计算预测日预测时刻的光伏发电功率。利用灰色关联度分析方法去除历史数据中关联度较低的辐射强度预测相关时刻的太阳辐射强度,提高了太阳辐射强度预测精度。采用BP神经网络良好的非线性函数逼近能力,对太阳辐射强度预测样本和光伏发电功率预测样本进行训练建立预测模型,提高了预测模型的预测精度,从而提高光伏并网***的稳定性,降低运行成本。
附图说明
图1为一实施例中光伏发电功率预测方法的流程图;
图2为一实施例中光伏发电功率预测***的结构图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图对本发明的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似改进,因此本发明不受下面公开的具体实施例的限制。
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在于限制本发明。本文所使用的术语“及/或”包括一个或多个相关的所列项目的任意的和所有的组合。
一种光伏发电功率预测方法,如图1所示,包括以下步骤:
步骤S110:根据历史数据,采用相关性分析方法获取辐射强度预测相关时刻和发电功率预测相关时刻。
历史数据包括太阳辐射强度、温度和光伏发电功率。具体地,既包括预测日之前各时刻的太阳辐射强度、温度和光伏发电功率,也包括预测日预测时刻之前各时刻的太阳辐射强度和光伏发电功率,以及预测时刻的温度。由于温度是唯一能够直接得到的预测时刻的非历史相关因素,因此也将其作为一个参考量,预测时刻的温度可通过天气预报获取。为便于理解,本发明以预测10月1日10点的光伏发电功率,历史数据包括9月1日至9月30日每天各时刻的太阳辐射强度和光伏发电功率,10月1日10点之前各时刻的太阳辐射强度和光伏发电功率,以及10月1日10点的温度为例,对本发明的具体实施方式进行解释说明。
辐射强度预测相关时刻指太阳辐射强度与预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻。发电功率预测相关时刻指太阳辐射强度和光伏发电功率与预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻。采用相关性分析方法对历史数据进行分析,确定辐射强度预测相关时刻和发电功率预测相关时刻,提高了各相关时刻的获取准确度。
在其中一个实施例中,步骤S110包括步骤111至步骤116。
步骤111:提取历史数据中预测日之前,预测时刻及预测时刻前设定时间段内的太阳辐射强度。
提取9月份每天10点以及10点之前设定时间段内的太阳辐射强度,设定时间段可以是固定的时间段,如2小时或3小时等,以设定时间段为2小时为例,步骤111即是提取9月份每天8至10点的太阳辐射强度。设定时间段也可根据实际情况调整。由于与预测时刻越接近,太阳辐射强度的相关性越高,本实施例中设定时间段为3小时,有效减少不必要的数据,提高后续步骤的处理效率。
步骤112:分别计算设定时间段内各时刻点的太阳辐射强度,与预测时刻的太阳辐射强度的相关系数。
cov(X,Y)=E[X-E(X)]E[Y-E(Y)]
其中,X为预测时刻的太阳辐射强度,Y为设定时间段内一时刻点的太阳辐射强度,ρXY表示X和Y的相关系数,cov(X,Y)表示X和Y的协方差,E(X)和E(Y)分别表示X和Y的期望,D(X)和D(Y)分别表示X和Y的方差。
同样以预测时刻为10月1日10点为例,则X包括9月份每天10点的太阳辐射强度,Y包括9月份每天设定时间段内一时刻点的太阳辐射强度。
步骤113:提取设定时间段内,太阳辐射强度与预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻,得到辐射强度预测相关时刻。
辐射相关阈值可以是根据经验设置的固定值,计算设定时间段内各时刻的太阳辐射强度与预测时刻的太阳辐射强度相关系数后,将相关系数高于辐射相关阈值的太阳辐射强度对应的时刻最为辐射强度预测相关时刻。也可以通过调整辐射相关阈值,使相关系数高于辐射相关阈值的太阳辐射强度对应时刻的数量为设定个数,如2个、3个、4个等。本实施例中提取相关系数最高的3个太阳辐射强度对应的时刻作为辐射强度预测相关时刻。
步骤114:提取历史数据中预测日之前,预测时刻的光伏发电功率及预测时刻前设定时间段内的太阳辐射强度和光伏发电功率。
提取9月份每天10点的光伏发电功率,及10点之前设定时间段内的太阳辐射强度和光伏发电功率,本实施例中设定时间段同样可为3小时。
步骤115:分别计算设定时间段内各时刻点的太阳辐射强度和光伏发电功率,与预测时刻的光伏发电功率的相关系数。
cov(X',Y')=E[X'-E(X')]E[Y'-E(Y')]
其中,X'为预测时刻的光伏发电功率,Y'为设定时间段内一时刻点的太阳辐射强度或光伏发电功率,ρX'Y'表示X'和Y'的相关系数,cov(X',Y')表示X'和Y'的协方差,E(X')和E(Y')分别表示X'和Y'的期望,D(X')和D(Y')分别表示X'和Y'的方差。
同样以预测时刻为10月1日10点为例,则X包括9月份每天10点的光伏发电功率,Y包括9月份每天设定时间段内一时刻点的太阳辐射强度,或者包括9月份每天设定时间段内一时刻点的光伏发电功率。
步骤116:提取设定时间范围内,太阳辐射强度和光伏发电功率与预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻,得到发电功率预测相关时刻。
同理,功率相关阈值可以是根据经验设置的固定值,也可以通过调整功率相关阈值。本实施例中提取相关系数最高的3个太阳辐射强度和光伏发电功率对应的时刻作为发电功率预测相关时刻。
在另一个实施例中,步骤S110包括步骤1至步骤6。
步骤1:提取历史数据中预测日之前,预测时刻及预测时刻前若干个等间隔时刻点的太阳辐射强度。
若干个等间隔时刻点之间的时间间隔可以是0.2小时、0.5小时、0.8小时等,本实施例中时间间隔为0.5小时。以提取预测时刻前5个等间隔时刻点的太阳辐射强度为例,步骤1即是提取9月份每天10点、9点半、9点、8点半、8点和7点半的太阳辐射强度。
步骤2:分别计算若干个等间隔时刻点的太阳辐射强度,与预测时刻的太阳辐射强度的相关系数。
cov(X,Y)=E[X-E(X)]E[Y-E(Y)]
其中,X为预测时刻的太阳辐射强度,Y为若干个等间隔时刻点中一时刻点的太阳辐射强度,ρXY表示X和Y的相关系数,cov(X,Y)表示X和Y的协方差,E(X)和E(Y)分别表示X和Y的期望,D(X)和D(Y)分别表示X和Y的方差。
同样以预测时刻为10月1日10点为例,则X包括9月份每天10点的太阳辐射强度,Y包括9月份每天若干个等间隔时刻点中一时刻点的太阳辐射强度。
步骤3:提取若干个等间隔时刻点中,太阳辐射强度与预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻点,得到辐射强度预测相关时刻。
辐射相关阈值可以是根据经验设置的固定值,也可以通过调整辐射相关阈值,使相关系数高于辐射相关阈值的太阳辐射强度对应时刻的数量为设定个数。本实施例中提取相关系数最高的3个太阳辐射强度对应的时刻作为辐射强度预测相关时刻。通过比较相关系数,以预测时刻前0.5小时、前1.0小时和前1.5小时作为辐射强度预测相关时刻。
步骤4:提取历史数据中预测日之前,预测时刻的光伏发电功率及预测时刻前若干个等间隔时刻点的太阳辐射强度和光伏发电功率。
若干个等间隔时刻点之间的时间间隔可以是0.2小时、0.5小时、0.8小时等,本实施例中时间间隔为0.5小时。以提取预测时刻前5个等间隔时刻点的太阳辐射强度和光伏发电功率为例,步骤4即是提取9月份每天10点的光伏发电功率,以及9点半、9点、8点半、8点和7点半的太阳辐射强度和光伏发电功率。
步骤5:根据分别计算若干个等间隔时刻点的太阳辐射强度和光伏发电功率,与预测时刻的光伏发电功率的相关系数。
cov(X',Y')=E[X'-E(X')]E[Y'-E(Y')]
其中,X'为预测时刻的光伏发电功率,Y'为若干个等间隔时刻点中一时刻点的太阳辐射强度或光伏发电功率,ρX'Y'表示X'和Y'的相关系数,cov(X',Y')表示X'和Y'的协方差,E(X')和E(Y')分别表示X'和Y'的期望,D(X')和D(Y')分别表示X'和Y'的方差。
同样以预测时刻为10月1日10点为例,则X包括9月份每天10点的光伏发电功率,Y包括9月份每天若干个等间隔时刻点中一时刻点的太阳辐射强度,或者包括9月份每天若干个等间隔时刻点中一时刻点的光伏发电功率。
步骤6:提取若干个等间隔时刻点中,太阳辐射强度和光伏发电功率与预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻,得到发电功率预测相关时刻。
同理,功率相关阈值可以是根据经验设置的固定值,也可以通过调整功率相关阈值。本实施例中提取相关系数最高的3个太阳辐射强度和光伏发电功率对应的时刻作为发电功率预测相关时刻。通过比较相关系数,以预测时刻前0.5小时、前1.0小时和前1.5小时作为发电功率预测相关时刻。
以上即是对步骤S110提供了两种具体实施方式,在采用相关性分析方法分析历史数据时,可以是提取预测时刻前设定时间段内的全部相关数据,确保计算准确度;也可以是提取预测时刻前若干个等间隔时刻点的相关数据,提高处理速度。在实际应用场景中可根据情况选择适合的方式,提高了光伏发电功率预测方法的适用性。
步骤S120:提取历史数据中辐射强度预测相关时刻的太阳辐射强度,并进行灰色关联度分析,得到太阳辐射强度预测样本。
以预测时刻为10月1日10点、辐射强度预测相关时刻为预测时刻前0.5小时、前1.0小时和前1.5小时为例,步骤S120即是提取9月份和10月1日所有8点半、9点和9点半的太阳辐射强度,然后进行灰色关联度分析,建立太阳辐射强度预测样本。利用灰色关联度分析方法去除历史数据中关联度较低的辐射强度预测相关时刻的太阳辐射强度,提高了太阳辐射强度预测精度。
步骤S120具体可包括步骤121至步骤126。
步骤121:获取预测日辐射强度预测相关时刻的太阳辐射强度,得到参考序列。
在本实施例中,参考序列即是包括10月1日8点半、9点和9点半的太阳辐射强度。
步骤122:获取预测日之前辐射强度预测相关时刻的太阳辐射强度,并对同一日的太阳辐射强度建立序列,得到多个比较序列。
获取9月份每天8点半、9点和9点半的太阳辐射强度,并将同一天的太阳辐射强度按时间顺序排成一个序列,本实施例中比较序列的数量为30个。
步骤123:对参考序列和比较序列进行无量纲化处理,得到无量纲化参考序列和无量纲化比较序列。
在进行灰色关联度分析时,一般都要对数据进行无量纲化处理,由于不同信息特征间数量级差别可能较大,因此在对数据进行无量纲化处理时,是将该信息特征值除以这个信息特征所有值的最大值。
其中表示序列Xi(k)中的最大值,且i=0时X0(k)为参考序列,i不为0时Xi(k)为比较序列,xi(k)表示序列Xi(k)对应的无量纲化序列。n为比较序列的个数,本实施例中为30个,m为参考序列或比较序列中的太阳辐射强度值的个数,本实施例中为3个。
步骤124:计算无量纲化参考序列和无量纲化比较序列的灰色关联系数。
其中x0(k)为无量纲化参考序列,xi(k)为无量纲化比较序列,ρ为分辨系数,本实施例中ρ=0.5,ξ0i(k)表示无量纲化参考序列x0(k)中第k个值,与无量纲化比较序列xi(k)中第k个值的灰色关联系数。以i=1,k=2为例,ξ01(2)即为无量纲化参考序列中的第2个值,与第1个无量纲化比较序列中的第2个值的灰色关联系数。
步骤125:计算比较序列与参考序列的关联度。
其中ξ0i(k)表示无量纲化参考序列x0(k)中第k个值,与无量纲化比较序列xi(k)中第k个值的灰色关联系数,r0i表示参考序列X0(k)与比较序列Xi(k)的关联度。
步骤126:根据比较序列与参考序列关联度的大小,对比较序列进行排序,提取前预设个数的比较序列,得到太阳辐射强度预测样本。
预设个数可以是5至15个,也可根据实际情况进行调整。本实施例中预设个数为7个,即将30个比较序列按与参考序列关联度的大小进行排列后,提取前7个比较序列,组成太阳辐射强度预测样本。
步骤S130:采用BP神经网络对太阳辐射强度预测样本进行训练,得到太阳辐射强度预测模型。
步骤S130具体可包括如下步骤1至步骤11。
步骤1:计算隐含层节点数。
其中l为隐含层节点数,n为输入层节点数,m为输出层节点数,a为调节常数。本实施例中太阳辐射强度预测模型的输入为预测时刻前0.5小时、前1.0小时和前1.5小时的太阳辐射强度,输出为预测时刻的太阳辐射强度,故n为3,m为1,取a=1,因此隐含层节点数l=3。
步骤2:初始化输入层与隐含层的连接权值、隐含层与输出层的连接权值、隐含层各节点的阈值以及输出层各节点的阈值。
步骤3:获取太阳辐射强度预测样本中多日的辐射强度预测相关时刻的太阳辐射强度,并将同一天的辐射强度预测相关时刻的太阳辐射强度作为一个学习向量,得到学习向量集。
为便于理解可结合步骤S120的具体实施例进行说明。本步骤即是提取太阳辐射强度预测样本中的若干个比较序列,如3个、4个或5个等,将提取的每一个比较序列作为一个学习向量,得到学习向量集。
步骤4:将学习向量集中的一个学习向量作为输入层的输入,计算隐含层的输出向量。具体为
输入层各节点对输入不进行处理,只是简单地把接收到的输入向量作为相应的输出向量传递给隐含层,即输入层的输出向量与输入向量相同。隐含层的神经元激励函数采用S型正切函数f1(x)。输入层的输入向量m为输入向量的个数。n为输入层节点数,本实施例中为3,隐含层的输入向量隐含层的输出向量p为隐含层节点数,本实施例中为3。wij为输入层与隐含层的连接权值,θj为隐含层各节点的阈值。
步骤5:根据隐含层的输出向量计算输出层的输出向量。具体为
f2(x)=purelin(x)=x
输出层的神经元激励函数采用纯线性函数f2(x)。隐含层的输出向量输出层的输入向量输出层的输出向量q为输出层节点数,本实施例中为1。vjt为隐含层与输出层的连接权值,γt为输出层各节点的阈值。本实施例中输出层的输出向量即是对预测时刻的太阳辐射强度的预测值。
步骤6:计算输出层各节点的校正误差。
其中表示输出层各节点的校正误差,为与对应的实际太阳辐射强度,表示对求导。举例说明,学习向量集中包括9月15日、18日和22日的8点半、9点和9点半的太阳辐射强度,是将9月15日的8点半、9点和9点半的太阳辐射强度作为输入层的输入而得到的9月15日10点太阳辐射强度的预测值。则取9月15日10点太阳辐射强度的实际值,通过上式计算输出层各节点的校正误差。
步骤7:计算隐含层各节点的校正误差。
步骤8:根据输出层各节点的校正误差,对隐含层与输出层的连接权值,及输出层各节点的阈值进行修正。具体为
其中α为学习速率,Δvjt表示对隐含层与输出层的连接权值vjt的修正值,Δγt表示对输出层各节点的阈值γt的修正值。根据步骤6得到的输出层各节点的校正误差,对隐含层与输出层的连接权值,及输出层各节点的阈值进行修正。
步骤9:根据隐含层各节点的校正误差,对输入层与隐含层的连接权值,及隐含层各节点的阈值进行修正。具体为
其中β为学习速率,Δwij表示对输入层与隐含层的连接权值wij的修正值,Δθj表示对隐含层各节点的阈值θj的修正值。根据步骤7得到的隐含层各节点的校正误差,对输入层与隐含层的连接权值,及隐含层各节点的阈值进行修正。
步骤10:判断学习向量集中的学习向量是否训练完毕。若否,则返回上述步骤4,将学习向量集中还未进行训练的学习向量代入输入层。若是,则进行步骤11。
步骤11:判断全局误差是否小于误差阈值。
将学习向量集中的所有学习向量都代入输入层进行训练后,对模型中的各连接权值和节点阈值进行多次修正。根据修正后的参数再次代入各学习向量计算预测时刻的太阳辐射强度的预测值,将预测值与实际值比较,计算全局误差。若全局误差不小于误差阈值,返回步骤3,在太阳辐射强度预测样本中获取其他日期的辐射强度预测相关时刻的太阳辐射强度,组建新的学习向量集,并将新获取的学习向量代入输入层进行训练,对模型中的各连接权值和节点阈值进行修正。
若全局误差小于误差阈值,则得到太阳辐射强度预测模型。
步骤S140:将预测日辐射强度预测相关时刻的太阳辐射强度作为太阳辐射强度预测模型的输入,计算预测日预测时刻的太阳辐射强度。
本实施例中即是将10月1日的8点半、9点和9点半的太阳辐射强度作为太阳辐射强度预测模型的输入,预测10月1日10点的太阳辐射强度。
由于步骤S120建立的太阳辐射强度预测样本只包括前预设个数的比较序列,因此本实施例中在每次对太阳辐射强度进行预测时,都重新建立太阳辐射强度预测样本和太阳辐射强度预测模型,进一步提高预测准确度。
步骤S150:提取历史数据中预测日之前发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及预测时刻的温度和太阳辐射强度,得到光伏发电功率预测样本。
提取9月每天8点半、9点和9点半的太阳辐射强度和光伏发电功率,以及10点的温度和太阳辐射强度,得到光伏发电功率预测样本。
步骤S160:采用BP神经网络对光伏发电功率预测样本进行训练,得到光伏发电功率预测模型。
同样也是先计算隐含层节点数。本实施例中光伏发电功率预测模型的输入为预测时刻前0.5小时、前1.0小时和前1.5小时的太阳辐射强度和光伏发电功率,以及预测时刻的温度和太阳辐射强度,输出为预测时刻的光伏发电功率。故输入层节点数为8,输出层节点数为1,取调节常数a=2,则根据可计算出隐含层节点数为5。
然后对输入层与隐含层的连接权值、隐含层与输出层的连接权值、隐含层各节点的阈值以及输出层各节点的阈值进行初始化。获取光伏发电功率预测样本中多日的发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及预测时刻的温度和太阳辐射强度,将同一日的相关数据作为一个学习向量,得到学习向量集。
将学习向量代入输入层进行训练,隐含层和输出层采用的神经元激励函数与步骤S130中的相同。计算输出层各节点的校正误差,隐含层各节点的校正误差,对隐含层与输出层的连接权值,及输出层各节点的阈值进行修正,对输入层与隐含层的连接权值,及隐含层各节点的阈值进行修正等步骤的原理也与步骤S130中类似,不再赘述。
步骤S170:将预测日发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及预测时刻的太阳辐射强度和温度作为光伏发电功率预测模型的输入,计算预测日预测时刻的光伏发电功率。
本实施例中将10月1日8点半、9点和9点半的太阳辐射强度和光伏发电功率,以及10点的太阳辐射强度和温度作为光伏发电功率预测模型的输入,预测10月1日10点的光伏发电功率。
由于步骤S150建立的光伏发电功率预测样本包括了历史数据中,预测日之前所有发电功率预测相关时刻的太阳辐射强度和光伏发电功率,数据比较多,能够保证光伏发电功率预测模型的预测准确度。本实施例中在光伏发电功率预测模型建好之后便固定不变,即是模型建好之后每次对光伏发电功率进行预测时,直接将预测日发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及预测时刻的太阳辐射强度和温度作为光伏发电功率预测模型的输入来进行预测,可节省时间。可以理解,在其他实施例中,也可以每次对光伏发电功率进行预测时,都重新建立光伏发电功率预测模型。
上述光伏发电功率预测方法,采用相关性分析方法对历史数据进行分析,确定辐射强度预测相关时刻和发电功率预测相关时刻。采用BP神经网络对太阳辐射强度预测样本和光伏发电功率预测样本进行训练,得到太阳辐射强度预测模型和光伏发电功率预测模型。利用太阳辐射强度预测模型计算预测日预测时刻的太阳辐射强度,利用光伏发电功率预测模型计算预测日预测时刻的光伏发电功率。利用灰色关联度分析方法去除历史数据中关联度较低的辐射强度预测相关时刻的太阳辐射强度,提高了太阳辐射强度预测精度。采用BP神经网络良好的非线性函数逼近能力,对太阳辐射强度预测样本和光伏发电功率预测样本进行训练建立预测模型,提高了预测模型的预测精度,从而提高光伏并网***的稳定性,降低运行成本。
本发明还提供了一种光伏发电功率预测***,如图2所示,包括相关时刻计算模块110、辐射强度预测样本建立模块120、辐射强度预测模型建立模块130、预测时刻辐射强度计算模块140、发电功率预测样本建立模块150、发电功率预测模型建立模块160和预测时刻发电功率计算模块170。
相关时刻计算模块110用于根据历史数据,采用相关性分析方法获取辐射强度预测相关时刻和发电功率预测相关时刻。
历史数据包括太阳辐射强度、温度和光伏发电功率,具体地,既包括预测日之前各时刻的太阳辐射强度、温度和光伏发电功率,也包括预测日预测时刻之前各时刻的太阳辐射强度和光伏发电功率,以及预测时刻的温度。由于温度是唯一能够直接得到的预测时刻的非历史相关因素,因此也将其作为一个参考量,预测时刻的温度可通过天气预报获取。为便于理解,本发明以预测10月1日10点的光伏发电功率,历史数据包括9月1日至9月30日每天各时刻的太阳辐射强度和光伏发电功率,10月1日10点之前各时刻的太阳辐射强度和光伏发电功率,以及10月1日10点的温度为例,对本发明的具体实施方式进行解释说明。
辐射强度预测相关时刻指太阳辐射强度与预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻,发电功率预测相关时刻指太阳辐射强度和光伏发电功率与预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻。采用相关性分析方法对历史数据进行分析,确定辐射强度预测相关时刻和发电功率预测相关时刻,提高了各相关时刻的获取准确度。
在其中一个实施例中,相关时刻计算模块110包括如下的第一提取单元、第一计算单元、第一处理单元、第二提取单元、第二计算单元和第二处理单元。
第一提取单元用于提取历史数据中预测日之前,预测时刻及预测时刻前设定时间段内的太阳辐射强度。
提取9月份每天10点以及10点之前设定时间段内的太阳辐射强度,设定时间段可以是固定的时间段,如2小时或3小时等,以设定时间段为2小时为例,第一提取单元即是提取9月份每天8至10点的太阳辐射强度。设定时间段也可根据实际情况调整。由于与预测时刻越接近,太阳辐射强度的相关性越高,本实施例中设定时间段为3小时,有效减少不必要的数据,提高后续步骤的处理效率。
第一计算单元用于分别计算设定时间段内各时刻点的太阳辐射强度,与预测时刻的太阳辐射强度的相关系数。
cov(X,Y)=E[X-E(X)]E[Y-E(Y)]
其中,X为预测时刻的太阳辐射强度,Y为设定时间段内一时刻点的太阳辐射强度,ρXY表示X和Y的相关系数,cov(X,Y)表示X和Y的协方差,E(X)和E(Y)分别表示X和Y的期望,D(X)和D(Y)分别表示X和Y的方差。
同样以预测时刻为10月1日10点为例,则X包括9月份每天10点的太阳辐射强度,Y包括9月份每天设定时间段内一时刻点的太阳辐射强度。
第一处理单元用于提取设定时间段内,太阳辐射强度与预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻,得到辐射强度预测相关时刻。
辐射相关阈值可以是根据经验设置的固定值,计算设定时间段内各时刻的太阳辐射强度与预测时刻的太阳辐射强度相关系数后,将相关系数高于辐射相关阈值的太阳辐射强度对应的时刻最为辐射强度预测相关时刻。也可以通过调整辐射相关阈值,使相关系数高于辐射相关阈值的太阳辐射强度对应时刻的数量为设定个数,如2个、3个、4个等。本实施例中提取相关系数最高的3个太阳辐射强度对应的时刻作为辐射强度预测相关时刻。
第二提取单元用于提取历史数据中预测日之前,预测时刻的光伏发电功率及预测时刻前设定时间段内的太阳辐射强度和光伏发电功率。
提取9月份每天10点的光伏发电功率,及10点之前设定时间段内的太阳辐射强度和光伏发电功率,本实施例中设定时间段同样可为3小时。
第二计算单元用于分别计算设定时间段内各时刻点的太阳辐射强度和光伏发电功率,与预测时刻的光伏发电功率的相关系数。
cov(X',Y')=E[X'-E(X')]E[Y'-E(Y')]
其中,X'为预测时刻的光伏发电功率,Y'为设定时间段内一时刻点的太阳辐射强度或光伏发电功率,ρX'Y'表示X'和Y'的相关系数,cov(X',Y')表示X'和Y'的协方差,E(X')和E(Y')分别表示X'和Y'的期望,D(X')和D(Y')分别表示X'和Y'的方差。
同样以预测时刻为10月1日10点为例,则X包括9月份每天10点的光伏发电功率,Y包括9月份每天设定时间段内一时刻点的太阳辐射强度,或者包括9月份每天设定时间段内一时刻点的光伏发电功率。
第二处理单元用于提取设定时间范围内,太阳辐射强度和光伏发电功率与预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻,得到发电功率预测相关时刻。
同理,功率相关阈值可以是根据经验设置的固定值,也可以通过调整功率相关阈值。本实施例中提取相关系数最高的3个太阳辐射强度和光伏发电功率对应的时刻作为发电功率预测相关时刻。
在另一个实施例中,相关时刻计算模块110包括如下的第一提取单元、第一计算单元、第一处理单元、第二提取单元、第二计算单元和第二处理单元。
第一提取单元用于提取历史数据中预测日之前,预测时刻及预测时刻前若干个等间隔时刻点的太阳辐射强度。
若干个等间隔时刻点之间的时间间隔可以是0.2小时、0.5小时、0.8小时等,本实施例中时间间隔为0.5小时。以提取预测时刻前5个等间隔时刻点的太阳辐射强度为例,第一提取单元即是提取9月份每天10点、9点半、9点、8点半、8点和7点半的太阳辐射强度。
第一计算单元用于分别计算若干个等间隔时刻点的太阳辐射强度,与预测时刻的太阳辐射强度的相关系数。
cov(X,Y)=E[X-E(X)]E[Y-E(Y)]
其中,X为预测时刻的太阳辐射强度,Y为若干个等间隔时刻点中一时刻点的太阳辐射强度,ρXY表示X和Y的相关系数,cov(X,Y)表示X和Y的协方差,E(X)和E(Y)分别表示X和Y的期望,D(X)和D(Y)分别表示X和Y的方差。
同样以预测时刻为10月1日10点为例,则X包括9月份每天10点的太阳辐射强度,Y包括9月份每天若干个等间隔时刻点中一时刻点的太阳辐射强度。
第一处理单元用于提取若干个等间隔时刻点中,太阳辐射强度与预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻点,得到辐射强度预测相关时刻。
辐射相关阈值可以是根据经验设置的固定值,也可以通过调整辐射相关阈值,使相关系数高于辐射相关阈值的太阳辐射强度对应时刻的数量为设定个数。本实施例中提取相关系数最高的3个太阳辐射强度对应的时刻作为辐射强度预测相关时刻。通过比较相关系数,以预测时刻前0.5小时、前1.0小时和前1.5小时作为辐射强度预测相关时刻。
第二提取单元用于提取历史数据中预测日之前,预测时刻的光伏发电功率及预测时刻前若干个等间隔时刻点的太阳辐射强度和光伏发电功率。
若干个等间隔时刻点之间的时间间隔可以是0.2小时、0.5小时、0.8小时等,本实施例中时间间隔为0.5小时。以提取预测时刻前5个等间隔时刻点的太阳辐射强度和光伏发电功率为例,第二提取单元即是提取9月份每天10点的光伏发电功率,以及9点半、9点、8点半、8点和7点半的太阳辐射强度和光伏发电功率。
第二计算单元用于分别计算若干个等间隔时刻点的太阳辐射强度和光伏发电功率,与预测时刻的光伏发电功率的相关系数。
cov(X',Y')=E[X'-E(X')]E[Y'-E(Y')]
其中,X'为预测时刻的光伏发电功率,Y'为若干个等间隔时刻点中一时刻点的太阳辐射强度或光伏发电功率,ρX'Y'表示X'和Y'的相关系数,cov(X',Y')表示X'和Y'的协方差,E(X')和E(Y')分别表示X'和Y'的期望,D(X')和D(Y')分别表示X'和Y'的方差。
同样以预测时刻为10月1日10点为例,则X包括9月份每天10点的光伏发电功率,Y包括9月份每天若干个等间隔时刻点中一时刻点的太阳辐射强度,或者包括9月份每天若干个等间隔时刻点中一时刻点的光伏发电功率。
第二处理单元用于提取若干个等间隔时刻点中,太阳辐射强度和光伏发电功率与预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻,得到发电功率预测相关时刻。
同理,功率相关阈值可以是根据经验设置的固定值,也可以通过调整功率相关阈值。本实施例中提取相关系数最高的3个太阳辐射强度和光伏发电功率对应的时刻作为发电功率预测相关时刻。通过比较相关系数,以预测时刻前0.5小时、前1.0小时和前1.5小时作为发电功率预测相关时刻。
以上即是在采用相关性分析方法分析历史数据时,提供了两种具体实施方式,可以是提取预测时刻前设定时间段内的全部相关数据,确保计算准确度;也可以是提取预测时刻前若干个等间隔时刻点的相关数据,提高处理速度。在实际应用场景中可根据情况选择适合的方式,提高了光伏发电功率预测***的适用性。
辐射强度预测样本建立模块120用于提取历史数据中辐射强度预测相关时刻的太阳辐射强度,并进行灰色关联度分析,得到太阳辐射强度预测样本。
以预测时刻为10月1日10点、辐射强度预测相关时刻为预测时刻前0.5小时、前1.0小时和前1.5小时为例,辐射强度预测样本建立模块120即是提取9月份和10月1日所有8点半、9点和9点半的太阳辐射强度,然后进行灰色关联度分析,建立太阳辐射强度预测样本。利用灰色关联度分析方法去除历史数据中关联度较低的辐射强度预测相关时刻的太阳辐射强度,提高了太阳辐射强度预测精度。
辐射强度预测样本建立模块120具体可包括参考序列建立单元、比较序列建立单元、无量纲化处理单元、灰色关联系数计算单元、关联度计算单元和关联度排序单元。
参考序列建立单元用于获取预测日辐射强度预测相关时刻的太阳辐射强度,得到参考序列。
在本实施例中,参考序列即是包括10月1日8点半、9点和9点半的太阳辐射强度。
比较序列建立单元用于获取预测日之前辐射强度预测相关时刻的太阳辐射强度,并对同一日的太阳辐射强度建立序列,得到多个比较序列。
获取9月份每天8点半、9点和9点半的太阳辐射强度,并将同一天的太阳辐射强度按时间顺序排成一个序列,本实施例中比较序列的数量为30个。
无量纲化处理单元用于对参考序列和比较序列进行无量纲化处理,得到无量纲化参考序列和无量纲化比较序列。
在进行灰色关联度分析时,一般都要对数据进行无量纲化处理,由于不同信息特征间数量级差别可能较大,因此在对数据进行无量纲化处理时,是将该信息特征值除以这个信息特征所有值的最大值。
其中表示序列Xi(k)中的最大值,且i=0时X0(k)为参考序列,i不为0时Xi(k)为比较序列,xi(k)表示序列Xi(k)对应的无量纲化序列。n为比较序列的个数,本实施例中为30个,m为参考序列或比较序列中的太阳辐射强度值的个数,本实施例中为3个。
灰色关联系数计算单元用于计算无量纲化参考序列和无量纲化比较序列的灰色关联系数。
其中x0(k)为无量纲化参考序列,xi(k)为无量纲化比较序列,ρ为分辨系数,本实施例中ρ=0.5,ξ0i(k)表示无量纲化参考序列x0(k)中第k个值,与无量纲化比较序列xi(k)中第k个值的灰色关联系数。以i=1,k=2为例,ξ01(2)即为无量纲化参考序列中的第2个值,与第1个无量纲化比较序列中的第2个值的灰色关联系数。
关联度计算单元用于计算比较序列与参考序列的关联度。
其中ξ0i(k)表示无量纲化参考序列x0(k)中第k个值,与无量纲化比较序列xi(k)中第k个值的灰色关联系数,r0i表示参考序列X0(k)与比较序列Xi(k)的关联度。
关联度排序单元用于根据比较序列与参考序列关联度的大小,对比较序列进行排序,提取前预设个数的比较序列,得到太阳辐射强度预测样本。
预设个数可以是5至15个,也可根据实际情况进行调整。本实施例中预设个数为7个,即将30个比较序列按与参考序列关联度的大小进行排列后,提取前7个比较序列,组成太阳辐射强度预测样本。
辐射强度预测模型建立模块130用于采用BP神经网络对太阳辐射强度预测样本进行训练,得到太阳辐射强度预测模型。
辐射强度预测模型建立模块130具体包括隐含层节点计算单元、初始化单元、学习向量集建立单元、隐含层输出向量计算单元、输出层输出向量计算单元、输出层节点校正误差计算单元、隐含层节点校正误差计算单元、第一修正单元、第二修正单元、判断单元和控制单元。
隐含层节点计算单元用于计算隐含层节点数。
其中l为隐含层节点数,n为输入层节点数,m为输出层节点数,a为调节常数。本实施例中太阳辐射强度预测模型的输入为预测时刻前0.5小时、前1.0小时和前1.5小时的太阳辐射强度,输出为预测时刻的太阳辐射强度,故n为3,m为1,取a=1,因此隐含层节点数l=3。
初始化单元用于初始化输入层与隐含层的连接权值、隐含层与输出层的连接权值、隐含层各节点的阈值以及输出层各节点的阈值。
学习向量集建立单元用于获取太阳辐射强度预测样本中多日的辐射强度预测相关时刻的太阳辐射强度,并将同一天的辐射强度预测相关时刻的太阳辐射强度作为一个学习向量,得到学习向量集。
本步骤即是提取太阳辐射强度预测样本中的若干个比较序列,如3个、4个或5个等,将提取的每一个比较序列作为一个学习向量,得到学习向量集。
隐含层输出向量计算单元用于将学习向量集中的一个学习向量作为输入层的输入,计算隐含层的输出向量。具体为
输入层各节点对输入不进行处理,只是简单地把接收到的输入向量作为相应的输出向量传递给隐含层,即输入层的输出向量与输入向量相同。隐含层的神经元激励函数采用S型正切函数f1(x)。输入层的输入向量m为输入向量的个数,n为输入层节点数,本实施例中为3,隐含层的输入向量隐含层的输出向量p为隐含层节点数,本实施例中为3。wij为输入层与隐含层的连接权值,θj为隐含层各节点的阈值。
输出层输出向量计算单元用于根据隐含层的输出向量计算输出层的输出向量。具体为
f2(x)=purelin(x)=x
输出层节点校正误差计算单元用于计算输出层各节点的校正误差。
其中表示输出层各节点的校正误差,为与对应的实际太阳辐射强度,表示对求导。举例说明,学习向量集中包括9月15日、18日和22日的8点半、9点和9点半的太阳辐射强度,是将9月15日的8点半、9点和9点半的太阳辐射强度作为输入层的输入而得到的9月15日10点太阳辐射强度的预测值。则取9月15日10点太阳辐射强度的实际值,通过上式计算输出层各节点的校正误差。
隐含层节点校正误差计算单用于计算隐含层各节点的校正误差。
第一修正单元用于根据输出层各节点的校正误差,对隐含层与输出层的连接权值,及输出层各节点的阈值进行修正。具体为
其中α为学习速率,Δvjt表示对隐含层与输出层的连接权值vjt的修正值,Δγt表示对输出层各节点的阈值γt的修正值。根据输出层节点校正误差计算单元得到的输出层各节点的校正误差,对隐含层与输出层的连接权值,及输出层各节点的阈值进行修正。
第二修正单元用于根据隐含层各节点的校正误差,对输入层与隐含层的连接权值,及隐含层各节点的阈值进行修正。具体为
其中β为学习速率,Δwij表示对输入层与隐含层的连接权值wij的修正值,Δθj表示对隐含层各节点的阈值θj的修正值。根据隐含层节点校正误差计算单得到的隐含层各节点的校正误差,对输入层与隐含层的连接权值,及隐含层各节点的阈值进行修正。
判断单元用于判断学习向量集中的学习向量是否训练完毕,并在学习向量集中的学习向量未训练完毕时,控制隐含层输出向量计算单元将学习向量集中还未进行训练的学习向量代入输入层。
控制单元用于在学习向量集中的学习向量训练完毕时,判断全局误差是否小于误差阈值。
将学习向量集中的所有学习向量都代入输入层进行训练后,对模型中的各连接权值和节点阈值进行多次修正。根据修正后的参数再次代入各学习向量计算预测时刻的太阳辐射强度的预测值,将预测值与实际值比较,计算全局误差。控制单元在全局误差不小于误差阈值时,控制学习向量集建立单元在太阳辐射强度预测样本中获取其他日期的发电功率预测相关时刻的太阳辐射强度,组建新的学习向量集,并将新获取的学习向量代入输入层进行训练,对模型中的各连接权值和节点阈值进行修正。
若全局误差小于误差阈值,则得到太阳辐射强度预测模型。
预测时刻辐射强度计算模块140用于将预测日辐射强度预测相关时刻的太阳辐射强度作为太阳辐射强度预测模型的输入,计算预测日预测时刻的太阳辐射强度。
本实施例中即是将10月1日的8点半、9点和9点半的太阳辐射强度作为太阳辐射强度预测模型的输入,预测10月1日10点的太阳辐射强度。
由于辐射强度预测样本建立模块120建立的太阳辐射强度预测样本只包括前预设个数的比较序列,因此本实施例中在每次对太阳辐射强度进行预测时,都重新建立太阳辐射强度预测样本和太阳辐射强度预测模型,进一步提高预测准确度。
发电功率预测样本建立模块150用于提取历史数据中预测日之前发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及预测时刻的温度和太阳辐射强度,得到光伏发电功率预测样本。
提取9月每天8点半、9点和9点半的太阳辐射强度和光伏发电功率,以及10点的温度和太阳辐射强度,得到光伏发电功率预测样本。。
发电功率预测模型建立模块160用于采用BP神经网络对光伏发电功率预测样本进行训练,得到光伏发电功率预测模型。
同样也是先计算隐含层节点数。本实施例中光伏发电功率预测模型的输入为预测时刻前0.5小时、前1.0小时和前1.5小时的太阳辐射强度和光伏发电功率,以及预测时刻的温度和太阳辐射强度,输出为预测时刻的光伏发电功率。故输入层节点数为8,输出层节点数为1,取调节常数a=2,则根据可计算出隐含层节点数为5。
然后对输入层与隐含层的连接权值、隐含层与输出层的连接权值、隐含层各节点的阈值以及输出层各节点的阈值进行初始化。获取光伏发电功率预测样本中多日的发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及预测时刻的温度和太阳辐射强度,将同一日的相关数据作为一个学习向量,得到学习向量集。
将学习向量代入输入层进行训练,隐含层和输出层采用的神经元激励函数与辐射强度预测模型建立模块130中的相同。计算输出层各节点的校正误差,隐含层各节点的校正误差,对隐含层与输出层的连接权值,及输出层各节点的阈值进行修正,对输入层与隐含层的连接权值,及隐含层各节点的阈值进行修正等步骤的原理也与辐射强度预测模型建立模块130中类似,不再赘述。
预测时刻发电功率计算模块170用于将预测日发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及预测时刻的太阳辐射强度和温度作为光伏发电功率预测模型的输入,计算预测日预测时刻的光伏发电功率。
本实施例中将10月1日8点半、9点和9点半的太阳辐射强度和光伏发电功率,以及10点的太阳辐射强度和温度作为光伏发电功率预测模型的输入,预测10月1日10点的光伏发电功率。
由于发电功率预测样本建立模块150建立的光伏发电功率预测样本包括了历史数据中,预测日之前所有发电功率预测相关时刻的太阳辐射强度和光伏发电功率,数据比较多,能够保证光伏发电功率预测模型的预测准确度。本实施例中在光伏发电功率预测模型建好之后便固定不变,即是模型建好之后每次对光伏发电功率进行预测时,直接将预测日发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及预测时刻的太阳辐射强度和温度作为光伏发电功率预测模型的输入来进行预测,可节省时间。可以理解,在其他实施例中,也可以每次对光伏发电功率进行预测时,都重新建立光伏发电功率预测模型。
上述光伏发电功率预测***,采用相关性分析方法对历史数据进行分析,确定辐射强度预测相关时刻和发电功率预测相关时刻。采用BP神经网络对太阳辐射强度预测样本和光伏发电功率预测样本进行训练,得到太阳辐射强度预测模型和光伏发电功率预测模型。利用太阳辐射强度预测模型计算预测日预测时刻的太阳辐射强度,利用光伏发电功率预测模型计算预测日预测时刻的光伏发电功率。利用灰色关联度分析方法去除历史数据中关联度较低的辐射强度预测相关时刻的太阳辐射强度,提高了太阳辐射强度预测精度。采用BP神经网络良好的非线性函数逼近能力,对太阳辐射强度预测样本和光伏发电功率预测样本进行训练建立预测模型,提高了预测模型的预测精度,从而提高光伏并网***的稳定性,降低运行成本。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种光伏发电功率预测方法,其特征在于,包括以下步骤:
根据历史数据,采用相关性分析方法获取辐射强度预测相关时刻和发电功率预测相关时刻,所述历史数据包括太阳辐射强度、温度和光伏发电功率,所述辐射强度预测相关时刻指太阳辐射强度与预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻,所述发电功率预测相关时刻指太阳辐射强度和光伏发电功率与所述预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻;
提取所述历史数据中所述辐射强度预测相关时刻的太阳辐射强度,并进行灰色关联度分析,得到太阳辐射强度预测样本;
采用BP神经网络对所述太阳辐射强度预测样本进行训练,得到太阳辐射强度预测模型;
将预测日所述辐射强度预测相关时刻的太阳辐射强度作为所述太阳辐射强度预测模型的输入,计算预测日所述预测时刻的太阳辐射强度;
提取所述历史数据中预测日之前所述发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及所述预测时刻的温度和太阳辐射强度,得到光伏发电功率预测样本;
采用BP神经网络对所述光伏发电功率预测样本进行训练,得到光伏发电功率预测模型;
将预测日所述发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及所述预测时刻的太阳辐射强度和温度作为所述光伏发电功率预测模型的输入,计算预测日所述预测时刻的光伏发电功率。
2.根据权利要求1所述的光伏发电功率预测方法,其特征在于,所述根据历史数据,采用相关性分析方法获取辐射强度预测相关时刻和发电功率预测相关时刻的步骤,包括以下步骤:
提取所述历史数据中所述预测日之前,所述预测时刻及预测时刻前设定时间段内的太阳辐射强度;
根据分别计算所述设定时间段内各时刻点的太阳辐射强度,与所述预测时刻的太阳辐射强度的相关系数,其中,cov(X,Y)=E[X-E(X)]E[Y-E(Y)],X为预测时刻的太阳辐射强度,Y为设定时间段内一时刻点的太阳辐射强度,ρXY表示X和Y的相关系数,cov(X,Y)表示X和Y的协方差,E(X)和E(Y)分别表示X和Y的期望,D(X)和D(Y)分别表示X和Y的方差;
提取所述设定时间段内,太阳辐射强度与所述预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻,得到所述辐射强度预测相关时刻;
提取所述历史数据中所述预测日之前,所述预测时刻的光伏发电功率及预测时刻前设定时间段内的太阳辐射强度和光伏发电功率;
根据分别计算所述设定时间段内各时刻点的太阳辐射强度和光伏发电功率,与所述预测时刻的光伏发电功率的相关系数,其中cov(X',Y')=E[X'-E(X')]E[Y'-E(Y')],X'为预测时刻的光伏发电功率,Y'为设定时间段内一时刻点的太阳辐射强度或光伏发电功率,ρX'Y'表示X'和Y'的相关系数,cov(X',Y')表示X'和Y'的协方差,E(X')和E(Y')分别表示X'和Y'的期望,D(X')和D(Y')分别表示X'和Y'的方差;
提取设定时间范围内,太阳辐射强度和光伏发电功率与所述预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻,得到所述发电功率预测相关时刻。
3.根据权利要求1所述的光伏发电功率预测方法,其特征在于,所述根据历史数据,采用相关性分析方法获取辐射强度预测相关时刻和发电功率预测相关时刻的步骤,包括以下步骤:
提取所述历史数据中所述预测日之前,所述预测时刻及预测时刻前若干个等间隔时刻点的太阳辐射强度;
根据分别计算所述若干个等间隔时刻点的太阳辐射强度,与所述预测时刻的太阳辐射强度的相关系数,其中,cov(X,Y)=E[X-E(X)]E[Y-E(Y)],X为预测时刻的太阳辐射强度,Y为若干个等间隔时刻点中一时刻点的太阳辐射强度,ρXY表示X和Y的相关系数,cov(X,Y)表示X和Y的协方差,E(X)和E(Y)分别表示X和Y的期望,D(X)和D(Y)分别表示X和Y的方差;
提取所述若干个等间隔时刻点中,太阳辐射强度与所述预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻点,得到所述辐射强度预测相关时刻;
提取所述历史数据中所述预测日之前,所述预测时刻的光伏发电功率及预测时刻前若干个等间隔时刻点的太阳辐射强度和光伏发电功率;
根据分别计算所述若干个等间隔时刻点的太阳辐射强度和光伏发电功率,与所述预测时刻的光伏发电功率的相关系数,其中cov(X',Y')=E[X'-E(X')]E[Y'-E(Y')],X'为预测时刻的光伏发电功率,Y'为若干个等间隔时刻点中一时刻点的太阳辐射强度或光伏发电功率,ρX'Y'表示X'和Y'的相关系数,cov(X',Y')表示X'和Y'的协方差,E(X')和E(Y')分别表示X'和Y'的期望,D(X')和D(Y')分别表示X'和Y'的方差;
提取所述若干个等间隔时刻点中,太阳辐射强度和光伏发电功率与所述预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻,得到所述发电功率预测相关时刻。
4.根据权利要求1所述的光伏发电功率预测方法,其特征在于,所述提取所述历史数据中所述辐射强度预测相关时刻的太阳辐射强度,并进行灰色关联度分析,得到太阳辐射强度预测样本的步骤,包括以下步骤:
获取预测日所述辐射强度预测相关时刻的太阳辐射强度,得到参考序列;
获取预测日之前所述辐射强度预测相关时刻的太阳辐射强度,并对同一日的太阳辐射强度建立序列,得到多个比较序列;
根据对所述参考序列和比较序列进行无量纲化处理,得到无量纲化参考序列和无量纲化比较序列,其中表示序列Xi(k)中的最大值,且i=0时X0(k)为参考序列,i不为0时Xi(k)为比较序列,xi(k)表示序列Xi(k)对应的无量纲化序列;
根据 计算所述无量纲化参考序列和无量纲化比较序列的灰色关联系数,其中x0(k)为无量纲化参考序列,xi(k)为无量纲化比较序列,ρ为分辨系数,ξ0i(k)表示无量纲化参考序列x0(k)中第k个值,与无量纲化比较序列xi(k)中第k个值的灰色关联系数;
根据计算所述比较序列与参考序列的关联度,其中ξ0i(k)表示无量纲化参考序列x0(k)中第k个值,与无量纲化比较序列xi(k)中第k个值的灰色关联系数,r0i表示参考序列X0(k)与比较序列Xi(k)的关联度;
根据比较序列与参考序列关联度的大小,对比较序列进行排序,提取前预设个数的比较序列,得到所述太阳辐射强度预测样本。
5.根据权利要求1所述的光伏发电功率预测方法,其特征在于,所述采用BP神经网络对所述太阳辐射强度预测样本进行训练,得到太阳辐射强度预测模型的步骤,包括以下步骤:
初始化输入层与隐含层的连接权值、隐含层与输出层的连接权值、隐含层各节点的阈值以及输出层各节点的阈值;
获取所述太阳辐射强度预测样本中多日的辐射强度预测相关时刻的太阳辐射强度,并将同一天的所述辐射强度预测相关时刻的太阳辐射强度作为一个学习向量,得到学习向量集;
将所述学习向量集中的一个学习向量作为所述输入层的输入,计算所述隐含层的输出向量;具体为
根据所述隐含层的输出向量计算所述输出层的输出向量;具体为
f2(x)=purelin(x)=x
根据所述输出层各节点的校正误差,对隐含层与输出层的连接权值,及输出层各节点的阈值进行修正;具体为
其中α为学习速率,Δvjt表示对所述隐含层与输出层的连接权值vjt的修正值,Δγt表示对所述输出层各节点的阈值γt的修正值;
根据所述隐含层各节点的校正误差,对输入层与隐含层的连接权值,及隐含层各节点的阈值进行修正;具体为
其中β为学习速率,Δwij表示对所述输入层与隐含层的连接权值wij的修正值,Δθj表示对所述隐含层各节点的阈值θj的修正值;
判断所述学习向量集中的学习向量是否训练完毕;
若否,则返回所述将所述学习向量集中的一个学习向量作为所述输入层的输入,计算所述隐含层的输出向量的步骤;
若是,则判断全局误差是否小于误差阈值;若否,则返回所述获取所述太阳辐射强度预测样本中多日的辐射强度预测相关时刻的太阳辐射强度,并将同一天的所述辐射强度预测相关时刻的太阳辐射强度作为一个学习向量,得到学习向量集的步骤,若是,则得到所述太阳辐射强度预测模型。
6.一种光伏发电功率预测***,其特征在于,包括:
相关时刻计算模块,用于根据历史数据,采用相关性分析方法获取辐射强度预测相关时刻和发电功率预测相关时刻,所述历史数据包括太阳辐射强度、温度和光伏发电功率,所述辐射强度预测相关时刻指太阳辐射强度与预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻,所述发电功率预测相关时刻指太阳辐射强度和光伏发电功率与所述预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻;
辐射强度预测样本建立模块,用于提取所述历史数据中所述辐射强度预测相关时刻的太阳辐射强度,并进行灰色关联度分析,得到太阳辐射强度预测样本;
辐射强度预测模型建立模块,用于采用BP神经网络对所述太阳辐射强度预测样本进行训练,得到太阳辐射强度预测模型;
预测时刻辐射强度计算模块,用于将预测日所述辐射强度预测相关时刻的太阳辐射强度作为所述太阳辐射强度预测模型的输入,计算预测日所述预测时刻的太阳辐射强度;
发电功率预测样本建立模块,用于提取所述历史数据中预测日之前所述发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及所述预测时刻的温度和太阳辐射强度,得到光伏发电功率预测样本;
发电功率预测模型建立模块,用于采用BP神经网络对所述光伏发电功率预测样本进行训练,得到光伏发电功率预测模型;
预测时刻发电功率计算模块,用于将预测日所述发电功率预测相关时刻的太阳辐射强度和光伏发电功率,以及所述预测时刻的太阳辐射强度和温度作为所述光伏发电功率预测模型的输入,计算预测日所述预测时刻的光伏发电功率。
7.根据权利要求6所述的光伏发电功率预测***,其特征在于,所述相关时刻计算模块包括:
第一提取单元,用于提取所述历史数据中所述预测日之前,所述预测时刻及预测时刻前设定时间段内的太阳辐射强度;
第一计算单元,用于根据分别计算所述设定时间段内各时刻点的太阳辐射强度,与所述预测时刻的太阳辐射强度的相关系数,其中,cov(X,Y)=E[X-E(X)]E[Y-E(Y)],X为预测时刻的太阳辐射强度,Y为设定时间段内一时刻点的太阳辐射强度,ρXY表示X和Y的相关系数,cov(X,Y)表示X和Y的协方差,E(X)和E(Y)分别表示X和Y的期望,D(X)和D(Y)分别表示X和Y的方差;
第一处理单元,用于提取所述设定时间段内,太阳辐射强度与所述预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻,得到所述辐射强度预测相关时刻;
第二提取单元,用于提取所述历史数据中所述预测日之前,所述预测时刻的光伏发电功率及预测时刻前设定时间段内的太阳辐射强度和光伏发电功率;
第二计算单元,用于根据分别计算所述设定时间段内各时刻点的太阳辐射强度和光伏发电功率,与所述预测时刻的光伏发电功率的相关系数,其中cov(X',Y')=E[X'-E(X')]E[Y'-E(Y')],X'为预测时刻的光伏发电功率,Y'为设定时间段内一时刻点的太阳辐射强度或光伏发电功率,ρX'Y'表示X'和Y'的相关系数,cov(X',Y')表示X'和Y'的协方差,E(X')和E(Y')分别表示X'和Y'的期望,D(X')和D(Y')分别表示X'和Y'的方差;
第二处理单元,用于提取设定时间范围内,太阳辐射强度和光伏发电功率与所述预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻,得到所述发电功率预测相关时刻。
8.根据权利要求6所述的光伏发电功率预测***,其特征在于,所述相关时刻计算模块包括:
第一提取单元,用于提取所述历史数据中所述预测日之前,所述预测时刻及预测时刻前若干个等间隔时刻点的太阳辐射强度;
第一计算单元,用于根据分别计算所述若干个等间隔时刻点的太阳辐射强度,与所述预测时刻的太阳辐射强度的相关系数,其中,cov(X,Y)=E[X-E(X)]E[Y-E(Y)],X为预测时刻的太阳辐射强度,Y为若干个等间隔时刻点中一时刻点的太阳辐射强度,ρXY表示X和Y的相关系数,cov(X,Y)表示X和Y的协方差,E(X)和E(Y)分别表示X和Y的期望,D(X)和D(Y)分别表示X和Y的方差;
第一处理单元,用于提取所述若干个等间隔时刻点中,太阳辐射强度与所述预测时刻的太阳辐射强度相关系数高于辐射相关阈值的时刻点,得到所述辐射强度预测相关时刻;
第二提取单元,用于提取所述历史数据中所述预测日之前,所述预测时刻的光伏发电功率及预测时刻前若干个等间隔时刻点的太阳辐射强度和光伏发电功率;
第二计算单元,用于根据分别计算所述若干个等间隔时刻点的太阳辐射强度和光伏发电功率,与所述预测时刻的光伏发电功率的相关系数,其中cov(X',Y')=E[X'-E(X')]E[Y'-E(Y')],X'为预测时刻的光伏发电功率,Y'为若干个等间隔时刻点中一时刻点的太阳辐射强度或光伏发电功率,ρX'Y'表示X'和Y'的相关系数,cov(X',Y')表示X'和Y'的协方差,E(X')和E(Y')分别表示X'和Y'的期望,D(X')和D(Y')分别表示X'和Y'的方差;
第二处理单元,用于提取所述若干个等间隔时刻点中,太阳辐射强度和光伏发电功率与所述预测时刻的光伏发电功率相关系数高于功率相关阈值的时刻,得到所述发电功率预测相关时刻。
9.根据权利要求6所述的光伏发电功率预测***,其特征在于,所述辐射强度预测样本建立模块包括:
参考序列建立单元,用于获取预测日所述辐射强度预测相关时刻的太阳辐射强度,得到参考序列;
比较序列建立单元,用于获取预测日之前所述辐射强度预测相关时刻的太阳辐射强度,并对同一日的太阳辐射强度建立序列,得到多个比较序列;
无量纲化处理单元,用于根据对所述参考序列和比较序列进行无量纲化处理,得到无量纲化参考序列和无量纲化比较序列,其中表示序列Xi(k)中的最大值,且i=0时X0(k)为参考序列,i不为0时Xi(k)为比较序列,xi(k)表示序列Xi(k)对应的无量纲化序列;
灰色关联系数计算单元,用于根据 计算所述无量纲化参考序列和无量纲化比较序列的灰色关联系数,其中x0(k)为无量纲化参考序列,xi(k)为无量纲化比较序列,ρ为分辨系数,ξ0i(k)表示无量纲化参考序列x0(k)中第k个值,与无量纲化比较序列xi(k)中第k个值的灰色关联系数;
关联度计算单元,用于根据计算所述比较序列与参考序列的关联度,其中ξ0i(k)表示无量纲化参考序列x0(k)中第k个值,与无量纲化比较序列xi(k)中第k个值的灰色关联系数,r0i表示参考序列X0(k)与比较序列Xi(k)的关联度;
关联度排序单元,用于根据比较序列与参考序列关联度的大小,对比较序列进行排序,提取前预设个数的比较序列,得到所述太阳辐射强度预测样本。
10.根据权利要求6所述的光伏发电功率预测***,其特征在于,所述辐射强度预测模型建立模块包括:
隐含层节点计算单元,用于根据计算隐含层节点数,其中l为隐含层节点数,n为输入层节点数,m为输出层节点数,a为调节常数;
初始化单元,用于初始化输入层与隐含层的连接权值、隐含层与输出层的连接权值、隐含层各节点的阈值以及输出层各节点的阈值;
学习向量集建立单元,用于获取所述太阳辐射强度预测样本中多日的辐射强度预测相关时刻的太阳辐射强度,并将同一天的所述辐射强度预测相关时刻的太阳辐射强度作为一个学习向量,得到学习向量集;
隐含层输出向量计算单元,用于将所述学习向量集中的一个学习向量作为所述输入层的输入,计算所述隐含层的输出向量;具体为
输出层输出向量计算单元,用于根据所述隐含层的输出向量计算所述输出层的输出向量;具体为
f2(x)=purelin(x)=x
第一修正单元,用于根据所述输出层各节点的校正误差,对隐含层与输出层的连接权值,及输出层各节点的阈值进行修正;具体为
其中α为学习速率,Δvjt表示对所述隐含层与输出层的连接权值vjt的修正值,Δγt表示对所述输出层各节点的阈值γt的修正值;
第二修正单元,用于根据所述隐含层各节点的校正误差,对输入层与隐含层的连接权值,及隐含层各节点的阈值进行修正;具体为
其中β为学习速率,Δwij表示对所述输入层与隐含层的连接权值wij的修正值,Δθj表示对所述隐含层各节点的阈值θj的修正值;
判断单元,用于判断所述学习向量集中的学习向量是否训练完毕,并在所述学习向量集中的学习向量未训练完毕时,控制所述隐含层输出向量计算单元将所述学习向量集中的一个学习向量作为所述输入层的输入;
控制单元,用于在所述学习向量集中的学习向量训练完毕时,判断全局误差是否小于误差阈值,并在全局误差不小于误差阈值时,控制所述学习向量集建立单元获取所述太阳辐射强度预测样本中多日的辐射强度预测相关时刻的太阳辐射强度,并将同一天的所述辐射强度预测相关时刻的太阳辐射强度作为一个学习向量,得到学习向量集;若全局误差小于误差阈值,则得到所述太阳辐射强度预测模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310430694.8A CN103500365B (zh) | 2013-09-18 | 2013-09-18 | 光伏发电功率预测方法和*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310430694.8A CN103500365B (zh) | 2013-09-18 | 2013-09-18 | 光伏发电功率预测方法和*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103500365A true CN103500365A (zh) | 2014-01-08 |
CN103500365B CN103500365B (zh) | 2016-05-04 |
Family
ID=49865571
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310430694.8A Active CN103500365B (zh) | 2013-09-18 | 2013-09-18 | 光伏发电功率预测方法和*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103500365B (zh) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103984988A (zh) * | 2014-05-06 | 2014-08-13 | 国家电网公司 | 测光网络实时校正arma模型光伏功率超短期预测方法 |
CN104182564A (zh) * | 2014-06-26 | 2014-12-03 | 内蒙古电力勘测设计院有限责任公司 | 光伏发电站设计专家*** |
CN104469810A (zh) * | 2014-12-01 | 2015-03-25 | 中国联合网络通信集团有限公司 | 一种移动终端无线上网的网速计算方法及装置 |
JP5736530B1 (ja) * | 2015-02-17 | 2015-06-17 | オーナンバ株式会社 | 太陽光発電システムの未来の電流値または発電量の低下の時期を予測する方法 |
CN104978611A (zh) * | 2015-07-06 | 2015-10-14 | 东南大学 | 一种基于灰色关联分析的神经网络光伏发电出力预测方法 |
CN105005825A (zh) * | 2015-07-03 | 2015-10-28 | 深圳供电局有限公司 | 一种基于动态神经网络预测光伏功率的方法及*** |
CN105303254A (zh) * | 2015-10-26 | 2016-02-03 | 国网浙江省电力公司电力科学研究院 | 一种对光伏电站接收辐射进行预测的方法及装置 |
CN105678397A (zh) * | 2015-11-18 | 2016-06-15 | 江苏省电力公司泰州市姜堰区供电公司 | 一种基于改进EMD算法及Elman算法的短期光伏功率预测方法 |
CN106446494A (zh) * | 2016-05-11 | 2017-02-22 | 新疆大学 | 基于小波包‑神经网络的风光功率预测方法 |
CN106897998A (zh) * | 2017-02-24 | 2017-06-27 | 深圳市昊睿智控科技服务有限公司 | 太阳能直射辐射强度信息预测方法和*** |
CN107800353A (zh) * | 2017-12-01 | 2018-03-13 | 尹希月 | 一种太阳能电池输出控制方法 |
CN107886445A (zh) * | 2017-11-09 | 2018-04-06 | 王钊 | 一种基于神经元大数据分析的功率调整方法 |
CN108564192A (zh) * | 2017-12-29 | 2018-09-21 | 河海大学 | 一种基于气象因子权重相似日的短期光伏功率预测方法 |
CN109060618A (zh) * | 2018-08-20 | 2018-12-21 | 国网冀北电力有限公司张家口供电公司 | 一种光伏发电功率预测装置及方法 |
CN109390976A (zh) * | 2018-09-21 | 2019-02-26 | 武汉大学 | 一种低压台区分布式光伏发电功率辨识方法 |
CN109447843A (zh) * | 2018-10-22 | 2019-03-08 | 北京中电普华信息技术有限公司 | 一种光伏出力预测方法及装置 |
CN109829572A (zh) * | 2019-01-14 | 2019-05-31 | 国网江苏省电力有限公司苏州供电分公司 | 雷电气候下的光伏发电功率预测方法 |
CN110598896A (zh) * | 2019-07-26 | 2019-12-20 | 陕西省水利电力勘测设计研究院 | 基于预测误差校正的光伏功率预测方法 |
CN112651537A (zh) * | 2019-10-10 | 2021-04-13 | 国网河北省电力有限公司 | 一种光伏发电超短期功率预测方法及*** |
DE102015104016B4 (de) * | 2015-03-18 | 2021-06-24 | Vse Ag | Verfahren zur Erstellung einer Leistungsprognose und Prognosevorrichtung |
CN114253302A (zh) * | 2021-12-15 | 2022-03-29 | 阳光电源(上海)有限公司 | 一种跟踪控制方法及装置 |
CN116418293A (zh) * | 2023-04-21 | 2023-07-11 | 攀枝花中电光伏发电有限公司 | 一种基于物联网大数据的光伏电站智慧诊断*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102999786A (zh) * | 2012-10-17 | 2013-03-27 | 浙江埃菲生能源科技有限公司 | 基于相似日特征分类与极端学习机的光伏发电功率短期预报方法 |
JP2013099143A (ja) * | 2011-11-01 | 2013-05-20 | Nippon Telegr & Teleph Corp <Ntt> | 予測モデル構築装置、方法、及びプログラム、並びに発電量予測装置、及び方法 |
-
2013
- 2013-09-18 CN CN201310430694.8A patent/CN103500365B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013099143A (ja) * | 2011-11-01 | 2013-05-20 | Nippon Telegr & Teleph Corp <Ntt> | 予測モデル構築装置、方法、及びプログラム、並びに発電量予測装置、及び方法 |
CN102999786A (zh) * | 2012-10-17 | 2013-03-27 | 浙江埃菲生能源科技有限公司 | 基于相似日特征分类与极端学习机的光伏发电功率短期预报方法 |
Non-Patent Citations (2)
Title |
---|
王军 等: "基于模式识别和BP神经网络的光伏发电功率检测", 《第四届电能质量及柔性输电技术研讨会》, 21 November 2012 (2012-11-21) * |
陈倩: "区域风电功率一体化预测方法研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》, 15 April 2012 (2012-04-15) * |
Cited By (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103984988A (zh) * | 2014-05-06 | 2014-08-13 | 国家电网公司 | 测光网络实时校正arma模型光伏功率超短期预测方法 |
CN103984988B (zh) * | 2014-05-06 | 2018-03-02 | 国家电网公司 | 测光网络实时校正arma模型光伏功率超短期预测方法 |
CN104182564B (zh) * | 2014-06-26 | 2017-09-12 | 内蒙古电力勘测设计院有限责任公司 | 光伏发电站设计专家*** |
CN104182564A (zh) * | 2014-06-26 | 2014-12-03 | 内蒙古电力勘测设计院有限责任公司 | 光伏发电站设计专家*** |
CN104469810A (zh) * | 2014-12-01 | 2015-03-25 | 中国联合网络通信集团有限公司 | 一种移动终端无线上网的网速计算方法及装置 |
CN104469810B (zh) * | 2014-12-01 | 2018-04-17 | 中国联合网络通信集团有限公司 | 一种移动终端无线上网的网速计算方法及装置 |
JP5736530B1 (ja) * | 2015-02-17 | 2015-06-17 | オーナンバ株式会社 | 太陽光発電システムの未来の電流値または発電量の低下の時期を予測する方法 |
EP3059856A1 (en) * | 2015-02-17 | 2016-08-24 | Onamba Co., Ltd. | Method for predicting a future timing of lowering of a current value or power generation quantity of a solar power generation system |
WO2016132569A1 (ja) * | 2015-02-17 | 2016-08-25 | オーナンバ株式会社 | 太陽光発電システムの未来の電流値または発電量の低下の時期を予測する方法 |
DE102015104016B4 (de) * | 2015-03-18 | 2021-06-24 | Vse Ag | Verfahren zur Erstellung einer Leistungsprognose und Prognosevorrichtung |
CN105005825B (zh) * | 2015-07-03 | 2018-12-07 | 深圳供电局有限公司 | 一种基于动态神经网络预测光伏功率的方法及*** |
CN105005825A (zh) * | 2015-07-03 | 2015-10-28 | 深圳供电局有限公司 | 一种基于动态神经网络预测光伏功率的方法及*** |
CN104978611A (zh) * | 2015-07-06 | 2015-10-14 | 东南大学 | 一种基于灰色关联分析的神经网络光伏发电出力预测方法 |
CN105303254B (zh) * | 2015-10-26 | 2019-05-31 | 国网浙江省电力公司电力科学研究院 | 一种对光伏电站接收辐射进行预测的方法及装置 |
CN105303254A (zh) * | 2015-10-26 | 2016-02-03 | 国网浙江省电力公司电力科学研究院 | 一种对光伏电站接收辐射进行预测的方法及装置 |
CN105678397A (zh) * | 2015-11-18 | 2016-06-15 | 江苏省电力公司泰州市姜堰区供电公司 | 一种基于改进EMD算法及Elman算法的短期光伏功率预测方法 |
CN106446494A (zh) * | 2016-05-11 | 2017-02-22 | 新疆大学 | 基于小波包‑神经网络的风光功率预测方法 |
CN106446494B (zh) * | 2016-05-11 | 2019-01-11 | 新疆大学 | 基于小波包-神经网络的风光功率预测方法 |
CN106897998A (zh) * | 2017-02-24 | 2017-06-27 | 深圳市昊睿智控科技服务有限公司 | 太阳能直射辐射强度信息预测方法和*** |
CN106897998B (zh) * | 2017-02-24 | 2020-09-04 | 深圳市微埃智能科技有限公司 | 太阳能直射辐射强度信息预测方法和*** |
CN107886445A (zh) * | 2017-11-09 | 2018-04-06 | 王钊 | 一种基于神经元大数据分析的功率调整方法 |
CN107800353A (zh) * | 2017-12-01 | 2018-03-13 | 尹希月 | 一种太阳能电池输出控制方法 |
CN107800353B (zh) * | 2017-12-01 | 2019-08-16 | 肇庆学院 | 一种太阳能电池输出控制方法 |
CN108564192A (zh) * | 2017-12-29 | 2018-09-21 | 河海大学 | 一种基于气象因子权重相似日的短期光伏功率预测方法 |
CN108564192B (zh) * | 2017-12-29 | 2021-06-08 | 河海大学 | 一种基于气象因子权重相似日的短期光伏功率预测方法 |
CN109060618A (zh) * | 2018-08-20 | 2018-12-21 | 国网冀北电力有限公司张家口供电公司 | 一种光伏发电功率预测装置及方法 |
CN109390976A (zh) * | 2018-09-21 | 2019-02-26 | 武汉大学 | 一种低压台区分布式光伏发电功率辨识方法 |
CN109447843A (zh) * | 2018-10-22 | 2019-03-08 | 北京中电普华信息技术有限公司 | 一种光伏出力预测方法及装置 |
CN109447843B (zh) * | 2018-10-22 | 2022-01-28 | 北京中电普华信息技术有限公司 | 一种光伏出力预测方法及装置 |
CN109829572A (zh) * | 2019-01-14 | 2019-05-31 | 国网江苏省电力有限公司苏州供电分公司 | 雷电气候下的光伏发电功率预测方法 |
CN109829572B (zh) * | 2019-01-14 | 2021-03-12 | 国网江苏省电力有限公司苏州供电分公司 | 雷电气候下的光伏发电功率预测方法 |
CN110598896A (zh) * | 2019-07-26 | 2019-12-20 | 陕西省水利电力勘测设计研究院 | 基于预测误差校正的光伏功率预测方法 |
CN112651537A (zh) * | 2019-10-10 | 2021-04-13 | 国网河北省电力有限公司 | 一种光伏发电超短期功率预测方法及*** |
CN114253302A (zh) * | 2021-12-15 | 2022-03-29 | 阳光电源(上海)有限公司 | 一种跟踪控制方法及装置 |
CN116418293A (zh) * | 2023-04-21 | 2023-07-11 | 攀枝花中电光伏发电有限公司 | 一种基于物联网大数据的光伏电站智慧诊断*** |
CN116418293B (zh) * | 2023-04-21 | 2024-02-27 | 攀枝花中电光伏发电有限公司 | 一种基于物联网大数据的光伏电站智慧诊断*** |
Also Published As
Publication number | Publication date |
---|---|
CN103500365B (zh) | 2016-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103500365A (zh) | 光伏发电功率预测方法和*** | |
Fouilloy et al. | Solar irradiation prediction with machine learning: Forecasting models selection method depending on weather variability | |
Harrou et al. | A data-driven soft sensor to forecast energy consumption in wastewater treatment plants: A case study | |
Zhang et al. | Air quality index forecast in Beijing based on CNN-LSTM multi-model | |
CN102930358B (zh) | 一种光伏电站发电功率的神经网络预测方法 | |
Chai et al. | PV power prediction based on LSTM with adaptive hyperparameter adjustment | |
CN105426956B (zh) | 一种超短期光伏预测方法 | |
CN104408913B (zh) | 一种考虑时空相关性的交通流三参数实时预测方法 | |
KR101703981B1 (ko) | 강우예측을 이용한 하수처리시설의 제어방법 | |
CN105913151A (zh) | 基于自适应变异粒子群与bp网络光伏电站发电量预测方法 | |
CN109376906B (zh) | 基于多维度轨迹的出行时间预测方法、***及电子设备 | |
CN104598765A (zh) | 一种基于弹性自适应神经网络的建筑物能耗预测方法 | |
CN103023065A (zh) | 一种基于相对误差熵值法的风电短期功率预测方法 | |
CN111652425A (zh) | 一种基于粗糙集和长短期记忆网络的河流水质预测方法 | |
CN105117975B (zh) | 一种光伏输出功率值的分频预测方法 | |
CN102867217A (zh) | 基于投影寻踪的设施农业气象灾害风险评价方法 | |
CN103793887A (zh) | 基于自适应增强算法的短期电力负荷在线预测方法 | |
JP2017135964A (ja) | 太陽光発電システムの動作 | |
CN105678397A (zh) | 一种基于改进EMD算法及Elman算法的短期光伏功率预测方法 | |
CN105225006A (zh) | 一种短期风电功率非参数概率预测方法 | |
CN103914740B (zh) | 一种基于数据驱动的输电线路覆冰预测与自校正方法 | |
CN110991776A (zh) | 一种基于gru网络实现水位预测的方法及*** | |
CN103886397A (zh) | 一种分治法与水质周期性结合的水质预测方法 | |
CN106570594A (zh) | 一种基于tmbp的相似日光伏发电短期预测方法 | |
Mellios et al. | A multivariate analysis of the daily water demand of Skiathos Island, Greece, implementing the artificial neuro-fuzzy inference system (ANFIS) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |