CN112117757B - Key line identification method under N-k fault of power system considering information physical coupling relation - Google Patents

Key line identification method under N-k fault of power system considering information physical coupling relation Download PDF

Info

Publication number
CN112117757B
CN112117757B CN202010674779.0A CN202010674779A CN112117757B CN 112117757 B CN112117757 B CN 112117757B CN 202010674779 A CN202010674779 A CN 202010674779A CN 112117757 B CN112117757 B CN 112117757B
Authority
CN
China
Prior art keywords
line
fault
physical coupling
information
model
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
CN202010674779.0A
Other languages
Chinese (zh)
Other versions
CN112117757A (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.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN202010674779.0A priority Critical patent/CN112117757B/en
Publication of CN112117757A publication Critical patent/CN112117757A/en
Application granted granted Critical
Publication of CN112117757B publication Critical patent/CN112117757B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S40/00—Systems for electrical power generation, transmission, distribution or end-user application management characterised by the use of communication or information technologies, or communication or information technology specific aspects supporting them
    • Y04S40/20—Information technology specific aspects, e.g. CAD, simulation, modelling, system security

Abstract

The invention discloses a key line identification method under an N-k fault of an electric power system considering an information physical coupling relation in the field of electric power system safety analysis research. The method comprises the steps of firstly establishing an information physical coupling relation mathematical model according to the information physical coupling relation between a power grid and an information grid, then establishing a key line identification model under the N-k fault by using a multi-stage double-layer programming theory, and finally converting the model into a single-layer mixed integer linear programming problem by using a dual theory to solve and determine a key line. The invention considers the information physical coupling relation between the power network and the information network, so that the key line identification is more effective.

Description

Key line identification method under N-k fault of power system considering information physical coupling relation
Technical Field
The invention relates to the field of safety analysis of power systems, in particular to a method for identifying a key line under an N-k fault of a power system.
Background
Static safety analysis is one of the important means of traditional power system safety analysis, and is used for verifying the running state of a power system under an N-k fault, and further providing preventive measures or corrective measures to be taken before and after the fault. The national standard of the state of China, the electric power system safety and stability guide, clearly requires that the electric power system must meet static safety analysis. In recent years, a plurality of major power failure accidents at home and abroad show that part of key lines play an important role in large-scale cascading failure propagation of a power system. Therefore, line fault analysis is one of the important links of static safety analysis.
Static security analysis requires a very large number of predicted incidents to be verified, resulting in a large computational effort. Aiming at the problem, the main solution method comprises the following steps: (1) carrying out fault sequencing according to the possibility of overload, undervoltage or overvoltage of a system caused by the disconnection of each line, and then sequentially checking the lines with higher possibility according to the sequence; (2) developing an efficient static security analysis and calculation method according to the characteristics of the power system; (3) and identifying the key line according to the system load loss after the N-k fault. Compared with the former two methods, the method has the advantages that the worst operation state of the power system under the fault of the key line can be checked, and meanwhile, the key line can be used as a constraint condition to optimize the operation state of the system. The method is widely applied to the aspects of system risk assessment, safety constraint optimal power flow calculation, cascading failure prevention control, unit combination optimization and the like. Various advanced information communication technologies are continuously introduced into the current power system, so that the traditional power system is gradually developed into an information physical system, and important technical support is provided for the characteristics of flexibility, considerable controllability, sustainability, interoperability and the like of a power grid. However, the introduction of information communication technology enables a close information physical coupling relationship to exist between a power grid (primary system) and a communication grid (secondary system), and a fault on one side can weaken the capability of elements or systems on the other side against other disturbances, and even directly cause malfunction or shutdown of corresponding elements. However, the above studies do not consider the influence of the physical coupling relationship between the power grid and the communication grid on the vulnerability of the system.
Disclosure of Invention
The invention aims to solve the technical problem of providing a method and a system for identifying a key line under an N-k fault of an electric power system by considering an information physical coupling relation, aiming at the defect that the information physical coupling relation between an electric power network and an information network is not considered in the prior art.
The technical scheme adopted by the invention for solving the technical problems is as follows:
the method for identifying the key line under the N-k fault of the power system by considering the information physical coupling relation comprises the following steps:
step 1: establishing an information physical coupling relation mathematical model according to an information physical coupling relation between a power grid and an information grid;
step 2: establishing a key line identification model under the N-k fault by using a multi-stage double-layer programming theory; assuming a variable v of 0-1lIndicating the open state of line l, if line l is open, vl0, otherwise vlThe key line under the fault of 1, N-k is the k line l which has the largest influence on the system;
and step 3: and converting the key line identification model under the N-k fault into a single-layer mixed integer linear programming problem by using a dual theory to solve, and determining the key line.
According to the technical scheme, the key line identification model under the N-k fault specifically comprises the following steps:
Figure GDA0002757524120000021
Figure GDA0002757524120000022
Figure GDA0002757524120000023
Figure GDA0002757524120000024
Figure GDA0002757524120000025
Figure GDA0002757524120000026
Figure GDA0002757524120000027
Figure GDA0002757524120000028
Figure GDA0002757524120000029
Figure GDA00027575241200000210
Figure GDA0002757524120000031
Figure GDA0002757524120000032
Figure GDA0002757524120000033
Figure GDA0002757524120000034
Figure GDA0002757524120000035
in the formula (I), the compound is shown in the specification,
Figure GDA0002757524120000036
wherein v islIs a variable 0-1, indicating that line l is open, if line l is open, vl0, otherwise vl＝1，Pg/Pg,s,Fl/Fl,s,Sd/ΔDd,s,Dd,s,θdThe variables are continuous variables and respectively represent the output of the generator g, the load flow of the line l, the load cut off from the bus d, the actual load on the bus d and the phase angle of the system in the initial stage (s is 1) and the stage s is more than or equal to 2. Variables of
Figure GDA0002757524120000037
And
Figure GDA0002757524120000038
represents the optimal generator output and shear load, x, determined by the underlying planning modellDenotes the reactance of the line l, Nd，NlAnd NsRespectively representing the number of busesNumber of line and cascading failures propagation, G and GdRespectively representing the set of generators connected to the bus bar d,
Figure GDA0002757524120000039
representing a node-branch incidence matrix, if the bus d is the origin of the line l, AdlIf bus d is the end of line l, a is 1dlNot being equal to-1, otherwise Adl0; and O represents a power primary system bus set, a generator and a load are connected, and L represents a transmission line set.
According to the technical scheme, the mathematical model of the information physical coupling relation is specifically as follows:
Figure GDA00027575241200000310
Figure GDA00027575241200000311
in which the variables are 0-1
Figure GDA00027575241200000312
Indicating the working state of the information node v in the cascading failure propagation stage s, and if the information node is working normally,
Figure GDA00027575241200000313
otherwise
Figure GDA00027575241200000314
The generator corresponding to the information node v is g,
Figure GDA00027575241200000315
represents the generator output, P, determined by the dispatching center according to the optimal power flow model after the line faultg,sRepresenting the actual output of the generator at stage s. The normal load on the bus D of the primary power system is DdAnd the minimum power supply requirement of the information node v corresponding to the bus D is represented as alpha DdI.e. only when on the busActual load greater than or equal to α DdThe information node can normally work, alpha represents an energy control coefficient, alpha is more than or equal to 0 and less than or equal to 1 and can be used for representing the information physical coupling strength, the larger the alpha is, the larger the control-energy information physical coupling strength is, otherwise, the smaller the information physical coupling strength is. Dd,sRepresents the actual load of the bus d in the phase s, BvRepresenting the set of busbars providing power to the information node v. The function I (f (x)) represents an indicator function, I (f (x)) being 1 if f (x) ≦ 0, otherwise I (f (x)) being 0.
According to the technical scheme, the more critical the line is, the more the load is
Figure GDA00027575241200000316
The greater the loss.
According to the technical scheme, the step 3 is specifically as follows:
the key line identification model under the N-k fault established by using the multi-stage double-layer programming theory is expressed by adopting a compact form as shown in the following:
Figure GDA0002757524120000041
s.t.g(γ,x,y*)≥0 (2)
γ∈{0,1} (3)
Figure GDA0002757524120000042
s.t.A+B(γ)y＝b1:ψ (5)
C+Dy≥b2:ζ (6)
in this model, f (-) and dTy represents upper and lower optimization model objective functions, γ and x represent upper optimization model 0-1 variables and continuous variables, y represents lower optimization model variables, ψ and ζ represent lagrangian multipliers corresponding to expressions (5) and (6), D, a, B (γ), C, D, B1And b2Respectively representing coefficient matrix and constant term of corresponding constraint, B (gamma) indicating the coefficient matrix and parameterGamma correlation; the formula (2) represents the constraint condition of the upper-layer optimization model, and the formulas (5) and (6) represent the equality and inequality constraints of the lower-layer optimization model respectively;
by using the dual theory, the model can be converted into a single-layer optimization model as shown in the following:
Figure GDA0002757524120000043
s.t.g(γ,x,y)≥0
γ∈{0,1}
B(γ)Tψ+DTζ＝d
ζ≥0
A+B(γ)y＝b1
C+Dy≥b2
dTy＝ψT[b1-A]+ζT[b2-C]。
the invention also provides a key line identification system under the N-k fault of the power system considering the information physical coupling relation, which comprises the following steps:
the information physical coupling relation mathematical model establishing module is used for establishing an information physical coupling relation mathematical model according to the information physical coupling relation between the power grid and the information grid;
the key line identification model establishing module is used for establishing a key line identification model under the N-k fault by utilizing a multi-stage double-layer planning theory; assuming a variable v of 0-1lIndicating the open state of line l, if line l is open, vl0, otherwise vlThe key line under the fault of 1, N-k is the k line l which has the largest influence on the system;
and the key line solving module is used for converting the key line identification model under the N-k fault into a single-layer mixed integer linear programming problem by utilizing a dual theory to solve and determine the key line.
According to the technical scheme, the key line identification model under the N-k fault specifically comprises the following steps:
Figure GDA0002757524120000051
Figure GDA0002757524120000052
Figure GDA0002757524120000053
Figure GDA0002757524120000054
Figure GDA0002757524120000055
Figure GDA0002757524120000056
Figure GDA0002757524120000057
Figure GDA0002757524120000058
Figure GDA0002757524120000059
Figure GDA00027575241200000510
Figure GDA00027575241200000511
Figure GDA00027575241200000512
Figure GDA00027575241200000513
Figure GDA00027575241200000514
Figure GDA00027575241200000515
in the formula (I), the compound is shown in the specification,
Figure GDA00027575241200000516
wherein v islIs a variable 0-1, indicating that line l is open, if line l is open, vl0, otherwise vl＝1，Pg/Pg,s,Fl/Fl,s,Sd/ΔDd,s,Dd,s,θdThe variables are continuous variables and respectively represent the output of the generator g, the load flow of the line l, the load cut off from the bus d, the actual load on the bus d and the phase angle of the system in the initial stage (s is 1) and the stage s is more than or equal to 2. Variables of
Figure GDA00027575241200000517
And
Figure GDA00027575241200000518
represents the optimal generator output and shear load, x, determined by the underlying planning modellDenotes the reactance of the line l, Nd，NlAnd NsRespectively representing the number of buses, lines and the number of cascading failure transmissions, G and GdRespectively representing the set of generators connected to the bus bar d,
Figure GDA00027575241200000519
representing node-branch incidence matrix if bus d is origin of line lThen A isdlIf bus d is the end of line l, a is 1dlNot being equal to-1, otherwise Adl0; and O represents a power primary system bus set, a generator and a load are connected, and L represents a transmission line set.
According to the technical scheme, the mathematical model of the information physical coupling relation is specifically as follows:
Figure GDA0002757524120000061
Figure GDA0002757524120000062
in which the variables are 0-1
Figure GDA0002757524120000063
Indicating the working state of the information node v in the cascading failure propagation stage s, and if the information node is working normally,
Figure GDA0002757524120000064
otherwise
Figure GDA0002757524120000065
The generator corresponding to the information node v is g,
Figure GDA0002757524120000066
represents the generator output, P, determined by the dispatching center according to the optimal power flow model after the line faultg,sRepresenting the actual output of the generator at stage s. The normal load on the bus D of the primary power system is DdAnd the minimum power supply requirement of the information node v corresponding to the bus D is represented as alpha DdI.e. only when the actual load on the bus is greater than or equal to α DdThe information node can normally work, alpha represents an energy control coefficient, alpha is more than or equal to 0 and less than or equal to 1 and can be used for representing the information physical coupling strength, the larger the alpha is, the larger the control-energy information physical coupling strength is, otherwise, the smaller the information physical coupling strength is. Dd,sIndicating bus d in phases actual load, BvRepresenting the set of busbars providing power to the information node v. The function I (f (x)) represents an indicator function, I (f (x)) being 1 if f (x) ≦ 0, otherwise I (f (x)) being 0.
According to the technical scheme, the key line solving module is specifically used for:
the key line identification model under the N-k fault established by using the multi-stage double-layer programming theory is expressed by adopting a compact form as shown in the following:
Figure GDA0002757524120000067
s.t.g(γ,x,y*)≥0 (2)
γ∈{0,1} (3)
Figure GDA0002757524120000068
s.t.A+B(γ)y＝b1:ψ (5)
C+Dy≥b2:ζ (6)
in this model, f (-) and dTy represents upper and lower optimization model objective functions, γ and x represent upper optimization model 0-1 variables and continuous variables, y represents lower optimization model variables, ψ and ζ represent lagrangian multipliers corresponding to expressions (5) and (6), D, a, B (γ), C, D, B1And b2Respectively representing a coefficient matrix and a constant item corresponding to the constraint, wherein B (gamma) indicates that the coefficient matrix is related to a parameter gamma; the formula (2) represents the constraint condition of the upper-layer optimization model, and the formulas (5) and (6) represent the equality and inequality constraints of the lower-layer optimization model respectively;
by using the dual theory, the model can be converted into a single-layer optimization model as shown in the following:
Figure GDA0002757524120000071
s.t.g(γ,x,y)≥0
γ∈{0,1}
B(γ)Tψ+DTζ＝d
ζ≥0
A+B(γ)y＝b1
C+Dy≥b2
dTy＝ψT[b1-A]+ζT[b2-C]。
the invention also provides a computer storage medium, which stores a computer program capable of being executed by a processor, wherein the computer program specifically executes the method for identifying the critical line under the N-k fault of the power system considering the information physical coupling relation in the technical scheme.
The invention has the following beneficial effects: according to the method, the key line identification is more effective due to the consideration of the information physical coupling relation between the power grid and the information grid, so that the key line under the N-k fault of the power system considering the information physical coupling relation can be effectively identified, and the N-k fault causing the maximum load loss is determined.
Drawings
The invention will be further described with reference to the accompanying drawings and examples, in which:
FIG. 1 is a flowchart of a method for identifying a critical line under an N-k fault of an electric power system in consideration of an information physical coupling relationship according to an embodiment of the present invention;
FIG. 2 is a multi-stage dual level programming diagram according to an embodiment of the present invention;
FIG. 3 is a schematic structural diagram of a critical line identification system under an N-k fault of an electrical power system in consideration of an information physical coupling relationship according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
As shown in fig. 1, the method for identifying a critical line under an N-k fault of an electric power system considering an information physical coupling relationship according to the embodiment of the present invention includes the following steps:
step 1: and establishing an information physical coupling relation mathematical model according to the information physical coupling relation between the power grid and the information grid. The information physical coupling relation between the power grid and the information grid means that the safe and reliable operation of the power primary system is controlled by the information node, and the normal operation of the information node in the power information communication system is powered by the power primary system. When the primary power system fails, the frequency of the power grid fluctuates. For a power primary system, when an information node fails, the power primary equipment may fail due to loss of control. In particular, for a generator, the generator loses control due to the failure of the corresponding information node, and finally exits operation because of low-frequency or over-frequency protection.
The physical coupling relationship of the information can be determined according to the following mathematical model:
Figure GDA0002757524120000081
Figure GDA0002757524120000082
in which the variables are 0-1
Figure GDA0002757524120000083
Indicating the working state of the information node v in the cascading failure propagation stage s, and if the information node is working normally,
Figure GDA0002757524120000084
otherwise
Figure GDA0002757524120000085
The generator corresponding to the information node v is g,
Figure GDA0002757524120000086
represents the generator output, P, determined by the dispatching center according to the optimal power flow model after the line faultg,sRepresenting the actual output of the generator at stage s. Electric power is onceThe normal load on the system bus D is DdAnd the lowest power supply requirement of the information node v corresponding to the bus D is represented as alpha DdI.e. only when the actual load on the bus is greater than or equal to α DdThe information node can normally work, alpha represents an energy control coefficient, alpha is more than or equal to 0 and less than or equal to 1 and can be used for representing the information physical coupling strength, the larger the alpha is, the larger the control-energy information physical coupling strength is, otherwise, the smaller the information physical coupling strength is. Dd,sRepresents the actual load of the bus d in the phase s, BvRepresenting the set of busbars providing power to the information node v. The function I (f (x)) represents an indicator function, I (f (x)) being 1 if f (x) ≦ 0, otherwise I (f (x)) being 0.
Step 2: considering cascading failures between a power grid and a communication grid, a key line identification model under the N-k failures is established by utilizing a multi-stage double-layer planning theory, and the objective function of the model is maximum load shedding.
Specifically, as shown in fig. 2, the upper layer model represents a line fault and determines a faulty line, and the lower layer model represents a dispatching center and determines an operation power flow of the power grid. Because the existing optimal power flow model does not consider the information physical coupling relation, cascading failures can occur between the primary power system and the power information communication system. The key line identification model mathematical model is as follows:
Figure GDA0002757524120000091
Figure GDA0002757524120000092
Figure GDA0002757524120000093
Figure GDA0002757524120000094
Figure GDA0002757524120000095
Figure GDA0002757524120000096
Figure GDA0002757524120000097
Figure GDA0002757524120000098
Figure GDA0002757524120000099
Figure GDA00027575241200000910
Figure GDA00027575241200000911
Figure GDA00027575241200000912
Figure GDA00027575241200000913
Figure GDA00027575241200000914
Figure GDA00027575241200000915
in the formula (I), the compound is shown in the specification,
Figure GDA00027575241200000916
wherein v islIs a variable 0-1, indicating that line l is open, if line l is open, vl0, otherwise vl＝1，Pg/Pg,s,Fl/Fl,s,Sd/ΔDd,s,Dd,s,θdThe variables are continuous variables and respectively represent the output of the generator g, the load flow of the line l, the load cut off from the bus d, the actual load on the bus d and the phase angle of the system in the initial stage (s is 1) and the stage s is more than or equal to 2. Variables of
Figure GDA00027575241200000917
And
Figure GDA00027575241200000918
represents the optimal generator output and shear load, x, determined by the underlying planning modellDenotes the reactance of the line l, Nd，NlAnd NsRespectively representing the number of buses, lines and the number of cascading failure transmissions, G and GdRespectively representing the set of generators connected to the bus bar d,
Figure GDA00027575241200000919
representing a node-branch incidence matrix, if the bus d is the origin of the line l, AdlIf bus d is the end of line l, a is 1dlNot being equal to-1, otherwise Adl0. O denotes a set of primary system buses (power nodes) for power, connecting generators and loads, and L denotes a set of transmission lines. The more critical the line, the load
Figure GDA00027575241200000920
The greater the loss.
And step 3: converting the model into a single-layer mixed integer linear programming problem by using a dual theory to solve, and determining a key line vl。
The key line identification model under the N-k fault established by using the multi-stage double-layer programming theory is expressed by adopting a compact form as shown in the following:
Figure GDA0002757524120000101
s.t.g(γ,x,y*)≥0 (2)
γ∈{0,1} (3)
Figure GDA0002757524120000102
s.t.A+B(γ)y＝b1:ψ (5)
C+Dy≥b2:ζ (6)
in this model, f (-) and dTy represents upper and lower optimization model objective functions, γ and x represent upper optimization model 0-1 variables and continuous variables, y represents lower optimization model variables, ψ and ζ represent lagrangian multipliers corresponding to expressions (5) and (6), D, a, B (γ), C, D, B1And b2Respectively representing a coefficient matrix and a constant item corresponding to the constraint, wherein B (gamma) indicates that the coefficient matrix is related to a parameter gamma; the formula (2) represents the constraint condition of the upper-layer optimization model, and the formulas (5) and (6) represent the equality and inequality constraints of the lower-layer optimization model respectively;
by using the dual theory, the model can be converted into a single-layer optimization model as shown in the following:
Figure GDA0002757524120000103
s.t.g(γ,x,y)≥0
γ∈{0,1}
B(γ)Tψ+DTζ＝d
ζ≥0
A+B(γ)y＝b1
C+Dy≥b2
dTy＝ψT[b1-A]+ζT[b2-C]。
the invention also provides a key line identification system under the N-k fault of the power system considering the information physical coupling relationship, as shown in fig. 3, the system comprises:
the information physical coupling relation mathematical model establishing module is used for establishing an information physical coupling relation mathematical model according to the information physical coupling relation between the power grid and the information grid;
the key line identification model establishing module is used for establishing a key line identification model under the N-k fault by utilizing a multi-stage double-layer planning theory; assume that a variable v of 0-1 is usedlIndicating the open state of line l, if line l is open, vl0, otherwise vlThe key line under the fault of 1, N-k is the k line l which has the largest influence on the system;
and the key line solving module is used for converting the key line identification model under the N-k fault into a single-layer mixed integer linear programming problem by utilizing a dual theory to solve and determine the key line.
The specific implementation manner of each module is as in the above embodiment, and the identification of the critical line is not repeated herein, when the N-k fault of the power system is considered in the information physical coupling relationship.
The computer storage medium of the embodiment of the present invention stores therein a computer program executable by a processor, where the computer program specifically executes the method for identifying a critical line in an N-k fault of an electrical power system considering an information physical coupling relationship according to the above embodiment, and the specific program is not described herein again.
It will be understood that modifications and variations can be made by persons skilled in the art in light of the above teachings and all such modifications and variations are intended to be included within the scope of the invention as defined in the appended claims.

Claims (8)

1. A method for identifying a key line under an N-k fault of an electric power system by considering an information physical coupling relation is characterized by comprising the following steps:
step 1: establishing an information physical coupling relation mathematical model according to an information physical coupling relation between a power grid and an information grid;
step 2: establishing a key line identification model under the N-k fault by using a multi-stage double-layer programming theory; assume that a variable v of 0-1 is usedlIndicating the open state of line l, if line l is open, vl0, otherwise vlThe key line under the fault of 1, N-k is the k line l which has the largest influence on the system;
and step 3: converting the key line identification model under the N-k fault into a single-layer mixed integer linear programming problem by using a dual theory to solve and determine a key line;
the key line identification model under the N-k fault is specifically as follows:
Figure FDA0003541342860000011
Figure FDA0003541342860000012
Figure FDA0003541342860000013
Figure FDA0003541342860000014
Figure FDA0003541342860000015
Figure FDA0003541342860000016
Figure FDA0003541342860000017
Figure FDA0003541342860000018
Figure FDA0003541342860000019
Figure FDA00035413428600000110
Figure FDA00035413428600000111
Figure FDA00035413428600000112
Figure FDA0003541342860000021
Figure FDA0003541342860000022
Figure FDA0003541342860000023
in the formula (I), the compound is shown in the specification,
Figure FDA0003541342860000024
wherein P isg/Pg,s,Fl/Fl,s,Sd/ΔDd,s,Dd,s,θdAll are continuous variables, which respectively represent the output of the generator g, the load flow of the line l and the bus d of the system at the initial stage (s is 1) and the stage s is more than or equal to 2The load cut off, the actual load on the bus d and the phase angle; variables of
Figure FDA0003541342860000025
And
Figure FDA0003541342860000026
represents the optimal generator output and shear load, x, determined by the underlying planning modellDenotes the reactance of the line l, Nd，NlAnd NsRespectively representing the number of buses, lines and the number of cascading failure transmissions, G and GdRespectively representing the set of generators connected to the bus bar d,
Figure FDA0003541342860000027
representing a node-branch incidence matrix, if the bus d is the origin of the line l, AdlIf bus d is the end of line l, a is 1dlNot being equal to-1, otherwise Adl0; and O represents a power primary system bus set, a generator and a load are connected, and L represents a transmission line set.
2. The method for identifying the key line under the N-k fault of the power system considering the information physical coupling relationship as claimed in claim 1, wherein the mathematical model of the information physical coupling relationship is specifically as follows:
Figure FDA0003541342860000028
Figure FDA0003541342860000029
in which the variables are 0-1
Figure FDA00035413428600000210
Indicating the working state of the information node v in the cascading failure propagation stage s, and if the information node is working normally,
Figure FDA00035413428600000211
otherwise
Figure FDA00035413428600000212
The generator corresponding to the information node v is g,
Figure FDA00035413428600000213
represents the generator output, P, determined by the dispatching center according to the optimal power flow model after the line faultg,sRepresenting the actual output of the generator at stage s; the normal load on the bus D of the primary power system is DdAnd the minimum power supply requirement of the information node v corresponding to the bus D is represented as alpha DdI.e. only when the actual load on the bus is greater than or equal to α DdThe information node can normally work, alpha represents an energy control coefficient, alpha is more than or equal to 0 and less than or equal to 1 and can be used for representing the information physical coupling strength, the larger the alpha is, the larger the control-energy information physical coupling strength is, otherwise, the smaller the information physical coupling strength is; dd,sRepresents the actual load of the bus d in the phase s, BvA set of busbars representing the supply of power to the information node v; the function I (f (x)) represents an indicator function, I (f (x)) being 1 if f (x) ≦ 0, otherwise I (f (x)) being 0.
3. The method for identifying key lines under N-k faults of electric power system in consideration of information physical coupling relation as claimed in claim 1, wherein the more critical a line is, the more load is
Figure FDA00035413428600000214
The greater the loss.
4. The method for identifying the key line under the N-k fault of the power system considering the information physical coupling relationship according to claim 1, wherein the step 3 specifically comprises:
the key line identification model under the N-k fault established by using the multi-stage double-layer programming theory is expressed by adopting a compact form as shown in the following:
Figure FDA0003541342860000031
s.t.g(γ,x,y*)≥0 (2)
γ∈{0,1} (3)
Figure FDA0003541342860000032
s.t.A+B(γ)y＝b1:ψ (5)
C+Dy≥b2:ζ (6)
in this model, f (-) and dTy represents upper and lower optimization model objective functions, γ and x represent upper optimization model 0-1 variables and continuous variables, y represents lower optimization model variables, ψ and ζ represent lagrangian multipliers corresponding to expressions (5) and (6), D, a, B (γ), C, D, B1And b2Respectively representing a coefficient matrix and a constant item corresponding to the constraint, wherein B (gamma) indicates that the coefficient matrix is related to a parameter gamma; the formula (2) represents the constraint condition of the upper-layer optimization model, and the formulas (5) and (6) represent the equality and inequality constraints of the lower-layer optimization model respectively;
by using the dual theory, the model can be converted into a single-layer optimization model as shown in the following:
Figure FDA0003541342860000034
s.t.g(γ,x,y)≥0
γ∈{0,1}
B(γ)Tψ+DTζ＝d
ζ≥0
A+B(γ)y＝b1
C+Dy≥b2
dTy＝ψT[b1-A]+ζT[b2-C]。
5. a critical line identification system under N-k faults of a power system considering information physical coupling relation, which is characterized by comprising:
the information physical coupling relation mathematical model establishing module is used for establishing an information physical coupling relation mathematical model according to the information physical coupling relation between the power grid and the information grid;
the key line identification model establishing module is used for establishing a key line identification model under the N-k fault by utilizing a multi-stage double-layer planning theory; assuming a variable v of 0-1lIndicating the open state of line l, if line l is open, vl0, otherwise vlThe key line under the fault of 1, N-k is the k line l which has the largest influence on the system;
the key line solving module is used for converting the key line identification model under the N-k fault into a single-layer mixed integer linear programming problem by utilizing a dual theory to solve and determine a key line;
the key line identification model under the N-k fault is specifically as follows:
Figure FDA0003541342860000041
Figure FDA0003541342860000042
Figure FDA0003541342860000043
Figure FDA0003541342860000044
Figure FDA0003541342860000045
Figure FDA0003541342860000046
Figure FDA0003541342860000047
Figure FDA0003541342860000048
Figure FDA0003541342860000049
Figure FDA00035413428600000410
Figure FDA00035413428600000411
Figure FDA00035413428600000412
Figure FDA00035413428600000413
Figure FDA00035413428600000414
Figure FDA00035413428600000415
in the formula (I), the compound is shown in the specification,
Figure FDA00035413428600000416
wherein P isg/Pg,s,Fl/Fl,s,Sd/ΔDd,s,Dd,s,θdThe variable values are continuous variables and respectively represent the output of a generator g, the load flow of a line l, the load cut off from a bus d, the actual load on the bus d and the phase angle of the system at an initial stage (s is 1) and a stage s is more than or equal to 2; variables of
Figure FDA00035413428600000417
And
Figure FDA00035413428600000418
represents the optimal generator output and shear load, x, determined by the underlying planning modellDenotes the reactance of the line l, Nd，NlAnd NsRespectively representing the number of buses, lines and the number of cascading failure transmissions, G and GdRespectively representing the set of generators connected to the bus bar d,
Figure FDA00035413428600000419
representing a node-branch incidence matrix, if the bus d is the origin of the line l, AdlIf bus d is the end of line l, a is 1dlNot being equal to-1, otherwise Adl0; and O represents a power primary system bus set, a generator and a load are connected, and L represents a transmission line set.
6. The key line identification system under the N-k fault of the electric power system considering the cyber-physical coupling relationship according to claim 5, wherein the cyber-physical coupling relationship mathematical model is specifically:
Figure FDA0003541342860000051
Figure FDA0003541342860000052
in which the variables are 0-1
Figure FDA0003541342860000053
Indicating the working state of the information node v in the cascading failure propagation stage s, and if the information node is working normally,
Figure FDA0003541342860000054
otherwise
Figure FDA0003541342860000055
The generator corresponding to the information node v is g,
Figure FDA0003541342860000056
represents the generator output, P, determined by the dispatching center according to the optimal power flow model after the line faultg,sRepresenting the actual output of the generator at stage s; the normal load on the bus D of the primary power system is DdAnd the minimum power supply requirement of the information node v corresponding to the bus D is represented as alpha DdI.e. only when the actual load on the bus is greater than or equal to α DdThe information node can normally work, alpha represents an energy control coefficient, alpha is more than or equal to 0 and less than or equal to 1 and can be used for representing the information physical coupling strength, the larger the alpha is, the larger the control-energy information physical coupling strength is, otherwise, the smaller the information physical coupling strength is; dd,sRepresents the actual load of the bus d in the phase s, BvA set of busbars representing the supply of power to the information node v; the function I (f (x)) represents an indicator function, I (f (x)) being 1 if f (x) ≦ 0, otherwise I (f (x)) being 0.
7. The key line identification system under the N-k fault of the power system considering the information physical coupling relation according to claim 5, wherein the key line solving module is specifically configured to:
the key line identification model under the N-k fault established by using the multi-stage double-layer programming theory is expressed by adopting a compact form as shown in the following:
Figure FDA0003541342860000057
s.t.g(γ,x,y*)≥0 (2)
γ∈{0,1} (3)
Figure FDA0003541342860000058
s.t.A+B(γ)y＝b1:ψ (5)
C+Dy≥b2:ζ (6)
in this model, f (-) and dTy represents the upper and lower optimization model objective functions, gamma and x represent the upper optimization model 0-1 variable and continuous variable, respectively, y represents the lower optimization model variable, psi and zeta represent the lagrangian multipliers corresponding to equations (5) and (6), D, a, B (gamma), C, D, B1And b2Respectively representing a coefficient matrix and a constant item corresponding to the constraint, wherein B (gamma) indicates that the coefficient matrix is related to a parameter gamma; the formula (2) represents the constraint condition of the upper-layer optimization model, and the formulas (5) and (6) represent the equality and inequality constraints of the lower-layer optimization model respectively;
by using the dual theory, the model can be converted into a single-layer optimization model as shown in the following:
Figure FDA0003541342860000061
s.t.g(γ,x,y)≥0
γ∈{0,1}
B(γ)Tψ+DTζ＝d
ζ≥0
A+B(γ)y＝b1
C+Dy≥b2
dTy＝ψT[b1-A]+ζT[b2-C]。
8. a computer storage medium, in which a computer program executable by a processor is stored, the computer program being specifically configured to perform the method for identifying a critical line under an N-k fault of a power system in consideration of an information physical coupling relationship according to any one of claims 1 to 4.
CN202010674779.0A 2020-07-14 2020-07-14 Key line identification method under N-k fault of power system considering information physical coupling relation Active CN112117757B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010674779.0A CN112117757B (en) 2020-07-14 2020-07-14 Key line identification method under N-k fault of power system considering information physical coupling relation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010674779.0A CN112117757B (en) 2020-07-14 2020-07-14 Key line identification method under N-k fault of power system considering information physical coupling relation

Publications (2)

Publication Number Publication Date
CN112117757A CN112117757A (en) 2020-12-22
CN112117757B true CN112117757B (en) 2022-05-20

Family

ID=73799035

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010674779.0A Active CN112117757B (en) 2020-07-14 2020-07-14 Key line identification method under N-k fault of power system considering information physical coupling relation

Country Status (1)

Country Link
CN (1) CN112117757B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113344743B (en) * 2021-07-26 2022-05-17 国网四川省电力公司经济技术研究院 Fault hazard index calculation and vulnerability assessment method for smart power grid
CN113836679B (en) * 2021-10-14 2024-02-23 国网湖南省电力有限公司 Method and device for identifying vulnerable line combination in N-K attack mode
CN114069724B (en) * 2021-12-06 2023-06-20 国网四川省电力公司电力科学研究院 Information link vulnerability analysis method for electric power information physical system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009052448A2 (en) * 2007-10-17 2009-04-23 V2Green, Inc. Electric resource power meter in a power aggregation system for distributed electric resources
CN107317328A (en) * 2017-08-17 2017-11-03 国网江苏省电力公司宿迁供电公司 Power system N K failure analysis methods based on bi-level optimization
CN109801183A (en) * 2018-12-24 2019-05-24 全球能源互联网研究院有限公司 A kind of the fault harm appraisal procedure and device of power information physics system
CN110135727A (en) * 2019-05-14 2019-08-16 华北电力大学 Power network and Information Network Fusion Model based on information physical fusion

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9960637B2 (en) * 2015-07-04 2018-05-01 Sunverge Energy, Inc. Renewable energy integrated storage and generation systems, apparatus, and methods with cloud distributed energy management services

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009052448A2 (en) * 2007-10-17 2009-04-23 V2Green, Inc. Electric resource power meter in a power aggregation system for distributed electric resources
CN107317328A (en) * 2017-08-17 2017-11-03 国网江苏省电力公司宿迁供电公司 Power system N K failure analysis methods based on bi-level optimization
CN109801183A (en) * 2018-12-24 2019-05-24 全球能源互联网研究院有限公司 A kind of the fault harm appraisal procedure and device of power information physics system
CN110135727A (en) * 2019-05-14 2019-08-16 华北电力大学 Power network and Information Network Fusion Model based on information physical fusion

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Cyber-Constrained Optimal Power Flow Model for Smart Grid Resilience Enhancement;Gang Huang等;《IEEE TRANSACTIONS ON SMART GRID》;20190930;全文 *
电力系统信息物理融合建模与综合安全评估：驱动力与研究构想;郭庆来等;《中 国 电 机 工 程 学 报》;20160320;全文 *

Also Published As

Publication number Publication date
CN112117757A (en) 2020-12-22

Similar Documents

Publication Publication Date Title
CN112117757B (en) Key line identification method under N-k fault of power system considering information physical coupling relation
US20110264276A1 (en) Interconnected electrical network and building management system and method of operation
CN110311427B (en) Two-stage N-K robust fault constraint unit combination method considering fault probability
CN107462761B (en) Distribution transformer anti-burning grading early warning method based on multivariable monitoring multi-threshold constraint
CN107658895B (en) Threshold value setting method for direct-current commutation failure prediction link
CN108493998B (en) Robust power transmission network planning method considering demand response and N-1 expected faults
CN105244865A (en) Power system safe and stable operation control method
CN108767848B (en) Electric power system fragile line identification and cascading failure prevention method
CN113141022B (en) Alternating current-direct current series-parallel connection receiving end power grid structure optimization method
CN109802417B (en) Emergency control method and device for power grid for coping with direct current fault impact weak alternating current channel
CN110932291B (en) Active splitting method for preserving stability of local power grid
Bourgin et al. Present practices and trends on the French power system to prevent voltage collapse
CN107067126B (en) Thermal stability key power transmission channel identification method based on power flow transfer ratio
CN109672175B (en) Power grid control method and device
CN106960262B (en) Active reserve capacity evaluation method and device for expected faults of extra-high voltage direct current line
CN110096767B (en) Cascading failure simulation method for alternating current-direct current series-parallel power grid
CN111130116A (en) Scheduling operation power flow checking method based on key topology change item identification
CN112615363B (en) Emergency coordination control method and system for dealing with multi-loop direct current blocking fault
CN111652469B (en) Weak link identification method and system for alternating current-direct current hybrid power grid
CN107317328A (en) Power system N K failure analysis methods based on bi-level optimization
CN114252735A (en) Monitoring and judging method for 110kV line broken branch positioning
CN111160746B (en) Intelligent power distribution network planning evaluation method based on risk assessment
CN110137926B (en) Power limit calculation method and system of power transmission section considering multiple faults
CN113270871A (en) Flexible interconnection device capacity configuration optimization method based on power distribution network N-1 safety assessment
CN104636607A (en) Evaluation method for refinery enterprise electric network static security features based on BPA

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