CN114362210B - Wind farm oscillation risk assessment test method, avoidance method and storage medium - Google Patents
Wind farm oscillation risk assessment test method, avoidance method and storage medium Download PDFInfo
- Publication number
- CN114362210B CN114362210B CN202210029824.6A CN202210029824A CN114362210B CN 114362210 B CN114362210 B CN 114362210B CN 202210029824 A CN202210029824 A CN 202210029824A CN 114362210 B CN114362210 B CN 114362210B
- Authority
- CN
- China
- Prior art keywords
- oscillation
- frequency
- risk
- mode
- signal
- 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
- 230000010355 oscillation Effects 0.000 title claims abstract description 185
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000012502 risk assessment Methods 0.000 title claims abstract description 40
- 238000010998 test method Methods 0.000 title claims abstract description 37
- 239000011159 matrix material Substances 0.000 claims description 52
- 238000000354 decomposition reaction Methods 0.000 claims description 28
- 238000013016 damping Methods 0.000 claims description 27
- 238000002347 injection Methods 0.000 claims description 14
- 239000007924 injection Substances 0.000 claims description 14
- 238000012216 screening Methods 0.000 claims description 11
- 230000009466 transformation Effects 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000009467 reduction Effects 0.000 claims description 5
- 238000011217 control strategy Methods 0.000 abstract description 5
- 230000015654 memory Effects 0.000 description 20
- 238000012360 testing method Methods 0.000 description 10
- 238000012545 processing Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 206010014357 Electric shock Diseases 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
Classifications
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Wind Motors (AREA)
Abstract
The invention discloses a wind farm oscillation risk assessment test method, an avoidance method and a storage medium, wherein the assessment test method comprises the following steps: collecting electrical signals at grid connection points of a wind farm; calculating an oscillation risk index of each mode according to the signal mode corresponding to each electric signal and the frequency of the mode; injecting disturbance current signals into a preset number of grid-connected points selected based on the oscillation risk indexes according to the frequency of the mode; and determining whether the wind power plant has oscillation risk according to the electric signals at the grid connection point after the disturbance current signals are injected. According to the wind power plant oscillation risk assessment test method provided by the embodiment of the invention, the risk weak grid-connected point is identified based on the oscillation risk index, and then the judgment of the oscillation risk of the wind power plant is realized by directly injecting the disturbance current signal at the physical side, so that the oscillation risk existing in the wind power plant can be effectively identified, and effective information is provided for formulating a stable control strategy of the wind power plant.
Description
Technical Field
The invention relates to the technical field of broadband oscillation of power systems, in particular to a wind power plant oscillation risk assessment test method, an avoidance method and a storage medium.
Background
According to the data published by the national energy agency, the newly-increased installed capacity of the power supply is 19087 kilowatts by the year 2020, wherein the installed capacity of the wind power grid-connected system reaches 7167 kilowatts, the occupied ratio reaches 37.5%, and the wind power integration installation breaks through 2.8 hundred million kilowatts. However, the problem of power oscillation caused by interaction of a large-scale wind power grid-connected system and an alternating current power grid poses a great threat to stable operation of the system, and a plurality of related oscillation accidents are generated at present: such as system power oscillation caused by the action of the doubly-fed wind turbine and the series compensation capacitor, system power oscillation caused by the action of the direct-driven wind turbine and the alternating-current power grid, and the like.
When the wind farm has the risk of broadband oscillation, the probability of attack is also greatly improved. However, damage to the wind farm may cause abnormal fluctuation of regional power supply, and cause destructive influence on the whole power grid, so that an oscillation risk assessment test for the wind farm is necessary. The oscillation risk frequency existing in the system can be effectively identified, and effective information is provided for formulating a stable control strategy of the wind power plant.
Disclosure of Invention
In view of the above, the embodiment of the invention provides a wind farm oscillation risk assessment test method, an avoidance method and a storage medium, so as to solve the technical problem that wind farm oscillation risk assessment is lack in the prior art.
The technical scheme provided by the invention is as follows:
the first aspect of the embodiment of the invention provides a wind farm oscillation risk assessment test method, which comprises the following steps: collecting electrical signals at grid connection points of a wind farm; calculating an oscillation risk index of each mode according to the signal mode corresponding to each electric signal and the frequency of the mode; injecting disturbance current signals into a preset number of grid-connected points selected based on the oscillation risk indexes according to the modal frequency; and determining whether the wind power plant has oscillation risk according to the electric signals at the grid connection point after the disturbance current signals are injected.
Optionally, calculating the oscillation risk index of each mode according to the signal mode corresponding to each electric signal and the frequency of the mode includes: performing modal analysis according to the electrical signals to obtain corresponding signal modes; calculating the frequency and damping ratio of each mode according to the signal modes; screening the damping ratio and the frequency according to a damping ratio threshold value and a preset broadband oscillation key frequency band range to obtain a corresponding mode; and calculating the vibration risk index of the mode obtained by screening according to the frequency and the damping ratio.
Optionally, injecting disturbance current signals into a preset number of grid-connected points selected based on the oscillation risk index according to the frequency of the mode, including: ordering the corresponding modes according to the oscillation risk indexes; selecting a preset number of modes according to the sorting result; and injecting disturbance current signals into grid-connected points corresponding to the preset number of modes according to the frequency of the modes.
Optionally, determining whether the wind farm has an oscillation risk according to the electrical signal at the grid-connected point after the disturbance current signal is injected includes: collecting the electric signals of each fan at the grid-connected point after the disturbance current signals are injected; judging whether the corresponding fans are off-grid according to the relation between the electrical signals and a preset threshold value; and determining whether the wind power plant has oscillation risk according to the number of fans which are off-grid.
A second aspect of the embodiment of the present invention provides a wind farm oscillation risk avoidance method, including: when the wind power plant oscillation risk assessment test method according to the first aspect of the embodiment of the invention determines that the wind power plant has an oscillation risk, acquiring participation indexes of each fan of the wind power plant during oscillation; determining an oscillating dominant fan according to participation indexes of each fan; modifying parameters of the control link corresponding to the dominant wind turbine, and repeating the wind power plant oscillation risk assessment test method according to the first aspect of the embodiment of the invention until the wind power plant has no oscillation risk.
Optionally, determining the oscillating dominant fan according to the participation index of each fan includes: singular value decomposition is carried out according to the electrical signal data matrix of the fan, so that a transformation matrix is obtained; performing singular value decomposition according to a low-dimensional matrix obtained by performing dimension reduction on the electrical signal data matrix based on the transformation matrix to obtain a decomposition matrix; performing eigenvalue decomposition according to the approximate state matrix obtained by calculation of the decomposition matrix and the low-dimensional matrix to obtain eigenvalues and eigenvectors; calculating according to the characteristic value and the characteristic vector to obtain a participation factor; and determining the position of the dominant fan according to the participation factors.
Optionally, modifying parameters of the control link corresponding to the dominant fan includes: calculating an oscillation frequency according to the characteristic value; and adjusting parameters of the corresponding control links of the leading fans according to the range of the oscillation frequency.
Optionally, according to the range of the oscillation frequency, adjusting parameters of a control link corresponding to the dominant fan includes: when the oscillation frequency is in the range of 1Hz-50Hz, PI parameters of a direct-current voltage ring of a leading fan and a phase-locked loop controller are adjusted, and the bandwidth of the direct-current voltage ring is increased or the bandwidth of the phase-locked loop is reduced; and when the oscillation frequency is within the range of 50Hz-2500Hz, PI parameters of the phase-locked loop of the leading fan and the current inner loop controller are adjusted, and the bandwidth of the phase-locked loop is reduced or the bandwidth of the current inner loop is reduced.
A third aspect of the embodiment of the present invention provides a computer-readable storage medium, where computer instructions are stored, where the computer instructions are configured to cause the computer to perform the wind farm oscillation risk assessment test method according to the first aspect of the embodiment of the present invention and any one of the first aspect of the embodiment of the present invention and the avoidance method according to the second aspect of the embodiment of the present invention and any one of the second aspect of the embodiment of the present invention.
A fourth aspect of an embodiment of the present invention provides an electronic device, including: the wind farm oscillation risk assessment test method according to the first aspect of the embodiment of the present invention and the avoidance method according to the second aspect of the embodiment of the present invention are performed by the processor through executing the computer instructions.
The technical scheme provided by the invention has the following effects:
according to the wind power plant oscillation risk assessment test method provided by the embodiment of the invention, the oscillation risk index is calculated by collecting the electrical signals at the grid-connected points and based on the frequency determined by the signal mode of the electrical signals, the disturbance current signals are injected into the grid-connected points selected according to the oscillation risk index, and whether the wind power plant has oscillation risk is judged based on the electrical signals after the disturbance current signals are injected. Therefore, the assessment test method identifies the weak risk grid-connected point based on the oscillation risk index, and then judges the oscillation risk of the wind power plant by directly injecting the disturbance current signal at the physical side, so that the oscillation risk existing in the wind power plant can be effectively identified, and effective information is provided for formulating a stable control strategy of the wind power plant.
According to the wind power plant oscillation risk avoiding method provided by the embodiment of the invention, when the wind power plant has oscillation risk, the participation index in the wind power plant is obtained, the vibrating dominant fan is determined according to the participation index, and the control link parameters of the dominant fan are adjusted until the wind power plant has no oscillation risk. Therefore, the risk avoiding method realizes risk avoiding of the wind power plant with oscillation risk, and avoids damage attack possibly suffered by the wind power plant.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are needed in the description of the embodiments or the prior art will be briefly described, and it is obvious that the drawings in the description below are some embodiments of the present invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a wind farm oscillation risk assessment test method according to an embodiment of the invention;
FIG. 2 is a flow chart of a wind farm oscillation risk assessment test method according to another embodiment of the invention;
FIG. 3 is a flow chart of a wind farm oscillation risk avoidance method according to an embodiment of the invention;
FIG. 4 is a flow chart of a wind farm oscillation risk avoidance method according to another embodiment of the invention;
FIG. 5 is a flow chart of a wind farm oscillation risk avoidance method according to another embodiment of the invention;
FIG. 6 is a block diagram of a wind farm oscillation risk assessment test apparatus according to an embodiment of the present invention;
FIG. 7 is a block diagram of a wind farm oscillation risk avoidance apparatus according to an embodiment of the invention;
FIG. 8 is a schematic diagram of a computer-readable storage medium provided according to an embodiment of the present invention;
fig. 9 is a schematic structural diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
In order that those skilled in the art will better understand the present invention, a technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in which it is apparent that the described embodiments are only some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the present invention without making any inventive effort, shall fall within the scope of the present invention.
The terms first, second, third, fourth and the like in the description and in the claims and in the above drawings are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used may be interchanged where appropriate such that the embodiments described herein may be implemented in other sequences than those illustrated or otherwise described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
According to an embodiment of the present invention, a wind farm oscillation risk assessment test method and an avoidance method are provided, it should be noted that the steps illustrated in the flowchart of the drawings may be performed in a computer system such as a set of computer executable instructions, and that although a logical order is illustrated in the flowchart, in some cases the steps illustrated or described may be performed in an order different from that illustrated herein.
In this embodiment, a wind farm oscillation risk assessment test method is provided, which may be used in electronic devices, such as a computer, a mobile phone, a tablet computer, etc., fig. 1 is a flowchart of a wind farm oscillation risk assessment test method according to an embodiment of the present invention, and as shown in fig. 1, the method includes the following steps:
the embodiment of the invention provides a wind farm oscillation risk assessment test method, as shown in fig. 1, comprising the following steps:
step S101: and collecting electrical signals at the grid connection point of the wind power plant. The grid connection point can be a connection point of the wind power plant and the power grid. Specifically, a broadband signal measurement device may be used to collect electrical signals at the point of connection. In one embodiment, the broadband signal measurement device has a measurement range of 1Hz-5000 Hz. The electrical signal collected may be a power signal or a current signal. When the collected electric signal is a power signal, high-pass filtering processing is carried out on the collected signal, and direct current components in the signal are filtered. When the collected electric signal is a current signal, notch filtering processing is carried out on the collected signal, and the power frequency component in the signal is filtered. Meanwhile, the collected electrical signals can be processed by abnormal values. The outlier processing comprises the steps of removing outliers in signals and filling the removed outliers by adopting average values of adjacent data.
Step S102: and calculating an oscillation risk index of each mode according to the signal mode corresponding to each electric signal and the frequency of the mode. Specifically, after the electrical signals are collected and preprocessed, the electrical signals can be subjected to modal analysis to obtain corresponding signal modes. And then calculating the frequency of each mode according to the calculated signal modes, and calculating a corresponding oscillation risk index based on the frequency.
Step S103: and injecting disturbance current signals into the preset number of grid-connected points selected based on the oscillation risk indexes according to the frequency of the mode. Specifically, after the oscillation risk index of the mode is determined, a preset number of grid-connected points can be selected based on a preset threshold. And then injecting current disturbance signals into the corresponding grid-connected points based on the corresponding frequencies. Wherein, a disturbance injection device can be adopted to inject disturbance current signals into the grid-connected point. The injection of the disturbance current signal may last for a preset time, after which the injection is stopped.
Step S104: and determining whether the wind power plant has oscillation risk according to the electric signals at the grid connection point after the disturbance current signals are injected. Specifically, after the disturbance current signal is injected, an electrical signal at a grid-connected point of the wind power plant is collected, for example, a voltage and a current signal at the grid-connected point are collected, and whether the voltage and the current are 0 is judged, so that whether the wind power plant has an oscillation risk is determined. The determination of whether the voltage and the current are 0 is one way to determine whether the oscillation disconnection occurs, or other ways may be adopted to determine whether the oscillation disconnection occurs.
According to the wind power plant oscillation risk assessment test method provided by the embodiment of the invention, the oscillation risk index is calculated by collecting the electrical signals at the grid-connected points and based on the frequency determined by the signal mode of the electrical signals, the disturbance current signals are injected into the grid-connected points selected according to the oscillation risk index, and whether the wind power plant has oscillation risk is judged based on the electrical signals after the disturbance current signals are injected. Therefore, the assessment test method identifies the weak risk grid-connected point based on the oscillation risk index, and then judges the oscillation risk of the wind power plant by directly injecting the disturbance current signal at the physical side, so that the oscillation risk existing in the wind power plant can be effectively identified, and effective information is provided for formulating a stable control strategy of the wind power plant.
In an embodiment, as shown in fig. 2, the oscillation risk indicator of each mode is calculated according to the signal mode corresponding to each electric signal and the frequency of the mode, and the method includes the following steps:
step S201: and carrying out modal analysis according to the electrical signals to obtain corresponding signal modes. In particular, the modal analysis may employ a global least squares method-rotation invariant algorithm (Estimating Signal Parameter via Rotational Invariance Techniques, TLS-ESPRIT), the resulting modal signal being denoted λ i =α i +jω i (i=1, 2, …, I). Wherein lambda is i Representing the computed feature root, alpha i Representing the real part, omega i Representing the imaginary part. i represents the number of corresponding electrical signals.
Step S202: and calculating the frequency and damping ratio of each mode according to the signal modes. Specifically, after the corresponding signal mode is calculated, the frequency and damping ratio of each mode of the signal are calculated based on the signal mode.
Wherein the frequency and damping ratio are expressed by the following formula:
wherein f i Represent frequency, ζ i Representing the damping ratio.
Step S203: and screening the damping ratio and the frequency according to a damping ratio threshold value and a preset broadband oscillation key frequency band range to obtain a corresponding mode. Wherein the preset broadband oscillation critical frequency band range is expressed as [ f ] min , f max ]In an embodiment, the preset wideband oscillation critical frequency band range is 2Hz-2500Hz, and in other embodiments, the range may be set to other ranges according to practical situations, and the embodiment of the present invention does not limit the specific range. Furthermore, the damping ratio threshold may be determined based on actual conditions.
And screening the damping ratio of the modes obtained by calculation according to the damping ratio threshold, selecting the mode with the damping ratio smaller than the damping ratio threshold, judging whether the frequency of the selected mode is in a preset broadband oscillation key frequency band range, and screening out the corresponding mode when the frequency of the selected mode is in the preset broadband oscillation key frequency band range. In an embodiment, during the screening of the modes, whether the mode is located in the preset broadband oscillation critical frequency band range or not may be determined first, and then the mode is compared with the damping ratio threshold value.
Step S204: and calculating the vibration risk index of the mode obtained by screening according to the frequency and the damping ratio. For the screened modes, the oscillation risk index is calculated according to the following formula:
wherein k is 1 ,k 2 Is a constant coefficient.
In an embodiment, injecting disturbance current signals into a preset number of grid-connected points selected based on the oscillation risk index according to the frequency of the mode includes the following steps: ordering the corresponding modes according to the oscillation risk indexes; selecting a preset number of modes according to the sorting result; and injecting disturbance current signals into grid-connected points corresponding to the preset number of modes according to the frequency of the modes.
Specifically, after the oscillation risk indexes of the screened modes are calculated, the corresponding oscillation risk indexes are ranked from large to small according to the sizes of the oscillation risk indexes, and then corresponding modes with the corresponding number are selected as test modes according to a preset oscillation risk index threshold, for example, the number of the selected test modes is S. After the test mode is selected, a disturbance current signal with preset time is injected into the grid-connected point corresponding to the selected test mode by adopting a preset disturbance injection device.
Before the disturbance current signal is injected, judging whether the signal in the modal analysis is a power signal or a current signal; if the power signal is the power signal, firstly acquiring the rated frequency of the system, and calculating the injection frequency of the disturbance current signal based on the rated frequency and the frequency of the corresponding mode, namely, the injection frequency is expressed as:
f ds =f 0 ±f s
wherein f 0 For the nominal frequency of the system, f s Representing the frequency of the corresponding modality.
If it isThe injection frequency of the disturbance current signal is determined directly based on the frequency of the corresponding mode, i.e. the injection frequency of the disturbance current signal is f ds =f s . In addition, the amplitude of the disturbance current signal is 5% -10% of the rated current when the disturbance current signal is injected.
In an embodiment, determining whether the wind farm has a risk of oscillation according to the electrical signal at the grid-connected point after the disturbance current signal is injected includes: collecting the electric signals of each fan at the grid-connected point after the disturbance current signals are injected; judging whether the corresponding fans are off-grid according to the relation between the electrical signals and a preset threshold value; and determining whether the wind power plant has oscillation risk according to the number of fans which are off-grid.
Specifically, after the disturbance current signal of the preset time is injected, the outlet voltage current signal of each fan at the grid-connected point of the wind power plant is collected, whether the voltage current signal is 0 is judged, and if the voltage current signal is 0, the corresponding fan is considered to be off-grid. After judging each fan, counting the total number of all fans which are off-grid, and judging whether the total number of the fans exceeds a limit value N lim If the wind power plant is out of network, the wind power plant is considered to have oscillation risk, wherein N lim The number of the fans of the wind power plant can be 10-20%.
The embodiment of the invention also provides a wind farm oscillation risk avoiding method, as shown in fig. 3, which comprises the following steps:
step S301: when the wind farm oscillation risk assessment test method according to any one of the embodiments determines that the wind farm has an oscillation risk, acquiring participation indexes of each fan of the wind farm during oscillation. Specifically, after the wind power plant is judged to have the oscillation risk, current signals of all fans of the wind power plant during oscillation are collected.
Step S302: and determining the oscillating dominant fan according to the participation index of each fan. Specifically, after the current parameters of each fan are collected, the oscillating dominant fan can be determined through singular value decomposition and other processes according to the data matrix constructed by the current parameters.
Step S303: modifying parameters of the control link corresponding to the dominant wind turbine, and repeating the wind power plant oscillation risk assessment test method according to any one of the embodiments until the wind power plant has no oscillation risk. After the main wind guiding machine is confirmed, parameters of control links such as a direct-current voltage ring, a phase-locked loop and a current inner ring of the main wind guiding machine are adjusted, after adjustment, whether the vibration risk exists is judged according to the wind power plant vibration risk assessment test method, and if the vibration risk exists, the adjustment and judgment process is continued until the wind power plant does not exist.
According to the wind power plant oscillation risk avoiding method provided by the embodiment of the invention, when the wind power plant has oscillation risk, the participation index in the wind power plant is obtained, the vibrating dominant fan is determined according to the participation index, and the control link parameters of the dominant fan are adjusted until the wind power plant has no oscillation risk. Therefore, the risk avoiding method realizes risk avoiding of the wind power plant with oscillation risk, and avoids damage attack possibly suffered by the wind power plant.
In one embodiment, as shown in fig. 4, determining the oscillating dominant blower according to the participation index of each blower includes the following steps:
step S401: and carrying out singular value decomposition according to the electrical signal data matrix of the fan to obtain a transformation matrix. Specifically, the current signals at a plurality of moments during oscillation can be acquired in real time to construct a data matrix. Wherein t can be selected 1 To t n Constructing a data matrix X at the moment, and selecting t 2 To t n+1 Time of day construction data matrix X 1 . The two data matrices are represented as:
wherein matrices X and X 1 The order of (m×n), W 1 To W m Representing a plurality of fans.
For the two constructed data matrixes, singular value decomposition is carried out on the data matrix X, and then a right singular vector matrix V obtained by singular value decomposition is stacked or intercepted according to a target dimension reduction multiple p to obtain a transformation matrix C. The singular value decomposition process adopts singular value decomposition in the prior art, and is not described herein.
Step S402: and performing singular value decomposition according to a low-dimensional matrix obtained by performing dimension reduction on the electric signal data matrix based on the transformation matrix to obtain a decomposition matrix. Specifically, after the transformation matrix is obtained, the data matrix X and the data matrix X are respectively mapped based on the transformation matrix 1 Performing dimension reduction processing to obtain a low-dimension matrix Y and a low-dimension matrix Y corresponding to the two data matrixes 1 . Wherein, the two low-dimensional matrices are calculated by adopting the following two formulas:
wherein: y and Y 1 The order of the matrix is (m×a), a=n/p;
after obtaining two low-dimensional matrixes, performing singular value decomposition on the low-dimensional matrix Y to obtain a decomposition matrix U y 、 S y And V y The three decomposition matrices satisfy the equation y=u y ×S y ×V y ’,V y ' is V y Is a transpose of (a).
Step S403: and carrying out eigenvalue decomposition according to the approximate state matrix obtained by calculation of the decomposition matrix and the low-dimensional matrix to obtain eigenvalues and eigenvectors. Specifically, for the three decomposition matrices obtained, they are summed with a low-dimensional matrix Y 1 An approximation state matrix of a computing system, the approximation state matrix being represented as a=u y ’× Y 1 ×V y ×S y The method comprises the steps of carrying out a first treatment on the surface of the Then, the approximate state matrix A is subjected to eigenvalue decomposition to obtain eigenvalue lambda y And feature vector W y 。
Step S404: calculating according to the characteristic value and the characteristic vector to obtain a participation factor; specifically, the present invention relates to a method for manufacturing a semiconductor device. The oscillation mode Φ=u of the actual system can be calculated using the calculated eigenvalues and eigenvectors y ×W y Eigenvalue λ=λ y The method comprises the steps of carrying out a first treatment on the surface of the The oscillation mode is then used to calculate the participation factor P f . The participation factor is calculated by the following formula:
the matrix ψ is an inverse matrix of the system oscillation mode Φ, and k and i respectively represent any two of the singular values reserved after singular value decomposition is performed on the data matrix X.
Step S405: and determining the position of the dominant fan according to the participation factors. After the participation factors corresponding to any two k and i are calculated based on the formula, determining the position of the dominant fan according to the sizes of the participation factors. Specifically, after a plurality of participation factors are calculated, the participation factors with orders of magnitude different from the other participation factors are determined, and then a dominant fan is determined according to fans corresponding to the participation factors.
In an embodiment, as shown in fig. 5, the modification of the parameters of the control link corresponding to the dominant fan includes the following steps:
step S501: and calculating the oscillation frequency according to the characteristic value. Specifically, the eigenvector W is obtained by eigenvalue decomposition of the approximated state matrix y And a characteristic value lambda y Then, the oscillation frequency and the damping ratio are calculated using the eigenvalues. The oscillation frequency and damping ratio are calculated using the following formula:
wherein: r=1, 2, …, R being the number of singular values remaining after singular value decomposition of the data matrix X; Δt is the oscillation data time interval.
Step S502: and adjusting parameters of the corresponding control links of the leading fans according to the range of the oscillation frequency. Specifically, after the oscillation frequency is calculated according to the formula, the parameters of the control link corresponding to the dominant fan are adjusted according to the range of the oscillation frequency. In one embodiment, when the oscillation frequency is within the range of 1Hz-50Hz, PI parameters of the direct-current voltage ring of the leading fan and the phase-locked loop controller are adjusted, and the bandwidth of the direct-current voltage ring is increased or the bandwidth of the phase-locked loop is reduced; and when the oscillation frequency is within the range of 50Hz-2500Hz, PI parameters of the phase-locked loop of the leading fan and the current inner loop controller are adjusted, and the bandwidth of the phase-locked loop is reduced or the bandwidth of the current inner loop is reduced.
The embodiment of the invention also provides a wind farm oscillation risk assessment test device, as shown in fig. 6, which comprises:
the signal acquisition module is used for acquiring electrical signals at the grid-connected point of the wind power plant; the specific content refers to the corresponding parts of the above method embodiments, and will not be described herein.
The index calculation module is used for calculating an oscillation risk index of each mode according to the signal mode corresponding to each electric signal and the frequency of the mode; the specific content refers to the corresponding parts of the above method embodiments, and will not be described herein.
The disturbance injection module is used for injecting disturbance current signals to the grid-connected points with the preset number selected based on the oscillation risk indexes according to the modal frequency; the specific content refers to the corresponding parts of the above method embodiments, and will not be described herein.
And the risk judging module is used for determining whether the wind power plant has oscillation risk or not according to the electric signals at the grid-connected point after the disturbance current signals are injected. The specific content refers to the corresponding parts of the above method embodiments, and will not be described herein.
According to the wind power plant oscillation risk assessment test device provided by the embodiment of the invention, the oscillation risk index is calculated by collecting the electrical signals at the grid-connected points and based on the frequency determined by the signal mode of the electrical signals, the disturbance current signals are injected into the grid-connected points selected according to the oscillation risk index, and whether the wind power plant has oscillation risk is judged based on the electrical signals after the disturbance current signals are injected. Therefore, the assessment test device identifies the weak risk grid-connected point based on the oscillation risk index, and then judges the oscillation risk of the wind power plant by directly injecting the disturbance current signal at the physical side, so that the oscillation risk existing in the wind power plant can be effectively identified, and effective information is provided for formulating a stable control strategy of the wind power plant.
The functional description of the wind farm oscillation risk assessment test device provided by the embodiment of the invention is described in detail by referring to the wind farm oscillation risk assessment test method in the embodiment.
The embodiment of the invention also provides a wind farm oscillation risk avoiding device, as shown in fig. 7, which comprises:
the oscillation parameter obtaining module is used for obtaining participation indexes of each fan of the wind power plant during oscillation when the wind power plant oscillation risk assessment test method according to any one of the embodiments determines that the wind power plant has an oscillation risk; the specific content refers to the corresponding parts of the above method embodiments, and will not be described herein.
The dominant fan determining module is used for determining an oscillating dominant fan according to the participation indexes of each fan; the specific content refers to the corresponding parts of the above method embodiments, and will not be described herein.
And the adjusting module is used for modifying the parameters of the control links corresponding to the dominant fans, and repeating the wind power plant oscillation risk assessment test method according to any one of the embodiments until the wind power plant has no oscillation risk. The specific content refers to the corresponding parts of the above method embodiments, and will not be described herein.
According to the wind power plant oscillation risk avoiding device provided by the embodiment of the invention, when the wind power plant has oscillation risk, the participation index in the wind power plant is obtained, the vibrating dominant fan is determined according to the participation index, and the control link parameters of the dominant fan are adjusted until the wind power plant has no oscillation risk. Therefore, the risk avoiding device realizes risk avoiding of the wind power plant with oscillation risk, and damage attack possibly suffered by the wind power plant is avoided.
The functional description of the wind farm oscillation risk avoiding device provided by the embodiment of the invention is detailed in the description of the wind farm oscillation risk avoiding method in the embodiment.
The embodiment of the present invention further provides a storage medium, as shown in fig. 8, on which a computer program 601 is stored, which when executed by a processor, implements the steps of the electric-shock-field oscillation risk assessment test method and the avoidance method of the above embodiment. The storage medium also stores audio and video stream data, characteristic frame data, interactive request signaling, encrypted data, preset data size and the like. The storage medium may be a magnetic Disk, an optical Disk, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), a Flash Memory (Flash Memory), a Hard Disk (HDD), or a Solid State Drive (SSD); the storage medium may also comprise a combination of memories of the kind described above.
It will be appreciated by those skilled in the art that implementing all or part of the above-described embodiment method may be implemented by a computer program to instruct related hardware, where the program may be stored in a computer readable storage medium, and the program may include the above-described embodiment method when executed. Wherein the storage medium may be a magnetic Disk, an optical Disk, a Read-Only Memory (ROM), a random access Memory (RandomAccessMemory, RAM), a Flash Memory (Flash Memory), a Hard Disk (HDD), a Solid State Drive (SSD), or the like; the storage medium may also comprise a combination of memories of the kind described above.
The embodiment of the present invention further provides an electronic device, as shown in fig. 9, where the electronic device may include a processor 51 and a memory 52, where the processor 51 and the memory 52 may be connected by a bus or other means, and in fig. 9, the connection is exemplified by a bus.
The processor 51 may be a central processing unit (Central Processing Unit, CPU). The processor 51 may also be other general purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), field programmable gate arrays (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or combinations thereof.
The memory 52 serves as a non-transitory computer readable storage medium that may be used to store non-transitory software programs, non-transitory computer-executable programs, and modules, such as corresponding program instructions/modules in embodiments of the present invention. The processor 51 executes various functional applications and data processing of the processor by running non-transitory software programs, instructions and modules stored in the memory 52, i.e. implementing the wind farm oscillation risk assessment test method and avoidance method in the above-described method embodiments.
The memory 52 may include a memory program area that may store an operating device, an application program required for at least one function, and a memory data area; the storage data area may store data created by the processor 51, etc. In addition, memory 52 may include high-speed random access memory, and may also include non-transitory memory, such as at least one magnetic disk storage device, flash memory device, or other non-transitory solid state storage device. In some embodiments, memory 52 may optionally include memory located remotely from processor 51, which may be connected to processor 51 via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
The one or more modules are stored in the memory 52, which when executed by the processor 51, perform the wind farm oscillation risk assessment test method and avoidance method in the embodiments shown in fig. 1-5.
The specific details of the electronic device may be understood correspondingly with reference to the corresponding related descriptions and effects in the embodiments shown in fig. 1 to 5, which are not repeated here.
Although embodiments of the present invention have been described in connection with the accompanying drawings, various modifications and variations may be made by those skilled in the art without departing from the spirit and scope of the invention, and such modifications and variations fall within the scope of the invention as defined by the appended claims.
Claims (8)
1. A wind farm oscillation risk assessment test method, comprising:
collecting electrical signals at grid connection points of a wind farm;
calculating an oscillation risk index of each mode according to the signal mode corresponding to each electric signal and the frequency of the mode;
injecting disturbance current signals into a preset number of grid-connected points selected based on the oscillation risk indexes according to the modal frequency;
determining whether the wind power plant has oscillation risk according to the electric signals at the grid connection point after the disturbance current signals are injected;
calculating an oscillation risk index of each mode according to the signal mode corresponding to each electric signal and the frequency of the mode, wherein the method comprises the following steps:
performing modal analysis according to the electrical signals to obtain corresponding signal modes;
calculating the frequency and damping ratio of each mode according to the signal modes;
screening the damping ratio and the frequency according to a damping ratio threshold value and a preset broadband oscillation key frequency band range to obtain a corresponding mode;
calculating the vibration risk index of the mode obtained by screening according to the frequency and the damping ratio;
injecting disturbance current signals into a preset number of grid-connected points selected based on the oscillation risk indexes according to the frequency of the mode, wherein the disturbance current signals comprise:
ordering the corresponding modes according to the oscillation risk indexes;
selecting a preset number of modes according to the sorting result;
injecting disturbance current signals into grid-connected points corresponding to the preset number of modes according to the frequency of the modes;
screening the damping ratio and the frequency according to a damping ratio threshold value and a preset broadband oscillation key frequency band range to obtain corresponding modes, wherein the method comprises the following steps:
selecting a mode with a damping ratio smaller than the damping ratio threshold, judging whether the frequency of the selected mode is in a preset broadband oscillation key frequency band range, and screening out the corresponding mode when the frequency of the selected mode is in the preset broadband oscillation key frequency band range;
injecting disturbance current signals into grid-connected points corresponding to the preset number of modes according to the frequency of the modes, wherein the disturbance current signals comprise:
before the disturbance current signal is injected, judging whether the signal is a power signal or a current signal when the modal analysis is carried out; if the power signal is the power signal, firstly acquiring the rated frequency of the system, and calculating the injection frequency of the disturbance current signal based on the rated frequency and the frequency of the corresponding mode, namely, the injection frequency is expressed as:
wherein,for the nominal frequency of the system, +.>Representing the frequency of the corresponding modality;
if the current signal is the current signal, determining the injection frequency of the disturbance current signal directly based on the frequency of the corresponding mode, namely the injection frequency of the disturbance current signal is。
2. The wind farm oscillation risk assessment test method according to claim 1, wherein determining whether the wind farm is at risk of oscillation based on the electrical signal at the grid-tie point after injection of the disturbance current signal, comprises:
collecting the electric signals of each fan at the grid-connected point after the disturbance current signals are injected;
judging whether the corresponding fans are off-grid according to the relation between the electrical signals and a preset threshold value;
and determining whether the wind power plant has oscillation risk according to the number of fans which are off-grid.
3. A wind farm oscillation risk avoidance method, comprising:
when the wind farm oscillation risk assessment test method according to claim 1 or 2 determines that the wind farm has an oscillation risk, acquiring participation indexes of each fan of the wind farm during oscillation;
determining an oscillating dominant fan according to participation indexes of each fan;
modifying parameters of the corresponding control link of the dominant wind turbine, and repeating the wind farm oscillation risk assessment test method according to claim 1 or 2 until the wind farm has no oscillation risk.
4. A wind farm oscillation risk avoidance method according to claim 3, wherein determining the dominant wind turbine for oscillation from the participation indicators of each wind turbine comprises:
singular value decomposition is carried out according to the electrical signal data matrix of the fan, so that a transformation matrix is obtained;
performing singular value decomposition according to a low-dimensional matrix obtained by performing dimension reduction on the electrical signal data matrix based on the transformation matrix to obtain a decomposition matrix;
performing eigenvalue decomposition according to the approximate state matrix obtained by calculation of the decomposition matrix and the low-dimensional matrix to obtain eigenvalues and eigenvectors;
calculating according to the characteristic value and the characteristic vector to obtain a participation factor;
and determining the position of the dominant fan according to the participation factors.
5. The method for avoiding risk of wind farm oscillation according to claim 4, wherein modifying parameters of the control link corresponding to the dominant wind turbine comprises:
calculating an oscillation frequency according to the characteristic value;
and adjusting parameters of the corresponding control links of the leading fans according to the range of the oscillation frequency.
6. The method for avoiding risk of wind farm oscillation according to claim 5, wherein adjusting parameters of a control link corresponding to a dominant wind turbine according to the range of oscillation frequencies comprises:
when the oscillation frequency is in the range of 1Hz-50Hz, PI parameters of a direct-current voltage ring of a leading fan and a phase-locked loop controller are adjusted, and the bandwidth of the direct-current voltage ring is increased or the bandwidth of the phase-locked loop is reduced;
and when the oscillation frequency is within the range of 50Hz-2500Hz, PI parameters of the phase-locked loop of the leading fan and the current inner loop controller are adjusted, and the bandwidth of the phase-locked loop is reduced or the bandwidth of the current inner loop is reduced.
7. A computer-readable storage medium storing computer instructions for causing the computer to perform the wind farm oscillation risk assessment test method according to claim 1 or 2 and the wind farm oscillation risk avoidance method according to any of claims 3 to 6.
8. An electronic device, comprising: the wind farm oscillation risk assessment test method according to claim 1 or 2 and the wind farm oscillation risk avoidance method according to any one of claims 3 to 6 are executed by the processor.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210029824.6A CN114362210B (en) | 2022-01-12 | 2022-01-12 | Wind farm oscillation risk assessment test method, avoidance method and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210029824.6A CN114362210B (en) | 2022-01-12 | 2022-01-12 | Wind farm oscillation risk assessment test method, avoidance method and storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114362210A CN114362210A (en) | 2022-04-15 |
CN114362210B true CN114362210B (en) | 2023-11-21 |
Family
ID=81109200
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210029824.6A Active CN114362210B (en) | 2022-01-12 | 2022-01-12 | Wind farm oscillation risk assessment test method, avoidance method and storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114362210B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109103903A (en) * | 2018-09-13 | 2018-12-28 | 华北电力大学 | A kind of judgment method causing sub-synchronous oscillation for straight drive blower |
CN109586273A (en) * | 2017-09-29 | 2019-04-05 | 中国电力科学研究院 | Sub-synchronous oscillation appraisal procedure and system based on wind energy HVDC delivery system |
CN110137948A (en) * | 2019-05-14 | 2019-08-16 | 中国电力科学研究院有限公司 | A kind of appraisal procedure and device of double-fed fan motor play synchronized oscillation risk |
CN110854884A (en) * | 2019-12-07 | 2020-02-28 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | Wind power collection region subsynchronous oscillation risk online assessment and early warning method |
CN111525611A (en) * | 2020-04-26 | 2020-08-11 | 西安热工研究院有限公司 | Frequency coupling effect-considering doubly-fed grid-connected system subsynchronous oscillation analysis method |
-
2022
- 2022-01-12 CN CN202210029824.6A patent/CN114362210B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109586273A (en) * | 2017-09-29 | 2019-04-05 | 中国电力科学研究院 | Sub-synchronous oscillation appraisal procedure and system based on wind energy HVDC delivery system |
CN109103903A (en) * | 2018-09-13 | 2018-12-28 | 华北电力大学 | A kind of judgment method causing sub-synchronous oscillation for straight drive blower |
CN110137948A (en) * | 2019-05-14 | 2019-08-16 | 中国电力科学研究院有限公司 | A kind of appraisal procedure and device of double-fed fan motor play synchronized oscillation risk |
CN110854884A (en) * | 2019-12-07 | 2020-02-28 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | Wind power collection region subsynchronous oscillation risk online assessment and early warning method |
CN111525611A (en) * | 2020-04-26 | 2020-08-11 | 西安热工研究院有限公司 | Frequency coupling effect-considering doubly-fed grid-connected system subsynchronous oscillation analysis method |
Non-Patent Citations (3)
Title |
---|
Analysis and control on ultra low frequency oscillation at seeding end of UHVDC power system;Gong Tao-Rong et al;2014 international conference on power system technology;全文 * |
基于关联机组风险系数的发电机振荡模式辨识;林波;吕林;钟继友;吴晓蓉;;电源技术(第07期);全文 * |
直驱风电场接入弱电网宽频带振荡机理与抑制 方法(一):宽频带阻抗特性与振荡机理分析;李光辉 等;中 国 电 机 工 程 学 报;第39卷(第22期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114362210A (en) | 2022-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10452845B2 (en) | Generic framework to detect cyber threats in electric power grid | |
Reddy et al. | Condition monitoring of 11 kV distribution system insulators incorporating complex imagery using combined DOST-SVM approach | |
CN116703183B (en) | Novel power distribution network system power load assessment method | |
US20160079758A1 (en) | Method and device for monitoring and suppressing resonance | |
CN112834224B (en) | Nuclear power steam turbine generator health state assessment method and system | |
CN116505653A (en) | Transmission line monitoring system, method, device, computer equipment and storage medium | |
CN114362210B (en) | Wind farm oscillation risk assessment test method, avoidance method and storage medium | |
CN110867889B (en) | Oscillation stability judging method and system for wind power plant/unit connected with alternating current power grid | |
CN109613324A (en) | A kind of detection method and device of Harmonics amplification | |
CN112198375A (en) | Phase identification method, device, equipment and storage medium for single-phase user | |
CN114691662A (en) | Data quality inspection rule self-adaption method, storage medium and system | |
CN111062590A (en) | Electricity abnormal behavior detection method and device, computer equipment and storage medium | |
CN107370149B (en) | Method and system for determining load reduction amount of power system | |
Rao et al. | Fault recognition and diagnosis based on a decision tree for power distribution systems | |
CN111414582B (en) | Photovoltaic theoretical power calculation method, device, equipment and storage medium | |
CN112731022B (en) | Photovoltaic inverter fault detection method, equipment and medium | |
CN116760193B (en) | Comprehensive data monitoring method and system for electric power system | |
CN109888847B (en) | Optimized operation method, device and system for multi-feed-in alternating current-direct current hybrid power grid | |
CN113884889B (en) | Method and device for battery safety pre-warning, storage medium and electronic equipment | |
CN114389241B (en) | Relay protection setting value setting method and device for active power distribution network | |
CN116169690A (en) | Subsynchronous oscillation suppression method for wind farm grid-connected system based on stability control device | |
CN107565578B (en) | The reactive output control method and device of power distribution network | |
CN112686503A (en) | Evaluation method and system for asynchronous power grid frequency regulation quality | |
CN114047372B (en) | Voltage characteristic-based platform region topology identification system | |
Yang et al. | Research on Attack Identification Method and Device Method Based on Random Forest Algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |