CN117748449A

CN117748449A - Regional electric energy substitution distributed optimization method and system based on electric carbon transaction

Info

Publication number: CN117748449A
Application number: CN202311439440.2A
Authority: CN
Inventors: 张思瑞; 卜凡鹏; 张静; 成岭; 苗淑平; 高子寒; 李静; 李文; 刘伟; 李建锋; 郭京超; 李德智; 刘畅; 李阳; 王勇; 林晶怡; 王占博; 屈博; 蒋利民; 李梦
Original assignee: State Grid Corp of China SGCC; China Electric Power Research Institute Co Ltd CEPRI; Jining Power Supply Co
Current assignee: State Grid Corp of China SGCC; China Electric Power Research Institute Co Ltd CEPRI; Jining Power Supply Co
Priority date: 2023-11-01
Filing date: 2023-11-01
Publication date: 2024-03-22

Abstract

Regional electric energy substitution distributed optimization method and system based on electric carbon transaction, comprising: calculating based on the obtained set value of the electric quantity of the regional electric energy electric transaction in combination with a pre-constructed upper model to obtain the marginal electricity price of a distribution node of a node connected with each energy device in the region; the method comprises the steps that based on the marginal electricity price of a distribution node and a pre-built lower model, distributed solving is carried out by utilizing an alternating direction multiplying sub-distributed optimization algorithm to coordinate electricity transaction and carbon transaction among electric energy replacement of each energy device in an area, and electric quantity of the electric energy replacement electricity transaction of the area and output of each energy device are obtained; re-inputting the electric quantity of the regional electric energy replacement electric transaction into an upper model for iterative calculation until the optimal electric quantity of the regional electric energy replacement electric transaction and the optimal output of each energy device are obtained; the upper model constructed based on the double-layer interactive framework for trading regional electric energy substitution and the power distribution network can provide proper price signals for regional electric energy substitution, and guide reasonable optimization.

Description

Regional electric energy substitution distributed optimization method and system based on electric carbon transaction

Technical Field

The invention relates to the field of regional electric energy substitution distributed optimization methods, in particular to a regional electric energy substitution distributed optimization method and system based on electric carbon transaction.

Background

In recent years, the problem of climate change which is increasingly focused by the whole society is faced, the low-carbon economic development is more and more emphasized in the society nowadays, the carbon emission is reduced under the background, and the realization of economic sustainable development is one of the hot problems focused by the society. The recent energy price, double control index, supply and demand situation and other factors influence the electric energy replacement development to enter a new development stage. With the help of the output of the related documents of regional electric energy replacement, the regional low-carbon electric energy replacement practice is developed deeply on the premise of ensuring the safe operation and the reliable energy supply of the electric power system, a novel electric power system taking new energy as a main body is constructed, and the high-quality development of regional electric energy replacement is promoted. Regional power replacement is a typical representative of a novel power system, and carbon transaction cost is brought into the process of optimizing operation of regional power replacement, so that carbon emission of regional power replacement projects can be reduced, and good economical efficiency can be ensured. However, if the carbon emission reduction effect of the regional electric energy replacement project is unknown and the economical efficiency is poor in the electric power market and the carbon market environment, the construction and operation costs of the electric energy replacement project are high, and the investment enthusiasm of users is affected. In addition, if the regional electric energy substitution has the problems of unclear trading capacity and multiple trading coupling mechanisms, the regional electric energy substitution is lost in proportion of participating in electric and carbon trading, the regional electric energy substitution is difficult to deal with complex scenes in the trading operation process, and the problems of power trading interruption, unbalanced supply and demand and the like are easily caused. Therefore, the regional electric energy replacement distributed optimization characteristic based on the electric carbon transaction is enhanced, so that the regional electric energy replacement can realize interaction between the electric transaction and the carbon transaction and generate a new optimization plan, and the method has important significance in promoting on-site consumption of new energy, reducing carbon dioxide emission and promoting energy conservation and emission reduction of an electric power system.

Disclosure of Invention

In order to solve the problem of how to strengthen the regional electric energy substitution distributed optimization characteristic based on electric carbon transaction in the prior art, so that the regional electric energy substitution realizes the interaction of the electric transaction and the carbon transaction and generates a new optimization plan, the invention provides a regional electric energy substitution distributed optimization method based on the electric carbon transaction, which comprises the following steps:

calculating based on the obtained set value of the electric quantity of the regional electric energy electric transaction in combination with a pre-constructed upper model to obtain the marginal electricity price of a distribution node of a node connected with each energy device in the region;

based on the power distribution node marginal electricity price and a pre-built lower model, utilizing an alternating direction multiplication sub-distributed optimization algorithm to coordinate electric transaction and carbon transaction among electric energy replacement of each energy device in the area to carry out distributed solving, so as to obtain electric quantity of the electric energy replacement electric transaction of the area and output of each energy device;

re-inputting the electric quantity of the regional electric energy replacement electric transaction into the upper model for iterative calculation until the optimal electric quantity of the regional electric energy replacement electric transaction and the optimal output of each energy device are obtained;

the upper model is constructed based on a double-layer interactive framework for trading with the power distribution network by replacing regional electric energy containing each energy source device;

The lower model is constructed based on electric energy substitution of each energy source device in the region to perform electric transaction and carbon transaction.

Preferably, the method further comprises a construction process of an upper model, wherein the construction process comprises the following steps:

minimizing the power purchase cost of the power distribution network based on the price of the power distribution network for purchasing power from the power transmission line, the power generation cost of each energy device in the area, the power purchase quantity of the power distribution network from the power transmission line and the cost calculation of the transaction between each energy device in the area and the power distribution network;

the branch balance constraint is constructed by combining an active power balance principle and a reactive power balance principle based on the active power purchased by the power distribution network on the power transmission line, the reactive power purchased by the power distribution network on the power transmission line, the active power traded with the power distribution network in place of the regional power, the reactive power traded with the power distribution network in place of the regional power, the active power loaded by the node i in the t period, the reactive power loaded by the node i in the t period, the total active power loss of the power distribution network and the total reactive power loss of the power distribution network; wherein t is a positive integer;

the method comprises the steps of constructing an active power output constraint and a reactive power output constraint by combining a rated power upper limit principle based on active power and reactive power traded by each electric energy substitution with an upper-level power distribution network in a t period, an active power minimum value and a reactive power minimum value traded by each electric energy substitution with the upper-level power distribution network in the t period and an active power maximum value and a reactive power maximum value traded by each electric energy substitution with the upper-level power distribution network in the t period; wherein t is a positive integer;

Based on the maximum value of the node voltage of the power distribution network and the minimum value of the node voltage of the power distribution network, constructing and obtaining node voltage amplitude constraint by combining a kirchhoff current law and a kirchhoff voltage law;

and constructing an upper model based on the objective function, the branch balance constraint, the active output constraint, the reactive output constraint and the node voltage amplitude constraint.

Preferably, the method further comprises a construction process of the lower model, wherein the construction process comprises the following steps:

constructing an objective function of a lower model based on the sum minimization of the transaction cost of each electric energy replacement carbon in the area, the operation cost of the electric energy replacement in the area and the point-to-point transaction cost of each electric energy replacement in the area;

based on the gas turbine output, the photovoltaic power generation output, the power sold to the upper-level distribution network, the power purchased from the upper-level distribution network, the storage battery discharging power, the storage battery charging power and the electric quantity and load traded among the electric energy substitution projects in the area in a unit time period, and combining the active power balance construction to obtain a power balance constraint;

based on the rated power of the gas turbine and the output of the gas turbine in a unit time period, combining the upper limit principle of the rated power to construct and obtain the output constraint of the gas turbine;

based on the minimum climbing value of the gas turbine, the maximum climbing value of the gas turbine, the output of the gas turbine in the t period, the output of the gas turbine in the t-1 period and the output of the gas turbine at the initial moment, the climbing constraint of the gas turbine is obtained by combining the concept of dynamic programming; wherein t is a positive integer;

Based on the real-time energy of the storage battery, the charging power of the storage battery in a unit period, the charging efficiency of the storage battery, the discharging power of the storage battery in a unit period and the initial energy of the storage battery, and combining the upper limit of the capacity of the storage battery to construct and obtain the capacity constraint of the storage battery;

based on the charge and discharge state of the storage battery, the charge and discharge power of the storage battery and the capacity of the storage battery, the charge and discharge constraint of the storage battery is constructed by combining the charge and discharge principle of the storage battery;

based on the maximum value of the buying electric power and the maximum value of the selling electric power, combining a power interaction principle to construct and obtain an interaction power constraint;

based on the regional electric energy to replace the buying electric power, selling electric power, buying electric power maximum value and selling electric power maximum value, and combining a power interaction principle to construct and obtain interaction power constraint;

based on the minimum value and the maximum value of the photovoltaic power, combining the upper limit principle of rated power to construct and obtain the output constraint of the photovoltaic device;

and constructing a lower model based on the objective function, the power balance constraint, the gas turbine output constraint, the gas turbine climbing constraint, the storage battery capacity constraint, the storage battery charge and discharge constraint, the interactive power constraint and the photovoltaic device output constraint.

Preferably, the lower model further includes a carbon transaction cost calculation model, and the construction process of the carbon transaction cost calculation model includes:

calculating based on the output of the gas turbine and the thermal power unit in a unit time period and the carbon emission distribution of the unit electric quantity, so as to obtain carbon emission quota of each electric energy substitution in the area;

calculating based on the carbon emission intensity of unit active power output of the gas turbine and the thermal power unit and the output of the gas turbine and the thermal power unit in a unit time period to obtain the actual carbon emission amount of regional electric energy substitution;

calculating by using a step carbon price method based on carbon emission quota of each electric energy substitution in the area and actual carbon emission amount of each electric energy substitution in the area to obtain trade carbon price of the electric energy substitution carbon in the area;

and constructing a carbon trade cost calculation formula based on the carbon emission quota of each electric energy substitution in the area, the actual carbon emission amount of each electric energy substitution in the area and the carbon trade price of each electric energy substitution in the area.

Preferably, the calculation formula of the marginal electricity price of the distribution node is as follows:

in the method, in the process of the invention,the marginal electricity price of the active power distribution node is set; lambda (lambda) _P A first dual factor that balances constraints for the branches; p (P) _loss The total active loss of the distribution network is calculated; />The active power of the load of the node i at the moment t; lambda (lambda) _q A second dual factor that balances the constraint for the branch; q (Q) _loss The total reactive power loss of the power distribution network; />Reactive power of the load of the node i at the moment t; />A first dual factor that is a node voltage magnitude constraint; v (V) _j.t The node voltage of the node i at the time t is obtained; />And (3) a second dual factor for node voltage amplitude constraint, j is node ordering, and N is the number of nodes.

Preferably, the power distribution node marginal electricity price is combined with a pre-built lower model to coordinate electric transaction and carbon transaction among electric energy replacement of each energy device in the area by using an alternate direction multiplier sub-distributed optimization algorithm to perform distributed solution, so as to obtain electric quantity of the electric energy replacement electric transaction of the area and output of each energy device, including:

inputting the marginal electricity price of the distribution node as the price of each energy device for replacing the electricity transaction with the distribution network into the lower model for solving, and obtaining the optimal output and decision variables of each energy device;

wherein the decision variables include: the regional electric energy replaces the output power of each energy power generation device in the region, the electric quantity of the regional electric energy replacement electric transaction and the transaction quantity between each electric energy replacement in the region.

Preferably, the re-inputting the electric quantity of the regional electric energy replacement electric transaction into the upper model for iterative calculation until the optimal electric quantity of the regional electric energy replacement electric transaction and the optimal output of each energy device are obtained, including:

inputting the electric quantity of the regional electric energy replacement electric transaction into the upper model for calculation, and outputting a new power distribution node marginal electricity price;

inputting the new power distribution node marginal electricity price into the lower model for calculation, and outputting new regional electric energy to replace the electric quantity of electric transaction and new output of each energy device;

judging based on the new power distribution node marginal electricity price, and if the new power distribution node marginal electricity price is smaller than a set error value, taking the electric quantity of the new regional electric energy replacement electric transaction and the new energy equipment output as the optimal electric quantity of the regional electric energy replacement electric transaction and the optimal energy equipment output; otherwise, inputting the new regional electric energy to the upper model for replacing the electric quantity of the electric transaction, and continuing to iterate the calculation until the optimal regional electric energy to replace the electric quantity of the electric transaction and the optimal output of each energy device are obtained.

Based on the same inventive concept, the invention also provides a regional electric energy substitution distributed optimization system based on electric carbon transaction, which comprises the following steps:

The upper model calculation module is used for calculating based on the acquired set value of the electric quantity of the regional electric energy electric transaction and combining with a pre-constructed upper model to obtain the marginal electricity price of the distribution node of the node connected with each energy source device in the power distribution network and the region;

the lower model calculation module is used for carrying out distributed solution on electric transaction and carbon transaction among electric energy replacement of each energy device in the area by utilizing an alternating direction multiplication sub-distributed optimization algorithm based on the power distribution node marginal electricity price and a pre-built lower model to obtain electric quantity of the electric transaction of the regional electric energy replacement and output of each energy device;

the iteration solving module is used for inputting the electric quantity of the regional electric energy replacement electric transaction into the upper model again for iteration calculation until the optimal electric quantity of the regional electric energy replacement electric transaction and the optimal output of each energy device are obtained;

Preferably, the method further comprises an upper layer model building module, wherein the upper layer model building module is specifically used for:

Preferably, the method further comprises a lower model building module, wherein the lower model building module is specifically used for:

Preferably, the lower layer model further includes a carbon transaction cost calculation model construction sub-module, and the carbon transaction cost calculation model construction sub-module is specifically configured to:

Preferably, the calculation formula of the marginal electricity price of the distribution node in the upper model calculation module is as follows:

in the method, in the process of the invention,the marginal electricity price of the active power distribution node is set; lambda (lambda) _P A first dual factor that balances constraints for the branches; p (P) _loss The total active loss of the distribution network is calculated; />With load at time t for node iA power; lambda (lambda) _q A second dual factor that balances the constraint for the branch; q (Q) _loss The total reactive power loss of the power distribution network; />Reactive power of the load of the node i at the moment t; />A first dual factor that is a node voltage magnitude constraint; v (V) _j.t The node voltage of the node i at the time t is obtained; />And (3) a second dual factor for node voltage amplitude constraint, j is node ordering, and N is the number of nodes.

Preferably, the lower model calculation module is specifically configured to:

Preferably, the iterative solving module is specifically configured to:

In yet another aspect, the present invention also provides a computer device, including:

one or more processors;

a processor for executing one or more programs;

the one or more programs, when executed by the one or more processors, implement an regional power replacement distributed optimization method based on an electrical carbon transaction as described above.

In yet another aspect, the present invention further provides a computer readable storage medium having a computer program stored thereon, wherein the computer program, when executed, implements a regional power replacement distributed optimization method based on an electric carbon transaction as described above.

Compared with the prior art, the invention has the beneficial effects that:

regional electric energy substitution distributed optimization method and system based on electric carbon transaction, comprising: calculating based on the obtained set value of the electric quantity of the regional electric energy electric transaction in combination with a pre-constructed upper model to obtain the marginal electricity price of a distribution node of a node connected with each energy device in the region; based on the power distribution node marginal electricity price and a pre-built lower model, utilizing an alternating direction multiplication sub-distributed optimization algorithm to coordinate electric transaction and carbon transaction among electric energy replacement of each energy device in the area to carry out distributed solving, so as to obtain electric quantity of the electric energy replacement electric transaction of the area and output of each energy device; re-inputting the electric quantity of the regional electric energy replacement electric transaction into the upper model for iterative calculation until the optimal electric quantity of the regional electric energy replacement electric transaction and the optimal output of each energy device are obtained; the upper model is constructed based on a double-layer interactive framework for trading with the power distribution network by replacing regional electric energy containing each energy source device; the lower model is constructed by carrying out electric transaction and carbon transaction based on the replacement of electric energy of each energy device in the area; the upper model constructed based on the double-layer interactive framework for trading the regional electric energy substitution containing each energy source device with the power distribution network can provide proper price signals for the regional electric energy substitution, and guide the regional electric energy substitution to be reasonably optimized; the optimal scheduling plan can be calculated under the condition that the data privacy among the electric energy substitution of different main bodies can be ensured by adopting the alternating direction multiplication sub-distributed optimization algorithm.

Drawings

FIG. 1 is a flow chart of a regional power replacement distributed optimization method based on electric carbon transaction according to the present invention;

FIG. 2 is a block diagram of a two-layer optimization framework of regional power replacement and distribution network in consideration of electric-to-carbon transactions in accordance with the present invention;

FIG. 3 is a flow chart of an area power substitution optimization method taking into account electric carbon transactions according to the present invention;

FIG. 4 is a graph of the daily load curve and photovoltaic predicted force curve of the present invention versus grid DLMP;

FIG. 5 is a graph of the upper and lower iterative processes of calculating DLMP at the upper layer and solving regional power substitution trading volume at the lower layer;

FIG. 6 is a graph of the active power price of a distribution node after the convergence of the considered regional power substitution transaction amount and the DLMP optimization iteration of the present invention;

FIG. 7 is a graph of the power substitution (1) day-ahead optimization result after convergence of the consideration area power substitution transaction amount and the DLMP optimization iterations according to the present invention;

FIG. 8 is a graph of the power substitution (2) day-ahead optimization result after convergence of the consideration area power substitution transaction amount and the DLMP optimization iterations according to the present invention;

FIG. 9 is a graph of the power substitution (3) day-ahead optimization results after convergence of the consideration area power substitution transaction amount and the DLMP optimization iterations of the present invention;

FIG. 10 is a graph of the power substitution (1) day before optimization results after convergence of the power substitution trade amount and the DLMP optimization iterations without consideration of the region in the present invention.

Detailed Description

Under a plurality of regional electric energy substitution cooperative operation modes, carrying out trade of electric energy and carbon emission quota by regional electric energy substitution; the main body with rich electric energy resources can sell surplus electric energy to other main bodies in a Peer-to-Peer (P2P) mode; carbon trade quotas there are remaining subjects that can sell carbon quotas to carbon emission rights markets, while subjects that do not have sufficient carbon trade quotas need to purchase carbon quotas from carbon emission rights markets, minimizing regional power replacement operating costs and carbon trade costs. The research of the electric carbon collaborative operation mechanism and the interval optimization method of the electric energy replacement project based on carbon tracing is carried out, the operation of the electric energy replacement project is guided, the economy of the electric energy replacement project is improved, the regional electric energy replacement electric carbon transaction interaction is improved, the regional electric energy replacement carbon emission is reduced, the regional electric energy replacement new energy consumption is improved, and the regional electric energy replacement scheduling level is optimized, so that the method has a certain guiding significance and has an extremely important significance for realizing the double-carbon target. For a better understanding of the present invention, reference is made to the following description, drawings and examples.

Example 1:

an regional electric energy substitution distributed optimization method based on electric carbon transaction specifically comprises the following steps as shown in fig. 1:

Step 1, calculating based on a set value of the electric quantity of the acquired regional electric energy electric transaction in combination with a pre-constructed upper model to obtain a marginal electricity price of a distribution node of a node connected with each energy device in a power distribution network and a region;

step 2, carrying out distributed solving by utilizing an alternating direction multiplying sub-distributed optimization algorithm to coordinate electric transaction and carbon transaction among electric energy replacement of each energy device in the area based on the marginal electricity price of the distribution node and a pre-built lower model, so as to obtain the electric quantity of the electric transaction of the electric energy replacement of the area and the output of each energy device;

step 3, re-inputting the electric quantity of the regional electric energy replacement electric transaction into the upper model for iterative calculation until the optimal electric quantity of the regional electric energy replacement electric transaction and the optimal output of each energy device are obtained;

The method further comprises the step of constructing an upper layer model before the step 1, wherein the construction process of the upper layer model specifically comprises the following steps:

The invention provides a double-layer interactive framework of a distribution node marginal electricity price distribution locational marginal pricing based on local market trading, a power distribution network of DLMP and a regional electric energy replacement project, as shown in figure 2, the distribution node marginal electricity price can be regarded as expansion of the electricity price of a transmission node, a distributed power supply is guided to reasonably access the power distribution network through price signals, and therefore, how the power distribution network provides proper price signals for regional electric energy replacement is analyzed to guide regional electric energy replacement optimization. The electric power market operators conduct day-ahead market clearing to achieve maximum social benefit. The marginal electricity price of the distribution node for market clearing calculation by the upper model is the electricity transaction price between regional electric energy substitution and the distribution network

The power purchase cost of the power distribution network, which is obtained based on the price of the power distribution network for purchasing power from a power transmission line, the power generation cost of each energy device in an area, the power purchase quantity of the power distribution network from the power transmission line and the cost calculation of each energy device in the area for trading with the power distribution network, is minimized as an objective function:

wherein T is the number of study period time periods, T is the time period ordering, c _grid For the distribution network to purchase the active electricity prices of electricity from the transmission lines,active electric quantity purchased on power transmission line for t-period power distribution network, d _grid Reactive power price for purchasing power from transmission line for distribution network,/-for power distribution network>For the reactive power purchased by the power distribution network in the period t on the power transmission line, N is the number of nodes, i is the node ordering, and M is the zoneNumber of intra-domain electric energy substitutions, x is electric energy substitution order, c _i Active power cost for replacing power generation for each electric energy inside the area, +.>Substituting active power traded by the power distribution network for regional power connected by nodes of the power distribution network at t time intervals, d _i Reactive power cost for replacing power generation for each electric energy in the area, < >>And substituting reactive power transacted with the power distribution network for the regional power connected with the power distribution network node at the t period.

The branch balance constraint is constructed by combining an active power balance principle and a reactive power balance principle to obtain the branch balance constraint based on the active power purchased by the power distribution network on the power transmission line, the reactive power purchased by the power distribution network on the power transmission line, the active power traded with the power distribution network by regional power substitution, the reactive power traded with the power distribution network by regional power substitution, the active power loaded by the node i in the t period, the reactive power loaded by the node i in the t period, the total active power loss of the power distribution network and the total reactive power loss of the power distribution network:

wherein lambda is _p A first dual factor that balances the constraint for the branch, For the active electric quantity purchased by the power distribution network in the t period on the power transmission line, i is node ordering, nmg is the number of electric energy substitutions connected with the node i, M is the number of electric energy substitutions in the area, x is the order of electric energy substitutions, and +>Substituting the active power traded with the distribution network for the regional electric energy connected by the nodes of the distribution network at the t period,/->Active power for node i to load in t period, P _loss For the total active loss of the distribution network, S _x Lambda is the set of decision variables _q A second dual factor for the branch balancing constraint, < ->For the reactive power purchased by the power distribution network in the t period on the power transmission line, N is the number of nodes, and +.>Reactive power for node i load in t period, Q _loss The total reactive power loss of the distribution network is obtained.

The method comprises the steps of constructing an active output constraint and a reactive output constraint by combining a rated power upper limit principle based on active power and reactive power traded by each electric energy substitution with an upper-level power distribution network in a t period, an active power minimum value and a reactive power minimum value traded by each electric energy substitution with the upper-level power distribution network in the t period and an active power maximum value and a reactive power maximum value traded by each electric energy substitution with the upper-level power distribution network in the t period:

in the method, in the process of the invention,first dual factor for active force constraint, +.>A second dual factor that is the active force constraint, Is a nodeEach electric energy in the connected area replaces the minimum value of active power traded with the upper-level distribution network in the period t,substituting active power traded with the upper level distribution network in t period for each electric energy in the area connected with the node, < >>And substituting each electric energy in the area connected with the node for the maximum active power which is traded with the upper-level distribution network in the period t.

The reactive power output constraint is calculated as follows:

in the method, in the process of the invention,first dual factor for reactive force constraint, < ->A second dual factor that is reactive force constraint,the reactive power minimum value traded with the upper-level distribution network in the period t is replaced by each electric energy in the area connected with the node,substituting reactive power traded with the upper level distribution network in period t for each electric energy in the area connected with the node, +.>And substituting the maximum reactive power value traded with the upper-level distribution network in the period t for each electric energy in the area connected with the node.

Based on the maximum value of the node voltage of the power distribution network and the minimum value of the node voltage of the power distribution network, the node voltage amplitude constraint is constructed by combining the kirchhoff current law and the kirchhoff voltage law:

in the method, in the process of the invention,first dual factor for node voltage amplitude constraint, +.>A second dual factor, V, being the node voltage magnitude constraint _i ^min Minimum node voltage for node i, V _i,t Node voltage of node i, S _x For the decision variable set, V _i ^max Is the maximum node voltage of node i.

Before the step 1, the method further comprises the step of constructing a lower model, wherein the construction process of the lower model specifically comprises the following steps:

regional power replacement lower layer models comprising gas turbines, photovoltaic power generation, storage batteries and loads are constructed, and a carbon transaction mechanism is added into the regional power replacement lower layer models, and the total cost is minimized as an objective function.

An objective function for constructing an underlying model based on the sum of the cost of each electrical energy replacement carbon transaction in the region, the cost of regional electrical energy replacement operation, and the cost of each electrical energy replacement point-to-point transaction in the region is minimized:

P1：

wherein F is _x As a total cost of the product,replacement of operating costs for regional power,/-)>Substitution of point-to-point transaction costs for each electrical energy in the area, < >>Replacing the electricity transaction cost for each electrical energy in the area.

Based on gas turbine output, photovoltaic power generation output, power sold to the upper level distribution grid, power purchased from the upper level distribution grid, battery discharge power, battery charging power, and power and load traded between the electric energy replacement projects in the area in a unit time period, a power balance constraint is obtained by combining active power balance construction:

Wherein j is the node order, N is the node number, i is the node order,power purchased for unit period to node j for replacing x for electric energy, +.>Substitution of the gas turbine unit time period output in x for node i electric energy, +.>Power output for node i power instead of photovoltaic unit time period in x, < >>For the power of discharging of the storage battery in unit time period in x of node i power substitution, y is power substitution order, x is power substitution order, M is the number of power substitution in the area, < >>When the unit of the regional power is replaced by x and the unit of the regional power is replaced by yBlock electric transaction amount, < >>Power sold to node j for power replacement x unit periods, +.>Substituting node i for battery unit period charging power in x, < >>The load in the power substitution of the inode.

Based on rated power of gas turbine and output of gas turbine in unit time period, and combined with upper limit principle of rated power, the gas turbine output constraint is constructed

The gas turbine output constraint is calculated as follows:

wherein P is _gn,x For the rated power of the gas turbine,the output is output per unit time period of the gas turbine in the inode electric energy substitution.

Based on the minimum climbing value of the gas turbine, the maximum climbing value of the gas turbine, the output of the gas turbine at the t period, the output of the gas turbine at the t-1 period and the output of the gas turbine at the initial moment, the climbing constraint of the gas turbine is obtained by combining the concept of dynamic programming:

When t is not 1, the number of times,

when the t is a number of times 1,

in vmt ^min Is the minimum climbing value of the gas turbine,substituting the power output of the internal gas turbine in the t period for each power of the area connected with the node i, +.>Power output of internal gas turbine at t-1 period for each power substitution in the region connected to node i vmt ^max For maximum climbing value of gas turbine, P _g0,x The regional power is used to replace the initial output of the internal gas turbine.

Based on the real-time energy of the storage battery, the charging power of the storage battery in unit time period, the charging efficiency of the storage battery, the discharging power of the storage battery in unit time period and the initial energy of the storage battery, and combining the upper limit of the storage battery capacity to construct and obtain the capacity constraint of the storage battery:

when t is not 1 =it is,

when the t is a number of times 1,

wherein ES1 _t ES1 is the real-time energy of the accumulator in the t period _t-1 Is the real-time energy of the storage battery in the t-1 period,substituting the power of the node i for the charging power of the storage battery in x in unit time period, eta _{_c} For battery charging efficiency, eta_d is battery discharging efficiency, < >>Substituting node i for electric energy in xES0 is the initial energy of the battery.

Based on the charge and discharge state of the storage battery, the charge and discharge power of the storage battery and the capacity of the storage battery, the charge and discharge constraint of the storage battery is constructed by combining the charge and discharge principle of the storage battery:

X _t ×Y _t ＝0

In the method, in the process of the invention,the charging power of the storage battery in unit time period in X is replaced by the electric energy _t Eb is the capacity of the battery for the battery state of charge, < > for the battery>Charging and discharging power of storage battery unit time period in x is replaced by electric energy, Y _t Is in a discharge state of the storage battery. />

Based on regional electric energy instead of buying electric power, selling electric power, buying electric power maximum value and selling electric power maximum value, and combining a power interaction principle to construct and obtain interaction power constraint:

in the method, in the process of the invention,substituting regional electric energy for buying electric power, +.>Instead of buying the maximum value of the electric power for the regional electric energy,replacing the power for the regional power, +.>The maximum value of electric power is replaced by the regional electric energy.

Based on the minimum value and the maximum value of the interaction power, the interaction power constraint is constructed by combining the power interaction principle:

in the method, in the process of the invention,maximum electric transaction amount for unit time period of electric energy substitution x and electric energy substitution y, +.>For the electric trade quantity of the unit time interval of the electric energy substitution x and the electric energy substitution y, < >>And replacing the electric transaction amount of x unit time periods for the electric energy by using the electric energy.

Based on the minimum photovoltaic power value and the maximum photovoltaic power value, the output constraint of the photovoltaic device is constructed and obtained by combining the rated power upper limit principle:

in the method, in the process of the invention,for the node i electric energy to replace the discharge power of the storage battery in x in unit time period, P _pv,max Is the maximum output value of the photovoltaic device.

The calculation formula of the regional power replacement operation cost is as follows:

in the method, in the process of the invention,to replace the running cost for regional electric energy, C _v The regional electric energy replaces the running cost of the internal photovoltaic power generation,c, substituting regional electric energy for internal photovoltaic power generation output in t period _b,t Time-of-use electricity price purchased from the upper distribution grid is replaced by regional electricity>Purchasing electric quantity from power distribution network for replacing x unit time period by electric energy, C _g Substitution of regional power for internal gas turbine operating costs, < >>To replace the output of the gas turbine in x in unit time period, C _x To replace the charge and discharge cost of the internal storage battery for regional electric energy, < >>C, substituting the regional electric energy for the discharge power of the internal storage battery in the t period _s,t Time-of-use electricity price for regional electricity to replace electricity sold to upper distribution grid, < >>Selling electric quantity to a distribution network for replacing x unit time periods by electric energy, < >>And replacing the charging power of the internal storage battery t period for regional electric energy.

The calculation formula of the point-to-point transaction cost of each electric energy in the area is as follows:

/>

in the method, in the process of the invention,for each electric energy replacement point-to-point transaction cost in the area, y is electric energy replacement ordering, x is electric energy replacement ordering, M is electric energy replacement number in the area, T is research period time period number, T is time period ordering, and>trade price for electricity of unit time period of electricity substitution x and y for electricity substitution +.>And substituting the electric energy for x and y for the electric trade amount of the unit time period.

The transaction cost of each electric energy replacing carbon in the area is calculated according to a carbon transaction cost calculation model; the method also comprises the construction of a carbon transaction cost calculation model, and the construction process specifically comprises the following steps:

the present embodiment recognizes that the carbon transaction operator uniformly manages carbon transactions for each electrical energy replacement in the area. Each electric energy replaces the electric power purchased from the distribution network to be generated by the thermal power generating unit. Thermal power generating units and gas turbines produce carbon emissions during power generation. And its carbon emission allowance was determined using a baseline method. In this model, each electrical energy replacement may purchase or sell carbon emission credits directly to the carbon trade operator according to its own needs to achieve its own maximization of benefits. The gratuitous carbon quota and carbon emissions for each electrical energy substitution can be expressed as:

Each electric energy in the region replaces gratuitous carbon quota:

wherein Y is _MG(x) For the carbon emission quota of each electric energy substitution in the region, T is the number of research period time periods, T is the time period sequencing,carbon emission allocation rate for industry uniform unit electric quantity, < >>Respectively, t is the active power of the gas turbine in the period of t,>the active power output of the thermal power generating unit in the t period is given, and delta t is the variation of the period.

Actual carbon emissions of regional power substitution:

wherein E is _MG(x) For the actual carbon emission of each electric energy substitution in the region, T is the number of research period time periods, T is the time period sequencing, a _gen Carbon emission intensity for the unit active power output of the gas turbine,respectively, the active power output of the gas turbine in the period t, a _in Carbon emission intensity of unit active power output of thermal power generating unit,/-for>The active power output of the thermal power generating unit in the t period.

The regional electric energy replacement carbon trade price is calculated by adopting a step carbon price method, and the regional electric energy replacement carbon trade price is defined as follows:

wherein C is _c For the price of the carbon trade,is the basic carbon price of the market; beta is the price increasing rate of each step carbon; a is the length of a carbon emission interval, Y _MG(x) Carbon emission allowance for each electric energy substitution in the region, E _MG(x) Actual carbon emissions for each electrical energy substitution in the region.

If each electric energy in the region replaces the actual carbon emission E _MG(x) Above its carbon emission allowance Y _MG(x) When the electric energy in the area is replaced, the carbon transaction cost is paid to obtain additional carbon emission quota, otherwise, if the electric energy in the area is replaced with the actual carbon emission E _MG(x) Below its carbon emission allowance Y _MG(x) When the electric energy replacement in the area can sell carbon quota to obtain benefits through carbon transaction, the calculation formula of the electric energy replacement carbon transaction cost in the area is as follows:

in the method, in the process of the invention,for the trade cost of substituting carbon for each electric energy in the area, T is the number of research period time, T is the time period sequencing, C _c Price for carbon trade->For replacing each electric energy in t time period regionActual carbon emissions of the generation, +.>Carbon emission quota for each electrical energy substitution in the t-slot region.

The lower model is regional electric energy substitution containing gas turbines, photovoltaic power generation, storage batteries and loads respectively, carbon transaction is introduced into the regional electric energy substitution, P2P point-to-point transaction can be directly carried out on each regional electric energy substitution, transaction prices are negotiated and determined by two transaction parties, and intervention of a power distribution network operator is not needed. Through the cooperative transaction of a plurality of electric energy substitutes in the area, clean energy utilization on the source side and green energy utilization behavior guidance on the load side can be effectively promoted, emission reduction potential on the user side is further exerted, the win-win situation of multiparty benefit bodies is realized, various flexible resources are stimulated to eliminate new energy, benign development of the new energy is promoted, and low-carbon operation in the area is further realized.

In step 1, calculating based on the obtained set value of the electric quantity of the regional electric energy electric transaction in combination with a pre-constructed upper model to obtain the marginal electricity price of the distribution node of the node connected with each energy device in the power distribution network and the region, specifically comprising:

the upper model comprehensively considers the voltage constraint, the network constraint and the power balance constraint of the power distribution network, and calculates the marginal electricity price of the power distribution node as the electricity price of regional electric energy replacement electricity transaction with the aim of minimizing the cost of the power distribution network.

By definition of the node marginal electricity price, the reactive power distribution node marginal electricity price is not considered. And carrying out linearization treatment on node voltage, power loss and branch power flow of the power distribution network, wherein the calculation formula of the marginal electricity price of the active power distribution node is as follows:

in the method, in the process of the invention,the marginal electricity price of the active power distribution node is set; lambda (lambda) _P A first dual factor that balances constraints for the branches; p (P) _loss The total active loss of the distribution network is calculated; />The active power of the load of the node i at the moment t; lambda (lambda) _q A second dual factor that balances the constraint for the branch; q (Q) _loss The total reactive power loss of the power distribution network; />Reactive power of the load of the node i at the moment t; />A first dual factor that is a node voltage magnitude constraint; v (V) _j.t The node voltage of the node i at the time t is obtained; / >And (3) a second dual factor for node voltage amplitude constraint, j is node ordering, and N is the number of nodes.

In step 2, the power transaction and the carbon transaction between the electric energy substitution of each energy device in the area are coordinated by using an alternating direction multiplier sub-distributed optimization algorithm based on the marginal electricity price of the distribution node and a pre-built lower model to perform distributed solving, so as to obtain the electric quantity of the electric energy substitution electric transaction of the area and the output of each energy device, and the method specifically comprises the following steps:

the marginal electricity price of the distribution node of the connection node of each energy device in the distribution network and the area calculated by the upper model is the price of the electricity transaction between the replacement of the electric energy of each energy device in the lower model and the distribution network, namely:

wherein C is _b,t The price of electric trade with the distribution network is replaced for the electric energy of each energy source device,and replacing the marginal electricity prices of the distribution nodes of each node connected with the power distribution network and the internal electric energy of the area.

And taking the marginal electricity price of the distribution node as the price of each energy device for replacing the electricity transaction with the distribution network to be input into the lower model, carrying out distributed solving by adopting an ADMM algorithm by using a matlab solver, wherein the optimal output of each energy device and the electric quantity of the regional electricity replacement electricity transaction are the output of the lower model, and the initial optimal control of each device in the power plant and the electric quantity of the electricity transaction can be obtained.

Decision variables include: the regional electric energy replaces the output power of each energy power generation device in the region, the electric quantity of the regional electric energy replacement electric transaction and the transaction quantity between each electric energy replacement in the region.

The regional power substitution system constructed by the invention contains a plurality of power substitution models, wherein different main bodies freely trade carbon emission and power in the models, and the regional power substitution system belongs to the multi-block optimization problem. The optimization problem is difficult to solve directly due to coupling constraints. The method separates coupling constraint, and utilizes the multiplication sub-distributed optimization algorithm based on the alternating direction to solve the model on the premise of guaranteeing data privacy between the electric energy substitution in different main body areas.

The solving steps of the alternating direction multiplication sub-distributed optimization algorithm are as follows:

s1: and reading in parameters and load requirements of each electric energy replacing internal distributed generator set, and entering S2.

S2: initializing an initial iteration index k=1 of an alternating direction multiplier sub-distributed optimization algorithm, and the maximum iteration number k _max =50, initial value of the augmented lagrangian multiplier Quadratic term penalty function initial value->Algorithm convergence accuracy epsilon=0.001,and (3) entering S3.

S3: solving by using the following calculation formula, and entering into S4;

in the method, in the process of the invention,replacement of operating costs for regional power,/-)>For regional electric energy replacement operation cost, M is the number of electric energy replacement in the region, x is the electric energy replacement sequence, T is the study period time period number, T is the time period sequence, and +.>For the trade marginal price between energy replacement x and energy replacement y, +.>For the electric trade quantity of the kth iteration unit time period of the electric energy substitution x and the electric energy substitution y,to augment Lagrangian quadratic term penalty function, +.>And replacing the electric transaction amount of the unit time period for the kth iteration of the y and x for the electric energy.

S4: each electric energy is replaced to mutually transmit expected transaction power according to the self optimization result, whether the expected transaction power is smaller than a preset value or not is judged, if so, stopping is carried out, and if not, S5 is carried out;

s5: setting k ₁ Returning to S3 after=k+1, until the desired trade power is greater than or equal to the preset value.

The upper model and the lower model are coupled, and the alternating direction multiplier sub-distributed optimization algorithm can coordinate privacy among the electric energy substitutes of different areas and coordinate electric transaction and carbon transaction among the electric energy substitutes of each area to obtain an optimal scheduling plan of the model.

The regional power replaces the power of the power transaction:

in the method, in the process of the invention,replacing the electric quantity of an electric transaction for regional electric energy, +.>For generating power of photovoltaic power within t period, < >>For t period of active power of the gas turbine, < >>For the discharge power of the accumulator in the t period, +.>For load in t period,/>And charging power of the storage battery in the t period.

In step 3, the electric quantity of the regional electric energy replacement electric transaction is input into the upper model again for iterative computation until the optimal electric quantity of the regional electric energy replacement electric transaction and the optimal output of each energy device are obtained, which specifically comprises the following steps:

inputting the electric quantity of the regional electric energy replacement electric transaction into the upper model for calculation, and outputting a new power distribution node marginal electricity price:

in order to strengthen the coupling of regional electric energy to replace electric carbon trade, an iteration mechanism is needed to coordinate the framework of electric trade and carbon trade, and the electric quantity of the regional electric energy to replace electric trade is used as a parameter electric power reference value and an adjustment quantity to be transmitted to an upper model according to a scheduling plan of each energy source device in the regional electric energy to replace electric trade, so that a new marginal electricity price of a distribution node is calculated and generated by the distribution network.

Inputting the new power distribution node marginal electricity price into the lower model for calculation, and outputting new regional electric energy to replace the electric quantity of electric transaction and new output of each energy device:

And the marginal electricity prices of the distribution nodes of the new distribution network and the area electric energy substitution connection nodes are used as parameters to be input into a lower model, and the electric quantity of the new area electric energy substitution electric transaction and the output of each new energy source device are obtained through re-solving.

Judging based on the new power distribution node marginal electricity price, and if the new power distribution node marginal electricity price is smaller than a set error value, taking the electric quantity of the new regional electric energy replacement electric transaction and the new energy equipment output as the optimal electric quantity of the regional electric energy replacement electric transaction and the optimal energy equipment output; otherwise, inputting the new regional electric energy to the upper model for iteration calculation until the optimal regional electric energy to replace the electric energy of the electric transaction and the optimal output of each energy device are obtained;

if the new marginal electricity price of the distribution node is larger than or equal to the set error value, continuously inputting the electric quantity of the new regional electric energy replacement electric transaction into the upper model for calculation until the marginal electricity price of the distribution node obtained by two continuous iterations is smaller than the allowable error value, and considering iteration convergence to obtain the optimal solution meeting the double-layer model.

Compared with the prior art, the invention provides a distribution network and regional electric energy replacing two-layer interactive framework based on the marginal electricity price of the distribution node;

The framework models the problem that regional electric energy substitution containing various energy devices participates in centralized bidding of the power distribution network layer, so that the power distribution network resource allocation can be fully guided, a proper price signal is provided for regional electric energy substitution, and reasonable optimization is guided.

The invention also constructs a lower model for combining regional electric energy replacement electric transaction and carbon transaction;

the model is a combination point of regional electric energy replacement electric transaction and carbon transaction, initial scheduling obtained by the model is solved, electric quantity of the regional electric energy replacement electric transaction is transmitted into an upper model of an interactive framework of the power distribution network and the regional electric energy replacement as parameters, and an optimal solution meeting regional electric energy replacement optimization is obtained through loop iteration. The invention comprehensively considers the transaction modes of the electric market and the carbon market, reflects that the change of the electricity price can influence the optimization control of each energy device in the regional electric energy replacement electricity and carbon transaction, and otherwise, the optimization of each energy device in the regional electric energy replacement can influence the change of the electricity price. So as to achieve the purposes of reducing the carbon emission of regional electric energy, improving new energy consumption and promoting the clean development of the power industry.

The iteration mechanism adopted by the invention can provide support for dynamic adjustment operation of regional electric energy replacement internal energy equipment in different time periods, so that flexible interaction among the regional electric energy replacement different energy equipment, double-layer distributed optimization of electric carbon transaction and regional electric energy replacement operation, optimal control of each energy equipment in the regional electric energy replacement and enhanced coupling of electric transaction and carbon transaction are realized.

The invention provides a regional electric energy substitution distributed optimization method based on electric carbon transaction based on the hot spot problems in the background of novel electric power system construction, carbon transaction market establishment and electric power market innovation, wherein the electric power generation side continuously improves the market competitiveness of clean energy through dynamic price adjustment according to the formed electric-carbon transaction price according to a distribution network based on the marginal electricity price of a distribution node and a regional electric energy substitution double-layer interactive frame, and promotes the use of clean energy; the energy utilization side actively bears the carbon emission cost according to the energy utilization equipment, so that the price advantage of clean electric energy to fossil energy is formed, and the energy utilization side is stimulated to reduce the use of fossil fuel. According to the proposed method for solving the upper model to obtain the marginal electricity price of the distribution node and solving the upper model to obtain the upper and lower iteration of the regional electric energy substitution transaction amount, support can be provided for dynamic adjustment operation of regional electric energy substitution internal energy equipment in different time periods, and flexible interaction, double-layer distributed optimization of electric carbon transaction and regional electric energy substitution operation among the regional electric energy substitution different energy equipment and optimal control of the regional electric energy substitution internal energy equipment are realized. The carbon emission of regional electric energy substitution is reduced while the good economical efficiency of regional electric energy substitution is ensured. The related achievement can be used in trade production work of regional electric energy substitution, and can effectively guide popularization, construction and planning operation of regional electric energy substitution.

According to the distribution network and regional electric energy replacement two-layer interactive framework based on the marginal electricity price of the distribution node, the distribution network can provide dynamic price signals for regional electric energy replacement to guide the regional electric energy replacement to carry out reasonable optimal scheduling; according to the carbon emission characteristics of different energy devices in regional electric energy substitution, what mode is adopted for regional electric energy substitution to set regional electric energy substitution carbon emission limit, and how to determine reasonable carbon trade price of regional electric energy substitution has a critical influence on regional electric energy substitution trade;

the regional electric energy substitution needs to fully show the characteristics of electric carbon transaction in the distributed optimization process, however, as the regional electric energy substitution is in conflict with objective functions of different electric energy substitution in a lower model, excessive coupling constraint exists, and the solving difficulty is high. Therefore, the coupling constraint in the model is decomposed, and the solution of the model is realized on the basis of guaranteeing the data privacy by using a distributed optimization method. The optimization method needs to combine the formed electricity trade price, the electricity market and the carbon trade market trade rules to determine whether the electricity price of the electricity trade can be updated when the regional electric energy substitution is subjected to distributed optimization, so that a new control mode is generated by the regional electric energy substitution.

Example 2:

taking an area electric energy substitution distributed optimization method based on electric carbon transaction as an example in the implementation, referring to fig. 3, the area electric energy substitution distributed optimization method based on electric carbon transaction provided by the invention is described;

an upper model of the regional electric energy to replace the electric price adopted for carrying out electric transaction is established, and the regional electric energy is calculated by adopting a method of marginal electric price of a distribution node to replace the electric price adopted for carrying out electric transaction. Taking an improved IEEE14 node system as an example, the node 1 is connected with a power transmission network, the power purchase price of the power transmission network to the upper-level power transmission network is the power price of the power transmission network at the node, and each power substitution is respectively connected with the node 7, the node 10 and the node 12 to form regional power substitution. 7. And the marginal electricity prices of the distribution nodes of the 10 nodes and the 12 nodes are the electricity prices of each electric energy substitution electricity transaction in the area.

Each electric energy substitution combination forms a regional electric energy substitution lower model containing a gas turbine, a photovoltaic power generation, a storage battery and a load. The curve of the daily load curve, the photovoltaic predicted force curve and the marginal electricity price of the power distribution node of the power transmission network is shown in fig. 4, the data of fig. 4 are imported into a lower model, a Yalmip is used for calling a Gurobi solver to solve, and simulation is carried out through a matlab software environment, so that regional electric energy can be obtained to replace the output.

Solving the lower model to form preliminary scheduling of each internal device, and the specific method is as follows: the calculated marginal electricity prices of the distribution nodes of the power distribution network and the area electric energy substitution connection nodes are used as electricity purchase prices of all electric energy substitution in the area from the upper distribution network, the total cost of the electric energy substitution operation of the lower-layer model area is used as an objective function, and the distributed solving is carried out by adopting a multiplication sub-distributed optimization algorithm in an alternating direction, so that only the electric quantity and the price information are interacted and traded among all the electric energy substitution in the area. The decision variables comprise the output power of each energy power generation device in the regional power substitution, the trade quantity between the regional power substitution and the trade electric quantity of the upper-level power distribution network and the trade quantity between each power substitution in the regional. And calling a Gurobi solver by using a Yalmip toolbox to solve the established regional power substitution model to form preliminary optimized scheduling of regional power substitution devices.

In order to better embody the characteristics of regional electric energy replacement electric carbon transaction, the electric quantity of regional electric energy replacement electric transaction is transmitted to an upper model as a parameter, as shown in fig. 5, the upper model recalculates to generate new power distribution node marginal electricity price, the regional electric energy replacement is solved according to the new obtained electricity price to obtain new control output, and the control output is circulated in sequence, if the power distribution node marginal electricity price obtained by two continuous iterations is smaller than an allowable error value, iteration convergence is considered, and the optimal solution meeting the double-layer model is obtained.

Under the condition of permitting light rejection in one control day, the active electricity price of the distribution node after the iteration convergence of the regional electricity substitution trading volume and the marginal electricity price optimization of the distribution node is considered is shown in figure 6; the power replacement (1) day-ahead optimization result after the consideration of the regional power replacement trading volume and the distribution node marginal power price optimization iteration convergence is shown in fig. 7; the power replacement (2) day-ahead optimization result after the consideration of the regional power replacement trading volume and the distribution node marginal power price optimization iteration convergence is shown in fig. 8; the power replacement (3) day-ahead optimization result after the consideration of the regional power replacement trading volume and the distribution node marginal power price optimization iteration convergence is shown in fig. 9; the power replacement (1) day-ahead optimization result after the iteration convergence of the regional power replacement transaction amount and the power distribution node marginal power price optimization is not considered is shown in fig. 10;

furthermore, in order to verify the effectiveness of the model proposed by the present invention in terms of carbon reduction and photovoltaic digestion, the calculation results of the following three cases were compared based on the above-described examples.

7-9, the sum of the generated energy of each internal electric energy replacement energy device in regional electric energy replacement and the interactive power of the power distribution network is positive, the electricity purchasing is negative, and the sum of the generated energy and the sold electricity is equal to the load value, so that the supply and demand balance is met. At this time, the total carbon emission of electric energy substitute 1 is 72.18t, the total carbon emission of electric energy substitute 2 is 56.05t, and the total carbon emission of electric energy substitute 3 is 104.62t. Regional power replacement purchased from the upper grid is significantly reduced due to the high carbon strength of the gas turbine compared to the situation where carbon trade is not a concern. The total capacity of the photovoltaic for regional power substitution was 14.85MW, the total capacity of the photovoltaic without regional power substitution transaction and DLMP iteration was 10MW, under conditions that allowed for disposal. This illustrates that regional power substitution enhances environmental friendliness while enhancing the electric carbon transaction.

Taking a 12:00 example analysis, the DLMP of FIG. 10 is lower than that of FIG. 7. Through the power interaction between the distribution network and the regional power replacement, the DLMP of the distribution network is higher, the power replacement 1 is prompted to reduce the power purchased from the distribution network, and the P2P interaction power is increased. The total carbon emission of the electric energy substitution 1 on one control day is 79.18t, and compared with fig. 7, the total carbon emission of the electric energy substitution 1 after the electric trade and the carbon trade are interacted is 72.18t, which verifies that the system further reduces the carbon emission while the electric-carbon trade market is further interacted.

After calculation, the corresponding control measures can be executed. There are a variety of control means that can be employed, including: adjusting the output of the gas turbine; regulating grid-connected power of photovoltaic power generation; the input power of the coupling equipment of the thermal power generating unit is regulated, so that the power conversion among different forms of energy sources is realized; and (3) regulating the charge and discharge behaviors of the storage battery to realize load translation. The coupling constraint in the lower model of the regional electric energy substitution is decomposed by adopting an alternate direction multiplication sub-distributed optimization algorithm when the regional electric energy substitution distributed transaction is calculated, in the algorithm, the model solving can be realized by only acquiring a small amount of non-transaction data by the platform, and the data privacy requirement between different electric energy substitution in the region is ensured. The invention provides the association interaction between the marginal electricity price of the distribution node and the upper and lower iteration of the distributed electricity transaction electric quantity of the regional electric energy substitution, and the regional electric energy substitution electricity-carbon transaction interaction is realized, so that the invention plays an important role in reducing the emission of the regional electric energy substitution carbon and improving the absorption of the photovoltaic device as compared with the situation.

In conclusion, the model provided by the invention has certain reference significance for improving the low carbon property and economy of regional electric energy substitution.

Example 3:

an regional power replacement distributed optimization system based on an electrical carbon transaction, comprising:

The system further comprises an upper model building module, wherein the upper model building module is specifically used for:

The system further comprises a lower model building module, wherein the lower model building module is specifically used for:

The lower layer model further comprises a carbon transaction cost calculation model construction sub-module, wherein the carbon transaction cost calculation model construction sub-module is specifically used for:

The calculation formula of the marginal electricity price of the distribution node in the upper model calculation module is as follows:

The lower model calculation module is specifically configured to:

The iteration solving module is specifically configured to:

Example 4:

in a further embodiment of the present invention, a computer device is provided, which includes a processor and a memory, where the memory is configured to store a computer program, the computer program includes program instructions, and the processor is configured to execute the program instructions stored in the computer storage medium. The processor may be a central processing unit (Central Processing Unit, CPU), but may also be other general purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), off-the-shelf Programmable gate arrays (FPGAs) or other Programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc., which are the computational core and control core of the terminal adapted to implement one or more instructions, in particular adapted to load and execute one or more instructions within a computer storage medium to implement the corresponding method flow or corresponding functions; the processor according to the embodiment of the invention can be used for executing the steps of a regional power replacement distributed optimization method based on electric carbon transaction.

Example 5:

in still another embodiment of the present invention, based on the same inventive concept, the present invention further provides a storage medium, in particular, a computer readable storage medium (Memory), which is a Memory device in a computer device, for storing programs and data. It is understood that the computer readable storage medium herein may include both built-in storage media in a computer device and extended storage media supported by the computer device. The computer-readable storage medium provides a storage space storing an operating system of the terminal. Also stored in the memory space are one or more instructions, which may be one or more computer programs (including program code), adapted to be loaded and executed by the processor. The computer readable storage medium herein may be a high-speed RAM memory or a non-volatile memory (non-volatile memory), such as at least one magnetic disk memory. One or more instructions stored in a computer-readable storage medium may be loaded and executed by a processor to implement the corresponding steps of an electric-carbon transaction-based regional power replacement distributed optimization method in the above-described embodiments.

It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.

The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

The foregoing is illustrative of the present invention and is not to be construed as limiting thereof, but rather as providing for the use of additional embodiments within the spirit and scope of the present invention.

Claims

1. An regional power substitution distributed optimization method based on electric carbon transaction, which is characterized by comprising the following steps:

2. The method of claim 1, further comprising a build process of an upper layer model, the build process comprising:

3. The method of claim 1, further comprising a build process of the underlying model, the build process comprising:

4. The method of claim 3, wherein the underlying model further comprises a carbon trade cost calculation model, and wherein the process of constructing the carbon trade cost calculation model comprises:

5. The method of claim 1, wherein the power distribution node marginal electricity price is calculated as follows:

6. The method according to claim 1, wherein the distributively solving the electric transaction and the carbon transaction between the electric energy replacement of each energy device in the area by using the alternate direction multiplier sub-distributively optimizing algorithm based on the distribution node marginal electricity price and the pre-built lower model to obtain the electric quantity of the electric energy replacement electric transaction of the area and the output of each energy device comprises:

7. The method according to claim 1, wherein the step of re-inputting the electric quantity of the regional electric energy replacement electric transaction into the upper model for iterative calculation until obtaining the optimal electric quantity of the regional electric energy replacement electric transaction and the optimal output of each energy device comprises the steps of:

8. An regional power replacement distributed optimization system based on an electrical carbon transaction, comprising:

9. The system of claim 8, further comprising an upper model building module, the upper model building module being specifically configured to:

10. The system of claim 8, further comprising a lower model building module, the lower model building module being specifically configured to:

11. The system of claim 10, wherein the underlying model further comprises a carbon trade cost calculation model building sub-module, the carbon trade cost calculation model building sub-module being specifically configured to:

12. The system of claim 8, wherein the calculation formula of the marginal electricity price of the distribution node in the upper model calculation module is as follows:

13. The system according to claim 8, wherein the lower model calculation module is specifically configured to:

14. The system according to claim 8, wherein the iterative solution module is specifically configured to:

15. A computer device, comprising:

one or more processors;

a processor for executing one or more programs;

an regional power replacement distributed optimization method based on an electrical carbon transaction as claimed in any one of claims 1 to 7, when said one or more programs are executed by said one or more processors.

16. A computer readable storage medium, characterized in that a computer program is stored thereon, which computer program, when executed, implements a regional power replacement distributed optimization method based on an electric carbon transaction as claimed in any one of claims 1 to 7.