申请/专利权人：中南大学;伊春鹿鸣矿业有限公司

申请日：2023-09-11

公开（公告）日：2024-04-16

公开（公告）号：CN117131644B

主分类号：G06F30/18

分类号：G06F30/18;G06F30/20;G06F17/18;G06F111/08;G06F111/06

优先权：

专利状态码：有效-授权

法律状态：2024.04.16#授权;2023.12.15#实质审查的生效;2023.11.28#公开

摘要：本公开实施例中提供了一种基于克里金插值及泰森多边形的冻土钻孔布置方法，属于计算技术领域，具体包括：获得每个初始钻孔点对应的泰森多边形的面积和影响范围，获得整体性指标；获取冻土分布重点区域的几何分布，并与预设的敏感区域对应几何分布合并，得到重点区域几何分布；进行蒙特卡洛均匀采样，随后进行Delaunay三角剖分，得到重点区域网格；选取Kriging经验协方差函数计算重点区域网格中的每个网格重元素的精度指标；将整体性指标和精度指标线性加权得到目标函数，并以钻孔点坐标位置为变量，利用遗传算法迭代求解目标函数获得最优解作为钻孔布置方案。通过本公开的方案，提高了钻孔布置的合理性和精准度。

主权项：1.一种基于克里金插值及泰森多边形的冻土钻孔布置方法，其特征在于，包括：步骤1，在目标区域内拟定初始钻孔点，获得每个初始钻孔点对应的泰森多边形的面积和影响范围，获得整体性指标，其中，所述整体性指标的表达式为 其中，EP为钻孔布置为P时的整体性指标，n为钻孔点总个数，Si为序号为i的泰森多边形网格面积大小，为所有泰森多边形面积的均值，即步骤2，根据历史冻土分布数据，获取冻土分布重点区域的几何分布，并与预设的敏感区域对应几何分布合并，得到重点区域几何分布；步骤3，针对重点区域几何分布中的每个重元素进行蒙特卡洛均匀采样，随后进行Delaunay三角剖分，得到重点区域网格；所述步骤3具体包括：步骤3.1，使用随机方法在重点区域几何分布中任一个重元素Wi内生成大量随机点作为初始采样点，对每个初始采样点，计算其与周围点的距离值，并根据距离阈值筛选出符合均匀分布要求的点；步骤3.2，对于筛选出的无效采样点，进行进一步的筛选，通过迭代计算每个点与最近邻点的距离，直到满足均匀分布条件，最终剩余的有效采样点形成均匀分布的点集I'i作为结果，同时将重元素Wi的边缘线按dmin长度进行等分获取点集I”i，合并点集，令最终结果Ii＝I'i+I”i；步骤3.3，针对点集Ii进行Delaunay三角剖分，获取重元素Wi的网格集信息{Mi|mi1,mi2,…mij,…,miN}，其中N为重元素Wi中的三角网格总个数，mij为重元素Wi中的第j个三角网格元素，N为重元素Wi中的网格元素总数，同时获取重元素Wi中网格面积信息{Si|si1,si2,…sij,…,siN}和中心点位置信息{Pi|pi1,pi2,…pij,…,piN}，得到重点区域网格，sij为区域Wi离散化后的网格集Mi中序号为i的网格元素面积，pij为网格集Mi中序号为i的网格元素中心点；步骤4，选取经验协方差函数计算重点区域网格中的每个网格重元素的精度指标；所述步骤4具体包括：步骤4.1，针对工程情况，在有往年插值数据的情况下，将往年插值时选用的变异函数转化为协方差函数，并以该协方差函数的参数作为经验协方差函数的参数，或者，在往年插值数据缺失的情况下，通过收集工程所在区域相关类型的工程数据获取经验协方差函数的基本参数；步骤4.2，获取了经验协方差函数后，根据钻孔布置P，针对重点区域中的任一网格元素，定义其中心点pj的精度值；步骤4.3，根据精度值计算该中心点对应网格元素的精度函数；步骤4.4，根据重元素Wi离散化后的网格集Mi，得到重点区域网格中的每个网格重元素的精度函数；步骤4.5，将重点区域网格中的每个网格重元素的精度函数累和形成重元素Wi的精度函数，获取精度指标集{RP|R1P,R2P,…,RiP,…,Rm+1P}，其中，RiP表示重元素Wi的精度函数，m为离散冻土总数量，l为敏感区域总数量；步骤5，将整体性指标和精度指标线性加权得到目标函数，并以钻孔点坐标位置为变量，利用遗传算法迭代求解目标函数获得最优解作为钻孔布置方案。

全文数据：

权利要求：

百度查询： 中南大学;伊春鹿鸣矿业有限公司 基于克里金插值及泰森多边形的冻土钻孔布置方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD