针对克里金(Kriging)算法在复杂地质构造应用中的局限性,提出了一种基于边界约束的复杂曲面插值方法。其基本思想是将断层多边形作为层面边界的约束条件,根据种子点与待插值点穿越多边形的关系为依据,判断待插值点与控制点之间的空间拓扑关系,并将满足条件的种子点利用克里金算法进行插值计算。通过实际数据的测试,解决了传统的网格化插值方案层位与断层无法严格相交以及多重逆掩断层构造的层面拟合等难题,为等值线绘制、地质块状模型构建等提供了新的思路。
According to the limitation of the Kriging algorithm,this paper proposed a new complex surface interpolation algo- rithm based on the boundary constraints of space. The basic idea was the fault polygon as the boundary constraints of the hori- zon, and according to the seed points and the relationship between the interpolation points across the polygon to determine the spatial topological relations between the interpolation points and control points, and would meet the conditions of the seed point be interpolated using the Kriging algorithm. Through the test of actual data to solve the problem of the traditional grid interpola- tion scheme layer and the fault can not be strictly intersect, as well as the horizon of the multi-level inverse fault structure. This paper provided a new way of contour mapping, geological block model.