Multi-energy storage device scheduling method and device

Abstract

The invention provides a multi-energy storage device scheduling method and device, and belongs to the technical field of comprehensive energy system operation regulation and control. The method comprises the following steps: constructing a comprehensive energy system optimization scheduling model comprising the multi-energy storage equipment, and solving the model to obtain a reference scheduling plan of the multi-energy storage equipment; adding the reference scheduling plan as an equality constraint condition into the optimization scheduling model, solving the updated optimization scheduling model, and obtaining a Lagrange multiplier of the equality constraint of the reference scheduling plan of the multi-energy storage device, so as to identify and obtain a rolling scheduling moment set; and then combining actual operation data, establishing a rolling scheduling model of the integrated energy system and solving to obtain a rolling scheduling plan of the multi-energy storage device. According to the method, the uncertainty of new energy output and load is considered, the storage battery, the heat storage device, the cold storage device, the hydrogen storage tank and other multi-energy storage devices are scheduled, the energy efficiency of the system is improved, the energy supply cost is reduced, and the new energy consumption capacity is improved.

Description

Technical Field

[0001] This invention belongs to the field of integrated energy system operation and control technology, and specifically relates to a method and apparatus for scheduling multi-energy storage devices. Background Technology

[0002] Integrated energy systems are energy systems that couple various heterogeneous energy flows, such as electricity, heat, cooling, and gas (natural gas / hydrogen). Optimizing scheduling by utilizing the synergistic characteristics of multiple energy flows in integrated energy systems is a powerful way to improve system energy efficiency, reduce energy supply costs, and absorb renewable energy. Considering the uncertainty of renewable energy and the volatility of multi-energy loads, the role of multi-energy storage devices is becoming increasingly important, necessitating research into scheduling and control methods for multi-energy storage devices under uncertain conditions. Since accurate forecasting information is difficult to obtain in advance, a scheduling plan for multi-energy storage devices is typically derived based on actual renewable energy output and load consumption. During real-time operation, the scheduling plan serves as a benchmark, and the operating strategy is adjusted based on actual renewable energy output and load consumption information.

[0003] Current research on multi-energy storage device scheduling mostly provides reference values ​​for the energy storage status at each point in time, while real-time scheduling aims to make the actual energy storage status as close as possible to the reference values. For example, invention patents with application numbers 201911275796.0 and 202510370199.5 indicate that when there are uncertainties in new energy sources and loads, insufficient energy storage regulation capacity can affect the overall operational safety and economy. Another example is the invention patent with application number 202311332849.4, which requires that the energy storage capacity during rolling scheduling should not be less than the planned value. This may lead to the curtailment of new energy power due to insufficient reserve capacity during periods of high new energy generation. Summary of the Invention

[0004] The purpose of this invention is to overcome the shortcomings of existing technologies and propose a multi-energy storage device scheduling method and apparatus. This invention considers the uncertainties of new energy output and load, and schedules multi-energy storage devices such as batteries, thermal storage devices, cold storage devices, and hydrogen storage tanks to improve system energy efficiency, reduce energy supply costs, and increase the capacity to absorb new energy sources.

[0005] A first aspect of the present invention provides a multi-functional storage device scheduling method, comprising:

[0006] An integrated energy system optimization scheduling model is constructed, which includes multi-energy storage devices, namely: batteries, thermal storage devices, cold storage devices, and hydrogen storage tanks.

[0007] Solving the optimized scheduling model yields the baseline scheduling plan for the multi-energy storage device;

[0008] The baseline scheduling plan is added as an equality constraint to the optimized scheduling model. The updated optimized scheduling model is solved to obtain the Lagrange multipliers of the equality constraint of the baseline scheduling plan for multi-energy storage devices, thereby identifying the set of rolling scheduling times.

[0009] Based on the set of rolling scheduling times and combined with actual operating data, the rolling scheduling plan of the multi-energy storage device is obtained by establishing and solving a rolling scheduling model of the integrated energy system.

[0010] In a specific embodiment of the present invention, the objective function of the integrated energy system optimization scheduling model is to minimize the total operating cost of the integrated energy system, as expressed below:

[0011]

[0012] S D =S G ∪S W ∪S CHP ∪S GB ∪S AC ∪S EB ∪S EC ∪S EH (1)

[0013] In the formula, t is the number of the scheduling period; Υ is the set of all scheduling periods, Υ={1,2,...,N T}, N T This represents the total number of scheduling periods; S represents the active power value of the tie line in the integrated energy system during dispatch period t; Δt represents the time interval between adjacent dispatch periods; i represents the number of any device in the integrated energy system; S D S is the collection of all equipment in an integrated energy system. G S is the set of all thermal power units in an integrated energy system. W S is the set of all new energy generator units in an integrated energy system. CHP S is the collection of all combined heat and power units in an integrated energy system. GB S is the collection of all gas-fired boilers in an integrated energy system. AC S is the collection of all absorption chillers in an integrated energy system. EB S is the collection of all electric boilers in an integrated energy system. EC S is a collection of electric chillers in an integrated energy system; EH It is a collection consisting of electric hydrogen production units in an integrated energy system; The electricity purchase price during the dispatch period t;

[0014] When i∈SG ∪S GB ∪S AC At that time, c i This is the operating cost coefficient for thermal power units, gas-fired boilers, or absorption turbine units; for thermal power units, For thermal power unit i, the active power generated during dispatch period t; for gas-fired boilers, Let be the thermal power generated by gas-fired boiler i during the dispatch period t; for absorption chillers, The thermal power generated by absorption chiller i during the dispatch period t;

[0015] When i∈S CHP At that time, c E,i and c H,i These are the power generation cost coefficient and heating cost coefficient of cogeneration unit i, respectively. and These represent the active power and thermal power generated by cogeneration unit i during the dispatch period t, respectively.

[0016] The constraints of the integrated energy system optimal scheduling model include:

[0017] Equipment operating constraints:

[0018]

[0019]

[0020] In the formula, for thermal power units, P i,min and P i,max These represent the lower and upper limits of the active power generated by thermal power unit i, respectively; for gas-fired boilers, P... i,min and P i,max These represent the lower and upper limits of the thermal power generated by gas-fired boiler i, respectively; for absorption chillers, P... i,min and P i,max These represent the lower and upper limits of the thermal power generated by absorption chiller i, respectively; for electric boilers, P... i,min and P i,max These represent the lower and upper limits of the electrical power consumed by electric boiler i, respectively; for electric chiller, P... i,min and P i,max These are the lower and upper limits of the electrical power consumed by the electric chiller i, respectively;

[0021] For thermal power units For thermal power unit i, the active power generated during the dispatch period t-1; for gas-fired boilers, Let i be the thermal power generated by gas-fired boiler i during the dispatch period t-1; for absorption chillers, For an absorption chiller unit i, the thermal power generated during the dispatch period t-1; for an electric boiler, For electric boiler i, the electrical power consumed during the scheduling period t-1; for electric chiller, Let i be the electrical power consumed by the electric chiller i during the scheduling period t-1;

[0022] R i,up and R i,down These are the upward ramp rate and downward ramp rate of device i, respectively;

[0023] Operating constraints of new energy generator sets:

[0024]

[0025] In the formula, P t i,pre This is the predicted value of the active power generated by the new energy generator unit i during the dispatch period t;

[0026] Operating constraints of combined heat and power units:

[0027]

[0028] In the formula, EP i Let P be the set of feasible domain endpoints of cogeneration unit i. i,k and H i,k These are the active power and thermal power values ​​at the k-th endpoint of the feasible region of cogeneration unit i, respectively. Let i be the k-th combination coefficient of cogeneration unit i during the scheduling period t; Let η be the mass flow rate of hydrogen consumed by the cogeneration unit at time t. CHP The operating efficiency of the combined heat and power unit; LHV is the calorific value of hydrogen combustion;

[0029] Battery operating constraints:

[0030]

[0031] In the formula, S ES It is the collection of all batteries in an integrated energy system; and These represent the charging power and discharging power of battery i during the scheduling period t, respectively. P represents the charge level of battery i during the scheduling period t. c,i,max and P dc,i,max These represent the maximum charging power and the maximum discharging power of battery i, respectively; E i,min and E i,max η represents the minimum and maximum values ​​of the battery capacity i, respectively; c,i and η dc,i These refer to the charging efficiency and discharging efficiency of the battery, respectively. Let i be the charge level of battery i during the scheduling period t-1; These represent the battery's charge level, charging power, and discharging power during the first scheduling period, respectively. The amount of electricity stored in battery i at the initial moment before scheduling;

[0032] Operating constraints of heat storage devices and cold storage devices:

[0033]

[0034] In the formula, S HS S is the collection of all thermal storage devices in an integrated energy system. CS It is the collection of all cold storage devices in an integrated energy system; for heat storage devices, and These represent the heat storage rate and heat release power of heat storage device i during the scheduling period t, respectively; for cold storage device, and These represent the cold storage power and cold release power of cold storage device i during the scheduling period t, respectively; for heat storage device, For heat storage device i, the amount of heat stored during the scheduling period t; for cold storage device, H represents the amount of cold storage device i that is stored during the scheduling period t; for heat storage devices, H represents the amount of cold storage device i that is stored during the scheduling period t. c,i,max and H dc,i,max These represent the maximum heat storage power and the maximum heat release power of heat storage device i, respectively; for cold storage device, H... c,i,max and H dc,i,max These represent the maximum value of the cold storage power and the maximum value of the cold release power of the cold storage device i, respectively; for the heat storage device, E HS,i,min and E HS,i,max These are the minimum and maximum values ​​of the heat storage capacity of heat storage device i, respectively; for cold storage device, E HS,i,min and E HS,i,max These are the minimum and maximum values ​​of the cold storage capacity of cold storage device i, respectively; for heat storage devices, For heat storage device i, the heat stored at the initial moment before scheduling; for cold storage device, For cold storage device i, the amount of cold energy stored at the initial moment before scheduling; for heat storage device i, s i For heat storage device i, s is the self-loss rate; for cold storage device i, s is the self-loss rate. i The self-loss rate of the cold storage device i;

[0035] Hydrogen storage tank operating constraints:

[0036]

[0037] In the formula, S GS It is the collection of all hydrogen storage tanks in an integrated energy system; and These represent the hydrogen storage and hydrogen release mass flow rates of hydrogen storage tank i during scheduling period t, respectively. V represents the amount of hydrogen stored in the thermal or cold storage device i during the scheduling period t; c,i,max and V dc,i,max These represent the maximum hydrogen storage mass flow rate and the maximum hydrogen release mass flow rate of hydrogen storage tank i, respectively; E GS,i,min and E GS,i,max These are the minimum and maximum values ​​of the hydrogen storage capacity of hydrogen storage tank i, respectively. The amount of hydrogen stored in hydrogen storage tank i at the initial moment before scheduling;

[0038] Energy balance constraints:

[0039]

[0040]

[0041] In the formula, η EB,i For the efficiency of electric boiler i, COP i η is the coefficient of performance (COP) of the electric chiller i. EH For the efficiency of the electric hydrogen production unit; P t L,i Let be the active power consumed by electrical load i during the dispatch period t. Let i be the heat power consumed by heat load i during the scheduling period t. S represents the cooling power consumed by cooling load i during the scheduling period t; L S HL S CL It is the collection of all electrical loads, heat loads, and cooling loads within a comprehensive energy system.

[0042] In one specific embodiment of the present invention, obtaining the baseline scheduling plan for the multi-energy storage device includes:

[0043] The predicted power of new energy P t i,pre Predicted power P for electrical load, heat load, and cooling load t L,i , Substituting the parameters into the optimized scheduling model, and then using the branch and bound method to solve the scheduling model with equation (1) as the objective function and equations (2) to (29) as constraints, the optimal solution of the model is the baseline scheduling plan of the integrated energy system. This baseline scheduling plan includes the baseline scheduling plans corresponding to each multi-energy storage device, i.e. The optimal values ​​are denoted as follows:

[0044] In one specific embodiment of the present invention, the identification of the rolling scheduling time set includes:

[0045] 1) Add the equality constraint corresponding to the baseline scheduling plan to the scheduling optimization model:

[0046]

[0047] 2) Solve the scheduling optimization model updated in step 1) using the branch and bound method. The model uses equation (1) as the objective function and equations (2) to (32) as constraints. The Lagrange multipliers corresponding to the constraints (30) to (32) are obtained based on the Lagrange function method, denoted as Let i represent the Lagrange multipliers corresponding to the equality constraints at time t for the storage battery, heat storage device, cold storage device, and hydrogen storage tank, respectively.

[0048] 3) Based on the results of step 2), the Lagrange multipliers at each time point are classified using a binary classification method. The time points corresponding to the group with larger absolute values ​​of the Lagrange multipliers form the rolling scheduling time set Υ'.

[0049] In a specific embodiment of the present invention, obtaining the rolling scheduling plan of the multi-functional storage device includes:

[0050] 1) Establish the objective function of the integrated energy system rolling dispatch model, expressed as follows:

[0051]

[0052] Where σ is a positive number, representing the equivalent weight of the baseline scheduling plan; These represent the deviations of the battery, thermal storage device, cold storage device, and hydrogen storage tank from the baseline scheduling plan at time t, respectively; Υ'={t|t∈Υ,t≥t now} represents the set of scheduling times covered by the rolling scheduling, i.e., the current time t. now Up to the last scheduling moment of the overall scheduling time set Υ;

[0053] 2) Constraints for establishing a rolling dispatch model for an integrated energy system;

[0054] Among them, the constraints of the integrated energy system rolling scheduling model include equations (2) to (29), wherein the scheduling time set Y in equations (2) to (29) is updated to the rolling scheduling time set Y';

[0055] The constraints of the integrated energy system rolling dispatch model also include constraints describing the deviation from the baseline dispatch plan, expressed as follows:

[0056]

[0057] 3) The branch and bound method is used to solve the rolling scheduling model with equation (33) as the objective function and equations (2) to (29) and (34) to (45) as constraints. The optimal solution of the rolling scheduling model is the rolling scheduling plan of the integrated energy system. This rolling scheduling plan includes the rolling scheduling plan of multi-energy storage devices, that is... The optimal value.

[0058] A second aspect of the present invention provides a multi-energy storage device scheduling apparatus, comprising:

[0059] An optimized scheduling model construction module is used to construct an optimized scheduling model for an integrated energy system that includes multiple energy storage devices, such as batteries, thermal storage devices, cold storage devices, and hydrogen storage tanks.

[0060] A baseline scheduling plan generation module is used to solve the optimized scheduling model to obtain the baseline scheduling plan for the multi-energy storage device;

[0061] The Lagrange multiplier calculation module is used to add the baseline scheduling plan as an equality constraint to the optimized scheduling model, solve the updated optimized scheduling model, obtain the Lagrange multipliers of the equality constraint of the multi-energy storage device baseline scheduling plan, and thus identify the rolling scheduling time set.

[0062] The rolling scheduling module is used to obtain the rolling scheduling plan of the multi-energy storage device by establishing and solving a rolling scheduling model of the integrated energy system based on the rolling scheduling time set and actual operating data.

[0063] In a specific embodiment of the present invention, the objective function of the integrated energy system optimization scheduling model is to minimize the total operating cost of the integrated energy system, as expressed below:

[0064]

[0065] S D =S G ∪S W ∪S CHP ∪S GB ∪S AC ∪S EB ∪S EC ∪S EH (1)

[0066] In the formula, t is the number of the scheduling period; Υ is the set of all scheduling periods, Υ={1,2,...,N T}, N T This represents the total number of scheduling periods; S represents the active power value of the tie line in the integrated energy system during dispatch period t; Δt represents the time interval between adjacent dispatch periods; i represents the number of any device in the integrated energy system; S D S is the collection of all equipment in an integrated energy system. G S is the set of all thermal power units in an integrated energy system. W S is the set of all new energy generator units in an integrated energy system. CHP S is the collection of all combined heat and power units in an integrated energy system. GB S is the collection of all gas-fired boilers in an integrated energy system. AC S is the collection of all absorption chillers in an integrated energy system. EB S is the collection of all electric boilers in an integrated energy system. EC S is a collection of electric chillers in an integrated energy system; EH It is a collection consisting of electric hydrogen production units in an integrated energy system; The electricity purchase price during the dispatch period t;

[0067] When i∈S G ∪S GB ∪S AC At that time, c i This is the operating cost coefficient for thermal power units, gas-fired boilers, or absorption turbine units; for thermal power units, For thermal power unit i, the active power generated during dispatch period t; for gas-fired boilers, Let be the thermal power generated by gas-fired boiler i during the dispatch period t; for absorption chillers, The thermal power generated by absorption chiller i during the dispatch period t;

[0068] When i∈S CHP At that time, c E,i and c H,i These are the power generation cost coefficient and heating cost coefficient of cogeneration unit i, respectively. and These represent the active power and thermal power generated by cogeneration unit i during the dispatch period t, respectively.

[0069] The constraints of the integrated energy system optimal scheduling model include:

[0070] Equipment operating constraints:

[0071]

[0072] In the formula, for thermal power units, P i,min and P i,max These represent the lower and upper limits of the active power generated by thermal power unit i, respectively; for gas-fired boilers, P...i,min and P i,max These represent the lower and upper limits of the thermal power generated by gas-fired boiler i, respectively; for absorption chillers, P... i,min and P i,max These represent the lower and upper limits of the thermal power generated by absorption chiller i, respectively; for electric boilers, P... i,min and P i,max These represent the lower and upper limits of the electrical power consumed by electric boiler i, respectively; for electric chiller, P... i,min and P i,max These are the lower and upper limits of the electrical power consumed by the electric chiller i, respectively;

[0073] For thermal power units For thermal power unit i, the active power generated during the dispatch period t-1; for gas-fired boilers, Let i be the thermal power generated by gas-fired boiler i during the dispatch period t-1; for absorption chillers, For an absorption chiller unit i, the thermal power generated during the dispatch period t-1; for an electric boiler, For electric boiler i, the electrical power consumed during the scheduling period t-1; for electric chiller, Let i be the electrical power consumed by the electric chiller i during the scheduling period t-1;

[0074] R i,up and R i,down These are the upward ramp rate and downward ramp rate of device i, respectively;

[0075] Operating constraints of new energy generator sets:

[0076]

[0077] In the formula, P t i,pre This is the predicted value of the active power generated by the new energy generator unit i during the dispatch period t;

[0078] Operating constraints of combined heat and power units:

[0079]

[0080] In the formula, EP i Let P be the set of feasible domain endpoints of cogeneration unit i. i,k and H i,k These are the active power and thermal power values ​​at the k-th endpoint of the feasible region of cogeneration unit i, respectively. Let i be the k-th combination coefficient of cogeneration unit i during the scheduling period t; Let η be the mass flow rate of hydrogen consumed by the cogeneration unit at time t. CHP The operating efficiency of the combined heat and power unit; LHV is the calorific value of hydrogen combustion;

[0081] Battery operating constraints:

[0082]

[0083] In the formula, S ES It is the collection of all batteries in an integrated energy system; and These represent the charging power and discharging power of battery i during the scheduling period t, respectively. P represents the charge level of battery i during the scheduling period t. c,i,max and P dc,i,max These represent the maximum charging power and the maximum discharging power of battery i, respectively; E i,min and E i,max η represents the minimum and maximum values ​​of the battery capacity i, respectively; c,i and η dc,i These refer to the charging efficiency and discharging efficiency of the battery, respectively. Let i be the charge level of battery i during the scheduling period t-1; These represent the battery's charge level, charging power, and discharging power during the first scheduling period, respectively. The amount of electricity stored in battery i at the initial moment before scheduling;

[0084] Operating constraints of heat storage devices and cold storage devices:

[0085]

[0086] In the formula, S HS S is the collection of all thermal storage devices in an integrated energy system. CS It is the collection of all cold storage devices in an integrated energy system; for heat storage devices, and These represent the heat storage rate and heat release power of heat storage device i during the scheduling period t, respectively; for cold storage device, and These represent the cold storage power and cold release power of cold storage device i during the scheduling period t, respectively; for heat storage device, For heat storage device i, the amount of heat stored during the scheduling period t; for cold storage device, H represents the amount of cold storage device i that is stored during the scheduling period t; for heat storage devices, H represents the amount of cold storage device i that is stored during the scheduling period t. c,i,max and H dc,i,max These represent the maximum heat storage power and the maximum heat release power of heat storage device i, respectively; for cold storage device, H... c,i,max and H dc,i,max These represent the maximum value of the cold storage power and the maximum value of the cold release power of the cold storage device i, respectively; for the heat storage device, E HS,i,min and E HS,i,maxThese are the minimum and maximum values ​​of the heat storage capacity of heat storage device i, respectively; for cold storage device, E HS,i,min and E HS,i,max These are the minimum and maximum values ​​of the cold storage capacity of cold storage device i, respectively; for heat storage devices, For heat storage device i, the heat stored at the initial moment before scheduling; for cold storage device, For cold storage device i, the amount of cold energy stored at the initial moment before scheduling; for heat storage device i, s i For heat storage device i, s is the self-loss rate; for cold storage device i, s is the self-loss rate. i The self-loss rate of the cold storage device i;

[0087] Hydrogen storage tank operating constraints:

[0088]

[0089] In the formula, S GS It is the collection of all hydrogen storage tanks in an integrated energy system; and These represent the hydrogen storage and hydrogen release mass flow rates of hydrogen storage tank i during scheduling period t, respectively. V represents the amount of hydrogen stored in the thermal or cold storage device i during the scheduling period t; c,i,max and V dc,i,max These represent the maximum hydrogen storage mass flow rate and the maximum hydrogen release mass flow rate of hydrogen storage tank i, respectively; E GS,i,min and E GS,i,max These are the minimum and maximum values ​​of the hydrogen storage capacity of hydrogen storage tank i, respectively. The amount of hydrogen stored in hydrogen storage tank i at the initial moment before scheduling;

[0090] Energy balance constraints:

[0091]

[0092] In the formula, η EB,i For the efficiency of electric boiler i, COP i η is the coefficient of performance (COP) of the electric chiller i. EH For the efficiency of the electric hydrogen production unit; P t L,i Let be the active power consumed by electrical load i during the dispatch period t. Let i be the heat power consumed by heat load i during the scheduling period t. S represents the cooling power consumed by cooling load i during the scheduling period t; L S HL S CL It is the collection of all electrical loads, heat loads, and cooling loads within a comprehensive energy system.

[0093] In one specific embodiment of the present invention, obtaining the baseline scheduling plan for the multi-energy storage device includes:

[0094] The predicted power of new energy P t i,pre Predicted power P for electrical load, heat load, and cooling load t L,i , Substituting the parameters into the optimized scheduling model, and then using the branch and bound method to solve the scheduling model with equation (1) as the objective function and equations (2) to (29) as constraints, the optimal solution of the model is the baseline scheduling plan of the integrated energy system. This baseline scheduling plan includes the baseline scheduling plans corresponding to each multi-energy storage device, i.e. The optimal values ​​are denoted as follows:

[0095] In one specific embodiment of the present invention, the identification of the rolling scheduling time set includes:

[0096] 1) Add the equality constraint corresponding to the baseline scheduling plan to the scheduling optimization model:

[0097]

[0098]

[0099] 2) Solve the scheduling optimization model updated in step 1) using the branch and bound method. The model uses equation (1) as the objective function and equations (2) to (32) as constraints. The Lagrange multipliers corresponding to the constraints (30) to (32) are obtained based on the Lagrange function method, denoted as Let i represent the Lagrange multipliers corresponding to the equality constraints at time t for the storage battery, heat storage device, cold storage device, and hydrogen storage tank, respectively.

[0100] 3) Based on the results of step 2), the Lagrange multipliers at each time point are classified using a binary classification method. The time points corresponding to the group with larger absolute values ​​of the Lagrange multipliers form the rolling scheduling time set Υ'.

[0101] In a specific embodiment of the present invention, obtaining the rolling scheduling plan of the multi-functional storage device includes:

[0102] 1) Establish the objective function of the integrated energy system rolling dispatch model, expressed as follows:

[0103]

[0104] Where σ is a positive number, representing the equivalent weight of the baseline scheduling plan; These represent the deviations of the battery, thermal storage device, cold storage device, and hydrogen storage tank from the baseline scheduling plan at time t, respectively; Υ'={t|t∈Υ,t≥tnow} represents the set of scheduling times covered by the rolling scheduling, i.e., the current time t. now Up to the last scheduling moment of the overall scheduling time set Υ;

[0105] 2) Constraints for establishing a rolling dispatch model for an integrated energy system;

[0106] Among them, the constraints of the integrated energy system rolling scheduling model include equations (2) to (29), wherein the scheduling time set Y in equations (2) to (29) is updated to the rolling scheduling time set Y';

[0107] The constraints of the integrated energy system rolling dispatch model also include constraints describing the deviation from the baseline dispatch plan, expressed as follows:

[0108]

[0109]

[0110] 3) The branch and bound method is used to solve the rolling scheduling model with equation (33) as the objective function and equations (2) to (29) and (34) to (45) as constraints. The optimal solution of the rolling scheduling model is the rolling scheduling plan of the integrated energy system. This rolling scheduling plan includes the rolling scheduling plan of multi-energy storage devices, that is... The optimal value.

[0111] A third aspect of the present invention provides an electronic device comprising:

[0112] At least one processor; and a memory communicatively connected to said at least one processor;

[0113] The memory stores instructions that can be executed by the at least one processor, and the instructions are configured to execute the aforementioned multi-functional storage device scheduling method.

[0114] A fourth aspect of the present invention provides a computer-readable storage medium storing computer instructions for causing the computer to execute the above-described multi-functional storage device scheduling method.

[0115] The features and beneficial effects of this invention are as follows:

[0116] This invention considers the uncertainties in renewable energy output and load, and when determining the baseline energy storage scheduling plan, only provides the energy storage plan for critical scheduling moments. Critical scheduling moments can be determined by adding equation constraints corresponding to the optimal energy plan for energy storage, and a method for selecting critical scheduling moments based on Lagrange multipliers is proposed to achieve optimal rolling scheduling of multi-energy storage devices, thereby increasing the capacity to absorb renewable energy. Attached Figure Description

[0117] Figure 1 This is an overall flowchart of a multi-energy storage device scheduling method according to an embodiment of the present invention. Detailed Implementation

[0118] This invention proposes a method and apparatus for scheduling multi-energy storage devices, which will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0119] A first aspect of the present invention provides a multi-functional storage device scheduling method, comprising:

[0120] An integrated energy system optimization scheduling model is constructed, which includes multi-energy storage devices, namely: batteries, thermal storage devices, cold storage devices, and hydrogen storage tanks.

[0121] Solving the optimized scheduling model yields the baseline scheduling plan for the multi-energy storage device;

[0122] The baseline scheduling plan is added as an equality constraint to the optimized scheduling model. The updated optimized scheduling model is solved to obtain the Lagrange multipliers of the equality constraint of the baseline scheduling plan for multi-energy storage devices, thereby identifying the set of rolling scheduling times.

[0123] Based on the set of rolling scheduling times and combined with actual operating data, the rolling scheduling plan of the multi-energy storage device is obtained by establishing and solving a rolling scheduling model of the integrated energy system.

[0124] In a specific embodiment of the present invention, the overall process of the multi-energy storage device scheduling method is as follows: Figure 1 As shown, the method includes the following steps:

[0125] 1) Construct an integrated energy system optimization scheduling model that includes multi-energy storage devices.

[0126] In this embodiment, the integrated energy system includes multi-energy storage devices, various loads, new energy power generation, and multi-energy coupling devices (such as gas boilers, absorption chillers, electric boilers, electric chillers, and electric hydrogen generators). The multi-energy storage devices include batteries, thermal storage devices, cold storage devices, and hydrogen storage tanks. The optimization scheduling model consists of an objective function and constraints. The objective function minimizes the total operating cost of the integrated energy system containing multi-energy storage devices. The constraints include unit / equipment constraints, multi-energy storage device constraints, and energy balance constraints. The specific steps are as follows:

[0127] 1-1) Construct the objective function of the integrated energy system optimization scheduling model.

[0128] In this embodiment, the objective function minimizes the total operating cost of the integrated energy system while satisfying various constraints.

[0129]

[0130] In the formula, t is the number of the scheduling period; Υ is the set of all scheduling periods, Υ={1,2,...,N T}, N T N represents the total number of scheduling periods. For example, for daytime scheduling, each scheduling period lasts 15 minutes. T It is 96; S represents the active power value of the tie line in the integrated energy system during dispatch period t; Δt represents the time interval between adjacent dispatch periods; i represents the number of any device in the integrated energy system; S D S is the collection of all equipment in an integrated energy system. G S is the set of all thermal power units in an integrated energy system. W S is the set of all new energy generator units in an integrated energy system. CHP S is the collection of all combined heat and power units in an integrated energy system. GB S is the collection of all gas-fired boilers in an integrated energy system. AC S is the collection of all absorption chillers in an integrated energy system. EB S is the collection of all electric boilers in an integrated energy system. EC S is a collection of electric chillers in an integrated energy system; EH It is a collection consisting of electric hydrogen production units in an integrated energy system; The electricity purchase price is for the dispatch period t.

[0131] When i∈S G ∪S GB ∪S AC At that time, c i This refers to the operating cost coefficient for traditional thermal power units, gas-fired boilers, or absorption turbine units; for traditional thermal power units, This refers to the active power generated by a traditional thermal power unit i during the dispatch period t; for a gas-fired boiler, Let be the thermal power generated by gas-fired boiler i during the dispatch period t; for absorption chillers, Let t be the thermal power generated by absorption chiller i during the scheduling period t.

[0132] When i∈S CHP At that time, c E,i and c H,i These are the power generation cost coefficient and heating cost coefficient of cogeneration unit i, respectively. and These represent the active power and thermal power generated by cogeneration unit i during the dispatch period t, respectively.

[0133] 1-2) Construct constraints for the integrated energy system optimization scheduling model, including:

[0134] Equipment operating constraints:

[0135]

[0136] In the formula, for thermal power units, P i,min and P i,max These represent the lower and upper limits of the active power generated by thermal power unit i, respectively; for gas-fired boilers, P... i,min and P i,max These represent the lower and upper limits of the thermal power generated by gas-fired boiler i, respectively; for absorption chillers, P... i,min and P i,max These represent the lower and upper limits of the thermal power generated by absorption chiller i, respectively; for electric boilers, P... i,min and P i,max These represent the lower and upper limits of the electrical power consumed by electric boiler i, respectively; for electric chiller, P... i,min and P i,max These are the lower and upper limits of the electrical power consumed by the electric chiller i, respectively.

[0137] For thermal power units For thermal power unit i, the active power generated during the dispatch period t-1; for gas-fired boilers, Let i be the thermal power generated by gas-fired boiler i during the dispatch period t-1; for absorption chillers, For an absorption chiller unit i, the thermal power generated during the dispatch period t-1; for an electric boiler, For electric boiler i, the electrical power consumed during the scheduling period t-1; for electric chiller, Let be the electrical power consumed by electric chiller i during the scheduling period t-1.

[0138] R i,up and R i,down These represent the upward ramp rate and downward ramp rate of device i, respectively.

[0139] Operating constraints of new energy generator sets:

[0140]

[0141] In the formula, P t i,pre This is the predicted value of the active power generated by the new energy generator unit i during the dispatch period t.

[0142] Operating constraints of combined heat and power units:

[0143]

[0144] In the formula, EP i Let P be the set of feasible domain endpoints of cogeneration unit i. i,k and H i,k These are the active power and thermal power values ​​at the k-th endpoint of the feasible region of cogeneration unit i, respectively. Let i be the k-th combination coefficient of cogeneration unit i during the scheduling period t; Let η be the mass flow rate of hydrogen consumed by the cogeneration unit at time t. CHP LHV represents the operating efficiency of the combined heat and power unit; LHV represents the calorific value of hydrogen combustion.

[0145] Battery operating constraints:

[0146]

[0147] In the formula, S ES It is the collection of all batteries in an integrated energy system; and These represent the charging power and discharging power of battery i during the scheduling period t, respectively. P represents the charge level of battery i during the scheduling period t. c,i,max and P dc,i,max These represent the maximum charging power and the maximum discharging power of battery i, respectively; E i,min and E i,max η represents the minimum and maximum values ​​of the battery capacity i, respectively; c,i and η dc,i These refer to the charging efficiency and discharging efficiency of the battery, respectively. Let i be the charge level of battery i during the scheduling period t-1; These represent the battery's charge level, charging power, and discharging power during the first scheduling period, respectively. This refers to the amount of electricity stored in battery i at the initial moment before scheduling.

[0148] Operating constraints of heat storage devices and cold storage devices:

[0149]

[0150] In the formula, S HS S is the collection of all thermal storage devices in an integrated energy system. CS It is the collection of all cold storage devices in an integrated energy system; for heat storage devices, and These represent the heat storage rate and heat release power of heat storage device i during the scheduling period t, respectively; for cold storage device, and These represent the cold storage power and cold release power of cold storage device i during the scheduling period t, respectively; for heat storage device, For heat storage device i, the amount of heat stored during the scheduling period t; for cold storage device, H represents the amount of cold storage device i that is stored during the scheduling period t; for heat storage devices, H represents the amount of cold storage device i that is stored during the scheduling period t. c,i,max and H dc,i,max These represent the maximum heat storage power and the maximum heat release power of heat storage device i, respectively; for cold storage device, H... c,i,max and H dc,i,max These represent the maximum value of the cold storage power and the maximum value of the cold release power of the cold storage device i, respectively; for the heat storage device, E HS,i,min and E HS,i,max These are the minimum and maximum values ​​of the heat storage capacity of heat storage device i, respectively; for cold storage device, E HS,i,min and E HS,i,max These are the minimum and maximum values ​​of the cold storage capacity of cold storage device i, respectively; for heat storage devices, For heat storage device i, the heat stored at the initial moment before scheduling; for cold storage device, For cold storage device i, the amount of cold energy stored at the initial moment before scheduling; for heat storage device i, s i For heat storage device i, s is the self-loss rate; for cold storage device i, s is the self-loss rate. i Let i be the self-loss rate of the cold storage device i.

[0151] Hydrogen storage tank operating constraints:

[0152]

[0153]

[0154] In the formula, S GS It is the collection of all hydrogen storage tanks in an integrated energy system; and These represent the hydrogen storage and hydrogen release mass flow rates of hydrogen storage tank i during scheduling period t, respectively. V represents the amount of hydrogen stored in thermal (cold) storage device i during the scheduling period t; c,i,max and V dc ,i,max These represent the maximum hydrogen storage mass flow rate and the maximum hydrogen release mass flow rate of hydrogen storage tank i, respectively; E GS,i,min and E GS,i,max These are the minimum and maximum values ​​of the hydrogen storage capacity of hydrogen storage tank i, respectively. The amount of hydrogen stored in hydrogen storage tank i at the initial moment before scheduling.

[0155] Energy balance constraints:

[0156]

[0157] In the formula, ηEB,i For the efficiency of electric boiler i, COP i η is the coefficient of performance (COP) of the electric chiller i. EH For the efficiency of the electric hydrogen production unit; P t L,i Let be the active power consumed by electrical load i during the dispatch period t. Let i be the heat power consumed by heat load i during the scheduling period t. S represents the cooling power consumed by cooling load i during the scheduling period t; L S HL S CL It is the collection of all electrical loads, heat loads, and cooling loads within a comprehensive energy system.

[0158] 2) Solve the optimized scheduling model established in step 1) to obtain the baseline scheduling plan for the multi-energy storage device.

[0159] The obtained new energy predicted power P t i,pre Predicted power P for electrical load, heat load, and cooling load t L,i , Substitute these parameters into the corresponding positions in the optimized scheduling model established in step 1) (i.e., equations (5), (26), (27), and (28)). Then, use the branch and bound method to solve the scheduling model with equation (1) as the objective function and equations (2) to (29) as constraints. The optimal solution of the model obtained is the baseline scheduling plan of the integrated energy system. This baseline scheduling plan includes the baseline scheduling plans corresponding to each multi-energy storage device, i.e. The optimal values ​​are denoted as follows:

[0160] 3) Add the baseline scheduling plan obtained in step 2) as an equality constraint to the optimized scheduling model established in step 1), solve the updated optimized scheduling model, and obtain the Lagrange multipliers of the equality constraints of the baseline scheduling plan for the multi-energy storage device, thereby identifying the rolling scheduling time set of the multi-energy storage device; the specific steps are as follows:

[0161] 3-1) Based on the results of step 2), add the equality constraints corresponding to the baseline scheduling plan to the scheduling optimization model established in step 1):

[0162]

[0163] 3-2) The branch and bound method is used to solve the scheduling optimization model updated in step 3-1). The model takes equation (1) as the objective function and equations (2) to (32) as constraints. The Lagrange multipliers corresponding to the constraints (30) to (32) are obtained based on the Lagrange function method, denoted as Let represent the Lagrange multipliers corresponding to the equation constraints at time t for the battery, thermal storage device, cold storage device, and hydrogen storage tank, respectively. 3-3) Based on the results of step 3-2), the Lagrange multipliers at each time are classified using typical binary classification methods such as thresholding, logistic regression, and support vector machines. The group with the larger absolute value of the Lagrange multipliers constitutes the critical scheduling time Υ'.

[0164] 4) Using the Lagrange multipliers obtained in step 3), combined with actual operating data, a rolling scheduling model of the integrated energy system is established and solved to achieve rolling correction of the scheduling plan for multi-energy storage devices.

[0165] In this embodiment, based on actual operating data of uncertainties such as predicted power of new energy sources and predicted power of load, the subsequent operating strategy is rolled over and corrected. The specific steps are as follows:

[0166] 4-1) Establish the objective function of the integrated energy system rolling dispatch model, expressed as follows:

[0167]

[0168] Wherein, σ is a positive number, meaning the equivalent weight of the baseline scheduling plan; the larger the value, the more the scheduling target values ​​the overall scheduling effect of the day, but the higher the requirement for prediction accuracy; the smaller the value, the more the scheduling target values ​​the current moment, and the value can be selected according to the actual user needs. These represent the deviations of the battery, thermal storage device, cold storage device, and hydrogen storage tank from the baseline scheduling plan at time t, respectively; Υ'={t|t∈Υ,t≥t now} represents the set of scheduling times covered by the rolling scheduling, i.e., the current time t. now The last scheduling moment of the overall scheduling time set Υ; t now Since the previous period has passed and the operational plan has been determined, it is not within the scope of rolling scheduling.

[0169] 4-2) Constraints for establishing a rolling dispatch model for an integrated energy system.

[0170] In this embodiment, the constraints of the integrated energy system rolling scheduling model include equations (2) to (29), wherein the scheduling time set Y in equations (2) to (29) is updated to the scheduling time set Y' during rolling scheduling. The constraints of the integrated energy system rolling scheduling model also include constraints describing the deviation from the baseline scheduling plan, expressed as follows:

[0171]

[0172] 4-3) The branch and bound method is used to solve the rolling scheduling model with equation (33) as the objective function and equations (2) to (29) and (34) to (45) as constraints. The optimal solution of the rolling scheduling model is the rolling scheduling plan of the integrated energy system. This rolling scheduling plan includes the rolling scheduling plan of multi-energy storage devices, that is... The optimal value is determined to achieve optimal rolling scheduling of multi-energy storage devices, thereby improving system energy efficiency, reducing energy supply costs, and increasing the capacity to absorb new energy sources.

[0173] To implement the above embodiments, a second aspect of the present invention provides a multi-energy storage device scheduling apparatus, comprising:

[0174] An optimized scheduling model construction module is used to construct an optimized scheduling model for an integrated energy system that includes multiple energy storage devices, such as batteries, thermal storage devices, cold storage devices, and hydrogen storage tanks.

[0175] A baseline scheduling plan generation module is used to solve the optimized scheduling model to obtain the baseline scheduling plan for the multi-energy storage device;

[0176] The Lagrange multiplier calculation module is used to add the baseline scheduling plan as an equality constraint to the optimized scheduling model, solve the updated optimized scheduling model, obtain the Lagrange multipliers of the equality constraint of the multi-energy storage device baseline scheduling plan, and thus identify the rolling scheduling time set.

[0177] The rolling scheduling module is used to obtain the rolling scheduling plan of the multi-energy storage device by establishing and solving a rolling scheduling model of the integrated energy system based on the rolling scheduling time set and actual operating data.

[0178] In a specific embodiment of the present invention, the objective function of the integrated energy system optimization scheduling model is to minimize the total operating cost of the integrated energy system, as expressed below:

[0179]

[0180] In the formula, t is the number of the scheduling period; Υ is the set of all scheduling periods, Υ={1,2,...,N T}, N T This represents the total number of scheduling periods; S represents the active power value of the tie line in the integrated energy system during dispatch period t; Δt represents the time interval between adjacent dispatch periods; i represents the number of any device in the integrated energy system; S D S is the collection of all equipment in an integrated energy system. G S is the set of all thermal power units in an integrated energy system. WS is the set of all new energy generator units in an integrated energy system. CHP S is the collection of all combined heat and power units in an integrated energy system. GB S is the collection of all gas-fired boilers in an integrated energy system. AC S is the collection of all absorption chillers in an integrated energy system. EB S is the collection of all electric boilers in an integrated energy system. EC S is a collection of electric chillers in an integrated energy system; EH It is a collection consisting of electric hydrogen production units in an integrated energy system; The electricity purchase price during the dispatch period t;

[0181] When i∈S G ∪S GB ∪S AC At that time, c i This is the operating cost coefficient for thermal power units, gas-fired boilers, or absorption turbine units; for thermal power units, For thermal power unit i, the active power generated during dispatch period t; for gas-fired boilers, Let be the thermal power generated by gas-fired boiler i during the dispatch period t; for absorption chillers, The thermal power generated by absorption chiller i during the dispatch period t;

[0182] When i∈S CHP At that time, c E,i and c H,i These are the power generation cost coefficient and heating cost coefficient of cogeneration unit i, respectively. and These represent the active power and thermal power generated by cogeneration unit i during the dispatch period t, respectively.

[0183] The constraints of the integrated energy system optimal scheduling model include:

[0184] Equipment operating constraints:

[0185]

[0186] In the formula, for thermal power units, P i,min and P i,max These represent the lower and upper limits of the active power generated by thermal power unit i, respectively; for gas-fired boilers, P... i,min and P i,max These represent the lower and upper limits of the thermal power generated by gas-fired boiler i, respectively; for absorption chillers, P... i,min and P i,max These represent the lower and upper limits of the thermal power generated by absorption chiller i, respectively; for electric boilers, P... i,min and P i,maxThese represent the lower and upper limits of the electrical power consumed by electric boiler i, respectively; for electric chiller, P... i,min and P i,max These are the lower and upper limits of the electrical power consumed by the electric chiller i, respectively;

[0187] For thermal power units For thermal power unit i, the active power generated during the dispatch period t-1; for gas-fired boilers, Let i be the thermal power generated by gas-fired boiler i during the dispatch period t-1; for absorption chillers, For an absorption chiller unit i, the thermal power generated during the dispatch period t-1; for an electric boiler, For electric boiler i, the electrical power consumed during the scheduling period t-1; for electric chiller, Let i be the electrical power consumed by the electric chiller i during the scheduling period t-1;

[0188] R i,up and R i,down These are the upward ramp rate and downward ramp rate of device i, respectively;

[0189] Operating constraints of new energy generator sets:

[0190]

[0191] In the formula, P t i,pre This is the predicted value of the active power generated by the new energy generator unit i during the dispatch period t;

[0192] Operating constraints of combined heat and power units:

[0193]

[0194] In the formula, EP i Let P be the set of feasible domain endpoints of cogeneration unit i. i,k and H i,k These are the active power and thermal power values ​​at the k-th endpoint of the feasible region of cogeneration unit i, respectively. Let i be the k-th combination coefficient of cogeneration unit i during the scheduling period t; Let η be the mass flow rate of hydrogen consumed by the cogeneration unit at time t. CHP The operating efficiency of the combined heat and power unit; LHV is the calorific value of hydrogen combustion;

[0195] Battery operating constraints:

[0196]

[0197] In the formula, S ES It is the collection of all batteries in an integrated energy system; and These represent the charging power and discharging power of battery i during the scheduling period t, respectively. P represents the charge level of battery i during the scheduling period t. c,i,max and P dc,i,max These represent the maximum charging power and the maximum discharging power of battery i, respectively; E i,min and E i,max η represents the minimum and maximum values ​​of the battery capacity i, respectively; c,i and η dc,i These refer to the charging efficiency and discharging efficiency of the battery, respectively. Let i be the charge level of battery i during the scheduling period t-1; These represent the battery's charge level, charging power, and discharging power during the first scheduling period, respectively. The amount of electricity stored in battery i at the initial moment before scheduling;

[0198] Operating constraints of heat storage devices and cold storage devices:

[0199]

[0200] In the formula, S HS S is the collection of all thermal storage devices in an integrated energy system. CS It is the collection of all cold storage devices in an integrated energy system; for heat storage devices, and These represent the heat storage rate and heat release power of heat storage device i during the scheduling period t, respectively; for cold storage device, and These represent the cold storage power and cold release power of cold storage device i during the scheduling period t, respectively; for heat storage device, For heat storage device i, the amount of heat stored during the scheduling period t; for cold storage device, H represents the amount of cold storage device i that is stored during the scheduling period t; for heat storage devices, H represents the amount of cold storage device i that is stored during the scheduling period t. c,i,max and H dc,i,max These represent the maximum heat storage power and the maximum heat release power of heat storage device i, respectively; for cold storage device, H... c,i,max and H dc,i,max These represent the maximum value of the cold storage power and the maximum value of the cold release power of the cold storage device i, respectively; for the heat storage device, E HS,i,min and E HS,i,max These are the minimum and maximum values ​​of the heat storage capacity of heat storage device i, respectively; for cold storage device, E HS,i,min and E HS,i,max These are the minimum and maximum values ​​of the cold storage capacity of cold storage device i, respectively; for heat storage devices, For heat storage device i, the heat stored at the initial moment before scheduling; for cold storage device, For cold storage device i, the amount of cold energy stored at the initial moment before scheduling; for heat storage device i, s i For heat storage device i, s is the self-loss rate; for cold storage device i, s is the self-loss rate. i The self-loss rate of the cold storage device i;

[0201] Hydrogen storage tank operating constraints:

[0202]

[0203] In the formula, S GS It is the collection of all hydrogen storage tanks in an integrated energy system; and These represent the hydrogen storage and hydrogen release mass flow rates of hydrogen storage tank i during scheduling period t, respectively. V represents the amount of hydrogen stored in the thermal or cold storage device i during the scheduling period t; c,i,max and V dc,i,max These represent the maximum hydrogen storage mass flow rate and the maximum hydrogen release mass flow rate of hydrogen storage tank i, respectively; E GS,i,min and E GS,i,max These are the minimum and maximum values ​​of the hydrogen storage capacity of hydrogen storage tank i, respectively. The amount of hydrogen stored in hydrogen storage tank i at the initial moment before scheduling;

[0204] Energy balance constraints:

[0205]

[0206] In the formula, η EB,i For the efficiency of electric boiler i, COP i η is the coefficient of performance (COP) of the electric chiller i. EH For the efficiency of the electric hydrogen production unit; P t L,i Let be the active power consumed by electrical load i during the dispatch period t. Let i be the heat power consumed by heat load i during the scheduling period t. S represents the cooling power consumed by cooling load i during the scheduling period t; L S HL S CL It is the collection of all electrical loads, heat loads, and cooling loads within a comprehensive energy system.

[0207] In one specific embodiment of the present invention, obtaining the baseline scheduling plan for the multi-energy storage device includes:

[0208] The predicted power of new energy P t i,pre Predicted power P for electrical load, heat load, and cooling load t L,i , Substituting the parameters into the optimized scheduling model, and then using the branch and bound method to solve the scheduling model with equation (1) as the objective function and equations (2) to (29) as constraints, the optimal solution of the model is the baseline scheduling plan of the integrated energy system. This baseline scheduling plan includes the baseline scheduling plans corresponding to each multi-energy storage device, i.e. The optimal values ​​are denoted as follows:

[0209] In one specific embodiment of the present invention, the identification of the rolling scheduling time set includes:

[0210] 1) Add the equality constraint corresponding to the baseline scheduling plan to the scheduling optimization model:

[0211]

[0212] 2) Solve the scheduling optimization model updated in step 1) using the branch and bound method. The model uses equation (1) as the objective function and equations (2) to (32) as constraints. The Lagrange multipliers corresponding to the constraints (30) to (32) are obtained based on the Lagrange function method, denoted as Let i represent the Lagrange multipliers corresponding to the equality constraints at time t for the storage battery, heat storage device, cold storage device, and hydrogen storage tank, respectively.

[0213] 3) Based on the results of step 2), the Lagrange multipliers at each time point are classified using a binary classification method. The time points corresponding to the group with larger absolute values ​​of the Lagrange multipliers form the rolling scheduling time set Υ'.

[0214] In a specific embodiment of the present invention, obtaining the rolling scheduling plan of the multi-functional storage device includes:

[0215] 1) Establish the objective function of the integrated energy system rolling dispatch model, expressed as follows:

[0216]

[0217] Where σ is a positive number, representing the equivalent weight of the baseline scheduling plan; These represent the deviations of the battery, thermal storage device, cold storage device, and hydrogen storage tank from the baseline scheduling plan at time t, respectively; Υ'={t|t∈Υ,t≥t now} represents the set of scheduling times covered by the rolling scheduling, i.e., the current time t. now Up to the last scheduling moment of the overall scheduling time set Υ;

[0218] 2) Constraints for establishing a rolling dispatch model for an integrated energy system;

[0219] Among them, the constraints of the integrated energy system rolling scheduling model include equations (2) to (29), wherein the scheduling time set Y in equations (2) to (29) is updated to the rolling scheduling time set Y';

[0220] The constraints of the integrated energy system rolling dispatch model also include constraints describing the deviation from the baseline dispatch plan, expressed as follows:

[0221]

[0222] 3) The branch and bound method is used to solve the rolling scheduling model with equation (33) as the objective function and equations (2) to (29) and (34) to (45) as constraints. The optimal solution of the rolling scheduling model is the rolling scheduling plan of the integrated energy system. This rolling scheduling plan includes the rolling scheduling plan of multi-energy storage devices, that is... The optimal value.

[0223] This allows for the scheduling of multiple energy storage devices, such as batteries, thermal storage devices, cold storage devices, and hydrogen storage tanks, taking into account the uncertainties of new energy output and load, thereby improving system energy efficiency, reducing energy supply costs, and increasing the capacity to absorb new energy.

[0224] To implement the above embodiments, a third aspect of the present invention provides an electronic device, comprising:

[0225] At least one processor; and a memory communicatively connected to said at least one processor;

[0226] The memory stores instructions that can be executed by the at least one processor, and the instructions are configured to execute the aforementioned multi-functional storage device scheduling method.

[0227] To implement the above embodiments, a fourth aspect of the present invention provides a computer-readable storage medium storing computer instructions for causing the computer to execute the above-described multi-functional storage device scheduling method.

[0228] It should be noted that the computer-readable medium described in this disclosure can be a computer-readable signal medium or a computer-readable storage medium, or any combination thereof. A computer-readable storage medium can be, for example,—but not limited to—an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of a computer-readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination thereof. In this disclosure, a computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in connection with an instruction execution system, apparatus, or device. In this disclosure, a computer-readable signal medium can include a data signal propagated in baseband or as part of a carrier wave, carrying computer-readable program code. Such propagated data signals can take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. A computer-readable signal medium can be any computer-readable medium other than a computer-readable storage medium, which can send, propagate, or transmit a program for use by or in connection with an instruction execution system, apparatus, or device. The program code contained on the computer-readable medium can be transmitted using any suitable medium, including but not limited to: wires, optical fibers, RF (radio frequency), etc., or any suitable combination thereof.

[0229] The aforementioned computer-readable medium may be included in the aforementioned electronic device; or it may exist independently and not assembled into the electronic device. The aforementioned computer-readable medium carries one or more programs, which, when executed by the electronic device, cause the electronic device to perform a multi-functional storage device scheduling method according to the above embodiments.

[0230] Computer program code for performing the operations of this disclosure can be written in one or more programming languages ​​or a combination thereof, including object-oriented programming languages ​​such as Java, Smalltalk, and C++, and conventional procedural programming languages ​​such as the "C" language or similar programming languages. The program code can be executed entirely on the user's computer, partially on the user's computer, as a standalone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving remote computers, the remote computer can be connected to the user's computer via any type of network—including a local area network (LAN) or a wide area network (WAN)—or can be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0231] In the description of this specification, the references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., refer to specific features, structures, materials, or characteristics described in connection with that embodiment or example, which are included in at least one embodiment or example of this application. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of different embodiments or examples.

[0232] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of that feature. In the description of this application, "multiple" means at least two, such as two, three, etc., unless otherwise explicitly specified.

[0233] Any process or method described in the flowchart or otherwise herein can be understood as representing a module, segment, or portion of code comprising one or more executable instructions for implementing a particular logical function or process, and the scope of the preferred embodiments of this application includes additional implementations in which functions may be performed not in the order shown or discussed, including substantially simultaneously or in reverse order depending on the function involved, as will be understood by those skilled in the art to which embodiments of this application pertain.

[0234] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (such as a computer-based system, a processor-included system, or other system that can fetch and execute instructions from, an instruction execution system, apparatus, or device). For the purposes of this specification, "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transmit programs for use by, or in conjunction with, an instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of computer-readable media include: an electrical connection having one or more wires (electronic device), a portable computer disk drive (magnetic device), random access memory (RAM), read-only memory (ROM), erasable and editable read-only memory (EPROM or flash memory), fiber optic devices, and portable optical disc read-only memory (CDROM). Furthermore, computer-readable media can even be paper or other suitable media on which programs can be printed, because programs can be obtained electronically, for example, by optically scanning the paper or other media, followed by editing, interpreting, or otherwise processing as necessary, and then stored in computer memory.

[0235] It should be understood that various parts of this application can be implemented using hardware, software, firmware, or a combination thereof. In the above embodiments, multiple steps or methods can be implemented using software or firmware stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, it can be implemented using any one or a combination of the following techniques known in the art: discrete logic circuits having logic gates for implementing logical functions on data signals, application-specific integrated circuits (ASICs) having suitable combinational logic gates, programmable gate arrays (PGAs), field-programmable gate arrays (FPGAs), etc.

[0236] Those skilled in the art will understand that all or part of the steps of the methods in the above embodiments can be implemented by a program instructing related hardware. The program can be stored in a computer-readable storage medium, and when executed, the program includes one or a combination of the steps of the method embodiments.

[0237] Furthermore, the functional units in the various embodiments of this application can be integrated into a processing module, or each unit can exist physically separately, or two or more units can be integrated into a module. The integrated module can be implemented in hardware or as a software functional module. If the integrated module is implemented as a software functional module and sold or used as an independent product, it can also be stored in a computer-readable storage medium.

[0238] The storage medium mentioned above can be a read-only memory, a disk, or an optical disk, etc. Although embodiments of this application have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting this application. Those skilled in the art can make changes, modifications, substitutions, and variations to the above embodiments within the scope of this application.

Claims

1. A method for scheduling multi-energy storage devices, characterized in that, include: An integrated energy system optimization scheduling model is constructed, which includes multi-energy storage devices, namely: batteries, thermal storage devices, cold storage devices, and hydrogen storage tanks. Solving the optimized scheduling model yields the baseline scheduling plan for the multi-energy storage device; The baseline scheduling plan is added as an equality constraint to the optimized scheduling model. The updated optimized scheduling model is solved to obtain the Lagrange multipliers of the equality constraint of the baseline scheduling plan for multi-energy storage devices, thereby identifying the set of rolling scheduling times. Based on the set of rolling scheduling times and combined with actual operating data, the rolling scheduling plan of the multi-energy storage device is obtained by establishing and solving a rolling scheduling model of the integrated energy system.

2. The method according to claim 1, characterized in that, The objective function of the integrated energy system optimization scheduling model is to minimize the total operating cost of the integrated energy system, as expressed below: In the formula, t is the number of the scheduling period; Υ is the set of all scheduling periods, Υ={1,2,...,N T }, N T This represents the total number of scheduling periods; Δt represents the active power value of the tie line in the integrated energy system during the dispatch period t; Δt is the time interval between adjacent dispatch periods. i represents the number of any device in the integrated energy system; S D S is the collection of all equipment in an integrated energy system. G S is the set of all thermal power units in an integrated energy system. W S is the set of all new energy generator units in an integrated energy system. CHP S is the collection of all combined heat and power units in an integrated energy system. GB S is the collection of all gas-fired boilers in an integrated energy system. AC S is the collection of all absorption chillers in an integrated energy system. EB S is the collection of all electric boilers in an integrated energy system. EC S is a collection of electric chillers in an integrated energy system; EH It is a collection consisting of electric hydrogen production units in an integrated energy system; The electricity purchase price during the dispatch period t; When i∈S G ∪S GB ∪S AC At that time, c i This is the operating cost coefficient for thermal power units, gas-fired boilers, or absorption turbine units; for thermal power units, For thermal power unit i, the active power generated during dispatch period t; for gas-fired boilers, Let be the thermal power generated by gas-fired boiler i during the dispatch period t; for absorption chillers, The thermal power generated by absorption chiller i during the dispatch period t; When i∈S CHP At that time, c E,i and c H,i These are the power generation cost coefficient and heating cost coefficient of cogeneration unit i, respectively. and These represent the active power and thermal power generated by cogeneration unit i during the dispatch period t, respectively. The constraints of the integrated energy system optimal scheduling model include: Equipment operating constraints: In the formula, for thermal power units, P i,min and P i,max These represent the lower and upper limits of the active power generated by thermal power unit i, respectively; for gas-fired boilers, P... i,min and P i,max These represent the lower and upper limits of the thermal power generated by gas-fired boiler i, respectively; for absorption chillers, P... i,min and P i,max These represent the lower and upper limits of the thermal power generated by absorption chiller i, respectively; for electric boilers, P... i,min and P i,max These represent the lower and upper limits of the electrical power consumed by electric boiler i, respectively; for electric chiller, P... i,min and P i,max These are the lower and upper limits of the electrical power consumed by the electric chiller i, respectively; For thermal power units For thermal power unit i, the active power generated during the dispatch period t-1; for gas-fired boilers, Let i be the thermal power generated by gas-fired boiler i during the dispatch period t-1; for absorption chillers, For an absorption chiller unit i, the thermal power generated during the dispatch period t-1; for an electric boiler, For electric boiler i, the electrical power consumed during the scheduling period t-1; for electric chiller, Let i be the electrical power consumed by the electric chiller i during the scheduling period t-1; R i,up and R i,down These are the upward ramp rate and downward ramp rate of device i, respectively; Operating constraints of new energy generator sets: In the formula, P t i,pre This is the predicted value of the active power generated by the new energy generator unit i during the dispatch period t; Operating constraints of combined heat and power units: In the formula, EP i Let P be the set of feasible domain endpoints of cogeneration unit i. i,k and H i,k These are the active power and thermal power values ​​at the k-th endpoint of the feasible region of cogeneration unit i, respectively. Let i be the k-th combination coefficient of cogeneration unit i during the scheduling period t; Let η be the mass flow rate of hydrogen consumed by the cogeneration unit at time t. CHP The operating efficiency of the combined heat and power unit; LHV is the calorific value of hydrogen combustion; Battery operating constraints: In the formula, S ES It is the collection of all batteries in an integrated energy system; and These represent the charging power and discharging power of battery i during the scheduling period t, respectively. P represents the charge level of battery i during the scheduling period t. c,i,max and P dc,i,max These represent the maximum charging power and the maximum discharging power of battery i, respectively; E i,min and E i,max η represents the minimum and maximum values ​​of the battery capacity i, respectively; c,i and η dc,i These refer to the charging efficiency and discharging efficiency of the battery, respectively. Let i be the charge level of battery i during the scheduling period t-1; These represent the battery's charge level, charging power, and discharging power during the first scheduling period, respectively. The amount of electricity stored in battery i at the initial moment before scheduling; Operating constraints of heat storage devices and cold storage devices: In the formula, S HS S is the collection of all thermal storage devices in an integrated energy system. CS It is the collection of all cold storage devices in an integrated energy system; for heat storage devices, and These represent the heat storage rate and heat release power of heat storage device i during the scheduling period t, respectively; for cold storage device, and These represent the cold storage power and cold release power of cold storage device i during the scheduling period t, respectively; for heat storage device, For heat storage device i, the amount of heat stored during the scheduling period t; for cold storage device, H represents the amount of cold storage device i that is stored during the scheduling period t; for heat storage devices, H represents the amount of cold storage device i that is stored during the scheduling period t. c,i,max and H dc,i,max These represent the maximum heat storage power and the maximum heat release power of heat storage device i, respectively; for cold storage device, H... c,i,max and H dc,i,max These represent the maximum value of the cold storage power and the maximum value of the cold release power of the cold storage device i, respectively; for the heat storage device, E HS,i,min and E HS,i,max These are the minimum and maximum values ​​of the heat storage capacity of heat storage device i, respectively; for cold storage device, E HS,i,min and E HS,i,max These are the minimum and maximum values ​​of the cold storage capacity of cold storage device i, respectively; for heat storage devices, For heat storage device i, the heat stored at the initial moment before scheduling; for cold storage device, For cold storage device i, the amount of cold energy stored at the initial moment before scheduling; for heat storage device i, s i For heat storage device i, s is the self-loss rate; for cold storage device i, s is the self-loss rate. i The self-loss rate of the cold storage device i; Hydrogen storage tank operating constraints: In the formula, S GS It is the collection of all hydrogen storage tanks in an integrated energy system; and These represent the hydrogen storage and hydrogen release mass flow rates of hydrogen storage tank i during scheduling period t, respectively. V represents the amount of hydrogen stored in the thermal or cold storage device i during the scheduling period t; c,i,max and V dc ,i,max These represent the maximum hydrogen storage mass flow rate and the maximum hydrogen release mass flow rate of hydrogen storage tank i, respectively; E GS,i,min and E GS,i,max These are the minimum and maximum values ​​of the hydrogen storage capacity of hydrogen storage tank i, respectively. The amount of hydrogen stored in hydrogen storage tank i at the initial moment before scheduling; Energy balance constraints: In the formula, η EB,i For the efficiency of electric boiler i, COP i η is the coefficient of performance (COP) of the electric chiller i. EH For the efficiency of the electric hydrogen production unit; P t L,i Let be the active power consumed by electrical load i during the dispatch period t. Let i be the heat power consumed by heat load i during the scheduling period t. S represents the cooling power consumed by cooling load i during the scheduling period t; L S HL S CL It is the collection of all electrical loads, heat loads, and cooling loads within a comprehensive energy system.

3. The method according to claim 2, characterized in that, The process of obtaining the baseline scheduling plan for the multi-energy storage device includes: The predicted power of new energy P t i,pre Predicted power P for electrical load, heat load, and cooling load t L,i , Substituting the parameters into the optimized scheduling model, and then using the branch and bound method to solve the scheduling model with equation (1) as the objective function and equations (2) to (29) as constraints, the optimal solution of the model is the baseline scheduling plan of the integrated energy system. This baseline scheduling plan includes the baseline scheduling plans corresponding to each multi-energy storage device, i.e. The optimal values ​​are denoted as follows:

4. The method according to claim 3, characterized in that, The identified set of rolling scheduling times includes: 1) Add the equality constraint corresponding to the baseline scheduling plan to the scheduling optimization model: 2) Solve the scheduling optimization model updated in step 1) using the branch and bound method. The model uses equation (1) as the objective function and equations (2) to (32) as constraints. The Lagrange multipliers corresponding to the constraints (30) to (32) are obtained based on the Lagrange function method, denoted as Let i represent the Lagrange multipliers corresponding to the equality constraints at time t for the storage battery, heat storage device, cold storage device, and hydrogen storage tank, respectively. 3) Based on the results of step 2), the Lagrange multipliers at each time point are classified using a binary classification method. The time points corresponding to the group with larger absolute values ​​of the Lagrange multipliers form the rolling scheduling time set Υ'.

5. The method according to claim 4, characterized in that, The process of obtaining the rolling scheduling plan for the multi-functional storage device includes: 1) Establish the objective function of the integrated energy system rolling dispatch model, expressed as follows: Where σ is a positive number, representing the equivalent weight of the baseline scheduling plan; These represent the deviations of the battery, thermal storage device, cold storage device, and hydrogen storage tank from the baseline scheduling plan at time t, respectively; Υ'={t|t∈Υ,t≥t now } represents the set of scheduling times covered by the rolling scheduling, i.e., the current time t. now Up to the last scheduling moment of the overall scheduling time set Υ; 2) Constraints for establishing a rolling dispatch model for an integrated energy system; Among them, the constraints of the integrated energy system rolling scheduling model include equations (2) to (29), wherein the scheduling time set Y in equations (2) to (29) is updated to the rolling scheduling time set Y'; The constraints of the integrated energy system rolling dispatch model also include constraints describing the deviation from the baseline dispatch plan, expressed as follows: 3) The branch and bound method is used to solve the rolling scheduling model with equation (33) as the objective function and equations (2) to (29) and (34) to (45) as constraints. The optimal solution of the rolling scheduling model is the rolling scheduling plan of the integrated energy system. This rolling scheduling plan includes the rolling scheduling plan of multi-energy storage devices, that is... The optimal value.

6. A multi-energy storage device scheduling apparatus, characterized in that, include: An optimized scheduling model construction module is used to construct an optimized scheduling model for an integrated energy system that includes multiple energy storage devices, such as batteries, thermal storage devices, cold storage devices, and hydrogen storage tanks. A baseline scheduling plan generation module is used to solve the optimized scheduling model to obtain the baseline scheduling plan for the multi-energy storage device; The Lagrange multiplier calculation module is used to add the baseline scheduling plan as an equality constraint to the optimized scheduling model, solve the updated optimized scheduling model, obtain the Lagrange multipliers of the equality constraint of the multi-energy storage device baseline scheduling plan, and thus identify the rolling scheduling time set. The rolling scheduling module is used to obtain the rolling scheduling plan of the multi-energy storage device by establishing and solving a rolling scheduling model of the integrated energy system based on the rolling scheduling time set and actual operating data.

7. The apparatus according to claim 6, characterized in that, The objective function of the integrated energy system optimization scheduling model is to minimize the total operating cost of the integrated energy system, as expressed below: In the formula, t is the number of the scheduling period; Υ is the set of all scheduling periods, Υ={1,2,...,N T }, N T This represents the total number of scheduling periods; Δt represents the active power value of the tie line in the integrated energy system during the dispatch period t; Δt is the time interval between adjacent dispatch periods. i represents the number of any device in the integrated energy system; S D S is the collection of all equipment in an integrated energy system. G S is the set of all thermal power units in an integrated energy system. W S is the set of all new energy generator units in an integrated energy system. CHP S is the collection of all combined heat and power units in an integrated energy system. GB S is the collection of all gas-fired boilers in an integrated energy system. AC S is the collection of all absorption chillers in an integrated energy system. EB S is the collection of all electric boilers in an integrated energy system. EC S is a collection of electric chillers in an integrated energy system; EH It is a collection consisting of electric hydrogen production units in an integrated energy system; The electricity purchase price during the dispatch period t; When i∈S G ∪S GB ∪S AC At that time, c i This is the operating cost coefficient for thermal power units, gas-fired boilers, or absorption turbine units; for thermal power units, For thermal power unit i, the active power generated during dispatch period t; for gas-fired boilers, Let be the thermal power generated by gas-fired boiler i during the dispatch period t; for absorption chillers, The thermal power generated by absorption chiller i during the dispatch period t; When i∈S CHP At that time, c E,i and c H,i These are the power generation cost coefficient and heating cost coefficient of cogeneration unit i, respectively. and These represent the active power and thermal power generated by cogeneration unit i during the dispatch period t, respectively. The constraints of the integrated energy system optimal scheduling model include: Equipment operating constraints: In the formula, for thermal power units, P i,min and P i,max These represent the lower and upper limits of the active power generated by thermal power unit i, respectively; for gas-fired boilers, P... i,min and P i,max These represent the lower and upper limits of the thermal power generated by gas-fired boiler i, respectively; for absorption chillers, P... i,min and P i,max These represent the lower and upper limits of the thermal power generated by absorption chiller i, respectively; for electric boilers, P... i,min and P i,max These represent the lower and upper limits of the electrical power consumed by electric boiler i, respectively; for electric chiller, P... i,min and P i,max These are the lower and upper limits of the electrical power consumed by the electric chiller i, respectively; For thermal power units For thermal power unit i, the active power generated during the dispatch period t-1; for gas-fired boilers, Let i be the thermal power generated by gas-fired boiler i during the dispatch period t-1; for absorption chillers, For an absorption chiller unit i, the thermal power generated during the dispatch period t-1; for an electric boiler, For electric boiler i, the electrical power consumed during the scheduling period t-1; for electric chiller, Let i be the electrical power consumed by the electric chiller i during the scheduling period t-1; R i,up and R i,down These are the upward ramp rate and downward ramp rate of device i, respectively; Operating constraints of new energy generator sets: In the formula, P t i,pre This is the predicted value of the active power generated by the new energy generator unit i during the dispatch period t; Operating constraints of combined heat and power units: In the formula, EP i Let P be the set of feasible domain endpoints of cogeneration unit i. i,k and H i,k These are the active power and thermal power values ​​at the k-th endpoint of the feasible region of cogeneration unit i, respectively. Let i be the k-th combination coefficient of cogeneration unit i during the scheduling period t; Let η be the mass flow rate of hydrogen consumed by the cogeneration unit at time t. CHP The operating efficiency of the combined heat and power unit; LHV is the calorific value of hydrogen combustion; Battery operating constraints: In the formula, S ES It is the collection of all batteries in an integrated energy system; and These represent the charging power and discharging power of battery i during the scheduling period t, respectively. P represents the charge level of battery i during the scheduling period t. c,i,max and P dc,i,max These represent the maximum charging power and the maximum discharging power of battery i, respectively; E i,min and E i,max η represents the minimum and maximum values ​​of the battery capacity i, respectively; c,i and η dc,i These refer to the charging efficiency and discharging efficiency of the battery, respectively. Let i be the charge level of battery i during the scheduling period t-1; These represent the battery's charge level, charging power, and discharging power during the first scheduling period, respectively. The amount of electricity stored in battery i at the initial moment before scheduling; Operating constraints of heat storage devices and cold storage devices: In the formula, S HS S is the collection of all thermal storage devices in an integrated energy system. CS It is the collection of all cold storage devices in an integrated energy system; for heat storage devices, and These represent the heat storage rate and heat release power of heat storage device i during the scheduling period t, respectively; for cold storage device, and These represent the cold storage power and cold release power of cold storage device i during the scheduling period t, respectively; for heat storage device, For heat storage device i, the amount of heat stored during the scheduling period t; for cold storage device, H represents the amount of cold storage device i that is stored during the scheduling period t; for heat storage devices, H represents the amount of cold storage device i that is stored during the scheduling period t. c,i,max and H dc,i,max These represent the maximum heat storage power and the maximum heat release power of heat storage device i, respectively; for cold storage device, H... c,i,max and H dc,i,max These represent the maximum value of the cold storage power and the maximum value of the cold release power of the cold storage device i, respectively; for the heat storage device, E HS,i,min and E HS,i,max These are the minimum and maximum values ​​of the heat storage capacity of heat storage device i, respectively; for cold storage device, E HS,i,min and E HS,i,max These are the minimum and maximum values ​​of the cold storage capacity of cold storage device i, respectively; for heat storage devices, For heat storage device i, the heat stored at the initial moment before scheduling; for cold storage device, For cold storage device i, the amount of cold energy stored at the initial moment before scheduling; for heat storage device i, s i For heat storage device i, s is the self-loss rate; for cold storage device i, s is the self-loss rate. i The self-loss rate of the cold storage device i; Hydrogen storage tank operating constraints: In the formula, S GS It is the collection of all hydrogen storage tanks in an integrated energy system; and These represent the hydrogen storage and hydrogen release mass flow rates of hydrogen storage tank i during scheduling period t, respectively. V represents the amount of hydrogen stored in the thermal or cold storage device i during the scheduling period t; c,i,max and V dc,i,max These represent the maximum hydrogen storage mass flow rate and the maximum hydrogen release mass flow rate of hydrogen storage tank i, respectively; E GS,i,min and E GS,i,max These are the minimum and maximum values ​​of the hydrogen storage capacity of hydrogen storage tank i, respectively. The amount of hydrogen stored in hydrogen storage tank i at the initial moment before scheduling; Energy balance constraints: In the formula, η EB,i For the efficiency of electric boiler i, COP i η is the coefficient of performance (COP) of the electric chiller i. EH For the efficiency of the electric hydrogen production unit; P t L,i Let be the active power consumed by electrical load i during the dispatch period t. Let i be the heat power consumed by heat load i during the scheduling period t. S represents the cooling power consumed by cooling load i during the scheduling period t; L S HL S CL It is the collection of all electrical loads, heat loads, and cooling loads within a comprehensive energy system.

8. The apparatus according to claim 7, characterized in that, The process of obtaining the baseline scheduling plan for the multi-energy storage device includes: The predicted power of new energy P t i,pre Predicted power P for electrical load, heat load, and cooling load t L,i , Substituting the parameters into the optimized scheduling model, and then using the branch and bound method to solve the scheduling model with equation (1) as the objective function and equations (2) to (29) as constraints, the optimal solution of the model is the baseline scheduling plan of the integrated energy system. This baseline scheduling plan includes the baseline scheduling plans corresponding to each multi-energy storage device, i.e. The optimal values ​​are denoted as follows:

9. The apparatus according to claim 8, characterized in that, The identified set of rolling scheduling times includes: 1) Add the equality constraint corresponding to the baseline scheduling plan to the scheduling optimization model: 2) Solve the scheduling optimization model updated in step 1) using the branch and bound method. The model uses equation (1) as the objective function and equations (2) to (32) as constraints. The Lagrange multipliers corresponding to the constraints (30) to (32) are obtained based on the Lagrange function method, denoted as Let i represent the Lagrange multipliers corresponding to the equality constraints at time t for the storage battery, heat storage device, cold storage device, and hydrogen storage tank, respectively. 3) Based on the results of step 2), the Lagrange multipliers at each time point are classified using a binary classification method. The time points corresponding to the group with larger absolute values ​​of the Lagrange multipliers form the rolling scheduling time set Υ'.

10. The apparatus according to claim 9, characterized in that, The process of obtaining the rolling scheduling plan for the multi-functional storage device includes: 1) Establish the objective function of the integrated energy system rolling dispatch model, expressed as follows: Where σ is a positive number, representing the equivalent weight of the baseline scheduling plan; These represent the deviations of the battery, thermal storage device, cold storage device, and hydrogen storage tank from the baseline scheduling plan at time t, respectively; Υ'={t|t∈Υ,t≥t now } represents the set of scheduling times covered by the rolling scheduling, i.e., the current time t. now Up to the last scheduling moment of the overall scheduling time set Υ; 2) Constraints for establishing a rolling dispatch model for an integrated energy system; Among them, the constraints of the integrated energy system rolling scheduling model include equations (2) to (29), wherein the scheduling time set Y in equations (2) to (29) is updated to the rolling scheduling time set Y'; The constraints of the integrated energy system rolling dispatch model also include constraints describing the deviation from the baseline dispatch plan, expressed as follows: 3) The branch and bound method is used to solve the rolling scheduling model with equation (33) as the objective function and equations (2) to (29) and (34) to (45) as constraints. The optimal solution of the rolling scheduling model is the rolling scheduling plan of the integrated energy system. This rolling scheduling plan includes the rolling scheduling plan of multi-energy storage devices, that is... The optimal value.

Images