CN105703369B - Optimal energy flow modeling and solving method for multi-energy coupling transmission and distribution network - Google Patents

Optimal energy flow modeling and solving method for multi-energy coupling transmission and distribution network Download PDF

Info

Publication number
CN105703369B
CN105703369B CN201610081125.0A CN201610081125A CN105703369B CN 105703369 B CN105703369 B CN 105703369B CN 201610081125 A CN201610081125 A CN 201610081125A CN 105703369 B CN105703369 B CN 105703369B
Authority
CN
China
Prior art keywords
energy
distribution network
power
network
optimal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610081125.0A
Other languages
Chinese (zh)
Other versions
CN105703369A (en
Inventor
马瑞
李晅
颜宏文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201610081125.0A priority Critical patent/CN105703369B/en
Publication of CN105703369A publication Critical patent/CN105703369A/en
Application granted granted Critical
Publication of CN105703369B publication Critical patent/CN105703369B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to a modeling and solving method for an optimal energy flow of a multi-energy coupling transmission and distribution network, belonging to the research field of a day-ahead scheduling plan of an electric power system under energy interconnection and comprising the following steps: acquiring basic data in a system scheduling period, and acquiring a random fuzzy time-space sequence model of large-scale wind power, distributed power supplies and multi-energy loads through historical data mining; the power and voltage of the power transmission network and each active distribution network at the common node are used as shared variables, and a multi-target system engineering theoretical dynamic optimal energy flow model which gives consideration to economy, low carbon and renewable energy consumption and loss reduction under the static safety constraint is constructed; and (3) obtaining multi-energy source load prediction by random fuzzy simulation, and obtaining a Pareto solution set, an optimal compromise solution and an energy flow result by adopting an improved system engineering theory hierarchical optimization algorithm based on approximate dynamic programming and NSGA-II. The invention adapts to the development trend of energy interconnection and realizes the day-ahead scheduling comprehensive coordination optimization of transmission and distribution parties on the premise of meeting the static safety and stability of the system.

Description

Optimal energy flow modeling and solving method for multi-energy coupling transmission and distribution network
Technical Field
The invention belongs to the field of research on day-ahead scheduling plans of power systems under energy interconnection, and relates to a method for modeling and solving optimal energy flows of a multi-energy coupling transmission and distribution network.
Background
Under the background of energy Internet, the development of distributed renewable energy sources and energy source hubs in active distribution networks is beneficial to realizing energy conservation and emission reduction and comprehensive and efficient utilization of multiple energy sources. The source-load side power injection of the multi-energy system has uncertainty and obvious space-time difference, and the traditional scheduling method prohibits the distribution network from transmitting power to the transmission network, so that part of renewable energy sources are difficult to be locally consumed and reduced, and therefore the feeder line power interaction of the active distribution network and the transmission network is a necessary trend. Based on the method, the power transmission network is fully utilized as a carrier, and the power transmission network and each active distribution network are cooperatively scheduled under safety constraint, so that the space-time complementation of comprehensive energy in a wide area range is promoted, and the method has certain significance for comprehensive utilization of multiple energy sources and consumption of renewable energy sources.
In the past scheduling model, a power transmission network or a distribution network is often used as research objects respectively, the equivalence of the distribution network is used as a PQ load node during analysis of the power transmission network, the equivalence of the power transmission network is used as an infinite power supply during analysis of the distribution network, and the competition of the power transmission network and the distribution network and the internal multi-energy coupling of the active distribution network are not considered. In recent years, researchers pay attention to the integrated scheduling problem research of power transmission and distribution networks, and describe the problem as a decomposition coordination unit combination optimization model and solve the problem through an iterative optimization algorithm based on Benders decomposition, a system engineering theory and the like. Multiple energy sources such as electricity, gas and heat are coupled, converted and supplied through an energy hub, so that corresponding energy hub modeling and scheduling research is developed by students. Also, based on the phenomenon that a natural gas network and a power network are coupled and interconnected through an energy hub, a student researches an analysis method of a multi-energy hybrid energy flow and analyzes network energy interaction information characteristics.
In summary, a multi-energy coupling transmission and distribution network optimal energy flow modeling and solving method under the energy interconnection is researched to comprehensively coordinate the multi-target optimization of all transmission and distribution parties on the premise of meeting the static safety and stability of the system, such as production cost, pollutant gas emission, renewable energy consumption, line loss and the like.
Disclosure of Invention
The technical problem to be solved is as follows: aiming at the defects of the prior art, the invention discloses a multi-energy coupling transmission and distribution network optimal energy flow modeling and solving method, which provides a phenomenon of multiple uncertain injection of multi-energy sources on the source and load sides of a power system under energy interconnection, constructs a multi-target optimal energy flow model of a multi-energy coupling active distribution network and power transmission network integrated system based on an SoS idea, researches a solving step method of the multi-target optimal energy flow model, and aims to provide theoretical reference for a wide area energy source space-time complementary absorption scheduling method under the energy interconnection background.
The technical scheme is as follows: a modeling and solving method for optimal energy flow of a multi-energy coupling transmission and distribution network comprises the following steps:
step 1: basic data of the system in a scheduling period are obtained, and a random fuzzy time-space sequence model of large-scale wind power output, distributed power supply output and multi-energy load is obtained through historical data mining.
Step 2: and taking the power and voltage of the power transmission network and each active distribution network in the common node set as shared variables, and constructing a multi-objective system engineering theoretical dynamic optimal energy flow model which gives consideration to economy, low carbon and renewable energy consumption and loss reduction under the constraint of static safety.
And step 3: the method comprises the steps of obtaining a predicted value of a multi-energy source load space-time sequence through random fuzzy simulation, judging different operation modes of an active distribution network, and obtaining a Pareto optimal solution set, an optimal compromise solution and a corresponding energy flow result by adopting an improved system engineering theory layered optimization algorithm based on approximate dynamic programming and NSGA-II.
Has the advantages that: the invention is suitable for the development trend of multi-energy access, and effectively realizes the daily scheduling comprehensive coordination optimization of transmission and distribution parties on the premise of meeting the static safety and stability of the system.
Drawings
FIG. 1: the invention relates to an implementation process of an optimal energy flow modeling and solving method of a multi-energy coupling transmission and distribution network;
FIG. 2: schematic diagram of an example energy system of the invention;
FIG. 3: the invention relates to a flow chart of a solving algorithm of a multi-target random fuzzy dynamic optimal energy flow.
Detailed Description
The present invention will be described in detail with reference to the accompanying drawings and specific embodiments.
The overall implementation flow of the optimal energy flow modeling and solving method for the multi-energy coupling transmission and distribution network provided by the invention is shown in fig. 1, and a certain energy system is taken as a specific embodiment to be explained in detail, and a schematic diagram of the energy flow modeling and solving method is shown in fig. 2. The examples are intended to illustrate but not to limit the invention.
Step 1: basic data of the system in a scheduling period are obtained, and a random fuzzy time-space sequence model of large-scale wind power output, distributed power output, multi-energy load and the like is obtained through historical data mining.
Acquiring basic information and data of basic topological structures, branch impedances and thermal power unit output limits of a power transmission network, each active distribution network and a natural gas network in a scheduling period; based on the mining and distribution parameter fitting of historical data of wind speed, illumination intensity and electricity/gas/heat load at different time/regions, a random fuzzy space-time sequence model and an opportunity measurement function thereof are obtained to describe a large-scale wind power plant, a distributed power supply, a multi-energy load and the like.
Step 2: and taking the power and voltage of the power transmission network and each active distribution network in the common node set as shared variables, and constructing a multi-objective system engineering theoretical dynamic optimal energy flow model which gives consideration to economy, low carbon and renewable energy consumption and loss reduction under the constraint of static safety.
The main designations, variables, subscripts, etc. appearing in the model are first described as follows:
t power transmission network;
d, actively distributing a network;
an H energy hub;
a period t;
i, numbering the nodes of the power transmission network;
j, actively distributing the number of the internal nodes of the network;
k, numbering natural gas network nodes;
ΦTa set of transmission grid nodes;
ΦT-Gthe nodes of the power transmission network and the thermal power generating unit are gathered,
Figure GDA0002451775160000021
ΦT-tDthe nodes connected with the traditional distribution network in the transmission network are gathered,
Figure GDA0002451775160000022
ΦT-Dthe nodes connected with the active distribution network in the power transmission network are gathered,
Figure GDA0002451775160000031
ΦDactively distributing an internal node set of the network;
ΦD-Hnodes connected with a natural gas network inside the active distribution network, namely an energy hub node set;
xtlocal decision variable vectors of the power transmission network, active power output of thermal power, voltage, wind power absorption power and the like;
ytactive distribution network local decision variable vectors, namely the absorbed power of a DG below the active distribution network local decision variable vectors, the natural gas inflow of an energy hub node, the electric power inflow, operation mode parameters and the like;
ztpower transmission network and active power distributionNet shared variable vector zt={(PLi,t,Ui,t)|i∈ΦT-DI.e. the active power, voltage of the connection node, is a subset of the decision variables.
The invention constructs an optimal energy flow model of a multi-energy coupling transmission and distribution network as follows:
(1) and (3) upper layer: power transmission network optimization model
1) An objective function:
a) economic objective
The economic objective of a power grid can be described as maximizing the benefit of the electricity purchasing/selling link, purchasing behavior including purchasing electricity from thermal power plants and active distribution networks that can reverse power, selling behavior including selling electricity to the active distribution networks and traditional distribution networks as loads, and cost of wind farms with additional consideration.
Figure GDA0002451775160000032
Wherein: a isi、bi、ciGenerating cost parameters of the generator set positioned at the node i; pGi,tThe active output of the generator set positioned at the node i, and the inflow node is positive; pLi,tFor the load power of the distribution network at the node i, the outflow node is positive, that is, the direction from the transmission network to the distribution network is positive, when the active distribution network is used as the power supply, PLi,t< 0, the second term in the above formula and the negative sign are positive after the second term and the negative sign indicate that the transmission network purchases electricity from the active distribution network side ξiAnd rhoiThe electricity selling (purchasing) fees of the active distribution network and the traditional distribution network are respectively represented and are positive numbers.
Figure GDA0002451775160000033
b) Low carbon goal
Described as thermal power generating unit CO2The emissions are minimized.
Figure GDA0002451775160000034
Wherein αi、βi、γiFor generating sets CO at node i2An emission parameter.
c) Renewable energy consumption target
Described as the maximum contribution of a scaled wind farm.
Figure GDA0002451775160000041
Wherein:
Figure GDA0002451775160000042
representing predicted output of the wind farm, PWG,tWhich is the actual force applied.
d) Loss reduction target
The network active loss described as the transmission grid is minimal.
Figure GDA0002451775160000043
Wherein: u shapei,tRepresents the voltage of node i; bi represents a node connected to the node i, and its set is ΓT,Ubi,tRepresents the voltage at node bi; gi,bi,tAnd Bi,bi,tRespectively representing the conductance and susceptance of a connecting branch of a node i and a node bi; deltai,bi,tRepresenting the phase angle difference.
2) Constraint conditions are as follows:
a) power balance constraint
For the joint node (namely the shared variable z) of the transmission network and the active distribution networktNodes involved):
Figure GDA0002451775160000044
for other nodes of the transmission network:
Figure GDA0002451775160000045
the difference between the formulas (5) and (6) is that P in the formula (5)Li,tAnd Ui,tBeing a decision variable of the transmission networkAnd (b) a portion.
b) Upper and lower limit constraints of node voltage
Figure GDA0002451775160000046
c) Branch transmission power constraint
Figure GDA0002451775160000047
d) Power constraint of feeder line for connecting power transmission network and active distribution network
Figure GDA0002451775160000048
e) Thermal power generating unit constraint
The method comprises the active output upper and lower limits, the reactive output upper and lower limits and the climbing constraint of the generator.
Figure GDA0002451775160000049
Figure GDA00024517751600000410
Figure GDA00024517751600000411
f) System spinning reserve
Per system maximum load
Figure GDA00024517751600000412
And a certain percentage upsilon of the system is considered for standby rotation and is borne by the thermal power generating unit.
Figure GDA0002451775160000051
g) Node static safety margin constraint
Figure GDA0002451775160000052
Wherein the content of the first and second substances,
Figure GDA0002451775160000053
is a reactive load of QLi,tVoltage breakdown point of time, Ai、Bi、CiRespectively, the P-Q-V quiescent voltage stability boundary expression coefficients of node i,
Figure GDA0002451775160000054
is the largest load margin and λ is the stability margin.
(2) The lower layer: active distribution network optimization model (with active distribution network at node i of power transmission network as research object)
1) An objective function:
a) economic objective
The economic objective of the active distribution network is described as minimizing the electricity purchase costs, including those from the main network and from its internal distributed energy sources.
Figure GDA0002451775160000055
Wherein, ξDG,jRepresenting the electricity price, P, of the electricity purchased from the distributed energy source in the active distribution networkDG,j,tRepresenting distributed energy contribution.
Also have
Figure GDA0002451775160000056
The active distribution network operation mode and its determination will be discussed in detail in the algorithm description of step 3.
b) Low carbon goal
CO for a gas turbine of a cogeneration unit described as an active distribution network lower energy hub node2The discharge amount is the smallest.
Figure GDA0002451775160000057
Wherein eEHCO of gas turbine of cogeneration unit representing energy hub node2Coefficient of emission, CgeRepresenting the conversion coefficient, P, of natural gas energy into electric powergj,tIndicating natural gas inflow
c) Renewable energy consumption target
Wind/light distributed renewable energy described as an active distribution network gives the greatest absorption.
Figure GDA0002451775160000058
Wherein, PDGj,tAnd
Figure GDA0002451775160000059
the actual internet power and the predicted power of the distributed power supply are respectively.
d) Loss reduction target
The network described as an active distribution network has minimal active loss.
Figure GDA0002451775160000061
2) Constraint conditions are as follows:
a) active distribution network node power balance constraint
For the node (i.e. shared variable z) connected with the transmission network in the active distribution networktNodes involved, i.e. active distribution network bus):
Figure GDA0002451775160000062
Figure GDA0002451775160000063
wherein P isESSj,tThe power of the energy storage device is positive with its charged state as a load. Then formula (17) means: adding the sum of all electric power loads in the active power distribution network to all energy storage device powers and all network losses, and deducting all distributed energy output to obtain bus exchange power PLi,tIts positive and negative blocksAnd determining the operation mode of the active distribution network.
For energy hub node, PLj,tSubstantially representing the electric power flowing from the grid into the energy hub, where the natural gas energy flows into Pgj,tOther energy sources such as distributed energy sources are introduced as
Figure GDA0002451775160000064
The electrical load under the energy hub is
Figure GDA0002451775160000065
Natural gas load of
Figure GDA0002451775160000066
A heat load of
Figure GDA0002451775160000067
Wherein the superscript "" indicates that the variable is a random ambiguity uncertainty. The multi-energy power balance constraint of the energy hub node is expressed as
Figure GDA0002451775160000068
For other nodes of the active distribution network:
Figure GDA0002451775160000069
b) active distribution network and transmission network connection feeder power constraint
Figure GDA00024517751600000610
c) Active distribution network node voltage constraints
Figure GDA00024517751600000611
d) Natural gas network node flow balance constraint
Figure GDA00024517751600000612
Wherein
Figure GDA0002451775160000071
Figure GDA0002451775160000072
The natural gas flow rate of the pipeline between two natural gas nodes k and bk;
Figure GDA0002451775160000073
e) natural gas network node pressure constraints
Figure GDA0002451775160000074
f) Natural gas network pipeline flow restriction
Figure GDA0002451775160000075
g) Energy storage device restraint
The energy storage device at the node j should satisfy the charge and discharge power constraint, the electric quantity state constraint and the like during operation scheduling.
Figure GDA0002451775160000076
Figure GDA0002451775160000077
And step 3: the method comprises the steps of obtaining a predicted value of a multi-energy source load space-time sequence through random fuzzy simulation, judging different operation modes of an active distribution network, and obtaining a Pareto solution set, an optimal compromise solution and a corresponding energy flow result by adopting an improved system engineering theory layered optimization algorithm based on approximate dynamic programming and NSGA-II.
The model created in step 3 is first transformed. In the optimization model of the upper-layer power transmission network,will share variable ztLocal variables η expressed as upper layer optimizationstIn the lower-layer active distribution network optimization model, z is settExpressed as a local variable mutThus ensuring ηt=μtThe optimization results of the transmission network and the active distribution network can be consistent. Defining a relaxation variable ct=ηttAnd its penalty function fπ(ct)=σctWhere σ is a dynamic penalty coefficient that increases gradually as the number of iterations increases, fπ(ct) Attached to each objective function. Thus, the double-layer optimization models are decoupled, optimized and mutually constrained, and iteration is carried out until ctIf the value is less than epsilon, a decomposition coordination optimization solution is obtained.
The algorithm flow is shown in fig. 3 and is set forth below.
1) Generating a time-space sequence of output of a large-scale wind power plant, distributed wind/light and multi-energy load of each node by adopting a random fuzzy simulation method;
2) respectively setting the iteration times of the inner layer and the outer layer to zero, namely setting l to 0 and w to 0, and initializing and setting a penalty coefficient sigma;
3) in order to enable the active distribution network to fully exert self-sufficient coordination capability and realize the consumption of internal distributed renewable energy sources, firstly, the problem of multi-target dynamic optimal energy flow of the lower active distribution network is solved, and the corresponding mu at the moment is obtained according to the optimization resultt(0)={(PLi,t,Ui,t)|i∈ΦT-D},PLi,tThe positive/negative/zero of (1) respectively corresponds to three operation modes of load/power supply/island of the active distribution network, if the operation mode is the load/power supply mode, the operation is switched to 5), and if the operation mode is the island mode, the operation is switched to 4);
4) correspondingly disconnecting the contact feeder of the power transmission network and the active distribution network in the island mode, completely decoupling the power transmission network and the active distribution network for operation, and performing independent optimization solution on the power transmission network and the active distribution network to 9);
5) w is w +1, in min { f }n[xtt(w)]+fπt(w)-μt(w-1)]1,2,3 and 4, solving the multi-target dynamic optimal power flow of the upper layer power transmission network to obtain ηt(w);
6) In min { fn[ytt(w)]+fπt(w)-μt(w)]1,2,3 and 4, solving the multi-target dynamic optimal energy flow of the lower active distribution network, and acquiring mut(w);
7) Determining η whether an inner iteration convergence criterion is satisfiedt(w)-ηt(w-1) < ε and μt(w)-μt(w-1) < epsilon, if not, turning to 4), and if yes, turning to 8);
8) determining η whether an outer iteration convergence criterion is satisfiedt(w)-μt(w) < epsilon, if not, l is l +1, updating the penalty coefficient sigma, and turning to the step 4), and if yes, turning to the step 9);
9) and outputting the optimized solution and the energy flow result, and finishing the algorithm.
In the above process:
(1) in step 3), in order to make the active distribution network fully exert self-sufficient coordination capability and realize the consumption of internal distributed renewable energy, three operation modes in the optimization of the active distribution network and the determination thereof are explained as follows.
After a multi-energy load space-time prediction sequence generated by random fuzzy simulation, when an active distribution network is optimized, firstly calculating a formula (19) in the current state, if P isLi,tIf the current value is less than 0, indicating that surplus is generated in distributed energy resources in the active distribution network, firstly considering to call an energy storage device for charging and storing, and if surplus still exists after the energy storage device is fully called on the premise of meeting safety constraint, transmitting power to the power transmission network on the premise of meeting safety and bus power constraint conditions, wherein the active distribution network is in a power supply state; if it is initially PLi,tIf the self-supply of the active distribution network cannot be realized, the discharge supply of the energy storage device is considered to be called, if the self-supply of the energy storage device cannot be realized after the self-supply is called fully, the electric power is absorbed from the power transmission network, and the active distribution network operates in a load state; if P is enabled in the initial state or by calling up the energy storage deviceLi,tAnd (5) considering the active distribution network as an isolated island and optimizing the decoupling operation with the power transmission network when the active distribution network is 0. In the above process, if complete consumption of renewable energy cannot be achieved under the constraint conditions, it is considered to reduce a partial output by means of active management.
(2) The adopted multi-target dynamic optimal energy flow solving method is a multi-target optimal energy flow algorithm based on approximate dynamic programming and NSGA-II.
The method comprises the steps of obtaining a Pareto optimal solution set by adopting NSGA-II, calculating the satisfaction degree of each target value of each non-dominated solution in the Pareto optimal solution set by adopting a partial small fuzzy satisfaction calculation formula, calculating the comprehensive satisfaction degree of each non-dominated solution, and selecting a path with the maximum accumulated satisfaction degree through strategy iterative approximate dynamic programming to form the solution of the multi-target dynamic optimal energy flow.
The above embodiments are merely illustrative, and not restrictive, and those skilled in the relevant art can make various changes and modifications without departing from the spirit and scope of the invention, and therefore all equivalent technical solutions also fall within the scope of the invention.

Claims (6)

1. A method for modeling and solving optimal energy flow of a multi-energy coupling transmission and distribution network is characterized by comprising the following steps:
step 1: acquiring basic data of the system in a scheduling period, and acquiring a random fuzzy time-space sequence model of large-scale wind power output, distributed power output and multi-energy load through historical data mining;
step 2: taking the power and voltage of the power transmission network and each active distribution network in the common node set as shared variables, and constructing a system engineering theoretical dynamic optimal energy flow model which gives consideration to economic, low-carbon and renewable energy consumption and loss reduction multiple targets under the constraint of static safety;
and step 3: the method comprises the steps of obtaining a predicted value of a multi-energy source load space-time sequence through random fuzzy simulation, judging different operation modes of an active distribution network, and obtaining a Pareto optimal solution set, an optimal compromise solution and a corresponding energy flow result by adopting an improved system engineering theory layered optimization algorithm based on approximate dynamic programming and NSGA-II.
2. The optimal energy flow modeling and solving method for the multi-energy coupling transmission and distribution network according to claim 1, wherein basic information and data of basic topological structures, branch impedances and thermal power generating unit output limits of the transmission network, each active distribution network and the natural gas network in a dispatching cycle are obtained in step 1; and acquiring a random fuzzy space-time sequence model and an opportunity measurement function thereof to describe large-scale wind power, distributed power and multi-energy loads based on mining and distribution parameter fitting of historical data of wind speed, illumination intensity, electricity/gas/heat load at different time/regions.
3. The method for modeling and solving the optimal energy flow of the multi-energy coupling transmission and distribution network according to claim 1, wherein step 2 considers static safety stability margins in static safety constraints, and considers the safety constraints of the stochastic fuzzy multi-energy supply and demand balance and the natural gas network in the active distribution network model, wherein the safety constraints are coupled with the natural gas network at energy hub nodes.
4. The method for modeling and solving the optimal energy flow of the multi-energy coupling transmission and distribution network according to claim 1, wherein in the step 3, a double-layer optimization model is converted to solve the coupling coordination, and a shared variable z is used in an upper-layer transmission network optimization modeltExpressed as its local variable ηtZ is set in the optimization model of the lower active distribution networktExpressed as a local variable mutDefine the relaxation variable ct=ηttAnd its penalty function fπ(ct)=σctWhere σ is a dynamic penalty coefficient that increases gradually as the number of iterations increases, fπ(ct) Adding the parameters to each objective function to enable the double-layer optimization model to be decoupled and optimized and mutually constrained, and iterating until ctIf the value is less than epsilon, a decomposition coordination optimization solution is obtained.
5. The method for modeling and solving the optimal energy flow of the multi-energy coupling transmission and distribution network according to claim 1, wherein the system engineering theory in the step 3 is optimized hierarchically to enable the active distribution network to fully exert self-sufficient coordination capability and realize the consumption of internal distributed renewable energy sources, and methods for determining three operation modes of the active distribution network are provided, namely, firstly, the method for calculating the initial state is usedActive distribution network and transmission network connecting feeder power PLi,tIf P isLi,tIf the distributed energy resources in the active distribution network have surplus in power generation, the energy storage device is considered to be called for charging and storing, on the premise that safety constraint is met, if surplus still exists after the energy storage device is called fully, power is transmitted to the power transmission network on the premise that safety and bus power constraint conditions are met, and at the moment, the active distribution network is in a power supply state; if it is initially PLi,tIf the self-supply of the active distribution network cannot be realized, the discharging supply of the energy storage device is called, if the self-supply of the energy storage device cannot be realized after the self-supply is fully called, the electric power is absorbed from the power transmission network, and the active distribution network operates in a load state; if P is enabled in the initial state or by calling up the energy storage deviceLi,tAnd (3) considering that the active distribution network is an isolated island, decoupling operation optimization is carried out on the active distribution network and the power transmission network, and in the process, if complete consumption of renewable energy sources cannot be realized under the constraint condition, partial output is reduced by adopting an active management means.
6. The optimal energy flow modeling and solving method for the multi-energy coupling transmission network according to claim 1, characterized in that a multi-objective optimal energy flow algorithm based on approximate dynamic programming and NSGA-II is adopted in step 3, a Pareto optimal solution set is solved by adopting NSGA-II, for each non-dominated solution in the Pareto optimal solution set, a partial small fuzzy satisfaction calculation formula is adopted to calculate the satisfaction degree of each target value, then the comprehensive satisfaction degree of each non-dominated solution is calculated, and a path with the maximum accumulated satisfaction degree is selected through strategy iterative approximate dynamic programming to form a solution of the multi-objective dynamic optimal energy flow.
CN201610081125.0A 2016-02-04 2016-02-04 Optimal energy flow modeling and solving method for multi-energy coupling transmission and distribution network Active CN105703369B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610081125.0A CN105703369B (en) 2016-02-04 2016-02-04 Optimal energy flow modeling and solving method for multi-energy coupling transmission and distribution network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610081125.0A CN105703369B (en) 2016-02-04 2016-02-04 Optimal energy flow modeling and solving method for multi-energy coupling transmission and distribution network

Publications (2)

Publication Number Publication Date
CN105703369A CN105703369A (en) 2016-06-22
CN105703369B true CN105703369B (en) 2020-06-09

Family

ID=56222147

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610081125.0A Active CN105703369B (en) 2016-02-04 2016-02-04 Optimal energy flow modeling and solving method for multi-energy coupling transmission and distribution network

Country Status (1)

Country Link
CN (1) CN105703369B (en)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107332249A (en) * 2017-08-25 2017-11-07 长沙理工大学 Consider the fuzzy optimal trend method of the multiple target dynamic random of power system containing wind-powered electricity generation of transmission & distribution collaboration
CN107591806B (en) * 2017-10-09 2019-10-25 华中科技大学 A kind of major network dispatching method a few days ago considering the active regulating power of active distribution network
CN109039660B (en) * 2018-08-05 2020-11-20 广西师范大学 Energy sharing and power distribution method of EH distributed base station system
CN109193666B (en) * 2018-10-12 2021-07-30 南方电网科学研究院有限责任公司 General energy bus-based comprehensive energy system time sequence energy flow calculation method
CN109858774B (en) * 2019-01-09 2020-10-27 燕山大学 Source network load planning method for improving system safety and coordination
CN110070213B (en) * 2019-03-28 2022-05-20 广东工业大学 Day-ahead scheduling method of electricity-gas comprehensive energy system
CN110210747B (en) * 2019-05-28 2022-07-29 河海大学 Flexible scheduling method for electric-heating-gas interconnection energy system
CN110797911A (en) * 2019-10-31 2020-02-14 上海电力大学 Multi-target active power distribution network switch state adjustment method based on state splitting
CN110929405B (en) * 2019-11-28 2023-08-04 国网辽宁省电力有限公司经济技术研究院 Electro-pneumatic dynamic analysis method considering wind turbine unit and gas turbine unit
CN111612248B (en) * 2020-05-20 2023-09-08 云南电网有限责任公司电力科学研究院 Power distribution network side source-load coordination method and system
CN112100871B (en) * 2020-11-19 2021-02-19 清华四川能源互联网研究院 Decoupling method and device of multi-energy coupling system, electronic device and storage medium
CN112396238A (en) * 2020-11-24 2021-02-23 华中科技大学 Real-time optimized scheduling method for thermoelectric combined system
CN113051752B (en) * 2021-03-22 2023-04-07 西安热工研究院有限公司 Method for determining optimal heat source of high-pressure air energy storage system electrically coupled with coal
CN115147016B (en) * 2022-09-01 2023-01-24 国网浙江省电力有限公司宁波供电公司 Distributed power supply and demand side resource combined planning method considering low-carbon benefits
CN115454614B (en) * 2022-11-14 2023-03-24 安徽大学 Intelligent scheduling method for robot cluster energy supply
CN116722657B (en) * 2023-08-10 2023-10-20 福建合众鑫智能科技有限公司 Distributed energy grid-connected management system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103887790A (en) * 2014-03-05 2014-06-25 广州供电局有限公司 Cooperative power flow control method and system for power transmission network and power distribution network
CN104866919A (en) * 2015-05-20 2015-08-26 天津大学 Multi-target planning method for power grid of wind farms based on improved NSGA-II

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103887790A (en) * 2014-03-05 2014-06-25 广州供电局有限公司 Cooperative power flow control method and system for power transmission network and power distribution network
CN104866919A (en) * 2015-05-20 2015-08-26 天津大学 Multi-target planning method for power grid of wind farms based on improved NSGA-II

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于NSGA-II的含风电场电力系统多目标调度计划研究;罗斌;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20131215;C042-507 *

Also Published As

Publication number Publication date
CN105703369A (en) 2016-06-22

Similar Documents

Publication Publication Date Title
CN105703369B (en) Optimal energy flow modeling and solving method for multi-energy coupling transmission and distribution network
Kong et al. Robust stochastic optimal dispatching method of multi-energy virtual power plant considering multiple uncertainties
Wang et al. Economic and efficient multi-objective operation optimization of integrated energy system considering electro-thermal demand response
CN111881616B (en) Operation optimization method of comprehensive energy system based on multi-main-body game
Moghaddam et al. Multi-operation management of a typical micro-grids using Particle Swarm Optimization: A comparative study
Li et al. A microgrids energy management model based on multi-agent system using adaptive weight and chaotic search particle swarm optimization considering demand response
Bidgoli et al. Multi-stage optimal scheduling of multi-microgrids using deep-learning artificial neural network and cooperative game approach
Li et al. A novel bi-level robust game model to optimize a regionally integrated energy system with large-scale centralized renewable-energy sources in Western China
CN107546743B (en) Distributed power flow optimization method for radial power distribution network
CN109980685A (en) It is a kind of to consider probabilistic active distribution network distributed optimization operation method
Ma et al. Decentralized and coordinated scheduling model of interconnected multi-microgrid based on virtual energy storage
Cavazzini et al. Techno-economic benefits deriving from optimal scheduling of a Virtual Power Plant: Pumped hydro combined with wind farms
Zhang et al. Multi-objective planning of energy hub on economic aspects and resources with heat and power sources, energizable, electric vehicle and hydrogen storage system due to uncertainties and demand response
CN113471976B (en) Optimal scheduling method based on multi-energy complementary micro-grid and active power distribution network
CN115204562B (en) Multi-energy sharing-based distributed collaborative optimization scheduling method and system for interconnected micro energy networks
Yang et al. Multi-Objective optimal scheduling of island microgrids considering the uncertainty of renewable energy output
Cao et al. Optimal operation of cold–heat–electricity multi-energy collaborative system based on price demand response
Zhao et al. Day-ahead robust optimal dispatch of integrated energy station considering battery exchange service
Si et al. Cost-effective operation of the urban energy system with variable supply and demand via coordination of multi-energy flows
Yu et al. Research on energy management of a virtual power plant based on the improved cooperative particle swarm optimization algorithm
Zhu et al. Optimal scheduling of a wind energy dominated distribution network via a deep reinforcement learning approach
Zhang et al. Research on two-level energy optimized dispatching strategy of microgrid cluster based on IPSO algorithm
Liang et al. Capacity configuration optimization of wind-solar combined power generation system based on improved grasshopper algorithm
Gao et al. Fuzzy decision-based optimal energy dispatch for integrated energy systems with energy storage
An et al. Real-time optimal operation control of micro energy grid coupling with electricity-thermal-gas considering prosumer characteristics

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant