一种基于元胞自动机的地表下沉仿真方法及系统技术方案

本申请公开了一种基于元胞自动机的地表下沉仿真方法及系统，该方法包括：获取地面空间数据以及开采工作面布置情形和采矿地质数据；将地面空间数据进行空间插值，形成预设分辨率的格网DEM面域数据；对格网DEM面域数据以及开采工作面布置情形、采矿地质数据进行分析以确定开采下沉影响区域，构建元胞自动机模型；利用构建的元胞自动机模型对地表下沉过程进行仿真。该方法能够仿真较为复杂地质采矿条件下的矿产开采引起的下沉过程，动态模拟开采沉陷的开采影响范围、下沉活跃的强度和面积等，并且得到接近真实情况的地表下沉结果，从而在很大程度上提高了对复杂地质采矿条件下矿产开采时的地表沉陷预测的准确性。

【技术实现步骤摘要】
一种基于元胞自动机的地表下沉仿真方法及系统
本申请涉及工业监测
，更具体地说，涉及一种基于元胞自动机的地表下沉仿真方法及系统。
技术介绍
大规模矿产资源的开发和利用，给人类带来巨大的经济和社会效益的同时，也对人类生存的环境带来了巨大的破坏影响。由于地下矿产资源在开采过程中，开采矿层的上覆岩层将出现位移和变形，引起的开采沉陷可导致一系列环境问题，甚至引发重大的地质灾害事故，造成巨大的经济损失及严重的社会问题。据不完全统计，我国自20世纪90年代以来，因煤炭开采形成的地表塌陷每年约2.2万hm2，到本世纪初，我国的采矿业每年占用和破坏的土地达3.4万hm2。而开采沉陷涉及到地质、采矿等多种因素，是一种复杂的三维、时空动态过程，地表移动是岩层移动传播到地表的沉陷现象，反映了岩层移动的传播方式和移动状况。因此，在开采前就对可能引起的开采沉陷状况进行一个合理的预测，并在开采过程中，实时对引起的沉陷过程进行监测，这样才便于对安全生产进行指导。目前，对矿区的地质分析以及沉陷过程的监测仅是借鉴地面观测站的监测与预计研究成果。由于其利用的数据资料是位于观测线上离散点集数据，即只利用了地面沉降观测站布设的观测线上的有限观测点数据，因此所得的预测结果可靠性差，从而导致开采时地表沉陷的无法准确预知，造成环境的破坏以及经济的巨大损失。
技术实现思路
有鉴于此，本申请提供了一种基于元胞自动机的地表下沉仿真方法及系统，可以实现对复杂地质采矿条件下矿产开采时的地表沉陷更加准确地预测。为了实现上述目的，现提出的方案如下：一种基于元胞自动机的地表下沉仿真方法，包括：获取地面空间数据以及开采工作面布置情形和采矿地质数据；将所述地面空间数据进行空间插值，形成预设分辨率的格网DEM面域数据；对所述格网DEM面域数据以及开采工作面布置情形、采矿地质数据进行分析以确定开采下沉影响区域，构建元胞自动机模型；利用构建的所述元胞自动机模型对地表下沉过程进行仿真。优选的，所述地面空间数据包括：地面三维激光扫描采集的数据和/或摄影测量获取的数据以及沉降观测站采集的数据。优选的，所述构建元胞自动机模型包括：A、将所述开采下沉影响区域进行划分，以构建元胞自动机模型运行的元胞空间，并确定所述元胞空间中的每个元胞的邻居结构；B、根据经典预计方法和/或地面观测站布置情形，在观测线上确定初始失稳点，将所述初始失稳点所在的元胞加入失稳元胞集合；C、判断所述失稳元胞集合中各元胞的稳定状态，若所述失稳元胞集合中的元胞没有达到下沉停止状态，则更新该元胞的下沉状态，并将该元胞的邻居结构中没有被引起失稳的元胞加入影响元胞集合中；D、判断所述影响元胞集合中的高危失稳元胞的稳定状态，若所述高危失稳元胞的失稳概率大于设定阈值，则更新所述高危失稳元胞的下沉状态，并将该元胞加入所述失稳元胞集合；E、迭代执行上述步骤C～D，直到所述失稳元胞集合为空。优选的，所述邻居结构具体为Moore型的邻居结构。优选的，还包括：绘制所述仿真结果对应的地面三维模型图和/或下沉等值线图和/或三维剖面图。一种基于元胞自动机的地表下沉仿真系统，包括：获取单元、空间插值单元、分析单元、构建单元以及仿真单元，其中：所述获取单元，用于获取地面空间数据以及开采工作面布置情形和采矿地质数据；与所述获取单元相连的所述空间插值单元，用于将所述地面空间数据进行空间插值，形成预设分辨率的格网DEM面域数据；与所述获取单元、空间插值单元相连的所述分析单元，用于对所述格网DEM面域数据以及开采工作面布置情形、采矿地质数据进行分析以确定开采下沉影响区域；与所述分析单元相连的所述构建单元，用于根据所述分析单元确定的开采下沉影响区域构建元胞自动机模型；与所述构建单元相连的所述仿真单元，用于利用构建的所述元胞自动机模型对地表下沉过程进行仿真。优选的，所述地面空间数据包括：地面三维激光扫描采集的数据和/或摄影测量获取的数据以及沉降观测站采集的数据。优选的，所述构建单元包括：划分单元、确定单元、失稳元胞集合判断单元、影响元胞集合判断单元以及迭代单元，其中：所述划分单元，用于将所述开采下沉影响区域进行划分，以构建元胞自动机模型运行的元胞空间，并确定所述元胞空间中的每个元胞的邻居结构；与所述划分单元相连的所述确定单元，用于根据经典预计方法和/或地面观测站布置情形，在观测线上确定初始失稳点，将所述初始失稳点所在的元胞加入失稳元胞集合；与所述确定单元相连的所述失稳元胞集合判断单元，用于判断所述失稳元胞集合中各元胞的稳定状态，若所述失稳元胞集合中的元胞没有达到下沉停止状态，则更新该元胞的下沉状态，并将该元胞的邻居结构中没有被引起失稳的元胞加入影响元胞集合中；与所述失稳元胞集合判断单元相连的所述影响元胞集合判断单元，用于判断所述影响元胞集合中的高危失稳元胞的稳定状态，若所述高危失稳元胞的失稳概率大于设定阈值，则更新所述高危失稳元胞的下沉状态，并将该元胞加入所述失稳元胞集合；与所述失稳元胞集合判断单元、影响元胞集合判断单元相连的所述迭代单元，用于控制所述失稳元胞集合判断单元以及所述影响元胞集合判断单元迭代执行判断过程，直到所述失稳元胞集合为空。优选的，所述邻居结构具体为Moore型的邻居结构。优选的，还包括：绘制单元，用于绘制所述仿真单元输出的仿真结果所对应的地面三维模型图和/或下沉等值线图和/或三维剖面图。从上述的技术方案可以看出，本申请公开了一种基于元胞自动机的地表下沉仿真方法及系统，该方法通过获取地面空间数据以及开采工作面布置情形和采矿地质数据，并对地面空间数据进行空间插值，形成预设分辨率的格网DEM面域数据，然后对格网DEM面域数据以及开采工作面布置情形、采矿地质数据进行分析以确定开采下沉影响区域，从而构建元胞自动机模型，以对地表下沉过程进行仿真。该方法能够仿真较为复杂地质采矿条件下的矿产开采引起的下沉过程，动态模拟开采沉陷的开采影响范围、下沉活跃的强度和面积等，并且得到接近真实情况的地表下沉结果，从而在很大程度上提高了对复杂地质采矿条件下矿产开采时的地表沉陷预测的准确性，而且能够更加精细的刻画地表某一时刻下沉盆地的。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。图1为本申请实施例一公开的一种基于元胞自动机的地表下沉仿真方法的流程图；图2为本申请实施例一公开的一种构建元胞自动机模型的流程图；图3为本申请实施例二公开的一种基于元胞自动机的地表下沉仿真方法的流程图；图4(a)为本申请实施例二公开的一种元胞自动机模型绘制的下沉等值线图；图4(b1)为本申请实施例二公开的一种元胞自动机模型仿真模拟的沉陷区三维模型图；图4(b2)为本申请实施例二所仿真模拟绘制的三维沉陷区模型图对应的同期地面三维激光扫描获取的真实的沉陷区三维模型图；图5为本申请实施例三公开的一种基于元胞自动机的地表下沉仿真系统的结构示意图；图6为本申请实施例三公开的一种构建单元的结构示意图；图7为本申请实施例四公开的一种基于元胞自动机的地表下沉仿真系统的结构示意图。本文档来自技高网...

【技术保护点】
一种基于元胞自动机的地表下沉仿真方法，其特征在于，包括：获取地面空间数据以及开采工作面布置情形和采矿地质数据；将所述地面空间数据进行空间插值，形成预设分辨率的格网DEM面域数据；对所述格网DEM面域数据以及开采工作面布置情形、采矿地质数据进行分析以确定开采下沉影响区域，构建元胞自动机模型；利用构建的所述元胞自动机模型对地表下沉过程进行仿真。

【技术特征摘要】
1.一种基于元胞自动机的地表下沉仿真方法，其特征在于，包括：获取地面空间数据以及开采工作面布置情形和采矿地质数据；将所述地面空间数据进行空间插值，形成预设分辨率的格网DEM面域数据；对所述格网DEM面域数据以及开采工作面布置情形、采矿地质数据进行分析以确定开采下沉影响区域，构建元胞自动机模型；利用构建的所述元胞自动机模型对地表下沉过程进行仿真；其中，所述构建元胞自动机模型包括：A、将所述开采下沉影响区域进行划分，以构建元胞自动机模型运行的元胞空间，并确定所述元胞空间中的每个元胞的邻居结构；B、根据经典预计方法和/或地面观测站布置情形，在观测线上确定初始失稳点，将所述初始失稳点所在的元胞加入失稳元胞集合；C、判断所述失稳元胞集合中各元胞的稳定状态，若所述失稳元胞集合中的元胞没有达到下沉停止状态，则更新该元胞的下沉状态，并将该元胞的邻居结构中没有被引起失稳的元胞加入影响元胞集合中；D、判断所述影响元胞集合中的高危失稳元胞的稳定状态，若所述高危失稳元胞的失稳概率大于设定阈值，则更新所述高危失稳元胞的下沉状态，并将该元胞加入所述失稳元胞集合；E、迭代执行上述步骤C～D，直到所述失稳元胞集合为空。2.根据权利要求1所述的地表下沉仿真方法，其特征在于，所述地面空间数据包括：地面三维激光扫描采集的数据和/或摄影测量获取的数据以及沉降观测站采集的数据。3.根据权利要求1所述的地表下沉仿真方法，其特征在于，所述邻居结构具体为Moore型的邻居结构。4.根据权利要求1所述的地表下沉仿真方法，其特征在于，还包括：绘制所述仿真模结果对应的地面三维模型图和/或下沉等值线图和/或三维剖面图。5.一种基于元胞自动机的地表下沉仿真系统，其特征在于，包括：获取单元、空间插值单元、分析单元、构建单元以及仿真单元，其中：所述获取单元，用于获取地面空间数据以及开采工作面布置情形和采矿地质数据；与所述获取单元相连的所述空间插值单元，用于将所述地面空间数据进行空间插值，形成预设分辨率的格网DEM面域数据；与所述获取单元、空间插值...

【专利技术属性】
技术研发人员：葛小三，丁安民，张彦，付治河，王石岩，拓海燕，
申请(专利权)人：河南理工大学，
类型：发明
国别省市：河南;41