CN109818346A - Alloy White Iron Method for Transient Calculation and device based on frequency-domain analysis - Google Patents
Alloy White Iron Method for Transient Calculation and device based on frequency-domain analysis Download PDFInfo
- Publication number
- CN109818346A CN109818346A CN201711173997.0A CN201711173997A CN109818346A CN 109818346 A CN109818346 A CN 109818346A CN 201711173997 A CN201711173997 A CN 201711173997A CN 109818346 A CN109818346 A CN 109818346A
- Authority
- CN
- China
- Prior art keywords
- pole
- frequency
- public pole
- matrix
- fdne
- 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.)
- Pending
Links
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The invention discloses a kind of Alloy White Iron Method for Transient Calculation and device based on frequency-domain analysis, wherein method includes: to obtain N2One group of identical initial public pole in a element；Vector fitting is carried out one by one to each element of FDNE matrix according to initial public pole, pole and its number to guarantee each element are identical；If pole is not converged, one group of newly public pole is obtained after vector fitting；Vector fitting is re-started using new public pole as initial public pole, until the difference of initial public pole and new public pole in predetermined threshold level, obtains the FDNE matrix with public pole；Solve cyclization transient current.This method can realize Alloy White Iron Transient calculation by the rapid vector fitting process based on frequency-domain analysis, improve and calculate accuracy and speed, and further increase power supply reliability, ensure the normal production and living of user.
Description
Technical field
The present invention relates to closed loop network current temporary state computing technique fields in electric system, in particular to a kind of based on frequency
The Alloy White Iron Method for Transient Calculation and device of domain analysis.
Background technique
Cyclization transient current is the principal element of cyclization success, therefore, quickly and accurately calculates cyclization transient current
Be related to can cyclization success, and for ensure the normal production and living of user, improve power supply reliability be of great significance.
In the related technology, the simplified Thevenin's equivalence circuit counting Alloy White Iron of some model uses, but actuarial precision
It is insufficient；By the whole network EM transient calculation Alloy White Iron, although precision meets, because calculation amount cannot achieve greatly very much；According to
Micro-judgment whether can cyclization, can only be directed to specific network, do not have generality and representativeness, it would be highly desirable to solve.
Summary of the invention
The present invention is directed to solve at least some of the technical problems in related technologies.
For this purpose, an object of the present invention is to provide a kind of Alloy White Iron Method for Transient Calculation based on frequency-domain analysis,
This method, which can be improved, calculates accuracy and speed, and further increases power supply reliability, ensures the normal production of user
Life.
It is another object of the present invention to propose a kind of Alloy White Iron Transient calculation device based on frequency-domain analysis.
In order to achieve the above objectives, one aspect of the present invention embodiment proposes a kind of Alloy White Iron transient state based on frequency-domain analysis
Calculation method, comprising the following steps: obtain N2One group of identical initial public pole in a element；According to the initial public pole
Point carries out vector fitting to each element of FDNE matrix one by one, and pole and its number to guarantee each element are identical；
If pole is not converged, one group of newly public pole is obtained after vector fitting；Using the new public pole as described in
Initial public pole re-starts vector fitting, until the difference of the initial public pole and the new public pole is default
In threshold value, the FDNE matrix with public pole is obtained；Solve cyclization transient current.
The Alloy White Iron Method for Transient Calculation based on frequency-domain analysis of the embodiment of the present invention proposes according to embodiments of the present invention
The Alloy White Iron Method for Transient Calculation based on frequency-domain analysis, can be realized by the rapid vector fitting process based on frequency-domain analysis
Alloy White Iron Transient calculation overcomes the shortcomings that accuracy and speed can not be fully met in Alloy White Iron Transient calculation, effectively improves
Accuracy and speed are calculated, and further increases power supply reliability, ensures the normal production and living of user.
In addition, the Alloy White Iron Method for Transient Calculation according to the above embodiment of the present invention based on frequency-domain analysis can also have
There is following additional technical characteristic:
Further, in one embodiment of the invention, the acquisition N2In a element one group it is identical initial public
Pole further comprises: if sampling number is N0, then the sampling value matrix of N number of N × N-dimensional is sought；By N number of N × N-dimensional
Sampling value matrix in N in each sampling value matrix2A element value is successively added, to obtain N0× 1 sampled value vector；To institute
It states sampled value vector and carries out vector fitting, obtain rational function；The pole for obtaining the rational function obtains described initial public
Pole.
Further, in one embodiment of the invention, each element of the FDNE matrix with public pole
All it is a frequency-domain function:
Wherein, pole aiWith residual ciIt for real number or is complex conjugate pair, d and h are real number, and n is pole number, and s is multiple
Frequency.
Further, in one embodiment of the invention, the solution cyclization transient current further comprises:
Establish FDNE (Frequency Dependent Network Equivalent, based on frequency-domain analysis technology equivalence
Model, to solve cyclization transient current, wherein power transmission network voltage and current relationship can indicate are as follows:
Wherein, IBFor boundary bus current, IEFor other bus currents, U in addition to the bus of boundaryEFor in addition to the bus of boundary its
His busbar voltage, UBFor boundary busbar voltage, f is frequency, and subscript B is boundary bus, and subscript E is that other in power transmission network are female
Line；
To obtain being retracted to FDNE matrix and Nuo Dun equal currents behind boundary, and the promise equal currents are as follows:
Wherein, f0For fundamental frequency.
Further, in one embodiment of the invention, further includes: convert three-phase electricity for the promise equal currents
Stream source；Difference is carried out to the three-phase current source by implicit trapezoid method.
In order to achieve the above objectives, it is temporary to propose a kind of Alloy White Iron based on frequency-domain analysis for another aspect of the present invention embodiment
State computing device, comprising: module is obtained, for obtaining N2One group of identical initial public pole in a element；First fitting mould
Block, it is described every to guarantee for carrying out vector fitting one by one to each element of FDNE matrix according to the initial public pole
The pole and its number of a element are identical；If judgment module obtains one after vector fitting not converged for pole
The new public pole of group；Second fitting module, for re-starting arrow for the new public pole as the initial public pole
Amount fitting, until the difference of the initial public pole and the new public pole in predetermined threshold level, obtains having public
The FDNE matrix of pole；Computing module solves cyclization transient current.
The Alloy White Iron Transient calculation device based on frequency-domain analysis of the embodiment of the present invention proposes according to embodiments of the present invention
The Alloy White Iron Method for Transient Calculation based on frequency-domain analysis, can be realized by the rapid vector fitting process based on frequency-domain analysis
Alloy White Iron Transient calculation overcomes the shortcomings that accuracy and speed can not be fully met in Alloy White Iron Transient calculation, effectively improves
Accuracy and speed are calculated, and further increases power supply reliability, ensures the normal production and living of user.
In addition, the Alloy White Iron Transient calculation device according to the above embodiment of the present invention based on frequency-domain analysis can also have
There is following additional technical characteristic:
Further, in one embodiment of the invention, the acquisition module further comprises: first acquisition unit,
If being N for sampling number0, then the sampling value matrix of N number of N × N-dimensional is sought；Computing unit is used for N number of N × N-dimensional
Sampling value matrix in N in each sampling value matrix2A element value is successively added, to obtain N0× 1 sampled value vector；Fitting
Unit obtains rational function for carrying out vector fitting to the sampled value vector；Second acquisition unit, it is described for obtaining
The pole of rational function obtains the initial public pole.
Further, in one embodiment of the invention, each element of the FDNE matrix with public pole
All it is a frequency-domain function:
Wherein, pole aiWith residual ciIt for real number or is complex conjugate pair, d and h are real number, and n is pole number, and s is multifrequency
Rate.
Further, in one embodiment of the invention, the computing module further comprises: establishing unit, is used for
FDNE Equivalent Model is established, to solve cyclization transient current, wherein power transmission network voltage and current relationship can indicate are as follows:
Wherein, IBFor boundary bus current, IEFor other bus currents, U in addition to the bus of boundaryEFor in addition to the bus of boundary its
His busbar voltage, UBFor boundary busbar voltage, f is frequency, and subscript B is boundary bus, and subscript E is that other in power transmission network are female
Line；
To obtain being retracted to FDNE matrix and Nuo Dun equal currents behind boundary, and the promise equal currents are as follows:
Wherein, f0For fundamental frequency.
Further, in one embodiment of the invention, further includes: convert three-phase electricity for the promise equal currents
Stream source；Difference is carried out to the three-phase current source by implicit trapezoid method.
The additional aspect of the present invention and advantage will be set forth in part in the description, and will partially become from the following description
Obviously, or practice through the invention is recognized.
Detailed description of the invention
Above-mentioned and/or additional aspect and advantage of the invention will become from the following description of the accompanying drawings of embodiments
Obviously and it is readily appreciated that, in which:
Fig. 1 is the flow chart according to the Alloy White Iron Method for Transient Calculation based on frequency-domain analysis of the embodiment of the present invention；
Fig. 2 is the process according to the Alloy White Iron Method for Transient Calculation based on frequency-domain analysis of one embodiment of the invention
Figure；
Fig. 3 is the Equivalent Model schematic diagram based on FDNE according to one embodiment of the invention；
Fig. 4 is the structural representation according to the Alloy White Iron Transient calculation device based on frequency-domain analysis of the embodiment of the present invention
Figure.
Specific embodiment
The embodiment of the present invention is described below in detail, examples of the embodiments are shown in the accompanying drawings, wherein from beginning to end
Same or similar label indicates same or similar element or element with the same or similar functions.Below with reference to attached
The embodiment of figure description is exemplary, it is intended to is used to explain the present invention, and is not considered as limiting the invention.
The Alloy White Iron Transient calculation based on frequency-domain analysis proposed according to embodiments of the present invention is described with reference to the accompanying drawings
Method and device describes the Alloy White Iron transient state based on frequency-domain analysis proposed according to embodiments of the present invention with reference to the accompanying drawings first
Calculation method.
Fig. 1 is the flow chart of the Alloy White Iron Method for Transient Calculation based on frequency-domain analysis of the embodiment of the present invention.
As shown in Figure 1, should Alloy White Iron Method for Transient Calculation based on frequency-domain analysis the following steps are included:
In step s101, N is obtained2One group of identical initial public pole in a element.
Further, in one embodiment of the invention, N is obtained2One group of identical initial public pole in a element,
Further comprise: if sampling number is N0, then the sampling value matrix of N number of N × N-dimensional is sought；By the sampled value square of N number of N × N-dimensional
N in each sampling value matrix in battle array2A element value is successively added, to obtain N0× 1 sampled value vector；To sampled value vector into
Row vector fitting, obtains rational function；The pole for obtaining rational function obtains initial public pole.
It is understood that the embodiment of the present invention can first be this N combined with Figure 1 and Figure 2,2A element seek one group it is identical
Initial public pole.If sampling number is N0, then the sampling value matrix of N number of N × N-dimensional can be sought, and by each sampled value square
N in battle array2A element value is successively added, for N0A N then can be obtained in a sampling value matrix0× 1 " sampled value vector ", by this
A " sampled value vector " carries out vector fitting, and the pole of this rational function is set as original by an available rational function
The initial public pole of FDNE matrix.
In step s 102, vector fitting is carried out to each element of FDNE matrix according to initial public pole one by one, with
Guarantee that pole and its number of each element are identical.
It is understood that combined with Figure 1 and Figure 2, based on the initial public pole in step S101, to this FDNE matrix
In each element carry out vector fitting one by one, pole and its number to guarantee each element of FDNE matrix are identical.
In step s 103, if pole is not converged, one group of newly public pole is obtained after vector fitting.
In step S104, vector fitting is re-started using new public pole as initial public pole, until initial public
The difference of copolar point and new public pole obtains the FDNE matrix with public pole in predetermined threshold level.
It is understood that combined with Figure 1 and Figure 2, by vector fitting method it is found that when pole is not converged, vector fitting method
After can find out one group of new pole, and this can be organized to the initial public pole that new pole is re-used as in step S102
Point, and vector fitting is re-started, the difference of the new pole after initial common node and algorithm is in certain threshold value
It is interior, it can restrain under normal circumstances and acquire the FDNE matrix with public pole for iteration 3-4 times.
It is understood that the embodiment of the present invention during vector fitting, can take following practical processing side
Method:
First in practical finding process, can generally enable h is 0, for real network, in infinitely great frequency
Under, the node admittance of real network can not be infinity.
Therefore, in order to guarantee that accurate frequency characteristic is all presented in the rational function of fitting under fundamental frequency and wideband, generally exist
Relatively large weight is arranged to the sampled value under fundamental frequency and wideband when fitting.
Rapid vector approximating method can guarantee that the pole of Y (s) each element is identical, therefore only residual item and constant
Item is not identical, then N × N-dimensional FDNE matrix Y (s) may be expressed as:
Above-mentioned matrix is write as to the form of transmission function:
Y (s)=C (sE-A)-1B+D,
It is E unit matrix in formula, dimension is identical as matrix A.
A=diag (A1 … Ak …AN),
B=diag (B1 … Bk … BN),
Bk=[(1 1 ... 1)(1×n)]T,
K, m=1,2 ..., N.
In step s105, cyclization transient current is solved.
Further, in one embodiment of the invention, cyclization transient current is solved, further comprises: establishing FDNE
Equivalent Model, to solve cyclization transient current, wherein power transmission network voltage and current relationship can indicate are as follows:
Wherein, IBFor boundary bus current, IEFor other bus currents, U in addition to the bus of boundaryEFor in addition to the bus of boundary its
His busbar voltage, UBFor boundary busbar voltage, f is frequency, and subscript B is boundary bus, and subscript E is that other in power transmission network are female
Line；
To obtain being retracted to FDNE matrix and Nuo Dun equal currents behind boundary, and promise equal currents are as follows:
Wherein, f0For fundamental frequency.
It is understood that FDNE Equivalent Model is as shown in figure 3, the embodiment of the present invention can use rapid vector fitting process
Solve cyclization transient current.The embodiment of the present invention considers FDNE, therefore power transmission network voltage and current relationship can indicate are as follows:
Wherein, subscript B indicates boundary bus, and subscript E indicates other buses in power transmission network.Gauss is carried out to formula to disappear
It goes, eliminates UE, to obtain being retracted to the FDNE matrix and Nuo Dun equal currents behind boundary.It is defeated due to 220kV voltage class
Farther out, therefore promise equal currents can only consider fundamental component to power grid clutch circling point distance, therefore the equivalent formula of promise are as follows:
Wherein, f0It is fundamental frequency.
Further, in one embodiment of the invention, the method for the embodiment of the present invention further include: by the equivalent electricity of promise
Circulation turns to three-phase current source；Difference is carried out to three-phase current source by implicit trapezoid method.
It should be noted that in order to allow to need to do in the equivalent result application Alloy White Iron Transient calculation of above-mentioned promise
Two-step pretreatment: three-phase current source is converted by promise equal currents first；Secondly difference is carried out using implicit trapezoid method.Due to this
The pole and its number of each element are all identical in the FDNE matrix of inventive embodiments, therefore can unify difference.And according to
Simulation result calculates cyclization transient current peak value and transient characterisitics.
Further, in one embodiment of the invention, each element of the FDNE matrix with public pole is
One frequency-domain function:
Wherein, pole aiWith residual ciIt for real number or is complex conjugate pair, d and h are real number, and n is pole number, and s is multifrequency
Rate.
It is understood that the essence of FDNE is one using frequency as the node admittance matrix of function, it can indicate frequency
Node admittance matrix under rate situation of change.Therefore, N × N-dimensional FDNE can be expressed as a frequency domain matrix:
A series of sampled value of each element of FDNE matrix under frequencies is fitted to one continuously by vector fitting method
Rational function, and each element in the FDNE matrix found out by vector fitting method may be expressed as a frequency-domain function:
Wherein, pole aiWith residual ciIt is real number or respectively with complex conjugate to appearance；D and h is real number；N is pole
It counts, a of different elements in FDNE matrixi、ci, d and h be different.
It is understood that working as N2When the pole of a rational function is not identical, the N in FDNE matrix2A element is answered one by one
Time-domain-simulation is used, calculation amount can sharply increase, while time domain realizes also inconvenience, and rapid vector fitting process can be effectively
Reduce calculation amount.
The Alloy White Iron Method for Transient Calculation based on frequency-domain analysis proposed according to embodiments of the present invention, can be by being based on
The rapid vector fitting process of frequency-domain analysis realizes Alloy White Iron Transient calculation, overcomes accuracy and speed in Alloy White Iron Transient calculation
The shortcomings that can not fully meeting, effectively improves and calculates accuracy and speed, and further increases power supply reliability, ensures and uses
The normal production and living at family.
The Alloy White Iron Transient calculation based on frequency-domain analysis proposed according to embodiments of the present invention is described referring next to attached drawing
Device.
Fig. 4 is the structural schematic diagram of the Alloy White Iron Transient calculation device based on frequency-domain analysis of the embodiment of the present invention.
As shown in figure 4, being somebody's turn to do the Alloy White Iron Transient calculation device 10 based on frequency-domain analysis includes: to obtain module 100, first
Fitting module 200, judgment module 300, the second fitting module 400 and computing module 500.
Wherein, mould 100 is obtained for obtaining N2One group of identical initial public pole in a element.First fitting module
200 for carrying out vector fitting to each element of FDNE matrix according to initial public pole one by one, to guarantee each element
Pole and its number are identical.If judgment module 300 for pole it is not converged, obtained after vector fitting one group it is newly public
Copolar point.Second fitting module 400 is used to re-start vector fitting for new public pole as initial public pole, until just
The difference of the public pole that begins and new public pole obtains the FDNE matrix with public pole in predetermined threshold level.Calculate mould
Block 500 solves cyclization transient current.The device 10 of the embodiment of the present invention can be fitted by the rapid vector based on frequency-domain analysis
Method realizes Alloy White Iron Transient calculation, improves and calculates accuracy and speed, and further increases power supply reliability, ensures and uses
The normal production and living at family.
Further, in one embodiment of the invention, obtaining module 100 further comprises: first acquisition unit, meter
Calculate unit, fitting unit and second acquisition unit.
Wherein, if it is N that first acquisition unit, which is used for sampling number,0, then the sampling value matrix of N number of N × N-dimensional is sought.Meter
Unit is calculated to be used for N in sampling value matrix each in the sampling value matrix of N number of N × N-dimensional2A element value is successively added, to obtain N0
× 1 sampled value vector.Fitting unit is used to carry out vector fitting to sampled value vector, obtains rational function.Second obtains list
Member obtains initial public pole for obtaining the pole of rational function.
Further, in one embodiment of the invention, each element of the FDNE matrix with public pole is
One frequency-domain function:
Wherein, pole aiWith residual ciIt for real number or is complex conjugate pair, d and h are real number, and n is pole number, and s is multifrequency
Rate.
Further, in one embodiment of the invention, computing module 500 further comprises: establishing unit.Wherein,
Unit is established for establishing FDNE Equivalent Model, to solve cyclization transient current, wherein power transmission network voltage and current relationship can
To indicate are as follows:
Wherein, IBFor boundary bus current, IEFor other bus currents, U in addition to the bus of boundaryEFor in addition to the bus of boundary its
His busbar voltage, UBFor boundary busbar voltage, f is frequency, and subscript B is boundary bus, and subscript E is that other in power transmission network are female
Line；
To obtain being retracted to FDNE matrix and Nuo Dun equal currents behind boundary, and promise equal currents are as follows:
Wherein, f0For fundamental frequency.
Further, in one embodiment of the invention, the device 10 of the embodiment of the present invention further include: by promise equivalence
Electric current is converted into three-phase current source；Difference is carried out to three-phase current source by implicit trapezoid method.
It should be noted that the aforementioned explanation to the Alloy White Iron Method for Transient Calculation embodiment based on frequency-domain analysis
It is also applied for the Alloy White Iron Transient calculation device based on frequency-domain analysis of the embodiment, details are not described herein again.
The Alloy White Iron Transient calculation device based on frequency-domain analysis proposed according to embodiments of the present invention, can be by being based on
The rapid vector fitting process of frequency-domain analysis realizes Alloy White Iron Transient calculation, overcomes accuracy and speed in Alloy White Iron Transient calculation
The shortcomings that can not fully meeting, effectively improves and calculates accuracy and speed, and further increases power supply reliability, ensures and uses
The normal production and living at family.
In the description of the present invention, it is to be understood that, term " center ", " longitudinal direction ", " transverse direction ", " length ", " width ",
" thickness ", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom" "inner", "outside", " up time
The orientation or positional relationship of the instructions such as needle ", " counterclockwise ", " axial direction ", " radial direction ", " circumferential direction " be orientation based on the figure or
Positional relationship is merely for convenience of description of the present invention and simplification of the description, rather than the device or element of indication or suggestion meaning must
There must be specific orientation, be constructed and operated in a specific orientation, therefore be not considered as limiting the invention.
In addition, term " first ", " second " are used for descriptive purposes only and cannot be understood as indicating or suggesting relative importance
Or implicitly indicate the quantity of indicated technical characteristic.Define " first " as a result, the feature of " second " can be expressed or
Implicitly include at least one this feature.In the description of the present invention, the meaning of " plurality " is at least two, such as two, three
It is a etc., unless otherwise specifically defined.
In the present invention unless specifically defined or limited otherwise, term " installation ", " connected ", " connection ", " fixation " etc.
Term shall be understood in a broad sense, for example, it may be being fixedly connected, may be a detachable connection, or integral；It can be mechanical connect
It connects, is also possible to be electrically connected；It can be directly connected, can also can be in two elements indirectly connected through an intermediary
The interaction relationship of the connection in portion or two elements, unless otherwise restricted clearly.For those of ordinary skill in the art
For, the specific meanings of the above terms in the present invention can be understood according to specific conditions.
In the present invention unless specifically defined or limited otherwise, fisrt feature in the second feature " on " or " down " can be with
It is that the first and second features directly contact or the first and second features pass through intermediary mediate contact.Moreover, fisrt feature exists
Second feature " on ", " top " and " above " but fisrt feature be directly above or diagonally above the second feature, or be merely representative of
First feature horizontal height is higher than second feature.Fisrt feature can be under the second feature " below ", " below " and " below "
One feature is directly under or diagonally below the second feature, or is merely representative of first feature horizontal height less than second feature.
In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", " specifically show
The description of example " or " some examples " etc. means specific features, structure, material or spy described in conjunction with this embodiment or example
Point is included at least one embodiment or example of the invention.In the present specification, schematic expression of the above terms are not
It must be directed to identical embodiment or example.Moreover, particular features, structures, materials, or characteristics described can be in office
It can be combined in any suitable manner in one or more embodiment or examples.In addition, without conflicting with each other, the skill of this field
Art personnel can tie the feature of different embodiments or examples described in this specification and different embodiments or examples
It closes and combines.
Although the embodiments of the present invention has been shown and described above, it is to be understood that above-described embodiment is example
Property, it is not considered as limiting the invention, those skilled in the art within the scope of the invention can be to above-mentioned
Embodiment is changed, modifies, replacement and variant.
Claims (10)
1. a kind of Alloy White Iron Method for Transient Calculation based on frequency-domain analysis, which comprises the following steps:
Obtain N2One group of identical initial public pole in a element；
Vector fitting is carried out to each element of FDNE matrix according to the initial public pole, one by one to guarantee each member
The pole and its number of element are identical；
If pole is not converged, one group of newly public pole is obtained after vector fitting；
Vector fitting is re-started using the new public pole as the initial public pole, until the initial public pole
Difference with the new public pole obtains the FDNE matrix with public pole in predetermined threshold level；And
Solve cyclization transient current.
2. the Alloy White Iron Method for Transient Calculation according to claim 1 based on frequency-domain analysis, which is characterized in that described to obtain
Take N2One group of identical initial public pole in a element further comprises:
If sampling number is N0, then the sampling value matrix of N number of N × N-dimensional is sought；
By N in sampling value matrix each in the sampling value matrix of N number of N × N-dimensional2A element value is successively added, to obtain N0×
1 sampled value vector；
Vector fitting is carried out to the sampled value vector, obtains rational function；
The pole for obtaining the rational function obtains the initial public pole.
3. the Alloy White Iron Method for Transient Calculation according to claim 1 based on frequency-domain analysis, which is characterized in that the tool
The each element for having the FDNE matrix of public pole is a frequency-domain function:
Wherein, pole aiWith residual ciIt for real number or is complex conjugate pair, d and h are real number, and n is pole number, and s is complex frequency.
4. the Alloy White Iron Method for Transient Calculation according to claim 1 based on frequency-domain analysis, which is characterized in that described to ask
Cyclization transient current is solved, further comprises:
FDNE Equivalent Model is established, to solve cyclization transient current, wherein power transmission network voltage and current relationship can indicate
Are as follows:
Wherein, IBFor boundary bus current, IEFor other bus currents, U in addition to the bus of boundaryEFor other mothers in addition to the bus of boundary
Line voltage, UBFor boundary busbar voltage, f is frequency, and subscript B is boundary bus, and subscript E is other buses in power transmission network；
To obtain being retracted to FDNE matrix and Nuo Dun equal currents behind boundary, and the promise equal currents are as follows:
Wherein, f0For fundamental frequency.
5. the Alloy White Iron Method for Transient Calculation according to claim 4 based on frequency-domain analysis, which is characterized in that also wrap
It includes:
Three-phase current source is converted by the promise equal currents；
Difference is carried out to the three-phase current source by implicit trapezoid method.
6. a kind of Alloy White Iron Transient calculation device based on frequency-domain analysis characterized by comprising
Module is obtained, for obtaining N2One group of identical initial public pole in a element；
First fitting module, it is quasi- for carrying out vector one by one to each element of FDNE matrix according to the initial public pole
It closes, pole and its number to guarantee each element are identical；
If judgment module obtains one group of newly public pole not converged for pole after vector fitting；
Second fitting module, for re-starting vector fitting for the new public pole as the initial public pole, directly
To the difference of the initial public pole and the new public pole in predetermined threshold level, the FDNE with public pole is obtained
Matrix；And
Computing module solves cyclization transient current.
7. the Alloy White Iron Transient calculation device according to claim 6 based on frequency-domain analysis, which is characterized in that described to obtain
Modulus block further comprises:
First acquisition unit, if being N for sampling number0, then the sampling value matrix of N number of N × N-dimensional is sought；
Computing unit is used for N in sampling value matrix each in the sampling value matrix of N number of N × N-dimensional2A element value successively phase
Add, to obtain N0× 1 sampled value vector；
Fitting unit obtains rational function for carrying out vector fitting to the sampled value vector；
Second acquisition unit obtains the initial public pole for obtaining the pole of the rational function.
8. the Alloy White Iron Transient calculation device according to claim 6 based on frequency-domain analysis, which is characterized in that the tool
The each element for having the FDNE matrix of public pole is a frequency-domain function:
Wherein, pole aiWith residual ciIt for real number or is complex conjugate pair, d and h are real number, and n is pole number, and s is complex frequency.
9. the Alloy White Iron Transient calculation device according to claim 6 based on frequency-domain analysis, which is characterized in that the meter
Module is calculated, further comprises:
Unit is established, for establishing FDNE Equivalent Model, to solve cyclization transient current, wherein power transmission network voltage and current
Relationship can indicate are as follows:
Wherein, IBFor boundary bus current, IEFor other bus currents, U in addition to the bus of boundaryEFor other mothers in addition to the bus of boundary
Line voltage, UBFor boundary busbar voltage, f is frequency, and subscript B is boundary bus, and subscript E is other buses in power transmission network；
To obtain being retracted to FDNE matrix and Nuo Dun equal currents behind boundary, and the promise equal currents are as follows:
Wherein, f0For fundamental frequency.
10. the Alloy White Iron Method for Transient Calculation according to claim 9 based on frequency-domain analysis, which is characterized in that also wrap
It includes:
Three-phase current source is converted by the promise equal currents；
Difference is carried out to the three-phase current source by implicit trapezoid method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711173997.0A CN109818346A (en) | 2017-11-22 | 2017-11-22 | Alloy White Iron Method for Transient Calculation and device based on frequency-domain analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711173997.0A CN109818346A (en) | 2017-11-22 | 2017-11-22 | Alloy White Iron Method for Transient Calculation and device based on frequency-domain analysis |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109818346A true CN109818346A (en) | 2019-05-28 |
Family
ID=66601412
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711173997.0A Pending CN109818346A (en) | 2017-11-22 | 2017-11-22 | Alloy White Iron Method for Transient Calculation and device based on frequency-domain analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109818346A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113051777A (en) * | 2021-04-28 | 2021-06-29 | 北京华大九天科技股份有限公司 | Method for correcting data by using vector fitting |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103106328A (en) * | 2012-11-02 | 2013-05-15 | 南方电网科学研究院有限责任公司 | Frequency correlation network equivalence generation method based on integral vector fitting method |
| CN106528934A (en) * | 2016-10-10 | 2017-03-22 | 清华大学 | Frequency-related network equivalence passivity correction method |
-
2017
- 2017-11-22 CN CN201711173997.0A patent/CN109818346A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103106328A (en) * | 2012-11-02 | 2013-05-15 | 南方电网科学研究院有限责任公司 | Frequency correlation network equivalence generation method based on integral vector fitting method |
| CN106528934A (en) * | 2016-10-10 | 2017-03-22 | 清华大学 | Frequency-related network equivalence passivity correction method |
Non-Patent Citations (1)
| Title |
|---|
| XU ZHANG等: ""Implicitly Coupled Electromechanical and Electromagnetic Transient Analysis Using a Frequency-Dependent Network Equivalent"", 《IEEE TRANSACTIONS ON POWER DELIVERY》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113051777A (en) * | 2021-04-28 | 2021-06-29 | 北京华大九天科技股份有限公司 | Method for correcting data by using vector fitting |
| CN113051777B (en) * | 2021-04-28 | 2021-08-31 | 北京华大九天科技股份有限公司 | Method for correcting data by using vector fitting |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lin et al. | Three-phase unbalanced distribution power flow solutions with minimum data preparation | |
| CN111216592B (en) | Power battery charging remaining time estimation method and device and electric vehicle | |
| CN102403724B (en) | Symmetrical obtaining method for nodal voltage sensitivity in alternating-direct current parallel-serial power grid | |
| CN105868513B (en) | Electromagnetic transient simulation newly-added model development method based on transient calculation network equivalence | |
| CN105388392B (en) | The single-ended online Fault Locating Method of DC distribution cable based on apparent impedance identification | |
| CN103544378A (en) | Method for calculating harmonic impedance of alternating current system for direct current transmission | |
| CN111308264B (en) | Power distribution network single-phase earth fault section positioning method based on cosine similarity | |
| CN110348038A (en) | A kind of battery pack electromagnetic interference modeling and simulating and test method | |
| CN102104252A (en) | Power system dynamic equivalence method suitable for electromagnetic transient analysis | |
| Li et al. | Fault location algorithm in distribution networks considering distributed capacitive current | |
| CN109738714A (en) | The method and device of secondary nuclear phase | |
| CN107395076B (en) | A kind of low parameter identification method for encouraging limiter of broken line type | |
| CN106291046B (en) | Mixed pressure common-tower double-return line is single-phase across single-phase across voltage failure current calculation method | |
| CN109494720A (en) | A kind of voltage dip Stochastic prediction method based on Internet communication characteristic | |
| CN106093705A (en) | A kind of computational methods of one-phase earthing failure in electric distribution network wavefront | |
| CN111668843A (en) | Low-voltage distribution network three-phase load flow calculation method based on phase sequence mixing method | |
| CN103412199A (en) | Computing method of unbalancedness of multi-circuit power transmission lines on same tower | |
| DE102015113771A1 (en) | Vehicle battery Ladeystem notification | |
| Cataliotti et al. | Simulation of a power line communication system in medium and low voltage distribution networks | |
| CN105305392B (en) | Short circuit calculation symmetrical component method suitable for the type IIDG power distribution network of control containing voltage | |
| CN102684229B (en) | Random neighborhood searching method for distribution network reconfiguration | |
| CN109818346A (en) | Alloy White Iron Method for Transient Calculation and device based on frequency-domain analysis | |
| CN104167733B (en) | A kind of power grid equivalence method based on load impedance | |
| Meerimatha et al. | A novel distribution system power flow algorithm using forward backward matrix method | |
| CN109655710B (en) | Correction method, device and system for phasor measurement at two ends of same-pole double-circuit transmission line |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190528 |