Preconceived fault scanning method, device and equipment applied to power system element

Abstract

The invention relates to an anticipated fault scanning method and device applied to a power system element, computer equipment, a storage medium and a computer program product, and relates to the technical field of power grid safety monitoring. The method can shorten the time of detailed load flow calculation and improve the efficiency of anticipated fault scanning. The method comprises the steps of obtaining a node susceptance matrix and an active power flow total increment of a to-be-detected element under an expected fault condition, and obtaining a voltage phase angle increment corresponding to each bus in the to-be-detected element according to the node susceptance matrix and the active power flow total increment; for each branch in the to-be-detected element, obtaining the active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch; and determining a weight coefficient and a power flow limit value of the current branch, and inputting the active power flow, the weight coefficient and the power flow limit value into a preset severity index operation model to obtain an anticipated fault severity index when the current branch is disconnected.

Description

Technical Field

[0001] The present application relates to the technical field of power grid security monitoring, and in particular to a method, device, computer equipment, storage medium and computer program product for scanning anticipated faults applied to power system components. Background Art

[0002] As the scale of the power grid continues to expand, the types of possible faults are also increasing, so it is necessary to classify the different forms of faults according to different conditions or criteria. In the era when there was no network analysis capability, the main safety measure of the power system was how to ensure sufficient spinning reserve, which made the emergency cost heavy.

[0003] In traditional technology, the ability of network security analysis is improved by calculating the power flow through computers, that is, the power flow calculated for a fault is a fault analysis. Early anticipated fault analysis used AC power flow and tributary power flow to continuously calculate the power flow solutions of various fault conditions in the anticipated fault set to determine its impact on the safe operation of the system. However, this approach has the defect of low fault scanning efficiency and is difficult to cope with the ever-expanding scale of power systems and the increasing number of anticipated faults. Summary of the invention

[0004] Based on this, it is necessary to provide a method, device, computer equipment, computer-readable storage medium and computer program product for anticipatory fault scanning of power system components in response to the above technical problems.

[0005] In a first aspect, the present application provides a method for scanning anticipated faults applied to power system components. The method comprises:

[0006] Obtaining a node susceptance matrix of the element to be detected under the expected fault condition, determining the total increment of active power flow of the element to be detected under the expected fault condition, and obtaining a voltage phase angle increment corresponding to each bus in the element to be detected according to the node susceptance matrix and the total increment of active power flow;

[0007] For each branch in the element to be detected, determining a first voltage phase angle increment, a second voltage phase angle increment and a reactance value of the current branch corresponding to the current branch;

[0008] Obtaining an active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch;

[0009] Obtaining the active power flow of the current branch according to the base power flow of the current branch and the active power flow increment;

[0010] The weight coefficient and power flow limit value of the current branch are determined, and the active power flow, weight coefficient and power flow limit value are input into a preset severity index calculation model to obtain the expected fault severity index when the current branch is disconnected.

[0011] In one embodiment, the step of obtaining a node susceptance matrix of the component to be detected under a predicted fault condition includes:

[0012] Determine the initial matrix corresponding to the base state power flow of the current branch; obtain the matrix of the base state power flow after the change under the expected fault condition based on the initial matrix; use the matrix correction method to correct the changed matrix to obtain the node susceptance matrix of the component to be detected under the expected fault condition.

[0013] In one embodiment, obtaining the voltage phase angle increment corresponding to each bus in the element to be detected according to the node susceptance matrix and the total increment of active power flow includes:

[0014] The node susceptance matrix and the total increment of active power flow are input into a preset voltage phase angle increment matrix operation model to obtain the voltage phase angle increment matrix of the component to be detected under the expected fault condition; based on the voltage phase angle increment matrix, the voltage phase angle increment corresponding to each bus in the component to be detected is determined.

[0015] In one embodiment, determining a first voltage phase angle increment, a second voltage phase angle increment, and a reactance value of the current branch corresponding to the current branch includes:

[0016] Determine the first bus and the second bus corresponding to the current branch; use the voltage phase angle increment of the first bus as the first voltage phase angle increment, and use the voltage phase angle increment of the second bus as the second voltage phase angle increment; based on the first bus and the second bus, obtain the reactance value of the current branch.

[0017] In one embodiment, obtaining the active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch includes:

[0018] According to the first voltage phase angle increment and the second voltage phase angle increment, a voltage phase angle increment difference corresponding to the first voltage phase angle increment and the second voltage phase angle increment is obtained; the voltage phase angle increment difference and the reactance value are input into a preset active power flow increment calculation model to obtain the active power flow increment of the current branch.

[0019] In one embodiment, the determining the total increment of active power flow of the component to be detected under the expected fault condition includes:

[0020] The vector of the voltage phase angle of each node in the element to be detected is obtained; the vector of the voltage phase angle of each node in the element to be detected is input into a preset active power flow total increment calculation model to obtain the active power flow total increment.

[0021] In a second aspect, the present application also provides a device for scanning anticipated faults applied to power system components. The device comprises:

[0022] A phase angle increment module is used to obtain the node susceptance matrix of the element to be detected under the expected fault condition, determine the total increment of active power flow of the element to be detected under the expected fault condition, and obtain the voltage phase angle increment corresponding to each bus in the element to be detected according to the node susceptance matrix and the total increment of active power flow;

[0023] A reactance calculation module, used to determine, for each branch in the element to be detected, a first voltage phase angle increment, a second voltage phase angle increment corresponding to a current branch, and a reactance value of the current branch;

[0024] A power flow increment module, used for obtaining the active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value;

[0025] An active power flow module, used to obtain the active power flow of the current branch according to the base power flow of the current branch and the active power flow increment;

[0026] The index determination module is used to determine the weight coefficient and power flow limit value of the current branch, input the active power flow, weight coefficient and power flow limit value into a preset severity index calculation model, and obtain the expected fault severity index when the current branch is disconnected.

[0027] In a third aspect, the present application further provides a computer device. The computer device includes a memory and a processor, the memory stores a computer program, and the processor implements the following steps when executing the computer program:

[0028] Obtain the node susceptance matrix of the component to be detected under the expected fault condition, determine the total increment of the active power flow of the component to be detected under the expected fault condition, and obtain the voltage phase angle increment corresponding to each bus in the component to be detected according to the node susceptance matrix and the total increment of the active power flow; for each branch in the component to be detected, determine the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch corresponding to the current branch; obtain the active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch; obtain the active power flow of the current branch according to the ground state power flow of the current branch and the active power flow increment; determine the weight coefficient and power flow limit value of the current branch, input the active power flow, weight coefficient and power flow limit value into a preset severity index calculation model, and obtain the expected fault severity index when the current branch is disconnected.

[0029] In a fourth aspect, the present application further provides a computer-readable storage medium. The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the following steps are implemented:

[0030] Obtain the node susceptance matrix of the component to be detected under the expected fault condition, determine the total increment of the active power flow of the component to be detected under the expected fault condition, and obtain the voltage phase angle increment corresponding to each bus in the component to be detected according to the node susceptance matrix and the total increment of the active power flow; for each branch in the component to be detected, determine the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch corresponding to the current branch; obtain the active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch; obtain the active power flow of the current branch according to the ground state power flow of the current branch and the active power flow increment; determine the weight coefficient and power flow limit value of the current branch, input the active power flow, weight coefficient and power flow limit value into a preset severity index calculation model, and obtain the expected fault severity index when the current branch is disconnected.

[0031] In a fifth aspect, the present application further provides a computer program product. The computer program product includes a computer program, and when the computer program is executed by a processor, the following steps are implemented:

[0032] Obtain the node susceptance matrix of the component to be detected under the expected fault condition, determine the total increment of the active power flow of the component to be detected under the expected fault condition, and obtain the voltage phase angle increment corresponding to each bus in the component to be detected according to the node susceptance matrix and the total increment of the active power flow; for each branch in the component to be detected, determine the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch corresponding to the current branch; obtain the active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch; obtain the active power flow of the current branch according to the ground state power flow of the current branch and the active power flow increment; determine the weight coefficient and power flow limit value of the current branch, input the active power flow, weight coefficient and power flow limit value into a preset severity index calculation model, and obtain the expected fault severity index when the current branch is disconnected.

[0033] The above-mentioned anticipated fault scanning method, device, computer equipment, storage medium and computer program product applied to power system components obtain the node susceptance matrix of the component to be detected under the anticipated fault condition, determine the total increment of the active power flow of the component to be detected under the anticipated fault condition, and obtain the voltage phase angle increment corresponding to each bus in the component to be detected according to the node susceptance matrix and the total increment of the active power flow; for each branch in the component to be detected, determine the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch corresponding to the current branch; obtain the active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch; obtain the active power flow of the current branch according to the base state power flow and the active power flow increment of the current branch; determine the weight coefficient and power flow limit value of the current branch, input the active power flow, weight coefficient and power flow limit value into a preset severity index calculation model, and obtain the anticipated fault severity index when the current branch is disconnected. Therefore, a fast scanning method for anticipated faults based on graph parallel computing and superposition principle is proposed, as well as a detailed power flow calculation method for anticipated faults that reuses the base state power flow bus numbering results and parallel calls of sub-queries on the graph database platform, thereby shortening the time of detailed power flow calculation and improving the efficiency of anticipated fault scanning. BRIEF DESCRIPTION OF THE DRAWINGS

[0034] FIG1 is an application environment diagram of a method for scanning anticipated faults applied to power system components in one embodiment;

[0035] Figure 2 A schematic flow chart of a method for scanning anticipated faults applied to power system components in one embodiment;

[0036] Figure 3 In one embodiment Schematic diagram of the linear network model;

[0037] Figure 4 In one embodiment, a branch circuit is opened Schematic diagram of applying the superposition principle when

[0038] Figure 5 A schematic flow chart of a method for scanning anticipated faults applied to power system components in a specific embodiment;

[0039] Figure 6 is a structural block diagram of a predicted fault scanning device applied to a power system component in one embodiment;

[0040] FIG. 7 is a diagram showing the internal structure of a computer device in one embodiment. DETAILED DESCRIPTION

[0041] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0042] The anticipated fault scanning method for power system components provided in the embodiment of the present application can be applied in the application environment shown in FIG1. ​​The terminal communicates with the server through a network. The data storage system can store data that the server needs to process. The data storage system can be integrated on the server, or placed on the cloud or other network servers.

[0043] Specifically, the anticipated fault scanning method for power system components provided in the embodiments of the present application can be executed by a server.

[0044] Exemplarily, the server obtains the node admittance matrix of the component to be detected under the expected fault condition, determines the total increment of the active power flow of the component to be detected under the expected fault condition, and obtains the voltage phase angle increment corresponding to each bus in the component to be detected based on the node admittance matrix and the total increment of the active power flow; for each branch in the component to be detected, the server determines the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch corresponding to the current branch; the server obtains the active power flow increment of the current branch based on the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch; obtains the active power flow of the current branch based on the base state power flow and the active power flow increment of the current branch; the server determines the weight coefficient and power flow limit value of the current branch, inputs the active power flow, weight coefficient and power flow limit value into a preset severity index calculation model, and obtains the expected fault severity index when the current branch is disconnected.

[0045] In the application environment shown in FIG1 , the terminal may be, but is not limited to, various personal computers, notebook computers, smart phones, and tablet computers. The server may be implemented as an independent server or a server cluster consisting of multiple servers.

[0046] In one embodiment, Figure 2 As shown, a method for scanning anticipated faults applied to power system components is provided, and the method is described by taking the server in FIG. 1 as an example, and includes the following steps:

[0047] Step S201, obtaining the node susceptance matrix of the component to be detected under the expected fault condition, determining the total increment of active power flow of the component to be detected under the expected fault condition, and obtaining the voltage phase angle increment corresponding to each bus in the component to be detected according to the node susceptance matrix and the total increment of active power flow.

[0048] Among them, the node susceptance matrix is ​​a matrix that describes the relationship between conductance and susceptance between nodes in the power system.

[0049] Among them, the active power flow increment refers to the amplitude of the change in the active power flow of other nodes in the system when the power generation or load of a node in the power system changes; ΔP is usually used to represent the active power flow increment, which indicates the change in power flow caused by the change in node power.

[0050] Specifically, the server obtains the node susceptance matrix of the component to be detected under the expected fault condition, determines the total increment of active power flow of the component to be detected under the expected fault condition, and obtains the voltage phase angle increment corresponding to each bus in the component to be detected according to the node susceptance matrix and the total increment of active power flow.

[0051] Step S202: for each branch in the element to be detected, determine the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch corresponding to the current branch.

[0052] Specifically, the server determines a first voltage phase angle increment, a second voltage phase angle increment, and a reactance value of the current branch corresponding to the current branch in the element to be detected.

[0053] Step S203, obtaining the active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch.

[0054] Specifically, the server obtains the active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch.

[0055] Step S204, obtaining the active power flow of the current branch according to the base power flow and the active power flow increment of the current branch.

[0056] Among them, active power flow refers to the transmission and distribution process of active power in the power system on the transmission line. It describes the flow of active power between the power source node and the load node.

[0057] Specifically, the server obtains the active power flow of the current branch according to the base power flow and the active power flow increment of the current branch.

[0058] For example, branch The flow of each branch after disconnection can be calculated by the following formula:

[0059]

[0060] In the above formula, For branch After disconnection The meritorious trend on is the ground state current, is the active power flow increment.

[0061] Step S205, determine the weight coefficient and power flow limit value of the current branch, input the active power flow, weight coefficient and power flow limit value into a preset severity index calculation model, and obtain the expected fault severity index when the current branch is disconnected.

[0062] Among them, the preset severity index calculation model can be:

[0063]

[0064] In the above formula, For branch The weight coefficient of is a branch set, For branch When the branch is disconnected The trend, For branch The flow limit value.

[0065] Specifically, the server determines the weight coefficient and power flow limit value of the current branch, and inputs the active power flow, weight coefficient and power flow limit value into a preset severity index calculation model to obtain the expected fault severity index when the current branch is disconnected.

[0066] In the above-mentioned anticipated fault scanning method applied to power system components, the node susceptance matrix of the component to be detected under the anticipated fault condition is obtained to determine the total increment of the active power flow of the component to be detected under the anticipated fault condition, and the voltage phase angle increment corresponding to each bus in the component to be detected is obtained according to the node susceptance matrix and the total increment of the active power flow; for each branch in the component to be detected, the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch corresponding to the current branch are determined; according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch, the active power flow increment of the current branch is obtained; according to the base state power flow and the active power flow increment of the current branch, the active power flow of the current branch is obtained; the weight coefficient and power flow limit value of the current branch are determined, and the active power flow, weight coefficient and power flow limit value are input into a preset severity index calculation model to obtain the anticipated fault severity index when the current branch is disconnected. Therefore, a fast scanning method for anticipated faults based on graph parallel computing and superposition principle is proposed, as well as a detailed power flow calculation method for anticipated faults that reuses the base state power flow bus numbering results and parallel calls of sub-queries on the graph database platform, thereby shortening the time of detailed power flow calculation and improving the efficiency of anticipated fault scanning.

[0067] It should be noted that in the analysis of anticipated faults in power systems, in order to improve the speed of fault scanning, the direct method is one of the commonly used fault scanning methods, that is, to quickly calculate the approximate power flow after the fault, usually using the DC power flow model, and its basic calculation formula is:

[0068]

[0069] in, is the node susceptance matrix of the currently running network, is the vector of voltage phase angles of each node in the network, Active power vector injected into the node.

[0070] Considering the linear network of the power system, the superposition principle can be applied, and the above The algorithm basically meets this condition. The following specifically describes the graph parallel computing process of the fast scanning method based on the superposition principle for an open circuit fault in a certain line in a power network.

[0071] For Figure 3 The network model shown, the branch When the circuit is open, it is equivalent to and Add an electric impedance between The branch, such as Figure 4 (a) as shown; branch The increase will make the bus and The net injected power changes, and the change is the power flowing through the newly added branch Since the anticipated fault analysis always starts from a converged base state power flow solution, Figure 4 The network in (b) shows the ground state power flow solution, which is known. According to the superposition principle, if we want to open the branch The network power flow solutions after , can be obtained respectively Figure 4 The network power flow solution and the base state power flow solution in (c) can be superimposed together.

[0072] against Figure 4 The network in (c) has

[0073]

[0074] in, is the node susceptance matrix after the fault, is the node voltage phase angle change matrix; according to the above formula, we first need to find and , then we can solve for .

[0075] In one embodiment, in the above step S201, obtaining the node susceptance matrix of the component to be detected under the expected fault condition specifically includes the following steps:

[0076] Determine the initial matrix corresponding to the base state flow of the current branch; obtain the changed matrix of the base state flow under the expected fault condition based on the initial matrix; use the matrix correction method to correct the changed matrix to obtain the node susceptance matrix of the component to be detected under the expected fault condition.

[0077] Among them, the base state power flow refers to the power flow distribution of the power system under steady-state operating conditions, which describes the power flow direction and power size between power sources, loads and transmission lines.

[0078] Specifically, the server determines the initial matrix corresponding to the base state flow of the current branch; based on the initial matrix, obtains the matrix after the base state flow changes under the expected fault condition; uses the matrix correction method to correct the changed matrix to obtain the node susceptance matrix of the component to be detected under the expected fault condition.

[0079] For example, the matrix correction method is a quick correction in the analysis of expected faults. A method for the matrix. When the dimension of the electrical coefficient matrix is ​​constant, the first-order correction technique can be used.

[0080] For the base state power flow, the matrix can be decomposed into:

[0081]

[0082] When the grid structure or parameters partially change, the changed matrix can be expressed as:

[0083]

[0084]

[0085] in, and They are dimensional vector, is a scalar

[0086]

[0087]

[0088]

[0089]

[0090] Therefore, the modified formula can be written as:

[0091]

[0092]

[0093]

[0094]

[0095]

[0096]

[0097] That is, the corrected matrix can be calculated based on the previous decomposed matrix, and the added corrected matrix satisfies:

[0098]

[0099]

[0100]

[0101]

[0102] The subscript 1 here represents the first fault. It can be seen that for multiple faults, the above process can be repeated until the entire matrix decomposition process is completed.

[0103] In this embodiment, the changed matrix is ​​corrected by using a matrix correction method, so as to quickly and accurately obtain the node susceptance matrix of the component to be detected under the expected fault condition.

[0104] In one embodiment, in the above step S201, the voltage phase angle increment corresponding to each bus in the element to be detected is obtained according to the node susceptance matrix and the total increment of active power flow, which specifically includes the following steps:

[0105] The node susceptance matrix and the total increment of active power flow are input into the preset voltage phase angle increment matrix operation model to obtain the voltage phase angle increment matrix of the component to be detected under the expected fault condition; according to the voltage phase angle increment matrix, the voltage phase angle increment corresponding to each bus in the component to be detected is determined.

[0106] Among them, the preset voltage phase angle increment matrix operation model can be

[0107]

[0108] In the above formula, is the node susceptance matrix after the fault, is the node voltage phase angle increment matrix, is the total increase of active power flow.

[0109] Specifically, the server inputs the node susceptance matrix and the total increment of active power flow into a preset voltage phase angle increment matrix operation model to obtain the voltage phase angle increment matrix of the component to be detected under the expected fault condition; based on the voltage phase angle increment matrix, the voltage phase angle increment corresponding to each bus in the component to be detected is determined.

[0110] In this embodiment, by using a preset voltage phase angle increment matrix operation model, the voltage phase angle increment matrix of the component to be detected under the expected fault condition is accurately obtained, thereby determining the voltage phase angle increment corresponding to each bus in the component to be detected.

[0111] In one embodiment, in the above step S202, determining the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch corresponding to the current branch specifically includes the following steps:

[0112] Determine the first bus and the second bus corresponding to the current branch; use the voltage phase angle increment of the first bus as the first voltage phase angle increment, and use the voltage phase angle increment of the second bus as the second voltage phase angle increment; based on the first bus and the second bus, obtain the reactance value of the current branch.

[0113] Among them, in the power system, busbar refers to the conductive rod or wire that connects multiple power equipment or circuits together.

[0114] Specifically, the server determines the first bus and the second bus corresponding to the current branch; uses the voltage phase angle increment of the first bus as the first voltage phase angle increment, and uses the voltage phase angle increment of the second bus as the second voltage phase angle increment; based on the first bus and the second bus, obtains the reactance value of the current branch.

[0115] In this embodiment, by first determining the first bus and the second bus corresponding to the current branch, the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch are correspondingly determined.

[0116] In one embodiment, in the above step S203, obtaining the active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch specifically includes the following steps:

[0117] According to the first voltage phase angle increment and the second voltage phase angle increment, a voltage phase angle increment difference corresponding to the first voltage phase angle increment and the second voltage phase angle increment is obtained; the voltage phase angle increment difference and the reactance value of the current branch are input into a preset active power flow increment calculation model to obtain the active power flow increment of the current branch.

[0118] The preset active power flow increment calculation model may be:

[0119]

[0120] In the above formula, For branch Active power flow increment; For bus The voltage phase angle increment; For branch reactance.

[0121] Specifically, the server obtains the voltage phase angle increment difference corresponding to the first voltage phase angle increment and the second voltage phase angle increment based on the first voltage phase angle increment and the second voltage phase angle increment; inputs the voltage phase angle increment difference and the reactance value of the current branch into a preset active power flow increment calculation model to obtain the active power flow increment of the current branch.

[0122] In this embodiment, the active power flow increment of the current branch is accurately obtained by using the preset active power flow increment calculation model.

[0123] In one embodiment, in the above step S201, determining the total increment of active power flow of the component to be detected under the expected fault condition specifically includes the following steps:

[0124] The vector of the voltage phase angle of each node in the element to be detected is obtained; the vector of the voltage phase angle of each node in the element to be detected is input into a preset active power flow total increment calculation model to obtain the active power flow total increment.

[0125] The preset total increment calculation model of active power flow can be:

[0126]

[0127] In the above formula, is the vector of the voltage phase angles of each node in the element to be detected.

[0128] Specifically, the server obtains the vector of the voltage phase angle of each node in the element to be detected; inputs the vector of the voltage phase angle of each node in the element to be detected into a preset active power flow total increment calculation model to obtain the active power flow total increment.

[0129] In this embodiment, the total increment of active power flow corresponding to the component to be detected is accurately obtained by using the preset total increment calculation model of active power flow.

[0130] In one embodiment, Figure 6 As shown, a method for scanning anticipated faults applied to power system components in a specific embodiment is provided, which specifically includes the following steps:

[0131] Step S501, determine the initial matrix corresponding to the base state power flow of the current branch; obtain the changed matrix of the base state power flow under the expected fault condition according to the initial matrix; use the matrix correction method to correct the changed matrix to obtain the node susceptance matrix of the component to be detected under the expected fault condition.

[0132] Step S502, obtaining the vector of the voltage phase angle of each node in the element to be detected; inputting the vector of the voltage phase angle of each node in the element to be detected into a preset active power flow increment calculation model to obtain the total active power flow increment of the element to be detected under the expected fault condition.

[0133] Step S503, input the node susceptance matrix and the total increment of active power flow into the preset voltage phase angle increment matrix operation model to obtain the voltage phase angle increment matrix of the component to be detected under the expected fault condition; according to the voltage phase angle increment matrix, determine the voltage phase angle increment corresponding to each bus in the component to be detected.

[0134] Step S504, for each branch in the element to be detected, determine the first bus and the second bus corresponding to the current branch; use the voltage phase angle increment of the first bus as the first voltage phase angle increment, and use the voltage phase angle increment of the second bus as the second voltage phase angle increment; based on the first bus and the second bus, obtain the reactance value of the current branch.

[0135] Step S505, obtaining a voltage phase angle increment difference corresponding to the first voltage phase angle increment and the second voltage phase angle increment based on the first voltage phase angle increment and the second voltage phase angle increment; inputting the voltage phase angle increment difference and the reactance value into a preset active power flow increment calculation model to obtain the active power flow increment of the current branch.

[0136] Step S506, obtaining the active power flow of the current branch according to the base power flow and the active power flow increment of the current branch.

[0137] Step S507, determine the weight coefficient and power flow limit value of the current branch, input the active power flow, weight coefficient and power flow limit value into the preset severity index calculation model, and obtain the expected fault severity index when the current branch is disconnected.

[0138] The beneficial effects brought by the above embodiments are as follows:

[0139] This scheme proposes a fast scanning method for anticipated faults based on graph parallel computing and superposition principles, as well as a detailed flow calculation method for anticipated faults that reuses the base state flow bus numbering results and parallel calls of graph database platform subqueries. By combining the above methods, the calculation time for fast anticipated fault scanning is increased by 13 times, and the calculation time for 100 detailed flows of a 10,000-level bus system is less than 3 seconds. Specifically: for the main grid 2749 bus system, the fast scanning calculation time for 2003 faults is less than 300 milliseconds, and the average calculation time for each fault is 0.15 milliseconds; for the main distribution network 8505 bus system, the fast scanning calculation time for 2128 faults is less than 800 milliseconds, and the average calculation time for each fault is 0.375 milliseconds; compared with the average time of 2.09 milliseconds for the current single fault fast scanning of the main grid (2749 bus), the calculation speed has increased by about 14 times.

[0140] For example, Table 1 and Table 2 respectively give the calculation performance of the main grid 2749 busbar system and the main distribution grid 8505 busbar system based on graph calculation.

[0141] Table 1 Fast scan calculation results of the expected fault analysis of the main grid 2749 bus system

[0142]

[0143] Table 2 Fast scan calculation results of the expected fault analysis of the main distribution network 8505 bus system

[0144]

[0145] It should be understood that, although the various steps in the flowcharts involved in the above-mentioned embodiments are displayed in sequence according to the indication of the arrows, these steps are not necessarily executed in sequence according to the order indicated by the arrows. Unless there is a clear explanation in this article, the execution of these steps does not have a strict order restriction, and these steps can be executed in other orders. Moreover, at least a part of the steps in the flowcharts involved in the above-mentioned embodiments can include multiple steps or multiple stages, and these steps or stages are not necessarily executed at the same time, but can be executed at different times, and the execution order of these steps or stages is not necessarily to be carried out in sequence, but can be executed in turn or alternately with other steps or at least a part of the steps or stages in other steps.

[0146] Based on the same inventive concept, the embodiment of the present application also provides a device for predicting fault scanning applied to power system components for implementing the above-mentioned method for predicting fault scanning applied to power system components. The implementation scheme for solving the problem provided by the device is similar to the implementation scheme recorded in the above-mentioned method, so the specific limitations in one or more embodiments of the device for predicting fault scanning applied to power system components provided below can refer to the limitations of the method for predicting fault scanning applied to power system components above, and will not be repeated here.

[0147] In one embodiment, Figure 6 As shown, a device for scanning anticipated faults applied to power system components is provided, comprising:

[0148] The phase angle increment module 601 is used to obtain the node susceptance matrix of the element to be detected under the expected fault condition, determine the total increment of the active power flow of the element to be detected under the expected fault condition, and obtain the voltage phase angle increment corresponding to each bus in the element to be detected according to the node susceptance matrix and the total increment of the active power flow;

[0149] The reactance calculation module 602 is used to determine, for each branch in the element to be detected, a first voltage phase angle increment, a second voltage phase angle increment and a reactance value of the current branch corresponding to the current branch;

[0150] The power flow increment module 603 is used to obtain the active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value;

[0151] Active power flow module 604, used to obtain the active power flow of the current branch according to the base power flow and active power flow increment of the current branch;

[0152] The index determination module 605 is used to determine the weight coefficient and power flow limit value of the current branch, input the active power flow, weight coefficient and power flow limit value into the preset severity index calculation model, and obtain the expected fault severity index when the current branch is disconnected.

[0153] In one embodiment, the phase angle increment module 601 is also used to determine the initial matrix corresponding to the base state flow of the current branch; based on the initial matrix, the matrix after the base state flow changes under the expected fault condition is obtained; and the changed matrix is ​​corrected using the matrix correction method to obtain the node susceptance matrix of the component to be detected under the expected fault condition.

[0154] In one embodiment, the phase angle increment module 601 is also used to input the node susceptance matrix and the total increment of active power flow into a preset voltage phase angle increment matrix operation model to obtain the voltage phase angle increment matrix of the component to be detected under the expected fault condition; according to the voltage phase angle increment matrix, determine the voltage phase angle increment corresponding to each bus in the component to be detected.

[0155] In one embodiment, the reactance calculation module 602 is also used to determine the first bus and the second bus corresponding to the current branch; use the voltage phase angle increment of the first bus as the first voltage phase angle increment, and use the voltage phase angle increment of the second bus as the second voltage phase angle increment; based on the first bus and the second bus, obtain the reactance value of the current branch.

[0156] In one embodiment, the power flow increment module 603 is also used to obtain the voltage phase angle increment difference corresponding to the first voltage phase angle increment and the second voltage phase angle increment based on the first voltage phase angle increment and the second voltage phase angle increment; input the voltage phase angle increment difference and the reactance value of the current branch into a preset active power flow increment calculation model to obtain the active power flow increment of the current branch.

[0157] In one embodiment, the phase angle increment module 601 is also used to obtain the vector of the voltage phase angle of each node in the element to be detected; the vector of the voltage phase angle of each node in the element to be detected is input into a preset active power flow total increment calculation model to obtain the active power flow total increment.

[0158] Each module in the above-mentioned anticipated fault scanning device for power system components can be implemented in whole or in part by software, hardware, or a combination thereof. Each of the above-mentioned modules can be embedded in or independent of a processor in a computer device in the form of hardware, or can be stored in a memory in a computer device in the form of software, so that the processor can call and execute operations corresponding to each of the above modules.

[0159] In one embodiment, a computer device is provided, which may be a server, and its internal structure diagram may be shown in FIG7. The computer device includes a processor, a memory, an input / output interface (Input / Output, referred to as I / O), and a communication interface. The processor, the memory, and the input / output interface are connected via a system bus, and the communication interface is connected to the system bus via the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program, and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal via a network connection. When the computer program is executed by the processor, a method for scanning anticipated faults applied to power system components is implemented.

[0160] Those skilled in the art will understand that the structure shown in FIG. 7 is merely a block diagram of a partial structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than shown in the figure, or combine certain components, or have a different arrangement of components.

[0161] In one embodiment, a computer device is further provided, including a memory and a processor, wherein a computer program is stored in the memory, and the processor implements the steps in the above method embodiments when executing the computer program.

[0162] In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the steps in the above-mentioned method embodiments are implemented.

[0163] In one embodiment, a computer program product is provided, including a computer program, which implements the steps in the above method embodiments when executed by a processor.

[0164] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with relevant laws, regulations and standards of relevant countries and regions.

[0165] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to the memory, database or other medium used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetoresistive random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. As an illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM). The database involved in each embodiment provided in this application may include at least one of a relational database and a non-relational database. Non-relational databases may include distributed databases based on blockchains, etc., but are not limited to this. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., but are not limited to this.

[0166] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0167] The above-described embodiments only express several implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the present application. It should be pointed out that, for a person of ordinary skill in the art, several variations and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the protection scope of the present application shall be subject to the attached claims.

Claims

1. A method for scanning anticipated faults of power system components, characterized in that: The method comprises: Obtaining a node susceptance matrix of the element to be detected under the expected fault condition, determining the total increment of active power flow of the element to be detected under the expected fault condition, and obtaining a voltage phase angle increment corresponding to each bus in the element to be detected according to the node susceptance matrix and the total increment of active power flow; For each branch in the element to be detected, determining a first voltage phase angle increment, a second voltage phase angle increment and a reactance value of the current branch corresponding to the current branch; Obtaining an active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch; Obtaining the active power flow of the current branch according to the base power flow of the current branch and the active power flow increment; The weight coefficient and power flow limit value of the current branch are determined, and the active power flow, weight coefficient and power flow limit value are input into a preset severity index calculation model to obtain the expected fault severity index when the current branch is disconnected.

2. The method according to claim 1, characterized in that The step of obtaining the node susceptance matrix of the component to be detected under the expected fault condition includes: Determine an initial matrix corresponding to the base state power flow of the current branch; According to the initial matrix, a matrix after the base state power flow changes under the expected fault condition is obtained; The changed matrix is ​​corrected by using a matrix correction method to obtain a node susceptance matrix of the element to be detected under the expected fault condition.

3. The method according to claim 1, characterized in that The step of obtaining the voltage phase angle increment corresponding to each bus in the element to be detected according to the node susceptance matrix and the total increment of active power flow includes: Inputting the node susceptance matrix and the total increment of active power flow into a preset voltage phase angle increment matrix operation model to obtain the voltage phase angle increment matrix of the element to be detected under the expected fault condition; The voltage phase angle increment corresponding to each bus in the element to be detected is determined according to the voltage phase angle increment matrix.

4. The method according to claim 3, characterized in that The determining of a first voltage phase angle increment, a second voltage phase angle increment and a reactance value of the current branch corresponding to the current branch includes: Determine a first bus and a second bus corresponding to the current branch; Using the voltage phase angle increment of the first bus as the first voltage phase angle increment, and using the voltage phase angle increment of the second bus as the second voltage phase angle increment; Based on the first bus and the second bus, the reactance value of the current branch is obtained.

5. The method according to claim 4, characterized in that The step of obtaining the active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value of the current branch comprises: According to the first voltage phase angle increment and the second voltage phase angle increment, a voltage phase angle increment difference corresponding to the first voltage phase angle increment and the second voltage phase angle increment is obtained; The voltage phase angle increment difference and the reactance value of the current branch are input into a preset active power flow increment calculation model to obtain the active power flow increment of the current branch.

6. The method according to any one of claims 1 to 5, characterized in that The determining of the total increment of active power flow of the component to be detected under the expected fault condition comprises: Obtaining a vector of voltage phase angles of each node in the element to be detected; The vector of the voltage phase angle of each node in the to-be-detected element is input into a preset active power flow total increment calculation model to obtain the active power flow total increment.

7. A device for scanning anticipated faults of power system components, characterized in that: The device comprises: A phase angle increment module is used to obtain the node susceptance matrix of the element to be detected under the expected fault condition, determine the total increment of active power flow of the element to be detected under the expected fault condition, and obtain the voltage phase angle increment corresponding to each bus in the element to be detected according to the node susceptance matrix and the total increment of active power flow; A reactance calculation module, used to determine, for each branch in the element to be detected, a first voltage phase angle increment, a second voltage phase angle increment corresponding to a current branch, and a reactance value of the current branch; A power flow increment module, used for obtaining the active power flow increment of the current branch according to the first voltage phase angle increment, the second voltage phase angle increment and the reactance value; An active power flow module, used to obtain the active power flow of the current branch according to the base power flow of the current branch and the active power flow increment; The index determination module is used to determine the weight coefficient and power flow limit value of the current branch, input the active power flow, weight coefficient and power flow limit value into a preset severity index calculation model, and obtain the expected fault severity index when the current branch is disconnected.

8. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 6 are implemented.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

10. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

Images