Low-carbon optimal dynamic regulation control method for wind energy storage hydrogen energy system

Abstract

The application discloses a low-carbon optimal dynamic regulation control method for a wind energy storage hydrogen energy system, aiming at solving the problems of low prediction accuracy, single optimization target, incomplete constraint consideration and lack of dynamic feedback in the existing wind energy-storage-hydrogen energy system control. Through real-time collection of multi-dimensional operation data of the system, an LSTM and XGBoost fusion model is used to predict the wind power, power grid load and carbon emission coefficient in the next 24 hours; a double-target function of carbon emission and cost minimization is constructed, and component and power grid constraints are set; a Pareto optimal solution set is obtained through an NSGA-II algorithm, and optimal parameters are screened to generate control instructions; after the instructions are executed, actual data is collected and compared with the optimal parameters, and when the deviation exceeds the limit, re-optimization regulation is performed to form a closed loop. The method improves the prediction accuracy, realizes low-carbon and economic synergy, guarantees the stability of the system and helps energy transformation.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of wind energy storage hydrogen energy system optimal control, and particularly relates to a low-carbon optimal dynamic regulation control method for a wind energy storage hydrogen energy system. BACKGROUND

[0002] Under the background of global double carbon target promotion and energy structure transformation, the large-scale application of renewable energy has become the core path to realize energy low carbonization. Wind energy, as a mature technology and abundant resource, its development and utilization scale continues to expand. However, wind energy has significant intermittency, volatility and randomness, such as wind speed being affected by season, weather, time period and other factors, which fluctuates sharply, resulting in unstable wind power output. Direct grid connection will impact the frequency and voltage stability of the power grid, and even cause power grid scheduling difficulties, frequent wind curtailment and other problems, which seriously restricts the efficient use of wind energy resources.

[0003] To alleviate the impact of wind power volatility on the power grid, energy storage technology and hydrogen energy technology are widely introduced into the wind energy utilization system to build a wind energy-energy storage-hydrogen energy integrated system. Energy storage devices can smooth wind power fluctuations through charging and discharging, store electricity when there is excess wind power, and release electricity when there is a shortage of wind power. The hydrogen energy subsystem can convert excess wind power into hydrogen through electrolysis tanks to realize long-term energy storage and cross-scenario application of electricity, further improving the flexibility of system energy utilization. However, the integrated system has a complex structure, involving wind turbine generators, energy storage devices, electrolytic tanks, hydrogen storage tanks and other components. The operating characteristics of each component are quite different, and the system needs to interact with the external power grid, while also needing to consider carbon emission control and economic operation. Traditional control methods have been difficult to meet the requirements of efficient, low-carbon and economic operation of the system.

[0004] Traditional control methods mostly use single target optimization, such as only pursuing the lowest operating cost or the highest wind power consumption rate, without considering the total carbon emissions as a core optimization target, which is out of line with the current low-carbon development needs. In the prediction stage, a single model is often used to predict wind power and grid load demand. Due to the limitations of the model itself, the prediction accuracy is insufficient, especially in extreme weather or load mutation scenarios, the prediction deviation is large, which leads to a lack of reliable basis for subsequent control instruction development. In the optimization solving stage, traditional algorithms are prone to local optimal solutions when dealing with multi-objective and multi-constraint problems, making it difficult to generate globally optimal system operating parameters, and lacking dynamic feedback adjustment mechanisms. When the actual operating state of the system deviates from the optimal parameters beyond the allowed range, the control strategy cannot be adjusted in time, resulting in decreased system operating efficiency, excessive carbon emissions or increased operating costs.

[0005] In addition, the existing system lacks systematicness in parameter integration and intermediate variable definition, does not fully consider the quantitative calculation of system full life cycle carbon emissions and multi-type costs, and is difficult to realize accurate management and control of the system running state. For example, some control methods do not dynamically adjust the electrolytic cell carbon emission coefficient according to the energy structure, resulting in carbon emission calculation deviation; the cost accounting ignores the energy storage charging and discharging loss cost or the hydrogen storage and transportation cost, resulting in the disconnection between the economic optimization goal and the actual operation cost. These problems together cause the wind energy-storage-hydrogen energy system to be difficult to realize the coordinated optimization of low carbonization, economy and stability in actual operation, which restricts the large-scale promotion and application of such system.

[0006] Therefore, the present application provides a low-carbon optimal dynamic regulation control method for a wind energy-storage-hydrogen energy system to solve the above technical problems. SUMMARY

[0007] The present application aims to provide a low-carbon optimal dynamic regulation control method for a wind energy-storage-hydrogen energy system to solve the technical problems existing in the background art.

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

[0009] The low-carbon optimal dynamic regulation control method for a wind energy-storage-hydrogen energy system comprises the following steps:

[0010] S1: Real-time acquisition of system running data, the system running data including wind turbine data, energy storage device data, hydrogen energy subsystem data, real-time electricity price of the power grid, carbon emission coefficient and load demand data;

[0011] S2: Construction of a multi-objective prediction model, with historical wind speed data, historical power grid load data and historical carbon emission coefficient data as input features, and outputting wind power prediction results, power grid load demand prediction results and carbon emission coefficient prediction results in a specified future time period;

[0012] S3: Based on the wind power prediction results, power grid load demand prediction results and carbon emission coefficient prediction results, a low-carbon optimal objective function is constructed, and the constraint conditions of system operation are set;

[0013] S4: Based on the prediction results, combined with the constraint conditions, through operations such as population initialization, crossover, mutation, non-dominated sorting and crowding degree calculation, an optimal solution set satisfying the low-carbon and economic double objectives is obtained, and the optimal operation parameters are selected from the optimal solution set to generate corresponding dynamic regulation control instructions;

[0014] S5: The dynamic regulation control instructions are transmitted to the actuators of each component of the wind energy-storage-hydrogen energy system, and each actuator executes the corresponding operation according to the control instructions;

[0015] S6: Real-time acquisition system components perform operation after the actual running state data; the actual running state data is transmitted to the central control unit of the comparative analysis module, compared with the optimal running parameters, and the deviation value is calculated; if the deviation value is within the allowable error range, the current control instruction is maintained; if the deviation value exceeds the allowable error range, the solving parameters of the multi-objective optimization algorithm are re-adjusted, and a new dynamic adjustment control instruction is generated again, and steps S5 and S6 are repeated to realize closed-loop feedback adjustment of the system running state.

[0016] Preferably, the multi-objective prediction model of step S2 includes an LSTM model and an XGBoost model.

[0017] The LSTM model adopts a 3-layer network architecture, including an input layer, 2 hidden layers and an output layer, 64 neurons are set for each hidden layer, the ReLU activation function is selected, and a Dropout layer is added.

[0018] The number of decision trees of the XGBoost model is 100, the maximum tree depth is 6, the learning rate is 0.1, and the loss function adopts the root mean square error.

[0019] Preferably, the specific process of wind power prediction in step S2 is as follows:

[0020] S21: Obtain at least 1 year of historical wind speed data, and match the historical output power data of the wind turbine generator set at the corresponding period, extract the time sequence characteristics including the mean, maximum value, fluctuation amplitude and trend slope of the wind speed in each window by using the sliding window method, form a historical wind speed time sequence characteristic matrix, screen non-time sequence characteristics including date type and seasonal factors, and convert them into numerical features by using one-hot encoding;

[0021] S22: Input the historical wind speed time sequence characteristic matrix into the LSTM model for training, and input the real-time collected wind turbine generator set data into the LSTM model, and output the preliminary prediction result of the wind power within 24 hours in the future;

[0022] S23: Input the preliminary prediction result of the wind power as the basic feature into the XGBoost model together with the non-time sequence characteristics, and output the corrected prediction result of the wind power;

[0023] S24: Weighted fusion is performed on the preliminary prediction result of the wind power and the corrected prediction result of the wind power to obtain the final wind power prediction value.

[0024] Preferably, the specific process of constructing the low-carbon optimal objective function based on the wind power prediction result, the power grid load demand prediction result and the carbon emission coefficient prediction result in step S3 is as follows:

[0025] S31: Integrate the final prediction results of wind power generation, power grid load demand, and carbon emission coefficient by time dimension:

[0026] S32: Obtain inherent parameters of wind energy storage hydrogen energy system as the basis for setting target function calculation and constraint conditions, including cost type parameters, technology type parameters, and efficiency type parameters;

[0027] S33: Based on the prediction results and system parameters, define intermediate variables needed in the target function and constraint conditions, and establish calculation relationships;

[0028] S34: Based on the preprocessed prediction results and intermediate variables, construct a low-carbon optimal double-objective function with the dual objectives of minimizing total carbon emissions and minimizing operating costs in the entire life cycle of the system.

[0029] Preferably, the specific process of step S34 is as follows:

[0030] S341: Construct the first target function: minimize total carbon emissions of the system minf 1:

[0031] The calculation formula is based on the predicted carbon emission coefficient and the actual power consumption:

[0032] ;

[0033] Wherein: P grid-buy (t) is the active power purchased from the power grid in the tth period, i.e., the power purchased from the power grid when the system's own power generation is insufficient, P grid-buy (t)= max [ P grid (t),0], P grid (t) is the net interaction power between the system and the external power grid at time t, which is a core variable reflecting the direction and size of energy flow between the system and the power grid. When P grid (t)>0, it means that the system's own power generation (wind power, energy storage discharge, etc.) is less than the total demand of the load and the power consumption equipment (such as electrolytic cells), and it needs to purchase power from the power grid, so P grid-buy (t)= P grid (t).

[0034] When P grid (t)≤0, it means that the system's own power generation is greater than the total demand, and the excess power will be sold to the power grid (or abandoned), so there is no need to purchase power from the power grid, and therefore P grid-buy(t)=0; α (t) is the carbon emission coefficient corresponding to the power grid purchase in time period t, which reflects the carbon emission per unit of electricity generated by the power grid in that time period, and varies with the power grid energy structure and time period; P elec (t) represents the power generation of the system's self-contained power source (such as a diesel generator or fuel cell) during time period t. α (t) represents the carbon emission coefficient corresponding to the power grid's electricity purchase in time period t, reflecting the carbon emission per unit of electricity generated by the power grid in that time period, which varies with the power grid's energy structure and time period. β This is the carbon emission coefficient per unit power of the electrolyzer, in kgCO2 / kWh. If the electricity for the electrolyzer comes from wind power or energy storage, it can be adjusted according to the actual energy structure. It is assumed to be included in the total carbon emissions of the system here. The summation range is the next 24 hours, that is, it covers the carbon emissions of the entire forecast period.

[0035] S342: Constructing the second objective function: minimizing the system operating cost. minf 2:

[0036] Operating costs encompass four categories: energy storage charging and discharging losses, electrolyzer energy consumption, hydrogen storage and transportation, and electricity purchase from the grid. All costs are quantified based on forecast results and intermediate variables.

[0037] ;

[0038] Energy storage charging and discharging loss costs: P bat-loss (t)× δ Based on real-time calculations of power loss and unit loss cost, P bat-loss (t) is the th t The power loss generated during the charging and discharging process of the time-limited energy storage device, and the energy loss caused by charging and discharging efficiency. δ This is the cost coefficient corresponding to the unit power loss of energy storage, reflecting the economic cost of energy storage losses.

[0039] Energy consumption cost of electrolytic cell operation: The cost is directly linked to the operating power of the electrolytic cell; the higher the power, the higher the cost. P elec (t) represents the operating power of the electrolytic cell in time period t. This is the energy consumption cost coefficient corresponding to the unit operating power of the electrolytic cell;

[0040] Hydrogen storage and transportation costs: Q H2 (t)× ζ It is produced only when hydrogen production exceeds consumption, and is directly proportional to the amount stored and transported. Q H2(t) represents the net increase in hydrogen that needs to be stored or transported in time period t, which is the production minus the consumption. ζ This is the storage and transportation cost coefficient per unit mass of hydrogen.

[0041] Electricity purchase cost from the grid: P grid-buy (t)×γ(t), calculated based on the predicted power purchase and the corresponding electricity price for the time period. P grid-buy γ(t) represents the active power purchased from the grid in time period t, and γ(t) represents the unit price of electricity purchased from the grid in time period t, which fluctuates with time.

[0042] S343: Integration of Dual Objective Functions: Integrating two single objective functions into a unified low-carbon optimal objective function, represented in vector form, clearly defining the direction of dual-objective optimization: minF =[ f 1, f 2).

[0043] Preferably, the constraints for system operation in step S3 include:

[0044] Wind turbine generator output power constraints:

[0045] P wind-min ≤ P wind-act (t)≤ P wind-rated , t=1,2,...,24;

[0046] in, P wind-act (t) represents the actual output power of the wind turbine generator during time period t. P wind-min This is the minimum operating power of the wind turbine generator set. When the power output is below this value, the wind turbine will enter a shutdown or standby state to avoid mechanical wear caused by frequent start-stop cycles. P wind-rated This is the rated output power of the wind turbine generator set, which is also the maximum safe operating power. Exceeding this value will cause the wind turbine to overload and trigger the protection mechanism.

[0047] SOC constraints:

[0048] SOC boundary constraints:

[0049] SOC min × C bat ≤ SOC (t)× C bat ≤ SOC max × Cbat ; t=1,2,...,24;

[0050] in, SOC min The lower limit of the safe state of charge of an energy storage device is usually 0.2~0.3, or 20%~30%. When it is below this value, the energy storage enters a protection state and deep discharge is prohibited. SOC max The upper limit of the safe state of charge of an energy storage device is usually 0.8~0.9, that is, 80%~90%. When it exceeds this value, the energy storage enters a protection state and charging is prohibited. C bat The rated capacity of the energy storage device, i.e., the maximum amount of electricity that can be stored, is SOC ( t )for t period The state of charge of the energy storage device is 0 ≤ SOC(t) ≤ 1. SOC (t)× C bat The actual stored electricity of the energy storage device in time period t;

[0051] SOC dynamic constraints:

[0052] SOC(t) = SOC(t-1) + [ P bat-ch ( t -1)× η bat - P bat-dis ( t- 1) / η bat ] / C bat ;

[0053] t=2,3,...,24;

[0054] Where SOC(t) is the state of charge of the energy storage device at time t, and SOC(t-1) is the state of charge of the energy storage device at time t-1; P bat-ch (t-1) represents the charging power of the energy storage device at time t-1, i.e., the power of the externally input energy storage. P bat-dis ( t -1) represents the discharge power of the energy storage device in the (t-1)th time period, that is, the power output by the energy storage device. η bat The charge / discharge efficiency of an energy storage device is typically taken as the average of the charge and discharge efficiencies, ranging from 0.85 to 0.95. C bat This refers to the rated capacity of the energy storage device.

[0055] Charging section: Pbat-ch ( t -1)× η bat Indicates the input power during charging. P bat-ch ( t -1) Not 100% can be converted into stored energy, multiplied by efficiency. η bat The actual effective power stored is obtained;

[0056] Discharge section: P bat-dis ( t- 1) / η bat This indicates that during discharge, in order to output externally... P bat-dis ( t- 1) The power required for energy storage is higher, which means the energy storage system needs to consume more electricity, divided by the efficiency. η bat To obtain the actual internal power consumed;

[0057] Based on the previous SOC and charging / discharging power, calculate the current SOC. The initial SOC (1) is set as the current actual SOC of the system.

[0058] Electrolytic cell operating power constraints:

[0059] ;

[0060] in, P elec (t) represents the actual operating power of the electrolytic cell in time period t; P elec (t-1) represents the actual operating power of the electrolytic cell during the (t-1)th time period; P elec-min This is the minimum starting power of the electrolytic cell, i.e., the minimum operating power. Below this value, the electrolysis efficiency will drop sharply, and the cell will usually be shut down. P elec-rated This is the rated operating power of the electrolytic cell, i.e., the maximum safe operating power. Exceeding this value will lead to equipment overload and shortened lifespan. P elec (t)- P elec (t-1)| represents the absolute change in the electrolytic cell's operating power between adjacent time periods; Δ P elec-max This is the maximum allowable power fluctuation range for the electrolyzer, to prevent sudden power changes from causing fluctuations in electrolyte concentration and damage to electrodes;

[0061] Hydrogen storage tank capacity constraints: Based on hydrogen storage tank capacity parameters, the real-time storage range is limited, while simultaneously considering the dynamic relationship between hydrogen production and storage.

[0062] Reserve boundary constraints: Q H2-min ≤ Q H2-act (t)≤ Q H2-max , t=1,2,...,24;

[0063] in, Q H2-min This is the minimum emergency reserve capacity of a hydrogen storage tank (the minimum amount of hydrogen that must be maintained) to cope with sudden hydrogen demand or system failures, and is typically 5% to 10% of the tank's maximum capacity. Q H2-act (t) represents the actual storage capacity of the hydrogen storage tank in time period t, reflecting the current quality of hydrogen stored in the tank; Q H2-max This is the maximum storage capacity limit for the hydrogen storage tank. Exceeding this value will cause the tank to overpressure and trigger the safety protection mechanism.

[0064] Dynamic constraints on reserves: Q H2-act (t)= Q H2-act (t-1)+ Q H2-prod (t-1)- Q H2-use (t-1), t=2,3,...,24;

[0065] in, Q H2-act (t) represents the actual storage capacity of the hydrogen storage tank during time period t. Q H2-act (t-1) represents the actual storage capacity of the hydrogen storage tank during the (t-1)th time period. Q H2-prod (t-1) represents the hydrogen production of the electrolyzer in the (t-1)th time period, which is determined by the electrolysis power and electrolysis efficiency. Q H2-use (t-1) The hydrogen consumption of the system during the t-1 time period, including internal use and external supply;

[0066] Power grid interaction power constraints: Based on the maximum permissible power interaction parameters of the power grid, the range of power purchased and sold is limited, while also considering the power balance of the associated system.

[0067] - P grid-sell-max ≤ P grid (t)≤ P grid-buy-max , t=1,2,...,24;

[0068] inP grid (t) represents the power exchange between the power grids during time period t. P grid-buy-max This is the maximum power that the power grid is allowed to purchase from the system, preventing the system from purchasing too much power from the grid and causing grid overload. P grid-sell-max This is the maximum power output allowed by the power grid system, preventing excessive power injection into the grid that could cause voltage and frequency fluctuations. P grid-sell-max This represents the lower limit of the power that can be sold to the grid, i.e., the negative maximum power that can be sold. It indicates that the power sold by the system to the grid cannot exceed [a certain limit]. P grid-sell-max .

[0069] Preferably, the specific process of solving the constructed low-carbon optimal objective function in step S4 is as follows:

[0070] S41: Define and encode decision variables: Set the operating parameters that the system needs to optimize as decision variables, use real number encoding, and bind the variable dimension with the time granularity to form a single-entity encoding vector:

[0071] X (t)=[ P wind-ac t(t), P bat-ch (t), P bat-dis (t), P elec (t), P grid [(t)];t=1,2,...,24;

[0072] P wind-act (t): Actual output power of the fan during time period t;

[0073] P bat-ch (t): Energy storage charging power during time period t;

[0074] P bat-dis (t): Energy storage discharge power during time period t;

[0075] P elec (t): Operating power of the electrolytic cell during time period t;

[0076] P grid (t): Power grid interaction power during time period t;

[0077] S42: Set core parameters:

[0078] Population size: N = 100~200;

[0079] Maximum number of iterations: G = 50~100;

[0080] Crossover probability: P c =0.8~0.9;

[0081] Mutation probability: P m =0.01~0.05;

[0082] Crowding threshold: C=0.1;

[0083] S43: Generating the initial population based on decision variable constraints and prediction data:

[0084] For each decision variable, an initial value is randomly generated within its constraint interval;

[0085] S44: Perform multi-objective optimization iterative solution.

[0086] Preferably, the specific process of S44 for multi-objective optimization iterative solution is as follows:

[0087] S441: Randomly select 2 parent individuals from the population. X 1(t), X 2(t), perform a crossover operation on each decision variable to generate offspring. X 1'(t), X 2'(t):

[0088] X 1'(t)=0.5[(1+ β 1) X 1(t)+(1- β 1) X 2(t)];

[0089] X 2'(t)=0.5[(1- β 1) X 1(t)+(1+ β 1) X 2(t)];

[0090] in β 1 represents the cross-distribution index;

[0091] For each offspring individual, the decision variable X'(t) variates according to the formula:

[0092] X''(t)=X'(t)+(X max (t)-X min (t))×δ 1;

[0093] Where X''(t) represents the offspring decision variable value after mutation, X max (t) represents the maximum allowable value of the decision variable in time period t, determined by system constraints such as the rated power of the wind turbine and the maximum charging and discharging power of the energy storage. min (t) represents the minimum allowable value of this decision variable in time period t, determined by system constraints, such as the minimum technical output of the wind turbine, the power corresponding to the lower limit of SOC, etc. δ 1 represents the variation distribution index;

[0094] S442: For the population after merging parents and offspring, calculate the two objective function values ​​for each individual. f 1. f 2:

[0095] ;

[0096] ;

[0097] Stratification by dominance: If individual A's... f 1. f If 2 is less than or equal to individual B, and at least one objective is strictly smaller, then A dominates B; individuals without domination are classified as level 1, individuals dominated by level 1 are classified as level 2, and so on.

[0098] For individuals within each layer, according to f 1. f 2. Sort them separately and calculate the objective function distance between each individual and its neighboring individuals:

[0099] ;

[0100] in, CD ( i (for individuals) i The level of congestion, f m max 、f m min The first m The maximum and minimum values ​​of each target are set, and the crowding degree of the boundary individuals is set to infinity (preferential retention). f m ( i+ 1) is the first i +1 individual in the first m Function values ​​on each objective function; f m ( i- 1) is the first i- One individual in the firstm Function values ​​on each objective function;

[0101] S443: Employing a hierarchical selection + crowding ranking strategy, N individuals are selected from the merged population to form the next generation population.

[0102] Select individuals based on the non-dominant sorting hierarchy, first selecting the first level, then the second level, until the number of selected individuals approaches N;

[0103] If the number of individuals in a certain layer exceeds the remaining quota, they will be selected in descending order of crowding.

[0104] Repeat the crossover, mutation, sorting, and selection operations until the number of iterations reaches G, and output the final Pareto optimal solution set;

[0105] S44: Pareto optimal solution set selection:

[0106] Construct a comprehensive evaluation function, calculate the comprehensive score for each solution, and select the solution with the lowest score as the optimal solution:

[0107] The formula for the comprehensive evaluation function is as follows: S = w x f 1 ′ + (1 - w) x f 2 ′ ;

[0108] in, S This represents the individual's overall evaluation score; a lower score indicates a better overall performance of the solution under the dual objectives. w The low-carbon target weighting coefficient reflects the system's priority in minimizing carbon emissions, and its value ranges from 0 to 1. w The larger the value, the more likely it is to choose a low-carbon solution. w The smaller the value, the more likely the economical solution will be chosen. To obtain the normalized carbon emission objective function value, the original f Mapping 1 to the [0,1] interval eliminates the influence of dimensions. To obtain the normalized operating cost objective function value, the original f 2 is mapped to the [0,1] interval to eliminate the influence of dimensions;

[0109] S45: Optimal Operating Parameter Output: Extract the decision variables of the optimal individual to obtain the optimal parameters for each time period over 24 hours, including:

[0110] Target output power of the wind turbine: P wind-obj (t), t=1,2,...,24;

[0111] Target charge / discharge power for energy storage: P bat-ch-obj (t), P bat-dis-obj (t);

[0112] Target operating power of electrolytic cell: P elec-obj (t);

[0113] Target power exchange in the power grid: P grid-obj (t).

[0114] Preferably, the specific process of generating the corresponding dynamic adjustment control command based on the selected optimal operating parameters in step S4 is as follows:

[0115] Wind turbine generator control commands:

[0116] Pitch angle adjustment command: based on target power P wind-obj (t) and current actual power P wind-act The deviation of (t) is used to calculate the pitch angle adjustment Δ. θ ;

[0117] Energy storage devices: mode switching and power correction based on state of charge / discharge and SOC;

[0118] SOC protection and power correction: Real-time acquisition of the current SOC value of energy storage, SOC(t), and correction of the target power if it exceeds the safe range;

[0119] Electrolyzer: Start-stop control and stepped power regulation based on power threshold;

[0120] Power adjustment command generation: Calculate the deviation between the current actual power and the target power to generate a power adjustment command.

[0121] The beneficial effects of this invention include:

[0122] 1. Improve prediction accuracy and provide a reliable basis for optimal control: In the prediction stage, a multi-objective prediction model integrating LSTM and XGBoost is adopted to overcome the limitations of insufficient prediction accuracy of traditional single models. For wind power prediction, the time-series features of historical wind speed are first extracted to output preliminary prediction results; then, non-time-series features are combined and corrected by the XGBoost model, and the two prediction results are weighted and fused, which effectively reduces the impact of wind speed fluctuations, seasonal changes and other factors on prediction accuracy, and significantly reduces the prediction error of wind power in the next 24 hours.

[0123] 2. Achieving Dual-Objective Synergistic Optimization, Balancing Low-Carbon and Economic Efficiency: An optimization function is constructed with the dual objectives of minimizing total carbon emissions and operating costs throughout the system's lifecycle. This comprehensively covers the system's carbon emissions and cost composition, overcoming the limitations of single-objective optimization in traditional control methods. In carbon emission calculations, multiple carbon emission sources are considered, including grid power purchases and electrolyzer operating power. The carbon emission coefficient per unit power of the electrolyzer can be dynamically adjusted according to the actual energy structure, ensuring accurate calculation of total carbon emissions. Cost accounting covers four core cost categories: energy storage charging and discharging losses, electrolyzer energy consumption, hydrogen storage and transportation, and grid power purchases, achieving comprehensive cost quantification. The dual-objective function is solved using NSGA-II to generate a Pareto optimal solution set, and the optimal operating parameters are selected through a comprehensive evaluation function, enabling the system to achieve an optimal balance between low-carbon and economic objectives.

[0124] 3. Comprehensive consideration of system constraints to ensure operational stability: This method systematically analyzes the operational constraints of each component of the wind-powered energy storage hydrogen energy system and the external power grid, including wind turbine output power constraints, energy storage SOC constraints, electrolyzer operating power constraints, and hydrogen storage tank capacity constraints. These constraints are designed to closely align with the operating characteristics of each component and the safety requirements of the power grid. During the optimization process, the optimal operating parameters generated are ensured not to exceed the system's safe operating boundaries, effectively preventing system failures caused by component overload, overcharging / over-discharging of energy storage, and hydrogen storage tank overpressure, thus guaranteeing long-term stable system operation.

[0125] 4. Dynamic Closed-Loop Feedback Adjustment for Enhanced Control Adaptability: A closed-loop feedback adjustment mechanism is introduced, overcoming the limitations of traditional one-time control methods. After the system executes control commands, real-time data on the actual operating status of each component is collected and compared with the optimal operating parameters to calculate the deviation. If the deviation exceeds the allowable error range, the solution parameters of the multi-objective optimization algorithm are immediately readjusted, generating new dynamic adjustment control commands, which are then executed and fed back, forming a dynamic cycle of "prediction-optimization-control-feedback-re-optimization." This ensures the system always operates in its optimal state, avoiding carbon emission exceedances, increased costs, or operational instability caused by initial prediction deviations or external disturbances, significantly improving the system's adaptability to complex operating conditions and changes in the external environment. Attached Figure Description

[0126] Fig. 1 This is a schematic flowchart of the low-carbon optimal dynamic adjustment and control method for wind energy storage hydrogen energy system of the present invention.

[0127] Fig. 2 This is a schematic diagram of the process for constructing the optimal low-carbon objective function according to the present invention.

[0128] Fig. 3 This is a schematic diagram of the architecture of the multi-objective prediction model of the present invention. Detailed Implementation

[0129] The following is in conjunction with the appendix Figs. 1 to 3 The present invention will be further described in detail below:

[0130] Example 1

[0131] See appendix Fig. 1 As shown, the low-carbon optimal dynamic regulation and control method for wind energy storage hydrogen energy systems includes the following steps:

[0132] S1: Real-time acquisition of system operation data, including wind turbine generator data, energy storage device data, hydrogen energy subsystem data, and real-time electricity price, carbon emission coefficient, and load demand data of the power grid. Wind turbine generator data includes real-time wind speed, output power, and pitch angle. Energy storage device data includes state of charge (SOC), charging and discharging power, and charging and discharging efficiency. Hydrogen energy subsystem data includes electrolysis efficiency of the electrolyzer, hydrogen production, and pressure and storage capacity of the hydrogen storage tank.

[0133] S2: Based on system operation data, a multi-objective prediction model is constructed. The multi-objective prediction model uses historical wind speed data, historical grid load data, and historical carbon emission coefficient data as input features. It uses a fusion algorithm of Long Short-Term Memory Network (LSTM) and Gradient Boosting Tree (XGBoost). The LSTM network is specifically designed to perform deep mining on wind speed data with strong time series characteristics to capture the pattern of wind speed changes over time and to make a preliminary prediction of wind power generation based on real-time wind speed. The preliminary wind power generation prediction result is then input into the XGBoost model along with non-time series features (date type, seasonal factors) to correct the preliminary prediction result and improve accuracy. Finally, through weighted fusion, the wind power generation prediction, grid load demand prediction, and carbon emission coefficient prediction are realized for the next 24 hours, respectively, to obtain the wind power generation prediction result, grid load demand prediction result, and carbon emission coefficient prediction result.

[0134] S3: Based on the wind power generation forecast, grid load demand forecast, and carbon emission coefficient forecast, a low-carbon optimal objective function is constructed. The low-carbon optimal objective function has two objectives: minimizing the total carbon emissions and minimizing the operating costs throughout the system's entire life cycle. The operating costs include the charging and discharging loss costs of the energy storage device, the energy consumption costs of the electrolyzer, the hydrogen storage and transportation costs, and the grid purchase costs. At the same time, constraints are set for the system operation, including the output power constraints of the wind turbine generator set, the state of charge constraints of the energy storage device, the operating power constraints of the electrolyzer, the storage capacity constraints of the hydrogen storage tank, and the power interaction constraints of the grid.

[0135] S4: The Non-Dominated Sorting Genetic Algorithm (NSGA-II) is used to solve the constructed low-carbon optimal objective function. Based on the prediction results and combined with the constraints, Pareto optimal solution set that satisfies both low-carbon and economic objectives is obtained through operations such as population initialization, crossover, mutation, non-dominated sorting, and crowding calculation. From the Pareto optimal solution set, the optimal operating parameters are selected, including the target output power of the wind turbine generator, the target charging and discharging power of the energy storage device, the target operating power of the electrolyzer, and the target interaction power with the power grid. According to the selected optimal operating parameters, corresponding dynamic adjustment and control commands are generated, including the pitch angle adjustment command of the wind turbine generator, the charging and discharging switching command and power adjustment command of the energy storage device, and the start-up and shutdown command, power adjustment command, and power control command of the electrolyzer.

[0136] S5: The generated dynamic adjustment and control commands are transmitted to the actuators of each component of the wind energy storage hydrogen energy system. The actuators include the pitch angle adjustment mechanism of the wind turbine generator, the charge and discharge controller of the energy storage device, the power controller of the electrolyzer, and the power controller of the grid interaction interface. Each actuator performs corresponding operations according to the control commands.

[0137] S6: Real-time acquisition of actual operating status data of each component of the system after execution of operations. The actual operating status data includes the actual output power of the wind turbine generator, the actual state of charge and charging / discharging power of the energy storage device, the actual operating power and hydrogen production of the electrolyzer, the actual storage capacity of the hydrogen storage tank, and the actual interaction power with the power grid. The actual operating status data is transmitted to the comparison and analysis module of the central control unit, compared with the optimal operating parameters, and the deviation value (the difference between the actual operating parameters and the optimal operating parameters) is calculated.

[0138] If the deviation is within the allowable error range (the absolute value of the deviation does not exceed 5% of the optimal operating parameters), the current control command is maintained; if the deviation exceeds the allowable error range, the solution parameters of the multi-objective optimization algorithm are readjusted, a new dynamic adjustment control command is generated again, and steps S5 and S6 are repeated to realize the closed-loop feedback adjustment of the system's operating state, ensuring that the wind energy storage hydrogen energy system always operates in a low-carbon optimal state.

[0139] In this embodiment, see Fig. 3 As shown, the multi-objective prediction model in step S2 includes an LSTM model and an XGBoost model;

[0140] The LSTM model uses a 3-layer network architecture, including an input layer, 2 hidden layers, and an output layer. Each hidden layer has 64 neurons, and the activation function is ReLU. A dropout layer (dropout rate=0.2) is added to prevent overfitting.

[0141] The XGBoost model has 100 decision trees, a maximum tree depth of 6, a learning rate of 0.1, and uses root mean square error (RMSE) as the loss function.

[0142] The XGBoost correction model is trained using "LSTM preliminary prediction results + non-time series features" from historical periods as input and the corresponding actual power as the label.

[0143] Example 2

[0144] Based on Example 1, the specific process of wind power generation prediction in step S2 is as follows:

[0145] S21: Obtain at least one year of historical wind speed data, with a time granularity of 15 minutes per data point. Simultaneously, match the historical output power data of wind turbine generators for the corresponding time period for model training labels. Remove invalid data with wind speeds <3m / s (generator not started) or >25m / s (generator shutdown protection). Use the sliding window method with a window length of 24 hours, i.e., 96 15-minute data points. Extract the mean, maximum, fluctuation range, trend slope, and other time-series features of wind speed within each window to form a historical wind speed time-series feature matrix. Filter non-time-series features, including date type (weekday / weekend / holiday) and seasonal factors (spring / summer / autumn / winter), and convert them into numerical features using one-hot encoding. Associate these features with the wind speed time-series feature matrix for later use.

[0146] S22: Input the historical wind speed time series feature matrix into the LSTM model for training, and input the real-time collected wind turbine generator data into the LSTM model to output the preliminary prediction results of wind power generation in the next 24 hours, with a time granularity of 15 minutes / data point and a total of 96 data points.

[0147] S23: Using the preliminary prediction results of wind power generation as the basic feature, input the non-time series features into the XGBoost model to output the corrected prediction results of wind power generation.

[0148] S24: The preliminary wind power generation forecast and the revised wind power generation forecast are weighted and fused to obtain the final wind power generation forecast value, as shown in the following formula:

[0149] Final wind power generation forecast = (Preliminary wind power generation forecast × 0.6) + (Corrected wind power generation forecast × 0.4);

[0150] The 96 merged data points are arranged in chronological order to form a wind power generation forecast for the next 24 hours, in kW.

[0151] The process of predicting grid load demand, carbon emission coefficient, and wind power generation in step S2 is obtained in the same way.

[0152] See Fig. 2 As shown, the specific process of constructing the optimal low-carbon objective function based on the wind power generation forecast results, grid load demand forecast results, and carbon emission coefficient forecast results in step S3 is as follows:

[0153] S31: Integrate the final forecast results of wind power generation, grid load demand, and carbon emission coefficients by time dimension:

[0154] Time alignment: Based on the time granularity of the grid load demand forecast results (1 hour / data point) and carbon emission coefficient forecast results (1 hour / data point), the wind power generation forecast results (15 minutes / data point) are summarized by hour, and the average of 4 data points within each hour is taken as the wind power generation forecast value for that hour, thus forming 3 sets of data sequences for the next 24 hours with a time granularity of 1 hour.

[0155] Wind power generation forecast P wind (t), t=1,2,...,24, unit: kW;

[0156] Forecast of power grid load demand P load (t), t=1,2,...,24, unit: kW;

[0157] Carbon emission coefficient prediction α (t), t=1,2,...,24, unit: kgCO2 / kWh;

[0158] S32: Obtain the inherent parameters of the wind energy storage hydrogen energy system as the basis for calculating the objective function and setting constraints, including:

[0159] Cost parameters: Unit charge / discharge loss cost of energy storage devices δ (Unit: Yuan / kWh) Unit operating energy consumption cost of electrolyzer ε (Unit: Yuan / kWh) Unit storage and transportation costs of hydrogen ζ (Unit: Yuan / kg), Power grid unit electricity purchase cost γ(t) (varying over time).

[0160] Technical parameters: Rated power of wind turbine generator set P wind-rated (Unit: kW); Minimum technical output P wind-min (Unit: kW) Rated capacity of energy storage device C bat (Unit: kWh), State of Charge (SOC) Safety Limit SOC min (e.g., 20%), safety limit SOCmax (e.g., 80%), rated operating power of the electrolytic cell P elec-rated (Unit: kW) Minimum starting power P elec-min (Unit: kW) Maximum storage capacity of hydrogen storage tank Q H2-max (Unit: kg) Emergency Reserves Q H2-min (Unit: kg) Maximum grid-connected power allowed by the power grid P grid-sell-max (Unit: kW) Maximum purchased power P grid-buy-max (Unit: kW)

[0161] Efficiency parameters: Charge and discharge efficiency of energy storage devices η bat (e.g., 90%) Electrolytic cell efficiency η elec (e.g., 75%, which is the amount of electricity consumed to produce 1 kg of hydrogen, unit: kWh / kg).

[0162] S33: Based on the prediction results and system parameters, define the intermediate variables needed in the objective function and constraints, and establish the calculation relationships:

[0163] Energy storage device charging and discharging power: (Assume) P bat-ch (t) represents the energy storage charging power during time period t. P bat-ch (t)≥0, unit: kW; P bat-dis (t) represents the energy storage discharge power during time period t. P bat-dis (t)≥0, unit: kW, satisfying P bat-ch (t)× P bat-dis (t)=0, meaning that simultaneous charging and discharging are not allowed at the same time.

[0164] Energy storage charging and discharging power loss:

[0165] P bat-loss (t)= P bat-ch (t)×(1- η bat )+ P bat-dis (t) / η bat - P bat-dis(t); that is, the charging loss is the difference between the input power and the actual stored power, and the discharging loss is the difference between the stored power and the output power.

[0166] Hydrogen production from the electrolyzer: Q H2-prod (t)= P elec (t)× η elec ;P elec (t) represents the operating power of the electrolyzer during time period t, in kW. Hydrogen production is directly proportional to the electrolysis power.

[0167] Hydrogen storage and transportation volume: Q H2 (t)= Q H2-prod (t)- Q H2-use (t); Q H2-use (t) represents the amount of hydrogen consumed during time period t. If there is no external supply or internal use of hydrogen, it can be set to 0. The unit is kg. That is, the amount of hydrogen that needs to be stored and transported is the difference between the production and consumption. Q H2 (t) represents the amount of hydrogen stored and transported during time period t. Q H2-prod (t) represents the hydrogen production of the electrolyzer during time period t.

[0168] Grid interaction power:

[0169] P grid (t)= P load (t)+ P elec (t)+ P bat-ch (t)- P wind (t)- P bat-dis (t);

[0170] If electricity needs to be purchased from the grid during time period t, the purchased power is: P grid (t); if it represents the electricity sold to the grid during time period t, the electricity sold is - P grid (t).

[0171] S34: With the dual objectives of minimizing total carbon emissions and minimizing operating costs throughout the system's entire life cycle, a low-carbon optimal dual objective function is constructed based on the preprocessed prediction results and intermediate variables.

[0172] In this embodiment, the specific process of step S34 is as follows:

[0173] S341: Construct the first objective function: Minimize the total carbon emissions of the system. minf 1:

[0174] The total carbon emissions of the system come from the carbon emissions from electricity purchases from the grid and the carbon emissions from the operation of the electrolyzers. The calculation formula is based on the predicted carbon emission coefficient and the actual power consumption.

[0175] ;

[0176] in: P grid-buy (t) represents the active power purchased from the grid in time period t, that is, the power purchased from the grid when the system's own power generation is insufficient. P grid-buy (t)= max [ P grid [(t),0], P grid (t) represents the net interaction power between the system and the external power grid at time t, which is a core variable reflecting the direction and magnitude of energy flow between the system and the power grid. P grid When (t)>0, it indicates that the system's own power generation (wind power, energy storage discharge, etc.) is less than the total demand of the load and electrical equipment (such as electrolytic cells), and it needs to purchase electricity from the grid. P grid-buy (t)= P grid (t).

[0177] when P grid When (t)≤0, it indicates that the system's own power generation exceeds the total demand, and the excess electricity will be sold to the grid (or abandoned). In this case, there is no need to purchase electricity from the grid. P grid-buy (t)=0; α (t) is the carbon emission coefficient corresponding to the power grid purchase in time period t, which reflects the carbon emission per unit of electricity generated by the power grid in that time period, and varies with the power grid energy structure and time period; P elec (t) represents the power generation of the system's self-contained power source (such as a diesel generator or fuel cell) during time period t. α (t) represents the carbon emission coefficient corresponding to the power grid's electricity purchase in time period t, reflecting the carbon emission per unit of electricity generated by the power grid in that time period, which varies with the power grid's energy structure and time period. βThis is the carbon emission coefficient per unit power of the electrolyzer, in kgCO2 / kWh. If the electricity for the electrolyzer comes from wind power or energy storage, it can be adjusted according to the actual energy structure. It is included in the total carbon emissions of the system by default. The summation range is the next 24 hours, that is, it covers the carbon emissions of the entire forecast period.

[0178] S342: Constructing the second objective function: minimizing the system operating cost. minf 2:

[0179] Operating costs encompass four categories: energy storage charging and discharging losses, electrolyzer energy consumption, hydrogen storage and transportation, and electricity purchase from the grid. All costs are quantified based on forecast results and intermediate variables.

[0180] ;

[0181] Energy storage charging and discharging loss costs: P bat-loss (t)× δ Based on real-time calculations of power loss and unit loss cost, P bat-loss (t) is the th t The power loss generated during the charging and discharging process of the time-limited energy storage device, and the energy loss caused by charging and discharging efficiency. δ This is the cost coefficient corresponding to the unit power loss of energy storage, reflecting the economic cost of energy storage losses.

[0182] Energy consumption cost of electrolytic cell operation: The cost is directly linked to the operating power of the electrolytic cell; the higher the power, the higher the cost. P elec (t) represents the operating power of the electrolytic cell in time period t. This is the energy consumption cost coefficient corresponding to the unit operating power of the electrolytic cell;

[0183] Hydrogen storage and transportation costs: Q H2 (t)× ζ It is produced only when hydrogen production exceeds consumption, and is directly proportional to the amount stored and transported. Q H2 (t) represents the net increase in hydrogen that needs to be stored or transported in time period t, which is the production minus the consumption. ζ This is the storage and transportation cost coefficient per unit mass of hydrogen.

[0184] Electricity purchase cost from the grid: P grid-buy (t)×γ(t), calculated based on the predicted power purchase and the corresponding electricity price for the time period. P grid-buy γ(t) represents the active power purchased from the grid during time period t, and γ(t) represents the unit price of electricity purchased from the grid during time period t, which fluctuates with time.

[0185] S343: Integration of Dual Objective Functions: Integrating two single objective functions into a unified low-carbon optimal objective function, represented in vector form, clearly defining the direction of dual-objective optimization: minF =[ f 1, f 2).

[0186] The constraints for system operation in step S3 include:

[0187] Wind turbine output power constraints: Based on wind power generation prediction results P wind (t) and the unit's technical parameters limit the actual output power range:

[0188] P wind-min ≤ P wind-act (t)≤ P wind-rated (t=1,2,...,24);

[0189] in, P wind-act (t) represents the actual output power of the wind turbine generator during time period t. P wind-min This is the minimum operating power of the wind turbine generator set. When the power output is below this value, the wind turbine will enter a shutdown or standby state to avoid mechanical wear caused by frequent start-stop cycles. P wind-rated This is the rated output power of the wind turbine generator set, i.e., the maximum safe operating power. Exceeding this value will cause the wind turbine to overload and trigger the protection mechanism.

[0190] Energy storage device state of charge (SOC) constraints: based on rated energy storage capacity C bat This is related to the safety threshold, which limits the real-time range of SOC, and also correlates the dynamic relationship between charge / discharge power and SOC:

[0191] SOC boundary constraints:

[0192] SOC min × C bat ≤ SOC (t)× C bat ≤ SOC max × C bat ; t=1,2,...,24;

[0193] in, SOC minThe lower limit of the safe state of charge of an energy storage device is usually 0.2~0.3, or 20%~30%. When it is below this value, the energy storage enters a protection state and deep discharge is prohibited. SOC max The upper limit of the safe state of charge of an energy storage device is usually 0.8~0.9, that is, 80%~90%. When it exceeds this value, the energy storage enters a protection state and charging is prohibited. C bat The rated capacity of the energy storage device, i.e., the maximum amount of electricity that can be stored, is SOC ( t )for t The state of charge (SOC) of the time-limited energy storage device is 0 ≤ SOC(t) ≤ 1. SOC (t)× C bat The actual stored electricity of the energy storage device during time period t is denoted as t.

[0194] SOC dynamic constraints:

[0195] SOC(t) = SOC(t-1) + [ P bat-ch ( t -1)× η bat - P bat-dis ( t- 1) / η bat ] / C bat ;

[0196] t=2,3,...,24.

[0197] Where SOC(t) is the state of charge of the energy storage device at time t, and SOC(t-1) is the state of charge of the energy storage device at time t-1; P bat-ch (t-1) represents the charging power of the energy storage device at time t-1, i.e., the power of the externally input energy storage. P bat-dis ( t -1) represents the discharge power of the energy storage device in the (t-1)th time period, that is, the power output by the energy storage device. η bat The charge / discharge efficiency of an energy storage device is typically taken as the average of the charge and discharge efficiencies, ranging from 0.85 to 0.95. C bat This refers to the rated capacity of the energy storage device.

[0198] Charging section: P bat-ch ( t -1)× η bat Indicates the input power during charging. Pbat-ch ( t -1) Not 100% can be converted into stored energy, multiplied by efficiency. η bat The actual effective power stored is obtained;

[0199] Discharge section: P bat-dis ( t- 1) / η bat This indicates that during discharge, in order to output externally... P bat-dis ( t- 1) The power required for energy storage is higher, which means the energy storage system needs to consume more electricity, divided by the efficiency. η bat To obtain the actual internal power consumed.

[0200] Based on the previous SOC and charging / discharging power, the current SOC is calculated to ensure that the change in energy storage capacity conforms to the charging / discharging efficiency law. The initial SOC (1) is set as the current actual SOC of the system (obtained from real-time data).

[0201] Electrolytic cell operating power constraints:

[0202] .

[0203] in, P elec (t) represents the actual operating power of the electrolytic cell in time period t; P elec (t-1) represents the actual operating power of the electrolytic cell during the (t-1)th time period; P elec-min This is the minimum starting power of the electrolytic cell, i.e., the minimum operating power. Below this value, the electrolysis efficiency will drop sharply, and the cell will usually be shut down. P elec-rated This is the rated operating power of the electrolytic cell, i.e., the maximum safe operating power. Exceeding this value will lead to equipment overload and shortened lifespan. P elec (t)- P elec (t-1)| represents the absolute change in the electrolytic cell's operating power between adjacent time periods; Δ P elec-max This is the maximum allowable power fluctuation range for the electrolyzer, preventing sudden power changes from causing fluctuations in electrolyte concentration and electrode damage.

[0204] Hydrogen storage tank capacity constraints: Based on hydrogen storage tank capacity parameters, the real-time storage range is limited, while simultaneously considering the dynamic relationship between hydrogen production and storage.

[0205] Reserve boundary constraints: QH2-min ≤ Q H2-act (t)≤ Q H2-max , t=1,2,...,24;

[0206] in, Q H2-min This is the minimum emergency reserve capacity of a hydrogen storage tank (the minimum amount of hydrogen that must be maintained) to cope with sudden hydrogen demand or system failures, and is typically 5% to 10% of the tank's maximum capacity. Q H2-act (t) represents the actual storage capacity of the hydrogen storage tank in time period t, reflecting the current quality of hydrogen stored in the tank; Q H2-max This is the maximum storage capacity limit for the hydrogen storage tank. Exceeding this value will cause the tank to overpressure and trigger the safety protection mechanism.

[0207] Dynamic constraints on reserves: Q H2-act (t)= Q H2-act (t-1)+ Q H2-prod (t-1)- Q H2-use (t-1), t=2,3,...,24.

[0208] in, Q H2-act (t) represents the actual storage capacity of the hydrogen storage tank during time period t. Q H2-act (t-1) represents the actual storage capacity of the hydrogen storage tank during the (t-1)th time period. Q H2-prod (t-1) represents the hydrogen production of the electrolyzer in the (t-1)th time period, which is determined by the electrolysis power and electrolysis efficiency. Q H2-use (t-1) The hydrogen consumption of the system in time period t-1, including internal consumption and external supply. Based on the storage capacity and hydrogen production and sales relationship of the previous time period, calculate the storage capacity at the current time. (Initial time period) Q H2-act (1) Set as the current actual storage of the system, which is obtained from real-time data collection.

[0209] Power grid interaction power constraints: Based on the maximum permissible power interaction parameters of the power grid, the range of power purchased and sold is limited, while also considering the power balance of the associated system.

[0210] - P grid-sell-max ≤ P grid (t)≤ P grid-buy-max, t=1,2,...,24;

[0211] in P grid (t) represents the power exchange between the power grids during time period t. P grid-buy-max This is the maximum power that the power grid is allowed to purchase from the system, preventing the system from purchasing too much power from the grid and causing grid overload. P grid-sell-max This is the maximum power output allowed by the power grid system, preventing excessive power injection into the grid that could cause voltage and frequency fluctuations. P grid-sell-max This represents the lower limit of the power that can be sold to the grid, i.e., the negative maximum power that can be sold. It indicates that the power sold by the system to the grid cannot exceed [a certain limit]. P grid-sell-max This constraint relates to the system power balance equation ( P grid (t)= P load (t)+ P elec (t)+ P bat-ch (t)- P wind (t)- P bat-dis (t) linkage ensures that the interaction power is within a safe range.

[0212] Example 3

[0213] Based on Example 1 or Example 2, the specific process of solving the constructed low-carbon optimal objective function using the Non-dominated Sorting Genetic Algorithm (NSGA-II) in step S4 is as follows:

[0214] S41: Define and encode decision variables: Set the operating parameters that the system needs to optimize as decision variables, use real number encoding (for easy continuous power adjustment), and bind the variable dimension with the time granularity (24 hours) to form a single-entity encoding vector:

[0215] X (t)=[ P wind-ac t(t), P bat-ch (t), P bat-dis (t), P elec (t), P grid [(t)];t=1,2,...,24;

[0216] P wind-act(t): Actual output power of the wind turbine (kW) during time period t, subject to P wind-min ≤ P wind-act (t)≤ P wind-rated Constraints, and with predicted values P wind (t) Deviation ≤±5%;

[0217] P bat-ch (t) and P bat-dis (t): Energy storage charging power and discharging power (kW) during time period t, respectively, satisfying:

[0218] P bat-ch (t)× P bat-dis (t)=0, and 0≤ P bat-ch (t)≤ P bat-ch-max ;

[0219] P elec (t): Electrolytic cell operating power (kW) during time period t, satisfying... P elec (t)=0 or P elec-min ≤ P elec (t)≤ P elec-rated And power fluctuation ≤ Δ P elec-max ;

[0220] P grid (t): Power exchanged between the power grid during time period t (kW), positive for power purchase and negative for power sale, received by - P grid-sell-max ≤ P grid (t)≤ P grid-buy-max constraint.

[0221] S42: Set NSGA-II core parameters:

[0222] Population size: N = 100~200;

[0223] Maximum number of iterations: G = 50~100;

[0224] Crossover probability: P c =0.8~0.9;

[0225] Mutation probability: P m =0.01~0.05;

[0226] Crowding threshold: C=0.1.

[0227] S43: Based on the decision variable constraints and predicted data, an initial population is generated using a "random sampling + constraint verification" method.

[0228] For each decision variable, an initial value is randomly generated within its constraint interval, such as Pwind-act(t) in [ max ( P wind-min 0.95 P wind (t)), min ( P wind-rated 1.05 P wind (t))] internal sampling;

[0229] S44: Perform multi-objective optimization iterative solution, the specific process is as follows:

[0230] S441: Randomly select 2 parent individuals from the population. X 1(t), X 2(t), for each decision variable, generate offspring according to the formula. X 1'(t), X 2'(t):

[0231] X 1'(t)=0.5[(1+ β 1) X 1(t)+(1- β 1) X 2(t)];

[0232] X 2'(t)=0.5[(1- β 1) X 1(t)+(1+ β 1) X 2(t)];

[0233] in β 1 is the crossover index, set to 20, which controls the similarity between offspring and parents. If the offspring exceeds the constraint boundary, it is truncated to the boundary value.

[0234] For each offspring individual, the decision variable X'(t) variates according to the formula:

[0235] X''(t)=X'(t)+(X max (t)-X min(t))× δ 1;

[0236] Where X''(t) represents the offspring decision variable value after mutation, X max (t) represents the maximum allowable value of the decision variable in time period t, determined by system constraints such as the rated power of the wind turbine and the maximum charging and discharging power of the energy storage. min (t) represents the minimum allowable value of this decision variable in time period t, determined by system constraints, such as the minimum technical output of the wind turbine, the power corresponding to the lower limit of SOC, etc. δ 1 is the variation distribution index, with a value range of [-1, 1].

[0237] S442: For the population after merging parents and offspring (size 2N), calculate the two objective function values ​​for each individual. f 1 (Total carbon emissions) f 2 (Operating Costs):

[0238] ;

[0239] ;

[0240] Stratification by dominance: If individual A's... f 1. f If 2 is less than or equal to individual B, and at least one objective is strictly smaller, then A dominates B; individuals without domination are classified as layer 1 (optimal layer), individuals dominated by layer 1 are classified as layer 2, and so on.

[0241] For individuals within each layer, according to f 1. f 2. Sort them separately and calculate the objective function distance between each individual and its neighboring individuals:

[0242] ;

[0243] in, CD ( i (for individuals) i The level of congestion, f m max 、f m min The first m The maximum and minimum values ​​of each target are set, and the crowding degree of the boundary individuals is set to infinity (preferential retention). f m ( i+ 1) is the first i +1 individual in the first m Function values ​​on each objective function; f m (i- 1) is the first i- One individual in the first m Function values ​​on the objective function.

[0244] S443: Employing a hierarchical selection + crowding ranking strategy, N individuals are selected from the merged population to form the next generation population.

[0245] Select individuals based on the non-dominant sorting hierarchy, first selecting the first level, then the second level, until the number of selected individuals approaches N;

[0246] If the number of individuals in a certain layer exceeds the remaining quota, select individuals in descending order of crowding, retaining those with high crowding to ensure that the solution set is evenly distributed in the target space.

[0247] Repeat the crossover, mutation, sorting, and selection operations until the number of iterations reaches G, and output the final Pareto optimal solution set (the set of individuals in the first layer).

[0248] S44: Pareto optimal solution set selection:

[0249] Construct a comprehensive evaluation function, calculate the comprehensive score for each solution, and select the solution with the lowest score as the optimal solution:

[0250] The formula for the comprehensive evaluation function is as follows: S = w x f 1 ′+ (1 -w ) × f 2 ′ ;

[0251] in, S This represents the individual's overall evaluation score; a lower score indicates a better overall performance of the solution under the dual objectives. w The low-carbon target weighting coefficient reflects the system's priority in minimizing carbon emissions, and its value ranges from 0 to 1. w The larger the value, the more likely it is to choose a low-carbon solution. w The smaller the value, the more likely the economical solution will be chosen. To obtain the normalized carbon emission objective function value, the original f Mapping 1 to the [0,1] interval eliminates the influence of dimensions. To obtain the normalized operating cost objective function value, the original f 2 is mapped to the [0,1] interval to eliminate the influence of dimensions.

[0252] S45: Optimal Operating Parameter Output: Extract the decision variables of the optimal individual to obtain the optimal parameters for each time period over 24 hours, including:

[0253] Target output power of the wind turbine: P wind-obj (t) (t=1,2,...,24);

[0254] Target charging power for energy storage: P bat-ch-obj (t);

[0255] Energy storage target discharge power: P bat-dis-obj (t);

[0256] Target operating power of electrolytic cell: P elec-obj (t);

[0257] Target power exchange in the power grid: P grid-obj (t) (including the direction and power of electricity purchase / sale).

[0258] The specific process of generating the corresponding dynamic adjustment and control command based on the selected optimal operating parameters in step S4 is as follows:

[0259] Wind turbine generator control commands:

[0260] Pitch angle adjustment command: based on target power P wind-obj (t) and current actual power P wind-act The deviation of (t) is used to calculate the pitch angle adjustment Δ. θ :

[0261] like P wind-act (t)< P wind-obj (t): Generates a command to reduce the pitch angle, with an adjustment amount Δ θ =-0.5°~-2° (the larger the deviation, the larger the adjustment amount);

[0262] The command format is: Wind turbine pitch angle adjustment: target pitch angle for time period t. θ = θ current +Δ θ Power Target P= P wind-obj (t), by reducing the pitch angle, the wind capture efficiency is improved and the output power is increased;

[0263] like P wind-act (t)> P wind-obj (t); Generates a command to increase the pitch angle, with an adjustment amount Δ θ =+0.5°~-2°, the command format is the same as above, by increasing the pitch angle to reduce the wind energy captured and reduce the output power;

[0264] like P wind-act(t)= P wind-obj (t): Generates a command to maintain the current pitch angle, avoiding frequent adjustments that could lead to mechanical wear.

[0265] Energy storage devices: Mode switching and power correction based on state of charge / discharge and SOC:

[0266] Charge / discharge mode determination:

[0267] based on P bat-ch-obj (t) and P bat-dis-obj The value of (t) determines the energy storage operating mode:

[0268] like P bat-ch-obj (t)>0 (requires charging): Generate a charging mode switching command, formatted as energy storage switching to charging mode, target charging power for time period t. P = P bat-ch-obj (t), and at the same time close the charging circuit switch and disconnect the discharging circuit;

[0269] like P bat-dis-obj (t)>0 (Discharge Required): Generate a discharge mode switching command, formatted as energy storage switching to discharge mode, target discharge power for time period t. P = P bat-dis-obj (t), and at the same time close the discharge circuit switch and disconnect the charging circuit;

[0270] like P bat-ch-obj (t)=P bat-dis-obj (t)=0 (no need for charging and discharging): Generate a "standby mode command" with the format of energy storage entering standby mode, maintaining the current SOC, disconnecting the charging and discharging circuit, and only retaining power supply to the control circuit.

[0271] SOC protection and power correction:

[0272] The current SOC value (SOC(t)) of the energy storage is collected in real time. If it exceeds the safe range, the target power is corrected.

[0273] When SOC(t)≥SOC max ×0.95 (near full charge): Generate a charging power reduction command, formatted as: Energy Storage Power Adjustment: Correct Charging Power for Time Period t P = P bat-ch-obj (t)×0.8, SOC protection to prevent overcharging;

[0274] When SOC(t) ≤ SOC min×1.05 (close to power depletion): Generate a discharge power reduction command, formatted as: Energy storage power adjustment: Correct discharge power P = P during time period t. bat-dis-obj (t)×0.8, SOC protection to prevent over-discharge;

[0275] When SOC(t) is in [SOC min SOC max Within the specified range: Generate a sustaining power command based on the target power to ensure stable charging and discharging.

[0276] Electrolyzer: Start-up / shutdown control and stepped power regulation based on power threshold:

[0277] Start / Stop Status Determination:

[0278] based on P elec-obj The value of (t) determines the operating status of the electrolytic cell:

[0279] like P elec-obj (t)≥ P elec-min Currently, the system is shut down: a start-up command is generated, formatted as "electrolytic cell start-up, power applied after 5 minutes of preheating," and the execution process is "power-on preheating → electrolyte melting → gradual loading to target power" (preheating is to avoid electrode damage caused by a sudden increase in power at low temperatures).

[0280] like P elec-obj (t)=0 and currently running: Generate a shutdown command in the format of "electrolyzer shutdown", execute the pressure relief procedure, and the execution flow is "reduce power to 0 → shut down power → depressurize hydrogen storage tank → inert gas protection".

[0281] like P elec-obj (t) remains at [ P elec-min , P elec-rated Within the specified range and currently running: Generate instructions to maintain operation without requiring state switching.

[0282] Power regulation command generation:

[0283] Calculate the current actual power P elec-act The deviation Δ(t) from the target power P elec = P elec-obj (t)- P elec-act (t):

[0284] If |Δ P elec |≤ΔP elec-max (Fluctuations within allowable range): Generate direct adjustment command, formatted as electrolytic cell power adjustment: target power P = Pelec - obj(t) for time period t, one-time adjustment.

[0285] If |Δ P elec |>Δ P elec-max (Fluctuation Exceeds Upper Limit): Generate a stepped adjustment command, formatted as: Electrolyzer Power Adjustment: Target Power for Time Period t P = P elec-obj (t), adjusted in 3 steps (each step is 1 minute apart, step size Δt). P =Δ P elec / 3, by adjusting in steps, we can avoid sudden changes in electrolyte concentration from affecting electrolysis efficiency.

Claims

1. A low-carbon optimal dynamic adjustment and control method for wind energy storage hydrogen energy systems, characterized in that, Includes the following steps: S1: Real-time acquisition of system operation data, including wind turbine generator data, energy storage device data, hydrogen energy subsystem data, real-time electricity price of the power grid, carbon emission coefficient, and load demand data; S2: Construct a multi-objective prediction model, using historical wind speed data, historical power grid load data, and historical carbon emission coefficient data as input features, and output the wind power generation prediction results, power grid load demand prediction results, and carbon emission coefficient prediction results for a specified future time period. S3: Based on the wind power generation forecast, grid load demand forecast, and carbon emission coefficient forecast, construct the low-carbon optimal objective function and set the system operation constraints. S4: Based on the prediction results and combined with the constraints, the optimal solution set that satisfies the dual objectives of low carbon and economy is obtained through operations such as population initialization, crossover, mutation, non-dominated sorting and crowding calculation. The optimal operating parameters are then selected from the set to generate the corresponding dynamic adjustment and control commands. S5: Transmit dynamic adjustment and control commands to the actuators of each component of the wind energy storage hydrogen energy system, and each actuator performs corresponding operations according to the control commands; S6: Real-time acquisition of actual operating status data of each component of the system after performing operations; The actual operating status data is transmitted to the comparison and analysis module of the central control unit, compared with the optimal operating parameters, and the deviation value is calculated. If the deviation value is within the allowable error range, the current control command is maintained. If the deviation value exceeds the allowable error range, the solution parameters of the multi-objective optimization algorithm are readjusted, a new dynamic adjustment control command is generated again, and steps S5 and S6 are repeated to realize the closed-loop feedback adjustment of the system operating status.

2. The low-carbon optimal dynamic adjustment and control method for wind energy storage hydrogen energy systems according to claim 1, characterized in that, The multi-objective prediction model in step S2 includes the LSTM model and the XGBoost model; The LSTM model uses a 3-layer network architecture, including an input layer, 2 hidden layers and an output layer. Each hidden layer has 64 neurons, the activation function is ReLU, and a Dropout layer is added. The XGBoost model has 100 decision trees, a maximum tree depth of 6, a learning rate of 0.1, and uses root mean square error as the loss function.

3. The low-carbon optimal dynamic adjustment and control method for wind energy storage hydrogen energy systems according to claim 2, characterized in that, The specific process of wind power generation prediction in step S2 is as follows: S21: Obtain at least one year of historical wind speed data, and match it with the historical output power data of wind turbine generators for the corresponding time period. Use the sliding window method to extract the time series features of wind speed, including mean, maximum, fluctuation range and trend slope within each window, to form a historical wind speed time series feature matrix. Filter out non-time series features including date type and seasonal factors, and use one-hot encoding to convert them into numerical features. S22: Input the historical wind speed time series feature matrix into the LSTM model for training, and input the real-time collected wind turbine generator data into the LSTM model to output the preliminary prediction results of wind power generation in the next 24 hours. S23: Using the preliminary prediction results of wind power generation as the basic feature, input the non-time series features into the XGBoost model and output the corrected prediction results of wind power generation. S24: The preliminary prediction results and the revised prediction results of wind power generation are weighted and fused to obtain the final prediction value of wind power generation.

4. The low-carbon optimal dynamic adjustment and control method for wind energy storage hydrogen energy systems according to claim 3, characterized in that, The specific process of constructing the optimal low-carbon objective function based on the wind power generation forecast, grid load demand forecast, and carbon emission coefficient forecast in step S3 is as follows: S31: Integrate the final forecast results of wind power generation, grid load demand, and carbon emission coefficients by time dimension: S32: Obtain the inherent parameters of the wind energy storage hydrogen energy system as the basis for objective function calculation and constraint setting, including cost parameters, technical parameters, and efficiency parameters; S33: Based on the prediction results and system parameters, define the intermediate variables required in the objective function and constraints, and establish the calculation relationship; S34: With the dual objectives of minimizing total carbon emissions and minimizing operating costs throughout the system's entire life cycle, a low-carbon optimal dual objective function is constructed based on the preprocessed prediction results and intermediate variables.

5. The low-carbon optimal dynamic adjustment and control method for wind energy storage hydrogen energy systems according to claim 4, characterized in that, The specific process of step S34 is as follows: S341: Construct the first objective function: Minimize the total carbon emissions of the system. minf 1: The calculation formula is based on the predicted carbon emission factor and actual power consumption: ； in: P grid-buy (t) represents the active power purchased from the grid in time period t. P grid-buy (t)= max [ P grid [(t),0], P grid (t) represents the net interaction power between the system and the external power grid at time t; P elec (t) represents the power generation of the system's self-contained power source during time period t. α (t) represents the carbon emission coefficient corresponding to the power grid's electricity purchase in time period t, reflecting the carbon emission per unit of electricity generated by the power grid in that time period, which varies with the power grid's energy structure and time period. β The carbon emission coefficient per unit power of the electrolytic cell; S342: Constructing the second objective function: minimizing the system operating cost. minf 2: Operating costs encompass four categories: energy storage charging and discharging losses, electrolyzer energy consumption, hydrogen storage and transportation, and electricity purchase from the grid. All costs are quantified based on forecast results and intermediate variables. ； Energy storage charging and discharging loss costs: P bat-loss (t)× δ Based on real-time calculations of power loss and unit loss cost, P bat-loss (t) is the th t Power loss generated during the charging and discharging process of a time-limited energy storage device. δ This is the cost coefficient corresponding to the unit power loss of energy storage. Energy consumption cost of electrolytic cell operation: , P elec (t) represents the operating power of the electrolytic cell in time period t. This is the energy consumption cost coefficient corresponding to the unit operating power of the electrolytic cell; Hydrogen storage and transportation costs: Q H2 (t)× ζ , Q H2 (t) represents the net increase in hydrogen that needs to be stored or transported during time period t. ζ This is the storage and transportation cost coefficient per unit mass of hydrogen. Electricity purchase cost from the grid: P grid-buy (t)×γ(t), P grid-buy γ(t) represents the active power purchased from the grid in time period t, and γ(t) represents the unit price of electricity purchased from the grid in time period t. S343: Integration of Dual Objective Functions: Integrating two single objective functions into a unified low-carbon optimal objective function, represented in vector form, clearly defining the direction of dual-objective optimization: minF =[ f 1, f 2).

6. The low-carbon optimal dynamic adjustment and control method for wind energy storage hydrogen energy systems according to claim 4, characterized in that, The constraints for system operation in step S3 include: Wind turbine generator output power constraints: P wind-min ≤ P wind-act (t)≤ P wind-rated ,t=1,2,...,24; in, P wind-act (t) represents the actual output power of the wind turbine generator during time period t. P wind-min This is the minimum operating power of the wind turbine generator set. P wind-rated This refers to the rated output power of the wind turbine generator set. SOC constraints: SOC boundary constraints: SOC min × C bat ≤ SOC (t)× C bat ≤ SOC max × C bat ;t=1,2,...,24; in, SOC min This is the lower limit of the safe state of charge for energy storage devices; SOC max This represents the safe upper limit for the state of charge of an energy storage device. C bat SOC (State of Charge) is the rated capacity of the energy storage device. t )for t The state of charge (SOC) of the time-limited energy storage device is 0 ≤ SOC(t) ≤ 1. SOC (t)× C bat The actual stored electricity of the energy storage device in time period t; SOC dynamic constraints: SOC(t)= SOC(t-1)+[ P bat-ch ( t -1)× η bat - P bat-dis ( t- 1) / η bat ] / C bat ; t=2,3,...,24； Where SOC(t) is the state of charge of the energy storage device at time t, and SOC(t-1) is the state of charge of the energy storage device at time t-1; P bat-ch (t-1) represents the charging power of the energy storage device at time t-1; P bat-dis ( t -1) represents the discharge power of the energy storage device during the (t-1)th time period; η bat The charging and discharging efficiency of energy storage devices; C bat This refers to the rated capacity of the energy storage device. Based on the previous SOC and charging / discharging power, calculate the current SOC. The initial SOC (1) is set as the current actual SOC of the system. Electrolytic cell operating power constraints: ； in, P elec (t) represents the actual operating power of the electrolytic cell in time period t; P elec (t-1) represents the actual operating power of the electrolytic cell during the (t-1)th time period; P elec-min This is the minimum starting power of the electrolytic cell; P elec-rated The rated operating power of the electrolytic cell; P elec (t)- P elec (t-1)| represents the absolute change in the electrolytic cell's operating power between adjacent time periods; Δ P elec-max This refers to the maximum allowable power fluctuation range of the electrolytic cell; Hydrogen storage tank capacity constraints: Based on hydrogen storage tank capacity parameters, the real-time storage range is limited, while simultaneously considering the dynamic relationship between hydrogen production and storage. Reserve boundary constraints: Q H2-min ≤ Q H2-act (t)≤ Q H2-max , t=1,2,...,24; in, Q H2-min This is the lower limit for the emergency reserve capacity of hydrogen storage tanks; Q H2-act (t) represents the actual storage capacity of the hydrogen storage tank during time period t; Q H2-max This represents the maximum storage capacity limit of the hydrogen storage tank. Dynamic constraints on reserves: Q H2-act (t)= Q H2-act (t-1)+ Q H2-prod (t-1)- Q H2-use (t-1), t=2,3,...,24; in, Q H2-act (t) represents the actual storage capacity of the hydrogen storage tank during time period t. Q H2-act (t-1) represents the actual storage capacity of the hydrogen storage tank during the (t-1)th time period. Q H2-prod (t-1) represents the hydrogen production of the electrolyzer in the (t-1)th time period. Q H2-use (t-1) The hydrogen consumption of the system during the t-1 time period, including internal use and external supply; Power grid interaction power constraints: Based on the maximum permissible power interaction parameters of the power grid, the range of power purchased and sold is limited, while also considering the power balance of the associated system. - P grid-sell-max ≤ P grid (t)≤ P grid-buy-max ,t=1,2,...,24; in, P grid (t) represents the power exchange between the power grids during time period t. P grid-buy-max The maximum power purchase capacity allowed by the power grid system. P grid-sell-max The maximum electricity sales capacity allowed by the power grid system, - P grid-sell-max This represents the lower limit of the power that can be sold to the grid, i.e., the negative maximum power that can be sold. It indicates that the power sold by the system to the grid cannot exceed [a certain limit]. P grid-sell-max .

7. The low-carbon optimal dynamic adjustment and control method for wind energy storage hydrogen energy systems according to claim 5, characterized in that, The specific process of solving the constructed low-carbon optimal objective function in step S4 is as follows: S41: Define and encode decision variables: Set the operating parameters that the system needs to optimize as decision variables, use real number encoding, and bind the variable dimension with the time granularity to form a single-entity encoding vector: X (t)=[ P wind-ac t(t), P bat-ch (t), P bat-dis (t), P elec (t), P grid (t)];t=1,2,...,24; P wind-act (t): Actual output power of the fan during time period t; P bat-ch (t): Energy storage charging power during time period t; P bat-dis (t): Energy storage discharge power during time period t; P elec (t): Operating power of the electrolytic cell during time period t; P grid (t): Power grid interaction power during time period t; S42: Set core parameters: Population size: N = 100~200; Maximum number of iterations: G = 50~100; Crossover probability: P c =0.8~0.9; Mutation probability: P m =0.01~0.05; Crowding threshold: C=0.1; S43: Generating the initial population based on decision variable constraints and prediction data: For each decision variable, an initial value is randomly generated within its constraint interval; S44: Perform multi-objective optimization iterative solution.

8. The low-carbon optimal dynamic adjustment and control method for wind energy storage hydrogen energy systems according to claim 7, characterized in that, The specific process of multi-objective optimization iterative solution for S44 is as follows: S441: Randomly select 2 parent individuals from the population. X 1(t), X 2(t), perform a crossover operation on each decision variable to generate offspring. X 1'(t), X 2'(t): X 1'(t)=0.5[(1+ β 1) X 1(t)+(1- β 1) X 2(t)]; X 2'(t)=0.5[(1- β 1) X 1(t)+(1+ β 1) X 2(t)]; in β 1 represents the cross-distribution index; For each offspring individual, the decision variable X'(t) variates according to the formula: X''(t)=X'(t)+(X max (t)-X min (t))× δ 1; Where X''(t) represents the offspring decision variable value after mutation, X max (t) represents the maximum allowable value of the decision variable in time period t, X min (t) represents the minimum allowable value of the decision variable in time period t. δ 1 represents the variation distribution index; S442: For the population after merging parents and offspring, calculate the two objective function values ​​for each individual. f 1. f 2: ； ； Stratification by dominance: If individual A's... f 1. f If 2 is less than or equal to individual B, and at least one objective is strictly smaller, then A dominates B; individuals without domination are classified as level 1, individuals dominated by level 1 are classified as level 2, and so on. For individuals within each layer, according to f 1. f 2. Sort them separately and calculate the objective function distance between each individual and its neighboring individuals: ； in, CD ( i (for individuals) i The level of congestion, f m max 、f m min The first m The maximum and minimum values ​​of each objective are set, and the crowding degree of boundary individuals is set to infinity. f m ( i+ 1) is the first i +1 individual in the first m Function values ​​on each objective function; f m ( i- 1) is the first i- One individual in the first m Function values ​​on each objective function; S443: Employing a hierarchical selection + crowding ranking strategy, N individuals are selected from the merged population to form the next generation population. Select individuals based on the non-dominant sorting hierarchy, first selecting the first level, then the second level, until the number of selected individuals approaches N; If the number of individuals in a certain layer exceeds the remaining quota, they will be selected in descending order of crowding. Repeat the crossover, mutation, sorting, and selection operations until the number of iterations reaches G, and output the final Pareto optimal solution set; S44: Pareto optimal solution set selection: Construct a comprehensive evaluation function, calculate the comprehensive score for each solution, and select the solution with the lowest score as the optimal solution: The formula for the comprehensive evaluation function is as follows: S=w×f 1 ′+(1-w)× f 2 ′ ; in, S This is the individual's overall evaluation score. w The weighting coefficients for low-carbon targets; This represents the normalized carbon emission objective function value. This represents the normalized operating cost objective function value. S45: Optimal Operating Parameter Output: Extract the decision variables of the optimal individual to obtain the optimal parameters for each time period over 24 hours, including: Target output power of the wind turbine: P wind-obj (t), t=1,2,...,24; Target charging power for energy storage: P bat-ch-obj (t); Energy storage target discharge power: P bat-dis-obj (t); Target operating power of the electrolyzer: P elec-obj (t); Target power exchange in the power grid: P grid-obj (t).

9. The low-carbon optimal dynamic adjustment and control method for wind energy storage hydrogen energy systems according to claim 1, characterized in that, The specific process of generating the corresponding dynamic adjustment and control command based on the selected optimal operating parameters in step S4 is as follows: Wind turbine generator control commands: Pitch angle adjustment command: based on target power P wind-obj (t) and current actual power P wind-act The deviation of (t) is used to calculate the pitch angle adjustment Δ. θ ; Energy storage devices: mode switching and power correction based on state of charge / discharge and SOC; SOC protection and power correction: Real-time acquisition of the current SOC value of energy storage, SOC(t), and correction of the target power if it exceeds the safe range; Electrolytic cells: Start-up and shutdown control and stepped power regulation based on power thresholds; Power adjustment command generation: Calculate the deviation between the current actual power and the target power to generate a power adjustment command.

Images