CN106451459A - Power flow calculation method for accessing of large number of distributed power flow controllers - Google Patents
Power flow calculation method for accessing of large number of distributed power flow controllers Download PDFInfo
- Publication number
- CN106451459A CN106451459A CN201610855875.9A CN201610855875A CN106451459A CN 106451459 A CN106451459 A CN 106451459A CN 201610855875 A CN201610855875 A CN 201610855875A CN 106451459 A CN106451459 A CN 106451459A
- Authority
- CN
- China
- Prior art keywords
- power flow
- flow controller
- power
- distributed power
- dpfc
- 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
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a power flow calculation method for accessing of a large number of distributed power flow controllers. Alternatively iterated power flow calculation is utilized, thereby keeping an iteration form of a power flow algorithm and a characteristic of a Jacobian matrix. Furthermore accessing modes are flexibly changed according to different control objects of the distributed power flow controllers. Compared with an existing algorithm which comprises a single power flow controller simultaneous iteration control variable, the power flow calculation method has an advantage of no requirement for modification of the Jacobian matrix. Furthermore through setting a sub-iteration process, the power flow calculation method can be easily combined with an original power flow algorithm just through performing corresponding modification in the sub-iteration process. The power flow calculation method is particularly suitable for power flow calculation on the condition that a large number of distributed power flow controllers are accessed and different control strategies are utilized.
Description
Technical field
The present invention relates to a kind of flexible tidal current computing method accessing for a large amount of Distributed Power Flow controllers, belong to flexible
UHV AC transmission technique field.
Background technology
According to national energy and power planning, trans-regional Large Copacity extra-high voltage direct-current transmission, concentrate access extensive
Generation of electricity by new energy will obtain significant progress in China's power system, and after its access, the impact to operation of power networks and trend are flexible
The problems such as control, is also needed further investigation badly and explores.
In the system of no any tidal current controller, the electric network swim mainly impedance function according to system power line road etc.
And NATURAL DISTRIBUTION.But this self-assembling formation, uncontrolled trend distribution is not typically the trend desired by system operation
Distribution, and conveying capacity is not enough, loss is big, voltage out-of-limit to bring to the operation of system, or even system stability breaks
Bad a series of problems, such as.For improving system load flow distribution, current power system adopts some jumbo Survey of Flexible AC Transmission Systems,
But jumbo centralized flow controller is due to involving great expense and floor space is big, therefore Distributed Power Flow controller starts because of it
Unique advantage starts to cause the concern of industry.
The Load flow calculation containing power flow control devices at present mainly passes through system state variables, including node voltage, phase
Angle, with the state variable simultaneous iteration of power flow control devices it is therefore desirable to the element correction in Jacobian matrix, and in trend
In calculating, the control targe of power flow control devices is difficult to change.When containing Distributed Power Flow controller in system, lead to Ya Ke
Greatly increase than matrix dimension, lead to Load flow calculation complicated.
Content of the invention
Goal of the invention：The technical problem to be solved is that a large amount of Distributed Power Flow controllers access Load flow calculation,
A kind of tidal current computing method that can conveniently realize and accessing containing a large amount of Distributed Power Flow controllers is provided.
The technical solution adopted for the present invention to solve the technical problems is as follows：
A kind of tidal current computing method that a large amount of distribution flow controller accesses it is characterised in that this computational methods include as
Lower step：
(1) system initialization, determines number, on-position and the control mode of Distributed Power Flow controller, described control
Mode processed includes constant circuit active power controller strategy and constant output control strategy, and arranges at the beginning of its state variable
Value；
(2) power deviation of each node of computing system, repaiies according to the equivalent power injection model of Distributed Power Flow controller
The power deviation of positive articulare；
Equivalent injecting power mathematical model according to Distributed Power Flow controller：
Wherein, Ui, θi, Pis, QisThe respectively voltage magnitude of i-node, phase angle, equivalent injection active power, equivalent injection
Reactive power；Uj、θj、Pjs、QjsIt is respectively the voltage magnitude of j node, phase angle, equivalent injection active power, equivalent injection is idle
Power g, b are respectively the conductance of circuit and admittance, b between i and j nodecFor the admittance over the ground of i and j node two ends, UDPFC、δDPoint
The voltage of Distributed Power Flow controller between the i that Wei not access, j two node and phase angle；
(3) main iterative process adopts Newton-Laphson method solving system state variable, obtains the Ya Ke of system state variables
Ratio passes through the running status of matrixing solving system after matrix, including voltage and phase angle, and carry out changing of system state variables
In generation, is revised；
(4) following sub- iterative process, the control targe equation according to Distributed Power Flow controller and internal constraint side are entered
Journey is modified to its state variable value, if the condition of convergence being not reaching in sub- iteration, repeats step (4), if
Reach, then enter step (5)；
The internal constraint equation of wherein Distributed Power Flow controller is：
PDPFCThe active power sending for Distributed Power Flow controller, its control strategy has constant line road active power controller
Strategy and constant output control strategy；Adopt suitable control strategy according to the needs that system is controlled, to control line
When road active power is steady state value, using constant circuit active power controller strategy, its control targe equation is：
PrefFor the circuit active power reference value setting, sub- iteration termination condition is max | Δ Pij,PDPFC|≤ε or iteration
Number of times k>100, wherein max | Δ Pij,PDPFC| represent active power of output and circuit wattful power in each Distributed Power Flow controller
Maximum in rate deviation；
If desired Distributed Power Flow controller is carried out when constant output controls using constant output as controlling
Strategy, its control targe is：
QrefCapacity for Distributed Power Flow controller sets reference value, and sub- iteration convergence condition is max | Δ QDPFC,PDPFC|
≤ ε or iterationses k>100, wherein max | Δ QDPFC,PDPFC| represent each Distributed Power Flow controller in output deviation and
Maximum in active power；
(5) judge whether to meet system convergence condition, if so, then terminate Load flow calculation；If it is not, then going to step (2)；
Wherein system convergence condition is max | Δ P, Δ Q | ＜ ε or main iterationses m >=100；Δ P, Δ Q are each node
Active power and the deviation of reactive power.
Further, in step (1), it is perunit value that the voltage magnitude initial value of Distributed Power Flow controller all takes 0.1p.u.
For 0.1, phase angle takesIn step (4), (5), ε takes 10-5.
Beneficial effect：The present invention proposes the tidal current computing method that a kind of a large amount of Distributed Power Flow controller accesses, the party
Method is by the state variable of Distributed Power Flow controller and system variable alternating iteration, and considers Distributed Power Flow controller not
Same control targe, realizes Load flow calculation when accessing containing a large amount of Distributed Power Flow controllers.With existing conventional simultaneous iteration tide
Flow algorithm is compared, and can retain traditional iteration form of power flow algorithm and the feature of Jacobian matrix using this method, not increase refined
Gram than matrix dimension it is only necessary to accordingly be changed it is easy to combine with original power flow algorithm in sub- iterative process.
Brief description
Fig. 1 is the general flow chart of the inventive method.
Fig. 2 is IEEE-30 node system structure chart.
Specific embodiment：
It is described in further detail with reference to the accompanying drawings and in conjunction with the embodiments to the present invention.But the invention is not restricted to
The example going out.
The present embodiment is to do corresponding Load flow calculation to test after Distributed Power Flow controller accesses only to 1 and 2 DPFC
The correctness of card algorithm, if desired additionally increases the quantity of Distributed Power Flow controller it is only necessary to increase distributed tide accordingly
The sub- iterative process of stream controller.
(1) system initialization, main iteration and sub- iterationses m, k are disposed as 0, determine the individual of Distributed Power Flow controller
Number, on-position and control mode, described control mode includes constant circuit active power controller strategy and constant output work(
Rate control strategy, and its state variable initial value voltage magnitude is set for 0.1p.u., phase angle takes
(2) power deviation of each node of computing system, repaiies according to the equivalent power injection model of Distributed Power Flow controller
The power deviation of positive articulare；
Equivalent injecting power mathematical model according to Distributed Power Flow controller：
Wherein, Ui, θi, Pis, QisThe respectively voltage magnitude of i-node, phase angle, equivalent injection active power, equivalent injection
Reactive power；Uj、θj、Pjs、QjsIt is respectively the voltage magnitude of j node, phase angle, equivalent injection active power, equivalent injection is idle
Power g, b are respectively the conductance of circuit and admittance, b between i and j nodecFor the admittance over the ground of i and j node two ends, UDPFC、δDPoint
The voltage of Distributed Power Flow controller between the i that Wei not access, j two node and phase angle；
(3) main iterative process adopts Newton-Laphson method solving system state variable, obtains the Ya Ke of system state variables
Ratio passes through the running status of matrixing solving system after matrix, including voltage and phase angle, and carry out changing of system state variables
In generation, is revised；
(4) following sub- iterative process, the control targe equation according to Distributed Power Flow controller and internal constraint side are entered
Journey is modified to its state variable value, if the condition of convergence being not reaching in sub- iteration, repeats step (4), if
Reach, then enter step (5)；
The internal constraint equation of wherein Distributed Power Flow controller is：
PDPFCThe active power sending for Distributed Power Flow controller, its control strategy has constant line road active power controller
Strategy and constant output control strategy；Adopt suitable control strategy according to the needs that system is controlled, to control line
When road active power is steady state value, using constant circuit active power controller strategy, its control targe equation is：
PrefFor the circuit active power reference value setting, sub- iteration termination condition is max | Δ Pij,PDPFC|≤ε or iteration
Number of times k>100, wherein max | Δ Pij,PDPFC| represent active power of output and circuit wattful power in each Distributed Power Flow controller
Maximum in rate deviation；；
If desired Distributed Power Flow controller is carried out when constant output controls using constant output as controlling
Strategy, its control targe is：
QrefCapacity for Distributed Power Flow controller sets reference value, and sub- iteration convergence condition is max | Δ QDPFC,PDPFC|
≤ ε or iterationses k>100, wherein max | Δ QDPFC,PDPFC| represent each Distributed Power Flow controller in output deviation and
Maximum in active power；
(5) judge whether to meet system convergence condition, if so, then terminate Load flow calculation；If it is not, then going to step (2).
Wherein system convergence condition is max | Δ P, Δ Q | ＜ ε or main iterationses m >=100；Δ P, Δ Q are each node
Active power and the deviation of reactive power.
In step (4), (5), ε takes 10-5.
As shown in Fig. 2 according to the network structure of IEEE-30 node power system model and network parameter, this network includes 6
Platform electromotor (node 1,2,5,8,11,13) and 41 branch roads.Table 1 sets forth the initial fortune of branch road 2-6 and branch road 9-10
Row result, table 2 sets forth the situation installing DPFC on branch road 2-6 and branch road 9-10, and DPFC is separately mounted to above-mentioned branch road
First section exit.
Table 1 IEEE-30 node system partial branch initial operating state
For verifying the correctness of Load flow calculation, this case study on implementation assumes three kinds of methods of operation.
The method of operation one：By installing DPFC in branch road 2,6 first section, controlling brancher transmits active reference value Pref=0.59,
Line transmission power is made to increase 50%.
The method of operation two：By installing DPFC in branch road 9,10 first section, controlling brancher transmits active reference value Pref=
0.25 so that line transmission power reduction 25%.
The method of operation three：By installing DPFC in branch road 2-6 and branch road 9-10 simultaneously, controlling brancher transmits active reference value
P26ref=0.4, P910ref=0.3.
Table 2 IEEE-30 node system associated branch adds the power flow solutions after DPFC
Calculation of tidal current in table 2 shows, under mode one, mode two and mode three, Distributed Power Flow controller is equal
Branch Power Flow can be made to reach the reference value of control, thus demonstrating the reasonability of this algorithm.
The above is only the preferred embodiment of the present invention it is noted that ordinary skill people for the art
For member, under the premise without departing from the principles of the invention, some improvements and modifications can also be made, these improvements and modifications also should
It is considered as protection scope of the present invention.
Claims (2)
1. a kind of a large amount of distribution flow controller accesses tidal current computing method it is characterised in that this computational methods include as follows
Step：
(1) system initialization, determines number, on-position and the control mode of Distributed Power Flow controller, described controlling party
Formula includes constant circuit active power controller strategy and constant output control strategy, and arranges its state variable initial value；
(2) power deviation of each node of computing system, according to the equivalent power injection model correction phase of Distributed Power Flow controller
The power deviation of articulare；
Equivalent injecting power mathematical model according to Distributed Power Flow controller：
Wherein, Ui, θi, Pis, QisIt is respectively the voltage magnitude of i-node, phase angle, equivalent injection active power, equivalent injection is idle
Power；Uj、θj、Pjs、QjsThe respectively voltage magnitude of j node, phase angle, equivalent injection active power, equivalent injection reactive power
G, b are respectively the conductance of circuit and admittance, b between i and j nodecFor the admittance over the ground of i and j node two ends, UDPFC、δDIt is respectively
I, the voltage of Distributed Power Flow controller between j two node and the phase angle accessing；
(3) main iterative process adopts Newton-Laphson method solving system state variable, obtains the Jacobi square of system state variables
By the running status of matrixing solving system after battle array, including voltage and phase angle, and carry out the iteration of system state variables and repair
Just；
(4) following sub- iterative process, the control targe equation according to Distributed Power Flow controller and internal constraint equation pair are entered
Its state variable value is modified, if the condition of convergence being not reaching in sub- iteration, repeats step (4), if reaching,
Then enter step (5)；
The internal constraint equation of wherein Distributed Power Flow controller is：
PDPFCThe active power sending for Distributed Power Flow controller, its control strategy has constant line road active power controller strategy
And constant output control strategy；Adopt suitable control strategy according to the needs that system is controlled, have to control circuit
When work(power is steady state value, using constant circuit active power controller strategy, its control targe equation is：
PrefFor the circuit active power reference value setting, sub- iteration termination condition is max | Δ Pij,PDPFC|≤ε or iterationses
k>100, wherein max | Δ Pij,PDPFC| represent that in each Distributed Power Flow controller, active power of output and circuit active power are inclined
Maximum in difference；
If desired Distributed Power Flow controller is carried out when constant output controls, adopting constant output to be control strategy,
Its control targe is：
QrefCapacity for Distributed Power Flow controller sets reference value, and sub- iteration convergence condition is max | Δ QDPFC,PDPFC|≤ε
Or iterationses k>100, wherein max | Δ QDPFC,PDPFC| represent output deviation and active in each Distributed Power Flow controller
Maximum in power；
(5) judge whether to meet system convergence condition, if so, then terminate Load flow calculation；If it is not, then going to step (2)；
Wherein system convergence condition is max | Δ P, Δ Q | ＜ ε or main iterationses m >=100；Δ P, Δ Q are the active of each node
Power and the deviation of reactive power.
2. the tidal current computing method that a large amount of distribution flow controller according to claim 1 accesses, is characterized in that：Step
(1), in, the voltage magnitude initial value of Distributed Power Flow controller all takes 0.1p.u., and phase angle takesIn step (4), (5), ε takes 10-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610855875.9A CN106451459B (en) | 2016-09-27 | 2016-09-27 | A kind of tidal current computing method of a large amount of Distributed Power Flow controller accesses |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610855875.9A CN106451459B (en) | 2016-09-27 | 2016-09-27 | A kind of tidal current computing method of a large amount of Distributed Power Flow controller accesses |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106451459A true CN106451459A (en) | 2017-02-22 |
| CN106451459B CN106451459B (en) | 2019-01-18 |
Family
ID=58170526
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610855875.9A Active CN106451459B (en) | 2016-09-27 | 2016-09-27 | A kind of tidal current computing method of a large amount of Distributed Power Flow controller accesses |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106451459B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107425526A (en) * | 2017-06-12 | 2017-12-01 | 武汉理工大学 | Distributed Power Flow controller one phase equivalent emulation mode |
| CN109245105A (en) * | 2018-09-06 | 2019-01-18 | 国电南瑞科技股份有限公司 | A kind of Distributed Power Flow control method and system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103107559A (en) * | 2013-02-06 | 2013-05-15 | 武汉理工大学 | Method of confirming parameters of distributed power flow controller system |
| CN105610158A (en) * | 2016-02-06 | 2016-05-25 | 武汉理工大学 | Distributed power flow controller and control method therefor |
| CN105896545A (en) * | 2016-06-16 | 2016-08-24 | 国网江苏省电力公司电力科学研究院 | Accident chain mode based distributed power flow controller location method |
-
2016
- 2016-09-27 CN CN201610855875.9A patent/CN106451459B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103107559A (en) * | 2013-02-06 | 2013-05-15 | 武汉理工大学 | Method of confirming parameters of distributed power flow controller system |
| CN105610158A (en) * | 2016-02-06 | 2016-05-25 | 武汉理工大学 | Distributed power flow controller and control method therefor |
| CN105896545A (en) * | 2016-06-16 | 2016-08-24 | 国网江苏省电力公司电力科学研究院 | Accident chain mode based distributed power flow controller location method |
Non-Patent Citations (2)
| Title |
|---|
| G. MADHUSUDHANA RAO 等: "Design of a neural network based distributed power flow controller (DPFC) for power system stability", 《SIGNAL PROCESSING, COMMUNICATION, POWER AND EMBEDDED SYSTEM (SCOPES), 2016 INTERNATIONAL CONFERENCE ON》 * |
| 唐爱红 等: "分布式潮流控制器对***功率控制的研究", 《电力***保护与控制》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107425526A (en) * | 2017-06-12 | 2017-12-01 | 武汉理工大学 | Distributed Power Flow controller one phase equivalent emulation mode |
| CN107425526B (en) * | 2017-06-12 | 2019-07-23 | 武汉理工大学 | Distributed Power Flow controller one phase equivalent emulation mode |
| CN109245105A (en) * | 2018-09-06 | 2019-01-18 | 国电南瑞科技股份有限公司 | A kind of Distributed Power Flow control method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106451459B (en) | 2019-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106374513B (en) | A kind of more microgrid dominant eigenvalues optimization methods based on leader-followers games | |
| CN106532711B (en) | Change the Newton load flow calculation method of Jacobian matrix with iteration and node type | |
| CN110504691A (en) | It is a kind of meter and VSC control mode alternating current-direct current power distribution network optimal load flow calculation method | |
| CN108429294A (en) | A kind of the AC/DC network tide model and method for solving of energy content router | |
| CN108493940A (en) | A kind of Ac/dc Power Systems trend alternating iteration computational methods containing flexible direct current | |
| CN104899396B (en) | A kind of algorithm quicksort tidal current computing method of correction factor matrix | |
| CN108988344B (en) | Load flow calculation optimization method of power system | |
| CN111969609B (en) | Second-order cone optimal power flow model construction method and device for alternating current-direct current power transmission network | |
| CN114336674B (en) | Distributed toughness frequency control method for alternating-current micro-grid | |
| CN106712044A (en) | Multi-objective optimization method of reactive voltage in wind power plant | |
| CN110535189A (en) | Distributed generation resource low voltage traversing control method and system | |
| CN106451459A (en) | Power flow calculation method for accessing of large number of distributed power flow controllers | |
| CN108233408B (en) | A kind of MMC-MTDC system self-adaption droop control method | |
| CN106099996A (en) | The conforming control method of virtual synchronous electromotor containing state predictor | |
| CN107332277B (en) | Active power distribution network island operation method considering source load storage operation characteristics | |
| CN109638871B (en) | Main network dividing method of large-scale alternating current-direct current hybrid system considering wind power access | |
| CN104901313B (en) | Microgrid Voltage unbalance compensation method based on the Gossip communications of algorithms | |
| CN108347057B (en) | LCC-MMC mixed direct-current power grid alternating iteration load flow calculation method | |
| CN106410811B (en) | Iteration small impedance branches endpoint changes the tidal current computing method of Jacobian matrix for the first time | |
| CN105896579B (en) | Consider the charging pile planing method that distribution voltage dip causes electric car tripping to cause overvoltage | |
| CN110829448A (en) | Distributed voltage interactive support control method and system for alternating current-direct current hybrid power distribution network | |
| Guo et al. | Multi-objective planning for voltage sag compensation of sparse distribution networks with unified power quality conditioner using improved NSGA-III optimization | |
| CN106682363B (en) | The one type sagging control isolated island micro-capacitance sensor tidal current computing method that this is decomposed to moral | |
| Aguirre et al. | Improved power flow program for unbalanced radial distribution systems including voltage dependent loads | |
| CN106712029B (en) | The Newton load flow calculation method of small impedance branches PQ endpoint change Jacobian matrix |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: 210009 No. 87 Dingjiaqiao, Gulou District, Nanjing City, Jiangsu Province Applicant after: Southeast University Address before: 210096 No. 2 Southeast University Road, Jiangning District, Nanjing City, Jiangsu Province Applicant before: Southeast University |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |