Similar simulation experiment of water loss and settlement in thick loose aquifer
-
摘要: 目前缺乏对厚松散含水层地质采矿条件下覆岩破断及变形规律的深入研究。以淮南矿区潘四东煤矿11111工作面为工程背景，构建相似材料模型，采用数字摄影测量提取位移法记录模型开挖过程中覆岩破断过程及覆岩变形情况。分析了含水层失水沉降原因：覆岩在W型剪切应力拱作用下形成2条纵向的主导水裂隙带，导水裂隙带的进一步发育引起含水层失水固结，在厚松散层重力作用下进一步压实，随着覆岩破断运动的加剧，在弯曲带和覆岩共同挤压下形成О型剪切应力拱，压缩薄层空间，导致地表下沉量增大。分析了失水状态下覆岩损伤情况：工作面开采工作完成且覆岩达到稳态后，前垮落角为57°，后垮落角为62°，导水裂隙带高度为63 m，开切眼及终采线上方覆岩在应力集中作用下断裂，产生纵向裂隙，开切眼及终采线上方垮落带区域内覆岩产生横向离层裂隙，纵向裂隙和横向离层裂隙加剧了覆岩与含水层间的水力联系。给出了失水状态下覆岩动态运动规律：随着开采工作面的推进，各观测线覆岩下沉量逐渐增大，接近开采工作面的观测线覆岩下沉量最大，工作面上方覆岩的观测线下沉量曲线走势基本类似且跳变一致，含水层上方的观测线下沉量曲线走势基本吻合且跳变同步，工作面上方与含水层上方的观测线下沉量跳变异步，表明含水层对覆岩移动变形具有重要作用。Abstract: The in-depth research on the breaking and deformation law of overburden rock under the geological and mining conditions of thick loose aquifer are lacking at present. Taking 11111 working face of Pansidong Coal Mine in Huainan mining area as the engineering background, the similar material model is constructed, and the digital photogrammetry extraction displacement method is used to record the overburden rock breaking process and overburden rock deformation during the model roadway heading. The causes of water loss and settlement of aquifer are analyzed. The overburden rocks form two main longitudinal diversion fissure zones under the action of W-type shear stress arch. The further development of the diversion fissure zone causes water loss and consolidation of the aquifer, and the aquifer is further compacted under the action of gravity of the thick loose layer. With the intensification of the overburden rock breaking movement, О type shear stress arch is formed under the joint extrusion of bending zone and overburden rock, which compresses the thin space and leads to the large amount of surface subsidence. The damage of overburden rock under water loss condition is analyzed. After the roadway heading work of the working face is completed and the overburden rock reaches a steady state, the front caving angle is 57°, the rear caving angle is 62°, and the height of the diversion fissure zone is 63 m. Under the action of stress concentration, the overburden rock above the open-cut hole and the stop-mining line is broken to produce longitudinal fissure, and the overburden rock in the area of the collapse zone above the open-cut hole and the stop-mining line produces lateral separation fissure. The longitudinal fissures and lateral separation fissures intensify the hydraulic connection between overburden rock and the aquifer. The dynamic movement law of overburden rock under water loss state is given. With the advance of mining face, the overburden settlement of each observation line increases gradually, and the overburden settlement of the observation line close to the working face is the largest. The trend of the subsidence curves of the observation lines in the overburden rock above the working face is basically similar, and the jump of the subsidence curves is consistent. The trend of the subsidence curves of the observation lines above the aquifer is basically consistent, and the jump of the subsidence curves is synchronous. The jump of the subsidence curves of observation lines in the overburden rock above the working face and the one of the observation line above the aquifer are asynchronous, indicating that the aquifer plays an important role in the movement and deformation of the overburden rock.
-
0. 引言
放顶煤开采技术[1-3]是高效开采厚煤层的重要方法，于1982年引入我国，并在我国得到快速发展。目前，我国放顶煤开采大都采用简单的单轮顺序或间隔放煤工艺，导致放煤过程中煤岩分界面下降极不均匀，矸石会提前窜入煤流，使放煤工误认为顶煤已经放净，造成大量的分界面处顶煤损失。即使某些工作面为尽可能多地回收顶煤，不采取“见矸关门”的准则，放出部分矸石以提高顶煤采出率，但会导致煤炭含矸率偏高，影响原煤质量，增加运输及洗选成本。
为提高顶煤采出率，降低含矸率，诸多新型放煤工艺被广泛提出，如多口同时放煤、分段逆序放煤、多轮顺序放煤[4]等。其中，多轮顺序放煤工艺通过增加放煤轮数，使煤岩分界面均匀下降，实现了综采工作面的顶煤采出率和含矸率的改善。但多轮顺序放煤工艺在现场应用时需要对每一轮放煤时间进行精确测定与控制，如何综合运用多方面的信息来实现智能化精准放煤对于放顶煤开采至关重要。现在大部分工作面还依赖于人工放煤，现场工人通过看和听来把握放煤过程、煤矸混放过程和放矸过程，无法精准掌控放煤口的开关时机。为精准掌控放煤口开关时机，科研人员以煤矸识别技术[5-8]为核心开展了大量研究。煤矸识别传感器大多安装在刮板输送带附近，通过传感器是否识别到刮板输送带上的矸石及矸石含量来把握放煤过程、煤矸混放过程和放矸过程，进而掌控放煤口的开关，但在现场应用时存在某些局限，例如基于图像识别的方法无法解决现场粉尘大导致图像质量低下无法识别的问题，基于煤矸撞击声音识别的方法易受到现场机械噪声的影响等问题，此类方法均易受到煤层夹矸的影响而提前终止放煤，导致大量资源浪费。
为解决现有技术难以掌控放煤口开关时机的问题，文献[9]提出了基于顶煤运移跟踪仪的自动化放煤技术，每个顶煤运移跟踪仪具有独立的ID编号及无线发送模块，并按固定频率向外发送无线信号，调度室根据接收到的无线信息来确定放煤进程。但在实际应用中，顶煤运移跟踪仪仅作为标志点安置于顶煤中，无法获得更多的顶煤运移信息。为此，本文在顶煤运移跟踪仪的基础上设计了一种基于加速度计的顶煤运移时间测量系统，通过加速度计获取标签的比力值信息，并调用时间测量算法分辨顶煤的放煤进程，得到顶煤运移时间。设计了顶煤运移时间测量标签（以下简称标签）的软硬件。由于加速度计采集比力值时误差较大，搭建了以3D转台为核心，Gauss-Newton方法[10-11]为标定算法的标定平台，完成加速度计的标定工作，使标定后的加速度计能够精确地采集比力值数据；设计完成基于阈值的时间测量算法与基于LSTM（Long Short-Term Memory，长短期记忆网络）[12]的时间测量算法，并通过时间测量实验完成2种时间测量算法的误差分析。
1. 顶煤运移时间测量系统
基于加速度计的顶煤运移时间测量系统包括标签、采集器及中心计算机3个部分，如图1所示。通过该系统能够测量得到每一轮次的顶煤运移时间，指导综采工作面现场实现多轮顺序放煤。
标签放置于顶煤内部，放煤过程中跟随顶煤运动，通过内置的加速度计实时采集比力值数据，并调用时间测量算法（基于阈值的时间测量算法与基于LSTM的时间测量算法），实现顶煤运移情况监测，进而确定不同放煤阶段，计算得到不同阶段的顶煤运移时间信息；当标签从放煤口放出后，与刮板输送带产生碰撞，通过射频信号将顶煤运移时间信息向外发送。采集器包括采集器从机和采集器主机，从机安装在刮板输送带两侧，负责接收顶煤运移时间信息的射频信号并通过现场总线发送给上一级采集器；采集器主机安装在巷道集控中心，通过现场总线将数据上传至巷道集控中心计算机；巷道集控中心计算机负责根据接收到的信息确定放煤进程，进而控制电液控制系统，调度室计算机负责同步演示放煤进程，二者通过工业环网相连接。
2. 顶煤运移时间测量标签设计
2.1 硬件设计
标签需具备实时采集比力值数据的功能，能够通过无线模块向外发送顶煤运移时间数据并被采集器接收；需防爆防潮，满足煤矿本质安全要求，能够被固定在顶煤中，并随顶煤向下运动；具备一定的数据存储功能。根据以上需求完成标签的设计，设计完成后的标签由三轴加速度计、电源模块、控制模块、通信模块、存储模块等构成，如图2所示。
控制模块选择微处理器STM32L051系列芯片，该芯片功耗低、体积小、反应速度快，具备2 kB的EEPROM。三轴加速度计选择ADXL345芯片，实时采集x，y，z三轴比力值。通信模块选择nRF24L01芯片[13]，用于实现标签与采集器间的无线通信。存储模块选择FM25W256芯片，与EEPROM共同完成参数、数据的存储。电源模块配备3.6 V锂亚电池，为标签提供安全且稳定的电压，搭载MCP1700芯片，提供稳定的3.3 V直流电压，保证电路的顺畅和其他模块的正常运作。经测试，标签在全功耗工作模式下可正常工作240 h以上。
标签外壳采用直径为4 cm的PVC圆管制成，长为16 cm，防爆防潮。标签外壳同时设有锚爪结构，用于将标签固定在煤层钻孔之中。装配好的标签如图3所示。
2.2 软件设计
标签控制模块采用STM32CubeMX（v6.0.1）和keil uVision5软件完成设计。STM32CubeMX软件以HAL库为基础，用于完成芯片选型、时钟配置、管脚配置等设计。在STM32CubeMX生成的工程文件基础上，采用keil uVision5软件对标签功能进行开发。标签能够实现实时采集三轴比力值数据、断电存储及无线通信、运动时间测量等功能，标签软件流程如图4所示。
（1） 各模块初始化。
（2） 读片上EEPROM，判断标签是否已标定，若未标定，给出标签未标定指示信号，配合标定平台进行标定。否则，读取片上EEPROM内参数并写入变量。
（3） 进入比力值数据采集状态并引入变量对原始数据进行修正。
（4） 调用时间测量算法确定顶煤运移时间。
（5） 待标签与刮板输送带产生碰撞，向外发送2.4 G无线信号。
2.3 加速度计标定
比力是载体相对惯性空间的绝对加速度和引力加速度之差。3个线性加速度计垂直安装可测量三轴比力值矢量，比力值矢量和F即为重力加速度g和运动加速度a的矢量和。
$$ {\boldsymbol{F}}={\boldsymbol{a}}+{\boldsymbol{g}} $$ (1) 未标定前，ADXL345芯片采集处于不同摆放姿势下的静态比力值矢量和F（单位为重力加速度g）曲线如图5所示。可看出加速度计在未进行标定的情况下，采集到的比力值矢量和存在失真的情况，影响运动阶段确定与顶煤运移时间测量。
![]()
图 5 未标定加速度计的比力值矢量和F测量曲线Figure 5. Specific force vector sum F measurement curve of non-calibrated accelerometer在比力值实际测量过程中存在轴偏移、零偏及刻度失准的情况。比力值矢量和F模长的计算公式为
$$ ||{\boldsymbol{F}}|| = \sqrt {f_x^{{}2} + f_y^{{}2} + f_z^{{}2}} $$ (2) 式中fx，fy，fz为修正后的三轴加速度计比力数据。
$$ \left( {\begin{array}{*{20}{c}} {{f_x}} \\ {{f_y}} \\ {{f_{\textit{z}}}} \end{array}} \right) = \left[ {\begin{array}{*{20}{c}} {{s_x}}&0&0 \\ 0&{{s_y}}&0 \\ 0&0&{{s_{\textit{z}}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} 1&{ - {\alpha _{y{\textit{z}}}}}&{{\alpha _{{\textit{z}}y}}} \\ 0&1&{ - {\alpha _{{\textit{z}}x}}} \\ 0&0&1 \end{array}} \right]\left( {\begin{array}{*{20}{c}} {{{\boldsymbol{F}}_x}} \\ {{{\boldsymbol{F}}_y}} \\ {{{\boldsymbol{F}}_{\textit{z}}}} \end{array}} \right) + \left( {\begin{array}{*{20}{c}} {{b_x}} \\ {{b_y}} \\ {{b_{\textit{z}}}} \end{array}} \right) $$ (3) 式中：sx，sy，sz为比例因子，用于调整三轴加速度计的灵敏度；
$ \alpha $ yz，$ \alpha $ zy，$ \alpha $ zx为失准角，用于调整三轴比力数据使其互相垂直；Fx，Fy，Fz为三轴加速度计原始数据；bx，by，bz为静态偏移量，用于调整静态三轴比力数据。因静止时，比力值矢量和F为重力加速度g，且ADXL345芯片采集到的比力值以g为单位，故静止时||F|| = 1g，则
$$ f\left( {\theta ,{\boldsymbol{F}}} \right) = {{\boldsymbol{F}}^2} - 1 = 0 $$ (4) 式中：
$\theta = \left( \begin{gathered} {\theta _1} \\ {\theta _2} \\ {\theta _3} \\ {\theta _4} \\ {\theta _5} \\ {\theta _6} \\ {\theta _7} \\ {\theta _8} \\ {\theta _9} \\ \end{gathered} \right) = \left( \begin{gathered} {s_x} \\ {s_y} \\ {s_{\textit{z}}} \\ {\alpha _{y{\textit{z}}}} \\ {\alpha _{{\textit{z}}y}} \\ - {\alpha _{{\textit{z}}x}} \\ {b_x} \\ {b_y} \\ {b_{\textit{z}}} \\ \end{gathered} \right)$ ；$f\left( \theta \right) = \left( \begin{gathered} f\left( {\theta ,{{\boldsymbol{F}}_1}} \right) \\ f\left( {\theta ,{{\boldsymbol{F}}_2}} \right) \\ \vdots \\ f\left( {\theta ,{{\boldsymbol{F}}_n}} \right) \\ \end{gathered} \right)$ 。为方便表示，将9个修正参数用θ进行替换，Fn代表采集到的第n个姿势下的静态比力值矢量和。
$f\left( {\theta ,{\boldsymbol{F}}} \right)$ 满足非线性最小二乘法标定方法，采用Gauss-Newton方法对加速度计进行标定。目标函数公式为$$ \hat \theta = \arg \mathop {\min }\limits_\theta \left\{ {{\boldsymbol{F}}\left( \theta \right)} \right\},{\boldsymbol{F}}\left( \theta \right) = \frac{1}{2}f{\left( \theta \right)^{\rm{T}}}f\left( \theta \right) $$ (5) 标定平台由3D转台、信号接收器和计算机3个部分构成，如图6所示。3D转台有3台电动机，可独立实现3个平面的自由旋转。信号接收器负责接收待标定标签发送的无线信号。计算机负责接收并记录比力值数据，调用Gauss-Newton算法计算得到标定后的参数值。
标定过程：① 控制3D转台使标签在多个位置下保持静止，在每一位置采集100组比力值矢量和F数据，求算数平均值作为该位置的比力值矢量和F。② 将多个位置的比力值矢量和F代入Gauss-Newton算法中计算得到标定后的参数值。③ 通过无线模块将标定后参数值返回标签中并进行记录。
标定参数设置见表1。未标定设备虽在数据采集过程中具备一定的采集误差，但数据仍在合理范围内，故比例因子sx，sy，sz初始值设为[1 1 1]，失准角
$ \mathrm{\alpha } $ yz，$ \mathrm{\alpha } $ zy，$ \mathrm{\alpha } $ zx用于弥补安装误差，初始值设为[0 0 0]，静态偏移量bx，by，bz用于静态校准，初始值设为[0 0 0]。标定结果见表2。结果表明，加速度计标定后平均值接近于1，方差接近于0，具有很好的标定效果。表 1 标定参数设置Table 1. Calibration parameter setting参数 数值 采样次数 100 迭代步数itmax 10 000 收敛阈值ε 1×10−10 初始值 [1,1,1,0,0,0,0,0,0] 表 2 标定前后比力值矢量和F对比Table 2. Comparison of specific force vector sum F before and after calibration比力值矢量和 标定前 标定后 平均值/g 1.016 739 0.999 988 方差 0.013 154 0.000 004 7 采用标定后的ADXL345芯片重新采集不同摆放姿势下的标签的静态比力值矢量和F，结果如图7所示。可看出标定后的加速计能够精确地采集比力值矢量和数据。
![]()
图 7 标定前后比力值矢量和F的测量曲线Figure 7. Specific force vector sum F measurement curve before and after calibration3. 顶煤运移时间测量算法
放顶煤作业开始后，标签主要受重力影响，随着顶煤向下运动，标签在每一轮放煤过程中运动分为4个阶段，模拟变化曲线如图8所示。
![]()
图 8 放煤过程中比力值矢量和F模拟曲线Figure 8. Specific force vector sum F simulating curve in the coal drawing process阶段①：放煤开始前，标签在顶煤中处于静止状态，运动加速度a为0，比力值矢量和F为g。 阶段②：放煤开始，标签随顶煤主要受重力影响向下运动，且比力值矢量和F处于持续变化的状态。阶段③：当前轮次放煤结束，标签再次固定在顶煤中，处于静止状态，运动加速度a为0，比力值矢量和F为g。阶段④：标签从放煤口放出，与刮板输送带产生反复碰撞，比力值矢量和F处于持续变化的状态。
将标签由阶段①变化为阶段②的突变时刻记为t0 1，将标签由阶段②变化为阶段③的突变时刻记为t1 r（r=1，2，…，R，R为放煤轮数），将标签由阶段③变化为阶段②的突变时刻记为t0 r+1，将标签由阶段②变化为阶段④的突变时刻记为t1 R，标签第r轮运动时间tr为t1 r−t0 r，标签运动时间也为顶煤运移时间。
根据顶煤在放煤过程中的运移特点，本文提出了基于阈值的时间测量算法与基于LSTM的时间测量算法。基于阈值的时间测量算法通过引入静态阈值识别每一轮次运动起始点时刻t0 r，引入最大阈值识别每一轮次运动终止点时刻t1 r；基于LSTM的时间测量算法通过识别时域下比力值矢量和的动态变化，捕捉突变点并记录突变时刻，进而得到每一轮次的运动时间。
3.1 基于阈值的时间测量算法
基于阈值的时间测量算法引入静态阈值G与最大阈值M识别运动起始点时刻t0 r与运动终止点时刻t1 r。时间测量流程如下：
（1） 标签以固定频率持续读取比力值矢量和F。
（2） 当
${\boldsymbol{F}} \notin G$ 时，记录当前时刻为t0 r。（3） 继续读取F，当F＞M时，记录当前时刻为t1 r。
（4） 若F≤M，且满足
${\boldsymbol{F}} \in G$ 时，认为标签没有发生运动，回到步骤（1），清除当前轮次的t0 r。（5） 记录到t1 r后，若连续采集到
${\boldsymbol{F}} \in G$ ，回到步骤（1），令r = r + 1，准备采集下一轮次运动时间，否则认为标签已从放煤口放出，启动通信模块，向外发送各轮次放煤时间。静态阈值的确定通过3个独立的标签采集4个不同摆放姿势下的多组比力值矢量和F，计算平均值得出不同摆放姿势下静止时标签的比力值矢量和F，实验结果见表3。
表 3 时间测量标签静止实验Table 3. Stationary experiment of time measurement label标签 标签
姿势比力值矢量和 平均值/g 方差 标签1 1 1.013228 4.05×10−5 2 1.019693 1.59×10−5 3 1.000189 3.46×10−5 4 1.000916 2.61×10−5 标签2 1 0.996145 2.82×10−5 2 1.008506 2.76×10−5 3 0.99687 1.52×10−5 4 1.006726 1.76×10−5 标签3 1 0.986678 2.81×10−5 2 0.966817 1.87×10−5 3 1.031543 2.57×10−5 4 0.996493 4.03×10−5 从表3可看出，标签处于静止状态时采集到的比力值矢量和的变化范围为0.966817~1.031543g。因此，静态阈值选择0.9~1.1g即可识别t0 r。
最大阈值的确定采用使标签掉落与地面发生碰撞的方式进行测量，记录标签与地面碰撞瞬间采集到的比力值矢量和F，得到多次碰撞下的比力值矢量和F的最小值，实验结果见表4。可看出标签与地面产生碰撞时，采集到的最小比力值矢量和F为2.06g。最大阈值设置为2g可识别t1 r。
表 4 时间测量标签自由落体实验Table 4. Free-fall experiment of time measurement labelg 标签 落地时比力值矢量和 标签1 12.87, 10.10, 11.52, 15.13, 2.26, 8.47, 16.74, 4.36, 8.80, 14.14, 5.14, 15.29, 13.13, 3.46, 6.36, 14.34, 3.25, 17.38, 5.37, 13.36 标签2 13.91, 8.34, 16.37, 2.06, 2.83, 19.30, 13.45, 13.22, 15.93, 5.46, 6.23, 11.71, 2.29, 11.97, 17.03, 16.92, 15.00, 9.41, 8.18, 7.30 标签3 14.52, 2.46, 21.63, 2.28, 21.25, 15.88, 12.28, 22.32, 9.59, 16.71, 13.22, 15.65, 5.02, 5.27, 6.05, 14.62, 20.55, 12.19, 18.69, 11.90 时间测量实验采用3个独立的标签，分别进行多组高度为1 m的自由落体实验，用于检测时间测量效果，理论下落时间为451.72 ms。实验结果见表5，通过表5数据计算可得下落时间平均值为510.8 ms，方差为0.000 6，误差为13.07%，满足现场顶煤运移时间测量需求。
表 5 基于阈值的时间测量算法的测量结果Table 5. Time measurement results of time measurement algorithm based on state thresholdms 标签 下落时间 标签1 510.9, 529.9, 509.0, 513.0, 487.9, 487.9, 529.9, 509.0, 510.9, 490.0, 487.9, 552.9, 532.0, 553.0, 510.9, 513.0, 490.0, 490.0, 487.9, 472.9 标签2 515.0, 515.0, 463.9, 509.0, 467.0, 529.9, 509.0, 509.0, 551.0, 490.0, 510.9, 469.0, 490.0, 487.9, 509.0, 552.0, 470.9, 490.0, 510.9, 490.0 标签3 529.9, 509.0, 536.0, 469.0, 534.0, 552.0, 529.9, 532.0, 555.0, 555.0, 509.0, 532.0, 487.9, 467.0, 509.0, 551.0, 529.9, 509.0, 533.9, 509.0 3.2 基于LSTM的时间测量算法
LSTM是一种时间循环的神经网络，解决了一般的RNN（Recurrent Neural Network，循环神经网络）中存在的长期依赖问题。所有的RNN都具有一种重复神经网络模块的链式形式，在标准RNN中，这个重复的结构模块仅有一个非常简单的结构，例如一个tanh层。而LSTM由于独特的设计结构能够更适合于预测和处理时间序列中间隔和延迟非常长的重要事件。
基于LSTM的时间测量原理通过调用训练好的模型寻找比力值矢量和曲线突变时的采样点，进而识别并记录运动起始点时刻t0 r与运动终止点时刻t1 r。
训练数据集为标签在不同姿势下采集到的共1 000组静态比力值矢量和，学习率为0.001，批大小为32，最大迭代次数设置为100，LSTM模型训练过程中损失值变化情况如图9所示。经过100次迭代，模型收敛，损失值收敛到1.87×10−4。
LSTM模型的测试结果如图10所示，蓝色曲线为原始采集比力值矢量和曲线，橙色曲线为调用LSTM模型后，根据比力值矢量和曲线得出的异常识别曲线。可看出，训练后的模型能够较好地在曲线突变处做出反应，能够及时检测出曲线的变化并用红色标志点标出。与此同时，程序输出状态切换位置：35，37，38，39，44，48，49，50，53，55，56，65，66，67，70，72，73，74，81，82，83，85，86，87，88，89，90，91，100，101，106，110，111，112，115，117，118，127，128，129，132，134，135，136，143，144，145，147，148，149，150，151，152，153，154，155，157，158，159，160，161，162，163，172，173，178。训练与测试结果表明，LSTM算法能够较好地判断标签是否发生运动及运动状态是否发生改变，进而识别标签在放煤过程中的运动阶段。
![]()
图 10 LSTM算法识别曲线突变点Figure 10. Mutation points of identification curve based on LSTM algorithm时间测量实验采用3个独立的标签，分别进行多组高度为1 m的自由落体实验，用于检测时间测量效果，实验结果见表6。理论下落时间为451.72 ms。经计算得到下落时间平均值为475.3 ms，方差为0.0002，误差为5.22%，满足现场顶煤运移时间测量需求。
表 6 基于LSTM的时间测量算法的时间测量结果Table 6. Time measurement results of time measurement algorithm based on LSTM algorithmms 标签 下落时间 标签1 466.4, 484.7, 493.6, 484.7, 492.1, 466.7, 484.7, 475.7, 493.7, 484.7, 475.8, 466.8, 475.8, 475.7, 493.7, 484.5, 493.9, 484.8, 484.9, 484.8 标签2 475.6, 493.4, 494.0, 457.5, 448.8, 448.8, 440.1, 457.8, 457.8, 457.9, 475.5, 458.0, 475.7, 475.7, 479.1, 457.7, 484.6, 475.8, 466.5, 484.9 标签3 466.4, 484.7, 493.6, 484.7, 492.1, 466.7, 484.7, 475.7, 493.7, 484.7, 475.8, 466.8, 475.8, 475.7, 493.7, 484.5, 493.9, 484.8, 484.9, 484.8 相较于基于阈值的时间测量算法，基于LSTM的时间测量算法具有以下优势：
（1） 基于阈值的时间测量算法标签与物体碰撞时的比力值矢量和波动较大，难以确定合适的最大阈值；基于LSTM的时间测量算法采用异常检测方式检测比力值矢量和曲线是否出现与静态时曲线产生差异，能够及时准确地定位曲线突变位置。
（2） 基于LSTM的时间测量算法在实际应用中，时间数据波动范围较小，误差较小。
（3） 基于阈值的时间测量算法在实际应用中，标签需按照预定好的运动方式运动才能够得到准确的时间数据；基于LSTM的时间测量算法能够从异常采样点观测标签的比力值矢量和的变化情况，进而得出标签的运动情况，根据异常采样点与采样时间能够确定多种运动时间。
4. 结论
（1） 设计了顶煤运移时间测量标签。标签采用微处理器STM32L051系列芯片，并搭载加速度计等多种模块，完成了各模块与微处理器的通信，实现了比力值实时采集、断电存储及无线通信等功能。
（2） 采用Gauss-Newton方法完成加速度计的标定工作，标定后的加速度计能够准确采集标签的比力值矢量和。
（3） 根据顶煤运移时间测量原理，提出了基于阈值的时间测量算法与基于LSTM的时间测量算法。通过标签的自由落体实验完成了2种时间测量算法的性能测试，其中时间测量方差分别为0.0006、0.0002，时间测量误差分别为13.07%、5.22%，满足现场顶煤运移时间测量需求。通过实验对比可知，基于LSTM的时间测量算法在顶煤运移时间测量方面时间数据波动范围较小，误差较小，可获取多种运动时间，具有明显的应用优势。
-
![]()
图 1 相似材料模拟实验方案及模型
Figure 1. Scheme and model of Similar material simulation experiment
![]()
图 7 含水层失水状态下各观测线的下沉量变化
Figure 7. The variation of subsidence of each observation line under water loss state of aquifer
表 1 研究区地层结构
Table 1 A histogram of stratigraphic structure in the study area
界 系 统 组 厚度(m) 主要岩性 新生界 第四系 全新统 — 40～130 浅黄、灰黄色粘土夹砂层 更新统 第三系上 上新统 — 0～152 灰绿、浅黄，多为粘土夹杂砂层 中新统 第三系下 渐新统 — >205 浅灰、棕色砂泥岩互层，夹杂砂砾岩 始新统 中生界 白垩系 上统 — >647 紫红色粉、细砂岩，砂砾岩 下统 — 844 棕红粉砂岩、泥岩及细中粒砂岩 侏罗系 上统 — >637 凝灰质砂砾岩，凝灰岩和安山岩 三叠系 下统 — 316～446 紫红色砂泥岩 古生界 二叠系 上统 石千峰组 114～260 紫红杂色砂泥岩，夹石英砂岩及砂砾岩 上石盒子组 316～566 灰绿色和浅灰色砂岩，底为石英砂岩且为含煤层 下统 下石盒子组 106～265 灰色砂泥岩及其互层，底含粗砂岩，含煤层 山西组 52～88 上部细砂岩、粗砂岩，下部深灰色泥岩，含煤层 石炭系 上统 太原组 102～148 灰岩为主，夹泥岩及砂岩，含薄煤层 奥陶系 中下统 — 400 中厚层为白云岩及白云质灰岩，夹灰岩 寒武系 上统 土坝组 170～220 硅质结核白云岩，产Heleionellasp.化石 固山组 9～78 白云岩，竹叶状灰岩，鲕状灰岩。 中统 张夏组 146 鲕状灰岩，白云岩产Dameselluasp.化石 徐庄组 190 棕黄砂岩，夹页岩及石灰岩 毛庄组 152 多为砾状灰岩，鲕状灰岩和页岩 下统 馒头组 215 紫色页岩夹灰岩，产Redlichasp.化石 猴家山组 100～150 鲕状灰岩，孔洞灰岩和砂灰岩 凤台组 10～100 页岩，砾岩 上元古界 震旦系 徐淮群 九顶山组 117 白云岩，底部夹竹叶状灰岩 倪园组 92 上部含泥白云岩，夹黄绿色钙质页岩，下部硅质条带白云岩 四顶山组 137 厚层白云岩为主，产蠕形动物化石 九里桥组 119 泥灰岩，砂灰岩 四十里长山组 93 石英岩及钙质砂岩 青白口系 八公
山群刘老碑组 1050 页岩，泥灰岩，石英砂岩，底部铁质砂砾岩，含藻及疑源类化石 伍山组 张店组 下元古界 — — 凤阳群 1 171 千枚岩，白云岩，大理岩，白云质石英片岩，石英岩，含藻化石 上太古界 — — 五河群 >6 422 片麻岩，浅粒岩，变粒岩，斜长角闪岩互层，夹少量大理岩及磁 铁矿层，岩石混合岩化 
下载: 导出CSV
-
[1] 谢和平, 吴立新, 郑德志. 2025 年中国能源消费及煤炭需求预测[J]. 煤炭学报,2019,44(7):1949-1960. XIE Heping, WU Lixin, ZHENG Dezhi. Prediction on the energy consumption and coal demand of China in 2025[J]. Journal of China Coal Society,2019,44(7):1949-1960.
[2] 高永格, 牛矗, 张强, 等. 厚松散层下采煤地表沉陷特征研究[J]. 煤炭科学技术,2019,47(6):192-198. GAO Yongge, NIU Chu, ZHANG Qiang, et al. Study on surface subsidence characteristics of coal mining under thick loose layer[J]. Coal Science and Technology,2019,47(6):192-198.
[3] 张海峰, 李文, 李少刚, 等. 浅埋深厚松散层综放工作面覆岩破坏监测技术[J]. 煤炭科学技术,2014,42(10):24-27. ZHANG Haifeng, LI Wen, LI Shaogang, et al. Research on technology of overlying strata failure monitoring in fully mechanized caving coal face with shallow and thick loose bed[J]. Coal Science and Technology,2014,42(10):24-27.
[4] 陈俊杰, 邹友峰, 郭文兵. 厚松散层下下沉系数与采动程度关系研究[J]. 采矿与安全工程学报,2012,29(2):250-254. DOI: 10.3969/j.issn.1673-3363.2012.02.018 CHEN Junjie, ZOU Youfeng, GUO Wenbing. Study on the relationship between subsidence coefficient and mining degree under a thick alluvium stratum[J]. Journal of Mining & Safety Engineering,2012,29(2):250-254. DOI: 10.3969/j.issn.1673-3363.2012.02.018
[5] 李新岭, 郭文兵, 赵高博. 巨厚松散层土体压缩特性对开采沉陷影响研究[J]. 中国安全科学学报,2018,28(7):135-141. LI Xinling, GUO Wenbing, ZHAO Gaobo. Study on influence of compression characteristics of super-thick alluvium on mining subsidence[J]. China Safety Science Journal,2018,28(7):135-141.
[6] 许国胜, 许胜军, 李德海, 等. 赵固矿区厚松散层下开采地表移动特征实测研究[J]. 中国矿业, 2021, 30(5): 168-172. XU Guosheng, XU Shengjun, LI Dehai, et al. Field research on surface movement characteristic by coal mining under thick alluvium in Zhaogu mine area[J]China Mining Magazine, 2021, 30(5): 168-172.
[7] 许国胜, 李德海, 侯得峰, 等. 厚松散层下开采地表动态移动变形规律实测及预测研究[J]. 岩土力学,2016,37(7):2056-2062. XU Guosheng, LI Dehai, HOU Defeng, et al. Measurement and prediction of the transient surface movement and deformation induced by mining under thick alluvium[J]. Rock and Soil Mechanics,2016,37(7):2056-2062.
[8] 崔希民, 方志海, 左红飞, 等. 开采引起的含水层失水对地表下沉的影响[J]. 煤田地质与勘探,2000,28(5):47-48. DOI: 10.3969/j.issn.1001-1986.2000.05.015 CUI Ximin, FANG Zhihai, ZUO Hongfei, et al. De-watering effect of aquifer to surface subsidence[J]. Coal Geology & Exploration,2000,28(5):47-48. DOI: 10.3969/j.issn.1001-1986.2000.05.015
[9] 崔希民, 邓喀中. 煤矿开采沉陷预计理论与方法研究评述[J]. 煤炭科学技术,2017,45(1):160-169. CUI Ximin, DENG Kazhong. Research review of predicting theory and method for coal mining subsidence[J]. Coal Science and Technology,2017,45(1):160-169.
[10] 李春意, 郭增长. 含水松散层固结对地表沉陷影响的数值分析[J]. 河南理工大学学报(自然科学版),2012,31(1):48-53. LI Chunyi, GUO Zengzhang. Numerical analysis of water-bearing alluvium consolidation influence on surface subsidence[J]. Journal of Henan Polytechnic University(Natural Science),2012,31(1):48-53.
[11] 张丁丁. 兖州矿区第四系厚松散层沉降特性研究[D]. 西安: 西安科技大学, 2015. ZHANG Dingding. Study on settlement characteristics of thick quaternary unconsolidated layer in Yanzhou mining area. [D]. Xi′an: Xi′an University of Science and Technology, 2015.
[12] 张劲满. 厚松散层含水层失水开采条件下覆岩运动及地表移动规律研究[D]. 淮南: 安徽理工大学, 2019. ZHANG Jinman. Study on overburden movement and surface movement under water loss mining in thick loose aquifer[D]. Huainan: Anhui University of Science and Technology, 2019.
[13] 杜锋, 白海波. 厚松散层薄基岩综放开采覆岩破断机理研究[J]. 煤炭学报,2012,37(7):1105-1110. DU Feng, BAI Haibo. Mechanism research of overlying strata activity with fully mechanized caving in thin bedrock with thick alluvium[J]. Journal of China Coal Society,2012,37(7):1105-1110.
[14] 梁庆华, 温兴林, 何刚, 等. 黏土体失水引起的地表下沉计算方法研究[J]. 采矿与安全工程学报,2007,24(1):105-108. DOI: 10.3969/j.issn.1673-3363.2007.01.023 LIANG Qinghua, WEN Xinglin, HE Gang, et al. Study on calculation methods for surface subsidence caused by water loss of clay[J]. Journal of Mining & Safety Engineering,2007,24(1):105-108. DOI: 10.3969/j.issn.1673-3363.2007.01.023
[15] 赵春虎, 靳德武, 王皓, 等. 榆神矿区中深煤层开采覆岩损伤变形与含水层失水模型构建[J]. 煤炭学报,2019,44(7):2227-2235. ZHAO Chunhu, JIN Dewu, WANG Hao, et al. Construction and application of overburden damage and aquifer water loss model in medium-deep buried coal seam mining in Yushen mining area[J]. Journal of China Coal Society,2019,44(7):2227-2235.
[16] 徐良骥, 朱楠, 马荣振, 等. 厚松散承压含水层失水沉降模拟实验研究[J]. 采矿与安全工程学报,2015,32(5):821-826. XU Liangji, ZHU Nan, MA Rongzhen, et al. Water loss settlement simulation of thick unconsolidated confined aquifer layer[J]. Journal of Mining & Safety Engineering,2015,32(5):821-826.
[17] 周大伟. 煤矿开采沉陷中岩土体的协同机理及预测[D]. 徐州: 中国矿业大学, 2014. ZHOU Dawei. The synergy mechanism between rock mass and soil in mining subsidence and its prediction[D]. Xuzhou: China University of Mining and Technology, 2014.
-
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 109
- HTML全文浏览量: 107
- PDF下载量: 18
- 被引次数: 1
