CN113346554A - Distributed cooperative regulation and control method for power distribution network - Google Patents

Distributed cooperative regulation and control method for power distribution network Download PDF

Info

Publication number
CN113346554A
CN113346554A CN202110529223.7A CN202110529223A CN113346554A CN 113346554 A CN113346554 A CN 113346554A CN 202110529223 A CN202110529223 A CN 202110529223A CN 113346554 A CN113346554 A CN 113346554A
Authority
CN
China
Prior art keywords
distribution network
microgrid
power distribution
power
cost
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.)
Granted
Application number
CN202110529223.7A
Other languages
Chinese (zh)
Other versions
CN113346554B (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.)
Shenyang Institute of Engineering
Original Assignee
Shenyang Institute of Engineering
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 Shenyang Institute of Engineering filed Critical Shenyang Institute of Engineering
Priority to CN202110529223.7A priority Critical patent/CN113346554B/en
Publication of CN113346554A publication Critical patent/CN113346554A/en
Application granted granted Critical
Publication of CN113346554B publication Critical patent/CN113346554B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48—Controlling the sharing of the in-phase component
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Abstract

The invention relates to a distributed cooperative regulation and control method for a power distribution network, which comprises the following steps: issuing an autonomous space plan by the power distribution network; the power distribution network judges the operation interval in which the actual operation point of each micro-grid falls according to the autonomous space plan; the operation point falls into a limit area, the boundary crossing behavior is stopped, the relevant information is counted and a punishment mechanism is incorporated; the operation point falls into an autonomous operation area, the power distribution network does not do forced interference, and the behavior information is collected into a reward mechanism after being counted; and the operation point falls into the cooperative control region, and the boundaries of the cooperative control region are set from outside to inside by utilizing a Benders decomposition algorithm, so that the boundaries of the cooperative control region are determined, and the upper limit and the lower limit of each-stage operation reward region of the multi-microgrid are obtained. The method and the system can improve the operation efficiency and safety of the intelligent power distribution network and protect the privacy information of the microgrid main body.

Description

Distributed cooperative regulation and control method for power distribution network
Technical Field
The invention relates to a power distribution network regulation and control method, in particular to a power distribution network distributed cooperative regulation and control method.
Background
The micro-grid provides a flexible and effective way for distributed grid-connected power generation and consumption and utilization of renewable energy sources, and also provides a friendly interface for operation and regulation of the power distribution network. The micro-grid connection enables a traditional passive power distribution network facing load regulation to step forward to an active intelligent power distribution system facing multiple micro-grids. The multi-microgrid high-permeability development situation inevitably leads to great changes of future intelligent power distribution systems in aspects of regulation and control object type diversification, regulation and control space-time dimension enlargement, regulation and control target diversification and the like, and the scientific and technical problems of trend management and control on a physical level, power distribution market construction on an economic level and the like appear. The conventional passive regulation and control mechanism of the power distribution network gradually shows obvious inadaptability, and a new cooperative and autonomous regulation and control mechanism of the power distribution system needs to be researched urgently to solve the contradiction between the standardization and the high autonomy of power generation and utilization behaviors of the multiple micro-grids, promote cooperative win-win among the multiple micro-grids while meeting diversified service requirements such as multiple micro-grid privacy protection, market bidding and the like, and improve the flexibility and the initiative of participating in the operation of the power distribution system.
Centralized cooperative regulation and control depends on high-density communication and control over all elements in the power distribution network and the microgrid, namely high perceptibility. When uncertainty or privacy exists in the multi-microgrid operation information, the centralized cooperative regulation and control method is difficult to seek the optimal operation point, so that the operation control of the power distribution network is trapped in trouble.
Disclosure of Invention
The invention aims to provide a distributed cooperative regulation and control method for a small-capacity distributed microgrid considering the protection requirement of private information, so that the private information of a microgrid main body is protected while the operation efficiency and the safety of an intelligent power distribution network are improved.
The technical scheme of the invention is that a distributed cooperative regulation and control method comprises the following steps:
step 1, issuing an autonomous space plan by a power distribution network;
step 2, the power distribution network judges the operation interval in which the actual operation point of each micro-grid falls according to the autonomous space plan;
step 3, for the micro-grid with the operating point falling into the limit area, when a forcible limit measure is adopted to terminate the border crossing behavior, relevant information is counted and incorporated into a punishment mechanism; for the micro-grid with the operating point falling into the autonomous operating area, the power distribution network does not do forced interference, and the behavior information is collected and then brought into a reward mechanism; and for the micro-grid with the operating points falling into the cooperative control area, setting the boundary of the cooperative control area from outside to inside by using a Benders decomposition algorithm, so as to determine the boundary of the cooperative control area and obtain the upper limit and the lower limit of each-level operation reward area of the multi-micro-grid.
In the step 3, determining the boundary of the cooperative control region comprises the following steps:
step a),
Figure BDA0003066598300000021
Figure BDA0003066598300000022
Iterate to tBNext total operating cost ($);
Cimpthe total cost of power input from the grid ($);
Cshedthe total load of the power distribution network reduces the cost ($);
Cvolthe cumulative voltage deviation cost ($) of the distribution network;
zmthe microgrid running cost ($) at the position m of the power distribution network bus;
ΩMa node set connected with the microgrid in the power distribution network;
step b), iterating to tBNext time, the feasibility checking sub-problem of the microgrid at the bus m is expressed as:
Figure BDA0003066598300000023
Figure BDA0003066598300000024
Figure BDA0003066598300000025
Figure BDA0003066598300000026
actual active power exchange capacity (kW) determined in the main question;
Figure BDA0003066598300000027
the actual amount of reactive power exchanged (kVar) determined in the main problem;
sP1,madditional non-negative relaxation variables
sP2,mAn additional non-negative slack variable;
sQ1,man additional non-negative slack variable;
sQ2,man additional non-negative slack variable;
Figure BDA0003066598300000031
minimum sum of relaxation variables of the microgrid;
Figure BDA0003066598300000032
the bivariate associated with equation (1-2);
Figure BDA0003066598300000033
the bivariate associated with equations (1-3);
step c), if the power exchange is feasible, solving an optimality check sub-problem by the multi-microgrid, as follows:
Figure BDA0003066598300000034
Figure BDA0003066598300000035
Figure BDA0003066598300000036
Figure BDA0003066598300000037
determining the minimum operation cost ($) of the microgrid under the power exchange schedule;
Figure BDA0003066598300000038
the sensitivity coefficients associated with equations (1-6), i.e., the Lagrangian multipliers;
Figure BDA0003066598300000039
the sensitivity coefficients associated with equations (1-7), i.e., the Lagrangian multipliers;
step d), after solving, utilizing the optimized target value
Figure BDA00030665983000000310
And the following valid inequality expressed as the optimality cut:
Figure BDA00030665983000000311
according to the weak duality of the convex optimization problem, an approximate value of the microgrid operation cost in the iteration process can be obtained through the formula;
step e), if the upper limit and the lower limit of the power exchange between the power distribution network and all the micro-grids are feasible, the upper limit of the operating cost of the power distribution network is expressed as follows:
Figure BDA00030665983000000312
Figure BDA00030665983000000313
at iteration tBAn upper bound ($) determined next time;
Figure BDA00030665983000000314
at iteration t, at the optimum value ofBDetermining ($) in a secondary main problem;
Figure BDA00030665983000000315
at iteration t, at the optimum value ofBDetermining ($) in a secondary main problem;
Figure BDA00030665983000000316
at iteration t, at the optimum value ofBDetermining ($) in a secondary main problem;
if the power exchange scheme is not feasible, its upper limit is set to infinity.
In said step b), if the determined power exchange is feasible, all slack variables will take zero, so that the target value in equation (1-1) is zero, otherwise, the non-zero target value in equation (1-1) indicates that the resulting power exchange scheme is not feasible, which will correct the main problem in the form of a feasibility cut, specifically, equation (1-1) will be included in the main problem to be solved again:
Figure BDA0003066598300000041
the method for formulating the upper and lower limits of the multi-microgrid operation reward area specifically comprises the following steps:
step 1), the total cost C of the power distribution network containing the non-private information protection microgrid is represented as follows:
C＝Cimp+Cshed+Cvol+Cmg
Figure BDA0003066598300000042
Figure BDA0003066598300000043
Figure BDA0003066598300000044
Figure BDA0003066598300000045
Figure BDA0003066598300000046
Cshed,m＝ρLSPLS,m
Figure BDA0003066598300000047
Cimpthe total cost of power input from the grid ($);
Cshedthe total load of the power distribution network reduces the cost ($);
Cvolthe cumulative voltage deviation cost ($) of the distribution network;
Cmgtotal microgrid operating cost ($);
Cgen,mconnected to bus-bar mTotal microgrid production cost ($);
Cshed,mthe total load of the microgrid connected to the bus m is reduced by a cost ($);
Cdeg,mtotal depreciation cost ($) for the microgrid connected to bus m;
ρE,mat boundary bus m, the price of the grid input power ($/kWh);
PB,mactive power (kW) is exchanged at a boundary bus m of the power distribution network and the power transmission network;
ρLScost ($) of load reduction in the distribution network;
PLS,mactive power demand reduction (kW) at bus m;
ρVpenalty factor (%) relating to the deviation of the bus voltage amplitude;
υmdeviation of the voltage amplitude at bus m ((kV)2)；
agConstant coefficient in the production cost function of the generator set ($/(kW)2)；
bgConstant coefficient ($/kW) in the generator set production cost function;
cga constant coefficient ($) in the generator set production cost function;
PG,gactive power output of g unit (kW)
ρLSCost ($) of load reduction in the distribution network;
PLS,mactive power demand reduction (kW) at bus m;
aea constant coefficient in a degradation cost function of the energy storage device e;
PCh,echarging power (kW) of the energy storage device e;
PDch,edischarge power (kW) of the energy storage means e;
step 2), after the total cost C of the power distribution network containing the non-privacy information protection microgrid is obtained, the unit power cost of the power distribution network containing the non-privacy information protection microgrid is expressed as follows:
Figure BDA0003066598300000051
lambda contains unit power cost ($/kW) of the power distribution network of the non-privacy information protection microgrid;
c Power distribution network total cost ($) with non-privacy information protection microgrid
PiThe non-privacy information protection type microgrid i has active power output (kW);
step 3), the total cost C' of the power distribution network with the privacy information protection microgrid is represented as follows:
C'＝C′imp+Cshed+Cvol+C'mg
Figure BDA0003066598300000052
Figure BDA0003066598300000053
Figure BDA0003066598300000054
step 4), after the total cost C' of the power distribution network containing the privacy information protection microgrid is obtained, the unit power cost of the power distribution network containing the privacy information protection microgrid is expressed as follows:
Figure BDA0003066598300000061
step 5), the unit power cost variation delta lambda of the power distribution network meets the preset setting value lambdaDIf the active output power of the privacy information protection type microgrid is within the allowable range, the cost deviation of other people caused by the active output power of the privacy information protection type microgrid is considered to be within the allowable range;
Δλ＝|λ-λ'|≤λD
the invention has the advantages that:
a power interaction model of a power distribution network in a Benders decomposition form and multi-microgrid operation is established, and a cooperative regulation and control method for a small-capacity microgrid with privacy information protection requirements is provided;
the operation efficiency and safety of the intelligent power distribution network can be improved, the privacy information of the microgrid main body is protected, and the solving process is accelerated while the optimality of the solution is ensured;
the cooperative regulation and control interval in the multi-microgrid autonomous space plan can be shortened, the constructed operation reward area is beneficial to seeking self benefit maximization of the multi-microgrid and meanwhile is also beneficial to electricity price compensation obtained by a long-acting reward mechanism, and therefore the multi-microgrid operation is promoted to approach to global optimization.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
Examples
The technical scheme of the invention is that a distributed cooperative regulation and control method comprises the following steps:
step 1, issuing an autonomous space plan by a power distribution network;
step 2, the power distribution network judges the operation interval in which the actual operation point of each micro-grid falls according to the autonomous space plan;
step 3, for the micro-grid with the operating point falling into the limit area, when a forcible limit measure is adopted to terminate the border crossing behavior, relevant information is counted and incorporated into a punishment mechanism; for the micro-grid with the operating point falling into the autonomous operating area, the power distribution network does not do forced interference, and the behavior information is collected and then brought into a reward mechanism; and for the micro-grid with the operating points falling into the cooperative control area, setting the boundary of the cooperative control area from outside to inside by using a Benders decomposition algorithm, so as to determine the boundary of the cooperative control area and obtain the upper limit and the lower limit of each-level operation reward area of the multi-micro-grid.
In the step 3, determining the boundary of the cooperative control region comprises the following steps:
step a),
Figure BDA0003066598300000071
Figure BDA0003066598300000072
Iterate to tBNext total operating cost ($);
Cimpthe total cost of power input from the grid ($);
Cshedthe total load of the power distribution network reduces the cost ($);
Cvolthe cumulative voltage deviation cost ($) of the distribution network;
zmthe microgrid running cost ($) at the position m of the power distribution network bus;
ΩMa node set connected with the microgrid in the power distribution network;
step b), iterating to tBNext time, the feasibility checking sub-problem of the microgrid at the bus m is expressed as:
Figure BDA0003066598300000073
Figure BDA0003066598300000074
Figure BDA0003066598300000075
Figure BDA0003066598300000076
actual active power exchange capacity (kW) determined in the main question;
Figure BDA0003066598300000077
the actual amount of reactive power exchanged (kVar) determined in the main problem;
sP1,madditional non-negative relaxation variables
sP2,mAn additional non-negative slack variable;
sQ1,madditional non-negative relaxation variables；
sQ2,mAn additional non-negative slack variable;
Figure BDA0003066598300000078
minimum sum of relaxation variables of the microgrid;
Figure BDA0003066598300000079
the bivariate associated with equation (1-2);
Figure BDA00030665983000000710
the bivariate associated with equations (1-3);
if the determined power exchange is feasible, all slack variables will take zero, so that the target value in equation (1-1) is zero, otherwise, the non-zero target value in equation (1-1) indicates that the resulting power exchange scheme is not feasible, which will correct the main problem in the form of a feasibility cut, in particular, equation (1-1) will be included in the main problem to be solved again:
Figure BDA0003066598300000081
step c), if the power exchange is feasible, solving an optimality check sub-problem by the multi-microgrid, as follows:
Figure BDA0003066598300000082
Figure BDA0003066598300000083
Figure BDA0003066598300000084
Figure BDA0003066598300000085
determining the minimum operation cost ($) of the microgrid under the power exchange schedule;
Figure BDA0003066598300000086
the sensitivity coefficients associated with equations (1-6), i.e., the Lagrangian multipliers;
Figure BDA0003066598300000087
the sensitivity coefficients associated with equations (1-7), i.e., the Lagrangian multipliers;
step d), after solving, utilizing the optimized target value
Figure BDA0003066598300000088
And the following valid inequality expressed as the optimality cut:
Figure BDA0003066598300000089
according to the weak duality of the convex optimization problem, an approximate value of the microgrid operation cost in the iteration process can be obtained through the formula;
step e), if the upper limit and the lower limit of the power exchange between the power distribution network and all the micro-grids are feasible, the upper limit of the operating cost of the power distribution network is expressed as follows:
Figure BDA00030665983000000810
Figure BDA00030665983000000811
at iteration tBAn upper bound ($) determined next time;
Figure BDA00030665983000000812
at iteration t, at the optimum value ofBDetermination of secondary main problem($)；
Figure BDA00030665983000000813
At iteration t, at the optimum value ofBDetermining ($) in a secondary main problem;
Figure BDA00030665983000000814
at iteration t, at the optimum value ofBDetermining ($) in a secondary main problem;
if the power exchange scheme is not feasible, its upper limit is set to infinity.
The method for formulating the upper and lower limits of the multi-microgrid operation reward area specifically comprises the following steps:
step 1), the total cost C of the power distribution network containing the non-private information protection microgrid is represented as follows:
C＝Cimp+Cshed+Cvol+Cmg
Figure BDA0003066598300000091
Figure BDA0003066598300000092
Figure BDA0003066598300000093
Figure BDA0003066598300000094
Figure BDA0003066598300000095
Cshed,m＝ρLSPLS,m
Figure BDA0003066598300000096
Cimpthe total cost of power input from the grid ($);
Cshedthe total load of the power distribution network reduces the cost ($);
Cvolthe cumulative voltage deviation cost ($) of the distribution network;
Cmgtotal microgrid operating cost ($);
Cgen,mtotal microgrid production cost ($) connected to bus m;
Cshed,mthe total load of the microgrid connected to the bus m is reduced by a cost ($);
Cdeg,mtotal depreciation cost ($) for the microgrid connected to bus m;
ρE,mat boundary bus m, the price of the grid input power ($/kWh);
PB,mactive power (kW) is exchanged at a boundary bus m of the power distribution network and the power transmission network;
ρLScost ($) of load reduction in the distribution network;
PLS,mactive power demand reduction (kW) at bus m;
ρVpenalty factor (%) relating to the deviation of the bus voltage amplitude;
υmdeviation of the voltage amplitude at bus m ((kV)2)；
agConstant coefficient in the production cost function of the generator set ($/(kW)2)；
bgConstant coefficient ($/kW) in the generator set production cost function;
cga constant coefficient ($) in the generator set production cost function;
PG,gactive power output of g unit (kW)
ρLSCost ($) of load reduction in the distribution network;
PLS,mactive power demand reduction (kW) at bus m;
aeenergy storageMeans e degenerates constant coefficients in the cost function;
PCh,echarging power (kW) of the energy storage device e;
PDch,edischarge power (kW) of the energy storage means e;
step 2), after the total cost C of the power distribution network containing the non-privacy information protection microgrid is obtained, the unit power cost of the power distribution network containing the non-privacy information protection microgrid is expressed as follows:
Figure BDA0003066598300000101
lambda contains unit power cost ($/kW) of the power distribution network of the non-privacy information protection microgrid;
c Power distribution network total cost ($) with non-privacy information protection microgrid
PiThe non-privacy information protection type microgrid i has active power output (kW);
step 3), the total cost C' of the power distribution network with the privacy information protection microgrid is represented as follows:
C'＝C′imp+Cshed+Cvol+C'mg
Figure BDA0003066598300000102
Figure BDA0003066598300000103
Figure BDA0003066598300000104
step 4), after the total cost C' of the power distribution network containing the privacy information protection microgrid is obtained, the unit power cost of the power distribution network containing the privacy information protection microgrid is expressed as follows:
Figure BDA0003066598300000105
step 5), the unit power cost variation delta lambda of the power distribution network meets the preset setting value lambdaDIf the active output power of the privacy information protection type microgrid is within the allowable range, the cost deviation of other people caused by the active output power of the privacy information protection type microgrid is considered to be within the allowable range;
Δλ＝|λ-λ'|≤λD
the privacy information and protection level of the invention mainly comprises:
(1) in the aspect of a power distribution network: the 1 st level protection information is power distribution network cost information (main network electricity purchasing cost and the like), the 2 nd level protection information is power distribution network allocable resource information (reactive compensation configuration, branch capacity margin, tie line state and the like), and the 3 rd level protection information is power distribution network topology information (branch power, node voltage, branch impedance and the like). The open information is power exchange request information and the like sent by the power distribution network to the microgrid.
(2) In the aspect of microgrid: the 1 st level protection information is microgrid cost information (distributed power generation cost, energy storage loss cost, microgrid operation and maintenance cost, load shedding cost and the like), the 2 nd level protection information is microgrid operation information (a daily power generation and utilization curve, an energy storage charging and discharging curve, load shedding information and the like), and the 3 rd level protection information is microgrid equipment and load information (distributed power capacity and parameters in the microgrid, energy storage capacity and parameters, interruptible load information and the like). The open information is power exchange request information sent by the microgrid to the power distribution network, and the like.
The above description is only for the preferred embodiment of the present invention, and is not intended to limit the present invention in any way. Any simple modification, change and equivalent changes of the above embodiments according to the technical essence of the invention are still within the protection scope of the technical solution of the invention.

Claims (4)

1. A distributed cooperative regulation and control method for a power distribution network is characterized by comprising the following steps:
step 1, issuing an autonomous space plan by a power distribution network;
step 2, the power distribution network judges the operation interval in which the actual operation point of each micro-grid falls according to the autonomous space plan;
step 3, for the micro-grid with the operating point falling into the limit area, when a forcible limit measure is adopted to terminate the border crossing behavior, relevant information is counted and incorporated into a punishment mechanism; for the micro-grid with the operating point falling into the autonomous operating area, the power distribution network does not do forced interference, and the behavior information is collected and then brought into a reward mechanism; and for the micro-grid with the operating points falling into the cooperative control area, setting the boundary of the cooperative control area from outside to inside by using a Benders decomposition algorithm, so as to determine the boundary of the cooperative control area and obtain the upper limit and the lower limit of each-level operation reward area of the multi-micro-grid.
2. The distributed cooperative control method for the power distribution network according to claim 1, wherein in the step 3, determining the boundary of the cooperative control region comprises the following steps:
step a),
Figure FDA0003066598290000011
Figure FDA0003066598290000012
Iterate to tBNext total operating cost ($);
Cimpthe total cost of power input from the grid ($);
Cshedthe total load of the power distribution network reduces the cost ($);
Cvolthe cumulative voltage deviation cost ($) of the distribution network;
zmthe microgrid running cost ($) at the position m of the power distribution network bus;
ΩMa node set connected with the microgrid in the power distribution network;
step b), iterating to tBNext time, the feasibility checking sub-problem of the microgrid at the bus m is expressed as:
Figure FDA0003066598290000013
Figure FDA0003066598290000014
Figure FDA0003066598290000015
Figure FDA0003066598290000021
actual active power exchange capacity (kW) determined in the main question;
Figure FDA0003066598290000022
the actual amount of reactive power exchanged (kVar) determined in the main problem;
sP1,madditional non-negative relaxation variables
sP2,mAn additional non-negative slack variable;
sQ1,man additional non-negative slack variable;
sQ2,man additional non-negative slack variable;
Figure FDA0003066598290000023
minimum sum of relaxation variables of the microgrid;
Figure FDA0003066598290000024
the bivariate associated with equation (1-2);
Figure FDA0003066598290000025
the bivariate associated with equations (1-3);
step c), if the power exchange is feasible, solving an optimality check sub-problem by the multi-microgrid, as follows:
Figure FDA0003066598290000026
Figure FDA0003066598290000027
Figure FDA0003066598290000028
Figure FDA0003066598290000029
determining the minimum operation cost ($) of the microgrid under the power exchange schedule;
Figure FDA00030665982900000210
the sensitivity coefficients associated with equations (1-6), i.e., the Lagrangian multipliers;
Figure FDA00030665982900000211
the sensitivity coefficients associated with equations (1-7), i.e., the Lagrangian multipliers;
step d), after solving, utilizing the optimized target value
Figure FDA00030665982900000212
And the following valid inequality expressed as the optimality cut:
Figure FDA00030665982900000213
according to the weak duality of the convex optimization problem, an approximate value of the microgrid operation cost in the iteration process can be obtained through the formula;
step e), if the upper limit and the lower limit of the power exchange between the power distribution network and all the micro-grids are feasible, the upper limit of the operating cost of the power distribution network is expressed as follows:
Figure FDA00030665982900000214
Figure FDA00030665982900000215
at iteration tBAn upper bound ($) determined next time;
Figure FDA0003066598290000031
——
Figure FDA0003066598290000032
at iteration t, at the optimum value ofBDetermining ($) in a secondary main problem;
Figure FDA0003066598290000033
——
Figure FDA0003066598290000034
at iteration t, at the optimum value ofBDetermining ($) in a secondary main problem;
Figure FDA0003066598290000035
——
Figure FDA0003066598290000036
at iteration t, at the optimum value ofBDetermining ($) in a secondary main problem;
if the power exchange scheme is not feasible, its upper limit is set to infinity.
3. The distributed coordinated control method for the power distribution network according to claim 2, wherein in the step b), if the determined power exchange is feasible, all the slack variables will be taken to be zero, so that the target value in the formula (1-1) is zero, otherwise, the non-zero target value in the formula (1-1) indicates that the obtained power exchange scheme is not feasible, which will correct the main problem in the form of a feasibility cut, specifically, the formula (1-1) will be included in the main problem to be solved again:
Figure FDA0003066598290000037
4. the distributed cooperative control method for the power distribution network according to claim 1, wherein the step of formulating the upper and lower limits of the multi-microgrid operation reward zone specifically comprises the following steps:
step 1), the total cost C of the power distribution network containing the non-private information protection microgrid is represented as follows:
C＝Cimp+Cshed+Cvol+Cmg
Figure FDA0003066598290000038
Figure FDA0003066598290000039
Figure FDA00030665982900000310
Figure FDA00030665982900000311
Figure FDA00030665982900000312
Cshed,m＝ρLSPLS,m
Figure FDA00030665982900000313
Cimpthe total cost of power input from the grid ($);
Cshedthe total load of the power distribution network reduces the cost ($);
Cvolthe cumulative voltage deviation cost ($) of the distribution network;
Cmgtotal microgrid operating cost ($);
Cgen,mtotal microgrid production cost ($) connected to bus m;
Cshed,mthe total load of the microgrid connected to the bus m is reduced by a cost ($);
Cdeg,mtotal depreciation cost ($) for the microgrid connected to bus m;
ρE,mat boundary bus m, the price of the grid input power ($/kWh);
PB,mactive power (kW) is exchanged at a boundary bus m of the power distribution network and the power transmission network;
ρLScost ($) of load reduction in the distribution network;
PLS,mactive power demand reduction (kW) at bus m;
ρVpenalty factor (%) relating to the deviation of the bus voltage amplitude;
υmdeviation of the voltage amplitude at bus m ((kV)2)；
agConstant coefficient in the production cost function of the generator set ($/(kW)2)；
bgConstant coefficient ($/kW) in the generator set production cost function;
cga constant coefficient ($) in the generator set production cost function;
PG,gactive power output of g unit (kW)
ρLSNegative in power distribution networkCost of load reduction ($);
PLS,mactive power demand reduction (kW) at bus m;
aea constant coefficient in a degradation cost function of the energy storage device e;
PCh,echarging power (kW) of the energy storage device e;
PDch,edischarge power (kW) of the energy storage means e;
step 2), after the total cost C of the power distribution network containing the non-privacy information protection microgrid is obtained, the unit power cost of the power distribution network containing the non-privacy information protection microgrid is expressed as follows:
Figure FDA0003066598290000041
lambda contains unit power cost ($/kW) of the power distribution network of the non-privacy information protection microgrid;
c Power distribution network total cost ($) with non-privacy information protection microgrid
PiThe non-privacy information protection type microgrid i has active power output (kW);
step 3), the total cost C' of the power distribution network with the privacy information protection microgrid is represented as follows:
C'＝C′imp+Cshed+Cvol+C'mg
Figure FDA0003066598290000051
Figure FDA0003066598290000052
Figure FDA0003066598290000053
step 4), after the total cost C' of the power distribution network containing the privacy information protection microgrid is obtained, the unit power cost of the power distribution network containing the privacy information protection microgrid is expressed as follows:
Figure FDA0003066598290000054
step 5), the unit power cost variation delta lambda of the power distribution network meets the preset setting value lambdaDIf the active output power of the privacy information protection type microgrid is within the allowable range, the cost deviation of other people caused by the active output power of the privacy information protection type microgrid is considered to be within the allowable range;
Δλ＝|λ-λ'|≤λD。
CN202110529223.7A 2021-05-14 2021-05-14 Distributed cooperative regulation and control method for power distribution network Active CN113346554B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110529223.7A CN113346554B (en) 2021-05-14 2021-05-14 Distributed cooperative regulation and control method for power distribution network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110529223.7A CN113346554B (en) 2021-05-14 2021-05-14 Distributed cooperative regulation and control method for power distribution network

Publications (2)

Publication Number Publication Date
CN113346554A true CN113346554A (en) 2021-09-03
CN113346554B CN113346554B (en) 2023-07-14

Family

ID=77470225

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110529223.7A Active CN113346554B (en) 2021-05-14 2021-05-14 Distributed cooperative regulation and control method for power distribution network

Country Status (1)

Country Link
CN (1) CN113346554B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116362523A (en) * 2023-06-01 2023-06-30 吉林大学 Coordinated optimization method for site selection and operation strategy of power exchange station considering temperature adaptability

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150019034A1 (en) * 2013-07-15 2015-01-15 Constantine Gonatas Device for Smoothing Fluctuations in Renewable Energy Power Production Cause by Dynamic Environmental Conditions
CN104617577A (en) * 2015-02-10 2015-05-13 清华大学 Active power distribution network full-distributed type self-discipline economic dispatch method
CN107292449A (en) * 2017-07-18 2017-10-24 广东双新电气科技有限公司 One kind is containing the scattered collaboration economic load dispatching method of many microgrid active distribution systems
CN109347095A (en) * 2018-10-29 2019-02-15 国网山东省电力公司枣庄供电公司 Consider the active distribution network layered distribution type tuning controller of characteristic of providing multiple forms of energy to complement each other
CN110011298A (en) * 2018-07-09 2019-07-12 东北林业大学 A kind of operation control strategy constructing the restructural microgrid group system of Autonomous Model
CN111769543A (en) * 2020-03-24 2020-10-13 绍兴大明电力设计院有限公司 Regional power distribution network autonomous cooperative operation optimization method containing multiple micro-grids
CN112542843A (en) * 2020-12-11 2021-03-23 南京南瑞继保工程技术有限公司 Independent micro-grid rapid frequency modulation method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150019034A1 (en) * 2013-07-15 2015-01-15 Constantine Gonatas Device for Smoothing Fluctuations in Renewable Energy Power Production Cause by Dynamic Environmental Conditions
CN104617577A (en) * 2015-02-10 2015-05-13 清华大学 Active power distribution network full-distributed type self-discipline economic dispatch method
CN107292449A (en) * 2017-07-18 2017-10-24 广东双新电气科技有限公司 One kind is containing the scattered collaboration economic load dispatching method of many microgrid active distribution systems
CN110011298A (en) * 2018-07-09 2019-07-12 东北林业大学 A kind of operation control strategy constructing the restructural microgrid group system of Autonomous Model
CN109347095A (en) * 2018-10-29 2019-02-15 国网山东省电力公司枣庄供电公司 Consider the active distribution network layered distribution type tuning controller of characteristic of providing multiple forms of energy to complement each other
CN111769543A (en) * 2020-03-24 2020-10-13 绍兴大明电力设计院有限公司 Regional power distribution network autonomous cooperative operation optimization method containing multiple micro-grids
CN112542843A (en) * 2020-12-11 2021-03-23 南京南瑞继保工程技术有限公司 Independent micro-grid rapid frequency modulation method

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116362523A (en) * 2023-06-01 2023-06-30 吉林大学 Coordinated optimization method for site selection and operation strategy of power exchange station considering temperature adaptability
CN116362523B (en) * 2023-06-01 2023-09-26 吉林大学 Coordinated optimization method for site selection and operation strategy of power exchange station considering temperature adaptability

Also Published As

Publication number Publication date
CN113346554B (en) 2023-07-14

Similar Documents

Publication Publication Date Title
CN108599373B (en) Cascade analysis method for transmission and distribution coordination scheduling target of high-proportion renewable energy power system
Liu et al. Uses, cost-benefit analysis, and markets of energy storage systems for electric grid applications
Rezaeimozafar et al. A review of behind-the-meter energy storage systems in smart grids
Mohammadkhani et al. Energy and emission management of CCHPs with electric and thermal energy storage and electric vehicle
Lv et al. Interactive energy management of networked microgrids-based active distribution system considering large-scale integration of renewable energy resources
Hu et al. Application of network-constrained transactive control to electric vehicle charging for secure grid operation
Dagdougui et al. Peak load reduction in a smart building integrating microgrid and V2B-based demand response scheme
Stenclik et al. Maintaining balance: the increasing role of energy storage for renewable integration
EP2688173B1 (en) Multi-service provision with energy storage system
Zhang et al. Grid-level application of electrical energy storage: Example use cases in the United States and China
CN110826880B (en) Active power distribution network optimal scheduling method for large-scale electric automobile access
Xu et al. Multi-objective chance-constrained optimal day-ahead scheduling considering BESS degradation
Zhang et al. Aggregating distributed energy storage: Cloud-based flexibility services from China
Gundogdu et al. Battery energy management strategies for UK firm frequency response services and energy arbitrage
Karandeh et al. Distributed control of behind-the-meter energy resources for multiple services
Ma et al. Energy storage capacity optimization for improving the autonomy of grid-connected microgrid
Singh et al. Profit maximization in ADN through voltage control and DR management with networked community micro-grids
CN113346554B (en) Distributed cooperative regulation and control method for power distribution network
Lu et al. Two-stage robust scheduling and real-time load control of community microgrid with multiple uncertainties
CN113762632A (en) Collaborative optimization operation method and system of electrical comprehensive energy system
Lee et al. Chance-constrained optimization for active distribution networks with virtual power lines
CN115660204A (en) Power grid congestion relief regulation and control method considering service type virtual power plant
Johnson et al. Utilization of stochastically located customer owned battery energy storage systems for violation management on UK LV residential feeders with varying renewables penetrations
Hu et al. Benefits Analysis of Energy Storage System in Power Systems with High Penetration of Renewable Energy
Brahmia et al. Efficient energy management of multi-microgrids systems based on distributed cooperative control in autonomous mode

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