An Energy Supply Restoration Coordination Method and System for an Electrical Integrated Energy System with Mobile Emergency Resources

Abstract

The present invention discloses a power supply restoration coordination method and system for an electrical integrated energy system containing mobile emergency resources. The method conducts state identification on the electrical integrated energy system to identify the electrical loads to be restored, the gas loads to be restored, and the set of gas turbines to be restored; makes an optimal scheduling decision for mobile emergency resources to determine the transportation paths, transportation costs, and corresponding connection points of each emergency mobile resource; optimizes the gas-to-electricity flow to determine the switch states of each gas turbine, the optimal gas supply system or quasi-gas supply system, the gas acquisition paths, and the gas acquisition flows; based on the above optimization results, conducts parallel partition division of the power distribution system and the gas distribution system to determine the loads to be restored in each partition; based on the partition results of the power distribution system and the gas distribution system, respectively optimize the restoration plans independently for each partition to determine the optimal restoration plan for each partition. The present invention is significantly helpful for effectively improving the utilization rate of energy supply resources, enhancing the restoration efficiency, and obtaining a large net restoration benefit.

Description

Technical Field

[0001] The present invention relates to a method and system for coordinating the energy supply restoration of an electrical integrated energy system containing mobile emergency resources, and belongs to the technical field of energy regulation. Background Art

[0002] The contradiction between the increasing energy demand and environmental ecological protection worldwide is becoming increasingly prominent. Developing clean energy and multi-energy complementarity to improve energy utilization efficiency has become the only way for energy transformation in countries around the world. Natural gas is rich in reserves, easy to store, highly efficient and environmentally friendly, and has a high calorific value and low production cost, making it an ideal source of urban gas. Therefore, at the user side, a gas supply system of a considerable scale and a power supply system together provide options for residents' terminal energy use. Gas-to-Power (G2P) is achieved through power generation equipment such as Gas-fired Turbines (GTs) to provide clean energy power generation. By the end of 2021, the installed capacity of gas-fired power in China exceeded 1×10 8 kW, and the installed capacity of gas-fired power accounted for approximately 4.5% of the total installed power generation capacity in the country. Power-to-Gas (P2G) technology can convert electrical energy into natural gas for storage and transportation, providing guarantee for the consumption of intermittent renewable energy such as solar energy and wind energy. With the continuous development of the natural gas system and the power system, the connection between the two is becoming increasingly close, which is manifested as the gradual deep coupling of the Power Distribution System (PDS) and the Gas Distribution System (GDS) at the user side.

[0003] Due to the high equipment investment and operating costs, and the fact that the electricity price is higher than the natural gas price at the present stage, the proportion of Power-to-Gas equipment is not large, and the electricity-gas integrated energy system considering Gas-to-Power is more common. With the large-scale grid connection of gas-fired power generation, the gas supply blockage in the natural gas system will seriously affect the safe and reliable power supply of the power system. How to allocate gas load under the limited gas supply resources in the natural gas system, and whether to supply the power system to generate electricity need to be dynamically decided according to the actual situation. Conversely, the change in the demand for natural gas caused by the power supply blockage in the power system also requires readjusting the flow distribution of the natural gas system. Therefore, the energy supply restoration of the IEGDS after outage requires coordinated decision-making. At present, the delay characteristics of the energy supply system are often ignored during the restoration of the electricity-gas integrated energy system, and the structural integrity of the natural gas system, the gas load restoration demand, and the gas source restoration ability are rarely considered.

[0004] As a supplement to the fixed power supply and gas source, the Mobile Emergency Power Source (MEPS) and the Mobile Emergency Gas Source (MEGS) can be flexibly connected to the system, playing the role of quickly supplying energy after power outage or ensuring reliable energy supply for the system. At present, the relevant research mainly focuses on considering the power supply restoration decision of MEPS. There is little research on considering both MEPS and MEGS for the restoration decision of the integrated energy system. In order to improve the practicality of the restoration decision, it is necessary to consider the coordinated decision-making of MEPS and MEGS. Summary of the Invention

[0005] The technical problem to be solved by the present invention is to overcome the defects of the prior art and provide a method and system for coordinated energy supply restoration of an electrical integrated energy system containing mobile emergency resources, aiming at maximizing the net income during the restoration process of the electrical integrated energy system, and coordinately optimizing the dispatching of each component in the electrical integrated energy system.

[0006] To solve the above technical problems, the technical solution adopted by the present invention is as follows:

[0007] The present invention provides a method for coordinated energy supply restoration of an electrical integrated energy system containing mobile emergency resources, including:

[0008] Identifying the state of the electrical integrated energy system containing mobile emergency resources, analyzing the states of the loads to be restored, energy supply resources, branches, and switch equipment in the system, obtaining the set of electrical loads to be restored, the gas loads to be restored, and the set of gas turbines to be restored, and determining the energy supply scope of each energy supply system; the mobile emergency resources include the Mobile Emergency Power Source (MEPS) and the Mobile Emergency Gas Source (MEGS); the energy supply systems include the power supply system on the distribution side and the gas supply system on the gas distribution side;

[0009] On the distribution side, optimizing the configuration of the partition strategy for the electrical loads to be restored, specifically: dividing partitions for the fixed power supply, and restoring the electrical loads to be restored within each partition based on the available power of the fixed power supply until the fixed power supply has no available power; taking the maximum net income of electrical load restoration as the goal, sequentially connecting the Mobile Emergency Power Source (MEPS) to each candidate connection point in the power grid and updating the partition until all the electrical loads to be restored are restored or the MEPS has no available power, obtaining the optimal partition strategy on the distribution side, the transportation paths, transportation costs, and corresponding connection points of each Mobile Emergency Power Source (MEPS);

[0010] On the gas distribution side, optimize the configuration of the zoning strategy for the gas load to be restored, specifically: calculate the initial revenue and required flow rate of the gas turbine, divide the fixed gas sources into zones, and sequentially connect the mobile emergency gas sources (MEGS) to each candidate connection point in the gas network with the goal of maximizing the net revenue of gas load restoration and update the zones until all the gas load to be restored is restored or the MEGS has no available power, so as to obtain the optimal zoning strategy on the gas distribution side, the transportation paths, transportation costs and corresponding connection points of each mobile emergency gas source (MEGS);

[0011] Optimize the G2P flow configuration of the gas turbine in the electrical coupling component, and with the goal of maximizing the net revenue of load restoration, allocate the flow rate supplied to the gas load to be restored and the flow rate supplied to the electrical load to be restored through the gas turbine, so as to obtain the switching states, gas supply systems, gas supply volumes and gas supply paths of each gas turbine;

[0012] According to the optimal zoning strategy on the power distribution side and the optimal zoning strategy on the gas distribution side, as well as the flow configuration of the gas turbine, perform parallel zoning on the power distribution side and the gas distribution side to determine the loads to be restored in each zone;

[0013] Perform load restoration for each zone on the power distribution side and the gas distribution side with the goal of maximizing the net revenue during the restoration process of the electrical integrated energy system, so as to obtain the optimal restoration plan for each load. The restoration plan includes the restoration status of each step of the load, the status of the energy supply path, the energy supply volume, the connection points and transportation paths of the mobile emergency resources;

[0014] Furthermore, perform state identification on the electrical integrated energy system containing mobile emergency resources, analyze the states of the loads to be restored, energy supply resources, branches and switch equipment in the system, obtain the set of electrical loads to be restored, the gas load to be restored and the set of gas turbines to be restored, and determine the energy supply scope of each energy supply system, including:

[0015] According to the topological structure of the electrical integrated energy system and the fault conditions of the equipment, identify the isolated power outage areas and non-isolated power outage areas, and then obtain the set of electrical loads to be restored The gas load to be restored Set and the set of gas turbines to be restored Ω GT ; The isolated power outage area refers to the power outage area that has no recoverable path to all energy supply systems or quasi-energy supply systems; The quasi-energy supply system refers to the system of a source that is in a power outage state but has self-starting ability; The quasi-energy supply system includes the quasi-power supply system on the power distribution side and the quasi-gas supply system on the gas distribution side;

[0016] Determine the power supply scope of the power supply system or quasi-power supply system according to the available power of the power supply system or quasi-power supply system and the voltage level on the power transmission path;

[0017] Determine the gas supply range of the gas supply system or quasi-gas supply system according to the available flow rate of the gas supply system or quasi-gas supply system and the natural gas pressure level on the gas transmission path.

[0018] Further, partition the fixed power supply, and within each partition, restore the electrical load to be restored based on the available power of the fixed power supply until there is no available power in the fixed power supply, including:

[0019] Step S11: Based on the power supply systems and quasi-power supply systems with available power, initially partition the distribution side, divide each power supply system or quasi-power supply system into an initial partition, and uniquely partition the electrical load nodes to be restored into the partition that can supply power to them and is the closest;

[0020] Step S12: Each partition preferentially supplies power to the electrical load to be restored with the highest net restoration benefit, and conducts simulation and deduction. Delete the electrical load that can be restored through deduction from and update Traverse the electrical load to be restored in this way until all the electrical load to be restored is restored or there is no available power in both the power supply system and the quasi-power supply system.

[0021] Further, taking the maximum net restoration benefit of the electrical load restoration as the goal, sequentially connect the mobile emergency power supply MEPS to each candidate connection point in the power grid and update the partition until all the electrical load to be restored is restored or the MEPS has no available power, including:

[0022] Step S13: If is not an empty set, then determine whether the set Ω MEPSJ of the candidate connection points of the power grid MEPS and the set Ω MEPS of the MEPS is an empty set. If any one of them is an empty set, end; otherwise, sort the elements in Ω MEPS in descending order of rated power, and go to step S14;

[0023] Step S14: Connect the first element Ω MEPS of Ω MEPS (1) to a certain candidate connection point in Ω MEPSJ , and update the distribution side partition with Ω MEPS (1) as a power supply system;

[0024] Calculate the net restoration benefit of the load and the shortest transportation path after the access of Ω MEPS (1) according to the following formula, and traverse all the candidate connection points in Ω MEPSJ to determine the maximum net restoration benefit of the load;

[0025] I M =β M I M,J -C M,J ;

[0026] Among them, I M is the net load recovery benefit obtained after the mobile emergency resource accesses the candidate connection point; β M is the conversion coefficient of the recovery benefit of the mobile emergency resource candidate connection point, I M,J , C M,J are respectively the recovery benefit and cost of the mobile emergency resource supplying energy to the load to be recovered at the candidate connection point;

[0027] If the maximum net load recovery benefit is less than or equal to 0, then Ω MEPS (1) is deleted from Ω MEPS , and go to step S13; otherwise, determine the candidate connection point with the maximum net load recovery benefit obtained as the optimal access point of Ω MEPS (1), and go to step S15;

[0028] Step S15: Connect Ω MEPS (1) to the determined optimal access point to obtain the connection point, the shortest transportation path and the transportation cost of the mobile emergency power supply MEPS;

[0029] Regard Ω MEPS (1) as a power supply system, update the distribution side partition, and delete Ω MEPS (1) from Ω MEPS , and go to step S12.

[0030] Furthermore, calculating the initial benefit and required flow rate of the gas turbine, partitioning the fixed gas source, and sequentially connecting the mobile emergency gas source MEGS to each candidate connection point in the gas network with the goal of maximizing the net gas load recovery benefit and updating the partition until all the gas loads to be recovered are restored or MEGS has no available power, including:

[0031] Step S21: For the elements in Ω GT , regard one gas turbine as a power supply system, uniquely partition the nodes in into the gas turbine partition that can supply power to it and is the closest in distance, take the maximum benefit that can be obtained from the rated power of the gas turbine in the partition as the initial benefit on the distribution side, and take the gas flow rate required by the gas turbine as its gas load demand on the gas distribution side; delete the GT with an initial benefit equal to zero from Ω GT , and update

[0032] Step S22: Divide each gas supply system or quasi - gas supply system into an initial partition, and uniquely partition the nodes in into the partition that can supply gas to it and is the closest in distance;

[0033] Step S23: Judge whether the set of candidate connection points Ω MEGSJ of the gas network MEGS is the same as the set Ω of MEGSMEGS Whether it is an empty set. If any of them is an empty set, end; otherwise, sort the elements in Ω MEGS in descending order according to the rated flow rate, and go to step S24;

[0034] Step S24: Take the first element Ω of Ω MEGS (1) and connect it to a certain candidate connection point, and use Ω MEGS (1) as a gas supply system to update the distribution side partition, calculate the net benefit of load restoration and the shortest transportation path, and traverse all candidate connection points in Ω MEGSJ to determine the maximum net benefit of load restoration; MEGS (1) If the maximum net benefit of load restoration is less than or equal to 0, then delete Ω MEGSJ (1) from Ω, and go to step S23; otherwise, determine the candidate connection point with the maximum net benefit of load restoration as the optimal access point of Ω

[0035] (1), and go to step S25; MEGS (1) and delete it from Ω MEGS MEGS MEGS (1), and go to step S25;

[0036] Step S25: Connect Ω MEGS (1) to the determined optimal access point to obtain the connection point, the shortest transportation path and the transportation cost of the mobile emergency gas source MEGS;

[0037] Use Ω MEGS (1) as a gas supply system to update the distribution side partition accordingly, delete Ω MEGS (1) from Ω MEGS and delete the determined optimal access point from Ω MEGSJ ;

[0038] Step S26: Based on the updated partition, re-optimize and simulate the restoration of the gas load to be restored in each partition, and then delete the gas load that can be restored after the simulation from and update Ω LG ;

[0039] Judge whether it is an empty set. If it is an empty set, end; otherwise, go to step S23.

[0040] Furthermore, the G2P flow optimization configuration of the electrical coupling component gas turbine includes:

[0041] Step S31: Take the maximum benefit that can be obtained by the rated power of the gas turbine in the partition as its initial benefit on the distribution side, update Ω in descending order according to the initial benefit GT , and judge whether Ω GT is an empty set. If Ω GTIf it is an empty set, end; otherwise, go to step S32;

[0042] Step S32: Determine whether there is still a load to be restored on the distribution side; if it is an empty set, end; otherwise, go to step S33;

[0043] Step S33: If there is a situation where a gas turbine in Ω GT cannot obtain stable gas supply, delete it from Ω GT and determine whether Ω GT is an empty set. If Ω GT is an empty set, end; otherwise, go to step S34;

[0044] Step S34: Obtain the first element Ω GT of Ω GT (1) The path with the shortest gas acquisition delay from the gas supply system or the quasi-gas supply system The maximum available gas flow where, is the available gas flow of the gas supply system or the quasi-gas supply system, is the rated gas flow of Ω GT (1), and go to step S35;

[0045] Step S35: Assume that the maximum net benefits obtained by the gas acquisition flow of Ω GT (1) in its distribution side partition and gas distribution side partition are respectively and Optimize to obtain the optimal gas acquisition flow F of Ω GT (1) that maximizes GT(1) ;

[0046] If then delete Ω GT (1) from Ω GT and go to step S31;

[0047] Otherwise, allocate the gas supply system or the quasi-gas supply system to plan to supply gas to Ω GT (1) Regard Ω GT (1) as a power supply system, update the distribution side partition, optimize the load restoration plan and simulate and deduce, update and the partitions and initial benefits of other gas turbines on the distribution side, update the remaining available gas flow on the gas distribution side, update the gas distribution side partition and Delete Ω GT (1) from Ω GT and go to step S31.

[0048] Furthermore, the parallel partitioning of the distribution side and the gas distribution side includes:

[0049] If there is no gas intake flow in all gas turbines, on the power distribution side, determine the partitions of the power supply system and the quasi-power supply system other than the gas turbines according to the method of partitioning by the power distribution side; on the gas distribution side, partition the gas supply system and the quasi-gas supply system according to the method of partitioning by the gas distribution side.

[0050] If there are gas turbines with gas intake flow, on the power distribution side, according to the method of partitioning by the power distribution side, uniquely partition the remaining unpartitioned electrical load nodes to the partition of the gas turbine that can supply power to them, is the closest in distance, and has gas intake flow; on the gas distribution side, for the gas supply system and the quasi-gas supply system that plan to supply gas to the gas turbines, partition the nodes related to the gas intake path of the gas turbines into their partitions, and uniquely partition the remaining unpartitioned nodes to the partition of the gas supply system or the quasi-gas supply system that can supply gas to them and has the shortest gas supply delay.

[0051] Furthermore, perform load restoration with the goal of maximizing the net restoration benefit during the restoration process of the electrical integrated energy system for each partition on the power distribution side and the gas distribution side, and obtain the optimal restoration plan for each load, including:

[0052] Step S41: In each partition on the power distribution side and the gas distribution side, according to the power outage scenario, arrange the net restoration benefits of the targets to be restored in descending order to generate a sequence Ω of the targets to be restored in the partition. TR ; The targets to be restored are the electrical loads to be restored and the gas loads to be restored.

[0053] Step S42: Select the first target to be restored Ω TR (1) from Ω TR , and evaluate the restoration benefit and restoration cost of restoring Ω TR (1) based on the following formula, and simulate and deduce its power outage scenario:

[0054]

[0055]

[0056]

[0057]

[0058] C T =∑χ T (t)u T (t);

[0059] Among them, and are respectively the electrical load at the k P th step and the gas load restoration benefit at the k G th step; and are respectively the kP The power supply at step k and the air supply at step k G The cost of supplying energy to the load; C T The transportation costs of MEPS and MEGS; t e Is the preset evaluation end time of the system, set as the larger value after the estimated execution completion time of the power distribution side and gas distribution side restoration plans; For step k P The power load restoration time at step k; Are respectively the unit restoration benefit and restoration volume of the power load restored at step k at time t; P Are respectively the unit power cost and power supply volume of the fixed power supply at step k providing power to the power load at time t; P Are respectively the unit power cost and power supply volume of MEPS at step k providing power to the power load at time t; P For step k G The gas load restoration time at step k; Are respectively the unit restoration benefit and restoration volume of the gas load restored at step k at time t; G Are respectively the unit gas cost and gas supply volume of the fixed gas source at step k providing gas flow to the gas load at time t; G Are respectively the unit gas cost and gas supply volume of MEGS at step k providing gas flow to the gas load at time t; u G T ψ(t) is the state of a certain section of the road passed by the mobile emergency resource at time t during transportation, 1 means passed, 0 means not passed; χ T φ(t) is the transportation cost of the mobile emergency resource passing through this section of the road at time t;

[0060] Step S43: Delete Ω TR (1) from Ω TR , update Ω TR , if the breadth deduction is completed in the current scenario or Ω TR is an empty set, then go to step S44; otherwise go to step S42 to generate a restoration plan for the next restoration target in the current scenario;

[0061] Step S44: If the deduction depth has not reached the set value, then switch to the simulation stop - energy scenario of the next step and go to step S41; otherwise go to step S45;

[0062] ​​​​​​Step S45: Calculate the net recovery benefit per unit capacity of all multi-step recovery plans for the current partition. With the goal of maximizing the net recovery benefit per unit capacity of all multi-step recovery plans in the partition during the recovery process of the integrated electrical energy system, determine the optimal recovery plan. Among them, define the ratio of the net recovery benefit of the recovery plan to its recovery capacity as the net recovery benefit per unit capacity.

[0063] The objective function for maximizing the net benefit of all multi-step recovery plans in the partition during the recovery process of the integrated electrical energy system is as follows:

[0064]

[0065] Among them, and are the total number of steps of the recovery plans on the power distribution side and the gas distribution side within the same decision period, respectively; k P and k G are the step numbers in the PDS and GDS recovery plans, respectively.

[0066] Furthermore, the objective function needs to satisfy the following constraint conditions:

[0067] A. Power distribution side constraints, including steady-state constraints and topological structure constraints:

[0068] A1. Steady-state constraints:

[0069]

[0070]

[0071]

[0072]

[0073]

[0074] U i,min ≤U i (t)≤U i,max , i∈Ω B ;

[0075] I ij,min ≤I ij (t)≤I ij,max , i,j∈Ω B ;

[0076]

[0077] Among them, Ω B is the set of energized nodes in the power supply area; Ω MEPS,i is the set of MEPSs connected to grid node i; Pi grid (t) and are the active and reactive powers injected by the external network at node i at time t; P i DG (t) and are the active and reactive powers injected by the distributed power source at node i at time t; P i PM (t) and are the active and reactive powers injected by the MEPS at node i at time t; P i L (t) and are the active and reactive powers of the load at node i at time t; U i (t), U j (t) are the voltage values of nodes i and j at time t; G ij and B ij are the real and imaginary parts of the admittance of line ij; θ ij (t) is the phase difference between the voltages of nodes i and j at time t; are the upper and lower limits of the active power of the distributed power source at node i; are the rated power and the active power output at time t of the mobile emergency power source c, c ∈ Ω MEPS,i ; is the energy supply time of the mobile emergency power source c; U i,max and U i,min are the upper and lower limits of the voltage of node i; I ij (t), I ij,max and I ij,min are the current of branch ij and its upper and lower limits at time t; are the transmission capacity of branch ij and its upper limit at time t.

[0078] A2. Transient constraint:

[0079] ΔP k,step ≤ ζP kN ;

[0080] where ΔP k,step is the active power capacity of the single-step restoration of the electrical load in the k-th step; ζ is the proportionality coefficient of the allowable single-step restoration capacity of the electrical load; P kN is the rated capacity in the energized area at the initial moment of the k-th step;

[0081] A3. Topological structure constraint:

[0082] g ∈ G;

[0083] where g is the topological structure of the energized area after the distribution network reconstruction; G is the set of all radial network topologies on the distribution side;

[0084] B. Gas distribution side constraints, including steady-state constraints and transient constraints:

[0085] B1. Steady-state constraints:

[0086]

[0087]

[0088]

[0089]

[0090]

[0091] Π m,min ≤ Π m (t) ≤ Π m,max , m ∈ Ω G ;

[0092] F mn,min ≤ F mn (t) ≤ F mn,max , m, n ∈ Ω G ;

[0093] where Ω G is the set of gas acquisition nodes in the gas supply area; are the gas flow rates injected by the external network and the gas source at node m at time t, respectively; is the gas flow rate injected by MEGS at node m at time t; is the gas flow rate of the load at node m at time t; F mn (t) is the pipe segment flow rate from node m to node n at time t; are the upper and lower limits of the gas injection flow rate at node m, respectively; the mobile emergency gas source d is the MEGS connected to node m of the gas network; are the rated flow rate and the output flow rate at time t of the mobile emergency gas source d, respectively; is the energy supply time of MEGS at node m; Π m (t), Π m,max , Π m,min are the pressure value at node m at time t and its upper and lower limits, respectively; F mn,max , F mn,min are the upper and lower limits of the flow rate of pipe segment mn, respectively;

[0094] B2. Transient constraints:

[0095]

[0096]

[0097] Among them, Ω GT is the set of gas turbines to be restored; are the natural gas pressure value injected into gas turbine a at time t and its upper and lower limits, respectively; are the natural gas flow value injected into gas turbine a at time t and its upper and lower limits, respectively;

[0098] C. Mobile emergency resource transportation constraints, including:

[0099] C1. Mobile emergency power supply constraints:

[0100]

[0101] Among them, Ω MEPSJ is the set of candidate connection points of the power grid MEPS; is the number of MEPS units connected to node i; N MEPS is the total number of MEPS units;

[0102] C2. Mobile emergency gas source constraints:

[0103]

[0104] Among them, Ω MEGSJ is the set of candidate connection points of the gas network MEGS; is the number of MEGS units connected to node m; N MEGS is the total number of MEGS units;

[0105] C3. Transportation path constraints:

[0106] l ∈ L;

[0107] Among them, l is the transportation path of the mobile emergency resource from its storage location to the candidate connection point; L is the set of all radial paths of the transportation network;

[0108] C4. Transportation time constraints:

[0109]

[0110] Among them, t 0 is the starting time of restoration; is the transportation time of mobile emergency power supply c passing through section pq; l c is the shortest transportation path of mobile emergency power supply c from its storage location to the candidate connection point; is the startup time consumption of mobile emergency power supply c;

[0111] The transportation time constraint of MEGS is the same as that of MEPS.

[0112] The present invention also provides an energy supply restoration coordination system for an electrical integrated energy system with mobile emergency resources, which is used to implement the energy supply restoration coordination method for the electrical integrated energy system with mobile emergency resources as described above. The system includes:

[0113] An initialization module, which is used to identify the state of the electrical integrated energy system with mobile emergency resources, analyze the states of the system's load to be restored, energy supply resources, branches, and switch equipment, obtain the set of electrical loads to be restored, the gas load to be restored, and the set of gas turbines to be restored, and determine the energy supply scope of each energy supply system; the mobile emergency resources include mobile emergency power sources (MEPS) and mobile emergency gas sources (MEGS); the energy supply systems include the power supply system on the power distribution side and the gas supply system on the gas distribution side;

[0114] A mobile emergency power source coordination module, which is used to optimize the configuration of the partition strategy for the electrical load to be restored on the power distribution side. Specifically: divide partitions for the fixed power sources, and restore the electrical load to be restored within each partition based on the available power of the fixed power sources until there is no available power for the fixed power sources; with the goal of maximizing the net benefit of electrical load restoration, connect the mobile emergency power source (MEPS) to each candidate connection point in the power grid in sequence and update the partitions until all the electrical load to be restored is restored or the MEPS has no available power, obtaining the optimal partition strategy on the power distribution side, the transportation paths, transportation costs, and corresponding connection points of each mobile emergency power source (MEPS);

[0115] A mobile emergency gas source coordination module, which is used to optimize the configuration of the partition strategy for the gas load to be restored on the gas distribution side. Specifically: calculate the initial benefit and required flow rate of the gas turbine, divide partitions for the fixed gas sources, and with the goal of maximizing the net benefit of gas load restoration, connect the mobile emergency gas source (MEGS) to each candidate connection point in the gas network in sequence and update the partitions until all the gas load to be restored is restored or the MEGS has no available power, obtaining the optimal partition strategy on the gas distribution side, the transportation paths, transportation costs, and corresponding connection points of each mobile emergency gas source (MEGS);

[0116] An electrical coupling element coordination module, which is used to optimize the G2P flow configuration of the electrical coupling element gas turbine. With the goal of maximizing the net benefit of load restoration, allocate the flow rate for supplying the gas load to be restored and the flow rate for supplying the electrical load to be restored through the gas turbine, obtaining the switch states, gas supply systems, gas supply volumes, and gas supply paths of each gas turbine;

[0117] A parallel recovery module is used to perform parallel zoning on the power distribution side and the gas distribution side according to the optimal zoning strategy on the power distribution side, the optimal zoning strategy on the gas distribution side, and the flow configuration of the gas turbine, determine the loads to be recovered in each zone; and perform load recovery on each zone of the power distribution side and the gas distribution side with the goal of maximizing the net recovery benefit during the recovery process of the electrical integrated energy system, to obtain the optimal recovery plan for each load. The recovery plan includes the recovery status of each step of the load, the status of the energy supply path, the energy supply quantity, the connection points of the mobile emergency resources, and the transportation path.

[0118] The beneficial effects achieved by the present invention:

[0119] The present invention comprehensively considers the recovery benefit and cost, and establishes an energy supply recovery optimization model for an electrical-gas integrated energy system containing emergency mobile resources; performs state identification on the system, analyzes the states of equipment such as the loads to be recovered, energy supply resources, branches, and switches in the system, and obtains the set Ω of electrical loads to be recovered LP , the gas load to be recovered and the set Ω of gas turbines to be recovered GT , and determines the energy supply scope of each energy supply system; makes an optimal scheduling decision for mobile emergency resources to determine the transportation path, transportation cost, and their corresponding connection points of each emergency mobile resource; makes an optimal decision on the gas-to-electricity flow to determine the optimal gas supply system (quasi-gas supply system), gas acquisition path, and gas acquisition flow of each gas turbine; based on the results of the optimal decision on the gas-to-electricity flow, performs parallel zoning division on the power distribution system and the gas distribution system to determine the loads to be recovered in each zone; based on the zoning results of the power distribution system and the gas distribution system, respectively optimize the recovery plan for each zone to determine the optimal recovery plan for each zone. The energy supply recovery coordination optimization decision method for the electrical-gas integrated energy system containing emergency mobile resources provided by the present invention is significantly helpful for effectively improving the utilization rate of energy supply resources, enhancing the recovery efficiency, and obtaining a larger net recovery benefit. Description of the Drawings

[0120] Figure 1 It is a flowchart of an energy supply recovery coordination method for an electrical integrated energy system containing mobile emergency resources provided by an embodiment of the present invention;

[0121] Figure 2 It is a schematic diagram of the G2P flow optimization principle in an embodiment of the present invention. Detailed Embodiment

[0122] The present invention will be further described below. The following embodiments are only used to more clearly illustrate the technical solutions of the present invention and cannot be used to limit the protection scope of the present invention.

[0123] An embodiment of the present invention provides an energy supply recovery coordination method for an electrical integrated energy system containing mobile emergency resources, including:

[0124] Perform state identification on the electrical integrated energy system with mobile emergency resources, analyze the states of equipment such as the load to be restored, energy supply resources, branches, and switches in the system, and obtain the set of electrical loads to be restored. Gas load to be restored and the set Ω of gas turbines to be restored GT , and determine the energy supply scope of each energy supply system;

[0125] On the distribution side, optimize the configuration of the partition strategy for the electrical load to be restored. Specifically: divide partitions for fixed power sources, and within each partition, restore the electrical load to be restored based on the available power of the fixed power source until the fixed power source has no available power; with the goal of maximizing the net revenue of electrical load restoration, sequentially connect the mobile emergency power source MEPS to each candidate connection point in the power grid and update the partition until all the electrical load to be restored is restored or the MEPS has no available power, to obtain the optimal partition strategy on the distribution side, the transportation paths, transportation costs, and corresponding connection points of each mobile emergency power source MEPS.

[0126] On the gas distribution side, optimize the configuration of the partition strategy for the gas load to be restored. Specifically: calculate the initial revenue and required flow rate of the gas turbine, divide partitions for fixed gas sources, and with the goal of maximizing the net revenue of gas load restoration, sequentially connect the mobile emergency gas source MEGS to each candidate connection point in the gas network and update the partition until all the gas load to be restored is restored or the MEGS has no available power, to obtain the optimal partition strategy on the gas distribution side, the transportation paths, transportation costs, and corresponding connection points of each mobile emergency gas source MEGS.

[0127] Perform G2P flow optimization configuration on the electrical coupling element GT. Under the condition that the total available flow rate of the GDS gas source is determined, with the goal of maximizing the net revenue of load restoration, allocate the flow rate supplied to the gas load to be restored and the flow rate supplied to the electrical load to be restored through GT, to obtain the switch state, gas supply system, gas supply volume, and gas supply path of each GT.

[0128] According to the optimal partition strategy on the distribution side and the optimal partition strategy on the gas distribution side, as well as the flow configuration of the gas turbine, perform parallel partitioning on the distribution side and the gas distribution side to determine the loads to be restored within each partition.

[0129] Perform load restoration for each partition on the distribution side and the gas distribution side with the goal of maximizing the net revenue during the restoration process of the electrical integrated energy system, to obtain the optimal restoration plan for each load, where the restoration plan includes the restoration state of each step of the load, the state of the energy supply path, the energy supply volume, the connection points of mobile emergency resources, and the transportation paths, etc.

[0130] In the embodiments of the present invention, the mobile emergency resources include a mobile emergency power source MEPS and a mobile emergency gas source MEGS.

[0131] In the embodiment of the present invention, the state of the electric-gas integrated energy system containing mobile emergency resources is identified to determine the energy supply range, and the set of electric loads to be restored is identified. The set of gas loads to be restored and the set Ω of gas turbines GT to be restored GT , and the specific implementation process is as follows:

[0132] The electric-gas integrated energy system containing mobile emergency resources is abbreviated as MG2PIES (Mobile emergency resource included Gas-to-Power Integrated Energy System). Among them, the stable energy supply part is called the energy supply system, and the source in the power-off state but with self-starting ability is called the quasi-energy supply system, corresponding to the power supply system and the quasi-power supply system in PDS, and the gas supply system and the quasi-gas supply system in GDS respectively.

[0133] According to the topological structure of the system and the fault conditions of the equipment, the isolated power-off areas and non-isolated power-off areas are identified. Among them, the isolated power-off area refers to the power-off area that has no recoverable path to all energy supply systems or quasi-energy supply systems, and then the set of electric loads to be restored The gas loads to be restored set and the set Ω of gas turbines to be restored GT are obtained.

[0134] Analyze the available energy of each energy supply system or quasi-energy supply system, corresponding to the available power and available electricity in PDS, and the available flow and available gas volume in GDS respectively. It should be noted that the available power in PDS is the standby power of the power supply system or the predicted maximum output of the quasi-power supply system; the available flow in GDS is the standby flow of the gas supply system or the predicted maximum flow of the quasi-gas supply system.

[0135] Determine the power supply range of the power supply system or the quasi-power supply system according to factors such as available power and voltage level on the power transmission path; determine the gas supply range of the gas supply system or the quasi-gas supply system according to factors such as available flow and natural gas pressure level on the gas transmission path.

[0136] In the embodiment of the present invention, the optimal scheduling strategy of the mobile emergency power supply MEPS can be referred to Figure 1 , including:

[0137] Step S11: Set the initial partition of PDS as the power supply system and the quasi-power supply system itself with available power. Partition the power supply system and the quasi-power supply system in PDS. Considering the voltage over-limit problem based on the power supply range of the initial partition of PDS, the nodes in are uniquely partitioned into the nearest partition that can supply power to them.

[0138] In this embodiment, each power supply system or quasi-power supply system is divided into an initial partition. The available power of the partition is the available power of its power supply system or quasi-power supply system. The power load nodes to be restored are uniquely divided into the partition that can supply power to them and is the closest.

[0139] Step S12: Calculate the net recovery benefit of the power loads to be restored in each partition considering the available power of the power supply system and the quasi-power supply system. Each partition preferentially supplies power to the power loads to be restored with a high net recovery benefit, optimizes the recovery plan for each partition and simulates and deduces. Delete the power loads that can be restored after deduction from and update Traverse the power loads to be restored in this way until all the power loads to be restored are restored or there is no available power in the power supply system and the quasi-power supply system;

[0140] Judge Whether it is an empty set. If it is an empty set, end this process; otherwise, go to step S13.

[0141] Step S13: Judge the set Ω of candidate connection points of the power grid MEPS MEPSJ and the set Ω of MEPS MEPS Whether it is an empty set. If any one of them is an empty set, end this process; otherwise, sort the elements in Ω MEPS in descending order of rated power and go to step S14.

[0142] Step S14: Connect the first element Ω SPEM of Ω MEPS (1) to a certain candidate connection point in Ω MEPSJ and regard it as a power supply system to update the PDS partition. The available power of the partition where this Ω MEPS (1) is located is the rated power of this Ω MEPS (1);

[0143] Calculate the net load recovery benefit and the shortest transportation path that can be obtained after the access of Ω MEPS (1) according to formula (1), traverse all candidate connection points in Ω MEPSJ to determine the maximum net load recovery benefit among them;

[0144] If the maximum net load recovery benefit is less than or equal to 0, delete Ω MEPS (1) from Ω MEPS and go to step S13; otherwise, determine the candidate connection point with the maximum net load recovery benefit obtained as the optimal access point of Ω MEPS (1) and go to step S15;

[0145] I M = β M I M,J - C M,J(1)

[0146] Wherein, I M is the net load restoration benefit obtained after the mobile emergency resource accesses the candidate connection point; β M (0 < β M ≤ 1) is the conversion coefficient of the restoration benefit of the mobile emergency resource candidate connection point; I M,J and C M,J are respectively the restoration benefit and cost of the mobile emergency resource for supplying energy to the load to be restored at the candidate connection point.

[0147] Step S15: Connect Ω MEPS (1) to the determined optimal access point to obtain the connection point, the shortest transportation path, and the transportation cost of the mobile emergency power supply MEPS;

[0148] At this time, regard Ω MEPS (1) as a power supply system, update the PDS partition, delete Ω MEPS (1) from Ω MEPS and go to step S12.

[0149] Those skilled in the art should know that when the black start power source point has been determined and there are multiple black start power source points in the system, in order to improve the grid restoration speed, according to the grid structure characteristics, the large-scale grid should be divided into two or more subsystems for independent restoration. After each subsystem operates stably for a period of time, they are interconnected to restore the entire system. Adopting a parallel restoration strategy can greatly shorten the grid restoration time. If any subsystem fails to restore due to some unforeseen factors, it will not affect the restoration process of other subsystems, which improves the system restoration speed. Therefore, the optimal scheduling strategy of the mobile emergency power supply MEPS in this embodiment is mainly used to partition the entire power system, that is, divide the entire large system into several subsystems and implement parallel restoration. The gas system can be solved by the same idea.

[0150] In the embodiment of the present invention, for the optimal scheduling strategy of the mobile emergency gas source MEGS, see Figure 1 , including:

[0151] Considering that each candidate connection point in the GDS can access at most one MEGS, the following steps are used for the optimal scheduling of MEGS:

[0152] Step S21: For each element in Ω GT , determine its initial partition according to the connection situation of each GT in the PDS, determine the power supply range according to the rated power of the GT, and The medium nodes are uniquely partitioned into the GT partitions that can supply power to them and are the closest in distance. The maximum benefit that can be obtained from the GT rated power within the partition is used as their initial benefit in the PDS, and the gas flow rate required by the GT is used as its gas load demand in the GDS. The GTs with zero initial benefit are deleted from Ω GT and updated

[0153] Those skilled in the art should know that the GT is a power supply system in the electrical system and belongs to the gas load to be restored in the gas system.

[0154] In this embodiment, one GT is regarded as a power supply system, each GT is divided into an initial partition, and the rated power of the GT is the available power of the partition.

[0155] Step S22: Set the initial partition of the GDS as the gas supply system with available flow and the quasi-gas supply system itself. Considering the problem of node gas pressure over-limit based on the gas supply range of the GDS initial partition, the medium nodes are uniquely partitioned into the partitions that can supply gas to them and are the closest in distance.

[0156] In this embodiment, each gas supply system or quasi-gas supply system is divided into an initial partition, the available power of the partition is the available power of its gas supply system or quasi-gas supply system, and the gas load nodes to be restored are uniquely partitioned into the partitions that can supply gas to them and are the closest in distance.

[0157] Step S23: Determine whether the set Ω MEGSJ and the set Ω MEGS of the MEGS are empty sets. If any of them is an empty set, end this process; otherwise, sort the elements in Ω MEGS in descending order of the rated flow rate, and go to Step S24.

[0158] Step S24: Connect the first element Ω MEGS of Ω MEGS (1) to a certain candidate connection point in Ω MEGSJ , regard it as a gas supply system, update the GDS partition, and calculate the net benefit of load restoration and the shortest transportation path that can be obtained after connecting according to Equation (1). Traverse all candidate connection points in Ω MEGS (1), determine the maximum net benefit of load restoration among them, MEGSJ If the maximum net benefit of load restoration is less than or equal to 0, delete Ω

[0159] (1) from Ω MEGS and go to Step S23; otherwise, determine the candidate connection point that can obtain the maximum net benefit of load restoration as the optimal access point of Ω MEGS (1), and go to Step S25. MEGS (1), and go to Step S25.

[0160] Step S25: Consider MEGS (1) Connect to the determined optimal access point to obtain the connection point, the shortest transportation path, and the transportation cost of the mobile emergency gas source MEGS;

[0161] At this time, it is equivalent to adding a gas supply system that can provide available flow to the system. Update the GDS partition with this, and MEGS (1) Delete it from MEGS and delete the corresponding candidate connection points from MEGSJ ;

[0162] Step S26: Based on the updated partition, re-optimize and simulate the restoration of the gas loads to be restored in each partition. After that, delete the gas loads that can be restored through simulation from , update

[0163] Judge whether it is an empty set. If it is an empty set, end this process; otherwise, go to Step S23.

[0164] As a preferred implementation, after obtaining the optimal scheduling decision of the mobile emergency resources, for the GDS partition containing the electrical coupling element GT, considering to obtain the maximum net benefit, it is necessary to compare the restoration benefits of restoring the electrical loads to be restored in the PDS partition with the gas flow and the restoration benefits of restoring the gas loads to be restored in the corresponding GDS partition, so as to decide whether GDS support is required during the PDS power supply restoration process.

[0165] As Figure 2 shown, are respectively the maximum net benefit curves that the gas supply volume F GT can obtain in the PDS and GDS. Under the condition that the total available gas flow of the GDS gas source is determined, reduce the gas supply volume of the gas source in the GDS from F 1 by ΔF GT to F 1 ', and increase the gas supply volume obtained in the PDS from F 2 by ΔF GT to F 2 '. In this way, the benefit of the GDS will be reduced from I 1 to I ', while the benefit of the PDS will increase from I 1 to I 2 by to I 2 '. As long as , a higher benefit can be obtained for the total available gas flow. Considering the efficiency and loss of gas-electricity conversion, the embodiment of the present invention only further considers the G2P flow optimization decision when the load to be restored in the PDS is greater than its available power.

[0166] The G2P flow optimization strategy is as follows:

[0167] Establish the objective function for G2P flow optimization as follows:

[0168]

[0169] Wherein, is the gas supply volume F of the gas source GT The difference between the maximum net income obtained in PDS and the maximum net income obtained in GDS. When GDS does not consider supplying gas to PDS for power generation.

[0170] The optimization process is as follows:

[0171] Step S31: Update Ω in descending order of the initial income GT , and judge whether Ω GT is an empty set. If Ω GT is an empty set, end this process; otherwise, go to step S32.

[0172] It should be noted that based on the initial partition in step S21, the maximum income that can be obtained by the GT rated power within the partition is used as its initial income in PDS.

[0173] Step S32: Judge whether there are still nodes to be restored in PDS; if Ω LP is an empty set, end this process; otherwise, go to step S33.

[0174] Step S33: If there is a situation in Ω GT where GT cannot obtain stable gas supply, that is, the available flow in the gas supply system or the quasi-gas supply system partition is less than the lower limit of the GT rated operating flow within the partition, then delete this GT from Ω GT , and judge whether Ω GT is an empty set. If Ω GT is an empty set, end this process; otherwise, go to step S34.

[0175] Step S34: Obtain the first element Ω GT of Ω GT (1) The path with the shortest gas acquisition delay from the gas supply system or the quasi-gas supply system The maximum available flow Wherein, is the available flow of the gas supply system or the quasi-gas supply system, is the rated flow of Ω GT (1), and go to step S35.

[0176] Step S35: Let Ω GT(1) The maximum net benefits that can be obtained for the gas acquisition flow rate in its PDS partition and GDS partition are and Optimize to obtain Ω that satisfies equation (2) GT The optimal gas acquisition flow rate of (1)

[0177] If Then remove Ω GT (1) from Ω GT and go to step S31;

[0178] Otherwise, this gas supply system or quasi - gas supply system plans to supply gas to Ω GT (1) Take Ω GT (1) as a power supply system, update the partition of PDS, optimize the restoration plan and simulate, update and the partitions and initial benefits of other GTs to be restored in PDS, update the remaining available gas flow rate in GDS, update the GDS partition and Take Ω GT (1) from Ω GT and go to step S31.

[0179] As a preferred implementation manner, to improve the energy utilization efficiency and the net restoration benefits of the electrical integrated energy system with mobile emergency resources, based on obtaining the states of PDS, GDS and their coupling components, perform parallel partition of PDS and GDS;

[0180] The process of parallel partition of PDS and GDS is as follows:

[0181] If there is no gas acquisition flow rate for all gas turbines, in PDS, determine the partitions of the power supply system and the quasi - power supply system except gas turbines according to the result of step S11 in the MEPS optimal scheduling decision; in GDS, divide the partitions for the gas supply system and the quasi - gas supply system according to step S22 in the MEGS optimal scheduling decision.

[0182] If there are gas turbines with gas acquisition flow rate, in PDS, according to the method in step S11 of the MEGS optimal scheduling decision, uniquely divide the remaining unpartitioned electrical load nodes into the partition of the gas turbine that can supply power to it, is the closest and has gas acquisition flow rate; in GDS, for the gas supply system and the quasi - gas supply system that plan to supply gas to the gas turbine, divide the nodes related to the gas acquisition path of the gas turbine into its partition, and uniquely divide the remaining unpartitioned nodes into the partition of the gas supply system or the quasi - gas supply system that can supply gas to it and has the shortest gas supply delay.

[0183] Given the transportation paths, transportation costs of the above-mentioned MEPS and MEGS, as well as the energy supply ranges, the order of nodes to be restored, and the restoration amounts of each partition of PDS and GDS, the optimal restoration plan with the maximum net restoration benefit per unit capacity can be determined through the following PDS and GDS parallel partition internal load restoration optimization decision-making.

[0184] The load restoration strategies within the PDS and GDS parallel partitions are shown in Figure 1 , including:

[0185] Step S41: According to the power outage scenario, screen the partition targets to be restored, and generate a sequence Ω of the partition targets to be restored in descending order of the net restoration benefit of the targets to be restored based on Equation (1). TR ;

[0186] Step S42: Select the first target to be restored Ω TR from Ω TR (1), optimize the restoration plan of Ω TR (1), evaluate the restoration benefit and restoration cost of the restoration plan of Ω TR (1) at this step, obtain the corresponding restoration benefit and restoration cost, and simulate and deduce its power outage scenario;

[0187]

[0188]

[0189]

[0190]

[0191] C T = ∑χ T (t)u T (t) (7)

[0192] Among them, and are the electrical load and the gas load restoration benefits at the k P th step respectively; G and are the costs of the power source and the gas source supplying energy to the load at the k P th step respectively; C G T is the transportation cost of MEPS and MEGS; t e is the preset evaluation end time of MG2PIES, set as the larger value after the estimated execution completion time of the PDS and GDS restoration plans; P is the electrical load restoration time at the k P th step; are respectively the unit recovery benefit and recovery volume of the electrical load restored at the k-th P step at time t; are respectively the unit electricity cost and power supply volume for the electrical load provided by the fixed power source at the k-th P step at time t; are respectively the unit electricity cost and power supply volume for the electrical load provided by the MEPS at the k-th P step at time t; is the gas load recovery time at the k-th G step; are respectively the unit recovery benefit and recovery volume of the gas load restored at the k-th G step at time t; are respectively the unit gas cost and gas supply volume for the gas load provided by the fixed gas source at the k-th G step at time t; are respectively the unit gas cost and gas supply volume for the gas load provided by the MEGS at the k-th G step at time t; u T (t) is the state of a certain section of the road passed by the mobile emergency resource at time t during the transportation process, 1 indicates passing, and 0 indicates not passing; χ T (t) is the transportation cost of the mobile emergency resource passing through this section of the road at time t;

[0193] It should be noted that when the unit recovery benefit, unit energy supply cost of electricity and gas loads, and transportation costs of MEPS and MEGS are known, the recovery plan such as the recovery state of each step of the load, the energy supply path state, the energy supply volume, the access state of the mobile emergency resource, and the transportation path can be optimized;

[0194] Step S43: Delete Ω TR (1) from Ω TR , update Ω TR . Since the scenario tree is deduced branch by branch in a breadth-first manner, if the breadth deduction is completed in the current scenario or Ω TR is an empty set, then go to step S44; otherwise, go to step S42 to generate a recovery plan for the next target to be restored in the current scenario;

[0195] Step S44: If the deduction depth of the scenario tree has not reached the set value, then switch to the simulation and power outage scenario of the next step, and go to step S41; otherwise, go to step S45;

[0196] Step S45: Calculate the unit capacity recovery net benefit of all multi-step recovery plans in the current partition, and determine the recovery plan with the largest unit capacity recovery net benefit as the optimal recovery plan. Among them, the ratio of the recovery net benefit of the recovery plan to its recovery capacity is defined as its unit capacity recovery net benefit.

[0197] In the embodiment of the present invention, an objective function is established with the maximum net income during the restoration process of the integrated electrical energy system as the goal, as follows:

[0198]

[0199] Among them, and are respectively the total number of steps of the restoration schemes on the power distribution side PDS and the gas distribution side GDS within the same decision period; k P and k G are respectively the step numbers in the PDS and GDS restoration schemes. One-step scheme refers to a series of operations taken to restore a certain object.

[0200] In the embodiment of the present invention, the constraint conditions that the objective function needs to satisfy are as follows:

[0201] A. PDS constraints, including steady-state constraints and topological structure constraints

[0202] A1. Steady-state constraints:

[0203]

[0204]

[0205]

[0206]

[0207]

[0208] U i,min ≤U i (t)≤U i,max ,i∈Ω B (14)

[0209] I ij,min ≤I ij (t)≤I ij,max ,i,j∈Ω B (15)

[0210]

[0211] Among them, Ω B is the set of energized nodes in the power supply area; Ω MEPS,i is the set of MEPSs connected to grid node i; P i grid (t), are respectively the active and reactive powers injected by the external network at node i at time t; P iDG (t) and are the active and reactive powers injected by the distributed power source at node i at time t; P i PM (t) and are the active and reactive powers injected by the MEPS at node i at time t; P i L (t) and are the active and reactive powers of the load at node i at time t; U i (t), U j (t) are the voltage values of nodes i and j at time t; G ij and B ij are the real and imaginary parts of the admittance of line ij; θ ij (t) is the phase difference between the voltages of nodes i and j at time t; are the upper and lower limits of the active power of the distributed power source at node i; are the rated power and the active power output at time t of the mobile emergency power source c, c ∈ Ω MEPS,i ; is the energy supply time of the mobile emergency power source c; U i,max and U i,min are the upper and lower limits of the voltage of node i; I ij (t), I ij,max , I ij,min are the current of branch ij and its upper and lower limits at time t; are the transmission capacity of branch ij and its upper limit at time t.

[0212] A2. Transient constraint:

[0213] To avoid misoperation of protection caused by transient impact caused by the connection of electrical loads, a maximum single-step restored electrical load constraint is set,

[0214] ΔP k,step ≤ ζP kN (17)

[0215] where ΔP k,step is the active power capacity of the single-step restored electrical load at the k-th step; ζ is the proportionality coefficient of the allowable single-step restored electrical load capacity; P kN is the rated capacity within the energized area at the initial moment of the k-th step.

[0216] A3. Topological structure constraint:

[0217] g ∈ G (18)

[0218] where g is the topological structure of the energized area after PDS reconstruction; G is the set of all radial network topologies of PDS.

[0219] B, GDS constraints, including steady-state constraints and transient constraints

[0220] B1. Steady-state constraints:

[0221]

[0222]

[0223]

[0224]

[0225]

[0226] Π m,min ≤ Π m (t) ≤ Π m,max , m ∈ Ω G (24)

[0227] F mn,min ≤ F mn (t) ≤ F mn,max , m, n ∈ Ω G (25)

[0228] where Ω G is the set of gas acquisition nodes in the gas supply area; are the gas flow rates injected by the external network and the gas source at node m at time t, respectively; is the gas flow rate injected by the MEGS at node m at time t; is the gas flow rate of the load at node m at time t; F mn (t) is the pipe segment flow rate from node m to node n at time t; are the upper and lower limits of the gas injection flow rate at node m, respectively; the mobile emergency gas source d is the MEGS connected to node m of the gas network; are the rated flow rate and the output flow rate at time t of the mobile emergency gas source d, respectively; is the energy supply time of the MEGS at node m; Π m (t), Π m,max , Π m,min are the pressure value at node m at time t and its upper and lower limits, respectively; F mn,max , F mn,min are the upper and lower limits of the pipe segment mn flow rate, respectively.

[0229] B2. Transient constraints:

[0230] The pressure and flow rate of the natural gas injected by the gas turbine should be within the regulation range of the regulating valve of the gas regulation system,

[0231]

[0232]

[0233] Among them, Ω GT is the set of gas turbines GT to be restored; are the natural gas pressure value injected into gas turbine a at time t and its upper and lower limits, respectively; are the natural gas flow rate value injected into gas turbine a at time t and its upper and lower limits, respectively.

[0234] C. Transportation constraints of mobile emergency resources, including:

[0235] C1. Constraint of mobile emergency power supply:

[0236]

[0237] Among them, Ω MEPSJ is the set of candidate connection points of the power grid MEPS; is the number of MEPS units connected to node i; N MEPS is the total number of MEPS units.

[0238] C2. Constraint of mobile emergency gas source:

[0239]

[0240] Among them, Ω MEGSJ is the set of candidate connection points of the gas network MEGS; is the number of MEGS units connected to node m; N MEGS is the total number of MEGS units.

[0241] C3. Transportation path constraint:

[0242] l ∈ L (30)

[0243] Among them, l is the transportation path of mobile emergency resources from their storage location to the candidate connection point; L is the set of all radial paths of the transportation network. In this embodiment, the transportation path of each mobile emergency resource is constrained to be open-loop.

[0244] C4. Transportation time constraint:

[0245]

[0246] Among them, t 0 is the starting time of restoration; is the transportation time of mobile emergency power supply c passing through section pq; l c is the shortest transportation path of mobile emergency power supply c from its storage location to the candidate connection point; is the startup time consumption of mobile emergency power supply c.

[0247] The transportation time constraint of MEGS is the same as that of MEPS.

[0248] The present invention also provides an energy supply restoration coordination system for an electrical integrated energy system with mobile emergency resources, which is used to implement the energy supply restoration coordination method for the electrical integrated energy system with mobile emergency resources as described above. The system includes:

[0249] An initialization module, which is used to identify the state of the electrical integrated energy system with mobile emergency resources, analyze the states of the system's load to be restored, energy supply resources, branches, and switch equipment, obtain the set of electrical loads to be restored, the gas load to be restored, and the set of gas turbines to be restored, and determine the energy supply scope of each energy supply system; the mobile emergency resources include mobile emergency power sources (MEPS) and mobile emergency gas sources (MEGS);

[0250] A mobile emergency power source coordination module, which is used to optimize the configuration of the partition strategy for the electrical load to be restored on the distribution side. Specifically: divide partitions for the fixed power sources, and restore the electrical load to be restored within each partition based on the available power of the fixed power sources until the fixed power sources have no available power; with the goal of maximizing the net benefit of electrical load restoration, connect the mobile emergency power source (MEPS) to each candidate connection point in the power grid in sequence and update the partition until all the electrical load to be restored is restored or the MEPS has no available power, and obtain the optimal partition strategy on the distribution side, the transportation paths, transportation costs, and corresponding connection points of each mobile emergency power source (MEPS);

[0251] A mobile emergency gas source coordination module, which is used to optimize the configuration of the partition strategy for the gas load to be restored on the gas distribution side. Specifically: calculate the initial benefit and required flow rate of the gas turbine, divide partitions for the fixed gas sources, and with the goal of maximizing the net benefit of gas load restoration, connect the mobile emergency gas source (MEGS) to each candidate connection point in the gas network in sequence and update the partition until all the gas load to be restored is restored or the MEGS has no available power, and obtain the optimal partition strategy on the gas distribution side, the transportation paths, transportation costs, and corresponding connection points of each mobile emergency gas source (MEGS);

[0252] An electrical coupling element coordination module, which is used to optimize the G2P flow configuration of the electrical coupling element gas turbine, and with the goal of maximizing the net benefit of load restoration, allocate the flow rate for supplying the gas load to be restored and the flow rate for supplying the electrical load to be restored through the gas turbine, and obtain the switch states, gas supply systems, gas supply volumes, and gas supply paths of each gas turbine;

[0253] A parallel recovery module is used to perform parallel zoning on the power distribution side and the gas distribution side according to the optimal zoning strategy on the power distribution side, the optimal zoning strategy on the gas distribution side, and the flow configuration of the gas turbine, determine the loads to be restored in each zone; and perform load restoration on each zone of the power distribution side and the gas distribution side with the goal of maximizing the net recovery benefit during the restoration process of the integrated electrical energy system, and obtain the optimal restoration plan for each load. The restoration plan includes the restoration status of each step of the load, the status of the energy supply path, the energy supply quantity, the connection points of mobile emergency resources, and the transportation path.

[0254] Those skilled in the art should understand that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can take the form of a computer program product implemented 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.

[0255] The present application is described with reference to the flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or block in the flowchart and / or block diagram can be implemented by computer program instructions, and the combination of the processes and / or blocks in the flowchart and / or block diagram can also be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate means for implementing the specified functions in Figure 1 one process or multiple processes and / or blocks Figure 1 one block or multiple blocks.

[0256] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory generate a manufactured article including instruction means, and the instruction means implements the specified functions in Figure 1 one process or multiple processes and / or blocks Figure 1 one block or multiple blocks.

[0257] These computer program instructions can also be loaded onto a computer or other programmable data processing device, so that a series of operation steps are executed on the computer or other programmable device to generate a computer-implemented process. Thus, the instructions executed on the computer or other programmable device provide steps for implementing the specified functions in Figure 1 one process or multiple processes and / or blocks Figure 1 one block or multiple blocks.

[0258] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit them. Although the present invention has been described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: it is still possible to modify the specific implementation manners of the present invention or make equivalent substitutions, and any modification or equivalent substitution that does not depart from the spirit and scope of the present invention shall be covered by the protection scope of the claims of the present invention.

Claims

1. A method for coordinating the energy supply restoration of an electrical integrated energy system with mobile emergency resources, characterized in that, it includes: Conduct state identification on the electrical integrated energy system with mobile emergency resources, analyze the states of the system's load to be restored, energy supply resources, branches, and switch equipment, obtain the set of electrical loads to be restored, the gas load to be restored, and the set of gas turbines to be restored, and determine the energy supply scope of each energy supply system; the mobile emergency resources include mobile emergency power sources (MEPS) and mobile emergency gas sources (MEGS); the energy supply systems include the power supply system on the distribution side and the gas supply system on the gas distribution side; On the distribution side, optimize the configuration of the partition strategy for the electrical load to be restored. Specifically: divide partitions for the fixed power sources, and within each partition, restore the electrical load to be restored based on the available power of the fixed power sources until there is no available power for the fixed power sources; With the goal of maximizing the net revenue of electrical load restoration, sequentially connect the mobile emergency power source (MEPS) to each candidate connection point in the power grid and update the partitions until all the electrical load to be restored is restored or the MEPS has no available power, obtaining the optimal partition strategy on the distribution side, the transportation paths, transportation costs, and corresponding connection points of each mobile emergency power source (MEPS); On the gas distribution side, optimize the configuration of the partition strategy for the gas load to be restored. Specifically: calculate the initial revenue and required flow rate of the gas turbine, divide partitions for the fixed gas sources, and with the goal of maximizing the net revenue of gas load restoration, sequentially connect the mobile emergency gas source (MEGS) to each candidate connection point in the gas network and update the partitions until all the gas load to be restored is restored or the MEGS has no available power, obtaining the optimal partition strategy on the gas distribution side, the transportation paths, transportation costs, and corresponding connection points of each mobile emergency gas source (MEGS); Conduct G2P flow optimization configuration for the gas turbine, an electrical coupling component. With the goal of maximizing the net revenue of load restoration, allocate the flow rate supplied to the gas load to be restored and the flow rate supplied to the electrical load to be restored through the gas turbine, obtaining the switch states, gas supply systems, gas supply volumes, and gas supply paths of each gas turbine; According to the optimal partition strategy on the distribution side and the optimal partition strategy on the gas distribution side, as well as the flow configuration of the gas turbine, conduct parallel partitioning of the distribution side and the gas distribution side to determine the load to be restored within each partition; Conduct load restoration for each partition on the distribution side and the gas distribution side with the goal of maximizing the net revenue during the restoration process of the electrical integrated energy system, obtaining the optimal restoration plan for each load. The restoration plan includes the restoration states of each step of the load, the states of the energy supply paths, the energy supply volumes, the connection points of the mobile emergency resources, and the transportation paths.

2. The method for coordinating the energy supply restoration of an electrical integrated energy system with mobile emergency resources according to claim 1, characterized in that, the conduct of state identification on the electrical integrated energy system with mobile emergency resources, analyzing the states of the system's load to be restored, energy supply resources, branches, and switch equipment, obtaining the set of electrical loads to be restored, the gas load to be restored, and the set of gas turbines to be restored, and determining the energy supply scope of each energy supply system includes: Identify isolated power outage areas and non-isolated power outage areas according to the topological structure of the electrical integrated energy system and the fault conditions of the loads to be restored, energy supply resources, branches, and switching equipment, and then obtain the set of electrical loads to be restored Gas load to be restored Set and the set Ω of gas turbines to be restored GT ; The isolated power outage area refers to the power outage area that has no recoverable path to all energy supply systems or quasi-energy supply systems; The quasi-energy supply system refers to the system of sources that are in a power outage state but have self-starting ability; The quasi-energy supply system includes the quasi-power supply system on the distribution side and the quasi-gas supply system on the gas distribution side; Determine the power supply scope of the power supply system or quasi-power supply system according to the available power of the power supply system or quasi-power supply system and the voltage level on the power transmission path; Determine the gas supply range of the gas supply system or quasi-gas supply system according to the available flow rate of the gas supply system or quasi-gas supply system and the natural gas pressure level on the gas transmission path.

3. A power supply restoration coordination method for an electrical integrated energy system with mobile emergency resources according to claim 2, characterized in that partitioning the fixed power supplies, and restoring the electrical loads to be restored within each partition based on the available power of the fixed power supplies until there is no available power in the fixed power supplies, including: Step S11: Based on the power supply systems and quasi-power supply systems with available power, initially partition the distribution side, divide each power supply system or quasi-power supply system into an initial partition, and uniquely partition the electrical load nodes to be restored into the partition that can supply power to them and is the closest; Step S12: Each partition preferentially powers the electricity load to be restored with the highest net recovery benefit, and conducts simulation and deduction. The electricity load that can be restored through deduction is deleted from and is updated. Traverse the electricity load to be restored in this way until all the electricity load to be restored is restored or there is no available power in both the power supply system and the quasi-power supply system.

4. A power supply restoration coordination method for an electrical integrated energy system with mobile emergency resources according to claim 3, characterized in that sequentially connecting the mobile emergency power source MEPS to each candidate connection point in the power grid and updating the partition with the goal of maximizing the net revenue of electrical load restoration until all the electrical loads to be restored are restored or the MEPS has no available power, including: Step S13: If is not an empty set, then determine whether the set Ω of candidate connection points of the power grid MEPS MEPSJ and the set Ω of MEPS MEPS are empty sets. If any one of them is an empty set, end; otherwise, sort the elements in Ω MEPS in descending order of rated power and go to Step S14; Step S14: Take Ω MEPS The first element of Ω MEPS (1) Connect to a certain candidate connection point of Ω MEPSJ and use Ω MEPS (1) to update the distribution side partition as a power supply system; Ω is calculated according to the following formula MEPS (1) After connection, the net revenue of load restoration and the shortest transportation path are traversed for Ω MEPSJ All candidate connection points in are determined to obtain the maximum net revenue of load restoration; I M = β M I M,J - C M,J ; Among them, I M is the net load recovery benefit obtained after the mobile emergency resource access candidate connection point; β M is the conversion coefficient of the recovery benefit of the mobile emergency resource candidate connection point, I M,J , C M,J are respectively the recovery benefit and cost of the mobile emergency resource for supplying energy to the load to be recovered at the candidate connection point; If the net revenue of maximum load restoration is less than or equal to 0, then Ω MEPS (1) Delete it from Ω MEPS and go to step S13; otherwise, determine the candidate connection point with the obtained net revenue of maximum load restoration as the optimal access point of MEPS (1), and go to step S15; Step S15: Take Ω MEPS (1) Connect to the determined optimal access point to obtain the connection point, the shortest transportation path, and the transportation cost of the mobile emergency power supply (MEPS); Take Ω MEPS (1) As a power supply system, update the distribution side partition and take Ω MEPS (1) Delete from Ω MEPS and go to step S12.

5. A power supply restoration coordination method for an electrical integrated energy system with mobile emergency resources according to claim 2, characterized in that calculating the initial revenue and required flow rate of the gas turbine, partitioning the fixed gas sources, and sequentially connecting the mobile emergency gas source MEGS to each candidate connection point in the gas network and updating the partition with the goal of maximizing the net revenue of gas load restoration until all the gas loads to be restored are restored or the MEGS has no available power, including: Step S21: For the elements in Ω GT , take a gas turbine as a power supply system, and uniquely partition the nodes in into the gas turbine partition that can supply power to them and is the closest in distance. Take the maximum benefit that can be obtained from the rated power of the gas turbine within the partition as the initial benefit on the distribution side, and take the gas flow required by the gas turbine as its gas load demand on the gas distribution side; Delete the GT with an initial benefit equal to zero from Ω GT , and update Step S22: Divide each gas supply system or quasi-gas supply system into an initial partition, and uniquely partition the nodes in into the partition that can supply gas to them and is the closest in distance; Step S23: Determine the candidate connection point set Ω of the gas network MEGS MEGSJ and the set Ω of MEGS MEGS whether it is an empty set. If any one of them is an empty set, end; otherwise, sort the elements in Ω MEGS in descending order of the rated flow rate and go to Step S24; Step S24: Take Ω MEGS The first element of Ω MEGS (1) Connect to a certain candidate connection point of Ω MEGSJ Take Ω MEGS (1) As a gas supply system, update the gas distribution side partition, calculate the net benefit of load restoration and the shortest transportation path, traverse all candidate connection points in Ω MEGSJ to determine the maximum net benefit of load restoration; If the net revenue of maximum load restoration is less than or equal to 0, then Ω MEGS (1) Delete from Ω MEGS and go to step S23; otherwise, determine the candidate connection point that obtains the net revenue of maximum load restoration as the optimal access point of Ω MEGS (1), and go to step S25; Step S25: Take Ω MEGS (1) Connect to the determined optimal access point to obtain the connection point, shortest transportation path and transportation cost of the mobile emergency gas source MEGS; Take Ω MEGS (1) As a gas supply system, update the partition on the gas distribution side and take Ω MEGS (1) Delete from Ω MEGS and delete the determined optimal access point from Ω MEGSJ ; Step S26: Based on the updated partitions, after re-optimizing and simulating the gas loads to be restored in each partition, delete the gas loads that can be restored after the simulation deduction from , and update Judge whether it is an empty set. If it is an empty set, end; otherwise, go to step S23.

6. A power supply restoration coordination method for an electrical integrated energy system with mobile emergency resources according to claim 5, characterized in that performing G2P flow optimization configuration on the electrical coupling component gas turbine, including: Step S31: Take the maximum benefit that can be obtained from the rated power of the gas turbine within the partition as its initial benefit on the distribution side, and update Ω by arranging it in descending order of the initial benefit GT , and determine Ω GT Whether it is an empty set. If Ω GT is an empty set, then end; otherwise, go to Step S32; Step S32: Determine whether there is still a load to be restored on the distribution side; if it is an empty set, then end; otherwise, go to step S33; Step S33: If there is a gas turbine in Ω GT that cannot obtain stable gas supply, then remove it from Ω GT , and determine whether Ω GT is an empty set. If Ω GT is an empty set, then end; otherwise, go to Step S34; Step S34: Obtain Ω GT The first element Ω GT (1) The path with the shortest gas acquisition delay from the gas supply system or the quasi-gas supply system The maximum available gas flow rate wherein is the available gas flow rate of the gas supply system or the quasi-gas supply system, is Ω GT (1) The rated gas flow rate, go to step S35; Step S35: Set Ω GT (1) The maximum net benefits obtained by the gas acquisition flow rate in its power distribution side partition and gas distribution side partition are respectively and Optimize to make The maximum Ω GT (1) The optimal gas acquisition flow rate If then delete Ω GT (1) from Ω GT and go to step S31; Otherwise, allocate the gas supply system or quasi-gas supply system to plan for Ω GT (1) Gas supply Take Ω GT (1) As a power supply system, update the distribution side zoning, optimize the load restoration plan and simulate and deduce, update The zoning and initial benefits of Ω and other gas turbines on the distribution side, update the remaining available gas flow on the gas distribution side, update the gas distribution side zoning and Take Ω GT (1) Delete from Ω GT and go to step S31.

7. A power supply restoration coordination method for an electrical integrated energy system with mobile emergency resources according to claim 6, characterized in that performing parallel partitioning on the distribution side and the gas distribution side, including: If there is no gas acquisition flow rate for all gas turbines, on the distribution side, determine the partitions of the power supply systems and quasi-power supply systems other than the gas turbines according to the optimal partitioning strategy of the distribution side; on the gas distribution side, divide the partitions of the gas supply systems and quasi-gas supply systems according to the optimal partitioning strategy of the gas distribution side; If there are gas turbines with gas acquisition flow rates, on the distribution side, uniquely partition the remaining electrical load nodes to be partitioned into the partition of the gas turbine that can supply power to them, is the closest, and has a gas acquisition flow rate according to the optimal partitioning strategy of the distribution side; on the gas distribution side, for the gas supply systems and quasi-gas supply systems planned to supply gas to the gas turbines, partition the nodes related to the gas acquisition path of the gas turbines into their partitions, and uniquely partition the remaining nodes to be partitioned into the partition of the gas supply system or quasi-gas supply system that can supply gas to them and has the shortest gas transmission delay.

8. A power supply restoration coordination method for an electrical integrated energy system with mobile emergency resources according to claim 7, characterized in that performing load restoration with the goal of maximizing the net revenue during the restoration process of the electrical integrated energy system within each partition of the distribution side and the gas distribution side to obtain the optimal restoration plan for each load. Step S41: In each partition of the power distribution side and the gas distribution side, according to the power outage scenario, the net recovery benefits of the targets to be restored are sorted in descending order to generate a sequence Ω of the targets to be restored in the partition. TR The targets to be restored are the electrical loads to be restored and the gas loads to be restored. Step S42: Select the first target to be restored Ω TR from Ω TR (1), and evaluate the restoration benefit and restoration cost of restoring Ω TR (1) based on the following formula, and simulate and deduce its power outage scenario: C T = ∑χ T (t)u T (t); Among them, and are the electrical load at the k P -th step and the gas load recovery benefit at the k G -th step respectively; and are the power supply at the k P -th step and the cost of the gas source supplying energy to the load at the k G -th step respectively; C T is the transportation cost of MEPS and MEGS; t e is the preset evaluation end time of the system, which is set to the larger value after the estimated execution completion time of the power distribution side and gas distribution side recovery plans; is the electrical load recovery time at the k P -th step; are the unit recovery benefit and recovery volume of the electrical load restored at the k P -th step at time t respectively; are the unit power cost and power supply of the fixed power supply providing power to the electrical load at time t at the k P -th step respectively; are the unit power cost and power supply of MEPS providing power to the electrical load at time t at the k P -th step respectively; is the gas load recovery time at the k G -th step; are the unit recovery benefit and recovery volume of the gas load restored at the k G -th step at time t respectively; are the unit gas cost and gas supply of the fixed gas source providing flow to the gas load at time t at the k G -th step respectively; are the unit gas cost and gas supply of MEGS providing flow to the gas load at time t at the k G -th step respectively; u T (t) is the state of a certain section of the road passed by the mobile emergency resource at time t during the transportation process, 1 means passed, 0 means not passed; χ T (t) is the transportation cost of the mobile emergency resource passing through this section of the road at time t; Step S43: Take Ω TR (1) Delete it from Ω TR and update Ω TR . If the breadth deduction is completed in the current scenario or Ω TR is an empty set, then go to Step S44; otherwise, go to Step S42 to generate a recovery plan for the next target to be recovered in the current scenario. Step S44: if the deduction depth does not reach the set value, switch to the next step of simulating the power outage scenario and go to step S41; otherwise, go to step S45; Step S45: Calculate the unit capacity recovery net benefits of all multi-step recovery schemes in the current partition, and determine the optimal recovery scheme with the goal of maximizing the unit capacity recovery net benefits of all multi-step recovery schemes in the partition during the recovery process of the electric integrated energy system; wherein the ratio of the recovery net benefit of the recovery scheme to its recovery capacity is defined as the unit capacity recovery net benefit; The objective function of maximizing the net benefits of all multi-step restoration schemes in the partition of the electrical integrated energy system restoration process is as follows: Among them, and are the total number of steps of the power distribution side and the gas distribution side restoration plans within the same decision-making period, respectively; k P and k G are the step numbers in the PDS and GDS restoration plans, respectively.

9. A method for coordinating energy supply restoration of an electrical integrated energy system including mobile emergency resources according to claim 8, It is characterized in that The objective function must satisfy the following constraints: A. Distribution side constraints, including steady-state constraints and topology constraints: A1. Steady-state constraints: U i,min ≤U i (t)≤U i,max , i ∈ Ω B ; I ij,min ≤I ij (t)≤I ij,max , i, j ∈ Ω B ; Among them, Ω B is the set of energized nodes in the power supply area; Ω MEPS,i is the set of MEPSs connected to grid node i; P i grid (t), are the active and reactive powers injected by the external grid at node i at time t, respectively; P i DG (t), Q i DG (t) are the active and reactive powers injected by the distributed power source at node i at time t, respectively; P i PM (t), are the active and reactive powers injected by the MEPS at node i at time t, respectively; P i L (t), are the active and reactive powers of the load at node i at time t, respectively; U i (t), U j (t) are the voltage values of nodes i and j at time t, respectively; G ij , B ij are the real and imaginary parts of the admittance of line ij, respectively; θ ij (t) is the phase difference between the voltages of nodes i and j at time t; are the upper and lower limits of the active power of the distributed power source at node i, respectively; are the rated power and the active power output at time t of the mobile emergency power source c, respectively, where c ∈ Ω MEPS,i ; is the energy supply time of the mobile emergency power source c; U i,max , U i,min are the upper and lower limits of the voltage of node i, respectively; I ij (t), I ij,max , I ij,min are the current of branch ij and its upper and lower limits at time t, respectively; are the transmission capacity of branch ij and its upper limit at time t, respectively; A2. Transient constraints: ΔP k,step ≤ζP kN ; Among them, ΔP k,step is the active power capacity for a single-step restoration of the electrical load; ζ is the proportionality coefficient for the allowable single-step restoration capacity of the electrical load; P kN is the rated capacity within the energized area at the initial moment of the k-th step; A3. Topological constraints: g∈G; Among them, g is the topological structure of the live domain after the distribution side is reconstructed; G is the set of all radial network topologies on the distribution side; B. Valve side constraints, including steady-state constraints and transient constraints: B1. Steady-state constraints: Π m,min ≤Π m (t)≤Π m,max ,m∈Ω G ; F mn,min ≤F mn (t)≤F mn,max , m, n ∈ Ω G ; Among them, Ω G is the set of gas acquisition nodes in the gas supply area; are the gas flow rates injected by the external network and the gas source at node m at time t, respectively; is the gas flow rate injected by MEGS at node m at time t; is the gas flow rate of the load at node m at time t; F mn (t) is the pipe segment flow rate from node m to node n at time t; are the upper and lower limits of the gas injection flow rate at the gas source at node m, respectively; The mobile emergency gas source d is the MEGS connected to node m of the gas network; are the rated flow rate of the mobile emergency gas source d and the output flow rate at time t, respectively; is the energy supply time of MEGS at node m; Π m (t), Π m,max , Π m,min are the pressure value of node m at time t and its upper and lower limits, respectively; F mn,max , F mn,min are the upper and lower limits of the flow rate of pipe segment mn, respectively; B2. Transient constraints: Among them, Ω GT is the set of gas turbines to be restored; are the natural gas pressure value injected into gas turbine a at time t and its upper and lower limits respectively; are the natural gas flow rate value injected into gas turbine a at time t and its upper and lower limits respectively; C. Mobile emergency resource transportation constraints, including: C1. Mobile emergency power supply restrictions: Among them, Ω MEPSJ is the set of candidate connection points of the power grid MEPS; is the number of units of node i accessing the MEPS; N MEPS is the total number of MEPS units; C2. Mobile emergency gas source constraints: Among them, Ω MEGSJ is the set of candidate connection points for the gas network MEGS; is the number of MEGS units connected to node m; N MEGS is the total number of MEGS units; C3. Transportation path constraints: l∈L； Among them, l is the transportation path of the mobile emergency resource from its storage location to the candidate connection point; L is the set of all radial paths in the transportation network; C4. Transportation time constraints: Among them, t 0 is the starting moment of restoration; is the transportation time for the mobile emergency power supply c to pass through the section pq; l c is the shortest transportation path for the mobile emergency power supply c from its storage location to the candidate connection point; is the startup time consumption of the mobile emergency power supply c; The transportation time constraints for MEGS are the same as those for MEPS.

10. An electrical integrated energy system energy supply recovery coordination system including mobile emergency resources, It is characterized in that A method for coordinating energy supply restoration of an electrical integrated energy system containing mobile emergency resources according to any one of claims 1 to 9, the system comprising: An initialization module is used to identify the state of the electrical integrated energy system containing mobile emergency resources, analyze the state of the system's loads to be restored, energy supply resources, branches and switchgear, obtain a set of electric loads to be restored, a set of gas loads to be restored and a set of gas turbines to be restored, and determine the energy supply range of each energy supply system; the mobile emergency resources include a mobile emergency power supply MEPS and a mobile emergency gas source MEGS; the energy supply system includes a power supply system on the power distribution side and a gas supply system on the gas distribution side; The mobile emergency power coordination module is used to optimize the configuration of the partition strategy for the electric load to be restored on the distribution side, specifically: divide the fixed power supply into zones, and restore the electric load to be restored based on the available power of the fixed power supply in each zone until the fixed power supply has no available power; with the goal of maximizing the net benefit of electric load restoration, the mobile emergency power supply MEPS is connected to each candidate connection point in the power grid in turn and the partition is updated until all the electric load to be restored is restored or the MEPS has no available power, and the optimal partition strategy on the distribution side, the transportation path, transportation cost and corresponding connection point of each mobile emergency power supply MEPS are obtained; The mobile emergency gas source coordination module is used to optimize the partition strategy configuration for the gas load to be restored on the gas distribution side. Specifically, it calculates the initial revenue and required flow rate of the gas turbine, divides partitions for the fixed gas sources, and sequentially connects the mobile emergency gas sources (MEGS) to each candidate connection point in the gas network with the goal of maximizing the net revenue of gas load restoration and updates the partitions until all the gas loads to be restored are restored or the MEGS has no available power, so as to obtain the optimal partition strategy on the gas distribution side, the transportation paths, transportation costs and corresponding connection points of each MEGS. The electrical coupling component coordination module is used to optimize the G2P flow rate configuration of the electrical coupling component gas turbine. With the goal of maximizing the net revenue of load restoration, it allocates the flow rate supplied to the gas load to be restored and the flow rate supplied to the electrical load to be restored through the gas turbine, so as to obtain the switching state, gas supply system, gas supply volume and gas supply path of each gas turbine. The parallel restoration module is used to perform parallel partitioning on the power distribution side and the gas distribution side according to the optimal partition strategy on the power distribution side and the optimal partition strategy on the gas distribution side, as well as the flow rate configuration of the gas turbine, to determine the loads to be restored in each partition; and perform load restoration for each partition on the power distribution side and the gas distribution side with the goal of maximizing the net revenue during the restoration process of the electrical integrated energy system, so as to obtain the optimal restoration plan for each load. The restoration plan includes the restoration state of each step of the load, the state of the energy supply path, the energy supply volume, the connection points of the mobile emergency resources and the transportation paths.

Images