CN109167354B - Power grid expected fault parallel analysis and calculation method based on file exchange - Google Patents
Power grid expected fault parallel analysis and calculation method based on file exchange Download PDFInfo
- Publication number
- CN109167354B CN109167354B CN201811166572.1A CN201811166572A CN109167354B CN 109167354 B CN109167354 B CN 109167354B CN 201811166572 A CN201811166572 A CN 201811166572A CN 109167354 B CN109167354 B CN 109167354B
- Authority
- CN
- China
- Prior art keywords
- calculation
- resource
- task
- computing
- file
- 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.)
- Active
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
Abstract
The invention relates to a power grid fault parallel analysis and calculation method based on file exchange, which is technically characterized by comprising the following steps of: the method comprises the following steps: step 1, establishing a task table, a computing resource table and a system fault set; establishing computing resource information and computing task information according to the participating computing resource server; step 2, running a total task scheduling process on a task scheduling server, adding a new power grid static safety analysis calculation total task in a task table when each round of static safety analysis function running period comes, decomposing the total task calculated according to expected faults of static safety analysis into a specific parallel calculation task of each calculation resource, and writing the parallel calculation task into a calculation resource table; step 3, on the computing resource server, computing resource daemon process; and 4, waiting for the next static security analysis calculation period to arrive on the task scheduling server. The invention can improve the online analysis and calculation speed of the static security analysis function.
Description
Technical Field
The invention belongs to the technical field of power systems, and relates to a power grid expected fault parallel analysis and calculation method, in particular to a power grid expected fault parallel analysis and calculation method based on file exchange.
Background
The static safety analysis of the power grid is one of main application software of a modern power grid operation control center, and the main aim of the static safety analysis is to preset faults which may occur in the power grid according to the needs of personnel, accurately give out the operation mode of a system after the faults by analyzing and calculating the power grid operation conditions after the faults occur on line, and visually and accurately display various fault results, so that the damage of the faults to the power grid operation is quickly judged, and operation on-duty personnel is prompted to collect and prevent errors in time. The main core function of the static safety analysis of the power grid is to perform rapid on-line calculation on various expected power grid faults, the calculation result gives an overload branch, a voltage abnormal bus, a pressure abnormal node, an out-of-limit generator, a thermal power plant and the like under the expected faults, the out-of-limit degree of the overload branch, the voltage abnormal bus, the pressure abnormal node, the out-of-limit generator, the thermal power plant and the like is given, and an analysis calculation basis is provided for guaranteeing the stable operation safety and reliability of the power system.
With the continuous and rapid increase of the power demand in China and the rapid development of the power industry, the scale of the power grid is continuously enlarged. The extensive research on the blackout accidents at home and abroad in recent years shows that the blackout accidents with great influence are gradually evolved from static safety constraint damage: power transmission elements such as power system lines, transformers and the like break down and quit operation, large-range power transfer is caused, overload of adjacent elements is caused, linkage overload tripping is formed, and finally large-area power failure accidents are caused. Therefore, the calculation of expected faults in the static safety analysis of the power system has become an important means for ensuring the safe operation of the power system.
A common power grid fault analysis and calculation method in static safety analysis is a load flow calculation method. The tidal current calculation is a calculation for determining steady-state operation state parameters of each part of the power system according to the given power grid structure, parameters and operation conditions of elements such as a generator and a load. Typically given operating conditions there are power at each source and load point in the system, pivot point voltage, voltage at the balance point and phase angle. The operating state parameters to be solved comprise voltage amplitude and phase angle of each bus node of the power grid, power distribution of each branch circuit, power loss of the network and the like. The fault analysis is that after initial operation conditions of load flow calculation are given, influences of expected faults on the operation conditions are given, for example, faults of lines, main transformers, generators and the like exit the operation, and then the load flow calculation is carried out, so that the voltage of each bus node and the power of each branch circuit after the expected faults occur are obtained.
At present, monitoring systems of power grid dispatching control centers above province level all have a static safety analysis function. The scanning period of the online static security analysis in the dispatching center is generally 1-15 minutes, and most of the required periods are 1 minute, so that the calculation speed of the online static security analysis is urgently needed to be improved; meanwhile, with the rapid development of power grids in China, the scale of the power grids monitored by each level of dispatching centers is larger and larger, the number of calculation nodes in a power grid calculation model is mostly more than 4000, the number of faults for performing on-line expected fault calculation is more than 8000, the large-scale calculation also provides higher requirements for the speed of an on-line expected fault analysis and calculation method, and the traditional existing method cannot meet the requirements.
Disclosure of Invention
The invention aims to overcome the defects of the prior art, provides a power grid expected fault parallel analysis and calculation method based on file exchange, and improves the online analysis and calculation speed of a static safety analysis function.
The invention solves the practical problem by adopting the following technical scheme:
a power grid fault parallel analysis and calculation method based on file exchange comprises the following steps:
step 1, before calculation, establishing a task table, a calculation resource table and a system fault set; establishing computing resource information and computing task information according to the participating computing resource server;
step 2, running a total task scheduling process on a task scheduling server, adding a new power grid static safety analysis calculation total task in a task table when each round of static safety analysis function running period comes, decomposing the total task calculated according to expected faults of static safety analysis into a specific parallel calculation task of each calculation resource, writing the parallel calculation task into a calculation resource table, and preparing calculation data;
step 3, on the computing resource server, the computing resource daemon starts a corresponding computing resource computing process according to the task allocation state of each computing resource updated in the parallel computing resource table;
and 4, on the task scheduling server, when waiting for the next static security analysis calculation period to arrive, returning to execute the step 2.
Further, the specific steps of step 1 include:
(1-1) establishing a task table, comprising: the system comprises a task name, a task type, task starting personnel, a task computing module, a task computing function, task overtime time, a task computing family, a task current state, a task state change time scale, task starting time, task completion time, a task completion ratio and task current operation information;
(1-2) building a computing resource table, comprising: the resource calculation method comprises the following steps of resource name, resource server name, resource module function name, whether to be put into operation, resource calculation module process, resource calculation task state and resource current day failure times, wherein the resource calculation task state can be as follows: "idle", "release complete", "compute failure", "exception";
(1-3) establishing records of a computing resource table according to the number n of servers participating in parallel computing, wherein the records are recorded as a set R ═ { Ri, i ═ 1, …, n }, and n is the number of configured servers adopting the technology, namely the number of parallel computing resources; setting the resource task state of each computing resource REi to be idle initially;
(1-4) establishing a grid fault set according to the grid condition, wherein the grid fault set is recorded as NF (NFx, x is 1, …, n 2), and n2 is the maximum fault number;
further, the specific steps of step 2 include:
(2-1) firstly, generating a model file Fz for power grid load flow calculation, wherein the file comprises the connection relation of each branch in the power grid and load flow injection data of each node, and the model file is stored in a compression mode;
(2-2) scanning and checking the calculation resource table, finding out all resources with the resource task state set as 'free', recording the resources as a set RE ═ { REi, i ═ 1, …, n1}, wherein the total number of the resources is n 1;
(2-3) automatically generating n2 calculation subtasks according to the power grid fault set NF: and TF, TFj, j is 1, …, n2, and the distribution method is that according to the average distribution mode, the number of calculation faults borne by each calculation subtask is as follows: n3 ═ n2/n1+ 1; n3 number of computing failures assumed for each computing subtask; in order to generate a task file TFile of each computing subtask, the file is recorded as TFF { TFFj, j ═ 1, …, n2}, each file TFFj contains n3 faults, the file is recorded as TFFj-k, k ═ 1, …, n3, and the file adopts a text file described by an E language and contains resource state information and expected fault set information;
(2-4) performing task allocation on each computation subtask TFj in the set TF;
(2-5) scanning and checking the calculation resource table R, and judging the resource task state of each resource Ri in the table;
moreover, the specific steps of the step 2 (2-3) include:
(2-3-1) generating resource state information in the TFFj file, wherein the class name is calsc _ info and the attribute is described in the following table; each subtask corresponds to a line of data in the data block calsc _ info;
(2-3-2) generating an expected failure information class name flt _ info in the TFFj file, wherein the attribute is described in the following table; each expected failure TFFj-k corresponds to a row of data in the data block flt _ info;
moreover, the specific steps of the (2-4) th step of the step 2 include:
(2-4-1) scanning and checking the calculation resource table R, finding out a resource with a certain resource task state set as 'idle', and recording the resource as Ri;
(2-4-2) allocating the subtasks TFj to the computing resources Ri by: sending a load flow calculation model file Fz and a calculation file TFFj to a server for calculating the resource Ri, and setting the resource task state of the calculation resource Ri as 'release completion';
moreover, the specific steps of the (2-5) th step of the step 2 include:
(2-5-1) if the resource task state is 'calculation completion', recycling a resource Ri calculation result file FRsti, namely transmitting the result file FRsti to a task scheduling server, and displaying the content in each result file on a human-computer interface; after the result file FRsti is recovered, the resource task state of the resource Ri is set as 'idle';
(2-5-2) if the resource task state is 'computing failure', re-executing the step (2-4) of the step 2 on the subtask TFj allocated to the resource to re-allocate the computing task; meanwhile, accumulating the calculation failure times m of the resource Ri on the same day, and setting the calculation task state of the resource Ri to be idle when m is smaller than the specified failure time threshold value N; when m is larger than or equal to N, setting the state of the calculation task of the resource Ri to be abnormal, and enabling the resource Ri not to participate in the subsequent calculation task in the day any more;
further, the specific steps of step 3 include:
(3-1) judging the resource task state of the resource record Ri corresponding to the current server in the computing resource table, if the computing task state of the resource is 'release completion', starting a fault analysis computing process to compute, and executing the step (3-2) of the subsequent step 3; otherwise, no step is executed;
(3-2) reading in a newly received model file Fz of power grid load flow calculation to form a power grid calculation model;
(3-3) reading in the newly received task file TFFj, and for each failure TFFj-k, executing the following steps:
(3-3-1) simulating and setting the fault of the branch or the node in the power grid model;
(3-3-2) carrying out power grid load flow calculation, if the calculation is successful, executing the step (3-3-3), and if the calculation is failed, directly executing the step (3-6) of the step 3 and terminating the calculation;
(3-3-3) counting and recording branch and node out-of-limit conditions of the power grid after the fault according to the load flow calculation result to form a calculation result TFFRsti, and storing the calculation result TFFRsti into a memory table of the computer;
(3-3-4) restoring the branch and node states set by the fault to initial states;
(3-3-5) returning to the step 3-3-1) until all fault calculations are completed;
(3-4) after all the fault calculation is finished, generating a calculation result, and writing the calculation result into a calculation result file FRsti of the current calculation resource Ri;
(3-4-1) generating out-of-limit information of the fault-caused voltage, class name in the calculation result file Frsti: a fall _ v _ rst including the following attribute information; each fault causes the out-of-limit of the voltage to correspond to one record of the fall _ v _ rst data block;
(3-4-2) generating out-of-limit information of the branch caused by the fault in the calculation result file Frstill, class name: fall _ q _ rst, which contains the following attribute information; the out-of-limit of each branch caused by each fault corresponds to a record of a fall _ q _ rst data block;
(3-5) setting the resource task state of the computing resource Ri corresponding to the current server as 'computing completed', and returning to the (3-1) of the step (3);
and (3-6) after the calculation is failed, writing the calculation resource Ri corresponding to the current server into the calculation resource Ri, wherein the task state is 'calculation failure', and returning to the (3-1) in the step (3).
The invention has the advantages and beneficial effects that:
1. the invention relates to a power grid fault parallel analysis and calculation method based on file exchange. The human-computer interface is responsible for configuring data of the computing resource table and the task table; the task management is operated on a main server, the computing resource can be a plurality of servers, a fault analysis computing function is installed on each server, and the plurality of servers can realize parallel computing. The task management layer and the computing resource layer transmit interactive data through files, and meanwhile, the task management layer can also receive messages sent by the human-computer interface and respond to human-computer interaction. For the 'failure analysis' module, a task table shared by a network is used for driving a computing task, and input and output files are used as data interaction media. The method can shorten the total time of fault analysis calculation and improve the online analysis calculation speed of the static safety analysis function.
2. The invention defines the task table and the computing resource table in the memory of the task scheduling server, and defines the E file as the input and output data interaction medium, the task scheduling server uniformly distributes and issues the tasks, and simultaneously, a plurality of computing resource servers perform parallel computing, thereby greatly shortening the computing time; the online analysis and calculation speed of the static security analysis function is improved, and the minute-level calculation requirement is met; and then the stability of the power grid is improved.
Drawings
FIG. 1 is a flow chart of task server execution of the present invention;
fig. 2 is a flow chart of the resource server implementation of the present invention.
Detailed Description
The embodiments of the invention will be described in further detail below with reference to the accompanying drawings:
the electric power system data markup language-E language specification (hereinafter referred to as "E language") is a novel and efficient electric power system data markup language which is formulated for the efficiency defect of a common data model (CIM) of an IEC 61970-301 electric power system on the basis of object-oriented abstraction and aims at the efficiency defect of the CIM when the CIM is described in an XML mode. The E language data file includes comments (optional), system declarations, data block start tags, data block header definitions, data blocks, data block end tags, and the like. Generally classified as class name definitions, attribute definitions, etc.; the national grid company enterprise standard Q-GDW-215-2008 power system data markup language-E language specification specifies the description mode of the E language, and the specific format is defined as follows:
1) note
Comments are an optional part of the E language, starting with a parallel-slash "//" guide. Annotations may appear anywhere in the data file, either as a separate line or as allowed at the back of the line.
2) System declarations
The system declaration of E language is located at the front end of the data block and is juxtaposed by a left angle bracket plus an exclamation mark "<! "guide Start, side by side with a sigh plus right braces! "end of mark". The name of the application system generating the data, the E language version, the character set and the data version are stated in the system declaration. The specific content within the tag is as follows:
a) system is an application System name. Refers to the name of the application system that generated the data.
b) Version ═ E language Version. Refers to the E language version used to encapsulate the data. Due to the development of the E language itself, grammar rules
And the symbol may be changed, this version information can be used for the difference of the program recognition language itself.
c) Code is the name of the character set. Refers to the character set used in this document.
d) Data is a version of the Data. The version of the data is referred to. Due to the change of application requirements, the packaging structure of the same type of data may be locally expanded, modified or simplified, and the version information can be used for program identification of slight differences of the data.
The above parts are divided by spaces and can be expanded according to the system condition. The following is a practical example of a system declaration: < | A! System 1.0 Code UTF-8 Data 1.0 | System SYS Version! < CHEM > A
3) Data block start and end
The start of the E language data block is represented by an angle bracket plus a horizontal slash in the angle bracket, wherein the angle bracket is used for indicating the start of the E language data block, and the angle bracket is used for indicating the end of the E language data block. In "class name:" entity name ", the former indicates the class of the data object in the data block, and the latter indicates the belonging object of the data in the data block.
4) Data block header definition
The data block header definition is marked with an address symbol "@" for indicating the basic structure of the data: i.e. a cross-table type, a single-row type or a multi-row type. The specific notation of the data block header definition is as follows: the single address character "@" represents that the basic structure of the data is in a horizontal table form, and the syntax is as follows: @ serial number attribute name 1 attribute name 2 attribute name 3 … "@", serial number, and space division among attribute names may be used, the serial number may be represented by a data ID (the same applies hereinafter), and the processing manner of the single-column and multi-column structure with respect to the serial number is similar. Each object in the data block of the horizontal table type structure occupies one row, and each attribute occupies one column, so that the data block is suitable for table type data or data with fewer attributes and more objects.
5) Data block
The data block is composed of a plurality of data lines, each of which is guided by a pound sign "#". Each data value in the data line is divided by a space. The expression mode of the data line corresponds to the definition of the data block head one by one. The first column of data blocks is always followed by a "#" followed by one or more spaces, which are followed by a sequence number, e.g., "# 10", to make the sequence number a single column. The sequence number indicates the order of the line data in the data block, starting with a 1. The sequence number is followed by a space before the sequence of data values begins.
The text interaction in the invention adopts E language description, and the data block definition adopts a horizontal table type, namely each class is a table, and each table has different attributes.
A power grid fault parallel analysis and calculation method based on file exchange comprises the following steps:
step 1, before calculation, establishing a task table, a calculation resource table and a system fault set; the task table and the computing resource table are established in a memory of the server, and other resource servers can access data in the task table and the computing resource table through network interfaces; establishing computing resource information and computing task information according to the participating computing resource server;
the specific steps of the step 1 comprise:
(1-1) creating a TASK table (DEF _ TASK) including: task name, task type (automatic calculation and manual calculation), task starting personnel, task calculation module, task calculation function, task overtime time, task calculation family, task current state, task state change time scale, task starting time, task completion ratio and task current operation information;
(1-2) building a computation resource table (DEF _ CALSC) comprising: the resource calculation method comprises the following steps of resource name, resource server name, resource module function name, whether to be put into operation, resource calculation module process, resource calculation task state and resource current day failure times, wherein the resource calculation task state can be as follows: "idle", "release complete", "compute failure", "exception";
(1-3) establishing records of a computing resource table according to the number n of servers participating in parallel computing, wherein the records are recorded as a set R ═ { Ri, i ═ 1, …, n }, and n is the number of configured servers adopting the technology, namely the number of parallel computing resources; setting the resource task state of each computing resource REi to be idle initially;
(1-4) establishing a grid fault set according to the grid condition, wherein the grid fault set is recorded as NF (NFx, x is 1, …, n 2), and n2 is the maximum fault number;
step 2, running a total task scheduling process on a task scheduling server, adding a new power grid static safety analysis calculation total task in a task table when each round of static safety analysis function running period comes, decomposing the total task calculated according to expected faults of static safety analysis into a specific parallel calculation task of each calculation resource, writing the parallel calculation task into a calculation resource table, and preparing calculation data;
as shown in fig. 1, the specific steps of step 2 include:
(2-1) firstly, generating a model file Fz for power grid load flow calculation, wherein the file comprises the connection relation of each branch in the power grid and load flow injection data of each node, and the model file is stored in a compression mode;
(2-2) scanning and checking the calculation resource table, finding out all resources with the resource task state set as 'free', recording the resources as a set RE ═ { REi, i ═ 1, …, n1}, wherein the total number of the resources is n 1;
(2-3) automatically generating n2 calculation subtasks according to the power grid fault set NF: and TF, TFj, j is 1, …, n2, and the distribution method is that according to the average distribution mode, the number of calculation faults borne by each calculation subtask is as follows: n3 ═ n2/n1+ 1; n3 number of computing failures assumed for each computing subtask; in order to generate a task file TFile of each computing subtask, the file is recorded as TFF { TFFj, j ═ 1, …, n2}, each file TFFj contains n3 faults, the file is recorded as TFFj-k, k ═ 1, …, n3, and the file adopts a text file described by an E language and contains resource state information and expected fault set information;
the step 2, the step (2-3) comprises the following specific steps:
(2-3-1) generating resource state information in the TFFj file, a class name (calsc _ info), and attributes described in the following table; each subtask corresponds to a line of data in the data block calsc _ info;
| serial number | Properties | Data type | Description of the invention |
| 1 | Resource module name | Character string | |
| 2 | Resource status | Shaping of | |
| 3 | Current computing task ID | Long plastic | |
| 4 | Current computing task function | Character string | |
| 5 | Current computing task start time | Shaping of | |
| 6 | Current computing task end time | Shaping of | |
| 7 | Current computing task completion ratio | Floating point type |
(2-3-2) generating an expected failure information class name (flt _ info) in the TFFj file, the attributes of which are described in the following table; each expected failure TFFj-k corresponds to a row of data in the data block flt _ info;
(2-4) performing task allocation on each computation subtask TFj in the set TF;
the specific steps of the (2-4) step of the step 2 comprise:
(2-4-1) scanning and checking the calculation resource table R, finding out a resource with a certain resource task state set as 'idle', and recording the resource as Ri;
(2-4-2) allocating the subtasks TFj to the computing resources Ri by: sending a load flow calculation model file Fz and a calculation file TFFj to a server for calculating the resource Ri, and setting the resource task state of the calculation resource Ri as 'release completion';
(2-5) scanning and checking the calculation resource table R, and judging the resource task state of each resource Ri in the table;
the specific steps of the (2) th step to the (5) th step of the step 2 comprise:
(2-5-1) if the resource task state is 'calculation completion', recycling a resource Ri calculation result file FRsti, namely transmitting the result file FRsti to a task scheduling server, and displaying the content in each result file on a human-computer interface (see step 3) for a result file format); after the result file FRsti is recovered, the resource task state of the resource Ri is set as 'idle';
(2-5-2) if the resource task state is 'computing failure', re-executing the step (2-4) of the step 2 on the subtask TFj allocated to the resource to re-allocate the computing task; meanwhile, accumulating the calculation failure times m of the resource Ri on the same day, and setting the calculation task state of the resource Ri to be idle when m is smaller than the specified failure time threshold value N; when m is larger than or equal to N, setting the state of the calculation task of the resource Ri to be abnormal, and enabling the resource Ri not to participate in the subsequent calculation task in the day any more;
and 3, on the computing resource server, starting the corresponding computing resource computing process by the computing resource daemon according to the task allocation state of each computing resource updated in the parallel computing resource table.
As shown in fig. 2, the specific steps of step 3 include:
(3-1) judging the resource task state of the resource record Ri corresponding to the current server in the computing resource table, if the computing task state of the resource is 'release completion', starting a fault analysis computing process to compute, and executing the step (3-2) of the subsequent step 3; otherwise, no step is executed;
(3-2) reading in a newly received model file Fz of power grid load flow calculation to form a power grid calculation model;
(3-3) reading in the newly received task file TFFj, and for each failure TFFj-k, executing the following steps:
(3-3-1) simulating and setting the fault of the branch or the node in the power grid model;
(3-3-2) carrying out power grid load flow calculation, if the calculation is successful, executing the step (3-3-3), and if the calculation is failed, directly executing the step (3-6) of the step 3 and terminating the calculation;
(3-3-3) counting and recording branch and node out-of-limit conditions of the power grid after the fault according to the load flow calculation result to form a calculation result TFFRsti, and storing the calculation result TFFRsti into a memory table of the computer;
(3-3-4) restoring the branch and node states set by the fault to initial states;
(3-3-5) returning to the step 3-3-1) until all fault calculations are completed.
(3-4) after all the fault calculation is finished, generating a calculation result, and writing the calculation result into a calculation result file FRsti of the current calculation resource Ri;
(3-4-1) generating out-of-limit information of the fault-caused voltage, class name in the calculation result file Frsti: a fall _ v _ rst including the following attribute information; each fault causes the out-of-limit of the voltage to correspond to one record of the fall _ v _ rst data block;
| serial number | Properties | Data type | Description of the invention |
| 1 | Fault ID | Long plastic | |
| 2 | Name of failure | Character string | |
| 3 | The upper limit of bus bar quantity caused by fault | Shaping of | |
| 4 | Maximum-over-limit bus name | Character string | |
| 5 | Maximum and maximum bus voltage | Floating point type | |
| 6 | Lower limit of bus number caused by fault | Shaping of | |
| 7 | Lowest-lower limit bus name | Character string | |
| 8 | Lowest lower limit bus voltage value | Floating point type |
(3-4-2) generating out-of-limit information of the branch caused by the fault in the calculation result file Frstill, class name: fall _ q _ rst, which contains the following attribute information; each failure causes the violation of a branch to correspond to a record of the fall _ q _ rst data block.
| Serial number | Properties | Data type | Description of the invention |
| 1 | Fault ID | Long plastic | |
| 2 | Name of failure | Character string | |
| 3 | Out-of-limit number of branches caused by fault | Shaping of | |
| 4 | Highest-over-bound branch name | Character string | |
| 5 | Maximum upper limit branch load flow active value | Floating point type | |
| 6 | Maximum upper limit branch load flow reactive value | Floating point type |
(3-5) setting the resource task state of the computing resource Ri corresponding to the current server as 'computing completed', and returning to the (3-1) of the step (3);
(3-6) after the calculation is failed, writing the calculation resource Ri task state of the calculation resource Ri corresponding to the current server into the calculation resource Ri, namely 'calculation failure', and returning to the (3-1) in the step (3);
and 4, on the task scheduling server, when waiting for the next static security analysis calculation period to arrive, returning to execute the step 2.
The working principle of the invention is as follows:
the invention provides a power grid fault parallel analysis and calculation method based on file exchange, which defines a task table and a calculation resource table in a task scheduling server memory, uniformly distributes and issues tasks by defining an E file as an input and output data interaction medium, and simultaneously performs parallel calculation by a plurality of calculation resource servers, thereby greatly shortening the calculation time; the online analysis and calculation speed of the static security analysis function is improved, and the minute-level calculation requirement is met; and then the stability of the power grid is improved.
It should be emphasized that the examples described herein are illustrative and not restrictive, and thus the present invention includes, but is not limited to, those examples described in this detailed description, as well as other embodiments that can be derived from the teachings of the present invention by those skilled in the art and that are within the scope of the present invention.
Claims (6)
1. A power grid fault parallel analysis and calculation method based on file exchange is characterized in that: the method comprises the following steps:
step 1, before calculation, establishing a task table, a calculation resource table and a system fault set; establishing computing resource information and computing task information according to the participating computing resource server;
step 2, running a total task scheduling process on a task scheduling server, adding a new power grid static safety analysis calculation total task in a task table when each round of static safety analysis function running period comes, decomposing the total task calculated according to expected faults of static safety analysis into a specific parallel calculation task of each calculation resource, writing the parallel calculation task into a calculation resource table, and preparing calculation data;
step 3, on the computing resource server, computing resource daemon, and starting the corresponding computing resource computing process according to the task allocation state of each computing resource updated in the parallel computing resource table;
step 4, on the task scheduling server, when waiting for the next static security analysis calculation period to arrive, returning to execute the step 2;
the specific steps of the step 2 comprise:
(2-1) firstly, generating a model file Fz for power grid load flow calculation, wherein the file comprises the connection relation of each branch in the power grid and load flow injection data of each node, and the model file is stored in a compression mode;
(2-2) scanning and checking the calculation resource table, finding out all resources with the resource task state set as 'free', recording the resources as a set RE ═ { REi, i ═ 1, …, n1}, wherein the total number of the resources is n 1;
(2-3) automatically generating n2 calculation subtasks according to the power grid fault set NF: and TF, TFj, j is 1, …, n2, and the distribution method is that according to the average distribution mode, the number of calculation faults borne by each calculation subtask is as follows: n3 ═ n2/n1+ 1; n3 number of computing failures assumed for each computing subtask; in order to generate a task file TFile of each computing subtask, the file is recorded as TFF { TFFj, j ═ 1, …, n2}, each file TFFj contains n3 faults, the file is recorded as TFFj-k, k ═ 1, …, n3, and the file adopts a text file described by an E language and contains resource state information and expected fault set information;
(2-4) performing task allocation on each computation subtask TFj in the set TF;
and (2-5) scanning and checking the calculation resource table R, and judging the resource task state of each resource Ri in the table.
2. The file exchange-based power grid fault parallel analysis and calculation method according to claim 1, characterized in that: the specific steps of the step 1 comprise:
(1-1) establishing a task table, comprising: the system comprises a task name, a task type, task starting personnel, a task computing module, a task computing function, task overtime time, a task computing family, a task current state, a task state change time scale, task starting time, task completion time, a task completion ratio and task current operation information;
(1-2) building a computing resource table, comprising: the resource calculation method comprises the following steps of resource name, resource server name, resource module function name, whether to be put into operation, resource calculation module process, resource calculation task state and resource current day failure times, wherein the resource calculation task state can be as follows: "idle", "release complete", "compute failure", "exception";
(1-3) establishing records of a computing resource table according to the number n of servers participating in parallel computing, wherein the records are recorded as a set R ═ { Ri, i ═ 1, …, n }, and n is the number of configured servers adopting the technology, namely the number of parallel computing resources; setting the resource task state of each computing resource REi to be idle initially;
and (1-4) establishing a grid fault set according to the grid condition, wherein NF is recorded as { NFx, x is 1, …, n2}, and n2 is the maximum fault number.
3. The file exchange-based power grid fault parallel analysis and calculation method according to claim 1, characterized in that: the step 2, the step (2-3) comprises the following specific steps:
(2-3-1) generating resource state information in the TFFj file, wherein the class name is calsc _ info and the attribute is described in the following table; each subtask corresponds to a line of data in the data block calsc _ info;
(2-3-2) generating an expected failure information class name flt _ info in the TFFj file, wherein the attribute is described in the following table; each expected failure TFFj-k corresponds to a row of data in the data block flt _ info.
4. The file exchange-based power grid fault parallel analysis and calculation method according to claim 1, characterized in that: the specific steps of the (2-4) step of the step 2 comprise:
(2-4-1) scanning and checking the calculation resource table R, finding out a resource with a certain resource task state set as 'idle', and recording the resource as Ri;
(2-4-2) allocating the subtasks TFj to the computing resources Ri by: and sending the load flow calculation model file Fz and the calculation file TFFj to a server of the calculation resource Ri, and setting the resource task state of the calculation resource Ri as 'release completion'.
5. The file exchange-based power grid fault parallel analysis and calculation method according to claim 1, characterized in that: the specific steps of the (2) th step to the (5) th step of the step 2 comprise:
(2-5-1) if the resource task state is 'calculation completion', recycling a resource Ri calculation result file FRsti, namely transmitting the result file FRsti to a task scheduling server, and displaying the content in each result file on a human-computer interface; after the result file FRsti is recovered, the resource task state of the resource Ri is set as 'idle';
(2-5-2) if the resource task state is 'computing failure', re-executing the step (2-4) of the step 2 on the subtask TFj allocated to the resource to re-allocate the computing task; meanwhile, accumulating the calculation failure times m of the resource Ri on the same day, and setting the calculation task state of the resource Ri to be idle when m is smaller than the specified failure time threshold value N; when m is larger than or equal to N, the state of the calculation task of the resource Ri is abnormal, and the resource Ri does not participate in the subsequent calculation task in the current day any more.
6. The file exchange-based power grid fault parallel analysis and calculation method according to claim 1, characterized in that: the specific steps of the step 3 comprise:
(3-1) judging the resource task state of the resource record Ri corresponding to the current server in the computing resource table, if the computing task state of the resource is 'release completion', starting a fault analysis computing process to compute, and executing the step (3-2) of the subsequent step 3; otherwise, no step is executed;
(3-2) reading in a newly received model file Fz of power grid load flow calculation to form a power grid calculation model;
(3-3) reading in the newly received task file TFFj, and for each failure TFFj-k, executing the following steps:
(3-3-1) simulating and setting the fault of the branch or the node in the power grid model;
(3-3-2) carrying out power grid load flow calculation, if the calculation is successful, executing the step (3-3-3), and if the calculation is failed, directly executing the step (3-6) of the step 3 and terminating the calculation;
(3-3-3) counting and recording branch and node out-of-limit conditions of the power grid after the fault according to the load flow calculation result to form a calculation result TFFRsti, and storing the calculation result TFFRsti into a memory table of the computer;
(3-3-4) restoring the branch and node states set by the fault to initial states;
(3-3-5) returning to the step (3-3-1) until all fault calculations are completed;
(3-4) after all the fault calculation is finished, generating a calculation result, and writing the calculation result into a calculation result file FRsti of the current calculation resource Ri;
(3-4-1) generating out-of-limit information of the fault-caused voltage, class name in the calculation result file Frsti: a fall _ v _ rst including the following attribute information; each fault causes the out-of-limit of the voltage to correspond to one record of the fall _ v _ rst data block;
(3-4-2) generating out-of-limit information of the branch caused by the fault in the calculation result file Frstill, class name: fall _ q _ rst, which contains the following attribute information; the out-of-limit of each branch caused by each fault corresponds to a record of a fall _ q _ rst data block;
(3-5) setting the resource task state of the computing resource Ri corresponding to the current server as 'computing completed', and returning to the (3-1) of the step (3);
and (3-6) after the calculation is failed, writing the calculation resource Ri corresponding to the current server into the calculation resource Ri, wherein the task state is 'calculation failure', and returning to the (3-1) in the step (3).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811166572.1A CN109167354B (en) | 2018-10-08 | 2018-10-08 | Power grid expected fault parallel analysis and calculation method based on file exchange |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811166572.1A CN109167354B (en) | 2018-10-08 | 2018-10-08 | Power grid expected fault parallel analysis and calculation method based on file exchange |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109167354A CN109167354A (en) | 2019-01-08 |
| CN109167354B true CN109167354B (en) | 2022-02-22 |
Family
ID=64877414
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811166572.1A Active CN109167354B (en) | 2018-10-08 | 2018-10-08 | Power grid expected fault parallel analysis and calculation method based on file exchange |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109167354B (en) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101441580A (en) * | 2008-12-09 | 2009-05-27 | 华北电网有限公司 | Distributed paralleling calculation platform system and calculation task allocating method thereof |
| CN101685480A (en) * | 2008-09-27 | 2010-03-31 | 国家电力调度通信中心 | Parallel computing method for security and stability analysis of large power grid and computing platform |
| CN103325012A (en) * | 2013-06-18 | 2013-09-25 | 国家电网公司 | Parallel computing dynamic task distribution method applicable to grid security correction |
| CN103873321A (en) * | 2014-03-05 | 2014-06-18 | 国家电网公司 | Distributed file system-based simulation distributed parallel computing platform and method |
| CN104239144A (en) * | 2014-09-22 | 2014-12-24 | 珠海许继芝电网自动化有限公司 | Multilevel distributed task processing system |
| CN104993477A (en) * | 2015-06-29 | 2015-10-21 | 中国电力科学研究院 | Static safety parallel analysis method for AC-DC power grid based on configuration mode |
| CN105576648A (en) * | 2015-11-23 | 2016-05-11 | 中国电力科学研究院 | Static security analysis double-layer parallel method based on GPU-CUP heterogeneous computing platform |
| JP2016224832A (en) * | 2015-06-03 | 2016-12-28 | 富士通株式会社 | Dump management device, dump management program, and dump management method |
| CN107132762A (en) * | 2017-04-17 | 2017-09-05 | 国电南瑞科技股份有限公司 | Online Static security assessment method based on forecast failure collection automatic screening |
| CN108256263A (en) * | 2018-02-07 | 2018-07-06 | 中国电力科学研究院有限公司 | A kind of electric system hybrid simulation concurrent computational system and its method for scheduling task |
| CN108521122A (en) * | 2018-04-13 | 2018-09-11 | 深圳供电局有限公司 | A kind of Contingency Analysis of Power Systems method considering prepared auto restart |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105610157B (en) * | 2016-01-27 | 2018-01-09 | 国网山东省电力公司烟台供电公司 | A kind of area power grid static security analysis system of failure driving |
| US10355966B2 (en) * | 2016-03-25 | 2019-07-16 | Advanced Micro Devices, Inc. | Managing variations among nodes in parallel system frameworks |
| CN106058856B (en) * | 2016-07-01 | 2019-07-23 | 中国南方电网有限责任公司 | A kind of method of quick analysis power grid static security |
| CN106874113A (en) * | 2017-01-19 | 2017-06-20 | 国电南瑞科技股份有限公司 | A kind of many GPU heterogeneous schemas static security analysis computational methods of CPU+ |
-
2018
- 2018-10-08 CN CN201811166572.1A patent/CN109167354B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101685480A (en) * | 2008-09-27 | 2010-03-31 | 国家电力调度通信中心 | Parallel computing method for security and stability analysis of large power grid and computing platform |
| CN101441580A (en) * | 2008-12-09 | 2009-05-27 | 华北电网有限公司 | Distributed paralleling calculation platform system and calculation task allocating method thereof |
| CN103325012A (en) * | 2013-06-18 | 2013-09-25 | 国家电网公司 | Parallel computing dynamic task distribution method applicable to grid security correction |
| CN103873321A (en) * | 2014-03-05 | 2014-06-18 | 国家电网公司 | Distributed file system-based simulation distributed parallel computing platform and method |
| CN104239144A (en) * | 2014-09-22 | 2014-12-24 | 珠海许继芝电网自动化有限公司 | Multilevel distributed task processing system |
| JP2016224832A (en) * | 2015-06-03 | 2016-12-28 | 富士通株式会社 | Dump management device, dump management program, and dump management method |
| CN104993477A (en) * | 2015-06-29 | 2015-10-21 | 中国电力科学研究院 | Static safety parallel analysis method for AC-DC power grid based on configuration mode |
| CN105576648A (en) * | 2015-11-23 | 2016-05-11 | 中国电力科学研究院 | Static security analysis double-layer parallel method based on GPU-CUP heterogeneous computing platform |
| CN107132762A (en) * | 2017-04-17 | 2017-09-05 | 国电南瑞科技股份有限公司 | Online Static security assessment method based on forecast failure collection automatic screening |
| CN108256263A (en) * | 2018-02-07 | 2018-07-06 | 中国电力科学研究院有限公司 | A kind of electric system hybrid simulation concurrent computational system and its method for scheduling task |
| CN108521122A (en) * | 2018-04-13 | 2018-09-11 | 深圳供电局有限公司 | A kind of Contingency Analysis of Power Systems method considering prepared auto restart |
Non-Patent Citations (4)
| Title |
|---|
| Large-Scale Transient Stability Simulation of Electrical Power Systems on Parallel GPUs;Vahid Jalili-Marandi 等;《IEEE TRANSACTIONS ON PARALLEL AND DISTRIBUTED SYSTEMS》;20120731;第23卷(第7期);全文 * |
| 基于GPU并行的静态安全分析设计与应用;陆娟娟 等;《计算机科学与探索》;20161231;全文 * |
| 基于多进程并行计算的日前调度计划安全校核技术;李虎成 等;《电网与清洁能源》;20130331;第29卷(第3期);全文 * |
| 电力***在线安全分析并行分布式计算;宗永升，张祎;《计算机工程与设计》;20101231;第31卷(第19期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109167354A (en) | 2019-01-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107330056B (en) | Wind power plant SCADA system based on big data cloud computing platform and operation method thereof | |
| CN102043719B (en) | Method for testing IEC61850 configuration description file | |
| CN112799321B (en) | DBC file rapid creating method following CAN bus communication protocol | |
| DE112018000125B4 (en) | Quantifying a combined impact of interdependent uncertain resources in an electric power grid | |
| US20200334562A1 (en) | Reserve optimization method and apparatus based on support outage event constrained unit commitment | |
| CN101364284A (en) | Production management system for electrical power system | |
| CN110837496A (en) | Data quality management method and system based on dynamic sql | |
| CN109167354B (en) | Power grid expected fault parallel analysis and calculation method based on file exchange | |
| CN107069742B (en) | Electric system continuous tide computing system based on Python and PSD-BPA | |
| CN105608637A (en) | Substation equipment model topology-based wiring diagram automatic generation method | |
| CN112541011A (en) | Power grid terminal fault processing method and device based on RPA and AI | |
| CN102393820A (en) | Logical implementation method for controlling building monitoring system | |
| CN112446619B (en) | Power distribution network rush-repair processing method and device | |
| CN107748701A (en) | A kind of analysis method for reliability of electric energy measurement automation system | |
| CN115471215A (en) | Business process processing method and device | |
| CN112532432B (en) | Automatic mapping method and system based on LTU terminal | |
| CN114035783A (en) | Software code knowledge graph construction method and tool | |
| CN107391133A (en) | A kind of application generator | |
| CN113537942A (en) | Method and system for increasing number of sample marks | |
| CN112184026A (en) | Purchase material quality feedback system based on cloud platform | |
| CN110990371A (en) | Operation and maintenance management system containing knowledge base | |
| CN110262869B (en) | Private cloud Paas resource supply and demand party information interaction method | |
| CN111176257B (en) | New energy control system parameter debugging method and system | |
| CN108491303A (en) | A method of based on the O＆M automation of acceptance of the bid kylin operating system software | |
| CN117220419B (en) | Power scheduling method, system 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 | ||
| CP02 | Change in the address of a patent holder |
Address after: No. 8, Haitai Huake 4th Road, Huayuan Industrial Zone, High tech Zone, Binhai New Area, Tianjin, 300384 Patentee after: ELECTRIC POWER SCIENCE & RESEARCH INSTITUTE OF STATE GRID TIANJIN ELECTRIC POWER Co. Patentee after: STATE GRID TIANJIN ELECTRIC POWER Co. Patentee after: STATE GRID CORPORATION OF CHINA Patentee after: BEIJING KING STAR HI-TECH SYSTEM CONTROL Co.,Ltd. Address before: No.8, Haitai Huake 4th Road, Xiqing District, Tianjin 300384 Patentee before: ELECTRIC POWER SCIENCE & RESEARCH INSTITUTE OF STATE GRID TIANJIN ELECTRIC POWER Co. Patentee before: STATE GRID TIANJIN ELECTRIC POWER Co. Patentee before: STATE GRID CORPORATION OF CHINA Patentee before: BEIJING KING STAR HI-TECH SYSTEM CONTROL Co.,Ltd. |
|
| CP02 | Change in the address of a patent holder |