A power system oscillation source positioning method, system, electronic device and medium

Abstract

The present application relates to a power system oscillation source positioning method, system, electronic equipment and medium, the oscillation source positioning method is through fast Fourier transform spectrum analysis to the instantaneous value waveform data of multiple key nodes is identified, obtains one or more leading interharmonic frequency, thereby obtaining interharmonic active power, based on interharmonic active power, the oscillation source of power system is identified and contribution degree sorting, directly and accurately locate the oscillation source under the target oscillation mode from multiple key nodes, without knowing the detailed mathematical model and parameter of system, especially suitable for modern power system containing a large number of inverter resources, and through contribution degree sorting, accurately quantifies the contribution degree of each key node under the target oscillation mode, provides accurate basis for subsequent suppression measures.

Description

Technical Field

[0001] This invention relates to the field of power system monitoring technology, specifically to a method, system, electronic device, and medium for locating power system oscillation sources. Background Technology

[0002] Power system oscillations are one of the main threats to small-signal stability. With the large-scale integration of inverter-type renewable energy sources, the system damping characteristics have changed, and broadband oscillation problems have become increasingly prominent. Quickly and accurately locating the oscillation source is a prerequisite for taking effective suppression measures.

[0003] Existing methods for locating oscillation sources have the following main shortcomings: Strong model dependence: The traditional participation factor method requires an accurate mathematical model and state equation of the system, which is computationally complex and not applicable to systems composed of inverters, making it difficult to apply online.

[0004] Limitations of the phasor method: Methods such as energy dissipation based on synchronous phasor measurement units are theoretically based on dynamic phasors. However, dynamic phasors are essentially a mixture of multiple frequency components. When the oscillation frequency is high (e.g., close to or exceeding the upper limit of the frequency reported by the PMU), aliasing effects can lead to distorted or even erroneous analysis results, failing to accurately reflect the oscillation energy of a specific mode.

[0005] Unable to quantify contribution: Most existing methods can only qualitatively determine the direction of oscillation energy, making it difficult to accurately quantify the "contribution" of each power generation device to a specific oscillation mode, and making it inconvenient to sort multiple oscillation sources and formulate precise control strategies. Summary of the Invention

[0006] To address the problems existing in the prior art, this invention proposes a method, system, electronic device, and medium for locating oscillation sources in power systems.

[0007] In a first aspect, the present invention provides a method for locating oscillation sources in a power system, the method comprising: Based on the instantaneous waveform data of three-phase voltage and three-phase current of multiple key nodes in the power system, fast Fourier transform spectrum analysis is used to identify one or more dominant interharmonic frequencies of each key node in the target oscillation mode. Based on one or more dominant interharmonic frequencies of each key node in the target oscillation mode, the interharmonic active power of each key node is obtained. Based on the interharmonic active power of each key node, the oscillation sources of the power system are identified and ranked by contribution, resulting in the oscillation source identification results and contribution ranking results. Based on the oscillation source identification results and the contribution ranking results, the oscillation source of the power system is located.

[0008] Optionally, the instantaneous waveform data of three-phase voltage and three-phase current based on multiple key nodes in the power system are identified using Fast Fourier Transform spectral analysis to obtain one or more dominant interharmonic frequencies of each key node in the target oscillation mode, including: Acquire instantaneous waveform data of three-phase voltage and three-phase current at multiple key nodes in the power system; Based on each of the key nodes, the instantaneous waveform data of the key nodes are preprocessed to obtain preprocessed instantaneous waveform data. Based on a preset time window, the preprocessed instantaneous waveform data is selected and truncated to obtain instantaneous waveform data of a fixed duration; Based on the instantaneous waveform data of the fixed duration, fast Fourier transform spectrum analysis is used to identify one or more dominant interharmonic frequencies of the key node in the target oscillation mode.

[0009] Optionally, the instantaneous waveform data based on the fixed duration is used to identify one or more dominant interharmonic frequencies of the key node in the target oscillation mode using fast Fourier transform spectral analysis, including: Perform a Fast Fourier Transform on the instantaneous waveform data of the fixed duration to obtain the complex voltage spectrum and complex current spectrum of each phase in the power system; Based on the complex voltage spectrum and the complex current spectrum, multiple interharmonic candidate frequencies of the key node in the target oscillation mode are identified; Interharmonic candidate frequencies whose difference between any interharmonic candidate frequency and the fundamental frequency is within the error range of the power oscillation frequency are determined as dominant interharmonic frequencies, so as to obtain one or more dominant interharmonic frequencies of the key node in the target oscillation mode.

[0010] Optionally, after determining the interharmonic candidate frequencies whose difference between any interharmonic candidate frequency and the fundamental frequency is within the error range of the power oscillation frequency as the dominant interharmonic frequency, the method further includes: If there are symmetrically distributed frequency components on both sides of the fundamental frequency, the symmetrically distributed frequency components are determined as a set of dominant interharmonic frequencies.

[0011] Optionally, the interharmonic active power includes: the positive and negative signs of the interharmonic active power and the magnitude of the interharmonic active power; The process of identifying and ranking the oscillation sources of the power system based on the interharmonic active power of each key node yields the oscillation source identification results and contribution ranking results, including: Based on the sign of the interharmonic active power of each key node, the oscillation source is identified for the multiple key nodes, and the oscillation source identification result is obtained. Based on the interharmonic active power of each oscillation source in the oscillation source identification results, the multiple oscillation sources in the oscillation source identification results are sorted by contribution to obtain the contribution ranking results.

[0012] Optionally, the step of identifying oscillation sources for the plurality of key nodes based on the sign of the interharmonic active power of each key node, and obtaining oscillation source identification results, includes: Key nodes whose interharmonic active power signs correspond to positive signs are identified as oscillation source nodes; All key nodes identified as oscillation source nodes, along with the interharmonic active power corresponding to those oscillation source nodes, are used as the oscillation source identification results.

[0013] Optionally, obtaining the interharmonic active power of each key node based on one or more dominant interharmonic frequencies of each key node in the target oscillation mode includes: Based on each of the key nodes, the voltage phasor and current phasor for each dominant interharmonic frequency are obtained; Based on the voltage phasor and current phasor of each dominant interharmonic frequency, the interharmonic active power of the key node is obtained.

[0014] Optionally, obtaining the voltage phasor and current phasor for each dominant interharmonic frequency based on each of the key nodes includes: Obtain the complex voltage spectrum and complex current spectrum for each dominant interharmonic frequency; Based on the complex voltage spectrum and complex current spectrum of each dominant interharmonic frequency, the voltage amplitude, voltage phase, current amplitude, and current phase of the dominant interharmonic frequency are calculated respectively. The voltage amplitude and voltage phase, as well as the current amplitude and current phase of each dominant interharmonic frequency, are respectively used as the voltage phasor and current phasor of each dominant interharmonic frequency.

[0015] Optionally, obtaining the interharmonic active power of the key node based on the voltage phasor and current phasor of each dominant interharmonic frequency includes: Based on the voltage phasor and current phasor of each dominant interharmonic frequency, the available interharmonic active power of each dominant interharmonic frequency is calculated. Based on each key node, the available interharmonic active power of all dominant interharmonic frequencies of the key node is summed to obtain the interharmonic active power of the key node.

[0016] Optionally, after locating the oscillation source of the power system based on the oscillation source identification result and the contribution ranking result, the method further includes: The key nodes corresponding to the located oscillation sources are used as monitoring points; Obtain the time-series variation sequence of the intermediate harmonic active power at the monitoring point; Based on the time-series variation sequence of the intermediate harmonic active power at the monitoring points, the trend of the interharmonic active power is obtained. Based on the interharmonic active power trend, it is determined whether there is a risk of increased oscillation at the monitoring point; When there is a risk of increased oscillation at the monitoring point, an oscillation warning signal is issued to provide early warning to the monitoring point.

[0017] Optionally, determining whether there is a risk of increased oscillation at the monitoring point based on the interharmonic active power trend includes: When the trend of the interharmonic active power is greater than a preset threshold, or when the interharmonic active power shows an increasing trend within a preset time period, it is determined that there is a risk of increased oscillation at the monitoring point.

[0018] A second aspect of the present invention also provides a power system oscillation source location system, the system comprising: The dominant interharmonic frequency determination unit is used to identify one or more dominant interharmonic frequencies of each key node in the target oscillation mode by using fast Fourier transform spectrum analysis based on the instantaneous waveform data of three-phase voltage and three-phase current of multiple key nodes in the power system. An interharmonic active power determination unit is used to obtain the interharmonic active power of each key node based on one or more dominant interharmonic frequencies of each key node in the target oscillation mode. The identification and ranking result determination unit is used to identify and rank the oscillation sources of the power system based on the interharmonic active power of each key node, and obtain the oscillation source identification result and the contribution ranking result. A power system oscillation source localization unit is used to locate the oscillation source of the power system based on the oscillation source identification result and the contribution ranking result.

[0019] Optionally, the dominant interharmonic frequency determination unit includes: The instantaneous waveform data acquisition module is used to acquire instantaneous waveform data of three-phase voltage and three-phase current at multiple key nodes in the power system. The instantaneous waveform data preprocessing module is used to preprocess the instantaneous waveform data of each key node to obtain preprocessed instantaneous waveform data. The instantaneous waveform data extraction module is used to select and extract preprocessed instantaneous waveform data based on a preset time window to obtain instantaneous waveform data of a fixed duration. The first dominant interharmonic frequency determination module is used to identify one or more dominant interharmonic frequencies of the key node in the target oscillation mode based on the instantaneous waveform data of the fixed duration using fast Fourier transform spectrum analysis.

[0020] Optionally, the first dominant interharmonic frequency determination module is specifically used for: Perform a Fast Fourier Transform on the instantaneous waveform data of the fixed duration to obtain the complex voltage spectrum and complex current spectrum of each phase in the power system; Based on the complex voltage spectrum and the complex current spectrum, multiple interharmonic candidate frequencies of the key node in the target oscillation mode are identified; Interharmonic candidate frequencies whose difference between any interharmonic candidate frequency and the fundamental frequency is within the error range of the power oscillation frequency are determined as dominant interharmonic frequencies, so as to obtain one or more dominant interharmonic frequencies of the key node in the target oscillation mode.

[0021] Optionally, the dominant interharmonic frequency determination unit further includes: The second dominant interharmonic frequency determination module is used to determine the symmetrically distributed frequency components as a set of dominant interharmonic frequencies if there are symmetrically distributed frequency components on both sides of the fundamental frequency.

[0022] Optionally, the interharmonic active power includes: the positive and negative signs of the interharmonic active power and the magnitude of the interharmonic active power; The identification and sorting result determination unit includes: The oscillation source identification result determination module is used to identify the oscillation source of the multiple key nodes based on the positive and negative signs of the interharmonic active power of each key node, and obtain the oscillation source identification result. The contribution ranking result determination module is used to rank multiple oscillators in the oscillator identification result based on the interharmonic active power of each oscillator in the oscillator identification result, and obtain the contribution ranking result.

[0023] Optionally, the oscillation source identification result determination module is specifically used for: Key nodes whose interharmonic active power signs correspond to positive signs are identified as oscillation source nodes; All key nodes identified as oscillation source nodes, along with the interharmonic active power corresponding to those oscillation source nodes, are used as the oscillation source identification results.

[0024] Optionally, the interharmonic active power determination unit includes: The phasor determination module is used to obtain the voltage phasor and current phasor for each dominant interharmonic frequency based on each of the key nodes. The interharmonic active power determination module is used to obtain the interharmonic active power of the key node based on the voltage phasor and current phasor of each dominant interharmonic frequency.

[0025] Optionally, the phasor determination module is specifically used for: Obtain the complex voltage spectrum and complex current spectrum for each dominant interharmonic frequency; Based on the complex voltage spectrum and complex current spectrum of each dominant interharmonic frequency, the voltage amplitude, voltage phase, current amplitude, and current phase of the dominant interharmonic frequency are calculated respectively. The voltage amplitude and voltage phase, as well as the current amplitude and current phase of each dominant interharmonic frequency, are respectively used as the voltage phasor and current phasor of each dominant interharmonic frequency.

[0026] Optionally, the interharmonic active power determination module is specifically used for: Based on the voltage phasor and current phasor of each dominant interharmonic frequency, the available interharmonic active power of each dominant interharmonic frequency is calculated. Based on each key node, the available interharmonic active power of all dominant interharmonic frequencies of the key node is summed to obtain the interharmonic active power of each key node.

[0027] Optionally, the system further includes: The detection point determination unit is used to identify the key nodes corresponding to the located oscillation source as monitoring points. A time-series change sequence determination unit is used to obtain the time-series change sequence of the intermediate harmonic active power at the monitoring point; The interharmonic active power trend determination unit is used to obtain the interharmonic active power trend based on the time-series change sequence of the interharmonic active power at the monitoring point. The oscillation amplification risk judgment unit is used to determine whether there is a risk of oscillation amplification at the monitoring point based on the interharmonic active power trend. The detection point early warning unit is used to issue an oscillation early warning signal when there is a risk of increased oscillation at the monitoring point, thereby providing an early warning to the monitoring point.

[0028] Optionally, the oscillation aggravation risk judgment unit is used for: When the trend of the interharmonic active power is greater than a preset threshold, or when the interharmonic active power shows an increasing trend within a preset time period, it is determined that there is a risk of increased oscillation at the monitoring point.

[0029] A third aspect of the present invention also provides an electronic device, comprising: one or more processors; A processor is used to execute one or more programs; When the one or more programs are executed by the one or more processors, a power system oscillation source localization method as described above is implemented.

[0030] In a fourth aspect, the present invention also provides a computer-readable storage medium having a computer program stored thereon, wherein when the computer program is executed, it implements a power system oscillation source localization method as described above.

[0031] Compared with the prior art, the beneficial effects of the present invention are as follows: This invention provides a method, system, electronic device, and medium for locating oscillation sources in a power system. The method identifies instantaneous waveform data from multiple key nodes using Fast Fourier Transform (FFT) spectral analysis to obtain one or more dominant interharmonic frequencies, thereby determining the interharmonic active power. Based on the interharmonic active power, the oscillation sources in the power system are identified and ranked by contribution. This allows for direct and precise location of the oscillation source within the target oscillation mode from multiple key nodes, without requiring knowledge of the system's detailed mathematical model and parameters. It is particularly suitable for modern power systems with a large number of inverter-type resources. Furthermore, by ranking the contribution, the method accurately quantifies the contribution of each key node within the target oscillation mode, providing a precise basis for subsequent suppression measures. Attached Figure Description

[0032] Figure 1 A flowchart of a method for locating oscillation sources in a power system provided by the present invention; Figure 2 A block diagram of a power system oscillation source localization system provided by the present invention; Figure 3 A block diagram of an electronic device provided by the present invention. Detailed Implementation

[0033] Example 1: Figure 1 A flowchart of a power system oscillation source localization method provided by the present invention is shown below. Figure 1 As shown, the method may include the following steps 101 to 104: In step 101, based on the instantaneous waveform data of three-phase voltage and three-phase current of multiple key nodes in the power system, fast Fourier transform spectrum analysis is used to identify one or more dominant interharmonic frequencies of each key node in the target oscillation mode.

[0034] In step 102, the interharmonic active power of each key node is obtained based on one or more dominant interharmonic frequencies of each key node in the target oscillation mode.

[0035] In step 103, based on the interharmonic active power of each key node, the oscillation sources of the power system are identified and ranked by contribution, resulting in oscillation source identification results and contribution ranking results.

[0036] In step 104, the oscillation source of the power system is located based on the oscillation source identification result and the contribution ranking result.

[0037] The interharmonic active power includes: the sign of the interharmonic active power and its magnitude. The sign of the interharmonic active power is used to characterize the energy flow direction of the node during oscillation, and its absolute value is used to characterize the intensity of the node's participation in oscillation. The dominant interharmonic frequency can be contained in a pair of frequency components symmetrically distributed on both sides of the fundamental frequency.

[0038] above Figure 1 Possible implementations of step 101 shown may include the following steps 1011 to 1014: In step 1011, instantaneous waveform data of three-phase voltage and three-phase current at multiple key nodes in the power system are acquired.

[0039] In step 1012, based on each key node, the instantaneous waveform data of the key node is preprocessed to obtain preprocessed instantaneous waveform data.

[0040] In step 1013, based on a preset time window, the preprocessed instantaneous waveform data is selected and truncated to obtain instantaneous waveform data of a fixed duration.

[0041] In step 1014, based on the instantaneous waveform data of the fixed duration, fast Fourier transform spectrum analysis is used to identify one or more dominant interharmonic frequencies of the key node in the target oscillation mode.

[0042] The preprocessing may include timestamp correction, removal of abnormal acquisition points, elimination of DC components, amplitude calibration, phase calibration, anti-aliasing filtering, noise suppression, and three-phase channel consistency verification.

[0043] It should be noted that after obtaining the spectrum by performing FFT, the following procedure is used to identify the dominant interharmonics corresponding to the target oscillation mode: Pattern matching method: 1) Determine the target oscillation frequency: The frequency of the oscillation mode existing in the system can be initially obtained by performing fast Fourier transform or Prony analysis on the power, frequency and other signals provided by the traditional phasor measurement unit (PMU). 2) Spectrum Search: Directly search the FFT spectrum of voltage or current. The search amplitude within the nearby narrow band is significantly higher than the spectral peak of the background noise. These two spectral peaks should be roughly symmetrical and have similar amplitudes (in the absence of severe distortion).

[0044] Automatic search method (when the target pattern is unknown): 1) Symmetrical Spectral Peak Detection: Performs peak detection on the FFT spectrum, automatically finding all paired peaks related to... Symmetrical spectral peak pairs and satisfy .

[0045] 2) Energy significance assessment: Calculate the sum of squares of the amplitudes of the pair of spectral peaks as the "energy index" of the potential mode. Sort the energy indices and eliminate spurious modes with excessively low energy (e.g., below a certain percentage of the total signal energy).

[0046] 3) Mode frequency calculation: Determined dominant interharmonic pairs Afterwards, the oscillation mode frequency is .

[0047] The possible implementation of step 1014 above may include the following steps S1 to S3: In step S1, a fast Fourier transform is performed on the instantaneous waveform data of the fixed duration to obtain the complex voltage spectrum and complex current spectrum of each phase in the power system.

[0048] In step S2, based on the complex voltage spectrum and the complex current spectrum, multiple interharmonic candidate frequencies of the key node in the target oscillation mode are identified.

[0049] In step S3, interharmonic candidate frequencies whose difference between any interharmonic candidate frequency and the fundamental frequency is within the error range of the power oscillation frequency are determined as dominant interharmonic frequencies, so as to obtain one or more dominant interharmonic frequencies of the key node in the target oscillation mode.

[0050] In some implementations, after step S3, the method may further include the following step S4: In step S4, if there are symmetrically distributed frequency components on both sides of the fundamental frequency, the symmetrically distributed frequency components are determined as a set of dominant interharmonic frequencies.

[0051] It should be noted that power oscillations in a power system typically manifest as a pair of symmetrical interharmonic components near the power frequency. Their physical essence is the amplitude-frequency modulation (AM-FM) effect of the oscillation frequency on the power frequency signal.

[0052] If the system has a frequency of f m A power oscillation mode (unit: Hz) will manifest in the voltage and current waveforms as a pair of oscillations oscillating at the system power frequency. f Interharmonic components symmetrically distributed with 0 (e.g., 50Hz) as the center: Top Sideband: f 0+ f m Lower sideband: f 0- f m For example, if a power oscillation of 1.8 Hz is observed, the corresponding dominant interharmonic frequencies in a 50 Hz system are 48.2 Hz and 51.8 Hz.

[0053] above Figure 1 Possible implementations of step 103 shown may include the following steps 1031 to 1032: In step 1031, based on the positive and negative signs of the interharmonic active power of each key node, the multiple key nodes are identified as oscillation sources to obtain the oscillation source identification results.

[0054] In step 1032, based on the interharmonic active power of each oscillation source in the oscillation source identification result, the multiple oscillation sources in the oscillation source identification result are sorted by contribution to obtain the contribution ranking result.

[0055] One possible implementation of step 1031 above may include: identifying key nodes whose interharmonic active power has a positive sign as oscillation source nodes. All key nodes identified as oscillation source nodes, and the corresponding interharmonic active power of each oscillation source node, are used as the oscillation source identification result.

[0056] above Figure 1 Possible implementations of step 102 shown may include the following steps 1021 to 1022: In step 1021, based on each of the key nodes, the voltage phasor and current phasor for each dominant interharmonic frequency are obtained.

[0057] In step 1022, the interharmonic active power of the key node is obtained based on the voltage phasor and current phasor of each dominant interharmonic frequency.

[0058] One possible implementation of step 1021 above may include: obtaining the complex voltage spectrum and complex current spectrum for each dominant interharmonic frequency. Based on the complex voltage spectrum and complex current spectrum for each dominant interharmonic frequency, calculating the voltage amplitude, voltage phase, current amplitude, and current phase for each dominant interharmonic frequency. Using the voltage amplitude and voltage phase, as well as the current amplitude and current phase, for each dominant interharmonic frequency as the voltage phasor and current phasor for each dominant interharmonic frequency, respectively.

[0059] Possible implementations of step 1022 above may include: calculating the available interharmonic active power for each dominant interharmonic frequency based on the voltage phasor and current phasor of each dominant interharmonic frequency; and summing the available interharmonic active power for all dominant interharmonic frequencies of each critical node to obtain the interharmonic active power of the critical node.

[0060] It should be noted that identifying the dominant interharmonic frequencies Next, the amplitude and phase of this frequency component need to be accurately extracted for power calculation. This invention recommends the following high-precision method, which surpasses simple FFT spectral interpolation: Windowed interpolation FFT method: 1) Apply a Hanning or Blackman window to data truncated to integer power frequency cycles to suppress spectral leakage.

[0061] 2) Perform an FFT on the windowed data to find... The spectral line with the largest amplitude in the vicinity and its adjacent spectral lines to the left and right.

[0062] 3) By using bi- or tri-spectral line interpolation algorithms, frequency, amplitude, and phase errors caused by the picket fence effect can be corrected, and frequency and phasor information accurate to 0.01 Hz can be obtained.

[0063] Processing methods based on analytic signals: 1) Perform Hilbert transform on the original signal to construct its analytic signal.

[0064] 2) Design a center frequency of A complex bandpass filter is used to filter analytic signals to directly obtain the signal. The complex envelope signal is a carrier wave.

[0065] 3) The amplitude of this complex envelope is the amplitude of the interharmonic voltage or current. or Its phase angle (relative to a fixed time reference point) is the phase. or This method has strong anti-interference capabilities and is especially suitable for noisy environments.

[0066] Optionally, in the above Figure 1 Following step 104, the method may further include the following steps 105 to 109: In step 105, the key node corresponding to the located oscillation source is used as a monitoring point.

[0067] In step 106, the time-series variation sequence of the intermediate harmonic active power at the monitoring point is obtained.

[0068] In step 107, the trend of interharmonic active power is obtained based on the time-series change sequence of the intermediate harmonic active power at the monitoring point.

[0069] In step 108, based on the interharmonic active power trend, it is determined whether there is a risk of increased oscillation at the monitoring point.

[0070] In step 109, when there is a risk of increased oscillation at the monitoring point, an oscillation warning signal is issued to warn the monitoring point.

[0071] One possible implementation of step 108 above is as follows: when the trend of the interharmonic active power is greater than a preset threshold, or when the interharmonic active power shows an increasing trend within a preset time period, it is determined that there is a risk of increased oscillation at the monitoring point.

[0072] It should be noted that in complex scenarios, a single node may participate in multiple oscillation modes simultaneously. This system possesses multi-mode parallel analysis capabilities: For each identified independent oscillation mode (corresponding to a pair of interharmonics), calculate the interharmonic active power P_IH of each node under that mode.

[0073] Contribution ranking is based on pattern. For example, in pattern A, node X contributes +100kW and node Y contributes +50kW; in pattern B, node Y may contribute +80kW and node X may contribute -10kW.

[0074] This analysis provides operators with a clear picture, distinguishing the sources and sinks of different oscillation modes, thereby enabling the development of more targeted control strategies.

[0075] In some implementations, the power system oscillation source location method can also be referred to as the power system oscillation source location and early warning method based on interharmonic power.

[0076] This invention provides a method for locating and warning of power system oscillation sources based on interharmonic power, comprising the following steps: S1: Simultaneously collect instantaneous waveform data of three-phase voltage and current at multiple key nodes in the power grid; S2: Preprocess the acquired waveform data.

[0077] S3: Select and extract data within a fixed-duration time window; the width of the time window is an integer multiple of the power frequency cycle.

[0078] S4: Perform fast Fourier transform spectrum analysis on the data within the time window to identify one or more dominant interharmonic frequencies corresponding to the target oscillation mode.

[0079] S5: For each of the dominant interharmonic frequencies, accurately extract its voltage phasor and current phasor.

[0080] S6: For each of the dominant interharmonic frequencies, calculate the corresponding interharmonic active power.

[0081] S7: Based on the direction and magnitude of the interharmonic active power, locate the oscillation source and rank its contribution. If the interharmonic active power is positive, the equipment at the monitoring point is determined to be injecting oscillation energy into the power grid and is thus an oscillation source; the larger the power value, the higher its contribution to the oscillation.

[0082] S8: Continue executing steps S3 to S7 to track the trend of the interharmonic active power. When the interharmonic active power value at a certain monitoring point exceeds the preset threshold or continues to increase, an oscillation warning signal is issued.

[0083] Furthermore, in step S3, the width of the time window is 5 to 20 power frequency cycles.

[0084] Furthermore, in step S4, the method for identifying the dominant interharmonic frequency is as follows: the difference between the dominant interharmonic frequency and the fundamental frequency is equal to the power oscillation frequency observed in the phasor domain.

[0085] Furthermore, in step S6, the formula for calculating the interharmonic active power is: P_IH = V_IH * I_IH * cos(θ_VI) Where V_IH and I_IH are the voltage amplitude and current amplitude at the interharmonic frequency, respectively, and θ_VI is the phase difference between the voltage phasor and the current phasor at the frequency.

[0086] The present invention has the following advantages: directness and physical clarity: it directly measures the essential physical quantity that sustains the oscillation—interharmonic active power, with clear physical meaning, avoiding the theoretical defects and aliasing problems of indirect methods based on phasors.

[0087] Model-independent: No need to know the detailed mathematical model and parameters of the system, making it particularly suitable for modern power systems with a large number of inverter resources, thus achieving true "data-driven" operation.

[0088] Precise quantification: It can not only determine the direction of oscillation energy, but also accurately quantify the contribution of each oscillation source through power values, providing precise basis for operators to take suppression measures (such as switching off the machine or adjusting controller parameters).

[0089] Highly implementable online: The calculation process is clear, the computational load is moderate, and it can be embedded in existing power quality monitoring devices or PMUs to achieve online, real-time oscillation monitoring and early warning.

[0090] Early warning: By tracking the trend of interharmonic power, signs of instability can be identified before the oscillation amplitude becomes dangerous, thus enabling early warning.

[0091] For example, in a certain regional power grid, a sustained power oscillation with a frequency of approximately 1.8 Hz was observed. To locate the source of the oscillation, the following steps were performed: 1) Data acquisition: Install or utilize existing high-performance power quality monitoring devices at the grid connection points of wind farm A, wind farm B and thermal power plant C in the region to synchronously acquire voltage and current waveforms at a sampling rate of 3.2kHz.

[0092] 2) Window processing: Extract a data window with a length of 10 power frequency cycles (approximately 200ms). Apply the Hanning window function to this data window to reduce spectral leakage.

[0093] 3) Spectrum Analysis: FFT was performed on the windowed data. Obvious interharmonic components were observed in the spectrum, with frequencies of 48.2 Hz and 51.8 Hz. The difference between these two frequencies and the fundamental frequency of 50 Hz is exactly 1.8 Hz, therefore they were identified as the dominant interharmonic frequencies corresponding to the 1.8 Hz oscillation mode.

[0094] 4) Power Calculation: Extract the voltage and current phasors (amplitude and phase) at 48.2Hz and 51.8Hz respectively. Calculate the interharmonic active power at these two frequencies respectively, and add them together to obtain the total interharmonic active power P_IH_total for each monitoring point.

[0095] 5) Location and Ranking: The calculation results are as follows: Wind farm A's P_IH_total = +85 kW, wind farm B's P_IH_total = -15 kW, and thermal power plant C's P_IH_total = +120 kW. Based on these results, it is determined that thermal power plant C and wind farm A are oscillation sources, with thermal power plant C contributing more than wind farm A; wind farm B is a participant in the oscillation, absorbing some of the oscillation energy.

[0096] 6) Early Warning: Continuous monitoring revealed that the P_IH_total of thermal power plant C increased from 120kW to 180kW within one minute. The system immediately issued a Level 1 early warning to the dispatcher, indicating a risk of escalating oscillations and recommending that the excitation system controller parameters of thermal power plant C be checked and adjusted first.

[0097] In a certain regional power grid, the dispatching master station detected a continuous power oscillation of approximately 1.8 Hz. To quickly determine the location of the oscillation source and assess the contribution of different grid-connected devices to the oscillation, high-performance power quality monitoring devices were deployed at the grid connection points of wind farm A, wind farm B, and thermal power plant C within the region. Alternatively, synchronous monitoring devices with waveform acquisition capabilities were activated to simultaneously collect instantaneous waveform data of the three-phase voltage and three-phase current at each node. The devices used a unified time synchronization method, with a sampling rate set to 3.2 kHz to ensure the accuracy of identifying interharmonic components near the target frequency.

[0098] After acquiring the raw waveforms, preprocessing is performed first, including outlier removal, DC bias elimination, channel calibration, and time alignment. Then, waveform data is truncated into 10 power frequency cycles as an analysis time window, and a Hanning window function is applied to each time window to reduce spectral leakage caused by finite data truncation. Fast Fourier Transforms are then performed on the windowed three-phase voltage and three-phase current data to obtain the corresponding spectral results.

[0099] Since the frequency of the target oscillation mode observed in the phasor domain is approximately 1.8 Hz, while the system fundamental frequency is 50 Hz, the focus is on searching for interharmonic components near 50 Hz that correspond to a 1.8 Hz offset. Spectral analysis results show significant frequency peaks at 48.2 Hz and 51.8 Hz, and the difference between these two frequencies and 50 Hz is 1.8 Hz. Therefore, 48.2 Hz and 51.8 Hz are identified as the dominant interharmonic frequencies corresponding to the target oscillation mode. Furthermore, voltage and current phasors are extracted from the corresponding frequencies to obtain the voltage amplitude, current amplitude, and phase difference at each monitoring point at 48.2 Hz and 51.8 Hz.

[0100] For each dominant interharmonic frequency, the interharmonic active power is calculated using the formula P_IH=V_IH×I_IH×cos(θ_VI). The interharmonic active power corresponding to 48.2Hz and 51.8Hz is then summed to obtain the total interharmonic active power P_IH_total for each monitoring point in this 1.8Hz oscillation mode. The calculation results show that P_IH_total for wind farm A is +85kW, for wind farm B it is -15kW, and for thermal power plant C it is +120kW. This indicates that both thermal power plant C and wind farm A inject energy related to this oscillation mode into the grid, acting as oscillation sources. Thermal power plant C injects more power and contributes more significantly. Wind farm B, however, shows a negative value, indicating that it is primarily a participating node in the oscillation and absorbs some of the oscillation energy.

[0101] During continuous monitoring, the system used a sliding time window to continuously track the interharmonic active power at each monitoring point. The results showed that the P_IH_total of thermal power plant C increased from 120kW to 180kW within one minute, and this upward trend was observed across multiple analysis windows. Based on this, the early warning module determined that there was a risk of increased oscillation at this node, issued a Level 1 early warning signal to the dispatchers, and simultaneously displayed the oscillation source ranking, power trend curve, and suggested priority inspection targets through the human-machine interface. The system indicated that relevant control parameters of thermal power plant C should be checked first, such as the excitation system or other control links related to low-frequency power oscillations.

[0102] Example 2: Figure 2 A block diagram of a power system oscillation source localization system provided by the present invention is shown below. Figure 2 As shown, the system includes: The dominant interharmonic frequency determination unit is used to identify one or more dominant interharmonic frequencies of each key node in the target oscillation mode by using fast Fourier transform spectrum analysis based on the instantaneous waveform data of three-phase voltage and three-phase current of multiple key nodes in the power system. An interharmonic active power determination unit is used to obtain the interharmonic active power of each key node based on one or more dominant interharmonic frequencies of each key node in the target oscillation mode. The identification and ranking result determination unit is used to identify and rank the oscillation sources of the power system based on the interharmonic active power of each key node, and obtain the oscillation source identification result and the contribution ranking result. A power system oscillation source localization unit is used to locate the oscillation source of the power system based on the oscillation source identification result and the contribution ranking result.

[0103] Optionally, the dominant interharmonic frequency determination unit includes: The instantaneous waveform data acquisition module is used to acquire instantaneous waveform data of three-phase voltage and three-phase current at multiple key nodes in the power system. The instantaneous waveform data preprocessing module is used to preprocess the instantaneous waveform data of each key node to obtain preprocessed instantaneous waveform data. The instantaneous waveform data extraction module is used to select and extract preprocessed instantaneous waveform data based on a preset time window to obtain instantaneous waveform data of a fixed duration. The first dominant interharmonic frequency determination module is used to identify one or more dominant interharmonic frequencies of the key node in the target oscillation mode based on the instantaneous waveform data of the fixed duration using fast Fourier transform spectrum analysis.

[0104] Optionally, the first dominant interharmonic frequency determination module is specifically used for: Perform a Fast Fourier Transform on the instantaneous waveform data of the fixed duration to obtain the complex voltage spectrum and complex current spectrum of each phase in the power system; Based on the complex voltage spectrum and the complex current spectrum, multiple interharmonic candidate frequencies of the key node in the target oscillation mode are identified; Interharmonic candidate frequencies whose difference from any interharmonic candidate frequency to the fundamental frequency is within the error range of the power oscillation frequency are determined as dominant interharmonic frequencies, so as to obtain one or more dominant interharmonic frequencies of the key node in the target oscillation mode.

[0105] Optionally, the dominant interharmonic frequency determination unit further includes: The second dominant interharmonic frequency determination module is used to determine the symmetrically distributed frequency components as a set of dominant interharmonic frequencies if there are symmetrically distributed frequency components on both sides of the fundamental frequency.

[0106] Optionally, the interharmonic active power includes: the positive and negative signs of the interharmonic active power and the magnitude of the interharmonic active power; The identification and sorting result determination unit includes: The oscillation source identification result determination module is used to identify the oscillation source of the multiple key nodes based on the positive and negative signs of the interharmonic active power of each key node, and obtain the oscillation source identification result. The contribution ranking result determination module is used to rank multiple oscillators in the oscillator identification result based on the interharmonic active power of each oscillator in the oscillator identification result, and obtain the contribution ranking result.

[0107] Optionally, the oscillation source identification result determination module is specifically used for: Key nodes whose interharmonic active power signs correspond to positive signs are identified as oscillation source nodes; All key nodes identified as oscillation source nodes, along with the interharmonic active power corresponding to those oscillation source nodes, are used as the oscillation source identification results.

[0108] Optionally, the interharmonic active power determination unit includes: The phasor determination module is used to obtain the voltage phasor and current phasor for each dominant interharmonic frequency based on each of the key nodes. The interharmonic active power determination module is used to obtain the interharmonic active power of the key node based on the voltage phasor and current phasor of each dominant interharmonic frequency.

[0109] Optionally, the phasor determination module is specifically used for: Obtain the complex voltage spectrum and complex current spectrum for each dominant interharmonic frequency; Based on the complex voltage spectrum and complex current spectrum of each dominant interharmonic frequency, the voltage amplitude, voltage phase, current amplitude, and current phase of the dominant interharmonic frequency are calculated respectively. The voltage amplitude and voltage phase, as well as the current amplitude and current phase of each dominant interharmonic frequency, are respectively used as the voltage phasor and current phasor of each dominant interharmonic frequency.

[0110] Optionally, the interharmonic active power determination module is specifically used for: Based on the voltage phasor and current phasor of each dominant interharmonic frequency, the available interharmonic active power of each dominant interharmonic frequency is calculated. Based on each key node, the available interharmonic active power of all dominant interharmonic frequencies of the key node is summed to obtain the interharmonic active power of each key node.

[0111] Optionally, the system further includes: The detection point determination unit is used to identify the key nodes corresponding to the located oscillation source as monitoring points. A time-series change sequence determination unit is used to obtain the time-series change sequence of the intermediate harmonic active power at the monitoring point; The interharmonic active power trend determination unit is used to obtain the interharmonic active power trend based on the time-series change sequence of the interharmonic active power at the monitoring point. The oscillation amplification risk judgment unit is used to determine whether there is a risk of oscillation amplification at the monitoring point based on the interharmonic active power trend. The detection point early warning unit is used to issue an oscillation early warning signal when there is a risk of increased oscillation at the monitoring point, thereby providing an early warning to the monitoring point.

[0112] Optionally, the oscillation aggravation risk judgment unit is used for: When the trend of the interharmonic active power is greater than a preset threshold, or when the interharmonic active power shows an increasing trend within a preset time period, it is determined that there is a risk of increased oscillation at the monitoring point.

[0113] It should be noted that the block diagram illustrates the four core modules of the system: data acquisition module, data processing and analysis module, positioning and sorting module, and early warning module, as well as the data flow relationships between them.

[0114] 1. Data Acquisition Module Composition: High-performance power quality monitoring devices or PMUs with waveform acquisition capabilities are deployed at key nodes of the power grid (such as power plants, wind farms, and substations).

[0115] Function: Synchronously acquires instantaneous waveform data of three-phase voltage and current, providing high-fidelity raw data for subsequent analysis. Supports GPS / IRIG-B and other time synchronization methods to ensure data synchronization.

[0116] 2. Data Processing and Analysis Module Function: This is the core of the system's calculations. It receives the raw waveform and outputs key interharmonic power parameters.

[0117] Sub-steps: Data synchronization and preprocessing: Align the timestamps of data from each monitoring point, and perform filtering, calibration, etc.

[0118] Windowing and FFT Spectrum Analysis: Apply an integer power frequency period window (such as the Hanning window) to the data and perform FFT to transform the signal from the time domain to the frequency domain.

[0119] Interharmonic frequency identification: Automatically searches for and locks onto the dominant interharmonic frequencies (such as 58.2Hz, 61.8Hz) associated with the target oscillation mode in the spectrum.

[0120] Interharmonic power calculation: Accurately extract the voltage and current phasors of the above frequencies, and calculate the active power P_IH for each frequency.

[0121] 3. Intelligent Analysis Layer 1) Positioning and sorting module: Input: P_IH of each monitoring point.

[0122] Logic: Determine the power direction based on the sign of P_IH (positive value is the source, negative value is the participant); rank the oscillation sources by contribution based on the magnitude of P_IH.

[0123] Output: A list of ranked oscillation sources.

[0124] 2) Early warning module: Input: Time-series historical data of P_IH.

[0125] Logic: Track the changing trend of P_IH, and trigger the corresponding level of warning when it exceeds the static threshold or shows a continuous growth trend.

[0126] Output: Early warning signal.

[0127] 4. Application Presentation Layer 1) Human-computer interaction interface: Functions: Visualizes the location results for operators (e.g., highlighting the oscillation source on a map), provides contribution ranking bar charts, power trend curves, and early warning pop-ups. This is the ultimate manifestation of the system's value.

[0128] 2) Control system (optional): Function: Receive instructions from the positioning and sorting module, and automatically or assistedly execute suppression measures, such as sending instructions to adjust controller parameters to the oscillation source with the largest contribution, or performing a machine switching operation in an emergency.

[0129] Example 3: Based on the same inventive concept, the present invention also provides an electronic device, such as... Figure 3 As shown, the electronic device includes a processor and a memory. The memory stores a computer program, which includes program instructions. The processor executes the program instructions stored in the computer storage medium. The processor may be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. It is the computing and control core of the terminal, and is suitable for implementing one or more instructions. Specifically, it is suitable for loading and executing one or more instructions in the computer storage medium to implement the corresponding method flow or corresponding function, so as to implement the steps of the power system oscillation source localization method in the above embodiment.

[0130] Example 4: Based on the same inventive concept, this invention also provides a storage medium, specifically a computer-readable storage medium (Memory), which is a memory device in an electronic device used to store programs and data. It is understood that the computer-readable storage medium here can include both built-in storage media in the electronic device and extended storage media supported by the electronic device. The computer-readable storage medium provides storage space that stores the terminal's operating system. Furthermore, this storage space also stores one or more instructions suitable for loading and execution by a processor. These instructions can be one or more computer programs (including program code). It should be noted that the computer-readable storage medium here can be high-speed RAM or non-volatile memory, such as at least one disk storage device. The processor can load and execute one or more instructions stored in the computer-readable storage medium to implement the steps of the power system oscillation source localization method in the above embodiments.

[0131] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0132] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0133] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1The function specified in one or more boxes.

[0134] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0135] The above are merely embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention are included within the scope of the claims of the present invention pending approval.

Claims

1. A method for locating oscillation sources in a power system, characterized in that, The method includes: Based on the instantaneous waveform data of three-phase voltage and three-phase current of multiple key nodes in the power system, fast Fourier transform spectrum analysis is used to identify one or more dominant interharmonic frequencies of each key node in the target oscillation mode. Based on one or more dominant interharmonic frequencies of each key node in the target oscillation mode, the interharmonic active power of each key node is obtained. Based on the interharmonic active power of each key node, the oscillation sources of the power system are identified and ranked by contribution, resulting in the oscillation source identification results and contribution ranking results. Based on the oscillation source identification results and the contribution ranking results, the oscillation source of the power system is located.

2. The method according to claim 1, characterized in that, The instantaneous waveform data of three-phase voltage and three-phase current based on multiple key nodes in the power system are identified using Fast Fourier Transform spectral analysis to obtain one or more dominant interharmonic frequencies of each key node in the target oscillation mode, including: Acquire instantaneous waveform data of three-phase voltage and three-phase current at multiple key nodes in the power system; Based on each of the key nodes, the instantaneous waveform data of the key nodes are preprocessed to obtain preprocessed instantaneous waveform data. Based on a preset time window, the preprocessed instantaneous waveform data is selected and truncated to obtain instantaneous waveform data of a fixed duration; Based on the instantaneous waveform data of the fixed duration, fast Fourier transform spectrum analysis is used to identify one or more dominant interharmonic frequencies of the key node in the target oscillation mode.

3. The method according to claim 2, characterized in that, The instantaneous waveform data based on the fixed duration is identified using Fast Fourier Transform spectral analysis to obtain one or more dominant interharmonic frequencies of the key node in the target oscillation mode, including: Perform a Fast Fourier Transform on the instantaneous waveform data of the fixed duration to obtain the complex voltage spectrum and complex current spectrum of each phase in the power system; Based on the complex voltage spectrum and the complex current spectrum, multiple interharmonic candidate frequencies of the key node in the target oscillation mode are identified; Interharmonic candidate frequencies whose difference from any interharmonic candidate frequency to the fundamental frequency is within the error range of the power oscillation frequency are determined as dominant interharmonic frequencies, so as to obtain one or more dominant interharmonic frequencies of the key node in the target oscillation mode.

4. The method according to claim 3, characterized in that, After determining the interharmonic candidate frequencies whose difference between any interharmonic candidate frequency and the fundamental frequency is within the error range of the power oscillation frequency as the dominant interharmonic frequency, the method further includes: If there are symmetrically distributed frequency components on both sides of the fundamental frequency, the symmetrically distributed frequency components are determined as a set of dominant interharmonic frequencies.

5. The method according to claim 1, characterized in that, The interharmonic active power includes: the positive and negative signs of the interharmonic active power and the magnitude of the interharmonic active power; The process of identifying and ranking the oscillation sources of the power system based on the interharmonic active power of each key node yields the oscillation source identification results and contribution ranking results, including: Based on the sign of the interharmonic active power of each key node, the oscillation source is identified for the multiple key nodes, and the oscillation source identification result is obtained. Based on the interharmonic active power of each oscillation source in the oscillation source identification results, the multiple oscillation sources in the oscillation source identification results are sorted by contribution to obtain the contribution ranking results.

6. The method according to claim 5, characterized in that, The oscillation source identification is performed on the multiple key nodes based on the sign of the interharmonic active power of each key node, and the oscillation source identification results are obtained, including: Key nodes whose interharmonic active power signs correspond to positive signs are identified as oscillation source nodes; All key nodes identified as oscillation source nodes, along with the interharmonic active power corresponding to those oscillation source nodes, are used as the oscillation source identification results.

7. The method according to claim 1, characterized in that, The process of obtaining the interharmonic active power of each key node based on one or more dominant interharmonic frequencies of each key node in the target oscillation mode includes: Based on each of the key nodes, the voltage phasor and current phasor for each dominant interharmonic frequency are obtained; Based on the voltage phasor and current phasor of each dominant interharmonic frequency, the interharmonic active power of the key node is obtained.

8. The method according to claim 7, characterized in that, The process of obtaining the voltage phasor and current phasor for each dominant interharmonic frequency based on each of the key nodes includes: Obtain the complex voltage spectrum and complex current spectrum for each dominant interharmonic frequency; Based on the complex voltage spectrum and complex current spectrum of each dominant interharmonic frequency, the voltage amplitude, voltage phase, current amplitude, and current phase of the dominant interharmonic frequency are calculated respectively. The voltage amplitude and voltage phase, as well as the current amplitude and current phase of each dominant interharmonic frequency, are respectively used as the voltage phasor and current phasor of each dominant interharmonic frequency.

9. The method according to claim 7, characterized in that, The process of obtaining the interharmonic active power of the key node based on the voltage phasor and current phasor of each dominant interharmonic frequency includes: Based on the voltage phasor and current phasor of each dominant interharmonic frequency, the available interharmonic active power of each dominant interharmonic frequency is calculated. Based on each key node, the available interharmonic active power of all dominant interharmonic frequencies of the key node is summed to obtain the interharmonic active power of the key node.

10. The method according to claim 1, characterized in that, After locating the oscillation source of the power system based on the oscillation source identification result and the contribution ranking result, the method further includes: The key nodes corresponding to the located oscillation sources are used as monitoring points; Obtain the time-series variation sequence of the intermediate harmonic active power at the monitoring point; Based on the time-series variation sequence of the intermediate harmonic active power at the monitoring points, the trend of the interharmonic active power is obtained. Based on the interharmonic active power trend, it is determined whether there is a risk of increased oscillation at the monitoring point; When there is a risk of increased oscillation at the monitoring point, an oscillation warning signal is issued to provide early warning to the monitoring point.

11. The method according to claim 10, characterized in that, The determination of whether there is a risk of increased oscillation at the monitoring point based on the interharmonic active power trend includes: When the trend of the interharmonic active power is greater than a preset threshold, or when the interharmonic active power shows an increasing trend within a preset time period, it is determined that there is a risk of increased oscillation at the monitoring point.

12. A power system oscillation source location system, characterized in that, The system includes: The dominant interharmonic frequency determination unit is used to identify one or more dominant interharmonic frequencies of each key node in the target oscillation mode by using fast Fourier transform spectrum analysis based on the instantaneous waveform data of three-phase voltage and three-phase current of multiple key nodes in the power system. An interharmonic active power determination unit is used to obtain the interharmonic active power of each key node based on one or more dominant interharmonic frequencies of each key node in the target oscillation mode. The identification and ranking result determination unit is used to identify and rank the oscillation sources of the power system based on the interharmonic active power of each key node, and obtain the oscillation source identification result and the contribution ranking result. A power system oscillation source localization unit is used to locate the oscillation source of the power system based on the oscillation source identification result and the contribution ranking result.

13. An electronic device, characterized in that, include: One or more processors; The processor is used to store one or more programs; When the one or more programs are executed by the one or more processors, a method for locating power system oscillation sources as described in any one of claims 1 to 11 is implemented.

14. A computer-readable storage medium, characterized in that, It contains a computer program, which, when executed, implements a method for locating a power system oscillation source as described in any one of claims 1 to 11.

Images