Climate-driven hybrid energy storage comprehensive energy optimization method, system and equipment and storage medium

Abstract

The invention discloses a climate-driven hybrid energy storage integrated energy optimization method, system and device and a storage medium, and is applied to the technical field of integrated energy system planning and optimization, and the method comprises the steps: generating a year-round wind and light hourly scene covering seasonal features based on a Wasserstein generative adversarial network fusing gradient penalty and small-batch discrimination terms; the method comprises the following steps of: constructing a hybrid energy storage integrated energy system capacity-scheduling joint optimization model by taking an annual wind and light hourly scene as random input, and solving by adopting a Monte Carlo decomposition strategy to obtain optimal capacity configuration and an hourly scheduling strategy of each scene under the optimal capacity configuration; and performing energy storage duration integral evaluation on the hourly scheduling strategy to obtain utilization distribution of battery energy storage and hydrogen energy storage in the hybrid energy storage integrated energy system under different time scales. According to the method, the robustness and the economical efficiency of the integrated energy system under the extreme climate condition are remarkably improved.

Description

Technical Field

[0001] This invention relates to the field of integrated energy system planning and optimization technology, and more specifically to a climate-driven hybrid energy storage integrated energy optimization method, system, equipment and storage medium. Background Technology

[0002] In recent years, the proportion of renewable energy in the energy system has risen rapidly. Clean energy sources such as solar and wind power have significant intermittency and volatility, posing challenges to the safe and economical operation of integrated energy systems.

[0003] Integrated energy systems, as a novel energy utilization model, enable the interconnection and coordinated optimization of multiple energy forms such as electricity, heat, and hydrogen. By integrating heterogeneous energy units such as wind power, photovoltaic power generation, combined heat and power, and energy storage, integrated energy systems can improve energy utilization efficiency and reduce operating costs. At the same time, the coupled operation of multiple energy units increases the complexity of the model and places higher demands on the optimal scheduling of the system.

[0004] In integrated energy systems, not only does the load exhibit significant volatility and randomness, but the renewable energy generation side is also subject to strong uncertainties due to climate conditions. This dual uncertainty exposes the system to the risk of supply-demand imbalances at different time scales. For example, intraday load fluctuations and changes in solar irradiance may lead to short-term power shortages, while seasonal variations in wind power output can result in long-term energy imbalances.

[0005] To address this uncertainty, academic and engineering practices typically employ typical day methods, deterministic optimization, or scenario generation methods. However, typical day methods struggle to cover the diverse climate characteristics throughout the year, while deterministic methods may overlook extreme situations caused by fluctuations, leading to system planning results deviating from actual needs. Scenario generation methods improve the model's adaptability to uncertainty to some extent, but still suffer from problems such as unstable generated samples and excessive computational burden due to a large number of scenarios.

[0006] Energy storage technology is considered an important means to mitigate the fluctuations in renewable energy and improve system flexibility. Battery energy storage systems are characterized by fast response and high efficiency, making them suitable for hourly-level energy balancing and frequency regulation / peak shaving; hydrogen energy storage systems, on the other hand, have the potential for large-scale and long-term energy storage, enabling energy transfer on weekly, monthly, and even seasonal scales. However, how to rationally allocate and schedule short-term and long-term energy storage to complement each other's advantages is a key issue currently facing integrated energy system planning.

[0007] Furthermore, with the expansion of integrated energy systems and the dramatic increase in climate and operational data, the dimensionality and complexity of optimization models are constantly rising. Traditional mixed-integer programming methods suffer from significantly increased computational costs when dealing with high-dimensional annual time series and large-scale scenarios, and their solution speed and stability are insufficient to meet the needs of practical engineering applications.

[0008] Therefore, how to provide a climate-driven hybrid energy storage integrated energy optimization method, system, equipment and storage medium to improve model calculation efficiency and the rationality of energy storage system configuration while ensuring the realism and diversity of climate scenarios is an urgent problem to be solved by those skilled in the art. Summary of the Invention

[0009] In view of this, the present invention provides a climate-driven integrated energy optimization method, system, device, and storage medium for hybrid energy storage. By introducing climate-driven scene generation technology, it is possible to generate time-series data of wind and solar resources that are both realistic and diverse. On this basis, the Monte Carlo decomposition optimization method is used to solve the hourly capacity configuration and scheduling problem throughout the year, solving the problem of excessive computational burden of traditional optimization models in high-dimensional scenarios. At the same time, an energy storage duration integral index is proposed to quantify the synergistic effect of short-term battery energy storage and long-term hydrogen energy storage, thereby improving the reliability and economy of the integrated energy system at different time scales.

[0010] To achieve the above objectives, the present invention adopts the following technical solution: A climate-driven integrated energy optimization approach for hybrid energy storage includes: Step 1: Generate hourly scenes of scenery throughout the year, covering seasonal features, based on a Wasserstein generative adversarial network that integrates gradient penalty and mini-batch discriminant terms; Step 2: Using the annual wind and solar hourly scenarios as random input, construct a joint optimization model of capacity-schedule for a hybrid energy storage integrated energy system with annualized total cost as the objective, including the coupling of electricity, hydrogen energy, and thermal energy, as well as equipment operation constraints. Then, use the Monte Carlo decomposition strategy to solve the joint optimization model of capacity-schedule for the hybrid energy storage integrated energy system, and obtain the optimal capacity configuration and the hourly scheduling strategy for each scenario under the optimal capacity configuration. Step 3: Perform an integral evaluation of the energy storage duration on the hourly scheduling strategy to obtain the utilization distribution of battery energy storage and hydrogen energy storage in the hybrid energy storage integrated energy system at different time scales.

[0011] Optionally, in step 1, the Wasserstein generative adversarial network that integrates gradient penalty and mini-batch discriminant terms is as follows: Using generators With critic The adversarial structure; Mini-batch discrimination is introduced into the WGAN-GP architecture; among which, the generator... Using wind and light as two-channel inputs respectively, and upsampling through 1D convolution, downsampling, GLU residuals, and deconvolution, the output is a two-channel annual sequence; Critics It is a three-layer fully connected network, with the middle layer splicing mini-batch distance features; The loss function is as follows: The basic goal of WGAN-GP:

[0012]

[0013] in, , These are the base loss values ​​for the critic and the generator, respectively. The distribution of the actual climate data sample is denoted as . ; The distribution of the fake samples generated by the generator is denoted as . ; Represents the mathematical expectation; This represents the critic's rating output for the input sample; This is the gradient penalty coefficient, used to adjust the strength of the gradient penalty term; The random interpolation sampling points between the real sample and the generated sample are defined as follows: ,in for Random numbers between; This indicates that the critiquer output is relative to the input. The gradient; express Norm; Small batch discrimination criteria: Set up a critic The middle layer The features are Define the intra-batch similarity measure :

[0014] in, and These represent the first and second digits in the current batch. The and the first The feature vector of each sample in the intermediate layer of the critique; Indicates the Euclidean distance between features; By adding generator and critique losses, the following is formed:

[0015]

[0016] in, and These represent the total losses of the generator and the critic after adding mini-batch discriminant terms, respectively. This refers to the batch size during model training. and These are the weight coefficients for the mini-batch items in the generator and the critique, respectively, used to adjust the strength of the diversity constraint; This is an exponential function used to map distance to similarity weights.

[0017] Optionally, in step 2, the capacity-dispatch joint optimization model for the hybrid energy storage integrated energy system is as follows: Decision variables and parameters: Key capacity decision: Photovoltaic installed capacity installed capacity of wind turbines Electrolytic cell installed capacity Fuel cell installed capacity Battery installed capacity Hydrogen storage equivalent energy ; Hourly variables in a climate scenario: , , , , , ;in, This refers to the hourly power generation of photovoltaic systems. This refers to the hourly power generation of wind power. This refers to the power consumption of the electrolytic cell; Power generation for fuel cells; The charging power for battery energy storage; This refers to the battery's energy storage discharge power. For time period numbering, ; This refers to the hourly state of charge of the battery. This refers to the hourly equivalent state of charge of hydrogen energy storage. , These refer to the power purchased from the power grid and the power sold to the grid, respectively. Objective function: Minimize annualized total cost:

[0018]

[0019]

[0020]

[0021]

[0022] in, This represents the system's annualized investment cost; Annual equipment maintenance costs; This represents the annual net cost of electricity purchases. It is a collection of equipment, including photovoltaics, wind turbines, electrolyzers, fuel cells, batteries, and hydrogen storage containers; For equipment The installed capacity; For equipment The unit investment price per unit capacity; For equipment The capital recovery coefficient; The discount rate; For equipment Service life; For equipment exist Operating power during a given time period; For equipment The unit power operation and maintenance cost per unit; for Time-of-use electricity pricing for specific time periods; for The on-grid electricity price for a given time period; Electricity and hydrogen balance and equipment constraints: Power balance:

[0023] in, for System base electrical load demand during a given period; Hydrogen energy balance, in terms of equivalent energy:

[0024]

[0025]

[0026]

[0027]

[0028] in, and They are respectively Hydrogen charging and discharging power of hydrogen storage during specific time periods; This represents the self-loss coefficient for hydrogen energy storage. and These represent the efficiencies of the hydrogen storage and release processes, respectively. and These are the maximum rated power limits for hydrogen storage and hydrogen release, respectively. and for The binary variables represent the hydrogen storage and hydrogen release states respectively, and are used to ensure that the two are mutually exclusive. Constraints of Electrolyzers and Fuel Cells:

[0029]

[0030]

[0031]

[0032] in, The efficiency of hydrogen production in the electrolyzer; For fuel cell power generation efficiency; This is the lower heating value of hydrogen. Constraints of lithium battery energy storage systems:

[0033]

[0034]

[0035]

[0036]

[0037] in, This is the battery self-discharge rate coefficient; and These refer to the battery charging and discharging efficiency, respectively. and These are the maximum charging and discharging power limits for the battery, respectively. and for A binary variable used to ensure that charging and discharging are mutually exclusive; Fan output constraints:

[0038] in, for Real-time wind speed during the specified time period; To cut in wind speed; To cut off the wind speed; Rated wind speed; This refers to the rated power of the fan. The overall efficiency of wind turbine power generation; Photovoltaic output constraints:

[0039] in, Light intensity under standard test conditions; for The actual solar radiation intensity during the time period; For the overall efficiency of photovoltaic systems; Power grid interconnection mutual exclusion constraints:

[0040]

[0041]

[0042] in: and These refer to the power purchased from the power grid and the power sold to the grid, respectively. and These are the maximum permissible power sales and power purchase capacity at the grid connection point, respectively. and for Binary variables ensure that electricity sales and electricity purchases do not occur simultaneously.

[0043] Optionally, in step 2, the Monte Carlo decomposition strategy is used to solve the capacity-scheduling joint optimization model of the hybrid energy storage integrated energy system to obtain the optimal capacity configuration and the hourly scheduling strategy for each scenario under the optimal capacity configuration, specifically: The generated hourly scenery scenes throughout the year are sampled multiple times to form multiple sub-models of the year's duration and solved independently. The capacity configuration distribution is then aggregated, and the distribution center is selected as the optimal capacity configuration based on the solution set distribution. Given a fixed optimal capacity configuration, solve for the hourly scheduling strategy for each scenario under the optimal capacity configuration.

[0044] Optionally, in step 3, the time-by-time scheduling strategy is evaluated by energy storage duration integration to obtain the utilization distribution of battery energy storage and hydrogen energy storage in the hybrid energy storage integrated energy system at different time scales, specifically: Based on the time-by-time scheduling strategy, the time-by-time energy curve of energy storage is generated; The hourly energy curve of energy storage is discretized along the vertical axis and scanned row by row. The product of energy storage and duration in each segment is recorded as the energy storage duration integral index. The energy storage duration distribution map is obtained by accumulating the data for the whole year and all scenarios.

[0045] This invention also provides a climate-driven hybrid energy storage integrated energy optimization system utilizing a climate-driven hybrid energy storage integrated energy optimization method, comprising: The year-round landscape scene generation module is used to generate year-round landscape scenes covering seasonal features based on a Wasserstein generative adversarial network that integrates gradient penalty and mini-batch discriminant terms. The hourly scheduling strategy solution module is used to construct a capacity-scheduling joint optimization model of a hybrid energy storage integrated energy system with annualized total cost as the objective, including the coupling of electricity, hydrogen energy, and thermal energy, as well as equipment operation constraints, using hourly wind and solar scenarios throughout the year as random input. The Monte Carlo decomposition strategy is then used to solve the capacity-scheduling joint optimization model of the hybrid energy storage integrated energy system to obtain the optimal capacity configuration and the hourly scheduling strategy for each scenario under the optimal capacity configuration. Energy storage duration integral evaluation module: used to evaluate the energy storage duration integral of the hourly scheduling strategy, and obtain the utilization distribution of battery energy storage and hydrogen energy storage in the hybrid energy storage integrated energy system at different time scales.

[0046] The present invention also provides an electronic device, comprising: Memory, used to store computer programs; A processor is used to implement a climate-driven hybrid energy storage integrated energy optimization method when executing computer programs.

[0047] The present invention also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of a climate-driven hybrid energy storage integrated energy optimization method.

[0048] As can be seen from the above technical solutions, compared with the prior art, this invention discloses a climate-driven hybrid energy storage integrated energy optimization method, system, device, and storage medium. By using an improved generative adversarial network to generate climate-driven wind and solar scenarios, it can maintain temporal relevance and seasonal characteristics, avoiding the problems of insufficient samples and model collapse in traditional scenario generation methods. Through Monte Carlo decomposition optimization, the computational complexity of hourly multi-scenario models throughout the year is effectively reduced, ensuring the robustness of the optimization results. Through the energy storage duration integral index, the division of labor and synergistic effects of battery energy storage and hydrogen energy storage at different time scales can be intuitively revealed, providing a quantitative basis for the configuration ratio of short-term regulation and long-term energy transfer, thereby improving the economy, flexibility, and reliability of the integrated energy system. Attached Figure Description

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

[0050] Figure 1 This is a schematic diagram of the method flow provided by the present invention.

[0051] Figure 2 This is a schematic diagram of the Wasserstein generative adversarial network structure that integrates gradient penalty and mini-batch discriminant terms provided by the present invention.

[0052] Figure 3 This is a schematic diagram of the hybrid energy storage integrated energy system provided by the present invention.

[0053] Figure 4 This is a schematic diagram of the ESDI evaluation process provided by the present invention. Detailed Implementation

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

[0055] Example 1: Embodiment 1 of this invention discloses a climate-driven hybrid energy storage integrated energy optimization method, such as... Figure 1 As shown, it includes: Step 1: Based on the Wasserstein generative adversarial network that integrates gradient penalty and mini-batch discriminant terms, train it for spring, summer, autumn and winter respectively to learn seasonal features and generate hourly wind and solar scenes covering seasonal features throughout the year (annual renewable energy power generation scene).

[0056] A Wasserstein generative adversarial network (WGAN-GP-minibatch) integrating gradient penalty and mini-batch discriminant terms is used to improve the diversity and stability of climate samples, such as... Figure 2 As shown, specifically: Using generators With critic The adversarial structure; Mini-batch discrimination is introduced into the WGAN-GP architecture to encourage intra-batch sample diversity; among which, the generator... As shown in Table 1, wind and light are used as 2-channel inputs respectively, and after 1D convolution, downsampling, GLU residual and deconvolution upsampling, a 2-channel annual sequence is output; Critics As shown in Table 2, it is a three-layer fully connected network with the middle layer splicing mini-batch distance features; Table 1 Generator structure

[0057] Table 2 Critics structure

[0058] The loss function is as follows: The basic goal of WGAN-GP:

[0059]

[0060] in, , These are the base loss values ​​for the critic and the generator, respectively. The distribution of the actual climate data sample is denoted as . ; The distribution of the fake samples generated by the generator is denoted as . ; Represents the mathematical expectation; This represents the critic's rating output for the input sample; This is the gradient penalty coefficient, used to adjust the strength of the gradient penalty term; The random interpolation sampling points between the real sample and the generated sample are defined as follows: ,in for Random numbers between; This indicates that the critiquer output is relative to the input. The gradient; express Norm (Euclidean norm); Small batch discrimination criteria: Set up a critic The middle layer The features are Define the intra-batch similarity measure :

[0061] in, and These represent the nth and nth elements in the current batch. The and the first The feature vector of each sample in the intermediate layer of the critique; Indicates the Euclidean distance between features; By adding generator and critique losses, the following is formed:

[0062]

[0063] in, and These represent the total losses of the generator and the critic after adding mini-batch discriminant terms, respectively. This refers to the batch size during model training. and These are the weight coefficients for the mini-batch items in the generator and the critique, respectively, used to adjust the strength of the diversity constraint; The exponential function is used to map distance to similarity weights. This design improves the critiquer's sensitivity to intra-batch diversity, prompting the generator... Generate more diverse sequences and mitigate pattern collapse.

[0064] It also includes comparing the generated wind and solar sequences with historical measurements using three types of indicators: autocorrelation (ACF), cumulative distribution function (CDF), and spectral density (SD), in order to assess the degree of matching between temporal correlation, distribution consistency, and frequency domain characteristics.

[0065] Step 2: Using the annual wind and solar hourly scenarios as random inputs, construct a joint optimization model of capacity-schedule for a hybrid energy storage integrated energy system (multi-scenario mixed integer linear programming model) with annualized total cost as the objective and including coupling of electricity, hydrogen energy, and thermal energy as well as equipment operation constraints. Then, use the Monte Carlo decomposition strategy to solve the joint optimization model of capacity-schedule for the hybrid energy storage integrated energy system to obtain the optimal capacity configuration and the hourly scheduling strategy for each scenario under the optimal capacity configuration.

[0066] Hybrid energy storage integrated energy system, such as Figure 3 As shown, it includes: Battery Energy Storage System (BESS) and Hydrogen Energy Storage System (P2H2P); wherein: BESS is mainly used for regulating intraday and daily-scale fluctuations; Hydrogen Energy Storage System is used for energy balance on weekly, monthly and seasonal scales.

[0067] Decision variables and parameters: Key capacity decision: Photovoltaic installed capacity installed capacity of wind turbines Electrolytic cell installed capacity Fuel cell installed capacity Battery installed capacity Hydrogen storage equivalent energy ; Hourly variables in a climate scenario: , , , , , ;in, This refers to the hourly power generation of photovoltaic systems. This refers to the hourly power generation of wind power. This refers to the power consumption of the electrolytic cell; Power generation for fuel cells; The charging power for battery energy storage; This refers to the battery's energy storage discharge power. For time period numbering, ; This refers to the hourly state of charge of the battery. This refers to the hourly equivalent state of charge of hydrogen energy storage. , These refer to the power purchased from the power grid and the power sold to the grid, respectively. Objective function: Minimizing annualized total cost (TAC) includes equipment investment costs, maintenance costs, and grid interaction costs:

[0068]

[0069]

[0070]

[0071]

[0072] in, This represents the system's annualized investment cost; Annual equipment maintenance costs; This represents the annual net cost of electricity purchases. It is a collection of equipment, including photovoltaics, wind turbines, electrolyzers, fuel cells, batteries, and hydrogen storage containers; For equipment The installed capacity; For equipment The unit investment price per unit capacity; For equipment The capital recovery coefficient; The discount rate; For equipment Service life; For equipment exist Operating power during a given time period; For equipment The unit power operation and maintenance cost per unit; for Time-of-use electricity pricing for specific time periods; for The on-grid electricity price for a given time period; Electricity and hydrogen balance and equipment constraints: Power balance:

[0073] in, for System base electrical load demand during a given period; Hydrogen energy balance, in terms of equivalent energy:

[0074]

[0075]

[0076]

[0077]

[0078] in, and They are respectively Hydrogen charging and discharging power (equivalent to energy flow) of hydrogen energy storage during a given time period; This is the self-loss coefficient (storage efficiency) for hydrogen energy storage. and These represent the efficiencies of the hydrogen storage and release processes, respectively. and These are the maximum rated power limits for hydrogen storage and hydrogen release, respectively. and for The binary variables represent the hydrogen storage and hydrogen release states respectively, and are used to ensure that the two are mutually exclusive. Constraints of Electrolyzers and Fuel Cells:

[0079]

[0080]

[0081]

[0082] in, The efficiency of hydrogen production in the electrolyzer; For fuel cell power generation efficiency; This is the lower heating value of hydrogen. Constraints of lithium battery energy storage systems:

[0083]

[0084]

[0085]

[0086]

[0087] in, This is the battery self-discharge rate coefficient; and These refer to the battery charging and discharging efficiency, respectively. and These are the maximum charging and discharging power limits for the battery, respectively. and for A binary variable used to ensure that charging and discharging are mutually exclusive; Fan output constraints:

[0088] in, for Real-time wind speed during the specified time period; To cut in wind speed; To cut off the wind speed; Rated wind speed; This refers to the rated power of the fan. The overall efficiency of wind turbine power generation; Photovoltaic output constraints:

[0089] in, The light intensity under standard test conditions (usually taken as...) ); for The actual solar radiation intensity during the time period; For the overall efficiency of photovoltaic systems; Power grid interconnection mutual exclusion constraints:

[0090]

[0091]

[0092] in: and These refer to the power purchased from the power grid and the power sold to the grid, respectively. and These are the maximum permissible power sales and power purchase capacity at the grid connection point, respectively. and for Binary variables ensure that electricity sales and electricity purchases do not occur simultaneously.

[0093] The Monte Carlo decomposition strategy is used to solve the capacity-scheduling joint optimization model of the hybrid energy storage integrated energy system, obtaining the optimal capacity configuration and the hourly scheduling strategy for each scenario under the optimal capacity configuration, specifically: The generated hourly landscape scenes throughout the year (generator) Output) Perform multiple samplings (random sampling) The system is divided into several scenarios, forming multiple sub-models for the entire year (each scenario constitutes an 8760-hour MILP sub-model), which are solved independently and then aggregated to obtain the capacity configuration distribution. Based on the solution set distribution, the distribution center (geometric centroid) is selected as the optimal capacity configuration. ; Fixed optimal capacity configuration Solve for the time-by-time scheduling strategy for each scenario under optimal capacity configuration. .

[0094] Step 3: Perform an Energy Storage Duration Integral (ESDI) evaluation on the hourly scheduling strategy to obtain the utilization distribution of battery energy storage and hydrogen energy storage in the hybrid energy storage integrated energy system at different time scales.

[0095] An integral evaluation of the energy storage duration is performed on the hourly dispatch strategy to obtain the utilization distribution of battery energy storage and hydrogen energy storage in the hybrid energy storage integrated energy system at different time scales, such as... Figure 4 As shown, specifically: Based on the time-by-time scheduling strategy, the time-by-time energy curve of energy storage is generated; The hourly energy curve of energy storage is discretized along the vertical axis and scanned row by row. The product of energy storage and duration in each segment is recorded as the energy storage duration integral index. The energy storage duration distribution map is obtained by accumulating the data for the whole year and all scenarios.

[0096] Specifically, set in the scene The hourly energy sequence of the lower energy storage is (Recommended) Discretize the vertical axis (SOC value) into several layers. Perform a lateral connectivity scan on each energy level and record the duration of each segment. With the energy amplitude of this layer The product is denoted as:

[0097] The ESDI duration distribution was obtained by accumulating it throughout the year and across all scenarios:

[0098] This method essentially calculates the distribution of the "energy-time area," which can separate and visualize the contributions of "short-term regulation (battery energy storage)" and "long-term balance (hydrogen energy storage)." The implementation steps are: "discrete the vertical axis → scan line by line → record ESDI → statistics / plotting."

[0099] Example 2: Embodiment 2 of the present invention discloses a climate-driven hybrid energy storage integrated energy optimization system utilizing a climate-driven hybrid energy storage integrated energy optimization method, comprising: The year-round landscape scene generation module is used to generate year-round landscape scenes covering seasonal features based on a Wasserstein generative adversarial network that integrates gradient penalty and mini-batch discriminant terms. The hourly scheduling strategy solution module is used to construct a capacity-scheduling joint optimization model of a hybrid energy storage integrated energy system with annualized total cost as the objective, including the coupling of electricity, hydrogen energy, and thermal energy, as well as equipment operation constraints, using hourly wind and solar scenarios throughout the year as random input. The Monte Carlo decomposition strategy is then used to solve the capacity-scheduling joint optimization model of the hybrid energy storage integrated energy system to obtain the optimal capacity configuration and the hourly scheduling strategy for each scenario under the optimal capacity configuration. Energy storage duration integral evaluation module: used to evaluate the energy storage duration integral of the hourly scheduling strategy, and obtain the utilization distribution of battery energy storage and hydrogen energy storage in the hybrid energy storage integrated energy system at different time scales.

[0100] Example 3: Embodiment 3 of the present invention also discloses an electronic device, comprising: Memory, used to store computer programs; A processor is used to implement a climate-driven hybrid energy storage integrated energy optimization method when executing computer programs.

[0101] Example 4: Embodiment 4 of the present invention discloses a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of a climate-driven hybrid energy storage integrated energy optimization method.

[0102] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the apparatus disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the description is relatively simple; relevant parts can be referred to the method section.

[0103] The above description of the disclosed embodiments enables those skilled in the art to make or use the invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the invention. Therefore, the invention is not to be limited to the embodiments shown herein, but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A climate-driven hybrid energy storage integrated energy optimization method, characterized in that, Comprise: Step 1: a Wasserstein generative adversarial network based on fusion gradient penalty and small batch discriminant term is used to generate a full-year light hourly scene covering seasonal characteristics; Step 2: the full-year light hourly scene is used as a random input to build a mixed energy storage comprehensive energy system capacity-scheduling joint optimization model with annual total cost as the target, including power, hydrogen energy, heat energy coupling and equipment operation constraints, and a Monte Carlo decomposition strategy is used to solve the mixed energy storage comprehensive energy system capacity-scheduling joint optimization model to obtain the optimal capacity configuration and hourly scheduling strategy of each scene under the optimal capacity configuration; Step 3: the energy storage time length integral evaluation is performed on the hourly scheduling strategy to obtain the utilization distribution of battery energy storage and hydrogen energy storage in the mixed energy storage comprehensive energy system at different time scales.

2. The climate-driven hybrid energy storage integrated energy optimization method according to claim 1, wherein, In step 1, the Wasserstein generative adversarial network based on fusion gradient penalty and small batch discriminant term is as follows: Adopting a generator Adversarial structure with a critic ; Introducing small batch discrimination on the WGAN-GP architecture; wherein the generator Respectively take wind and light as 2 channel input, and output 2 channel annual sequence through 1D convolution, down sampling, GLU residual and deconvolution up sampling; the discriminator Is a three-layer fully connected network, and the middle layer splices small batch distance features; The loss function is as follows: The basic goal of WGAN-GP is: wherein, , are the base loss values of the critic and the generator, respectively; denotes the real climate data samples, whose distribution is denoted as ; denotes the fake samples generated by the generator, whose distribution is denoted as ; denotes the mathematical expectation; denotes the score output of the critic for the input sample; is the gradient penalty coefficient, used to adjust the intensity of the gradient penalty term; is a random interpolation sampling point between the real sample and the generated sample, defined as , wherein is a random number between ; denotes the gradient of the critic output with respect to the input ; denotes the norm; Mini-batch discriminative term: let the discriminator be the middle layer of the first feature , define the intra-batch similarity measure : wherein, and denote the feature vector of the i-th and j-th sample in the current batch at the intermediate layer of the critic, respectively; denote the feature vector of the i-th and j-th sample in the current batch at the intermediate layer of the critic, respectively; denote the feature vector of the i-th and j-th sample in the current batch at the intermediate layer of the critic, respectively; denote the Euclidean distance between features. The generator and critic loss are added to form: where, and are the total loss of generator and discriminator with small batch discrimination term added respectively; is the batch size during model training; and are the weight coefficients of small batch term in generator and discriminator respectively, used to adjust the strength of diversity constraint; is the exponential function, used to map the distance to the similarity weight.

3. The climate-driven hybrid energy storage integrated energy optimization method of claim 1, wherein, In step 2, the mixed energy storage comprehensive energy system capacity-scheduling joint optimization model is as follows: Decision variables and parameters: Main capacity decisions: PV installed capacity , wind turbine installed capacity , electrolyzer installed capacity , fuel cell installed capacity , battery installed capacity , hydrogen storage energy equivalent energy ; Hourly operating variables under climate scenarios: , , , , , ; wherein, is the photovoltaic hourly power production; is the wind power hourly power production; is the electrolyzer power consumption; is the fuel cell power production; is the battery storage charging power; is the battery storage discharging power; is the time period number, ; is the battery hourly state of charge; is the hydrogen storage hourly equivalent state of charge; , are the grid hourly power purchase and power sale, respectively; Objective function: Annual total cost minimization: wherein, is the annualized investment cost of the system; is the annual operation and maintenance cost of the equipment; is the annual net electricity purchase cost; is the set of equipment, including photovoltaic, wind turbine, electrolyzer, fuel cell, battery and hydrogen storage container; is the installed capacity of the equipment ; is the unit capacity investment unit price of the equipment ; is the capital recovery coefficient of the equipment ; is the discount rate; is the service life of the equipment ; is the operating power of the equipment in the time period; is the unit power operation and maintenance cost unit price of the equipment ; is the time-of-use electricity purchase price in the time period; is the grid-connected electricity price in the time period; Electricity, hydrogen balance and equipment constraints: Electricity balance: wherein, is system base electrical load demand for the period; Hydrogen energy balance in equivalent energy: wherein, and are respectively charging and discharging power of hydrogen storage during the period; is the self-loss coefficient of hydrogen storage; and are respectively the efficiencies of hydrogen storage and release processes; and are respectively the upper limits of maximum rated power of hydrogen storage and release; and are binary variables, respectively representing the hydrogen storage and release states, for ensuring mutual exclusivity. Electrolytic cell and fuel cell constraints: wherein, is the hydrogen production efficiency of the electrolyzer; is the power generation efficiency of the fuel cell; is the lower heating value of the hydrogen gas; Lithium battery energy storage system constraints: wherein, is the battery self-discharge rate coefficient; and are the battery charge and discharge efficiencies, respectively; and are the maximum charge and discharge power limits of the battery, respectively; and are binary variables to guarantee that charging and discharging are mutually exclusive. Fan output constraints: wherein, is real-time wind speed of the period; is the cut-in wind speed; is the cut-out wind speed; is the rated wind speed; is the fan rated power; is the fan power generation overall efficiency; Photovoltaic output constraints: wherein is the light intensity under standard test conditions; is is the actual solar radiation intensity for the period; is the overall efficiency of the photovoltaic system; Grid interconnection and exclusion constraints: wherein: and are the grid purchase and sell power, respectively; and are the maximum allowed sell and purchase power at the point of common coupling, respectively; and are binary variables that guarantee that sell and purchase actions do not occur simultaneously.

4. The climate-driven hybrid energy storage integrated energy optimization method of claim 1, wherein, In step 2, the Monte Carlo decomposition strategy is used to solve the mixed energy storage comprehensive energy system capacity-scheduling joint optimization model to obtain the optimal capacity configuration and hourly scheduling strategy of each scene under the optimal capacity configuration, which is as follows: The generated full-year light hourly scene is sampled multiple times to form multiple full-year length sub-models and is independently solved to obtain a capacity configuration distribution, and the distribution center is selected as the optimal capacity configuration based on the solution set distribution; The optimal capacity configuration is fixed, and the hourly scheduling strategy of each scene under the optimal capacity configuration is solved.

5. The climate-driven hybrid energy storage integrated energy optimization method of claim 1, wherein, In step 3, the energy storage time length integral evaluation is performed on the hourly scheduling strategy to obtain the utilization distribution of battery energy storage and hydrogen energy storage in the mixed energy storage comprehensive energy system at different time scales, which is as follows: Based on the hourly scheduling strategy, an energy storage hourly energy curve is generated; The energy storage hourly energy curve is discretized in the vertical axis and scanned row by row, and the product of each energy storage energy x duration is recorded as an energy storage time length integral indicator, and the energy storage time length distribution diagram is obtained by accumulating the full year and full scene.

6. A climate-driven hybrid energy storage integrated energy optimization system using the climate-driven hybrid energy storage integrated energy optimization method of any one of claims 1-5, characterized in that, Comprise: The full-year light hourly scene generation module is used to generate a full-year light hourly scene covering seasonal characteristics based on a Wasserstein generative adversarial network based on fusion gradient penalty and small batch discriminant term; a hourly scheduling strategy solving module, configured to take the hourly wind and solar scene as a random input, construct a capacity-scheduling joint optimization model of a hybrid energy storage comprehensive energy system, and solve the capacity-scheduling joint optimization model by using a Monte Carlo decomposition strategy, to obtain an optimal capacity configuration and a hourly scheduling strategy of each scene under the optimal capacity configuration; an energy storage time length integral evaluation module, configured to perform energy storage time length integral evaluation on the hourly scheduling strategy, to obtain utilization distribution of battery energy storage and hydrogen energy storage in the hybrid energy storage comprehensive energy system at different time scales.

7. An electronic device, comprising: The method comprises the following steps: a memory, configured to store a computer program; a processor, configured to execute the computer program to implement the steps of the climate-driven hybrid energy storage comprehensive energy optimization method according to any one of claims 1-5.

8. A computer-readable storage medium, characterized in that, The computer program is stored on the computer readable storage medium, and the computer program is executed by the processor to implement the steps of the climate-driven hybrid energy storage comprehensive energy optimization method according to any one of claims 1-5.

Images