A three-dimensional geological voxel model modeling method based on drilling data

Abstract

The application discloses a three-dimensional geological voxel model modeling method based on drilling data, and comprises the following steps: obtaining geological stratum drilling data of a target area, preprocessing the geological stratum drilling data, and obtaining independent layered geological stratum drilling point sets; for each layered geological stratum drilling point set, generating a continuous stratum surface; based on the upper and lower boundary constraint relationship of the multi-layer continuous stratum surface, constructing a three-dimensional voxel grid, and judging the stratum to which each voxel unit belongs to obtain a geological stratum voxel model; performing spatial interpolation assignment on the multidimensional geological attribute parameters of each voxel; taking voxel three-dimensional integer indexes as unique associated primary keys to establish a one-to-one mapping relationship between the three-dimensional geological voxel model and a structured database. Based on sparse drilling data, the stratum surface is continuously reconstructed, the stratum is voxelized, and multidimensional geological attribute mapping is realized, so that the spatial integrity and attribute expression capability of geological information modeling are improved, and the reliability of subsequent site selection analysis, disaster assessment and engineering decision is improved.

Description

Technical Field

[0001] This invention relates to the field of three-dimensional geological modeling technology, and in particular to a method for modeling three-dimensional geological voxel models based on borehole data. Background Technology

[0002] In engineering fields such as underground space development and geological assessment, the spatial distribution patterns of underground stratigraphic structure and geological properties are the core decision-making basis for project site selection, construction risk assessment, and long-term operation and maintenance safety. Borehole data, as a direct in-situ survey result of underground stratigraphic information, has outstanding advantages such as strong geological authenticity and traceable parameters, and is currently the most core basic data source for 3D geological modeling.

[0003] However, in actual engineering projects, borehole layout is constrained by multiple factors such as surveying costs, terrain conditions, surrounding environment, and construction conditions. Borehole data generally exhibits sparse spatial distribution and uneven density of points, making it impossible to form high-density, full-coverage continuous observation data. Traditional geological modeling methods, limited by sparse and discrete data foundations, can usually only rely on borehole points to conduct two-dimensional profile geological interpretation or complete the fitting and reconstruction of stratigraphic layers in local areas. They are unable to restore the continuous undulating morphology of the stratigraphic interface across the entire region and cannot form a complete and continuous three-dimensional geological spatial representation.

[0004] At the same time, existing 3D geological modeling technology has two significant shortcomings: First, it lacks the integration and utilization of high-dimensional geological attribute parameters. Most modeling work only focuses on the modeling and expression of stratigraphic interfaces and lithological categories, and cannot achieve a unified mapping of multi-dimensional geological parameters and 3D space, resulting in limited information carrying capacity of the model. Second, the model expression form has inherent limitations. Traditional 3D geological models mostly use surface models or solid boundary models, which can only represent the boundary morphology of strata and cannot perform fine-grained attribute control of each spatial unit inside the model, making it difficult to support subsequent spatial data query, statistical analysis, risk quantification assessment, and 3D spatial calculation. Summary of the Invention

[0005] The purpose of this invention is to provide a three-dimensional geological voxel modeling method based on borehole data. This method addresses the shortcomings of existing three-dimensional geological modeling methods in engineering geological survey scenarios with sparse borehole data. These methods typically focus on stratigraphic interface reconstruction, only constructing surface morphology-based stratigraphic boundary models. They lack a quantifiable carrier for irregular subsurface units, cannot achieve refined voxel representation of subsurface space, struggle to accurately couple multidimensional geological attributes with spatial units, and have poor compatibility with general volumetric data systems, making it difficult to support high-precision spatial analysis and quantitative engineering evaluation. This invention achieves integrated three-dimensional representation of stratigraphic structure and multidimensional geological attributes through continuous stratigraphic reconstruction, three-dimensional voxel mesh construction, multi-attribute spatial interpolation, and database mapping. This comprehensively improves the spatial integrity, attribute refinement, and engineering application value of geological modeling.

[0006] In a first aspect, the present invention provides a method for modeling a three-dimensional geological voxel model based on borehole data, comprising the following steps: Obtain geological strata borehole data for the target area, preprocess the geological strata borehole data to obtain an independent set of layered geological strata borehole points; For each geological stratum borehole set, based on the regionalized variable theory, the smooth and continuous ground plane of each stratum is reconstructed; Based on the upper and lower boundary constraints of the multi-layer continuous stratigraphic interface, the target 3D modeling space is uniformly and regularly divided into voxels to construct a 3D voxel mesh covering the entire domain and generate a set of 3D voxel units. The stratigraphic affiliation of all voxel units is determined by spatial location matching to form a 3D geological voxel skeleton model. A hierarchical adaptive interpolation algorithm is adopted. Based on the geological attribute parameters in the borehole data of each geological stratum, the multidimensional geological attribute fine interpolation is performed on each voxel unit in the three-dimensional geological voxel skeleton model to obtain the complete geological attribute parameters corresponding to each voxel. A one-to-one mapping relationship between the three-dimensional geological voxel model and the structured database is established using the three-dimensional integer index of the voxel as the unique associated primary key.

[0007] According to the three-dimensional geological voxel modeling method based on borehole data provided by the present invention, the geological strata borehole data includes spatial coordinate information under a unified geographic coordinate system, strata layering parameters, and multi-dimensional geological attribute parameters. The spatial coordinate information is the original collected spherical latitude, longitude, and elevation coordinates, and the strata layering parameters are the strata numbers corresponding to each stratum.

[0008] According to the three-dimensional geological voxel modeling method based on borehole data provided by the present invention, the preprocessing of the geological strata borehole data includes spatial coordinate normalization, outlier removal, missing parameter completion, and clustering and aggregating the geological strata borehole data by stratum number. The method for spatial coordinate normalization of geological strata borehole data includes: The original collected spherical latitude and longitude coordinates are converted into three-dimensional Euclidean engineering rectangular coordinates under a unified projection coordinate system through the Gauss-Kruger projection algorithm, thus constructing a unified spatial reference system for the entire domain. The method for outlier removal from geological strata borehole data includes: Using the same observation parameter from all borehole data within the same geological stratum as a statistical sample, a single-parameter observation sequence is constructed, and the sample mean of the single-parameter observation sequence is calculated. with standard deviation When a single observation parameter satisfy When determining single observation parameters If the entire parameter corresponding to the observation point is an outlier, the outlier will be removed directly. The method for completing missing parameters in geological strata borehole data includes: After outlier removal from the geological strata borehole data, the geological strata borehole data with complete spatial coordinate information and stratigraphic number information are defined as valid geological strata borehole data. In cases where one or more geological attribute parameters are missing in effective geological formation borehole data, a parameter-by-parameter independent completion strategy is adopted to complete the missing parameters. The method for clustering and classifying geological strata borehole data by strata number includes: Based on the stratigraphic numbering of the geological strata borehole data, the entire dataset is hierarchically clustered according to the stratigraphic number to obtain multiple independent sets of hierarchical geological strata borehole points.

[0009] According to the three-dimensional geological voxel modeling method based on borehole data provided by the present invention, the method for completing missing parameters using a parameter-by-parameter independent completion strategy includes: Based on the principle of spatial proximity within the same stratum, for each missing geological attribute parameter, the same geological attribute parameter that is spatially adjacent within the same stratum is selected from the effective geological stratum borehole data, and the missing geological attribute parameter is supplemented by the measured mean value of the corresponding geological attribute parameter.

[0010] According to the three-dimensional geological voxel modeling method based on borehole data provided by the present invention, the method for reconstructing smooth and continuous strata for each geological stratum based on regionalized variable theory for each set of borehole points includes: For each stratum, the stratum interface elevation is defined as a stratum variable. A spherical semi-variogram is used to fit the spatial variation law of the stratum interface elevation to construct a continuous stratum interface function. The optimal interpolation weights are calculated based on the constructed continuous stratigraphic interface function and the Lagrange multipliers obtained by constraining the Kriging equations. For any point to be valued, the estimated elevation of the strata corresponding to the plane where the point is located is obtained by weighted fitting of the elevations of the effectively occurring strata borehole data in the neighborhood according to the optimal interpolation weights.

[0011] According to the three-dimensional geological voxel modeling method based on borehole data provided by the present invention, the method for uniformly and regularly partitioning the target three-dimensional modeling space into voxels includes: Based on the preset modeling accuracy, the target 3D modeling space is uniformly divided along three axes to construct a 3D voxel mesh that covers the entire domain. This generates an ordered set of 3D voxel units with uniform size. Each voxel unit corresponds to a unique 3D integer index, providing a standardized spatial primary key for standardized format file storage and database mapping.

[0012] According to the three-dimensional geological voxel modeling method based on borehole data provided by the present invention, the method for determining the stratigraphic affiliation of all voxel units through spatial location matching includes: Using the elevation of the center point of the voxel unit as the core criterion, the elevation range of the upper and lower strata corresponding to the plane where the center point of the voxel unit is located is matched, and the stratigraphic affiliation of all voxel units is determined one by one to construct a three-dimensional geological voxel skeleton model.

[0013] According to the three-dimensional geological voxel modeling method based on borehole data provided by the present invention, the hierarchical adaptive interpolation algorithm includes the Kriging interpolation algorithm and the inverse distance weighted interpolation algorithm. The Kriging interpolation algorithm is used for geological parameters with strong spatial continuity and structural stability, including mechanical and structural parameters.

[0014] According to the three-dimensional geological voxel modeling method based on borehole data provided by the present invention, inverse distance weighted interpolation is used for geological parameters with strong spatial abrupt changes and significant randomness. Geological parameters with strong spatial abrupt changes and significant randomness include hydrological and topographic parameters.

[0015] The three-dimensional geological voxel modeling method based on borehole data provided by the present invention includes the following method for establishing a one-to-one mapping relationship between the three-dimensional geological voxel model and the structured database using the three-dimensional integer index of the voxel as the unique associative primary key: Using the voxel three-dimensional integer index as the globally unique association key, a field-level one-to-one mapping is established between each voxel unit and the geological strata borehole data in the structured database; A hierarchical storage architecture is adopted. The three-dimensional geological voxel model only stores the three-dimensional integer index of voxels and stratigraphic information in a lightweight manner. All multi-dimensional geological attribute parameters are stored independently in the structured database as external plug-ins, forming a three-dimensional geological voxel model with decoupled structure and attributes and the ability to be extended to high-dimensional geological data.

[0016] Compared with the prior art, the beneficial effects of the present invention are as follows: Compared to traditional stratigraphic interface models, which can only represent stratigraphic boundary morphology through curved surfaces and cannot discretize the interior of three-dimensional space, nor support spatial unit-level calculation and analysis, this invention provides a three-dimensional geological voxel modeling method based on borehole data. Through uniform and regular voxel partitioning, the entire underground space is transformed into a standardized set of voxel units, supporting unit-by-unit spatial query, statistical analysis, and numerical calculation. It has computable and analyzable spatial quantization capabilities, and the constructed three-dimensional voxel model can realize a refined unit-based quantitative expression of underground space.

[0017] Existing modeling methods often attach geological attributes to interfaces or discrete measurement points, which cannot accurately correspond to spatial units. This invention can carry multiple types of geological attribute parameters on a single voxel unit, and achieve deep integration of stratigraphic structure and geological attributes at the smallest spatial unit. This enables the integrated coupling and bearing of stratigraphic structure and multi-dimensional geological attributes, effectively improving the information carrying capacity of the three-dimensional geological model and supporting refined engineering geological evaluation and risk analysis.

[0018] Regular voxel grids and GIS raster data models share a natural structural similarity, facilitating data conversion and algorithm reuse. The three-dimensional integer index structure of voxels is consistent with the organizational logic of two-dimensional raster data and three-dimensional voxel raster data, enabling convenient data interaction and format conversion with geographic information systems. It also allows for the direct reuse of mature spatial analysis, visualization, and numerical calculation algorithms, reducing the adaptation costs for subsequent engineering applications.

[0019] Existing solutions often embed all attribute data into the model file, leading to bloated files, low reading efficiency, and poor scalability. This invention uses a voxel 3D integer index as the unique primary key to construct a mapping mechanism between the voxel model and the structured database, enabling efficient management and flexible access to geological data. It adopts a hierarchical storage model of "lightweight structural model + structured attribute database": stratigraphic structure information is stored in the voxel model file, while multi-dimensional geological attribute parameters are stored in the backend database. This supports rapid storage, retrieval, statistical analysis, and visualization of geological data, and is adaptable to the needs of online web systems.

[0020] The model output is highly standardized and has strong engineering adaptability. It supports direct export of voxel models in the .vti standard format, which is compatible with mainstream geological visualization, numerical simulation and 3D modeling platforms without secondary format conversion, effectively improving the reusability and scene adaptability of the modeling results.

[0021] The entire method is adapted to engineering survey scenarios with sparse borehole data. Through the optimization of the entire process, including layered preprocessing, ground plane reconstruction, voxel structured subdivision, adaptive attribute interpolation, and database mapping and storage, it can complete the construction of high-precision three-dimensional geological voxel models under the conditions of sparse borehole distribution and uneven point location, effectively improving the integrity and accuracy of modeling under complex geological conditions. Attached Figure Description

[0022] To more clearly illustrate the technical solutions in this invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of this invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0023] Figure 1This is a flowchart illustrating a three-dimensional geological voxel modeling method based on borehole data provided by the present invention. Detailed Implementation

[0024] The embodiments of the present invention will be described in further detail below with reference to the accompanying drawings and examples. The following examples are for illustrative purposes only and should not be construed as limiting the scope of the invention.

[0025] In the description of this invention, it should be noted that the terms "center," "longitudinal," "lateral," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," and "outer," etc., indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings. They are used only for the convenience of describing the invention and for simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, they should not be construed as limitations on the invention. Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance.

[0026] In the description of this invention, it should be noted that, unless otherwise explicitly specified and limited, the terms "connected" and "linked" should be interpreted broadly. For example, they can refer to a fixed connection, a detachable connection, or an integral connection; they can refer to a mechanical connection or an electrical connection; they can refer to a direct connection or an indirect connection through an intermediate medium. Those skilled in the art can understand the specific meaning of the above terms in this invention based on the specific circumstances.

[0027] like Figure 1 As shown in the figure, an embodiment of the present invention provides a method for modeling a three-dimensional geological voxel model based on borehole data, which includes the following steps: S1, Obtain geological strata borehole data of the target area, preprocess the geological strata borehole data to obtain an independent set of layered geological strata borehole points; S2, for each geological stratum borehole point set, based on the regional variable theory, reconstructs the smooth and continuous ground plane of each stratum; S3, based on the upper and lower boundary constraints of the multi-layer continuous stratigraphic interface, uniformly and regularly divides the target three-dimensional modeling space into voxels, constructs a three-dimensional voxel mesh covering the entire domain, and generates a set of three-dimensional voxel units; the stratigraphic attribution of all voxel units is determined by spatial location matching, forming a three-dimensional geological voxel skeleton model. S4 employs a hierarchical adaptive interpolation algorithm. Based on the geological attribute parameters in the borehole data of each geological stratum, it performs refined interpolation of multidimensional geological attributes on each voxel unit in the three-dimensional geological voxel skeleton model to obtain the complete geological attribute parameters corresponding to each voxel. S5 establishes a one-to-one mapping relationship between the three-dimensional geological voxel model and the structured database, using the voxel three-dimensional integer index as the unique associated primary key.

[0028] This invention sequentially performs continuous stratigraphic reconstruction, 3D voxel mesh construction, multi-attribute spatial interpolation, and database mapping and association steps to achieve an integrated 3D representation of stratigraphic structure and multi-dimensional geological attributes. This addresses the technical problems encountered in engineering geological survey scenarios with sparse borehole data. Existing 3D geological modeling generally focuses on stratigraphic interface reconstruction, only capable of constructing surface morphology-based stratigraphic boundary models. It lacks a quantifiable carrier for irregular subsurface units, cannot achieve refined voxel representation of subsurface space, struggles to accurately couple multi-dimensional geological attributes with spatial units, and suffers from poor compatibility between model formats and general volumetric data systems, making it difficult to support high-precision spatial analysis and quantitative engineering evaluation. On one hand, it can completely reconstruct the spatial distribution of subsurface strata based on sparse discrete borehole data, improving the spatial integrity of the geological model. On the other hand, it utilizes voxel interpolation to achieve refined assignment of geological parameters across the entire domain. Combined with the mapping and binding of voxel units to the attribute database, this ensures the accuracy of geological attribute characterization while facilitating dynamic management of geological data. Ultimately, it effectively enhances the detail representation of 3D geological models and their adaptability to various engineering geological analysis scenarios.

[0029] In step S1 of this embodiment of the invention, geological strata borehole data of the target area are obtained, and the geological strata borehole data are preprocessed to obtain an independent set of layered geological strata borehole points.

[0030] In step S1 of this embodiment of the invention, the geological strata borehole data includes spatial coordinate information under a unified geographic coordinate system, strata layering parameters, and multidimensional geological attribute parameters.

[0031] In real-world engineering geological survey scenarios, the geological strata borehole data are non-uniformly discrete survey data, lacking continuous profile observation data, and characterized by low sample size and sparse distribution.

[0032] The spatial coordinate information is the original spherical latitude, longitude, and elevation coordinates collected.

[0033] The stratigraphic stratification parameter is the stratigraphic number corresponding to each stratum.

[0034] The borehole data of a single geological stratum are constructed into a multi-dimensional parameter vector to achieve a standardized and quantitative definition of the geological stratum borehole data. The expression for a single geological stratum borehole data is as follows:

[0035] In the formula: For the first Complete observation vector for each borehole; These are the three-dimensional spatial coordinates of the borehole survey points; Stratigraphic strata are numbered for use in stratigraphic unit classification; For the first Each drill hole corresponds to The class can be extended with geological attribute parameters. This represents the total number of geological attribute parameters, which can be adaptively adjusted according to engineering modeling needs.

[0036] In one specific implementation, five strata exist within the target area, and several sparse geological strata borehole data were collected for each stratum. Each geological stratum borehole data includes the actual geographical coordinates, stratum number, and eight attributes: wave height, seaward slope, reef limestone strength, overburden thickness, stress, water inflow, deformation, and rock permeability.

[0037] Preprocessing the geological strata borehole data yields an independent set of layered geological strata borehole points, which can eliminate noise errors in discrete geological strata borehole data and interference from cross-strata data.

[0038] Raw field geological borehole datasets commonly suffer from technical defects such as inconsistent spatial benchmarks, abnormal observation noise, missing parameters, and mixed stratigraphic samples. Direct interpolation can lead to stratigraphic interface distortion, inaccurate attribute fitting, and disordered stratigraphic logic. This step preprocesses the geological borehole data to standardize, normalize, and reconstruct the original dataset, eliminating systematic and random errors at the data source and ensuring the accuracy of subsequent modeling and consistency with stratigraphic logic.

[0039] In step S1 of this embodiment of the invention, the preprocessing of the geological strata borehole data includes spatial coordinate normalization, outlier removal, missing parameter completion, and clustering and aggregating the geological strata borehole data by stratum number.

[0040] In step S1 of this embodiment of the invention, the method for spatial coordinate normalization of geological strata borehole data includes: converting the original collected spherical latitude and longitude coordinates into three-dimensional Euclidean engineering rectangular coordinates under a unified projection coordinate system through the Gauss-Kruger projection algorithm, constructing a unified spatial reference system for the entire domain, and eliminating spherical coordinate distortion and spatial reference deviation.

[0041] The original survey coordinates are mostly spherical latitude and longitude coordinates, which suffer from spherical distortion and cannot be directly used for 3D Euclidean space modeling. A coordinate transformation must be performed using the Gauss-Kruger projection algorithm to achieve an accurate mapping from spherical geographic coordinates to planar engineering rectangular coordinates. The transformation model is as follows:

[0042] In the formula: For projection transformation operators; For original longitude, For the original latitude, Original elevation depth; This represents the three-dimensional Euclidean space coordinates under the unified projected coordinate system after transformation. After transformation, all borehole samples are uniformly incorporated into the same spatial reference system, eliminating coordinate distortion and reference deviation.

[0043] In step S1 of this embodiment of the invention, the method for outlier removal from geological strata borehole data includes: Outlier identification was performed using the 3σ criterion combined with geological constraint thresholds. A single-parameter observation sequence was constructed by using the same observation parameter from all borehole data within the same geological stratum as a statistical sample, and the sample mean and standard deviation of the single-parameter observation sequence were calculated.

[0044] When the single observation parameter satisfies When this happens, the entire parameter corresponding to that observation point is determined to be an outlier and is directly removed. This represents the number of samples in a single-parameter observation sequence. This represents the sample mean of a single-parameter observation sequence. Standard deviation of a single-parameter observation sequence.

[0045] When the single observation parameter satisfies When an observation point is identified as having slight anomalies consistent with local geological variations, a weighted correction is applied using geological and stratigraphic borehole data from neighboring observation points in the same stratum to eliminate the interference of random observation errors. Simultaneously, duplicate and redundant geological and stratigraphic borehole data are removed to ensure the uniqueness of measurement points and the validity of parameters within the dataset.

[0046] In step S1 of this embodiment of the invention, the method for completing missing parameters in geological strata borehole data includes: After outlier removal, geological and stratigraphic borehole data with complete spatial coordinate information and stratigraphic number information are defined as valid geological and stratigraphic borehole data.

[0047] To address the issue of missing one or more geological attribute parameters in effective geological formation borehole data, a parameter-by-parameter independent completion strategy is adopted to complete the missing parameters.

[0048] In step S1 of this embodiment of the invention, the method of using a parameter-by-parameter independent completion strategy to complete missing parameters includes: Based on the principle of spatial proximity within the same stratum, for each missing geological attribute parameter, the same geological attribute parameter spatially adjacent within the same stratum is selected from the effective geological stratum borehole data, and the missing parameter is filled in using the measured mean of the corresponding geological attribute parameter. The formula for filling in a single missing parameter is as follows:

[0049] In the formula: To fill in missing parameters with values; The number of borehole data for adjacent geological strata in the same formation that participated in this supplementary calculation; The measured value of the parameter corresponding to the borehole data of the j-th adjacent geological stratum.

[0050] By using the above method, each missing parameter is completed independently according to the parameter dimension. Multiple missing parameters of borehole data of the same geological stratum can be calculated independently without interference. This can preserve the original survey sample size to the maximum extent while ensuring the geological rationality of the data.

[0051] In step S1 of this embodiment of the invention, the method for clustering and collecting geological strata borehole data by stratum number includes: Based on the embedded stratigraphic numbers in the geological strata borehole data, hierarchical clustering was performed on the entire dataset according to the stratigraphic numbers, resulting in multiple independent stratigraphic point sets. The clustering result expression is as follows:

[0052]

[0053] In the formula: D is the global borehole dataset; Let m be the set of independent boreholes corresponding to the t-th stratum; m is the total number of strata. For the first Complete observation vector for each borehole; Stratigraphic layers are numbered. By aggregating data into layers, physical isolation of data from different stratigraphic layers is achieved, avoiding interference from cross-strata geological differences on interface reconstruction.

[0054] The borehole data are divided according to the stratigraphic number:

[0055] Following the previous example, the Gauss-Kruger projection algorithm was used to uniformly convert the original spherical latitude and longitude coordinates of all boreholes into the same engineering rectangular coordinate system, thus completing the spatial reference normalization. The 3σ criterion combined with the geological constraint threshold was used to identify parameter outliers, and outlier samples were directly removed. For minor local anomalies, the weighted correction of the neighboring measurement points in the same stratum was used. For the missing geological attribute parameters in the valid borehole samples, the mean value of the corresponding parameters of the spatially adjacent measured boreholes in the same stratum was selected to complete the adaptive completion. After data cleaning, all borehole points were hierarchically clustered and classified into independent borehole point sets corresponding to 5 strata based on the initial stratum number labels of the boreholes.

[0056] In step S2 of this embodiment of the invention, for each geological stratum borehole point set, the smooth and continuous strata of each stratum are reconstructed based on the regional variable theory.

[0057] Even after preprocessing the geological strata borehole data, the borehole measurement points corresponding to each stratum are still sparse and discrete. Relying solely on these scattered measurement points is insufficient to fully reconstruct the continuous undulating shape of the underground strata interface.

[0058] In step S2 of this embodiment of the invention, the method for reconstructing smooth and continuous strata for each geological stratum based on regionalized variable theory includes: For each stratum, the stratum interface elevation is defined as a stratum variable. A spherical semi-variogram is used to fit the spatial variation law of the stratum interface elevation to construct a continuous stratum interface function. The optimal interpolation weights are calculated based on the constructed continuous stratigraphic interface function and the Lagrange multipliers obtained by constraining the Kriging equations. For any point to be valued, the estimated elevation of the strata corresponding to the plane where the point is located is obtained by weighted fitting of the elevations of the effectively occurring strata borehole data in the neighborhood according to the optimal interpolation weights.

[0059] In step S2 of this embodiment of the invention, the Kriging interpolation algorithm is used. By utilizing the spatial autocorrelation characteristics of the stratum elevation, the spatial variation law of the elevation is quantified by fitting a semi-variogram function, and the optimal interpolation weight is solved. The measured elevation data of discrete boreholes are extrapolated into a continuous and smooth stratum, which solves the problem that the stratum is prone to faults and distortions under sparse borehole conditions, and ensures the continuity and accuracy of the stratum reconstruction.

[0060] For the k-th stratum, the stratum interface elevation is defined as a regional variable, and a continuous stratum interface function is constructed:

[0061] In the formula: For the first Arbitrary plane coordinates of strata The interface elevation value at the location enables a continuous mapping from a two-dimensional plane to an elevation dimension.

[0062] The spatial variation of formation elevation is fitted using a spherical semivariogram to accurately characterize the spatial structure and randomness of formation variables.

[0063] In the formula: h is the Euclidean distance between the two points in space; denoted as nugget value, representing micro-level random variation; C represents structural variance, representing regional structural variation; and a represents range, representing the scope of the variable's spatial autocorrelation.

[0064] To satisfy the interpolation unbiasedness and minimum variance, the optimal interpolation weights are calculated by solving the Kriging equations based on the semivariogram and Lagrange multipliers. The constructed Kriging equations are as follows:

[0065] In the formula: Given the known half-variogram values ​​among the samples; denoted as the semivariogram value between the known sample and the point to be estimated; μ is the Lagrange multiplier, ensuring that the unbiased interpolation constraint holds.

[0066] For any point to be estimated in the plane The optimal unbiased estimate of the formation elevation is obtained by weighted fitting of the effective borehole elevations in the neighborhood:

[0067] In the formula: This is the optimal estimate of the elevation of the point to be estimated; The number of valid borehole samples participating in the interpolation for the k-th stratum; The measured elevation values ​​of the known borehole points; The optimal interpolation weights are for each known sample.

[0068] For each set of geological strata borehole points, the geological strata borehole data is continuously extrapolated based on geostatistical interpolation algorithms to reconstruct smooth and continuous ground planes for each stratum, thus making up for the defects of discrete and discontinuous sparse geological strata borehole data.

[0069] The above algorithm completes interpolation calculations layer by layer, ultimately generating a continuous, smooth, and realistic geological stratum model. By supplementing the bottom constraint interface, a complete upper and lower boundary constraint system of the strata can be formed, effectively solving the problems of stratum faults, distortions, and discontinuities in sparse discrete borehole data.

[0070] When conventional methods rely on borehole interpolation to generate stratigraphic interfaces, the borehole exploration depth is limited, there is no effective stratigraphic interface data in the deep underground part of the model, and the bottom of the model is in a completely open state.

[0071] The underlying constraint interface is essentially an artificially defined lowest reference plane with a fixed elevation, locking the maximum downward range of the three-dimensional mesh along the Z-axis, forming a closed three-dimensional modeling space.

[0072] Following the previous example, in a specific implementation, for each set of borehole points completed in groups, the ordinary Kriging interpolation algorithm is used one by one. Combined with the spatial distribution law of stratum elevation, the corresponding semi-variogram model is fitted, the optimal interpolation weight is solved, and the point-by-point extrapolation calculation of stratum elevation is completed on the global plane grid. A smooth and continuous single-layer stratum surface that fits the real undulation of the stratum is generated. The interface reconstruction of all strata is completed in sequence to form a multi-layer continuous stratum interface system.

[0073] In step S3 of this embodiment of the invention, based on the upper and lower boundary constraints of the multi-layer continuous stratigraphic interface, the target three-dimensional modeling space is uniformly and regularly divided into voxels to construct a three-dimensional voxel grid covering the entire domain and generate a set of three-dimensional voxel units; the stratigraphic attribution of all voxel units is determined by spatial location matching to form a three-dimensional geological voxel skeleton model.

[0074] Traditional stratigraphic interface models can only express the morphology of stratigraphic boundaries through curved surfaces, and cannot perform unitized quantitative representation of the subsurface three-dimensional space. This step, based on the upper and lower boundary constraints of multiple continuous stratigraphic interfaces, performs uniform and regular voxel subdivision of the target three-dimensional modeling space to construct a three-dimensional voxel mesh covering the entire domain; through spatial location matching, the stratigraphic affiliation of all voxel units is determined, forming a three-dimensional geological voxel skeleton model, realizing the transformation from a two-dimensional stratigraphic interface to a three-dimensional solid geological model.

[0075] In step S3 of this embodiment of the invention, the following step is also included: Based on the global borehole coordinates and formation interface extreme parameters, a closed 3D modeling spatial domain is defined. The spatial domain constraint formula is as follows:

[0076] In the formula: These are the extreme values ​​of the three-dimensional coordinates of the modeling area, ensuring that the modeling space completely covers the entire geological distribution range of the target area.

[0077] In step S3 of this embodiment of the invention, the method for uniformly and regularly partitioning the target three-dimensional modeling space into voxels includes: uniformly partitioning the target three-dimensional modeling space along three axes according to a preset modeling precision, constructing a three-dimensional voxel mesh that covers the entire domain, generating an ordered set of three-dimensional voxel units with uniform size, and each voxel corresponding to a unique three-dimensional integer index, providing a standardized spatial primary key for standardized format file storage and database mapping.

[0078] Standardized format files include common structured voxel data formats such as .vti format files, which are compatible with mainstream geological visualization, numerical simulation and 3D modeling platforms. No secondary format conversion is required, which effectively improves the reusability and scene adaptability of modeling results.

[0079] Based on the complexity of stratigraphic undulations and the requirements for modeling accuracy, the triaxial uniform partition size is set. The closed spatial domain is uniformly partitioned in three dimensions to determine the voxel resolution of the voxel model. The calculation formula is as follows:

[0080] In the formula This represents the average distance between adjacent boreholes. This represents the proportionality coefficient, which is usually taken as... to .

[0081] Based on a defined voxel resolution (Δx, Δy, Δz), the number of triaxial mesh elements is calculated by uniformly partitioning the closed spatial domain in three dimensions:

[0082] In the formula: For the floor operator, , , These represent the number of voxel mesh elements in the x, y, and z axes, respectively.

[0083] Finally, an ordered set of voxel units is generated:

[0084] In the formula: These are the voxel triaxial index numbers, which together constitute the three-dimensional integer index of the voxel, uniquely identifying the spatial position of the voxel in the three-dimensional mesh, consistent with the native indexing system of the subsequent .vti format.

[0085] In step S3 of this embodiment of the invention, the method of determining the stratigraphic affiliation of all voxel units by spatial location matching includes: using the elevation of the center point of the voxel unit as the core criterion, matching the elevation range of the upper and lower ground planes corresponding to the plane where the center point of the voxel unit is located, determining the stratigraphic affiliation of all voxel units one by one, and constructing a three-dimensional geological voxel skeleton model with complete topological relationships and clear layering logic.

[0086] It supports direct export of standardized .vti format 3D geological voxel model files. The .vti format 3D geological voxel model has built-in 3D integer indexes and stratigraphic numbers for voxel, enabling lightweight and structured stratigraphic model output.

[0087] In step S3 of this embodiment of the invention, the step of: accurately solving the three-dimensional coordinates of each voxel center is performed using the voxel center point as the unit feature representation point.

[0088]

[0089] In the formula The x-axis represents the engineering rectangular coordinates of the center point of the voxel with index number (a,b,c). The representation of the engineering rectangular coordinates of the voxel center point with the 3D index number (a,b,c) along the Y-axis. The Z-axis represents the engineering rectangular coordinates of the center point of the voxel with index number (a,b,c).

[0090] In step S3 of this embodiment of the invention, the formula for determining the stratigraphic affiliation of each voxel unit by matching the elevation range of the ground plane above and below the voxel center point elevation as the core criterion is as follows:

[0091] Voxel elements that meet the above criteria are uniformly assigned the value of the k-th stratum. The stratum labeling formula is as follows:

[0092] By traversing all voxel units, the stratigraphic classification and labeling of the entire region are completed, ultimately forming a three-dimensional geological voxel skeleton model with accurate topological relationships and clear layering logic.

[0093] After constructing the voxel skeleton model, this method supports the direct output of standardized 3D voxel stratigraphic model files, and can adaptively export .vti format model files. The .vti format, as a common structured voxel data format for visualization tools, can completely preserve the model's 3D mesh topology, stratigraphic boundary constraints, and spatial coordinate reference information without secondary format conversion. It can directly adapt to the import and parsing requirements of various geological visualization, numerical simulation, and 3D modeling platforms, significantly improving model reusability and scene adaptability flexibility. It can independently output pure stratigraphic structure models, adapting to engineering application scenarios that only require the stratigraphic skeleton structure and do not need geological attributes, thus enriching the output formats and usage dimensions of modeling results.

[0094] Continuing from the previous example, in a specific implementation, based on the upper and lower boundary constraints of multi-layered continuous stratigraphic interfaces and combined with preset voxel precision dimensions, the target 3D modeling space is uniformly and regularly partitioned to generate a set of 3D voxel units that cover the entire domain. By comparing the elevation of the voxel center point with the elevation of the stratigraphic interface using a threshold comparison, the stratigraphic affiliation of each voxel is accurately determined, completing the stratigraphic classification and labeling of all voxel units, and constructing a structurally complete and precisely boundaryed 3D geological voxel skeleton model.

[0095] In step S4 of this embodiment of the invention, a hierarchical adaptive interpolation algorithm is used to perform fine interpolation of multidimensional geological attributes on each voxel unit in the three-dimensional geological voxel skeleton model based on the geological attribute parameters in the borehole data of each geological stratum, so as to obtain the complete geological attribute parameters corresponding to each voxel.

[0096] Voxel skeleton models only possess spatial coordinate information and stratigraphic layering information, lacking quantitative parameters for geological attributes, thus failing to achieve digital and refined characterization of geological features. Addressing the geological characteristics of varying spatial continuity, variability sensitivity, and abrupt change features across different geological attributes, this invention employs a stratigraphically independent interpolation + dual-algorithm adaptive adaptation technique. This avoids interference from cross-strata parameters, matches the optimal interpolation algorithm to the spatial variation characteristics of different attributes, and performs high-precision assignment of multidimensional geological attributes on a voxel-by-voxel basis, constructing voxel-level multidimensional geological attribute vectors to achieve deep coupling between the model's geometric structure and geological attributes.

[0097] In step S4 of this embodiment of the invention, the following step is also included: A multi-dimensional geological attribute vector is constructed for each voxel unit to achieve integrated encapsulation of multiple parameters:

[0098] In the formula A multidimensional geological attribute vector for voxel units; For each voxel center point, there are n types of expandable geological attribute parameters, and the parameter types and quantities correspond one-to-one with the borehole observation vectors.

[0099] In step S4 of this embodiment of the invention, the hierarchical adaptive interpolation algorithm includes the Kriging interpolation algorithm and the inverse distance weighted interpolation algorithm. The Kriging interpolation algorithm is used for geological parameters with strong spatial continuity and structural stability. These geological parameters include mechanical and structural parameters. The Kriging interpolation algorithm used is as follows:

[0100] In the formula: Assign attribute values ​​to voxel interpolation; These are measured attribute values ​​from boreholes in the same stratum. This is the optimal interpolation weight for this attribute.

[0101] In step S4 of this embodiment of the invention, an inverse distance weighted interpolation algorithm is used for geological parameters with strong spatial abrupt changes and significant randomness. This algorithm takes into account local geological abrupt change characteristics and uses only borehole data from the same stratum for interpolation calculation, thus avoiding interference from geological parameters across strata. Geological parameters with strong spatial abrupt changes and significant randomness include hydrological and topographical parameters.

[0102] The inverse distance weighted interpolation algorithm is as follows:

[0103]

[0104] In the formula: The three-dimensional Euclidean distance between the voxel and the borehole location; It is the distance decay power exponent, used to adjust the decay rate of spatial influence weights.

[0105] After interpolation is performed on each voxel unit, an expandable attribute vector containing multiple geological parameters is generated, which fully records the multidimensional geological feature information of a single voxel unit and supports subsequent geological analysis and engineering evaluation.

[0106] Through the above-mentioned hierarchical adaptive interpolation algorithm, the voxel attributes of the whole domain are assigned layer by layer and parameter by parameter, which accurately adapts to the spatial distribution of different geological parameters, avoids cross-layer parameter interference, and enables each voxel unit to have complete and accurate multi-dimensional geological attribute parameters, and finally completes the construction of a high-precision three-dimensional geological voxel model integrating structure and attributes.

[0107] Following the previous example, in a specific implementation, based on a stratigraphic adaptive interpolation strategy, appropriate spatial interpolation algorithms are used for the spatial distribution characteristics of different types of geological attributes. Valid borehole attribute data for the corresponding strata are called upon to perform refined interpolation and assignment of the eight geological attribute parameters for each voxel unit, ensuring that each voxel unit carries complete multidimensional geological attribute information, thus achieving the integrated fusion of voxel model structure and geological attributes. Furthermore, after the voxel skeleton model is constructed, a .vti format stratigraphic model containing only stratigraphic structure and spatial topology information can be independently exported, achieving lightweight and standardized output of the structured stratigraphic model.

[0108] In step S5 of this embodiment of the invention, the three-dimensional integer index corresponding to the voxel unit is used as the unique associated primary key to construct the voxel global code, establish a one-to-one mapping relationship between the three-dimensional geological voxel model and the structured database, and store the voxel spatial coordinate information, stratigraphic number, and multi-dimensional geological attribute parameters in the database in layers to realize the two-way linkage between the visualization of the lightweight voxel model and the structured storage of geological attributes.

[0109] To address the technical shortcomings of fragmented, weakly correlated, and unstructured data management of spatial coordinate information, stratigraphic layering information, and multidimensional geological attribute information in 3D voxel models, this step constructs a globally unique voxel encoding mechanism and a one-to-one mapping model for multidimensional parameters. This enables precise association between 3D spatial entity units and structured database fields, and completes the standardization, structuring, and integrated storage of all modeling data, achieving digital, quantifiable, and traceable management of 3D geological models.

[0110] In step S5 of this embodiment of the invention, the method for constructing global voxel encoding and establishing a one-to-one mapping relationship between the three-dimensional geological voxel model and the structured database includes: Based on three-dimensional voxel grid coordinates, a mapping from continuous spatial coordinates to voxel indices is established to achieve precise voxel positioning. Through this mapping model, the continuous three-dimensional geological space is discretized into a regular integer index raster system. Each unique integer index (a, b, c) precisely corresponds to a voxel unit at a fixed spatial location in the work area, realizing the binding between physical spatial location and discrete raster index. This provides a standardized spatial carrier for subsequent stratigraphic numbering and built-in storage of multi-dimensional geological attributes, and is fully compatible with the structured storage logic of .vti voxel files.

[0111] In step S5 of this embodiment of the invention, the structured database adopts a hierarchical storage mapping mechanism. The vti format voxel model only stores the three-dimensional integer index and stratigraphic number corresponding to the voxel unit in a lightweight manner. All multi-dimensional geological attribute parameters are bound to the three-dimensional integer index corresponding to the voxel unit and then stored in the database. The database storage fields can be expanded as needed to achieve bidirectional adaptation between lightweight model visualization and refined structured management of geological attributes.

[0112] Based on the underlying storage characteristics of .vti format voxel models, a single voxel unit is suitable for storing structured, lightweight basic identifier parameters, but it is not suitable for directly storing massive multidimensional geological attribute parameters in batches. Therefore, this invention constructs a hierarchical coupling and association system of lightweight built-in voxel storage and refined attribute storage in a database: using the converted voxel 3D integer index (a, b, c) as the unique association key, only the spatial index and stratigraphic number are natively written into the .vti voxel, ensuring the lightweight and highly compatible output characteristics of the voxel model; while all refined and extensible multidimensional geological attribute parameters are uniformly bound to voxel coordinates and unique numbers and stored in a structured database, achieving hierarchical adaptation between model visualization and refined data management, avoiding the problems of model redundancy, reading lag, and poor compatibility caused by embedding multiple types of parameters within voxels.

[0113] Corresponding extended multidimensional geological attribute vector Using voxel 3D integer indices as a link, each voxel is bound to its spatial coordinates and stratigraphic number, and the data is fully entered into a backend structured database for storage and management. This hierarchical storage mechanism strictly adapts to the native storage specification of the .vti format, ensuring that the exported .vti stratigraphic models are lightweight, highly compatible, and can be used independently for 3D visualization and stratigraphic structure analysis. Furthermore, by storing multi-dimensional geological attributes in the database, it enables queryable, statistical, and scalable management of geological parameters, solving the technical drawbacks of traditional models that suffer from mixed storage of structure and attributes, bloated models, and poor parameter scalability.

[0114] Meanwhile, by directly using the voxel 3D integer index as the bidirectional association primary key, the precise data linkage between the .vti lightweight structural model and the structured attribute database is established, realizing an integrated adaptive modeling system of "voxel spatial structure visualization + geological attribute refined management". This eliminates the need to generate additional global code, simplifies the modeling process, and reduces data redundancy.

[0115] In an optional implementation, in step S5 of this embodiment, the structured database employs a hierarchical storage mapping mechanism to build a web-based online geological analysis system. The .vti format stratigraphic voxel models are deployed on the front-end visualization engine for 3D display, spatial sectioning, and interactive querying of stratigraphic structures. Multidimensional geological attribute parameters are stored in the back-end database. The front-end can retrieve the complete geological attribute parameters of the corresponding voxel in real time by sending a request to the back-end via the voxel's 3D integer index. This deployment method eliminates the need to load all attribute data, significantly improving front-end loading and interaction efficiency, and supporting multi-terminal online access and collaborative analysis.

[0116] Based on the aforementioned integrated voxel indexing system, a structured database mapping operator is constructed. Using the voxel 3D integer index (a, b, c) as the unique primary key, a precise four-dimensional mapping is established for .vti voxel raster information, spatial coordinate information, stratigraphic number, and multi-dimensional geological attributes. This achieves a bidirectional and precise correspondence between the voxel model visualization data and the structured database data, completely resolving the problems of fragmented and weakly correlated data in traditional models. The mapping model is as follows:

[0117] In the formula: M is the structured mapping function; a, b, c are the voxel three-axis integer indices; The voxel space coordinate vector; The stratigraphic number to which the voxel belongs; It is a voxel-based multidimensional extensible geological attribute vector.

[0118] This invention abandons the traditional voxel global encoding method and directly uses the voxel's native three-dimensional integer index as the unique associated primary key, adapting to the .vti format native spatial identification system, simplifying the mapping logic, reducing redundant calculations, and ensuring the nativeness and accuracy of data association.

[0119] Based on the mapping relationship, a standardized database table is constructed, transforming the complete parameter information of each voxel unit into an independent structured data record. This achieves a one-to-one correspondence between the 3D model entity and the structured data in the database, completing the entire high-precision modeling process for sparse borehole 3D geological voxels. It realizes the complete transformation from raw discrete survey data to continuous interface models, 3D entity voxel models, and structured attribute databases. The fields included in the established voxel database table are shown in Table 1. Table 1. Voxel Database

[0120] By constructing a high-precision mapping relationship, the structured and integrated management of all unit data of the three-dimensional geological voxel model was realized. The entire process of sparse borehole data was accurately transformed from raw discrete survey data to continuous ground planes, three-dimensional solid voxel models, and multi-attribute structured databases, achieving a comprehensive improvement in the spatial integrity, attribute refinement, and data structuring of geological modeling.

[0121] Following the previous example, a globally unique number is assigned to each voxel unit, and a precise mapping relationship is constructed between voxel spatial coordinates, stratigraphic affiliation, multidimensional geological attributes and database fields. The complete parameter information of all voxel units is then entered into the database in batches, forming a structured, traceable, and integrated three-dimensional geological voxel model database, thus completing the entire modeling process.

[0122] The following describes a three-dimensional geological voxel modeling system based on borehole data provided by the present invention. The three-dimensional geological voxel modeling system based on borehole data described below and the three-dimensional geological voxel modeling method based on borehole data described above can be referred to and correspond to each other.

[0123] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0124] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments.

[0125] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for modeling a three-dimensional geological voxel model based on borehole data, characterized in that, Includes the following steps: Obtain geological strata borehole data for the target area, preprocess the geological strata borehole data to obtain an independent set of layered geological strata borehole points; For each geological stratum borehole set, based on the regionalized variable theory, the smooth and continuous ground plane of each stratum is reconstructed; Based on the upper and lower boundary constraints of the multi-layer continuous stratigraphic interface, the target 3D modeling space is uniformly and regularly divided into voxels to construct a 3D voxel mesh covering the entire domain and generate a set of 3D voxel units. The stratigraphic affiliation of all voxel units is determined by spatial location matching to form a 3D geological voxel skeleton model. A hierarchical adaptive interpolation algorithm is adopted. Based on the geological attribute parameters in the borehole data of each geological stratum, the multidimensional geological attribute fine interpolation is performed on each voxel unit in the three-dimensional geological voxel skeleton model to obtain the complete geological attribute parameters corresponding to each voxel. A one-to-one mapping relationship between the three-dimensional geological voxel model and the structured database is established using the three-dimensional integer index of the voxel as the unique associated primary key.

2. The three-dimensional geological voxel modeling method based on borehole data according to claim 1, characterized in that: The geological and stratigraphic borehole data includes spatial coordinate information under a unified geographic coordinate system, stratigraphic layering parameters, and multidimensional geological attribute parameters. The spatial coordinate information is the original spherical latitude, longitude, and elevation coordinates, and the stratigraphic layering parameters are the stratigraphic numbers corresponding to each stratum.

3. The three-dimensional geological voxel modeling method based on borehole data according to claim 1, characterized in that: The preprocessing of geological strata borehole data includes spatial coordinate normalization, outlier removal, missing parameter completion, and clustering by stratum number.

4. The three-dimensional geological voxel modeling method based on borehole data according to claim 3, characterized in that: The method for spatial coordinate normalization of geological strata borehole data includes: The original collected spherical latitude and longitude coordinates are converted into three-dimensional Euclidean engineering rectangular coordinates under a unified projection coordinate system through the Gauss-Kruger projection algorithm, thus constructing a unified spatial reference system for the entire domain. The method for outlier removal from geological strata borehole data includes: Using the same observation parameter from all borehole data within the same geological stratum as a statistical sample, a single-parameter observation sequence is constructed, and the sample mean of the single-parameter observation sequence is calculated. with standard deviation When a single observation parameter satisfy When determining single observation parameters If the entire parameter corresponding to the observation point is an outlier, the outlier will be removed directly. The method for completing missing parameters in geological strata borehole data includes: After outlier removal from the geological strata borehole data, the geological strata borehole data with complete spatial coordinate information and stratigraphic number information are defined as valid geological strata borehole data. To address the issue of missing one or more geological attribute parameters in effective geological formation borehole data, a parameter-by-parameter independent completion strategy is adopted to complete the missing parameters.

5. The three-dimensional geological voxel modeling method based on borehole data according to claim 3, characterized in that: The method for clustering and classifying geological strata borehole data by strata number includes: Based on the stratigraphic numbering of the geological strata borehole data, the entire dataset is hierarchically clustered according to the stratigraphic number to obtain multiple independent sets of hierarchical geological strata borehole points.

6. The method for modeling a three-dimensional geological voxel model based on borehole data according to claim 1, characterized in that, The method for completing missing parameters using a parameter-by-parameter independent completion strategy includes: Based on the principle of spatial proximity within the same stratum, for each missing geological attribute parameter, the same geological attribute parameter that is spatially adjacent within the same stratum is selected from the effective borehole data of the geological stratum, and the missing geological attribute parameter is supplemented by the measured mean value of the corresponding geological attribute parameter.

7. The method for modeling a three-dimensional geological voxel model based on borehole data according to claim 1, characterized in that, The method for reconstructing smooth and continuous ground planes for each geological stratum based on regionalized variable theory, using the borehole point set for each stratum, includes: For each stratum, the stratum interface elevation is defined as a stratum variable. A spherical semi-variogram is used to fit the spatial variation law of the stratum interface elevation to construct a continuous stratum interface function. The optimal interpolation weights are calculated based on the constructed continuous stratigraphic interface function and the Lagrange multipliers obtained by constraining the Kriging equations. For any point to be valued, the estimated elevation of the strata corresponding to the plane where the point is located is obtained by weighted fitting of the elevations of the effectively occurring strata borehole data in the neighborhood according to the optimal interpolation weights.

8. The method for modeling a three-dimensional geological voxel model based on borehole data according to claim 1, characterized in that, The method for uniformly and regularly voxelizing the target 3D modeling space includes: Based on the preset modeling accuracy, the target 3D modeling space is uniformly divided along three axes to construct a 3D voxel mesh that covers the entire domain. This generates an ordered set of 3D voxel units with uniform size. Each voxel unit corresponds to a unique 3D integer index, providing a standardized spatial primary key for standardized format file storage and database mapping.

9. The method for modeling a three-dimensional geological voxel model based on borehole data according to claim 1, characterized in that, The method for determining the stratigraphic affiliation of all voxel units through spatial location matching includes: Using the elevation of the center point of the voxel unit as the core criterion, the elevation range of the upper and lower strata corresponding to the plane where the center point of the voxel unit is located is matched, and the stratigraphic affiliation of all voxel units is determined one by one to construct a three-dimensional geological voxel skeleton model.

10. The method for modeling a three-dimensional geological voxel model based on borehole data according to claim 1, characterized in that, The method for establishing a one-to-one mapping relationship between a three-dimensional geological voxel model and a structured database using a voxel three-dimensional integer index as the unique associative primary key includes: Using the voxel three-dimensional integer index as the globally unique association key, a field-level one-to-one mapping is established between each voxel unit and the geological strata borehole data in the structured database; A hierarchical storage architecture is adopted. The three-dimensional geological voxel model only stores the three-dimensional integer index of the voxel and the stratigraphic information in a lightweight manner, while all multi-dimensional geological attribute parameters are stored independently in the structured database as an external component.

Images