CN111245013B - DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling - Google Patents

DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling Download PDF

Info

Publication number
CN111245013B
CN111245013B CN202010123388.XA CN202010123388A CN111245013B CN 111245013 B CN111245013 B CN 111245013B CN 202010123388 A CN202010123388 A CN 202010123388A CN 111245013 B CN111245013 B CN 111245013B
Authority
CN
China
Prior art keywords
rsc
gsc
impedance
axis
current
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.)
Expired - Fee Related
Application number
CN202010123388.XA
Other languages
Chinese (zh)
Other versions
CN111245013A (en
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.)
North China Electric Power University
Original Assignee
North China Electric Power 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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN202010123388.XA priority Critical patent/CN111245013B/en
Publication of CN111245013A publication Critical patent/CN111245013A/en
Application granted granted Critical
Publication of CN111245013B publication Critical patent/CN111245013B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/24Arrangements for preventing or reducing oscillations of power in networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

The invention discloses a DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling, which comprises a filter for filtering rotor current and base frequency components of GSC current, a voltage instruction virtual impedance generation device for generating a voltage instruction corresponding to virtual impedance, and a suppression device for superposing the voltage instruction corresponding to virtual impedance and an original voltage given value to obtain a new voltage given value; the equivalent impedance of a doubly-fed wind power generation system is changed by adding virtual impedance control in a Rotor Side Converter (RSC) and a Grid Side Converter (GSC) of the doubly-fed wind power generator, so that subsynchronous oscillation is quickly and effectively suppressed, stability and economy are guaranteed without adding other equipment, the suppression effect is stronger compared with the traditional additional damping control, the suppression effect is realized in the whole secondary same-frequency band, and the suppression effect is improved.

Description

DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling
Technical Field
The invention relates to the technical field of wind power generation, in particular to a DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling.
Background
With the increase of wind power capacity, the phenomenon of subsynchronous oscillation of a plurality of large wind power plant grid-connected systems occurs at home and abroad, so that large-area wind power sets are separated from a power grid, and the safety of the power grid is seriously influenced. To avoid this effect, fast and effective measures need to be taken to suppress subsynchronous oscillations.
The existing suppression strategies for subsynchronous oscillation of doubly-fed wind generators (DFIGs) are mainly classified into 3 types: the flexible power transmission system and an additional damping controller, a subsynchronous filter and an additional subsynchronous damping controller are added on a fan converter control loop. The first two types of suppression strategies require additional hardware devices in the system and are not cost-effective. And the subsynchronous damping controller strategy added on the control loop of the existing fan converter only adds a virtual resistor on a Rotor Side Converter (RSC), so that the suppression effect is limited. Particularly, for a wind power generation system with high series compensation degree, the subsynchronous oscillation degree is relatively serious, and the subsynchronous oscillation cannot be inhibited only by adding virtual damping in RSC. Therefore, a new technical solution is needed to solve the above problems.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling, which improves the capability of suppressing subsynchronous oscillation, reduces the equipment cost and ensures the stability of a system.
The invention provides a DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling, which comprises a filter for filtering rotor current and a fundamental frequency component of GSC current, a voltage instruction virtual impedance generation device for generating a virtual impedance and a suppression device for obtaining a new voltage given value by superposing a voltage instruction corresponding to the virtual impedance and an original voltage given value, and comprises the following steps:
step S1: the DQ axis component i of the rotor current and GSC currentrd、irq、igd、igqRespectively inputting into a filter to obtain a current i 'containing only subsynchronous oscillation frequency components'rd、i′rq、i′gd、i′gq
Step S2: i'rd、i′rqMultiplying the gain Kr _ RSC to obtain DQ axis components Urd _ RSC and Urq _ RSC of the RSC side virtual resistor; i'rd、i′rqMultiplying the gain Kx _ RSC to obtain DQ axis components Uxd _ RSC and Uxq _ RSC of the virtual reactance on the RSC side; i'gd、i′gqMultiplying the gain Kr _ gsc to obtain DQ axis components Urd _ gsc and Urq _ gsc of the RSC side virtual resistor;
step S3: the D-axis voltage set values of Urd _ RSC, -Uxd _ RSC and the original RSC are superposed to be used as a new D-axis voltage ring set value of the RSC
Figure BDA0002393678120000021
The Q-axis voltage given values of Urq _ RSC, Uxq _ RSC and the original RSC are superposed to be used as a new Q-axis voltage ring given value of RSC
Figure BDA0002393678120000022
Superposing the specified values of the DQ axis voltages of the Urd _ GSC, the Urq _ GSC and the original GSC to be used as the new DQ axis voltage ring specified value of the GSC
Figure BDA0002393678120000023
The invention has the beneficial effects that: the equivalent impedance of the doubly-fed wind power generation system is changed by adding virtual impedance control in the RSC and grid-side converter (GSC) of the doubly-fed wind power generator, so that subsynchronous oscillation is quickly and effectively suppressed, other equipment is not needed to be added, the stability and the economy of the doubly-fed wind power generation system are guaranteed, the suppression effect is stronger compared with the traditional additional damping control, the suppression effect is realized in the whole secondary same frequency band, and the suppression effect is improved.
Drawings
FIG. 1 is a schematic structural diagram of a double-fed wind power generation system according to the present invention;
FIG. 2 is a diagram of an original RSC control block;
FIG. 3 is a block diagram of the original GSC control;
FIG. 4 is an equivalent circuit diagram of the doubly-fed wind power generation system of the present invention;
FIG. 5 is an equivalent circuit diagram of multi-branch impedance remodeling of the doubly-fed wind power generation system in the invention;
fig. 6 is a control block diagram of a sub-synchronous oscillation suppression strategy of a doubly-fed wind turbine generator based on multi-branch impedance remodeling.
Detailed Description
The following examples are only for illustrating the technical solutions of the present invention more clearly, and therefore are only examples, and the protection scope of the present invention is not limited thereby.
It is to be noted that, unless otherwise specified, technical or scientific terms used herein shall have the ordinary meaning as understood by those skilled in the art to which the invention pertains.
The double-fed wind power generation system is shown in fig. 1 and comprises a double-fed fan, a GSC, a RSC, a grid-side filter, a transformer, a power transmission line, a series compensation capacitor and the like. An equivalent circuit of the impedance model of the doubly-fed wind power generation system established based on the original RSC impedance model (shown in FIG. 2), the GSC impedance model (shown in FIG. 3) and the DFIG impedance model is shown in FIG. 4, and the impedance model is equivalently aggregated into the RLC circuit by utilizing the Noton theorem. The equivalent impedance of the system is:
Figure BDA0002393678120000031
in the formula, Req、Leq、CeqRespectively, the equivalent resistance, the equivalent inductance and the equivalent capacitance of the system.
The natural oscillation frequency of the system is:
Figure BDA0002393678120000032
at natural oscillation frequency, the system equivalent reactance is 0 (i.e., the system equivalent reactance is 0)
Figure BDA0002393678120000033
). At this time, if the equivalent resistance Req<0, the system will oscillate. Therefore, the equivalent resistance of the system needs to be changed to change the equivalent resistance under the natural oscillation frequency from negative to positive, thereby achieving the purpose of inhibiting subsynchronous oscillation; the DFIG impedance model expression is
Figure BDA0002393678120000034
In the formula, RrIs rotor resistance, RsIs stator resistance, Lls、LlrRespectively, the leakage inductance of the stator and the rotor, and slip(s) is the slip ratio; GSC impedance model expression:
Figure BDA0002393678120000035
in the formula, Kgp、KgiRespectively are GSC current inner loop proportion and integral coefficient, omega is synchronous rotation speed, LfilterA network side filter inductor; RSC impedance model expression:
Figure BDA0002393678120000041
in the formula, Krp、KriRespectively, RSC inner ring proportion and integral coefficient, LrFor self-inductance of the rotor winding, LsFor self-inductance of the stator winding, LmThe stator and rotor windings are mutually inducted.
Specifically, the method comprises a filter for filtering fundamental frequency components of rotor current and GSC current, a virtual impedance generating device for generating a voltage command corresponding to virtual impedance, and a suppression device for superposing the voltage command corresponding to the virtual impedance and an original voltage given value to obtain a new voltage given value, and comprises the following steps of:
step S1: the DQ axis component i of the rotor current and GSC currentrd、irq、igd、igqSeparately fed to filters to obtain components containing subsynchronous oscillation frequency onlyCurrent i'rd、i′rq、i′gd、i′gq
Step S2: i'rd、i′rqMultiplying the gain Kr _ RSC to obtain DQ axis components Urd _ RSC and Urq _ RSC of the RSC side virtual resistor; i'rd、i′rqMultiplying the gain Kx _ RSC to obtain DQ axis components Uxd _ RSC and Uxq _ RSC of the virtual reactance on the RSC side; i'gd、i′gqMultiplying the gain Kr _ gsc to obtain DQ axis components Urd _ gsc and Urq _ gsc of the RSC side virtual resistor;
step S3: the D-axis voltage set values of Urd _ RSC, -Uxd _ RSC and the original RSC are superposed to be used as a new D-axis voltage ring set value of the RSC
Figure BDA0002393678120000042
The Q-axis voltage given values of Urq _ RSC, Uxq _ RSC and the original RSC are superposed to be used as a new Q-axis voltage ring given value of RSC
Figure BDA0002393678120000043
Superposing the specified values of the DQ axis voltages of the Urd _ GSC, the Urq _ GSC and the original GSC to be used as the new DQ axis voltage ring specified value of the GSC
Figure BDA0002393678120000044
The new voltage instruction value changes the equivalent impedance of the double-fed wind power generation system, so that the equivalent resistance under the oscillation frequency is changed from negative to positive, and subsynchronous oscillation is effectively inhibited.
The fundamental frequency filter is used for filtering the fundamental frequency components of the rotor current and the GSC current; the DQ axis component i of the rotor current and GSC currentrd、irq、igd、igqRespectively inputting the current into a filter, calculating a direct current component in the current through Fourier transform (a fundamental frequency component in an ABC coordinate system corresponds to a direct current component in a DQ coordinate system), and subtracting the direct current component from the input current to obtain a current i 'only containing subsynchronous oscillation frequency components'rd、i′rq、i′gd、i′gq
The virtual impedance generating device is used for generating a voltage command corresponding to the virtual impedance; i'rd、i′rqMultiplying the gain Kr _ RSC to obtain DQ axis components Urd _ RSC and Urq _ RSC of the RSC side virtual resistor; i'rd、i′rqMultiplying the gain Kx _ RSC to obtain DQ axis components Uxd _ RSC and Uxq _ RSC of the virtual reactance on the RSC side; i'gd、i′gqMultiplying by a gain Kr _ gsc respectively to obtain DQ axis components Urd _ gsc and Urq _ gsc of the RSC side virtual resistor;
the suppression device is used for superposing a voltage command corresponding to the virtual impedance and an original voltage given value to obtain a new voltage given value; the D-axis voltage set values of Urd _ RSC, -Uxd _ RSC and the original RSC are superposed to be used as a new D-axis voltage ring set value of the RSC
Figure BDA0002393678120000051
The Q-axis voltage given values of Urq _ RSC, Uxq _ RSC and the original RSC are superposed to be used as a new Q-axis voltage ring given value of RSC
Figure BDA0002393678120000052
Superposing the specified values of the DQ axis voltages of the Urd _ GSC, the Urq _ GSC and the original GSC to be used as the new DQ axis voltage ring specified value of the GSC
Figure BDA0002393678120000053
The new voltage instruction value changes the equivalent impedance of the double-fed wind power generation system, so that the equivalent resistance under the oscillation frequency is changed from negative to positive, and subsynchronous oscillation is effectively inhibited.
As shown in FIG. 6, the DQ-axis component i of the rotor current, GSC currentrd、irq、igd、igqThe current can be obtained by current sampling of the current transformer through Park conversion. Will ird、irq、igd、igqRespectively inputting the current into a filter, calculating a direct current component in the current through Fourier transform (a fundamental frequency component in an ABC coordinate system corresponds to a direct current component in a DQ coordinate system), and subtracting the direct current component from the input current to obtain a current i 'only containing subsynchronous oscillation frequency components'rd、i′rq、i′gd、i′gq;i′rd、i′rqMultiplying the gain Kr _ RSC to obtain DQ axis components Urd _ RSC and Urq _ RSC of the RSC side virtual resistor; i'rd、i′rqMultiplying the gain Kx _ RSC to obtain DQ axis components Uxd _ RSC and Uxq _ RSC of the virtual reactance on the RSC side; i'gd、i′gqMultiplying the gain Kr _ gsc to obtain DQ axis components Urd _ gsc and Urq _ gsc of the RSC side virtual resistor; the D-axis voltage set values of Urd _ RSC, -Uxd _ RSC and the original RSC are superposed to be used as a new D-axis voltage ring set value of the RSC
Figure BDA0002393678120000061
The Q-axis voltage given values of Urq _ RSC, Uxq _ RSC and the original RSC are superposed to be used as a new Q-axis voltage ring given value of RSC
Figure BDA0002393678120000062
Superposing the specified values of the DQ axis voltages of the Urd _ GSC, the Urq _ GSC and the original GSC to be used as the new DQ axis voltage ring specified value of the GSC
Figure BDA0002393678120000063
The new voltage command value is equivalent to adding virtual impedance Z to each of the RSC and GSC branchesrsc_ssrAnd Zgsc_ssrThereby changing the equivalent impedance of the doubly-fed wind power generation system, the schematic circuit diagram of which is shown in fig. 5. And the inhibition method utilizes a semi-physical simulation experiment to verify the effectiveness of the invention.
Simulation experiments show that when no additional control is selected, the system generates subsynchronous oscillation with oscillation frequency of 32Hz, the maximum amplitude is 20MW, and the impact on a power grid is severe; when a machine side resistance remodeling control parameter Kr _ RSC with the fastest convergence rate is added in the RSC control strategy, the maximum amplitude is 16MW, constant amplitude oscillation with the amplitude of about +/-4 MW still exists after 6.18s, and oscillation cannot be inhibited; and (3) selecting a multi-branch impedance remodeling inhibition strategy, namely when three impedance remodeling control parameters of Kr _ rsc, Kx _ rsc and Kx _ gsc are added at the same time, the maximum amplitude is 14.1MW, and after 6.2s, the stable operation state is recovered.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; such modifications and substitutions do not depart from the spirit and scope of the present invention, and they should be construed as being included in the following claims and description.

Claims (1)

1. The DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling is characterized by comprising the following steps: the method comprises a filter for filtering the rotor current and the fundamental frequency component of the GSC current, a voltage instruction virtual impedance generating device for generating a virtual impedance, and a restraining device for superposing the voltage instruction corresponding to the virtual impedance and the original voltage given value to obtain a new voltage given value, and comprises the following steps:
step S1: the DQ axis component i of the rotor current and GSC currentrd、irq、igd、igqRespectively inputting into a filter to obtain a current i 'containing only subsynchronous oscillation frequency components'rd、i′rq、i′gd、i′gq
Step S2: i'rd、i′rqMultiplying the gain Kr _ RSC to obtain DQ axis components Urd _ RSC and Urq _ RSC of the RSC side virtual resistor; i'rd、i′rqMultiplying the gain Kx _ RSC to obtain DQ axis components Uxd _ RSC and Uxq _ RSC of the virtual reactance on the RSC side; i'gd、i′gqMultiplying the gain Kr _ gsc to obtain DQ axis components Urd _ gsc and Urq _ gsc of the RSC side virtual resistor;
step S3: the D-axis voltage given values of Urd _ RSC, -Uxd _ RSC and the original RSC are superposed to be used as a new D-axis voltage ring given value of the RSC
Figure FDA0002393678110000011
The Q-axis voltage given values of Urq _ RSC, Uxq _ RSC and the original RSC are superposed to be used as a new Q-axis voltage ring given value of RSC
Figure FDA0002393678110000012
Mixing Urd _ gsc and Urq _ gsc is superposed with the DQ axis voltage given value of the original GSC to be used as a new DQ axis voltage ring given value of the GSC
Figure FDA0002393678110000013
CN202010123388.XA 2020-02-27 2020-02-27 DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling Expired - Fee Related CN111245013B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010123388.XA CN111245013B (en) 2020-02-27 2020-02-27 DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010123388.XA CN111245013B (en) 2020-02-27 2020-02-27 DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling

Publications (2)

Publication Number Publication Date
CN111245013A CN111245013A (en) 2020-06-05
CN111245013B true CN111245013B (en) 2022-05-13

Family

ID=70875046

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010123388.XA Expired - Fee Related CN111245013B (en) 2020-02-27 2020-02-27 DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling

Country Status (1)

Country Link
CN (1) CN111245013B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112039112A (en) * 2020-06-23 2020-12-04 湖南大学 Subsynchronous oscillation suppression method and control system for series compensation grid-connected system of virtual synchronous machine
CN114006387B (en) * 2021-11-19 2024-02-02 华北电力大学 Self-adaptive frequency subsynchronous oscillation suppression method and system based on multi-branch impedance
CN116544964B (en) * 2023-05-06 2023-12-01 燕山大学 Impedance optimization method of wind power generation system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107017646A (en) * 2017-05-25 2017-08-04 东南大学 The double-fed blower fan sub-synchronous oscillation suppression method controlled based on virtual impedance
CN108321841A (en) * 2018-02-26 2018-07-24 国电南瑞科技股份有限公司 The double-fed fan motor unit sub-synchronous oscillation suppression method of wideband rotor additional damping
CN108631337A (en) * 2018-05-23 2018-10-09 哈尔滨工业大学 The double-fed fan motor unit subsynchronous resonance suppressing method of damping control is cooperateed with based on converter
CN109120001A (en) * 2018-09-29 2019-01-01 华中科技大学 Double-fed fan motor field grid-connected system sub-synchronous oscillation suppression method based on virtual resistance

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109546664B (en) * 2017-09-21 2022-05-24 通用电气公司 Power generation system, system for suppressing subsynchronous oscillations, and method for controlling operation of power system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107017646A (en) * 2017-05-25 2017-08-04 东南大学 The double-fed blower fan sub-synchronous oscillation suppression method controlled based on virtual impedance
CN108321841A (en) * 2018-02-26 2018-07-24 国电南瑞科技股份有限公司 The double-fed fan motor unit sub-synchronous oscillation suppression method of wideband rotor additional damping
CN108631337A (en) * 2018-05-23 2018-10-09 哈尔滨工业大学 The double-fed fan motor unit subsynchronous resonance suppressing method of damping control is cooperateed with based on converter
CN109120001A (en) * 2018-09-29 2019-01-01 华中科技大学 Double-fed fan motor field grid-connected system sub-synchronous oscillation suppression method based on virtual resistance

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
双馈风机的次同步振荡特性与基于虚拟电阻的抑制策略研究;王玉鹏等;《可再生能源》;20180120;第36卷(第01期);全文 *

Also Published As

Publication number Publication date
CN111245013A (en) 2020-06-05

Similar Documents

Publication Publication Date Title
CN109120001B (en) Doubly-fed wind power plant grid-connected system subsynchronous oscillation suppression method based on virtual resistor
CN108631338B (en) Method for inhibiting grid-connected subsynchronous oscillation of doubly-fed wind power plant
Xu et al. Dynamic modeling and control of DFIG-based wind turbines under unbalanced network conditions
CN111245013B (en) DFIG subsynchronous oscillation suppression method based on multi-branch impedance remodeling
CN107732939B (en) Subsynchronous oscillation suppression control method based on voltage source type converter decoupling control
CN107681698A (en) Double-fed fan motor rotor string resistance low voltage traversing control method based on power optimization
CN107147132A (en) A kind of control method for suppressing double-fed fan motor unit subsynchronous resonance
Kerrouche et al. Fractional-order sliding mode control for D-STATCOM connected wind farm based DFIG under voltage unbalanced
CN106329571A (en) Operation control method of DFIG grid side and rotor side PWM converter
CN111130142B (en) Virtual capacitance control method for double-fed wind power generation grid-connected system
CN112448398A (en) Stator side analog resistance-based doubly-fed wind power plant subsynchronous oscillation suppression method
CN106972511B (en) Direct-drive wind power subsynchronous oscillation suppression method based on network side control loop optimization
CN111130136A (en) Subsynchronous oscillation suppression method based on additional virtual impedance control
CN112838589A (en) Subsynchronous oscillation suppression method for voltage source type doubly-fed wind turbine generator under series compensation power grid
Meng et al. A sub-synchronous oscillation suppression strategy for doubly fed wind power generation system
Mishra et al. Flicker attenuation using FACTS device for DFIG-based WECS connected to distribution network
CN112332427B (en) Double-fed wind power plant subsynchronous oscillation suppression method based on integral control
CN112993991A (en) Fan two-channel damping control low-frequency oscillation wide-area damping control method
CN115995828A (en) Subsynchronous oscillation suppression method and system for virtual synchronous doubly-fed fan through series compensation grid-connected system
CN104883109A (en) Control method for restraining harmonic current of doubly-fed wind generator stator side
Zhang et al. Coordinated power oscillation damping from a VSC-HVDC grid integrated with offshore wind farms: using capacitors energy
Shah et al. Impedance-based prediction of SSR-generated harmonics in doubly-fed induction generators
Li et al. VSG virtual inertial control strategy based on lead-lag link and fuzzy logic control
Mahvash et al. A look-up table based approach for fault ride-through capability enhancement of a grid connected DFIG wind turbine
CN115459300A (en) Doubly-fed wind power plant subsynchronous oscillation suppression method based on linear active disturbance rejection control

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: 20220513

CF01 Termination of patent right due to non-payment of annual fee