CN117665445A - New energy unit electromechanical model checking method based on whole group fault ride-through test - Google Patents
New energy unit electromechanical model checking method based on whole group fault ride-through test Download PDFInfo
- Publication number
- CN117665445A CN117665445A CN202311629258.3A CN202311629258A CN117665445A CN 117665445 A CN117665445 A CN 117665445A CN 202311629258 A CN202311629258 A CN 202311629258A CN 117665445 A CN117665445 A CN 117665445A
- Authority
- CN
- China
- Prior art keywords
- voltage
- disturbance
- deviation
- whole group
- new energy
- 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.)
- Pending
Links
- 238000012360 testing method Methods 0.000 title claims abstract description 83
- 238000000034 method Methods 0.000 title claims abstract description 71
- 238000004088 simulation Methods 0.000 claims abstract description 39
- 238000004364 calculation method Methods 0.000 claims abstract description 19
- 230000001052 transient effect Effects 0.000 claims description 41
- 230000011218 segmentation Effects 0.000 claims description 8
- 238000011084 recovery Methods 0.000 claims description 3
- 230000000630 rising effect Effects 0.000 claims description 3
- 238000012795 verification Methods 0.000 description 5
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention relates to a new energy unit electromechanical model checking method based on a whole group of fault ride-through test, and belongs to the technical field of power system simulation. The invention utilizes the whole group of fault ride-through curves of high voltage and low voltage output by the power electronic fault ride-through test equipment to carry out the whole group of fault ride-through test on the new energy unit, calculates the voltage, the active power, the reactive power, the active current and the reactive current of the fundamental wave positive sequence component by utilizing the full-wave Fourier coefficient based on the test result, carries out deviation calculation with model simulation data, and judges that the electromechanical model of the new energy unit meets the requirement if the deviation is within the allowable range; otherwise, adjusting the model parameters, and checking again; if the model parameters are adjusted for multiple times and the deviation requirement cannot be completely met all the time, selecting a group of optimal model parameters as a final result. The invention can reduce the workload of checking the new energy unit and accelerate the checking speed; a set of optimal model parameters can be screened, and the method is convenient to use in the stable calculation of the power system.
Description
Technical Field
The invention belongs to the technical field of power system simulation, and particularly relates to a new energy unit electromechanical model checking method based on a whole group of fault ride-through tests.
Background
Under the background of realizing ' carbon peak, carbon neutralization ' and ' target, pushing to construct a novel power system taking new energy as a main body, new large-scale wind power and photovoltaic grid connection are added, the wind power and the photovoltaic occupy higher and higher installed proportion of the power system, and the stability influence analysis of the power system of the wind power and the photovoltaic is not negligible.
GB 38755-2019 "safety and stability guidance of electric power system" requires that electric power system stability calculation analysis should be performed according to specific conditions and requirements of the system, but most electric power system stability calculation analysis is not added with models of wind power plants and photovoltaic power stations at present. On the one hand, as the model checking speeds of a wind power plant and a photovoltaic power station are slower, according to the requirements of the existing standards GB/T32826-2016 on the modeling guideline of a photovoltaic power generation system, GB/T32892-2016 on the model of the photovoltaic power generation system, the test of a parameter test procedure of an electrical simulation model of a wind turbine generator, NB/T31053-2021 on the verification procedure of the electrical simulation model of the wind turbine generator and the like, the new energy machine set model checking of one-by-one voltage working condition (0% Un, 20% Un, 35% Un, 50% Un, 75% Un, 115% Un, 120% Un, 125% Un and 130% Un) needs to be carried out, the checking data source test is slower, and the checking workload is large; on the other hand, the traditional electromechanical transient simulation supports the identification of four sets of model parameters (large-load symmetrical faults, large-load asymmetrical faults, small-load symmetrical faults and small-load asymmetrical fault control strategy parameters), but because 'black box models' of wind turbines and photovoltaic array inverters are required to adopt different model parameters aiming at different fault working conditions, the four sets of parameters can not completely meet deviation requirements, and different model parameters need to be manually adjusted to form different stable calculation simulations, so that the stable calculation of a power grid cannot be used.
Therefore, how to overcome the defects of the prior art is a problem to be solved in the current power system simulation technical field.
Disclosure of Invention
The invention aims to solve the defects of the prior art and provides a new energy unit electromechanical model checking method based on a whole group of fault ride-through tests. The invention perfects the existing new energy unit model checking method, improves checking speed, screens an optimal model parameter and is easy to popularize and apply.
In order to achieve the above purpose, the technical scheme adopted by the invention is as follows:
the new energy unit electromechanical model checking method based on the whole group of fault ride-through tests comprises the following steps:
step (1), a whole group of fault ride-through curves of high voltage and low voltage are output by using power electronic fault ride-through test equipment, and a whole group of fault ride-through tests are carried out on a new energy unit;
step (2), calculating voltage, active power, reactive power, active current and reactive current of a fundamental wave positive sequence component by using a full-wave Fourier coefficient based on the new energy unit test result obtained in the step (1);
step (3), carrying out per unit on the calculation result in the step (2), and segmenting the per unit data to obtain segmented test data; the specific segmentation method comprises the following steps: segmenting data according to three periods of time before disturbance, during disturbance and after disturbance;
step (4), building a new energy unit electromechanical transient simulation model in the electromechanical transient simulation software, and obtaining a whole group of fault ride-through simulation test results of the new energy unit electromechanical transient simulation model at high voltage and low voltage; the simulation test result comprises voltage, active power, reactive power, active current and reactive current of fundamental wave positive sequence component
Step (5), carrying out per unit on the simulation test result in the step (4), and segmenting the per unit data to obtain segmented simulation data; the specific segmentation method comprises the following steps: segmenting data according to three periods of time before disturbance, during disturbance and after disturbance;
step (6), performing deviation calculation on the test data segmented in the step (3) and the simulation data corresponding to the segmented test data segmented in the step (5);
step (7), if all the deviations calculated in the step (6) are within the allowable range, judging that the new energy unit electromechanical model meets the requirements; otherwise, readjusting and adjusting the model parameters, and repeating the checking process from the step (4) to the step (7).
Further, it is preferable that the whole set of fault ride-through curves includes a high voltage whole set of ride-through curves including 120% un, 125% un, 130% un voltage rise points, and a low voltage whole set of ride-through curves including 20% un, 35% un, 50% un, 75% un voltage drop points; where Un represents the rated voltage.
Further, it is preferable that the full-group fault ride-through test includes 8 working condition tests of a high-voltage full-group traversing large-load symmetrical fault, a large-load asymmetrical fault, a small-load symmetrical fault, a small-load asymmetrical fault, and a low-voltage full-group traversing large-load symmetrical fault, a large-load asymmetrical fault, a small-load symmetrical fault, and a small-load asymmetrical fault, respectively; the large load is that the output active power of the unit is larger than 0.7 rated active power; the small load is the rated active power of which the unit output active power is more than or equal to 0.1 and less than or equal to 0.3.
Further, it is preferable that in the step (3) and the step (5), the disturbance period is divided into a disturbance period steady-state section and a disturbance period steady-state section when the disturbance period is segmented; the post-disturbance steady-state section and the post-disturbance transient-state section are divided into the following concrete:
the voltage drops to a rated voltage value of 0.9 or 2 seconds before the voltage rises to a rated voltage of 1.1 is taken as a starting point of a section before disturbance; the 20 milliseconds before the voltage drops to 0.9 rated voltage value or the voltage rises to 1.1 rated voltage is the end point of the section before disturbance and is also the start point of the transient section during disturbance;
the transient state interval end point in the disturbance period is 40ms after the transient state interval end point is reduced and the active current and the reactive current are stabilized, and is also the steady state interval start point in the disturbance period;
the 20 milliseconds before the normal time of the voltage recovery is the end point of the steady-state section in the disturbance period and is also the start point of the transient-state section after the disturbance;
the end point of the transient state section after disturbance is 20 milliseconds after the voltage is recovered to be normal and the active power of the new energy unit is recovered to the point before the fault, and is also the start point of the steady state section after disturbance;
and the end point of the steady-state interval after disturbance is the end of a period of 2 seconds after steady-state output of the active power.
Further, it is preferable that the deviation type includes average deviation, average absolute deviation, maximum deviation, and weighted average absolute deviation;
calculating average deviation, average absolute deviation and maximum deviation in a steady-state interval before disturbance;
calculating average deviation, average absolute deviation and maximum deviation in a steady state period in a disturbance period, and calculating average deviation and average absolute deviation in a transient state in the disturbance period;
calculating average deviation, average absolute deviation and maximum deviation of the steady-state interval after disturbance, and calculating average deviation and average absolute deviation of the transient-state interval after disturbance;
and carrying out weighted average on the average absolute deviation of each disturbance interval to obtain the weighted average absolute deviation of the whole disturbance process.
Further, it is preferable that the specific method for calculating the voltage, active power, reactive power, active current and reactive current of the fundamental wave positive sequence component using the full wave fourier coefficient is:
firstly, calculating Fourier coefficients of fundamental wave components in a fundamental wave period;
wherein: f (f) 1 Is the fundamental frequency, namely 50 Hz; u (u) a,cos 、u b,cos 、u c,cos The voltage fundamental wave cosine components of a phase, b phase and c phase are respectively; u (u) a,sin 、u b,sin 、u c,sin The voltage fundamental wave sinusoidal components of a phase, b phase and c phase; i.e a,cos 、i b,cos 、i c,cos The current fundamental wave cosine components of a phase, b phase and c phase are respectively; i.e a,sin 、i b,sin 、i c,sin The current fundamental wave sinusoidal components of a phase, b phase and c phase; u (u) a 、u b 、u c The voltages of a phase, b phase and c phase; i.e a 、i b 、i c The current is a phase current, b phase current and c phase current; t is the fundamental wave period; t is a time variable;
the voltage and current vector components of the fundamental positive sequence component are calculated using:
wherein: u (u) 1+,cos Cosine component of positive sequence of fundamental wave voltage; u (u) 1+,sin Sinusoidal components of the fundamental voltage positive sequence; i.e 1+,cos Cosine component of positive sequence of fundamental wave voltage; i.e 1+,sin Sinusoidal components of the fundamental voltage positive sequence;
voltage U of fundamental positive sequence component 1+ The method comprises the following steps:
active power P of fundamental positive sequence component 1+ The method comprises the following steps:
reactive power Q of fundamental positive sequence component 1+ The method comprises the following steps:
active current I of fundamental positive sequence component P1+ The method comprises the following steps:
reactive current I of fundamental positive sequence component Q1+ The method comprises the following steps:
further, it is preferable that, if the deviation allowance requirement is not completely satisfied all the time after the model parameters are adjusted a plurality of times, a set of model parameters having the smallest weighted average absolute deviation is selected as the final result.
Further, it is preferable that the number of times is not less than ten.
The power electronic fault ride-through test equipment is existing equipment.
In the invention, the method for identifying the parameters of the electromechanical transient simulation model of the new energy unit adopts the existing mature identification method, and the invention does not have special requirements on the method.
The invention divides the per-unit test data and the per-unit simulation data sequence into three periods before disturbance, during disturbance and after disturbance; dividing the disturbance period and the disturbance period into a transient state section and a steady state section according to the response characteristics of the active power and the reactive power,
the specific calculation of the deviation of the invention is shown in table 1;
TABLE 1
And carrying out weighted average on the average absolute deviation of each disturbance interval to obtain the weighted average absolute deviation of the whole disturbance process.
Wherein the calculation of the average deviation, the average absolute deviation, the maximum deviation and the weighted average absolute deviation refers to the NB/T31053-2021 standard.
In the invention, the test data segmented in the step (3) and the simulation data corresponding to the segmented test data in the step (5) are subjected to deviation calculation; the correspondence here means that the working conditions correspond, the segments correspond, and the data correspond.
In step (7) of the present invention, all deviations are within the allowable range, i.e., all deviation results meet the NB/T31053-2021 standard.
Compared with the prior art, the invention has the beneficial effects that:
(1) The invention reduces the workload of checking the electromechanical model of the new energy unit and accelerates the checking speed; the existing checking at least checks 32 high and low voltage ride through test conditions, and the method only needs to check 8 high and low voltage ride through test conditions, so that the checking efficiency is greatly improved.
(2) The method is characterized in that a plurality of new energy units adopt different control strategies under different load working conditions and at different voltage drop or rise moments, and the electromechanical transient simulation software needs to adopt different parameters for fitting; the method can screen out the optimal set of model parameters under each fault and load working condition, and meets the use requirement of the existing electromechanical transient simulation software.
Drawings
In order to more clearly illustrate the technical solutions of the present application, the drawings that are needed in the embodiments will be briefly described below, and it will be obvious to those skilled in the art that other drawings can be obtained from these drawings without inventive effort.
FIG. 1 is a flow chart of a new energy unit electromechanical model checking method based on a whole group of fault ride-through tests;
FIG. 2 is a schematic diagram of a high voltage full-set ride-through test curve provided by the present invention;
FIG. 3 is a schematic diagram of a low voltage full-set ride-through test curve provided by the present invention;
fig. 4 is a graph of the high voltage whole group crossing test result of a certain photovoltaic array test provided in application example 1 of the present invention;
FIG. 5 is a graph comparing the results of a high voltage full-set ride-through test and an electromechanical simulation for a photovoltaic array test provided in application example 1 of the present invention;
FIG. 6 is a graph of the results of a high voltage full-set ride-through test for a photovoltaic array test provided in application example 2 of the present invention;
fig. 7 is a graph comparing the results of the high voltage full set of pass through tests and the electromechanical simulation of a certain photovoltaic array test provided in application example 2 of the present invention.
Detailed Description
The present invention will be described in further detail with reference to examples.
It will be appreciated by those skilled in the art that the following examples are illustrative of the present invention and should not be construed as limiting the scope of the invention. The specific techniques or conditions are not identified in the examples and are performed according to techniques or conditions described in the literature in this field or according to the product specifications. The materials or equipment used are conventional products available from commercial sources, not identified to the manufacturer.
The new energy unit electromechanical model checking method based on the whole group of fault ride-through tests comprises the following steps:
step (1), a whole group of fault ride-through curves of high voltage and low voltage are output by using power electronic fault ride-through test equipment, and a whole group of fault ride-through tests are carried out on a new energy unit;
step (2), calculating voltage, active power, reactive power, active current and reactive current of a fundamental wave positive sequence component by using a full-wave Fourier coefficient based on the new energy unit test result obtained in the step (1);
step (3), carrying out per unit on the calculation result in the step (2), and segmenting the per unit data to obtain segmented test data; the specific segmentation method comprises the following steps: segmenting data according to three periods of time before disturbance, during disturbance and after disturbance;
step (4), building a new energy unit electromechanical transient simulation model in the electromechanical transient simulation software, and obtaining a whole group of fault ride-through simulation test results of the new energy unit electromechanical transient simulation model at high voltage and low voltage; the simulation test result comprises voltage, active power, reactive power, active current and reactive current of fundamental wave positive sequence component
Step (5), carrying out per unit on the simulation test result in the step (4), and segmenting the per unit data to obtain segmented simulation data; the specific segmentation method comprises the following steps: segmenting data according to three periods of time before disturbance, during disturbance and after disturbance;
step (6), performing deviation calculation on the test data segmented in the step (3) and the simulation data corresponding to the segmented test data segmented in the step (5);
step (7), if all the deviations calculated in the step (6) are within the allowable range, judging that the new energy unit electromechanical model meets the requirements; otherwise, readjusting and adjusting the model parameters, and repeating the checking process from the step (4) to the step (7).
Specifically, the whole group of fault ride-through curves comprises a high-voltage whole group of ride-through curves and a low-voltage whole group of ride-through curves, wherein the high-voltage whole group of ride-through curves comprise 120% Un, 125% Un and 130% Un voltage rising points, and the low-voltage ride-through curves comprise 20% Un, 35% Un, 50% Un and 75% Un voltage falling points; where Un represents the rated voltage.
Specifically, the whole group fault ride-through test respectively comprises 8 working condition tests of high-voltage whole group passing large-load symmetrical faults, high-load asymmetrical faults, small-load symmetrical faults and small-load asymmetrical faults, and low-voltage whole group passing large-load symmetrical faults, large-load asymmetrical faults, small-load symmetrical faults and small-load asymmetrical faults; the large load is that the output active power of the unit is larger than 0.7 rated active power; the small load is the rated active power of which the unit output active power is more than or equal to 0.1 and less than or equal to 0.3.
Specifically, in the step (3) and the step (5), during segmentation, the disturbance period is divided into a stable period of the disturbance period and a stable period of the disturbance period; the post-disturbance steady-state section and the post-disturbance transient-state section are divided into the following concrete:
the voltage drops to a rated voltage value of 0.9 or 2 seconds before the voltage rises to a rated voltage of 1.1 is taken as a starting point of a section before disturbance; the 20 milliseconds before the voltage drops to 0.9 rated voltage value or the voltage rises to 1.1 rated voltage is the end point of the section before disturbance and is also the start point of the transient section during disturbance;
the transient state interval end point in the disturbance period is 40ms after the transient state interval end point is reduced and the active current and the reactive current are stabilized, and is also the steady state interval start point in the disturbance period;
the 20 milliseconds before the normal time of the voltage recovery is the end point of the steady-state section in the disturbance period and is also the start point of the transient-state section after the disturbance;
the end point of the transient state section after disturbance is 20 milliseconds after the voltage is recovered to be normal and the active power of the new energy unit is recovered to the point before the fault, and is also the start point of the steady state section after disturbance;
and the end point of the steady-state interval after disturbance is the end of a period of 2 seconds after steady-state output of the active power.
Specifically, the deviation types include average deviation, average absolute deviation, maximum deviation and weighted average absolute deviation;
calculating average deviation, average absolute deviation and maximum deviation in a steady-state interval before disturbance;
calculating average deviation, average absolute deviation and maximum deviation in a steady state period in a disturbance period, and calculating average deviation and average absolute deviation in a transient state in the disturbance period;
calculating average deviation, average absolute deviation and maximum deviation of the steady-state interval after disturbance, and calculating average deviation and average absolute deviation of the transient-state interval after disturbance;
and carrying out weighted average on the average absolute deviation of each disturbance interval to obtain the weighted average absolute deviation of the whole disturbance process.
Specifically, the specific method for calculating the voltage, the active power, the reactive power, the active current and the reactive current of the fundamental wave positive sequence component by using the full-wave Fourier coefficient is as follows:
firstly, calculating Fourier coefficients of fundamental wave components in a fundamental wave period;
wherein: f (f) 1 Is the fundamental frequency, namely 50 Hz; u (u) a,cos 、u b,cos 、u c,cos The voltage fundamental wave cosine components of a phase, b phase and c phase are respectively; u (u) a,sin 、u b,sin 、u c,sin The voltage fundamental wave sinusoidal components of a phase, b phase and c phase; i.e a,cos 、i b,cos 、i c,cos The current fundamental wave cosine components of a phase, b phase and c phase are respectively; i.e a,sin 、i b,sin 、i c,sin The current fundamental wave sinusoidal components of a phase, b phase and c phase; u (u) a 、u b 、u c The voltages of a phase, b phase and c phase; i.e a 、i b 、i c The current is a phase current, b phase current and c phase current; t is the fundamental wave period; t is a time variable;
the voltage and current vector components of the fundamental positive sequence component are calculated using:
wherein: u (u) 1+,cos Cosine component of positive sequence of fundamental wave voltage; u (u) 1+,sin Sinusoidal components of the fundamental voltage positive sequence; i.e 1+,cos Cosine component of positive sequence of fundamental wave voltage; i.e 1+,sin Sinusoidal components of the fundamental voltage positive sequence;
voltage U of fundamental positive sequence component 1+ The method comprises the following steps:
active power P of fundamental positive sequence component 1+ The method comprises the following steps:
reactive power Q of fundamental positive sequence component 1+ The method comprises the following steps:
active current I of fundamental positive sequence component P1+ The method comprises the following steps:
reactive current I of fundamental positive sequence component Q1+ The method comprises the following steps:
specifically, if the model parameters are adjusted multiple times and the deviation allowance requirement cannot be met completely all the time, a group of model parameters with the smallest weighted average absolute deviation is selected as a final result.
Specifically, the number of times is not less than ten.
Application example 1
The voltage, active power, reactive power, active current and reactive current of the fundamental wave positive sequence component of the high-voltage whole group crossing test of the field large-load symmetrical fault of a certain 3.15MW photovoltaic array are shown in fig. 4, and the curves are subjected to per unit treatment, specifically: the rated active power of the photovoltaic array is 3.15MW, and the voltage class of the whole group of high voltage penetrating through the test points is 35kV, so that the voltage, active power, reactive power, active current and reactive current curves of the fundamental wave positive sequence component are subjected to per unit treatment based on the rated voltage of 35kV, the rated active power of 3.15MW and the rated current of 51.96A; the photovoltaic array high voltage ride through result simulated by the new energy unit electromechanical transient simulation model in the electromechanical transient simulation software is subjected to per unit processing, and specifically comprises the following steps: the fundamental wave positive sequence component voltage, active current and reactive current obtained by the electromechanical transient software are per unit values, and the fundamental wave positive sequence component active power and reactive power which are simulated according to the rated active power of 3.15MW are subjected to per unit treatment, as shown in fig. 5.
Dividing the test data after per unit and the simulation data sequence after per unit into three periods before disturbance, during disturbance and after disturbance; according to the response characteristics of the active power and the reactive power, the disturbance period and the post-disturbance period are divided into a transient state section and a steady state section, and the sections divided by the dotted line in fig. 5 are a pre-disturbance period transient state section, a disturbance period steady state section, a post-disturbance transient state section and a post-disturbance steady state section in sequence. And respectively calculating deviation of voltage, active power, reactive power, active current and reactive current of fundamental wave positive sequence components of field actual measurement test data and simulation data, wherein after the electromechanical transient model parameters are adjusted for ten times, deviation calculation results are shown in a table 1, and most deviation is within a meeting range, namely the electromechanical model parameters of the photovoltaic array meet requirements. The verification of the high-voltage whole group crossing test of the large-load asymmetrical faults, the small-load symmetrical faults and the small-load asymmetrical faults is the same as that of the verification, and the repeated description is omitted.
Assuming that the existing method is adopted, 4 sets of electromechanical transient model parameters are needed in order to ensure that the deviation meets the NB/T31053-2021 standard requirement, but the voltage rise amplitude is estimated in advance in the power grid stability calculation, and then the proper model parameters are switched manually, so that the method is inconvenient to use and the estimated voltage rise amplitude is difficult; the method provided by the invention has the advantages that although the deviation of reactive power and reactive current can not meet the standard requirement, only one set of model parameters are adopted as the optimal result, and the auxiliary voltage rising and manual switching model parameters are not needed to be estimated in the power grid stability calculation, so that the method is convenient to use.
TABLE 1 high Voltage Whole set penetration test heavy load, symmetrical Fault check results
Application example 2
The voltage, active power, reactive power, active current and reactive current of the fundamental wave positive sequence component of the low-voltage whole group crossing test of the field large-load symmetrical fault of a certain 3.15MW photovoltaic array are subjected to per unit treatment on the curve as shown in fig. 6; and the whole group of low-voltage traversing test results of the photovoltaic array simulated by the new energy unit electromechanical model in the electromechanical transient simulation software are subjected to per unit processing, as shown in fig. 7.
Dividing the test data after per unit and the simulation data sequence after per unit into three periods before disturbance, during disturbance and after disturbance; according to the response characteristics of the active power and the reactive power, the disturbance period and the post-disturbance period are divided into a transient state section and a steady state section, and the sections divided by the dotted line in fig. 7 are a pre-disturbance period transient state section, a disturbance period steady state section, a post-disturbance transient state section and a post-disturbance steady state section in sequence. And respectively calculating deviation of voltage, active power, reactive power, active current and reactive current of fundamental wave positive sequence components of field actual measurement test data and simulation data, wherein the deviation calculation results are shown in table 2, and the deviation is in a meeting range, namely the electromechanical model parameters of the photovoltaic array meet the requirements. The verification of the high-voltage whole group crossing test of the large-load asymmetrical faults, the small-load symmetrical faults and the small-load asymmetrical faults is the same as that of the verification, and the repeated description is omitted.
The whole group of fault ride-through test checking method of the wind turbine generator is consistent with the photovoltaic array checking method.
TABLE 2 high load and symmetrical fault check results for low voltage full set ride through test
The foregoing has shown and described the basic principles, principal features and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, and that the above embodiments and descriptions are merely illustrative of the principles of the present invention, and various changes and modifications may be made without departing from the spirit and scope of the invention, which is defined in the appended claims. The scope of the invention is defined by the appended claims and equivalents thereof.
Claims (8)
1. The new energy unit electromechanical model checking method based on the whole group of fault ride-through tests is characterized by comprising the following steps:
step (1), a whole group of fault ride-through curves of high voltage and low voltage are output by using power electronic fault ride-through test equipment, and a whole group of fault ride-through tests are carried out on a new energy unit;
step (2), calculating voltage, active power, reactive power, active current and reactive current of a fundamental wave positive sequence component by using a full-wave Fourier coefficient based on the new energy unit test result obtained in the step (1);
step (3), carrying out per unit on the calculation result in the step (2), and segmenting the per unit data to obtain segmented test data; the specific segmentation method comprises the following steps: segmenting data according to three periods of time before disturbance, during disturbance and after disturbance;
step (4), building a new energy unit electromechanical transient simulation model in the electromechanical transient simulation software, and obtaining a whole group of fault ride-through simulation test results of the new energy unit electromechanical transient simulation model at high voltage and low voltage; the simulation test result comprises voltage, active power, reactive power, active current and reactive current of fundamental wave positive sequence component
Step (5), carrying out per unit on the simulation test result in the step (4), and segmenting the per unit data to obtain segmented simulation data; the specific segmentation method comprises the following steps: segmenting data according to three periods of time before disturbance, during disturbance and after disturbance;
step (6), performing deviation calculation on the test data segmented in the step (3) and the simulation data corresponding to the segmented test data segmented in the step (5);
step (7), if all the deviations calculated in the step (6) are within the allowable range, judging that the new energy unit electromechanical model meets the requirements; otherwise, readjusting and adjusting the model parameters, and repeating the checking process from the step (4) to the step (7).
2. The new energy unit electromechanical model checking method based on the whole group of fault ride through tests as claimed in claim 1, wherein the method is characterized by comprising the following steps: the whole group of fault ride-through curves comprise a high-voltage whole group of ride-through curves and a low-voltage whole group of ride-through curves, the high-voltage whole group of ride-through curves comprise 120% Un, 125% Un and 130% Un voltage rising points, and the low-voltage ride-through curves comprise 20% Un, 35% Un, 50% Un and 75% Un voltage falling points; where Un represents the rated voltage.
3. The new energy unit electromechanical model checking method based on the whole group of fault ride through tests as claimed in claim 1, wherein the method is characterized by comprising the following steps: the whole group fault ride-through test comprises 8 working condition tests of high-voltage whole group passing large-load symmetrical faults, high-load asymmetrical faults, small-load symmetrical faults and small-load asymmetrical faults, and low-voltage whole group passing large-load symmetrical faults, high-load asymmetrical faults, small-load symmetrical faults and small-load asymmetrical faults respectively; the large load is that the output active power of the unit is larger than 0.7 rated active power; the small load is the rated active power of which the unit output active power is more than or equal to 0.1 and less than or equal to 0.3.
4. The new energy unit electromechanical model checking method based on the whole group of fault ride through tests as claimed in claim 1, wherein the method is characterized by comprising the following steps: in the step (3) and the step (5), during segmentation, the disturbance period is divided into a disturbance period steady-state interval and a disturbance period steady-state interval; the post-disturbance steady-state section and the post-disturbance transient-state section are divided into the following concrete:
the voltage drops to a rated voltage value of 0.9 or 2 seconds before the voltage rises to a rated voltage of 1.1 is taken as a starting point of a section before disturbance; the 20 milliseconds before the voltage drops to 0.9 rated voltage value or the voltage rises to 1.1 rated voltage is the end point of the section before disturbance and is also the start point of the transient section during disturbance;
the transient state interval end point in the disturbance period is 40ms after the transient state interval end point is reduced and the active current and the reactive current are stabilized, and is also the steady state interval start point in the disturbance period;
the 20 milliseconds before the normal time of the voltage recovery is the end point of the steady-state section in the disturbance period and is also the start point of the transient-state section after the disturbance;
the end point of the transient state section after disturbance is 20 milliseconds after the voltage is recovered to be normal and the active power of the new energy unit is recovered to the point before the fault, and is also the start point of the steady state section after disturbance;
and the end point of the steady-state interval after disturbance is the end of a period of 2 seconds after steady-state output of the active power.
5. The new energy unit electromechanical model checking method based on the whole group of fault ride through tests as claimed in claim 1, wherein the method is characterized by comprising the following steps:
the deviation type comprises average deviation, average absolute deviation, maximum deviation and weighted average absolute deviation;
calculating average deviation, average absolute deviation and maximum deviation in a steady-state interval before disturbance;
calculating average deviation, average absolute deviation and maximum deviation in a steady state period in a disturbance period, and calculating average deviation and average absolute deviation in a transient state in the disturbance period;
calculating average deviation, average absolute deviation and maximum deviation of the steady-state interval after disturbance, and calculating average deviation and average absolute deviation of the transient-state interval after disturbance;
and carrying out weighted average on the average absolute deviation of each disturbance interval to obtain the weighted average absolute deviation of the whole disturbance process.
6. The new energy unit electromechanical model checking method based on the whole group of fault ride through tests as claimed in claim 1, wherein the method is characterized by comprising the following steps: the specific method for calculating the voltage, the active power, the reactive power, the active current and the reactive current of the fundamental wave positive sequence component by using the full wave Fourier coefficient comprises the following steps:
firstly, calculating Fourier coefficients of fundamental wave components in a fundamental wave period;
wherein: f (f) 1 Is the fundamental frequency, namely 50 Hz; u (u) a,cos 、u b,cos 、u c,cos The voltage fundamental wave cosine components of a phase, b phase and c phase are respectively; u (u) a,sin 、u b,sin 、u c,sin For phase aFundamental sinusoidal components of b-phase and c-phase voltages; i.e a,cos 、i b,cos 、i c,cos The current fundamental wave cosine components of a phase, b phase and c phase are respectively; i.e a,sin 、i b,sin 、i c,sin The current fundamental wave sinusoidal components of a phase, b phase and c phase; u (u) a 、u b 、u c The voltages of a phase, b phase and c phase; i.e a 、i b 、i c The current is a phase current, b phase current and c phase current; t is the fundamental wave period; t is a time variable;
the voltage and current vector components of the fundamental positive sequence component are calculated using:
wherein: u (u) 1+,cos Cosine component of positive sequence of fundamental wave voltage; u (u) 1+,sin Sinusoidal components of the fundamental voltage positive sequence; i.e 1+,cos Cosine component of positive sequence of fundamental wave voltage; i.e 1+,sin Sinusoidal components of the fundamental voltage positive sequence;
voltage U of fundamental positive sequence component 1+ The method comprises the following steps:
active power P of fundamental positive sequence component 1+ The method comprises the following steps:
reactive power Q of fundamental positive sequence component 1+ The method comprises the following steps:
active current I of fundamental positive sequence component P1+ The method comprises the following steps:
reactive current I of fundamental positive sequence component Q1+ The method comprises the following steps:
7. the new energy unit electromechanical model checking method based on the whole group of fault ride through tests as claimed in claim 5, wherein the method is characterized by comprising the following steps: if the deviation allowance requirement cannot be completely met all the time after the model parameters are adjusted for a plurality of times, a group of model parameters with the minimum weighted average absolute deviation is selected as a final result.
8. The new energy unit electromechanical model checking method based on the whole group of fault ride through tests as claimed in claim 7, wherein the method is characterized by comprising the following steps: the times are not less than ten times.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311629258.3A CN117665445A (en) | 2023-11-30 | 2023-11-30 | New energy unit electromechanical model checking method based on whole group fault ride-through test |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311629258.3A CN117665445A (en) | 2023-11-30 | 2023-11-30 | New energy unit electromechanical model checking method based on whole group fault ride-through test |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117665445A true CN117665445A (en) | 2024-03-08 |
Family
ID=90065508
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311629258.3A Pending CN117665445A (en) | 2023-11-30 | 2023-11-30 | New energy unit electromechanical model checking method based on whole group fault ride-through test |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117665445A (en) |
-
2023
- 2023-11-30 CN CN202311629258.3A patent/CN117665445A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103454521B (en) | A kind of wind power plant power grid operation simulator | |
| CN104215904B (en) | A kind of full power convertor low-voltage ride-through testing system of wind generating set and method | |
| Li et al. | Transient voltage control of sending-end wind farm using a synchronous condenser under commutation failure of HVDC transmission system | |
| Liu et al. | A dynamic equivalent model for DFIG-based wind farms | |
| CN107766596B (en) | Low-voltage ride-through capability simulation evaluation method based on typical fault working condition set | |
| Liu et al. | Research on short-circuit currents calculation method considering dynamic reactive power support of renewable energy systems | |
| CN110311415B (en) | Transient stability margin evaluation method for renewable energy power generation base | |
| Qi et al. | Generic EMT modeling method of Type-4 wind turbine generators based on detailed FRT studies | |
| CN104678302A (en) | Low voltage ride-through test system and method for wind turbine generator of Boost current transformer | |
| CN103501010B (en) | The wind energy turbine set reactive power support method of a kind of pair of Hysteresis control | |
| Li et al. | DFIG-based wind farm electromagnetic dynamic model and impact on protection relay of transmission network | |
| CN203479929U (en) | Grid operation simulating device for wind power plant | |
| Wang et al. | Development of wind-energy modeling technology and standards | |
| CN117665445A (en) | New energy unit electromechanical model checking method based on whole group fault ride-through test | |
| CN116187082A (en) | Single-machine equivalent modeling method for wind power plant | |
| Liu et al. | A new calculation method of short-circuit currents contributed by doubly-fed wind turbines cluster | |
| Li et al. | Structure preserving aggregation method for doubly-fed induction generators in wind power conversion | |
| CN108448595B (en) | Photovoltaic-series compensation system time domain simulation small disturbance analysis method and system | |
| Li et al. | Wind farm electromagnetic dynamic model and outgoing line protection relay RTDS testing | |
| CN106972510B (en) | Directly driven wind-powered sub-synchronous oscillation analysis method based on the delay of net side control loop | |
| Gevorgian et al. | Wgrid-49 GMLC project report: Understanding the role of short-term energy storage and large motor loads for active power controls by wind power | |
| Dong et al. | Group equivalent modeling and analysis of group subsynchronous oscillation of direct-drive wind turbine | |
| Zeng et al. | Research on Wind Turbine Control Hardware in Loop Real-time Simulation and its Application | |
| Joshi et al. | Performance Enhancement of Doubly Fed Induction Generator Based Wind Farm for Grid Voltage Dip and Harmonics Mitigation | |
| Ma et al. | Research on Transient Control Characteristic Testing Technology for New Energy Connected to Weak Grid |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |