Energy supply restoration method and system for integrated energy system considering demand response compensation evaluation

Abstract

A method and system for energy supply restoration of an integrated energy system considering demand response compensation evaluation, which evaluate the compensations for the reducible energy consumption loads in each zone in different energy supply restoration schemes; use the revenues, costs, evaluation values of the compensations for the reducible energy consumption loads, transportation costs of mobile energy, and risk costs in each zone in different energy supply restoration schemes to form the net revenues of energy supply restoration in each zone; establish an optimization module for system energy supply restoration. The upper-layer model determines the zone scheme with the goal of maximizing the total sum of the net revenues of energy supply restoration in all zones, and optimizes the zone scheme according to the real-time scheduling of mobile energy; the lower-layer model, based on the optimized zone scheme, re-evaluates the compensations for the optimized reducible energy consumption loads in each zone in different energy supply restoration schemes, and determines the energy supply restoration scheme for each zone with the goal of maximizing the net revenues of energy supply restoration in each updated zone; iteratively obtain the optimal zones of the integrated energy system and the energy supply restoration schemes for each optimal zone, realizing the energy supply restoration considering demand response compensation evaluation.

Description

Technical Field

[0001] The present invention belongs to the technical field of integrated energy demand response. Specifically, it relates to a power supply restoration method and system for an electrical integrated energy system considering the compensation evaluation of integrated load demand response. Background Art

[0002] With the changes in energy demand and energy structure, the integrated energy system (IES) has become the main body to meet user needs and energy supply guarantee. The integrated energy system is usually based on electricity, with multiple energy sources collaborating, covering various energy forms such as electricity, gas, cold, and heat. The coupling between each energy network is tight. However, the integrated energy system also faces problems such as flexibility during operation and stability after being affected by disasters. Demand-side response is an effective way to participate in grid interaction. Demand-side response usually refers to the power department taking various incentive subsidy measures to stimulate the initiative of users to participate. By changing the user's electricity consumption behavior through increasing or decreasing the electricity load, it promotes the dynamic balance between system supply and demand. According to different classification criteria of the current demand response mechanism, the demand response mechanism can be divided into two categories according to the guiding method and the evaluation method of the contribution made by users to the system. Among them, the guiding method can be divided into price-based and incentive-based. The response effect of the price-based type depends on the formulation of electricity prices and has a certain degree of unpredictability. The incentive-based demand response does not involve the formulation of electricity prices. Through incentive compensation, the enthusiasm of a large number of users can be mobilized, and it has the advantages of strong response ability and fast response speed. In summary, demand response, as a form of flexible participation on the load side, can effectively enhance the stability of the system.

[0003] In the prior art, most of the demand responses aim at minimizing the system operation cost, formulating a comprehensive demand response method for a multi-energy system, and improving the reliability of system operation. However, factors such as the cuttable capacity of the contracted demand response electrical load affected by the power outage scenario and the recovery process and the system recovery adequacy are not considered. The comprehensive load demand response compensation evaluation method and its effect in the recovery scenario are not involved. Moreover, currently, demand response usually provides auxiliary services in normal operation modes such as peak shaving or standby for the system. And whether it is time-scheduled or real-time demand response service, a preset fixed unit capacity compensation price is usually adopted to participate in the operation. Summary of the Invention

[0004] To solve the deficiencies in the prior art, the present invention provides a power supply restoration method and system for an electrical integrated energy system considering the compensation evaluation of integrated load demand response, realizing the demand response compensation evaluation during the power outage recovery process, solving the evaluation of aspects such as the cuttable capacity of the demand response electrical load affected by the power outage scenario and the recovery process and the system recovery adequacy. The evaluation result represents the urgency of the system for demand response during the recovery process, which is beneficial to timely taking measures to effectively improve the stability of the system.

[0005] The present invention adopts the following technical solutions.

[0006] The present invention proposes a method for restoring energy supply in an integrated energy system considering demand response compensation evaluation. The reducible energy loads of the system include: reducible electricity loads and reducible gas loads; the method includes:

[0007] Establish a unit capacity compensation cost model for reducible energy loads to evaluate the compensation for reducible energy loads in each partition within the system in different energy supply restoration schemes; use the benefits, costs, evaluation values of reducible energy load compensation, transportation costs of mobile energy, and risk costs in each partition in different energy supply restoration schemes to form the net energy supply restoration benefits of each partition; based on a two-layer model structure, establish an optimization module for system energy supply restoration; among them, the upper-layer model determines the partition scheme with the goal of maximizing the sum of the net energy supply restoration benefits of all partitions, and optimizes the partition scheme according to the real-time scheduling of mobile energy; the lower-layer model, based on the optimized partition scheme, uses the unit capacity compensation cost model of reducible energy loads to re-evaluate the compensation for reducible energy loads in each optimized partition in different energy supply restoration schemes, and determines the energy supply restoration scheme for each partition with the goal of maximizing the updated net energy supply restoration benefits of each partition; the lower-layer model feeds back the updated net energy supply restoration benefits of each partition to the upper-layer model; through multiple iterations, the optimal partition of the integrated energy system and the energy supply restoration schemes for each optimal partition are obtained.

[0008] Preferably, collect data on reducible electricity loads to determine the proportion and reduction rate of reducible electricity loads, and determine the cumulative reduction times and cumulative reduction amount proportions of reducible electricity loads at each node within the system; collect data on distributed power sources to determine the smoothness of the predicted output values of distributed power sources; collect and use data on electricity loads and distributed power sources to determine the adequacy of energy supply restoration in the area to be restored; use the proportion and reduction rate of reducible electricity loads, the cumulative reduction times and cumulative reduction amount proportions of reducible electricity loads at each node within the system, the smoothness of the predicted output values of distributed power sources, and the adequacy of energy supply restoration in the area to be restored to establish a unit capacity compensation cost model for reducible electricity loads;

[0009] Among them, the data on reducible electricity loads includes: reduction amount, rated power, power.

[0010] Preferably, the proportion of reducible electricity loads satisfies the following relational expression:

[0011]

[0012] In the formula, is the proportion of reducible electricity loads at time in the system, and its value range is [0, 1], is the total number of electrical nodes, is the rated power of the electrical load at node ; is the rated power of the reducible electrical load at node at time ;

[0013] The reduction rate of the reducible electrical load satisfies the following relation:

[0014]

[0015] wherein is the reduction rate of the reducible electrical load at time in the system, and the value range is [0, 1]. 1 means that all reducible electrical loads are operating normally at time , and 0 means that all reducible electrical loads are out of service at time ; is the power of the reducible electrical load at node

[0016] The cumulative reduction times of the reducible electrical load satisfy the following relation:

[0017]

[0018] wherein is the cumulative reduction times of the reducible electrical load at node at time , ; is the sign function; is the evaluation time step of the cumulative reduction times of the reducible electrical load; is the reduction amount of the reducible electrical load at node at time ;

[0019] Among them, the observation interval of the cumulative reduction times of the reducible electrical load is from the starting restoration time of the area to be restored in the system to the set evaluation end time ;

[0020] The proportion of the cumulative reduction amount of the reducible electrical load satisfies the following relation:

[0021]

[0022] wherein is the proportion of the cumulative reduction amount of the reducible electrical load at node at time in the interval , and the value range is [0, 1], is 1, indicating that the reducible electricity load is fully reduced within the interval and is 0, indicating that the reducible electricity load is not reduced within the interval is the number of time periods within the interval and is the reduction rate of the reducible electricity load at node during the time period within the interval is the node at which the reduction rate of the reducible electricity load during the time period

[0023] Preferably, the data of the distributed power source includes: rated power, predicted output value;

[0024] The smoothness of the predicted output value of the distributed power source satisfies the following relational expression:

[0025]

[0026] wherein, is the smoothness of the predicted output value of the distributed power source in the system at time, and the value range is [0, 1]. Being 1 indicates that the predicted output value is constant, and being 0 indicates that the predicted output value fluctuates; is the evaluation time step of the smoothness of the predicted output value; is from time to the set evaluation end time the total number of time periods in between; is the time period serial number; is the node at which the distributed power source at time is the predicted output value; is the rated power output by the distributed power source at node is the set of nodes corresponding to the distributed power source in the system.

[0027] Preferably, the data of the electricity load includes: rated power of the electricity load, power of the electricity load restored with power supply;

[0028] The adequacy of power supply restoration in the area to be restored with power supply in the system satisfies the following relational expression:

[0029]

[0030] wherein, is the adequacy of power supply restoration in the area to be restored with power supply in the system, , and the smaller the value of the lower the adequacy of power supply restoration in the area to be restored with power supply, and is the node at the power output by the distributed power source at the moment; is the node at the power of the electrical load restored at the moment.

[0031] Preferably, the unit capacity compensation cost model of the electrical load that can be curtailed satisfies the following relationship:

[0032]

[0033] In the formula, is the node at the unit capacity compensation cost of the electrical load that can be curtailed at the moment; is the node the quotation when all the electrical load that can be curtailed at the place is curtailed at one time; is the proportion of the distributed power source in the system, and the value range is [0,1], being 0 means that there is no distributed power source in the system; is the node at the set end moment of the evaluation the cumulative curtailment times of the electrical load that can be curtailed.

[0034] Preferably, collect the data of the gas load that can be curtailed to determine the proportion and curtailment rate of the gas load that can be curtailed, and determine the cumulative curtailment times and the proportion of the cumulative curtailment volume of the gas load that can be curtailed at each node in the system; collect and utilize the data of the gas load in the system and the data of the gas source in the system to determine the adequacy of gas supply restoration in the area to be restored;

[0035] Use the proportion and curtailment rate of the gas load that can be curtailed, the cumulative curtailment times and the proportion of the cumulative curtailment volume of the gas load that can be curtailed at each node in the system, and the adequacy of gas supply restoration in the area to be restored to establish a unit capacity compensation cost model for the gas load that can be curtailed;

[0036] Among them, the data of the gas load that can be curtailed includes: curtailment volume, rated power, power.

[0037] Preferably, the proportion of the gas load that can be curtailed satisfies the following relationship:

[0038]

[0039] In the formula, is the proportion of the gas load that can be curtailed in the system at the moment, and the value range is [0,1], is the total number of gas nodes, is the node the rated power of the gas load at the place; is the rated power of the gas consumption load that can be reduced at node at time;

[0040] The reduction rate of the gas consumption load that can be reduced satisfies the following relationship:

[0041]

[0042] where is the reduction rate of the gas consumption load that can be reduced in the system at time, and its value range is [0, 1]. A value of 1 means all the gas consumption loads that can be reduced in the system are operating normally at time, and a value of 0 means all the gas consumption loads that can be reduced in the system are shut down at time; is the power of the gas consumption load that can be reduced at node

[0043] The cumulative reduction times of the gas consumption loads that can be reduced at each node in the system satisfy the following relationship:

[0044]

[0045] where is the cumulative reduction times of the gas consumption load that can be reduced at node at time, ; is the sign function; is the evaluation time step of the cumulative reduction times of the gas consumption load that can be reduced; is the reduction amount of the gas consumption load that can be reduced at node at time;

[0046] Among them, the observation interval of the cumulative reduction times of the gas consumption load that can be reduced is from the starting restoration time of the area to be restored in the system to the set evaluation end time ;

[0047] The proportion of the cumulative reduction amount of the gas consumption loads that can be reduced at each node in the system satisfies the following relationship:

[0048]

[0049] where is the proportion of the cumulative reduction amount of the gas consumption load that can be reduced at node at time within the interval , and its value range is [0, 1], It is 1, indicating that the reducible gas load is completely reduced within the interval and it is 0, indicating that the reducible gas load is not reduced within the interval ; is the number of time periods within the interval , and is the reduction rate of the reducible gas load at node at time period;

[0050] The gas supply restoration adequacy of the area to be restored in the system satisfies the following relationship:

[0051]

[0052] In the formula, is the gas supply restoration adequacy of the area to be restored in the system, , and the smaller the value of , the lower the gas supply restoration adequacy of the area to be restored, and the larger the value of , the higher the gas supply restoration adequacy of the area to be restored; is the power output by the distributed power source at node at time; is the power of the electricity load restored at node at

[0053] Preferably, the unit capacity compensation cost model of the reducible gas load satisfies the following relationship:

[0054]

[0055] In the formula, is the unit capacity compensation cost of the reducible gas load at node at time; is the quotation when the reducible gas load at node is completely reduced at one time; is the cumulative reduction times of the reducible gas load at node at the set evaluation end time ;

[0056] Preferably, based on the unit capacity compensation cost model of the reducible electricity load, evaluate the reducible electricity load compensation of each partition in the system in different energy supply restoration plans; based on the unit capacity compensation cost model of the reducible gas load, evaluate the reducible gas load compensation of each partition in the system in different energy supply restoration plans.

[0057] Preferably, the compensation cost for the load that can be curtailed in the th step of the power restoration plan satisfies the following relationship:

[0058]

[0059] In the formula, is the compensation cost for the load that can be curtailed in the th step of the power restoration plan, is the moment when the power load is restored in the th step of the power restoration plan, and are respectively the unit capacity compensation cost and the compensation capacity of the user for the load that can be curtailed at the moment in the th step of the power restoration plan, is the set evaluation end moment.

[0060] Preferably, the compensation cost for the load that can be curtailed in the th step of the gas supply restoration plan satisfies the following relationship:

[0061]

[0062] In the formula, is the compensation cost for the load that can be curtailed in the th step of the power supply restoration plan, is the moment when the gas load is restored in the th step of the gas supply restoration plan, and are respectively the unit capacity compensation cost and the compensation capacity of the user for the load that can be curtailed at the moment in the th step of the gas supply restoration plan, is the set evaluation end moment.

[0063] Preferably, within each zone, the net income of the power restoration plan is composed of the income of the power restoration plan, the power supply cost, the compensation cost for the load that can be curtailed, the transportation cost of the mobile power source, and the risk cost; the net income of the gas supply restoration plan is composed of the income of the gas supply restoration plan, the gas supply cost, the compensation cost for the load that can be curtailed, the transportation cost of the mobile gas source, and the risk cost.

[0064] Preferably, the sum of the net income of the power restoration plan and the net income of the gas supply restoration plan within a zone is maximized as the optimization objective function of the lower-level model, satisfying the following relationship:

[0065]

[0066] Wherein, is the sum of the net benefits of the power supply restoration plan and the net benefits of the gas supply restoration plan in the partition within the system, and are respectively the number of steps of the power supply restoration plan and the number of steps of the gas supply restoration plan in the partition during the same decision-making period; and are respectively the benefit of restoring power supply for the electrical load and the benefit of restoring gas supply for the gas load in the th step of the power supply restoration plan and the th step of the gas supply restoration plan; and are respectively the power supply cost of the power source and the gas supply cost of the gas source in the th step of the power supply restoration plan and the th step of the gas supply restoration plan; and are respectively the compensation cost for the cuttable electrical load and the compensation cost for the cuttable gas load in the th step of the power supply restoration plan and the th step of the gas supply restoration plan; and are respectively the transportation cost of the mobile power source and the transportation cost of the mobile gas source in the th step of the power supply restoration plan and the th step of the gas supply restoration plan; and are respectively the risk cost of the electrical load and the risk cost of the gas load in the th step of the power supply restoration plan and the th step of the gas supply restoration plan.

[0067] Preferably, the risk cost of the electrical load in the th step of the power supply restoration plan satisfies the following relational expression:

[0068]

[0069] Wherein, is the unit capacity restoration risk of the electrical load at the moment in the th step of the power supply restoration plan; and are respectively the unit restoration benefit and the restoration amount of the restored electrical load at time t in the th step of the power supply restoration plan;

[0070] The risk cost of the gas load in the th step of the gas supply restoration plan satisfies the following relational expression:

[0071]

[0072] In the formula, is the unit capacity restoration risk of the gas load at the th step of the gas supply restoration plan; The unit capacity restoration risk of the gas load at the , are respectively the unit restoration revenue and restoration volume of the gas load restored at the th step of the gas supply restoration plan at time t.

[0073] Preferably, the constraint conditions of the lower-layer objective function include: power supply system operation and demand response constraints, gas supply system operation and demand response constraints.

[0074] Preferably, the optimization objective function of the upper-layer model is that the sum of the net revenues of the power supply restoration plan and the gas supply restoration plan in each partition of the system is maximized, and the following relationship is satisfied:

[0075]

[0076] In the formula, is the net revenue of the system's energy supply restoration; is the total number of partitions containing electrical loads, and each partition containing an electrical load has at least one power source under the V / F control mode; is the sum of the net revenue of power supply restoration of the distribution system and the net revenue of gas supply restoration of the gas distribution system in partition of the system, .

[0077] Preferably, the constraint conditions of the upper-layer objective function include: radial topology constraints, system operation constraints, mobile power source and mobile gas source constraints.

[0078] Preferably, when solving the lower-layer model, a delay value characterizing the delay characteristics of gas load transmission is set, and the following relationship is satisfied:

[0079]

[0080] In the formula, is the delay value, , are respectively the length and diameter of the natural gas pipeline in the partition , is the gas load supply at time is the delay value at the inlet end where the gas flows into the natural gas pipeline from the previous pipeline. When there is gas flowing into the inlet end of the natural gas pipeline from multiple previous pipelines, the average value of each delay value is taken.

[0081] The present invention also provides an energy supply restoration system for an integrated energy system considering demand response compensation evaluation, including:

[0082] A compensation evaluation module, configured to establish a unit capacity compensation cost model for the load of reducible energy consumption to evaluate the compensation for the load of reducible energy consumption in each partition within the system in different energy supply restoration schemes;

[0083] An energy supply restoration net income module, configured to use the income, cost, evaluation value of the compensation for the load of reducible energy consumption, transportation cost of mobile energy, and risk cost in each partition in different energy supply restoration schemes to form the net income of energy supply restoration for each partition;

[0084] A partition and energy supply restoration optimization module, configured to be a system energy supply restoration optimization module based on a two-layer model structure; wherein, the upper-layer model determines the partition scheme with the goal of maximizing the total sum of the net income of energy supply restoration for all partitions, and optimizes the partition scheme according to the real-time scheduling of mobile energy; the lower-layer model, based on the optimized partition scheme, uses the unit capacity compensation cost model for the load of reducible energy consumption to re-evaluate the compensation for the load of reducible energy consumption in each optimized partition in different energy supply restoration schemes, and determines the energy supply restoration scheme for each partition with the goal of maximizing the updated net income of energy supply restoration for each partition; the lower-layer model feeds back the updated net income of energy supply restoration for each partition to the upper-layer model; through multiple iterations, the optimal partition of the integrated energy system and the energy supply restoration scheme for each optimal partition are obtained.

[0085] The present invention is also a terminal, including a processor and a storage medium; the storage medium is used to store instructions; the processor is used to operate according to the instructions to execute the steps of the method.

[0086] The present invention is also a computer-readable storage medium, on which a computer program is stored, and when the program is executed by a processor, the steps of the method are implemented.

[0087] The beneficial effects of the present invention are at least as follows compared with the prior art. The present invention proposes a demand response for the power restoration process. Considering the large-scale access of new energy to the distribution system at the present stage, various factors affecting power restoration are taken into account on the basis of traditional demand response, so as to obtain a unit capacity compensation cost model for the reducible electricity load. The present invention comprehensively considers the proportion of the reducible electricity load, the reduction rate, the cumulative reduction times, and the proportion of the cumulative reduction amount to characterize the real situation of the cumulative response of the reducible electricity load during the power supply restoration process in the actual power outage scenario, and reflects the user's own demand response ability and the actual participation willingness of the user; at the same time, the smoothness of the predicted output value of the distributed power source and the adequacy of the restored energy supply in the area to be restored are introduced to characterize the influence of the external environment of the system. The established unit capacity compensation cost model for the reducible electricity load covers factors such as the reducible ability of the contracted demand response electricity load affected by the power outage scenario and the restoration process, and the system restoration adequacy, and can characterize the effect of the comprehensive load demand response in the restoration scenario.

[0088] The objective function established by the two-layer model includes a compensation evaluation characterized by the compensation cost. When facing the energy supply restoration scenario, the effects of different energy supply restoration schemes under the demand response compensation evaluation are fully considered; moreover, the present invention not only proposes an energy supply restoration scheme with the optimal net income for each partition (the restoration time, restoration amount of the load within each partition, the response time and response amount of the demand response, etc.), but also gives the optimal partition scheme to support the best energy supply restoration effect for each partition, so as to jointly form an integrated energy supply restoration scheme for the comprehensive energy system with the system optimal partition scheme and the optimal energy supply restoration scheme for each partition. It not only meets the overall management of the states of each partition from the perspective of the whole system and the mobile energy scheduling instructions issued by the superior dispatching in real time, but also details the implementation steps of the energy supply restoration within each partition. Brief Description of the Drawings

[0089] Figure 1 is a flowchart of the method for restoring the energy supply of the integrated energy system considering the demand response compensation evaluation proposed by the present invention;

[0090] Figure 2 is a schematic structural diagram of the two-layer optimization model in the embodiment of the present invention;

[0091] Figure 3 is a flowchart of the iterative solution of the lower-layer model in the embodiment of the present invention. Detailed Embodiment

[0092] To make the objectives, technical solutions, and advantages of the present invention clearer, the following will clearly and completely describe the technical solutions of the present invention with reference to the accompanying drawings in the embodiments of the present invention. The embodiments described in this application are only a part of the embodiments of the present invention, rather than all embodiments. Based on the spirit of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts belong to the protection scope of the present invention.

[0093] The present invention proposes a power supply restoration method for an integrated electrical energy system considering the compensation evaluation of comprehensive load demand response. The demand response is embodied as the reducible energy load. In the embodiment, the incentive-based demand response participates in the interaction with the power grid in the form of reducible electricity load.

[0094] As Figure 1 shown, the method includes:

[0095] Step 1: Establish a unit capacity compensation cost model for the reducible energy load to evaluate the compensation of the reducible energy load in each partition of the system in different power supply restoration schemes.

[0096] The reducible energy load of the integrated electrical energy system includes: reducible electricity load and reducible gas load; therefore, Step 1 includes:

[0097] Step 1.1: Collect data on the reducible electricity load to determine the proportion and reduction rate of the reducible electricity load, and determine the cumulative reduction times and the proportion of the cumulative reduction amount of the reducible electricity load at each node in the system; collect data on distributed power sources to determine the smoothness of the predicted output value of the distributed power sources; collect and utilize the data of the electricity load and the distributed power sources to determine the adequacy of power supply restoration in the area to be restored; use the proportion and reduction rate of the reducible electricity load, the cumulative reduction times and the proportion of the cumulative reduction amount of the reducible electricity load at each node in the system, the smoothness of the predicted output value of the distributed power sources, and the adequacy of power supply restoration in the area to be restored to establish a unit capacity compensation cost model for the reducible electricity load.

[0098] Specifically, Step 1.1 includes:

[0099] Step 1.1.1: Collect data on the reducible electricity load in the system, including but not limited to: the reduction amount of the reducible electricity load and the rated power of the reducible electricity load.

[0100] The data of the reducible electricity load satisfies the following relational expression:

[0101]

[0102]

[0103] In the formula, is the power load reduction status at node at moment. 1 indicates being reduced, and 0 indicates not being reduced; is the reduction amount of the reducible power load at node at moment; is the rated power of the reducible power load at node at moment; is the reduction rate of the reducible power load at node at moment;

[0104] In the embodiment, the power load at node at moment is a contracted power load, which is a non - restrictive and better choice;

[0105] Step 1.1.2, the proportion of the reducible power load in the system satisfies the following relational expression:

[0106]

[0107] In the formula, is the proportion of the reducible power load in the system at moment, and its value range is [0, 1], is the total number of electrical nodes, is the rated power of the power load at node ;

[0108] In the embodiment, the proportion of the reducible power load is used to represent the proportion of the rated power of all reducible power loads with which the demand response agreement has been signed in the total rated power of all power loads in the system;

[0109] Step 1.1.3, collect data of the reducible power load in the system, including but not limited to: the power of the reducible power load;

[0110] The reduction rate of the reducible power load in the system satisfies the following relational expression:

[0111]

[0112] In the formula, is the reduction rate of the reducible power load in the system at moment, and its value range is [0, 1]. 1 indicates that at moment, all reducible power loads in the system are operating normally, and 0 indicates that at moment, all reducible power loads in the system are out of service; is at node at The power of the electricity consumption load that can be reduced at a moment;

[0113] In the embodiment, the reduction rate of the electricity consumption load that can be reduced in the system is used to characterize the proportion of the actually reducible electricity consumption load in the system among the rated powers of all the reducible electricity consumption loads for which the demand response agreement has been signed;

[0114] Step 1.1.4, the cumulative reduction times of the electricity consumption load that can be reduced at each node in the system satisfy the following relational expression:

[0115]

[0116] In the formula, is the cumulative reduction times of the electricity consumption load that can be reduced at node at moment, ; is the sign function; is the evaluation time step of the cumulative reduction times of the electricity consumption load that can be reduced; is node at the reduction amount of the electricity consumption load that can be reduced at moment;

[0117] In the embodiment, the electricity consumption load that can be reduced at node at moment is a signed reducible electricity consumption load and is a non - restrictive better choice;

[0118] The observation interval of the cumulative reduction times of the electricity consumption load that can be reduced is from the starting restoration moment of the area to be restored in the system to the set evaluation end moment ;

[0119] Step 1.1.5, the proportion of the cumulative reduction amount of the electricity consumption load that can be reduced at each node in the system satisfies the following relational expression:

[0120]

[0121] In the formula, is the proportion of the cumulative reduction amount of the electricity consumption load that can be reduced at node at moment within the interval , and the value range is [0, 1]. being 1 means that the electricity consumption load that can be reduced is completely reduced within the interval , being 0 means that the electricity consumption load that can be reduced is not reduced within the interval ; is the number of time periods within the interval , is node at Reduction rate of the electricity load that can be curtailed during a time period

[0122] In the embodiment, the proportion of the cumulative reduction amount of the electricity load that can be curtailed is used to characterize the duration of power resupply after the reduction of the electricity load that can be curtailed and the influence of each reduction amount

[0123] Step 1.1.6: Collect data of distributed power sources in the system, including but not limited to: rated power and predicted output value of the distributed power sources

[0124] The smoothness of the predicted output value of the distributed power source satisfies the following relationship

[0125]

[0126] In the formula is the smoothness of the predicted output value of the distributed power source in the system at time, and the value range is [0, 1]. 1 means the predicted output value is constant, and 0 means the predicted output value fluctuates is the evaluation time step of the smoothness of the predicted output value is from time to the set evaluation end time total number of time periods between is the time period serial number is the node distributed power source at predicted output value at time is the node rated power output by the distributed power source at is the set of nodes corresponding to the distributed power sources in the system

[0127] In the embodiment, the smoothness of the predicted output value of the distributed power source is used to characterize the output fluctuation degree of the new energy units in the system within the observation interval

[0128] Step 1.1.7: Collect data of the electricity load in the system, including but not limited to: rated power of the electricity load, power of the electricity load restored with power supply

[0129] The adequacy of power restoration for the area to be restored with power supply in the system satisfies the following relationship

[0130]

[0131] In the formula is the adequacy of power restoration for the area to be restored with power supply in the system , and the smaller the value of, the lower the adequacy of power restoration for the area to be restored with power supply The larger the value, the higher the adequacy of power restoration in the area to be restored with power supply. is the node at the power output of the distributed power source at time is the node at the power of the electricity consumption load restored at time

[0132] In the embodiment, the adequacy of power restoration in the area to be restored with power supply in the system is used to characterize the ratio relationship between the callable power of the new energy units with power supply capacity in the area to be restored with power supply in the system and the power demand of the electricity consumption load to be restored with power supply.

[0133] Step 1.1.8, establish a unit capacity compensation cost model for the electricity consumption load that can be curtailed, satisfying the following relational expression:

[0134]

[0135] In the formula, is the unit capacity compensation cost of the electricity consumption load that can be curtailed at the node at time; is the quotation when all the electricity consumption load that can be curtailed at the node is curtailed at one time; is the proportion of the distributed power source in the system, and the value range is [0,1], being 0 means that there is no distributed power source in the system;

[0136] is at time the proportion of the electricity consumption load that can be curtailed in the system, is at time the curtailment rate of the electricity consumption load that can be curtailed in the system; due to the intermittency and volatility of the new energy output, which will affect the available capacity of the system and thus the power restoration ability on the power source side of the system. In some power outage scenarios and the power restoration process, even if the electricity load has signed a demand response agreement, it will not have the curtailment ability due to its outage. Only the signed electricity consumption load that can be curtailed and is operating normally is the actual electricity consumption load that can be curtailed. Therefore, and are used to be corrected to obtain a quotation that matches the power restoration ability on the power source side of the system and the actual electricity consumption load that can be curtailed;

[0137] is the smoothness of the predicted output value of the distributed power source in the system at time; the output volatility of the new energy is also a factor affecting the power restoration process. Therefore, is used to Make corrections to obtain quotes that match the output volatility of new energy;

[0138] For the adequacy of energy restoration in the area to be restored to power in the system; therefore, utilize For Make corrections to obtain quotes that match the adequacy of power grid energy restoration in the area to be restored to power;

[0139] For node At The cumulative reduction times of the load that can be curtailed at a certain moment; For node At the set end moment of the assessment The cumulative reduction times of the load that can be curtailed; For node At The proportion of the cumulative reduction amount of the load that can be curtailed at a certain moment;

[0140] The present invention proposes a demand response for the power restoration process. Considering the large-scale access of new energy to the distribution system at the present stage, various factors affecting power restoration are taken into account on the basis of traditional demand response, so as to obtain a unit capacity compensation cost model for the load that can be curtailed. The present invention comprehensively considers the proportion, reduction rate, cumulative reduction times, and proportion of cumulative reduction amount of the load that can be curtailed to characterize the real situation of the cumulative response of the load that can be curtailed in the power supply restoration process under the actual power outage scenario, and reflects the user's own demand response ability and the user's actual participation willingness; at the same time, the smoothness of the predicted value of the distributed power generation output and the adequacy of energy restoration in the area to be restored to power are introduced to characterize the influence of the external environment of the system. The established unit capacity compensation cost model for the load that can be curtailed covers factors such as the curtailment ability of the contracted demand response load affected by the power outage scenario and the restoration process, and the system restoration adequacy, and can characterize the effect of the comprehensive load demand response in the restoration scenario.

[0141] Step 1.2, collect data on the load that can be curtailed in the system to determine the proportion and reduction rate of the load that can be curtailed in the system, and determine the cumulative reduction times and proportion of cumulative reduction amount of the load that can be curtailed at each node in the system; collect and utilize the data of the gas load in the system and the data of the gas source in the system to determine the adequacy of gas supply restoration in the area to be restored to gas supply; use the proportion and reduction rate of the load that can be curtailed in the system, the cumulative reduction times and proportion of cumulative reduction amount of the load that can be curtailed at each node in the system, and the adequacy of gas supply restoration in the area to be restored to gas supply to establish a unit capacity compensation cost model for the load that can be curtailed;

[0142] Specifically, step 2 includes:

[0143] Step 1.2.1, collect the data of the gas load that can be reduced in the system, including but not limited to: the reduction amount of the gas load that can be reduced, the rated power of the gas load that can be reduced, and the power of the gas load that can be reduced;

[0144] The data of the gas load that can be reduced satisfies the following relational expression:

[0145]

[0146]

[0147] In the formula, is the gas load reduction state at node at time , where 1 means reduction and 0 means no reduction; is the reduction amount of the gas load that can be reduced at node at time ; is the rated power of the gas load that can be reduced at node at time ; is the reduction rate of the gas load that can be reduced at node at time ;

[0148] Step 1.2.2, the proportion of the gas load that can be reduced in the system satisfies the following relational expression:

[0149]

[0150] In the formula, is the proportion of the gas load that can be reduced in the system at time , and its value range is [0, 1]. is the total number of gas nodes, is the rated power of the gas load at node ; is the rated power of the gas load that can be reduced at node at time ;

[0151] Step 1.2.3, the reduction rate of the gas load that can be reduced in the system satisfies the following relational expression:

[0152]

[0153] In the formula, is the reduction rate of the gas load that can be reduced in the system at time , and its value range is [0, 1]. 1 means all the gas loads that can be reduced in the system are operating normally at time All the reducible gas loads in the time system are shut down; For node at the power of the reducible gas load at time;

[0154] Step 1.2.4, the cumulative reduction times of the reducible gas loads at each node in the system satisfy the following relationship:

[0155]

[0156] In the formula, is the cumulative reduction times of the reducible gas load at node at time, ; is the sign function; is the evaluation time step of the cumulative reduction times of the reducible gas load; is for node at the reduction amount of the reducible gas load at time;

[0157] Among them, the observation interval of the cumulative reduction times of the reducible gas load is from the starting restoration time of the area to be restored in the system to the set evaluation end time ;

[0158] Step 1.2.5, the proportion of the cumulative reduction amount of the reducible gas loads at each node in the system satisfies the following relationship:

[0159]

[0160] In the formula, is the proportion of the cumulative reduction amount of the reducible gas load at node at time within the interval , and the value range is [0,1], being 1 means that the reducible gas load is completely reduced within the interval , being 0 means that the reducible gas load is not reduced within the interval ; is the number of time periods within the interval , is the reduction rate of the reducible gas load at node at time period;

[0161] Step 1.2.6, the restoration gas supply adequacy of the area to be restored with gas supply in the system satisfies the following relationship:

[0162]

[0163] In the formula, is the restoration gas supply adequacy of the area to be restored in the system, and The smaller the value of, the lower the restoration gas supply adequacy of the area to be restored, The larger the value of, the higher the restoration gas supply adequacy of the area to be restored; is the power output of the distributed power source at node at time; is the power of the electricity load restored at node at time; is the set of nodes corresponding to the gas sources in the system;

[0164] Step 1.2.7, the unit capacity compensation cost model of the gas load that can be curtailed, satisfies the following relationship:

[0165]

[0166] In the formula, is the unit capacity compensation cost of the gas load that can be curtailed at node at time; is the quotation when all the gas load that can be curtailed at node is curtailed at once; is the cumulative curtailment times of the gas load that can be curtailed at node at the set evaluation end time ;

[0167] The characteristics of the compensation cost model of the gas load that can be curtailed are the same as those of the compensation cost model of the electricity load that can be curtailed.

[0168] is the proportion of the gas load that can be curtailed in the system at time, is the curtailment rate of the gas load that can be curtailed in the system at time. In some gas outage scenarios and restoration processes, even if the gas load has signed a demand response agreement, it will not have the curtailment ability due to its outage. Only the signed and curtailable gas load in normal operation is the actual curtailable gas load. Therefore, and are used to for correction to obtain a quotation that matches the actual curtailable gas load of the system;

[0169] is the restoration gas supply adequacy of the area to be restored in the system. Therefore, is used to Make corrections to obtain an offer that matches the adequacy of energy restoration in the area where gas supply is to be restored;

[0170] The present invention proposes a demand response for the gas supply restoration process. Considering the large-scale access of new energy to the distribution system at the present stage, various factors affecting gas supply restoration are taken into account on the basis of traditional demand response, so as to obtain a unit capacity compensation cost model for the reducible gas load. The present invention comprehensively considers the proportion of the reducible gas load, the reduction rate, the cumulative reduction times, and the proportion of the cumulative reduction amount to characterize the real situation of the cumulative response of the reducible gas load during the gas supply restoration process in the actual gas outage scenario, and reflects the user's own demand response ability and the user's actual participation willingness; at the same time, the smoothness of the predicted output of distributed power sources and the adequacy of energy restoration in the area where gas supply is to be restored are introduced to characterize the influence of the external system environment. The established unit capacity compensation cost model for the reducible gas load covers factors such as the reducible capacity of the contracted demand response gas load affected by the gas outage scenario and the restoration process, and the system restoration adequacy, and can characterize the effect of the comprehensive load demand response in the restoration scenario.

[0171] Step 1.3: Based on the unit capacity compensation cost model for the reducible electricity load, evaluate the compensation for the reducible electricity load in each partition of the system in different energy restoration plans; based on the unit capacity compensation cost model for the reducible gas load, evaluate the compensation for the reducible gas load in each partition of the system in different energy restoration plans.

[0172] Evaluating the compensation for the reducible energy load in each partition of the system in different energy restoration plans includes: calculating the compensation cost for the reducible energy load in each partition of the system in different energy restoration plans.

[0173] In the th step of the power restoration plan, the compensation cost for the reducible electricity load satisfies the following relational expression:

[0174]

[0175] In the formula, is the compensation cost for the reducible electricity load in the th step of the power restoration plan, is the moment when the electricity load is restored in the th step of the power restoration plan, and are respectively the unit capacity compensation cost for the reducible electricity load and the compensation capacity of the user at the moment of in the th step of the power restoration plan, is the set evaluation end moment.

[0176] In the In step , the compensation cost for the reducible gas load satisfies the following relationship:

[0177]

[0178] In the formula, is the compensation cost for the reducible gas load in the th step of the power restoration plan, is the moment when the gas load resumes gas supply in the th step of the gas restoration plan, and are respectively the unit capacity compensation cost and the compensation capacity of the user for the reducible gas load at the th step in the gas restoration plan moment, is the set evaluation end moment.

[0179] Step 2: Use the benefits, costs, evaluation values of the compensation for the reducible energy load, transportation costs of mobile energy, and risk costs in each partition of different energy restoration plans to form the net energy restoration benefits of each partition.

[0180] Specifically, within each partition, use the benefits, power supply costs, compensation costs for the reducible power load, transportation costs of mobile power sources, and risk costs of the power restoration plan to form the net benefits of the power restoration plan; use the benefits, gas supply costs, compensation costs for the reducible gas load, transportation costs of mobile gas sources, and risk costs of the gas restoration plan to form the net benefits of the gas restoration plan.

[0181] Step 3: Based on the two-layer model structure, establish an optimization module for the system energy restoration; among them, the upper-layer model determines the partition plan with the goal of maximizing the sum of the net energy restoration benefits of all partitions, and optimizes the partition plan according to the real-time scheduling of mobile energy; the lower-layer model, based on the optimized partition plan, uses the unit capacity compensation cost model of the reducible energy load to re-evaluate the compensation for the reducible energy load in each optimized partition of different energy restoration plans, and determines the energy restoration plan for each partition with the goal of maximizing the updated net energy restoration benefits of each partition; the lower-layer model feeds back the updated net energy restoration benefits of each partition to the upper-layer model; through multiple iterations, obtain the optimal partition of the integrated energy system and the energy restoration plans for each optimal partition.

[0182] Specifically, step 3 includes:

[0183] Step 3.1: Maximize the sum of the net benefits of the power restoration plan and the net benefits of the gas restoration plan within a partition as the optimization objective function of the lower-layer model, which satisfies the following relationship:

[0184]

[0185] In the formula, and are respectively the number of steps of the power supply restoration plan and the number of steps of the gas supply restoration plan in the same decision-making period for the sub-region ; and are respectively the benefit of restoring the power load in the -th step of the power supply restoration plan and the benefit of restoring the gas load in the -th step of the gas supply restoration plan; and are respectively the power supply cost of the power source in the -th step of the power supply restoration plan and the gas supply cost of the gas source in the -th step of the gas supply restoration plan; and are respectively the compensation cost for the reducible power load in the -th step of the power supply restoration plan and the compensation cost for the reducible gas load in the -th step of the gas supply restoration plan; and are respectively the transportation cost of the mobile power source in the -th step of the power supply restoration plan and the transportation cost of the mobile gas source in the -th step of the gas supply restoration plan; and are respectively the risk cost of the power load in the -th step of the power supply restoration plan and the risk cost of the gas load in the -th step of the gas supply restoration plan.

[0186] Among them, the risk cost of the power load in the -th step of the power supply restoration plan satisfies the following relational expression:

[0187]

[0188] In the formula, is the unit capacity restoration risk of the power load at the -th step of the power supply restoration plan at the moment; , are respectively the unit restoration benefit and the restoration volume of the restored power load at the -th step of the power supply restoration plan at the t moment.

[0189] The risk cost of the gas load in the -th step of the gas supply restoration plan satisfies the following relational expression:

[0190]

[0191] In the formula, is the In step the unit capacity restoration risk of gas load at a certain moment; and are respectively the unit restoration benefit and restoration volume of the gas load restored in the step at a certain moment.

[0192] Step 3.2. Establish the constraint conditions of the lower - layer objective function, including: power supply system operation and demand response constraints, gas supply system operation and demand response constraints;

[0193] 1) Power supply system operation and demand response constraints

[0194] For the lower - layer objective function, the system power balance constraint needs to be adjusted, and at the same time, the constraints of the demand response load need to be satisfied, satisfying the following relational expressions:

[0195]

[0196]

[0197]

[0198] In the formula, is the set of energized node of the distribution system in the sub - region; and are respectively the active response amount and reactive response amount of the demand response load reduced at node at a certain moment; is the maximum value of the cumulative reduction times of the demand response electrical load at node ; and are respectively the actual reducible electrical load at node at the th reduction during the observation interval and the lower limit of the duration of power supply restoration after the reduction, , is the cumulative reduction times of the reducible electrical load at node at the set evaluation end moment;

[0199] 2) Gas supply system operation and demand response constraints

[0200] In GDS, in addition to adjusting the node flow balance constraint, the constraints of the demand response load need to be satisfied, satisfying the following relational expressions:

[0201]

[0202]

[0203]

[0204] In the formula, is the set of gas acquisition nodes of the gas supply system in sub-region a; is the node at time the demand response amount of the load at; is the node the maximum cumulative reducible times of the signed demand response gas load at; , are respectively the duration and its lower limit of the resumed gas supply after the actual reducible gas load at the node during the observation interval at the time for the th reduction, .

[0205] Step 3.3, obtain the net income of the power supply restoration plan of the distribution system and the net income of the gas supply restoration plan of the gas distribution system in each sub-region of the system. Taking the maximum total net income as the upper-layer objective function, the following relationship is satisfied:

[0206]

[0207] In the formula, is the net income of the system's energy supply restoration; is the total number of sub-regions with electrical loads, and each sub-region with an electrical load has at least one power source under the V / F control mode; is the sum of the net income of the power supply restoration of the distribution system and the net income of the gas supply restoration of the gas distribution system in sub-region in the system, .

[0208] Step 3.4, establish the constraint conditions of the upper-layer objective function, including: radial topology constraint, system operation constraint, mobile power source and mobile gas source constraint;

[0209] 1), Radial topology constraint

[0210]

[0211] In the formula, , are respectively the number of nodes and the number of lines in sub-region a.

[0212] 2), System operation constraint

[0213]

[0214]

[0215] In the formula, and are respectively the sets of nodes where the electrical load flows into and out of node ; The active power and reactive power injected by the distributed power source at node i are respectively and ; and The active power and reactive power injected by the mobile power source (MPS) from the MPS candidate connection point are respectively and ; The active power and reactive power flowing from node to node are respectively and ; The active power and reactive power flowing from node to node are respectively and ; The active power and reactive power restored at node are respectively ;

[0216]

[0217] In the formula, and are respectively the voltages of the two ends of the line at nodes and ; and are respectively the resistance and reactance of the line ;

[0218]

[0219] In the formula, and are respectively the lower limit and upper limit of the voltage at node . In the embodiment, they are respectively taken as 0.9 p.u. and 1.1 p.u.

[0220]

[0221] In the formula, and are respectively The inflow gas load and the node set of the outflow gas load ; The gas flow rate injected by the gas source at the node ; The gas flow rate injected by the mobile gas source (MGS) from the MGS candidate connection point ; The gas flow rate from the node flowing to the node ; The gas flow rate from the node flowing to the node ; The gas flow rate recovered at the node ;

[0222]

[0223]

[0224] In the formula, , , are the pressure value, upper limit, and lower limit of node m respectively; , are the upper limit and lower limit of the flow rate of pipeline segment mn respectively.

[0225] 3), Constraints of mobile power sources and mobile gas sources

[0226]

[0227]

[0228]

[0229] In the formula, is the total number of MPSs, is the total number of MGSs; is the set of MPS candidate connection points in the power grid, is the set of MGS candidate connection points in the gas grid; , are the access status variables of the mobile power source to the MPS candidate connection point , the gas source to the MGS candidate connection point , 1 means access, 0 means no access; is the shortest transportation path of the mobile power source from its storage location to its candidate connection point; is the length of the road section , It is the set of all radial paths of the transportation network.

[0230] In the embodiment, the transportation paths of the mobile power source and the mobile gas source are constrained to be open loops.

[0231] Step 3.5: The upper-layer model determines the partition scheme with the goal of maximizing the total net benefit of energy supply recovery for all partitions, and optimizes the partition scheme according to the real-time scheduling of mobile energy; the lower-layer model, based on the optimized partition scheme, uses the unit capacity compensation cost model of the load that can be curtailed to re-evaluate the compensation for the load that can be curtailed in each optimized partition in different energy supply recovery schemes, and determines the energy supply recovery scheme for each partition with the goal of maximizing the net benefit of energy supply recovery for each updated partition; the lower-layer model feeds back the net benefit of energy supply recovery for each updated partition to the upper-layer model; through multiple iterations, the optimal partition of the integrated energy system and the energy supply recovery scheme for each optimal partition are obtained.

[0232] In the present invention, as Figure 2 shown, the objective functions established by the two-layer model include the compensation evaluation characterized by the compensation cost. When facing the energy supply recovery scenario, the effects of different energy supply recovery schemes under the demand response compensation evaluation are fully considered; moreover, the present invention not only proposes the energy supply recovery scheme with the optimal net benefit for each partition (the recovery time, recovery volume, response time and response volume of demand response, etc. within each partition), but also gives the optimal partition scheme to support the best energy supply recovery effect for each partition. Thus, the energy supply recovery comprehensive scheme of the integrated energy system is jointly composed of the system-optimal partition scheme and the optimal energy supply recovery scheme for each partition, which not only satisfies the overall management of the states of each partition from the perspective of the whole system and the real-time mobile energy scheduling instructions issued by the superior dispatching, but also details the implementation steps of energy supply recovery within each partition.

[0233] When solving the lower-layer model, according to the gas load transmission delay characteristics of the integrated energy system, in cooperation with the iteration step size and iteration node control parameters, the decoupling of the power supply recovery scheme and the gas supply recovery scheme is realized; among them, a delay value representing the gas load transmission delay characteristics is set to satisfy the following relational expression:

[0234]

[0235] In the formula, is the delay value, , are the length and diameter of the natural gas pipeline within the partition respectively, is the gas load supply at time, is the delay value when the gas flows into the inlet end of the natural gas pipeline from the previous pipeline, and when there is gas flowing into the inlet end of the natural gas pipeline from multiple previous pipelines, Take the average value of each delay value;

[0236] As Figure 3 shown, at the beginning of the lower-layer model solution, obtain the initialization partition parameters sent by the upper-layer model, and let , , observation time ; determine the simulation step size of the power distribution network PDS, and update the available capacity of the gas distribution network GDS; when there is no available capacity in the GDS and the delay value , then perform electrical coupling information update; when there is available capacity in the GDS, optimize the GDS by one time step , solve the gas supply recovery plan and the delay value , simulate and deduce the gas supply recovery plan, and let , perform electrical coupling information update;

[0237] Set , update the available power of the PDS. When there is no available capacity in the PDS and the delay value , let the observation time , compare with . When , end the solution, otherwise update the available capacity of the GDS again; when there is available capacity in the PDS, optimize the PDS by one time step , solve the power supply recovery plan, simulate and deduce the power supply recovery plan, and let , . If and , let the observation time . If , then update the available power of the PDS again.

[0238] The present invention also proposes an integrated energy system power supply recovery system considering demand response compensation evaluation, including:

[0239] A compensation evaluation module for establishing a unit capacity compensation cost model for the load of energy that can be curtailed to evaluate the compensation for the load of energy that can be curtailed in each partition of the system in different power supply recovery plans;

[0240] A net power supply recovery benefit module for using the benefits, costs, evaluation values of the compensation for the load of energy that can be curtailed, transportation costs of mobile energy, and risk costs in each partition in different power supply recovery plans to form the net power supply recovery benefit of each partition;

[0241] The partition and energy supply recovery optimization module is used for the system energy supply recovery optimization module based on a two-layer model structure. Among them, the upper-layer model determines the partition scheme with the goal of maximizing the total net income of energy supply recovery for all partitions, and optimizes the partition scheme according to the real-time scheduling of mobile energy. The lower-layer model, based on the optimized partition scheme, uses the unit capacity compensation cost model of the load with reducible energy consumption to re-evaluate the compensation of the load with reducible energy consumption in each optimized partition in different energy supply recovery schemes, and determines the energy supply recovery scheme for each partition with the goal of maximizing the net income of energy supply recovery for each updated partition. The lower-layer model feeds back the net income of energy supply recovery for each updated partition to the upper-layer model. Through multiple iterations, the optimal partition of the integrated energy system and the energy supply recovery scheme for each optimal partition are obtained.

[0242] The present disclosure may be a system, a method, and / or a computer program product. The computer program product may include a computer-readable storage medium having thereon computer-readable program instructions for causing a processor to implement various aspects of the present disclosure.

[0243] The computer-readable storage medium may be a tangible device that can retain and store instructions for use by an instruction execution device. The computer-readable storage medium may be, for example, but not limited to, an electrical storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. More specific examples (a non-exhaustive list) of the computer-readable storage medium include: a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), a static random access memory (SRAM), a portable compact disc read-only memory (CD-ROM), a digital versatile disc (DVD), a memory stick, a floppy disk, a mechanically encoded device such as a punch card or raised structures in grooves having instructions stored thereon, and any suitable combination of the foregoing. The computer-readable storage medium as used herein is not construed as an instantaneous signal itself, such as a radio wave or other freely propagating electromagnetic wave, an electromagnetic wave propagated through a waveguide or other transmission medium (e.g., an optical pulse through an optical fiber cable), or an electrical signal transmitted through a wire.

[0244] The computer-readable program instructions described herein may be downloaded from a computer-readable storage medium to various computing / processing devices, or may be downloaded to an external computer or an external storage device through a network, such as the Internet, a local area network, a wide area network, and / or a wireless network. The network may include a copper transmission cable, an optical fiber transmission, a wireless transmission, a router, a firewall, a switch, a gateway computer, and / or an edge server. The network adapter or network interface in each computing / processing device receives the computer-readable program instructions from the network and forwards the computer-readable program instructions for storage in the computer-readable storage medium in each computing / processing device.

[0245] The computer program instructions for performing the operations of the present disclosure may be assembly instructions, instruction set architecture (ISA) instructions, machine instructions, machine-related instructions, microcode, firmware instructions, state-setting data, or source code or object code written in any combination of one or more programming languages, including object-oriented programming languages such as Smalltalk, C++, etc., and conventional procedural programming languages such as the "C" language or similar programming languages. The computer-readable program instructions may be executed entirely on the user's computer, partially on the user's computer, executed as a stand-alone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., through the Internet using an Internet service provider). In some embodiments, by using the state information of the computer-readable program instructions to customize an electronic circuit, such as a programmable logic circuit, a field-programmable gate array (FPGA), or a programmable logic array (PLA), the electronic circuit can execute the computer-readable program instructions to implement various aspects of the present disclosure.

[0246] 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 replacements, and any modification or equivalent replacement 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 energy supply restoration of an integrated energy system considering demand response compensation evaluation, characterized in that, Including: The reducible energy consumption load of the system includes: reducible electricity load and reducible gas load; establishing a unit capacity compensation cost model for the reducible energy consumption load to evaluate the compensation of the reducible energy consumption load in each zone in different energy supply restoration schemes, including: collecting data on the reducible electricity load to determine the proportion and reduction rate of the reducible electricity load, determining the cumulative reduction times and the proportion of the cumulative reduction amount of the reducible electricity load at each node in the system; collecting data on distributed power sources to determine the smoothness of the predicted output value of the distributed power sources; collecting and using the data of the electricity load and the data of the distributed power sources to determine the adequacy of the energy supply restoration in the area to be restored with power supply; using the proportion and reduction rate of the reducible electricity load, the cumulative reduction times and the proportion of the cumulative reduction amount of the reducible electricity load at each node in the system, the smoothness of the predicted output value of the distributed power sources, and the adequacy of the energy supply restoration in the area to be restored with power supply to establish a unit capacity compensation cost model for the reducible electricity load; collecting data on the reducible gas load to determine the proportion and reduction rate of the reducible gas load, determining the cumulative reduction times and the proportion of the cumulative reduction amount of the reducible gas load at each node in the system; collecting and using the data of the gas load in the system and the data of the gas sources in the system to determine the adequacy of the energy supply restoration in the area to be restored with gas supply; using the proportion and reduction rate of the reducible gas load, the cumulative reduction times and the proportion of the cumulative reduction amount of the reducible gas load at each node in the system, and the adequacy of the energy supply restoration in the area to be restored with gas supply to establish a unit capacity compensation cost model for the reducible gas load; Using the benefits, costs, compensation for the reducible energy consumption load, transportation costs and risk costs of mobile energy in each zone in different energy supply restoration schemes to form the net energy supply restoration benefit of each zone; establishing an energy supply restoration optimization module for the system based on a two-layer model structure; the upper model determines the zoning scheme with the goal of maximizing the total sum of the net energy supply restoration benefits of all zones and optimizes the zoning scheme with real-time scheduling of mobile energy; the lower model, based on the optimized zoning scheme, uses the unit capacity compensation cost model for the reducible energy consumption load to re-evaluate the compensation of the reducible energy consumption load in each optimized zone in different energy supply restoration schemes, and determines the energy supply restoration scheme for each zone with the goal of maximizing the updated net energy supply restoration benefit of each zone; the lower model feeds back the updated net energy supply restoration benefit of each zone to the upper model; through iteration, the optimal zoning of the integrated energy system and the energy supply restoration schemes for each optimal zone are obtained.

2. The method for restoring the energy supply of an integrated energy system considering demand response compensation evaluation according to claim 1, characterized in that The data of the reducible electricity load includes: reduction amount, rated power, power.

3. The method for restoring the energy supply of an integrated energy system considering demand response compensation evaluation according to claim 2, characterized in that The proportion of the reducible electricity load satisfies the following relational expression: In the formula, is the proportion of the electricity load that can be reduced by the system at moment, and its value range is [0, 1]. is the total number of electrical nodes. is the rated power of the electricity load at node . is the rated power of the electricity load that can be reduced at node at moment. The reduction rate of the reducible electricity load satisfies the following relational expression: In the formula, is the reduction rate of the electricity load that can be reduced by the system at time. The value range is [0, 1]. When it is 1, it means that all the electricity loads that can be reduced are operating normally at time. When it is 0, it means that at node at time, the power of the electricity load that can be reduced The cumulative reduction times of the reducible electricity load satisfies the following relational expression: Wherein, is the node at the cumulative reduction times of the electricity load that can be reduced at the moment, ; is the sign function; is the evaluation time step of the cumulative reduction times of the electricity load that can be reduced; is the node at the reduction amount of the electricity load that can be reduced at the moment; Among them, the observation interval of the cumulative reduction times of the electricity load that can be reduced is from the starting restoration moment of the area to be restored in the system to the set evaluation end moment ; The proportion of the cumulative reduction amount of the reducible electricity load satisfies the following relational expression: In the formula, is the node at the cumulative reduction ratio of the electricity load that can be curtailed at time in the interval, and the value range is [0, 1]. When it is 1, it means that the electricity load that can be curtailed is completely curtailed in the interval ; When it is 0, it means that the electricity load that can be curtailed is not curtailed in the interval ; is the number of time periods in the interval ; is the node at the curtailment rate of the electricity load that can be curtailed in the time period.

4. The method for restoring the energy supply of an integrated energy system considering demand response compensation evaluation according to claim 3, wherein Data of distributed power sources, including: rated power, predicted output value; The smoothness of the predicted output value of the distributed power source satisfies the following relational expression: Wherein, is the smoothness of the predicted output value of the distributed power source in the system at time, and its value range is [0, 1]. A value of 1 indicates that the predicted output value is constant, and a value of 0 indicates that the predicted output value fluctuates; is the evaluation time step of the smoothness of the predicted output value; is from time to the set end time of evaluation the total number of time periods in between; is the time period serial number; is the node at the distributed power source at time predicted output value; is the rated power output by the distributed power source at node ; is the set of nodes corresponding to the distributed power sources in the system.

5. The method for restoring the energy supply of an integrated energy system considering demand response compensation evaluation according to claim 4, wherein Data of the electrical load, including: rated power of the electrical load, power of the electrical load for which power supply is restored; The adequacy of power supply restoration in the area to be restored in the system satisfies the following relational expression: In the formula, is the adequacy of power restoration for the area to be restored in the system, , and The smaller the value of, the lower the adequacy of power restoration for the area to be restored, The larger the value of, the higher the adequacy of power restoration for the area to be restored; is the node at The power output by the distributed power source at time; is the node at The power of the electrical load restored at time.

6. The method for restoring the energy supply of an integrated energy system considering demand response compensation evaluation according to claim 5, wherein The unit capacity compensation cost model of the electrical load that can be curtailed satisfies the following relational expression: In the formula, is the unit capacity compensation cost of the load that can be curtailed at node at time ; is the quotation when all the load that can be curtailed at node is curtailed at once; is the proportion of distributed power sources in the system, and its value range is [0, 1]. being 0 means that there is no distributed power source in the system; is the cumulative curtailment times of the load that can be curtailed at node at the set evaluation end time .

7. The method for restoring the energy supply of an integrated energy system considering demand response compensation evaluation according to claim 1, wherein Data of the gas load that can be curtailed, including: curtailment volume, rated power, power.

8. The method for restoring the energy supply of an integrated energy system considering demand response compensation evaluation according to claim 7, wherein The proportion of the gas load that can be curtailed satisfies the following relational expression: wherein, is the proportion of the gas load that can be reduced by the system at time, and its value range is [0, 1], is the total number of gas nodes, is the rated power of the gas load at node ; is the rated power of the gas load that can be reduced at node at time; The curtailment rate of the gas load that can be curtailed satisfies the following relational expression: In the formula, is the reduction rate of the gas load that can be reduced by the system at moment, and the value range is [0, 1]. A value of 1 means that all the gas loads that can be reduced in the system are operating normally at the moment, and a value of 0 means that all the gas loads that can be reduced in the system are out of service at the moment; is the power of the gas load that can be reduced at node at moment; The cumulative curtailment times of the gas load that can be curtailed at each node in the system satisfies the following relational expression: In the formula, is the node at the cumulative reduction times of the gas load that can be reduced at the moment, ; is the sign function; is the evaluation time step of the cumulative reduction times of the gas load that can be reduced; is the node at the reduction amount of the gas load that can be reduced at the moment; Among them, the observation interval of the cumulative reduction times of the gas load that can be reduced is from the starting restoration moment of the area to be restored in the system to the set evaluation end moment ; The proportion of the cumulative curtailment volume of the gas load that can be curtailed at each node in the system satisfies the following relational expression: In the formula, is the node at the cumulative reduction ratio of the reducible gas load within the interval at time , and the value range is [0, 1]. When it is 1, it means that the reducible gas load is completely reduced within the interval . When it is 0, it means that the reducible gas load is not reduced within the interval . is the number of time periods within the interval , is the node at the reduction rate of the reducible gas load in the time period. The adequacy of gas supply restoration in the area to be restored in the system satisfies the following relational expression: Wherein, is the restoration gas supply adequacy of the area to be restored with gas supply in the system, , and the smaller the value of means the lower the restoration gas supply adequacy of the area to be restored with gas supply, the larger the value of means the higher the restoration gas supply adequacy of the area to be restored with gas supply; is the power output by the distributed power source at node at time; is the power of the electricity consumption load restored with power supply at node at time; is the set of nodes corresponding to the gas sources in the system.​​ 9. The method for restoring the energy supply of an integrated energy system considering demand response compensation evaluation according to claim 8, wherein The unit capacity compensation cost model of the gas load that can be curtailed satisfies the following relational expression: Wherein, is the unit capacity compensation cost of the gas load that can be curtailed at node at time ; the quotation when the gas load that can be curtailed at node is completely curtailed at one time; and the cumulative curtailment times of the gas load that can be curtailed at node at the set evaluation end time is .​ 10. The method for restoring the energy supply of an integrated energy system considering demand response compensation evaluation according to claim 1, wherein Based on the unit capacity compensation cost model of the electrical load that can be curtailed, evaluate the compensation for the electrical load that can be curtailed in each partition of the system in different energy supply restoration schemes; based on the unit capacity compensation cost model of the gas load that can be curtailed, evaluate the compensation for the gas load that can be curtailed in each partition of the system in different energy supply restoration schemes.

11. The method for restoring the energy supply of an integrated energy system considering demand response compensation evaluation according to claim 10, wherein In the step of the power restoration plan, the compensation cost for the load that can be curtailed satisfies the following relationship: Wherein, is the compensation cost for the curtailable electricity load in the th step of the power restoration plan, is the moment when the electricity load is restored in the th step of the power restoration plan, and are respectively the unit capacity compensation cost for the curtailable electricity load and the compensation capacity of the user at the th step in the power restoration plan, and is the set evaluation end time.

12. The method for restoring the energy supply of an integrated energy system considering demand response compensation evaluation according to claim 10, wherein In the step of the gas supply resumption plan, the compensation cost for the reducible gas load satisfies the following relationship: In the formula, is the compensation cost for the reducible gas load in the th step of the power supply restoration plan, is the moment when the gas load resumes gas supply in the th step of the gas supply restoration plan, and are respectively the unit capacity compensation cost for the reducible gas load and the compensation capacity of the user at the th step in the gas supply restoration plan, and is the set evaluation end time.

13. The method for restoring the energy supply of an integrated energy system considering demand response compensation evaluation according to claim 1, wherein In each partition, use the revenue of the power supply restoration scheme, power supply cost, compensation cost for the electrical load that can be curtailed, transportation cost of mobile power sources, and risk cost to constitute the net revenue of the power supply restoration scheme; use the revenue of the gas supply restoration scheme, gas supply cost, compensation cost for the gas load that can be curtailed, transportation cost of mobile gas sources, and risk cost to constitute the net revenue of the gas supply restoration scheme.

14. The energy supply restoration method for an integrated energy system considering demand response compensation evaluation according to claim 13, characterized in that Maximizing the sum of the net benefits of the power supply restoration plan and the net benefits of the gas supply restoration plan within a partition is used as the optimization objective function of the lower-layer model, satisfying the following relationship: wherein, is the sum of the net benefits of the power supply restoration plan and the net benefits of the gas supply restoration plan in the partition in the system, and and are respectively the number of steps of the power supply restoration plan and the number of steps of the gas supply restoration plan in the partition during the same decision-making period; ; and are respectively the benefit of restoring power supply to the electrical load in the -th step of the power supply restoration plan and the benefit of restoring gas supply to the gas load in the -th step of the gas supply restoration plan; and are respectively the power supply cost of the power source in the -th step of the power supply restoration plan and the gas supply cost of the gas source in the -th step of the gas supply restoration plan; and are respectively the compensation cost for the reducible electrical load in the -th step of the power supply restoration plan and the compensation cost for the reducible gas load in the -th step of the gas supply restoration plan; and are respectively the transportation cost of the mobile power source in the -th step of the power supply restoration plan and the transportation cost of the mobile gas source in the -th step of the gas supply restoration plan; and are respectively the risk cost of the electrical load in the -th step of the power supply restoration plan and the risk cost of the gas load in the -th step of the gas supply restoration plan.

15. The energy supply restoration method for an integrated energy system considering demand response compensation evaluation according to claim 14, characterized in that In the step, the risk cost of the electrical load satisfies the following relationship: In the formula, is the unit capacity restoration risk of the electrical load at the th step of the power restoration plan at time; and are respectively the unit restoration benefit and the restoration amount of the electrical load restored in the th step of the power restoration plan at time t; The risk cost of the gas load in the step of the gas supply restoration plan satisfies the following relationship: In the formula, is the unit capacity restoration risk of gas load at the th step of the gas supply restoration plan; The unit restoration income and restoration volume of the gas load restored at the - th step of the gas supply restoration plan at time t are and 16. The energy supply restoration method for an integrated energy system considering demand response compensation evaluation according to claim 14, characterized in that The constraint conditions of the lower-layer objective function include: power supply system operation and demand response constraints, and gas supply system operation and demand response constraints.

17. The energy supply restoration method for an integrated energy system considering demand response compensation evaluation according to claim 14, characterized in that Maximizing the sum of the net benefits of the power supply restoration plan and the net benefits of the gas supply restoration plan within each partition of the system is used as the optimization objective function of the upper-layer model, satisfying the following relationship: Wherein, is the net benefit of the system for restoring power supply; is the total number of zones with electrical loads, and each zone with electrical loads has at least one power supply under the V / F control mode; is the zone in the system is the sum of the net benefit of restoring power supply of the distribution system and the net benefit of restoring gas supply of the gas distribution system within; .

18. The energy supply restoration method for an integrated energy system considering demand response compensation evaluation according to claim 17, characterized in that The constraint conditions of the upper-layer objective function include: radial topology constraints, system operation constraints, mobile power source and mobile gas source constraints.

19. The energy supply restoration method for an integrated energy system considering demand response compensation evaluation according to claim 1, characterized in that When solving the lower-layer model, a delay value characterizing the delay characteristics of gas load transmission is set, satisfying the following relationship: Wherein, is the delay value, , are respectively the length and diameter of the natural gas pipeline within the partition , is the gas load supply at time is the delay value when the gas flows from the previous pipeline into the inlet end of the natural gas pipeline . When there is gas flowing from multiple previous pipelines into the inlet end of the natural gas pipeline , the average value of each delay value is taken.

20. A comprehensive energy system power supply restoration system considering demand response compensation evaluation, characterized in that, Including: A compensation evaluation module for establishing a unit capacity compensation cost model for the reducible energy consumption load to evaluate the compensation for the reducible energy consumption load in each partition of the system in different energy supply restoration plans; The energy load that can be reduced in the system includes: the electricity load that can be reduced and the gas load that can be reduced; collect data on the electricity load that can be reduced to determine the proportion and reduction rate of the electricity load that can be reduced, and determine the cumulative reduction times and the proportion of the cumulative reduction amount of the electricity load that can be reduced at each node in the system; collect data on distributed power sources to determine the smoothness of the predicted output value of distributed power sources; collect and utilize the data of the electricity load and the data of distributed power sources to determine the adequacy of energy supply restoration in the area to be restored with power supply; use the proportion and reduction rate of the electricity load that can be reduced, the cumulative reduction times and the proportion of the cumulative reduction amount of the electricity load that can be reduced at each node in the system, the smoothness of the predicted output value of distributed power sources, and the adequacy of energy supply restoration in the area to be restored with power supply to establish a unit capacity compensation cost model for the electricity load that can be reduced; collect data on the gas load that can be reduced to determine the proportion and reduction rate of the gas load that can be reduced, and determine the cumulative reduction times and the proportion of the cumulative reduction amount of the gas load that can be reduced at each node in the system; collect and utilize the data of the gas load in the system and the data of the gas sources in the system to determine the adequacy of gas supply restoration in the area to be restored with gas supply; use the proportion and reduction rate of the gas load that can be reduced, the cumulative reduction times and the proportion of the cumulative reduction amount of the gas load that can be reduced at each node in the system, and the adequacy of gas supply restoration in the area to be restored with gas supply to establish a unit capacity compensation cost model for the gas load that can be reduced; The energy supply restoration net income module is used to use the income, cost, compensating evaluation value of the energy load that can be reduced, transportation cost of mobile energy, and risk cost of each sub-region in different energy supply restoration plans to constitute the energy supply restoration net income of each sub-region; The sub-region and energy supply restoration optimization module is used to, based on a two-layer model structure, the system energy supply restoration optimization module; among them, the upper-layer model determines the sub-region plan with the goal of maximizing the sum of the energy supply restoration net incomes of all sub-regions, and optimizes the sub-region plan according to the real-time scheduling of mobile energy; the lower-layer model, based on the optimized sub-region plan, uses the unit capacity compensation cost model of the energy load that can be reduced to re-evaluate the compensating of the energy load that can be reduced in the optimized each sub-region in different energy supply restoration plans, and determines the energy supply restoration plan of each sub-region with the goal of maximizing the updated energy supply restoration net income of each sub-region; the lower-layer model feeds back the updated energy supply restoration net income of each sub-region to the upper-layer model; through multiple iterations, the optimal sub-regions of the integrated energy system and the energy supply restoration plans of each optimal sub-region are obtained.

21. A terminal, including a processor and a storage medium; characterized in that: The storage medium is used to store instructions; The processor is used to operate according to the instructions to execute the steps of the method described in any one of claims 1-19.

22. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the program is executed by the processor, it realizes the steps of the method described in any one of claims 1-19.

Images