Energy system capacity optimization and mobile resource space-time decoupling cooperative scheduling method, system and device under extreme disaster disturbance and medium

Abstract

The invention relates to the technical field of energy scheduling, in particular to an energy system capacity optimization and mobile resource space-time decoupling cooperative scheduling method, system and device under extreme disaster disturbance and a medium. Load change parameters are coupled to form a disaster-causing response database; renewable energy output fluctuation and user load demand deviation are converted into fuzzy variables, and a capacity configuration optimization model is established in combination with a disaster-causing response database to generate multi-energy complementary capacity configuration; mobile resources are introduced to establish an elastic lifting scheduling model, and an elastic lifting scheme is generated by implementing damaged road network topology reconstruction of the mobile resources, space-time transfer chain modeling and action state coordinated regulation and control. According to the method, collaborative decision-making of capacity optimization of equipment and dynamic scheduling of mobile resources is driven through the disaster-causing response database, and dual elastic gains of fixed resource optimization and after-disaster quick response of the mobile resources under extreme disaster events are achieved.

Description

Technical Field

[0001] The present invention relates to the technical field of energy dispatch, and particularly to a method, system, device and medium for optimizing the capacity of an energy system and decoupling and coordinating the spatio-temporal scheduling of mobile resources under extreme disaster disturbances. Background Art

[0002] To enhance the resilience of the energy system, resistance and recovery are two key aspects, that is, a resilient energy system should be able to resist the impact of destructive events, especially catastrophic events, and recover quickly. On the one hand, the energy system needs to provide safer and more reliable power services for users under normal weather conditions; on the other hand, it also needs to overcome various disturbances during extreme weather events to ensure the continuity and stability of energy supply. Therefore, it is particularly urgent to develop new energy system planning methods and mobile resource scheduling schemes. By optimizing equipment capacity and realizing multi-energy complementarity using energy conversion equipment, the economy and sustainability of the energy system can be effectively improved. Especially during extreme weather or disasters, the flexible scheduling ability of mobile resources will become the key to meeting user needs and enhancing system resilience.

[0003] However, in the existing technical research on the optimization of the energy system configuration under extreme weather, the impact on renewable energy is often only considered, while the dynamic changes in the demand-side load and the collaborative scheduling potential of multiple types of mobile resources (such as emergency generators and in-vehicle energy storage) are ignored, resulting in insufficient system resilience and difficulty in effectively coping with the complex disturbances on both the supply and demand sides in disaster scenarios. Summary of the Invention

[0004] (I) Technical Problems to be Solved In view of the above-mentioned shortcomings and deficiencies of the existing technology, the present invention provides a method, system, device and medium for optimizing the capacity of an energy system and decoupling and coordinating the spatio-temporal scheduling of mobile resources under extreme disaster disturbances, which solves the technical problems that the existing energy system planning and mobile resource scheduling methods do not fully consider the dynamic change characteristics of user loads under extreme weather and insufficiently explore the collaborative scheduling potential of multiple types of power sources, resulting in limited optimization of the energy system.

[0005] (II) Technical Solutions To achieve the above object, the main technical solutions adopted by the present invention include: In a first aspect, an embodiment of the present invention provides a method for optimizing the capacity of an energy system and coordinately scheduling spatio-temporal decoupling of mobile resources under extreme disaster disturbances, including: according to the obtained extreme disaster event data, linking the equipment failure probability and line failure probability obtained through two-state resilience adaptive Markov and Monte Carlo processing, and coupling the load change parameters driven by extreme disaster events to form a disaster response database; converting the fluctuations in renewable energy output and the deviation of user load demand into fuzzy variables through fuzzy chance constraints, and establishing a capacity configuration optimization model in combination with the disaster response database, with the minimum power levelized cost as the optimization goal to generate a multi-energy complementary capacity configuration; based on the multi-energy complementary capacity configuration, introducing mobile resources to establish an elastic enhancement scheduling architecture, and generating an elastic enhancement plan adaptive to extreme disaster events through implementing damaged road network topology reconstruction, spatio-temporal transfer chain modeling, and action state coordinated control of mobile resources.

[0006] Optionally, according to the obtained extreme disaster event data, linking the equipment failure probability and line failure probability obtained through two-state resilience adaptive Markov and Monte Carlo processing, and coupling the load change parameters driven by extreme disaster events to form a disaster response database includes: According to the obtained extreme disaster event data of multiple observation points, a multi-dimensional extreme disaster event database is formed; wherein, the extreme disaster event data is a set of disaster events exceeding a preset extreme threshold, including earthquakes, storms, extremely high temperature and drought, and extremely low temperature and ice disasters. Using a two-state resilience adaptive Markov switching model to analyze the state transition chain of equipment, generating the equipment failure probability through the dynamic game of failure rate and repair rate, and generating a time-series equipment failure state set through Monte Carlo method for multi-scenario sampling simulation. According to the spatial distribution intensity and time accumulation parameters under the current extreme disaster event, simulating the dynamic diffusion process of the extreme disaster event through parameter space gradient mapping and time accumulation effect analysis, and generating a set of line failure probabilities including energy network lines and pipe networks. Constructing a dynamic load growth evolution model, and quantifying the load dynamic increment of extreme low temperature events or extreme high temperature events by introducing humidity correction terms, wind speed correction terms, and irradiance correction terms to form a time-series load dynamic parameter set. Aligning the equipment failure state set, the line failure probability set, and the load dynamic parameter set in space and time, and generating a disaster response database in combination with the multi-dimensional extreme disaster event database.

[0007] Optionally, the two-state resilience adaptive Markov switching model is: ; ; ; ; ; In the formula, is the probability of the equipment j operating normally, is the probability of the equipment j failing. 1 indicates that the operating state of the equipment j is normal, and 0 indicates that the equipment j has a fault during operation. , is the time-varying repair rate; for the power grid and gas network, , is the number of emergency repair teams, is the baseline repair rate without additional emergency repair resources, is the emergency repair baseline quantity; for the energy system, , and are the time-varying normal operating time of the equipment and the time-varying repair time of the equipment respectively, is a number randomly generated from the uniform distribution [0, 1], is a number randomly generated from the uniform distribution [0, 1], is the state of the equipment j at t time, is the failure rate, , K ={earthquake, storm, extreme high temperature and drought, extreme low temperature and ice disaster}, K is the set of disaster types, is the real-time intensity index of the disaster K , is the weight coefficient of the disaster K , satisfying w k = 1, is the non-linear influence function of the disaster K , which is used to standardize the physical dimensions of different disasters, is the quantified maintenance resource input disaster intensity maintenance and disaster game intensity that suppresses the ability of R ( t ) ≥ D ( t ) when γ ( t ) ≥ 0.5, indicating that maintenance is dominant and the probability of the normal state of the equipment is increased; R ( t ) < D ( t ) when γ ( t ) < 0.5, indicating that the disaster suppresses and increases the risk of failure.θ is the reference threshold value of the resource disaster ratio; The dynamic load growth evolution model is: ; In the formula, is the predicted load at time t when extreme weather occurs, is t the predicted load at time when extreme weather does not occur, is the load increment at time caused by extreme weather, related to the perceived temperature t ; is the neural network load prediction model, represents the actual temperature at time, t ; represents t the air relative humidity correction term at time, represents the wind speed correction term, represents the solar radiation correction term.

[0008] Optionally, convert the fluctuations in renewable energy output and the deviation of user load demand into fuzzy variables through fuzzy chance constraints, and establish a capacity configuration optimization model in combination with the disaster response database, with the goal of minimizing the levelized cost of electricity, and generate a multi-energy complementary capacity configuration including: Construct a membership function using fuzzy chance constraints to convert the random fluctuations in wind and solar power generation and the deviation of user load demand into a set of fuzzy variables under credible probability constraints; In response to regional energy demand, construct a system topology including energy supply equipment, energy storage equipment, and energy conversion units based on the disaster response database, and configure capacity configuration elements including equipment output sets, equipment status sets, and efficiency parameter sets; Establish a capacity configuration optimization model with the goal of minimizing the levelized cost of electricity, generate model constraint conditions based on the capacity configuration elements, and embed the set of fuzzy variables into the model constraint conditions to generate a fuzzy feasible solution space that integrates uncertainties; Iteratively solve the capacity configuration optimization model in the fuzzy feasible solution space, and output a multi-energy complementary capacity configuration adapted to the system topology.

[0009] Optionally, the equipment output set is: ; In the formula, is t the electricity price at time, is the electricity generation of equipment i at t time, is the cooling capacity of equipment i at t time, is the equipmenti At t the heating capacity at time is the hydrogen production amount of the electro-hydrogen conversion device i at t time are respectively t the building's electrical load demand, heat load demand, and cooling load demand at time is t the capacity of the energy storage device i at time The set of device states is expressed as: ; In the formula, represents the operating state of device i at t time. When is 1, it means that device i is in a normal operating state. When is 0, it means that device i is in a running fault state are respectively the charging state and discharging state of the energy storage device i at t time. When the energy storage device i is in a fault state, or is 0. When in a normal operating state, or is 1 ; The capacity configuration optimization model is: ; ; ; ; ; In the formula, TNPC is the total net present value, CRF is the capital recovery factor, is t the electrical load at time is t the charging load of the electric vehicle at time is t the reduced electrical load at time is the annual investment cost of the system, is the discount rate, is the device life, is the device i unit investment cost, is the device iThe unit investment cost, is the annual system operation and maintenance cost, is the unit operation and maintenance cost of the equipment i , and i is the operating power of the equipment.

[0010] Optionally, based on the multi - energy complementary capacity configuration, mobile resources are introduced to establish an elastic improvement scheduling architecture. By implementing the damaged road network topology reconstruction, spatio - temporal transfer chain modeling, and action - state collaborative regulation of mobile resources, an elastic improvement plan adaptable to extreme disaster events is generated, including: Based on the multi - energy complementary capacity configuration, mobile resources including mobile emergency power vehicles, electric vehicles, and diesel generators are introduced as elastic supplementary energy supply units, and an elastic improvement scheduling architecture for the coordination of fixed equipment and mobile resources is established; According to the obtained disaster damage information, damaged road network nodes and edge sets are removed, the road network topology is reconstructed, and a set of feasible scheduling paths that minimize the distance from the starting point to the target energy system for mobile resources is generated; The departure timestamp, expected arrival timestamp, driving duration interval, and charging stop time period of mobile resources are analyzed to generate time - chain features, and the starting position coordinates, target energy supply site coordinates, and path node coordinate sets of mobile resources are extracted to generate space - chain features; The time - chain features and space - chain features are verified. Based on the spatio - temporal feature combinations that pass the verification, the maximum service radius, shortest response time, and continuous power supply duration threshold of mobile resources are calculated to form a spatio - temporal schedulable boundary; According to the energy supply characteristics of the obtained mobile resources, an action set representing the energy supply behavior parameters of mobile resources is obtained, and a state set representing the operating state of mobile resources is constructed based on the mutual exclusion constraints of equipment operation; In the elastic improvement scheduling architecture, with the goals of minimizing the total cost of mobile resource scheduling, the lowest load curtailment, and the minimum path energy consumption, based on the spatio - temporal schedulable boundary, the optimal solutions of the parameters of the action set of mobile resources, the flag sequence of the state set, and the priority of the set of feasible scheduling paths are dynamically solved through a multi - layer collaborative optimization algorithm, and an elastic improvement plan adaptable to extreme disaster events is generated.

[0011] Optionally, in the elastic improvement scheduling architecture, with the goals of minimizing the total cost of mobile resource scheduling, the lowest load curtailment, and the minimum path energy consumption, based on the spatio - temporal schedulable boundary, the optimal solutions of the parameters of the action set of mobile resources, the flag sequence of the state set, and the priority of the set of feasible scheduling paths are dynamically solved through a multi - layer collaborative optimization algorithm, and an elastic improvement plan adaptable to extreme disaster events is generated, including: The parameters of the action set and the flag sequence of the state set of mobile resources are encoded as the position vectors of honey badger individuals to generate an initial honey badger population that meets the spatio - temporal schedulable boundary; Encode the priority of the set of feasible scheduling paths into the migration path chain of the individuals in the bird flock, generate an initial bird flock population that meets the spatio-temporal schedulable boundary, and load it into the path memory bank; Construct an inner and outer layer collaborative optimization framework. In the inner layer architecture, the upper layer honey badger population performs a two-objective optimization of the global scheduling cost and the load reduction amount. The lower layer bird flock population optimizes the priority of the set of feasible scheduling paths with the goal of minimizing the path energy consumption through the path memory bank sharing mechanism. In the outer layer architecture, use the simulated annealing algorithm model to receive the optimal solution output by the inner layer architecture as the initial solution and perform perturbation optimization; For each honey badger individual, match the lowest energy consumption sequence of the corresponding path chain from the initial bird flock population, calculate the total cost of mobile resource scheduling, the load reduction amount, and the path energy consumption index, aggregate them into a comprehensive performance index through normalized weights, and generate a dual-driven scale of prey odor concentration - migration gravity; Sort the dual populations according to the dual-driven scale, select the current global optimal solution, and combine the introduced dynamic density factor to control the honey badger and bird flock individuals to perform two-stage optimization within the neighborhood of the current global optimal solution; In the wide-area mining stage, expand the search step of the prey odor concentration in the upper layer honey badger population and traverse the combinations of the action set parameters and the state set flag sequences. In the lower layer bird flock population, update the topological weights in the path memory bank according to the path requirements of the honey badger individuals to generate a wide-area feasible path set; In the fine-grained mining stage, shrink the search step in the upper layer honey badger population to lock the local optimum of the combinations of the action set parameters and the state set flag sequences. In the lower layer bird flock population, perform a priority sorting of the wide-area feasible path set driven by the migration gravity; Perform constraint verification on the updated honey badger and bird flock individuals. If the updated honey badger individuals and bird flock individuals violate the spatio-temporal schedulable boundary, perform topological re-routing repair, adjust the node order and passing weights of the path chain, and apply bird flock gravity calibration to reassign priorities based on the path memory bank, and feedback the repaired individuals to the corresponding populations to update the distribution of the honey badger and the bird flock; After the iteration in the inner layer architecture, input the current optimal scheduling scheme into the outer layer simulated annealing algorithm model, generate a new solution through Gaussian perturbation and calculate the fitness difference, and jump out of the local optimum with a probabilistic acceptance strategy, and output a scheme that meets the spatio-temporal schedulable boundary conditions and is jointly optimal for multiple objectives, including the mobile resource action set parameters, the state set flag sequences, and the priority of the set of feasible scheduling paths.

[0012] In a second aspect, an energy system capacity optimization and mobile resource spatio-temporal decoupling collaborative scheduling system under extreme disaster disturbances provided by an embodiment of the present invention includes: a database construction module, configured to link the device failure probability and line failure probability obtained through dual-state resilience adaptive Markov and Monte Carlo processing according to the acquired extreme disaster event data, and couple the load change parameters driven by extreme disaster events to form a disaster response database; a configuration generation module, configured to convert the fluctuations in renewable energy output and the deviation of user load demand into fuzzy variables through fuzzy chance constraints, and establish a capacity configuration optimization model in combination with the disaster response database, with the goal of minimizing the levelized cost of electricity, to generate a multi-energy complementary capacity configuration; a configuration improvement module, configured to introduce mobile resources based on the multi-energy complementary capacity configuration to establish an elastic improvement scheduling architecture, and generate an elastic improvement plan adaptive to extreme disaster events by implementing damaged road network topology reconstruction, spatio-temporal transfer chain modeling, and action state collaborative regulation of mobile resources.

[0013] In a third aspect, an automatic annotation device for high-quality drainage pipeline image data based on a deep learning model with multi-feature fusion provided by an embodiment of the present invention includes: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to execute the method for optimizing the capacity of the energy system and spatio-temporally decoupling and collaboratively scheduling mobile resources under extreme disaster disturbances as described above.

[0014] In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium, on which computer-executable instructions are stored, and when the executable instructions are executed by a processor, the method for optimizing the capacity of the energy system and spatio-temporally decoupling and collaboratively scheduling mobile resources under extreme disaster disturbances as described above is implemented.

[0015] (III) Beneficial effects The beneficial effects of the present invention are as follows: First, by integrating multi-source extreme disaster event observation data, a dual-state resilience adaptive Markov process is used to dynamically simulate the device failure probability and line failure risk, and a dynamic incremental model of user cooling / heating / electricity load driven by meteorological parameters is coupled to construct a disaster response database. This database breaks through the limitations of the decoupling of meteorological data from device failure models and user load characteristics in traditional methods. This database not only covers multi-dimensional failure mechanisms on the device side, but also couples the time-varying incremental model of cooling / heating / electricity load driven by meteorological factors, breaking through the limitations of one-sided disaster risk assessment or static load prediction in traditional planning, and deeply exploring the impacts of extreme weather and natural disasters on the supply side and demand side.

[0016] On this basis, based on fuzzy chance constraints, the fluctuations in renewable energy output and user load ranges are converted into credibility fuzzy variables. Combining with the set of fault scenarios generated by the disaster response database, a capacity configuration optimization model with the goal of minimizing the levelized cost of electricity is constructed. Through multi-energy coupling constraints, embedding of equipment failure probabilities, and definition of fuzzy chance constraint boundaries, this model realizes the coordinated planning of the capacities of energy supply equipment, energy storage devices, and energy conversion equipment under extreme weather conditions. This method reduces the total system cost while ensuring the reliability of energy supply in scenarios with a high proportion of renewable energy penetration, and solves the problem of insufficient resilience optimization in existing planning caused by simplifying meteorological disaster risk factors into fixed thresholds. Through this comprehensive optimization, the adaptability and resilience of the system in the face of extreme climate conditions or natural disasters are significantly improved, and while pursuing system economy, the ability of the system to respond to sudden extreme events is enhanced.

[0017] Furthermore, on the basis of multi-energy complementary capacity optimization, three types of mobile resources, namely mobile emergency power vehicles, electric vehicles, and diesel generators, are innovatively introduced. Elasticity enhancement is achieved through the implementation of damaged road network topology reconstruction, spatio-temporal transfer chain modeling, and action state coordination mechanisms for mobile resources.

[0018] Thus, through the collaborative decision-making of device capacity optimization driven by the disaster response database and dynamic scheduling of mobile resources, the present invention realizes a dual elastic gain of strengthening the disaster resistance ability of fixed devices and rapid post-disaster response of mobile resources under extreme weather conditions, providing an innovative solution for the design and operation of highly resilient energy systems. Description of the Drawings

[0019] Figure 1 Schematic diagram of the composition of the energy system for the method provided by the embodiment of the present invention; Figure 2 Schematic diagram of the process for the method provided by the embodiment of the present invention; Figure 3 Specific process schematic diagram of step S1 of the method provided by the embodiment of the present invention; Figure 4 Schematic diagram of device failure transfer of the dual-state resilience adaptive Markov for the method provided by the embodiment of the present invention; Figure 5 Schematic diagram of the Monte Carlo-based system device failure extraction process for the method provided by the embodiment of the present invention; Figure 6 Specific process schematic diagram of step S2 of the method provided by the embodiment of the present invention; Figure 7 Specific process schematic diagram of step S3 of the method provided by the embodiment of the present invention; Figure 8 Road network path diagram for the method provided by the embodiment of the present invention; Figure 9 It is a schematic flowchart of the specific process of step S35 of the method provided by the embodiment of the present invention; Figure 10 It is a schematic diagram of the elastic effect of the energy system of the method provided by the example of the present invention; Figure 11 It is a schematic overall flowchart of the method provided by the embodiment of the present invention. Detailed implementation manners

[0020] In order to better explain the present invention for easy understanding, the present invention will be described in detail below in conjunction with the accompanying drawings through specific implementation manners.

[0021] Prior to this, in order to facilitate the understanding of the technical solution provided by the present invention, the concept of the energy system will be introduced first below.

[0022] It should be clear that, as shown in Figure 1 The energy system in the embodiment of the present invention includes: an input end, a system end and a demand end; among them, the input end includes a power grid, solar energy, wind energy, geothermal energy and natural gas, and the energy conversion equipment (electrolyzer, photovoltaic panel, fan, ground source heat pump, gas turbine, electric chiller, fuel cell and waste heat boiler, etc.) in the system end is combined to realize the conversion between energies to meet the cooling, heating, power and hydrogen loads of each building at the demand end. At the same time, relying on energy storage equipment such as hydrogen storage tanks, storage batteries, chilled water storage tanks and heat storage tanks to achieve dynamic balance between supply and demand; using mobile resources such as electric vehicles, mobile emergency vehicles, mobile diesel generators, etc. to ensure the load supply when the power supply is insufficient due to extreme weather and other reasons, and improve the system elasticity.

[0023] The energy system equipment has a strong coupling relationship. The operation strategy of the established energy system is proposed, which is to determine heat based on cooling, determine hydrogen based on heat, and determine electricity based on hydrogen; electricity is used as the bottom-layer energy, and the consumption of renewable energy is given priority, followed by power supply from gas turbines. The excess electricity is stored in the storage battery, and the insufficient electricity is purchased from the power grid; the building cooling load is mainly supplied by the electric chiller and the ground source heat pump. When the rated power of the ground source heat pump is insufficient, the electric chiller is started, and the excess cooling capacity is stored in the chilled water storage tank; the building heating load is mainly supplied by the waste heat recovery of the gas turbine and the ground source heat pump. When the waste heat of the gas turbine cannot meet the demand, the ground source heat pump is started to supply heat, and the excess heat is stored in the heat storage tank; when the renewable energy and the gas turbine cannot meet the electricity load, the electricity consumed by the electrolyzer for hydrogen production is purchased from the power grid.

[0024] On this basis, as shown in Figure 2As shown in the figure, a method for optimizing the capacity of an energy system and coordinately scheduling the spatio-temporal decoupling of mobile resources under extreme disaster disturbances proposed by an embodiment of the present invention includes: According to the obtained extreme disaster event data, linking the device failure probability and line failure probability obtained through the dual-state resilience adaptive Markov and Monte Carlo processing, and coupling the load change parameters driven by extreme disaster events to form a disaster response database; Converting the fluctuations in the output of renewable energy and the deviation of user load demand into fuzzy variables through fuzzy chance constraints, and establishing a capacity configuration optimization model in combination with the disaster response database, with the minimum power levelized cost as the optimization goal, to generate a multi-energy complementary capacity configuration; Based on the multi-energy complementary capacity configuration, introducing mobile resources to establish an elastic enhancement scheduling architecture, and generating an elastic enhancement plan adaptable to extreme disaster events through implementing the topological reconstruction of the damaged road network of mobile resources, spatio-temporal transfer chain modeling, and coordinated regulation of action states.

[0025] First, by integrating multi-source extreme disaster event observation data, using the dual-state resilience adaptive Markov process to dynamically simulate the device failure probability and line failure risk, and coupling the dynamic incremental models of user cooling / heating / electricity loads driven by meteorological parameters, a disaster response database is constructed. This database breaks through the limitations of the decoupling of meteorological data, device failure models, and user load characteristics in traditional methods. This database not only covers the multi-dimensional failure mechanisms on the device side, but also couples the time-varying incremental models of cooling / heating / electricity loads driven by meteorological factors, breaking through the limitations of one-sided disaster risk assessment or static load prediction in traditional planning, and deeply exploring the impacts of extreme weather and natural disasters on the supply side and demand side.

[0026] On this basis, based on fuzzy chance constraints, converting the fluctuations in the output of renewable energy and the fluctuation range of user loads into credibility fuzzy variables, and combining with the set of failure scenarios generated by the disaster response database, a capacity configuration optimization model with the minimum power levelized cost as the goal is constructed. This model realizes the coordinated planning of the capacities of energy supply devices, energy storage devices, and energy conversion devices under extreme weather through multi-energy coupling constraints, embedding of device failure probabilities, and definition of fuzzy chance constraint boundaries. This method reduces the total system cost while ensuring the energy supply reliability in the scenario of high penetration of renewable energy, and solves the problem of insufficient elastic optimization caused by the simplification of meteorological disaster risk factors to fixed thresholds in existing planning. Through this comprehensive optimization, the adaptability and resilience of the system in the face of extreme climate conditions or natural disasters are significantly improved, and while pursuing the system economy, the ability of the system to cope with sudden extreme events is enhanced.

[0027] Furthermore, on the basis of multi-energy complementary capacity optimization, three types of mobile resources, namely mobile emergency power vehicles, electric vehicles, and diesel generators, are innovatively introduced, and elastic enhancement is achieved through implementing the topological reconstruction of the damaged road network of mobile resources, spatio-temporal transfer chain modeling, and coordinated action state mechanism.

[0028] Thus, through the collaborative decision-making of disaster response database-driven equipment capacity optimization and dynamic scheduling of mobile resources, the present invention realizes the dual elastic gain of strengthening the disaster resistance of fixed equipment and the rapid post-disaster response of mobile resources under extreme weather, providing an innovative solution for the design and operation of a highly resilient energy system.

[0029] To better understand the above technical solution, the exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings.

[0030] Specifically, the embodiment of the present invention provides a method for collaborative scheduling of energy system capacity optimization and spatio-temporal decoupling of mobile resources under extreme disaster disturbances, which includes: S1. According to multi-source extreme disaster event data, link the equipment failure probability and line failure probability under extreme disaster events obtained through two-state resilience adaptive Markov and Monte Carlo processing, and couple the weather-driven load change parameters to form a disaster response database.

[0031] Further, as Figure 3 shown, step S1 includes: S11. According to the extreme disaster event data collected at multiple observation points, form a multi-dimensional extreme disaster event database; wherein, the extreme disaster event data is a set of meteorological events including earthquakes, storms, extremely high temperature droughts, and extremely low temperature ice disasters that exceed the corresponding extreme thresholds.

[0032] In this step, historical extreme weather data of multiple meteorological observation stations in different regions is selected, and the collected extreme weather includes earthquakes, typhoons, droughts caused by extremely high temperatures, and ice disasters caused by extremely low temperatures.

[0033] Specifically, for earthquake disasters, relevant parameters such as earthquake time, magnitude, and epicenter area intensity need to be collected, and the collected earthquake data is mainly medium-level earthquakes above magnitude 5; for typhoon disasters, relevant parameters such as solar radiation intensity and real-time wind speed need to be collected; for high-temperature drought weather, relevant meteorological parameters such as the temperature and solar irradiance intensity on the day when the highest temperature of the day exceeds the extreme threshold need to be collected; for low-temperature ice disasters, relevant meteorological parameters such as the temperature of the frozen weather, real-time wind speed, freezing rain time, and freezing rain rate need to be collected.

[0034] Furthermore, analyze the impacts of extreme weather on the supply side and load side of the energy system: When an earthquake disaster occurs, the power lines are damaged and there is no power supply from the external power grid. The ground source heat pumps are damaged and unable to provide heating and cooling. The photovoltaic and wind turbines are damaged and unable to supply power. Some gas pipelines fail, and some roads are damaged or congested, resulting in changes in the dispatching routes of mobile emergency vehicles and affecting the restoration of load supply. Typhoon weather mainly affects overhead power lines, wind turbines, and photovoltaics. Extreme high-temperature weather often causes drought, leading to a decrease in power supply from the external power grid, limited output of photovoltaics and wind turbines. When the temperature is higher than 15°C, the efficiency of gas turbines decreases by 0.7% - 1% for every 1°C increase in temperature, and the cooling load of users increases when high-temperature weather occurs. Extreme low-temperature weather will result in ice and snow disasters. At this time, some overhead lines fail, the output of photovoltaics and wind turbines is limited, and the heating load of users increases when low-temperature weather occurs.

[0035] S12. Use the two-state resilient adaptive Markov process to analyze the state transition chain of the equipment, generate the equipment failure probability through the dynamic game of the failure rate and repair rate, and perform multi-scenario sampling simulation through the Monte Carlo method to generate a time-series set of equipment failure states.

[0036] In a specific example, referring to Figure 4 , to improve the resilience and reliability of the energy system, a method combining the two-state resilient adaptive Markov process and Monte Carlo simulation is used to generate the time-sequence scenarios of equipment failures. As shown in the two-state resilient adaptive Markov state transition model in Figure 4 , its core lies in constructing a time-varying transition mechanism for the operating state (normal / fault) of the equipment through the dynamic interaction of the failure rate and repair rate, and generalizing it to the external power grid circuit and natural gas underground pipeline. The two-state resilient adaptive Markov switching model is specifically as follows: ; ; ; ; ; In the formula, is the probability that the equipment j operates normally, is the probability that the equipment j fails. 1 represents that the operating state of the equipment j is normal, and 0 represents that the equipment j has a failure in operation. , is the time-varying repair rate; for the power grid and gas network, , is the number of repair teams, is the benchmark repair rate without additional repair resources, is the emergency repair reference quantity; for the energy system, , or are the normal operation time of the time-varying equipment and the time-varying repair time of the equipment respectively, is a number randomly generated from the uniform distribution [0, 1], is a number randomly generated from the uniform distribution [0, 1], is the state of the equipment j at t moment, is the failure rate, , K = {earthquake, storm, extreme high temperature and drought, extreme low temperature and ice disaster}, K is the set of disaster types, is the disaster K real-time intensity index, is the weight coefficient of the disaster K satisfying w k = 1, is the non-linear influence function of the disaster K used to standardize the physical dimensions of different disasters, is the quantified maintenance resource input for the suppression ability of the disaster intensity maintenance and disaster game intensity, R ( t ) ≥ D ( t ) when, γ ( t ) ≥ 0.5, indicating that maintenance is dominant and the probability of the normal state of the equipment is increased; R ( t ) < D ( t ) when, γ ( t ) < 0.5, indicating that the disaster suppresses and increases the failure risk, is the reference threshold of the resource-disaster ratio.

[0037] Furthermore, referring to Figure 5 , the steps of fault extraction based on Monte Carlo are as follows: (1) Define the model parameters: First, define the state. State 0 indicates that the equipment is in the normal working state, and state 1 indicates that the equipment is in the operating failure state; Second, calculate the probability j of the normal working state of the equipment and the probability j of the operating failure state of the equipment ; Determine the time interval. In this embodiment, the time step is set to 1 hour. Usually under normal circumstances, the time step of the operation cycle is shorter than the maintenance time of the equipment. Therefore, within this time step, the equipment is considered irreparable and in a failure state. Finally, obtain the j interruption model of the equipment.

[0038] (2) Initialize the simulation: First, set the initial state of the equipment. Usually, the initial state of the equipment is considered to be the working state; second, determine the number of Monte Carlo simulations. Initialize 100 simulations to obtain sufficient data for analysis.

[0039] (2) Conduct the simulation: Each simulation starts from the initial state. At each time step, randomly determine whether the equipment undergoes a state transition according to the probability of the current equipment and the transition state, and record the state of the equipment and the time of the transition.

[0040] (3) Data collection: Collect the time when the equipment enters the failure state in each simulation, and calculate key indicators such as the failure rate and the mean time to failure.

[0041] In this embodiment, by breaking through the limitations of the traditional static reliability model through the dynamic failure rate - repair rate game mechanism, combined with Monte Carlo iteration under the constraint of the disaster duration, high-fidelity simulations of complex scenarios such as earthquake-induced cascading failures of gas pipelines and hourly line breaks caused by typhoons are realized, providing high-confidence input boundary conditions for subsequent capacity configuration optimization.

[0042] S13. According to the spatial distribution intensity and time accumulation parameters under the current extreme disaster event, simulate the dynamic diffusion process of the extreme disaster event through parameter space gradient mapping and time accumulation effect analysis, and generate a set of line failure probabilities including energy network lines and pipe networks.

[0043] Specifically, for different types of extreme disaster events, the following steps are used to dynamically generate a set of failure probabilities of energy network lines and pipe networks: (1) In the earthquake scenario, based on the spatial gradient distribution of the intensity in the epicenter area, combined with the correlation parameters fitted from historical earthquake data, analyze the seismic damage risk gradient of different sections along the gas pipeline as follows: ; In the formula, is the earthquake occurrence frequency, is the earthquake magnitude, and its size depends on the intensity in the earthquake epicenter area ; a and b are fitting parameters, and the fitting parameters in different regions need to be obtained by fitting local historical earthquake data.

[0044] Judge the earthquake hazard rate according to the intensity in the epicenter area, and thus calculate the failure probability of the gas pipeline: ; ; In the formula, the intensity of the earthquake epicenter area, is the earthquake hazard rate, is the failure probability of the gas pipeline, is the pipeline length.

[0045] In the earthquake scenario, the energy supply system has suffered serious damage, presenting a complex and severe energy supply and demand situation. The characteristics of the energy supply of the energy system are as follows: 1) Interruption and limitation of energy supply: The external power grid lines and natural gas pipelines are damaged, resulting in a complete interruption of the external power grid power supply, a significant reduction in the underground natural gas pipelines, and directly restricting the power generation of gas turbines. At the same time, renewable energy sources such as photovoltaic and wind turbines are damaged by the earthquake and cannot supply power normally, and the energy supply method has been severely hit. 2) Associated impact of the energy supply system: The damage to the power supply system not only leads to insufficient power supply but also has a chain reaction on building heating and cooling. The ground source heat pump is damaged and cannot supply heating and cooling normally. The electric chiller still needs to consume electricity to operate under the condition of tight cooling supply, further exacerbating the contradiction between power supply and demand. 3) Adjustment of the energy supply strategy: In the extremely tight power supply situation, in order to give priority to meeting the energy demand of buildings, the electrolytic cell originally used for hydrogen production stops working, and the limited resources are centrally allocated. The electrical load on the demand side is mainly met by the discharge of batteries and mobile resources. At the same time, some non-critical electrical loads are restricted or reduced.

[0046] (2) In the typhoon scenario, extract the ratio of the real-time wind speed to the designed wind speed under the typhoon movement path, identify the designed wind speed threshold and the duration of the typhoon weather. The failure probability of the wire under the typhoon is: ; ;

[0047] In the formula, is the wind speed at time t ; is the designed wind speed; is the maximum wind speed that the line can withstand; is the length of the line; is the duration of the typhoon weather.

[0048] In the typhoon scenario, the characteristics of energy supply among energy system equipment are as follows: 1) Fault simulation and detection: Using the Batts typhoon model, integrating the typhoon wind field model and the typhoon movement attenuation model, and obtaining typhoon data such as wind speed and wind direction information to dynamically simulate and characterize the transmission line fault problem, providing a basis for the subsequent response of power supply equipment. 2) Energy supply interruption: When the typhoon passes through the external lines of the energy system and causes damage, it leads to a complete interruption of the external power grid supply. At the same time, renewable energy sources such as photovoltaic and wind turbines are damaged due to excessive wind speed and cannot supply power normally. At this time, the gas turbine undertakes the main power supply task, and the battery and mobile resources cooperate as emergency power supply equipment to jointly ensure power supply in the typhoon situation. 3) Adjustment of energy supply strategy: The main force of power supply changes to the gas turbine, which undertakes the task of meeting the main power load. The battery and mobile resources cooperate as emergency power supply equipment to jointly ensure power supply in the typhoon situation and relieve the contradiction between power supply and demand.

[0049] (3) In the ice storm scenario, according to the spatial heterogeneity of the freezing rain rate and the persistence of the freezing rain time, calculate the spatio-temporal evolution trajectory of the ice thickness on the transmission line as follows: ; ;

[0050] In the formula, is the thickness of the ice; is the freezing rain time; is the freezing rain rate; is the density of the ice; is the density of the freezing rain; is the liquid water content, calculated by the formula ; is the ice load per unit wire line; is the conductor diameter; is the wind load; is a constant coefficient, in the unit of ; is the cross-boundary factor; is the combined ice ball load, calculated from the ice load and the wind load.

[0051] Based on the spatio-temporal evolution trajectory of the ice thickness, establish the following probability model of the transmission line under the ice storm: ; ; In the formula, is the fault probability per unit line; and are the first threshold and the second threshold of the ice ball load respectively.

[0052] In the extreme weather scenario of ice storms, the heating load increases, and the energy supply operation and linkage in the equipment room present the following specific states: 1) Energy supply is blocked: The ice storm causes the photovoltaic and wind turbines to be affected by the ice thickness, and their power generation capacity is severely limited. At the same time, the external lines are overloaded due to icing, which limits the power supply capacity of the power grid. The normal power supply system is impacted, and various original energy supply channels are ineffective or weakened to varying degrees, and the overall energy supply is in trouble; 2) Energy supply dominance and cooperation: The gas turbine becomes the main power supply equipment during the ice storm. When the power generation of the gas turbine still cannot meet the power demand, mobile resources are quickly put into use to supplement the power.

[0053] S14. Construct a temperature-driven load growth model. By introducing humidity correction terms, wind speed correction terms, and irradiance correction terms, quantify the load dynamic increment of extreme low temperature events or extreme low temperature events, and form a time-series load dynamic parameter set.

[0054] Specifically, the dynamic load growth evolution models for drought weather caused by high temperature and ice storm weather caused by low temperature are as follows: ; In the formula, is t the predicted load at the moment when extreme weather occurs, is t the predicted load at the moment when extreme weather does not occur, is the load increment at the moment t after extreme weather occurs, which is related to the perceived temperature ; is the neural network load prediction model, and the predicted load is related to the perceived temperature. represents t the actual air temperature at the moment, represents t the air relative humidity correction term at the moment, represents the wind speed correction term, represents the solar radiation correction term. Among them, the neural network load prediction model is a data-driven mathematical function, and its core function is to predict the load under specific meteorological conditions by learning the hidden laws between historical meteorological data and power load. The predicted load output by the model reflects the quantitative impact of the combined effects of temperature, humidity, wind speed, and radiation under normal meteorological conditions on electricity consumption behavior. The neural network load prediction model can adopt network structures such as LSTM (Long Short-Term Memory Network), CNN (Convolutional Neural Network), or Transformer.

[0055] In extremely high temperature scenarios, all types of power generation equipment are significantly affected, as follows: 1) Limited renewable energy power generation: Excessively high temperatures will reduce the conversion efficiency of photovoltaic modules, resulting in reduced power generation. At the same time, high temperatures may affect the normal operation of mechanical parts and electronic components inside the wind turbine, thereby limiting the power generation capacity of the wind turbine, making it difficult for photovoltaics and wind turbines, which were originally important sources of renewable energy power generation, to fully exert their efficiency under extremely high temperatures. 2) Limited traditional energy power generation: Under extremely high temperatures, gas turbine power generation is also limited. Excessively high temperatures will reduce the density of air entering the gas turbine, resulting in a decrease in combustion efficiency, thereby limiting the power generation of the gas turbine and making it difficult to maintain normal power generation levels.

[0056] S15. The equipment fault state set, line failure probability set and load dynamic parameter set are aligned in time and space, and combined with the multi-dimensional extreme disaster event database to generate a disaster response database. Through the time and space grid alignment technology, the equipment fault state set, line failure probability set and load dynamic parameter set are multi-dimensionally coupled to generate a disaster response database.

[0057] S2. Through fuzzy opportunity constraints, the fluctuation of renewable energy output and the deviation of user load demand are converted into fuzzy variables. Combined with the disaster response database, a capacity configuration optimization model is established. Taking cost minimization as the optimization goal, multi-energy complementary capacity configuration is generated.

[0058] Furthermore, if Figure 6 As shown, step S2 includes: S21. In response to regional energy demand, a system topology including energy supply equipment, energy storage equipment and energy conversion units is constructed based on a disaster response database, and capacity configuration elements including equipment output set, equipment status set and efficiency parameter set are configured.

[0059] Specifically, the equipment output set is used to characterize the dynamic balance state of energy production and demand of the energy system at time t. It includes the energy supply efficiency of energy supply equipment such as wind turbines and photovoltaics in the energy system, building electrical load demand, heat load demand and cooling load demand, hydrogen fuel cell vehicle hydrogen, load demand, energy storage equipment SOC, etc., to fully reflect the output information of the energy system at the current moment. Therefore, the output of each device in the energy system at time t is expressed as: ; In the formula, for t Electricity price at any time, For equipment i exist t The amount of electricity generated at the time; For equipment i exist t Cooling capacity at any time; For equipment i existt The heating capacity at a certain moment; For the electro-hydrogen conversion device i At t The hydrogen production capacity at a certain moment; They are respectively the building's electrical load demand, heating load demand, and cooling load demand at time t; Is the capacity of energy storage device i at time t.

[0060] The device status set is used to describe the operation reliability of system devices and the charging and discharging dynamics of energy storage. It includes the operation status (normal operation or fault status) of all devices in the energy system and the charging and discharging status of energy storage devices, etc., to comprehensively reflect the operation status information of each device in the energy system at the current moment. Therefore, t The status of each device in the energy system at a certain moment is expressed as: ; In the formula, Indicates the operation status of device i At t A certain moment, obtained from formula (14). When it is 1, it means that device i is in a normal operation state. When Is 1, it means that device i is in a normal operation state. When Is 0, it means that device i Is in a running fault state; They are respectively the charging state and discharging state of energy storage device i At t A certain moment. When the energy storage device i Is in a fault state, Or Is 0. When in a normal operation state, Or Is 1, and .

[0061] The energy efficiency parameter set is used to quantify the energy conversion and transmission losses of devices, and it includes: the efficiency of energy supply devices, the energy storage loss coefficient, and the energy conversion efficiency.

[0062] S22. Establish a capacity configuration optimization model with the goal of minimizing the levelized cost of electricity, generate model constraint conditions based on capacity configuration elements, and embed the fuzzy variable set into the model constraint conditions to generate a fuzzy feasible solution space that integrates uncertainties.

[0063] In a specific embodiment, collect the unit investment cost, discount rate, and investment period of each device in the energy system, collect the unit operation and maintenance cost of each device, determine the power supply and electricity load reduction at each moment, and construct a capacity configuration optimization model for improving the resilience of the energy system considering extreme weather and natural disasters; the capacity configuration optimization model aims to minimize the Levelized Cost of Electricity (LCOE), which represents the average electricity price per kilowatt-hour and is calculated from the Total net present cost (TNPC), Capital recovery factor (CRF), and total electricity consumption, and evaluate the annual total cost of the system over a 20-year service life, including initial investment, operation, and maintenance costs, as follows: ; ; ; ; ; In the formula, TNPC is the total net present cost, CRF is the capital recovery factor, is t the electricity load at time is t the charging load of electric vehicles at time is t the reduced electricity load at time is the annual investment cost of the system, is the discount rate, is the equipment life, is the equipment i unit investment cost, is the equipment i installed capacity, is the annual operation and maintenance cost of the system, is the equipment i unit operation and maintenance cost, is the equipment i operating power.

[0064] S23. Iteratively solve the capacity configuration optimization model in the fuzzy feasible solution space, and output the multi-energy complementary capacity configuration adapted to the system topology.

[0065] S3. Based on the multi - energy complementary capacity configuration, introduce mobile resources to establish an elastic improvement scheduling framework. Through the implementation of damaged road network topology reconstruction, spatio - temporal transfer chain modeling, and action - state collaborative regulation of mobile resources, generate an elastic improvement plan adaptable to extreme disaster events.

[0066] Further, as Figure 7 shown, step S3 includes: S31. On the basis of the multi - energy complementary capacity configuration, introduce mobile resources including mobile emergency power vehicles, electric vehicles, and diesel generators as elastic supplementary energy supply units, and establish an elastic improvement scheduling framework for the coordination of fixed equipment and mobile resources.

[0067] S32. According to the obtained disaster damage information, remove the damaged road network nodes and edge sets, reconstruct the road network topology, and generate a set of feasible scheduling paths that minimize the distance from the starting point to the target energy system for mobile resources.

[0068] In this embodiment, as Figure 8 shown, the road topology network includes nodes A - M and edge sets, where the edge weights represent the distances between adjacent nodes. When a road is damaged due to an earthquake disaster, the scheduling path of mobile emergency resources is blocked, resulting in a delay in load restoration. To address this problem, the following steps are executed: First, based on the road damage information collected after the earthquake (such as bridge collapse, road surface cracking), identify the damaged nodes and edge sets (such as the road interruption between nodes D - E). Second, remove the damaged edge sets and update the road network adjacency matrix to generate a sparse road network structure under the post - disaster reachability constraint; furthermore, use Dijkstra's algorithm to search for the optimal scheduling path of the mobile emergency power vehicle from the starting point (such as node A) to the target energy system (such as node M) in the reconstructed road network, and use Dijkstra's algorithm to search for the shortest driving route of mobile resources, and output a set of feasible paths that minimize the distance (such as A→B→F→L→M, A→B→G→K→M, A→C→H→J→M, etc.).

[0069] S33. Analyze the departure timestamp, estimated arrival timestamp, driving duration interval, and charging stop times of mobile resources to generate time - chain features including time availability constraints, and extract the starting position coordinates, target energy supply station coordinates, and path node coordinate sets of mobile resources to generate space - chain features including space connectivity constraints.

[0070] In this step, the travel chain of electric vehicles is divided into two parts: the time chain and the space chain; the time chain mainly describes the time distribution of electric vehicle travel, including: departure time ( ), arrival time ( ), driving time ( ), and parking time ( ); the space chain mainly describes the spatial transfer of electric vehicle travel, including: starting point ( ), end point ( ), driving path ( ); The time chain and space chain are expressed as: ; ; S34. Verify the time chain characteristics and space chain characteristics. Based on the spatio-temporal feature combinations that pass the verification, calculate the maximum service radius, shortest response time, and continuous power supply duration threshold of the mobile resources, and form the spatio-temporal schedulable boundary. Specifically, establish the association rules between the time chain and the space chain, and implement the spatio-temporal consistency verification rules and spatio-temporal conflict detection to screen out the scheduling schemes with overlapping time window conflicts, path node failures, or unavailable charging periods at stops.

[0071] S35. Obtain the action set representing the power supply behavior parameters of the mobile resources according to the energy supply characteristics of the obtained mobile resources, and construct the state set representing the operating state of the mobile resources according to the mutual exclusion constraints of the device operation.

[0072] In this step, the mobile resource action set mainly includes the charging and discharging power of electric vehicles ( ), and the electric power consumed during driving ( ), the power supply power of mobile emergency power supply vehicles ( ), and the power generation power of mobile diesel generators ( ), so as to comprehensively represent the situation of scheduling mobile resources at a certain moment to restore the building energy load after extreme weather or natural disasters. The action sets of mobile resources such as electric vehicles (Electric vehicles, EVs), mobile emergency vehicles (Mobile energystorage system, MESSs), and mobile diesel generators (Diesel Generator, DGs) are expressed as: ; The mobile resource state set mainly includes the charging state of electric vehicles ( ), and the discharging state ( ), the start-stop state of mobile emergency vehicles ( ), the start-stop state of mobile diesel generators ( ), so as to comprehensively represent the scheduling states of mobile resources such as electric vehicles; the mobile resource state set is expressed as: ; In the formula, when the i th electric vehicle is in the charging state the value is 1, when the electric vehicle is in the discharging state the value is 1, and it cannot be charged and discharged simultaneously ; When the mobile emergency power supply vehicle arrives at the energy system and discharges to supply the building energy demand The value is 1; conversely, when there is no need to move the emergency vehicle for discharging or the mobile emergency power vehicle has not arrived at the energy system The value is 0; when the mobile diesel generator is in the power supply state The value is 1; conversely, it is 0.

[0073] S36. In the elastic promotion scheduling architecture, with the goals of minimizing the total cost of mobile resource scheduling, minimizing the load curtailment amount, and minimizing the path energy consumption, based on the spatio-temporal schedulable boundary, through a multi-layer collaborative optimization algorithm, dynamically solve the parameters of the action set, the flag sequence of the state set, and the optimal solution of the priority of the feasible scheduling path set of mobile resources, and generate an elastic promotion plan adaptable to extreme disaster events.

[0074] In order to quickly restore the load supply and further improve the system elasticity, with the core optimization goals of minimizing the operating cost of mobile resources, maximizing the load restoration amount, and minimizing the path energy consumption within the scheduling period, construct the following multi-objective function system: ; ; ; ; ; In the formula, represents the number of days of continuous occurrence of extreme weather or natural disasters; , , and are the scheduling costs of EV, DG, and MESS and the amount of unsupplied load respectively; is the path energy consumption during the scheduling of electric vehicles; , and are the unit energy supply prices respectively; and are the discharge amounts of EV and MESS respectively; is the amount of diesel consumed by DG for power supply; , , are the amounts of unsupplied electric load, unsupplied cooling load, and unsupplied heating load respectively.

[0075] Furthermore, as shown in Figure 9 , step S35 includes: S351. Encode the action set parameters and the flag sequence of the state set of mobile resources into the position vector of honey badger individuals, and generate an initial honey badger population that meets the spatio-temporal schedulable boundary. And set the algorithm parameters and the maximum number of iterations for each layer.

[0076] S352. Encode the priority of the set of feasible scheduling paths into the migration path chain of the bird flock individuals, generate the initial bird flock population that meets the spatio-temporal schedulable boundary, and load it into the path memory bank.

[0077] S353. Construct an inner and outer layer collaborative optimization framework. In the inner layer architecture, the upper layer honey badger population performs the dual-objective optimization of the global scheduling cost and the load reduction amount. The lower layer bird flock population optimizes the priority of the set of feasible scheduling paths with the goal of minimizing the path energy consumption through the path memory bank sharing mechanism. In the outer layer architecture, the simulated annealing algorithm model receives the optimal solution output by the inner layer architecture as the initial solution and performs perturbation optimization.

[0078] S354. For each honey badger individual, match the lowest energy consumption sequence of the corresponding path chain from the initial bird flock population, calculate the total cost of mobile resource scheduling, the load reduction amount, and the path energy consumption index, and aggregate them into a comprehensive performance index through normalized weights to generate a dual-driven scale of prey odor concentration - migration attraction.

[0079] S355. Sort the dual populations according to the dual-driven scale, select the current global optimal solution, and combine the introduced dynamic density factor to control the two-stage optimization of the honey badger and bird flock individuals within the neighborhood of the current global optimal solution.

[0080] S356. In the wide-area mining stage, expand the search step of the prey odor concentration in the upper layer honey badger population, traverse the combinations of the action set parameters and the state set flag sequences. In the lower layer bird flock population, update the topological weights in the path memory bank according to the path requirements of the honey badger individuals to generate a wide-area feasible path set.

[0081] S357. In the fine-grained mining stage, shrink the search step in the upper layer honey badger population, lock the local optimum of the combination of the action set parameters and the state set flag sequences. In the lower layer bird flock population, perform the priority sorting driven by migration attraction on the wide-area feasible path set.

[0082] S358. Perform constraint verification on the updated honey badger and bird flock individuals. If the updated honey badger individuals and bird flock individuals violate the spatio-temporal schedulable boundary, perform topological re-routing repair, adjust the node order and passing weights of the path chain, and apply bird flock gravity calibration to reassign priorities based on the path memory bank. Feed the repaired individuals back to the corresponding populations to update the distribution of the honey badger and the bird flock.

[0083] S359. After the inner-layer architecture is iterated, the current optimal scheduling scheme is input into the outer-layer simulated annealing algorithm model. New solutions are generated through Gaussian perturbation, and the fitness differences are calculated. A probabilistic acceptance strategy is used to jump out of the local optimum, and a solution that satisfies the spatio-temporal schedulable boundary conditions and is jointly optimal for the dual objectives is output, including the parameter set of mobile resource actions, the state set flag sequence, and the priority of the feasible scheduling path set. Specifically, after the inner-layer dual-model algorithm is iterated, the currently found optimal scheduling scheme is obtained, and the currently found optimal scheduling scheme is output and used as the initial solution for the outer-layer simulated annealing optimization stage to further optimize and adjust the scheduling scheme. In the simulated annealing optimization stage, new solutions for the parameter set of mobile resource actions, the state set flag sequence, and the priority of the feasible scheduling path set that are jointly optimal for the dual objectives are generated by perturbing the initial solution, and the fitness differences are calculated. If the difference is less than 0, the new solution is accepted. If the difference is greater than 0, the acceptance probability of the new solution is calculated, and the new solution is selectively accepted according to the acceptance probability. When the simulated annealing algorithm layer reaches the maximum number of iterations at different temperatures according to the cooling coefficient, the multi-layer collaborative optimization is terminated.

[0084] Reference Figure 10 , which shows the trend of system performance changing over time after adopting the above scheme, and the impact of different strategies adopted in different stages on system performance. Under ideal conditions, the energy supply level of the system remains at , when extreme weather or natural disasters occur at a certain moment, the equipment of the energy system is affected, and the system energy supply level begins to decline until the end time of the extreme weather or natural disaster, the system energy supply level drops to . And at the end time of the extreme weather or natural disaster, the system recovery strategy begins to be implemented. At time , the system energy supply begins to recover until time when the system performance is fully restored. After considering the energy system capacity configuration optimization model provided in the embodiments of the present invention, the ability of the energy system to overcome extreme weather and natural disasters is improved. When extreme weather or natural disasters occur, the system energy supply level drops to , and the system performance is improved by compared with that before the capacity configuration optimization. After considering the mobile resource scheduling model provided in the embodiments of the present invention, such as electric vehicles, the recovery time can be advanced from to , quickly restoring the load supply.

[0085] Additionally, an energy system resilience improvement system considering mobile resource load restoration under extreme weather provided by an embodiment of the present invention includes: a database construction module, configured to link the equipment failure probability and line failure probability under extreme disaster events obtained through dual-effect Markov and Monte Carlo processing based on multi-source extreme disaster event data, and couple the load change parameters driven by meteorology to construct a disaster response database. A configuration generation module, configured to convert the fluctuations in renewable energy output and the deviation of user load demand into fuzzy variables through fuzzy chance constraints, establish a capacity configuration optimization model in combination with the disaster response database, and generate a multi-energy complementary capacity configuration with the minimization of the levelized cost of electricity as the optimization objective. A configuration improvement module, configured to introduce mobile resources based on the multi-energy complementary capacity configuration to establish a resilience improvement scheduling architecture, and generate a resilience improvement plan adaptable to extreme disaster events by implementing the damaged road network topology reconstruction, spatio-temporal transfer chain modeling, and action state collaborative regulation of mobile resources.

[0086] Moreover, an automatic annotation device for high-quality drainage pipeline image data based on a deep learning model with multi-feature fusion provided by an embodiment of the present invention includes: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, and when the instructions are executed by the at least one processor, the at least one processor is enabled to execute the method for optimizing the energy system capacity and decoupling and coordinating the spatio-temporal scheduling of mobile resources under extreme disaster disturbances as described above.

[0087] Furthermore, an embodiment of the present invention provides a computer-readable storage medium, on which computer-executable instructions are stored, and when the executable instructions are executed by a processor, the method for optimizing the energy system capacity and decoupling and coordinating the spatio-temporal scheduling of mobile resources under extreme disaster disturbances as described above is implemented.

[0088] In summary, an embodiment of the present invention provides a method, system, device, and medium for optimizing the energy system capacity and decoupling and coordinating the spatio-temporal scheduling of mobile resources under extreme disaster disturbances, aiming to implement a three-level linkage architecture of the energy system under disaster scenarios, including disaster data-driven modeling, multi-energy capacity collaborative optimization, and dynamic scheduling of mobile resources. Refer to Figure 11 , and its core logic is divided into the following steps: Step 1: Obtain the historical extreme weather data of multiple different meteorological station observation points. The extreme weather and natural disasters collected mainly include earthquakes, typhoons, extreme high-temperature weather, and extreme low-temperature weather. And establish an extreme weather database based on the historical extreme weather data; Step 2: Simulate the device state transition probability through a two-state resilience adaptive Markov process, combine Monte Carlo random sampling to generate device cascading failure scenarios, construct a weather disaster energy system response database including device output limitation, load dynamic disturbance, and energy supply-demand imbalance, and generate a disaster response database based on the extreme weather database; Step 3: Define the uncertainty of wind-solar power generation as a fuzzy variable using fuzzy chance constraints and establish a fuzzy chance constraint model, and generate an energy system configuration optimization model under extreme weather through the disaster response database. With the goal of minimizing the levelized electricity cost, call solvers such as Gurobi to optimize the device capacity of the energy system, determine the optimal energy system capacity configuration, and transfer the optimized results to the dispatching layer. The optimized energy system can effectively reduce the interruption of energy supply under extreme events and improve the resilience and reliability of the system; Step 4: To further enhance the resilience of the energy system, on the basis of capacity configuration, couple mobile resources such as mobile emergency power vehicles, mobile diesel generators, and electric vehicles. With the goal of jointly minimizing the total dispatching cost, load shedding amount, and path energy consumption, adopt a multi-layer collaborative optimization framework to dynamically solve the power distribution, charging and discharging time sequence, and path priority, and generate an adaptive resilience plan for extreme disaster events, thereby dispatching mobile resources to supply power to restore load supply at the fastest speed.

[0089] The present invention realizes accurate risk quantification through the drive of a disaster database, and the fuzzy constraint model ensures the economic-resilience balance of multi-energy capacity coordination, and finally breaks through the upper limit of the response ability of fixed devices by means of dynamic scheduling of mobile resources.

[0090] Since the system / device described in the above embodiments of the present invention is the system / device adopted for implementing the method in the above embodiments of the present invention, based on the method described in the above embodiments of the present invention, those skilled in the art can understand the specific structure and deformation of the system / device, so it will not be repeated here. Any system / device adopted by the method in the above embodiments of the present invention belongs to the scope of protection of the present invention.

[0091] The present invention is described with reference to the flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowcharts and / or block diagrams can be implemented by computer program instructions, as well as the combination of processes and / or blocks in the flowcharts and / or block diagrams.

[0092] Although the preferred embodiments of the present invention have been described, those skilled in the art can make additional changes and modifications after learning the basic creative concepts. Therefore, the technical solutions should be interpreted as including the preferred embodiments and all changes and modifications falling within the scope of the present invention.

[0093] Obviously, those skilled in the art can make various modifications and variations to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the technical solutions of the present invention and their equivalent technologies, the present invention should also include these modifications and variations.

Claims

1. A method for energy system capacity optimization and mobile resource spatiotemporal decoupling coordinated scheduling under extreme disaster disturbances, characterized in that: include: Based on the acquired extreme disaster event data, the equipment failure probability and line failure probability obtained through two-state resilient adaptive Markov and Monte Carlo processing are linked, and the load change parameters driven by extreme disaster events are coupled to form a disaster response database; The output fluctuation of renewable energy and the deviation of user load demand are converted into fuzzy variables through fuzzy chance constraints. A capacity configuration optimization model is established in combination with the disaster response database. The optimization goal is to minimize the levelized cost of electricity and generate multi-energy complementary capacity configuration. Based on the multi-energy complementary capacity configuration, mobile resources are introduced to establish a flexible scheduling architecture. By implementing the damaged road network topology reconstruction, space-time transfer chain modeling and coordinated control of action states of mobile resources, an adaptive resilience enhancement plan for extreme disaster events is generated.

2. The method for energy system capacity optimization and mobile resource spatiotemporal decoupling coordinated scheduling under extreme disaster disturbances as claimed in claim 1 is characterized in that: According to the acquired extreme disaster event data, the equipment failure probability and line failure probability obtained by two-state resilient adaptive Markov and Monte Carlo processing are linked, and the load change parameters driven by extreme disaster events are coupled to form a disaster response database including: A multi-dimensional extreme disaster event database is constructed based on the extreme disaster event data obtained from multiple observation points; wherein the extreme disaster event data is a collection of disaster events that exceed the preset extreme threshold, including earthquakes, storms, extremely high temperature droughts, and extremely low temperature ice disasters; The two-state toughness adaptive Markov switching model is used to analyze the state transition chain of the equipment, and the equipment failure probability is generated through the dynamic game of failure rate and maintenance rate. The Monte Carlo method is used to perform multi-scenario sampling simulation to generate a time-series equipment failure state set. According to the spatial distribution intensity and time accumulation parameters of current extreme disaster events, the dynamic diffusion process of extreme disaster events is simulated through parameter spatial gradient mapping and time accumulation effect analysis, and a set of line failure probabilities including energy network lines and pipelines is generated; A dynamic load growth evolution model is constructed to quantify the dynamic load increment of extreme low temperature events or extreme low temperature events by introducing humidity correction terms, wind speed correction terms, and radiation correction terms, thus forming a time-series load dynamic parameter set; The equipment fault state set, line failure probability set and load dynamic parameter set are aligned in time and space, and combined with the multi-dimensional extreme disaster event database to generate a disaster response database.

3. The method for energy system capacity optimization and mobile resource spatiotemporal decoupling coordinated scheduling under extreme disaster disturbances as claimed in claim 2 is characterized in that: The two-state resilient adaptive Markov switching model is: ; ; ; ; ; In the formula, For equipment j The probability of normal operation, For equipment j The probability of failure, 1 means the device j The running status is normal, 0 means the device j Failure in operation, , is the time-varying maintenance rate; for the power grid and gas grid, , For the number of emergency repair teams, is the baseline maintenance rate when there are no additional repair resources. is the base number for emergency repairs; for energy systems, , and They are respectively the normal operation time of time-varying equipment and the time-varying repair time of equipment. is a randomly generated number from a uniform distribution [0, 1]. is a randomly generated number from a uniform distribution [0, 1]. For equipment j exist t The state of the moment, is the failure rate, , K ={earthquake, storm, extreme heat and drought, extreme cold and ice disaster}, K is a set of disaster types, For disasters K Real-time strength indicator, For disasters K The weight coefficient satisfies w k =1, For disasters K The nonlinear influence function is used to standardize the physical dimensions of different disasters. To quantify the amount of maintenance resources invested Disaster intensity The maintenance and disaster game intensity of the suppression capacity, R ( t )≥ D ( t )hour, γ ( t )≥0.5, indicating that maintenance is dominant and the probability of normal state of equipment is improved; R ( t )< D ( t )hour, γ ( t )<0.5, indicating disaster suppression and increased failure risk, θ is the benchmark threshold of resource disaster ratio; The dynamic load growth evolution model is: ; In the formula, is the predicted load when extreme weather occurs at time t, for t The predicted load when there is no extreme weather at all times, Caused by extreme weather t Load increment at any moment and body temperature Related, is a neural network load forecasting model, express t The actual temperature at the moment, express t The relative humidity correction term of the air at that moment, represents the wind speed correction term, Represents the solar radiation correction term.

4. The method for energy system capacity optimization and mobile resource spatiotemporal decoupling coordinated scheduling under extreme disaster disturbances as claimed in claim 1 is characterized in that: The output fluctuation of renewable energy and the deviation of user load demand are converted into fuzzy variables through fuzzy opportunity constraints. The capacity configuration optimization model is established in combination with the disaster response database. The optimization goal is to minimize the levelized cost of electricity. The multi-energy complementary capacity configuration is generated, including: The membership function is constructed by using fuzzy chance constraints to convert the random fluctuation of wind and solar power generation and the deviation of user load demand into a set of fuzzy variables under credible probability constraints. In response to regional energy demand, a system topology including energy supply equipment, energy storage equipment and energy conversion units is constructed based on a disaster response database, and capacity configuration elements including equipment output sets, equipment status sets and efficiency parameter sets are configured; Establish a capacity configuration optimization model with the goal of minimizing the levelized cost of electricity, generate model constraints based on capacity configuration factors, and embed the fuzzy variable set into the model constraints to generate a fuzzy feasible solution space that integrates uncertainty; The capacity configuration optimization model is iteratively solved in the fuzzy feasible solution space, and the multi-energy complementary capacity configuration that is compatible with the system topology is output.

5. The method for energy system capacity optimization and mobile resource spatiotemporal decoupling coordinated scheduling under extreme disaster disturbances as claimed in claim 4 is characterized in that: The equipment output is: ; In the formula, for t Electricity price at any time, For equipment i exist t The power generation at the time, For equipment i exist t The cooling capacity at the time, For equipment i exist t The heating amount at the time, Electricity-hydrogen conversion equipment i exist t The amount of hydrogen produced at the time, They are t The building electricity load demand, heating load demand and cooling load demand at the moment, for t Momentary energy storage equipment i Capacity; The device state set is represented as: ; In the formula, Indicates the device i exist t When the running state When it is 1, it means the device i It is a normal operating state. When it is 0, it means the device i It is a running fault state. Energy storage devices i exist t Always in charging and discharging state, when energy storage equipment i In a fault state, or When it is 0, it is in normal operation. or is 1, and ; The capacity configuration optimization model is: ; ； ； ； ； In the formula, TNPC is the total net present value, CRF is the capital recovery factor, for t The electrical load at any time, for t The charging load of electric vehicles at any moment, for t The amount of electricity load reduced at any time, is the annual investment cost of the system, is the discount rate, For equipment life, For equipment i The unit investment cost, For equipment i The installed capacity of is the annual operation and maintenance cost of the system, For equipment i Unit operation and maintenance cost, For equipment i operating power.

6. The method for energy system capacity optimization and mobile resource spatiotemporal decoupling coordinated scheduling under extreme disaster disturbances as claimed in claim 1, characterized in that: Based on the configuration of multi-energy complementary capacity, mobile resources are introduced to establish a flexible scheduling framework. By implementing the damaged road network topology reconstruction, spatiotemporal transfer chain modeling and action state coordinated regulation of mobile resources, an adaptive elastic improvement plan for extreme disaster events is generated, including: On the basis of multi-energy complementary capacity configuration, mobile resources including mobile emergency power supply vehicles, electric vehicles and diesel generators are introduced as elastic supplementary energy supply units to establish an elastic scheduling framework that coordinates fixed equipment and mobile resources; According to the disaster damage information obtained, the damaged road network nodes and edges are removed, the road network topology is reconstructed, and a feasible scheduling path set is generated to minimize the distance of mobile resources from the starting point to the target energy system; Parse the departure timestamp, estimated arrival timestamp, travel time interval and charging stop period of the mobile resource to generate time chain features, and extract the starting position coordinates, target energy supply station coordinates and path node coordinate sets of the mobile resource to generate space chain features; Verify the time chain characteristics and space chain characteristics, and calculate the maximum service radius, shortest response time and continuous power supply duration threshold of mobile resources based on the verified time and space characteristics combination to form the time and space dispatchable boundary; An action set representing the energy supply behavior parameters of the mobile resource is obtained according to the energy supply characteristics of the mobile resource, and a state set representing the operation state of the mobile resource is constructed according to the mutual exclusion constraint of the device operation; In the elasticity enhancement scheduling architecture, with the goal of minimizing the total cost of mobile resource scheduling, the lowest load reduction and the smallest path energy consumption, based on the spatiotemporal schedulable boundaries, a multi-layer collaborative optimization algorithm is used to dynamically solve the optimal solution for the parameters of the action set of mobile resources, the flag sequence of the state set and the priority of the feasible scheduling path set, to generate an adaptive elasticity enhancement plan for extreme disaster events.

7. The method for energy system capacity optimization and mobile resource spatiotemporal decoupling coordinated scheduling under extreme disaster disturbances as claimed in claim 6 is characterized in that: In the elasticity enhancement scheduling architecture, with the goal of minimizing the total cost of mobile resource scheduling, the lowest load reduction, and the minimum path energy consumption, based on the spatiotemporal schedulable boundary, a multi-layer collaborative optimization algorithm is used to dynamically solve the optimal solution of the action set parameters of the mobile resource, the flag sequence of the state set, and the priority of the feasible scheduling path set, and generate an adaptive elasticity enhancement plan for extreme disaster events, including: The action set parameters and state set flag sequences of mobile resources are encoded into the position vectors of honey badgers, and the initial honey badger population that meets the spatiotemporal schedulable boundaries is generated. The priority of the feasible scheduling path set is encoded into the migration path chain of individual bird flocks, and the initial bird flock population that meets the spatiotemporal scheduling boundary is generated and loaded into the path memory library; An inner-outer layer collaborative optimization framework is constructed. In the inner layer, the upper honey badger population performs dual-objective optimization of global scheduling cost and load reduction. The lower bird population optimizes the priority of feasible scheduling path sets through the path memory library sharing mechanism with the goal of minimizing path energy consumption. In the outer layer, the simulated annealing algorithm model is used to receive the optimal solution output by the inner layer as the initial solution to perform disturbance optimization. For each honey badger individual, the minimum energy consumption sequence of the corresponding path chain is matched from the initial bird population, and the total cost of mobile resource scheduling, load reduction and path energy consumption index are calculated. The normalized weights are aggregated into comprehensive performance indicators to generate a dual-driven scale of prey odor concentration-migration gravity; The dual populations are sorted according to the dual-drive scale, the current global optimal solution is selected, and combined with the introduced dynamic density factor, the honey badgers and bird individuals are controlled to perform two-stage optimization in the neighborhood of the current global optimal solution; In the wide-area mining stage, the search step length is expanded in the upper honey badger population and the concentration of prey odor, and the combination of action set parameters and state set flag sequences is traversed. In the lower bird population, the topological weights in the path memory are updated according to the path requirements of individual honey badgers to generate a wide-area feasible path set; In the fine mining stage, the search step size is shrunk in the upper honey badger population to lock the local optimum of the combination of action set parameters and state set flag sequence, and the migration gravity-driven priority sorting is performed on the wide-area feasible path set in the lower bird population; Constraint checks are performed on the updated honey badgers and bird flocks. If the updated honey badgers and bird flocks violate the spatiotemporal schedulable boundaries, topological rerouting repair is performed, the node order and pass weight of the path chain are adjusted, and the bird flock gravity calibration is applied. Priorities are reallocated based on the path memory library, and the repaired individuals are fed back to the corresponding populations to update the distribution of honey badgers and bird flocks. After the iteration of the inner architecture is completed, the current optimal scheduling plan is input into the outer simulated annealing algorithm model, and a new solution is generated through Gaussian perturbation and the fitness difference is calculated. The local optimum is jumped out with a probabilistic acceptance strategy, and the mobile resource action set parameters, state set flag sequence and feasible scheduling path set priority plan that meet the spatiotemporal schedulable boundary conditions and are jointly optimal for multiple objectives are output.

8. A system for optimizing energy system capacity and decoupling mobile resources in time and space under extreme disaster disturbances, characterized in that: include: A database construction module is used to link the equipment failure probability and line failure probability obtained through two-state resilience adaptive Markov and Monte Carlo processing based on the acquired extreme disaster event data, and couple the load change parameters driven by extreme disaster events to form a disaster response database; The configuration generation module is used to convert the output fluctuation of renewable energy and the deviation of user load demand into fuzzy variables through fuzzy opportunity constraints, establish a capacity configuration optimization model in combination with the disaster response database, and generate multi-energy complementary capacity configuration with the minimization of the levelized cost of electricity as the optimization goal; The configuration enhancement module is used to introduce mobile resources to establish an elastic enhancement scheduling architecture based on multi-energy complementary capacity configuration. By implementing the damaged road network topology reconstruction, space-time transfer chain modeling and coordinated control of action states of mobile resources, an adaptive elastic enhancement plan for extreme disaster events is generated.

9. A high-quality drainage pipe image data automatic annotation device based on a deep learning model with multi-feature fusion, characterized in that: include: at least one processor; And a memory communicatively connected to at least one processor; wherein the memory stores instructions executable by at least one processor, and the instructions are executed by at least one processor so that at least one processor can execute the method for energy system capacity optimization and spatiotemporal decoupling and coordinated scheduling of mobile resources under extreme disaster disturbances as described in any one of claims 1-7.

10. A computer-readable storage medium having computer-executable instructions stored thereon, characterized in that: When the executable instructions are executed by the processor, the method for energy system capacity optimization and mobile resource spatiotemporal decoupling and coordinated scheduling under extreme disaster disturbances as described in any one of claims 1 to 7 is implemented.

Images