CN106896267A - A kind of Distribution Network Harmonics resonance improves modal analysis method - Google Patents
A kind of Distribution Network Harmonics resonance improves modal analysis method Download PDFInfo
- Publication number
- CN106896267A CN106896267A CN201710067984.9A CN201710067984A CN106896267A CN 106896267 A CN106896267 A CN 106896267A CN 201710067984 A CN201710067984 A CN 201710067984A CN 106896267 A CN106896267 A CN 106896267A
- Authority
- CN
- China
- Prior art keywords
- node
- frequency
- resonance
- mode
- crucial
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
- G01R23/163—Spectrum analysis; Fourier analysis adapted for measuring in circuits having distributed constants
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
Modal analysis method is improved the invention discloses a kind of Distribution Network Harmonics resonance, is comprised the following steps：Determine resonant frequency using frequency spectrum analysis method, then obtain the eigenvalue matrix of the bus admittance matrix under resonant frequency, the characteristic value to bus admittance matrix carries out ascending order arrangement, records the position of each characteristic value to it by size；Take feasible solution of the preceding m minimal eigenvalue as crucial mode under resonant frequency, determine the corresponding crucial right characteristic vector of each feasible solution, maximum right eigenvalue and its position of these crucial right characteristic vector feasible solutions are can obtain, position is the feasible solution that highest excites node.The beneficial effects of the invention are as follows：The present invention has the advantages such as computational efficiency is high, highest excites node locating accurate, is more applicable for Distribution Network Harmonics harmonic analysis.
Description
Technical field
The present invention relates to Power System Analysis and control technology field, more particularly to a kind of Distribution Network Harmonics resonance improves mould
State analysis method.
Background technology
In recent years, being continuously increased with need for electricity, the extensive application of converter technique, the nonlinear-load in power distribution network
Increase year by year, the harmonic problem of power network is on the rise, the stable operation to power network causes potential hazard.On the one hand, harmonic wave is deposited
The quality of power supply of power network is being reduced, the safe and reliable operation of electrical equipment is being endangered, is increasing plant capacity loss, or even damage sets
It is standby, cause power-off fault.On the other hand, when harmonic frequency and close system resonance frequencies in power network, system may be induced
Resonance, produces resonance overvoltage, has a strong impact on the safe and reliable operation of power network.
Harmonic resonance is the important harm that harmonic wave causes, and parallel resonance causes overvoltage, series resonance to cause electricity
Stream, causes harm to power system safety and stability operation, so that suppressing the harmonic resonance of power system, this is accomplished by carrying out
Harmonic resonance is analyzed, and obtains resonant frequency and resonance range, so as to take corresponding measure to avoid the generation of harmonic resonance, is reduced humorous
Shake the breakdown loss for causing, and this stable operation to ensureing power system has important practical significance.
Harmonic resonance analysis method mainly has Spectral Analysis Method, modal analysis method etc., and Spectral Analysis Method can determine resonance
Frequency, but can not provide the more information of resonance, such as resonance location and each element are to the influence degree of resonance, traditional modal analysis
Method is wide concerned due to efficiently convenient, by the characteristic root of analysis node admittance matrix, it is possible to obtain the various letters of resonance
Breath, such as resonant frequency, highest shot point, highest observation station and participation factors, but its computational efficiency is relatively low, and cannot accurately sentence
Highest in disconnected Distribution Network Harmonics resonance excites node.Therefore, mode is improved in the urgent need to studying a kind of Distribution Network Harmonics resonance
Analysis method, the various information of Distribution Network Harmonics resonance are determined with efficiently and accurately.
The content of the invention
It is contemplated that overcoming the shortcomings of traditional modal analysis method, propose that a kind of Distribution Network Harmonics resonance improves mode point
Analysis method, the present invention determines resonant frequency first with Spectral Analysis Method analysis power distribution network node impedance extreme point position, eliminates
The step of traditional modal analysis method interior joint admittance matrix is decomposed, eigenvalue matrix is inverted, improves computational efficiency；It is simultaneously right
Traditional modal analysis method is improved, by increase judge to excite node can occur at the resonant frequency fx the condition of resonance come
Determine that highest excites node and optimal observer nodes, the positioning produced when reducing traditional modal analysis method for power distribution network is wrong
By mistake, it is capable of various information (such as resonant frequency, highest shot point, the highest sight for determining Distribution Network Harmonics resonance of efficiently and accurately
Measuring point and participation factors), the harmonic resonance analysis method is applicable not only to power distribution network, is equally applicable to power transmission network.
In order to realize foregoing invention purpose, modal analysis method is improved the invention provides a kind of Distribution Network Harmonics resonance,
Specifically include following steps：
S1, determine resonant frequency using frequency spectrum analysis method, specifically, when in system for distribution network of power i-node occur it is humorous
When vibration frequency is the parallel resonance of f, i-node impedance shows maximum, and each node maximum impedance is by formula (1) in power system
Obtain：
In formula, ziifIt is the node impedance under a certain frequency, takes famous value；
The frequency corresponding to each node impedance maximum is obtained by according to each node maximum impedance obtained, i.e., it is described
Resonant frequency；
Wherein, the frequency spectrum analysis method uses self adaptation sweep frequency technique, specifically includes following steps：
Step S101, setting original frequency f0, final frequency fmax, initial step length h0=0.2Hz, maximum step-length hmax=
10Hz (increasing step-length during frequency sweep can improve computational efficiency, but step-length is too big may be interval across some unimodal interval, single paddy,
Therefore the maximum of hmax restricted steps is set)；
Step S102, h1=h2=h0 is set gradually, f1=f0, f2=f1+h1, f3=f2+h2；
Step S103, Z is node self-impedance, is calculated
Step S104, the product by calculating the slope between twice sweep, i.e. A=k1*k2 control the size of step-length；
Step S105, work as A>When 0, i.e., frequency sweep, with incremental interval or with interval of successively decreasing, now increases step-length (h1 twice
=h2, h2=h2+0.2), advance (f1=f2, f2=f3, f3=f2+h2, k1=k2), until encounter one it is unimodal interval or
(Rule of judgment is A untill single paddy interval<0), turn 6.；During frequency sweep, if reaching interval right margin fmax, program end
Only.
Work as A<When 0, i.e., frequency sweep is interval in unimodal interval or single paddy twice, and now node impedance has in the interval of frequency sweep twice
Extreme point, [f1, f3] interval domestic demand frequency sweep again, step-length is changed into the half (h1=h1/2, h2=h12+0.2) of original steps,
Turn 3., to continue frequency sweep and find out Impedance Peak position, you can determine resonant frequency.
S2, the eigenvalue matrix that the bus admittance matrix under the resonant frequency that step S1 draws is obtained using formula (2),
Y=L Λ T (2)
Wherein, Λ is diagonal eigenmatrix, and L and T is respectively left eigenvector matrix and right eigenvectors matrix, L=T-1；
Mode voltage vector U=TV, mode current vector J=TI are defined, then U=Λ-1J, can be represented by the formula：
In formula, λ-1Referred to as mode impedance Zm, work as λi=0 or during very little, even if mode i Injection Currents JiVery little will also be produced
Very big mode i voltages Ui, claim minimum λiIt is the crucial mode of resonance, corresponding [Ti1,Ti2,…Tin]、[L1i,L2i,…Lni
]TRespectively crucial right characteristic vector and crucial left eigenvector；
S3, the characteristic value to bus admittance matrix carry out ascending order arrangement to it by size, record the position of each characteristic value；Such as
Shown in lower:
S4, the feasible solution of preceding m minimal eigenvalue in step S3 as crucial mode under resonant frequency is taken, preferably m=10,
The corresponding crucial right characteristic vector of each feasible solution is determined according to above-mentioned formula (2), you can obtain these crucial right characteristic vectors
The maximum right eigenvalue of feasible solution and its position, the position are the feasible solution that highest excites node；
S5, by mode electric current JiIt is expressed as：
Ji=Ti1I1+Ti2I2+Ti3I3+…+TinIn (4)
Work as TijWhen maximum, its corresponding IjTo mode electric current JiContribution it is maximum, illustrate that j nodes are most easily generation mode i
The position of resonance, works as TijWhen=0, no matter IjIt is much all without exciting resonance, therefore with crucial right characteristic vector [Ti1,Ti2,…
Tin] represent that each node current excites degree to mode i resonance；The excitability of crucial mode can use right key feature vector
Represent, the node with maximum right eigenvalue is that highest excites node in mode.
Further, due to correspondence certain pattern have highest can incentive node simultaneously also with highest Observable
Property, thus highest observer nodes can excite node identical with above-mentioned highest, therefore be that can obtain highest to observe by the above method
Node；It is that can determine that corresponding participation factors by formula (5)：
PFbm=LbmTmb (5)
In formula, PF is participation factors, and L is crucial left eigenvector, and T is crucial right characteristic vector, and b is bus number, and m is mould
State number.
Beneficial effects of the present invention are：Compared with prior art, the present invention has that computational efficiency is high, highest excites node to determine
Level really etc. advantage, be more applicable for Distribution Network Harmonics harmonic analysis, by improve modal analysis method by spectrum analysis with tradition
Model analysis is combined, and resonant frequency is determined using Spectral Analysis Method analysis node Impedance Peak point position, has been eliminated node and has been led
Matrix decomposition, eigenvalue matrix are received the step of invert, computational efficiency is improve；Traditional modal analysis are utilized at the resonant frequency fx
Can method, node be excited determining highest while increased and judge to excite node that the condition of resonance occur at the resonant frequency fx, is subtracted
Lack the Wrong localization of traditional modal analysis method, be capable of the various information for determining Distribution Network Harmonics resonance of efficiently and accurately, such as
Resonant frequency, highest shot point, highest observation station and participation factors.
Brief description of the drawings
Fig. 1 is the flow chart of method in the embodiment of the present invention.
Fig. 2 is the flow chart of the self adaptation frequency sweep of step S1 in the embodiment of the present invention.
Fig. 3 is the schematic diagram of embodiment of the present invention contrast test Zhong Bonanyou areas electricity distribution network model.
Fig. 4 is 62, No. 66 nodes correspondence determined using traditional modal analysis method in embodiment of the present invention contrast test
Impedance frequency characteristic curve map.
Fig. 5 is No. 62 nodes pair determined using present invention improvement modal analysis method in embodiment of the present invention contrast test
The impedance frequency characteristic curve map answered.
Specific embodiment
The present invention for mention in background technology existing harmonic resonance analysis method exist it is less efficient, and cannot be accurate
Judge the problems such as highest in Distribution Network Harmonics resonance excites node, there is provided a kind of Distribution Network Harmonics resonance improves model analysis
Method, resonant frequency is determined first with Spectral Analysis Method, obtains the eigenvalue matrix of resonant frequency lower node admittance matrix, and
Ascending order arrangement is carried out to it by size, the position of each characteristic value is recorded, preceding m minimal eigenvalue is taken as crucial under resonant frequency
The feasible solution of mode, determines the corresponding crucial right characteristic vector of each feasible solution, draws crucial right characteristic vector feasible solution most
Big right eigenvalue and its position, these positions are the feasible solution that highest excites node.By find out at the resonant frequency fx occur it is humorous
The highest shaken excites the feasible solution of node, and takes and wherein possess the feasible solution of maximum impedance value (perunit value) and excite section as highest
Point, its corresponding mode is crucial mode, and highest observer nodes are obtained using the crucial mode and traditional modal analysis method that determine
And corresponding participation factors.
Technical scheme is described in detail with reference to specific embodiment.
Embodiment 1
The embodiment of the invention provides a kind of Distribution Network Harmonics resonance and improve modal analysis method, specifically include following step
Suddenly：
S1, determine resonant frequency using frequency spectrum analysis method, specifically, when in system for distribution network of power i-node occur it is humorous
When vibration frequency is the parallel resonance of f, i-node impedance shows maximum, and each node maximum impedance is by formula (1) in power system
Obtain：
In formula, ziifIt is the node impedance under a certain frequency, takes famous value；
The frequency corresponding to each node impedance maximum is obtained by according to each node maximum impedance obtained, i.e., it is described
Resonant frequency；
Wherein, the frequency spectrum analysis method uses self adaptation sweep frequency technique, specifically includes following steps：
Step S101, setting original frequency f0, final frequency fmax, initial step length h0=0.2Hz, maximum step-length hmax=
10Hz (increasing step-length during frequency sweep can improve computational efficiency, but step-length is too big may be interval across some unimodal interval, single paddy,
Therefore the maximum of hmax restricted steps is set)；
Step S102, h1=h2=h0 is set gradually, f1=f0, f2=f1+h1, f3=f2+h2；
Step S103, Z is node self-impedance, is calculated
Step S104, the product by calculating the slope between twice sweep, i.e. A=k1*k2 control the size of step-length；
Step S105, work as A>When 0, i.e., frequency sweep, with incremental interval or with interval of successively decreasing, now increases step-length (h1 twice
=h2, h2=h2+0.2), advance (f1=f2, f2=f3, f3=f2+h2, k1=k2), until encounter one it is unimodal interval or
(Rule of judgment is A untill single paddy interval<0), turn 6.；During frequency sweep, if reaching interval right margin fmax, program end
Only.
Work as A<When 0, i.e., frequency sweep is interval in unimodal interval or single paddy twice, and now node impedance has in the interval of frequency sweep twice
Extreme point, [f1, f3] interval domestic demand frequency sweep again, step-length is changed into the half (h1=h1/2, h2=h12+0.2) of original steps,
Turn 3., to continue frequency sweep and find out Impedance Peak position, you can determine resonant frequency.
S2, the eigenvalue matrix that the bus admittance matrix under the resonant frequency that step S1 draws is obtained using formula (2),
Y=L Λ T (2)
Wherein, Λ is diagonal eigenmatrix, and L and T is respectively left eigenvector matrix and right eigenvectors matrix, L=T-1；
Mode voltage vector U=TV, mode current vector J=TI are defined, then U=Λ-1J, can be represented by the formula：
In formula, λ-1Referred to as mode impedance Zm, work as λi=0 or during very little, even if mode i Injection Currents JiVery little will also be produced
Very big mode i voltages Ui, claim minimum λiIt is the crucial mode of resonance, corresponding [Ti1,Ti2,…Tin]、[L1i,L2i,…Lni
]TRespectively crucial right characteristic vector and crucial left eigenvector；
S3, the characteristic value to bus admittance matrix carry out ascending order arrangement to it by size, record the position of each characteristic value；Such as
Shown in lower:
S4, the feasible solution of preceding m minimal eigenvalue in step S3 as crucial mode under resonant frequency is taken, preferably m=10,
The corresponding crucial right characteristic vector of each feasible solution is determined according to above-mentioned formula (2), you can obtain these crucial right characteristic vectors
The maximum right eigenvalue of feasible solution and its position, the position are the feasible solution that highest excites node；
S5, by mode electric current JiIt is expressed as：
Ji=Ti1I1+Ti2I2+Ti3I3+…+TinIn (4)
Work as TijWhen maximum, its corresponding IjTo mode electric current JiContribution it is maximum, illustrate that j nodes are most easily generation mode i
The position of resonance, works as TijWhen=0, no matter IjIt is much all without exciting resonance, therefore with crucial right characteristic vector [Ti1,Ti2,…
Tin] represent that each node current excites degree to mode i resonance；The excitability of crucial mode can use right key feature vector
Represent, the node with maximum right eigenvalue is that highest excites node in mode.
Further, due to correspondence certain pattern have highest can incentive node simultaneously also with highest Observable
Property, thus highest observer nodes can excite node identical with above-mentioned highest, therefore be that can obtain highest to observe by the above method
Node；It is that can determine that corresponding participation factors by formula (5)：
PFbm=LbmTmb (5)
In formula, PF is participation factors, and L is crucial left eigenvector, and T is crucial right characteristic vector, and b is bus number, and m is mould
State number.
Contrast experiment
This contrast experiment by taking Shengli Oil Field Bo Nanyou areas power distribution network as an example, respectively using traditional modal analysis method and the present invention
Improved modal analysis method analyzes Distribution Network Harmonics resonance respectively, contrasts the analysis result of two methods, and checking improves mode point
The superiority of analysis method.
Bo Nanyou areas electricity distribution network model as shown in figure 3, carry branched line under wherein every circuit, 733 nodes altogether,
749 branch roads；Two methods analysis result and simulation time are as shown in table 1.
The Bo Nanyou areas power distribution network modal analysis result of table 1 and simulation time
By the data of table 1, it can be seen that No. 62 node on behalf Bo Nan note -1250 nodes, No. 66 node on behalf oil plants -
1250 nodes, two resonant frequencies respectively 905.2Hz, 1545.2Hz are determined using traditional modal analysis method, resonant frequency
It is respectively 493, No. 625 nodes that highest excites node, adds the frequency filtering to be at No. 493 nodes in Bo Nanyou areas electricity distribution network model
The resonance filter of 905.2Hz, adds the resonance filter that frequency filtering is 1545.2Hz, 62, No. 66 nodes at No. 625 nodes
Impedance frequency characteristic curve is as shown in Figure 4.
62,66 node impedance frequency characteristics when red line is represented and do not add wave filter in Fig. 4, green line adds wave filter deutomerite
Point impedance frequency characteristic, as seen from the figure, adds the resonance filter that frequency filtering is 905.2Hz not at No. 493 nodes
Harmonic resonance can be suppressed, the resonance filter for adding frequency filtering to be 1545.2Hz at No. 625 nodes can effectively suppress harmonic wave
Resonance.Above analysis shows traditional modal analysis method is likely to result in highest and excites node locating when Distribution Network Harmonics resonance is analyzed
Mistake, the inapplicable Distribution Network Harmonics harmonic analysis of traditional modal analysis method.
By the data of table 1, it can be seen that No. 62 node on behalf Bo Nan note -1250 nodes, No. 66 node on behalf oil plants -
1250 nodes, improve modal analysis method and determine that the highest of resonant frequency 905.2Hz, 1545.2Hz excites node point using the present invention
It is not 63, No. 625 nodes, No. 625 nodes have verified that highest during for 1545.2Hz excites node, in Bo Nanyou areas power distribution network mould
Add the resonance filter that frequency filtering is 905.2Hz, No. 62 node impedance frequency characteristics such as Fig. 5 institutes in type at No. 63 nodes
Show.
62,66 node impedance frequency characteristics when red line is represented and do not add wave filter in Fig. 5, green line adds wave filter deutomerite
Point impedance frequency characteristic, as seen from the figure, the resonance filter for adding frequency filtering to be 905.2Hz at No. 63 nodes can have
Effect suppresses harmonic resonance.
By above-mentioned contrast, it can be seen that compared with traditional modal analysis method, improving modal analysis method in analysis power distribution network
Highest can be accurately positioned during harmonic resonance excites node, and traditional modal analysis method simulation time 38120.5s to be much larger than and change
Enter modal analysis method simulation time 6545.4s, show the computational efficiency for improving modal analysis method apparently higher than traditional modal analysis
Method computational efficiency.
In sum, although traditional modal analysis method can determine that the resonant frequency of power network, highest excite node, highest
Observer nodes and bus participation factors, but it needs to carry out bus admittance matrix decomposition, and computational efficiency is relatively low, and applies in distribution
Cannot accurately judge that highest excites node in the analysis of net harmonic resonance, improvement modal analysis method proposed by the present invention is by spectrum analysis
It is combined with traditional modal analysis method, excites the node bar of resonance can occur at the resonant frequency fx while increased judgement
Part, can efficiently, accurately determine that Distribution Network Harmonics resonant frequency, highest excite node.Improvement model analysis proposed by the present invention
Method is applied to Distribution Network Harmonics harmonic analysis.
The foregoing is only presently preferred embodiments of the present invention, be not intended to limit the invention, it is all it is of the invention spirit and
Within principle, any modification, equivalent substitution and improvements made etc. should be included within the scope of the present invention.
Claims (4)
1. a kind of Distribution Network Harmonics resonance improves modal analysis method, it is characterised in that the analysis method is comprised the following steps that：
S1, determine resonant frequency using frequency spectrum analysis method, specifically, when in system for distribution network of power i-node there is resonance frequently
When rate is the parallel resonance of f, i-node impedance shows maximum, and each node maximum impedance is asked by formula (1) in power system
Go out：
In formula, ziifIt is the node impedance under a certain frequency, takes famous value；
The frequency corresponding to each node impedance maximum, i.e., described resonance are obtained by according to each node maximum impedance obtained
Frequency；
S2, the eigenvalue matrix that the bus admittance matrix under the resonant frequency that step S1 draws is obtained using formula (2),
Y=L Λ T (2)
Wherein, Λ is diagonal eigenmatrix, and L and T is respectively left eigenvector matrix and right eigenvectors matrix, L=T-1；
Mode voltage vector U=TV, mode current vector J=TI are defined, then U=Λ-1J, can be represented by the formula：
In formula, λ-1Referred to as mode impedance Zm, work as λi=0 or during very little, even if mode i Injection Currents JiVery little will also produce very big
Mode i voltages Ui, claim minimum λiIt is the crucial mode of resonance, corresponding [Ti1,Ti2,…Tin]、[L1i,L2i,…Lni]TPoint
Wei not crucial right characteristic vector and crucial left eigenvector；
S3, the characteristic value to bus admittance matrix carry out ascending order arrangement to it by size, record the position of each characteristic value；
S4, the feasible solution of preceding m minimal eigenvalue in step S3 as crucial mode under resonant frequency is taken, according to above-mentioned formula
(2) the corresponding crucial right characteristic vector of each feasible solution is determined, you can obtain the maximum of these crucial right characteristic vector feasible solutions
Right eigenvalue and its position, the position are the feasible solution that highest excites node；
S5, by mode electric current JiIt is expressed as：
Ji=Ti1I1+Ti2I2+Ti3I3+…+TinIn (4)
Work as TijWhen maximum, its corresponding IjTo mode electric current JiContribution it is maximum, illustrate that j nodes are most easy generation mode i resonance
Position, works as TijWhen=0, no matter IjIt is much all without exciting resonance, therefore with crucial right characteristic vector [Ti1,Ti2,…Tin] represent
Each node current excites degree to mode i resonance；The excitability of crucial mode can use right key feature vector representation, tool
The node for having maximum right eigenvalue is that highest excites node, namely highest observer nodes in mode.
2. Distribution Network Harmonics resonance according to claim 1 improves modal analysis method, it is characterised in that by formula
(5) participation factors are calculated：
PFbm=LbmTmb (5)
In formula, PF is participation factors, and L is crucial left eigenvector, and T is crucial right characteristic vector, and b is bus number, and m is mode
Number.
3. Distribution Network Harmonics resonance according to claim 1 and 2 improves modal analysis method, it is characterised in that the step
Frequency spectrum analysis method in rapid S1 uses self adaptation frequency sweeping method, specially：
Step S101, setting original frequency f0, final frequency fmax, initial step length h0=0.2Hz, maximum step-length hmax=10Hz
(increasing step-length during frequency sweep can improve computational efficiency, but step-length is too big may be interval across some unimodal interval, single paddy, therefore
The maximum of hmax restricted steps is set)；
Step S102, h1=h2=h0 is set gradually, f1=f0, f2=f1+h1, f3=f2+h2；
Step S103, Z is node self-impedance, is calculated
Step S104, the product by calculating the slope between twice sweep, i.e. A=k1*k2 control the size of step-length；
Step S105, work as A>When 0, i.e., frequency sweep, with incremental interval or with interval of successively decreasing, now increases step-length (h1=twice
H2, h2=h2+0.2), advance (f1=f2, f2=f3, f3=f2+h2, k1=k2), it is unimodal interval or single until encountering one
(Rule of judgment is A untill paddy interval<0), turn 6.；During frequency sweep, if reaching interval right margin fmax, program determination.
Work as A<When 0, i.e., frequency sweep is interval in unimodal interval or single paddy twice, and now node impedance has extreme value in the interval of frequency sweep twice
Point, [f1, f3] interval domestic demand frequency sweep again, step-length is changed into the half (h1=h1/2, h2=h12+0.2) of original steps, turns 3.,
Continue frequency sweep and find out Impedance Peak position, you can determine resonant frequency.
4. the Distribution Network Harmonics resonance according to claim any one of 1-3 improves modal analysis method, it is characterised in that institute
In stating step S4, the feasible solution of preceding m=10 minimal eigenvalue in step S3 as crucial mode under resonant frequency is taken.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710067984.9A CN106896267B (en) | 2017-02-07 | 2017-02-07 | A kind of Distribution Network Harmonics resonance improvement modal analysis method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710067984.9A CN106896267B (en) | 2017-02-07 | 2017-02-07 | A kind of Distribution Network Harmonics resonance improvement modal analysis method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106896267A true CN106896267A (en) | 2017-06-27 |
| CN106896267B CN106896267B (en) | 2019-05-10 |
Family
ID=59198113
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710067984.9A Active CN106896267B (en) | 2017-02-07 | 2017-02-07 | A kind of Distribution Network Harmonics resonance improvement modal analysis method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106896267B (en) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107453358A (en) * | 2017-09-15 | 2017-12-08 | 浙江大学 | A kind of electric network natural resonance structure analysis method based on complex frequency domain lower node admittance matrix |
| CN108268999A (en) * | 2017-12-12 | 2018-07-10 | 国网河南省电力公司电力科学研究院 | For more pivot modal analysis methods of electric power networks node Harmonic characteristic analysis |
| CN109613324A (en) * | 2018-11-07 | 2019-04-12 | 全球能源互联网研究院有限公司 | A kind of detection method and device of Harmonics amplification |
| CN110718917A (en) * | 2019-11-18 | 2020-01-21 | 国网青海省电力公司 | Implementation method for constructing generation of power electronic harmonic resonance source |
| CN110889081A (en) * | 2019-11-25 | 2020-03-17 | 国网河南省电力公司驻马店供电公司 | Method for quickly acquiring harmonic resonance of new energy grid-connected system |
| CN112034251A (en) * | 2020-08-28 | 2020-12-04 | 武汉大学 | Method for evaluating parallel resonance excitation capability of Nonton type inter-harmonic source access point |
| CN112072662A (en) * | 2020-08-28 | 2020-12-11 | 武汉大学 | Method for avoiding inter-harmonic parallel resonance of one frequency band |
| CN112701689A (en) * | 2020-12-16 | 2021-04-23 | 深圳供电局有限公司 | Wide-area resonance evaluation and early warning method based on limited distribution points |
| CN112858784A (en) * | 2021-04-03 | 2021-05-28 | 国网四川省电力公司电力科学研究院 | Traction power supply system-regional power grid parallel harmonic resonance frequency identification method |
| CN113612231A (en) * | 2021-07-28 | 2021-11-05 | 合肥工业大学 | Impedance adapter address selecting method based on modal analysis method |
| CN115173417A (en) * | 2022-06-20 | 2022-10-11 | 中国石油化工股份有限公司 | Comprehensive identification and treatment method for resonance overvoltage of shale gas field power distribution network |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103474966A (en) * | 2013-09-12 | 2013-12-25 | 国家电网公司 | Capacitor harmonic resonance overvoltage detection method and protection device |
| CN103580025A (en) * | 2013-08-31 | 2014-02-12 | 西南交通大学 | Adaptive filtering method for resonant frequency iteration transfer based on modal frequency sensitivity |
| CN105529727A (en) * | 2016-02-22 | 2016-04-27 | 哈尔滨工业大学 | Node admittance matrix eigenvalue analysis method applied to grid-connected-inverter-included parallel resonance situation |
| CN106250660A (en) * | 2016-08-31 | 2016-12-21 | 上海交通大学 | A kind of marine wind electric field harmonic analysis method |
| CN106295231A (en) * | 2016-08-31 | 2017-01-04 | 上海交通大学 | A kind of marine wind electric field resonance suppressing method based on modal calculation |
-
2017
- 2017-02-07 CN CN201710067984.9A patent/CN106896267B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103580025A (en) * | 2013-08-31 | 2014-02-12 | 西南交通大学 | Adaptive filtering method for resonant frequency iteration transfer based on modal frequency sensitivity |
| CN103474966A (en) * | 2013-09-12 | 2013-12-25 | 国家电网公司 | Capacitor harmonic resonance overvoltage detection method and protection device |
| CN105529727A (en) * | 2016-02-22 | 2016-04-27 | 哈尔滨工业大学 | Node admittance matrix eigenvalue analysis method applied to grid-connected-inverter-included parallel resonance situation |
| CN106250660A (en) * | 2016-08-31 | 2016-12-21 | 上海交通大学 | A kind of marine wind electric field harmonic analysis method |
| CN106295231A (en) * | 2016-08-31 | 2017-01-04 | 上海交通大学 | A kind of marine wind electric field resonance suppressing method based on modal calculation |
Non-Patent Citations (2)
| Title |
|---|
| 徐文远 等: "基于模态分析的谐波谐振评估方法", 《中国电机工程学报》 * |
| 蒋平 等: "谐波分析及谐振评估软件的开发与应用", 《电力自动化设备》 * |
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107453358B (en) * | 2017-09-15 | 2019-10-15 | 浙江大学 | A kind of electric network natural resonance structure analysis method based on complex frequency domain lower node admittance matrix |
| CN107453358A (en) * | 2017-09-15 | 2017-12-08 | 浙江大学 | A kind of electric network natural resonance structure analysis method based on complex frequency domain lower node admittance matrix |
| CN108268999B (en) * | 2017-12-12 | 2021-05-25 | 国网河南省电力公司电力科学研究院 | Multi-principal-element modal analysis method for harmonic characteristic analysis of power network nodes |
| CN108268999A (en) * | 2017-12-12 | 2018-07-10 | 国网河南省电力公司电力科学研究院 | For more pivot modal analysis methods of electric power networks node Harmonic characteristic analysis |
| CN109613324A (en) * | 2018-11-07 | 2019-04-12 | 全球能源互联网研究院有限公司 | A kind of detection method and device of Harmonics amplification |
| CN109613324B (en) * | 2018-11-07 | 2021-10-08 | 全球能源互联网研究院有限公司 | Harmonic amplification detection method and device |
| CN110718917A (en) * | 2019-11-18 | 2020-01-21 | 国网青海省电力公司 | Implementation method for constructing generation of power electronic harmonic resonance source |
| CN110718917B (en) * | 2019-11-18 | 2022-03-25 | 国网青海省电力公司 | Implementation method for constructing generation of power electronic harmonic resonance source |
| CN110889081A (en) * | 2019-11-25 | 2020-03-17 | 国网河南省电力公司驻马店供电公司 | Method for quickly acquiring harmonic resonance of new energy grid-connected system |
| CN112034251A (en) * | 2020-08-28 | 2020-12-04 | 武汉大学 | Method for evaluating parallel resonance excitation capability of Nonton type inter-harmonic source access point |
| CN112072662A (en) * | 2020-08-28 | 2020-12-11 | 武汉大学 | Method for avoiding inter-harmonic parallel resonance of one frequency band |
| CN112034251B (en) * | 2020-08-28 | 2021-11-23 | 武汉大学 | Method for evaluating parallel resonance excitation capability of Nonton type inter-harmonic source access point |
| CN112072662B (en) * | 2020-08-28 | 2022-02-22 | 武汉大学 | Method for avoiding inter-harmonic parallel resonance of one frequency band |
| CN112701689A (en) * | 2020-12-16 | 2021-04-23 | 深圳供电局有限公司 | Wide-area resonance evaluation and early warning method based on limited distribution points |
| CN112858784A (en) * | 2021-04-03 | 2021-05-28 | 国网四川省电力公司电力科学研究院 | Traction power supply system-regional power grid parallel harmonic resonance frequency identification method |
| CN112858784B (en) * | 2021-04-03 | 2023-05-26 | 国网四川省电力公司电力科学研究院 | Traction power supply system-regional power grid parallel harmonic resonance frequency identification method |
| CN113612231A (en) * | 2021-07-28 | 2021-11-05 | 合肥工业大学 | Impedance adapter address selecting method based on modal analysis method |
| CN113612231B (en) * | 2021-07-28 | 2022-12-02 | 合肥工业大学 | Impedance adapter addressing method based on modal analysis method |
| CN115173417A (en) * | 2022-06-20 | 2022-10-11 | 中国石油化工股份有限公司 | Comprehensive identification and treatment method for resonance overvoltage of shale gas field power distribution network |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106896267B (en) | 2019-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106896267A (en) | A kind of Distribution Network Harmonics resonance improves modal analysis method | |
| CN106250626B (en) | A kind of extra-high voltage transformer zero load D.C. magnetic biasing quick calculation method based on L-I curve | |
| Lin | Power harmonics and interharmonics measurement using recursive group-harmonic power minimizing algorithm | |
| CN105160058B (en) | A kind of network boundary equivalence method calculated for electric network electromagnet transient | |
| CN110221237A (en) | A kind of capacitance type potential transformer harmonic voltage measurement method and apparatus | |
| Gustavsen | A hybrid measurement approach for wideband characterization and modeling of power transformers | |
| CN106226661A (en) | Use detection structure and the method for succusion detection reactor turn-to-turn insulation fault | |
| CN103544377A (en) | Calculating and analyzing method for resonance characteristics of direct current circuit of high-voltage direct current power transmission | |
| CN103018534A (en) | Method and system for determining harmonic voltage | |
| CN206594257U (en) | Air reactor turn-to-turn insulation intelligent detection device | |
| CN110554274A (en) | adaptive weight grounding line selection method based on wavelet singular information | |
| CN108490379B (en) | Self-excited oscillation wave-based transformer winding wave process calibration method | |
| CN102185314A (en) | Power quality comprehensive evaluation method | |
| CN106340876A (en) | Multi-direct-current feed-in city power transmission network harmonic wave amplification characteristic analysis method | |
| CN103675560A (en) | Online identifying method for detuning of direct current side filter of high-voltage direct current transmission system | |
| CN105277913A (en) | Capacitive-type voltage transformer testing method | |
| CN104297706A (en) | Magnetic field generation device based on resonance and design method thereof | |
| CN105024352B (en) | The method and system of overvoltage between suppression converter bridge | |
| CN107390025A (en) | Power system method for distinguishing multiple harmonic sources based on blind source separating | |
| Mwaniki et al. | Characterizing power transformer frequency responses using bipolar pseudo-random current impulses | |
| CN108735488A (en) | A kind of design of transformer method based on self-oscillation wave | |
| CN206960539U (en) | The harmonic elimination apparatus of capacitance type potential transformer | |
| CN108037424A (en) | Power capacitor Condition assessment of insulation method and system based on multiple data origin | |
| CN109856506B (en) | Single-phase earth fault area positioning method based on adjacent point difference method | |
| CN107345983A (en) | Multi-harmonic Sources system harmonicses transmitting appraisal procedure based on subharmonic source correlation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |