时间序列预测实战：数据清洗与集成策略-编程实验室

1. 时间序列预测的挑战与机遇

刚接手时间序列预测项目时，我常常被两个问题困扰：为什么同样的算法在不同数据集上表现差异巨大？为什么论文里的漂亮指标在实际业务中总是不尽如人意？经过多年实战，我发现数据预处理的质量和预测策略的设计才是决定成败的关键因素。超市销售额预测、电力负荷预测、股票价格波动这些看似迥异的场景，本质上都在处理带时间戳的连续数据流，而优秀的预处理和集成策略能让预测效果提升30%以上。

时间序列数据就像一位固执的谈判对手，它带着季节性、趋势性、噪声和各种异常值与你周旋。我曾见过某电商平台的GMV预测因为忽视春节效应导致备货严重失误，也经历过医疗设备故障预测因未处理传感器漂移而产生大量误报。这些教训让我意识到：预测模型只是工具，真正的高手功夫都在数据准备阶段。

2. 时间序列数据深度清洗实战

2.1 异常值检测与修复方案

在电力负荷预测项目中，我总结出异常值处理的"三级防御体系"：

统计阈值法：用滑动窗口计算均值±3σ范围，标记超出阈值的点

def sigma_filter(data, window=24, n_sigma=3): rolling_mean = data.rolling(window).mean() rolling_std = data.rolling(window).std() upper = rolling_mean + n_sigma * rolling_std lower = rolling_mean - n_sigma * rolling_std return np.where((data > upper) | (data < lower), np.nan, data)

孤立森林检测：对具有多重周期性的数据（如同时含日周期和周周期），传统统计方法容易误判，这时需要无监督学习算法
业务规则校验：与领域专家确认异常区间，比如电力数据中零值持续超过2小时明显不符合常理

特别注意：不要简单删除异常值！我常用线性插值处理短暂缺失，对于连续异常段则采用历史同期数据加权平均填充。某次处理风电功率数据时，盲目删除"异常低值"导致低估了台风天气的影响，教训深刻。

2.2 趋势与季节性分解技巧

STL分解（Seasonal-Trend decomposition using Loess）是我的首选工具，相比传统的X11或SEATS方法，它对异常值更鲁棒。在零售销售预测中，我这样调整参数：

from statsmodels.tsa.seasonal import STL result = STL(sales_data, period=7, seasonal=13, trend=21).fit()

seasonal参数应小于周期长度（这里是周数据所以period=7）
trend窗口建议取周期长度的3倍左右
对明显非线性趋势，先用对数变换压缩动态范围

实践中发现，很多教程忽略的残差项其实包含宝贵信息。某次设备故障预测中，正是残差中的特定模式帮助我们提前2周发现了轴承磨损问题。

3. 特征工程黄金法则

3.1 时序特征构造

除了常规的滞后项(lagged features)，这些特征常被忽视但效果显著：

周期位置特征：将时间戳映射到正弦/余弦空间，比独热编码更能保持周期性
```
hour_sin = np.sin(2 * np.pi * df['hour']/24) hour_cos = np.cos(2 * np.pi * df['hour']/24)
```
变化速率特征：一阶差分往往不够，可以添加二阶差分、百分比变化等
事件窗口统计量：最近3/7/30天的均值、方差、偏度等

3.2 外部特征融合

在预测餐厅客流量时，我发现这些外部特征至关重要：

天气数据（温度、降水量、风速）
节假日标记（包括节前节后特殊日期）
周边竞品促销信息（通过爬虫获取）
宏观经济指标（适用于中长期预测）

血泪教训：外部特征一定要检查数据延迟！某次用天气预报数据训练模型，却在实际预测时误用了实况数据，导致线上线下效果严重不一致。

4. 集成预测策略精要

4.1 分层建模架构

我的标准工作流包含三个层级：

基础模型层：ARIMA、ETS、Prophet等传统时序模型
机器学习层：LightGBM、XGBoost、随机森林
深度学习层：N-BEATS、Temporal Fusion Transformer

关键技巧是让不同模型聚焦不同成分：

ARIMA捕捉线性关系
LightGBM处理外部特征
N-BEATS学习复杂非线性模式

4.2 动态权重集成方案

静态加权平均（如按验证集表现分配权重）效果有限，我开发了两种动态策略：

残差感知加权：

def dynamic_weight(models, X, y): recent_errors = [] for model in models: pred = model.predict(X[-30:]) # 最近30个样本 error = mean_absolute_error(y[-30:], pred) recent_errors.append(error) weights = softmax(-np.array(recent_errors)) return weights

概率集成：用分位数回归生成预测区间，当各模型预测差异较大时自动降低置信度，触发人工复核机制。这套方案在金融风控场景中帮助我们将误判率降低了42%。

5. 生产环境部署陷阱

5.1 实时预测的冷启动问题

新上线模型常遇到历史数据不足的情况，我的解决方案是：

构建"数据生成器"模块，用类似场景数据做迁移学习
设计渐进式置信度机制，初期更多依赖业务规则
实现自动回填功能，用新数据持续更新特征库

5.2 概念漂移检测

某电商平台的预测模型上线3个月后效果逐渐变差，后来我们实现了滑动窗口KS检验：

from scipy import stats def detect_drift(new_data, baseline, window=100): p_values = [] for i in range(0, len(new_data), window): _, p = stats.ks_2samp(baseline, new_data[i:i+window]) p_values.append(p) return np.mean(p_values) < 0.01

当检测到漂移时自动触发模型再训练，并将新旧模型预测结果做对比分析。