CN109948281B - 基于天气类型有效识别和组合预测的直散分离建模方法 - Google Patents
基于天气类型有效识别和组合预测的直散分离建模方法 Download PDFInfo
- Publication number
- CN109948281B CN109948281B CN201910249101.5A CN201910249101A CN109948281B CN 109948281 B CN109948281 B CN 109948281B CN 201910249101 A CN201910249101 A CN 201910249101A CN 109948281 B CN109948281 B CN 109948281B
- Authority
- CN
- China
- Prior art keywords
- weather
- model
- weather type
- prediction
- optimal
- 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
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于天气类型有效识别和组合预测的直散分离建模方法,Sandia方法考虑了气象要素的综合变化特征,挑选出长序列气象数据中具有气候特征代表性的典型年资料,使直散分离模型的研究和分析结果更具有典型性和说服力;依据修正后的清晰度指数进行天气类型的划分;每种天气类型下的最优经典小时模型选定后进行本地化修正,得到适用于所有天气类型的本地组合预测模型;通过偏相关分析剔除一部分对散射比影响较小的因子,再通过主成分分析提取主成分建立线性模型,有效提取原有数据携带的隐含信息;从每种天气类型下的线性模型和本地化模型中选择该种天气类型下的最优模型,再进行组合得到组合预测模型进行天气预测。提高预测结果准确度。
Description
技术领域
本发明涉及一种光伏技术,特别涉及一种基于天气类型有效识别和组合预测的直散分离建模方法。
背景技术
自20世纪70年代以来,全世界大部分地区能源消耗成倍增加,气候变化明显,随之而来的是一系列环境问题,太阳能等清洁能源的开发利用越来越受到关注。据世界能源委员会的研究报告称,到21世纪下半叶,太阳能将成为能源利用中非常重要的一种。光伏发电是目前太阳能利用中技术较成熟的一种太阳能发电系统,具有不消耗燃料、不排放污染物、规模灵活、安全可靠、维护简单等优点。太阳能光伏发电将是今后太阳能利用的主要利用形式,有离网和并网两种形式,而并网光伏发电是目前乃至将来的主流趋势。太阳能光伏发电技术的应用在欧美国家早已完成初期开发和示范阶段,目前正在大批量生产和规模应用。
在当下光伏装机快速增长的大环境里,以提高电站收益率为出发点的设计方案已引起全行业的关注。在系统优化设计及功率预报中,产生了进行直散分离获得水平面直接辐射和散射辐射的新需求,以此为固定式、跟踪式阵列斜面辐射计算、光伏功率预测以及聚光热发电设计等提供依据,同时对建筑能耗动态模拟研究具有一定参考价值。
目前,国内外研究人员对于直散分离预测模型的研究多以单一预测模型为主,然而单一预测模型本身具有局限性,预测精度提高很有限,且各单一预测模型的准确度有所不同,预测存在不确定性。
常见的直散分离预测模型大多只考虑少量的气象因素如清晰度指数和日照百分率等,这将影响模型预测精度;然而如果考虑大量的气象因素,不仅会增加模型复杂度,预测效果也会受各气象因素之间的多重共线性关系的影响。所以在建立预测模型时设法降低气象因素间的多重共线性关系对于模型预测精度的提高显得尤其重要。
由于太阳能资源评估缺少详细的技术规范,目前基于长期辐射序列获取典型气象年辐射资料的方法除上述的气候平均法之外,频率(数)最大法也在光伏电站资源评估中有了初步的应用。这两种算法均具有明确的数学统计意义,但在典型天气形势的刻画方面不具有代表性。因此,太阳能资源评估领域中需要进一步提出切实可行的、能较好地代表评估区域长期气候特征的典型年辐射数据计算方法,以有效提高太阳能资源评估结果的科学代表性。
发明内容
本发明是针对使用单一预测模型进行直散分离建模,来对太阳辐射进行预测存在的问题,提出了一种基于天气类型有效识别和组合预测的直散分离建模方法,首先采用Sandia法,对北京地区多年小时尺度原始数据选取典型气象年;其次根据修正后的清晰度指数划分天气类型,在每种天气类型下,对比选择最优经典小时预测模型进行本地化修正,并根据每种天气类型所占比例得到适用于所有天气类型的组合预测模型;同时在每种天气类型下对与散射比相关的PM2.5、总云量、能见度等影响因子进行偏相关分析,利用主成分分析法对偏相关强的因子提取主成分从而建立线性模型;然后将本地化的最优经典小时模型和线性模型进行比较,得到每种天气类型下的最优模型;最后根据每种天气类型所占总样本的比例,将各个天气类型下对应的最优模型进行组合,得到适用于所有天气类型的最终组合预测模型。
本发明的技术方案为:一种基于天气类型有效识别和组合预测的直散分离建模方法,具体包括如下步骤:
1)收集整理数据,并对数据进行样本筛选和数据质量检查,根据Sandia法对北京地区多年气象、辐射数据选取典型气象月组成典型气象年;
2)为了降低太阳高度角对清晰度指数的影响,对天气类型中清晰度指数进行修正,依据修正后的清晰度指数k′T进行天气类型的划分;根据日出、日落时角,对北京地区总样本日照时间内的小时辐射数据进行筛选,计算散射比,并将其作为散射比的实测值,选取三种典型直散分离模型计算水平面上散射辐照度与水平面上太阳总辐照度比值作为预测值,选用平均绝对百分比误差、相对均方根误差及相关系数作为评估模型预测效果的指标,分别计算并选取每种天气类型下误差较小、相关系数大的模型为对应天气类型下的最优经典小时模型;
其中kT是修正前的清晰度指数;m是大气质量;
3)对步骤2)得到的每种天气类型下的最优经典小时模型进行本地化修正,即再利用大样本数据中以频率代概率的方法,根据本地每类天气所占本地总样本的比例得到该种天气类型下最优模型的权重,最终得到适用于所有天气类型的本地组合预测模型;
4)在每种天气类型下,对与散射比相关的各种气象环境影响因子分别进行偏相关分析,并根据主成分分析对偏相关强的因子进行主成分选取,利用选取的主成分建立线性模型;
5)在每种天气类型下,将线性模型和本地化修正后的最优经典小时模型做比较选择对应最优模型;
6)对步骤5)选择的最优模型,根据不同天气类型在总样本中的占比情况,确定每种天气类型模型的权重系数,得到该地区的最终预测模型。
本发明的有益效果在于:本发明基于天气类型有效识别和组合预测的直散分离建模方法,直散分离组合模型使单个模型预测值过于偏大或偏小的极端情况得以削弱,降低了单一预测模型出现极端误差的概率,有效地解决了单一模型预测散射比不稳定的问题;Sandia方法考虑了气象要素的综合变化特征,可以挑选出长序列气象数据中具有气候特征代表性的典型年资料,从而代表典型大气环境的演变特点,使直散分离模型的研究和分析结果更具有典型性和说服力;直散分离预测模型若只考虑少量气象因素,将会影响模型预测精度;如果考虑大量气象因素,不仅会增加模型复杂度,预测效果也会受各气象因素之间的多重共线性关系的影响,因而通过偏相关分析剔除一部分对散射比影响较小的因子,再通过主成分分析提取主成分进行建模,可有效提取原有数据携带的隐含信息并使模型更简洁、方便使用;此建模方法给单一模型预测结果一定权值,可以优化直散分离预测模型最终的输出结果,降低均方根误差和平均相对误差,从而提高预测结果准确度。
附图说明
图1为本发明基于天气类型有效识别和组合预测的直散分离建模方法流程示意图。
具体实施方式
Sandia方法选取典型气象年:
Sandia方法可挑选出长序列气象数据中具有气候特征代表性的典型月份和典型年资料,首先采用Finkelstein-Schafer统计算法,挑选出12个均具有代表性的典型月组成一个典型气象年。而典型月的选取需要考虑各个气象要素在大气环境中所占的权重,被分析的气象要素应尽可能代表大气环境的整体变化特征,且相互独立。
根据以下技术路线实现典型气象年资料的挑选:
1、计算各个气象要素分布的长期累积分布函数值和逐年月累积分布函数值。
式中:Sn(x)为要素x处的长期累积分布值;n为样本总数;k为要素x在增序时间序列中的排序。然后,按照年份顺序分别在目标月份内按照增序排列,再按上式计算该要素逐年的月累积分布值。
2、对每个气象要素的分布计算其Finkelstein-Schafer统计值Cfs,即:
式中:δi为各要素长期累积分布值与逐年各月累积分布值的绝对差值;nd为各分析月内的天数。在获得各气象要素每个月份的Cfs后,按一定权重系数WFi把各Cfs汇总成一个参数Ws,该系数代表了各气象要素在大气环境中的作用比重。
式中:KK为气象要素的个数。
天气类型分类:
清晰度指数表示大气的透明程度,与天气状况及太阳辐射密切相关,其公式为:
式中,I为水平面上太阳总辐照度,其值是到达地表水平面上的太阳直射辐照度和散射辐照度的总和;
大气层外水平面上太阳辐射量I0:
ESC=1367±7W/m2
式中Tn为一年中的日期序号,如元旦为Tn=1,春分日为Tn=81,12月31日为Tn=365。
然而清晰度指数不仅与气象条件相关,而且与天空中太阳位置有关。为了降低太阳高度角对清晰度指数的影响,对其进行修正如下:
其中kT是修正前的清晰度指数,k′T是修正后的清晰度指数,m是大气质量。
按照修正后的清晰度指数k′T将不同时次对应的天气状态进行划分归类。类型1:当k′T≥0.5时,定义为晴、晴转多云、多云转晴;类型2:当0.5>k′T≥0.2时,定义为多云、阴转多云、多云转阴;类型3:当0.2>k′T≥0.1时,定义为小雨、阵雨、小雪、轻雾、霾等;类型4:当k′T<0.1,定义为中雨及以上、中雪及以上,从而深度探讨不同天气类型下逐时散射比的预测效果及与各气象因子的关系。
不同模型计算公式:
Erbs模型:
Orgill and Hollands模型:
Reindl(1990)模型:
式中,Id为水平面上散射辐照度。
偏相关分析:
多个变量间的关系往往是比较复杂的,简单相关分析分析的可能只是表面的非本质的相互关系,没有把其他变量的影响考虑在内。偏相关系数就是在多元回归分析中,在消除其他变量影响的条件下,所计算的某2个变量之间的相关系数,它所反映的变量之间的相关性会更准确。考虑到光伏影响因子的实际情况,本专利取r≥0.3的因子为强相关因子。
主成分分析:
主成分分析是最常用的一种降维方法,其算法描述如下:
1、对偏相关分析中得到的对散射比影响较强的气象因子的标准化的矩阵X*求解其相关系数矩阵R,由相关系数矩阵R计算其特征值λ1≥λ2≥…≥λN及对应的特征向量α1,α2,…,αN。
2、计算贡献率ek和累计贡献率Et。
通过累计方差贡献率确定主成分个数,一般选取累计方差贡献率位于75%~95%范围内的特征值λ1,λ2,…,λp所对应的第1,第2,…,第p个主成分。
3、P个主成分构成的主成分矩阵为:
U=X*T
其中
T=[α1,α2,^αp]
式中,T为p个主成分对应特征向量。
如图1所示基于天气类型有效识别和组合预测的直散分离建模方法流程示意图,方法步骤如下:
S1:收集整理数据,并对数据进行样本筛选和数据质量检查,根据Sandia法对北京地区多年气象、辐射数据选取典型气象月组成典型气象年;
S2:依据修正的清晰度指数k′T进行天气类型的划分,根据日出、日落时角,对北京地区总样本日照时间内的小时辐射数据进行筛选,计算散射比,并将其作为散射比的实测值,上述三种典型直散分离模型计算结果(Id/I)作为预测值。选用平均绝对百分比误差(MAPE)、相对均方根误差(NRMSE)及相关系数(CORR)作为评估模型预测效果的指标,分别计算并选取每种天气类型下误差较小、相关系数较大的模型为对应天气类型下的最优经典小时模型;
S3:分别本地化修正每种天气类型下的最优经典小时模型,即采用1stOpt软件中Levenberg-Marquardt法(LM)+通用全局优化算法(UGO)对最优经典小时模型进行分段函数拟合,然后根据每类天气所占总样本的比例得到该种天气类型下最优模型的权重,最终得到适用于所有天气类型的组合预测模型;
S4:在每种天气类型下,对与散射比相关的PM2.5、总云量、能见度等因子进行偏相关分析,并根据主成分分析对偏相关强的因子(考虑到光伏影响因子的实际情况,本专利取r≥0.3的因子为强相关因子)进行主成分选取,利用选取的主成分建立线性模型;
S5:在每种天气类型下,将线性模型和本地化修正后的最优经典小时模型做比较选择最优模型;
S6:根据不同天气类型在总样本中的占比情况,确定每种天气类型模型的权重系数,得到该地区的最终预测模型。
Claims (1)
1.一种基于天气类型有效识别和组合预测的直散分离建模方法,其特征在于,具体包括如下步骤:
1)收集整理数据,并对数据进行样本筛选和数据质量检查,根据Sandia法对北京地区多年气象、辐射数据选取典型气象月组成典型气象年;
2)为了降低太阳高度角对清晰度指数的影响,对天气类型中清晰度指数进行修正,依据修正后的清晰度指数kT′进行天气类型的划分;根据日出、日落时角,对北京地区总样本日照时间内的小时辐射数据进行筛选,计算散射比,并将其作为散射比的实测值,选取三种典型直散分离模型计算水平面上散射辐照度与水平面上太阳总辐照度比值作为预测值,选用平均绝对百分比误差、相对均方根误差及相关系数作为评估模型预测效果的指标,分别计算并选取每种天气类型下误差小、相关系数大的模型为对应天气类型下的最优经典小时模型;
其中kT是修正前的清晰度指数;m是大气质量;
式中,I为水平面上太阳总辐照度,其值是到达地表水平面上的太阳直射辐照度和散射辐照度的总和;
大气层外水平面上太阳辐射量I0:
ESC=1367±7W/m2
式中Tn为一年中的日期序号,元旦为Tn=1,春分日为Tn=81,12月31日为Tn=365;所述三种典型直散分离模型的水平面上散射辐照度Id与水平面上太阳总辐照度I比值如下:
Erbs模型:
Orgill and Hollands模型:
Reindl(1990)模型:
3)对步骤2)得到的每种天气类型下的最优经典小时模型进行本地化修正,即采用1stOpt软件中Levenberg-Marquardt法+通用全局优化算法对最优经典小时模型进行分段函数拟合;再根据本地每类天气所占本地总样本的比例得到该种天气类型下最优模型的权重,最终得到适用于所有天气类型的本地组合预测模型;4)在每种天气类型下,对与散射比相关的各种气象环境影响因子分别进行偏相关分析,并根据主成分分析对偏相关强的因子进行主成分选取,利用选取的主成分建立线性模型;
5)在每种天气类型下,将线性模型和本地化修正后的最优经典小时模型做比较选择对应最优模型;
6)对步骤5)选择的最优模型,根据不同天气类型在总样本中的占比情况,确定每种天气类型模型的权重系数,得到该地区的最终预测模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910249101.5A CN109948281B (zh) | 2019-03-29 | 2019-03-29 | 基于天气类型有效识别和组合预测的直散分离建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910249101.5A CN109948281B (zh) | 2019-03-29 | 2019-03-29 | 基于天气类型有效识别和组合预测的直散分离建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109948281A CN109948281A (zh) | 2019-06-28 |
CN109948281B true CN109948281B (zh) | 2023-05-05 |
Family
ID=67013011
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910249101.5A Active CN109948281B (zh) | 2019-03-29 | 2019-03-29 | 基于天气类型有效识别和组合预测的直散分离建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109948281B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111539846B (zh) * | 2020-04-22 | 2022-10-25 | 上海电力大学 | 基于天气类型细分的光伏功率预测方法 |
CN111596384B (zh) * | 2020-05-11 | 2022-05-17 | 国网浙江省电力有限公司电力科学研究院 | 一种基于天气类型有效识别的斜面辐射预测方法 |
CN111815020A (zh) * | 2020-06-04 | 2020-10-23 | 上海电力大学 | 一种基于太阳辐射气候特征识别的南墙辐射预测方法 |
CN111815021B (zh) * | 2020-06-04 | 2024-05-31 | 上海电力大学 | 一种基于太阳辐射气候特征识别的光伏功率预测方法 |
CN111967655A (zh) * | 2020-07-28 | 2020-11-20 | 中国南方电网有限责任公司 | 一种短期负荷预测方法及系统 |
CN113379143A (zh) * | 2021-06-23 | 2021-09-10 | 阳光电源股份有限公司 | 一种典型气象年的构建方法、发电量预测方法及相关装置 |
CN113722375B (zh) * | 2021-08-10 | 2023-10-17 | 国网浙江省电力有限公司电力科学研究院 | 一种基于多算法优化的双层混合光伏预测方法及计算机设备 |
KR102526361B1 (ko) * | 2022-11-09 | 2023-04-26 | 한국외국어대학교 연구산학협력단 | 지역단위의 폭염 및 한파에 대한 지표면효과 분석 시스템 및 방법 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104021427A (zh) * | 2014-06-10 | 2014-09-03 | 上海电力学院 | 一种基于因子分析的并网光伏电站日发电量预测方法 |
CN105184678A (zh) * | 2015-09-18 | 2015-12-23 | 齐齐哈尔大学 | 基于多种神经网络组合算法的光伏电站发电量短期预测模型的构建方法 |
CN107991721A (zh) * | 2017-11-21 | 2018-05-04 | 上海电力学院 | 一种基于天文和气象环境因子的逐时散射比预测方法 |
CN109299820A (zh) * | 2018-09-20 | 2019-02-01 | 国网河南省电力公司电力科学研究院 | 基于斜面辐射修正的光伏发电功率预测方法及装置 |
-
2019
- 2019-03-29 CN CN201910249101.5A patent/CN109948281B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104021427A (zh) * | 2014-06-10 | 2014-09-03 | 上海电力学院 | 一种基于因子分析的并网光伏电站日发电量预测方法 |
CN105184678A (zh) * | 2015-09-18 | 2015-12-23 | 齐齐哈尔大学 | 基于多种神经网络组合算法的光伏电站发电量短期预测模型的构建方法 |
CN107991721A (zh) * | 2017-11-21 | 2018-05-04 | 上海电力学院 | 一种基于天文和气象环境因子的逐时散射比预测方法 |
CN109299820A (zh) * | 2018-09-20 | 2019-02-01 | 国网河南省电力公司电力科学研究院 | 基于斜面辐射修正的光伏发电功率预测方法及装置 |
Non-Patent Citations (4)
Title |
---|
R.Perez.Making Full Use of The Clearness Index for Parameterizing Hourly Insolation Conditions.《Solar Energy》.1990,第45卷(第45期),111-114. * |
常蕊.太阳能资源典型年挑选方法的适用性对比研究.《高原气象》.2017,第36卷(第36期),1713-1720. * |
李芬.基于GRA-BPNN时变权重的光伏短期出力组合预测.《可再生能源》.2018,第36卷(第36期),1605-1611. * |
李芬.基于PCA-LMBP神经网络的北京地区直散分离预测.《水电能源科学》.2017,第35卷(第35期),208-212. * |
Also Published As
Publication number | Publication date |
---|---|
CN109948281A (zh) | 2019-06-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109948281B (zh) | 基于天气类型有效识别和组合预测的直散分离建模方法 | |
CN107991721B (zh) | 一种基于天文和气象环境因子的逐时散射比预测方法 | |
Zeng et al. | Daily global solar radiation in China estimated from high‐density meteorological observations: a random forest model framework | |
Yohanna et al. | A model for determining the global solar radiation for Makurdi, Nigeria | |
Zhao et al. | Optimal PV panel tilt angle based on solar radiation prediction | |
CN109726867B (zh) | 一种基于多模式集合的高分辨率电力气象预报方法 | |
Besson et al. | Long-term soiling analysis for three photovoltaic technologies in Santiago region | |
CN111815021B (zh) | 一种基于太阳辐射气候特征识别的光伏功率预测方法 | |
Durre et al. | Factors influencing the cold-season diurnal temperature range in the United States | |
Ogbulezie et al. | A review of regression models employed for predicting diffuse solar radiation in North-Western Africa | |
CN113837426A (zh) | 一种基于天气分型的光伏功率预测方法 | |
Salazar et al. | Practical models to estimate horizontal irradiance in clear sky conditions: Preliminary results | |
Pattnayak et al. | How good is regional climate model version 4 in simulating the monsoon onset over Kerala? | |
Czekalski et al. | Parameterization of daily solar irradiance variability | |
CN106355264B (zh) | 一种参考作物腾发量组合预测方法 | |
CN111815020A (zh) | 一种基于太阳辐射气候特征识别的南墙辐射预测方法 | |
Rincón et al. | Assessment of short-term irradiance forecasting based on post-processing tools applied on WRF meteorological simulations | |
Lin et al. | Assessing forecasting performance of daily reference evapotranspiration: A comparative analysis of updated temperature penman-monteith and penman-monteith forecast models | |
Sawadogo et al. | Hourly global horizontal irradiance over West Africa: A case study of one-year satellite-and reanalysis-derived estimates vs. in situ measurements | |
CN111596384B (zh) | 一种基于天气类型有效识别的斜面辐射预测方法 | |
CN111950813B (zh) | 一种气象干旱监测与预测方法 | |
Park et al. | An improved photovoltaic system output prediction model under limited weather information | |
Pereira et al. | Prediction of solar resource and photovoltaic energy production through the generation of a typical meteorological year and Meso-NH simulations: Application to the south of Portugal | |
Ghasemi Mobtaker et al. | Estimation of monthly mean daily global solar radiation in Tabriz using empirical models and artificial neural networks | |
Wilson et al. | Evaluation of solar capacity factor of~ 2000 solar plants across the united states using multilayer perceptron regressor models |
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 |