CN112531715B - 基于虚拟电阻的下垂控制多端直流微电网潮流计算方法 - Google Patents

基于虚拟电阻的下垂控制多端直流微电网潮流计算方法 Download PDF

Info

Publication number
CN112531715B
CN112531715B CN202011293340.XA CN202011293340A CN112531715B CN 112531715 B CN112531715 B CN 112531715B CN 202011293340 A CN202011293340 A CN 202011293340A CN 112531715 B CN112531715 B CN 112531715B
Authority
CN
China
Prior art keywords
converter valve
node
vsc
power
load
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.)
Active
Application number
CN202011293340.XA
Other languages
English (en)
Other versions
CN112531715A (zh
Inventor
王守相
刘琪
赵倩宇
王振宇
田云峰
沈宇
董文琦
马步云
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Xinyuan Zhangjiakou Scenery Storage Demonstration Power Plant Co ltd
Tianjin University
Original Assignee
State Grid Xinyuan Zhangjiakou Scenery Storage Demonstration Power Plant Co ltd
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Xinyuan Zhangjiakou Scenery Storage Demonstration Power Plant Co ltd, Tianjin University filed Critical State Grid Xinyuan Zhangjiakou Scenery Storage Demonstration Power Plant Co ltd
Priority to CN202011293340.XA priority Critical patent/CN112531715B/zh
Publication of CN112531715A publication Critical patent/CN112531715A/zh
Application granted granted Critical
Publication of CN112531715B publication Critical patent/CN112531715B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/60Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及基于虚拟电阻的下垂控制多端直流微电网潮流计算方法,包括如下步骤:获取直流微电网拓扑结构、换流阀设备参数、线路参数、负荷参数及分布式电源参数;根据换流阀下垂控制参数构建虚拟电阻模型;得到直流微电网潮流计算的雅可比矩阵;采用牛顿‑拉夫逊迭代法计算节点电压,当前后两次迭代计算的节点电压差满足收敛条件时,计算完成,输出潮流计算结果。本发明通过在直流潮流模型中引入虚拟电阻,将下垂控制模型转化为精确的潮流模型,避免了直流潮流计算中松弛母线电压的更新,提高了计算速度。

Description

基于虚拟电阻的下垂控制多端直流微电网潮流计算方法
技术领域
本发明属于微电网计算技术领域,尤其涉及基于虚拟电阻的下垂控 制多端直流微电网潮流计算方法。
背景技术
随着直流型分布式电源的大规模接入和电动汽车、数据中心等直流 负荷的不断增加,配电系统的直流特征日益突出。发展直流微电网能有 效提高对光伏等直流型分布式电源的接纳能力,减少换流环节,降低线 路损耗,提高能源利用效率。
目前,潮流计算是电力系统的基本工具,前推回代法和牛顿迭代法 是配电网最常用的潮流算法:前推回代法计算速度快,适用于辐射状配 电网,虽然牛顿迭代法比前推回代法需要更多的执行时间来收敛,但是 牛顿迭代法的有效性不受网络拓扑结构的影响。因此,基于牛顿迭代法 的潮流算法更适合于直流配电网的潮流计算。
在直流微电网中,主从控制和下垂控制是两种常见的控制方式。在 下垂控制的直流微电网中,交/直流变换器的输出电压与交/直流变换器 的输出功率有关。因此,在直流微电网中没有松弛节点维持母线节点电 压恒定。
在已有的基于下垂控制的直流微电网潮流计算方法中,通常首先假 设一个最大容量的交/直流变换器的输出电压是恒定的。然后,根据每次 迭代的结果,得到该交/直流变换器的输出功率,并更新相应的输出电压。 经过多次迭代,满足了交/直流变换器输出电压的收敛条件,完成潮流计 算。在上述算法中,潮流计算多次迭代增加了计算负担,因此,基于这 些问题,提供一种更为快速准确的直流微电网潮流计算方法,具有重要 的现实意义。
发明内容
本发明的目的在于克服现有技术的不足,引入虚拟电阻模型来表示 直流微电网下垂控制对潮流的影响,通过向下垂节点加入注入功率来补 偿虚拟电阻的功率损耗,保证了虚拟电阻模型的精度,在此模型基础上 提出了基于牛顿-拉夫逊迭代法的潮流计算方法。
本发明解决其技术问题是采取以下技术方案实现的:
基于虚拟电阻的下垂控制多端直流微电网潮流计算方法,包括如下 步骤:
获取直流微电网拓扑结构、换流阀设备参数、线路参数、负荷参数 及分布式电源参数;
根据换流阀下垂控制参数构建虚拟电阻模型;
得到直流微电网潮流计算的雅可比矩阵;
采用牛顿-拉夫逊迭代法计算节点电压,当前后两次迭代计算的节点 电压差满足收敛条件时,计算完成,输出潮流计算结果。
进一步的,根据换流阀下垂控制参数构建虚拟电阻模型的具体方法 为:
在下垂控制的直流微电网中,第n个换流阀的直流侧电压UVSC,n和输出 功率PVSC,n满足:
式中,UVSC,n,ref和PVSC,n,ref分别是第n个换流阀的运行参考电压和运行参考 功率;kVSC,n是第n个换流阀的下垂系数;
在所构造的虚拟电阻模型中,换流阀下垂控制模型相当于一个理想 的电压源和一条线路,下垂节点相当于负荷节点;其中,理想电压源的 输出电压即为换流阀的运行参考电压,同时,通过向下垂节点增加负的 负载来补偿虚拟电阻的功率损耗,从而保证下垂节点电压和虚拟线路的 传输功率满足下垂控制方程;
对于第n个换流阀,抵消等效线路损耗的等效荷载Ploss,n可表示为:
式中,RVSC,n为第n个换流阀的等效线路电阻;UVSC,n,sou为第n个换流阀 的理想电压源的电压;
第n个换流阀下垂节点视为负荷节点,负荷大小为:
Pload,n=-Ploss,n-PVSC,n,ref (3)
第n个换流阀的虚拟电阻的潮流模型为:
其中,Pn是第n个换流阀的虚拟线路的传输功率。
根据(1),当换流阀的输出功率变化ΔP时,换流阀的输出电压变 化量ΔUVSC,n可以表示为:
ΔUVSC,n=kVSC,nΔP (5)
由(4)可知,第n个换流阀由于换流阀输出功率的变化而产生的输 出电压变化可以表示为:
因此,第n个换流阀下垂控制对应的虚拟电阻大小为:
RVSC,n=-kVSC,nUVSC,n,sou (7) 。
进一步的,令负荷节点序号为1,2,…,h,下垂节点的序号为h+1, h+2,…,h+m,通过虚拟电阻的应用,在直流配电网中引入了m个理 想电压源,分布式电源输出功率假定为负的负载,第i条负载母线的节 点功率Pi可表示为:
Pi=PDG,i-PL,i i=1,2,L,h (8)
其中,PDG,i和PL,i是位于节点i的分布式电源输出功率和负载功率;
设第n个换流阀的下垂节点和理想电压源节点分别为节点i和节点 k,下垂节点的节点功率可以表示为:
Pi=PVSC,n,ref+gVSC,n(UVSC,n,ref-Un)2 i=h+1,h+2,L,h+m (9)
式中,gVSC,n是第n个换流阀的虚拟线路的导纳;
直流微电网模型为:
在每次迭代中,负荷节点和下垂节点的节点功率变化可以表示为:
牛顿-拉夫逊迭代方法的迭代形式为:
其中,ΔF=[ΔP1,ΔP2,L,ΔPi],x=[U1,U2,L,Ui],k是迭代次数,J是雅可比矩阵;
根据(11),雅可比矩阵中的元素可以表示为:
进一步的,采用牛顿-拉夫逊迭代法计算节点电压,当前后两次迭 代计算的节点电压差小于10-8p.u.时,计算完成,输出潮流计算结果。
本发明的优点和积极效果是:
本发明通过在直流潮流模型中引入虚拟电阻,将下垂控制模型转化 为精确的潮流模型,避免了直流潮流计算中松弛母线电压的更新,提高 了计算速度。
具体实施方式
首先,需要说明的是,以下将以示例方式来具体说明本发明的具体 结构、特点和优点等,然而所有的描述仅是用来进行说明的,而不应将 其理解为对本发明形成任何限制。此外,在本文所提及各实施例中予以 描述或隐含的任意单个技术特征,仍然可在这些技术特征(或其等同物) 之间继续进行任意组合或删减,从而获得可能未在本文中直接提及的本 发明的更多其他实施例。
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中 的特征可以相互组合。
本实施例提供的基于虚拟电阻的下垂控制多端直流微电网潮流计算 方法,包括如下步骤:
获取直流微电网拓扑结构、换流阀设备参数、线路参数、负荷参数 及分布式电源参数;
根据换流阀下垂控制参数构建虚拟电阻模型;
得到直流微电网潮流计算的雅可比矩阵;
采用牛顿-拉夫逊迭代法计算节点电压,当前后两次迭代计算的节点 电压差满足收敛条件时,计算完成,输出潮流计算结果。
具体的,在本实施例中,根据换流阀下垂控制参数构建虚拟电阻模 型的具体方法为:
在下垂控制的直流微电网中,第n个换流阀的直流侧电压UVSC,n和输出 功率PVSC,n满足:
式中,UVSC,n,ref和PVSC,n,ref分别是第n个换流阀的运行参考电压和运行参考 功率;kVSC,n是第n个换流阀的下垂系数;
在虚拟电阻模型中,换流阀下垂控制模型相当于一个理想的电压源 和一条线路。其中,理想电压源的输出电压即为换流阀的运行参考电压。
通过向下垂节点增加负的负载来补偿虚拟电阻的功率损耗,从而保 证下垂节点电压和虚拟线路的传输功率满足下垂控制方程。对于第n个 换流阀,抵消等效线路损耗的等效荷载Ploss,n可表示为:
式中,RVSC,n为第n个换流阀的等效线路电阻;UVSC,n,sou为第n个换流阀 的理想电压源的电压。
第n个换流阀下垂节点视为负荷节点,负荷大小为:
Pload,n=-Ploss,n-PVSC,n,ref (3)
第n个换流阀的虚拟电阻的潮流模型为:
其中,Pn是第n个换流阀的虚拟线路的传输功率。
根据(1),当换流阀的输出功率变化ΔP时,换流阀的输出电压变 化量ΔUVSC,n可以表示为:
ΔUVSC,n=kVSC,nΔP (5)
由(4)可知,第n个换流阀由于换流阀输出功率的变化而产生的输 出电压变化可以表示为:
因此,第n个换流阀下垂控制对应的虚拟电阻大小为:
RVSC,n=-kVSC,nUVSC,n,sou (7)
令负荷节点序号为1,2,…,h,下垂节点的序号为h+1,h+2,…, h+m,通过虚拟电阻的应用,在直流配电网中引入了m个理想电压源, 分布式电源输出功率假定为负的负载,第i条负载母线的节点功率Pi可 表示为:
Pi=PDG,i-PL,i i=1,2,L,h (8)
其中,PDG,i和PL,i是位于节点i的分布式电源输出功率和负载功率;
设第n个换流阀的下垂节点和理想电压源节点分别为节点i和节点 k,下垂节点的节点功率可以表示为:
Pi=PVSC,n,ref+gVSC,n(UVSC,n,ref-Un)2 i=h+1,h+2,L,h+m (9)
式中,gVSC,n是第n个换流阀的虚拟线路的导纳;
直流微电网模型为:
在每次迭代中,负荷节点和下垂节点的节点功率变化可以表示为:
牛顿-拉夫逊迭代方法的迭代形式为:
其中,ΔF=[ΔP1,ΔP2,L,ΔPi],x=[U1,U2,L,Ui],k是迭代次数,J是雅可比矩阵;
根据(11),雅可比矩阵中的元素可以表示为:
采用牛顿-拉夫逊迭代法计算节点电压:当前后两次迭代计算的节点 电压差小于10-8p.u.时,计算完成,输出潮流计算结果。
相比于传统的针对下垂控制直流配电网的潮流计算,所提出的潮流 计算方法不需要多次潮流迭代修正下垂节点的节点电压,针对多个直流 微电网的仿真对比表明,本发明所提的针对下垂控制的直流微电网的潮 流计算方法的计算速度是已有的潮流计算方法的3倍。
以上实施例对本发明进行了详细说明,但所述内容仅为本发明的较 佳实施例,不能被认为用于限定本发明的实施范围。凡依本发明申请范 围所作的均等变化与改进等,均应仍归属于本发明的专利涵盖范围之 内。

Claims (2)

1.基于虚拟电阻的下垂控制多端直流微电网潮流计算方法,其特征在于,包括如下步骤:
获取直流微电网拓扑结构、换流阀设备参数、线路参数、负荷参数及分布式电源参数;
根据换流阀下垂控制参数构建虚拟电阻模型;
得到直流微电网潮流计算的雅可比矩阵;
采用牛顿-拉夫逊迭代法计算节点电压,当前后两次迭代计算的节点电压差满足收敛条件时,计算完成,输出潮流计算结果;
根据换流阀下垂控制参数构建虚拟电阻模型的具体方法为:
在下垂控制的直流微电网中,第n个换流阀的直流侧电压UVSC,n和输出功率PVSC,n满足:
式中,UVSC,n,ref和PVSC,n,ref分别是第n个换流阀的运行参考电压和运行参考功率;kVSC,n是第n个换流阀的下垂系数;
在所构造的虚拟电阻模型中,换流阀下垂控制模型相当于一个理想的电压源和一条线路,下垂节点相当于负荷节点;其中,理想电压源的输出电压即为换流阀的运行参考电压,同时,通过向下垂节点增加负的负载来补偿虚拟电阻的功率损耗,从而保证下垂节点电压和虚拟线路的传输功率满足下垂控制方程;
对于第n个换流阀,抵消等效线路损耗的等效荷载Ploss,n可表示为:
式中,RVSC,n为第n个换流阀的等效线路电阻;UVSC,n,sou为第n个换流阀的理想电压源的电压;
第n个换流阀下垂节点视为负荷节点,负荷大小为:
Pload,n=-Ploss,n-PVSC,n,ref (3)
第n个换流阀的虚拟电阻的潮流模型为:
其中,Pn是第n个换流阀的虚拟线路的传输功率;
根据(1),当换流阀的输出功率变化ΔP时,换流阀的输出电压变化量ΔUVSC,n可以表示为:
ΔUVSC,n=kVSC,nΔP (5)
由(4)可知,第n个换流阀由于换流阀输出功率的变化而产生的输出电压变化可以表示为:
因此,第n个换流阀下垂控制对应的虚拟电阻大小为:
RVSC,n=-kVSC,nUVSC,n,sou (7);
令负荷节点序号为1,2,…,h,下垂节点的序号为h+1,h+2,…,h+m,通过虚拟电阻的应用,在直流配电网中引入了m个理想电压源,分布式电源输出功率假定为负的负载,第i条负载母线的节点功率Pi可表示为:
Pi=PDG,i-PL,ii=1,2,…,h (8)
其中,PDG,i和PL,i是位于节点i的分布式电源输出功率和负载功率;
设第n个换流阀的下垂节点和理想电压源节点分别为节点i和节点k,下垂节点的节点功率可以表示为:
Pi=PVSC,n,ref+gVSC,n(UVSC,n,ref-Un)2 i=h+1,h+2,…,h+m (9)
式中,gVSC,n是第n个换流阀的虚拟线路的导纳;
直流微电网模型为:
在每次迭代中,负荷节点和下垂节点的节点功率变化可以表示为:
牛顿-拉夫逊迭代方法的迭代形式为:
其中,ΔF=[ΔP1,ΔP2,…,ΔPi],x=[U1,U2,…,Ui],k是迭代次数,J是雅可比矩阵;
根据(11),雅可比矩阵中的元素可以表示为:
2.根据权利要求1所述的基于虚拟电阻的下垂控制多端直流微电网潮流计算方法,其特征在于:采用牛顿-拉夫逊迭代法计算节点电压,当前后两次迭代计算的节点电压差小于10-8p.u.时,计算完成,输出潮流计算结果。
CN202011293340.XA 2020-11-18 2020-11-18 基于虚拟电阻的下垂控制多端直流微电网潮流计算方法 Active CN112531715B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011293340.XA CN112531715B (zh) 2020-11-18 2020-11-18 基于虚拟电阻的下垂控制多端直流微电网潮流计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011293340.XA CN112531715B (zh) 2020-11-18 2020-11-18 基于虚拟电阻的下垂控制多端直流微电网潮流计算方法

Publications (2)

Publication Number Publication Date
CN112531715A CN112531715A (zh) 2021-03-19
CN112531715B true CN112531715B (zh) 2023-07-25

Family

ID=74981498

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011293340.XA Active CN112531715B (zh) 2020-11-18 2020-11-18 基于虚拟电阻的下垂控制多端直流微电网潮流计算方法

Country Status (1)

Country Link
CN (1) CN112531715B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112952906B (zh) * 2021-03-25 2023-06-30 天津大学 下垂控制直流微电网参考运行功率快速获取方法及装置
CN113964842B (zh) * 2021-11-22 2024-04-19 上海交通大学 直流微电网潮流计算及其可靠收敛域求取方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105514971A (zh) * 2015-11-02 2016-04-20 湘潭大学 一种适用于各种运行模式微电网的潮流计算方法
CN109522676A (zh) * 2018-12-24 2019-03-26 国网辽宁省电力有限公司电力科学研究院 一种基于牛顿拉夫逊算法的分布式光伏潮流计算方法
CN109638839A (zh) * 2019-01-21 2019-04-16 东南大学 一种双极柔性直流输电系统潮流计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105514971A (zh) * 2015-11-02 2016-04-20 湘潭大学 一种适用于各种运行模式微电网的潮流计算方法
CN109522676A (zh) * 2018-12-24 2019-03-26 国网辽宁省电力有限公司电力科学研究院 一种基于牛顿拉夫逊算法的分布式光伏潮流计算方法
CN109638839A (zh) * 2019-01-21 2019-04-16 东南大学 一种双极柔性直流输电系统潮流计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Power Flow Analysis for Low-Voltage AC and DC Microgrids Considering Droop Control and Virtual Impedance;Chendan Li等;《IEEE TRANSACTIONS ON SMART GRID》;第8卷(第6期);2754-2764 *
Power sensitivity models with wide adaptability in active distribution networks considering loops and DC networks;Qi Liu等;《WILEY》;1-24 *

Also Published As

Publication number Publication date
CN112531715A (zh) 2021-03-19

Similar Documents

Publication Publication Date Title
JP6025575B2 (ja) 三相電力潮流解析のためのシステムおよび方法
Pan et al. An impedance-based stability assessment methodology for DC distribution power system with multivoltage levels
CN112531715B (zh) 基于虚拟电阻的下垂控制多端直流微电网潮流计算方法
CN106549396B (zh) 一种配电网多目标概率无功优化方法
CN110265991B (zh) 一种直流微电网的分布式协调控制方法
CN107591807B (zh) 一种新能源接入下的输电网规划的优化方法
CN109638839B (zh) 一种双极柔性直流输电系统潮流计算方法
CN108429294B (zh) 一种含能量路由器的交直流网络潮流模型及求解方法
CN110323979B (zh) 一种考虑电压稳定的发电机励磁调差系数优化整定方法
CN109774525B (zh) 基于交替方向乘子法的电动汽车优化调度方法
CN106786606B (zh) 一种基于多种随机变量的电力系统概率潮流的计算方法
CN112467748A (zh) 三相不平衡主动配电网双时标分布式电压控制方法及系统
Kim et al. Power flow calculation method of DC distribution network for actual power system
Tagore et al. Harmonic load flow analysis of radial distribution system in presence of distributed generation
CN118333219A (zh) 一种考虑风光储充接入的双层协同优化配置方法及装置
Liu et al. Function controller design in tertiary level for DC microgrid clusters
Deshmukh et al. A control algorithm for energy management and transient mitigation in DC microgrid
CN114725920B (zh) 一种适用于多分布式电源虚拟惯量直流并网系统的功率精确分配方法及装置
Javid et al. Efficient graph theory based load flow solver for DC distribution networks considering DC/DC converter models
CN114172149A (zh) 一种分布式发电集群自治控制方法及系统
CN109449968B (zh) 一种电力电子变压器与交直流源网荷多变流设备集成方法
CN109494747B (zh) 一种基于交替梯度算法的电网概率潮流计算方法
CN111313430A (zh) 一种基于局部灵敏度信息分布式电压控制方法
Miao et al. Recent advances in distributed cooperative droop control of DC microgrids: A brief survey
Kumar et al. A load flow algorithm for radial systems having renewable energy generation sources

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