CN107359616A - One kind solves the not convergent method of Load flow calculation after large scale electric network mode adjusts - Google Patents
One kind solves the not convergent method of Load flow calculation after large scale electric network mode adjusts Download PDFInfo
- Publication number
- CN107359616A CN107359616A CN201710583770.7A CN201710583770A CN107359616A CN 107359616 A CN107359616 A CN 107359616A CN 201710583770 A CN201710583770 A CN 201710583770A CN 107359616 A CN107359616 A CN 107359616A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- subnet
- unbalance
- newton
- 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
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- 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]
Abstract
The invention discloses one kind to solve the not convergent method of Load flow calculation after large scale electric network mode adjusts, and belongs to technical field of electric power system control.Present invention meter and mode adjust the emergent control measure after operation, are shared using imbalance power, and it is inferior that meter and upper wheel do not restrain newton pressgang in Load flow calculation（NR）The optimal solution of iteration simultaneously adjusts lower whorl calculating iterative initial value on this basis, and the trend mode self-optimizing model of the interior point method optimal load flow according to ill feature is proposed to not restraining trend mode, obtains the convergence trend after large scale electric network mode adjusts.The present invention is efficiently solved after mode adjusts operation not to be restrained routinely by PQ decomposition with the flow calculation program based on Newton-Laphson method, mode can not directly be calculated and adjust the trend after operation, quickly meter can not be provided for dispatcher and mode adjusts the safety and stability evaluation of operation and the problem of early warning result.
Description
Technical field
The invention belongs to technical field of electric power system control, and more precisely, the present invention relates to one kind to solve extensive electricity
The not convergent method of Load flow calculation after net mode adjustment.
Background technology
The offline research submodule of the online dynamic secure estimation of power system and early warning system, it usually needs analog form is adjusted
The whole influence to power system security stability characteristic (quality).Mode, which adjusts the operation being related to, to be included：Generator power adjustment, load power
Setting, load tap changer are set, circuit/bus/plant stand is stopped transport etc..With extra-high voltage alternating current-direct current and Large Copacity transformer station
Build, after mode adjusts operation, it is understood that there may be the imbalance power of Large Copacity, it may occur that large-scale power flow transfer, for example (,) it is special
HVDC stoppage in transit etc., its result can cause the conventional flow calculation program based on PQ decomposition and Newton-Raphson NR methods
Do not restrain, the mode that can not directly calculate adjusts the trend after operation, quickly can not provide meter for dispatcher and mode adjusts behaviour
The safety and stability evaluation and early warning result of work.
Load flow calculation, which is not restrained, mainly two aspect reasons：(1) power flow equation has solution, but due to tidal current computing method not
Improve and can not seek solution, be convergence of algorithm sex chromosome mosaicism.Flow solution can be tried to achieve by improving the convergence of tidal current computing method.
(2) power flow equation real solution is not present.Equally show as calculating not restraining when Load flow calculation.Now need to pass through adjustment
Systematic parameter such as generator output, load tap changer, cut-out load etc. change power flow equation, make the trend side after change
Journey has solution.
At present, two class technique studies have a lot.Have in terms of Research on Calculation：The improvement of Newton-like method；Newton
The Optimal Multiplier method that method is combined with minimizing trend；Using the Nonlinear Programming Method of optimization thought；There to be global convergence to calculate
Method such as Homotopy Method is applied in electric power system tide calculating.Research in terms of flow data adjustment also has a lot, is mostly
From the relation analyzing between the solvability of power flow equation and voltage stability, by adjust the reactive power of key node come
Obtain trend convergence solution.
Obviously, to solve power system operation mode adjustment after Load flow calculation not convergence problem, it is necessary to carry out Load flow calculation
Adjustment in terms of method and calculating data.How reasonably to reflect the method for operation adjustment after system dynamic changing process and provide
Mode calculates convergent steady-state load flow after adjusting, lack paractical research at present.
The content of the invention
The present invention seeks to：In view of the shortcomings of the prior art, propose that the mode of a set of Practical adjusts injustice after operation
Weighing apparatus power share simulation, flow data inspection, trend initial value give, trend convergence automatic identification, trend mode adjust automatically
Method, for solve the problems, such as power system operation mode adjustment after Load flow calculation do not restrain.
Specifically, the present invention is realized using following technical scheme, is comprised the following steps：
1) remember that power system initial launch mode is S0, to S0Ways of carrying out adjustment operation, which adjustment operation is designated as
F, the equipment composition of each subnet and node injection after F operations are determined, calculates active amount of unbalance caused by the operation of each subnet, note
Active amount of unbalance caused by j-th of subnet operation is Δ Pj,c, retain each subnet that can be connected by direct current, into step 2)；
2) whether judgment mode adjustment operation F can start emergent control measure, if emergent control measure will not be started,
The active amount of unbalance for determining each subnet is active amount of unbalance caused by the operation of each subnet, that is, has Δ Pj=Δ Pj,c, wherein Δ
PjFor the active amount of unbalance of j-th of subnet；
Otherwise, obtain respectively each subnet due to when caused active amount of unbalance is implemented in value measure, remember j-th of subnet by
It is Δ P in implementing caused active amount of unbalance when value measurej,a, each subnet when value measure caused by implementing
Active amount of unbalance is when machine amount is cut in selection to be negative when selecting cutting load amount just；Determine the active amount of unbalance of each subnet
Add up for active amount of unbalance caused by the operation of each subnet and due to the active amount of unbalance caused by value measure implementation, i.e.,
There is Δ Pj=Δ Pj,c+ΔPj,a；
If 3) the active amount of unbalance of each subnet is equal to 0, the unallocated amount of unbalance for counting each subnet is 0, remembers j-th of subnet
Unallocated amount of unbalance be Δ Pj,t, into step 4)；
If the active amount of unbalance of each subnet is not equal to 0, to the active amount of unbalance of each subnet, according to what is be previously set
Imbalance power simulation, which shares rule and carries out simulation, shares, the power balancing procedure after simulated operation, obtains each after simulation is shared
The unallocated amount of unbalance of subnet, into step 4)；
4) method of operation that the active amount of unbalance simulation of each subnet is shared after adjustment is designated as S1, to S1Carry out conventional tide
The inspection of the small branch impedance data of stream calculation, it is 1 that juxtaposition Newton-Raphson NR iteration round, which counts NT, when small branch impedance
Less than setting minimum impedance value when, then limit its as set minimum impedance value；
If 5) it is 1 that Newton-Raphson NR iteration round, which counts NT, using the fixation PQ decomposition method numbers of setting, opened by flat
The dynamic trend iterative initial value V for obtaining Newton-Raphson NR iteration(0), as Newton-Raphson NR iteration in epicycle Load flow calculation
Trend iterative initial value, the flat startup refers to take all PQ node voltage amplitudes to be equal to 1.0, and PV node voltage sets for given
Determine magnitude of voltage, phase angle is 0；
If Newton-Raphson NR iteration round, which counts NT, is more than 1, remember in previous Newton-Raphson NR iterative process most
Excellent solution is designated as V(k), judge corresponding optimal solution V(k)Newton-Raphson NR iterative powers deviation whether meet setting threshold value：
If satisfied, then using adjustment trend initial value strategy restarting Load flow calculation, the trend iterative initial value of Newton-Raphson NR iteration
The initial value V obtained is solved using the flat startup after change PQ decomposition method numbers(0)With in previous Newton-Raphson NR iterative process
Optimal solution V(k)Weighting, i.e. V=α V(0)+βV(k), wherein α and β are the weight coefficient set and have alpha+beta=1, using V as epicycle
The trend iterative initial value of Newton-Raphson NR iteration in Load flow calculation；If not satisfied, then directly using change PQ decomposition method number plans
Slightly, and as Newton-Raphson NR iteration wheel numbers NT increases automatic increase PQ decomposition method numbers, initial value V is obtained by flat start(0),
Trend iterative initial value as Newton-Raphson NR iteration in epicycle Load flow calculation；
After the trend iterative initial value for determining Newton-Raphson NR iteration, the iteration of NT wheel Newton-Raphson NR iteration is put
Number K=1, into step 6)；
6) candidate's balancing machine is simulated in Newton-Raphson NR iteration, adjustable space or capacity are pressed between candidate's balancing machine
The unallocated amount of unbalance or balancing machine that size shares each subnet are more limited the quantity, and voltage magnitude when calculating iteration and phase angle change
For attenuation ratio index；
7) if Newton-Raphson NR iteration trend restrains, into step 8)；Otherwise, if epicycle Newton-Raphson NR
The iterations K of iteration then makes K increase 1, continues epicycle newton-pressgang into step 6) no more than the maximum times Kmax of setting
Inferior NR iteration, if the iterations K of epicycle Newton-Raphson NR iteration is equal to the maximum iteration Kmax of setting, judge
Whether epicycle wheel number NT is more than the parameter NTmax of setting, made if being not more than NT increase 1 enter step 5) progress next round newton-
The inferior NR iteration of pressgang, such as larger than then according to the voltage magnitude and phase angle in Newton-Raphson NR iterative process at index catastrophe point
Iterative attenuation determines Load Flow Solution morbid state feature than index size, enters step 8) after carrying out automatic load flow rectification；
8) mode Load flow calculation terminates.
Above-mentioned technical proposal is further characterized by, and the mode of the step 1) adjusts operation F, including generator power
Operation including adjustment, load power are set, load tap changer is set and circuit/bus/plant stand is stopped transport.
Above-mentioned technical proposal is further characterized by, to the active amount of unbalance of each subnet in the step 3), according to
The imbalance power simulation being previously set, which shares rule and carries out simulation, shares, and the power balancing procedure after simulated operation, obtains mould
Intend the unallocated amount of unbalance for sharing rear each subnet, including procedure below：
If selection is shared by simulation primary frequency modulation response, into step 3-1), if selection is by simulation frequency modulation frequency modulation response
Share, then into step 3-2), if neither selecting to share by simulation primary frequency modulation response nor selecting to ring by simulation frequency modulation frequency modulation
It should share, then the unallocated amount of unbalance for counting each subnet is the active amount of unbalance of each subnet；
3-1) become according to the difference coefficient of each unit, the current loads of each load bus and static load frequency in each subnet
Change 1% and cause active change percentage, each unit and each load bus are shared by formula (1) and (2) respectively：
Wherein：Ri,jFor the difference coefficient of i-th of unit in j-th of subnet, Plk0,jFor k-th of load section in j-th of subnet
The current loads of point, lk,jStatic load frequency change 1% for k-th of load bus in j-th of subnet causes active change hundred
Divide ratio, Δ PGi,jWith Δ Plk,jThe imbalance that the simulation of i-th of unit and k-th of load bus is shared respectively in j-th of subnet
Amount；
If after sharing by formula (1), some unit output is out-of-limit, then fixes the unit output, represents the unit primary frequency modulation
Ability is failed, and remainder is continued to share inside the unit and load bus still with Primary frequency control ability；
After sharing end, the unallocated amount of unbalance for putting each subnet is 0；
3-2) the active amount of unbalance according to each subnet, Automation generation control unit allocation logic is simulated, according to first area
AGC shares in the control deviation ACE areas of domain, and rear upper level scheduling administration AGC participates in sharing, finally AGC machines in other ACE areas in net
Group participates in the mode shared, and shares as follows：
By the active amount of unbalance of each subnet according to ACE regional extents, its affiliated ACE region is defined；
If for the active amount of unbalance of some subnet, power shortage such as formula (3) caused by some ACE region belonging to it
It is described, represent that the power shortage can be maintained to share balance by the AGC units in the ACE regions, then by the power shortage at this
Share between AGC units in ACE regions, now frequency difference is 0：
Or
Wherein, Δ PiPower shortage caused by i-th of ACE region belonging to the active amount of unbalance of some subnet is represented,
Pj0、WithRepresent respectively in i-th of ACE region under current output and the output upper limit and output of j-th of AGC unit
Limit；
If for the active amount of unbalance of some subnet, power shortage such as formula (4) caused by some ACE region belonging to it
It is described, represent that the power shortage exceedes the adjustable extent of AGC units in the ACE regions, but still AGC machines are administered in upper level degree
The adjustable extent of group, then share the upper limit or lower limit, overage is dispatched by upper level by AGC units in the ACE regions
Administration AGC participates in sharing, and it is 0 finally to realize frequency difference：
Or
Wherein, S represents the set of upper level degree administration AGC units, Ps0、WithThe administration of upper level degree is represented respectively
Current output and the output upper limit and output lower limit of s-th of AGC unit in AGC units；
If for the active amount of unbalance of some subnet, power shortage such as formula (5) caused by some ACE region belonging to it
It is described, represent that the power shortage exceedes the adjustable extent of AGC units and the AGC units of schedule level one administration thereon in this ACE areas,
AGC units in the ACE regions and the AGC units of schedule level one administration thereon are then shared into the upper limit or lower limit, overage by
AGC units participate in droop control together in other ACE regions in net：
Or
After sharing end, the unallocated amount of unbalance for putting each subnet is 0.
Above-mentioned technical proposal is further characterized by, and the step 6), specifically includes following steps：
If 6-1) quantity of candidate's balancing machine is 0, into step 6-3),
If the quantity of candidate's balancing machine is not 0, the unallocated amount of unbalance of such as each subnet is not equal to 0 and epicycle newton-drawing
The iterations K of the inferior NR iteration of husband is equal to 1, then the unallocated amount of unbalance of each subnet is pressed each subnet candidate balancing machine can
Space or amount of capacity is adjusted to share out equally, into step 6-2)；
If 6-2) meet that the accumulative active deviation of each node is less than setting value Psh and school in Newton-Raphson NR iteration
The balancing machine of core, which is contributed, has out-of-limit deviation, then shares out equally by adjustable space or amount of capacity and more limit the quantity, into step 6-3)；
Voltage magnitude and phase angle iterative attenuation when 6-3) calculating Newton-Raphson NR iteration by the following method compare index：
If the correction [Δ x] during epicycle Newton-Raphson NR iteration kths(k)Middle voltage magnitude and phase angle greatest member
Absolute value be respectively Δ U(k)With Δ δ(k), by formulaWithVoltage magnitude and phase angle iterative attenuation ratio are calculated respectively
Index, wherein Δ U(k-1)With Δ δ(k-1)Respectively epicycle Newton-Raphson NR iteration the K-1 times when correction [Δ x](k-1)
The absolute value of middle voltage magnitude and phase angle greatest member.
Above-mentioned technical proposal is further characterized by, according in Newton-Raphson NR iterative process in the step 7)
Voltage magnitude and phase angle iterative attenuation at indexes suddenly changed point determine Load Flow Solution morbid state feature than index size, carry out automatic tide
Stream adjustment, specifically includes procedure below：
If each indexes suddenly changed point Δ U(k)/ΔU(k-1)The size of sum is more than Δ δ(k)/Δδ(k-1)The size of sum, then it is damp
Stream solves ill feature and shows as local reactive power support deficiency, into step 7-1)；Otherwise Load Flow Solution morbid state feature is shown as
Partial cut-away trend is overweight, into step 7-2)：
The trend mode of local reactive power support deficiency 7-1) is shown as Load Flow Solution morbid state feature, is generated electricity by adjusting
Machine is idle to contribute, throws/cut electric capacity, reactor group, on the premise of the various constraints of system are met, makes to give unit in total system
Active power output adjustment is minimum, and its mathematical modeling is described such as formula (6)：
Wherein：SGRepresent the adjustable unit output set of total system, PG0iAnd PGiThe respectively adjustable unit output set of total system
In the front and rear active power of i-th unit adjustment, wPGFor the active weight coefficient of corresponding unit；H (x) is decision variable x functions
Equality constraint, refer to the network and power-balance constraint of power system；Inequality constraints g (x) is by bound gDetermine
Operation constraint, the bound gFor the idle bound contributed, throw/cut electric capacity and reactor group of controllable electric generator；
The overweight trend mode of section tidal current 7-2) is shown as Load Flow Solution morbid state feature, is had by adjusting generator
Work(and idle output, on the premise of the various constraints of system are met, make the active power output adjustment for giving unit in total system minimum,
Mathematical modeling describes same formula (6), and now bound g, g is the bound of controllable electric generator active power output and idle output.
Beneficial effects of the present invention are as follows：The present invention is directed to the Load flow calculation after large scale electric network mode adjusts, it is proposed that
After the operation of a set of Practical imbalance power share simulation, flow data inspection, trend initial value be given, trend convergence from
Dynamic identification, the method for trend mode adjust automatically.The present invention is efficiently solved after mode adjusts operation and decomposed and newton-drawing with PQ
Conventional Load Flow calculation procedure based on Fu Xunfa is not restrained, and the mode that can not directly calculate adjusts the trend after operation, can not
Quickly meter is provided for dispatcher and mode adjusts the safety and stability evaluation of operation and the problem of early warning result.
Brief description of the drawings
Fig. 1 is the flow chart of the inventive method.
Embodiment
The present invention is described in further detail with reference to the accompanying drawings.
Step 1 describes in Fig. 1, and note power system initial launch mode is S0, to S0Ways of carrying out adjustment operation (hair
Power of motor adjustment, load power are set, load tap changer is set and circuit/bus/plant stand is stopped transport etc.), which is adjusted
Whole operation is designated as F, determines the equipment composition of each subnet and node injection after F operations, calculates active caused by the operation of each subnet
Amount of unbalance, remember that active amount of unbalance caused by j-th of subnet operation is Δ Pj,c, retain each subnet that can be connected by direct current.
Step 2 describes in Fig. 1, and whether judgment mode adjustment operation F can start emergent control measure, if will not open
Dynamic emergent control measure, it is determined that the active amount of unbalance of each subnet is active amount of unbalance caused by each subnet operates, that is, is had
ΔPj=Δ Pj,c, wherein Δ PjFor the active amount of unbalance of j-th of subnet；
Otherwise, obtain respectively each subnet due to when caused active amount of unbalance is implemented in value measure, remember j-th of subnet by
It is Δ P in implementing caused active amount of unbalance when value measurej,a, each subnet is due to active caused by value measure implementation
Amount of unbalance is when machine amount is cut in selection to be negative when selecting cutting load amount just；The active injustice of operation of each subnet is considered simultaneously
Weigh and start emergent control and work as the active amount of unbalance that value measure is implemented, add up each subnet because operation causes and starts urgent control
System cuts machine-cut load when value measure, that is, the active amount of unbalance for determining each subnet is active injustice caused by the operation of each subnet
Weigh and add up due to the active amount of unbalance caused by value measure implementation, that is, have Δ Pj=Δ Pj,c+ΔPj,a。
Step 3 describes in Fig. 1, and active amount of unbalance is shared in simulation.If the specifically, active amount of unbalance of each subnet
Equal to 0, the unallocated amount of unbalance for counting each subnet is 0, and the unallocated amount of unbalance of j-th of subnet of note is Δ Pj,t, into step
4；If the active amount of unbalance of each subnet is not equal to 0, to the active amount of unbalance of each subnet, according to the imbalance being previously set
Power analog, which shares rule and carries out simulation, to be shared, the power balancing procedure after simulated operation, is obtained simulation and is shared rear each subnet
Unallocated amount of unbalance, into step 4.
If selection is shared by simulation primary frequency modulation response, into step 3-1, if selection is by simulation frequency modulation frequency modulation response point
Stand, then into step 3-2, if neither selecting to share by simulation primary frequency modulation response nor selecting by simulation frequency modulation frequency modulation response point
Stand, then the unallocated amount of unbalance for counting each subnet is the active amount of unbalance of each subnet.
Step 3-1 is that simulation primary frequency modulation response is shared.According to the difference coefficient of each unit, each load bus in each subnet
Current loads and the change of static load frequency 1% cause active change percentage, each unit and each load bus are pressed respectively
Formula (1) and (2) are shared：
Wherein：Ri,jFor the difference coefficient of i-th of unit in j-th of subnet, Plk0,jFor k-th of load section in j-th of subnet
The current loads of point, lk,jStatic load frequency change 1% for k-th of load bus in j-th of subnet causes active change hundred
Divide ratio, Δ PGi,jWith Δ Plk,jThe imbalance that the simulation of i-th of unit and k-th of load bus is shared respectively in j-th of subnet
Amount；
If after sharing by formula (1), some unit output is out-of-limit, then fixes the unit output, represents the unit primary frequency modulation
Ability is failed, and remainder is continued to share inside the unit and load bus still with Primary frequency control ability；
After sharing end, the unallocated amount of unbalance for putting each subnet is 0；
Step 3-2 is that simulation frequency modulation frequency modulation response is shared.According to the active amount of unbalance of each subnet, automatic generation control is simulated
AGC unit allocations logic processed, share according to AGC in first district control deviation ACE areas, rear upper level scheduling administration AGC participates in dividing
Stand, finally AGC units participate in the mode shared in other ACE areas in net, share as follows：
By the active amount of unbalance of each subnet according to ACE regional extents, its affiliated ACE region is defined (because subnet j is not
Aequum Δ PjCause for ACE control areas i operations and emergent control when value acts, Δ PjMultiple ACE control zones may be corresponded to
Domain, now each regional imbalance amount is simulated respectively)；
If for the active amount of unbalance of some subnet, power shortage such as formula (3) caused by some ACE region belonging to it
It is described, represent that the power shortage can be maintained to share balance by the AGC units in the ACE regions, then by the power shortage at this
Share between AGC units in ACE regions, now frequency difference is 0：
Or
Wherein, Δ PiPower shortage caused by i-th of ACE region belonging to the active amount of unbalance of some subnet is represented,
Pj0、WithRepresent respectively in i-th of ACE region under current output and the output upper limit and output of j-th of AGC unit
Limit；
If for the active amount of unbalance of some subnet, power shortage such as formula (4) caused by some ACE region belonging to it
It is described, represent that the power shortage exceedes the adjustable extent of AGC units in the ACE regions, but still AGC machines are administered in upper level degree
The adjustable extent of group, then share the upper limit or lower limit, overage is dispatched by upper level by AGC units in the ACE regions
Administration AGC participates in sharing, and it is 0 finally to realize frequency difference：
Or
Wherein, S represents the set of upper level degree administration AGC units, Ps0、WithThe administration of upper level degree is represented respectively
Current output and the output upper limit and output lower limit of s-th of AGC unit in AGC units；
If for the active amount of unbalance of some subnet, power shortage such as formula (5) caused by some ACE region belonging to it
It is described, represent that the power shortage exceedes the adjustable extent of AGC units and the AGC units of schedule level one administration thereon in this ACE areas,
AGC units in the ACE regions and the AGC units of schedule level one administration thereon are then shared into the upper limit or lower limit, overage by
AGC units participate in droop control together in other ACE regions in net：
Or
After sharing end, the unallocated amount of unbalance for putting each subnet is 0.
Step 4 describes in Fig. 1, the method for operation S that the amount of unbalance simulation to each subnet is shared after adjustment1Carry out normal
The inspection of the small branch impedance data of Load flow calculation is advised, it is 1 that juxtaposition Newton-Raphson NR iteration round, which counts NT,.When small branch road
Impedance is too small, less than the minimum impedance value (taking 1e-5 when such as single precision calculates) of setting, then limits its minimum impedance as setting
Value.
Step 5 is described in Fig. 1, and strategy is given to the trend initial value of Newton-Raphson NR iteration.If newton-pressgang
It is 1 that inferior NR iteration round, which counts NT, using the fixation PQ decomposition method numbers of setting, is changed by the flat acquisition Newton-Raphson NR that starts
The trend iterative initial value V in generation(0), it is described flat as the trend iterative initial value of Newton-Raphson NR iteration in epicycle Load flow calculation
Start and refer to take all PQ node voltage amplitudes to be equal to 1.0, PV node voltage is given setting magnitude of voltage, and phase angle is 0.
If Newton-Raphson NR iteration round, which counts NT, is more than 1, remember in previous Newton-Raphson NR iterative process most
Excellent solution is designated as V(k), judge corresponding optimal solution V(k)Newton-Raphson NR iterative powers deviation whether meet setting threshold value
(such as integral power deviation<Node total number × convergence precision)：If satisfied, then using adjustment trend initial value strategy restarting trend
Calculate, the trend iterative initial value of Newton-Raphson NR iteration solves what is obtained using the flat startup after change PQ decomposition method numbers
Initial value V(0)With the optimal solution V in previous Newton-Raphson NR iterative process(k)Weighting, i.e. V=α V(0)+βV(k), wherein α and β are
The weight coefficient of setting simultaneously has alpha+beta=1, using V as at the beginning of the trend iteration of Newton-Raphson NR iteration in epicycle Load flow calculation
Value；If not satisfied, then directly using change PQ decomposition method number strategies, and with Newton-Raphson NR iteration wheel numbers NT increases certainly
Dynamic increase PQ decomposition method numbers, initial value V is obtained by flat start(0), as Newton-Raphson NR iteration in epicycle Load flow calculation
Trend iterative initial value.
After the trend iterative initial value for determining Newton-Raphson NR iteration, the iteration of NT wheel Newton-Raphson NR iteration is put
Number K=1.
Step 6 describes in Fig. 1, candidate's balancing machine is simulated in Newton-Raphson NR iteration, between candidate's balancing machine
The unallocated amount of unbalance or balancing machine for sharing each subnet by adjustable space or amount of capacity are more limited the quantity, and when calculating iteration
Voltage magnitude and phase angle iterative attenuation compare index.It is divided into three parts：
Step 6-1：If the quantity of candidate's balancing machine is 0, into step 6-3；
If the quantity of candidate's balancing machine is not 0, the unallocated amount of unbalance of such as each subnet is not equal to 0 and epicycle newton-drawing
The iterations K of the inferior NR iteration of husband is equal to 1 (when i.e. epicycle iteration starts), then to the unallocated amount of unbalance of each subnet in each son
Shared out equally between net candidate's balancing machine by adjustable space or amount of capacity, into step 6-2；
Step 6-2：If meeting in Newton-Raphson NR iteration, the accumulative active deviation of each node is less than setting value Psh
(the system reference capacity for such as taking 0.2 times) and the balancing machine checked are contributed and out-of-limit deviation be present, then big by adjustable space or capacity
Small share out equally more is limited the quantity, into step 6-3；
Step 6-3：Voltage magnitude and phase angle iterative attenuation ratio when calculating Newton-Raphson NR iteration by the following method refer to
Mark：
If the correction [Δ x] during epicycle Newton-Raphson NR iteration kths(k)Middle voltage magnitude and phase angle greatest member
Absolute value be respectively Δ U(k)With Δ δ(k), by formulaWithVoltage magnitude and phase angle iterative attenuation ratio are calculated respectively
Index, wherein Δ U(k-1)With Δ δ(k-1)Respectively epicycle Newton-Raphson NR iteration the K-1 times when correction [Δ x](k-1)
The absolute value of middle voltage magnitude and phase angle greatest member.
Step 7 describes in Fig. 1, if Newton-Raphson NR iteration trend restrains, into step 8；Otherwise, if
The iterations K of epicycle Newton-Raphson NR iteration no more than setting maximum times Kmax, then make K increase 1, into step 6 after
Continuous epicycle Newton-Raphson NR iteration, if the iterations K of epicycle Newton-Raphson NR iteration is equal to the greatest iteration of setting
Number Kmax, then judge whether epicycle wheel number NT is more than the parameter NTmax (as set NTmax=3) of setting, is made as no more than if
NT increases 1 and enters step 5 progress next round Newton-Raphson NR iteration, such as larger than then according in Newton-Raphson NR iterative process
Voltage magnitude and phase angle iterative attenuation at indexes suddenly changed point determine Load Flow Solution morbid state feature than index size, carry out automatic tide
Enter step 8 after stream adjustment.
If each indexes suddenly changed point Δ U(k)/ΔU(k-1)The size of sum is more than Δ δ(k)/Δδ(k-1)The size of sum, then it is damp
Stream solves ill feature and shows as local reactive power support deficiency, into step 7-1；Otherwise Load Flow Solution morbid state feature shows as portion
Divide section tidal current overweight, into step 7-2：
Step 7-1：The trend mode of local reactive power support deficiency is shown as Load Flow Solution morbid state feature, passes through adjustment
Generator reactive contributes, throws/cut the various control variables such as electric capacity, reactor group, on the premise of the various constraints of system are met, makes
The active power output adjustment minimum of unit is given in total system, its mathematical modeling is described such as formula (6)：
Wherein：SGRepresent the adjustable unit output set of total system, PG0iAnd PGiThe respectively adjustable unit output set of total system
In the front and rear active power of i-th unit adjustment, wPGFor the active weight coefficient of corresponding unit；H (x) is decision variable x functions
Equality constraint, refer to the network and power-balance constraint of power system；Inequality constraints g (x) is by bound gDetermine
Operation constraint, the bound gFor the idle bound contributed, throw/cut electric capacity and reactor group of controllable electric generator；
Step 7-2：The overweight trend mode of section tidal current is shown as Load Flow Solution morbid state feature, is generated electricity by adjusting
The control variables such as the active and idle output of machine, on the premise of the various constraints of system are met, make to give having for unit in total system
Work(is contributed, and adjustment is minimum, and mathematical modeling describes same formula (6), now bound g, g be controllable electric generator active power output and it is idle go out
The bound of power.
Step 8 is that mode Load flow calculation terminates in Fig. 1.
Although the present invention is disclosed as above with preferred embodiment, embodiment is not for limiting the present invention's.Not
In the spirit and scope for departing from the present invention, any equivalence changes done or retouching, the protection domain of the present invention is also belonged to.Cause
This protection scope of the present invention should be using the content that claims hereof is defined as standard.
Claims (5)
1. one kind solves the not convergent method of Load flow calculation after large scale electric network mode adjusts, it is characterised in that including following step
Suddenly：
1) remember that power system initial launch mode is S0, to S0Ways of carrying out adjustment operation, which adjustment operation is designated as F, really
Determine equipment composition and the node injection of each subnet after F is operated, calculate active amount of unbalance caused by each subnet operation, note jth
Active amount of unbalance caused by individual sub- net operation is Δ Pj,c, retain each subnet that can be connected by direct current, into step 2)；
2) whether judgment mode adjustment operation F can start emergent control measure, if will not start emergent control measure, it is determined that
The active amount of unbalance of each subnet is active amount of unbalance caused by each subnet operates, that is, has Δ Pj=Δ Pj,c, wherein Δ PjFor
The active amount of unbalance of j-th of subnet；
Otherwise, active amount of unbalance of each subnet caused by implementing when value measure is obtained respectively, and j-th of subnet of note is due to working as
It is Δ P that caused active amount of unbalance is implemented in value measurej,a, each subnet is due to active caused by value measure implementation
Amount of unbalance is when machine amount is cut in selection to be negative when selecting cutting load amount just；Determine that the active amount of unbalance of each subnet is each
Active amount of unbalance caused by subnet operation and due to implementing adding up for caused active amount of unbalance when value measure, that is, have Δ
Pj=Δ Pj,c+ΔPj,a；
If 3) the active amount of unbalance of each subnet is equal to 0, the unallocated amount of unbalance for counting each subnet is 0, and j-th of subnet of note is not
Distribution amount of unbalance is Δ Pj,t, into step 4)；
If the active amount of unbalance of each subnet is not equal to 0, to the active amount of unbalance of each subnet, according to the injustice being previously set
Weighing apparatus power analog, which shares rule and carries out simulation, to be shared, and the power balancing procedure after simulated operation, is obtained simulation and is shared rear each subnet
Unallocated amount of unbalance, into step 4)；
4) method of operation that the active amount of unbalance simulation of each subnet is shared after adjustment is designated as S1, to S1Carry out conventional Load Flow meter
The inspection for the small branch impedance data calculated, it is 1 that juxtaposition Newton-Raphson NR iteration round, which counts NT, when small branch impedance is less than
During the minimum impedance value of setting, then its minimum impedance value as setting is limited；
5) it is 1 that if Newton-Raphson NR iteration round, which counts NT, using the fixation PQ decomposition method numbers of setting, obtained by flat startup
Take the trend iterative initial value V of Newton-Raphson NR iteration(0), the tide as Newton-Raphson NR iteration in epicycle Load flow calculation
Iterative initial value is flowed, the flat startup refers to take all PQ node voltage amplitudes to be equal to 1.0, and PV node voltage is given setting electricity
Pressure value, phase angle are 0；
If Newton-Raphson NR iteration round, which counts NT, is more than 1, the optimal solution in previous Newton-Raphson NR iterative process is remembered
It is designated as V(k), judge corresponding optimal solution V(k)Newton-Raphson NR iterative powers deviation whether meet setting threshold value：It is if full
Foot, then used using adjustment trend initial value strategy restarting Load flow calculation, the trend iterative initial value of Newton-Raphson NR iteration
The flat startup changed after PQ decomposition method numbers solves the initial value V obtained(0)With in previous Newton-Raphson NR iterative process most
Excellent solution V(k)Weighting, i.e. V=α V(0)+βV(k), wherein α and β are the weight coefficient set and have alpha+beta=1, using V as epicycle trend
The trend iterative initial value of Newton-Raphson NR iteration in calculating；Become PQ decomposition method number strategies if not satisfied, then directly using,
And as Newton-Raphson NR iteration wheel numbers NT increases automatic increase PQ decomposition method numbers, initial value V is obtained by flat start(0), make
For the trend iterative initial value of Newton-Raphson NR iteration in epicycle Load flow calculation；
After the trend iterative initial value for determining Newton-Raphson NR iteration, the iterations of NT wheel Newton-Raphson NR iteration is put
K=1, into step 6)；
6) candidate's balancing machine is simulated in Newton-Raphson NR iteration, adjustable space or amount of capacity are pressed between candidate's balancing machine
The unallocated amount of unbalance or balancing machine for sharing each subnet are more limited the quantity, and voltage magnitude when calculating iteration and phase angle iteration decline
Subtract and compare index；
7) if Newton-Raphson NR iteration trend restrains, into step 8)；Otherwise, if epicycle Newton-Raphson NR iteration
Iterations K no more than the maximum times Kmax of setting, then make K increase 1, continue epicycle Newton-Raphson NR into step 6)
Iteration, if the iterations K of epicycle Newton-Raphson NR iteration is equal to the maximum iteration Kmax of setting, judge epicycle
The parameter NTmax whether number NT is more than setting is taken turns, next round newton-pressgang is carried out as made NT increasings 1 enter step 5) no more than if
Inferior NR iteration, such as larger than then according to the voltage magnitude and phase angle iteration in Newton-Raphson NR iterative process at index catastrophe point
Attenuation ratio index size determines Load Flow Solution morbid state feature, enters step 8) after carrying out automatic load flow rectification；
8) mode Load flow calculation terminates.
2. solving the not convergent method of Load flow calculation after large scale electric network mode adjusts as claimed in claim 1, its feature exists
In the mode of the step 1) adjusts operation F, including generator power adjustment, load power are set, load tap changer is set
And the operation including circuit/bus/plant stand stoppage in transit.
3. solving the not convergent method of Load flow calculation after large scale electric network mode adjusts as claimed in claim 1, its feature exists
In to the active amount of unbalance of each subnet in the step 3), sharing rule according to the imbalance power simulation being previously set and enter
Row simulation is shared, the power balancing procedure after simulated operation, obtains the unallocated amount of unbalance that rear each subnet is shared in simulation, including
Procedure below：
If selection is shared by simulation primary frequency modulation response, into step 3-1), if selection is shared by simulation frequency modulation frequency modulation response,
Then enter step 3-2), if neither selecting to share by simulation primary frequency modulation response nor selecting by simulation frequency modulation frequency modulation response point
Stand, then the unallocated amount of unbalance for counting each subnet is the active amount of unbalance of each subnet；
3-1) according to the difference coefficient of each unit, the current loads of each load bus and the change of static load frequency in each subnet
1% causes active change percentage, and each unit and each load bus are shared by formula (1) and (2) respectively：
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mi>j</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mi>j</mi>
</mrow>
</munder>
<msub>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>&Element;</mo>
<mi>j</mi>
</mrow>
</munder>
<msub>
<mi>P</mi>
<mrow>
<mi>l</mi>
<mi>k</mi>
<mn>0</mn>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>l</mi>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mi>j</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>l</mi>
<mi>k</mi>
<mn>0</mn>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mi>j</mi>
</mrow>
</munder>
<msub>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>&Element;</mo>
<mi>j</mi>
</mrow>
</munder>
<msub>
<mi>P</mi>
<mrow>
<mi>l</mi>
<mi>k</mi>
<mn>0</mn>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein：Ri,jFor the difference coefficient of i-th of unit in j-th of subnet, Plk0,jFor k-th load bus in j-th of subnet
Current loads, lk,jStatic load frequency change 1% for k-th of load bus in j-th of subnet causes active change percentage
Than Δ PGi,jWith Δ Plk,jThe imbalance that the simulation of i-th of unit and k-th of load bus is shared respectively in j-th of subnet
Amount；
If after sharing by formula (1), some unit output is out-of-limit, then fixes the unit output, represents the unit Primary frequency control ability
Failure, remainder is continued to share inside the unit and load bus still with Primary frequency control ability；
After sharing end, the unallocated amount of unbalance for putting each subnet is 0；
3-2) the active amount of unbalance according to each subnet, Automation generation control unit allocation logic is simulated, controlled according to first region
AGC shares in deviation ACE areas processed, and rear upper level scheduling administration AGC is participated in sharing, and finally AGC units are joined in other ACE areas in net
With the mode shared, share as follows：
By the active amount of unbalance of each subnet according to ACE regional extents, its affiliated ACE region is defined；
If for the active amount of unbalance of some subnet, power shortage such as formula (3) institute caused by some ACE region belonging to it
State, represent that the power shortage can be maintained to share balance by the AGC units in the ACE regions, then by the power shortage in the ACE
Share between AGC units in region, now frequency difference is 0：
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msub>
<mi>&Delta;P</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<mi>i</mi>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>j</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
2
Or
Wherein, Δ PiRepresent power shortage caused by i-th of ACE region belonging to the active amount of unbalance of some subnet, Pj0、WithCurrent output and the output upper limit and output lower limit of j-th of AGC unit in i-th of ACE region are represented respectively；
If for the active amount of unbalance of some subnet, power shortage such as formula (4) institute caused by some ACE region belonging to it
State, represent that the power shortage exceedes the adjustable extent of AGC units in the ACE regions, but still AGC units are administered in upper level degree
Adjustable extent, then AGC units in the ACE regions are shared into the upper limit or lower limit, overage is dispatched by upper level and managed
Linchpin AGC participates in sharing, and it is 0 finally to realize frequency difference：
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<mi>i</mi>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>j</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>s</mi>
<mo>&Element;</mo>
<mi>S</mi>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>s</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<msub>
<mi>&Delta;P</mi>
<mi>i</mi>
</msub>
<mo>></mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<mi>i</mi>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>j</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
Or
Wherein, S represents the set of upper level degree administration AGC units, Ps0、WithUpper level degree administration AGC machines are represented respectively
Current output and the output upper limit and output lower limit of s-th of AGC unit in group；
If for the active amount of unbalance of some subnet, power shortage such as formula (5) institute caused by some ACE region belonging to it
State, represent that the power shortage exceedes the adjustable extent of AGC units and the AGC units of schedule level one administration thereon in this ACE areas, then
AGC units in the ACE regions and the AGC units of schedule level one administration thereon are shared into the upper limit or lower limit, overage is by net
AGC units participate in droop control together in interior other ACE regions：
<mrow>
<msub>
<mi>&Delta;P</mi>
<mi>i</mi>
</msub>
<mo>></mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<mi>i</mi>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>j</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>j</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>s</mi>
<mo>&Element;</mo>
<mi>S</mi>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>s</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
Or
After sharing end, the unallocated amount of unbalance for putting each subnet is 0.
4. solving the not convergent method of Load flow calculation after large scale electric network mode adjusts as claimed in claim 3, its feature exists
In the step 6), specifically including following steps：
6-1) if the quantity of candidate's balancing machine is 0, into step 6-3)；
If the quantity of candidate's balancing machine is not 0, the unallocated amount of unbalance of such as each subnet is not equal to 0 and epicycle Newton-Raphson
The iterations K of NR iteration is equal to 1, then adjustable sky is pressed each subnet candidate balancing machine the unallocated amount of unbalance of each subnet
Between or amount of capacity share out equally, into step 6-2)；
If 6-2) meet what the accumulative active deviation of each node was less than setting value Psh and checked in Newton-Raphson NR iteration
Balancing machine, which is contributed, has out-of-limit deviation, then shares out equally by adjustable space or amount of capacity and more limit the quantity, into step 6-3)；
Voltage magnitude and phase angle iterative attenuation when 6-3) calculating Newton-Raphson NR iteration by the following method compare index：
If the correction [Δ x] during epicycle Newton-Raphson NR iteration kths(k)Middle voltage magnitude and phase angle greatest member it is exhausted
It is respectively Δ U to value(k)With Δ δ(k), by formulaWithCalculate voltage magnitude respectively and phase angle iterative attenuation ratio refers to
Mark, wherein Δ U(k-1)With Δ δ(k-1)Respectively epicycle Newton-Raphson NR iteration the K-1 times when correction [Δ x](k-1)In
The absolute value of voltage magnitude and phase angle greatest member.
5. solving the not convergent method of Load flow calculation after large scale electric network mode adjusts as claimed in claim 4, its feature exists
In being declined in the step 7) according to the voltage magnitude in Newton-Raphson NR iterative process at index catastrophe point and phase angle iteration
Subtract and determine Load Flow Solution morbid state feature than index size, carry out automatic load flow rectification, specifically include procedure below：
If each indexes suddenly changed point Δ U(k)/ΔU(k-1)The size of sum is more than Δ δ(k)/Δδ(k-1)The size of sum, then Load Flow Solution
Ill feature shows as local reactive power support deficiency, into step 7-1)；Otherwise it is disconnected to show as part for Load Flow Solution morbid state feature
Face trend is overweight, into step 7-2)：
The trend mode of local reactive power support deficiency 7-1) is shown as Load Flow Solution morbid state feature, by adjust generator without
Work(contributes, throws/cut electric capacity, reactor group, on the premise of the various constraints of system are met, makes to give the active of unit in total system
Adjustment of contributing is minimum, and its mathematical modeling is described such as formula (6)：
<mrow>
<mtable>
<mtr>
<mtd>
<mi>min</mi>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<msub>
<mi>S</mi>
<mi>G</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>w</mi>
<mrow>
<mi>P</mi>
<mi>G</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mn>0</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<munder>
<mi>g</mi>
<mo>&OverBar;</mo>
</munder>
<mo>&le;</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mover>
<mi>g</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein：SGRepresent the adjustable unit output set of total system, PG0iAnd PGiRespectively in the adjustable unit output set of total system
The front and rear active power of i platforms unit adjustment, wPGFor the active weight coefficient of corresponding unit；H (x) be decision variable x functions etc.
Formula constrains, and refers to the network and power-balance constraint of power system；Inequality constraints g (x) is by bound gThe operation of decision
Constraint, the bound gFor the idle bound contributed, throw/cut electric capacity and reactor group of controllable electric generator；
7-2) show as the overweight trend mode of section tidal current for Load Flow Solution morbid state feature, by adjust generated power and
Idle output, on the premise of the various constraints of system are met, make the active power output adjustment for giving unit in total system minimum, mathematics
Model describes same formula (6), now bound gFor controllable electric generator active power output and the bound of idle output.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710583770.7A CN107359616B (en) | 2017-07-18 | 2017-07-18 | Method for solving problem of non-convergence of load flow calculation after large-scale power grid mode adjustment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710583770.7A CN107359616B (en) | 2017-07-18 | 2017-07-18 | Method for solving problem of non-convergence of load flow calculation after large-scale power grid mode adjustment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107359616A true CN107359616A (en) | 2017-11-17 |
| CN107359616B CN107359616B (en) | 2020-05-19 |
Family
ID=60284245
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710583770.7A Active CN107359616B (en) | 2017-07-18 | 2017-07-18 | Method for solving problem of non-convergence of load flow calculation after large-scale power grid mode adjustment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107359616B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108281967A (en) * | 2018-01-04 | 2018-07-13 | 国网浙江省电力有限公司电力科学研究院 | A kind of algorithm startup method improving electric power system tide computational convergence |
| CN108365609A (en) * | 2018-03-26 | 2018-08-03 | 国网山东省电力公司淄博供电公司 | The iteration tidal current computing method of meter and interconnected network subregion Automatic Generation Control |
| CN108808683A (en) * | 2018-06-12 | 2018-11-13 | 国网福建省电力有限公司 | Alternating current-direct current power grid load flow calculation method containing multi-balancing machine and multimachine joint pressure regulation |
| CN110460061A (en) * | 2019-08-29 | 2019-11-15 | 国电南瑞科技股份有限公司 | Trend quick calculation method and system after a kind of failure |
| CN112381401A (en) * | 2020-11-13 | 2021-02-19 | 中国南方电网有限责任公司 | Power market power generation dispatching plan safety checking method |
| CN114239252A (en) * | 2021-12-06 | 2022-03-25 | 佛山电力设计院有限公司 | Power system operation mode generation method, system, computer and storage medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103425878A (en) * | 2013-08-01 | 2013-12-04 | 哈尔滨工业大学 | Method for rapidly calculating electrical power system quasi dynamic trend and power grid operation situation |
| CN103632046A (en) * | 2013-11-20 | 2014-03-12 | 国网河南省电力公司南阳供电公司 | Power grid load flow calculation method |
| CN104156609A (en) * | 2014-08-19 | 2014-11-19 | 河海大学 | Power grid flow solvability identification and initial value generating method based on continuous power flow |
| CN104732459A (en) * | 2015-03-31 | 2015-06-24 | 上海交通大学 | Large-scale power system ill-condition load flow analysis system |
-
2017
- 2017-07-18 CN CN201710583770.7A patent/CN107359616B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103425878A (en) * | 2013-08-01 | 2013-12-04 | 哈尔滨工业大学 | Method for rapidly calculating electrical power system quasi dynamic trend and power grid operation situation |
| CN103632046A (en) * | 2013-11-20 | 2014-03-12 | 国网河南省电力公司南阳供电公司 | Power grid load flow calculation method |
| CN104156609A (en) * | 2014-08-19 | 2014-11-19 | 河海大学 | Power grid flow solvability identification and initial value generating method based on continuous power flow |
| CN104732459A (en) * | 2015-03-31 | 2015-06-24 | 上海交通大学 | Large-scale power system ill-condition load flow analysis system |
Non-Patent Citations (1)
| Title |
|---|
| 张明广 等: "潮流计算中牛顿一拉夫逊法与PQ分解法结合使用探讨", 《机电信息》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108281967A (en) * | 2018-01-04 | 2018-07-13 | 国网浙江省电力有限公司电力科学研究院 | A kind of algorithm startup method improving electric power system tide computational convergence |
| CN108365609A (en) * | 2018-03-26 | 2018-08-03 | 国网山东省电力公司淄博供电公司 | The iteration tidal current computing method of meter and interconnected network subregion Automatic Generation Control |
| CN108365609B (en) * | 2018-03-26 | 2020-01-14 | 国网山东省电力公司淄博供电公司 | Unified iterative power flow calculation method considering automatic power generation control of interconnected power grid partitions |
| CN108808683A (en) * | 2018-06-12 | 2018-11-13 | 国网福建省电力有限公司 | Alternating current-direct current power grid load flow calculation method containing multi-balancing machine and multimachine joint pressure regulation |
| CN110460061A (en) * | 2019-08-29 | 2019-11-15 | 国电南瑞科技股份有限公司 | Trend quick calculation method and system after a kind of failure |
| CN110460061B (en) * | 2019-08-29 | 2022-08-23 | 国电南瑞科技股份有限公司 | Method and system for rapidly calculating power flow after fault |
| CN112381401A (en) * | 2020-11-13 | 2021-02-19 | 中国南方电网有限责任公司 | Power market power generation dispatching plan safety checking method |
| CN114239252A (en) * | 2021-12-06 | 2022-03-25 | 佛山电力设计院有限公司 | Power system operation mode generation method, system, computer and storage medium |
| CN114239252B (en) * | 2021-12-06 | 2022-10-11 | 佛山电力设计院有限公司 | Method and system for generating operation mode of power system, computer and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107359616B (en) | 2020-05-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107359616A (en) | One kind solves the not convergent method of Load flow calculation after large scale electric network mode adjusts | |
| CN104362648B (en) | A kind of idle phase modulation method of photovoltaic plant | |
| CN110970906B (en) | Method and system for coordinated control of reactive voltage of wind power plant | |
| CN108599259B (en) | Micro-grid active operation decision method based on sensitivity analysis | |
| CN110932281B (en) | Multi-section cooperative correction method and system based on quasi-steady-state sensitivity of power grid | |
| CN105490304A (en) | Coordination control method for a multiterminal flexible DC power grid to stabilize wind power fluctuations in a cross-regional manner | |
| CN109687497B (en) | Prevention control method and system for inhibiting extra-high voltage direct current continuous commutation failure | |
| CN110970907B (en) | Method and system for coordinated control of reactive voltage of photovoltaic power station | |
| CN105515012A (en) | Power storage participation power system voltage control method and device | |
| CN106487042A (en) | A kind of Multiple Time Scales micro-capacitance sensor voltage power-less optimized controlling method | |
| CN108711868A (en) | It is a kind of meter and islet operation voltage security GA for reactive power optimization planing method | |
| CN104767199A (en) | Directional control method for cut set cross section load flow in electric power system | |
| CN108390387A (en) | A kind of source lotus peak regulation control method of dynamic self-discipline decentralized coordinating | |
| CN110970905B (en) | Multi-power-supply reactive power coordination method and system for optimizing voltage control capability of wind power plant | |
| Morales-Paredes et al. | Cooperative control of power quality compensators in microgrids | |
| CN107451743B (en) | Active power distribution network island operation method considering source load operation characteristics | |
| CN113097995A (en) | Alternating current-direct current hybrid power distribution network optimized operation method considering multiple time scales | |
| CN109301833B (en) | AC/DC series-parallel power grid power transmission section overload control method based on path stripping | |
| Zhang et al. | Distributed predefined-time control for hybrid AC/DC microgrid | |
| CN111313461B (en) | Black start partitioning method based on step-by-step solving | |
| CN113131531B (en) | Adjustment standby sharing method and system suitable for different operation conditions of power grid | |
| CN107069699A (en) | A kind of optimal emergency control policy method for fast searching of power system | |
| CN114977316A (en) | Photovoltaic power station active power control method and system based on adjustable capacity margin | |
| CN112271724A (en) | Virtual power plant partition construction model and construction method based on voltage regulation | |
| CN110854853A (en) | Generating method and system for starting scheme of generator set |
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 |