CN106202815B - Double-feed wind power field Dynamic Equivalence based on active response - Google Patents
Double-feed wind power field Dynamic Equivalence based on active response Download PDFInfo
- Publication number
- CN106202815B CN106202815B CN201610594669.7A CN201610594669A CN106202815B CN 106202815 B CN106202815 B CN 106202815B CN 201610594669 A CN201610594669 A CN 201610594669A CN 106202815 B CN106202815 B CN 106202815B
- Authority
- CN
- China
- Prior art keywords
- power
- wind
- active
- equivalent
- 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.)
- Active
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)
- Wind Motors (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
Double-feed wind power field Dynamic Equivalence based on active response, is related to electric system simulation modeling technique field.The present invention is to solve the problem of that existing wind power plant equivalence method cannot be considered in terms of calculation amount and equivalent precision while can not complete offline blower and divide group.It builds Wind turbines electrical-magnetic model and carries out three-phase shortcircuit emulation and on-the-spot test Experimental comparison, transient state active power response curve of the Wind turbines in full blast speed operation area is obtained if precision complies with standard, divide group's index according to its characteristic acquisition of bunching；Active power error threshold value is set according to equivalent required precision, the equivalent example for carrying out typical wind speed distribution calculates, if precision is unsatisfactory for requiring, single machine equivalence wind speed is as supplement cut-point in the corresponding wind speed section of addition；Divide group's index according to offline acquisition, be one by unit equivalence in each group of practical wind power plant, the equivalent parameters for calculating each equivalent unit sum aggregate electric network obtains wind power plant Equivalent Model.For establishing wind power plant Equivalent Model.
Description
Technical field
The present invention relates to a kind of Dynamic Equivalences of Double-feed wind power field, belong to electric system simulation modeling technique neck
Domain.
Background technique
It is increasingly serious due to energy and environmental problem, wind energy with the advantage in terms of its technology and cost in the power system
It is widely applied.However, quickling increase with wind energy permeability, the operational reliability and stability of electric system by
Certain influence and impact, therefore access to wind-powered electricity generation electric system on a large scale and carry out the accuracy requirement of analytical calculation increasingly
It is high.One wind power plant often has tens even up to a hundred Wind turbines, each Wind turbines again by wind energy conversion system, transmission shaft,
The modules composition such as generator, current transformer and its controller, protective device and its controller, if all independent to every Wind turbines
Modeling can then greatly increase the complexity of simulation model and calculate the time.Therefore, the research of wind power plant Equivalent Model is carried out very
It is necessary.
In wind power plant equivalent modeling, wanted since the fitting precision of large-scale wind power field single machine equivalence is generally difficult to meet
It asks, how Wind turbines is carried out that group is rationally effectively divided to be the matter of utmost importance that wind power plant equivalent modeling research needs to solve.Mesh
Before until, the grouping method of wind power plant mainly has following two categories:
(1) grouping method of machine group cluster is carried out according to the similitude of wind speed.This method is made according to wake effect and wind direction
At wind speed difference, wind power plant is divided into multiple regions, respectively corresponding one equivalent unit.Or direct basis wind speed is close
Property, divide group according to the working region of Wind turbines.This kind of grouping method only considered Index For Steady-state, when wind power plant inner blower quantity
When numerous, operating condition differs greatly, systematic error can be larger.
(2) the characteristic quantity of operating states of the units can be characterized as the grouping method of index.It is linearized with Wind turbines
The variable that the influence factor of state dominant characteristics root or revolving speed, pitch angle etc. can characterize unit working condition, which is used as, divides group to refer to
Mark, although this kind of grouping method can achieve higher equivalent precision, but generally require with more complex algorithm, calculation amount
Greatly.Operating condition difference can still result in the increase of equivalent wind turbine number of units when big, and when inputting wind speed variation, can not be divided in real time
Group's index, engineering, which uses, has biggish limitation.
Summary of the invention
The present invention be in order to solve existing wind power plant equivalence method and cannot be considered in terms of calculation amount and equivalent precision, while can not
It is offline to complete the problem of blower divides group.Double-feed wind power field Dynamic Equivalence based on active response is now provided.
Double-feed wind power field Dynamic Equivalence based on active response, it the following steps are included:
Step 1: the transient Model of Double-feed wind power unit is built, respectively in the small-power output of the Wind turbines and greatly
Under two kinds of operating conditions of power output, three phase short circuit fault emulation is carried out, the supplemental characteristic which is obtained and on-the-spot test are tested
The supplemental characteristic of acquisition compares, and obtains the weighted average absolute error of simulation result and field experiment test result, and the weighting is flat
Equal absolute error complies with standard required precision, thens follow the steps two；
Step 2: obtaining in full blast speed operation area, that three-phase occurs at grid entry point is short for Double-feed wind power unit transient Model
Transient state active power response curve when the failure of road, according to the characteristic of bunching on the curve, acquisition divides group's index；
Step 3: according to the active power error threshold value of practical equivalent required precision setting, progress typical wind speed distribution etc.
Be worth example to calculate, obtain the active opposite systematic error of each example, by the active opposite systematic error of each example respectively with setting
Active power error threshold compares, and verifies to cut-point, is set with if the active opposite systematic error of some example is greater than
Function power error threshold value then adds the single machine equivalence wind speed in corresponding wind speed section as supplement cut-point to the example；
Step 4: taking the union of cut-point to be used as according to the cut-point that step 2 and step 3 obtain offline and dividing group's index
Wind turbines are divided, are one by the Wind turbines equivalence in each group, calculate the equivalent parameters and corresponding collection of each equivalent unit
Electric network equivalent parameters obtains wind power plant Equivalent Model.
The invention has the benefit that the present invention is 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 conditions of high-power output, three phase short circuit fault is set, Wind turbines exit is measured
Active power, reactive power, voltage, the response curve of electric current and field test data compare, and are weighted by each period
Mean absolute error formula obtain simulation parameter and field test data error of fitting, the error required precision range it
It is interior, illustrate that simulation model is consistent with the dynamic characteristic of practical LVRT Capability of Wind Turbine Generator；Then active power is selected to analyze
Its transient response characteristic, blower work at different wind speed regions, and the transient response curve of active power has apparent poor
It is different, according to the characteristic of bunching on the active power curve, it can be seen that work in promoter region, maximal power tracing area and perseverance and turn
The Wind turbines in fast invariable power area, which have, significantly bunches characteristic, can identify 7m/s (starting point in maximal power tracing area) and
Two cut-points of 12m/s (terminal in maximal power tracing area).According to the stronger roll-off characteristic of active power curves, setting is typical
Wind speed profile scene and active opposite systematic error threshold value have the active opposite systematic error of each scene with setting respectively
Function power relative error allows threshold value comparison, verifies to the cut-point of each scene, if some scene is active relatively equivalent
The active power relative error that error is greater than setting allows threshold value, then the single machine equivalence wind in corresponding wind speed section is added to the scene
Speed is as supplement cut-point；Finally, obtaining the cut-point of each scene, takes the union of cut-point to be used as and group's index is divided to divide wind turbine
Wind turbines equivalence in each group is one by group, calculates the equivalent parameters and corresponding current collection network equivalent of each equivalent unit
Parameter obtains wind power plant Equivalent Model, the wind power plant Equivalent Model obtained using step 4, using actual wind speed data verification etc.
It is worth the accuracy of model.
It cannot be considered in terms of calculation amount and equivalent precision the method overcome current equivalence method and wind-powered electricity generation can not be completed offline
The problem of machine component group, can be obviously improved the accuracy of conventional individual equivalence, and principle is simple, and physical meaning is clear, not need
Complicated calculating has good adaptability to air speed data and position of failure point.
Detailed description of the invention
The simulation model for the Double-feed wind power unit that Fig. 1 is built for the present invention；
Fig. 2 is Double-feed wind power unit under the operating condition that small-power exports, having when generation three-phase shortcircuit at fan outlet
The comparison diagram of function Simulation of SAR power image and measured result；
Fig. 3 is Double-feed wind power unit under the operating condition that small-power exports, and nothing when three-phase shortcircuit occurs at fan outlet
The comparison diagram of function Simulation of SAR power image and measured result；
Fig. 4 is Double-feed wind power unit under the operating condition that small-power exports, and electricity when three-phase shortcircuit occurs at fan outlet
The comparison diagram of pressure emulation and measured result；
Fig. 5 is Double-feed wind power unit under the operating condition that small-power exports, and electricity when three-phase shortcircuit occurs at fan outlet
The comparison diagram of stream emulation and measured result；
Fig. 6 is Double-feed wind power unit under the operating condition of high-power output, having when generation three-phase shortcircuit at fan outlet
The comparison diagram of function Simulation of SAR power image and measured result；
Fig. 7 is Double-feed wind power unit under the operating condition of high-power output, and nothing when three-phase shortcircuit occurs at fan outlet
The comparison diagram of function Simulation of SAR power image and measured result；
Fig. 8 is Double-feed wind power unit under the operating condition of high-power output, and electricity when three-phase shortcircuit occurs at fan outlet
The comparison diagram of pressure emulation and measured result；
Fig. 9 is Double-feed wind power unit under the operating condition of high-power output, and electricity when three-phase shortcircuit occurs at fan outlet
The comparison diagram of stream emulation and measured result；
Figure 10 is the definition schematic diagram of error calculation period, at the time of map parameter t indicates that failure occurs, parameter t0Generation
The time of table crowbar circuit investment, parameter t1Representing fault duration, parameter t3-t1When needed for representing fault recovery
Between, T1For the parameter reciprocal for testing mains frequency, t+t3At the time of representative enters stable state；
Figure 11 is the active power response curve in the full wind speed range of Wind turbines, and appended drawing reference 1 indicates promoter region, attached drawing
Label 2 indicates maximal power tracing area, and appended drawing reference 3 indicates permanent revolving speed invariable power area；
Figure 12 is the partial enlarged view at the 70s to 70.2s in Figure 11；
Figure 13 is the air speed data of l-G simulation test；
When Figure 14 is the short circuit of wind power plant exit, detailed model, single machine Equivalent Model and the use of the 6th group of air speed data are originally
The active power simulation result comparison diagram that inventive method obtains；
When Figure 15 is the short circuit of wind power plant exit, detailed model, single machine Equivalent Model and the use of the 6th group of air speed data are originally
The reactive power simulation result comparison diagram that inventive method obtains；
When Figure 16 is the short circuit of wind power plant exit, detailed model, single machine Equivalent Model and the use of the 6th group of air speed data are originally
The voltage simulation result comparison diagram that inventive method obtains；
When Figure 17 is the short circuit of wind power plant exit, detailed model, single machine Equivalent Model and the use of the 6th group of air speed data are originally
The current simulations comparative result figure that inventive method obtains；
When Figure 18 is the short circuit of wind power plant exit, detailed model, single machine Equivalent Model and the use of the 30th group of air speed data
The active power simulation result comparison diagram that the method for the present invention obtains；
When Figure 19 is the short circuit of wind power plant exit, detailed model, single machine Equivalent Model and the use of the 30th group of air speed data
The reactive power simulation result comparison diagram that the method for the present invention obtains；
When Figure 20 is the short circuit of wind power plant exit, detailed model, single machine Equivalent Model and the use of the 30th group of air speed data
The voltage simulation result comparison diagram that the method for the present invention obtains；
When Figure 21 is the short circuit of wind power plant exit, detailed model, single machine Equivalent Model and the use of the 30th group of air speed data
The current simulations comparative result figure that the method for the present invention obtains；
Figure 22 is the detailed model of the 6th group of air speed data, single machine Equivalent Model when pointing out short circuit in one loop line of supply lines
With the active power simulation result comparison diagram obtained using equivalence method of the present invention；
Figure 23 is the detailed model of the 6th group of air speed data, single machine Equivalent Model when pointing out short circuit in one loop line of supply lines
With the reactive power simulation result comparison diagram obtained using equivalence method of the present invention；
Figure 24 is the detailed model of the 6th group of air speed data, single machine Equivalent Model when pointing out short circuit in one loop line of supply lines
With the voltage simulation result comparison diagram obtained using equivalence method of the present invention；
Figure 25 is the detailed model of the 6th group of air speed data, single machine Equivalent Model when pointing out short circuit in one loop line of supply lines
With the current simulations comparative result figure obtained using equivalence method of the present invention；
Figure 26 is detailed model, the single machine equivalence mould of the 30th group of air speed data when pointing out short circuit in one loop line of supply lines
Type and the active power simulation result comparison diagram obtained using equivalence method of the present invention；
Figure 27 is detailed model, the single machine equivalence mould of the 30th group of air speed data when pointing out short circuit in one loop line of supply lines
Type and the reactive power simulation result comparison diagram obtained using equivalence method of the present invention；
Figure 28 is detailed model, the single machine equivalence mould of the 30th group of air speed data when pointing out short circuit in one loop line of supply lines
Type and the voltage simulation result comparison diagram obtained using equivalence method of the present invention；
Figure 29 is detailed model, the single machine equivalence mould of the 30th group of air speed data when pointing out short circuit in one loop line of supply lines
Type and the current simulations comparative result figure obtained using equivalence method of the present invention；
Figure 30 is the stream of the Double-feed wind power field Dynamic Equivalence based on active response described in specific embodiment one
Cheng Tu.
Specific embodiment
Specific embodiment 1: present embodiment is illustrated referring to Figure 30, based on active sound described in present embodiment
The Double-feed wind power field Dynamic Equivalence answered, it the following steps are included:
Step 1: building the electrical-magnetic model of Double-feed wind power unit, exported respectively in the small-power of the Wind turbines
Under two kinds of operating conditions of high-power output, three phase short circuit fault emulation, the supplemental characteristic which is obtained and on-the-spot test are carried out
The supplemental characteristic comparison that experiment obtains, obtains the weighted average absolute error of simulation result and field experiment test result, this is flat
Equal absolute error complies with standard required precision, thens follow the steps two；
Step 2: obtaining in full blast speed operation area, that three-phase occurs at grid entry point is short for Double-feed wind power unit transient Model
Transient state active power response curve when the failure of road, according to the characteristic of bunching on the curve, acquisition divides group's index；
Step 3: according to the active power error threshold value of practical equivalent required precision setting, progress typical wind speed distribution etc.
Be worth example to calculate, the adequacy and necessity of two cut-point of observation of steps, by the active opposite systematic error of each example respectively with
The relative error that active power is arranged allows threshold value comparison, verifies to each cut-point, if some example it is active relatively etc.
It is worth the relative error permission threshold value that error is greater than setting active power, then it is equivalent to add the single machine in corresponding wind speed section to the example
Wind speed is as supplement cut-point；
Step 4: taking the union of cut-point to be used as according to the cut-point that step 2 and step 3 obtain offline and dividing group's index
Wind turbines are divided, are one by the Wind turbines equivalence in each group, calculate the equivalent parameters and corresponding collection of each equivalent unit
Electric network equivalent parameters obtains wind power plant Equivalent Model.
In present embodiment, Wind turbines refer to incision wind speed to cutting wind speed, between wind speed in full blast speed operation area
It is divided into 0.1m/s.
Embodiment:
It is one embodiment of the present of invention below, by taking 3 × 6 wind power plant of certain Double-feed wind power unit composition as an example,
Specific steps and result are as follows:
One, the electromagnetic transient simulation model for building Double-feed wind power unit, according to " Wind turbines access power system technology
Regulation ", under the output of Wind turbines small-power and two kinds of operating conditions of high-power output, carry out three phase short circuit fault emulation, and scene
Test experiments comparison, if precision meets the requirement of " Wind turbines electrical simulation model modeling directive/guide ", carries out step 2.
Two, Wind turbines are obtained full blast speed operation area (incision wind speed is divided into 0.1m/s between wind speed to cut-out wind speed)
It is interior, transient state active power response curve when three phase short circuit fault occurs at grid entry point, according to its characteristic of bunching, acquisition divides group to refer to
Mark.
Three, active power error threshold value is set according to required precision, successively attempts typical wind speed and is distributed example, observation segmentation
The adequacy and necessity of point.If the single machine equivalence wind speed that precision is unsatisfactory for requirement and can add in corresponding wind speed section divides as supplement
Cutpoint.
Four, group's index is divided to divide Wind turbines according to what second step and third step off-line calculation obtained, by the wind in each group
Motor group equivalence is one, the equivalent parameters of each equivalent unit of calculating and corresponding current collection network equivalent parameter.
Five, accuracy is carried out to equivalence method and is tested using actual wind speed data using the wind power plant Equivalent Model obtained
Card.
The Double-feed wind power unit electromagnetic transient simulation model that step 1 is built is as shown in Figure 1, unit passes through generator terminal transformation
Device, transmission line of electricity, step-up transformer are connected with equivalent power grid, and simulation parameter is as shown in table 1.When 70s at blower grid entry point B690
Three-phase ground short trouble, fault clearance after 0.625s occurs.Measurement parameter is 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 conditions of high-power output, when three phase short circuit fault, wind
Active power, reactive power, voltage, the response curve of electric current and the field test data that motor group exit measures to score
Not not as shown in Fig. 2 to Fig. 9.
1 simulation parameter of table
By taking active power as an example, the definition of A1, B1, B2, C1, C2 period are as shown in Figure 10.Each period weighted average
Shown in the calculation formula of absolute error such as formula (1), T1For the inverse for testing mains frequency.Wherein, XM(i)、Xs(It i) is i-th point
Measured value and simulation value.Kend, KstartFor the start sequence number of each time hop counts evidence.
The error of fitting that simulation model and measured data can be obtained by Fig. 2 to Fig. 9 is as shown in table 2.
Table 2 is weighted and averaged absolute error (%)
The result shows that the electrical-magnetic model built is consistent with the dynamic characteristic of practical LVRT Capability of Wind Turbine Generator,
Precision meets " the Wind turbines electrical simulation model modeling directive/guide " in China.
In step 2, Wind turbines are the important components of wind power plant, study its in the voltage of grid entry point, electric current, have
Function, idle response characteristic are to carry out the basis of wind power plant dynamic equivalent.Voltage responsive 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 power plant
Equipment is repaid, therefore is suitable for selecting active power to analyze in its transient response characteristic.
Transient state wattful power when three phase short circuit fault occurs in full blast speed operation area at Wind turbines grid entry point to obtain
Rate response curve, the present invention is based on the single machine simulation models built, and blower is connected with Infinite bus system, wind speed from 4.5m/s,
25m/s is increased to every 0.1m/s, fault setting and main simulation parameter are the same as above.The active response measured at B690 is bent
Line is as shown in figure 11.Each curve corresponds to the simulation result under a kind of wind speed in figure.
As shown in Figure 11, the wattful power of blower takes the lead in after undergoing the downward overshoot of certain time and overshooting upwards, according to oblique
Rate is restored to initial value.Furthermore it can also find, blower works at different wind speed regions, and the transient response of active power is bent
Line has apparent difference, works in promoter region (appended drawing reference 1 in Figure 11), maximal power tracing area (appended drawing reference 2 in Figure 11)
There is characteristic of significantly bunching with the Wind turbines of permanent revolving speed invariable power area (appended drawing reference 3 in Figure 11), can identify 7m/s
(starting point in maximal power tracing area) and two cut-points of 12m/s (terminal in maximal power tracing area).
Be not added with it is low wear control module when the phenomenon that it is different, show the difference of active amplitude and its degree of fluctuation compared with
Small, active power curves when different wind speed have very strong 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 restored, it is low in the case where slope sets unanimous circumstances
The Wind turbines of wind speed can quickly reach steady-state value, and the Wind turbines of high wind speed take a long time.
In step 3, further to explore and verifying the selection principle that wind power plant divides group's index, wind power plant emulation mould is built
Type.According to the stronger roll-off characteristic of active power curves, the typical wind speed distribution scene of setting and active opposite systematic error
As shown in table 3, the calculation formula of error is as follows:
Wherein, JstartAnd JendRespectively active power starts to stablize after Voltage Drop previous 1 period and fault clearance
1 period corresponding data segment serial number after output.
Systematic error (%) under the different wind speed scenes of table 3
Assuming that the relative error of active power allows threshold value to be 5%.As shown in Table 3, under wind speed scene 1 and scene 2, single machine
The opposite systematic error of equivalent active power is respectively less than the threshold value set, it was demonstrated that works in promoter region, permanent revolving speed invariable power area
Wind turbines can equivalent be respectively one.
As shown in Table 3, under scene 3, the precision of single machine equivalence is unsatisfactory for threshold requirement.If scene 3 using 7m/s as cut-point,
Precision when by wind power plant equivalence being two machines is doubled left and right compared with single machine, it was demonstrated that with maximal power tracing area beginning wind speed
Make the necessity of cut-point.
Equally, under scene 4, the precision of single machine equivalence is unsatisfactory for threshold requirement, and only bent according to the response of single machine fault power
Line can not identify the cut-point in maximal power tracing region.Since the linearity of maximal power tracing area response curve is preferable,
A large amount of l-G simulation test show with its single machine equivalent wind speed (9.6m/s) it is equivalent be two groups when two machine equivalence effects have compared with single machine
It is apparent to improve, precision is doubled or so.
As shown in Table 3, under scene 5, the error of single machine equivalence seriously exceeds threshold range, punishes in 12m/s into two groups
Two machine equivalence effects improve significantly compared with single machine equivalence, and error is about 4.78%, are closer to threshold requirement.For there are one
Fixed nargin, the study found that the response base of three machine equivalence effects and detailed model when equivalence is three groups at 7m/s, 12m/s
Originally being completely coincident error is about 2.75%.It demonstrates using maximal power tracing area end wind speed as the necessity of cut-point.
The simulation result of the above example is as shown in table 4.Discovery is summarized, if the relative error permission threshold value of active power is
5%, with the starting point in maximal power tracing area, terminal and work in maximal power tracing area Wind turbines single machine equivalence wind speed
Required precision can be met for cut-point.
Cut-point situation under the different wind speed ranges of table 4
It is worth noting that, if the Wind turbines of some producers are in strong wind and small wind (with power output for the corresponding wind of 0.5p.u.
Speed be boundary) when using different power slope recovery policies, then there is natural cut-point, the as line of demarcation of strong wind and small wind
(generally between 9m/s-10m/s), generally by the equivalent wind speed replacement in maximal power tracing area, segmentation wind speed can expire thus
Sufficient required precision.
In step 4, in full blast speed operation area, with temporary when generation three phase short circuit fault at Wind turbines grid entry point
The characteristic of bunching of state active power response curve, the cut-point of acquisition；Active power error threshold value is set according to required precision, according to
The secondary typical examples for attempting wind speed profile, observe the adequacy and necessity of cut-point, obtain cut-point.Take the union of cut-point
Wind turbines equivalence in each group is one, calculated by (7m/s, 9.6m/s, 12m/s) as dividing group's index to divide Wind turbines
The equivalent parameters of each equivalence unit and corresponding current collection network equivalent parameter.Equivalent wind speed is all Wind turbines power outputs in the group
The sum of mean power corresponding to equivalent wind speed.The parameters of equivalent Wind turbines are calculated according to capacity weighting method, equivalent
Cable data is calculated according to the equal principle of current collection network power losses before and after equivalence, and the detailed process present invention will not be described in great detail.
Step 5 carries out Accuracy Verification to equivalence method, guarantees the accuracy of wind power plant Equivalent Model in controlled range
It is interior.Specifically, can according at least one set detect wind speed, respectively to wind power plant Equivalent Model, conventional individual Equivalent Model and in detail
Thin model carries out fault simulation test respectively, compares wind power plant Equivalent Model, conventional individual Equivalent Model and detailed model same
Sample wind speed, the active power under same failure, reactive power, voltage and current response curve, to judge wind power plant equivalence mould
The simulated effect of type.
This example has randomly selected in 2017 group air speed datas of certain 3 × 6 wind power plant May 19 to May 26 such as figure
30 groups of air speed datas shown in 13 carry out equivalent experiment.Position of failure point is wind power plant exit, wind power plant and Infinite bus system
It is connected, main simulation parameter is the same as above.For the equivalent effect for showing the method for the present invention, the 6th group, the 30th group of wind speed have been randomly selected
It is compared with the effect of single machine equivalence, as a result respectively as shown in Figure 14 to Figure 17 and Figure 18 to Figure 21.
By Figure 14 to Figure 17 and Figure 18 to Figure 21 it is found that dividing group's strategy that can significantly improve conventional individual using of the invention
The equivalent precision of polymerization, the active power in wind power plant exit, reactive power, the tracking effect of voltage and current are preferable.Cause
This, equivalence strategy of the invention has good adaptability to air speed data.
For the adaptability for verifying this point of group's strategy on fault point position, fault point has been moved to double circuit power supply line by this example
One time midpoint has re-started simulation calculation to above 30 groups of air speed datas using grouping method above.Still choose
6 groups and the 30th group of simulation result is shown, respectively as shown in Figure 22 to Figure 25 and Figure 26 to Figure 29.
It is of the invention it can be seen from Figure 22 to Figure 25 and Figure 26 to Figure 29 to divide group's strategy that significantly improve conventional individual
Equivalent equivalent precision, the active power in wind power plant exit, reactive power, the tracking effect of voltage and current are preferable.Cause
This, also there is good adaptability to position of failure point.
Specific embodiment 2: present embodiment be to described in specific embodiment one based on the double-fed of active response
Wind power plant Dynamic Equivalence is described further, and in present embodiment, in step 1, the supplemental characteristic for emulating acquisition is wind-powered electricity generation
Voltage, electric current, active power and reactive power at machine set outlet, trouble duration 0.625s；
The supplemental characteristic that the supplemental characteristic for obtaining the emulation in step 1 and on-the-spot test experiment obtain compares, and obtains
The detailed process of the weighted average absolute error of simulation result and field experiment test result are as follows:
Active power response curve and the active power data of on-the-spot test are compared, by reactive power response curve
It is compared with the reactive power data of on-the-spot test, the response curve of voltage and the voltage data of on-the-spot test is carried out pair
Than, the current data of the response curve of electric current and on-the-spot test is compared,
According to formula:
The parameter number of the supplemental characteristic that emulation obtains in period A1, B1, B2, C1 and C2 and on-the-spot test is obtained respectively
According to weighted average absolute error Fii,
In formula, Kend, KstartFor the start sequence number of each time hop counts evidence, XM(i) and Xs(i) it is respectively i-th point of actual measurement
Value and simulation value, F are voltage, electric current, active power or the total weighted average absolute error of reactive power,
By weighted average absolute error F value compared with the accuracy value of standard, the transient state for building Double-feed wind power unit is determined
Model is consistent with the dynamic characteristic of actually measured Wind turbines.
Specific embodiment 3: present embodiment be to described in specific embodiment one based on the double-fed of active response
Wind power plant Dynamic Equivalence is described further, and in present embodiment, the full blast speed operation area in step 2 is incision wind
Speed arrives cut-out wind speed, is divided into 0.1m/s, trouble duration 0.1s between wind speed.
Specific embodiment 4: present embodiment be to described in specific embodiment one based on the double-fed of active response
Wind power plant Dynamic Equivalence is described further, in present embodiment, in step 2, and according to the characteristic of bunching on the curve,
Obtain the detailed process for dividing group's index are as follows:
Blower works at different wind speed regions, and the transient response curve of active power is different, work in promoter region,
The Wind turbines in maximal power tracing area and permanent revolving speed invariable power area have characteristic of significantly bunching, and can recognize that maximum power
The starting point of tracking area and two cut-points of terminal in maximal power tracing area.
Specific embodiment 5: present embodiment be to described in specific embodiment one based on the double-fed of active response
Wind power plant Dynamic Equivalence is described further, in present embodiment, trouble duration 0.1s, and in step 3, according to
The active detailed process with respect to systematic error for dividing group's index to obtain each example in step 2 are as follows:
According to formula:
Obtain the active opposite systematic error E of each examplep,
In formula, JstartAnd JendRespectively active power starts to stablize after Voltage Drop previous 1 period and fault clearance
1 period corresponding data segment serial number after output, XM(i) and Xs(i) it is respectively i-th point of measured value and simulation value.
Specific embodiment 6: present embodiment be to described in specific embodiment one based on the double-fed of active response
Wind power plant Dynamic Equivalence is described further, and in present embodiment, calculates the equivalent parameters of each equivalent unit and corresponding
The detailed process of current collection network equivalent parameter are as follows:
Each parameter of equivalent Wind turbines is calculated according to capacity weighting method, and the parameter of equivalence collection electric network is according to before and after equivalence
The equal principle of current collection network power losses calculates.
Specific embodiment 7: present embodiment be to described in specific embodiment one based on the double-fed of active response
Wind power plant Dynamic Equivalence is described further, in present embodiment, in step 3, and single machine equivalence wind speed veqCalculation formula
Are as follows:
In formula, m is the number of units of Wind turbines, and j is the number of unit, PjFor the active power output of jth platform unit.Pj=f (vj)
Represent wind speed-power characteristic of unit.
Specific embodiment 8: present embodiment be to described in specific embodiment one based on the double-fed of active response
Wind power plant Dynamic Equivalence is described further, further comprising the steps of after step 4 in present embodiment,
The wind power plant Equivalent Model obtained using step 4, using the accurate of the actual wind speed data verification Equivalent Model
Property.
Claims (4)
1. the Double-feed wind power field Dynamic Equivalence based on active response, which is characterized in that it the following steps are included:
Step 1: the electrical-magnetic model of Double-feed wind power unit is built, respectively in the small-power output of the Wind turbines and greatly
Under two kinds of operating conditions of power output, three phase short circuit fault emulation is carried out, the supplemental characteristic which is obtained and on-the-spot test are tested
The supplemental characteristic of acquisition compares, and obtains the weighted average absolute error of simulation result and field experiment test result, and the weighting is flat
Equal absolute error complies with standard required precision, thens follow the steps two；
Step 2: obtain Double-feed wind power unit transient Model in full blast speed operation area occurs three-phase shortcircuit event at grid entry point
Transient state active power response curve when barrier, according to the characteristic of bunching on the curve, acquisition divides group's index；
Step 3: setting active power error threshold value according to practical equivalent required precision, carries out the equivalent of typical wind speed distribution and calculate
Example calculates, and obtains the active opposite systematic error of each example, and the active opposite systematic error of each example is active with setting respectively
Power error threshold value comparison, verifies cut-point, if the active opposite systematic error of some example is greater than setting wattful power
Rate error threshold then adds the single machine equivalence wind speed in corresponding wind speed section as supplement cut-point to the example；
Step 4: taking the union of cut-point to be used as according to the cut-point that step 2 and step 3 obtain offline and group's index being divided to divide
Wind turbines equivalence in each group is one by Wind turbines, calculates the equivalent parameters and corresponding collection power grid of each equivalent unit
Network equivalent parameters obtains wind power plant Equivalent Model；
In step 1, the supplemental characteristic for emulating acquisition is voltage, electric current, active power and the idle function in Wind turbines exit
Rate；
The supplemental characteristic that the supplemental characteristic for obtaining the emulation in step 1 and on-the-spot test experiment obtain compares, and is emulated
As a result with the detailed process of the weighted average absolute error of field experiment test result are as follows:
Active power response curve and the active power data of on-the-spot test are compared, by reactive power response curve and now
The reactive power data of field test compare, and the voltage data of the response curve of voltage and on-the-spot test is compared, will
The response curve of electric current and the current data of on-the-spot test compare,
According to formula:
The supplemental characteristic of the supplemental characteristic that emulation obtains in period A1, B1, B2, C1 and C2 and on-the-spot test is obtained respectively
It is weighted and averaged absolute error Fii,
In formula, Kend, KstartFor the start sequence number of each time hop counts evidence, XM(i) and Xs(i) be respectively i-th point measured value and
Simulation value, F are voltage, electric current, active power or the total weighted average absolute error of reactive power, FA1To be imitated in period A1
The weighted average absolute error of the supplemental characteristic of the supplemental characteristic and on-the-spot test that really obtain；FB1It is obtained to be emulated in period B1
The weighted average absolute error of the supplemental characteristic of the supplemental characteristic and on-the-spot test that obtain；FB2It is obtained for emulation in period B2
The weighted average absolute error of this of supplemental characteristic and on-the-spot test supplemental characteristic；FC1For the parameter for emulating acquisition in period C1
The weighted average absolute error of this of data and on-the-spot test supplemental characteristic；FC2For the supplemental characteristic for emulating acquisition in period C2
With the weighted average absolute error of the supplemental characteristic of on-the-spot test；
By weighted average absolute error F value compared with the accuracy value of standard, the transient Model for building Double-feed wind power unit is determined
It is consistent with the dynamic characteristic of actually measured Wind turbines.
2. the Double-feed wind power field Dynamic Equivalence according to claim 1 based on active response, which is characterized in that step
Full blast speed operation area in rapid two is incision wind speed to cut-out wind speed, is divided into 0.1m/s between wind speed.
3. the Double-feed wind power field Dynamic Equivalence according to claim 1 based on active response, which is characterized in that step
In rapid two, according to the characteristic of bunching on the curve, the detailed process for dividing group's index is obtained are as follows:
Blower works at different wind speed regions, and the transient response curve of active power is different, works in promoter region, maximum
The Wind turbines in power tracking area and permanent revolving speed invariable power area have characteristic of significantly bunching, and can recognize that maximal power tracing
The starting point wind speed in area and two cut-points of terminal wind speed in maximal power tracing area.
4. the Double-feed wind power field Dynamic Equivalence according to claim 1 based on active response, which is characterized in that step
In rapid three, the detailed process of the active opposite systematic error of each example is obtained are as follows:
According to formula:
Obtain the active opposite systematic error E of each examplep,
In formula, JstartAnd JendRespectively active power starts to stablize output after Voltage Drop previous 1 period and fault clearance
1 period corresponding data segment serial number afterwards, XM(i) and Xs(i) it is respectively i-th point of measured value and simulation value.
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 CN106202815A (en) | 2016-12-07 |
| CN106202815B true 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) |
Families Citing this family (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 |
| CN107272643B (en) * | 2017-06-19 | 2019-03-29 | 哈尔滨工业大学 | The single machine equivalence method of new energy station |
| 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 |
| CN110032778A (en) * | 2019-03-26 | 2019-07-19 | 河海大学 | A kind of wind power plant single machine equivalent modeling verification method |
| CN110119570B (en) * | 2019-05-10 | 2023-08-29 | 东北电力大学 | Actually measured data driven wind farm model parameter checking method |
| CN110175933B (en) * | 2019-05-29 | 2021-06-22 | 哈尔滨工业大学 | Direct-drive wind power plant dynamic equivalence method suitable for frequency modulation control |
| CN113627035B (en) * | 2021-09-01 | 2024-01-19 | 国网黑龙江省电力有限公司电力科学研究院 | Method for establishing single-machine equivalent model of wind power plant |
| CN114188943B (en) * | 2021-12-15 | 2023-09-08 | 河海大学 | Dynamic equivalent modeling method and system for photovoltaic field group |
| CN114357787B (en) * | 2022-01-10 | 2024-05-28 | 华北电力大学 | Offshore wind farm equivalent modeling method and system |
| CN115882526B (en) * | 2023-02-13 | 2023-05-05 | 华北电力科学研究院有限责任公司 | Active and reactive power regulation capability data processing method and device facing wind power plant group |
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 |
|---|
| Investigation on Fundamental Issues of Dynamic Equivalent Modeling of Wind Farms;Pupu Chao等;《Preprints of the 5th international conference on Electric Utility Deregulation and Restructuring and Power Technologies》;20151129;第192-196页 |
| 基于聚类–判别分析的风电场概率等值建模研究;朱乾龙等;《中国电机工程学报》;20141005;第34卷(第28期);第4770-4780页 |
| 永磁直驱同步风电场多机动态等值模型;蒙晓航;《电力系统保护与控制》;20130716;第41卷(第14期);第25-32页 |
| 风电场动态等值建模方法研究;张洋溢;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20160215;第2016年卷(第02期);第C042-623页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106202815A (en) | 2016-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106202815B (en) | Double-feed wind power field Dynamic Equivalence based on active response | |
| Larsson | Flicker emission of wind turbines during continuous operation | |
| CN107453401B (en) | Doubly-fed wind generator parameter identification method | |
| 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 | |
| CN106058937B (en) | A kind of power distribution network broad sense load modeling method of the wind power plant containing direct-drive permanent-magnetism | |
| CN103700036B (en) | A kind of transient stability projecting integral method being suitable to power system Multiple Time Scales | |
| CN102819221B (en) | LVRT Capability of Wind Turbine Generator characteristic associative simulation model and joint simulation method thereof | |
| CN104392056B (en) | A kind of directly driven wind-powered unit modeling method and device | |
| CN105678033B (en) | A kind of wind-powered electricity generation station equivalent modeling method suitable for electromagnetic transient simulation | |
| CN103219725A (en) | Wind power plant equivalent modeling method based on real-time operation data | |
| CN105741188B (en) | Doubly-fed wind turbine group's electro-magnetic transient coherence grouping method | |
| CN105576654B (en) | Directly driven wind-powered field equivalence method and system | |
| Mahmoud et al. | A developed integrated scheme based approach for wind turbine intelligent control | |
| CN109522607A (en) | A kind of double-fed fan motor field electromechanical transient equivalent modeling method | |
| CN106526347A (en) | Digital-analog hybrid simulation-based photovoltaic inverter low voltage ride through evaluation method | |
| CN108845230A (en) | A kind of sub-synchronous oscillation random time-dependent modal identification method | |
| CN104820741A (en) | Wind power plant dynamic equivalence method combining wind field disperstiveness and unit difference | |
| CN102830692A (en) | Low-voltage ride through consistency testing system and method of main control system of wind turbine generator | |
| CN109190320A (en) | Parallel Heterogeneous simulation method suitable for ADPSS double-fed blower closed-loop test | |
| CN109921426A (en) | Wind-electricity integration system probability load flow calculation method based on CV-KDE | |
| 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 | |
| CN104538979B (en) | The emulation modelling method that wind energy turbine set dynamic passive compensation capacity optimizes |
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 |