A Triangulated Irregular Network Constrained Ordinary Kriging Method for Three-Dimensional Modeling of Faulted Geological Surfaces

Abstract

Three-dimensional (3D) modeling of geological surfaces, such as coal seams and strata horizons, from sparsely sampled data collected in the field, is a crucial task in geological modeling. Interpolation is a common approach for this task to construct continuous geological surface models. However, this problem becomes challenging considering the impact of the faults on geological surfaces. Existing methods tend to solve this problem through three steps, including interpolating stratum and fault surface, applying a fault modeling method to modify the geological surface, and optimizing the modified surface to pass sample points fallen into the fault displacement zone. This paper presents a more concise method to generate a faulted geological surface, in which 1) a constrained Delaunay triangulated irregular network (CD-TIN) is constructed to facilitate the neighborhood search process of the ordinary kriging (OK) interpolation, 2) the CD-TIN is also directly constrained by horizon cutoff lines formed from theoretical fault displacement profiles, and 3) subsequently, neighbors of the location to be estimated are selected effectively in the CD-TIN considering the fault topology. The proposed method significantly improves the time efficiency of the OK interpolation by utilizing the CD-TIN and incorporates fault effects directly into the interpolation process by inserting fault horizontal cutoff lines into CD-TIN. Moreover, by integrating the fault effects directly into the interpolation process, the surface modeling process is accomplished in a single stage instead of two separate stages of interpolation first and then modifying the surface in the fault area. By this strategy, the proposed method significantly improves the time efficiency of the OK interpolation algorithm and achieves more accurate modeling of the faulted geological surface. Experiments were designed to compare the performance of our method with several commonly used approaches, and the results indicate that the proposed TINconstrained OK method achieves better accuracy and efficiency in modeling faulted geological surfaces than other methods. This method could also be used in geospatial interpolation studies, such as meteorological data interpolation. INDEX TERMS 3D geological modeling, fault modeling, interpolation, ordinary kriging, neighborhood search.

Three-Dimensional Information Modeling Based on Incomplete Data for Anchor Engineering

“…Geostatistical analysis methods [13], [14], such as the inverse distance weighting and kriging methods, have been widely used in stratum interpolation for model construction. Jia, QR [15] improved the kriging method, which can be applied to generate a faulted geological surface. Topological information from stratum modeling is another common modeling method that can be used to build a solid model.…”

Bandung City Traffic Classification Map With Machine Learning and Ordinary Kriging

“…. , n} is y(xi), the estimated value of y(x0) at the unsampled position x0 is the weighted sum of value X [23]:…”

Enhancing Efficiency in Explicit 3D Geological Modeling and Accuracy of Attribute Interpolation for Urban Underground Spaces: A Case Study in Shanghai

“…Over the past three decades, various modeling theories and techniques have emerged to address the need for 3D models of sedimentary systems. 3D geological modeling methods are generally divided into two categories: explicit and implicit modeling (Jia et al, 2020). Explicit modeling directly defines each object in the model by obtaining the coordinates of key nodes, line segments, and surfaces to establish geological layers (Jessell et al, 2014).…”

“…This method demonstrated superior fitting effectiveness over the general inverse distance weighting and ordinary kriging methods, especially when there was significant spatial variability in the borehole data. Jia et al (2020) introduced an ordinary kriging method with Triangulated Irregular Networks (TIN) constraints to construct fault surfaces in geological settings. This method integrates the influence of faults directly into the interpolation process, completing surface modeling in a single stage instead of the traditional two-stage process of interpolating and modifying the surface in the fault area.…”