CN107506868A - 一种短时电力负荷预测的方法及装置 - Google Patents
一种短时电力负荷预测的方法及装置 Download PDFInfo
- Publication number
- CN107506868A CN107506868A CN201710770943.6A CN201710770943A CN107506868A CN 107506868 A CN107506868 A CN 107506868A CN 201710770943 A CN201710770943 A CN 201710770943A CN 107506868 A CN107506868 A CN 107506868A
- Authority
- CN
- China
- Prior art keywords
- data
- forecast model
- load
- electric power
- forecast
- 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
- 230000005611 electricity Effects 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000012549 training Methods 0.000 claims description 55
- 230000006870 function Effects 0.000 claims description 30
- 238000004458 analytical method Methods 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 7
- 239000012141 concentrate Substances 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 3
- 238000009826 distribution Methods 0.000 abstract description 18
- 238000005457 optimization Methods 0.000 abstract description 8
- 230000002159 abnormal effect Effects 0.000 abstract description 4
- 239000002245 particle Substances 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000013277 forecasting method Methods 0.000 description 9
- 230000008859 change Effects 0.000 description 6
- 241001269238 Data Species 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 238000010606 normalization Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 210000002569 neuron Anatomy 0.000 description 3
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000002354 daily effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- BULVZWIRKLYCBC-UHFFFAOYSA-N phorate Chemical compound CCOP(=S)(OCC)SCSCC BULVZWIRKLYCBC-UHFFFAOYSA-N 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000013432 robust analysis Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 238000005303 weighing 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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/048—Activation functions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- 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)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Software Systems (AREA)
- Marketing (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Public Health (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)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种短时电力负荷预测的方法,通过综合分位回归和鲁棒极限学习机,并利用混合粒子群算法(PSOGSA)优化后所建立的混合预测模型来对电力负荷进行预测,分位回归利用历史电力数据影响因素的多个分位数来得到未来某时刻电力负荷预测数据的条件分布的相应的分位数方程,分位回归中输入的电力数据的随机扰动不需要做任何分布上的假定,就可以详细的描述预测负荷值的统计分布,使得整个预测模型具有很强的稳健性;而鲁棒极限学习机对异常负荷值的鲁棒性更强,将上述两种方法结合在一起,并通过PSOGSA优化后所形成的混合模型可以准确的对电力负荷进行预测;本发明还提供了一种短时电力负荷预测的装置,同样具有上述有益效果。
Description
技术领域
本发明涉及数据分析领域,特别是涉及一种短时电力负荷预测的方法及装置。
背景技术
由于电不能直接储存的特点,当今社会在进行供电时都是随时发电随时使用的状态。但是由于在日常生活中,人们使用电量总体是波动的,是存在波峰波谷的,但是人们的用电状况可以进行预测,以此使得电力公司可以准确的进行供电。
短期电力负荷预测对电力系统的运行和规划具有重要的经济意义,准确的电力负荷预测有助于电力公司做出合理的发电计划,可以有效减少电力公司的运营成本。
在现有技术中,通常是使用人工神经网络结合最小二乘法来对电力负荷进行预测。在具体操作中需要根据电力公司的历史电力负荷数据,预先建立相应的预测模型,之后再根据输入的实时电力负荷数据对后续一段时间内的电力负荷数据进行预测。
但是在现有技术中,普遍存在预测精度不高的情况,这会增加电力公司的运营成本。
发明内容
本发明的目的是提供一种短时电力负荷预测的方法,可以准确对电力系统中的电力负荷进行预测;本发明的另一目的在于提供一种短时电力负荷预测的装置,可以准确对电力系统中的电力负荷进行预测。
为解决上述技术问题,本发明提供一种短时电力负荷预测的方法,所述方法包括:
获取历史电力数据;
调用预先建立的预测模型,根据所述历史电力数据计算出电力负荷预测数据,所述预测模型为综合分位回归和鲁棒极限学习机,并利用PSOGSA优化后所建立的预测模型;所述预测模型为通过预先获取的训练数据集训练得到的预测模型。
可选的,所述的历史电力数据为历史电力负荷数据。
可选的,所述预测模型具体为:
其中Yt为所述预测模型的输出值,Qt为t时刻所述预测模型在τ分位下的输出值,为所述分位回归的影响参数,为所述鲁棒极限学习机的影响参数;相应的,训练所述预测模型时所用的损失函数为:
其中C是惩罚系数,β是所述鲁棒极限学习机中隐含层与输出层之间的权值,yt是训练所述预测模型时所使用的训练数据中t时刻的电力负荷值,ρτ为检验函数。
可选的,所述训练数据集具体为:
将预先获取的历史电力数据进行归一化处理,以生成归一化的样本数据集;
通过自相关分析,提取所述归一化的样本数据集中各变量间的自相关系数;
根据所述自相关系数和预先设定的选取规则,从所述归一化的样本数据集中选取可用的样本变量数据,以构成的所述训练数据集。
本发明还提供了一种短时电力负荷预测的装置,所述装置包括:
获取模块:用于获取历史电力数据;
调用模块:用于调用预先建立的预测模型,根据所述历史电力数据计算出电力负荷预测数据,所述预测模型为综合分位回归模型和鲁棒极限学习机,并利用PSOGSA优化后所建立的预测模型;所述预测模型为通过预先获取的训练数据集训练得到的预测模型。
可选的,所述获取模块具体用于:
获取历史电力数据,其中历史电力数据为历史电力负荷数据。
可选的,所述调用模块具体用于:
调用预先建立的预测模型,根据所述历史电力负荷数据计算出电力负荷预测数据,所述预测模型具体为:
其中Yt为所述预测模型的输出值,Qt为t时刻所述预测模型在τ分位下的输出值,为所述分位回归的影响参数,为所述鲁棒极限学习机的影响参数;相应的,训练所述预测模型时所用的损失函数为:
其中C是惩罚系数,β是所述鲁棒极限学习机中隐含层与输出层之间的权值,yt是训练所述预测模型时所使用的训练数据中t时刻的电力负荷值,ρτ为检验函数。
可选的,所述调用模块具体用于:
调用预先建立的预测模型,根据所述历史电力数据计算出电力负荷预测数据,所述预测模型为通过预先获取的训练数据集训练得到的预测模型,所述训练数据集具体为将预先获取的历史电力数据进行归一化处理,以生成归一化的样本数据集;通过自相关分析,提取所述归一化的样本数据集中各变量间的自相关系数;根据所述自相关系数和预先设定的选取规则,从所述归一化的样本数据集中选取可用的样本变量数据,以构成的所述训练数据集。
本发明所提供的一种短时电力负荷预测的方法,可以通过综合分位回归和鲁棒极限学习机,并利用PSOGSA优化后所建立的混合预测模型来对电力负荷进行预测,分位回归利用历史电力数据影响因素的多个分位数来得到未来某时刻电力负荷预测数据的条件分布的相应的分位数方程,相较传统的最小二乘法而言,分位回归中输入的电力数据的随机扰动不需要做任何分布上的假定,就可以详细的描述预测负荷值的统计分布,使得整个预测模型具有很强的稳健性;而鲁棒极限学习机在继承极限学习机高泛化性能的同时,该模型对异常负荷值的鲁棒性更强,将上述两种方法结合在一起形成的混合模型可以准确的对电力负荷进行预测;本发明还提供了一种短时电力负荷预测的装置,同样具有上述有益效果,在此不再进行赘述。
附图说明
为了更清楚的说明本发明实施例或现有技术的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例所提供的一种短时电力负荷预测方法的流程图;
图2为本发明实施例所提供的一种短时电力负荷预测方法中预测模型的训练流程图;
图3输入的异常值为10个时预测模型输出的数据分布图;
图4输入的异常值为20个时预测模型输出的数据分布图;
图5输入的异常值为50个时预测模型输出的数据分布图;
图6输入的异常值为80个时预测模型输出的数据分布图;
图7为自相关系数分布图;
图8为MAE指标随候选数据个数变化的趋势图;
图9为MAPE指标随候选数据个数变化的趋势图;
图10为RMSE指标随候选数据个数变化的趋势图;
图11为本发明实施例提供的一种短时电力负荷预测装置的结构框图。
具体实施方式
本发明的核心是提供一种短时电力负荷预测的方法。在现有技术中,通常是使用人工神经网络结合最小二乘法来对电力负荷进行预测。在具体操作中需要根据电力公司的历史电力数据,预先建立相应的预测模型,之后再根据输入的实时电力负荷数据对后续一段时间内的电力负荷数据进行预测。作为估计回归系数最基础的方法,最小二乘法描述了自变量X对因变量Y的均值影响,但是最小二乘法对于预测模型中输入的电力数据的随机扰动需要提前做出分布假定,例如预先假定所述随机扰动为均值分布、正态分布等等。但是在实际生活当中,上述假设通常是不能满足的,数据经常会出现尖峰,厚尾以及存在显著的异方差等情况,此时使用最小二乘法预测出来的电力负荷数据往往存在较大误差。
而本发明所提供的一种短时电力负荷预测的方法,可以通过综合分位回归和鲁棒极限学习机,并利用PSOGSA优化后所建立的混合预测模型来对电力负荷进行预测,分位回归利用历史电力数据影响因素的多个分位数来得到未来某时刻电力负荷预测数据的条件分布的相应的分位数方程,相较传统的最小二乘法而言,分位回归中输入的电力数据的随机扰动不需要做任何分布上的假定,就可以详细的描述预测负荷值的统计分布,使得整个预测模型具有很强的稳健性;而鲁棒极限学习机在继承极限学习机高泛化性能的同时,该模型对异常负荷值的鲁棒性更强,将上述两种方法结合在一起形成的混合模型可以准确的对电力负荷进行预测。
为了使本技术领域的人员更好地理解本发明方案,下面结合附图和具体实施方式对本发明作进一步的详细说明。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参考图1,图1为本发明实施例所提供的第一种短时电力负荷预测方法的流程图,该方法包括:
S101:获取历史电力数据。
在本步骤中,会先获取历史电力数据,其中有关具体获取数据的过程可以是操作人员手动输入,也可以是自动将当前的电力数据导入下述的预测模型,当然还可以选择其他获取数据的方法,在此不做具体限定。
其中,所述电力数据为历史电力负荷数据,可能也会考虑实时电价数据,详细内容将在下述实施例中做详细说明。
S102:调用预先建立的预测模型,根据历史电力数据计算出电力负荷预测数据。
在本发明实施例中,所述预测模型为综合分位回归和鲁棒极限学习机,并利用PSOGSA优化后所建立的预测模型;所述预测模型为通过预先获取的训练数据集训练得到的预测模型,所述训练数据集为根据历史电力数据得到的所述训练数据集。
其中,鲁棒极限学习机是一种前馈神经网络的算法,基于鲁棒极限学习机所建立的预测模型通常具有输入层,隐藏层和输出层,其中隐藏层通常具有多个。在确定所述预测模型中的具体参数时,需要先对所述预测模型进行训练,而训练所用的数据通常是预先存储的用于训练预测模型的历史电力负荷数据。
对于给定的N个训练数据集其中,对于第i个训练数据集,xi=[xi1,xi2,…,xin]为输入向量,yi为相应的期望输出值。此时,对于包含个隐含层,激励函数为g(x)的鲁棒极限学习机的数学模型可以表示为:
其中,激励函数g(x)的形式可以有多种,例如sigmoid、cos、threshold等等,在此不做具体限定,本发明实施例中所使用的激励函数为sigmoid函数:g(x)=1/(1+e-x)。
上述鲁棒极限学习机的数学模型中,wi=[wi1,wi1,…,w1n]T为连接隐含层中第i个神经元和输入层之间的权值;为上述第i个神经元的偏差;βi=[βi1,βi2,…,βin]T为连接所述第i个神经元和输出层之间的权值。
上述鲁棒极限学习机的数学模型可以用矩阵的形式来进行表达:
Hβ=y
其中:
上式中,H为隐含层的输出矩阵。
为了方便理解,在本实施例接下来的叙述中,将在一维情况下进行详细说明。上述输出权值的求解是要保证损失函数取得最小值,所述损失函数如下:
其中,γ为调节系数,用于权衡训练误差和输出权值。
而在本发明中,是将分位回归与上述鲁棒极限学习机结合在一起,所提供的预测模型为:
其中,Yt为所述预测模型的输出值,Qt为t时刻所述预测模型在τ分位下的输出值,υ是线性分位回归模型的参数向量,ut=xt=(lt,pt)T是上述线性分位回归模型的影响因子,为所述分位回归的影响参数;h(xt)是上述H矩阵在只考虑一个样本集合时候的情况,即在一维时所能出现的情况,为所述鲁棒极限学习机的影响参数。
相应的,训练所述预测模型时所用的损失函数为:
其中C是惩罚系数,β是所述鲁棒极限学习机中隐含层与输出层之间的权值,yt是训练所述预测模型时所使用的训练数据中t时刻的电力负荷值,ρτ为检验函数。在本发明实施例中,是通过最小化上述损失函数来得出最优的参数,其中检验函数ρτ定义如下:
其中qτ(xt)是在一维情况时,上述预测模型所输出的电力负荷预测数据。
在确定出上述损失函数之后,构建拉格朗日函数去解决上述损失函数的优化问题,即:
其中,λ∈Rn是拉格朗日向量乘子,μ=2N/||y||1,μ是惩罚系数。
利用上式对参数{βτ,υτ}求偏导为零的函数的最优解,具体过程如下:
其中,
在求得上述参数{βτ,υτ}之后,就得到了经过训练的预测模型,在输入S101中的电力数据来计算得出本发明实施例中的电力负荷预测数据。
有关混合粒子群算法(PSOGSA)的详细内容,将在下述事实例中做详细描述。
本发明所提供的一种短时电力负荷预测的方法,可以通过综合分位回归和鲁棒极限学习机,并利用PSOGSA优化后所建立的混合预测模型来对电力负荷进行预测,分位回归利用历史电力数据影响因素的多个分位数来得到未来某时刻电力负荷预测数据的条件分布的相应的分位数方程,相较传统的最小二乘法而言,分位回归中输入的电力数据的随机扰动不需要做任何分布上的假定,就可以详细的描述预测负荷值的统计分布,使得整个预测模型具有很强的稳健性;而鲁棒极限学习机在继承极限学习机高泛化性能的同时,该模型对异常负荷值的鲁棒性更强,将上述两种方法结合在一起形成的混合模型可以准确的对电力负荷进行预测。
在对所述预测模型进行训练时,可以进一步的提高预测模型的训练速度以及提高所述预测模型的精度,具体的优化步骤请参照下述实施例。
请参考图2,图2为本发明实施例所提供的一种短时电力负荷预测方法中预测模型的训练流程图,该方法包括:
S201:获取历史电力数据。
本步骤与上述实施例中S101大体相同,详细情况请参照上述实施例,在此不再进行赘述。
在本发明实施例中,所述历史电力数据会包括有历史电力负荷数据,详细内容将在下述实施例中做详细说明。
S202:历史电力数据进行归一化处理,以生成归一化的样本数据集。
在本步骤中,会先对历史电力数据进行归一化,具体步骤如下:
对于所述历史电力数据,假定其中第i个样本集包括j个历史电力负荷数据,用L表示历史电力负荷数据,则得到的样本集为:
xi={l1,l2,l3...lj};
对上述样本集进行归一化处理,需要使其中的数据分布在[0,1],可以通过下式进行归一化:
其中,表示归一化后的历史电力负荷数据,lmax为上述样本集xi={l1,l2,l3...lj}中的最大值。此时可以得到所述归一化数据
在对数据进行归一化之后,使得输入的数据之间的差异得到了控制,同时还保持了数据独立的类别信息。在对数据归一化之后,可以提高上述发明实施例中的预测模型的计算速度和预测精度。
S203:通过自相关分析,提取归一化的样本数据集中各变量间的自相关系数。
在本步骤中,由于电力负荷数据是跟时间有关的数据,所以电力负荷数据可以按照时间的先后循序排列成时间序列样本集合xt,t是给定的时刻,是所述样本集合xt中数据的平均值。根据实际情况预先设置延迟时刻k,由于电力负荷数据呈现周期性,通常是以一周为单位进行循环,假如输入的电力负荷数据是每30分钟获取一次电力负荷数据,那么一周就能获取336个电力负荷数据,此时延迟时刻k可选为335。
在t和t-k时刻的电力负荷数据相关性定义如下:
通过上式计算得到的rk为自相关系数,rk的值在[-1,1]之间变换,若其值越接近1,则说明上述电力负荷数据之间的正相关性越强;若其值越接近-1,则说明上述电力负荷数据之间的负相关性越强;若越接近0,则代表着缺乏相关性。
S204:根据自相关系数和预先设定的选取规则,从归一化的样本数据集中选取可用的样本变量数据,以构成的训练数据集。
在本步骤中,会先根据上述延迟时刻的先后顺序对S203中计算得到的自相关系数进行排序,此时会得到包含多个波峰与波谷的自相关系数分布图,之后会根据所述自相关系数分布图选取相应的电力负荷数据,例如选取上述自相关系数分布图中最大波峰值所对应的延迟时刻的电力负荷数据以及该延迟时刻前面10个延迟时刻所对应的电力负荷数据作为可用的样本数据同样选取第二大和第三大的波峰值所对应的延迟时刻的电力负荷数据以及该延迟时刻前后各3个延迟时刻所对应的电力负荷数据作为可用的样本数据对于后续波峰,选取其延迟时刻以及前后各一个延迟时刻所对应的电力负荷数据作为可用的样本数据
将上述选取的可用的样本数据作为训练数据集VD={V1,V2.....,Vb}。当然,除了上述选取相关性较大的电力负荷数据的方法之外,还可以通过其他的方法选取电力负荷数据,例如在选取每个波峰的峰值所对应的延迟时刻以及该延迟时刻前后各5个延迟时刻所对应的电力负荷数据,还可以通过其他的方式选取电力负荷数据,在此不做具体限定。
S205:根据训练数据集对预测模型进行训练。
根据所述训练数据集对预测模型进行训练的相关内容以在上述实施例中所详细说明,在此不再进行赘述。
S206:在训练过程中,利用PSOGSA算法对所述预测模型中的参数进行优化。
在本发明实施例中,对于预测模型中的参数,可以进一步的通过PSOGSA算法对预测模型中的参数进行优化。所述PSOGSA算法综合了粒子群算法(particle swarmoptimization,简称:PSO)的社会认知能力和万有引力搜索法(gravitational searchalgorithm,简称:GSA)的局部搜索能力,所述PSOGSA算法具体步骤如下:
其中,Vi(k)是第i个单元在进行第k次迭代时的速度,Xi(k)是第i个单元在进行第k次迭代时的位置,c'1和c'2是在上述算法的迭代过程中的加速度的系数,r'1和r'2是在[0,1]之间的随机数,gbest是在第k次迭代中的最优解,w为权值函数。
在上式中,aci(k)是第i个单元在第k此迭代时的加速度,所述加速度可以通过下式进行计算:
其中,是单元i(i≠j)对单元j的引力,rj是分布在[0,1]之间的随机数,时在t时刻单元i与单元j之间的引力,其中Maj为单元i的自重,Mpj为单元j的自重,G(t)是在t时刻的重力常数,ε是一个数值很小的常数,Rij(k)是单元i与单元j之间的欧氏距离,所述Rij(k)可以通过下述式子进行计算:
G(t)=G0×exp(-α×k/maxiter);
Rij(t)=dis(Xi(t),Xi(t));
其中α是衰减常数,G0是初始重力常数,maxiter是预先设置的最大迭代次数。
通过上述PSOGSA算法可以对预测模型中经过训练得到的参数进行进一步的优化。在优化过程中,需要先给出需要优化的参数的初始值,之后经过上述PSOGSA算法进行更新,直到达到最大迭代次数或者是达到要求的精度时,迭代停止,此时就得到了参数的最优值,此时预测模型为最优模型。
在本发明实施例中,可以通过对电力负荷数据进行归一化以及通过自相关分析选取相关性较大的电力负荷数据来提高模型的计算速度和预测精度。进一步的,本发明实施例中还可以通过PSOGSA算法对预测模型中的参数进行进一步优化,以此来提高所述预测模型计算出电力负荷预测数据的准确性。
在本发明实施例中,可以通过一个非线性函数对本发明实施例中的预测模型进行鲁棒性分析,所用的非线性函数如下:
在本发明实施例中,根据上述非线性函数,选取200个点作为训练数据对所述预测模型进行训练,其中每个训练数据的x都服从在[-10,10]之间的均匀分布。由于现在需要对预测模型的鲁棒性进行验证,需要对上述训练数据的y值进行[-1,1]之间的扰动;之后再取11001个点作为验证数据,用于验证预测模型所输出的数据是否准确,其中验证数据的x都服从在[-11,11]之间的均匀分布,y值为上述非线性函数的输出值。
有关预测模型的参数设置如下:
惩罚系数C=2^30,分位数τ=0.95,为保证预测可靠性,进行100次实验取平均值,预测模型的隐含层节点个数从5至80每递增5取一个值,输入的数据的异常值分别为10、20、50和80个。通过上述数据得到的预测模型输出的数据分布图如图3至图6所示,其中CQR-ORELM为本发明实施例所提供的预测模型的输出值的曲线。当输入的异常值为10个时、预测模型的隐含层节点个数最优为25个;当输入的异常值为20个时、预测模型的隐含层节点个数最优为70个;当输入的异常值为50个时、预测模型的隐含层节点个数最优为75个;当输入的异常值为80个时、预测模型的隐含层节点个数最优为35个。从图3至图6可以看出,本发明实施例所提供的预测模型对于异常值的鲁棒性较强,随着输入数据的异常值的个数从10增加到80时,本发明实施例所提供的预测模型依旧能保持良好的预测性能。
对于上述预测模型,我们还可以进一步的通过下述3个误差指标:MAE(平均绝对误差)、MAPE(平均值绝对百分误差)、RMSE(均方根误差) 在对本发明实施例所提供的预测模型进行评估,有关上述3个误差指标的计算公式如下:
其中,为t时刻的数据的真实值,为所述预测模型输出的数据在t时刻的计算值。有关计算结果参见下表:
通过上表同样可以看出,本发明实施例所提供的预测模型具有较好的鲁棒性,其中RMSE误差指标值随着输入数据中异常值的增加而呈现增长的趋势,说明异常值对模型的精度有影响,但是变化幅度不大,也进一步验证了模型具有较强的鲁棒性,即预测模型在异常值的干扰下仍然能够保持较高的预测精度。
上述在对预测模型鲁棒性的验证过程中,所出现的具体数据仅仅是对本发明起到说明的作用,仅为了可以更加方便的理解本发明,而不是对本发明实施例做具体的限定,上述出现的数值均可以有其他的替代方案,在此不做具体限定。
除了仅通过电力负荷数据进行预测,还可以增加实时电价来对电力负荷数据进行预测。电价的高低也直接影响到用户的用电习惯,所以在本方案中,还可以进一步的考虑实时电价对电力负荷预测数据所造成的影响。
请参考图7、图8、图9和图10,图7为自相关系数分布图;图8为MAE指标随候选变量个数变化的趋势图;图9为MAPE指标随候选变量个数变化的趋势图;图10为RMSE指标随候选变量个数变化的趋势图;
在本发明实施例中,预先收集了某电力公司公布于2016年7月的电力负荷数据以及实时电价,其中,数据的采样频率为30分钟,首先对上述数据进行归一化处理,处理方式参照上文电力负荷数据的归一化,再选取6至12日(周一至周日)的336个电力负荷样本数据进行自相关分析,其延迟时刻为335,得到如图7所示的自相关系数分布图。如图7所示,自相关系数出现了7个正相关的峰值点,对应每一天的样本中第一个采样的数据,其延迟时刻、自相关系数和所选取的候选数据如下表所示:
在本发明实施例中,由于添加了实时电价作为对电力负荷数据的影响条件,将分为两种情况进行说明,一种是考虑电价,一种是不考虑电价。在上述37个候选数据中,向预测模型输入不同个数的候选数据,并通过上述MAE、MAPE、RMSE三个误差指标对输出的电力负荷数据进行评价,得到图8至图10三个变化趋势图,其中横坐标为候选数据个数,纵坐标分别为三个误差指标的数值。
参见图8至图10,候选数据个数从1增加到37的过程中,各误差指标均呈现先迅速减少再趋于平缓的趋势,综合考虑上述3个误差指标,本发明提供的预测模型在考虑电价时,选取21个候选数据;在不考虑电价时,选取16个候选数据。
可以将本发明实施例所提供的预测模型CQR-ORELM(综合分位回归-鲁棒极限学习机)模型与现有的几种模型,例如:QR(分位回归)模型,ELM(极限学习机)模型和QR-SVM(分位回归-支持向量机)模型进行对比,通过上述方法先确定各个模型的最优候选数据个数:
再将上述候选数据输入到上述各个模型之中以生成电力负荷预测数据,其中对于QR模型,QR-SVM模型,CQR-ORELM模型取分位值τ分别为0.1、0.5和0.9,在上述3种分位值下分别通过各个模型计算电力负荷预测数据,并计算各组数据的MAE、MAPE、RMSE三个误差指标。在考虑电价时,计算结果如下表所示:
未考虑电价时,计算结果如下表所示:
通过上述两个表格,可以清晰的看出本发明所提供的预测模型对于其他现有的预测模型而言,在考虑电价和未考虑电价的情况下,均有较好的预测效果,并且相比于QR模型,QR-SVM模型和ELM模型,均显著提高了预测精度。
考虑电价时,本发明实施例所提供的预测模型,当τ=9.0时预测精度最好,其MAPE,MAE,RMSE三种误差指标值分别为0.3151%,16.7821MWh和21.0694MWh。与ELM模型预测结果相比较,三种误差指标值分别下降了2.597%,2.616%和3.662%。与QR-SVM模型的最优结果(即τ=0.9时)相比较,MAPE指标值下降了13.243%,而与QR模型的最优结果(即τ=0.9时)相比,RMSE指标值下降了19.264%。
未考虑电价时,本发明实施例所提供的CQR-ORELM模型的预测精度分别为0.2726%,14.5344MWh和19.8202MWh,在所有模型中,预测性能最优。
下面对本发明实施例提供的一种短时电力负荷预测的装置进行介绍,下文描述的短时电力负荷预测装置与上文描述的短时电力负荷预测方法可相互对应参照。
图11为本发明实施例提供的一种短时电力负荷预测装置的结构框图,参照图11,短时电力负荷预测装置可以包括:
获取模块100:用于获取历史电力数据;
调用模块200:用于调用预先建立的预测模型,根据所述历史电力数据计算出电力负荷预测数据,所述预测模型为综合分位回归和鲁棒极限学习机,并利用PSOGSA优化后所建立的混合预测模型;所述预测模型为通过预先获取的训练数据集训练得到的预测模型,所述训练数据集为根据历史电力数据得到的所述训练数据集。
在本发明实施例中,所述获取模块100可以具体用于:
获取电力数据,其中电力数据为历史电力负荷数据
在本发明实施例中,所述调用模块200可以具体用于:
调用预先建立的预测模型,根据所述电力数据计算出电力负荷预测数据,所述预测模型具体为:
其中Yt为所述预测模型的输出值,Qt为t时刻所述预测模型在τ分位下的输出值,为所述分位回归的影响参数,为所述鲁棒极限学习机的影响参数;相应的,训练所述预测模型时所用的损失函数为:
其中C是惩罚系数,β是所述鲁棒极限学习机中隐含层与输出层之间的权值,yt是训练所述预测模型时所使用的训练数据中t时刻的电力负荷值,ρτ为检验函数。
在本发明实施例中,所述调用模块200具体用于:
调用预先建立的预测模型,根据所述历史电力数据计算出电力负荷预测数据,所述预测模型为通过预先获取的训练数据集训练得到的预测模型,所述训练数据集具体为将预先获取的历史电力数据进行归一化处理,以生成归一化的样本数据集;通过自相关分析,提取所述归一化的样本数据集中各变量间的自相关系数;根据所述自相关系数和预先设定的选取规则,从所述归一化的样本数据集中选取可用的样本变量数据,以构成的所述训练数据集。
在本发明实施例中,所述调用模块200可以具体用于:
调用预先建立的预测模型,根据所述电力数据计算出电力负荷预测数据,所述预测模型中的参数是根据PSOGSA算法进行优化的参数。
本实施例的短时电力负荷预测装置用于实现前述的短时电力负荷预测方法,因此短时电力负荷预测装置中的具体实施方式可见前文中的短时电力负荷预测方法的实施例部分,例如,获取模块100,调用模块200,分别用于实现上述短时电力负荷预测方法中步骤S101和S102,所以,其具体实施方式可以参照相应的各个部分实施例的描述,在此不再赘述。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其它实施例的不同之处,各个实施例之间相同或相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
专业人员还可以进一步意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
以上对本发明所提供的一种短时电力负荷预测方法以及装置进行了详细介绍。本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以对本发明进行若干改进和修饰,这些改进和修饰也落入本发明权利要求的保护范围内。
Claims (8)
1.一种短时电力负荷预测的方法,其特征在于,所述方法包括:
获取历史电力数据;
调用预先建立的预测模型,根据所述历史电力数据计算出电力负荷预测数据,所述预测模型为综合分位回归和鲁棒极限学习机所建立的预测模型;所述预测模型为通过预先获取的训练数据集训练并利用PSOGSA进行优化得到的预测模型,所述训练数据集为根据所述历史电力数据得到的所述训练数据集。
2.根据权利要求1所述的方法,其特征在于,所述历史电力数据为历史电力负荷数据。
3.根据权利要求1所述的方法,其特征在于,所述预测模型具体为:
其中Yt为所述预测模型的输出值,Qt为t时刻所述预测模型在τ分位下的输出值,为所述分位回归的影响参数,为所述鲁棒极限学习机的影响参数;相应的,训练所述预测模型时所用的损失函数为:
其中C是惩罚系数,β是所述鲁棒极限学习机中隐含层与输出层之间的权值,yt是训练所述预测模型时所使用的训练数据中t时刻的电力负荷值,ρτ为检验函数。
4.根据权利要求3所述的方法,其特征在于,所述训练数据集具体
为:
将预先获取的所述历史电力数据进行归一化处理,以生成归一化的样本数据集;
通过自相关分析,提取所述归一化的样本数据集中各变量间的自相关系数;
根据所述自相关系数和预先设定的选取规则,从所述归一化的样本数据集中选取可用的样本变量数据,以构成所述训练数据集。
5.一种短时电力负荷预测的装置,其特征在于,所述装置包括:
获取模块:用于获取历史电力数据;
调用模块:用于调用预先建立的预测模型,根据所述历史电力数据计算出电力负荷预测数据,所述预测模型为综合分位回归和鲁棒极限学习机,并利用PSOGSA优化后所建立的预测模型;所述预测模型为通过预先获取的训练数据集训练得到的预测模型,所述训练数据集为根据所述历史电力数据得到的所述训练数据集。
6.根据权利要求5所述的装置,其特征在于,所述获取模块具体用于:
获取历史电力数据,其中历史电力数据为历史电力负荷数据。
7.根据权利要求5所述的装置,其特征在于,所述调用模块具体用于:
调用预先建立的预测模型,根据所述历史电力数据计算出电力负荷预测数据,所述预测模型具体为:
其中Yt为所述预测模型的输出值,Qt为t时刻所述预测模型在τ分位下的输出值,为所述分位回归的影响参数,为所述鲁棒极限学习机的影响参数;相应的,训练所述预测模型时所用的损失函数为:
其中C是惩罚系数,β是所述鲁棒极限学习机中隐含层与输出层之间的权值,yt是训练所述预测模型时所使用的训练数据中t时刻的电力负荷值,ρτ为检验函数。
8.根据权利要求7所述的装置,其特征在于,所述调用模块具体用
于:
调用预先建立的预测模型,根据所述历史电力数据计算出电力负荷预测数据,所述预测模型为通过预先获取的训练数据集训练得到的预测模型,所述训练数据集具体为将预先获取的历史电力数据进行归一化处理,以生成归一化的样本数据集;通过自相关分析,提取所述归一化的样本数据集中变量间的自相关系数;根据所述自相关系数和预先设定的选取规则,从所述归一化的样本数据集中选取可用的样本变量数据,以构成的所述训练数据集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710770943.6A CN107506868B (zh) | 2017-08-31 | 2017-08-31 | 一种短时电力负荷预测的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710770943.6A CN107506868B (zh) | 2017-08-31 | 2017-08-31 | 一种短时电力负荷预测的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107506868A true CN107506868A (zh) | 2017-12-22 |
CN107506868B CN107506868B (zh) | 2021-01-26 |
Family
ID=60694510
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710770943.6A Active CN107506868B (zh) | 2017-08-31 | 2017-08-31 | 一种短时电力负荷预测的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107506868B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108846517A (zh) * | 2018-06-12 | 2018-11-20 | 清华大学 | 一种分位数概率性短期电力负荷预测集成方法 |
CN109146063A (zh) * | 2018-08-27 | 2019-01-04 | 广东工业大学 | 一种基于重要点分割的多分段短期负荷预测方法 |
CN109656208A (zh) * | 2018-12-17 | 2019-04-19 | 七彩安科智慧科技有限公司 | 一种楼宇电力预警物联网的构建方法 |
CN109768989A (zh) * | 2019-02-27 | 2019-05-17 | 重庆邮电大学 | 基于lahp-igfnn的网络安全态势评估模型 |
CN111240282A (zh) * | 2019-12-31 | 2020-06-05 | 联想(北京)有限公司 | 一种工艺优化方法、装置、设备及计算机可读存储介质 |
CN111612139A (zh) * | 2019-02-22 | 2020-09-01 | 浙江宇视科技有限公司 | 一种神经网络模型训练方法、系统、设备及计算机介质 |
CN111697570A (zh) * | 2020-05-25 | 2020-09-22 | 国网河北省电力有限公司石家庄市藁城区供电分公司 | 一种电网负荷预测方法 |
CN112036081A (zh) * | 2020-08-26 | 2020-12-04 | 北京科技大学 | 基于收得率预测的转炉出钢硅锰合金加入量确定方法 |
CN112465250A (zh) * | 2020-12-08 | 2021-03-09 | 深圳供电局有限公司 | 电力负荷预测方法、装置、计算机设备和存储介质 |
CN113391256A (zh) * | 2021-05-28 | 2021-09-14 | 国网河北省电力有限公司营销服务中心 | 一种现场作业终端的电能表计量故障分析方法及系统 |
WO2022126809A1 (zh) * | 2020-12-15 | 2022-06-23 | 电子科技大学中山学院 | 一种以错因强化方式动态优化损失函数的方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105046374A (zh) * | 2015-08-25 | 2015-11-11 | 华北电力大学 | 一种基于核极限学习机模型的功率区间预测方法 |
-
2017
- 2017-08-31 CN CN201710770943.6A patent/CN107506868B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105046374A (zh) * | 2015-08-25 | 2015-11-11 | 华北电力大学 | 一种基于核极限学习机模型的功率区间预测方法 |
Non-Patent Citations (2)
Title |
---|
彭显刚等: "考虑负荷自适应检测和修复的鲁棒极限学习机短期负荷预测方法", 《中国电机工程学报》 * |
秦玥等: "基于分位数回归的电力负荷特性预测面板数据模型", 《电力需求侧管理》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108846517B (zh) * | 2018-06-12 | 2021-03-16 | 清华大学 | 一种分位数概率性短期电力负荷预测集成方法 |
CN108846517A (zh) * | 2018-06-12 | 2018-11-20 | 清华大学 | 一种分位数概率性短期电力负荷预测集成方法 |
CN109146063A (zh) * | 2018-08-27 | 2019-01-04 | 广东工业大学 | 一种基于重要点分割的多分段短期负荷预测方法 |
CN109656208A (zh) * | 2018-12-17 | 2019-04-19 | 七彩安科智慧科技有限公司 | 一种楼宇电力预警物联网的构建方法 |
CN111612139B (zh) * | 2019-02-22 | 2024-04-09 | 浙江宇视科技有限公司 | 一种神经网络模型训练方法、系统、设备及计算机介质 |
CN111612139A (zh) * | 2019-02-22 | 2020-09-01 | 浙江宇视科技有限公司 | 一种神经网络模型训练方法、系统、设备及计算机介质 |
CN109768989A (zh) * | 2019-02-27 | 2019-05-17 | 重庆邮电大学 | 基于lahp-igfnn的网络安全态势评估模型 |
CN111240282B (zh) * | 2019-12-31 | 2021-12-24 | 联想(北京)有限公司 | 一种工艺优化方法、装置、设备及计算机可读存储介质 |
CN111240282A (zh) * | 2019-12-31 | 2020-06-05 | 联想(北京)有限公司 | 一种工艺优化方法、装置、设备及计算机可读存储介质 |
CN111697570A (zh) * | 2020-05-25 | 2020-09-22 | 国网河北省电力有限公司石家庄市藁城区供电分公司 | 一种电网负荷预测方法 |
CN112036081A (zh) * | 2020-08-26 | 2020-12-04 | 北京科技大学 | 基于收得率预测的转炉出钢硅锰合金加入量确定方法 |
CN112036081B (zh) * | 2020-08-26 | 2021-06-01 | 北京科技大学 | 基于收得率预测的转炉出钢硅锰合金加入量确定方法 |
CN112465250A (zh) * | 2020-12-08 | 2021-03-09 | 深圳供电局有限公司 | 电力负荷预测方法、装置、计算机设备和存储介质 |
WO2022126809A1 (zh) * | 2020-12-15 | 2022-06-23 | 电子科技大学中山学院 | 一种以错因强化方式动态优化损失函数的方法 |
CN113391256A (zh) * | 2021-05-28 | 2021-09-14 | 国网河北省电力有限公司营销服务中心 | 一种现场作业终端的电能表计量故障分析方法及系统 |
CN113391256B (zh) * | 2021-05-28 | 2022-07-22 | 国网河北省电力有限公司营销服务中心 | 一种现场作业终端的电能表计量故障分析方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107506868B (zh) | 2021-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107506868A (zh) | 一种短时电力负荷预测的方法及装置 | |
US11043808B2 (en) | Method for identifying pattern of load cycle | |
Macedo et al. | Demand side management using artificial neural networks in a smart grid environment | |
Ghadimi et al. | PSO based fuzzy stochastic long-term model for deployment of distributed energy resources in distribution systems with several objectives | |
CN103092699A (zh) | 一种云计算资源预分配实现方法 | |
KR102346188B1 (ko) | 전력 수요 예측 방법 및 장치 | |
CN109636010A (zh) | 基于相关因素矩阵的省级电网短期负荷预测方法及系统 | |
CN104537433A (zh) | 基于存量容量与业扩报装特征的售电量预测方法 | |
CN108596242A (zh) | 基于小波神经网络和支持向量机的电网气象负荷预测方法 | |
CN110119948A (zh) | 基于时变权重动态组合的电力用户信用评价方法及系统 | |
CN105005711A (zh) | 获取统计线损的方法及装置 | |
CN110210670A (zh) | 一种基于电力系统短期负荷的预测方法 | |
CN106503833A (zh) | 基于支持向量机算法的光伏发电短期功率滚动预测方法 | |
CN106548285A (zh) | 计及小水电出力的趸售电量预测方法 | |
Kim et al. | Extracting baseline electricity usage using gradient tree boosting | |
Kasule et al. | Forecasting Uganda’s net electricity consumption using a hybrid PSO-ABC Algorithm | |
Karthik et al. | OELF: short term load forecasting for an optimal electrical load forecasting using hybrid whale optimization based convolutional neural network | |
CN109858667A (zh) | 一种基于雷电气候对负荷影响的短期负荷聚类方法 | |
CN109214610A (zh) | 一种基于长短期记忆神经网络的饱和电力负荷预测方法 | |
Bompard et al. | Dynamic price forecast in a competitive electricity market | |
Khosravi et al. | Wind farm power uncertainty quantification using a mean-variance estimation method | |
Ponoćko | Data analytics-based demand profiling and advanced demand side management for flexible operation of sustainable power networks | |
CN110232478A (zh) | 一种社会经济分析预测方法 | |
Wen et al. | Optimal operation framework of customer-premise battery storage for energy charge reduction and primary frequency regulation | |
CN108615091A (zh) | 基于聚类筛选和神经网络的电力气象负荷数据预测方法 |
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 |