A low-carbon operation optimization method and device for a multi-energy park

Abstract

The present invention discloses a method and device for optimizing low-carbon operation of a multi-energy park. The method for optimizing low-carbon operation of a multi-energy park comprises: calculating the carbon potential of the grid nodes of the multi-energy park according to the prediction curve of the power grid input to the multi-energy park and the carbon flow; constructing a storage energy charging and discharging state judgment model according to the grid node carbon potential; constructing an objective function with the minimum sum of the power generation cost, carbon emission cost and energy storage operation cost, solving the operation curves of various types of generator sets and the energy storage charging and discharging curves in the multi-energy park in combination with the energy storage charging and discharging state judgment model and the constructed constraints; calculating the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system according to the operation curves of various types of generator sets and the energy storage charging and discharging curves. The present invention can optimize the low-carbon operation of the multi-energy park and effectively reduce the carbon emissions of the multi-energy park.

Description

Technical Field

[0001] The present invention relates to the technical field of low-carbon operation of power systems, and in particular to a method and device for optimizing low-carbon operation of a multi-energy park. Background Art

[0002] As energy and climate change issues become increasingly prominent, achieving low-carbon development and reducing excessive consumption of fossil energy have gradually become the common goals of human society. The low-carbon transformation of the power industry is considered to be the key to low-carbon development and the core measure to achieve the goals of "carbon peak and carbon neutrality".

[0003] In recent years, the carbon emission flow theory has brought new directions to low-carbon electricity. The carbon emission flow is regarded as a virtual network flow analogous to the active power flow in the power grid, which is used to characterize the carbon emissions corresponding to the power flowing through a branch. After the introduction of the carbon emission flow, the real-time carbon potential of any node in the system can be obtained according to any known system flow distribution, and the distributed power source can be guided to adjust the output to achieve the minimum carbon emission of the system. At present, the research on low-carbon multi-energy systems based on carbon emission flows is mainly concentrated on the macro-discussion and planning level, and there is relatively little research on optimized operation. Summary of the invention

[0004] In order to overcome the defects of the prior art, the present invention provides a method and device for optimizing the low-carbon operation of a multi-energy park, which can optimize the low-carbon operation of the multi-energy park and effectively reduce the carbon emissions of the multi-energy park.

[0005] In order to solve the above technical problems, in a first aspect, an embodiment of the present invention provides a low-carbon operation optimization method for a multi-energy park, comprising:

[0006] Calculate the carbon potential of the grid nodes in the multi-energy park based on the forecast curve of the power grid’s input of electricity and carbon flow to the multi-energy park;

[0007] According to the carbon potential of the grid node, a model for judging the energy storage charging and discharging state is constructed;

[0008] The objective function is constructed by minimizing the sum of power generation cost, carbon emission cost and energy storage operation cost, and the operation curves of various power generation units and energy storage charge and discharge curves in the multi-energy park are solved by combining the energy storage charge and discharge state judgment model and the constructed constraints;

[0009] According to the operating curves of the various types of generator sets and the energy storage charging and discharging curves, the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system is calculated.

[0010] Furthermore, the grid node carbon potential is:

[0011]

[0012] Among them, ek N is the equivalent carbon emission value on the power generation side caused by the unit electricity consumption of the kth node in the multi-energy park; + is the set of all branches connected to the kth node that have carbon flow flowing into them; l is the branch number; R l is the carbon flow rate of the lth branch, F is the carbon flow of the lth branch, t m is the unit time; P l is the injected power of the lth branch at the kth node.

[0013] Furthermore, the energy storage charging and discharging state judgment model is:

[0014]

[0015] Where M is of order of magnitude 10 4 A positive integer, y1 and y2 are both 0 / 1 variables, is the discharge carbon potential of the energy storage at time t, t∈[0,T], T is the preset period, e k,t is the charging carbon potential of the kth node in the multi-energy park at time t; when y1=0, y2=1, and When y1=1、y2=0, and When y1=y2=0, the energy storage is in the charging state. When the energy storage is off-grid.

[0016] Furthermore, the objective function is:

[0017]

[0018] Wherein, T is the preset period, t∈[0,T]; are the unit operating costs of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set in the multi-energy park, respectively. are the power generation of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set at time t; c emi For carbon tax, is the carbon potential of the nth gas-fired generator set at time t, C ESS is the unit operating cost of energy storage, P t is the charging or discharging power of the energy storage at time t.

[0019] Furthermore, the constraints include line flow constraints, unit operation constraints and multi-energy load balance constraints.

[0020] Further, the line power flow constraint includes a line power constraint and a node voltage constraint;

[0021] The line power constraint is:

[0022] |P l |≤P lmax ,|Q l |≤Q lmax ;

[0023] Among them, P l is the active power of the lth branch, P lmax is the maximum active power of the lth branch; Q l is the reactive power of the lth branch, Q lmax is the maximum reactive power of the lth branch;

[0024] The node voltage constraints are:

[0025] U k,min ≤U k ≤U k,max ;

[0026] Among them, U k is the voltage of the kth node in the multi-energy park, U k,min is the minimum voltage of the kth node, U k,max is the maximum voltage of the kth node.

[0027] Furthermore, the unit operation constraints include unit power constraints and unit climbing constraints;

[0028] The unit power constraint is:

[0029] P jGkmin ≤P jGk,t ≤P jGkmax ;

[0030] Among them, P jGk,t is the output power of the jth unit in the multi-energy park at time t, P jGkmin is the minimum output power of the jth unit, P jGkmax is the maximum output power of the jth unit;

[0031] The unit climbing constraint is:

[0032]

[0033] Among them, P jGk,t-1 is the output power of the jth unit at time t-1, is the maximum output power that can be adjusted downward by the j-th unit in the period [t-1, t], It is the maximum output power that can be adjusted upward by the j-th unit in the time period [t-1, t].

[0034] Furthermore, the multi-energy load balancing constraint is:

[0035]

[0036] Among them, L t is the load demand of the multi-energy park at time t; P t grid is the power transmitted from the power grid to the multi-energy park at time t, They are respectively the power generation capacities of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set in the multi-energy park at time t.

[0037] Furthermore, the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system is:

[0038]

[0039] Among them, E 0 and E' are the carbon emissions of the multi-energy park before and after the construction of the energy storage system, which are composed of the carbon amount transmitted from the power grid to the multi-energy park and the carbon emissions of the units in the multi-energy park; P G0,t and P G0,t ' are the active power transmitted from the power grid to the multi-energy park at time t before and after the construction of the energy storage system; and are the output of the nth gas unit at time t in the multi-energy park before and after the construction of the energy storage system; e G0,t is the carbon potential of the power grid at time t; is the carbon potential of the nth gas generator set at time t; t∈[0,T], T is the preset period.

[0040] In a second aspect, an embodiment of the present invention provides a multi-energy park low-carbon operation optimization device, comprising:

[0041] A grid node carbon potential calculation module is used to calculate the grid node carbon potential of the multi-energy park based on the prediction curve of the power and carbon flow input by the grid to the multi-energy park;

[0042] A state judgment model construction module is used to construct an energy storage charging and discharging state judgment model according to the carbon potential of the grid node;

[0043] The energy storage system operation construction module is used to construct the objective function with the minimum sum of power generation cost, carbon emission cost and energy storage operation cost, and solve the operation curves of various generator sets and energy storage charge and discharge curves in the multi-energy park in combination with the energy storage charge and discharge state judgment model and the constructed constraints;

[0044] The carbon emission reduction calculation module is used to calculate the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system based on the operating curves of the various types of generator sets and the energy storage charging and discharging curves.

[0045] Compared with the prior art, the embodiments of the present invention have the following beneficial effects:

[0046] The carbon potential of the grid nodes of the multi-energy park is calculated based on the prediction curve of the power grid's input of electricity and carbon flow to the multi-energy park; a storage energy charging and discharging state judgment model is constructed based on the carbon potential of the grid nodes; the objective function is constructed with the minimum sum of the power generation cost, carbon emission cost and energy storage operation cost, and the operation curves of various generator sets and energy storage charging and discharging curves in the multi-energy park are solved by combining the energy storage charging and discharging state judgment model and the constructed constraints; the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system is calculated based on the operation curves of various generator sets and energy storage charging and discharging curves, which can optimize the low-carbon operation of the multi-energy park and effectively reduce the carbon emissions of the multi-energy park. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] Figure 1 A schematic flow chart of a multi-energy park low-carbon operation optimization method in the first embodiment of the present invention;

[0048] Figure 2 This is a schematic structural diagram of a multi-energy park low-carbon operation optimization device in the second embodiment of the present invention. DETAILED DESCRIPTION

[0049] The technical solution of the present invention will be described clearly and completely below in conjunction with the accompanying drawings of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0050] It should be noted that the step numbers in the text are only for the convenience of explaining the specific embodiment and do not limit the order of execution of the steps. The method provided in this embodiment can be executed by a related terminal device, and the following description takes the controller as the execution subject as an example.

[0051] like Figure 1 As shown, the first embodiment provides a multi-energy park low-carbon operation optimization method, including steps S1 to S4:

[0052] S1. Calculate the carbon potential of the grid nodes of the multi-energy park based on the prediction curve of the power grid input to the multi-energy park and the carbon flow;

[0053] S2. Construct a model for judging the energy storage charging and discharging status based on the carbon potential of the grid nodes;

[0054] S3. Construct an objective function by minimizing the sum of power generation cost, carbon emission cost and energy storage operation cost, and solve the operation curves of various power generation units and energy storage charge and discharge curves in the multi-energy park by combining the energy storage charge and discharge state judgment model and the constructed constraints;

[0055] S4. Calculate the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system based on the operating curves of various types of generating units and the energy storage charging and discharging curves.

[0056] It should be noted that the multi-energy park combines distributed power sources such as thermal power generating units, photovoltaic power generating units, gas power generating units, energy storage, loads and other related monitoring and protection devices. It is an effective model for aggregating intermittent low-carbon electricity, multi-energy coordinated control, and demand-side load control.

[0057] As an example, in step S1, a prediction curve of the amount of electricity and carbon flow input from the power grid to the multi-energy park is obtained, and the carbon potential of the power grid nodes of the multi-energy park is calculated based on the prediction curve of the amount of electricity and carbon flow input from the power grid to the multi-energy park.

[0058] In step S2, according to the carbon potential of the grid node, an energy storage charge and discharge state judgment model is constructed to judge the energy storage charge and discharge state based on the energy storage charge and discharge state judgment model.

[0059] In step S3, the objective function is constructed by minimizing the sum of the power generation cost, carbon emission cost and energy storage operation cost of the multi-energy park, and the constraint conditions such as line flow constraints, unit operation constraints and multi-energy load balance constraints are considered. According to the objective function, the energy storage charge and discharge state judgment model and the constraint conditions, the operation curves of various generator sets and the energy storage charge and discharge curves in the multi-energy park are solved.

[0060] In step S4, the energy storage charging and discharging state is determined according to the operating curves of various types of generating units and the energy storage charging and discharging curves, the energy storage is connected to the multi-energy park or the energy storage charging and discharging state is adjusted to build an energy storage system, the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system is calculated, and the optimization effect on the low-carbon operation of the multi-energy park is evaluated.

[0061] This embodiment can introduce carbon emission flows to optimize the low-carbon operation of the multi-energy park and effectively reduce the carbon emissions of the multi-energy park.

[0062] In a preferred embodiment, the grid node carbon potential is:

[0063]

[0064] Among them, e k N is the equivalent carbon emission value on the power generation side caused by the unit electricity consumption of the kth node in the multi-energy park; +is the set of all branches connected to the kth node that have carbon flow flowing into them; l is the branch number; R l is the carbon flow rate of the lth branch, F is the carbon flow of the lth branch, t m is the unit time; P l is the injected power of the lth branch at the kth node.

[0065] As an example, the prediction curve of the power grid's input electricity and carbon flow to the multi-energy park is read, and the carbon potential of each node in the grid is calculated using the following formula.

[0066] The branch carbon flow rate is defined as the carbon flow rate that passes through a single branch with the tidal current in unit time. The calculation formula for the branch carbon flow rate is:

[0067]

[0068] In formula (2), R is the branch carbon flow rate, in kgCO2 / h; F is the branch carbon flow rate, in kgCO2; t m is the unit time, in h.

[0069] The physical meaning of the carbon potential of a grid node is the carbon emission value equivalent to the power generation side caused by the consumption of a unit of electricity at the node. The calculation formula for the carbon potential of a grid node is:

[0070]

[0071] In formula (1), e k is the equivalent carbon emission value on the power generation side caused by the unit electricity consumption of the kth node in the multi-energy park, in kgCO2 / kWh; N + is the set of all branches connected to the kth node that have carbon flow flowing into them; l is the branch number; R l is the carbon flow rate of the lth branch; P l is the injected power of the lth branch at the kth node, in kW.

[0072] In a preferred embodiment, the energy storage charging and discharging state judgment model is:

[0073]

[0074] Where M is of order of magnitude 10 4 A positive integer, y1 and y2 are both 0 / 1 variables, is the discharge carbon potential of the energy storage at time t, t∈[0,T], T is the preset period, e k,t is the charging carbon potential of the kth node in the multi-energy park at time t; when y1=0, y2=1, and When y1=1、y2=0, and When y1=y2=0, the energy storage is in the charging state. When the energy storage is off-grid.

[0075] As an example, according to the carbon potential calculation method of the carbon emission flow node, the charging carbon potential e of the kth node in the multi-energy park at time t can be calculated. k,t , the charging carbon potential of each node at time t is obtained, from which the discharging carbon potential of the energy storage at time t can be calculated. The calculation formula for the discharging carbon potential of the energy storage at time t is:

[0076]

[0077] In formula (4), is the discharge carbon potential of the energy storage at time t; η is the charge and discharge efficiency; P i 、e i are the charging power and charging carbon potential of the node at time i, respectively; Q0 ​​and F0 are the remaining power and carbon emission flow of the energy storage at time T0, respectively; Δt is the length of the time period, T is the preset cycle, and its value is 24, which represents the working time of the energy storage throughout the day.

[0078] Since formula (1) is a nonlinear constraint and cannot be solved using the optimization toolbox, this embodiment uses steps (4)-(11) to convexify the constraints of the energy storage element, without repeatedly updating the energy storage boundary conditions, so that the model is converted into a large optimization model, which can be directly solved by calling the optimization solver. Decisions are made from the perspective of multi-period coupling to avoid possible charge and discharge oscillation problems and improve the efficiency of calculations.

[0079] Take the upper bound of the charging carbon potential at each node at each time, that is, The value of is considered according to the carbon potential curve of the power grid node, and the average value of the minimum and maximum carbon potential of the multi-energy park power grid node is taken, then formula (4) can be changed to:

[0080]

[0081]

[0082] By transposing the term, formula (5) can be transformed into:

[0083]

[0084] By converting the product term into a square term, equation (7) can be transformed into:

[0085]

[0086] Formula (8) satisfies the form of second-order cone optimization:

[0087] When F'0≥Q0, (8) is transformed into:

[0088]

[0089] When F'0<Q0, (8) is transformed into:

[0090]

[0091] When judging the charge and discharge states of the energy storage element, 0 / 1 variables y1 and y2 are introduced:

[0092]

[0093] In formula (11), M is a positive integer of the order of 10 4 , y1 and y2 are both 0 / 1 variables, is the discharge carbon potential of the energy storage at time t, t∈[0, T], T is the preset period, e i is the charge carbon potential of the node at time i; when y1 = 0, y2 = 1, and when, the energy storage is in the discharge state, when y1 = 1, y2 = 0, and when, the energy storage is in the charge state, when y1 = y2 = 0, and when, the energy storage is in the off-grid state.

[0094] In a preferred embodiment, the objective function is:

[0095]

[0096] where, T is the preset period, t∈[0, T]; are respectively the unit operating costs of the vth wind turbine, the mth photovoltaic generator, and the nth gas generator in the multi-energy park, are respectively the power generation powers of the vth wind turbine, the mth photovoltaic generator, and the nth gas generator at time t; c emi is the carbon tax, is the carbon potential of the nth gas generator at time t, C ESS is the unit operating cost of the energy storage, P t is the charging or discharging power of the energy storage at time t.

[0097] As an example, assuming that a multi-energy park combines wind turbines (WT), photovoltaic generators (PV), and gas generators (GPP), where the gas generators produce carbon emissions, the objective function is constructed based on the power generation costs of wind turbines, photovoltaic generators, and gas generators, as well as the carbon emission costs of gas generators:

[0098]

[0099] In formula (12), T is the preset period, t∈[0,T]; are the unit operating costs of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set in the multi-energy park, respectively. are the power generation of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set at time t; c emi For carbon tax, is the carbon potential of the nth gas-fired generator set at time t, C ESS is the unit operating cost of energy storage, P t is the charging or discharging power of the energy storage at time t.

[0100] In a preferred embodiment, the constraints include line flow constraints, unit operation constraints and multi-energy load balance constraints.

[0101] In a preferred embodiment, the line power flow constraint includes a line power constraint and a node voltage constraint;

[0102] Line power constraints:

[0103] |P l |≤P lmax ,|Q l |≤Q lmax (13);

[0104] Among them, P l is the active power of the lth branch, P lmax is the maximum active power of the lth branch; Q l is the reactive power of the lth branch, Q lmax is the maximum reactive power of the lth branch;

[0105] Node voltage constraints:

[0106] U k,min ≤U k ≤U k,max (14);

[0107] Among them, U k is the voltage of the kth node in the multi-energy park, U k,min is the minimum voltage of the kth node, Uk,max is the maximum voltage of the kth node.

[0108] In a preferred embodiment, the unit operation constraints include unit power constraints and unit ramp constraints;

[0109] Unit power constraints:

[0110] P jGkmin ≤P jGk,t ≤P jGkmax (15);

[0111] Among them, P jGk,t is the output power of the jth unit in the multi-energy park at time t, P jGkmin is the minimum output power of the jth unit, P jGkmax is the maximum output power of the jth unit;

[0112] The unit climbing constraint is:

[0113]

[0114] Among them, P jGk,t-1 is the output power of the jth unit at time t-1, is the maximum output power that can be adjusted downward by the j-th unit in the period [t-1, t], It is the maximum output power that can be adjusted upward by the j-th unit in the time period [t-1, t].

[0115] It is understandable that coal-fired units include unit power constraints and unit ramp constraints, while gas, photovoltaic and wind turbine units only have unit power constraints.

[0116] In a preferred embodiment, the multi-energy load balancing constraint is:

[0117]

[0118] Among them, L t is the load demand of the multi-energy park at time t; P t grid is the power transmitted from the power grid to the multi-energy park at time t, They are respectively the power generation capacities of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set in the multi-energy park at time t.

[0119] In a preferred embodiment, the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system is:

[0120]

[0121] Among them, E 0and E' are the carbon emissions of the multi-energy park before and after the construction of the energy storage system, which are composed of the carbon amount transmitted from the power grid to the multi-energy park and the carbon emissions of the units in the multi-energy park; P G0,t and P G0,t ' are the active power transmitted from the power grid to the multi-energy park at time t before and after the construction of the energy storage system; and are the output of the nth gas unit at time t in the multi-energy park before and after the construction of the energy storage system; e G0,t is the carbon potential of the power grid at time t; is the carbon potential of the nth gas generator set at time t; t∈[0,T], T is the preset period.

[0122] It is understandable that based on the node power balance, the output curves of the power grid and the gas-fired generator sets in different periods will be different before and after the multi-energy park is connected to energy storage, and therefore the corresponding carbon emissions will be different.

[0123] Based on the same inventive concept as the first embodiment, the second embodiment provides Figure 2 A multi-energy park low-carbon operation optimization device shown includes: a grid node carbon potential calculation module 21, which is used to calculate the grid node carbon potential of the multi-energy park according to the prediction curve of the power grid input to the multi-energy park and the carbon flow; a state judgment model construction module 22, which is used to construct a storage energy charging and discharging state judgment model according to the grid node carbon potential; an energy storage system operation construction module 23, which is used to construct an objective function with the minimum sum of power generation cost, carbon emission cost and energy storage operation cost, and solve the operation curves and energy storage charging and discharging curves of various types of generator sets in the multi-energy park in combination with the energy storage charging and discharging state judgment model and the constructed constraints; a carbon emission reduction calculation module 24, which is used to calculate the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system according to the operation curves of various types of generator sets and the energy storage charging and discharging curves.

[0124] In a preferred embodiment, the grid node carbon potential is:

[0125]

[0126] Among them, e k N is the carbon emission value equivalent to the power generation side caused by the unit power consumption of the kth node in the power grid; + is the set of all branches connected to the kth node that have carbon flow flowing into them; l is the branch number; R l is the carbon flow rate of the lth branch, F is the carbon flow of the lth branch, t is the time; P l is the injected power of the lth branch at the kth node.

[0127] In a preferred embodiment, the energy storage charging and discharging state judgment model is:

[0128]

[0129] Where M is of order of magnitude 10 4 A positive integer, y1 and y2 are both 0 / 1 variables, is the discharge carbon potential of the energy storage at time t, t∈[0,T], T is the preset period, e k,t is the charging carbon potential of the kth node in the multi-energy park at time t; when y1=0, y2=1, and When y1=1、y2=0, and When y1=y2=0, the energy storage is in the charging state. When the energy storage is off-grid.

[0130] In a preferred embodiment, the objective function is:

[0131]

[0132] Wherein, T is the preset period, t∈[0,T]; are the unit operating costs of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set in the multi-energy park, respectively. are the power generation of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set at time t; c emi For carbon tax, is the carbon potential of the nth gas-fired generator set at time t, C ESS is the unit operating cost of energy storage, P t is the charging or discharging power of the energy storage at time t.

[0133] In a preferred embodiment, the constraints include line flow constraints, unit operation constraints and multi-energy load balance constraints.

[0134] In a preferred embodiment, the line power flow constraint includes a line power constraint and a node voltage constraint;

[0135] Line power constraints:

[0136] |P l |≤P lmax ,|Q l |≤Q lmax (twenty two);

[0137] Among them, P l is the active power of the lth branch, P lmax is the maximum active power of the lth branch; Ql is the reactive power of the lth branch, Q lmax is the maximum reactive power of the lth branch;

[0138] Node voltage constraints:

[0139] U k,min ≤U k ≤U k,max (twenty three);

[0140] Among them, U k is the voltage of the kth node in the multi-energy park, U k,min is the minimum voltage of the kth node, U k,max is the maximum voltage of the kth node.

[0141] In a preferred embodiment, the unit operation constraints include unit power constraints and unit ramp constraints;

[0142] Unit power constraints:

[0143] P jGkmin ≤P jGk,t ≤P jGkmax (twenty four);

[0144] Among them, P jGk,t is the output power of the jth unit in the multi-energy park at time t, P jGkmin is the minimum output power of the jth unit, P jGkmax is the maximum output power of the jth unit;

[0145] The unit climbing constraint is:

[0146]

[0147] Among them, P jGk,t-1 is the output power of the jth unit at time t-1, is the maximum output power that can be adjusted downward by the j-th unit in the period [t-1, t], It is the maximum output power that can be adjusted upward by the j-th unit in the time period [t-1, t].

[0148] In a preferred embodiment, the multi-energy load balancing constraint is:

[0149]

[0150] Among them, L t is the load demand of the multi-energy park at time t; P t grid is the power transmitted from the power grid to the multi-energy park at time t, They are respectively the power generation capacities of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set in the multi-energy park at time t.

[0151] In a preferred embodiment, the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system is:

[0152]

[0153] Among them, E 0 and E' are the carbon emissions of the multi-energy park before and after the construction of the energy storage system, which are composed of the carbon amount transmitted from the power grid to the multi-energy park and the carbon emissions of the units in the multi-energy park; P G0,t and P G0,t ' are the active power transmitted from the power grid to the multi-energy park at time t before and after the construction of the energy storage system; and are the output of the nth gas unit at time t in the multi-energy park before and after the construction of the energy storage system; e G0,t is the carbon potential of the power grid at time t; is the carbon potential of the nth gas generator set at time t; t∈[0,T], T is the preset period.

[0154] In summary, the implementation of the embodiments of the present invention has the following beneficial effects:

[0155] The carbon potential of the grid nodes of the multi-energy park is calculated based on the prediction curve of the power grid's input of electricity and carbon flow to the multi-energy park; a storage energy charging and discharging state judgment model is constructed based on the carbon potential of the grid nodes; the objective function is constructed with the minimum sum of the power generation cost, carbon emission cost and energy storage operation cost, and the operation curves of various generator sets and energy storage charging and discharging curves in the multi-energy park are solved by combining the energy storage charging and discharging state judgment model and the constructed constraints; the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system is calculated based on the operation curves of various generator sets and energy storage charging and discharging curves, which can optimize the low-carbon operation of the multi-energy park and effectively reduce the carbon emissions of the multi-energy park.

[0156] The above is a preferred embodiment of the present invention. It should be pointed out that a person skilled in the art can make several improvements and modifications without departing from the principle of the present invention. These improvements and modifications are also considered to be within the scope of protection of the present invention.

[0157] Those skilled in the art can understand that all or part of the processes in the above embodiments can be implemented by instructing related hardware through a computer program, and the program can be stored in a computer-readable storage medium, and when the program is executed, it can include the processes in the above embodiments. The storage medium can be a disk, an optical disk, a read-only memory (ROM) or a random access memory (RAM), etc.

Claims

1. A low-carbon operation optimization method for a multi-energy park, characterized in that: include: According to the prediction curve of the power grid input electricity and carbon flow to the multi-energy park, the carbon potential of the grid node of the multi-energy park is calculated; the carbon potential of the grid node is: Among them, e k N is the equivalent carbon emission value on the power generation side caused by the unit electricity consumption of the kth node in the multi-energy park; + is the set of all branches connected to the kth node that have carbon flow flowing into them; l is the branch number; R l is the carbon flow rate of the lth branch, F is the carbon flow of the lth branch, t m is the unit time; P l is the injected power of the lth branch at the kth node; According to the carbon potential of the grid node, an energy storage charging and discharging state judgment model is constructed; the energy storage charging and discharging state judgment model is: Where M is of order of magnitude 10 4 A positive integer, y1 and y2 are both 0 / 1 variables, is the discharge carbon potential of the energy storage at time t, t∈[0,T], T is the preset period, e k,t is the charging carbon potential of the kth node in the multi-energy park at time t; when y1=0, y2=1, and When y1=1、y2=0, and When y1=y2=0, the energy storage is in the charging state. When , the energy storage is in an off-grid state; The objective function is constructed by minimizing the sum of power generation cost, carbon emission cost and energy storage operation cost. Combined with the energy storage charge and discharge state judgment model and the constructed constraints, the operation curves of various generator sets and energy storage charge and discharge curves in the multi-energy park are solved; the objective function is: Wherein, T is the preset period, t∈[0,T]; are the unit operating costs of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set in the multi-energy park, respectively. are the power generation of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set at time t; c emi For carbon tax, is the carbon potential of the nth gas-fired generator set at time t, C ESS is the unit operating cost of energy storage, P t is the charging or discharging power of the energy storage at time t; the constraints include line flow constraints, unit operation constraints and multi-energy load balance constraints; the line flow constraints include line power constraints and node voltage constraints; the line power constraints: |P l |≤P lmax ,|Q l |≤Q lmax ; Among them, P l is the active power of the lth branch, P lmax is the maximum active power of the lth branch; Q l is the reactive power of the lth branch, Q lmax is the maximum reactive power of the lth branch; The node voltage constraints are: IN k,min ≤U k ≤U k,max ; Among them, U k is the voltage of the kth node in the multi-energy park, U k,min is the minimum voltage of the kth node, U k,max is the maximum voltage of the kth node; the unit operation constraints include unit power constraints and unit climbing constraints; The unit power constraint is: P jGkmin ≤P jGk,t ≤P jGkmax ; Among them, P jGk,t is the output power of the jth unit in the multi-energy park at time t, P jGkmin is the minimum output power of the jth unit, P jGkmax is the maximum output power of the jth unit; The unit climbing constraint is: Among them, P jGk,t-1 is the output power of the jth unit at time t-1, is the maximum output power that can be adjusted downward by the j-th unit in the period [t-1, t], is the maximum output power that can be adjusted upward by the j-th unit in the period [t-1, t]; the multi-energy load balance constraint is: Among them, L t is the load demand of the multi-energy park at time t; P t grid is the power transmitted from the power grid to the multi-energy park at time t, are the power generation of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set in the multi-energy park at time t; the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system is: Among them, E 0 and E' are the carbon emissions of the multi-energy park before and after the construction of the energy storage system, which are composed of the carbon amount transmitted from the power grid to the multi-energy park and the carbon emissions of the units in the multi-energy park; P G0,t and P G0,t ' are the active power transmitted from the power grid to the multi-energy park at time t before and after the construction of the energy storage system; and are the output of the nth gas unit at time t in the multi-energy park before and after the construction of the energy storage system; e G0,t is the carbon potential of the power grid at time t; is the carbon potential of the nth gas generator set at time t; t∈[0,T], T is the preset period; According to the operating curves of the various types of generator sets and the energy storage charging and discharging curves, the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system is calculated.

2. A multi-energy park low-carbon operation optimization device, characterized in that: include: The grid node carbon potential calculation module is used to calculate the grid node carbon potential of the multi-energy park according to the prediction curve of the power grid input to the multi-energy park and the carbon flow; the grid node carbon potential is: Among them, e k N is the equivalent carbon emission value on the power generation side caused by the unit electricity consumption of the kth node in the multi-energy park; + is the set of all branches connected to the kth node that have carbon flow flowing into them; l is the branch number; R l is the carbon flow rate of the lth branch, F is the carbon flow of the lth branch, t m is the unit time; P l is the injected power of the lth branch at the kth node; The state judgment model construction module is used to construct an energy storage charge and discharge state judgment model according to the carbon potential of the grid node; the energy storage charge and discharge state judgment model is: Where M is of order of magnitude 10 4 A positive integer, y1 and y2 are both 0 / 1 variables, is the discharge carbon potential of the energy storage at time t, t∈[0,T], T is the preset period, e k,t is the charging carbon potential of the kth node in the multi-energy park at time t; when y1=0, y2=1, and When y1=1、y2=0, and When y1=y2=0, the energy storage is in the charging state. When , the energy storage is in an off-grid state; The energy storage system operation construction module is used to construct an objective function with the minimum sum of power generation cost, carbon emission cost and energy storage operation cost, and solve the operation curves of various generator sets and energy storage charge and discharge curves in the multi-energy park in combination with the energy storage charge and discharge state judgment model and the constructed constraints; the objective function is: Wherein, T is the preset period, t∈[0,T]; are the unit operating costs of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set in the multi-energy park, respectively. are the power generation of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set at time t; c emi For carbon tax, is the carbon potential of the nth gas-fired generator set at time t, C ESS is the unit operating cost of energy storage, P t is the charging or discharging power of the energy storage at time t; the constraints include line flow constraints, unit operation constraints and multi-energy load balance constraints; the line flow constraints include line power constraints and node voltage constraints; the line power constraints: |P l |≤P lmax ,|Q l |≤Q lmax ; Among them, P l is the active power of the lth branch, P lmax is the maximum active power of the lth branch; Q l is the reactive power of the lth branch, Q lmax is the maximum reactive power of the lth branch; The node voltage constraints are: IN k,min ≤U k ≤U k,max ; Among them, U k is the voltage of the kth node in the multi-energy park, U k,min is the minimum voltage of the kth node, U k,max is the maximum voltage of the kth node; the unit operation constraints include unit power constraints and unit climbing constraints; The unit power constraint is: P jGkmin ≤P jGk,t ≤P jGkmax ; Among them, P jGk,t is the output power of the jth unit in the multi-energy park at time t, P jGkmin is the minimum output power of the jth unit, P jGkmax is the maximum output power of the jth unit; The unit climbing constraint is: Among them, P jGk,t-1 is the output power of the jth unit at time t-1, is the maximum output power that can be adjusted downward by the j-th unit in the period [t-1, t], is the maximum output power that can be adjusted upward by the j-th unit in the period [t-1, t]; the multi-energy load balance constraint is: Among them, L t is the load demand of the multi-energy park at time t; P t grid is the power transmitted from the power grid to the multi-energy park at time t, are the power generation of the vth wind turbine generator set, the mth photovoltaic generator set, and the nth gas generator set in the multi-energy park at time t; the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system is: Among them, E 0 and E' are the carbon emissions of the multi-energy park before and after the construction of the energy storage system, which are composed of the carbon amount transmitted from the power grid to the multi-energy park and the carbon emissions of the units in the multi-energy park; P G0,t and P G0,t ' are the active power transmitted from the power grid to the multi-energy park at time t before and after the construction of the energy storage system; and are the output of the nth gas unit at time t in the multi-energy park before and after the construction of the energy storage system; e G0,t is the carbon potential of the power grid at time t; is the carbon potential of the nth gas generator set at time t; t∈[0,T], T is the preset period; The carbon emission reduction calculation module is used to calculate the carbon emission reduction of the multi-energy park before and after the construction of the energy storage system based on the operating curves of the various types of generator sets and the energy storage charging and discharging curves.

Images