CN111342450B - Robust energy management method considering uncertain photovoltaic and load for traction power supply system - Google Patents

Robust energy management method considering uncertain photovoltaic and load for traction power supply system Download PDF

Info

Publication number
CN111342450B
CN111342450B CN202010154115.1A CN202010154115A CN111342450B CN 111342450 B CN111342450 B CN 111342450B CN 202010154115 A CN202010154115 A CN 202010154115A CN 111342450 B CN111342450 B CN 111342450B
Authority
CN
China
Prior art keywords
power
traction
photovoltaic
robust
supply system
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
CN202010154115.1A
Other languages
Chinese (zh)
Other versions
CN111342450A (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN202010154115.1A priority Critical patent/CN111342450B/en
Publication of CN111342450A publication Critical patent/CN111342450A/en
Application granted granted Critical
Publication of CN111342450B publication Critical patent/CN111342450B/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
    • 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/008Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
    • 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/26Arrangements for eliminating or reducing asymmetry in polyphase networks
    • 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/28Arrangements for balancing of the load in a network by storage of energy
    • 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/28Arrangements for balancing of the load in a network by storage of energy
    • H02J3/32Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
    • 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/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J7/00Circuit arrangements for charging or depolarising batteries or for supplying loads from batteries
    • H02J7/34Parallel operation in networks using both storage and other dc sources, e.g. providing buffering
    • H02J7/345Parallel operation in networks using both storage and other dc sources, e.g. providing buffering using capacitors as storage or buffering devices
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J7/00Circuit arrangements for charging or depolarising batteries or for supplying loads from batteries
    • H02J7/34Parallel operation in networks using both storage and other dc sources, e.g. providing buffering
    • H02J7/35Parallel operation in networks using both storage and other dc sources, e.g. providing buffering with light sensitive cells
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/50Arrangements for eliminating or reducing asymmetry in polyphase networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a robust energy management method for a traction power supply system considering uncertain photovoltaic and load, which comprises the following steps of: 1. acquiring photovoltaic output data in a load process of a traction substation, and constructing an uncertain set; 2. establishing an objective function of a robust optimization model; 3. establishing a constraint condition of a robust optimization model, and linearizing the constraint condition of the optimization model; 4. establishing a two-stage robust optimization model according to the objective function obtained in the step (2) and the constraint condition obtained in the step (3); 5. solving the model obtained in the step 4 by using a column and constraint generation algorithm to obtain a traction substation power flow controller, a hybrid energy storage device and a photovoltaic optimal operation plan in the worst scene, namely completing robust energy management of a traction power supply system; the method overcomes the influence of uncertainty of photovoltaic and traction load on the operation of the traction power supply system, improves the operation economy of the traction substation in severe operation environment, and is closer to the reality.

Description

Robust energy management method considering uncertain photovoltaic and load for traction power supply system
Technical Field
The invention belongs to the field of energy management optimization of a traction power supply system, and particularly relates to a robust energy management method of the traction power supply system considering uncertain photovoltaic and load.
Background
Renewable energy sources in China develop rapidly in recent years, and photovoltaic power generation has the advantages of no pollution, no noise, small regional limitation and the like. Meanwhile, electrified railways represented by high-speed railways and heavy haul railways are constructed and operated in large scale in China, and the operation mileage of high-speed railways in China exceeds 3.5 kilometers by 2019. The electrified railway has the following features: 1) the regional distribution is wide, and the traffic network and the renewable energy network have more geographical intersections, such as Lanzhou New (Wulumangoqi) high-speed rails and Sichuan (Chengdu) Tibetan (Lasa) railways crossing southwest and northwest areas with rich wind and light resources; 2) the traction load is high in demand and has high consumption potential, and particularly, in weak areas of an external power supply, the access of renewable energy sources can also play an important supporting role. Therefore, the photovoltaic-based renewable energy is connected into the traction power supply system to promote local consumption and improve the permeability of the renewable energy.
On the other hand, the electric power department in China charges the electric charges for the electrified railway by adopting two large-scale industrial electricity-consumption power rates, wherein the basic electric charges are charged according to the basic capacity or the maximum demand. Maximum demand charging also sets a "threshold" for minimum installed capacity. From the statistics of railway operation departments, the electricity fee has become one of the main operation payment fees. If the hybrid energy storage system is accessed, the peak clipping and valley filling of the traction load can be effectively realized, the required electric charge and the electric charge are greatly reduced, the regenerative braking energy utilization of a traction power supply system is promoted, and the system operation efficiency is improved. In addition, the requirement of the system on the capacity of the power supply equipment can be reduced. Therefore, the photovoltaic access and hybrid energy storage system in the traction power supply system of the electrified railway are increasingly concerned.
However, the random fluctuating nature of the traction load and the distributed power output, represented by photovoltaic, present challenges to the operation of the traction power supply system. How to effectively deal with uncertain factors in a traction power supply system and realize safe, high-quality and high-efficiency operation is the key for solving the problem of energy management of the traction power supply system.
Disclosure of Invention
In order to overcome the influence of photovoltaic and load uncertainty on a traction power supply system, the electricity charge expenditure cost of a railway department is reduced in the worst scene, and the energy management optimization method is closer to reality. The invention provides a robust energy management method of a traction power supply system considering uncertain photovoltaic and load.
The invention discloses a robust energy management method of a traction power supply system considering uncertain photovoltaic and load, which comprises the following steps:
step 1: acquiring predicted values and fluctuation intervals of traction substation load process data and photovoltaic output data, and constructing uncertain sets of traction load and photovoltaic output;
and 2, step: establishing a target function of a robust optimization model according to the electric charge parameters and the intervals of the load process data and the photovoltaic output data of the traction substation obtained in the step 1;
and step 3: according to power capacity parameters and three-phase voltage unbalance international parameters of a hybrid energy storage system and a photovoltaic system, establishing constraint conditions of a robust optimization model based on the traction substation load process data and the photovoltaic output data interval obtained in the step 1, and linearizing the constraint conditions of the robust optimization model;
and 4, step 4: establishing a robust energy management method of the traction power supply system based on a two-stage robust optimization model according to the objective function obtained in the step 2 and the constraint condition obtained in the step 3;
and 5: and (4) solving the model obtained in the step (4) by using a column and constraint generation algorithm to obtain the optimal charging and discharging power of the hybrid energy storage device, the optimal photovoltaic grid-connected power and the optimal power flow power of a back-to-back converter in the power flow controller under the worst scene, namely completing the robust energy management optimization of the traction power supply system.
Further, in the step 1, the load process data of the traction substation in the existing line can be obtained from the historical load data, and the traction substation in the newly-built line can be obtained by calculating load process simulation software, such as elbase/WEBANET, according to the high-speed railway line, the train and the schedule.
Further, the uncertain set of photovoltaic output and load in step 1 is:
Figure BDA0002403459970000021
Figure BDA0002403459970000022
Figure BDA0002403459970000023
in the formula: p, A and R are the uncertain set of photovoltaic power output, active load and reactive load, p, respectively t 、a t And r t For uncertain variables of photovoltaic contribution, active load and reactive load, Δ p t 、Δa t And Δ r t Respectively represent predicted values p t f 、a t f And r t f Deviation value of (d), parameter Γ p 、Γ a And Γ r Representing uncertainty margin, i.e. the total relative deviation of the uncertainty variable from the predicted value, which ranges from 0 to within a dayTotal time period N T In the meantime.
Further, the objective function in step 2 is:
Figure BDA0002403459970000024
Figure BDA0002403459970000025
in the formula: f is an objective function and represents daily operation cost of the traction substation, t is a time period and P is t grid,buy Electric power purchased from the grid for traction substations, P t grid,fed Electric power P fed back to the grid by the traction substation t PV For photovoltaic output, P t b,dis For discharging power of the battery, P t b,ch Charging power for batteries, P t u,dis Discharging power, P, for the super-capacitor t u,ch Charging power to the super capacitor, P t dem In order to demand the amount of power,
Figure BDA0002403459970000031
in order to achieve the cost price of photovoltaic operation and maintenance,
Figure BDA0002403459970000032
in order to keep the cost price of the battery running and maintaining,
Figure BDA0002403459970000033
cost price for operating and maintaining super capacitor, c dem The price of the electricity charge of the required quantity,
Figure BDA0002403459970000034
the price of the electricity purchased by the traction substation is,
Figure BDA0002403459970000035
charge price levied for feedback back to the grid power, N day For the number of days of operation per month,
Figure BDA0002403459970000036
and
Figure BDA0002403459970000037
are all binary variables;
wherein:
Figure BDA0002403459970000038
further, the constraint conditions in the step 3 comprise power balance constraint, hybrid energy storage system constraint, public power grid power constraint, photovoltaic power generation constraint, back-to-back converter constraint and three-phase voltage unbalance constraint; the method specifically comprises the following steps:
power balance constraint conditions:
P t grid,buy -P t grid,fed =P t T +P t α (6)
P t α +P t b,dis +P t u,dis +P t PV =P t β +P t b,ch +P t u,ch (7)
P t T +P t β =a t (8)
Figure BDA0002403459970000039
Figure BDA00024034599700000310
Figure BDA00024034599700000311
in the formula: p t T For active power of traction transformers, P t α Is alpha-phase active power, P, of a back-to-back converter t β Is a back-to-back converterActive power of beta phase, Q t β Is beta-phase reactive power, P, of a back-to-back converter t PV The photovoltaic output is obtained;
Figure BDA00024034599700000312
is the maximum limit value of the interaction power between the traction substation and the power grid,
Figure BDA00024034599700000313
is a binary variable representing the power direction of interaction between a traction substation and a power grid,
Figure BDA00024034599700000314
indicating that the interactive power flows from the grid to the traction substation,
Figure BDA00024034599700000315
and the representation interactive power is fed back to the power grid by the traction substation.
Hybrid energy storage system constraint conditions:
Figure BDA00024034599700000316
Figure BDA00024034599700000317
Figure BDA00024034599700000318
Figure BDA00024034599700000319
Figure BDA00024034599700000320
Figure BDA0002403459970000041
Figure BDA0002403459970000042
in the formula: epsilon b Is the self-discharge rate of the battery, epsilon b Is the self-discharge rate of the super capacitor, eta b,dis Eta, discharge efficiency of the cell b,ch Is the charging efficiency of the battery, eta u,dis Is the discharge efficiency of the super capacitor, eta u,ch To the charging efficiency of the super capacitor, Δ t is the unit time period,
Figure BDA0002403459970000043
for the electric energy stored in the battery during the time period t +1, E t b Storing the electric energy for the battery in the time period t;
Figure BDA0002403459970000044
the electric energy stored by the super capacitor in the time period of t +1,
Figure BDA0002403459970000045
the electric energy stored for the super capacitor in the time period t;
Figure BDA0002403459970000046
is the rated power of the battery and is,
Figure BDA0002403459970000047
the power of the super capacitor is rated,
Figure BDA0002403459970000048
is the minimum state of charge of the battery,
Figure BDA0002403459970000049
the maximum state of charge of the battery is,
Figure BDA00024034599700000410
is the rated capacity of the battery,
Figure BDA00024034599700000411
the capacity of the super capacitor is the rated capacity,
Figure BDA00024034599700000412
the electrical energy stored by the battery for the time period t-1,
Figure BDA00024034599700000413
the electric energy stored by the super capacitor in the t-1 period,
Figure BDA00024034599700000414
the minimum state of charge of the super capacitor is obtained,
Figure BDA00024034599700000415
the maximum charge state of the super capacitor;
Figure BDA00024034599700000416
the electrical energy stored in the battery for the initial period of time each day,
Figure BDA00024034599700000417
the electrical energy stored in the battery for the last period of the day,
Figure BDA00024034599700000418
for the purpose of the initial state of charge per day,
Figure BDA00024034599700000419
the super capacitor stores electric energy for the initial time period every day,
Figure BDA00024034599700000420
the electric energy stored by the super capacitor for the last time period of each day,
Figure BDA00024034599700000421
the initial charge state of the super capacitor every day;
Figure BDA00024034599700000422
and
Figure BDA00024034599700000423
are all binary variables.
Photovoltaic power generation constraint:
0≤P t PV ≤p t (19)
in the formula: p is a radical of t The photovoltaic output uncertain variable is the solar photovoltaic output upper limit value.
Back-to-back converter constraint:
Figure BDA00024034599700000424
Figure BDA00024034599700000425
in the formula:
Figure BDA00024034599700000426
the capacity of the alpha phase of the back-to-back converter,
Figure BDA00024034599700000427
the capacity of the beta-phase of the back-to-back converter.
And (3) three-phase voltage unbalance degree constraint:
Figure BDA00024034599700000428
Figure BDA00024034599700000429
Figure BDA00024034599700000430
in the formula: epsilon U For traction substation power grid side three-phase voltage unbalance degree, U S For the grid side line voltage, S cap For the short-circuit capacity of the side line of the power grid,
Figure BDA0002403459970000051
is the upper limit value of the unbalance degree of the three-phase voltage in the national standard,
Figure BDA0002403459970000052
for grid side negative sequence current, U T For the voltage at the outlet of the traction transformer, U α Is the voltage at the alpha-phase outlet of the back-to-back converter, N 1 For single-phase traction transformer transformation ratio, N 2 For high voltage matching transformer transformation ratio, a is complex operator e j120°
Figure BDA0002403459970000053
Is the voltage-current phase angle difference of the single-phase traction transformer,
Figure BDA0002403459970000054
is the voltage-current phase angle difference of alpha phases of the back-to-back converter, I T For drawing the transformer current, I α Is the current of the alpha phase of the back-to-back converter.
Further, the constraint condition linearization method in step 3 is as follows:
the max function in equation (5) is linearized by the following equation:
max(P t dem )=P dem,max (25)
Figure BDA0002403459970000055
in the formula: p is dem,max Is an auxiliary variable representing the maximum demand value during the day.
The formula (21) after linearization is given by:
Figure BDA0002403459970000056
Figure BDA0002403459970000057
Figure BDA0002403459970000058
Figure BDA0002403459970000059
in the formula: n is a radical of p The number of the sectors is divided into equal parts of a PQ semicircle, and the Q is more than or equal to 0; delta theta is the sector angle (P) k ,Q k ) The coordinates of the division points of the fan shape and the PQ semicircle.
Equation (24) is linearized as follows:
Figure BDA00024034599700000510
Figure BDA00024034599700000511
Figure BDA00024034599700000512
in the formula:
Figure BDA00024034599700000513
and
Figure BDA00024034599700000514
are all auxiliary variables, and are all the auxiliary variables,
Figure BDA00024034599700000515
is a binary variable.
Further, the two-stage robust optimization model for robust energy management of the traction power supply system established in step 4 is as follows:
Figure BDA00024034599700000516
Figure BDA00024034599700000517
in the formula: x represents the binary decision variable vector of the first stage,
Figure BDA0002403459970000061
y denotes a second stage continuous decision variable vector,
Figure BDA0002403459970000062
Figure BDA0002403459970000063
c. b, D, D, E, E, F, F, G, H and I are all parameter matrixes or parameter vectors.
Further, in step 5, based on the robust energy management model of the traction power supply system obtained in step 4, a main model and a sub model of the robust energy management model of the traction power supply system are formed through a column and constraint generation algorithm, and finally the main model and the sub model are solved in an iterative loop manner until a convergence standard is met.
Wherein, the expression of the main model is as follows:
Figure BDA0002403459970000064
s.t.x∈{0,1} (37)
Figure BDA0002403459970000065
Figure BDA0002403459970000066
Figure BDA0002403459970000067
Figure BDA0002403459970000068
Figure BDA0002403459970000069
Figure BDA00024034599700000610
Figure BDA00024034599700000611
in the formula: k is the number of iterative solution times, y l The decision variables added to the main model at the first loop,
Figure BDA00024034599700000612
to solve for the worst photovoltaic output obtained by the submodel,
Figure BDA00024034599700000613
for solving the worst active load scenario, r, obtained for the submodel l * And the worst reactive load scene obtained by solving the submodel is shown.
The sub-model expression is:
Figure BDA00024034599700000614
s.t.By≤d,(γ 1 ) (46)
Dy=e,(γ 2 ) (47)
Fy≤f-Ex * ,(γ 3 ) (48)
Gy≤p,(γ 4 ) (49)
Hy=a,(γ 5 ) (50)
Iy=r,(γ 6 ) (51)
in the formula: x is the number of * Is mainly composed ofModel optimal solution, { gamma { 123456 Is a dual variable of the constraint.
The submodel equivalent representation method is as follows:
worst scenario p in sub-model * 、a * And r * To not determine the extreme values in sets P, A and R, equations (1) - (3) are therefore equivalent to:
Figure BDA0002403459970000071
Figure BDA0002403459970000072
Figure BDA0002403459970000073
in the formula: g t 、m t
Figure BDA0002403459970000074
And
Figure BDA0002403459970000075
are all binary variables.
Based on strong dual theory, submodels (45) - (51) are equivalent to:
Figure BDA0002403459970000076
-B T γ 1 +D T γ 2 -F T γ 3 -G T γ 4 +H T γ 5 +I T γ 6 =c T (56)
γ 1 ≥0,γ 3 ≥0,γ 4 ≥0,γ 256 is a free variable (57)
Figure BDA0002403459970000077
Figure BDA0002403459970000078
Figure BDA0002403459970000081
Figure BDA0002403459970000082
In the formula: lambda, mu, omega 1 And ω 2 Are all auxiliary variables, and are all the auxiliary variables,
Figure BDA0002403459970000083
the invention has the beneficial effects that:
1. the robust energy management optimization method for the traction power supply system can overcome the influence of uncertainty of photovoltaic and traction load on safe and efficient operation of the traction power supply system, effectively reduces the electric charge expenditure of railway departments, and accords with the development trend of green and intelligent traffic.
2. The robust energy management optimization model of the traction power supply system is established based on the uncertain set of photovoltaic output and traction load, the robustness of the energy management model of the traditional traction power supply system is improved, and the optimality of the traction power supply system in the worst scene is ensured;
3. according to the method, the nonlinear elements in the objective function and the constraint condition are subjected to linearization processing, the mixed integer linear programming model is established, direct solving by using an optimization solver is facilitated, and the complexity of solving the mixed integer nonlinear model is avoided.
Drawings
Fig. 1 is a schematic diagram of a robust energy management structure of a traction power supply system in the invention.
FIG. 2 is a schematic flow chart of the present invention.
Detailed Description
The invention is further described with reference to the following figures and specific embodiments.
The structure of a traction power supply system aimed at by the invention is shown in fig. 1, and the flow of the robust energy management method of the traction power supply system considering uncertain photovoltaic and load is shown in fig. 2, specifically:
step 1: based on actual measurement historical data, obtaining predicted values and fluctuation intervals of load process data and photovoltaic output data of the traction substation, and establishing an uncertain set of traction load and photovoltaic output:
Figure BDA0002403459970000084
Figure BDA0002403459970000085
Figure BDA0002403459970000091
in the formula: p, A and R are the uncertain set of photovoltaic power output, active load and reactive load, p, respectively t 、a t And r t For uncertain variables of photovoltaic output, active and reactive loads, Δ p t 、Δa t And Δ r t Respectively represent predicted values
Figure BDA0002403459970000092
And
Figure BDA00024034599700000913
deviation value of, parameter Γ p 、Γ a And Γ r Representing uncertainty margin, i.e. the total relative deviation of the uncertainty variable from the predicted value, which ranges from 0 to the total time period N of the day T In between.
And 2, step: establishing an objective function of an optimization model according to the electric charge parameters and the intervals of the load process data and the photovoltaic output data of the traction substation obtained in the step 1;
the target function expression is:
Figure BDA0002403459970000093
Figure BDA0002403459970000094
in the formula: f is an objective function and represents daily operation cost of the traction substation, t is a time period and P is t grid,buy Electric power purchased from the grid for traction substations, P t grid,fed Electric power P fed back to the grid by the traction substation t PV For photovoltaic output, P t b,dis For discharging power of the battery, P t b,ch Charging power for batteries, P t u,dis Discharging power, P, for the super-capacitor t u,ch Charging power to the super capacitor, P t dem In order to demand the power, the power supply is,
Figure BDA0002403459970000095
in order to achieve the cost price of photovoltaic operation and maintenance,
Figure BDA0002403459970000096
in order to maintain the cost price of the battery operation,
Figure BDA0002403459970000097
cost price for operating and maintaining super capacitor, c dem The price of the required electricity charge is calculated,
Figure BDA0002403459970000098
the price of the electricity purchased by the traction substation is low,
Figure BDA0002403459970000099
charge price levied for feedback back to the grid power, N day For the number of days of operation per month,
Figure BDA00024034599700000910
and
Figure BDA00024034599700000911
are all binary variables;
wherein:
Figure BDA00024034599700000912
and 3, step 3: establishing constraint conditions of an optimization model based on the traction substation load process data and the photovoltaic output data interval obtained in the step 1 according to power capacity parameters and three-phase voltage unbalance international parameters of the hybrid energy storage system and the photovoltaic system, and linearizing the constraint conditions of the optimization model;
the constraint conditions comprise power balance constraint, hybrid energy storage system constraint, public power grid power constraint, photovoltaic power generation constraint, back-to-back converter constraint in a power flow controller and three-phase unbalance constraint.
The constraints are as follows:
power balance constraint conditions:
P t grid,buy -P t grid,fed =P t T +P t α (6)
P t α +P t b,dis +P t u,dis +P t PV =P t β +P t b,ch +P t u,ch (7)
P t T +P t β =a t (8)
Figure BDA0002403459970000101
Figure BDA0002403459970000102
Figure BDA0002403459970000103
in the formula: p t T For active power of traction transformers, P t α Is alpha-phase active power, P, of a back-to-back converter t β Is the beta-phase active power of the back-to-back converter,
Figure BDA0002403459970000104
is beta-phase reactive power P of a back-to-back converter t PV The photovoltaic output is obtained;
Figure BDA0002403459970000105
is the maximum limit value of the interaction power between the traction substation and the power grid,
Figure BDA0002403459970000106
is a binary variable representing the direction of power interaction between the traction substation and the grid,
Figure BDA0002403459970000107
indicating that the interactive power flows from the grid to the traction substation,
Figure BDA0002403459970000108
and the representation interactive power is fed back to the power grid by the traction substation.
Hybrid energy storage system constraint conditions:
Figure BDA0002403459970000109
Figure BDA00024034599700001010
Figure BDA00024034599700001011
Figure BDA00024034599700001012
Figure BDA00024034599700001013
Figure BDA00024034599700001014
Figure BDA00024034599700001015
in the formula: epsilon b Is the self-discharge rate of the battery, epsilon b Is the self-discharge rate of the super capacitor, eta b,dis Is the discharge efficiency of the battery, η b,ch Is the charging efficiency of the battery, eta u,dis Is the discharge efficiency of the super capacitor, eta u,ch To the charging efficiency of the super capacitor, Δ t is the unit time period,
Figure BDA00024034599700001016
for the electrical energy stored by the battery during the time period t +1,
Figure BDA00024034599700001017
storing the electric energy for the battery in the time period t;
Figure BDA00024034599700001018
the electric energy stored by the super capacitor in the time period of t +1,
Figure BDA00024034599700001019
the electric energy stored for the super capacitor in the time period t;
Figure BDA00024034599700001020
the power of the battery is rated for the rated power,
Figure BDA00024034599700001021
the power of the super capacitor is rated,
Figure BDA00024034599700001022
is the minimum state of charge of the battery,
Figure BDA00024034599700001023
the maximum state of charge of the battery is,
Figure BDA00024034599700001024
is the rated capacity of the battery and is,
Figure BDA00024034599700001025
the capacity of the super capacitor is the rated capacity,
Figure BDA00024034599700001026
the electrical energy stored by the battery for the time period t-1,
Figure BDA00024034599700001027
the electric energy stored by the super capacitor in the t-1 period,
Figure BDA00024034599700001028
the minimum state of charge of the super capacitor is obtained,
Figure BDA00024034599700001029
the maximum charge state of the super capacitor;
Figure BDA00024034599700001030
the electrical energy stored in the battery for the initial period of time each day,
Figure BDA0002403459970000111
the electrical energy stored in the battery for the last period of the day,
Figure BDA0002403459970000112
for the purpose of the initial state of charge per day,
Figure BDA0002403459970000113
the electric energy stored by the super capacitor for the initial time period every day,
Figure BDA0002403459970000114
the electric energy stored by the super capacitor for the last time period of each day,
Figure BDA0002403459970000115
the initial charge state of the super capacitor every day;
Figure BDA0002403459970000116
and
Figure BDA0002403459970000117
are all binary variables.
Photovoltaic power generation constraint:
0≤P t PV ≤p t (19)
in the formula: p is a radical of t The photovoltaic output uncertain variable is the solar photovoltaic output upper limit value.
Back-to-back converter constraint:
Figure BDA0002403459970000118
Figure BDA0002403459970000119
in the formula:
Figure BDA00024034599700001110
the capacity of the alpha phase of the back-to-back converter,
Figure BDA00024034599700001111
is the capacity of the beta-phase of the back-to-back converter.
And (3) three-phase voltage unbalance degree constraint:
Figure BDA00024034599700001112
Figure BDA00024034599700001113
Figure BDA00024034599700001114
in the formula: epsilon U For traction substation power grid side three-phase voltage unbalance degree, U S For the grid side line voltage, S cap For the short-circuit capacity of the side line of the power grid,
Figure BDA00024034599700001115
is the upper limit value of the unbalance degree of the three-phase voltage in the national standard,
Figure BDA00024034599700001116
for grid side negative sequence current, U T For the voltage at the outlet of the traction transformer, U α Is the voltage at the alpha-phase outlet of the back-to-back converter, N 1 For single-phase traction transformer transformation ratio, N 2 For high voltage matching transformer transformation ratio, a is complex operator e j120°
Figure BDA00024034599700001117
Is the voltage-current phase angle difference of the single-phase traction transformer,
Figure BDA00024034599700001118
is the voltage-current phase angle difference of the alpha phase of the back-to-back converter, I T For drawing the transformer current, I α Is the current of the alpha phase of the back-to-back converter.
The constraint linearization method is as follows:
the max function in equation (5) is linearized by the following equation:
max(P t dem )=P dem,max (25)
Figure BDA00024034599700001119
in the formula: p is dem,max Is an auxiliary variable representing the maximum demand value during the day.
The formula (21) is linearized by the following equation:
Figure BDA0002403459970000121
Figure BDA0002403459970000122
Figure BDA0002403459970000123
Figure BDA0002403459970000124
in the formula: n is a radical of p The number of the fan-shaped PQ semicircles is equal to that of the fan-shaped PQ semicircles, and the Q is more than or equal to 0; delta theta is the sector angle, (P) k ,Q k ) Is the division point coordinate of the fan shape and the PQ semicircle.
Equation (24) is linearized as follows:
Figure BDA0002403459970000125
Figure BDA0002403459970000126
Figure BDA0002403459970000127
in the formula:
Figure BDA0002403459970000128
and
Figure BDA0002403459970000129
are all auxiliary variables, and are all the auxiliary variables,
Figure BDA00024034599700001210
is a binary variable.
And 4, step 4: establishing a two-stage robust optimization model of robust energy management of the traction power supply system according to the objective function obtained in the step (2) and the constraint condition obtained in the step (3);
the two-stage robust optimization model for robust energy management of the traction power supply system comprises the following steps:
Figure BDA00024034599700001211
Figure BDA00024034599700001212
in the formula: x represents the binary decision variable vector of the first stage,
Figure BDA00024034599700001213
y denotes a second stage continuous decision variable vector,
Figure BDA00024034599700001214
Figure BDA00024034599700001215
c. b, D, D, E, E, F, F, G, H and I are parameter matrixes or parameter vectors.
And 5: and forming a main model and a sub model of the robust energy management model of the traction power supply system through a column and constraint generation algorithm, and finally, circularly and iteratively solving the main model and the sub model by using optimization software (such as a hybrid integer optimization solver CPLEX in a Matlab environment) to obtain the optimal charging and discharging power of the hybrid energy storage device, the optimal photovoltaic grid-connected power and the optimal power flow power of a back-to-back converter in a power flow controller under the worst scene, namely completing the robust energy management optimization of the traction power supply system.
The main model expression is:
Figure BDA0002403459970000131
s.t.x∈{0,1} (37)
Figure BDA0002403459970000132
Figure BDA0002403459970000133
Figure BDA0002403459970000134
Figure BDA0002403459970000135
Figure BDA0002403459970000136
Figure BDA0002403459970000137
Figure BDA0002403459970000138
in the formula: k is the number of iterative solution times, y l The decision variables added to the main model at the first loop,
Figure BDA0002403459970000139
to solve for the worst photovoltaic output obtained by the submodel,
Figure BDA00024034599700001310
for solving the worst active load scenario, r, obtained for the submodel l * Representing the worst reactive load scene obtained by solving the submodel;
the sub-model expression is:
Figure BDA00024034599700001311
s.t.By≤d,(γ 1 ) (46)
Dy=e,(γ 2 ) (47)
Fy≤f-Ex * ,(γ 3 ) (48)
Gy≤p,(γ 4 ) (49)
Hy=a,(γ 5 ) (50)
Iy=r,(γ 6 ) (51)
in the formula: x is the number of * For the optimal solution of the main model, { gamma { 123456 Is a constraint dual variable;
the submodel equivalent representation method is as follows:
worst scene p in sub-model * 、a * And r * To not determine the extreme values in sets P, A and R, equations (1) - (3) are therefore equivalent to:
Figure BDA0002403459970000141
Figure BDA0002403459970000142
Figure BDA0002403459970000143
in the formula: g t 、m t
Figure BDA0002403459970000144
And
Figure BDA0002403459970000145
are all binary variables;
based on strong dual theory, submodels (45) - (51) are equivalent to:
Figure BDA0002403459970000146
-B T γ 1 +D T γ 2 -F T γ 3 -G T γ 4 +H T γ 5 +I T γ 6 =c T (56)
γ 1 ≥0,γ 3 ≥0,γ 4 ≥0,γ 256 is a free variable (57)
Figure BDA0002403459970000147
Figure BDA0002403459970000148
Figure BDA0002403459970000149
Figure BDA00024034599700001410
In the formula: lambda, mu, omega 1 And omega 2 Are all auxiliary variables, and are all the auxiliary variables,
Figure BDA00024034599700001411
examples
The robust energy management structure of the electrified railway traction power supply system considering uncertain photovoltaic and load is shown in figure 1, the parameters of the energy storage system are shown in table 1, the parameters of a traction substation are shown in table 2, and the power price parameters of a power grid are shown in table 3.
TABLE 1 hybrid energy storage System parameters
Figure BDA0002403459970000151
TABLE 2 traction substation parameters
Figure BDA0002403459970000152
TABLE 3 Electricity price parameter
Figure BDA0002403459970000153
The operating and maintenance costs of other parameters, photovoltaic, battery and super capacitor are all 0.1 rmw/kWh, and the charging c of feedback electric energy fed =-0.8c buy In linearization of PFC converter, the fan angle Δ θ is 30 °, C&The CG algorithm iteration convergence precision epsilon is 0.01. The uncertain interval of the photovoltaic output and the traction load is expressed as [ 1-lambda ] p ,1+λ p ]p f 、[1-λ a ,1+λ a ]a f And [ 1-lambda ] r ,1+λ r ]r f Wherein λ is p 、λ a And λ r For predicting the deviation coefficients, 0.1 is taken.
The robust energy management optimization method of the traction power supply system is compared with the traditional deterministic energy management method, and the results are shown in tables 4 and 5 after simulation calculation.
TABLE 4 comparison of results of deterministic method and robust method under different uncertainty margins
Figure BDA0002403459970000154
TABLE 5 comparison of deterministic and robust method results under different prediction biases
Figure BDA0002403459970000161
As can be seen from tables 4 and 5, with the increase of the uncertainty margin and the prediction deviation coefficient, compared with the conventional deterministic energy management method, the robust energy management optimization method for the traction power supply system has the advantages that the total operating cost is lower in the worst scenario, and the saving rate is continuously increased. Particularly, when the prediction deviation coefficient reaches 0.6, the traditional deterministic energy management method has the situation of being unsolvable, so the robustness is low, and the robust energy management optimization method can be solvated in all scenes, so the robustness is high.
The invention provides a double-stage robust model-based energy management robust optimization method for a traction power supply system, which aims at the traction power supply system connected with a photovoltaic power generation system and an energy storage device, takes photovoltaic output and traction load uncertainty into consideration, and provides the double-stage robust model-based energy management robust optimization method for the traction power supply system. According to the method, in the first stage, based on photovoltaic output and traction load prediction information, a charging and discharging strategy of an energy storage device and an electric quantity trading scheme with a power grid are formulated, and in the second stage, the worst scene of the photovoltaic output and traction load and the corresponding optimal power flow are found. The column and constraint generation algorithm is adopted to solve the two-stage robust optimization model, and the result shows that compared with the traditional deterministic energy management method, the robust energy management optimization method has the advantages of optimality and robustness, especially under the condition that the uncertainty of the operating environment is increased. The method makes the energy management method of the traction power supply system more practical and provides a foundation for the access and engineering application of renewable energy sources and energy storage systems in future electrified railways.

Claims (7)

1. The method for managing the robust energy of the traction power supply system in consideration of uncertain photovoltaic and load is characterized by comprising the following steps of:
step 1: acquiring predicted values and fluctuation intervals of traction substation load process data and photovoltaic output data, and constructing uncertain sets of traction load and photovoltaic output;
step 2: establishing a target function of a robust optimization model according to the electric charge parameters and the intervals of the load process data and the photovoltaic output data of the traction substation obtained in the step 1;
and step 3: according to power capacity parameters of a hybrid energy storage system and a photovoltaic system and national standard limit values of three-phase voltage unbalance degrees, establishing constraint conditions of a robust optimization model based on traction substation load process data and photovoltaic output data intervals obtained in the step 1, and linearizing the constraint conditions of the robust optimization model;
and 4, step 4: establishing a robust energy management method of the traction power supply system based on a two-stage robust optimization model according to the objective function obtained in the step 2 and the constraint condition obtained in the step 3;
and 5: and (4) solving the model obtained in the step (4) by using a column and constraint generation algorithm to obtain the optimal charging and discharging power of the hybrid energy storage device, the optimal photovoltaic grid-connected power and the optimal power flow power of a back-to-back converter in the power flow controller under the worst scene, namely completing the robust energy management optimization of the traction power supply system.
2. The method for robust energy management of a traction power supply system considering uncertain photovoltaic and load according to claim 1, wherein the uncertain set of photovoltaic output and load in step 1 is:
Figure FDA0003637731690000011
Figure FDA0003637731690000012
Figure FDA0003637731690000013
in the formula: p, A and R are the uncertain sets of photovoltaic power output, active load and reactive load, p, respectively t 、a t And r t For uncertain variables of photovoltaic output, active and reactive loads, Δ p t 、Δa t And Δ r t Respectively represent predicted values
Figure FDA0003637731690000014
And r t f Deviation value of, parameter Γ p 、Γ a And Γ r Representing uncertainty margin, i.e. the total relative deviation of the uncertainty variable from the predicted value, which ranges from 0 to the total time period N of the day T In the meantime.
3. The method for robust energy management of a traction power supply system considering uncertain photovoltaic and load according to claim 2, wherein the objective function in the step 2 is as follows:
Figure FDA0003637731690000015
Figure FDA0003637731690000021
in the formula: f is an objective function and represents daily operation cost of the traction substation, t is a time period and P is t grid,buy Electric power purchased from the grid for traction substations, P t grid,fed Electric power P fed back to the grid by the traction substation t PV For photovoltaic output, P t b,dis For discharging power of the battery, P t b,ch Charging power for batteries, P t u,dis Discharging power, P, for the super-capacitor t u,ch Charging power to the super capacitor, P t dem In order to demand the amount of power,
Figure FDA0003637731690000022
in order to reduce the cost price of photovoltaic operation and maintenance,
Figure FDA0003637731690000023
in order to maintain the cost price of the battery operation,
Figure FDA0003637731690000024
cost price for operating and maintaining super capacitor, c dem The price of the electricity charge of the required quantity,
Figure FDA0003637731690000025
the price of the electricity purchased by the traction substation is,
Figure FDA0003637731690000026
charge price levied for feedback back to the grid power, N day For the number of days of operation per month,
Figure FDA0003637731690000027
and
Figure FDA0003637731690000028
are all binary variables; p t T For active power of traction transformers, P t α For back-to-back converter alpha-phase active power, P t β Is beta-phase active power of a back-to-back converter,
Figure FDA0003637731690000029
is beta-phase reactive power of a back-to-back converter; p dem,max Is an auxiliary variable representing the maximum demand value during the day;
Figure FDA00036377316900000210
and
Figure FDA00036377316900000211
are all auxiliary variables;
wherein:
Figure FDA00036377316900000212
4. the method for robust energy management of a traction power supply system considering uncertain photovoltaic and load according to claim 3, wherein the constraint conditions in the step 3 are as follows:
power balance constraint conditions:
P t grid,buy -P t grid,fed =P t T +P t α (6)
P t α +P t b,dis +P t u,dis +P t PV =P t β +P t b,ch +P t u,ch (7)
P t T +P t β =a t (8)
Figure FDA00036377316900000213
Figure FDA00036377316900000214
Figure FDA00036377316900000215
in the formula: p t T For active power of traction transformers, P t α For back-to-back converter alpha-phase active power, P t β Is the beta-phase active power of the back-to-back converter,
Figure FDA00036377316900000216
is beta-phase reactive power, P, of a back-to-back converter t PV The magnitude of photovoltaic output is obtained;
Figure FDA00036377316900000217
is the maximum limit value of the interaction power between the traction substation and the power grid,
Figure FDA00036377316900000218
is a binary variable representing the power direction of interaction between a traction substation and a power grid,
Figure FDA00036377316900000219
indicating that the interactive power flows from the grid to the traction substation,
Figure FDA00036377316900000220
representing that the interactive power is fed back to the power grid by the traction substation;
hybrid energy storage system constraint conditions:
Figure FDA00036377316900000221
Figure FDA0003637731690000031
Figure FDA0003637731690000032
Figure FDA0003637731690000033
Figure FDA0003637731690000034
Figure FDA0003637731690000035
Figure FDA0003637731690000036
in the formula: epsilon b Is the self-discharge rate of the battery, epsilon u Is the self-discharge rate of the super capacitor, eta b,dis Is the discharge efficiency of the battery, η b,ch Is the charging efficiency of the battery, eta u,dis Is the discharge efficiency of the super capacitor, eta u,ch To the charging efficiency of the super capacitor, Δ t is the unit time period,
Figure FDA0003637731690000037
for the electrical energy stored by the battery during the time period t +1,
Figure FDA0003637731690000038
storing the electric energy for the battery in the time period t;
Figure FDA0003637731690000039
the electric energy stored by the super capacitor in the time period of t +1,
Figure FDA00036377316900000330
the electric energy stored for the super capacitor in the time period t;
Figure FDA00036377316900000310
is the rated power of the battery and is,
Figure FDA00036377316900000311
the power of the super capacitor is rated, b SOCis the minimum state of charge of the battery,
Figure FDA00036377316900000312
the maximum state of charge of the battery is,
Figure FDA00036377316900000313
is the rated capacity of the battery,
Figure FDA00036377316900000314
the capacity of the super capacitor is rated,
Figure FDA00036377316900000315
the electrical energy stored by the battery for the time period t-1,
Figure FDA00036377316900000316
the electric energy stored by the super capacitor in the period of t-1, u SOCis the minimum state of charge of the super capacitor,
Figure FDA00036377316900000317
the maximum charge state of the super capacitor;
Figure FDA00036377316900000318
the electrical energy stored in the battery for the initial period of the day,
Figure FDA00036377316900000319
the electrical energy stored in the battery for the last period of the day,
Figure FDA00036377316900000320
for the purpose of the initial state of charge per day,
Figure FDA00036377316900000321
the electric energy stored by the super capacitor for the initial time period every day,
Figure FDA00036377316900000322
the stored energy for the super capacitor for the last period of the day,
Figure FDA00036377316900000323
the initial charge state of the super capacitor every day;
Figure FDA00036377316900000324
and
Figure FDA00036377316900000325
are all binary variables;
photovoltaic power generation constraint:
0≤P t PV ≤p t (19)
in the formula: p is a radical of t The photovoltaic output uncertain variable is the solar photovoltaic output upper limit value;
back-to-back converter constraint:
Figure FDA00036377316900000326
Figure FDA00036377316900000327
in the formula:
Figure FDA00036377316900000328
the capacity of the alpha phase of the back-to-back converter,
Figure FDA00036377316900000329
the capacity of the beta phase of the back-to-back converter;
and (3) three-phase voltage unbalance degree constraint:
Figure FDA0003637731690000041
Figure FDA0003637731690000042
Figure FDA0003637731690000043
in the formula:ε U for the unbalance of three-phase voltage at the side of the power grid of the traction substation, U S For the grid side line voltage, S cap For the short-circuit capacity on the grid side,
Figure FDA0003637731690000044
is the upper limit value of the unbalance degree of the three-phase voltage in the national standard,
Figure FDA0003637731690000045
for grid side negative sequence current, U T For the voltage at the outlet of the traction transformer, U α Is the voltage at the alpha-phase outlet of the back-to-back converter, N 1 For single-phase traction transformer transformation ratio, N 2 For high voltage matching transformer transformation ratio, a is complex operator e j120°
Figure FDA00036377316900000413
Is the voltage-current phase angle difference of the single-phase traction transformer,
Figure FDA0003637731690000046
is the voltage-current phase angle difference of the alpha phase of the back-to-back converter, I T For drawing transformer currents, I α Is the current of the alpha phase of the back-to-back converter.
5. The method for robust energy management of a traction power supply system considering uncertain photovoltaic and load according to claim 4, wherein the constraint condition linearization method in the step 3 is as follows:
the max function in equation (5) is linearized by the following equation:
max(P t dem )=P dem,max (25)
Figure FDA0003637731690000047
in the formula: p is dem,max Is an auxiliary variable representing the maximum demand value during the day;
the formula (21) is linearized by the following equation:
Figure FDA0003637731690000048
Figure FDA0003637731690000049
Figure FDA00036377316900000410
Figure FDA00036377316900000411
in the formula: n is a radical of p The number of the fan-shaped PQ semicircles is equal to that of the fan-shaped PQ semicircles, and the Q is more than or equal to 0; delta theta is the sector angle, (P) k ,Q k ) The coordinates of the division points of the fan shape and the PQ semicircle;
equation (24) is linearized as follows:
Figure FDA00036377316900000412
Figure FDA0003637731690000051
Figure FDA0003637731690000052
in the formula:
Figure FDA0003637731690000053
and
Figure FDA0003637731690000054
are all auxiliary variables, and are all the auxiliary variables,
Figure FDA0003637731690000055
is a binary variable.
6. The method for robust energy management of a traction power supply system considering uncertain photovoltaic and load according to claim 5, wherein the two-stage robust optimization model for robust energy management of the traction power supply system established in the step 4 is as follows:
Figure FDA0003637731690000056
Figure FDA0003637731690000057
in the formula: x represents the binary decision variable vector of the first stage,
Figure FDA0003637731690000058
y represents the second stage continuous decision variable vector,
Figure FDA0003637731690000059
Figure FDA00036377316900000510
c. b, D, D, E, E, F, F, G, H and I are parameter matrixes or parameter vectors.
7. The method for managing the robust energy of the traction power supply system in consideration of the uncertain photovoltaic and load according to claim 6, wherein a main model and a sub model of the robust energy management model of the traction power supply system are formed through a column and constraint generation algorithm in the step 5, the main model and the sub model are solved in a circulating iteration mode to obtain the optimal charging and discharging power of the hybrid energy storage device in the worst scene, the optimal photovoltaic grid-connected power and the optimal power flow power of a back-to-back converter in a power flow controller, and therefore the robust energy management optimization of the traction power supply system is completed;
wherein, the expression of the main model is as follows:
Figure FDA00036377316900000511
s.t.x∈{0,1} (37)
Figure FDA00036377316900000512
Figure FDA00036377316900000513
Figure FDA00036377316900000514
Figure FDA00036377316900000515
Figure FDA00036377316900000516
Figure FDA0003637731690000061
Figure FDA0003637731690000062
in the formula: k is the number of iterative solution times, y l Is a first circulation time direction main modelThe decision variables that are added are the variables of the decision,
Figure FDA0003637731690000063
to solve for the worst photovoltaic output obtained by the submodel,
Figure FDA0003637731690000064
for solving the worst active load scenario, r, obtained for the submodel l * Representing the worst reactive load scene obtained by solving the submodel;
the sub-model expression is:
Figure FDA0003637731690000065
s.t.By≤d,(γ 1 ) (46)
Dy=e,(γ 2 ) (47)
Fy≤f-Ex * ,(γ 3 ) (48)
Gy≤p,(γ 4 ) (49)
Hy=a,(γ 5 ) (50)
Iy=r,(γ 6 ) (51)
in the formula: x is the number of * For the optimal solution of the main model, { gamma { 123456 Is a constraint dual variable;
the submodel equivalent representation method is as follows:
worst scenario p in sub-model * 、a * And r * To not determine the extreme values in sets P, A and R, equations (1) - (3) are therefore equivalent to:
Figure FDA0003637731690000066
Figure FDA0003637731690000067
Figure FDA0003637731690000068
in the formula: g is a radical of formula t 、m t
Figure FDA0003637731690000069
And
Figure FDA00036377316900000610
are all binary variables;
based on strong dual theory, submodels (45) - (51) are equivalent to:
Figure FDA0003637731690000071
-B T γ 1 +D T γ 2 -F T γ 3 -G T γ 4 +H T γ 5 +I T γ 6 =c T (56)
γ 1 ≥0,γ 3 ≥0,γ 4 ≥0,γ 256 is a free variable (57)
Figure FDA0003637731690000072
Figure FDA0003637731690000073
Figure FDA0003637731690000074
Figure FDA0003637731690000075
In the formula: lambda, mu, omega 1 And ω 2 Are all auxiliary variables, and are all the auxiliary variables,
Figure FDA0003637731690000076
CN202010154115.1A 2020-03-07 2020-03-07 Robust energy management method considering uncertain photovoltaic and load for traction power supply system Active CN111342450B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010154115.1A CN111342450B (en) 2020-03-07 2020-03-07 Robust energy management method considering uncertain photovoltaic and load for traction power supply system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010154115.1A CN111342450B (en) 2020-03-07 2020-03-07 Robust energy management method considering uncertain photovoltaic and load for traction power supply system

Publications (2)

Publication Number Publication Date
CN111342450A CN111342450A (en) 2020-06-26
CN111342450B true CN111342450B (en) 2022-08-23

Family

ID=71184310

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010154115.1A Active CN111342450B (en) 2020-03-07 2020-03-07 Robust energy management method considering uncertain photovoltaic and load for traction power supply system

Country Status (1)

Country Link
CN (1) CN111342450B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113422364B (en) * 2021-05-06 2023-04-28 华翔翔能科技股份有限公司 Full-buried variable load management method considering two electricity rates
CN114977182B (en) * 2022-04-25 2024-05-28 西南交通大学 Flexible traction power supply system interval optimal power flow optimization method considering photovoltaic access
CN117060491B (en) * 2023-10-11 2024-01-30 中国电建集团西北勘测设计研究院有限公司 Operation optimization method, system, medium and equipment of wind-solar hybrid energy storage system
CN117937487B (en) * 2024-03-22 2024-05-28 华东交通大学 Distributed energy management method for traction power supply system of heavy haul railway

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107887903A (en) * 2017-10-31 2018-04-06 深圳供电局有限公司 Consider the micro-capacitance sensor robust Optimization Scheduling of element frequency characteristic
CN107979111A (en) * 2017-07-21 2018-05-01 天津大学 A kind of energy management method for micro-grid based on the optimization of two benches robust
CN108108846A (en) * 2017-12-28 2018-06-01 东南大学 A kind of alternating current-direct current mixing microgrid robust optimizes coordinated scheduling method
CN108233415A (en) * 2018-01-15 2018-06-29 合肥工业大学 Two-stage type photovoltaic DC-to-AC converter virtual synchronous generator control method
CN108258695A (en) * 2018-02-28 2018-07-06 东南大学 A kind of random robust coupled mode Optimization Scheduling of alternating current-direct current series-parallel connection microgrid
CN108321854A (en) * 2018-03-13 2018-07-24 燕山大学 It is a kind of to consider that the micro-capacitance sensor event of communication time lag triggers hierarchical control method
CN109274134A (en) * 2018-11-08 2019-01-25 东南大学 A kind of active distribution network robust active reactive coordination optimizing method based on time series scene analysis
CN109659980A (en) * 2019-01-22 2019-04-19 西南交通大学 The tractive power supply system energy management optimization method of integrated hybrid energy-storing and photovoltaic devices
CN110390467A (en) * 2019-06-25 2019-10-29 河海大学 A kind of random ADAPTIVE ROBUST Optimization Scheduling of virtual plant distinguished based on key scenes
CN110474367A (en) * 2019-08-05 2019-11-19 广东工业大学 A kind of micro-capacitance sensor capacity configuration optimization method considering risk of loss
CN110829433A (en) * 2019-11-14 2020-02-21 青岛特锐德高压设备有限公司 Novel electric energy interaction system

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105140917B (en) * 2015-09-06 2017-05-10 清华大学 Robust restoration control method suitable for active power distribution network under uncertain environment

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107979111A (en) * 2017-07-21 2018-05-01 天津大学 A kind of energy management method for micro-grid based on the optimization of two benches robust
CN107887903A (en) * 2017-10-31 2018-04-06 深圳供电局有限公司 Consider the micro-capacitance sensor robust Optimization Scheduling of element frequency characteristic
CN108108846A (en) * 2017-12-28 2018-06-01 东南大学 A kind of alternating current-direct current mixing microgrid robust optimizes coordinated scheduling method
CN108233415A (en) * 2018-01-15 2018-06-29 合肥工业大学 Two-stage type photovoltaic DC-to-AC converter virtual synchronous generator control method
CN108258695A (en) * 2018-02-28 2018-07-06 东南大学 A kind of random robust coupled mode Optimization Scheduling of alternating current-direct current series-parallel connection microgrid
CN108321854A (en) * 2018-03-13 2018-07-24 燕山大学 It is a kind of to consider that the micro-capacitance sensor event of communication time lag triggers hierarchical control method
CN109274134A (en) * 2018-11-08 2019-01-25 东南大学 A kind of active distribution network robust active reactive coordination optimizing method based on time series scene analysis
CN109659980A (en) * 2019-01-22 2019-04-19 西南交通大学 The tractive power supply system energy management optimization method of integrated hybrid energy-storing and photovoltaic devices
CN110390467A (en) * 2019-06-25 2019-10-29 河海大学 A kind of random ADAPTIVE ROBUST Optimization Scheduling of virtual plant distinguished based on key scenes
CN110474367A (en) * 2019-08-05 2019-11-19 广东工业大学 A kind of micro-capacitance sensor capacity configuration optimization method considering risk of loss
CN110829433A (en) * 2019-11-14 2020-02-21 青岛特锐德高压设备有限公司 Novel electric energy interaction system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Energy Management of Connected Co-Phase Traction Power System Considering HESS and PV;Yuanli Liu et al;《2019 14th IEEE Conference on Industrial Electronics and Applications (ICIEA)》;20190916;1408-1412 *
考虑碳排放交易的日前调度双阶段鲁棒优化模型;张刚等;《中国电机工程学报》;20180920;第38卷(第18期);5490-5499 *

Also Published As

Publication number Publication date
CN111342450A (en) 2020-06-26

Similar Documents

Publication Publication Date Title
CN111342450B (en) Robust energy management method considering uncertain photovoltaic and load for traction power supply system
CN109659980B (en) Energy management optimization method for traction power supply system integrating hybrid energy storage and photovoltaic device
CN110739725B (en) Optimal scheduling method for power distribution network
CN111313465B (en) Energy management method for flexible traction power supply system containing photovoltaic and hybrid energy storage device
CN112734098B (en) Power distribution network power dispatching method and system based on source-load-network balance
CN110148969B (en) Active power distribution network optimized operation method based on model predictive control technology
CN110135662B (en) Energy storage site selection constant volume multi-objective optimization method considering reduction of peak-valley difference
CN109978240B (en) Electric automobile ordered charging optimization method and system
CN114977320A (en) Power distribution network source-network charge-storage multi-target collaborative planning method
CN107732937B (en) Peak clipping and valley filling method for grid-connected micro-grid containing wind-light-storage-electric automobile
CN110994697B (en) Optimal operation control method and system for alternating current-direct current distribution network containing light storage combination
Chen et al. Dynamic voltage unbalance constrained economic dispatch for electrified railways integrated energy storage
CN115102160A (en) Day-ahead energy optimization scheduling method for traction power supply system under weak power grid condition
Zhang et al. Multi-objective optimal scheduling strategy of microgrid based on v2g technology
Zhou et al. Research review on energy storage technology in power grid
CN114977182A (en) Photovoltaic access-considered optimal power flow optimization method for flexible traction power supply system interval
Chang et al. A coordinated charging strategy for PV-Assisted charging station of electric vehicles based on charging service price
Hafiz et al. Network constraints consideration for grid-edge energy management system
CN113690925A (en) Energy interaction optimization method and system based on micro-grid
Xu et al. Smart grid dispatching strategy considering the difference of electric vehicle demand
Zhang et al. Research on optimal configuration of energy storage in wind-solar microgrid considering real-time electricity price
Thakre et al. Potentially affect of a vehicle to grid on the electricity system
Cuiqing et al. A hybrid ac/dc microgrid energy management strategy based on neural network
Xiao et al. Optimization of Distribution Network Considering Distributed Renewable Energy and V2G of Electric Vehicles
Wu et al. A study on photovoltaic charging strategy of electric vehicles with multi-objective constraints

Legal Events

Date Code Title Description
PB01 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