一种岩石物性结构三维建模的方法技术

本发明专利技术公开了一种岩石物性结构三维建模的方法，其包括取地球物理测井数据和钻孔空间信息数据，建立地下岩石物性的三维数据场；选择岩石物性参数值或参数值组合，确定物性参数值的值域范围；将三维数据场中的数据映射为网格节点状态值；将相邻八个网格节点状态值组合成六面体，并根据六面体的网格节点状态值选择对应体构型的索引值；计算岩石物性参数值组合等值面与六面体棱边的交点；根据对应的体构型连接棱边交点，创建岩石物性三维结构模型轮廓表面。本发明专利技术能够解决现有技术中岩石物性结构三维建模的准确性和效率较低的问题，准确率高、测定速度快、范围大。

A method of three-dimensional modeling of rock physical structure

【技术实现步骤摘要】
一种岩石物性结构三维建模的方法
本专利技术涉及地质工程领域，具体涉及一种岩石物性结构三维建模的方法。
技术介绍
三维建模一直是工程、地质、矿产行业研究的技术热点和难点问题。在研究地下各种复杂的地质结构时，根据有限的岩石物性参数，构建地下地质体三维空间属性模型是岩石结构分析的一种重要手段，已经有很多的学者和研究文献对其相关技术问题进行了研究和讨论。其中，岩石物性结构三维建模主要是针对各种地质数据、地球物理数据、测井数据及其他方式获取的地下岩石物性参数，建立三维物性空间物性数据场，并从物性数据场构建不同物性分布的三维几何结构模型，为地下空间的三维工程地质分析提供重要的数据基础。现有岩石物性结构三维建模主要包括两个步骤，其一是利用各种方法采集的有限的物性数据，采用空间插值算法，建立三维地下空间的物性数据场。其次，在三维物性数据场的基础上，生成不同岩石物性的三维空间结构模型。现有生成不同岩石物性的三维空间结构模型的方法，主要包括通过移动立方体网格方法MarchingCubes算法，即MC算法，提取三维数据场的等值面。MC算法来源于人体CT医学成像，每次计算获取特定属性值的三维等值面。若要显示完整的岩石物性结构或一定的地质构造，则需要进行多次的等值面提取计算。而每一次的计算都需要遍历整个网格数据，计算效率比较低，完成特定岩石物性结构分析需要大量时间，因此，在实际研究中，常常会为了保证效率减少计算的等值面数量，从而影响建模准确性。另外，对于地质模型来说，按照视觉观察原理，根据观察角度不同，每次观察到的只是地层(构造)的表面，内部或者深部进一步信息不作为视觉信息进行处理，而通过多次等值面提取计算，则需要会将大量地层内部或者深部信息纳入处理中，占用计算设备时间，严重影响建模效率。现有对于MC算法的研究，主要问题集中在等值面连接的二义性的消除、优化等值面存储和计算效率的提高等方面。僧德文等(2006)提出了采用线性八叉树的方式，对有效体元进行编码，从而实现对计算效率进行优化。李程(2014)将六面体分解程四面体的方法来避免二义性的问题。帅仁俊等(2016)为了提高计算速度，提出了基于黄金分割的等值面与棱边的交点算法。王旭初(2012)等提出了构建体元邻接查找子表的方法，使面绘制时只处理有效体元，从而提高计算效率的方法。王铮，李瑞明(2017)提出了基于插值点连线交点的解决方法来解决等值面连接的歧义问题。在地质工程领域中，该方法也被引入到地层建模、矿体建模及地下岩石属性的建模中。杨鸿翼等(2008)将MC算法用于矿化等值面提取及三维模拟；夏艳华等(2012)采用MC算法实现了三维地层建模；赵玲等(2018)采用MC算法用于建立三维数字岩心模型；王铮(2018)采用MC算法，利用地质勘探剖面数据建立矿体的表面模型。综合分析以上对于MC方法的应用上，大多数的文献只是讨论了该方法在特定领域的应用上，对于方法原理、构型设计讨论不多。而且均是讨论针对单一等值面的构建问题，现有技术中缺乏对于利用岩石的物性数据场来构建多个属性等值面的情况的研究。
技术实现思路
本专利技术针对现有技术中的上述不足，提供了一种能够解决现有技术中岩石物性结构三维建模的准确性和效率较低的问题的岩石物性结构三维建模的方法。为解决上述技术问题，本专利技术采用了下列技术方案：提供了一种岩石物性结构三维建模的方法，其包括如下步骤：S1、获取地球物理测井数据和钻孔空间信息数据，建立地下岩石物性的三维数据场；S2、选择岩石物性参数值或参数值组合，确定物性参数值的值域范围；S3、将满足物性参数值域范围的三维数据场中的网格节点值映射为网格节点状态值；S4、将相邻八个网格节点状态值组合成六面体，并根据六面体的网格节点状态值计算索引值并选择对应的体构型；S5、根据所述索引值，建立体构型和构型关系表；S6、根据构型关系表，计算体构型上的物性等值面与六面体棱边的交点坐标；S7、遍历所有网格节点，重复S6，得到物性等值面与所有六面体棱边的交点坐标；S8、根据构型关系表，将所有交点连接成三角网，得到岩石物性结构的三维模型。进一步地，地球物理测井数据包括波速、密度、孔隙度、电阻率；钻孔空间信息数据包括钻孔坐标，孔深、测斜。进一步地，建立地下岩石物性的三维数据场的方法包括基于地层约束的线性插值算法和反距离加权插值算法。进一步地，将满足物性参数值域范围的三维数据场中的网格节点值映射为网格节点状态值的方法包括：S1、统计三维数据场中所有网格节点物性值的分布范围，并将其映射到设定的有限颜色空间中；S2、从物性参数值中选取若干物性值作为比较值；S3、将比较值对应到有限颜色空间中，确定比较值的邻域范围；S4、遍历整个网格节点，将所有节点的物性值分别与比较值进行比较，当节点的物性值位于任一比较值所在有限颜色空间的邻域内时，则将物性值对应的节点状态值设置为1；当节点的物性值位于所有比较值所在有限颜色空间的邻域外时，则将物性值对应的节点状态值设置为0。进一步地，有限颜色空间为包含有限个离散的物性值点的空间和点对应的RGB颜色值。进一步地，索引值为包含六面体的所有网格节点状态值的数组。进一步地，构型关系表包括与索引值对应的体构型类型、每个构型中的三角面连接顺序表；构型关系表与索引值一一对应。进一步地，建立体构型的方法包括如下步骤：s1、对六面体中的相邻两个节点的状态值进行判断，当两个节点状态值不同时，则两个节点对应的棱边存在一个交点；S2、将棱边的交点按构型关系表上的连接顺序进行三角网连接，构成封闭的三角形面；S3、对生成的若干三角形面进行比较，取使整个体构型体积最大化的三角形面，为体构型结构面；S4、将所有体构型结构面连接，得到目标体构型。进一步地，根据构型关系表，确定体构型上的物性等值面与六面体棱边的交点坐标的方法包括如下步骤：S1、建立空间直角坐标系，选择一个网格节点作为坐标原点，按照右手螺旋法则，选择六面体的三个相交棱边为坐标系的三个坐标轴；S2、选择一个等值面，将等值面上与棱边相交的直线投影到坐标系中，以计算等值面对应直线与坐标轴的交点坐标；S3、每个棱边上选取最靠近显示状态为0的节点的交点，作为等值面与六面体的棱边的交点；S4、重复步骤S2和S3，对体构型的每个等值面均进行计算，直到得到体构型所有棱边的交点坐标。进一步地，计算等值面对应直线与坐标轴的交点的坐标的计算公式为：其中，i为网格节点组合成的六面体单元编号，px(i)、py(i)、pz(i)分别为等值面与六面体单元i的三个棱边的交点坐标，n为组合的等值面数量，x、y、z分别为六面体单元的三条棱边组成的三个坐标轴；dx、dy、dz为棱边方向，沿着坐标轴方向，如果棱边节点状态从1到0，则棱边方向为正，否则为负。本专利技术提供的上述岩本文档来自技高网...

【技术保护点】
1.一种岩石物性结构三维建模的方法，其特征在于，包括如下步骤：/nS1、获取地球物理测井数据和钻孔空间信息数据，建立地下岩石物性的三维数据场；/nS2、选择岩石物性参数值或参数值组合，确定物性参数值的值域范围；/nS3、将满足物性参数值域范围的三维数据场中的网格节点值映射为网格节点状态值；/nS4、将相邻八个网格节点状态值组合成六面体，并根据六面体的网格节点状态值计算索引值并选择对应的体构型；/nS5、根据所述索引值，建立体构型和构型关系表；/nS6、根据构型关系表，计算体构型上的物性等值面与六面体棱边的交点坐标；/nS7、遍历所有网格节点，重复S6，得到物性等值面与所有六面体棱边的交点坐标；/nS8、根据构型关系表，将所有交点连接成三角网，得到岩石物性结构的三维模型。/n

【技术特征摘要】
1.一种岩石物性结构三维建模的方法，其特征在于，包括如下步骤：
S1、获取地球物理测井数据和钻孔空间信息数据，建立地下岩石物性的三维数据场；
S2、选择岩石物性参数值或参数值组合，确定物性参数值的值域范围；
S3、将满足物性参数值域范围的三维数据场中的网格节点值映射为网格节点状态值；
S4、将相邻八个网格节点状态值组合成六面体，并根据六面体的网格节点状态值计算索引值并选择对应的体构型；
S5、根据所述索引值，建立体构型和构型关系表；
S6、根据构型关系表，计算体构型上的物性等值面与六面体棱边的交点坐标；
S7、遍历所有网格节点，重复S6，得到物性等值面与所有六面体棱边的交点坐标；
S8、根据构型关系表，将所有交点连接成三角网，得到岩石物性结构的三维模型。


2.根据权利要求1所述的岩石物性结构三维建模的方法，其特征在于，所述地球物理测井数据包括波速、密度、孔隙度、电阻率；钻孔空间信息数据包括钻孔坐标，孔深、测斜。


3.根据权利要求2所述的岩石物性结构三维建模的方法，其特征在于，所述建立地下岩石物性的三维数据场的方法包括基于地层约束的线性插值算法和反距离加权插值算法。


4.根据权利要求1所述的岩石物性结构三维建模的方法，其特征在于，所述将满足物性参数值域范围的三维数据场中的网格节点值映射为网格节点状态值的方法包括：
S1、统计三维数据场中所有网格节点物性值的分布范围，并将其映射到设定的有限颜色空间中；
S2、从物性参数值中选取若干物性值作为比较值；
S3、将比较值对应到有限颜色空间中，确定比较值的邻域范围；
S4、遍历整个网格节点，将所有节点的物性值分别与比较值进行比较，当节点的物性值位于任一比较值所在有限颜色空间的邻域内时，则将物性值对应的节点状态值设置为1；
当节点的物性值位于所有比较值所在有限颜色空间的邻域外时，则将物性值对应的节点状态值设置为0。


5.根据权利要求4所述的岩石物性结构三维建模的方法，其特征在于，所述有限颜色空间为包含有限个离散的物性值点的空间和点对应的RGB颜色值。


...

【专利技术属性】
技术研发人员：简兴祥，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：四川;51