CN111106612B - Energy storage type charging pile participating power grid demand side response combined operation optimization model and solving algorithm - Google Patents

Energy storage type charging pile participating power grid demand side response combined operation optimization model and solving algorithm Download PDF

Info

Publication number
CN111106612B
CN111106612B CN201811239625.8A CN201811239625A CN111106612B CN 111106612 B CN111106612 B CN 111106612B CN 201811239625 A CN201811239625 A CN 201811239625A CN 111106612 B CN111106612 B CN 111106612B
Authority
CN
China
Prior art keywords
energy storage
charging pile
storage type
type charging
response
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
CN201811239625.8A
Other languages
Chinese (zh)
Other versions
CN111106612A (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.)
Guizhou Power Grid Co Ltd
Original Assignee
Guizhou Power Grid Co Ltd
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 Guizhou Power Grid Co Ltd filed Critical Guizhou Power Grid Co Ltd
Priority to CN201811239625.8A priority Critical patent/CN111106612B/en
Publication of CN111106612A publication Critical patent/CN111106612A/en
Application granted granted Critical
Publication of CN111106612B publication Critical patent/CN111106612B/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/28Arrangements for balancing of the load in a network by storage of energy

Landscapes

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

Abstract

The invention discloses a method for optimizing participation of an energy storage type charging pile in power grid demand side response joint operation, which is characterized by comprising the steps of initializing, establishing an energy storage type charging pile participation power grid demand side response joint operation optimization model, solving distributed energy joint optimization based on original-dual decomposition and the like. And establishing an energy storage type charging pile participation power grid demand side response combined operation optimization model which comprises an objective function of a maximized benefit function and constraint conditions such as energy storage type charging pile charging and discharging capacity constraint, capacity constraint and demand response electric quantity constraint. The distributed energy joint optimization solution based on the original-dual decomposition comprises the steps of constructing a Lagrangian function, constructing a Lagrangian dual function, applying a KKT condition, updating an iterative operator based on the original-dual decomposition and the like. The method provided by the invention can guide an energy storage type charging pile operator to adjust the power consumption according to the real-time power price change signal, improve the resource efficiency of the power, realize the maximization of social welfare, improve the economy of the energy storage type charging pile and fully exert the advantage of the energy storage type charging pile in the power grid for participating in the response of the demand side.

Description

Energy storage type charging pile participating power grid demand side response combined operation optimization model and solving algorithm
Technical Field
The invention belongs to the field of operation optimization of charging facilities, and relates to a method for optimizing the operation of an energy storage type charging pile participating in the response joint of the demand side of a power grid by combining with real-time electricity price incentive and comprehensively considering the influence of the performance of an energy storage system on the management compensation of the demand side.
Background
In recent years, electric automobiles in China are rapidly developed. The development of the electric automobile needs to invest in the construction of matched charging infrastructures including charging piles, so that the construction of the charging infrastructures of the electric automobile is greatly promoted, and the development of the electric automobile is an urgent task for accelerating the popularization and the application of the electric automobile. The mobile/distributed energy storage system is installed in the charging pile and the nearby charging pile in a matched mode, impact on a power grid caused by a random charging mode of an electric automobile can be relieved, and the energy storage type charging pile system plays an important role in future charging pile design.
To energy storage system and fill electric pile complex optimization operation problem, the existing research mainly focuses on how to utilize energy storage system to reduce electric automobile charging fluctuation. However, as a large number of electric vehicles are connected to the power grid under the condition of no adjustment and control, peak-to-peak superposition of a load curve of the power distribution network may occur, so that the peak-to-valley difference is further enlarged, the load risk of a line transformer is increased, and the effect and the economy of the research on the charging harmonic suppression on the optimal operation of the energy storage type charging pile and the power grid are not obvious enough.
At present, the cost of an energy storage system is high, and the investment cost of energy storage is difficult to recover in a mode of singly carrying out price arbitrage on a charging and discharging integrated energy storage charging pile. And energy storage formula fills electric pile and participates in power grid demand side response service and both can exert energy storage system power rapid adjustment's advantage, can obtain considerable income again, simultaneously, also can reduce the economic cost to energy storage system through the management of ordered charge-discharge.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: aiming at the existing problems, the method for optimizing the response joint operation of the energy storage type charging pile participating in the power grid demand side is provided by combining the real-time electricity price incentive and comprehensively considering the influence of the performance of an energy storage system on the management compensation of the demand side.
The invention discloses a method for optimizing the participation of an energy storage type charging pile in power grid demand side response combined operation, which is characterized by comprising the following steps of:
step (1): initialization
Step (2): establishing energy storage type charging pile participation power grid demand side response combined operation optimization model
Step (2-1): the objective function is determined as follows:
Figure BSA0000172610710000021
represents: the energy storage type charging pile participates in a power grid demand side response joint operation optimization model, and a maximum benefit function is taken as a target.
Step (2-2): determining constraints
And (3): solving energy storage type charging pile participation power grid demand side response joint operation optimization model by adopting distributed energy joint optimization solving algorithm based on original-dual decomposition
Still further, the step (1) includes:
a) leading in energy storage formula fills electric pile parameter: and importing parameters such as energy storage capacity, maximum charge and discharge power, initial charge state and the like.
b) Importing power grid demand side response excitation power price information: and importing relevant excitation factors and electricity price information of the power grid for response excitation at the demand side.
Furthermore, in the step (2-1), in the objective function, the energy storage type charging pile benefit function is composed of weighted corresponding profit of the demand side and economic consumption of the energy storage type charging pile. The energy storage type charging pile operator determines power dispatching of the energy storage type charging pile participating in power grid demand response according to the change of the power price of the power grid operator demand response
Figure BSA0000172610710000031
Charging and discharging scheduling of energy storage type charging pile
Figure BSA0000172610710000032
As shown in the following formula:
Figure BSA0000172610710000033
wherein the content of the first and second substances,
Figure BSA0000172610710000034
electric quantity for representing participation of energy storage type charging pile i in power grid demand response in t time slot
Figure BSA0000172610710000035
The gain in the amount of gain that is obtained,
Figure BSA0000172610710000036
and indicating the electric quantity of the energy storage type charging pile i participating in the response of the power grid demand side in the t time slot.
Figure BSA0000172610710000037
Representing cost-effectiveness of energy storage systems, mainly by charging/discharging current
Figure BSA0000172610710000038
And (6) determining.
Under V2G technical condition, electric current can flow between energy storage formula fills electric pile and electric wire netting bidirectionally, and energy storage formula fills electric pile and can purchase the electric quantity from the electric wire netting and be used for real-time electric automobile demand and energy storage system charge-discharge demand, and the electric quantity that electric pile participated in the electric wire netting demand response is filled to the energy storage formula promptly
Figure BSA0000172610710000039
Real-time charging demand by electric vehicle
Figure BSA00001726107100000310
Charging and discharging scheduling of energy storage type charging pile
Figure BSA00001726107100000311
Consists of the following components:
Figure BSA00001726107100000312
furthermore, the two contents included in the energy storage charging pile benefit function are specifically as follows:
a) energy storage type charging pile participation demand side response benefit
When the energy storage type charging pile participates in the power grid demand side response combined operation, the energy storage type charging pile is adjusted to participate in the power grid demand electric quantity according to the change of the power grid side demand price
Figure BSA00001726107100000313
And the demand response benefits obtained. According to the response benefit of the energy storage system on the participation demand side, storing energyThe benefit that electric pile participated in the power grid demand side response is filled to the formula generally includes: incentive rewards based on participation demand response power and incentive rewards based on system benefits. The grid operator gives P the two demand response stimulitAnd mutThe unit electric quantity reward, the demand response income that fills electric pile operator and obtain promptly is as follows:
Figure BSA00001726107100000314
in the formula (I), the compound is shown in the specification,
Figure BSA00001726107100000315
show and fill electric pile average electric quantity demand in T time, show as:
Figure BSA00001726107100000316
b) energy storage economical loss of energy storage type charging pile
Considering the dynamic change of the State of charge (Soc) caused by the charging and discharging operations of the energy storage charging pile in the operation process, the Soc dynamic model can be expressed as:
Figure BSA0000172610710000041
in the formula (I), the compound is shown in the specification,
Figure BSA0000172610710000048
represents the Soc state of the charging post i at time t of the battery,
Figure BSA0000172610710000042
the initial Soc state of the battery is shown, the unit of the initial Soc state is battery energy storage capacity (Ah), and T represents the time period of the combined operation of the energy storage system and the charging pile. Different charging/discharging operations of the energy storage system can generate economic expenses of different degrees for the energy storage type charging pile, and the expenses are mainly generated by two charging/discharging behavior indexes, namely, charging/based on whether the charging/based on the charging status is on the basis of the charging status of the energy storage system or notDepth of discharge and charge/discharge frequency. The energy storage economy loss can be expressed as:
Figure BSA0000172610710000043
in the formula, deltai,αiAnd gammaiAre all that the normal number weight value represents the influence of different charging/discharging behaviors on the economic loss of the battery,
Figure BSA0000172610710000044
representing the energy storage system maximum capacity.
Still further, the constraint conditions of the step (2-2) include:
a) energy storage type charging pile charge and discharge amount constraint
The constraint condition of the charging and discharging amount of the energy storage type charging pile in any unit time slot represents the electric quantity limit constraint in the charging/discharging process of the charging pile, and is shown as the following formula:
Figure BSA0000172610710000045
wherein, the first and the second end of the pipe are connected with each other,
Figure BSA0000172610710000046
the maximum charging/discharging power of the charging pile. The formula represents the charge/discharge amount of any energy storage type charging pile in any time slot
Figure BSA0000172610710000047
The maximum charge/discharge capacity limit is satisfied and is limited within the rated charge/discharge capacity range.
b) Energy storage type charging pile capacity constraint
Because the energy storage formula fills electric pile has the restriction of certain capacity, energy storage system Soc value is too high or low will influence and fill the normal operation requirement of electric pile, for example, energy storage system Soc deviates from the median seriously, may lead to the energy storage because the electric quantity is too high or too low unable response fills electric pile system's instruction, and Soc is a continuous dynamic model, and its constraint condition can be expressed as:
Figure BSA0000172610710000051
according to the above formula, the capacity constraint of the energy storage system will charge/discharge the charging pile Ci tCreating an indirect effect.
c) Demand response power constraints
The electric quantity that electric pile participated in the power grid demand response should satisfy the maximum electric quantity restriction at the energy storage formula in arbitrary unit time slot, and its response electric quantity restraint is:
Figure BSA0000172610710000052
in the formula (I), the compound is shown in the specification,
Figure BSA0000172610710000053
and the maximum electric quantity limit of the energy storage type charging pile i participating in demand response in the unit time slot t is shown.
Furthermore, in the step (3), a solution algorithm of distributed energy joint optimization based on the original-dual decomposition is adopted to solve the optimization model obtained in the step (2). Through obtaining the operation function of each charging pile energy management, the KKT condition is applied to realize the decomposition of the problem, so that each charging pile can be iteratively updated according to the operation state of the charging pile, and finally, the optimal solution is obtained. The solving algorithm is divided into the following four steps:
a) construction of Lagrangian functions
Constructing a Lagrangian function based on the optimization objective function as shown in the following formula:
Figure BSA0000172610710000054
where λ, v are lagrange multipliers greater than 0.
b) Constructing a Lagrangian dual function:
on the basis of the Lagrangian function constructed in the previous step, taking the dual form of the Lagrangian function to construct and obtain the Lagrangian dual function, wherein the Lagrangian dual function is shown as the following formula:
Figure BSA0000172610710000061
according to the dual theorem, g (lambda, v) is less than or equal to P on the feasible fields of lambda and v*And P is the optimal solution of the constraint of the original problem in the constraint limiting condition.
c) Applying KKT conditions
When the feasible domain of solving the original problem is satisfied, the feasible domain of the dual problem is as follows:
λ>0,v>0
the complementary relaxation amounts are:
Figure BSA0000172610710000062
Figure BSA0000172610710000063
because the optimization problem model is a convex function, the feasible solution is a convex set, and the solution obtained by applying the KKT condition is the solution of the original problem, which is a sufficient necessary condition. Therefore, the problem decomposition method can be applied to obtain the operator of each charging pile needing to be solved and updated.
d) Updating operators based on primitive-dual decomposition
According to the principle of original-dual decomposition, each energy storage charging pile is according to P collectedt、μtAfter the signal, according to self restrictive condition distributed update energy storage formula fill electric pile operation dispatch, the update operator as follows:
Figure BSA0000172610710000064
Figure BSA0000172610710000065
Figure BSA0000172610710000066
in the formula (I), the compound is shown in the specification,
Figure BSA0000172610710000067
representing the gradient of the function L.
No parameter is coupled between the charging piles in the solving process, so that the splitting and solving can be carried out. The complexity of the algorithm mainly depends on the scale of an applicable scene, is related to the total number N of the charging piles participating in the joint optimization operation, is also related to the iteration number k in the operation process of each algorithm, and is referred to as o (kN).
In conclusion, due to the adoption of the technical scheme, the energy storage type charging pile has the advantages that an energy storage type charging pile operator can be guided to adjust the power consumption according to the real-time power price change signal, the resource efficiency of the power is improved, the social welfare maximization is realized, the economy of the energy storage type charging pile is improved, the advantage of the energy storage type charging pile in participating in the response of a demand side in a power grid is fully exerted, and the healthy development of electric vehicles and energy storage industries is facilitated.
Drawings
Fig. 1 is a flowchart of a method for optimizing the joint operation of energy storage charging piles participating in the response of the demand side of a power grid.
Detailed Description
All of the features disclosed in this specification, or all of the steps of any method or process so disclosed, may be combined in any combination, except combinations of features and/or steps that are mutually exclusive.
Any feature disclosed in this specification may be replaced by alternative features serving an equivalent or similar purpose, unless expressly stated otherwise. That is, unless expressly stated otherwise, each feature is only an example of a generic series of equivalent or similar features.
The energy storage type charging pile participating power grid demand side response combined operation optimization method provided by the invention is described in more detail below with reference to the accompanying drawings.
The invention discloses a method for optimizing the participation of an energy storage type charging pile in power grid demand side response combined operation, which is characterized by comprising the following steps of:
step (1): initialization
Step (2): establishing energy storage type charging pile participation power grid demand side response combined operation optimization model
Step (2-1): the objective function is determined as follows:
Figure BSA0000172610710000081
represents: the energy storage type charging pile participates in a power grid demand side response joint operation optimization model, and a maximum benefit function is taken as a target.
Step (2-2): determining constraints
And (3): and solving a power storage type charging pile participation power grid demand side response joint operation optimization model by adopting a distributed energy joint optimization solving algorithm based on original-dual decomposition.
Still further, the step (1) includes:
a) importing energy storage type charging pile parameters: and importing parameters such as energy storage capacity, maximum charge and discharge power, initial charge state and the like.
b) Importing power grid demand side response excitation power price information: and importing relevant excitation factors and electricity price information of the power grid for response excitation at the demand side.
Furthermore, in the step (2-1), in the objective function, the energy storage type charging pile benefit function is composed of weighted demand side corresponding income and energy storage type charging pile economic consumption. The energy storage type charging pile operator determines power dispatching of participation of the energy storage type charging pile in power grid demand response according to change of demand response power price of the power grid operator
Figure BSA0000172610710000082
Charging and discharging scheduling of energy storage type charging pile
Figure BSA0000172610710000083
As shown in the following formula:
Figure BSA0000172610710000084
wherein the content of the first and second substances,
Figure BSA0000172610710000085
electric quantity for representing participation of energy storage type charging pile i in power grid demand response in t time slot
Figure BSA0000172610710000086
The gain in the amount of gain that is obtained,
Figure BSA0000172610710000087
and indicating the electric quantity of the energy storage type charging pile i participating in response of the power grid demand side in the t time slot.
Figure BSA0000172610710000088
Representing cost-effectiveness of energy storage systems, mainly by charging/discharging current
Figure BSA0000172610710000089
And (6) determining.
Under V2G technical condition, electric current can flow between energy storage formula fills electric pile and electric wire netting bidirectionally, and energy storage formula fills electric pile and can purchase the electric quantity from the electric wire netting and be used for real-time electric automobile demand and energy storage system charge-discharge demand, and the electric quantity that electric pile participated in the electric wire netting demand response is filled to the energy storage formula promptly
Figure BSA00001726107100000810
Real-time charging demand by electric vehicle
Figure BSA00001726107100000811
Charging and discharging scheduling of energy storage type charging pile
Figure BSA00001726107100000812
Consists of the following components:
Figure BSA0000172610710000091
furthermore, the two contents included in the energy storage type charging pile benefit function are specifically as follows:
a) energy storage type charging pile participation demand side response benefit
When the energy storage type charging pile participates in power grid demand side response combined operation, the energy storage type charging pile is adjusted to participate in power grid demand electric quantity according to the change of power grid side demand price
Figure BSA0000172610710000092
And the demand response benefits obtained. According to the benefit that energy storage system responded in participating in demand side, the benefit that energy storage formula fills electric pile and participates in electric wire netting demand side response generally includes: incentive rewards based on participation demand response power and incentive rewards based on system benefits. The grid operator gives P the two demand response stimulitAnd mutThe unit electric quantity of reward, the demand response income that fills electric pile operator and obtain promptly is as follows:
Figure BSA0000172610710000093
in the formula (I), the compound is shown in the specification,
Figure BSA0000172610710000094
show and fill electric pile average electric quantity demand in T time, show as:
Figure BSA0000172610710000095
b) energy storage economical loss of energy storage type charging pile
Considering the dynamic change of the State of charge (Soc) caused by the charging and discharging operations of the energy storage charging pile in the operation process, the Soc dynamic model can be expressed as:
Figure BSA0000172610710000096
in the formula (I), the compound is shown in the specification,
Figure BSA0000172610710000097
represents the Soc state of the charging post i at time t of the battery,
Figure BSA0000172610710000098
the initial Soc state of the battery is shown, the unit of the initial Soc state is battery energy storage capacity (Ah), and T represents the time period of the combined operation of the energy storage system and the charging pile. Different charging/discharging operations of the energy storage system can generate economic expenses of different degrees for the energy storage type charging pile, and the expenses are mainly generated by two charging/discharging behavior indexes, namely charging/discharging depth and charging/discharging frequency. The energy storage economy loss can be expressed as:
Figure BSA0000172610710000099
in the formula, deltai,αiAnd gammaiAre all that the normal number weight value represents the influence of different charging/discharging behaviors on the economic loss of the battery,
Figure BSA0000172610710000101
representing the energy storage system maximum capacity.
Further, the constraints of the step (2-2) include:
a) energy storage type charging pile charge and discharge amount constraint
The constraint condition of the charging and discharging amount of the energy storage type charging pile in any unit time slot represents the electric quantity limit constraint in the charging/discharging process of the charging pile, and is shown as the following formula:
Figure BSA0000172610710000102
wherein the content of the first and second substances,
Figure BSA0000172610710000103
the maximum charging/discharging power of the charging pile. The above formula shows the charge/discharge amount of any energy storage type charging pile in any time slot
Figure BSA0000172610710000104
The maximum charge/discharge capacity limit is satisfied and is limited within the rated charge/discharge capacity range.
b) Energy storage type charging pile capacity constraint
Because the energy storage formula fills electric pile has the restriction of certain capacity, energy storage system Soc value is too high or low will influence and fill the normal operation requirement of electric pile, for example, energy storage system Soc deviates from the median seriously, may lead to the energy storage because the electric quantity is too high or too low unable response fills electric pile system's instruction, and Soc is a continuous dynamic model, and its constraint condition can be expressed as:
Figure BSA0000172610710000105
according to the above formula, the capacity constraint of the energy storage system will charge/discharge the charging pile Ci tCreating an indirect effect.
c) Demand response power constraints
The electric quantity that electric pile participated in the power grid demand response should satisfy the maximum electric quantity restriction at the energy storage formula in arbitrary unit time slot, and its response electric quantity restraint is:
Figure BSA0000172610710000106
in the formula (I), the compound is shown in the specification,
Figure BSA0000172610710000107
and the maximum electric quantity limit of the energy storage type charging pile i participating in demand response in the unit time slot t is shown.
Furthermore, in the step (3), a solution algorithm of distributed energy joint optimization based on primitive-dual decomposition is adopted to solve the optimization model obtained in the step (2). Through obtaining the operation function of each charging pile energy management, the KKT condition is applied to realize the decomposition of the problem, so that each charging pile can be iteratively updated according to the operation state of the charging pile, and finally, the optimal solution is obtained. The solving algorithm is divided into the following four steps:
a) construction of Lagrangian functions
Constructing a Lagrangian function based on the optimization objective function as shown in the following formula:
Figure BSA0000172610710000111
where λ, v are lagrange multipliers greater than 0.
b) Constructing a Lagrangian dual function:
on the basis of the Lagrangian function constructed in the previous step, taking the dual form of the Lagrangian function to construct and obtain the Lagrangian dual function, wherein the Lagrangian dual function is shown as the following formula:
Figure BSA0000172610710000112
according to the dual theorem, g (lambda, v) is less than or equal to P on the feasible fields of lambda and v*P is the optimal solution of the original problem constrained within the constraint constraints.
c) Applying KKT conditions
When the feasible domain of solving the original problem is satisfied, the feasible domain of the dual problem is as follows:
λ>0,v>0
the complementary relaxation amounts are:
Figure BSA0000172610710000113
Figure BSA0000172610710000114
because the optimization problem model is a convex function, the feasible solution is a convex set, and the solution obtained by applying the KKT condition is the solution of the original problem, which is a sufficient necessary condition. Therefore, the problem decomposition method can be applied to obtain the operator of each charging pile needing to be solved and updated.
d) Updating operators based on primitive-dual decomposition
According to the principle of original-dual decomposition, each energy storage charging pile is according to P collectedt、μtAfter the signal, the operation scheduling of the energy storage type charging pile is updated in a distributed mode according to the self limiting condition, and the updating operator is as follows:
Figure BSA0000172610710000121
Figure BSA0000172610710000122
Figure BSA0000172610710000123
in the formula (I), the compound is shown in the specification,
Figure BSA0000172610710000124
representing the gradient of the function L.
No parameter is coupled between the charging piles in the solving process, so that the splitting and solving can be carried out. The complexity of the algorithm mainly depends on the scale of an applicable scene, is related to the total number N of the charging piles participating in the joint optimization operation, is also related to the iteration number k in the operation process of each algorithm, and is referred to as o (kN).
The invention is not limited to the foregoing embodiments. The invention extends to any novel feature or any novel combination of features disclosed in this specification, and to any novel method or process steps or any novel combination of steps disclosed.

Claims (4)

1. A method for optimizing the joint operation of energy storage type charging piles participating in power grid demand side response is characterized by comprising the following steps:
step (1): initialization
The step (1) comprises the following steps:
a) importing parameters of the energy storage type charging pile, including importing energy storage capacity, maximum charging and discharging power and an initial charge state;
b) importing power grid demand side response excitation power price information: importing relevant excitation factors and electricity price information of a power grid for response excitation at a demand side;
step (2): establishing energy storage type charging pile participation power grid demand side response combined operation optimization model
Step (2-1): the objective function is determined as follows:
Figure FDA0003638497660000011
represents: the energy storage type charging pile participates in a power grid demand side response joint operation optimization model, and a maximum benefit function is taken as a target; qi tAnd C, representing the electric quantity of the energy storage type charging pile i participating in response of the power grid demand side in the t time sloti tThe charging/discharging amount of any energy storage type charging pile in any time slot is represented;
step (2-2): determining constraints
And (3): and solving a power storage type charging pile participation power grid demand side response joint operation optimization model by adopting a distributed energy joint optimization solving algorithm based on original-dual decomposition.
2. The method for optimizing the response joint operation of the energy storage type charging pile participating in the power grid on the demand side according to claim 1, characterized by comprising the following steps: the algorithm is used for realizing the decomposition of the problem by solving the operation function of energy management of each charging pile and then applying the KKT condition, so that each charging pile can be iteratively updated according to the operation state of the charging pile, and finally the optimal solution is obtained, and the algorithm comprises the following steps:
a) construction of Lagrangian functions
Constructing a Lagrangian function based on the optimization objective function as shown in the following formula:
Figure FDA0003638497660000021
wherein λ, v are lagrange multipliers greater than 0;
b) constructing a Lagrangian dual function:
on the basis of the Lagrangian function constructed in the previous step, taking the dual form of the Lagrangian function to construct and obtain the Lagrangian dual function, wherein the Lagrangian dual function is shown as the following formula:
Figure FDA0003638497660000022
according to the dual theorem, g (lambda, v) is less than or equal to P on the feasible fields of lambda and v*P is the optimal solution of the original problem constrained within constraint limiting conditions;
c) applying KKT conditions
When the feasible domain of solving the original problem is satisfied, the feasible domain of the dual problem is as follows:
λ>0,v>0
the complementary relaxation amounts are:
Figure FDA0003638497660000023
Figure FDA0003638497660000024
because the optimization problem model is a convex function, the feasible solution is a convex set, the solution obtained by applying the KKT condition is the solution of the original problem and is a sufficient necessary condition, the operator required to be updated by each charging pile can be obtained by applying a problem decomposition method;
d) updating operators based on primitive-dual decomposition
According to the principle of original-dual decomposition, each energy storage charging pile is according to P collectedt、μtAfter the signal, the energy storage type charging is updated in a distributed way according to the self limiting conditionAnd (3) electric pile operation scheduling, wherein an update operator is as follows:
Figure FDA0003638497660000031
Figure FDA0003638497660000032
Figure FDA0003638497660000033
in the formula (I), the compound is shown in the specification,
Figure FDA0003638497660000034
represents the gradient of the function L;
in the solving process, parameters are not coupled among the charging piles, so that the splitting solving can be carried out; the complexity of the algorithm mainly depends on the scale of an applicable scene, is related to the total number N of the charging piles participating in the joint optimization operation, is also related to the iteration times k in each algorithm operation process, and is referred to as o (kN);
Coi t(Ci t) Represents an economic cost of the energy storage system; qi tRepresenting the electric quantity of the energy storage type charging pile i participating in response of the power grid demand side in a t time slot; ci tThe charging/discharging amount of any energy storage type charging pile in any time slot is represented; drt t(Qi t) Representing energy storage type charging pile i participating in power grid demand response electric quantity Q in t time sloti tThe gain obtained; ci t,maxThe maximum charging/discharging power of the charging pile; qi t,maxRepresenting the maximum electric quantity limit of the energy storage type charging pile i participating in demand response in the unit time slot t; and T represents the time period of the combined operation of the energy storage system and the charging pile.
3. The method of claim 1The energy storage type charging pile participation power grid demand side response combined operation optimization method is characterized by comprising the following steps: in the step (2-1), in the objective function, the energy storage type charging pile benefit function is composed of weighted demand side corresponding income and energy storage type charging pile economic consumption; the energy storage type charging pile operator determines the power dispatching Q of the energy storage type charging pile participating in power grid demand response according to the change of the demand response power price of the power grid operatori tCharging and discharging scheduling C for energy storage type charging pilei tAs shown in the following formula:
Figure FDA0003638497660000041
wherein, Drt t(Qi t) Representing energy storage type charging pile i participating in power grid demand response electric quantity Q in t time sloti tGain obtained, Qi tRepresenting the electric quantity of the energy storage type charging pile i participating in response of the power grid demand side in a t time slot; coi t(Ci t) Shows the economic cost of the energy storage system, and mainly comprises the charge/discharge capacity C of any energy storage type charging pile in any time sloti tDetermining;
under V2G technical condition, electric current can flow between energy storage formula fills electric pile and electric wire netting bidirectionally, and energy storage formula fills electric pile and can purchase the electric quantity from the electric wire netting and be used for real-time electric automobile demand and energy storage system charge-discharge demand, and the electric quantity Q that electric pile participated in the electric wire netting demand response is filled to the energy storage formula promptlyi tReal-time charging demand D of electric vehiclei tCharging and discharging scheduling C for energy storage type charging pilei tConsists of the following components:
Figure FDA0003638497660000042
furthermore, the two contents included in the energy storage charging pile benefit function are specifically as follows:
a) energy storage type charging pile participation demand side response benefit
When the energy storage type charging pile participates in the power grid demand side response combined operation, the energy storage type charging pile is adjusted to participate in the power grid demand electric quantity Q according to the change of the power grid side demand pricei tAnd the obtained demand response benefit; according to the benefit that energy storage system responded in participating in demand side, the benefit that energy storage formula fills electric pile and participates in electric wire netting demand side response generally includes: incentive rewards based on participation demand response power and incentive rewards based on system benefits; the electric network operator gives P two demand response stimulitAnd mutThe unit electric quantity reward, the demand response income that fills electric pile operator and obtain promptly is as follows:
Figure FDA0003638497660000051
in the formula (I), the compound is shown in the specification,
Figure FDA0003638497660000052
show and fill electric pile average electric quantity demand in the T time, show and be:
Figure FDA0003638497660000053
b) energy storage economical loss of energy storage type charging pile
Considering the dynamic change of the State of charge (Soc) caused by the charging and discharging operations of the energy storage charging pile in the operation process, the Soc dynamic model can be expressed as:
Figure FDA0003638497660000054
in the formula, Soci tRepresents the Soc state of the charging pile i at time t of the batteryi t0The initial Soc state of the battery is represented, the unit of the initial Soc state is the energy storage capacity (Ah) of the battery, and T represents the time period of the combined operation of the energy storage system and the charging pile; due to energy storageDifferent charging/discharging operations of the system can generate economic expenses of different degrees for the energy storage type charging pile, and the expenses are mainly generated by two charging/discharging behavior indexes, namely charging/discharging depth and charging/discharging frequency; the energy storage economy loss can be expressed as:
Figure FDA0003638497660000055
in the formula, deltai,αiAnd gammaiAll represent the influence of different charging/discharging behaviors on the economic loss of the battery by a normal number weight value, Soci maxRepresenting the maximum capacity of the energy storage system.
4. The method for optimizing the joint operation of the energy storage type charging pile participating in the response of the power grid demand side according to claim 1, wherein the method comprises the following steps: the constraint conditions of the step (2-2) comprise:
a) energy storage type charging pile charge and discharge amount constraint
The constraint condition of the charging and discharging amount of the energy storage type charging pile in any unit time slot represents the electric quantity limit constraint in the charging/discharging process of the charging pile, and is shown as the following formula:
Figure FDA0003638497660000061
wherein, Ci t,maxThe maximum charging/discharging power of the charging pile; the formula represents the charge/discharge capacity C of any energy storage type charging pile in any time sloti tThe maximum charge/discharge capacity limit is met and limited in the rated charge/discharge capacity range;
b) energy storage type charging pile capacity constraint
Because the energy storage formula fills electric pile has the restriction of certain capacity, energy storage system Soc value is too high or low will influence and fill the normal operation requirement of electric pile, for example, energy storage system Soc deviates from the median seriously, may lead to the energy storage because the electric quantity is too high or too low unable response fills electric pile system's instruction, and Soc is a continuous dynamic model, and its constraint condition can be expressed as:
Figure FDA0003638497660000062
according to the above formula, the capacity constraint of the energy storage system will charge/discharge C to/from the charging pilei tGenerating indirect influence; soci maxRepresenting the maximum capacity of the energy storage system;
c) demand response power constraints
The electric quantity that electric pile participated in the power grid demand response should satisfy the maximum electric quantity restriction at the energy storage formula in arbitrary unit time slot, and its response electric quantity restraint is:
Figure FDA0003638497660000063
in the formula, Qi t,maxRepresenting the maximum electric quantity limit of the energy storage type charging pile i participating in demand response in the unit time slot t;
furthermore, in the step (3), a solution algorithm of distributed energy joint optimization based on the original-dual decomposition is adopted to solve the optimization model obtained in the step (2).
CN201811239625.8A 2018-10-25 2018-10-25 Energy storage type charging pile participating power grid demand side response combined operation optimization model and solving algorithm Active CN111106612B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811239625.8A CN111106612B (en) 2018-10-25 2018-10-25 Energy storage type charging pile participating power grid demand side response combined operation optimization model and solving algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811239625.8A CN111106612B (en) 2018-10-25 2018-10-25 Energy storage type charging pile participating power grid demand side response combined operation optimization model and solving algorithm

Publications (2)

Publication Number Publication Date
CN111106612A CN111106612A (en) 2020-05-05
CN111106612B true CN111106612B (en) 2022-07-05

Family

ID=70417572

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811239625.8A Active CN111106612B (en) 2018-10-25 2018-10-25 Energy storage type charging pile participating power grid demand side response combined operation optimization model and solving algorithm

Country Status (1)

Country Link
CN (1) CN111106612B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111724080A (en) * 2020-06-29 2020-09-29 南京工程学院 Mobile charging pile group scheduling method considering battery pack health state balance
CN111738611B (en) * 2020-06-29 2024-04-26 南京工程学院 Intelligent scheduling method for mobile charging pile group based on Sarsa algorithm
CN112636330B (en) * 2020-12-02 2023-01-10 上海交通大学 Comprehensive energy system peak regulation method based on user contribution behaviors

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106253315A (en) * 2016-08-02 2016-12-21 万马联合新能源投资有限公司 A kind of energy accumulation capacity configuration considering electric automobile charging station schedulability
CN106364344A (en) * 2016-09-19 2017-02-01 北京慧峰聚能科技有限公司 Energy management method of electric automobile charging demand side distributed system

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104966127B (en) * 2015-06-03 2018-06-01 东南大学 A kind of electric vehicle economic load dispatching method based on demand response
CN106026152B (en) * 2016-05-19 2017-06-06 合肥工业大学 A kind of electric automobile accesses the discharge and recharge dispatching method of micro-capacitance sensor
CN106058904B (en) * 2016-06-30 2022-04-29 全球能源互联网研究院 Capacity matching method for charging pile of energy storage type electric automobile
CN106299511B (en) * 2016-08-02 2019-02-26 万马联合新能源投资有限公司 Electric automobile charging station energy storage capacity optimization method
CN106779487A (en) * 2017-01-17 2017-05-31 国网江西省电力公司电力科学研究院 A kind of virtual plant dynamic economic dispatch model based on prim al- dual interior point m ethod
CN107689634B (en) * 2017-08-24 2021-02-26 国网新疆电力公司电力科学研究院 Demand side bidirectional interaction control method based on energy storage type charging pile

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106253315A (en) * 2016-08-02 2016-12-21 万马联合新能源投资有限公司 A kind of energy accumulation capacity configuration considering electric automobile charging station schedulability
CN106364344A (en) * 2016-09-19 2017-02-01 北京慧峰聚能科技有限公司 Energy management method of electric automobile charging demand side distributed system

Also Published As

Publication number Publication date
CN111106612A (en) 2020-05-05

Similar Documents

Publication Publication Date Title
Martinez-Laserna et al. Battery second life: Hype, hope or reality? A critical review of the state of the art
Ahmad et al. A cost-efficient energy management system for battery swapping station
Heydarian-Forushani et al. Robust scheduling of variable wind generation by coordination of bulk energy storages and demand response
Saber et al. One million plug-in electric vehicles on the road by 2015
White et al. Using vehicle-to-grid technology for frequency regulation and peak-load reduction
Amirhosseini et al. Scheduling charging of hybrid-electric vehicles according to supply and demand based on particle swarm optimization, imperialist competitive and teaching-learning algorithms
CN111106612B (en) Energy storage type charging pile participating power grid demand side response combined operation optimization model and solving algorithm
CN105046371A (en) Electric vehicle charge-discharge scheduling method based on demand side bidding
Debnath et al. Energy storage model with gridable vehicles for economic load dispatch in the smart grid
Alfaverh et al. Optimal vehicle-to-grid control for supplementary frequency regulation using deep reinforcement learning
CN107394802A (en) Distributed energy storage participates in the control method for coordinating of Automatic Generation Control
Panwar et al. Strategic Energy Management (SEM) in a micro grid with modern grid interactive electric vehicle
CN108062619B (en) Rail vehicle-ground integrated capacity configuration method and device
CN116001624A (en) Ordered charging method for one-pile multi-connected electric automobile based on deep reinforcement learning
Rahman et al. Comprehensive impact analysis of electric vehicle charging scheduling on load-duration curve
Bhandari et al. The profitability of vehicle to grid for system participants-A case study from the Electricity Reliability Council of Texas
Bayram et al. A probabilistic capacity planning methodology for plug-in electric vehicle charging lots with on-site energy storage systems
Galus Agent-based modeling and simulation of large scale electric mobility in power systems
CN113887032A (en) Electric automobile ordered charging and discharging control method based on Lagrange distributed algorithm
Rezania et al. Business models for the integration of electric vehicles into the Austrian energy system
Gupta et al. Voltage profile improvement using demand side management in distribution networks under frequency linked pricing regime
KR101581685B1 (en) Apparatus and method for charge and discharge scheduling in energy storage device
Anand et al. Optimal day ahead scheduling of distributed EVs in a smart distribution network
Foley et al. Electric vehicles and displaced gaseous emissions
Konda et al. Binary fireworks algorithm application for optimal schedule of electric vehicle reserve in traditional and restructured electricity markets

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20201123

Address after: 550002 No. 17 Binhe Road, Nanming District, Guizhou, Guiyang

Applicant after: GUIZHOU POWER GRID Corp.

Address before: 550002 No. 32, Nanming, Guiyang District, Guizhou, Jiefang Road

Applicant before: ELECTRIC POWER RESEARCH INSTITUTE, GUIZHOU POWER GRID Co.,Ltd.

SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant