CN115986768A - Method and device for identifying dominant oscillation mode of power system - Google Patents

Method and device for identifying dominant oscillation mode of power system Download PDF

Info

Publication number
CN115986768A
CN115986768A CN202211611779.1A CN202211611779A CN115986768A CN 115986768 A CN115986768 A CN 115986768A CN 202211611779 A CN202211611779 A CN 202211611779A CN 115986768 A CN115986768 A CN 115986768A
Authority
CN
China
Prior art keywords
power system
oscillation mode
dominant
determinant
transfer function
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211611779.1A
Other languages
Chinese (zh)
Inventor
雷傲宇
谢小荣
梅勇
汪林光
高永强
许诘翊
黄方能
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
China Southern Power Grid Co Ltd
Original Assignee
Tsinghua University
China Southern Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsinghua University, China Southern Power Grid Co Ltd filed Critical Tsinghua University
Priority to CN202211611779.1A priority Critical patent/CN115986768A/en
Publication of CN115986768A publication Critical patent/CN115986768A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The application relates to the technical field of power system stability analysis, in particular to a method and a device for identifying a dominant oscillation mode of a power system, wherein the method comprises the following steps: establishing an impedance model of each power device in the power system, converting the power system into a multi-input multi-output system based on the impedance model of each power device, calculating a determinant frequency characteristic curve of a transfer function matrix of the multi-input multi-output system, identifying a leading oscillation mode of the power system according to the curve, and quantitatively analyzing the oscillation stability of the power system to obtain a final identification result. According to the method and the device, the leading oscillation mode of the power system is identified according to the curve, so that the oscillation stability of the power system is quantitatively analyzed, the calculated amount is small, and the quantitative analysis of the system oscillation stability of large-scale equipment with the black and gray box can be realized.

Description

Method and device for identifying dominant oscillation mode of power system
Technical Field
The present disclosure relates to the field of power system stability analysis technologies, and in particular, to a method and an apparatus for identifying a dominant oscillation mode of a power system.
Background
With the access of a large number of power electronic equipment in the current power grid, the problem of negative damping electromagnetic oscillation caused by the interaction between the power electronic equipment and the power grid is increasingly prominent, and the safe and stable operation of the power grid is seriously threatened.
In the related art, the existing power system oscillation stability analysis method commonly includes a characteristic value method and an impedance method. The eigenvalue method firstly establishes a differential algebraic equation system of the system, then converts the system into a state space equation form, obtains each oscillation mode of the system by calculating matrix eigenvalues, and can accurately and quantitatively analyze the oscillation stability. The existing impedance method is mainly based on a (generalized) Nyquist method to analyze stability, impedance models of an equipment side and a system side are established, then a system open-loop transfer function is calculated, finally, stability judgment is achieved through a (generalized) Nyquist criterion, and the calculation is simple and easy to use.
However, in the related art, it is difficult for the eigenvalue method to model a large number of "black and gray box" models existing in the system, and in the actual engineering, the number of power electronic devices of the system is large and the topology is complex, resulting in a high dynamic order of the whole system and a large calculation amount, while the impedance method is often only used for qualitative judgment and stability, and it is difficult to accurately obtain the quantitative analysis results of damping, and for a system with a high dimension, the Nyquist criterion is not intuitive enough, and the two methods cannot satisfy the conditions of easy modeling and low calculation complexity while accurately and quantitatively analyzing the oscillation stability, resulting in a small calculation amount of the method, and cannot realize the quantitative analysis of the oscillation stability of the system with the "black and gray box" devices on a large scale, and thus a solution is urgently needed.
Disclosure of Invention
The application provides a method and a device for identifying a dominant oscillation mode of a power system, which are used for solving the technical problems that in the related technology, a characteristic value method and an impedance method cannot accurately and quantitatively analyze oscillation stability, and simultaneously meet the conditions of easiness in modeling and low calculation complexity, so that the method has small calculation amount, and the quantitative analysis of the system oscillation stability of large-scale equipment with a black and grey box cannot be realized.
An embodiment of a first aspect of the present application provides a method for identifying a dominant oscillation mode of a power system, including the following steps: establishing an impedance model of each power device of the power system; converting the power system into a multi-input multi-output system based on the impedance model of each power device; calculating a determinant frequency characteristic curve of a transfer function matrix of the multi-input multi-output system; identifying a dominant oscillation mode of the power system according to the determinant frequency characteristic curve of the transfer function matrix; and quantitatively analyzing the oscillation stability of the power system according to the dominant oscillation mode to obtain a final identification result.
Optionally, in an embodiment of the present application, the establishing an impedance model of each power device of the power system includes: and constructing an impedance model of each power device by using an external characteristic identification method of preset injection disturbance.
Optionally, in an embodiment of the present application, the converting the power system into a multiple-input multiple-output system based on the impedance model of each power device includes: determining input information of voltage of each node of the power system and output information of output current of each node; and constructing the multi-input multi-output system based on the input information, the output information, the impedance models of the electric power devices and the Kirchoff KCL and KVL equations.
Optionally, in an embodiment of the present application, the identifying a dominant oscillation mode of the power system according to the determinant frequency characteristic of the transfer function matrix includes: generating a determinant frequency characteristic curve graph of a transfer function matrix according to the determinant frequency characteristic curve of the transfer function matrix; and searching zero-crossing points of the two curves in the determinant frequency characteristic curve of the transfer function matrix, wherein the dominant oscillation mode exists if an equivalent reactance curve zero-crossing point and an equivalent resistance curve zero-crossing point with the frequency reaching a preset condition exist.
Optionally, in an embodiment of the present application, the quantitatively analyzing the oscillation stability of the power system according to the dominant oscillation mode to obtain a final identification result includes: if a dominant oscillation mode with a real part larger than 0 exists, the dominant oscillation mode is an unstable oscillation mode, the electric power system is judged to be unstable, the oscillation frequency is an imaginary part of the dominant oscillation mode, and the divergence rate is a real part of the dominant oscillation mode; and if the dominant oscillation mode with the real part larger than 0 does not exist, judging that the oscillation of the power system is stable.
An embodiment of a second aspect of the present application provides an apparatus for identifying a dominant oscillation mode of a power system, including: the system comprises a building module, a control module and a control module, wherein the building module is used for building an impedance model of each power device of the power system; the conversion module is used for converting the power system into a multi-input multi-output system based on the impedance model of each power device; the calculation module is used for calculating a determinant frequency characteristic curve of a transfer function matrix of the multi-input multi-output system; the identification module is used for identifying a dominant oscillation mode of the power system according to the determinant frequency characteristic curve of the transfer function matrix; and the analysis module is used for quantitatively analyzing the oscillation stability of the power system according to the dominant oscillation mode to obtain a final identification result.
Optionally, in an embodiment of the present application, the building module includes: the first construction unit is used for constructing the impedance model of each power device by using an external characteristic identification method of preset injection disturbance.
Optionally, in an embodiment of the present application, the conversion module includes: a determination unit for determining input information of each node voltage and output information of each node output current of the power system; a second construction unit configured to construct the multiple-input multiple-output system based on the input information, the output information, the impedance models of the respective power devices, and kirchhoff KCL and KVL equations.
Optionally, in an embodiment of the present application, the identification module includes: the generating unit is used for generating a determinant frequency characteristic curve graph of a transfer function matrix according to the determinant frequency characteristic curve of the transfer function matrix; and the inspection unit is used for searching zero-crossing points of the two curves in the determinant frequency characteristic curve of the transfer function matrix, wherein the dominant oscillation mode exists if an equivalent reactance curve zero-crossing point and an equivalent resistance curve zero-crossing point with the frequency reaching a preset condition exist.
Optionally, in an embodiment of the present application, the analysis module includes: a first determination unit, configured to determine that the power system is unstable, an oscillation frequency is an imaginary part of a dominant oscillation mode, and a divergence rate is a real part of the dominant oscillation mode when the dominant oscillation mode with a real part greater than 0 exists; and the second determination unit is used for determining that the power system oscillation is stable when the dominant oscillation mode with the real part larger than 0 does not exist.
An embodiment of a third aspect of the present application provides an electronic device, including: a memory, a processor and a computer program stored on the memory and executable on the processor, the processor executing the program to implement the power system dominant oscillation pattern recognition method as described in the above embodiments.
An embodiment of a fourth aspect of the present application provides a computer-readable storage medium storing a computer program which, when executed by a processor, implements the power system dominant oscillation pattern recognition method as above.
According to the method and the device, the impedance model of each power device of the power system can be established, the power system is converted into the multi-input multi-output system based on the impedance model of each power device, the determinant frequency characteristic curve of the transfer function matrix of the multi-input multi-output system is calculated, the leading oscillation mode of the power system is identified according to the curve, the oscillation stability of the power system is quantitatively analyzed, and the final identification result is obtained. Therefore, the problems that in the related technology, a characteristic value method and an impedance method cannot accurately and quantitatively analyze the oscillation stability, and simultaneously meet the conditions of easiness in modeling and low calculation complexity, so that the method is small in calculation amount, and the system oscillation stability quantitative analysis of large-scale equipment with a black and gray box cannot be realized are solved.
Additional aspects and advantages of the present application will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the present application.
Drawings
The foregoing and/or additional aspects and advantages of the present application will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:
fig. 1 is a flowchart of a method for identifying a dominant oscillation mode of a power system according to an embodiment of the present application;
FIG. 2 is a schematic diagram of a dominant oscillation mode identification device of a power system according to an embodiment of the present disclosure;
fig. 3 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Detailed Description
Reference will now be made in detail to embodiments of the present application, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the drawings are exemplary and intended to be used for explaining the present application and should not be construed as limiting the present application.
The following describes a method and an apparatus for identifying dominant oscillation mode of a power system according to an embodiment of the present application with reference to the drawings. In order to solve the problems that in the related technologies mentioned in the background art, a characteristic value method and an impedance method cannot accurately and quantitatively analyze oscillation stability, and simultaneously, the conditions of easiness in modeling and low calculation complexity are met, so that the method is small in calculation amount, and quantitative analysis of system oscillation stability of large-scale equipment including a black and gray box cannot be realized, the method for identifying the dominant oscillation mode of the power system is provided. Therefore, the problems that in the related technology, a characteristic value method and an impedance method cannot accurately and quantitatively analyze the oscillation stability, and simultaneously meet the conditions of easiness in modeling and low calculation complexity, so that the method is small in calculation amount, and the system oscillation stability quantitative analysis of large-scale equipment with a black and gray box cannot be realized are solved.
Specifically, fig. 1 is a schematic flowchart of a method for identifying a dominant oscillation mode of a power system according to an embodiment of the present disclosure.
As shown in fig. 1, the method for identifying a dominant oscillation mode of a power system includes the following steps:
in step S101, impedance models of the respective electric devices of the electric power system are established.
It can be understood that, in the embodiment of the present application, the impedance model of each electrical device of the electrical power system may be established through the following steps, so as to obtain the impedance model of the relevant electrical device, and further obtain the relevant information about the impedance influence of the electrical device in the electrical power system, thereby performing the following steps.
Optionally, in an embodiment of the present application, establishing an impedance model of each power device of the power system includes: and constructing an impedance model of each power device by using an external characteristic identification method of preset injection disturbance.
According to the method and the device for identifying the external characteristics of the injected disturbance, the impedance model of each power device can be constructed by using the method for identifying the external characteristics of the preset injected disturbance, and in the actual execution process, the method for identifying the external characteristics of the injected disturbance can be replaced by other impedance modeling methods, for example, the impedance model of each power device is constructed by using a state space equation modeling method.
It should be noted that the method for identifying the external characteristic of the preset injection disturbance is set by a person skilled in the art according to actual situations, and is not limited herein.
By adopting the external characteristic impedance identification, the frequency domain impedance model of the equipment can be obtained under the condition that the internal topology and parameters of the equipment are unknown, so that the method and the device can be suitable for an actual system containing the equipment with the black and gray box, and the application level is wider.
In step S102, the power system is converted into a multiple-input multiple-output system based on the impedance model of each power device.
It can be understood that, in the embodiment of the present application, the power system may be converted into a Multiple input Multiple Output system MIMO (Multiple input Multiple Output) system based on the impedance model of each power device, and the Multiple input Multiple Output system converted by the power system is obtained by using the established impedance model of each power device of the power system, so that the analysis process of the situation that the number of electronic devices in the power system is large and the topology is complex can be simplified, and further calculation is facilitated.
Optionally, in an embodiment of the present application, converting the power system into a multiple-input multiple-output system based on an impedance model of each power device includes: determining input information of voltage of each node of a power system and output information of output current of each node; and constructing a multi-input multi-output system based on the input information, the output information, the impedance model of each power device and the kirchhoff KCL and KVL equations.
As a possible implementation manner, in the embodiment of the present application, it may be determined that the input of the system is the voltage of each node, the output is the output current of each node, and system establishment is performed according to each device impedance model and kirchhoff KCL and KVL equations, so as to convert the power system into a multiple-input multiple-output system
I(ω)=(ω)U(ω)
Wherein, I (omega) is the injection current of each node, U (omega) is the voltage of each node, and Y (omega) is the function matrix of the constructed multi-input multi-output system.
According to the embodiment of the application, the multi-input multi-output system can be constructed by determining the input information of the voltage of each node of the power system and the output information of the output current of each node, based on the input information, the output information, the impedance model of each power device and the Kirchoff KCL and KVL equations, the complexity degree of system construction is reduced, and the applicability in practical engineering is enhanced.
In step S103, a determinant frequency characteristic curve of a transfer function matrix of the mimo system is calculated.
It can be understood that, in the embodiments of the present application, a determinant frequency characteristic curve of a transfer function matrix of a MIMO system can be calculated, and MIMO stability of the MIMO system can be determined by a characteristic equation thereof according to a stability criterion of MIMO of the MIMO system
Y(s)=0
Judging that Y(s) is a transfer function matrix of the multi-input multi-output system, solving the characteristic equation into each oscillation mode of the system, and when the oscillation mode with the real part larger than 0 exists, the oscillation mode is an unstable oscillation mode, and the system generates an oscillation phenomenon.
In the actual implementation process, since a manufacturer does not provide detailed equipment parameters, the impedance model of the equipment can only be obtained by an external characteristic identification method of injection disturbance, so that the transfer function matrix Y(s) of the multi-input multi-output system in the s-domain is difficult to establish, and only the transfer function matrix Y (omega) in the frequency domain can be obtained, therefore, the formula for establishing the determinant frequency characteristic of the transfer function matrix is
D(ω)=etY(ω)
D (omega) is the determinant frequency characteristic of the transfer function matrix, Y (omega) is the transfer function matrix in the frequency domain, and then a graph of the determinant frequency characteristic of the transfer function matrix is drawn, so that the method can be applied to electronic systems with unknown parameters and topological structures of electronic equipment, and the practical application range is widened.
In step S104, a dominant oscillation mode of the power system is identified according to the determinant frequency characteristic curve of the transfer function matrix.
It can be understood that the embodiment of the present application may identify the dominant oscillation mode of the power system according to a transfer function matrix determinant frequency characteristic curve, where in a plotted transfer function matrix determinant frequency characteristic curve, a real part Re (D (ω)) of the transfer function matrix determinant frequency characteristic D (ω) is an equivalent resistance, and an imaginary part Im (D (ω)) of the transfer function matrix determinant frequency characteristic D (ω) is an equivalent reactance.
The dominant oscillation mode is identified and obtained based on the frequency characteristic curve, so that the process of solving equation solutions of traditional characteristic value analysis and s-domain analysis is avoided, the consumed calculation amount is small, and the method and the device are suitable for large-scale system analysis.
Optionally, in an embodiment of the present application, identifying the dominant oscillation mode of the power system according to the determinant frequency characteristic curve of the transfer function matrix includes: generating a determinant frequency characteristic curve graph of a transfer function matrix according to the determinant frequency characteristic curve of the transfer function matrix; and searching zero crossing points of the two curves in a determinant frequency characteristic curve of the transfer function matrix, wherein if an equivalent reactance curve zero crossing point and an equivalent resistance curve zero crossing point with the frequency reaching a preset condition exist, a dominant oscillation mode exists.
It can be understood that the equivalent reactance curve when the frequency reaches the preset condition in the embodiment of the present application may be an equivalent reactance curve with a similar frequency, and if there are an equivalent reactance curve zero crossing point and an equivalent resistance curve zero crossing point with a similar frequency in the determinant frequency characteristic curve of the transfer function matrix, there is a dominant oscillation mode
Figure BDA0003999652740000061
Wherein, ω is R Frequency, ω, at zero crossing of the equivalent resistance Re (D (ω)) curve X Is the frequency, k, at the zero crossing of the equivalent reactance Im (D (ω)) curve R Is the slope, k, of the equivalent resistance Re (D (ω)) curve at the zero crossing point X J is an imaginary unit, which is the slope at the zero crossing of the equivalent reactance Im (D (ω)) curve.
The preset conditions are set by those skilled in the art according to actual conditions, and are not specifically limited herein.
According to the method and the device, a determinant frequency characteristic curve of the transfer function matrix can be generated according to the determinant frequency characteristic curve of the transfer function matrix, and the zero crossing points of two curves are searched in the determinant frequency characteristic curve of the transfer function matrix, wherein if the zero crossing point of an equivalent reactance curve and the zero crossing point of an equivalent resistance curve with the frequency reaching a preset condition exist, a dominant oscillation mode exists, the dominant oscillation mode is judged through the zero crossing point of the frequency characteristic curve, the complexity of a calculation process is reduced, and the execution capacity of the oscillation stability analysis of a large-scale equipment system is improved.
In step S105, the oscillation stability of the power system is quantitatively analyzed according to the dominant oscillation mode, and a final identification result is obtained.
It can be understood that, in the embodiment of the present application, each dominant oscillation mode of the system in the full frequency band is obtained through identification, and the result is obtained by performing quantitative analysis on the oscillation stability of the power system according to the identified dominant oscillation mode.
According to the method and the device, the oscillation stability of the power system can be quantitatively analyzed according to the leading oscillation mode, and the final identification result is obtained, so that the oscillation stability degree of the power system is judged, and the safe and stable operation of a power grid is further guaranteed.
Optionally, in an embodiment of the present application, the quantitatively analyzing the oscillation stability of the power system according to the dominant oscillation mode to obtain a final identification result includes: if the dominant oscillation mode with the real part larger than 0 exists, the dominant oscillation mode is an unstable oscillation mode, the power system is judged to be unstable, the oscillation frequency is the imaginary part of the dominant oscillation mode, and the divergence rate is the real part of the dominant oscillation mode; and if the dominant oscillation mode with the real part larger than 0 does not exist, judging that the oscillation of the power system is stable.
It can be understood that, in the embodiment of the present application, if there is a dominant oscillation mode whose real part is greater than 0, the dominant oscillation mode is an unstable oscillation mode, and thus it is determined that the power system is unstable, and the oscillation frequency is an imaginary part of the dominant oscillation mode, and the divergence rate is a real part of the dominant oscillation mode, and if there is no dominant oscillation mode whose real part is greater than 0, it is determined that the power system is stable in oscillation, and the criterion is more intuitive and accurate by accurately and quantitatively analyzing the oscillation stability of the power system.
According to the method for identifying the dominant oscillation mode of the power system, the power system can be converted into a multi-input multi-output system by establishing impedance models of all power equipment of the power system and based on the impedance models of all the power equipment, a determinant frequency characteristic curve of a transfer function matrix of the multi-input multi-output system is calculated, the dominant oscillation mode of the power system is identified according to the curve, the oscillation stability of the power system is quantitatively analyzed, and a final identification result is obtained. Therefore, the problems that in the related technology, a characteristic value method and an impedance method cannot accurately and quantitatively analyze the oscillation stability, and simultaneously meet the conditions of easiness in modeling and low calculation complexity, so that the method is small in calculation amount, and the system oscillation stability quantitative analysis of large-scale equipment with a black and gray box cannot be realized are solved.
Next, a dominant oscillation pattern recognition apparatus of a power system according to an embodiment of the present application is described with reference to the drawings.
Fig. 2 is a block diagram of a dominant oscillation mode identification device of a power system according to an embodiment of the present application.
As shown in fig. 2, the apparatus 10 for identifying dominant oscillation mode of power system includes: a construction module 100, a transformation module 200, a calculation module 300, a recognition module 400, and an analysis module 500.
The building module 100 is configured to build an impedance model of each power device of the power system.
The conversion module 200 is configured to convert the power system into a multiple-input multiple-output system based on the impedance model of each power device.
The calculation module 300 is configured to calculate a determinant frequency characteristic curve of a transfer function matrix of a mimo system.
An identification module 400 for identifying a dominant oscillation mode of the power system according to the determinant frequency characteristic curve of the transfer function matrix.
The analysis module 500 is configured to quantitatively analyze the oscillation stability of the power system according to the dominant oscillation mode to obtain a final identification result.
Optionally, in an embodiment of the present application, the building module 100 includes: a first building element.
The first construction unit is used for constructing an impedance model of each power device by using an external characteristic identification method of preset injection disturbance.
Optionally, in an embodiment of the present application, the conversion module 200 includes: a validation unit and a second building unit.
The confirming unit is used for confirming input information of voltage of each node of the power system and output information of output current of each node.
And the second construction unit is used for constructing the multi-input multi-output system based on the input information, the output information, the impedance model of each power device and the Kirchoff KCL and KVL equations.
Optionally, in an embodiment of the present application, the recognition module 400 includes: generating unit and checking unit
The generating unit is used for generating a determinant frequency characteristic curve graph of the transfer function matrix according to the determinant frequency characteristic curve of the transfer function matrix.
And the inspection unit is used for searching zero crossing points of the two curves in a determinant frequency characteristic curve of the transfer function matrix, wherein if an equivalent reactance curve zero crossing point and an equivalent resistance curve zero crossing point with the frequency reaching a preset condition exist, a dominant oscillation mode exists.
Optionally, in an embodiment of the present application, the analysis module 500 includes: a first determination unit and a second determination unit.
The first determination unit is used for determining that the dominant oscillation mode is an unstable oscillation mode when the dominant oscillation mode with the real part larger than 0 exists, the oscillation frequency is the imaginary part of the dominant oscillation mode, and the divergence rate is the real part of the dominant oscillation mode.
And the second determination unit is used for determining that the power system oscillation is stable when the dominant oscillation mode with the real part larger than 0 does not exist.
It should be noted that the above explanation of the embodiment of the method for identifying a dominant oscillation mode of an electric power system is also applicable to the device for identifying a dominant oscillation mode of an electric power system of the embodiment, and is not repeated herein.
According to the dominant oscillation mode identification device for the power system, the impedance model of each power device of the power system is established, the power system is converted into the multi-input multi-output system based on the impedance model of each power device, the determinant frequency characteristic curve of the transfer function matrix of the multi-input multi-output system is calculated, the dominant oscillation mode of the power system is identified according to the curve, the oscillation stability of the power system is quantitatively analyzed, and a final identification result is obtained. Therefore, the problems that in the related technology, a characteristic value method and an impedance method cannot accurately and quantitatively analyze the oscillation stability, and simultaneously meet the conditions of easiness in modeling and low calculation complexity, so that the method is small in calculation amount, and the system oscillation stability quantitative analysis of large-scale equipment with a black and gray box cannot be realized are solved.
Fig. 3 is a schematic structural diagram of an electronic device according to an embodiment of the present application. The electronic device may include:
a memory 301, a processor 302, and a computer program stored on the memory 301 and executable on the processor 302.
The processor 302, when executing the program, implements the method for identifying dominant oscillation mode of power system provided in the above embodiments.
Further, the electronic device further includes:
a communication interface 303 for communication between the memory 301 and the processor 302.
A memory 301 for storing computer programs executable on the processor 302.
The memory 301 may comprise high-speed RAM memory, and may also include non-volatile memory, such as at least one disk memory.
If the memory 301, the processor 302 and the communication interface 303 are implemented independently, the communication interface 303, the memory 301 and the processor 302 may be connected to each other through a bus and perform communication with each other. The bus may be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, an Extended ISA (EISA) bus, or the like. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of illustration, only one thick line is shown in FIG. 3, but this does not mean only one bus or one type of bus.
Alternatively, in practical implementation, if the memory 301, the processor 302 and the communication interface 303 are integrated on one chip, the memory 301, the processor 302 and the communication interface 303 may complete communication with each other through an internal interface.
The processor 302 may be a Central Processing Unit (CPU), an Application Specific Integrated Circuit (ASIC), or one or more Integrated circuits configured to implement embodiments of the present Application.
The present embodiment also provides a computer-readable storage medium on which a computer program is stored, which when executed by a processor implements the power system dominant oscillation pattern recognition method as above.
In the description of the present specification, reference to the description of "one embodiment," "some embodiments," "an example," "a specific example," or "some examples" or the like means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present application. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or N embodiments or examples. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction.
Furthermore, the terms "first", "second" and "first" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of the feature. In the description of the present application, "N" means at least two, e.g., two, three, etc., unless specifically limited otherwise.
Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or N executable instructions for implementing steps of a custom logic function or process, and alternate implementations are included within the scope of the preferred embodiment of the present application in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the embodiments of the present application.
The logic and/or steps represented in the flowcharts or otherwise described herein, such as an ordered listing of executable instructions that can be considered to implement logical functions, can be embodied in any computer-readable medium for use by or in connection with an instruction execution system, apparatus, or device, such as a computer-based system, processor-containing system, or other system that can fetch the instructions from the instruction execution system, apparatus, or device and execute the instructions. For the purposes of this description, a "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of the computer-readable medium would include the following: an electrical connection (electronic device) having one or N wires, a portable computer diskette (magnetic device), a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). Further, the computer readable medium could even be paper or another suitable medium upon which the program is printed, as the program can be electronically captured, via optical scanning of the paper or other medium, then compiled, interpreted or otherwise processed in a suitable manner if necessary, and then stored in a computer memory.
It should be understood that portions of the present application may be implemented in hardware, software, firmware, or a combination thereof. In the above embodiments, the N steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. If implemented in hardware, as in another embodiment, any one or combination of the following techniques, which are known in the art, may be used: a discrete logic circuit having a logic gate circuit for implementing a logic function on a data signal, an application specific integrated circuit having an appropriate combinational logic gate circuit, a Programmable Gate Array (PGA), a Field Programmable Gate Array (FPGA), or the like.
It will be understood by those skilled in the art that all or part of the steps carried out in the method of implementing the above embodiments may be implemented by hardware related to instructions of a program, which may be stored in a computer readable storage medium, and the program, when executed, includes one or a combination of the steps of the method embodiments.
In addition, functional units in the embodiments of the present application may be integrated into one processing module, or each unit may exist alone physically, or two or more units are integrated into one module. The integrated module can be realized in a hardware mode, and can also be realized in a software functional module mode. The integrated module, if implemented in the form of a software functional module and sold or used as a stand-alone product, may also be stored in a computer readable storage medium.
The storage medium mentioned above may be a read-only memory, a magnetic or optical disk, etc. Although embodiments of the present application have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present application, and that variations, modifications, substitutions and alterations may be made to the above embodiments by those of ordinary skill in the art within the scope of the present application.

Claims (12)

1. A method for identifying a dominant oscillation mode of a power system is characterized by comprising the following steps:
establishing an impedance model of each power device of the power system;
converting the power system into a multi-input multi-output system based on the impedance model of each power device;
calculating a determinant frequency characteristic curve of a transfer function matrix of the multi-input multi-output system;
identifying a dominant oscillation mode of the power system according to the determinant frequency characteristic curve of the transfer function matrix; and
and quantitatively analyzing the oscillation stability of the power system according to the dominant oscillation mode to obtain a final identification result.
2. The method of claim 1, wherein the establishing an impedance model of each power device in the power system comprises:
and constructing an impedance model of each power device by using an external characteristic identification method of preset injection disturbance.
3. The method of claim 1, wherein transforming the power system into a multiple-input multiple-output system based on the impedance models of the respective power devices comprises:
determining input information of voltage of each node of the power system and output information of output current of each node;
and constructing the multi-input multi-output system based on the input information, the output information, the impedance models of the electric power devices and the Kirchoff KCL and KVL equations.
4. The method of claim 1, wherein identifying the dominant oscillatory pattern of the power system according to the transfer function matrix determinant frequency characteristic curve comprises:
generating a determinant frequency characteristic curve graph of a transfer function matrix according to the determinant frequency characteristic curve of the transfer function matrix;
and searching zero-crossing points of the two curves in the determinant frequency characteristic curve of the transfer function matrix, wherein the dominant oscillation mode exists if an equivalent reactance curve zero-crossing point and an equivalent resistance curve zero-crossing point with the frequency reaching a preset condition exist.
5. The method of claim 1, wherein the quantitatively analyzing the oscillation stability of the power system according to the dominant oscillation mode to obtain a final identification result comprises:
if a dominant oscillation mode with a real part larger than 0 exists, the dominant oscillation mode is an unstable oscillation mode, the electric power system is judged to be unstable, the oscillation frequency is an imaginary part of the dominant oscillation mode, and the divergence rate is a real part of the dominant oscillation mode;
and if the dominant oscillation mode with the real part larger than 0 does not exist, judging that the oscillation of the power system is stable.
6. An apparatus for identifying dominant oscillation modes of a power system, comprising:
the system comprises a building module, a control module and a control module, wherein the building module is used for building an impedance model of each power device of the power system;
the conversion module is used for converting the power system into a multi-input multi-output system based on the impedance model of each power device;
the calculation module is used for calculating a determinant frequency characteristic curve of a transfer function matrix of the multi-input multi-output system;
the identification module is used for identifying a dominant oscillation mode of the power system according to the determinant frequency characteristic curve of the transfer function matrix; and
and the analysis module is used for quantitatively analyzing the oscillation stability of the power system according to the dominant oscillation mode to obtain a final identification result.
7. The apparatus of claim 6, wherein the building module comprises:
the first construction unit is used for constructing the impedance model of each power device by using an external characteristic identification method of preset injection disturbance.
8. The apparatus of claim 6, wherein the conversion module comprises:
a confirmation unit for determining input information of each node voltage and output information of each node output current of the power system;
a second construction unit configured to construct the multiple-input multiple-output system based on the input information, the output information, the impedance models of the respective power devices, and kirchhoff KCL and KVL equations.
9. The apparatus of claim 6, wherein the recognition module comprises:
the generating unit is used for generating a determinant frequency characteristic curve graph of a transfer function matrix according to the determinant frequency characteristic curve of the transfer function matrix;
and the inspection unit is used for searching zero-crossing points of the two curves in the determinant frequency characteristic curve of the transfer function matrix, wherein the dominant oscillation mode exists if an equivalent reactance curve zero-crossing point and an equivalent resistance curve zero-crossing point with the frequency reaching a preset condition exist.
10. The apparatus of claim 6, wherein the analysis module comprises:
a first determination unit, configured to determine that, when a dominant oscillation mode with a real part greater than 0 exists, the dominant oscillation mode is an unstable oscillation mode, and that the power system is unstable, the oscillation frequency is an imaginary part of the dominant oscillation mode, and the divergence rate is a real part of the dominant oscillation mode;
and the second determination unit is used for determining that the power system oscillation is stable when the dominant oscillation mode with the real part larger than 0 does not exist.
11. An electronic device, comprising: memory, a processor and a computer program stored on the memory and executable on the processor, the processor executing the program to implement the power system dominant oscillation pattern recognition method as claimed in any one of claims 1 to 5.
12. A computer-readable storage medium, on which a computer program is stored, the program being executable by a processor for implementing a method for dominant oscillation pattern recognition of a power system as claimed in any one of claims 1 to 5.
CN202211611779.1A 2022-12-14 2022-12-14 Method and device for identifying dominant oscillation mode of power system Pending CN115986768A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211611779.1A CN115986768A (en) 2022-12-14 2022-12-14 Method and device for identifying dominant oscillation mode of power system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211611779.1A CN115986768A (en) 2022-12-14 2022-12-14 Method and device for identifying dominant oscillation mode of power system

Publications (1)

Publication Number Publication Date
CN115986768A true CN115986768A (en) 2023-04-18

Family

ID=85969257

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211611779.1A Pending CN115986768A (en) 2022-12-14 2022-12-14 Method and device for identifying dominant oscillation mode of power system

Country Status (1)

Country Link
CN (1) CN115986768A (en)

Similar Documents

Publication Publication Date Title
US9964980B2 (en) Method and apparatus for optimal power flow with voltage stability for large-scale electric power systems
CN113037589B (en) Pressure testing method and device of gateway equipment, testing platform and storage medium
EP3861615A1 (en) Distribution grid admittance estimation with limited nonsynchronized measurements
CN112330488A (en) Power grid frequency situation prediction method based on transfer learning
CN110532731B (en) Rapid calculation method for voltage sag concave domain
US20150331062A1 (en) Failure Detection Method and Detection Device for Inverter
Colbert et al. Using trajectory measurements to estimate the region of attraction of nonlinear systems
Ramezani et al. Multiple objective functions for falsification of cyber-physical systems
Ma et al. Two-stage bayesian sequential change diagnosis
Tostado-Véliz et al. A Novel Power Flow Solution Paradigm for Well and Ill-Conditioned Cases
CN113011055B (en) Motor high-temperature fault simulation method and system based on finite element technology
CN115986768A (en) Method and device for identifying dominant oscillation mode of power system
WO2009120190A1 (en) Stable equilibrium point (sep) calculation apparatus of power system
Faraji-Niri et al. Stochastic stability and stabilization of semi-Markov jump linear systems with uncertain transition rates
CN115864502A (en) Method and device for analyzing stability and designing stability of small signal of alternating current and direct current power system
US8326589B2 (en) Stable equilibrium point (SEP) calculation apparatus of power system
CN113933718B (en) Retired battery capacity sorting method, retired battery capacity sorting device, retired battery capacity sorting equipment and storage medium
CN108964282A (en) Dispositions method, device and the computer-readable medium of equipment for monitoring power quality
CN114784800A (en) Method and system for analyzing resonance stability of power distribution network
CN111045407A (en) Method and device for controlling production of corrosion foil
CN109436980A (en) The condition detection method and system of elevator components
KR101330526B1 (en) Method, Apparatus and Recording Medium for Migration of Prediction Systems
Theiler et al. Kalman Filter Tuning Using Multi-Objective Genetic Algorithm for State and Parameter Estimation of Lithium-Ion Cells. Batteries 2022, 8, 104
CN114552599A (en) Oscillation stability analysis method based on polymerization impedance determinant logarithmic derivative calculation
Colbert et al. Estimating the region of attraction using stable trajectory measurements

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination