CN108717593B - Micro-siting power generation capacity evaluation method based on equivalent wind speed of wind wheel surface - Google Patents
Micro-siting power generation capacity evaluation method based on equivalent wind speed of wind wheel surface Download PDFInfo
- Publication number
- CN108717593B CN108717593B CN201810335465.0A CN201810335465A CN108717593B CN 108717593 B CN108717593 B CN 108717593B CN 201810335465 A CN201810335465 A CN 201810335465A CN 108717593 B CN108717593 B CN 108717593B
- Authority
- CN
- China
- Prior art keywords
- wind
- rews
- power plant
- shear
- power generation
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
Abstract
A micro-siting power generation capacity evaluation method based on wind wheel surface equivalent wind speed comprises the following steps: a) calculating wind shear at each wind turbine location in the wind farm: when wind shear of each machine position point of the wind power plant is calculated, a set of CFD simulation calculation method is established according to the air flow theory of the wind power plant, and meanwhile wind shear of each machine position point of the wind power plant is calculated more accurately by combining with actually measured data of a wind measuring tower; b) the influence of wind shear at the position of each wind turbine generator on the generated energy is considered, and the generated energy of the wind power plant is evaluated: converting the conventional dynamic power curve into a dynamic power curve represented by REWS when evaluating the power generation capacity of the wind farm; according to the concept of the REWS, calculating the REWS aiming at wind shear of each locomotive position point of the wind power plant; and evaluating the generated energy of each wind turbine generator and the wind power plant according to the dynamic power curve based on the REWS, the REWS and the Weibull distribution. The method has the advantages of high result precision, practical theory, concise calculation flow and high efficiency.
Description
Technical Field
The invention relates to the technical field of micro site selection of wind power plants, in particular to a micro site selection power generation amount evaluation method based on equivalent wind speed of a wind wheel surface.
Background
The micro-site selection power generation evaluation is a core link of wind power plant development and plays an important role in the success or failure of wind power plant benefits and wind power plant investment. With the gradual reduction of the wind power price on the internet in China, the development of wind power plants gradually turns to regions in the southeast mountains with low wind speed and complex terrain, and the profit and loss of the wind power plants almost touch the balance point, so that more rigorous requirements are provided for the accuracy of the generated energy evaluation method.
The current generated energy evaluation method is mainly obtained by calculating the Hub Height Wind Speed (HHWS), Wind frequency distribution and a dynamic power curve. In calculating the dynamic power curve, a wind condition input condition of wind shear (vertical shear of horizontal wind) needs to be considered. Typically, wind shear is assumed to be an exponential function, and the model index defaults to 0.2 (exponential wind shear). However, due to the effects of terrain and thermal variations, the actual wind shear of a wind farm typically varies widely and does not follow an exponential law. In addition, the current dynamic power curve is defined by the wind speed at the height of the hub, and the wind speed at the height of the hub cannot represent the change of the wind speed of the whole wind wheel surface, so that the dynamic power curve and the annual energy production amount calculated by the current method have larger difference compared with the actual situation. As the gravity center of wind power development in China turns to a region with low wind speed and complex mountainous regions, the diameter of an impeller of a wind turbine generator is continuously increased, and the influence of wind shear on the wind sweeping surface is increasingly and non-negligibly avoided.
Disclosure of Invention
In order to overcome the defects of lack of credibility of theory and large error of result caused by not considering the actual wind shear of a wind power plant in the existing micro-siting power generation evaluation technology, the invention provides a power generation evaluation method with high result precision, practical theory, simple calculation flow and high efficiency based on the concepts of CFD technology and REWS.
The technical scheme adopted by the invention for solving the technical problems is as follows:
a micro-siting power generation capacity evaluation method based on wind wheel surface equivalent wind speed comprises the following steps:
a) calculating wind shear at each wind turbine location within a wind farm
When wind shear of each machine position point of the wind power plant is calculated, a set of CFD simulation calculation method is established according to the air flow theory of the wind power plant, and meanwhile wind shear of each machine position point of the wind power plant is calculated more accurately by combining with actually measured data of a wind measuring tower;
b) the influence of wind shear at the position of each wind turbine generator on the generated energy is considered, and the generated energy of the wind power plant is evaluated
Converting the conventional dynamic power curve into a dynamic power curve represented by REWS when evaluating the power generation capacity of the wind farm; according to the concept of the REWS, calculating the REWS aiming at wind shear of each locomotive position point of the wind power plant; and evaluating the generated energy of each wind turbine generator and the wind power plant according to the dynamic power curve of the REWS, the REWS and the Weibull distribution.
Further, in the step a), a wind power plant air flow theory and strategy are given, and CFD software is used for carrying out numerical simulation on a wind power plant flow field with a plurality of inlet wind directions to obtain wind shear of each machine site of the wind power plant of each sector; and obtaining the actual comprehensive wind shear of each machine site of the wind power plant in a set time period by a weighted average method according to the measured data of the wind measuring tower, the wind speed acceleration factor, the sector interpolation theory and the wind direction distribution frequency. The wind shear calculation method is verified to be more consistent with an actual test result.
In the step b), an exponential type wind shear with an index of 0.2, a specific turbulence intensity and an air density are used as input conditions to calculate a conventional dynamic power curve, and the dynamic power curve is represented by HHWS; and finally obtaining a dynamic power curve expressed by each machine position point through the REWS according to the relation between the HHWS and the REWS.
In the step b), according to the concept of REWS, combining the wind shear of each machine position of the wind wheel rotating surface and the wind farm to calculate the REWS, wherein the calculation formula is
Wherein i represents the height, n represents the number of height layers, UiRepresents the wind speed at altitude i, AiRepresents the area of the i-th section, and a represents the area of the rotor rotation plane.
In the step b), the annual generating capacity of each wind turbine generator and the wind farm considering the wind shear influence is calculated according to the dynamic power relation curve of the REWS, the REWS and the Weibull distribution.
The technical conception of the invention is as follows: the method is characterized in that a power generation amount evaluation method is improved in an innovative mode, Wind shear of each Wind turbine of a Wind power plant is accurately calculated by adopting a CFD (computational Fluid dynamics) technology, and a dynamic power curve is defined by adopting Rotor Equivalent Wind Speed (REWS), so that the influence of actual Wind shear is considered when the power generation amount of the Wind power plant is evaluated. The method is high in calculation accuracy, simple in calculation process and significant in improvement of the micro-siting power generation amount evaluation accuracy.
The invention has the following beneficial effects: the wind power plant wind shear can be obtained by calculating by using the CFD technology, and a scientific and accurate tool is provided for the evaluation of the generated energy, so that the cost of actual evaluation is greatly reduced. The REWS innovation theory is used for expressing the dynamic power curve and evaluating the power generation amount, the influence of actual wind shear can be considered during power generation amount evaluation, and evaluation precision is greatly improved. The calculation flow is simple, and compared with the traditional calculation method, the extra calculation cost is not increased.
Drawings
FIG. 1 is a schematic representation of the calculation of REWS modeling in conjunction with wind turbine face and wind shear in the present invention;
fig. 2 is a flow chart of a method for evaluating micro-siting power generation based on REWS.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
Referring to fig. 1 and 2, a micro-siting power generation evaluation method based on wind turbine equivalent wind speed includes the following steps:
a) calculating wind shear at each wind turbine location within a wind farm
When wind shear of each machine position point of the wind power plant is calculated, a set of CFD simulation calculation method is established according to the air flow theory of the wind power plant, and meanwhile, the actual wind shear of each machine position point of the wind power plant is calculated more accurately by combining the wind speed acceleration factor and the wind direction sector interpolation theory;
b) the influence of wind shear at the position of each wind turbine generator on the generated energy is considered, and the generated energy of the wind power plant is evaluated
Converting the conventional dynamic power curve into a dynamic power curve represented by REWS when evaluating the power generation capacity of the wind farm; according to the concept of the REWS, calculating the REWS aiming at wind shear of each locomotive position point of the wind power plant; and evaluating the generated energy of each wind turbine generator and the wind power plant according to the dynamic power curve of the REWS, the REWS and the Weibull distribution.
Further, in the step a), a wind power plant air flow theory and strategy are given, and CFD software is used for carrying out numerical simulation on a wind power plant flow field with a plurality of inlet wind directions to obtain wind shear of each machine site of the wind power plant of each sector; and obtaining the actual comprehensive wind shear of each machine site of the wind power plant in a set time period by a weighted average method according to the measured data of the wind measuring tower, the wind speed acceleration factor, the sector interpolation theory and the wind direction distribution frequency.
In the step b), an exponential type wind shear with an index of 0.2, a specific turbulence intensity and an air density are used as input conditions to calculate a conventional dynamic power curve, and the dynamic power curve is represented by HHWS; and finally obtaining a dynamic power curve expressed by each machine position point through the REWS according to the relation between the HHWS and the REWS.
In the step b), according to the concept of the REWS, combining the wind shear of each machine position of the wind wheel rotating surface and the wind farm, and calculating the REWS.
In the step b), the annual generating capacity of each wind turbine generator and the wind farm considering the wind shear influence is calculated according to the dynamic power relation curve of the REWS, the REWS and the Weibull distribution.
The method for evaluating the micro addressing power generation amount based on the REWS comprises the following steps: 1. preparing in an early stage; 2. CFD orientation calculation (multi-sector); 3. calculating wind shear and REWS of each machine position; 4. converting the conventional dynamic power curve into a dynamic power curve based on REWS; 5. and evaluating the power generation capacity.
Firstly, CFD directional calculation is carried out on the premise that a geometric file and a grid file are prepared and various parameter conditions are set on a simulation platform; secondly, obtaining the REWS on the premise that the actual wind shear of each machine position point is obtained; the premise of obtaining the dynamic power curve based on the REWS is that a conventional dynamic power curve based on the HHWS and the relationship between the HHWS and the REWS are obtained by calculation; the anemometer tower anemometry data processing can provide a basis for verification of the correctness of the simulation technology and evaluation of the generated energy; some parameters, eventually derived from the weibull distribution, may provide a necessary input condition for evaluating the power generation.
The early-stage preparation comprises a terrain geometric file, a grid file and parameter setting, wherein the terrain file is mainly a geometric parameter modeling file, the grid file is obtained through grid generation software, and the parameter setting is mainly carried out in numerical simulation software correspondingly according to the air flow theory of the wind power plant.
The CFD directional calculation is as shown in FIG. 2, numerical simulation is carried out on a wind power plant flow field of a plurality of inlet wind directions by using a CFD numerical calculation method, and wind shear of each machine position point of the wind power plant in each sector is obtained; and obtaining the actual wind shear of each machine position point of the wind power plant in the specific time period by a weighted average method according to the actually measured data of the wind measuring tower, the wind speed acceleration factor, the sector interpolation theory and the wind direction distribution frequency.
The calculation REWS is calculated according to the REWS theory for the wind turbine surface and the wind shear of the wind turbine generator, as shown in fig. 1. As shown in fig. 1, the distribution of the wind speed in the vertical direction of the wind wheel plane can be determined by U ═ UhIndicates that (h is the height from the ground), and the area infinitesimal corresponding to each height passes through AhAnd (4) showing. According to the REWS formula, multiplying the ratio of the area infinitesimal of the area of the wind wheel face at different heights in the area of the wind wheel face by the wind speed at the infinitesimal, summing the items at different heights in the area, and then opening the sum to obtain the REWS. The formula for calculating REWS is
Wherein i represents the height, n represents the number of height layers, UiRepresents the wind speed at altitude i, AiRepresents the area of the i-th section, and a represents the area of the rotor rotation plane.
The conventional dynamic power curve is shown in fig. 2 and is calculated by using an exponential wind shear with an index of 0.2, a specific turbulence intensity and an air density as input conditions, and is expressed by HHWS.
The REWS-based dynamic power curve is shown in fig. 2 and is obtained by converting a conventional dynamic power curve. And converting the relation curve of the power and the HHWS into a dynamic power curve based on the REWS according to the relation between the HHWS and the REWS, and finally obtaining the dynamic power curve based on the REWS of each locomotive position of the wind power plant.
The anemometer tower anemometry data processing module can be combined with wind shear of each sector of each machine site of the wind power plant obtained by CFD simulation as shown in FIG. 2, and finally comprehensive wind shear of each machine site is obtained through calculation.
According to the technical scheme, the method for evaluating the power generation amount of the micro site selection based on the REWS has clear logic, considers the influence of the actual wind shear of each unit in the wind power plant on the power generation amount, and can effectively meet the requirement for evaluating the accuracy of the power generation amount of the wind power plant. The wind flow field CFD simulation technology is combined with the measured data of the anemometer tower to obtain the actual wind shear of each machine site, the influence of more actual factors on the generated energy of the wind power plant is considered, and the wind shear is closer to the actual wind; the REWS innovation theory is used for expressing a dynamic power curve and evaluating the generated energy, the influence of actual wind shear can be considered during the evaluation of the generated energy, and the evaluation precision of the generated energy is greatly improved; the method can be applied to units with different models and different capacities, can be applied to wind power plants with different landforms, and has universality; the calculation process is simple, the theory is practical, the calculation process can be directly processed through a program, the problems of long calculation time and low efficiency are avoided, and the method can be widely applied to the field of wind power.
While the invention has been described in conjunction with specific embodiments thereof, it is not intended to limit the invention to the particular embodiments described, and alternatives, modifications, and variations, or equivalents may be substituted for those skilled in the art based on the foregoing description. Any substitutions, modifications, equivalents and the like which are within the spirit and principle of the present invention shall be included within the scope of the present invention.
Claims (1)
1. A micro-siting power generation capacity evaluation method based on wind wheel surface equivalent wind speed is characterized by comprising the following steps: the method comprises the following steps:
a) calculating wind shear at each wind turbine location within a wind farm
When wind shear of each machine position of the wind power plant is calculated, weighted average solving is carried out on wind shear of each sector and wind direction distribution frequency obtained by CFD simulation, and wind shear of each machine position of the wind power plant with high precision is obtained;
b) the influence of wind shear at the position of each wind turbine generator on the generated energy is considered, and the generated energy of the wind power plant is evaluated
Converting the conventional dynamic power curve into a dynamic power curve represented by REWS when evaluating the power generation capacity of the wind farm; calculating an REWS aiming at wind shear of each machine position point of a wind power plant according to the concept of equivalent wind speed REWS of a wind wheel surface; according to the dynamic power curve of the REWS, the REWS and the Weibull distribution, the power generation amount of each wind turbine generator and the power generation amount of the wind power plant are evaluated;
in the a), a wind power plant air flow theory and strategy are given, and CFD software is used for carrying out numerical simulation on a wind power plant flow field with a plurality of inlet wind directions to obtain wind shear of each machine site of the wind power plant of each sector; according to measured data of a wind measuring tower, a wind speed acceleration factor, a sector interpolation theory and wind direction distribution frequency, obtaining actual comprehensive wind shear of each machine site of a wind power plant in a set time period by a weighted average method;
in the b), an exponential type wind shear with an index of 0.2, turbulence intensity and air density are used as input conditions to calculate a conventional dynamic power curve, and the dynamic power curve is represented by a hub height wind speed HHWS; according to the relation between HHWS and REWS, finally obtaining a dynamic power curve expressed by each machine position point through REWS;
in the b), according to the concept of REWS, combining the wind shear of each machine position of the wind wheel rotating surface and the wind farm, calculating the REWS by the formula
Wherein i represents heightN represents the number of layers of height, UiRepresents the wind speed at altitude i, AiRepresents the area of the i-th part, and A represents the area of the rotating surface of the wind wheel;
and in the step b), according to the dynamic power relation curve of the REWS, the REWS and the Weibull distribution, calculating by a program to obtain annual power generation quantities of each wind turbine and the wind power plant considering wind shear influence.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810335465.0A CN108717593B (en) | 2018-04-16 | 2018-04-16 | Micro-siting power generation capacity evaluation method based on equivalent wind speed of wind wheel surface |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810335465.0A CN108717593B (en) | 2018-04-16 | 2018-04-16 | Micro-siting power generation capacity evaluation method based on equivalent wind speed of wind wheel surface |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108717593A CN108717593A (en) | 2018-10-30 |
| CN108717593B true CN108717593B (en) | 2021-07-16 |
Family
ID=63899171
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810335465.0A Active CN108717593B (en) | 2018-04-16 | 2018-04-16 | Micro-siting power generation capacity evaluation method based on equivalent wind speed of wind wheel surface |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108717593B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111859662B (en) * | 2020-07-17 | 2023-09-22 | 泰州市都瑞堡船舶设备有限公司 | Layout method of uniform velocity tube flowmeter in tail end unit of variable-air-volume air conditioner |
| CN111967205A (en) * | 2020-08-14 | 2020-11-20 | 中国华能集团清洁能源技术研究院有限公司 | Anemometer tower micro-addressing method based on wind acceleration factor |
| CN113792430B (en) * | 2021-09-14 | 2023-08-15 | 华能陇东能源有限责任公司 | Wind shear index calculation method and system for wind power plant in complex terrain |
| CN115238344A (en) * | 2022-07-13 | 2022-10-25 | 北京瑞科同创能源科技有限公司 | Method for determining optimal hub height according to provinces based on wind shear |
| CN116415838B (en) * | 2023-02-14 | 2024-01-19 | 中节能风力发电股份有限公司 | Old wind power plant reconstruction project power generation capacity assessment method based on operation data |
| CN116169727B (en) * | 2023-02-14 | 2023-11-17 | 中节能风力发电股份有限公司 | Old wind farm reconstruction project generating capacity assessment method and system without anemometry data |
| CN117521282B (en) * | 2023-11-07 | 2024-04-12 | 国家气候中心 | Density-dependent fan parameterization method for wind farm climate characteristic simulation |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105587474A (en) * | 2015-12-10 | 2016-05-18 | 中国电力科学研究院 | Method and device for testing equivalent wind speed of wind wheels of wind turbine generator systems |
| EP2719893B1 (en) * | 2012-10-09 | 2017-01-04 | ALSTOM Renewable Technologies | Method of operating a variable speed wind turbine |
| CN107885959A (en) * | 2017-12-06 | 2018-04-06 | 华北电力大学 | A kind of wind-powered electricity generation modeling and performance estimating method based on confidence equivalent power curve belt |
-
2018
- 2018-04-16 CN CN201810335465.0A patent/CN108717593B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2719893B1 (en) * | 2012-10-09 | 2017-01-04 | ALSTOM Renewable Technologies | Method of operating a variable speed wind turbine |
| CN105587474A (en) * | 2015-12-10 | 2016-05-18 | 中国电力科学研究院 | Method and device for testing equivalent wind speed of wind wheels of wind turbine generator systems |
| CN107885959A (en) * | 2017-12-06 | 2018-04-06 | 华北电力大学 | A kind of wind-powered electricity generation modeling and performance estimating method based on confidence equivalent power curve belt |
Non-Patent Citations (1)
| Title |
|---|
| 风力发电场微观选址与机型优化综合研究;孙永岗;《中国优秀硕士学位论文全文数据库（电子期刊）》;20101215;摘要、第4、22、28-29、32-33、36-37、50、54-56页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108717593A (en) | 2018-10-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108717593B (en) | Micro-siting power generation capacity evaluation method based on equivalent wind speed of wind wheel surface | |
| CN107194097B (en) | Analysis method based on wind power plant aerodynamic simulation and wind speed and direction data | |
| CN104699936A (en) | Sector management method based on CFD short-term wind speed forecasting wind power plant | |
| CN107292514B (en) | Site selection method and device for anemometer tower in production operation of wind power plant | |
| CN112818590A (en) | Method for evaluating power generation capacity of wind power plant and establishing micro-site selection model | |
| Schepers et al. | Analysis of Mexico wind tunnel measurements: Final report of IEA Task 29, Mexnext (Phase 1) | |
| CN113205210B (en) | Wind power plant wind speed and power prediction method, system and equipment for complex terrain and storage medium | |
| CN114169614B (en) | Wind power plant optimal scheduling method and system based on wind turbine wake model optimization | |
| CN106777499B (en) | Complete machine dynamic modeling method for double-fed asynchronous wind generating set | |
| Ali et al. | Effects of instantaneous tangential velocity on the aerodynamic performance of an H-Darrieus wind turbine | |
| CN114169197A (en) | Wind turbine generator on-loop simulation system and method | |
| CN106951977B (en) | Construction method of wind speed prediction model based on wake effect | |
| CN113627101A (en) | Wind turbine wake flow simulation method based on improved AD/RSM model | |
| CN109636019B (en) | Wind measuring tower arrangement scheme determination method based on neural network algorithm | |
| CN113657055A (en) | Inflow condition generation method and system for numerical simulation of wind field in complex terrain | |
| WO2023165159A1 (en) | Online indirect measurement system and method for pitching and yawing moments of wind energy or tidal current energy generator set | |
| CN109446548B (en) | Marine wind field automation machine position arrangement software and operation method thereof | |
| CN102288159A (en) | Method and device for determining the deflection of a tower | |
| Yi et al. | Airfoil design for vertical axis wind turbine operating at variable tip speed ratios | |
| Francis et al. | Development of a novel airfoil for low wind speed vertical axis wind turbine using QBlade simulation tool | |
| CN109325248A (en) | Three-dimensional wind speed profile associated diagram method for building up in wind power plant | |
| CN115310375A (en) | Wind turbine generator cabin transfer function fitting method based on fluid dynamics model | |
| CN114186407A (en) | Wind power plant wake velocity field calculation method and system with parameters capable of being adjusted in self-adaptive mode | |
| CN112861301A (en) | Wind power plant theoretical power intelligent calculation method based on real-time data of fans | |
| CN111222233A (en) | Angular motion based motion simulation method for offshore wind turbine generator system platform |
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 |