CN107886575A - 一种露天矿采场三角形网格裁剪煤层四边形网格的方法 - Google Patents
一种露天矿采场三角形网格裁剪煤层四边形网格的方法 Download PDFInfo
- Publication number
- CN107886575A CN107886575A CN201711120978.1A CN201711120978A CN107886575A CN 107886575 A CN107886575 A CN 107886575A CN 201711120978 A CN201711120978 A CN 201711120978A CN 107886575 A CN107886575 A CN 107886575A
- Authority
- CN
- China
- Prior art keywords
- pol
- polygon
- intersection
- chained list
- mesh
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
- G06T17/205—Re-meshing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/32—Image data format
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明提供一种露天矿采场三角形网格裁剪煤层四边形网格的方法，涉及煤层地质技术领域。该方法首先生成煤层顶底板界面的四边形网格和露天矿采场的三角形网格，然后生成由三角形网格和四边形网格相交得到的交线多边形，并对其按逆时针排序，得到裁剪多边形，最后使用裁剪多边形对煤层顶底板界面的四边形网格进行裁剪，实现对煤层地质模型的更新。本发明提供的露天矿采场三角形网格裁剪煤层四边形网格的方法，有效解决了四边形曲面与三角形面的交线是空间曲线的问题，使得实际应用得以实现，这种直接利用三角形网格对四边形网格进行裁剪的方法，以不增加存储量为前提，高效地实现了三角形网格对四边形网格的裁剪。
Description
技术领域
本发明涉及煤层地质技术领域，尤其涉及一种露天矿采场三角形网格裁剪煤层四边形网格的方法。
背景技术
煤矿床是由一系列在时间上连续或间断沉积的沉积岩组成，是典型的沉积矿床。地质界面是在规定的范围内划分规则的网格，以勘探钻孔数据作为样本，经过插值而形成的。煤矿床中除了岩层接触面外，还有一类重要的界面就是地表面，地表面在矿山未开采前，指的是原始地表面，在开采后，由于原始地表面被采矿工程所破坏，这时的地表面是由未遭破坏的原始地表面部分和采场表面组成。露天矿的地表面有两类：其一是自然地表面，即原始地表面，是自然形成的；其二是工程地表面，是人工形成的，主要有采剥工程开挖形成的采剥工程面和排土工程堆筑形成的排土工程面。但无论是自然形成的还是人工形成的，地表面都存在一个共同特点，即裸露于地上或曾经裸露于地上，其上的样本数据一般通过地表测量获得，取样点密度高，且常常充分考虑地貌特点因此，地表面上的测量数据既有点特征数据，又有线特征数据，构造的地表面模型不仅要反映点状特征信息，还要反映线状特征信息，规则网格模型不能对地表面上的点和线特征给予显式的表达，而对于工程地表面来讲是非常重要的。而不规则网格模型是沿着所提供的地学曲面上的特征点和特征线划分网格的，可以对地表面上的点和线特征数据给予明确、显式的表达。因此，地表面模型一般均采用不规则三角网格模型。
随着露天矿开采的进行，煤层和岩层被不断开采出去，形成采场表面与地质界面相交。为了获得更新的矿床模型，需要利用地表面模型对地质界面模型进行裁剪，即地表面三角形网格裁剪煤层顶底板界面四边形网格。但由于四边形网格是曲面，它与三角形面的交线是空间曲线，在实际应用中难以实现，现有的裁剪方法是将四边形网格转换成三角形网格进行处理，这些方法既增加了存储量，也降低了效率。
发明内容
针对现有技术的缺陷，本发明提供一种露天矿采场三角形网格裁剪煤层四边形网格的方法，利用露天矿采场的三角形网格对煤层的顶、底板的四边形网格进行裁切，实现露天矿矿床地质模型的更新。
一种露天矿采场三角形网格裁剪煤层四边形网格的方法，包括以下步骤：
步骤1：生成煤层界面的四边形网格，具体方法为：
在煤层界面范围内，按照一定的间距划分网格，对煤层顶底板界面网格的节点的高程进行插值，插值点为P(x,y)，共有n个样本点Pi(xi,yi)，其中，i＝1,2,…,n，其样本值为zi，采用加权最小二乘拟合法计算插值点P的高程值z＝f(x,y)，其中，f(x,y)为加权最小二乘拟合曲面，生成煤层界面的四边形网格；
生成的煤层界面的四边形网格，由如下所示的四边形链表存取：
pQu＝{pQu0(v00,v01,v02,v03)→…→pQuk(vk0,vk1,vk2,vk3)→…→pQum(vm0,vm1,vm2,vm3)}
其中，pQuk(vk0,vk1,vk2,vk3)，为空间四边形网格中的第k个四边形，vk0,vk1,vk2,vk3分别为第k个四边形的四个顶点，k＝0、1、…m，m+1为煤层界面生成的四边形网格总数；
链表pQu中的四边形的顶点都按逆时针排列；
步骤2：根据露天矿采场数据的特点，建立露天矿采场的三角形网格模型，具体方法为：
步骤2.1：根据露天矿采场数据的来源，将露天矿采场数据分为两部分，一部分为原地表数据，另一部分为露天矿采场数据；原地表数据包括地表测点、地形等高线及地貌特性线数据；露天矿采场数据包括台阶点、台阶线和采场散点数据；根据露天矿采场数据的几何形态，将这些露天矿采场数据分为点状数据和线状数据两类；
步骤2.2：获取线状数据的顶点，将其与点状数据一起形成进行三角剖分所需的数据点；
步骤2.3：对步骤2.2形成的数据点进行Delaunay三角剖分，生成剖分三角网；
步骤2.4：将地形等高线、地貌特性线和台阶线这些线状数据作为约束边插入到步骤2.3生成的三角网中，实现对露天矿采场数据的约束三角剖分(Constrained DelaunayTriangulation,CDT)，形成露天矿采场的三角网模型，由如下所示的三角形链表存取：
pTr＝{pTr0(p00,p01,p02)→…→pTrl(pl0,pl1,pl2)→…→pTrs(ps0,ps1,ps2)}
其中，pTrl(pl0,pl1,pl2)为空间三角形网格中的第l个三角形，pl0,pl1,pl2分别为第l个三角形的三个顶点，l＝0、1、…s，s+1为生成的三角网模型中三角形的总个数；
链表pTr中的三角形的顶点按逆时针排列；
步骤3：通过使用以直线代替曲线求解空间三角形与空间四边形相交的交线的方法，得到露天矿采场形成的不规则的三角形空间网格与煤层顶底板界面形成的规则的空间四边形网格相交构成的多边形，具体方法为：
步骤3.1：选取三角形链表pTr中一个三角形，判断该三角形是否有两条直线段与煤层界面的四边形网格相交，如果有则执行步骤3.2，否则重新选取三角形链表pTr中其他的三角形进行判断；
步骤3.2：将三角形链表中选取的三角形与煤层界面的四边形的两个交点的连线作为该三角形与四边形网格的交线；
步骤3.3：依次求出三角形链表中所有三角形与四边形链表中各四边形网格的交线，建立交线链表，如下所示：
pL＝{pL0(l00,l01)→pL1(l10,l11)…→pLr(lr0,lr1)→…→pLo(lo0,lo1)}
其中，pLr(lr0,lr1)为第r条交线，lr0,lr1分别为第r条交线的两个端点，r＝0、1、…、o，o+1为交线总数，交线链表中的所有交线组成一个闭合的多边形；
步骤4：使用交线构成的闭合多边形对煤层顶底板界面的四边形网格进行裁剪，具体方法为：
步骤4.1：将闭合的多边形按逆时针排序；
交线链表pL是一个线段间无序连接的多边形，将这些无序线段重新按逆时针进行排序的具体方法为：
步骤4.1.1：将交线多边形链表pL中的所有线段按首、尾端点相连顺序排列，形成一个顺序连接的交线多边形链表pOL，具体方法为：
将交线多边形链表pL的首结点pL0(l00,l01)作为pOL的首结点，即pOL0(l00,l01)；根据端点l01的坐标值，在交线链表pL其余的交线中找到与端点l01有相同端点的线段，建立首结点的相邻节点pOL1(l10,l11)，使得这两个节点满足端点l01＝l10；重复同样的操作，将各端点依次相连，构成一个顺序多边形链表，如下所示：
pOL＝{pOL0(l00,l01)→pOL1(l10,l11)…→pOLr′(lr′0,lr′1)→…→pOLo(lo0,lo1)}
其中，pOLr′(lr0,lr′1)为顺序连接的交线多边形的第r′条线段，lr′0,lr′1分别为第r′条线段的两个端点,且lr′1＝l(r′+1)0，r′＝0、1、…o；
步骤4.1.2：遍历交线多边形链表pOL，得到关于x,y坐标的最小和最大值端点的顺序号Nxmin、Nxmax、Nymin和Nymax，如果这四个顺序号的大小关系满足下列四个大小关系式中任意一个，则该交线多边形为逆时针排列的多边形，执行步骤4.2，否则，该交线多边形为顺时针排列的多边形，执行步骤4.1.3；
步骤4.1.3：对顺时针排列的交线多边形重新排序，得到逆时针排序的交线多边形，具体方法为：
将多边形链表pOL进行反向，并将线段的首点和末点互换，即完全将顺时针多边形链表pOL转变成按逆时针排列的交线多边形链表；
步骤4.2：求逆时针排列的交线多边形与顶底板四边形网格的交点，构成裁剪多边形链表；
求交线多边形链表pOL中的每条线段与四边形链表pQu中的每个四边形的边的交点，如果交线多边形链表pOL中第r′条线段pOLr′(lr′0,lr′1)与四边形链表pQu中的某个四边形的边相交，交点是ptr′，则ptr′将线段pOLr′(lr′0,lr′1)分成两条线段pOLr′(lr′0,ptr′)和pOLr′+1(ptr′,lr′1)，在交线多边形链表pOL中，使用新生成的线段pOLr′(lr′0,ptr′)替换线段pOLr′(lr′0,lr′1)，并将另一条新生成的线段pOLr′+1(ptr′,lr′1)插入到pOLr′(lr′0,ptr′)之后，更新后的交线多边形链表pOL′如下所示：
pOL′＝{pOL0(l00,l01)→pOL1(l10,l11)…→pOLr′(lr′0,ptr′)→pOLr′+1(ptr′,lr′1)→…→pOLo(lo0,lo1)}
得到的多边形链表pOL′为裁剪多边形链表；
步骤4.3：用裁剪多边形对煤层界面四边形网格进行裁剪，实现对露天矿矿床地质模型的更新，具体方法为：
步骤4.3.1：根据空间关系，裁剪多边形将煤层顶底板界面四边形网格分成三部分，即外部四边形网格、内部四边形网格和相交四边形网格；
步骤4.3.2：将外部四边形网格和相交四边形网格转存入新的网格单元链表中，实现去掉内部四边形网格；
步骤4.3.3：遍历所有相交四边形网格，分别对每个相交四边形网格进行处理，生成不规则的多边形，具体处理方法为：
第一步，选取位于相交四边形内部的裁剪多边形的线段；
第二步，选取第一步所选取四边形的位于裁剪多边形外部的顶点，并与第一步得到的线段组成逆时针方向的多边形，由这些多边形组成裁剪后的不规则网格；
步骤4.3.4：由生成的不规则的多边形网格代替相应的相交四边形网格，完成对相交四边形网格的处理，实现露天矿采场三角形网格对所有煤层顶底板四边形网格的裁剪，从而实现对露天矿矿床地质模型的更新。
由上述技术方案可知，本发明的有益效果在于：本发明提供的一种露天矿采场三角形网格裁剪煤层四边形网格的方法，使用通过三角剖分得到的露天矿采场三角形网格和插值得到的煤层顶底板界面四边形网格在空间相交的交线多边形，对煤层顶底板界面网格进行裁剪，实现对煤层界面模型进行更新。该方法有效解决了四边形曲面与三角形面的交线是空间曲线的问题，使得实际应用得以实现，这种直接利用三角形网格对四边形网格进行裁剪的方法，以不增加存储量为前提，高效地实现了三角形网格对四边形网格的裁剪。
附图说明
图1为本发明实施例提供的采场表面与地质界面相交的示意图；
图2为本发明实施例提供的一种露天矿采场三角形网格裁剪煤层四边形网格的方法的流程图；
图3为本发明实施例提供的样本钻孔的分布俯视示意图；
图4为本发明实施例提供的16#煤的顶板界面四边形网格的示意图；
图5为本发明实施例提供的露天矿采场的测量示意图；
图6为本发明实施例提供的露天矿采场约束三角剖分示意图；
图7为本发明实施例提供的空间三角形与四边形相交的示意图；
图8为本发明实施例提供的交线闭合多边形示意图；
图9为本发明实施例提供的对裁剪多边形进行逆时针排序的示意图；
图10为本发明实施例提供的经过裁剪后的采场表面与地质界面相交的示意图。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
本发明实施例以如图1所示的某露天煤矿的采场勘测数据作为采场表面模型的数据源，以钻孔的位置坐标(x,y)、煤层的底板标高z作为煤层界面的插值样本数据，通过本发明方法实现露天矿采场的三角形网格对煤层界面的四边形网格裁切。
一种露天矿采场三角形网格裁剪煤层四边形网格的方法，如图2所示，包括以下步骤：
步骤1：生成煤层界面的四边形网格，具体方法为：
在煤层界面范围内，按照一定的间距划分网格，对煤层顶底板界面网格的节点的高程进行插值，插值点为P(x,y)，共有n个样本点Pi(xi,yi)，其中，i＝1,2,…,n，其样本值为zi，采用加权最小二乘拟合法计算插值点P的高程值z＝f(x,y)，其中，f(x,y)为加权最小二乘拟合曲面，生成煤层界面的四边形网格；
生成的煤层界面的四边形网格，由如下所示的四边形链表存取：
pQu＝{pQu0(v00,v01,v02,v03)→…→pQuk(vk0,vk1,vk2,vk3)→…→pQum(vm0,vm1,vm2,vm3)}
其中，pQuk(vk0,vk1,vk2,vk3)，为空间四边形网格中的第k个四边形，vk0,vk1,vk2,vk3分别为第k个四边形的四个顶点，k＝0、1、…m，m+1为煤层界面生成的四边形网格总数；
链表pQu中的四边形的顶点都按逆时针排列；
所述采用加权最小二乘拟合法进行计算的具体计算方法为：
首先选择如下式所求的二次多项式曲面作为拟合曲面，
f(x,y)＝c1+c2x+c3y+c4xy
其中，c1,c2,c3,c4均为曲面参数；
将满足公式的二次曲面f(x,y)称为最小二乘拟合曲面；
在上述最小二乘拟合曲面中引入一个距离权，如下公式所示：
其中，di为第i个样本点Pi到插值点的距离，m′为距离的幂指数；
实现对拟合曲面的控制作用，得到二次曲面f(x,y)在插值点P的加权最小二乘拟合曲面，如下式所示：
根据最小二乘原理，得到其中，j＝1,2,3,4；
进而得到加权最小二乘拟合曲面f(x,y)的四个系数为变元的四元一次方程组，解此方程组，获得加权最小二乘拟合曲面f(x,y)，即得到插值点P(x,y)的高程值z＝f(x,y)；
本实施例以16#煤的顶板数据为样本，共有638个钻孔为样本数据，部分钻孔样本数据如表1所示，样本钻孔的分布俯视图如图3所示，X轴和Y轴分别表示钻孔数据水平方向的数值和竖直方向的数值。在矿山境界范围内，对这个范围区域按照30m间隔划分为单元规则网格，形成如图4所示的顶板界面四边形网格，为了显示方便，图中只显示了矿山境界范围内的四边形网格。
表1钻孔样本数据表(部分)
| 序号 | 钻孔号 | x | y | 16#顶板标高(m) |
| 1 | G21012 | -18746.94 | 85475.66 | 639.44 |
| 2 | G21022 | -18145.9 | 85115.86 | 632.94 |
| 3 | G21029 | -17866.29 | 84949.03 | 623.39 |
| 4 | G20110 | -18147.67 | 84813.04 | 619.8 |
| … | … | … | … | … |
| 638 | D01006 | -24649.15 | 80613.89 | 537.03 |
步骤2：根据露天矿采场数据的特点，建立露天矿采场的三角形网格模型，具体方法为：
步骤2.1：根据露天矿采场数据的来源，将露天矿采场数据分为两部分，一部分为原地表数据，另一部分为露天矿采场数据；原地表数据包括地表测点、地形等高线及一些地貌特性线数据；露天矿采场数据包括台阶点、台阶线和采场散点数据；根据露天矿采场数据的几何形态，将这些数据分为点状数据和线状数据两类；
步骤2.2：获取线状数据的顶点，将其与点状数据一起形成进行三角剖分所需的数据点；
步骤2.3：对步骤2.2形成的数据点进行Delaunay三角剖分，生成剖分三角网；
步骤2.4：将地形等高线、地貌特性线和台阶线这些线状数据作为约束边插入到步骤2.3生成的三角网中，实现对露天矿采场数据的约束三角剖分(Constrained DelaunayTriangulation,CDT)，形成露天矿采场的三角网模型，由如下所示的三角形链表存取：
pTr＝{pTr0(p00,p01,p02)→…→pTrl(pl0,pl1,pl2)→…→pTrs(ps0,ps1,ps2)}
其中，pTrl(pl0,pl1,pl2)为空间三角形网格中的第l个三角形，pl0,pl1,pl2分别为第l个三角形的三个顶点，l＝0、1、…s，s+1为生成的三角网模型中三角形的总个数；
链表pTr中的三角形的顶点按逆时针排列；
本实施例中，采集如图5所示的露天矿采场的点状数据和线状数据作为Delaunay三角剖分的输入数据，最终形成如图6所示的露天矿采场的三角网模型，即露天矿采场的三维几何模型，该三角网模型共生成1226个三角形。
步骤3：通过使用以直线代替曲线求解空间三角形与空间四边形相交的交线的方法得到露天矿采场形成的不规则的三角形空间网格与煤层顶底板界面形成的规则的空间四边形网格相交构成的多边形；
由于煤层顶底板的四边形网格采用的是线性插值，所以空间四边为双线性曲面，其曲面方程如下式所示：
z＝f(x,y)＝c1+c2x+c3y+c4xy
式中，曲面参数c1,c2,c3,c4由四边形的四个顶点确定；
露天矿采场空间三角形的空间平面方程，如下式所示：
z＝ax+by+c
式中，a、b、c为方程参数，由三角形的三个顶点确定；
空间三角形与空间四边形求交线，直接的方法是将上述空间四边形的曲面方程和空间三角形的空间平面方程联立，得到空间三角形与空间四边形相交的交线，两个联立的方程如下所示：
上述联立后的方程为一空间曲线方程，这种理论上的交线方程，在实际操作中很难求解，所以，以直线来近似代替曲线的方法求解空间三角形与空间四边形相交的交线，具体求解方法为：
步骤3.1：选取三角形链表pTr中一个三角形，判断该三角形是否有两条直线段与煤层界面的四边形网格相交，如果有则执行步骤3.2，否则重新选取三角形链表pTr中其他的三角形；
本实施例以图7为例，描述判断三角形中的直线段与四边形网格是否相交的方法：
在三角形T0T1T2中，直线段T0T1的方程为：
其中，x0、y0和z0分别为直线段中一个端点T0的坐标值，x1、y1和z1分别为直线段的另一个端点T1的坐标值，t为比例参数，0≤t≤1；
直线段T0T1的参数方程为：
令m1＝(x1-x0),n1＝(y1-y0),p＝(z1-z0)，将直线段的参数方程代入四边形的曲面方程，整理得如下方程：
c4m1n1t2+[c4(y0m1+x0n1)+c3n1+c2m1-p]t+(c4x0y0+c3y0+c2x0+c1-z0)＝0
将上式化简为：
at2+bt+c＝0
其中，a＝c4m1n1，b＝c4(y0m1+x0n1)+c3n1+c2m1-p，c＝c4x0y0+c3y0+c2x0+c1-z0；
根据该标准一元二次方程根的情况判断三角形的直线段与四边形网格是否相交，如果该方程至少有一个实数根，则三角形的直线段与四边形网格相交，否则三角形的直线段与四边形网格不相交；
步骤3.2：将三角形链表中三角形与煤层界面的四边形的两个交点的连线作为该三角形与四边形网格的交线；
步骤3.3：依次求出三角形链表中所有三角形与四边形链表中各四边形网格的交线，建立交线链表，如下所示：
pL＝{pL0(l00,l01)→pL1(l10,l11)…→pLr(lr0,lr1)→…→pLo(lo0,lo1)}
其中，pLr(lr0,lr1)为第r条交线，lr0,lr1分别为第r条交线的两个端点，r＝0、1、…、o，o+1为交线总数，交线链表中的所有交线组成一个闭合的多边形；
本实施例中，三角形链表中所有三角形与四边形链表中各四边形网格的交线所构成的闭合的多边形如图8所示，该闭合多边形由89条交线组成。
步骤4：使用交线构成的闭合多边形对煤层顶底板界面的四边形网格进行裁剪，具体方法为：
步骤4.1：将闭合的多边形按逆时针排序；
交线链表pL是一个线段间无序连接的多边形，将这些无序线段重新进行逆时针排序，具体方法为：
步骤4.1.1：将交线多边形链表pL中的所有线段按首、尾端点相连顺序排列，形成一个顺序连接的交线多边形链表pOL；
将交线多边形链表pL的首结点pL0(l00,l01)作为pOL的首结点，即pO0(l00,l01)；根据端点l01的坐标值，在交线链表pL其余的线段中找到与端点l01有相同端点的线段，建立首结点的相邻节点pO1(l10,l11)，使得这两个节点满足端点l01＝l10；重复同样的操作，将各端点依次相连，构成一个顺序多边形链表：
pOL＝{pOL0(l00,l01)→pOL1(l10,l11)…→pOLr′(lr′0,lr′1)→…→pOLo(lo0,lo1)}
其中，pOLr′(lr′0,lr′1)为顺序连接的交线多边形中的第r′条线段，r′＝0、1、…、o，lr′0,lr′1分别为线段的两个端点,且lr′1＝l(r′+1)0；
步骤4.1.2：按如图9所示遍历交线多边形链表pOL，得到关于x,y坐标的最小和最大值端点的顺序号Nxmin、Nxmax、Nymin和Nymax，如果这四个顺序号的大小关系满足下列四个大小关系式中任意一个，则该交线多边形为逆时针排列的多边形，执行步骤4.2，否则，该交线多边形为顺时针排列的多边形，执行步骤4.1.3；
步骤4.1.3：对顺时针排列的交线多边形重新排序，得到逆时针排序的交线多边形，具体方法为：
将多边形链表pOL进行反向，并将线段的首点和末点互换，即完全将顺时针多边形链表pOL转变成按逆时针排列的交线多边形链表。
步骤4.2：求逆时针排列的交线多边形与顶底板四边形网格的交点，构成裁剪多边形链表；
求交线多边形链表pOL中的每条线段与四边形链表pQu中的每个四边形的边的交点，如果交线多边形链表pOL中第r′条线段pOLr′(lr′0,lr′1)与四边形链表pQu中的某个四边形的边相交，交点是ptr′，则ptr′将线段pOLr′(lr′0,lr′1)分成两条线段pOLr′(lr′0,ptr′)和pOLr′+1(ptr′,lr′1)，在交线多边形链表pOL中，使用新生成的线段pOLr′(lr′0,ptr′)替换线段pOLr′(lr′0,lr′1)，并将另一条新生成的线段pOLr′+1(ptr′,lr′1)插入到pOLr′(lr′0,ptr′)之后，更新后的交线多边形链表pOL′如下所示：
pOL′＝{pOL0(l00,l01)→pOL1(l10,l11)…→pOLr′(lr′0,ptr′)→pOLr′+1(ptr′,lr′1)→…→pOLo(lo0,lo1)}得到的多边形链表pOL′为裁剪多边形链表；
步骤4.3：用裁剪多边形对煤层界面四边形网格进行裁剪，具体方法为：
步骤4.3.1：根据空间关系，裁剪多边形将煤层顶底板界面四边形网格分成三部分，即外部四边形网格、内部四边形网格和相交四边形网格；
步骤4.3.2：将外部四边形网格和相交四边形网格转存入新的网格单元链表中，实现去掉内部四边形网格；
步骤4.3.3：遍历所有相交四边形网格，分别对每个相交四边形网格进行处理，生成不规则的多边形，具体处理方法为：
第一步，选取位于相交四边形内部的裁剪多边形的线段；
第二步，选取第一步所选取四边形的位于裁剪多边形外部的顶点，并与第一步得到的线段组成逆时针方向的多边形，由这些多边形组成裁剪后的不规则网格；
步骤4.3.4：由生成的不规则的多边形网格单元代替相应的相交四边形网格单元，完成对相交四边形网格的处理，实现露天矿采场三角形网格对所有煤层顶底板四边形网格的裁剪，从而实现对露天矿矿床地质模型的更新。
本实施的最终裁剪结果如图10所示，采场三角形网内的四边形网格已经全部裁剪掉。通过相同的方法，可以对矿床中其他地质界面进行裁剪，最终实现对整个矿床模型的更新。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明权利要求所限定的范围。
Claims (6)
1.一种露天矿采场三角形网格裁剪煤层四边形网格的方法，其特征在于：包括以下步骤：
步骤1：生成煤层界面的四边形网格，具体方法为：
在煤层界面范围内，按照一定的间距划分网格，对煤层顶底板界面网格的节点的高程进行插值，插值点为P(x，y)，共有n个样本点Pi(xi，yi)，其中，i＝1，2，...，n，其样本值为zi，采用加权最小二乘拟合法计算插值点P的高程值z＝f(x，y)，其中，f(x，y)为加权最小二乘拟合曲面，生成煤层界面的四边形网格；
生成的煤层界面的四边形网格，由如下所示的四边形链表存取：
pQu＝{pQu0(v00，v01，v02，v03)→…→pQuk(vk0，vk1，vk2，vk3)→…→pQum(vm0，Vm1，Vm2，Vm3)}
其中，pQuk(Vk0，vk1，vk2，vk3)，为空间四边形网格中的第k个四边形，vk0，vk1，vk2，vk3分别为第k个四边形的四个顶点，k＝0、1、…m，m+1为煤层界面生成的四边形网格总数；
链表pQu中的四边形的顶点都按逆时针排列；
步骤2：根据露天矿采场数据的特点，建立露天矿采场的三角形网格模型，具体方法为：
步骤2.1：根据露天矿采场数据的几何形态，将露天矿采场数据分为点状数据和线状数据两类；
步骤2.2：获取线状数据的顶点，将其与点状数据一起形成进行三角剖分所需的数据点；
步骤2.3：对步骤2.2形成的数据点进行Delaunay三角剖分，生成剖分三角网；
步骤2.4：将线状数据作为约束边插入到步骤2.3生成的三角网中，实现对露天矿采场数据的约束三角剖分(Constrained Delaunay Triangulation，CDT)，形成露天矿采场的三角网模型，由如下所示的三角形链表存取：
pTr＝{pTr0(p00，p01，p02)→…→pTrl(pl0，pl1，pl2)→…→pTrs(ps0，ps1，ps2)}
其中，pTrl(pl0，pl1，pl2)为空间三角形网格中的第l个三角形，pl0，pl1，pl2分别为第l个三角形的三个顶点，l＝0、1、…s，s+1为生成的三角网模型中三角形的总个数；
链表pTr中的三角形的顶点按逆时针排列；
步骤3：通过使用以直线代替曲线求解空间三角形与空间四边形相交的交线的方法，得到露天矿采场形成的不规则的三角形空间网格与煤层顶底板界面形成的规则的空间四边形网格相交构成的多边形，具体方法为：
步骤3.1：选取三角形链表pTr中一个三角形，判断该三角形是否有两条直线段与煤层界面的四边形网格相交，如果有则执行步骤3.2，否则重新选取三角形链表pTr中其他的三角形进行判断；
步骤3.2：将三角形链表中选取的三角形与煤层界面的四边形的两个交点的连线作为该三角形与四边形网格的交线；
步骤3.3：依次求出三角形链表中所有三角形与四边形链表中各四边形网格的交线，建立交线链表，如下所示：
pL＝{pL0(l00，l01)→pL1(l10，l11)…→pLr(lr0，lr1)→…→pLo(lo0，lo1)}
其中，pLr(lr0，lr1)为第r条交线，lr0，lr1分别为第r条交线的两个端点，r＝0、1、…、o，o+1为交线总数，交线链表中的所有交线组成一个闭合的多边形；
步骤4：使用交线构成的闭合多边形对煤层顶底板界面的四边形网格进行裁剪，具体方法为：
步骤4.1：交线链表pL是一个线段间无序连接的多边形，将这些无序线段重新按逆时针进行排序，得到逆时针排序的交线多边形链表pOL；
步骤4.2：求逆时针排列的交线多边形与顶底板四边形网格的交点，构成裁剪多边形链表pOL′；
步骤4.3：用裁剪多边形对煤层界面四边形网格进行裁剪，实现对露天矿矿床地质模型的更新。
2.根据权利要求1所述的一种露天矿采场三角形网格裁剪煤层四边形网格的方法，其特征在于：步骤2.1所述点状数据包括地表测点、台阶平盘点、台阶点，所述线状数据包括地形等高线、地貌特征线、台阶线。
3.根据权利要求2所述的一种露天矿采场三角形网格裁剪煤层四边形网格的方法，其特征在于：步骤4.1所述将交线构成的闭合多边形按逆时针排序的具体方法为：
步骤4.1.1：将交线多边形链表pL中的所有线段按首、尾端点相连顺序排列，形成一个顺序连接的交线多边形链表pOL，具体方法为：
将交线多边形链表pL的首结点pL0(l00，l01)作为pOL的首结点，即pOL0(l00，l01)；根据端点l01的坐标值，在交线链表pL其余的交线中找到与端点l01有相同端点的线段，建立首结点的相邻节点pOL1(l10，l11)，使得这两个节点满足端点l01＝l10；重复同样的操作，将各端点依次相连，构成一个顺序多边形链表，如下所示：
pOL＝{pOL0(l00，l01)→pOL1(l10，l11)…→pOLr′(lr′0，lr′1)→…→pOLo(lo0，lo1)}
其中，pOLr′(lr0，lr′1)为顺序连接的交线多边形的第r′条线段，lr′0，lr′1分别为第r′条线段的两个端点，且lr′1＝l(r′+1)0，r′＝0、1、…o；
步骤4.1.2：遍历交线多边形链表pOL，得到关于x，y坐标的最小和最大值端点的顺序号Nxmin、Nxmax、Nymin和Nymax，如果这四个顺序号的大小关系满足下列四个大小关系式中任意一个，则该交线多边形为逆时针排列的多边形，执行步骤4.2，否则，该交线多边形为顺时针排列的多边形，执行步骤4.1.3；
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mi>N</mi>
<mi>y</mi>
<mi> </mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo><</mo>
<mi>N</mi>
<mi>x</mi>
<mi> </mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo><</mo>
<mi>N</mi>
<mi>y</mi>
<mi> </mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo><</mo>
<mi>N</mi>
<mi>x</mi>
<mi> </mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>N</mi>
<mi>x</mi>
<mi> </mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo><</mo>
<mi>N</mi>
<mi>y</mi>
<mi> </mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo><</mo>
<mi>N</mi>
<mi>x</mi>
<mi> </mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo><</mo>
<mi>N</mi>
<mi>y</mi>
<mi> </mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>N</mi>
<mi>y</mi>
<mi> </mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo><</mo>
<mi>N</mi>
<mi>x</mi>
<mi> </mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo><</mo>
<mi>N</mi>
<mi>y</mi>
<mi> </mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo><</mo>
<mi>N</mi>
<mi>x</mi>
<mi> </mi>
<mi>max</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>N</mi>
<mi>x</mi>
<mi> </mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo><</mo>
<mi>N</mi>
<mi>y</mi>
<mi> </mi>
<mi>min</mi>
<mo><</mo>
<mi>N</mi>
<mi>x</mi>
<mi> </mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo><</mo>
<mi>N</mi>
<mi>y</mi>
<mi> </mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
步骤4.1.3：对顺时针排列的交线多边形重新排序，得到逆时针排序的交线多边形，具体方法为：
将多边形链表pOL进行反向，并将线段的首点和末点互换，即完全将顺时针多边形链表pOL转变成按逆时针排列的交线多边形链表。
4.根据权利要求3所述的一种露天矿采场三角形网格裁剪煤层四边形网格的方法，其特征在于：所述步骤4.2的具体方法为：
求交线多边形链表pOL中的每条线段与四边形链表pQu中的每个四边形的边的交点，如果交线多边形链表pOL中第r′条线段pOLr′(lr′0，lr′1)，与四边形链表pQu中的某个四边形的边相交，交点是ptr′，则ptr′将线段pOLr′(lr′0，lr′1)分成两条线段pOLr′(lr′0，ptr′)和pOLr′+1(ptr′，lr′1)，在交线多边形链表pOL中，使用新生成的线段pOLr′(lr′0，ptr′)替换线段pOLr′(lr′0，lr′1)，并将另一条新生成的线段pOLr′+1(ptr′，lr′1)插入到pOLr′(lr′0，ptr′)之后，更新后的交线多边形链表pOL′如下所示：
pOL′＝{pOL0(l00，l01)→pOL1(l10，l11)…→pOLr′(lr′0，ptr′)→pOLr′+1(ptr′，lr′1)→…→pOLo(lo0lo1)}
其中，r′＝0、1、…、o，得到的多边形链表pOL′为裁剪多边形链表。
5.根据权利要求3所述的一种露天矿采场三角形网格裁剪煤层四边形网格的方法，其特征在于：所述步骤4.3的具体方法为：
步骤4.3.1：根据空间关系，裁剪多边形将煤层顶底板界面四边形网格分成三部分，即外部四边形网格、内部四边形网格和相交四边形网格；
步骤4.3.2：将外部四边形网格和相交四边形网格转存入新的网格单元链表中，实现去掉内部四边形网格；
步骤4.3.3：遍历所有相交四边形网格，分别对每个相交四边形网格进行处理，生成不规则的多边形；
步骤4.3.4：由生成的不规则的多边形网格代替相应的相交四边形网格，完成对相交四边形网格的处理，实现露天矿采场三角形网格对所有煤层顶底板四边形网格的裁剪，从而实现对露天矿矿床地质模型的更新。
6.根据权利要求4所述的一种露天矿采场三角形网格裁剪煤层四边形网格的方法，其特征在于：所述步骤4.3.3的具体处理方法为：
第一步，选取位于相交四边形内部的裁剪多边形的线段；
第二步，选取第一步所选取四边形的位于裁剪多边形外部的顶点，并与第一步得到的线段组成逆时针方向的多边形，由这些多边形组成裁剪后的不规则网格。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711120978.1A CN107886575A (zh) | 2017-11-14 | 2017-11-14 | 一种露天矿采场三角形网格裁剪煤层四边形网格的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711120978.1A CN107886575A (zh) | 2017-11-14 | 2017-11-14 | 一种露天矿采场三角形网格裁剪煤层四边形网格的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107886575A true CN107886575A (zh) | 2018-04-06 |
Family
ID=61777079
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711120978.1A Pending CN107886575A (zh) | 2017-11-14 | 2017-11-14 | 一种露天矿采场三角形网格裁剪煤层四边形网格的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107886575A (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108564635A (zh) * | 2018-04-14 | 2018-09-21 | 深圳市图智能科技有限公司 | 一种基于面模型的断面网格生成方法 |
| CN109754469A (zh) * | 2019-01-10 | 2019-05-14 | 上海联影医疗科技有限公司 | 表面网格切割和切口边界优化处理的方法和系统 |
| CN110032771A (zh) * | 2019-03-18 | 2019-07-19 | 辽宁工程技术大学 | 一种顾及局部细节特征的露天煤矿dem精确裁剪方法 |
| CN110516930A (zh) * | 2019-08-12 | 2019-11-29 | 山东东山新驿煤矿有限公司 | 一种基于Voronoi图的煤层稳定性定量评价方法 |
| CN110673227A (zh) * | 2019-10-31 | 2020-01-10 | 中国石油集团东方地球物理勘探有限责任公司 | 地层不整合交切的处理方法及处理装置 |
| CN111354086A (zh) * | 2018-12-24 | 2020-06-30 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种适用于dsmc方法网格位置属性判断的双向三维扫描方法 |
| CN111723503A (zh) * | 2020-06-04 | 2020-09-29 | 中国飞机强度研究所 | 一种曲面石墨加热器设计方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001076167A (ja) * | 1999-08-31 | 2001-03-23 | Square Co Ltd | コンピュータ読み取り可能な記録媒体、画像処理方法および画像処理装置 |
| CN101515372A (zh) * | 2009-02-04 | 2009-08-26 | 北京石油化工学院 | 基于虚拟地质模型的可视化分析预测方法 |
| CN102360507A (zh) * | 2011-10-19 | 2012-02-22 | 浙江大学 | 一种基于交点排序的多边形裁剪的方法 |
| CN102496179A (zh) * | 2011-11-08 | 2012-06-13 | 中国测绘科学研究院 | 一种三维场景的裁剪方法 |
-
2017
- 2017-11-14 CN CN201711120978.1A patent/CN107886575A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001076167A (ja) * | 1999-08-31 | 2001-03-23 | Square Co Ltd | コンピュータ読み取り可能な記録媒体、画像処理方法および画像処理装置 |
| CN101515372A (zh) * | 2009-02-04 | 2009-08-26 | 北京石油化工学院 | 基于虚拟地质模型的可视化分析预测方法 |
| CN102360507A (zh) * | 2011-10-19 | 2012-02-22 | 浙江大学 | 一种基于交点排序的多边形裁剪的方法 |
| CN102496179A (zh) * | 2011-11-08 | 2012-06-13 | 中国测绘科学研究院 | 一种三维场景的裁剪方法 |
Non-Patent Citations (4)
| Title |
|---|
| 刘光伟: "数字化露天矿虚拟开采模型构建技术及应用研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
| 王志宏 等: "一种基于地质界面的三维模型裁剪法", 《辽宁工程技术大学学报》 * |
| 陈应显 等: "基于约束剖分三维露天矿采场模型的建立", 《煤炭学报》 * |
| 陈应显: "露天矿矿床三维建模技术及可视化研究", 《中国优秀博硕士学位论文全文数据库 (硕士) 工程科技Ⅰ辑》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108564635A (zh) * | 2018-04-14 | 2018-09-21 | 深圳市图智能科技有限公司 | 一种基于面模型的断面网格生成方法 |
| CN111354086B (zh) * | 2018-12-24 | 2023-04-14 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种适用于dsmc方法网格位置属性判断的双向三维扫描方法 |
| CN111354086A (zh) * | 2018-12-24 | 2020-06-30 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种适用于dsmc方法网格位置属性判断的双向三维扫描方法 |
| CN109754469A (zh) * | 2019-01-10 | 2019-05-14 | 上海联影医疗科技有限公司 | 表面网格切割和切口边界优化处理的方法和系统 |
| CN109754469B (zh) * | 2019-01-10 | 2023-08-22 | 上海联影医疗科技股份有限公司 | 表面网格切割和切口边界优化处理的方法和系统 |
| CN110032771A (zh) * | 2019-03-18 | 2019-07-19 | 辽宁工程技术大学 | 一种顾及局部细节特征的露天煤矿dem精确裁剪方法 |
| CN110032771B (zh) * | 2019-03-18 | 2023-04-07 | 辽宁工程技术大学 | 一种顾及局部细节特征的露天煤矿dem精确裁剪方法 |
| CN110516930B (zh) * | 2019-08-12 | 2022-11-04 | 山东东山新驿煤矿有限公司 | 一种基于Voronoi图的煤层稳定性定量评价方法 |
| CN110516930A (zh) * | 2019-08-12 | 2019-11-29 | 山东东山新驿煤矿有限公司 | 一种基于Voronoi图的煤层稳定性定量评价方法 |
| CN110673227B (zh) * | 2019-10-31 | 2021-07-09 | 中国石油集团东方地球物理勘探有限责任公司 | 地层不整合交切的处理方法及处理装置 |
| CN110673227A (zh) * | 2019-10-31 | 2020-01-10 | 中国石油集团东方地球物理勘探有限责任公司 | 地层不整合交切的处理方法及处理装置 |
| CN111723503B (zh) * | 2020-06-04 | 2021-12-28 | 中国飞机强度研究所 | 一种曲面石墨加热器设计方法 |
| CN111723503A (zh) * | 2020-06-04 | 2020-09-29 | 中国飞机强度研究所 | 一种曲面石墨加热器设计方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107886575A (zh) | 一种露天矿采场三角形网格裁剪煤层四边形网格的方法 | |
| CN100530167C (zh) | 固体矿床三维可视化储量计算系统及计算方法 | |
| CN100478996C (zh) | 水利水电工程地质信息的三维统一模型构建方法 | |
| CN102254349B (zh) | 一种使用钻孔数据构建沉积地层系统三维实体模型的方法 | |
| CN101582173B (zh) | 复杂地质构造块状模型构建方法 | |
| CN106529755A (zh) | 一种矿山地质资源储量管理方法 | |
| CN103699751A (zh) | 一种基于空间矢量的砂体储层构型建模方法及系统 | |
| WO2018222331A1 (en) | Constructing structural models of the subsurface | |
| CN103514631A (zh) | 基于钻孔数据的海量三维地质模型网格式并行构建方法 | |
| CN102609982B (zh) | 空间地质数据非结构化模式的拓扑发现方法 | |
| CN102129712A (zh) | 基于多地层及三维的土石方数量的三角网模型构建方法 | |
| CN106023313A (zh) | 一种三维地质数字勘察报告生成方法 | |
| CN109102564B (zh) | 一种复杂地质体数值模型的耦合建模方法 | |
| CN102222365A (zh) | 复杂空间的曲面重构方法 | |
| CN104318618B (zh) | 一种广义三棱柱空间数据模型的三维剖切方法 | |
| Guo et al. | Explicit-implicit-integrated 3-D geological modelling approach: A case study of the Xianyan Demolition Volcano (Fujian, China) | |
| CN102194252A (zh) | 一种基于地质层面结构的三角形格架网格生成方法 | |
| CN104635262A (zh) | 一种基于增强型矩形网格的正逆断层等值线自动生成方法 | |
| CN105321204A (zh) | 一种三维地质栅格体模型构建方法 | |
| CN112150582A (zh) | 一种面向多模态数据的地质剖面图近似表达方法 | |
| CN104866682B (zh) | 一种基于地面形变进行页岩气勘探区构造应力场反演的方法 | |
| CN110705168A (zh) | 构造应力场的模拟方法 | |
| Turner et al. | A review of geological modeling | |
| CN105719342A (zh) | 一种地裂缝地质体的三维建模可视化方法及装置 | |
| CN117197377A (zh) | 矿山开采区三维动态可视化建模方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180406 |