一种面向智能开采的地层界面插值建模方法技术

本发明专利技术公开了一种面向智能开采的地层界面插值建模方法：步骤1：将工作面测量的数据整理为集合α；步骤2：将集合α进行网格化处理形成网格点的集合β1；步骤3：得到插入点集φ1中插入点坐标；步骤4：设置点集φ1中每个点的初始z坐标；步骤5：设置点集φ1和集合β1中每个点的权重；步骤6：定义点集φ1和集合β1中各点的势η

【技术实现步骤摘要】
一种面向智能开采的地层界面插值建模方法


[0001]本专利技术属于煤田地质
，涉及一种地层界面插值建模的方法，用于实现地层界面网格剖分、地层界面动态更新的功能。

技术介绍

[0002]煤矿智能化、少人化是实现煤矿安全的重要途经之一。煤矿智能化经历了可视远程干预、工作面自动找直、基于透明工作面智能割煤、全智能化自适应开采四个阶段，目前正处于透明工作面智能割煤技术关键阶段。透明工作面智能割煤技术通过三维地震、井巷测量、钻孔测量等探测手段获取工作面的实际展布情况，在分析上述数据后，通过不同的插值方法建立工作面模型，然后使用工作面模型指导采煤机开采工作。其中，高精度采煤工作面煤层模型是其中的关键所在，而插值方法是工作面模型实现的必经途径。
[0003]目前常用的空间插值方法有函数插值、克里金插值和光滑离散插值。由于光滑离散插值算法的原理，并不能在X、Y方向上形成大小统一的网格，但是对于工作面来说工作面每一刀大小都是基本固定。所以在模型的使用过程中，需要网格的大小与每一刀的大小相匹配。克里金插值算法由于其算法原理，模型的更新必须全部更新，难以实现局部更新的目的，而对于工作面模型来说，模型通常较大，模型完全更新对系统提出了更高的要求，还浪费资源，如果能够实现局部动态更新便可解决上述问题。函数插值包括距离加权平均插值法(IDW)、趋势面法、样条函数法等，其插值方法通过拟合已知点的数据，通过函数预测未知点的数据，插值得到模型的精度较低。

技术实现思路

[0004]本专利技术的目的是提供一种地层界面插值建模的方法，以克服现有技术存在的不能形成大小统一的网格以及难以局部更新的问题。
[0005]为实现上述目的，本专利技术公开了一种地层界面插值建模的方法，至少包括以下步骤：
[0006]一种面向智能开采的地层界面插值建模方法，具体包括以下步骤：
[0007]步骤1：将工作面测量获取的数据进行整理，形成集合α，集合中包含三列数据(x,y,z)，x表示横坐标、y表示纵坐标，z表示高程；
[0008]步骤2：将集合α进行网格化处理形成网格点的集合β1；
[0009]步骤3：根据集合β1的x、y坐标，得到插入点集φ1中插入点的坐标；
[0010]步骤4：设置点集φ1中每个点的初始z坐标；
[0011]具体的，将集合α的平均z坐标赋值给点集φ1中各点的初始z坐标；
[0012]步骤5：设置点集φ1和集合β1中每个点的权重大小；
[0013]步骤6：定义点集φ1和集合β1中各点的势η
a
；
[0014]步骤7：令η
a
＝0，根据步骤6的定义，计算得到当前的点集φ1和集合β1中各网格点的z
a
值，得到点集φ1和点集β1的并集
[0015]步骤8，循环执行步骤7，循环次数不小于10次，得到当前的点集φ1和点集β1的并集
[0016]步骤9：采用步骤2得到的集合β1替换步骤8中得到的当前的点集φ1和点集β1的并集中的集合β1，得到新的并集
[0017]步骤10：计算步骤9得到的并集的网格间距，判断网格间距是否达到目标网格间距，若已达到，则输出并集作为最终的插值结果；否则将并集作为集合β1，返回步骤3；
[0018]进一步的，所述步骤2中的网格化处理采用函数拟合计算或者克里金插值。
[0019]进一步的，所述步骤3中，插入点的坐标的计算公式为：
[0020]x
n
＝(x
n+1
+x
n
‑1)/2
[0021]y
n
＝(y
n+1
+y
n
‑1)/2
[0022]其中，x
n
、y
n
为插入点集φ1中一个插入点的坐标值；x
n
‑1、y
n
‑1、x
n+1
、y
n+1
分别为集合β1中两个相邻的点的坐标值。
[0023]进一步的，所述步骤5中，将每个点的权重设置为1。
[0024]进一步的，所述步骤6中，点集φ1和集合β1中各点的势η
a
定义为：
[0025][0026]其中，Ο(a)表示与点a有共同边界的上、下、左、右四个网格点的集合，i为上、下、左、右四个点的序号；λ
i
为点i的权重；
[0027]若a点为中间网格，则与其有共同边的a1、a2、a3、a4四个网格点的z值参与计算；
[0028]若a点为边界网格，则与其有共同边的a1、a2、a3三个网格点的z值参与计算；
[0029]若a点为角网格，则与其有共同边的a1、a2两个网格点的z值参与计算。
[0030]进一步的，所述步骤8中，循环次数为10次。
[0031]相较于现有技术，本专利技术的方法具有如下技术效果：
[0032]1、能够根据工作面的具体情况，设置网格的大小，能够实现固定网格大小。在智能开采过程中，能根据采煤机的截深，无需经过数据处理，直接提供截割曲线。
[0033]2、能够进行局部动态更新，有效利用实时揭露的数据。随着工作面回采，不断揭露新数据，将揭露的数据一定范围内进行动态更新，得到更加精确的模型。
附图说明
[0034]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本专利技术的一些实施案例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0035]图1是网格a的位置示意图，其中，(a)中间网格；(b)边网格；(c)角网格；
[0036]图2是通过本专利技术的方法生成的工作面的效果图；
[0037]图3是本专利技术在工作面进行动态更新后的效果图。
具体实施方式
[0038]本专利技术给出的面向智能开采的地层界面插值建模方法具体包括以下步骤：
[0039]步骤1：将工作面测量获取的数据进行整理，形成集合α，集合中包含三列数据(x,y,z)，这三列数据(x,y,z)为工作面测量写实得到的煤层的坐标数据。其中，x表示横坐标、y表示纵坐标(x、y常用经纬度表示)，z表示高程；
[0040]表1集合α中点的格式
[0041]序号xyz1x1y1z12x2y2z23x3y3z3............
[0042]步骤2：将集合α进行网格化处理形成网格点的集合β1；
[0043]具体的，使步骤1得到的离散数据转换等间距的网格数据，能够实现离散数据网格化的方法有很多，其中包括函数拟合计算、克里金插值等。其中克里金插值可用PyKrige工具箱。
[0044]步骤3：根据集合β1的x、y坐标，得到插入点集φ1中插入点的坐标；
[0045]x
n
＝(x
本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种面向智能开采的地层界面插值建模方法，其特征在于，具体包括以下步骤：步骤1：将工作面测量获取的数据进行整理，形成集合α，集合中包含三列数据(x,y,z)，x表示横坐标、y表示纵坐标，z表示高程；步骤2：将集合α进行网格化处理形成网格点的集合β1；步骤3：根据集合β1的x、y坐标，得到插入点集φ1中插入点的坐标；步骤4：设置点集φ1中每个点的初始z坐标；具体的，将集合α的平均z坐标赋值给点集φ1中各点的初始z坐标；步骤5：设置点集φ1和集合β1中每个点的权重大小；步骤6：定义点集φ1和集合β1中各点的势η
a
；步骤7：令η
a
＝0，根据步骤6的定义，计算得到当前的点集φ1和集合β1中各网格点的z
a
值，得到点集φ1和点集β1的并集步骤8，循环执行步骤7，循环次数不小于10次，得到当前的点集φ1和点集β1的并集步骤9：采用步骤2得到的集合β1替换步骤8中得到的当前的点集φ1和点集β1的并集中的集合β1，得到新的并集步骤10：计算步骤9得到的并集的网格间距，判断网格间距是否达到目标网格间距，若已达到，则输出并集作为最终的插值结果；否则将并集作为集合β1，返回步骤3。2.如权利要求1所述的面向智能开采的地层界面插值建模方法，其特征在于，所述步骤2中的网格化处理采用函数拟合计算或者克里金插值。3.如权利要求1所述的面向智能开采的地层界面插值建模方法，其特征在于，所述步骤3中，插入点的坐标的计...

【专利技术属性】
技术研发人员：安林，刘再斌，王少龙，李鹏，王新苗，刘晨光，高珺，
申请(专利权)人：中煤科工西安研究院集团有限公司，
类型：发明
国别省市：