CN107947246B - Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power - Google Patents

Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power Download PDF

Info

Publication number
CN107947246B
CN107947246B CN201711126400.7A CN201711126400A CN107947246B CN 107947246 B CN107947246 B CN 107947246B CN 201711126400 A CN201711126400 A CN 201711126400A CN 107947246 B CN107947246 B CN 107947246B
Authority
CN
China
Prior art keywords
wind power
wind
provincial
output
power
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
CN201711126400.7A
Other languages
Chinese (zh)
Other versions
CN107947246A (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.)
Tsinghua University
North China Grid Co Ltd
Original Assignee
Tsinghua University
North China Grid Co Ltd
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 Tsinghua University, North China Grid Co Ltd filed Critical Tsinghua University
Priority to CN201711126400.7A priority Critical patent/CN107947246B/en
Publication of CN107947246A publication Critical patent/CN107947246A/en
Application granted granted Critical
Publication of CN107947246B publication Critical patent/CN107947246B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48Controlling the sharing of the in-phase component
    • H02J3/386
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Wind Motors (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

The invention discloses a wind power generation index distribution and power increase evaluation method considering frequency modulation power increase, and belongs to the technical field of active power control of power systems. The method comprises the steps of firstly, calculating the maximum available output and the actual output of each provincial dispatching wind power virtual machine set through related information in an automatic power generation control information system; the regional dispatching center calculates a wind power regulation instruction of the whole network according to the wind power consumption space of the whole network and the frequency modulation instruction condition; the regional dispatching center distributes wind power regulation instructions of each provincial dispatching center according to the maximum available output of each provincial wind power regulation virtual unit and by considering the historical output characteristics of each provincial wind power regulation virtual unit; and evaluating the peak-load power increasing amount and the frequency modulation power increasing amount of each provincial wind power modulation virtual machine set according to the regulation instruction data and the measurement data. According to the method, the wind turbine generator is introduced to participate in power grid frequency modulation so as to reduce wind abandon and improve system regulation performance, fairness and punishment mechanisms of regulation instruction distribution are considered, and reliability of participation of wind turbine generator in power grid regulation is improved.

Description

Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power
Technical Field
The invention belongs to the technical field of active power control of a power system, and particularly relates to a wind power generation index distribution and power increase evaluation method considering frequency modulation and power increase.
Background
As a renewable energy source with mature technology and outstanding economic benefits, wind power generation develops rapidly over the past two decades and is gradually becoming one of the main sources of electric energy supply. With the continuous improvement of the installation and grid-connected scale of the wind power generation, the influence of randomness and intermittence on the safe and stable operation of the power system is gradually shown: for the small-scale scheduling of the power system, the intermittency and uncertainty of wind power bring difficulties in the aspects of unit combination, tie line planning, climbing standby and the like; for the scheduling of the minute level and the second level, the fluctuation of the wind power brings difficulties in the aspects of system frequency control, tie line fluctuation adjustment and the like.
With the continuous development of wind power generation control technology, wind power plants gradually have the capability of receiving power system scheduling and participating in power system scheduling operation, and it has become possible to bring wind power generation into a power grid for operation as one of controlled resources. At present, the main control means of grid-connected wind power in China still adopts a wind power plant to track a day-ahead plan, which means that the actual operation of the wind power plant is still open-loop operation. However, the formulation of the day-ahead plan depends on wind power output prediction and load prediction, and prediction errors of the day-ahead plan have great influence on the actual operation of a wind power plant and a power system, which may cause poor utilization effect of wind power resources and exhaustion of the regulated reserve capacity of the power system, and reduce the economy and safety of power grid operation.
In order to solve the problem, the patent 'a wind power generation increasing control method with two-stage coordination of split-dispatching-provincial-dispatching' (patent application number: CN201610626975.4) discloses a wind power closed-loop control method, and the closed-loop operation of wind power generation participating in power system peak regulation is realized through the coordination of split-dispatching-provincial-dispatching and the calculation and distribution of wind power generation indexes. Furthermore, under the condition of wind limitation caused by difficult peak regulation, if the available output of wind power and the corresponding section constraint are still surplus, a proper control strategy can be adopted to enable the wind power generation unit to participate in the automatic power generation control of the power system, and the wind power generation increase is realized while the frequency modulation performance of the system is improved. However, the method does not consider the historical output characteristic and the regulation performance of the wind turbine generator, and the regulation performance is not rewarded or punished.
The patent 'a coordinated control method for wind generation sets to participate in automatic power generation control of an electric power system' (patent application number: CN201610248302.X) discloses a wind power frequency modulation and power increase control method under a wind limit condition. However, the calculation method of the method based on model predictive control in the regional dispatching center is based on the optimization thought, and the calculation real-time performance and the engineering availability are not high; when the instruction allocation is adjusted, the fairness and the reward and punishment mechanism of the allocation are not specifically considered.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a wind power generation index distribution and generation increase evaluation method considering frequency modulation and generation increase. The method can realize reasonable calculation, distribution and post evaluation of the wind power generation indexes by combining the running state of the wind power generation unit, the regional control constraint and the historical performance of the wind power generation unit.
The invention provides a wind power generation index distribution method considering frequency modulation and power increase, which is characterized by comprising the following steps of:
1) the maximum available output and the actual output of all wind power plants in the whole network and section constraint information are collected through an automatic generation control information system AGC, and the maximum available output and the actual output of each provincial dispatching wind power virtual machine set are calculated; the calculation formula is as follows:
Figure BDA0001468435580000023
wherein,
Figure BDA0001468435580000024
in order to save and adjust the maximum available output of the wind turbine virtual unit j,adjusting the actual output of the wind turbine virtual unit j for provincial wind turbine virtual unit number;PA,iThe maximum available output of the wind power plant i; pW,iThe actual output of the wind power plant i;and SkRespectively representing the maximum available wind power output under a section k and the capacity constraint of the section k, wherein k is 1,2, the. N is a radical ofW,kAnd NG,kRespectively the total number of the wind power plant under the section k and the total number of the conventional power plant under the section k; pG,iReal-time output of a conventional unit i; l is the total number of provincial dispatching wind power virtual units;
2) the regional dispatching center calculates a wind power regulation instruction of the whole network according to the wind power consumption space of the whole network and the frequency modulation instruction condition; the method comprises the following specific steps:
2-1) calculating the whole network wind power generation index
Figure BDA0001468435580000027
2-2) judging whether wind power in the area participates in frequency modulation; the method comprises the following specific steps:
2-2-1) judging the maximum available output of the wind power of the whole networkWhether the power generation index is larger than the full-grid wind power generation index
Figure BDA0001468435580000029
If it is
Figure BDA00014684355800000210
The wind power output in the area is limited, and the step 2-2-2) is carried out; otherwise, the wind-electricity output in the region is not limited, the wind-electricity in the region does not participate in frequency modulation, and the method is ended;
2-2-2) judging whether the wind power in the area is in a climbing stage: if the wind power in the region is in the output climbing stage, the wind power in the region does not participate in frequency modulation, and the method is ended; if the wind power in the region is not in the output climbing stage, the wind power in the region participates in frequency modulation, and the step 2-3 is carried out;
2-3) calculating a full-network wind power regulation instruction according to the full-network wind power consumption space and the frequency modulation instruction condition, wherein the calculation formula is as follows:
Figure BDA0001468435580000031
wherein,
Figure BDA0001468435580000032
for regulating the wind power in the whole network, PARRThe method comprises the steps that a full-network frequency modulation command is adopted, and alpha is a sharing coefficient of wind power participating in frequency modulation;
3) the regional dispatching center distributes wind power regulation instructions of each provincial dispatching center according to the maximum available output of each provincial wind power regulation virtual machine set and by considering the historical output characteristics of each provincial wind power regulation virtual machine set; the method comprises the following specific steps:
3-1) calculating a sharing weight coefficient of the adjusting power distribution of each provincial wind power adjusting virtual machine set;
new energy ultrashort-term prediction precision c of provincial dispatching wind power virtual unit jp,jThe definition is as follows:
Figure BDA0001468435580000033
in the formula, Ppredict,jUltra-short-term predicted output, P, of wind turbine virtual unit j for provincial regulation of non-wind-limiting timeactual,jAdjusting the actual output of the wind power virtual unit j for the non-wind-limiting time period province, and D is the predicted time point number of the non-wind-limiting time period in the day, ICjAdjusting the installed capacity of the wind power virtual machine set j for province;
adjusting performance score s of provincial wind power virtual unit jp,jSelecting an adjusting performance evaluation index of a provincial wind power virtual unit j controlled by a region;
respectively calculating the ultra-short-term prediction precision reference value of the provincial dispatching wind turbine virtual unit j
Figure BDA0001468435580000034
And adjusting the performance score benchmark value
Figure BDA0001468435580000035
The expression is as follows:
Figure BDA0001468435580000036
wherein,
Figure BDA0001468435580000038
is cp,jThe average of the degree of the week or month,
Figure BDA0001468435580000039
is s isp,jA weekly or monthly average of;
sharing weight coefficient W for adjusting power distribution of provincial wind power virtual unit jjThe calculation expression is as follows:
Figure BDA00014684355800000310
wherein, a and b are a prediction precision index weight coefficient and an adjustment performance index weight coefficient respectively;
3-2) distributing wind power adjusting instructions of the wind power virtual units of each provincial dispatching center;
adjusting instruction I of provincial wind power virtual unit jW,jThe calculation formula is as follows:
Figure BDA0001468435580000041
the invention provides a wind power generation index increase evaluation method considering frequency modulation increase, which is characterized by comprising the following steps of:
1) the maximum available output and the actual output of all wind power plants in the whole network and section constraint information are collected through an automatic generation control information system AGC, and the maximum available output and the actual output of each provincial dispatching wind power virtual machine set are calculated; the calculation formula is as follows:
Figure BDA0001468435580000042
Figure BDA0001468435580000043
Figure BDA0001468435580000044
wherein,
Figure BDA0001468435580000045
in order to save and adjust the maximum available output of the wind turbine virtual unit j,
Figure BDA0001468435580000046
the actual output of the provincial dispatching wind power virtual machine set j is obtained, and j is the serial number of the provincial dispatching wind power virtual machine set; pA,iThe maximum available output of the wind power plant i; pW,iThe actual output of the wind power plant i;
Figure BDA0001468435580000047
and SkRespectively representing the maximum available wind power output under a section k and the capacity constraint of the section k, wherein k is 1,2, the. N is a radical ofW,kAnd NG,kRespectively the total number of the wind power plant under the section k and the total number of the conventional power plant under the section k; pG,iReal-time output of a conventional unit i; l is the total number of provincial dispatching wind power virtual units;
2) the regional dispatching center calculates a wind power regulation instruction of the whole network according to the wind power consumption space of the whole network and the frequency modulation instruction condition; the method comprises the following specific steps:
2-1) calculating the whole network wind power generation index
Figure BDA0001468435580000048
2-2) judging whether wind power in the area participates in frequency modulation; the method comprises the following specific steps:
2-2-1) judging the maximum available output of the wind power of the whole network
Figure BDA0001468435580000049
Whether the power generation index is larger than the full-grid wind power generation index
Figure BDA00014684355800000410
If it is
Figure BDA00014684355800000411
The wind power output in the region is limited, and the step 2-2-2) is carried out; otherwise, the wind-electricity output in the region is not limited, the wind-electricity in the region does not participate in frequency modulation, and the method is ended;
2-2-2) judging whether the wind power in the area is in a climbing stage: if the wind power in the region is in the output climbing stage, the wind power in the region does not participate in frequency modulation, and the method is ended; if the wind power in the region is not in the output climbing stage, the wind power in the region participates in frequency modulation, and the step 2-3 is carried out;
2-3) calculating a full-network wind power regulation instruction according to the full-network wind power consumption space and the frequency modulation instruction condition, wherein the calculation formula is as follows:
Figure BDA0001468435580000051
wherein,
Figure BDA0001468435580000052
for regulating the wind power in the whole network, PARRThe method comprises the steps that a full-network frequency modulation command is adopted, and alpha is a sharing coefficient of wind power participating in frequency modulation;
3) the regional dispatching center distributes wind power regulation instructions of each provincial dispatching center according to the maximum available output of each provincial wind power regulation virtual machine set and by considering the historical output characteristics of each provincial wind power regulation virtual machine set; the method comprises the following specific steps:
3-1) calculating a sharing weight coefficient of the adjusting power distribution of each provincial wind power adjusting virtual machine set;
new energy ultrashort-term prediction precision c of provincial dispatching wind power virtual unit jp,jThe definition is as follows:
Figure BDA0001468435580000053
in the formula, Ppredict,jUltra-short-term predicted output, P, of wind turbine virtual unit j for provincial regulation of non-wind-limiting timeactual,jAdjusting the actual output of the wind power virtual unit j for the non-wind-limiting time period province, and D is the predicted time point number of the non-wind-limiting time period in the day, ICjAdjusting the installed capacity of the wind power virtual machine set j for province;
adjusting performance score s of provincial wind power virtual unit jp,jSelecting an adjusting performance evaluation index of a provincial wind power virtual unit j controlled by a region;
respectively calculating the ultra-short-term prediction precision reference value of the provincial dispatching wind turbine virtual unit j
Figure BDA0001468435580000054
And adjusting the performance score benchmark value
Figure BDA0001468435580000055
The expression is as follows:
Figure BDA0001468435580000057
wherein,
Figure BDA0001468435580000058
is cp,jThe average of the degree of the week or month,
Figure BDA0001468435580000059
is s isp,jA weekly or monthly average of;
sharing weight coefficient W for adjusting power distribution of provincial wind power virtual unit jjThe calculation expression is as follows:
Figure BDA00014684355800000510
wherein, a and b are a prediction precision index weight coefficient and an adjustment performance index weight coefficient respectively;
3-2) distributing wind power adjusting instructions of the wind power virtual units of each provincial dispatching center;
adjusting instruction I of provincial wind power virtual unit jW,jThe calculation formula is as follows:
Figure BDA0001468435580000061
4) evaluating peak-load capacity and frequency-modulation capacity of each provincial wind power modulation virtual machine set according to the regulation instruction data and the measurement data; the method comprises the following specific steps:
4-1) matching the adjusting instruction and the actual output of the provincial dispatching wind power virtual unit in time sequence to obtain the corresponding response delay time when the correlation degree of the actual output of the provincial dispatching wind power virtual unit and the adjusting instruction of the provincial dispatching wind power virtual unit is maximum;
setting the whole continuous regulation and control time period as a time window, and transversely translating the actual output of the provincial wind power virtual unit j along a time axis under the time window to obtain the translated actual output of the provincial wind power virtual unit j and the response delay time corresponding to the maximum correlation degree of the regulation instruction of the provincial wind power virtual unit j; wherein the maximum degree of correlation rjThe calculation formula is as follows:
Figure BDA0001468435580000062
wherein, deltajFor response delay time corresponding to the maximum correlation degree of the provincial dispatching wind power virtual unit j, R (x, y) is a function for solving the correlation degree of the vector x and the vector y; i isW,j(t) is an adjusting instruction of the provincial wind power virtual unit j at the moment t,adjusting the actual output of the wind turbine virtual unit j for the province at the moment t;
4-2) calculating the power increasing quantity of the provincial wind power and the wind power virtual machine group participating in peak shaving and the power increasing quantity of frequency modulation;
will economize on the power of wind turbine generator system j actual outputSimulated artificial wind limiting instruction P of provincial dispatching wind power virtual unit jC,jResponse delay time delta calculated according to step 4-1)jTranslating to respectively obtain the actual output of the translated provincial dispatching wind turbine virtual unit j
Figure BDA0001468435580000065
And the translated simulation artificial wind limiting instruction of the provincial dispatching wind turbine virtual unit j
Figure BDA0001468435580000066
Peak regulation and power generation increasing quantity of provincial wind power generation virtual machine set j
Figure BDA0001468435580000067
Sum frequency modulation power generation quantity increase
Figure BDA0001468435580000068
Respectively as follows:
Figure BDA0001468435580000069
where Δ t is an AGC command cycle duration.
The invention has the characteristics and beneficial effects that:
the method considers that wind power participates in frequency modulation control under the condition of wind limitation to realize power generation increase, distributes wind power generation indexes by combining the running state of the wind power generation set, regional control constraint and historical performance of the wind power generation set, and reasonably evaluates the power generation increase amount in two stages of peak regulation and frequency modulation. By the aid of wind power participating in frequency modulation, wind power generation index distribution considering a reward and punishment mechanism can be realized while the amount of abandoned wind power is reduced, and the wind power regulation performance is promoted.
Drawings
Fig. 1 is a flow chart of a wind power generation index distribution and generation increase evaluation method considering frequency modulation and generation increase according to the present invention.
Detailed Description
The invention provides a wind power generation index distribution and generation increase evaluation method considering frequency modulation and generation increase, which is further described in detail below by combining the accompanying drawings and specific embodiments.
The invention provides a wind power generation index distribution and generation increase evaluation method considering frequency modulation and generation increase, the whole flow is shown as figure 1, and the method comprises the following steps:
1) the maximum available output and the actual output of all wind power plants in the whole network and section constraint information from an automatic generation control information (AGC) system are collected by the AGC system, and the maximum available output and the actual output of each provincial dispatching wind power virtual unit are calculated; the calculation formula is as follows:
Figure BDA0001468435580000071
Figure BDA0001468435580000073
wherein,
Figure BDA0001468435580000074
in order to save and adjust the maximum available output of the wind turbine virtual unit j,the actual output of the provincial dispatching wind power virtual machine set j is obtained, and j is the serial number of the provincial dispatching wind power virtual machine set; pA,iIs the maximum available output (in MW) of the wind farm i; pW,iIs the actual contribution (in MW) of the wind farm i;
Figure BDA0001468435580000076
and SkThe maximum available power of the wind power under a section k and the capacity constraint of the section k are respectively set, wherein k is 1,2, the section number, M and k are the total number of the sections (province regulations may govern a plurality of wind power output bases, and each wind power output base is generally output through one section); n is a radical ofW,kAnd NG,kRespectively the total number of the wind power plant under the section k and the total number of the conventional power plant under the section k; pG,iReal-time output of a conventional unit i; and L is the total number of the provincial dispatching wind power virtual units. In the invention, the whole network refers to all power networks administered by a regional dispatching center and consists of a plurality of control areas; provincial dispatching is the mechanism responsible for the dispatching of each control area.
2) And the regional dispatching center calculates the wind power regulation instruction of the whole network according to the wind power consumption space of the whole network and the frequency modulation instruction condition. The method comprises the following specific steps:
2-1) calculating and considering peak regulation constraint full-grid wind power generation index
Figure BDA0001468435580000077
The calculation method of the index (in the prior art) should be selected according to the scheduling mechanism and the regulation and control resources of the regional scheduling center, and the main considered factors include: the method comprises the steps of regional load prediction, regional tie line planning, minimum adjustability of regional overall thermal power, reserved lower rotation standby, pumped storage power station pumped power generation power and the like.
Specifically, in the present embodiment, the following formula is used for calculation
Figure BDA0001468435580000081
Figure BDA0001468435580000082
Wherein, VfFor ultra-short term load prediction (ultra-short term load prediction generally refers to prediction of load change condition of 5 to 1 hour in the future), PtieScheduling plans for links, Vh-regRegulating the minimum output of electricity, V, for a regional systemr-regIn order to reserve the lower part for rotation for standby,Vr-plantfor self-contained power plant output, VpumpAnd pumping water to generate power for the pumped storage power station.
2-2) judging whether the wind power in the region participates in frequency modulation, wherein the wind power in the region consists of all provincial and tonal wind power virtual units in the whole network; the method comprises the following specific steps:
2-2-1) judging whether the wind power output in the area is limited, namely judging the maximum available wind power output of the whole networkWhether the power generation index is larger than the full-grid wind power generation indexIf it is
Figure BDA0001468435580000085
The wind power output in the region is limited, and the step 2-2-2) is carried out; otherwise, the wind-electricity output in the region is not limited, the wind-electricity in the region does not participate in frequency modulation, and the method is ended;
2-2-2) judging whether the wind power in the area is in a climbing stage: if the wind power in the region is in the output climbing stage, the wind power in the region does not participate in frequency modulation, and the method is ended; and if the wind power in the region is not in the output climbing stage, the wind power in the region can participate in frequency modulation, and the step 2-3) is carried out.
In particular, the embodiment adopts the full-grid wind power generation indexFluctuation ratio VR ofWAs a basis for judging whether the wind power is in the stage of climbing by exerting force in the region, VRWThe following formula is used for calculation:
Figure BDA0001468435580000087
wherein,and (4) predicting the power generation index of the whole network at the next time point (namely after 5 minutes). If VRWIf the wind power is less than or equal to 5%, the wind power in the region is not in the output climbing stage; otherwise, the wind power in the region is considered to be in the output climbing stage.
2-3) calculating a wind power regulation instruction of the whole network according to the wind power consumption space of the whole network and the frequency modulation instruction condition (all AGC systems have the calculation function of calculating the instruction), wherein the calculation formula is as follows:
Figure BDA0001468435580000089
wherein,
Figure BDA00014684355800000810
for regulating the wind power in the whole networkARRThe command is a whole network frequency modulation command (unit: MW), and alpha is a sharing coefficient (the value range is 0-1, and the value of the embodiment is 0.6) of wind power participating in frequency modulation. The formula shows that the adjustment instruction of the wind power of the whole network is to superpose a frequency modulation signal on the existing total wind power output; meanwhile, the condition that the adjustment instruction after superposition is larger than the full-network wind power generation index needs to be ensured
Figure BDA00014684355800000811
The maximum available output of the wind power of the whole network is smaller than the maximum available output of the wind power of the whole network, so that the regulation performance is ensured.
3) And the regional dispatching center distributes the wind power regulation instruction of each provincial dispatching center according to the maximum available output of each provincial wind power regulation virtual machine set and by considering the historical output characteristics of each provincial wind power regulation virtual machine set.
In order to ensure fairness in the new energy power generation index distribution process and consider concrete performances of each control area in the links of day-ahead prediction, real-time response and the like, evaluation indexes of new energy ultra-short-term prediction and real-time response can be introduced, the sharing weight of the adjustment power distribution of each provincial wind power regulation virtual machine set is obtained by combining the maximum available output of the wind power virtual machine set, and then wind power adjustment instructions of each provincial level scheduling center are distributed. The method comprises the following specific steps:
and 3-1) calculating a sharing weight coefficient of the adjusting power distribution of each province wind power adjusting virtual machine set. The weight coefficient is obtained by correcting two indexes of ultra-short prediction precision and adjustment performance score of new energy based on the maximum available output of each provincial wind power adjusting virtual machine set.
New energy ultrashort-term prediction precision c of provincial dispatching wind power virtual unit jp,jThe definition is as follows:
in the formula, Ppredict,jUltra-short-term predicted output, P, of wind turbine virtual unit j for provincial regulation of non-wind-limiting timeactual,jAdjusting the actual output of the wind power virtual unit j for the non-wind-limited time-period province, D is the number (generally 288 points) of predicted time points of the non-wind-limited time-period province in the day, and ICjAnd adjusting the installed capacity of the wind power virtual machine set j for the province.
Adjusting performance score s of provincial wind power virtual unit jp,jUsually, selecting an adjusting performance evaluation index of a provincial wind power virtual unit j controlled by a region; particularly, reference can be made to relevant indexes of AGC examination compensation, such as various regulation indexes specified by 'two detailed rules';
in order to ensure fairness and real-time performance, the two indexes can be calculated in a rolling mode, and the average value of the period or month of the two indexes is selected
Figure BDA0001468435580000092
For subsequent calculations. In order to more fairly take the two indexes into account in the distribution process of the new energy power generation indexes and correct the two indexes, each index can be divided by the average value of the indexes of each control area to serve as a calculation index, and the calculation result can be respectively expressed as an ultra-short-term prediction precision reference value of the provincial dispatching wind power virtual unit j
Figure BDA0001468435580000093
And adjusting the performance score benchmark valueThe expression is as follows:
Figure BDA0001468435580000096
the calculation process is to perform per unit on the two indexes respectively, so that the influence of the two indexes on the distribution result is closer, and the selection of the subsequent weight coefficient is facilitated.
Sharing weight coefficient W for adjusting power distribution of provincial wind power virtual unit jjThe calculation expression is as follows:
Figure BDA0001468435580000101
and a and b are respectively a prediction precision index weight coefficient and an adjustment performance index weight coefficient, and the value ranges are 0-1. In particular, in this embodiment, a is 0.5.
And 3-2) distributing wind power adjusting instructions of the wind power virtual units of each provincial dispatching center.
Adjusting instruction I of provincial wind power virtual unit jW,jThe calculation formula is as follows:
Figure BDA0001468435580000102
4) and evaluating the peak-load power increasing amount and the frequency modulation power increasing amount of each provincial wind power modulation virtual machine set according to the regulation instruction data and the measurement data. The method comprises the following specific steps:
4-1) matching the adjusting instruction and the actual output of the provincial dispatching wind power virtual unit in time sequence to obtain the corresponding response delay time when the correlation degree of the actual output of the provincial dispatching wind power virtual unit and the adjusting instruction of the provincial dispatching wind power virtual unit is maximum; considering that the whole continuous regulation and control time period is a time window (continuous regulation and control comprises only peak regulation and frequency modulation), the regulation instruction of the provincial wind power regulation virtual machine set issued in the time window is executed by wind power after a certain fixed delay. To determine the fixed delay, the actual output of provincial wind turbine generator set j may be advanced along the time axisAnd horizontally translating the line to obtain the actual output of the provincial dispatching wind turbine virtual unit j after translation and the response delay time corresponding to the maximum correlation degree of the regulating instruction of the provincial dispatching wind turbine virtual unit j. Wherein the maximum degree of correlation rjThe calculation formula is as follows:
Figure BDA0001468435580000103
wherein, deltajAnd R (x, y) is a function for solving the correlation degree of the vector x and the vector y in order to adjust the response delay time corresponding to the maximum correlation degree of the wind turbine virtual unit j. I isW,j(t) is an adjusting instruction of the provincial wind power virtual unit j at the moment t,
Figure BDA0001468435580000104
adjusting the actual output of the wind power virtual machine set for each province at the moment t;
and 4-2) calculating the power increasing quantity of the provincial wind power and the wind power virtual set participating in peak shaving and the power increasing quantity of the frequency modulation.
Will economize on the power of wind turbine generator system j actual output
Figure BDA0001468435580000105
Simulated artificial wind limiting instruction P of provincial dispatching wind power virtual unit jC,j(PC,jThe wind power generation index with the lowest wind limit time period can be selected) according to the calculated response delay time deltajTranslating to respectively obtain the actual output of the translated provincial dispatching wind turbine virtual unit jAnd the translated simulation artificial wind limiting instruction of the provincial dispatching wind turbine virtual unit j
Figure BDA0001468435580000107
And the frequency modulation power increasing quantity is calculated.
Peak regulation and power generation increasing quantity of provincial wind power generation virtual machine set jSum frequency modulation power generation quantity increaseRespectively as follows:
Figure BDA00014684355800001010
Figure BDA0001468435580000111
wherein,
Figure BDA0001468435580000112
andthe actual output of the translated provincial wind power virtual unit j and the simulated artificial wind limiting instruction of the translated provincial wind power virtual unit j are respectively, and delta t is the period duration of one AGC instruction cycle.
By the aid of the increased power generation evaluation, peak regulation and frequency modulation increased power generation quantity distribution and splitting calculation can be achieved to serve as a basis for compensation of the wind turbine generator.

Claims (6)

1. A wind power generation index distribution method considering frequency modulation and power increase is characterized by comprising the following steps:
1) the maximum available output and the actual output of all wind power plants in the whole network and section constraint information are collected through an automatic generation control information system AGC, and the maximum available output and the actual output of each provincial dispatching wind power virtual machine set are calculated; the calculation formula is as follows:
Figure FDA0002194351090000011
Figure FDA0002194351090000012
wherein,
Figure FDA0002194351090000014
in order to save and adjust the maximum available output of the wind turbine virtual unit j,
Figure FDA0002194351090000015
the actual output of the provincial dispatching wind power virtual machine set j is obtained, and j is the serial number of the provincial dispatching wind power virtual machine set; pA,iThe maximum available output of the wind power plant i; pW,iThe actual output of the wind power plant i;and SkRespectively representing the maximum available wind power output under a section k and the capacity constraint of the section k, wherein k is 1,2, the. N is a radical ofW,kAnd NG,kRespectively the total number of the wind power plant under the section k and the total number of the conventional power plant under the section k; pG,iReal-time output of a conventional unit i; l is the total number of provincial dispatching wind power virtual units;
2) the regional dispatching center calculates a wind power regulation instruction of the whole network according to the wind power consumption space of the whole network and the frequency modulation instruction condition; the method comprises the following specific steps:
2-1) calculating the whole network wind power generation index
Figure FDA0002194351090000017
2-2) judging whether wind power in the area participates in frequency modulation; the method comprises the following specific steps:
2-2-1) judging the maximum available output of the wind power of the whole network
Figure FDA0002194351090000018
Whether the power generation index is larger than the full-grid wind power generation index
Figure FDA0002194351090000019
If it is
Figure FDA00021943510900000110
The wind power output in the region is limited, and the step 2-2-2) is carried out; otherwise, the wind-electricity output in the region is not limited, the wind-electricity in the region does not participate in frequency modulation, and the method is ended;
2-2-2) judging whether the wind power in the area is in a climbing stage: if the wind power in the region is in the output climbing stage, the wind power in the region does not participate in frequency modulation, and the method is ended; if the wind power in the region is not in the output climbing stage, the wind power in the region participates in frequency modulation, and the step 2-3 is carried out;
2-3) calculating a full-network wind power regulation instruction according to the full-network wind power consumption space and the frequency modulation instruction condition, wherein the calculation formula is as follows:
wherein,for regulating the wind power in the whole network, PARRThe method comprises the steps that a full-network frequency modulation command is adopted, and alpha is a sharing coefficient of wind power participating in frequency modulation;
3) the regional dispatching center distributes wind power regulation instructions of each provincial dispatching center according to the maximum available output of each provincial wind power regulation virtual machine set and by considering the historical output characteristics of each provincial wind power regulation virtual machine set; the method comprises the following specific steps:
3-1) calculating a sharing weight coefficient of the adjusting power distribution of each provincial wind power adjusting virtual machine set;
new energy ultrashort-term prediction precision c of provincial dispatching wind power virtual unit jp,jThe definition is as follows:
Figure FDA0002194351090000023
in the formula, Ppredict,jUltra-short-term predicted output, P, of wind turbine virtual unit j for provincial regulation of non-wind-limiting timeactual,jAdjusting the actual output of the wind power virtual unit j for the provincial wind-free time period, D is the number of predicted time points for the non-wind-limited time period in the day,ICjadjusting the installed capacity of the wind power virtual machine set j for province;
calculating the regulation performance score s of the provincial dispatching wind turbine virtual unit j according to the regulation performance evaluation index of the provincial dispatching wind turbine virtual unit j controlled by the regionp,j
Respectively calculating the ultra-short-term prediction precision reference value of the provincial dispatching wind turbine virtual unit jAnd adjusting the performance score benchmark value
Figure FDA0002194351090000025
The expression is as follows:
Figure FDA0002194351090000026
wherein,is cp,jThe average of the degree of the week or month,
Figure FDA0002194351090000029
is s isp,jA weekly or monthly average of;
sharing weight coefficient W for adjusting power distribution of provincial wind power virtual unit jjThe calculation expression is as follows:
Figure FDA0002194351090000031
wherein, a and b are a prediction precision index weight coefficient and an adjustment performance index weight coefficient respectively;
3-2) distributing wind power adjusting instructions of the wind power virtual units of each provincial dispatching center;
province regulationAdjusting instruction I of wind power virtual unit jW,jThe calculation formula is as follows:
2. the method of claim 1, wherein the full grid wind power generation indicator in step 2-1)
Figure FDA0002194351090000033
The calculation method is as follows:
Figure FDA0002194351090000034
wherein, VfFor ultra-short term load prediction, PtieScheduling plans for links, Vh-regRegulating the minimum output of electricity, V, for a regional systemr-regFor reserve rotation, Vr-plantFor self-contained power plant output, VpumpAnd pumping water to generate power for the pumped storage power station.
3. The method according to claim 1, wherein in step 2-2-2), whether the wind power in the area is in the output climbing stage is determined by the following specific method:
calculating the whole network wind power generation index
Figure FDA0002194351090000035
The expression of fluctuation ratio of (a) is as follows:
wherein,
Figure FDA0002194351090000037
the predicted value of the power generation index of the whole network at the next time point is obtained;
if VRWIf the wind power is less than or equal to 5%, the wind power in the region is not in the stage of output climbing; otherwise, the wind power in the region is in the output climbing stage.
4. A wind power generation index increase evaluation method considering frequency modulation increase is characterized by comprising the following steps:
1) the maximum available output and the actual output of all wind power plants in the whole network and section constraint information are collected through an automatic generation control information system AGC, and the maximum available output and the actual output of each provincial dispatching wind power virtual machine set are calculated; the calculation formula is as follows:
Figure FDA0002194351090000038
Figure FDA0002194351090000039
wherein,
Figure FDA0002194351090000042
in order to save and adjust the maximum available output of the wind turbine virtual unit j,
Figure FDA0002194351090000043
the actual output of the provincial dispatching wind power virtual machine set j is obtained, and j is the serial number of the provincial dispatching wind power virtual machine set; pA,iThe maximum available output of the wind power plant i; pW,iThe actual output of the wind power plant i;
Figure FDA0002194351090000044
and SkRespectively representing the maximum available wind power output under a section k and the capacity constraint of the section k, wherein k is 1,2, the. N is a radical ofW,kAnd NG,kTotal number of wind power plants under section k and conventional power plants under section k respectivelyTotal number; pG,iReal-time output of a conventional unit i; l is the total number of provincial dispatching wind power virtual units;
2) the regional dispatching center calculates a wind power regulation instruction of the whole network according to the wind power consumption space of the whole network and the frequency modulation instruction condition; the method comprises the following specific steps:
2-1) calculating the whole network wind power generation index
Figure FDA0002194351090000045
2-2) judging whether wind power in the area participates in frequency modulation; the method comprises the following specific steps:
2-2-1) judging the maximum available output of the wind power of the whole network
Figure FDA0002194351090000046
Whether the power generation index is larger than the full-grid wind power generation indexIf it is
Figure FDA0002194351090000048
The wind power output in the region is limited, and the step 2-2-2) is carried out; otherwise, the wind-electricity output in the region is not limited, the wind-electricity in the region does not participate in frequency modulation, and the method is ended;
2-2-2) judging whether the wind power in the area is in a climbing stage: if the wind power in the region is in the output climbing stage, the wind power in the region does not participate in frequency modulation, and the method is ended; if the wind power in the region is not in the output climbing stage, the wind power in the region participates in frequency modulation, and the step 2-3 is carried out;
2-3) calculating a full-network wind power regulation instruction according to the full-network wind power consumption space and the frequency modulation instruction condition, wherein the calculation formula is as follows:
wherein,for regulating the wind power in the whole network, PARRThe method comprises the steps that a full-network frequency modulation command is adopted, and alpha is a sharing coefficient of wind power participating in frequency modulation;
3) the regional dispatching center distributes wind power regulation instructions of each provincial dispatching center according to the maximum available output of each provincial wind power regulation virtual machine set and by considering the historical output characteristics of each provincial wind power regulation virtual machine set; the method comprises the following specific steps:
3-1) calculating a sharing weight coefficient of the adjusting power distribution of each provincial wind power adjusting virtual machine set;
new energy ultrashort-term prediction precision c of provincial dispatching wind power virtual unit jp,jThe definition is as follows:
Figure FDA0002194351090000051
in the formula, Ppredict,jUltra-short-term predicted output, P, of wind turbine virtual unit j for provincial regulation of non-wind-limiting timeactual,jAdjusting the actual output of the wind power virtual unit j for the non-wind-limiting time period province, and D is the predicted time point number of the non-wind-limiting time period in the day, ICjAdjusting the installed capacity of the wind power virtual machine set j for province;
calculating the regulation performance score s of the provincial dispatching wind turbine virtual unit j according to the regulation performance evaluation index of the provincial dispatching wind turbine virtual unit j controlled by the regionp,j
Respectively calculating the ultra-short-term prediction precision reference value of the provincial dispatching wind turbine virtual unit j
Figure FDA0002194351090000052
And adjusting the performance score benchmark value
Figure FDA0002194351090000053
The expression is as follows:
Figure FDA0002194351090000055
wherein,
Figure FDA0002194351090000056
is cp,jThe average of the degree of the week or month,
Figure FDA0002194351090000057
is s isp,jA weekly or monthly average of;
sharing weight coefficient W for adjusting power distribution of provincial wind power virtual unit jjThe calculation expression is as follows:
Figure FDA0002194351090000058
wherein, a and b are a prediction precision index weight coefficient and an adjustment performance index weight coefficient respectively;
3-2) distributing wind power adjusting instructions of the wind power virtual units of each provincial dispatching center;
adjusting instruction I of provincial wind power virtual unit jW,jThe calculation formula is as follows:
4) evaluating peak-load capacity and frequency-modulation capacity of each provincial wind power modulation virtual machine set according to the regulation instruction data and the measurement data; the method comprises the following specific steps:
4-1) matching the adjusting instruction and the actual output of the provincial dispatching wind power virtual unit in time sequence to obtain the corresponding response delay time when the correlation degree of the actual output of the provincial dispatching wind power virtual unit and the adjusting instruction of the provincial dispatching wind power virtual unit is maximum;
setting the whole continuous regulation and control time period as a time window, and transversely translating the actual output of the provincial wind power virtual unit j along a time axis under the time window to obtain the translated actual output of the provincial wind power virtual unit j and the response delay time corresponding to the maximum correlation degree of the regulation instruction of the provincial wind power virtual unit j; wherein the maximum phaseDegree of closeness rjThe calculation formula is as follows:
wherein, deltajFor response delay time corresponding to the maximum correlation degree of the provincial dispatching wind power virtual unit j, R (x, y) is a function for solving the correlation degree of the vector x and the vector y; i isW,j(t) is an adjusting instruction of the provincial wind power virtual unit j at the moment t,
Figure FDA0002194351090000062
adjusting the actual output of the wind turbine virtual unit j for the province at the moment t;
4-2) calculating the power increasing quantity of the provincial wind power and the wind power virtual machine group participating in peak shaving and the power increasing quantity of frequency modulation;
will economize on the power of wind turbine generator system j actual output
Figure FDA0002194351090000063
Simulated artificial wind limiting instruction P of provincial dispatching wind power virtual unit jC,jResponse delay time delta calculated according to step 4-1)jTranslating to respectively obtain the actual output of the translated provincial dispatching wind turbine virtual unit j
Figure FDA0002194351090000064
And the translated simulation artificial wind limiting instruction of the provincial dispatching wind turbine virtual unit j
Figure FDA0002194351090000065
Peak regulation and power generation increasing quantity of provincial wind power generation virtual machine set j
Figure FDA0002194351090000066
Sum frequency modulation power generation quantity increase
Figure FDA0002194351090000067
Respectively as follows:
Figure FDA0002194351090000068
Figure FDA0002194351090000069
where Δ t is an AGC command cycle duration.
5. The method of claim 4, wherein the full-grid wind power generation index in step 2-1)
Figure FDA00021943510900000610
The calculation method is as follows:
Figure FDA00021943510900000611
wherein, VfFor ultra-short term load prediction, PtieScheduling plans for links, Vh-regRegulating the minimum output of electricity, V, for a regional systemr-regFor reserve rotation, Vr-plantFor self-contained power plant output, VpumpAnd pumping water to generate power for the pumped storage power station.
6. The method according to claim 4, wherein the step 2-2-2) of determining whether the wind power in the area is in the output climbing stage is as follows:
calculating the whole network wind power generation indexThe expression of fluctuation ratio of (a) is as follows:
Figure FDA0002194351090000072
wherein,
Figure FDA0002194351090000073
the predicted value of the power generation index of the whole network at the next time point is obtained;
if VRWIf the wind power is less than or equal to 5%, the wind power in the region is not in the stage of output climbing; otherwise, the wind power in the region is in the output climbing stage.
CN201711126400.7A 2017-11-15 2017-11-15 Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power Active CN107947246B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711126400.7A CN107947246B (en) 2017-11-15 2017-11-15 Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711126400.7A CN107947246B (en) 2017-11-15 2017-11-15 Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power

Publications (2)

Publication Number Publication Date
CN107947246A CN107947246A (en) 2018-04-20
CN107947246B true CN107947246B (en) 2020-01-10

Family

ID=61932208

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711126400.7A Active CN107947246B (en) 2017-11-15 2017-11-15 Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power

Country Status (1)

Country Link
CN (1) CN107947246B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110429666B (en) * 2019-08-22 2021-08-27 国网冀北电力有限公司电力科学研究院 Power generation control method and system based on electricity abandoning balance
CN112821469B (en) * 2021-03-09 2022-10-25 中国南方电网有限责任公司 Day-ahead power generation scheduling optimization method and device based on frequency modulation absorption domain analysis
CN113346475B (en) * 2021-04-16 2022-09-02 国电南瑞科技股份有限公司 Configurable wind-solar absorption index distribution method and device based on newly added object

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104037805A (en) * 2014-06-19 2014-09-10 清华大学 Power generation margin distribution method for photovoltaic power station based on power grid security constraints
CN104600747A (en) * 2015-01-21 2015-05-06 西安交通大学 Operation optimizing method capable of coordinating operation risk and wind energy consumption of power system
CN105914779A (en) * 2016-04-20 2016-08-31 国家电网公司 Coordination control method for wind turbine generator set participating in automatic power generation control of power system
CN106160000A (en) * 2016-08-02 2016-11-23 国电南瑞科技股份有限公司 Adjust for a kind of point and save the wind-powered electricity generation additional issue control method adjusting two-stage to coordinate

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10230242B2 (en) * 2014-10-10 2019-03-12 Edge Electrons Limited Maximizing energy savings by utilizing conservation voltage reduction with adaptive voltage control and peak demand reduction at point of use
CN106471695B (en) * 2014-11-24 2019-05-14 Abb瑞士股份有限公司 Black starting-up blower, wind power plant and the method for restoring wind power plant and blower, and the blower using this method, wind power plant

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104037805A (en) * 2014-06-19 2014-09-10 清华大学 Power generation margin distribution method for photovoltaic power station based on power grid security constraints
CN104600747A (en) * 2015-01-21 2015-05-06 西安交通大学 Operation optimizing method capable of coordinating operation risk and wind energy consumption of power system
CN105914779A (en) * 2016-04-20 2016-08-31 国家电网公司 Coordination control method for wind turbine generator set participating in automatic power generation control of power system
CN106160000A (en) * 2016-08-02 2016-11-23 国电南瑞科技股份有限公司 Adjust for a kind of point and save the wind-powered electricity generation additional issue control method adjusting two-stage to coordinate

Also Published As

Publication number Publication date
CN107947246A (en) 2018-04-20

Similar Documents

Publication Publication Date Title
CN108092324B (en) AGC control system and control method for wind power participating in peak shaving frequency modulation
CN105046395B (en) Method for compiling day-by-day rolling plan of power system containing multiple types of new energy
CN110782281B (en) Day-ahead market clearing method for multi-owner cascade power station basin electric quantity transfer
CN109245183A (en) A kind of honourable permeability area power grid peak regulating method of height based on load control system
CN109636674B (en) Large-scale hydropower station group monthly transaction electric quantity decomposition and checking method
CN107248751A (en) A kind of energy storage station dispatch control method for realizing distribution network load power peak load shifting
CN108039737B (en) Source-grid-load coordinated operation simulation system
CN105006843A (en) Multi-time-scale flexible load scheduling method for handling wind power uncertainties
CN104915737A (en) Coordinating economic dispatching method for power system containing large-scale wind power based on risk control
CN105162149A (en) Fuzzy adaptive control based method for tracking output of power generation plan of light storage system
CN107947246B (en) Wind power generation index distribution and increased power evaluation method considering frequency modulation and increased power
CN109447405A (en) A kind of library multi-stag step library group's short-term plan formulating method undertaking peak regulation task
CN111832824A (en) Electric power frequency modulation market trading clearing and settlement method, device and system
CN113054669B (en) Block chain technology-based distribution network peak-shifting valley-leveling self-adaptive self-balancing method
CN114648195B (en) Intelligent cooperative regulation and control method for energy storage power station in multi-level auxiliary service of power grid
CN110829408A (en) Multi-domain scheduling method considering energy storage power system based on power generation cost constraint
CN112909933B (en) Intraday rolling optimization scheduling method containing pumped storage unit under spot market environment
CN110535119A (en) Adapt to the energy-accumulating power station period electricity scroll control method of peak load regulation network frequency modulation
CN117709682B (en) Reservoir water level disturbance interval assessment and control method under cascade hydropower response flexibility demand mode
CN109726894A (en) Ensure the new energy active command calculation method of spot exchange and medium-term and long-term electricity
CN105678394A (en) Multi-source and multi-cycle generation schedule formulation method
CN103248050A (en) Competition game-based cooperative wind power ramp control system and competition game-based cooperative wind power ramp control method
CN114498773A (en) Control method and device for participating in APC frequency modulation of master station by adjustable load
CN117767349A (en) Optimal scheduling and coordination control method and system for wind-solar new energy power station
Hong-Tao et al. Local Consumption model of wind and photovoltaic power based on demand side response

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