CN109256770A - A kind of Distributed power net jamming control method based on Demand Side Response - Google Patents
A kind of Distributed power net jamming control method based on Demand Side Response Download PDFInfo
- Publication number
- CN109256770A CN109256770A CN201811168489.8A CN201811168489A CN109256770A CN 109256770 A CN109256770 A CN 109256770A CN 201811168489 A CN201811168489 A CN 201811168489A CN 109256770 A CN109256770 A CN 109256770A
- Authority
- CN
- China
- Prior art keywords
- power
- family
- moment
- branch
- node
- 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
Links
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
-
- 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/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- 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/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/14—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
-
- 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/18—Arrangements for adjusting, eliminating or compensating reactive power in networks
-
- 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/24—Arrangements for preventing or reducing oscillations of power in networks
-
- 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]
-
- 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/003—Load forecast, e.g. methods or systems for forecasting future load demand
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B70/00—Technologies for an efficient end-user side electric power management and consumption
- Y02B70/30—Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
- Y02B70/3225—Demand response systems, e.g. load shedding, peak shaving
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/30—Reactive power compensation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S20/00—Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
- Y04S20/20—End-user application control systems
- Y04S20/222—Demand response systems, e.g. load shedding, peak shaving
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The present invention proposes a kind of Distributed power net jamming control method based on Demand Side Response, belongs to Operation of Electric Systems and control technology field.This method comprises: establishing the electric appliance utility models in Demand Side Response;It establishes power distribution network network model and considers that congestion constrains;Establish load control system model;It establishes the Distributed power net Congestion Control Model objective function based on Demand Side Response and carries out distributed solution, real-time congestion control is carried out to power distribution network using every step solving result as control amount, until reaching the condition of convergence, the Distributed power net congestion control based on Demand Side Response terminates.The method of the present invention has given full play to Demand Side Response bring and has adjusted nargin, and Demand-side user is made to participate in the congestion control of network, is suitable for quick power distribution network congestion control, low in cost, is suitble to large-scale promotion.
Description
Technical field
The invention belongs to Operation of Electric Systems and control technology field, in particular to a kind of distribution based on Demand Side Response
Formula power distribution network jamming control method.
Background technique
With pushing forward comprehensively for smart grid construction, occurs a large amount of elastic load in power distribution network.This type load is one
Power distribution network scheduling can be responded under fixed condition, be increasingly becoming the important supplier of power distribution network ancillary service.Therefore, Demand Side Response
It is increasingly becoming engineering reality, just plays key player in each control algolithm of active distribution network.
Meanwhile as renewable energy permeability increasingly improves, the gradually development of power distribution network form, by traditional centralization
Power distribution network is gradually converted to multiagent active distribution network, and the development of distribution system is faced with new opportunities and challenges.If control
Method is improper, and power supply, load deteriorate the trend form of power distribution network, limited distribution line is competed, so that in power distribution network
The power delivery value of partial branch approaches the upper limit, is easy to cause distribution system failure, causes economic loss.Therefore, how to avoid
Above-mentioned this congestion phenomenon carries out reasonable congestion control to power distribution network, is a problem to be solved.
The congestion control of multiagent active distribution network, it is critical that, in the security constraint for meeting global system
Under the premise of (including node voltage and Branch Power Flow), by timely control strategy, flexibilities various in power distribution network are provided
Source is controlled, and realizes the minimum of power distribution network Congestion Level SPCC.Current jamming control method method, often failing to fully consider needs
Side response bring is asked to influence and act on, so that a large amount of controllable resources of Demand-side cannot give full play to the work of ancillary service
With.
Existing technology, such as a kind of active distribution network branch power congestion real-time control method, to active distribution network into
It has gone modeling, and has used and offline on Branch Power Flow power gathered around as constraint using the optimizer of centralization to model solution
Fill in control program.However, it in this method modeling process, there is no consideration Demand Side Responses, thus can not play Demand-side
The participation effect in active distribution network is responded, causes resulting congestion control scheme performance more low.
On the other hand, Demand-side number of resources is various, while power distribution network scale is growing also with development, traditional collection
Chinese style power distribution network uses the controller of centralization.Measurement information in distribution is all collected and is gone forward side by side by the controller of centralization
The centralization of row model solves, and result is then issued to each execution unit by communication network.Using centralized controller
The jamming control method of implementation is faced with the problem that response speed is slow, transinformation is big, can not carry out in time to network congestion
Control.Meanwhile in multiagent active distribution network, centralized algorithm will cause the trust between multiagent, coordination problem, difficult
With competent.Therefore, control is carried out to power distribution network congestion using distributed algorithm to be necessary.
Summary of the invention
The purpose of the present invention is the shortcomings to overcome prior art, propose a kind of distribution based on Demand Side Response
Power distribution network jamming control method.The present invention is based on Demand Side Responses, and it is important in power grid ancillary service sufficiently to excavate Demand-side
Effect, can reach better power distribution network control effect.
The present invention proposes a kind of Distributed power net jamming control method based on Demand Side Response, which is characterized in that should
Method the following steps are included:
1) the electric appliance utility models in Demand Side Response are established;
It enables the family in active distribution network collect and is combined into H={ h1,h2,...,hH, H is family's sum, remembers h-th of family
Electric appliance collection is combined into Ah={ ah,1,ah,2,...,ah,A, A is the electric appliance sum of h-th of family;
Define utility function Uh,a(ph,a) it is electric appliance a in h family in output power ph,aWhen bring effectiveness;
All electric appliances of h-th of family are divided into critical load and non-key load, then the output power of critical load meets
Following expression:
Wherein, ph,aIt (t) is the output power of a-th of electric appliance in h-th of family of t moment,It is t moment h-th
The necessary power that a-th of electric appliance works normally in front yard;
For non-key load, it is divided into interruptible load, transferable load and thermic load three classes:
Wherein, the effectiveness expression formula of interruptible load is as follows:
Wherein, T represents total moment number of operation;
The effectiveness expression formula of transferable load is as follows:
Meanwhile total electric energy of transferable load consumption meets following operation constraint:
Wherein, Eh,aThe total electric energy for representing a-th of the electric appliance of h-th of family needs lower limit,Represent the of h-th of family
Total electric energy of a electric appliance needs the upper limit;
The effectiveness expression formula of thermic load is as follows:
Wherein,For h-th of domestic consumer's the most comfortable temperature of t moment,It is the interior of h-th of family of t moment
Temperature;Meet following thermodynamical equilibrium equation:
Wherein, αh,βhIt is h-th of family thermal parameter relevant to thermal environment and thermal parameter relevant with electric appliance respectively;
Meet following constraint:
Wherein,It is the lower limit of h-th of home interior temperature of t moment,It is h-th of home interior temperature of t moment
The upper limit of degree;
There is following constraint for all electric appliances of h-th of family:
Power factor constraint:
Wherein, ηh,aFor the power factor of h-th of family, a-th of electric appliance, sh,aIt (t) is h-th of family, a-th of electric appliance in t
The apparent energy at moment;
Electric operation constraint:
Wherein, ph,a(t) h-th of family, a-th of electric appliance is represented in the lower limit of the power of t moment,Represent h-th of family
The upper limit of the power of a-th of the electric appliance in t moment;
Electrical work time-constrain:
Wherein,For the working time range of h-th of family, a-th of electric appliance;
2) power distribution network network model is established;
Active distribution network is modeled as corresponding connected graph (N, E), wherein N is node collection, and E is branch collection;Node 0 indicates
Feeder line root node, voltage V0It is set as the reference voltage controlled by transformer substation side, p0It (t) is t moment power distribution network from feeder line root section
The power that point obtains, each load node i ∈ N { 0 } to be connected to the family of the node provide electrical power, and t moment family
Active power required for h are as follows:
Reactive power required for t moment family h are as follows:
Wherein, qh,a(t) for a-th of electric appliance in h family t moment reactive power;
Net injection active power of the t moment in i-th of node φ phase:
Wherein, G is the set of power generation node, pGjIt (t) is active power of j-th of power generation node in t moment,It is
The injection active power of i critical load φ phase;
Net injection reactive power of the t moment in i-th of node φ phase:
Wherein, qGjIt (t) is reactive power of j-th of power generation node in t moment,For i-th critical load φ phase
Inject reactive power;
Three-phase distribution net power flow equation:
vj=vi-2(RijPij+XijQij)+Δvij(P,Q) (16)
WhereinFor the active power vector of Three-phase Power Flow ij branch,It is three
Phase trend ij branch reactive power vector,Active vector power is injected for the three-phase of j-th of node,Reactive power vector is injected for the three-phase of j-th of node,For j-th node
Voltage quantities,Node voltage, V are referred to for φ phaserefFor reference voltage virtual value;
For the active power loss of ij branch,For the idle network loss of ij branch, Δ vij(P, Q) is ij
The voltage drop of branch, specific linearized expression are as follows:
Wherein, P0For the initial active power vector of power distribution network, Q0For the initial reactive power vector of power distribution network, Pij0For power distribution network
The initial active power vector of ij branch, Qij0For the initial reactive power vector of power distribution network ij branch;
Wherein " .* " respectively represents two with "/" and is multiplied or is divided by with corresponding element between dimensional vector;
It is angularly measured at each node:
In formula, i is nodal scheme, and j is imaginary unit;
It defines the first kind and corrects impedance matrixImpedance matrix is corrected with the second classIn the i-th row j column element express difference
Are as follows:
Meanwhile the real and imaginary parts of first kind amendment impedance matrix respectively correspond first kind amendment resistor matrixWith first
Class corrects reactance matrixThe real and imaginary parts of second class amendment impedance matrix respectively correspond the second class amendment resistor matrixWith
Second class corrects reactance matrixExpression formula is as follows:
Wherein, partial derivative fpij、fqij、gpij、gqij、lpij、lqijExpression formula difference is as follows:
And
hxx(A, x)=diag (Ax)+diag (x) A (31)
hxy(A, x, y)=diag (Ay) (32)
hyx(A, x, y)=diag (y) A (33)
The constraint of the congestion of branch power and voltage magnitude:
Wherein,For the active power lower limit of ij branch φ phase,For the active power upper limit of ij branch φ phase,
For the active power of ij branch φ phase t moment;For the reactive power lower limit of ij branch φ phase,For the nothing of ij branch φ phase
The function upper limit of the power,For the reactive power of ij branch φ phase t moment;V iFor the lower voltage limit of ij branch,For ij branch
Upper voltage limit,For the voltage of ij branch φ phase t moment;
3) load control system model is established, expression formula is as follows:
Wherein,It is the total injection complex power upper limit of the power distribution network t moment in φ phase,At t moment node j
Injection active power,For the injection reactive power at t moment node j;
4) the Distributed power net Congestion Control Model objective function based on Demand Side Response is established, expression formula is as follows:
Under conditions of linear cost, CtFor power distribution company cost parameter, G is the set of power generation node, pGjIt (t) is jth
Active power of a power generation node in t moment;κ is weight of the power distribution company totle drilling cost in social welfare;
5) to model built (39), (1)-(38) carry out distributed solution, using every step solving result as control amount to matching
Power grid carries out real-time congestion control, until reaching the condition of convergence, the Distributed power net congestion control knot based on Demand Side Response
Beam;Specific step is as follows:
5-1) enable the initial step number k=1 of iteration;Choose Lagrange multiplier λh(t)、μh(t) initial value, give one it is positive
Penalty factor ρ > 0;
It 5-2) solves and matches net side subproblem:
s.t.(1)-(38)
Wherein L function follows following expression:
5-3) solve family subproblem:
s.t.(1)-(38)
Wherein, subscript k represents kth time iteration;
5-4) after kth time iteration, j-th of power station solves to obtainRepresent j-th of power station active power of output;The
H family solves to obtainAndRespectively represent family h distributed generation resource output active power and
Reactive power, family h each electric appliance output active power and reactive power, each node according to solving result control equipment pair
It should export, carry out the Distributed power net congestion control based on Demand Side Response;
5-5) update Lagrange multiplier:
5-6) judge whether iteration restrains according to the following formula:
Wherein ε is positive real number, represents convergence error;
If formula (44) is set up, iteration convergence, the Distributed power net congestion control based on Demand Side Response terminates;If formula
(44) invalid, then iteration does not restrain, and enables k=k+1, then returns to step 5-2).
Advantages of the present invention and beneficial effect are:
1. the method for the present invention has given full play to Demand Side Response bring and has adjusted nargin compared with traditional control method,
So that Demand-side user is participated in the congestion control of network, reduces the Congestion Level SPCC of network, and using distributed method to mould
Type carries out solution calculating, can quickly reduce local congestion's degree.
2. Distributed power net jamming control method of the invention, the independence of each family and hidden is maintained in distribution
It is private.In kth time iteration, the data volume exchanged between each family and distribution regulation center is that the general power of plan and boundary multiply
Sonλ0,μ0;And the private data of family, the details and customer priorities of each electrical equipment of example, all protection is each
In the household energy management system of family local.
3. the present invention needs not rely on the controller equiment of centralization, centralized communication network construction and dilatation had both been saved
Bring cost, and the reduction of the raising of control speed, information interaction amount is brought, it is suitable for quick power distribution network congestion control
System, it is low in cost, it is suitble to large-scale promotion.
Specific embodiment
The present invention proposes a kind of Distributed power net jamming control method based on Demand Side Response, below with reference to specific reality
Example is applied to be further described below.
The present invention proposes a kind of Distributed power net jamming control method based on Demand Side Response, comprising the following steps:
1) the electric appliance utility models in Demand Side Response are established;
It enables the family in active distribution network collect and is combined into H={ h1,h2,...,hH, H is family's sum, in each family hi,i
There is household energy management system (HEMS) in ∈ [1, H] to manage the household electrical appliance of this family, remembers the electric appliance set of h-th of family
For Ah={ ah,1,ah,2,...,ah,A, A is the electric appliance sum of h-th of family.In the load control system project based on contract, family
Front yard Energy Management System receives load control order from distribution regulation center.
Define utility function Uh,a(ph,a) it is electric appliance a in h family in output power ph,aWhen bring effectiveness.
All electric appliances of h-th of family, which are divided into critical load, (continuous work or must must ensure the symbol of power supply
Close, such as crucial illumination, refrigerator, cooker etc.) and non-key load, then the output power of critical load should meet following expression
Formula:
Wherein, ph,aIt (t) is the output power of a-th of electric appliance in h-th of family of t moment,It is t moment h-th
The necessary power that a-th of electric appliance works normally in front yard.
For remaining non-key load, it is divided into interruptible load, transferable load and thermic load three classes:
For interruptible load, such as optional illumination and plug and play load, bring effectiveness and t moment use electric work
Rate is related, expression formula:
Wherein T represent operation total moment number, for using hour as the stage, using day as the scene in period, often select T=
24。
For transferable load, such as washing machine and dryer, bring effectiveness depends on total electric energy of consumption, form
It is as follows:
Meanwhile total electric energy of transferable load consumption needs to meet certain operation constraint:
Wherein,E h,aThe total electric energy for representing a-th of the electric appliance of h-th of family needs lower limit,Represent the of h-th of family
Total electric energy of a electric appliance needs the upper limit.
For thermic load, such as air-conditioning and heating, room temperature is closer to h-th of domestic consumer's the most comfortable temperature of t moment
DegreeThen user utility is bigger, and bring effectiveness form is as follows:
Meanwhile h-th of home interior temperature of t moment meets following thermodynamical equilibrium equation:
Wherein, αh,βhIt is that h-th of family thermal parameter relevant to thermal environment and thermal parameter relevant with electric appliance are (usual respectively
Service manual by inquiring electric appliance obtains),It is the room temperature of h-th of family of t moment.SimultaneouslyIt should be in human body
Within the scope of tolerable:
WhereinIt is the lower limit of h-th of home interior temperature of t moment,It is h-th of home interior temperature of t moment
The upper limit.
There is following constraint for all electric appliances of h-th of family:
Power factor constraint:
Wherein, ηh,aFor the power factor of h-th of family, a-th of electric appliance, can be obtained from the service manual of electric appliance, sh,a
It (t) is apparent energy of h-th of family, a-th of the electric appliance in t moment.
Electric operation constraint:
Wherein, ph,a(t) h-th of family, a-th of electric appliance is represented in the lower limit of the power of t moment,Represent h-th of family
The upper limit of the power of a-th of the electric appliance in t moment.
Electrical work time-constrain:
Wherein,For the working time range of h-th of family, a-th of electric appliance.
2) power distribution network network model is established;
Active distribution network is modeled as corresponding connected graph (N, E), wherein N is node collection, and E is branch collection.Node 0 indicates
Feeder line root node, voltage V0It is set as the reference voltage controlled by transformer substation side, p0It (t) is t moment power distribution network from feeder line root section
The power that point obtains.Each load node i ∈ N { 0 } to be connected to the family of the node provide electrical power, and t moment family
Active power required for h are as follows:
Reactive power required for t moment family h are as follows:
Wherein qh,a(t) for a-th of electric appliance in h family t moment reactive power.
Net injection active power of the t moment in i-th of node φ phase:
Wherein, G is the set of power generation node, pGjIt (t) is active power of j-th of power generation node in t moment,It is
The injection active power of i critical load φ phase.
Net injection reactive power of the t moment in i-th of node φ phase:
Wherein, qGjIt (t) is reactive power of j-th of power generation node in t moment,For i-th critical load φ phase
Inject reactive power.
Three-phase distribution net power flow equation:
vj=vi-2(RijPij+XijQij)+Δvij(P,Q) (16)
WhereinFor the active P of Three-phase Power Flow ij branchij0Vector power,
For Three-phase Power Flow ij branch reactive power vector,For j-th node three-phase inject active power to
Amount,Reactive power vector is injected for the three-phase of j-th of node,For j-th of node
Voltage quantities,Node voltage, V are referred to for φ phaserefFor reference voltage virtual value.
For the active power loss of ij branch,For the idle network loss of ij branch, Δ vij(P, Q) is ij
The voltage drop of branch, specific linearized expression are as follows:
Wherein, P0The initial active power vector for being power distribution network before implementing this method, Q0Implementing this method for power distribution network
Preceding initial reactive power vector is initial active power vector of the power distribution network ij branch before implementing this method, Qij0For distribution
Net initial reactive power vector of the ij branch before implementing this method.
Wherein " .* " respectively represents two with "/" and is multiplied or is divided by with corresponding element between dimensional vector;
It is angularly measured at each node:
In formula, i is nodal scheme, and j is imaginary unit.
It defines the first kind and corrects impedance matrixImpedance matrix is corrected with the second classIn the i-th row j column element express difference
Are as follows:
Meanwhile the real and imaginary parts of first kind amendment impedance matrix respectively correspond first kind amendment resistor matrixWith first
Class corrects reactance matrixThe real and imaginary parts of second class amendment impedance matrix respectively correspond the second class amendment resistor matrixWith
Second class corrects reactance matrixExpression formula is as follows:
Wherein, partial derivative fpij、fqij、gpij、gqij、lpij、lqijExpression formula difference is as follows:
And
hxx(A, x)=diag (Ax)+diag (x) A (31)
hxy(A, x, y)=diag (Ay) (32)
hyx(A, x, y)=diag (y) A (33)
The constraint of the congestion of branch power and voltage magnitude:
Wherein,For the active power lower limit of ij branch φ phase,For the active power upper limit of ij branch φ phase,
For the active power of ij branch φ phase t moment;For the reactive power lower limit of ij branch φ phase,For the nothing of ij branch φ phase
The function upper limit of the power,For the reactive power of ij branch φ phase t moment;V iFor the lower voltage limit of ij branch,For ij branch
Upper voltage limit,For the voltage of ij branch φ phase t moment;
3) the load control system model i.e. restricted model of Demand Side Response is established, expression formula is as follows:
WhereinIt is the total injection complex power upper limit of the power distribution network t moment in φ phase,For the note at t moment node j
Enter active power,For the injection reactive power at t moment node j.
4) the Distributed power net Congestion Control Model objective function based on Demand Side Response is established, expression formula is as follows:
Under conditions of linear cost, CtFor power distribution company cost parameter, usually in positive real number domain value.
The objective function as shown in formula (39) is social welfare maximization, i.e. the total utility of family side user subtracts distribution public affairs
The totle drilling cost of department.Wherein κ is weight (usually in positive real number domain value) of the power distribution company totle drilling cost in social welfare, if it
Value is bigger, it is meant that compared with promoting user utility, power distribution company can more preferably press the operation cost for subtracting oneself.
5) to model built (39), (1)-(38) carry out distributed solution, using every step solving result as control amount to matching
Power grid carries out real-time congestion control, until reaching the condition of convergence, the Distributed power net congestion control knot based on Demand Side Response
Beam;Specific step is as follows:
5-1) enable the initial step number k=1 of iteration;One group of suitable parameter is chosen as Lagrange multiplier λh(t)、μh(t)
Initial value is usually initial value with 0.Give a positive penalty factor ρ > 0.
It 5-2) solves and matches net side subproblem:
s.t.(1)-(38)
Wherein L function follows following expression:
Wherein λh(t)、μhIt (t) is Lagrange multiplier, value changes with iteration.
5-3) solve family subproblem:
s.t.(1)-(38)
Wherein, subscript k represents kth time iteration.
5-4) after kth time iteration, j-th of power station solves to obtainRepresent j-th of power station active power of output;The
H family solves to obtainAndRespectively represent family h distributed generation resource output active power and
Reactive power, family h each electric appliance output active power and reactive power.Each node controls pair of equipment according to solving result
It should export, carry out the Distributed power net congestion control based on Demand Side Response.
5-5) update Lagrange multiplier:
5-6) according to convergence criterion, judge whether iteration restrains:
Wherein ε is a given positive real number, represents convergence error, should rationally adjust in practical applications, usually positive real
Number field value, typical value such as 1e-2.If formula (44) is set up, iteration convergence, the Distributed power net based on Demand Side Response
Congestion control terminates, and has reached control purpose;If formula (44) is invalid, iteration does not restrain, then enables k=k+1, then return again
Return step 5-2).
Claims (1)
1. a kind of Distributed power net jamming control method based on Demand Side Response, which is characterized in that this method includes following
Step:
1) the electric appliance utility models in Demand Side Response are established;
It enables the family in active distribution network collect and is combined into H={ h1,h2,...,hH, H is family's sum, remembers the electric appliance of h-th of family
Collection is combined into Ah={ ah,1,ah,2,...,ah,A, A is the electric appliance sum of h-th of family;
Define utility function Uh,a(ph,a) it is electric appliance a in h family in output power ph,aWhen bring effectiveness;
All electric appliances of h-th of family are divided into critical load and non-key load, then the output power of critical load meets as follows
Expression formula:
Wherein, ph,aIt (t) is the output power of a-th of electric appliance in h-th of family of t moment,For in h-th of family of t moment
The necessary power that a-th of electric appliance works normally;
For non-key load, it is divided into interruptible load, transferable load and thermic load three classes:
Wherein, the effectiveness expression formula of interruptible load is as follows:
Wherein, T represents total moment number of operation;
The effectiveness expression formula of transferable load is as follows:
Meanwhile total electric energy of transferable load consumption meets following operation constraint:
Wherein,E h,aThe total electric energy for representing a-th of the electric appliance of h-th of family needs lower limit,Represent a-th of h-th of family
Total electric energy of electric appliance needs the upper limit;
The effectiveness expression formula of thermic load is as follows:
Wherein,For h-th of domestic consumer's the most comfortable temperature of t moment,It is the Indoor Temperature of h-th of family of t moment
Degree;Meet following thermodynamical equilibrium equation:
Wherein, αh,βhIt is h-th of family thermal parameter relevant to thermal environment and thermal parameter relevant with electric appliance respectively;
Meet following constraint:
Wherein,It is the lower limit of h-th of home interior temperature of t moment,It is the upper of h-th of home interior temperature of t moment
Limit;
There is following constraint for all electric appliances of h-th of family:
Power factor constraint:
Wherein, ηh,aFor the power factor of h-th of family, a-th of electric appliance, sh,aIt (t) is h-th of family, a-th of electric appliance in t moment
Apparent energy;
Electric operation constraint:
Wherein,p h,a(t) h-th of family, a-th of electric appliance is represented in the lower limit of the power of t moment,Represent h-th of family a
The upper limit of the power of a electric appliance in t moment;
Electrical work time-constrain:
Wherein,For the working time range of h-th of family, a-th of electric appliance;
2) power distribution network network model is established;
Active distribution network is modeled as corresponding connected graph (N, E), wherein N is node collection, and E is branch collection;Node 0 indicates feeder line
Root node, voltage V0It is set as the reference voltage controlled by transformer substation side, p0It (t) is that t moment power distribution network is obtained from feeder line root node
The power taken, each load node i ∈ N { 0 } to be connected to the family of the node provide electrical power, and h institute of t moment family
The active power needed are as follows:
Reactive power required for t moment family h are as follows:
Wherein, qh,a(t) for a-th of electric appliance in h family t moment reactive power;
Net injection active power of the t moment in i-th of node φ phase:
Wherein, G is the set of power generation node, pGjIt (t) is active power of j-th of power generation node in t moment,It is i-th
The injection active power of critical load φ phase;
Net injection reactive power of the t moment in i-th of node φ phase:
Wherein, qGjIt (t) is reactive power of j-th of power generation node in t moment,For the injection of i-th of critical load φ phase
Reactive power;
Three-phase distribution net power flow equation:
vj=vi-2(RijPij+XijQij)+Δvij(P,Q) (16)
WhereinFor the active power vector of Three-phase Power Flow ij branch,For three-phase tide
Ij branch reactive power vector is flowed,Active vector power is injected for the three-phase of j-th of node,Reactive power vector is injected for the three-phase of j-th of node,For j-th node
Voltage quantities,Node voltage, V are referred to for φ phaserefFor reference voltage virtual value;
For the active power loss of ij branch,For the idle network loss of ij branch, Δ vij(P, Q) is ij branch
Voltage drop, specific linearized expression is as follows:
Wherein, P0For the initial active power vector of power distribution network, Q0For the initial reactive power vector of power distribution network, Pij0For power distribution network ij branch
The initial active power vector in road, Qij0For the initial reactive power vector of power distribution network ij branch;
Wherein " .* " respectively represents two with "/" and is multiplied or is divided by with corresponding element between dimensional vector;
It is angularly measured at each node:
In formula, i is nodal scheme, and j is imaginary unit;
It defines the first kind and corrects impedance matrixImpedance matrix is corrected with the second classIn the i-th row j column element expression be respectively as follows:
Meanwhile the real and imaginary parts of first kind amendment impedance matrix respectively correspond first kind amendment resistor matrixIt is repaired with the first kind
Positive reactive matrixThe real and imaginary parts of second class amendment impedance matrix respectively correspond the second class amendment resistor matrixWith second
Class corrects reactance matrixExpression formula is as follows:
Wherein, partial derivative fpij、fqij、gpij、gqij、lpij、lqijExpression formula difference is as follows:
And
hxx(A, x)=diag (Ax)+diag (x) A (31)
hxy(A, x, y)=diag (Ay) (32)
hyx(A, x, y)=diag (y) A (33)
The constraint of the congestion of branch power and voltage magnitude:
Wherein,For the active power lower limit of ij branch φ phase,For the active power upper limit of ij branch φ phase,For ij
The active power of branch φ phase t moment;For the reactive power lower limit of ij branch φ phase,For the idle function of ij branch φ phase
The rate upper limit,For the reactive power of ij branch φ phase t moment;V iFor the lower voltage limit of ij branch,For the voltage of ij branch
The upper limit,For the voltage of ij branch φ phase t moment;
3) load control system model is established, expression formula is as follows:
Wherein,It is the total injection complex power upper limit of the power distribution network t moment in φ phase,For the injection at t moment node j
Active power,For the injection reactive power at t moment node j;
4) the Distributed power net Congestion Control Model objective function based on Demand Side Response is established, expression formula is as follows:
Under conditions of linear cost, CtFor power distribution company cost parameter, G is the set of power generation node, pGjIt (t) is j-th of hair
Active power of the electrical nodes in t moment;κ is weight of the power distribution company totle drilling cost in social welfare;
5) to model built (39), (1)-(38) carry out distributed solution, using every step solving result as control amount to power distribution network
Real-time congestion control is carried out, until reaching the condition of convergence, the Distributed power net congestion control based on Demand Side Response terminates;Tool
Steps are as follows for body:
5-1) enable the initial step number k=1 of iteration;Choose Lagrange multiplier λh(t)、μh(t) initial value, given one positive penalize because
Sub- ρ > 0;
It 5-2) solves and matches net side subproblem:
s.t. (1)-(38)
Wherein L function follows following expression:
5-3) solve family subproblem:
s.t.(1)-(38)
Wherein, subscript k represents kth time iteration;
5-4) after kth time iteration, j-th of power station solves to obtainRepresent j-th of power station active power of output;H-th
Front yard solves to obtainAndRespectively represent the active power and idle function that the distributed generation resource of family h exports
Rate, family h each electric appliance output active power and reactive power, each node according to solving result control equipment correspondence it is defeated
Out, the Distributed power net congestion control based on Demand Side Response is carried out;
5-5) update Lagrange multiplier:
5-6) judge whether iteration restrains according to the following formula:
Wherein ε is positive real number, represents convergence error;
If formula (44) is set up, iteration convergence, the Distributed power net congestion control based on Demand Side Response terminates;If formula (44)
Invalid, then iteration does not restrain, and enables k=k+1, then returns to step 5-2).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811168489.8A CN109256770B (en) | 2018-10-08 | 2018-10-08 | Distributed power distribution network congestion control method based on demand side response |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811168489.8A CN109256770B (en) | 2018-10-08 | 2018-10-08 | Distributed power distribution network congestion control method based on demand side response |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109256770A true CN109256770A (en) | 2019-01-22 |
CN109256770B CN109256770B (en) | 2020-09-25 |
Family
ID=65045086
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811168489.8A Active CN109256770B (en) | 2018-10-08 | 2018-10-08 | Distributed power distribution network congestion control method based on demand side response |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109256770B (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103715687A (en) * | 2014-01-10 | 2014-04-09 | 贵州电网公司信息通信分公司 | Method for controlling branch power congestion of active power distribution network in real time |
US20140244792A1 (en) * | 2007-06-21 | 2014-08-28 | Jeffrey G. Fedders | Distributing intelligence across networks |
CN104734154A (en) * | 2014-12-18 | 2015-06-24 | 国家电网公司 | Multilevel hierarchical control method of transmission congestion management based on multi-source active distribution network |
CN107092985A (en) * | 2017-04-13 | 2017-08-25 | 大连理工大学 | A kind of active distribution network congestion Dispatching Method towards the flexible charge and discharge electroresponse of electric automobile group |
US20170324243A1 (en) * | 2016-05-09 | 2017-11-09 | Utopus Insights, Inc. | Congestion control in electric power system under load and uncertainty |
CN107359649A (en) * | 2017-08-31 | 2017-11-17 | 华中科技大学 | A kind of transmission & distribution Integrated optimization dispatching method considered under more distribution access situations |
CN107734020A (en) * | 2017-09-30 | 2018-02-23 | 国网青海省电力公司 | A kind of coordinated operation method of multiple photo-voltaic power generation station data transfer congestions |
US20180226800A1 (en) * | 2017-02-08 | 2018-08-09 | Battelle Energy Alliance, Llc | System, method, and apparatus for peak shaving and load shaping using power storage with variable power generation |
-
2018
- 2018-10-08 CN CN201811168489.8A patent/CN109256770B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140244792A1 (en) * | 2007-06-21 | 2014-08-28 | Jeffrey G. Fedders | Distributing intelligence across networks |
CN103715687A (en) * | 2014-01-10 | 2014-04-09 | 贵州电网公司信息通信分公司 | Method for controlling branch power congestion of active power distribution network in real time |
CN104734154A (en) * | 2014-12-18 | 2015-06-24 | 国家电网公司 | Multilevel hierarchical control method of transmission congestion management based on multi-source active distribution network |
US20170324243A1 (en) * | 2016-05-09 | 2017-11-09 | Utopus Insights, Inc. | Congestion control in electric power system under load and uncertainty |
US20180226800A1 (en) * | 2017-02-08 | 2018-08-09 | Battelle Energy Alliance, Llc | System, method, and apparatus for peak shaving and load shaping using power storage with variable power generation |
CN107092985A (en) * | 2017-04-13 | 2017-08-25 | 大连理工大学 | A kind of active distribution network congestion Dispatching Method towards the flexible charge and discharge electroresponse of electric automobile group |
CN107359649A (en) * | 2017-08-31 | 2017-11-17 | 华中科技大学 | A kind of transmission & distribution Integrated optimization dispatching method considered under more distribution access situations |
CN107734020A (en) * | 2017-09-30 | 2018-02-23 | 国网青海省电力公司 | A kind of coordinated operation method of multiple photo-voltaic power generation station data transfer congestions |
Non-Patent Citations (3)
Title |
---|
J. HAZRA ET AL.: "Smart Grid Congestion Management Through Demand Response", 《IEEE》 * |
JOHN W. HERON ET AL.: "Demand-Response Round-Trip Latency of IoT SmartGrid Network Topologies", 《IEEE》 * |
邰雪 等: "能源互联网中基于区块链的电力交易和阻塞管理方法", 《电网技术》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109256770B (en) | 2020-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bahrami et al. | Deep reinforcement learning for demand response in distribution networks | |
CN104734168B (en) | Microgrid running optimization system and method based on power and heat combined dispatching | |
CN105762841B (en) | A kind of parallel virtual synchronous generator distributed collaboration progress control method and system | |
Wang et al. | Privacy-preserving energy scheduling in microgrid systems | |
CN108734350A (en) | A kind of independent method for solving with combined dispatching of the power distribution network containing micro-capacitance sensor | |
CN109740786A (en) | A kind of multiterminal flexible interconnection distribution network planning method and system | |
CN103903073A (en) | Planning method and system for optimizing micro-grid containing distributed power sources and stored energy | |
CN108599154A (en) | A kind of three-phase imbalance power distribution network robust dynamic reconfiguration method considering uncertain budget | |
CN105870949B (en) | A kind of micro-capacitance sensor energy-storage units optimal control method based on distributed gradient algorithm | |
Wang et al. | Load curve smoothing strategy based on unified state model of different demand side resources | |
CN106786753A (en) | The system and its adjusting method of the Regional Energy internet of multi-user | |
CN104659810B (en) | A kind of micro-capacitance sensor cooperative control method for not knowing communication topology | |
CN107608218A (en) | Household electrical appliance intelligent power control method | |
CN105071410A (en) | Power distribution network reactive power optimized dispatching method and system based on virtual power plant | |
Wang et al. | A test bed for self-regulating distribution systems: Modeling integrated renewable energy and demand response in the GridLAB-D/MATLAB environment | |
CN106816879B (en) | Active power distribution network voltage coordination control method based on emergency demand response | |
CN111476394A (en) | Robust operation optimization method suitable for multi-energy systems such as electric, thermal and the like | |
CN109598429A (en) | A kind of island microgrid intelligent dispatching system and dispatching method | |
CN111144642B (en) | Comprehensive energy system optimization method and system | |
CN111461417A (en) | Comprehensive energy system collaborative optimization integrated modeling method based on large system theory | |
CN109256770A (en) | A kind of Distributed power net jamming control method based on Demand Side Response | |
Lu et al. | Planning of Soft Open Point Considering Demand Response | |
CN107171333A (en) | The real-time control method of virtual battery model based on intelligent load | |
Li et al. | Electricity ancillary service market model considering demand response of multi-energy load aggregators | |
Jiang et al. | Electricity optimal scheduling strategy considering multiple parks shared energy in the absence of grid power supply |
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 |