矿井 瓦斯浓度 时序往往同时含 线性趋势 与 非线性波动,单一 ARIMA 或 LSTM 难以兼顾。Li 等(2023)在 Processes 提出 ARIMA-LSTM 组合模型:ARIMA 给出线性预报与残差,LSTM 学习残差中的非线性,再叠加得到最终预测。下文按「背景—数据—方法—结果—局限—工程启示」归纳。
1. 问题背景
1.1 单模型的不足
| 模型 | 优势 | 在瓦斯序列上的局限 | |------|------|-------------------| | ARIMA | 线性结构、可解释 | 非线性模式捕捉弱 | | LSTM | 非线性、长依赖 | 线性成分利用不足 |
既有 SVM、蚁群、随机森林、纯 LSTM 等工作提升了精度,但较少 显式分解线性/非线性。
1.2 组合思路
将序列分解为 y = 线性 + 非线性:
- ARIMA → 线性预报
l与残差δ = y − l - LSTM → 对残差建模得非线性项
nl - 叠加 →
ŷ = l + nl
2. 数据与预处理
- 测点: 某矿 803 工作面上隅角
- 时段: 2021 年 3 月 1–6 日;30 min 采样,共 288 条
- 划分: 训练 3 月 1–4 日(192 条);测试 3 月 5–6 日(96 条)
- 清洗: 3σ(拉依达)准则 剔除粗差;异常值以邻点均值替代
3. 方法流程
原始序列 → 3σ 清洗 → ADF + 差分 → ARIMA(3,1,0) → 残差 → LSTM → 叠加 → 预报
3.1 ARIMA 定阶
- ADF 检验: 原序列非平稳 → 一阶差分
- ACF/PACF + BIC → 最优 ARIMA (3, 1, 0)
- 残差 Ljung–Box(p ≈ 0.55)、Q–Q 图支持白噪声假设
3.2 LSTM 拟合残差
- 输入为 ARIMA 残差(非线性成分)
- ReLU 激活;batch_size = 1;最多 100 轮;损失约 第 30 轮 收敛
- 数据经归一化后送入网络
4. 结果(原文 Table 2,测试集)
| 模型 | R² | MAPE | RMSE | |------|-----|------|------| | 仅 ARIMA | 0.3648 | 1.4135 | 1.5769 | | 仅 LSTM | 0.5244 | 0.4253 | 0.7823 | | ARIMA-LSTM | 0.9825 | 0.0124 | 0.0830 |
组合模型在报告指标上明显优于两种单模型。作者指出,误差仍可能来自 通风、采速等未建模外生因素。
5. 局限与复现注意
- 样本短: 单点上隅角 6 天——部署前宜在多工作面、更长时段验证。
- 单变量: 尚未纳入外生协变量;作者计划扩展多因素预测。
- 指标量纲: RMSE/MAPE 与浓度尺度相关,不宜当作通用验收门槛。
- 高 R² 风险: 需排查泄漏、小测试集过拟合,并做 时间外推 检验。
- 安全链路: 预报服务预警,不能替代甲烷传感器、通风联锁与规程限值。
6. 工程启示
- 瓦斯浓度序列可优先尝试 线性模型 + 残差深度网络,再对比端到端黑箱。
- 保留 ADF、ACF/PACF、Ljung–Box 等 ARIMA 诊断,便于审计与运维。
- 与 同一测点阈值/报警策略 对接,形成「预测 → 告警 → 人工确认」闭环。
- 下一步可引入 风量、推进速度、抽采量 等外生变量,与本系列 RFECV-Bi-LSTM 涌出量工作形成对照。
参考文献
Li, C.; Fang, X.; Yan, Z.; Huang, Y.; Liang, M. Research on Gas Concentration Prediction Based on the ARIMA-LSTM Combination Model. Processes 2023, 11 (1), 174.