CN110569577A - Frequency domain mode method-based sub-synchronous oscillation tracing method and system for new energy system - Google Patents
Frequency domain mode method-based sub-synchronous oscillation tracing method and system for new energy system Download PDFInfo
- Publication number
- CN110569577A CN110569577A CN201910789949.7A CN201910789949A CN110569577A CN 110569577 A CN110569577 A CN 110569577A CN 201910789949 A CN201910789949 A CN 201910789949A CN 110569577 A CN110569577 A CN 110569577A
- Authority
- CN
- China
- Prior art keywords
- oscillation
- matrix
- mode
- frequency domain
- subsynchronous
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000010355 oscillation Effects 0.000 title claims abstract description 206
- 238000000034 method Methods 0.000 title claims abstract description 63
- 239000011159 matrix material Substances 0.000 claims abstract description 101
- 230000035945 sensitivity Effects 0.000 claims abstract description 46
- 238000004458 analytical method Methods 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000013016 damping Methods 0.000 claims description 11
- 239000013598 vector Substances 0.000 description 7
- 239000000463 material Substances 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
the invention discloses a frequency domain mode method-based sub-synchronous oscillation tracing method and system for a new energy system, wherein the method comprises the following steps: establishing an impedance network model of a target system, forming a node admittance matrix and a loop impedance matrix of the network, and obtaining an oscillation mode of the system according to determinant zero points of the matrixes to obtain a subsynchronous oscillation mode; calculating participation factors of the nodes and the branches to the subsynchronous oscillation mode based on frequency domain mode analysis, and positioning an oscillation influence area by comparing the participation factors of the nodes and the branches; and establishing a sensitivity matrix according to the influence area of the subsynchronous oscillation, calculating the sensitivity of each element, and positioning the key elements influencing the oscillation by comparing the sensitivities of the elements. And finally, carrying out oscillation tracing according to the affected area and the key element of the subsynchronous oscillation. The method can position the main influence area of the oscillation and the key elements influencing the oscillation, thereby realizing the identification of the oscillation source of the subsynchronous oscillation.
Description
Technical Field
The invention relates to the technical field of power system stability analysis, in particular to a new energy system subsynchronous oscillation tracing method and system based on a frequency domain mode method.
Background
With the continuous improvement of new energy permeability, the power electronic converter is widely connected to a power system, and the dynamic behavior of the power system is obviously changed. In recent years, many novel oscillation accidents occur at home and abroad. Compared with typical subsynchronous/supersynchronous oscillation accidents caused by variable-current power supplies such as wind power and the like, the subsynchronous/supersynchronous oscillation accidents reduce the quality of electric energy, destroy the safety of equipment and endanger the reliable operation of a power grid. In order to improve the stability of the system, the oscillation source of the novel subsynchronous oscillation needs to be identified, so that important information is provided for the prevention and suppression of the oscillation.
The existing oscillation source identification research is mainly developed aiming at forced oscillation in low-frequency oscillation, and common methods include an energy function method, a simulation method and the like. Because the mechanism of the subsynchronous oscillation and the forced oscillation generated by the new energy system is different, the subsynchronous oscillation and the forced oscillation are the result of the interaction of the power electronic converter and the alternating current and direct current power grid, and a uniform oscillation source definition is not formed yet, the existing method is not suitable for the source tracing of the novel subsynchronous oscillation. Considering that new energy systems are generally complex, detailed modeling of target systems in the time domain is difficult.
disclosure of Invention
The present invention is directed to solving, at least to some extent, one of the technical problems in the related art.
therefore, an object of the present invention is to provide a method for tracing to a source of a subsynchronous oscillation of a new energy system based on a frequency domain mode method, which can locate a main influence area of the oscillation and key elements influencing the oscillation, so as to identify an oscillation source of the subsynchronous oscillation.
The invention further aims to provide a new energy system subsynchronous oscillation traceability system based on the frequency domain mode method.
In order to achieve the above object, the present invention provides a sub-synchronous oscillation tracing method for a new energy system based on a frequency domain mode method, including: establishing an impedance network model of a target system, and forming a node admittance matrix and a loop impedance matrix of a network according to the impedance network model; calculating determinant zeros of the node admittance matrix and the loop impedance matrix to obtain an oscillation mode of a target system so as to position a subsynchronous oscillation mode; calculating participation factors of each node and each branch to the subsynchronous oscillation mode based on frequency domain mode analysis, and positioning an influence area of subsynchronous oscillation by comparing the participation factors of each node and each branch; establishing a sensitivity matrix according to the affected area of the subsynchronous oscillation, calculating the sensitivity of each element, and positioning the key element affecting the oscillation by comparing the sensitivities of the elements; and tracing the subsynchronous oscillation according to the affected area of the subsynchronous oscillation and the key element.
According to the new energy system subsynchronous oscillation tracing method based on the frequency domain mode method, an impedance network model of a target system is established, and an oscillation mode of the system is obtained based on a node admittance matrix and a loop impedance matrix of the network, so that oscillation damping and frequency information of subsynchronous oscillation are obtained; the participation factor indexes of the nodes/branches and the sensitivity indexes of the elements are defined, and the solution is carried out based on a frequency domain mode analysis method; according to the calculation results of the participation factors and the sensitivity, the main influence area of the oscillation and the key elements influencing the oscillation can be positioned, so that the oscillation source of the subsynchronous oscillation is identified.
in addition, the new energy system subsynchronous oscillation tracing method based on the frequency domain mode method according to the embodiment of the invention may further have the following additional technical features:
in an embodiment of the present invention, the establishing an impedance network model of the target system further includes: establishing a frequency domain impedance and admittance model of each element in a target system; and connecting the frequency domain impedance of each element according to the topological structure of the target system to obtain the impedance network model.
in an embodiment of the present invention, the calculating determinant zero of the node admittance matrix and the loop impedance matrix to obtain the oscillation mode of the target system further includes: and calculating the numerator polynomial of the determinant of the node admittance matrix and the loop impedance matrix to obtain a plurality of pairs of complex zeros, wherein each pair of complex zeros corresponds to an oscillation mode.
in one embodiment of the invention, the oscillation mode is represented as:
si＝σi±jωi
Wherein s isiis an oscillation mode, σiFor damping of oscillations, omegaiJ is the unit of an imaginary number for the oscillation angular frequency.
In an embodiment of the present invention, the calculating the participation factor of each node and branch to the subsynchronous oscillation mode based on the frequency domain pattern analysis further includes: calculating participation factors of each node to subsynchronous oscillation according to the node admittance matrix; and calculating the participation factor of each branch to subsynchronous oscillation according to the loop impedance matrix.
In an embodiment of the present invention, the calculating the sensitivity of each element further includes: dividing the elements into series elements and parallel elements according to their positions in the nodal admittance matrix; calculating the sensitivity formula of the series element as follows: sen ═ Sjj-Sjk-Skj+Skkwhere S is the sensitivity matrix, j and k are connection nodes, SjkElement of jth row and kth column of S, Sjjis the jth row and jth column element of S, SkjIs the jth row and jth column element of S, SkkThe kth row and the kth column of S; calculating the sensitivity formula of the parallel element as follows: sen ═ Siiwhere S is the sensitivity matrix, i is a connection node, SiiIs the ith row and ith column element of S.
In order to achieve the above object, another aspect of the present invention provides a sub-synchronous oscillation traceability system of a new energy system based on a frequency domain mode method, including: the system comprises an establishing module, a judging module and a judging module, wherein the establishing module is used for establishing an impedance network model of a target system and forming a node admittance matrix and a loop impedance matrix of a network according to the impedance network model; the calculation module is used for calculating determinant zeros of the node admittance matrix and the loop impedance matrix, and solving an oscillation mode of a target system so as to position a subsynchronous oscillation mode; the first comparison module is used for calculating participation factors of each node and each branch to the subsynchronous oscillation mode based on frequency domain mode analysis, and positioning an influence area of subsynchronous oscillation by comparing the participation factors of each node and each branch; the second comparison module is used for establishing a sensitivity matrix according to the influence area of the subsynchronous oscillation, calculating the sensitivity of each element and positioning the key element influencing the oscillation by comparing the sensitivity of each element; and the tracing module is used for tracing the subsynchronous oscillation according to the affected area of the subsynchronous oscillation and the key element.
according to the new energy system subsynchronous oscillation tracing system based on the frequency domain mode method, an impedance network model of a target system is established, and the oscillation mode of the system is obtained based on a node admittance matrix and a loop impedance matrix of the network, so that the information of oscillation damping and frequency of subsynchronous oscillation is obtained; the participation factor indexes of the nodes/branches and the sensitivity indexes of the elements are defined, and the solution is carried out based on a frequency domain mode analysis method; according to the calculation results of the participation factors and the sensitivity, the main influence area of the oscillation and the key elements influencing the oscillation can be positioned, so that the oscillation source of the subsynchronous oscillation is identified.
in addition, the sub-synchronous oscillation traceability system of the new energy system based on the frequency domain mode method according to the embodiment of the invention may further have the following additional technical features:
Further, in an embodiment of the present invention, the establishing module further includes: the establishing unit is used for establishing frequency domain impedance and admittance models of all elements in the target system; and the connecting unit is used for connecting the frequency domain impedance of each element according to the topological structure of the target system so as to obtain the impedance network model.
further, in an embodiment of the present invention, the calculation module further includes: and calculating the numerator polynomial of the determinant of the node admittance matrix and the loop impedance matrix to obtain a plurality of pairs of complex zeros, wherein each pair of complex zeros corresponds to an oscillation mode, and a subsynchronous oscillation mode is found in the oscillation mode.
Further, in one embodiment of the present invention, the oscillation mode is represented as:
si＝σi±jωi
wherein s isiis an oscillation mode, σifor damping of oscillations, omegaiJ is the unit of an imaginary number for the oscillation angular frequency.
Additional aspects and advantages of the invention 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 invention.
Drawings
The foregoing and/or additional aspects and advantages of the present invention 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 sub-synchronous oscillation tracing method for a new energy system based on a frequency domain mode method according to an embodiment of the present invention;
Fig. 2 is a flowchart of an embodiment of a sub-synchronous oscillation tracing method for a new energy system based on a frequency domain mode method according to the present invention;
Fig. 3 is a schematic structural diagram of a sub-synchronous oscillation traceability system of a new energy system based on a frequency domain mode method according to an embodiment of the present invention.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, 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 illustrative and intended to be illustrative of the invention and are not to be construed as limiting the invention.
The method and system for tracing the source of the sub-synchronous oscillation of the new energy system based on the frequency domain mode method provided by the embodiment of the invention are described below with reference to the accompanying drawings.
fig. 1 is a flowchart of a sub-synchronous oscillation tracing method for a new energy system based on a frequency domain mode method according to an embodiment of the present invention.
As shown in fig. 1, the frequency domain mode method-based sub-synchronous oscillation tracing method for the new energy system includes the following steps:
in step S1, an impedance network model of the target system is established, and a node admittance matrix and a loop impedance matrix of the network are formed according to the impedance network model.
that is, as shown in fig. 2, an impedance network model of the target system is established, forming a node admittance matrix and a loop impedance matrix of the network.
in one embodiment of the present invention, establishing an impedance network model of the target system further comprises: establishing a frequency domain impedance and admittance model of each element in a target system; and connecting the frequency domain impedance of each element according to the topological structure of the target system to obtain an impedance network model.
For example, step S1 in the embodiment of the present invention may specifically be:
101, establishing a frequency domain impedance and an admittance model of each element in the system, which are respectively marked as Z(s) and Y(s), wherein Z(s) and Y(s) are reciprocal, as follows:
Wherein j and k are the order of a polynomial, ai(i ═ 0,1, …, k) and bi(i ═ 0,1, …, j) is the coefficient of the ith order of the polynomial, and s is the laplacian.
And 102, connecting the impedances of the elements according to the topological structure of the system to obtain an impedance network model of the system.
Each element impedance in the network is treated as a branch, connected by a node 103. Each node and branch in the network is numbered. Let n be the number of independent nodes of the network, b be the number of branches, and l be the number of independent loops (l ═ b-n).
104, respectively recording a node admittance matrix and a loop impedance matrix of the network as Y(s) and Z(s), and calculating the formula as follows:
Y(s)＝Ay(s)AT
Z(s)＝Bz(s)BT
A and B are respectively a node-branch incidence matrix and a loop-branch incidence matrix; y(s) and z(s) are respectively a branch admittance matrix and a branch impedance matrix, which are diagonal matrices, and diagonal elements are admittance and impedance of each branch.
In step S2, determinant zeros of the node admittance matrix and the loop impedance matrix are calculated, and the oscillation mode of the target system is obtained to locate the subsynchronous oscillation mode.
That is, as shown in fig. 2, in step S2, a subsynchronous oscillation pattern is found out of the obtained oscillation patterns.
In an embodiment of the present invention, calculating determinant zeros of a node admittance matrix and a loop impedance matrix to obtain an oscillation mode of a target system, further includes: and calculating the numerator polynomial of the determinant of the node admittance matrix and the loop impedance matrix to obtain a plurality of pairs of complex zeros, wherein each pair of complex zeros corresponds to an oscillation mode.
Wherein the oscillation mode is represented as:
si＝σi±jωi
Wherein s isiIs an oscillation mode, σiFor damping of oscillations, omegaij is the unit of an imaginary number for the oscillation angular frequency.
specifically, after step S1, step S2 calculates the process as:
201, calculating a numerator polynomial of determinants of a node admittance matrix Y(s) and a loop impedance matrix Z(s), and marking the numerator polynomial as P(s);
202, zero point of the molecular polynomial P(s) is found. Each pair of complex zeros corresponds to an oscillation mode, denoted si(i ═ 1,2, …), as follows:
si＝σi±jωi
wherein σiFor damping of oscillations, omegaiIs the oscillation angular frequency. Let the oscillation frequency be fi，fi＝ωi/2π。
203, finding out subsynchronous oscillation mode from the obtained oscillation modes, and recording the subsynchronous oscillation mode as sk。
In step S3, based on the frequency domain pattern analysis, the participation factors of each node and branch to the subsynchronous oscillation mode are calculated, and the influence region of the subsynchronous oscillation is located by comparing the participation factors of each node and branch
In an embodiment of the present invention, calculating the participation factor of each node and branch to the sub-synchronous oscillation mode based on the frequency domain pattern analysis further includes:
Calculating participation factors of each node to subsynchronous oscillation according to the node admittance matrix;
and calculating the participation factor of each branch to subsynchronous oscillation according to the loop impedance matrix.
As shown in fig. 2, step S3 is to calculate the participation factor of each node in subsynchronous oscillation based on the node admittance matrix Y (S). The method specifically comprises the following steps:
Substituting subsynchronous oscillation mode into node admittance matrix, i.e. making s equal to skObtaining Y(s)k)；
mixing Y(s)k) Diagonalization to obtain Y(s)k) The eigenvalue matrix and left and right eigenvector matrix of (a), are as follows:
Y(sk)＝LΛR
Λ＝diag(λ1,λ2,…,λn)
L＝[L1,L2,…,Ln]T
R＝[R1,R2,…,Rn]
Where Λ is the eigenvalue matrix and L and R are the left and right eigenvector matrices, respectively. Lambda [ alpha ]tIs Y(s)k) Is determined by the characteristic value of (a),And Rtare each lambdatCorresponding left and right feature vectors (t ═ 1,2, …, n).
In Y(s)k) There is one zero feature root among all the feature roots of (1), called a key modality. Suppose the key modality is λp。
Left eigenvector corresponding to key modalitythe considerable degree of subsynchronous oscillation mode of each node is characterized. Defining node observability vector, denoted as oN，oNthe ith element ofFor the observability index of node # i (i ═ 1,2, …, n), the calculation formula is as follows:
Wherein L istpIs the p-th row and t-th column element of the left eigenvector matrix L (t ═ 1,2, …, n).
Right eigenvector R corresponding to key modepthe controllable degree of each node to the subsynchronous oscillation mode is characterized. Defining a node controllability vector, denoted as cN，cNthe ith element ofFor the controllability index of the node # i (i ═ 1,2, …, n), the calculation formula is as follows:
Wherein R istpIs the t-th row and p-th column element of the right eigenvector matrix R (t ═ 1,2, …, n).
The participation factor of the node is defined as the product of the observability and the controllability of the node, and the calculation formula is as follows:
Wherein the content of the first and second substances,Is the participation factor of node # i (i ═ 1,2, …, n).
then, as shown in fig. 2, step S3 calculates the participation factor of each branch in the subsynchronous oscillation mode according to the loop impedance matrix. The method specifically comprises the following steps:
Substituting subsynchronous oscillation mode into loop impedance matrix, i.e. making s equal to skobtaining Z(s)k)。
A reaction of Z(s)k) Diagonalizing to obtain an eigenvalue matrix and a left and right eigenvector matrix, as follows:
Z(sk)＝VΛW
Λ＝diag(μ1,μ2,…,μl)
V＝[V1,V2,…,Vn]T
W＝[W1,W2,…,Wn]
Where Λ is the eigenvalue matrix and V and W are the left and right eigenvector matrices, respectively. Mu.st(t ═ 1,2, …, l) is Z(s)k) Is determined by the characteristic value of (a),And WtAre respectively mutcorresponding left and right feature vectors.
At Z(s)k) There is one zero feature root among all the feature roots of (1), called a key modality. Suppose the key modality is μq。
Vqthe considerable extent of the subsynchronous oscillation mode of each loop is characterized. Define branch observability vector, denoted as oB，oBThe ith element ofFor the observability index (i ═ 1,2, …, b) of branch # i, the calculation formula is as follows:
V'＝BTV
Wherein B is a loop-branch correlation matrix, Vq'tThe qth row of matrix V' and the tth column element (t ═ 1,2, …, b).
WqThe controllable degree of the subsynchronous oscillation mode of each loop is represented. Defining branch controllability vector, denoted as cB，cBThe ith element offor the controllability index of branch # i (i ═ 1,2, …, b), the calculation formula is as follows:
W＝WB
Wherein, WtqThe t-th row and the q-th column of the matrix W' (t ═ 1,2, …, b).
the participation factor of the branch is defined as the product of the observability and the controllability of the branch, and the calculation formula is as follows:
Wherein the content of the first and second substances,Is the participation factor of branch # i (i ═ 1,2, …, b).
Finally, as shown in fig. 2, the calculated participation factors of each node and branch are compared, and the region formed by the nodes and branches with high participation factors is the main influence region of oscillation.
In step S4, a sensitivity matrix is established based on the influence region of the subsynchronous oscillation, the sensitivity of each element is calculated, and the sensitivity of each element is compared to locate the key element influencing the oscillation.
Specifically, in step 401, a sensitivity matrix is defined, denoted as S, and the calculation formula is as follows:
Wherein L ispAnd RpIs a key mode lambdapCorresponding left and right feature vectors.
in step 402, the sensitivity of the element is calculated based on the sensitivity matrix.
The elements are divided into series elements and parallel elements according to their position in the nodal admittance matrix. The parallel elements are connected by non-reference nodes and reference nodes, the admittance of which only appears at the diagonal elements of the matrix. The series elements are connected by two non-reference nodes, the admittance of which occurs at diagonal and non-diagonal elements of the matrix.
For a parallel element, assuming that its connection node is node # i, its sensitivity calculation formula is:
Sen＝Sii
wherein S isiiIs the ith row and ith column element of S.
For a series element, assuming it is connected by node # j and node # k, its sensitivity calculation formula is:
Sen＝Sjj-Sjk-Skj+Skk
Wherein S isjkIs the jth row and kth column element of S.
In step 403, as shown in fig. 2, the sensitivities of the elements are compared, and the element with high sensitivity is a key element for influencing oscillation.
In step S5, subsynchronous oscillation tracing is performed according to the affected area and the key element of the subsynchronous oscillation.
According to the new energy system subsynchronous oscillation tracing method based on the frequency domain mode method, which is provided by the embodiment of the invention, an impedance network model of a target system is established, and the oscillation mode of the system is obtained based on a node admittance matrix and a loop impedance matrix of the network, so that the information of the oscillation damping and the frequency of subsynchronous oscillation is obtained; the participation factor indexes of the nodes/branches and the sensitivity indexes of the elements are defined, and the solution is carried out based on a frequency domain mode analysis method; according to the calculation results of the participation factors and the sensitivity, the main influence area of the oscillation and the key elements influencing the oscillation can be positioned, so that the oscillation source of the subsynchronous oscillation is identified.
next, a sub-synchronous oscillation traceability system of a new energy system based on a frequency domain mode method according to an embodiment of the present invention is described with reference to the drawings.
Fig. 3 is a subsynchronous oscillation tracing system of a new energy system based on a frequency domain mode method according to an embodiment of the invention.
As shown in fig. 3, the frequency domain mode method-based sub-synchronous oscillation tracing system 10 of the new energy system includes: the system comprises a building module 100, a calculating module 200, a first comparing module 300, a second comparing module 400 and a tracing module 500.
The establishing module 100 is configured to establish an impedance network model of a target system, and form a node admittance matrix and a loop impedance matrix of a network according to the impedance network model.
In one embodiment of the present invention, the establishing module 100 further comprises:
The establishing unit is used for establishing frequency domain impedance and admittance models of all elements in the target system;
And the connecting unit is used for connecting the frequency domain impedance of each element according to the topological structure of the target system so as to obtain an impedance network model.
And the calculation module 200 is configured to calculate determinant zero points of the node admittance matrix and the loop impedance matrix, and obtain an oscillation mode of the target system to locate the subsynchronous oscillation mode.
In one embodiment of the present invention, the calculation module 200 further comprises:
And calculating the numerator polynomial of the determinant of the node admittance matrix and the loop impedance matrix to obtain a plurality of pairs of complex zeros, wherein each pair of complex zeros corresponds to an oscillation mode, and a subsynchronous oscillation mode is found in the oscillation mode.
Wherein the oscillation mode is represented as:
si＝σi±jωi
Wherein s isiis an oscillation mode, σiFor damping of oscillations, omegaij is the unit of an imaginary number for the oscillation angular frequency.
And the first comparison module 300 is configured to calculate participation factors of each node and branch to the sub-synchronous oscillation mode based on frequency domain mode analysis, and locate an influence area of the sub-synchronous oscillation by comparing the participation factors of each node and branch.
And the second comparison module 400 is used for establishing a sensitivity matrix according to the influence area of the subsynchronous oscillation, calculating the sensitivity of each element and positioning the key element influencing the oscillation by comparing the sensitivity of each element.
And a tracing module 500, configured to perform subsynchronous oscillation tracing according to the affected area and the key element of the subsynchronous oscillation.
It should be noted that the explanation of the embodiment of the source tracing method of sub-synchronous oscillation of the new energy system based on the frequency domain mode method is also applicable to the system, and is not repeated here.
According to the new energy system subsynchronous oscillation tracing system based on the frequency domain mode method, which is provided by the embodiment of the invention, an impedance network model of a target system is established, and the oscillation mode of the system is obtained based on a node admittance matrix and a loop impedance matrix of the network, so that the information of the oscillation damping and the frequency of subsynchronous oscillation is obtained; the participation factor indexes of the nodes/branches and the sensitivity indexes of the elements are defined, and the solution is carried out based on a frequency domain mode analysis method; according to the calculation results of the participation factors and the sensitivity, the main influence area of the oscillation and the key elements influencing the oscillation can be positioned, so that the oscillation source of the subsynchronous oscillation is identified.
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 such feature. In the description of the present invention, "a plurality" means at least two, e.g., two, three, etc., unless specifically limited otherwise.
in the present invention, unless otherwise expressly stated or limited, the terms "mounted," "connected," "secured," and the like are to be construed broadly and can, for example, be fixedly connected, detachably connected, or integrally formed; can be mechanically or electrically connected; they may be directly connected or indirectly connected through intervening media, or they may be connected internally or in any other suitable relationship, unless expressly stated otherwise. The specific meanings of the above terms in the present invention can be understood by those skilled in the art according to specific situations.
in the present invention, unless otherwise expressly stated or limited, the first feature "on" or "under" the second feature may be directly contacting the first and second features or indirectly contacting the first and second features through an intermediate. Also, a first feature "on," "over," and "above" a second feature may be directly or diagonally above the second feature, or may simply indicate that the first feature is at a higher level than the second feature. A first feature being "under," "below," and "beneath" a second feature may be directly under or obliquely under the first feature, or may simply mean that the first feature is at a lesser elevation than the second feature.
In the description herein, references to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., mean 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 invention. 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 more 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.
Although embodiments of the present invention have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present invention, and that variations, modifications, substitutions and alterations can be made to the above embodiments by those of ordinary skill in the art within the scope of the present invention.
Claims (10)
1. A new energy system subsynchronous oscillation tracing method based on a frequency domain mode method is characterized by comprising the following steps:
establishing an impedance network model of a target system, and forming a node admittance matrix and a loop impedance matrix of a network according to the impedance network model;
Calculating determinant zero points of the node admittance matrix and the loop impedance matrix to obtain an oscillation mode of a target system so as to position a subsynchronous oscillation mode;
calculating participation factors of each node and each branch to the subsynchronous oscillation mode based on frequency domain mode analysis, and positioning an influence area of subsynchronous oscillation by comparing the participation factors of each node and each branch;
Establishing a sensitivity matrix according to the affected area of the subsynchronous oscillation, calculating the sensitivity of each element, and positioning the key element affecting the oscillation by comparing the sensitivities of the elements;
and tracing the subsynchronous oscillation according to the affected area of the subsynchronous oscillation and the key element.
2. The method for tracing the subsynchronous oscillation of the new energy system based on the frequency domain mode method as claimed in claim 1, wherein the establishing of the impedance network model of the target system further comprises:
Establishing a frequency domain impedance and admittance model of each element in a target system;
And connecting the frequency domain impedance of each element according to the topological structure of the target system to obtain the impedance network model.
3. The method for tracing the subsynchronous oscillation of the new energy system based on the frequency domain mode method as recited in claim 1, wherein the calculating determinant zero points of the node admittance matrix and the loop impedance matrix to obtain the oscillation mode of the target system further comprises:
and calculating the numerator polynomial of the determinant of the node admittance matrix and the loop impedance matrix to obtain a plurality of pairs of complex zeros, wherein each pair of complex zeros corresponds to an oscillation mode.
4. The method for tracing the subsynchronous oscillation of the new energy system based on the frequency domain mode method as claimed in claim 3, wherein the oscillation mode is represented as:
si＝σi±jωi
wherein s isiis an oscillation mode, σiFor damping of oscillations, omegaiJ is the unit of an imaginary number for the oscillation angular frequency.
5. The method for tracing the subsynchronous oscillation of the new energy system based on the frequency domain pattern method according to claim 1, wherein the method for calculating the participation factors of each node and branch to the subsynchronous oscillation mode based on the frequency domain pattern analysis further comprises:
Calculating participation factors of each node to subsynchronous oscillation according to the node admittance matrix;
And calculating the participation factor of each branch to subsynchronous oscillation according to the loop impedance matrix.
6. The method for tracing the subsynchronous oscillation of the new energy system based on the frequency domain mode method as claimed in claim 1, wherein the calculating the sensitivity of each element further comprises:
Dividing the elements into series elements and parallel elements according to their positions in the nodal admittance matrix;
Calculating the sensitivity formula of the series element as follows: sen ═ Sjj-Sjk-Skj+SkkWhere S is the sensitivity matrix, j and k are connection nodes, SjkElement of jth row and kth column of S, Sjjof SRow jth column jth element, SkjIs the jth row and jth column element of S, Skkthe kth row and the kth column of S;
Calculating the sensitivity formula of the parallel element as follows: sen ═ SiiWhere S is the sensitivity matrix, i is a connection node, SiiIs the ith row and ith column element of S.
7. A subsynchronous oscillation traceability system of a new energy system based on a frequency domain mode method is characterized by comprising the following steps:
The system comprises an establishing module, a judging module and a judging module, wherein the establishing module is used for establishing an impedance network model of a target system and forming a node admittance matrix and a loop impedance matrix of a network according to the impedance network model;
The calculation module is used for calculating determinant zero points of the node admittance matrix and the loop impedance matrix, and solving an oscillation mode of a target system so as to position a subsynchronous oscillation mode;
The first comparison module is used for calculating participation factors of each node and each branch to the subsynchronous oscillation mode based on frequency domain mode analysis, and positioning an influence area of subsynchronous oscillation by comparing the participation factors of each node and each branch;
The second comparison module is used for establishing a sensitivity matrix according to the influence area of the subsynchronous oscillation, calculating the sensitivity of each element and positioning the key element influencing the oscillation by comparing the sensitivity of each element;
and the tracing module is used for tracing the subsynchronous oscillation according to the affected area of the subsynchronous oscillation and the key element.
8. The system according to claim 7, wherein the establishing module further includes:
The establishing unit is used for establishing frequency domain impedance and admittance models of all elements in the target system;
And the connecting unit is used for connecting the frequency domain impedance of each element according to the topological structure of the target system so as to obtain the impedance network model.
9. The system according to claim 7, wherein the calculation module further includes:
And calculating the numerator polynomial of the determinant of the node admittance matrix and the loop impedance matrix to obtain a plurality of pairs of complex zeros, wherein each pair of complex zeros corresponds to an oscillation mode, and a subsynchronous oscillation mode is found in the oscillation mode.
10. the frequency domain mode method-based sub-synchronous oscillation traceability system of the new energy system, as claimed in claim 9, wherein the oscillation mode is represented as:
si＝σi±jωi
Wherein s isiIs an oscillation mode, σifor damping of oscillations, omegaij is the unit of an imaginary number for the oscillation angular frequency.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910789949.7A CN110569577A (en) | 2019-08-26 | 2019-08-26 | Frequency domain mode method-based sub-synchronous oscillation tracing method and system for new energy system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910789949.7A CN110569577A (en) | 2019-08-26 | 2019-08-26 | Frequency domain mode method-based sub-synchronous oscillation tracing method and system for new energy system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110569577A true CN110569577A (en) | 2019-12-13 |
Family
ID=68775960
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910789949.7A Pending CN110569577A (en) | 2019-08-26 | 2019-08-26 | Frequency domain mode method-based sub-synchronous oscillation tracing method and system for new energy system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110569577A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111679124A (en) * | 2020-06-04 | 2020-09-18 | 北京交通大学 | Oscillation identification method and device for power system |
| CN112163532A (en) * | 2020-09-30 | 2021-01-01 | 国网冀北电力有限公司电力科学研究院 | Oscillation mode identification method and device of power system |
| CN112202211A (en) * | 2020-10-10 | 2021-01-08 | 华北电力大学 | Multi-VSG grid-connected power frequency oscillation evaluation method based on modal analysis method |
| CN114966296A (en) * | 2022-08-01 | 2022-08-30 | 锦浪科技股份有限公司 | Method and device for positioning oscillation source |
| CN116780502A (en) * | 2023-05-15 | 2023-09-19 | 华北电力大学 | Method and system for determining influence of power generation energy on low-frequency oscillation of power system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107769172A (en) * | 2017-11-09 | 2018-03-06 | 清华大学 | The system-level protection system and method for secondary/supersynchronous resonance/vibration |
| CN108462192A (en) * | 2018-04-28 | 2018-08-28 | 清华大学 | The pattern analysis method of electric system broadband oscillation |
-
2019
- 2019-08-26 CN CN201910789949.7A patent/CN110569577A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107769172A (en) * | 2017-11-09 | 2018-03-06 | 清华大学 | The system-level protection system and method for secondary/supersynchronous resonance/vibration |
| CN108462192A (en) * | 2018-04-28 | 2018-08-28 | 清华大学 | The pattern analysis method of electric system broadband oscillation |
Non-Patent Citations (1)
| Title |
|---|
| YING ZHAN等: "Frequency-Domain Modal Analysis of the", 《IEEE TRANSACTIONS ON SUSTAINABLE ENERGY》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111679124A (en) * | 2020-06-04 | 2020-09-18 | 北京交通大学 | Oscillation identification method and device for power system |
| CN111679124B (en) * | 2020-06-04 | 2021-07-23 | 北京交通大学 | Oscillation identification method and device for power system |
| CN112163532A (en) * | 2020-09-30 | 2021-01-01 | 国网冀北电力有限公司电力科学研究院 | Oscillation mode identification method and device of power system |
| CN112202211A (en) * | 2020-10-10 | 2021-01-08 | 华北电力大学 | Multi-VSG grid-connected power frequency oscillation evaluation method based on modal analysis method |
| CN114966296A (en) * | 2022-08-01 | 2022-08-30 | 锦浪科技股份有限公司 | Method and device for positioning oscillation source |
| CN114966296B (en) * | 2022-08-01 | 2022-12-02 | 锦浪科技股份有限公司 | Method and device for positioning oscillation source |
| CN116780502A (en) * | 2023-05-15 | 2023-09-19 | 华北电力大学 | Method and system for determining influence of power generation energy on low-frequency oscillation of power system |
| CN116780502B (en) * | 2023-05-15 | 2024-03-12 | 华北电力大学 | Method and system for determining influence of power generation energy on low-frequency oscillation of power system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110569577A (en) | Frequency domain mode method-based sub-synchronous oscillation tracing method and system for new energy system | |
| Aoki et al. | 2+ 1 flavor lattice QCD toward the physical point | |
| CN106443587B (en) | A kind of high-resolution quick deconvolution sound source imaging algorithm | |
| CN108462192B (en) | Mode analysis method for broadband oscillation of power system | |
| CN109657313A (en) | Missile performance appraisal procedure based on improved AHP method | |
| CN107871216A (en) | A kind of recognition methods of power distribution network fragility node | |
| CN105160058B (en) | A kind of network boundary equivalence method calculated for electric network electromagnet transient | |
| Ferro et al. | Anisotropic mesh adaptation for crack propagation induced by a thermal shock in 2D | |
| US20230048665A1 (en) | Two-stage frequency selection method and device for microwave frequency sweep data | |
| CN109951468A (en) | A kind of network attack detecting method and system based on the optimization of F value | |
| CN106410791B (en) | A kind of bad data is pushed forward back substitution method for tracing | |
| CN108319776A (en) | Simulation Parameters decision-making technique based on the soft collection of group's generalized interval intuitionistic fuzzy | |
| Guney et al. | Resonant frequency calculation for circular microstrip antennas with a dielectric cover using adaptive network-based fuzzy inference system optimized by various algorithms | |
| TWI452529B (en) | Combined with the system equivalent model of the system and its computer program products | |
| CN109308360A (en) | Channel parameter finding method, device, storage medium and computer equipment | |
| Hattori et al. | Theoretical and evolutionary parameter tuning of neural oscillators with a double-chain structure for generating rhythmic signals | |
| CN114977222A (en) | Analysis method for broadband oscillation influence factors and propagation paths of power system | |
| CN106803233B (en) | The optimization method of perspective image transformation | |
| CN113139295A (en) | Method and system for estimating comprehensive state of power system | |
| Melchior et al. | System identification of thermal transfers inside the lungs using fractional models | |
| CN109063928B (en) | Transformer oil density prediction model modeling method and device | |
| CN105633982B (en) | Power system stabilizer parameter tuning method and device based on mobile terminal | |
| He et al. | Optimization of renewable energy big data transactions based on vector evaluated genetic algorithm | |
| CN108921827A (en) | Stereo image quality evaluation method is referred to based on independent component analysis and the complete of binocular fusion characteristic | |
| Anh | Solving inverse problem in structural health monitoring by differential evolution algorithm with nearest neighbor comparison method |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191213 |
|
| RJ01 | Rejection of invention patent application after publication |