<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Energy Res.</journal-id>
<journal-title>Frontiers in Energy Research</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Energy Res.</abbrev-journal-title>
<issn pub-type="epub">2296-598X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1252067</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2023.1252067</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>An ensemble model for short-term wind power prediction based on EEMD-GRU-MC</article-title>
<alt-title alt-title-type="left-running-head">Wang et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenrg.2023.1252067">10.3389/fenrg.2023.1252067</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Peilin</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Su</surname>
<given-names>Chengguo</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2103421/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Li</surname>
<given-names>Li</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Yuan</surname>
<given-names>Wenlin</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Guo</surname>
<given-names>Chaoyu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>School of Water Conservancy and Transportation</institution>, <institution>Zhengzhou University</institution>, <addr-line>Zhengzhou</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Yellow River Laboratory</institution>, <institution>Zhengzhou University</institution>, <addr-line>Zhengzhou</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1879202/overview">Pedro Haro</ext-link>, Universidad de Sevilla, Spain</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2058200/overview">Bo Ming</ext-link>, Xi&#x2019;an University of Technology, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1525708/overview">Ling-Ling Li</ext-link>, Hebei University of Technology, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Chengguo Su, <email>suchguo@163.com</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>23</day>
<month>01</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>11</volume>
<elocation-id>1252067</elocation-id>
<history>
<date date-type="received">
<day>03</day>
<month>07</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>27</day>
<month>12</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Wang, Su, Li, Yuan and Guo.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Wang, Su, Li, Yuan and Guo</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>As a kind of clean and renewable energy, wind power is of great significance for alleviating energy crisis and environmental pollution. However, the strong randomness and large volatility of wind power bring great challenges to the dispatching and safe operation of the power grid. Hence, accurate and reliable short-term prediction of wind power is crucial for the power grid dispatching department arranging reasonable day-ahead generation schedules. Targeting the problem of low model prediction accuracy caused by the strong intermittency and large volatility of wind power, this paper develops a novel ensemble model for short-term wind power prediction which integrates the ensemble empirical mode decomposition (EEMD) algorithm, the gated recurrent unit (GRU) model and the Markov chain (MC) technique. Firstly, the EEMD algorithm is used to decompose the historical wind power sequence into a group of relatively stationary subsequences to reduce the influence of random fluctuation components and noise. Then, the GRU model is employed to predict each subsequence, and the predicted values of each subsequence are aggregated to get the preliminary prediction results. Finally, to further enhance the prediction accuracy, the MC is used to modified the prediction results. A large number of numerical examples indicates that the proposed EEMD-GRU-MC model outperforms the six benchmark models (i.e., LSTM, GRU, EMD-LSTM, EMD-GRU, EEMD-LSTM and EEMD-GRU) in terms of multiple evaluation indicators. Taking the spring dataset of the ZMS wind farm, for example, the MAE, RMSE and MAPE of the EEMD-GRU-MC model is 1.37 MW, 1.97 MW, and from 1.76%, respectively. Moreover, the mean prediction error of the developed model in all scenarios is less than or close to 2%. After 30 iterations, the proposed model uses an average of about 35&#xa0;min to accurately predict the wind power of the next day, proving its high computation efficiency. It can be concluded that the ensemble model based on EEMD-GRU-MC is a promising prospect for short-term wind power prediction.</p>
</abstract>
<kwd-group>
<kwd>ensemble empirical mode decomposition</kwd>
<kwd>gated recurrent unit</kwd>
<kwd>Markov chain</kwd>
<kwd>short-term wind power prediction</kwd>
<kwd>ensemble forecasting models</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Process and Energy Systems Engineering</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>In order to cope with the global energy crisis and climate change, renewable energy has become the focus of the development of countries around the world. As an important part of renewable energy, wind power has developed rapidly in recent years due to its low cost and mature technology (<xref ref-type="bibr" rid="B14">Chen et al., 2017</xref>; <xref ref-type="bibr" rid="B57">Yuan et al., 2022</xref>). According to statistics from the International Energy Agency, wind electricity generation reached 1,870&#xa0;TWh in 2021 and it remains the leading non-hydro renewable technology. To achieve the goal of net-zero emissions by 2050, which is to generate around 7,900&#xa0;TWh of wind power by 2030, it will be necessary to increase average annual electricity generation to almost 250&#xa0;GW (<xref ref-type="bibr" rid="B22">International Energy Agency, 2021</xref>). However, with the increasing penetration of wind power into the power grids, the randomness, volatility and intermittency of wind power bring great challenges to the safe and stable operation of the power grids (<xref ref-type="bibr" rid="B37">Shafiullah et al., 2013</xref>; <xref ref-type="bibr" rid="B16">Dai et al., 2019</xref>). Accurate and reliable wind power forecasting is an effective way to cope with this problem and has therefore become quite a hot topic of research (<xref ref-type="bibr" rid="B40">Tascikaraoglu and Uzunoglu, 2014</xref>; <xref ref-type="bibr" rid="B48">Wang et al., 2021</xref>).</p>
<p>According to the length of the foresight period, wind power forecasting can be divided into: ultra-short-term forecasting (0&#x2013;4&#xa0;h) for real-time load balancing, short-term forecasting (4&#x2013;72&#xa0;h) for unit commitment and flexibility reserve, and medium and long-term forecasting (several days, weeks or months) for unit maintenance scheduling and generation capacity evaluation. This study only focuses on the short-term wind power prediction. In recent years, many short-term wind power forecasting methods have been proposed. These can be summarized into three categories: physical methods, statistical methods, and ensemble forecasting models.</p>
<p>Based on the meteorological conditions of the underlying surface of the wind farms and the output curve of the fans, the physical prediction methods can establish the mapping relationship between wind power output and meteorological information using micro-meteorology to realize the wind power prediction. Numerical weather prediction (NWP) is the most commonly used physical method. <xref ref-type="bibr" rid="B10">Charabi et al. (2011)</xref> evaluated the performance of NWP model data for wind energy applications in Oman and demonstrated that NWP data has better accuracy than satellite data compared to ground measurements. <xref ref-type="bibr" rid="B30">Liu et al. (2022)</xref> proposed a novel NWP-enhanced wind power prediction method based on rank ensemble and probabilistic fluctuation awareness. <xref ref-type="bibr" rid="B36">Pr&#xf3;sper et al. (2019)</xref> focused on production prediction and validation of actual onshore wind farms using high horizontal and vertical resolution Weather Research and Prediction (WRF) model simulations. <xref ref-type="bibr" rid="B53">Ye et al. (2017)</xref> proposed a short-term wind power prediction model based on physical methods and spatial correlations to characterize the uncertainty and dependency structure of turbine&#x2019;s output in wind farms. However, physical methods rely on very precise meteorological and geographic data, which are sometimes difficult to obtain. In addition, the physical methods usually need significant computational time, making their application to short-term wind power forecasting difficult.</p>
<p>The statistical methods do not usually consider the complex physical mechanism of wind power generation, and only construct a statistical model based on the historical operational data of wind farms in order to achieve future wind power prediction. Compared with physical methods, the statistical methods have simpler calculation and can directly predict wind power by mapping the relationship between historical wind power data and the prediction target. Statistical models can be further divided into time series models, other machine learning models and deep learning models: 1) The typical time series models include the autoregressive moving average model (ARMA) (<xref ref-type="bibr" rid="B42">Torres et al., 2005</xref>), the autoregressive integrated moving average model (ARIMA) (<xref ref-type="bibr" rid="B12">Chen et al., 2010</xref>; <xref ref-type="bibr" rid="B7">Barbosa et al., 2017</xref>), the exponential smoothing method (<xref ref-type="bibr" rid="B8">Cadenas et al., 2010</xref>), and the generalized autoregressive conditional heteroscedasticity (GARCH) model (<xref ref-type="bibr" rid="B23">Jeon and Taylor, 2016</xref>). Nevertheless, time series models only analyze the potential relationship of time series variables, which makes it difficult for them to mine the nonlinear relationship between data, hence the prediction accuracy of this kind of model is poor. 2) Machine learning models can adaptively learn to make decisions and predict future data based on given historical data (<xref ref-type="bibr" rid="B31">Liu et al., 2019</xref>). Commonly used machine learning models, such as support vector machine (SVM) (<xref ref-type="bibr" rid="B32">Liu et al., 2017</xref>; <xref ref-type="bibr" rid="B2">Abedinia et al., 2022</xref>), random forest (RF) (<xref ref-type="bibr" rid="B26">Lahouar and Slama, 2017</xref>; <xref ref-type="bibr" rid="B39">Shi et al., 2018</xref>), and Bayesian additive regression tree (BART) (<xref ref-type="bibr" rid="B13">Chen et al., 2018</xref>), are widely used in wind power output prediction, wind speed prediction and other fields. However, the effect of SVM is closely related to the selection of kernel function and its parameters, which is strongly dependent on the user&#x2019;s experience. RF is prone to overfitting, and the BART method requires a long computation time. 3) With the rapid development of deep learning, artificial intelligence (AI) technology has also been applied to wind power prediction. The AI models, back-propagation (BP) neural network (<xref ref-type="bibr" rid="B59">Zhang et al., 2018</xref>), artificial neural network (ANN) (<xref ref-type="bibr" rid="B9">Carolin and Fernandez, 2008</xref>), convolution neural network (CNN) (<xref ref-type="bibr" rid="B43">Wang et al, 2017a</xref>; <xref ref-type="bibr" rid="B4">Afrasiabi et al., 2019</xref>) and recursive neural network (RNN) (<xref ref-type="bibr" rid="B27">Li et al., 2019</xref>) have been the focus of previous research on prediction models. These models have higher prediction accuracy than other machine learning models but have the same problem with difficulty in model training. Hence improved RNN and CNN models, such as long short-term memory (LSTM) (<xref ref-type="bibr" rid="B58">Zhang et al., 2019a</xref>; <xref ref-type="bibr" rid="B62">Zhang et al, 2019b</xref>; <xref ref-type="bibr" rid="B49">Wu et al., 2019</xref>), GRU (<xref ref-type="bibr" rid="B17">Ding et al., 2019</xref>; <xref ref-type="bibr" rid="B15">Chen et al., 2022</xref>), and temporal convolutional network (TCN) (<xref ref-type="bibr" rid="B19">Gan et al., 2021</xref>; <xref ref-type="bibr" rid="B20">He et al., 2022</xref>) have been widely used in wind power prediction. In recent years, the generative adversarial network (GAN) has attracted a lot of attention (<xref ref-type="bibr" rid="B56">Yuan et al., 2021</xref>; <xref ref-type="bibr" rid="B63">Zhou et al., 2021</xref>; <xref ref-type="bibr" rid="B51">Xia et al., 2022</xref>). Its generative model maps noise variables to multi-layer perceptron networks to make the generated data as close as possible to the distribution of training samples. In general, the AI models can better mine the hidden feature information of the wind power series, improve the overall prediction accuracy, and have strong learning ability and robustness.</p>
<p>Due to the high randomness and volatility of wind power, the prediction abilities of a single model often do not meet actual needs. In recent years, ensemble forecasting models which combine the advantages of multiple single models have become a popular direction for wind power prediction research. Current research on ensemble forecasting models can be summarized in four categories. 1) Ensemble forecasting models based on multi-model weighting. In these models, multiple single models, such as SVM and RNN (<xref ref-type="bibr" rid="B55">Yu et al., 2018</xref>), extreme learning machine (ELM), Elman neural network (ENN) and LSTM (<xref ref-type="bibr" rid="B1">Abedinia and Bagheri, 2022</xref>), least square SVM (LSSVM) and radial basis function neural network (RBFNN) (<xref ref-type="bibr" rid="B38">Shi et al., 2013</xref>), outlier robust ELM (ORELM), ENN, and bidirectional LSTM (BiLSTM) (<xref ref-type="bibr" rid="B11">Chen and Liu, 2020</xref>), are used to predict wind power series, and the prediction results are weighted to improve the prediction accuracy. 2) Ensemble forecasting models based on data preprocessing. To cope with the non-stationary wind power sequence, these methods use signal decomposition and denoising algorithms to decompose the original wind power data into multiple stationary subsequences, and use the prediction model to predict each subsequence separately. Commonly used mode decomposition algorithms include empirical mode decomposition (EMD) (<xref ref-type="bibr" rid="B6">Amjady and Abedinia, 2017</xref>; <xref ref-type="bibr" rid="B3">Abedinia et al., 2020</xref>), variational mode decomposition (VMD) (<xref ref-type="bibr" rid="B54">Yin et al., 2019</xref>; <xref ref-type="bibr" rid="B18">Duan et al., 2021</xref>), singular value decomposition (<xref ref-type="bibr" rid="B44">Wang et al., 2020</xref>), ensemble empirical mode decomposition (EEMD) (<xref ref-type="bibr" rid="B45">Wang et al, 2017b</xref>), wavelet transform (WT) (<xref ref-type="bibr" rid="B64">Zucatelli et al., 2021</xref>; <xref ref-type="bibr" rid="B25">Khazaei et al., 2022</xref>), and complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN) (<xref ref-type="bibr" rid="B34">Lu et al., 2020</xref>) Ensemble forecasting models based on optimization techniques. In order to improve the prediction accuracy, the parameters of the forecasting model are optimized by using optimization techniques. These models include the Multilayer Perceptron (MLP) neural network optimized by Non-dominated Sorting Genetic Algorithm II (NSGA- &#x4c0;&#x4c0;) (<xref ref-type="bibr" rid="B25">Khazaei et al., 2022</xref>), SVM optimized by cuckoo search algorithm (SVM-CSA) (<xref ref-type="bibr" rid="B28">Li et al., 2021</xref>), ENN optimized by multi-objective grey wolf optimization (ENN-MOGWO) (<xref ref-type="bibr" rid="B46">Wang et al., 2019a</xref>), ELM optimized by Particle Swarm Optimization (ELM-PSO) (<xref ref-type="bibr" rid="B41">Tian et al., 2019</xref>), Echo State Network optimized by MOGWO (ESN-MOGWO) (<xref ref-type="bibr" rid="B47">Wang et al., 2019b</xref>) Ensemble forecasting models based on error correction. In order to further reduce the prediction error, error correction technology has been widely used in wind power prediction, usually by predicting the error extracted from the initial prediction result as a secondary prediction. The Markov chain (MC) model (<xref ref-type="bibr" rid="B61">Zhang et al., 2014</xref>; <xref ref-type="bibr" rid="B60">Zhang et al., 2021</xref>), the GARCH (<xref ref-type="bibr" rid="B24">Jiang and Huang, 2017</xref>), the temporally local moving window technique (<xref ref-type="bibr" rid="B52">Yan et al., 2015</xref>), and machine learning methods (<xref ref-type="bibr" rid="B29">Liang et al., 2016</xref>) are commonly used to deal with the error component.</p>
<p>Although many advances have been made in wind power forecasting methods, wind power forecasting remains challenging due to the high instability of wind power output. Moreover, few prediction methods combine data decomposition, model prediction, and error correction techniques to further improve the prediction accuracy. Based on the above analysis, this research is driven by the following concepts: The EEMD method is an improved and robust decomposition technique, and can effectively discover the potential characteristics of wind power output; The GRU model shows good performance in extracting temporal correlation hidden features from time series, hence is making a figure in short-term power prediction of new energy sources; The MC approach is a very popular error correction technique because it is easy to understand and implement. Hence in this paper, following the concept of &#x201c;data decomposition - model prediction - error correction&#x201d;, a novel ensemble forecasting model for short-term wind power sequences based on EEMD-GRU-MC is developed. The proposed model consists of three important steps: Firstly, the EEMD method is employed to decompose the original wind power output sequence into a set of relatively stationary subsequences and denoise the data sequence. Secondly, the GRU model is used to individually forecast each subsequence, and the predicted value of each subsequence is superimposed to obtain the predicted result of the original data. Finally, to further enhance the prediction accuracy, the MC is applied to correct the preliminary prediction results. Extensive numerical experiments are conducted to test the performance of the proposed forecasting model when applied to different wind farms and in different seasons. This testing indicates that the proposed hybrid model outperforms the benchmark models in terms of multiple evaluation indicators. Moreover, the mean prediction error of the developed model in all scenarios is less than or close to 2%, proving that it is a promising prospect for short-term wind power prediction.</p>
<p>The rest of this paper is structured as follows: <xref ref-type="sec" rid="s2">Section 2</xref> introduces the proposed ensemble forecasting method for short-term wind power sequences. Case studies are presented and discussed in <xref ref-type="sec" rid="s3">Section 3</xref>. Finally, conclusions are drawn in <xref ref-type="sec" rid="s4">Section 4</xref>.</p>
</sec>
<sec sec-type="methods" id="s2">
<title>2 Methodology</title>
<sec id="s2-1">
<title>2.1 Data decomposition based on EEMD</title>
<p>EMD is a signal preprocessing analysis method proposed by <xref ref-type="bibr" rid="B50">Wu and Huang (2009)</xref>, which is widely used in non-stationary and nonlinear signal processing. It progressively breaks down fluctuations or trends in different frequencies in the signal, and finally obtains a set of intrinsic mode functions (IMFs), where each decomposed IMF represents the characteristic signals of different frequencies in the original signal. However, mode mixing may occur in EMD signal processing, which prevents the IMFs from being separated effectively. The EEMD method introduces Gaussian white noise into the original signal and realizes the automatic distribution of the signal for the appropriate timescale after several averaging calculations, which effectively solves the mode mixing problem. Wind power output is easily affected by wind direction, wind speed and other factors, and presents large random fluctuations, which result in a large number of outliers in the wind power sequence. Therefore, the EEMD algorithm is applied to decompose and denoise the wind power sequence, extract the main trend component in the sequence, and eliminate the random fluctuation component. The decomposition of the wind power sequence by EEMD can be summarized as the following steps:</p>
<p>
<statement content-type="step" id="Step_1">
<label>Step 1</label>
<p>The white noise signal <inline-formula id="inf1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is added to the original wind power sequence <inline-formula id="inf2">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, and the new power sequence <inline-formula id="inf3">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is obtained using Eq. <xref ref-type="disp-formula" rid="e1">1</xref>:<disp-formula id="e1">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2026;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>K</mml:mi>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
</statement>
</p>
<p>
<statement content-type="step" id="Step_2">
<label>Step 2</label>
<p>The new sequence <inline-formula id="inf4">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> (see Eq. <xref ref-type="disp-formula" rid="e2">2</xref>) is decomposed into a series of IMFs using the EMD algorithm (<xref ref-type="bibr" rid="B35">Naik et al., 2018</xref>):<disp-formula id="e2">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:msubsup>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>m</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
</statement>
</p>
<p>
<statement content-type="step" id="Step_3">
<label>Step 3</label>
<p>Steps (1) and (2) are repeated <italic>K</italic> times, and white noise with different amplitude is added each time.</p>
</statement>
</p>
<p>
<statement content-type="step" id="Step_4">
<label>Step 4</label>
<p>Since the mean value of the white noise spectrum is 0, the mean value of all IMFs calculated for K iterations is the final IMF obtained by the EEMD method (see Eq. <xref ref-type="disp-formula" rid="e3">3</xref>):<disp-formula id="e3">
<mml:math id="m7">
<mml:mrow>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>m</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>K</mml:mi>
</mml:msubsup>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>m</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:mo>/</mml:mo>
</mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
</statement>
</p>
<p>
<statement content-type="step" id="Step_5">
<label>Step 5</label>
<p>The original wind power sequence can be reconstructed as Eq. <xref ref-type="disp-formula" rid="e4">4</xref>:<disp-formula id="e4">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:msubsup>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>m</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>r</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>M</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>The amplitude of <inline-formula id="inf5">
<mml:math id="m9">
<mml:mrow>
<mml:msubsup>
<mml:mi>r</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>M</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is so small that it can be ignored in wind power prediction.</p>
</statement>
</p>
</sec>
<sec id="s2-2">
<title>2.2 Model prediction based on GRU neural network</title>
<p>LSTM is an enhanced type of RNN, which effectively solves the problem of the vanishing collateral gradient of traditional RNNs. GRU is an improved version of LSTM, simplifying the number of gating units and improving the computational efficiency of the model while ensuring the output accuracy. The GRU neuron is the basic unit of the GRU neural network (GRUNN) model and its structure is shown in <xref ref-type="fig" rid="F1">Figure 1</xref>. The GRU neuron includes reset gate <inline-formula id="inf6">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and update gate <inline-formula id="inf7">
<mml:math id="m11">
<mml:mrow>
<mml:msub>
<mml:mi>z</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. The update gate receives the current state <inline-formula id="inf8">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and the previously hidden state <inline-formula id="inf9">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. After receiving the input and the matrix operation, the sigmoid function &#x3c3; determines whether the neuron is activated. The reset gate receives <inline-formula id="inf10">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf11">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, and the result determines how much past information needs to be forgotten. The current memory <inline-formula id="inf12">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is a summary of the input and output of the previous hidden layer. <inline-formula id="inf13">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>h</mml:mi>
<mml:mo>&#x223c;</mml:mo>
</mml:mover>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf14">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> determine the final output <inline-formula id="inf15">
<mml:math id="m19">
<mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> by dynamic control of the update gate and transmit <inline-formula id="inf16">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> to the next GRU neuron. The mathematical model of GRU is shown as Eqs (<xref ref-type="disp-formula" rid="e5">5</xref>&#x2013;<xref ref-type="disp-formula" rid="e8">8</xref>).<disp-formula id="e5">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mi>z</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="normal">&#x3c3;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>z</mml:mi>
</mml:msub>
<mml:mo>&#x2219;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
<disp-formula id="e6">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>z</mml:mi>
</mml:msub>
<mml:mo>&#x2219;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
<disp-formula id="e7">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>h</mml:mi>
<mml:mo>&#x223c;</mml:mo>
</mml:mover>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="italic">tanh</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
<mml:mo>&#x2219;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x2299;</mml:mo>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
<disp-formula id="e8">
<mml:math id="m24">
<mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>z</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2299;</mml:mo>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mover accent="true">
<mml:mi>h</mml:mi>
<mml:mo>&#x223c;</mml:mo>
</mml:mover>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x2299;</mml:mo>
<mml:msub>
<mml:mi>z</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Structure of the GRU neuron.</p>
</caption>
<graphic xlink:href="fenrg-11-1252067-g001.tif"/>
</fig>
<p>Based on the GRU neuron, the time series prediction of GRUNN is shown in <xref ref-type="fig" rid="F2">Figure 2</xref>.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Schematic diagram of GRUNN prediction.</p>
</caption>
<graphic xlink:href="fenrg-11-1252067-g002.tif"/>
</fig>
</sec>
<sec id="s2-3">
<title>2.3 Detailed description of error correlation based on MC</title>
<sec id="s2-3-1">
<title>2.3.1 Basic theory of Markov chain</title>
<p>The Markov process is a typical stochastic process proposed by the famous mathematician Markov, which is applicable to both time series and interval sequences. The main content of the Markov process research is the state of a given stochastic process and its transition law. The MC refers to the Markov process with discrete time and state, and it can predict the changing trend of each state according to the initial probability of each state and the transition probability between each state. Hence the preliminary prediction results are corrected by MC to make up for the prediction error caused by the elimination of some components in the data decomposition process and the corrected wind power output is therefore closer to the actual value.</p>
<p>Assuming that <inline-formula id="inf17">
<mml:math id="m25">
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1,2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2026;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> is a random sequence where <italic>t</italic> represents any time period, if for any state <inline-formula id="inf18">
<mml:math id="m26">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#x2026;</mml:mo>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> state <italic>i</italic> and <italic>j</italic> satisfy Eq. <xref ref-type="disp-formula" rid="e9">9</xref>, then <inline-formula id="inf19">
<mml:math id="m27">
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1,2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2026;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> is a MC. <italic>i</italic> and <italic>j</italic> represent the possible states of the system at present and in future time, respectively:<disp-formula id="e9">
<mml:math id="m28">
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>j</mml:mi>
<mml:mo>&#x7c;</mml:mo>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#x2026;</mml:mo>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>j</mml:mi>
<mml:mo>&#x7c;</mml:mo>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>Supposing there are <italic>n</italic> states in state space <italic>I</italic> and, since each state can turn to itself, each state has <italic>n</italic> turns. So, the one-step transition probability from state <italic>i</italic> to state <italic>j</italic> can be expressed as Eq. <xref ref-type="disp-formula" rid="e10">10</xref>:<disp-formula id="e10">
<mml:math id="m29">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mi>M</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>M</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
<p>The matrix composed of the one-step transition probability set of all states is called the one-step transition probability matrix, and is expressed as Eq. <xref ref-type="disp-formula" rid="e11">11</xref>:<disp-formula id="e11">
<mml:math id="m30">
<mml:mrow>
<mml:msup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msup>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mn>11</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mn>12</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mn>21</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mn>22</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2026;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2026;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mo>&#x22ee;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:mo>&#x22ee;</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2026;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:mo>&#x22ee;</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2026;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
</p>
<p>Accordingly, the matrix composed of the <italic>k</italic>-step transition probabilities of all states is called the <italic>k</italic>-step transition probability matrix of the system. According to the homogeneity of MC, the <italic>k</italic>-step state transition probability matrix is expressed as Eq. <xref ref-type="disp-formula" rid="e12">12</xref>:<disp-formula id="e12">
<mml:math id="m31">
<mml:mrow>
<mml:msup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msup>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mn>11</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mn>12</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mn>21</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mn>22</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2026;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2026;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mo>&#x22ee;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:mo>&#x22ee;</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2026;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:mo>&#x22ee;</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2026;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>k</mml:mi>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
<p>In the process of MC error correction, the classification of states is very important. In this paper, the mean-standard deviation classification method, which is simple in theory and widely used, is employed to divide the state space according to the mean and standard deviation of the samples. Let the sample sequence be <inline-formula id="inf20">
<mml:math id="m32">
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>n</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1,2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2026;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>N</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>, the sample mean is <inline-formula id="inf21">
<mml:math id="m33">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula>, and the standard deviation is <inline-formula id="inf22">
<mml:math id="m34">
<mml:mrow>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>. According to the central limit theorem in mathematical statistics, the sample sequence can be divided into five intervals: <inline-formula id="inf23">
<mml:math id="m35">
<mml:mrow>
<mml:mfenced open="" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close="" separators="|">
<mml:mrow>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="italic">min</mml:mi>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#x2200;</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf24">
<mml:math id="m36">
<mml:mrow>
<mml:mfenced open="" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close="" separators="|">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.5</mml:mn>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf25">
<mml:math id="m37">
<mml:mrow>
<mml:mfenced open="" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close="" separators="|">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.5</mml:mn>
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.5</mml:mn>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf26">
<mml:math id="m38">
<mml:mrow>
<mml:mfenced open="" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close="" separators="|">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.5</mml:mn>
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf27">
<mml:math id="m39">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="italic">max</mml:mi>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#x2200;</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
</sec>
<sec id="s2-3-2">
<title>2.3.2 Basic theory of Markov chain</title>
<p>Based on the above analysis, the correction process for wind power prediction error is as follows:</p>
<p>
<statement content-type="step" id="step_1">
<label>Step 1</label>
<p>Calculate the historical wind power error sequence using Eq. <xref ref-type="disp-formula" rid="e13">13</xref>:<disp-formula id="e13">
<mml:math id="m40">
<mml:mrow>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>s</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:mo>&#x2200;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1,2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2026;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>where <inline-formula id="inf28">
<mml:math id="m41">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> denotes the actual historical wind power output value of sample point <italic>s</italic>; <inline-formula id="inf29">
<mml:math id="m42">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>s</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> denotes the predicted value of sample point <italic>s</italic> obtained by the EEMD-GRU model; <italic>S</italic> is the total number of sample points.</p>
</statement>
</p>
<p>
<statement content-type="step" id="step_2">
<label>Step 2</label>
<p>Calculate the mean and standard deviation of the error sequence and divide the error sequence into five intervals <inline-formula id="inf30">
<mml:math id="m43">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf31">
<mml:math id="m44">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> using the mean-standard deviation classification method.</p>
</statement>
</p>
<p>
<statement content-type="step" id="step_3">
<label>Step 3</label>
<p>The number of sample points belonging to different intervals is counted, and then the one-step and <italic>k</italic>-step transition probability matrices of each error state are calculated by using Eqs. <xref ref-type="disp-formula" rid="e11">11</xref> and <xref ref-type="disp-formula" rid="e2">(2)</xref>, respectively.</p>
</statement>
</p>
<p>
<statement content-type="step" id="step_4">
<label>Step 4</label>
<p>The states of the error sequence 5&#xa0;days before the forecast days are taken as the initial states. In the transition matrix <inline-formula id="inf32">
<mml:math id="m45">
<mml:mrow>
<mml:msup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula>, the row vectors corresponding to each initial state <inline-formula id="inf33">
<mml:math id="m46">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:mo>&#x22ef;</mml:mo>
<mml:mo>,</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>5</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2026;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> are taken to form a new probability matrix (see Eq. <xref ref-type="disp-formula" rid="e14">14</xref>):<disp-formula id="e14">
<mml:math id="m47">
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2026;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>5</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2026;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>5</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mo>&#x22ee;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:mo>&#x22ee;</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2026;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:mo>&#x22ee;</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2026;</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>5</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>
</p>
</statement>
</p>
<p>
<statement content-type="step" id="step_5">
<label>Step 5</label>
<p>The state corresponding to <inline-formula id="inf34">
<mml:math id="m48">
<mml:mrow>
<mml:mi mathvariant="italic">max</mml:mi>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mn>5</mml:mn>
</mml:munderover>
</mml:mstyle>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mn>5</mml:mn>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mn>1,5</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, i.e., the state to which the error is most likely to be transferred in the future, is taken as the state of the modified error. Thus the modified error <inline-formula id="inf35">
<mml:math id="m49">
<mml:mrow>
<mml:mover accent="true">
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>&#x223c;</mml:mo>
</mml:mover>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.5</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, where <inline-formula id="inf36">
<mml:math id="m50">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf37">
<mml:math id="m51">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the lower and upper bound of the state interval of the error to be modified.</p>
</statement>
</p>
<p>
<statement content-type="step" id="step_6">
<label>Step 6</label>
<p>Correct the predicted wind power sequence on the forecast days using Eq. <xref ref-type="disp-formula" rid="e15">15</xref>:<disp-formula id="e15">
<mml:math id="m52">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mover accent="true">
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>&#x223c;</mml:mo>
</mml:mover>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>
</p>
</statement>
</p>
</sec>
</sec>
<sec id="s2-4">
<title>2.4 Overall model prediction process</title>
<p>The overall flowchart of the proposed EEMD-GRU-MC model is depicted in <xref ref-type="fig" rid="F3">Figure 3</xref>, and the main steps are as follows:</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Flowchart of the proposed forecasting model.</p>
</caption>
<graphic xlink:href="fenrg-11-1252067-g003.tif"/>
</fig>
<p>
<statement content-type="step" id="STEP_1">
<label>Step 1</label>
<p>Use the EEMD method to decompose the historical power data into <italic>K</italic> IMF components (subsequences) and one RES component.</p>
</statement>
</p>
<p>
<statement content-type="step" id="STEP_2">
<label>Step 2</label>
<p>Divide each subsequence into a training set and a test set, and then use the GRU model to predict each component. The prediction results of each subsequence are aggregated as the preliminary prediction results of the EEMD-GRU model.</p>
</statement>
</p>
<p>
<statement content-type="step" id="STEP_3">
<label>Step 3</label>
<p>Calculate the prediction error between the historical wind power and the predicted power, then use the MC to correct the preliminary prediction results to get the final predicted wind power sequence.</p>
</statement>
</p>
</sec>
</sec>
<sec id="s3">
<title>3 Case studies</title>
<sec id="s3-1">
<title>3.1 Data description</title>
<p>To verify the effectiveness and practicability of the EEMD-GRU-MC model, it was applied to the short-term wind power prediction of two wind farms, ZMS and YMC, which are located in Yunnan Province, China. For each wind farm, four datasets were collected to test the forecasting performance of the proposed model in different seasons. The datasets were collected from 1 February 2021 to 22 March 2021, from 1 June 2021 to 20 July 2021, from 1 September 2021 to 20 October 2021 and from 1 December 2021 to 19 January 2022, representing the wind power data in spring, summer, autumn and winter, respectively. Each dataset was recorded for a time period of 15&#xa0;min. As shown in <xref ref-type="fig" rid="F4">Figures 4</xref>, <xref ref-type="fig" rid="F5">5</xref>, there is a total of 50 days, representing 4800 sample points, included in each dataset. The first 3840 sample points are used for the training set, the middle 864 are used for the validation set to avoid modeling over fitting, and the last 96 are used for the test set. This study focuses on short-term wind power prediction 24&#xa0;h in advance to assist day-ahead dispatching of power grids (<xref ref-type="bibr" rid="B5">Alham et al., 2016</xref>; <xref ref-type="bibr" rid="B21">Hu et al., 2019</xref>; <xref ref-type="bibr" rid="B33">Liu et al., 2021</xref>). Therefore, 96 sample points are selected as the test set in this study. The ratio of training set, verification set and test set is usually 6:2:2. In order to improve the prediction accuracy of the model, a longer training set and verification set were selected in this study, which made the data volume of the whole data set reach 4800. <xref ref-type="table" rid="T1">Table 1</xref> lists the statistical information for the datasets, including the maximum, minimum, mean, median and standard deviation. It can be observed that the power variation of each wind farm in all seasons is close to the installed capacity, showing strong volatility and non-stationarity. The power output of all wind farms is larger in spring and winter, but smaller in summer.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Wind power dataset for the ZMS wind farm.</p>
</caption>
<graphic xlink:href="fenrg-11-1252067-g004.tif"/>
</fig>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Wind power dataset for the YMS wind farm.</p>
</caption>
<graphic xlink:href="fenrg-11-1252067-g005.tif"/>
</fig>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Statistical information of the datasets used in this study.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Wind farm</th>
<th align="center">Dataset</th>
<th align="center">Maximum</th>
<th align="center">Minimum</th>
<th align="center">Range</th>
<th align="center">Mean</th>
<th align="center">Median</th>
<th align="center">Std.</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="4" align="center">ZMS</td>
<td align="center">Spring</td>
<td align="center">145.2</td>
<td align="center">0</td>
<td align="center">145.2</td>
<td align="center">101.4</td>
<td align="center">114.4</td>
<td align="center">37.2</td>
</tr>
<tr>
<td align="center">Summer</td>
<td align="center">147.6</td>
<td align="center">0</td>
<td align="center">147.6</td>
<td align="center">66.8</td>
<td align="center">61.9</td>
<td align="center">41.9</td>
</tr>
<tr>
<td align="center">Autumn</td>
<td align="center">131.8</td>
<td align="center">0</td>
<td align="center">131.8</td>
<td align="center">29.4</td>
<td align="center">21.4</td>
<td align="center">27.9</td>
</tr>
<tr>
<td align="center">Winter</td>
<td align="center">143.3</td>
<td align="center">0</td>
<td align="center">143.3</td>
<td align="center">68.5</td>
<td align="center">67.3</td>
<td align="center">38.7</td>
</tr>
<tr>
<td rowspan="4" align="center">YMS</td>
<td align="center">Spring</td>
<td align="center">83.0</td>
<td align="center">0.1</td>
<td align="center">82.9</td>
<td align="center">53.7</td>
<td align="center">58.5</td>
<td align="center">20.9</td>
</tr>
<tr>
<td align="center">Summer</td>
<td align="center">81.8</td>
<td align="center">0</td>
<td align="center">81.8</td>
<td align="center">18.9</td>
<td align="center">11.5</td>
<td align="center">19.7</td>
</tr>
<tr>
<td align="center">Autumn</td>
<td align="center">88.6</td>
<td align="center">0</td>
<td align="center">88.6</td>
<td align="center">22.6</td>
<td align="center">20.8</td>
<td align="center">17.4</td>
</tr>
<tr>
<td align="center">Winter</td>
<td align="center">98.5</td>
<td align="center">0</td>
<td align="center">98.5</td>
<td align="center">48.5</td>
<td align="center">51.7</td>
<td align="center">21.6</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The root mean square error (RMSE), mean absolute error (MAE) and mean absolute percentage error (MAPE) are used as indexes to evaluate the predictive performance of the forecasting models (see Eqs <xref ref-type="disp-formula" rid="e16">16</xref>&#x2013;<xref ref-type="disp-formula" rid="e18">18</xref>). The smaller the RMSE, MAE and MAPE, the better the predictive performance of the models:<disp-formula id="e16">
<mml:math id="m53">
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>M</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>E</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:msubsup>
</mml:mstyle>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>a</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo>/</mml:mo>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
<label>(16)</label>
</disp-formula>
<disp-formula id="e17">
<mml:math id="m54">
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi>A</mml:mi>
<mml:mi>E</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:msubsup>
</mml:mstyle>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>a</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mo>/</mml:mo>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:math>
<label>(17)</label>
</disp-formula>
<disp-formula id="e18">
<mml:math id="m55">
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi>A</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>E</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:msubsup>
</mml:mstyle>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>a</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mo>/</mml:mo>
</mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>a</mml:mi>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mo>/</mml:mo>
</mml:mrow>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:math>
<label>(18)</label>
</disp-formula>
</p>
<p>Where <italic>L</italic> is the number of sample points in the test set, which is 96 in this paper; <inline-formula id="inf38">
<mml:math id="m56">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>a</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the actual wind power value of sample point <italic>l</italic>; <inline-formula id="inf39">
<mml:math id="m57">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the predicted value.</p>
<p>The EEMD decomposition of the wind power sequence and MC error correction of the preliminary prediction results were realized by Matlab 2020a, and the training and prediction of the GRU model were realized by Python programming language. All numerical experiments were conducted on a Dell workstation equipped with an Intel Xeon Gold processor, with 20 cores and 40 threads, 2.1G main frequency and 64G memory.</p>
</sec>
<sec id="s3-2">
<title>3.2 Case 1: short-term wind power prediction for the ZMS wind farm</title>
<p>Aiming to solve the problem of poor model robustness caused by the randomicity and intermittent nature of wind power, the EEMD algorithm was introduced to decompose wind power data into a set of subsequences. Due to space limitations, <xref ref-type="sec" rid="s10">Supplementary Figure S1</xref> only displays the EEMD decomposition results of the spring dataset. It can be seen that EEMD decomposes the wind power sequence into seven IMF subsequences and one residual subsequence with different frequency characteristics, which facilitates the analysis of the hidden information in the data and overcomes the shortcomings of the original wind power sequence with its high volatility and non-stationarity.</p>
<p>In order to verify the superiority, reliability and stability of the proposed model, six other forecasting models based on LSTM, GRU, EMD-LSTM, EMD-GRU, EEMD-LSTM and EEMD-GRU were constructed as comparison models. It should be mentioned that the LSTM and GRU methods are adopted by (<xref ref-type="bibr" rid="B18">Duan et al., 2021</xref>; <xref ref-type="bibr" rid="B15">Chen et al., 2022</xref>), respectively. The prediction results of the seven models in different seasons are shown in <xref ref-type="fig" rid="F6">Figure 6</xref>, and regression analysis of the prediction results for the ZMS wind farm using different models is presented in <xref ref-type="fig" rid="F7">Figure 7</xref>. For further quantitative comparison, the evaluation indicators of various prediction models, including MAE, RMSE and MAPE are listed in <xref ref-type="table" rid="T2">Table 2</xref>. The detailed analyses are summarized as follows: (1) The predicted wind power curves of all models are generally consistent with the trend of the actual power curve. However, it is clear that the predicted wind power curve obtained by the proposed EEMD-GRU-MC model is very close to the actual power curve, and has the smallest RMSE, MAE and MAPE among the seven models for all seasons. In addition, the correlation between the observed data and the predicted data generated by the proposed model is greater than that of comparison models. Therefore, the proposed model connecting EEMD and MC to the GRU model has the ability to capture the dynamic characteristics of wind power output data series. (2) The LSTM has the largest prediction error and the prediction effect of EEMD-LSTM is also inferior to that of EEMD-GRU for different seasons, which proves that GRU has more advantages in predicting short-term wind power time series data compared to LSTM. (3) EEMD-LSTM and EEMD-GRU models are superior to LSTM and GRU respectively in various performance evaluation indexes. Taking the summer dataset with the strongest stochastic wind power volatility as an example, the MAE, RSME and MAPE of the EEMD-GRU model are 1.90 MW, 3.36&#xa0;MW and 1.58%, which decreased by 69.3%, 48.7% and 71.2% compared with the GRU model. Similar results also appear in the comparison of the EEMD-LSTM and LSTM model, whose MAE, RMSE and MAPE decreased by 77.9%, 47.3% and 79.2%, respectively. This proves that EEMD can separate the noise information from the complex wind power data and facilitate the prediction model to extract the hidden information in the data. (4) Compared with EMD-LSTM and EMD-GRU, EEMD-LSTM and EEMD-GRU have better predictive performance, showing EEMD technique is more helpful for improving the prediction accuracy than EMD technique. (5) The prediction accuracy of the EEMD-GRU model can be further improved after MC error correlation. Taking the spring dataset, for example, after MC correction, the MAE, RMSE and MAPE of the EEMD-GRU model decreased from 1.87&#xa0;MW to 1.37 MW, from 2.37&#xa0;MW to 1.97 MW, and from 2.18% to 1.76%, respectively.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Short-term wind power forecasting results for the ZMS wind farm.</p>
</caption>
<graphic xlink:href="fenrg-11-1252067-g006.tif"/>
</fig>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>
<bold>(A)</bold> Spring Regression analysis of prediction results for the ZMS wind farm using different models. <bold>(B)</bold> Summer Regression analysis of prediction results for the ZMS wind farm using different models. <bold>(C)</bold> Autumn Regression analysis of prediction results for the ZMS wind farm using different models. <bold>(D)</bold> Winter Regression analysis of prediction results for the ZMS wind farm using different models.</p>
</caption>
<graphic xlink:href="fenrg-11-1252067-g007.tif"/>
</fig>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Statistical indexes of short-term power prediction for the ZMS wind farm using different models.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Dataset</th>
<th align="center">Model</th>
<th align="center">MAE(MW)</th>
<th align="center">RMSE(MW)</th>
<th align="center">MAPE (%)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="7" align="center">Spring</td>
<td align="center">LSTM</td>
<td align="center">9.82</td>
<td align="center">11.50</td>
<td align="center">10.74</td>
</tr>
<tr>
<td align="center">GRU</td>
<td align="center">7.49</td>
<td align="center">8.61</td>
<td align="center">8.32</td>
</tr>
<tr>
<td align="center">EMD-LSTM</td>
<td align="center">6.13</td>
<td align="center">7.26</td>
<td align="center">7.27</td>
</tr>
<tr>
<td align="center">EMD-GRU</td>
<td align="center">3.61</td>
<td align="center">4.45</td>
<td align="center">4.38</td>
</tr>
<tr>
<td align="center">EEMD-LSTM</td>
<td align="center">4.03</td>
<td align="center">5.15</td>
<td align="center">4.59</td>
</tr>
<tr>
<td align="center">EEMD-GRU</td>
<td align="center">1.87</td>
<td align="center">2.37</td>
<td align="center">2.18</td>
</tr>
<tr>
<td align="center">EEMD-GRU-MC</td>
<td align="center">1.37</td>
<td align="center">1.97</td>
<td align="center">1.76</td>
</tr>
<tr>
<td rowspan="7" align="center">Summer</td>
<td align="center">LSTM</td>
<td align="center">9.62</td>
<td align="center">11.11</td>
<td align="center">8.35</td>
</tr>
<tr>
<td align="center">GRU</td>
<td align="center">6.19</td>
<td align="center">6.55</td>
<td align="center">5.49</td>
</tr>
<tr>
<td align="center">EMD-LSTM</td>
<td align="center">3.58</td>
<td align="center">7.97</td>
<td align="center">2.90</td>
</tr>
<tr>
<td align="center">EMD-GRU</td>
<td align="center">2.27</td>
<td align="center">6.03</td>
<td align="center">1.84</td>
</tr>
<tr>
<td align="center">EEMD-LSTM</td>
<td align="center">2.13</td>
<td align="center">5.86</td>
<td align="center">1.74</td>
</tr>
<tr>
<td align="center">EEMD-GRU</td>
<td align="center">1.90</td>
<td align="center">3.36</td>
<td align="center">1.58</td>
</tr>
<tr>
<td align="center">EEMD-GRU-MC</td>
<td align="center">1.43</td>
<td align="center">1.41</td>
<td align="center">1.28</td>
</tr>
<tr>
<td rowspan="7" align="center">Autumn</td>
<td align="center">LSTM</td>
<td align="center">5.45</td>
<td align="center">7.41</td>
<td align="center">9.69</td>
</tr>
<tr>
<td align="center">GRU</td>
<td align="center">4.07</td>
<td align="center">5.52</td>
<td align="center">7.43</td>
</tr>
<tr>
<td align="center">EMD-LSTM</td>
<td align="center">4.37</td>
<td align="center">5.72</td>
<td align="center">8.15</td>
</tr>
<tr>
<td align="center">EMD-GRU</td>
<td align="center">3.95</td>
<td align="center">5.67</td>
<td align="center">7.00</td>
</tr>
<tr>
<td align="center">EEMD-LSTM</td>
<td align="center">2.68</td>
<td align="center">4.32</td>
<td align="center">4.67</td>
</tr>
<tr>
<td align="center">EEMD-GRU</td>
<td align="center">2.24</td>
<td align="center">3.22</td>
<td align="center">4.49</td>
</tr>
<tr>
<td align="center">EEMD-GRU-MC</td>
<td align="center">0.76</td>
<td align="center">0.84</td>
<td align="center">1.69</td>
</tr>
<tr>
<td rowspan="7" align="center">Winter</td>
<td align="center">LSTM</td>
<td align="center">6.08</td>
<td align="center">6.33</td>
<td align="center">11.53</td>
</tr>
<tr>
<td align="center">GRU</td>
<td align="center">2.64</td>
<td align="center">2.90</td>
<td align="center">4.94</td>
</tr>
<tr>
<td align="center">EMD-LSTM</td>
<td align="center">2.48</td>
<td align="center">3.01</td>
<td align="center">4.62</td>
</tr>
<tr>
<td align="center">EMD-GRU</td>
<td align="center">1.60</td>
<td align="center">2.52</td>
<td align="center">3.00</td>
</tr>
<tr>
<td align="center">EEMD-LSTM</td>
<td align="center">1.50</td>
<td align="center">2.07</td>
<td align="center">2.85</td>
</tr>
<tr>
<td align="center">EEMD-GRU</td>
<td align="center">0.75</td>
<td align="center">1.15</td>
<td align="center">1.33</td>
</tr>
<tr>
<td align="center">EEMD-GRU-MC</td>
<td align="center">0.03</td>
<td align="center">0.07</td>
<td align="center">0.68</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-3">
<title>3.3 Case 2: Short-term wind power prediction of the YMS wind farm</title>
<p>In order to verify its robustness, the proposed model was applied to the short-term power prediction of the YMS wind farm, whose output characteristics are quite different from those of ZMS. The experiments were also conducted using six comparison models as mentioned in <xref ref-type="sec" rid="s3-1">Section 3.1</xref>. The forecasting results from these models and regression analysis are shown in <xref ref-type="fig" rid="F8">Figures 8</xref>, <xref ref-type="fig" rid="F9">9</xref>, while the evaluation indicators of the forecasting results are illustrated in <xref ref-type="table" rid="T3">Table 3</xref>. It can be found that the LSTM and GRU based models without data preprocessing fail to obtain satisfactory forecasting results. Especially in spring and summer, the MAPE values of both models exceed 10%. The proposed EEMD-GRU-MC model achieves the smallest MAE, RMSE and MAPE among the models for the four seasons, and the forecasted wind power curves closely match the trend of the actual power curves. Except for in spring, the MAPE values of the predicted results of the developed model are all within 2%. Although the developed model&#x2019;s prediction accuracy dropped slightly in the spring dataset where the wind power fluctuation is more severe, the proposed model performs best in the wind power prediction for the YMS wind farm. It can be concluded that the developed ensemble forecasting model has more outstanding potential and is better able to capture valuable information in complex and non-stationary wind power data.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Short-term wind power forecasting results for the YMS wind farm.</p>
</caption>
<graphic xlink:href="fenrg-11-1252067-g008.tif"/>
</fig>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>
<bold>(A)</bold> Spring Regression analysis of prediction results for the YMS wind farm using different models. <bold>(B)</bold> Summer Regression analysis of prediction results for the YMS wind farm using different models. <bold>(C)</bold> Autumn Regression analysis of prediction results for the YMS wind farm using different models. <bold>(D)</bold> Winter Regression analysis of prediction results for the YMS wind farm using different models.</p>
</caption>
<graphic xlink:href="fenrg-11-1252067-g009.tif"/>
</fig>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Statistical indexes of short-term power prediction for YMS wind farm using different models.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Dataset</th>
<th align="center">Model</th>
<th align="center">MAE(MW)</th>
<th align="center">RMSE(MW)</th>
<th align="center">MAPE (%)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="7" align="center">Spring</td>
<td align="center">LSTM</td>
<td align="center">3.59</td>
<td align="center">4.87</td>
<td align="center">11.58</td>
</tr>
<tr>
<td align="center">GRU</td>
<td align="center">2.55</td>
<td align="center">3.22</td>
<td align="center">10.07</td>
</tr>
<tr>
<td align="center">EMD-LSTM</td>
<td align="center">2.38</td>
<td align="center">3.14</td>
<td align="center">12.87</td>
</tr>
<tr>
<td align="center">EMD-GRU</td>
<td align="center">1.95</td>
<td align="center">2.68</td>
<td align="center">10.61</td>
</tr>
<tr>
<td align="center">EEMD-LSTM</td>
<td align="center">1.49</td>
<td align="center">1.96</td>
<td align="center">6.94</td>
</tr>
<tr>
<td align="center">EEMD-GRU</td>
<td align="center">1.14</td>
<td align="center">1.52</td>
<td align="center">5.86</td>
</tr>
<tr>
<td align="center">EEMD-GRU-MC</td>
<td align="center">0.32</td>
<td align="center">0.39</td>
<td align="center">2.55</td>
</tr>
<tr>
<td rowspan="7" align="center">Summer</td>
<td align="center">LSTM</td>
<td align="center">2.36</td>
<td align="center">2.47</td>
<td align="center">16.89</td>
</tr>
<tr>
<td align="center">GRU</td>
<td align="center">1.71</td>
<td align="center">1.90</td>
<td align="center">12.05</td>
</tr>
<tr>
<td align="center">EMD-LSTM</td>
<td align="center">1.34</td>
<td align="center">1.60</td>
<td align="center">9.90</td>
</tr>
<tr>
<td align="center">EMD-GRU</td>
<td align="center">1.00</td>
<td align="center">1.33</td>
<td align="center">6.57</td>
</tr>
<tr>
<td align="center">EEMD-LSTM</td>
<td align="center">1.41</td>
<td align="center">1.57</td>
<td align="center">11.55</td>
</tr>
<tr>
<td align="center">EEMD-GRU</td>
<td align="center">0.92</td>
<td align="center">1.12</td>
<td align="center">7.29</td>
</tr>
<tr>
<td align="center">EEMD-GRU-MC</td>
<td align="center">0.15</td>
<td align="center">0.32</td>
<td align="center">1.12</td>
</tr>
<tr>
<td rowspan="7" align="center">Autumn</td>
<td align="center">LSTM</td>
<td align="center">4.57</td>
<td align="center">5.10</td>
<td align="center">7.63</td>
</tr>
<tr>
<td align="center">GRU</td>
<td align="center">3.20</td>
<td align="center">3.61</td>
<td align="center">5.16</td>
</tr>
<tr>
<td align="center">EMD-LSTM</td>
<td align="center">3.31</td>
<td align="center">3.86</td>
<td align="center">6.14</td>
</tr>
<tr>
<td align="center">EMD-GRU</td>
<td align="center">1.81</td>
<td align="center">2.31</td>
<td align="center">3.40</td>
</tr>
<tr>
<td align="center">EEMD-LSTM</td>
<td align="center">1.68</td>
<td align="center">2.03</td>
<td align="center">2.78</td>
</tr>
<tr>
<td align="center">EEMD-GRU</td>
<td align="center">1.26</td>
<td align="center">1.70</td>
<td align="center">2.47</td>
</tr>
<tr>
<td align="center">EEMD-GRU-MC</td>
<td align="center">0.61</td>
<td align="center">0.77</td>
<td align="center">1.15</td>
</tr>
<tr>
<td rowspan="7" align="center">Winter</td>
<td align="center">LSTM</td>
<td align="center">4.52</td>
<td align="center">5.57</td>
<td align="center">6.35</td>
</tr>
<tr>
<td align="center">GRU</td>
<td align="center">3.52</td>
<td align="center">4.69</td>
<td align="center">4.99</td>
</tr>
<tr>
<td align="center">EMD-LSTM</td>
<td align="center">3.63</td>
<td align="center">5.28</td>
<td align="center">5.29</td>
</tr>
<tr>
<td align="center">EMD-GRU</td>
<td align="center">2.85</td>
<td align="center">5.06</td>
<td align="center">4.11</td>
</tr>
<tr>
<td align="center">EEMD-LSTM</td>
<td align="center">2.77</td>
<td align="center">4.36</td>
<td align="center">4.15</td>
</tr>
<tr>
<td align="center">EEMD-GRU</td>
<td align="center">2.61</td>
<td align="center">3.95</td>
<td align="center">4.00</td>
</tr>
<tr>
<td align="center">EEMD-GRU-MC</td>
<td align="center">0.98</td>
<td align="center">1.15</td>
<td align="center">1.44</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-4">
<title>3.3 Analysis of computational efficiency of the proposed model</title>
<p>Computational complexity is an important index to evaluate the efficiency of a wind power forecasting method, which describes how the execution time of a method changes with the increase of the input size. The computational complexity of EEMD-GRU-MC model mainly depends on the respective complexity of EEMD, GRU and MC. EEMD is an adaptive signal processing method for nonlinear and non-stationary data. The time complexity of EEMD is actually equivalent to the time complexity of Fourier transform. This means that although EEMD is considered computationally intensive, it is actually a computationally efficient method. A GRU is a recurrent neural network used to process sequential data. The computational complexity of a GRU depends on several factors, including sequence length, number of network layers, and number of hidden units per layer. In general, the computational complexity of a GRU is proportional to these factors. An MC is a statistical model that describes random changes in the state of a system. The computational complexity of MC depends on the number of states. If the number of states is fixed, then the computational complexity of MC can be considered constant. In general, computational complexity is related to the parameters of the model and the amount of data.</p>
<p>In order to ensure that the proposed prediction model is realizable in practical applications, the prediction efficiency of the proposed model is analyzed, which is shown in <xref ref-type="table" rid="T4">Table 4</xref>. This model is designed to predict the day-ahead wind power output, rather than real-time prediction. After 30 iterations, the model uses an average of about 35&#xa0;min to accurately predict the wind power of the next day. Considering that the accuracy of the model prediction is very high, the prediction time is completely acceptable and can meet the timeliness requirement of the short-term wind power prediction.</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Prediction time for different wind farms in different seasons.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="center">Wind farm</th>
<th rowspan="2" align="center">Dataset</th>
<th colspan="3" align="center">Time (s)</th>
</tr>
<tr>
<th align="center">Decomposition</th>
<th align="center">Total prediction</th>
<th align="center">Average</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="4" align="center">ZMS</td>
<td align="center">Spring</td>
<td align="center">10.135</td>
<td align="center">2,135</td>
<td rowspan="4" align="center">2,140</td>
</tr>
<tr>
<td align="center">Summer</td>
<td align="center">10.158</td>
<td align="center">2,124</td>
</tr>
<tr>
<td align="center">Autumn</td>
<td align="center">10.206</td>
<td align="center">2,138</td>
</tr>
<tr>
<td align="center">Winter</td>
<td align="center">10.319</td>
<td align="center">2,162</td>
</tr>
<tr>
<td rowspan="4" align="center">YMS</td>
<td align="center">Spring</td>
<td align="center">10.159</td>
<td align="center">2,141</td>
<td rowspan="4" align="center">2,154</td>
</tr>
<tr>
<td align="center">Summer</td>
<td align="center">10.171</td>
<td align="center">2,136</td>
</tr>
<tr>
<td align="center">Autumn</td>
<td align="center">10.238</td>
<td align="center">2,150</td>
</tr>
<tr>
<td align="center">Winter</td>
<td align="center">10.338</td>
<td align="center">2,170</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Few prediction methods combine data decomposition, model prediction, and error correction techniques to further improve the prediction accuracy. The EEMD-GRU-MC model proposed in this paper provides a novel method for wind power prediction. By integrating EEMD algorithm, GRU model and MC technology, this model effectively deals with the strong intermittency and large volatility of wind power, thus improving the accuracy of model prediction. This novel method provides a new perspective and idea for the theoretical research of wind power prediction. For the research community, the EEMD-GRU-MC model has enriched the theoretical research of wind power prediction, and provided a new reference and inspiration for the subsequent research. For practitioners, especially power grid dispatching departments, the research results of this paper can help them make more accurate and reliable short-term wind power forecasts, so as to arrange more reasonable day-ahead generation plans. Moreover, the wind farms&#x2019; historical operation data and the source code of the proposed EEMD-GRU-MC model will be made available on request.</p>
</sec>
</sec>
<sec sec-type="conclusion" id="s4">
<title>4 Conclusion</title>
<p>Accurate and reliable short-term prediction of wind power is of important reference value for the power grid dispatching department to arrange reasonable day-ahead generation plans. This study innovatively combines a data decomposition technique, an AI-based prediction model and an error correction technique, and proposes a short-term wind power prediction method based on EMD-GRU-MC. Two case studies, including two wind farms and eight datasets, are used to verify the performance of the proposed forecasting model when applied to different wind farms and in different seasons. The conclusions can be summarized as follows:<list list-type="simple">
<list-item>
<p>(1) Compared with LSTM, GRU, EMD-LSTM, EMD-GRU, EEMD-LSTM and EEMD-GRU models, the proposed EEMD-GRU-MC model achieves the smallest MAE, RMSE and MAPE for all datasets, and the forecasted wind power curves very closely match the trend of the actual power curves. Taking the spring dataset of the ZMS wind farm for example, the MAE, RMSE and MAPE of the EEMD-GRU-MC model is 1.37 MW, 1.97 MW, and from 1.76%, respectively. Moreover, except for the YMS wind farm in spring, the mean forecasting error of the proposed model is always within 2%. This demonstrates that the proposed model has excellent forecasting performance and generalization ability, and can be used as an effective tool for short-term wind power prediction.</p>
</list-item>
<list-item>
<p>(2) After 30 iterations, the proposed model uses an average of about 35&#xa0;min to accurately predict the wind power of the next day, proving its high computation efficiency.</p>
</list-item>
<list-item>
<p>(3) GRU has more advantages in predicting short-term wind power sequences than LSTM. EEMD-LSTM and EEMD-GRU models are also achieve better prediction performance than LSTM and GRU respectively in various scenarios, indicating that the EEMD algorithm can overcome the shortcomings of the original wind power sequence with its high volatility and non-stationarity and facilitate the prediction model to extract the hidden information in the data. Taking the summer dataset of the ZMS wind farm as an example, the MAE, RSME and MAPE of the EEMD-GRU model are 1.90 MW, 3.36&#xa0;MW and 1.58%, which decreased by 69.3%, 48.7% and 71.2%, respectively, compared with the GRU model.</p>
</list-item>
<list-item>
<p>(4) For the spring dataset of the ZMS wind farm, the MAE, RMSE and MAPE of the EEMD-GRU model decreased by 26.73%, 16.88% and 19.27%, respectively, after MC correction. Similar results also appear in other datasets. This proves the effectiveness and applicability of the MC error correlation technique in short-term wind power forecasting.</p>
</list-item>
</list>
</p>
<p>The proposed EEMD-GRU-MC model is a deterministic wind power forecasting model, and does not take into account the complex meteorological factors. Moreover, the process of dissecting and projecting all of the data is not online forecasting, resulting in its temporary can not be applied to real-time forecasting. In future studies, the uncertainty of wind power prediction error will be considered and meteorological factors will be embedded to build a multi-feature interval prediction model, so as to obtain more comprehensive wind power prediction results. Moreover, the number of decomposed IMFs of EEMD and CEEMDAN, along with the standard EMD and its upgraded algorithms, is uncertain for different data characteristics. The developed model and the commercial solver Matlab 2020a can be integrated into the wind farm short-term power prediction support systems. In this case, the system program can automatically identify each IMF decomposed by the EEMD method, and then give it to the GRU model one by one for training and prediction. The support system would cope with the problem that the number of decomposed IMFs is uncertain for different data characteristics and realize online wind power forecasting. Hence, how to integrate the hybrid prediction model based on EEMD-GRU-MC into the decision support system for the short-term power prediction of wind farms to realize online and real-time prediction is our next research direction.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s5">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="sec" rid="s10">Supplementary Material</xref>, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s6">
<title>Author contributions</title>
<p>PW: Software, Writing&#x2013;Original Draft, Conceptualization. CS: Methodology, Writing&#x2013;Review and Editing, Funding acquisition, Supervision. LL: Validation, Writing&#x2013;Review and Editing. WY: Investigation, Validation. CG: Data Curation, Formal analysis. All authors contributed to manuscript revision, read, and approved the submitted version.</p>
</sec>
<sec sec-type="funding-information" id="s7">
<title>Funding</title>
<p>The research presented in this paper was supported by the National Natural Science Foundation of China (No. 52109041), China Postdoctoral Science Foundation (No. 2021M690139) and First-class Project Special Funding of Yellow River Laboratory (No. YRL22LT08).</p>
</sec>
<sec sec-type="COI-statement" id="s8">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s9">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s10">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fenrg.2023.1252067/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fenrg.2023.1252067/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Image1.TIF" id="SM1" mimetype="application/TIF" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abedinia</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Bagheri</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Execution of synthetic Bayesian model average for solar energy forecasting</article-title>. <source>IET Renew. Power Gener.</source> <volume>16</volume> (<issue>6</issue>), <fpage>1134</fpage>&#x2013;<lpage>1147</lpage>. <pub-id pub-id-type="doi">10.1049/rpg2.12389</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Abedinia</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Ghasemi-Marzbali</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Shafiei</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Sobhani</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Gharehpetian</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Bagheri</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2022</year>). &#x201c;<article-title>A multi-level model for hybrid short term wind forecasting based on SVM, wavelet transform and feature selection</article-title>,&#x201d; in <conf-name>Proceedings of the 2022 IEEE International Conference on Environment and Electrical Engineering and 2022 IEEE Industrial and Commercial Power Systems Europe (EEEIC/I&#x26;CPS Europe)</conf-name>, <conf-loc>Prague, Czech Republic</conf-loc>, <conf-date>July 2022</conf-date> (<publisher-name>IEEE</publisher-name>), <fpage>1</fpage>&#x2013;<lpage>6</lpage>. <pub-id pub-id-type="doi">10.1109/EEEIC/ICPSEUROPE54979.2022.9854519</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abedinia</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Lotfi</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Bagheri</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Sobhani</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Shafie-khah</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Catalao</surname>
<given-names>J. P. S.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Improved EMD-based complex prediction model for wind power forecasting</article-title>. <source>IEEE Trans. Sustain. Energy</source> <volume>11</volume> (<issue>4</issue>), <fpage>2790</fpage>&#x2013;<lpage>2802</lpage>. <pub-id pub-id-type="doi">10.1109/TSTE.2020.2976038</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Afrasiabi</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Mohammadi</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Rastegar</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Kargarian</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Probabilistic deep neural network price forecasting based on residential load and wind speed predictions</article-title>. <source>IET Renew. POWER Gen.</source> <volume>13</volume> (<issue>11</issue>), <fpage>1840</fpage>&#x2013;<lpage>1848</lpage>. <pub-id pub-id-type="doi">10.1049/iet-rpg.2018.6257</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alham</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Elshahed</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Ibrahim</surname>
<given-names>D. K.</given-names>
</name>
<name>
<surname>Zahab</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>A dynamic economic emission dispatch considering wind power uncertainty incorporating energy storage system and demand side management</article-title>. <source>Renew. energy</source> <volume>96</volume>, <fpage>800</fpage>&#x2013;<lpage>811</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2016.05.012</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Amjady</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Abedinia</surname>
<given-names>O.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Short term wind power prediction based on improved Kriging interpolation, empirical mode decomposition, and closed-loop forecasting engine</article-title>. <source>Sustainability</source> <volume>9</volume> (<issue>11</issue>), <fpage>2104</fpage>. <pub-id pub-id-type="doi">10.3390/su9112104</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Barbosa</surname>
<given-names>A. D.</given-names>
</name>
<name>
<surname>de Mattos Affonso</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Lima&#x2dc;o</surname>
<given-names>O. R.</given-names>
</name>
<name>
<surname>Moya</surname>
<given-names>R. J.</given-names>
</name>
<name>
<surname>Leite</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Reston</surname>
<given-names>F. J.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Different models for forecasting wind power generation: case study</article-title>. <source>ENERGIES</source> <volume>10</volume> (<issue>12</issue>), <fpage>1976</fpage>. <pub-id pub-id-type="doi">10.3390/en10121976</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cadenas</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Jaramillo</surname>
<given-names>O. A.</given-names>
</name>
<name>
<surname>Rivera</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Analysis and forecasting of wind velocity in chetumal, quintana roo, using the single exponential smoothing method</article-title>. <source>Renew. ENERGY</source> <volume>35</volume> (<issue>5</issue>), <fpage>925</fpage>&#x2013;<lpage>930</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2009.10.037</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Carolin</surname>
<given-names>M. M.</given-names>
</name>
<name>
<surname>Fernandez</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Analysis of wind power generation and prediction using ANN: a case study</article-title>. <source>Renew. ENERGY</source> <volume>33</volume> (<issue>5</issue>), <fpage>986</fpage>&#x2013;<lpage>992</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2007.06.013</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Charabi</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Al-Yahyai</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Gastli</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Evaluation of NWP performance for wind energy resource assessment in Oman</article-title>. <source>Renew. Sust. ENERG Rev.</source> <volume>15</volume> (<issue>3</issue>), <fpage>1545</fpage>&#x2013;<lpage>1555</lpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2010.11.055</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Medium-term wind power forecasting based on multi-resolution multi-learner ensemble and adaptive model selection</article-title>. <source>ENERG Convers. MANAGE</source> <volume>206</volume>, <fpage>112492</fpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2020.112492</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Pedersen</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Birgitte</surname>
<given-names>B. J.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>ARIMA-based time series model of stochastic wind power generation</article-title>. <source>IEEE T POWER Syst.</source> <volume>25</volume> (<issue>2</issue>), <fpage>667</fpage>&#x2013;<lpage>676</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRS.2009.2033277</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Lehr</surname>
<given-names>J. M.</given-names>
</name>
<name>
<surname>Lavrova</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Martinez-Ramon</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Distribution feeder-level day-ahead peak load forecasting methods and comparative study</article-title>. <source>IET Gener. Transm. Dis.</source> <volume>12</volume> (<issue>13</issue>), <fpage>3270</fpage>&#x2013;<lpage>3278</lpage>. <pub-id pub-id-type="doi">10.1049/iet-gtd.2017.1745</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Mcelroy</surname>
<given-names>M. B.</given-names>
</name>
<name>
<surname>Kang</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Integrated energy systems for higher wind penetration in China: formulation, Implementa-tion and Impacts</article-title>. <source>IEEE T POWER Syst.</source> <volume>33</volume> (<issue>2</issue>), <fpage>1</fpage>&#x2013;<lpage>1319</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRS.2017.2736943</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>A review of ultra-short-term forecasting of wind power based on data decomposition-forecasting technology combination model</article-title>. <source>ENERGY Rep.</source> <volume>8</volume>, <fpage>14200</fpage>&#x2013;<lpage>14219</lpage>. <pub-id pub-id-type="doi">10.1016/j.egyr.2022.10.342</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dai</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Wen</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>W.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Security region of renewable energy integration: characterization and flexibility</article-title>. <source>ENERGY</source> <volume>187</volume>, <fpage>115975</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2019.115975</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ding</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Nakanishi</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yokoyama</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>A gated recurrent unit neural networks based wind speed error correction model for short-term wind power forecasting</article-title>. <source>NEUROCOMPUTING</source> <volume>365</volume>, <fpage>54</fpage>&#x2013;<lpage>61</lpage>. <pub-id pub-id-type="doi">10.1016/j.neucom.2019.07.058</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Duan</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Tian</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Fang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Short-term wind power forecasting using the hybrid model of improved variational mode decomposition and Correntropy Long Short -term memory neural network</article-title>. <source>ENERGY</source> <volume>214</volume>, <fpage>118980</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2020.118980</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gan</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Temporal convolutional networks interval prediction model for wind speed forecasting</article-title>. <source>Electr. Pow. Syst. Res.</source> <volume>191</volume>, <fpage>106865</fpage>. <pub-id pub-id-type="doi">10.1016/j.epsr.2020.106865</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>He</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Tian</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>RUL prediction of wind turbine gearbox bearings based on self-calibration temporal convolutional network</article-title>. <source>IEEE T Instrum. Meas.</source> <volume>71</volume>, <fpage>1</fpage>&#x2013;<lpage>12</lpage>. <pub-id pub-id-type="doi">10.1109/TIM.2022.3143881</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hu</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Hughes</surname>
<given-names>K. J.</given-names>
</name>
<name>
<surname>Ingham</surname>
<given-names>D. B.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Pourkashanian</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Dynamic economic and emission dispatch model considering wind power under Energy Market Reform: a case study</article-title>. <source>Int. J. Electr. power &#x26; energy Syst.</source> <volume>110</volume>, <fpage>184</fpage>&#x2013;<lpage>196</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2019.03.004</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="web">
<collab>International Energy Agency</collab> (<year>2021</year>). <article-title>Wind generation increased by a record amount in 2021, but even faster growth is needed to get on the Net Zero Scenario trajectory</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.iea.org/reports/wind-electricity">https://www.iea.org/reports/wind-electricity</ext-link>.</comment>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jeon</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Taylor</surname>
<given-names>J. W.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Short-term density forecasting of wave energy using ARMA- GARCH models and kernel density estimation</article-title>. <source>Int. J. Forecast</source> <volume>32</volume> (<issue>3</issue>), <fpage>991</fpage>&#x2013;<lpage>1004</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijforecast.2015.11.003</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jiang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Short-term wind speed prediction: hybrid of ensemble empirical mode decomposition, feature selection and error correction</article-title>. <source>ENERGY Convers. MANAG.</source> <volume>144</volume>, <fpage>340</fpage>&#x2013;<lpage>350</lpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2017.04.064</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Khazaei</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Ehsan</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Soleymani</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Mohammadnezhad-Shourkaei</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>A high-accuracy hybrid method for short-term wind power forecasting</article-title>. <source>ENERGY</source> <volume>238</volume>, <fpage>122020</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2021.122020</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lahouar</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Slama</surname>
<given-names>J. B. H.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Hour-ahead wind power forecast based on random forests</article-title>. <source>Renew. ENERGY</source> <volume>109</volume>, <fpage>529</fpage>&#x2013;<lpage>541</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2017.03.064</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Ren</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Multi-step wind speed prediction based on turbulence intensity and hybrid deep neural networks</article-title>. <source>ENERG Convers. MANAGE</source> <volume>186</volume>, <fpage>306</fpage>&#x2013;<lpage>322</lpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2019.02.045</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Cen</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Tseng</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Shen</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Ali</surname>
<given-names>M. H.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Improving short-term wind power prediction using hybrid improved cuckoo search arithmetic - support vector regression machine</article-title>. <source>J. Clean. Prod.</source> <volume>279</volume>, <fpage>123739</fpage>. <pub-id pub-id-type="doi">10.1016/j.jclepro.2020.123739</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Liang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Miao</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Short-term wind power combined forecasting based on error forecast correction</article-title>. <source>ENERGY Convers. MANAG.</source> <volume>119</volume>, <fpage>215</fpage>&#x2013;<lpage>226</lpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2016.04.036</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Mei</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhen</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Jia</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Numerical weather prediction enhanced wind power forecasting: rank ensemble and probabilistic fluctuation awareness</article-title>. <source>Appl. ENERG</source> <volume>313</volume>, <fpage>118769</fpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2022.118769</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Lv</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Deterministic wind energy forecasting: a review of intelligent predictors and auxiliary methods</article-title>. <source>ENERGY Convers. MANAG.</source> <volume>195</volume>, <fpage>328</fpage>&#x2013;<lpage>345</lpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2019.05.020</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Infield</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Han</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>A hybrid forecasting method for wind power ramp based on orthogonal test and support vector machine (OT-SVM)</article-title>. <source>IEEE T SUSTAIN ENERG</source> <volume>8</volume> (<issue>2</issue>), <fpage>451</fpage>&#x2013;<lpage>457</lpage>. <pub-id pub-id-type="doi">10.1109/TSTE.2016.2604852</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Shen</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Dynamic economic emission dispatch considering renewable energy generation: a novel multi-objective optimization approach</article-title>. <source>Energy (Oxford)</source> <volume>235</volume>, <fpage>121407</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2021.121407</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Azimi</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Prediction of offshore wind farm power using a novel two-stage model combining kernel-based nonlinear extension of the Arps decline model with a multi-objective grey wolf optimizer</article-title>. <source>Renew. SUSTAIN ENERGY Rev.</source> <volume>127</volume>, <fpage>109856</fpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2020.109856</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Naik</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Satapathy</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Dash</surname>
<given-names>P. K.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Short-term wind speed and wind power prediction using hybrid empirical mode decomposition and kernel ridge regression</article-title>. <source>Appl. SOFT Comput.</source> <volume>70</volume>, <fpage>1167</fpage>&#x2013;<lpage>1188</lpage>. <pub-id pub-id-type="doi">10.1016/j.asoc.2017.12.010</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pr&#xf3;sper</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Otero-Casal</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Fern&#xe1;ndez</surname>
<given-names>F. C.</given-names>
</name>
<name>
<surname>Miguez-Macho</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Wind power forecasting for a real onshore wind farm on complex terrain using WRF high resolution simulations</article-title>. <source>Renew. ENERG</source> <volume>135</volume>, <fpage>674</fpage>&#x2013;<lpage>686</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2018.12.047</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shafiullah</surname>
<given-names>G. M.</given-names>
</name>
<name>
<surname>Oo</surname>
<given-names>A. M. T.</given-names>
</name>
<name>
<surname>Ali</surname>
<given-names>A. B. M. S.</given-names>
</name>
<name>
<surname>Wolfs</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Potential challenges of integrating large-scale wind energy into the power grid-A review</article-title>. <source>Renew. Sust. ENERG Rev.</source> <volume>20</volume>, <fpage>306</fpage>&#x2013;<lpage>321</lpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2012.11.057</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ding</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>W. J.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Hybrid forecasting model for very-short term wind power forecasting based on grey relational analysis and wind speed distribution features</article-title>. <source>IEEE T SMART GRID</source> <volume>5</volume> (<issue>1</issue>), <fpage>521</fpage>&#x2013;<lpage>526</lpage>. <pub-id pub-id-type="doi">10.1109/TSG.2013.2283269</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Qiao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>An improved random forest model of short&#x2010;term wind&#x2010;power forecasting to enhance accuracy, efficiency, and robustness</article-title>. <source>WIND ENERGY</source> <volume>21</volume> (<issue>12</issue>), <fpage>1383</fpage>&#x2013;<lpage>1394</lpage>. <pub-id pub-id-type="doi">10.1002/we.2261</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tascikaraoglu</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Uzunoglu</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>A review of combined approaches for prediction of short-term wind speed and power</article-title>. <source>Renew. Sust. ENERG Rev.</source> <volume>34</volume>, <fpage>243</fpage>&#x2013;<lpage>254</lpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2014.03.033</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tian</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Ren</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Short-term wind speed prediction based on improved PSO algorithm optimized EM-ELM</article-title>. <source>ENERG SOURCE PART A</source> <volume>41</volume>, <fpage>26</fpage>&#x2013;<lpage>46</lpage>. <pub-id pub-id-type="doi">10.1080/15567036.2018.1495782</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Torres</surname>
<given-names>J. L.</given-names>
</name>
<name>
<surname>Garc&#xed;aa</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>De Blas</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>De Francisco</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Forecast of hourly average wind speed with ARMA models in Navarre (Spain)</article-title>. <source>Sol. ENERGY</source> <volume>79</volume> (<issue>1</issue>), <fpage>65</fpage>&#x2013;<lpage>77</lpage>. <pub-id pub-id-type="doi">10.1016/j.solener.2004.09.013</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2017a</year>). <article-title>A new chaotic time series hybrid prediction method of wind power based on EEMD-SE and full-parameters continued fraction</article-title>. <source>ENERGY</source> <volume>138</volume>, <fpage>977</fpage>&#x2013;<lpage>990</lpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2017.07.112</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Wind power forecasting based on singular spectrum analysis and a new hybrid Laguerre neural network</article-title>. <source>Appl. ENERGY</source> <volume>259</volume>, <fpage>114139</fpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2019.114139</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2017b</year>). <article-title>Deep learning based ensemble approach for probabilistic wind power forecasting</article-title>. <source>Appl. ENERG</source> <volume>188</volume>, <fpage>56</fpage>&#x2013;<lpage>70</lpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2016.11.111</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2019a</year>). <article-title>A no<italic>vel non</italic>-linear combination system for short-term wind speed forecast</article-title>. <source>Renew. ENERG</source> <volume>143</volume>, <fpage>1172</fpage>&#x2013;<lpage>1192</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2019.04.154</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Niu</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2019b</year>). <article-title>A novel system for wind speed forecasting based on multi-objective optimization and Echo state network</article-title>. <source>SUSTAINABILITY-BASEL</source> <volume>11</volume> (<issue>2</issue>), <fpage>526</fpage>. <pub-id pub-id-type="doi">10.3390/su11020526</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zou</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A review of wind speed and wind power forecasting with deep neural networks</article-title>. <source>Appl. ENERG</source> <volume>304</volume>, <fpage>117766</fpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2021.117766</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Data-driven wind speed forecasting using deep feature extraction and LSTM</article-title>. <source>IET Renew. POWER Gen.</source> <volume>13</volume> (<issue>12</issue>), <fpage>2062</fpage>&#x2013;<lpage>2069</lpage>. <pub-id pub-id-type="doi">10.1049/iet-rpg.2018.5917</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>N. E.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Ensemble empirical mode decomposition: a noise-assisted data analysis method</article-title>. <source>Adv. DATA Sci. ADAPT</source> <volume>1</volume>, <fpage>1</fpage>&#x2013;<lpage>41</lpage>. <pub-id pub-id-type="doi">10.1142/s1793536909000047</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xia</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Ying</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Wind power prediction in view of ramping events based on classified spatiotemporal network</article-title>. <source>Front. ENERGY Res.</source> <volume>9</volume>, <fpage>754274</fpage>. <pub-id pub-id-type="doi">10.3389/fenrg.2021.754274</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yan</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Bai</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Deng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Foley</surname>
<given-names>A. M.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Hybrid probabilistic wind power forecasting using temporally local Gaussian process</article-title>. <source>IEEE T SUSTAIN ENERG</source> <volume>7</volume> (<issue>1</issue>), <fpage>87</fpage>&#x2013;<lpage>95</lpage>. <pub-id pub-id-type="doi">10.1109/TSTE.2015.2472963</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ye</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zeng</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Short-term wind power prediction based on spatial model</article-title>. <source>Renew. ENERG</source> <volume>101</volume>, <fpage>1067</fpage>&#x2013;<lpage>1074</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2016.09.069</pub-id>
</citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yin</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Ou</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Meng</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>A cascaded deep learning wind power prediction approach based on a two-layer of mode decomposition</article-title>. <source>ENERGY</source> <volume>189</volume>, <fpage>116316</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2019.116316</pub-id>
</citation>
</ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Bao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhai</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>A novel framework for wind speed prediction based on recurrent neural networks and support vector machine</article-title>. <source>ENERG Convers. MANAGE</source> <volume>178</volume>, <fpage>137</fpage>&#x2013;<lpage>145</lpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2018.10.008</pub-id>
</citation>
</ref>
<ref id="B56">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yuan</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Mao</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Watada</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Multi-objective wind power scenario forecasting based on PG-GAN</article-title>. <source>ENERGY</source> <volume>226</volume>, <fpage>120379</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2021.120379</pub-id>
</citation>
</ref>
<ref id="B57">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yuan</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Xin</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Cross-regional integrated transmission of wind power and pumped-storage hydropower considering the peak shaving demands of multiple power grids</article-title>. <source>Renew. ENERG</source> <volume>190</volume>, <fpage>1112</fpage>&#x2013;<lpage>1126</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2021.10.046</pub-id>
</citation>
</ref>
<ref id="B58">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Infield</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Lien</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2019a</year>). <article-title>Short-term forecasting and uncertainty analysis of wind turbine power based on long short-term memory network and Gaussian mixture model</article-title>. <source>Appl. ENERG</source> <volume>241</volume>, <fpage>229</fpage>&#x2013;<lpage>244</lpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2019.03.044</pub-id>
</citation>
</ref>
<ref id="B59">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Pan</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Wind speed prediction of IPSO-BP neural network based on lorenz disturbance</article-title>. <source>IEEE ACCESS</source> <volume>6</volume>, <fpage>53168</fpage>&#x2013;<lpage>53179</lpage>. <pub-id pub-id-type="doi">10.1109/ACCESS.2018.2869981</pub-id>
</citation>
</ref>
<ref id="B60">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Han</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Pan</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A multi-stage predicting methodology based on data decomposition and error correction for ultra-short-term wind energy prediction</article-title>. <source>J. Clean. Prod.</source> <volume>292</volume>, <fpage>125981</fpage>. <pub-id pub-id-type="doi">10.1016/j.jclepro.2021.125981</pub-id>
</citation>
</ref>
<ref id="B61">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Review on probabilistic forecasting of wind power generation</article-title>. <source>Renew. SUSTAIN ENERGY Rev.</source> <volume>32</volume>, <fpage>255</fpage>&#x2013;<lpage>270</lpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2014.01.033</pub-id>
</citation>
</ref>
<ref id="B62">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yao</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<etal/>
</person-group> (<year>2019b</year>). <article-title>Long Short-Term Memory Network based on Neighborhood Gates for processing complex causality in wind speed prediction</article-title>. <source>ENERG Convers. MANAGE</source> <volume>192</volume>, <fpage>37</fpage>&#x2013;<lpage>51</lpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2019.04.006</pub-id>
</citation>
</ref>
<ref id="B63">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Duan</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Or</surname>
<given-names>S. W.</given-names>
</name>
<name>
<surname>Chan</surname>
<given-names>K. W.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Short-term prediction of wind power and its ramp events based on semi-supervised generative adversarial network</article-title>. <source>Int. J. ELEC POWER</source> <volume>125</volume>, <fpage>106411</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2020.106411</pub-id>
</citation>
</ref>
<ref id="B64">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zucatelli</surname>
<given-names>P. J.</given-names>
</name>
<name>
<surname>Nascimento</surname>
<given-names>E. G. S.</given-names>
</name>
<name>
<surname>Santos</surname>
<given-names>A. &#xc1;. B.</given-names>
</name>
<name>
<surname>Arce</surname>
<given-names>A. M. G.</given-names>
</name>
<name>
<surname>Moreira</surname>
<given-names>D. M.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>An investigation on deep learning and wavelet transform to nowcast wind power and wind power ramp: a case study in Brazil and Uruguay</article-title>. <source>ENERGY</source> <volume>230</volume>, <fpage>120842</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2021.120842</pub-id>
</citation>
</ref>
</ref-list>
<sec id="s11">
<title>Nomenclature</title>
<table-wrap id="udT1" position="float">
<table>
<tbody valign="top">
<tr>
<td align="left">
<bold>Sets and indices</bold>
</td>
<td/>
</tr>
<tr>
<td align="left">
<bold>
<italic>K</italic>
</bold>
</td>
<td align="left">Total number of times white noise is added</td>
</tr>
<tr>
<td align="left">
<bold>
<italic>k</italic>
</bold>
</td>
<td align="left">Index of times white noise is added</td>
</tr>
<tr>
<td align="left">
<bold>
<italic>M</italic>
</bold>
</td>
<td align="left">Total number of decomposed IMFs</td>
</tr>
<tr>
<td align="left">
<bold>
<italic>m</italic>
</bold>
</td>
<td align="left">Index of decomposed IMFs</td>
</tr>
<tr>
<td align="left">
<bold>
<italic>t</italic>
</bold>
</td>
<td align="left">Index of time periods</td>
</tr>
<tr>
<td align="left">
<bold>
<italic>S</italic>
</bold>
</td>
<td align="left">Total number of sample points</td>
</tr>
<tr>
<td align="left">
<bold>
<italic>s</italic>
</bold>
</td>
<td align="left">Index of sample points</td>
</tr>
<tr>
<td align="left">
<bold>Constants</bold>
</td>
<td/>
</tr>
<tr>
<td align="left">
<inline-formula id="inf40">
<mml:math id="m58">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">P</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">Original wind power sequence</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf41">
<mml:math id="m59">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mi mathvariant="bold-italic">Z</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">Weight matrixes of the update gate</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf42">
<mml:math id="m60">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mi mathvariant="bold-italic">r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">Weight matrixes of the reset gate</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf43">
<mml:math id="m61">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">W</mml:mi>
<mml:mi mathvariant="bold-italic">h</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">Weight matrixes of the intermediate state</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf44">
<mml:math id="m62">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mi mathvariant="bold-italic">h</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">Actual historical wind power output value of sample point <italic>s</italic>
</td>
</tr>
<tr>
<td align="left">
<bold>Variables</bold>
</td>
<td/>
</tr>
<tr>
<td align="left">
<inline-formula id="inf45">
<mml:math id="m63">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">Wind power sequence after adding white noise for the <italic>k</italic>th time</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf46">
<mml:math id="m64">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">The <italic>m</italic>th IMF obtained by the EMD method for the <italic>k</italic>th time</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf47">
<mml:math id="m65">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mi mathvariant="bold-italic">m</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">The <italic>m</italic>th IMF obtained by the EEMD method</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf48">
<mml:math id="m66">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:mrow>
<mml:mi mathvariant="bold-italic">M</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">RES after EMD decomposition for the <italic>k</italic>th time</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf49">
<mml:math id="m67">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
<mml:mi mathvariant="bold-italic">M</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">RES after EEMD decomposition</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf50">
<mml:math id="m68">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">r</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">Reset gate</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf51">
<mml:math id="m69">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">White noise signal added at the <italic>k</italic>th time</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf52">
<mml:math id="m70">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">z</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">Update gate</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf53">
<mml:math id="m71">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">hidden state and load data of GRU neuron at time <italic>t</italic>
</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf54">
<mml:math id="m72">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">h</mml:mi>
<mml:mo>&#x223c;</mml:mo>
</mml:mover>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">Intermediate state</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf55">
<mml:math id="m73">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">h</mml:mi>
<mml:mi mathvariant="bold-italic">t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">Output of GRU neuron</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf56">
<mml:math id="m74">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">k</mml:mi>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mi mathvariant="bold-italic">j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">Number of times state <italic>i</italic> turns into state <italic>j</italic> after one step</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf57">
<mml:math id="m75">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mi mathvariant="bold-italic">j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">Total number of occurrences of state <italic>j</italic>
</td>
</tr>
<tr>
<td align="left">
<inline-formula id="inf58">
<mml:math id="m76">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mi mathvariant="bold-italic">f</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">Predicted value of sample point <italic>s</italic> obtained by the EEMD-GRU model</td>
</tr>
<tr>
<td align="left">
<bold>Functions</bold>
</td>
<td/>
</tr>
<tr>
<td align="left">
<bold>&#x3c3;</bold>
</td>
<td align="left">Sigmoid function</td>
</tr>
<tr>
<td align="left">
<bold>&#x2299;</bold>
</td>
<td align="left">Element-wise multiplication (Hadamard product)</td>
</tr>
<tr>
<td align="left">
<bold>Abbreviations</bold>
</td>
<td/>
</tr>
<tr>
<td align="left">
<bold>MC</bold>
</td>
<td align="left">Markov chain</td>
</tr>
<tr>
<td align="left">
<bold>GRU</bold>
</td>
<td align="left">Gated recurrent unit</td>
</tr>
<tr>
<td align="left">
<bold>EEMD</bold>
</td>
<td align="left">Ensemble empirical mode decomposition</td>
</tr>
<tr>
<td align="left">
<bold>NWP</bold>
</td>
<td align="left">Numerical weather prediction</td>
</tr>
<tr>
<td align="left">
<bold>WRF</bold>
</td>
<td align="left">Weather Research and Forecasting</td>
</tr>
<tr>
<td align="left">
<bold>ARMA</bold>
</td>
<td align="left">Autoregressive moving average model</td>
</tr>
<tr>
<td align="left">
<bold>ARIMA</bold>
</td>
<td align="left">Autoregressive integrated moving average model</td>
</tr>
<tr>
<td align="left">
<bold>GARCH</bold>
</td>
<td align="left">Generalized autoregressive conditional heteroscedasticity</td>
</tr>
<tr>
<td align="left">
<bold>SVM</bold>
</td>
<td align="left">Support vector machine</td>
</tr>
<tr>
<td align="left">
<bold>RF</bold>
</td>
<td align="left">Random forest</td>
</tr>
<tr>
<td align="left">
<bold>BART</bold>
</td>
<td align="left">Bayesian additive regression tree</td>
</tr>
<tr>
<td align="left">
<bold>AI</bold>
</td>
<td align="left">Artificial intelligence</td>
</tr>
<tr>
<td align="left">
<bold>BP</bold>
</td>
<td align="left">Back-propagation</td>
</tr>
<tr>
<td align="left">
<bold>CNN</bold>
</td>
<td align="left">Convolution neural network</td>
</tr>
<tr>
<td align="left">
<bold>RNN</bold>
</td>
<td align="left">Recursive neural network</td>
</tr>
<tr>
<td align="left">
<bold>LSTM</bold>
</td>
<td align="left">Long short-term memory</td>
</tr>
<tr>
<td align="left">
<bold>TCN</bold>
</td>
<td align="left">Temporal convolutional network</td>
</tr>
<tr>
<td align="left">
<bold>ELM</bold>
</td>
<td align="left">Extreme learning machine</td>
</tr>
<tr>
<td align="left">
<bold>ENN</bold>
</td>
<td align="left">Elman neural network</td>
</tr>
<tr>
<td align="left">
<bold>LSSVM</bold>
</td>
<td align="left">Least square SVM</td>
</tr>
<tr>
<td align="left">
<bold>RBFNN</bold>
</td>
<td align="left">Radial basis function neural network</td>
</tr>
<tr>
<td align="left">
<bold>ORELM</bold>
</td>
<td align="left">Outlier robust ELM</td>
</tr>
<tr>
<td align="left">
<bold>BiLSTM</bold>
</td>
<td align="left">Bidirectional LSTM</td>
</tr>
<tr>
<td align="left">
<bold>EMD</bold>
</td>
<td align="left">Empirical mode decomposition</td>
</tr>
<tr>
<td align="left">
<bold>VMD</bold>
</td>
<td align="left">Variational mode decomposition</td>
</tr>
<tr>
<td align="left">
<bold>WT</bold>
</td>
<td align="left">Wavelet transform</td>
</tr>
<tr>
<td align="left">
<bold>CEEMDAN</bold>
</td>
<td align="left">Complete ensemble empirical mode decomposition with adaptive noise</td>
</tr>
<tr>
<td align="left">
<bold>MLP</bold>
</td>
<td align="left">Multilayer Perceptron</td>
</tr>
<tr>
<td align="left">
<bold>NSGA- &#x4c0;&#x4c0;</bold>
</td>
<td align="left">Non-dominated Sorting Genetic Algorithm II</td>
</tr>
<tr>
<td align="left">
<bold>MOGWO</bold>
</td>
<td align="left">Multi-objective grey wolf optimization</td>
</tr>
<tr>
<td align="left">
<bold>PSO</bold>
</td>
<td align="left">Particle Swarm Optimization</td>
</tr>
<tr>
<td align="left">
<bold>ESN</bold>
</td>
<td align="left">Echo State Network</td>
</tr>
<tr>
<td align="left">
<bold>IMF</bold>
</td>
<td align="left">Intrinsic mode function</td>
</tr>
<tr>
<td align="left">
<bold>RES</bold>
</td>
<td align="left">Residual</td>
</tr>
<tr>
<td align="left">
<bold>GRUNN</bold>
</td>
<td align="left">GRU neural network</td>
</tr>
<tr>
<td align="left">
<bold>RMSE</bold>
</td>
<td align="left">Root mean square error</td>
</tr>
<tr>
<td align="left">
<bold>MAE</bold>
</td>
<td align="left">Mean absolute error</td>
</tr>
<tr>
<td align="left">
<bold>MAPE</bold>
</td>
<td align="left">Mean absolute percentage error</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</back>
</article>