Tab. 1
A part of the HDD data from 1974 to 2017.
| Jan | Feb | Mar | Apr | May | Jun | Jul | Aug | Sep | Oct | Nov | Dec | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 1974 | 434.2 | 385.8 | 327.9 | 248.8 | 128.5 | 24.71 | 4.658 | 2.940 | 30.13 | 150.3 | 297.2 | 398.4 |
| 1975 | 368.1 | 345.4 | 386.7 | 270.5 | 186.6 | 77.73 | 17.10 | 11.21 | 40.27 | 215.6 | 331.3 | 469.7 |
| . | ||||||||||||
| . | ||||||||||||
| . | ||||||||||||
| 2015 | 414.6 | 402.0 | 320.5 | 200.5 | 117.4 | 31.28 | 9.428 | 8.295 | 77.75 | 184.5 | 222.4 | 279.9 |
| 2016 | 363.3 | 341.3 | 347.9 | 254.1 | 135.3 | 41.40 | 12.62 | 8.048 | 21.08 | 180.9 | 289.1 | 390.3 |
| 2017 | 486.8 | 321.9 | 256.0 | 236.1 | 109.4 | 24.20 | 6.902 | 10.01 | 71.63 | 113.2 | 300.3 | 400.1 |

Fig. 1
Time series graph of the HHD series.

Fig. 2
The correlogram of the HDD series for k=36 month delays.
Tab. 2
Monthly rate averages and seasonal indices.
| Jan | Feb | Mar | Apr | May | Jun | Jul | Aug | Sep | Oct | Nov | Dec | Total | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Arithmetic mean of ratios | 2.09 | 1.78 | 1.54 | 1.13 | 0.59 | 0.21 | 0.06 | 0.06 | 0.27 | 0.79 | 1.50 | 1.97 | 11.99 |
| Seasonal index | 2.09 | 1.78 | 1.54 | 1.13 | 0,59 | 0.21 | 0.06 | 0.06 | 0.27 | 0.79 | 1.50 | 1.97 | 12.00 |
| Seasonal index(%) | 209 | 178 | 154 | 113 | 59 | 21 | 6 | 6 | 27 | 79 | 150 | 197 | 1200 |

Fig. 3
The plot of actual and seasonal adjusted HDD data sets.

Fig. 4
The plot of seasonally adjusted data and trend line.

Fig. 5
The graph of trend corrected data.

Fig. 6
Distribution of cyclic indices.

Fig. 7
Actual and forecasted values obtained by the decomposition method.
Tab. 3
Tentative models generated in this study.
| Models | |||
|---|---|---|---|
| *SARIMA (1,0,1)(1,0,1)12 | SARIMA (1,0,0)(1,1,1)12 | SARIMA (1,1,1)(1,1,0)12 | *SARIMA (1,0,2)(0,1,1)12 |
| SARIMA (0,0,1)(1,0,1)12 | *SARIMA (1,0,1)(0,1,1)12 | SARIMA (0,1,0)(1,1,1)12 | SARIMA (1,0,1)(0,1,2)12 |
| SARIMA (1,0,0)(1,0,1)12 | SARIMA (1,0,1)(1,1,0)12 | SARIMA (1,1,0)(1,1,0)12 | *SARIMA (0,0,2)(0,1,1)12 |
| SARIMA (1,0,1)(0,0,1)12 | SARIMA (1,0,0)(1,1,0)12 | SARIMA (0,1,1)(1,1,0)12 | SARIMA (0,0,1)(0,1,2)12 |
| SARIMA (1,0,1)(1,0,0)12 | SARIMA (0,0,1)(1,1,0)12 | *SARIMA (0,1,1)(0,1,1)12 | SARIMA (2,0,0)(0,1,1)12 |
| SARIMA (1,0,0)(1,0,0)12 | *SARIMA (0,0,1)(0,1,1)12 | SARIMA (1,1,0)(0,1,1)12 | SARIMA (1,0,0)(0,1,2)12 |
| SARIMA (0,0,1)(1,0,0)12 | *SARIMA (1,0,0)(0,1,1)12 | *SARIMA (2,0,1)(1,0,1)12 | SARIMA (2,1,1)(0,1,1)12 |
| SARIMA (1,0,0)(0,0,1)12 | SARIMA (1,1,1)(1,1,1)12 | SARIMA (1,0,2)(1,0,1)12 | SARIMA (1,1,2)(0,1,1)12 |
| SARIMA (0,0,1)(0,0,1)12 | SARIMA (0,1,1)(1,1,1)12 | SARIMA (1,0,1)(2,0,1)12 | SARIMA (1,1,1)(0,1,2)12 |
| SARIMA (1,0,1)(1,1,1)12 | SARIMA (1,1,0)(1,1,1)12 | SARIMA (1,0,1)(1,0,2)12 | SARIMA (0,1,2)(0,1,1)12 |
| SARIMA (0,0,1)(1,1,1)12 | *SARIMA (1,1,1)(0,1,1)12 | SARIMA (2,0,1)(0,1,1)12 | SARIMA (0,1,1)(0,1,2)12 |
Tab. 4
The selected models in the first step and second step selection criteria.
| Model | Adjusted R2 | Residual sum of squares | AIC | SIC | |
|---|---|---|---|---|---|
| 1 | SARIMA (1,0,1)(1,0,1)12 | 0.936654 | 828032.9 | 10.33716 | 10.37759 |
| 2 | SARIMA (1,0,1)(0,1,1)12 | 0.475825 | 812361.4 | 10.26529 | 10.29820 |
| 3 | SARIMA (0,0,1)(0,1,1)12 | 0.464422 | 831655.0 | 10.27557 | 10.30026 |
| 4 | SARIMA (1,0,0)(0,1,1)12 | 0.615053 | 1090103.0 | 10.57437 | 10.59910 |
| 5 | SARIMA (1,1,1)(0,1,1)12 | 0.721348 | 787551.5 | 10.27085 | 10.30382 |
| 6 | SARIMA (0,1,1)(0,1,1)12 | 0.715896 | 804533.5 | 10.28096 | 10.30569 |
| 7 | SARIMA (2,0,1)(1,0,1)12 | 0.936911 | 816848.7 | 10.32561 | 10.37412 |
| 8 | SARIMA (1,0,2)(0,1,1)12 | 0.487184 | 793205.5 | 10.26017 | 10.30132 |
| 9 | SARIMA (0,0,2)(0,1,1)12 | 0.467373 | 825461.2 | 10.27325 | 10.30617 |
Tab. 5
The actual and forecasted HDD values for 2017.
| Months | Actual (kWh) | Forecasted (°C*day) | Forecasted (°C*day) |
|---|---|---|---|
| 2017M01 | 497.52 | 490.2 | 417.986 |
| 2017M02 | 303.62 | 315.16 | 361.840 |
| 2017M03 | 248.77 | 251.1 | 308.607 |
| 2017M04 | 233.89 | 228.96 | 217.635 |
| 2017M05 | 95.36 | 97.75 | 111.711 |
| 2017M06 | 14.73 | 16.71 | 32.316 |
| 2017M07 | 4.98 | 4.95 | 4.335 |
| 2017M08 | 10.36 | 9.84 | 5.148 |
| 2017M09 | 78.48 | 75 | 48.420 |
| 2017M010 | 117.86 | 127.85 | 151.317 |
| 2017M011 | 321.16 | 313.33 | 294.122 |
| 2017M012 | 411.22 | 411.29 | 395.733 |

Fig. 8
Plots of the HDD series’ residuals according to the selected model.

Fig. 9
ACF plot of residuals.

Fig. 10
PACF plot of residuals.

Fig. 11
The actual values, the forecasted values, and the confidence intervals of the HDD data according to model SARIMA (2,0,1) (1,0,1) 12.