Next Article in Journal
Review and Evaluation of Multi-Agent Control Applications for Energy Management in Buildings
Next Article in Special Issue
Evaluation and Long-Term Prediction of Annual Wind Farm Energy Production
Previous Article in Journal
Electromagnetic Energy Harvester Using Pulsating Airflows—Reeds Waving in the Wind
Previous Article in Special Issue
Comprehensive Review of Building Energy Management Models: Grid-Interactive Efficient Building Perspective
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Advanced Optimal System for Electricity Price Forecasting Based on Hybrid Techniques

1
School of Economics and Finance, Shanghai International Studies University, Shanghai 201620, China
2
School of Statistics, Dongbei University of Finance and Economics, Dalian 116025, China
*
Author to whom correspondence should be addressed.
Energies 2024, 17(19), 4833; https://doi.org/10.3390/en17194833
Submission received: 9 August 2024 / Revised: 8 September 2024 / Accepted: 24 September 2024 / Published: 26 September 2024
(This article belongs to the Special Issue Advanced Forecasting Methods for Sustainable Power Grid)

Abstract

:
In the context of the electricity sector’s liberalization and deregulation, the accurate forecasting of electricity prices has emerged as a crucial strategy for market participants and operators to minimize costs and maximize profits. However, their effectiveness is hampered by the variable temporal characteristics of real-time electricity prices and a wide array of influencing factors. These challenges hinder a single model’s ability to discern the regularity, thereby compromising forecast precision. This study introduces a novel hybrid system to enhance forecast accuracy. Firstly, by employing an advanced decomposition technique, this methodology identifies different variation features within the electricity price series, thus bolstering feature extraction efficiency. Secondly, the incorporation of a novel multi-objective intelligent optimization algorithm, which utilizes two objective functions to constrain estimation errors, facilitates the optimal integration of multiple deep learning models. The case study uses electricity market data from Australia and Singapore to validate the effectiveness of the algorithm. The forecast results indicate that the hybrid short-term electricity price forecasting system proposed in this paper exhibits higher prediction accuracy compared to traditional single-model predictions, with MAE values of 7.3363 and 4.2784, respectively.

1. Introduction

In recent years, renewable energy generation, including wind and solar photovoltaics, has seen rapid expansion and active integration into the electricity market. Due to the uncontrollability of their power generation, it breaks the real-time supply and demand balance of electric energy, which causes the fluctuation of electricity prices in the power market to become more severe [1,2]. Electricity prices are considered an important reference for decision making in the market. Accurate short-term electricity price forecasting can help market participants and power suppliers flexibly adjust pricing strategies according to changes in electricity prices, thereby increasing profits [3]. Moreover, electricity purchasers can use the difference between peak and off-peak electricity prices to carry out dynamic cost management to save expenses. At the same time, it also provides an important scientific basis for the real-time supervision of regulatory authorities, ensuring the normal operation of electricity [4].
Our study focuses on enhancing the precision and robustness of hybrid models for forecasting electricity prices, employing cutting-edge optimization techniques. The motivation for our research lies in the potential impact of reliable short-term electricity price forecasts on the energy market, fostering intelligent data-driven strategies beneficial to energy managers, suppliers, and consumers, especially in an era of rapid technological advancements and severe resource wastage. The essence of our research lies in developing time-series forecasting models suitable for further academic investigation, thereby underscoring the goal of the rigorous evaluation and enhancement of current predictive methodologies.
Electricity prices exhibit non-stationary and variable characteristics within the energy market, coupled with notable anomalies. Such factors introduce increased uncertainty and complexity to the electricity market, complicating the task of short-term price forecasting [5,6]. Generally, the prevailing approaches to forecasting encompass statistical methods, along with models based on machine learning or deep learning, and those that integrate hybrid methodologies.
Models based on statistical analysis usually analyze historical data mathematically and statistically to build probabilistic models to predict future trends. Examples include the hidden Markov model (HMM) [7], autoregressive integrated generalized autoregressive conditional heteroskedasticity (Autoregressive-GARCH) model [8], and autoregressive integrated moving average (ARIMA) model [9]. These models are simple to implement, require fewer parameters for adjustment, and demand fewer computing resources. They can predict the trend in electricity price changes to a certain extent. However, their effectiveness is diminished by their reliance on assumptions, such as data smoothness and autocorrelation, which render them less capable of accommodating a nonlinear series.
To address the limitations inherent in statistical analysis models, scholars have increasingly shifted their focus toward machine learning and deep learning methodologies [10]. These techniques have proven to be highly effective in identifying and processing complex, nonlinear characteristics, thus offering superior performances in applications. Machine learning models mainly include artificial neural networks (ANN) [11], support vector machine (SVM) [12], and random forest (RF) [13]. Machine learning can adjust weights to approximate multivariate functions and capture complex, dynamic, nonlinear features of electricity prices. However, these methods require a significant amount of historical data and effective features. Additionally, the algorithm’s robustness needs improvement. More researchers opt for deep learning models when it comes to forecasting tasks [14]. Typical studies include the following: Rezaei, Rajabi, & Estebsar (2022) utilized a Stacked Auto Encoder (SAE) to extract features from denoised electricity price data efficiently. To further refine the accuracy of electricity price forecasts and minimize the potential for overfitting, this model incorporated a Gated Recurrent Unit (GRU) [15]. And B. Wang, Wei, & Su (2022) employed Long Short-Term Memory (LSTM) networks, which led to the selection of more accurate forecast samples, resulting in a reduction in both relative and absolute forecasting errors [16]. Furthermore, three different artificial intelligence algorithms were applied to forecast day-ahead electricity prices in the Turkish electricity market [17]. The experimental findings demonstrated that, in the context of error metrics, the model based on Extreme Gradient Boosting Decision Trees (XGBoost) exhibited superior performance compared to its counterparts. Deep learning approaches typically outperform traditional machine learning algorithms and encounter notable difficulties [18]. However, the process of directly incorporating exogenous variables through partial models risks overwhelming the system and obscuring recognizable patterns. Furthermore, the conventional method of initializing weights randomly in deep learning frameworks might underemphasize vital features and moments, thereby reducing the effectiveness of feature extraction. This necessitates a refined strategy to ensure balanced weight distribution and optimal pattern discernment, enhancing the overall efficacy of models.
Researchers have employed a variety of hybrid strategies to address the limitations inherent in singular models by leveraging the unique advantages of different methodologies to “learn from the strong to offset one’s weaknesses effect” [19,20]. These efforts can be categorized into three distinct groups.
The first group focuses on optimizing from a singular model perspective by combining singular models with optimization models or by integrating singular models with variable selection models [21]. For instance, Qu et al. (2024) introduced an improved Wavelet Neural Network (WNN) model, building upon the foundation laid by ELM to construct classification prediction models for various electricity pricing schemes, effectively overcoming the traditional WNN issues of slow or even non-converging rates [22]. Similarly, J. Wang, Liu, Song, & Zhao (2016) updated the dynamic choice artificial neural network (DCANN) using the Cuckoo Search algorithm (CS) to address the initial parameter shortcomings of DCANN [23]. This approach also involved the identification of poor training sample sets through the Index of Bad Samples Matrix (IBSM) and the selection of appropriate inputs for each desired output. While such models have achieved commendable predictive performance, they struggle to capture the complex characteristics of electricity price series.
To address this issue, a second group of models was introduced, incorporating signal decomposition strategies to forecast different trend variables independently and subsequently integrate the final predictions. Typical studies include the following: H. Zhang et al. (2021) introduced an innovative hybrid model that combined Temporal Convolutional Networks (TCNs) and ARIMA with empirical mode decomposition (EMD) [24]. This method targets high- and low-frequency data parts and significantly enhances accuracy. However, EMDs may experience modal aliasing at different frequencies when processing signals, especially when the signals contain similar frequency components. To address this challenge, P. Jiang, Nie, Wang, & Huang (2023) implemented the Improved Complete Ensemble EMD with Adaptive Noise (ICEEMDAN) technique in their research [25]. This approach ameliorates the mode mixing problem by introducing additional noise signals. However, the selection of noise levels within the algorithm significantly impacts the outcome, and the lack of a unified standard for choosing the optimal noise level complicates the utilization of this method. Conversely, variational mode decomposition (VMD) operates through predetermined bandwidth constraints to decompose signals, effectively reducing mode mixing between different modes [26]. The decomposition results are less susceptible to signal noise, enhancing the precision and interpretability of the decomposition [27]. Compared to EMD and its variants, VMD requires fewer parameters, making the algorithm easier to adjust and implement [28]. The predictive accuracy of second-group hybrid models demonstrates an improvement, yet they struggle to surmount the intrinsic limitations of singular models and lack the capacity for adaptive parameter adjustments in response to variations in input samples.
Consequently, contemporary research has expanded to encompass a third group of models that advocate for the integration of multiple predictive models, employing optimization algorithms to facilitate a weighted amalgamation of diverse models. Z. Yang, Ce, & Lian (2017) introduced the hybrid methodology aimed at enhancing predictive accuracy by leveraging the distinct characteristics of each model [29]. Specifically, our approach integrated the wavelet transform with the kernel extreme learning machine (KELM), which was further optimized through self-adapting particle swarm optimization techniques and the autoregressive moving average (ARMA) model. X. Zhang, Wang, & Gao (2019) assessed a groundbreaking hybrid forecasting framework in the context of the New South Wales electricity market [30]. This framework innovatively melds a variety of predictive methodologies, encompassing seasonal ARIMA, Backpropagation (BP) neural networks, and SVM, aiming to bolster the precision of forecasts. Y. Xu, Li, Wang, & Du (2024) adeptly utilized an advanced version of the Multi-Objective Tuna Swarm Optimization (MOTSO) algorithm [31]. This innovative approach was pivotal in determining the combined weights assigned to each model. The third group of models had significantly improved predictive accuracy; however, the selection of the model and its parameter configuration greatly influenced the operational efficiency and predictive performance. Thus, it is imperative to engage in additional simulation experiments to identify more optimal methods for short-term electricity price forecasts. Moreover, conducting multi-step forecasting experiments is crucial for maintaining the integrity of the research framework. The representative articles, models, and findings discussed in the preceding literature review are systematically outlined and summarized in Table 1.
In reviewing the literature, it is evident that hybrid forecasting methods are gaining acknowledgment for their significant advantages in predicting electricity prices. However, despite their rising importance, these models have limitations. The following gaps exist in current research: (1) The electricity price dataset is extensive and rife with extraneous information. Conventional time-series decomposition methods frequently generate mixed-frequency components, compromising the precision of predictive analyses. (2) In the realm of model combination, the absence of a scientifically rigorous weighting scheme can yield less-than-ideal forecasting outcomes. Moreover, the utilization of antiquated weighting strategies may impede optimal model parameter tuning and prolong the computational runtime.
Therefore, this study employed an innovative approach to address the existing gaps in research by integrating trend decomposition techniques with multi-objective intelligent optimization algorithms, effectively optimizing prediction accuracy. The principal contributions of the study are as follows:
(1) Introduce an innovative trend decomposition technique that enhances the efficiency of feature extraction for predictive analysis. By employing an improved variational mode decomposition, we can dissect electricity price sequences into trend components, cyclical components, and irregular fluctuating components. This allows for the separate prediction of each component, facilitating a more effective capture of the distinct characteristics inherent in electricity price data.
(2) Apply a novel multi-objective intelligent optimization technique to further refine the predictive outcomes of deep learning models through scientific weighting, thereby increasing the forecasting accuracy for non-stationary sequences. By constraining the estimation error through an objective function and seeking Pareto optimal solutions, we promote the integration of multiple deep learning models, leveraging the unique predictive strengths of each model.
(3) The proposed hybrid short-term electricity price forecasting system demonstrates adaptability and scalability, making it applicable to various types of power markets. Through cases and discussions conducted across 14 comparative models in markets experiencing different levels of volatility, we validated the effectiveness of the proposed model’s forecasting capabilities.
This paper is structured as follows: Section 1 is the introduction presenting the research motivation, related literature, current research gaps, and innovations; Section 2 introduces the basic concepts of model preparation; Section 3 demonstrates the main framework of this study; Section 4 gives the details of the experimental setup, data presentation, numerical results, and case study; Section 5 discusses the results of the present study; and Section 6 presents the conclusion and future research plans.

2. Basic Concepts of Model Preparation

2.1. Trend Decomposition Technique

Electricity price series exhibit complexities such as instability, nonlinearity, and unpredictability, posing a significant challenge to accurate forecasting [32]. Decomposition algorithms become essential in addressing prediction inaccuracies due to noise, elucidating crucial electricity price attributes, and augmenting the efficacy of predictive models in various domains [33]. The essence of selecting apt decomposition algorithms lies in their ability to discern and utilize the principal characteristics of electricity prices efficiently. However, traditional signal decomposition methods, such as Fourier transform and wavelet transform, may have limitations in some cases. For example, they usually assume that the signal is linear and stationary, which is less effective when dealing with nonlinear and non-stationary signals. In order to overcome these limitations, K. Dragomiretskiy & K. Zosso (2014) proposed the variational mode decomposition (VMD) algorithm [34].
VMD stands out as a premier technique for the iterative separation of multi-component signals into the quasi-orthogonal intrinsic mode function. VMD is notably superior to the adaptive decomposition method EMD as it circumvents the issues of recursive computation errors and premature termination. Therefore, for the analysis of electricity price data, VMD emerges as the methodology of choice [35].
The price sequence γ ( t ) can be decomposed into P sub-signals (modes) { E p ( t ) } = { E 1 ( t ) , , E P ( t ) } , where each mode is a set of band-limited data with a center frequency. { ω p } = { ω 1 , , ω P } and E p ( t ) ae able to be transformed into their analytical form by E A , p ( t ) = E p ( t ) j H E p ( t ) = ( t ) + j / π t * E p ( t ) , where H is the Hilbert transform operator, δ ( ) represents the use of the Dirac function, and * embodies the convolution.
To derive the E p ( t ) in the frequency domain, we can recenter the spectrum of E p ( t ) and apply a low-pass filter. The formulation of the recenter process is E A , p ( t ) e j ω p t = δ ( t ) + j / π t * E p ( t ) e j ω p t .
Following Fourier transform, e j ω p t should be converted to δ ( ω + ω p ) in the frequency domain. The bandwidth of E p ( t ) can be determined by optimizing the demodulated spectrum. Thus, the objective and constraint are formed in (1).
min t δ ( t ) + j / π t * E p ( t ) e j ω p t 2 2 s . t .   E p ( t ) = γ ( t )
The utilization of the augmented Lagrangian method ( L ) assists in transforming the aforementioned constrained optimization problem to an unconstrained form:
L = α t δ ( t ) + j / π t * E p ( t ) e j ω p t 2 2 + γ ( t ) E p ( t ) 2 2 + λ [ γ ( t ) E p ( t ) ]
where α is the balancing parameter of the data-fidelity constraint, and λ is the Lagrangian multiplier.
To simplify (2), a more feasible approach would be to transpose L into the frequency domain using the Fourier transform as follows (3):
min α j ω 1 + sgn ( ω + ω p ) E p ( ω + ω p ) 2 2 + γ ( ω ) E p ( ω ) + λ ( ω ) / 2 2 2
where the multiplication of the exponential term becomes the frequency spectrum shift operation, the convolution becomes multiplication, the derivative operation becomes multiplying the objective function with j ω , and sgn means the sign function.
Employing the Alternating Direction Method of Multipliers (ADMMs) [36], we computed the spectral properties of mode E P ( ω ) by executing a series of iterative refinement with:
E P r + 1 ( ω ) = argmin E P α j ω ( 1 + sgn ( ω + ω P r ) ) E P r ( ω + ω P r ) 2 2 + γ ( ω ) l = 1 p 1 E l r + 1 ( ω ) E l r ( ω ) + λ r ( ω ) / 2 2 2
The iteration number is denoted by r, where l is mode index, P is the total number of modes, and Hermitian symmetry is applied to simplify the following:
E P r + 1 ( ω ) = argmin E P 0 4 α ω ω p r 2 E p r ( ω ) 2 + 2 γ ( ω ) l = 1 p 1 E l r + 1 ( ω ) l = p P E l r ( ω ) + λ r ( ω ) / 2 2 d ω
The optimal solution is letting the derivative equal 0. Thus,
E p r + 1 ( ω ) = γ ( ω ) l = 1 p 1 E l r + 1 ( ω ) l = p P E l r ( ω ) + λ r ( ω ) / 2 1 + 2 α ω ω p r 2
The center frequency is ω p r + 1 = argmin ω p { 0 ( ω ω p r ) 2 | E P r + 1 ( ω ) 2 | d ω } . The solution is found by the following equations: ω p r + 1 = ( 0 ω E p r + 1 ( ω ) 2 d ω ) / ( 0 E p r + 1 ( ω ) 2 d ω ) and λ r + 1 ( ω ) = λ r ( ω ) + τ ( γ ( ω ) E p r + 1 ( ω ) ) , where τ is the penalty parameter. And they run until the convergence criterion v E p r + 1 ( ω ) E p r ( ω ) 2 2 / E p r ( ω ) 2 2 ε is met.
Following a thorough examination of the data, this study employed a K-means clustering methodology to categorize dataset E p r ( ω ) into n distinct groups. Each chronological day within the dataset was assigned a specific pattern label to facilitate analysis. The Euclidean distance metric served as the basis for assessing the similarity between patterns. In the process of de-duplication and revision, the d-dimensional vector space of data points E i ,   D = E i i = 1 , 2 , , t was partitioned by minimizing the objective function: F = ( arg min | | E i c j | | 2 ) . Here, c j represents the k-cluster centroids in the central set C = { c j j = 1 , 2 , , k } . This method helps to identify the optimal number of clusters for the dataset, ensuring the efficient and accurate partitioning of the data points [37]. The final result obtained contained several trend components; hence, the method is called tVMD.

2.2. Multi-Objective Wild Horse Optimization Algorithm

The living behavior of stallions, several mares, and foals in the wild horse society inspired the wild horse optimization algorithm (WHO) [38]. The decency behavior of horses, such as grazing, chasing, dominance, leadership, and mating, has been used to construct a new method. In comparison to other intelligence algorithms, the optimization process of horses departing from the group and mating with horses from other groups efficiently avoids the algorithm slipping into local optimization and balances the algorithm’s exploration and development. On this basis, this paper introduces the concepts of the Pareto optimal solution and external archiving mechanism into WHO and creatively proposes a multi-objective wild horse swarm optimization algorithm (MOWHO) to solve complex forecasting problems. Figure 1 illustrates the key link and the flowchart of the algorithm.
(1) 
Creating an initial population.
An initialized random population is constructed X = x 1 , x 2 , , x n and determines the objective value Ő i = Ő i 1 , Ő i 2 , , Ő i n . Then, the groups are divided. Let the population size be N ¯ ¯ ; then, the number of groups is Ğ = [ N ¯ ¯ * P S ] , where P S indicates the proportion of stallions. The remaining populations [ N ¯ ¯ Ğ ] are divided equally into these groups according to the number of stallions, S t a l l i o n (Leader). Figure 1a shows an example of this population division.
(2) 
Grazing behavior.
The algorithm positions the stallion Ğ as the central point of the grazing area, with group members moving and searching around the S t a l l i o n ’s perimeter by X i , Ğ j = 2 Ź * COS ( 2 π R Ź ) × ( S t a l l i o n j X i , G j ) + S t a l l i o n j , where X i , Ğ j is the position of the mare and foal and Ź is an adaptive mechanism: Ź = R 2 Θ Ї + R 2 Θ ( Ї ) . R i [ 2 , 2 ] is a uniform random number. The COS function enables the population to move in a 360-degree direction. The Ї index satisfies the condition ( P = = 0 ) and this results in P = R 1 < TDR   , Ї = P = = 0 and Ź = R 2 Θ Ї + R 2 Θ ( Ї ) , where TDR is the parameter regulator by TDR = 1 Iter × 1 / IMax , which gradually decreases to 0 with the number of iterations.
(3) 
Improved nonlinear update strategy.
To overcome the shortage of the only local search, a new nonlinear TDR ¯ ¯ n o n was proposed by TDR ¯ ¯ n o n = 2 × ( 1 ( Iter / IMax 1 / 3 ) 1 / 3 × ( Rand 1 ) , where Rand [ 0 , 1 ] is a random number. Through the nonlinear updating mode, the algorithm performs global and local searches in the whole iterative process. Therefore, TDR = TDR ¯ ¯ n o n in the paper.
(4) 
Horse mating behavior.
Horses have unique forms of decent behavior, i.e., the foal will leave the original population when it is about to mature, the male horse will join the single horse group, and the female horse will join another family group. Such decent behavior avoids inbreeding among horses. WHO uses the mean crossover operator to simulate this behavior in horses: X Ğ , k p ¯ = C r o s s o v e r ( X Ğ , i q ¯ , X Ğ , j z ¯ ) , where X Ğ , k p ¯ is the distance from group k; X Ğ , i q ¯ and X Ğ , j z ¯ are two maturing foals from i group and j group, respectively. Figure 1b shows this mating and departure behavior.
(5) 
Group leadership.
Each group has a leader ( S t a l l i o n Ğ i ) who leads the group to occupy the water hole Њ and occupy a dominant position. ( Њ is the current optimal solution).
S t a l l i o n Ğ i = X Ğ , i i f   C O S [ t ( X Ğ , i ) ] > C O S [ t ( S t a l l i o n Ğ i ) ] S t a l l i o n Ğ i i f   C O S [ t ( X Ğ , i ) ] < C O S [ t ( S t a l l i o n Ğ i ) ]
S t a l l i o n ¯ ¯ Ğ i = 2 Ź * COS ( 2 π R Ź ) × ( Њ S t a l l i o n Ğ i ) + Њ i f   R 3 > 0.5 2 Ź * COS ( 2 π R Ź ) × ( Њ S t a l l i o n Ğ i ) Њ i f   R 3 0.5
where S t a l l i o n ¯ ¯ Ğ i and S t a l l i o n Ğ i are the next and current positions of the leader of the i group. In addition, for the selection of leaders, S t a l l i o n ¯ ¯ Ğ i , in each group, we randomly select leaders S t a l l i o n ¯ ¯ Ğ 1 in the initial stage to maintain the diversity of algorithms. In the later stage, the position of the leader S t a l l i o n ¯ ¯ Ğ i , i = 2 , 3 , is exchanged in the group according to the Pareto dominance relationship.
(6) 
Adaptive T-distribution variation strategy.
The variation operator in the algorithm uses the T-distribution with the number of iterations (Iter) serving as the degree of freedom. This aims to balance exploration and development by S t a l l i o n ¯ ¯ Ğ i = S t a l l i o n ¯ ¯ Ğ i + T ( I t e r ) * S t a l l i o n ¯ ¯ Ğ i . The probability density function is as follows:
t n ( x ) = Γ n + 1 2 Γ n 2 n π 1 + x 2 n n + 1 2 , < x <
where Γ ( ) is the Gamma function.
(7) 
Multi-objective version introduction.
This paper innovatively introduces the Pareto optimal solution A * , and external Archiving mechanism Å r c into WHO, and constructs a multi-objective wild horse optimization algorithm (MOWHO) suitable for time-series forecasts. In the current iteration, Å r c serves primarily as a repository for the non-dominated solution A , thus establishing an upper limit Ň N . Upon acquiring a new non-dominated solution A n d via the iterative method, it undergoes comparison with the archived member { A 1 , A 2 , , A N } . If F ( A nd ) F ( A i ) , i [ 1 , N ] , it shall incorporate A n d in Å r c . To refine Å r c , the optimization algorithm will retain the superior solution within Å r c while discarding the least effective one. The selection of the most advantageous starting point for the subsequent iteration is conducted through the application of a roulette-based mechanism þ i r = ξ / Д ¯ i r , ξ > 1 , where ξ is considered a constant, while Д ¯ i r represents the number of solutions that are closely related to the i solution. þ i r assigns higher weights to solutions located in regions with fewer samples at the Pareto optimality boundary. This strategy aims to expand the coverage to provide a more diverse and comprehensive solution. It is worth noting that the two optimization objectives [ O b f 1 ¯ ¯ M , O b f 2 ¯ ¯ M ] and constraints of this study are as follows: O b f 1 ¯ ¯ M = Min M A E and O b f 2 ¯ ¯ M = Min s t d ( E P E P ) .

3. The Main Structure

The principal architecture is depicted in Figure 2, with the specific workflow delineated as follows. Initially, trend decomposition techniques are employed to extract trends, seasonal, and irregular components from raw electricity price data. Subsequently, the three most accurate forecasting models are selected for multiple simulations and parameter adjustments to ensure the stability of results. Ultimately, the multi-objective optimization algorithm is utilized to execute a weighted combination of the three models. And the comparative results are obtained by analyzing the error evaluation indexes.

4. Case Study and Evaluations

4.1. Datasets and Preprocessing

In this comparative analysis, we leveraged datasets derived from authentic electricity markets. The Singapore electricity market (SGP-EP), inaugurated in 2003, represents Asia’s pioneering deregulated market. Recent years have seen the Singaporean government adapt to the global movement towards energy transition, enacting long-term strategies for the decarbonization of its electricity sector. This initiative has culminated in the development of a sophisticated framework for electricity generation and retailing. Conversely, the Australian National Electricity Market (AUS-EP), initiated in 1998, extends over 5000 km, making it the world’s most extensive alternating current electricity system. Characterized by its leading position in global price volatility, it ranks among the most unstable electricity markets internationally.
In the context of electricity markets, the data on electricity prices exhibit pronounced volatility and are frequently characterized by a surplus of peak values. Utilizing unadjusted electricity price series as inputs for predictive models inherently compromises the precision of future price forecasts. To mitigate this challenge, our approach encompasses the identification and subsequent adjustment of aberrant outlier data points. For example, within the Singaporean electricity market, outliers are ascertained utilizing the 3σ criterion. Conversely, the Australian market’s outlier data, distinguished by their considerable deviation resulting in atypical standard deviations, necessitates outlier identification through box plot analysis. Following outlier adjustment, data imputation is conducted employing the piecewise cubic Hermite interpolation polynomial (PCHIP) method. It aims to approximate the original data by a series of cubic polynomials, each of which is responsible for interpolating between two points in the dataset. The polynomial is expressed as follows: H 3 ( x ) = 1 + 2 χ y 0 + x x 0 y 0 ψ 2 + 1 + 2 ψ y 1 + x x 1 y 1 χ , where x 0 , x 1 are the positions of the two neighboring points of the point to be interpolated, χ = ( x x 0 ) / ( x 1 x 0 ) , ψ = ( x x 1 ) / ( x 0 x 1 ) , y 0 , y 1 corresponds to the dependent variables of the independent variable x 0 , x 1 , and y 0 , y 1 are the corresponding derivatives. This methodology ensures the interpolation curve’s monotonicity and smoothness, effectively eliminating the introduction of supplementary local extrema. Such a feature renders it exceptionally applicable to contexts necessitating the mitigation of undue oscillations and the retention of inherent data trends and attributes.
Our dataset encompasses 31 days from 1 January to 31 January 2024 (sample interval: 15 min). We adopted a stratified split for each dataset into training, validation, and test sets in a 6:2:2 ratio, ensuring rigorous methodology and comprehensive analysis. The specific features of the datasets are shown in Table 2.

4.2. Evaluation Metrics

In assessing the predictive model, three essential metrics were utilized: the Mean Absolute Error (MAE), Mean Absolute Percentage Error (MAPE), and Root Mean Square Error (RMSE) [39,40]. The equation and indicator strengths are shown in Table 3.

4.3. Experimental Setup

In this segment, we test the proposed hybrid short-term electricity price forecast (HSEPF) model through experimental studies. We begin by outlining the model’s complex operational settings, followed by a thorough examination of the results obtained from these processes. The first case demonstrates the model’s significant ability to integrate various forecasting elements, showing a clear advantage over traditional methods. The second case highlights an improvement in performance compared to models based on different decomposition techniques. The third case emphasizes the model’s exceptional accuracy in predictions thanks to its novel approach of intelligently weighting multi-objectives.

4.3.1. Operational Settings

In the field of deep learning, model parameters play a pivotal role in determining both the computational efficiency and the accuracy of experimental simulations. This study has rigorously tested these parameters across a series of simulation experiments, with the definitive configurations enumerated in Table 4.

4.3.2. Case I: Comparison with Predictive Classical Models

Case I involved conducting a comparative experiment with a single model. In this scenario, ARIMA, BP, ELM, LSTM, and RF single models were first used for direct prediction. Subsequently, considering that combining multiple models would increase the complexity of the prediction framework, this paper intended to use the three best-performing single models for decomposition-ensemble [41]. The prediction results are shown in Table 5.
(a) For AUS-EP forecasting, the developed HSEFP system was evaluated alongside five benchmark models. This comparison provides a fundamental insight into the advanced capabilities of the proposed system. Specifically, when making a 1-step prediction—which involves forecasting the upcoming half-hour—the HSEFP system demonstrates superior performance with parameter settings of α tVMD = 2000, Κ tVMD = 5, C tVMD = 3, N P MOWHO = 100, and I t e r max MOWHO = 100. It achieved M A E one - step HSEPF = 7.3362, M A P E one - step HSEPF = 14.9100%, and R M S E one - step HSEPF = 9.7486. In contrast, the benchmarks reported MAE values above 10, MAPE values ranging from 22% to 27%, and RMSE values exceeding 15. The prediction accuracies are significantly better relative to all the base single models. Figure 3 shows the forecast trends and error indicators of each model.
(b) For SGP-EP forecasting, the overall forecast error was smaller because this market is more stable than the AUS-EP. When 1-step forecasting was performed, the best forecasting performance in the benchmark model was ELM, with parameters I E l m = 5, H E l m = 15, O E l m = 1, and T f Elm ¯ ¯ = sigmoid resulting in M A E one - step ELM = 5.9693, M A P E one - step ELM = 6.0769%, and R M S E one - step ELM = 9.3607. HSEFP has a better forecasting performance (with parameters N P MOWHO = 100, I t e r max MOWHO = 100, H E l m = 15, T f Elm ¯ ¯ = sigmoid, H L s t m = 150, I t e r max Lstm = 500, D e p t h max RF = 50, and L e a f max RF = 5), where the prediction is expressed as M A E one - step HSEPF = 4.2784, M A P E one - step HSEPF = 4.1918%, and R M S E one - step HSEPF = 5.9268. There is still a gap in the forecasting effectiveness of the single model.
(c) To compare the capability of the data decomposition–integration strategy, the three models ELM, LSTM, and RF, which performed better in the benchmark model, were trend decomposed, and each trend variable was predicted individually. As shown in Table 5, the prediction accuracy improved after integration. From the forecasting results for the AUS-EP, to ensure objectivity in the model comparisons, the ELM was unified with parameters I E l m = 5, H E l m = 15, O E l m = 1, and T f Elm ¯ ¯ = sigmoid; the changes in the indicators at 1-step prediction were M A E { ELM , tVMD - ELM } one - step = {10.8203→7.3713}, M A P E { ELM , tVMD - ELM } one - step = {22.6797%→15.3919%}, and R M S E { ELM , tVMD - ELM } one - step = {15.4311→9.9353}. The LSTM uniform parameters were H L s t m = 150, Γ t a u L s t m = 2, and I t e r max Lstm = 500; the changes in the indicators at one-step prediction were as follows: M A E { LSTM , tVMD - LSTM } one - step = {10.7900→7.5255}, M A P E { LSTM , tVMD - LSTM } one - step = {22.5088%→15.6030%}, and R M S E { LSTM , tVMD - LSTM } one - step = {15.4000→9.8997}. The RF model unification parameters were D e p t h max RF = 50 and L e a f max RF = 5; the changes in the indicators were M A E { RF , tVMD - RF } one - step = {11.1553→8.8530}, M A P E { RF , tVMD - RF } one - step = {24.9916%→18.8079%}, and R M S E { RF , tVMD - RF } one - step = {15.4319→11.8251}. From the performance of the metrics, the tVMD decomposition technique effectively improved the accuracy prediction. In the multi-step simulation results, it seemed that 88.89% of the indicators could reflect the superiority of the HSEPF system.
(d) The fitting effect of the SGP-EP models was satisfactory. Additionally, decomposing the model using tVMD significantly enhanced the performance. Under the parameters I E l m = 5, H E l m = 15, O E l m = 1, and T f Elm ¯ ¯ = sigmoid, there were M A E { ELM , tVMD - ELM } one - step = {5.9693→4.3112}, M A P E { ELM , tVMD - ELM } one - step = {6.0769%→4.2667%}, and R M S E { ELM , tVMD - ELM } one - step = {9.3607→5.9489}. The LSTM uniform parameters were I L s t m = 5, H L s t m = 150, O L s t m = 1, Γ t a u L s t m = 2, and I t e r max Lstm = 500; the results were M A E { LSTM , tVMD - LSTM } one - step = {6.0180→4.5092}, M A P E { LSTM , tVMD - LSTM } one - step = {6.0780%→4.4401%}, and R M S E { LSTM , tVMD - LSTM } one - step = {9.4193→6.2313}. The errors of the RF model were reduced (parameters were D e p t h max RF = 50 and L e a f min RF = 5) and the results were M A E { RF , tVMD - RF } one - step = {6.3716→5.4662}, M A P E { RF , tVMD - RF } one - step = {6.2949%→5.5410%}, and R M S E { RF , tVMD - RF } one - step = {9.4881→7.3617}. The prediction errors all decreased significantly, indicating that the introduction of the tVMD decomposition technique significantly enhanced the hybrid model.
Remark 1.
The effectiveness of the HSEPF system in predicting future electricity prices was validated through experimental results. This system surpassed other single models in forecasting accuracy, as evidenced by the results that support the validity and effectiveness of the proposed modeling theory and methodology. Additionally, the use of the tVMD decomposition technique significantly improved forecasting accuracy. Crucially, this method supports automatic adjustments based on incoming data, enhancing their applicability across various forecasting domains and its practical implementation.

4.3.3. Case ΙΙ: Comparison with Other Data Decomposition Methods

This experiment validated the ability of the tVMD method proposed in this paper to decompose trends in electricity price data, comparing it with the following three models: one without decomposition and two with other varieties of decomposition methods (EMD and EEMD models). The research results for the AUS-EP and SGP-EP markets are summarized in Table 6 and Figure 4, with specific comparative analyses as follows.
(a) Due to various uncertainties in the electricity market, electricity price series are typical examples of nonlinear and non-stationary data. Decomposition techniques can significantly improve predictive performance. For the AUS-EP, which is characterized by significant price fluctuations, the multi-step forecasting results of a combined deep learning and machine learning approach without data decomposition are M A E ELsR step = {10.7941, 13.8730, 14.6527}, M A P E ELsR step = {22.4224%, 31.1126%, 34.2981%}, and R M S E ELsR step = {15.3988, 18.7547, 19.3794}. And for the relatively stable SGP-EP, the combined model ELsR without data decomposition processing can obtain M A E ELsR step = {5.8931, 8.7458, 10.9286}, M A P E ELsR step = {5.8611%, 8.7855%, 11.0888%}, and R M S E ELsR step = {9.0542, 12.3424, 14.7562} in the calculation of error indicators for multi-step prediction. However, the use of decomposition techniques in the forecasting model reduces the MAPE, MAE, and RMSE values.
(b) For the 1-step forecasting analysis using two market datasets, the accuracy improvements were most pronounced following decomposition through the tVMD method. Notably, there was a significant reduction in all error metrics for the AUS-EP: M A P E AUS - EP H S E P F = 7.5124%, R M A E AUS - EP H S E P F = 5.6502 USD/MWh, and M A E AUS - EP H S E P F = 3.4579 USD/MWh. In the context of SGP-EP, forecasts derived from tVMD decomposition were most satisfactory, showing improved prediction precision with R M S E SGP - EP H S E P F = 3.12174 USD/MWh, M A P E SGP - EP H S E P F = 1.6693%, and M A E SGP - EP H S E P F = 1.6147 USD/MWh.
(c) In multi-step forecasting performances across two major market datasets, models decomposed via the tVMD method consistently maintained excellent predictive performance. There are significant improvements in the 2-step and 3-step forecasting accuracy in AUS-EP: M A P E AUS - EP { 2 , 3 } H S E P F = {10.1376%, 12.3727%}. The improvement effects of EEMD and EMD decomposed models are similar in the 2-step, but EEMD performs better at the 3-step, though it still lags behind the proposed HSEFP model by about 3% in terms of the MAPE index. For SGP-EP, the forecasting errors decreased by M A P E SGP - EP { 2 , 3 } H S E P F = {3.8889%, 4.5777%} for 2-step and 3-step forecasting, where the improvement effects of EEMD and EMD were closer.
Remark 2.
Before integration into predictive models, it is crucial to stabilize electricity price data and to adaptively isolate several components that exhibit periodicity. Comparative analyses robustly validate that forecasting precision enhancements afforded by tVMD outperform those achieved with alternative decomposition methodologies. Furthermore, the benefits of this decomposition approach intensify progressively with the increase in forecast steps. The potential of tVMD as a superior tool for improving the reliability and accuracy of forecasting models is underscored by this finding.

4.3.4. Case ΙΙΙ: Comparison with Other Combinatorial Weighting Methods

In this case, we compared the HSEPF model proposed in this paper with a combinatorial weighting strategy based on three different principles. Among them, simple average (SA) weightings were the most common method of model mixing, which reflects the average performance of multiple model combinations. The root mean square (RMS) weighting strategy amplifies large deviations by squaring them and emphasizes the volatility or difference in these values. The multi-objective Dragonfly optimization algorithm (MODA) is widely used in the fields of energy and financial forecasting. Compared with the model prediction results in Table 7, it was found that the HSEPF method achieves higher prediction accuracy and a smaller prediction error. Figure 5 shows the forecast trends and error indicators of each method for the SGP-EP market.
(a) In 1-step forecasting, the HSEPF technique outperformed the other weighting strategies in the forecasting performance in Singapore and Australian markets. tVMD-ELsR-SA, tVMD-ELsR-RMS, and tVMD-ELsR-MODA all performed well when applied to AUS-EP prediction, with tVMD-ELsR-SA being the strategy with the smaller forecasting errors, M A E AUS - EP SA = 7.7246, M A P E AUS - EP SA = 16.2421%, and R M S E AUS - EP SA = 10.4957. However, there was a gap relative to HSEPF, M A E AUS - EP HSEPF = 7.3362, M A P E AUS - EP HSEPF = 14.9100%, and R M S E AUS - EP HSEPF = 9.7486. On the other hand, among the comparison models in Singapore, tVMD-ELsR-RMS was the strategy with less prediction error with M A E SGP - RP RMS = 4.4409, M A P E SGP - RP RMS = 4.3885%, and R M S E SGP - RP RMS = 6.1485. Similarly, the prediction performance lost out to the proposed HSEPF system with M A E SGP - RP HSEPF = 4.2784, M A P E SGP - RP HSEPF = 4.1918%, and R M S E SGP - RP HSEPF = 5.9268.
(b) In the 2-step prediction of the AUS-EP, the prediction accuracy of the tVMD-ELsR-SA in MAPE index was better: M A P E AUS - EP SA = 20.8713%. At this time, with M A P E AUS - EP HSEPF = 20.9750%, M A E AUS - EP HSEPF = 9.8568, and R M S E AUS - EP HSEPF = 12.9024, the HSEPF model was slightly better than the tVMD-ELsR-SA model in the performance of MAE and MAPE prediction errors. In SGP-EP, the HSEPF model still maintained a good forecasting ability, where M A E SGP - RP HSEPF = 5.0703, M A P E SGP - RP HSEPF = 4.8966%, R M S E SGP - RP HSEPF = 6.7367. The overall performance was better than that of the best-performing comparative models: tVMD-ELsR-MODA with M A E SGP - RP MODA = 5.1215, M A P E SGP - RP MODA = 4.9873%, and R M S E SGP - RP MODA = 6.8570. Therefore, in general, the overall performance of the HSEPF model was better and can be widely used in different electricity price markets.
(c) For a 3-step prediction with a longer prediction period, the HSEPF model had excellent performance in model comparison. In the AUS-EP market, the forecast errors were M A E AUS - EP HSEPF = 10.2954, M A P E AUS - EP HSEPF = 21.9254%, and R M S E AUS - EP HSEPF = 13.2731. The tVMD-ELsR-SA, tVMD-ELsR-RMS, and tVMD-ELsR-MODA models had closer prediction accuracies, and the MAPE errors were all above 22%. The tVMD-ELsR-RMS model had a better RMSE in the SGP-EP market, which was lower than the HSEPF model with R M S E SGP - EP RMS = 0.1586 USD/MWh. However, the HSEPF model of the other two indicators still maintained low prediction errors as M A E SGP - EP HSEPF = 6.7639, M A P E SGP - EP HSEPF = 6.5111%.
Remark 3.
In this comprehensive study, we substantiated the robust theoretical foundation of the HSEPF model and its superiority in enhancing the predictive accuracy of complex systems. Comparative analyses revealed that other models failed to guarantee stable electricity price forecasts across varied market conditions. In contrast, the model proposed herein consistently demonstrates heightened precision in numerous forecasting scenarios, proficiently responding to the challenges presented by the significant fluctuations in electricity prices.

5. Discussion

To further demonstrate the accuracy of the proposed HSEPF, we perform two additional sets of model tests.

5.1. DM Test

From the comparison of error indicators, we can understand the model prediction performance but cannot discern whether the model is statistically significant or not. The Diebold–Mariano (DM) test can solve this problem, where DM applies to the comparison of predictive ability between two or more models [42]. We can assume that the true value is EP t ; t = 1 , , T , the two predicted values in the different models are E P 1 , t ; t = 1 , , T and E P 2 , t ; t = 1 , , T , respectively, and the prediction errors of the two models are i , t ; i = 1 , 2 , where i , t = E P i E P i , t i = 1 , 2 . We can assume that the loss associated with prediction i is a function of the prediction error, i , t , and denotes the loss function by Ψ ( i , t ) , i = 1 , 2 . δ t Δ = Ψ 1 , t Ψ 2 , t , where δ t Δ is a first-order difference function of the loss function, if and only if G ( δ t Δ ) = 0, E P 1 , t and E P 2 , t have the same prediction accuracy, so we construct the null hypothesis H 0 : G ( δ t Δ ) = 0   t and the alternative hypothesis H 1 : G ( δ t Δ ) 0   t . Here, H 0 means that the two prediction models have the same prediction accuracy, and H 1 states that the two do not have related prediction accuracies. The DM statistic can be constructed as T ( δ Δ ¯ μ ) δ Δ N ( 0 , 2 π f δ ( 0 ) ) , where δ Δ ¯ is the mean of δ t Δ , μ is the overall mean of the loss function, and f δ ( 0 ) is the spectral line density function when the loss difference is 0. The expression of f δ ( 0 ) is as follows: f δ ( 0 ) = + γ δ ( ϕ ) , where γ δ ( ϕ ) is the order self-covariance of δ t Δ , and the specific functional equation is γ δ ( ϕ ) = E [ ( δ t Δ μ ) ( δ t ϕ Δ μ ) ] , where D M = δ Δ ¯ / 2 π f δ Δ ( 0 ) / T ; N ( 0 , 1 ) , f δ Δ ( 0 ) and f δ ( 0 ) are the consistent estimators when the null hypothesis is valid; the DM statistic obeys N(0,1), HSEPF is used as the benchmark model for the DM test, and the other 14 control models are selected for the DM test. Table 8 presents the outcomes of the tests.
(a) In the comparison of single-mode models, whether for 1-step, 2-step, or 3-step predictions, the results uniformly achieved statistical significance at the 1% level. This outcome indicates a disparity in predictive capabilities between the proposed hybrid system and single models. Further analysis revealed that the predictive strength of models diminishes as the forward prediction steps increase. Moreover, the DM test values exhibited an increasing trend in 85% of the cases. This trend signifies that in 85% of instances, the decline in the predictive ability of single models is notably more pronounced than that of the HSEPF model, underscoring the superior multi-step forecasting stability of the hybrid model.
(b) In the comparison of machine learning models decomposed by tVMD, undecomposed deep learning integration models, and hybrid deep learning models based on other data decomposition techniques, it was found that there were cases where individual models performed better in prediction, but the prediction performance was not stable in multi-step prediction. Among them, the hybrid model ELsR without data decomposition had the highest DM values: D M ELsR AUS - EP = {7.1362, 8.5006, 8.1303} and D M ELsR SGP - EP = {3.5945, 6.4505, 6.8781}, which are both greater than z 0.01 / 2 = 2.576, suggesting that the prediction without data decomposition is the weakest. Further analysis revealed that the decomposition technique helped improve the prediction accuracy, whether it was a single model or a hybrid model. For example, in the 2-step prediction in AUS-EP, tVMD-ELM has D M t V M D E L M AUS - EP = 0.5748, and in the 3-step prediction model, tVMD-LSTM has a statistically non-significant difference with the model proposed in this paper, D M t V M D L S T M AUS - EP = 0.1518. It suggests that the forecasting performance of deep learning models, after undergoing trend decomposition through tVMD, closely matches that of the optimal model in specific instances.
(c) In the optimization comparison of different weighting strategies, the DM results in the 1-step prediction for the Australian and Singaporean markets, which both exceed z 0.01 / 2 = 2.576, indicating a significant difference between the HSEPF model and others at a 99% confidence level. This suggests the effectiveness of the HSEPF model. In the multi-step prediction performance, the majority of models have DM values exceeding z 0.1 / 2 = 1.645, signifying a significant difference between models at a 90% confidence level. The HSEPF model demonstrates a leading advantage in electricity price prediction. These findings underscore the robustness and superiority of the HSEPF model in forecasting electricity prices, positioning them as a valuable tool in the energy market.
Remark 4.
According to the results in Table 8, 90.48% of the models passed the test in the AUS-EP, while 88% of the models passed the test in the SGP-EP. This demonstrates a significant difference between the comparative models, confirming the predictive effectiveness of the HSEPF model.

5.2. Performance Improvement Analysis

In predictive research, the improvement percentage of predictive indicators serves as a crucial metric for assessing the predictive capabilities of models. This is quantified as p i m p r o v e m e n t ¯ ¯ = | ( M e t r i c compared M e t r i c HSEPF ) / M e t r i c compared | , where M e t r i c compared denotes the evaluation metric of the reference model, and M e t r i c HSEPF signifies the predictive indicator from the composite system developed in this study. p i m p r o v e m e n t ¯ ¯ quantifies the enhancement in predictive performance, with higher values reflecting the superior predictive accuracy of HSEFP. The comparative outcomes of models in AUS-EP and SGP-EP are summarized in Table 9.
(a) Upon analyzing the results of multi-step forecasting, it was evident that the HSEFP prediction system had significantly improved the accuracy and stability of electricity price forecasting in Australia. Particularly, for the single prediction model, the minimum values of the model average improvement rates were p L S T M Avg . MAE ¯ ¯ = 30.30%, p A R I M A Avg . MAPE ¯ ¯ = 32.78%, and p R F Avg . RMSE ¯ ¯ = 33.38%. However, when predicting using the single model after tVMD model decomposition, the model improvement rates were all below 10%, indicating the significant impact of tVMD technology on enhancing model prediction accuracy. The average improvement rates of models decomposed by EMD and EEMD are the minimal values of p E M D E L s R Avg . RMSE ¯ ¯ = 13.68% and p E E M D E L s R Avg . MAE ¯ ¯ = 15.93%, respectively, suggesting that these two decomposition techniques have weaker capabilities in decomposing electricity price sequences compared to the tVMD technology proposed in this study. Among the three different combination weighted models, only p t V M D E L s R M O D A Avg . MAPE ¯ ¯ exceeded 10%, indicating that the hybrid weighting strategy effectively enhanced model prediction capability, yet the HSEFP proposed in this study still outperformed the other comparative strategies.
(b) The electricity price sequence in Singapore exhibited relative stability, with small prediction errors. However, the improvement in forecast accuracy with the HSEFP model was also significant when considering the enhanced percentage. Compared to ARIMA, the HSEFP model demonstrates a superiority of p A R I M A Avg . RMSE ¯ ¯ = 51.10%. In the models decomposed by tVMD, the improvement rates were all above 10%, with p t V M D R F Avg . MAE ¯ ¯ = 18.33%, p t V M D R F Avg . MAPE ¯ ¯ = 22.91%, and p t V M D R F Avg . RMSE ¯ ¯ = 15.31%, showcasing the excellent predictive performance of the MOWHO-weighted strategy. This finding highlights the capability of tVMD to handle electricity price data, significantly enhancing the predictive accuracy of deep learning algorithms and positioning it as a crucial tool for advancing predictive analytics in the field.
Remark 5.
The improvement in percentage through model enhancement further confirms the predictive ability of the HSEFP system. It validates the significant effect of the tVMD technique in enhancing model prediction accuracy, and compared with other weighting strategies, the MOWHO algorithm has clear advantages, which can be further extended to the analysis of time-series prediction.

6. Conclusions and Future Work

The HSEPF system, introduced in this study, is a novel approach for forecasting short-term electricity prices. Initially, acknowledging the high volatility and irregular movements characteristic of electricity price sequences, the HSEPF system applied a stabilizing transformation to the raw data before forecasting. This involves the adaptive extraction of several feature components with distinct central frequencies. Specifically, this study employed the tVMD technique to isolate the trend, cyclical, and irregular components of the dataset. This approach significantly enhanced the efficiency of feature extraction in the forecasting model. Comparative analysis and discussion provided compelling evidence that the application of tVMD exerts a scalable impact on forecasting accuracy, warranting further exploration and validation in a broader data context.
Subsequently, the model inputs the decomposed sequence components into both linear and nonlinear machine learning models for separate predictions before integrating the results from each component. The final stage involved selecting the top three performing single models and employing the MOWHO model to obtain a weighted amalgamation of their outcomes. The experimental results demonstrated that the proposed HSEPF model was capable of predicting electricity prices across different market types. The forecast results in the Singapore market were as follows: M A E one - step HSEPF = 4.2784, M A P E one - step HSEPF = 4.1918%, and R M S E one - step HSEPF = 5.9268. Australian market data fluctuated dramatically but still obtained smaller errors as follows: M A E one - step HSEPF = 7.3362, M A P E one - step HSEPF = 14.9100%, and R M S E one - step HSEPF = 9.7486. This system not only enhances prediction accuracy through innovative decomposition and optimization techniques but also demonstrates its universality and effectiveness under different market conditions. This study highlights the potential of combining advanced analytical methods with deep learning models to improve predictive performance in the non-stationary, complex field of electricity price forecasting, providing valuable insights for academia and industry stakeholders.
In future research, two aspects can be explored. Firstly, a generalized and precise electricity price prediction model can be developed by considering various factors, such as user behavior, holidays, and weather. Secondly, future work can focus on peak electricity price prediction based on the current research to achieve higher prediction accuracy.

Author Contributions

Conceptualization, methodology, funding acquisition, supervision, writing—review and editing, H.L.; writing—original draft preparation, data curation, software, validation, visualization, Y.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the HSSF of the Chinese Ministry of Education (No. 23YJA790053) and the NSFC (No. 12071302).

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Acknowledgments

The authors extend their gratitude to the editors and reviewers for their invaluable comments and suggestions.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Jiang, H.; Dong, Y.; Wang, J. Electricity Price Forecasting Using Quantile Regression Averaging with Nonconvex Regularization. J. Forecast. 2024, 43, 1859–1879. [Google Scholar] [CrossRef]
  2. Indira, G.; Bhavani, M.; Brinda, R.; Zahira, R. Electricity Load Demand Prediction for Microgrid Energy Management System Using Hybrid Adaptive Barnacle-Mating Optimizer with Artificial Neural Network Algorithm. Energy Technol. 2024, 12, 2301091. [Google Scholar] [CrossRef]
  3. Iftikhar, H.; Turpo-Chaparro, J.E.; Canas Rodrigues, P.; López-Gonzales, J.L. Forecasting Day-Ahead Electricity Prices for the Italian Electricity Market Using a New Decomposition—Combination Technique. Energies 2023, 16, 6669. [Google Scholar] [CrossRef]
  4. Tan, Y.Q.; Shen, Y.X.; Yu, X.Y.; Lu, X. Day-Ahead Electricity Price Forecasting Employing a Novel Hybrid Frame of Deep Learning Methods: A Case Study in NSW, Australia. Electr. Power Syst. Res. 2023, 220, 109300. [Google Scholar] [CrossRef]
  5. Sai, W.; Pan, Z.; Liu, S.; Jiao, Z.; Zhong, Z.; Miao, B.; Chan, S.H. Event-Driven Forecasting of Wholesale Electricity Price and Frequency Regulation Price Using Machine Learning Algorithms. Appl. Energy 2023, 352, 121989. [Google Scholar] [CrossRef]
  6. Zhang, T.; Tang, Z.; Wu, J.; Du, X.; Chen, K. Short Term Electricity Price Forecasting Using a New Hybrid Model Based on Two-Layer Decomposition Technique and Ensemble Learning. Electr. Power Syst. Res. 2022, 205, 107762. [Google Scholar] [CrossRef]
  7. Liu, D.; Wei, Y.; Yang, S.; Guan, Z. Electricity Price Forecast Using Combined Models with Adaptive Weights Selected and Errors Calibrated by Hidden Markov Model. Math. Probl. Eng. 2013, 2013, 648101. [Google Scholar] [CrossRef]
  8. Girish, G.P. Spot Electricity Price Forecasting in Indian Electricity Market Using Autoregressive-GARCH Models. Energy Strategy Rev. 2016, 11–12, 52–57. [Google Scholar] [CrossRef]
  9. Zhao, Z.; Wang, C.; Nokleby, M.; Miller, C.J. Improving Short-Term Electricity Price Forecasting Using Day-Ahead LMP with ARIMA Models. In Proceedings of the 2017 IEEE Power & Energy Society General Meeting, Chicago, IL, USA, 16–20 July 2017; pp. 1–5. [Google Scholar]
  10. Imani, M.H.; Bompard, E.; Colella, P.; Huang, T. Forecasting Electricity Price in Different Time Horizons: An Application to the Italian Electricity Market. IEEE Trans. Ind. Appl. 2021, 57, 5726–5736. [Google Scholar] [CrossRef]
  11. Singh, A.; Sahay, K.B. Short-Term Demand Forecasting by Using ANN Algorithms. In Proceedings of the 2018 International Electrical Engineering Congress, Krabi, Thailand, 7–9 March 2018; pp. 1–4. [Google Scholar] [CrossRef]
  12. Wu, S.; He, L.; Zhang, Z.; Du, Y. Forecast of Short-Term Electricity Price Based on Data Analysis. Math. Probl. Eng. 2021, 2021, 6637183. [Google Scholar] [CrossRef]
  13. Wang, P.; Xu, K.; Ding, Z.; Du, Y.; Liu, W.; Sun, B.; Zhu, Z.; Tang, H. An Online Electricity Market Price Forecasting Method Via Random Forest. IEEE Trans. Ind. Appl. 2022, 58, 7013–7021. [Google Scholar] [CrossRef]
  14. Huang, S.; Shi, J.; Wang, B.; Lyu, J.; Nie, N.; Dong, X.; Li, H.; Zhang, S.; Ren, X. A Meta-Learning Based Method for Day Ahead Electricity Price Forecasting in Markets with Renewable Energy Resources. In Proceedings of the 2023 IEEE 12th Global Conference on Consumer Electronics, Nara, Japan, 10–13 October 2023; pp. 288–290. [Google Scholar] [CrossRef]
  15. Rezaei, N.; Rajabi, R.; Estebsari, A. Electricity Price Forecasting Model Based on Gated Recurrent Units. In Proceedings of the 2022 IEEE International Conference on Environment and Electrical Engineering and 2022 IEEE Industrial and Commercial Power Systems Europe, Prague, Czech Republic, 28 June–1 July 2022; pp. 1–5. [Google Scholar] [CrossRef]
  16. Wang, B.; Wei, W.; Su, W. Short-Term Electricity Price Forecasting Based on Data Mining. In Proceedings of the 2022 2nd International Conference on Algorithms, High Performance Computing and Artificial Intelligence, Guangzhou, China, 21–23 October 2022; pp. 394–397. [Google Scholar] [CrossRef]
  17. Yorat, E.; Ozbek, N.S.; Zor, K.; Saribulut, L. Day-Ahead Electricity Price Forecasting Using Artificial Intelligence-Based Algorithms. In Proceedings of the 2023 International Conference on Innovation and Intelligence for Informatics, Computing, and Technologies, Sakheer, Bahrain, 20–21 November 2023; pp. 121–126. [Google Scholar] [CrossRef]
  18. Xu, H.; Hu, F.; Liang, X.; Gunmi, M.A. Attention Mechanism Multi-Size Depthwise Convolutional Long Short-Term Memory Neural Network for Forecasting Real-Time Electricity Prices. IEEE Trans. Power Syst. 2024, 39, 6277–6289. [Google Scholar] [CrossRef]
  19. Yang, W.; Sun, S.; Hao, Y.; Wang, S. A Novel Machine Learning-Based Electricity Price Forecasting Model Based on Optimal Model Selection Strategy. Energy 2022, 238, 121989. [Google Scholar] [CrossRef]
  20. Bozlak, Ç.B.; Yaşar, C.F. An Optimized Deep Learning Approach for Forecasting Day-Ahead Electricity Prices. Electr. Power Syst. Res. 2024, 229, 110129. [Google Scholar] [CrossRef]
  21. Najafi, A.; Homaee, O.; Jasinski, M.; Golshan, M.; Leonowicz, Z. Application of Extreme Learning Machine-Autoencoder to Medium Term Electricity Price Forecasting. In Proceedings of the 2022 IEEE International Conference on Environment and Electrical Engineering and 2022 IEEE Industrial and Commercial Power Systems Europe, Prague, Czech Republic, 28 June–1 July 2022; pp. 1–4. [Google Scholar] [CrossRef]
  22. Qu, Z.; He, L.; Ge, X.; Wang, F.; Xu, F.; Lu, J. A Two-Stage Forecasting Approach for Day-Ahead Electricity Price Based on Improved Wavelet Neural Network with ELM Initialization. IEEE Trans. Ind. Appl. 2024, 60, 5061–5073. [Google Scholar] [CrossRef]
  23. Wang, J.; Liu, F.; Song, Y.; Zhao, J. A Novel Model: Dynamic Choice Artificial Neural Network (DCANN) for an Electricity Price Forecasting System. Appl. Soft Comput. J. 2016, 48, 281–297. [Google Scholar] [CrossRef]
  24. Zhang, H.; Hu, W.; Member, S.; Cao, D.; Member, S.; Huang, Q.; Chen, Z.; Blaabjerg, F. A Temporal Convolutional Network Based Hybrid Model of Short-Term Electricity Price Forecasting. CSEE J. Power Energy Syst. 2021, 10, 1119–1130. [Google Scholar] [CrossRef]
  25. Jiang, P.; Nie, Y.; Wang, J.; Huang, X. Multivariable Short-Term Electricity Price Forecasting Using Artificial Intelligence and Multi-Input Multi-Output Scheme. Energy Econ. 2023, 117, 106471. [Google Scholar] [CrossRef]
  26. Wang, J.; Yang, W.; Du, P.; Niu, T. Outlier-Robust Hybrid Electricity Price Forecasting Model for Electricity Market Management. J. Clean. Prod. 2020, 249, 119318. [Google Scholar] [CrossRef]
  27. Zhang, J.; Tan, Z.; Wei, Y. An Adaptive Hybrid Model for Short Term Electricity Price Forecasting. Appl. Energy 2020, 258, 114087. [Google Scholar] [CrossRef]
  28. Huang, C.J.; Shen, Y.; Chen, Y.H.; Chen, H.C. A Novel Hybrid Deep Neural Network Model for Short-Term Electricity Price Forecasting. Int. J. Energy Res. 2021, 45, 2511–2532. [Google Scholar] [CrossRef]
  29. Yang, Z.; Ce, L.; Lian, L. Electricity Price Forecasting by a Hybrid Model, Combining Wavelet Transform, ARMA and Kernel-Based Extreme Learning Machine Methods. Appl. Energy 2017, 190, 291–305. [Google Scholar] [CrossRef]
  30. Zhang, X.; Wang, J.; Gao, Y. A Hybrid Short-Term Electricity Price Forecasting Framework: Cuckoo Search-Based Feature Selection with Singular Spectrum Analysis and SVM. Energy Econ. 2019, 81, 899–913. [Google Scholar] [CrossRef]
  31. Xu, Y.; Li, J.; Wang, H.; Du, P. A Novel Probabilistic Forecasting System Based on Quantile Combination in Electricity Price. Comput. Ind. Eng. 2024, 187, 109834. [Google Scholar] [CrossRef]
  32. Yang, W.; Wang, J.; Niu, T.; Du, P. A Novel System for Multi-Step Electricity Price Forecasting for Electricity Market Management. Appl. Soft Comput. J. 2020, 88, 106029. [Google Scholar] [CrossRef]
  33. Kavakci, G.; Cicekdag, B.; Ertekin, S. Time Series Prediction of Solar Power Generation Using Trend Decomposition. Energy Technol. 2024, 12, 2300914. [Google Scholar] [CrossRef]
  34. Dragomiretskiy, K.; Zosso, D. Variational Mode Decomposition. IEEE Trans. Signal Process. 2014, 62, 531–544. [Google Scholar] [CrossRef]
  35. Zhang, C.; Fu, Y.; Gong, L. Short-Term Electricity Price Forecast Using Frequency Analysis and Price Spikes Oversampling. IEEE Trans. Power Syst. 2023, 38, 4739–4751. [Google Scholar] [CrossRef]
  36. Rockafellar, R.T. A Dual Approach to Solving Nonlinear Programming Problems by Unconstrained Optimization. Math. Program. 1973, 5, 354–373. [Google Scholar] [CrossRef]
  37. Shi, C.; Wei, B.; Wei, S.; Wang, W.; Liu, H.; Liu, J. A Quantitative Discriminant Method of Elbow Point for the Optimal Number of Clusters in Clustering Algorithm. EURASIP J. Wirel. Commun. Netw. 2021, 2021, 31. [Google Scholar] [CrossRef]
  38. Naruei, I.; Keynia, F. Wild Horse Optimizer: A New Meta-Heuristic Algorithm for Solving Engineering Optimization Problems; Springer: London, UK, 2022; Volume 38, ISBN 0123456789. [Google Scholar]
  39. Wang, X.; An, Y.; Wang, J.; Shao, Y. A Novel Prediction-Integration Forecasting System for Short Wind Speed Based on Combined Data Preprocessing Technique and Weight Determination Strategy. Energy Technol. 2024, 12, 2300889. [Google Scholar] [CrossRef]
  40. Zhang, W.; Kou, M.; Lv, M.; Shao, Y. Improved Combined System and Application to Precipitation Forecasting Model. Alex. Eng. J. 2022, 61, 12739–12757. [Google Scholar] [CrossRef]
  41. Wang, J.; Guo, H.; Song, A. Photovoltaic Power Combination Prediction System Based on Improved Multi-Objective Optimization Algorithm and Nonlinear Weighting Strategy. Expert Syst. 2023, 40, e13209. [Google Scholar] [CrossRef]
  42. Diebold, F.X.; Mariano, R.S. Comparing Predictive Accuracy. J. Bus. Econ. Stat. 1995, 13, 253–263. [Google Scholar] [CrossRef]
Figure 1. Key link and the flowchart of MOWHO.
Figure 1. Key link and the flowchart of MOWHO.
Energies 17 04833 g001
Figure 2. The main structure of this study.
Figure 2. The main structure of this study.
Energies 17 04833 g002
Figure 3. The forecasting performances in Case I for the AUS market.
Figure 3. The forecasting performances in Case I for the AUS market.
Energies 17 04833 g003
Figure 4. The forecasting performances in Case ΙΙ.
Figure 4. The forecasting performances in Case ΙΙ.
Energies 17 04833 g004
Figure 5. The forecasting performances in Case ΙΙΙ for the SGP market.
Figure 5. The forecasting performances in Case ΙΙΙ for the SGP market.
Energies 17 04833 g005
Table 1. A comprehensive overview based on the representative literature.
Table 1. A comprehensive overview based on the representative literature.
CategorizationRelated StudiesMain TechnologyDiscovery
Statistical analysis-based categorizationLiu, Wei, Yang, & Guan, 2013 [7]HMMUsing the firing matrix of the HMM to correct the residuals of the combined prediction.
Girish, 2016 [8]Autoregressive-GARCHA study on modeling and forecasting prices in a market with demand exceeding supply.
Zhao, Wang, Nokleby, & Miller, 2017 [9]ARIMAIncorporating an exogenous time series into the model enhances precision.
Machine learning modelsSingh & Sahay, 2018 [11]ANNAccuracy can be maintained even in the presence of high volatility in the raw data.
Wu, He, Zhang, & Du, 2021 [12]DE-SVMThe differential evolutionary (DE) enhanced the precision of the SVM model while concurrently diminishing the computational time required.
P. Wang et al., 2022 [13]RFAn adaptive online prediction methodology that employed RF algorithms was designed to dynamically adjust the sizes of its training datasets.
Deep learning modelsRezaei, Rajabi, & Estebsari, 2022 [15]GRUThis approach resulted in stable and robust forecasting performance.
B. Wang, Wei, & Su, 2022 [16]LSTMLSTM networks led to the selection of more accurate forecast samples.
Yorat, Ozbek, Zor, & Saribulut, 2023 [17].XGBoostThe test results indicated that the XGBoost model outperformed the others in terms of error metrics.
Hybrid modelsQu et al., 2024 [22]ELM-WNN (IWNN)A two-stage forecasting algorithm that integrates pattern recognition and classification-based forecasting was proposed.
P. Jiang, Nie, Wang, & Huang, 2023 [25]ICEEMDAN-MSSAAnchoring electricity price forecasts to a loading series marked by minimal variability.
Y. Xu, Li, Wang, & Du, 2024 [31]IMOTSOThis approach yields probability density estimation curves that closely align with the actual values, demonstrating minimal deviation.
Table 2. Statistical features of the datasets.
Table 2. Statistical features of the datasets.
DatasetSamplesNumbersStatistical Indicator (USD/MWh)
MaxMinMedianMEANStd.
AUS-EPAll samples1488127.834.1559.9861.5425.03
Training set912127.834.1558.0060.1026.53
Validation set288125.7210.0061.0265.9321.83
Testing set288127.819.5460.7761.6922.55
SGP-EPAll samples1488183.6247.48112.17109.2325.75
Training set912183.5429.08107.68101.3833.22
Validation set288177.1749.43116.27107.5526.19
Testing set288183.6229.08113.12107.3827.58
Table 3. Three indicators.
Table 3. Three indicators.
EquationStrengths
M A E = 1 n t = 1 n E P t E P t The average difference between predicted and actual values is directly provided, relatively unaffected by extreme values.
M A P E = 1 n t = 1 n E P t E P t E P t * 100 % Expression of the error as a percentage allows for fair comparisons between datasets of different scales.
R M S E = 1 n t = 1 n E P t E P t 2 As the error is squared and then averaged, a larger penalty is given to large prediction errors.
Note:  E P t and E P t represent actual and predicted prices, t represents the sample index in chronological order, and t = 1, 2,…, n, and n denotes the total number of samples.
Table 4. Parameter settings.
Table 4. Parameter settings.
ModelParameterSymbolValueReason
ELMInput-layer nodes I E l m 5Trial and error
Hidden-layer nodes H E l m 15Trial and error
Output-layer nodes O E l m 1Target outcome
Transfer function T f Elm ¯ ¯ sigmoidCommon preset
LSTMInput-layer nodes I L s t m 5Trial and error
Hidden-layer nodes H L s t m 150Trial and error
Output-layer nodes O L s t m 1Target outcome
Gradient threshold Γ t a u L s t m 2Common preset
Maximum iteration I t e r max Lstm 500Trial and error
RFMax_depth D e p t h max RF 50Trial and error
Min_samples_leaf L e a f max RF 5Trial and error
tVMDModerate bandwidth constraint α tVMD 2000Trial and error
Modes K Κ tVMD 5Trial and error
categories C tVMD 3outcome
MOWHOPopulation number N P MOWHO 100Trial and error
Maximum iteration I t e r max MOWHO 100Trial and error
MODAPopulation number N P MODA 100Parameter consistency
Maximum iteration I t e r max MODA 100Parameter consistency
Table 5. Results of the compared models for different markets.
Table 5. Results of the compared models for different markets.
Electricity
Markets
Model1-Step2-Step3-Step
MAE
(USD/MWh)
MAPE
(%)
RMSE
(USD/MWh)
MAE
(USD/MWh)
MAPE
(%)
RMSE
(USD/MWh)
MAE
(USD/MWh)
MAPE
(%)
RMSE
(USD/MWh)
AUS-EPARIMA11.346021.887716.050315.186930.116920.313017.061434.312922.6320
BP11.628226.098315.867314.869232.616719.617115.667737.422920.5253
ELM10.820322.679715.431113.887330.937018.947214.812934.431919.7708
LSTM10.790022.508815.400013.908930.911318.924214.657834.282919.5256
RF11.155324.991615.431914.408532.761318.738515.215236.823019.5678
tVMD-ELM7.371315.39199.93539.914720.875312.935610.490822.434313.5512
tVMD-LSTM7.525515.603010.22629.899720.765213.077010.295822.046013.2776
tVMD-RF8.853018.807911.825110.154321.782913.306310.718823.179613.5515
Proposed HSEPF7.336214.91009.74869.856820.975012.902410.295421.925413.2731
SGP-EPARIMA6.57996.609210.370310.591610.651315.531213.982214.291719.7375
BP6.22456.39449.55379.36449.549913.075311.435911.954215.4864
ELM5.96936.07699.36078.94509.055812.748211.191011.550615.3115
LSTM6.01806.07809.41939.11829.260712.960911.193711.650215.4007
RF6.37166.29499.48819.33939.426912.613510.931711.461914.7077
tVMD-ELM4.31124.26675.94895.27125.18736.93896.87676.84019.0130
tVMD-LSTM4.50924.44016.23135.12335.02266.83306.79556.71798.8900
tVMD-RF5.46625.54107.36176.32406.49098.24567.81418.12129.9183
Proposed HSEPF4.27844.19185.92685.07034.89666.73676.76396.51119.1110
Note: This table shows the prediction results for five single models, three models based on the tVMD trend decomposition afterward and the proposed HSEPF, where M A E = 1 n t = 1 n E P t E P t , M A P E = 1 n t = 1 n ( | E P t E P t | ) / ( E P t ) * 100 % , and R M S E = 1 n t = 1 n E P t E P t 2 . Bold numbers indicate that the model performed better in all comparison results.
Table 6. Results of the compared models for data decomposition methods.
Table 6. Results of the compared models for data decomposition methods.
Electricity
Markets
Model1-Step2-Step3-Step
MAE
(USD/MWh)
MAPE
(%)
RMSE
(USD/MWh)
MAE
(USD/MWh)
MAPE
(%)
RMSE
(USD/MWh)
MAE
(USD/MWh)
MAPE
(%)
RMSE
(USD/MWh)
AUS-EPELsR10.794122.422415.398813.873031.112618.754714.652734.298119.3794
EMD-ELsR8.143217.392810.225710.850423.284714.128613.925830.545618.3599
EEMD-ELsR9.689420.206813.269611.220923.732215.364011.614924.757515.7821
Proposed HSEPF7.336214.91009.74869.856820.975012.902410.295421.925413.2731
SGP-EPELsR5.89315.86119.05428.74588.785512.342410.928611.088814.7562
EMD-ELsR4.90544.91677.30405.58055.64817.56678.47508.785211.2588
EEMD-ELsR4.37534.57516.42195.80975.87728.66738.19028.171811.7028
Proposed HSEPF4.27844.19185.92685.07034.89666.73676.76396.51119.1110
Note: ELsR represents the combination model of ELM, LSTM, and RF, which is modeled without any decomposition processing. This table shows the prediction results for different decompositions afterward, where M A E = 1 n t = 1 n E P t E P t , M A P E = 1 n t = 1 n ( | E P t E P t | ) / ( E P t ) * 100 % , and R M S E = 1 n t = 1 n E P t E P t 2 . Bold numbers indicate that the model performed better in all comparison results.
Table 7. Results of the compared models for other combinatorial weighting methods.
Table 7. Results of the compared models for other combinatorial weighting methods.
Electricity
Markets
Model1-Step2-Step3-Step
MAE
(USD/MWh)
MAPE
(%)
RMSE
(USD/MWh)
MAE
(USD/MWh)
MAPE
(%)
RMSE
(USD/MWh)
MAE
(USD/MWh)
MAPE
(%)
RMSE
(USD/MWh)
AUS-EPtVMD-ELsR-SA7.724616.242110.49579.859520.871312.984010.421822.429013.3598
tVMD-ELsR-RMS7.731616.259610.50439.859520.873212.984310.423222.433213.3602
tVMD-ELsR-MODA8.248117.793610.762810.364323.216413.436710.490622.899813.3782
Proposed HSEPF7.336214.91009.74869.856820.975012.902410.295421.925413.2731
SGP-EPtVMD-ELsR-SA4.61014.57746.38665.40945.38097.17796.93547.00199.0398
tVMD-ELsR-RMS4.44094.38856.14855.24885.18297.00806.85466.85288.9524
tVMD-ELsR-MODA4.53254.43056.32155.12154.98736.85706.86316.69238.9942
Proposed HSEPF4.27844.19185.92685.07034.89666.73676.76396.51119.1110
Note: ELsR represents the combination model of ELM, LSTM, and RF, which is modeled without any decomposition processing; SA is the simple average method; and RMS is the root mean square. M A E = 1 n t = 1 n E P t E P t , M A P E = 1 n t = 1 n ( | E P t E P t | ) / ( E P t ) * 100 % , and R M S E = 1 n t = 1 n E P t E P t 2 . Bold numbers indicate that the model performed better in all comparison results.
Table 8. DM test results.
Table 8. DM test results.
Electricity MarketAUS-EPSGP-EP
Model1-Step2-Step3-Step1-Step2-Step3-Step
ARIMA6.8043 ***7.4900 ***8.0509 ***4.8452 ***8.4343 ***8.8529 ***
BP7.1934 ***9.0850 ***8.8016 ***3.9281 ***7.5013 ***8.2048 ***
ELM7.1672 ***8.3486 ***8.6239 ***3.8637 ***7.2651 ***8.1112 ***
LSTM7.1313 ***8.3954 ***8.2850 ***3.9831 ***7.3610 ***8.2098 ***
RF7.1672 ***8.7916 ***8.3099 ***4.2884 ***6.2514 ***6.4664 ***
tVMD-ELM1.8208 *0.57482.3776 **0.42412.1572 **0.7557
tVMD-LSTM3.5418 ***1.8477 *0.15183.8626 ***1.3480 *1.8091 *
tVMD-RF5.6180 ***2.3662 **1.8415 *6.1373 ***5.3650 ***3.1619 ***
ELsR7.1362 ***8.5006 ***8.1303 ***3.5945 ***6.4505 ***6.8781 ***
EMD-ELsR0.99412.0821 **7.2362 ***3.0622 ***2.5157 **6.1797 ***
EEMD-ELsR5.8997 ***3.9738 ***3.5646 ***1.11192.4545 **3.2436 ***
tVMD-ELsR-SA3.7340 ***1.5210 *1.5346 *4.1060 ***3.0546 ***0.5165
tVMD-ELsR-RMS3.7501 ***1.5189 *1.5313 *3.0160 ***2.4181 **1.2954 *
tVMD-ELsR-MODA3.3949 ***2.3408 **1.11994.5819 ***1.6640 *0.9010
Note: ***, **, and * represent the significance levels under 99%, 95%, and 90% confidence level, respectively. D M = δ Δ ¯ / 2 π f δ Δ ( 0 ) / T ; N ( 0 , 1 ) . Here, z 0.01 / 2 = 2.576, z 0.05 / 2 = 1.96, and z 0.1 / 2 = 1.645.
Table 9. Improvement ratio results (%).
Table 9. Improvement ratio results (%).
Electricity
Markets
Model1-Step2-Step3-StepAvg.
MAEMAPERMSEMAEMAPERMSEMAEMAPERMSEMAEMAPERMSE
AUS-EPARIMA35.3431.8839.2635.1030.3536.4839.6636.1041.3536.7032.7839.03
BP36.9142.8738.5633.7135.6934.2334.2941.4135.3334.9739.9936.04
ELM32.2034.2636.8229.0232.2031.9030.5036.3232.8730.5734.2633.86
LSTM32.0133.7636.7029.1332.1431.8229.7636.0532.0230.3033.9833.51
RF34.2440.3436.8331.5935.9831.1432.3340.4632.1732.7238.9233.38
tVMD-ELM0.483.131.880.580.480.261.862.272.050.971.961.40
tVMD-LSTM2.524.444.670.431.011.340.000.550.030.982.002.01
tVMD-RF17.1320.7217.562.933.713.043.955.412.058.009.957.55
ELsR32.0433.5036.6928.9532.5831.2029.7436.0731.5130.2434.0533.14
EMD-ELsR9.9114.274.679.169.928.6826.0728.2227.7115.0517.4713.68
EEMD- ELsR24.2926.2126.5312.1611.6216.0211.3611.4415.9015.9316.4219.48
tVMD-ELsR-SA5.038.207.120.030.500.631.212.250.652.093.652.80
tVMD-ELsR-RMS5.118.307.190.030.490.631.232.260.652.123.682.83
tVMD-ELsR-MODA11.0616.219.424.909.653.981.864.260.795.9410.044.73
SGP-EPARIMA34.9836.5842.8552.1354.0356.6251.6254.4453.8446.2448.3551.10
BP31.2734.4537.9645.8648.7348.4840.8545.5341.1739.3242.9042.54
ELM28.3331.0236.6843.3245.9347.1639.5643.6340.5037.0740.1941.45
LSTM28.9131.0337.0844.3947.1248.0239.5744.1140.8437.6240.7641.98
RF32.8533.4137.5345.7148.0646.5938.1343.1938.0538.9041.5540.73
tVMD-ELM0.761.760.373.815.602.911.644.811.092.074.061.46
tVMD-LSTM5.125.594.891.032.511.410.473.082.492.213.732.93
tVMD-RF21.7324.3519.4919.8224.5618.3013.4419.838.1418.3322.9115.31
ELsR27.4028.4834.5442.0344.2645.4238.1141.2838.2635.8438.0139.41
EMD-ELsR12.7814.7418.869.1413.3110.9720.1925.8919.0814.0417.9816.30
EEMD- ELsR2.218.387.7112.7316.6822.2717.4120.3222.1510.7915.1317.38
tVMD-ELsR-SA7.208.427.206.279.006.152.477.010.795.318.144.71
tVMD-ELsR-RMS3.664.483.613.405.523.871.324.991.772.795.003.08
tVMD-ELsR-MODA5.615.396.241.001.821.751.452.711.302.683.303.10
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Luo, H.; Shao, Y. Advanced Optimal System for Electricity Price Forecasting Based on Hybrid Techniques. Energies 2024, 17, 4833. https://doi.org/10.3390/en17194833

AMA Style

Luo H, Shao Y. Advanced Optimal System for Electricity Price Forecasting Based on Hybrid Techniques. Energies. 2024; 17(19):4833. https://doi.org/10.3390/en17194833

Chicago/Turabian Style

Luo, Hua, and Yuanyuan Shao. 2024. "Advanced Optimal System for Electricity Price Forecasting Based on Hybrid Techniques" Energies 17, no. 19: 4833. https://doi.org/10.3390/en17194833

APA Style

Luo, H., & Shao, Y. (2024). Advanced Optimal System for Electricity Price Forecasting Based on Hybrid Techniques. Energies, 17(19), 4833. https://doi.org/10.3390/en17194833

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop