A method, apparatus, equipment and medium for assessing the resilience of a multi-energy power distribution system

Abstract

This application discloses a method, apparatus, equipment, and medium for assessing the resilience of a multi-energy distribution system. The method constructs cyber-physical coupling constraints for the multi-energy distribution system and combines them with operational constraints. It then performs optimal power flow calculations using a first objective function to obtain the retained load of the multi-energy distribution system during the disaster response period. Next, it constructs a cyber-physical collaborative recovery model and, based on this model and real-time updated operational constraints, obtains the retained load of the multi-energy distribution system at each moment during the disaster recovery period using a second objective function. Finally, it determines the resilience index based on the retained load during the disaster response period and the retained load at each moment during the disaster recovery period. This effectively improves the accuracy of predicting the scope of fault impact and minimizes recovery time to reduce power outage losses.

Description

Technical Field

[0001] This application relates to the field of multi-energy power distribution system technology, and in particular to a method, apparatus, equipment and medium for evaluating the resilience of multi-energy power distribution systems. Background Technology

[0002] The safe and reliable operation of multi-energy distribution systems plays a crucial role in sustainable development. In recent years, frequent extreme natural disasters have threatened the normal operation of multi-energy distribution systems. Furthermore, the strong cyber-physical coupling relationship allows faults to propagate between information and physical systems, leading to a continuous expansion of disaster impacts. Existing research largely focuses on assessing the independent resilience of single energy systems or emphasizes the analysis of physical systems, neglecting the cyber-physical coupling fault mechanism. This makes it difficult for assessment models to accurately capture the propagation path of information faults causing abnormal physical equipment operation, resulting in significant discrepancies between predicted fault impact ranges and actual conditions. Consequently, subsequent dispatching schemes cannot efficiently and accurately allocate resources to reduce fault losses. Moreover, the research fails to consider cyber-physical collaborative recovery during the post-disaster recovery phase, making it impossible to minimize recovery time and reduce power outage losses. Summary of the Invention

[0003] This application provides a method, apparatus, equipment, and medium for assessing the resilience of a multi-energy power distribution system. This addresses the shortcomings of existing technologies where significant discrepancies between the predicted impact range of a fault and the actual situation prevent efficient and accurate resource allocation to mitigate losses. Furthermore, the lack of consideration for cyber-physical collaborative recovery during the post-disaster recovery phase hinders efforts to minimize recovery time and reduce power outage losses. This approach achieves minimal load shedding and rapid load recovery, improving the accuracy of fault impact range prediction and minimizing recovery time to reduce power outage losses.

[0004] In a first aspect, embodiments of this application provide a method for assessing the resilience of a multi-energy power distribution system, including: The information substations of the multi-energy power distribution system are used as information nodes, and the working status of each information node is obtained as the node status. By obtaining the initial opening and closing state of the power line and combining it with the node state of the corresponding information node, the cyber-physical coupling constraints of the multi-energy power distribution system are constructed. The operational constraints of the multi-energy distribution system are obtained, and combined with the cyber-physical coupling constraints, the optimal power flow is calculated through the first objective function to obtain the retained load of the multi-energy distribution system during the disaster response period. Construct physical system recovery constraints and information system recovery constraints, and combine the aforementioned cyber-physical coupling constraints to generate a cyber-physical collaborative recovery model; Based on the aforementioned cyber-physical collaborative recovery model and real-time updated operational constraints, the retained load of the multi-energy distribution system at each moment during the disaster recovery period is obtained through the second objective function; The resilience index is determined based on the retained load during the disaster response period and the retained load at each moment during the disaster recovery period.

[0005] In conjunction with the first aspect, in one possible implementation, the cyber-physical coupling constraints of the multi-energy power distribution system include power supply node output constraints and topology control constraints; The power output constraint of the power supply node is as follows:

[0006]

[0007] In the formula, Indicates the power supply node The corresponding information nodes Indicates the power supply node Corresponding information nodes The node status, Indicates the power supply node Minimum active power output, Indicates the power supply node Maximum active power output Indicates the power supply node Output active power, This represents the set of all power supply nodes in a multi-energy power distribution system. This represents the set of information nodes corresponding to all power supply nodes in a multi-energy power distribution system. The topology control constraints are:

[0008] In the formula, Endpoints and endpoints The initial open / closed state of the power lines between them. Endpoints Corresponding information nodes The node status, Endpoints Corresponding information nodes The node status, This indicates the permissible opening and closing states of power lines during network reconfiguration in a multi-energy distribution system. It represents the set of all power lines in a multi-energy distribution system. This represents the constraint control coefficient. In power lines, Endpoints and endpoints The power lines between them.

[0009] In conjunction with the first aspect, in one possible implementation, the physical system restoration constraints include power line restoration constraints and load restoration constraints; wherein, The power line restoration constraints include:

[0010] In the formula, Endpoints With endpoints The power lines between The state of recovery at any time, Endpoints With endpoints The power lines between The state of recovery at any time, Representation and endpoints , The adjacent upstream or downstream power lines between the power lines The state of recovery at any time, Indicates the restored status of the power line. and Indicates the endpoints on an upstream or downstream power line. Endpoints With endpoints The power lines between The opening and closing state at any moment, It represents the set of all power lines in a multi-energy distribution system. Representation and endpoints , The set of all power lines adjacent to each other. This indicates the number of power line maintenance teams. Indicates the total number of disaster recovery periods; in power lines, Endpoints and endpoints The power lines between them; The load recovery constraints include:

[0011] In the formula, Indicates load node exist The recovery state at any moment, u D Indicates the restored status of power lines. Indicates load node The set of adjacent power lines, Endpoints With endpoints The power lines between The state of recovery at any time, Indicates the restored status of power lines. Indicates load node exist Path to all controllable power sources at all times The opening and closing state, Indicates the path from the load node to all controllable power sources. The first A power line in The opening and closing state at any moment, Represents the set of all load nodes. This represents the set of all power lines along the path. This represents the set of all paths from a load node to a controllable power source; in power lines, Endpoints and endpoints The power lines between them.

[0012] In conjunction with the first aspect, in one possible implementation, the information system recovery constraints include communication link repair constraints, information node operational state recovery constraints, and information physical coupling constraints; wherein, The communication link repair constraints include:

[0013] In the formula, Endpoints With endpoints The recovery status of the communication link between them. Endpoints With endpoints The power lines between The state of recovery at any time, Endpoints With endpoints The communication link between them The state of recovery at any time, This represents the set of all communication links in a multi-energy power distribution system. This indicates the number of communication link maintenance teams. Indicates the restored status of power lines. Indicates the recovery status of the communication link. Indicates the total number of disaster recovery periods; in communication links, Endpoints and endpoints Communication links between them; The constraints for restoring the working state of the information node include:

[0014] In the formula, Endpoints With endpoints The recovery status of the communication link between them. Represents the endpoint at one end of the communication link. The corresponding information nodes are in The state of recovery at any time, Indicates the endpoint at the other end of the communication link. The corresponding information nodes are in The state of recovery at any time, Endpoints The corresponding information nodes are in The state of recovery at any time, Indicates the recovery status of the communication link. Indicates the recovery status of the information node. This represents the set of all communication links in a multi-energy power distribution system. Represents the set of all information nodes. Indicates the total number of disaster recovery periods; in communication links, Endpoints and endpoints Communication links between them; The cyber-physical coupling constraints include:

[0015] In the formula, Indicates the power supply node exist The active power output at all times. Indicates the power supply node exist The active power output at all times. Indicates active power. Indicates the power supply node The corresponding active power regulation speed of the equipment Indicates the duration of the disaster recovery period. Indicates the power supply node The corresponding information nodes are in The state of recovery at any time, This indicates the time required for manual adjustment of output. Indicates the power supply node Minimum active power output, Indicates the power supply node Maximum active power output Endpoints With endpoints The power lines between The opening and closing state at any moment, Represents the endpoint at one end of the communication link. The corresponding information nodes are in The state of recovery at any time, Indicates the endpoint at the other end of the communication link. The corresponding information nodes are in The state of recovery at any time, Endpoints With endpoints The power lines between The opening and closing state at any moment, Represents the set of power supply nodes. Indicates the total number of disaster recovery periods. This represents the set of endpoints of all communication links in a multi-energy distribution system. This represents the constraint control coefficient. In the communication link, Endpoints and endpoints Communication links between them.

[0016] In conjunction with the first aspect, in one possible implementation, the expression of the first objective function is specifically as follows:

[0017] In the formula, This indicates the retained load of the multi-energy power distribution system during each disaster response period; These represent the weighting coefficients for electricity, gas, and heat loads, respectively. and These represent the coefficients for the reduction of gas load and heat load to electrical load, respectively. These represent the retained values ​​of electricity, gas, and heat loads during the disaster response period, respectively; among them, Indicates electrical load, Indicates gas load, Indicates heat load, Represents the first in the power system Each node Represents the first in the natural gas system Each node Represents the first in the thermodynamic system 1 node Represents the set of nodes in a power system. Represents the set of nodes in a natural gas system. This represents the set of nodes in a heating system.

[0018] In conjunction with the first aspect, in one possible implementation, the expression for the second objective function is specifically as follows:

[0019] In the formula, This indicates the retained load of the multi-energy power distribution system during each disaster recovery period; These represent the weighting coefficients for electricity, gas, and heat loads, respectively. and These represent the coefficients for the reduction of gas load and heat load to electrical load, respectively. Let represent the retained values ​​of electricity, gas, and heat loads at time t during the disaster recovery period; where, Indicates electrical load, Indicates gas load, Indicates heat load, Represents the first in the power system Each node Represents the first in the natural gas system Each node Represents the first in the thermodynamic system 1 node Represents the set of nodes in a power system. Represents the set of nodes in a natural gas system. This represents the set of nodes in a heating system.

[0020] In conjunction with the first aspect, in one possible implementation, determining the resilience index based on the retention load during the disaster response period and the retention load at each time point during the disaster recovery period includes: Based on the retained load during the disaster response period and the retained load at each moment during the disaster recovery period, a load retention rate curve is plotted, and the load loss area is determined based on the closed area constructed by the load retention rate curve and the constant 100%. The area of ​​the fault curve is obtained and combined with the area of ​​load loss, and the resilience index is determined by the resilience index calculation formula. The formula for calculating the resilience index is:

[0021] In the formula, Indicators of resilience Represents a set of fault scenarios. Indicates the fault scenario The probability of occurrence Indicates the area of ​​load deficiency. This represents the area of ​​the fault curve.

[0022] Secondly, embodiments of this application provide a multi-energy power distribution system resilience assessment device, comprising: The node status acquisition module is used to take the information substation of the multi-energy power distribution system as information nodes and acquire the working status of each information node as the node status. The cyber-physical coupling constraint construction module is used to obtain the initial opening and closing state of the power line and combine it with the node state of the corresponding information node to construct the cyber-physical coupling constraints of the multi-energy power distribution system. The optimal power flow calculation module is used to obtain the operating constraints of the multi-energy distribution system, and in combination with the information-physical coupling constraints, to perform optimal power flow calculation through the first objective function to obtain the retained load of the multi-energy distribution system during the disaster response period; The collaborative recovery model generation module is used to construct physical system recovery constraints and information system recovery constraints, and to generate a cyber-physical collaborative recovery model by combining the cyber-physical coupling constraints. The module for determining the retained load during the recovery period is used to obtain the retained load of the multi-energy distribution system at each moment during the disaster recovery period based on the cyber-physical collaborative recovery model and real-time updated operating constraints, through a second objective function. The resilience index calculation module is used to determine the resilience index based on the retained load during the disaster response period and the retained load at each time during the disaster recovery period.

[0023] Thirdly, embodiments of this application provide an apparatus comprising: a processor; a memory for storing processor-executable instructions; wherein, when the processor executes the executable instructions, it implements the method as described in the first aspect or any possible implementation of the first aspect.

[0024] Fourthly, embodiments of this application provide a non-volatile computer-readable storage medium, the non-volatile computer-readable storage medium including storage for storing a computer program or instructions that, when executed, cause the method described in the first aspect or any possible implementation of the first aspect to be implemented.

[0025] One or more technical solutions provided in the embodiments of this application have at least the following technical effects or advantages: This application embodiment acquires the node status of each information node and constructs cyber-physical coupling constraints for a multi-energy distribution system by combining the initial opening and closing status of the power line. Then, it acquires the operational constraints of the multi-energy distribution system and, combined with the cyber-physical coupling constraints, performs optimal power flow calculations using a first objective function to obtain the retained load of the multi-energy distribution system during the disaster response period. Next, it constructs physical system recovery constraints and information system recovery constraints, and generates a cyber-physical collaborative recovery model by combining the cyber-physical coupling constraints. Based on the cyber-physical collaborative recovery model and the real-time updated operational constraints, it processes the updated load retention value using a second objective function to obtain the retained load of the multi-energy distribution system during each disaster recovery period. Finally, it determines the resilience index based on the retained load during the disaster response period and the retained load at each moment during the disaster recovery period, effectively improving the accuracy of fault impact range prediction and minimizing recovery time to reduce power outage losses. Attached Figure Description

[0026] To more clearly illustrate the technical solutions of the embodiments of this application, the drawings used in the description of the embodiments of this application or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0027] Figure 1 A flowchart illustrating a method for assessing the resilience of a multi-energy power distribution system, provided as an embodiment of this application; Figure 2 A diagram illustrating the operational state requirements of information nodes at both ends of a power line under different initial states during power line opening and closing, as provided in embodiments of this application. Figure 3 An example diagram illustrating the change of active power output by the power supply node over time, provided in an embodiment of this application. Figure 4 The load retention rate variation curve provided in the embodiments of this application; Figure 5 This is a schematic diagram of a multi-energy power distribution system resilience assessment device provided in an embodiment of this application. Detailed Implementation

[0028] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this invention. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without creative effort are within the scope of protection of this invention.

[0029] The following description of some technologies involved in the embodiments of this application is provided to aid understanding and should be considered merely exemplary. Therefore, those skilled in the art should recognize that various changes and modifications can be made to the embodiments described herein without departing from the scope and spirit of this application. Similarly, for clarity and brevity, some descriptions of well-known functions and structures are omitted in the following description.

[0030] Figure 1 This is a flowchart of a method for assessing the resilience of a multi-energy power distribution system according to an embodiment of this application, including steps S10 to S50. Figure 1 This is merely one execution order shown in the embodiments of this application and does not represent the only execution order for a multi-energy power distribution system resilience assessment method. Where the final result can be achieved, Figure 1 The steps shown can be performed in parallel or in reverse order.

[0031] S10: Use the information substation of the multi-energy power distribution system as information nodes, and obtain the working status of each information node as the node status.

[0032] Specifically, the multi-energy power distribution system is abstracted into a node system, and the physical layer uses... This indicates that the information layer uses It indicates. Among them, Let be the set of nodes in a power system, representing the set of all busbars in a multi-energy distribution system; The set of power branches represents the collection of all power lines in a multi-energy distribution system; The information node set represents the collection of all information substations in a multi-energy power distribution system. It is an information branch set, representing the set of all communication links in a multi-energy power distribution system.

[0033] For cyber-physical coupled multi-energy power distribution systems, physical equipment is controlled by information substations. In disaster situations (such as typhoon disasters), the working status of physical equipment is affected by the working status of information substations. It is necessary to analyze the working status of the power communication system, that is, to analyze the status of each information node in the information layer based on the fault scenario.

[0034] Generally, whether an information node can function properly depends on two factors: 1) whether the information node has a valid transmission path connected to the control center; and 2) whether the information node has sufficient power supply.

[0035] Assume that all information nodes are equipped with uninterruptible power supplies to ensure sufficient power supply to the information substations. In this case, the operational status of the information nodes depends on the working condition of the communication link on their transmission path to the control center.

[0036] Assuming information nodes There is a connection with the control center As long as at least one transmission path is valid, the information node can function normally. work status Represented as,

[0037] In the formula, Represents information nodes The corresponding number The validity of the transmission path. , representing information node The corresponding number The transmission path is valid. , representing information node The corresponding number The transmission path is invalid.

[0038] Whether information can be accurately transmitted along a path depends on path topology, communication delay, and communication errors. Therefore, the effectiveness of each transmission path depends on three requirements: topology requirements, delay requirements, and error-free requirements, as detailed below.

[0039] In the formula, Represents information nodes The corresponding number Topology requirements for each transmission path Represents information nodes The corresponding number Transmission path delay requirements Represents information nodes The corresponding number Error-free requirements for each transmission path.

[0040] Among them, topology requirements specifically refer to the requirement that the transmission path can function normally, assuming that the information node... The corresponding number The transmission path contains A communication link, only when The transmission path can only function normally when all communication links are working properly. ,otherwise .

[0041]

[0042] In the formula, Indicates the first The operational status of each communication link. Indicates the first The communication link is working normally. Indicates the first One communication link is malfunctioning.

[0043] The delay requirement specifically refers to ensuring that the delay in path transmission information is within the allowable range, so as to prevent information substations from being unable to receive control commands and adjustments in a timely manner. This indicates that the delay requirement has been met. This indicates that the delay requirement is not met.

[0044]

[0045] In the formula, Represents information nodes The corresponding number The latency value of each transmission path, This indicates the preset latency threshold.

[0046] The path transmission delay in this embodiment includes node delay and line transmission delay. This path transmission delay depends on the number of information nodes traversed and the path length. Specifically, the delay for each information node is... (e.g. 25) The line transmission delay is (e.g. 5) ),but

[0047] In the formula, Represents information nodes The corresponding number The latency value of each transmission path, Represents information nodes The corresponding number The number of information nodes traversed by the transmission path Represents information nodes The corresponding number The path length of each transmission path.

[0048] The error-free requirement specifically refers to the requirement that information transmission be free of communication errors. Since the information system in this embodiment uses synchronous digital layer communication, checks are performed at both ends of the transmission path during information transmission to correct error codes. Therefore, the transmitted information can be considered completely reliable. .

[0049] S20: Obtain the initial opening and closing state of the power line and combine it with the node state of the corresponding node to construct the cyber-physical coupling constraints of the multi-energy power distribution system.

[0050] Specifically, when an information substation fails, it cannot receive control commands from the control center, and the corresponding physical equipment cannot operate in a timely manner, causing coupled failures in the physical equipment, including power supply equipment tripping and remote control switches failing to operate. During the disaster response process, optimal power flow calculations are performed, taking into account the constraints of cyber-physical coupled failures, and load shedding is minimized through network reconfiguration and mutual support between different energy sources.

[0051] The cyber-physical coupling constraints of the multi-energy power distribution system in this embodiment include power output constraints of the power supply nodes and topology control constraints. Specifically, 1) Power supply equipment tripped In cyber-physical systems, the output of a power supply node (considering a power supply device as a power supply node) requires monitoring and control by the corresponding information node.

[0052] The control center optimizes the calculation and outputs control commands to adjust the power output. When the information node is operating normally, the power output of the power supply node can be adjusted between the minimum and maximum power output limits.

[0053] When an information node fails, the output of the power supply node cannot be controlled by the control center. Although small load changes can be handled through local control, the adjustment range is limited.

[0054] When a power outage causes significant disturbance, the system experiences power imbalance, leading to power supply equipment tripping and output dropping to zero. The power supply node output constraint set in this embodiment is as follows:

[0055] In the formula, Indicates the power supply node The corresponding information nodes Indicates the power supply node Corresponding information nodes The node status, This indicates that the information node has failed. This indicates that the information node is functioning normally; Indicates the power supply node Minimum active power output, Indicates the power supply node Minimum active power output, Indicates the power supply node Maximum active power output Indicates the power supply node Output active power, This represents the set of all power supply nodes in a multi-energy power distribution system. It represents the set of information nodes corresponding to all power supply nodes in a multi-energy power distribution system.

[0056] 2) Remote switch fails to operate In a multi-energy distribution system with multiple distributed power sources, there is a situation where the two ends of the line supply power to each other. Assume the power lines... Busbar nodes at both ends and Each location is equipped with a remote control switch to control the opening and closing of the circuit, and the information nodes corresponding to the busbar nodes. and The system receives instructions from the control center to dynamically adjust the opening and closing status of the remote control switch, thereby enabling network reconfiguration.

[0057] When the information node fails, the remote control switch fails to operate, making it impossible to control the opening and closing of the power line.

[0058] The operational requirements of the information nodes at both ends of a power line differ depending on its initial state. Specifically, for a power line initially open, the line can only be closed when both information nodes at both ends are functioning normally; for a power line initially closed, the line can be opened as long as one information node is functioning normally. Figure 2 As shown.

[0059] The control relationship between the initial state of the power line and the working state of the information nodes at both ends of the power line on the opening and closing state of the power line is shown in Table 1.

[0060] In this system, 0 represents the power line being open, and 1 represents the power line being closed. If the information node can function normally, it is in a normal state, represented by 1; if the information node cannot function normally, it is in a failed state, represented by 0.

[0061] Table 1

[0062] Based on the relationships in Table 1, the topology control constraints are constructed as follows:

[0063] In the formula, Endpoints endpoints The initial open / closed state of the power lines between them. Endpoints Corresponding information nodes The node status, Endpoints Corresponding information nodes The node status, This indicates the permissible opening and closing states of power lines during network reconfiguration in a multi-energy distribution system. It represents the set of all power lines in a multi-energy distribution system. This represents the constraint control coefficient. In power lines, Endpoints and endpoints The power lines between them.

[0064] It should be noted that in power lines, the endpoints... and endpoints This represents the nodes at both ends of a power line.

[0065] S30: Obtain the operating constraints of the multi-energy distribution system, and combine them with cyber-physical coupling constraints to perform optimal power flow calculation through the first objective function, thereby obtaining the retained load of the multi-energy distribution system during the disaster response period.

[0066] The disaster response phase refers to the period during which a disaster occurs.

[0067] Specifically, a multi-energy distribution system includes power systems, natural gas systems, heating systems, and energy conversion components. When performing optimal power flow calculations, it is necessary to conform to actual operating laws and meet the basic operating constraints of each energy system.

[0068] For ease of reference, the variable definitions under various constraints in the multi-energy power distribution system in this embodiment are shown in Table 2.

[0069] Table 2

[0070] I. Power System 1) Power line state constraints Under the damage of disasters (such as typhoons), some power lines may be out of service due to faults and remain disconnected until they are repaired. Because some information nodes fail, corresponding remote switches may malfunction, causing some power lines that are not faulty to be unable to open or close according to optimal power flow calculations, thus losing control. Power lines initially disconnected (such as tie lines) cannot be closed and remain disconnected; power lines initially closed cannot be disconnected and remain closed. The constraints are as follows:

[0071] In the formula, This indicates the permissible opening and closing states of power lines during network reconfiguration in a multi-energy distribution system. This indicates that the power line is disconnected. This indicates that the power line is in a closed state. This refers to a set of power lines that remain disconnected due to faults or loss of control. This refers to a set of power lines that remain closed due to loss of control.

[0072] 2) Radial topological constraints When reconfiguring a multi-energy power distribution system, it is necessary to prevent the network from forming a ring topology, thereby avoiding circulating currents, increasing line losses, and maintaining the network in a radial topology to facilitate voltage control and operation management.

[0073] When some lines fail and disconnect, the system reconstructs islands through network reconfiguration. Islanding islands containing controllable power sources powering loads are considered operational islands, ensuring each operational island has only one power source. Islanding islands without power supply paths to loads are considered outage islands, where any node within an island is designated as a virtual power source to meet unified constraints. The number of controllable power sources represents the number of operational islands, and the number of virtual power sources represents the number of outage islands. The total number of closed power lines must equal the total number of nodes minus the number of operational and outage islands.

[0074] Specify the upstream and downstream relationships of nodes. For each power line, the positive direction of power flow is from the node with the smaller number to the node with the larger number. The node with the smaller number is designated as the upstream node, and the node with the larger number is designated as the downstream node. Represents a node The set of downstream nodes, Represents a node The set of upstream nodes, nodes The power line connected to the upstream node is called a node. The upstream branch, and the power line connecting the downstream node, are the nodes. The downstream branches of the load node. Using the single-commodity flow method, virtual power is defined to flow in the distribution network topology. A requirement is that the difference between the sum of the virtual power of the downstream branches and the sum of the virtual power of the upstream branches of a load node is -1, meaning the node meets the unit virtual load demand, indicating the presence of a power source or virtual power source supplying power to that node. The virtual power emitted by controllable power sources and virtual power sources can be any value. The magnitude of the virtual power of a power line is related to its open / closed state. When the power line is closed, the virtual power flowing through is unrestricted; when the power line is open, the virtual power flowing through is zero. The open / closed state of the power line is constrained by the node's virtual power demand and the relationship between the power line's virtual power and its open / closed state. Specific constraints are as follows.

[0075]

[0076]

[0077]

[0078] In the formula, Indicates a set of power branch circuits; Represents the set of nodes in a power system; Represents the set of controllable power supply points; Represents a set of virtual power points; It is a large number.

[0079] 3) Voltage constraints at both ends of the power line For a closed power line, the power flow through the line must satisfy the voltage drop formula, while the voltage difference between the two ends of an open power line is unrestricted, and the constraint is expressed as follows.

[0080] In the formula, Indicates power lines The resistance; Indicates power lines The reactance; This indicates the system's rated voltage.

[0081] 4) Power line capacity constraints The power flow magnitude of a closed power line is limited by the line's capacity. When the line is open, both active and reactive power flow is zero. The constraints are expressed as follows:

[0082] In the formula, Indicates power lines The maximum apparent power.

[0083] 5) Node power balance constraints The power flowing into the node is equal to the power flowing out of the node, under the following constraint:

[0084] 6) Node voltage constraints The voltage at normal load points is between the system's minimum and maximum allowable voltages. The voltage at controllable power supply points is equal to the system's rated voltage. The voltage at virtual power supply points is 0. The voltage at all load points within the power outage island is 0. The constraints are as follows:

[0085] In the formula: Indicates the minimum allowable voltage of the system; This indicates the maximum allowable voltage of the system.

[0086] 7) Node supply and demand constraints The output power constraints of the controllable power supply point are as follows.

[0087] In the formula, , Indicates controllable power point The minimum and maximum active power output; Indicates controllable power point Maximum apparent power output; Controllable power supply point The minimum operating power factor.

[0088] After the response, the active and reactive power of the load are retained proportionally within a range less than the pre-fault load. The load retention rate within the power outage island is 0, subject to the following constraints.

[0089] In the formula, Indicates the load point Normal active power value; Indicates the load point The normal value for reactive power.

[0090] (2) Natural gas system 1) Nodal pressure and supply and demand constraints The constraints on natural gas node pressure and supply and demand are as follows.

[0091]

[0092]

[0093] In the formula, Represents a node The square of the lower limit of air pressure Represents a node The square of the upper limit of air pressure; Represents the set of natural gas nodes; This represents the air pressure value at a fixed pressure node; Represents the set of fixed pressure nodes; Represents a node Minimum supply Represents a node The maximum supply; Represents the set of gas source nodes; Represents the set of natural gas load nodes; Represents a node Minimum demand Represents a node The maximum required value.

[0094] The demand value after the response is retained within a range that is less than the pre-fault load, subject to the following constraints:

[0095] In the formula, Indicates the load point Normal demand value Indicates the load point after response The demand value, Indicates the load point Retention rate.

[0096] 2) Nodal airflow balance constraints When the inflow airflow at a node equals the outflow airflow, it is represented as follows:

[0097] In the formula, Indicates pipeline airflow size The amount of airflow through a pipe is limited by the pipe's capacity, as shown below.

[0098] In the formula: , Indicates pipeline Minimum and maximum airflow values; This refers to a collection of natural gas pipelines.

[0099] 3) Pipeline gas flow constraints The flow of natural gas through pipelines causes pressure losses, which satisfy the pipeline gas flow equation. A piecewise linearization method is used to handle this.

[0100]

[0101]

[0102]

[0103]

[0104]

[0105] In the formula, Indicates pipeline The Weymouth constant; Indicates pipeline The first of the airflow magnitude Segmentation point; Indicates the total number of segments; and Indicates pipeline The right-hand side of the airflow equation and its first term Segmentation point; Indicates pipeline gas flow solution in the first The proportion of values ​​within a segment; It is a 0-1 variable, which controls the subsequent segment to take a value only when the preceding segment is full.

[0106] 4) Compressor operating constraints The airflow through the compressor is limited by the pipeline capacity. For a working compressor, the gas pressure increases after the natural gas flows through it. Assuming the pressure increase is constant and that the compressor does not consume natural gas and the gas flows in one direction, it can be represented as follows.

[0107]

[0108] In the formula, This represents the set of power lines for the compressor; This represents the increase in air pressure squared.

[0109] (3) Heating system The heating system consists of a heating network and a regenerative network, which have the same topology but opposite hot water flow directions. Therefore, the analysis only requires analyzing the heating network when considering network topology. The heating system employs quality regulation. The operating model of the thermal system includes a hydraulic model and a thermal model. It is assumed that the water flow rate and heat load of each pipe are fixed, the hydraulic model yields a feasible solution, and the heat load can be adjusted.

[0110] 1) Pipeline Flow-Mass Balance Model Using a quality-controlled approach, the flow rates of the thermal pipelines can be calculated given the fluid-mass requirements of each heat load. When no pipelines in the thermal network are out of service, the calculation is performed based on the initial flow rates of the pipelines, with the following constraints:

[0111] In the formula: Represents a node The magnitude of the heat load fluid is assumed to be constant. Indicates pipeline The flow rate; This represents the set of thermal nodes.

[0112] 2) Nodal heat load model and temperature constraints Water flowing into a heat load will lower its temperature, as shown below.

[0113] In the formula, This indicates the specific heat capacity of water.

[0114] The required heat power after the response load is retained within a range less than the pre-fault load, subject to the following constraints:

[0115] In the formula, This represents the normal value of the required heat power at load point j.

[0116] The temperature limits for each node are as follows:

[0117] In the formula: These represent the nodes of the heating network. The lowest and highest temperatures; These represent the minimum and maximum temperatures of node j in the regenerating network, respectively.

[0118] 3) Pipeline heat transfer model The heat transfer of hot water in a pipe will incur losses and delays. Ignoring the heat transfer time, the temperatures at the beginning and end of the pipe must satisfy the following constraints.

[0119]

[0120] In the formula, This represents the ambient temperature, assumed to be 0℃. For pipelines Heat transfer coefficient per unit length; This refers to the length of the pipe. This refers to the collection of pipes in a heating system.

[0121] The temperature limits at the beginning and end of the pipeline are as follows.

[0122] In the formula: These represent the pipes in the heating pipeline and the regeneration pipeline, respectively. Minimum and maximum water temperature limits at the beginning and end.

[0123] 4) Nodal thermal mixing model When hot water mixes at a node, assuming no heat loss, the total heat of the inflow fluid equals the total heat of the outflow fluid. After the hot water in the heating network and the regenerating network mixes at the confluence, the initial temperature of subsequent pipes is equal to the mixed temperature. The constraints are as follows:

[0124]

[0125]

[0126]

[0127]

[0128] In the formula, This represents the set of heat source nodes.

[0129] (4) Energy conversion element

[0130] In the formula, Indicates the energy converted; This represents the conversion coefficient, which includes losses such as efficiency and calorific value. Indicates the energy that has been converted; Represents a y-type conversion element The Middle One component; This represents the overall conversion element.

[0131] Combining the operational constraints and cyber-physical coupling constraints of the multi-energy distribution system, the optimal power flow calculation aims to minimize load outages. Since electrical, gas, and heat loads are not equivalent, loads of different energy types are classified as electrical loads. The first objective function is constructed to minimize the load cut off by the system, thus obtaining the maximum load remaining in the multi-energy distribution system during the disaster response period. The first objective function is specifically as follows:

[0132] In the formula, This indicates the retained load of the multi-energy power distribution system during each disaster response period; These represent the weighting coefficients for electricity, gas, and heat loads, respectively. and These represent the coefficients for the reduction of gas load and heat load to electrical load, respectively. These represent the retained values ​​of electricity, gas, and heat loads during the disaster response period, respectively; among them, Indicates electrical load, Indicates gas load, Indicates heat load, Represents the first in the power system Each node Represents the first in the natural gas system Each node Represents the first in the thermodynamic system 1 node Represents the set of nodes in a power system. Represents the set of nodes in a natural gas system. This represents the set of nodes in a heating system.

[0133] S40: Construct physical system recovery constraints and information system recovery constraints, and combine them with information-physical coupling constraints to generate an information-physical collaborative recovery model.

[0134] Specifically, due to the limited number of maintenance teams, it is impossible to repair all faulty lines simultaneously under extreme circumstances. Therefore, this embodiment considers resource constraints based on the actual situation, mainly including physical system recovery constraints and information system recovery constraints.

[0135] Furthermore, the physical system recovery constraints include: After a disaster (such as after a typhoon), damaged power lines and communication links need to be repaired manually. Because the repairs rely on manual labor, the recovery status is affected by the available resources. Assuming that the time required to repair each power line is considered a time period, and the number of power lines that can be restored in each time period cannot exceed the available resources, a first constraint is established.

[0136] The expression for the first constraint is as follows:

[0137] Since the recovery status of a power line is affected by the recovery status of adjacent power lines, if both the upstream and downstream power lines of a power line are in a fault state, the power line cannot be restored. Therefore, a second constraint is constructed.

[0138] The expression for the second constraint is as follows:

[0139] Since power line restoration is a prerequisite for power lines to participate in opening and closing regulation during the reconfiguration of a multi-energy distribution system, and if the power line is not restored, it remains disconnected, a third constraint is constructed.

[0140] The expression for the third constraint is as follows:

[0141] In the formula, Endpoints With endpoints The power lines between The state of recovery at any time, Endpoints With endpoints The power lines between The state of recovery at any time, Representation and endpoints , The adjacent upstream or downstream power lines between the power lines The state of recovery at any time, Indicates the restored status of power lines. and Indicates the endpoints on an upstream or downstream power line. Endpoints With endpoints The power lines between The opening and closing state at any moment, It represents the set of all power lines in a multi-energy distribution system. Representation and endpoints , The set of all power lines adjacent to each other. This indicates the number of power line maintenance teams. Indicates the total number of disaster recovery periods; in power lines, Endpoints and endpoints The power lines between them.

[0142] Furthermore, load restoration constraints include: Since some loads are isolated after a disaster (such as a typhoon), load restoration is achieved by restoring the path between the de-energized load and the controllable power source. Therefore, the load restoration status is related to the restoration status of the power lines adjacent to the load. If all the power lines adjacent to the load are faulty, the load cannot be connected to the controllable power source and cannot be restored. Therefore, a fourth constraint is constructed.

[0143] The expression for the fourth constraint is as follows:

[0144] The recovery status of a load is determined by whether a valid path exists between the load and a controllable power source at a given moment. This involves iterating through all paths from the load to all controllable power sources. If all paths are broken, it indicates that there is no power supply path to the load, meaning the load has not recovered. Therefore, a fifth constraint is constructed.

[0145] The expression for the fifth constraint is as follows:

[0146]

[0147] In the formula, Indicates load node exist The recovery state at any moment, u D Indicates the restored status of power lines. Indicates load node The set of adjacent power lines, Endpoints With endpoints The power lines between The state of recovery at any time, Indicates the restored status of power lines. Indicates load node exist Path to all controllable power sources at all times The opening and closing state, Indicates the path from the load node to all controllable power sources. The first A power line in The opening and closing state at any moment, Represents the set of all load nodes. This represents the set of all power lines along the path. This represents the set of all paths from a load node to a controllable power source; in power lines, Endpoints and endpoints The power lines between them.

[0148] It should be noted that in power lines, the endpoints... and endpoints This represents the nodes at both ends of a power line.

[0149] Furthermore, information system recovery constraints include: Since fiber optic composite overhead ground wires are mostly constructed along with the power lines, the corresponding power lines need to be restored before the communication links can be restored. In other words, the communication links cannot be restored if the power lines are not restored. Therefore, a sixth constraint is constructed.

[0150] The expression for the sixth constraint is as follows:

[0151] Since communication links rely on manual repair, the number of communication links repaired at the same time cannot exceed the recovery resources. Therefore, a seventh constraint is constructed.

[0152] The expression for the seventh constraint is as follows:

[0153] In the formula, Endpoints With endpoints The recovery status of the communication link between them. Endpoints With endpoints The power lines between The state of recovery at any time, Endpoints With endpoints The communication link between them The state of recovery at any time, This represents the set of all communication links in a multi-energy power distribution system. This indicates the number of communication link maintenance teams. Indicates the restored status of power lines. Indicates the recovery status of the communication link. Indicates the total number of disaster recovery periods; in communication links, Endpoints and endpoints Communication links between them.

[0154] It should be noted that in a communication link, the endpoint and endpoints This refers to the nodes at both ends of the communication link.

[0155] Information node working status recovery refers to restoring a valid path between the information node and the control center. Therefore, priority is given to restoring the communication link between normal information nodes and failed information nodes. If the information nodes at both ends of a communication link are in a failed state, the communication link cannot be restored. Therefore, the eighth constraint is constructed.

[0156] The expression for the eighth constraint is as follows:

[0157] Similarly, when all the communication links adjacent to an information node are in a faulty state, the information node does not recover. Therefore, a ninth constraint is constructed.

[0158] The expression for the ninth constraint is as follows:

[0159] In the formula, Endpoints With endpoints The recovery status of the communication link between them. Represents the endpoint at one end of the communication link. The corresponding information nodes are in , Indicates the endpoint at the other end of the communication link. The corresponding information nodes are in The state of recovery at any time, Endpoints The corresponding information nodes are in The state of recovery at any time, Indicates the recovery status of the communication link. Indicates the recovery status of the information node. This represents the set of all communication links in a multi-energy power distribution system. Represents the set of all information nodes. Indicates the total number of disaster recovery periods; in communication links, Endpoints and endpoints Communication links between them.

[0160] It should be noted that in a communication link, the endpoint and endpoints This refers to the nodes at both ends of the communication link.

[0161] Furthermore, during the recovery phase, system network reconstruction and the adjustment of different energy outputs are still affected by cyber-physical coupling. Power supply nodes with normally functioning information nodes can adjust their output according to commands from the control center; if an information node fails, output adjustment must be done manually. The communication time and manual adjustment time slow down the equipment response speed, reducing the output range. In other words, the uncontrollability of the power supply equipment in the previous phase leads to a reduction in the output range in this phase. Therefore, a tenth constraint is constructed. Figure 3 As shown, These represent the time intervals of the power supply node when communication is normal. The lower and upper limits of output.

[0162] The expression for the tenth constraint is as follows:

[0163]

[0164] The allowed opening and closing state of the line is related to the line opening and closing state at the previous moment and the working state of the information nodes at both ends of the line at this moment. As the failed information node is restored, the allowed opening and closing state of the line is updated at different times. Therefore, the eleventh constraint is constructed.

[0165] The expression for the eleventh constraint is as follows:

[0166] In the formula, Indicates the power supply node The active power output at all times. Indicates the power supply node exist The active power output at all times. Indicates active power. Indicates the power supply node The corresponding active power regulation speed of the equipment Indicates the duration of the disaster recovery period. Indicates the power supply node The corresponding information nodes are in The state of recovery at any time, This indicates the time required for manual adjustment of output. Indicates the power supply node Minimum active power output, Indicates the power supply node Maximum active power output Endpoints With endpoints The power lines between The opening and closing state at any moment, Represents the endpoint at one end of the communication link. The corresponding information nodes are in The state of recovery at any time, Indicates the endpoint at the other end of the communication link. The corresponding information nodes are in The state of recovery at any time, Endpoints With endpoints The power lines between The opening and closing state at any moment, Represents the set of power supply nodes. Indicates the total number of disaster recovery periods. This represents the set of endpoints of all communication links in a multi-energy distribution system. This represents the constraint control coefficient. In the communication link, Endpoints and endpoints Communication links between them.

[0167] S50: Based on the cyber-physical collaborative recovery model and real-time updated operational constraints, the retained load of the multi-energy distribution system at each moment during the disaster recovery period is obtained through the second objective function.

[0168] The second objective function is as follows:

[0169] In the formula, This indicates the retained load of the multi-energy power distribution system during each disaster recovery period; These represent the weighting coefficients for electricity, gas, and heat loads, respectively. and These represent the coefficients for the reduction of gas load and heat load to electrical load, respectively. Let represent the retained values ​​of electricity, gas, and heat loads at time t during the disaster recovery period; where, Indicates electrical load, Indicates gas load, Indicates heat load, Represents the first in the power system Each node Represents the first in the natural gas system Each node Represents the first in the thermodynamic system 1 node Represents the set of nodes in a power system. Represents the set of nodes in a natural gas system. This represents the set of nodes in a heating system.

[0170] S60: Determine resilience indicators based on the retained load during the disaster response period and the retained load at each time point during the disaster recovery period.

[0171] Specifically, at each stage of a disaster, resilience is mainly manifested in its ability to defend against the disaster and its ability to recover after the disaster; that is, the ability to minimize load loss during a disaster and to quickly restore normal operation after the disaster. During a disaster, load retention rate is chosen to represent defense capability, and after a disaster, load recovery time is chosen to represent recovery capability. The load retention rate change curve is used to reflect the system's resilience, such as... Figure 4 As shown, This represents the missing area corresponding to the closed region constructed by the load retention rate curve and the constant 100%. The area of ​​the fault curve during the disaster response period.

[0172] The resilience index is determined based on the fault curve area and the load-loss area using the resilience index calculation formula. The resilience index calculation formula is:

[0173] In the formula, Indicators of resilience Represents a set of fault scenarios. Indicates the fault scenario The probability of occurrence Indicates the area of ​​load deficiency. This represents the area of ​​the fault curve.

[0174] In the above formula for calculating the resilience index, the area of ​​the fault curve is used. As a molecule, the larger the molecule, the greater the system load retention rate during a disaster. (The missing area is then considered.) Area of ​​fault curve The smaller the denominator, the shorter the post-disaster load recovery time, corresponding to the resilience index. The larger the value, the better. Since multiple fault scenarios exist when using simulation methods to assess resilience, the resilience index in this embodiment takes the expected value under multiple scenarios.

[0175] Furthermore, to compare the impact of cyber-physical coupling on system resilience, four case studies were set up for comparative analysis. Case 1 did not consider the cyber-physical coupling relationship and assumed that the information system operated normally during the fault. Case 2 only considered whether the remote control switch was controllable. Case 3 only considered the operating status of the information node of the power supply node. Case 4 considered both the controllability of the remote control switch and the operating status of the information node of the power supply node. Scenarios where the information system's operating status was affected by the fault were selected, and the equivalent load deficit curves for different cases were obtained as follows: Figure 4 As shown in Table 3, the parameters under different cases are compared.

[0176] Table 3

[0177] Case 1 did not consider information system failures, and the simulated fault scale was significantly smaller than that of Case 4, resulting in a higher resilience index than other cases. Due to information system failures, the resilience index of Case 4 was 1.96% lower than that of Case 1. This demonstrates that ignoring the operational status of the information layer and only modeling and evaluating the physical layer leads to optimistic estimates of the system's load deficit during a disaster and post-disaster recovery time, significantly deviating from the evaluation results under cyber-physical coupling conditions. In Case 3, the information failure of the power supply node limited output adjustment. With other normally operating power supply equipment and controllable remote switches, network reconfiguration could optimize the allocation of limited resources, transmitting power to critical load areas and minimizing the deficit. In Case 2, the remote switch malfunctioned, preventing power lines from being opened and closed according to optimal power flow calculations for network reconfiguration. Even with adjustable power supply node output, power could not be effectively transmitted to demand areas; therefore, the load deficit in Case 3 was smaller than in Case 2, and the resilience index was greater. The curves for Case 3 and Case 2 are located between Case 1 and Case 4, closer to Case 4, indicating that any factor in cyber-physical coupling will have a significant impact on the system resilience assessment results. Improving the resilience of cyber-physical coupling systems requires attention to the synergy between power supply and output regulation and network reconfiguration capabilities.

[0178] In this embodiment, by establishing a working state model of the power communication system, considering the cyber-physical coupling fault mechanism, and modeling it as a cyber-physical coupling constraint to perform optimal power flow calculation during the disaster response period, a cyber-physical collaborative recovery strategy is proposed in the post-disaster recovery stage to achieve minimum load shedding and rapid load recovery, thereby improving the accuracy of fault impact range prediction and minimizing recovery time to reduce power outage losses.

[0179] While this application provides method operation steps as shown in the embodiments or flowcharts, more or fewer operation steps may be included based on conventional or non-inventive labor. The order of steps listed in this embodiment is merely one possible execution order among many and does not represent the only execution order. In actual device or client product execution, the method can be executed sequentially according to this embodiment or the accompanying drawings, or in parallel (e.g., in a parallel processor or multi-threaded processing environment).

[0180] like Figure 5 As shown in the illustration, this application also provides a device for assessing the resilience of a multi-energy power distribution system. The device includes: The node status acquisition module 10 is used to take the information substation of the multi-energy power distribution system as information nodes and acquire the working status of each information node as the node status.

[0181] The cyber-physical coupling constraint construction module 20 is used to obtain the initial opening and closing state of the power line and combine it with the node state of the corresponding information node to construct the cyber-physical coupling constraints of the multi-energy power distribution system.

[0182] The optimal power flow calculation module 30 is used to obtain the operating constraints of the multi-energy distribution system and, in combination with the cyber-physical coupling constraints, perform optimal power flow calculation through the first objective function to obtain the retained load of the multi-energy distribution system during the disaster response period.

[0183] The collaborative recovery model generation module 40 is used to construct physical system recovery constraints and information system recovery constraints, and to generate a cyber-physical collaborative recovery model by combining cyber-physical coupling constraints.

[0184] The module 50 for determining the retained load during the recovery period is used to obtain the retained load of the multi-energy distribution system at each moment during the disaster recovery period based on the cyber-physical collaborative recovery model and real-time updated operational constraints, through a second objective function.

[0185] The resilience index calculation module 60 is used to determine the resilience index based on the retained load during the disaster response period and the retained load at each time during the disaster recovery period.

[0186] Some modules in the apparatus described in this application can be described in the general context of computer-executable instructions that are executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, classes, etc., that perform a specific task or implement a specific abstract data type. This application can also be practiced in distributed computing environments where tasks are performed by remote processing devices connected via a communication network. In distributed computing environments, program modules can reside in local and remote computer storage media, including storage devices.

[0187] The apparatus or module described in the above embodiments can be implemented by a computer chip or physical entity, or by a product with a certain function. For ease of description, the above apparatus is described by dividing it into various modules according to their functions. When implementing the embodiments of this application, the functions of each module can be implemented in one or more software and / or hardware. Of course, a module that implements a certain function can also be implemented by combining multiple sub-modules or sub-units.

[0188] The methods, apparatus, or modules described in this application can be implemented in a computer-readable program code manner. The controller can be implemented in any suitable manner, such as a microprocessor or processor and a computer-readable medium storing computer-readable program code (e.g., software or firmware) executable by the (micro)processor, logic gates, switches, application-specific integrated circuits (ASICs), programmable logic controllers, and embedded microcontrollers. Examples of controllers include, but are not limited to, the following microcontrollers: ARC 625D, Atmel AT91SAM, Microchip PIC18F26K20, and Silicon Labs C8051F320. A memory controller can also be implemented as part of the control logic of a memory. Those skilled in the art will also recognize that, in addition to implementing the controller in purely computer-readable program code manner, the same functionality can be achieved by logically programming the method steps to make the controller take the form of logic gates, switches, application-specific integrated circuits, programmable logic controllers, and embedded microcontrollers. Therefore, such a controller can be considered a hardware component, and the means included within it for implementing various functions can also be considered as structures within the hardware component. Alternatively, the device used to implement various functions can be viewed as either a software module that implements the method or a structure within a hardware component.

[0189] This application also provides an apparatus, the apparatus comprising: a processor; a memory for storing processor-executable instructions; wherein, when the processor executes the executable instructions, it implements the method described in this application.

[0190] This application also provides a non-volatile computer-readable storage medium storing a computer program or instructions thereon, which, when executed, enables the method described in this application embodiment to be implemented.

[0191] Furthermore, in the various embodiments of the present invention, each functional module can be integrated into a processing module, or each module can exist independently, or two or more modules can be integrated into a single module.

[0192] The aforementioned storage media include, but are not limited to, Random Access Memory (RAM), Read-Only Memory (ROM), Cache, Hard Disk Drive (HDD), or Memory Card. The memory can be used to store computer program instructions.

[0193] As can be seen from the above description of the embodiments, those skilled in the art can clearly understand that this application can be implemented by means of software plus necessary hardware. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, can be embodied in the form of a software product, or it can be embodied in the process of data migration. The computer software product can be stored in a storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, mobile terminal, server, or network device, etc.) to execute the methods described in various embodiments or some parts of the embodiments of this application.

[0194] The various embodiments described in this specification are presented in a progressive manner. Similar or identical parts between embodiments can be referred to interchangeably. Each embodiment focuses on its differences from other embodiments. All or part of this application can be used in numerous general-purpose or special-purpose computer system environments or configurations. Examples include: personal computers, server computers, handheld or portable devices, tablet devices, mobile communication terminals, multiprocessor systems, microprocessor-based systems, programmable electronic devices, network PCs, minicomputers, mainframe computers, and distributed computing environments including any of the above systems or devices, etc.

[0195] The above embodiments are only used to illustrate the technical solutions of this application, and are not intended to limit this application. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features therein. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of this application.

Claims

1. A method for assessing the resilience of a multi-energy power distribution system, characterized in that, include: The information substations of the multi-energy power distribution system are used as information nodes, and the working status of each information node is obtained as the node status. By obtaining the initial opening and closing state of the power line and combining it with the node state of the corresponding information node, the cyber-physical coupling constraints of the multi-energy power distribution system are constructed. The operational constraints of the multi-energy distribution system are obtained, and combined with the cyber-physical coupling constraints, the optimal power flow is calculated through the first objective function to obtain the retained load of the multi-energy distribution system during the disaster response period. Construct physical system recovery constraints and information system recovery constraints, and combine the aforementioned cyber-physical coupling constraints to generate a cyber-physical collaborative recovery model; Based on the aforementioned cyber-physical collaborative recovery model and real-time updated operational constraints, the retained load of the multi-energy distribution system at each moment during the disaster recovery period is obtained through the second objective function; The resilience index is determined based on the retained load during the disaster response period and the retained load at each moment during the disaster recovery period.

2. The method for assessing the resilience of a multi-energy power distribution system according to claim 1, characterized in that, The cyber-physical coupling constraints of the multi-energy power distribution system include power supply node output constraints and topology control constraints. The power output constraint of the power supply node is as follows: In the formula, Indicates the power supply node The corresponding information nodes Indicates the power supply node Corresponding information nodes The node status, Indicates the power supply node Minimum active power output, Indicates the power supply node Maximum active power output Indicates the power supply node Output active power, This represents the set of all power supply nodes in a multi-energy power distribution system. This represents the set of information nodes corresponding to all power supply nodes in a multi-energy power distribution system. The topology control constraints are: In the formula, Endpoints and endpoints The initial open / closed state of the power lines between them. Endpoints Corresponding information nodes The node status, Endpoints Corresponding information nodes The node status, This indicates the permissible opening and closing states of power lines during network reconfiguration in a multi-energy distribution system. It represents the set of all power lines in a multi-energy distribution system. Represents the constraint control coefficient. In power lines, Endpoints and endpoints The power lines between them.

3. The method for assessing the resilience of a multi-energy power distribution system according to claim 1, characterized in that, The physical system recovery constraints include power line recovery constraints and load recovery constraints; wherein... The power line restoration constraints include: In the formula, Endpoints With endpoints The power lines between The state of recovery at any time, Endpoints With endpoints The power lines between The state of recovery at any time, Representation and endpoints , The adjacent upstream or downstream power lines between the power lines The state of recovery at any time, Indicates the restored status of the power line. and Indicates the endpoints on an upstream or downstream power line. Endpoints With endpoints The power lines between The opening and closing state at any moment, It represents the set of all power lines in a multi-energy distribution system. Representation and endpoints , The set of all power lines adjacent to each other. This indicates the number of power line maintenance teams. Indicates the total number of disaster recovery periods; in power lines, Endpoints and endpoints The power lines between them; The load recovery constraints include: In the formula, Indicates load node exist The recovery state at any moment, u D Indicates the restored status of the power line. Indicates load node The set of adjacent power lines, Endpoints With endpoints The power lines between The state of recovery at any time, Indicates the restored status of the power line. Indicates load node exist Paths to all controllable power sources at all times The opening and closing state, Indicates the path from the load node to all controllable power sources. The first A power line in The opening and closing state at any moment, Represents the set of all load nodes. This represents the set of all power lines along the path. This represents the set of all paths from a load node to a controllable power source; in power lines, Endpoints and endpoints The power lines between them.

4. The method for assessing the resilience of a multi-energy power distribution system according to claim 1, characterized in that, The information system recovery constraints include communication link repair constraints, information node operational state recovery constraints, and information physical coupling constraints; among which... The communication link repair constraints include: In the formula, Endpoints With endpoints The recovery status of the communication link between them. Endpoints With endpoints The power lines between The state of recovery at any time, Endpoints With endpoints The communication link between them The state of recovery at any time, This represents the set of all communication links in a multi-energy power distribution system. This indicates the number of communication link maintenance teams. Indicates the restored status of the power line. Indicates the recovery status of the communication link. Indicates the total number of disaster recovery periods; in communication links, Endpoints and endpoints Communication links between them; The constraints for restoring the working state of the information node include: In the formula, Endpoints With endpoints The recovery status of the communication link between them. Represents the endpoint at one end of the communication link. The corresponding information nodes are in The state of recovery at any time, Indicates the endpoint at the other end of the communication link. The corresponding information nodes are in The state of recovery at any time, Endpoints The corresponding information nodes are in The state of recovery at any time, Indicates the recovery status of the communication link. Indicates the recovery status of the information node. This represents the set of all communication links in a multi-energy power distribution system. Represents the set of all information nodes. Indicates the total number of disaster recovery periods; in communication links, Endpoints and endpoints Communication links between them; The cyber-physical coupling constraints include: In the formula, Indicates the power supply node exist The active power output at all times. Indicates the power supply node exist The active power output at all times. Indicates active power. Indicates the power supply node The corresponding active power regulation speed of the equipment Indicates the duration of the disaster recovery period. Indicates the power supply node The corresponding information nodes are in The state of recovery at any time, This indicates the time required for manual adjustment of output. Indicates the power supply node Minimum active power output, Indicates the power supply node Maximum active power output Endpoints With endpoints The power lines between The opening and closing state at any moment, Represents the endpoint at one end of the communication link. The corresponding information nodes are in The state of recovery at any time, Indicates the endpoint at the other end of the communication link. The corresponding information nodes are in The state of recovery at any time, Endpoints With endpoints The power lines between The opening and closing state at any moment, Represents the set of power supply nodes. Indicates the total number of disaster recovery periods. This represents the set of endpoints of all communication links in a multi-energy distribution system. Represents the constraint control coefficient. In the communication link, Endpoints and endpoints Communication links between them.

5. The method for assessing the resilience of a multi-energy power distribution system according to claim 1, characterized in that, The expression for the first objective function is as follows: In the formula, This indicates the retained load of the multi-energy power distribution system during each disaster response period; These represent the weighting coefficients for electricity, gas, and heat loads, respectively. and These represent the coefficients for the reduction of gas load and heat load to electrical load, respectively. These represent the retained values ​​of electricity, gas, and heat loads during the disaster response period, respectively; among them, Indicates electrical load, Indicates gas load, Indicates heat load, Represents the first in the power system Each node Represents the first in the natural gas system Each node Represents the first in the thermodynamic system 1 node Represents the set of nodes in a power system. Represents the set of nodes in a natural gas system. This represents the set of nodes in a heating system.

6. The method for assessing the resilience of a multi-energy power distribution system according to claim 1, characterized in that, The expression for the second objective function is as follows: In the formula, This indicates the retained load of the multi-energy power distribution system during each disaster recovery period; These represent the weighting coefficients for electricity, gas, and heat loads, respectively. and These represent the coefficients for the reduction of gas load and heat load to electrical load, respectively. Let represent the retained values ​​of electricity, gas, and heat loads at time t during the disaster recovery period; where, Indicates electrical load, Indicates gas load, Indicates heat load, Represents the first in the power system Each node Represents the first in the natural gas system Each node Represents the first in the thermodynamic system 1 node Represents the set of nodes in a power system. Represents the set of nodes in a natural gas system. This represents the set of nodes in a heating system.

7. The method for assessing the resilience of a multi-energy power distribution system according to claim 1, characterized in that, The determination of resilience indicators based on the retention load during the disaster response period and the retention load at each moment during the disaster recovery period includes: Based on the retained load during the disaster response period and the retained load at each moment during the disaster recovery period, a load retention rate curve is plotted, and the load loss area is determined based on the closed area constructed by the load retention rate curve and the constant 100%. The area of ​​the fault curve is obtained and combined with the area of ​​load loss, and the resilience index is determined by the resilience index calculation formula. The formula for calculating the resilience index is: In the formula, Indicators of resilience Represents a set of fault scenarios. Indicates the fault scenario The probability of occurrence Indicates the area of ​​load deficiency. This represents the area of ​​the fault curve.

8. A device for assessing the resilience of a multi-energy power distribution system, characterized in that, include: The node status acquisition module is used to take the information substation of the multi-energy power distribution system as information nodes and acquire the working status of each information node as the node status. The cyber-physical coupling constraint construction module is used to obtain the initial opening and closing state of the power line and combine it with the node state of the corresponding information node to construct the cyber-physical coupling constraints of the multi-energy power distribution system. The optimal power flow calculation module is used to obtain the operating constraints of the multi-energy distribution system, and in combination with the information-physical coupling constraints, to perform optimal power flow calculation through the first objective function to obtain the retained load of the multi-energy distribution system during the disaster response period; The collaborative recovery model generation module is used to construct physical system recovery constraints and information system recovery constraints, and to generate a cyber-physical collaborative recovery model by combining the cyber-physical coupling constraints. The module for determining the retained load during the recovery period is used to obtain the retained load of the multi-energy distribution system at each moment during the disaster recovery period based on the cyber-physical collaborative recovery model and real-time updated operating constraints, through a second objective function. The resilience index calculation module is used to determine the resilience index based on the retained load during the disaster response period and the retained load at each time during the disaster recovery period.

9. An apparatus for performing a method for assessing the resilience of a multi-energy power distribution system, characterized in that, include: processor; Memory used to store processor-executable instructions; When the processor executes the executable instructions, it implements the method as described in any one of claims 1 to 7.

10. A non-volatile computer-readable storage medium, characterized in that, Includes storage of computer programs or instructions that, when executed, cause the method as described in any one of claims 1 to 7 to be implemented.

Images