Electricity-gas integrated energy system recovery decision-making method and device considering recovery process charged domain fusion

Abstract

The invention discloses an electricity-gas integrated energy system recovery decision-making method and device considering recovery process charged domain fusion in the technical field of electricity-gas integrated energy systems, and aims to solve the problem of recovery process charged domain fusion. The method comprises the following steps: constructing a sequence recovery problem constraint, an equipment constraint and a system operation constraint; based on the sequential recovery problem constraint, the equipment constraint and the system operation constraint, establishing a power-gas integrated energy system recovery decision model with the goal of maximizing the recovery net income of the power-gas integrated energy system; and linearizing the sequence recovery problem constraint, the equipment constraint and the system operation constraint, carrying out optimization solution on a target function of the power-gas integrated energy system recovery decision model based on the linearized constraints, and determining a power-gas integrated energy system recovery decision scheme. According to the method, the power-gas comprehensive energy system recovery decision model is established, and the mutual supporting capacity of the power system and the gas system in the recovery process can be fully exerted.

Description

Technical Field

[0001] This invention relates to a method and apparatus for making recovery decisions for an integrated electric-gas energy system that takes into account the fusion of charged domains during the recovery process, and belongs to the field of integrated energy system technology. Background Technology

[0002] In recent years, major power outages caused by natural disasters, cyberattacks, and human error have occurred globally. Extreme cold weather has led to a significant reduction in local natural gas supply, forcing numerous gas turbine units to shut down and causing widespread power outages. Simultaneously, the power outages have caused the shutdown of many electrically driven natural gas production and transportation facilities, further exacerbating the system's natural gas shortage. Therefore, researching power restoration decisions for Electricity-Gas Integrated Energy Systems (EGIES) will help reduce losses from electricity-gas outages.

[0003] During distribution network restoration, power sources with self-starting capability (SC) typically form energized domains. Multiple energized domains restore power independently in parallel to accelerate the restoration process. When available capacity is exhausted or there are no restorable loads, the restoration process terminates with each energized domain operating independently. This not only confines surplus available capacity and loads to be restored to each energized domain but is clearly not the desired final state of grid restoration. Considering energized domain fusion can not only expand the solution space for restoration decisions but also make the restoration scheme more complete. Therefore, EGIES power restoration decisions need to consider energized domain fusion and the bidirectional electro-pneumatic coupling characteristics to coordinate the allocation of electrical and pneumatic energy during the restoration process, thereby breaking down the barriers of multi-energy coupling in the power system and further expanding the solution space for restoration decisions. Summary of the Invention

[0004] The purpose of this invention is to propose a method and apparatus for decision-making in the recovery of an integrated electric-gas energy system that takes into account the fusion of charged domains during the recovery process, in order to address the problem of fusion of charged domains during the recovery process.

[0005] To achieve the above objectives, the present invention is implemented using the following technical solution:

[0006] In a first aspect, the present invention proposes a restoration decision-making method for an integrated electric-gas energy system that takes into account the fusion of charged domains during the restoration process, comprising:

[0007] Based on the energized domain fusion strategy and the number of steps during line operation, a sequential recovery problem constraint that satisfies the electro-pneumatic coupling constraint is constructed. The energized domain fusion strategy is used for step-by-step collaborative recovery of the power outage distribution network. The number of steps during line operation refers to the number of line switching operation steps that must be performed to energize critical loads.

[0008] Based on the physical characteristics of the integrated electric-gas energy system and its equipment, equipment constraints and system operation constraints are constructed.

[0009] Based on the constraints of the sequential recovery problem, equipment constraints, and system operation constraints, a recovery decision model for the integrated electric-gas energy system is established with the objective of maximizing the net recovery benefit of the integrated electric-gas energy system.

[0010] The constraints of the sequential recovery problem, equipment constraints, and system operation constraints are linearized. Based on the linearized constraints, the objective function of the electric-gas integrated energy system recovery decision model is optimized and solved to determine the electric-gas integrated energy system recovery decision scheme.

[0011] Furthermore, the sequential recovery problem constraint introduces an initial source node, which is defined as the first controllable power source node to be put into operation to form an energy recovery area during the recovery process of the integrated electric-gas energy system. This initial source node is used to select the best location based on competition and merge neighboring areas based on cooperation.

[0012] Furthermore, the electric-gas integrated energy system consists of a power grid system and a gas grid system. The power grid system includes a fixed power source, a gas turbine, and several power grid nodes. The gas grid system includes a fixed gas source, an electric-driven compressor, and several gas grid nodes. The power grid nodes supply power to the electric-driven compressor to support its operation, and the gas grid nodes supply gas to the gas turbine to support its power generation and supply the power grid node demand.

[0013] Furthermore, the constraints of the sequential recovery problem include:

[0014] (1)

[0015]

[0016]

[0017]

[0018]

[0019]

[0020]

[0021]

[0022]

[0023]

[0024]

[0025] (12)

[0026] (13)

[0027] (14)

[0028] (15)

[0029] (16)

[0030] (17)

[0031] (18)

[0032] (19)

[0033] in, A collection of branches in a power distribution network; It is the set of electrical nodes in the distribution network that have self-starting capability; For the power supply nodes of the distribution network; For the set of non-power source nodes in the distribution network; For the set of SC power supply nodes in the distribution network; This refers to the set of NSC power supply nodes in the distribution network that lack self-starting capability. It is a set of FSC power supply nodes in the distribution network that are fully self-starting. Let PSC be the set of power supply nodes in the distribution network that do not have full self-starting capability. The relationships between the sets are as follows:

[0034] (20)

[0035] The recovery state of branch ij at step t is represented by 1, indicating that it has been recovered and 0 indicates that it has not been recovered. The recovery state of step node i at time t, where 1 indicates that it has recovered and 0 indicates that it has not recovered; For time step t, whether SC power node i is a primary power node, 1 indicates that it is a primary power node, and 0 indicates that it is not a primary power node; and These are the sets of injection and outflow branches for electrical node i, respectively. and These are the virtual flows of the injection and outflow branches of electrical node i, respectively; This represents the virtual output of the SC power node; H is a positive number. This represents the recovery state of branch ij at step t+1; The recovery state of step node i at time t+1; T is the total number of observation steps; Let be the number of time steps elapsed from when branch ij is energized at the terminal node until restoration is complete. The recovery state of step node j at time t, where 1 indicates that it has recovered and 0 indicates that it has not recovered; The number of operation steps required to restore branch ij is the sum of the operation steps of all line switches on the branch; The recovery status of the step circuit branch ij at time T is 1, indicating that it has recovered and 0 indicates that it has not recovered. Is power node i of SC at time step t+1 the initial power node? The start / stop status of power supply m under step node i at time t, 1 indicates that it has been started, and 0 indicates that it has not been started; The set of fixed power supplies configured under node i; The number of operation steps required to start the fixed power supply m at node i; Let i be the set of FSC power sources at node i; Let i be the set of NSC power sources at node i.

[0036] Furthermore, the equipment constraints include gas turbine allocation constraints, mobile emergency power supply (MEPS) allocation constraints, fixed power supply output constraints, fixed gas source constraints, electric compressor start / stop status constraints, mobile emergency power supply (MEPS) output constraints, mobile emergency gas source allocation constraints, electrochemical energy storage constraints, and gas storage tank constraints. The equipment includes a gas turbine, mobile emergency power supply, fixed power supply, fixed gas source, electric compressor, mobile emergency gas source, electrochemical energy storage, and gas storage tank.

[0037] Furthermore, the system operation constraints satisfy the load recovery state constraints:

[0038] (twenty one)

[0039] (twenty two)

[0040] (twenty three)

[0041] in, and These represent the recovery states of node i's load at time step t and time step t+1, respectively. and These represent the recovery states of node g's load at time step t and time step t+1, respectively. Let i be the recovery state of the step node i at time t. and These are the sets of load nodes for the power distribution network and the gas distribution network, respectively.

[0042] Furthermore, the system operation constraints also include power flow constraints of the distribution network and dynamic power flow constraints of the gas distribution network;

[0043] The power flow constraints of the distribution network include:

[0044] (twenty four)

[0045] (25)

[0046] (26)

[0047] (27)

[0048] (28)

[0049] (29)

[0050] (30)

[0051] (31)

[0052] (32)

[0053] (33)

[0054] (34)

[0055] in, For the gas node g, which is connected to the electro-pneumatic coupling device at electrical node i, A collection of branches in a power distribution network; It is the set of electrical nodes in the distribution network that have self-starting capability; and These are the sets of injection and outflow branches for electrical node i, respectively. and These are the active power and reactive power of branch ij at time step t, respectively. and These represent the active power and reactive power flowing out of step node i to the outside at time t, respectively. and These represent the active power and reactive power flowing into step node i at time t, respectively. and Let be the resistance and reactance of branch ij, respectively; Let be the voltage at step node i at time t; This serves as the reference voltage in the power system. As an auxiliary variable; Let H represent the recovery state of branch ij at step t; H is a positive number. Let be the maximum apparent power that can pass through branch ij; and These are the maximum and minimum allowable voltages for electrical node i, respectively; The power demand of step-electric node i at time t after the load participates in IDR, excluding power supply and gas supply equipment; This refers to the set of NSC power supply nodes in the distribution network that lack self-starting capability. Let m be the start / stop state of power supply m at step node i at time t. The active power required for the N-SC power supply m at node i to start; Let i be the electrochemical energy storage collection at electrical node i. Let m be the charging state of the electrochemical energy storage at step node i at time t. The active power required to charge the electrochemical energy storage m at the step node i at time t; The fixed electrically driven gas source set at gas node g. This represents the start / stop state of the fixed gas source n at the gas stepping node g at time t. The power consumed by the fixed electric drive air source n at step node i at time t; This refers to the collection of P2G devices at gas node g. Let t be the power consumed by P2G device n at step node i; The start / stop state of the electrically driven compressor at step gas node g at time t is given. Let t be the power consumed by the compressor driven by step node i when step node i is powered on; Let i be the reactive power demand factor of electrical node i; and These are the active power and reactive power of branch ki at time step t, respectively. Let be the voltage at step node j at time t; Let be the maximum active power that can pass through branch ij; Let be the maximum reactive power that can pass through branch ij;

[0056] The dynamic flow constraints of the gas distribution network include dynamic flow constraints of natural gas, node flow balance and node pressure constraints, flow and pressure constraints in pipelines, and node gas pressure constraints.

[0057] The dynamic power flow constraint for natural gas is:

[0058] (35)

[0059] (36)

[0060] in, A collection of pipes for the gas distribution network; Let gh be the cross-sectional area of ​​the pipe; For time step; This refers to the length of the pipe. The average flow velocity of natural gas in pipeline gh; Where gh is the inner diameter of the pipe; and These represent the flow rates at the end and beginning of the pipeline at time step t, respectively. and These are the end pressure and beginning pressure of the pipeline at time t, respectively. This refers to the pipe friction coefficient; The speed of sound in a gas;

[0061] The node flow balancing and node pressure constraints are as follows:

[0062] (37)

[0063] (38)

[0064] (39)

[0065] (40)

[0066] in, Let t be the flow rate of gas injected into node g from the outside; Let g be the gas flow rate of the object to be restored at the gas node g at time t; and These are the sets of injection and outflow branches for gas node g, respectively. Let fg be the flow rate at the end of pipe at time step t; The recovery state of the load at node g at time step t; To calculate the flow demand of the gas node g at time t after adjustment via IDR, excluding power supply and gas supply equipment. Let be the set of GT at electrical node i; Let be the gas flow rate of the gas turbine m connected to the gas stepping node g at time t; Let g be the set of GT at the gas node g; The gas storage state of gas storage tank n at gas step node g at time t; Let n be the gas storage flow rate at gas storage tank n at gas step node g at time t; Let g be the air pressure at step node g at time t; A collection of air nodes equipped with electrically driven compressors; This refers to the set of gas nodes in the gas distribution network. Let fg be the flow rate at the end of pipe at time step t; , where are the flow rates at the beginning of pipe gh at time step t;

[0067] The flow and pressure constraints in the pipeline are as follows:

[0068] (41)

[0069] (42)

[0070] (43)

[0071] (44)

[0072] in, The maximum permissible mass flow rate for pipeline gh; and These are the upper and lower limits of the allowable pressure for pipeline gh, respectively.

[0073] The node air pressure constraint is:

[0074] (45)

[0075] (46)

[0076] (47)

[0077] (48)

[0078] in, Electrical node i is connected to the electro-pneumatic coupling device at gas node g. and These are the minimum and maximum air pressures allowed at node g, respectively; Let g be the air pressure at the stepping node; The power consumed by the compressor driven by the pipeline gh at time step t; K1, K2, and K3 are the outlet pressure setpoint of the electrically driven compressor at node g; K1, K2, and K3 are the compressor equation coefficients and exponents. Let t be the power consumed by the compressor driven by step node i at time t.

[0079] Furthermore, the objective function of the electric-gas integrated energy system recovery decision model is:

[0080] (49)

[0081] in, The net benefit of obtaining electricity or gas supply in advance for an integrated electric-gas energy system is calculated using the following formula:

[0082] (50)

[0083] in, and These are the sets of load nodes for the power distribution network and the gas distribution network, respectively. Let be the unit recovery revenue of load node i at time step t; This represents the recovery status of the load at node i at time step t, where 1 indicates recovery and 0 indicates no recovery. The electrical energy to be restored at node i at time step t, excluding power supply and gas supply equipment; Let g be the unit recovery revenue of the gas load node g at time t; This represents the recovery status of the load at node g at time step t, where 1 indicates recovery and 0 indicates no recovery. Let g be the volume of natural gas to be restored at time step t, excluding power and gas supply equipment.

[0084] Furthermore, the linearized sequential recovery problem constraints, device constraints, and system operation constraints include:

[0085] For the nonlinear terms appearing in the constraint system, equipment constraints, and system operation constraints of the sequence recovery problem, auxiliary variable constraints are introduced for linearization, transforming the original mixed integer nonlinear programming problem into a mixed integer linear programming problem.

[0086] Secondly, this invention proposes a power-gas integrated energy system recovery decision-making device that considers the fusion of charged domains during the recovery process, comprising:

[0087] The constraint module is used to construct sequential recovery problem constraints, device constraints, and system operation constraints;

[0088] The objective function construction module for the recovery decision model is used to establish the decision model and objective function of the integrated electric-gas energy system with the goal of maximizing the net recovery benefit of the integrated electric-gas energy system.

[0089] The recovery decision scheme determination module is used to linearize the nonlinear constraints in the constraint module and optimize the objective function of the recovery decision model to determine the recovery decision scheme of the integrated electric-gas energy system.

[0090] Compared with the prior art, the beneficial effects achieved by the present invention are as follows:

[0091] (1) For the problem of fusion of charged domains, this invention proposes the concept of the initial node and its competitive and cooperative capabilities, and improves the constraint on the branch recovery delay in the traditional sequential recovery model. On this basis, considering the impact of the natural gas system and the transportation delay of mobile emergency power supply, a distribution network sequential recovery model that takes into account the fusion of charged domains, the order of equipment operation and the operation time is constructed.

[0092] (2) Considering the bidirectional electro-pneumatic coupling characteristics of EGIES, this invention takes into account various types of electro-pneumatic and pneumatic-electrical equipment on the source and load sides, and takes into account the dynamic characteristics of the pneumatic system and the fusion of charged domains. An EGIES recovery decision model is established and converted into a mixed integer linear programming model, which can be solved using existing commercial solvers. The proposed model can fully leverage the mutual support capabilities of the electro-pneumatic system during the recovery process. Attached Figure Description

[0093] Figure 1 A flowchart illustrating the electric-gas integrated energy system recovery decision-making method based on the integration of the electric domain during the recovery process, as provided in an embodiment of the present invention;

[0094] Figure 2 This is a schematic diagram of the charged domain fusion strategy provided in an embodiment of the present invention;

[0095] Figure 3 This is a schematic diagram of an integrated electric-gas energy system provided in an embodiment of the present invention. Detailed Implementation

[0096] The technical solution of the present invention will be described in detail below with reference to the accompanying drawings and specific embodiments. The same reference numerals in the drawings indicate the same or similar components or parts. Those skilled in the art should understand that these drawings are not necessarily drawn to scale. The embodiments and specific features within the embodiments of this application are detailed descriptions of the technical solution of this application, and not limitations thereof. Where there is no conflict, the embodiments and technical features within the embodiments of this application can be combined with each other.

[0097] In this article, the term "and / or" is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A existing alone, A and B existing simultaneously, or B existing alone. Additionally, the character " / " in this article generally indicates that the preceding and following related objects have an "or" relationship.

[0098] Example 1:

[0099] This embodiment proposes a restoration decision-making method for an integrated electric-gas energy system that considers the fusion of charged domains during the restoration process. Figure 1 As shown, it includes:

[0100] Based on the energized domain fusion strategy and the number of steps during line operation, a sequential recovery problem constraint that satisfies the electro-pneumatic coupling constraint is constructed. The energized domain fusion strategy is used for the step-by-step collaborative recovery of the power outage distribution network. The number of steps during line operation refers to the number of "steps" or "batches" of line switches that need to be operated (such as closing) in order to supply power to critical loads during the process of gradually and orderly restoring the power grid from a large-scale power outage to a normal operating state.

[0101] Based on the physical characteristics of the integrated electric-gas energy system and its equipment, equipment constraints and system operation constraints are constructed.

[0102] Based on the constraints of sequential recovery problem, equipment constraints, and system operation constraints, a recovery decision model for an integrated electric-gas energy system is established with the objective of maximizing the net recovery benefit of the integrated electric-gas energy system.

[0103] The constraints of sequential recovery, equipment, and system operation are linearized. Based on the linearized constraints, the objective function of the integrated power-gas energy system recovery decision model is optimized to determine the integrated power-gas energy system recovery decision scheme. This scheme covers the recovery status and action arrangements of the grid side and the gas grid side at each time step within the observation period. On the grid side, this includes the power supply status of each node and load at each time step, the commissioning status of each line, the recovery progress of various power sources except for fully self-starting FSC power sources, the nodes and time steps of mobile emergency power source access and their subsequent output plans, as well as the start-up and shutdown status and power supply of electrochemical energy storage. On the gas grid side, this includes the gas supply status of each node and load at each time step, the pipeline on / off status, the gas supply recovery status of various gas sources, the nodes and time steps of mobile emergency gas source access and their subsequent gas supply arrangements, as well as the start-up and shutdown status and gas supply of gas storage tanks.

[0104] In this embodiment, the integrated electric-gas energy system consists of a power grid system and a gas grid system, and its structure is as follows: Figure 3 As shown, the power grid system includes the main grid power source, gas turbine, and several grid nodes, while the gas grid system includes the gas source, several gas grid nodes, and an electrically driven compressor. The two systems are coupled together. The grid nodes supply power to the compressor to support its operation, and the gas grid nodes supply gas to the gas turbine to support the gas turbine's power generation and supply the grid node's demand.

[0105] In this embodiment, a self-starting capability (SC) source that does not rely on other energy systems for power supply during the recovery process is defined as a fully self-starting capability (FSC) source, such as a thermal power unit with self-starting capability and no fuel limitations, a non-electrically driven gas source, and a mobile emergency source (MES) pre-installed at a designated location. An SC source whose power supply during the recovery process is constrained by other energy systems is defined as a partially self-starting capability (PSC) source, such as a gas turbine (GT), an electrically driven gas source, a power-to-gas (P2G) device, and a MES that relies on traffic network status adjustments during the recovery process. The nodes are further classified. Nodes with fixed power / gas supply and candidate MES connection points are designated as source nodes. Based on the self-starting capability of the power / gas supply, they can be further divided into FSC source nodes, PSC source nodes, and source nodes without self-starting capability (NSC). Other nodes are designated as non-source nodes.

[0106] In this embodiment, the constraint system for the sequential recovery problem introduces the initial source node. Let the SC source node initially deployed to form the energy recovery region be the initial source node of that region, and the other source nodes be non-initial source nodes. The ability of the initial source node to be deployed first, before non-initial source nodes, is called the competitive capability of the initial source node. Assuming that there is only one initial source node in an energy recovery region, during the merging of multiple energy recovery regions, only one initial source node is retained, and the remaining original initial source nodes will be converted into non-initial source nodes. This conversion process is called the cooperative capability of the initial source nodes. All SC source nodes possess the competitive capability to become initial source nodes, and all initial source nodes possess the cooperative capability. For EGIES with many and densely distributed source nodes, if the competitive and cooperative capabilities of SC source nodes to become initial source nodes are ignored, then all SC source nodes will generate independent energy recovery regions, which is neither conducive to the aggregation of available energy nor to its rapid and full utilization.

[0107] In this embodiment, the constraints of the sequential recovery problem include:

[0108] (1)

[0109]

[0110]

[0111]

[0112]

[0113]

[0114]

[0115]

[0116]

[0117]

[0118]

[0119] (12)

[0120] (13)

[0121] (14)

[0122] (15)

[0123] (16)

[0124] (17)

[0125] (18)

[0126] (19)

[0127] in, A collection of branches in a power distribution network; It is the set of electrical nodes in the distribution network that have self-starting capability; For the power supply nodes of the distribution network; For the set of non-power source nodes in the distribution network; For the set of SC power supply nodes in the distribution network; This refers to the set of NSC power supply nodes in the distribution network that lack self-starting capability. It is a set of FSC power supply nodes in the distribution network that are fully self-starting. Let PSC be the set of power supply nodes in the distribution network that do not have full self-starting capability. The relationships between the sets are as follows:

[0128] (20)

[0129] The recovery state of branch ij at step t is represented by 1, indicating that it has been recovered and 0 indicates that it has not been recovered. The recovery state of step node i at time t, where 1 indicates that it has recovered and 0 indicates that it has not recovered; For time step t, whether SC power node i is a primary power node, 1 indicates that it is a primary power node, and 0 indicates that it is not a primary power node; and These are the sets of injection and outflow branches for electrical node i, respectively. and These are the virtual flows of the injection and outflow branches of electrical node i, respectively; This represents the virtual output of the SC power node; H is a positive number. This represents the recovery state of branch ij at step t+1; The recovery state of step node i at time t+1; T is the total number of observation steps; Let be the number of time steps elapsed from when branch ij is energized at the terminal node until restoration is complete. The recovery state of step node j at time t, where 1 indicates that it has recovered and 0 indicates that it has not recovered; The number of operation steps required to restore branch ij is the sum of the operation steps of all line switches on the branch; The recovery status of the step circuit branch ij at time T is 1, indicating that it has recovered and 0 indicates that it has not recovered. Is power node i of SC at time step t+1 the initial power node? The start / stop status of power supply m under step node i at time t, 1 indicates that it has been started, and 0 indicates that it has not been started; The set of fixed power supplies configured under node i; The number of operation steps required to start the fixed power supply m at node i; Let i be the set of FSC power sources at node i; Let i be the set of NSC power sources at node i.

[0130] During the system power restoration process, the topology of the restored network should remain connected. Equations (2)-(4) constrain the connectivity of the system power restoration process; Equation (5) constrains the virtual output of the SC power nodes; Equations (6)-(7) indicate that nodes and branches in the system should remain in the restored state after restoration; Equations (8)-(9) indicate that the nodes at both ends of a branch should also be restored after restoration; Equations (10)-(11) indicate that a sufficient condition for the restoration of a branch in the system is that at least one of its nodes has been restored; Equation (12) indicates that a branch needs to go through a certain number of operation steps to become a restored state; Equation (13) indicates that the number of operation steps required for the restoration of a branch must be greater than or equal to the number of operation steps of the branch; Equation (14) indicates that, except for mobile emergency power (Mobile Emergency Power) The necessary condition for the recovery of a power node outside the candidate connection point of the Source (MEPS) is that its connected node is recovered; Equation (15) indicates that the necessary but not sufficient condition for an SC power node to be a primary source node is that a power source is recovered; Equation (16) indicates that an SC power node is allowed to switch from a non-primary source node to a primary source node, but once a primary source node switches to a non-primary source node, it cannot switch back to a primary source node. Equation (17) constrains the continuous operation of a fixed power supply after startup; Equation (18) constrains the self-starting time of the FSC power supply; Equation (20) indicates that an NSC power supply can only recover after its adjacent node has recovered and after its startup time has elapsed.

[0131] Charged domain fusion strategy such as Figure 2 As shown, this strategy uses the SC power node as the initial source and takes into account the number of steps during line operation to gradually restore the power supply of the distribution network. During this process, all restored paths are connected, and the restoration of all nodes except the initial source node depends on at least one restored adjacent node. That is, the power supply sequence restoration path in the distribution network power restoration process should satisfy both connectivity and operation order conditions. Step S0 represents the initial power outage scenario. In step S1, nodes 1 and 10 become the initial source nodes of two electrified domains due to the first startup of DG1 and DG3, respectively. Node 9, due to the longer startup time of DG2, failed to compete with the initial source node 10 and became a non-initial source node in step S2. In step S3, the gray and orange electrified domains in step S2 are merged, and the initial source node 10 becomes a non-initial source node.

[0132] In this embodiment, the equipment constraints and system operation constraints are mathematically modeled based on the physical characteristics of the equipment and system. Taking the gas turbine as an example in the equipment operation mode, its physical operation characteristics can be referenced for modeling. For example, the output of the gas turbine has upper and lower limits, which can be modeled as an inequality: lower limit <= output <= upper limit. The modeling of other equipment constraints and system operation constraints is also based on physical characteristics.

[0133] Equipment constraints include gas turbine allocation constraints, mobile emergency power supply (MEPS) allocation constraints, stationary power supply output constraints, stationary gas source constraints, electric compressor start / stop status constraints, mobile emergency power supply (MEPS) output constraints, mobile emergency gas source allocation constraints, electrochemical energy storage constraints, and gas storage tank constraints. Among these, the equipment includes gas turbines, mobile emergency power supplies, stationary power supplies, stationary gas sources, electric compressors, mobile emergency gas sources, electrochemical energy storage, and gas storage tanks.

[0134] Assuming that all gas turbines (GTs) in this invention have self-starting capability; the competitiveness of an electric node equipped with a gas turbine (GT) to become a primary energy source node depends on whether the gas system can supply it with energy and whether the gas pressure at its gas acquisition node reaches a specified value, with the following allocation constraints:

[0135] (twenty one)

[0136] (twenty two)

[0137] in, For a set of electric nodes configured with GT; Let be the set of GT at electrical node i; For electrical node i, gas node g is connected via GT / electrically driven compressor / electrically driven gas source / P2G device; The recovery status of the gas node g at time t, where 1 indicates recovery and 0 indicates non-recovery; Let g be the air pressure at step node g at time t; and These are the minimum and maximum allowable gas pressures for the gas turbine m at gas node g, respectively;

[0138] The competitiveness of a candidate connection point equipped with a Mobile Emergency Power Supply (MEPS) to become a primary source node depends on whether it has been successfully connected to the power grid system. Its allocation constraints are as follows:

[0139] (twenty three)

[0140] (twenty four)

[0141] (25)

[0142] (26)

[0143] (27)

[0144] (28)

[0145] (29)

[0146] (30)

[0147] (31)

[0148] (32)

[0149] in, For the MEPS set; Let m be the set of candidate connection points for MEPS. To observe the startup status of MEPS m at its candidate connection point i at the end of the time step, 1 indicates that it has started and 0 indicates that it has not started; and These are the start time and initiation time of MEPS m to its candidate connection point i, respectively; The number of time steps required for MEPS m to be transported to the candidate connection point i and started; The start / stop status of MEPS m at candidate connection point i at time step t is 1, indicating that it has started and 0 indicates that it has not started.

[0150] During the process of transporting MEPS from its storage location to the candidate connection point and starting the grid connection after a preparation period, the departure time and start time involved can be calculated by equations (23)-(26), where the start time is the departure time plus the number of time steps required for the MEPS to be transported to the candidate connection point and start. For candidate connection points that have not been connected to MEPS, both the departure time and start time must be set to T. Equation (27) indicates that MEPS can only be connected to one candidate connection point at most; Equation (328) indicates that MEPS candidate connection points can only be connected to one MEPS at most; Equation (29) indicates that the necessary but not sufficient condition for a MEPS candidate connection point to be a source node is that MEPS is started; Equation (30) is the consistency constraint between the time step state of each MEPS candidate connection point and the start time of the MEPS it is connected to; Equation (31) indicates that MEPS must continue to run after it is started; Equation (32) indicates that the necessary but not sufficient condition for the recovery of a MEPS candidate connection point is that MEPS is started.

[0151] Fixed power supply output constraints include:

[0152] (33)

[0153] (34)

[0154] (35)

[0155] in, and These are the active power and reactive power output by the fixed power supply m connected to electrical node i, respectively. , , , These are the upper and lower limits of the active and reactive power output of the fixed power supply m connected to electrical node i, respectively. , These are the ramping constraints for the fixed power supply m connected to electrical node i.

[0156] Equations (33)-(34) provide the output constraints of the fixed power source; Equation (35) is the ramp constraint of the fixed power source.

[0157] The energy conversion relationship constraint of GT is shown in Equation (36).

[0158] (36)

[0159] in, For the energy conversion parameters of the gas turbine m connected to electrical node i; Let be the gas flow rate of the gas turbine m connected to the gas stepping node g at time t; Electrical node i is connected to the GT / electrically driven compressor / electrically driven air source / P2G device at air node g.

[0160] Fixed gas source constraints include:

[0161] (37)

[0162] (39)

[0163] (39)

[0164] (40)

[0165] (41)

[0166] in, This refers to the set of SC gas source nodes in the gas distribution network. This refers to the set of NSC gas source nodes in the gas distribution network. Let g be the set of gas sources at gas node g; For the collection of P2G devices at gas node g; A collection of air nodes equipped with a fixed electrically driven air source; A collection of gas nodes equipped with P2G devices; A fixed set of electrically driven air sources at air node g; The start / stop status of the fixed gas source n at the gas stepping node g at time t, 1 indicates that it has been started, and 0 indicates that it has not been started; Let n be the air flow rate output by the fixed air source n at the step air node g at time t; and These are the upper and lower limits of the output flow rate of the fixed gas source n at gas node g, respectively; Set the outlet pressure of the fixed gas source at gas node g; The number of operation steps required to start the fixed gas source n at gas node g.

[0167] Equation (37) indicates that the stationary gas source other than the P2G device can continue to operate after startup; Equations (38)-(39) constrain the flow rate and pressure output of the stationary gas source; Equation (40) indicates that the stationary electric drive gas source and the P2G device can only start after the power supply to their connected electrical nodes is restored; Equation (41) indicates that the electric drive gas source can only start after the number of time steps required for its self-starting.

[0168] The energy conversion relationship of P2G devices is shown in equation (42).

[0169] (42)

[0170] in, For the energy conversion parameters of P2G device n at gas node g; Let t be the power consumed by P2G device n at step node i.

[0171] The start-stop state constraints of the electrically driven compressor include:

[0172] (43)

[0173] (44)

[0174] in, A collection of air nodes equipped with electrically driven compressors; The start / stop status of the electrically driven compressor at the step gas node g at time t is 1, which indicates that it has started and 0 indicates that it has not started.

[0175] Equation (43) indicates that the electric compressor can only start after the power supply to its connected electrical node is restored; Equation (44) indicates that the electric compressor must remain in the started state after it has started.

[0176] MEPS output constraints include:

[0177] (45)

[0178] (46)

[0179] in, , , and These are the upper and lower limits of active and reactive power for MEPS m, respectively. and Let be the active and reactive power output of MEPS m at candidate connection point i at time step t, respectively.

[0180] Equations (45) and (46) are the output constraints of MEPS.

[0181] The preparation process for mobile emergency gas sources (MEGS) is similar to that for MEPS, including:

[0182] (47)

[0183] (48)

[0184] (49)

[0185] (50)

[0186] (51)

[0187] (52)

[0188] (53)

[0189] (54)

[0190] (55)

[0191] in, For the MEGS set; For the set of candidate join points in MEGS; and These are the start and initiation steps of MEGSn for its candidate connection point g, respectively. The number of time steps required for MEGS n to be transported to the candidate connection point g and started; To observe the access status of MEGS n at its candidate connection point g at the end of the time step, 1 indicates that it has been started and 0 indicates that it has not been started; The access status of MEGS n at time step t at candidate connection point g is 1, indicating that it has been started and 0 indicates that it has not been started. , These are the upper and lower limits of the gas flow rate for MEGS n, respectively; Let t be the gas flow rate of MEGS n at the candidate connection point g.

[0192] Equations (47)-(50) are constraints on the start time and activation time of MEGS; Equation (51) indicates that MEGS can only connect to one candidate connection point at most; Equation (52) indicates that a candidate connection point of MEGS can only connect to one MEGS at most; Equation (53) is a constraint on the consistency of the state of each time step of the candidate connection point of MEGS with the activation time of the MEGS it connects to; Equation (54) indicates that MEGS must maintain the connected state after it is connected; Equation (55) is a constraint on the output of MEGS.

[0193] Electrochemical energy storage constraints include:

[0194] (56)

[0195] (57)

[0196] (58)

[0197] (59)

[0198] (60)

[0199] (61)

[0200] (62)

[0201] (63)

[0202] in, A collection of electrical nodes equipped with energy storage; Let i be the electrochemical energy storage set at electrical node i; and These represent the active power and reactive power required to charge the electrochemical energy storage m at step node i at time t, respectively. and These represent the active power and reactive power released by the electrochemical energy storage m at step node i at time t, respectively. and These represent the charging and discharging states of the electrochemical energy storage m at step node i at time t, where 1 indicates a charging / discharging state and 0 indicates that it is not in a charging / discharging state. Let m be the initial state of charge of the electrochemical energy storage m at electrical node i. Let m be the state of charge of the electrochemical energy storage at step node i at time t. and , respectively, are the charge and discharge efficiencies of the electrochemical energy storage m at electrical node i; Let m be the capacity of the electrochemical energy storage at electrical node i; and These represent the upper and lower limits of the allowable state of charge (SOC) for electrochemical energy storage m at electrical node i, respectively.

[0203] Equations (56)-(59) are the output constraints of energy storage; Equation (60) is the charging and discharging state constraints of energy storage; Equations (61) and (62) are used to calculate the state of charge of energy storage; Equation (63) is the state of charge constraints of energy storage.

[0204] Gas storage tank constraints include:

[0205] (64)

[0206] (65)

[0207] (66)

[0208] (67)

[0209] (68)

[0210] (69)

[0211] in, A set of gas nodes equipped with gas storage tanks; The collection of gas storage tanks at gas node g; and These represent the gas storage flow rate and the gas release flow rate at gas node g at time t, respectively; and These represent the gas storage and venting states of gas storage tank n at gas step node g at time t, where 1 indicates gas storage / venting state and 0 indicates that it is not in gas storage / venting state. Let g be the initial volume of natural gas in gas storage tank n at gas node g; Let n be the volume of natural gas in storage tank n at step gas node g at time t; The density of natural gas; and These represent the upper and lower limits of the permissible natural gas volume for gas storage tank n at gas node g, respectively.

[0212] Equations (64)-(65) are constraints on the output of the gas storage tank; Equation (66) is a constraint on the storage and release state of energy storage; Equations (67)-(68) are used to calculate the natural gas volume of the gas storage tank; Equation (69) is a constraint on the natural gas volume of the gas storage tank.

[0213] The system operation constraints satisfy the load recovery state constraints:

[0214] (70)

[0215] (71)

[0216] (72)

[0217] in, and These represent the recovery states of node i's load at time step t and time step t+1, respectively. and These represent the recovery states of node g's load at time step t and time step t+1, respectively. Let i be the recovery state of the step node i at time t. and These are the sets of load nodes for the power distribution network and the gas distribution network, respectively.

[0218] In this embodiment, the system operation constraints also include power flow constraints of the distribution network and dynamic power flow constraints of the gas distribution network;

[0219] This embodiment uses linearized Distflow power flow constraints, which include:

[0220] (73)

[0221] (74)

[0222] (75)

[0223] (76)

[0224] (77)

[0225] (78)

[0226] (79)

[0227] (80)

[0228] (81)

[0229] (82)

[0230] (83)

[0231] in, For the gas node g, which is connected to the electro-pneumatic coupling device at electrical node i, A collection of branches in a power distribution network; It is the set of electrical nodes in the distribution network that have self-starting capability; and These are the sets of injection and outflow branches for electrical node i, respectively. and These are the active power and reactive power of branch ij at time step t, respectively. and These represent the active power and reactive power flowing out of step node i to the outside at time t, respectively. and These represent the active power and reactive power flowing into step node i at time t, respectively. and Let be the resistance and reactance of branch ij, respectively; Let be the voltage at step node i at time t; This serves as the reference voltage in the power system. As an auxiliary variable; Let H represent the recovery state of branch ij at step t; H is a positive number. Let be the maximum apparent power that can pass through branch ij; and These are the maximum and minimum allowable voltages for electrical node i, respectively; The power demand of step-electric node i at time t after the load participates in IDR, excluding power supply and gas supply equipment; This refers to the set of NSC power supply nodes in the distribution network that lack self-starting capability. Let m be the start / stop state of power supply m at step node i at time t. The active power required for the N-SC power supply m at node i to start; Let i be the electrochemical energy storage collection at electrical node i. Let m be the charging state of the electrochemical energy storage at step node i at time t. The active power required to charge the electrochemical energy storage m at the step node i at time t; The fixed electrically driven gas source set at gas node g. This represents the start / stop state of the fixed gas source n at the gas stepping node g at time t. The power consumed by the fixed electric drive air source n at step node i at time t; This refers to the collection of P2G devices at gas node g. Let t be the power consumed by P2G device n at step node i; The start / stop state of the electrically driven compressor at step gas node g at time t is given. Let t be the power consumed by the compressor driven by step node i when step node i is powered on; Let i be the reactive power demand factor of electrical node i; and These are the active power and reactive power of branch ki at time step t, respectively. Let be the voltage at step node j at time t; Let be the maximum active power that can pass through branch ij; Let be the maximum reactive power that can pass through branch ij;

[0232] Equations (73)-(74) are the active and reactive power balance constraints of the electrical node; Equations (75)-(77) are the voltage relationship constraints of the restored branch ij; Equation (78) is the apparent power constraint of the branch; Equations (79)-(81) are the node voltage and branch power constraints; Equations (82)-(83) are the active and reactive power constraints flowing out of the node.

[0233] Dynamic flow constraints in gas distribution networks include dynamic flow constraints for natural gas, node flow balance and node pressure constraints, flow and pressure constraints in pipelines, and node gas pressure constraints.

[0234] To illustrate the mathematical expression of these constraints, it is necessary to first clarify the topology of the gas network and the definition of the positive direction of the parameters. The assumed positive direction rule is similar to that of the power system. Therefore, pipelines fg and gh are composed of nodes f, g, and h, where node f is the first node of pipeline fg, pipeline fg is the outflow pipeline from node f, node g is the last node of pipeline fg, and pipeline fg is the inflow pipeline from node g. In this embodiment, the assumed positive direction of the mass flow rate in pipeline gh is assumed to be the same as the assumed positive direction of the gas system topology. Regarding the pressure and flow rate in the pipeline, this embodiment defines the mass flow rate flowing out of the node g partition along the assumed positive direction at node g as the first-end flow rate of pipeline gh, and the corresponding pressure as the first-end pressure of pipeline gh. The mass flow rate flowing into the node h partition along the assumed positive direction at the node is called the last-end flow rate of the pipeline, and the corresponding pressure as the last-end pressure of the pipeline.

[0235] The dynamic flow constraint for natural gas is:

[0236] (84)

[0237] (85)

[0238] in, A collection of pipes for the gas distribution network; Let gh be the cross-sectional area of ​​the pipe; For time step; This refers to the length of the pipe. The average flow velocity of natural gas in pipeline gh; Where gh is the inner diameter of the pipe; and These represent the flow rates at the end and beginning of the pipeline at time step t, respectively. and These are the end pressure and beginning pressure of the pipeline at time t, respectively. This refers to the pipe friction coefficient; The speed of sound in a gas;

[0239] The node traffic balancing and node pressure constraints are as follows:

[0240] (86)

[0241] (87)

[0242] (88)

[0243] (89)

[0244] in, Let t be the flow rate of gas injected into node g from the outside; Let g be the gas flow rate of the object to be restored at the gas node g at time t; and These are the sets of injection and outflow branches for gas node g, respectively. Let fg be the flow rate at the end of pipe at time step t; The recovery state of the load at node g at time step t; To calculate the flow demand of the gas node g at time t after adjustment via IDR, excluding power supply and gas supply equipment. Let be the set of GT at electrical node i; Let be the gas flow rate of the gas turbine m connected to the gas stepping node g at time t; Let g be the set of GT at the gas node g; The gas storage state of gas storage tank n at gas step node g at time t; Let n be the gas storage flow rate at gas storage tank n at gas step node g at time t; Let g be the air pressure at step node g at time t; A collection of air nodes equipped with electrically driven compressors; This refers to the set of gas nodes in the gas distribution network. Let fg be the flow rate at the end of pipe at time step t; , where are the flow rates at the beginning of pipe gh at time step t;

[0245] The flow and pressure constraints in the pipeline are as follows:

[0246] (90)

[0247] (91)

[0248] (92)

[0249] (93)

[0250] in, The maximum permissible mass flow rate for pipeline gh; and These are the upper and lower limits of the allowable pressure for pipeline gh, respectively.

[0251] The nodal pressure constraint is:

[0252] (94)

[0253] (95)

[0254] (96)

[0255] (97)

[0256] in, Electrical node i is connected to the electro-pneumatic coupling device at gas node g. and These are the minimum and maximum air pressures allowed at node g, respectively; Let g be the air pressure at the stepping node; The power consumed by the compressor driven by the pipeline gh at time step t; K1, K2, and K3 are the outlet pressure setpoint of the electrically driven compressor at node g; K1, K2, and K3 are the compressor equation coefficients and exponents. Let t be the power consumed by the compressor driven by step node i at time t.

[0257] Equation (94) indicates that gas load can only be put into the node when the node gas pressure reaches the specified value; Equation (95) indicates that the outlet pressure of the electric drive compressor is the set value when it starts; Equation (96) indicates that the outlet pressure and inlet pressure of the electric drive compressor are equal when it is not started; Equation (97) indicates the relationship between the flow rate, inlet pressure and outlet pressure of the electric drive compressor and its power consumption.

[0258] In this embodiment, the objective function of the integrated electric-gas energy system recovery decision model is:

[0259] (98)

[0260] in, The net benefit of obtaining electricity or gas supply in advance for an integrated electric-gas energy system is calculated using the following formula:

[0261] (99)

[0262] in, and These are the sets of load nodes for the power distribution network and the gas distribution network, respectively. Let be the unit recovery revenue of load node i at time step t; This represents the recovery status of the load at node i at time step t, where 1 indicates recovery and 0 indicates no recovery. The electrical energy to be restored at node i at time step t, excluding power supply and gas supply equipment; Let g be the unit recovery revenue of the gas load node g at time t; This represents the recovery status of the load at node g at time step t, where 1 indicates recovery and 0 indicates no recovery. Let g be the volume of natural gas to be restored at time step t, excluding power and gas supply equipment.

[0263] In this embodiment, the objective function and constraints of the linearized electric-gas integrated energy system recovery decision model include:

[0264] For the nonlinear terms appearing in the sequential recovery problem constraints, equipment constraints, and system operation constraints, auxiliary variable constraints are introduced for linearization, transforming the original mixed-integer nonlinear programming problem into a mixed-integer linear programming problem.

[0265] For example, for the nonlinear term max(a,b) in equation (13), an auxiliary variable is introduced. , , Substitute max(a,b) as shown in equation (100):

[0266] (100)

[0267] For the multiplication of a 0-1 variable with a continuous variable, the formula (82) is used. For example, introduce auxiliary variables. replace After linearization, it is shown in equation (101):

[0268] (101)

[0269] The quadratic constraint (79) is approximated by two square constraints (102) at a 45-degree angle:

[0270] (102)

[0271] For nonlinear terms (97), use fixed values. The estimated power required by the compressor, after linearization, is shown in equation (103):

[0272] (103)

[0273] Example 2

[0274] This embodiment proposes a power-gas integrated energy system recovery decision-making device that considers the fusion of the charged domain during the recovery process, including:

[0275] The constraint module is used to construct constraints for sequential recovery problems, equipment constraints, and system operation constraints;

[0276] The objective function construction module for the recovery decision model is used to establish the decision model and objective function of the integrated electric-gas energy system with the goal of maximizing the net recovery benefit of the integrated electric-gas energy system.

[0277] The recovery decision scheme determination module is used to linearize the nonlinear constraints in the constraint module and optimize the objective function of the recovery decision model to determine the recovery decision scheme of the integrated electric-gas energy system.

[0278] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0279] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0280] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0281] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0282] The embodiments of the present invention have been described above with reference to the accompanying drawings. However, the present invention is not limited to the specific embodiments described above. The specific embodiments described above are merely illustrative and not restrictive. Those skilled in the art can make many other forms under the guidance of the present invention without departing from the spirit and scope of the claims. All of these forms are within the protection scope of the present invention.

Claims

1. A method for decision-making in the recovery of an integrated electric-gas energy system, considering the fusion of charged domains during the recovery process, characterized in that, include: Based on the energized domain fusion strategy and the number of steps during line operation, a sequential recovery problem constraint that satisfies the electro-pneumatic coupling constraint is constructed. The energized domain fusion strategy is used for step-by-step collaborative recovery of the power outage distribution network. The number of steps during line operation refers to the number of line switching operations that must be performed to energize critical loads. Based on the physical characteristics of the integrated electric-gas energy system and its equipment, equipment constraints and system operation constraints are constructed. Based on the constraints of the sequential recovery problem, equipment constraints, and system operation constraints, a recovery decision model for the integrated electric-gas energy system is established with the objective of maximizing the net recovery benefit of the integrated electric-gas energy system. The constraints of the sequential recovery problem, equipment constraints, and system operation constraints are linearized. Based on the linearized constraints, the objective function of the electric-gas integrated energy system recovery decision model is optimized and solved to determine the electric-gas integrated energy system recovery decision scheme.

2. The electric-gas integrated energy system restoration decision-making method considering the fusion of charged domains during the restoration process, as described in claim 1, is characterized in that... The sequential recovery problem constraint introduces an initial source node, which is defined as the first controllable power source node to be put into operation to form an energy recovery area during the recovery process of the integrated electric-gas energy system. It is used to select the best location based on competition and merge the neighboring areas based on cooperation.

3. The electric-gas integrated energy system restoration decision-making method considering the fusion of charged domains during the restoration process, as described in claim 2, is characterized in that... The integrated electric-gas energy system includes a coupled power grid system and a gas grid system. The power grid system includes a fixed power source, a gas turbine, and several power grid nodes. The gas grid system includes a fixed gas source, an electric-driven compressor, and several gas grid nodes. The power grid nodes supply power to the electric-driven compressor to support its operation, and the gas grid nodes supply gas to the gas turbine to support its power generation and supply the power grid node demand.

4. The electric-gas integrated energy system restoration decision-making method considering the fusion of charged domains during the restoration process, as described in claim 2, is characterized in that... The constraints of the sequential recovery problem include: （1） (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) in, A collection of branches in a power distribution network; It is the set of electrical nodes in the distribution network that have self-starting capability; For the power supply nodes of the distribution network; For the set of non-power source nodes in the distribution network; For the set of SC power supply nodes in the distribution network; This refers to the set of NSC power supply nodes in the distribution network that lack self-starting capability. It is a set of FSC power supply nodes in the distribution network that are fully self-starting. Let PSC be the set of power supply nodes in the distribution network that do not have full self-starting capability. The relationships between the sets are as follows: (20) The recovery state of branch ij at step t is represented by 1, indicating that it has been recovered and 0 indicates that it has not been recovered. The recovery state of step node i at time t, where 1 indicates that it has recovered and 0 indicates that it has not recovered; For time step t, whether SC power node i is a primary power node, 1 indicates that it is a primary power node, and 0 indicates that it is not a primary power node; and These are the sets of injection and outflow branches for electrical node i, respectively. and These are the virtual flows of the injection and outflow branches of electrical node i, respectively; This represents the virtual output of the SC power node; H is a positive number. This represents the recovery state of branch ij at step t+1; The recovery state of step node i at time t+1; T is the total number of observation steps; Let be the number of time steps elapsed from when branch ij is energized at the terminal node until restoration is complete. The recovery state of step node j at time t, where 1 indicates that it has recovered and 0 indicates that it has not recovered; The number of operation steps required to restore branch ij is the sum of the operation steps of all line switches on the branch; The recovery status of the step circuit branch ij at time T is 1, indicating that it has recovered and 0 indicates that it has not recovered. Is power node i of SC at time step t+1 the initial power node? The start / stop status of power supply m under step node i at time t, 1 indicates that it has been started, and 0 indicates that it has not been started; The set of fixed power supplies configured under node i; The number of operation steps required to start the fixed power supply m at node i; Let i be the set of FSC power sources at node i; Let i be the set of NSC power sources at node i.

5. The electric-gas integrated energy system restoration decision-making method considering the fusion of charged domains during the restoration process, as described in claim 1, is characterized in that... The equipment constraints include gas turbine allocation constraints, mobile emergency power supply (MEPS) allocation constraints, fixed power supply output constraints, fixed gas source constraints, electric compressor start / stop status constraints, mobile emergency power supply (MEPS) output constraints, mobile emergency gas source allocation constraints, electrochemical energy storage constraints, and gas storage tank constraints. The equipment includes a gas turbine, mobile emergency power supply, fixed power supply, fixed gas source, electric compressor, mobile emergency gas source, electrochemical energy storage, and gas storage tank.

6. The electric-gas integrated energy system restoration decision-making method considering the fusion of charged domains during the restoration process, as described in claim 1, is characterized in that... The system operation constraints satisfy the load recovery state constraints: (21) (22) (23) in, and These represent the recovery states of node i's load at time step t and time step t+1, respectively. and These represent the recovery states of node g's load at time step t and time step t+1, respectively. Let i be the recovery state of the step node at time t. and These are the sets of load nodes for the power distribution network and the gas distribution network, respectively.

7. The electric-gas integrated energy system restoration decision-making method considering the fusion of charged domains during the restoration process, as described in claim 6, is characterized in that... The system operation constraints also include power flow constraints of the distribution network and dynamic power flow constraints of the gas distribution network. The power flow constraints of the distribution network include: (24) (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) in, For the gas node g, which is connected to the electro-pneumatic coupling device at electrical node i, A collection of branches in a power distribution network; It is the set of electrical nodes in the distribution network that have self-starting capability; and These are the sets of injection and outflow branches for electrical node i, respectively. and These are the active power and reactive power of branch ij at time step t, respectively. and These represent the active power and reactive power flowing out of step node i to the outside at time t, respectively. and These represent the active power and reactive power flowing into step node i at time t, respectively. and Let be the resistance and reactance of branch ij, respectively; Let be the voltage at step node i at time t; This serves as the reference voltage in the power system. As an auxiliary variable; Let H represent the recovery state of branch ij at step t; H is a positive number. Let be the maximum apparent power that can pass through branch ij; and These represent the maximum and minimum allowable voltages for electrical node i, respectively; The power demand of step-electric node i at time t after the load participates in IDR, excluding power supply and gas supply equipment; This refers to the set of NSC power supply nodes in the distribution network that lack self-starting capability. Let m be the start / stop state of power supply m at step node i at time t. The active power required for the N-SC power supply m at node i to start; Let i be the electrochemical energy storage collection at electrical node i. Let m be the charging state of the electrochemical energy storage at step node i at time t. The active power required to charge the electrochemical energy storage m at the step node i at time t; The fixed electrically driven gas source set at gas node g. This represents the start / stop state of the fixed gas source n at the gas stepping node g at time t. Let n be the power consumed by the fixed electric drive air source n at step node i at time t. This refers to the collection of P2G devices at gas node g. Let t be the power consumed by P2G device n at step node i; The start / stop state of the electrically driven compressor at step gas node g at time t is given. Let t be the power consumed by the compressor driven by step node i when step node i is powered on; Let i be the reactive power demand factor of electrical node i; and These are the active power and reactive power of branch ki at time step t, respectively. Let be the voltage at step node j at time t; Let be the maximum active power that can pass through branch ij; Let be the maximum reactive power that can pass through branch ij; The dynamic flow constraints of the gas distribution network include dynamic flow constraints of natural gas, node flow balance and node pressure constraints, flow and pressure constraints in pipelines, and node gas pressure constraints. The dynamic power flow constraint for natural gas is: (35) (36) in, A collection of pipes for the gas distribution network; Let gh be the cross-sectional area of ​​the pipe; For time step; This refers to the length of the pipe. The average flow velocity of natural gas in pipeline gh; Where gh is the inner diameter of the pipe; and These represent the flow rates at the end and beginning of the pipeline at time step t, respectively. and These are the end pressure and beginning pressure of the pipeline at time t, respectively. This refers to the pipe friction coefficient; The speed of sound in a gas; The node flow balancing and node pressure constraints are as follows: (37) (38) (39) (40) in, Let t be the flow rate of gas injected into node g from the outside; Let g be the gas flow rate of the object to be restored at the gas node g at time t; and These are the sets of injection and outflow branches for gas node g, respectively. Let fg be the flow rate at the end of pipe at time step t; The recovery state of the load at node g at time step t; To calculate the flow demand of the gas node g at time t after adjustment via IDR, excluding power supply and gas supply equipment. Let be the set of GT at electrical node i; Let be the gas flow rate of the gas turbine m connected to the gas stepping node g at time t; Let g be the set of GT at the gas node g; The gas storage state of gas storage tank n at gas node g at time t; Let n be the gas storage flow rate at gas storage tank n at gas step node g at time t; Let g be the air pressure at step node g at time t; A collection of air nodes equipped with electrically driven compressors; This refers to the set of gas nodes in the gas distribution network. Let fg be the flow rate at the end of pipe at time step t; , where are the flow rates at the beginning of pipe gh at time step t; The flow and pressure constraints in the pipeline are as follows: (41) (42) (43) (44) in, The maximum permissible mass flow rate for pipeline gh; and These are the upper and lower limits of the allowable pressure for the pipeline gh, respectively. The node air pressure constraint is: (45) (46) (47) (48) in, Electrical node i is connected to the electro-pneumatic coupling device at gas node g. and These are the minimum and maximum air pressures allowed at node g, respectively; Let g be the air pressure at the stepping node; The power consumed by the compressor driven by the pipeline gh at time step t; K1, K2, and K3 are the outlet pressure setpoint of the electrically driven compressor at node g; K1, K2, and K3 are the compressor equation coefficients and exponents. Let t be the power consumed by the compressor driven by step node i at time t.

8. The method for decision-making on the restoration of an integrated electric-gas energy system considering the fusion of charged domains during the restoration process, as described in claim 1, is characterized in that... The objective function of the electric-gas integrated energy system recovery decision model is: (49) in, The net benefit of obtaining electricity or gas supply in advance for an integrated electric-gas energy system is calculated using the following formula: (50) in, and These are the sets of load nodes for the power distribution network and the gas distribution network, respectively. Let be the unit recovery revenue of load node i at time step t; This represents the recovery status of the load at node i at time step t, where 1 indicates recovery and 0 indicates no recovery. The electrical energy to be restored at node i at time step t, excluding power supply and gas supply equipment; Let g be the unit recovery revenue of the gas load node g at time t; This represents the recovery status of the load at node g at time step t, where 1 indicates recovery and 0 indicates no recovery. Let g be the volume of natural gas to be restored at time step t, excluding power and gas supply equipment.

9. The electric-gas integrated energy system restoration decision-making method considering the fusion of charged domains during the restoration process, as described in claim 1, is characterized in that... The linearized sequential recovery problem constraints, device constraints, and system operation constraints include: For the nonlinear terms appearing in the constraint system, equipment constraints, and system operation constraints of the sequence recovery problem, auxiliary variable constraints are introduced for linearization, transforming the original mixed integer nonlinear programming problem into a mixed integer linear programming problem.

10. A power-electric integrated energy system recovery decision-making device considering the fusion of charged domains during the recovery process, comprising: The constraint module is used to construct constraints for sequential recovery problems, equipment constraints, and system operation constraints; The objective function construction module for the recovery decision model is used to establish the decision model and objective function of the integrated electric-gas energy system with the goal of maximizing the net recovery benefit of the integrated electric-gas energy system. The recovery decision scheme determination module is used to linearize the nonlinear constraints in the constraint module and optimize the objective function of the recovery decision model to determine the recovery decision scheme of the integrated electric-gas energy system.

Images