CN109861214A - Judge method, the system of regional power grid transient rotor angle stability weakness route - Google Patents

Judge method, the system of regional power grid transient rotor angle stability weakness route Download PDF

Info

Publication number
CN109861214A
CN109861214A CN201910152065.0A CN201910152065A CN109861214A CN 109861214 A CN109861214 A CN 109861214A CN 201910152065 A CN201910152065 A CN 201910152065A CN 109861214 A CN109861214 A CN 109861214A
Authority
CN
China
Prior art keywords
route
power grid
matrix
regional power
branch
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
Application number
CN201910152065.0A
Other languages
Chinese (zh)
Other versions
CN109861214B (en
Inventor
孙玉娇
周勤勇
张健
代倩
陈典
刘新元
施浩波
贺海磊
张立波
赵珊珊
吴俊玲
覃琴
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Shanxi Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Shanxi Electric Power 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 State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, State Grid Shanxi Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201910152065.0A priority Critical patent/CN109861214B/en
Publication of CN109861214A publication Critical patent/CN109861214A/en
Priority to PCT/CN2019/117585 priority patent/WO2020173132A1/en
Application granted granted Critical
Publication of CN109861214B publication Critical patent/CN109861214B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks

Abstract

A kind of method for judging regional power grid transient rotor angle stability weakness route of disclosure of the invention, system, utilize the initial trend of power grid, pass through simple acquisition of information and matrix operation, route fragility factor LFF matrix after different types of faults occurs for the different routes that can calculate certain regional power grid, and most weak route and most serious failure under the transient rotor angle stability of system can be directly judged by fragile factor matrix, method is simply direct, the influence of weak link and All other routes to weak link to system is very clear, simulation time can be greatlyd save, it can be used in a variety of analyses of electric system.

Description

Judge method, the system of regional power grid transient rotor angle stability weakness route
Technical field
The invention belongs to technical field of power systems more particularly to a kind of judgement regional power grid transient rotor angle stability line of weakness Method, the system on road.
Background technique
When finding electrical power system transient angle stability weakness route at present, often more concerned be system jam after Influence to system, the method that analysis failure influences mainly use emulation mode, expend the time, existing judgement weakness section or The analysis or method that the method or shortage failure of weak route influence are complex.
Summary of the invention
One of the technical issues of the present invention is directed to solve in the related technology to a certain extent.
For this purpose, the first purpose of this invention is to propose the side for judging regional power grid transient rotor angle stability weakness route Method.This method can directly acquire weak route according only to a certain status information of electric system and by tight to the weakness line influence The system failure of weight degree sequence, reduces and calculates time and simulation time, has to the transient rotor angle stability analysis of electric system Significance.
Second object of the present invention is the system for proposing to judge regional power grid transient rotor angle stability weakness route.
In order to achieve the above objectives, the judgement regional power grid transient rotor angle stability that first aspect present invention embodiment proposes is weak The method of route, comprising:
S1: the bus information for obtaining each line information of the regional power grid and being connected with each route of the regional power grid.
S2: it is connected according to each line information of the regional power grid obtained in S1 and with each route of the regional power grid Bus information calculates the power of each route of the regional power grid before failure and calculates the event of each route of the regional power grid The matrix B of each branch admittance composition before barrierL, before failure the node-branch incident matrix A of each route of the regional power grid therefore Hinder front nodal point admittance matrix B0
S3: increasing associated with the connection bus outside regional power grid branch, and by the corresponding admittance of the branch Value is set as infinitely great, and the matrix B of the corresponding each branch admittance composition modified in the S2L, the regional power grid before failure The node-branch incident matrix A of each route, the node admittance matrix B of reactance is only considered before failure0
S4: according to the different line faults of setting, branch admittance matrix, the node-branch in the S3 are accordingly modified Incidence matrix, node admittance matrix, the admittance matrix B ' of the branch after forming failureL, node-branch incident matrix A ', Node admittance matrix B '0
S5: according to the admittance matrix B ' of branch described in the S4L, node-branch incident matrix A ', node admittance matrix B′0Calculate the power of the faulty line, the fragile factor of the faulty line and the faulty line fragility factor matrix LEF。
S6: the faulty line fragility factor matrix absolute value LEF ' in the S5 is calculated, and descending arrangement is described crisp Weak factor matrix absolute value LEF ', the corresponding failure (the corresponding row of LEF ' maximum value) of maximum value is in the LEF ' sequence Most serious failure, the corresponding route (the corresponding column of LEF ' maximum value) are most weak route.
According to one embodiment of present invention, each line information of regional power grid described in the S1 includes: the route two The reactance value of the bus title at end, the route both end voltage grade and the route;With each route phase of the regional power grid Bus information even, including the bus title being connected with each route of the regional power grid, described each with the regional power grid The connected voltage class of route, the busbar voltage being connected with each route of the regional power grid, the described and regional power grid The connected bus phase angle of each route.
According to one embodiment of present invention, the matrix B of each branch admittance composition described in the S2L, node admittance square Battle array B0Acquisition only consider reactance factor.
The system for judging regional power grid transient rotor angle stability weakness route that second aspect of the present invention embodiment proposes, it is described System includes:
Module is obtained, for obtaining each line information of the regional power grid and being connected with each route of the regional power grid Bus information.
First computing module, for according to it is described first obtain module obtain each line information of the regional power grid and The bus information being connected with each route of the regional power grid, calculate failure before the regional power grid each route power and meter Calculate the matrix B of each branch admittance composition before the failure of each route of the regional power gridL, before failure the regional power grid each line Node-branch incident matrix A, the failure front nodal point admittance matrix B on road0
First modified module, increases associated with the connection bus outside regional power grid branch, and by the branch Corresponding admittance value is set as infinitely great, and the matrix B of the corresponding each branch admittance composition modified in the S2L, institute before failure State the node admittance matrix B that reactance is only considered before the node-branch incident matrix A of each route of regional power grid, failure0
Second modified module accordingly modifies second computing module for the different line faults according to setting Branch admittance matrix, node-branch incident matrix, node admittance matrix in described, the branch after forming failure are led Receive matrix B 'L, node-branch incident matrix A ', node admittance matrix B '0
Second computing module, for the admittance matrix B ' according to the branch of second modified moduleL, node- Road incidence matrix A ', node admittance matrix B '0Calculate the power of the faulty line, the fragile factor of the faulty line and institute State faulty line fragility factor matrix LEF.
Third computing module, it is absolute for calculating the faulty line fragility factor matrix in second computing module Value LEF ', and descending arranges the fragile factor matrix absolute value LEF ', the corresponding event of maximum value in the LEF ' sequence Hindering (the corresponding row of LEF ' maximum value) is most serious failure, and the corresponding route (the corresponding column of LEF ' maximum value) is most weak Route.
According to one embodiment of present invention, each line information of the regional power grid for obtaining module includes: described The reactance value of the bus title at route both ends, the route both end voltage grade and the route;It is each with the regional power grid The connected bus information of route, including the bus title being connected with each route of the regional power grid, the described and region The connected voltage class of each route of power grid, the busbar voltage being connected with each route of the regional power grid, the described and area The connected bus phase angle of each route of domain power grid.
According to one embodiment of present invention, the matrix B of each branch admittance composition described in first computing moduleL、 Node admittance matrix B0Acquisition only consider reactance factor.
The present invention can calculate certain region by simple acquisition of information and matrix operation using the initial trend of power grid Route fragility factor LFF matrix after different types of faults occurs for the different routes of power grid, and can be straight by fragile factor matrix Most weak route and the most serious failure under the transient rotor angle stability of judgement system are connect, method is simply direct, to the weakness of system The influence of link and All other routes to weak link is very clear, can greatly save simulation time, can be used for electric system In a variety of analyses.
Detailed description of the invention
By reference to the following drawings, exemplary embodiments of the present invention can be more fully understood by:
Fig. 1 is the method for disclosing the judgement regional power grid transient rotor angle stability weakness route of embodiment according to the present invention Flow chart;
Fig. 2 (a) is system maximum generation machine generator rotor angle difference simulation curve under the different faults for disclosing embodiment according to the present invention Figure;
Fig. 2 (b) is system maximum generation machine generator rotor angle difference simulation curve under the different faults for disclosing embodiment according to the present invention Figure;
Fig. 3 is the system of the disclosed judgement regional power grid transient rotor angle stability weakness route according to embodiments of the present invention Block diagram.
Specific embodiment
In order to enable those skilled in the art to better understand the solution of the present invention, below in conjunction in the embodiment of the present invention Attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is only The embodiment of a part of the invention, instead of all the embodiments.Based on the embodiments of the present invention, ordinary skill people The model that the present invention protects all should belong in member's every other embodiment obtained without making creative work It encloses.
It should be noted that description and claims of this specification and term " first " in above-mentioned attached drawing, " Two " etc. be to be used to distinguish similar objects, without being used to describe a particular order or precedence order.It should be understood that using in this way Data be interchangeable under appropriate circumstances, so as to the embodiment of the present invention described herein can in addition to illustrating herein or Sequence other than those of description is implemented.In addition, term " includes " and " having " and their any deformation, it is intended that cover Cover it is non-exclusive include, for example, the process, method, system, product or equipment for containing a series of steps or units are not necessarily limited to Step or unit those of is clearly listed, but may include be not clearly listed or for these process, methods, product Or other step or units that equipment is intrinsic.
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, right combined with specific embodiments below The present invention is described in further detail.
The present invention proposes the method for obtaining regional power grid transient rotor angle stability weakness route and failure, based on a certain of system State, computing system occur the route weakness factor matrix after different types of faults, judge power grid weakness route, and provide difference Failure sorts to the index of the weakness route influence degree, can directly give the weak link and the event of system most serious of power grid Barrier can reduce and calculate time and simulation time, have very apparent directive function to the analysis of electrical power system transient angle stability.
Under normal circumstances, to certain regional power grid, after route k breaks down, it is assumed that the power flow transfer value on route i is Δ Pki, then the ratio between the power of route i and maximum power after trend turns are as follows:
According to DC power flow calculation formula:
PSP=B0θ (2)
Order matrix A is node-branch incident matrix, then has
PL=BLATθ (4)
PSP=B0θ (5)
PL=BLATB0 -1PSP (6)
B0(AAT)-1ABL -1PL=B0(AAT)-1ABL -1BLATB0 -1PSP=PSP (7)
When route is cut-off:
Then:
In formula, PLFor the vector power of branch each before failure, BLOnly to consider what each branch admittance of reactance formed before failure Matrix, A are failure front nodal point-branch incidence matrix, B0For the node admittance matrix for only considering reactance before failure, P 'LAfter failure The power of each branch, B 'LFor only the matrix of each branch admittance composition of consideration reactance, A ' are failure posterior nodal point-branch after failure Incidence matrix, B '0For the node admittance matrix for only considering reactance after failure.
If it is assumed that each route δmaxIt is all the same, then can be | P 'lk·X′lk| as each route weakness after route k failure because Sub- index, the index is bigger, then illustrates that route is weaker.Wherein P 'lk、X′lkIt is per unit value.
When the regional power grid investigated and exterior domain power grid are there are when electrical link, before carrying out index calculating, need by Exterior domain power grid carries out virtual link by zero impedance, so that index calculating is more reasonable.
The premise of this method application thinks that power grid is substantially layering and zoning.
To the power grid of a certain level, Fig. 1 is that the judgement failure rear region stabilization of power grids of the invention for disclosing embodiment is thin The method flow diagram on weak line road, detailed process are as follows:
Step S1: each line information of the regional power grid and the bus being connected with each route of regional power grid letter are obtained Breath obtains each line information of regional power grid, including route both ends bus title, voltage class that is, by the operating status of system And the reactance value of route;The bus information being connected with route is obtained, including the mother under bus title, voltage class and operating status Line voltage, phase angle.
S2: it is connected according to each line information of the regional power grid obtained in S1 and with each route of the regional power grid Bus information calculates the power of each route of the regional power grid before failure and calculates the event of each route of the regional power grid The matrix B of each branch admittance composition before barrierL, before failure the node-branch incident matrix A of each route of the regional power grid therefore Hinder front nodal point admittance matrix B0
S3: to the connection bus outside investigated regional power grid, increase associated with the connection bus outside the regional power grid Branch, and set infinitely great for the corresponding admittance value of the branch, and the corresponding each branch admittance group modified in the S2 At matrix BL, reactance is only considered before failure before the node-branch incident matrix A of each route of the regional power grid, failure Node admittance matrix B0
S4: according to the different line faults of setting, branch admittance matrix, the node-branch in the S3 are accordingly modified Incidence matrix, node admittance matrix, the admittance matrix B ' of the branch after forming failureL, node-branch incident matrix A ', Node admittance matrix B '0
S5: according to the admittance matrix B ' of branch described in the S4L, node-branch incident matrix A ', node admittance matrix B′0And it is crisp by the power of formula (9) calculating faulty line, the fragile factor of the faulty line and the faulty line Weak factor matrix LEF.
S6: the faulty line fragility factor matrix absolute value LEF ' in the S5 is calculated, and descending arrangement is described crisp Weak factor matrix absolute value LEF ', the corresponding failure (the corresponding row of LEF ' maximum value) of maximum value is in the LEF ' sequence Most serious failure, the corresponding route (the corresponding column of LEF ' maximum value) are most weak route.
Application verification has been carried out so that power grid scheme is planned in East China Jiangsu as an example.Under the programme, East China Power Grid The 500kV interconnection of formation extra-high voltage looped network, Jiangsu Power Grid and external electrical network has 4.Under certain extreme manner, other provinces of East China are enabled Part power transmission Jiangsu, and Jiangsu Power Grid reaches and most receives electric scale greatly.The route fragility factor under different route catastrophe failures is calculated at this time Matrix selects the corresponding route of the biggish fragility factor of absolute value, and the results are shown in Table 1.
Vulnerable line and its corresponding fragile factor under 1 extreme manner of table
Vulnerable line The fragile factor
' Su rotation peak 51 --- three 51' of Anhui grassland 0.682418
' Anhui ought painting 51 --- revive Liyang 51' 0.558321
' Suzhou south 51 --- Soviet Union Liyang 51' 0.358903
' Su Anlan 51 --- Su Shuansi 51' 0.351139
' Su Anlan 51 --- Su Yudong 51' 0.329789
' Su Minzhu 51 --- Soviet Union rotation peak 51' 0.30947
Table 1 shows: the most fragile route of system is " Su rotation peak 51 --- Anhui grassland 3 51 " route, and followed by " Anhui is when painting 51 --- Soviet Union Liyang 51 ".The All other routes being followed successively by inside Jiangsu later.
It is as follows by the vulnerable line and corresponding error listing of fragile factor size sequence:
Vulnerable line and corresponding failure under 2 extreme manner of table by fragile factor sequence
In order to check the failure being affected to vulnerable line, list " Su rotation peak 51 --- Anhui grassland 3 51 ", " Anhui is when applying 51 --- the faulty line that several routes such as Soviet Union Liyang 51 " sort by the factor fragile after failure is respectively as shown in table 3 to table 5.
" Soviet Union rotation peak --- Anhui grassland three " route and fragile factor when other line faults under 3 extreme manner of table
Failure in table 3 is emulated, the results showed that, come front two " Anhui when apply 51 --- Soviet Union Liyang 51 ", " Soviet Union State south 51 --- system unstability after Soviet Union's 51 " line fault of Liyang, system is stablized after remaining line fault.Demonstrate the correct of index Property.
" Anhui when apply 51 --- Soviet Union Liyang 51 " route and fragile factor when other line faults under 4 extreme manner of table
" Suzhou south 51 --- Soviet Union Liyang 51 " route and fragile factor when other line faults under 5 extreme manner of table
Failure in table 4, table 5 is emulated, the results showed that, come the " Su rotation peak 51 of foremost --- Anhui grassland 3 51 " System unstability after line fault, system is stablized after remaining line fault.Demonstrate the correctness of index.
Shown in simulation curve Fig. 2 (a), Fig. 2 (b) of partial fault.
The system for judging regional power grid transient rotor angle stability weakness route that second aspect of the present invention embodiment proposes, Fig. 3 For the system block diagram for judging regional power grid transient rotor angle stability weakness route, system 300 includes:
Obtain module 301, for obtain each line information of the regional power grid and with each route phase of the regional power grid Bus information even.
First computing module 302, for obtaining each line information of the regional power grid that module obtains according to described first And the bus information being connected with each route of the regional power grid, calculate failure before the regional power grid each route power with And the matrix B that each branch admittance forms before the failure of each route of the calculating regional power gridL, the regional power grid before failure Node-branch incident matrix A, the failure front nodal point admittance matrix B of each route0
First modified module 303, for the connection bus outside investigated regional power grid, increase with outside the regional power grid The associated branch of connection bus, and set infinitely great for the corresponding admittance value of the branch, and corresponding modify the S2 In each branch admittance composition matrix BL, before failure the node-branch incident matrix A of each route of the regional power grid therefore The node admittance matrix B of reactance is only considered before barrier0
Second modified module 304 accordingly modifies second computing module for the different line faults according to setting It is described in branch admittance matrix, node-branch incident matrix, node admittance matrix, the branch after forming failure Admittance matrix B 'L, node-branch incident matrix A ', node admittance matrix B '0
Second computing module 305, for the admittance matrix B ' according to the branch of second modified moduleL, node- Branch incidence matrix A ', node admittance matrix B '0Calculate the power of the faulty line, the fragile factor of the faulty line and The faulty line fragility factor matrix LEF.
Third computing module 306, for calculating the faulty line fragility factor matrix in second computing module Absolute value LEF ', and descending arranges the fragile factor matrix absolute value LEF ', the corresponding institute of maximum value in the LEF ' sequence Stating failure (the corresponding row of LEF ' maximum value) is most serious failure, and the corresponding route (the corresponding column of LEF ' maximum value) is most Weak route.
According to one embodiment of present invention, each line information of the regional power grid for obtaining module 301 includes: institute State the reactance value of the bus title at route both ends, the route both end voltage grade and the route;With the regional power grid The connected bus information of each route, including the bus title being connected with each route of the regional power grid, the described and area The connected voltage class of each route of domain power grid, the busbar voltage being connected with each route of the regional power grid, it is described with it is described The connected bus phase angle of each route of regional power grid.
According to one embodiment of present invention, the matrix of each branch admittance composition described in first computing module 302 BL, node admittance matrix B0Acquisition only consider reactance factor.
Here the system for judging regional power grid transient rotor angle stability weakness route and the judgement area in first aspect embodiment The method that domain power grid transient rotor angle stability weakness route uses is same or similar, and which is not described herein again.
Above-mentioned the embodiment of the present application serial number is for illustration only, does not represent the advantages or disadvantages of the embodiments.
In above-described embodiment of the application, all emphasizes particularly on different fields to the description of each embodiment, do not have in some embodiment The part of detailed description, reference can be made to the related descriptions of other embodiments.
In several embodiments provided herein, it should be understood that disclosed technology contents can pass through others Mode is realized.Wherein, the apparatus embodiments described above are merely exemplary, such as the division of the unit or module, Only a kind of logical function partition, there may be another division manner in actual implementation, such as multiple units or module or group Part can be combined or can be integrated into another system, or some features can be ignored or not executed.Another point, it is shown Or the mutual coupling, direct-coupling or communication connection discussed can be through some interfaces, module or unit it is indirect Coupling or communication connection can be electrical or other forms.
The unit as illustrated by the separation member or module may or may not be physically separated, as The component that unit or module are shown may or may not be physical unit or module, it can and it is in one place, or Person may be distributed in multiple network units or module.Some or all of list therein can be selected according to the actual needs Member or module achieve the purpose of the solution of this embodiment.
In addition, each functional unit or module in each embodiment of the application can integrate in a processing unit or mould In block, it is also possible to each unit or module physically exists alone, can also be integrated in two or more units or module In one unit or module.Above-mentioned integrated unit or module both can take the form of hardware realization, can also use software The form of functional unit or module is realized.
If the integrated unit is realized in the form of SFU software functional unit and sells or use as independent product When, it can store in a computer readable storage medium.Based on this understanding, the technical solution of the application is substantially The all or part of the part that contributes to existing technology or the technical solution can be in the form of software products in other words It embodies, which is stored in a storage medium, including some instructions are used so that a computer Equipment (can for personal computer, server or network equipment etc.) execute each embodiment the method for the application whole or Part steps.And storage medium above-mentioned includes: that USB flash disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited Reservoir (RAM, Random Access Memory), mobile hard disk, magnetic or disk etc. be various to can store program code Medium.
The above is only the preferred embodiment of the application, it is noted that for the ordinary skill people of the art For member, under the premise of not departing from the application principle, several improvements and modifications can also be made, these improvements and modifications are also answered It is considered as the protection scope of the application.

Claims (6)

1. judging the method for regional power grid transient rotor angle stability weakness route, which is characterized in that the described method includes:
S1: the bus information for obtaining each line information of the regional power grid and being connected with each route of the regional power grid;
S2: according to each line information of the regional power grid obtained in S1 and the bus being connected with each route of the regional power grid Information, before calculating the power of each route of the regional power grid before failure and calculating the failure of each route of the regional power grid The matrix B of each branch admittance compositionL, before failure before the node-branch incident matrix A of each route of the regional power grid, failure Node admittance matrix B0
S3: increase branch associated with the connection bus outside the regional power grid, and the corresponding admittance value of the branch is set It is set to infinity, and the matrix B of the corresponding each branch admittance composition modified in the S2L, the regional power grid is each before failure The node admittance matrix B of reactance is only considered before the node-branch incident matrix A of route, failure0
S4: according to the different line faults of setting, branch admittance matrix, the node-branch association in the S3 are accordingly modified Matrix, node admittance matrix, the admittance matrix B ' of the branch after forming failureL, node-branch incident matrix A ', node Admittance matrix B '0
S5: according to the admittance matrix B ' of branch described in the S4L, node-branch incident matrix A ', node admittance matrix B '0, Calculate the power of the faulty line, the fragile factor of the faulty line and the faulty line fragility factor matrix LEF;
S6: calculating the faulty line fragility factor matrix absolute value LEF ' in the S5, and descending arrangement it is described it is fragile because Submatrix absolute value LEF ', (LEF ' is maximum for the corresponding failure in the LEF ' matrix for maximum value in the LEF ' sequence It is worth corresponding row) it is most serious failure, the corresponding route (the corresponding column of LEF ' maximum value) is most weak route.
2. the method according to claim 1, wherein each line information of regional power grid described in the S1 includes: The reactance value of the bus title at the route both ends, the route both end voltage grade and the route;With the region electricity Net the connected bus information of each route, including the bus title being connected with each route of the regional power grid, it is described with it is described The connected voltage class of each route of regional power grid, the busbar voltage being connected with each route of the regional power grid, described and institute State the connected bus phase angle of each route of regional power grid.
3. the method according to claim 1, wherein the matrix B of each branch admittance composition described in the S2L, section Point admittance matrix B0Acquisition only consider reactance factor.
4. judging the system of regional power grid transient rotor angle stability weakness route, which is characterized in that the system comprises:
Obtain module, the bus for obtaining each line information of the regional power grid and being connected with each route of the regional power grid Information;
First computing module, for according to it is described first obtain module obtain each line information of the regional power grid and with institute The connected bus information of each route of regional power grid is stated, the power of each route of the regional power grid and calculating institute before calculating failure State the matrix B of each branch admittance composition before the failure of each route of regional power gridL, each route of the regional power grid before failure Node-branch incident matrix A, failure front nodal point admittance matrix B0
First modified module, for increasing branch associated with the connection bus outside the regional power grid, and by the branch Corresponding admittance value is set as infinitely great, and the matrix B of the corresponding each branch admittance composition modified in the S2L, institute before failure State the node admittance matrix B that reactance is only considered before the node-branch incident matrix A of each route of regional power grid, failure0
Second modified module accordingly modifies the described of second computing module for the different line faults according to setting In branch admittance matrix, node-branch incident matrix, node admittance matrix, formed failure after the branch admittance square Battle array B 'L, node-branch incident matrix A ', node admittance matrix B '0
Second computing module, for the admittance matrix B ' according to the branch of second modified moduleL, the association of node-branch Matrix A ', node admittance matrix B '0Calculate the power of the faulty line, the fragile factor of the faulty line and the failure Route fragility factor matrix LEF;
Third computing module, for calculating the faulty line fragility factor matrix absolute value in second computing module LEF ', and descending arranges the fragile factor matrix absolute value LEF ', the corresponding failure of maximum value in the LEF ' sequence (the corresponding row of LEF ' maximum value) is most serious failure, and the corresponding route (the corresponding column of LEF ' maximum value) is most line of weakness Road.
5. system according to claim 4, which is characterized in that each line information of the regional power grid for obtaining module It include: the reactance value of the bus title at the route both ends, the route both end voltage grade and the route;With it is described The connected bus information of each route of regional power grid, including the bus title being connected with each route of the regional power grid, described Voltage class, the busbar voltage being connected with each route of the regional power grid, institute being connected with each route of the regional power grid State the bus phase angle being connected with each route of the regional power grid.
6. system according to claim 4, which is characterized in that each branch admittance composition described in first computing module Matrix BL, node admittance matrix B0Acquisition only consider reactance factor.
CN201910152065.0A 2019-02-28 2019-02-28 Method and system for judging weak line with stable transient power angle of regional power grid Active CN109861214B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201910152065.0A CN109861214B (en) 2019-02-28 2019-02-28 Method and system for judging weak line with stable transient power angle of regional power grid
PCT/CN2019/117585 WO2020173132A1 (en) 2019-02-28 2019-11-12 Method and system for determining transient power angle stability weak line of regional power grid

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910152065.0A CN109861214B (en) 2019-02-28 2019-02-28 Method and system for judging weak line with stable transient power angle of regional power grid

Publications (2)

Publication Number Publication Date
CN109861214A true CN109861214A (en) 2019-06-07
CN109861214B CN109861214B (en) 2022-10-04

Family

ID=66899350

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910152065.0A Active CN109861214B (en) 2019-02-28 2019-02-28 Method and system for judging weak line with stable transient power angle of regional power grid

Country Status (2)

Country Link
CN (1) CN109861214B (en)
WO (1) WO2020173132A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020173132A1 (en) * 2019-02-28 2020-09-03 中国电力科学研究院有限公司 Method and system for determining transient power angle stability weak line of regional power grid

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005160257A (en) * 2003-11-27 2005-06-16 Japan Science & Technology Agency Single phase/three-phase converter, and its control method
WO2006098703A1 (en) * 2005-02-21 2006-09-21 Domel Elektromotorji In Gospodinjski Aparati, D.D. A circuit and a method for controlling a reluctance motor
CN102403721A (en) * 2011-11-23 2012-04-04 昆明理工大学 Method for adjusting super real-time coincidence time sequence based on transient power angle stability margin
CN103094901A (en) * 2013-01-09 2013-05-08 中国电力科学研究院 Method for identification of transient state power angle instability in real time based on actual measurement response message
CN103762579A (en) * 2013-12-22 2014-04-30 南京南瑞集团公司 Direct-current power emergency control performance index calculation method for improving transient power angle stability of power grid
WO2015044803A1 (en) * 2013-09-26 2015-04-02 Abb Technology Method for connecting subsystems of an electrical power system and an intelligent electronic device therefor
US20150112498A1 (en) * 2013-09-23 2015-04-23 Sureshchandra B Patel Methods of Patel Loadflow Computation for Electrical Power System
CN105656039A (en) * 2016-03-21 2016-06-08 国网宁夏电力公司电力科学研究院 Recognition method for vulnerable lines of power system
WO2017072593A1 (en) * 2015-10-26 2017-05-04 Rakesh Goel System and method for improved performance and synchronization of a grid tie inverter
CN106684855A (en) * 2015-11-10 2017-05-17 中国电力科学研究院 Transient stability emergency control method based on key branch identification
CN106786529A (en) * 2016-12-23 2017-05-31 国网山西省电力公司电力科学研究院 A kind of distributed static security analysis method
CN107317354A (en) * 2017-06-30 2017-11-03 天津大学 A kind of multi-computer system Transient angle stability analysis method containing wind power plant
WO2018006499A1 (en) * 2016-07-06 2018-01-11 南方电网科学研究院有限责任公司 Dominant instability mode identification method and system for power system
CN107785889A (en) * 2017-10-18 2018-03-09 中国电力科学研究院有限公司 A kind of minimum start optimization method of weak receiving end power network for supporting high-capacity direct current feed-in

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10090777B2 (en) * 2011-05-08 2018-10-02 Koolbridge Solar, Inc. Inverter with independent current and voltage controlled outputs
CN103795058B (en) * 2014-02-12 2015-10-28 国家电网公司 The air extract analysis of electric power system and system failure sort method
CN107895947B (en) * 2017-12-20 2023-07-07 广东电网有限责任公司惠州供电局 Power distribution network fragile line identification method
CN109861214B (en) * 2019-02-28 2022-10-04 中国电力科学研究院有限公司 Method and system for judging weak line with stable transient power angle of regional power grid

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005160257A (en) * 2003-11-27 2005-06-16 Japan Science & Technology Agency Single phase/three-phase converter, and its control method
WO2006098703A1 (en) * 2005-02-21 2006-09-21 Domel Elektromotorji In Gospodinjski Aparati, D.D. A circuit and a method for controlling a reluctance motor
CN102403721A (en) * 2011-11-23 2012-04-04 昆明理工大学 Method for adjusting super real-time coincidence time sequence based on transient power angle stability margin
CN103094901A (en) * 2013-01-09 2013-05-08 中国电力科学研究院 Method for identification of transient state power angle instability in real time based on actual measurement response message
US20150112498A1 (en) * 2013-09-23 2015-04-23 Sureshchandra B Patel Methods of Patel Loadflow Computation for Electrical Power System
WO2015044803A1 (en) * 2013-09-26 2015-04-02 Abb Technology Method for connecting subsystems of an electrical power system and an intelligent electronic device therefor
CN103762579A (en) * 2013-12-22 2014-04-30 南京南瑞集团公司 Direct-current power emergency control performance index calculation method for improving transient power angle stability of power grid
WO2017072593A1 (en) * 2015-10-26 2017-05-04 Rakesh Goel System and method for improved performance and synchronization of a grid tie inverter
CN106684855A (en) * 2015-11-10 2017-05-17 中国电力科学研究院 Transient stability emergency control method based on key branch identification
CN105656039A (en) * 2016-03-21 2016-06-08 国网宁夏电力公司电力科学研究院 Recognition method for vulnerable lines of power system
WO2018006499A1 (en) * 2016-07-06 2018-01-11 南方电网科学研究院有限责任公司 Dominant instability mode identification method and system for power system
CN106786529A (en) * 2016-12-23 2017-05-31 国网山西省电力公司电力科学研究院 A kind of distributed static security analysis method
CN107317354A (en) * 2017-06-30 2017-11-03 天津大学 A kind of multi-computer system Transient angle stability analysis method containing wind power plant
CN107785889A (en) * 2017-10-18 2018-03-09 中国电力科学研究院有限公司 A kind of minimum start optimization method of weak receiving end power network for supporting high-capacity direct current feed-in

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
BIN LIU 等: "Recognition and Vulnerability Analysis of Key Nodes in Power Grid Based on Complex Network Centrality", 《IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—II: EXPRESS BRIEFS》 *
YUJIAO SUN 等: "Study on the transmission line vulnerability index and its application based on the transient power angle stability mechanism", 《2018 IEEE 2ND INTERNATIONAL ELECTRICAL AND ENERGY CONFERENCE (CIEEC)》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020173132A1 (en) * 2019-02-28 2020-09-03 中国电力科学研究院有限公司 Method and system for determining transient power angle stability weak line of regional power grid

Also Published As

Publication number Publication date
WO2020173132A1 (en) 2020-09-03
CN109861214B (en) 2022-10-04

Similar Documents

Publication Publication Date Title
AlHajri et al. Exploiting the radial distribution structure in developing a fast and flexible radial power flow for unbalanced three-phase networks
CN106302574B (en) A kind of service availability management method, device and its network function virtualization architecture
CN104901308B (en) Power system critical circuits discrimination method
CN102932479B (en) Virtual network mapping method for realizing topology awareness based on historical data
CN106934246A (en) The computational methods and device of power network line fragility
CN109256970B (en) MMC-MTDC transmission system monopolar grounding fault current calculation method
CN103326351A (en) Distribution type parallel load flow calculation method based on hierarchical clustering automatic partition
CN107947151A (en) A kind of recognition methods of electric system key node
CN106779269A (en) Reliability assessment method for secondary system of intelligent substation based on VC Method
CN103956735B (en) The harmonic trend of a kind of distributed generation system analyzes method
Abdelkader et al. Network usage determination using a transformer analogy
CN103823998B (en) Weak cross section determination method taking influence of network topology changes on power transmission capacity into consideration
CN105974271A (en) Distribution line fault port node impedance matrix calculation method
CN109861214A (en) Judge method, the system of regional power grid transient rotor angle stability weakness route
CN104142419B (en) Power grid short-circuit current acquiring method with influences of loads considered
CN109375034A (en) A kind of electric system branch, short-circuit current dc component acquisition methods
CN106532727A (en) Configuration method and device of installation points of dynamic reactive power compensator
CN106645868A (en) Short-circuit current scanning method and device
CN106532712A (en) Rectangular coordinate Newton method load flow calculation method for small-impedance-branch-containing power grid based on compensation method
CN111900720B (en) Transmission network fragile line identification method based on double-layer webpage sorting algorithm
Talaq Fault calculations using three terminal Thevenin’s equivalent circuit
Fan et al. Optimal fault location for power distribution systems with distributed generations using synchronized measurements
CN107045596A (en) The method and device that analysis atmospheric conditions parameter influences on the air gap discharge voltage
CN107436995A (en) Consider outer net to ground leg and the equivalent three-phase shortcircuit computational methods of sensitivity information
CN105842516B (en) A kind of computational methods of the Power System Shortcuts electric current containing a large amount of wind turbines

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