CN110232475A - A kind of distribution power transmission network power distribution network collaboration economic load dispatching method - Google Patents

A kind of distribution power transmission network power distribution network collaboration economic load dispatching method Download PDF

Info

Publication number
CN110232475A
CN110232475A CN201910458453.1A CN201910458453A CN110232475A CN 110232475 A CN110232475 A CN 110232475A CN 201910458453 A CN201910458453 A CN 201910458453A CN 110232475 A CN110232475 A CN 110232475A
Authority
CN
China
Prior art keywords
power
network
distribution network
power transmission
power distribution
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
CN201910458453.1A
Other languages
Chinese (zh)
Other versions
CN110232475B (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.)
Guangdong Power Grid Co Ltd
Zhuhai Power Supply Bureau of Guangdong Power Grid Co Ltd
Original Assignee
Guangdong Power Grid Co Ltd
Zhuhai Power Supply Bureau of Guangdong Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangdong Power Grid Co Ltd, Zhuhai Power Supply Bureau of Guangdong Power Grid Co Ltd filed Critical Guangdong Power Grid Co Ltd
Priority to CN201910458453.1A priority Critical patent/CN110232475B/en
Publication of CN110232475A publication Critical patent/CN110232475A/en
Application granted granted Critical
Publication of CN110232475B publication Critical patent/CN110232475B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0637Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06Electricity, gas or water supply
    • H02J3/383
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48Controlling the sharing of the in-phase component
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A30/00Adapting or protecting infrastructure or their operation
    • Y02A30/60Planning or developing urban green infrastructure
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Abstract

The present invention relates to the technical fields of electric system O&M, cooperate with economic load dispatching method more particularly, to a kind of distributed power transmission network power distribution network, comprising: establish network structure, and establish transmission & distribution Joint economics scheduling model based on network structure;Network structure includes power transmission network and n power distribution network, and scheduling model includes objective function and constraint condition;Outer circulation iterative solution, while F are carried out to F (x) based on Newton iterative method calculation methodTWithMiddle variable recycles the solution of initial value distributed iterative in carrying out;Until power transmission network and power distribution network n export respective scheduling result respectively when outer circulation iteration meets termination condition.The Newton iterative method that the present invention introduces outer ring in optimization calculates calculates and inner ring distributed iterative calculates, two layers of iterative problem is converted by power transmission network power distribution network collaboration dynamic economic dispatch optimization problem, since external iteration is Newton iterative method, with twice convergence, the accurate Efficient Solution of distribution of transmission & distribution cooperative scheduling optimization computational problem may be implemented.

Description

A kind of distribution power transmission network power distribution network collaboration economic load dispatching method
Technical field
The present invention relates to the technical fields of electric system O&M, more particularly, to a kind of distributed power transmission network power distribution network Cooperate with economic load dispatching method.
Background technique
As the permeability that distributed generation resource accesses the renewable energy such as power distribution network, especially photovoltaic, wind-powered electricity generation extensively improves, Influence of the power distribution network to power transmission network can not be ignored, and bring security risk for it.Meanwhile power transmission network is also to distributed generation resource Power distribution network combined operating has important coordinative role.The combined optimization of power transmission network power distribution network mostly uses greatly power transmission network power distribution network integrated The method for changing model splicing completes combined calculation, not only computationally intensive, solution efficiency is low, while being also required on distribution network operation business Such as model, operating parameter and real time data mass data are passed to transmission of electricity network operation business, communications burden is big.In addition, due to transmission of electricity Net and power distribution network adhere to different scheduling system independent operatings separately, it is difficult to directly analyze to the joint model, the privacy information of each main body It is unable to get protection.Therefore, it is necessary to design the combined optimization decoupling computation frame of efficient power transmission network power distribution network, it is also required to simultaneously More suitable algorithm meets quick iteration demand.
Currently, there is scholar by power transmission network power distribution network combined optimization problem based on distributed method solution, i.e. power transmission network and crowd Only need exchange boundary node information that calculating can be completed between polygamy power grid, each body interior solves respective subproblem, is not necessarily to Share all privacy informations.However, former PROBLEM DECOMPOSITION is passed through side between multiagent after multiple subproblems by most of these methods Boundary's information exchange carries out interior optimization respectively, but does not close optimization problem from global angle distich and analyzed, and not can guarantee and obtains And the centralized consistent globally optimal solution of model, there is the possibility for further promoting economy;And these methods are most of Some distributed methods are based on gradient method and realize iteration, and the number of iterations is higher, and power transmission network power distribution network can not be adapted in practical application The quick calculating demand of combined optimization.
Summary of the invention
It is an object of the invention to overcome the deficiencies of the prior art and provide a kind of distributed power transmission network power distribution network collaboration is economical Dispatching method can be quickly obtained global optimum's scheduling strategy, and can improve renewable energy in power distribution network using collaboration optimization The consumption in source is horizontal, improves the economic benefit of total system, guarantees the safe operation of power transmission network and power distribution network.
In order to solve the above technical problems, the technical solution adopted by the present invention is that:
A kind of distributed power transmission network power distribution network collaboration economic load dispatching method is provided, comprising the following steps:
S10. network structure is established, and transmission & distribution Joint economics scheduling model is established based on the network structure;The network Structure includes power transmission network and n power distribution network connecting with power transmission network, and the scheduling model indicates are as follows:
S.t.Ax-b=0
g(x)≤0
In formula, F (x) is objective function, FTIndicate the objective function of power transmission network,Indicate the power distribution network n of access power transmission network Objective function,For the distribution net collective for accessing power transmission network, Ax-b indicates active power balance equality constraint, and g (x) is indicated Run inequality constraints;
S20. outer circulation iterative solution, while F are carried out to F (x) based on Newton iterative method calculation methodTWithMiddle variable into Circulation initial value distributed iterative solves in row;Until power transmission network and power distribution network n difference are defeated when outer circulation iteration meets termination condition Respective scheduling result out.
Distributed power transmission network power distribution network of the invention cooperates with economic load dispatching method, in power transmission network power distribution network pool Optimization calculate in introduce the Newton iterative method of outer ring and calculate and inner ring distributed iterative calculates, power transmission network power distribution network is cooperateed with dynamic State economic load dispatching optimization problem is converted into two layers of iterative problem, since external iteration is Newton iterative method, has twice convergence, The accurate Efficient Solution of distribution that transmission & distribution cooperative scheduling optimization computational problem may be implemented, in addition without power distribution network and power transmission network All privacy informations are shared, the demand of data exchange is reduced, improve the stability of calculating.
Preferably, pass through interconnection l between n power distribution network and power transmission network respectively1、l2、l3、…、lnConnection, interconnection l1、 l2、l3、…、lnIt is t in the node of power distribution network1、t2、t3、…、tn, interconnection l1、l2、l3、…、lnIn power distribution network 1, power distribution network 2 ..., the endpoint of power distribution network n is expressed as d1、d2、d3、…、dn
Preferably, the objective function F of the power transmission networkTWith the objective function of the power distribution network n of access power transmission networkTable respectively It is shown as:
In formula,For the distribution net collective for accessing power transmission network, K is that power transmission network power distribution network combines the total of dynamic economic dispatch When number of segment,It is the set of power transmission network generating set, ai、biAnd ciThe respectively power generation consumption cost coefficient of power transmission network unit, ai、biAnd ciValue be respectively 0.05~1,10~11 and 400~2000,It is i-th unit of power transmission network moment t's Plan active power;an,i、bn,iAnd cn,iI-th conventional distributed power generation unit consumption cost coefficient in respectively power distribution network n, an,i、bn,iAnd cn,iValue be respectively 0.01~1,30~35 and 550~1400, αn,iFor i-th in power distribution network n distributed light Overhead generator abandons light penalty coefficient, and value is 700~1000, pn,tNode for power distribution network n and power transmission network boundary in moment t is electric Valence,The active power for being the j node that is connect in power distribution network n with power transmission network in moment t to power transmission network planned purchases,For In power distribution network n i-th distributed conventional generator moment t plan active power,For i-th distribution in power distribution network n Formula photovoltaic generator moment t plan active power,It is i-th distributed photovoltaic power generation machine in power distribution network n in moment t Prediction active power,For conventional generator set distributed in power distribution network n,For distributed photovoltaic in power distribution network n Generating equipment set.
Preferably, the active power balance equality constraint includes power transmission network active power balance constraint and power distribution network Active power balance constraint, the operation inequality constraints includes that power transmission network operation constraint and power distribution network operation constrain, described defeated Operation of power networks constraint includes power grids circuits power constraint and power transmission network generator active power output constraint, the power distribution network fortune Row constraint includes the constraint of power distribution network conventional generator active power output, the constraint of power distribution network photovoltaic power generation active power output, distribution network line Tie-line power transmission constrains between power constraint and power distribution network and power transmission network.
Preferably, the power transmission network active power balance constraint indicates are as follows:
In formula,For the node set being connect in power transmission network with i-node,For i-th of node in power transmission network when Carve the prediction load active power of t;
The power distribution network active power balance constraint indicates are as follows:
In formula,WithT moment node i flows to node u in respectively power distribution network n and node j flows to node i Active power,The line loss of node i is flowed to for t moment node j in power distribution network n,WithRespectively power distribution network In n t moment node i flow to node i history run is active and reactive power flow,It is gone through for t moment node j in power distribution network n History working voltage,(i, j) branch resistance in power distribution network n,I-node to connect in power distribution network n with power transmission network exists Active power from moment t to power transmission network planned purchases,Distributed power generation plan for t moment node i in power distribution network n is active Power,For the load active power predicted value of t moment node i in power distribution network n.
Preferably, Ax-b=0 described in step S20 is by power transmission network active power balance constraint and power distribution network active power Unknown variable concentrates the matrix expression rewritten and obtained in Constraints of Equilibrium.
Preferably, power grids circuits powerConstraint representation are as follows:
In formula,WithRespectively branch active power lower and upper limit of the power transmission network from node i to node j;
The output of power transmission network generator active powerConstraint representation are as follows:
In formula,WithThe respectively active power bottoming and the upper limit of i-th generating set of power transmission network.
Preferably, power distribution network conventional generator active power outputConstraint representation are as follows:
In formula,WithI-th conventional generator active power output bound in respectively power distribution network n;
The power distribution network photovoltaic power generation active power outputConstraint representation are as follows:
In formula,For the history power generation active power output of t moment i-node in power distribution network n;
The distribution network line powerConstraint representation are as follows:
In formula,WithRespectivelyTo flow to the active power of node i branch most from node j in power distribution network n Small value and maximum value;
Tie-line power transmission between the power distribution network and power transmission networkConstraint representation are as follows:
In formula,WithRespectively power distribution network n interacts active power with power transmission network minimum and maximum.
Preferably, the power transmission network is run into constraint and power distribution network operation constraint is written as about the abstract of column vector x Function representation form are as follows:
G (x)=[gtrans(x),gdist,n(x)]T≤0
In formula, gtransIt (x) is inequality constraints column vector abstract function expression formula in power transmission network, gdist,nIt (x) is power distribution network Inequality constraints column vector abstract function expression formula in n.
Preferably, solution procedure follows these steps to carry out in step S20:
S21. setting subscript trans represents the correlated variables in power transmission network, and subscript dist, n represent the correlation in power distribution network n Variable;Set subscriptkFor the number of iterations variable, the iteration subscript of functionkIt then indicates to bring x intokValue afterwards;Column vector is set to become Amount v represents the Lagrange multiplier variable of equality constraint, and column vector variable λ represents the corresponding Lagrange multiplier of inequality constraints Variable;Set yk=[xkk,vk]TWith Δ yk=[Δ xk,Δλk,Δvk]T;Feasible initial point y is set0=[x00,v0]TWith And Δ y0Whole element values are 0;Setting the number of iterations is k=0;
S22. in outer circulation kth time iteration, power transmission network variableInitial value is recycled in power distribution network nIt is distributed Iterative solution;
S23. power transmission network calculates
Power distribution network n is calculatedIn formula
S24. power transmission network calculates Power distribution network n is calculated
S25. power transmission network updatesPower distribution network n updates
S26. power distribution network updatesWithAnd after being sent to power transmission network, power transmission network calculates convergence signal W And power distribution network n is passed to, in formulaε2For outer circulation convergence Standard, generally desirable 10-4;If W=1, terminate outer circulation, it is active that power transmission network and power distribution network n export each self generating sets respectively Power output plan;If W=0, step S22 is gone to.
Compared with prior art, the beneficial effects of the present invention are:
Distributed power transmission network power distribution network of the invention cooperates with economic load dispatching method, in power transmission network power distribution network pool Optimization calculate in introduce the Newton iterative method of outer ring and calculate and inner ring distributed iterative calculates, power transmission network power distribution network is cooperateed with dynamic State economic load dispatching optimization problem is converted into two layers of iterative problem, since external iteration is Newton iterative method, has twice convergence, The accurate Efficient Solution of distribution that transmission & distribution cooperative scheduling optimization computational problem may be implemented, in addition without power distribution network and power transmission network All privacy informations are shared, the demand of data exchange is reduced, improve the stability of calculating;And method of the invention can be fast Speed obtains global optimum's scheduling strategy, and the consumption that renewable energy in power distribution network can be improved using collaboration optimization is horizontal, mentions The economic benefit of high total system guarantees the safe operation of power transmission network and power distribution network.
Detailed description of the invention
Fig. 1 is the flow chart that distributed power transmission network power distribution network of the invention cooperates with economic load dispatching method;
Fig. 2 is the schematic diagram that distributed power transmission network power distribution network cooperates with network structure in economic load dispatching method;
Fig. 3 is that distributed power transmission network power distribution network cooperates with solution procedure flow chart in economic load dispatching method and step S20.
Specific embodiment
The present invention is further illustrated With reference to embodiment.
Embodiment
It is as shown in Figure 1 the embodiment of distributed power transmission network power distribution network collaboration economic load dispatching method of the invention, including with Lower step:
S10. network structure is established, and transmission & distribution Joint economics scheduling model is established based on network structure;Network structure includes Power transmission network and n power distribution network connecting with power transmission network, scheduling model indicate are as follows:
S.t.Ax-b=0
g(x)≤0
In formula, F (x) is objective function, and F indicates the objective function of power transmission network,Indicate the power distribution network n of access power transmission network Objective function,For the distribution net collective for accessing power transmission network, Ax-b indicates active power balance equality constraint, and g (x) is indicated Run inequality constraints;
S20. outer circulation iterative solution, while F are carried out to F (x) based on Newton iterative method calculation methodTAnd FDIn n variable into Circulation initial value distributed iterative solves in row;Until power transmission network and power distribution network n difference are defeated when outer circulation iteration meets termination condition Respective scheduling result out.
As shown in Fig. 2, the network structure in the present embodiment includes power transmission network and n power distribution network connecting with power transmission network, n is a Pass through interconnection l between power distribution network and power transmission network respectively1、l2、l3、…、lnConnection, interconnection l1、l2、l3、…、lnIn power distribution network Node be t1、t2、t3、…、tn, interconnection l1、l2、l3、…、lnPower distribution network 1, power distribution network 2 ..., the endpoint of power distribution network n point D is not expressed as it1、d2、d3、…、dn
Objective function in the present embodiment indicates are as follows:
In formula,For the distribution net collective for accessing power transmission network, K is that power transmission network power distribution network combines the total of dynamic economic dispatch When number of segment,It is the set of power transmission network generating set, ai、biAnd ciThe respectively power generation consumption cost coefficient of power transmission network unit, ai、biAnd ciValue be respectively 0.05~1,10~11 and 400~2000,It is i-th unit of power transmission network moment t's Plan active power.an,i、bn,iAnd cn,iI-th conventional distributed power generation unit consumption cost coefficient in respectively power distribution network n, an,i、bn,iAnd cn,iValue be respectively 0.01~1,30~35 and 550~1400, αn,iFor i-th in power distribution network n distributed light Overhead generator abandons light penalty coefficient, value 700-1000, pn,tNode for power distribution network n and power transmission network boundary in moment t is electric Valence,The active power for being the j node that is connect in power distribution network n with power transmission network in moment t to power transmission network planned purchases,For In power distribution network n i-th distributed conventional generator moment t plan active power,For i-th distribution in power distribution network n Formula photovoltaic generator moment t plan active power,It is i-th distributed photovoltaic power generation machine in power distribution network n in moment t Prediction active power,For conventional generator set distributed in power distribution network n,For distributed photovoltaic in power distribution network n Generating equipment set.
In the present embodiment, active power balance equality constraint includes power transmission network active power balance constraint and power distribution network Active power balance constraint, operation inequality constraints include power transmission network operation constraint and power distribution network operation constraint, in which:
Power transmission network active power balance constraint based on DC flow model indicates are as follows:
In formula,For the node set being connect in power transmission network with i-node,It is i-th of node in power transmission network at the moment The prediction load active power of t.
Consider the power distribution network active power balance constraint of via net loss
In formula,WithT moment node i flows to node u in respectively power distribution network n and node j flows to node i Active power,The line loss of node i is flowed to for t moment node j in power distribution network n,WithRespectively power distribution network In n t moment node i flow to node i history run is active and reactive power flow,It is gone through for t moment node j in power distribution network n History working voltage,(i, j) branch resistance in power distribution network n,I-node to connect in power distribution network n with power transmission network exists Active power from moment t to power transmission network planned purchases,It is (including normal for the distributed power generation of t moment node i in power distribution network n Advise generator and photovoltaic power generation equipment) plan active power,Load active power for t moment node i in power distribution network n is pre- Measured value.
The unknown variable of above-mentioned power transmission network active power balance constraint and power distribution network active power balance constraint is concentrated, and It is rewritten into following matrix form:
Ax=b
In formula, A represents expanding node branch matrix,Define Atrans= [Ctrans,Itrans], Adist,n=[Bdist,n,Cdist,n,Idist,n], [Bdist,n,Cdist,n,Idist,n] it is [Bdist Cdist Idist] The matrix intersection of the middle part corresponding power distribution network n, Matrix CtransIn arbitrary element beMatrix ItransIn arbitrary element beMatrix BtransIn arbitrary element beMatrix BdistIn arbitrary element beMatrix CdistIn arbitrary element beMatrix IdistMiddle arbitrary element are as follows:
Column vector x, which is represented, to be needed Variables collection is sought, column vector x is by power transmission network column vector variable x to be askedtransWith power distribution network n column vector variable x to be askeddist,nIt constitutes, Column vector xtransBy power transmission network branch active power column vector variable to be askedWith Power transmission network generator active power column vector variable to be askedConstitute, arrange to Measure xdistBy power distribution network interconnection active power column vector variable to be askedPower distribution network branch active power column vector variable to be askedWith power transmission network generator active power column vector variable to be askedIt constitutes, column vector b Represent extension matrix of loadings, b=[(btrans)T, (bdist)T]T, subscriptTRepresent transposition, column vector btransMiddle arbitrary element isColumn vector bdistBy the column vector b of each power distribution networkdist,nCollection is combined into, column vector bdist,nIn any member Element is
Power transmission network operation constraint includes power grids circuits power constraint and power transmission network generator active power output constraint, is matched Operation of power networks constraint includes the constraint of power distribution network conventional generator active power output, the constraint of power distribution network photovoltaic power generation active power output, distribution Tie-line power transmission constrains between cable road power constraint and power distribution network and power transmission network;Specifically:
Power grids circuits power constraint indicates are as follows:
In formula,WithRespectively branch active power lower and upper limit of the power transmission network from node i to node j;
Power transmission network generator active power output constraint indicates are as follows:
In formula,WithThe respectively active power bottoming and the upper limit of i-th generating set of power transmission network;
Power distribution network conventional generator active power output constraint representation are as follows:
In formula,WithI-th conventional generator active power output bound in respectively power distribution network n;
Power distribution network photovoltaic power generation active power output constraint representation are as follows:
In formula,For the history power generation active power output of t moment i-node in power distribution network n;
Distribution network line power constraint indicates are as follows:
In formula,WithRespectivelyTo flow to the active power of node i branch most from node j in power distribution network n Small value and maximum value.
Tie-line power transmission constrains between power distribution network and power transmission network
In formula,WithRespectively power distribution network n interacts active power with power transmission network minimum and maximum.
Power transmission network is run into constraint and power distribution network operation constrains the abstract function expression-form being written as about column vector x are as follows:
G (x)=[gtrans(x),gdist,n(x)]T≤0
In formula, gtransIt (x) is inequality constraints column vector abstract function expression formula in power transmission network, gdist,nIt (x) is power distribution network Inequality constraints column vector abstract function expression formula in n.
In addition, solution procedure follows these steps to carry out in step S20, as shown in Figure 3:
S21. setting subscript trans represents the correlated variables in power transmission network, and subscript dist, n represent the correlation in power distribution network n Variable;Set subscriptkFor the number of iterations variable, the iteration subscript of functionkIt then indicates to bring x intokValue afterwards;Column vector is set to become Amount v represents the Lagrange multiplier variable of equality constraint, and column vector variable λ represents the corresponding Lagrange multiplier of inequality constraints Variable;Set yk=[xkk,vk]TWith Δ yk=[Δ xk,Δλk,Δvk]T;Feasible initial point y is set0=[x00,v0]TWith And Δ y0Whole element values are 0;Setting the number of iterations is k=0;
S22. in outer circulation kth time iteration, power transmission network variableInitial value is recycled in power distribution network nDistribution Formula iterative solution;The following steps are included:
(1) initial value is recycled in setting power transmission networkInitial value is recycled in power distribution network nLoop iteration in being arranged Number is z=0;
(2) it recycles in the z times iteration inside, power transmission network receives each power distribution networkIt calculatesPower distribution network n receives power transmission network simultaneouslyIt calculatesIn formulaRepresent power transmission network variableInside follow Value of the ring at the z times,Represent power distribution network n variableInside circulate in the z times value, subscript-1Represent matrix inversion It calculates, matrix Mtrans、Mdist,nWithAnd column vector stransAnd sdist,nCalculation method narration such as Under:
Order matrix Y=AH-1AT, in formula, H=▽2Fk-(Dgk)Tdiag(gk)-1diag(λk)Dgk, ▽2Representative is found a function Hessian matrix, that is, secondary local derviation matrix, Dg represent Jacobian matrix operator i.e.▽ represents gradient Operator, gqFor q-th of daughter element function of column vector function g, diag is diagonal matrix function, that is, is generated with each member of independent variable Element is the diagonal matrix of diagonal element.
It can be found that H is symmetrical diagonal matrix, H diagonal element can be acquired by power transmission network and power distribution network local computing:
Y is symmetrical sparse matrix, and the quantity of row or column is power transmission network and the total nodal point number amount of power distribution network association system, in formula The element being not zero shows that the node of this element node of the row and column has in power transmission network power distribution network association system Physical couplings.
ThenIn formulaFor diagonal matrix, only power transmission network interior joint i is correspondingElement, εiRepresent all nodes being connected with node i in power transmission network power distribution network association system, ω be less than 1 is greater than 0.5 decomposition coefficient, and generally desirable 2/3,
In formulaFor diagonal matrix, only power distribution network n interior joint i is correspondingElement,
For diagonal matrix, only power transmission network interior joint i is correspondingElement.For Diagonal matrix, only power distribution network n node i are correspondingElement.
Column vectorFor power transmission network kth time iteration raw residual, in formulaMiddle member Element isVectorFor power distribution network n kth time iteration raw residual.
Column vectorFor power transmission network kth time iteration antithesis residual error.To AmountFor distribution n net kth time iteration antithesis residual error, in formula,
Column vectorFor power transmission network kth time iteration center residual error, τ in formula=μm/η, μ For parameter, general value is the inequality quantity that 10, m is inequality column vector function g, ηk=-(λk)Tgk, 1 is and column vector The consistent unit vector of function g line number, after τ informs power transmission network by power distribution network interface boundary information, power transmission network is completed to calculate anti-again It feeds power transmission network.VectorFor distribution n net kth time iteration center residual error.
Column vectorIn formula Column vector
(3) loop convergence criterion epsilon in being arranged1, generally desirable 10-6.Power transmission network calculatesMatch Power grid n is calculatedIfOrZ=z+1 is then set and goes to step 2-2- 2), ifAndThen terminate interior circulation to enableWithAnd go to step S23.
S23. power transmission network calculatesPower distribution network n is calculatedIn formula
S24. power transmission network calculates Power distribution network n is calculated
S25. power transmission network updatesPower distribution network n updates
S26. power distribution network updatesWithAnd after being sent to power transmission network, power transmission network calculates convergence signal W And power distribution network n is passed to, in formulaε2For outer circulation convergence Standard, generally desirable 10-4;If W=1, terminate outer circulation, it is active that power transmission network and power distribution network n export each self generating sets respectively Power output plan;If W=0, step S22 is gone to.
By above step, the present embodiment introduces outer ring in the optimization of power transmission network power distribution network pool calculates Newton iterative method calculates and inner ring distributed iterative calculates, by the collaboration dynamic economic dispatch optimization problem conversion of power transmission network power distribution network For two layers of iterative problem, the accurate Efficient Solution of distribution of transmission & distribution cooperative scheduling optimization computational problem may be implemented, and without matching Power grid and power transmission network share all privacy informations, reduce the demand of data exchange, improve the stability of calculating.
Obviously, the above embodiment of the present invention be only to clearly illustrate example of the present invention, and not be pair The restriction of embodiments of the present invention.For those of ordinary skill in the art, may be used also on the basis of the above description To make other variations or changes in different ways.There is no necessity and possibility to exhaust all the enbodiments.It is all this Made any modifications, equivalent replacements, and improvements etc., should be included in the claims in the present invention within the spirit and principle of invention Protection scope within.

Claims (10)

1. a kind of distribution power transmission network power distribution network cooperates with economic load dispatching method, which comprises the following steps:
S10. network structure is established, and transmission & distribution Joint economics scheduling model is established based on the network structure;The network structure Including power transmission network and n power distribution network connecting with power transmission network, the scheduling model is indicated are as follows:
S.t.Ax-b=0
g(x)≤0
In formula, F (x) is objective function, FTIndicate the objective function of power transmission network,Indicate the mesh of the power distribution network n of access power transmission network Scalar functions,For the distribution net collective for accessing power transmission network, Ax-b=0 indicates active power balance equality constraint, g (x) ≤ 0 indicates operation inequality constraints condition;
S20. outer circulation iterative solution, while F are carried out to F (x) based on Newton iterative method calculation methodTWithIn middle variable carries out Initial value distributed iterative is recycled to solve;Until power transmission network and power distribution network n are exported respectively respectively when outer circulation iteration meets termination condition Self scheduling result.
2. distribution power transmission network power distribution network according to claim 1 cooperates with economic load dispatching method, which is characterized in that n are matched Pass through interconnection l between power grid and power transmission network respectively1、l2、l3、…、lnConnection, interconnection l1、l2、l3、…、lnIn power distribution network Node is t1、t2、t3、…、tn, interconnection l1、l2、l3、…、lnPower distribution network 1, power distribution network 2 ..., the endpoint of power distribution network n difference It is expressed as d1、d2、d3、…、dn
3. distribution power transmission network power distribution network according to claim 1 cooperates with economic load dispatching method, which is characterized in that described defeated The objective function F of power gridTWith the objective function of the power distribution network n of access power transmission networkIt respectively indicates are as follows:
In formula,For the distribution net collective for accessing power transmission network, K is total period that power transmission network power distribution network combines dynamic economic dispatch Number,It is the set of power transmission network generating set, ai、biAnd ciThe respectively power generation consumption cost coefficient of power transmission network unit, ai、bi And ciValue be respectively 0.05~1,10~11 and 400~2000,For i-th unit of power transmission network moment t plan Active power;an,i、bn,iAnd cn,iI-th conventional distributed power generation unit consumption cost coefficient, a in respectively power distribution network nn,i、 bn,iAnd cn,iValue be respectively 0.01~1,30~35 and 550~1400, αn,iFor i-th distributed photovoltaic hair in power distribution network n Motor abandons light penalty coefficient, and value is 700~1000, pn,tFor power distribution network n and power transmission network boundary moment t node electricity price,The active power for being the j node that is connect in power distribution network n with power transmission network in moment t to power transmission network planned purchases,For with In power grid n i-th distributed conventional generator moment t plan active power,For i-th distribution in power distribution network n Photovoltaic generator moment t plan active power,It is i-th distributed photovoltaic power generation machine in power distribution network n moment t's Predict active power,For conventional generator set distributed in power distribution network n,It is sent out for distributed photovoltaic in power distribution network n Electric equipment set.
4. distribution power transmission network power distribution network according to claim 1 cooperates with economic load dispatching method, which is characterized in that described to have Function power-balance equality constraint includes power transmission network active power balance constraint and power distribution network active power balance constraint, described Operation inequality constraints includes power transmission network operation constraint and power distribution network operation constraint, and the power transmission network operation constraint includes power transmission network Line power constraint and power transmission network generator active power output constraint, the power distribution network operation constraint include that power distribution network is routinely sent out The constraint of motor active power output, the constraint of power distribution network photovoltaic power generation active power output, distribution network line power constraint and power distribution network and transmission of electricity Tie-line power transmission constrains between net.
5. distribution power transmission network power distribution network according to claim 4 cooperates with economic load dispatching method, which is characterized in that described defeated Network re-active power Constraints of Equilibrium indicates are as follows:
In formula,For the node set being connect in power transmission network with i-node,It is i-th of node in power transmission network moment t's Predict load active power;
The power distribution network active power balance constraint indicates are as follows:
In formula,WithT moment node i flows to node u in respectively power distribution network n and node j flows to the active of node i Power,The line loss of node i is flowed to for t moment node j in power distribution network n,WithT in respectively power distribution network n Moment node i flow to node i history run is active and reactive power flow,For the history fortune of t moment node j in power distribution network n Row voltage,(i, j) branch resistance in power distribution network n,It is the i-node that is connect in power distribution network n with power transmission network in moment t To the active power of power transmission network planned purchases,For the distributed power generation plan active power of t moment node i in power distribution network n,For the load active power predicted value of t moment node i in power distribution network n.
6. distribution power transmission network power distribution network according to claim 5 cooperates with economic load dispatching method, which is characterized in that step Ax-b=0 described in S20 is by becoming in power transmission network active power balance constraint and power distribution network active power balance constraint wait ask The matrix expression rewritten in quantity set.
7. cooperateing with economic load dispatching method according to the described in any item distributed power transmission network power distribution networks of claim 4 to 6, feature exists In:
Power grids circuits powerConstraint representation are as follows:
In formula,WithRespectively branch active power lower and upper limit of the power transmission network from node i to node j;
The output of power transmission network generator active powerConstraint representation are as follows:
In formula,WithThe respectively active power bottoming and the upper limit of i-th generating set of power transmission network.
8. distribution power transmission network power distribution network according to claim 7 cooperates with economic load dispatching method, it is characterised in that:
Power distribution network conventional generator active power outputConstraint representation are as follows:
In formula,WithI-th conventional generator active power output bound in respectively power distribution network n;
The power distribution network photovoltaic power generation active power outputConstraint representation are as follows:
In formula,For the history power generation active power output of t moment i-node in power distribution network n;
The distribution network line powerConstraint representation are as follows:
In formula,WithRespectivelyFor the active power minimum value for flowing to node i branch in power distribution network n from node j And maximum value;
Tie-line power transmission between the power distribution network and power transmission networkConstraint representation are as follows:
In formula,WithRespectively power distribution network n interacts active power with power transmission network minimum and maximum.
9. distribution power transmission network power distribution network according to claim 8 cooperates with economic load dispatching method, which is characterized in that will be described Power transmission network operation constraint and power distribution network operation constraint are written as the abstract function expression-form about column vector x are as follows:
G (x)=[gtrans(x),gdist,n(x)]T≤0
In formula, gtransIt (x) is inequality constraints column vector abstract function expression formula in power transmission network, gdist,n(x) in power distribution network n Inequality constraints column vector abstract function expression formula.
10. distribution power transmission network power distribution network according to claim 1 cooperates with economic load dispatching method, which is characterized in that step Solution procedure follows these steps to carry out in S20:
S21. subscript is settransRepresent the correlated variables in power transmission network, subscriptdist,nRepresent the correlated variables in power distribution network n;If Determine subscriptkFor the number of iterations variable, the iteration subscript of functionkIt then indicates to bring x intokValue afterwards;Column vector variable v is set to represent The Lagrange multiplier variable of equality constraint, column vector variable λ represent the corresponding Lagrange multiplier variable of inequality constraints;If Determine yk=[xkk,vk]TWith Δ yk=[Δ xk,Δλk,Δvk]T;Feasible initial point y is set0=[x00,v0]TAnd Δ y0Entirely Portion's element value is 0;Setting the number of iterations is k=0;
S22. in outer circulation kth time iteration, power transmission network variableInitial value is recycled in power distribution network nDistributed iterative It solves;
S23. power transmission network calculates
Power distribution network n is calculatedIn formula
S24. power transmission network calculatesDistribution N is netted to calculate
S25. power transmission network updatesPower distribution network n updates
S26. power distribution network updatesWithAnd after being sent to power transmission network, power transmission network calculates convergence signal W and transmits Power distribution network n is given, in formula ε2For outer circulation convergence, generally desirable 10-4;If W=1, terminate outer circulation, power transmission network and power distribution network n are exported respectively respectively The power output plan of self generating sets active power;If W=0, step S22 is gone to.
CN201910458453.1A 2019-05-29 2019-05-29 Distributed transmission network power distribution network coordinated economic dispatching method Active CN110232475B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910458453.1A CN110232475B (en) 2019-05-29 2019-05-29 Distributed transmission network power distribution network coordinated economic dispatching method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910458453.1A CN110232475B (en) 2019-05-29 2019-05-29 Distributed transmission network power distribution network coordinated economic dispatching method

Publications (2)

Publication Number Publication Date
CN110232475A true CN110232475A (en) 2019-09-13
CN110232475B CN110232475B (en) 2021-07-13

Family

ID=67858679

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910458453.1A Active CN110232475B (en) 2019-05-29 2019-05-29 Distributed transmission network power distribution network coordinated economic dispatching method

Country Status (1)

Country Link
CN (1) CN110232475B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111555370A (en) * 2020-05-20 2020-08-18 云南电网有限责任公司电力科学研究院 Power distribution network layered coordination scheduling method and device based on cloud edge coordination
CN111612248A (en) * 2020-05-20 2020-09-01 云南电网有限责任公司电力科学研究院 Power distribution network side source-load coordination method and system
CN111654025A (en) * 2020-06-12 2020-09-11 国网河北省电力有限公司经济技术研究院 Transmission and distribution integrated coordination scheduling method considering interconnection between distribution networks
CN112003277A (en) * 2020-08-21 2020-11-27 山东大学 Transmission and distribution cooperative load recovery optimization control method and system
CN112165086A (en) * 2020-08-12 2021-01-01 华南理工大学 Online optimization system of active power distribution network
CN112308331A (en) * 2020-11-09 2021-02-02 中国南方电网有限责任公司 Optimal power flow model in power market separated from power distribution network and iterative solution method thereof
CN116722605A (en) * 2023-08-08 2023-09-08 武汉创星空间科技发展有限公司 Power distribution network scheduling optimization method based on Internet of things

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150310366A1 (en) * 2012-11-09 2015-10-29 Tianjin University Security region based security-constrained economic dispatching method
CN107069706A (en) * 2017-02-17 2017-08-18 清华大学 A kind of dynamic economic dispatch method that transmission and distribution network based on multi-parametric programming is coordinated

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150310366A1 (en) * 2012-11-09 2015-10-29 Tianjin University Security region based security-constrained economic dispatching method
CN107069706A (en) * 2017-02-17 2017-08-18 清华大学 A kind of dynamic economic dispatch method that transmission and distribution network based on multi-parametric programming is coordinated

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张旭 等: ""高比例可再生能源电力系统的输配协同优化调度方法"", 《电力系统自动化》 *
林青 等: ""输配电网一体化分布式潮流计算方法"", 《广东电力》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111555370A (en) * 2020-05-20 2020-08-18 云南电网有限责任公司电力科学研究院 Power distribution network layered coordination scheduling method and device based on cloud edge coordination
CN111612248A (en) * 2020-05-20 2020-09-01 云南电网有限责任公司电力科学研究院 Power distribution network side source-load coordination method and system
CN111612248B (en) * 2020-05-20 2023-09-08 云南电网有限责任公司电力科学研究院 Power distribution network side source-load coordination method and system
CN111654025A (en) * 2020-06-12 2020-09-11 国网河北省电力有限公司经济技术研究院 Transmission and distribution integrated coordination scheduling method considering interconnection between distribution networks
CN111654025B (en) * 2020-06-12 2021-10-08 国网河北省电力有限公司经济技术研究院 Transmission and distribution integrated coordination scheduling method considering interconnection between distribution networks
CN112165086A (en) * 2020-08-12 2021-01-01 华南理工大学 Online optimization system of active power distribution network
CN112003277A (en) * 2020-08-21 2020-11-27 山东大学 Transmission and distribution cooperative load recovery optimization control method and system
CN112308331A (en) * 2020-11-09 2021-02-02 中国南方电网有限责任公司 Optimal power flow model in power market separated from power distribution network and iterative solution method thereof
CN116722605A (en) * 2023-08-08 2023-09-08 武汉创星空间科技发展有限公司 Power distribution network scheduling optimization method based on Internet of things
CN116722605B (en) * 2023-08-08 2023-10-17 武汉创星空间科技发展有限公司 Power distribution network scheduling optimization method based on Internet of things

Also Published As

Publication number Publication date
CN110232475B (en) 2021-07-13

Similar Documents

Publication Publication Date Title
CN110232475A (en) A kind of distribution power transmission network power distribution network collaboration economic load dispatching method
Li et al. Multiagent-based distributed state of charge balancing control for distributed energy storage units in AC microgrids
Zhao et al. Coordinated restoration of transmission and distribution system using decentralized scheme
Wu et al. Distributed optimal coordination for distributed energy resources in power systems
CN109950907B (en) Dispatching method and system for alternating current-direct current hybrid power distribution network containing power electronic transformer
CN110266038B (en) Distributed coordination regulation and control method for multiple virtual power plants
CN110504685A (en) A kind of Polynuclear complex system control parameters optimization method
CN109474017A (en) A kind of real-time distributed economic load dispatching method of power distribution network
CN107706921B (en) Micro-grid voltage regulation method and device based on Nash game
CN105703368B (en) Multi-uncertainty energy flow modeling integrated with active power distribution network and transmission network
CN102801790A (en) Microgrid energy management system based on IEC (International Electrotechnical Commission) 61850 and multi-agent method
CN105406459B (en) Utilize the method for the optimal load flow estimated based on uniformity distributed treatment in power network
CN107769213A (en) A kind of more transverters alternating current-direct current distribution power system load flow calculation method in parallel
Zhang et al. Distributed optimal active power dispatch under constraints for smart grids
CN105119279B (en) A kind of distributed power source planing method and its system
Yang et al. Deep learning-based distributed optimal control for wide area energy Internet
CN106684932B (en) A kind of interconnected network economic load dispatching method considering electric network composition adjustment
CN106549376B (en) The branch comprehensive stability appraisal procedure of power distribution network containing DG based on equivalent node method
CN110120673A (en) Distributed transmission & distribution collaboration idle work optimization method and system based on Thevenin's equivalence parameter identification
CN106877309A (en) Stability judgment method for direct-current micro-grid in master-slave control mode
CN111654025B (en) Transmission and distribution integrated coordination scheduling method considering interconnection between distribution networks
CN107196339A (en) Microgrid stratified sampling distributed and coordinated control method based on event trigger mechanism
CN110535123A (en) A kind of more microgrid economic load dispatching optimization methods of quick analytic expression distributed multiple target
CN110048407A (en) Distributed energy power generation plan feasible zone method for optimization analysis
CN113673912A (en) Distribution-gas network distributed collaborative planning method and system considering influence of power transmission network

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