【技术实现步骤摘要】
本专利技术涉及地球物理的三维地质建模，尤其涉及一种不连续封闭矿体、岩体的建模方法、装置、设备及介质。
技术介绍
1、在固体矿、油气矿产勘查中，不连续封闭矿体、岩体，体，其建模方法在储量评估以及地震、重磁电正反演上有重要的应用。固体矿多发育于火成岩沉积环境，地质构造多以非层状的多值地质体为主，地质体形状及相互之间的不整合接触关系比较复杂，所获取的数据比如钻孔、重磁电二维剖面数据都非常稀疏，导致复杂岩体或固体矿体的三维地质模型建立比较困难。
2、现有技术中，一般固体矿或者岩体的建模数据主要通过重磁电震反演剖面、钻孔联井数据以及地表地图测斜数据得来，主要包括以下数据：重磁电震剖面解释的地质体边界线；地表露头地图的不整合接触线及测斜数据；大量钻孔地层岩性段信息数据。这些数据在剖面上和钻孔上分布比较密，在线间和井间分布非常稀疏，类似地震的二维测线工区构造解释数据和测井数据；同时由于这些数据表示的是复杂的多值地质体而非层状地层体，而此种地质体的形状比简单的盐丘要复杂，所以需要此种复杂地质体的建模算法可以满足支持其稀疏数据的特点，并且尽量可以引入各种地质约束条件约束模型符合地质预期。
技术实现思路
1、有鉴于此，本专利技术实施例的目的在于提出一种不连续封闭矿体、岩体的建模方法、装置、设备及可读介质，其中，所述方法是基于离散隐函数的不连续封闭矿体、岩体的建模方法，对于不连续封闭矿体、岩体的数据进行预处理，建立其表面的点集以及梯度信息，以此为根据建立断层面约束下的计算网格，在计算网格上建
2、基于上述目的，本专利技术实施例的一方面提供了一种不连续封闭矿体、岩体的建模方法，包括以下步骤：获取不连续封闭矿体、岩体的断层数据和非断层数据并对其进行预处理，得到所述不连续封闭矿体、岩体的表面点集和表面梯度信息；根据所述断层数据建立断面模型，根据所述不连续封闭矿体、岩体的建模需求和所述断面模型设置所述不连续封闭矿体、岩体的建模的计算网格；基于所述不连续封闭矿体、岩体的表面点集和表面梯度信息在所述计算网格上进行插值，生成所述不连续封闭矿体、岩体的符号距离场，并且提取所述符号距离场中的等值面片，得到所述不连续封闭矿体、岩体的模型。
3、在一些实施方式中，所述获取不连续封闭矿体、岩体的断层数据和非断层数据并对其进行预处理，得到所述不连续封闭矿体、岩体的表面点集和表面梯度信息包括：根据获取的不连续封闭矿体、岩体的断层数据和非断层数据建立所述不连续封闭矿体、岩体的二维模型，在所述二维模型上得到断层和非断层对应的折线和不连续多边形；基于最小张力样条函数插值所述折线和不连续多边形，得到所述断层和非断层对应的曲线段和曲线多边形；遍历计算所述不连续封闭矿体、岩体的所有曲线多边形上点的梯度信息，得到所述不连续封闭矿体、岩体的表面梯度信息。
4、在一些实施方式中，所述获取不连续封闭矿体、岩体的断层数据和非断层数据并对其进行预处理，得到所述不连续封闭矿体、岩体的表面点集和表面梯度信息还包括：响应于在同一断层的不同剖面方向存在相交的断层折线，获取相交的交点并且将所述交点插入到相交的断层折线上，以降低断层折线相交产生的闭合差；响应于在同一非断层的不同剖面方向存在相交的多边形，获取相交的交点并且将所述交点插入到相交的多边形上，以降低多边形相交产生的闭合差。
5、在一些实施方式中，所述遍历计算所述不连续封闭矿体、岩体的所有曲线多边形上点的梯度信息，得到所述不连续封闭矿体、岩体的表面梯度信息包括：计算当前曲线多边形所在的平面法矢并且归一化；依次计算所述当前曲线多边形上的点在所述当前曲线多边形上的切向量并且归一化；按序遍历所述当前曲线多边形上的点的切向量并且分别与所述平面法矢叉积运算，得到所述当前曲线多边形的梯度信息；按照以上步骤遍历所述不连续封闭矿体、岩体的所有曲线多边形对应的梯度信息，得到所述不连续封闭矿体、岩体的表面梯度信息。
6、在一些实施方式中，所述根据所述断层数据建立断面模型，根据所述不连续封闭矿体、岩体的建模需求和所述断面模型设置所述不连续封闭矿体、岩体的建模的计算网格包括：根据所述断层数据计算断面的点集的凸包边界多边形；根据所述断面的加密折线段将所述凸包边界多边形切割为块区域的集合；对所述块区域的集合中的所有块区域进行三角剖分，得到所述断面的三角网格面模型；根据所述不连续封闭矿体、岩体的建模需求和所述三角网格面模型设置所述不连续封闭矿体、岩体的建模的计算网格。
7、在一些实施方式中，所述根据所述不连续封闭矿体、岩体的建模需求和所述三角网格面模型设置所述不连续封闭矿体、岩体的建模的计算网格包括：按照所述不连续封闭矿体、岩体的建模精度和复杂程度设置建模的网格尺度，并且基于所述网格尺度得到所述不连续封闭矿体、岩体的建模的初步计算网格；将所述三角网格面模型引入所述初步计算网格中，对所述初步计算网格进行约束以得到所述不连续封闭矿体、岩体的建模的计算网格。
8、在一些实施方式中，所述将所述三角网格面模型引入所述初步计算网格中，对所述初步计算网格进行约束以得到所述不连续封闭矿体、岩体的建模的计算网格包括：计算所述三角网格面模型的三角网格与所述初步计算网格的交点，并且根据所述交点判断相邻初步计算网格是否连通，通过对所述初步计算网格是否连通进行标记，得到所述不连续封闭矿体、岩体的建模的计算网格。
9、在一些实施方式中，所述基于所述不连续封闭矿体、岩体的表面点集和表面梯度信息在所述计算网格上进行插值，生成所述不连续封闭矿体、岩体的符号距离场，并且提取所述符号距离场中的等值面片，得到所述不连续封闭矿体、岩体的模型包括：根据所述不连续封闭矿体、岩体的表面点集和表面梯度信息在所述计算网格上插值建立离散最小曲率插值方程组，对所述离散最小曲率差值方程组求解得到所述不连续封闭矿体、岩体的符号距离场；遍历所有计算网格，提取所述符号距离场中的等值面片，得到所述不连续封闭矿体、岩体的模型。
10、在一些实施方式中，所述根据所述不连续封闭矿体、岩体的表面点集和表面梯度信息在所述计算网格上插值建立离散最小曲率插值方程组，对所述离散最小曲率差值方程组求解得到所述不连续封闭矿体、岩体的符号距离场包括：将所述表面点集的每个点按照所在的计算网格进行分配，得到所述每个点对应的计算网格；根据所述不连续封闭矿体、岩体的线性约束信息在所述每个点对应的计算网格中建立约束方程组；通过集合所述约束方程组和引入所述不连续封闭矿体、岩体的不连续性约束信息建立离散最小曲率插值方程组；求解所述离散最小曲率插值方程组，得到所述不连续封闭矿体、岩体的符号距离场。
11、在一些实施方式中，所述根据所述不连续封闭矿体、岩体的线性约束信息在所述每个点对应的计算网格中建立约束方程组包括：建立所述每个点在对应的计算网格中的线性约束局部方程组；建立所述每个点的梯度信息在对应的计算网格中的线性约束局部方程组；基于能量最小几何偏微分方程建立所述计本文档来自技高网...
【技术保护点】
1.一种不连续封闭矿体、岩体的建模方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的方法，其特征在于，所述获取不连续封闭矿体、岩体的断层数据和非断层数据并对其进行预处理，得到所述不连续封闭矿体、岩体的表面点集和表面梯度信息包括：
3.根据权利要求2所述的方法，其特征在于，所述获取不连续封闭矿体、岩体的断层数据和非断层数据并对其进行预处理，得到所述不连续封闭矿体、岩体的表面点集和表面梯度信息还包括：
4.根据权利要求2所述的方法，其特征在于，所述遍历计算所述不连续封闭矿体、岩体的所有曲线多边形上点的梯度信息，得到所述不连续封闭矿体、岩体的表面梯度信息包括：
5.根据权利要求1所述的方法，其特征在于，所述根据所述断层数据建立断面模型，根据所述不连续封闭矿体、岩体的建模需求和所述断面模型设置所述不连续封闭矿体、岩体的建模的计算网格包括：
6.根据权利要求5所述的方法，其特征在于，所述根据所述不连续封闭矿体、岩体的建模需求和所述三角网格面模型设置所述不连续封闭矿体、岩体的建模的计算网格包括：
7.根据权利要求
8.根据权利要求1所述的方法，其特征在于，所述基于所述不连续封闭矿体、岩体的表面点集和表面梯度信息在所述计算网格上进行插值，生成所述不连续封闭矿体、岩体的符号距离场，并且提取所述符号距离场中的等值面片，得到所述不连续封闭矿体、岩体的模型包括：
9.根据权利要求8所述的方法，其特征在于，所述根据所述不连续封闭矿体、岩体的表面点集和表面梯度信息在所述计算网格上插值建立离散最小曲率插值方程组，对所述离散最小曲率差值方程组求解得到所述不连续封闭矿体、岩体的符号距离场包括：
10.根据权利要求9所述的方法，其特征在于，所述根据所述不连续封闭矿体、岩体的线性约束信息在所述每个点对应的计算网格中建立约束方程组包括：
11.根据权利要求9所述的方法，其特征在于，所述求解所述离散最小曲率插值方程组，得到所述不连续封闭矿体、岩体的符号距离场包括：
12.根据权利要求8所述的方法，其特征在于，所述遍历所述计算网格，提取所述符号距离场中的等值面片，得到所述不连续封闭矿体、岩体的模型包括：
13.根据权利要求12所述的方法，其特征在于，所述遍历所有计算网格，并且根据所述计算网格是否被所述断面模型切割匹配对应的提取算法提取所述符号距离场中的等值面片，以得到所述不连续封闭矿体、岩体的模型包括：
14.一种不连续封闭矿体、岩体的建模装置，其特征在于，包括：
15.一种计算机设备，其特征在于，包括：
16.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1-13任意一项所述方法的步骤。
...【技术特征摘要】
1.一种不连续封闭矿体、岩体的建模方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的方法，其特征在于，所述获取不连续封闭矿体、岩体的断层数据和非断层数据并对其进行预处理，得到所述不连续封闭矿体、岩体的表面点集和表面梯度信息包括：
3.根据权利要求2所述的方法，其特征在于，所述获取不连续封闭矿体、岩体的断层数据和非断层数据并对其进行预处理，得到所述不连续封闭矿体、岩体的表面点集和表面梯度信息还包括：
4.根据权利要求2所述的方法，其特征在于，所述遍历计算所述不连续封闭矿体、岩体的所有曲线多边形上点的梯度信息，得到所述不连续封闭矿体、岩体的表面梯度信息包括：
5.根据权利要求1所述的方法，其特征在于，所述根据所述断层数据建立断面模型，根据所述不连续封闭矿体、岩体的建模需求和所述断面模型设置所述不连续封闭矿体、岩体的建模的计算网格包括：
6.根据权利要求5所述的方法，其特征在于，所述根据所述不连续封闭矿体、岩体的建模需求和所述三角网格面模型设置所述不连续封闭矿体、岩体的建模的计算网格包括：
7.根据权利要求6所述的方法，其特征在于，所述将所述三角网格面模型引入所述初步计算网格中，对所述初步计算网格进行约束以得到所述不连续封闭矿体、岩体的建模的计算网格包括：
8.根据权利要求1所述的方法，其特征在于，所述基于所述不连续封闭矿体、岩体的表面点集和表面梯度信息在所述计算网格上进行插值，生成所述不连续...
【专利技术属性】
技术研发人员：于海生，崔京彬，陈茂山，王文涛，郝彦国，李咏梅，
申请(专利权)人：中国石油天然气集团有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。