CN102801180A - 一种微电网并网等效模型的参数辨识方法 - Google Patents
一种微电网并网等效模型的参数辨识方法 Download PDFInfo
- Publication number
- CN102801180A CN102801180A CN2012102104627A CN201210210462A CN102801180A CN 102801180 A CN102801180 A CN 102801180A CN 2012102104627 A CN2012102104627 A CN 2012102104627A CN 201210210462 A CN201210210462 A CN 201210210462A CN 102801180 A CN102801180 A CN 102801180A
- Authority
- CN
- China
- Prior art keywords
- prime
- equivalent
- axle
- power
- motor
- 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
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P80/00—Climate change mitigation technologies for sector-wide applications
- Y02P80/10—Efficient use of energy, e.g. using compressed air or pressurized fluid as energy carrier
- Y02P80/14—District level solutions, i.e. local energy networks
Landscapes
- Control Of Eletrric Generators (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种微电网并网等效模型的参数辨识方法,利用微电网并网接入点电压、功率数据,对等效模型进行特征研究和建模,得到适用于微电网并网仿真的等效模型参数,同时利用优化辨识方法建立等效模型的数学模型,通过最小化误差的平方和寻找数据的最佳匹配,得到适用于微电网并网仿真的等效模型参数。
Description
技术领域
发明涉及一种微电网并网等效建模方法,属于电力系统建模和控制技术领域。
背景技术
为了充分利用分布式发电,发挥其经济、环境效益,解决分布式发电高渗透给电网带来的威胁,微电网作为分布式发电高效利用的一种电网组织形式,为分布发电安全、高效利用提出了一种解决思路。微电网本质是小型电力系统(10MW以下),集分布式发电系统、储能系统、负荷于一体,具备孤岛运行和并网运行的能力,并且可以在两种运行状态之间平滑过渡。
并网模式下,由于容量小,单个微电网对于大电网的影响有限,而随着微电网数量、容量以及电压等级的提高,其对于配电网的影响将逐渐显现,大量微电网的接入必然给电力系统仿真分析和实时控制带来影响。
微电网中元件包括各类分布式发电、储能系统、负荷以及电力电子设备。相对大电网,微电网可以看成一整体可控的“电源”或者“负荷”,表现出一定的电源或者负荷特性。微电网内部的运行状态和运行方式相对较复杂,分析微电网和大电网相互作用机理时,并不需要关心微电网内部运行状态,而主要从微电网和大电网的接口处电压、频率和功率特性入手。因此,在分析微电网并网影响时,可以对微电网进行等效处理,利用微电网的等效模型来反映微电网接入时的动态特性。
发明内容
本发明所要解决的技术问题是针对将来大量微电网接入微电网和大电网相互作用问题,提供一种适用于微电网并网接入仿真分析的微电网等效模型并且得到适用于微电网并网仿真的等效模型参数。
为解决上述技术问题,本发明提供了一种微电网并网等效模型的参数辨识方法,其特征在于,包括以下步骤:
1)建立等效电机动态元件模型:
ud=rsid-X′qiq+E′d
(2)
uq=rsiq+X′did+E′q
Pm=udid+uqiq
(3)
Qm=udiq-uqid
建立等效静态元件模型:
(4)
Pm、Qm为等效电机元件功率,ud、uq接入点电压d、q轴分量,id、iq接入点电流d、q轴分量,rs为等效电机定子电阻,E′d,E′q为d、q轴暂态电势,Xd,Xq为d、q轴稳态电抗,X′d,Xq′为d、q轴暂态电抗,T′d0,T′q0为d、q轴开路时间常数,Ef、Eg为d、q轴转子电压,Tj为转子惯性时间常数,Xf为转子d轴电抗,Xad为d轴互抗,Xg为转子q轴电抗,Xaq为q轴互抗,Tm为机械转矩,Ps、Qs为静态元件功率,Ps0、Qs0为静态元件稳态功率,U微电网接入点母线电压幅值,U0为微电网接入点母线电压稳态幅值,pu、qu为静态元件电压有功和无功指数;t为时间,s为滑差,ωr为转子转速;ωs为系统公共参考轴的转速,均为标么值;
上述模型参数中包括等效电机元件模型基本参数、等效静态元件模型参数及关联参数,等效电机元件基本参数包括:定子电阻rs、定子漏抗Xsl、d轴互感抗Xad、q轴互感抗Xaq、转子d轴电阻rf、转子d轴漏抗Xfl,转子q轴电阻rg、转子q轴漏抗Xgl、转子惯性时间常数Tj;等效静态元件参数包括:静态元件电压有功指数pu、静态元件电压无功指数qu;关联参数包括:s0和Kmp,s0表示等效电机的初始滑差,Kmp表示为等效电机元件初始功率比例;
等效电机模型实用参数(实用参数就是实用模型的参数)和基本参数间存在如下对应关系:
Xf=Xfl+Xad
Xg=Xgl+Xaq
Xd=Xsl+Xad
Xq=Xsl+Xaq
2)输入量测数据,包括微电网接入点电压ud+juq和接入点有功功率和无功功率PM+jQM;
3)根据等效电机有功功率比例计算得到稳态时等效电机有功功率Pm=Kmp*PM,根据粒子群优化算法初始给定的模型参数计算等效电机动态方程状态变量的初始值:d、q轴暂态电势E′d0、E′q0,等效电机的初始滑差s0;
4)利用龙格库塔法求解等效电机动态元件模型方程,获得等效电机的动态输出功率值;
5)根据电压值和等效静态元件的初试功率,计算等效静态元件的输出功率值;
6)计算目标函数即微电网总体等效模型输出功率与量测功率的偏差
其中,E(θ)表示量测数据和辨识数据的功率偏差,PM(k,θ)、QM(k,θ)表示测量得到的微电网接入点有功功率和无功功率值,k表示测量的数据点数,θ表示优化变量集,P∑(k)、Q∑(k)表示模型计算得到的有功功率和无功功率值:包括等效电机和等效静态元件的功率;
7)利用粒子群算法调整模型参数进行优化计算,判断迭代次数最终获得模型参数值。
本发明所达到的有益效果:
本发明利用微电网并网接入点电压、功率数据,对等效模型进行特征研究和建模,得到适用于微电网并网仿真的等效模型参数,同时利用优化辨识方法建立等效模型的数学模型,通过最小化误差的平方和寻找数据的最佳匹配,得到适用于微电网并网仿真的等效模型参数。
等效模型能够反映微电网内部元件的机理特性,利用等效电机元件反映微电网中同步电机、异步风力发电机和异步感应电动机的电机旋转特性;利用等效静态元件反映微电网中静态负荷、光伏发电等具有静态特性元件的分布式发电系统外部动态特性。
附图说明
图1为微电网等效模型结构图;
图2为等效电机坐标矢量图;
图3为等效电机d、q轴等效电路;
图4为建模方程流程图。
具体实施方式
电力系统仿真计算是电力系统动态分析和安全控制的基本工具,也是电力生产部门用于指导电网运行的基本依据。不恰当的模型会使计算结果与实际情况不一致,从而构成系统的潜在威胁或造成不必要的浪费。电力系统四大模型的建模方法已经得到了广泛的研究,随着分布式可再生能源发电技术的发展和成熟,微电网作为分布式发电开发利用的重要电网组织形式,在未来智能电网建设框架下必然得到更广泛的认同和重视。因此,关于微电网的建模将作为一重要分支而需要得到深入的研究,包括微电网内部元件的详细建模以及微电网并网仿真研究的等效建模。
本发明在微电网元件建模的基础上,提炼微电网元件模型之间的共性,提出来一种适合于微电网并网仿真研究的微电网等效模型,等效模型包括等效电机元件和等效静态元件。模型的推导如下所示:
定义等效电机各坐标关系如图2所示,其中abc为电机定子坐标,ABC为电机转子坐标,dq为转速为同步速ωs的旋转坐标,xy为转速为同步速ωs的系统公共坐标。假定在0时刻,定子a轴和d轴重合,则t时刻d轴领先于a轴角度为θs=ωst,转子绕组领先于定子绕组的角度为θ=ωrt,d轴领先A轴θr=(ωs-ωr)t。
经过Park变换得到在dq坐标中的定子电压方程
uds=pψds-ωsψqs+rsids (6)
uqs=pψqs+ωsψds+rsiqs
转子电压方程
uf=pψf-ωssψg+rrif (7)
ug=pψg+ωssψf+rrig
此时转子电路也分为d,q轴,其中d就是传统上的f,q轴视为g,为了和同步发电机方程对应,在此用下标用f轴和g表示。其磁链与电流的关系如下所示
定子磁链
ψds=Ldid+Ladif
(8)
ψqs=Lqiq+Laqig
转子磁链
ψfr=Ladid+Lfif
(9)
ψgr=Laqiq+Lgig
经过Park变化后的等效电机d,q轴等效电路如图3所示。
定义,并忽略定子暂态pψds=pψqs=0
此时转子电压方程(7)可以写成
(11)
定子磁链方程带入定子电压方程(6),并忽略定子暂态得到
(12)
将电感转化为电抗X′q=ωsL′q,X′d=ωsL′d,(12)式可以表示为
ud=rsid-X′aiq+ωsE′d
(13)
uq=rsiq+X′did+ωsE′q
转子运动方程
等效静态元件理解为其输出功率为电压和频率的代数方程,主要表示微电网中静态负荷和静态电源。利用幂指数来描述静态部分,如下所示。
由于微电网并网运行,此时微电网的频率完全由大电网来维持,因此,模型中频率的标幺值始终为1,并且电感和电抗的标幺值相等。此时模型方程可以写成:
电压方程
ud=rsid-X′qiq+E′d
(17)
uq=rsiq+X′did+E′q
Pm=udid+uqiq
(18)
Qm=udiq-uqid
静态元件模型
模型参数及辨识
由模型(16)、(17)、(18)、(19)构成的微电网等效模型方程可知,模型包含各种参数,因此需要通过相应的建模方法来确定模型参数。等效模型的参数确定为基本参数:Kmp,s0,rs,Xsl,Xad,Xaq,rf,Xfl,rg,Xgl,Tj,pu,qu,s0、Kmp。
模型参数辨识是利用电压、频率和功率等量测数据,通过辨识方法获得模型参数值。将参数辨识转化为参数优化问题,利用粒子群优化方法来求解。用每一个粒子代表一组模型参数X,将每组参数代入模型中,求解微分方程,计算出每一点的输出量,把这些计算出的输出量与测量值相比较,得到相应的误差平方和指标E(X)。具体辨识过程如下:
(1)根据稳态的母线电压ud0+juq0和测量功率PM+jQM、每一个粒子表示的等效模型参数由此进行初始化计算求得方程各状态变量的初始值E′d0,E′q0,s0。
初始电动势E′d0,E′q0计算过程
(20)
ud0=rsid0-X′qiq0+ωs0E′d0
(21)
uq0=rsiq0+X′did0+ωs0E′q0
可得
令
最终得到状态变量的初始解
id0=Hrs[ud0-E′d0]+HX′q[uq0-E′q0]
(24)
iq0=Hrs[uq0-E′q0]-HX′d[ud0-E′d0]
由于电机参数标幺值的基准值为电机容量,而系统电压功率的标幺值为系统容量,因此模型求解过程中所有参数需要统一到基准值下面进行求解。电机容量和系统容量之间的转换公式为:
其中SBM为电机自身容量基准值,SBS为系统容量基准值,KmpP∑为电机系统容量下初始有功功率,PMM0为电机在自身容量下的稳态有功功率,可以通过电机参数和端口电压计算得到。
(2)根据初始化计算过程中求得的各状态量的初值,运用四阶龙格库塔法求解微分方程组得到每一时步的E′di,E′qi,si和电动机消耗的功率PMi+jQMi,1≤i≤N,N为测量数据的长度。
龙格库塔法是求解微分方程的有效方法,利用四阶龙格库塔法对上述模型的微分方程进行求解
龙格库塔法没迭代一次需要四次计算函数值,其精度很高,但其步长可以变化,可以设计成变步长微分方程求解方法。在微分方程求解过程中所涉及的变量和参数都应该是标幺值,因此,步长时间h做标幺化处理。h*=kωB,步长的标幺值为采样时间间隔和角速度基准值的乘积,角速度的基准值ωB=2πfB=100π。
(3)根据公式(18)计算等效电机有功功率Pm和无功功率Qm。
(4)根据公式(19)计算等效静态元件有功功率Ps和无功功率Qs
(5)目标函数计算:将微电网并网时与大电网的功率交换为观测量,在此以有功功率和无功功率的偏差为目标函数,偏差最小作为等效模型与实际微电网之间的匹配关系。
目标函数定义为各时步观测功率和计算功率之间的偏差平方和 其中PM(k)、Pm(k)、Ps(k)分别为第k时刻的量测有功功率、等效电机元件有功功率和等效静态元件有功功率;QM(k)、Qm(k)、Qs(k)分别为第k时刻的量测无功功率、等效电机元件无功功率和等效静态元件无功功率。
(6)用粒子群优化算法进行参数辨识。直到满足收敛条件,得到参数的辨识结果。参数辨识过程如图4所示。
本发明按照优选实施例进行了说明,应当理解,但上述实施例不以任何形式限定本发明,凡采用等同替换或等效变换的形式所获得的技术方案,均落在本发明的保护范围之内。
Claims (1)
1.一种微电网并网等效模型的参数辨识方法,其特征在于,包括以下步骤:
1)建立等效电机动态元件模型:
ud=rsid-X′qiq+E′d
(2)
uq=rsiq+X′did+E′q
Pm=udid+uqiq
(3)
Qm=udiq-uqid
建立等效静态元件模型:
Pm、Qm为等效电机元件功率,ud、uq接入点电压d、q轴分量,id、iq接入点电流d、q轴分量,rs为等效电机定子电阻,E′d,E′q为d、q轴暂态电势,Xd,Xq为d、q轴稳态电抗,X′d,X′q为d、q轴暂态电抗,T′d0,T′q0为d、q轴开路时间常数,Ef、Eg为d、q轴转子电压,Tj为转子惯性时间常数,Xf为转子d轴电抗,Xad为d轴互抗,Xg为转子q轴电抗,Xaq为q轴互抗,Tm为机械转矩,Ps、Qs为静态元件功率,Ps0、Qs0为静态元件稳态功率,U微电网接入点母线电压幅值,U0为微电网接入点母线电压稳态幅值,pu、qu为静态元件电压有功和无功指数;t为时间,s为滑差,ωr为转子转速;ωs为系统公共参考轴的转速,均为标么值;
上述模型参数中包括等效电机元件模型基本参数、等效静态元件模型参数及关联参数,等效电机元件基本参数包括:定子电阻rs、定子漏抗Xsl、d轴互感抗Xad、q轴互感抗Xaq、转子d轴电阻rf、转子d轴漏抗Xfl,转子q轴电阻rg、转子q轴漏抗Xgl、转子惯性时间常数Tj;等效静态元件参数包括:静态元件电压有功指数pu、静态元件电压无功指数qu;关联参数包括:s0和Kmp,s0表示等效电机的初始滑差,Kmp表示为等效电机元件初始功率比例;
等效电机模型实用参数和基本参数间存在如下对应关系:
Xf=Xfl+Xad
Xg=Xgl+Xaq
Xd=Xsl+Xad
Xq=Xsl+Xaq
2)输入量测数据,包括微电网接入点电压ud+juq和接入点有功功率和无功功率PM+jQM;
3)根据等效电机有功功率比例计算得到稳态时等效电机有功功率Pm=Kmp*PM。根据粒子群优化算法初始给定的模型参数计算等效电机动态方程状态变量的初始值:d、q轴暂态电势E′d0、E′q0,等效电机的初始滑差s0;
4)利用龙格库塔法求解等效电机动态元件模型方程,获得等效电机的动态输出功率值;
5)根据电压值和等效静态元件的初试功率,计算等效静态元件的输出功率值;
6)计算目标函数即微电网总体等效模型输出功率与量测功率的偏差
其中,E(θ)表示量测数据和辨识数据的功率偏差,PM(k,θ)、QM(k,θ)表示测量得到的微电网接入点有功功率和无功功率值,k表示测量的数据点数,θ表示优化变量集,P∑(k)、Q∑(k)表示模型计算得到的有功功率和无功功率值:包括等效电机和等效静态元件的功率;
7)利用粒子群算法调整模型参数进行优化计算,判断迭代次数最终获得模型参数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210210462.7A CN102801180B (zh) | 2012-06-25 | 2012-06-25 | 一种微电网并网等效模型的参数辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210210462.7A CN102801180B (zh) | 2012-06-25 | 2012-06-25 | 一种微电网并网等效模型的参数辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102801180A true CN102801180A (zh) | 2012-11-28 |
CN102801180B CN102801180B (zh) | 2014-11-26 |
Family
ID=47200195
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210210462.7A Active CN102801180B (zh) | 2012-06-25 | 2012-06-25 | 一种微电网并网等效模型的参数辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102801180B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103023027A (zh) * | 2012-12-06 | 2013-04-03 | 湖南大学 | Dfig风电场外特性等效模型的构建方法 |
CN103529698A (zh) * | 2013-10-17 | 2014-01-22 | 广东电网公司电力科学研究院 | 发电机调速系统参数辨识方法 |
CN103631991A (zh) * | 2013-11-05 | 2014-03-12 | 国家电网公司 | 原动机调速系统参数辨识系统及方法 |
CN104537196A (zh) * | 2014-09-25 | 2015-04-22 | 西安交通大学 | 一种获取基于vsc的多端直流系统可靠性的方法 |
CN106777636A (zh) * | 2016-12-07 | 2017-05-31 | 天津大学 | 面向微电网电磁暂态实时仿真的分块并行方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011087372A (ja) * | 2009-10-14 | 2011-04-28 | Nec Engineering Ltd | 太陽光発電システムの発電量予測方法及び予測装置 |
CN102403731A (zh) * | 2011-11-01 | 2012-04-04 | 湖南大学 | 微型燃气轮机发电系统的模拟方法 |
CN102412596A (zh) * | 2011-11-28 | 2012-04-11 | 湖南大学 | 一种并网光伏发电系统的等效模拟方法 |
-
2012
- 2012-06-25 CN CN201210210462.7A patent/CN102801180B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011087372A (ja) * | 2009-10-14 | 2011-04-28 | Nec Engineering Ltd | 太陽光発電システムの発電量予測方法及び予測装置 |
CN102403731A (zh) * | 2011-11-01 | 2012-04-04 | 湖南大学 | 微型燃气轮机发电系统的模拟方法 |
CN102412596A (zh) * | 2011-11-28 | 2012-04-11 | 湖南大学 | 一种并网光伏发电系统的等效模拟方法 |
Non-Patent Citations (3)
Title |
---|
孙黎霞等: "基于Park模型的同步发电机参数辨识", 《中国电机工程学报》 * |
程颖等: "负荷模型参数辨识的粒子群优化法及其与基因算法比较", 《电力系统自动化》 * |
鞠平等: "基于物理背景的微电网总体模型", 《电力自动化设备》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103023027A (zh) * | 2012-12-06 | 2013-04-03 | 湖南大学 | Dfig风电场外特性等效模型的构建方法 |
CN103529698A (zh) * | 2013-10-17 | 2014-01-22 | 广东电网公司电力科学研究院 | 发电机调速系统参数辨识方法 |
CN103529698B (zh) * | 2013-10-17 | 2016-09-14 | 广东电网公司电力科学研究院 | 发电机调速系统参数辨识方法 |
CN103631991A (zh) * | 2013-11-05 | 2014-03-12 | 国家电网公司 | 原动机调速系统参数辨识系统及方法 |
CN104537196A (zh) * | 2014-09-25 | 2015-04-22 | 西安交通大学 | 一种获取基于vsc的多端直流系统可靠性的方法 |
CN104537196B (zh) * | 2014-09-25 | 2018-01-19 | 西安交通大学 | 一种获取基于vsc的多端直流系统可靠性的方法 |
CN106777636A (zh) * | 2016-12-07 | 2017-05-31 | 天津大学 | 面向微电网电磁暂态实时仿真的分块并行方法 |
CN106777636B (zh) * | 2016-12-07 | 2020-08-07 | 天津大学 | 面向微电网电磁暂态实时仿真的分块并行方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102801180B (zh) | 2014-11-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103887815B (zh) | 基于运行数据的风电场参数辨识及动态等值方法 | |
CN103973203B (zh) | 一种适用于安全稳定分析的大型光伏电站在线等值建模方法 | |
CN102801180B (zh) | 一种微电网并网等效模型的参数辨识方法 | |
CN107679769A (zh) | 含风电的电力系统频率响应模型建立方法和频率特性指标计算方法 | |
CN103414212B (zh) | 考虑双馈式风电机接入的配网系统的模拟方法 | |
CN102856899B (zh) | 一种微电网降低网损的方法 | |
CN102842917B (zh) | 一种通用的并网式光伏发电系统机电暂态模型 | |
CN103248049A (zh) | 含dfig风电场电力系统电压稳定域切平面的获取方法 | |
CN103986189A (zh) | 一种储能型双馈风电场简化模型建模方法 | |
CN103208817B (zh) | 一种基于二阶滑模的dfig控制方法 | |
CN105680473A (zh) | 面向光伏发电系统的通用机电暂态信息物理融合建模方法 | |
CN104617578A (zh) | 一种含风电场电力系统的可用输电能力的获取方法 | |
CN113378347B (zh) | 一种基于模块化多端口的风电机组频域阻抗建模方法 | |
CN113721476A (zh) | 100mw级可变速海水抽水蓄能机组与可再生能源联合运行系统硬件在环仿真平台及方法 | |
CN106294959A (zh) | 模型参考自适应控制与双馈风机降阶模型相结合的建模仿真方法 | |
CN203218889U (zh) | 一种通用的并网式光伏发电系统机电暂态模型 | |
CN103956767B (zh) | 一种考虑尾流效应的风电场并网稳定性分析方法 | |
CN103199524B (zh) | 一种适应多种分布式电源接入的潮流计算方法 | |
CN110336299B (zh) | 一种考虑综合能源系统小干扰稳定的配网重构方法 | |
CN102055197B (zh) | 一种建立可控串补线性化模型的方法 | |
Yu et al. | Small-signal modeling of wind farm with direct-drive PMSG using the component connection method | |
CN103532167A (zh) | 基于fpga的并网光伏系统小干扰稳定性判断系统及方法 | |
CN103248011B (zh) | 机端对称性故障双馈发电机阻抗电压源等效电路设计方法 | |
CN105140957A (zh) | 基于风电场和光伏电站聚合模型的机电振荡模式估算方法 | |
CN102707161B (zh) | 一种基于动态转移阻抗的短路容量在线监测方法及其装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |