CN109861214B - Method and system for judging weak line with stable transient power angle of regional power grid - Google Patents

Method and system for judging weak line with stable transient power angle of regional power grid Download PDF

Info

Publication number
CN109861214B
CN109861214B CN201910152065.0A CN201910152065A CN109861214B CN 109861214 B CN109861214 B CN 109861214B CN 201910152065 A CN201910152065 A CN 201910152065A CN 109861214 B CN109861214 B CN 109861214B
Authority
CN
China
Prior art keywords
line
matrix
power grid
regional power
fault
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910152065.0A
Other languages
Chinese (zh)
Other versions
CN109861214A (en
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.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Shanxi Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Shanxi Electric Power 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 State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, State Grid Shanxi Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201910152065.0A priority Critical patent/CN109861214B/en
Publication of CN109861214A publication Critical patent/CN109861214A/en
Priority to PCT/CN2019/117585 priority patent/WO2020173132A1/en
Application granted granted Critical
Publication of CN109861214B publication Critical patent/CN109861214B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks

Abstract

The invention discloses a method and a system for judging weak lines of a regional power grid in transient state power angle stability, which can calculate line fragile factor LFF matrixes after different types of faults occur to different lines of a regional power grid by utilizing initial tide of the power grid through simple information acquisition and matrix operation, can directly judge weakest lines and severest faults of the system in transient state power angle stability through the fragile factor matrixes, have simple and direct methods, have clear influence on weak links and other lines of the system on the weak links, can greatly save simulation time, and can be used for various analyses of a power system.

Description

Method and system for judging weak line with stable transient power angle of regional power grid
Technical Field
The invention belongs to the technical field of power systems, and particularly relates to a method and a system for judging a weak line with a stable transient power angle of a regional power grid.
Background
At present, when a weak line with a stable transient power angle of a power system is searched, the influence on the system after the system is in fault is often more concerned, a simulation method is mainly adopted in a method for analyzing the influence of the fault, time is consumed, and the existing method for judging the weak section or the weak line is lack of analysis of the influence of the fault or is complex.
Disclosure of Invention
The present invention is directed to solving, to some extent, one of the technical problems in the related art.
Therefore, the first purpose of the invention is to provide a method for judging a weak line with a stable transient power angle of a regional power grid. The method can directly acquire the weak line and system faults sorted according to the influence severity of the weak line only according to certain state information of the power system, reduces the calculation time and the simulation time, and has important significance for transient power angle stability analysis of the power system.
The second purpose of the invention is to provide a system for judging the weak line with the stable transient power angle of the regional power grid.
In order to achieve the above object, a method for determining a weak line with a stable transient power angle of a regional power grid according to an embodiment of the first aspect of the present invention includes:
s1: and acquiring information of each line of the regional power grid and information of buses connected with each line of the regional power grid.
S2: calculating the power of each line of the regional power grid before the fault and calculating a matrix B consisting of admittance of each branch of each line of the regional power grid before the fault according to the information of each line of the regional power grid and the information of the bus connected with each line of the regional power grid acquired in the S1 L Node-branch incidence matrix A of each line of the regional power grid before fault and node admittance matrix B before fault 0
S3: adding branch circuits related to the connecting buses outside the regional power grid, setting admittance values corresponding to the branch circuits to be infinite, and correspondingly modifying a matrix B formed by admittance of each branch circuit in the S2 L Node-branch incidence matrix A of each line of the regional power grid before fault, and node admittance matrix B only considering reactance before fault 0
S4: correspondingly modifying the branch admittance matrix, the node-branch incidence matrix and the node admittance matrix in the S3 according to different set line faults to form an admittance matrix B 'of the branch after the fault' L Node-branch incidence matrix A ', node admittance matrix B' 0
S5: according to admittance matrix B 'of the branch in the S4' L Node-branch incidence matrix A ', node admittance matrix B' 0 And calculating the power of the fault line, the vulnerability factor of the fault line and the vulnerability factor matrix LEF of the fault line.
S6: calculating the vulnerability factor matrix absolute value LEF ' of the fault line in the S5, and arranging the vulnerability factor matrix absolute value LEF ' in a descending order, wherein the fault (the row corresponding to the maximum value of LEF ') corresponding to the maximum value in the LEF ' ordering is the most serious fault, and the line (the column corresponding to the maximum value of LEF ') corresponding to the line is the weakest line.
According to an embodiment of the present invention, the information of each line of the regional power grid in S1 includes: the name of a bus at two ends of the line, the voltage level at two ends of the line and the reactance value of the line; the bus information connected with each line of the regional power grid comprises the name of the bus connected with each line of the regional power grid, the voltage grade connected with each line of the regional power grid, the bus voltage connected with each line of the regional power grid and the phase angle of the bus connected with each line of the regional power grid.
According to an embodiment of the present invention, the matrix B of admittance of each branch in S2 L Node admittance matrix B 0 Is obtained considering only the reactance factor.
The system for judging the weak line with the stable transient power angle of the regional power grid provided by the embodiment of the second aspect of the invention comprises:
and the acquisition module is used for acquiring the information of each line of the regional power grid and the bus information connected with each line of the regional power grid.
A first calculating module, configured to calculate, according to the information of each line of the regional power grid and the information of the bus connected to each line of the regional power grid, the power of each line of the regional power grid before the fault, and calculate a matrix B formed by admittances of each branch of each line of the regional power grid before the fault L Node-branch incidence matrix A of each line of the regional power grid before fault and node admittance matrix B before fault 0
A first modification module, configured to add a branch associated with a connection bus outside the regional power grid, set an admittance value corresponding to the branch to infinity, and modify a matrix B formed by admittances of each branch in the S2 correspondingly L Node-branch incidence matrix A of each line of the regional power grid before fault, and node admittance matrix B only considering reactance before fault 0
A second modifying module, configured to modify the branch admittance matrix, the node-branch association matrix, and the node admittance matrix of the second computing module according to different set line faults, so as to form a faulted admittance matrix B 'of the branch' L Node-branch incidence matrix A ', node admittance matrix B' 0
A second calculation module for modifying the model according to the second modification modelAdmittance matrix B 'of the branches of a block' L Node-branch incidence matrix A ', node admittance matrix B' 0 And calculating the power of the fault line, the vulnerability factor of the fault line and the fault line vulnerability factor matrix LEF.
And the third calculation module is used for calculating the fault line vulnerability factor matrix absolute value LEF 'in the second calculation module, and sequencing the vulnerability factor matrix absolute values LEF' in a descending order, wherein the fault (the row corresponding to the maximum value of LEF ') corresponding to the maximum value in the LEF' sequencing is the most serious fault, and the line (the column corresponding to the maximum value of LEF ') corresponding to the maximum value of LEF' is the weakest line.
According to an embodiment of the present invention, the obtaining module obtains information of each line of the regional power grid, where the information includes: the name of a bus at two ends of the line, the voltage level at two ends of the line and the reactance value of the line; the bus information connected with each line of the regional power grid comprises the name of the bus connected with each line of the regional power grid, the voltage grade connected with each line of the regional power grid, the bus voltage connected with each line of the regional power grid and the bus phase angle connected with each line of the regional power grid.
According to an embodiment of the present invention, the matrix B formed by the admittance of each branch in the first calculation module L Node admittance matrix B 0 Is obtained considering only the reactance factor.
The method can calculate the line fragile factor LFF matrix after different types of faults occur to different lines of the power grid in a certain area by utilizing the initial tide of the power grid and through simple information acquisition and matrix operation, can directly judge the weakest line and the severest fault under the condition of stable transient power angle of the system through the fragile factor matrix, is simple and direct, has clear influence on the weak link and other lines of the system on the weak link, can greatly save simulation time, and can be used for various analyses of a power system.
Drawings
Exemplary embodiments of the invention may be more completely understood in consideration of the following drawings:
fig. 1 is a flowchart of the method for determining a weak line with a transient power angle stability of a regional power grid according to an embodiment of the disclosure;
FIG. 2 (a) is a graph of simulation of the maximum generator power angle difference for a system under different faults according to a disclosed embodiment of the invention;
FIG. 2 (b) is a graph of simulation of the maximum generator power angle difference for different faults in a system according to an embodiment of the disclosure;
fig. 3 is a block diagram of a system for determining a weak line in transient power angle stability of a regional power grid according to an embodiment of the present invention.
Detailed Description
In order to make the technical solutions of the present invention better understood, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be noted that the terms "first," "second," and the like in the description and claims of the present invention and in the drawings described above are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used is interchangeable under appropriate circumstances such that the embodiments of the invention described herein are capable of operation in sequences other than those illustrated or described herein. Moreover, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
The present invention will be described in further detail with reference to specific embodiments in order to make the above objects, features and advantages of the present invention more comprehensible.
The invention provides a method for acquiring weak lines and faults of transient power angle stability of a regional power grid, which is characterized in that based on a certain state of a system, a line weak factor matrix after different types of faults occur in the system is calculated, the weak lines of the power grid are judged, indexes of influence degrees of different faults on the weak lines are sequenced, the weak links of the power grid and the most serious faults of the system can be directly given, the calculation time and the simulation time can be reduced, and the method has a very clear guiding effect on transient power angle stability analysis of a power system.
In general, for a regional power grid, when a line k has a fault, a power flow transfer value on a line i is assumed to be Δ P ki Then, the ratio of the power of the line i to the maximum power after the tidal current is converted is:
Figure BDA0001981836690000051
according to a direct current power flow calculation formula:
P SP =B 0 θ (2)
Figure BDA0001981836690000052
let matrix A be a node-branch incidence matrix, then have
P L =B L A T θ (4)
P SP =B 0 θ (5)
P L =B L A T B 0 -1 P SP (6)
B 0 (AA T ) -1 AB L -1 P L =B 0 (AA T ) -1 AB L -1 B L A T B 0 -1 P SP =P SP (7)
When the line is disconnected:
Figure BDA0001981836690000053
then:
Figure BDA0001981836690000054
in the formula, P L As power vectors of branches before failure, B L A matrix formed by admittance of each branch considering reactance only before fault, A is a node-branch incidence matrix before fault, B 0 Node admittance matrix, P ', considering only reactances before failure' L Is the power of each branch after failure, B' L A matrix of branch admittances for which only the reactance is considered after a fault, A 'being a node-branch correlation matrix after a fault, B' 0 Only the reactive node admittance matrix is considered after a fault.
If it is assumed that each line delta max All are the same, then | P 'can be prepared' lk ·X′ lk And l is used as the weak factor index of each line after the k fault of the line, and the larger the index is, the weaker the line is. Wherein P' lk 、X′ lk Are each per unit values.
When the investigated regional power grid is electrically connected with the external regional power grid, the external regional power grid needs to be virtually connected through zero impedance before index calculation, so that the index calculation is more reasonable.
The premise of the application of the method is that the power grid is basically divided into layers and zones.
For a certain level of power grid, fig. 1 is a flowchart of the method for determining a stable weak line of a regional power grid after a fault according to the embodiment of the present disclosure, and the specific flow is as follows:
step S1: acquiring information of each line of the regional power grid and information of buses connected with each line of the regional power grid, namely acquiring information of each line of the regional power grid through the running state of a system, wherein the information comprises names of buses at two ends of the line, voltage levels and reactance values of the line; and acquiring bus information connected with the line, wherein the bus information comprises a bus name, a voltage grade, and bus voltage and phase angle in an operating state.
S2: calculating the power of each line of the regional power grid before the fault and calculating a matrix B consisting of admittance of each branch circuit before the fault of each line of the regional power grid according to the information of each line of the regional power grid and the information of the bus connected with each line of the regional power grid, which is acquired in the S1 L Node-branch incidence matrix A and node admittance matrix B of each line of the regional power grid before fault 0
S3: for a connecting bus outside an inspected regional power grid, adding branch circuits associated with the connecting bus outside the regional power grid, setting admittance values corresponding to the branch circuits to be infinite, and correspondingly modifying a matrix B formed by admittance of each branch circuit in the S2 L Node-branch incidence matrix A of each line of the regional power grid before fault, and node admittance matrix B only considering reactance before fault 0
S4: correspondingly modifying the branch admittance matrix, the node-branch incidence matrix and the node admittance matrix in the S3 according to different set line faults to form an admittance matrix B 'of the branch after the fault' L Node-branch incidence matrix A ', node admittance matrix B' 0
S5: admittance matrix B 'according to the branch in the S4' L Node-branch incidence matrix A ', node admittance matrix B' 0 And calculating the power of the fault line, the vulnerability factor of the fault line and the vulnerability factor matrix LEF of the fault line according to the formula (9).
S6: calculating the vulnerability factor matrix absolute value LEF ' of the fault line in the S5, and arranging the vulnerability factor matrix absolute value LEF ' in a descending order, wherein the fault (the row corresponding to the maximum value of LEF ') corresponding to the maximum value in the LEF ' ordering is the most serious fault, and the line (the column corresponding to the maximum value of LEF ') corresponding to the line is the weakest line.
The application verification is carried out by taking a certain planning power grid scheme of Jiangsu in the east China area as an example. Under the planning scheme, the east China power grid forms an extra-high voltage ring network, and there are 4 500kV connecting lines between Jiangsu power grids and external power grids. Under a certain extreme mode, other provinces in east China transmit power to Jiangsu, and the Jiangsu power grid reaches the maximum power receiving scale. At this time, the line vulnerability factor matrix under different line serious faults is calculated, and the line corresponding to the vulnerability factor with a larger absolute value is selected, and the result is shown in table 1.
TABLE 1 vulnerable line and its corresponding vulnerability factor in extreme mode
Vulnerable line Fragility factor
'Sujianfeng 51-Wan radix et rhizoma Rhei three 51' 0.682418
'Wan Dang Tu 51-Su Li Yang 51' 0.558321
'Suzhou south 51-SuLiyang 51' 0.358903
'Su Anlan 51-suduobin 51' 0.351139
'Suansuanlan 51-Suyudong 51' 0.329789
'Su Min Zhu 51-Su Jifeng 51' 0.30947
Table 1 shows that: the weakest line of the system is a 'suzufeng 51-Wanacu three 51' line, and then 'Wandangtu 51-SuLiyang 51'. And then the other lines in Jiangsu.
The fragile lines and corresponding fault lists sorted by the size of the fragile factors are as follows:
TABLE 2 vulnerable line sorted by vulnerable factor in extreme manner and corresponding fault
Figure BDA0001981836690000071
Figure BDA0001981836690000083
In order to check faults with large influence on fragile lines, fault lines in which several lines such as 'suzufeng 51' -wan umei 51',' wan dangui 51 '-suyuyang 51' and the like are sorted according to fragile factors after faults are respectively shown in tables 3 to 5.
TABLE 3 Fragile factor in extreme failure of SUZUIFENG-WanWU line and other lines
Figure BDA0001981836690000081
The faults in the table 3 are simulated, and the results show that the system is unstable after the line faults of the first two rows of 'Wan Dantu 51' -SuLiyang 51',' Suzhou south 51 '-SuLiyang 51', and the system is stable after the faults of the other lines. The correctness of the index is verified.
TABLE 4 vulnerable factor in extreme modes of "Wan Dang Tu 51-Su Li Yang 51" line and other line failures
Figure BDA0001981836690000082
Figure BDA0001981836690000091
TABLE 5 vulnerable factors in extreme modes of "Suzhou south 51-Su Li Yang 51" line and other line failures
Figure BDA0001981836690000092
The faults in the tables 4 and 5 are simulated, and the result shows that the system is unstable after the line fault of the first 'sujianfeng 51-wan turni three 51', and the system is stable after the fault of the other lines. The correctness of the index is verified.
Fig. 2 (a) and 2 (b) show simulation graphs of partial faults.
Fig. 3 is a block diagram of a system for determining a weak line of transient power angle stability of a regional power grid, where the system 300 includes:
an obtaining module 301, configured to obtain information of each line of the regional power grid and information of a bus connected to each line of the regional power grid.
A first calculating module 302, configured to calculate, according to the information of each line of the regional power grid and the information of the bus connected to each line of the regional power grid, which is obtained by the first obtaining module, the power of each line of the regional power grid before the fault, and calculate a matrix B formed by admittances of each branch of each line of the regional power grid before the fault L Node-branch incidence matrix A and node admittance matrix B of each line of the regional power grid before fault 0
A first modifying module 303, configured to add, to a connection bus outside an observed area power grid, a branch associated with the connection bus outside the area power grid, set an admittance value corresponding to the branch to infinity, and correspondingly modify a matrix B formed by admittances of each branch in S2 L Node-branch incidence matrix A of each line of the regional power grid before fault, and node admittance matrix B only considering reactance before fault 0
Second modification dieA block 304, configured to modify the branch admittance matrix, the node-branch association matrix, and the node admittance matrix of the second computing module according to different set line faults, to form a faulted admittance matrix B 'of the branch' L Node-branch incidence matrix A ', node admittance matrix B' 0
A second calculation module 305 for calculating an admittance matrix B 'from the branches of the second modification module' L Node-branch incidence matrix A ', node admittance matrix B' 0 And calculating the power of the fault line, the vulnerability factor of the fault line and the fault line vulnerability factor matrix LEF.
A third calculation module 306, configured to calculate the vulnerability factor matrix absolute values LEF ' of the faulty lines in the second calculation module, and arrange the vulnerability factor matrix absolute values LEF ' in a descending order, where the fault corresponding to the maximum value in the LEF ' ordering (the row corresponding to the maximum value of LEF ') is the most serious fault and the corresponding line (the column corresponding to the maximum value of LEF ') is the weakest line.
According to an embodiment of the present invention, the obtaining module 301 obtains information of each line of the regional power grid, including: the names of the buses at the two ends of the line, the voltage levels at the two ends of the line and the reactance value of the line; the bus information connected with each line of the regional power grid comprises the name of the bus connected with each line of the regional power grid, the voltage grade connected with each line of the regional power grid, the bus voltage connected with each line of the regional power grid and the bus phase angle connected with each line of the regional power grid.
According to an embodiment of the present invention, the matrix B formed by the admittance of each branch in the first calculation module 302 L Node admittance matrix B 0 Is obtained considering only the reactance factor.
The system for determining the weak line with the transient power angle stability of the regional power grid is the same as or similar to the method for determining the weak line with the transient power angle stability of the regional power grid in the embodiment of the first aspect, and is not described herein again.
The above-mentioned serial numbers of the embodiments of the present application are merely for description and do not represent the merits of the embodiments.
In the embodiments of the present application, the descriptions of the respective embodiments have respective emphasis, and for parts that are not described in detail in a certain embodiment, reference may be made to the related descriptions of other embodiments.
In the embodiments provided in the present application, it should be understood that the disclosed technology can be implemented in other ways. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units or modules is only one logical division, and there may be other divisions when the actual implementation is performed, for example, a plurality of units or modules or components may be combined or integrated into another system, or some features may be omitted or not executed. In addition, the shown or discussed coupling or direct coupling or communication connection between each other may be through some interfaces, indirect coupling or communication connection between modules or units, and may be in an electrical or other form.
The units or modules described as separate parts may or may not be physically separate, and parts displayed as units or modules may or may not be physical units or modules, may be located in one place, or may be distributed on a plurality of network units or modules. Some or all of the units or modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment.
In addition, functional units or modules in the embodiments of the present application may be integrated into one processing unit or module, or each unit or module may exist alone physically, or two or more units or modules are integrated into one unit or module. The integrated unit or module may be implemented in the form of hardware, or may also be implemented in the form of a software functional unit or module.
The integrated unit, if implemented in the form of a software functional unit and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application may be substantially implemented or contributed to by the prior art, or all or part of the technical solution may be embodied in a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a removable hard disk, a magnetic or optical disk, and other various media capable of storing program codes.
The foregoing is only a preferred embodiment of the present application and it should be noted that those skilled in the art can make several improvements and modifications without departing from the principle of the present application, and these improvements and modifications should also be considered as the protection scope of the present application.

Claims (4)

1. The method for judging the weak line with the stable transient power angle of the regional power grid is characterized by comprising the following steps:
s1: acquiring information of each line of a regional power grid and information of buses connected with each line of the regional power grid, wherein the information of each line of the regional power grid comprises names of buses at two ends of the line, voltage levels at two ends of the line and reactance values of the line; the bus information connected with each line of the regional power grid comprises a bus name connected with each line of the regional power grid, a voltage grade connected with each line of the regional power grid, a bus voltage connected with each line of the regional power grid and a bus phase angle connected with each line of the regional power grid;
s2: calculating the power of each line of the regional power grid before the fault and calculating a matrix B consisting of admittance of each branch circuit before the fault of each line of the regional power grid according to the information of each line of the regional power grid and the information of the bus connected with each line of the regional power grid obtained in the step S1 L Node-branch incidence matrix A and node admittance matrix B of each line of regional power grid before fault 0
S3: adding branches associated with connecting busbars outside a regional gridSetting the admittance value corresponding to the branch circuit to infinity, and correspondingly modifying the matrix B formed by the admittance of each branch circuit before the fault in S2 L Node-branch incidence matrix A and node admittance matrix B of each line of regional power grid before fault 0
S4: according to different set line faults, correspondingly modifying the matrix B formed by admittance of each branch circuit before the fault in S3 L Node-branch incidence matrix A and node admittance matrix B of each line of regional power grid before fault 0 Admittance matrix B 'of the branch after failure is formed' L Node-branch incidence matrix A ', node admittance matrix B' 0
S5: admittance matrix B 'according to branch in S4' L Node-branch incidence matrix A ', node admittance matrix B' 0 Calculating the power of the fault line, the fragile factor of the fault line and a fault line fragile factor matrix LEF;
s6: and calculating the fault line vulnerability factor matrix absolute value LEF 'in the S5, and arranging the vulnerability factor matrix absolute values LEF' in a descending order, wherein the fault corresponding to the maximum value in the LEF 'matrix in the LEF' ordering is the most serious fault, and the corresponding line is the weakest line.
2. The method of claim 1, wherein the admittance of each branch in S2 forms a matrix B L Node admittance matrix B 0 Is obtained considering only the reactance factor.
3. Judge system of regional electric wire netting transient state merit angle stable weak circuit, its characterized in that, the system includes:
the system comprises an acquisition module, a control module and a control module, wherein the acquisition module is used for acquiring information of each line of a regional power grid and information of buses connected with each line of the regional power grid, and the information of each line of the regional power grid comprises names of buses at two ends of the line, voltage levels at two ends of the line and reactance values of the line; the bus information connected with each line of the regional power grid comprises a bus name connected with each line of the regional power grid, a voltage grade connected with each line of the regional power grid, a bus voltage connected with each line of the regional power grid and a bus phase angle connected with each line of the regional power grid;
the first calculation module is used for calculating the power of each line of the regional power grid before the fault and calculating a matrix B consisting of admittance of each branch circuit before the fault of each line of the regional power grid according to the information of each line of the regional power grid acquired by the acquisition module and the information of the bus connected with each line of the regional power grid L Node-branch incidence matrix A and node admittance matrix B of each line of regional power grid before fault 0
A first modification module, configured to add a branch associated with a connection bus outside the regional power grid, set an admittance value corresponding to the branch to infinity, and correspondingly modify a matrix B formed by admittances of the branches in the first calculation module L Node-branch incidence matrix A and node admittance matrix B of each line of regional power grid before fault 0
A second modification module for correspondingly modifying the matrix B formed by the admittance of each branch before the fault in the first modification module according to the set different line faults L Node-branch incidence matrix A and node admittance matrix B of each line of regional power grid before fault 0 Admittance matrix B 'of the branch after failure is formed' L Node-branch incidence matrix A ', node admittance matrix B' 0
A second calculation module for calculating admittance matrix B 'according to branches of the second modification module' L Node-branch incidence matrix A ', node admittance matrix B' 0 Calculating the power of the fault line, the fragile factor of the fault line and a fault line fragile factor matrix LEF;
and the third calculation module is used for calculating the fault line vulnerability factor matrix absolute values LEF 'in the second calculation module, arranging the vulnerability factor matrix absolute values LEF' in a descending order, wherein the fault corresponding to the maximum value in the LEF 'matrix in the LEF' sorting is the most serious fault, and the corresponding line is the weakest line.
4. The system of claim 3, wherein each branch guide in the first computing moduleA nano-constituent matrix B L Node admittance matrix B 0 Is obtained considering only the reactance factor.
CN201910152065.0A 2019-02-28 2019-02-28 Method and system for judging weak line with stable transient power angle of regional power grid Active CN109861214B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201910152065.0A CN109861214B (en) 2019-02-28 2019-02-28 Method and system for judging weak line with stable transient power angle of regional power grid
PCT/CN2019/117585 WO2020173132A1 (en) 2019-02-28 2019-11-12 Method and system for determining transient power angle stability weak line of regional power grid

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910152065.0A CN109861214B (en) 2019-02-28 2019-02-28 Method and system for judging weak line with stable transient power angle of regional power grid

Publications (2)

Publication Number Publication Date
CN109861214A CN109861214A (en) 2019-06-07
CN109861214B true CN109861214B (en) 2022-10-04

Family

ID=66899350

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910152065.0A Active CN109861214B (en) 2019-02-28 2019-02-28 Method and system for judging weak line with stable transient power angle of regional power grid

Country Status (2)

Country Link
CN (1) CN109861214B (en)
WO (1) WO2020173132A1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109861214B (en) * 2019-02-28 2022-10-04 中国电力科学研究院有限公司 Method and system for judging weak line with stable transient power angle of regional power grid

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102403721A (en) * 2011-11-23 2012-04-04 昆明理工大学 Method for adjusting super real-time coincidence time sequence based on transient power angle stability margin
CN105656039A (en) * 2016-03-21 2016-06-08 国网宁夏电力公司电力科学研究院 Recognition method for vulnerable lines of power system

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005160257A (en) * 2003-11-27 2005-06-16 Japan Science & Technology Agency Single phase/three-phase converter, and its control method
SI21983A (en) * 2005-02-21 2006-08-31 DOMEL, Elketromotorji in gospodinjski aparati, d.d. Circuit and procedure for the control of reluctance motors
US10090777B2 (en) * 2011-05-08 2018-10-02 Koolbridge Solar, Inc. Inverter with independent current and voltage controlled outputs
CN103094901B (en) * 2013-01-09 2014-12-03 中国电力科学研究院 Method for identification of transient state power angle instability in real time based on actual measurement response message
CA2827701A1 (en) * 2013-09-23 2015-03-23 Sureshchandra B. Patel Methods of patel decoupled loadlow computation for electrical power system
WO2015044803A1 (en) * 2013-09-26 2015-04-02 Abb Technology Method for connecting subsystems of an electrical power system and an intelligent electronic device therefor
CN103762579B (en) * 2013-12-22 2015-11-11 国家电网公司 A kind of direct-current power emergency control performance index calculation method improving transient power angle stability of power grid
CN103795058B (en) * 2014-02-12 2015-10-28 国家电网公司 The air extract analysis of electric power system and system failure sort method
WO2017072593A1 (en) * 2015-10-26 2017-05-04 Rakesh Goel System and method for improved performance and synchronization of a grid tie inverter
CN106684855A (en) * 2015-11-10 2017-05-17 中国电力科学研究院 Transient stability emergency control method based on key branch identification
CN106130039A (en) * 2016-07-06 2016-11-16 南方电网科学研究院有限责任公司 The leading instability mode recognition method of power system and system
CN106786529B (en) * 2016-12-23 2019-05-03 国网山西省电力公司电力科学研究院 A kind of distribution static security analysis method
CN107317354B (en) * 2017-06-30 2019-12-17 天津大学 Method for analyzing transient power angle stability of multi-machine system with wind power plant
CN107785889B (en) * 2017-10-18 2020-08-21 中国电力科学研究院有限公司 Minimum startup optimization method for weak receiving-end power grid supporting high-capacity direct current feed-in
CN107895947B (en) * 2017-12-20 2023-07-07 广东电网有限责任公司惠州供电局 Power distribution network fragile line identification method
CN109861214B (en) * 2019-02-28 2022-10-04 中国电力科学研究院有限公司 Method and system for judging weak line with stable transient power angle of regional power grid

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102403721A (en) * 2011-11-23 2012-04-04 昆明理工大学 Method for adjusting super real-time coincidence time sequence based on transient power angle stability margin
CN105656039A (en) * 2016-03-21 2016-06-08 国网宁夏电力公司电力科学研究院 Recognition method for vulnerable lines of power system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Recognition and Vulnerability Analysis of Key Nodes in Power Grid Based on Complex Network Centrality;Bin Liu 等;《IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—II: EXPRESS BRIEFS》;20170518;第346-350页 *
Study on the transmission line vulnerability index and its application based on the transient power angle stability mechanism;Yujiao Sun 等;《2018 IEEE 2nd International Electrical and Energy Conference (CIEEC)》;20181106;第434-439页 *

Also Published As

Publication number Publication date
CN109861214A (en) 2019-06-07
WO2020173132A1 (en) 2020-09-03

Similar Documents

Publication Publication Date Title
Swetapadma et al. Improved fault location algorithm for multi‐location faults, transforming faults and shunt faults in thyristor controlled series capacitor compensated transmission line
Feng et al. Fault location using wide-area measurements and sparse estimation
WO2017016021A1 (en) State enumeration reliability evaluation method based on influence increment and device therefor
Barman et al. Detection and location of faults in large transmission networks using minimum number of phasor measurement units
CN106099891B (en) The analysis method that shunt reactor is distributed rationally is carried out to marine wind electric field leading-out terminal extra large cable
CN107271852A (en) Complicated Distribution Network Fault Locating Method based on voltage dip information
Soltani et al. Simultaneous robust state estimation, topology error processing, and outage detection for unbalanced distribution systems
Hou et al. Impact‐increment based decoupled reliability assessment approach for composite generation and transmission systems
CN106779269A (en) Reliability assessment method for secondary system of intelligent substation based on VC Method
CN112649642A (en) Electricity stealing position judging method, device, equipment and storage medium
Fan et al. Fault identification and location for distribution network with distributed generations
CN109861214B (en) Method and system for judging weak line with stable transient power angle of regional power grid
Yang et al. Three‐phase power‐flow solutions using decomposed quasi‐Newton method for unbalanced radial distribution networks
Alwash et al. New impedance‐based fault location method for unbalanced power distribution systems
CN103647272A (en) Alternating-current and direct-current power grid static equivalence method suitable for cascading failures
de Oliveira et al. Voltage sags: Validating short-term monitoring by using long-term stochastic simulation
Ghaedi et al. Modified WLS three-phase state estimation formulation for fault analysis considering measurement and parameter errors
Fan et al. Optimal fault location for power distribution systems with distributed generations using synchronized measurements
CN111339624B (en) PSASP and EMTP/ATP-based short-circuit current direct-current component calculation method
Wijeweera et al. Derivation of an equivalent circuit for real‐time security assessment
CN115983634A (en) Power distribution network operation risk assessment method and device
Kesherwani et al. Synchrophasor measurement‐based approach for online available transfer capability evaluation
CN110829413B (en) Radial distribution network reliability assessment method based on road matrix
Chen et al. Real‐time risk assessment of cascading failure in power system with high proportion of renewable energy based on fault graph chains
Mortazavi et al. An analytical fault location method based on minimum number of installed PMUs

Legal Events

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