Forecasting Voltage Time Series under Anomalies Using a WGAN Model With Lipschitz Control and Progressive Kernel Critic Strategy
Organizations: Institut für Maschinelle Sprachverarbeitung, University Of Stuttgart, Stuttgart, Germany
Abstract
In time-series forecasting tasks, the use of recursive forecasting on high-dimensional data substantially increases the complexity of the prediction process. The main challenge is that long-term dynamics becomes harder to model, since the sequences show greater temporal variability. To counteract this limitation and improve the modeling of long-term dynamics, and more importantly, address high-dimensional recursive forecasting, this paper uses a Wasserstein GAN equipped with a progressive Lipschitz penalty and an adaptive kernel. This model is capable of taking a window of the data stream and forecasting the next values, subsequently reconstructing the time-series stream and dynamically adapting its kernel to each new window. The proposed model is a deep learning framework based on a Wasserstein GAN (WGAN-CLP) with an adaptive kernel to model alterations in streaming voltage time series. To establish the methodological foundation, we review the WGAN-LP formulation another variant model in this class and extend its framework through an adaptive kernel critic model. We then evaluate this extended model using power ratios, correlation coefficient (CC), and relative root mean squared error (RRMSE). The experimental results demonstrate that the proposed model outperforms benchmark methods in both accuracy and robustness. The evaluation shows promising performance of the WGAN-CLP with an adaptive kernel framework in forecasting sequence alterations across the voltage and energy time-series datasets.
Figures & tables
| Model | Year | Distinguishing feature |
|---|---|---|
| WGAN-DIV | 2018 | Uses Wasserstein divergence instead of an explicit Lipschitz constraint [ 57 ] . |
| CWGAN | 2017 | Conditions the generator and critic on auxiliary input information [ 43 , 2 ] . |
| WGAN-QC | 2019 | Uses quadratic transport cost and an optimal-transport regularizer [ 41 ] . |
| WGAN-LP | 2018 | Applies a one-sided gradient penalty for the Lipschitz constraint [ 47 ] . |
| WGAN-CT | 2018 | Adds a consistency term penalizing Lipschitz violations between critic inputs [ 53 ] . |
| Dataset | Rows | Features | Anomalies | Anom. Rate | Quality | Key Issue |
|---|---|---|---|---|---|---|
| Millisecond-PMU (Forced) | 5.2K | 156 | 2,156 | 41.5% | 76.2% | Transient spikes |
| Millisecond-PMU (Natural) | 6.8K | 156 | 2,131 | 31.3% | 80.4% | Frequency drift |
| Minute-Load-Renewable (CAISO) | 3.1K | 89 | 742 | 23.9% | 84.1% | Solar ramps |
| Minute-Load-Renewable (ERCOT) | 2.4K | 89 | 621 | 25.9% | 82.5% | Wind variability |
| Minute-Load-Renewable (NYISO) | 2.5K | 89 | 560 | 22.4% | 88.3% | Load forecasting |
| Minute-PMU (Zone 1-5) | 4.1K | 124 | 1,089 | 26.5% | 81.2% | Voltage sag |
| Anomaly Type | Count | % of Total | Top Var. | Severity | Pattern | Mitigation |
|---|---|---|---|---|---|---|
| Voltage transients | 2,847 | 34.0% | VOLT_BUS_650 | High | Spike | Filter |
| Power oscillations | 1,923 | 23.0% | POWR_SOURCEBUS | High | Ring | Damping |
| Frequency deviations | 1,056 | 12.6% | FREQ_SYS | Medium | Drift | Regulation |
| Reactive power swings | 894 | 10.7% | VARS_RENEWABLE | Medium | Oscillate | Compensation |
| Load forecast error | 742 | 8.9% | LOAD_CAISO | Low | Bias | Recalibrate |
| Sensor saturation | 562 | 6.7% | VOLT_ZONE | High | Clipping | Calibrate |
| Model | Gen. Loss | Critic Loss | Div. | Stability | Epochs | Time (min) |
|---|---|---|---|---|---|---|
| M1a: Baseline (lr=1e-3) | 0.512 | -0.678 | High | Poor | 100 | 11.2 |
| M1b: Baseline (lr=1e-4) | 0.342 | -0.456 | Medium | Fair | 100 | 12.3 |
| M2a: WGAN-WC (clip=0.01) | 0.267 | -0.234 | Low | Good | 100 | 14.1 |
| M2b: WGAN-WC (clip=0.001) | 0.215 | -0.189 | Low | Good | 100 | 15.7 |
| M3a: WGAN-GP ( =5) | 0.198 | -0.145 | Very Low | V. Good | 100 | 16.8 |
| M3b: WGAN-GP ( =10) | 0.156 | -0.098 | Very Low | V. Good | 100 | 18.2 |
| Model | MSE | KL Div. | Wasserstein | FID | Inception | JS Div. |
|---|---|---|---|---|---|---|
| M1a (Baseline, high-lr) | 0.0421 | 0.567 | 0.234 | 28.34 | 0.612 | 0.289 |
| M1b (Baseline, low-lr) | 0.0234 | 0.342 | 0.127 | 18.45 | 0.456 | 0.178 |
| M2a (WGAN-WC, clip=0.01) | 0.0212 | 0.298 | 0.112 | 16.23 | 0.398 | 0.156 |
| M2b (WGAN-WC, clip=0.001) | 0.0198 | 0.281 | 0.098 | 14.23 | 0.345 | 0.134 |
| M3a (WGAN-GP, =5) | 0.0167 | 0.215 | 0.078 | 11.12 | 0.287 | 0.098 |
| M3b (WGAN-GP, =10) | 0.0156 | 0.198 | 0.076 | 9.87 | 0.267 | 0.087 |
| Config | Architecture | Regularization | Batch Size | Opt. | Init. Scale |
|---|---|---|---|---|---|
| M1a | 256 128 | None | 64 | Adam (1e-3) | Glorot |
| M1b | 256 128 | None | 64 | Adam (1e-4) | Glorot |
| M2a | 256 128 | WC (0.01) | 64 | RMSprop (1e-4) | Glorot |
| M2b | 256 128 | WC (0.001) | 64 | RMSprop (1e-4) | Glorot |
| M3a | 256 128 | GP ( =5) | 64 | Adam (1e-4) | Glorot |
| M3b | 256 128 | GP ( =10) | 64 | Adam (1e-4) | Glorot |
| Model | Mode Coverage | Throughput | GPU Mem | Peak RAM | Convergence | Cost/Sample |
|---|---|---|---|---|---|---|
| M1a | 72% | 9,234 s/s | 0.98 GB | 1.2 GB | Slow | $0.00012 |
| M1b | 78% | 8,652 s/s | 1.0 GB | 1.3 GB | Medium | $0.00014 |
| M2a | 80% | 7,891 s/s | 1.15 GB | 1.45 GB | Fast | $0.00018 |
| M2b | 82% | 6,802 s/s | 1.2 GB | 1.52 GB | Fast | $0.00023 |
| M3a | 86% | 5,934 s/s | 1.38 GB | 1.68 GB | V. Fast | $0.00028 |
| M3b | 87% | 5,867 s/s | 1.4 GB | 1.71 GB | V. Fast | $0.00031 |
| Feature | Variance | Anom. Freq | Correlation | Sensitivity | Impact |
|---|---|---|---|---|---|
| VOLT_BUS_650 | 0.342 | 2,847 | 0.892 | High | Critical |
| POWR_SOURCEBUS | 0.287 | 1,923 | 0.856 | High | Critical |
| FREQ_SYS | 0.156 | 1,056 | 0.743 | Medium | Important |
| VARS_RENEWABLE | 0.098 | 894 | 0.621 | Medium | Important |
| LOAD_CAISO | 0.089 | 742 | 0.534 | Low | Moderate |
| VOLT_ZONE | 0.067 | 562 | 0.478 | Low | Moderate |
| Model | Clean Data | Anom. Data | Delta | Robustness | Bias |
|---|---|---|---|---|---|
| M1b (Baseline) | 0.0234 | 0.0521 | 0.0287 | 55.0% | +2.23pp |
| M2b (WGAN-WC) | 0.0198 | 0.0412 | 0.0214 | 67.8% | +1.89pp |
| M3b (WGAN-GP) | 0.0156 | 0.0289 | 0.0133 | 81.8% | +0.92pp |
| M4b (AK-LP) | 0.0089 | 0.0156 | 0.0067 | 94.2% | +0.34pp |
| Statistic | Real Data | M1b Synth. | M3b Synth. | M4b Synth. | Fidelity (M4b) |
|---|---|---|---|---|---|
| Mean | 2.341 | 2.156 | 2.298 | 2.334 | 99.7% |
| Std Dev | 1.287 | 0.987 | 1.156 | 1.267 | 98.4% |
| Skewness | 0.456 | 0.234 | 0.401 | 0.448 | 98.2% |
| Kurtosis | 3.214 | 2.456 | 2.987 | 3.189 | 99.2% |
| Min | -3.456 | -2.123 | -3.201 | -3.412 | 98.7% |
| Max | 9.876 | 8.234 | 9.543 | 9.812 | 99.3% |
| Epoch | M1b Loss | M2b Loss | M3b Loss | M4b Loss | Gradient Norm | Step Size |
|---|---|---|---|---|---|---|
| 10 | 0.567 | 0.421 | 0.356 | 0.287 | 2.341 | 1e-4 |
| 20 | 0.478 | 0.334 | 0.267 | 0.198 | 1.876 | 1e-4 |
| 30 | 0.423 | 0.298 | 0.218 | 0.134 | 1.234 | 1e-4 |
| 50 | 0.345 | 0.245 | 0.167 | 0.098 | 0.876 | 1e-4 |
| 75 | 0.287 | 0.203 | 0.121 | 0.067 | 0.456 | 1e-4 |
| 100 | 0.234 | 0.198 | 0.089 | 0.045 | 0.234 | 1e-4 |
| Train Set | Test Set | MSE | KL Div | Domain Gap | Generalization |
|---|---|---|---|---|---|
| Millisecond-PMU | Millisecond-PMU | 0.0089 | 0.087 | 0.0% | Excellent |
| Millisecond-PMU | Minute-Load | 0.0156 | 0.234 | 11.2% | Good |
| Millisecond-PMU | Minute-PMU | 0.0198 | 0.301 | 24.5% | Fair |
| Minute-Load | Minute-Load | 0.0091 | 0.095 | 0.0% | Excellent |
| Minute-Load | Millisecond-PMU | 0.0267 | 0.456 | 34.8% | Poor |
| Minute-Load | Minute-PMU | 0.0145 | 0.189 | 12.3% | Good |
| Use Case | Baseline | M4b | Effi. Gain | ROI | Payback Period |
|---|---|---|---|---|---|
| Data augmentation (100K samples) | $4,200 | $1,890 | 55.0% | +$2,310 | 2.1 months |
| Anomaly detector training | $6,800 | $2,450 | 63.9% | +$4,350 | 1.8 months |
| Real-time monitoring (1M samples/day) | $18,500 | $7,234 | 60.9% | +$11,266 | 1.3 months |
| Forecasting model development | $9,200 | $3,120 | 66.1% | +$6,080 | 1.9 months |
| Scenario analysis (50K scenarios) | $8,900 | $3,670 | 58.8% | +$5,230 | 2.4 months |
| Grid resilience testing | $12,300 | $4,980 | 59.5% | +$7,320 | 2.0 months |
Time Series Forecasting
Time series forecasting predicts future values from ordered observations using statistical, neural, and probabilistic models to support decisions under uncertainty.
See all 839 papers