CN112821383A - 一种基于深度学习的电力系统自然频率特性系数区间预测方法 - Google Patents
一种基于深度学习的电力系统自然频率特性系数区间预测方法 Download PDFInfo
- Publication number
- CN112821383A CN112821383A CN202110003572.5A CN202110003572A CN112821383A CN 112821383 A CN112821383 A CN 112821383A CN 202110003572 A CN202110003572 A CN 202110003572A CN 112821383 A CN112821383 A CN 112821383A
- Authority
- CN
- China
- Prior art keywords
- frequency characteristic
- natural frequency
- characteristic coefficient
- prediction
- beta
- 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 73
- 238000013135 deep learning Methods 0.000 title claims abstract description 19
- 238000012549 training Methods 0.000 claims abstract description 73
- 238000013507 mapping Methods 0.000 claims abstract description 5
- 230000006870 function Effects 0.000 claims description 33
- 210000002569 neuron Anatomy 0.000 claims description 33
- 230000004913 activation Effects 0.000 claims description 17
- 230000008859 change Effects 0.000 claims description 11
- 230000003068 static effect Effects 0.000 claims description 10
- 150000001875 compounds Chemical class 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 2
- 238000011217 control strategy Methods 0.000 abstract description 3
- 230000000875 corresponding effect Effects 0.000 description 13
- 238000004364 calculation method Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 239000013598 vector Substances 0.000 description 7
- 238000010606 normalization Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000000605 extraction Methods 0.000 description 5
- 108090000672 Annexin A5 Proteins 0.000 description 4
- 238000013528 artificial neural network Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 239000000126 substance Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000012952 Resampling Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000003292 diminished effect Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000007786 learning performance Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 238000002759 z-score normalization Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- 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/045—Combinations of networks
-
- 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
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/003—Load forecast, e.g. methods or systems for forecasting future load demand
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/241—The oscillation concerning frequency
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Power Engineering (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- General Business, Economics & Management (AREA)
- Biomedical Technology (AREA)
- Tourism & Hospitality (AREA)
- Artificial Intelligence (AREA)
- Computing Systems (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种基于深度学习的电力系统自然频率特性系数区间预测方法,步骤为:1)确定自然频率特性系数β的主要影响因素;2)利用DNN模型建立系统功率扰动、备用容量、机组启停方式与自然频率特性系数β间的映射关系;3)建立自然频率特性系数β的区间预测模型;4)获取Bootstrap训练集,并对DNN模型进行训练;5)估计自然频率特性系数β预测误差的方差,利用区间预测模型计算自然频率特性系数β预测结果在给定置信水平下的置信区间。本发明可作为AGC控制策略中频率偏差系数B整定的参考依据。
Description
技术领域
本发明涉及电力系统及其自动化领域,具体是一种基于深度学习的电力系统自然频率特性系数区间预测方法。
背景技术
当今社会,电力是人们生活中不可或缺的一部分。维持电力系统的安全、稳定、经济运行对于整个社会的发展而言至关重要。而电力系统的安全、经济运行则依赖于良好的频率质量,过高或过低的频率都将对整个系统带来严重危害,甚至引发大停电事故。自动发电控制(Automatic Generation Control,AGC)系统是控制互联电网频率的重要手段。在AGC控制过程中,一般通过计算区域控制偏差(Area Control Error,ACE)来确定系统有功调节量。AGC系统包含众多控制模式,不同控制模式具有不同的ACE计算方式。其中,联络线频率偏差控制(Tie-line frequency Bias Control)是目前应用最为广泛的控制模式。该模式的控制目标为将互联电网频率与联络线功率均控制在计划值附近。该模式下的ACE计算中,B系数是确定ACE的关键参数,其整定值决定了系统调整量对联络线功率偏差和频率偏差的偏向程度,将直接影响到频率控制性能的优劣,如式(1)所示。
ACE=(Ptie,a-Ptie,s)+B(fa-fs)=ΔPtie+BΔf (1)
式中,Ptie,a、Ptie,s为联络线功率实际值和计划值;fa、fs为频率实际值与计划值;B为人为整定的频率偏差系数;ΔPtie、Δf分别表示联络线功率偏差和频率偏差。
理想情况下,应将B系数整定的与电力系统自然频率特性β系数相等,这样可使得互联电网中每个控制区的AGC仅负责本区域的负荷扰动。当B系数略大于β系数时,系统有功调节量将大于系统功率偏差,可使得系统频率更快地恢复到计划值,但过大的B系数将使有功调节量及相关运行费用大幅增加。当B系数小于β系数时,将引起系统欠调甚至反调,危害系统频率质量。综上可见,B系数的整定原则是使其近似跟踪β系数,可略大于β系数,但应避免其小于β系数的的情形。
目前,我国工业界主要采用固定系数法来整定B系数,即取为年最大负荷的1%~2%,并每年调整一次。然而,β系数受负荷波动、备用容量、机组启停方式等系统运行状态的影响,是非线性时刻变化的。因此,固定系数法难以满足B≈β,将引起系统的超调、欠调甚至反调。为使B系数更好地跟踪β系数,国内外学者提出了分段B系数整定法,包括:以系统频率偏差的大小为分段依据的两段式B系数整定,以保障B>β;将频率偏差量划分为4个等级,随着频率偏差量的增大,不断增大B系数的整定值,以确保B>β;根据系统运行方式建立了全天分时段B系数整定模型,使其更好地近似β系数;根据火电、水电机组一次调频死区设置情况而建立的三段式B系数整定方法。上述分段B系数整定法可在一定程度上缓解系统欠调、反调现象,但仍不能较好地满足B≈β。为更好地跟踪β系数的变化,有学者提出时变B系数整定方法,即通过实时在线计算β系数来整定B系数。然而,该方法对数据的实时采集处理以及调度自动化水平提出了高要求。此外,由于β系数的在线计算需要在扰动发生后方可进行,故该方法控制时延较高,不利于频率的快速恢复。对此,有文献提出了β系数点预测方法,可根据系统功率扰动量对β系数的取值进行预测,为B系数的整定提供了参考。然而该方法仅考虑了功率扰动对β系数的影响,且预测误差的存在可能使得B<β的情况发生。
发明内容
本发明的目的是提供一种基于深度学习的电力系统自然频率特性系数区间预测方法,包括以下步骤:
1)确定自然频率特性系数β的主要影响因素。
所述自然频率特性系数β用于反映稳定状态下电力系统有功功率和频率的关系,包括负荷静态频率特性KL和发电机静态频率特性KG。
自然频率特性系数β如下所示:
式中,ΔP为系统功率偏差。Δf为频率偏差。
所述自然频率特性系数β的主要影响因素包括系统功率扰动、备用容量和机组启停方式。
2)利用DNN模型建立影响因素与自然频率特性系数β间的映射关系。
所述DNN模型的输入为自然频率特性系数β的影响因素,输出为自然频率特性系数β。
所述DNN模型包括输入层、若干隐藏层和输出层。层与层之间以全连接的方式连接在一起。
层与层之间的数据传递公式如下所示:
al=s(Wlal-1+bl) (2)
式中,al代表第l层神经元的输出。al-1代表第l-1层神经元的输出。l=1,2,…n。n为DNN的总层数。a0表示输入数据x。Wl为连接第l-1层和第l层神经元的权值矩阵,bl为第l层神经元的偏置项。函数s(·)为激活函数。
式中,an代表第n层神经元的输出。an-1代表第n-1层神经元的输出。Wn为连接第n-1层和第n层神经元的权值矩阵。bn为第n层神经元的偏置项。
隐藏层的激活函数如下所示:
s(t)=max(0.01t,t) (4)
输出层的激活函数如下所示:
s(t)=t (5)
DNN模型的参数θ={W,b}迭代更新如下:
θ(k)=θ(k-1)-Δθ(k) (6)
式中,θ(k-1)是在第k次更新之前的参数。Δθ(k)是在第k次更新时参数的变化量。θ(k)是第k次后的参数。k为迭代更新次数。
其中,第k次更新时参数的变化量Δθ(k)如下所示:
式中,dθ(k)是第k次更新时θ的梯度。r(k)是在第k次更新时参数的平方梯度的移动平均数。⊙代表Hadamard乘积。η(k)是在第k次更新时的学习率。δ为常数。
第k次更新时θ的梯度dθ(k)、第k次更新时参数的平方梯度的移动平均数r(k)、第k次更新时的学习率η(k)分别如下所示:
r(k)=ρr(k-1)+(1-ρ)dθ(k)⊙dθ(k) (9)
3)建立自然频率特性系数的区间预测模型。
建立自然频率特性系数β区间预测模型的步骤包括:
3.1)建立自然频率特性系数β的误差分解方程,即:
其中,模型误差εm和数据误差εd相互独立且服从期望值为0的正态分布,即:
自然频率特性系数β的总预测误差ε=εm+εd服从期望为0的正态分布,即:
3.2)给定置信水平100(1-α)%,其中α为显著性水平。
3.3)计算自然频率特性系数β在给定置信水平下的置信区间,即:
式中,Lα、Uα分别为置信水平为100(1-α)%时,自然频率特性系数β的下限和上限。z1-α/2为标准正态分布的(1-α/2)分位数。
4)获取Bootstrap训练集,并对DNN模型进行训练。
所述区间预测方法包括N个通过Bootstrap训练集训练得到的DNN模型。每个DNN模型的Bootsrap训练集相互独立。所述Bootstrap训练集为利用有放回的方式从m个训练样本中进行m次抽样得到的训练集。所述训练样本的数据包括系统功率扰动、备用容量、机组启停方式及自然频率特性系数。
计算自然频率特性系数β预测结果置信区间的步骤包括:
5.4)计算第i个输入数据对应的自然频率特性系数β的上下限,即:
式中,Lα(xi)为第i个输入数据对应的自然频率特性系数β的预测下限。Uα(xi)为第i个输入数据对应的自然频率特性系数β的预测上限。
值得说明的是,本发明首先分析自然频率特性系数β的主要影响因素;然后利用DNN强大的非线性特征提取能力建立系统功率扰动、备用容量、机组启停方式与自然频率特性系数β间的映射关系,实现β系数的点预测;最后结合Bootstrap方法得到β系数预测结果的置信区间。
本发明的技术效果是毋庸置疑的,本发明提出的基于DNN的自然频率特性系数点预测模型利用了DNN模型强大的非线性特征提取能力,可实现β系数的准确预测,能够有效追踪β系数的变化情况。
本发明所提结合Bootstrap的自然频率特性系数区间预测方法能够得到给定置信度下的β系数预测上下限,可为频率偏差B系数的整定提供参考依据。若将B系数整定为β系数的预测上限值,则可极大程度地避免B系数小于β系数,大大提高电力系统的安全性和经济性。
本发明可广泛应用于电力系统自然频率特性系数β区间预测,可作为AGC控制策略中频率偏差系数B整定的参考。
附图说明
图1为自然频率特性系数随功率扰动的变化情况;
图2为自然频率特性系数随备用容量的变化情况;
图3为DNN的结构;
图4为PowerFactory仿真界面。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
一种基于深度学习的电力系统自然频率特性系数区间预测方法,包括以下步骤:
1)确定自然频率特性系数β的主要影响因素。
所述自然频率特性系数β用于反映稳定状态下电力系统有功功率和频率的关系,包括负荷静态频率特性KL和发电机静态频率特性KG。
自然频率特性系数β如下所示:
式中,ΔP为系统功率偏差。Δf为频率偏差。
所述自然频率特性系数β的主要影响因素包括系统功率扰动、备用容量和机组启停方式。在分析过程中,上述主要因素可由如下方式获得:系统功率扰动量可通过分析负荷实时变化获得,负荷实时变化可通过广义相量测量系统(WAMS)得到,机组备用容量和启停方式也可由WAMS系统实时得到。自然频率特性系数β也可基于WAMS系统的频率数据、负荷数据等计算得到。在对自然频率特性系数β进行预测的过程中,需将系统未来功率扰动作为输入,可通过预测得到,后续实施方案将验证本发明所提方法对功率扰动预测误差的鲁棒性。而对于机组备用容量和启停方式而言,短时间内不会发生较大变动,可用当前时刻实测数据作为自然频率特性系数β预测的输入。
系统功率扰动的值和自然频率特性系数β的值正相关;当系统功率扰动的值大于扰动阈值后,自然频率特性系数β稳定;
备用容量的值和自然频率特性系数β的值正相关;当备用容量的值大于备用容量阈值后,自然频率特性系数β稳定;
表1列出了固定功率扰动下(20MW),机组启停方式对自然频率特性系数的影响。由表1可见,当有机组停运时,由于在线机组所能提供的一次调频能力减弱,β系数将减小。且由于不同机组的一次调频能力存在差异,不同机组的停运将对β系数带来不同的影响。表1中,机组1和机组2为同类型机组,故机组1/机组2停运时,β系数相同。机组3和机组4同理。
表1机组启停方式对自然频率特性系数的影响
2)利用DNN模型建立影响因素与自然频率特性系数β间的映射关系。
所述DNN模型的输入为自然频率特性系数β的主要影响因素,输出为自然频率特性系数β。
所述DNN模型包括输入层、若干隐藏层和输出层。层与层之间以全连接的方式连接在一起。
层与层之间的数据传递公式如下所示:
al=s(Wlal-1+bl) (2)
式中,al代表第l层神经元的输出。al-1代表第l-1层神经元的输出。l=1,2,…n。n为DNN的总层数。a0表示输入数据x。Wl为连接第l-1层和第l层神经元的权值矩阵,bl为第l层神经元的偏置项。函数s(·)为激活函数。
式中,an代表第n层神经元的输出。an-1代表第n-1层神经元的输出。Wn为连接第n-1层和第n层神经元的权值矩阵。bn为第n层神经元的偏置项。
隐藏层的激活函数如下所示:
s(t)=max(0.01t,t) (4)
输出层的激活函数如下所示:
s(t)=t (5)
当激活函数应用于公式2时,激活函数中的t指代Wlal-1+bl,应用于公式3时,激活函数中的t指代Wnan-1+bn。
DNN模型的参数θ={W,b}迭代更新如下:
θ(k)=θ(k-1)-Δθ(k) (6)
式中,θ(k-1)是在第k次更新之前的参数。Δθ(k)是在第k次更新时参数的变化量。θ(k)是第k次后的参数。k为迭代更新次数。W表示权重,b表示偏置。
其中,第k次更新时参数的变化量Δθ(k)如下所示:
式中,dθ(k)是第k次更新时θ的梯度。r(k)是在第k次更新时参数的平方梯度的移动平均数。⊙代表Hadamard乘积。η(k)是在第k次更新时的学习率。δ为常数。
第k次更新时θ的梯度dθ(k)、第k次更新时参数的平方梯度的移动平均数r(k)、第k次更新时的学习率η(k)分别如下所示:
r(k)=ρr(k-1)+(1-ρ)dθ(k)⊙dθ(k) (9)
3)建立自然频率特性系数的区间预测模型。
建立自然频率特性系数β区间预测模型的步骤包括:
3.1)建立自然频率特性系数β的预测方程,即:
其中,模型误差εm和数据误差εd相互独立且服从期望值为0的正态分布,即:
β系数的总预测误差ε=εm+εd服从期望为0的正态分布,即:
3.2)给定置信水平100(1-α)%,其中α为显著性水平。
3.3)计算自然频率特性系数β在给定置信水平下的置信区间,即:
式中,Lα、Uα分别为置信水平为100(1-α)%时,自然频率特性系数β的下限和上限。z1-α/2为标准正态分布的(1-α/2)分位数。
4)获取Bootstrap训练集,并对DNN模型进行训练。
所述区间预测方法包括N个通过Bootsrap训练集训练得到的DNN模型。每个DNN模型的Bootsrap训练集相互独立。所述Bootsrap训练集为利用有放回的方式从m个训练样本中进行m次抽样得到的训练集。所述训练样本的数据包括系统功率扰动、备用容量、机组启停方式和自然频率特性系数。机组启停方式为0-1变量,1表示机组开启,0表示机组关停。
计算自然频率特性系数β预测结果置信区间的步骤包括:
5.4)计算第i个输入数据对应的自然频率特性系数β的上下限,即:
式中,Lα(xi)为第i个输入数据对应的自然频率特性系数β的预测下限。Uα(xi)为第i个输入数据对应的自然频率特性系数β的预测上限。
6)将AGC控制策略中频率偏差B系数整定为自然频率特性系数β的预测上限,据此可在系统实际运行过程中计算区域控制偏差ACE:
ACE=(Ptie,a-Ptie,s)+B(fa-fs)=ΔPtie+BΔf (20)
式中,Ptie,a、Ptie,s为联络线功率实际值和计划值;fa、fs为频率实际值与计划值;B为频率偏差系数,其值整定为自然频率特性系数β的上限Uα(xi);ΔPtie、Δf分别表示联络线功率偏差和频率偏差。
根据ACE的计算值可确定电力系统有功调节量,实现电力系统有功功率的合理调节。
采用本方法整定频率偏差B系数的优势在于能有效避免B<β的情况发生,从而避免出现欠调或反调现象,有利于系统频率稳定性。此外,由于本发明所提区间预测方法所得自然频率特性系数β的置信区间较窄,将频率偏差B系数整定为自然频率特性系数β的预测上限后,所得系统ACE可较好地反映系统实际功率缺额,进而确定合理的系统功率调节量,避免出现过调现象。
实施例2:
一种基于深度学习的电力系统自然频率特性系数区间预测方法,包括以下步骤:
1)建立DNN模型,对电力系统自然频率特性系数进行预测;
2)结合Bootstrap方法,计算电力系统自然频率特性系数预测结果的置信区间。
实施例3:
一种基于深度学习的电力系统自然频率特性系数区间预测方法,主要步骤见实施例2,其中,计算电力系统自然频率特性系数预测结果置信区间的步骤包括:
1)确定自然频率特性系数及其影响因素
电力系统自然频率特性是电网的固有特性,反映了稳定状态下系统有功功率和频率的关系,它包括负荷静态频率特性及发电机静态频率特性,常用自然频率特性系数β来表示,如式(1)所示。
其中,β为系统自然频率特性系数;ΔP为系统功率偏差;KL、KG分别为负荷及发电机静态频率特性。
在电力系统运行过程中,β系数并非固定不变,它受系统功率扰动、备用容量、机组启停方式、负荷性质等因素的影响而表现出非线性和时变性。由于负荷静态频率特性对β系数的影响较小,本发明主要考虑了系统功率扰动、备用容量及机组启停方式对β系数的影响。上述相关因素对β系数的影响如图1、图2和表1所示,相关数据由PowerFactory仿真软件生成。
由图1可见,随着功率扰动的增大,β系数也随之增大。但当扰动增大到一定程度后,发电机组的一次调频能力趋于饱和,β系数的变化曲线也趋于平缓。
图2展示了固定功率扰动(60MW)下,β系数随备用容量的变化情况。其中,R表示系统备用容量。由图2可见,当功率扰动与备用容量较为接近时,由于扰动大小逼近了机组调节能力的上限,此时β系数较小。随着备用容量的增大,β系数随之增大。当备用容量增大到一定程度后,系统所留备用容量已足以补偿系统功率扰动,故β系数基本保持不变。
表1列出了固定功率扰动下(20MW),机组启停方式对自然频率特性系数的影响。由表1可见,当有机组停运时,由于在线机组所能提供的一次调频能力减弱,β系数将减小。且由于不同机组的一次调频能力存在差异,不同机组的停运将对β系数带来不同的影响。其中,机组1和机组2为同类型机组,故机组1/机组2停运时,β系数相同。机组3和机组4同理。
表1机组启停方式对自然频率特性系数的影响
2)特征向量选择
特征向量选择即确定DNN的输入向量及输出向量。本发明构建DNN模型的目的在于实现β系数的预测。因此,DNN的输出即为实际的β系数。为实现β系数的准确预测,DNN的输入应为影响β系数的关键因素。由1)的分析可见,β系数的关键影响因素包括功率扰动量、系统备用容量以及机组开机方式。其中,功率扰动量和系统备用容量可直接用相应的数值表示。对于机组开机方式,本发明则采用0-1向量来表示,1表示机组运行,0表示机组停运。另外,需要说明的是,作为DNN输入的功率扰动量一般需要在扰动发生后才能获得。为提前预测未来某时刻β系数的取值,可采用该时刻功率扰动量的预测值作为DNN的输入。
确定好特征向量后,便可收集相应数据,并形成训练样本用于DNN模型的训练。
3)基于DNN模型的自然频率特性系数点预测
深度神经网络(DNN)是传统浅层神经网络的拓展,它是多隐藏层结构,具备强大的特征提取能力,可自动挖掘隐含在训练数据中的复杂非线性关系。在分析β系数影响因素的基础上,本发明进一步利用DNN构建了β系数点预测模型。
DNN的结构如图3所示。
DNN由输入层、多个隐藏层以及输出层构成。同时,层与层之间以全连接的方式连接在一起。层与层之间的数据传递公式如式(2)所示。
al=s(Wlal-1+bl) (2)
其中,al代表第l层神经元的输出,l=1,2,…n;n为DNN的总层数;数据的输入x则表示为a0;Wl为连接第l-1层和第l层神经元的权值矩阵,bl为第l层神经元的偏置项。则最后一层第n层神经元的输出an(即DNN的预测输出)可由下式得到。
式(2)和式(3)中,函数s(·)为激活函数,在隐藏层,本发明选用LReLU(LeakyRectified Linear Unit)作为激活函数,可有效提高DNN模型的学习性能。由式(4)可知,LReLU激活函数是一个分段函数,当输入值t小于或等于零时,输出值等于0.01t;当输入值t大于零时,输出值等于输入值t。
s(z)=max(0.01t,t) (4)
在输出层选用线性激活函数作为激活函数,可以很大程度上减少计算量,其表达式如下。
s(t)=t (5)
确定了DNN的基本框架后,便可通过训练样本对DNN进行训练。DNN的训练是利用有标签样本来调整参数θ={W,b},使得损失函数L最小,本发明选用的损失函数为均方差损失函数,如式(6)所示。
其中,m为样本数量,y为训练样本实际输出,即真实的β系数。
本发明选用的训练算法为结合学习率衰减策略的均方根传播算法。均方根传播算法可以保持参数的平方梯度的移动平均数不变,使每个参数的步长自适应地更新,从而加快了参数的收敛速度。此外,学习率衰减策略随着训练的进行逐渐降低学习率,避免了算法来回振荡。该算法的具体步骤如式(7)-(11)所示:
r(k)=ρr(k-1)+(1-ρ)dθ(k)⊙dθ(k) (8)
θ(k)=θ(k-1)-Δθ(k) (10)
式中,dθ(k)是第k次更新时θ的梯度,θ(k-1)是在第k次更新之前的系数,是在第k次更新时θ的目标函数的偏导;r(k)是在第k次更新时参数的平方梯度的移动平均数;ρ是衰减率,其值设为0.9;⊙代表Hadamard乘积;Δθ(k)是在第k次更新时参数的变化量;η(k)是在第k次更新时的学习率;δ是为了防止出现而设置的常数,其值设为10-8;θ(k)是在第k次更新后的参数;η0代表最开始的学习率,其值设为0.001;τ表示学习率的衰减率,其值设为0.95。
训练完成后,便形成了基于DNN的自然频率特性系数点预测模型。
4)自然频率特性系数区间预测
自然频率特性系数点预测模型可预测得到β系数的取值,为B系数的整定提供参考依据。然而,该预测模型必然存在误差,若直接将B系数整定为β系数的预测值,则可能出现B<β这一不利于频率恢复的情况。
对β系数的预测值而言,其误差包含两个部分:模型误差以及数据误差。模型误差指模型训练效果不佳而引起的误差,引起训练效果不佳的因素包括训练过程中陷入局部最优、训练样本有限等。数据误差指因训练数据存在噪声等数据质量问题引起的误差。故存在误差的预测值可表示为:
本发明假设模型误差与数据误差相互独立且服从期望值为0的正态分布:
现有文献已证明,即使误差的实际分布不服从正态分布,基于正态分布假设的误差分析仍可取得良好效果。故本发明关于误差服从正态分布的假设是合理可行的。
β系数的总预测误差ε为模型误差εm与数据误差εd的叠加,基于上述假设可知,总预测误差也服从期望为0的正态分布:
其中,Lα、Uα分别为置信水平为100(1-α)%时,β系数的预测下限和预测上限;z1-α/2为标准正态分布的(1-α/2)分位数。
5)Bootstrap方法及误差方差计算
为避免B<β的情况,本发明进一步结合Bootstrap方法,得到了给定置信度下β系数的预测上下限。若将B系数整定为β系数的上限值,则可极大程度地避免B<β的情况发生。
本发明采用Bootstrap方法来计算预测误差的方差。Bootstrap是一种重采样技术,其基本原理如下:假设共有m个训练样本,采用有放回的方式从m个训练样本中进行抽样,抽样m次,便可得到一个新的包含m个训练样本的训练集,称为Bootstrap训练集。这样重复N次,便可得到N个Bootsrap训练集。通过每一个Bootstrap训练集都可以训练一个DNN模型,一共可得到N个DNN模型。通过N个训练好的DNN模型,可分别对模型误差与数据误差的方差进行计算。
对于模型误差方差的计算,首先,将N个DNN所得β系数预测结果的平均值作为最终预测输出,则第i个样本的预测输出可由下式得到。
由此,可得到模型误差的方差为:
其中,μεd为数据误差的均值。
由模型误差方差和数据误差方差的计算结果可得到总预测误差方差,如式(19)所示。
于是,由式(15)便可得到第i个样本对应的β系数上下限,如式(20)所示。
实施例4:
一种基于深度学习的电力系统自然频率特性系数区间预测方法的验证试验,包括以下步骤:
1)训练样本获取
为验证本发明所提自然频率特性系数区间预测方法的有效性,本发明利用PowerFactory仿真软件,通过调整自然频率特性系数的各影响因素共生成了685个样本。其中600个作为训练样本,余下85个用作测试样本。PowerFactory仿真界面如图4所示。
2)基于Bootstrap方法的重采样
采用Bootstrap方法对训练集进行有放回的抽样,共抽样600次,便可得到一个新的包含600个训练样本的训练集,称为Bootstrap训练集。重复上述过程100次,最终形成了100个Bootstrap训练集(即Bootstrap训练集个数N设为100)。
3)样本预处理
样本中包含的系统功率扰动、备用容量、机组启停方式以及自然频率特性系数具有不同的量纲且数值上具有较大差别,这不利于DNN模型的训练,为此,需要对样本进行归一化预处理。本发明采用z-score方法对样本的输入数据和输出数据进行归一化预处理,如式(1)所示。z-score标准化方法利用样本均值和标准差来进行归一化处理,经过z-score处理的数据满足均值为0,标准差为1。
式中,x*为归一化后的样本值,x为待归一化的样本值,μ为样本均值,σ为样本标准差。
另外,值得说明的是,本发明中机组启停方式采用0-1向量来表示,无需进行归一化处理。因此,本发明仅对系统功率扰动、备用容量和自然频率特性系数进行归一化预处理。
4)DNN模型训练
对100个Bootstrap训练集进行归一化预处理后,便利用每个训练集对DNN进行训练,最终可训练得到100个DNN模型。DNN模型的训练过程中,迭代次数设为2000次。
5)预测误差方差计算及β系数区间预测
根据100个DNN模型的训练结果,可计算得到自然频率特性系数预测结果的模型误差方差、数据误差方差和总预测误差方差。将置信水平设定为95%,便可得到β系数的预测区间。
具体仿真结果如下:
I)训练方法对比
本发明将对比如下区间预测方法:
M1:基于DNN的Bootstrap方法,其中,DNN包含4个隐藏层,每层50个神经元;
M2:基于浅层反向传播(Backpropagation,BP)神经网络的Bootstrap方法,其中BP神经网络包含1个隐藏层,隐藏层神经元数为200;
M3:基于极限学习机(Extreme Machine Learning,ELM)的Bootstrap方法,其中ELM包含1个隐藏层,隐藏层神经元数为200。
本发明所用计算机硬件环境为Intel(R)Core(TM)i7-7500U CPU@2.70GHz 8GBRAM。
II)自然频率特性系数区间预测的准确性
为验证本发明所提自然频率特性系数区间预测方法M1的有效性,将通过如下几个指标来衡量M1-M3的区间预测结果的准确性。
β系数的点预测误差可用平均绝对百分比误差(Mean Absolute PercentageError,MAPE)来衡量,如式(2)所示。
其中,mt为测试样本的个数。
所得预测区间是否有效涵盖了β系数的真实值由预测区间覆盖率(PredictionInterval Coverage Probability,PICP)来衡量,如式(3)所示。
所得预测区间的宽度一般由平均预测区间宽度(Mean Prediction IntervalWidth,MPIW)来量化,如式(4)所示。为更直观展示预测区间宽度的大小,本发明进一步定义相对平均预测区间宽度(Relative Mean Prediction Interval Width,RMPIW),如式(5)所示。
一般而言,较好的预测结果应具有较小的MAPE,较小的RMPIW以及较高的PICP。M1-M3所得区间预测结果对比如表2所示。
表2区间预测结果对比
方法 | MAPE | RMPIW | PICP |
M1 | 0.86% | 6.22% | 97.65% |
M2 | 1.85% | 8.89% | 91.76% |
M3 | 3.33% | 13.86% | 92.94% |
由表2可见,在隐藏层神经元总数均相同的情况下,具有深层结构的M1能够取得最优的预测结果,其各项指标均优于仅具有浅层结构的M2及M3,验证了DNN在特征提取方面的有效性。
3)自然频率特性系数区间预测结果对功率扰动量预测误差的鲁棒性
本发明所用特征向量包含系统功率扰动量。由于在功率扰动发生前,该扰动量是未知的,故在对未来β系数进行预测时,所用输入为功率扰动量的预测值。由于功率扰动量的预测存在误差,将预测值作为输入必然影响β系数区间预测的准确性。为验证本发明所提区间预测方法M1对功率扰动量预测误差的鲁棒性,将通过在真实功率扰动量的基础上添加服从正态分布的预测误差来生成功率扰动量的预测值,并测试了M1在不同功率扰动量误差情况下的预测效果,如表3所示。
表3 M1对功率扰动量预测误差的鲁棒性
扰动量误差大小 | MAPE | RMPIW | PICP |
1% | 1.10% | 6.68% | 96.47% |
2% | 1.57% | 9.28% | 96.47% |
3% | 1.61% | 10.30% | 95.29% |
4% | 1.99% | 11.32% | 94.12% |
5% | 2.21% | 12.26% | 95.29% |
由表3可见,随着功率扰动量预测误差的增大,区间预测的准确性呈下降趋势。但在误差为1%-5%的范围内,M1仍能保持较高的预测区间覆盖率,且可将相对平均预测区间宽度保持在13%以内,预测精度能够满足实际需求,可为B系数的整定提供参考。由此可见,本发明所提自然频率特性系数区间预测方法M1对于功率扰动量预测误差具备良好的鲁棒性。
实验结果表明,本发明所提基于DNN和Bootstrap的电力系统自然频率特性系数区间预测方法能够以高精度追踪β系数的变化,所得预测区间能够有效涵盖β系数的真实值,且该方法对功率扰动量预测误差具备良好的鲁棒性。
本发明在总结分析β系数主要影响因素的基础上,提出了基于DNN及Bootstrap的β系数区间预测方法。本发明首先确定了自然频率特性系数β的三个主要影响因素:系统功率扰动、备用容量、机组启停方式;然后利用DNN强大的非线性特征提取能力,实现β系数的准确点预测;最后,基于DNN点预测模型及Bootstrap方法得到了β系数的置信区间。β系数的预测上限值可为B系数的整定提供参考,在很大程度上避免B系数整定值小于β系数的情况发生。此外,由算例仿真可见,本发明所提区间预测方法能够取得良好的预测精度,且对功率扰动量预测误差具备良好的鲁棒性,验证了本发明所提方法的有效性。
Claims (10)
3.根据权利要求1所述的一种基于深度学习的电力系统自然频率特性系数区间预测方法,其特征在于:所述自然频率特性系数β的主要影响因素包括系统功率扰动、备用容量和机组启停方式。
4.根据权利要求1所述的一种基于深度学习的电力系统自然频率特性系数区间预测方法,其特征在于:所述DNN模型的输入数据包括自然频率特性系数β的影响因素,输出为自然频率特性系数β。
5.根据权利要求1所述的一种基于深度学习的电力系统自然频率特性系数区间预测方法,其特征在于:所述DNN模型包括输入层、若干隐藏层和输出层;层与层之间以全连接的方式连接在一起;
层与层之间的数据传递公式如下所示:
al=s(Wlal-1+bl) (2)
式中,al代表第l层神经元的输出;al-1代表第l-1层神经元的输出;l=1,2,…n;n为DNN的总层数;a0表示输入数据x;Wl为连接第l-1层和第l层神经元的权值矩阵,bl为第l层神经元的偏置项;函数s(·)为激活函数;
式中,an代表第n层神经元的输出;an-1代表第n-1层神经元的输出;Wn为连接第n-1层和第n层神经元的权值矩阵;bn为第n层神经元的偏置项。
6.根据权利要求5所述的一种基于深度学习的电力系统自然频率特性系数区间预测方法,其特征在于:隐藏层的激活函数如下所示:
s(t)=max(0.01t,t) (4)
输出层的激活函数如下所示:
s(t)=t (5)。
7.根据权利要求5所述的一种基于深度学习的电力系统自然频率特性系数区间预测方法,其特征在于:DNN模型的参数θ={W,b}迭代更新如下:
θ(k)=θ(k-1)-Δθ(k) (6)
式中,θ(k-1)是在第k次更新之前的参数;Δθ(k)是在第k次更新时参数的变化量;θ(k)是第k次后的参数;k为迭代更新次数;
其中,第k次更新时参数的变化量Δθ(k)如下所示:
式中,dθ(k)是第k次更新时θ的梯度;r(k)是在第k次更新时参数的平方梯度的移动平均数;⊙代表Hadamard乘积;η(k)是在第k次更新时的学习率;δ为常数;
第k次更新时θ的梯度dθ(k)、第k次更新时参数的平方梯度的移动平均数r(k)、第k次更新时的学习率η(k)分别如下所示:
r(k)=ρr(k-1)+(1-ρ)dθ(k)⊙dθ(k) (9)
8.根据权利要求1所述的一种基于深度学习的电力系统自然频率特性系数区间预测方法,其特征在于,建立自然频率特性系数β区间预测模型的步骤包括:
1)建立自然频率特性系数β的误差分解方程,即:
其中,模型误差εm和数据误差εd相互独立且服从期望值为0的正态分布,即:
自然频率特性系数β的总预测误差ε=εm+εd服从期望为0的正态分布,即:
2)给定置信水平100(1-α)%,其中α为显著性水平;
3)计算自然频率特性系数β在给定置信水平下的置信区间,即:
9.根据权利要求1所述的一种基于深度学习的电力系统自然频率特性系数区间预测方法,其特征在于,所述区间预测方法包括N个通过Bootstrap训练集训练得到的DNN模型;每个DNN模型的Bootstrap训练集相互独立;所述Bootstrap训练集为利用有放回的方式从m个训练样本中进行m次抽样得到的训练集;所述训练样本的数据包括系统功率扰动、备用容量、机组启停方式以及对应自然频率特性系数。
10.根据权利要求1所述的一种基于深度学习的电力系统自然频率特性系数区间预测方法,其特征在于,利用区间预测模型计算自然频率特性系数β预测结果置信区间的步骤包括:
4)计算第i个输入数据对应的自然频率特性系数β的上下限,即:
式中,Lα(xi)为第i个输入数据对应的自然频率特性系数β的预测下限;Uα(xi)为第i个输入数据对应的自然频率特性系数β的预测上限。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110003572.5A CN112821383A (zh) | 2021-01-04 | 2021-01-04 | 一种基于深度学习的电力系统自然频率特性系数区间预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110003572.5A CN112821383A (zh) | 2021-01-04 | 2021-01-04 | 一种基于深度学习的电力系统自然频率特性系数区间预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112821383A true CN112821383A (zh) | 2021-05-18 |
Family
ID=75857019
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110003572.5A Pending CN112821383A (zh) | 2021-01-04 | 2021-01-04 | 一种基于深度学习的电力系统自然频率特性系数区间预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112821383A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114066089A (zh) * | 2021-11-25 | 2022-02-18 | 中国工商银行股份有限公司 | 批量作业运行耗时区间确定方法及装置 |
CN114530889A (zh) * | 2022-03-28 | 2022-05-24 | 四川大学 | 基于大扰动区间识别的光伏系统双模式频率控制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105808889A (zh) * | 2016-04-25 | 2016-07-27 | 中国电力科学研究院 | 一种频率偏差系数仿真配置方法 |
CN110232409A (zh) * | 2019-05-30 | 2019-09-13 | 国网上海市电力公司 | 一种配网跳闸故障类型自动识别方法 |
-
2021
- 2021-01-04 CN CN202110003572.5A patent/CN112821383A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105808889A (zh) * | 2016-04-25 | 2016-07-27 | 中国电力科学研究院 | 一种频率偏差系数仿真配置方法 |
CN110232409A (zh) * | 2019-05-30 | 2019-09-13 | 国网上海市电力公司 | 一种配网跳闸故障类型自动识别方法 |
Non-Patent Citations (3)
Title |
---|
CAN WAN,等: "Probabilistic Forecasting of Wind Power Generation Using Extreme Learning Machine", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 * |
MINGXU XIANG,等: "Probabilistic power flow with topology changes based on deep neural network", 《ELECTRICAL POWER AND ENERGY SYSTEMS》 * |
董飞飞: "基于静态频率稳定性分析的电网容量交换能力评估", 《陕西电力》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114066089A (zh) * | 2021-11-25 | 2022-02-18 | 中国工商银行股份有限公司 | 批量作业运行耗时区间确定方法及装置 |
CN114530889A (zh) * | 2022-03-28 | 2022-05-24 | 四川大学 | 基于大扰动区间识别的光伏系统双模式频率控制方法 |
CN114530889B (zh) * | 2022-03-28 | 2022-12-13 | 四川大学 | 基于大扰动区间识别的光伏系统双模式频率控制方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11409347B2 (en) | Method, system and storage medium for predicting power load probability density based on deep learning | |
CN110659779B (zh) | 一种基于长短期记忆网络的配电系统网损预测方法 | |
CN110826791A (zh) | 一种基于长短时记忆神经网络的混合风功率预测方法 | |
CN112821383A (zh) | 一种基于深度学习的电力系统自然频率特性系数区间预测方法 | |
CN112018758A (zh) | 基于数字孪生的含高比例新能源交直流混联系统建模方法 | |
CN110766320A (zh) | 一种机场智能电网运行安全评价方法及装置 | |
LU500551B1 (en) | Virtual load dominant parameter identification method based on incremental learning | |
CN111415068A (zh) | 基于改造措施与失负荷量指标关联性的配电决策建模方法 | |
CN116796644A (zh) | 基于多智能体sac深度强化学习的风电场参数辨识方法 | |
CN115358437A (zh) | 基于卷积神经网络的供电负荷预测方法 | |
CN112821424A (zh) | 一种基于数据-模型融合驱动的电力系统频率响应分析方法 | |
CN115329669A (zh) | 基于深度贝叶斯主动学习的电力系统暂态稳定评估方法 | |
CN111080000A (zh) | 基于psr-dbn的超短期母线负荷预测方法 | |
CN112821420B (zh) | 一种基于XGBoost的ASFR模型中动态阻尼因子、多维频率指标的预测方法及系统 | |
CN113890017A (zh) | 基于关键量测的配电网电压自适应控制方法 | |
WO2024082836A1 (zh) | 一种水风光储多能互补系统综合效益评估方案寻优方法 | |
CN110580578A (zh) | 一种智能变电站二次系统运行质量多分层评价方法 | |
CN114268102B (zh) | 基于解析式概率潮流模型的电力系统运行状态量化方法 | |
CN116823008A (zh) | 一种园区能源利用效率评估方法、系统、设备和存储介质 | |
CN114759579A (zh) | 一种基于数据驱动的电网有功优化控制系统、方法和介质 | |
CN115313422A (zh) | 主动支撑主网频率的电-气综合能源系统运行控制方法 | |
CN115099448A (zh) | 一种基于VMD-Prophet的短期负荷预测方法 | |
Yang et al. | Research on power system small signal stability analysis and correction based on LightGBM algorithm | |
Hu et al. | An online deep reinforcement learning based parameter identification method for HVDC system | |
Huang et al. | Probabilistic prediction intervals of wind speed based on explainable neural network |
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: 20210518 |
|
RJ01 | Rejection of invention patent application after publication |