CN-121998372-A - Low-carbon scheduling-oriented layered distributed robust operation method for power distribution network-comprehensive energy system

CN121998372ACN 121998372 ACN121998372 ACN 121998372ACN-121998372-A

Abstract

The invention discloses a layered distributed robust operation method of a distribution network-integrated energy system for low-carbon scheduling, which comprises the steps of 1, establishing a multi-layer low-carbon cooperative operation framework of an IES-distribution network based on a carbon coupling electricity price mechanism for realizing low-carbon interaction of different operation bodies of the distribution network and the IES, 2, establishing a low-carbon operation model of the distribution network containing a reconfigurable intelligent soft switch R-SOP for effectively supporting layered optimization decisions on a physical execution level of the distribution network, 3, propelling a system carbon reduction target from a multi-energy flow coupling and cooperation angle, establishing an IES operation model containing various energy devices, and 4, taking the uncertainty of new energy output into consideration, and establishing an IES-distribution network layered distributed robust optimization strategy. The invention ensures interaction efficiency by constructing a multi-layer low-carbon cooperative operation framework of the power distribution network-IES and utilizing the feeder-port dynamic recombination capability of the R-SOP, and also aims at solving privacy problems of different subjects and randomness of new energy output by constructing an ATC-DRO algorithm.

Inventors

YANG XIAODONG
ZHAO XUELI
ZHANG CHENGJIA
YANG YUE
YANG ZHIQING
ZHAO SHUANG
LI HELONG
DING LIJIAN
LOU WEI

Assignees

合肥工业大学

Dates

Publication Date: 20260508
Application Date: 20260211

Claims (10)

1. The hierarchical distributed robust operation method of the power distribution network-comprehensive energy system for low-carbon scheduling is characterized by comprising the following steps of: Step 1, constructing a carbon coupling electric energy interaction mechanism based on carbon potential; Step 2, constructing a recombinant intelligent soft switch R-SOP model; step 3, constructing an operation objective function of the power distribution network ; Step 4, constructing a power distribution network operation model containing R-SOP; step 5, constructing an IES operation model containing various energy devices; step 6, constructing a decoupled power distribution network and IES operation model; Step 7, constructing a distribution network-IES layered distribution robust operation model based on the decoupled distribution network and IES operation model, and converting the distribution network-IES layered distribution robust operation model into a double-layer main model and a double-layer sub-model; and 8, solving the double-layer main model and the double-layer sub model to obtain a final operation scheme comprising R-SOP actions of the reconfigurable soft switch in the power distribution network, prediction actions and adjustment actions of each device in the IES and interaction actions of the power distribution network and the IES.
2. The hierarchical distributed robust operation method for the low-carbon scheduling-oriented power distribution network-integrated energy system according to claim 1, wherein in the step 1, a carbon-coupled electric energy interaction mechanism is obtained by a formula (1) -a formula (3): (1) (2) (3) In the formulas (1) - (3), For nodes in distribution networks Carbon potential of (2); For inflow nodes Is a set of branches; For inflow nodes Is the first of (2) Active power of the branch circuit; For inflow nodes Is the first of (2) Carbon potential of the branch; To the node of the generator G The active power delivered; Is a node The carbon emission intensity of the generator G is connected; Is that Carbon potential matrix of time period distribution network; Is that Node active flux matrix of time period distribution network; Is that A line power flow distribution matrix of the time period distribution network; Is that The output matrix of the generator set of the time period distribution network; Is that Carbon emission factor matrix of the time period generator set; Access node for IES of any integrated energy system Is a carbon coupling unit interaction characteristic quantity; outputting power for a distribution network a unit interaction feature quantity; Receiving power for a distribution network a unit interaction feature quantity; and T represents the transpose.
3. The hierarchical distributed robust operation method for a power distribution network-integrated energy system for low-carbon scheduling according to claim 2, wherein in the step 2, a formula (4) -formula (11) reconfigurable intelligent soft switch R-SOP model is utilized: (4) (5) (6) (7) (8) (9) (10) (11) in the formulae (4) - (11), 、 Respectively representing a node set communicated with the converters and a set of all converters; the number of the nodes is the number of the nodes communicated with the converter; 、、 nodes respectively accessing R-SOP All converters connected are in Total loss of period, dc side power, and ac side power; for accessing R-SOP nodes All converters connected are in Reactive power of the time period; for accessing R-SOP nodes All converters connected are in Apparent power of the period; Is an inverter Is a capacity of (2); is the loss factor of R-SOP; Is indicative of an inverter At the position of Node for communicating with voltage source converter at moment Are connected; for the R-SOP to the node communicating with the voltage source converter Carbon potential when delivering electric energy; for the R-SOP to the node communicating with the voltage source converter The active power delivered; For connecting with voltage source converter Carbon emission intensity of the R-SOP is accessed; From a node in communication with a voltage source converter for R-SOP Carbon potential upon receiving electrical energy.
4. The hierarchical distributed robust operation method for power distribution network-integrated energy system for low-carbon dispatch according to claim 3, wherein in step 3, the operation objective function of the power distribution network is obtained by the formula (12) -formula (15) : (12) (13) (14) (15) In the formulae (12) to (15), Is the first Probability values for the individual discrete scenes; To the power distribution network at the first Running objective functions in discrete scenes; a cost function for receiving power from a superior power grid for the power distribution network; the method is a power distribution network line loss function; is the loss function of R-SOP; interaction functions for the distribution network and all IES; Participating in a carbon market cost function for the power distribution network; the total number of time periods running under a fixed period; output electric energy for upper power grid Unit interactive feature quantity of time period; to be in line with the upper power grid Interactive power of the time period; The node set is a power distribution network node set; The total number of all IES accessed for the distribution network; Penalty coefficients for line loss; Is that Time period node And node Branches between A current flowing therethrough; Is that Time period node And node Branches between Resistance value of (2); Is that Time period and m-th IES' unit power interaction feature; Is that The interactive power of the time interval distribution network and the mth IES; Is a unit carbon emission characteristic quantity; carbon emission generated for the load side of the distribution network; Carbon quota for the distribution network; the total carbon quota for the power distribution network area; Is that Energy for load of the power distribution network in a period of time; Is that The interaction power of the time interval distribution network and the mth IES.
5. The hierarchical distributed robust operation method for a low-carbon scheduling-oriented power distribution network-integrated energy system according to claim 4, wherein in the step 4, an R-SOP-containing power distribution network operation model is constructed by a formula (16) -a formula (22): (16) (17) (18) (19) (20) (21) (22) in the formulae (16) to (22), 、 Is that Time period node And node Branches between Active and reactive power of (a); 、 Is that Time period node And node Branches between Active and reactive power of (a); Is that Time period node And node Branches between Reactance value of (2); 、 Is that Time period node Net injected active and reactive power of (2); 、 Is that Time period node Node Voltage magnitude of (2); Is that Time period distribution network node And interaction power of the mth IES; Is that Time period distribution network node Is set to the load power of (1); 、 The upper and lower limit values of the voltage.
6. The hierarchical distributed robust operation method for a low-carbon dispatch-oriented power distribution network-integrated energy system of claim 5, wherein said step 5 comprises: Step 5-1, constructing the operation objective function of IES using the formulas (24) - (27) : (24) (25) In the formula (24) -formula (25): A day-ahead running objective function for IES; The first of IES Adjusting an objective function under a plurality of discrete scenes; A purchase energy cost function that is IES; device operation and maintenance functions of IES; A wind-discarding and light-punishing function for IES; Is the first An IES's energy purchasing cost adjustment function in a discrete scenario; Is the first An equipment operation cost adjustment function of the IES under discrete scenes; Is the first An IES wind-discarding and light-punishing adjustment function under discrete scenes; Is the first Carbon trade cost function of IES in discrete scenarios; 、 the interactive characteristic quantity is the unit interactive characteristic quantity of electric energy and natural gas; Is that The power distribution network interacts with the mth IES in time period; Is that Natural gas amount received from the upper gas network in the mth IES of the period; Is that Period m natural gas production of P2G devices in IES and D represents the set of all devices in IES, Is a device A unit operation and maintenance feature quantity of (1); Is that Time period device Output power of (a); the set of new energy units comprises a photovoltaic and a fan; is a new energy unit in the mth IES Maximum active power output; Is that New energy unit in mth IES of time period The output actual active power; Is the first In a plurality of discrete scenes The interaction power adjusted by the mth IES of the period and the power distribution network; Is the first In a plurality of discrete scenes The adjusted natural gas amount received in the mth IES of the time period; Is the first In a plurality of discrete scenes The electricity-gas conversion equipment P2G in the m-th IES of the period adjusts the synthesized natural gas quantity; Is the first In a plurality of discrete scenes Time period device The adjusted output power; Is the first In a plurality of discrete scenes New energy unit in mth IES of time period The output actual regulated active power; Is the first Actual carbon emissions of the mth IES in the discrete scenario; Is the first Carbon emission quota of mth IES in the discrete scenes, K represents the total number of discrete scenes, and has: (26) (27) In the formulas (26) - (27): Is that Node carbon potential of the mth IES of the period; 、 Is that Carbon emission of the cogeneration unit CHP and the gas boiler GB of the m-th IES in the period; Is that Carbon capture amount of the mth IES of the period; And 5-2, constructing a constraint model of equipment in the IES, wherein the constraint model comprises a carbon capture system CCS, a constraint model of electricity-to-gas equipment P2G, a constraint model of a cogeneration unit CHP, a constraint model of a gas boiler GB and a constraint model of an energy storage device, the CHP comprises a gas turbine GT and a waste heat boiler WHB, and the energy storage equipment comprises a storage battery BT and a heat storage tank HST.
7. The hierarchical distributed robust operation method for a low-carbon dispatch-oriented power distribution network-integrated energy system of claim 6, wherein said step 5-2 comprises: Step 5-2-1, construction using formula (28) and formula (31) Constraint model of the period CCS: (28) (29) (30) (31) formula (28) -formula (31): carbon emission intensity for CCS; Is the heat value of natural gas; Is that Natural gas amount consumed by CCS in the mth IES of period; Is that Fixed energy consumption of CCS in the mth IES of the period; The energy consumption for the CCS operation in the mth IES; The energy consumption required for carbon dioxide treatment units of CCS; Step 5-2-2, construction using formula (31) and formula (32) Operational model of period P2G: (32) (33) in the formula (31) -formula (32): the heat energy can be converted by unit electric energy; The working efficiency of P2G; Is that Power consumed by P2G in the mth IES of the period; Is that Carbon dioxide amount consumed by P2G in the mth IES of the period; the amount of carbon dioxide required to generate a unit power of natural gas; step 5-2-3, construction using formula (31) and formula (32) Constraint model of period CHP unit: (34) (35) formula (34) -formula (35): 、、 3 carbon emission coefficients for the gas plant; Is that The generated power of the period CHP; Is that Heating power of the period CHP; Is the power generation efficiency of CHP; Is that Natural gas consumed by CHP during periods; Is that Exhaust residual heat of the period CHP; the heat dissipation loss rate is the heat dissipation loss rate; The waste heat recovery efficiency is achieved; the flue gas recovery rate of CHP; Step 5-2-4, construction using formula (36) and formula (37) A constraint model of the time period gas boiler GB; (36) (37) In the formula (36) -formula (37): Is that Heating power of the period GB; heating efficiency is GB; Is that Natural gas consumed in period GB; step 5-2-5, construction using formula (38) and formula (39) Constraint model of time period energy storage device: (38) (39) Formula (38) -formula (39): Is that The state of charge of the period BT, Representation of The state of charge of the period BT; 、 Is that Charge and discharge power of the period BT; 、 Is that Charging and discharging efficiency of the period BT; Is that The heat storage energy of the time period HST; Is that The heat storage energy of the time period HST; 、 Is that The heat storage and release power of the time period HST; Is the energy self-loss rate of HST; 、 Is the heat storage and release efficiency of HST.
8. The hierarchical distributed robust operation method for a low carbon dispatch-oriented power distribution network-integrated energy system of claim 7, wherein said step 6 comprises: Step 6-1, construction of the first Using formula (40) Penalty function of decoupled IES operation model under secondary internal loop : (40) In formula (40): 、 Is that Lagrangian penalty function first term and second term multipliers corresponding to the mth IES of the period; Is the first Under the secondary internal circulation The mth IES of the period expects interactive power with the distribution network; Is the first -1 Internal circulation The period distribution network expects interactive power with the mth IES; step 6-2, construction of the first Using formula (41) Punishment function of power distribution network operation model after decoupling under secondary internal circulation : (41) In the formula (41), Is the first -1 Internal circulation The mth IES of the period expects interactive power with the distribution network; Is the first Under the secondary internal circulation The period distribution network expects interactive power with the mth IES; step 6-3, construction of the first Using formula (42) Decoupling constraints under the secondary inner loop: (42) In formula (42): Is that The interaction power value of the mth IES in the period and the power distribution network; is the maximum value of the interaction power of the IES and the power distribution network.
9. The hierarchical distributed robust operation method for a low carbon dispatch-oriented power distribution network-integrated energy system of claim 8, wherein said step 8 comprises the steps of: step 7-1, constructing a double-layer main model by using the formula (45) -formula (46): (45) (46) Formula (45) -formula (46): an IES objective function that is a bilayer master; the power distribution network target function is a double-layer main model; An IES constraint goal for a bilayer submodel; A constraint target of the power distribution network of the double-layer sub-model is adopted; coefficients that are objective functions; Predicting a force value for the new energy; Is the first New energy output values of the discrete scenes; Is the first Probability values for the individual discrete scenes; Step 7-2, constructing a double-layer submodel by using the formula (47) -formula (48): (47) (48) in the formulae (47) - (48), An IES objective function that is a bilayer submodel; The power distribution network objective function is a double-layer sub-model; a feasible domain which is a probability constraint, and comprises: (49) in formula (49): The initial probability value of the new energy output scene is set; 、 the bias values are allowed for probabilities under the 1-norm and + -norm constraints, respectively.
10. The hierarchical distributed robust operation method for a low carbon dispatch-oriented power distribution network-integrated energy system of claim 9, wherein said step 8 comprises the steps of: step 8-1, defining CCG outer layer circulation times as Defining the circulating times of the ATC inner layer as Initializing =1; Initialize the first Lower distribution network target upper bound of secondary outer loop First, the Lower boundary of power distribution network target under secondary external circulation First, the IES target upper bound under minor outer loop First, the Lower IES target bound under secondary outer circulation =0 And the first Under the secondary external circulation Probability values for discrete scenes Definition of The external circulation convergence precision is obtained; step 8-2 initializing =1, Initialize the first Second outer circulation Lagrange penalty function once term corresponding to mth IES in t period under secondary internal circulation Initializing the first to a fixed value between (0,0.1) Second outer circulation Lagrange penalty function quadratic term corresponding to mth IES of t period under secondary internal circulation Initializing the first to a fixed value between (0,0.1) Second outer circulation IES objective function of double-layer main model under secondary internal circulation =0, Initialize the first Second outer circulation Power distribution network objective function of double-layer main model under secondary internal circulation =0; Step 8-3 calculating the first according to the formulas (45) to (46) The second external circulation IES total objective function of double-layer main model under secondary internal circulation Day-ahead objective function of IES And the first thereof Second outer circulation Under the secondary internal circulation Period mth IES expects interaction power with power distribution network General objective function of power distribution network of double-layer main model And the first thereof Second outer circulation Under the secondary internal circulation Period distribution network expects interactive power with mth IES Thereby obtaining the first Second outer circulation IES action scheme set under secondary internal circulation Comprises the following steps of Second outer circulation Each device in the IES predicts actions and the power distribution network and the IES interact actions under the secondary internal circulation; Step 8-4, judging whether the formula (50) -formula (51) is satisfied, if so, making the first step Lower boundary of power distribution network target under secondary external circulation First, the Lower IES target bound under secondary outer circulation After that, step 8-5 is executed, otherwise Assignment to After updating the penalty function according to the formula (52), returning to the step 8-3 for sequential execution; (50) (51) (52) in the formulae (52) - (54), The accuracy value of the convergence condition 1 of the ATC algorithm is obtained; the accuracy value of the convergence condition 2 of the ATC algorithm is obtained; Is the first Second outer circulation A Lagrangian penalty function primary term corresponding to the mth IES in the t period under the secondary internal circulation; Is the first Second outer circulation A Lagrange penalty function quadratic term corresponding to the mth IES in the t period under the secondary internal circulation; Is the first The second external circulation Under the secondary internal circulation The mth IES of the period expects interactive power with the distribution network; Is the first The second external circulation Under the secondary internal circulation The period distribution network expects interactive power with the mth IES; step 8-5 initializing =1, Initialize the first Second outer circulation Lagrange penalty function once term corresponding to mth IES in t period under secondary internal circulation Initializing the first to a fixed value between (0,0.1) Second outer circulation Lagrange penalty function quadratic term corresponding to mth IES of t period under secondary internal circulation Initializing the first to a fixed value between (0,0.1) Second outer circulation IES objective function of double-layer submodel under secondary internal circulation =0, Initialize the first Second outer circulation Power distribution network objective function of double-layer sub-model under secondary internal circulation =0; Step 8-6 based on Calculate the first according to the formulas (53) and (54) Under the secondary external circulation Probability values for discrete scenes First, the The second external circulation IES total objective function of bilayer submodel under secondary internal circulation And the first thereof The second external circulation Under the secondary internal circulation Period mth IES expects interaction power with power distribution network General objective function of power distribution network of double-layer main model And the first thereof Second outer circulation Under the secondary internal circulation Period distribution network expects interactive power with mth IES Thereby obtaining the first Second outer circulation IES adjustment action scheme set under secondary internal circulation First of all The second external circulation Action scheme set of power distribution network under secondary internal circulation , wherein, Comprises the first step of Second outer circulation Individual device tuning actions in IES and power distribution network and IES interaction tuning actions under the secondary internal loop, Includes the first Second outer circulation R-SOP action under secondary internal circulation; (53) (54) in the formulas (53) and (54), Is the first Second outer circulation Optimal objective function of IES in double-layer submodel under secondary internal circulation under kth discrete scene; Is the first Second outer circulation An optimal objective function of the power distribution network in the double-layer sub-model under the secondary internal circulation under the kth discrete scene; the new energy output set is the new energy output set in the kth discrete scene; 8-7, judging whether the formula (55) -formula (56) is satisfied, and if so, obtaining the first step Lower distribution network target upper bound of secondary outer loop First, the IES target upper bound under minor outer loop After that, step 8-8 is executed, otherwise Assignment to Then, updating the penalty function according to the formula (52), and returning to the step 8-6 for sequential execution; (55) (56) Step 8-8, judging And Whether or not it is true, if so, then the final operation scheme is obtained Otherwise, will Assignment to And then, returning to the step 8-2 for sequential execution.

Description

Low-carbon scheduling-oriented layered distributed robust operation method for power distribution network-comprehensive energy system Technical Field The invention belongs to the field of operation optimization of flexible power distribution networks containing comprehensive energy systems, and particularly relates to a hierarchical distributed robust operation method of a power distribution network-comprehensive energy system oriented to low-carbon scheduling. Background In recent years, integrated energy systems (INTEGRATED ENERGY SYSTEM, IES) have been widely deployed in power distribution networks (distribution network, DN) as autonomous bodies capable of coordinating heterogeneous source-load resources of multiple classes. By means of source-network-charge-storage coordination control and a multi-energy complementary mechanism, the IES actively promotes new energy grid connection and integrates various energy forms such as electricity, heat, gas and the like while realizing economical, efficient and clean operation, and becomes an important direction of current energy system research. The overall flexibility and economy can be effectively balanced by flexibly planning distributed resources such as intelligent soft point (SOP) and energy storage. And the recombined intelligent soft switch (reconfigurable soft open point, R-SOP) can switch feeder lines, adjust port capacity and improve active power transmission efficiency compared with multi-port SOP with the same size capacity. It can be seen that R-SOP is more advantageous than conventional SOP in helping the interactive regulation of the distribution network. However, the traditional R-SOP relies on a mechanical switch to dynamically switch port capacity, is limited by action speed and response precision, and is difficult to meet the requirement of rapid power flow regulation in a complex power distribution scene. Therefore, improvement of R-SOP is needed to improve the action speed and control accuracy, respond quickly and adjust the power transmission path dynamically, and guide the low-carbon electric energy to flow to the high-demand area effectively, so as to optimize the allocation mode of carbon emission responsibility and further improve the low-carbon benefit and flexibility of the operation of the distribution network. Meanwhile, in order to solve the problem of multi-layer collaborative operation of the power distribution network and the IES, a distributed algorithm is often adopted to realize optimal solution among different main bodies. Distributed optimization algorithms such as an alternate direction multiplier method, a consistency algorithm, a target cascading method and the like all have different advantages in the field, wherein the target cascading method (ANALYTICAL TARGET CASCADING, ATC) has greater flexibility in the aspect of coordination of sub-problems, and various types of penalty functions can be applied to different levels of a power distribution network and therefore various types of distributed solutions. However, when the distribution network and the IES interaction strategy are in the face of the new energy uncertainty scene, the phenomena of overlarge scheduling deviation and the like occur. Therefore, the existing deterministic analysis strategy is difficult to be applied to the interactive analysis of the IES-distribution network considering the uncertainty of new energy. Disclosure of Invention The invention provides a layered and distributed robust operation method of a power distribution network-integrated energy system for low-carbon scheduling, aiming at solving the defects existing in the prior art. The method has the advantages of improving flexibility and stability of power flow regulation and control of the power distribution network through optimizing a system operation mode, improving interaction efficiency of the power distribution network and the IES, standardizing carbon emission responsibility allocation standards, promoting deep collaborative carbon reduction of the power distribution network and the comprehensive energy system, converting the power-assisted energy system into low carbonization and clean energy, guaranteeing economy of system operation and reliability of coping with random characteristics of new energy output, providing technical support for low-carbon scheduling of the power distribution network-IES in the background of a novel power system, conforming to energy low-carbon development strategy, and having important engineering application value and environmental protection significance. In order to achieve the aim of the invention, the invention adopts the following technical scheme: The invention relates to a layered distributed robust operation method of a distribution network-integrated energy system oriented to low-carbon scheduling, which is characterized by comprising the following steps: Step 1, constructing a carbon coupling electric energy interaction mechanism based on car