CN106099946A - The collocation method of electrical network dynamic reactive capacity and system - Google Patents

The collocation method of electrical network dynamic reactive capacity and system Download PDF

Info

Publication number
CN106099946A
CN106099946A CN201610615753.2A CN201610615753A CN106099946A CN 106099946 A CN106099946 A CN 106099946A CN 201610615753 A CN201610615753 A CN 201610615753A CN 106099946 A CN106099946 A CN 106099946A
Authority
CN
China
Prior art keywords
dynamic reactive
reactive capacity
destination node
configuration mode
risk factor
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
CN201610615753.2A
Other languages
Chinese (zh)
Other versions
CN106099946B (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.)
Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd
Original Assignee
Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd
Beijing King Star Hi Tech System Control 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 Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd, Beijing King Star Hi Tech System Control Co Ltd filed Critical Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd
Priority to CN201610615753.2A priority Critical patent/CN106099946B/en
Publication of CN106099946A publication Critical patent/CN106099946A/en
Application granted granted Critical
Publication of CN106099946B publication Critical patent/CN106099946B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/18Arrangements for adjusting, eliminating or compensating reactive power in networks
    • 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 collocation method and the system of a kind of electrical network dynamic reactive capacity, described method includes: according to each destination node voltage after fault occurs and default normal voltage range, calculate the risk factor that multiple dynamic reactive capacity configuration mode is respectively corresponding;The dynamic reactive power of each destination node that multiple dynamic reactive capacity configuration mode described in matching is corresponding and the risk factor of correspondence, obtain relationship by objective (RBO) formula;According to the first default optimal conditions, described relationship by objective (RBO) formula is optimized calculating, obtains the optimal solution of the dynamic reactive capacity of each destination node;The dynamic reactive capacity of each destination node described in the optimal solution corresponding configuration of the dynamic reactive capacity according to described each destination node.The collocation method of above-mentioned electrical network dynamic reactive capacity and system, the overall situation considers multiple dynamic reactive capacity configuration mode and tackles the risk factor of various fault, the various faults being likely to occur network system have general applicability, can be effectively improved stability and the reliability of electrical network.

Description

The collocation method of electrical network dynamic reactive capacity and system
Technical field
The present invention relates to distribution technique field, particularly relate to the collocation method of a kind of electrical network dynamic reactive capacity and be System.
Background technology
In the operation and scheduling of electrical network, implement Reactive Power Control, keep certain reactive reserve to be to ensure that electrical network supplies Electricity quality, prevent electrical network from the necessary means of collapse of voltage accident occurring after being impacted by big disturbance or fault.Therefore, at electricity During network operation, need some nodes in network system are configured reactive capability.
In prior art, the configuration of the reactive capability of electrical network is static, i.e. mainly for electrical network certain time point therefore Barrier situation determines a kind of reactive capability configuration mode, and according to the reactive capability configuration mode determined to configurable in network system The node of reactive capability carries out reactive capability configuration.But, in practical situation, disturbance or fault that electrical network is subject to have undulatory property And unpredictability, the reactive capability collocation method of this static state is difficult to persistently protect when being impacted by different faults at electrical network The stable operation of card electrical network.
Summary of the invention
Based on this, it is necessary to be difficult to when electrical network is by different disturbances or fault hold for existing reactive capability configuration mode The defect of continuation of insurance card power grid operation, it is provided that the collocation method of a kind of electrical network dynamic reactive capacity and system.
Embodiment of the present invention first aspect provides the collocation method of a kind of electrical network dynamic reactive capacity, comprising:
According to each destination node fault occur after voltage and default normal voltage range, calculate multiple dynamic reactive The risk factor that capacity configuration mode is the most corresponding;
The dynamic reactive power of each destination node that multiple dynamic reactive capacity configuration mode described in matching is corresponding and correspondence Risk factor, obtain relationship by objective (RBO) formula;
According to the first default optimal conditions, described relationship by objective (RBO) formula is optimized calculating, obtains the dynamic of each destination node The optimal solution of state reactive capability;
Described in the optimal solution corresponding configuration of the dynamic reactive capacity according to described each destination node, each destination node is dynamic Reactive capability.
In one embodiment, each destination node that multiple dynamic reactive capacity configuration mode described in described matching is corresponding Dynamic reactive power and the risk factor of correspondence, obtain relationship by objective (RBO) formula, including:
According to gaussian kernel function and the second optimal conditions calculation optimization multiplier based on described risk factor, wherein said height The variable of this kernel function is included in the dynamic reactive power of each destination node in the case of described multiple dynamic reactive capacity configuration;
Using described optimization multiplier as the coefficient of described gaussian kernel function, it is thus achieved that described relationship by objective (RBO) formula.
In one embodiment, described gaussian kernel function isDescribed second excellent Change condition isDescribed optimization multiplier is α=(α1, α2..., αk);Described relationship by objective (RBO) formula is
y k = f ( x k ) = K ( x k , x k ) - 2 Σ p = 1 10000 α p K ( x k , x p ) + Σ q = 1 10000 Σ p = 1 10000 α p α q K ( x p , x q ) ;
Wherein, N1Quantity for described multiple dynamic reactive capacity configuration mode;Described D represents the element that pth row q arranges For Dpq=K (xp, xq)ypyqMatrix;Described c=(-1 ... ,-1)T, expression line number is N1, columns is the matrix of 1;Described Aα =(1,1 ..., 1), representing that line number is 1, columns is N1Matrix;DescribedykRepresent that kth kind is moved The risk factor that state reactive capability configuration mode is corresponding; Represent n-th destination node dynamic reactive merit after fault occurs under kth kind dynamic reactive capacity configuration mode The maximum of rate;K=1,2 ..., N1, n=1,2 ..., Nsvc, NsvcQuantity for described destination node.
In one embodiment, described first optimal conditions isWherein QsvcTotal dynamic reactive capacity for network system.
In one embodiment, described according to each destination node fault occur after voltage and default normal voltage model Enclose, calculate the risk factor that multiple dynamic reactive capacity configuration mode is the most corresponding, including:
According to default constraints, accounting equationObtain The risk factor that described multiple dynamic reactive capacity configuration mode is corresponding
Wherein, described default constraints includes system load flow constraints, each busbar voltage constraints, transmission line The constraints such as or not phase angle difference constraints and control variable;
Described risk (k) is the risk factor that kth kind dynamic reactive capacity configuration mode is corresponding, described NsvcFor described mesh The quantity of mark node;DescribedFor under kth kind dynamic reactive capacity configuration mode, the n-th destination node is sent out the s fault Voltage after life;DescribedIt it is the upper voltage limit in the normal voltage range of the n-th destination node;DescribedV n Represent the n-th mesh Lower voltage limit in the normal voltage range of mark node;DescribedWith In characterizing kth kind dynamic reactive capacity configuration mode, whereinRepresent under kth kind dynamic reactive capacity configuration mode n-th The maximum of individual destination node dynamic reactive power after fault occurs, n=1,2 ..., Nsvc, NsvcFor described destination node Quantity.
Embodiment of the present invention second aspect provides the configuration system of a kind of electrical network dynamic reactive capacity, comprising:
First computing module, for according to each destination node fault occur after voltage and default normal voltage model Enclose, calculate the risk factor that multiple dynamic reactive capacity configuration mode is the most corresponding;
Fitting module, for the dynamic nothing of each destination node corresponding to dynamic reactive capacity configuration mode multiple described in matching Merit power and the risk factor of correspondence, obtain relationship by objective (RBO) formula;
Second computing module, for described relationship by objective (RBO) formula being optimized calculating according to the first optimal conditions preset, Obtain the optimal solution of the dynamic reactive capacity of each destination node;And
Configuration module, for according to each mesh described in the optimal solution corresponding configuration of the dynamic reactive capacity of described each destination node The dynamic reactive capacity of mark node.
In one embodiment, described fitting module includes:
Computing unit, for taking advantage of according to gaussian kernel function and the second optimal conditions calculation optimization based on described risk factor Son, the variable of wherein said gaussian kernel function is included in the dynamic of each destination node in the case of described multiple dynamic reactive capacity configuration State reactive power;And
Obtain unit, for the coefficient using described optimization multiplier as described gaussian kernel function, it is thus achieved that described relationship by objective (RBO) Formula.
In one embodiment, described gaussian kernel function isDescribed second excellent Change condition isDescribed optimization multiplier is α=(α1, α2..., αk);Described relationship by objective (RBO) formula is
y k = f ( x k ) = K ( x k , x k ) - 2 Σ p = 1 10000 α p K ( x k , x p ) + Σ q = 1 10000 Σ p = 1 10000 α p α q K ( x p , x q ) ;
Wherein, N1Quantity for described multiple dynamic reactive capacity configuration mode;Described D represents the element that pth row q arranges For Dpq=K (xp, xq)ypyqMatrix;Described c=(-1 ... ,-1)T, expression line number is N1, columns is the matrix of 1;Described Aα =(1,1 ..., 1), representing that line number is 1, columns is N1Matrix;DescribedykRepresent that kth kind is moved The risk factor that state reactive capability configuration mode is corresponding; Represent n-th destination node dynamic reactive merit after fault occurs under kth kind dynamic reactive capacity configuration mode The maximum of rate;K=1,2 ..., N1, n=1,2 ..., Nsvc, NsvcQuantity for described destination node.
In one embodiment, described first optimal conditions isWherein QsvcTotal dynamic reactive capacity for network system.
In one embodiment, described first computing module is used for:
According to default constraints, accounting equationObtain institute State the risk factor that multiple dynamic reactive capacity configuration mode is corresponding
Wherein, described default constraints includes system load flow constraints, each busbar voltage constraints, transmission line The constraints such as or not phase angle difference constraints and control variable;
Described risk (k) is the risk factor that kth kind dynamic reactive capacity configuration mode is corresponding, described NsvcFor described mesh The quantity of mark node;DescribedFor under kth kind dynamic reactive capacity configuration mode, the n-th destination node is sent out the s fault Voltage after life;DescribedIt it is the upper voltage limit in the normal voltage range of the n-th destination node;DescribedV n Represent the n-th mesh Lower voltage limit in the normal voltage range of mark node;DescribedWith In characterizing kth kind dynamic reactive capacity configuration mode, whereinRepresent under kth kind dynamic reactive capacity configuration mode n-th The maximum of individual destination node dynamic reactive power after fault occurs, n=1,2 ..., Nsvc, NsvcFor described destination node Quantity.
The collocation method of above-mentioned electrical network dynamic reactive capacity and system, corresponding according to multiple dynamic reactive capacity configuration mode The dynamic reactive power of each destination node and risk factor when tackling various fault, fit object relational expression also asks optimum Solve, owing to this optimal solution is that the overall situation considers the risk factor that multiple dynamic reactive capacity configuration mode tackles various fault and obtains , the various faults being likely to occur network system have general applicability, therefore according to this optimal solution to each in network system The dynamic reactive capacity of destination node configures, and can at utmost maintain electrical network electricity when electrical network is by different disturbances or fault Stablizing of pressure, the fluctuation of suppression line voltage, thus it is effectively improved stability and the reliability of electrical network.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of the collocation method of the electrical network dynamic reactive capacity of one embodiment of the invention;
Fig. 2 is the schematic flow sheet of the collocation method of the electrical network dynamic reactive capacity of another embodiment of the present invention;
Fig. 3 is the modular structure schematic diagram of the configuration system of the electrical network dynamic reactive capacity of one embodiment of the invention;And
Fig. 4 is the modular structure schematic diagram of the fitting module of one embodiment of the invention.
Detailed description of the invention
Understandable, below in conjunction with the accompanying drawings to the present invention for enabling the above-mentioned purpose of the present invention, feature and advantage to become apparent from Detailed description of the invention be described in detail.Elaborate a lot of detail in the following description so that fully understanding this Bright.But the present invention can implement to be much different from alternate manner described here, and those skilled in the art can be not Doing similar improvement in the case of running counter to intension of the present invention, therefore the present invention is not limited by following public specific embodiment.
In describing the invention, it is to be understood that term " first ", " second " are only used for describing purpose, and can not It is interpreted as instruction or hint relative importance or the implicit quantity indicating indicated technical characteristic.Thus, define " the One ", the feature of " second " can express or implicitly include at least one this feature.In describing the invention, " multiple " It is meant that at least two, such as two, three etc., unless otherwise expressly limited specifically.
Refer to the flow process signal of the collocation method of the electrical network dynamic reactive capacity that Fig. 1, Fig. 1 are one embodiment of the invention Figure.As it is shown in figure 1, the collocation method of this electrical network dynamic reactive capacity comprises the steps that
Step 110, according to each destination node fault occur after voltage and default normal voltage range, calculate multiple The risk factor that dynamic reactive capacity configuration mode is the most corresponding.
In the present embodiment, destination node is the node that can be configured dynamic reactive capacity, such as, include in network system Reactive power source.
Specifically, the normal voltage range of each destination node in network system can be pre-set, each in monitoring network system Destination node voltage before and after various faults occur.Under various dynamic reactive capacity configuration modes, by each target is saved The normal voltage range putting the voltage after every kind of fault occurs corresponding compares, and can calculate every kind of dynamic reactive capacity The risk factor that configuration mode is corresponding.
Above-mentioned risk factor can characterize dynamic reactive capacity configuration mode and maintain when tackling various fault line voltage steady Fixed, the ability of suppression voltage ripple of power network.Risk factor is the biggest, represents that dynamic reactive capacity configuration mode suppresses line voltage ripple Dynamic ability is the most weak.Wherein, the most repeatedly data of fault, if certain dynamic reactive capacity configuration mode is after fault occurs, The voltage of each destination node more deviates normal voltage range, then the risk factor that this kind of dynamic reactive capacity configuration mode is corresponding is more Greatly.
Step 130, the dynamic reactive merit of each destination node that multiple dynamic reactive capacity configuration mode described in matching is corresponding Rate and the risk factor of correspondence, obtain relationship by objective (RBO) formula.
Wherein, the dynamic reactive power of each destination node that multiple dynamic reactive capacity configuration mode is corresponding and step 110 The risk factor that calculated multiple dynamic reactive capacity configuration mode is the most corresponding, organizes discrete numerical value, by right for more These discrete numerical value are fitted, and can obtain between expression trend reactive capability configuration mode and risk factor is corresponding The relationship by objective (RBO) formula of relation.
In being embodied as, suitable core can be selected based on SVC (Support Vector Machines, support vector machine) Function and the dynamic reactive power of the optimal conditions each destination node corresponding to above-mentioned multiple dynamic reactive capacity configuration mode and Corresponding risk factor is fitted.For example, it is possible to save with each target that above-mentioned multiple dynamic reactive capacity configuration mode is corresponding The dynamic reactive power of point is as the variable of gaussian kernel function, with risk corresponding to above-mentioned multiple dynamic reactive capacity configuration mode Coefficient is fitted as the variable of risk optimization condition.
Step 150, is optimized calculating according to the first default optimal conditions to described relationship by objective (RBO) formula, obtains each target The optimal solution of the dynamic reactive capacity of node.
In being embodied as, above-mentioned first optimal conditions can be joining based on above-mentioned multiple dynamic reactive capacity of pre-setting The optimal conditions of the risk factor that mode of putting is corresponding.Such as, the first optimal conditions can be that above-mentioned multiple dynamic reactive capacity is joined The risk factor sum minimum that mode of putting is corresponding.Such as, the first optimal conditions can be to be multiple dynamic reactive capacity configuration mode The meansigma methods of corresponding risk factor is minimum.Such as, the first optimal conditions can be above-mentioned multiple dynamic reactive capacity configuration side The variance of the risk factor that formula is corresponding is minimum.
By seeking above-mentioned relationship by objective (RBO) formula optimal solution under the first optimal conditions, the dynamic nothing of available each destination node The optimal solution of merit capacity configuration, is designated asWhereinTable Show the optimal solution of the dynamic reactive capacity of destination node 1,Represent the optimal solution of the dynamic reactive capacity of destination node 2,Represent destination node NsvcThe optimal solution of dynamic reactive capacity.
Step 170, saves according to each target described in the optimal solution corresponding configuration of the dynamic reactive capacity of described each destination node The dynamic reactive capacity of point.
After obtaining the optimal solution of dynamic reactive capacity of each destination node, can be according in this optimal solution corresponding configuration network system The dynamic reactive capacity of each destination node.Such as, according to this optimal solution By the dynamic reactive capacity configuration of destination node 1 it isBy the dynamic reactive capacity configuration of destination node 2 it is By the dynamic reactive capacity configuration of destination node n it isBy destination node NsvcDynamic reactive capacity configuration be
The collocation method of above-mentioned electrical network dynamic reactive capacity, according to each mesh that multiple dynamic reactive capacity configuration mode is corresponding The mark dynamic reactive power of node and risk factor when tackling various fault, fit object relational expression also seeks optimal solution, due to This optimal solution is that the overall situation considers the risk factor that multiple dynamic reactive capacity configuration mode tackles various fault and obtains, to electricity The various faults that net system is likely to occur have general applicability, therefore according to this optimal solution to destination node each in network system Dynamic reactive capacity configure, can at utmost maintain the steady of line voltage when electrical network is by different disturbances or fault Fixed, the fluctuation of suppression line voltage, thus it is effectively improved stability and the reliability of electrical network.
In one embodiment, before step 110, can be according to the configurable total dynamic reactive capacity of network system and electricity The quantity of destination node in net system, the dynamic reactive capacity of each destination node of random arrangement, obtain multiple dynamic reactive capacity Configuration mode.Specifically, the dynamic nothing that every kind of dynamic reactive capacity configuration mode can be sent by destination node each under which Merit power represents.
Such as, if with QsvcRepresent the configurable total dynamic reactive capacity of network system, with NsvcRepresent mesh in network system The quantity of mark node, with N1Represent the quantity of above-mentioned multiple dynamic reactive capacity configuration mode, then kth kind dynamic reactive capacity is joined Putting under mode, maximum and the minima of destination node n dynamic reactive power after fault occurs can be respectivelyAndWherein k=1,2 ..., N1, n=1,2 ..., Nsvc.Wherein, quantity N of above-mentioned multiple dynamic reactive capacity configuration mode1Can enter according to the operational capability of system or calculation resources Row is arranged, N1The biggest, i.e. the discrete data for fit object relational expression is the most, and the relationship by objective (RBO) formula finally given the most more has Representativeness, but computing simultaneously is more complicated, and the calculation resources of system consumes the most.Such as, N1May be provided at 8000 to 20000 Between.
In one embodiment, step 110 is particularly as follows: according to default constraints, accounting equationObtain described multiple dynamic reactive capacity configuration mode corresponding Risk factor
Wherein, default constraints includes system load flow constraints, each busbar voltage constraints, transmission line phase angle Difference constraints and the constraints such as or not control variable;Risk (k) is the wind that kth kind dynamic reactive capacity configuration mode is corresponding Danger coefficient, NsvcQuantity for destination node;For under kth kind dynamic reactive capacity configuration mode, the n-th destination node exists Voltage after the s fault generation;It it is the upper voltage limit in the normal voltage range of the n-th destination node;V n Represent n-th Lower voltage limit in the normal voltage range of destination node;For Dynamic reactive capacity configuration mode is planted, wherein after characterizing theRepresent under kth kind dynamic reactive capacity configuration mode n-th The maximum of individual destination node dynamic reactive power after fault occurs, n=1,2 ..., Nsvc, NsvcNumber for destination node Amount.
Specifically, said system trend constraints includes network system trend constraints under non-faulting stateAnd network system nonserviceable under trend constraintsWherein N is the quantity of network system interior joint.WithRepresent respectively under kth kind dynamic reactive capacity configuration mode, node i during electrical network properly functioning (being i.e. in non-faulting state) Meritorious and reactive power;WithRepresent respectively under kth kind dynamic reactive capacity configuration mode, after fault s occurs Meritorious and the reactive power of node i.WithBeing illustrated respectively under kth kind dynamic reactive capacity configuration mode, electrical network is just Node i and the voltage of node j when often running;WithIt is illustrated respectively under kth kind dynamic reactive capacity configuration mode, There is posterior nodal point i and the voltage of node j in fault s.GijAnd BijRepresent the conductance between node i and node j and susceptance respectively.WithRepresent respectively under kth kind dynamic reactive capacity configuration mode, after electrical network is properly functioning and fault s occurs, node i And the phase angle difference between node j.
Specifically, above-mentioned each busbar voltage constraints includesWherein,WithV i Represent joint respectively The upper voltage limit of some i and lower limit,Represent under kth kind dynamic reactive capacity configuration mode, node i when electrical network is properly functioning Voltage.
Specifically, above-mentioned transmission line phase angle difference constraints includesWherein,Withδ ij Respectively The upper voltage limit of the phase angle difference between expression node i and node j and lower limit,Represent in kth kind dynamic reactive capacity configuration Under mode, phase angle difference between node i and node j when electrical network is properly functioning.
Specifically, the constraints such as or not above-mentioned control variable includesWherein,WithQ i Represent the reactive power upper limit and the reactive power lower limit of node i respectively,WithP i Represent the active power of node i respectively The upper limit and active power lower limit,WithBeing illustrated respectively under kth kind dynamic reactive capacity configuration mode, electrical network is normally transported The active power of node i and reactive power during row.WithIt is illustrated respectively in kth kind dynamic reactive capacity configuration mode Under, the reactive power of destination node n when electrical network is properly functioning and after fault s occurs.Represent in the dynamic nothing of kth kind Under power capacity amount configuration mode, the dynamic reactive power value that destination node n sends after fault s occurs.With Being illustrated respectively under kth kind dynamic reactive capacity configuration mode, the dynamic reactive power that destination node n sends after a failure is maximum Value and minima.
In one embodiment, as in figure 2 it is shown, step 130 comprises the steps that
Step 131, according to gaussian kernel function and the second optimal conditions calculation optimization multiplier based on described risk factor, its Described in the variable of gaussian kernel function be included in the dynamic nothing of each destination node in the case of described multiple dynamic reactive capacity configuration Merit power.
Specifically, described gaussian kernel function isDescribed second optimal conditions isCalculated optimization multiplier can be designated as α=(α1, α2..., αk)。
Step 133, using described optimization multiplier as the coefficient of described gaussian kernel function, it is thus achieved that described relationship by objective (RBO) formula.
Specifically, by above-mentioned optimization multiplier α=(α1, α2..., αk) as the coefficient of above-mentioned gaussian kernel function, obtain Relationship by objective (RBO) formula is:
y k = f ( x k ) = K ( x k , x k ) - 2 Σ p = 1 10000 α p K ( x k , x p ) + Σ q = 1 10000 Σ p = 1 10000 α p α q K ( x p , x q ) ;
Wherein, N1Quantity for described multiple dynamic reactive capacity configuration mode;Described D represents the element that pth row q arranges For Dpq=K (xp, xq)ypyqMatrix;Described c=(-1 ... ,-1)T, expression line number is N1, columns is the matrix of 1;Described Aα =(1,1 ..., 1), representing that line number is 1, columns is N1Matrix;DescribedykRepresent that kth kind is moved The risk factor that state reactive capability configuration mode is corresponding; Represent n-th destination node dynamic reactive merit after fault occurs under kth kind dynamic reactive capacity configuration mode The maximum of rate;K=1,2 ..., N1, n=1,2 ..., Nsvc, NsvcQuantity for described destination node.
In one embodiment, step S150 can be particularly as follows: ask relationship by objective (RBO) formula in the first optimal conditionsUnder optimal solution Wherein ykFor the risk factor that kth kind dynamic reactive capacity configuration mode is corresponding;QsvcTotal dynamic reactive for network system holds Amount; Represent under kth kind dynamic reactive capacity configuration mode the The maximum of n destination node dynamic reactive power after fault occurs, NsvcFor the quantity of destination node, k in network system =1,2 ..., N1, n=1,2 ..., Nsvc
In the present embodiment, minimum as the first optimal conditions using risk factor, try to achieve the complete of dynamic reactive capacity configuration Office's optimal solution, configures the dynamic reactive capacity of destination node each in network system according to this optimal solution, can be subject at electrical network At utmost maintain stablizing of line voltage, the fluctuation of suppression line voltage during to different disturbances or fault, thus be effectively improved The stability of electrical network and reliability.
Fig. 3 is the structural representation of the configuration system of the electrical network dynamic reactive capacity of one embodiment of the invention.Such as Fig. 3 institute Showing, the configuration system of this electrical network dynamic reactive capacity comprises the steps that
First computing module 310, for according to each destination node fault occur after voltage and default normal voltage Scope, calculates the risk factor that multiple dynamic reactive capacity configuration mode is the most corresponding;
Fitting module 330, dynamic for each destination node corresponding to dynamic reactive capacity configuration mode multiple described in matching State reactive power and the risk factor of correspondence, obtain relationship by objective (RBO) formula;
Second computing module 350, based on being optimized described relationship by objective (RBO) formula according to the first optimal conditions preset Calculate, obtain the optimal solution of the dynamic reactive capacity of each destination node;
Configuration module 370, for according to described in the optimal solution corresponding configuration of the dynamic reactive capacity of described each destination node The dynamic reactive capacity of each destination node.
In one embodiment, the first computing module 310 can be specifically for according to presetting constraints, accounting equationObtain described multiple dynamic reactive capacity configuration side The risk factor that formula is corresponding
Wherein, described default constraints includes system load flow constraints, each busbar voltage constraints, transmission line The constraints such as or not phase angle difference constraints and control variable;
Described risk (k) is the risk factor that kth kind dynamic reactive capacity configuration mode is corresponding, described NsvcFor described mesh The quantity of mark node;DescribedFor under kth kind dynamic reactive capacity configuration mode, the n-th destination node is sent out the s fault Voltage after life;DescribedIt it is the upper voltage limit in the normal voltage range of the n-th destination node;DescribedV n Represent the n-th mesh Lower voltage limit in the normal voltage range of mark node;DescribedWith In characterizing kth kind dynamic reactive capacity configuration mode, whereinRepresent under kth kind dynamic reactive capacity configuration mode n-th The maximum of individual destination node dynamic reactive power after fault occurs, n=1,2 ..., Nsvc, NsvcFor described destination node Quantity.
Specifically, said system trend constraints includes network system trend constraints under non-faulting stateAnd network system nonserviceable under trend constraintsWherein N is the quantity of network system interior joint.WithRepresent respectively under kth kind dynamic reactive capacity configuration mode, node i during electrical network properly functioning (being i.e. in non-faulting state) Meritorious and reactive power;WithRepresent respectively under kth kind dynamic reactive capacity configuration mode, after fault s occurs Meritorious and the reactive power of node i.WithBeing illustrated respectively under kth kind dynamic reactive capacity configuration mode, electrical network is just Node i and the voltage of node j when often running;WithIt is illustrated respectively under kth kind dynamic reactive capacity configuration mode, There is posterior nodal point i and the voltage of node j in fault s.GijAnd BijRepresent the conductance between node i and node j and susceptance respectively.WithRepresent respectively under kth kind dynamic reactive capacity configuration mode, after electrical network is properly functioning and fault s occurs, node i And the phase angle difference between node j.
Specifically, above-mentioned each busbar voltage constraints includesWherein,WithV i Represent joint respectively The upper voltage limit of some i and lower limit,Represent under kth kind dynamic reactive capacity configuration mode, node i when electrical network is properly functioning Voltage.
Specifically, above-mentioned transmission line phase angle difference constraints includesWherein,Withδ ij Respectively The upper voltage limit of the phase angle difference between expression node i and node j and lower limit,Represent in kth kind dynamic reactive capacity configuration Under mode, phase angle difference between node i and node j when electrical network is properly functioning.
Specifically, the constraints such as or not above-mentioned control variable includesWherein,WithQ i Represent the reactive power upper limit and the reactive power lower limit of node i respectively,WithP i Represent the active power of node i respectively The upper limit and active power lower limit,WithBeing illustrated respectively under kth kind dynamic reactive capacity configuration mode, electrical network is normally transported The active power of node i and reactive power during row.WithIt is illustrated respectively in kth kind dynamic reactive capacity configuration mode Under, the reactive power of destination node n when electrical network is properly functioning and after fault s occurs.Represent in the dynamic nothing of kth kind Under power capacity amount configuration mode, the dynamic reactive power value that destination node n sends after fault s occurs.With Being illustrated respectively under kth kind dynamic reactive capacity configuration mode, the dynamic reactive power that destination node n sends after a failure is maximum Value and minima.
In one embodiment, as shown in Figure 4, fitting module 330 can include computing unit 331 and obtain unit 333, its In:
Computing unit 331 is for according to gaussian kernel function and the second optimal conditions calculation optimization based on described risk factor Multiplier, the variable of wherein said gaussian kernel function is included in each destination node in the case of described multiple dynamic reactive capacity configuration Dynamic reactive power.
Obtain unit 333 for the coefficient using described optimization multiplier as described gaussian kernel function, it is thus achieved that described target is closed It it is formula.
Specifically, above-mentioned gaussian kernel function can beSecond optimal conditions isOptimization multiplier is α=(α1, α2..., αk);Relationship by objective (RBO) formula is
y k = f ( x k ) = K ( x k , x k ) - 2 Σ p = 1 10000 α p K ( x k , x p ) + Σ q = 1 10000 Σ p = 1 10000 α p α q K ( x p , x q ) ;
Wherein, N1Quantity for multiple dynamic reactive capacity configuration mode;D represents that the element that pth row q arranges is Dpq=K (xp, xq)ypyqMatrix;C=(-1 ... ,-1)T, expression line number is N1, columns is the matrix of 1;Aα=(1,1 ..., 1), table Showing that line number is 1, columns is N1Matrix;ykRepresent that kth kind dynamic reactive capacity configuration mode is corresponding Risk factor; Represent kth kind dynamic reactive capacity configuration The maximum of n-th destination node dynamic reactive power after fault occurs under mode;K=1,2 ..., N1, n=1, 2 ..., Nsvc, NsvcQuantity for destination node.
In one embodiment, above-mentioned first optimal conditions isWherein QsvcTotal dynamic reactive capacity for network system.
The configuration system of above-mentioned electrical network dynamic reactive capacity, according to each mesh that multiple dynamic reactive capacity configuration mode is corresponding The mark dynamic reactive power of node and risk factor when tackling various fault, fit object relational expression also seeks optimal solution, due to This optimal solution is that the overall situation considers the risk factor that multiple dynamic reactive capacity configuration mode tackles various fault and obtains, to electricity The various faults that net system is likely to occur have general applicability, therefore according to this optimal solution to destination node each in network system Dynamic reactive capacity configure, can at utmost maintain the steady of line voltage when electrical network is by different disturbances or fault Fixed, the fluctuation of suppression line voltage, thus it is effectively improved stability and the reliability of electrical network.
Should be noted that in said system embodiment, included modules simply carries out drawing according to function logic Point, but it is not limited to above-mentioned division, as long as being capable of corresponding function;It addition, each functional module is concrete Title also only to facilitate mutually distinguish, is not limited to protection scope of the present invention.
It addition, one of ordinary skill in the art will appreciate that all or part of step realizing in the various embodiments described above method The program that can be by completes to instruct relevant hardware, and corresponding program can be stored in read/write memory medium.
Each technical characteristic of embodiment described above can combine arbitrarily, for making description succinct, not to above-mentioned reality The all possible combination of each technical characteristic executed in example is all described, but, as long as the combination of these technical characteristics is not deposited In contradiction, all it is considered to be the scope that this specification is recorded.
Embodiment described above only have expressed the several embodiments of the present invention, and it describes more concrete and detailed, but also Can not therefore be construed as limiting the scope of the patent.It should be pointed out that, come for those of ordinary skill in the art Saying, without departing from the inventive concept of the premise, it is also possible to make some deformation and improvement, these broadly fall into the protection of the present invention Scope.Therefore, the protection domain of patent of the present invention should be as the criterion with claims.

Claims (10)

1. the collocation method of an electrical network dynamic reactive capacity, it is characterised in that including:
According to each destination node fault occur after voltage and default normal voltage range, calculate multiple dynamic reactive capacity The risk factor that configuration mode is the most corresponding;
The dynamic reactive power of each destination node that multiple dynamic reactive capacity configuration mode described in matching is corresponding and the wind of correspondence Danger coefficient, obtains relationship by objective (RBO) formula;
According to the first default optimal conditions, described relationship by objective (RBO) formula is optimized calculating, obtains the dynamic nothing of each destination node The optimal solution of power capacity amount;And
The dynamic reactive of each destination node described in the optimal solution corresponding configuration of the dynamic reactive capacity according to described each destination node Capacity.
The collocation method of electrical network dynamic reactive capacity the most according to claim 1, it is characterised in that many described in described matching Plant dynamic reactive power and the risk factor of correspondence of each destination node corresponding to dynamic reactive capacity configuration mode, obtain target Relational expression, including:
According to gaussian kernel function and the second optimal conditions calculation optimization multiplier based on described risk factor, wherein said gaussian kernel The variable of function is included in the dynamic reactive power of each destination node in the case of described multiple dynamic reactive capacity configuration;And
Using described optimization multiplier as the coefficient of described gaussian kernel function, it is thus achieved that described relationship by objective (RBO) formula.
The collocation method of electrical network dynamic reactive capacity the most according to claim 2, its feature exists In, described gaussian kernel function isDescribed second optimal conditions isDescribed optimization multiplier is α=(α1, α2..., αk);Described relationship by objective (RBO) formula is
Wherein, N1Quantity for described multiple dynamic reactive capacity configuration mode;Described D represents that the element that pth row q arranges is Dpq =K (xp, xq)ypyqMatrix;Described c=(-1 ... ,-1)T, expression line number is N1, columns is the matrix of 1;Described Aα=(1, 1 ..., 1), representing that line number is 1, columns is N1Matrix;DescribedykRepresent kth kind dynamic reactive The risk factor that capacity configuration mode is corresponding; Represent the maximum of n-th destination node dynamic reactive power after fault occurs under kth kind dynamic reactive capacity configuration mode Value;K=1,2 ..., N1, n=1,2 ..., Nsvc, NsvcQuantity for described destination node.
The collocation method of electrical network dynamic reactive capacity the most according to claim 3, it is characterised in that described first optimizes bar Part isWherein QsvcTotal dynamic reactive capacity for network system.
The collocation method of electrical network dynamic reactive capacity the most according to claim 1, it is characterised in that described according to each target Node fault occur after voltage and default normal voltage range, calculate multiple dynamic reactive capacity configuration mode the most right The risk factor answered, including:
According to default constraints, accounting equationObtain described The risk factor that multiple dynamic reactive capacity configuration mode is corresponding
Wherein, described default constraints includes system load flow constraints, each busbar voltage constraints, transmission line phase angle Difference constraints and the constraints such as or not control variable;
Described risk (k) is the risk factor that kth kind dynamic reactive capacity configuration mode is corresponding, described NsvcSave for described target The quantity of point;DescribedFor under kth kind dynamic reactive capacity configuration mode, the n-th destination node is after the s fault occurs Voltage;Described VnIt it is the upper voltage limit in the normal voltage range of the n-th destination node;DescribedV n Represent the n-th destination node Normal voltage range in lower voltage limit;DescribedFor characterizing Kth kind dynamic reactive capacity configuration mode, whereinRepresent the n-th target under kth kind dynamic reactive capacity configuration mode The maximum of node dynamic reactive power after fault occurs, n=1,2 ..., Nsvc, NsvcNumber for described destination node Amount.
6. the configuration system of an electrical network dynamic reactive capacity, it is characterised in that including:
First computing module, for according to each destination node fault occur after voltage and default normal voltage range, meter Calculate the risk factor that multiple dynamic reactive capacity configuration mode is the most corresponding;
Fitting module, for the dynamic reactive merit of each destination node corresponding to dynamic reactive capacity configuration mode multiple described in matching Rate and the risk factor of correspondence, obtain relationship by objective (RBO) formula;
Second computing module, for described relationship by objective (RBO) formula being optimized calculating according to the first optimal conditions preset, obtains The optimal solution of the dynamic reactive capacity of each destination node;And
Configuration module, for saving according to each target described in the optimal solution corresponding configuration of the dynamic reactive capacity of described each destination node The dynamic reactive capacity of point.
The configuration system of electrical network dynamic reactive capacity the most according to claim 6, it is characterised in that described fitting module bag Include:
Computing unit, is used for according to gaussian kernel function and the second optimal conditions calculation optimization multiplier based on described risk factor, The variable of wherein said gaussian kernel function is included in the dynamic of each destination node in the case of described multiple dynamic reactive capacity configuration Reactive power;And
Obtain unit, for the coefficient using described optimization multiplier as described gaussian kernel function, it is thus achieved that described relationship by objective (RBO) formula.
The configuration system of electrical network dynamic reactive capacity the most according to claim 7, its feature exists In, described gaussian kernel function isDescribed second optimal conditions isDescribed optimization multiplier is α=(α1, α2..., αk);Described relationship by objective (RBO) formula is
Wherein, N1Quantity for described multiple dynamic reactive capacity configuration mode;Described D represents that the element that pth row q arranges is Dpq =K (xp, xq)ypyqMatrix;Described c=(-1 ... ,-1)T, expression line number is N1, columns is the matrix of 1;Described Aα=(1, 1 ..., 1), representing that line number is 1, columns is N1Matrix;DescribedykRepresent kth kind dynamic reactive The risk factor that capacity configuration mode is corresponding; Represent the maximum of n-th destination node dynamic reactive power after fault occurs under kth kind dynamic reactive capacity configuration mode Value;K=1,2 ..., N1, n=1,2 ..., Nsvc, NsvcQuantity for described destination node.
The configuration system of electrical network dynamic reactive capacity the most according to claim 8, it is characterised in that described first optimizes bar Part isWherein QsvcTotal dynamic reactive capacity for network system.
The configuration system of electrical network dynamic reactive capacity the most according to claim 6, it is characterised in that described first calculates Module is used for:
According to default constraints, accounting equationObtain institute State the risk factor that multiple dynamic reactive capacity configuration mode is corresponding
Wherein, described default constraints includes system load flow constraints, each busbar voltage constraints, transmission line phase angle Difference constraints and the constraints such as or not control variable;
Described risk (k) is the risk factor that kth kind dynamic reactive capacity configuration mode is corresponding, described NsvcSave for described target The quantity of point;DescribedFor under kth kind dynamic reactive capacity configuration mode, the n-th destination node is after the s fault occurs Voltage;DescribedIt it is the upper voltage limit in the normal voltage range of the n-th destination node;DescribedV n Represent the n-th target joint Lower voltage limit in the normal voltage range of point;DescribedFor table Levy kth kind dynamic reactive capacity configuration mode, whereinRepresent the n-th mesh under kth kind dynamic reactive capacity configuration mode The maximum of mark node dynamic reactive power after fault occurs, n=1,2 ..., Nsvc, NsvcNumber for described destination node Amount.
CN201610615753.2A 2016-07-28 2016-07-28 The configuration method and system of power grid dynamic reactive capacity Active CN106099946B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610615753.2A CN106099946B (en) 2016-07-28 2016-07-28 The configuration method and system of power grid dynamic reactive capacity

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610615753.2A CN106099946B (en) 2016-07-28 2016-07-28 The configuration method and system of power grid dynamic reactive capacity

Publications (2)

Publication Number Publication Date
CN106099946A true CN106099946A (en) 2016-11-09
CN106099946B CN106099946B (en) 2019-04-09

Family

ID=57478779

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610615753.2A Active CN106099946B (en) 2016-07-28 2016-07-28 The configuration method and system of power grid dynamic reactive capacity

Country Status (1)

Country Link
CN (1) CN106099946B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109038572A (en) * 2018-08-30 2018-12-18 国家电网公司华东分部 Consider the voltage security range determining method of AC-DC hybrid power grid failure
CN106099946B (en) * 2016-07-28 2019-04-09 广东电网有限责任公司电力调度控制中心 The configuration method and system of power grid dynamic reactive capacity
CN111260248A (en) * 2020-02-13 2020-06-09 东方电子股份有限公司 Distribution network fault self-healing scheduling method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2093855A2 (en) * 2008-02-25 2009-08-26 General Electric Company Systems and methods of dynamic reactive voltage support for a power transmission system
CN103124075A (en) * 2013-03-20 2013-05-29 东南大学 Reactive power configuration method for wind power base
CN103795068A (en) * 2014-03-05 2014-05-14 广东电网公司电力调度控制中心 Optimal configuration method for high-voltage distribution network dynamic reactive power compensation equipment capacity

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106099946B (en) * 2016-07-28 2019-04-09 广东电网有限责任公司电力调度控制中心 The configuration method and system of power grid dynamic reactive capacity

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2093855A2 (en) * 2008-02-25 2009-08-26 General Electric Company Systems and methods of dynamic reactive voltage support for a power transmission system
CN103124075A (en) * 2013-03-20 2013-05-29 东南大学 Reactive power configuration method for wind power base
CN103795068A (en) * 2014-03-05 2014-05-14 广东电网公司电力调度控制中心 Optimal configuration method for high-voltage distribution network dynamic reactive power compensation equipment capacity

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106099946B (en) * 2016-07-28 2019-04-09 广东电网有限责任公司电力调度控制中心 The configuration method and system of power grid dynamic reactive capacity
CN109038572A (en) * 2018-08-30 2018-12-18 国家电网公司华东分部 Consider the voltage security range determining method of AC-DC hybrid power grid failure
CN109038572B (en) * 2018-08-30 2021-05-07 国家电网公司华东分部 Voltage safety range determination method considering faults of alternating current-direct current hybrid power grid
CN111260248A (en) * 2020-02-13 2020-06-09 东方电子股份有限公司 Distribution network fault self-healing scheduling method

Also Published As

Publication number Publication date
CN106099946B (en) 2019-04-09

Similar Documents

Publication Publication Date Title
EP2891221B1 (en) Power distribution system loss reduction with distributed energy resource control
CN107317334B (en) A kind of electric system rack reconstructing method and device
CN106099946A (en) The collocation method of electrical network dynamic reactive capacity and system
CN104753061B (en) Distributed type power supply accessed into power distribution network and microgrid group zone control method and microgrid group zone control system
CN106384162B (en) Consider electric system electric energy-spare joint optimal operation method of zonal reserve and N-1 verification
Adhikari et al. Recovery risk mitigation of wind integrated bulk power system with flywheel energy storage
CN105048483A (en) Energy-storage capacity configuration method and system for PV-containing power distribution network
CN107346267A (en) A kind of cpu performance optimization method and device based on NUMA architecture
CN104123683A (en) Electrical power system black-start scheme generation method based on dynamic programming
CN103310296A (en) Operation ticket security check method based on disturbance evaluation and trend analysis
CN107565550A (en) A kind of power distribution network partition method and system
CN109962485A (en) A kind of composite energy storing device addressing constant volume method towards source net lotus close friend interaction
CN107516909B (en) It can access the optimization method and device of wind power output in a kind of rack restructuring procedure
CN109560552A (en) Pricing method, device, computer equipment and the storage medium of power grid distribution
Shahgholian et al. Impact of PSS and STATCOM on dynamic parameters of power system based on neuro-fuzzy controllers
CN116384212B (en) Collaborative loss reduction method for DG-containing area, computer equipment and storage medium
CN116093995B (en) Multi-target network reconstruction method and system for power distribution system
CN105279297A (en) PSS/E power system equivalent method based on engineering application
CN107179688A (en) Consider the Power System Reliability Analysis method of Monte Carlo state sampling truncation
CN116316607A (en) Combined optimization scheduling method, system, equipment and medium
CN108306302B (en) Voltage control method and device and storage medium
CN109286204A (en) A kind of distribution network black starting-up reconstructing method based on minimum expectation loss of outage
CN104218607A (en) New energy access system application based on artificial fish swarm and tabu search algorithms
Liu et al. A novel acceleration strategy for nl contingency screening in distribution system
CN107565578B (en) The reactive output control method and device of power distribution network

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210726

Address after: 510699 No. 75 Meihua Road, Yuexiu District, Guangzhou City, Guangdong Province

Patentee after: ELECTRIC POWER DISPATCHING CONTROL CENTER OF GUANGDONG POWER GRID Co.,Ltd.

Address before: 510699 No. 75 Meihua Road, Yuexiu District, Guangzhou City, Guangdong Province

Patentee before: ELECTRIC POWER DISPATCHING CONTROL CENTER OF GUANGDONG POWER GRID Co.,Ltd.

Patentee before: BEIJING KING STAR HI-TECH SYSTEM CONTROL Co.,Ltd.