CN102073766A - Method for accelerating asynchronous parallel computation of transient stability of power grid based on single-computer multiple-core mode - Google Patents
Method for accelerating asynchronous parallel computation of transient stability of power grid based on single-computer multiple-core mode Download PDFInfo
- Publication number
- CN102073766A CN102073766A CN2011100046124A CN201110004612A CN102073766A CN 102073766 A CN102073766 A CN 102073766A CN 2011100046124 A CN2011100046124 A CN 2011100046124A CN 201110004612 A CN201110004612 A CN 201110004612A CN 102073766 A CN102073766 A CN 102073766A
- Authority
- CN
- China
- Prior art keywords
- transient stability
- data
- fault
- computation
- batch
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000001052 transient effect Effects 0.000 title claims abstract description 191
- 238000000034 method Methods 0.000 title claims abstract description 86
- 230000008569 process Effects 0.000 claims abstract description 54
- 238000004458 analytical method Methods 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims description 79
- 230000015572 biosynthetic process Effects 0.000 claims description 44
- 238000005755 formation reaction Methods 0.000 claims description 44
- 230000005283 ground state Effects 0.000 claims description 26
- 238000004088 simulation Methods 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 14
- 238000000354 decomposition reaction Methods 0.000 claims description 12
- 238000011160 research Methods 0.000 claims description 10
- 238000012797 qualification Methods 0.000 claims description 7
- 238000007726 management method Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000000717 retained effect Effects 0.000 claims description 5
- 230000003068 static effect Effects 0.000 claims description 4
- 230000005587 bubbling Effects 0.000 claims description 3
- 238000010224 classification analysis Methods 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 230000006641 stabilisation Effects 0.000 claims description 3
- 238000011105 stabilization Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 2
- 230000002452 interceptive effect Effects 0.000 abstract description 2
- 230000002265 prevention Effects 0.000 abstract 1
- 238000012360 testing method Methods 0.000 description 17
- 238000010205 computational analysis Methods 0.000 description 10
- 238000012545 processing Methods 0.000 description 10
- 230000005611 electricity Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 238000013439 planning Methods 0.000 description 3
- 239000003381 stabilizer Substances 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000005059 dormancy Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000011049 filling Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003012 network analysis Methods 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 230000008278 dynamic mechanism Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005662 electromechanics Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000011068 loading method Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000010223 real-time analysis Methods 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a method for quickly accelerating asynchronous parallel computation of transient stability based on a single-machine multiple-core hardware architecture, comprising the following steps: (1) realizing the batch, parallel batch, automatic load balancing batch computation of transient stability faults in fault sets, carrying out discriminatory analysis on the stability of each fault computation result, and drawing up a statistic analysis statement for fault computation results in all the fault sets; (2) carrying out batch computation and analysis on the transient stability work; and (3) providing a transient stability batch computation dead prevention method, realizing a process speed control method for parallel batch computation speed, and providing and realizing an implementation technique method combining the transient stability interactive computation and the transient stability automatic computation. The parallel batch computation is realized on the basis of the task decoupling mode without changing the transient stability computation algorithm and implementation form, the transient stability computation independence is guaranteed and the transient stability computation efficiency is improved, while the massive fault computation and the parallel computation are realized.
Description
Technical field:
The present invention relates to a kind of asynchronous parallel computing method, be specifically related to based on unit multinuclear pattern speed-raising electrical network transient stability asynchronous parallel computing method.
Background technology:
Electric system simulation is electric system scientific research, planning, design, build, put into operation, the most basic analysis and research method in the operation, maintenance work.There is multiple different manifestations type in electric system simulation by different sorting techniques, shows as off-line and in-circuit emulation by the real-time analysis of emulation, divides by the means that realize emulation to show as dynamic similation, digital simulation, digital-to-analogue hybrid analog-digital simulation emulation.Press the fineness of realistic model description and the difference of simulation step length, can show as electromechanics and electromagnetic transient simulation.What the present invention relates to is the electromechanical simulation analysis of being undertaken by computing machine.
Electromechanical simulation generally is divided into stable state computational analysis and transient state computational analysis.Stable state is calculated main pointer the trend that the digital model electrical network set up of research electrical network carries out is calculated, and is for other calculating and analyzes the ground state that makes up an operation of power networks, shows as a kind of running status that electrical network may exist.It is by equipment such as generator, excitation system, debug system, power system stabilizer, PSS, quiet dynamic load, line transformer are carried out digital modeling that transient stability calculates, calculate by the electrical network that is in certain running status is applied certain disturbance, can all generators keep the process of synchronous operation in the calculating electrical network.
Transient stability calculates the analysis means that is absolutely necessary in Power System Planning and operating analysis, have very important status.It not only for the power generation configuration of planning system and rationally send, the rationality of channel layout and the network optimization, reactive-load compensation and Control of Voltage, protection configuration and steady control means provides the check of power system transient stability; for the power system safety and stability reliability service provides quantitative basis; and can be used for studying the measure of various raising transient stabilities, calculate relay protection and aut.eq. parameter tuning definite value foundation is provided.
In the electric system, for the transmission of electricity major network, the transient stability computational tool of employing is mainly the PSD-BPA and the PSASP of China Electric Power Research Institute's exploitation at home.The two computation schema mainly is to carry out single transient stability calculating or set up the fault cascaded computation that simple fault collection is realized single mode by single fault cooperation transient stability digital model is set.For the combined fault computing of magnanimity fault and a plurality of methods of operation, adopted the parallel computation mode of computer cluster.
And the fast development of electrical network, new forms of energy are built, new electric pressure occurs, multiple complex AC/DC transmission of electricity, little operation of power networks, various power electronic devices are used, the fault collection that causes the computational analysis of electrical network transient stability, the mode amount is the trend of exploding, and a large amount of computational analysiss often has more than in fixing workplace and carries out, a large amount of computational analysiss may be in construction site, debugging is on-the-spot, meeting scene etc., all set up parallel computer group at each magnanimity calculation of fault scene, be unpractical also be uneconomic, how to solve the quick computational problem of magnanimity fault that power network development brings, be a reality also be the problem that presses for solution.
According to " PSD-BPA transient stability program user handbook is introduced, and typical transient stability calculation procedure is as follows:
1) target grid of emulation is received money, set up basic device parameter and network topology model data, utilize the trend program to set up the basic method of operation of target grid.
2) receive money target grid unit, excitation, debugging, PSS (power system stabilizer, PSS), line transformer Zero sequence parameter etc., set up basic stable model parameter, carry out tentative calculation by transient stability computation program, guarantee the rationality and the correctness of transient stability model data.
3) according to the parameter information in the tide model, set up a fault data by certain call format, fault data can be bus-bar fault, transformer fault, line fault, direct current fault, or even cascading failure.Failure process can comprise a plurality of processes, comprises that fault begins, fault disappearance, switch tripping, switch are returned, switch is jumped again etc.
4) tape jam carries out transient stability calculating, and result of calculation is carried out the judgement of stability analysis, by calculating curve of output and calculating output data and carry out the electrical network specificity analysis.
Above step is suitable for the interactive calculating of electrical network transient stability, can realize the transient stability detailed calculated of the hand-guided of single failure form.But can not adapt to fast automatic computational analysis under the magnanimity fault multimode, not utilize and calculate the multinuclear processing power, the different application mode of engineering calculation is lacked the selectivity support.
Summary of the invention:
At the deficiencies in the prior art, the present invention seeks to utilize the powerful processing power of stand-alone multinuclear, do not changing transient stability computation program itself and authorizing under the judgment mode condition, the parallel batch of realizing the operation of magnanimity fault multimode calculates, improve grid disturbance computational analysis speed, reduce grid simulation and analyze labour intensity.By setting up good a plurality of electrical networks basis method of operation in advance, by flow calculation program it is carried out precomputation and form a plurality of basic methods of operation.It is limited that trend is calculated iterations, computing velocity is very fast, and must obtain just carrying out transient stability calculating after the convergent trend method of operation, therefore need not implement the trend batch or the parallel computation of a plurality of modes in advance, here after not getting rid of the trend automatic adjustment technologies and breaking through, adopt parallel calculating method for obtaining a large amount of trend tentative calculations that the ground state method of operation carried out.Also must be according to the tide model data of setting up, by different analyze demands, be ready to magnanimity fault collection according to the form of appointment.In addition, also need according to the electrical network analysis purpose and analyze the level of detail, make up the transient stability digital model, and need to guarantee that the transient stability digital model calculates holds water under the non-fault noisy data.
After the data ready, realizing on the unit that according to the present invention the magnanimity fault collection transient stability speed-raising asynchronous parallel under the multi-operating condition calculates.
Provided by the invention based on unit multinuclear pattern speed-raising electrical network transient stability asynchronous parallel computing method, its improvements are: described method comprises the steps:
(1) determines the ground state method of operation that grid simulation calculates, set up the electrical network transient stability model parameter that ground state trend level adapts, and the electrical network that needs are studied is carried out the ground state mode prepare; Set up different classes of failure collection according to the difference of analysis purpose, and adjust electrical network transient stability model;
(2) determine ground state trend, fault collection and transient stability model corresponding relation, make up transient stability computational tasks task groups, with the TU task unit form, form the job patch task groups of whole evaluating objects, be job task batch queue data, the job batch formation is carried out by the serial sequence pattern, and it is right promptly to form transient stability compute associations data successively by the priority level of calculating;
(3) integrality of judgment data and validity are opened the transient stability operation and are calculated;
(4) be single fault data formation with current operation fault collection data decomposition;
(5) implementing selected fault collection transient stability according to the computation schema of appointment calculates in batches;
(6) transient stability that forms whole fault collection calculates the comprehensive summary form;
(7) finishing whole operations after all the transient stability computational tasks is finished calculates.
The computing method of first preferred version provided by the invention, its improvements are,
Described integrality judgement comprises the existence of ground state data existence, fault collection data, the existence of transient Model data;
Described validity is judged rationality and the validity judgement that comprises single failure in the failure collection.
The computing method of second preferred version provided by the invention, its improvements are that described computation schema is divided into:
Basic cascade serial is account form in batches;
Band is specified the parallel batch computation schema of check figure;
Dynamic auto load balance computation schema.
The computing method of the 3rd preferred version provided by the invention, its improvements are that the step of described basic cascade serial batch account form is as follows:
(a) from the fault data formation that described step (4) decomposition obtains, get one group of single failure data, form the computational data group with mode data, the pairing of transient stability model data;
(b) start transient stability with (a) computational data group and calculate nuclear;
(c) transient stability jumps to (e) after calculating and finishing;
(d) each judgement detect cycle is added up, and compares, surpass predefined threshold values, just stop transient stability computation process by force with the maximum transient stability Calculation threshold values of setting consuming time, and by force operation process recording in the calculation of fault result; Wherein, end up being a total detection time section from step (3) to All Jobs, in a total detection time section according to several detect cycles that how much are divided into of workload.
(e) the extraction result carries out the judgement of stability analysis, differentiates to be stabilization result, and promptly primary fault calculates and finishes, and related conclusions is placed in the results list, withdraws from this calculation of fault;
(f) extract fault data queuing data item, after the fault in the current computational tasks content has been calculated and finished, gather all result of calculations in this subjob, form statistical abstract, and withdraw from this subjob and calculate.
The computing method of the 4th preferred version provided by the invention, its improvements are that described parallel batch computation schema step is as follows:
1) multinuclear existence and setting value jump to 2 greater than 1);
2) set up the data queue that equates with it according to the check figure of appointment;
When 3) having idle queues in the computational data formation, from the fault data formation that described step (4) decomposition obtains, get one group of single failure data and be filled in the idle data formation, specify to be retained in the formation before the transient stability computation program application always;
4) data have been filled up in all formations, change the 5th over to) step;
5) calculating of appointment nuclear has vacant calculating nuclear, gets one group of parameter from data queue, starts transient stability computation process on the nuclear of free time, and the transient stability that carries out the fault data of getting calculates;
6) in a detect cycle, calculate end when detecting transient stability, just start resolving as a result immediately, implement power grid accident result of calculation and carry out classification analysis, interpretation of result is finished, and promptly a calculation of fault finishes, withdraw from this calculation of fault, jump to 8);
7) all parallel transient stabilities calculate the ossified qualification threshold values that duration consuming time surpasses appointment, stop surpassing the transient stability computation process of threshold values by force, jump to 8);
8) finish to calculate the calculating nuclear of release busy immediately and data parameters formation;
9) fill in the parallel data formation process,, just fill it in the vacant data queue when the prior fault operation also has unfilled fault data;
The computing method of the 5th preferred version provided by the invention, its improvements are that described ALB calculation procedure is:
I. set up the data queue that equates with it according to the check figure of appointment;
When II. having idle queues in the computational data formation, from the fault data formation that described step (4) decomposition obtains, get one group of single failure data and be filled in the idle data formation, specify to be retained in the formation before the transient stability computation program application always;
III. data have been filled up in all formations, change the V step over to;
IV. the calculating of appointment nuclear has vacant calculating nuclear, gets one group of parameter from data queue, starts transient stability computation process according to the load judgment result on the nuclear of free time, and the transient stability that carries out the fault data of getting calculates, and the dynamic load judgment formula is as follows:
The computing method of the 6th preferred version provided by the invention, its improvements are: the substep of described step (2) comprising:
1. select to calculate convergent ground state mode data;
2. select to be used for the failure collection data that transient stability calculates;
3. select to be matched with the transient stability data of ground state trend mode;
4. will be 1.-data message that 3. obtains is integrated into an independently computational tasks T
i
5. repeat 1.-4. process, set up independent computational tasks more than two, P sorts to task priority by weight factor, forms to calculate priority sequence, according to calculating priority level, arranges the sequencing of calculation task; For complete operation, with weight factor P zero setting, when new operation begins to carry out, will be again by weight factor P to task ranking, because of as calculated operation weight factor P is 0, be placed to the descending end;
Set up deque's administrative power repeated factor, the weight factor that one of them queue management is provided with, the variation of the weight factor when another is used for Management Calculation.
The computing method of the 7th preferred version provided by the invention, its improvements are, described step 5. in according to dynamic change priority weight system, realize the serializing cascade dynamic execution of computing operation in batches, it carries out weight by the [p of deque
1, p
2...., p
N]
TColumn vector control, static reaction is manual sets control, in dynamic response execution and the computation process to the not dynamic execution sequence change of computational tasks, the performance coeffcient p that complete transient stability is calculated
iZero setting to unenforced operation, when determining operation next time, realizes that by the column vector descending bubbling of computational tasks is selected.
The computing method of the 8th preferred version provided by the invention, its improvements are: the transient stability that provides in the described step (c) calculates end judgement substep and judges that the empirical value of detect cycle is 1 second.
Compared with the prior art, beneficial effect of the present invention is:
Realize three kinds of computation schemas that the magnanimity fault is calculated in batches, can seamless switching between three kinds of patterns, do not influence the transient stability computational item of having carried out after switching.
No matter transient stability fault item calculates which kind of pattern of employing, adopt the serial sequence computing method between a plurality of operations, by introducing the calculating priority level weight factor, dynamically adjust the not calculating priority of computational tasks, make the calculating of dynamic control operation successively become possibility.
The supervision judgement cycle for transient stability computation process provides the engineering experience value, and the too short or long conclusion that all will influence counting yield of theoretical analysis.The batch calculation of fault pattern of the asynchronous starting executed in parallel of the multinuclear of design can promote the efficient that transient stability calculates rapidly.
In transient stability calculates, introduce the load factor factor, the load factor that designs a calculating machine judgment model, to be fit under the routine office work condition, carry out transient stability Calculation Method in batches, can make full use of computing machine idle interval and maximum its processing power of performance, improve the transient stability counting yield.
By changing transient stability computational algorithm and way of realization, the parallel batch of realization does not calculate the task decoupling zero mode that is based in the present invention, when realizing the parallel computation of magnanimity calculation of fault, guarantees the independence of transient stability calculating itself.
The basic cascade serial serial sequenceization of account form realization transient stability fault is in batches calculated in batches, is basic form of calculation; Band specifies the parallel batch computation schema of check figure to realize asynchronous starting, the parallel computation of the down a plurality of transient stabilities calculating of many calculating nuclears, can satisfy unattended full-automatic high-level efficiency computational analysis; The batch computation schema of dynamic auto load balance is to consider to coordinate with other non-calculation task, satisfy and calculate the comfort and the optionally automatic operation mode of working, the automatic conversion of under this mode, can realize calculating dormancy, cascade being calculated in batches and parallel batch calculates.In physical fault operation computation process, can not interrupt under the whole calculation process condition of whole operations by manual intervention, realize that three kinds of computation schemas freely switch.
Description of drawings
Fig. 1: the asynchronous parallel computing method performing step block diagram that the present invention proposes based on the speed-raising electrical network transient stability counting yield of unit multinuclear;
Fig. 2: unit serial batch processing transient stability calculates subprocess step block diagram;
Fig. 3: unit multi-core parallel concurrent batch processing transient stability calculates subprocess step block diagram;
Fig. 4: the batch processing transient stability of unit multinuclear bringing onto load balance calculates subprocess step block diagram;
Fig. 5: transient stability calculates three kinds of computation schema principle schematic in batches, mainly shows and calculates the calculation processes of checking fault operation and fault item.
Embodiment
The block diagram of the specific embodiment of the invention such as Fig. 2 comprise the steps: based on the asynchronous parallel computing method of the speed-raising electrical network transient stability counting yield of unit multinuclear
(1) according to the electrical network analysis calculation requirement, from electrical network network trend topological model data, comprise generator, bus, alternating current circuit, transformer, transient stability device, DC equipment etc. according to the grid equipment composition, form the fault collection that the electrical network transient stability analysis requires.Realization is calculated based on the trend of electrical network steady-state model, forms the ground state method of operation that grid simulation calculates.Build on the electrical network transient stability model parameter that ground state trend level adapts, comprise that mainly generator prime machine model, wind-powered electricity generation unit and power electronic equipment models such as control system model, field regulator model, speed regulator model, DC control system model, power system stabilizer, PSS model, static state and dynamic load model, SVC/STATCOM/FLC, electrical network zero sequence comprise transformer zero sequence and circuit zero sequence model etc. in the electrical network.And needs research electrical network is carried out multiple ground state mode prepare, comprise withered big, withered little, rich big, rich for a short time etc., guarantee the convergence of ground state trend mode; According to analysis purpose for example the difference of N-1 fault scanning, the screening of N-2 fault etc. set up different classes of failure collection, and adjust electrical network transient stability model.
(2) determine ground state trend, fault collection and transient stability model corresponding relation, make up transient stability computational tasks task groups.With the TU task unit form, form the job patch task groups of whole evaluating objects, i.e. job task batch queue data.The job batch formation is carried out by the serial sequence pattern, and it is right therefore should to form transient stability compute associations data successively by the priority level of calculating.The operation of this step subitem comprises:
1. select to calculate convergent ground state mode data;
2. select to be used for the failure collection data that transient stability calculates;
3. select to be matched with the transient stability data of ground state trend mode;
4. will be 1.-data message that 3. obtains is integrated into an independently computational tasks T
i, wherein T represents operation, and the sequence number that on behalf of this subjob, i produce must add up successively.
5. repeat 1.-4. process, set up a plurality of independent computational tasks,, arrange the sequencing of calculation task, set up a complete batch jobs set of tasks according to calculating priority level
Whole computational tasks numbers that the N representative is provided with.In single homework, being provided with the job priority weight factor is P=[p
1, p
2...., p
N]
T, weight factor P
iShown job priority, for greater than zero natural number.When the operation sequence is carried out, by weight factor task priority is sorted, form and calculate priority sequence.For complete operation, with weight factor zero setting, when new operation begins to carry out, to press weight factor again to task ranking, cause operation weight factor as calculated is 0, is placed to the descending end, if new operation weight factor is 0, show that all computational tasks finish, calculate otherwise just start this operation.In order to preserve weight factor, set up deque's administrative power repeated factor, the weight factor that one of them queue management is provided with, the variation of the weight factor when another is used for Management Calculation for using next time.
Step 5. in according to dynamic change priority weight system, realize the serializing cascade dynamic execution of computing operation in batches.It carries out weight by the [p of deque
1, p
2...., p
N]
TColumn vector control, static reaction is manual sets control, in dynamic response execution and the computation process to the not dynamic execution sequence change of computational tasks.Performance coeffcient p to complete transient stability calculating
iZero setting to unenforced operation, when determining operation next time, realizes that by the column vector descending bubbling of computational tasks is selected.
(3) starting the operation of single group transient stability calculates.Before calculating, consider to rejudge the integrality and the validity of trend ground state mode, fault collection file and transient Model data, because in grid simulation computation process, data modeling and modification, mode calculate and the setting of dynamic model parameter adjustment, concrete fault data all may be changed by manual, have only the modification of considering data to change, just may satisfy up-to-date transient stability computational analysis demand, therefore starting the calculating previous moment, the integrality and the validity of data are being judged reliability and the high efficiency that helps whole parallel batch.The integrality judgement comprises the existence of ground state data existence, fault collection data, three judgements of existence of transient Model data; Validity is judged rationality and the validity judgement that comprises single failure in the failure collection; To check also in addition whether the ground state trend restrains, whether the setting of stable model parameter is effective, automatically judge whether the result of calculation judgement of stability exists and effectively, the stable definite value of judgement is set manually if there is no, just get default judgement and stablize definite value in conjunction with expertise.
(4) at current operation fault collection data, it is decomposed into single fault data formation.During decomposition, by the format string operation, in fault collection data, each failed row begins by spaced rows, comprise between-failures indications, calculation of fault title and three data fields of fail result data name in the spaced rows, wherein fault title and fail result data name territory allow not exist, and this method can produce automatically according to the title formation rule.
(5) implementing selected fault collection transient stability according to the computation schema of appointment calculates in batches.This method is divided into three types with computation schema: basic cascade serial in batches account form, band specify check figure the parallel batch computation schema, take the batch computation schema of dynamic auto load balance.Three kinds of computation schemas can switch in computation process at any time.Divide three branches to carry out single failure set transient stability according to mode designated and calculate, step is respectively:
Basic cascade serial is account form in batches.The cascade pattern is that the transient stability fault mode is calculated in a kind of serial automatically.It takes the concentrated fault one by one of fault apart, and fault data, transient stability model data, mode data are constantly matched, and passes to transient stability and calculates nuclear, finishes transient stability calculating by calculating nuclear.The subprocess step of this mode is as follows:
(a) from the fault data formation that the decomposition of (4) step obtains, get one group of single failure data, form the computational data group with mode data, the pairing of transient stability model data.
(b) band (a) computational data group starts transient stability calculating nuclear, and the transient stability of beginning single failure form calculates.
(c) by monitoring that constantly transient stability calculates progress, judge to calculate and whether finish.If finish, just jump to (e).If do not finish, just wait for a judgement cycle.The judgement cycle shows as time-delay T
d, implied value T
dIt is 1 second.The appointment of delay time directly influences the efficient of cascaded computation, if be provided with long, so last judge calculate the time waiting that finishes in theory the limit near T
d, because if bide one's time entering at the beginning etc., calculate and just finish, the time of returning to judge just will postpone a T once more
dIf the consideration worst case, with an average operation N fault, M operation altogether, the limit time because of repeating query inspection waste is so: N * M * T
d, unit is second.Therefore, time-delay T
dN * M second is just amplified in every increase by 1 second, whole cascaded computation time-delay waste.The time-delay definite value can not be too short, because the too short frequency that will cause judging increases, also may cause counting yield to reduce.
Arrange each transient stability and calculate consuming time equating.Have following expression formula to set up:
(N
0+1)×T
0＝(T
0+ΔT
ε)×N
ε
T
0Be the initial time delay definite value of acquiescence, N
0For evenly judging the number of times that the single transient stability calculates to be needed (supposition is judged after being uniformly, is not subjected to the machine working conditions change), consider the consumed time of judgement own so, minimum judgement single transient stability calculates the number of times that finishes and is at least (N
0+ 1).If the variable quantity Δ T of time-delay definite value
εExpression, and cause judging that because of the time-delay definite value changes total degree is N
ε, because of changing the judgement number of times that causes, the time-delay definite value is changed to so:
Consider discrepancy process expense switching time Δ T
ω, so, as delay time shortening, i.e. Δ T
εLess than 0 o'clock, Δ N was on the occasion of, available absolute value representation, had following expression formula to set up because of delay value shortens the judgement increase consuming time that causes:
Because of time of the deterministic process of coming in and going out generally delicate to a millisecond one-level, therefore Δ T is arranged
ω＜＜T
0So, reduce T to a certain extent
0Can improve cascaded computation efficient, if but setting is too little, causes
Situation will be opposite, and it is also inadvisable that therefore too little delay judgement time value is set.
In single transient stability computation process, whether finish for judging transient stability computation process, the empirical value of suggestion detect cycle is set at 1 second, prevents the long or too short influence to whole computing time that causes.
(d) each judgement detection time is added up, and compare with the maximum transient stability Calculation threshold values of setting consuming time, if surpass predefined threshold values, just stop transient stability computation process by force, and by force operation process recording in the calculation of fault result.Maximum transient stability Calculation threshold values consuming time is set, and is in order to prevent the ossified of transient stability computation program, may have problems because of finding the solution matrix, data integration, coupling, statistics, will cause transient stability to calculate deadlock, makes that cascaded computation can't continue to carry out.Because transient stability is set calculates ossified threshold values, efficiently solve the whole probability of failure of cascaded computation technically, improved whole computation's reliability.
(e) after transient stability calculates end, need to extract the result and carry out the judgement of stability analysis, the judgement of stability conclusion comprises stable and the instability differentiation.Unstable character loses stable kind according to electrical network and comprises classifications such as angle of attack unstability, dynamic buckling, voltage unstability, frequency unstability again.Stabilization result is differentiated and to be finished, and shows that primary fault calculates to finish, and related conclusions is placed in the results list, withdraws from this calculation of fault.
(f) extract fault data queuing data item,, enter next fault item computation process, get back to (a) step if also there is fault data.If all faults in the current computational tasks content have been calculated finish, gather all results in this subjob, form statistical abstract, withdraw from this subjob then and calculate.
Band is specified the parallel batch computation schema of check figure.This method is supported the concurrent calculating of multi-core computer transient stability, before realizing parallel computation, specify to be used for the CPU check figure that transient stability calculates in computing machine allows maximum check figure scope.Calculating is only relevant with the check figure of computing machine, and is irrelevant with the type and the CPU number of computing machine, supports single personal desktop machine, notebook and server.Must specify more than one check figure, just start parallel computation in batches theoretically.Its subprocess step is as follows:
1) judges that multinuclear existence and setting value are greater than 1.If existence and setting value directly jump to step 2 greater than 1); Otherwise, enter basic cascade serial computation process in batches.
2) set up the data queue that equates with it according to the check figure of appointment.In parallel batch computation process, be followed successively by data importing derivation, transient stability calculating, the decomposition of fault collection, transient stability interpretation of result, the preparation of transient stability computational data by ordering consuming time.Wait for and raising overall calculation efficient for avoiding transient stability to calculate, must be to have computing machine to examine when vacant one, calculate with regard to starting transient stability immediately, therefore always before the startup transient stability calculates, the whole calculating parameter data that go up requirement of filling of all parallel data formations.Hinder data relative transient stability of setup time for some reason and calculate and will lack, therefore data computation can not occur examines vacant appearance, and data also do not have complete ready situation.Under the normal operating condition, after parallel data formation always was ready, the wait transient stability calculates to be finished.
3) judge whether there is idle queues in the computational data formation.If there is idle queues, from the fault data formation that the decomposition of (4) step obtains, get one group of single failure data immediately and be filled in the idle data formation, specify to be retained in the formation before the transient stability computation program application always.If there is no idle data formation, the task that just enters wait for, wait for that any one has started transient stability computation process and finish after, discharge data queue.
4) judge that whether all formations have filled up data, if fill up, just change the 5th over to) step; Otherwise change the 9th over to) step.
5) judge whether that the calculating nuclear of appointment is all busy,, from data queue, get one group of parameter immediately, on the nuclear of free time, start transient stability computation process, carry out the transient stability calculating of the fault data of getting if vacant calculating nuclear is arranged.Calculate as long as start transient stability, will set up transient stability and calculate monitor thread, so that pass through polling mode with certain detection time, done state is calculated in real-time detection.
6) in a detect cycle, calculate end when detecting transient stability, just start resolving as a result immediately, implement power grid accident result of calculation and carry out classification analysis, comprise stable the qualitative judgement, the unstable classification judge, and minimum voltage point statistical study, the statistical study of merit angle variation amplitude peak etc.Interpretation of result is finished, and shows finishing of a calculation of fault, withdraws from this calculation of fault, jumps to 8).If in the detect cycle, do not detect any transient stability and calculate end, the record transient stability is calculated the duration consuming time that experiences.
7) judge that all parallel transient stabilities calculate the ossified qualification threshold values whether duration consuming time surpasses appointment,, just stop surpassing the transient stability computation process of threshold values by force, jump to 8) if surpass ossified qualification threshold values.Do not limit threshold values if surpass, just enter next detect cycle, jump to 6).When each transient stability calculates startup, all open a counter T
i, i represents the transient stability computation process sequence number that starts.There is a T in each computation process
i=K * Δ T, each computation process all has the K value of oneself, have only after the startup, K just begin the counting, when K * Δ T greater than T
Setting, T
SettingBe ossified threshold values, it is overtime to illustrate that transient stability calculates, and will start and calculate termination procedure.There is requirement in this method to ossified threshold values is set, not getting rid of transient stability calculating normally carries out but overlength consuming time, therefore may exist to extract and stop the normal possibility of calculating, therefore in time that is provided with during threshold values, can carry out a transient stability to the research electrical network and calculate tentative calculation, calculate the time of calculating, increase certain nargin then, as ossified qualification threshold values.Generally speaking, desirable 5 to 10 minutes of time margin, ossified qualification threshold values is generally less than 30 minutes.Ossified qualification threshold values can not be set too little, if be provided with too smallly, will cause a large amount of transient stabilities to calculate failure, is not that this method will be recommended, and this method fundamental purpose is to calculate a large amount of electric network faults for rapid batch.
8) no matter be that fault is calculated normal termination, or over time limit being forced to finished.Have only and finish to calculate, the just calculating nuclear and the data parameters formation of release busy immediately.
9) connect 4 on), in filling parallel data formation process, also have unfilled fault data if work as the prior fault operation, just fill it in the vacant data queue.If the fault data that will calculate does not just enter and judges that whether whole transient stabilities calculating finish, if finish to finish this subjob, otherwise wait for that all calculating is examined vacant.
The ALB computation schema.The ALB pattern is to determine whether starting transient stability computation process according to the computer dynamic loading condition, in transient stability computation process, may relate to other initiation of services consuming time, if this service priority is higher than computation schema, if lack the load balance dynamic mechanism, the mean allocation that will cause resource certainly will influence operator's use feeling and job schedule.The ALB pattern specifies the parallel batch computation schema of check figure very similar with band, also specifies the parallel batch computation schema step of check figure to carry out by band substantially on execution flow process and step, and main difference step is:
I. under the ALB pattern, do not judge in advance and specify check figure whether must as long as the appointment check figure is arranged,, be made as the minimum of computation requirement, will set up data queue yet greater than 1 even specify 1 nuclear.This pattern of that is to say does not carry out the 1st) the step judgement.
II. the 5th) in the step, add the dynamic load Rule of judgment, be to start transient stability computation process according to the load judgment result, this is than only considering that simply calculating a nuclear free time more scientifically utilizes the computing system resource.The principle that dynamic load is judged is as follows:
In the following formula, R
CpuNuclear occupancy, R are calculated in expression
MemBe memory usage, be two key indexs that ALB monitors, because in transient stability computation process, what influence computing velocity and efficient mainly is dominant frequency, the EMS memory occupation amount of calculating nuclear.T
PeriodExtremely short monitoring period when calculating for automatic load, T
IdleBe to calculate in the load monitoring period free time of nuclear, their ratio is the idleness of current calculating nuclear, after representing by normalization, is subtracted by 1 and to be calculated the nuclear load factor timely.M
TotalMemory amount when calculating for automatic load, M
UseEMS memory occupation total amount when calculating for automatic load, their ratio is exactly memory usage.R
C-settingTake the definite value upper limit, R for calculating nuclear
M-settingBe the EMS memory occupation definite value upper limit,, and have idle calculating nuclear, just start single transient stability and calculate if calculate nuclear occupancy and memory usage smaller or equal to its higher limit.General R
C-settingValue is 1/N (N is the total check figure of computing machine), R
M-settingGenerally get 0.8.The framework of above-mentioned definite value and computer CPU and type have necessarily gets in touch, and its definite value can suitably be adjusted according to concrete type.
(6) judge in the computational tasks that all whether fault datas all calculate and finish, finish that the single failure Calculation results is gathered, and the transient stability that forms whole fault collection calculates the comprehensive summary form, and preserves and be shown in the form if calculate.Do not finish if calculate, will continue (5) step process.
(7) judge whether whole transient stability computational tasks are finished, if do not have, it is return to (3) step just entering next operation calculating, otherwise, finish whole operations and calculate.
2, require the asynchronous parallel computing method of the speed-raising electrical network transient stability counting yield of the 1 unit multinuclear of being narrated according to the power book, it is characterized in that: realized three kinds of computation schemas that a plurality of faults of many group jobs are calculated in batches be basic cascade serial in batches account form, band specify check figure the parallel batch computation schema, take the batch computation schema of dynamic auto load balance.The basic cascade serial serial sequenceization of account form realization transient stability fault is in batches calculated in batches, is basic form of calculation; Band specifies the parallel batch computation schema of check figure to realize asynchronous starting, the parallel computation of the down a plurality of transient stabilities calculating of many calculating nuclears, can satisfy unattended full-automatic high-level efficiency computational analysis; The batch computation schema of dynamic auto load balance is to consider to coordinate with other non-calculation task, satisfy and calculate the comfort and the optionally automatic operation mode of working, the automatic conversion of under this mode, can realize calculating dormancy, cascade being calculated in batches and parallel batch calculates.In physical fault operation computation process, can not interrupt under the whole calculation process condition of whole operations by manual intervention, realize that three kinds of computation schemas freely switch.
The present invention obtains the exploitation checking in the grid simulation analysis platform of new generation based on the PSD-BPA kernel, the independence that forms is set up and has been integrated in the grid simulation analysis platform, obtains widely applying in certain electrical network extra-high voltage safety and stability adaptability research, certain interconnected big electric network security in-depth research, certain wind-powered electricity generation factory atmosphere electricity connecting system and operation Control Study, electricity net safety stable control strategy Research of Calculation and Analysis.For checking its function and speed-raising efficient, randomly draw a working computer and carry out transient stability calculating test.Its process is as follows:
One. the picked at random working computer.
Dell OPTIPLEX 960 desktop computers are chosen in test, and operating system adopts 32
Basic (EM) 32-Bit, processor are two nuclears that calculate
2Quad CPU Q9400, predominant frequency is 2.66GHz, saves as 2.00GB in the installation, hard disk 450G.
Two. determine test cases.
Choose certain regional power grid model data, total 2332 bus nodes of node have circuit and prop up 1731 of ways, have transformer and prop up 1738 of ways.Prepare the flow data file, stablize data file, 4 fault collection files comprise 10 faults, 50 faults, 100 faults and 300 faults respectively.
Three. test and logging test results.
Test divides two periods to carry out, and for each period, adopts identical test mode.Promptly undertaken by the batch processing of single task serial batch processing, the parallel batch processing of double-core, bringing onto load balance respectively, and the time of ending is opened in operation calculating carry out record.When carrying out the on-load balance test, by normal working method operation, for example write report on computers, browse data, receiving and dispatching mail etc.Test result sees the following form:
Table 1 test duration 1 test result
Table 2 test durations 2 test result
Four. the analytical test result.
By to the test result data analysis, the inventive method can the satisfying magnanimity electric network fault calculation requirement in batches, when implementing the electrical network transient stability and calculate, can obtain obvious speed-raising effect, test process and data show:
Can realize that unmanned transient stability parallel batch calculates;
Further improve the transient stability computing velocity, can under computing hardware framework support situation, realize that by increasing the computing machine check figure parallel batch calculates.
Can realize carrying out synchronously of transient stability calculating and statistical study as a result;
The parallel batch transient stability calculates can promote calculation of fault speed significantly, and batch calculating that walks abreast is criticized than serial and calculated the speed-raising ratio more than or equal to 3.1 times;
Parallel batch processing speed-raising good stability, in the test of twice repetition example, same fault collection is consuming time approaching;
The batch of on-load balance calculates, its speed-raising is subjected to the restriction of computer operation operating mode, when heavy computational burden, influential to whole counting yield, but by the test cases analysis, even under heavy load amount situation, also can make full use of computer resource, improve the transient stability computing velocity, can adapt to the batch that transient stability calculates under the various working and calculate automatically.
Claims (9)
1. based on unit multinuclear pattern speed-raising electrical network transient stability asynchronous parallel computing method, it is characterized in that described method comprises the steps:
(1) determines the ground state method of operation that grid simulation calculates, set up the electrical network transient stability model parameter that ground state trend level adapts, and needs research electrical network is carried out the ground state mode prepare; Set up different classes of failure collection according to the difference of analysis purpose, and adjust electrical network transient stability model;
(2) determine ground state trend, fault collection and transient stability model corresponding relation, make up transient stability computational tasks task groups, with the TU task unit form, form the job patch task groups of whole evaluating objects, be job task batch queue data, the job batch formation is carried out by the serial sequence pattern, and it is right promptly to form transient stability compute associations data successively by the priority level of calculating;
(3) integrality of judgment data and validity are opened the transient stability operation and are calculated;
(4) be single fault data formation with current operation fault collection data decomposition;
(5) implementing selected fault collection transient stability according to the computation schema of appointment calculates in batches;
(6) transient stability that forms whole fault collection calculates the comprehensive summary form;
(7) finishing whole operations after all the transient stability computational tasks is finished calculates.
2. computing method as claimed in claim 1 is characterized in that, described integrality judgement comprises the existence of ground state data existence, fault collection data, the existence of transient Model data;
Described validity is judged rationality and the validity judgement that comprises single failure in the failure collection.
3. computing method as claimed in claim 1 is characterized in that, described computation schema is divided into:
Basic cascade serial is account form in batches;
Band is specified the parallel batch computation schema of check figure;
Dynamic auto load balance computation schema.
4. computing method as claimed in claim 3 is characterized in that, the step of described basic cascade serial batch account form is as follows:
(a) from the fault data formation that described step (4) decomposition obtains, get one group of single failure data, form the computational data group with mode data, the pairing of transient stability model data;
(b) start transient stability with (a) computational data group and calculate nuclear;
(c) transient stability jumps to (e) after calculating and finishing;
(d) each judgement detect cycle is added up, and compares, surpass predefined threshold values, just stop transient stability computation process by force with the maximum transient stability Calculation threshold values of setting consuming time, and by force operation process recording in the calculation of fault result;
(e) the extraction result carries out the judgement of stability analysis, differentiates to be stabilization result, and promptly primary fault calculates and finishes, and related conclusions is placed in the results list, withdraws from this calculation of fault;
(f) extract fault data queuing data item, after the fault in the current computational tasks content has been calculated and finished, gather all result of calculations in this subjob, form statistical abstract, and withdraw from this subjob and calculate.
5. computing method as claimed in claim 3 is characterized in that, described parallel batch computation schema step is as follows:
1) multinuclear existence and setting value jump to 2 greater than 1);
2) set up the data queue that equates with it according to the check figure of appointment;
When 3) having idle queues in the computational data formation, from the fault data formation that described step (4) decomposition obtains, get one group of single failure data and be filled in the idle data formation, specify to be retained in the formation before the transient stability computation program application always;
4) data have been filled up in all formations, change the 5th over to) step;
5) calculating of appointment nuclear has vacant calculating nuclear, gets one group of parameter from data queue, starts transient stability computation process on the nuclear of free time, and the transient stability that carries out the fault data of getting calculates;
6) in detect cycle, calculate end when detecting transient stability, just start resolving as a result immediately, to implement power grid accident result of calculation and carry out classification analysis, interpretation of result is finished, promptly finishing of a calculation of fault withdrawed from this calculation of fault, jumps to 8);
7) all parallel transient stabilities calculate the ossified qualification threshold values that duration consuming time surpasses appointment, stop surpassing the transient stability computation process of threshold values by force, jump to 8);
8) finish to calculate the calculating nuclear of release busy immediately and data parameters formation;
9) fill in the parallel data formation process,, just fill it in the vacant data queue when the prior fault operation also has unfilled fault data;
6. computing method as claimed in claim 3 is characterized in that, described ALB calculation procedure is:
I. set up the data queue that equates with it according to the check figure of appointment;
When II. having idle queues in the computational data formation, from the fault data formation that described step (4) decomposition obtains, get one group of single failure data and be filled in the idle data formation, specify to be retained in the formation before the transient stability computation program application always;
III. data have been filled up in all formations, change the V step over to;
IV. the calculating of appointment nuclear has vacant calculating nuclear, gets one group of parameter from data queue, starts transient stability computation process according to the load judgment result on the nuclear of free time, and the transient stability that carries out the fault data of getting calculates, and the dynamic load judgment formula is as follows:
7. computing method as claimed in claim 1 is characterized in that, the substep of described step (2) comprising:
1. select to calculate convergent ground state mode data;
2. select to be used for the failure collection data that transient stability calculates;
3. select to be matched with the transient stability data of ground state trend mode;
4. will be 1.-data message that 3. obtains is integrated into an independently computational tasks T
i
5. repeat 1.-4. process, set up independent computational tasks more than two, P sorts to task priority by weight factor, forms to calculate priority sequence, according to calculating priority level, arranges the sequencing of calculation task; For complete operation, with weight factor P zero setting, when new operation begins to carry out, will be again by weight factor P to task ranking, because of as calculated operation weight factor P is 0, be placed to the descending end;
Set up deque's administrative power repeated factor, the weight factor that one of them queue management is provided with, the variation of the weight factor when another is used for Management Calculation.
8. computing method as claimed in claim 7 is characterized in that, described step 5. in according to dynamic change priority weight system, realize the serializing cascade dynamic execution of computing operation in batches, it carries out weight by the [p of deque
1, p
2...., p
N]
TColumn vector control, static reaction is manual sets control, in dynamic response execution and the computation process to the not dynamic execution sequence change of computational tasks, the performance coeffcient p that complete transient stability is calculated
iZero setting to unenforced operation, when determining operation next time, realizes that by the column vector descending bubbling of computational tasks is selected.
9. computing method as claimed in claim 4 is characterized in that, the transient stability that provides in the described step (c) calculates and finishes to judge that the time set of single detect cycle is 1 second.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110004612.4A CN102073766B (en) | 2011-01-11 | 2011-01-11 | Method for accelerating asynchronous parallel computation of transient stability of power grid based on single-computer multiple-core mode |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110004612.4A CN102073766B (en) | 2011-01-11 | 2011-01-11 | Method for accelerating asynchronous parallel computation of transient stability of power grid based on single-computer multiple-core mode |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102073766A true CN102073766A (en) | 2011-05-25 |
| CN102073766B CN102073766B (en) | 2014-03-12 |
Family
ID=44032304
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110004612.4A Active CN102073766B (en) | 2011-01-11 | 2011-01-11 | Method for accelerating asynchronous parallel computation of transient stability of power grid based on single-computer multiple-core mode |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102073766B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102332709A (en) * | 2011-09-13 | 2012-01-25 | 中国电力科学研究院 | Equivalent high-voltage direct-current transmission direct-current protection digital simulation method and system |
| CN102629220A (en) * | 2012-03-08 | 2012-08-08 | 北京神州数码思特奇信息技术股份有限公司 | Dynamic task allocation and management method |
| CN102709953A (en) * | 2012-05-17 | 2012-10-03 | 中国电力科学研究院 | Online quantization evaluation method of transient state stability of power grid based on WAMS (wide area measurement system) and unit pair |
| CN103198188A (en) * | 2013-04-03 | 2013-07-10 | 国家电网公司 | Static synchronous compensator (STATCOM) modeling method in power system analysis software package |
| CN103365715A (en) * | 2013-06-17 | 2013-10-23 | 广东电网公司电力调度控制中心 | Method and system for quickly calling BPA transient calculation program |
| CN104899397A (en) * | 2015-06-19 | 2015-09-09 | 国家电网公司 | Fault simulation result processing method and device |
| CN109344117A (en) * | 2018-10-10 | 2019-02-15 | 四川新网银行股份有限公司 | A kind of risk detecting system based on concurrent |
| CN110348066A (en) * | 2019-06-17 | 2019-10-18 | 国网西藏电力有限公司电力科学研究院 | Access the protective relaying device check method of weak grid on a large scale for new energy |
| CN118094841A (en) * | 2024-04-25 | 2024-05-28 | 中国电力科学研究院有限公司 | Method and device for carrying out multi-tide parallel calculation on power grid mode |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101169743A (en) * | 2007-11-27 | 2008-04-30 | 南京大学 | Method for implementing parallel power flow calculation based on multi-core computer in electric grid |
-
2011
- 2011-01-11 CN CN201110004612.4A patent/CN102073766B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101169743A (en) * | 2007-11-27 | 2008-04-30 | 南京大学 | Method for implementing parallel power flow calculation based on multi-core computer in electric grid |
Non-Patent Citations (2)
| Title |
|---|
| 宋墩文 等: "实时动态安全评估中在线和离线仿真互访过程研究", 《电网技术》, vol. 34, no. 2, 28 February 2010 (2010-02-28) * |
| 郑超 等: "在线动态安全评估与预警系统的功能设计与实现", 《电网技术》, vol. 34, no. 3, 31 March 2010 (2010-03-31) * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102332709B (en) * | 2011-09-13 | 2014-10-01 | 中国电力科学研究院 | Equivalent high-voltage direct-current transmission direct-current protection digital simulation method and system |
| CN102332709A (en) * | 2011-09-13 | 2012-01-25 | 中国电力科学研究院 | Equivalent high-voltage direct-current transmission direct-current protection digital simulation method and system |
| CN102629220A (en) * | 2012-03-08 | 2012-08-08 | 北京神州数码思特奇信息技术股份有限公司 | Dynamic task allocation and management method |
| CN102709953B (en) * | 2012-05-17 | 2016-05-25 | 中国电力科学研究院 | A kind of based on WAMS and the right online quantitative estimation method of Power Network Transient Stability of unit |
| CN102709953A (en) * | 2012-05-17 | 2012-10-03 | 中国电力科学研究院 | Online quantization evaluation method of transient state stability of power grid based on WAMS (wide area measurement system) and unit pair |
| CN103198188A (en) * | 2013-04-03 | 2013-07-10 | 国家电网公司 | Static synchronous compensator (STATCOM) modeling method in power system analysis software package |
| CN103365715A (en) * | 2013-06-17 | 2013-10-23 | 广东电网公司电力调度控制中心 | Method and system for quickly calling BPA transient calculation program |
| CN104899397A (en) * | 2015-06-19 | 2015-09-09 | 国家电网公司 | Fault simulation result processing method and device |
| CN104899397B (en) * | 2015-06-19 | 2018-04-06 | 国家电网公司 | The treating method and apparatus of fault simulation result |
| CN109344117A (en) * | 2018-10-10 | 2019-02-15 | 四川新网银行股份有限公司 | A kind of risk detecting system based on concurrent |
| CN110348066A (en) * | 2019-06-17 | 2019-10-18 | 国网西藏电力有限公司电力科学研究院 | Access the protective relaying device check method of weak grid on a large scale for new energy |
| CN118094841A (en) * | 2024-04-25 | 2024-05-28 | 中国电力科学研究院有限公司 | Method and device for carrying out multi-tide parallel calculation on power grid mode |
| CN118094841B (en) * | 2024-04-25 | 2024-07-09 | 中国电力科学研究院有限公司 | Method and device for carrying out multi-tide parallel calculation on power grid mode |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102073766B (en) | 2014-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102073766B (en) | Method for accelerating asynchronous parallel computation of transient stability of power grid based on single-computer multiple-core mode | |
| Irisarri et al. | Economic dispatch with network and ramping constraints via interior point methods | |
| CN102063336B (en) | Distributed computing multiple application function asynchronous concurrent scheduling method | |
| CN101630273B (en) | Small interference stability simulation method of electric power system | |
| CN101447671B (en) | Automatic integrating and adjusting method for flow data | |
| CN107861606A (en) | A kind of heterogeneous polynuclear power cap method by coordinating DVFS and duty mapping | |
| CN108599151A (en) | A kind of Model in Reliability Evaluation of Power Systems dynamic parallel computational methods | |
| CN109685380A (en) | A kind of large complicated comprehensive performance evaluation method for repairing device | |
| CN106410789A (en) | Power system black-start portioning method for isolated network | |
| Huang | Modeling of software partition for distributed real-time applications | |
| CN113193881A (en) | Intelligent functional verification detection method based on HPLC (high performance liquid chromatography) advanced application | |
| CN102945198A (en) | Method for characterizing application characteristics of high performance computing | |
| CN106056305A (en) | Power generation system reliability rapid assessment method based on state clustering | |
| Wang et al. | PV hosting capacity estimation: Experiences with scalable framework | |
| CN117913857A (en) | Thermal power generating unit frequency modulation and peak shaving method and system | |
| CN101477168B (en) | Parallelization test system and method for transient stability of electric power system | |
| CN105391597A (en) | Multi-server-simulation-based avalanche testing method, apparatus, and system | |
| CN104578050A (en) | Transient stability strongly-correlated power transmission section identification method for power grid | |
| CN104578048B (en) | Method for quickly evaluating transient security and stability of enumerated combined successive fault set on line | |
| CN103973793A (en) | Distributed aviation electronic testing system | |
| CN101364286B (en) | Prediction method for equivalent forced outage rate of thermoelectric generating set | |
| CN113517693B (en) | Weak voltage transient stability node analysis method and device | |
| CN103577685A (en) | Electric system running state online and offline estimation hybrid scheduling system | |
| Pohl et al. | Integrating an autonomous agent-based power flow control system into control center software | |
| CN103236696B (en) | Planning load flow calculation-aided system of power transmission network |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| ASS | Succession or assignment of patent right |
Owner name: STATE ELECTRIC NET CROP. Effective date: 20130424 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20130424 Address after: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15 Applicant after: China Electric Power Research Institute Applicant after: State Grid Corporation of China Address before: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15 Applicant before: China Electric Power Research Institute |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |