一种基于地统计学和符号矩阵的三维地质建模方法技术

本发明专利技术公开了一种基于地统计学和符号矩阵的三维地质建模的方法，本方法包括以下步骤：预建模处理：确定建模区域的边界范围，初始化三维地质模型；数据采集：获得所述边界范围内各钻孔数据的对应的煤层数据和符号距离函数；计算边界范围区内的煤厚和距离函数估值；模型建立：对栅格单元的点进行赋值，对模型进行可视化分析。本发明专利技术基于地统计学和符号矩阵的三维地质建模的方法适用于尖灭、分叉煤层地质区域建模，其自动化程度高，可应用于煤层开采设计。采设计。采设计。

【技术实现步骤摘要】
一种基于地统计学和符号矩阵的三维地质建模方法


[0001]本专利技术属于地质建模领域，特别涉及基于地统计学和符号矩阵的三维地质建模的方法。

技术介绍

[0002]近些年来，继数字地球和智慧城市之后，数字孪生的提出进一步的推动了各个行业的数字化、信息化和智能化。在地质行业中，也在数字孪生之后，提出了数字孪生矿山和石油地质孪生等行业或国家级研究计划，而三维地质建模技术是数字孪生矿山和石油数字孪生的基础和关键技术，也是数字地球和智慧城市的重要组成部分。与传统的二维地质数据相比，三维地质模型能够直观地展示地层的地质结构、空间拓扑、岩性等信息，不仅为研究人员分析地质构造、断层分布等工作提供准确的信息，也能为地下空间的开发利用提供可靠依据。三维地质建模的完成一般依赖于钻孔数据，而钻孔数据往往都是一系列离散的、空间上分布不均匀的有限的反映地理现象特征的勘查数据，因此要完成地层的建模，离不开对这些离散点进行合理的插值。其中，常用的空间插值方法有反距离加权、克里金、趋势面、样条函数、自然邻域法，而IDW以其插值原理简单、计算简便且符合地理学第一定律而被广泛使用。而各插值方法的前提就是插值点的确立，因此，如何将不均匀的、离散的钻孔数据，通过虚拟钻孔点的插入，形成均匀的点集，是执行插值的首要任务。同时由于成煤时期的泥炭沼泽基底不平、地壳运动、地理环境、植物堆积及沉积后受到地质构造破坏等情况使得煤层出现分叉、尖灭现象，煤矿多数呈现为多个薄煤层。准确详实的三维地质模型可以将地下煤层的分布更直观的展示出来，为煤矿的开采设计提供可靠的依据。/>[0003]传统对煤层进行地质建模的方法有也有许多，但它们的实用性都有一定的局限性，如实体模型和三维数字化，需要大量人工解释，只适用于小范围的建模；基于网格化的指示模拟和截断高斯模拟，所建模型连通性较差；基于层面设计的三维地质模型，搭建模型时需要界面切割，所建模型效果不好，自动化程度低；基于有限的地质勘测数据，准确智能的建立三维可视化精细地质模型，也依旧是地质建模领域的一个不小的挑战。

技术实现思路

[0004]本专利技术的主要目的在于提供一种针对尖灭、分叉煤层的通过地统计学和符号矩阵的三维地质建模的方法，可解决上述技术问题。
[0005]为实现上述目的，本专利技术采取的技术方案为：
[0006]本专利技术实施例提供一种针对尖灭、分叉煤层的通过地统计学和符号矩阵的三维地质建模的方法，包括以下步骤：
[0007]S1、预建模处理：确定建模区域的边界范围，初始化三维地质模型；
[0008]S2、数据采集：获得所述边界范围内各钻孔数据的对应的煤层数据和符号距离函数；
[0009]S3、计算边界范围区内的煤厚和距离函数估值；
[0010]S4、模型建立：对栅格单元的点进行赋值，对模型进行可视化分析。
[0011]进一步地，所述步骤S1包括：
[0012]确定建模区域的边界范围，利用GOCAD软件对边界范围进行角点网格预建模处理，初始化三维地质模型。
[0013]进一步地，所述步骤S2包括：
[0014]S21、获得所述目标范围内各钻孔数据的对应的煤层数据，包括煤层顶底板位置及煤厚；
[0015]S22、获得所述目标范围内各钻孔数据对应的符号距离函数。
[0016]进一步地，所述步骤S3包括：
[0017]S31、通过步骤S21的总煤厚，利用全局克里金插值法求待插点煤厚及权重值；
[0018]S32、通过步骤S22的符号距离函数值及S31的权重值求垂直方向下各栅格点的符号距离函数估值。
[0019]进一步地，所述步骤S4包括：
[0020]S41、将计算出的二维坐标下的煤厚与三维坐标下的符号距离函数估值进行数据融合处理，并将数据融合后的属性值赋值到符号矩阵中。
[0021]S42、对模型进行可视化分析，利用符号矩阵的值进行空间建模。
[0022]与现有技术相比，本专利技术具有如下有益效果：
[0023]本专利技术提供的一种针对尖灭、分叉煤层的通过地统计学和符号矩阵的三维地质建模的方法，通过结合显性建模和隐形建模的优点，在煤层构建过程中，针对尖灭、分叉情况，对地质体进行一体化建模，保证了煤层的连续性。且由于网格剖分的方向与拟合方向垂直，赋予属性时也会依据构造的煤层情况进行赋予，不易发生属性混乱。
[0024]在本方法实施前期的数据准备上，只需要将所需拟合的煤层数据和符号距离函数数据输入依据本方法，可直接生成模型中煤层的位置关系，以达到所需要的工程要求。
附图说明
[0025]图1为本专利技术实施例提供的基于地统计学和符号矩阵的三维地质建模的方法流程图；
[0026]图2为本专利技术实施例提供的专利技术步骤功能结构示意图；
[0027]图3为本专利技术实施例提供的符号距离函数计算结构示意图；
[0028]图4为本专利技术实施例提供的拟合后所得结果示意图；
具体实施方式
[0029]为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施方式，对本专利技术进行进一步详细说明，应当理解，此处所描述的具体实施方式仅仅用以解释本专利技术，并不用于限定本专利技术。
[0030]本专利技术提供了基于地统计学和符号矩阵的三维地质建模的方法，参阅图1至图4所示，所述方法包括以下步骤：
[0031]S1、预建模处理：确定建模区域的边界范围，初始化三维地质模型；
[0032]S2、数据采集：获得所述边界范围内各钻孔数据的对应的煤层数据和符号距离函
数；
[0033]S3、计算边界范围区内的煤厚和距离函数估值；
[0034]S4、模型建立：对栅格单元的点进行赋值，对模型进行可视化分析。
[0035]下面分别对上述各个步骤进行详细的说明：
[0036]步骤S1中，确定建模区域的边界范围，利用GOCAD软件对边界范围进行角点网格预建模处理，初始化三维地质模型，可获得模型中各栅格的所有交点坐标矩阵，该模型中含I
×
J
×
K个栅格点；将预建立的网格数据的所有交点坐标储存于G矩阵中，G为一个(xm
×
ym
×
zm
×
S)
×
4的二维矩阵，横坐标共有xm
×
ym
×
zm
×
S个，其中xm为整个角点网格剖分的纵向栅格I的个数，它的序号为1～xm，ym为横向栅格J的个数，序号为1～ym，zm为垂直方向下栅格K的个数，序号为1～zm，S为各交点坐标的属性值，初始设置为0。
[0037]步骤S2包括：
[0038]步骤S21中，获得所述目标范围内各钻孔数据的对应的煤层数据；确定钻孔数据的各坐标(X
i
，Y
i
)，及煤层的顶底板位置(d
顶j
，d
底j
)，将各钻孔数据包括X，Y，M储存至文本文档meih本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于地统计学和符号矩阵的三维地质建模的方法，其特征在于：包括以下步骤：S1、预建模处理：确定建模区域的边界范围，初始化三维地质模型；S2、数据采集：获得所述边界范围内各钻孔数据的对应的煤层数据和符号距离函数；S3、计算边界范围区内的煤厚和距离函数估值；S4、模型建立：对栅格单元的点进行赋值，对模型进行可视化分析。2.根据权利要求1所述的一种基于地统计学和符号矩阵的三维地质建模的方法，其特征在于：所述步骤S1包括：确定建模区域的边界范围，利用GOCAD软件对边界范围进行角点网格预建模处理，初始化三维地质模型。3.根据权利要求1所述的一种基于地统计学和符号矩阵的三维地质建模的方法，其特征在于：所述步骤S2包括：S21、获得所述目标范围内各钻孔数据的对...

【专利技术属性】
技术研发人员：王艳美，柳博文，刘星，郭杨，许沐伟，
申请(专利权)人：安徽理工大学，
类型：发明
国别省市：