Abstract

A large variety of high-value substation relay protection equipment occupies a considerable amount of inventory space and capital in electric power companies. To improve this problem, this study proposes an inventory prediction model based on the remaining useful life (RUL) of equipment. The model acquires the RUL data of equipment by using the support vector regression (SVR) algorithm, and then, by taking this data as the main factor and the environmental factors and human factors during the operation of equipment as secondary factors, the model can realize the prediction of relay protection equipment in the substation. At the same time, the nature of the enterprise and the requirements for safety inventory are considered. The comparison of calculation results and error analysis, as well as the calculation time, all indicate that the RUL-based inventory forecasting is the best one. This model not only has high prediction accuracy but also has strong stability and portability. The model can provide a strong decision basis for improving the inventory management of the enterprise, enhancing the resource allocation capability, and formulating the spare parts procurement plan under the condition that the spare parts inventory reaches the safety stock.

1. Introduction

Safety stock is an important topic in inventory management [1, 2], and its main role is to meet the uncertainty of supply and demand. The definition of safety stock varies from industry to industry. In power systems, the safety stock of spare parts means that the stock level should always be slightly redundant to the actual demand, in case there are no spare parts available when an unexpected event occurs. Spare parts inventory levels are often calculated based on certain historical data and forecasting models. The authors of [3, 4] have summarized the inventory research results of spares, and the forecasting techniques cover everything from early expert system models to the current artificial intelligence models.

The expert meeting method and the Delphi method are less practical and reliable due to their reliance on expert experience. However, this qualitative model, combined with other methods, has an excellent performance in forecasting [57]. With the development of computers and the demand for short-term and medium-term forecasting techniques, quantitative forecasting models have taken the dominant position, such as gray forecasting models [8, 9], linear regression [10] and nonlinear regression [11], and autoregressive moving average (ARMA) models [1214]. These quantitative forecasting approaches tend to have certain data requirements and also have obvious advantages and disadvantages, such as an emphasis on model improvement to improve the prediction accuracy, but less consideration is given to the factors affecting the inventory, and the information contained in the data is not sufficiently mined. With the development of artificial intelligence, forecasting models based on AI techniques began to play an important role, such as machine learning [15, 16], deep learning [1719], support vector machine (SVM) [20, 21], backpropagation (BP) neural networks [22, 23], and long short-term memory (LSTM) [24, 25]. These intelligent prediction techniques have played an important role in various industries with the help of excellent data processing capabilities.

There is also rich research on inventory forecasting for spare parts in power systems. Zhang [26] established a smart meter inventory demand forecasting model based on analyzing its fault characteristic and installing requirements, and an optimal management strategy [27] was developed. Ding [28] improved the BP neural network model based on the Adam optimization method to forecast the demand for materials in the Guizhou power grid, and the method can significantly reduce the error. In [25], an effective model based on the long short-term memory (LSTM) recurrent neural network was put forward to predict the requirement for maintenance of spare parts. Hamoud and Yiu [29] described a practical reliability model based on a stationary Markov process for assessing the number of spare parts. Yang [30] constructed a deep convolutional neural network based on the graph to realize multicriteria classification for spare parts, which could supply good decision support to control inventory.

As an integral part of the power grid, relay protection equipment has the main characteristics of being diverse, small in size, difficult to obtain monitoring data, and its service life is affected by various factors. Therefore, the advanced prediction techniques mentioned in the above literature, according to our review, found that few prediction models have been applied to this type of equipment. For this reason, in this study, on the basis of analyzing the characteristics of the equipment and its elimination records and safety inventory requirements of electric power enterprises, we make an innovative model of an inventory prediction algorithm according to the RUL of the equipment and applied it to the inventory management system of relay protection equipment in the Changzhou substation. The model implements the following functions:(i)To improve the prediction accuracy of the model, this study fully considers the factors that have direct and indirect effects on the RUL of equipment; these factors are shown in Section 3.1(ii)To reduce the computational difficulty, we define the new variable of the life course as the main variable to obtain the RUL of equipment(iii)To meet the requirement of saving inventory cost and storage space under the condition of achieving inventory safety and to provide decision support for Changzhou Power Supply Company to realize spare parts procurement on a quarterly cycle

Section 1 of this study introduces the current research status of inventory forecasting models and inventory of spare parts in power companies. Section 2 is the inventory forecasting algorithm according to the remaining useful life of equipment, including the algorithm about obtaining the RUL of equipment with the aid of support vector regression and the prediction algorithm of spare parts based on the RUL. Section 3 is a case study for inventory prediction of the central processing unit (CPU) of a substation in Changzhou City. The portability of the model is illustrated by using the predictions for the liquid crystal display (LCD) and DC 220 V Power Supply (DCPS). The validity and stability of the model are indicated by comparing the RUL model with the LSTM, ARMA, support vector machine (SVM), and SVM + BP models. At the same time, the error analyses are given. Section 4 is the conclusions and the perspective.

2. Inventory Prediction Model on the Basis of RUL

The remaining useful life [31] plays an important role in many fields. It can provide strong support for upper-level decision-making, scientifically reduce the operating cost of the system, and enhance the reliability of the system. Especially, with the gradual maturity of AI technology, RUL has been successfully applied to many practical problems [32]. Figure 1 is the main ideation for this study; this framework can demonstrate our research steps.

In this section, we will discuss how to obtain the RUL of substation equipment by using SVR and the construction of an inventory prediction model based on the RUL.

2.1. Principle and Algorithm of RUL Based on SVR

In machine learning [33, 34], SVR is a nonparametric regression model that determines the regression hyperplane by optimizing the distance to nearby support vectors. We first consider the regression hyperplane that achieves the remaining useful life.

For the given training set , since it is difficult to determine whether the training set is linearly divisible in the original space, let the division hyperplane bewhere , is the eigenvector after mapping , is the weight vector, and b is the bias.

The problem of solving the regression hyperplane can be transformed into the following optimization problem:where the penalty coefficient and is the insensitivity loss function and its expression is

Introducing the slack variables and to normalize (2), we can obtain

Introducing the Lagrange multipliers , we construct the Lagrange function of (4) as follows:

Considering (1), let the derivatives of with respect to be zeros; we can obtain

Substituting equation (5)–(8) into (5), we obtain the dual programming problem as follows:

With the KKT condition satisfied, the optimal solution of (1) can be obtained as

(11) is the regression hyperplane used to calculate the remaining useful life.

2.2. Variables and Data Processing

Let the data matrix of a certain type of equipment be , where the th row of denotes the data vector of the th device of that class of devices and N is the number of devices that have been eliminated. The name and type of the variables in are shown in Table 1.

To facilitate the description, several variables are introduced, and the notation and description of the variables are shown in Table 2.

To improve the accuracy of the prediction results, firstly, the data Z need to be cleaned, and let be the days of usage of equipment, and the mean μ and variance of are calculated, and the singular values are screened out using the Z-score method; then, the elimination records containing the singular values are deleted. The formula of the Z-score is as follows:

Generally, the threshold . The data satisfying the above equation are retained, and the new data set matrix is still denoted as Z, where .

Next, we make certain variables labeled and also define a new variable named life course.

The textual variables to are labeled with the regular numerical values and are denoted as to . The numerical variables and are denoted as and . To reduce the dimensionality of the data and simplify the calculation, the life history course are constructed by using the installation date , damage date , and design to life as follows:(i)% M is the number of elimination records of a class of components after cleaning(ii)% I is the number of interval days(iii)(iv)for i from 1 to M(v)j = 1(vi)  Dt =  + I % Date after installation I days(vii)  while now < % now Indicates the current date(viii)   (ix)     =  % Design life course of equipment(x)     =  % Actual life course of equipment(xi)    = () % update variable symbols(xii)   j = j + 1(xiii)   now=now + Iwhere is expressed by the number of days. If the th device is active, take means the current date.

Then, the processed data matrix for this type of device is obtained as . The data of the th device are

2.3. Computation of RUL

Under the above theory and its data characteristics, Algorithm 1 is constructed to predict the RUL of equipment.

Input:  =  .
Output: the RUL value of the kth component .
Step 1. The data matrix is assigned to training set and test set in the ratio of 4 : 1.
Step 2. Get the prediction model by the training data. The regression hyperplane expression is obtained by inputting a certain class of components processed data into equation (2), i.e.,
Step 3. Validity test, let be the number of rows of test data, for the given , if
go to Step4, otherwise, go to (2), and adjust the penalty coefficients and until they are satisfied.
Step 4. Substitute the kth data of the in-use device into (14) and output . This is the value of the life course of the actual life of the kth component.
Step 5. Calculate the damage date of the kth in-service spare part:
Step 6. Compute the number of days left in the life of the kth in-service spare part:

With Algorithm 1, the number of days of the RUL of the kth in-service device can be calculated as

2.4. Algorithm for Inventory Forecasting of Spare Parts by RUL

Once the RUL of an in-service device is obtained, we can construct Algorithm 2 to predict the inventory. The main idea of Algorithm 2 is to consider how many components have an RUL less than T in a prediction period T. Also, a correction function is added for tuning to prevent overprediction and to protect the safety stock.

Input: RUL of the kth equipment
Output: the demand for spare parts Q in the next cycle T.
Step 7. Calculate the number of days remaining in the design life of the kth component
Step 8. Set the correction function to
Step 9. Predicting the inventory level of a certain type of component spare parts, the algorithm pseudocode is as follows:(i)% Q is the demand for a class of components.(ii)% Sum is the number of components of this class being operated in the system(iii)% T is the forecast period(iv)Q= 0(v)for k = 1: sum(vi)If (vii)Q = Q + 1(viii)If and (ix)Q = Q + (x)Output Q % the demand for spare parts in the next cycle.

3. Experiment and Analyses

3.1. Elimination Record

The data on the defect elimination records stored in Changzhou Power Supply Company contain 10 field variables, which are listed in Table 1. The spare parts in this research are installed in relay protection equipment of the substation. They share the following common features:(i)A wide variety and distribution, high value, and integration.(ii)It is difficult to obtain monitoring data since they are board-like devices of small size. The operation and maintenance (O&M) processes require a lot of manpower. It has a long procurement cycle.(iii)Service life is influenced by many kinds of factors listed in Table 1.

3.2. Prediction Results Based on RUL, Comparison, and Its Error

First, we make forecasts for the CPU inventory. Combining with the proposed procurement cycle of Changzhou Power Supply Company, we set the prediction period days and the interval days days and make in the Z-score during data screening and parameter in the correction function . Figure 1 shows the prediction results of CPU spare parts based on the RUL.

In Figure 2, the horizontal coordinate indicates the quarter and the vertical coordinate indicates the prediction result about CPU. From it, we can see that the prediction curve of the CPU with the quarterly cycle has almost the same trend as the real value, which means that the RUL-based prediction method is effective. In the graph, we will also notice that the predicted value is slightly larger than the true value, which is intended to make the spare parts slightly more available in case there are no spare parts available when an unexpected event occurs, i.e., it is beneficial to the stability of the power system.

Second, to illustrate the effectiveness of the algorithm, we employ LSTM, ARMA, SVM, and a combined model of BP and SVM to predict the CPU inventory, and a comparison of the predicted results with the true values is shown in Figure 3.

As we can see from Figure 3, although the prediction results of these four forecasting models are the same as the true values in some quarters, the overall prediction results are not suitable for electric power companies mainly because the difference between the predicted and true values is positive and negative, which leads to the inability of Changzhou Power Supply Company to achieve the goal of purchasing by quarters, and the safe and stable operation of the electric power system cannot be guaranteed under such prediction results.

Finally, the errors of five forecasting models are compared, the root mean square error (RMSE) and mean absolute error (MAE) corresponding to the forecast values under different models are discussed, and the formulas are as follows:where and are the predicted and true values, respectively. If the values of RMSE and MAE of an algorithm are smaller, it means that the algorithm is more effective. The RMSE and MAE of the five method prediction methods are shown in Table 3.

From Table 3, it can be seen that the RUL-based prediction results have the smallest RMSE and MAE among the above five prediction models. It further illustrates the effectiveness of RUL-based CPU prediction.

3.3. Portability and Stability

There are many other CPU-like devices in substations that have similar characteristics to CPUs. To exemplify the portability of the model, we apply the model to the LCD and DCPS inventory predictions, respectively. The parameter settings are unchanged. Figure 4 shows the prediction results based on the RUL model for LCD. Figure 5 presents the inventory prediction results for LCD based on LSTM, ARMA, SVM, and BP-SVM. Figure 6 illustrates the prediction results of DCPS based on the RUL model. Figure 7 demonstrates the inventory prediction results for DCPS based on LSTM, ARMA, SVM, and BP + SVM. Table 4 is the error comparison of the responses.

From Figures 4 and 6, we can notice that the trend of the predicted and true values is almost the same, and the predicted values are also slightly larger than the true values, thus indicating that the RUL-based prediction model has strong portability and also strong stability. This prediction model is suitable for the corporate requirements of Changzhou Power Supply Company and the requirements for power system stability.

From Figures 5 and 7, the predicted and true value curves intersect several times. Although, in some quarters, the predicted results are almost the same as the true ones, this does not satisfy the enterprise requirements.

In Table 4, it can be found that the RMSE and MAE of the forecast results based on RUL are the smallest, except for the MAE of the forecast based on the ARMA model for LCD. For this exception, the result of our analysis is that LCDs are durable equipment in substations with a long service life, resulting in small data for elimination records and data with a strong linear nature, while the ARMA model has a high advantage in dealing with this type of prediction problem.

The purpose of forecasting on a quarterly cycle is to provide a decision basis for the procurement of spare parts for Changzhou Power Supply Company to save storage space as well as improve capital utilization.

3.4. Computation Time

The elapsed time of a model is affected by several aspects, such as the complexity of the algorithm, the dimensionality of the variables, and the code, and therefore, it is often used as a metric to assess the efficiency of the model. In this study, the elapsed time of the five algorithms involved is statistically measured, and the results are shown in Figure 8.

Figure 8shows the time consumption for the computation of the RUL-based prediction model is not the most or the least. Combined with the previous prediction results, we can consider this time consumption to be acceptable for the following reasons: (1) For this study, we are pursuing the accuracy of the prediction and the conformity of the prediction results to the nature of the enterprise; (2) the engineering environment in which the algorithm is applied is not in a time-critical circumstance but has a very sufficient amount of time.

4. Conclusions and Perspective

In this study, a spare parts inventory forecasting model based on the remaining useful life of in-service equipment is proposed for the nature of electric power company enterprises, spare parts characteristics, and data. From the results of the case study, the model is suitable for inventory forecasting in electric power companies. The comparison of the model, the analysis of the error, and the computation time illustrate the validity, stability, and portability of the model.

According to the forecasting algorithm, the quarterly cycle is used for forecasting, which avoids the randomness of data brought by a small cycle and the large forecast error caused by a large cycle. At the same time, the quarterly cycle does not affect the annual procurement plan, while inspiring us to provide a feasible procurement proposal, i.e., to pay and supply according to the quarter, which can save inventory space and reduce holding costs.

The value of the forecasting model used in this study is reflected in three aspects: (1) accurate calculation, portability, and stability; (2) improving inventory management capability and saving inventory cost and storage space; (3) the model idea has certain generalizations, such as applying to the primary variable equipment of power system and enriching the whole life cycle management of power equipment.

The inability to obtain monitoring data makes it difficult to operate and maintain (O&M) relay protection equipment, and usually, routine maintenance does not reduce the failure rate of the device. Currently, the best maintenance strategy is predictive O&M [3537]. An essential prerequisite for predictive O&M is to have information about the status of the health of the equipment, and the remaining useful life is an influential parameter for assessing the operational status of equipment. Therefore, we can utilize RUL to assess the health of equipment and provide a reliable basis for predictive maintenance.

Data Availability

All data included in this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

This paper was supported by the Science and Technology Project of the State Grid Corporation Limited (SGJSCZ00KJJS2100747).