CN116191478B - Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors - Google Patents

Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors Download PDF

Info

Publication number
CN116191478B
CN116191478B CN202310439751.2A CN202310439751A CN116191478B CN 116191478 B CN116191478 B CN 116191478B CN 202310439751 A CN202310439751 A CN 202310439751A CN 116191478 B CN116191478 B CN 116191478B
Authority
CN
China
Prior art keywords
power system
equivalent model
frequency response
dynamic equivalent
test
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
CN202310439751.2A
Other languages
Chinese (zh)
Other versions
CN116191478A (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.)
State Grid Electric Power Research Institute Of Sepc
Xian Jiaotong University
Original Assignee
State Grid Electric Power Research Institute Of Sepc
Xian Jiaotong 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 Electric Power Research Institute Of Sepc, Xian Jiaotong University filed Critical State Grid Electric Power Research Institute Of Sepc
Priority to CN202310439751.2A priority Critical patent/CN116191478B/en
Publication of CN116191478A publication Critical patent/CN116191478A/en
Application granted granted Critical
Publication of CN116191478B publication Critical patent/CN116191478B/en
Active 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
    • H02J3/241The oscillation concerning frequency
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F3/00Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
    • G06F3/01Input arrangements or combined input and output arrangements for interaction between user and computer
    • G06F3/02Input arrangements using manually operated switches, e.g. using keyboards or dials
    • G06F3/023Arrangements for converting discrete items of information into a coded form, e.g. arrangements for interpreting keyboard generated codes as alphanumeric codes, operand codes or instruction codes
    • G06F3/027Arrangements for converting discrete items of information into a coded form, e.g. arrangements for interpreting keyboard generated codes as alphanumeric codes, operand codes or instruction codes for insertion of the decimal point
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/12Computing arrangements based on biological models using genetic models
    • G06N3/126Evolutionary algorithms, e.g. genetic algorithms or genetic programming
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/06Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
    • 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]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Biophysics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Evolutionary Biology (AREA)
  • General Physics & Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Biomedical Technology (AREA)
  • Data Mining & Analysis (AREA)
  • Physiology (AREA)
  • Genetics & Genomics (AREA)
  • Artificial Intelligence (AREA)
  • Human Computer Interaction (AREA)
  • Computational Linguistics (AREA)
  • Power Engineering (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention provides a method for evaluating equivalent inertia of a plurality of asynchronous motors and modeling frequency response, belonging to the technical field of frequency stability analysis of an electric power system; the problem that when a plurality of asynchronous motors are connected into a power system, inertia evaluation and frequency response modeling of the power system are difficult is solved; the method comprises the following steps: step one: acquiring frequency response data of a power system to be analyzed through a power grid SCADA system, and dividing historical data into a fitting set and a testing set; step two: establishing a dynamic equivalent model of the power system to be analyzed, wherein the dynamic equivalent model adopts a low-order power system model in a transfer function form; step three: fitting parameters of a dynamic equivalent model from the test set through a genetic algorithm; step four: testing whether the parameters of the established dynamic equivalent model meet the test indexes or not on the test set; the invention is applied to stability analysis of an electric power system comprising a plurality of asynchronous motors.

Description

Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors
Technical Field
The invention provides a method for evaluating equivalent inertia of a plurality of asynchronous motors and modeling frequency response, and belongs to the technical field of frequency stability analysis of power systems.
Background
In recent years, more and more new energy power generation systems are connected into a power system through a power electronic converter, and the ratio of the photovoltaic power station to the wind turbine generator in the new energy is more than 80%. The photovoltaic power station does not have rotational inertia, and the power electronic converter makes the inertia of the wind turbine generator difficult to effectively utilize. Therefore, a large number of conventional generator sets are replaced by new energy sources, which causes the moment of inertia of the power system to gradually decrease, causing concerns about frequency stability.
Under the above circumstances, the asynchronous motor directly connected to the power grid has the capability of participating in the frequency response, and can improve primary frequency modulation performance, so that the asynchronous motor has an increasingly important role in the frequency response.
However, when evaluating the inertia of an asynchronous motor and establishing a mechanism model of the frequency response of an electric power system, detailed parameters of the asynchronous motor are often required, for example, an equivalent inertia evaluation method of the asynchronous motor in an inertia response phase disclosed in the patent application No. 2022105164271, but it is very difficult to implement the equivalent inertia evaluation method for a large system having hundreds of asynchronous motors.
Disclosure of Invention
The invention provides a method for evaluating equivalent inertia and modeling frequency response of a plurality of asynchronous motors, which aims to solve the problem that the inertia evaluation and the frequency response modeling of a power system are difficult when the plurality of asynchronous motors are connected into the power system.
In order to solve the technical problems, the invention adopts the following technical scheme: a method for evaluating equivalent inertia and modeling frequency response of a plurality of asynchronous motors comprises the following steps:
step one: acquiring frequency response data of a power system to be analyzed through a power grid SCADA system, and dividing historical data into a fitting set and a testing set;
step two: establishing a dynamic equivalent model of the power system to be analyzed, wherein the dynamic equivalent model adopts a low-order power system model in a transfer function form;
step three: fitting parameters of a dynamic equivalent model from the test set through a genetic algorithm;
step four: and testing whether the parameters of the established dynamic equivalent model meet the test indexes or not on the test set.
The change rule of the frequency deviation of the dynamic equivalent model established in the second step is as follows:
Figure SMS_1
in the above formula:Din order to be a damping coefficient,His an inertial time constant of the power system;
ΔP d for load disturbance, the expression of load disturbance is:
Figure SMS_2
in the above formula:P d,mag the amplitude of the load disturbance is a step function;
ΔP g for the output variation of the synchronous generator, the expression of the output variation of the synchronous generator is as follows:
Figure SMS_3
in the above formula:T R for the reheat time constant, the temperature of the refrigerant is,K m as a function of the mechanical power gain factor,F H for the high-pressure turbine coefficient,Ris a governor coefficient.
In the third step, parameters of the dynamic equivalent model are fitted from the test set through a genetic algorithm, wherein the parameters of the dynamic equivalent model comprise: speed regulator coefficientRInertial time constant of power systemHDamping coefficientDCoefficient of mechanical power gainK m High pressure turbine coefficientF H And reheat time constantT R
The expression of the optimization target of the genetic algorithm is as follows:
Figure SMS_4
in the above formula:Jin order to adapt the function of the degree of adaptation,mas the number of sets of data to be used for fitting,n i for the number of samples of the power system frequency in each set of data,
Figure SMS_5
data for true frequency offset, +.>
Figure SMS_6
Is the frequency offset of the dynamic equivalent model.
Constraint conditions to be met by the optimization target of the genetic algorithm are as follows:
Figure SMS_7
the expression of the test index in the fourth step is as follows:
Figure SMS_8
in the above formula:J test for the test index, used to characterize the average fitting error,lfor the number of test sets, the standard for test passing is that the test index is less than 1%.
Compared with the prior art, the invention has the following beneficial effects: the method for evaluating the equivalent inertia of the plurality of asynchronous motors and modeling the frequency response aims to establish a dynamic equivalent model, adopts a genetic algorithm to fit parameters of the dynamic equivalent model from historical data, is convenient and quick, does not need detailed asynchronous motor parameters, and can be suitable for evaluating the inertia and modeling the frequency response of a large-scale power system comprising the plurality of asynchronous motors.
Drawings
The invention is further described below with reference to the accompanying drawings:
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a schematic diagram of an IEEE9 node system with three asynchronous motors in accordance with an embodiment of the present invention;
FIG. 3 is a schematic illustration of a Simulink modeling architecture of the IEEE9 node system with three asynchronous motors of FIG. 2;
FIG. 4 is a schematic diagram of a transfer function of a dynamic equivalent model of the frequency response of the power system according to the present invention;
FIG. 5 is a schematic diagram of the power system frequency response under three different sets of perturbations used for fitting in accordance with the present invention;
FIG. 6 is a graph showing the frequency response of the power system under four different sets of perturbations used in the test according to the present invention.
Detailed Description
The invention aims to provide a method for evaluating equivalent inertia and modeling frequency response of a plurality of asynchronous motors driven by data, which comprises the following implementation steps of:
step one: acquiring frequency response data of a power system to be analyzed through a power grid SCADA system, and dividing historical data into a fitting set and a testing set;
step two: establishing a dynamic equivalent model of the power system to be analyzed, wherein the dynamic equivalent model adopts a low-order power system model in a transfer function form;
step three: fitting parameters of a dynamic equivalent model from the test set through a genetic algorithm;
step four: and testing the established dynamic equivalent model parameters on a test set, wherein the standard of test passing is that the test index is less than 1%.
The method of the invention is further described below with reference to specific examples.
Embodiment one: the schematic diagram of the power system to be tested is shown in fig. 2, the Simulink model of the power system to be tested is shown in fig. 3, and the Simulink model is a frequency response historical data source of the power system to be tested.
Historical frequency response data of a power system comprising a plurality of asynchronous motors is obtained, and the historical data is divided into three fitting sets and four testing sets.
The dynamic equivalent model of the frequency response of the power system is shown in fig. 4, a dynamic equivalent model of the power system to be analyzed is established, and the dynamic equivalent model adopts a low-order power system model in a transfer function form. In this model, the law of variation of the frequency deviation is:
Figure SMS_9
in the method, in the process of the invention,Dthe damping coefficient is generally 0 to 2,Hthe inertia time constant of the power system is generally 3-9 s.
ΔP d For load disturbance, the expression of load disturbance is:
Figure SMS_10
in the method, in the process of the invention,P d,mag the amplitude of the load disturbance is a step function.
ΔP g For the output variation of the synchronous generator, the expression of the output variation of the synchronous generator is
Figure SMS_11
In the method, in the process of the invention,T R for the reheat time constant, it is usually between 6 to 14s,K m the mechanical power gain factor is mainly influenced by the power factor and the rotation reserve, and is generally between 0.4 and 1,F H the high-pressure turbine coefficient is generally between 0.15 and 0.4,Rthe speed regulator coefficient is generally 0.04-1.
Fitting parameters of a dynamic equivalent model from the test set through a genetic algorithm, wherein the parameters of the dynamic equivalent model comprise: speed regulator coefficientRInertial time constant of power systemHDamping coefficientDCoefficient of mechanical power gainK m High pressure turbine coefficientF H And reheat time constantT R
The optimization targets of the genetic algorithm are:
Figure SMS_12
in the method, in the process of the invention,Jin order to adapt the function of the degree of adaptation,mas the number of sets of data to be used for fitting,n i for the number of samples of the power system frequency in each set of data,
Figure SMS_13
data for true frequency offset, +.>
Figure SMS_14
Is the frequency offset of the dynamic equivalent model. The constraint conditions to be satisfied by the above optimization problem are:
Figure SMS_15
i.e. all parameters are guaranteed to be within reasonable limits.
The genetic algorithm fits on three fitting sets, and the frequency response of the obtained equivalent dynamic model and the equivalent dynamic model pair of the power system to be analyzed are as shown in fig. 5.
And testing the established dynamic equivalent model parameters on a test set, wherein the test indexes are as follows:
Figure SMS_16
in the method, in the process of the invention,J test for the test index, used to characterize the average fitting error,lfor the number of test sets, the standard for test passing is that the test index is less than 1%.
The equivalent dynamic model and the test set result of the power system to be analyzed are shown in fig. 6.
According to the invention, the inertia of the power system can be evaluated and a frequency response model can be established under the condition that the detailed parameters of the asynchronous motor are unknown, so that the frequency stability of the power system can be effectively analyzed.
The specific structure of the invention needs to be described that the connection relation between the component modules adopted by the invention is definite and realizable, and besides the specific description in the embodiment, the specific connection relation can bring about corresponding technical effects, and on the premise of not depending on execution of corresponding software programs, the technical problems of the invention are solved, the types of the components, the modules and the specific components, the connection modes of the components and the expected technical effects brought by the technical characteristics are clear, complete and realizable, and the conventional use method and the expected technical effects brought by the technical characteristics are all disclosed in patents, journal papers, technical manuals, technical dictionaries and textbooks which can be acquired by a person in the field before the application date, or the prior art such as conventional technology, common knowledge in the field, and the like, so that the provided technical scheme is clear, complete and the corresponding entity products can be reproduced or obtained according to the technical means.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and not for limiting the same; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some or all of the technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit of the invention.

Claims (2)

1. A method for evaluating equivalent inertia and modeling frequency response of a plurality of asynchronous motors is characterized by comprising the following steps: the method comprises the following steps:
step one: acquiring frequency response data of a power system to be analyzed through a power grid SCADA system, and dividing historical data into a fitting set and a testing set;
step two: establishing a dynamic equivalent model of the power system to be analyzed, wherein the dynamic equivalent model adopts a low-order power system model in a transfer function form;
step three: fitting parameters of the dynamic equivalent model from the fitting set through a genetic algorithm, wherein the parameters of the dynamic equivalent model comprise: speed regulator coefficientRInertial time constant of power systemHDamping coefficientDCoefficient of mechanical power gainK m High pressure turbine coefficientF H And reheat time constantT R
The expression of the optimization target of the genetic algorithm is as follows:
Figure QLYQS_1
in the above formula:Jin order to adapt the function of the degree of adaptation,mas the number of sets of data to be used for fitting,n i for the number of samples of the power system frequency in each set of data,
Figure QLYQS_2
data for true frequency offset, +.>
Figure QLYQS_3
Frequency offset for dynamic equivalent model;
constraint conditions to be met by the optimization target of the genetic algorithm are as follows:
Figure QLYQS_4
step four: testing whether the parameters of the established dynamic equivalent model meet the test indexes or not on the test set;
the expression of the test index in the fourth step is as follows:
Figure QLYQS_5
in the above formula:J test for the test index, used to characterize the average fitting error,lfor the number of test sets, the standard for test passing is that the test index is less than 1%.
2. A method for evaluating equivalent inertia and modeling frequency response of a plurality of asynchronous motors according to claim 1, wherein: the change rule of the frequency deviation of the dynamic equivalent model established in the second step is as follows:
Figure QLYQS_6
in the above formula:Din order to be a damping coefficient,His an inertial time constant of the power system;
ΔP d for load disturbance, the expression of load disturbance is:
Figure QLYQS_7
in the above formula:P d,mag is the amplitude of the load disturbance;
ΔP g for the output variation of the synchronous generator, the expression of the output variation of the synchronous generator is as follows:
Figure QLYQS_8
in the above formula:T R for the reheat time constant, the temperature of the refrigerant is,K m as a function of the mechanical power gain factor,F H for the high-pressure turbine coefficient,Ris a governor coefficient.
CN202310439751.2A 2023-04-23 2023-04-23 Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors Active CN116191478B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310439751.2A CN116191478B (en) 2023-04-23 2023-04-23 Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310439751.2A CN116191478B (en) 2023-04-23 2023-04-23 Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors

Publications (2)

Publication Number Publication Date
CN116191478A CN116191478A (en) 2023-05-30
CN116191478B true CN116191478B (en) 2023-07-11

Family

ID=86452381

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310439751.2A Active CN116191478B (en) 2023-04-23 2023-04-23 Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors

Country Status (1)

Country Link
CN (1) CN116191478B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113972654A (en) * 2021-11-01 2022-01-25 国网江苏省电力有限公司电力科学研究院 Multipoint-access area load inertia time constant identification method
CN115483707A (en) * 2022-10-12 2022-12-16 东南大学溧阳研究院 Novel power system frequency situation prediction method considering photovoltaic frequency modulation

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111276973A (en) * 2020-03-09 2020-06-12 国网江苏省电力有限公司 Method for evaluating inertia requirement of power system considering wind power fluctuation
CN114638074A (en) * 2021-12-31 2022-06-17 国网辽宁省电力有限公司电力科学研究院 Inertia evaluation method based on quantum derivation algorithm
CN114792055A (en) * 2022-05-05 2022-07-26 山西大学 Asynchronous motor equivalent inertia evaluation method based on transient reactance post-potential
CN114841005A (en) * 2022-05-12 2022-08-02 南京理工大学 Method for evaluating equivalent inertia of asynchronous motor in inertia response stage
CN115622149A (en) * 2022-09-06 2023-01-17 华中科技大学 System frequency response modeling method and system for double-fed fan participating in primary frequency modulation

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113972654A (en) * 2021-11-01 2022-01-25 国网江苏省电力有限公司电力科学研究院 Multipoint-access area load inertia time constant identification method
CN115483707A (en) * 2022-10-12 2022-12-16 东南大学溧阳研究院 Novel power system frequency situation prediction method considering photovoltaic frequency modulation

Also Published As

Publication number Publication date
CN116191478A (en) 2023-05-30

Similar Documents

Publication Publication Date Title
CN103887815B (en) Based on wind energy turbine set parameter identification and the Dynamic Equivalence of service data
Nadour et al. Comparative analysis between PI & backstepping control strategies of DFIG driven by wind turbine
Leithead et al. Control of variable speed wind turbines: Design task
Zhong et al. A novel frequency regulation strategy for a PV system based on the curtailment power-current curve tracking algorithm
CN110492479B (en) Method for identifying rotational inertia and damping of distributed grid-connected equipment
Mensou et al. An efficient nonlinear Backstepping controller approach of a wind power generation system based on a DFIG
CN112542855B (en) Modeling and simulation method for phasor model of double-fed wind power generation system
CN110829487A (en) Dynamic frequency prediction method for power system
CN111725840B (en) Parameter identification method for direct-drive wind generating set controller
CN107947228B (en) Stochastic stability analysis method for power system containing wind power based on Markov theory
CN111756038B (en) New energy power system equal frequency difference inertia estimation method considering frequency modulation characteristics
Li et al. Research on clustering equivalent modeling of large-scale photovoltaic power plants
CN111311021A (en) Theoretical power prediction method, device, equipment and storage medium for wind power plant
CN116738636A (en) Multi-machine equivalent method of doubly-fed wind power station considering impedance characteristics and synchronization mechanism
CN115578016A (en) Online evaluation method for frequency modulation capability of wind power plant with incomplete model
Weijie et al. Investigating instability of the wind turbine simulator with the conventional inertia emulation scheme
CN110829491A (en) Grid-connected photovoltaic power generation system parameter identification method based on transient disturbance
Beniss et al. Improvement of Power Quality Injected into the Grid by Using a FOSMC-DPC for Doubly Fed Induction Generator.
CN104617578A (en) Method for acquiring available power transmission capability of power system with wind power plant
CN116191478B (en) Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors
Yan et al. Transient modelling of doubly‐fed induction generator based wind turbine on full operation condition and rapid starting period based on low voltage ride‐through testing
Jin et al. Maximum power control of wind turbines with practical prescribed time stability based on wind estimation
CN109657380A (en) A kind of double-fed fan motor field Dynamic Equivalence based on Extended Kalman filter
Bouguerra Comparative study between PI, FLC, SMC and Fuzzy sliding mode controllers of DFIG wind turbine
Manjeera et al. Design and Implementation of Fuzzy logic-2DOF controller for Emulation of wind turbine System

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