CN105023058B - Power distribution network intelligence Sofe Switch running optimizatin method that is a kind of while considering switch motion - Google Patents

Power distribution network intelligence Sofe Switch running optimizatin method that is a kind of while considering switch motion Download PDF

Info

Publication number
CN105023058B
CN105023058B CN201510395612.XA CN201510395612A CN105023058B CN 105023058 B CN105023058 B CN 105023058B CN 201510395612 A CN201510395612 A CN 201510395612A CN 105023058 B CN105023058 B CN 105023058B
Authority
CN
China
Prior art keywords
node
switch
sofe switch
power
snop
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510395612.XA
Other languages
Chinese (zh)
Other versions
CN105023058A (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.)
Tianjin University
Hainan Power Grid Co Ltd
Original Assignee
Tianjin University
Hainan Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University, Hainan Power Grid Co Ltd filed Critical Tianjin University
Priority to CN201510395612.XA priority Critical patent/CN105023058B/en
Publication of CN105023058A publication Critical patent/CN105023058A/en
Application granted granted Critical
Publication of CN105023058B publication Critical patent/CN105023058B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

Power distribution network intelligence Sofe Switch running optimizatin method that is a kind of while considering switch motion:According to distribution system incoming line parameter, load level and network topology connection relation, system operation voltage level and branch current limitation, distributed generation resource on-position, type and capacity and parameter, intelligent Sofe Switch on-position and capacity and parameter, running optimizatin period internal loading and distributed generation resource operation characteristic prediction curve, system reference voltage and reference power initial value;Establish the timing optimization model of power distribution network interconnection switch and intelligent Sofe Switch synthetic operation;According to the canonical form of cone optimization in the timing optimization model of power distribution network interconnection switch and intelligent Sofe Switch synthetic operation object function and nonlinear restriction carry out Based On The Conic Model conversion;Object function is linearized respectively, nonlinear restriction is converted into linear restriction, second order cone constraint or rotating cone constraint, and obtained mathematical model is solved using CPLEX solvers.The invention avoids cumbersome iteration and a large amount of test.

Description

Power distribution network intelligence Sofe Switch running optimizatin method that is a kind of while considering switch motion
Technical field
The present invention relates to a kind of timing optimization methods of distribution network operation.Considering switch motion simultaneously more particularly to a kind of Power distribution network intelligence Sofe Switch running optimizatin method.
Background technology
Highest attention to the energy and environment so that the development of power distribution network is faced with new pressure and challenge, these pressure with It challenges while being also the opportunity for pushing conventional electrical distribution net to develop to intelligent distribution network.In intelligent distribution network, controllable device is increasingly Increase, network structure and the method for operation are more flexible changeable, and advanced Distribution Automation Technology, advanced Information and Communication Technology are able to Extensive use, distributed generation resource, energy storage, Demand-side resource etc. begin participating in the operation and optimization of power distribution network.Intelligent distribution network Development and the extensive use of new energy power generation technology are pushing the deep reform of adapted power mode and running fluidization air flow.
Intelligent Sofe Switch (Soft Normally Open Point, SNOP) device is exactly to derive in the above context A kind of novel intelligent power distribution equipment of the traditional interconnection switch of substitution.Compared with switching manipulation, the power control of SNOP is safer, Reliably, it might even be possible to realize real-time optimization, randomness and fluctuation that distributed generation resource and load are brought can be successfully managed.But It is that the realization of SNOP is based primarily upon full-control type power electronic device, and the cost of these devices itself is higher, in a short time in power distribution network Interconnection switch can not possibly be replaced completely by SNOP.This running optimizatin for allowing for power distribution network needs overall thinking contact to open It closes and SNOP and the case where deposit, optimal operation model will be one and need while solve discrete magnitude (on off state) and continuous quantity The mixed integer nonlinear programming problem of (SNOP transimission powers).
In intelligent distribution network, the distributed generation resource accessed extensively exacerbates the uncertainty of system operation, it is contemplated that opens The factors such as loss and dash current are closed, interconnection switch can not possibly frequently be cut-off, and traditional network reconfiguration is difficult to accomplish power distribution network Real-time adjustment, and SNOP can then change transimission power in real time, adjust operating status, to cope with band after distributed generation resource access Come the problems such as a series of voltage is out-of-limit, circuit overload, it is therefore desirable to from the angle of time series to power distribution network optimization problem into Row modeling, and need to consider interconnection switch and the coordination optimization problem of SNOP.After considering its temporal aspect, distribution network operation In optimization process, making rational planning for for switch motion just becomes a urgent problem to be solved.Meanwhile considering switch motion expense SNOP operation timing optimization problem can with when discontinuity surface number increase solve dimension sharply increase, become mix on a large scale it is whole Number nonlinear programming problem, causes its solution to become more difficult or even infeasible.
For solving this kind of extensive mixed integer nonlinear programming problem, it is also difficult to find at present a kind of quick, effective Method for solving.The solution of the problem has been proposed at present and has developed a variety of optimization methods, mainly have including:1) traditional Mathematics Optimization Method, including analytic method, successive elimination method etc.;2) heuritic approach, including Sensitivity Analysis Method, specially Family's system etc.;3) randomized optimization process, including genetic algorithm, particle cluster algorithm etc..
Although the above method or technology have certain application, also all there is clearly disadvantageous, as traditional mathematics are excellent Although change method can theoretically carry out global optimizing, it is inevitably present " dimension calamity " problem in practical application, calculates Explosive surge is often presented in time;Heuritic approach requires, there are one polynomial time, to calculate in terms of time complexity Speed is fast, but obtained optimal solution either lacks optimality in mathematical meaning or only locally optimal solution;Although random excellent The last solution that change method is searched is unrelated with initial solution, but the power distribution network of different scales is needed to reset its control ginseng Number, population quantity, iterations etc., to ensure to find globally optimal solution with larger probability.Heuristic and random device It is suitable for solving integer programming problem, but for considering that the power distribution network timing optimization of interconnection switch and SNOP synthetic operations is asked more Topic, mathematics is substantially extensive mixed integer nonlinear programming problem, so traditional mathematics optimization method, heuritic approach pair In on such issues that solve, speed or precision cannot mostly be met the requirements simultaneously.Therefore, it is necessary to one kind accurately, rapid solving is above-mentioned The model and algorithm of optimization problem.
Invention content
The technical problem to be solved by the invention is to provide one kind can consider switch motion expense it is equivalent loss, The power distribution networks running wastage such as SNOP running wastages and via net loss, while determining that rational switch motion sequential and SNOP are run Consider the power distribution network intelligence Sofe Switch running optimizatin method of switch motion.
The technical solution adopted in the present invention is:Power distribution network intelligence Sofe Switch operation that is a kind of while considering switch motion is excellent Change method, includes the following steps:
1) according to selected distribution system, incoming line parameter, load level and network topology connection relation, system operation Voltage level and branch current limitation, distributed generation resource on-position, type and capacity and parameter, intelligent Sofe Switch on-position With capacity and parameter, running optimizatin period internal loading and distributed generation resource operation characteristic prediction curve and system reference voltage With reference power initial value;
2) the distribution system structure and parameter provided according to step 1), while considering the switch of via net loss, network reconfiguration The running wastage of the equivalent loss of action expense and intelligent Sofe Switch, establishes power distribution network interconnection switch and intelligent Sofe Switch synthetic operation Timing optimization model, including:Selection root node is balance nodes, sets the minimum object function of distribution system running wastage, Network topology constraint, system load flow constraint, system operation constraint, intelligent Sofe Switch operation constraint are considered respectively;
3) according to the canonical form of cone optimization to the power distribution network interconnection switch and intelligence Sofe Switch collaboration fortune described in step 2) Object function and nonlinear restriction in capable timing optimization model carry out Based On The Conic Model conversion;
4) conversion Jing Guo step 3), respectively linearizes object function, and nonlinear restriction is converted into linear restriction, second order Cone constraint or rotating cone constraint, obtained mathematical model is solved using CPLEX solvers;
5) solving result of step 4), including on off state, intelligent Sofe Switch optimal transmission performance number, network trend are exported And target function value as a result.
The minimum object function of distribution system running wastage described in step 2), is expressed as
Min f=ES, loss+EL, loss+ESNOP, loss
In formula, the equivalent loss E of switch motion expenseS, loss, via net loss EL, lossWith the running wastage of intelligent Sofe Switch ESNOP, lossIt is indicated respectively with following formula
In formula, CSFor switch motion expense equivalent conversion coefficient;Δ t is the period interval that optimization calculates;NTIt is calculated for optimization When hop count, NNFor the node total number in system, NSNOPTo access the number of intelligent Sofe Switch in system;Ω (i) is node i The set of adjacent node;αij(t) on off state for being t period branches ij;rijFor the resistance of branch ij, Iij(t) it is the t periods Point i flows to the current amplitude of node j;PM, 1(t) and PM, 2(t) it is the active of two transverters of m-th of intelligent Sofe Switch of t periods Output power, AM, 1And AM, 2For the active loss coefficient of two transverters of m-th of intelligent Sofe Switch.
Network topology constraint representation described in step 2) is
αij(t)=βij(t)+βji(t)
αij(t) { 0,1 } ∈
βij(t) { 0,1 } ∈
In formula, NSFor the source node number in system;βij(t) indicate that t periods node i and node j relationships, node j are section It is 1 when the parent node of point i, is otherwise 0.
System load flow constraint representation described in step 2) is
Pi(t)=PDG, i(t)+PSNOP, i(t)-PLOAD, i(t)
Qi(t)=QDG, i(t)+QSNOP, i(t)-QLOAD, i(t)
In formula, Φ (i) is using node i as the branch headend node set of endpoint node, and Ψ (i) is using node i as headend node Branch set of end nodes;Ui(t) it is the voltage magnitude of t period node is, xijFor the reactance of branch ij;Pij(t) it is the t periods Node i flows to the active power of node j, Qij(t) reactive power of node j is flowed to for t period node is;Ω (i) is node i The set of adjacent node;Pi(t) it is the sum of the active power injected in t period node is, PDG, i(t)、PSNOP, i(t)、PLOAD, i(t) The active power of distributed generation resource injection, the active power of SNOP transmission, load consume active respectively in t period node is Power, Qi(t) it is the sum of the active power injected in t period node is, QDG, i(t)、QSNOP, i(t)、QLOAD, i(t) when being respectively t The reactive power of reactive power, load consumption that the reactive power of Duan Jiediani upper distributed generation resource injections, SNOP are sent out;M is one A maximum.
Intelligent Sofe Switch described in step 2) runs constraint representation
PM, 1(t)+PM, 2(t)+AM, 1|PM, 1(t)|+AM, 2|PM, 2(t) |=0
-QM, 1, max≤QM, 1(t)≤QM, 1, max
-QM, 2, max≤QM, 2(t)≤QM, 2, max
In formula, QM, 1(t) and QM, 2(t) it is the reactive power of two transverters output of m-th of intelligent Sofe Switch of t periods; SM, 1, max、SM, 2, max、QM, 1, max、QM, 2, maxIt the access capacity of respectively m-th intelligent two transverter of Sofe Switch and can be output The reactive power upper limit.
The power distribution network intelligence Sofe Switch running optimizatin method that is a kind of while considering switch motion of the present invention, based on solution Discontinuity surface interconnection switch and SNOP and the distribution network operation timing optimization deposited is asked when considering multiple under the premise of switch motion expense Topic, mathematics essence is extensive mixed integer nonlinear programming problem (MINLP), and current existing method is difficult to quick standard Really solve.Basic principle of the present invention according to cone optimization algorithm, bores the object function and constraints of Optimized model Conversion, converts former problem to MIXED INTEGER Second-order cone programming problem (MISOCP), greatly reduces solution difficulty, is easy to use Solution tool is solved.Cone optimization method of the present invention can carry out network reconfiguration and SNOP running optimizatin problems Unify legislation so that the problem of complicated mixed integer nonlinear programming, which solves, to be achieved, and cumbersome iteration and big is avoided The test of amount has in calculating speed and is significantly promoted.Also, because of graceful geometry and special place possessed by cone Reason mode can ensure the optimality of the solution of institute's Solve problems, apply it in the timing optimization problem of SNOP operations, Optimal system operation scheme can be quickly obtained.
Description of the drawings
Fig. 1 is the present invention while considering the flow chart of the power distribution network intelligence Sofe Switch running optimizatin method of switch motion;
Fig. 2 is 33 node examples of IEEE and distributed generation resource, the on-positions SNOP figure;
Fig. 3 is distributed generation resource and load operation Predicting Performance Characteristics curve;
Fig. 4 be switch motion expense conversion coefficient be 5kWh/ when it is each when discontinuity surface switch motion situation;
Fig. 5 a are that different switch motion expense conversion coefficients correspond to SNOP transmitting active power situations of change;
Fig. 5 b are that different switch motion expense conversion coefficients correspond to SNOP and send out reactive power situation of change;
Fig. 6 is different 18 voltage change situations of switch motion expense conversion coefficient corresponding node.
Specific implementation mode
With reference to embodiment and attached drawing to the power distribution network intelligence Sofe Switch that is a kind of while considering switch motion of the present invention Running optimizatin method is described in detail.
As shown in Figure 1, the power distribution network intelligence Sofe Switch running optimizatin method that is a kind of while considering switch motion of the present invention, Include the following steps:
1) according to selected distribution system, incoming line parameter, load level and network topology connection relation, system operation Voltage level and branch current limitation, distributed generation resource on-position, type and capacity and parameter, intelligent Sofe Switch (SNOP) connect Enter position and capacity and parameter, running optimizatin period internal loading and distributed generation resource operation characteristic prediction curve and system base The initial values such as quasi- voltage and reference power;
2) the distribution system structure and parameter provided according to step 1), while considering the switch of via net loss, network reconfiguration The running wastage of the equivalent loss of action expense and intelligent Sofe Switch, establishes power distribution network interconnection switch and intelligent Sofe Switch synthetic operation Timing optimization model, including:Selection root node is balance nodes, sets the minimum object function of distribution system running wastage, Network topology constraint, system load flow constraint, system operation constraint, intelligent Sofe Switch operation constraint are considered respectively;Wherein:
(1) the minimum object function of distribution system running wastage described in, is expressed as
Min f=ES, loss+EL, loss+ESNOP, loss (1)
In formula, the equivalent loss E of switch motion expenseS, loss, via net loss EL, lossWith the running wastage of intelligent Sofe Switch ESNOP, lossIt is indicated respectively with following formula
In formula, CSFor switch motion expense equivalent conversion coefficient;Δ t is the period interval that optimization calculates;NTIt is calculated for optimization When hop count, NNFor the node total number in system, NSNOPTo access the number of intelligent Sofe Switch in system;Ω (i) is node i The set of adjacent node;αij(t) on off state for being t period branches ij;rijFor the resistance of branch ij, Iij(t) it is the t periods Point i flows to the current amplitude of node j;PM, 1(t) and PM, 2(t) it is the active of two transverters of m-th of intelligent Sofe Switch of t periods Output power, AM, 1And AM, 2For the active loss coefficient of two transverters of m-th of intelligent Sofe Switch.
(2) the network topology constraint representation described in is
αij(t)=βij(t)+βji(t) (5)
αij(t) { 0,1 } (8) ∈
βij(t) { 0,1 } (9) ∈
In formula, NSFor the source node number in system;βij(t) indicate that t periods node i and node j relationships, node j are section It is 1 when the parent node of point i, is otherwise 0.
(3) the system load flow constraint representation described in is
Pi(t)=PDG, i(t)+PSNOP, i(t)-PLOAD, i(t) (12)
Qi(t)=QDG, i(t)+QSNOP, i(t)-QLOAD, i(t) (13)
In formula, Φ (i) is using node i as the branch headend node set of endpoint node, and Ψ (i) is using node i as headend node Branch set of end nodes;Ui(t) it is the voltage magnitude of t period node is, xijFor the reactance of branch ij;Pij(t) it is the t periods Node i flows to the active power of node j, Qij(t) reactive power of node j is flowed to for t period node is;Ω (i) is node i The set of adjacent node;Pi(t) it is the sum of the active power injected in t period node is, PDG, i(t)、PSNOP, i(t)、PLOAD, i(t) The active power of distributed generation resource injection, the active power of SNOP transmission, load consume active respectively in t period node is Power, Qi(t) it is the sum of the active power injected in t period node is, QDG, i(t)、QSNOP, i(t)、QLOAD, i(t) when being respectively t The idle work(of reactive power, load consumption that the reactive power of Duan Jiediani upper distributed generation resource injections, intelligent Sofe Switch are sent out Rate;M is a maximum.
(4) the system operation constraint representation described in is
-Mαij(t)≤Pii(t)≤Mαij(t) (19)
-Mαij(t)≤Qij(t)≤Mαij(t) (20)
In formula, UI, minAnd UI, maxThe respectively minimum allowable voltage value and maximum allowable voltage of node i;IIj, maxFor this The maximum allowed current value of branch.
(5) the intelligent Sofe Switch described in runs constraint representation
PM, 1(t)+PM, 2(t)+AM, 1|PM, 1(t)|+AM, 2|PM, 2(t) |=0 (22)
-QM, 1, max≤QM, 1(t)≤QM, 1, max (25)
-QM, 2, max≤QM, 2(t)≤QM, 2, max (26)
In formula, QM, 1(t) and QM, 2(t) it is the reactive power of two transverters output of m-th of intelligent Sofe Switch of t periods; SM, 1, max、SM, 2, max、QM, 1, max、QM, 2, maxIt the access capacity of respectively m-th intelligent two transverter of Sofe Switch and can be output The reactive power upper limit.
3) according to the canonical form of cone optimization to the power distribution network interconnection switch and intelligence Sofe Switch collaboration fortune described in step 2) Object function and nonlinear restriction in capable timing optimization model carry out Based On The Conic Model conversion, and specific method for transformation is as follows:
M0(t)=| αii(t)-αij(t-1) |=max { αij(t)-αij(t-1), αij(t-1)-αij(t) }, and increase constraint
M0(t)≥0 (27)
M0(t)≥αij(t)-αij(t-1) (28)
M0(t)≥αij(t-1)-αij(t) (29)
(2) object function via net lossWith constraints (10), (11), contain quadratic term in (14)~(18) and (21)WithUsing U2, i(t) and I2, ij(t) quadratic term is replacedWithIt is linearized.
(3) running wastage of object function intelligence Sofe Switch
Intelligent Sofe Switch operation constraints PM, 1(t)+PM, 2(t)+AM, 1|PM, 1(t)|+AM, 2|PM, 2(t) | in=0 (22) Contain absolute value term | PM, 1(t) | and | PM, 2(t) |, introduce auxiliary variable M1(t)=| PM, 1(t) |=max { PM, 1(t) ,-PM, 1 And M (t) }2(t)=| pM, 2(t) |=max { PM, 2(t) ,-PM, 2(t) }, and increase constraint
M1(t)≥0 (29)
M2(t)≥0 (30)
M1(t)≥PM, 1(t) (31)
M1(t)≥-PM, 1(t) (32)
M2(t)≥PM, 2(t) (33)
M2(t)≥-PM, 2(t) (34)
(4) system load flow constrainsFor the nonlinear restriction after above-mentioned steps are replaced, by it It is loose to be constrained for second order cone
||[2Pij(t) 2Qij(t) I2, ij(t)-U2, i(t)]T||2≤I2, ij(t)-U2, i(t) (35)
For nonlinear restriction, it is converted into rotating cone constraint
4) conversion Jing Guo step 3), respectively linearizes object function, and nonlinear restriction is converted into linear restriction, second order Cone constraint or rotating cone constraint, obtained mathematical model is solved using CPLEX solvers;
5) solving result of step 4), including on off state, intelligent Sofe Switch optimal transmission performance number, network trend are exported And target function value as a result.
While on off state, intelligent Sofe Switch optimal transmission power and trend being realized the present invention is based on cone optimization algorithm It solves.Establish the timing optimization problem of the power distribution network interconnection switch for considering switch motion expense and intelligent Sofe Switch synthetic operation Mathematical model, not only from it is single when discontinuity surface consider switch motion and the operation constraint of intelligent Sofe Switch, but also consider phase When adjacent between discontinuity surface switch change continuity and sequential relationship.
Specific example is given below:
For the present embodiment, the impedance value of circuit element in 33 node systems of IEEE as shown in Figure 2 is inputted first, is born Active power, the reactive power of lotus element, network topology connection relation, detail parameters are shown in Tables 1 and 2;Then 5 typhoons electricity is set The on-position of unit be node 10,16,17,30,33, access capacity be respectively 500kVA, 300kVA, 200kVA, 200kVA, 300kVA, the on-positions of 3 photovoltaic systems are node 7,13,27, and access capacity is respectively 500kVA, 300kVA, 400kVA, Power factor is 1.0;One group of SNOP access power distribution network is set again, replaces interconnection switch TS1, the capacity of two transverters equal For 500kVA, it is 200kVar that reactive power, which exports the upper limit,;Then, as unit of day, with 1 hour for time interval, using negative Lotus prediction technique simulates the day operation curve of load and wind-powered electricity generation, photovoltaic, as shown in Figure 3;The finally benchmark electricity of setting system Pressure is 12.66kV, reference power 1MVA, and maximum M takes 9999.
1 IEEE33 nodes example load on-position of table and power
2 IEEE33 node example line parameter circuit values of table
Discontinuity surface distich network switch and intelligent Sofe Switch and when the power distribution network deposited carries out when the present embodiment with 1 hour was one Sequence optimizes, and optimum results are shown in Table 3, and the switch motion situation of scheme 2 is as shown in Figure 4.
The different switch fare paths optimum results of table 3
It is Intel (R) Xeon (R) CPU E5-1620 to execute the computer hardware environment that optimization calculates, and dominant frequency is 3.70GHz inside saves as 32GB;Software environment is 7 operating systems of Windows.
Prioritization scheme considers different switch motion expense conversion coefficients, and to interconnection switch and intelligent Sofe Switch and that deposits match Operation of power networks carries out timing optimization, and considers the loss generated during intelligent Sofe Switch transimission power, intelligent Sofe Switch transmission Active power and go out the reactive power sent out in node 22 and see Fig. 5.With the increase of switch motion expense conversion coefficient, switch Action frequency significantly reduces, and helps to improve the service life of switch.On the other hand, network reconfiguration and intelligent Sofe Switch are run excellent The operation voltage level of system can be improved to a certain extent by changing, as shown in fig. 6, further improving power quality, improving confession Electric reliability.
The mathematics essence of power distribution network interconnection switch and the timing optimization problem of intelligent Sofe Switch synthetic operation is extensive mixed Integral nonlinear program-ming problem is closed, current existing optimization method can not solve this problem mostly.One kind proposed by the present invention The power distribution network intelligence Sofe Switch running optimizatin method for considering switch motion simultaneously, can fast and accurately solve problems, and It can guarantee the optimality of solution.For scheme two, at the same using a kind of hybrid solution method based on simulated annealing and cone optimization into Row solves, and is compared to the optimality and calculated performance of solution, and comparison result is shown in Table 4.
The different method for solving calculated performances of table 4 compare
/ The method of the present invention Method for mixing and optimizing
Target function value (member) 463.08 485.31
Solve the time (s) 27.12 1429.68

Claims (5)

1. power distribution network intelligence Sofe Switch running optimizatin method that is a kind of while considering switch motion, which is characterized in that including as follows Step:
1) according to selected distribution system, incoming line parameter, load level and network topology connection relation, system operation voltage The limitation of horizontal and branch current, distributed generation resource on-position, type and capacity and parameter, intelligent Sofe Switch on-position and appearance Amount and parameter, running optimizatin period internal loading and distributed generation resource operation characteristic prediction curve and system reference voltage and base Quasi- power initial value;
2) the distribution system structure and parameter provided according to step 1), while considering the switch motion of via net loss, network reconfiguration The running wastage of the equivalent loss of expense and intelligent Sofe Switch, establish power distribution network interconnection switch and intelligent Sofe Switch synthetic operation when Sequence Optimized model, including:Selection root node is balance nodes, sets the minimum object function of distribution system running wastage, respectively Consider network topology constraint, system load flow constraint, system operation constraint, intelligent Sofe Switch operation constraint;
3) according to the canonical form of cone optimization to described in step 2) power distribution network interconnection switch and intelligent Sofe Switch synthetic operation Object function and nonlinear restriction in timing optimization model carry out Based On The Conic Model conversion;
4) conversion Jing Guo step 3), respectively linearizes object function, and nonlinear restriction is converted into linear restriction, second order cone about Beam or rotating cone constraint, obtained mathematical model is solved using CPLEX solvers;
5) solving result of step 4), including on off state, intelligent Sofe Switch optimal transmission performance number, network power flow solutions are exported And target function value.
2. power distribution network intelligence Sofe Switch running optimizatin method that is according to claim 1 a kind of while considering switch motion, It is characterized in that, the minimum object function of distribution system running wastage described in step 2), is expressed as
Min f=ES, loss+EL, loss+ESNOP, loss
In formula, the equivalent loss E of switch motion expenseS, loss, via net loss EL, lossWith the running wastage E of intelligent Sofe SwitchSNOP, loss It is indicated respectively with following formula
In formula, CSFor switch motion expense equivalent conversion coefficient;Δ t is the period interval that optimization calculates;NTFor optimization calculate when Hop count, NNFor the node total number in system, NSNOPTo access the number of intelligent Sofe Switch in system;Ω (i) is the adjacent of node i The set of node;αij(t) on off state for being t period branches ij;rijFor the resistance of branch ij, Iij(t) it is t period node i streams To the current amplitude of node j;PM, 1(t) and PM, 2(t) it is the active output of two transverters of m-th of intelligent Sofe Switch of t periods Power, AM, 1And AM, 2For the active loss coefficient of two transverters of m-th of intelligent Sofe Switch.
3. power distribution network intelligence Sofe Switch running optimizatin method that is according to claim 1 a kind of while considering switch motion, It is characterized in that, the network topology constraint representation described in step 2) is
αij(t)=βij(t)+βji(t)
αij(t) { 0,1 } ∈
βij(t) { 0,1 } ∈
In formula, αij(t) on off state for being t period branches ij;Ω (i) is the set of the adjacent node of node i;NNFor in system Node total number;NSFor the source node number in system;βij(t) indicate that t periods node i and node j relationships, node j are node i Parent node when be 1, be otherwise 0.
4. power distribution network intelligence Sofe Switch running optimizatin method that is according to claim 1 a kind of while considering switch motion, It is characterized in that, the system load flow constraint representation described in step 2) is
Pi(t)=PDG, i(t)+PSNOP, i(t)-PLOAD, i(t)
Qi(t)=QDG, i(t)+QSNOP, i(t)-QLOAD, i(t)
In formula, Φ (i) is using node i as the branch headend node set of endpoint node, and Ψ (i) is using node i as the branch of headend node Road set of end nodes;Ui(t) it is the voltage magnitude of t period node is, xijFor the reactance of branch ij;Pij(t) it is t period node is Flow to the active power of node j, Qij(t) reactive power of node j is flowed to for t period node is;Ω (i) is the adjacent segments of node i The set of point;Pi(t) it is the sum of the active power injected in t period node is, PDG, i(t)、PSNOP, i(t)、PLOAD, i(t) it is respectively The active power of distributed generation resource injects in t period node is active power, SNOP transmission, the active power of load consumption, Qi (t) it is the sum of the active power injected in t period node is, QDG, i(t)、QSNOP, i(t)、QLOAD, i(t) it is respectively t period node is The reactive power of reactive power, load consumption that the reactive power of upper distributed generation resource injection, SNOP are sent out;M is one very big Value;rijFor the resistance of branch ij;Iij(t) current amplitude of node j is flowed to for t period node is;αij(t) it is t period branches ij On off state.
5. power distribution network intelligence Sofe Switch running optimizatin method that is according to claim 1 a kind of while considering switch motion, It is characterized in that, the intelligent Sofe Switch operation constraint representation described in step 2) is
PM, 1(t)+PM, 2(t)+AM, 1|PM, 1(t)|+AM, 2|PM, 2(t) |=0
-QM, 1, max≤QM, 1(t)≤QM, 1, max
-QM, 2, max≤QM, 2(t)≤QM, 2, max
In formula, QM, 1(t) and QM, 2(t) it is the reactive power of two transverters output of m-th of intelligent Sofe Switch of t periods; SM, 1, max、SM, 2, max、QM, 1, max、QM, 2, maxIt the access capacity of respectively m-th intelligent two transverter of Sofe Switch and can be output The reactive power upper limit;PM, 1(t) and PM, 2(t) it is the active output work of two transverters of m-th of intelligent Sofe Switch of t periods Rate;AM, 1And AM, 2For the active loss coefficient of two transverters of m-th of intelligent Sofe Switch.
CN201510395612.XA 2015-07-07 2015-07-07 Power distribution network intelligence Sofe Switch running optimizatin method that is a kind of while considering switch motion Active CN105023058B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510395612.XA CN105023058B (en) 2015-07-07 2015-07-07 Power distribution network intelligence Sofe Switch running optimizatin method that is a kind of while considering switch motion

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510395612.XA CN105023058B (en) 2015-07-07 2015-07-07 Power distribution network intelligence Sofe Switch running optimizatin method that is a kind of while considering switch motion

Publications (2)

Publication Number Publication Date
CN105023058A CN105023058A (en) 2015-11-04
CN105023058B true CN105023058B (en) 2018-09-18

Family

ID=54413011

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510395612.XA Active CN105023058B (en) 2015-07-07 2015-07-07 Power distribution network intelligence Sofe Switch running optimizatin method that is a kind of while considering switch motion

Country Status (1)

Country Link
CN (1) CN105023058B (en)

Families Citing this family (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105449713B (en) * 2015-12-11 2017-11-03 天津大学 Consider the intelligent Sofe Switch planing method of active power distribution network of distributed power source characteristic
CN105719196B (en) * 2016-01-18 2020-01-07 天津大学 Active power distribution network voltage reactive power control method based on intelligent soft switch
CN105896537B (en) * 2016-06-21 2018-05-01 中国南方电网有限责任公司电网技术研究中心 Power distribution network power supply recovery method based on intelligent soft switch
CN105977934B (en) * 2016-06-21 2018-04-06 天津大学 A kind of power distribution network intelligence Sofe Switch service restoration method for considering load importance
CN106329523B (en) * 2016-11-19 2018-09-25 中国南方电网有限责任公司电网技术研究中心 Active power distribution network intelligent soft switch robust optimization modeling method considering uncertainty
CN106655177B (en) * 2017-01-18 2018-12-28 中国南方电网有限责任公司电网技术研究中心 Distributed power supply maximum access capacity calculation method based on extended second-order cone programming
CN106655227B (en) * 2017-01-18 2019-02-15 天津大学 A kind of active power distribution network feeder line balancing method of loads based on intelligent Sofe Switch
CN107644118B (en) * 2017-08-04 2020-08-14 天津大学 Intelligent power distribution soft switch time sequence optimization method integrating energy storage
CN107451743B (en) * 2017-08-04 2021-03-19 天津大学 Active power distribution network island operation method considering source load operation characteristics
CN107341623A (en) * 2017-08-31 2017-11-10 国家电网公司 A kind of active power distribution network source storage lotus islet operation method of meter and network reconfiguration
CN107394775A (en) * 2017-08-31 2017-11-24 国家电网公司 A kind of active power distribution network sequential islet operation method of meter and load priority
CN107994577B (en) * 2017-12-29 2020-09-01 天津大学 Power distribution network dynamic reconstruction method considering switch short-time continuous action constraint
CN108649582A (en) * 2018-06-11 2018-10-12 中国石油大学(华东) A kind of distribution looped network optimal power flow control strategy based on SNOP
CN110783911B (en) * 2018-07-31 2021-04-02 上海交通大学 Soft intelligent switch configuration method and system for medium and low voltage distribution network
CN109492794B (en) * 2018-10-08 2021-08-31 国网天津市电力公司电力科学研究院 Electric vehicle charging station optimal configuration modeling method based on two-step equivalence method
CN109449945B (en) * 2018-11-22 2020-11-24 国网四川省电力公司天府新区供电公司 10kV loop closing current constraint power distribution network transfer optimization method based on DG (distributed generation) adjustment
US10971934B2 (en) 2018-12-31 2021-04-06 Abb Schweiz Ag Distribution networks with flexible direct current interconnection system
US11121543B2 (en) 2018-12-31 2021-09-14 Abb Schweiz Ag Fault mitigation in medium voltage distribution networks
CN109713670B (en) * 2019-01-28 2022-03-15 合肥工业大学 Fault recovery control optimization method based on three-port flexible multi-state switch
CN110021966A (en) * 2019-03-07 2019-07-16 华中科技大学 A kind of active distribution network Optimization Scheduling considering dynamic network reconfiguration
US10819112B1 (en) 2019-03-27 2020-10-27 Abb Schweiz Ag Feeder line fault response using direct current interconnection system
US11031773B2 (en) 2019-03-27 2021-06-08 Abb Power Grids Switzerland Ag Transformer isolation response using direct current link
CN110034560A (en) * 2019-04-23 2019-07-19 南京工程学院 Power distribution network multi-target optimum operation method based on intelligent Sofe Switch
CN110391673B (en) * 2019-07-03 2020-12-15 浙江大学 Multi-period active distribution network reconstruction method considering demand side response under high-proportion wind power
CN110490376B (en) * 2019-08-05 2023-04-07 天津大学 Intelligent soft switch planning method for improving reliability and economy of power distribution network
CN110445180A (en) * 2019-09-03 2019-11-12 国网河南省电力公司信阳供电公司 A kind of new-energy grid-connected method
CN110932259B (en) * 2019-11-21 2020-11-10 国网江苏省电力有限公司南通供电分公司 Distribution network reconstruction method considering multi-scene model and demand side response strategy
CN110991727A (en) * 2019-11-28 2020-04-10 海南电网有限责任公司 Power grid planning method based on power flow network loss model and line constraint model
CN111523213A (en) * 2020-04-15 2020-08-11 南京清然能源科技有限公司 Artificial intelligence-based optimized energy supply method for electricity core type comprehensive energy system
CN111541248B (en) * 2020-06-11 2021-10-15 南方电网科学研究院有限责任公司 Intelligent soft switch and energy storage system combined optimization method and device
CN113078678B (en) * 2021-04-13 2022-09-13 湖南大学 DG and SOP active-reactive collaborative planning method and device
CN118281858B (en) * 2023-03-17 2024-09-13 山东华科信息技术有限公司 Method for obtaining optimization result of active power distribution network
CN116780638A (en) * 2023-08-22 2023-09-19 国网天津市电力公司城南供电分公司 Snowflake power distribution network operation optimization method and device with soft switch and distributed energy storage

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103400208A (en) * 2013-08-01 2013-11-20 天津大学 Power distribution network distributive power supply optimal access capacity determining method based on cone optimization

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103400208A (en) * 2013-08-01 2013-11-20 天津大学 Power distribution network distributive power supply optimal access capacity determining method based on cone optimization

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"An evaluation method for power supply capacity of distribution system";Wang cheng shan,etc..;《Electric Utility Deregulation and Restructuring and Power Technologies》;20080409;1244-1248 *
"智能软开关非接触电源研究";周成虎;《计算机测量与控制》;20141130;630-633 *
"智能配电系统的发展与挑战";王成山;《供用电》;20150228;43-44 *

Also Published As

Publication number Publication date
CN105023058A (en) 2015-11-04

Similar Documents

Publication Publication Date Title
CN105023058B (en) Power distribution network intelligence Sofe Switch running optimizatin method that is a kind of while considering switch motion
CN104794541B (en) A kind of power distribution network running optimizatin method based on simulated annealing and cone optimization
CN103400207B (en) Operation optimization method for power distribution network comprising schedulable distributed power supply
CN107579518B (en) Power system environment economic load dispatching method and apparatus based on MHBA
CN108134401B (en) Multi-target power flow optimization and control method for alternating current-direct current hybrid system
CN108599154A (en) A kind of three-phase imbalance power distribution network robust dynamic reconfiguration method considering uncertain budget
CN112039069A (en) Double-layer collaborative planning method and system for power distribution network energy storage and flexible switch
CN105449713A (en) Active power distribution network soft normally open point (SNOP) planning method considering distributed generator characteristics
CN107563550A (en) A kind of Optimal Configuration Method of the real-time distributed state estimation of power distribution network based on PMU and PMU
CN103279661B (en) Substation capacity Optimal Configuration Method based on Hybrid quantum inspired evolution algorithm
CN111276968A (en) Singular perturbation-based distributed convergence control method and system for comprehensive energy system
CN112615368A (en) Distributed sensing-based platform area topology generation method and sensing terminal device
CN106684932B (en) A kind of interconnected network economic load dispatching method considering electric network composition adjustment
CN109888817B (en) Method for carrying out position deployment and capacity planning on photovoltaic power station and data center
Pedersen et al. Network condition based adaptive control and its application to power balancing in electrical grids
CN110675276B (en) Method and system for inversion droop control of direct current power transmission system
CN104821580A (en) Three-phase reactive power control partitioning method based on reactive power source clustering analysis
CN115360768A (en) Power scheduling method and device based on muzero and deep reinforcement learning and storage medium
CN106452837A (en) Intelligent substation switch dynamic configuration method
CN107425519B (en) Method for calculating maximum power supply capacity of three-phase power distribution network containing distributed power supply
CN109038640B (en) Capacity determining method for UPFC serial side converter based on phase diagram method
CN112835647A (en) Application software configuration method, configuration device, storage medium and electronic device
Luo et al. TanGo: A cost optimization framework for tenant task placement in geo-distributed clouds
Wu et al. Research on the Influence of Communication Time-Delay on Multi-agent System Based Distributed Energy Dispatching Control Consensus
CN115086435B (en) Intelligent energy comprehensive distribution control system

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