Oscillation source localization method, device, equipment and medium for multiple new energy power stations area

Abstract

The present invention discloses a method and device for locating oscillation sources in a multi-new energy power station area, which are used to solve the technical problem of low accuracy in locating oscillation sources in a multi-new energy power station area. The present invention includes: establishing a state space model of the multi-new energy power station area; the multi-new energy power station area includes a plurality of new energy power stations; based on the state space model, calculating the relative amplitude and predicted flow direction of the oscillation component in the power transmission network corresponding to the multi-new energy power station area; according to the relative amplitude and the predicted flow direction, selecting a new energy power station where the oscillation component meets a preset condition as the first new energy power station; calculating the participation factor of each new energy power station in the oscillation mode; determining the second new energy power station according to the participation factor; calculating the intersection of the first new energy power station and the second new energy power station, and locating the oscillation source of the multi-new energy power station area according to the intersection.

Description

Technical Field

[0001] The present invention relates to the technical field of oscillation source location, and particularly to a method, device, equipment and medium for locating an oscillation source in a multi-new energy power station area. Background Art

[0002] In recent years, the problem of broadband oscillation in new energy power stations has been frequent, seriously affecting the safe and stable operation of the power system and the development of new energy. In a certain area, a serious sub / supersynchronous oscillation occurred in the sending system of the conventional DC transmission with multiple wind farms and thermal power bundled. This oscillation was caused by the interaction between the wind farms and the weak AC grid to which they were connected, and further propagated to the vicinity of the thermal power units, triggering shaft torsional vibration and resulting in a significant reduction in the transmission power. At the same time, due to the large number of wind farms in the accident area, it was difficult for on-site operators to determine the dominant wind farm in the accident. After the accident, the operators cut off all the power stations in the area to suppress the oscillation in the system, further leading to a reduction in the outgoing power and affecting the safe and stable operation of the receiving and sending end systems. Therefore, for new energy-intensive areas such as wind power and photovoltaic, an effective oscillation source location method is needed to assist operators in determining the power stations to be cut off during an accident and to assist in the analysis and rectification measure research after the accident.

[0003] At present, for oscillation source location, on the one hand, it is considered to filter and identify the oscillation waveform extracted online or identify the transient energy flow, which requires the aid of artificial intelligence algorithms. The overall scheme is relatively complex and the accuracy of actual application is relatively low. On the other hand, it is considered to judge whether a certain power station / equipment emits or absorbs oscillation by the current direction. However, in the multi-new energy power station area, it is possible that a single new energy power station is the dominant oscillation source, causing the common oscillation of the remaining power stations, or there are two new energy power stations that are both the dominant oscillation sources, but show the mutual oscillation characteristics that one power station emits oscillation and the other power station absorbs oscillation. Therefore, it is impossible to simply locate the power station / equipment that emits oscillation as the dominant oscillation source by the oscillation current direction. Summary of the Invention

[0004] The present invention provides a method and device for locating an oscillation source in a multi-new energy power station area to solve the technical problem of low accuracy in locating an oscillation source in a multi-new energy power station area.

[0005] The present invention provides a method for locating an oscillation source in a multi-new energy power station area, including:

[0006] Establish a state space model of the multi-new energy power station area; the multi-new energy power station area includes a plurality of new energy power stations;

[0007] Based on the state space model, calculate the relative amplitude and predicted flow direction of the oscillation component in the transmission network corresponding to the multi-new energy power station area;

[0008] Select a new energy power station whose oscillation component meets the preset conditions as the first new energy power station according to the relative amplitude and the predicted flow direction;

[0009] Calculate the participation factors of each new energy power station in the oscillation mode;

[0010] Determine the second new energy power station according to the participation factors;

[0011] Calculate the intersection of the first new energy power station and the second new energy power station, and locate the oscillation source of the multi-new energy power station area according to the intersection.

[0012] Optionally, the step of establishing the state space model of the multi-new energy power station area includes:

[0013] Collect the state variables of each device in the multi-new energy power station area;

[0014] Use the device and the corresponding state variables to construct the mathematical model of each device;

[0015] Generate the interface equation between the coordinate systems of each device according to the preset rotating coordinate system;

[0016] Use the mathematical models of each device and the interface equation between the coordinate systems to generate the state space model of the multi-new energy power station area.

[0017] Optionally, the step of calculating the relative amplitude and the predicted flow direction of the oscillation component in the transmission network corresponding to the multi-new energy power station area based on the state space model includes:

[0018] Obtain the preset coefficient matrix, and calculate the eigenvalues of the state space model through the coefficient matrix;

[0019] Determine the eigenvalues with a real part not less than 0 as the unstable mode eigenvalues;

[0020] Obtain the current state variables of each line of the transmission network;

[0021] Calculate the right eigenvector of the line corresponding to the current state variable and the unstable mode eigenvalue;

[0022] Calculate the amplitude and phase of each right eigenvector of the line;

[0023] Take the amplitude of the right eigenvector of the line with the largest amplitude as the reference amplitude, and take the phase of the right eigenvector of the line with the largest amplitude as the reference phase, and normalize the amplitude and phase of each right eigenvector of the line to obtain the normalized amplitude and the normalized phase;

[0024] Determine the relative amplitude of the oscillation component in each line of the transmission network according to the normalized amplitude;

[0025] Determine the predicted flow direction of the oscillation component in the power transmission network according to the normalized phase.

[0026] Optionally, the step of calculating the participation factor of each new energy power station in the oscillation mode includes:

[0027] Calculate the state variables of each new energy power station and the right eigenvector of the power station corresponding to the eigenvalue of the instability mode.

[0028] Calculate the state variables of each new energy power station and the left eigenvector of the power station corresponding to the eigenvalue of the instability mode.

[0029] Calculate the participation factor of each new energy power station in the oscillation mode by using the right eigenvector and the left eigenvector of the power station.

[0030] The present invention also provides an oscillation source positioning device for a multi-new energy power station area, including:

[0031] A state space model establishment module for establishing a state space model of a multi-new energy power station area; the multi-new energy power station area includes a plurality of new energy power stations;

[0032] A relative amplitude and predicted flow direction calculation module for calculating the relative amplitude and predicted flow direction of the oscillation component in the power transmission network corresponding to the multi-new energy power station area based on the state space model;

[0033] A first new energy power station determination module for selecting a new energy power station whose oscillation component meets a preset condition as the first new energy power station according to the relative amplitude and the predicted flow direction;

[0034] A participation factor calculation module for calculating the participation factor of each new energy power station in the oscillation mode;

[0035] A second new energy power station determination module for determining the second new energy power station according to the participation factor;

[0036] An oscillation source positioning module for calculating the intersection of the first new energy power station and the second new energy power station, and positioning the oscillation source of the multi-new energy power station area according to the intersection.

[0037] Optionally, the state space model establishment module includes:

[0038] An equipment parameter collection sub-module for collecting the equipment parameters of each equipment in the multi-new energy power station area, the equipment includes the new energy power stations, traditional units, DC power transmission systems and each connection line of the power transmission network, and the equipment parameters include state variables;

[0039] A mathematical model construction sub-module, which is used to construct the mathematical models of each device by using the device and the corresponding state variables;

[0040] An interface equation generation sub-module between coordinate systems, which is used to generate the interface equations between the coordinate systems of each device according to a preset rotating coordinate system;

[0041] A state space model generation sub-module, which is used to generate the state space model of the multi-new energy power station area by using the mathematical models of each device and the interface equations between the coordinate systems;

[0042] Optionally, the relative amplitude and predicted flow direction calculation module includes:

[0043] An eigenvalue calculation sub-module, which is used to obtain a preset coefficient matrix and calculate the eigenvalues of the state space model through the coefficient matrix;

[0044] A judgment sub-module, which is used to judge whether the real parts of all the eigenvalues are less than 0;

[0045] An unstable mode eigenvalue determination sub-module, which is used to, if not, determine the eigenvalues with real parts not less than 0 as the unstable mode eigenvalues;

[0046] A current state variable acquisition sub-module, which is used to acquire the current state variables of each line of the power transmission network;

[0047] A line right eigenvector calculation sub-module, which is used to calculate the line right eigenvectors corresponding to the current state variables and the unstable mode eigenvalues;

[0048] An amplitude and phase calculation sub-module, which is used to calculate the amplitudes and phases of each of the line right eigenvectors;

[0049] A normalization sub-module, which is used to use the amplitude of the line right eigenvector with the largest amplitude as the reference amplitude and the phase of the line right eigenvector with the largest amplitude as the reference phase to normalize the amplitudes and phases of each of the line right eigenvectors to obtain the normalized amplitude and the normalized phase;

[0050] A relative amplitude determination sub-module, which is used to determine the relative amplitudes of the oscillation components in each line of the power transmission network according to the normalized amplitudes;

[0051] A predicted flow direction determination sub-module, which is used to determine the predicted flow direction of the oscillation components in the power transmission network according to the normalized phases;

[0052] Optionally, the participation factor calculation module includes:

[0053] A station right eigenvector calculation sub-module, which is used to calculate the station right eigenvectors corresponding to the state variables of each new energy power station and the unstable mode eigenvalues;

[0054] The sub-module for calculating the left eigenvector of the power station is used to calculate the left eigenvector of the power station corresponding to the instability mode eigenvalue for each state variable of the new energy power station;

[0055] The participation factor calculation sub-module is used to calculate the participation factor of each new energy power station in the oscillation mode by using the right eigenvector of the power station and the left eigenvector of the power station.

[0056] The present invention also provides an electronic device, which includes a processor and a memory:

[0057] The memory is used to store program codes and transmit the program codes to the processor;

[0058] The processor is used to execute the oscillation source localization method for the multi-new energy power station area as described in any one of the above according to the instructions in the program codes.

[0059] The present invention also provides a computer-readable storage medium, which is used to store program codes, and the program codes are used to execute the oscillation source localization method for the multi-new energy power station area as described in any one of the above.

[0060] It can be seen from the above technical solutions that the present invention has the following advantages: The present invention establishes a state space model for a multi-new energy power station area; the multi-new energy power station area includes multiple new energy power stations; based on the state space model, the relative amplitude and predicted flow direction of the oscillation component in the corresponding transmission network of the multi-new energy power station area are calculated; according to the relative amplitude and predicted flow direction, the new energy power station whose oscillation component meets the preset conditions is selected as the first new energy power station; the participation factor of each new energy power station in the oscillation mode is calculated; the second new energy power station is determined according to the participation factor; the intersection of the first new energy power station and the second new energy power station is calculated, and the oscillation source of the multi-new energy power station area is located according to the intersection. Thereby improving the accuracy of oscillation source localization. Description of the Drawings

[0061] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained according to these drawings without creative efforts.

[0062] Figure 1 It is a step flow chart of an oscillation source localization method for a multi-new energy power station area provided by an embodiment of the present invention;

[0063] Figure 2The structural block diagram of a multi-new energy power station area oscillation source positioning device provided by an embodiment of the present invention. Detailed implementation manners

[0064] An embodiment of the present invention provides a method and device for positioning an oscillation source in a multi-new energy power station area, which is used to solve the technical problem of low accuracy in positioning the oscillation source in the multi-new energy power station area.

[0065] In order to make the invention purpose, features, and advantages of the present invention more obvious and understandable, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the embodiments described below are only a part of the embodiments of the present invention, rather than all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts belong to the scope of protection of the present invention.

[0066] Please refer to Figure 1 , Figure 1 The step flow chart of a method for positioning an oscillation source in a multi-new energy power station area provided by an embodiment of the present invention.

[0067] A method for positioning an oscillation source in a multi-new energy power station area provided by the present invention may specifically include the following steps:

[0068] Step 101, establish a state space model of the multi-new energy power station area; the multi-new energy power station area includes multiple new energy power stations;

[0069] A new energy power station refers to the collection of all equipment below the grid connection point of a new energy field (wind farm, solar power station, etc.) that is centrally connected to the power system, including transformers, busbars, converters, energy storage, wind turbines, photovoltaic power generation equipment, power regulation equipment, and auxiliary equipment, etc.

[0070] In the embodiment of the present invention, by establishing a state space model of the multi-new energy power station area, data support for subsequent oscillation source positioning can be obtained.

[0071] In specific implementation, the steps of establishing a state space model of the multi-new energy power station area may include the following sub-steps:

[0072] S11, collect the state variables of each device in the multi-new energy power station area;

[0073] S12, use the device and the corresponding state variables to construct a mathematical model of each device;

[0074] S13, generate an interface equation between coordinate systems of each device according to a preset rotation coordinate system;

[0075] S14. Generate the state space model of the multi-new energy power station area by using the mathematical models of each device and the interface equations between coordinate systems.

[0076] In the specific implementation, first, the device parameters of new energy power stations, traditional units, DC transmission systems, and each connection line of the transmission network in the multi-new energy power station area can be collected. The device parameters of each device are reduced to a unified base power, and according to the corresponding device parameters, the mathematical models of each new energy power station, traditional unit, DC transmission system, and each connection line of the transmission network are established. The formulas are as follows:

[0077]

[0078] Among them, ΔX i is the state variable of each device X i (such as new energy power stations, traditional units, DC transmission systems, and each section of the connection lines of the transmission network), Δi i is the device port current variable, Δu i is the device port voltage variable, A ii , B ii , C ii , D ii are the state space coefficient matrices.

[0079] Then, a unified rotating coordinate system is selected to convert the port voltage and current components of each new energy power station, traditional unit, DC transmission system, and each section of the connection line to the unified rotating coordinate system, and the interface equations between different coordinate systems are completed. The interface equations between coordinate systems are as follows:

[0080]

[0081]

[0082] Among them, u xi / u yi is the x / y axis component of the port voltage of each device (such as new energy power stations, traditional units, DC transmission systems, and each section of the connection lines of the transmission network), i xi / i yi is the x / y axis component of the port current of each device, u di / u qi is the d / q axis component of the port voltage of each device, i di / i qi is the d / q axis component of the port current of each device, δ i is the angle between the dq coordinate system controlled by each device and the unified xy rotating coordinate system.

[0083] Next, combining the mathematical models of new energy power stations, traditional units, DC transmission systems, and transmission networks within the multi-new energy power station area and the interface equations between different coordinate systems, a state space model of the multi-new energy power station area is established, and the formula is as follows:

[0084]

[0085] Among them, A is the coefficient matrix of the state space model, ΔX = [X 1 , X 2 , X 3 , …, X n is the state variable of the equipment, and the current state variable of the line current is ΔX L = [X L1 , X L2 , X L3 , …, X Lm .

[0086] Step 102: Based on the state space model, calculate the relative amplitude and predicted flow direction of the oscillation component in the transmission network corresponding to the multi-new energy power station area;

[0087] After constructing the state space model of the multi-new energy power station, the relative amplitude and predicted flow direction of the oscillation component in the transmission network corresponding to the multi-new energy power station area can be calculated based on this state space model.

[0088] In one example, the steps of calculating the relative amplitude and predicted flow direction of the oscillation component in the transmission network corresponding to the multi-new energy power station area based on the state space model may specifically include the following sub-steps:

[0089] S21: Obtain the preset coefficient matrix and calculate the eigenvalues of the state space model through the coefficient matrix;

[0090] S22: Determine the eigenvalues with a real part not less than 0 as the unstable mode eigenvalues;

[0091] S23: Obtain the current state variables of each line in the transmission network;

[0092] S24: Calculate the right eigenvector of the line corresponding to the current state variable and the unstable mode eigenvalue;

[0093] S25: Calculate the amplitude and phase of each line's right eigenvector;

[0094] S26: Using the amplitude of the right eigenvector of the line with the largest amplitude as the reference amplitude and the phase of the right eigenvector of the line with the largest amplitude as the reference phase, normalize the amplitude and phase of each line's right eigenvector to obtain the normalized amplitude and normalized phase;

[0095] S27. Determine the relative amplitude of the oscillation component in each line of the transmission network according to the normalized amplitude;

[0096] S28. Determine the predicted flow direction of the oscillation component in the transmission network according to the normalized phase.

[0097] In specific implementation, first, the coefficient matrix in the state - space model can be obtained to calculate the eigenvalue λ in the state - space model:

[0098] |λI - A| = 0

[0099] where I is the identity matrix.

[0100] Then, judge whether the real part real(λ i ) of all the eigenvalues in the state - space model is less than zero. If so, the system has no oscillation problem and the calculation can be ended; if not, the system has an oscillation problem, and the unstable - mode eigenvalue λ i with real(λ i )≥0 can be calculated through the following formula for the right - eigenvector U i :

[0101] AU i =λ i U i

[0102] where the current state variable ΔX L =[X L1 ,X L2 ,X L3 ,…,X Lm of the transmission - network line with respect to the unstable - mode eigenvalue λ i has the line right - eigenvector U Li =[U Li1 ,U Li2 ,…,U Lim . T .

[0103] The amplitude M Li and phase P i of each element in U i can be calculated through the following formula:

[0104]

[0105]

[0106] Next, select the amplitude of the element with the largest amplitude in U Li as the reference amplitude of 1.0 pu and the phase as the reference phase of 0°, and normalize the amplitudes and phases of the remaining elements:

[0107] P' j = P j - P M

[0108] where M M is the amplitude of the maximum amplitude element, P M is the phase of the maximum amplitude element, M j and P j are the amplitudes and phases of the remaining elements respectively, M' j and P' j are the amplitudes and phases of the remaining elements after normalization respectively.

[0109] If the normalization result of each element is greater than 0.35 pu, it is determined that the corresponding line contains obvious oscillation components. If the normalization result of the phase of each element is less than 60°, it is determined that the oscillation components in the corresponding line are in phase with the line with the maximum amplitude oscillation. Thus, the relative amplitude and predicted flow direction of the oscillation components in the transmission network corresponding to the multi-new energy power station area are obtained.

[0110] Step 103, select the new energy power stations whose oscillation components meet the preset conditions as the first new energy power stations according to the relative amplitude and predicted flow direction;

[0111] After determining the relative amplitude and predicted flow direction of the oscillation components in the transmission network corresponding to the multi-new energy power station area, the new energy power stations with larger oscillation components in the line can be marked as the first new energy power stations.

[0112] Step 104, calculate the participation factors of each new energy power station in the oscillation mode;

[0113] In the embodiment of the present invention, the step of calculating the participation factors of each new energy power station in the oscillation mode may include the following sub-steps:

[0114] S41, calculate the right eigenvector of the station corresponding to the state variable of each new energy power station and the instability mode eigenvalue;

[0115] S42, calculate the left eigenvector of the station corresponding to the state variable of each new energy power station and the instability mode eigenvalue;

[0116] S43, calculate the participation factors of each new energy power station in the oscillation mode by using the right eigenvector and left eigenvector of the station.

[0117] In the embodiment of the present invention, based on the state space model, the state variable ΔX of each energy power station can be selected WF = [X WF1 , X WF2 , X WF3 , …, X WFz about the instability mode eigenvalue λ iThe right eigenvector U of the power station WFi = [U WF1 , U WF2 , U WF3 , …, U WFz T , and by calculating the state variable ΔX of each new energy power station WF = [X WF1 , X WF2 , X WF3 , …, X WFz , the left eigenvector V of the power station with respect to the eigenvalue λ i of the instability mode WFi = [V WF1 , V WF2 , V WF3 , …, V WFz T .

[0118]

[0119] Next, the participation factor PF of each new energy power station in the oscillation mode is calculated using the right eigenvector and left eigenvector of the power station WFi :

[0120] PF WFi = U WF1 V WF1 + U WF2 V WF2 + … + U WFz V WFz

[0121] Step 105, determine the second new energy power station according to the participation factor;

[0122] Take the amplitude of the participation factor with the largest amplitude among the participation factors of all new energy power stations as the reference amplitude, normalize the assignment of the participation factor of each new energy power station, and regard the new energy power station with the participation factor amplitude greater than 0.5 pu as the second new energy power station.

[0123] Step 106, calculate the intersection of the first new energy power station and the second new energy power station, and locate the oscillation source of the multi-new energy power station area according to the intersection.

[0124] Take the intersection of the first new energy power station and the second new energy power station and determine it as the oscillation source of the multi-new energy power station area.

[0125] ​​The present invention improves the accuracy of oscillation source localization by establishing a state space model for a multi-new energy power station area; the multi-new energy power station area includes multiple new energy power stations; based on the state space model, calculating the relative amplitude and predicted flow direction of the oscillation component in the transmission network corresponding to the multi-new energy power station area; selecting, according to the relative amplitude and predicted flow direction, the new energy power station whose oscillation component meets the preset conditions as the first new energy power station; calculating the participation factor of each new energy power station in the oscillation mode; determining the second new energy power station according to the participation factor; calculating the intersection of the first new energy power station and the second new energy power station, and locating the oscillation source of the multi-new energy power station area according to the intersection.

[0126] Please refer to Figure 2 , Figure 2 which is the structural block diagram of an oscillation source localization device for a multi-new energy power station area provided by an embodiment of the present invention.

[0127] An embodiment of the present invention provides an oscillation source localization device for a multi-new energy power station area, including:

[0128] A state space model establishment module 201, configured to establish a state space model for a multi-new energy power station area; the multi-new energy power station area includes multiple new energy power stations;

[0129] A relative amplitude and predicted flow direction calculation module 202, configured to calculate the relative amplitude and predicted flow direction of the oscillation component in the transmission network corresponding to the multi-new energy power station area based on the state space model;

[0130] A first new energy power station determination module 203, configured to select, according to the relative amplitude and predicted flow direction, the new energy power station whose oscillation component meets the preset conditions as the first new energy power station;

[0131] A participation factor calculation module 204, configured to calculate the participation factor of each new energy power station in the oscillation mode;

[0132] A second new energy power station determination module 205, configured to determine the second new energy power station according to the participation factor;

[0133] An oscillation source localization module 206, configured to calculate the intersection of the first new energy power station and the second new energy power station, and locate the oscillation source of the multi-new energy power station area according to the intersection.

[0134] In an embodiment of the present invention, the state space model establishment module 201 includes:

[0135] An equipment parameter collection sub-module, configured to collect the state variables of each equipment in the multi-new energy power station area;

[0136] A mathematical model construction sub-module, configured to construct a mathematical model of each equipment by using the equipment and the corresponding state variables;

[0137] The interface equation generation sub-module between coordinate systems is used to generate the interface equations between the coordinate systems of each device according to a preset rotated coordinate system;

[0138] The state space model generation sub-module is used to generate the state space model of the multi-new energy power station area by using the mathematical models of each device and the interface equations between the coordinate systems.

[0139] In the embodiment of the present invention, the relative amplitude and predicted flow direction calculation module 202 includes:

[0140] The eigenvalue calculation sub-module is used to obtain a preset coefficient matrix and calculate the eigenvalues of the state space model through the coefficient matrix;

[0141] The unstable mode eigenvalue determination sub-module is used to determine the eigenvalues with a real part not less than 0 as the unstable mode eigenvalues;

[0142] The current state variable acquisition sub-module is used to acquire the current state variables of each line of the transmission network;

[0143] The line right eigenvector calculation sub-module is used to calculate the line right eigenvectors corresponding to the current state variables and the unstable mode eigenvalues;

[0144] The amplitude and phase calculation sub-module is used to calculate the amplitude and phase of each line right eigenvector;

[0145] The normalization sub-module is used to normalize the amplitude and phase of each line right eigenvector with the amplitude of the line right eigenvector with the largest amplitude as the reference amplitude and the phase of the line right eigenvector with the largest amplitude as the reference phase to obtain the normalized amplitude and the normalized phase;

[0146] The relative amplitude determination sub-module is used to determine the relative amplitude of the oscillation component in each line of the transmission network according to the normalized amplitude;

[0147] The predicted flow direction determination sub-module is used to determine the predicted flow direction of the oscillation component in the transmission network according to the normalized phase.

[0148] In the embodiment of the present invention, the participation factor calculation module 204 includes:

[0149] The power station right eigenvector calculation sub-module is used to calculate the power station right eigenvectors corresponding to the state variables of each new energy power station and the unstable mode eigenvalues;

[0150] The power station left eigenvector calculation sub-module is used to calculate the power station left eigenvectors corresponding to the state variables of each new energy power station and the unstable mode eigenvalues;

[0151] The participation factor calculation sub-module is used to calculate the participation factors of each new energy power station in the oscillation mode by using the right eigenvector and left eigenvector of the power station.

[0152] An embodiment of the present invention also provides an electronic device, which includes a processor and a memory:

[0153] The memory is used to store program codes and transmit the program codes to the processor;

[0154] The processor is used to execute the oscillation source localization method for the multi-new energy power station area according to the instructions in the program codes.

[0155] An embodiment of the present invention also provides a computer-readable storage medium, which is used to store program codes, and the program codes are used to execute the oscillation source localization method for the multi-new energy power station area.

[0156] Those skilled in the art can clearly understand that for the convenience and brevity of description, the specific working processes of the above-described systems, devices, and units can refer to the corresponding processes in the foregoing method embodiments, and will not be described herein again.

[0157] Each embodiment in this specification is described in a progressive manner. Each embodiment focuses on the differences from other embodiments. The same or similar parts among the various embodiments can be referred to each other.

[0158] Those skilled in the art should understand that the embodiments of the present invention can be provided as a method, a device, or a computer program product. Therefore, the embodiments of the present invention can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the embodiments of the present invention can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.

[0159] The embodiments of the present invention are described with reference to the flowcharts and / or block diagrams of the methods, terminal devices (systems), and computer program products according to the embodiments of the present invention. It should be understood that each process and / or block in the flowcharts and / or block diagrams, and the combination of processes and / or blocks in the flowcharts and / or block diagrams can be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing terminal devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing terminal devices generate for implementing in the process Figure 1 each process or multiple processes and / or blocks Figure 1apparatus for the functions specified in one or more blocks.

[0160] These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing terminal device to work in a specific manner, such that the instructions stored in the computer-readable memory produce a manufacture including an instruction device that implements the functions specified in one Figure 1 one or more processes and / or blocks Figure 1 one or more blocks.

[0161] These computer program instructions may also be loaded onto a computer or other programmable data processing terminal device, such that a series of operational steps are executed on the computer or other programmable terminal device to produce a computer-implemented process, and thus the instructions executed on the computer or other programmable terminal device provide steps for implementing the functions specified in one Figure 1 one or more processes and / or blocks Figure 1 one or more blocks.

[0162] Although the preferred embodiments of the embodiments of the present invention have been described, those skilled in the art can make additional changes and modifications once they learn the basic creative concept. Therefore, the appended claims are intended to be construed to include the preferred embodiments as well as all changes and modifications falling within the scope of the embodiments of the present invention.

[0163] Finally, it should also be noted that in this text, relational terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the term "comprising", "including" or any other variant thereof is intended to cover non-exclusive inclusion, such that a process, method, article or terminal device comprising a series of elements not only includes those elements, but also includes other elements not expressly listed, or also includes elements inherent to such process, method, article or terminal device. Without further limitation, an element defined by the statement "comprising an..." does not exclude the presence of additional identical elements in the process, method, article or terminal device comprising the element.

[0164] As described above, the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements on some of the technical features; and these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for locating oscillation sources in a multi-new energy power station area, characterized in that, it includes: Establish a state space model of the multi-new energy power station area; The multi-new energy power station area includes multiple new energy power stations; Based on the state space model, calculate the relative amplitude and predicted flow direction of the oscillation component in the power transmission network corresponding to the multi-new energy power station area; According to the relative amplitude and the predicted flow direction, select the new energy power station where the oscillation component meets the preset conditions as the first new energy power station; Calculate the participation factor of each new energy power station in the oscillation mode; Determine the second new energy power station according to the participation factor; Calculate the intersection of the first new energy power station and the second new energy power station, and locate the oscillation source of the multi-new energy power station area according to the intersection; Among them, the step of calculating the relative amplitude and predicted flow direction of the oscillation component in the power transmission network corresponding to the multi-new energy power station area based on the state space model includes: Obtain a preset coefficient matrix, and calculate the eigenvalues of the state space model through the coefficient matrix; Determine the unstable mode eigenvalues with the real part not less than 0; Obtain the current state variables of each line in the power transmission network; Calculate the right eigenvector of the line corresponding to the current state variable and the unstable mode eigenvalue; Calculate the amplitude and phase of each line right eigenvector; Taking the amplitude of the line right eigenvector with the largest amplitude as the reference amplitude and the phase of the line right eigenvector with the largest amplitude as the reference phase, normalize the amplitude and phase of each line right eigenvector to obtain the normalized amplitude and normalized phase; Determine the relative amplitude of the oscillation component in each line of the power transmission network according to the normalized amplitude; Determine the predicted flow direction of the oscillation component in the power transmission network according to the normalized phase.

2. The method according to claim 1, characterized in that, The step of establishing the state space model of the multi-new energy power station area includes: Collect the state variables of each device in the multi-new energy power station area; Use the device and the corresponding state variables to construct the mathematical model of each device; Generate the interface equation between coordinate systems of each device according to a preset rotation coordinate system; Use the mathematical model of each device and the interface equation between coordinate systems to generate the state space model of the multi-new energy power station area.

3. The method according to claim 1, characterized in that, The step of calculating the participation factor of each new energy power station in the oscillation mode includes: Calculate the right eigenvector of the power station corresponding to the state variable of each new energy power station and the unstable mode eigenvalue; Calculate the left eigenvector of the power station corresponding to the state variable of each new energy power station and the unstable mode eigenvalue; Calculate the participation factor of each new energy power station in the oscillation mode by using the right eigenvector of the power station and the left eigenvector of the power station.

4. An oscillation source location device for a multi-new energy power station area, characterized in that, it includes: A state space model establishment module for establishing a state space model of a multi-new energy power station area; The multi-new energy power station area includes multiple new energy power stations; The relative amplitude and predicted flow calculation module is used to calculate the relative amplitude and predicted flow of the oscillation component in the transmission network corresponding to the multi-new energy power station area based on the state space model; The first new energy power station determination module is used to select the new energy power station whose oscillation component meets the preset conditions as the first new energy power station according to the relative amplitude and the predicted flow; The participation factor calculation module is used to calculate the participation factors of each new energy power station in the oscillation mode; The second new energy power station determination module is used to determine the second new energy power station according to the participation factor; The oscillation source positioning module is used to calculate the intersection of the first new energy power station and the second new energy power station, and locate the oscillation source in the multi-new energy power station area according to the intersection; Among them, the relative amplitude and predicted flow calculation module includes: The eigenvalue calculation sub-module is used to obtain a preset coefficient matrix and calculate the eigenvalues of the state space model through the coefficient matrix; The unstable mode eigenvalue determination sub-module is used to determine the eigenvalues with a real part not less than 0 as the unstable mode eigenvalues; The current state variable acquisition sub-module is used to acquire the current state variables of each line of the transmission network; The line right eigenvector calculation sub-module is used to calculate the line right eigenvector corresponding to the current state variable and the unstable mode eigenvalue; The amplitude and phase calculation sub-module is used to calculate the amplitude and phase of each line right eigenvector; The normalization sub-module is used to normalize the amplitude and phase of each line right eigenvector with the amplitude of the line right eigenvector with the largest amplitude as the reference amplitude and the phase of the line right eigenvector with the largest amplitude as the reference phase to obtain the normalized amplitude and the normalized phase; The relative amplitude determination sub-module is used to determine the relative amplitude of the oscillation component in each line of the transmission network according to the normalized amplitude; The predicted flow determination sub-module is used to determine the predicted flow of the oscillation component in the transmission network according to the normalized phase.

5. The device according to claim 4, wherein, The state space model establishment module includes: The device parameter collection sub-module is used to collect the state variables of each device in the multi-new energy power station area; The mathematical model construction sub-module is used to construct the mathematical models of each device by using the device and the corresponding state variables; The coordinate system interface equation generation sub-module is used to generate the coordinate system interface equations of each device according to a preset rotation coordinate system; The state space model generation sub-module is used to generate the state space model of the multi-new energy power station area by using the mathematical models of each device and the coordinate system interface equations.

6. The device according to claim 4, wherein, The participation factor calculation module includes: The station right eigenvector calculation sub-module is used to calculate the station right eigenvector corresponding to the state variable of each new energy power station and the unstable mode eigenvalue; The station left eigenvector calculation sub-module is used to calculate the station left eigenvector corresponding to the state variable of each new energy power station and the unstable mode eigenvalue; The participation factor calculation sub-module is used to calculate the participation factors of each new energy power station in the oscillation mode by using the right eigenvector and the left eigenvector of the power station.

7. An electronic device, characterized in that, the device includes a processor and a memory: the memory is used to store program codes and transmit the program codes to the processor; the processor is used to execute the oscillation source location method for the multi-new energy power station area according to the instructions in the program codes of any one of claims 1-3.

8. A computer-readable storage medium, characterized in that, the computer-readable storage medium is used to store program codes, and the program codes are used to execute the oscillation source location method for the multi-new energy power station area according to any one of claims 1-3.

Images