CN114065536A - Modeling method for layered collaborative recovery model of electricity-gas interconnection system - Google Patents

Modeling method for layered collaborative recovery model of electricity-gas interconnection system Download PDF

Info

Publication number
CN114065536A
CN114065536A CN202111385987.XA CN202111385987A CN114065536A CN 114065536 A CN114065536 A CN 114065536A CN 202111385987 A CN202111385987 A CN 202111385987A CN 114065536 A CN114065536 A CN 114065536A
Authority
CN
China
Prior art keywords
gas
power
constraint
transmission
distribution
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111385987.XA
Other languages
Chinese (zh)
Inventor
周健
叶洪波
凌晓波
杨心刚
郭灵瑜
时珊珊
顾承昱
贠靖洋
周云
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiaotong University
State Grid Shanghai Electric Power Co Ltd
Original Assignee
Shanghai Jiaotong University
State Grid Shanghai 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 Shanghai Jiaotong University, State Grid Shanghai Electric Power Co Ltd filed Critical Shanghai Jiaotong University
Priority to CN202111385987.XA priority Critical patent/CN114065536A/en
Publication of CN114065536A publication Critical patent/CN114065536A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/04—Power grid distribution networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/14—Pipes

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Human Resources & Organizations (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Public Health (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to a modeling method of a layered collaborative recovery model of an electricity-gas interconnection system, which comprises the following steps: establishing a layered collaborative recovery framework model of the electric-gas interconnection system; establishing an electric-gas interconnection system cooperative recovery model; adding power-natural gas network coupling constraints in an electric-gas interconnection system collaborative recovery model, wherein the power-natural gas network coupling constraints comprise power transmission-gas transmission system coupling constraints, power distribution-gas distribution system coupling constraints and boundary connection constraints of a power transmission-gas transmission system and a power distribution-gas distribution system; and finishing the establishment of the electric-gas interconnection system cooperative recovery model. Compared with the prior art, the method has the advantages of high cooperative recovery speed, full consideration of the coupling characteristic of the power-natural gas network and the like.

Description

Modeling method for layered collaborative recovery model of electricity-gas interconnection system
Technical Field
The invention relates to the technical field of power network recovery, in particular to a modeling method of a layered collaborative recovery model of an electricity-gas interconnection system.
Background
The traditional power system recovery means that after a local power failure or a large-area power failure of a system, a gas unit, a hydroelectric unit and other units with black start capability in the system or an external power supply provide starting power for other units without self-start capability in the system, a main network frame and load of the system are gradually recovered, and the system is safely and effectively recovered to a new normal running state.
However, the conventional power system recovery model has the following defects: the recovery of the power transmission network is taken as a core, the self recovery capability of the power distribution network is less considered, and the power distribution network needs to be recovered after the main unit is connected to the grid and the main power transmission network is recovered, so that the continuous power supply of important loads is influenced; secondly, after a large-scale energy supply interruption accident occurs, the recovery of the power and natural gas energy network has a close interdependence characteristic, and the traditional power network black start model cannot consider the influence of the power-natural gas network coupling characteristic on the recovery of the power network, so that the black start efficiency of the power network is low and the speed is slow.
In large-scale power outage, the starting time of black start is the most important consideration index of the black start problem, because a power system cooperative recovery model capable of considering the coupling characteristics of the power-natural gas network is urgently needed.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provide a modeling method of a layered cooperative recovery model of an electric-gas interconnection system, which has high cooperative recovery speed and fully considers the coupling characteristic of an electric power-natural gas network.
The purpose of the invention can be realized by the following technical scheme:
a modeling method of a layered collaborative recovery model of an electric-gas interconnection system comprises the following steps:
step 1: establishing a layered collaborative recovery framework model of the electric-gas interconnection system;
step 2: establishing an electric-gas interconnection system cooperative recovery model based on the framework model in the step 1;
and step 3: adding power-natural gas network coupling constraints in an electric-gas interconnection system collaborative recovery model, wherein the power-natural gas network coupling constraints comprise power transmission-gas transmission system coupling constraints, power distribution-gas distribution system coupling constraints and boundary connection constraints of a power transmission-gas transmission system and a power distribution-gas distribution system;
and 4, step 4: and finishing the establishment of the electric-gas interconnection system cooperative recovery model.
Preferably, the step 1 specifically comprises:
suppose that an energy supply interruption accident occurs at t0Time of day from tkStarting to recover the electric and gas loads on the level of the power transmission-gas transmission system from the moment, and formally starting to recover the load of the power distribution-gas distribution system in a large scale;
hierarchical collaborative recovery from t for an electrical-to-electrical interconnection systemkStarting at the moment, namely the recovery preparation stage is completed, starting each power supply and air source in the power transmission-gas transmission system, and putting the power transmission line, the gas transmission channel, the distribution network transformer substation and the distribution network voltage regulating station into operation;
and establishing an electric-gas interconnection system layered collaborative recovery framework model based on the electric-gas interconnection system at the moment.
Preferably, the electric-gas interconnection system layered collaborative recovery model specifically comprises:
the method comprises the following steps of maximizing the weighted supply quantity of the electric load and the gas load in the power transmission-gas transmission system and the power distribution-gas distribution system as an objective function of a model, and specifically comprises the following steps:
Figure BDA0003367040120000021
wherein, each variable superscript T represents a transmission network or a transmission network, and superscript D represents a distribution network or a distribution network; subscript E represents the power system, subscript G represents the natural gas system; superscript L represents load; the subscript T denotes the respective recovery time, TRA set representing each recovery time; subscript m represents the transmission network node number, subscript n represents the transmission network node number, subscript i represents the distribution network node number, and subscript j represents the distribution network node number;
Figure BDA0003367040120000022
Figure BDA0003367040120000023
and
Figure BDA0003367040120000024
respectively representing the sets of all nodes of a transmission network, a gas transmission network, a distribution network and a distribution network;
Figure BDA0003367040120000025
Figure BDA0003367040120000026
and
Figure BDA0003367040120000027
respectively representing the weight coefficients of each node of the corresponding system;
Figure BDA0003367040120000028
Figure BDA0003367040120000029
the load recovered by each system at each moment is respectively.
Preferably, the electric-gas interconnection system layered collaborative recovery model is provided with a power transmission-gas transmission system recovery constraint and a power distribution-gas distribution system recovery constraint;
the power transmission-gas transmission system recovery constraints comprise power transmission network flow constraints and safety constraints, gas transmission network flow constraints and safety constraints and power transmission-gas transmission system load recovery constraints;
the power distribution-gas distribution system recovery constraints comprise power distribution network flow constraints and safety constraints, gas distribution network flow constraints and safety constraints and power distribution-gas distribution system load recovery constraints.
More preferably, the power transmission network power flow constraint specifically includes:
and (3) carrying out active and reactive balance constraint on nodes of the power transmission system:
Figure BDA0003367040120000031
Figure BDA0003367040120000032
wherein, each variable superscript G represents a power supply; m' is a node adjacent to m, KTThe node m is a set of power transmission lines connected with the node m;
Figure BDA0003367040120000033
active power injection of a node m and active power flowing through a line mm' are respectively carried out;
Figure BDA0003367040120000034
respectively reactive power injection, reactive load and line mm' flowing reactive power of a node m;
equation Right side
Figure BDA0003367040120000035
And
Figure BDA0003367040120000036
respectively representing the active and reactive total losses on the transmission line connected with the node m;
the power transmission system line active and reactive balance constraint:
Figure BDA0003367040120000037
Figure BDA0003367040120000038
wherein the content of the first and second substances,
Figure BDA0003367040120000039
respectively mm' conductance and susceptance of the line;
Figure BDA00033670401200000310
voltage amplitudes of the node m and the node m' respectively;
Figure BDA00033670401200000311
is the phase angle difference of the voltage at the two ends of the line mm'; wherein the voltage amplitude and cosine terms
Figure BDA00033670401200000312
And sine term
Figure BDA00033670401200000313
The cross product term of the voltage is processed by a secondary relaxation technology, and the voltage amplitude is a square term
Figure BDA00033670401200000314
Processing by adopting an improved piecewise linearization method, thereby obtaining a convex expression form of the line power flow constraint and realizing convex relaxation of the constraint;
and (3) recovering safety constraint of the power transmission network:
including voltage amplitude restraint, power active, reactive power upper and lower limits restraint and climbing restraint, do respectively:
Figure BDA00033670401200000315
Figure BDA00033670401200000316
Figure BDA00033670401200000317
Figure BDA00033670401200000318
wherein the content of the first and second substances,
Figure BDA00033670401200000319
and
Figure BDA00033670401200000320
the upper limit and the lower limit of the active power output of the power supply are respectively;
Figure BDA00033670401200000321
and
Figure BDA00033670401200000322
respectively an upper limit and a lower limit of the reactive power output of the power supply,
Figure BDA00033670401200000323
the upper limit of the climbing capacity of the power supply in the adjacent recovery time period;
and (3) constraint between the variable quantity of the total quantity of the active load of the system between the adjacent recovery moments and the maximum total quantity of the active load recovery:
Figure BDA0003367040120000041
more preferably, the power flow constraint and the safety constraint of the gas transmission network are specifically as follows:
and (3) airflow balance constraint of a node n in the gas transmission network system:
Figure BDA0003367040120000042
wherein the content of the first and second substances,
Figure BDA0003367040120000043
and
Figure BDA0003367040120000044
the natural gas flow injected by the source at node n and the natural gas flow consumed by the gas load at time t,
Figure BDA0003367040120000045
the flow rate of the natural gas flowing through the gas transmission pipeline nn';
the relationship between the gas flow of the gas transmission pipeline and the pressure difference of two ends of the pipeline is restricted:
Figure BDA0003367040120000046
wherein the content of the first and second substances,
Figure BDA0003367040120000047
the gas pressures at node n and node n', respectively,
Figure BDA0003367040120000048
is the Weymouth constant for pipe nn';
wherein, the expression of the symbolic function of the air pressure difference between the two ends of the pipeline is as follows:
Figure BDA0003367040120000049
for the nonlinear terms in the above constraints, processing by adopting an improved piecewise linearization method;
and (3) recovering safety constraint of the gas transmission network:
the method comprises the following steps of pipeline airflow upper and lower limit restraint, node air pressure upper and lower limit restraint, air source output upper and lower limit restraint and climbing restraint, which are respectively as follows:
Figure BDA00033670401200000410
Figure BDA00033670401200000411
Figure BDA00033670401200000412
Figure BDA00033670401200000413
more preferably, the load recovery constraint of the power transmission and gas transmission system is specifically:
and active load recovery constraint:
Figure BDA00033670401200000414
reactive load recovery constraint:
Figure BDA00033670401200000415
and (3) upper limit constraint of the load of the power transmission network:
Figure BDA0003367040120000051
and (4) gas source safe operation constraint:
Figure BDA0003367040120000052
Figure BDA0003367040120000053
more preferably, the power flow constraint and the safety constraint of the power distribution network are specifically as follows:
the active and reactive power balance constraint at a point i at a time t in the mth power distribution network:
Figure BDA0003367040120000054
Figure BDA0003367040120000055
wherein the content of the first and second substances,
Figure BDA0003367040120000056
and
Figure BDA0003367040120000057
active power flow and reactive power flow flowing through a line ii' connected with the node i respectively;
Figure BDA0003367040120000058
is the square of the magnitude of the current flowing through line ii';
Figure BDA0003367040120000059
and
Figure BDA00033670401200000510
the active power and reactive power injected for node i respectively,
Figure BDA00033670401200000511
and
Figure BDA00033670401200000512
respectively an active load and a reactive load which are actually recovered at a node i;
the relationship of the square of the voltage amplitude between node i and node i' constrains:
Figure BDA00033670401200000513
wherein r isii′And xii′Resistance and reactance of line ii', respectively;
Figure BDA00033670401200000514
is the square of the line impedance modulus;
and (3) the node voltage amplitude and the corresponding line current amplitude, active power and reactive power current value are subjected to relation constraint:
Figure BDA00033670401200000515
the constraint is a non-convex constraint, for cross product terms
Figure BDA00033670401200000516
Adopting McCormick convex relaxation method to treat, and obtaining a variable square term
Figure BDA00033670401200000517
Processing by adopting a piecewise linearization method;
safety restraint of the power distribution network:
including voltage assignment page number, the restraint of the upper and lower limits of power supply active and reactive power output and climbing restraint, specifically do:
Figure BDA00033670401200000518
Figure BDA00033670401200000519
Figure BDA00033670401200000520
Figure BDA00033670401200000521
more preferably, the load recovery constraint of the distribution-distribution system is specifically:
and active load recovery constraint:
Figure BDA00033670401200000522
reactive load recovery constraint:
Figure BDA0003367040120000061
and (3) terminal power load recovery continuity constraint:
Figure BDA0003367040120000062
the active and reactive power constraints accepted by the distribution substation:
Figure BDA0003367040120000063
Figure BDA0003367040120000064
and (3) recovering and constraining the gas load for the terminal of the gas distribution network:
Figure BDA0003367040120000065
Figure BDA0003367040120000066
the gas distribution pressure regulating station receives natural gas flow constraint:
Figure BDA0003367040120000067
preferably, the power transmission-gas transmission system coupling constraint and the power distribution-gas distribution system coupling constraint in step 3 specifically include:
and (3) energy conversion relation constraint of the gas turbine set:
Figure BDA0003367040120000068
wherein the content of the first and second substances,
Figure BDA0003367040120000069
the natural gas consumption of a gas unit connected between a power distribution network node i and a gas distribution network node j at a time t;
Figure BDA00033670401200000610
the active power output of the gas unit;
Figure BDA00033670401200000611
to the efficiency of the power generation; g is the low heating value of natural gas;
energy conversion relation constraint of a gas turbine unit connected with the power transmission network and the gas transmission network;
P2G facility energy conversion relationship constraints:
Figure BDA00033670401200000612
wherein the content of the first and second substances,
Figure BDA00033670401200000613
as electric load for P2G facilityActive power consumed;
Figure BDA00033670401200000614
for the purpose of a corresponding natural gas flow rate,
Figure BDA00033670401200000615
for gas production efficiency;
the relation constraint of the gas flow rate of the natural gas transmitted by the power-driven pressure regulating station of the gas transmission network and the consumed electric power is as follows:
Figure BDA00033670401200000616
wherein the content of the first and second substances,
Figure BDA00033670401200000617
the natural gas flow which flows through the electric drive pressure regulating station between the gas transmission network node n and the gas distribution network node j;
Figure BDA00033670401200000618
electric power consumed as an electric load for the electrically driven regulator station; beta is areg,TThe energy consumption coefficient;
the boundary connection constraint of the power distribution-gas distribution system specifically comprises the following steps:
Figure BDA0003367040120000071
Figure BDA0003367040120000072
Figure BDA0003367040120000073
the above constraints represent the power and voltage correlations of the transmission system and the transmission system at the substation, respectively, where Vt D,m,sub,sqrIndicating connection to node m of the transmission networkThe square of the voltage amplitude of the substation;
the boundary connection constraint of the gas transmission-distribution system is as follows:
Figure BDA0003367040120000074
Figure BDA0003367040120000075
Figure BDA0003367040120000076
and
Figure BDA0003367040120000077
for the boundary connection constraint of the gas transmission-distribution system,
Figure BDA0003367040120000078
indicating the pressure at the boost side of the distribution pressure regulating station connected to the gas transmission network node n,
Figure BDA0003367040120000079
representing the pressure on the pressure-reducing side, gammacIs the voltage boosting ratio;
will be provided with
Figure BDA00033670401200000710
Figure BDA00033670401200000711
And
Figure BDA00033670401200000712
and setting boundary connection variables of the power transmission-gas transmission system and the power distribution-gas distribution system in the recovery process.
Compared with the prior art, the invention has the following beneficial effects:
firstly, the cooperative recovery speed is faster: the modeling method of the layered collaborative restoration model of the electric-gas interconnection system considers the influence of the coupling characteristic of the electric power-natural gas network on the restoration of the electric power network, adds the force-natural gas network coupling constraint in the model, and fully utilizes the close interdependence characteristic of the restoration of the electric power and natural gas energy networks to ensure that the restoration speed is higher when the model is used for carrying out the collaborative restoration of the electric power system, wherein the force-natural gas network coupling constraint comprises the coupling constraint of a power transmission-gas transmission system, the coupling constraint of a power distribution-gas distribution system and the boundary connection constraint of the power transmission-gas transmission system and the power distribution-gas distribution system.
And secondly, fully considering the coupling characteristic of the power-natural gas network: according to the modeling method of the layered collaborative recovery model of the electric-gas interconnection system, the layered collaborative recovery frame model of the electric-gas interconnection system is provided before formal modeling, the influence of the coupling characteristic of the electric power-natural gas network on the recovery of the electric power network is fully considered, modeling is carried out according to the coupling characteristic of the electric power-natural gas network, and the reliability of the model is high.
Drawings
FIG. 1 is a schematic flow chart of a modeling method of a layered cooperative recovery model of an electrical-gas interconnection system according to the present invention;
FIG. 2 is a schematic diagram of an overall framework mechanism for hierarchical cooperative recovery of an electrical-to-electrical interconnection system in an embodiment of the present invention;
FIG. 3 is a schematic diagram of the interaction between the power system and the natural gas system in the embodiment of the invention;
fig. 4 is a schematic diagram of energy conversion and information interaction in the collaborative recovery process of the model in the embodiment of the present invention.
Detailed Description
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 some, not all, embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, shall fall within the scope of protection of the present invention.
The embodiment relates to a modeling method of a layered collaborative recovery model of an electric-gas interconnection system, the flow of which is shown in fig. 1, and the modeling method comprises the following steps:
step 1: establishing a layered collaborative recovery framework model of the electric-gas interconnection system;
step 2: establishing an electric-gas interconnection system cooperative recovery model based on the framework model in the step 1;
and step 3: adding power-natural gas network coupling constraints in an electric-gas interconnection system collaborative recovery model, wherein the power-natural gas network coupling constraints comprise power transmission-gas transmission system coupling constraints, power distribution-gas distribution system coupling constraints and boundary connection constraints of a power transmission-gas transmission system and a power distribution-gas distribution system;
the following describes each step in detail:
firstly, establishing an electric-gas interconnection system layered cooperative recovery framework
The electrical-to-gas interconnection system recovery problem can be decomposed into a power transmission-gas transmission system recovery problem and a power distribution-gas distribution system recovery problem. The main tasks of the recovery of the transmission-transmission system mainly comprise: 1) starting the energy source (black starting of the generator set and energy recovery of the air source); 2) recovering the power transmission net rack and the gas transmission net rack; 3) and (4) recovering the electric load of the distribution network transformer substation and the gas load of the distribution network transformer substation. The restoration of a distribution-distribution system usually begins with a load restoration phase of a transmission-distribution system, and the tasks mainly include: 1) optimizing the input sequence and the size of the load according to the importance degree of the electric and gas load, and 2) determining the scheduling decision of a distributed power supply and a gas source (comprising a distributed gas storage facility and a P2G device) in the power distribution-gas distribution system.
The overall framework mechanism for hierarchical cooperative recovery of an electro-pneumatic interconnection system is shown in fig. 2. Suppose that a large-scale power supply interruption accident occurs at t0Time of day from t1Starting from the moment, the recovery of the electric-gas interconnection system enters a recovery preparation stage, the power transmission-gas transmission system starts the recovery of the energy source and the recovery of the power transmission and gas transmission net rack, at the moment, the net rack and the load nodes of the power distribution-gas distribution system start to recover, partial electric energy from the power transmission-gas transmission system can be received to start an internal power supply and a P2G device and recover partial terminal loads, and the unit output in the power transmission-gas transmission system is balanced; if the power distribution-distribution system has an energy source which is started quickly, the power transmission-gas transmission system can be injected with power or natural gas to start the energy source. From tkAt the beginning of timeThe electric and gas loads on the level of the power transmission-gas transmission system are recovered (the nodes of the transformer substation and the voltage regulating station of the power distribution-gas distribution system are the electric and gas load nodes), and the loads of the power distribution-gas distribution system are formally recovered in a large scale. Hierarchical cooperative recovery slave t of electric-gas interconnection system in this embodimentkStarting from moment, namely the recovery preparation stage is completed, starting each power supply and air source in the power transmission-gas transmission system, putting the power transmission line, the gas transmission channel, the distribution network transformer substation and the distribution network voltage regulating station into operation, and performing important research tkAnd the physical process and mechanism model of information interaction and mutual support of the power transmission-gas transmission system and the power distribution-gas distribution system in the load recovery process after the moment.
The framework mechanism for layered cooperative recovery of the electric-gas interconnection system specifically comprises the following implementation steps:
(1) and solving the load recovery problem of the power transmission-gas transmission system.
The output of each power supply and each air source in the power transmission-gas transmission system at each moment, the node voltage and the air pressure of each power distribution transformer station and each gas distribution pressure regulating station, and the load recovery supply quantity are determined.
(2) And solving the load recovery problem of the power distribution-gas distribution system.
From tkAnd determining the load recovery sequence and size in the power distribution-gas distribution system, the scheduling decision of each distributed power supply and gas source, and boundary connection variables such as voltage, air pressure, active and reactive power flow, gas flow and the like of the power distribution transformer station and the gas distribution and pressure regulation station.
(3) And (4) information interaction and updating.
The power transmission-gas transmission system load recovery problem and the power distribution-gas distribution load system recovery problem which are solved in the two steps are at tkAnd (4) interacting the information of the corresponding boundary contact variables of the power transmission-gas transmission system and the power distribution-gas distribution system and iteratively updating the corresponding Lagrange multipliers according to the results of all moments after the moment.
In summary, through a layered cooperative recovery mechanism of an electric-gas interconnection system, the solution of the original centralized electric-gas interconnection system recovery problem can be decomposed into the solution of a load recovery subproblem of a power transmission-gas transmission system and a plurality of load recovery subproblems of a power distribution-gas distribution system, and the interaction of limited boundary link variable information is performed after each solution, so as to ensure the global optimality of the solution.
Establishing an electric-gas interconnection system cooperative recovery model
Defining an objective function to maximize the weighted supply of electrical and gas loads within the transmission and distribution systems:
Figure BDA0003367040120000091
wherein, each variable superscript T represents a transmission network or a transmission network (transmission layer), and a superscript D represents a distribution network or a distribution network (distribution layer); subscript E represents the power system, subscript G represents the natural gas system; the superscript L represents the load. The subscript T denotes the respective recovery time, TRA set representing each recovery time; subscript m represents the transmission network node number, subscript n represents the transmission network node number, subscript i represents the distribution network node number, and subscript j represents the distribution network node number. In this embodiment, it is assumed that each load node of the transmission network is connected to the distribution network, and is not connected to the terminal power load, and the transmission network is the same.
Figure BDA0003367040120000101
And
Figure BDA0003367040120000102
respectively representing the sets of all nodes of a transmission network, a gas transmission network, a distribution network and a distribution network;
Figure BDA0003367040120000103
and
Figure BDA0003367040120000104
respectively representing the weight coefficients of the nodes of the corresponding system.
Figure BDA0003367040120000105
The load recovered by each system at each moment is respectively.
It should be noted that, in the above objective function, both the electrical load and the natural gas load are normalized, that is, in the objective function, both the recovered electrical load and the recovered gas load are percentages of the recovered load in the original total load. The present embodiment, for convenience of description, will still write the objective function in the above form.
The layered collaborative recovery model of the electric-gas interconnection system further comprises power transmission network flow constraint and safety constraint and power distribution-gas distribution system recovery constraint, which are respectively described as follows:
1. transmission-transmission system recovery constraints
(1) Power transmission network flow and safety constraints
The problem of recovery of the power system needs to consider reactive and voltage problems, so that an alternating current power flow model of the power transmission system needs to be established. The duration of the electromechanical transient process of the power system is in the order of seconds, and the duration of the electromagnetic transient process is in the order of milliseconds or even microseconds. The power transmission network power flow constraint can be established based on an algebraic power flow equation of a steady-state power system.
And (3) carrying out active and reactive balance constraint on nodes of the power transmission system:
Figure BDA0003367040120000106
Figure BDA0003367040120000107
wherein, each variable superscript G represents the power supply. m' is a node adjacent to m, KTIs the collection of the transmission lines connected with the node m.
Figure BDA0003367040120000108
Active power injection of a node m and active power flowing through a line mm' are respectively carried out;
Figure BDA0003367040120000109
respectively the reactive power injection of the node m, the reactive load and the reactive power flowing through the line mm'. Equation Right side
Figure BDA00033670401200001010
And
Figure BDA00033670401200001011
respectively representing the total active and reactive losses on the transmission line to which node m is connected.
The power transmission system line active and reactive balance constraint:
Figure BDA0003367040120000111
Figure BDA0003367040120000112
wherein the content of the first and second substances,
Figure BDA0003367040120000113
respectively mm' conductance and susceptance of the line;
Figure BDA0003367040120000114
voltage amplitudes of the node m and the node m' respectively;
Figure BDA0003367040120000115
is the phase angle difference of the voltage across the mm' line. Wherein the voltage amplitude and cosine terms
Figure BDA0003367040120000116
And sine term
Figure BDA0003367040120000117
The cross product term of the voltage is processed by a secondary relaxation technology, and the voltage amplitude is a square term
Figure BDA0003367040120000118
And processing by adopting an improved piecewise linearization method, thereby obtaining a convex expression form of the line power flow constraint and realizing convex relaxation of the constraint.
And (3) recovering safety constraint of the power transmission network:
including voltage amplitude restraint, power active, reactive power upper and lower limits restraint and climbing restraint, do respectively:
Figure BDA0003367040120000119
Figure BDA00033670401200001110
Figure BDA00033670401200001111
Figure BDA00033670401200001112
wherein the content of the first and second substances,
Figure BDA00033670401200001113
and
Figure BDA00033670401200001114
the upper limit and the lower limit of the active power output of the power supply are respectively;
Figure BDA00033670401200001115
and
Figure BDA00033670401200001116
respectively an upper limit and a lower limit of the reactive power output of the power supply,
Figure BDA00033670401200001117
the upper limit of the climbing capacity of the power supply in the adjacent recovery period.
Meanwhile, in order to ensure the frequency stability of the power transmission system in the recovery process, the variation of the total active load recovery amount of the system between adjacent recovery moments is lower than the maximum total active load recovery amount:
Figure BDA00033670401200001118
(2) power flow and safety constraints for gas transmission network
In the embodiment, a physical process and a mechanism model for information interaction and mutual support of a power transmission-gas transmission system and a power distribution-gas distribution system in a recovery process are mainly researched, and a steady-state Weymouth flow equation with neglected dynamic characteristics is adopted for flow constraint modeling of a gas transmission network and a gas distribution network.
And (3) airflow balance constraint of a node n in the gas transmission network system:
Figure BDA00033670401200001119
wherein the content of the first and second substances,
Figure BDA00033670401200001120
and
Figure BDA00033670401200001121
the natural gas flow injected by the source at node n and the natural gas flow consumed by the gas load at time t,
Figure BDA00033670401200001122
the flow of natural gas through the gas transmission line nn'.
The relationship between the gas flow of the gas transmission pipeline and the pressure difference of two ends of the pipeline is restricted:
Figure BDA0003367040120000121
wherein the content of the first and second substances,
Figure BDA0003367040120000122
the gas pressures at node n and node n', respectively,
Figure BDA0003367040120000123
is the Weymouth constant for the pipe nn'.
Wherein, the expression of the symbolic function of the air pressure difference between the two ends of the pipeline is as follows:
Figure BDA0003367040120000124
for the non-linear terms in the above constraints, an improved piecewise linearization method is adopted for processing.
And (3) recovering safety constraint of the gas transmission network:
the method comprises the following steps of pipeline airflow upper and lower limit restraint, node air pressure upper and lower limit restraint, air source output upper and lower limit restraint and climbing restraint, which are respectively as follows:
Figure BDA0003367040120000125
Figure BDA0003367040120000126
Figure BDA0003367040120000127
Figure BDA0003367040120000128
(3) transmission-transmission system load recovery constraints
The method comprises the following steps of active load recovery constraint, reactive load recovery constraint, transmission grid load upper limit constraint and gas source safe operation constraint of a gas transmission system, wherein the active load recovery constraint, the reactive load recovery constraint, the transmission grid load upper limit constraint and the gas source safe operation constraint are respectively as follows:
Figure BDA0003367040120000129
Figure BDA00033670401200001210
wherein the content of the first and second substances,
Figure BDA00033670401200001211
and
Figure BDA00033670401200001212
respectively the maximum active and reactive loads to be restored at node m. The recovery of the active load of the transmission network should be limited to the maximum load of its nodes, while the reactive load should be recovered in proportion to the active load according to a specific power factor.
Since this embodiment assumes that the transmission network load node is a distribution network substation node, the transmission network load upper limit is the distribution network substation capacity:
Figure BDA00033670401200001213
similar to the power transmission system, the gas source safe operation constraint of the power transmission system is as follows:
Figure BDA00033670401200001214
Figure BDA00033670401200001215
2. distribution-distribution system recovery constraint
(1) Power flow constraint and safety constraint of power distribution network
The influence of unbalanced three-phase load of the power distribution network on recovery is ignored in the embodiment, and the power distribution network is assumed to be a three-phase balance system. And modeling the power distribution network Flow constraint by adopting a single-phase Closed-form Branch Flow Model. For the mth distribution network (the number of the distribution network is the number of the load node of the transmission network, and the number superscript m of the distribution network is omitted here for the sake of simplicity), the mathematical expression form of the power flow constraint is as follows:
the active and reactive power balance constraint at a point i at a time t in the mth power distribution network:
Figure BDA0003367040120000131
Figure BDA0003367040120000132
wherein the content of the first and second substances,
Figure BDA0003367040120000133
and
Figure BDA0003367040120000134
active power flow and reactive power flow flowing through a line ii' connected with the node i respectively;
Figure BDA0003367040120000135
is the square of the magnitude of the current flowing through line ii'. It is worth noting that since the current magnitude in the distribution network flow constraint occurs only in the form of a square term, it will be
Figure BDA0003367040120000136
Treated as independent variables rather than secondary variables.
Figure BDA0003367040120000137
And
Figure BDA0003367040120000138
the active power and reactive power injected for node i respectively,
Figure BDA0003367040120000139
and
Figure BDA00033670401200001310
the active load and the reactive load actually restored at the node i are respectively.
The relationship of the square of the voltage amplitude between node i and node i' constrains:
Figure BDA00033670401200001311
wherein r isii′And xii′Resistance and reactance of line ii', respectively;
Figure BDA00033670401200001312
is the square of the module value of the line impedance. And
Figure BDA00033670401200001313
similarly, the voltage amplitude appears only in the form of a squared term, so it is considered as an independent variable.
And (3) the node voltage amplitude and the corresponding line current amplitude, active power and reactive power current value are subjected to relation constraint:
Figure BDA00033670401200001314
the constraint is a non-convex constraint, for cross product terms
Figure BDA00033670401200001315
Adopting McCormick convex relaxation method to treat, and obtaining a variable square term
Figure BDA00033670401200001316
And (4) processing by adopting a piecewise linearization method.
Safety restraint of the power distribution network:
including voltage assignment page number, the restraint of the upper and lower limits of power supply active and reactive power output and climbing restraint, specifically do:
Figure BDA00033670401200001317
Figure BDA00033670401200001318
Figure BDA00033670401200001319
Figure BDA00033670401200001320
(2) power flow constraint and safety constraint of gas distribution network
The power flow constraint of the gas distribution network is the same as the mathematical expression form of the safety constraint of the gas transmission network.
(3) Load recovery constraints
The load recovery constraints of the distribution network are similar to those of the transmission network:
Figure BDA0003367040120000141
Figure BDA0003367040120000142
to ensure the recovery effect, the recovery of the electrical load of the terminal should have continuity, that is, the load recovered at the subsequent time must not be less than the previous time:
Figure BDA0003367040120000143
in addition, the distribution substation should not accept any more active/reactive power than the power flowing out of its connected grid nodes:
Figure BDA0003367040120000144
Figure BDA0003367040120000145
similarly, the terminal air load recovery constraint of the air distribution network is as follows:
Figure BDA0003367040120000146
Figure BDA0003367040120000147
the natural gas flow rate accepted by the distribution pressure regulating station should also satisfy the following constraints:
Figure BDA0003367040120000148
third, coupling constraint of power transmission-gas transmission system and coupling constraint of power distribution-gas distribution system
Including power transmission-gas transmission system coupling constraints, power distribution-gas distribution system coupling constraints, and boundary connection constraints for power transmission-gas transmission systems and power distribution-gas distribution systems.
Fig. 3 is a schematic diagram of the interaction between the power system and the natural gas system. It can be seen that the operation of the power system requires the fuel supply of the gas system to the gas turbine unit and the micro gas turbine; pipeline transmission of natural gas in natural gas systems requires the supply of electrical energy from the power system to electrically driven pressure regulating stations. Therefore, the gas turbine set in the power transmission system and the electric drive pressure regulating station in the gas transmission system are main coupling points of the power transmission-gas transmission system, and the micro gas turbine, the electric gas conversion facility and the gas distribution system pressure reduction pressure regulating station in the power distribution system are main coupling points of the power distribution-gas distribution system.
The power transmission-gas transmission system coupling constraint and the power distribution-gas distribution system coupling constraint comprise:
and (3) energy conversion relation constraint of the gas turbine set:
Figure BDA0003367040120000149
wherein the content of the first and second substances,
Figure BDA0003367040120000151
the natural gas consumption of a gas unit connected between a power distribution network node i and a gas distribution network node j at a time t;
Figure BDA0003367040120000152
the active power output of the gas unit;
Figure BDA0003367040120000153
to the efficiency of the power generation; g is the low heat value of natural gas, and 35590kJ/m is taken3。
And the energy conversion relation constraint of the gas turbine set connected with the power transmission network and the gas transmission network is the same as the energy conversion relation constraint of the gas turbine set in expression, only the subscript on the variable symbol is different, and the constraint is omitted.
P2G facility energy conversion relationship constraints:
Figure BDA0003367040120000154
wherein the content of the first and second substances,
Figure BDA0003367040120000155
active power consumed as a load for the P2G facility;
Figure BDA0003367040120000156
for the purpose of a corresponding natural gas flow rate,
Figure BDA0003367040120000157
for gas production efficiency;
after an extreme event occurs, if the electric drive pressure regulating station loses the power supply, the residual natural gas drive pressure regulating station in the available state can maintain the basic operation of the network, but the transmission capacity of the gas transmission network is seriously influenced. The goal of restoration of the gas delivery network is therefore to restore the delivery capacity of the network as soon as possible, restoring the supply capacity to the distribution network.
The relation constraint of the gas flow rate of the natural gas transmitted by the power-driven pressure regulating station of the gas transmission network and the consumed electric power is as follows:
Figure BDA0003367040120000158
wherein the content of the first and second substances,
Figure BDA0003367040120000159
the natural gas flow which flows through the electric drive pressure regulating station between the gas transmission network node n and the gas distribution network node j;
Figure BDA00033670401200001510
electric power consumed as an electric load for the electrically driven regulator station; beta is areg,TIs the energy consumption coefficient.
The boundary connection constraint of the power transmission-gas transmission system and the power distribution-gas distribution system is as follows:
the transformer substation and the gas distribution and pressure regulation station are key hubs respectively connected with a power transmission system, a power distribution system, a gas transmission system and a gas distribution system. The transformer substation can be used as a load for a transmission network and a power supply for a power distribution network, and the distribution and voltage regulation station is the same. The boundary connection constraints of the power transmission-gas transmission system and the power distribution-gas distribution system are defined as follows:
Figure BDA00033670401200001511
Figure BDA00033670401200001512
Figure BDA00033670401200001513
the above constraints represent the power sum of the power transmission system and the gas transmission system at the substationA voltage relation relationship in which Vt D,m,sub,sqrRepresents the square of the voltage amplitude of the substation connected to the transmission network node m; it should be noted that all voltage values in this embodiment are per unit.
Meanwhile, considering that the boosting ratio of the gas distribution pressure regulating station is fixed and the natural gas flow at two sides is the same, the boundary relation constraint of the gas transmission-gas distribution system can be obtained:
Figure BDA0003367040120000161
Figure BDA0003367040120000162
Figure BDA0003367040120000163
and
Figure BDA0003367040120000164
for the boundary connection constraint of the gas transmission-distribution system,
Figure BDA0003367040120000165
indicating the pressure at the boost side of the distribution pressure regulating station connected to the gas transmission network node n,
Figure BDA0003367040120000166
representing the pressure on the pressure-reducing side, gammacIs the boost ratio.
This embodiment will be described
Figure BDA0003367040120000167
Figure BDA0003367040120000168
And
Figure BDA0003367040120000169
and setting boundary connection variables of the power transmission-gas transmission system and the power distribution-gas distribution system in the recovery process.
So far, the establishment of the power transmission-gas transmission system and the power distribution-gas distribution system in the cooperative recovery model can be completed, and the energy conversion and information interaction relationship in the process of performing the cooperative recovery on the model is summarized as shown in fig. 4.
While the invention has been described with reference to specific embodiments, the invention is not limited thereto, and various equivalent modifications and substitutions can be easily made by those skilled in the art within the technical scope of the invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (10)

1. A modeling method of a layered collaborative recovery model of an electric-gas interconnection system is characterized by comprising the following steps:
step 1: establishing a layered collaborative recovery framework model of the electric-gas interconnection system;
step 2: establishing an electric-gas interconnection system cooperative recovery model based on the framework model in the step 1;
and step 3: adding power-natural gas network coupling constraints in an electric-gas interconnection system collaborative recovery model, wherein the power-natural gas network coupling constraints comprise power transmission-gas transmission system coupling constraints, power distribution-gas distribution system coupling constraints and boundary connection constraints of a power transmission-gas transmission system and a power distribution-gas distribution system;
and 4, step 4: and finishing the establishment of the electric-gas interconnection system cooperative recovery model.
2. The modeling method of the layered collaborative restoration model of the electrical-electrical interconnection system according to claim 1, wherein the step 1 specifically comprises:
suppose that an energy supply interruption accident occurs at t0Time of day from tkStarting to recover the electric and gas loads on the level of the power transmission-gas transmission system from the moment, and formally starting to recover the load of the power distribution-gas distribution system in a large scale;
hierarchical collaborative recovery from t for an electrical-to-electrical interconnection systemkThe moment begins, namely the recovery preparation stage is completed, and the output is carried outStarting power supplies and gas sources in the electricity-gas transmission system, and putting the power transmission line, the gas transmission channel, the distribution network transformer substation and the distribution network voltage regulation station into operation;
and establishing an electric-gas interconnection system layered collaborative recovery framework model based on the electric-gas interconnection system at the moment.
3. The modeling method of the layered collaborative restoration model of the electrical-electrical interconnection system according to claim 1, wherein the layered collaborative restoration model of the electrical-electrical interconnection system is specifically:
the method comprises the following steps of maximizing the weighted supply quantity of the electric load and the gas load in the power transmission-gas transmission system and the power distribution-gas distribution system as an objective function of a model, and specifically comprises the following steps:
Figure FDA0003367040110000011
wherein, each variable superscript T represents a transmission network or a transmission network, and superscript D represents a distribution network or a distribution network; subscript E represents the power system, subscript G represents the natural gas system; superscript L represents load; the subscript T denotes the respective recovery time, TRA set representing each recovery time; subscript m represents the transmission network node number, subscript n represents the transmission network node number, subscript i represents the distribution network node number, and subscript j represents the distribution network node number;
Figure FDA0003367040110000012
Figure FDA0003367040110000013
and
Figure FDA0003367040110000014
respectively representing the sets of all nodes of a transmission network, a gas transmission network, a distribution network and a distribution network;
Figure FDA0003367040110000015
Figure FDA0003367040110000021
and
Figure FDA0003367040110000022
respectively representing the weight coefficients of each node of the corresponding system;
Figure FDA0003367040110000023
Figure FDA0003367040110000024
the load recovered by each system at each moment is respectively.
4. The modeling method of the layered collaborative restoration model of the electric-gas interconnection system according to claim 1, wherein the layered collaborative restoration model of the electric-gas interconnection system is provided with a restoration constraint of a power transmission-gas transmission system and a restoration constraint of a power distribution-gas distribution system;
the power transmission-gas transmission system recovery constraints comprise power transmission network flow constraints and safety constraints, gas transmission network flow constraints and safety constraints and power transmission-gas transmission system load recovery constraints;
the power distribution-gas distribution system recovery constraints comprise power distribution network flow constraints and safety constraints, gas distribution network flow constraints and safety constraints and power distribution-gas distribution system load recovery constraints.
5. The modeling method of the layered collaborative recovery model of the electrical-electrical interconnection system according to claim 4, wherein the power transmission network power flow constraint is specifically:
and (3) carrying out active and reactive balance constraint on nodes of the power transmission system:
Figure FDA0003367040110000025
Figure FDA0003367040110000026
wherein, each variable superscript G represents a power supply; m' is a node adjacent to m, KTThe node m is a set of power transmission lines connected with the node m;
Figure FDA0003367040110000027
active power injection of a node m and active power flowing through a line mm' are respectively carried out;
Figure FDA0003367040110000028
respectively reactive power injection, reactive load and line mm' flowing reactive power of a node m;
equation Right side
Figure FDA0003367040110000029
And
Figure FDA00033670401100000210
respectively representing the active and reactive total losses on the transmission line connected with the node m;
the power transmission system line active and reactive balance constraint:
Figure FDA00033670401100000211
Figure FDA00033670401100000212
wherein the content of the first and second substances,
Figure FDA00033670401100000213
respectively mm' conductance and susceptance of the line;
Figure FDA00033670401100000214
voltage amplitudes of the node m and the node m' respectively;
Figure FDA00033670401100000215
is the phase angle difference of the voltage at the two ends of the line mm'; wherein the voltage amplitude and cosine terms
Figure FDA00033670401100000216
And sine term
Figure FDA00033670401100000217
The cross product term of the voltage is processed by a secondary relaxation technology, and the voltage amplitude is a square term
Figure FDA00033670401100000218
Processing by adopting an improved piecewise linearization method, thereby obtaining a convex expression form of the line power flow constraint and realizing convex relaxation of the constraint;
and (3) recovering safety constraint of the power transmission network:
including voltage amplitude restraint, power active, reactive power upper and lower limits restraint and climbing restraint, do respectively:
Figure FDA0003367040110000031
Figure FDA0003367040110000032
Figure FDA0003367040110000033
Figure FDA0003367040110000034
wherein the content of the first and second substances,
Figure FDA0003367040110000035
and
Figure FDA0003367040110000036
the upper limit and the lower limit of the active power output of the power supply are respectively;
Figure FDA0003367040110000037
and
Figure FDA0003367040110000038
respectively an upper limit and a lower limit of the reactive power output of the power supply,
Figure FDA0003367040110000039
the upper limit of the climbing capacity of the power supply in the adjacent recovery time period;
and (3) constraint between the variable quantity of the total quantity of the active load of the system between the adjacent recovery moments and the maximum total quantity of the active load recovery:
Figure FDA00033670401100000310
6. the modeling method of the layered collaborative recovery model of the electrical-electrical interconnection system according to claim 4, wherein the power flow constraint and the safety constraint of the power transmission network are specifically as follows:
and (3) airflow balance constraint of a node n in the gas transmission network system:
Figure FDA00033670401100000311
wherein the content of the first and second substances,
Figure FDA00033670401100000312
and
Figure FDA00033670401100000313
the flow rate of the natural gas injected into the source at node n at time t andthe amount of natural gas flow consumed by the gas load,
Figure FDA00033670401100000314
the flow rate of the natural gas flowing through the gas transmission pipeline nn';
the relationship between the gas flow of the gas transmission pipeline and the pressure difference of two ends of the pipeline is restricted:
Figure FDA00033670401100000315
wherein the content of the first and second substances,
Figure FDA00033670401100000316
the gas pressures at node n and node n', respectively,
Figure FDA00033670401100000317
is the Weymouth constant for pipe nn';
wherein, the expression of the symbolic function of the air pressure difference between the two ends of the pipeline is as follows:
Figure FDA00033670401100000318
for the nonlinear terms in the above constraints, processing by adopting an improved piecewise linearization method;
and (3) recovering safety constraint of the gas transmission network:
the method comprises the following steps of pipeline airflow upper and lower limit restraint, node air pressure upper and lower limit restraint, air source output upper and lower limit restraint and climbing restraint, which are respectively as follows:
Figure FDA0003367040110000041
Figure FDA0003367040110000042
Figure FDA0003367040110000043
Figure FDA0003367040110000044
7. the modeling method of the layered cooperative recovery model of the electrical-gas interconnection system according to claim 4, wherein the load recovery constraint of the electrical transmission-gas transmission system is specifically as follows:
and active load recovery constraint:
Figure FDA0003367040110000045
reactive load recovery constraint:
Figure FDA0003367040110000046
and (3) upper limit constraint of the load of the power transmission network:
Figure FDA0003367040110000047
and (4) gas source safe operation constraint:
Figure FDA0003367040110000048
Figure FDA0003367040110000049
8. the modeling method of the layered collaborative recovery model of the electrical-electrical interconnection system according to claim 4, wherein the power flow constraint and the safety constraint of the power distribution network are specifically as follows:
the active and reactive power balance constraint at a point i at a time t in the mth power distribution network:
Figure FDA00033670401100000410
Figure FDA00033670401100000411
wherein the content of the first and second substances,
Figure FDA00033670401100000412
and
Figure FDA00033670401100000413
active power flow and reactive power flow flowing through a line ii' connected with the node i respectively;
Figure FDA00033670401100000414
is the square of the magnitude of the current flowing through line ii';
Figure FDA00033670401100000415
and
Figure FDA00033670401100000416
the active power and reactive power injected for node i respectively,
Figure FDA00033670401100000417
and
Figure FDA00033670401100000418
respectively an active load and a reactive load which are actually recovered at a node i;
the relationship of the square of the voltage amplitude between node i and node i' constrains:
Figure FDA00033670401100000419
wherein r isii′And xii′Resistance and reactance of line ii', respectively;
Figure FDA00033670401100000420
is the square of the line impedance modulus;
and (3) the node voltage amplitude and the corresponding line current amplitude, active power and reactive power current value are subjected to relation constraint:
Figure FDA00033670401100000421
the constraint is a non-convex constraint, for cross product terms
Figure FDA00033670401100000422
Adopting McCormick convex relaxation method to treat, and obtaining a variable square term
Figure FDA0003367040110000051
Processing by adopting a piecewise linearization method;
safety restraint of the power distribution network:
including voltage assignment page number, the restraint of the upper and lower limits of power supply active and reactive power output and climbing restraint, specifically do:
Figure FDA0003367040110000052
Figure FDA0003367040110000053
Figure FDA0003367040110000054
Figure FDA0003367040110000055
9. the modeling method of the layered cooperative recovery model of the electrical-electrical interconnection system according to claim 4, wherein the load recovery constraint of the distribution-distribution system is specifically as follows:
and active load recovery constraint:
Figure FDA0003367040110000056
reactive load recovery constraint:
Figure FDA0003367040110000057
and (3) terminal power load recovery continuity constraint:
Figure FDA0003367040110000058
the active and reactive power constraints accepted by the distribution substation:
Figure FDA0003367040110000059
Figure FDA00033670401100000510
and (3) recovering and constraining the gas load for the terminal of the gas distribution network:
Figure FDA00033670401100000511
Figure FDA00033670401100000512
the gas distribution pressure regulating station receives natural gas flow constraint:
Figure FDA00033670401100000513
10. the modeling method of the layered collaborative restoration model of the electrical-electrical interconnection system according to claim 1, wherein the coupling constraint of the transmission-transmission system and the coupling constraint of the distribution-distribution system in step 3 are specifically:
and (3) energy conversion relation constraint of the gas turbine set:
Figure FDA0003367040110000061
wherein the content of the first and second substances,
Figure FDA0003367040110000062
the natural gas consumption of a gas unit connected between a power distribution network node i and a gas distribution network node j at a time t;
Figure FDA0003367040110000063
the active power output of the gas unit;
Figure FDA0003367040110000064
to the efficiency of the power generation; g is the low heating value of natural gas;
energy conversion relation constraint of a gas turbine unit connected with the power transmission network and the gas transmission network;
P2G facility energy conversion relationship constraints:
Figure FDA0003367040110000065
wherein the content of the first and second substances,
Figure FDA0003367040110000066
active power consumed as a load for the P2G facility;
Figure FDA0003367040110000067
for the purpose of a corresponding natural gas flow rate,
Figure FDA0003367040110000068
for gas production efficiency;
the relation constraint of the gas flow rate of the natural gas transmitted by the power-driven pressure regulating station of the gas transmission network and the consumed electric power is as follows:
Figure FDA0003367040110000069
wherein the content of the first and second substances,
Figure FDA00033670401100000610
the natural gas flow which flows through the electric drive pressure regulating station between the gas transmission network node n and the gas distribution network node j;
Figure FDA00033670401100000611
electric power consumed as an electric load for the electrically driven regulator station; beta is areg,TThe energy consumption coefficient;
the boundary connection constraint of the power distribution-gas distribution system specifically comprises the following steps:
Figure FDA00033670401100000612
Figure FDA00033670401100000613
Figure FDA00033670401100000614
the above constraints represent the power and voltage correlations of the transmission system and the transmission system at the substation, respectively, where Vt D ,m,sub,sqrRepresents the square of the voltage amplitude of the substation connected to the transmission network node m;
the boundary connection constraint of the gas transmission-distribution system is as follows:
Figure FDA00033670401100000615
Figure FDA00033670401100000616
Figure FDA00033670401100000617
and
Figure FDA00033670401100000618
for the boundary connection constraint of the gas transmission-distribution system,
Figure FDA00033670401100000619
indicating the pressure at the boost side of the distribution pressure regulating station connected to the gas transmission network node n,
Figure FDA00033670401100000620
representing the pressure on the pressure-reducing side, gammacIs the voltage boosting ratio;
will Pt sub,D,m、
Figure FDA00033670401100000621
Vt D,m,sub,sqr、
Figure FDA00033670401100000622
Figure FDA00033670401100000623
And
Figure FDA00033670401100000624
and setting boundary connection variables of the power transmission-gas transmission system and the power distribution-gas distribution system in the recovery process.
CN202111385987.XA 2021-11-22 2021-11-22 Modeling method for layered collaborative recovery model of electricity-gas interconnection system Pending CN114065536A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111385987.XA CN114065536A (en) 2021-11-22 2021-11-22 Modeling method for layered collaborative recovery model of electricity-gas interconnection system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111385987.XA CN114065536A (en) 2021-11-22 2021-11-22 Modeling method for layered collaborative recovery model of electricity-gas interconnection system

Publications (1)

Publication Number Publication Date
CN114065536A true CN114065536A (en) 2022-02-18

Family

ID=80278868

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111385987.XA Pending CN114065536A (en) 2021-11-22 2021-11-22 Modeling method for layered collaborative recovery model of electricity-gas interconnection system

Country Status (1)

Country Link
CN (1) CN114065536A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115345391A (en) * 2022-10-20 2022-11-15 广东电网有限责任公司 Post-disaster recovery method and device for electric-gas energy system and storage medium

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115345391A (en) * 2022-10-20 2022-11-15 广东电网有限责任公司 Post-disaster recovery method and device for electric-gas energy system and storage medium

Similar Documents

Publication Publication Date Title
Van Hertem et al. Usefulness of DC power flow for active power flow analysis with flow controlling devices
Gil et al. Hierarchical frequency control scheme for islanded multi-microgrids operation
Lin et al. A restorative self-healing algorithm for transmission systems based on complex network theory
CN103730908B (en) Method for controlling energy storage converters in large-scale off-grid micro-grid
CN110034587B (en) Optimized scheduling method
CN105529701B (en) A kind of method for optimizing route of power up containing DC converter station based on artificial bee colony algorithm
CN112271726B (en) Power distribution system fault recovery method considering electricity-water-gas coupling relation
Chen et al. Impact of smart transformer voltage and frequency support in a high renewable penetration system
CN106992537A (en) A kind of power grid"black-start" method and device
CN109038532B (en) Distributed transmission and distribution cooperative unit combination method based on direct current power flow equivalence
CN114065536A (en) Modeling method for layered collaborative recovery model of electricity-gas interconnection system
El-Bahay et al. Computational methods to mitigate the effect of high penetration of renewable energy sources on power system frequency regulation: a comprehensive review
CN104124700A (en) Method and system for generating power distribution network black start scheme
Fahad et al. A coordination control strategy for power quality enhancement of an active distribution network
CN112018756A (en) Day-ahead robust coordinated optimization scheduling method for gas-electricity combined system
CN103997035A (en) Power transformer and line economic operation closed loop control system and method
Hemmati et al. Single-phase control of three-phase fuelcell-battery under unbalanced conditions considering off-grid and grid-tied states
CN110880791A (en) Coordination optimization method for hybrid alternating current-direct current power distribution network
CN105098774A (en) Method for generating power supply startup sequence of powered-off power grid containing direct-current converter station
CN117833273A (en) Power system look-ahead scheduling model solving method considering frequency constraint based on neural network assistance
CN110377973B (en) Construction method of standard linear comprehensive energy system model
CN112270096A (en) Flexible resource contribution degree solving method for promoting new energy consumption
Khoa et al. Experimental study on fault ride-through capability of VSC-based HVDC transmission system
Faraji et al. Centralized Multifunctional Control of Integrated Electricity–Gas Network in District Energy System Considering Resilience
CN103996149B (en) Regional power grid wind-powered electricity generation topological analysis method based on static electric voltage stability

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