CN108183487A - 一种基于线性数学模型的配电网潮流快速分析方法 - Google Patents
一种基于线性数学模型的配电网潮流快速分析方法 Download PDFInfo
- Publication number
- CN108183487A CN108183487A CN201810061899.6A CN201810061899A CN108183487A CN 108183487 A CN108183487 A CN 108183487A CN 201810061899 A CN201810061899 A CN 201810061899A CN 108183487 A CN108183487 A CN 108183487A
- Authority
- CN
- China
- Prior art keywords
- node
- power
- power flow
- node voltage
- line
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 25
- 238000013178 mathematical model Methods 0.000 title description 2
- 238000004364 calculation method Methods 0.000 claims abstract description 28
- 238000000205 computational method Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 28
- 239000013598 vector Substances 0.000 claims description 18
- 238000012545 processing Methods 0.000 claims description 8
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 230000005611 electricity Effects 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000005755 formation reaction Methods 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 239000004744 fabric Substances 0.000 claims 1
- 238000000034 method Methods 0.000 description 11
- 238000006467 substitution reaction Methods 0.000 description 4
- 230000006978 adaptation Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 206010068052 Mosaicism Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 210000003765 sex chromosome Anatomy 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- 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)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种配电网潮流快速分析方法。本发明包括节点相角计算方法、节点电压偏移量计算方法，线路潮流计算三部分，其中：配电网节点相角计算方法实现含R/X高比值支路配电网的节点电压相角计算；节点电压偏移量计算方法实现PQ节点电压幅值的线性化计算；线路潮流计算实现配电网支路潮流的快速分析计算。与基于经典直流潮流算法的潮流快速分析方法相比，本发明的配电网潮流快速分析方法能够适应配电网线路参数阻抗比较高的特点，并可进一步计算节点电压、线路无功潮流，且算法方程保留线性化结构，不须迭代收敛，在满足工程应用精度的前提下提高计算效率。
Description
技术领域
本发明创造涉及一种配电网潮流快速分析方法，具体地说，涉及一种适应配电网线路参数R/X高比值特点、具有线性化方程结构、不须迭代的配电网潮流快速分析计算技术。
背景技术
潮流快速分析是电力系统潮流分析的重要分支，准确的潮流快速分析方法对电网实时、准实时运行具有重要意义。
目前，电网潮流快速分析方法主要分为两类，一类方法利用计算机并行处理方式，结合稀疏矩阵技术、分解方法等实现电力系统的潮流快速计算，其中主要方法为迭代法。在使用迭代法时，需要考虑算法的收敛性能，特别对于非线性方程组，需要在每一步迭代中对矩阵进行处理，实现预条件的更新，方法复杂、计算量大，存在算法收敛性问题。
由于直流潮流算法模型简单、方程完全线性化、求解速度快，因此另一类方法利用直流潮流算法在电力系统的诸多场合中进行潮流快速分析与计算，如在电网静态安全分析中快速筛选过载线路；在电网调度工作中，对机组最优组合做出快速校核；在电力现货市场中为实时电价计算提供参考。
由于经典直流潮流算法对于线路参数较为敏感，在应用时，线路参数R/X的比值应尽可能的小，至少在0.25以下，而配电网线路参数R/X的比值通常较高，一般均大于1，此时应用经典直流潮流算法带来的误差不可忽略。且经典直流潮流算法将节点电压幅值均视为额定电压幅值，忽略了节点电压偏移对线路潮流的影响。目前文献对经典直流潮流算法的改进或利用等值网损负荷模型，计及电阻参数引起的有功网损，但需要反复迭代，计算效率较低；或基于初始潮流断面，需要启动条件；或虽考虑节点电压幅值计算，但在节点电压幅值计算中，沿用经典直流潮流算法的节点相角计算值，未考虑在配电网中应用时线路参数引入的较大误差。为此，本发明提出一种适应配电网线路参数R/X高比值特点、计及节点电压偏移、具有线性化方程结构、不须迭代收敛的配电网潮流快速分析方法。
发明内容
本发明的目的在于提供一种配电网潮流快速分析方法。
本发明的技术方案为：一种配电网潮流快速分析方法包括节点相角计算方法、节点电压偏移量计算方法、线路潮流计算。利用经典直流潮流算法计算各条线路的有功潮流初始分布；结合线路有功潮流初始分布结果，利用节点相角计算方法计算配电网各节点的节点相角值；结合节点相角计算结果，利用节点电压偏移量计算方法计算PQ节点的节点电压偏移量，并修正各节点的电压幅值；结合节点相角计算结果、节点电压偏移量计算结果，利用无功潮流方程计算线路无功潮流，利用计及节点电压和无功潮流耦合的直流潮流方程更新线路有功潮流分布和节点相角计算值。
其中，节点电压相角计算方法包括如下步骤：
S1：对配电网节点进行编号，获取配电网的网络参数，利用经典直流潮流算法计算配电网各条线路的有功潮流初步分布；
S2：构建有功潮流耦合向量Pcp：
式中：Pcp为N-1维列向量，不包括参考节点所在列；xij为节点i、j之间线路的线路电抗参数；rij为节点i、j之间线路的线路电阻参数；P0(i,j)为节点i、j之间线路的有功潮流初步分布；
S3：利用线路电阻参数构成支路电导，建立节点电导矩阵R0：
式中：R0为N-1阶方阵，不包括参考节点所在的行和列；R0(i,j)为节点电导矩阵R0第i行第j列的元素；
S4：构建计算节点电压相角θ向量的线性矩阵方程：
QSP-Pcp＝R0θ；
式中：QSP为节点净注入无功列向量，各元素为相应节点的节点净注入无功功率；θ为节点电压相角列向量，各元素为相应节点的节点电压相角计算结果。
节点电压偏移量计算方法包括：
S1：计算所述配电网的N-1阶节点导纳矩阵，根据节点导纳矩阵的实部和虚部构造中间变量：
ωi＝Qicosθ(i)-Pisinθ(i)；
ηij＝Bijcosθ(j)+Gijsinθ(j)；
式中：Gij、Bij分别为节点导纳矩阵第i行第j列元素Y0(i,j)的实部与虚部；θ(i)为节点i的节点电压相角；
S2：取节点潮流方程的虚部进行恒等变换，并进行近似处理1/(1+ΔVi)≈1-ΔVi，得到关于PQ节点i的节点电压偏移量的线性化方程：
式中：ΔVi为节点i的节点电压与额定电压之间的节点偏移量；Vj为节点j的节点电压；PQ为所述配电网的PQ节点集；PV为所述配电网的PV节点集；
S3:利用节点电压幅值偏移量修正各节点的电压幅值：
Vi＝1+ΔVi。
线路潮流计算部分，包括：
S1.根据上述节点电压相角计算结果和节点电压幅值计算结果，利用线路交流潮流方程计算各条线路的无功潮流：
Qij＝-Vi Vj(Gijsinθij-Bijcosθij)-Vi 2Bij；
式中：Qij为节点i与节点j之间线路的线路无功潮流；Vi为节点i的节点电压幅值；θij为节点i的节点电压相角与节点j的节点电压相角之间的相角差；
由上述线路无功潮流Qij，构建线路无功潮流耦合向量Qcp：
S2.在经典直流潮流方程的基础上，计及线路无功潮流、修正后的节点电压等，得到线路有功线性潮流方程：
式中：Pij为节点i与节点j之间线路的线路有功潮流；
在经典直流潮流算法的节点导纳矩阵的基础上，建立计及节点电压幅值的的节点导纳矩阵B：
S3.由上述线路有功线性潮流方程和计及节点电压幅值的节点导纳矩阵构建矩阵线性方程组：
PSP-Qcp＝Bθ；
式中：PSP为节点净注入有功功率列向量，各元素为相应节点的节点净注入有功功率；
求解矩阵线性方程组，更新节点相角值，并根据线路有功线性潮流方程求解配电网的线路有功潮流。
本发明的有益效果为：
(1)本发明所提一种配电网潮流快速分析方法以1/rij为支路电导，构建节点电导矩阵，计及了线路电阻参数对线路潮流的影响，适用于线路参数R/X比值较高的网络。因此将基于直流潮流算法的潮流快速分析方法的适用范围从线路参数R/X比值小于0.25的高压输电网拓展到线路参数R/X比值大于1的中低压配电网；
(2)本发明通过对节点潮流方程进行恒等变换和近似处理，得到关于PQ节点电压偏移量的线性方程，由于仅做一次1/(1+ΔVi)≈1-ΔVi的近似处理，且正常运行的网络，其节点电压偏移量不会超过额定电压幅值的-10％～10％，因此所做近似处理引起的误差在-0.901％～0.9901％之间，引入的误差较小。并且，关于节点电压偏移量的线性方程的计算使用本发明所提节点相角计算方法的计算结果，其计算精度比使用经典直流潮流算法的计算结果时更高；
(3)本发明所提一种配电网潮流快速分析方法保留了直流潮流算法的线性方程结构，不存在迭代收敛问题，计算量较小，计算速度更快，更加适用于配电网的实时和准实时运行。
附图说明
图1：本发明实验所用的IEEE33节点配电网网络结构示意图；
图2：本发明实验所用的IEEE33节点配电网络线路参数与节点功率参数；
图3：本发明所提配电网潮流快速分析方法与经典前推回代潮流算法的潮流计算结果对比；
图4：本发明所提配电网潮流快速分析方法与经典前推回代潮流算法的节点电压计算结果对比；
图5：本发明所提配电网潮流快速分析方法与经典前推回代潮流算法的线路潮流计算结果对比。
具体实施方式
下面结合附图，对本发明的具体实施方式作进一步详细说明。应该强调的是，下述说明仅仅是示例性的，而不是为了限制本发明的范围及其应用。
本发明提供了一种配电网潮流快速分析方法，其包括以下步骤：
S1.对配电网进行节点编号，并选取首端节点为参考节点，获取各节点的功率参数；
S2.以节点编号小的节点作为线路的首节点，以节点编号大的节点作为线路的末节点，定义首端节点到末端节点为线路正方向，以首末节点的编号标识线路，并获取各条线路的线路电阻参数和线路电抗参数；
S3.由经典直流潮流算法计算线路有功潮流的初步分布；
S4.利用上述线路有功潮流初步分布，构建N-1维线路有功潮流耦合向量Pcp：
式中：Pcp为N-1维列向量，不包括参考节点所在列；xij为节点i、j之间线路的线路电抗参数；rij为节点i、j之间线路的线路电阻参数；P0(i,j)为节点i、j之间线路的有功潮流初步分布；
S5.由线路电阻参数构成支路电导，建立节点电导矩阵R0：
式中：R0为N-1阶方阵，不包括参考节点所在的行和列；R0(i,j)为节点导纳矩阵R0第i行第j列的元素；
S6.构建计算节点电压相角θ的线性矩阵方程组：
QSP-Pcp＝R0θ；
式中：QSP为节点净注入无功列向量，各元素为相应节点的节点净注入无功功率；θ为节点电压相角列向量，各元素为相应节点的节点电压相角计算结果。
S7.计算所述配电网的N-1阶节点导纳矩阵，根据节点导纳矩阵的实部和虚部构造中间变量：
ωi＝Qicosθ(i)-Pisinθ(i)；
ηij＝Bijcosθ(j)+Gijsinθ(j)；
式中：Gij、Bij分别为节点导纳矩阵第i行第j列元素Y0(i,j)的实部与虚部；θ(i)为节点i的节点电压相角；
S8.取节点潮流方程的虚部进行恒等变换，并进行近似处理1/(1+ΔVi)≈1-ΔVi，得到关于PQ节点i的节点电压偏移量的线性化方程：
式中：ΔVi为节点i的节点电压与额定电压之间的节点偏移量；Vj为节点j的节点电压；PQ为所述配电网的PQ节点集；PV为所述配电网的PV节点集；
S9.利用节点电压幅值偏移量修正各节点的电压幅值：
Vi＝1+ΔVi。
S10.根据上述节点电压相角计算结果和节点电压幅值计算结果，利用线路潮流方程计算各条线路的无功潮流：
Qij＝-Vi Vj(Gijsinθij-Bijcosθij)-Vi 2Bij；
由上述线路无功潮流Qij，构建线路无功潮流耦合列向量Qcp：
在经典直流潮流算法节点导纳矩阵的基础上，建立计及节点电压幅值的的节点导纳矩阵B：
S11.由上述线路有功线性潮流方程和计及节点电压幅值的节点导纳矩阵构建矩阵线性方程组：
PSP-Qcp＝Bθ；
式中：PSP为节点净注入有功功率列向量，各元素为相应节点的节点净注入有功功率；
求解矩阵线性方程组，更新节点相角值，并根据线路有功线性潮流方程求解配电网的线路有功潮流：
根据图1所示的经典IEEE33节点配电网结构图和图2所示的配电网线路参数和节点负荷参数，利用经典的前推回代算法计算IEEE33节点配电网的节点电压幅值、节点电压相角、线路潮流和线路有功网损。对比利用本发明方法得到的IEEE33节点配电网潮流分析结果，得到图3、图4、图5所示对比结果，可以验证本发明所提出的一种配电网潮流快速分析方法及其验证过程的有效性。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。
Claims (3)
1.一种配电网潮流快速分析方法，包括节点相角计算方法、节点电压偏移量计算方法、线路潮流计算三部分。
2.如权利要求1所述的节点相角计算方法，其特征在于，通过线路电阻参数构建节点电导矩阵，适应配电网线路参数R/X高比值特点，R指配电网线路的线路电阻参数，X指配电网线路的线路电抗参数，包括如下步骤：
S1：对配电网节点进行编号，配电网节点数为N，利用经典直流潮流算法计算配电网线路有功潮流的初步分布；
S2：构建有功潮流耦合向量Pcp：
式中：Pcp为N-1维列向量，不包括参考节点所在列；xij为节点i、j之间线路的线路电抗参数；rij为节点i、j之间线路的线路电阻参数；P0(i,j)为节点i、j之间线路的有功潮流初步分布；
S3：利用线路电阻参数构成支路电导，建立节点电导矩阵R0：
式中：R0为N-1阶方阵，不包括参考节点所在的行和列；R0(i,j)为节点电导矩阵R0第i行第j列的元素；
S4：构建计算节点电压相角θ向量的线性矩阵方程组：
QSP-Pcp＝R0θ；
式中：QSP为节点净注入无功向量，各元素为相应节点的节点净注入无功功率；θ为节点电压相角向量，各元素为相应节点的节点电压相角计算结果。
3.根据权利要求1所述的节点电压偏移量计算方法，其特征在于，通过对潮流方程进行恒等变换和近似处理，构建线性方程求解PQ节点的节点电压偏移，包括如下步骤：
S1：计算所述配电网的N-1阶节点导纳矩阵，根据节点导纳矩阵的实部和虚部构造中间变量：
ωi＝Qicosθ(i)-Pisinθ(i)；
ηij＝Bijcosθ(j)+Gijsinθ(j)；
式中：Gij、Bij分别为节点导纳矩阵第i行第j列元素Y0(i,j)的实部与虚部；θ(i)为节点i的节点电压相角；
S2：取节点潮流方程的虚部进行恒等变换，并进行近似处理1/(1+ΔVi)≈1-ΔVi，得到关于PQ节点i的节点电压偏移量的线性化方程：
式中：ΔVi为节点i的节点电压与额定电压之间的节点电压偏移量；Vj为节点j的节点电压幅值；PQ为所述配电网的PQ节点集；PV为所述配电网的PV节点集；
S3:利用节点电压幅值偏移量修正各节点的电压幅值：
Vi＝1+ΔVi。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810061899.6A CN108183487B (zh) | 2018-01-23 | 2018-01-23 | 一种基于线性数学模型的配电网潮流快速分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810061899.6A CN108183487B (zh) | 2018-01-23 | 2018-01-23 | 一种基于线性数学模型的配电网潮流快速分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108183487A true CN108183487A (zh) | 2018-06-19 |
| CN108183487B CN108183487B (zh) | 2022-03-15 |
Family
ID=62551199
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810061899.6A Expired - Fee Related CN108183487B (zh) | 2018-01-23 | 2018-01-23 | 一种基于线性数学模型的配电网潮流快速分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108183487B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111864741A (zh) * | 2020-07-27 | 2020-10-30 | 国网安徽省电力有限公司电力科学研究院 | 一种线路参数误差对功率分布影响的定量分析方法及系统 |
| CN114912068A (zh) * | 2022-06-01 | 2022-08-16 | 西南交通大学 | 一种dc牵引网络的潮流分析方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104573384A (zh) * | 2015-01-23 | 2015-04-29 | 东南大学 | 一种适用于配电网在线分析的改进直流潮流算法 |
| CN107230982A (zh) * | 2017-07-07 | 2017-10-03 | 广西大学 | 一种微电网线性化潮流计算方法 |
-
2018
- 2018-01-23 CN CN201810061899.6A patent/CN108183487B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104573384A (zh) * | 2015-01-23 | 2015-04-29 | 东南大学 | 一种适用于配电网在线分析的改进直流潮流算法 |
| CN107230982A (zh) * | 2017-07-07 | 2017-10-03 | 广西大学 | 一种微电网线性化潮流计算方法 |
Non-Patent Citations (1)
| Title |
|---|
| 陈醒等: "基于双解耦的配电网三相不平衡快速潮流算法", 《电力自动化设备》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111864741A (zh) * | 2020-07-27 | 2020-10-30 | 国网安徽省电力有限公司电力科学研究院 | 一种线路参数误差对功率分布影响的定量分析方法及系统 |
| CN111864741B (zh) * | 2020-07-27 | 2022-09-30 | 国网安徽省电力有限公司电力科学研究院 | 一种线路参数误差对功率分布影响的定量分析方法及系统 |
| CN114912068A (zh) * | 2022-06-01 | 2022-08-16 | 西南交通大学 | 一种dc牵引网络的潮流分析方法 |
| CN114912068B (zh) * | 2022-06-01 | 2023-03-10 | 西南交通大学 | 一种dc牵引网络的潮流分析方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108183487B (zh) | 2022-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ju et al. | An extension of FBS three-phase power flow for handling PV nodes in active distribution networks | |
| CN106532711B (zh) | 随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法 | |
| CN104037764B (zh) | 一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法 | |
| CN105656031B (zh) | 基于高斯混合分布特征的含风电电力系统安全风险评估方法 | |
| CN103904643B (zh) | 一种考虑网损的直流潮流计算方法 | |
| CN107332240A (zh) | 基于优化模型的电力系统静态电压稳定域边界搜索的方法 | |
| CN104899396B (zh) | 一种修正系数矩阵的快速分解法潮流计算方法 | |
| CN106385035B (zh) | 一种计及充电电容的直流潮流计算方法 | |
| CN104578159A (zh) | 一种含分布式电源的配电网三相潮流修正方法 | |
| CN103956735B (zh) | 一种分布式发电系统的谐波潮流分析方法 | |
| CN103632046A (zh) | 一种电网潮流计算方法 | |
| CN104600697A (zh) | 一种计及温度影响的拟直流最优潮流方法 | |
| CN108183487A (zh) | 一种基于线性数学模型的配电网潮流快速分析方法 | |
| CN106410811B (zh) | 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法 | |
| CN106856327B (zh) | 一种含小阻抗支路电网串联补偿快速分解法潮流计算方法 | |
| CN106532712A (zh) | 含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法 | |
| CN114629122A (zh) | 一种交直流混合配电网线性潮流计算方法 | |
| CN105470948A (zh) | 潮流控制器在直流潮流模型下功率注入模型及选址方法 | |
| CN107196307A (zh) | 一种变压器故障后电网有功潮流快速估算的方法 | |
| CN105811394B (zh) | 一种用于电力系统潮流近似计算的电压估算方法 | |
| CN103632052A (zh) | 一种配电网潮流计算方法 | |
| CN104092213B (zh) | 一种基于优化方法的不确定潮流支路功率分析方法 | |
| CN106529089B (zh) | 用于含小阻抗支路电网的补偿法快速分解法潮流计算方法 | |
| CN106712029B (zh) | 小阻抗支路pq端点变雅可比矩阵的牛顿法潮流计算方法 | |
| CN105977967B (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220315 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |