CN107527116B - 基于支持向量回归的短期负荷预测方法 - Google Patents
基于支持向量回归的短期负荷预测方法 Download PDFInfo
- Publication number
- CN107527116B CN107527116B CN201710749300.3A CN201710749300A CN107527116B CN 107527116 B CN107527116 B CN 107527116B CN 201710749300 A CN201710749300 A CN 201710749300A CN 107527116 B CN107527116 B CN 107527116B
- Authority
- CN
- China
- Prior art keywords
- day
- load
- max
- daily
- similarity
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000013507 mapping Methods 0.000 claims abstract description 9
- 230000002354 daily effect Effects 0.000 claims description 37
- 230000008859 change Effects 0.000 claims description 7
- 238000012417 linear regression Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000012706 support-vector machine Methods 0.000 claims description 5
- 230000009471 action Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000003203 everyday effect Effects 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 16
- 238000010801 machine learning Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000000611 regression analysis Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- 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
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Primary Health Care (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (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)
- Feedback Control In General (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种基于支持向量回归的短期负荷预测方法,步骤如下:一、建立增量最大负荷模型和最小负荷模型:选取日特征量,包括类型日、最高温度、平均温度、平均相对湿度、平均风速;二、利用日特征量相似度,选取与预测日的日特征相似度较高的日期,并计算相似日的权重系数,最后获得预测日各点负荷;本发明针对夏冬两季气象因素修正了传统的休息日归一化数据,选取前四日最大(最小)负荷、最高温度、平均温度、平均风速、平均相对湿度作为最大(最小)负荷预测回归模型的输入,使用LS‑SVM建立了输入输出的映射关系,并取得了较好的预测效果。
Description
技术领域
本发明涉及的是一种基于支持向量回归的短期负荷预测方法。
背景技术
短期负荷预测是能量管理系统(EMS)的一个重要模块,在电力系统经济调度中具有重要的意义。日前经济调度根据第二天的负荷预测曲线,分配各台机组的出力计划。因此,提高短期负荷预测的准确度,一方面能够提高供电企业和发电厂的经济效益,另一方面够减轻日内计划的平衡压力,有助于系统电力系统安全稳定运行。
传统的负荷预测方法有时间序列法、回归分析法等。时间序列法对影响负荷的因素考虑较少,如气象因素突然变化时预测的准确度会降低。回归分析法考虑了气象因素、历史负荷对预测日负荷的影响,由于负荷与其影响因素之间强的非线性,该方法预测精度通常不高。随着机器学习算法研究的深入,神经网络和支持向量机等技术逐步引入到电力系统的短期负荷预测模型中。机器学习算法能够对训练数据进行学习,获得负荷的预测模型,然而机器学习算法容易陷入过拟合。支持向量机的结构风险最小化原理,使得其可以较好的防止过拟合。
发明内容
本发明的目的在于针对现有技术的缺陷和不足,提供了一种基于支持向量回归的短期负荷预测方法。
为实现上述目的,本发明采用的技术方案是:
一种基于支持向量回归的短期负荷预测方法,步骤如下:
一、增量最大负荷支持向量回归模型:
首先选取与气象数据不相关月份的日平均负荷(工作日与休息日分别计算)作为基准,然后预测日负荷减去日平均负荷为日增量负荷,可认为日增量负荷最大值与气象因素具有强相关性。
选取日特征量,包括类型日、最高温度、平均温度、平均相对湿度、平均风速,电力系统的日最大负荷曲线受该日及前几日的日特征量影响,建立日增量最大负荷回归模型如下:
xk=[ymax(k-1),ymax(k-2),ymax(k-3),ymax(k-4),u(k),u(k-1),u(k-2),u(k-3)]
(2)
则式(1)变为
采用最小二乘支持向量机(LS-SVM)来辨识回归模型,LS-SVM是把输入数据xk映射到一个高维空间,在高维空间构造一个线性回归函数,可以近似表示为:
其中训练数据集N为训练数据集个数,xk为式(2)表示的第k日输入向量,ymax,k为第k日增量最大负荷,另外w为权向量,表示某种非线性映射,可以将输入空间映射到高维特征空间,b为偏置,LS-SVM将回归转换成优化问题:
其中ek为拟合误差,γ为规则化参数,起到调节对拟合误差的惩罚程度的作用;采用Lagrange乘数法,求解一系列偏微分方程,可以构建下列线性方程组:
其中α=[α1 α2 L αN]T称为支持向量,αk=γek,γ为式(5)的规则化参数,1N=[1,1,L,1]T,ymax=[ymax,1,ymax,2,Lymax,N]T。IN为单位矩阵,核矩阵如下:
LS-SVM方法采用等式约束代替传统SVM的不等式约束,将回归简化为可以求解的一系列线性方程式,计算出式(6)中的b和α,回归模型变为:
K为满足Mercer条件的核函数,xk为训练数据集,x为新的输入数据,选用常用的高斯核函数
其中σ为宽度参数,控制函数的径向作用范围,||g||表示欧氏距离。这样增量最大负荷回归模型就变为:
二、负荷变化系数:
采用日特征相似度对日特征量进行衡量,日特征相似度是描述两日的日特征量的相似程度,其定义为:假设每日考虑了H个负荷相关因素,i,j两日的日特征量向量分别为:(ui1,ui2 L uiH)T,(uj1,uj2 L ujN)T,i,j两日的日特征相似度为:
根据公式(11)计算,选择与预测日的日特征相似度较高的日期,求取该日各点负荷归一化值,计算方法为:
其中k为选取的日特征相似度较高的日期个数,L(k,i)为第k日,第i点负荷数据,Lmax,k为第k日最大负荷数据,Lmin,k为第k日最小负荷数据;
按照相似度的大小,计算预测日的日负荷变化系数:
其中μk为第k个相似类型日所占的权重系数,
Okj为第k个相似类型日与预测日的相似度;
预测日的各点负荷模型为:
采用上述结构后,本发明的有益效果为:本发明针对夏冬两季气象因素修正了传统的休息日归一化数据,选取前四日最大(最小)负荷、最高温度、平均温度、平均风速、平均相对湿度作为最大(最小)负荷预测回归模型的输入,使用LS-SVM建立了输入输出的映射关系,能够提高对负荷的预测能力。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合具体实施方式,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施方式仅用以解释本发明,并不用于限定本发明。
实施例:
一种基于支持向量回归的短期负荷预测方法,步骤如下:
一、增量最大负荷支持向量回归模型:
首先选取日特征量,包括类型日、最高温度、平均温度、平均相对湿度、平均风速,电力系统的日最大负荷曲线受该日及前几日的日特征量影响,建立日增量最大负荷回归模型如下:
xk=[ymax(k-1),ymax(k-2),ymax(k-3),ymax(k-4),u(k),u(k-1),u(k-2),u(k-3)]
(2)
则式(1)变为
采用最小二乘支持向量机(LS-SVM)来辨识回归模型,LS-SVM是把输入数据xk映射到一个高维空间,在高维空间构造一个线性回归函数,可以近似表示为:
其中训练数据集N为训练数据集个数,xk为式(2)表示的第k日输入向量,ymax,k为第k日增量最大负荷,另外w为权向量,表示某种非线性映射,可以将输入空间映射到高维特征空间,b为偏置,LS-SVM将回归转换成优化问题:
其中ek为拟合误差,γ为规则化参数,起到调节对拟合误差的惩罚程度的作用;采用Lagrange乘数法,求解一系列偏微分方程,可以构建下列线性方程组:
其中α=[α1 α2 L αN]T称为支持向量,αk=γek,γ为式(5)的规则化参数,1N=[1,1,L,1]T,ymax=[ymax,1,ymax,2,L ymax,N]T。IN为单位矩阵,核矩阵如下:
LS-SVM方法采用等式约束代替传统SVM的不等式约束,将回归简化为可以求解的一系列线性方程式,计算出式(6)中的b和α,回归模型变为:
K为满足Mercer条件的核函数,xk为训练数据集,x为新的输入数据,选用常用的高斯核函数
其中σ为宽度参数,控制函数的径向作用范围,||g||表示欧氏距离。这样增量最大负荷回归模型就变为:
二、负荷变化系数:
采用日特征相似度对日特征量进行衡量,日特征相似度是描述两日的日特征量的相似程度,其定义为:假设每日考虑了H个负荷相关因素,i,j两日的日特征量向量分别为:(ui1,ui2 L uiH)T,(uj1,uj2 L ujN)T,i,j两日的日特征相似度为:
根据公式(11)计算,选择与预测日的日特征相似度较高的日期,求取该日各点负荷归一化值,计算方法为:
其中k为选取的日特征相似度较高的日期个数,L(k,i)为第k日,第i点负荷数据,Lmax,k为第k日最大负荷数据,Lmin,k为第k日最小负荷数据;
按照相似度的大小,计算预测日的日负荷变化系数:
其中μk为第k个相似类型日所占的权重系数,
Okj为第k个相似类型日与预测日的相似度;
预测日的各点负荷模型为:
本具体实施方式针对夏冬两季气象因素修正了传统的休息日归一化数据,选取前四日最大(最小)负荷、最高温度、平均温度、平均风速、平均相对湿度作为最大(最小)负荷预测回归模型的输入,使用LS-SVM建立了输入输出的映射关系,并取得了较好的预测效果,考虑温度的积分作用和重大社会活动,能够提高对负荷的预测能力。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (1)
1.一种基于支持向量回归的短期负荷预测方法,其特征在于:步骤如下:
一、增量最大负荷支持向量回归模型:
首先选取日特征量,包括类型日、最高温度、平均温度、平均相对湿度和平均风速,电力系统的日最大负荷曲线受某日及前几日的日特征量影响,建立日增量最大负荷回归模型如下:
xk=[ymax(k-1),ymax(k-2),ymax(k-3),ymax(k-4),u(k),u(k-1),u(k-2),u(k-3)] (2)
则式(1)变为
采用最小二乘支持向量机(LS-SVM)来辨识回归模型,LS-SVM是把输入数据xk映射到一个高维空间,在高维空间构造一个线性回归函数,可以表示为:
其中训练数据集N为训练数据集个数,xk为式(2)表示的第k日输入向量,ymax,k为第k日增量最大负荷,另外w为权向量,表示某种非线性映射,可以将输入空间映射到高维特征空间,b为偏置,LS-SVM将回归转换成优化问题:
其中ek为拟合误差,γ为规则化参数,起到调节对拟合误差的惩罚程度的作用;采用Lagrange乘数法,求解一系列偏微分方程,可以构建下列线性方程组:
其中α=[α1 α2 L αN]T称为支持向量,αk=γek,γ为式(5)的规则化参数,1N=[1,1,L,1]T,ymax=[ymax,1,ymax,2,L ymax,N]T,IN为单位矩阵,核矩阵如下:
LS-SVM方法采用等式约束代替传统SVM的不等式约束,将回归简化为可以求解的一系列线性方程式,计算出式(6)中的b和α,回归模型变为:
K为满足Mercer条件的核函数,xk为训练数据集,x为新的输入数据,选用常用的高斯核函数
其中σ为宽度参数,控制函数的径向作用范围,||g||表示欧氏距离,这样增量最大负荷回归模型就变为:
二、负荷变化系数:
采用日特征相似度对日特征量进行衡量,日特征相似度是描述两日的日特征量的相似程度,其定义为:假设每日考虑了H个负荷相关因素,i,j两日的日特征量向量分别为:(ui1,ui2L uiH)T,(uj1,uj2L ujN)T,i,j两日的日特征相似度为:
根据公式(11)计算,选择与预测日的日特征相似度较高的日期,求取该日各点负荷归一化值,计算方法为:
其中k为选取的日特征相似度较高的日期个数,L(k,i)为第k日、第i点负荷数据,Lmax,k为第k日最大负荷数据,Lmin,k为第k日最小负荷数据;
按照相似度的大小,计算预测日的日负荷变化系数:
其中μk为第k个相似类型日所占的权重系数,
Okj为第k个相似类型日与预测日的相似度;
预测日的各点负荷模型为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710749300.3A CN107527116B (zh) | 2017-08-28 | 2017-08-28 | 基于支持向量回归的短期负荷预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710749300.3A CN107527116B (zh) | 2017-08-28 | 2017-08-28 | 基于支持向量回归的短期负荷预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107527116A CN107527116A (zh) | 2017-12-29 |
CN107527116B true CN107527116B (zh) | 2021-02-02 |
Family
ID=60682741
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710749300.3A Expired - Fee Related CN107527116B (zh) | 2017-08-28 | 2017-08-28 | 基于支持向量回归的短期负荷预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107527116B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108416366B (zh) * | 2018-02-06 | 2021-09-21 | 武汉大学 | 一种基于气象指数的加权ls-svm的电力系统短期负荷预测方法 |
CN108879656B (zh) * | 2018-06-09 | 2021-08-10 | 宝鸡文理学院 | 一种基于子抽样svr集成的短期电力负荷预测方法 |
CN109002937A (zh) * | 2018-09-07 | 2018-12-14 | 深圳供电局有限公司 | 电网负荷预测方法、装置、计算机设备和存储介质 |
CN109460917A (zh) * | 2018-11-08 | 2019-03-12 | 中国南方电网有限责任公司 | 一种基于分布因子和支持向量机的母线负荷预测方法 |
CN109711621A (zh) * | 2018-12-27 | 2019-05-03 | 浙江华云信息科技有限公司 | 基于判别分析和支持向量机结合的工业园区负荷预测方法 |
CN111311025B (zh) * | 2020-03-17 | 2023-08-08 | 南京工程学院 | 一种基于气象相似日的负荷预测方法 |
CN111860979B (zh) * | 2020-07-01 | 2022-08-12 | 广西大学 | 一种基于tcn与ipso-lssvm组合模型的短期负荷预测方法 |
CN112330077B (zh) * | 2021-01-04 | 2021-09-24 | 南方电网数字电网研究院有限公司 | 电力负荷预测方法、装置、计算机设备和存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102270279A (zh) * | 2011-07-27 | 2011-12-07 | 华北电力大学 | 一种短期电力负荷预测方法 |
CN103488869A (zh) * | 2013-08-23 | 2014-01-01 | 上海交通大学 | 一种最小二乘支持向量机的风力发电短期负荷预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9020874B2 (en) * | 2011-10-31 | 2015-04-28 | Siemens Aktiengesellschaft | Short-term load forecast using support vector regression and feature learning |
-
2017
- 2017-08-28 CN CN201710749300.3A patent/CN107527116B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102270279A (zh) * | 2011-07-27 | 2011-12-07 | 华北电力大学 | 一种短期电力负荷预测方法 |
CN103488869A (zh) * | 2013-08-23 | 2014-01-01 | 上海交通大学 | 一种最小二乘支持向量机的风力发电短期负荷预测方法 |
Non-Patent Citations (1)
Title |
---|
基于改进最小二乘支持向量机的短期负荷预测;孙薇等;《电力科学与工程》;20151231;第31卷(第12期);16-21 * |
Also Published As
Publication number | Publication date |
---|---|
CN107527116A (zh) | 2017-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107527116B (zh) | 基于支持向量回归的短期负荷预测方法 | |
Gao et al. | Day-ahead power forecasting in a large-scale photovoltaic plant based on weather classification using LSTM | |
CN104636822B (zh) | 一种基于elman神经网络的居民负荷预测方法 | |
CN111932402B (zh) | 一种基于相似日和lstm的短期电力负荷双向组合预测方法 | |
CN111079980B (zh) | 基于自适应分类策略和混合优化svr的光功率预测方法 | |
CN108063456B (zh) | 分布式光伏发电并网规划方法及终端设备 | |
CN111460001B (zh) | 一种配电网理论线损率评估方法及系统 | |
CN108596242A (zh) | 基于小波神经网络和支持向量机的电网气象负荷预测方法 | |
CN111626473A (zh) | 一种考虑误差修正的两阶段光伏功率预测方法 | |
CN113705929B (zh) | 一种基于负荷特征曲线及典型特征值融合的春节假期负荷预测方法 | |
CN109299814B (zh) | 一种气象影响电量分解预测方法 | |
CN109034500A (zh) | 一种多时序协同的中期电力负荷预测方法 | |
CN112926772A (zh) | 一种基于lstm-gpr混合模型的光能预测方法 | |
CN109460893B (zh) | 一种光伏电站天气类型相关性指标计算方法和系统 | |
CN117332880A (zh) | 基于Shapley值的负荷预测模型组合方法 | |
CN112598181A (zh) | 负荷预测方法、装置、设备和存储介质 | |
CN110276478B (zh) | 基于分段蚁群算法优化svm的短期风电功率预测方法 | |
CN109409604B (zh) | 一种基于遗传算法-支持向量机的冷负荷预测方法 | |
CN116760006A (zh) | 一种风电功率预测方法、装置、设备及存储介质 | |
CN116542385A (zh) | 一种针对风、光功率预测的集成学习方法及系统 | |
Chen et al. | Research on medium-long term power load forecasting method based on load decomposition and big data technology | |
CN112348235A (zh) | 风光母线负荷自适应预测方法、装置、计算机设备 | |
CN114925931A (zh) | 台区负荷预测方法及系统 | |
CN114676815A (zh) | 一种太阳辐射预测方法、系统、存储介质及电子设备 | |
He et al. | Research on Photovoltaic Power Forecasting Based on SOM Weather Clustering |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210202 Termination date: 20210828 |