Source network load storage integrated optimization operation method and system considering dynamic carbon emission factor

Abstract

According to the method, a power generation side and a power utilization side of a power grid are regarded as two macroscopic whole bodies, and a power generation side outer layer model and a power utilization side inner layer model considering dynamic carbon emission factors are respectively constructed through mutual restriction of industrial and commercial time-of-use electricity prices and the dynamic carbon emission factors; and reconstructing and solving the power generation side outer layer model and the power utilization side inner layer model which consider the dynamic carbon emission factor to obtain a source-grid-load-storage integrated double-layer optimization operation strategy which considers the dynamic carbon emission factor. According to the method, the master-slave game based on the electricity price and the carbon neutralization cost between the main network and the distribution network is considered, a traditional fixed carbon emission factor is improved into a time-varying dynamic carbon emission factor, the electricity utilization side flexible load is guided to deeply participate in the demand side response through the time-of-use electricity price and the dynamic carbon emission factor, the time sequence characteristics of power generation and electricity utilization are optimized, and the power utilization efficiency is improved. Therefore, the comprehensive power utilization cost of the power utilization side is reduced, and the power generation income of the power generation side is improved.

Description

Technical Field

[0001] The present invention belongs to the field of integrated optimization operation of power grid, load and storage, and specifically relates to an integrated optimization operation method and system of power grid, load and storage taking into account dynamic carbon emission factors. Background Art

[0002] With the continuous advancement of the construction of new power systems, the proportion of renewable energy generation in the global energy generation structure has gradually increased. The rapid development of clean energy such as wind power and solar energy has effectively reduced dependence on traditional fossil fuels and promoted the transformation of the power system towards a low-carbon and environmentally friendly direction. However, due to the intermittent and volatile nature of renewable energy generation, which poses a challenge to the stability of the power system, the balance and dispatch of the power grid face higher requirements.

[0003] As a key indicator to measure the carbon emission intensity in the electricity production process, the electricity carbon emission factor plays an important role in the carbon reduction process. By calculating the carbon emissions per unit of electricity, it can help decision makers identify the sources of high-emission electricity, thereby promoting the use of more low-carbon energy. However, with the increase in the proportion of renewable energy, the real-time changes in the load, energy structure and grid dispatch in the power system, the carbon emission intensity of electricity also fluctuates. The fixed electricity carbon emission factor can no longer accurately reflect the actual carbon emission changes in the power grid. If the carbon emission factor cannot be adjusted dynamically, the optimization of power generation and electricity consumption behavior will be limited. Therefore, improving the flexibility of the carbon emission factor and promoting more sophisticated carbon emission management are important tasks in the construction of future low-carbon power systems. Summary of the invention

[0004] The purpose of the present invention is to address the above-mentioned problems existing in the prior art and to provide a source-grid-load-storage integrated optimization operation method and system that takes into account dynamic carbon emission factors.

[0005] To achieve the above objectives, the technical solution of the present invention is as follows:

[0006] In the first aspect, the present invention proposes a source-grid-load-storage integrated optimization operation method considering dynamic carbon emission factors, including:

[0007] S1. The power generation side and the power consumption side of the power grid are regarded as two macroscopic entities. Through the mutual constraints of industrial and commercial time-of-use electricity prices and dynamic carbon emission factors, the outer model of the power generation side and the inner model of the power consumption side considering the dynamic carbon emission factors are constructed respectively;

[0008] S2. Reconstruct and solve the outer model of the power generation side and the inner model of the power consumption side considering the dynamic carbon emission factors, and obtain a two-layer optimized operation strategy of source, grid, load and storage considering the dynamic carbon emission factors.

[0009] The outer layer model of the power generation side in S1 aims to maximize the power generation profit, and establishes the objective function of the outer layer model of the power generation side as follows:

[0010] max I EC T -C S ;

[0011]

[0012] In the above formula, I E is the total power generation revenue of the grid, C T is the operating cost of the thermal power generating unit, C S is the startup cost of the thermal power generating unit, is the power transmitted from the power generation side to the power consumption side at time period t, T U is the unit time, is the electricity price in time period t, d is the number of typical days in a year, is the output power of the kth thermal power generating unit at time period t, PR F is the unit operating cost of the thermal power generating unit, δ T is the power generation efficiency of the thermal power generating unit, is the startup operation variable of the kth thermal power generating unit at time period t, PR S It is the single startup cost of the thermal power generating unit.

[0013] In S1, the constraints of the external model on the power generation side include:

[0014]

[0015]

[0016] In the above formula, is the startup operation variable of the kth thermal power generating unit at time period t, δ M is the large M constant, is the operating state variable of the kth thermal power generating unit at time period t, is the shutdown operation variable of the kth thermal power generating unit at time period t, α M , β M are the upper limit of the number of start-up and shutdown operations of a single thermal power generating unit in a typical day, TR t，k is the cumulative operating time of the kth thermal power generating unit in time period t, tx is the auxiliary time period variable, T LR is the minimum operating time of the thermal power generating unit, TS t,k is the cumulative downtime of the kth thermal power generating unit in time period t, T LS is the minimum downtime of the thermal power generating unit, are the lower and upper limits of the output of the kth thermal power generating unit, respectively. is the output power of the kth thermal power generating unit at time period t, are the upper limits of the down-slope and up-slope output of the kth thermal power generating unit, respectively. is the output power of the centralized photovoltaic system at time period t, is the system's spare capacity requirement at time period t, is the power transmitted from the power generation side to the power consumption side at time period t, T U is the unit time, is the dynamic carbon emission factor at time period t, ψ T is the carbon emission caused by the unit power generation of thermal power generating units, δ T It is the power generation efficiency of thermal power generating units.

[0017] The inner layer model of the power consumption side in S1 aims to minimize the comprehensive power consumption cost, and establishes the objective function of the inner layer model of the power consumption side as follows:

[0018] min X E +X C ;

[0019]

[0020] In the above formula, X E is the electricity purchase cost, X C To achieve carbon neutrality in electricity consumption, is the output power of the kth thermal power generating unit at time period t, is the output power of the centralized photovoltaic system at time period t, T U is the unit time, is the electricity price in time period t, d is the number of typical days in a year, is the power transmitted from the power generation side to the power consumption side at time period t, is the dynamic carbon emission factor at time period t, PR C is the unit carbon neutrality cost.

[0021] In S1, the constraints of the inner model on the power consumption side include:

[0022]

[0023] In the above formula, is the power transmitted from the power generation side to the power consumption side at time period t, is the base load of the station i at time period t, is the flexibility load of the station area i at time period t, is the energy storage output of station area i at time period t, is the distributed photovoltaic output of area i at time period t, is the real-time maximum power of distributed photovoltaic connected to the grid area i at time period t, is the lower limit of the discharge power of the energy storage system connected to the station area i, is the upper limit of the charging power of the energy storage system connected to the station area i, η L , η H are the lower and upper limits of the energy storage state, ES i is the installed capacity of the energy storage system connected to the grid area i, is the initial capacity of the energy storage system connected to the station area i, T U is the unit time, tn is any period in a typical day, is the flexible load energy demand of zone i in a typical day.

[0024] In the second aspect, the present invention proposes a source-grid-load-storage integrated optimization operation system considering dynamic carbon emission factors, including a model building module and a model solving module;

[0025] The model building module is used to regard the power generation side and the power consumption side of the power grid as two macroscopic entities, and to respectively build an outer model of the power generation side and an inner model of the power consumption side that consider the dynamic carbon emission factor through the mutual constraints of the industrial and commercial time-of-use electricity price and the dynamic carbon emission factor;

[0026] The model solving module is used to reconstruct and solve the outer model of the power generation side and the inner model of the power consumption side that consider the dynamic carbon emission factors, and obtain a two-layer optimized operation strategy of source, grid, load and storage that considers the dynamic carbon emission factors.

[0027] The model building module includes an objective function building unit of the outer model of the power generation side;

[0028] The objective function establishment unit of the power generation side outer layer model is used to establish the objective function of the power generation side outer layer model as follows:

[0029] max I E -C T -C S ;

[0030]

[0031]

[0032] In the above formula, I E is the total power generation revenue of the grid, C T is the operating cost of the thermal power generating unit, C S is the startup cost of the thermal power generating unit, is the power transmitted from the power generation side to the power consumption side at time period t, T U is the unit time, is the electricity price in time period t, d is the number of typical days in a year, is the output power of the kth thermal power generating unit at time period t, PR F is the unit operating cost of the thermal power generating unit, δ T is the power generation efficiency of the thermal power generating unit, is the startup operation variable of the kth thermal power generating unit at time period t, PR S It is the single startup cost of the thermal power generating unit.

[0033] The model building module also includes a constraint condition establishment unit for the outer model of the power generation side;

[0034] The constraint condition establishment unit of the power generation side outer layer model is used to establish the following constraint conditions of the power generation side outer layer model:

[0035]

[0036]

[0037] In the above formula, is the startup operation variable of the kth thermal power generating unit at time period t, δ M is the large M constant, is the operating state variable of the kth thermal power generating unit at time period t, is the shutdown operation variable of the kth thermal power generating unit at time period t, α M , β M are the upper limit of the number of start-up and shutdown operations of a single thermal power generating unit in a typical day, TR t,k is the cumulative operating time of the kth thermal power generating unit in time period t, tx is the auxiliary time period variable, T LR is the minimum operating time of the thermal power generating unit, TS t,k is the cumulative downtime of the kth thermal power generating unit in time period t, T LS is the minimum downtime of the thermal power generating unit, are the lower and upper limits of the output of the kth thermal power generating unit, respectively. is the output power of the kth thermal power generating unit at time period t, are the upper limits of the down-slope and up-slope output of the kth thermal power generating unit, respectively. is the output power of the centralized photovoltaic system at time period t, is the system's spare capacity requirement at time period t, is the power transmitted from the power generation side to the power consumption side at time period t, T U is the unit time, is the dynamic carbon emission factor at time period t, ψ Tis the carbon emission caused by the unit power generation of thermal power generating units, δ T It is the power generation efficiency of thermal power generating units.

[0038] The model building module also includes an objective function building unit for the inner layer model of the power consumption side;

[0039] The objective function establishment unit of the inner layer model of the power consumption side is used to establish the objective function of the inner layer model of the power consumption side as follows:

[0040] min X E +X C ;

[0041]

[0042] In the above formula, X E is the electricity purchase cost, X C To achieve carbon neutrality in electricity consumption, is the output power of the kth thermal power generating unit at time period t, is the output power of the centralized photovoltaic system at time period t, T U is the unit time, is the electricity price in time period t, d is the number of typical days in a year, is the power transmitted from the power generation side to the power consumption side at time period t, is the dynamic carbon emission factor at time period t, PR C is the unit carbon neutrality cost.

[0043] The model building module also includes a constraint condition establishment unit for the inner layer model of the power consumption side;

[0044] The constraint condition establishment unit of the inner layer model of the power consumption side is used to establish the following constraint conditions of the inner layer model of the power consumption side:

[0045]

[0046] In the above formula, is the power transmitted from the power generation side to the power consumption side at time period t, is the base load of the station i at time period t, is the flexibility load of the station area i at time period t, is the energy storage output of station area i at time period t, is the distributed photovoltaic output of area i at time period t, is the real-time maximum power of distributed photovoltaic connected to the grid area i at time period t, is the lower limit of the discharge power of the energy storage system connected to the station area i, is the upper limit of the charging power of the energy storage system connected to the station area i, η L , ηH are the lower and upper limits of the energy storage state, ES i is the installed capacity of the energy storage system connected to the grid area i, is the initial capacity of the energy storage system connected to the station area i, T U is the unit time, tn is any period in a typical day, is the flexible load energy demand of zone i in a typical day.

[0047] Compared with the prior art, the present invention has the following beneficial effects:

[0048] The present invention proposes a source-grid-load-storage integrated optimization operation method and system considering dynamic carbon emission factors. The method regards the power generation side and the power consumption side of the power grid as two macroscopic wholes, and through the mutual constraints of industrial and commercial time-of-use electricity prices and dynamic carbon emission factors, respectively constructs an outer model of the power generation side and an inner model of the power consumption side that consider dynamic carbon emission factors. By reconstructing and solving the outer model of the power generation side and the inner model of the power consumption side that consider dynamic carbon emission factors, a source-grid-load-storage integrated two-layer optimization operation strategy considering dynamic carbon emission factors is obtained. On the one hand, this method takes into account the time-series changes in the power generation and proportion of thermal power generation units and centralized photovoltaic power generation on the power generation side of the power grid system in various time periods, and improves the traditional fixed carbon emission factor into a time-varying dynamic carbon emission factor, which more accurately reflects the carbon emission characteristics in the power generation process and provides a new incentive method for the adjustment of flexible loads on the power consumption side; on the other hand, this method regards the main grid and distribution network, that is, the power generation side and the power consumption side, as two macro-entities, and considers the master-slave game between the two based on electricity prices and carbon neutrality costs. Through time-of-use electricity prices and dynamic carbon emission factors, the flexible loads on the power consumption side are guided to deeply participate in the demand-side response, make full use of centralized and distributed photovoltaic output, and optimize the time-series characteristics of power generation and power consumption, thereby reducing the comprehensive electricity consumption cost on the power consumption side while increasing the power generation income on the power generation side. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] Figure 1 This is a topological diagram of the power generation side and the power consumption side described in Example 1.

[0050] Figure 2 This is the load characteristic curve diagram described in Example 1.

[0051] Figure 3 This is the industrial and commercial time-of-use electricity price diagram described in Example 1.

[0052] Figure 4 This is the photovoltaic output characteristic curve diagram described in Example 1.

[0053] Figure 5 The figure is an overall flow chart of the method of the present invention.

[0054] Figure 6The structure diagram of the system of the present invention is shown in FIG. DETAILED DESCRIPTION

[0055] The present invention is further described in detail below in conjunction with specific implementations and drawings.

[0056] The present invention proposes an integrated optimization operation method and system for source-grid-load-storage considering dynamic carbon emission factors, constructs an outer model on the power generation side with the goal of maximizing power generation profits, and constructs an inner model on the power consumption side with the goal of minimizing the comprehensive electricity consumption costs. The inner and outer layers are mutually constrained by industrial and commercial time-of-use electricity prices and dynamic carbon emission factors, forming a master-slave game relationship, and the KKT conditional method is used to reconstruct the double-layer optimization problem of the outer model on the power generation side and the inner model on the power consumption side into a single-layer optimization problem, and the single-layer optimization problem model is solved to obtain a double-layer optimization operation strategy for source-grid-load-storage considering dynamic carbon emission factors. This method reduces the comprehensive electricity consumption costs on the power consumption side while improving the power generation income on the power generation side.

[0057] Embodiment 1:

[0058] This embodiment uses Figure 1 As an example of the topological structure of the power generation side and the power consumption side, the characteristic curves of the flexible load and the basic load are as follows: Figure 2 As shown in the figure, the industrial and commercial time-of-use electricity price is as follows Figure 3 As shown in Figure 2, the photovoltaic output characteristic curve is as follows: Figure 4 As shown. Parameters: the unit time is 1 hour, the investment period is 10 years, the typical days in a year are 365 days, the power generation efficiency of the thermal power generating unit is 45%, the unit operation cost of the thermal power generating unit is 0.4 yuan per kilowatt-hour, the single startup cost of the thermal power generating unit is 100,000 yuan, the number of thermal power generating units is 3, the maximum power of a single thermal power generating unit is 100 megawatts, the capacity of centralized photovoltaic is 100 megawatts, the minimum operating time of the thermal power generating unit is 4 hours, the minimum shutdown time of the thermal power generating unit is 4 hours, the lower and upper limits of the output of the thermal power generating unit are 30% and 100% of the maximum power of the unit, respectively, the upper limits of the ramp-down and ramp-up rates of the output of the thermal power generating unit are 2 megawatts and 5 megawatts, respectively, the carbon emissions generated by the unit power generation of the thermal power generating unit are 1.2 kilograms of carbon dioxide per kilowatt-hour, the unit carbon neutrality cost is 400 yuan per ton of carbon dioxide, and the lower and upper limits of the energy storage energy state are 10% and 90% of the power of the energy storage system, respectively.

[0059] like Figure 5 As shown in the figure, the source-grid-load-storage integrated optimization operation method considering the dynamic carbon emission factor is carried out in the following steps:

[0060] 1. The power generation side and the power consumption side of the power grid are regarded as two macroscopic entities. Through the mutual constraints of industrial and commercial time-of-use electricity prices and dynamic carbon emission factors, the outer model of the power generation side and the inner model of the power consumption side considering the dynamic carbon emission factors are constructed respectively;

[0061] The main grid and distribution grid of the power grid system, namely the power generation side and the power consumption side, are regarded as two macroscopic entities. The outer model of the power generation side aims to maximize the power generation profit. It takes into account the time series changes of the power and proportion of traditional thermal power generating units and centralized photovoltaic power generation in various time periods, and improves the traditional fixed carbon emission factor into a time-varying dynamic carbon emission factor, which more accurately reflects the carbon emission characteristics and provides a new incentive method for the adjustment of flexible loads on the power consumption side.

[0062] The inner model on the electricity consumption side aims to minimize the comprehensive electricity cost. It constrains the outer model on the power generation side through the industrial and commercial time-of-use electricity price and dynamic carbon emission factors. It takes into account the master-slave game relationship between the two based on electricity price and carbon neutrality cost, guides the flexible load on the electricity consumption side to deeply participate in the demand side response, makes full use of centralized and distributed photovoltaic output, and optimizes the timing characteristics of power generation and electricity consumption, thereby reducing the comprehensive electricity cost on the electricity consumption side while increasing the power generation income on the power generation side;

[0063] It should be noted that the present invention adopts the “consumer responsibility method” for apportionment, that is, the electricity user bears the carbon neutrality cost of electricity consumption;

[0064] The objective functions of the external model on the power generation side include:

[0065] max I E -C T -C S ;

[0066]

[0067] In the above formula, I E is the total power generation revenue of the grid, C T is the operating cost of the thermal power generating unit, C S is the startup cost of the thermal power generating unit, is the power transmitted from the power generation side to the power consumption side at time period t, T U is the unit time, is the electricity price in time period t, d is the number of typical days in a year, is the output power of the kth thermal power generating unit at time period t, PR F is the unit operating cost of the thermal power generating unit, δ T is the power generation efficiency of the thermal power generating unit, is the startup operation variable of the kth thermal power generating unit at time period t, PR S is the single startup cost of the thermal power generating unit;

[0068] The constraints of the external model on the power generation side include:

[0069] The startup and shutdown operation constraints of thermal power generating units are:

[0070]

[0071] The upper limit of the number of starts and stops of thermal power generating units is:

[0072]

[0073]

[0074] The minimum operating time constraint of the thermal power generating unit is:

[0075]

[0076] The minimum shutdown time constraint of the thermal power generating unit is:

[0077]

[0078] The output power constraint of the thermal power generating unit is:

[0079]

[0080] The climbing rate constraint of the thermal power generating unit is:

[0081]

[0082] The system's spare capacity constraint is:

[0083]

[0084] The power balance constraint for transmission to the distribution network is:

[0085]

[0086] The dynamic carbon emission factor and the actual carbon emission constraints are:

[0087]

[0088] In the above formula, is the startup operation variable of the kth thermal power generating unit at time period t. When it is 1, the unit has a startup operation, otherwise it does not. M is the large M constant, is the operating state variable of the kth thermal power generating unit at time period t. When it is 1, the unit is in operation, otherwise it is in shutdown state. is the shutdown operation variable of the kth thermal power generating unit at time period t. When it is 1, the unit has a shutdown operation, otherwise it does not. M , β M are the upper limit of the number of start-up and shutdown operations of a single thermal power generating unit in a typical day, TR t，k is the cumulative operating time of the kth thermal power generating unit in time period t, tx is the auxiliary time period variable, T LR is the minimum operating time of the thermal power generating unit, TS t，k is the cumulative downtime of the kth thermal power generating unit at time period t, where TR t，k and TS t，k The multiples of a single period are used for counting, T LS is the minimum downtime of the thermal power generating unit, are the lower and upper limits of the output of the kth thermal power generating unit, respectively. is the output power of the kth thermal power generating unit at time period t, are the upper limits of the down-slope and up-slope output of the kth thermal power generating unit, respectively. is the output power of the centralized photovoltaic system at time period t, is the system's spare capacity requirement at time period t, is the power transmitted from the power generation side to the power consumption side at time period t, T U is the unit time, is the dynamic carbon emission factor at time period t, ψ T is the carbon emission caused by the unit power generation of thermal power generating units, δ T The power generation efficiency of the thermal power generating unit;

[0089] The objective functions of the inner layer model on the electricity consumption side include:

[0090] min X E +X C ;

[0091]

[0092] In the above formula, X E is the electricity purchase cost, X C To achieve carbon neutrality in electricity consumption, is the output power of the kth thermal power generating unit at time period t, is the output power of the centralized photovoltaic system at time period t, T U is the unit time, is the electricity price in time period t, d is the number of typical days in a year, is the power transmitted from the power generation side to the power consumption side at time period t, is the dynamic carbon emission factor at time period t, PRC is the unit carbon neutrality cost;

[0093] The constraints of the inner model on the power consumption side include:

[0094] The power balance constraint on the power consumption side is:

[0095]

[0096] The actual output constraints of distributed photovoltaics connected to the distribution network area are:

[0097]

[0098] The energy storage operation constraints are:

[0099]

[0100] The energy demand constraint of the flexible load in a typical day is:

[0101]

[0102] In the above formula, is the power transmitted from the power generation side to the power consumption side at time period t, is the base load of the station i at time period t, is the flexible load of the station area i at time period t. The flexible load refers to the power load that can be met at any time of the day, such as the electric vehicle charging load. is the energy storage output of station area i at time period t, is the distributed photovoltaic output of area i at time period t, is the real-time maximum power of distributed photovoltaic connected to the grid area i at time period t, is the lower limit of the discharge power of the energy storage system connected to the station area i, is the upper limit of the charging power of the energy storage system connected to the station area i, η L , η H are the lower and upper limits of the energy storage state, ES i is the installed capacity of the energy storage system connected to the grid area i, is the initial capacity of the energy storage system connected to the station area i, T U is the unit time, tn is any period in a typical day, is the flexible load energy demand of zone i in a typical day.

[0103] 2. The KKT conditional method is used to reconstruct the double-layer optimization problem of the outer model on the power generation side and the inner model on the power consumption side into a single-layer optimization problem, and the single-layer optimization problem model is solved to obtain a double-layer optimization operation strategy of source-grid-load-storage integration that considers dynamic carbon emission factors;

[0104] The single-layer optimization problem model of the KKT conditional method substitutes the dynamic carbon emission factor and actual carbon emission constraint of the outer model on the power generation side and the power balance constraint of the inner model on the power consumption side into the objective function of the inner model on the power consumption side. The variables of the inner model on the power consumption side are:

[0105] The reconstructed single-layer optimization problem model takes the actual output constraints of distributed photovoltaics connected to the distribution network area of ​​the inner layer model on the power consumption side, the energy storage operation constraints, and the flexible load energy demand constraints within a typical day as the original feasible conditions;

[0106] The dual feasible conditions of the reconstructed single-layer optimization problem model include:

[0107]

[0108] In the above formula, v t，i , All are Lagrange multipliers;

[0109] The complementary relaxation conditions of the reconstructed single-layer optimization problem model include:

[0110]

[0111] The Big-M method is used to linearize the complementary relaxation condition and the following formula is obtained:

[0112] In the above formula, All are binary auxiliary variables.

[0113] In order to verify the effectiveness of the source-grid-load-storage integrated dual-layer optimization operation strategy considering dynamic carbon emission factors proposed in this scheme, referred to as Strategy 1, the source-grid-load-storage integrated dual-layer optimization operation strategy based on fixed carbon emission factors, referred to as Strategy 2, is used to compare with Strategy 1 proposed in this scheme. The comparison results of the economic indicators of the two strategies on the power generation side and the power consumption side are shown in Table 1 and Table 2 respectively:

[0114] Table 1 Economic comparison table of power generation side

[0115]

[0116] Table 2 Economic comparison of electricity consumption side

[0117]

[0118] Strategy 2 uses a fixed carbon emission factor, so the incentive for flexible loads on the user side to participate in demand-side response is relatively low. As can be seen from the above table, at the economic level of the power generation side, when Strategy 1 is adopted, the power generation income increases by 1.21%, the operating cost of thermal power generating units decreases by 7.83%, the startup cost of thermal power generating units decreases by 4.02%, and the profit of the power generation side increases by 6.67% compared with Strategy 2; at the economic level of the power consumption side, when Strategy 1 is adopted, the power purchase cost increases by 1.21%, the carbon neutrality cost decreases by 15.37%, and the comprehensive power consumption cost of the power consumption side decreases by 2.36%. The above results show that the source-grid-load-storage integrated dual-layer optimization operation strategy proposed in this scheme considering the dynamic carbon emission factor has achieved the effect of reducing the comprehensive power consumption cost of the power consumption side while increasing the power generation income of the power generation side.

[0119] Embodiment 2:

[0120] like Figure 6 As shown, the source-grid-load-storage integrated optimization operation system considering dynamic carbon emission factors includes a model building module and a model solving module;

[0121] The model building module is used to regard the power generation side and the power consumption side of the power grid as two macroscopic entities, and to respectively build an outer model of the power generation side and an inner model of the power consumption side that consider the dynamic carbon emission factor through the mutual constraints of the industrial and commercial time-of-use electricity price and the dynamic carbon emission factor;

[0122] The model solving module is used to reconstruct and solve the outer model of the power generation side and the inner model of the power consumption side that consider the dynamic carbon emission factors, and obtain a two-layer optimized operation strategy of source, grid, load and storage that considers the dynamic carbon emission factors.

[0123] The model building module includes an objective function building unit of the outer model of the power generation side;

[0124] The objective function establishment unit of the power generation side outer layer model is used to establish the objective function of the power generation side outer layer model as follows:

[0125] max I E -C T -C S ;

[0126]

[0127] In the above formula, I E is the total power generation revenue of the grid, C T is the operating cost of the thermal power generating unit, C S is the startup cost of the thermal power generating unit, is the power transmitted from the power generation side to the power consumption side at time period t, T U is the unit time, is the electricity price in time period t, d is the number of typical days in a year, is the output power of the kth thermal power generating unit at time period t, PR F is the unit operating cost of the thermal power generating unit, δ T is the power generation efficiency of the thermal power generating unit, is the startup operation variable of the kth thermal power generating unit at time period t, PR S It is the single startup cost of the thermal power generating unit.

[0128] The model building module also includes a constraint condition establishment unit for the outer model of the power generation side;

[0129] The constraint condition establishment unit of the power generation side outer layer model is used to establish the following constraint conditions of the power generation side outer layer model:

[0130]

[0131]

[0132] In the above formula, is the startup operation variable of the kth thermal power generating unit at time period t, δ M is the large M constant, is the operating state variable of the kth thermal power generating unit at time period t, is the shutdown operation variable of the kth thermal power generating unit at time period t, α M , β M are the upper limit of the number of start-up and shutdown operations of a single thermal power generating unit in a typical day, TR t，k is the cumulative operating time of the kth thermal power generating unit in time period t, tx is the auxiliary time period variable, T LR is the minimum operating time of the thermal power generating unit, TS t，k is the cumulative downtime of the kth thermal power generating unit in time period t, T LS is the minimum downtime of the thermal power generating unit, are the lower and upper limits of the output of the kth thermal power generating unit, respectively. is the output power of the kth thermal power generating unit at time period t, are the upper limits of the down-slope and up-slope output of the kth thermal power generating unit, respectively. is the output power of the centralized photovoltaic system at time period t, is the system's spare capacity requirement at time period t, is the power transmitted from the power generation side to the power consumption side at time period t, T U is the unit time, is the dynamic carbon emission factor at time period t, ψ Tis the carbon emission caused by the unit power generation of thermal power generating units, δ T It is the power generation efficiency of thermal power generating units.

[0133] The model building module also includes an objective function building unit for the inner layer model of the power consumption side;

[0134] The objective function establishment unit of the inner layer model of the power consumption side is used to establish the objective function of the inner layer model of the power consumption side as follows:

[0135] min X E +X C ;

[0136]

[0137]

[0138] In the above formula, X E is the electricity purchase cost, X C To achieve carbon neutrality in electricity consumption, is the output power of the kth thermal power generating unit at time period t, is the output power of the centralized photovoltaic system at time period t, T U is the unit time, is the electricity price in time period t, d is the number of typical days in a year, is the power transmitted from the power generation side to the power consumption side at time period t, is the dynamic carbon emission factor at time period t, PR C is the unit carbon neutrality cost.

[0139] The model building module also includes a constraint condition establishment unit for the inner layer model of the power consumption side;

[0140] The constraint condition establishment unit of the inner layer model of the power consumption side is used to establish the following constraint conditions of the inner layer model of the power consumption side:

[0141]

[0142] In the above formula, is the power transmitted from the power generation side to the power consumption side at time period t, is the base load of the station i at time period t, is the flexibility load of the station area i at time period t, is the energy storage output of station area i at time period t, is the distributed photovoltaic output of area i at time period t, is the real-time maximum power of distributed photovoltaic connected to the grid area i at time period t, is the lower limit of the discharge power of the energy storage system connected to the station area i, is the upper limit of the charging power of the energy storage system connected to the station area i, η L , η H are the lower and upper limits of the energy storage state, ES i is the installed capacity of the energy storage system connected to the grid area i, is the initial capacity of the energy storage system connected to the station area i, T U is the unit time, tn is any period in a typical day, is the flexible load energy demand of zone i in a typical day.

Claims

1. A source-grid-load-storage integrated optimization operation method considering dynamic carbon emission factors, characterized in that: The method comprises: S1. The power generation side and the power consumption side of the power grid are regarded as two macroscopic entities. Through the mutual constraints of industrial and commercial time-of-use electricity prices and dynamic carbon emission factors, the outer model of the power generation side and the inner model of the power consumption side considering the dynamic carbon emission factors are constructed respectively; S2. Reconstruct and solve the outer model of the power generation side and the inner model of the power consumption side considering the dynamic carbon emission factors, and obtain a two-layer optimized operation strategy of source, grid, load and storage considering the dynamic carbon emission factors.

2. The source-grid-load-storage integrated optimization operation method considering dynamic carbon emission factors according to claim 1 is characterized in that: The outer layer model of the power generation side in S1 aims to maximize the power generation profit, and establishes the objective function of the outer layer model of the power generation side as follows: max I E -C T -C S ; In the above formula, I E is the total power generation revenue of the grid, C T is the operating cost of the thermal power generating unit, C S is the startup cost of the thermal power generating unit, is the power transmitted from the power generation side to the power consumption side at time period t, T U is the unit time, is the electricity price in time period t, d is the number of typical days in a year, is the output power of the kth thermal power generating unit at time period t, PR F is the unit operating cost of the thermal power generating unit, δ T is the power generation efficiency of the thermal power generating unit, is the startup operation variable of the kth thermal power generating unit at time period t, PR S It is the single startup cost of the thermal power generating unit.

3. The source-grid-load-storage integrated optimization operation method considering dynamic carbon emission factors according to claim 1 is characterized in that: In S1, the constraints of the external model on the power generation side include: In the above formula, is the startup operation variable of the kth thermal power generating unit at time period t, δ M is the large M constant, is the operating state variable of the kth thermal power generating unit at time period t, is the shutdown operation variable of the kth thermal power generating unit at time period t, α M , β M are the upper limit of the number of start-up and shutdown operations of a single thermal power generating unit in a typical day, TR t，k is the cumulative operating time of the kth thermal power generating unit in time period t, tx is the auxiliary time period variable, T LR is the minimum operating time of the thermal power generating unit, TS t，k is the cumulative downtime of the kth thermal power generating unit in time period t, T LS is the minimum downtime of the thermal power generating unit, are the lower and upper limits of the output of the kth thermal power generating unit, respectively. is the output power of the kth thermal power generating unit at time period t, are the upper limits of the down-slope and up-slope output of the kth thermal power generating unit, respectively. is the output power of the centralized photovoltaic system at time period t, is the system's spare capacity requirement at time period t, is the power transmitted from the power generation side to the power consumption side at time period t, T U is the unit time, is the dynamic carbon emission factor at time period t, ψ T is the carbon emission caused by the unit power generation of thermal power generating units, δ T It is the power generation efficiency of thermal power generating units.

4. The source-grid-load-storage integrated optimization operation method considering dynamic carbon emission factors according to claim 1 is characterized in that: The inner layer model of the power consumption side in S1 aims to minimize the comprehensive power consumption cost, and establishes the objective function of the inner layer model of the power consumption side as follows: my X E +X C ; In the above formula, X E is the electricity purchase cost, X C To achieve carbon neutrality in electricity consumption, is the output power of the kth thermal power generating unit at time period t, is the output power of the centralized photovoltaic system at time period t, T U is the unit time, is the electricity price in time period t, d is the number of typical days in a year, is the power transmitted from the power generation side to the power consumption side at time period t, is the dynamic carbon emission factor at time period t, PR C is the unit carbon neutrality cost.

5. The source-grid-load-storage integrated optimization operation method considering dynamic carbon emission factors according to claim 1 is characterized in that: In S1, the constraints of the inner model on the power consumption side include: In the above formula, is the power transmitted from the power generation side to the power consumption side at time period t, is the base load of the station i at time period t, is the flexibility load of the station area i at time period t, is the energy storage output of station area i at time period t, is the distributed photovoltaic output of area i at time period t, is the real-time maximum power of distributed photovoltaic connected to the grid area i at time period t, is the lower limit of the discharge power of the energy storage system connected to the station area i, is the upper limit of the charging power of the energy storage system connected to the station area i, η L , η H are the lower and upper limits of the energy storage state, ES i is the installed capacity of the energy storage system connected to the grid area i, is the initial capacity of the energy storage system connected to the station area i, T U is the unit time, tn is any period in a typical day, is the flexible load energy demand of zone i in a typical day.

6. A source-grid-load-storage integrated optimization operation system considering dynamic carbon emission factors, characterized in that: The system includes a model building module and a model solving module; The model building module is used to regard the power generation side and the power consumption side of the power grid as two macroscopic entities, and to respectively build an outer model of the power generation side and an inner model of the power consumption side that consider the dynamic carbon emission factor through the mutual constraints of the industrial and commercial time-of-use electricity price and the dynamic carbon emission factor; The model solving module is used to reconstruct and solve the outer model of the power generation side and the inner model of the power consumption side that consider the dynamic carbon emission factors, and obtain a two-layer optimized operation strategy of source, grid, load and storage that considers the dynamic carbon emission factors.

7. The source-grid-load-storage integrated optimization operation system considering dynamic carbon emission factors according to claim 6 is characterized in that: The model building module includes an objective function building unit of the outer model of the power generation side; The objective function establishment unit of the power generation side outer layer model is used to establish the objective function of the power generation side outer layer model as follows: max I E -C T -C S ; In the above formula, I E is the total power generation revenue of the grid, C T is the operating cost of the thermal power generating unit, C S is the startup cost of the thermal power generating unit, is the power transmitted from the power generation side to the power consumption side at time period t, T U is the unit time, is the electricity price in time period t, d is the number of typical days in a year, is the output power of the kth thermal power generating unit at time period t, PR F is the unit operating cost of the thermal power generating unit, δ T is the power generation efficiency of the thermal power generating unit, is the startup operation variable of the kth thermal power generating unit at time period t, PR S It is the single startup cost of the thermal power generating unit.

8. The source-grid-load-storage integrated optimization operation system considering dynamic carbon emission factors according to claim 6 is characterized in that: The model building module also includes a constraint condition establishment unit for the outer model of the power generation side; The constraint condition establishment unit of the power generation side outer layer model is used to establish the following constraint conditions of the power generation side outer layer model: In the above formula, is the startup operation variable of the kth thermal power generating unit at time period t, δ M is the large M constant, is the operating state variable of the kth thermal power generating unit at time period t, is the shutdown operation variable of the kth thermal power generating unit at time period t, α M , β M are the upper limit of the number of start-up and shutdown operations of a single thermal power generating unit in a typical day, TR t，k is the cumulative operating time of the kth thermal power generating unit in time period t, tx is the auxiliary time period variable, T LR is the minimum operating time of the thermal power generating unit, TS t，k is the cumulative downtime of the kth thermal power generating unit in time period t, T LS is the minimum downtime of the thermal power generating unit, are the lower and upper limits of the output of the kth thermal power generating unit, respectively. is the output power of the kth thermal power generating unit at time period t, are the upper limits of the down-slope and up-slope output of the kth thermal power generating unit, respectively. is the output power of the centralized photovoltaic system at time period t, is the system's spare capacity requirement at time period t, is the power transmitted from the power generation side to the power consumption side at time period t, T U is the unit time, is the dynamic carbon emission factor at time period t, ψ T is the carbon emission caused by the unit power generation of thermal power generating units, δ T It is the power generation efficiency of thermal power generating units.

9. The source-grid-load-storage integrated optimization operation system considering dynamic carbon emission factors according to claim 6 is characterized in that: The model building module also includes an objective function building unit for the inner layer model of the power consumption side; The objective function establishment unit of the inner layer model of the power consumption side is used to establish the objective function of the inner layer model of the power consumption side as follows: my X E +X C ; In the above formula, X E is the electricity purchase cost, X C To achieve carbon neutrality in electricity consumption, is the output power of the kth thermal power generating unit at time period t, is the output power of the centralized photovoltaic system at time period t, T U is the unit time, is the electricity price in time period t, d is the number of typical days in a year, is the power transmitted from the power generation side to the power consumption side at time period t, is the dynamic carbon emission factor at time period t, PR C is the unit carbon neutrality cost.

10. The source-grid-load-storage integrated optimization operation system considering dynamic carbon emission factors according to claim 6, characterized in that: The model building module also includes a constraint condition establishment unit for the inner layer model of the power consumption side; The constraint condition establishment unit of the inner layer model of the power consumption side is used to establish the following constraint conditions of the inner layer model of the power consumption side: In the above formula, is the power transmitted from the power generation side to the power consumption side at time period t, is the base load of the station i at time period t, is the flexibility load of the station area i at time period t, is the energy storage output of station area i at time period t, is the distributed photovoltaic output of area i at time period t, is the real-time maximum power of distributed photovoltaic connected to the grid area i at time period t, is the lower limit of the discharge power of the energy storage system connected to the station area i, is the upper limit of the charging power of the energy storage system connected to the station area i, η L , η H are the lower and upper limits of the energy storage state, ES i is the installed capacity of the energy storage system connected to the grid area i, is the initial capacity of the energy storage system connected to the station area i, T U is the unit time, tn is any period in a typical day, is the flexible load energy demand of zone i in a typical day.

Images