一种面向地表沉陷预测的凹多边形倾斜采煤工作面剖分方法技术

本发明专利技术公开一种面向地表沉陷预测的凹多边形倾斜采煤工作面剖分方法，适用于凹多边形倾斜煤矿开采评估领域。获取待预测工作面顶点坐标数据，在多边形边线及内部以固定间隔内插点，建立两个覆盖多边形工作面的大三角形，依据三角形外接圆准则将内插点逐个插入到大三角形中，将所有数据点构建剖分三角形，接着采集各剖分三角形重心坐标，判断各重心是否在多边形工作面边线或内部，以是留否去原则得到凹多边形采煤工作面所有剖分三角形。考虑到倾斜煤层不同剖分单元对应煤层参数不同。其基于内插点的构网算法理论严密，几何逻辑恰当，适用于凹多边形倾斜采煤工作面，剖分精度完善，有助于提高凹多边形倾斜采煤工作面地表移动变形的预测精度。形的预测精度。形的预测精度。

【技术实现步骤摘要】
一种面向地表沉陷预测的凹多边形倾斜采煤工作面剖分方法


[0001]本专利技术涉及一种面向地表沉陷预测的凹多边形倾斜采煤工作面剖分方法，属于多边形倾斜采煤工作面开采地表移动与变形预测

技术背景
[0002]煤矿开采打破原有地下岩层的应力平衡，采空存在使得应力传达到地表产生地表沉陷，威胁地面耕地及建构筑物安全和缩小使用寿命。在煤矿开采前进行采煤工作引起地表移动变形值的预测，可以提前采取地表次生灾害影响的预防和规划措施，实现采煤和环境保护共赢。
[0003]目前围绕煤炭开采地表沉陷预测方法进行了大量的研究工作，建立了满足常规地质采矿条件的地表沉陷预测方法，尤其是建立的概率积分法已成为现今最成熟的、应用最为广泛的矿山开采沉陷预计方法。在应用概率积分法预测采煤地表移动变形值时，因考虑到工作面尺寸较大而往往采取对工作面进行剖分处理的操作以提高预测精度。常规开采工作面多呈条状矩形，利用矩形剖分可以实现工作面的细小单元划分，而不规则工作面在应用矩形剖分时出现的边缘锯齿使得工作面剖分不完善，预测精度难以满足实践需求。由此三角剖分被引入到工作面剖分处理。三角剖分可以实现在工作面内插点的情况下将不规则工作面划分成细小单元三角形，其形状的多变可以实现与不规则多边形边缘的贴合，减少边缘精度的缺失。然而三角剖分算法应用到凹多边形采煤工作面时使得多边形凹陷区域产生了非法三角形，这一现象使得三角剖分算法应用到任意多边形剖分处理领域的适用性大大的降低，而如何解决凹多边形三角剖分处理产生非法三角形成为亟待解决的问题，对于提高煤矿开采地表沉陷预测精度具有重要的意义。

技术实现思路

[0004]针对上述技术不足之处，提供一种面向地表沉陷预测的凹多边形倾斜采煤工作面剖分方法，该方法理论严密，几何逻辑恰当，编程实现时间复杂度低，剖分精度完善，有助于提高不规则采煤工作面地表移动变形的预测精度。
[0005]为了实现上述技术目的，本专利技术的一种面向地表沉陷预测的凹多边形倾斜采煤工作面剖分方法，针对平面图为凹多边形的倾斜采煤工作面，步骤如下：
[0006]通过已有矿区开采资料获取该凹多边形倾斜煤层信息，包括开采工作面多边形空间的各顶点坐标及各顶点开采深度，即采深，将顶点坐标以工作面平面图逆时针分布的顺序排列并存储；
[0007]默认倾斜煤层的开采厚度均一，将多边形工作面简化为平面图，依据先边线后内部原则在多边形工作面的平面图中边线及内部标注内插点，内插点的间距根据工作经验或者固定比例进行，所述的固定比例为工作面上顶板至地表的深度即开采深度的五分之一；
[0008]标注所有内插点的坐标，然后将所有内插点坐标依X坐标从小至大原则排序存储入数组Vertices中，数组Vertices记录各内插点的坐标值(X,Y)；
[0009]提取凹多边形采煤工作面顶点坐标中最大的X
max
与最大的Y
max
坐标、最小的X
min
与Y
min
坐标，建立两个覆盖凹多边形的采煤工作面平面图的大三角形，并存储入三角形列表Triangles中，列表Triangles记录各三角形的三个顶点坐标；
[0010]对凹多边形采煤工作面的平面图进行三角剖分，具体过程如下：
[0011]遍历数组Vertices中存储的所有内插点坐标，若内插点坐标位于三角形列表Triangles中某三角形的外接圆内或外接圆上时，则将该三角形从列表中删除，然后将当前点连接刚刚删除的三角形的三个顶点，从而形成三个剖分三角形，并将三个剖分三角形的坐标加入三角形列表Triangles中；若该内插点坐标位于三角形列表Triangles中某三角形的外接圆外时，则直接将该三角形的坐标加入三角形列表Triangles中；
[0012]遍历完成数组Vertices中记载的所有内插点坐标之后，将两个大三角形的坐标从三角形列表Triangles中删除后余下的三角形即为所有初步三角剖分的剖分三角形；
[0013]由于凹多边形凹陷处初步剖分完成后存在非法三角形，因此再利用三角形重心定理计算得到各剖分三角形的重心坐标，若重心位于工作面坐标范围内部，则保留该三角形；若重心位于多边形工作面坐标范围外部或边线，则去除该三角形；最后得到不规则工作面所有剖分三角形。
[0014]进一步，由于倾斜煤层倾角的存在，后续地表移动变形预测过程使用平均采深将使得预测值脱离实际；根据各剖分三角形重心至多边形工作面各顶点距离作为权重，多边形工作面各顶点采深作为计算值，加权平均计算即可得到各剖分三角形对应的采深。
[0015]创建类UDelTris，该类用以记录各剖分单元信息，包括三角形三个顶点坐标、三角形序号、三角形重心点坐标、该三角形对应的开采深度。
[0016]进一步，判断内插点M与三角形外接圆圆心位置关系步骤为：
[0017]设某三角形三个顶点坐标分别为：A(x1,y1),B(x2,y2),C(x3,y3)，则该三角形重心坐标为P(X/D,Y/D)，然后比较内插点M点到三角形重心坐标P点的距离MP与三角形的A点到三角形重心坐标P点的距离AP两者的大小，即可确定该内插点与三角形外接圆圆心位置关系；
[0018]其中：
[0019]D＝(x1*(y2‑
y3)+x2*(y3‑
y1)+x3*(y1‑
y2))*2，
[0020]X＝(x
12
+y
12
)*(y2‑
y3)+(x
22
+y
22
)*(y3‑
y1)+(x
32
+y
32
)*(y1‑
y2),
[0021]Y＝(x
12
+y
12
)*(x3‑
x2)+(x
22
+y
22
)*(x1‑
x3)+(x
32
+y
32
)*(x2‑
x1)；
[0022]若MP＝AP，则判断内插点M在该三角形外接圆上；若MP>AP，则判断内插点M在该三角形外接圆外；若MP<AP，则判断内插点M在该三角形外接圆内。
[0023]进一步，计算各剖分三角形重心坐标方法为：某三角形三个顶点坐标分别为：A(x1,y1),B(x2,y2),C(x3,y3)，则该三角形重心坐标为P((x1+x2+x3)/3,(y1+y2+y3)/3)。
[0024]进一步，判断三角形重心点与凹多边形工作面位置关系步骤为：记某三角形重心点坐标为P，从P点作水平向左或向右的射线，若该射线与多边形工作面交点个数为奇数，则P点位于多边形工作面内部，若该射线与多边形交点个数为偶数或0，则P点位于多边形工作面外部。
[0025]进一步，计算各剖分三角形对应的采深的方法为：设多边形各顶点采深为H1、H2、
…
、H...

【技术保护点】

【技术特征摘要】
1.一种面向地表沉陷预测的凹多边形倾斜采煤工作面剖分方法，其特征在于：针对平面图为凹多边形的倾斜采煤工作面，步骤如下：通过已有矿区开采资料获取该凹多边形倾斜煤层信息，包括开采工作面多边形空间的各顶点坐标及各顶点开采深度，即采深，将顶点坐标以工作面平面图逆时针分布的顺序排列并存储；默认倾斜煤层的开采厚度均一，将多边形工作面简化为平面图，依据先边线后内部原则在多边形工作面的平面图中边线及内部标注内插点，内插点的间距根据工作经验或者固定比例进行，所述的固定比例为工作面上顶板至地表的深度即开采深度的五分之一；标注所有内插点的坐标，然后将所有内插点坐标依X坐标从小至大原则排序存储入数组Vertices中，数组Vertices记录各内插点的坐标值(X,Y)；提取凹多边形采煤工作面顶点坐标中最大的X
max
与最大的Y
max
坐标、最小的X
min
与Y
min
坐标，建立两个覆盖凹多边形的采煤工作面平面图的大三角形，并存储入三角形列表Triangles中，列表Triangles记录各三角形的三个顶点坐标；对凹多边形采煤工作面的平面图进行三角剖分，具体过程如下：遍历数组Vertices中存储的所有内插点坐标，若内插点坐标位于三角形列表Triangles中某三角形的外接圆内或外接圆上时，则将该三角形从列表中删除，然后将当前点连接刚刚删除的三角形的三个顶点，从而形成三个剖分三角形，并将三个剖分三角形的坐标加入三角形列表Triangles中；若该内插点坐标位于三角形列表Triangles中某三角形的外接圆外时，则直接将该三角形的坐标加入三角形列表Triangles中；遍历完成数组Vertices中记载的所有内插点坐标之后，将两个大三角形的坐标从三角形列表Triangles中删除后余下的三角形即为所有初步三角剖分的剖分三角形；由于凹多边形凹陷处初步剖分完成后存在非法三角形，因此再利用三角形重心定理计算得到各剖分三角形的重心坐标，若重心位于工作面坐标范围内部，则保留该三角形；若重心位于多边形工作面坐标范围外部或边线，则去除该三角形；最后得到不规则工作面所有剖分三角形。2.根据权利要求1所述的一种面向地表沉陷预测的凹多边形倾斜采煤工作面剖分方法，其特征在于，由于倾斜煤层倾角的存在，后续地表移动变形预测过程使用平均采深将使得预测值脱离实际；根据各剖分三角形重心至多边形工作面各顶点距离作为权重，多边形工作面各顶点采深作为计算值，加权平均计算即可得到各剖分三角形对应的采深；考虑到倾斜煤层不同剖分单元对应煤层参数不同，创建UDelTris类整合三角剖分单元对应坐标及参数信息。3.根据权利要求1所述的一种面向地表沉陷预测...

【专利技术属性】
技术研发人员：戴广礼，周华安，李怀展，黄建勇，孙兢超，唐路，李红佳，刘超，常友伟，霍文奇，
申请(专利权)人：中国矿业大学，
类型：发明
国别省市：