CN107590604B - Coherent unit grouping method and system combining S transformation and 2DPCA - Google Patents
Coherent unit grouping method and system combining S transformation and 2DPCA Download PDFInfo
- Publication number
- CN107590604B CN107590604B CN201710821583.8A CN201710821583A CN107590604B CN 107590604 B CN107590604 B CN 107590604B CN 201710821583 A CN201710821583 A CN 201710821583A CN 107590604 B CN107590604 B CN 107590604B
- Authority
- CN
- China
- Prior art keywords
- clustering
- neuron
- matrix
- sample
- 2dpca
- 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
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000001427 coherent effect Effects 0.000 title claims abstract description 40
- 230000009466 transformation Effects 0.000 title claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims abstract description 78
- 238000013528 artificial neural network Methods 0.000 claims abstract description 14
- 210000002569 neuron Anatomy 0.000 claims description 45
- 210000004205 output neuron Anatomy 0.000 claims description 25
- 238000005259 measurement Methods 0.000 claims description 16
- 239000006185 dispersion Substances 0.000 claims description 15
- 230000001360 synchronised effect Effects 0.000 claims description 12
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 210000002364 input neuron Anatomy 0.000 claims description 9
- 230000009467 reduction Effects 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims 4
- 238000005192 partition Methods 0.000 claims 2
- 238000000926 separation method Methods 0.000 claims 1
- 230000010355 oscillation Effects 0.000 abstract description 3
- 238000000605 extraction Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 206010033799 Paralysis Diseases 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
The invention relates to a coherent unit grouping method and system combining S transformation and 2 DPCA. The method comprises the steps of obtaining a real-time power angle of a system unit by using a WAMS phasor measuring device, converting a power angle signal of each generator into a time-frequency characteristic module value matrix by adopting rapid S transformation, reducing the dimension of the time-frequency characteristic module value matrix by using 2DPCA, converting the time-frequency characteristic module value matrix into a low-dimension characteristic index matrix, and inputting the characteristic index matrix into a self-organizing neural network for clustering identification. The method can divide the units of the system on line in real time along with different faults and different fault places under different operation modes, comprehensively considers the time-frequency domain information of the power angle, has high identification accuracy and stable clustering result, provides necessary premise for simplifying the power grid, determining the system oscillation center and controlling the splitting, and ensures the safe, efficient and stable operation of the power grid.
Description
Technical Field
The invention relates to the technical field of coherent unit identification, in particular to a coherent unit grouping method and system combining S transformation and 2 DPCA.
Background
With the implementation of the engineering project of 'west-east power transmission, south-north interconnection', national power grids are developing towards large-area interconnection power grids. However, while economic benefits are improved, the problem of grid safety is more and more prominent, and small disturbance may cause cascading failure, resulting in large-area power failure. The active splitting control is the last ring for avoiding the power grid paralysis, and the searching of the splitting section needs to be based on accurate unit coherent grouping. In addition, along with the expansion of the scale of the power grid, the system has various fault modes, and the system equivalence is a single-machine infinite system or a two-machine system during fault analysis, so that the requirement of transient stability analysis cannot be completely met. The dynamic equivalence method is a good equivalence simplification method of a system model, is beneficial to quickly analyzing the dynamic stability of the system and provides reference for active splitting, and the unit coherent grouping is the core technology of the dynamic equivalence method. Therefore, the quick and accurate unit coherent grouping has important significance.
When a power system fails, some generators in the system exhibit similar dynamic characteristics, referred to as unit coherence. The coherent grouping method of the units mainly comprises three types: identifying a coherent machine group according to a maneuvering parameter of a generator rotor; analyzing the unit coherence characteristics according to the disturbance system state matrix trajectory characteristic root; and thirdly, performing unit homodyne identification according to the time-frequency characteristics of the power angle swing curve of the disturbance generator. The first method needs to define the parameters of the system generator in advance, the accuracy of the unit coherent grouping depends on the accuracy of the generator parameters and the like to a great extent, but the accurate generator parameters are not easy to obtain sometimes. In addition, the method cannot consider the influence of the system operation condition and the accident type on the grouping result of the unit. In the second method, before the unit coherence identification, the power grid system needs to be simplified in a linearization manner, which may result in low accuracy of the unit clustering result and is not suitable for a large-scale interconnection system. The third method is originally restricted by obtaining the power angle information of the generator in real time, but with the development of a wide area measurement system and an optical fiber communication technology, the real-time power angle information of the generator can be easily obtained, and then the coherent cluster is identified according to the effective time-frequency characteristics of the power angle information. The method not only can realize real-time grouping by considering the system operation condition and the fault type, but also is suitable for a large-scale interconnection system.
And identifying the coherent cluster according to the power angle information of the generator, wherein effective characteristics of the power angle need to be extracted. At present, the power angle information feature extraction method comprises the following steps: fourier transform, wavelet transform, empirical mode decomposition, hilbert-yellow transform, and the like. However, these methods have disadvantages, such as poor effect of fourier transform on unstable signal processing, difficulty in selecting wavelet basis functions in wavelet transform, modal aliasing phenomenon in inherent modal components in empirical mode decomposition, and complex recursive operation required for hilbert-yellow transform.
Disclosure of Invention
The invention aims to provide a coherent unit grouping method and system combining S transformation and 2DPCA, so as to overcome the defects of complex operation, low efficiency and accuracy and the like in the prior art.
In order to achieve the purpose, the technical scheme of the invention is as follows: a coherent unit grouping method combining S transformation and 2DPCA,
the method comprises the following steps:
step S1: acquiring power angle information increment after the M generators have faults through a synchronous phasor measuring device;
step S2: obtaining corresponding power angle rocking curve information according to the power angle information increment of each generator obtained in the step S1, and performing rapid S transformation to obtain an mxn-dimensional time-frequency domain information matrix;
step S3: reducing the dimension of the m × n-dimensional time-frequency domain information matrix of the step S2 by adopting a 2DPCA algorithm to obtain a q × p-dimensional characteristic index matrix Ts; wherein s is a generator label;
step S4: and clustering M samples in the characteristic index matrix obtained in the step S3 to obtain a coherent unit.
In an embodiment of the present invention, in the step S1, the increment dimension of the power angle information is N, and an nxm-dimensional power angle information database is obtained.
In an embodiment of the present invention, in the step S4, the M samples in the characteristic index matrix obtained in the step S3 are clustered by using a self-organizing neural network.
In an embodiment of the present invention, the clustering by using the self-organizing neural network is implemented according to the following steps:
step S41: the number of input neurons is the number of elements of the characteristic index matrix Ts, and the input neurons are used as a sample characteristic index matrix; the output layer neuron is initially one; will connect the weight WijInitially setting an arbitrary sample characteristic index matrix, wherein i is an input layer neuron label, and j is an output layer neuron label;
step S42: recording a division precision threshold used for determining the accuracy degree of clustering as lambda;
step S43: inputting a new sample X, calculating Euclidean distance between the new sample and a corresponding clustering center of each output neuron, and recording the minimum distance as DjJ is the output layer neuron label corresponding to the minimum distance, and the formula is as follows:
in the formula: xiIs an input vector; connection weight WijA cluster center of a jth class mode;
step S44: if D isj<λ, then the jth output neuron wins, the sample is taken into the jth neuron, and the neuron connection weight is updated as follows:
Wij(n+1)=Wij(n)+η(n)(Xi-Wij(n))
wherein η (n) is the reciprocal of the total number of samples contained in the neuron;
if D isj>Lambda, the sample does not belong to any existing output neuron, an output neuron is newly added, and the initial connection weight of the output neuron and the input layer is a sample characteristic index matrix;
step S45: and sequentially inputting all samples according to the steps, wherein the number of final output layer neurons is the number of clustering clusters, each output layer neuron comprises a sample which is a same cluster sample, and the connection weight of the sample and the input layer is the clustering center of the cluster.
In an embodiment of the present invention, the method further includes a step S5: and introducing an internal effective index-CH index to determine the optimal cluster number.
In one embodiment of the present invention, in the step S5, the internal valid indicator — CH indicator is based on a measure of intra-class dispersion matrix and inter-class dispersion matrix of all samples, as follows:
wherein k is the number of clusters, n is the number of all samples, tr (B (k)) is the trace of the inter-class dispersion matrix, and tr (W (k)) is the trace of the intra-class dispersion matrix; and k corresponding to the maximum value obtained by the CH index is the optimal cluster number.
In an embodiment of the present invention, the method further includes a step S6: the synchrophasor measurement device acquires the power angle information increment of the M generators after the fault once again every certain time, refreshes the power angle information database, and repeats the steps S2, S3, and S4.
Further, a coherent unit grouping system combining S transformation and 2DPCA is also provided, including:
the synchronous phasor measurement unit is used for acquiring power angle information increment after the M generators are in fault through the synchronous phasor measurement device;
the S transformation unit is used for carrying out rapid S transformation on the power angle swing curve information of each generator acquired by the synchronous phasor measurement unit and acquiring an m multiplied by n dimensional time-frequency domain information matrix;
a dimension reduction unit based on a 2DPCA algorithm, which is used for reducing the dimension of the m multiplied by n dimension time-frequency domain information matrix through the 2DPCA algorithm to obtain a q multiplied by p dimension characteristic index matrix Ts, wherein s is the generator label;
and the clustering unit is used for clustering the characteristic index matrix Ts to obtain the coherent unit.
In an embodiment of the present invention, the method further includes: and an optimal clustering number calculating unit for determining the optimal clustering number according to an internal effective index-CH index.
Compared with the prior art, the invention has the following beneficial effects: the invention provides a coherent unit grouping method and a coherent unit grouping system combining S transformation and 2DPCA, which aim at the problems of feature extraction one-sided and complex and mutual restriction of calculation processes faced by coherent group identification. The multi-resolution time-frequency information of the power angle swing curve of the generator is obtained by utilizing the rapid S transformation and the 2DPCA, the characteristics are comprehensively extracted, and effective characteristics can be effectively mined. And the self-organizing neural network is used as a classifier to perform cluster analysis on the generator characteristic index matrix, so that a cluster result is stable. The method can acquire the grouping state of the system unit in real time, provide necessary premises for simplifying the power grid, determining the system oscillation center and controlling the splitting, and ensure the safe, efficient and stable operation of the power grid; the invention is easy to operate and suitable for practical engineering application.
Drawings
Fig. 1 is a flow chart of a coherent unit identification method using fast S-transform and 2DPCA according to the present invention.
FIG. 2 is a graph illustrating power angle incremental swing of an IEEE-39 generator in accordance with an embodiment of the present invention.
Fig. 3 is a schematic diagram illustrating the information content of the index matrix of different dimensional characteristics according to an embodiment of the present invention.
Fig. 4 is a graph comparing the trend of the power angle swing curve of the generator according to an embodiment of the present invention.
Detailed Description
The technical scheme of the invention is specifically explained below with reference to the accompanying drawings.
The invention provides a coherent unit grouping method combining S transformation and 2DPCA, and the invention is explained in detail with the accompanying drawings and the specific embodiments.
As shown in fig. 1, a flow chart of a coherent unit identification method combining fast S-transform and 2DPCA includes the following steps:
step S1: and obtaining the power angle information increment after the M generators have faults by the synchronous phasor measuring device, wherein the dimensionality is N, and obtaining an NxM dimensional power angle information database.
Step S2: and obtaining corresponding power angle swing curve information according to the power angle information increment of each generator obtained in the step S1, and performing S transformation to obtain an m × n-dimensional time-frequency domain information matrix.
S3, adopting 2DPCA algorithm to reduce the dimension of the time-frequency domain information matrix of the step S2, and finally obtaining a q × p-dimension characteristic index matrix TsAnd s is the generator number.
Step S4: and clustering the M samples obtained in the step S3, namely the characteristic index matrix, by using a self-organizing neural network, wherein the samples clustered into the same cluster have high similarity, and the change rules of the corresponding power angle swing curves are similar, so that the unit coherence is obtained.
Step S5: and introducing an internal effective index-CH index to determine the optimal cluster number.
Step S6: and the PMU (phasor measurement Unit) acquires the power angle data of the system unit once again every specific time, refreshes the database, and repeats the steps S2, S3 and S4. Preferably, the specific time period is 10 ms.
Furthermore, in this embodiment, the adopted S transform can perform multi-resolution analysis on the unstable signal, and the width of the gaussian window is adjustable, so that the method has good time-frequency transformation capability. The power angle data information is converted into a high-dimensional matrix through S conversion, and the matrix contains a large amount of redundant information. Yang et al propose two-dimensional principal component analysis (2 DPCA for short) can extract the main information characteristic of two-dimensional matrix, have high, advantage such as being simple and easy of accuracy with high efficiency.
Further, in this embodiment, the internal valid indicator — CH indicator is based on the measure of the intra-class dispersion matrix and the inter-class dispersion matrix of all samples, and is defined by the following formula:
wherein k is the number of clusters, n is the total number of samples, tr (B (k)) is the trace of the inter-class dispersion matrix, and tr (W (k)) is the trace of the intra-class dispersion matrix. And k corresponding to the maximum value obtained by the CH index is the optimal cluster number.
Further, in this embodiment, the specific steps of using the self-organizing neural network algorithm are as follows:
step S41: and (5) initializing the network. The number of input neurons is the number of elements of the characteristic index matrix. The output layer neurons are initially one. Connection weight Wij(i is the input layer neuron label, j is the output layer neuron label) is initially an arbitrary sample characteristic index matrix.
Step S42: a division accuracy threshold λ is set. The accuracy degree of clustering is determined by a threshold lambda, and the smaller the threshold lambda is, the more accurate the clustering result is and the more the classification is.
Step S43: and (5) identifying clusters. Inputting a new sample X, calculating Euclidean distance between the new sample and a corresponding clustering center of each output neuron, and recording the minimum distance as DjAnd j is the output layer neuron label corresponding to the minimum distance. The formula is as follows:
in the formula: xiIs an input vector; connection weight WijThe cluster center of the jth class pattern.
Step S44: and updating the network. If D isj<λ, then the jth output neuron wins, the sample is taken into the jth neuron, and the neuron connection weight is updated as follows:
Wij(n+1)=Wij(n)+η(n)(Xi-Wij(n)) (3)
η (n) is the reciprocal of the total number of samples contained in the neuron.
If D isj>And lambda, if the sample does not belong to any existing output neuron, adding a new output neuron, and taking the initial connection weight of the new output neuron and the input layer as a sample characteristic index matrix.
Step S45: and (5) performing cyclic training. And sequentially inputting all samples according to the steps, wherein the number of final output layer neurons is the number of clustering clusters, each output layer neuron comprises a sample which is a same cluster sample, and the connection weight of the sample and the input layer is the clustering center of the cluster.
Further, in this embodiment, a coherent group grouping system combining S transform and 2DPCA is further provided, including:
the synchronous phasor measurement unit is used for acquiring power angle information increment after the M generators are in fault through the synchronous phasor measurement device;
the S transformation unit is used for carrying out rapid S transformation on the power angle swing curve information of each generator acquired by the synchronous phasor measurement unit and acquiring an m multiplied by n dimensional time-frequency domain information matrix;
a dimension reduction unit based on a 2DPCA algorithm, which is used for reducing the dimension of the m multiplied by n dimension time-frequency domain information matrix through the 2DPCA algorithm to obtain a q multiplied by p dimension characteristic index matrix Ts, wherein s is the generator label;
and the clustering unit is used for clustering the characteristic index matrix Ts to obtain the coherent unit.
In this embodiment, the synchrophasor measurement unit obtains, through the synchrophasor measurement device, power angle information increments after failures of M generators, where the dimensionality is N, and obtains an nxm-dimensional power angle information database. After the grouping of the coherent units is completed, the synchrophasor measurement unit acquires the power angle information increment of the M generators after the fault once again every 10ms, refreshes the power angle information database, and repeats the step S2, the step S3 and the step S4.
And the clustering unit clusters and groups the M samples in the characteristic index matrix by adopting a self-organizing neural network to obtain a coherent unit. The self-organizing neural network adopted for clustering is realized according to the following steps:
step A: the number of input neurons is the number of elements of the characteristic index matrix Ts, and the input neurons are used as a sample characteristic index matrix; the output layer neuron is initially one; will connect the weight WijInitially setting an arbitrary sample characteristic index matrix, wherein i is an input layer neuron label, and j is an output layer neuron label;
and B: recording a division precision threshold used for determining the accuracy degree of clustering as lambda;
and C: inputting a new sample X, calculating Euclidean distance between the new sample and a corresponding clustering center of each output neuron, and recording the minimum distance as DjJ is the output layer neuron label corresponding to the minimum distance, and the formula is as follows:
in the formula: xiIs an input vector; connection weight WijA cluster center of a jth class mode;
step D: if D isj<λ, then the jth output neuron wins, the sample is taken into the jth neuron, and the neuron connection weight is updated as follows:
Wij(n+1)=Wij(n)+η(n)(Xi-Wij(n))
wherein η (n) is the reciprocal of the total number of samples contained in the neuron;
if D isj>Lambda, the sample does not belong to any existing output neuron, an output neuron is newly added, and the initial connection weight of the output neuron and the input layer is a sample characteristic index matrix;
step E: and sequentially inputting all samples according to the steps, wherein the number of final output layer neurons is the number of clustering clusters, each output layer neuron comprises a sample which is a same cluster sample, and the connection weight of the sample and the input layer is the clustering center of the cluster.
Further, in this embodiment, the method further includes: and an optimal clustering number calculating unit for determining the optimal clustering number according to an internal effective index-CH index. The internal effective index, the CH index, is based on the measure of the intra-class dispersion matrix and the inter-class dispersion matrix of all samples, as follows:
wherein k is the number of clusters, n is the number of all samples, tr (B (k)) is the trace of the inter-class dispersion matrix, and tr (W (k)) is the trace of the intra-class dispersion matrix; and k corresponding to the maximum value obtained by the CH index is the optimal cluster number.
Furthermore, in order to make those skilled in the art further understand the method proposed in the present invention, the IEEE-39 classic system is used as an example, and G is used as an example30The generator is used as a reference, a three-phase open circuit fault occurs on a line bus3-bus18, the fault is removed after the fault lasts for 0.14s, the front and the back of a power grid structure are kept consistent, a power angle change curve of a system generator at a time period of 4.0s from the fault occurrence time is monitored, simulated power angle data are derived from PSDEdit, power angle increment data are obtained through Matlab processing, a power angle swing curve is drawn, and the power angle increment swing curve of an IEEE-39 generator is shown in an attached figure 2.
Furthermore, rapid S conversion is carried out on the power angle increment time sequence number sequence of each generator, the sampling frequency is 100Hz, the total number of samples is 400, the corresponding frequency range of S conversion time frequency information matrix rows is 0 Hz-50 Hz, the interval is 0.25Hz, the corresponding time range of the columns is 0S-4S, and the interval is 0.01S. And 2DPCA dimension reduction processing is carried out on the time-frequency information matrix after S transformation, when the main component numbers p and q of the rows and the columns take different values respectively, the information content of the characteristic index matrix is different, and the information content of the characteristic index matrix with different dimensions is shown in an attached figure 3.
According to analysis, when p is 2 and q is 1 in 2DPCA, the principal component already contains 81% of the information content of the time-frequency information matrix, so that the self-organizing neural network input sample is selected as a 2 × 1 dimensional characteristic index matrix, the samples are clustered, the number of clusters changes along with the change of the threshold value λ, and the clustering result is shown in table 1.
TABLE 1 ST-2DPCA coherent unit grouping results
Number of classification | ST-2DPCA clustering results |
9 | G31/G32/G33/G34/G35/G36/G37/G38/G39 |
8 | G31/G32/G33/G34/G35G36/G37/G38/G39 |
7 | G31/G32/G33G34/G35G36/G37/G38/G39 |
6 | G31/G32/G33G34G35G36/G37/G38/G39 |
5 | G31/G32/G33G34G35G36/G37/G38G39 |
4 | G31G37/G32/G33G34G35G36/G38G39 |
3 | G31G37/G32G38G39/G33G34G35G36 |
2 | G31G32G37G38G39/G33G34G35G36 |
1 | G31G32G33G34G35G36G37G38G39 |
The CH index varies with the number of clusters, as shown in table 2. When k is 3, the CH index takes the maximum value, and therefore the optimal cluster number is 3.
TABLE 2 index values for different cluster numbers CH
Number of clusters clustered | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
CH index | 89.7 | 119.6 | 110.2 | 91.2 | 62.6 | 32.9 | 7.9 |
After the failure, the database is refreshed every 10ms, the grouping result is continuously changed along with the updating of the database,
when the units are divided into 3 groups, the change process of the grouping result at different times is shown in table 3:
TABLE 3 comparison of clustering results at different times
Time t(s) | Clustering results |
0.2 | G32/G31G37G38G39/G33G34G35G36 |
0.5 | G31G32G38G39/G33G34G35G36/G37 |
1.0 | G31G37/G32G38G39/G33G34G35G36 |
1.5 | G31G37/G32G38G39/G33G34G35G36 |
2.0 | G31G37/G32G38G39/G33G34G35G36 |
2.5 | G31G37/G32G38G39/G33G34G35G36 |
It can be seen from the table that after the fault occurs for 1s, the clustering result tends to be stable, because the power angle increment information in the database contains more and more power angle increment information along with the increase of time, and the actual characteristics of the unit can be reflected more and more.
The characteristic indexes extracted by the ST-2DPCA method and the characteristic indexes extracted by the PCA, ICA and NMF methods provided by the text are respectively input into the self-organizing neural network for clustering, and the optimal clustering result obtained according to the CH index is shown in Table 4.
TABLE 4 comparison of coherent unit grouping results
Feature extraction method | Clustering results |
ST-2DPCA | G31G37/G32G38G39/G33G34G35G36 |
PCA | G31G32G37/G33G34G35G36/G38G39 |
ICA | G31G32G37/G33G34G35G36/G38G39 |
NMF | G31G32/G33G34G35G36/G37G38G39 |
Defining a unit coherence criterion formula as follows:
k may take any constant such that SijMinimum, let:
the following can be found:
and (5) verifying the effect of the 4 clustering results by utilizing a coherence criterion formula. The maximum criterion among all clusters of samples in the ST-2DPCA feature extraction method grouping result is 0.6721, and the maximum criterion among all clusters of samples in the PCA, ICA and NMF feature extraction method grouping result is 0.9652, 0.9652 and 1.2745 respectively.
The four unit grouping methods all group the generators {33,34,35,36}, and it can be seen from fig. 4(2) that the power angle swing curves are similar and have strong coherence. As can be seen from FIG. 4(1) and FIG. 4(3), the oscillation frequencies of the generators {32,38,39} are similar, and the homologies of the generators {32} and {38,39} are stronger than those of the generators {32} and {31,37 }. Looking at fig. 4(4), it is clear that the generators {37} are not in tune with the generators {38,39 }.
Furthermore, the criterion and the trend of the power angle curve all show that the feature extraction method in the embodiment can effectively extract the time-frequency features of the power angle signals compared with the PCA, ICA and NMF methods. The method provided by the invention is verified to have higher accuracy.
Further, in this embodiment, a 30 # generator is used as a reference, different fault types are set at different positions, the fault is removed after the duration lasts for 0.12s, the power angle of the system generator after the fault starts is monitored, the coherent unit clustering is performed according to the method, and the clustering result is shown in table 5.
TABLE 5 comparison of clustering results for different accident types
Type of accident | Clustering results |
Three-phase open circuit of line 4-14 | G31G37/G33G34G35G36/G32G38G39 |
Line 25-26 three-phase short circuit | G31G32/G33G34G35G36/G37/G38G39 |
Three-phase short circuit of lines 16-24 | G31G32/G33G34G35G36/G37G38G39 |
Three-phase break of line 9-39 | G31G32/G33G35G36G37/G34G38/G39 |
The grouping results in table 5 show that the method can perform unit grouping on the system in real time along with the change of different faults and different fault locations in different operation modes, and the model-based grouping method (performing unit grouping according to system parameters) cannot realize the function.
The above are preferred embodiments of the present invention, and all changes made according to the technical scheme of the present invention that produce functional effects do not exceed the scope of the technical scheme of the present invention belong to the protection scope of the present invention.
Claims (7)
1. A coherent unit grouping method combining S transformation and 2DPCA is characterized by comprising the following steps:
step S1: acquiring power angle information increment after the M generators have faults through a synchronous phasor measuring device;
step S2: obtaining corresponding power angle swing curve information according to the power angle information increment of each generator obtained in the step S1, and performing rapid S transformation to obtain an m × n-dimensional time-frequency domain information matrix;
step S3: for said step S2 by using 2DPCA algorithmm×nReducing the dimension of the dimension time-frequency domain information matrix to obtainq ×pFeature index matrix of dimensionTs(ii) a Wherein the content of the first and second substances,sthe generator is labeled;
step S4: clustering M samples in the characteristic index matrix obtained in the step S3 by using a self-organizing neural network to obtain a coherent unit;
the clustering by adopting the self-organizing neural network is realized according to the following steps:
step S41: the number of input neurons is the number of elements of the characteristic index matrix Ts, and the input neurons are used as a sample characteristic index matrix; the output layer neuron is initially one; will connect the weightW ij Initially, an arbitrary sample characteristic index matrix is obtained,ifor the input layer neuron labels,jis the output layer neuron label;
step S42: the partition accuracy threshold, which is used to determine the accuracy of the clustering, isλ;
Step S43: inputting a new sample X, calculating Euclidean distance between the new sample and a corresponding clustering center of each output neuron, and recording the minimum distance asD j ,jFor the output layer neuron label corresponding to the minimum distance, the formula is as follows:
in the formula:X i is an input vector; connection weightW ij Is as followsjA cluster center of class group patterns;
step S44: if it isD j <λ, thenjAn output neuron wins, the sample is includedjIn each neuron, updating the neuron connection weight value, wherein the formula is as follows:
wherein the content of the first and second substances,η(n)is the reciprocal of the total number of samples contained in the neuron;
if it isD j >λIf the sample does not belong to any existing output neuron, adding an output neuron, and taking the initial connection weight of the output neuron and the input layer as a sample characteristic index matrix;
step S45: and sequentially inputting all samples according to the steps, wherein the number of final output layer neurons is the number of clustering clusters, each output layer neuron comprises a sample which is a same cluster sample, and the connection weight of the sample and the input layer is the clustering center of the cluster.
2. The coherent unit clustering method combining S-transform and 2DPCA according to claim 1, wherein in step S1, the power angle information increment dimension is N, and a power angle information database with N × M dimensions is obtained.
3. The coherent group clustering method combining the S-transform and the 2DPCA according to claim 1, further comprising a step S5: and introducing an internal effective index-CH index to determine the optimal cluster number.
4. The coherent group clustering method combining S-transform and 2DPCA according to claim 3, wherein in the step S5, the internal valid index-CH index is based on the measure of the intra-class dispersion matrix and the inter-class dispersion matrix of all samples, as follows:
wherein the content of the first and second substances,kto be the number of clusters to be clustered,nfor the purpose of the total number of samples,tr(B(k))is the trace of the inter-class separation difference matrix,tr(W(k))is the trace of the intra-class dispersion matrix; the CH index corresponding to the maximum valuekNamely the optimal clustering number.
5. The coherent group clustering method combining the S-transform and the 2DPCA according to claim 1, further comprising a step S6: the synchrophasor measurement device acquires the power angle information increment of the M generators after the fault once again every certain time, refreshes the power angle information database, and repeats the steps S2, S3, and S4.
6. A coherent unit clustering system combining S-transform and 2DPCA, comprising:
the synchronous phasor measurement unit is used for acquiring power angle information increment after the M generators are in fault through the synchronous phasor measurement device;
the S transformation unit is used for carrying out rapid S transformation on the power angle swing curve information of each generator acquired by the synchronous phasor measurement unit and acquiring an m multiplied by n dimensional time-frequency domain information matrix;
a dimension reduction unit based on a 2DPCA algorithm, which is used for reducing the dimension of the m multiplied by n dimension time-frequency domain information matrix through the 2DPCA algorithm to obtain a q multiplied by p dimension characteristic index matrix Ts, wherein s is the generator label;
a clustering unit, which clusters the characteristic index matrix Ts by using a self-organizing neural network to obtain a coherent unit;
the clustering by adopting the self-organizing neural network is realized according to the following steps:
step S41: the number of input neurons is the number of elements of the characteristic index matrix Ts, and the input neurons are used as a sample characteristic index matrix; the output layer neuron is initially one; will connect the weightW ij Initially, an arbitrary sample characteristic index matrix is obtained,ifor the input layer neuron labels,jis the output layer neuron label;
step S42: the partition accuracy threshold, which is used to determine the accuracy of the clustering, isλ;
Step S43: inputting a new sample X, calculating Euclidean distance between the new sample and a corresponding clustering center of each output neuron, and recording the minimum distance asD j ,jFor the output layer neuron label corresponding to the minimum distance, the formula is as follows:
in the formula:X i is an input vector; connection weightW ij Is as followsjA cluster center of class group patterns;
step S44: if it isD j <λ, thenjAn output neuron wins, the sample is includedjIn each neuron, updating the neuron connection weight value, wherein the formula is as follows:
wherein the content of the first and second substances,η(n)is the reciprocal of the total number of samples contained in the neuron;
if it isD j >λIf the sample does not belong to any existing output neuron, adding an output neuron, and taking the initial connection weight of the output neuron and the input layer as a sample characteristic index matrix;
step S45: and sequentially inputting all samples according to the steps, wherein the number of final output layer neurons is the number of clustering clusters, each output layer neuron comprises a sample which is a same cluster sample, and the connection weight of the sample and the input layer is the clustering center of the cluster.
7. The coherent team clustering system combining S-transform and 2DPCA according to claim 6, further comprising: and an optimal clustering number calculating unit for determining the optimal clustering number according to an internal effective index-CH index.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710821583.8A CN107590604B (en) | 2017-09-13 | 2017-09-13 | Coherent unit grouping method and system combining S transformation and 2DPCA |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710821583.8A CN107590604B (en) | 2017-09-13 | 2017-09-13 | Coherent unit grouping method and system combining S transformation and 2DPCA |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107590604A CN107590604A (en) | 2018-01-16 |
CN107590604B true CN107590604B (en) | 2020-09-04 |
Family
ID=61051260
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710821583.8A Active CN107590604B (en) | 2017-09-13 | 2017-09-13 | Coherent unit grouping method and system combining S transformation and 2DPCA |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107590604B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110110673B (en) * | 2019-05-10 | 2020-11-27 | 杭州电子科技大学 | Face recognition method based on bidirectional 2DPCA and cascade forward neural network |
CN110348540B (en) * | 2019-07-24 | 2021-06-01 | 国电南瑞科技股份有限公司 | Clustering-based method and device for screening transient power angle stability faults of power system |
CN112531706B (en) * | 2020-12-18 | 2023-04-18 | 东北电力大学 | Coherent cluster identification method based on complex invariance and deep neural network |
CN115754469B (en) * | 2023-01-05 | 2023-05-09 | 湖南红船科技有限公司 | Unmanned aerial vehicle weak signal detection and extraction method, system, equipment, medium and terminal |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105071380A (en) * | 2015-07-28 | 2015-11-18 | 天津大学 | Power system wide-area coherency identification method and device |
CN105429138B (en) * | 2015-12-11 | 2019-02-05 | 清华大学 | Coherent Generator Group recognition methods and device based on wavelet transformation |
CN106487026B (en) * | 2016-08-25 | 2019-04-23 | 天津大学 | A kind of Coherent Generator Group method for quickly identifying based on wide area measurement information |
-
2017
- 2017-09-13 CN CN201710821583.8A patent/CN107590604B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN107590604A (en) | 2018-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cai et al. | Classification of power quality disturbances using Wigner-Ville distribution and deep convolutional neural networks | |
CN107590604B (en) | Coherent unit grouping method and system combining S transformation and 2DPCA | |
CN108241873B (en) | A kind of intelligent failure diagnosis method towards pumping plant main equipment | |
He et al. | A novel data-driven situation awareness approach for future grids—Using large random matrices for big data modeling | |
Liu et al. | High-precision identification of power quality disturbances under strong noise environment based on FastICA and random forest | |
CN111738309B (en) | Multi-scale analysis and integrated learning gas sensor fault mode identification method | |
CN103245907A (en) | Artificial circuit fault diagnosis pattern sorting algorithm | |
CN104951787B (en) | The electrical energy power quality disturbance recognition methods of dictionary learning is differentiated under a kind of SRC frame | |
CN107895202B (en) | Short-term wind speed prediction method based on multi-view wind speed mode mining | |
CN106503279B (en) | A kind of modeling method for transient stability evaluation in power system | |
CN102364490B (en) | Automatic synchronization recognition method based on hierarchical analyzing model | |
Luo et al. | Graph convolutional network-based interpretable machine learning scheme in smart grids | |
CN105071380A (en) | Power system wide-area coherency identification method and device | |
Ma et al. | Underdetermined blind source separation based on source number estimation and improved sparse component analysis | |
CN113256443B (en) | Nuclear power water pump guide bearing fault detection method, system, equipment and readable storage medium | |
CN112882899B (en) | Log abnormality detection method and device | |
CN104834816A (en) | Short-term wind speed prediction method | |
Zhou et al. | Bearing fault diagnosis based on RF-PCA-LSTM model | |
CN105116323B (en) | A kind of electrical fault detection method based on RBF | |
Wang et al. | Fuzzy c-means clustering for power system coherency | |
CN110943473A (en) | Generator coherence identification method based on wide area measurement system and clustering theory | |
Wang et al. | Cluster division in wind farm through ensemble modelling | |
Dabou et al. | Supervised learning of overcomplete dictionaries for rapid response-based dynamic stability prediction | |
Zhu et al. | Aiming to Complex Power Quality Disturbances: A Novel Decomposition and Detection Framework | |
CN115713032A (en) | Power grid prevention control method, device, equipment and medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |