CN106600038A - 一种基于马尔科夫模型的负荷区间预测方法 - Google Patents
一种基于马尔科夫模型的负荷区间预测方法 Download PDFInfo
- Publication number
- CN106600038A CN106600038A CN201611013832.2A CN201611013832A CN106600038A CN 106600038 A CN106600038 A CN 106600038A CN 201611013832 A CN201611013832 A CN 201611013832A CN 106600038 A CN106600038 A CN 106600038A
- Authority
- CN
- China
- Prior art keywords
- interval
- load
- prediction
- curve
- trend
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000012706 support-vector machine Methods 0.000 claims abstract description 7
- 230000007704 transition Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 240000002853 Nelumbo nucifera Species 0.000 claims description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 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)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于马尔科夫模型的负荷区间预测方法,首先,提取负荷数据的趋势曲线。接着用日期加趋势项负荷作为输入,趋势项负荷作为输出,在支持向量机上实现负荷的趋势项预测。从负荷曲线中减去趋势项,得到随机项。对随机项利用马尔科夫模型进行区间预测。将趋势项预测值与随机项预测区间相加,即得到负荷的预测区间值。此外,引入误差反馈环节:利用三次样条插值法提取出历史负荷曲线的上下包络线,取样本中的点进行马尔科夫区间预测,将预测区间的上、下限与上、下包络线对比,计算误差,最后将该误差加入到原始的预测区间中,即得到最终的预测区间。本发明预测结果稳定,预测精度高。
Description
技术领域
本发明涉及能源预测技术领域,尤其是指一种基于马尔科夫模型的负荷区间预测方法。
背景技术
由于电力系统中蕴含了各种不确定因素,使得决策工作必然面临一定程度的风险,所以在决策时必须考虑电力需求的不确定性。传统确定性预测方法的结果不能反映需求的不确定性,而区间预测可满足这种客观要求。区间预测的结果不是一个简单的确定性数值,而是一个区间,并且这个区间对应了一定水平的概率置信水平,能描述未来预测结果的可能范围。根据区间预测结果,电力系统决策人员在进行生产计划、系统安全分析等工作时能够更好地认识到未来负荷可能存在的不确定性和面临的风险因素,从而及时作出更为合理的决策。因此,分析电力系统负荷的变化规律,研究电力负荷区间预测方法,实现电力负荷的不确定性预测具有重要的理论意义和实用价值。
目前区间预测方法存在以下不足:
1、计算复杂;
2、假设性强;
3、计算时间长。
发明内容
本发明的目的在于克服现有技术的缺点和不足,提供一种基于马尔科夫模型的负荷区间预测方法,实现对负荷的中长期准确预报。
为实现上述目的,本发明所提供的技术方案为:一种基于马尔科夫模型的负荷区间预测方法,首先,提取负荷数据的趋势曲线,然后用日期加趋势项负荷作为输入,趋势项负荷作为输出,在支持向量机上实现负荷的趋势项预测;接着从负荷曲线中减去趋势项,得到随机项,对随机项利用马尔科夫模型进行区间预测;接着将趋势项预测值与随机项预测区间相加,即得到负荷的预测区间值;此外,为了提高预测的精度,引入误差反馈环节:首先,利用三次样条插值法提取出历史负荷曲线的上、下包络线,取样本中的点进行马尔科夫区间预测,将预测区间的上、下限与上、下包络线对比,计算误差,最后将该误差加入到原始的预测区间中,即得到最终的预测区间;其包括以下步骤:
1)找出历史负荷曲线的趋势曲线,历史曲线减去趋势项,得到随机项;
2)对预测日的趋势曲线部分进行预测,预测模型为支持向量机,模型输入为日期和历史负荷的趋势项,输出为负荷的趋势项;
3)将随机项曲线按幅值最大最小为上下限平均划分成3个区间,计算转移概率矩阵,转移概率最大的区间,即为随机项的预测区间;
4)负荷的趋势项预测值与随机项的预测区间之和为原始预测区间;
5)通过三次样条插值的方法找出负荷曲线的上、下包络线;
6)计算上、下包络线与原始预测区间上、下限之间的误差;
7)原始预测区间与误差之和为最终的预测区间。
本发明与现有技术相比,具有如下优点与有益效果:
1、预测结果稳定,预测精度高;
2、本发明预测方法是数据驱动、自适应的方法,其预测结果不依赖于使用者的先验知识;
3、本发明思路简单、直观、易操作;
4、本发明对电网的调度、规划等实际情况具有很强的指导意义。
附图说明
图1为本发明基于马尔科夫模型的负荷区间预测方法的流程图之一。
图2为本发明基于马尔科夫模型的负荷区间预测方法的流程图之二。
图3为理想预报结果与实际数据的曲线图。
具体实施方式
下面结合具体实施例对本发明做进一步的说明。
如图1和图2所示,本实施例所述的基于马尔科夫模型的负荷区间预测方法,具体是:首先,提取负荷数据的趋势曲线,然后用日期加趋势项负荷作为输入,趋势项负荷作为输出,在支持向量机上实现负荷的趋势项预测;接着从负荷曲线中减去趋势项,得到随机项,对随机项利用马尔科夫模型进行区间预测;接着将趋势项预测值与随机项预测区间相加,即得到负荷的预测区间值;此外,为了提高预测的精度,引入误差反馈环节:首先,利用三次样条插值法提取出历史负荷曲线的上、下包络线,取样本中的点进行马尔科夫区间预测,将预测区间的上、下限与上、下包络线对比,计算误差,最后将该误差加入到原始的预测区间中,即得到最终的预测区间;其包括以下步骤:
1)找出历史负荷曲线的趋势曲线,历史曲线减去趋势项,得到随机项;
假设历史曲线f(t),t=1,2,…,N,其趋势曲线为F(t),则
S(t)=f(t)-F(t)
其中,f(t)为历史负荷曲线,F(t)为趋势曲线,S(t)为随机项。
2)对预测日的趋势曲线部分进行预测,预测模型为支持向量机,模型输入为日期和历史负荷的趋势项,输出为负荷的趋势项;
假设负荷曲线上第t个点代表Y年M月D日的负荷,则支持向量机模型的输入可表示为:
Xtra=[Y M D F(t-3) F(t-2) F(t-1)]
其中,Y,M,D分别代表与负荷F(t)对应的年、月、日。F(t-3),F(t-2),F(t-1)为第t-3,t-2,t-1天的负荷趋势项。
输出可表示为
Ytra=F(t)
其中,F(t)为第t天的负荷趋势项。
3)将随机项曲线按幅值最大最小为上下限平均划分成3个区间,计算转移概率矩阵,转移概率最大的区间,即为随机项的预测区间;
其中,Ω1,Ω2,Ω3代表区间1,区间2,区间3。
转移概率矩阵可表示为
其中,Aij(i=1,2,3;j=1,2,3)表示上一时刻的数值S(t-1)位于i区间,下一时刻的数值S(t)位于j区间的概率。即:
其中,Bij(t)为判断是否满足上一时刻的数值S(t-1)位于i区间,下一时刻的数值S(t)位于j区间的变量,若满足,则为1,否则为0。
若预测前一天的负荷值位于第i个区间,则预测日的负荷随机项区间SΩ(t)为Ai1,Ai2,Ai3中的最大值所在的区间。
4)负荷的趋势项预测值与随机项的预测区间之和为原始预测区间;
其中,为趋势项预测值,SΩ(t)预测日的负荷随机项区间,为原始预测区间。
5)通过三次样条插值的方法找出负荷曲线的上、下包络线fh和fl;
6)计算上、下包络线与原始预测区间上、下限之间的误差;
首先对预测日之前10天的负荷执行步骤(1)至步骤(4),得到这10天的区间预测值
其中,errorh为上限误差,fh为上包络线的值,为预测区间的上限值;errorl为下限误差,fl为下包络线的值,为预测区间的下限值。
7)原始预测区间与误差之和为最终的预测区间。
假设原始日最高负荷曲线为x,预测区间上限构成的上曲线为f,预测区间下限构成的下曲线为g。在该预测模型中,从准确度和精确度两方面评估本预测方法的预测精度,提出五个评估指标:
1、准确度
若预测点的实际数值落在预测结果的上限和下限之间,为预测准确。
其中,x为原始日最高负荷曲线,f为预测区间上限构成的上曲线,g为预测区间下限构成的下曲线。
2、精确度
2.1、上调平均距离与上调平均比例
上调平均距离:
上调平均比例:
其中,x为原始日最高负荷曲线,f为预测区间上限构成的上曲线,上划线表示取平均值。d1为上调平均距离,D1为上调平均比例。
2.2、下调平均距离与下调平均比例
下调平均距离:
下调平均比例:
其中,x为原始日最高负荷曲线,g为预测区间下限构成的下曲线。上划线表示取平均值。d2为下调平均距离,D2为下调平均比例。
上(下)调平均距离体现预测区间上(下)限与真实值之间的整体差距,为具体数值;上(下)调平均比例体现该整体差距在真实值中所占的比例,为百分值,可实现不同数据间的指标对比。在这四个指标中,值越小,表示预测值与真实值之间的阈值越小,体现了预测的精确度。
由以上分析可以看出,准确度与精确度之间存在相互制约的关系,准确度高的模型必将导致精确度偏低,精确度高的模型则准确度将受制约。如何在两者中进行权衡选择,则需要对具体情况进行分析。
如图3所示,为理想预报结果与实际数据的曲线图。由图可看出,在绝大多数情况下,真实值位于预测区间上限与预测区间下限之间,准确度极高。预测区间上限与预测区间下限较好地贴合了真实曲线的波动情况,误差较小,精确度较高。
以上所述实施例只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。
Claims (1)
1.一种基于马尔科夫模型的负荷区间预测方法,其特征在于:首先,提取负荷数据的趋势曲线,然后用日期加趋势项负荷作为输入,趋势项负荷作为输出,在支持向量机上实现负荷的趋势项预测;接着从负荷曲线中减去趋势项,得到随机项,对随机项利用马尔科夫模型进行区间预测;接着将趋势项预测值与随机项预测区间相加,即得到负荷的预测区间值;此外,为了提高预测的精度,引入误差反馈环节:首先,利用三次样条插值法提取出历史负荷曲线的上、下包络线,取样本中的点进行马尔科夫区间预测,将预测区间的上、下限与上、下包络线对比,计算误差,最后将该误差加入到原始的预测区间中,即得到最终的预测区间;其包括以下步骤:
1)找出历史负荷曲线的趋势曲线,历史曲线减去趋势项,得到随机项;
2)对预测日的趋势曲线部分进行预测,预测模型为支持向量机,模型输入为日期和历史负荷的趋势项,输出为负荷的趋势项;
3)将随机项曲线按幅值最大最小为上下限平均划分成3个区间,计算转移概率矩阵,转移概率最大的区间,即为随机项的预测区间;
4)负荷的趋势项预测值与随机项的预测区间之和为原始预测区间;
5)通过三次样条插值的方法找出负荷曲线的上、下包络线;
6)计算上、下包络线与原始预测区间上、下限之间的误差;
7)原始预测区间与误差之和为最终的预测区间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611013832.2A CN106600038A (zh) | 2016-11-18 | 2016-11-18 | 一种基于马尔科夫模型的负荷区间预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611013832.2A CN106600038A (zh) | 2016-11-18 | 2016-11-18 | 一种基于马尔科夫模型的负荷区间预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106600038A true CN106600038A (zh) | 2017-04-26 |
Family
ID=58591489
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611013832.2A Pending CN106600038A (zh) | 2016-11-18 | 2016-11-18 | 一种基于马尔科夫模型的负荷区间预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106600038A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108764528A (zh) * | 2018-04-23 | 2018-11-06 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种基于影响因子分析的日最大负荷区间预测方法 |
CN112651545A (zh) * | 2020-12-07 | 2021-04-13 | 华电郑州机械设计研究院有限公司 | 一种基于特性分解的工业热负荷预测方法 |
CN115310212A (zh) * | 2022-10-12 | 2022-11-08 | 中汽研(天津)汽车工程研究院有限公司 | 一种汽车减震器特性数据抽样方法 |
-
2016
- 2016-11-18 CN CN201611013832.2A patent/CN106600038A/zh active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108764528A (zh) * | 2018-04-23 | 2018-11-06 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种基于影响因子分析的日最大负荷区间预测方法 |
CN112651545A (zh) * | 2020-12-07 | 2021-04-13 | 华电郑州机械设计研究院有限公司 | 一种基于特性分解的工业热负荷预测方法 |
CN115310212A (zh) * | 2022-10-12 | 2022-11-08 | 中汽研(天津)汽车工程研究院有限公司 | 一种汽车减震器特性数据抽样方法 |
CN115310212B (zh) * | 2022-10-12 | 2023-03-24 | 中汽研(天津)汽车工程研究院有限公司 | 一种汽车减震器特性数据抽样方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Machine learning for solar irradiance forecasting of photovoltaic system | |
Hu et al. | A single-loop Kriging surrogate model method by considering the first failure instant for time-dependent reliability analysis and safety lifetime analysis | |
CN110135635A (zh) | 一种区域电力饱和负荷预测方法及系统 | |
CN103218675A (zh) | 一种基于聚类和滑动窗口的短期负荷预测方法 | |
CN103488869A (zh) | 一种最小二乘支持向量机的风力发电短期负荷预测方法 | |
CN109800898A (zh) | 一种智能短期负荷预测方法及系统 | |
CN107239849A (zh) | 一种基于改进的指数平滑灰色模型电力负荷预测方法 | |
CN110737876B (zh) | 区域电网光伏功率预测优化方法和装置 | |
CN102509173A (zh) | 一种基于马尔科夫链的电力系统负荷准确预测方法 | |
CN106600038A (zh) | 一种基于马尔科夫模型的负荷区间预测方法 | |
CN110705771A (zh) | 区域电网新能源发电功率预测优化方法和装置 | |
CN102831064A (zh) | 一种面向可靠性评估的软件自适应测试方法 | |
Yang et al. | Sub-minute probabilistic solar forecasting for real-time stochastic simulations | |
CN108764528A (zh) | 一种基于影响因子分析的日最大负荷区间预测方法 | |
CN115081307A (zh) | 一种基于神经网络的气象数值模式同化方法及装置 | |
Han et al. | Economic dispatch considering the wind power forecast error | |
CN105676645A (zh) | 一种基于函数型权rbf-arx模型的双回路水箱液位预测控制方法 | |
CN107358059A (zh) | 短期光伏能量预测方法和装置 | |
KR102570248B1 (ko) | 태양광 발전량 예측 시스템 및 방법 | |
Hayajneh et al. | Intelligent Solar Forecasts: Modern Machine Learning Models & TinyML Role for Improved Solar Energy Yield Predictions | |
CN111476437A (zh) | 一种光伏发电功率超短期区间预测的方法及设备 | |
CN110543942A (zh) | 一种多时空长短记忆深度网络精准预测方法 | |
Dong et al. | Very short-term photovoltaic power forecasting using uncertain basis function | |
Ling et al. | A thermal load forecasting algorithm based on trajectory tracking | |
Silva et al. | Application of a hybrid neural fuzzy inference system to forecast solar intensity |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170426 |
|
RJ01 | Rejection of invention patent application after publication |