CN110718919A - GPU acceleration-based large power grid static safety analysis fault screening method - Google Patents

GPU acceleration-based large power grid static safety analysis fault screening method Download PDF

Info

Publication number
CN110718919A
CN110718919A CN201910911216.6A CN201910911216A CN110718919A CN 110718919 A CN110718919 A CN 110718919A CN 201910911216 A CN201910911216 A CN 201910911216A CN 110718919 A CN110718919 A CN 110718919A
Authority
CN
China
Prior art keywords
branch
solving
factor
generator
threads
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
CN201910911216.6A
Other languages
Chinese (zh)
Other versions
CN110718919B (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.)
Tsinghua University
State Grid Corp of China SGCC
Beijing Jiaotong University
State Grid Jiangsu Electric Power Co Ltd
Global Energy Interconnection Research Institute
Original Assignee
Tsinghua University
State Grid Corp of China SGCC
Beijing Jiaotong University
State Grid Jiangsu Electric Power Co Ltd
Global Energy Interconnection Research Institute
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 Tsinghua University, State Grid Corp of China SGCC, Beijing Jiaotong University, State Grid Jiangsu Electric Power Co Ltd, Global Energy Interconnection Research Institute filed Critical Tsinghua University
Priority to CN201910911216.6A priority Critical patent/CN110718919B/en
Publication of CN110718919A publication Critical patent/CN110718919A/en
Application granted granted Critical
Publication of CN110718919B publication Critical patent/CN110718919B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • H02H—EMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H1/00—Details of emergency protective circuit arrangements
    • H02H1/0092—Details of emergency protective circuit arrangements concerning the data processing means, e.g. expert systems, neural networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02H—EMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H7/00—Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
    • H02H7/26—Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured

Abstract

The method for screening the static safety analysis faults of the large power grid based on the GPU acceleration establishes models of branch circuit breaking and generator breaking according to the expected faults of the branch circuit breaking and the generator breaking and solves the models to obtain the workflow and data flow required by calculation, further designs a fine-grained parallel algorithm for static safety analysis fault identification based on the GPU acceleration, and obtains the hidden fault line by generating the active power of the line under each working condition in batch. The static safety analysis fault screening method provided by the invention has higher calculation precision and obtains a great acceleration effect, and meets the requirements of online application.

Description

GPU acceleration-based large power grid static safety analysis fault screening method
Technical Field
The invention relates to the field of power grid safety analysis, in particular to a method for screening static safety analysis faults of a large power grid based on GPU acceleration.
Background
With the continuous increase of the scale of the power grid, the failure forecast is more comprehensive, and higher requirements are provided for the real-time and accurate dispatching control capability of the power grid. The static safety analysis is used for ensuring the steady-state safety of the system under the condition of single element loss, and plays an indispensable important role in the real-time operation of the power system.
The conventional static safety analysis refers to the calculation of the ac power flow under each working condition, and the analysis of modern power systems [ M ] wang xi fan, beijing: scientific publishing agency, 2003:98-104, states that the conventional static security analysis has too low solving speed, the alternating current load flow solving calculation amount for all expected accidents is huge, and a Graphic Processing Unit (GPU) has strong parallel computing capability due to the superior performance of floating point operation and memory bandwidth, and is successfully applied to a plurality of scientific computing fields including the field of power systems. The core problem of load flow calculation is to solve a Large sparse linear equation set, Jalivianadi V, Zhou Z, Dinavahi V, IEEE Transactions on Parallel and distributed Systems,2012,23(7):1255 + 1266. the obtained acceleration effect is far insufficient to support the calculation cost of all expected accidents. Therefore, for static safety analysis, an expected fault set needs to be screened out, and then alternating current power flow calculation needs to be carried out. A Novel GPU-accessed Stratgy for containment screening of Static Security Analysis [ J ], Gan Zhou, Xu Zhang, Yangsheng Lang, RuiBo, Yupei Jia, Jinghuai Lin and Yanjun Feng. Fast Algorithms for sensing Analysis of Large dimensions in power systems b.sun, n.xiang, s.wang. journal of Tsinghua university.vol.28, No.1, pp.1-9,1988.
Disclosure of Invention
The embodiment of the invention provides a method for screening static safety analysis faults of a large power grid based on GPU acceleration, which is used for solving the technical problems in the prior art.
In order to achieve the purpose, the invention adopts the following technical scheme.
The method for screening the static safety analysis faults of the large power grid based on GPU acceleration comprises the following steps:
aiming at an application scene of static safety analysis of a power grid, models of branch circuit disconnection and generator disconnection are established and solved, and a workflow and a data flow are obtained;
designing and establishing a fine-grained parallel solving algorithm based on a workflow and a data flow, and generating branch circuit breaking factors and generator output power transfer distribution factors in batches;
and generating active power of each branch under each working condition by using the branch disconnection factor and the generator output power transfer distribution factor, judging whether the branch is out of range according to the active power of each branch, and screening out the fault hidden trouble branches.
Preferably, for an application scenario of static security analysis, establishing a model of branch circuit disconnection and generator disconnection and solving, and obtaining a workflow and a data stream includes the following sub-steps:
under the steady state working condition, 1/x is taken as a branch, and a susceptance matrix B without a balance node is established0Solving the susceptance matrix B0Obtaining basic data;
based on basic data, and according to branch circuits and generator information under various working conditions, M is generated in batcheslAnd ei；
Wherein x is the reactance of a branch, branch l is a line break branch, MlA node branch correlation vector of the branch l, corresponding positions at the break point of the branch are +1 and-1, and the rest are 0; e.g. of the typeiIs a unit column vector, is +1 at the corresponding position of the generator i, and is 0 for the rest.
Preferably, M is generated in batches based on basic data and according to branch and generator information of each working conditionlAnd eiThe method comprises the following steps:
creating m + s threads;
according to the thread number, indexing a branch number, and determining outflow inflow nodes i and j of the broken branch;
generating node-branch association vector M of branch l in batchl；
Wherein m is the number of branches and s is the number of generators.
Preferably, a fine-grained parallel solving algorithm is established based on the workflow and the data flow, and the step of generating branch circuit breaking factors and generator output power transfer distribution factors in batches comprises the following substeps:
solving η in batchesl＝X0MlAnd Xi＝X0eiConverting the formula into B0ηl＝MlAnd B0Xl＝eiCarrying out LU decomposition once and mass pre-generation and back-generation;
generating self-impedance between node pairs corresponding to each branch and disconnected branch node pairs in batch
Figure BDA0002214758390000031
And mutual impedance
Figure BDA0002214758390000032
Batch generation of branch disconnection factor Dk-lAnd generator output power transfer distribution factor Gk-i；
Wherein, X0Is B0The inverse matrix of (c).
Preferably, the self-impedance between the node pair corresponding to each branch and the broken branch node pair is generated in batch
Figure BDA0002214758390000033
And mutual impedance
Figure BDA0002214758390000034
Further comprising:
creating m (m + s) threads, and solving X of each branch circuit disconnection by the first m (m) threadsl-lAnd Xk-lAnd then m multiplied by s threads are used for solving X under various working conditions of power change of the generatork-i(ii) a Wherein m is the number of branches and s is the number of generators.
Preferably, the branch breaking factor D is generated in batchesk-lAnd generator output power transfer distribution factor Gk-iFurther comprising:
branch cut-off factor Dk-lBy the mathematical formula
Figure BDA0002214758390000035
Generating;
generator output power transfer distribution factor Gk-iBy the mathematical formula
Figure BDA0002214758390000036
Generating;
wherein x iskIs the reactance of line k, xlIs the reactance of the line l, τkAnd τlAnd respectively representing the standard transformation ratio of the transformer of the branch k and the standard transformation ratio of the transformer of the branch l, and when the branch does not contain the transformer, the value of the transformation ratio is 1.
Preferably, the branch disconnection factors D are generated simultaneouslyk-lAnd generator output power transfer distribution factor Gk-iFurther comprising:
creating mx (m + s) threads, solving branch disconnection factors in parallel by the first mxm threads, and solving generator output power transfer distribution factors in parallel by the last mxs threads; wherein m is the number of branches and s is the number of generators.
Preferably, the branch active power under each working condition is generated in batch by using the branch disconnection factor and the generator output power transfer distribution factor, and identifying the hidden fault line according to the branch active power further includes:
according to the formula
Figure BDA0002214758390000037
Solving the power flow of each branch; wherein
Figure BDA0002214758390000038
Active power, P, of branch k after opening branch lkActive power of ground-state line k, PlActive power for the ground state line l;
and judging whether the boundary is crossed according to the active power of each branch obtained by solving, and screening out the fault hidden trouble branches.
Preferably, the method for generating the line active power under each working condition by using the branch circuit disconnection factor and the generator output power transfer distribution factor further includes:
creating mx (m + s) threads, and solving the line active power under each broken line scene by the first mxm threads; wherein m is the number of branches and s is the number of generators.
According to the technical scheme provided by the embodiment of the invention, the method for screening the static safety analysis faults of the large power grid based on the GPU acceleration establishes the models of the branch circuit disconnection and the generator disconnection and solves the models to obtain the workflow and the data stream aiming at the expected faults of the branch circuit disconnection and the generator disconnection, and further designs the fine-grained parallel algorithm for static safety analysis fault identification based on the GPU acceleration, so that the calculation precision is high, the acceleration effect is great, and the requirements of online application are met.
Additional aspects and advantages of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
FIG. 1 is a processing flow chart of a method for static safety analysis fault screening of a large power grid based on GPU acceleration according to the present invention;
FIG. 2 is a processing flow diagram of a first batch solution embodiment of the method for static safety analysis fault screening of a large power grid based on GPU acceleration according to the present invention;
FIG. 3 is a flowchart illustrating a first batch solving embodiment of the method for static security analysis fault screening of a large power grid based on GPU acceleration according to the present invention;
fig. 4 is a display diagram of an efficiency test of the method for static safety analysis fault screening of a large power grid based on GPU acceleration provided by the present invention.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the accompanying drawings are illustrative only for the purpose of explaining the present invention, and are not to be construed as limiting the present invention.
As used herein, the singular forms "a", "an", "the" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof. It will be understood that when an element is referred to as being "connected" or "coupled" to another element, it can be directly connected or coupled to the other element or intervening elements may also be present. Further, "connected" or "coupled" as used herein may include wirelessly connected or coupled. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
It will be understood by those skilled in the art that, unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. It will be further understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the prior art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein.
For the convenience of understanding the embodiments of the present invention, the following description will be further explained by taking several specific embodiments as examples in conjunction with the drawings, and the embodiments are not to be construed as limiting the embodiments of the present invention.
Referring to fig. 1, the method for screening faults of static safety analysis of large power grid based on GPU acceleration provided by the present invention is used for static safety analysis of large power grid,
the method comprises the following steps:
aiming at an application scene of static safety analysis of a power grid, models of branch circuit disconnection and generator disconnection are established and solved, and a workflow and a data flow required by calculation are obtained;
designing a fine-grained parallel solving algorithm based on the workflow and the data flow, and generating branch circuit breaking factors and generator output power transfer distribution factors in batches;
and generating branch active power under each working condition by using the branch disconnection factor and the generator output power transfer distribution factor, judging whether each branch has power out-of-range according to the branch active power, and screening out a fault hidden line.
The method for screening the static safety analysis faults of the large power grid based on the GPU acceleration establishes models of the branch circuit disconnection and the generator disconnection and solves the models to obtain a workflow and a data stream aiming at the expected faults of the branch circuit disconnection and the generator disconnection, further designs a fine-grained parallel algorithm for static safety analysis fault identification based on the GPU acceleration, and meets the precision requirement of fault screening, wherein the calculation result is consistent with a direct current power flow method.
Further, in some preferred embodiments, the application scenario for static safety analysis may be a preset fault scenario, such as an anticipated fault for branch circuit disconnection and generator disconnection; the aforementioned first step comprises the following substeps:
under the steady state working condition, 1/x is taken as a branch, and a susceptance matrix B without a balance node is established0Solving the susceptance matrix B0Obtaining basic data;
based on basic data, and according to branch circuits and generator information under various working conditions, M is generated in batcheslAnd ei；
Wherein x is the reactance of a branch, branch l is a line break branch, MlA node branch correlation vector of the branch l, corresponding positions at the break point of the branch are +1 and-1, and the rest are 0; e.g. of the typeiIs a unit column vector, is +1 at the corresponding position of the generator i, and the rest is 0;
the first substep described above is preferably performed in the CPU; further in the second substep, there is also a process of first executing the CPU to transmit the base data to the GPU.
Further, in some preferred embodiments, the second step further includes the following sub-steps:
solving η in batchesl＝X0MlAnd Xi＝X0eiConverting the formula into B0ηl＝MlAnd B0Xl＝eiCarrying out LU decomposition once and batch forward and backward substitution, wherein the function of the LU decomposition is to avoid huge calculation amount cost paid by inversion operation;
generating self-impedance between node pairs corresponding to each branch and disconnected branch node pairs in batch
Figure BDA0002214758390000061
And mutual impedance
Figure BDA0002214758390000062
Batch generation of branch disconnection factor Dk-lAnd generator output power transfer distribution factor Gk-i；
Wherein, X0Is B0The inverse matrix of (d); batch generation of branch breaking factors Dk-lBatch-generated generator output power transfer distribution factor G for describing active power change on branch k caused by opening of branch lk-iThe method is used for describing the power flow change of a branch k caused by the active output power change of a generator i;
the three sub-step execution sequences described above are preferably performed simultaneously.
Furthermore, in order to reduce the branch of warp during calculation, the case that the generator exists in the line and the case that the generator does not exist in the line are processed simultaneously, and the third substep further comprises:
branch cut-off factor Dk-lBy the mathematical formula
Figure BDA0002214758390000063
Generating;
generator output power transfer distribution factor Gk-iBy the mathematical formula
Figure BDA0002214758390000064
Generating;
wherein x iskIs the reactance of line k, xlIs the reactance of the line l, τkAnd τlThe standard transformer ratios of line k and line l are shown, respectively, and when the branch does not contain a transformer, the value of the ratio is 1.
Further, in some preferred embodiments, the third step further includes:
according to the formulaSolving the power flow of each branch; wherein
Figure BDA0002214758390000072
Active power, P, of branch k after opening branch lkActive power of ground-state line k, PlActive power for the ground state line l;
obtaining the active power and branch power limit of each branch under each working condition according to the formula
Figure BDA0002214758390000073
And comparing the branch power with the branch limit value, judging whether each branch has power out-of-limit or not, and screening out the branch with fault hidden trouble.
The fault identification method provided by the invention can be stored in a medium in the form of a program to realize corresponding functions, and the invention provides an embodiment for exemplary illustration, which comprises the following steps:
s1 is branched by 1/x under the steady state working condition to establish a susceptance matrix B without a balance node0Solving the susceptance matrix B0Obtaining basic data; this step is done in the CPU;
s2, transmitting the basic data to the GPU through the CPU, and further executing the following steps according to the branch circuit and the generator information of each working condition:
creating m + s threads;
according to the thread number, indexing a branch number, and determining outflow inflow nodes i and j of the broken branch;
generating node-branch association vector M of branch l in batchlBatch Generation of ei；
Pseudo code for implementing the function of step S2 is shown in table 1;
TABLE 1 batch Generation node-Branch association vector pseudocode S3 batch solving ηl＝X0MlAnd Xi＝X0eiConverting the formula into B0ηl＝MlAnd B0Xl＝eiCarrying out LU decomposition once and mass pre-generation and back-generation;
creating m (m + s) threads, and generating self-impedance between the node pair corresponding to each branch and the broken branch node pair in batches by the first m (m) threads
Figure BDA0002214758390000075
And mutual impedance
Figure BDA0002214758390000076
Solving the mutual impedance of the power change of the generator under various working conditions by using the m multiplied by s threads
Figure BDA0002214758390000077
By breaking branch 1 and branch 2, solve for Xl-lAnd Xk-lFor example, as shown in FIG. 2;
creating m (m + s) threads, the first m (m) threads being defined by the mathematical expression
Figure BDA0002214758390000081
Generating a branch breaking factor Dk-lThe last m × s threads are expressed by the following mathematical formula
Figure BDA0002214758390000082
Generating a motor output power transfer profile factor Gk-i(ii) a As shown in appendix FIG. 3, each m threads process a scene, and the denominator is the same in the same scene, so thatAnd multiplexing is carried out, only one time of calculation is needed, and the numerator needs to index to the corresponding line number to obtain the reactance of the line and the transformer transformation ratio to carry out corresponding calculation. When solving the distribution factor, X needs to be readl-l、Xk-lAnd Xk-iPseudo codes for realizing the functions are shown in an appendix table 2, and since the read memories are continuous, the merged access of the memories under the same warp is ensured, so that the algorithm has high efficiency;
Figure BDA0002214758390000083
TABLE 2 solving the distribution factor in batches pseudocode S4 creates m (m + S) threads, the first m threads according to the formula
Figure BDA0002214758390000084
Solving the trend of each branch, and further obtaining the active power of the branch k after the branch l is disconnected
Figure BDA0002214758390000085
Active power P of ground state line kkActive power P of ground state line ll(ii) a Taking a broken line scenario as an example, the pseudo code for implementing the above function is shown in appendix table 3.
Figure BDA0002214758390000086
Figure BDA0002214758390000091
TABLE 3 batch Generation of line active Power pseudo-codes under various line breaking scenarios
The present embodiment also provides an example test:
under the test platform shown in table 4, the actual examples in the IEEE standard examples and matpower were tested;
Figure BDA0002214758390000092
TABLE 4 test platform
As shown in table 2, the cases 118 and 300 are IEEE standard examples, and the cases 1354, 2869, 9241, and 13659 are all from matpower, where the cases 9241 and 13659 are from real examples of european power grid.
Figure BDA0002214758390000093
TABLE 5 correctness test
As shown in fig. 4, for the 13659 node system, the computation time for batch computation of the LODF and the GSDF by the GPU only needs about 420ms, while the computation time for batch computation of the LODF and the GSDF by the CPU needs about 29900ms, which reaches about 70 times of the speed-up ratio, and the efficiency improvement is very significant.
It should be understood by those skilled in the art that the above-described fault identification method is exemplary only, and other existing or future pseudo code application types for implementing the above-described fault identification function may be applied to the embodiments of the present invention, and are included in the scope of the present invention and are incorporated herein by reference.
In summary, the method for screening the static safety analysis fault of the large power grid based on the GPU acceleration provided by the invention has the following advantages:
(1) the defect that the existing static safety analysis fault recognition time is long is overcome;
(2) the calculation result is consistent with a direct current power flow method, and the accuracy requirement of fault screening is met;
(3) can be adapted to fault analysis software;
(4) the requirement of online application is met.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (9)

1. The method for screening the static safety analysis faults of the large power grid based on GPU acceleration is characterized by comprising the following steps:
aiming at an application scene of static safety analysis of a power grid, models of branch circuit disconnection and generator disconnection are established and solved, and a workflow and a data flow are obtained;
designing and establishing a fine-grained parallel solving algorithm based on the workflow and the data flow, and generating branch circuit breaking factors and generator output power transfer distribution factors in batches;
and generating the active power of each branch under each working condition by using the branch disconnection factor and the generator output power transfer distribution factor, and judging whether the branch is out of range according to the active power of each branch so as to screen out the fault hidden trouble branches.
2. The method according to claim 1, wherein the modeling and solving for branch circuit disconnection and generator disconnection for the application scenario of static security analysis comprises the following sub-steps:
under the steady state working condition, 1/x is taken as a branch, and a susceptance matrix B without a balance node is established0Solving the susceptance matrix B0Obtaining basic data;
based on basic data, and according to branch circuits and generator information under various working conditions, M is generated in batcheslAnd ei；
Wherein x is the reactance of a branch, branch l is a line break branch, MlA node branch correlation vector of the branch l, corresponding positions at the break point of the branch are +1 and-1, and the rest are 0; e.g. of the typeiIs a unit column vector, is +1 at the corresponding position of the generator i, and is 0 for the rest.
3. The method of claim 2, wherein M is generated in batches based on the basic data and based on the branch and generator information of each operating conditionlAnd eiThe method comprises the following steps:
creating m + s threads;
according to the thread number, indexing a branch number, and determining outflow inflow nodes i and j of the broken branch;
generating node-branch association vector M of branch l in batchl；
Wherein m is the number of branches and s is the number of generators.
4. The method of claim 1, wherein the establishing a fine-grained parallel solution algorithm based on workflow and dataflow, and the batch generation of branch cutoff factors and generator output power transfer distribution factors comprises the sub-steps of:
solving η in batchesl＝X0MlAnd Xi＝X0eiConverting the formula into B0ηl＝MlAnd B0Xl＝eiCarrying out LU decomposition once and mass pre-generation and back-generation;
generating self-impedance between node pairs corresponding to each branch and disconnected branch node pairs in batchAnd mutual impedance
Figure FDA0002214758380000022
Batch generation of branch disconnection factor Dk-lAnd generator output power transfer distribution factor Gk-i；
Wherein, X0Is B0The inverse matrix of (c).
5. The method of claim 4, wherein the step of generating the self-impedance between the pair of nodes for each branch and the pair of disconnected branch nodes is performed in a batch manner
Figure FDA0002214758380000023
And mutual impedance
Figure FDA0002214758380000024
Further comprising:
creating m (m + s) threads, the first m linesSolving X of each branch circuit breakl-lAnd Xk-lAnd then m multiplied by s threads are used for solving X under various working conditions of power change of the generatork-i(ii) a Wherein m is the number of branches and s is the number of generators.
6. The method of claim 4, wherein the batch-generated branch cutoff factor D isk-lAnd generator output power transfer distribution factor Gk-iFurther comprising:
the branch breaking factor Dk-lBy the mathematical formulaGenerating;
the generator output power transfer distribution factor Gk-iBy the mathematical formula
Figure FDA0002214758380000026
Generating;
wherein x iskIs the reactance of line k, xlIs the reactance of the line l, τkAnd τlAnd respectively representing the standard transformation ratio of the transformer of the branch k and the standard transformation ratio of the transformer of the branch l, and when the branch does not contain the transformer, the value of the transformation ratio is 1.
7. Method according to claim 6, characterized in that said simultaneous generation of the branch breaking factor Dk-lAnd generator output power transfer distribution factor Gk-iFurther comprising:
creating m (m + s) threads, solving the branch circuit breaking factors in parallel by the first m (m) threads, and solving the generator output power transfer distribution factors in parallel by the last m(s) threads; wherein m is the number of branches and s is the number of generators.
8. The method according to claim 1, wherein the branch active power under each operating condition is generated in batch by using the branch disconnection factor and the generator output power transfer distribution factor, and identifying the hidden fault line according to the branch active power further comprises:
according to the formula
Figure FDA0002214758380000027
Solving the power flow of each branch; wherein
Figure FDA0002214758380000028
Active power, P, of branch k after opening branch lkActive power of ground-state line k, PlActive power for the ground state line l;
and judging whether the boundary is crossed according to the active power of each branch obtained by solving, and screening out the fault hidden trouble branches.
9. The method for identifying a power grid fault according to claim 8, wherein the generating the line active power under each working condition by using the branch circuit breaking factor and the generator output power transfer distribution factor further comprises:
creating mx (m + s) threads, and solving the line active power under each broken line scene by the first mxm threads; wherein m is the number of branches and s is the number of generators.
CN201910911216.6A 2019-09-25 2019-09-25 GPU acceleration-based large power grid static safety analysis fault screening method Active CN110718919B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910911216.6A CN110718919B (en) 2019-09-25 2019-09-25 GPU acceleration-based large power grid static safety analysis fault screening method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910911216.6A CN110718919B (en) 2019-09-25 2019-09-25 GPU acceleration-based large power grid static safety analysis fault screening method

Publications (2)

Publication Number Publication Date
CN110718919A true CN110718919A (en) 2020-01-21
CN110718919B CN110718919B (en) 2021-06-01

Family

ID=69210838

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910911216.6A Active CN110718919B (en) 2019-09-25 2019-09-25 GPU acceleration-based large power grid static safety analysis fault screening method

Country Status (1)

Country Link
CN (1) CN110718919B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105955712A (en) * 2016-04-26 2016-09-21 东南大学 Direct current fault screening method based on GPU acceleration
CN106026107A (en) * 2016-07-26 2016-10-12 东南大学 QR decomposition method of power flow Jacobian matrix for GPU acceleration
CN106354479A (en) * 2016-08-12 2017-01-25 东南大学 GPU acceleration QR decomposition method for a large number of isomorphic sparse matrixes
CN106874113A (en) * 2017-01-19 2017-06-20 国电南瑞科技股份有限公司 A kind of many GPU heterogeneous schemas static security analysis computational methods of CPU+

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105955712A (en) * 2016-04-26 2016-09-21 东南大学 Direct current fault screening method based on GPU acceleration
CN106026107A (en) * 2016-07-26 2016-10-12 东南大学 QR decomposition method of power flow Jacobian matrix for GPU acceleration
CN106354479A (en) * 2016-08-12 2017-01-25 东南大学 GPU acceleration QR decomposition method for a large number of isomorphic sparse matrixes
CN106874113A (en) * 2017-01-19 2017-06-20 国电南瑞科技股份有限公司 A kind of many GPU heterogeneous schemas static security analysis computational methods of CPU+
WO2018133348A1 (en) * 2017-01-19 2018-07-26 国电南瑞科技股份有限公司 Static security analysis computation method, apparatus, and computer storage medium

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
NORBERTO GARCIA等: "Parallel power flow solutions using a biconjugate gradient algorithm and a Newton method: A GPU-based approach", 《IEEE PES GENERAL MEETING》 *
张伯明等: "《高等电力网络分析》", 31 August 1996 *
温柏坚等: "匹配GPU编程架构的直流故障筛选算法设计", 《电力自动化设备》 *
龚正毅编著: "《线性电路理论》", 31 December 1987 *

Also Published As

Publication number Publication date
CN110718919B (en) 2021-06-01

Similar Documents

Publication Publication Date Title
CN101957745B (en) Method for converting ladder diagram language into structure text language
CN107862103B (en) Simulation model generation method and system based on ADPSS
CN115688488B (en) Criterion determination method, criterion determination device, criterion determination equipment, computer storage medium and program product
Lei et al. Developing a benchmark test system for electric power grid cyber-physical reliability studies
CN116595395B (en) Inverter output current prediction method and system based on deep learning
CN110718919B (en) GPU acceleration-based large power grid static safety analysis fault screening method
CN111695238A (en) System and method for visualizing transient simulation data of power system
CN112345996B (en) Energy controller loop event generation method based on edge calculation
CN116382214A (en) Script automatic triggering flow early warning method and device and computer equipment
CN112800637B (en) Intelligent substation simulation test device and method
CN109861214B (en) Method and system for judging weak line with stable transient power angle of regional power grid
Zhang et al. Batch computing method for sensitivity analysis of large power grids based on GPU acceleration
CN110391935B (en) Measurement degree evaluation method and system considering information physical coupling characteristics and information disturbance
CN110148065B (en) Short-circuit current calculation processing method and system considering flexible and direct influence
CN112803428A (en) Receiving-end main network frame dynamic reactive power supply configuration node selection method and terminal
CN112906169A (en) Computer medium based on virtual stream fault positioning and application thereof
CN110783897A (en) GPU acceleration-based power grid short-circuit current parallel computing method and system
Doherr et al. Fieldbus material take-off estimation: Towards an automated cost estimation of fieldbus installations
CN116109108B (en) Distribution network automatic switch position distribution generation method, device and storage medium
CN114421439B (en) Power grid protection fixed value data changing method based on visualization and expert system
Parekh et al. Reliability Evaluation in Physical Power Systems and Communication Network for Protection Failures
CN117290688A (en) Virtual loop verification method and device, electronic equipment and storage medium
Yang et al. A unified scheme for power system transient stability simulation and transient stability assessment
CN113191579A (en) Risk identification and evaluation method and system for 500 kV 2/3 power grid
CN116973668A (en) Power grid line fault diagnosis method and device, electronic equipment and storage medium

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