CN109038551A - The MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution - Google Patents
The MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution Download PDFInfo
- Publication number
- CN109038551A CN109038551A CN201810804522.5A CN201810804522A CN109038551A CN 109038551 A CN109038551 A CN 109038551A CN 201810804522 A CN201810804522 A CN 201810804522A CN 109038551 A CN109038551 A CN 109038551A
- Authority
- CN
- China
- Prior art keywords
- power
- period
- node
- transmitting electricity
- transmission
- 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
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses the MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution, this method is on the basis of having formulated the whole network generating set Active Generation operation plan and DC power transmission line active power transmission plan, by boring the processing methods such as relaxation converts and large M is equivalent, to with nonlinearity objective function and constraint condition carry out convexification processing, convert the Optimized model of foundation to the MIXED INTEGER Convex Programming Model for being easy to solve, the AC and DC transmission channel power optimized allocation of acquisition can be effectively reduced the loss of AC and DC passway for transmitting electricity, with apparent economic benefit, and it ensure that the optimality of gained solution mathematically, improve the solution efficiency that optimization calculates, to realize the performance driving economy of AC-DC interconnecting power network.
Description
Technical field
The present invention relates to technical field of electric power, and in particular to a kind of AC/DC parallel passway for transmitting electricity power optimization distribution mixes
Close integer method of convex programming.
Background technique
Due to being unevenly distributed for the energy and load, remote, large capacity, ultra-high/extra-high voltage transmission of electricity becomes the hair of China's power grid
Exhibition trend, the regional power grids such as south electric network, East China Power Grid have all formd large-scale AC-DC interconnecting power network pattern.Currently, south
Power grid " transferring electricity from the west to the east " channel covers southern five provinces (area), since Yunnan, Guizhou Province are within the border, passes through Guangxi whole area, reaches wide
Eastern domestic, the quantity of electricity in a large amount of Yunnan and Guizhou is sent to Guangdong every year by the nearly 2000km of west and east span, to southern five provinces and regions electricity
The optimization of power resource and reasonable disposition and the socio-economic development of each provinces and regions play an important role.With ± 800 kilovolts of Yunnan
Northwest extra-high voltage direct-current transmission engineering is open to traffic, and south electric network transferring electricity from the west to the east main channel has formd the network of " eight friendships ten are straight "
Structure.Since there is some difference for the transmission loss rate of each communication channel and direct current channel, pass through Optimum operation of power networks
Transmission power and other operation of power networks states of mode, each direct current of reasonable arrangement and communication channel etc., can be effectively reduced
The loss of electricity of entire transmission system, brings significant economic benefit.Therefore, meeting system safety in operation and power quality
Under conditions of it is required that, how according to each province net between total electricity electricity exchange plan, reasonable distribution save net between each direct current with exchange
The power transmission power planning of passway for transmitting electricity is the key technology of a urgent need to resolve to reduce the power loss of entire passway for transmitting electricity
Problem.
AC/DC parallel passway for transmitting electricity power optimization assignment problem is that the corresponding unit inside the province worked out in known each province net is sent out
Under the premise of always exchanging power planning between electric power output plan and each province net, the transmission on AC-DC tie line road between each province is determined
Power planning, to ensure under the premise of the transimission power security constraint of alternating current interconnection section between meeting each province, so that whole
The sum of loss of each transmission line of alternation current and DC power transmission line is minimum on a passway for transmitting electricity.Currently, AC/DC parallel is transmitted electricity
The power-balance constraint of each node generally uses conventional ac and dc systems Load flow calculation mould in channel power model of optimizing allocation
Type, however conventional ac and dc systems power flow algorithm is nonlinear model;Furthermore, it is contemplated that the change of current of DC converter station becomes
Depressor is the on-load regulator transformer of adjustable no-load voltage ratio, is related to the integer variable of discrete gear, and in practical engineering applications,
Realize that DC line power planning is completed by adjusting power electronic equipment, too frequent operation is unfavorable to equipment, because
What the adjusting number of this DC line power was limited by within a certain period of time, and each transmission power shape must be limited
The minimum duration of state, thus require obtain direct current transportation channel transmission of electricity it is contemplated that ladder shape curve.These
Restrictive condition causes AC/DC parallel passway for transmitting electricity power optimization distribution model to belong to mixed-integer nonlinear programming model, solves
Difficulty is big, it is difficult to obtain globally optimal solution, and solving speed is slow.
In practical engineering application, it is desirable that the dc power Plan Curve finally made should be stepped, and mesh
The optimizing distribution method of preceding existing AC/DC parallel passway for transmitting electricity power is directly logical to direct current transportation not in optimization calculates
The ladderization processing of road transmission of electricity plan, does not consider that first solving optimization model obtains AC/DC parallel passway for transmitting electricity power often
Optimize plan of distribution, then ladder processing, the power optimization obtained in this way are carried out to obtained direct current transportation channel transmission of electricity plan
The degree of optimization of allocation plan is often poor, and the economy of operation is also poor.
In addition, the power-balance constraint of each node generally uses in AC/DC parallel passway for transmitting electricity power optimization distribution model
Conventional ac and dc systems power flow algorithm, is related to the operation of trigonometric function and quadratic function, is Nonlinear programming Model, examines
The adjusting for considering the converter power transformer no-load voltage ratio of DC converter station is the integer variable of discrete gear, and it is logical to need to introduce reflection direct current
The integer variable whether the transmission power state in road is adjusted.Therefore, the optimization distribution of AC/DC parallel passway for transmitting electricity power
Belong to mixed-integer nonlinear programming model, solution difficulty is big, using the common mixed integer nonlinear programming solver such as SBB
It is solved and often hardly results in optimal solution and be difficult to ensure the optimality of gained solution mathematically, and solving speed is slow.
Above-mentioned technology has the drawback that: conventional AC/DC parallel passway for transmitting electricity power optimization distributes calculation method, does not have
Have in view of directly being handled the ladderization of direct current transportation channel transmission of electricity plan in optimization calculates, and it is right again after optimization calculates
The transmission of electricity plan in direct current transportation channel carries out ladder processing, and the degree of optimization of the power optimization allocation plan obtained in this way is often
It is poor;It is introduced directly into the optimization distribution of the AC/DC parallel passway for transmitting electricity power of direct current transportation channel transmission of electricity plan ladder requirement
Model is related to the operation of trigonometric function and quadratic function, and is related to the discrete gear of converter power transformer and dc power adjustment state
Integer variable, belong to mixed-integer nonlinear programming model, it is big to solve difficulty, non-thread using the common MIXED INTEGER such as SBB
Property programming evaluation device solved and often hardly result in optimal solution and be difficult to ensure gained solution optimality mathematically, Er Qieqiu
It is slow to solve speed.
Summary of the invention
It is an object of the invention to overcome above-mentioned the deficiencies in the prior art, a kind of AC/DC parallel passway for transmitting electricity power is provided
Optimize the MIXED INTEGER method of convex programming of distribution, to realize the performance driving economy of AC-DC interconnecting power network.
To achieve the above object, the technical scheme is that
The MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution, including
Total exchange power planning between sending power data and each province to net under each exchange node is obtained, is handed over directly with establishing
Passway for transmitting electricity power optimization distribution model in parallel is flowed, the AC/DC parallel passway for transmitting electricity power optimization distribution model is entirely to transmit electricity
The minimum target of the sum of loss of each transmission line of alternation current and DC power transmission line on channel, and constraint condition is set;
By cone relaxation conversion and large M equivalent way, mould is distributed to AC/DC parallel passway for transmitting electricity power optimization
Non-linear partial in the objective function and constraint condition of type carries out convexification processing, and the AC/DC parallel passway for transmitting electricity power is excellent
Change distribution model and is converted into MIXED INTEGER Convex Programming Model;
MIXED INTEGER Convex Programming Model is solved, AC/DC parallel passway for transmitting electricity power optimization point between each province is obtained
With scheme;
AC/DC parallel passway for transmitting electricity power optimization allocation plan between obtained each province is handed down to each converter station
Operations staff executes.
Compared with prior art, the present invention the beneficial effect is that:
This method is formulating the whole network generating set Active Generation operation plan and the active power transmission meter of DC power transmission line
On the basis of drawing, by boring the processing methods such as relaxation converts and large M is equivalent, to the objective function peace treaty with nonlinearity
Beam condition carries out convexification processing, converts the Optimized model of foundation to the MIXED INTEGER Convex Programming Model for being easy to solve, acquisition
AC and DC transmission channel power optimized allocation can be effectively reduced the loss of alternating current circuit and DC line, have apparent
Economic benefit, and ensure that the optimality of gained solution mathematically, the solution efficiency that optimization calculates is improved, to realize alternating current-direct current
The performance driving economy of interconnected network.
Detailed description of the invention
Fig. 1 is AC-DC interconnecting power network main grid structure structure chart;
Fig. 2 is the total losses change curve in optimization front and back direct current transportation channel;
Fig. 3 is the total losses change curve in optimization front and back ac transmission channel;
Fig. 4 a is optimized allocation figure compared with the corresponding 4 direct current channel power of case in front of optimization, in figure, dotted line generation
Before table optimization, after solid line representing optimized;
Fig. 4 b is optimized allocation figure compared with the corresponding other 3 direct current channel power of optimization front case, in figure, void
Before line representing optimized, after solid line representing optimized.
Specific embodiment
The contents of the present invention are described in further details with reference to the accompanying drawings and detailed description.
Embodiment:
The MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution, specifically includes:
It is obtained from the O&M automated system of EHV transmission company and send power data and each under each exchange node
Total exchange power planning between net is saved, to establish AC/DC parallel passway for transmitting electricity power optimization distribution model, the AC/DC parallel
Passway for transmitting electricity power optimization distribution model is with the loss of each transmission line of alternation current and DC power transmission line on entire passway for transmitting electricity
The sum of minimum target, constraint condition contain always the sending of each province's net and outer net/by electrical power Constraints of Equilibrium, main channel exchange connection
Winding thread transmission cross-section security constraint, the operation constraint of DC line, DC line rectification side power and inverter side power relation side
Journey exchanges the active reactive equilibrium equation of node, and considers that the dc power of actual motion adjusts restriction etc..
The objective function of AC/DC parallel passway for transmitting electricity power optimization distribution model is as follows:
In formula, T is to formulate the corresponding period sum of channel power plan of distribution, and Δ T is the time span of each period, with
1 hour is a period,Indicate active loss of the l articles transmission line of alternation current in the t period,Indicate that kth direct current is defeated
Active loss of the electric line in the t period, NacAnd NdcThe AC and DC that respectively entire AC/DC parallel passway for transmitting electricity includes is defeated
The total number of electric line.
For transmission line of alternation current, power loss can be expressed as follows:
In formula,WithThe respectively voltage magnitude of the head end and endpoint node of the l articles transmission line of alternation current period t,WithThe respectively voltage phase angle of the head end and endpoint node of the l articles transmission line of alternation current period t,It is the l articles exchange
The line conductance of transmission line of electricity.
For DC power transmission line, active loss can be expressed as follows:
N in formulakFor the operation number of poles of DC power transmission line k,WithRespectively list of the DC power transmission line k in period t
Pole running current and Unipolar resistance.
Constraint condition is as follows:
(1) always the sending of each province's net and outer net/by electrical power Constraints of Equilibrium
In formula,Transmission power for communication channel l in the t period,Transmission power for direct current channel k in the t period,
If the rectification side of DC link k is saving net s,For kth DC link period t rectification side
Power, if the inverter side of DC link k is saving net s,It is kth DC link in the inverse of period t
Become side power,To save total electricity exchange plan value of the net s in period t and outer net, ΩsaAnd ΩsdRespectively save net s and outer net
The alternating current interconnection set and DC link set of connection, ε % are the tolerance percentage of exchange of electric power plan between saving net.
The transmission power of communication channel calculates as follows:
In formula,WithThe respectively voltage magnitude of the head end and endpoint node of the l articles transmission line of alternation current period t,It is the susceptance of the l articles transmission line of alternation current;For the phase difference of voltage at the l articles transmission line of alternation current of period t both ends.
The transmission power of direct current channel calculates as follows:
In formula,For the monopolar operation voltage of period t direct current channel k.
(2) alternating current interconnection transmission cross-section security constraint
In formula,WithThe lower safety limit and upper safety limit of the safe section transimission power of respectively m-th of ac transmission.
(3) the operation constraint of DC line: it includes dc power equation (such as formula (6)) that the operation of DC line, which constrains, straight
Voltage equation, inverter fundamental equation and DC network equation are flowed, specific as follows:
In formula, Vi,tFor with DC node k be connected the period t voltage for exchanging node i,There is load to adjust for converter power transformer
It is pressed in the no-load voltage ratio of period t,For inverter period t pilot angle,For the equivalent reactance of converter power transformer;kγFor constant,
0.995 generally is taken,For direct current station period t power factor;WithRespectively period t DC power transmission line k rectification
The DC voltage of side and inverter side.
(4) power balance equation of node is exchanged:
In formula, PGi,tIt contributes for the generated power of node i, PLi,tFor the load active power of node i,For with exchange
The dc power of the connected DC node k of node i, converting plant are negative, and Inverter Station is positive;QGi,tFor the generator reactive of node i
Power output, QLi,tFor the reactive load power of node i;For the power-factor angle at direct current station;Vi,tThe electricity of node i is exchanged for period t
Pressure amplitude value, GijAnd BijTransconductance and mutual susceptance between node i and j,Phase difference of voltage between node i and j.
(5) the adjusting limitation of direct current transportation channel transmission power:
In formula,For kth article DC line the t period power regulation;WithRespectively kth item
DC line single regulation power allows the maximum value and minimum value changed;For 0/1 integer variable, kth AC line is indicated
The state whether road transmission power is adjusted in the t period indicates that the period maintains and a upper period power phase when equal to 0
Same numerical value, indicates that the period will increase or reduce on the basis of a upper period power when equal to 1.For kth direct current
Transmission line of electricity power in one day allows the number adjusted.
(6) the bound constraint of variable: the bound including generator node active power output and idle power output exchanges node
The bound of voltage magnitude, the bound of DC node dc power, the bound of DC node DC voltage and DC current,
The no-load voltage ratio of converter power transformer and the bound of pilot angle, it may be assumed that
xmin≤x≤xmax (11)
In formula,
It can be seen that due to the transmission power calculation formula (5) of communication channel, the first of inverter fundamental equation (8)
With second equation, the AC and DC line loss meter in the power balance equation formula (9) and objective function of node is exchanged
It calculates formula (2) and (3), the nonlinear operation containing trigonometric function and quadratic term.Therefore, above-mentioned formula (1)~(11) describe
The optimization distribution of AC/DC parallel passway for transmitting electricity power is mixed-integer nonlinear programming model, which is that planning field is most difficult to
One of the problem of solution, belongs to NP- difficulty problem, if carried out using the common mixed integer nonlinear programming solver such as SBB
It solves, calculating speed is very slow.
It therefore, in the method, will be by boring the processing methods such as relaxation converts and large M is equivalent, to nonlinear target
Function and constraint condition carry out convexification processing, convert the Optimized model of foundation to the MIXED INTEGER convex programming mould for being easy to solve
Type, specific convexification process are as follows:
1) convexification of the operation constraint of DC line
For most of scenes in practical power systems, busbar voltage amplitude is about 1p.u., therefore can be direct current
The voltage magnitude of node k connected exchange node i is similar to 1.0, then inverter fundamental equation is reduced toIt enablesI-th of pressure regulation shelves whether is operated in for k-th of converter station converter power transformer
Position, is one 0/1 binary system discrete variable,It is then the corresponding no-load voltage ratio value of i-th of pressure regulation gear,It is whole as one
The continuous variable of body, thenThen it is the product of a discrete variable and a continuous variable, line can be converted into using large M
Property constraint, enableThe equation can be as follows with large M conversion linear inequality:
Then inverter fundamental equation is converted into linear restriction:
For the transmission power equation of direct current channel, based in operation generally using the control for determining converting plant side DC voltage
Rectification side DC voltage is fixed as its voltage rating, i.e., by meansThe then transmission power of rectification sideIt is converted into linear restriction, the transmission power of inverter side It can be converted by second order cone relaxation are as follows:
It enablesThe reactive power of the absorption of DC node k has when disregarding the power loss of inverter:
By Vi,tSimilar to 1.0, enableThenTwo can be passed through
Rank cone relaxation conversion are as follows:
With large M pairIt linearizes:
(2) the second order cone relaxation of the active reactive equilibrium equation of node is exchanged
It enablesYij=Vi,tVj,tcosθij, Zij=Vi,tVj,tsinθij, then the power balance equation formula of node is exchanged
(9) it converts are as follows:
It can be converted by second order cone relaxation are as follows:
(3) linearisation of the power loss and transmission power of communication channel
The power loss formula (2) of communication channel can convert in objective function are as follows:
Save always sending/being converted by the transmission power formula (5) of communication channel in electrical power Constraints of Equilibrium for net and outer net are as follows:
To sum up, after carrying out convexification processing to the non-linear partial in Optimized model, the optimization point of AC and DC transmission channel power
With model conversation be objective function be Convex quadratic function, constraint condition is linear or second order cone MIXED INTEGER Second-order cone programming mould
Type is as follows:
Above-mentioned MIXED INTEGER Second-order cone programming model can be used mature optimization solver such as GUROBI and carry out high efficient and reliable
Solution, to obtain the power optimization allocation plan of direct current and ac transmission channel between each province;
Finally, the power allocation scheme for optimizing resulting AC/DC parallel passway for transmitting electricity to be handed down to the fortune of each converter station
Administrative staff executes.
In addition, in order to further verify the validity of this method, below with south electric network EHV transmission company O&M
The data instance on the 1st of August in 2017 derived from automated system carries out simulation calculation, and it is defeated to analyze AC/DC parallel proposed by the present invention
The application effect of the MIXED INTEGER Second-order cone programming method of electric channel power optimization distribution.
The interconnected network main grid structure structure chart such as Fig. 1, entire AC/DC parallel passway for transmitting electricity include that passway for transmitting electricity is flowed in 8 backcrossings
(" ox from " common-tower double-return DC power transmission line is merged into 1 time) with 8 times direct current transportation channels, wherein it is total that section is sent out in Yunnan
Have 5 DC power transmission lines (including western Shandong back-to-back DC), Guizhou sends out section and shares 2 DC power transmission lines and 4 friendships
Transmission line of electricity is flowed, Guangdong is sent into section and shares 7 DC power transmission lines and 8 transmission line of alternation currents.Entire main grid structure shares 147
A node, 183 branches.The fundamental characteristics parameter of the DC line between each province of table 1.
1 direct current channel of table (except western Shandong is back-to-back) fundamental characteristics parameter
The mixing distributed using the GUROBI solver in GAMS software AC/DC parallel passway for transmitting electricity power optimization is whole
Number Second-order cone programming model is solved, and it is as follows to obtain result: the alternating current circuit loss of whole day day part passway for transmitting electricity and direct current
Comparison such as Fig. 2 and Fig. 3 before line loss and optimization, it can be seen that the friendship of resulting day part after channel power optimization distribution
There is decline before all opposite optimization of the loss of Flow Line and DC line, wherein before the total losses of DC line are than optimization after optimization
25.5% is reduced, reduces by 12.8% when the total losses of alternating current circuit are before optimization.The network loss of passway for transmitting electricity before the optimization of this day
Electricity is 454.06MWh, and Energy loss is reduced to 350.35MWh after optimization, is that 0.3 yuan/degree calculates with network loss electricity price, the Ke Jie
About 31.113 ten thousand yuan of power purchase expense, economic benefit is obvious.It follows that using AC/DC parallel passway for transmitting electricity proposed by the present invention
The MIXED INTEGER Second-order cone programming method of power optimization distribution can effectively reduce via net loss, improve the economic effect of system operation
Benefit.The power transportation program in optimization front and back direct current transportation channel compares as shown in Fig. 4 a-4b, it can be seen that for defeated after optimization
The biggish DC line of electrical loss rate, as it is wide and Xingan's direct current, the transmission power after optimization are declined;And two extra-high voltages
Direct current such as Chu Sui and general emigrant and an other DC line ox are from transmission loss rate is smaller, and the transmission power after optimization has
Increased, the reduction that the two is combined with conducive to DC line loss.
Compare MIXED INTEGER Second-order cone programming method proposed by the present invention and conventional mixed integer nonlinear programming method
Optimal solution and calculate the time, as shown in table 2.Wherein, mixed-integer nonlinear programming model formula (1)~(11) are asked using SBB
Solution device is solved.It can be seen that three different periods optimum results, calculate time-consuming aspect, MIXED INTEGER Second-order cone programming side
The calculating time of method is all obviously reduced than mixed integer nonlinear programming;Also, the target of MIXED INTEGER Second-order cone programming method
Function is all smaller than mixed integer nonlinear programming, illustrates the Xie Gengyou of MIXED INTEGER Second-order cone programming method.With optimization mould
The increase of type calculation scale, such as optimization for one day 24 period, since the number of integer variable is excessive, MIXED INTEGER is non-linear
Planning is unable to get qualified integer solution, and MIXED INTEGER Second-order cone programming can still obtain within a short period of time it is optimal
Solution;If the integer variable value in the optimal solution that GUROBI solution MIXED INTEGER Second-order cone programming model is obtained is as known quantity
Substitute into mixed-integer nonlinear programming model in, the objective function for using SBB solver to solve again for 353.26MWh,
Target function value with MIXED INTEGER Second-order cone programming model solution is all very close to calculating time-consuming 992.8s, this also indicates that mixing
The calculating speed of integer Second-order cone programming model greatly improves, and obtained optimum results accuracy is high.
The Comparative result of 2 Different Optimization model difference solver of table
It follows that the MIXED INTEGER convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution provided by the present invention
Method has the following advantages compared with prior art:
1) compared to mixed integer nonlinear programming method, MIXED INTEGER method of convex programming can be obtained in mathematical meaning
Optimal solution, it is logical that obtained AC/DC parallel passway for transmitting electricity power optimization allocation plan can be effectively reduced AC and DC transmission of electricity
The loss in road improves the economic benefit of system operation;
2) compared to mixed integer nonlinear programming method, MIXED INTEGER method of convex programming can be substantially reduced Optimized model
Difficulty in computation, reduce the solution time of Optimized model, improve the computational efficiency of model solution, and can be improved Optimized model and ask
The reliability of solution.
That is, the present invention has having formulated the whole network generating set Active Generation operation plan and DC power transmission line
On the basis of function power transmission plan, by boring the processing methods such as relaxation converts and large M is equivalent, to the target with nonlinearity
Function and constraint condition carry out convexification processing, convert the Optimized model of foundation to the MIXED INTEGER convex programming mould for being easy to solve
Type, the AC and DC transmission channel power optimized allocation of acquisition can be effectively reduced the loss of alternating current circuit and DC line,
It with apparent economic benefit, and ensure that the optimality of gained solution mathematically, improve the solution efficiency that optimization calculates, with
Realize the performance driving economy of AC-DC interconnecting power network.
Simply to illustrate that technical concepts and features of the invention, its purpose is allows in the art above-described embodiment
Those of ordinary skill cans understand the content of the present invention and implement it accordingly, and it is not intended to limit the scope of the present invention.It is all
It is the equivalent changes or modifications that the essence of content according to the present invention is made, should be covered by the scope of protection of the present invention.
Claims (10)
1. the MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution, which is characterized in that including
Total exchange power planning between sending power data and each province to net under each exchange node is obtained, to establish alternating current-direct current simultaneously
Join passway for transmitting electricity power optimization distribution model, the AC/DC parallel passway for transmitting electricity power optimization distribution model is with entire passway for transmitting electricity
The minimum target of the sum of the loss of upper each transmission line of alternation current and DC power transmission line, and constraint condition is set;
By cone relaxation conversion and large M equivalent way, to AC/DC parallel passway for transmitting electricity power optimization distribution model
Non-linear partial in objective function and constraint condition carries out convexification processing, by the AC/DC parallel passway for transmitting electricity power optimization
Distribution model is converted into MIXED INTEGER Convex Programming Model;
MIXED INTEGER Convex Programming Model is solved, AC/DC parallel passway for transmitting electricity power optimization distribution side between each province is obtained
Case;
AC/DC parallel passway for transmitting electricity power optimization allocation plan between obtained each province is handed down to the operation of each converter station
Personnel execute.
2. the MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution as described in claim 1,
It is characterized in that, the objective function of AC/DC parallel passway for transmitting electricity power optimization distribution model are as follows:
In formula, T is to formulate the corresponding period sum of channel power plan of distribution, and Δ T is the time span of each period,Table
Show active loss of the l articles transmission line of alternation current in the t period,Indicate kth DC power transmission line in the active damage of t period
Consumption, NacAnd NdcThe total number for the AC and DC transmission line of electricity that respectively entire AC/DC parallel passway for transmitting electricity includes;
For transmission line of alternation current, power loss is expressed as follows:
In formula,WithThe respectively voltage magnitude of the head end and endpoint node of the l articles transmission line of alternation current period t,WithThe respectively voltage phase angle of the head end and endpoint node of the l articles transmission line of alternation current period t,It is the l articles ac transmission
The line conductance of route;
For DC power transmission line, active loss is expressed as follows:
N in formulakFor the operation number of poles of DC power transmission line k,WithThe respectively monopolar operation of DC power transmission line k period t
Electric current and Unipolar resistance.
3. the MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution as claimed in claim 2,
It is characterized in that, constraint condition includes always sending/being pacified by electrical power Constraints of Equilibrium, alternating current interconnection transmission cross-section for each province's net and outer net
Staff cultivation, the operation constraint of DC line, the adjusting of the power balance equation of exchange node, direct current transportation channel transmission power limit
System and the constraint of the bound of variable.
4. the MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution as claimed in claim 3,
It is characterized in that, always the sending of each province's net and outer net/by electrical power Constraints of Equilibrium are as follows:
In formula,Transmission power for communication channel l in the t period,Transmission power for direct current channel k in the t period, if directly
The rectification side for flowing interconnection k is saving net s, then For kth DC link period t rectification side power,
If the inverter side of DC link k is saving net s, For kth DC link period t inverter side function
Rate,To save total electricity exchange plan value of the net s in period t and outer net, ΩsaAnd ΩsdRespectively save what net s was connect with outer net
Alternating current interconnection set and DC link set, ε % are the tolerance percentage of exchange of electric power plan between saving net;
The transmission power of communication channel calculates as follows:
In formula,WithThe respectively voltage magnitude of the head end and endpoint node of the l articles transmission line of alternation current period t,It is
The susceptance of the l articles transmission line of alternation current;For the phase difference of voltage at the l articles transmission line of alternation current of period t both ends;
The transmission power of direct current channel calculates as follows:
In formula,For the monopolar operation voltage of period t direct current channel k.
5. the MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution as claimed in claim 3,
It is characterized in that, alternating current interconnection transmission cross-section security constraint are as follows:
In formula,WithThe lower safety limit and upper safety limit of the safe section transimission power of respectively m-th of ac transmission;
The operation of DC line constrains
DC voltage equation, inverter fundamental equation and DC network equation, specific as follows:
In formula, Vi,tFor with DC node k be connected the period t voltage for exchanging node i,Exist for converter power transformer on-load voltage regulation
The no-load voltage ratio of period t,For inverter period t pilot angle,For the equivalent reactance of converter power transformer;kγFor constant,For
Power factor of the direct current station in period t;WithThe respectively direct current of period t DC power transmission line k rectification side and inverter side
Voltage.
6. the MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution as claimed in claim 3,
It is characterized in that, exchanges the power balance equation of node are as follows:
In formula, PGi,tIt contributes for the generated power of node i, PLi,tFor the load active power of node i,For with exchange node
The dc power of i connected DC node k, converting plant are negative, and Inverter Station is positive;QGi,tIt contributes for the generator reactive of node i,
QLi,tFor the reactive load power of node i;For the power-factor angle at direct current station;Vi,tThe voltage amplitude of node i is exchanged for period t
Value, GijAnd BijTransconductance and mutual susceptance between node i and j,Phase difference of voltage between node i and j;
The adjusting of direct current transportation channel transmission power limits are as follows:
In formula,For kth article DC line the t period power regulation;WithRespectively kth direct current
Route single regulation power allows the maximum value and minimum value changed;For 0/1 integer variable, indicate that kth DC line is defeated
The state whether electrical power is adjusted in the t period indicates that the period remains identical with a upper period power when equal to 0
Numerical value indicates that the period will increase or reduce on the basis of a upper period power when equal to 1;For kth direct current transportation
Route power in one day allows the number adjusted.
7. the MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution as claimed in claim 3,
It is characterized in that, the bound constraint of variable includes: the bound of generator node active power output and idle power output, exchange node electricity
The bound of pressure amplitude value, the bound of DC node dc power, the bound of DC node DC voltage and DC current are changed
The no-load voltage ratio of convertor transformer and the bound of pilot angle, it may be assumed that
xmin≤x≤xmax(11)
In formula,
8. the MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution as claimed in claim 5,
It is characterized in that, convexification processing is carried out to the operation constraint of DC line:
Inverter fundamental equation is reduced toIt enables For k-th of converter station change of current change
Whether depressor operates in i-th of pressure regulation gear, is one 0/1 binary system discrete variable, Ti kIt is then corresponding for i-th of pressure regulation gear
No-load voltage ratio value,As the continuous variable of an entirety, thenIt is then a discrete variable and a continuous variable
Product, linear restriction is converted into using large M, enableThe equation is as follows with large M conversion linear inequality:
Then inverter fundamental equation is converted into linear restriction:
For the transmission power equation of direct current channel, rectification side DC voltage is fixed as its voltage rating, i.e.,Then
The transmission power of rectification sideIt is converted into linear restriction, the transmission power of inverter side
It is converted by second order cone relaxation are as follows:
It enablesThe reactive power of the absorption of DC node k has when disregarding the power loss of inverter:
By Vi,tSimilar to 1.0, enableThenPass through second order cone pine
Relax conversion are as follows:
With large M pairIt linearizes:
9. the MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution as claimed in claim 6,
It is characterized in that, enablesYij=Vi,tVj,tcosθij, Zij=Vi,tVj,tsinθij, then the power balance equation of node is exchanged
Formula (9) conversion are as follows:
It is converted by second order cone relaxation are as follows:
10. the MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution as claimed in claim 4,
It is characterized in that, the power loss formula (2) of transmission line of alternation current in objective function is converted are as follows:
Always the sending of net and outer net/by the transmission power formula (5) turn of the communication channel of communication channel in electrical power Constraints of Equilibrium will be saved
It turns to:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810804522.5A CN109038551B (en) | 2018-07-20 | 2018-07-20 | Hybrid integer convex programming method for power optimized distribution of alternating current-direct current parallel power transmission channel |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810804522.5A CN109038551B (en) | 2018-07-20 | 2018-07-20 | Hybrid integer convex programming method for power optimized distribution of alternating current-direct current parallel power transmission channel |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109038551A true CN109038551A (en) | 2018-12-18 |
CN109038551B CN109038551B (en) | 2020-12-08 |
Family
ID=64643878
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810804522.5A Active CN109038551B (en) | 2018-07-20 | 2018-07-20 | Hybrid integer convex programming method for power optimized distribution of alternating current-direct current parallel power transmission channel |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109038551B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109617113A (en) * | 2018-12-19 | 2019-04-12 | 南京南瑞继保电气有限公司 | The device and method mutually converted for dc power and AC power control target |
CN111598718A (en) * | 2020-04-03 | 2020-08-28 | 广东电网有限责任公司电力调度控制中心 | STATCOM configuration planning method, device and equipment |
CN112909959A (en) * | 2021-01-20 | 2021-06-04 | 广东电网有限责任公司 | Reactive power optimization method, device, equipment and storage medium |
CN114896794A (en) * | 2022-05-17 | 2022-08-12 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | AC/DC parallel power transmission channel maintenance plan and power distribution cooperative processing method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130304266A1 (en) * | 2012-04-13 | 2013-11-14 | Regents Of The University Of Minnesota | State estimation of electrical power networks using semidefinite relaxation |
CN104376378A (en) * | 2014-11-14 | 2015-02-25 | 浙江工商大学 | Distributed-power-source-contained power distribution network reactive power optimization method based on mixed integer cone optimization |
CN104915724A (en) * | 2015-04-22 | 2015-09-16 | 华南理工大学 | AC-DC parallel power transmission channel power optimization distribution method and system |
CN105119280A (en) * | 2015-08-31 | 2015-12-02 | 天津大学 | Conic optimization-based AC/DC hybrid structure active power distribution network operation optimization method |
-
2018
- 2018-07-20 CN CN201810804522.5A patent/CN109038551B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130304266A1 (en) * | 2012-04-13 | 2013-11-14 | Regents Of The University Of Minnesota | State estimation of electrical power networks using semidefinite relaxation |
CN104376378A (en) * | 2014-11-14 | 2015-02-25 | 浙江工商大学 | Distributed-power-source-contained power distribution network reactive power optimization method based on mixed integer cone optimization |
CN104915724A (en) * | 2015-04-22 | 2015-09-16 | 华南理工大学 | AC-DC parallel power transmission channel power optimization distribution method and system |
CN105119280A (en) * | 2015-08-31 | 2015-12-02 | 天津大学 | Conic optimization-based AC/DC hybrid structure active power distribution network operation optimization method |
Non-Patent Citations (4)
Title |
---|
SHAHAB BAHRAMI等: "Semidefinite Relaxation of Optimal Power Flow for AC–DC Grids", 《IEEE TRANSACTIONS ON POWER SYSTEMS 》 * |
余秀月等: "含电压源换流器的交直流混联电网无功优化模型", 《电力系统保护与控制》 * |
周举: "特高压交直流系统的可用输电能力计算", 《红水河》 * |
张艺镨等: "基于极限场景的两阶段含分布式电源的配网无功优化", 《电工技术学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109617113A (en) * | 2018-12-19 | 2019-04-12 | 南京南瑞继保电气有限公司 | The device and method mutually converted for dc power and AC power control target |
CN109617113B (en) * | 2018-12-19 | 2022-02-18 | 南京南瑞继保电气有限公司 | Device and method for mutual conversion between direct current power and alternating current power control targets |
CN111598718A (en) * | 2020-04-03 | 2020-08-28 | 广东电网有限责任公司电力调度控制中心 | STATCOM configuration planning method, device and equipment |
CN112909959A (en) * | 2021-01-20 | 2021-06-04 | 广东电网有限责任公司 | Reactive power optimization method, device, equipment and storage medium |
CN114896794A (en) * | 2022-05-17 | 2022-08-12 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | AC/DC parallel power transmission channel maintenance plan and power distribution cooperative processing method |
CN114896794B (en) * | 2022-05-17 | 2024-05-28 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | AC/DC parallel power transmission channel maintenance plan and power distribution cooperative processing method |
Also Published As
Publication number | Publication date |
---|---|
CN109038551B (en) | 2020-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2021135332A1 (en) | Edge computing-based multi-agent load regulation and control method | |
CN109038551A (en) | The MIXED INTEGER method of convex programming of AC/DC parallel passway for transmitting electricity power optimization distribution | |
CN104915724B (en) | AC/DC parallel passway for transmitting electricity optimized power allocation method and system | |
CN108599157B (en) | Alternating current-direct current cooperation power grid optimal scheduling method considering electricity price type demand response | |
CN110263391B (en) | Opportunity constraint-based intelligent soft switch planning method for active power distribution network | |
CN106936152B (en) | Consider the ac and dc systems voltage and reactive power coordinated control method of converter station loss characteristic | |
CN109950907A (en) | The dispatching method and system of alternating current-direct current mixing power distribution network containing electric power electric transformer | |
CN109617092A (en) | A kind of alternating current-direct current mixing power grid dynamic reactive power optimization method and system | |
CN107749635A (en) | A kind of closed loop network running optimizatin method | |
CN108023364A (en) | Power distribution network distributed generation resource maximum access capability computational methods based on convex difference planning | |
CN105610158A (en) | Distributed power flow controller and control method therefor | |
CN104269847A (en) | Flexible looped network control system operation and power flow optimization method | |
CN104362650B (en) | A kind of Method for Reactive Power Optimization in Power for considering cost factor | |
CN110350600A (en) | A kind of flexible multimode switch regulation method promoting distributed generation resource consumption | |
CN106655234A (en) | Method for analyzing influence of line impedance and call wire power on generalized short-circuit ratio | |
CN107528321B (en) | Random optimal distribution method for power of interconnected power grid alternating current-direct current parallel power transmission channel | |
CN109377020B (en) | Power transmission network planning method considering load transfer capacity of power distribution network | |
CN105226649B (en) | One kind predicting improved provincial power network power generation dispatching optimization method based on bus load | |
CN110299713A (en) | A kind of tractive power supply system imbalance of three-phase voltage compensation method counted and wind power plant influences | |
CN116154849A (en) | Power grid distributed power supply bearing capacity assessment method | |
CN110198053A (en) | It is a kind of to concentrate with the micro-capacitance sensor real-time voltage control method and system combined on the spot | |
CN104201697A (en) | Optimal reactive power compensation method for 110KV and 35KV power grids | |
CN117371154B (en) | Power grid planning system of industrial frequency division interconnection system with M3C frequency conversion station for site selection and volume fixation | |
Ziyat et al. | Real-time Reactive Power Compensation by Distributed Generation Simulated using GridLAB-D and PSIM | |
CN209056932U (en) | A kind of integrated energy system topological structure |
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 |