CN109687439A - A kind of quick calculation method of power transmission network Alloy White Iron - Google Patents

A kind of quick calculation method of power transmission network Alloy White Iron Download PDF

Info

Publication number
CN109687439A
CN109687439A CN201811597769.0A CN201811597769A CN109687439A CN 109687439 A CN109687439 A CN 109687439A CN 201811597769 A CN201811597769 A CN 201811597769A CN 109687439 A CN109687439 A CN 109687439A
Authority
CN
China
Prior art keywords
cyclization
node
subnet
white iron
impedance
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.)
Pending
Application number
CN201811597769.0A
Other languages
Chinese (zh)
Inventor
郭英杰
张东院
武振华
冯长忠
朱英刚
杨锋
吕明超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Integrated Electronic Systems Lab Co Ltd
Original Assignee
Integrated Electronic Systems Lab Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Integrated Electronic Systems Lab Co Ltd filed Critical Integrated Electronic Systems Lab Co Ltd
Priority to CN201811597769.0A priority Critical patent/CN109687439A/en
Publication of CN109687439A publication Critical patent/CN109687439A/en
Pending legal-status Critical Current

Links

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/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 present invention relates to a kind of quick calculation methods of power transmission network Alloy White Iron, obtain electric network model and presence estimated data from area power grid EMS system, carry out electric network swim calculating, obtain the flow data of power grid before cyclization；According to switch state in power grid and equipment connecting relation, search obtains the All Paths between two node of cyclization；Using electric network model parameter and cyclization All Paths, forms cyclization and calculate the cyclization calculating subnet model needed；Subnet model is calculated according to cyclization, node admittance matrix is calculated, invert to obtain subnet nodal impedance matrix Z；Using subnet nodal impedance matrix Z and cyclization two o'clock position, loop circuit impedance is calculated；Using the flow data and loop circuit impedance of power grid before cyclization, Alloy White Iron is calculated.It is related that impedance matrix calculating is narrowed down to only close-loop path by the present invention, has greatly reduced calculation scale, solves the problems, such as when the calculating of fairly large power transmission network Alloy White Iron that speed is slow, low efficiency, improves work efficiency.

Description

A kind of quick calculation method of power transmission network Alloy White Iron
Technical field
The invention belongs to dispatching of power netwoks control technology field more particularly to a kind of quick calculating sides of power transmission network Alloy White Iron Method.
Background technique
To improve power supply reliability, area power grid is frequently necessary to carry out Loop Closing Operation.Because there is electricity in cyclization switch two sides The reasons such as pressure difference, phase difference can generate Alloy White Iron in operating process, influence the stable operation of power grid.Therefore, Loop Closing Operation it Before, dispatcher will do it Alloy White Iron calculating.
Traditional approach calculates Alloy White Iron, and main flow includes: Load flow calculation, the calculating of the whole network impedance matrix, loop circuit impedance It calculates, Alloy White Iron calculates.It is the power grid rule by inverting to obtain to the whole network admittance matrix when conventional method computing impedance matrix Mould is bigger, and admittance matrix order is higher, and calculating speed is slower.When power grid is larger, amount of calculation is very big, and speed is very Slowly, working efficiency is seriously affected.
Summary of the invention
For the above-mentioned problems in the prior art, the invention discloses a kind of quick calculating of power transmission network Alloy White Iron Method solves the problems, such as when the calculating of current Alloy White Iron that speed is slow, low efficiency.The technical solution adopted in the present invention is as follows:
A kind of quick calculation method of power transmission network Alloy White Iron, comprising the following steps:
Step 1 obtains electric network model and presence estimated data from area power grid EMS system, carries out electric network swim meter It calculates, obtains the flow data of power grid before cyclization；
Step 2, according to switch state in power grid and equipment connecting relation, search obtains whole roads between two node of cyclization Diameter；
Step 3 utilizes electric network model parameter and cyclization All Paths, the cyclization calculating subnet mould of composition cyclization calculating needs Type；
Node admittance matrix is calculated according to cyclization calculating subnet model in step 4, inverts to obtain subnet node impedance square Battle array Z；
Step 5, using subnet nodal impedance matrix Z and cyclization two o'clock position, calculate loop circuit impedance；
Step 6, flow data and loop circuit impedance using power grid before cyclization, calculate Alloy White Iron.
The invention has the following advantages that
It is related that impedance matrix calculating is narrowed down to only close-loop path by the present invention, has greatly reduced calculation scale, solve compared with The problem of speed is slow when extensive power transmission network Alloy White Iron calculates, low efficiency, improves work efficiency, and has and high applies valence Value.
Detailed description of the invention
Fig. 1 is logical procedure diagram of the invention.
Specific embodiment
With reference to the accompanying drawing, specific implementation step of the invention is described in detail.The present embodiment is based on area power grid EMS system, it can be achieved that power transmission network Alloy White Iron quick computing function.
Step 1 obtains electric network model and presence estimated data from area power grid EMS system, carries out electric network swim meter It calculates, obtains the flow data of power grid before cyclization；
Area power grid EMS system provides whole interfaces, includes the parameter of generator, transformer, route etc., and online Data of State Estimation etc..
Step 2, according to switch state in power grid and equipment connecting relation, search obtains whole roads between two node of cyclization Diameter.
Regard the whole network electric network model as non-directed graph, tie point regards the vertex of figure as, and equipment connecting relation regards the branch of figure as, Using graph theory search depth priority algorithm, search obtains the All Paths between two node of cyclization.Specific step is as follows:
S2.1, from cyclization starting point according to model connection relationship the conductive equipment of all connections is found first；
S2.2, again from conductive equipment, find the node that the equipment does not access, so use Depth Priority Algorithm Continue to traverse；
S2.3, when traversing cyclization terminal, an as paths.
During deep search, when to encounter conductive equipment be switch tool, judge the division state of switch, if close, Another side gusset of switch tool is then continued searching, is otherwise retracted.So search finishes, and can find All Paths.
Step 3 utilizes electric network model parameter and cyclization All Paths, the cyclization calculating subnet mould of composition cyclization calculating needs Type.
Electric network model mainly includes transformer, route etc., and parameter includes impedance etc..Model on path all counts subnet Model, connection relationship are constant.Forming cyclization calculating subnet model, specific step is as follows:
S3.1, by cyclization All Paths equipment and node separately constitute set of fingers and node set, member in set Element guarantees uniquely, duplicate to be added without；
All switch tools in S3.2, deletion set of fingers, the node contraction of switch tool two sides is at a node, together Shi Gengxin set of fingers and node set；
S3.3, node set have n element, then all nodes redefine section from 1 to n according to the sequence in node set Point number；Set of fingers is updated simultaneously.
Node admittance matrix is calculated according to cyclization calculating subnet model in step 4, inverts to obtain subnet node impedance square Battle array Z.
Subnet model connection relationship and impedance parameter are calculated according to cyclization and calculates admittance matrix, using Gaussian elimination method to leading It receives matrix inversion, calculates subnet nodal impedance matrix Z, the specific steps are as follows:
S4.1, subnet node admittance matrix Y [n, n] is established, ranks number is all n；
S4.2, traversal all branches of set of fingers, branch admittance calculate self-admittance and transadmittance, are separately added into Y matrix；
S4.3, Y matrix inversion obtain subnet nodal impedance matrix Z.
Step 5, using subnet nodal impedance matrix Z and cyclization two o'clock position, calculate loop circuit impedance.
Subnet nodal impedance matrix Z is that the relationship of impedance is presented in a manner of matrix, ZijIt is between node i and j Mutual impedance.According to ranks value of the cyclization point in impedance matrix, the loop circuit impedance Z between two node of cyclization can be calculatedeq=Zii +Zjj-2×Zij。
Step 6, flow data and loop circuit impedance using power grid before cyclization, calculate Alloy White Iron.
Electric network swim data before existing cyclization, the voltage difference including cyclization point both ends With And loop circuit impedance Zeq, it is directly substituted into Alloy White Iron calculation formula, obtains Alloy White Iron.
Cyclization dash current imax=2um/|Zeq|；
Cyclization steady-state current
Wherein it isImpedance angle, ω are 2 π f, f 50Hz.
The ground branch of area power grid is capacitive ground branch, and when engineering calculation is usually ignored.By nodal impedance matrix member The physical significance of element can be seen that subnet cyclization impedance and the whole network cyclization impedance computation result is almost the same, therefore subnet cyclization The calculated result that electric current calculated result and the whole network Alloy White Iron calculate is almost the same, and calculated result of the invention can be fully met Using needs.

Claims (8)

1. a kind of quick calculation method of power transmission network Alloy White Iron, which comprises the following steps:
Step 1 obtains electric network model and presence estimated data from area power grid EMS system, carries out electric network swim calculating, obtains The flow data of power grid before to cyclization；
Step 2, according to switch state in power grid and equipment connecting relation, search obtains the All Paths between two node of cyclization；
Step 3 utilizes electric network model parameter and cyclization All Paths, the cyclization calculating subnet model of composition cyclization calculating needs；
Node admittance matrix is calculated according to cyclization calculating subnet model in step 4, inverts to obtain subnet nodal impedance matrix Z；
Step 5, using subnet nodal impedance matrix Z and cyclization two o'clock position, calculate loop circuit impedance；
Step 6, flow data and loop circuit impedance using power grid before cyclization, calculate Alloy White Iron.
2. a kind of quick calculation method of power transmission network Alloy White Iron according to claim 1, which is characterized in that step 2 institute The method that the search stated obtains the All Paths between two node of cyclization is: regarding the whole network electric network model as non-directed graph, tie point is seen The vertex of figure is done, equipment connecting relation regards the branch of figure as, and using graph theory search depth priority algorithm, search obtains two section of cyclization All Paths between point.
3. a kind of quick calculation method of power transmission network Alloy White Iron according to claim 2, which is characterized in that step 2 institute The search stated obtains the All Paths between two node of cyclization, and specific step is as follows:
S2.1, from cyclization starting point according to model connection relationship the conductive equipment of all connections is found first；
S2.2, again from conductive equipment, find the node that the equipment does not access, so continued using Depth Priority Algorithm Traversal；
S2.3, when traversing cyclization terminal, an as paths；
During deep search, when to encounter conductive equipment be switch tool, judge the division state of switch, if close, then after Continuous another side gusset of search switch disconnecting link, is otherwise retracted.
4. a kind of quick calculation method of power transmission network Alloy White Iron according to claim 1-3, feature exist In the cyclization that composition cyclization calculates needs described in step 3 calculates subnet model, and specific step is as follows:
S3.1, by cyclization All Paths equipment and node separately constitute set of fingers and node set, in set element protect Card is unique, duplicate to be added without；
S3.2, all switch tools in set of fingers are deleted, the node contraction of switch tool two sides is at a node, while more New set of fingers and node set；
S3.3, node set have n element, then all nodes redefine node from 1 to n according to the sequence in node set Number；Set of fingers is updated simultaneously.
5. a kind of quick calculation method of power transmission network Alloy White Iron according to claim 4, which is characterized in that step 4 institute The method for obtaining subnet nodal impedance matrix Z stated is: calculating subnet model connection relationship according to cyclization and impedance parameter calculates Admittance matrix inverts to admittance matrix using Gaussian elimination method, calculates subnet nodal impedance matrix Z.
6. a kind of quick calculation method of power transmission network Alloy White Iron according to claim 5, which is characterized in that step 4 institute Specific step is as follows by the subnet nodal impedance matrix Z that obtains stated:
S4.1, subnet node admittance matrix Y [n, n] is established, ranks number is all n；
S4.2, traversal all branches of set of fingers, branch admittance calculate self-admittance and transadmittance, are separately added into Y matrix；
S4.3, Y matrix inversion obtain subnet nodal impedance matrix Z.
7. a kind of quick calculation method of power transmission network Alloy White Iron according to claim 6, which is characterized in that step 5 institute The method for calculating loop circuit impedance stated is:
In subnet nodal impedance matrix Z, ZijIt is the mutual impedance between node i and j, according to ranks of the cyclization point in impedance matrix Value, calculates the loop circuit impedance Z between two node of cyclizationeq=Zii+Zjj-2×Zij。
8. a kind of quick calculation method of power transmission network Alloy White Iron according to claim 7, which is characterized in that step 6 institute That states, which calculates the specific method of Alloy White Iron, is:
Electric network swim data before existing cyclization, the voltage difference including cyclization point both ends And loop Impedance Zeq, it is directly substituted into Alloy White Iron calculation formula, obtains Alloy White Iron:
Cyclization dash current imax=2um/|Zeq|；
Cyclization steady-state current
Wherein it isImpedance angle, ω are 2 π f, f 50Hz.
CN201811597769.0A 2018-12-25 2018-12-25 A kind of quick calculation method of power transmission network Alloy White Iron Pending CN109687439A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811597769.0A CN109687439A (en) 2018-12-25 2018-12-25 A kind of quick calculation method of power transmission network Alloy White Iron

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811597769.0A CN109687439A (en) 2018-12-25 2018-12-25 A kind of quick calculation method of power transmission network Alloy White Iron

Publications (1)

Publication Number Publication Date
CN109687439A true CN109687439A (en) 2019-04-26

Family

ID=66189458

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811597769.0A Pending CN109687439A (en) 2018-12-25 2018-12-25 A kind of quick calculation method of power transmission network Alloy White Iron

Country Status (1)

Country Link
CN (1) CN109687439A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110994610A (en) * 2019-12-24 2020-04-10 国网天津市电力公司 Monitoring information assistance-based loop closing method and device for power distribution network line
CN111697571A (en) * 2020-05-26 2020-09-22 云南电网有限责任公司玉溪供电局 Closed loop current calculation method during closed loop operation of power distribution network
CN112531717A (en) * 2020-12-09 2021-03-19 深圳供电局有限公司 Method and device for calculating closed loop power flow of power system
CN114221342A (en) * 2021-11-30 2022-03-22 国网北京市电力公司 Closed loop current calculation method and device, storage medium and electronic equipment

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101431255A (en) * 2008-11-17 2009-05-13 国电南瑞科技股份有限公司 Closed loop operation risk analysis method for power system
CN103324744A (en) * 2013-07-03 2013-09-25 国家电网公司 Topology search method based on distribution looped network automatic mapping
CN107834545A (en) * 2017-11-13 2018-03-23 国网四川省电力公司成都供电公司 A kind of city 110kV power network cyclization methods based on transfer load method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101431255A (en) * 2008-11-17 2009-05-13 国电南瑞科技股份有限公司 Closed loop operation risk analysis method for power system
CN103324744A (en) * 2013-07-03 2013-09-25 国家电网公司 Topology search method based on distribution looped network automatic mapping
CN107834545A (en) * 2017-11-13 2018-03-23 国网四川省电力公司成都供电公司 A kind of city 110kV power network cyclization methods based on transfer load method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
裴星宇等: "基于公共信息模型拓扑收缩的合环电流计算方法", 《黑龙江电力》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110994610A (en) * 2019-12-24 2020-04-10 国网天津市电力公司 Monitoring information assistance-based loop closing method and device for power distribution network line
CN111697571A (en) * 2020-05-26 2020-09-22 云南电网有限责任公司玉溪供电局 Closed loop current calculation method during closed loop operation of power distribution network
CN112531717A (en) * 2020-12-09 2021-03-19 深圳供电局有限公司 Method and device for calculating closed loop power flow of power system
CN112531717B (en) * 2020-12-09 2023-07-21 深圳供电局有限公司 Method and device for calculating loop closing power flow of power system
CN114221342A (en) * 2021-11-30 2022-03-22 国网北京市电力公司 Closed loop current calculation method and device, storage medium and electronic equipment
CN114221342B (en) * 2021-11-30 2024-04-26 国网北京市电力公司 Loop closing current calculation method and device, storage medium and electronic equipment

Similar Documents

Publication Publication Date Title
CN109687439A (en) A kind of quick calculation method of power transmission network Alloy White Iron
CN105260516B (en) A kind of electromagnetical transient emulation method of the sub-network containing switching characteristic
CN104462661B (en) Electronic power switch interpolation real-time emulation method
CN108879664B (en) AC/DC system online voltage stability evaluation method based on wide area measurement
CN102403720B (en) Hyper-real-time setting method for superposed time sequences based on transient voltage safety margin
CN108923401B (en) A kind of calculation method and device of superconductive current limiter fault resstance and cut-off current
CN108923396B (en) A kind of short circuit current quick calculation method of multiterminal flexible direct current power grid
CN108229021A (en) Modularization multi-level converter locking modeling method based on Real Time Digital Simulator
CN102750416A (en) Topological subnetting method of electromagnetic transient simulation containing switching characteristic circuit
CN103902814A (en) Electric power system operation state detecting method based on dynamic partitioning
CN106549377B (en) A kind of handover method of alternating current-direct current combined hybrid system electromechanics-electromagnetic transient hybrid simulation
CN108649554A (en) A kind of power network topology analysis method based on Gaussian elimination algorithm
CN107482633A (en) A kind of non-iterative Interval Power Flow algorithm suitable for radial distribution networks
CN115563921A (en) Method and system for determining fault transient electric quantity of flexible direct current transmission line
CN103729502B (en) A kind of method improving electromagnetic transient in power system simulation velocity
CN112684290B (en) Efficient calculation method for interelectrode short-circuit fault current of flexible direct-current power grid
CN107968399A (en) A kind of method of fast search Static Voltage Stability Region Boundary
CN105375475B (en) A kind of power system fast transient stability simulation method based on Real-time Dynamic Equivalence
CN108521344B (en) Method for evaluating reliability of two ends of command control network based on edge expansion diagram
CN109635374B (en) Short-circuit current calculation method and device based on dynamic phasor time domain method
CN102142682B (en) Method for calculating sensitivity of branch breakage based on direct-current flow model
Hou et al. Parameter estimation method of distribution network based on PMU measurement data
CN105305392A (en) Symmetrical component method for short circuit calculation of voltage-controlled type IIDG included power distribution network
CN110427635A (en) LCC-HVDC optimizes electromagnetical transient emulation method
CN108551178A (en) A kind of electromechanical transient simulation model of flexible direct current power transmission system

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190426