CN106202815A - Double-feed wind power field based on meritorious response Dynamic Equivalence - Google Patents
Double-feed wind power field based on meritorious response Dynamic Equivalence Download PDFInfo
- Publication number
- CN106202815A CN106202815A CN201610594669.7A CN201610594669A CN106202815A CN 106202815 A CN106202815 A CN 106202815A CN 201610594669 A CN201610594669 A CN 201610594669A CN 106202815 A CN106202815 A CN 106202815A
- Authority
- CN
- China
- Prior art keywords
- wind
- power
- equivalent
- meritorious
- equivalence
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- 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
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Wind Motors (AREA)
Abstract
Double-feed wind power field based on meritorious response Dynamic Equivalence, relates to electric system simulation modeling technique field.The present invention is cannot to take into account amount of calculation and equivalent precision to solve existing wind energy turbine set equivalence method, cannot complete the problem that blower fan hives off by off-line simultaneously.Build Wind turbines electrical-magnetic model and carry out three-phase shortcircuit emulation and on-the-spot test Experimental comparison, if precision conformance with standard, obtain Wind turbines transient state active power response curve in full blast speed operation area, obtain, according to its characteristic of bunching, index of hiving off；Setting meritorious power error threshold value according to equivalent required precision, the equivalent example carrying out typical wind speed distribution calculates, if precision is unsatisfactory for requirement, in adding corresponding wind speed section, unit equivalence wind speed is as supplementary cut-point；Obtaining, according to off-line, index of hiving off, be one by unit equivalence in each for actual wind energy turbine set group, the equivalent parameters calculating each equivalent unit and collection electric network obtains wind energy turbine set Equivalent Model.For setting up wind energy turbine set Equivalent Model.
Description
Technical field
The present invention relates to the Dynamic Equivalence of a kind of Double-feed wind power field, belong to electric system simulation modeling technique neck
Territory.
Background technology
Increasingly serious due to the energy and environmental problem, wind energy with its technology with the advantage becoming present aspect in power system
It is widely applied.But, along with the quick increase of wind energy permeability, operational reliability and the stability of power system are subject to
Certain impact and impact, therefore access power system extensive to wind-powered electricity generation is analyzed the accuracy requirement that calculates increasingly
High.One wind energy turbine set often has tens typhoon group of motors the most up to a hundred, each typhoon group of motors again by wind energy conversion system, power transmission shaft,
The module compositions such as electromotor, current transformer and controller thereof, protection device and controller thereof, if the most independent to every typhoon group of motors
Modeling, then can be greatly increased complexity and the calculating time of phantom.Therefore, the research of wind energy turbine set Equivalent Model is carried out very
Necessary.
In wind energy turbine set equivalent modeling, owing to the fitting precision that large-scale wind power field unit is equivalent is generally difficult to satisfied wanting
Asking, how carrying out Wind turbines rationally effectively hives off is the matter of utmost importance of wind energy turbine set equivalent modeling research needs solution.Mesh
Till before, the grouping method of wind energy turbine set mainly has a two categories below:
(1) grouping method of unit cluster is carried out according to the similarity of wind speed.The method is made according to wake effect and wind direction
The wind speed difference become, is divided into multiple region by wind energy turbine set, the most corresponding equivalent unit.Or direct basis wind speed is close
Property, hive off according to the working region of Wind turbines.This kind of grouping method only considered Index For Steady-state, when wind energy turbine set inner blower quantity
When numerous, operating condition differs greatly, systematic error can be bigger.
(2) so that the characteristic quantity grouping method as index of operating states of the units can be characterized.With Wind turbines linearisation
The influence factor of state dominant characteristics root or rotating speed, pitch angle etc. can characterize the variable of unit duty as the finger that hives off
Mark, although this kind of grouping method can reach higher equivalent precision, but it is generally required to uses more complicated algorithm, amount of calculation
Greatly.Operating mode difference can still result in the increase of equivalent wind turbine number of units time big, and when inputting wind speed change, it is impossible to divided in real time
Group's index, engineering uses has bigger limitation.
Summary of the invention
The present invention is cannot to take into account amount of calculation and equivalent precision to solve existing wind energy turbine set equivalence method, simultaneously cannot
Off-line completes the problem that blower fan hives off.Double-feed wind power field based on meritorious response Dynamic Equivalence is now provided.
Double-feed wind power field based on meritorious response Dynamic Equivalence, it comprises the following steps:
Step one, building the transient Model of Double-feed wind power unit, the small-power at this Wind turbines exports and big respectively
Under power two kinds of operating modes of output, carrying out three phase short circuit fault emulation, the supplemental characteristic this emulation obtained is tested with on-the-spot test
The supplemental characteristic contrast obtained, it is thus achieved that simulation result and the weighted average absolute error of field experimentation test result, this weighting is put down
All absolute error conformance with standard required precisions, then perform step 2；
In step 2, acquisition full blast speed operation area, Double-feed wind power unit transient Model occurs three-phase short at also site
Transient state active power response curve during the fault of road, according to the characteristic of bunching on this curve, it is thus achieved that index of hiving off；
Step 3, set meritorious power error threshold value according to actual equivalent required precision, carry out typical wind speed distribution etc.
Value example calculate, it is thus achieved that the meritorious relative systematic error of each example, by the meritorious relative systematic error of each example respectively with setting
Active power error threshold compares, and verifies cut-point, if the meritorious relative systematic error of certain example is more than being set with
Merit power error threshold value, then add the unit equivalence wind speed in corresponding wind speed section as supplementary cut-point to this example；
Step 4, the cut-point obtained according to step 2 and step 3 off-line, take the union of cut-point as index of hiving off
Divide Wind turbines, be one by the Wind turbines equivalence in each group, calculate the equivalent parameters of each equivalent unit and corresponding collection
Electric network equivalent parameters, it is thus achieved that wind energy turbine set Equivalent Model.
The invention have the benefit that the present invention by building the electromagnetic transient simulation model of Double-feed wind power unit,
Under the output of unit small-power and two kinds of operating modes of high-power output, three phase short circuit fault is set, Wind turbines exit is recorded
Active power, reactive power, voltage, the response curve of electric current and field test data contrast, and are weighted by each time period
Mean absolute error formula obtains simulation parameter and the error of fitting of field test data, this error required precision scope it
In, illustrate that phantom is consistent with the dynamic characteristic of actual LVRT Capability of Wind Turbine Generator；Then active power is selected to analyze
Its transient response characteristic, when blower fan is operated in different wind speed regions, it is the poorest that the transient response curve of its active power has
Different, according to the characteristic of bunching on this curve of active power, it can be seen that to work in promoter region, maximal power tracing district and permanent turn
The Wind turbines in speed invariable power district has characteristic of significantly bunching, and can identify 7m/s (starting point in MPPT district) and 12m/s
(terminal in MPPT district) two cut-points.According to the roll-off characteristic that active power curves is stronger, typical wind speed distribution scene is set
With meritorious relative systematic error threshold value, miss relative with the active power arranged respectively for the meritorious relative systematic error of each scene
Difference allows threshold ratio relatively, verifies the cut-point of each scene, if the meritorious relative systematic error of certain scene is more than arranging
Active power relative error allow threshold value, then this scene is added unit equivalence wind speed in corresponding wind speed section as supplementing point
Cutpoint；Finally, it is thus achieved that the cut-point of each scene, the union taking cut-point divides Wind turbines, by each group as index of hiving off
Wind turbines equivalence be one, calculate the equivalent parameters of each equivalent unit and corresponding current collection network equivalent parameter, it is thus achieved that wind
Electric field Equivalent Model, utilizes the wind energy turbine set Equivalent Model that step 4 draws, uses the standard of actual wind speed data verification Equivalent Model
Really property.
The method overcome current equivalence method to take into account amount of calculation and equivalent precision and wind-powered electricity generation cannot be completed by off-line
The problem of machine component group, can be obviously improved the accuracy that conventional individual is equivalent, and principle is simple, and physical meaning is clear and definite, it is not necessary to
Complicated calculating, has good adaptability to air speed data and position of failure point.
Accompanying drawing explanation
The phantom of the Double-feed wind power unit that Fig. 1 is built by the present invention；
Fig. 2 be Double-feed wind power unit under the operating mode that small-power exports, at fan outlet occur three-phase shortcircuit time having
Merit Simulation of SAR power image and the comparison diagram of measured result；
Fig. 3 be Double-feed wind power unit under the operating mode that small-power exports, at fan outlet occur three-phase shortcircuit time nothing
Merit Simulation of SAR power image and the comparison diagram of measured result；
Fig. 4 be Double-feed wind power unit under the operating mode that small-power exports, at fan outlet occur three-phase shortcircuit time electricity
Pressure emulation and the comparison diagram of measured result；
Fig. 5 be Double-feed wind power unit under the operating mode that small-power exports, at fan outlet occur three-phase shortcircuit time electricity
Stream emulation and the comparison diagram of measured result；
Fig. 6 be Double-feed wind power unit under the operating mode of high-power output, at fan outlet occur three-phase shortcircuit time having
Merit Simulation of SAR power image and the comparison diagram of measured result；
Fig. 7 be Double-feed wind power unit under the operating mode of high-power output, at fan outlet occur three-phase shortcircuit time nothing
Merit Simulation of SAR power image and the comparison diagram of measured result；
Fig. 8 be Double-feed wind power unit under the operating mode of high-power output, at fan outlet occur three-phase shortcircuit time electricity
Pressure emulation and the comparison diagram of measured result；
Fig. 9 be Double-feed wind power unit under the operating mode of high-power output, at fan outlet occur three-phase shortcircuit time electricity
Stream emulation and the comparison diagram of measured result；
Figure 10 is the definition schematic diagram of Error Calculation time period, and map parameter t represents the moment that fault occurs, parameter t0Generation
The time that table crowbar circuit puts into, parameter t1The representing fault persistent period, parameter t3-t1Representing fault recover needed for time
Between, T1For testing the parameter reciprocal of mains frequency, t+t3Represent the moment entering stable state；
Figure 11 is the active power response curve in the range of Wind turbines full blast speed, and reference 1 represents promoter region, accompanying drawing
Labelling 2 represents that maximal power tracing district, reference 3 represent permanent rotating speed invariable power district；
Figure 12 is the partial enlarged drawing at 70s to the 70.2s place in Figure 11；
Figure 13 is the air speed data of l-G simulation test；
When Figure 14 is the short circuit of wind energy turbine set exit, the 6th group of detailed model of air speed data, unit Equivalent Model and use this
The active power simulation result comparison diagram that inventive method obtains；
When Figure 15 is the short circuit of wind energy turbine set exit, the 6th group of detailed model of air speed data, unit Equivalent Model and use this
The reactive power simulation result comparison diagram that inventive method obtains；
When Figure 16 is the short circuit of wind energy turbine set exit, the 6th group of detailed model of air speed data, unit Equivalent Model and use this
The voltage simulation result comparison diagram that inventive method obtains；
When Figure 17 is the short circuit of wind energy turbine set exit, the 6th group of detailed model of air speed data, unit Equivalent Model and use this
The current simulations Comparative result figure that inventive method obtains；
When Figure 18 is the short circuit of wind energy turbine set exit, the 30th group of detailed model of air speed data, unit Equivalent Model and employing
The active power simulation result comparison diagram that the inventive method obtains；
When Figure 19 is the short circuit of wind energy turbine set exit, the 30th group of detailed model of air speed data, unit Equivalent Model and employing
The reactive power simulation result comparison diagram that the inventive method obtains；
When Figure 20 is the short circuit of wind energy turbine set exit, the 30th group of detailed model of air speed data, unit Equivalent Model and employing
The voltage simulation result comparison diagram that the inventive method obtains；
When Figure 21 is the short circuit of wind energy turbine set exit, the 30th group of detailed model of air speed data, unit Equivalent Model and employing
The current simulations Comparative result figure that the inventive method obtains；
When Figure 22 is to point out short circuit in supply lines one loop line, the 6th group of detailed model of air speed data, unit Equivalent Model
With the active power simulation result comparison diagram using equivalence method of the present invention to obtain；
When Figure 23 is to point out short circuit in supply lines one loop line, the 6th group of detailed model of air speed data, unit Equivalent Model
With the reactive power simulation result comparison diagram using equivalence method of the present invention to obtain；
When Figure 24 is to point out short circuit in supply lines one loop line, the 6th group of detailed model of air speed data, unit Equivalent Model
With the voltage simulation result comparison diagram using equivalence method of the present invention to obtain；
When Figure 25 is to point out short circuit in supply lines one loop line, the 6th group of detailed model of air speed data, unit Equivalent Model
With the current simulations Comparative result figure using equivalence method of the present invention to obtain；
When Figure 26 is to point out short circuit in supply lines one loop line, the 30th group of detailed model of air speed data, unit equivalence mould
The active power simulation result comparison diagram that type and employing equivalence method of the present invention obtain；
When Figure 27 is to point out short circuit in supply lines one loop line, the 30th group of detailed model of air speed data, unit equivalence mould
The reactive power simulation result comparison diagram that type and employing equivalence method of the present invention obtain；
When Figure 28 is to point out short circuit in supply lines one loop line, the 30th group of detailed model of air speed data, unit equivalence mould
The voltage simulation result comparison diagram that type and employing equivalence method of the present invention obtain；
When Figure 29 is to point out short circuit in supply lines one loop line, the 30th group of detailed model of air speed data, unit equivalence mould
The current simulations Comparative result figure that type and employing equivalence method of the present invention obtain；
Figure 30 is the stream of Double-feed wind power field based on the meritorious response Dynamic Equivalence described in detailed description of the invention one
Cheng Tu.
Detailed description of the invention
Detailed description of the invention one: illustrate present embodiment with reference to Figure 30, described in present embodiment based on meritorious sound
The Double-feed wind power field Dynamic Equivalence answered, it comprises the following steps:
Step one, building the electrical-magnetic model of Double-feed wind power unit, the small-power at this Wind turbines exports respectively
With under two kinds of operating modes of high-power output, carry out three phase short circuit fault emulation, the supplemental characteristic this emulation obtained and on-the-spot test
The supplemental characteristic contrast that experiment obtains, it is thus achieved that simulation result and the weighted average absolute error of field experimentation test result, this is put down
All absolute error conformance with standard required precisions, then perform step 2；
In step 2, acquisition full blast speed operation area, Double-feed wind power unit transient Model occurs three-phase short at also site
Transient state active power response curve during the fault of road, according to the characteristic of bunching on this curve, it is thus achieved that index of hiving off；
Step 3, set meritorious power error threshold value according to actual equivalent required precision, carry out typical wind speed distribution etc.
Value example calculates, and observes the adequacy of step 2 cut-point and necessity, by the meritorious relative systematic error of each example respectively with
The relative error arranging active power allows threshold ratio relatively, verifies each cut-point, if the meritorious equity mutually of certain example
Value error allows threshold value more than the relative error arranging active power, then this example adds the unit in corresponding wind speed section equivalent
Wind speed is as supplementary cut-point；
Step 4, the cut-point obtained according to step 2 and step 3 off-line, take the union of cut-point as index of hiving off
Divide Wind turbines, be one by the Wind turbines equivalence in each group, calculate the equivalent parameters of each equivalent unit and corresponding collection
Electric network equivalent parameters, it is thus achieved that wind energy turbine set Equivalent Model.
In present embodiment, Wind turbines refers in full blast speed operation area cut wind speed to cutting wind speed, between wind speed
It is divided into 0.1m/s.
Embodiment:
The following is one embodiment of the present of invention, as a example by the wind energy turbine set of the 3 × 6 of certain Double-feed wind power unit composition,
Concrete steps and result are as follows:
One, the electromagnetic transient simulation model of Double-feed wind power unit is built, according to " Wind turbines accesses power system technology
Regulation ", under the output of Wind turbines small-power and two kinds of operating modes of high-power output, carry out three phase short circuit fault emulation, and on-the-spot
Test experiments contrasts, if precision meets the requirement of " Wind turbines electrical simulation model modeling directive/guide ", then carries out step 2.
Two, Wind turbines is obtained full blast speed operation area (incision wind speed is spaced apart 0.1m/s to cut-out wind speed, wind speed)
In, and at site occur three phase short circuit fault time transient state active power response curve, according to its characteristic of bunching, it is thus achieved that hive off finger
Mark.
Three, set meritorious power error threshold value according to required precision, attempt typical wind speed distribution example successively, observe segmentation
The adequacy of point and necessity.If precision is unsatisfactory for unit equivalence wind speed that requirement can add in corresponding wind speed section as supplementing point
Cutpoint.
Four, the index of hiving off obtained according to second step and the 3rd step calculated off line divides Wind turbines, by the wind in each group
Group of motors equivalence is one, calculates the equivalent parameters of each equivalent unit and corresponding current collection network equivalent parameter.
Five, utilize the wind energy turbine set Equivalent Model drawn, use actual wind speed data, equivalence method is carried out accuracy and tests
Card.
The Double-feed wind power unit electromagnetic transient simulation model that step one is built is as it is shown in figure 1, unit is by machine end transformation
Device, transmission line of electricity, booster transformer are connected with equivalence electrical network, and simulation parameter is as shown in table 1.During 70s at blower fan site B690
There is three-phase ground short trouble, fault clearance after 0.625s.Measure parameter and be mainly the voltage in Wind turbines exit, electricity
Stream, active power and reactive power.Under the output of unit small-power and two kinds of operating modes of high-power output, during three phase short circuit fault, wind
Active power, reactive power, voltage, the response curve of electric current and the field test data that group of motors exit records to score
Not as shown in Fig. 2 to Fig. 9.
Table 1 simulation parameter
As a example by active power, the definition of A1, B1, B2, C1, C2 time period is as shown in Figure 10.Each time period weighted average
Shown in the computing formula of absolute error such as formula (1), T1For testing the inverse of mains frequency.Wherein, XM(i)、Xs(I) it is i-th point
Measured value and simulation value.Kend, KstartStart sequence number for each time period data.
The error of fitting that can be obtained phantom and measured data by Fig. 2 to Fig. 9 is as shown in table 2.
Table 2 weighted average absolute error (%)
Result shows, the electrical-magnetic model built is consistent with the dynamic characteristic of actual LVRT Capability of Wind Turbine Generator,
Precision meets " the Wind turbines electrical simulation model modeling directive/guide " of China.
In step 2, Wind turbines is the important component part of wind energy turbine set, study its and the voltage of site, electric current, have
Merit, idle response characteristic are by the basis of wind energy turbine set dynamic equivalent.Voltage response depends primarily on the Voltage Drop of setting
The computational accuracy of situation, current-responsive and impedance is closely related, and reactive power depends on the idle benefit of machine unit characteristic and wind energy turbine set
Repay equipment, the most suitably select active power to analyze in its transient response characteristic.
For there is transient state wattful power during three phase short circuit fault at the grid-connected point of Wind turbines in acquisition full blast speed operation area
Rate response curve, blower fan, based on the unit phantom built, is connected by the present invention with Infinite bus system, wind speed from 4.5m/s,
Increasing to 25m/s every 0.1m/s, fault is arranged and main simulation parameter is with above.The meritorious response song measured at B690
Line is as shown in figure 11.Simulation result under each the corresponding a kind of wind speed of curve in figure.
As shown in Figure 11, the wattful power of blower fan take the lead in experience certain time downward cross punching and upwards cross punching after, according to tiltedly
Rate returns to initial value.In addition can also find, when blower fan is operated in different wind speed regions, the transient response of its active power is bent
Line has obvious difference, works in promoter region (reference 1 in Figure 11), maximal power tracing district (reference 2 in Figure 11)
With the Wind turbines of permanent rotating speed invariable power district (reference 3 in Figure 11), there is characteristic of significantly bunching, 7m/s can be identified
(starting point in MPPT district) and two cut-points of 12m/s (terminal in MPPT district).
Be not added with low wear control module time phenomenon different, it shows the difference of meritorious amplitude and degree of fluctuation thereof relatively
Little, active power curves during different wind speed has the strongest roll-off characteristic, and difference is mainly manifested in its transient response speed side
Face.This is because low control module of wearing is provided with the characteristic that active power slope is recovered, in the case of slope sets unanimously, low
The Wind turbines of wind speed can quickly reach steady-state value, and the Wind turbines of high wind speed needs the longer time.
In step 3, the selection principle of index of hiving off for exploration further and checking wind energy turbine set, build wind energy turbine set emulation mould
Type.According to the roll-off characteristic that active power curves is stronger, typical wind speed distribution scene and the meritorious relative systematic error of setting
As shown in table 3, the computing formula of error is as follows:
Wherein, JstartAnd JendAfter being respectively Voltage Drop previous 1 cycle and fault clearance, active power starts stable
The data segment sequence number that 1 cycle after output is corresponding.
Systematic error (%) under the different wind speed scene of table 3
Assume that the relative error of active power allows threshold value to be 5%.As shown in Table 3, wind speed scene 1 and scene 2 times, unit
The relative systematic error of equivalent active power is respectively less than the threshold value set, it was demonstrated that work in promoter region, permanent rotating speed invariable power district
Wind turbines can equivalence be each one.
As shown in Table 3, scene 3 times, the precision of unit equivalence is unsatisfactory for threshold requirement.If scene 3 is with 7m/s as cut-point,
It is that the relatively unit of precision during two machines is doubled left and right by wind energy turbine set equivalence, it was demonstrated that make cut-point with MPPT district top wind speed
Necessity.
Equally, scene 4 times, the precision of unit equivalence is unsatisfactory for threshold requirement, and only according to unit fault power response song
Line None-identified goes out the cut-point in MPPT region.Owing to the linearity of MPPT district response curve is preferable, substantial amounts of l-G simulation test table
Bright, it is that two machine equivalence effect relatively units when two groups improve significantly with its unit equivalent wind speed (9.6m/s) equivalence, precision
Be doubled left and right.
As shown in Table 3, scene 5 times, the error of unit equivalence, seriously beyond threshold range, is divided into two groups at 12m/s
Two machine equivalence effect relatively unit equivalences improve significantly, and its error is about 4.78%, is closer to threshold requirement.For leaving one
Fixed nargin, research finds, at 7m/s, 12m/s, equivalence is three machine equivalence effects when three groups and the response base of detailed model
Originally it is completely superposed error and is about 2.75%.Demonstrate the necessity with MPPT district end wind speed as cut-point.
The simulation result of above example is as shown in table 4.Sum up and find, if the relative error permission threshold value of active power is
5%, with the starting point in maximal power tracing district, terminal and the unit equivalence wind speed of the Wind turbines working in maximal power tracing district
Required precision can be met for cut-point.
Cut-point situation under the different wind speed range of table 4
If it should be noted that the Wind turbines of some producers at strong wind with little wind (to exert oneself as wind corresponding for 0.5p.u.
Speed is for boundary) time use different power slope recovery policies, then there is natural cut-point, be strong wind and the demarcation line of little wind
(generally between 9m/s-10m/s), the equivalent wind speed in maximal power tracing district typically replaces with this segmentation wind speed can expire
Foot required precision.
In step 4, in full blast speed operation area, with at Wind turbines site occur three phase short circuit fault time temporary
The characteristic of bunching of state active power response curve, it is thus achieved that cut-point；Set meritorious power error threshold value according to required precision, depend on
The typical examples of secondary trial wind speed profile, observes adequacy and the necessity of cut-point, it is thus achieved that cut-point.Take the union of cut-point
(7m/s, 9.6m/s, 12m/s) divides Wind turbines as index of hiving off, and is one by the Wind turbines equivalence in each group, calculates
The equivalent parameters of each equivalent unit and corresponding current collection network equivalent parameter.Equivalent wind speed is that all Wind turbines are exerted oneself in this group
Equivalent wind speed corresponding to the mean power of sum.The parameters of equivalent Wind turbines calculates according to capacity weighting method, equivalent
Cable data calculates according to the principle that current collection network power losses before and after equivalence is equal, and detailed process is the present invention will not be described in great detail.
Step 5, carries out Accuracy Verification to equivalence method, it is ensured that the accuracy of wind energy turbine set Equivalent Model is in controlled range
In.Specifically, wind speed can be detected according to least one set, respectively to wind energy turbine set Equivalent Model, conventional individual Equivalent Model and in detail
Thin model carries out fault simulation test respectively, compares wind energy turbine set Equivalent Model, conventional individual Equivalent Model and detailed model same
Active power under sample wind speed, same fault, reactive power, voltage and current response curve, thus judge wind energy turbine set equivalence mould
The simulated effect of type.
This example, in certain 3 × 6 wind energy turbine set May 19 to 2017 groups of air speed datas on May 26, has randomly selected such as figure
Shown in 13 30 group air speed data carries out equivalent experiment.Position of failure point is wind energy turbine set exit, wind energy turbine set and Infinite bus system
Being connected, main simulation parameter is with above.For showing the equivalent effect of the inventive method, the 6th group, the 30th group of wind speed are randomly selected
The effect equivalent with unit is contrasted, and result is respectively as shown in Figure 14 to Figure 17 and Figure 18 to Figure 21.
From Figure 14 to Figure 17 and Figure 18 to Figure 21, use the strategy that hives off of the present invention can significantly improve conventional individual
Polymerization equivalent precision, the active power in wind energy turbine set exit, reactive power, voltage and current tracking effect preferable.Cause
This, the equivalent strategy of the present invention has good adaptability to air speed data.
For verifying the adaptability of this strategy on fault point position of hiving off, trouble point has been moved to double circuit power supply line by this example
The midpoint of one time, utilizes grouping method above, above 30 groups of air speed datas has been re-started simulation calculation.Still choose
The simulation result of 6 groups and the 30th group is shown, respectively as shown in Figure 22 to Figure 25 and Figure 26 to Figure 29.
By Figure 22 to Figure 25 and Figure 26 to Figure 29 it can be seen that the strategy that hives off of the present invention can significantly improve conventional individual
Equivalent equivalent precision, the active power in wind energy turbine set exit, reactive power, voltage and current tracking effect preferable.Cause
This, it also has good adaptability to position of failure point.
Detailed description of the invention two: present embodiment is based on meritorious response double-fed to described in detailed description of the invention one
Wind energy turbine set Dynamic Equivalence is described further, and in present embodiment, in step one, the supplemental characteristic that emulation obtains is wind-powered electricity generation
Voltage, electric current, active power and reactive power at machine set outlet, trouble duration is 0.625s；
In step one, the supplemental characteristic contrast that the supplemental characteristic this emulation obtained and on-the-spot test experiment obtain, it is thus achieved that
Simulation result with the detailed process of the weighted average absolute error of field experimentation test result is:
The active power data of active power response curve with on-the-spot test are contrasted, by reactive power response curve
Contrasting with the reactive power data of on-the-spot test, it is right the response curve of voltage and the voltage data of on-the-spot test to be carried out
Ratio, contrasts the current data of the response curve of electric current with on-the-spot test,
According to formula:
The weighted average of the active power of active power and on-the-spot test in acquisition time period A1, B1, B2, C1 and C2 respectively
Absolute error Fii, reactive power is absolute with the weighted average of the reactive power of on-the-spot test in time period A1, B1, B2, C1 and C2
Error Fii, weighted average absolute error F of voltage and the voltage of on-the-spot test in time period A1, B1, B2, C1 and C2iiAnd the time
The weighted average absolute error of the response curve of the electric current of response curve and the on-the-spot test of electric current in section A1, B1, B2, C1 and C2
Fii,
In formula, Kend, KstartFor the start sequence number of each time period data, XM(i) and Xs(I) actual measurement of i-th point it is respectively
Value and simulation value, F is the weighted average absolute error that voltage, electric current, active power or reactive power are total,
Weighted average absolute error F value is compared with the accuracy value of standard, determines the transient state building Double-feed wind power unit
Model is consistent with the dynamic characteristic of actually measured Wind turbines.
Detailed description of the invention three: present embodiment is based on meritorious response double-fed to described in detailed description of the invention one
Wind energy turbine set Dynamic Equivalence is described further, and in present embodiment, the full blast speed operation area in step 2 is incision wind
Speed is to cut-out wind speed, and wind speed is spaced apart 0.1m/s, and trouble duration is 0.1s.
Detailed description of the invention four: present embodiment is based on meritorious response double-fed to described in detailed description of the invention one
Wind energy turbine set Dynamic Equivalence is described further, in present embodiment, in step 2, according to the characteristic of bunching on this curve,
The detailed process obtaining index of hiving off is:
When blower fan is operated in different wind speed regions, the transient response curve of its active power is different, work in promoter region,
The Wind turbines in maximal power tracing district and permanent rotating speed invariable power district has characteristic of significantly bunching, it is possible to identify MPPT district
Two cut-points of the terminal in starting point and MPPT district.
Detailed description of the invention five: present embodiment is based on meritorious response double-fed to described in detailed description of the invention one
Wind energy turbine set Dynamic Equivalence is described further, and in present embodiment, trouble duration is 0.1s, in step 3, according to
The detailed process of the meritorious relative systematic error that the index of hiving off in step 2 obtains each example is:
According to formula:
Obtain the meritorious relative systematic error E of each examplep,
In formula, JstartAnd JendAfter being respectively Voltage Drop previous 1 cycle and fault clearance, active power starts stable
The data segment sequence number that 1 cycle after output is corresponding, XM(i) and XsI () is respectively measured value and the simulation value of i-th point.
Detailed description of the invention six: present embodiment is based on meritorious response double-fed to described in detailed description of the invention one
Wind energy turbine set Dynamic Equivalence is described further, and in present embodiment, calculates the equivalent parameters of each equivalent unit with corresponding
The detailed process of current collection network equivalent parameter is:
Each parameter of equivalent Wind turbines calculates according to capacity weighting method, before and after the parameter of equivalent collection electric network is according to equivalence
The principle that current collection network power losses is equal calculates.
Detailed description of the invention seven: present embodiment is based on meritorious response double-fed to described in detailed description of the invention one
Wind energy turbine set Dynamic Equivalence is described further, in present embodiment, in step 3, and unit equivalence wind speed veqComputing formula
For:
In formula, m is the number of units of Wind turbines, and j is the numbering of unit, PjMeritorious exerting oneself for jth platform unit.Pj=f (vj)
Represent the wind speed-power characteristic of unit.
Detailed description of the invention eight: present embodiment is based on meritorious response double-fed to described in detailed description of the invention one
Wind energy turbine set Dynamic Equivalence is described further, in present embodiment, further comprising the steps of after step 4,
Utilize the wind energy turbine set Equivalent Model that step 4 draws, use the accurate of actual wind speed this Equivalent Model of data verification
Property.
Claims (5)
1. Double-feed wind power field based on meritorious response Dynamic Equivalence, it is characterised in that it comprises the following steps:
Step one, building the electrical-magnetic model of Double-feed wind power unit, the small-power at this Wind turbines exports and big respectively
Under power two kinds of operating modes of output, carrying out three phase short circuit fault emulation, the supplemental characteristic this emulation obtained is tested with on-the-spot test
The supplemental characteristic contrast obtained, it is thus achieved that simulation result and the weighted average absolute error of field experimentation test result, this weighting is put down
All absolute error conformance with standard required precisions, then perform step 2；
In step 2, acquisition full blast speed operation area there is three-phase shortcircuit event in Double-feed wind power unit transient Model at also site
Transient state active power response curve during barrier, according to the characteristic of bunching on this curve, it is thus achieved that index of hiving off；
Step 3, power error threshold value of gaining merit according to actual equivalent required precision setting, carry out the equivalent calculation of typical wind speed distribution
Example calculates, it is thus achieved that the meritorious relative systematic error of each example, is gained merit with setting respectively by the meritorious relative systematic error of each example
Cut-point relatively, is verified by power error threshold ratio, if the meritorious relative systematic error of certain example is more than setting wattful power
Rate error threshold, then add the unit equivalence wind speed in corresponding wind speed section as supplementary cut-point to this example；
Step 4, the cut-point obtained according to step 2 and step 3 off-line, the union taking cut-point divides as index of hiving off
Wind turbines, is one by the Wind turbines equivalence in each group, calculates the equivalent parameters of each equivalent unit and collects electrical network accordingly
Network equivalent parameters, it is thus achieved that wind energy turbine set Equivalent Model.
Double-feed wind power field based on meritorious response the most according to claim 1 Dynamic Equivalence, it is characterised in that
In step one, the supplemental characteristic that emulation obtains is the voltage in Wind turbines exit, electric current, active power and idle merit
Rate；
In step one, the supplemental characteristic contrast that the supplemental characteristic this emulation obtained and on-the-spot test experiment obtain, it is thus achieved that emulation
Result with the detailed process of the weighted average absolute error of field experimentation test result is:
The active power data of active power response curve with on-the-spot test are contrasted, by reactive power response curve with existing
The reactive power data of field test contrast, and are contrasted by the voltage data of the response curve of voltage with on-the-spot test, will
The response curve of electric current contrasts with the current data of on-the-spot test,
According to formula:
In obtaining time period A1, B1, B2, C1 and C2 respectively, active power is absolute with the weighted average of the active power of on-the-spot test
Error Fii, the weighted average absolute error of reactive power and the reactive power of on-the-spot test in time period A1, B1, B2, C1 and C2
Fii, weighted average absolute error F of voltage and the voltage of on-the-spot test in time period A1, B1, B2, C1 and C2iiAnd the time period
The weighted average absolute error of the response curve of the electric current of the response curve of electric current and on-the-spot test in A1, B1, B2, C1 and C2
Fii,
In formula, Kend, KstartFor the start sequence number of each time period data, XM(i) and Xs(I) be respectively i-th point measured value and
Simulation value, F is the weighted average absolute error that voltage, electric current, active power or reactive power are total,
Weighted average absolute error F value is compared with the accuracy value of standard, determines the transient Model building Double-feed wind power unit
Consistent with the dynamic characteristic of actually measured Wind turbines.
Double-feed wind power field based on meritorious response the most according to claim 1 Dynamic Equivalence, it is characterised in that step
Full blast speed operation area in rapid two is spaced apart 0.1m/s for incision wind speed to cut-out wind speed, wind speed.
Double-feed wind power field based on meritorious response the most according to claim 1 Dynamic Equivalence, it is characterised in that step
In rapid two, according to the characteristic of bunching on this curve, it is thus achieved that the detailed process of index of hiving off is:
When blower fan is operated in different wind speed regions, the transient response curve of its active power is different, works in promoter region, maximum
The Wind turbines in power tracking district and permanent rotating speed invariable power district has characteristic of significantly bunching, it is possible to identify the starting point in MPPT district
Two cut-points of terminal with MPPT district.
Double-feed wind power field based on meritorious response the most according to claim 1 Dynamic Equivalence, it is characterised in that step
In rapid three, it is thus achieved that the detailed process of the meritorious relative systematic error of each example is:
According to formula:
Obtain the meritorious relative systematic error E of each examplep,
In formula, JstartAnd JendAfter being respectively Voltage Drop previous 1 cycle and fault clearance, active power starts stable output
After data segment sequence number corresponding to 1 cycle, XM(i) and XsI () is respectively measured value and the simulation value of i-th point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610594669.7A CN106202815B (en) | 2016-07-26 | 2016-07-26 | Double-feed wind power field Dynamic Equivalence based on active response |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610594669.7A CN106202815B (en) | 2016-07-26 | 2016-07-26 | Double-feed wind power field Dynamic Equivalence based on active response |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106202815A true CN106202815A (en) | 2016-12-07 |
| CN106202815B CN106202815B (en) | 2019-04-26 |
Family
ID=57495236
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610594669.7A Active CN106202815B (en) | 2016-07-26 | 2016-07-26 | Double-feed wind power field Dynamic Equivalence based on active response |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106202815B (en) |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107154624A (en) * | 2017-04-25 | 2017-09-12 | 华北电力大学 | A kind of appraisal procedure of wind power plant Equivalent Model accuracy |
| CN107272643A (en) * | 2017-06-19 | 2017-10-20 | 哈尔滨工业大学 | The practical unit equivalence method in new energy station |
| CN109407543A (en) * | 2018-01-29 | 2019-03-01 | 中国电力科学研究院有限公司 | A kind of verification method and device of Wind turbines electrical model voltage responsive characteristic |
| CN109522607A (en) * | 2018-10-22 | 2019-03-26 | 国网江西省电力有限公司电力科学研究院 | A kind of double-fed fan motor field electromechanical transient equivalent modeling method |
| CN109521693A (en) * | 2018-12-03 | 2019-03-26 | 中国电力科学研究院有限公司 | A kind of wind power plant electro-magnetic transient real-time emulation method and device |
| CN110032778A (en) * | 2019-03-26 | 2019-07-19 | 河海大学 | A kind of wind power plant single machine equivalent modeling verification method |
| CN110119570A (en) * | 2019-05-10 | 2019-08-13 | 东北电力大学 | A kind of farm model parameters validation method of measured data driving |
| CN110175933A (en) * | 2019-05-29 | 2019-08-27 | 哈尔滨工业大学 | Direct-driving type wind power plant Dynamic Equivalence suitable for frequency modulation control |
| CN113627035A (en) * | 2021-09-01 | 2021-11-09 | 国网黑龙江省电力有限公司电力科学研究院 | Method for establishing wind power plant single-machine equivalent model |
| CN114188943A (en) * | 2021-12-15 | 2022-03-15 | 河海大学 | Photovoltaic field group dynamic equivalent modeling method and system integrating voltage drop depth and voltage fluctuation similarity |
| CN114357787A (en) * | 2022-01-10 | 2022-04-15 | 华北电力大学 | Equivalent modeling method and system for offshore wind plant |
| CN115882526A (en) * | 2023-02-13 | 2023-03-31 | 华北电力科学研究院有限责任公司 | Wind power plant group-oriented active and reactive power regulation capacity data processing method and device |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104036073A (en) * | 2014-05-23 | 2014-09-10 | 电子科技大学 | Double-fed wind power plant dynamic equivalence modeling method suitable for active power characteristic analysis |
| CN105303319A (en) * | 2015-11-01 | 2016-02-03 | 国网山东省电力公司聊城供电公司 | Wind power plant dynamic clustering equivalence method based on operating data |
| CN105447255A (en) * | 2015-12-03 | 2016-03-30 | 哈尔滨工业大学 | Dynamic equivalence and modeling methods and systems for double-fed wind power plant |
| CN105576654A (en) * | 2016-01-08 | 2016-05-11 | 南方电网科学研究院有限责任公司 | Equivalent method and system of direct-driven wind power plant |
| CN105678033A (en) * | 2016-03-04 | 2016-06-15 | 国家电网公司 | Wind farm equivalent modeling method suitable for electromagnetism transient simulation |
| CN105703396A (en) * | 2016-03-21 | 2016-06-22 | 广西电网有限责任公司电力科学研究院 | Multi-turbine characterization equivalent modeling method of double-fed wind power plant |
| CN105703364A (en) * | 2016-04-18 | 2016-06-22 | 哈尔滨工业大学 | Photovoltaic power station equivalent modeling method |
-
2016
- 2016-07-26 CN CN201610594669.7A patent/CN106202815B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104036073A (en) * | 2014-05-23 | 2014-09-10 | 电子科技大学 | Double-fed wind power plant dynamic equivalence modeling method suitable for active power characteristic analysis |
| CN105303319A (en) * | 2015-11-01 | 2016-02-03 | 国网山东省电力公司聊城供电公司 | Wind power plant dynamic clustering equivalence method based on operating data |
| CN105447255A (en) * | 2015-12-03 | 2016-03-30 | 哈尔滨工业大学 | Dynamic equivalence and modeling methods and systems for double-fed wind power plant |
| CN105576654A (en) * | 2016-01-08 | 2016-05-11 | 南方电网科学研究院有限责任公司 | Equivalent method and system of direct-driven wind power plant |
| CN105678033A (en) * | 2016-03-04 | 2016-06-15 | 国家电网公司 | Wind farm equivalent modeling method suitable for electromagnetism transient simulation |
| CN105703396A (en) * | 2016-03-21 | 2016-06-22 | 广西电网有限责任公司电力科学研究院 | Multi-turbine characterization equivalent modeling method of double-fed wind power plant |
| CN105703364A (en) * | 2016-04-18 | 2016-06-22 | 哈尔滨工业大学 | Photovoltaic power station equivalent modeling method |
Non-Patent Citations (4)
| Title |
|---|
| PUPU CHAO等: "Investigation on Fundamental Issues of Dynamic Equivalent Modeling of Wind Farms", 《PREPRINTS OF THE 5TH INTERNATIONAL CONFERENCE ON ELECTRIC UTILITY DEREGULATION AND RESTRUCTURING AND POWER TECHNOLOGIES》 * |
| 张洋溢: "风电场动态等值建模方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
| 朱乾龙等: "基于聚类–判别分析的风电场概率等值建模研究", 《中国电机工程学报》 * |
| 蒙晓航: "永磁直驱同步风电场多机动态等值模型", 《电力系统保护与控制》 * |
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107154624A (en) * | 2017-04-25 | 2017-09-12 | 华北电力大学 | A kind of appraisal procedure of wind power plant Equivalent Model accuracy |
| CN107272643A (en) * | 2017-06-19 | 2017-10-20 | 哈尔滨工业大学 | The practical unit equivalence method in new energy station |
| CN107272643B (en) * | 2017-06-19 | 2019-03-29 | 哈尔滨工业大学 | The single machine equivalence method of new energy station |
| CN109407543A (en) * | 2018-01-29 | 2019-03-01 | 中国电力科学研究院有限公司 | A kind of verification method and device of Wind turbines electrical model voltage responsive characteristic |
| CN109407543B (en) * | 2018-01-29 | 2023-11-14 | 中国电力科学研究院有限公司 | Verification method and device for voltage response characteristics of electrical model of wind turbine generator |
| CN109522607A (en) * | 2018-10-22 | 2019-03-26 | 国网江西省电力有限公司电力科学研究院 | A kind of double-fed fan motor field electromechanical transient equivalent modeling method |
| CN109521693B (en) * | 2018-12-03 | 2023-05-16 | 中国电力科学研究院有限公司 | Wind power plant electromagnetic transient real-time simulation method and device |
| CN109521693A (en) * | 2018-12-03 | 2019-03-26 | 中国电力科学研究院有限公司 | A kind of wind power plant electro-magnetic transient real-time emulation method and device |
| CN110032778A (en) * | 2019-03-26 | 2019-07-19 | 河海大学 | A kind of wind power plant single machine equivalent modeling verification method |
| CN110119570A (en) * | 2019-05-10 | 2019-08-13 | 东北电力大学 | A kind of farm model parameters validation method of measured data driving |
| CN110119570B (en) * | 2019-05-10 | 2023-08-29 | 东北电力大学 | Actually measured data driven wind farm model parameter checking method |
| CN110175933A (en) * | 2019-05-29 | 2019-08-27 | 哈尔滨工业大学 | Direct-driving type wind power plant Dynamic Equivalence suitable for frequency modulation control |
| CN113627035A (en) * | 2021-09-01 | 2021-11-09 | 国网黑龙江省电力有限公司电力科学研究院 | Method for establishing wind power plant single-machine equivalent model |
| CN113627035B (en) * | 2021-09-01 | 2024-01-19 | 国网黑龙江省电力有限公司电力科学研究院 | Method for establishing single-machine equivalent model of wind power plant |
| CN114188943A (en) * | 2021-12-15 | 2022-03-15 | 河海大学 | Photovoltaic field group dynamic equivalent modeling method and system integrating voltage drop depth and voltage fluctuation similarity |
| CN114188943B (en) * | 2021-12-15 | 2023-09-08 | 河海大学 | Dynamic equivalent modeling method and system for photovoltaic field group |
| CN114357787A (en) * | 2022-01-10 | 2022-04-15 | 华北电力大学 | Equivalent modeling method and system for offshore wind plant |
| CN114357787B (en) * | 2022-01-10 | 2024-05-28 | 华北电力大学 | Offshore wind farm equivalent modeling method and system |
| CN115882526A (en) * | 2023-02-13 | 2023-03-31 | 华北电力科学研究院有限责任公司 | Wind power plant group-oriented active and reactive power regulation capacity data processing method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106202815B (en) | 2019-04-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106202815A (en) | Double-feed wind power field based on meritorious response Dynamic Equivalence | |
| Ellis et al. | Description and technical specifications for generic WTG models—A status report | |
| CN101937483B (en) | Dynamic equivalence modeling method for double-fed set wind power station based on action condition of pitch angle | |
| CN106383947B (en) | The fast acquiring method of wind power plant current collection network dynamic equivalent parameters | |
| CN103208798B (en) | Method for calculating probability power flow of power system containing wind farm | |
| CN102760191B (en) | Double-fed machine set wind power station equivalent modeling system and method based on rotation speed grouping | |
| CN103761349B (en) | A kind of wind energy turbine set equivalent modeling method based on Wind turbines probability same tone | |
| CN103700036B (en) | A kind of transient stability projecting integral method being suitable to power system Multiple Time Scales | |
| CN104392056B (en) | A kind of directly driven wind-powered unit modeling method and device | |
| CN103219725A (en) | Wind power plant equivalent modeling method based on real-time operation data | |
| CN105678033B (en) | A kind of wind-powered electricity generation station equivalent modeling method suitable for electromagnetic transient simulation | |
| CN112818491A (en) | Wind power plant aggregation equivalent modeling method based on principal component analysis and clustering algorithm | |
| CN102819221B (en) | LVRT Capability of Wind Turbine Generator characteristic associative simulation model and joint simulation method thereof | |
| CN106526347A (en) | Digital-analog hybrid simulation-based photovoltaic inverter low voltage ride through evaluation method | |
| CN105576654B (en) | Directly driven wind-powered field equivalence method and system | |
| CN105243230B (en) | A kind of equivalent modeling method of the directly driven wind-powered field of meter and discharging circuit | |
| CN104820741A (en) | Wind power plant dynamic equivalence method combining wind field disperstiveness and unit difference | |
| CN109583069A (en) | Wind speed modeling method based on temporal correlation | |
| CN104504285A (en) | Doubly-fed wind power farm equivalent modeling method for crowbar protection | |
| CN106410862A (en) | Wind power plant single machine equivalent method based on active recovery slope correction | |
| CN105701265B (en) | A kind of double-fed fan motor unit modeling method and device | |
| CN102510072A (en) | Power grid system transient destabilization differentiation method | |
| CN104268316B (en) | Off-line calculation and online matching based doubly-fed wind power station probability equivalent modeling method | |
| CN103853888B (en) | Equivalent modeling method for offshore wind power system suitable for reactive compensation | |
| CN107167733B (en) | A kind of acquisition methods of the basal evaluation data of synchronous generator excited system performance |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |