CN109003330A - 一种基于基岩边界约束的三维地层建模方法 - Google Patents

Abstract

本发明公开了一种基于基岩边界约束的三维地层建模方法，主要处理步骤包括：(1)处理原始钻孔数据；(2)读取工作区域边界点数据，并根据点的坐标构建边界虚拟钻孔；(3)离散工作区域范围；(4)构建初始三角网；(5)地层尖灭处理；(6)符合沉积规律的地层边界尖灭处理；(7)构建工程地质体GTP模型，本发明方法构建的模型精度较高，更加符合真实地层情况。

Description

一种基于基岩边界约束的三维地层建模方法

技术领域

本发明地理信息技术应用领域，尤其涉及一种基于基岩边界约束的三维地层建模方法。

背景技术

三维地层建模技术在城市工程地质建设中有着极其重要的作用，近年来，城市工程地质体三维地层建模技术得到长足发展，研究出多种三维数据模型和适用于不同工程地质环境下的建模方法，并且取得了良好的效果。目前，三维地层建模的数据源主要是工程地质钻孔，通常是先根据工程建设需要人为设定建模工作区范围，然后基于钻孔进行三维建模。这类方法建模效率较高，但其边界通常是工作区边界，而不是工程地层的自然边界，当设定的工作区内有基岩出露时，构建出来的地层三维模型与真实地层相比就会存在一定差异，同时，也难于实现与基岩地质体的有效融合。

因此，当工作区边界包含基岩埋深零值等值线时，各地层范围受基岩顶面的限制，从上到下各地层范围不断缩小，需要考虑地层范围问题，即如何确定每个地层的范围，并根据地层范围去构建符合真实地层形态的地层。

发明内容

发明目的：本发明针对现有技术存在的问题，提供一种基于基岩边界约束的三维地层建模方法。

技术方案：本发明所述的基于基岩边界约束的三维地层建模方法包括：

(1)将建模区域内的原始钻孔数据读入原始钻孔数据集合中，并按照标准地层顺序表，对缺失地层的原始钻孔添加缺失地层；

(2)读取建模区域边界点数据，并根据边界点的坐标构建边界虚拟钻孔，并将边界虚拟钻孔添加到原始钻孔数据集合中，形成第一钻孔数据集合；

(3)对建模区域范围进行格网离散，并根据建模区域离散点的坐标构建离散点虚拟钻孔，并将离散点虚拟钻孔添加到第一钻孔数据集合中，形成第二钻孔数据集合；

(4)以建模区域范围为约束条件，以原始钻孔、离散点虚拟钻孔和边界虚拟钻孔的坐标为基准点，采用标准的D-TIN构建算法，建立初始三角网；

(5)根据第二钻孔数据集合，采用地层尖灭处理规则对初始三角网进行地层尖灭处理，生成地层尖灭处理三角网，地层尖灭处理后的钻孔数据形成第三钻孔数据集合；

(6)根据第三钻孔数据集合，采用地层边界尖灭处理规则，对地层尖灭处理三角网进行处理，得到地层边界尖灭处理三角网集合，地层边界尖灭处理后的钻孔数据形成第四钻孔数据集合；

(7)根据地层三角网集合，构建工程地质体GTP模型。

进一步的，步骤(1)具体包括：

(1-1)读取建模区域内的原始钻孔数据和对应的钻孔地层数据，并将数据存入原始钻孔数据集合OD＝{od_i|i＝1,2,...,m}，其中，od_i表示第i个原始钻孔，od_i＝{lc_i,stglist_i}，lc_i为ad_i信息，stglist_i为ad_i对应的钻孔地层列表，存储有每个钻孔地层的序号、层顶深度、层底深度和地层深度，m为原始钻孔数量；

(1-2)根据标准地层顺序表，得到标准地层数stmcnt，对AD中钻孔地层总数小于标准地层数stmcnt钻孔进行如下处理：对于该钻孔的钻孔地层中的缺失地层，构建虚拟地层，并设置其序号为缺失地层的序号，再将缺失地层的层顶深度和层底深度均设置为其上邻接地层的层底深度，之后，将虚拟地层信息添加到该钻孔的钻孔地层列表中。

进一步的，步骤(2)具体包括：

(2-1)读取建模区域边界点数据，并将边界点数据存入边界点数据集合BP；

(2-2)遍历边界点数据集合BP中所有的边界点，并判断点的类型，若为基岩边界点，则根据点的坐标，构建基岩边界虚拟钻孔；若为普通边界点，则根据点的坐标，构建普通边界虚拟钻孔；之后，将边界虚拟钻孔存入边界虚拟钻孔集合BD；

(2-3)计算边界虚拟钻孔集合BD中所有钻孔的地层信息，具体处理步骤为：利用原始钻孔数据集合OD，采用IDW空间插值算法，拟合虚拟钻孔地层的层顶深度和层底深度，再计算虚拟钻孔高程和层底深度的差值，得到地层深度，之后，将地层信息添加到该钻孔的钻孔地层列表中；

(2-4)将边界虚拟钻孔集合BD中所有的边界虚拟钻孔，添加到原始钻孔数据集合OD中，形成第一钻孔数据集合OD₁。

进一步的，步骤(3)具体包括：

(3-1)读取建模区域范围数据roi；

(3-2)根据设置的离散间距d对roi进行格网离散处理，并将得到的离散点存入离散点集合DP；

(3-3)遍历离散点集合DP中所有的离散点，并根据点的坐标，构建离散点虚拟钻孔，并存入离散点虚拟钻孔集合DD；

(3-4)计算离散点虚拟钻孔集合DD中所有钻孔的地层信息，具体处理步骤为：利用原始钻孔数据集合OD，采用IDW空间插值算法，拟合虚拟钻孔地层的层顶深度和层底深度，再计算虚拟钻孔高程和层底深度的差值，得到地层深度，之后，将地层信息添加到该钻孔的钻孔地层列表中；

(3-5)将离散点虚拟钻孔集合DD中所有离散点虚拟钻孔，添加到第一钻孔数据集合OD₁中，形成第二钻孔数据集合OD₂。

进一步的，步骤(5)具体包括：

(5-1)遍历初始三角网T₀中所有的三角形的边，根据地层尖灭处理规则，在尖灭处添加尖灭点，得到地层尖灭处理三角网集合T_a，并根据尖灭点的坐标构建地层尖灭处理虚拟钻孔；

(5-2)从第二钻孔数据集合OD₂中获取初始三角网T₀中三角形边两端钻孔的地层信息，用于计算地层尖灭处理虚拟钻孔的地层信息，并将尖灭处理虚拟钻孔添加到第二钻孔数据集合OD₂中，形成第三钻孔数据集合OD₃。

进一步的，步骤(6)具体包括：

(6-1)根据标准地层数stmcnt，构建地层三角网集合TL＝{tl_i|i＝1,…,stmcnt}，并将tl_i的初始值都设为地层尖灭处理三角网T_a，其中，tl_i表示第i个地层三角网；

(6-2)将第三钻孔数据集合OD₃中基岩边界虚拟钻孔中基岩地层的层顶深度和层底深度，设置为该基岩边界虚拟钻孔的高程值；

(6-3)地层边界尖灭处理，具体处理步骤为：

A、按地层顺序，遍历地层三角网集合TL中除基岩地层三角网外的其余地层三角网，并从地层三角网中获取每一个三角形的端点坐标；

B、根据端点坐标，从第三钻孔数据集合OD₃中获得每一个三角形的钻孔数据；

C、统计每个三角形中各钻孔的地层层底深度小于基岩顶面深度的地层个数，记为Cnt；

D、针对Cnt>0的三角形进行如下处理：求取该三角形与基岩顶面的交点，根据交点的坐标，构建地层边界尖灭处理虚拟钻孔，之后根据三角形边两端钻孔的地层信息，计算地层边界尖灭处理虚拟钻孔的地层信息，最后将地层边界尖灭处理虚拟钻孔存入第三钻孔数据集合OD₃中，形成第四钻孔数据集合OD₄；

E、根据三角形与基岩顶面的交点以及三角形的顶点构建新的地层三角形，并将构建的地层三角形添加到当前所在地层三角网中，将原来的三角形舍弃，最终，形成地层边界尖灭处理三角网集合TL₁；

(6-4)处理钻孔数据，具体处理步骤如下：

A、遍历第四钻孔数据集合OD₄中所有的钻孔；

B、除基岩地层外，遍历各钻孔中其余地层，判断其层顶深度h_u和层底深度h_d，分别与基岩地层层顶深度h的大小关系，若关系下述三种情形之一，则将地层层底深度h_d的值更新为h_u：

I、h_u>h且h_d>h，

II、h_u＝h且h_d<h

III、h_u<h且h_d<h。

进一步的，步骤(7)具体包括：

(7-1)读取地层边界尖灭处理三角网集合TL₁中的一个地层三角网；

(7-2)遍历该地层三角网中所有三角形，获取三角形的三个端点坐标，并根据坐标，从第四钻孔数据集合OD₄中获得位于三角形端点的三个钻孔；

(7-3)获取三个钻孔的地层数据，并以该地层的层顶深度为GTP上三角形的顶点坐标、以该地层的层底深度为GTP下三角形顶点坐标、以地层岩性编号为GTP属性构建GTP；

(7-4)循环步骤(7-1)至(7-3)，直到处理完集合TL₁中所有的地层三角网。

有益效果：本发明与现有技术相比，其显著优点是：本发明通过构建基岩顶面，求取各地层面与基岩顶面交点，并根据交点构建地层边界虚拟钻孔，然后再进行地层处理等步骤，形成了一种基于基岩边界约束的工程地质体三维建模方法，模型精度较高，更加符合真实地层情况。

附图说明

图1是本发明提供的基于基岩边界约束的三维地层建模方法的流程图；

图2是钻孔与地层数据结构示意图；

图3是原始钻孔空间分布图；

图4是建模边界点空间分布图；

图5是合并后钻孔空间分布图；

图6是初始三角网示意图；

图7是地层边界尖灭处理虚拟钻孔示意图；

图8是地层边界尖灭处理后三角网示意图；

图9是GTP模型示意图。

具体实施方式

本实施例提供了一种基于基岩边界约束的三维地层建模方法，如图1所示，包括：

(1)将建模区域内的原始钻孔数据读入原始钻孔数据集合中，并按照标准地层顺序表，对原始钻孔添加缺失地层。

该步骤具体包括：

(1-1)读取建模区域内的原始钻孔数据和对应的钻孔地层数据，并将数据存入原始钻孔数据集合OD＝{od_i|i＝1,2,...,m}，其中，od_i表示第i个原始钻孔，od_i＝{lc_i,stglist_i}，lc_i为ad_i信息，stglist_i为ad_i对应的钻孔地层列表，存储有每个钻孔地层的序号、层顶深度、层底深度和地层深度，m为原始钻孔数量。其中，钻孔信息和钻孔地层信息的具体数据结构和关系如图2所示。

例如，以含有基岩零值等值线南京某区域为建模工作区域，从数据库中读取工作区域内29个原始钻孔位置数据和59条钻孔地层数据，得到集合OD＝{od_i|i＝1,2,...,29}，原始钻孔空间分布如图3所示。

(1-2)根据标准地层顺序表，得到标准地层数stmcnt，对AD中钻孔地层总数小于标准地层数stmcnt钻孔进行如下处理：对于该钻孔的钻孔地层中的缺失地层，构建虚拟地层，并设置其序号为缺失地层的序号，再将缺失地层的层顶深度和层底深度均设置为其上邻接地层的层底深度，之后，将虚拟地层信息添加到该钻孔的钻孔地层列表中。其中，标准地层顺序如表1所示，

表1

例如，接上例，标准地层数stmcnt＝7，因此，共处理了23个地层缺失的原始钻孔。

(2)读取建模区域放入边界点数据，并根据边界点的坐标构建边界虚拟钻孔，并将边界虚拟钻孔添加到原始钻孔数据集合中，形成第一钻孔数据集合。

该步骤具体包括：

(2-1)读取建模区域边界点数据，并将边界点数据存入边界点数据集合BP；

(2-2)遍历边界点数据集合BP中所有的边界点，并判断点的类型，若为基岩边界点，则根据点的坐标，构建基岩边界虚拟钻孔；若为普通边界点，则根据点的坐标，构建普通边界虚拟钻孔；之后，将边界虚拟钻孔存入边界虚拟钻孔集合BD；

(2-3)计算边界虚拟钻孔集合BD中所有钻孔的地层信息，具体处理步骤为：利用原始钻孔数据集合OD，采用IDW空间插值算法，拟合虚拟钻孔地层的层顶深度和层底深度，再计算虚拟钻孔高程和层底深度的差值，得到地层深度，之后，将地层信息添加到该钻孔的钻孔地层列表中；

(2-4)将边界虚拟钻孔集合BD中所有的边界虚拟钻孔，添加到原始钻孔数据集合OD中，形成第一钻孔数据集合OD₁。

例如，接上例，从数据库中读取了150个建模区域边界点数据，其空间分布如图4所示，得到集合BP＝{bp_i|i＝1,…,150}，再根据边界点构建边界虚拟钻孔，得到集合BD＝{bd_i|i＝1,…,150}，其中，有90个基岩边界虚拟钻孔和60个普通边界虚拟钻孔，将边界虚拟钻孔添加到集合集合OD中，得到OD₁＝{od_i|i＝1,…,179}。

(3)对建模区域范围进行格网离散，并根据建模区域离散点的坐标构建离散点虚拟钻孔，并将离散点虚拟钻孔添加到第一钻孔数据集合中，形成第二钻孔数据集合。

该步骤具体包括：

(3-1)读取建模区域范围数据roi；

(3-2)根据设置的离散间距d对roi进行格网离散处理，并将得到的离散点存入离散点集合DP；

(3-3)遍历离散点集合DP中所有的离散点，并根据点的坐标，构建离散点虚拟钻孔，并存入离散点虚拟钻孔集合DD；

(3-4)计算离散点虚拟钻孔集合DD中所有钻孔的地层信息，具体处理步骤为：利用原始钻孔数据集合OD，采用IDW空间插值算法，拟合虚拟钻孔地层的层顶深度和层底深度，再计算虚拟钻孔高程和层底深度的差值，得到地层深度，之后，将地层信息添加到该钻孔的钻孔地层列表中；

(3-5)将离散点虚拟钻孔集合DD中所有离散点虚拟钻孔，添加到第一钻孔数据集合OD₁中，形成第二钻孔数据集合OD₂。

例如，接上例，设置d＝200，则对roi离散得到110个离散点，形成离散点集合DP＝{dp_i|i＝1,…,110}；对离散点构建离散虚拟钻孔，得到集合DD＝{dd_i|i＝1,…,110}，将集合DD中所有离散点虚拟钻孔，添加到第一钻孔数据集合OD₁中，形成第二钻孔数据集合OD₂＝{od_i|i＝1,…,289}，边界虚拟钻孔和离散虚拟钻孔空间分布如图5所示。

(4)以建模区域范围为约束条件，以原始钻孔、离散点虚拟钻孔和边界虚拟钻孔的坐标为基准点，采用标准的D-TIN构建算法，建立初始三角网。上例中初始三角网如图6所示。

(5)根据第二钻孔数据集合，采用地层尖灭处理规则对初始三角网进行地层尖灭处理，生成地层尖灭处理三角网，地层尖灭处理后的钻孔数据形成第三钻孔数据集合。

该步骤具体包括：

(5-1)遍历初始三角网T₀中所有的三角形的边，根据地层尖灭处理规则，在尖灭处添加尖灭点，得到地层尖灭处理三角网集合T_a，并根据尖灭点的坐标构建地层尖灭处理虚拟钻孔；

(5-2)从第二钻孔数据集合OD₂中获取初始三角网T₀中三角形边两端钻孔的地层信息，用于计算地层尖灭处理虚拟钻孔的地层信息，并将尖灭处理虚拟钻孔添加到第二钻孔数据集合OD₂中，形成第三钻孔数据集合OD₃。

例如，接上例，对初始三角网T₀尖灭处理，构建了419个地层尖灭虚拟钻孔，得到集合OD₃＝{od_i|i＝1,…,997}。

(6)根据第三钻孔数据集合，采用地层边界尖灭处理规则，对地层尖灭处理三角网进行处理，得到地层边界尖灭处理三角网集合，地层边界尖灭处理后的钻孔数据形成第四钻孔数据集合。

该步骤具体包括

(6-1)根据标准地层数stmcnt，构建地层三角网集合TL＝{tl_i|i＝1,…,stmcnt}，并将tl_i的初始值都设为地层尖灭处理三角网T_a，其中，tl_i表示第i个地层三角网。例如，接上例，stmcnt＝10，则TL＝{tl_i|i＝1,…,7}。

(6-2)将第三钻孔数据集合OD₃中基岩边界虚拟钻孔中基岩地层的层顶深度和层底深度，设置为该基岩边界虚拟钻孔的高程值；接上例，共处理了OD₃中90个基岩边界虚拟钻孔。

(6-3)地层边界尖灭处理，具体处理步骤为：

A、按地层顺序，遍历地层三角网集合TL中除基岩地层三角网外的其余地层三角网，并从地层三角网中获取每一个三角形的端点坐标；

B、根据端点坐标，从第三钻孔数据集合OD₃中获得每一个三角形的钻孔数据；

C、统计每个三角形中各钻孔的地层层底深度小于基岩顶面深度的地层个数，记为Cnt；

D、针对Cnt>0的三角形进行如下处理：求取该三角形与基岩顶面的交点，根据交点的坐标，构建地层边界尖灭处理虚拟钻孔，之后根据三角形边两端钻孔的地层信息，计算地层边界尖灭处理虚拟钻孔的地层信息，最后将地层边界尖灭处理虚拟钻孔存入第三钻孔数据集合OD₃中，形成第四钻孔数据集合OD₄；

E、根据三角形与基岩顶面的交点以及三角形的顶点构建新的地层三角形，并将构建的地层三角形添加到当前所在地层三角网中，将原来的三角形舍弃，最终，形成地层边界尖灭处理三角网集合TL₁；

例如，接上例，地层三角网集合TL除基岩地层三角网外有6个地层，共计840个地层三角形，共构建了如图7所示的690个地层边界尖灭处理虚拟钻孔，得到集合OD₄OD₄＝{od_i|i＝1,…,2395}，将地层边界尖灭处理后的地层三角网集合，记为TL₁，其中一个地层三角网如图8所示。

(6-4)处理钻孔数据，具体处理步骤如下：

A、遍历第四钻孔数据集合OD₄中所有的钻孔；

B、除基岩地层外，遍历各钻孔中其余地层，判断其层顶深度h_u和层底深度h_d，分别与基岩地层层顶深度h的大小关系，若关系下述三种情形之一，则将地层层底深度h_d的值更新为h_u：

I、h_u>h且h_d>h，

II、h_u＝h且h_d<h

III、h_u<h且h_d<h。

(7)根据地层三角网集合，构建工程地质体GTP模型。

该步骤具体包括：

(7-1)读取地层边界尖灭处理三角网集合TL₁中的一个地层三角网；

(7-2)遍历该地层三角网中所有三角形，获取三角形的三个端点坐标，并根据坐标，从第四钻孔数据集合OD₄中获得位于三角形端点的三个钻孔；

(7-3)获取三个钻孔的地层数据，并以该地层的层顶深度为GTP上三角形的顶点坐标、以该地层的层底深度为GTP下三角形顶点坐标、以地层岩性编号为GTP属性构建GTP；

(7-4)循环步骤(7-1)至(7-3)，直到处理完集合TL₁中所有的地层三角网。

接上例，最终构建出如图9中图(a)、(c)所示的包含7个地层的GTP模型。

以上所揭露的仅为本发明一种较佳实施例而已，不能以此来限定本发明之权利范围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。

Claims (7)

1.一种基于基岩边界约束的三维地层建模方法，其特征在于该方法包括：

(1)将建模区域内的原始钻孔数据读入原始钻孔数据集合中，并按照标准地层顺序表，对缺失地层的原始钻孔添加缺失地层；

(2)读取建模区域边界点数据，并根据边界点的坐标构建边界虚拟钻孔，并将边界虚拟钻孔添加到原始钻孔数据集合中，形成第一钻孔数据集合；

(3)对建模区域范围进行格网离散，并根据建模区域离散点的坐标构建离散点虚拟钻孔，并将离散点虚拟钻孔添加到第一钻孔数据集合中，形成第二钻孔数据集合；

(4)以建模区域范围为约束条件，以原始钻孔、离散点虚拟钻孔和边界虚拟钻孔的坐标为基准点，采用标准的D-TIN构建算法，建立初始三角网；

(5)根据第二钻孔数据集合，采用地层尖灭处理规则对初始三角网进行地层尖灭处理，生成地层尖灭处理三角网，地层尖灭处理后的钻孔数据形成第三钻孔数据集合；

(6)根据第三钻孔数据集合，采用地层边界尖灭处理规则，对地层尖灭处理三角网进行处理，得到地层边界尖灭处理三角网集合，地层边界尖灭处理后的钻孔数据形成第四钻孔数据集合；

(7)根据地层三角网集合，构建工程地质体GTP模型。

2.根据权利要求1所述的基于基岩边界约束的三维地层建模方法，其特征在于：步骤(1)具体包括：

(1-1)读取建模区域内的原始钻孔数据和对应的钻孔地层数据，并将数据存入原始钻孔数据集合OD＝{od_i|i＝1,2,...,m}，其中，od_i表示第i个原始钻孔，od_i＝{lc_i,stglist_i}，lc_i为ad_i信息，stglist_i为ad_i对应的钻孔地层列表，存储有每个钻孔地层的序号、层顶深度、层底深度和地层深度，m为原始钻孔数量；

(1-2)根据标准地层顺序表，得到标准地层数stmcnt，对AD中钻孔地层总数小于标准地层数stmcnt钻孔进行如下处理：对于该钻孔的钻孔地层中的缺失地层，构建虚拟地层，并设置其序号为缺失地层的序号，再将缺失地层的层顶深度和层底深度均设置为其上邻接地层的层底深度，之后，将虚拟地层信息添加到该钻孔的钻孔地层列表中。

3.根据权利要求1所述的基于基岩边界约束的三维地层建模方法，其特征在于：步骤(2)具体包括：

(2-1)读取建模区域边界点数据，并将边界点数据存入边界点数据集合BP；

(2-2)遍历边界点数据集合BP中所有的边界点，并判断点的类型，若为基岩边界点，则根据点的坐标，构建基岩边界虚拟钻孔；若为普通边界点，则根据点的坐标，构建普通边界虚拟钻孔；之后，将边界虚拟钻孔存入边界虚拟钻孔集合BD；

(2-3)计算边界虚拟钻孔集合BD中所有钻孔的地层信息，具体处理步骤为：利用原始钻孔数据集合OD，采用IDW空间插值算法，拟合虚拟钻孔地层的层顶深度和层底深度，再计算虚拟钻孔高程和层底深度的差值，得到地层深度，之后，将地层信息添加到该钻孔的钻孔地层列表中；

(2-4)将边界虚拟钻孔集合BD中所有的边界虚拟钻孔，添加到原始钻孔数据集合OD中，形成第一钻孔数据集合OD₁。

4.根据权利要求1所述的基于基岩边界约束的三维地层建模方法，其特征在于：步骤(3)具体包括：

(3-1)读取建模区域范围数据roi；

(3-2)根据设置的离散间距d对roi进行格网离散处理，并将得到的离散点存入离散点集合DP；

(3-3)遍历离散点集合DP中所有的离散点，并根据点的坐标，构建离散点虚拟钻孔，并存入离散点虚拟钻孔集合DD；

(3-4)计算离散点虚拟钻孔集合DD中所有钻孔的地层信息，具体处理步骤为：利用原始钻孔数据集合OD，采用IDW空间插值算法，拟合虚拟钻孔地层的层顶深度和层底深度，再计算虚拟钻孔高程和层底深度的差值，得到地层深度，之后，将地层信息添加到该钻孔的钻孔地层列表中；

(3-5)将离散点虚拟钻孔集合DD中所有离散点虚拟钻孔，添加到第一钻孔数据集合OD₁中，形成第二钻孔数据集合OD₂。

5.根据权利要求1所述的基于基岩边界约束的三维地层建模方法，其特征在于：步骤(5)具体包括：

(5-1)遍历初始三角网T₀中所有的三角形的边，根据地层尖灭处理规则，在尖灭处添加尖灭点，得到地层尖灭处理三角网集合T_a，并根据尖灭点的坐标构建地层尖灭处理虚拟钻孔；

(5-2)从第二钻孔数据集合OD₂中获取初始三角网T₀中三角形边两端钻孔的地层信息，用于计算地层尖灭处理虚拟钻孔的地层信息，并将尖灭处理虚拟钻孔添加到第二钻孔数据集合OD₂中，形成第三钻孔数据集合OD₃。

6.根据权利要求1所述的基于基岩边界约束的三维地层建模方法，其特征在于：步骤(6)具体包括：

(6-1)根据标准地层数stmcnt，构建地层三角网集合TL＝{tl_i|i＝1,…,stmcnt}，并将tl_i的初始值都设为地层尖灭处理三角网T_a，其中，tl_i表示第i个地层三角网；

(6-2)将第三钻孔数据集合OD₃中基岩边界虚拟钻孔中基岩地层的层顶深度和层底深度，设置为该基岩边界虚拟钻孔的高程值；

(6-3)地层边界尖灭处理，具体处理步骤为：

A、按地层顺序，遍历地层三角网集合TL中除基岩地层三角网外的其余地层三角网，并从地层三角网中获取每一个三角形的端点坐标；

B、根据端点坐标，从第三钻孔数据集合OD₃中获得每一个三角形的钻孔数据；

C、统计每个三角形中各钻孔的地层层底深度小于基岩顶面深度的地层个数，记为Cnt；

D、针对Cnt>0的三角形进行如下处理：求取该三角形与基岩顶面的交点，根据交点的坐标，构建地层边界尖灭处理虚拟钻孔，之后根据三角形边两端钻孔的地层信息，计算地层边界尖灭处理虚拟钻孔的地层信息，最后将地层边界尖灭处理虚拟钻孔存入第三钻孔数据集合OD₃中，形成第四钻孔数据集合OD₄；

E、根据三角形与基岩顶面的交点以及三角形的顶点构建新的地层三角形，并将构建的地层三角形添加到当前所在地层三角网中，将原来的三角形舍弃，最终，形成地层边界尖灭处理三角网集合TL₁；

(6-4)处理钻孔数据，具体处理步骤如下：

A、遍历第四钻孔数据集合OD₄中所有的钻孔；

B、除基岩地层外，遍历各钻孔中其余地层，判断其层顶深度h_u和层底深度h_d，分别与基岩地层层顶深度h的大小关系，若关系下述三种情形之一，则将地层层底深度h_d的值更新为h_u：

I、h_u>h且h_d>h，

II、h_u＝h且h_d<h，

III、h_u<h且h_d<h。

7.根据权利要求1所述的基于基岩边界约束的三维地层建模方法，其特征在于：步骤(7)具体包括：

(7-1)读取地层边界尖灭处理三角网集合TL₁中的一个地层三角网；

(7-2)遍历该地层三角网中所有三角形，获取三角形的三个端点坐标，并根据坐标，从第四钻孔数据集合OD₄中获得位于三角形端点的三个钻孔；

(7-3)获取三个钻孔的地层数据，并以该地层的层顶深度为GTP上三角形的顶点坐标、以该地层的层底深度为GTP下三角形顶点坐标、以地层岩性编号为GTP属性构建GTP；

(7-4)循环步骤(7-1)至(7-3)，直到处理完集合TL₁中所有的地层三角网。