CN106208041B - A kind of Distribution Network Harmonics current estimation method based on piecewise constant independent component analysis - Google Patents
A kind of Distribution Network Harmonics current estimation method based on piecewise constant independent component analysis Download PDFInfo
- Publication number
- CN106208041B CN106208041B CN201610564399.5A CN201610564399A CN106208041B CN 106208041 B CN106208041 B CN 106208041B CN 201610564399 A CN201610564399 A CN 201610564399A CN 106208041 B CN106208041 B CN 106208041B
- Authority
- CN
- China
- Prior art keywords
- harmonic
- current
- matrix
- section
- harmonic current
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The Distribution Network Harmonics current estimation method based on piecewise constant independent component analysis that the invention discloses a kind of.The piecewise constant model for initially setting up harmonic current estimation, determines harmonic current order using the correlation matrix correlation method of overlay segment, determines harmonic current using the method that minimizes the error of overlay segment, reduces the needs to priori harmonic current data.Then the criterion for providing electrical network parameter variation determines the data interval of electrical network parameter variation finally by dichotomy, corrects harmonic current estimated result in electrical network parameter variation, improves estimated accuracy when electrical network parameter variation.The method of the present invention can reduce harmonic current evaluated error when electrical network parameter variation, and estimation procedure does not need electrical network parameter data, and the harmonic current estimation being suitable under power distribution network dynamic state of parameters variation operating condition.
Description
Technical field
The Distribution Network Harmonics current estimation method based on piecewise constant independent component analysis that the present invention relates to a kind of, belongs to electricity
Force system field.
Background technique
In recent years, with the extensive application of power electronic equipment, harmonic pollution in electric power net aggravation has seriously affected the fortune of power grid
Row.Therefore, in order to provide a user good electric energy, guarantee the safety and stability of power distribution network, it is necessary to administer to harmonic pollution.
Installation filter is the important means of harmonic wave control, the position of installation and specification require according to harmonic injection levels of current into
Row is made rational planning for；The estimation of harmonic current injection level also provides accurate, real-time harmonic wave status number for distribution system simultaneously
According to, be power distribution network Situation Awareness significant data basis.
Existing harmonic current estimation method common are minimum mostly with reference to the method for Power system state estimation
Square law, singular value decomposition method, neural network etc..But these methods require complete electrical network parameter data, as power grid is opened up
Structure, harmonic impedance, method of operation etc. are flutterred, in the actual motion of power grid, is generally difficult to obtain complete electrical network parameter data,
The application of existing harmonic current estimation method is very limited.Meanwhile power distribution network changes of operating modes, opening and closing ring transformation,
Load such as cuts out at the variation that can all lead to electrical network parameter, to obtain accurate harmonic injection current estimation as a result, will
Adapt to the variation of electrical network parameter.And " He Zhengyou, Zang Tianlei, Yang Yuan, Sun Zhongmin one kind are based on quickly independent existing patent of invention
The harmonic source location method of Component Analysis and mutual information, the harmonic current estimation method in CN103412242A [P] .2013. " exist
Estimated accuracy is lower under the operating condition of electrical network parameter variation, and practical application is restricted.
Therefore, under conditions of electrical network parameter dynamic change, accurate harmonic injection current estimation is provided as a result, being current
A harmonic detecting important topic urgently to be resolved with identification field, is of great significance to the comprehensive treatment of harmonic wave.
Summary of the invention
Distribution Network Harmonics current estimation can be carried out under electrical network parameter change condition the object of the present invention is to provide a kind of
Method estimates model using piecewise constant, more accurately estimates harmonic injection electric current.
The present invention is to realize its goal of the invention, used technical solution are as follows:
A kind of Distribution Network Harmonics current estimation method based on piecewise constant independent component analysis, first to each section of harmonic wave electricity
It presses data to carry out piecewise constant independent component analysis, then is determined using the related coefficient correlation method of overlay segment with the method that minimizes the error
The order and amplitude of harmonic current finally obtain revised harmonic current estimated value, include following key step:
A, data acquire
Acquire the h subharmonic voltage data V=[V of all buses of power distribution network1,V2,…,Vi…,VN] and priori harmonic injection
Electric currentWherein Vi(1≤i≤N) be harmonic voltage on i-th bus to
Amount, 1 × L rank；(1≤i≤N) is the priori harmonic injection current vector on i-th bus, 1 × LpRank；L is total sampling
Points, LpFor harmonic voltage data segment, length, Lp≤ L, N are bus number.
B, the segmentation estimation of harmonic current
B1, harmonic current is estimated using piecewise constant independent component analysis method:
B11, to harmonic voltage data V, be segmented according to sampled point, every segment data length be Lp, weight between any two sections
Folded length is Lc, Lc≤Lp, m sections of harmonic voltage data V are obtained1,V2,…,Vi…,VN, m is according to L >=mLp+(m-1)LcIt solves
It obtains；
B12, to m sections of harmonic voltage data, establish piecewise constant Independent Component Analysis Model, Vk=AkIk, IkFor kth section
Harmonic injection electric current, AkFor the hybrid matrix of kth section.To kth section harmonic voltage dataIndependent component analysis calculating is carried out, seeks solving mixed matrix Wk, according toIt obtains
Estimate harmonic current Matrix W is mixed to solvingkIt inverts to obtain AkEstimation
Wherein(1≤i≤N) is harmonic voltage vector ViThe kth section of (1≤i≤N),(1≤i≤n) is that kth section estimates harmonic wave note
Enter current vector.
B2, related coefficient correlation method determine harmonic current order:
B21, as k=1, i.e., harmonic injection electric current first segment estimatedIt is infused using known priori harmonic wave
Enter electric currentBy seeking order matrix Fk, calculate the harmonic injection electric current for determining orderIt asksWithRelated coefficient
Matrix C1, whereinFor C1In i row j column element, indicateIn i-th of harmonic injection current vectorWithIn j-th
Harmonic injection current vectorRelated coefficient, find the maximum value of the i-th row element of correlation matrix C, which is located at
NiColumn, by order matrix FkThe row corresponding position element1 is set, remaining element zero setting:
According toAcquire the paragraph 1 harmonic injection electric current of determining order
B22, work as k > 1, i.e., to the estimation harmonic current other than first segmentOrder has been determined using -1 section of kth
Estimate harmonic currentBy seeking order matrix Fk, calculate the harmonic current for determining orderIt asksWithSampled point
LapWithCorrelation matrix C, wherein cijFor the i row j column element in C, indicateIn i-th harmonic wave electricity
Flow vectorWithIn j-th of harmonic current vectorRelated coefficient, find the maximum value of the i-th row element of C, the maximum
Value is located at NiColumn, by order matrix FkThe row corresponding position element1 is set, remaining element zero setting:
According toAcquire the kth section harmonic current of determining order
B3, harmonic current is determined using the method for minimizing the error:
B31, the kth section harmonic current as k=1, i.e., to determining orderUtilize known priori harmonic injection electric currentDetermine its amplitude.It is right respectivelyIn each harmonic current vectorAccording to formulaWith
Calculate proportionality coefficient akWith deviation ratio bk。
B32, as k > 1, i.e., the harmonic current data of order has been determined to other than first segment, has utilized -1 section of kth
The estimation harmonic current of order and amplitude has been determinedWith the kth section harmonic current for determining orderLapWithIt is rightWithCorrespondence vector according to formulaWithCalculate proportionality coefficient akWith deviation ratio bk；
B33, according to formulaAcquire the kth section harmonic current vector of determining amplitudeMerge
It obtains that the kth section harmonic current matrix of order and amplitude has been determinedIt willIt presses
Data label save toLater, the two lap takesEstimated value.
C, the amendment of the judgement of electrical network parameter variation and harmonic current estimated result
C1, determine whether electrical network parameter changes: to estimated mixing matrixEach row vector place is normalized
Reason,Obtain normalization estimated mixing matrixIt calculatesWithDeviation
Degree dk,WhereinFor kth section estimated mixing matrixIn the i-th row jth column member
Element,Hybrid matrix is normalized for kth sectionIn the i-th row jth column element.Work as dkWhen > λ, determine that electrical network parameter becomes
Change, wherein λ is constant threshold.
C2, dichotomy determine electrical network parameter delta data section:
Divided equally using data segment of the dichotomy to the electrical network parameter variation that previous step determines, repeats abovementioned steps, diminution is sentenced
The data interval of fixed electrical network parameter variation, obtains the data interval of electrical network parameter variation
The amendment of C3, harmonic current estimated result: the delta data section of judgement is removedHarmonic wave electricity
Estimated result is flowed, fromRestart segmentation estimation harmonic current.
Further, the concrete operations of harmonic current and hybrid matrix are asked in heretofore described step B12 is:
B121, pretreatment: with kth section harmonic voltage dataFor known quantity, to it
Go mean value and albefaction.To eachMean value is gone respectively, is enabledObtain the harmonic wave electricity of mean value
The amount of pressing toComposition removes the harmonic voltage matrix after mean value Correlation matrix Rx
Carry out Eigenvalues Decomposition Rx=Q Σ2QT, whitening matrix T=Σ-1QT, rightAlbefaction is carried out, after obtaining mean value and albefaction
Harmonic voltage data
B122, it extracts isolated component: by solving following optimization problem, obtaining solving mixed matrix Wk。
It maximizes:
Constraint:
Wherein, wi,wjIt is to solve mixed matrix WkColumn vector, E { } expression average,HIndicate to ask the hermitian of vector to turn
It sets,
Estimation harmonic current is obtained according to formulaMatrix W is mixed to solvingkIt inverts to obtain estimated mixing matrix
So merely with harmonic voltage data, while extracting its isolated component and estimating harmonic currentIt is with hybrid matrix
Electrical network parameter dataThe applicable elements of algorithm are relaxed significantly.
Further, in heretofore described step C2 dichotomy determine electrical network parameter delta data section specific behaviour
Work is:
C21, whether electrical network parameter, which changes, is judged to kth section estimated mixing matrix, if so, variation occurs in L1、L2、L3
Or L4Between.
C22, with L2、L3Midpoint is as harmonic data segment endpoint to be estimated, with LpFor data sectional length, LcTo be overlapped number
It according to length, reevaluates and determines, if it is determined that being, then constant interval is L1Or L2Section, on the contrary determine that constant interval is L3Or L4
Section.
C23, as shown in Figure 1, if upper step is judged as L1Or L2Section, then with L1、L2Midpoint is whole as harmonic data section to be estimated
Point continues to estimate and determine, if it is determined that being result L1Section, on the contrary result is L2Section；If upper step judges section for L3Or L4Section,
Determination method is similarly.
The electrical network parameter refers to the impedance class parameter of monitored power grid, although in the methods of the invention primarily directed to electricity
The resistance and reactance index of net, but also without loss of generality to the topological structure of electric parameter in addition to both.
Compared with prior art, the beneficial effects of the present invention are:
The invention firstly uses piecewise constant Independent Component Analysis Models, extract independent point to each section of harmonic voltage data
Amount, estimates harmonic current, does not need the topology or component parameters data of power grid.When determining harmonic current order and amplitude,
Using overlapping phase method, the dependence to priori harmonic current data is avoided, the independence of algorithm is improved.Finally by power grid
The judgement of Parameters variation and the amendment of estimated result improve the harmonic current estimated accuracy in electrical network parameter variation.This hair
Bright method can reduce the brings estimated bias such as grid operation mode variation, switch operation, and computational efficiency is higher, is applicable in
Change the harmonic injection current estimation under operating condition in electrical network parameter.
Detailed description of the invention
Fig. 1 is the schematic diagram that dichotomy determines electrical network parameter delta data section in C2 step.
Fig. 2 is artificial circuit schematic diagram of the invention.
Specific embodiment
The present invention is described in further detail with reference to the accompanying drawings and detailed description.
Embodiment
Fig. 2 shows a kind of Distribution Network Harmonics current segmenting situation of the invention, Distribution Network Harmonics current segmenting constant is estimated,
Using the step of are as follows:
A, data acquire
Acquire the h subharmonic voltage data V=[V of all buses of power distribution network1,V2,…,Vi…,VN] and priori harmonic injection
Current dataWherein Vi(1≤i≤N) is the harmonic wave on i-th bus
Voltage vector, 1 × L rank；(1≤i≤N) is the priori harmonic injection current vector on i-th bus, 1 × LpRank；L is
Total sampling number, LpFor harmonic voltage data segment, length (Lp≤ L), N is bus number.
B, the segmentation estimation of harmonic current
B1, harmonic current is estimated using piecewise constant independent component analysis method:
B11, to harmonic voltage data V, be segmented according to sampled point, every segment data length be Lp, weight between any two sections
Folded length is Lc(Lc≤Lp), m sections of harmonic voltage data V are obtained1,V2,…,Vi…,VN, m is according to L >=mLp+(m-1)LcIt solves
It obtains；
B12, to m sections of harmonic voltage data, establish piecewise constant Independent Component Analysis Model, Vk=AkIk, IkFor kth section
Harmonic injection electric current, AkFor the hybrid matrix of kth section.To kth section harmonic voltage data
Independent component analysis calculating is carried out, seeks solving mixed matrix Wk, according toObtain estimation harmonic current Matrix W is mixed to solvingkIt inverts to obtain AkEstimationWherein(1≤i≤N)
For harmonic voltage vector ViThe kth section of (1≤i≤N),(1≤i≤n) is that kth section estimates harmonic injection current vector.
B2, related coefficient correlation method determine harmonic current order:
B21, as k=1, i.e., harmonic injection electric current first segment estimatedIt is infused using known priori harmonic wave
Enter electric currentBy seeking order matrix Fk, calculate the harmonic injection electric current for determining orderSeek I0WithRelated coefficient
Matrix C, wherein cijFor the i row j column element in C, indicateIn i-th of harmonic injection current vectorWithIn j-th it is humorous
Wave Injection Current vectorRelated coefficient, find the maximum value of the i-th row element of correlation matrix C, which is located at the
NiColumn, by order matrix FkThe row corresponding position element1 is set, remaining element zero setting:
According toAcquire the paragraph 1 harmonic injection electric current of determining order
B22, work as k > 1, i.e., to the estimation harmonic current other than first segmentOrder has been determined using -1 section of kth
Estimate harmonic currentBy seeking order matrix Fk, calculate the harmonic current for determining orderIt asksWithSampled point
LapWithCorrelation matrix C, wherein cijFor the i row j column element in C, indicateIn i-th harmonic wave electricity
Flow vectorWithIn j-th of harmonic current vectorRelated coefficient, find the maximum value of the i-th row element of C, the maximum
Value is located at NiColumn, by order matrix FkThe row corresponding position element1 is set, remaining element zero setting:
According toAcquire the kth section harmonic current of determining order
B3, harmonic current is determined using the method for minimizing the error:
B31, as k=1, i.e., the estimation harmonic current of order has been determined to first segmentUtilize known priori harmonic wave
Injection CurrentDetermine its amplitude.It is right respectivelyIn each harmonic current to
Amount(1 < i < n), according to formulaWithCalculate proportionality coefficient akWith deviation ratio bk。
B32, as k > 1, i.e., the harmonic current data of order has been determined to other than first segment, has utilized -1 section of kth
The estimation harmonic current of order and amplitude has been determinedWith the kth section harmonic current for determining orderLapWithIt is rightWithCorrespondence vector according to formulaWithCalculate proportionality coefficient akWith deviation ratio bk；
B33, according to formulaAcquire the kth section harmonic current vector of determining amplitudeMerge
It obtains that the kth section harmonic current matrix of order and amplitude has been determinedIt willIt presses
Data label save toLater, the two lap takesEstimated value.
C, the amendment of the judgement of electrical network parameter variation and harmonic current estimated result
C1, determine whether electrical network parameter changes: to estimated mixing matrixEach row vector be normalized,Obtain normalization estimated mixing matrixIt calculatesWithDeviation journey
Spend dk,WhereinFor kth section estimated mixing matrixIn the i-th row jth column member
Element,Hybrid matrix is normalized for kth sectionIn the i-th row jth column element.Work as dkWhen > λ, determine that electrical network parameter becomes
Change, wherein λ is constant threshold.
C2, dichotomy determine electrical network parameter delta data section:
Divided equally using data segment of the dichotomy to the electrical network parameter variation that previous step determines, repeats abovementioned steps, diminution is sentenced
The data interval of fixed electrical network parameter variation, obtains the data interval of electrical network parameter variation
The amendment of C3, harmonic current estimated result: the delta data section of judgement is removedHarmonic wave electricity
Estimated result is flowed, fromRestart segmentation estimation harmonic current.
It is accurate under the variation of power distribution network dynamic state of parameters for a kind of harmonic current piecewise constant estimation method of the verifying present invention
Property, simulating, verifying is carried out on illustrating 34 node feeder systems, the voltage rating of system is 10kV, rated power 10MVA, list
Position harmonic current benchmark is 1000A.3 harmonic current sources are injected separately into node 11,29 and 33, are imitated by taking 5 subharmonic as an example
Very.
Collecting sample points N=1400, and the white Gaussian noise that signal-to-noise ratio is 120dB is added, line impedance is at the 900th
It changes after data point, using the method for the present invention, takes harmonic voltage data sectional length Lp=250, it is overlapped segment length Lc=
10, electrical network parameter changes decision threshold λ=0.05, judges that electrical network parameter constant interval process is as shown in table 1, estimates to harmonic current
Meter result is modified front and back and estimates harmonic currentWith practical harmonic current IkError it is as shown in table 2.
1 electrical network parameter constant interval deterministic process of table
2 harmonic current estimated result error of table
By data in table as it can be seen that the method for the present invention can effectively determine the variation of electrical network parameter, harmonic current estimated accuracy exists
It is improved after carrying out the judgement and the amendment of harmonic current estimated result of electrical network parameter variation, illustrates that the method for the present invention can adapt to
The variation of distribution network parameters reduces the harmonic current evaluated error under the operating condition.
Claims (3)
1. a kind of Distribution Network Harmonics current estimation method based on piecewise constant independent component analysis, first to each section of harmonic voltage
Data carry out piecewise constant independent component analysis, then using overlay segment related coefficient correlation method and minimize the error method determine it is humorous
The order and amplitude of wave electric current finally obtain revised harmonic current estimated value, include following key step:
A, data acquire
Acquire the h subharmonic voltage data V=[V of all buses of power distribution network1,V2,…,Vi…,VN] and priori harmonic injection electric currentWherein Vi(1≤i≤N) is the harmonic voltage vector on i-th bus, 1
× L rank；For the priori harmonic injection current vector on i-th bus, 1 × LpRank；L is total sampled point
Number, LpFor harmonic voltage data segment, length, Lp≤ L, N are bus number；
B, the segmentation estimation of harmonic current
B1, harmonic current is estimated using piecewise constant independent component analysis method:
B11, to harmonic voltage data V, be segmented according to sampled point, every segment data length be Lp, overlap length between any two sections
For Lc, Lc≤Lp, m sections of harmonic voltage data V are obtained1,V2,…,Vi…,VN, m is according to L >=mLp+(m-1)LcSolution obtains；
B12, to m sections of harmonic voltage data, establish piecewise constant Independent Component Analysis Model, Vk=AkIk, IkFor the harmonic wave of kth section
Injection Current, AkFor the hybrid matrix of kth section；To kth section harmonic voltage dataIt carries out
Independent component analysis calculates, and seeks solving mixed matrix Wk, according toObtain estimation harmonic current Matrix W is mixed to solvingkIt inverts to obtain AkEstimationWherein
For harmonic voltage vector ViThe kth section of (1≤i≤N),Harmonic injection current vector is estimated for kth section；
B2, related coefficient correlation method determine harmonic current order:
B21, as k=1, i.e., harmonic injection electric current first segment estimatedUtilize known priori harmonic injection electricity
StreamBy seeking order matrix Fk, calculate the harmonic injection electric current for determining orderIt asksWithCorrelation matrix
C1, whereinFor C1In i row j column element, indicateIn i-th of harmonic injection current vectorWithIn j-th of harmonic wave
Injection Current vectorRelated coefficient, find the maximum value of the i-th row element of correlation matrix C, which is located at Ni
Column, by order matrix FkThe row corresponding position element1 is set, remaining element zero setting:
According toAcquire the paragraph 1 harmonic injection electric current of determining order
B22, work as k > 1, i.e., to the estimation harmonic current other than first segmentThe estimation of order has been determined using -1 section of kth
Harmonic currentBy seeking order matrix Fk, calculate the harmonic current for determining orderIt asksWithSampled point overlapping
PartWithCorrelation matrix C, wherein cijFor the i row j column element in C, indicateIn i-th of harmonic current to
AmountWithIn j-th of harmonic current vectorRelated coefficient, find the maximum value of the i-th row element of C, the maximum value position
In NiColumn, by order matrix FkThe row corresponding position element1 is set, remaining element zero setting:
According toAcquire the kth section harmonic current of determining order
B3, harmonic current is determined using the method for minimizing the error:
B31, the kth section harmonic current as k=1, i.e., to determining orderUtilize known priori harmonic injection electric currentDetermine its amplitude；It is right respectivelyIn each harmonic current vectorAccording to formulaWithCalculate proportionality coefficient akWith deviation ratio bk；B32, as k > 1, i.e., pair
The harmonic current data of order has been determined other than first segment, has been determined that the estimation of order and amplitude is humorous using -1 section of kth
Wave electric currentWith the kth section harmonic current for determining orderLapWithIt is rightWithCorrespondence to
Amount is according to formulaWithCalculating ratio
Example coefficient akWith deviation ratio bk；
B33, according to formulaAcquire the kth section harmonic current vector of determining amplitudeMerging obtains really
The kth section harmonic current matrix of order and amplitude is determinedIt willBy data label
It saves extremelyLater, the two lap takesEstimated value；
C, the amendment of the judgement of electrical network parameter variation and harmonic current estimated result
C1, determine whether electrical network parameter changes: to estimated mixing matrixEach row vector be normalized,Obtain normalization estimated mixing matrixIt calculatesWithDeparture degree
dk,WhereinFor kth section estimated mixing matrixIn the i-th row jth column element,Hybrid matrix is normalized for kth sectionIn the i-th row jth column element；Work as dkWhen > λ, determine that electrical network parameter changes,
Wherein λ is constant threshold；
C2, dichotomy determine electrical network parameter delta data section:
Divided equally using data segment of the dichotomy to the electrical network parameter variation that C1 step determines, repeats A-C1 step, reduce judgement
The data interval of electrical network parameter variation, obtains the data interval of electrical network parameter variation
The amendment of C3, harmonic current estimated result: the delta data section of judgement is removedHarmonic current estimation
As a result, fromRestart segmentation estimation harmonic current.
2. Distribution Network Harmonics current estimation method according to claim 1, which is characterized in that
The concrete operations of harmonic current and hybrid matrix are asked in the step B12 is:
B121, pretreatment: with kth section harmonic voltage dataFor known quantity, it is gone
Value and albefaction；To eachMean value is gone respectively, is enabledObtain the harmonic voltage of mean value to
AmountComposition removes the harmonic voltage matrix after mean value Correlation matrix RxIt carries out
Eigenvalues Decomposition Rx=Q Σ2QT, whitening matrix T=Σ-1QT, rightCarry out albefaction, the harmonic wave after obtaining mean value and albefaction
Voltage data
B122, it extracts isolated component: by solving following optimization problem, obtaining solving mixed matrix Wk；
It maximizes:
Constraint:
Wherein, wi,wjIt is to solve mixed matrix WkColumn vector, E { } expression average,HThe hermitian transposition of vector is sought in expression,
Estimation harmonic current is obtained according to formulaMatrix W is mixed to solvingkIt inverts to obtain estimated mixing matrix
3. Distribution Network Harmonics current estimation method according to claim 1, which is characterized in that
Dichotomy determines that the concrete operations in electrical network parameter delta data section are in the step C2:
C21, whether electrical network parameter, which changes, is judged to kth section estimated mixing matrix, if so, variation occurs in L1、L2、L3Or L4It
Between；
C22, with L2、L3Midpoint is as harmonic data segment endpoint to be estimated, with LpFor data sectional length, LcIt is long for overlapped data
Degree, reevaluates and determines, if it is determined that being, then constant interval is L1Or L2Section, on the contrary determine that constant interval is L3Or L4Section；
If C23, upper step are judged as L1Or L2Section, then with L1、L2Continue to estimate and sentence as harmonic data segment endpoint to be estimated in midpoint
It is fixed, if it is determined that being result L1Section, on the contrary result is L2Section；If upper step judges section for L3Or L4Section, determination method is similarly.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610564399.5A CN106208041B (en) | 2016-07-18 | 2016-07-18 | A kind of Distribution Network Harmonics current estimation method based on piecewise constant independent component analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610564399.5A CN106208041B (en) | 2016-07-18 | 2016-07-18 | A kind of Distribution Network Harmonics current estimation method based on piecewise constant independent component analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106208041A CN106208041A (en) | 2016-12-07 |
| CN106208041B true CN106208041B (en) | 2019-03-01 |
Family
ID=57492956
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610564399.5A Active CN106208041B (en) | 2016-07-18 | 2016-07-18 | A kind of Distribution Network Harmonics current estimation method based on piecewise constant independent component analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106208041B (en) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106980044B (en) * | 2017-03-22 | 2019-08-20 | 西南交通大学 | A kind of Harmonious Waves in Power Systems current estimation method adapting to wind power integration |
| CN107576866B (en) * | 2017-08-22 | 2019-08-20 | 西南交通大学 | It is a kind of to approach the sparse method for distinguishing multiple harmonic sources reconstructed with interior point method based on smooth |
| CN107991532B (en) * | 2017-11-21 | 2020-07-24 | 厦门理工学院 | Harmonic threshold value merging method based on multiple operation modes |
| CN110896218B (en) * | 2019-11-08 | 2021-09-03 | 山东大学 | Harmonic modeling method and system for establishing collective residential load |
| CN112834821B (en) * | 2020-12-31 | 2023-08-11 | 国网上海市电力公司 | Multi-DC feed-in full-cable power supply area multi-harmonic source harmonic responsibility quantification method |
| CN112858782B (en) * | 2021-01-07 | 2022-04-26 | 国网河南省电力公司电力科学研究院 | Harmonic responsibility quantification method and system under influence of multi-user interaction of power system |
| CN114660362B (en) * | 2022-05-19 | 2022-08-16 | 四川大学 | System side harmonic impedance estimation method and device based on mutual information data optimization |
| CN115453193B (en) * | 2022-09-15 | 2023-04-18 | 四川大学 | Power distribution network harmonic state estimation method based on cooperation of PQM, TTU and SM measurement data |
| CN116388190B (en) * | 2023-06-06 | 2023-08-18 | 山东大学 | Distributed photovoltaic harmonic quantitative evaluation method, system, terminal equipment and medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103412242A (en) * | 2013-08-01 | 2013-11-27 | 西南交通大学 | Method for locating harmonic source based on independent rapid component analysis and mutual information |
| CN105071387A (en) * | 2015-08-13 | 2015-11-18 | 河海大学 | Power distribution network rapid decomposition state estimation method based on complex number domain per unit theory |
| CN105548809A (en) * | 2015-12-18 | 2016-05-04 | 南瑞（武汉）电气设备与工程能效测评中心 | Harmonic wave responsibility determining method of power system |
| CN105548718A (en) * | 2015-12-14 | 2016-05-04 | 南瑞（武汉）电气设备与工程能效测评中心 | System harmonic impedance calculation method based on mixed overall least square method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101661069B (en) * | 2009-09-25 | 2011-07-20 | 北京四方继保自动化股份有限公司 | Dynamic process real-time estimation method of weak observable non-PMU measuring point independent of state matrix |
-
2016
- 2016-07-18 CN CN201610564399.5A patent/CN106208041B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103412242A (en) * | 2013-08-01 | 2013-11-27 | 西南交通大学 | Method for locating harmonic source based on independent rapid component analysis and mutual information |
| CN105071387A (en) * | 2015-08-13 | 2015-11-18 | 河海大学 | Power distribution network rapid decomposition state estimation method based on complex number domain per unit theory |
| CN105548718A (en) * | 2015-12-14 | 2016-05-04 | 南瑞（武汉）电气设备与工程能效测评中心 | System harmonic impedance calculation method based on mixed overall least square method |
| CN105548809A (en) * | 2015-12-18 | 2016-05-04 | 南瑞（武汉）电气设备与工程能效测评中心 | Harmonic wave responsibility determining method of power system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106208041A (en) | 2016-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106208041B (en) | A kind of Distribution Network Harmonics current estimation method based on piecewise constant independent component analysis | |
| CN103454559B (en) | A kind of one-phase earthing failure in electric distribution network Section Location and locating device | |
| Hu et al. | Measurement-based real-time voltage stability monitoring for load areas | |
| CN111428355B (en) | Modeling method for intelligent comprehensive statistics of power load numbers | |
| CN107843810B (en) | Active power distribution network fault section online positioning method based on state estimation | |
| CN103872681A (en) | Online real-time loop closing method based on integration of major network and distribution network | |
| CN107453357A (en) | A kind of State Estimation for Distribution Network based on hierarchical solving | |
| CN107944612B (en) | Bus net load prediction method based on ARIMA and phase space reconstruction SVR | |
| CN106208049A (en) | The practical approach that a kind of power distribution network simple state is estimated | |
| CN109767054A (en) | Efficiency cloud appraisal procedure and edge efficiency gateway based on deep neural network algorithm | |
| CN106443244A (en) | Method and system for identifying type of electric devices | |
| CN106295160B (en) | AC-DC interconnecting power network Thevenin's equivalence parameter on-line calculation method | |
| CN110289613A (en) | The identification of distribution net topology and line parameter circuit value discrimination method based on sensitivity matrix | |
| CN107947169B (en) | Information flow modeling method for power grid energy management system | |
| CN110797874A (en) | State estimation method for alternating current-direct current hybrid power distribution network containing power electronic transformer | |
| CN111965484B (en) | Power distribution network harmonic contribution calculation method and system based on continuous harmonic state estimation | |
| CN109524982A (en) | A kind of alternating current-direct current Power Network Transient Stability methods of risk assessment | |
| CN106877309B (en) | Stability judgment method for direct-current micro-grid in master-slave control mode | |
| CN106383280B (en) | The voltage transformer Model test Method of line model is maked somebody a mere figurehead based on two nodes | |
| CN109271925A (en) | A kind of electrical equipment kind identification method on intelligent socket | |
| Yi et al. | Phase identification of low-voltage distribution network based on stepwise regression method | |
| CN108233364A (en) | Network equivalence modeling method based on boundary node interconnection structure | |
| Hu et al. | A data-driven method of users-transformer relationship identification in the secondary power distribution system | |
| CN115455802A (en) | Low-voltage distribution network modeling based on data driving and visualization method thereof | |
| CN113408617A (en) | XGboost and Stacking model fusion-based non-invasive load identification method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |