CN104868480B - The active distribution network robust voltage control method that a kind of active reactive is coordinated - Google Patents

The active distribution network robust voltage control method that a kind of active reactive is coordinated Download PDF

Info

Publication number
CN104868480B
CN104868480B CN201510226565.6A CN201510226565A CN104868480B CN 104868480 B CN104868480 B CN 104868480B CN 201510226565 A CN201510226565 A CN 201510226565A CN 104868480 B CN104868480 B CN 104868480B
Authority
CN
China
Prior art keywords
distribution network
active
active distribution
power
voltage
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510226565.6A
Other languages
Chinese (zh)
Other versions
CN104868480A (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.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN201510226565.6A priority Critical patent/CN104868480B/en
Publication of CN104868480A publication Critical patent/CN104868480A/en
Application granted granted Critical
Publication of CN104868480B publication Critical patent/CN104868480B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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/30Reactive power compensation

Abstract

The present invention relates to the active distribution network robust voltage control method that a kind of active reactive is coordinated, belong to electric power system optimization operation control field.This method establishes the object function of optimal control, propose the interval scheduling method of photovoltaic generation, Three-phase Power Flow equation in being run using active distribution network under voltage highest and minimum scene, and subdivided capacitor is handled using binary coding mode, by discrete control problem continuous treatment.Finally, in Comprehensive Control active distribution network photovoltaic, Continuous Var Compensation device, the reactive power of subdivided capacitor and to issue permission active power to photovoltaic plant interval, to adapt to the uncertainty of photovoltaic generation and load power in active distribution network.The Optimal Control Strategy that this method is drawn ensure that voltage is safe when photovoltaic generation and load power are fluctuated in indeterminacy section, be conducive to the voltage of active distribution network to control.

Description

The active distribution network robust voltage control method that a kind of active reactive is coordinated
Technical field
The present invention relates to the active distribution network robust voltage control method that a kind of active reactive is coordinated, belong to power system excellent Change operation control field.
Background technology
The power distribution network accessed for high permeability photovoltaic, what the size of photovoltaic generation power was brought to the safe operation of power network Influence can not be ignored.Photovoltaic generation power is sized depending on intensity of illumination, therefore with randomness and intermittence.Distribution network planning Mould is huge, is limited to investment and maintenance cost, seldom installs real-time measurement device, the measuring value for state estimation is generally load work( Rate puppet is measured, and therefore, the ornamental of power distribution network is poor, and state estimation result particularly load estimation error is very big.It is traditional true Qualitative voltage optimization control is all based on imaginary " accurate " data, does not consider the uncertainty of photovoltaic generation and estimating for load Error is counted, then its obtained optimisation strategy does not ensure that power network when photovoltaic generation and load power are fluctuated in indeterminacy section Voltage is still safe.
The content of the invention
The purpose of the present invention is to propose to the active distribution network robust voltage control method that a kind of active reactive is coordinated, according to every The active power that control centre is sent on individual photovoltaic plant is interval interval with load apparent energy, and proposition can adapt to photovoltaic generation With the probabilistic active reactive robust voltage control method of load power, it is ensured that photovoltaic generation and load power are in uncertain region Between fluctuate when voltage it is safe.
The active distribution network robust voltage control method that active reactive proposed by the present invention is coordinated, comprises the following steps:
(1) the generating active power interval for presetting the photovoltaic plant in an access active distribution network isAnd Active distribution network control centre will be sent on default active power interval, active distribution network control centre is by photovoltaic plant The permission of active power is intervalPhotovoltaic plant is issued to, when the active power of photovoltaic plant is more than above-mentioned permission interval The upper limitWhen, the active power for making photovoltaic plant is the upper limit
(2) the voltage-controlled object function of active distribution network robust for setting up active reactive coordination is as follows:
Wherein, M abandons light penalty coefficient for photovoltaic plant, and M span is 4~6,P pvWithRespectively photovoltaic electric The lower and upper limit of the active power of active distribution network control centre are sent on standing,WithRespectively active distribution network is dispatched Center is issued to the lower and upper limit of the active power of photovoltaic plant,For each node active power composition in active distribution network Column vector, footmark u represents to make active distribution network magnitude of voltage to be in maximum, footmarkFor the ABC three-phases of active distribution network, For the column vector of each node active power composition in active distribution network, footmark l represents to make active distribution network magnitude of voltage be in minimum Value, sum () represents to sum to the element in column vector;
(3) set active distribution network voltage be in maximum when trend constraint as:
Wherein, footmarkFor the ABC three-phases of active distribution network,WithThe active respectively formed by topology search is matched somebody with somebody The three-phase node branch road incidence matrix and branch impedance matrix of power network,ForIntermediate value be -1 element constitute withDimension phase Same end branch incidence matrix,ForIntermediate value be 1 element constitute withDimension identical top branch road incidence matrix, Real () and imag () represent to take real and imaginary part respectively, × and × it is respectively the sum of products dot product of matrix,WithThe column vector of the column vector of each node active power composition and reactive power composition, angle respectively in active distribution network Mark u represents to make active distribution network magnitude of voltage in maximum,
Wherein,The column vector constituted for the active power lower limit of each load bus of active distribution network,Actively to match somebody with somebody Each photovoltaic plant allows the active power upper limit in power networkThe column vector of composition,For the nothing of each load bus of active distribution network The column vector of work(lower limit of the power composition, QinFor by the reactive power of all photovoltaic plants in active distribution networkStatic reacance The reactive power of compensatorWith the reactive power of subdivided capacitorThe column vector of composition;
For the column vector of active distribution network voltage value three phases active power composition of branch road when maximum:
Wherein, N is the circuitry number in active distribution network,WithRepresent active distribution network magnitude of voltage most respectively When big on branch road 1 A, B and C three-phase active power, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value three phase reactive power composition of branch road when maximum:
Wherein, N is the circuitry number in active distribution network,Active distribution network magnitude of voltage is represented respectively When maximum on branch road 1 A, B and C three-phase reactive power, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value electric current modulus value square composition of branch road when maximum:
Wherein, N is the circuitry number in active distribution network,WithRepresent that active distribution network magnitude of voltage is maximum respectively Square of A, B and C three-phase current modulus value on Shi Zhilu 1, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value three-phase voltage modulus value square composition of node when maximum:
Wherein, n is the nodes in active distribution network,Represent that active distribution network magnitude of voltage is maximum respectively Square of Shi Jiedian 1 A, B and C three-phase voltage modulus value, subscript T represents the transposition computing of vector;
(4) set active distribution network voltage be in minimum value when trend constraint as:
Wherein, footmarkFor the ABC three-phases of active distribution network,WithThe active respectively formed by topology search is matched somebody with somebody The three-phase node branch road incidence matrix and branch impedance matrix of power network,ForIntermediate value be -1 element constitute withDimension phase Same end branch incidence matrix,ForIntermediate value be 1 element constitute withDimension identical top branch road incidence matrix, Real () and imag () represent to take real and imaginary part respectively, × and × it is respectively the sum of products dot product of matrix,WithThe row arrow of the column vector of each node active power composition and each node reactive power composition respectively in active distribution network Amount, footmark l represents to make active distribution network magnitude of voltage in minimum value,
Wherein,The column vector constituted for the active power upper limit of each load bus of active distribution network,Actively to match somebody with somebody Each photovoltaic plant allows active power lower limit in power networkThe column vector of composition,For each load bus of active distribution network The column vector of reactive power upper limit composition, QinFor by the reactive power of all photovoltaic plants in active distribution networkStatic nothing The reactive power of work(compensatorWith the reactive power of subdivided capacitorThe column vector of composition;
For the column vector of active distribution network voltage value three phases active power composition of branch road when minimum:
Wherein, N is the circuitry number in active distribution network,WithRepresent that active distribution network magnitude of voltage is minimum respectively The active power of A, B and C three-phase on Shi Zhilu 1, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value three phase reactive power composition of branch road when minimum:
Wherein, N is the circuitry number in active distribution network,Active distribution network magnitude of voltage is represented respectively When minimum on branch road 1 A, B and C three-phase reactive power, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value electric current modulus value square composition of branch road when minimum:
Wherein, N is the circuitry number in active distribution network,WithRepresent that active distribution network magnitude of voltage is minimum respectively Square of A, B and C three-phase current modulus value on Shi Zhilu 1, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value three-phase voltage modulus value square composition of node when minimum:
Wherein, n is the nodes in active distribution network,Represent that active distribution network magnitude of voltage is minimum respectively Square of Shi Jiedian 1 A, B and C three-phase voltage modulus value, subscript T represents the transposition computing of vector;
(5) active power of photovoltaic plant is constrained in setting active distribution network:
Wherein,P pvWithRespectively be sent on photovoltaic plant active distribution network control centre active power lower limit and The upper limit,WithRespectively active distribution network control centre is issued to the lower and upper limit of the active power of photovoltaic plant;
(6) voltage security of setting active distribution network interior joint is constrained to:
Wherein,WithThe three-phase voltage modulus value of root node is put down when respectively active distribution network voltage is maximum The three-phase voltage modulus value square of root node when side and minimum active distribution network voltage,It is maximum for active distribution network voltage When other nodes in addition to root node three-phase voltage modulus value square,Root section is removed during for active distribution network voltage minimum The three-phase voltage modulus value square of other nodes outside point, UmaxFor the upper safety limit of active distribution network interior joint three-phase voltage, Umin For the lower safety limit of active distribution network interior joint three-phase voltage;
(7) active power and reactive power between setting power transmission network and active distribution network on transmission line are constrained to:
Wherein,WithActive distribution network root node three-phase when respectively active distribution network magnitude of voltage is maximum The three-dimensional column vector of active power composition and the three-dimensional column vector of three phase reactive power composition,WithRespectively The three-dimensional column vector and three-phase of active distribution network root node three phases active power composition are idle when active distribution network magnitude of voltage is minimum The three-dimensional column vector of power composition,WithRespectively between power transmission network and power distribution network on transmission line active power peace Full lower and upper limit,WithRespectively between power transmission network and power distribution network on transmission line the lower safety limit of reactive power and The upper limit;
(8) reactive capability of photovoltaic plant and SVC is constrained in setting active distribution network:
Wherein,WithThe reactive power of photovoltaic plant and SVC respectively in active distribution network,The respectively reactive power maximum size and lower limit of photovoltaic plant,Respectively static reacance The reactive power maximum size and lower limit of compensator;
(9) reactive power of subdivided capacitor is constrained in setting active distribution network:
T=α0*201*21
Wherein,For the reactive power of subdivided capacitor,For the idle work(in subdivided capacitor per group capacitor Rate, the group number that t puts into operation for subdivided capacitor in active distribution network, α0、α1Respectively variable, value is 0 or 1;
(10) Primal-dual Interior Point is used, according to above-mentioned steps (3), the trend constraint of (4) and step (5)~step Suddenly the operation constraint of the active distribution network of (9), the object function of solution procedure (2), calculating obtains photovoltaic electric in active distribution network Reactive power value, the reactive power value of SVC, the reactive power value of subdivided capacitor and the photovoltaic plant stood have The permission of work(power is interval, realizes that active distribution network robust voltage is controlled.
The active distribution network robust voltage control method that active reactive proposed by the present invention is coordinated, its advantage is:
1st, the inventive method considers the randomness and intermittence that photovoltaic generation is accessed in active distribution network, considers simultaneously The error that deficiency causes load power estimation to exist is measured in active distribution network, the active distribution network robust electricity that active reactive is coordinated The control strategy that pressure control method is drawn ensure that photovoltaic generation and the load power voltage in uncertain interval interior fluctuation It is safe.The Primal-dual Interior Point used in the inventive method, can effectively handle nonconvex nonlinear programming problems.
2nd, the inventive method proposes photovoltaic generation Robust Interval scheduling method, for the active distribution of high permeability photovoltaic Net, it is simple not necessarily to can guarantee that system voltage safety, it is necessary to be realized with reference to a certain amount of light of abandoning by adjustment reactive power Voltage security.Active power forecast interval is sent out on photovoltaic plant, control centre's optimization issues permission active power after calculating interval To photovoltaic plant.
3rd, the inventive method realizes that simply, model solution speed is fast, with very strong practicality.
Embodiment
The active distribution network robust voltage control method that active reactive proposed by the present invention is coordinated, comprises the following steps:
(1) the generating active power interval for presetting the photovoltaic plant in an access active distribution network isAnd Active distribution network control centre will be sent on default active power interval, active distribution network control centre is by photovoltaic plant Active power permission it is intervalPhotovoltaic plant is issued to, the active power of photovoltaic plant is allowing intervalIt is interior with maximum power tracing mode operation;When the active power of photovoltaic plant is more than the interval upper limit of above-mentioned permissionWhen, the active power for making photovoltaic plant is the upper limit
(2) the voltage-controlled object function of active distribution network robust for setting up active reactive coordination is as follows:
Wherein, M abandons light penalty coefficient for photovoltaic plant, and M span is 4~6,P pvWithRespectively photovoltaic electric The lower and upper limit of the active power of active distribution network control centre are sent on standing,WithRespectively active distribution network is dispatched Center is issued to the lower and upper limit of the active power of photovoltaic plant,For each node active power composition in active distribution network Column vector, footmark u represents to make active distribution network magnitude of voltage to be in maximum, footmarkFor the ABC three-phases of active distribution network, For the column vector of each node active power composition in active distribution network, footmark l represents to make active distribution network magnitude of voltage be in minimum Value, sum () represents to sum to the element in column vector;
(3) set active distribution network voltage be in maximum when trend constraint as:
Wherein, footmarkFor the ABC three-phases of active distribution network,WithThe active respectively formed by topology search is matched somebody with somebody The three-phase node branch road incidence matrix and branch impedance matrix of power network,ForIntermediate value be -1 element constitute withDimension phase Same end branch incidence matrix,ForIntermediate value be 1 element constitute withDimension identical top branch road incidence matrix, Real () and imag () represent to take real and imaginary part respectively, × and × it is respectively the sum of products dot product of matrix,WithThe column vector of the column vector of each node active power composition and reactive power composition, angle respectively in active distribution network Mark u represents to make active distribution network magnitude of voltage in maximum,
Wherein,The column vector constituted for the active power lower limit of each load bus of active distribution network,Actively to match somebody with somebody Each photovoltaic plant allows the active power upper limit in power networkThe column vector of composition,For each load bus of active distribution network The column vector of reactive power lower limit composition, QinFor by the reactive power of all photovoltaic plants in active distribution networkStatic nothing The reactive power of work(compensatorWith the reactive power of subdivided capacitorThe column vector of composition;
For the column vector of active distribution network voltage value three phases active power composition of branch road when maximum:
Wherein, N is the circuitry number in active distribution network,WithRepresent active distribution network magnitude of voltage most respectively When big on branch road 1 A, B and C three-phase active power, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value three phase reactive power composition of branch road when maximum:
Wherein, N is the circuitry number in active distribution network,Active distribution network magnitude of voltage is represented respectively When maximum on branch road 1 A, B and C three-phase reactive power, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value electric current modulus value square composition of branch road when maximum:
Wherein, N is the circuitry number in active distribution network,WithRepresent that active distribution network magnitude of voltage is maximum respectively Square of A, B and C three-phase current modulus value on Shi Zhilu 1, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value three-phase voltage modulus value square composition of node when maximum:
Wherein, n is the nodes in active distribution network,Represent that active distribution network magnitude of voltage is maximum respectively Square of Shi Jiedian 1 A, B and C three-phase voltage modulus value, subscript T represents the transposition computing of vector;
(4) set active distribution network voltage be in minimum value when trend constraint as:
Wherein, footmarkFor the ABC three-phases of active distribution network,WithThe active respectively formed by topology search is matched somebody with somebody The three-phase node branch road incidence matrix and branch impedance matrix of power network,ForIntermediate value be -1 element constitute withDimension phase Same end branch incidence matrix,ForIntermediate value be 1 element constitute withDimension identical top branch road incidence matrix, Real () and imag () represent to take real and imaginary part respectively, × and × it is respectively the sum of products dot product of matrix,WithThe row arrow of the column vector of each node active power composition and each node reactive power composition respectively in active distribution network Amount, footmark l represents to make active distribution network magnitude of voltage in minimum value,
Wherein,The column vector constituted for the active power upper limit of each load bus of active distribution network,For active distribution Each photovoltaic plant allows active power lower limit in netThe column vector of composition,For the nothing of each load bus of active distribution network The column vector of work(upper limit of the power composition, QinFor by the reactive power of all photovoltaic plants in active distribution networkStatic reacance The reactive power of compensatorWith the reactive power of subdivided capacitorThe column vector of composition;
For the column vector of active distribution network voltage value three phases active power composition of branch road when minimum:
Wherein, N is the circuitry number in active distribution network,WithRepresent that active distribution network magnitude of voltage is minimum respectively The active power of A, B and C three-phase on Shi Zhilu 1, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value three phase reactive power composition of branch road when minimum:
Wherein, N is the circuitry number in active distribution network,Active distribution network magnitude of voltage is represented respectively When minimum on branch road 1 A, B and C three-phase reactive power, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value electric current modulus value square composition of branch road when minimum:
Wherein, N is the circuitry number in active distribution network,WithRepresent that active distribution network magnitude of voltage is minimum respectively Square of A, B and C three-phase current modulus value on Shi Zhilu 1, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value three-phase voltage modulus value square composition of node when minimum:
Wherein, n is the nodes in active distribution network,Represent that active distribution network magnitude of voltage is minimum respectively Square of Shi Jiedian 1 A, B and C three-phase voltage modulus value, subscript T represents the transposition computing of vector;
(5) active power of photovoltaic plant is constrained in setting active distribution network:
Wherein,P pvWithRespectively be sent on photovoltaic plant active distribution network control centre active power lower limit and The upper limit,WithRespectively active distribution network control centre is issued to the lower and upper limit of the active power of photovoltaic plant;
(6) voltage security of setting active distribution network interior joint is constrained to:
Wherein,WithThe three-phase voltage modulus value of root node is put down when respectively active distribution network voltage is maximum The three-phase voltage modulus value square of root node when side and minimum active distribution network voltage,It is maximum for active distribution network voltage When other nodes in addition to root node three-phase voltage modulus value square,Root section is removed during for active distribution network voltage minimum The three-phase voltage modulus value square of other nodes outside point, UmaxFor the upper safety limit of active distribution network interior joint three-phase voltage, Umin For the lower safety limit of active distribution network interior joint three-phase voltage;
(7) active power and reactive power between setting power transmission network and active distribution network on transmission line are constrained to:
Wherein,WithActive distribution network root node three when respectively active distribution network magnitude of voltage is maximum The three-dimensional column vector of phase active power composition and the three-dimensional column vector of three phase reactive power composition,WithRespectively For active distribution network voltage value active distribution network root node three phases active power composition when minimum three-dimensional column vector and three-phase without The three-dimensional column vector of work(power composition,WithActive power on transmission line respectively between power transmission network and power distribution network Lower safety limit and the upper limit,WithRespectively between power transmission network and power distribution network on transmission line reactive power lower safety limit And the upper limit;
(8) reactive capability of photovoltaic plant and SVC is constrained in setting active distribution network:
Wherein,WithThe reactive power of photovoltaic plant and SVC respectively in active distribution network,The respectively reactive power maximum size and lower limit of photovoltaic plant,Respectively static reacance The reactive power maximum size and lower limit of compensator;
(9) reactive power of subdivided capacitor is constrained in setting active distribution network:
T=α0*201*21
Wherein,For the reactive power of subdivided capacitor,For the idle work(in subdivided capacitor per group capacitor Rate, the group number that t puts into operation for subdivided capacitor in active distribution network, α0、α1Respectively variable, value is 0 or 1;
(10) Primal-dual Interior Point is used, according to above-mentioned steps (3), the trend constraint of (4) and step (5)~step Suddenly the operation constraint of the active distribution network of (9), the object function of solution procedure (2), calculating obtains photovoltaic electric in active distribution network Reactive power value, the reactive power value of SVC, the reactive power value of subdivided capacitor and the photovoltaic plant stood have The permission of work(power is interval, realizes that active distribution network robust voltage is controlled.

Claims (1)

1. the active distribution network robust voltage control method that a kind of active reactive is coordinated, it is characterised in that this method includes following step Suddenly:
(1) the generating active power interval for presetting the photovoltaic plant in an access active distribution network isAnd this is pre- If active power interval on be sent to active distribution network control centre, active distribution network control centre is by the wattful power of photovoltaic plant The permission of rate is intervalPhotovoltaic plant is issued to, when the active power of photovoltaic plant is more than the interval upper limit of above-mentioned permissionWhen, the active power for making photovoltaic plant is the upper limit
(2) the voltage-controlled object function of active distribution network robust for setting up active reactive coordination is as follows:
Wherein, M abandons light penalty coefficient for photovoltaic plant, and M span is 4~6,WithSent out respectively on photovoltaic plant To the lower and upper limit of the active power of active distribution network control centre,WithRespectively under active distribution network control centre The lower and upper limit of the active power of photovoltaic plant are sent to,For the row arrow of each node active power composition in active distribution network Amount, footmark u represents to make active distribution network magnitude of voltage be in maximum, footmarkFor the ABC three-phases of active distribution network,For actively The column vector of each node active power composition in power distribution network, footmark l represents to make active distribution network magnitude of voltage be in minimum value, sum () represents to sum to the element in column vector;
(3) set active distribution network voltage be in maximum when trend constraint as:
Wherein, footmarkFor the ABC three-phases of active distribution network,WithThe active distribution network respectively formed by topology search Three-phase node branch road incidence matrix and branch impedance matrix,ForIntermediate value be -1 element constitute withDimension identical End branch incidence matrix,ForIntermediate value be 1 element constitute withDimension identical top branch road incidence matrix, real () and imag () represent to take real and imaginary part respectively, × and × it is respectively the sum of products dot product of matrix,WithThe column vector of the column vector of each node active power composition and reactive power composition, footmark u tables respectively in active distribution network Showing makes active distribution network magnitude of voltage be in maximum,
Wherein,The column vector constituted for the active power lower limit of each load bus of active distribution network,For in active distribution network Each photovoltaic plant allows the active power upper limitThe column vector of composition,For the idle work(of each load bus of active distribution network The column vector of rate lower limit composition, QinFor by the reactive power of all photovoltaic plants in active distribution networkStatic reactive The reactive power of deviceWith the reactive power of subdivided capacitorThe column vector of composition;
For the column vector of active distribution network voltage value three phases active power composition of branch road when maximum:
Wherein, N is the circuitry number in active distribution network,WithBranch when active distribution network magnitude of voltage is maximum is represented respectively The active power of A, B and C three-phase on road 1, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value three phase reactive power composition of branch road when maximum:
Wherein, N is the circuitry number in active distribution network,Represent that active distribution network magnitude of voltage is maximum respectively The reactive power of A, B and C three-phase on Shi Zhilu 1, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value electric current modulus value square composition of branch road when maximum:
Wherein, N is the circuitry number in active distribution network,WithBranch when active distribution network magnitude of voltage is maximum is represented respectively Square of A, B and C three-phase current modulus value on road 1, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value three-phase voltage modulus value square composition of node when maximum:
Wherein, n is the nodes in active distribution network,The active distribution network magnitude of voltage maximum time is represented respectively Square of A, B and C three-phase voltage modulus value of point 1, subscript T represents the transposition computing of vector;
(4) set active distribution network voltage be in minimum value when trend constraint as:
Wherein, footmarkFor the ABC three-phases of active distribution network,WithThe active distribution network respectively formed by topology search Three-phase node branch road incidence matrix and branch impedance matrix,ForIntermediate value be -1 element constitute withDimension identical End branch incidence matrix,ForIntermediate value be 1 element constitute withDimension identical top branch road incidence matrix, real () and imag () represent to take real and imaginary part respectively, × and × it is respectively the sum of products dot product of matrix,WithThe column vector of the column vector of each node active power composition and each node reactive power composition respectively in active distribution network, Footmark l represents to make active distribution network magnitude of voltage in minimum value,
Wherein,The column vector constituted for the active power upper limit of each load bus of active distribution network,For in active distribution network Each photovoltaic plant allows active power lower limitThe column vector of composition,For the idle work(of each load bus of active distribution network The column vector of rate upper limit composition, QinFor by the reactive power of all photovoltaic plants in active distribution networkStatic reactive The reactive power of deviceWith the reactive power of subdivided capacitorThe column vector of composition;
For the column vector of active distribution network voltage value three phases active power composition of branch road when minimum:
Wherein, N is the circuitry number in active distribution network,WithBranch when active distribution network magnitude of voltage is minimum is represented respectively The active power of A, B and C three-phase on road 1, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value three phase reactive power composition of branch road when minimum:
Wherein, N is the circuitry number in active distribution network,Represent that active distribution network magnitude of voltage is minimum respectively The reactive power of A, B and C three-phase on Shi Zhilu 1, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value electric current modulus value square composition of branch road when minimum:
Wherein, N is the circuitry number in active distribution network,WithBranch when active distribution network magnitude of voltage is minimum is represented respectively Square of A, B and C three-phase current modulus value on road 1, subscript T represents the transposition computing of vector;
For the column vector of active distribution network voltage value three-phase voltage modulus value square composition of node when minimum:
Wherein, n is the nodes in active distribution network,The active distribution network magnitude of voltage minimum time is represented respectively Square of A, B and C three-phase voltage modulus value of point 1, subscript T represents the transposition computing of vector;
(5) active power of photovoltaic plant is constrained in setting active distribution network:
P ‾ ^ pv ≤ P ‾ pv
P ^ ‾ pv ≤ P ‾ pv
Wherein,WithThe lower and upper limit of the active power of active distribution network control centre are respectively sent on photovoltaic plant,WithRespectively active distribution network control centre is issued to the lower and upper limit of the active power of photovoltaic plant;
(6) voltage security of setting active distribution network interior joint is constrained to:
Wherein,WithThe three-phase voltage modulus value quadratic sum of root node when respectively active distribution network voltage is maximum The three-phase voltage modulus value square of root node when active distribution network voltage is minimum,Removed during for active distribution network voltage maximum The three-phase voltage modulus value square of other nodes outside root node,During for active distribution network voltage minimum in addition to root node Other nodes three-phase voltage modulus value square, UmaxFor the upper safety limit of active distribution network interior joint three-phase voltage, UminBased on The lower safety limit of dynamic power distribution network interior joint three-phase voltage;
(7) active power and reactive power between setting power transmission network and active distribution network on transmission line are constrained to:
Wherein,WithActive distribution network root node three-phase is active when respectively active distribution network magnitude of voltage is maximum The three-dimensional column vector of power composition and the three-dimensional column vector of three phase reactive power composition,WithRespectively actively The three-dimensional column vector and three phase reactive power of active distribution network root node three phases active power composition when distribution network voltage value is minimum The three-dimensional column vector of composition,WithRespectively between power transmission network and power distribution network on transmission line under the safety of active power Limit and the upper limit,WithRespectively between power transmission network and power distribution network on transmission line reactive power lower safety limit and the upper limit;
(8) reactive capability of photovoltaic plant and SVC is constrained in setting active distribution network:
Wherein,WithThe reactive power of photovoltaic plant and SVC respectively in active distribution network,The respectively reactive power maximum size and lower limit of photovoltaic plant,Respectively static reacance The reactive power maximum size and lower limit of compensator;
(9) reactive power of subdivided capacitor is constrained in setting active distribution network:
Q in CB = t * Q in CB , step
T=α0*201*21
Wherein,For the reactive power of subdivided capacitor,For the reactive power in subdivided capacitor per group capacitor, t The group number put into operation for subdivided capacitor in active distribution network, α0、α1Respectively variable, value is 0 or 1;
(10) Primal-dual Interior Point is used, according to above-mentioned steps (3), the trend constraint of (4) and step (5)~step (9) Active distribution network operation constraint, the object function of solution procedure (2) calculates and obtains the nothing of photovoltaic plant in active distribution network Work(performance number, the reactive power value of SVC, the reactive power value of subdivided capacitor and photovoltaic plant active power Permission it is interval, realize that active distribution network robust voltage is controlled.
CN201510226565.6A 2015-05-06 2015-05-06 The active distribution network robust voltage control method that a kind of active reactive is coordinated Active CN104868480B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510226565.6A CN104868480B (en) 2015-05-06 2015-05-06 The active distribution network robust voltage control method that a kind of active reactive is coordinated

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510226565.6A CN104868480B (en) 2015-05-06 2015-05-06 The active distribution network robust voltage control method that a kind of active reactive is coordinated

Publications (2)

Publication Number Publication Date
CN104868480A CN104868480A (en) 2015-08-26
CN104868480B true CN104868480B (en) 2017-07-14

Family

ID=53914113

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510226565.6A Active CN104868480B (en) 2015-05-06 2015-05-06 The active distribution network robust voltage control method that a kind of active reactive is coordinated

Country Status (1)

Country Link
CN (1) CN104868480B (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105322534B (en) * 2015-10-08 2016-11-09 南京邮电大学 A kind of microgrid Optimization Scheduling based on bounded-but-unknown uncertainty
CN106655148B (en) * 2015-10-29 2019-02-05 中国电力科学研究院 A kind of generation of electricity by new energy base power grid static security risk online evaluation method
CN105337289B (en) * 2015-11-19 2018-07-24 河海大学 A kind of area power grid power-less optimized controlling method based on incidence matrix
CN105607480A (en) * 2016-01-22 2016-05-25 梅照付 Gateway flow control system for intelligent power grid
CN105607479A (en) * 2016-01-22 2016-05-25 梅照付 Gateway flow control method for intelligent power grid
CN105488357B (en) * 2016-01-26 2018-02-02 清华大学 A kind of active power rolling scheduling method of photo-thermal power station wind power plant association system
CN105870962B (en) * 2016-03-31 2019-03-12 国家电网公司 A kind of Robust Interval wind-powered electricity generation dispatching method considering power system frequency response characteristic
CN106159974B (en) * 2016-08-02 2019-01-15 清华大学 A kind of distributed reactive Voltage Optimum method that transmission & distribution are coordinated
CN108092321B (en) * 2016-11-22 2023-02-07 中国电力科学研究院 Active power and reactive power coordinated control method considering uncertainty for active power distribution network
CN108334696B (en) * 2018-02-02 2021-06-08 国网江西省电力有限公司经济技术研究院 Day-ahead network reconstruction method for power distribution network considering power randomness
CN108539757B (en) * 2018-05-09 2020-05-12 北京交通大学 Reactive power scheduling method for optical storage cluster of power distribution network
CN113328435B (en) * 2021-05-26 2022-12-09 国网河北省电力有限公司 Active and reactive power combined control method for active power distribution network based on reinforcement learning

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11289661A (en) * 1998-04-06 1999-10-19 Kansai Electric Power Co Inc:The Power flow-calculating device for power distribution system control system
CN103701132A (en) * 2013-12-20 2014-04-02 清华大学 Method for controlling overvoltage of active power distribution network on basis of active-reactive coordination
CN103701142A (en) * 2014-01-10 2014-04-02 贵州电网公司信息通信分公司 Active power distributing network real-time reactive power-voltage control method considering discrete control variable

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11289661A (en) * 1998-04-06 1999-10-19 Kansai Electric Power Co Inc:The Power flow-calculating device for power distribution system control system
CN103701132A (en) * 2013-12-20 2014-04-02 清华大学 Method for controlling overvoltage of active power distribution network on basis of active-reactive coordination
CN103701142A (en) * 2014-01-10 2014-04-02 贵州电网公司信息通信分公司 Active power distributing network real-time reactive power-voltage control method considering discrete control variable

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A Robust Wind Power Optimization Method for Look-Ahead Power Dispatch;Wenchuan Wu et al.;《IEEE Transactions on Sustainable Energy》;20140430;第5卷(第2期);第507-515页 *
Robust Nonlinear Controller Design for Three-Phase Grid-Connected Photovoltaic Systems Under Structured Uncertainties;M. A. Mahmud et al.;《IEEE Transactions on Power Delivery》;20140630;第29卷(第3期);第1221-1230页 *
基于混合整数二阶锥规划的三相有源配电网无功优化;刘一兵等;《电力系统自动化》;20140810;第38卷(第15期);第58-64页 *
自动电压控制中连续变量与离散变量的协调方法(一)变电站内协调电压控制;郭庆来等;《电力系统自动化》;20080425;第32卷(第8期);第39-42、107页 *

Also Published As

Publication number Publication date
CN104868480A (en) 2015-08-26

Similar Documents

Publication Publication Date Title
CN104868480B (en) The active distribution network robust voltage control method that a kind of active reactive is coordinated
CN105119280B (en) Alternating current-direct current mixed structure active power distribution network running optimizatin method based on cone optimization
CN104533714B (en) A kind of wind turbine generator system controller parameter identification and setting method
CN110690732B (en) Photovoltaic reactive power partition pricing power distribution network reactive power optimization method
CN107679769B (en) Method for establishing frequency response model of power system containing wind power and method for calculating frequency characteristic index
CN106410835A (en) Power grid simulation realization method of multi-terminal flexible direct-current transmission system
Yun et al. Online Thevenin equivalent parameter identification method of large power grids using LU factorization
CN105870968B (en) A kind of three-phase imbalance reactive voltage control method of meter and system negative sequence voltage
CN103701132B (en) A kind of superpotential method of control active distribution network based on meritorious-idle coordination
CN105322535A (en) Two-stage optimal power flow calculation method for power supply containing unified power flow controller
CN106026113A (en) Micro-grid system monitoring method having reactive automatic compensation function
CN101882785A (en) Flow analysis system and method for intelligent power distribution networks containing distributed generation
CN104779609B (en) A kind of trend cooperative control method for interconnected network
CN106779444B (en) The active plan load flow rectification method and apparatus extended out based on electric network model
CN107332252A (en) A kind of power distribution network low-voltage comprehensive processing method of meter and broad sense reactive source
CN105162129B (en) Meter and the distribution reactive voltage control method of distributed power source allocation optimum
CN102769292B (en) Method for enhancing voltage safety margin of wind farm
CN104993525A (en) Active power distribution network coordination optimization control method considering ZIP loads
CN107732926A (en) A kind of AC-DC hybrid power grid idle work optimization method containing voltage source converter
CN107039981A (en) One kind intends direct current linearisation probability optimal load flow computational methods
CN103715687B (en) The congested real-time control method of a kind of active distribution network branch power
CN107196343A (en) A kind of voltage scheduling method a few days ago of multiterminal flexible direct current island-grid sending end
CN103606949A (en) Control method for improving voltage stability of power grid of large-scale photovoltaic access area
CN105958530A (en) Microgrid system with reactive power automatic compensation function
CN109149566A (en) A kind of modeling method of the simulation model of the high-power minimum point prediction of missing lower frequency

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant