CN106896267B - A kind of Distribution Network Harmonics resonance improvement modal analysis method - Google Patents
A kind of Distribution Network Harmonics resonance improvement modal analysis method Download PDFInfo
- Publication number
- CN106896267B CN106896267B CN201710067984.9A CN201710067984A CN106896267B CN 106896267 B CN106896267 B CN 106896267B CN 201710067984 A CN201710067984 A CN 201710067984A CN 106896267 B CN106896267 B CN 106896267B
- Authority
- CN
- China
- Prior art keywords
- node
- frequency
- resonance
- mode
- analysis method
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
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
Abstract
The invention discloses a kind of Distribution Network Harmonics resonance to improve modal analysis method, the following steps are included: determining resonance frequency using frequency spectrum analysis method, the eigenvalue matrix of the node admittance matrix under resonance frequency is found out again, ascending order arrangement is carried out to it by size to the characteristic value of node admittance matrix, records the position of each characteristic value;Feasible solution of the m minimal eigenvalue as mode crucial under resonance frequency before taking, determine the corresponding right feature vector of key of each feasible solution, maximum right eigenvalue and its position of these crucial right feature vector feasible solutions, position, that is, highest excitation node feasible solution can be obtained.The advantages such as computational efficiency is high, highest excitation node locating is accurate the beneficial effects of the present invention are: the present invention has, are more applicable for Distribution Network Harmonics harmonic analysis.
Description
Technical field
The present invention relates to Power System Analysis and control technology field, in particular to a kind of Distribution Network Harmonics resonance improves mould
State analysis method.
Background technique
In recent years, being continuously increased with power demand, the extensive application of converter technique, the nonlinear-load in power distribution network
Increase year by year, the harmonic problem of power grid is on the rise, and causes potential hazard to the stable operation of power grid.On the one hand, harmonic wave is deposited
In the power quality for reducing power grid, the safe and reliable operation of electrical equipment is endangered, increases plant capacity loss, or even damage is set
It is standby, cause power-off fault.On the other hand, when the harmonic frequency and close system resonance frequencies in power grid, it may induce system
Resonance generates resonance overvoltage, seriously affects the safe and reliable operation of power grid.
Harmonic resonance is an important harm caused by harmonic wave, and parallel resonance causes overvoltage, and series resonance caused electricity
Stream, causes damages to safe and stable operation of power system, thus needs to inhibit the harmonic resonance of electric system, this just needs to carry out
Harmonic resonance analysis, obtains resonance frequency and resonance range, so that corresponding measure be taken to avoid the generation of harmonic resonance, reduces humorous
Breakdown loss caused by shaking, this has important practical significance to the stable operation for ensureing electric system.
Harmonic resonance analysis method mainly has Spectral Analysis Method, modal analysis method etc., and Spectral Analysis Method can determine resonance
Frequency, but the more information of resonance cannot be provided, the influence degree such as resonance location and each element to resonance, traditional modal analysis
Method is widely noticed due to efficiently convenient, by the characteristic root of analysis node admittance matrix, so that it may obtain the various letters of resonance
Breath, such as resonance frequency, highest excitation point, highest observation point and participates in the factor, but its computational efficiency is lower, and can not accurately sentence
Highest in disconnected Distribution Network Harmonics resonance excites node.Therefore, there is an urgent need to study a kind of Distribution Network Harmonics resonance improvement mode
Analysis method determines the various information of Distribution Network Harmonics resonance with efficiently and accurately.
Summary of the invention
The present invention is directed to overcome the shortcomings of traditional modal analysis method, propose that a kind of Distribution Network Harmonics resonance improves mode point
Analysis method, the present invention analyze power distribution network node impedance extreme point position first with Spectral Analysis Method, determine resonance frequency, eliminate
The step of traditional modal analysis method interior joint admittance matrix decomposes, eigenvalue matrix is inverted, improves computational efficiency;It is right simultaneously
Traditional modal analysis method is improved, by increase judge excitation node can occur at the resonant frequency fx the condition of resonance come
It determines highest excitation node and best observer nodes, it is wrong to reduce positioning of the traditional modal analysis method for generating when power distribution network
Accidentally, it is capable of various information (such as resonance frequency, highest excitation point, the highest sight of determining Distribution Network Harmonics resonance of efficiently and accurately
Measuring point and the participation factor), which is applicable not only to power distribution network, is equally applicable to power transmission network.
In order to achieve the above-mentioned object of the invention, the present invention provides a kind of Distribution Network Harmonics resonance to improve modal analysis method,
Specifically includes the following steps:
S1, resonance frequency is determined using frequency spectrum analysis method, specifically, it is humorous to work as i-node generation in system for distribution network of power
When vibration frequency is the parallel resonance of f, i-node impedance shows maximum value, and each node maximum impedance is by formula (1) in electric system
It finds out:
In formula, ziifFor the node impedance under a certain frequency, famous value is taken;
Frequency corresponding to each node impedance maximum value can be found out according to each node maximum impedance for finding out, i.e., described
Resonance frequency;
Wherein, the frequency spectrum analysis method uses adaptive sweep frequency technique, specifically includes the following steps:
Step S101, original frequency f0, final frequency fmax, initial step length h0=0.2Hz, maximum step-length hmax=are set
10Hz (increasing step-length when frequency sweep can be improved computational efficiency, but step-length too it is big may across some unimodal sections, the section Dan Gu,
Therefore the maximum value of hmax restricted step is set);
Step S102, h1=h2=h0, f1=f0, f2=f1+h1, f3=f2+h2 are set gradually;
Step S103, Z is node self-impedance, is calculated
Step S104, by calculating the product of the slope between twice sweep, i.e. A=k1*k2 controls the size of step-length;
Step S105, as A > 0, i.e., frequency sweep increases step-length (h1 with increasing section or in section of successively decreasing at this time twice
=h2, h2=h2+0.2), advance (f1=f2, f2=f3, f3=f2+h2, k1=k2), until encounter a unimodal section or
Until the section Dan Gu (Rule of judgment is A < 0), turn 6.;During frequency sweep, if reaching section right margin fmax, program is whole
Only.
As A < 0, i.e., frequency sweep is in unimodal section or the section Dan Gu twice, and node impedance has in frequency sweep section twice at this time
Extreme point, frequency sweep, step-length become the half (h1=h1/2, h2=h12+0.2) of original steps to the section [f1, f3] domestic demand again,
Turn 3., to continue frequency sweep and find out Impedance Peak position, that is, can determine resonance frequency.
S2, the eigenvalue matrix that the node admittance matrix under the resonance frequency that step S1 is obtained is found out using formula (2),
Y=L Λ T (2)
Wherein, Λ is diagonal eigenmatrix, and L and T are 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 when very little, even if mode i Injection Current JiVery little will also generate
Very big mode i voltage Ui, claim the smallest λiFor the crucial mode of resonance, corresponding [Ti1,Ti2,…Tin]、[L1i,L2i,…Lni]T
Respectively crucial right feature vector and crucial left eigenvector;
S3, ascending order arrangement is carried out to it by size to the characteristic value of node admittance matrix, records the position of each characteristic value;Such as
Shown in lower:
S4, take preceding m minimal eigenvalue in step S3 as the feasible solution of mode crucial under resonance frequency, preferably m=10,
The corresponding right feature vector of key of each feasible solution is determined according to above-mentioned formula (2), these crucial right feature vectors can be obtained
The maximum right eigenvalue of feasible solution and its position, the position, that is, highest excitation node feasible solution;
S5, by mode electric current JiIt indicates are as follows:
Ji=Ti1I1+Ti2I2+Ti3I3+…+TinIn (4)
Work as TijWhen maximum, corresponding IjTo mode electric current JiContribution it is maximum, illustrate that j node is that mode i most easily occurs
The position of resonance, works as TijWhen=0, no matter IjIt is much all to excite resonance, therefore with crucial right feature vector [Ti1,Ti2,…
Tin] indicate each node current to the excitation degree of mode i resonance;The excitability of crucial mode can use right key feature vector
It indicates, having the node of maximum right eigenvalue is that highest excites node in mode.
Further, due to certain corresponding mode have highest can incentive node simultaneously also with highest Observable
Property, thus highest observer nodes can excite node identical with above-mentioned highest, therefore highest observation can be obtained by the above method
Node;It can determine the corresponding participation factor by formula (5):
PFbm=LbmTmb (5)
In formula, PF is to participate in the factor, and L is crucial left eigenvector, and T is crucial right feature vector, and b is bus number, and m is mould
State number.
The invention has the benefit that compared with prior art, the present invention has computational efficiency height, highest excitation node fixed
Level really etc. advantages, be more applicable for Distribution Network Harmonics harmonic analysis, by improve modal analysis method by spectrum analysis and tradition
Model analysis combines, and determines resonance frequency using Spectral Analysis Method analysis node Impedance Peak point position, has eliminated node and led
It receives matrix decomposition, eigenvalue matrix the step of inverting, improves computational efficiency;Traditional modal analysis is utilized at the resonant frequency fx
Method, while increasing judgement and node is excited the condition of resonance can to occur at the resonant frequency fx to determine that highest excites node, subtract
The positioning mistake for having lacked traditional modal analysis method is capable of the various information for determining Distribution Network Harmonics resonance of efficiently and accurately, such as
Resonance frequency, highest excitation point, highest observation point and the participation factor.
Detailed description of the invention
Fig. 1 is the flow chart of method in the embodiment of the present invention.
Fig. 2 is the flow chart of the adaptive frequency sweep of step S1 in the embodiment of the present invention.
Fig. 3 is the schematic diagram of the area comparative test of embodiment of the present invention Zhong Bonanyou electricity distribution network model.
Fig. 4 is to be corresponded in comparative test of the embodiment of the present invention using 62, No. 66 nodes that traditional modal analysis method determines
Impedance frequency characteristic curve graph.
Fig. 5 is to improve No. 62 nodes pair that modal analysis method determines using the present invention in comparative test of the embodiment of the present invention
The impedance frequency characteristic curve graph answered.
Specific embodiment
The present invention is lower for mentioning existing harmonic resonance analysis method in background technique there are efficiency, and can not be accurate
The problems such as judging the highest excitation node in Distribution Network Harmonics resonance, provides a kind of Distribution Network Harmonics resonance improvement model analysis
Method determines resonance frequency first with Spectral Analysis Method, finds out the eigenvalue matrix of resonance frequency lower node admittance matrix, and
Ascending order arrangement is carried out to it by size, records the position of each characteristic value, m minimal eigenvalue is as crucial under resonance frequency before taking
The feasible solution of mode determines the corresponding right feature vector of key of each feasible solution, obtains crucial right feature vector feasible solution most
Big right eigenvalue and its position, these positions, that is, highest excitation node feasible solution.By find out occur at the resonant frequency fx it is humorous
The feasible solution of the highest excitation node of vibration, and take the feasible solution for wherein possessing maximum impedance value (per unit value) to excite as highest and save
Point, corresponding mode are crucial mode, obtain highest observer nodes using determining crucial mode and traditional modal analysis method
And the factor is participated in accordingly.
Technical solution of the present invention is described in detail combined with specific embodiments below.
Embodiment 1
The embodiment of the invention provides a kind of Distribution Network Harmonics resonance to improve modal analysis method, specifically includes following step
It is rapid:
S1, resonance frequency is determined using frequency spectrum analysis method, specifically, it is humorous to work as i-node generation in system for distribution network of power
When vibration frequency is the parallel resonance of f, i-node impedance shows maximum value, and each node maximum impedance is by formula (1) in electric system
It finds out:
In formula, ziifFor the node impedance under a certain frequency, famous value is taken;
Frequency corresponding to each node impedance maximum value can be found out according to each node maximum impedance for finding out, i.e., described
Resonance frequency;
Wherein, the frequency spectrum analysis method uses adaptive sweep frequency technique, specifically includes the following steps:
Step S101, original frequency f0, final frequency fmax, initial step length h0=0.2Hz, maximum step-length hmax=are set
10Hz (increasing step-length when frequency sweep can be improved computational efficiency, but step-length too it is big may across some unimodal sections, the section Dan Gu,
Therefore the maximum value of hmax restricted step is set);
Step S102, h1=h2=h0, f1=f0, f2=f1+h1, f3=f2+h2 are set gradually;
Step S103, Z is node self-impedance, is calculated
Step S104, by calculating the product of the slope between twice sweep, i.e. A=k1*k2 controls the size of step-length;
Step S105, as A > 0, i.e., frequency sweep increases step-length (h1 with increasing section or in section of successively decreasing at this time twice
=h2, h2=h2+0.2), advance (f1=f2, f2=f3, f3=f2+h2, k1=k2), until encounter a unimodal section or
Until the section Dan Gu (Rule of judgment is A < 0), turn 6.;During frequency sweep, if reaching section right margin fmax, program is whole
Only.
As A < 0, i.e., frequency sweep is in unimodal section or the section Dan Gu twice, and node impedance has in frequency sweep section twice at this time
Extreme point, frequency sweep, step-length become the half (h1=h1/2, h2=h12+0.2) of original steps to the section [f1, f3] domestic demand again,
Turn 3., to continue frequency sweep and find out Impedance Peak position, that is, can determine resonance frequency.
S2, the eigenvalue matrix that the node admittance matrix under the resonance frequency that step S1 is obtained is found out using formula (2),
Y=L Λ T (2)
Wherein, Λ is diagonal eigenmatrix, and L and T are 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 when very little, even if mode i Injection Current JiVery little will also generate
Very big mode i voltage Ui, claim the smallest λiFor the crucial mode of resonance, corresponding [Ti1,Ti2,…Tin]、[L1i,L2i,…Lni]T
Respectively crucial right feature vector and crucial left eigenvector;
S3, ascending order arrangement is carried out to it by size to the characteristic value of node admittance matrix, records the position of each characteristic value;Such as
Shown in lower:
S4, take preceding m minimal eigenvalue in step S3 as the feasible solution of mode crucial under resonance frequency, preferably m=10,
The corresponding right feature vector of key of each feasible solution is determined according to above-mentioned formula (2), these crucial right feature vectors can be obtained
The maximum right eigenvalue of feasible solution and its position, the position, that is, highest excitation node feasible solution;
S5, by mode electric current JiIt indicates are as follows:
Ji=Ti1I1+Ti2I2+Ti3I3+…+TinIn (4)
Work as TijWhen maximum, corresponding IjTo mode electric current JiContribution it is maximum, illustrate that j node is that mode i most easily occurs
The position of resonance, works as TijWhen=0, no matter IjIt is much all to excite resonance, therefore with crucial right feature vector [Ti1,Ti2,…
Tin] indicate each node current to the excitation degree of mode i resonance;The excitability of crucial mode can use right key feature vector
It indicates, having the node of maximum right eigenvalue is that highest excites node in mode.
Further, due to certain corresponding mode have highest can incentive node simultaneously also with highest Observable
Property, thus highest observer nodes can excite node identical with above-mentioned highest, therefore highest observation can be obtained by the above method
Node;It can determine the corresponding participation factor by formula (5):
PFbm=LbmTmb (5)
In formula, PF is to participate in the factor, and L is crucial left eigenvector, and T is crucial right feature vector, and b is bus number, and m is mould
State number.
Comparative experiments
This comparative experiments applies traditional modal analysis method and the present invention by taking the area Shengli Oil Field Bo Nanyou power distribution network as an example respectively
Improved modal analysis method analyzes Distribution Network Harmonics resonance respectively, compares the analysis of two methods as a result, verifying improves mode point
The superiority of analysis method.
The area Bo Nanyou electricity distribution network model as shown in figure 3, wherein under every route have branched line, amount to 733 nodes,
749 branches;Two methods analyze result and simulation time is as shown in table 1.
1 area Bo Nanyou power distribution network modal analysis result of table and simulation time
By the data of table 1, it can be seen that No. 62 node on behalf Bo Nan infuse -1250 nodes, No. 66 node on behalf oil plants -
1250 nodes determine that two resonance frequencies are respectively 905.2Hz, 1545.2Hz using traditional modal analysis method, resonance frequency
It is 493, No. 625 nodes respectively that highest, which excites node, adds the frequency filtering to be at No. 493 nodes in the area Bo Nanyou electricity distribution network model
The resonance filter of 905.2Hz, it is the resonance filter of 1545.2Hz, 62, No. 66 nodes that frequency filtering is added at No. 625 nodes
Impedance frequency characteristic curve is as shown in Figure 4.
Filter deutomerite is added in 62,66 node impedance frequency characteristics when red line indicates not add filter in Fig. 4, green line
Point impedance frequency characteristic adds the resonance filter that frequency filtering is 905.2Hz not as seen from the figure at No. 493 nodes
It can inhibit harmonic resonance, it is that the resonance filter of 1545.2Hz can effectively inhibit harmonic wave that frequency filtering is added at No. 625 nodes
Resonance.Above analysis shows traditional modal analysis method may cause highest excitation node locating when analyzing Distribution Network Harmonics resonance
Mistake, traditional modal analysis method are not suitable for Distribution Network Harmonics harmonic analysis.
By the data of table 1, it can be seen that No. 62 node on behalf Bo Nan infuse -1250 nodes, No. 66 node on behalf oil plants -
1250 nodes improve the highest excitation node point that modal analysis method determines resonance frequency 905.2Hz, 1545.2Hz using the present invention
It is not 63, No. 625 nodes, highest when No. 625 nodes have verified that as 1545.2Hz excites node, in the area Bo Nanyou power distribution network mould
Adding frequency filtering in type at No. 63 nodes is the resonance filter of 905.2Hz, No. 62 node impedance frequency characteristics such as Fig. 5 institute
Show.
Filter deutomerite is added in 62,66 node impedance frequency characteristics when red line indicates not add filter in Fig. 5, green line
Point impedance frequency characteristic, as seen from the figure, it is that the resonance filter of 905.2Hz can have that frequency filtering is added at No. 63 nodes
Effect inhibits harmonic resonance.
Through above-mentioned comparison, it can be seen that compared with traditional modal analysis method, modal analysis method is improved in analysis power distribution network
Highest excitation node can be accurately positioned when harmonic resonance, and traditional modal analysis method simulation time 38120.5s is much larger than and changes
Into modal analysis method simulation time 6545.4s, show that the computational efficiency for improving modal analysis method is apparently higher than traditional modal analysis
Method computational efficiency.
In conclusion although traditional modal analysis method can determine the resonance frequency of power grid, highest excitation node, highest
Observer nodes and bus participate in the factor, but it needs to carry out node admittance matrix decomposition, and computational efficiency is lower, and applies in distribution
Highest excitation node is unable to judge accurately 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, while increasing judgement excitation node the item of resonance can occur at the resonant frequency fx
Part can efficiently, accurately determine that Distribution Network Harmonics resonance frequency, highest excite node.Improvement model analysis proposed by the present invention
Method is suitable for Distribution Network Harmonics harmonic analysis.
The foregoing is merely presently preferred embodiments of the present invention, is not intended to limit the invention, it is all in spirit of the invention and
Within principle, any modification, equivalent replacement, improvement and so on be should all be included in the protection scope of the present invention.
Claims (3)
1. a kind of Distribution Network Harmonics resonance improves modal analysis method, which is characterized in that specific step is as follows for the analysis method:
S1, resonance frequency is determined using frequency spectrum analysis method, specifically, when i-node generation resonance frequency in system for distribution network of power
When rate is the parallel resonance of f, i-node impedance shows maximum value, and each node maximum impedance is asked by formula (1) in electric system
Out:
In formula, ziifFor the node impedance under a certain frequency, famous value is taken;
Frequency corresponding to each node impedance maximum value, the i.e. resonance can be found out according to each node maximum impedance found out
Frequency;
S2, the eigenvalue matrix that the node admittance matrix under the resonance frequency that step S1 is obtained is found out using formula (2),
Y=L Λ T (2)
Wherein, Λ is diagonal eigenmatrix, and L and T are 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 when very little, even if mode i Injection Current JiVery little will also generate very big
Mode i voltage Ui, claim the smallest λiFor the crucial mode of resonance, corresponding [Ti1,Ti2,…Tin]、[L1i,L2i,…Lni]TPoint
It Wei not crucial right feature vector and crucial left eigenvector;
S3, ascending order arrangement is carried out to it by size to the characteristic value of node admittance matrix, records the position of each characteristic value;
S4, feasible solution of the preceding m minimal eigenvalue as mode crucial under resonance frequency in step S3 is taken, according to above-mentioned formula
(2) it determines the corresponding right feature vector of key of each feasible solution, the maximum of these crucial right feature vector feasible solutions can be obtained
Right eigenvalue and its position, the position, that is, highest excitation node feasible solution;
S5, by mode electric current JiIt indicates are as follows:
Ji=Ti1I1+Ti2I2+Ti3I3+…+TinIn (4)
Work as TijWhen maximum, corresponding IjTo mode electric current JiContribution it is maximum, illustrate that j node is that mode i resonance most easily occurs
T is worked as in positionijWhen=0, no matter IjIt is much all to excite resonance, therefore with crucial right feature vector [Ti1,Ti2,…Tin] indicate
Excitation degree of each node current to mode i resonance;The excitability of crucial mode can indicate have with right key feature vector
The node of maximum right eigenvalue is highest excitation node namely highest observer nodes in mode;
Wherein, the frequency spectrum analysis method in the step S1 uses adaptive frequency sweeping method, specifically:
Step S101, original frequency f0, final frequency fmax, initial step length h0=0.2Hz, maximum step-length hmax=are set
10Hz, increasing step-length when frequency sweep can be improved computational efficiency, but step-length too it is big may across some unimodal sections, the section Dan Gu,
Therefore the maximum value of hmax restricted step is set;
Step S102, h1=h2=h0, f1=f0, f2=f1+h1, f3=f2+h2 are set gradually;
Step S103, Z is node self-impedance, is calculated
Step S104, by calculating the product of the slope between twice sweep, i.e. A=k1*k2 controls the size of step-length;
Step S105, as A > 0, i.e., frequency sweep increases step-length (h1=with increasing section or in section of successively decreasing at this time twice
H2, h2=h2+0.2), advance (f1=f2, f2=f3, f3=f2+h2, k1=k2), until encountering a unimodal section or list
Until paddy section (Rule of judgment is A < 0), S106 is gone to step;During frequency sweep, if reaching section right margin fmax, journey
Sequence terminates;
Step S106, as A < 0, i.e., frequency sweep is in unimodal section or the section Dan Gu twice, and node impedance is in frequency sweep area twice at this time
Interior to have extreme point, frequency sweep, step-length become half (h1=h1/2, the h2=h12+ of original steps to the section [f1, f3] domestic demand again
0.2) S103, is gone to step, continues frequency sweep and finds out Impedance Peak position, that is, can determine resonance frequency.
2. Distribution Network Harmonics resonance according to claim 1 improves modal analysis method, which is characterized in that pass through formula
(5) the participation factor is calculated:
In formula, PF is to participate in the factor, and L is crucial left eigenvector, and T is crucial right feature vector, and b is bus number, and m is mode
Number.
3. Distribution Network Harmonics resonance according to claim 1 or 2 improves modal analysis method, which is characterized in that the step
In rapid S4, feasible solution of the preceding m=10 minimal eigenvalue as mode crucial under resonance frequency in step S3 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 CN106896267A (en) | 2017-06-27 |
CN106896267B true 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) |
Families Citing this family (11)
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 |
CN108268999B (en) * | 2017-12-12 | 2021-05-25 | 国网河南省电力公司电力科学研究院 | Multi-principal-element modal analysis method for harmonic characteristic analysis of power network nodes |
CN109613324B (en) * | 2018-11-07 | 2021-10-08 | 全球能源互联网研究院有限公司 | Harmonic amplification detection method and device |
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 |
CN112072662B (en) * | 2020-08-28 | 2022-02-22 | 武汉大学 | 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 |
CN112701689B (en) * | 2020-12-16 | 2023-01-20 | 深圳供电局有限公司 | Wide-area resonance evaluation and early warning method based on limited distribution points |
CN112858784B (en) * | 2021-04-03 | 2023-05-26 | 国网四川省电力公司电力科学研究院 | Traction power supply system-regional power grid parallel harmonic resonance frequency identification 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 |
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 |
---|
基于模态分析的谐波谐振评估方法;徐文远 等;《中国电机工程学报》;20051130;第25卷(第22期);第1-2节 |
谐波分析及谐振评估软件的开发与应用;蒋平 等;《电力自动化设备》;20090430;第29卷(第4期);第126-130页 |
Also Published As
Publication number | Publication date |
---|---|
CN106896267A (en) | 2017-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106896267B (en) | A kind of Distribution Network Harmonics resonance improvement modal analysis method | |
CN108667048B (en) | Frequency domain stability judging method and device for oscillation stability of new energy grid-connected system | |
Hasan et al. | Harmonic compensation analysis in offshore wind power plants using hybrid filters | |
WO2019154135A1 (en) | Microgrid reliability evaluation method based on network survivability | |
CN112803377B (en) | Single-ended electric quantity protection method suitable for hybrid bipolar direct current transmission line | |
CN107271869B (en) | High-altitude area 750kV parallel reactor partial discharge test method | |
Li et al. | Probabilistic harmonic resonance assessment considering power system uncertainties | |
CN103544377A (en) | Calculating and analyzing method for resonance characteristics of direct current circuit of high-voltage direct current power transmission | |
Au et al. | Planning approaches for the strategic placement of passive harmonic filters in radial distribution networks | |
CN106226661A (en) | Use detection structure and the method for succusion detection reactor turn-to-turn insulation fault | |
CN105576662A (en) | Multiple-tuned filter design method for HVDC system | |
CN102185314A (en) | Power quality comprehensive evaluation method | |
CN110632394A (en) | Circuit impedance characteristic simulation system | |
Monem | Harmonic mitigation for power rectifier using passive filter combination | |
CN108983047A (en) | A kind of insulation imitative experimental appliance | |
CN110896218B (en) | Harmonic modeling method and system for establishing collective residential load | |
CN110336283B (en) | Method and system for analyzing influence of converter station alternating current filter group on power grid harmonic wave | |
Kovalenko et al. | Calculation of Resonant Modes for Power Supply Systems and Development of Measures on Higher Harmonics Filtering | |
CN110233476A (en) | Voltage stability assessment method and relevant apparatus during a kind of black starting-up | |
Ma et al. | Harmonie evaluation of grid with multiple harmonic sources based on DIgSILENT | |
Teng et al. | Resonance Risk in Infirm-interconnected Grid Considering Photovoltaic Devices | |
CN108347055A (en) | A kind of grid-connected filter inductor parameter evaluation circuit and its control method | |
CN107450024A (en) | Battery charging and discharging tester ripple absorption plant and method | |
CN114677031B (en) | Harmonic wave treatment demand evaluation method based on data envelope analysis | |
Zengyao et al. | Voltage fluctuation and ripple analysis method for DC grid connected PV converter |
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 |