回采工作面煤层动态标定预测技术应用研究

2022-02-22刘文明高耀全蒋必辞薛悟强

煤田地质与勘探 2022年1期

刘文明，高耀全，蒋必辞 ,2，李鹏 ,3，薛悟强

（1.中煤科工集团西安研究院有限公司，陕西西安 710077；2.煤炭科学研究总院，北京 100013；3.西安科技大学地质与环境学院，陕西西安 710054；4.鄂尔多斯市华兴能源有限责任公司，内蒙古鄂尔多斯 017000）

我国煤矿智能化开采基本上实现了从单机自动化向“三机”协同自动化以至于局部系统智能化的技术跨越[1-2]，但智能开采对于复杂地质条件的适应性远远不足[3]，究其原因，一是煤矿开采重型装备的自动化与智能化水平、远程控制系统的自主决策能力有待提高[4]；二是工作面基础地质条件查明精度不足，导致智能开采遭遇地质条件较大变化时无法迅速调整适应等[5]。

智能化开采对煤层厚度、煤层起伏形态提出了更高的精度要求[6]。构建高精度三维地质模型成为实现智能化的关键技术问题之一[7]。由于工作面内部已知数据较少，采用传统的插值算法无法有效预测工作面煤层起伏。三维地震勘探能够对煤层起伏、埋深、厚度和构造等进行高精度探测[8-10]，但由于其标定数据的稀缺性，导致反演得到的煤层底板精度受限[11]。回采工作面两侧巷道具有大量的探煤厚数据，同时随着切眼不断推进，会不断产生探煤厚数据，利用这些数据对三维地震进行动态解释和标定，可不断提高煤层底板反演的精度，从而从根本上提高地质建模精度。同时煤层厚度较为稳定的工作面，可通过插值算法迭代更新煤厚不断降低预测误差[12]。

前人对三维地震勘探技术开展了大量的研究工作，但是对于面向智能化开采的地震动态解释技术还没有进行系统的研究。本文以三维地震数据为基础，利用回采工作面巷道和切眼得到的探煤厚数据，重新标定三维地震数据，更新平均速度场，更新回采工作面煤层底板模型，同时采用离散光滑插值不断迭代工作面煤层厚度，利用假设的未回采处探煤厚数据进行精度对比，研究更新过程中精度变化及存在的问题。

1 三维地震动态解释原理

三维地震动态解释是在常规的地质分析、反射波标定、波组对比、层位构造解释的基础上进行的，由于矿井在采掘过程中可不断得到煤层顶底的空间坐标，从而能够对该煤层的底板反射波和煤层厚度进行重新标定校正。

工作面煤层底板反射波为密度较高的网格化数据，若其采集处理方法得当，其成果精度能够控制煤层起伏。由于反射波记录为时间域数据，需要利用已知煤层底板高程值构建速度场，进行时深转换。由于煤田声波测井资料稀缺，常采用平均速度场方法，其精度主要受标定点的密度和均匀度影响。

计算标定点处的速度值公式为：

式中：v(i)为第i个标定点处平均速度；Wb(i)为第i个标定点处煤层底板高程；t(i)为第i个标定点处煤层底板反射波双程旅行时；H为地震数据处理基准面。

对所有标定点计算得到的速度值进行内插，得到平均速度网格数据v(i,j)，按照时深转换公式可计算出目标层位底板高程：

式中：H为三维地震处理基准面；v(i,j)为第(i,j)个内插点的速度；t(i,j)为第(i,j)个解释点处目标层反射波地震双程旅行时；Wb(i,j)为第(i,j)个点目标层煤层底板高程。

采用离散光滑插值算法，不断融合采掘揭露的数据，更新工作面煤厚值。离散光滑插值是一种无维数的内插方法，它不以空间坐标为参数，主要依赖于网格结点的拓扑关系，通过使网格点满足特定约束条件，求解一个线性方程得到未知结点上的值[13]。

首先建立网格结点φ最优解目标函数：

式中：R(φ)为全局粗糙度函数；ρ(φ)为线性约束违反度函数，其约束点为插值点。

使式（3）无限逼近最小值，达到φ在任意结点的函数值逼近该点领域结点φ值的均值，使插值点处的结点的预测值尽可能逼近插值数据。

通过求解全局粗糙度函数使网格结点尽可能平滑：

式中：μ(k)为在结点k∈Ω上的权系数，可调整局部光滑度；a为Λ(k)的结点；Λ(k)为结点k邻域内不包含k的结点的集合；va(k)为在结点k的邻域内结点a的权系数；φ(a)为结点a处的真实值。

通过违反度函数进行结点约束：

式中：a∈Ω ；为给定的正数，Ai(a)φ(a)-bi(a)为结点a关于φ的第i个结点的线性约束；M为结点数量；Ai(a)和bi(a)为给定常数。

具体插值步骤如下。

①准备散点数据的三维空间坐标（x,y,z）。

② 以散点数据的平面坐标为边界范围，构建平面网格，网格间距分别为dx和dy，网格点初始纵向坐标值取所有散点数据z值的平均值。

③设置插值约束条件为：已知点之间曲率最小、已知点处的解算值小于设置的阈值；通过插值算法迭代求解各网格节点处的z值。

④ 将网格之间进行填充可视化。

2 三维地震动态解释实例

本次试验工作面为TJH304 工作面，主采6 号煤层，煤层埋深为469.7～533.1 m，煤层厚度为11.6～17.5 m，煤层倾角为0～6°，煤层整体为一单斜构造，呈东南高西北低，局部有褶曲。工作面倾向长度为240 m，走向长度为2 141 m，工作面辅运巷、运输巷和切眼共有72 个探煤厚点，工作面内部有1 个钻孔，采煤面共有44 个探煤厚点（图1）。

图1 TJH304 工作面分布Fig.1 Sketch map of TJH304 longwall panel

三维地震数据覆盖工作面全区，CDP 网格5 m×5 m，覆盖次数80 次。图2 为工作面推进方向地震剖面，图中T1 层位为6 号煤层底板反射波，反射波能量强、信噪比高，煤层局部有起伏。图3 为运输巷、辅运巷和切眼探煤后得到的6 号煤层底板高程和其对应的双程旅行时交会图，两者相关性可达62.26%。相关性较好时，说明平均速度变化不大，采用平均速度法进行时深转换比较稳定，可作为检验数据质量的一个判断准则。当煤层上覆地层存在速度异常体或较大的横向速度变化时，该处平均速度出现突变，会对时深转换精度影响较大。

图2 TJH304 工作面推进方向地震剖面（A-A’）Fig.2 Seismic section map of TJH304 working face along mining direction（A-A’）

图3 巷道及切眼煤层底板测量点高程与双程旅行时相关关系Fig.3 The correlation of coal bottom elevation value and two-way time value of measure points in tunnel and initial mining surface

在工作面形成后，还未开采前，利用辅运巷、运输巷和切眼得到的煤层底板高程作为新的时深转换标定点，对T1 反射层位平均速度场进行更新（图4a），时深转换得到新的6 号煤层底板高程数据，利用工作面开采后回采面获得的煤层底板高程，对验证点处更新前后高程误差进行统计分析；在工作面开始回采后，分别选取回采进尺300 m 和329 m 以内的探煤厚数据，继续更新T1 反射层位平均速度场（图4b），更新6 号煤层底板高程数据，分别利用回采进尺300 m 和329 m 以后的探煤厚点验证点，对更新前后高程误差进行统计分析。

同理，对于煤层厚度，采用以上过程，利用离散光滑插值对6 号煤层厚度数据迭代更新（图5），并进行误差统计分析。

由图4 可知，平均速度场由于采煤面标定点的增多，在0～300 m 范围速度变化更加精细。由图5 可知，煤层厚度由于采煤面厚度点的增多，在0～300 m 范围煤层厚度变化更加明显。

图4 不同开采阶段煤层底板高程数据更新前后的平均速度场平面图Fig.4 The average velocity comparison map updated by coal bottom elevation value of mining working face in different mining stages

图5 不同开采阶段煤层厚度数据更新前后的煤层厚度平面图Fig.5 The coal seam thickness comparison map updated by coal seam thickness value of mining working face in different mining stages

3 面向智能开采的解释精度分析

基于三维地震解释成果和探煤厚数据开展回采工作面煤层底板和厚度动态预测，其精度主要受以下因素影响：地表障碍物、复杂地表、表浅层地层非均质性影响，地层速度异常、速度横向突变，煤层局部起伏较大，煤层厚度局部变化较大，探煤厚点数量等因素。

表1 为不同开采阶段煤层底板与煤层厚度预测绝对误差统计。可以看出在回采前，利用速度场刷新和煤厚迭代插值，煤层底板绝对误差可控制在0.42～3.63 m，与常规三维地震解释精度相比具有较大提高，煤层厚度误差可控制在0.27～1.89 m 内。在进入回采后300 m 处，利用速度场刷新和煤厚迭代插值，煤层底板绝对误差可控制在0.34～3.26 m，煤层厚度误差可控制在0.11～1.76 m 内；在回采标定点最前方30 m 和29 m处的4 个验证点，煤层底板高程误差范围为0.37～0.58 m，煤层厚度误差范围为0.32～0.44 m。