融合离散元模拟与集成学习方法的
开采沉陷时间序列预测模型

研究背景与技术挑战

地下矿产资源的大规模开采会破坏原岩应力平衡,诱发覆岩垮落并导致地表沉陷。当前的沉陷防控主要依赖实时监测,但监测数据仅能反映已发生的沉陷,难以精准预测未来趋势。

本研究结合离散元模拟(DEM)的数据扩展优势与集成学习的数据驱动优势,提出LSTM-TST模型,实现采动区域地表沉陷监测与预测的一体化。

现有方法的四大局限

理论模型局限

概率积分法、典型曲线法在处理复杂地质条件和非线性变形时具有局限性,难以精确预测。

监测技术瓶颈

传统点状监测站受地形影响大、工作量繁重,且获取的数据在空间和时间上具有离散性。

深度学习短板

LSTM在长序列预测中易出现数据挖掘不充分,Transformer对短期动态变化的敏感性不足。

数据瓶颈

深度学习模型的性能高度依赖于大规模、高质量的数据集,而现场实测数据往往有限。

技术路线与方法创新

1

参数化建模(Rhino-Grasshopper)

利用Rhino-Grasshopper平台开发自定义电池组,实现综采工作面的高精度参数化建模,具有极高的灵活性和自动化批处理能力。

输入钻孔柱状图数据 → 自动生成三维地质体 → 快速适应不同地质条件

2

离散元模拟与自动监测(3DEC)

基于3DEC开展模拟,编写FISH语言程序实现监测点的自动布设,在不同推进距离下实时提取沉陷数据。

数据生成能力

  • • 模型规模:185,773个块体单元
  • • 开挖步数:60步(5m/步)
  • • 每个测点:240组数据
3

集成学习模型(LSTM-TST)

LSTM提取短期波动和局部时序模式,TST(时序Transformer)捕获长程依赖关系,通过集成实现优势互补。

优化机制

  • • 正则化准则调整权重
  • • 早停机制(Early Stopping)
  • • 粒子群算法(PSO)超参数寻优

工程案例:冯家塔煤矿1208工作面

地质与开采参数

煤层特征

  • 山西组2#煤层,平均厚度3.17m
  • 倾角约4°(近水平),平均埋深127m
  • 表土层厚度达78.42m

工作面尺寸

  • 走向长度:240m
  • 倾向长度:216.9m
  • 采高:3m

模型验证结果

初次来压步距

模拟值:59.60m

实测值:60.32m

周期来压步距

模拟值:17~20m

实测值:18~22m

模型性能突破

0.989
R²决定系数
倾向测点预测精度
0.0008
RMSE均方根误差
显著优于单一模型
2.06%
R²提升幅度
相比单一TST模型

LSTM-TST集成模型在倾向和走向测点上的R²均超过0.98,RMSE控制在0.001以内,实现了高精度的沉陷预测。

跨矿区验证结果

验证案例1:冯家塔煤矿1209工作面

走向方向

监测点:16个(269m)

最大误差:69.2mm

倾向方向

监测点:21个(374m)

最大误差:47.0mm

验证案例2:榆家梁煤矿52309工作面

走向方向

监测点:28个(432m)

最大误差:81.2mm

倾向方向

监测点:30个(468m)

最大误差:74.5mm

泛化能力分析

模型在不同矿区均表现出较高的预测精度。榆家梁煤矿的误差略高,这是由于地质条件与开采工艺的差异所致。通过迁移学习(Transfer Learning)对模型进行微调可进一步提高精度。

研究价值与应用前景

方法论创新

成功融合了参数化建模、离散元模拟数据增强与集成学习技术,实现了地表沉陷"监测-建模-预测"的一体化。

数据驱动优势

利用3DEC自动监测FISH程序生成的高频率时序数据集,有效解决了深度学习在开采沉陷领域训练样本不足的问题。

工程实用性

通过两个不同工作面的实测数据验证,证明模型具有良好的泛化能力和较高的工程应用价值,可为矿区沉陷防控提供科学依据。