CN109755942A - Extension trend method and device based on optimization method - Google Patents
Extension trend method and device based on optimization method Download PDFInfo
- Publication number
- CN109755942A CN109755942A CN201711064949.8A CN201711064949A CN109755942A CN 109755942 A CN109755942 A CN 109755942A CN 201711064949 A CN201711064949 A CN 201711064949A CN 109755942 A CN109755942 A CN 109755942A
- Authority
- CN
- China
- Prior art keywords
- bus
- busbar voltage
- optimized model
- equation
- optimization
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 93
- 238000005457 optimization Methods 0.000 title claims abstract description 76
- 239000011159 matrix material Substances 0.000 claims abstract description 45
- 238000002347 injection Methods 0.000 claims abstract description 42
- 239000007924 injection Substances 0.000 claims abstract description 42
- 238000011990 functional testing Methods 0.000 claims abstract description 11
- 230000021615 conjugation Effects 0.000 claims description 16
- 238000013507 mapping Methods 0.000 claims description 16
- 239000013598 vector Substances 0.000 claims description 16
- 230000005611 electricity Effects 0.000 claims description 9
- 230000009021 linear effect Effects 0.000 claims description 9
- 235000008331 Pinus X rigitaeda Nutrition 0.000 claims description 4
- 235000011613 Pinus brutia Nutrition 0.000 claims description 4
- 241000018646 Pinus brutia Species 0.000 claims description 4
- 238000009415 formwork Methods 0.000 claims 3
- 230000035945 sensitivity Effects 0.000 description 7
- 230000002040 relaxant effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000009022 nonlinear effect Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The extension trend method and device based on optimization method that the invention discloses a kind of, wherein method includes: to limit equation according to the busbar voltage that limitation formula establishes relaxation;The tide optimization problem for having partial voltage control device is established by optimization formula;The Optimized model, the Optimized model of PV bus and the Optimized model of PQ bus of zero injection bus are established by constraint equation;Optimized model, the Optimized model of PV bus and the Optimized model of PQ bus of zero injection bus are solved by the IPM of constant Hessian matrix, to obtain integer tap position, and carry out functional test under the conditions of three-phase equilibrium and imbalance.The convergence of nonlinear-sensitivity method under fully loaded transportation condition can be improved in this method, suitable for infeasible voltage limitation is detected and relaxed under fully loaded transportation condition, it is ensured that Power flow solvability.
Description
Technical field
The present invention relates to Power System Expansion band wagon technique field, in particular to a kind of extension trend based on optimization method
Method and device.
Background technique
In the related technology partial voltage control device, this method application constant impedance square are coped with based on the method for sensitivity
Battle array exports the nonlinear-sensitivity matrix of partial voltage control device.Since the non-linear property of trend is very tight under fully loaded transportation condition
Weight, the distributed network trend based on sensitivity method reply partial voltage control device need more iteration under case of heavy load
Number.
The nonlinear-sensitivity matrix of voltage control distributed generator can also be derived based on loop analysis, it was demonstrated that
Nonlinear-sensitivity method has better convergence than linear sensitivity method, however, under case of heavy load, nonlinear-sensitivity
Method might have bad convergence, and especially after case of heavy load has adjusted the position of tap switch, busbar voltage can
It can be limited still above voltage, it would be highly desirable to solve.
Summary of the invention
The present invention is directed to solve at least some of the technical problems in related technologies.
For this purpose, an object of the present invention is to provide a kind of extension trend method based on optimization method, this method can
To improve the convergence of nonlinear-sensitivity method under fully loaded transportation condition, suitable for detecting and relaxing infeasible under fully loaded transportation condition
Voltage limitation.
It is another object of the present invention to propose a kind of extension trend device based on optimization method.
In order to achieve the above objectives, one aspect of the present invention embodiment proposes a kind of extension trend side based on optimization method
Method, comprising the following steps: equation is limited according to the busbar voltage that limitation formula establishes relaxation;It is established by optimization formula and has office
The tide optimization problem of portion's voltage controller;The optimization of the Optimized model, PV bus of zero injection bus is established by constraint equation
The Optimized model of model and PQ bus;Pass through IPM (Intelligent Power Module, the intelligent function of constant Hessian matrix
Rate module) to Optimized model, the Optimized model of the PV bus and the Optimized model of the PQ bus of the zero injection bus
It is solved, to obtain integer tap position, and carries out functional test under the conditions of three-phase equilibrium and imbalance.
The extension trend method based on optimization method of the embodiment of the present invention, by the IPM of constant Hessian matrix to zero note
The Optimized model for entering the Optimized model of bus, the Optimized model of PV bus and PQ bus is solved, to obtain integer tap
Position, and functional test is carried out under the conditions of three-phase equilibrium and imbalance, realize the mesh of the extension trend based on optimization method
, it solves the non-linear of distributed network trend based on sensitive method reply partial voltage control device under fully loaded transportation condition and asks
The problem of detecting under topic and fully loaded transportation condition and relaxing the limitation of infeasible voltage, improves nonlinear-sensitivity method under fully loaded transportation condition
Convergence, suitable for the limitation of infeasible voltage is detected and relaxed under fully loaded transportation condition.
In addition, the extension trend method according to the above embodiment of the present invention based on optimization method can also have it is following attached
The technical characteristic added:
Further, in one embodiment of the invention, the busbar voltage limitation equation is established by limiting formula,
The limitation formula are as follows:
Wherein,For the lower limit of busbar voltage square,For the upper limit of the busbar voltage square, l and u are indicated
Non-negative slack variable, UreAnd UimThe respectively real and imaginary parts of busbar voltage complex variable;guFor the function upper limit, g is mapping, gl
For function lower limit,For busbar voltage square.
Further, in one embodiment of the invention, the tide optimization problem is established by following optimization formula,
The optimization formula are as follows:
min∑(u+v)
Wherein, t is transformer tap positions, tminAnd tmaxThe respectively described transformer tap positions minimum value limitation and
Maximum value limitation, h are mapping, and v is non-negative slack variable, and s.t. is constraint condition.
Further, in one embodiment of the invention, the Optimized model for establishing zero injection bus, PV bus
The Optimized model of Optimized model and PQ bus further comprises:
For the zero injection bus, the first constraint equation h is indicated with pairs of linear equation:
Wherein, IreAnd IimThe respectively real and imaginary parts of Injection Current, YzeroFor corresponding to zero injection bus admittance row,
U is multiple busbar voltage,WithThe respectively real and imaginary parts of equation.
For the PV bus, the second constraint equation h is indicated in pairs are as follows:
Wherein,Target is controlled for the voltage of PV bus, subscript PV is the respective vectors or row matrix of PV bus, diag
For diagonal matrix, * is the conjugation of complex variable, SreFor inject complex power real part,For the real part of equation, YPVFor corresponding to PQ
The admittance row of bus, UPVFor PV busbar voltage, U is multiple busbar voltage,For the real part of PV busbar voltage square,For
The imaginary part of PV busbar voltage square.
For the PQ bus, third constraint equation h is indicated in pairs are as follows:
Wherein, SreAnd SimThe real and imaginary parts of complex power are respectively injected,WithThe respectively real and imaginary parts of equation,
Diag is diagonal matrix, and * is the conjugation of complex variable, and subscript PQ is the respective vectors or row matrix of PQ bus, UPQFor PQ bus electricity
Pressure, YPQFor the admittance row corresponding to PQ bus, U is multiple busbar voltage.
Further, in one embodiment of the invention, described in the fractional value acquisition by the tap position that rounds up
Integer tap position.
In order to achieve the above objectives, another aspect of the present invention embodiment proposes a kind of extension trend dress based on optimization method
It sets, comprising: limitation establishing equation module, the busbar voltage for establishing relaxation according to limitation formula limit equation;Optimization problem
Module is established, for establishing the tide optimization problem for having partial voltage control device by optimization formula;Optimized model establishes mould
Block, for establishing the Optimized model, the Optimized model of PV bus and the optimization mould of PQ bus of zero injection bus by constraint equation
Type;Trend expansion module, for by the IPM of constant Hessian matrix to the Optimized model of the zero injection bus, PV bus
The Optimized model of Optimized model and PQ bus is solved, to obtain integer tap position, and in three-phase equilibrium and imbalance
Under the conditions of carry out functional test.
The extension trend device based on optimization method of the embodiment of the present invention, by the IPM of constant Hessian matrix to zero note
The Optimized model for entering the Optimized model of bus, the Optimized model of PV bus and PQ bus is solved, to obtain integer tap
Position, and functional test is carried out under the conditions of three-phase equilibrium and imbalance, realize the mesh of the extension trend based on optimization method
, it solves the non-linear of distributed network trend based on sensitive method reply partial voltage control device under fully loaded transportation condition and asks
The problem of detecting under topic and fully loaded transportation condition and relaxing the limitation of infeasible voltage, improves nonlinear-sensitivity method under fully loaded transportation condition
Convergence, suitable for the limitation of infeasible voltage is detected and relaxed under fully loaded transportation condition.
In addition, the extension trend device according to the above embodiment of the present invention based on optimization method can also have it is following attached
The technical characteristic added:
Further, in one embodiment of the invention, formula is limited described in the limitation establishing equation module are as follows:
Wherein,For the lower limit of busbar voltage square,For the upper limit of the busbar voltage square, l and u difference
For non-negative slack variable, UreAnd UimThe respectively real and imaginary parts of busbar voltage complex variable, guFor the function upper limit, g is mapping,
glFor function lower limit,For busbar voltage square.
Further, in one embodiment of the invention, the optimization problem, which is established, optimizes formula described in module are as follows:
min∑(u+v)
Wherein, t is transformer tap positions, tminAnd tmaxThe respectively described transformer tap positions minimum value limitation and
Maximum value limitation, h are mapping, and v is non-negative slack variable, and s.t. is constraint condition.
Further, in one embodiment of the invention, the Optimized model establishes module and further comprises:
For the zero injection bus, the first constraint equation h is indicated with pairs of linear equation:
Wherein, IreAnd IimThe respectively real and imaginary parts of Injection Current, YzeroFor corresponding to zero injection bus admittance row,
U is multiple busbar voltage,WithThe respectively real and imaginary parts of equation.
For the PV bus, the second constraint equation h is indicated in pairs are as follows:
Wherein,Target is controlled for the voltage of PV bus, subscript PV is the respective vectors or row matrix of PV bus,
Diag is diagonal matrix, and * is the conjugation of complex variable, SreFor inject complex power real part,For the real part of equation, YPVFor correspondence
In the admittance row of PV bus, UPVFor PV busbar voltage, U is multiple busbar voltage,For the real part of PV busbar voltage square,For the imaginary part of PV busbar voltage square.
For the PQ bus, third constraint equation h is indicated in pairs are as follows:
Wherein, SreAnd SimThe real and imaginary parts of complex power are respectively injected,WithThe respectively real and imaginary parts of equation,
Diag is diagonal matrix, and * is the conjugation of complex variable, and subscript PQ is the respective vectors or row matrix of PQ bus, UPQFor PQ bus electricity
Pressure, YPQFor the admittance row corresponding to PQ busbar voltage, U is multiple busbar voltage.
Further, in one embodiment of the invention, described in the fractional value acquisition by the tap position that rounds up
Integer tap position.
The additional aspect of the present invention and advantage will be set forth in part in the description, and will partially become from the following description
Obviously, or practice through the invention is recognized.
Detailed description of the invention
Above-mentioned and/or additional aspect and advantage of the invention will become from the following description of the accompanying drawings of embodiments
Obviously and it is readily appreciated that, in which:
Fig. 1 is the flow chart according to the extension trend method based on optimization method of the embodiment of the present invention;
Fig. 2 is the structural schematic diagram according to the extension trend device based on optimization method of the embodiment of the present invention.
Specific embodiment
The embodiment of the present invention is described below in detail, examples of the embodiments are shown in the accompanying drawings, wherein from beginning to end
Same or similar label indicates same or similar element or element with the same or similar functions.Below with reference to attached
The embodiment of figure description is exemplary, it is intended to is used to explain the present invention, and is not considered as limiting the invention.
The extension trend method and dress based on optimization method proposed according to embodiments of the present invention is described with reference to the accompanying drawings
It sets, describes the extension trend method based on optimization method proposed according to embodiments of the present invention with reference to the accompanying drawings first.
Fig. 1 is that the extension trend method based on optimization method of the embodiment of the present invention is flow chart.
As shown in Figure 1, should extension trend method based on optimization method the following steps are included:
In step s101, equation is limited according to the busbar voltage that limitation formula establishes relaxation.
Further, in one embodiment of the invention, busbar voltage limitation equation, limitation are established by limiting formula
Formula are as follows:
Wherein,For the lower limit of busbar voltage square,For the upper limit of busbar voltage square, l and u are non-negative pine
Relaxation variable, UreAnd UimThe respectively real and imaginary parts of busbar voltage complex variable, guFor the function upper limit, g is mapping, glFor function
Lower limit,For square of busbar voltage.
In step s 102, the tide optimization problem for having partial voltage control device is established by optimization formula.
Further, in one embodiment of the invention, tide optimization problem is established by following optimization formula, optimized
Formula are as follows:
min∑(u+v)
Wherein, t is transformer tap positions, tminAnd tmaxThe respectively minimum value limitation and maximum of transformer tap positions
Value limitation, h are mapping, and v is non-negative slack variable, and s.t. is constraint condition.
In step s 103, by constraint equation establish the Optimized model of zero injection bus, PV bus Optimized model and
The Optimized model of PQ bus.
Further, in one embodiment of the invention, the optimization of the Optimized model, PV bus of zero injection bus is established
The Optimized model of model and PQ bus further comprises:
For zero injection bus, the first constraint equation h is indicated with pairs of linear equation:
Wherein, IreAnd IimThe respectively real and imaginary parts of Injection Current, YzeroFor corresponding to zero injection bus admittance row,
U is multiple busbar voltage,WithThe respectively real and imaginary parts of equation.
For PV bus, the second constraint equation h is indicated in pairs are as follows:
Wherein,Target is controlled for the voltage of PV bus;Subscript PV is the respective vectors or row matrix of PV bus, diag
For diagonal matrix, * is the conjugation of complex variable, SreFor inject complex power real part,For the real part Y of equationPVFor corresponding to PV mother
The admittance row of line, UPVFor PV busbar voltage, U is multiple busbar voltage,For the real part of PV busbar voltage square,For
The imaginary part of PV busbar voltage square.
For PQ bus, third constraint equation h is indicated in pairs are as follows:
Wherein, SreAnd SimThe real and imaginary parts of complex power are respectively injected,WithThe respectively real and imaginary parts of equation,
Diag is diagonal matrix, and * is the conjugation of complex variable, and subscript PQ is the respective vectors or row matrix of PQ bus, UPQFor PQ bus electricity
Pressure, YPQFor the admittance row corresponding to PQ busbar voltage, U is multiple busbar voltage.
In step S104, pass through the optimization of the IPM of the constant Hessian matrix Optimized model, PV bus for injecting bus to zero
The Optimized model of model and PQ bus is solved, to obtain integer tap position, and in three-phase equilibrium and uneven condition
Lower carry out functional test.
Further, in one embodiment of the invention, integer is obtained by the fractional value for the tap position that rounds up
Tap position.
In one particular embodiment of the present invention, under condition of three-phase balance, the embodiment of the present invention propose based on excellent
The extension trend method of change method the following steps are included:
S1: the distributed network based on sensitivity method reply partial voltage control device not in view of non-linear property is collected
Data of the trend under underloading or heavy duty.
S2: loose busbar voltage limitation equation is established.
Wherein,For the lower limit of busbar voltage square,For the upper limit of busbar voltage square, l and u are non-negative pine
Relaxation variable, UreAnd UimThe respectively real and imaginary parts of busbar voltage complex variable, guFor the function upper limit, g is mapping, glFor function
Lower limit,For square of busbar voltage.
S3: the optimization problem that the trend with partial voltage control device extends is established.
min∑(u+v)
Wherein, UreAnd UimThe respectively real and imaginary parts of busbar voltage complex variable, t are transformer tap positions, tmin
And tmaxRespectively the minimum value limitation and maximum value limitation of transformer tap positions, h are mapping, and v is non-negative slack variable,
It s.t. is constraint condition.
S4: zero injection bus constraint equation is established.
Wherein, IreAnd IimThe respectively real and imaginary parts of Injection Current, YzeroFor corresponding to zero injection bus admittance row,
U is multiple busbar voltage,WithThe respectively real and imaginary parts of equation.
S5: PV bus constraint equation is established.
Wherein,Target is controlled for the voltage of PV bus, subscript PV is the respective vectors or row matrix of PV bus, diag
For diagonal matrix, * is the conjugation of complex variable, SreFor the real part for injecting complex power;For the real part of equation, YPVFor corresponding to PV
The admittance row of bus, UPVFor PV busbar voltage, U is multiple busbar voltage,For the real part of PV busbar voltage square,For
The imaginary part of PV busbar voltage square.
S6: PQ bus constraint equation is established.
Wherein, SreAnd SimThe real and imaginary parts of complex power are respectively injected,WithThe respectively real and imaginary parts of equation,
Diag is diagonal matrix, and * is the conjugation of complex variable, and subscript PQ is the respective vectors or row matrix of PQ bus, UPQFor PQ bus electricity
Pressure, YPQFor the admittance row corresponding to PQ busbar voltage, U is multiple busbar voltage.
S7: Matlab is applied, and solves above-mentioned optimization problem with the IPM with constant Hessian matrix.
S8: the number of iterations and the time needed for contrast sensitivity method and extension trend method, the embodiment of the present invention is verified
Method functionality.
In another specific embodiment of the invention, under non-equilibrium condition, the embodiment of the present invention propose based on excellent
The extension trend method of change method the following steps are included:
S1: loose busbar voltage limitation equation is established.
Wherein,For the lower limit of busbar voltage square,For the upper limit of busbar voltage square, l and u are non-negative pine
Relaxation variable UreAnd UimThe respectively real and imaginary parts of busbar voltage complex variable, guFor the function upper limit, g is mapping, glFor function
Lower limit,For square of busbar voltage.
S2: the optimization problem that the trend with partial voltage control device extends is established.
min∑(u+v)
Wherein, UreAnd UimThe respectively real and imaginary parts of busbar voltage complex variable, t are transformer tap positions, tmin
And tmaxRespectively the minimum value limitation and maximum value limitation of transformer tap positions, h are mapping, and v is non-negative slack variable,
It s.t. is constraint condition.
S3: zero injection bus constraint equation is established.
Wherein, IreAnd IimThe respectively real and imaginary parts of Injection Current, YzeroFor corresponding to zero injection bus admittance row,
U indicates multiple busbar voltage,WithThe respectively real and imaginary parts of equation.
S4: PV bus constraint equation is established.
Wherein,Target is controlled for the voltage of PV bus, subscript PV is the respective vectors or row matrix of PV bus, diag
For diagonal matrix, * is the conjugation of complex variable, SreFor the real part for injecting complex power;For the real part of equation, YPVFor corresponding to PV
The admittance row of bus, UPVFor PV busbar voltage, U is multiple busbar voltage,For the real part of PV busbar voltage square,For
The imaginary part of PV busbar voltage square.
S5: PQ bus constraint equation is established.
Wherein, SreAnd SimThe real and imaginary parts of complex power are respectively injected,WithThe respectively real and imaginary parts of equation,
Diag is diagonal matrix, and * is the conjugation of complex variable, and subscript PQ is the respective vectors or row matrix of PQ bus, UPQFor PQ bus electricity
Pressure, YPQFor the admittance row corresponding to PQ busbar voltage, U is multiple busbar voltage.
S6: applying Matlab, solves above-mentioned optimization problem with the IPM with constant Hessian matrix.
S7: according to [I.Kocar, J.Mahseredjian, U.Karaagac, G.Soykan, and O.Saad,
“Multiphase Load-Flow Solution for Large-Scale Distribution Systems Using
MANA, " IEEE Transactions on Power Delivery, vol.29, no.2, pp.908-915, Apr.2014.] it can
Know, step voltage adjuster model is indicated by Voltage Control Voltage Source and Current Controlled Current Source.
S8: [W.H.Kersting, Distribution System Modeling and Analysis, Third are collected
Edition.Boca Raton, Florida:CRC Press, 2012.] the IEEE4 bus in tests feeder line data.
S9: the number of iterations of the power flow solutions of step voltage regulator adjustment and present invention method is compared
Right, the method for verifying the embodiment of the present invention has the advantages that the number of iterations is few.
S10: by the C phase load on IEEE4 bus by multiplying 2, to obtain the infeasible electricity of present invention method
Pressure limitation, verify the embodiment of the present invention method have detection can not voltage limitation and relaxation capabilities functionality.
The extension trend method based on optimization method proposed according to embodiments of the present invention, passes through constant Hessian matrix
IPM solves Optimized model, the Optimized model of PV bus and the Optimized model of PQ bus of zero injection bus, to obtain
Integer tap position, and functional test is carried out under the conditions of three-phase equilibrium and imbalance, realize the extension based on optimization method
The purpose of trend solves the distributed network trend based on sensitive method reply partial voltage control device under fully loaded transportation condition
The problem of detecting under nonlinear problem and fully loaded transportation condition and relaxing the limitation of infeasible voltage, improves non-thread personality under fully loaded transportation condition
The convergence of sensitivity method, suitable for infeasible voltage limitation is detected and relaxed under fully loaded transportation condition.
The extension trend device based on optimization method proposed according to embodiments of the present invention is described referring next to attached drawing.
Fig. 2 is the structural schematic diagram of the extension trend device based on optimization method of the embodiment of the present invention.
As shown in Fig. 2, being somebody's turn to do the extension trend device 10 based on optimization method includes: limitation establishing equation module 100, optimization
Problem establishes module 200, Optimized model establishes module 300 and trend expansion module 400.
Wherein, the busbar voltage that limitation establishing equation module 100 is used to establish relaxation according to limitation formula limits equation.It is excellent
Change problem establish module 200 for by optimization formula establish with partial voltage control device tide optimization problem.Optimize mould
Type establishes module 300 and is used to establish the Optimized model of zero injection bus, the Optimized model of PV bus and PQ mother by constraint equation
The Optimized model of line.Trend expansion module 400 be used for by the IPM of constant Hessian matrix to the Optimized model of zero injection bus,
The Optimized model of PV bus and the Optimized model of PQ bus are solved, to obtain integer tap position, and three-phase equilibrium with
And functional test is carried out under the conditions of imbalance.Nonlinear-sensitivity under fully loaded transportation condition can be improved in the device 10 of the embodiment of the present invention
The convergence of method, suitable for infeasible voltage limitation is detected and relaxed under fully loaded transportation condition, it is ensured that Power flow solvability.
Further, in one embodiment of the invention, it limits in establishing equation module 100 and limits formula are as follows:
Wherein,For the lower limit of busbar voltage square,For the upper limit of busbar voltage square, l and u indicate non-negative
Slack variable, UreAnd UimThe respectively real and imaginary parts of busbar voltage complex variable, guFor the function upper limit, g is mapping, glFor letter
Number lower limit,For square of busbar voltage.
Further, in one embodiment of the invention, optimization problem, which is established in module 200, optimizes formula are as follows:
min∑(u+v)
Wherein, t is transformer tap positions, tminAnd tmaxThe respectively minimum value limitation and maximum of transformer tap positions
Value limitation, h are mapping, and v is non-negative slack variable, and s.t. is constraint condition.
Further, in one embodiment of the invention, Optimized model establishes module 300 and further comprises:
For zero injection bus, the first constraint equation h is indicated with pairs of linear equation:
Wherein, IreAnd IimThe respectively real and imaginary parts of Injection Current, YzeroFor corresponding to zero injection bus admittance row,
U is multiple busbar voltage,WithThe respectively real and imaginary parts of equation.
For PV bus, the second constraint equation h is indicated in pairs are as follows:
Wherein,Target is controlled for the voltage of PV bus, subscript PV is the respective vectors or row matrix of PV bus, diag
For diagonal matrix, * is the conjugation of complex variable, SreFor inject complex power real part,For the real part of equation, YPVFor corresponding to PV
The admittance row of bus, UPVFor PV busbar voltage, U is multiple busbar voltage,For the real part of PV busbar voltage square,For
The imaginary part of PV busbar voltage square.
For PQ bus, third constraint equation h is indicated in pairs are as follows:
Wherein, SreAnd SimThe real and imaginary parts of complex power are respectively injected,WithThe respectively real and imaginary parts of equation,
Diag is diagonal matrix, and * is the conjugation of complex variable, and subscript PQ is the respective vectors or row matrix of PQ bus, UPQFor PQ bus electricity
Pressure, YPQFor the admittance row corresponding to PQ busbar voltage, U is multiple busbar voltage.
Further not in one embodiment of the invention, integer is obtained by the fractional value for the tap position that rounds up
Tap position.
It should be noted that the aforementioned explanation to the extension trend embodiment of the method based on optimization method is also applied for
The extension trend device based on optimization method of the embodiment, details are not described herein again.
The extension trend device based on optimization method proposed according to embodiments of the present invention, passes through constant Hessian matrix
IPM solves Optimized model, the Optimized model of PV bus and the Optimized model of PQ bus of zero injection bus, to obtain
Integer tap position, and functional test is carried out under the conditions of three-phase equilibrium and imbalance, realize the extension based on optimization method
The purpose of trend solves the distributed network trend based on sensitive method reply partial voltage control device under fully loaded transportation condition
The problem of detecting under nonlinear problem and fully loaded transportation condition and relaxing the limitation of infeasible voltage, improves non-thread personality under fully loaded transportation condition
The convergence of sensitivity method, suitable for infeasible voltage limitation is detected and relaxed under fully loaded transportation condition.
In the description of the present invention, it is to be understood that, term " center ", " longitudinal direction ", " transverse direction ", " length ", " width ",
" thickness ", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom" "inner", "outside", " up time
The orientation or positional relationship of the instructions such as needle ", " counterclockwise ", " axial direction ", " radial direction ", " circumferential direction " be orientation based on the figure or
Positional relationship is merely for convenience of description of the present invention and simplification of the description, rather than the device or element of indication or suggestion meaning must
There must be specific orientation, be constructed and operated in a specific orientation, therefore be not considered as limiting the invention.
In addition, term " first ", " second " are used for descriptive purposes only and cannot be understood as indicating or suggesting relative importance
Or implicitly indicate the quantity of indicated technical characteristic.Define " first " as a result, the feature of " second " can be expressed or
Implicitly include at least one this feature.In the description of the present invention, the meaning of " plurality " is at least two, such as two, three
It is a etc., unless otherwise specifically defined.
In the present invention unless specifically defined or limited otherwise, term " installation ", " connected ", " connection ", " fixation " etc.
Term shall be understood in a broad sense, for example, it may be being fixedly connected, may be a detachable connection, or integral;It can be mechanical connect
It connects, is also possible to be electrically connected;It can be directly connected, can also can be in two elements indirectly connected through an intermediary
The interaction relationship of the connection in portion or two elements, unless otherwise restricted clearly.For those of ordinary skill in the art
For, the specific meanings of the above terms in the present invention can be understood according to specific conditions.
In the present invention unless specifically defined or limited otherwise, fisrt feature in the second feature " on " or " down " can be with
It is that the first and second features directly contact or the first and second features pass through intermediary mediate contact.Moreover, fisrt feature exists
Second feature " on ", " top " and " above " but fisrt feature be directly above or diagonally above the second feature, or be merely representative of
First feature horizontal height is higher than second feature.Fisrt feature can be under the second feature " below ", " below " and " below "
One feature is directly under or diagonally below the second feature, or is merely representative of first feature horizontal height less than second feature.
In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", " specifically show
The description of example " or " some examples " etc. means specific features, structure, material or spy described in conjunction with this embodiment or example
Point is included at least one embodiment or example of the invention.In the present specification, schematic expression of the above terms are not
It must be directed to identical embodiment or example.Moreover, particular features, structures, materials, or characteristics described can be in office
It can be combined in any suitable manner in one or more embodiment or examples.In addition, without conflicting with each other, the skill of this field
Art personnel can tie the feature of different embodiments or examples described in this specification and different embodiments or examples
It closes and combines.
Although the embodiments of the present invention has been shown and described above, it is to be understood that above-described embodiment is example
Property, it is not considered as limiting the invention, those skilled in the art within the scope of the invention can be to above-mentioned
Embodiment is changed, modifies, replacement and variant.
Claims (10)
1. a kind of extension trend method based on optimization method, which comprises the following steps:
Equation is limited according to the busbar voltage that limitation formula establishes relaxation;
The tide optimization problem for having partial voltage control device is established by optimization formula;
The Optimized model, the Optimized model of PV bus and the Optimized model of PQ bus of zero injection bus are established by constraint equation;
And
By the IPM of constant Hessian matrix to the Optimized model of the zero injection bus, the Optimized model of the PV bus and institute
The Optimized model for stating PQ bus is solved, to obtain integer tap position, and under the conditions of three-phase equilibrium and imbalance into
Row functional test.
2. the extension trend method according to claim 1 based on optimization method, which is characterized in that built by limiting formula
Found the busbar voltage limitation equation, the limitation formula are as follows:
Wherein,For the lower limit of busbar voltage square,For the upper limit of the busbar voltage square, l and u are non-negative pine
Relaxation variable, UreAnd UimThe respectively real and imaginary parts of busbar voltage complex variable, guFor the upper limit of function, g is mapping, glFor letter
Several lower limits,For busbar voltage square.
3. the extension trend method according to claim 2 based on optimization method, which is characterized in that public by following optimization
Formula establishes the tide optimization problem, the optimization formula are as follows:
Wherein, t is transformer tap positions, tminAnd tmaxThe minimum value of the respectively described transformer tap positions limits and maximum
Value limitation, h are mapping, and v is non-negative slack variable, and s.t. is constraint condition.
4. the extension trend method according to claim 1 based on optimization method, which is characterized in that described to establish zero injection
The Optimized model of the Optimized model of bus, the Optimized model of PV bus and PQ bus further comprises:
For the zero injection bus, the first constraint equation h is indicated with pairs of linear equation:
Wherein, IreAnd IimThe respectively real and imaginary parts of Injection Current, YzeroFor the admittance row corresponding to zero injection bus, U is
Multiple busbar voltage,WithThe respectively real and imaginary parts of equation.
For the PV bus, the second constraint equation h is indicated in pairs are as follows:
Wherein,Target is controlled for the voltage of PV bus, subscript PV is the respective vectors or row matrix of PV bus, and diag is pair
Angular moment battle array, * are the conjugation of complex variable, SreFor inject complex power real part,For the real part of equation, YPVFor corresponding to PV bus
Admittance row, UPVFor PV busbar voltage, U is multiple busbar voltage,For the real part of PV busbar voltage,For PV bus electricity
The imaginary part of pressure.
For the PQ bus, third constraint equation h is indicated in pairs are as follows:
Wherein, SreAnd SimThe real and imaginary parts of complex power are respectively injected,WithThe respectively real and imaginary parts of equation, diag
For diagonal matrix, * is the conjugation of complex variable, and subscript PQ is the respective vectors or row matrix of PQ bus, UPQFor PQ busbar voltage,
YPQFor the admittance row corresponding to PQ bus, U is multiple busbar voltage.
5. the extension trend method according to claim 1 based on optimization method, which is characterized in that taken out by rounding up
The fractional value of head position obtains the integer tap position.
6. a kind of extension trend device based on optimization method characterized by comprising
Establishing equation module is limited, the busbar voltage for establishing relaxation according to limitation formula limits equation;
Optimization problem establishes module, for establishing the tide optimization problem for having partial voltage control device by optimization formula;
Optimized model establishes module, for establishing the Optimized model of zero injection bus, the optimization mould of PV bus by constraint equation
The Optimized model of type and PQ bus;And
Trend expansion module, for by the IPM of constant Hessian matrix to the Optimized model of the zero injection bus, PV bus
The Optimized model of Optimized model and PQ bus is solved, to obtain integer tap position, and in three-phase equilibrium and imbalance
Under the conditions of carry out functional test.
7. the extension trend device according to claim 6 based on optimization method, it is characterised in that: the limitation equation is built
Formula is limited described in formwork erection block are as follows:
Wherein,For the lower limit of busbar voltage square,For the upper limit of the busbar voltage square, l and u indicate non-negative
Slack variable, UreAnd UimThe respectively real and imaginary parts of busbar voltage complex variable, guFor the function upper limit, g is mapping, glFor letter
Number lower limit,For busbar voltage square.
8. the extension trend device according to claim 7 based on optimization method, which is characterized in that the optimization problem is built
Optimize formula described in formwork erection block are as follows:
Wherein, t is transformer tap positions, tminAnd tmaxThe minimum value of the respectively described transformer tap positions limits and maximum
Value limitation, h are mapping, and v is non-negative slack variable, and s.t. is constraint condition.
9. the extension trend device according to claim 6 based on optimization method, which is characterized in that the Optimized model is built
Formwork erection block further comprises:
For the zero injection bus, the first constraint equation h is indicated with pairs of linear equation:
Wherein, IreAnd IimThe respectively real and imaginary parts of Injection Current, YzeroFor the admittance row corresponding to zero injection bus, U is
Multiple busbar voltage,WithThe respectively real and imaginary parts of equation.
For the PV bus, the second constraint equation h is indicated in pairs are as follows:
Wherein,Target is controlled for the voltage of PV bus, subscript PV is the respective vectors or row matrix of PV bus, and diag is pair
Angular moment battle array, * are the conjugation of complex variable, SreFor inject complex power real part,For the real part of equation, YPVFor corresponding to PV bus
Admittance row, UPVFor PV busbar voltage, U is multiple busbar voltage,For the real part of PV busbar voltage,For PV bus electricity
The imaginary part of pressure.
For the PQ bus, third constraint equation h is indicated in pairs are as follows:
Wherein, SreAnd SimThe real and imaginary parts of complex power are respectively injected,WithThe respectively real and imaginary parts of equation, diag
For diagonal matrix, * is the conjugation of complex variable, and subscript PQ is the respective vectors or row matrix of PQ bus, UPQFor PQ busbar voltage,
YPQFor the admittance row corresponding to PQ bus, U is multiple busbar voltage.
10. the extension trend device according to claim 6 based on optimization method, which is characterized in that by rounding up
The fractional value of tap position obtains the integer tap position.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711064949.8A CN109755942B (en) | 2017-11-02 | 2017-11-02 | Tidal current expanding method and device based on optimization method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711064949.8A CN109755942B (en) | 2017-11-02 | 2017-11-02 | Tidal current expanding method and device based on optimization method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109755942A true CN109755942A (en) | 2019-05-14 |
CN109755942B CN109755942B (en) | 2021-05-07 |
Family
ID=66397931
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711064949.8A Expired - Fee Related CN109755942B (en) | 2017-11-02 | 2017-11-02 | Tidal current expanding method and device based on optimization method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109755942B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2482417A2 (en) * | 2011-01-27 | 2012-08-01 | General Electric Company | Systems, Methods, and Apparatus for Accelerating Volt/VAR Load Flow Optimization |
CN103150606A (en) * | 2013-01-22 | 2013-06-12 | 中国电力科学研究院 | Optimal power flow optimization method of distributed power supplies |
CN103151784A (en) * | 2013-04-01 | 2013-06-12 | 国家电网公司 | AVC (automatic voltage control) system-based reactive voltage optimization method and device |
CN104092213A (en) * | 2014-07-30 | 2014-10-08 | 东南大学 | Power analyzing method for indeterminate power flow branches based on optimization method |
CN106253288A (en) * | 2016-07-19 | 2016-12-21 | 河海大学 | A kind of optimal load flow algorithm containing THE UPFC based on automatic differential |
CN106532710A (en) * | 2016-11-04 | 2017-03-22 | 国网山东省电力公司 | Microgrid power flow optimization method considering voltage stability constraint |
-
2017
- 2017-11-02 CN CN201711064949.8A patent/CN109755942B/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2482417A2 (en) * | 2011-01-27 | 2012-08-01 | General Electric Company | Systems, Methods, and Apparatus for Accelerating Volt/VAR Load Flow Optimization |
CN103150606A (en) * | 2013-01-22 | 2013-06-12 | 中国电力科学研究院 | Optimal power flow optimization method of distributed power supplies |
CN103151784A (en) * | 2013-04-01 | 2013-06-12 | 国家电网公司 | AVC (automatic voltage control) system-based reactive voltage optimization method and device |
CN104092213A (en) * | 2014-07-30 | 2014-10-08 | 东南大学 | Power analyzing method for indeterminate power flow branches based on optimization method |
CN106253288A (en) * | 2016-07-19 | 2016-12-21 | 河海大学 | A kind of optimal load flow algorithm containing THE UPFC based on automatic differential |
CN106532710A (en) * | 2016-11-04 | 2017-03-22 | 国网山东省电力公司 | Microgrid power flow optimization method considering voltage stability constraint |
Non-Patent Citations (1)
Title |
---|
雷婧婷 等: "含直流配电网的交直流潮流计算", 《中国电机工程学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109755942B (en) | 2021-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106233557B (en) | The control of microgrid | |
JP6140803B2 (en) | Photovoltaic power plant control technology | |
US20130099800A1 (en) | Algorithm and implementation system for measuring impedance in the d-q domain | |
CN103532367B (en) | interleaved PFC control method and device | |
CN106953317B (en) | The discrimination method of secondary/supersynchronous coupled impedance (admittance) model of power equipment | |
CN104137377B (en) | For by method Tong Bu with line voltage for supplying voltage | |
CN106130025B (en) | A kind of linear method that three-phase distribution net trend is calculated based on loop analysis | |
CN108365609B (en) | Unified iterative power flow calculation method considering automatic power generation control of interconnected power grid partitions | |
CN109802381A (en) | A kind of direct-current grid multi-source dynamic coordinate control method based on fuzzy control | |
CN109061338A (en) | Energy storage PCS is to dragging test platform and method | |
CN109462237A (en) | The method and system of meter and direct current layered structure alternating current-direct current mixed connection Load flow calculation | |
CN105046588B (en) | A kind of calculation method of the improvement direct current Dynamic Optimal Power Flow Problem based on network loss iteration | |
Macana et al. | Accurate proportional power sharing with minimum communication requirements for inverter-based islanded microgrids | |
CN104601029B (en) | Inverter parallel-control system droop control method | |
CN110927519A (en) | Active power distribution network fault positioning method based on mu PMU measurement value | |
CN109193633A (en) | Optimal configuration method for DC rated power in multi-DC feed-in system | |
CN110850175A (en) | Direct current micro-grid line impedance detection method and device based on low-frequency current injection | |
Cvetkovic et al. | Un-terminated, low-frequency terminal-behavioral dq model of three-phase converters | |
CN109755942A (en) | Extension trend method and device based on optimization method | |
Xue et al. | A comprehensive study on impedance models of grid-tied voltage-source converters | |
Marius et al. | Modeling and simuling power active filter using method of generalized reactive power theory | |
CN109245097B (en) | Method and system for calculating voltage coupling action factor of alternating current-direct current system based on extended Jacobian matrix | |
CN111256275B (en) | Operation control method and system, compressor and air conditioner | |
US20120016530A1 (en) | Method and system for managing power quality and efficiency | |
Wang et al. | Line impedance compensation control strategy for multiple interlinking converters in hybrid AC/DC microgrid |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210507 |
|
CF01 | Termination of patent right due to non-payment of annual fee |