Construction of EMD-SVR-QGA Model for Electricity Consumption: Case of University Dormitory

: In the context of the nationwide call for “energy savings” in China, it is desirable to establish a more accurate forecasting model to manage the electricity consumption from the university dormitory, to provide a suitable management approach, and eventually, to achieve the “green campus” policy. This paper applies the empirical mode decomposition (EMD) method and the quantum genetic algorithm (QGA) hybridizing with the support vector regression (SVR) model to forecast the daily electricity consumption. Among the decomposed intrinsic mode functions (IMFs), define three meaningful items: item A contains the terms but the residual term; item B contains the terms but without the top two IMFs (with high randomness); and item C contains the terms without the first two IMFs and the residual term, where the first two terms imply the first two high-frequency part of the electricity consumption data, and the residual term is the low-frequency part. These three items are separately modeled by the employed SVR-QGA model, and the final forecasting values would be computed as A + B − C. Therefore, this paper proposes an effective electricity consumption forecasting model, namely EMD-SVR-QGA model, with these three items to forecast the electricity consumption of a university dormitory, China. The forecasting results indicate that the proposed model outperforms other compared models.


Introduction
The university dormitory is the place where college students spend most of their daily life. The electricity consumption is often huge without any saving actions. In the context of the nationwide call for "energy savings" in China, it is necessary to provide some effective education courses, guidance, and management approach to cultivate "energy savings" behavior (habit), and eventually to achieve the "green campus" policy. Therefore, it is desirable to establish a more accurate forecasting model to monitor electricity consumption particular from the university dormitory. However, the electricity consumption data of a university dormitory are influenced by lots of factors, such as weather conditions, holidays, population of the dormitory, social activities of university, energy policy of university, and so on [1,2], the electricity consumption data demonstrate the non-linearity in nature, which lead the electricity consumption forecasting work be more complicate [3,4].
In the past decades, many researchers proposed lots of electricity consumption forecasting models to receive more accurate forecasting results. These forecasting models are often divided into two categories [5]: the statistical methods and the artificial intelligent (AI) models. The statistical methods only employ historical data to find out the linear relationships among time periods. There are various statistical models that contain the ARIMA models [6][7][8], regression models [9][10][11], exponential smoothing models [12,13], Kalman filtering models [14,15], Bayesian estimation models [16,17], and so on. However, the inherent shortcomings of these statistical models are that they are only defined to deal with the linear relationships among the electricity consumption and other influenced factors mentioned above, eventually, only receiving unsatisfied forecasting results [18].
As shown in the powerful nonlinear modeling performance, the support vector regression (SVR) model [45] has already been with abundant research results to deal with electric load forecasting [4,[46][47][48][49]. Based on the conclusions of these papers, the well determination of parameters for an SVR model by any meta-heuristic algorithm could guarantee more satisfied forecasting accuracy. Among those employed meta-heuristic algorithms in the literature, the quantum genetic algorithm (QGA) proposed in [50] is considered to be employed in this paper, not only it is empowered by the quantum computing mechanism (QCM), but also it is a classic algorithm with mature applications to solve real problems [51,52]. On the other hand, it is also concluded in the literature that decomposing the time series by the empirical mode decomposition (EMD) method [53] into several single and apparent components (namely intrinsic mode functions (IMFs)), then, separately forecasting each component, it would eventually receive higher forecasting accuracy. The EMD method had also been applied in many fields [53][54][55]. Therefore, this paper would also apply the EMD method to extract the electricity consumption time series data into several IMFs, then, each IMF is forecasted by an SVR model with QGA, and the final forecasting results are summed by the forecasted value of each IMF.
However, as mentioned in [40], the decomposed IMFs contain IMFs with different randomness and the residual IMF; then, these IMFs should be separately modeled by the proposed forecasting model to receive superior forecasting performance. Inspired by the different randomness of IMFs, these two kinds of IMFs could have many combinations to obtain more accurate forecasting results. Thus, based on the EMD method, the QGA, and the SVR model, the authors propose a new hybrid model, namely the EMD-SVR-QGA model, to compare the forecasting performance of different IMF combinations and different algorithms with the H-EMD-SVR-PSO model [40]. The electricity consumption data of a university dormitory is decomposed by the EMD method into six IMFs. Then, these decomposed IMFs are combined into three items: item A contains the terms minus the residual term, item B contains the terms minus the top two IMFs (with high randomness), and item C only contains the terms minus the top two IMFs and the residual term contents. Items A, B, and C would be separately modeled by the SVR-QGA model proposed in [50]. For item A, there are many IMFs with different randomness, which implies that it can have neutral volatility while the SVR-QGA modeling. For item B, the residual term could be fine-tuned under the remaining IMFs (with low randomness), it is also fine while the SVR-QGA modeling. For item C, it is complete with low randomness, which is excellent during the SVR-QGA modeling. Finally, the electricity consumption forecasting results are calculated by the forecasting values of A + B − C. The forecasting details would be presented in the following sections.
To compare the forecasting performance among the proposed EMD-SVR-QGA model and other compared models, the original SVR model, the SVR-QGA model (hybridizing the QGA algorithm with an SVR model), and the H-EMD-SVR-PSO model [40], the electricity consumption data are collected from a university dormitory, Quanshan Campus, Jiangsu Normal University, China, with daily type, from 1 September 2018 to 31 March 2019. The experimental results indicate the superiority of the proposed EMD-SVR-QGA model in terms of forecasting accuracy.
The remainder of this paper is organized as follows. The proposed EMD-SVR-QGA model is introduced in Section 2. The numerical example against other alternative models is illustrated in Section 3. The conclusion is shown in Section 4.

The Empirical Mode Decomposition (EMD) Method
The EMD [53] method is an adaptive time series decomposition method to decompose the original time series with its inherent characteristics into several different stationary IMFs. Each IMF should satisfy the following two criteria: (1) in the whole time range, each IMF at most has one extreme value over; and (2) the mean value of two envelopes by the local maxima and local minima should be zero at any point. Therefore, the EMD method can theoretically overcome the premature convergent problem. For the original time series, ( ), the decomposition procedures of the EMD method are demonstrated as follows.
Step 1: calculation of the mean envelope. Apply the cubic spline interpolation to connect all maxima and minima identified in ( ), to find out the upper envelope and lower envelope, respectively. The mean envelope, , is calculated by the average value of the upper envelope and the lower envelope.
Step 2: decomposing. The first IMF candidate, , is the difference between and the original signal ( ), as shown in Equation (1), Check whether meets the two criteria of an IMF mentioned above. If does not satisfy, then, would be viewed as the original data set, and would be zero. Then, repeat the above evolution procedure for k times, the k-th component, , is shown as by Equation (2), where and ( ) represent the data set after k times and k − 1 times evolution procedures, respectively.
Step 3: identify IMF. If meets the two criteria of an IMF, then, can be identified as the first IMF component, . A new series, , by subtracting IMF1 (as Equation (3)), would continue the decomposition procedure.
Step 4: IMF composition. Repeat steps 1 to 3, until the remained ( ) becomes a monotone function, i.e., there are no new IMFs decomposed from . The decomposition details of these IMFs are shown in Equation (4). As shown in Equation (5), the series, , is the remainder of x(t), i.e., the residual of ( ).

Support Vector Regression Model
The modeling details of the SVR-QGA model are briefly illustrated as follows. For a given nonlinear electricity consumption data set, ( , ) (where represents the actual electricity consumption data) is mapped into a high dimensional feature space by a non-linear mapping function, ( ), where the so-called SVR function, ( ), as shown in Equation (6), is theoretically to formulate the nonlinear relationship between the input training data ( ) and the output data ( ).
where ( ) represents the forecasted values; the weight, w, and the intercept, b, can be computed during the SVR modeling process. Based on the SVR theory, the empirical risk minimization problem can be transformed to solve the quadratic optimization problem with inequality constraints as shown in Equation (7), with the constraints: where ‖w‖ is used to maximize the distance of two separated training data, it also determines the steepness and the flatness of ( ); C is a parameter to penalize the training errors; ε is the width of the ε-insensitive loss function, which defines the loss is zero only if the forecasting value is within the range of ε; two positive slack variables, and ξ * , are used to represent the training statuses, training error above ε, denotes as * , training error below -ε, denotes as .
with the constraints: After solving the quadratic problem, Equation (7), the solution of the weight, w, in Equation (6) is computed by Equation (14), where and * are the Lagrangian multipliers. Eventually, the SVR forecasting function is calculated as Equation (15), where , is the so-called kernel function, and its value could be computed by the inner product of ( ) and ( ) , i.e., , = ( ) • ( ) . Due to superior ability to mapping nonlinear data into high dimensional space, the Gaussian function (as shown in Equation (16)), is employed in this study. Therefore, there are totally three parameters, ε, σ and C, to be optimized in the modeling processes of the Gaussian kernel-based SVR model, to successfully improve the forecasting accuracy.
For comparison with Hong and Fan [40], this paper uses the QGA algorithm to look for suitable parameters of the SVR model.

Quantum Genetic Algorithm (QGA)
Genetic algorithm (GA) is a well-known effective evolutionary algorithm with its advanced operations, including selection, crossover, and mutation operations. It has been applied to solve many optimization problems. However, the worst shortcoming is that its population diversity would decrease along with repeated iterative computation and leads to some problems such as time-consuming, slow convergence, and trapping into local optima. Recently, the quantum computing mechanism [50,56,57] has been applied to be hybridized with the genetic algorithm [50]. The principal quantum computing mechanism, such as qubit, quantum superposition, and quantum entanglement are used to represent the chromosome of QGA into the qubit coding; eventually, the quantum rotation gate is also employed to renew the chromosomes in the modeling process. Therefore, it is with much superiority in faster convergent speed, time-saving, high diversity of population, and robustness.
The quantum computing mechanism is briefly introduced. The smallest information unit in the quantum system is called as qubit (the quantum bit), which demonstrates its state "0", or "1", or any superposition of these two states. The state of a qubit can be as Equation (17), where |0⟩ and |1⟩ represent the values of bits 0 and 1, respectively; and are the probability of the corresponding states and satisfy the normalization criterion, as illustrated in Equation (18), The probability of a qubit individual as a string with qubits is presented as Equation (19), where | | + | | = 1, = 1,2, … , . Then, the chromosome with qubits is presented as, = ( , , … , ) , where ( = 1,2, … , ) denotes an individual qubit as shown in Equation (19). The quantum gate is an operator used for qubits to make the unitary transformations. In general, this operation is often presented by matrices. To obtain more results, it is feasible to use the trigonometric function with a rotation angle , which is so-called quantum rotation gate, as shown in Equation (20). It is used to update as a better state.
where is the updated chromosome; is the rotation angle.

Modeling Processes of QGA
The modeling processes of the SVR-QGA model are completely stated as follows.
Step 1: initialization. Three parameters of an SVR model are initialized as three genes to form a quantum chromosome with a string of m qubits (i.e., the superposition of all possible states), as shown in Figure 1. Based on the practical experiences in [50], a gene contains more qubits could be associated with better partitioning around the space, i.e., could obtain higher forecasting performances. Thus, in this paper, to receive more accurate forecasting results, the chromosome is initialized with 150 qubits (for each gene is 50 bits, i.e., m = 150). The population is initialized by setting all superposition states of qubits with equal amplitude as √ [51]. Step 2: objective fitness evaluation. The objective fitness is evaluated by the quantum chromosome. For convenience, the principal forecasting accuracy index, the mean absolute percentage error (MAPE), as shown in Equation (21), is used to evaluate the objective fitness.
where N is the number of forecasts; is the actual value for forecast point i; is the forecast of point i.
Step 3: selection. The elitist selection is employed to select the best chromosome, which is with the smallest MAPE value during iterations of a generation. The selected elitist is set as the prior chromosome for the next generation.
Step 4: quantum crossover. To maintain the diversity of the population, the quantum crossover operation is used. This paper also uses the single-point-crossover to implement the crossover operation by two randomly selected chromosomes with crossover probability, Pcr, which is suggested to be set as 0.9 [50].
Step 5: quantum mutation. To efficiently enrich the diversity of the population, the quantum mutation operation is also employed. Each selected quantum chromosome will be changed by other suitable oscillation amplitudes ( , ) of the corresponding state with the mutation probability, Pm, which is suggested to be set as 0.1 [50].
Step 6: quantum rotation. To avoid premature convergence, it is useful to change the quantum location of a qubit, therefore, the quantum rotation is used to modify the oscillation ranges of the qubit. The quantum rotation gate (shown in Equation (20)) with the rotation angle, , is practically applied to perform the quantum rotation. The rotation angle, , is corresponding to ( , ) to modify the quantum location of a qubit. The values of are suggested to be between 0.1π and 0.005π [52].
Step 7: stop criteria. If the required number of generations is reached, then, the best solution is identified as the last quantum chromosome, else, go back to Step 3 to continue the three operations for the next generation.

The Complete Processes of the Proposed EMD-SVR-QGA Model
The complete processes of the proposed EMD-SVR-QGA model are presented in Figure 2 and stated as following steps. Step 1: conduct decomposition of the training data. The electricity consumption data from the university dormitory would be decomposed by the EMD method into several characterized IMFs. These decomposed IMFs would be combined into three items: item A includes all IMFs but without the residual term; item B includes IMFs without top two (the first randomness and the second randomness which are with highly random, as shown in Figure 2) and the residual term; item C includes IMFs without top twos.
Step 2: defined items modeled by the SVR-QGA model. These three defined items, A, B, and C, are separately modeled by the SVR-QGA model. The modeling processes of the SVR-QGA model, such as quantum chromosome setting, superposition states of qubits, population initialization, selection operation, quantum crossover operation, quantum mutation operation, and quantum rotation operation, could refer Section 2.2.3 to learn about more details of the SVR-QGA model.
Step 3: final forecasting results. The final electricity consumption forecast is calculated by the values of these three items, A, B, and C, modeled separately by their corresponding SVR-QGA models. The calculation formula is easily given by the forecasts of A + B − C.

Data Set of Experimental Examples
The electricity consumption data are collected from a real-world case, two representative university dormitories, Quanshan Campus, Jiangsu Normal University, China. The detailed background of these data is stated as follows. There are about 12 persons per dormitory room. On the campus, each week about 10 seminars and associated lectures, each seminar or lecture would attract about 200 persons. The relative weight of Dormitories' consumption in the corresponding university campus total consumption is about 14.78%. Based on the time period of the collected data, the associated weather is in fall and winter, as the campus is located in northern China, the official heating system is available. However, if suffering from heavy snow or intense cold, the electrical heater is also used. Thus, the electricity consumption would be moderate fluctuated along with the weather conditions. Finally, due to "energy savings" activities in China, the administration department also announces some active policy to encourage "energy savings", eventually, to achieve the "green campus", such as record merits, award, and give honorary title, etc.
It would be further processed by the proposed EMD-SVR-QGA model to receive the forecasting electricity consumption of these two selected dormitories in the future for one month. For these two dormitories, the electricity consumption data are both collected, with daily type, from 1 September 2018 to 31 March 2019, thus, each data set contains 212 electricity consumption data. However, during the vacation period, such as the Chinese Lunar Year Festival and the winter vacation, the electricity consumption of each dormitory is almost zero, i.e., these data should be removed. Therefore, after pre-processing operation, these two data sets from two dormitories are arranged to 110 and 134 electricity consumption data, respectively.
For clear denotation of these two dormitory cases, the first dormitory with 110 electricity consumption data is denoted as dormitory I, the other one with 134 electricity consumption is denoted as dormitory II. For the division of the training and the testing sets, Dormitory I is with 88 and 22 electricity consumption data in the training and the testing sets, respectively; dormitory II is with 107 and 27 electricity consumption data in the training and the testing sets, respectively.
Besides two dormitory cases, the electric load data set, collected from New South Wales (NSW) market in Australia also used in [40], is also employed to demonstrate the superiority of the proposed EMD-SVR-QGA model. The employed electric load data is collected from 1 to 13 August 2019, based on half-hour format, i.e., there are 48 data per day. The training set is from 1 to 10 August 2019 (in total 480 load data points), and the testing data is from 10 to 13 August 2019 (in total 144 load data points).
During the training stage of the SVR-QGA model, the rolling-based process [59,60], is employed to help QGA to systematic look for suitable parameter combination of an SVR model with smaller forecasting error. Specifically, the first n electricity consumption data is fed into the SVR-QGA model to conduct the model training process, then, the (n + 1)-th forecast value of electricity consumption and the feasible parameter combination are both obtained by the SVR-QGA model. Repeat this operation to receive the (n + 2)-th, the (n + 3)-th, … forecast values of electricity consumption till the required forecasting values are received. Finally, the parameter combination with the smallest testing errors are verified as the most appropriate parameter combination of an SVR model. Eventually, the testing data set will be used to compute the forecasting accuracy of the proposed SVR-QGA model.

Parameters Setting of the EMD-SVR-QGA Model
The parameters of QGA for these two dormitory cases are set practically or by referencing other experiences, the population scale (Pscale) is set to be 200 [50]; the generations of the population (qmax) are no larger than 200; the qubit string length of a quantum chromosome (m) is set as 150; the probabilities of quantum crossover (Pcr) and quantum mutation (Pm) are set as 0.5 and 0.1 [50], respectively. The rotation angle, , is between 0.1π and 0.005π [52]. For SVR-QGA modeling, the maximal iteration for each example is all set as 10,000 in each generation; the searching range of the three parameters are set as followings, σ ∈ 0,10 , ε ∈ 0,1 , and ∈ 0,200 in both two dormitory cases; δ is fixed as 0.001 [50].

Forecasting Accuracy Indexes and Forecasting Performance Superiority Test
To compare the forecasting accuracy among the proposed model and other compared models, four forecasting accuracy indexes, the mean absolute percentage error (MAPE), the root mean squared error (RMSE), mean squared error (RMSE), and the mean absolute error (MAE) is employed to calculate the forecasting errors, they are defined as Equations (21)-(24).
where N is the number of forecasts; is the actual value for forecast point i is the forecast of point i.
On the other hand, it is necessary to conduct the statistical test to verify the statistical significance of the forecasting superiority of the proposed model. Based on Derrac et al. [61] suggestions, only Wilcoxon signed-rank test is implemented in this paper. The statistic of Wilcoxon signed-rank, W, is defined as Equation (25) where (+) is the total number of positive rank that the forecasting error of the proposed EMD-SVR-QGA model larger than other compared model; (−) is the total number of negative rank that the proposed EMD-SVR-QGA model smaller than the compared model. If the value of W satisfies the required value of Wilcoxon distribution with an associated degree of freedom, then, the null hypothesis of equal performance of these two compared models is not accepted. It also reveals the outstanding forecasting results of the proposed EMD-SVR-QGA model outperforms significantly the compared model.

Decomposition Results and Forecasting Results for Dormitory I
For Dormitory I case, the electricity consumption data is firstly decomposed by the EMD method into six terms. These six decomposed terms are demonstrated in Figure 3, in which the first IMF and the second IMF are obviously the random terms, the 6th IMF is also obviously the residual term. The three items A, B, and C, defined by the six IMFs are shown in Figure 4.
The three items A, B, and C, defined by the six IMFs are separately modeled by the SVR-QGA model, as mentioned in Section 3.1, the rolling-based process [59] is applied to collaborate with the QGA to determine the most appropriate parameter combination of their associated SVR-based models in the training phase. The appropriate parameter combination for each item is shown in Table 1. In addition, the appropriate parameter combination for other compared models, including the original SVR model, the SVR-QGA model, the H-EMD-SVR-PSO model [40], are also illustrated in Table 1. The forecasting performances for the three items are demonstrated in Figure 5.  The forecasting results of the original SVR model, the SVR-QGA model [50], the H-EMD-SVR-PSO model [40], and the proposed EMD-SVR-QGA model are shown in Figure 6. It illustrates the forecasting curve of the proposed EMD-SVR-QGA model has better fitting effects than other compared models. It also implies that the proposed EMD-SVR-QGA model has the potentials to capture the fluctuation variation of the electricity consumption from this dormitory. Thus, the consumption behaviors can be effectively monitored and some valuable activities regarding improving electricity consumption habits can be held to successfully save electricity consumption.  The proposed EMD-SVR-QGA model also demonstrates better generalization ability than other models. The forecasting performance of these models is listed in Table 2. In which, it demonstrates clearly that the proposed EMD-SVR-QGA model outperforms the original SVR model, the SVR-QGA model [50], and the H-EMD-SVR-PSO model [40] in terms of four employed forecasting accuracy indexes. Particularly, for the H-EMD-SVR-PSO model, it also receives outstanding forecasting results, which indicates the forecasting advantages from this kind of EMD-SVR-based model. Comparing the forecasting accuracy of the EMD-SVR-QGA model and the H-EMD-SVR-PSO model, their contributions are different due to different combination types of these decomposed IMFs and the residual term. Therefore, looking for a suitable IMFs combination approach is an interesting research topic. In addition, in Table 2, it also can see that the forecasting accuracy of the SVR-QGA model [50] is not superior to the EMD-SVR-based models. This is because the nonlinear effects are inherent in the electricity consumption data itself, and demonstrate the interactions among the IMFs and the residual term. After decomposition operation by the EMD method, the EMD-SVR-based models are able to well deal with this inherent non-linear data pattern by the IMFs and these defined items, A, B, and C. Therefore, the proposed EMD-SVR-QGA model can be a useful forecasting model to deal with the electricity consumption forecasting work of a university dormitory. Finally, to verify the significant contribution of the forecasting improvement from the proposed EMD-SVR-QGA model, in this paper, the Wilcoxon signed-rank test is suggested [60]. Wilcoxon signed-rank test is conducted on one-tail-test, under two significant levels, α = 0.025 and α = 0.05. The test results are shown in Table 3, which obviously indicate that the proposed EMD-SVR-QGA model is significantly superior to other compared models.

Decomposition Results and Forecasting Results for Dormitory II
Similar to the dormitory I case, the electricity consumption data in the dormitory II case is also firstly decomposed by the EMD method into six terms. These six decomposed terms are demonstrated in Figure 7, in which the 1st IMF and the 2nd IMF are obviously the random terms, the 6th IMF is also obviously the residual term. The three items A, B, and C, defined by the six IMFs are shown in Figure 8. The decomposed IMFs and the residual term for the dormitory II case. IMF (the most random term); second IMF (the second most random term); thirrd IMF (general term 1); fourth IMF (general term 2); fiftth IMF (general term 3); sixth IMF (the residual term). In the dormitory II case, the three items A, B, and C, defined by the six IMFs are also separately modeled by the SVR-QGA model to determine the most appropriate parameter combination of their associated SVR-based models. The appropriate parameter combination for each item for other compared models, including the original SVR model, the SVR-QGA model, the H-EMD-SVR-PSO model [40], are illustrated in Table 4. The forecasting performances for the three items are demonstrated in Figure 9.
For the dormitory II case, the forecasting results of the original SVR model, the SVR-QGA model [50], the H-EMD-SVR-PSO model [40], and the proposed EMD-SVR-QGA model are shown in Figure 10. It also demonstrates the forecasting curve of the proposed EMD-SVR-QGA model fits better with the actual electricity consumption than other compared models. It also implies that the proposed EMD-SVR-QGA model is qualified to well deal with the fluctuation characteristics of the electricity consumption from this dormitory. Thus, it could contribute to the campus managers providing more insight governance policy for the electricity-saving of the dormitory.   The forecasting performance of these models is listed in Table 5. As the same results in the dormitory I case, the proposed EMD-SVR-QGA model also shows better generalization ability than other models, the original SVR model, the SVR-QGA model [50], and the H-EMD-SVR-PSO model [40] in terms of four employed forecasting accuracy indexes. The forecasting advantages from this kind of EMD-SVR-based model are revealed once again. Different combination types of these decomposed IMFs and the residual term could obtain different forecasting accuracy. It is deserved to look for a suitable IMFs combination approach.
Furthermore, from Table 5, the forecasting results support the same conclusion in Dormitory I case once again: the EMD-SVR-based models able to well deal with this inherent non-linear data pattern by the IMFs and these defined items, A, B, and C. Therefore, the proposed EMD-SVR-QGA model can be a useful forecasting model to deal with the electricity consumption forecasting work of a university dormitory. For the significance test of the forecasting performance from the proposed EMD-SVR-QGA model, the results of Wilcoxon signed-rank test, based on one-tail-test and is under two significance levels, α = 0.025 and α = 0.05, are shown in Table 6. Obviously, it indicates that the proposed EMD-SVR-QGA model outperforms significantly other compared models.

Decomposition Results and Forecasting Results for New South Wales (NSW, Australia) Market
Similar to the above two dormitory cases, the electric load data in the New South Wales (NSW, Australia) market is also firstly decomposed by the EMD method into six terms. These six decomposed terms are demonstrated in Figure 11, in which the first IMF and the second IMF are obviously the random terms, the sixth IMF is also obviously the residual term. The three items A, B, and C, defined by the six IMFs are shown in Figure 12. Figure 11. The decomposed IMFs and the residual term for NSW (Australia) case. IMF (the most random term); second IMF (the second most random term); third IMF (general term 1); fouth IMF (general term 2); fifth IMF (general term 3); sixth IMF (residual term). In the dormitory II case, the three items A, B, and C, defined by the six IMFs are also separately modeled by the SVR-QGA model to determine the most appropriate parameter combination of their associated SVR-based models. The appropriate parameter combination for each item for other compared models, including the original SVR model, the SVR-QGA model, the H-EMD-SVR-PSO model [40], are illustrated in Table 7. The forecasting performances for the three items are demonstrated in Figure 13.
For NSW (Australia) case, the forecasting results of the original SVR model, the SVR-QGA model [50], the H-EMD-SVR-PSO model [40], and the proposed EMD-SVR-QGA model are shown in Figure 14. It also demonstrates the forecasting curve of the proposed EMD-SVR-QGA model fits better with the actual electric load than other compared models. It also implies that the proposed model is qualified to well deal with the fluctuation characteristics of the electric load in real-world applications, such as the NSW market in Australia. Thus, it supports that the proposed model is useful to the campus managers to receive accurate forecasting performance.   The forecasting performance of these models is listed in Table 5. As the same results in dormitory I case, the proposed EMD-SVR-QGA model also shows better generalization ability than other models, the original SVR model, the SVR-QGA model [50], and the H-EMD-SVR-PSO model [40] in terms of four employed forecasting accuracy indexes. It once again reveals the superiority of EMD-SVR-based models.
Furthermore, from Table 8, the forecasting results also support the same conclusion in NSW (Australia) case: the EMD-SVR-based models able to well deal with this inherent non-linear data pattern by the IMFs and these defined items, A, B, and C. Therefore, the proposed EMD-SVR-QGA model is appropriate to be applied to deal with the electricity consumption forecasting work. For the significance test of the forecasting performance from the proposed EMD-SVR-QGA model, the Wilcoxon signed-rank test is conducted as in the previous two cases, also based on one-tail-test and two significance levels, α = 0.025 and α = 0.05, the tested results are shown in Table  9. Obviously, it indicates that the proposed EMD-SVR-QGA model outperforms significantly other compared models.

Discussions
As shown in [40], they defined three items (A, B, and C) based on the excellent combination of these different decomposed IMFs, such as the random term, the middle terms, and the residual term. Therefore, a suitable combination of the decomposed IMFs, i.e., how to make good use of these IMFs, could be an interesting issue to help to improve the forecasting accuracy of the proposed model. However, due to the inherent complexity of the employed data set, the decomposition results are varied. For example, in [40], there is only one decomposed IMF with random volatility and the best ones are not, thus, the best combination of these IMFs is demonstrated in [40]. In addition, there may be two or more decomposed IMFs with random volatility, such as in this paper there are two terms with different randomness, therefore, we define different three items, i.e., different combinations of these decomposed IMFs (the proposed EMD-SVR-QGA model). Eventually, we receive more accurate forecasting results than the results using the combination approach (the H-EMD-SVR-PSO model) in [40]. We would like to claim that the superiority of our proposed EMD-SVR-QGA model than the H-EMD-SVR-PSO model. Therefore, it could be remarked that an effective combination of the decomposed IMFs could provide a clearer cue to understand the complex structure of the employed data set, then, make good use of these IMFs can help to receive more satisfying analysis results, such as forecasting accuracy improvement or classification accuracy enhancement. The exploration of the effective combination of the decomposed IMFs also inspires some interesting future research.

Conclusions
This paper proposes a novel EMD-SVR-QGA electricity consumption forecasting model to provide the campus managers with more accurate electricity consumption forecasting from the university dormitory. It is superior in capturing the fluctuation variation of the electricity consumption and reveals its potentials to indicate the daily patterns of the electricity consumption of the dormitory which is useful to take some valuable activities regarding improving electricity consumption habits. The proposed EMD-SVR-QGA model uses different combinations approach of the IMFs decomposed by the EMD method, to well deal with the interaction among those decomposed IMFs and the residual term and the inherent complexity in the electricity consumption data. Furthermore, this paper also applies the QGA to comprehensively search for the appropriate values of the three parameters of an SVR model, then, used the trained SVR-QGA model to forecast the defined items (A, B, and C), separately. It ensures receiving more satisfied forecasting results. Finally, two representative university dormitories' electricity consumption data sets are used to conduct modeling, the experimental results demonstrate that the proposed model has outperformed significantly other compared models in terms of forecasting accuracy indexes, such as the original SVR model, the SVR-QGA model, and the H-EMD-SVR-PSO model.
The experimental results also indicate the very issues that the different combination of these decomposed IMFs to deal with the inherent nonlinear or other complex interactions among the data set. This issue also asks researchers to carefully study the characteristics of the IMFs, then, determine an effective combination type of these IMFs and the associated SVR-based model. For future research, we will consider hybridizing different intelligent technologies to combine the IMFs, and other useful meta-heuristic algorithm with quantum computing mechanism, to enrich the research contents to receive higher forecasting accuracy.