Abstract
With the large-scale development of renewable energy power, China has faced with the challenges of the reverse regional distribution of wind and solar resources and power load, as well as the intermittency and randomness of renewable energy power. Therefore, China is vigorously developing ultra-high voltage direct current (UHVDC) transmission technology to solve the problem of insufficient flexibility caused by the uncertainty of renewable energy and the fluctuation of multi-energy loads in integrated community energy systems. UHVDC plays an increasingly pivotal role in the west-east transmission system in China’s power system due to its high transmission capacity and long transmission distance. Once the fault occurs in the ultra-high voltage direct (UHVD) transmission line, quick and accurate fault location identification is of great significance. Hence, this paper proposes a UHVDC transmission line diagnosis method based on wavelet analysis for integrated community energy systems. Wavelet transform (WT) is used to decompose the transient signal on a multi-scale, and then power systems computer-aided design (PSCAD) software is utilized for simulation calculation to obtain the singular spectrum entropy of each layer and facilitate wavelet transformations for signal denoising with advanced tools such as MATLAB. The prediction results can distinguish outside the rectification side fault, within the rectification side fault, and outside the inverter fault with an accuracy of 100%. A large number of simulations demonstrate that combining singular spectrum entropy with support vector machines (SVM) has emerged as a robust technique for integrated community energy systems, suggesting its potential as a standard method in UHVDC transmission line diagnosis. This study is of significant reference for realizing the complementarity of multiple types of power supply and ensuring a reliable power supply.
1 Introduction
As global energy demand rises and the public’s awareness of environmental protection increases, there is a growing demand for clean energy sources. China’s energy resource centers and load centers are distributed in the reverse direction, and the energy flow is in the pattern of “sending coal from the west to the east and from the north to the south” as well as “sending electricity from the west to the east and electricity from the north to the south.” China’s wind and solar energy resources have great potential for development, with the most abundant distribution in the northwest, Xinjiang, and Inner Mongolia regions. Although these regions have abundant wind and solar resources, the load power is low due to factors such as productivity and population, which is exactly the opposite of the central and eastern regions, leading to serious wind and solar abandonment in the northwest region (Zhou, 2019). In the future, China’s energy production centers are expected to move west and north, while the gravity center of demand may remain in the central and eastern regions for a long time, and the scale and distance of energy flow will further increase. For future large-scale access to new energy under the conditions of a large-capacity “west to east” scenario, there is a need to plan ahead with the appropriate transmission mode (; ), posing higher requirements for electric energy transmission technology. Among various energy transmission methods, UHVDC technology has played an increasingly important role in power transmission and has well solved the problem of unbalanced distribution between energy centers and load centers due to its advantages of large transmission power, low line cost, and good control performance (Wang et al., 2007; ). From an economic point of view, UHVDC is more suitable for long-distance power supply; the longer the distance, the higher the economy. From the existing UHV lines in our country, UHVDC transmission is generally chosen for ultra long-distance power transmission. The main way to fully implement the new national energy development plan is to transmit wind and photovoltaic power from large energy bases in the “Three Northern Regions” to the eastern and central load centers through the UHVDC transmission system.
Numerous studies have been dedicated to efficiently accommodating excess renewable generation, reducing renewable curtailment, and improving overall energy efficiency. According to the development research of the State Grid, clean energy, such as wind power located in the energy base of North China (including western Inner Mongolia), can be transmitted through UHVDC systems. Similarly, clean energy from energy bases in the northeast and northwest (including Xinjiang) regions can only be transported through UHVDC systems (Zhen et al., 2008; ). and Wang L. et al. (2018) established an planning model to minimize investment, construction, and operating costs. Zeng et al. (2018) studied methods for handling multiple uncertainties in the programming models. Zhou et al. (2020) studied the planning method based on maximum utilization hours and established to minimize construction and operating costs throughout the entire cycle. However, the above studies have overlooked the complex geographical environment along the UHVDC transmission line. The fault rate of the transmission line is high in the event of extreme events, making it difficult to ensure continuous energy supply to important users and rapid recovery from faults, which compromises the safety and reliability of the direct current transmission systems (). Therefore, the protection of UHV transmission lines becomes paramount, highlighting the importance of fast and accurate fault diagnosis for the safe and stable operation of the UHVDC transmission systems. Based on this phenomenon, , , Wang Y. et al. (2018), and conducted certain research. At present, transient protection utilizing the boundary to the attenuation characteristics of high-frequency quantities is the development direction of UHVDC transmission line protection. Studies by Yong et al. (2009), Yan et al. (2017), , and have made significant progress in this realm. However, these studies overlook the complex regional environment along UHVDC transmission lines. In the case of extreme events, it is difficult to ensure the continuous energy supply to important users and the rapid recovery of faults due to the high transmission line failure rate seriously affecting the safety and reliability of the HVDC transmission system.
Support Vector Machines (SVM) is a machine learning method based on statistical learning developed in the 1990s (), aimed at classifying samples by finding an optimal classification hyperplane between them (Xue et al., 2015). Given that SVM is a superior statistical learning method known for recognizing patterns in small sample data, the objective of distinguishing fault types is achieved through training and testing the sample sets (; ; ; Zhen et al., 2019). Zhu et al. (2011) introduced a new identification method combining SVM and multi-resolution singular spectral entropy to address the problem of classifying islanding and grid interference. Considering the practical difficulty in obtaining a large number of typical samples of UHVDC line boundary and fault transient signals, the combination of multi-resolution singular spectrum entropy and support vector machine is applied to the problem of fault identification inside and outside the UHVDC line transient protection. This approach demonstrates the effectiveness in classifying fault categories with minimal sample data, yielding the desired outcomes (; ). Hence, this study mainly proposes a transient protection method for UHVDC lines based on SVM and multi-resolution singular spectrum entropy. Section 2 provides a brief overview of the methodology designed in the study. Section 3 discusses the modeling of verification. Section 4 elaborates on the analysis of frequency domain waveform and singular spectral entropy to classify signals in three different positions. Finally, Section 5 concludes key insights derived from this study.
2 Methodology
2.1 Singular spectral entropy
Singular spectrum analysis represents an effective time-domain analysis method that transforms the embedding space into an equivalent orthogonal coordinate system. This transformation facilitates obtaining signal trajectories in the subspace with the minimum embedding dimension, thereby eliminating linear dependencies and artificial symmetry between delay coordinates. As a result, this process enhances the signal-to-noise ratio and sharpens signal singularity. However, as a time-domain analysis method, traditional singular spectrum analysis is not conducive to multi-scale monitoring of signal singularity and fault localization. Thus, the multi-resolution singular spectrum entropy combines the idea of multi-resolution analysis and information entropy to characterize the singular state of signal energy distribution and the probability of geometric feature distribution at different scales (Wu et al., 2012).
2.1.1 Wavelet singular spectral entropy extraction
1) A discrete sampling sequence (N is the sampling point) is provided. If represents the approximate value at scale , denoted as , then the formula for discrete dyadic wavelet transform can be expressed in Eq. 1.
where
and
are low-pass filters and high-pass filters respectively.
In addition,
and
represent the approximate and detailed parts of the signal at scale, respectively. The discrete signal
is decomposed into
by scale
which contains information from different frequency bands from the high-frequency to the low-frequency.
2) Perform wavelet reconstruction on each decomposed layer of the signal. Suppose only the wavelet coefficients of a certain frequency band are retained and the wavelet coefficients of other frequency bands are set to 0. In that case, a singular spectrum analysis of the reconstructed time-domain signal is performed. The reconstruction formula is in Eq. 2.
where
and
are respectively the dual operators of
and
.
3) Set on layer , the discrete wavelet reconstructed signal from the multi-resolution analysis is , and the reconstructed signal is mapped to the embedding space. For the sequence of point , the reconstruction state at discrete time is: , where represents the reconstruction delay, denotes the embedding dimension, and the corresponding reconstruction trajectory is in Eq. 3.
Among them, the element of
is the correlation function, in Eq.
4.
4) Perform singular value decomposition on each layer of matrix , assuming the singular value is . Then, forms the singular spectral value of the signal on that layer. Suppose is the number of non-zero singular values. In that case, the value of reflects the number of different patterns in each column of the feature matrix . The magnitude of singular value reflects the proportion of corresponding patterns in the total pattern. Therefore, the singular value is a partition of vibration signals in the time-frequency domain of wavelet signals. According to information entropy theory, the singular spectral entropy of the reconstructed signal at this level can be calculated in Eq. 5.
where
.
2.1.2 Wavelet singular spectrum entropy eigenvector extraction
Specific steps of feature vector extraction:
1) Perform WT on the selected voltage signal for analysis, use phase space reconstruction, equivalently exchange the embedding space into an orthogonal coordinate system, and construct a matrix in Eq. 6.
2) Perform singular value decomposition of matrix for each layer.
3) Perform singular spectral entropy operation on the singular values of each layer, combine the singular spectral entropy values of each layer, and obtain the eigenvector in Eq. 7.
where
represent the entropy value of each layer, and the corresponding frequency band is from high frequency to low frequency.
2.2 Support vector machines (SVM)
The singular spectral entropy value of the transient signal at the initial moment of a fault in a DC transmission line can reflect the information regarding the fault in the UHVDC transmission line. However, the values are all less than 1. Therefore, due to the presence of abnormal data and a limited number of faults, along with a small number of training samples, it is proposed to use the singular spectral entropy as the input for the SVM classification model to determine the type of fault.
The learning strategy of SVM is to find the optimal hyperplane, maximizing the interval between hyperplanes and transforming it into a convex quadratic problem. Nonlinear mapping is used to map the sample space to a high-dimensional feature space, and linear learning machines are applied in the feature space to solve nonlinear classification and regression problems in the sample space (Zhang et al., 2020). Figure 1 shows the network structure of the SVM, with the detailed definitions of variables and functions provided in the following text.
FIGURE 1
2.2.1 Principles of SVM
For a non-linear separable sample set, introducing the relaxation variable ( ≥0) and penalty factor , the objective function is in Eq. 8:where denotes the normal direction vector for dividing the hyperplane. The decision function of the optimal hyperplane is in Eq. 9.where indicates that and are the parameters for determining the optimal partition hyperplane. and represent points in the sample space, that is, support vectors. Moreover, symbolizes the value of the prediction factor, and indicates a signed function.
2.2.2 Selection of kernel function and parameter optimization methods
This paper uses the Gauss kernel function, whose formula is in Eq. 10:where represents the distance between two vectors and denotes the width of the kernel function.
The bilinear search method and grid search method are methods used to determine the SVM classifier. Unlike the bilinear search method, the grid search method offers the advantage of parallel SVM training for and as they are independent of each other. Therefore, this paper adopts the grid search method, taking M and N separately to form combination , which are trained separately. Subsequently, the learning accuracy is estimated, obtaining the combination with the highest learning accuracy among all combinations as the optimal parameter.
2.3 Failure diagnosis method
, Wang et al. (2015), and analyzed the attenuation characteristics of UHVDC transmission lines and their boundaries for high-frequency quantities of fault transient signals. According to the attenuation characteristics of the fault transient signal caused by the lines and their boundaries, wavelet decomposition applied to the fault signals in three different positions, including outside the rectification side, within the rectification side, and outside the inverter side to obtain different reconstructed high-frequency coefficient amplitudes (; ). Subsequently, phase space reconstruction is performed on the high-frequency coefficients of each layer, followed by singular value decomposition on the phase space of each layer to obtain singular values of different sizes (; ; Xiang et al., 2018).
Since the magnitude of singular values reflects the differences between various fault states, singular spectral entropy can quantitatively describe the degree of change (; ; ). Singular spectral entropy serves as a reflection of the uncertainty of signal energy. The simpler the signal component, the smaller the singular spectral entropy value, indicating more concentrated signal energy. Conversely, the more complex the signal components, the larger the singular spectral entropy value, indicating that the energy is more dispersed, and the signal is more evenly distributed throughout the entire frequency component (; Wang and Zheng, 2014).
Therefore, the singular spectral entropy of the fault signal at different decomposition levels is calculated, and then the fault signal is diagnosed through the singular spectral entropy (Xing et al., 2016). Singular spectral entropy can be used to reflect the different changes in fault signals at three different locations, including outside the rectification side, within the rectification side, and outside the inverter side after boundary and line attenuation.
3 Simulation verification
The simulation model is established based on the actual parameters of the Yunguang ±800 kV UHVDC transmission system. For training, 3 external grounding faults are considered on the rectifier side, 3 external grounding faults on the inverter side, and 9 internal grounding faults. The training set is shown in Table 1. For testing, 3 external grounding faults on the rectifier side, 3 external grounding faults on the inverter side, and 9 internal grounding faults are used. The testing set is shown in Table 2, and the identification results of the grounding fault are shown in Table 3.
TABLE 1
| Number | e1 | e2 | e3 | e4 | e5 | e6 | y |
|---|---|---|---|---|---|---|---|
| 1 | 4.234 851 | 4.007 698 | 3.530 038 | 2.426 193 | 2.153 935 | 1.620 465 | 1 |
| 2 | 4.230 357 | 4.007 169 | 3.541 980 | 2.470 981 | 2.151 344 | 1.416 290 | 1 |
| 3 | 4.226 877 | 4.006 925 | 3.542 546 | 2.491 254 | 2.498 893 | 1.239 845 | 1 |
| 4 | 3.992 142 | 3.968 133 | 3.608 143 | 2.954 392 | 2.368 172 | 1.700 132 | 2 |
| 5 | 4.001 413 | 3.968 919 | 3.610 928 | 2.967 437 | 2.371 527 | 1.709 861 | 2 |
| 6 | 4.009 175 | 3.964 558 | 3.609 616 | 2.903 287 | 2.325 60 | 1.731 467 | 2 |
| 7 | 4.032 883 | 3.958 139 | 3.604 873 | 2.964 006 | 2.355 154 | 1.554 912 | 2 |
| 8 | 3.966 016 | 3.901 362 | 3.610 336 | 2.880 234 | 2.321 510 | 1.699 400 | 2 |
| 9 | 3.977 207 | 3.900 431 | 3.609 016 | 2.895 399 | 2.342 394 | 1.715 578 | 2 |
| 10 | 3.988 813 | 3.897 810 | 3.610 564 | 2.886 855 | 2.349 761 | 1.766 450 | 2 |
| 11 | 4.022 361 | 3.904 598 | 3.609 334 | 2.880 222 | 2.334 955 | 1.632 571 | 2 |
| 12 | 4.023 365 | 3.912 566 | 3.609 254 | 2.881 215 | 2.335 145 | 1.653 458 | 2 |
| 13 | 4.260 101 | 4.006 756 | 3.550 898 | 2.533 301 | 2.185 532 | 1.645 652 | 3 |
| 14 | 4.254 768 | 3.998 255 | 3.586 552 | 2.478 692 | 2.184 178 | 1.564 385 | 3 |
| 15 | 4.253 896 | 3.989 569 | 3.612 551 | 2.434 569 | 2.183 158 | 1.512 548 | 3 |
Training set of line short circuit.
TABLE 2
| Number | e1 | e2 | e3 | e4 | e5 | e6 | Actual classification y | Output results y |
|---|---|---|---|---|---|---|---|---|
| 1 | 4.233 234 | 4.011 294 | 3.545 682 | 2.444 781 | 2.175 692 | 1.731 926 | 1 | 1 |
| 2 | 4.237 742 | 4.010 300 | 3.533 309 | 2.420 985 | 2.160 731 | 1.695 838 | 1 | 1 |
| 3 | 4.239 541 | 4.009 894 | 3.527 514 | 2.400 612 | 2.155 879 | 1.667 456 | 1 | 1 |
| 4 | 4.127 146 | 3.950 120 | 3.587 153 | 2.939 680 | 2.372 009 | 1.761 636 | 2 | 2 |
| 5 | 4.124 025 | 3.955 669 | 3.590 703 | 2.945 490 | 2.363 065 | 1.710 628 | 2 | 2 |
| 6 | 4.124 349 | 3.965 404 | 3.594 802 | 2.874 286 | 2.328 265 | 1.644 996 | 2 | 2 |
| 7 | 4.162 143 | 3.973 289 | 3.606 325 | 2.872 658 | 2.303 787 | 1.570 906 | 2 | 2 |
| 8 | 4.073 362 | 3.655 916 | 3.584 639 | 2.880 284 | 2.339 811 | 1.800 524 | 2 | 2 |
| 9 | 4.083 346 | 3.671 308 | 3.591 093 | 2.826 990 | 2.307 579 | 1.781 310 | 2 | 2 |
| 10 | 4.074 555 | 3.685 678 | 3.595 680 | 2.729 995 | 2.274 955 | 1.766 197 | 2 | 2 |
| 11 | 4.131 489 | 3.746 499 | 3.611 569 | 2.720 045 | 2.257 295 | 1.647 411 | 2 | 2 |
| 12 | 4.121 654 | 3.745 412 | 3.612 125 | 2.720 001 | 2.256 748 | 1.667 845 | 2 | 2 |
| 13 | 4.256 340 | 3.984 575 | 3.615 606 | 2.705 496 | 2.239 537 | 1.682 504 | 3 | 3 |
| 14 | 4.264 150 | 3.994 711 | 3.607 863 | 2.674 488 | 2.227 026 | 1.671 522 | 3 | 3 |
| 15 | 4.268 754 | 3.998 785 | 3.600 988 | 2.645 456 | 2.225 625 | 1.658 958 | 3 | 3 |
Testing set of line short circuit.
TABLE 3
| C | δ | Number of training samples | Number of testing samples | Testing accuracy/% | ||
|---|---|---|---|---|---|---|
| External grounding fault on the rectifier side | Internal grounding fault | External grounding fault on the inverter side | ||||
| 2 | 2 | 15 | 15 | 100 | 100 | 100 |
Identification results.
As can be seen in Table 6, electing appropriate C and can achieve the best testing accuracy.
When = 2 and = 2, the testing accuracy of the external grounding fault on the rectifier side, the internal grounding fault, and the external grounding fault on the inverter side is 100%.
4 Results and discussions
4.1 Frequency domain waveform analysis
This paper uses PSCAD to establish a model, with a sampling frequency of 40 kHz, a sampling time of 0.05 s, and a total of 2000 sampling points. The WT uses the db4 wavelet with a decomposition level of 6.
The reconstructed waveform of the high-frequency coefficients obtained from the WT reflecting ground faults at three different locations, outside the rectification side, within the rectification side, and outside the inverter side, is shown in Figures 2–4.
FIGURE 2
FIGURE 3
FIGURE 4
It can be seen from the waveform that the waveform changes of fault information at different positions are noticeable in the frequency domain, and the amplitude in the high-frequency range also varies.
4.2 Singular spectral entropy analysis
Three different fault signals, namely, external ground fault on the rectification side, internal ground fault, and external ground fault on the inverter side, are decomposed at three scales, and the wavelet singular spectral entropy is computed for each layer to form a feature vector. Three sets of characteristic vectors are taken for the external ground fault signals on the rectifier side, internal ground fault signals, and external ground fault signals on the inverter side, as listed in Table 4.
TABLE 4
| Number | Fault location | e1 | e2 | e3 | y |
|---|---|---|---|---|---|
| 1 | External on the rectifier side | 4.192 7 | 3.966 0 | 3.527 6 | 1 |
| 2 | External on the rectifier side | 4.192 5 | 3.958 6 | 3.501 9 | 1 |
| 3 | External on the rectifier side | 4.192 3 | 3.951 7 | 3.4886 | 1 |
| 4 | Internal | 3.959 7 | 3.901 1 | 3.609 3 | 2 |
| 5 | Internal | 3.974 4 | 3.896 8 | 3.607 8 | 2 |
| 6 | Internal | 3.999 2 | 3.889 2 | 3.60 63 | 2 |
| 7 | External on the inverter side | 4.248 8 | 3.983 4 | 3.606 9 | 3 |
| 8 | External on the inverter side | 4.253 8 | 4.011 8 | 3.572 0 | 3 |
| 9 | External on the inverter side | 4.255 7 | 4.035 4 | 3.542 5 | 3 |
Feature vectors of three different locations.
The in Table 4 denotes the category, and the meaning of in the following tables is the same as that in Table 4.
According to the presented data in Table 4, in the low frequency band , there is not much difference in the singular entropy values of the external grounding fault on the rectifier side, the internal grounding fault, and the external grounding fault on the inverter side. This phenomenon arises due to the attenuation of the low-frequency signal along the line and line boundary can be approximated to be zero. Consequently, the distribution of energy in this frequency band is relatively concentrated, leading to relatively low uncertainty in energy distribution. In the high-frequency range, the singular spectral entropy values of and are relatively large, exhibiting a significant difference, with values of the internal < the external on the rectifier side < the external on the inverter side. This is because when there is an external grounding fault on the inverter side, the high-frequency signal needs to pass through the double attenuation effect of the line and the line boundary to reach the installation site of the rectifier side protection. This severe high-frequency signal attenuation and the relatively uniform energy distribution in this frequency band increase the uncertainty of energy distribution.
In addition, Table 4 demonstrates the relative stability of wavelet singular spectral entropy for fault information at the same location. Therefore, this study uses the singular spectral entropy of faults outside the rectification side, within the rectification side, and outside the inverter side as input variables for SVM to classify faults.
4.3 SVM fault recognition algorithm
SVM algorithm process:
1) Perform wavelet multi-resolution decomposition on the fault signals at three different positions to obtain the corresponding reconstructed high-frequency coefficients;
2) Calculate the singular spectral entropy of the wavelet reconstruction coefficients of the fault signals at three different positions to obtain the feature vectors;
3) Utilize a portion of the feature vectors as the training set and apply the grid search method to determine the SVM classifier parameter ;
4) Input the remaining part of the feature vector as the test set into the SVM classifier for testing and obtaining the classification results.
The SVM algorithm process is shown in Figure 5.
FIGURE 5
The nine samples in Table 4 are employed as a training set, and the parameters of the SVM classifier are obtained using the grid search method.
Then, three groups of feature vectors are selected to form a testing set and input to the SVM for testing. The test set samples are shown in Table 5.
TABLE 5
| Number | Fault location | e1 | e2 | e3 | y |
|---|---|---|---|---|---|
| 1 | External on the rectifier side | 4.181 4 | 3.972 4 | 3.540 4 | 1 |
| 2 | External on the rectifier side | 4.182 2 | 3.976 8 | 3.536 2 | 1 |
| 3 | External on the rectifier side | 4.183 0 | 3.978 2 | 3.521 2 | 1 |
| 4 | Internal | 4.000 1 | 3.897 8 | 3.608 2 | 2 |
| 5 | Internal | 4.005 1 | 3.900 8 | 3.610 3 | 2 |
| 6 | Internal | 4.009 8 | 3.901 7 | 3.612 9 | 2 |
| 7 | External on the inverter side | 4.258 5 | 4.031 0 | 3.551 5 | 3 |
| 8 | External on the inverter side | 4.255 8 | 4.010 0 | 3.538 5 | 3 |
| 9 | External on the inverter side | 4.251 3 | 3.998 9 | 3.525 9 | 3 |
Testing set of three different locations.
The testing set is tested in
Table 5, and the following results are obtained:
1) When = 1024, = 1024, the correct rate is very low, which is 13.33%. Only the external grounding fault on the inverter side is correctly classified, and the faults in the external grounding on the rectifier side and the internal grounding are not correctly classified.
2) When = 512, = 512, the correct rate is 33.33%, and there are two samples with the external grounding fault on the rectifier side that is not correctly classified.
3) When = 128, = 128, the correct rate is 83.33%, and there is one sample with the external grounding fault on the rectifier side that is not correctly classified.
4) When = 2, = 2, the correct rate is 100%, and all samples are correctly classified. The classification results are shown in Table 6.
TABLE 6
| Number | e1 | e2 | e3 | Actual classification y | Output results y |
|---|---|---|---|---|---|
| 1 | 4.181 4 | 3.972 4 | 3.540 4 | 1 | 1 |
| 2 | 4.182 2 | 3.976 8 | 3.536 2 | 1 | 1 |
| 3 | 4.183 0 | 3.978 2 | 3.521 2 | 1 | 1 |
| 4 | 4.000 1 | 3.897 8 | 3.608 2 | 2 | 2 |
| 5 | 4.005 1 | 3.900 8 | 3.610 3 | 2 | 2 |
| 6 | 4.009 8 | 3.901 7 | 3.612 9 | 2 | 2 |
| 7 | 4.258 5 | 4.031 0 | 3.551 5 | 3 | 3 |
| 8 | 4.255 8 | 4.010 0 | 3.538 5 | 3 | 3 |
| 9 | 4.251 3 | 3.998 9 | 3.525 9 | 3 | 3 |
SVM diagnosis of three different locations.
The data in Table 6 indicate that an ideal accuracy can be achieved for small sample training and learning by selecting appropriate kernel parameters.
When = 2, = 2, the above external grounding fault on the rectifier side, the internal grounding fault, and the external grounding fault on the inverter side are classified by SVM, resulting in the classification diagram is obtained, as shown in Figure 6.
FIGURE 6
As can be seen from the SVM diagnosis results in Table 6 and the SVM classification diagram in Figure 6, the trained SVM classifier parameters can correctly classify the grounding fault signals in three different positions.
5 Conclusion
This study presented a fault diagnosis method for UHVDC transmission lines in a regional integrated energy system, establishing an accurate mathematical model and making a multi-scale analysis of fault signals using wavelet analysis. Then, the method for fault identification inside and outside the UHVDC line area using a multi-resolution singular spectrum entropy and SVM was proposed by analyzing the UHVDC transmission line boundary and the attenuation effect of lines on the high-frequency of fault transient signals. For this purpose, MATLAB software was utilized to facilitate wavelet transformations for signal denoising and various machine learning techniques for predictive modeling. Compared with the traditional artificial neural network, the revealed mathematical model offers the advantages of fewer training samples, shorter training time, less overfitting, and higher recognition accuracy.
This study revealed the fault diagnosis method for UHVDC transmission lines within a regional integrated energy system. The fault identification problem of transient protection for UHVDC transmission lines inside and outside the region was transformed into a pattern classification problem. As a result, the sample identification accuracy improved from 13.3% to 100% by establishing a mathematical model and selecting a suitable SVM classifier parameter. Thus, the external fault on the rectifier side, the internal fault, and the external fault on the inverter side could be distinguished at the same time, and automation in the whole classification process could be realized. Given the challenges of acquiring a large number of typical samples of UHVDC line boundary and fault transient signals of the line, the combination of multi-resolution singular spectrum entropy and SVM was applied to the problem of fault identification inside and outside the UHVDC line transient protection. With only small sample data, the fault categories could be effectively classified, and the expected results could be achieved. The accuracy reached 100%, indicating the correctness and effectiveness of the method. This study not only makes a great contribution to the fault diagnosis of UHVDC transmission lines in regional integrated energy systems but also lays a new theoretical foundation for the fault diagnosis of other regional integrated energy systems in the process of energy signal transmission.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions
TY: Conceptualization, Investigation, Methodology, Resources, Software, Validation, Writing–original draft, Writing–review and editing. JL: Conceptualization, Data curation, Methodology, Writing–original draft, Writing–review and editing. XZ: Investigation, Methodology, Writing–original draft, Writing–review and editing. KL: Investigation, Writing–original draft, Writing–review and editing. LF: Software, Writing–original draft, Writing–review and editing. ZD: Validation, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The authors gratefully acknowledge funding from the Foundation of Liaoning Provincial Key Laboratory of Energy Storage and Utilization (Grant No.CNNK202318) and the Natural Science Foundation of Liaoning Province of China (Grant No. 2023-MSLH-314).
Acknowledgments
The authors would like to acknowledge State Grid Yingkou Electric Power Supply Company for providing experimental data.
Conflict of interest
Author XZ was employed by the company State Grid Yingkou Electric Power Supply Company.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AbdullahA. (2017). Ultrafast transmission line fault detection using a dwt based ann. IEEE Trans. Industry Appl.54 (99), 1182–1193. 10.1109/TIA.2017.2774202
2
AdlyA. R.AleemS. H. E. A.AlgabalawyM. A.JuradoF.AliZ. M. (2020). A novel protection scheme for multi-terminal transmission lines based on wavelet transform. Electr. Power Syst. Res.183, 106286. 10.1016/j.epsr.2020.106286
3
BauerS.NolteL. P.ReyesM. (2011). Fully automatic segmentation of brain tumor images using support vector machine classification in combination with hierarchical conditional random field regularization. Springer Berlin Heidelberg. 10.1007/978-3-642-23626-6_44
4
ChenX.LiH.LiangY.WangG. (2020). A protection scheme for hybrid multi-terminal hvdc networks utilizing a time-domain transient voltage based on fault-blocking converters. Int. J. Electr. Power and Energy Syst.118, 105825. 10.1016/j.ijepes.2020.105825
5
DingJ.FangX. S.SongY. T.ZhengC.TangX. J.YaoW. (2021). Conception of electricity and hydrogen integrated energy network for renewableenergy transmission in western China under background of carbon neutralization. Automation Electr. Power Syst.45 (24), 1–9. 10.7500/AEPS20210211002
6
LiC.GoleA. M.ZhaoC. (2018a). A fast dc fault detection method using dc reactor voltages in hvdc grids. IEEE Trans. Power Deliv.33, 2254–2264. 10.1109/TPWRD.2018.2825779
7
LiS. L.ChenW.YinX.ChenD.TengY. (2019). A novel integrated protection for VSC-HVDC transmission line based on current limiting reactor power. IEEE Trans. Power Deliv.35 (1), 226–233. 10.1109/TPWRD.2019.2945412
8
LiY.WuL.LiJ.XiongL.ZhangX.SongG.et al (2018b). Dc fault detection in meshed mtdc systems based on transient average value of current. IEEE Trans. Power Deliv.67, 1932–1943. 10.1109/TIE.2019.2907499
9
LiuF.ChenH. (2017). Fault diagnosis of commutation failure based on wavelet packet decomposition and generalized regression neural network. Automation Instrum.6, 22–25. 10.14016/j.cnki.1001-9227.2017.06.022
10
MaX. Y.GuoX. B.LeiJ. Y. (2018). Capacity planning method of distributed PV and P2G in muti-energy coupled system. Automation Electr. Power Syst.42, 55–63. 10.7500/AEPS20170602002
11
MeghwaniA.SrivastavaS. C.ChakrabartiS. (2017). A non-unit protection scheme for dc microgrid based on local measurements. IEEE Trans. Power Deliv.32 (1), 172–181. 10.1109/TPWRD.2016.2555844
12
MetidjiB.TaibN.BaghliL.RekiouaT.BachaS. (2013). Phase current reconstruction using a single current sensor of three-phase AC motors fed by SVM-controlled direct matrix converters. IEEE Trans. Industrial Electron.60 (12), 5497–5505. 10.1109/TIE.2012.2232252
13
MohammadiF.RouzbehiK.HajianM.NiayeshK.GharehpetianG. B.SaadH.et al (2021). Hvdc circuit breakers: a comprehensive review. IEEE Trans. Power Electron.36 (36-12), 13726–13739. 10.1109/TPEL.2021.3073895
14
MorenoR.VisairoN.NunezC.RodriguezE. (2014). A novel algorithm for voltage transient detection and isolation for power quality monitoring. Electr. Power Syst. Res.114 (SEP), 110–117. 10.1016/j.epsr.2014.04.009
15
MuzzammelR. (2020). Restricted Boltzmann machines based fault estimation in multi terminal HVDC transmission system, 772, 790. 10.1007/978-981-15-5232-8_66
16
NiazyI.SadehJ. (2013). A new single ended fault location algorithm for combined transmission line considering fault clearing transients without using line parameters. Int. J. Electr. Power and Energy Syst.44 (1), 816–823. 10.1016/j.ijepes.2012.08.007
17
PanY.YinX. G.HuJ. B.HeJ. J.et al (2016). Centralized exploitation and large-scale delivery of wind and solar energies in west China based on flexible DC grid. Power Syst. Technol.40 (12), 3621–3629. 10.13335/j.1000-3673.pst.2016.12.001
18
ShengL.ChenM.ChenD. (2019). A novel hvdc double-terminal non-synchronous fault location method based on convolutional neural network. IEEE Trans. Power Deliv.34 (3), 848–857. 10.1109/TPWRD.2019.2901594
19
ShuH.HanY.HuangR.ZhangY. (2020). Fault model and travelling wave matching based single terminal fault location algorithm for t-connection transmission line: a yunnan power grid study. Energies, 13 (6). 10.3390/en13061506
20
ShuH.TianX.DaiY. (2010). The identification of internal and external faults for and#x00B1;800kV UHVDC transmission line based on S-transform of the polarity wave. Image and Signal Process. Int. Congr.7, 3060–3063. 10.1109/CISP.2010.5646181
21
ShuH. C.ZhangB.ZhangG. B.DuanR. M. (2012). Identification of lightning disturbance in uhvdc transmission lines using correlation degree based on short time window data. Adv. Mater. Res.433-440, 3787–3791. 10.4028/www.scientific.net/AMR.433-440.3787
22
Singh BainsT. P.SidhuT. S.XuZ.VolohI.Dadash ZadehM. R. (2017). Impedance-based fault location algorithm for ground faults in series capacitor compensated transmission lines. IEEE Trans. Power Deliv.33, 189–199. 10.1109/TPWRD.2017.2711358
23
SongG. B.CaiX. L.GaoS. P.SuoN. J. L.YangZ. L. (2011). New pilot protection for vsc-hvdc transmission lines using natural frequency characteristic of current. Gaodianya Jishu/High Volt. Eng.37 (8), 1989–1996. 10.1016/B978-0-444-53599-3.10005-8
24
TangG.LuoX.WeiX. (2013). Multi-terminal hvdc and dc-grid technology. Proceedings of the CSEE33 (10), 8–17.
25
VerraxP.BertinatoA.KiefferM.RaisonB. (2020). Transient-based fault identification algorithm using parametric models for meshed hvdc grids. Electr. Power Syst. Res.185, 106387. 10.1016/j.epsr.2020.106387
26
WangC.WangZ.WangJ.ZhaoD. (2018). Svm-based parameter identification for composite zip and electronic load modeling. IEEE Trans. Power Syst.34 (1), 182–193. 10.1109/TPWRS.2018.2865966
27
WangG.LiZ. K.LiH. F. (2007). Transient based protection for ±800kV UHVDC transmission lines. Automation Electr. Power Syst.31 (21), 40–43. 10.1002/jrs.1570
28
WangL.LiX.LiK.HuY.HuangR. (2015). Research of pseudo bipolar lcc-vsc hybrid hvdc system supplying passive network. Power Syst. Prot. Control.
29
WangL.LiuH.DaiL. V.LiuY. W. (2018). Novel method for identifying fault location of mixed lines. Energies11 (6), 1529–1547. 10.3390/en11061529
30
WangX.ZhengY. (2014). Vibration fault diagnosis for wind turbine based on enhanced supervised locally linear embedding. Adv. Mater. Res.1008-1009, 983–987. 10.4028/www.scientific.net/AMR.1008-1009.983
31
WangY.ZhangN.ZhuoZ. Y.KangC.KirschenD. (2018). Mixed-integer linear programming-based optimal configuration planning for energy hub: starting from scratch. Appl. Energy210, 1141–1150. 10.1016/j.apenergy.2017.08.114
32
WuS. D.WuP. H.WuC. W.DingJ. J.WangC. C. (2012). Bearing fault diagnosis based on multiscale permutation entropy and support vector machine. Entropy14 (8), 1343–1356. 10.3390/e14081343
33
XiangW.YangS.XuL.ZhangJ.LinW.WenJ. (2018). A transient voltage based dc fault line protection scheme for mmc based dc grid embedding dc breakers. IEEE Trans. Power Deliv.34, 334–345. 10.1109/tpwrd.2018.2874817
34
XingL.YangD.ZhangD.LiW.ChenQ. (2016). Fault transient simulation analysis for hvdc transmission lines. IEEE. 10.1109/DRPT.2015.7432248
35
XueH.ZhangK.LiB.PengC. (2015). Fault diagnosis of transformer based on the cuckoo search and support vector machine. Dianli Xit. Baohu yu Kongzhi/Power Syst. Prot. Control43 (8), 8–13.
36
YanY.ShengG.QiuR. C.JiangX. (2017). Big data modeling and analysis for power transmission equipment: a novel random matrix theoretical approach. IEEE Access6, 7148–7156. 10.1109/ACCESS.2017.2784841
37
YongJ.Zheng-YouH. E.JingZ. (2009). A method to identify voltage sag sources in distribution network based on wavelet entropy and probability neural network. Power Syst. Technol.33 (16), 63–69. 10.1016/j.apm.2007.10.019
38
ZengX. Q.ZhangJ. W.WangX. L. (2023). Optimal configuration of regional integrated energy system taking into account multiple uncertainties and the participation of concentrating solar power stations. High. Volt. Eng.49 (1), 353–363. 10.13336/j.1003-6520.hve.20211326
39
ZhangX.HanP.XuL.ZhangF.GaoL. (2020). Research on bearing fault diagnosis of wind turbine gearbox based on 1DCNN-PSO-SVM turbine gearbox based on 1dcnn-pso-svm. IEEE Access8, 192248–192258. 10.1109/ACCESS.2020.3032719
40
ZhenY.YangJ.ZhangB. (2019). Study on numerical sampling stability of traveling wave protection based on wavelet transform. Power Syst. Prot. Control47 (9), 42–48. 10.7667/PSPC180633
41
ZhenY. L.YinB. S.WenL. Z.YunZ. Z. (2008). Study on voltage class series for hvdc transmission system. Proc. CSEE. 10.13334/j.0258-8013.pcsee.2008.10.001
42
ZhouC. C.MaX. Y.GuoZ. G.HuY.et al (2020). User-level integrated energy system planning for engineering applications. Trans. China Electrotech. Soc.35 (13), 2843–2854. 10.19595/j.cnki.1000-6753.tces.191056
43
ZhouX. (2019). New generation power system and energy Internet. Electrotech. Appl.38 (1), 4–6.
44
ZhuY.ShiX.LiP. (2011). Classification of islanding and grid disturbance based on multi-resolution singular spectrum entropy and svm. Proc. CSEE31 (7), 64–70. 10.3788/gzxb20114002.0199
Summary
Keywords
integrated community energy system, wavelet analysis, UHVDC power transmission, singular spectrum entropy, support vector machines, characteristic vector
Citation
Yuan T, Liang J, Zhang X, Liang K, Feng L and Dong Z (2024) UHVDC transmission line diagnosis method for integrated community energy system based on wavelet analysis. Front. Energy Res. 12:1401285. doi: 10.3389/fenrg.2024.1401285
Received
15 March 2024
Accepted
11 April 2024
Published
10 May 2024
Volume
12 - 2024
Edited by
Xiaochen Hou, Shandong University of Technology, China
Reviewed by
Fubin Yang, Beijing University of Technology, China
Yonghong Xu, Beijing University of Technology, China
Genyu Xu, Yunnan University, China
Updates
Copyright
© 2024 Yuan, Liang, Zhang, Liang, Feng and Dong.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jinguang Liang, liangjingunag@yku.edu.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.