A tetrahedron classification method and system based on a geological body surface grid

Abstract

The application provides a tetrahedron classification method and system based on a geological body surface grid, and relates to the technical field of geological modeling. The method comprises the following steps: acquiring closed surface grid data of a plurality of geological bodies; performing tetrahedron subdivision on the closed surface grid to obtain corresponding tetrahedron grids; performing category division on the subdivided tetrahedrons based on a zoning strategy and a priority strategy according to spatial geometric information and attribute labels of the geological bodies to obtain tetrahedron classification information; identifying contact surfaces in the generated tetrahedron grids, and generating a contact surface category matrix based on the boundary relationship between the plurality of geological bodies; and outputting a tetrahedron grid data structure containing the tetrahedron classification information and the contact surface category matrix. The application can significantly reduce manual operation, improve classification accuracy and simulation efficiency, and is suitable for complex geological numerical simulation tasks.

Description

A tetrahedral classification method and system based on geological body surface mesh Technical Field

[0001] This invention relates to the field of geological modeling technology, and in particular to a tetrahedral classification method and system based on a geological body surface grid. Background Technology

[0002] In geological modeling and numerical simulation, the closed surface meshes of multiple geological bodies usually need to be converted into tetrahedral meshes for subsequent simulation calculations.

[0003] However, when dealing with complex geological formations and multi-contact surface problems, related segmentation techniques often require significant manual intervention, which increases workload and the likelihood of errors. Furthermore, the increasing scale of simulations further exacerbates the technical bottlenecks for multi-classification requirements, making it impossible for related technologies to achieve efficient and accurate automatic classification. Summary of the Invention

[0004] This invention provides a tetrahedral classification method and system based on geological body surface grids, aiming to solve the problems existing in the background art.

[0005] To solve the above-mentioned technical problems, the present invention is implemented as follows:

[0006] In a first aspect, embodiments of the present invention provide a tetrahedral classification method based on a geological body surface grid, the method comprising:

[0007] Acquire closed surface mesh data of multiple geological bodies, wherein the closed surface mesh data includes spatial geometric information and attribute labels of different geological bodies;

[0008] The closed surface mesh is tetrahedralized to obtain the corresponding tetrahedral mesh;

[0009] Based on the spatial geometric information and attribute labels of geological bodies, the segmented tetrahedrons are classified according to partitioning and priority strategies to obtain tetrahedron classification information. The classification ensures that the tetrahedron category is accurately identified and allows the tetrahedron to have multiple categories.

[0010] The contact surfaces in the generated tetrahedral mesh are identified, and a contact surface category matrix is ​​generated based on the boundary relationship between multiple geological bodies. The contact surface category matrix is ​​used to classify the contact surface categories between different geological bodies to ensure accurate expression of contact surface information under multiple geological body conditions.

[0011] The output is a tetrahedral mesh data structure containing the tetrahedral classification information and the contact surface category matrix. The tetrahedral mesh data structure meets the input requirements of a specific numerical simulation software, including FLAC3D and 3DEC.

[0012] Optionally, acquiring the closed surface mesh data of each of the multiple geological bodies includes:

[0013] Obtain the three-dimensional surface information and geological attribute data of each of the multiple geological bodies;

[0014] The three-dimensional surface information of each geological body is meshed, and a corresponding spatial geometric model is constructed.

[0015] Based on the geological attribute data, each geological body is assigned a unique attribute label, which is used to distinguish different geological body categories;

[0016] Based on the spatial geometric model and attribute labels of each geological body, corresponding closed surface mesh data is generated.

[0017] Optionally, the step of tetrahedralizing the closed surface mesh to obtain a corresponding tetrahedral mesh includes:

[0018] The closed surface mesh data is tetrahedralized using a tetrahedral subdivision algorithm to obtain the corresponding tetrahedral mesh. The tetrahedral subdivision algorithm generates a tetrahedral mesh suitable for geological simulation needs based on the boundary and structure of the geological body.

[0019] Adjust the shape and size of the tetrahedral mesh generated during the meshing process, wherein the difference between the ratio of the longest side and the shortest side of each tetrahedron and 1 is less than a preset first threshold, and the volume of each tetrahedron is not less than a preset second threshold.

[0020] The control parameters of the tetrahedral subdivision algorithm are adjusted to adapt to the shapes and boundaries of different geological bodies and to support the automatic subdivision requirements of multiple types of geological body scenarios.

[0021] Optionally, the step of classifying the segmented tetrahedrons based on the spatial geometric information and attribute labels of the geological body, using a partitioning strategy and a priority strategy, to obtain tetrahedron classification information includes:

[0022] Based on the spatial geometric information and attribute labels of the geological bodies, each tetrahedron is classified into the corresponding geological body category according to its spatial location.

[0023] Based on the preset priority of geological bodies, tetrahedrons belonging to multiple geological body categories are classified into the geological body category with higher priority.

[0024] Based on the geological body category to which the tetrahedron belongs, the corresponding tetrahedron classification information is obtained.

[0025] Optionally, the step of identifying contact surfaces in the generated tetrahedral mesh and generating a contact surface category matrix based on the boundary relationships between multiple geological bodies includes:

[0026] Identify the contact surfaces in a tetrahedral mesh at the interface of multiple geological bodies;

[0027] Based on the spatial geometric information and attribute labels of the geological body, each contact surface is classified into the corresponding contact surface category according to its spatial location;

[0028] Based on the contact surface categories of the contact surfaces in the tetrahedral mesh, a corresponding contact surface category matrix is ​​constructed. The contact surface category matrix defines the contact relationship between each pair of geological bodies and the contact surface category to which each contact surface belongs. The contact surface category to which each contact surface belongs is stored in the contact surface category matrix in matrix form.

[0029] Optionally, the method further includes:

[0030] The at least one geological body category to which each tetrahedron belongs, and the contact surface category to which each contact surface belongs, are output as multi-classification information to the tetrahedral mesh data structure. The multi-classification information is used to identify the boundary relationships between different geological bodies and the geological body category to which each tetrahedron belongs during numerical simulation.

[0031] Optionally, the output includes a tetrahedral mesh data structure comprising the tetrahedral classification information and the contact surface category matrix, including:

[0032] The output is a tetrahedral mesh data structure containing the tetrahedral mesh, the tetrahedral classification information of the tetrahedral mesh, and the contact surface category matrix, in a file format that conforms to the input format of the specific numerical simulation software. The specific numerical simulation software directly calls the tetrahedral mesh data in the simulation analysis.

[0033] Optionally, the classification of the partitioned tetrahedrons based on partitioning and priority strategies includes:

[0034] Based on partitioning and priority strategies, OpenCL technology is used to perform parallel computation on the class division of the tetrahedrons on the GPU to improve the efficiency of class division of the tetrahedrons.

[0035] Optionally, the step of using a tetrahedral subdivision algorithm to tetrahedronize the closed surface mesh data to obtain the corresponding tetrahedral mesh includes:

[0036] Load closed surface mesh data of multiple geological bodies, and repair duplicate vertices, duplicate faces, zero-area faces and holes in the closed surface mesh data;

[0037] The BW algorithm is used to calculate the convex hull of all input points in the repaired closed surface mesh data;

[0038] The boundary surfaces of the original closed surface mesh data are inserted into the convex hull to obtain the corresponding tetrahedral mesh, wherein the tetrahedral mesh contains the geometric information of the original closed surface mesh data.

[0039] Secondly, embodiments of the present invention provide a tetrahedral classification system based on a geological body surface grid, the system comprising:

[0040] The input module is used to acquire closed surface mesh data of multiple geological bodies, wherein the closed surface mesh data includes spatial geometric information and attribute labels of different geological bodies;

[0041] The partitioning module is used to perform tetrahedral partitioning on the closed surface mesh to obtain a corresponding tetrahedral mesh. The tetrahedral mesh includes multiple tetrahedra, and the difference between the ratio of the longest side and the shortest side of each tetrahedral and 1 is less than a preset first threshold, and the volume of each tetrahedral is not less than a preset second threshold.

[0042] The category division module is used to classify the segmented tetrahedrons according to the spatial geometric information and attribute labels of the geological body, based on the partitioning strategy and the priority strategy, to obtain tetrahedron classification information. The category division ensures that the tetrahedron category is accurately identified and allows the tetrahedron to have multiple categories.

[0043] The contact surface recognition module is used to identify the contact surfaces in the generated tetrahedral mesh and generate a contact surface category matrix based on the boundary relationship between multiple geological bodies. The contact surface category matrix is ​​used to classify the contact surface categories between different geological bodies to ensure accurate expression of contact surface information under multiple geological body conditions.

[0044] The output module is used to output a tetrahedral mesh data structure containing the tetrahedral classification information and the contact surface category matrix. The tetrahedral mesh data structure meets the input requirements of a specific numerical simulation software, including FLAC3D and 3DEC.

[0045] Optionally, the input module includes:

[0046] The acquisition unit is used to acquire the three-dimensional surface information and geological attribute data of each of the multiple geological bodies;

[0047] Mesh processing units are used to mesh the three-dimensional surface information of each geological body and construct the corresponding spatial geometric model;

[0048] A tagging unit is used to assign a unique attribute tag to each geological body based on the geological attribute data, and the attribute tag is used to distinguish different geological body categories;

[0049] The generation unit is used to generate corresponding closed surface mesh data based on the spatial geometric model and attribute labels of each geological body.

[0050] Optionally, the segmentation module includes:

[0051] The subdivision unit is used to perform tetrahedral subdivision on the closed surface mesh data using a tetrahedral subdivision algorithm to obtain the corresponding tetrahedral mesh. The tetrahedral subdivision algorithm generates a tetrahedral mesh suitable for geological simulation needs based on the boundary and structure of the geological body.

[0052] The first adjustment unit is used to adjust the shape and size of the tetrahedral mesh generated during the meshing process, wherein the difference between the ratio of the longest side and the shortest side of each tetrahedron and 1 is less than a preset first threshold, and the volume of each tetrahedron is not less than a preset second threshold.

[0053] The second adjustment unit is used to adjust the control parameters of the tetrahedral subdivision algorithm to adapt to the shape and boundary of different geological bodies and support the automatic subdivision requirements of multiple types of geological body scenarios.

[0054] Optionally, the category classification module includes:

[0055] The first classification unit is used to classify each tetrahedron into the corresponding geological body category according to its spatial geometric information and attribute labels.

[0056] The second classification unit is used to classify tetrahedrons belonging to multiple geological body categories into the category with higher priority based on the preset priority of geological bodies.

[0057] The information generation unit is used to obtain the corresponding tetrahedral classification information based on the geological body category to which the tetrahedron belongs.

[0058] Optionally, the contact surface recognition module includes:

[0059] The identification unit is used to identify the contact surfaces in a tetrahedral mesh at the interface of multiple geological bodies.

[0060] The third classification unit is used to classify each contact surface into the corresponding contact surface category according to its spatial location, based on the spatial geometric information and attribute labels of the geological body.

[0061] A construction unit is used to construct a contact surface category matrix based on the contact surface categories of the contact surfaces in the tetrahedral mesh. The contact surface category matrix defines the contact relationship between each pair of geological bodies and the contact surface category to which each contact surface belongs. The contact surface category to which each contact surface belongs is stored in the contact surface category matrix in matrix form.

[0062] Optionally, the system further includes:

[0063] The category information output module is used to output the at least one geological body category to which each tetrahedron belongs, and the contact surface category to which each contact surface belongs, as multi-classification information to the tetrahedral mesh data structure. The multi-classification information is used to identify the boundary relationships between different geological bodies and the geological body category to which each tetrahedron belongs during numerical simulation.

[0064] Optionally, the output module includes:

[0065] The output unit is used to output a tetrahedral mesh data structure containing the tetrahedral mesh, the tetrahedral classification information of the tetrahedral mesh, and the contact surface category matrix in a file format that conforms to the input format of the specific numerical simulation software. The specific numerical simulation software directly calls the tetrahedral mesh data in the simulation analysis.

[0066] Optionally, the category classification module includes:

[0067] The parallel computing unit is used to perform parallel computing on the GPU for the class division of the tetrahedron based on partitioning and priority strategies and using OpenCL technology, so as to improve the class division efficiency of the tetrahedron.

[0068] Optionally, the partitioning unit includes:

[0069] The loading sub-unit is used to load closed surface mesh data of multiple geological bodies and repair duplicate vertices, duplicate faces, zero-area faces and holes in the closed surface mesh data;

[0070] The convex hull calculation sub-element is used to calculate the convex hull of all input points of the repaired closed surface mesh data using the BW algorithm.

[0071] An insertion sub-unit is used to insert the boundary surface of the original closed surface mesh data into the convex hull to obtain a corresponding tetrahedral mesh, wherein the tetrahedral mesh contains the geometric information of the original closed surface mesh data.

[0072] The technical solutions provided by the embodiments of the present invention bring at least the following beneficial effects:

[0073] This invention employs a partitioning and priority strategy to automatically classify tetrahedral meshes, enabling each tetrahedron to be accurately categorized based on its spatial location and the attribute labels of the geological body it belongs to. This significantly reduces the complexity of manual operations and flexibly handles situations involving the boundaries and overlaps of multiple geological bodies, achieving efficient and accurate category identification and meeting the needs for multi-category classification in complex geological scenarios. Furthermore, when dealing with contact surfaces between multiple geological bodies, this invention can automatically identify contact surfaces and generate a contact surface category matrix, accurately representing the boundary relationships between different geological bodies. It supports multi-classification of contact surfaces between complex geological bodies, thus ensuring the accurate representation of contact surface information during numerical simulation. Through automated tetrahedral mesh generation and category division, this invention reduces manual intervention and operational errors, significantly improving the efficiency and accuracy of geological modeling and numerical simulation processes. Especially when facing large-scale, complex scenarios with multiple geological bodies and multiple contact surfaces, it greatly reduces technical bottlenecks, providing more efficient and accurate technical support for numerical simulation in fields such as geological resource exploration and underground engineering design. Attached Figure Description

[0074] Various other advantages and benefits will become apparent to those skilled in the art upon reading the following detailed description of preferred embodiments. The accompanying drawings are for illustrative purposes only and are not intended to limit the invention. Furthermore, the same reference numerals denote the same parts throughout the drawings. In the drawings:

[0075] Figure 1 is a schematic diagram of the steps of a tetrahedral classification method based on a geological body surface grid provided in an embodiment of the present invention;

[0076] Figure 2 is a flowchart of the overall process of a tetrahedral classification method based on geological body surface mesh provided in an embodiment of the present invention;

[0077] Figure 3 is a structural block diagram of a tetrahedral classification system based on a geological body surface grid provided in an embodiment of the present invention. Detailed Implementation

[0078] To make the above-mentioned objects, features, and advantages of the present invention more apparent and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present invention, and are only some, not all, embodiments of the present invention, and are not intended to limit the present invention.

[0079] In geological modeling and numerical simulation, especially when dealing with the intersections of multiple geological bodies and complex contact surfaces, it is often necessary to convert the closed surface mesh of the geological body into a tetrahedral mesh for subsequent simulation calculations. However, traditional mesh generation techniques often require significant manual intervention when facing complex multi-geological-body scenarios, increasing workload and the risk of errors. Particularly when dealing with overlapping areas and interfaces of multiple geological bodies, related techniques cannot achieve efficient and accurate automatic classification and contact surface identification, resulting in low simulation accuracy and efficiency. This invention aims to achieve automated tetrahedral mesh classification, automatic identification of contact surfaces and generation of contact surface category matrices, and, combined with an optimized tetrahedral mesh generation algorithm, achieve efficient tetrahedral mesh generation and classification. This significantly reduces manual intervention, improves the ability to handle complex geological body scenarios, and substantially enhances simulation accuracy and efficiency.

[0080] Figure 1 is a schematic diagram of the steps of a tetrahedral classification method based on a geological body surface grid provided by an embodiment of the present invention. Referring to Figure 1, the method includes the following steps:

[0081] Step 101: Obtain closed surface mesh data of multiple geological bodies, wherein the closed surface mesh data includes spatial geometric information and attribute labels of different geological bodies.

[0082] In embodiments of this invention, closed surface mesh data for each independent geological body is first extracted from a multi-geological-body scene. These multiple geological bodies are independent structures, either underground or on the surface, such as rock strata, ore bodies, and groundwater bodies. These independent geological bodies have specific spatial distributions and geometric shapes. The closed surface mesh is a mesh structure used to represent the outer surface of the geological body and can be obtained through 3D scanning, geological exploration, or geological model generation software. The closed surface mesh data for each geological body represents its outer boundary, ensuring accurate representation of the shape of each geological body during subsequent tetrahedral meshing. The closed surface mesh data for each geological body includes spatial geometric information and its unique attribute labels, used to identify and characterize different geological bodies. The spatial geometric information refers to the geometric features of each geological body in three-dimensional space, such as its position, shape, size, and boundaries. Spatial geometric information can be expressed in the form of point cloud data, triangular meshes, or polyhedral meshes, accurately describing the shape and boundaries of the geological body. The attribute labels refer to the classification and identification of different geological bodies, used to distinguish different types of geological bodies or different geological structures. For example, in mineral exploration, it is necessary to identify different types of rock strata, ore bodies, or water layers. In this invention, each geological body is assigned a unique attribute label so that, during subsequent tetrahedral mesh generation and classification, each tetrahedron can be accurately assigned to its corresponding geological body category based on the label information. Attribute labels help distinguish between different geological body categories and provide a basis for subsequent classification work.

[0083] Step 102: Perform tetrahedral subdivision on the closed surface mesh to obtain the corresponding tetrahedral mesh.

[0084] In embodiments of the present invention, dividing the closed surface mesh into tetrahedral meshes is to transform a complex geometry in three-dimensional space into a data structure suitable for numerical simulation, enabling subsequent simulation calculations to be performed on an accurate three-dimensional discrete model. The present invention utilizes an optimized tetrahedral meshing algorithm to tetrahedronize the closed surface mesh of a geological body, generating a high-quality tetrahedral mesh suitable for numerical simulation analysis. In an optional embodiment, step 102 includes substeps S1021–S1023:

[0085] Step S1021: The closed surface mesh data is tetrahedralized using a tetrahedral subdivision algorithm to obtain the corresponding tetrahedral mesh. The tetrahedral subdivision algorithm generates a tetrahedral mesh suitable for geological simulation requirements based on the boundary and structure of the geological body.

[0086] Given that geological body boundaries are typically complex and irregular, tetrahedral meshing algorithms need to generate adaptive tetrahedral shapes based on these boundaries, ensuring a balanced mesh shape and avoiding elongated, distorted, or anomalously angular tetrahedra. For example, if the geological body's boundaries are very complex or curved, the meshing algorithm needs to ensure that the generated tetrahedral mesh closely fits these boundaries, avoiding large geometric deviations or mismatches. Furthermore, different geological bodies may have different structural characteristics, such as irregular boundaries, abrupt rock strata, or multiple alternating geological structures. Therefore, the tetrahedral meshing algorithm used in this invention not only needs to handle complex boundaries but also needs to address the different structural features between geological bodies. Ultimately, the generated tetrahedral mesh exhibits balanced tetrahedral size and shape, suitable for geological simulation requirements.

[0087] Step S1022: Adjust the shape and size of the tetrahedral mesh generated during the meshing process, wherein the difference between the ratio of the longest side and the shortest side of each tetrahedron and 1 is less than a preset first threshold, and the volume of each tetrahedron is not less than a preset second threshold.

[0088] Because excessively long sides can make tetrahedrons overly elongated, thus affecting the numerical stability of the mesh, to achieve the "mesh shape balance" described in step S1021—that is, the side lengths of tetrahedrons should be as close as possible to avoid excessively long or short sides—this step stipulates that the difference between the ratio of the longest and shortest sides of each tetrahedron and 1 is less than a preset first threshold. Specifically, the ratio of the longest to the shortest side should be close to 1, meaning the shape of the tetrahedron should be as close as possible to a regular tetrahedron. If the ratio difference is too large, it indicates that the tetrahedron shape is not uniform enough, which will affect the calculation accuracy in numerical simulation. Similarly, to ensure that there are no "holes" or excessively large gaps in the mesh, the volume of each tetrahedron must meet a minimum value, that is, not less than a preset second threshold. Tetrahedrons with too small a volume may lead to excessive refinement in local areas, while tetrahedrons with too large a volume may not be able to effectively capture subtle changes in geological bodies. By setting a minimum volume threshold, it is possible to ensure that the tetrahedrons are of reasonable size throughout the mesh, avoiding elements that are too large or too small.

[0089] Step S1023: Adjust the control parameters of the tetrahedral subdivision algorithm to adapt to the shape and boundary of different geological bodies and support the automatic subdivision requirements of multiple types of geological body scenarios.

[0090] Adjusting the control parameters of the tetrahedral meshing algorithm allows the meshing results to adapt to the boundary characteristics of different geological bodies, meeting the requirements of numerical simulation for mesh shape, thereby enhancing computational stability and result accuracy. Based on the geometric features and physical properties of geological bodies, the tetrahedral meshing algorithm can automatically identify areas requiring refinement (such as fractures, faults, veins, seismic zones, etc.) and generate finer tetrahedral meshes in these areas; in simple, uniform areas, it can generate coarser meshes to improve computational efficiency. For complex boundaries, the tetrahedral meshing algorithm can automatically adjust the mesh generation strategy, such as using more nodes or specific boundary alignment algorithms, to make the mesh more accurately fit the boundary of the geological body. Possible implementation methods include using boundary stretching techniques or adaptive interpolation methods to optimize the mesh. The tetrahedral meshing algorithm can preset different meshing strategies and control parameters according to different types of geological bodies. Through flexible parameter adjustment, the algorithm can quickly switch and adapt to the automatic meshing needs of various geological scenarios.

[0091] Step S103: Based on the spatial geometric information and attribute labels of the geological body, the segmented tetrahedrons are classified according to the partitioning strategy and the priority strategy to obtain tetrahedron classification information. The classification ensures that the tetrahedron category is accurately identified and allows the tetrahedron to have multiple categories.

[0092] Based on the generated tetrahedral mesh, the tetrahedra are automatically classified using the spatial geometric information and attribute labels of the geological body. This tetrahedral classification employs a partitioning and priority strategy to ensure efficiency and accuracy. Specifically, for each tetrahedron, based on its location (coordinates in space), combined with the geological characteristics of the region (e.g., which rock layer the region belongs to, fracture characteristics, etc.), and the tetrahedron's attribute labels (e.g., lithology, porosity, etc.), the algorithm determines its category and automatically assigns each tetrahedron to the appropriate category. The entire geological body is divided into multiple sub-regions using a partitioning strategy. These regions can be divided based on different criteria, such as spatial location (e.g., upper and lower rock layers). Each sub-region has unique geological characteristics. For each tetrahedron, its geological region is examined, and its category is assigned based on its attribute labels. For example, if a tetrahedron is located in an aquifer region, it is labeled as the "aquifer" category. In some cases, a tetrahedron may lie within the intersection of multiple regions, or its attribute labels may contain multiple different features. In this situation, a priority strategy is used to determine which category the tetrahedron should be assigned to. For example, if a tetrahedron is located within the intersection of both the water layer and the rock layer, and the rock layer is pre-defined as having a higher priority, then the "rock layer" category is given priority.

[0093] It should be noted that in actual geological bodies, especially in the presence of complex interfaces, faults, and alternating mineral layers, a tetrahedron may involve multiple geological features. Therefore, embodiments of the present invention allow tetrahedra to have multiple categories. For example, a tetrahedron may be located in both "aquifer" and "rock" regions simultaneously. The category information of the tetrahedron as "aquifer" and "rock" should be recorded so that the characteristics of the tetrahedron can be accurately described during subsequent simulations.

[0094] Step S104: Identify the contact surfaces in the generated tetrahedral mesh and generate a contact surface category matrix based on the boundary relationship between multiple geological bodies. The contact surface category matrix is ​​used to classify the contact surface categories between different geological bodies to ensure accurate expression of contact surface information under multiple geological body conditions.

[0095] In this embodiment of the invention, the contact surface category matrix is ​​a data structure used to store and represent the boundary relationships between multiple geological bodies. In the contact surface category matrix, each element corresponds to a contact surface and is labeled with the geological body category involved (for example, a contact surface may involve two geological body categories: "water layer" and "rock layer"), ensuring that each contact surface can accurately express its contact information in numerical simulation. The contact surface category matrix, as contact surface category information, is used in subsequent analysis and simulation to accurately distinguish the boundary characteristics between different geological bodies.

[0096] Step S105: Output a tetrahedral mesh data structure containing the tetrahedral classification information and the contact surface category matrix. The tetrahedral mesh data structure meets the input requirements of a specific numerical simulation software, including FLAC3D and 3DEC.

[0097] In this embodiment of the invention, the tetrahedral mesh data structure with tetrahedral classification information and contact surface category information is exported as a data structure usable by specific numerical simulation software, which can be directly applied to simulation analysis.

[0098] In one optional implementation, the output includes a tetrahedral mesh data structure containing the tetrahedral classification information and the contact surface category matrix, comprising: outputting a tetrahedral mesh data structure containing the tetrahedral mesh, the tetrahedral classification information of the tetrahedral mesh, and the contact surface category matrix in a file format conforming to the input format of the specific numerical simulation software, wherein the specific numerical simulation software directly calls the tetrahedral mesh data in the simulation analysis.

[0099] The output data includes: spatial geometric information of the tetrahedrons, tetrahedron category information, and contact surface category information (presented as a contact surface category matrix). The output tetrahedral mesh data structure needs to meet the input format requirements of specific numerical simulation software, including FLAC3D (Fast Lagrangian Analysis of Continua in 3 Dimensions) and 3DEC (Three-Dimensional Distinct Element Code). FLAC3D and 3DEC are simulation tools for geological engineering and structural mechanics, used in numerical simulations of underground engineering, geotechnical engineering, and geological bodies. The input formats and requirements of each software may differ, therefore the data structure needs to be specifically designed and converted for these software programs. FLAC3D is a computational software for three-dimensional numerical simulation, suitable for geotechnical engineering, structural mechanics, fluid dynamics, and other fields. FLAC3D input includes mesh data, nodes, elements, material properties, boundary conditions, and loads. To ensure the output data meets FLAC3D's input requirements, the tetrahedral mesh needs to be represented as nodes, elements (tetrahedrons), contact surface information, etc., and each element needs to be associated with material and mechanical property information. 3DEC is another numerical software for simulating geotechnical and structural systems, particularly for simulating faults, cracks, and particle systems. 3DEC input files typically include the model's geometry, contact surface characteristics, material models, and physical properties. In 3DEC, the interaction between the tetrahedral mesh and contact surfaces can be represented as different types of mechanical elements, such as contact surfaces and particles.

[0100] When outputting data, the spatial geometry, tetrahedral classification, and contact surface category information need to be converted into a format that the specific numerical simulation software can understand, according to the software's specific input format. The generated tetrahedral mesh data structure can be directly applied to simulation analysis, ensuring data format compatibility and reducing information loss during data processing. In other words, during numerical simulation, the specific numerical simulation software will directly load and use the output tetrahedral mesh data. Based on the structure and classification information of the tetrahedral mesh, the specific numerical simulation software can assign appropriate physical properties for further calculation and analysis.

[0101] This invention employs a partitioning and priority strategy to automatically classify tetrahedral meshes, enabling each tetrahedron to be accurately categorized based on its spatial location and the attribute labels of the geological body it belongs to. This significantly reduces the complexity of manual operations and flexibly handles situations involving the boundaries and overlaps of multiple geological bodies, achieving efficient and accurate category identification and meeting the needs for multi-category classification in complex geological scenarios. Furthermore, when dealing with contact surfaces between multiple geological bodies, this invention can automatically identify contact surfaces and generate a contact surface category matrix, accurately representing the boundary relationships between different geological bodies. It supports multi-classification of contact surfaces between complex geological bodies, thus ensuring the accurate representation of contact surface information during numerical simulation. Through automated tetrahedral mesh generation and category division, this invention reduces manual intervention and operational errors, significantly improving the efficiency and accuracy of geological modeling and numerical simulation processes. Especially when facing large-scale, complex scenarios with multiple geological bodies and multiple contact surfaces, it greatly reduces technical bottlenecks, providing more efficient and accurate technical support for numerical simulation in fields such as geological resource exploration and underground engineering design.

[0102] In one optional implementation, acquiring the closed surface mesh data of each of the multiple geological bodies includes:

[0103] Step S1011: Obtain the three-dimensional surface information and geological attribute data of each of the multiple geological bodies.

[0104] The three-dimensional surface information refers to the geometric shape of each geological body, including its boundaries and external contours in three-dimensional space, which can be obtained through geological exploration data (such as seismic surveys, geological drilling, and remote sensing data). The geological attribute data refers to the physical, chemical, or structural characteristics associated with each geological body. For example, the mineral composition, porosity, density, elastic modulus, and permeability of rock strata can be obtained through laboratory analysis or field measurements.

[0105] Step S1012: The three-dimensional surface information of each geological body is meshed and the corresponding spatial geometric model is constructed.

[0106] Meshing is the process of converting three-dimensional surface information into a mesh structure. The purpose of meshing is to divide the surface of a geological body into smaller tetrahedral units. Meshing makes complex three-dimensional surfaces easier to process and represent in a computer.

[0107] Next, a spatial geometric model is constructed. A spatial geometric model refers to a complete three-dimensional geometric representation formed based on gridded three-dimensional surface information, including elements such as points (vertices), edges, faces, and volumes, defining the precise shape of each geological body.

[0108] Step S1013: Based on the geological attribute data, assign a unique attribute label to each geological body. The attribute label is used to distinguish different geological body categories.

[0109] Based on the geological attribute data, a unique attribute label is assigned to each geological body. Each geological body's attribute label reflects its main characteristics. As mentioned earlier, attribute labels ensure the distinction between different geological bodies and allow the behavior of each geological body to be tracked within the geological model.

[0110] Step S1014: Generate corresponding closed surface mesh data based on the spatial geometric model and attribute labels of each geological body.

[0111] The three-dimensional surface information and attribute labels of each geological body are combined to generate the final closed surface mesh data. Here, "closed" means that the surface of each geological body must be continuous and uninterrupted, that is, there are no holes or cracks on the surface.

[0112] In one optional implementation, the step of classifying the segmented tetrahedrons based on the spatial geometric information and attribute labels of the geological body, using a partitioning strategy and a priority strategy, to obtain tetrahedron classification information includes:

[0113] Step S1031: Based on the spatial geometric information and attribute labels of the geological bodies, each tetrahedron is classified into the corresponding geological body category according to its spatial location.

[0114] Each tetrahedron is automatically assigned to a corresponding geological body category based on its spatial location. If a tetrahedron is entirely located within a single geological body, it will be assigned to that category; if it lies within the intersection of multiple geological bodies, it will be assigned to multiple categories. In this way, each tetrahedron is assigned to a specific geological body category after this step.

[0115] Step S1032: Based on the preset priority of geological bodies, tetrahedrons belonging to multiple geological body categories are classified into the geological body category with higher priority.

[0116] In practical applications, certain geological bodies have higher priority than others. Priority is predefined for each geological body and can be determined based on its geological importance, engineering requirements, or other physical characteristics. For example, if a shale layer and a sandstone layer alternate, and the shale layer has a higher priority (because it has a greater impact on groundwater flow), then the tetrahedron within the intersection should be preferentially classified as a shale layer.

[0117] Step S1033: Obtain the corresponding tetrahedron classification information based on the geological body category to which the tetrahedron belongs.

[0118] Tetrahedral classification information includes the attribute label of the geological body to which the tetrahedron belongs, as well as the final category after the application of a priority strategy. Tetrahedral classification information identifies the geological body category to which each tetrahedron belongs, for use in subsequent analysis, simulation, and visualization.

[0119] The above steps ensured efficient and accurate tetrahedral classification in complex geological models, providing an accurate foundation for subsequent numerical simulations.

[0120] In one optional implementation, the step of identifying contact surfaces in the generated tetrahedral mesh and generating a contact surface category matrix based on the boundary relationships between multiple geological bodies includes:

[0121] Step S1041: Identify the contact surfaces in the tetrahedral mesh at the interfaces of multiple geological bodies.

[0122] In a tetrahedral mesh, contact surfaces can be identified by the shared faces of adjacent tetrahedra in the mesh. A contact surface is defined as a face shared by two tetrahedra that belong to different geological bodies. The identification process can include examining each pair of adjacent tetrahedra in the mesh to determine if they belong to different geological bodies; if so, their shared face is marked as a contact surface.

[0123] Step S1042: Based on the spatial geometric information and attribute labels of the geological body, each contact surface is classified into the corresponding contact surface category according to its spatial location.

[0124] By using the spatial geometry of the tetrahedron and the spatial location of the geological bodies, the positional relationship of each contact surface can be determined, and it can be identified which two geological bodies each contact surface belongs to. Furthermore, the attribute tags of the geological bodies can further confirm which two geological bodies the contact surface belongs to. For example, if the contact surface is located at the boundary between sandstone and shale, then this contact surface will be labeled as a "sandstone-shale contact surface." Based on the location of the contact surface and the category of the geological body it belongs to, a unique contact surface category can be assigned to each contact surface.

[0125] Step S1043: Based on the contact surface categories of the contact surfaces in the tetrahedral mesh, a corresponding contact surface category matrix is ​​constructed. The contact surface category matrix defines the contact relationship between each pair of geological bodies and the contact surface category to which each contact surface belongs. The contact surface category to which each contact surface belongs is stored in the contact surface category matrix in matrix form.

[0126] Based on the identified contact surface categories, a contact surface category matrix is ​​constructed. This matrix is ​​a two-dimensional symmetric structure, where the row and column indices correspond to different geological body types, and the matrix elements represent different contact surface categories. The contact surface category matrix is ​​stored using a compact data structure (such as a symmetric matrix or dictionary mapping), supporting fast lookup of the contact surface category between any two geological bodies. For example, when sandstone (ID:1) and shale (ID:2) have a boundary contact surface, the (1,2) and (2,1) positions in the matrix both store category codes 1-2.

[0127] The above steps ensured that the contact surfaces of geological bodies were efficiently and accurately classified in the complex geological body model, providing an accurate foundation for subsequent numerical simulations.

[0128] In an optional implementation, the method further includes: outputting at least one geological body category to which each tetrahedron belongs, and the contact surface category to which each contact surface belongs, as multi-classification information to the tetrahedral mesh data structure. The multi-classification information is used to identify the boundary relationships between different geological bodies and the geological body category to which each tetrahedron belongs during numerical simulation.

[0129] Multi-classification information refers to classifying each tetrahedron and contact surface into its corresponding category so that the relationships between them can be identified during subsequent simulations. That is, each tetrahedron and contact surface is simultaneously associated with a corresponding category, representing either the geological body category or the contact surface category, respectively.

[0130] A tetrahedral mesh data structure consists of several tetrahedral cells and their contact surfaces. Each tetrahedron has a specific spatial location, defined by the coordinates of its four vertices. The mesh structure can be represented as a complex data structure containing elements such as nodes, edges, and faces; it can be a three-dimensional mesh composed of multiple tetrahedral cells. In a tetrahedral mesh data structure, tetrahedra are stored using nodes (vertices) and cells (tetrahedrons) as basic elements. For each tetrahedron, the mesh data structure records the geological body category to which the tetrahedron belongs and its adjacent contact surfaces.

[0131] The tetrahedral classification information for each tetrahedron is stored in the attribute field of that tetrahedron, and the contact surface category for each contact surface is also marked in the tetrahedral mesh data structure.

[0132] It is understood that the embodiments of the present invention support multi-classification of tetrahedrons and contact surfaces based on multi-classification information. Multi-classification information allows each tetrahedron and contact surface to simultaneously reflect the geological body category to which it belongs. In this embodiment, the numerical models for each tetrahedron and contact surface can more flexibly adapt to complex geological conditions and better reflect the diversity and complexity of geological body boundaries in the real world, thereby improving the accuracy of simulation results.

[0133] In one optional implementation, classifying the partitioned tetrahedrons based on partitioning and priority strategies includes: using OpenCL technology on a GPU to perform parallel computation on the classification of the tetrahedrons based on partitioning and priority strategies, so as to improve the classification efficiency of the tetrahedrons.

[0134] OpenCL (Open Computing Language) is a framework for writing programs that can execute in parallel on various hardware platforms. It supports parallel computing on CPUs, GPUs, and other hardware, and can effectively allocate computing tasks to improve computational efficiency.

[0135] In one embodiment of the present invention, OpenCL technology is used to perform a parallel tetrahedral classification task on a GPU. Compared to a traditional CPU, a GPU has more computing cores, enabling it to handle a large number of computational tasks simultaneously. Especially for tasks like tetrahedral classification, where the classification process is typically independent, it is well-suited for parallel computation. A GPU can process multiple tetrahedral classification tasks in parallel on each computing core, thereby significantly improving computational efficiency. Compared to using multi-threaded parallelism (such as TBB: Threading Building Blocks) on a CPU, GPU acceleration can achieve a much higher speedup. According to actual tests, the computational speed using GPU acceleration can be up to 40 times faster than using TBB for parallel computation on a CPU. This significantly reduces computation time when handling a large number of tetrahedral classification tasks, especially with large-scale grid data, where the advantages of GPUs are even more pronounced.

[0136] In one optional implementation, the step of using a tetrahedral subdivision algorithm to tetrahedronize the closed surface mesh data to obtain a corresponding tetrahedral mesh includes:

[0137] Load closed surface mesh data of multiple geological bodies, and repair duplicate vertices, duplicate faces, zero-area faces, and holes in the closed surface mesh data.

[0138] In this embodiment of the invention, a tetrahedralization process is performed on closed surface meshes based on an optimized tetrahedral partitioning algorithm. Closed surface mesh data often exhibits geometric problems such as duplicate vertices, duplicate faces, zero-area faces, and holes. The steps to repair these geometric problems ensure that the input mesh data is clean and ready for subsequent geometric processing and tetrahedral partitioning. Specifically, duplicate vertices are merged into a single vertex to avoid redundancy and unnecessary computation; duplicate faces and zero-area faces are deleted to simplify the mesh; and the process of repairing holes involves filling in missing portions to ensure the mesh is closed.

[0139] The convex hull of all input points in the repaired closed surface mesh data is calculated using the Bowyer-Watson algorithm.

[0140] The Browsing and Woodruff algorithm (BW) is an algorithm for calculating the convex hull of a 3D point set. The convex hull is the smallest convex polyhedron that encloses a set of points. For a given set of points, the convex hull is a closed shape that contains all the input points and is the smallest convex boundary of those points. Given a set of points scattered in 3D space, the BW algorithm can compute a smallest polyhedron that encloses these points; the surface of this polyhedron is the convex hull of the point set. After fixing problems in closed-surface meshes (such as duplicate vertices, holes, etc.), the BW algorithm is applied to all input points to compute their convex hulls, which serve as the basis for generating tetrahedral meshes.

[0141] The boundary surface of the original closed surface mesh data is inserted into the convex hull to obtain the corresponding tetrahedral mesh, wherein the tetrahedral mesh has the same geometry as the original closed surface mesh data.

[0142] Boundary surfaces refer to the outer surfaces of the original mesh. By inserting these boundary surfaces (i.e., the surfaces of the original mesh) into the convex hull obtained through the BW algorithm, it is ensured that the final tetrahedral mesh not only conforms to the shape of the convex hull but also retains the geometry of the original mesh. Assuming the original closed-surface mesh is a complex rock strata model, its outer surface may consist of multiple faces. When generating the tetrahedral mesh, these outer surfaces need to be preserved and inserted into the calculated convex hull to ensure that the generated tetrahedral mesh is precisely aligned with the shape of the original mesh. The final generated tetrahedral mesh has the same geometry as the original closed-surface mesh. That is, the shape and boundaries of the generated tetrahedral mesh will precisely retain the geometric features of the original mesh (such as boundary shape, size, etc.). However, the internal structure will be divided into tetrahedral elements, which are the basis for subsequent calculations and simulation analysis.

[0143] Figure 2 is a flowchart illustrating the overall process of a tetrahedral classification method based on geological body surface mesh provided in this embodiment of the invention. Steps S101 to S105 can be summarized by the flowchart in Figure 2: Step 1: Obtain closed surface mesh data for multiple geological bodies, and establish a geometric model and attribute labels for each geological body. The result is a preliminary construction of the three-dimensional mesh and attribute information for each geological body, ensuring the clarity of geological body boundaries and categories during subsequent subdivision.

[0144] Step 2 involves tetrahedralizing the closed surface mesh using an optimized tetrahedral meshing algorithm. The result is a high-quality tetrahedral mesh with a balanced shape, suitable for numerical simulation of complex geological structures, and avoids the generation of elongated or distorted tetrahedra.

[0145] Step 3 automatically categorizes the generated tetrahedral meshes. Through partitioning and priority strategies, it ensures that the tetrahedral meshes are accurately classified according to their respective geological body categories. The result is categorized data of tetrahedral meshes with category labels, which can be directly applied in numerical simulation analysis.

[0146] Step 4 automatically identifies the contact surfaces between multiple geological bodies and classifies them based on a contact surface category matrix. The result is a classification matrix containing contact surface category information within a tetrahedral mesh, ensuring accurate representation of the contact relationships between geological bodies in numerical simulations.

[0147] Step 5 outputs a tetrahedral mesh data structure conforming to the input format of simulation software (such as FLAC3D, 3DEC), including tetrahedral geometric information, tetrahedral classification information, and contact surface category matrix. The result is a highly compatible mesh data file that can be directly applied to simulation analysis, reducing the complexity of data conversion.

[0148] In summary, the tetrahedral classification method based on geological body surface mesh of the present invention realizes the automated tetrahedral mesh generation and classification process, and is suitable for numerical simulation needs of complex interface scenarios of multiple geological bodies.

[0149] The technical solution of this invention has good versatility and adaptability, and is particularly suitable for numerical simulation systems that require efficient and accurate processing of contact relationships between multiple geological bodies. It can significantly reduce the involvement of manual operation and improve the accuracy and stability of simulation data.

[0150] Figure 3 is a structural block diagram of a tetrahedral classification system based on a geological body surface grid provided in an embodiment of the present invention. As shown in Figure 3, the system includes:

[0151] Input module 201 is used to acquire closed surface mesh data of multiple geological bodies, wherein the closed surface mesh data includes spatial geometric information and attribute labels of different geological bodies;

[0152] The subdivision module 202 is used to perform tetrahedral subdivision on the closed surface mesh to obtain a corresponding tetrahedral mesh. The tetrahedral mesh includes multiple tetrahedra. The difference between the ratio of the longest side and the shortest side of each tetrahedral and 1 is less than a preset first threshold, and the volume of each tetrahedral is not less than a preset second threshold.

[0153] The category division module 203 is used to classify the segmented tetrahedrons according to the spatial geometric information and attribute labels of the geological body, based on the partitioning strategy and the priority strategy, to obtain tetrahedron classification information. The category division ensures that the tetrahedron category is accurately identified and allows the tetrahedron to have multiple categories.

[0154] The contact surface recognition module 204 is used to identify the contact surfaces in the generated tetrahedral mesh and generate a contact surface category matrix based on the boundary relationship between multiple geological bodies. The contact surface category matrix is ​​used to classify the contact surface categories between different geological bodies to ensure accurate expression of contact surface information under multiple geological body conditions.

[0155] Output module 205 is used to output a tetrahedral mesh data structure containing the tetrahedral classification information and the contact surface category matrix. The tetrahedral mesh data structure meets the input requirements of a specific numerical simulation software, including FLAC3D and 3DEC.

[0156] In one optional implementation, the input module includes:

[0157] The acquisition unit is used to acquire the three-dimensional surface information and geological attribute data of each of the multiple geological bodies;

[0158] Mesh processing units are used to mesh the three-dimensional surface information of each geological body and construct the corresponding spatial geometric model;

[0159] A tagging unit is used to assign a unique attribute tag to each geological body based on the geological attribute data, and the attribute tag is used to distinguish different geological body categories;

[0160] The generation unit is used to generate corresponding closed surface mesh data based on the spatial geometric model and attribute labels of each geological body.

[0161] In one optional implementation, the partitioning module includes:

[0162] The subdivision unit is used to perform tetrahedral subdivision on the closed surface mesh data using a tetrahedral subdivision algorithm to obtain the corresponding tetrahedral mesh. The tetrahedral subdivision algorithm generates a tetrahedral mesh suitable for geological simulation needs based on the boundary and structure of the geological body.

[0163] The first adjustment unit is used to adjust the shape and size of the tetrahedral mesh generated during the meshing process, wherein the difference between the ratio of the longest side and the shortest side of each tetrahedron and 1 is less than a preset first threshold, and the volume of each tetrahedron is not less than a preset second threshold.

[0164] The second adjustment unit is used to adjust the control parameters of the tetrahedral subdivision algorithm to adapt to the shape and boundary of different geological bodies and support the automatic subdivision requirements of multiple types of geological body scenarios.

[0165] In one optional implementation, the category classification module includes:

[0166] The first classification unit is used to classify each tetrahedron into the corresponding geological body category according to its spatial geometric information and attribute labels.

[0167] The second classification unit is used to classify tetrahedrons belonging to multiple geological body categories into the category with higher priority based on the preset priority of geological bodies.

[0168] The information generation unit is used to obtain the corresponding tetrahedral classification information based on the geological body category to which the tetrahedron belongs.

[0169] In one optional implementation, the contact surface recognition module includes:

[0170] The identification unit is used to identify the contact surfaces in a tetrahedral mesh at the interface of multiple geological bodies.

[0171] The third classification unit is used to classify each contact surface into the corresponding contact surface category according to its spatial location, based on the spatial geometric information and attribute labels of the geological body.

[0172] A construction unit is used to construct a contact surface category matrix based on the contact surface categories of the contact surfaces in the tetrahedral mesh. The contact surface category matrix defines the contact relationship between each pair of geological bodies and the contact surface category to which each contact surface belongs. The contact surface category to which each contact surface belongs is stored in the contact surface category matrix in matrix form.

[0173] In one alternative implementation, the system further includes:

[0174] The category information output module is used to output the at least one geological body category to which each tetrahedron belongs, and the contact surface category to which each contact surface belongs, as multi-classification information to the tetrahedral mesh data structure. The multi-classification information is used to identify the boundary relationships between different geological bodies and the geological body category to which each tetrahedron belongs during numerical simulation.

[0175] In one optional implementation, the output module includes:

[0176] The output unit is used to output a tetrahedral mesh data structure containing the tetrahedral mesh, the tetrahedral classification information of the tetrahedral mesh, and the contact surface category matrix in a file format that conforms to the input format of the specific numerical simulation software. The specific numerical simulation software directly calls the tetrahedral mesh data in the simulation analysis.

[0177] In one optional implementation, the category classification module includes:

[0178] The parallel computing unit is used to perform parallel computing on the GPU for the class division of the tetrahedron based on partitioning and priority strategies and using OpenCL technology, so as to improve the class division efficiency of the tetrahedron.

[0179] In one optional implementation, the partitioning unit includes:

[0180] The loading sub-unit is used to load closed surface mesh data of multiple geological bodies and repair duplicate vertices, duplicate faces, zero-area faces and holes in the closed surface mesh data;

[0181] The convex hull calculation sub-element is used to calculate the convex hull of all input points of the repaired closed surface mesh data using the BW algorithm.

[0182] An insertion sub-unit is used to insert the boundary surface of the original closed surface mesh data into the convex hull to obtain a corresponding tetrahedral mesh, wherein the tetrahedral mesh has the same geometry as the original closed surface mesh data.

[0183] Although preferred embodiments of the present invention have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of the embodiments of the present invention.

[0184] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or terminal device that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or terminal device. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or terminal device that includes said element.

[0185] The embodiments of the present invention have been described above with reference to the accompanying drawings. However, the present invention is not limited to the specific embodiments described above. The specific embodiments described above are merely illustrative and not restrictive. Those skilled in the art can make many other forms under the guidance of the present invention without departing from the spirit and scope of the claims. All of these forms are within the protection scope of the present invention.

Claims

1. A tetrahedral classification method based on geological body surface mesh, characterized in that, The method includes: acquiring closed surface mesh data of multiple geological bodies, the closed surface mesh data including spatial geometric information and attribute labels of different geological bodies; performing tetrahedral subdivision on the closed surface mesh to obtain corresponding tetrahedral meshes; classifying the subdivided tetrahedrals according to the spatial geometric information and attribute labels of the geological bodies, based on a partitioning strategy and a priority strategy, to obtain tetrahedral classification information, the classification ensuring accurate tetrahedral category identification and allowing the tetrahedrals to have multiple categories; identifying the contact surfaces in the generated tetrahedral meshes, and generating a contact surface category matrix based on the boundary relationships between multiple geological bodies, the contact surface category matrix being used to classify the contact surface categories between different geological bodies, ensuring accurate expression of contact surface information under multi-geological-body conditions; outputting a tetrahedral mesh data structure containing the tetrahedral classification information and the contact surface category matrix, the tetrahedral mesh data structure meeting the input requirements of specific numerical simulation software, the specific numerical simulation software including FLAC3D and 3DEC; wherein, the partitioning strategy includes dividing the entire geological body into multiple sub-regions, the sub-regions being divided based on different spatial locations.

2. The tetrahedral classification method based on geological body surface mesh according to claim 1, characterized in that, The process of acquiring closed surface mesh data for multiple geological bodies includes: acquiring the three-dimensional surface information and geological attribute data of each of the multiple geological bodies; performing meshing processing on the three-dimensional surface information of each geological body and constructing a corresponding spatial geometric model; assigning a unique attribute label to each geological body based on the geological attribute data, the attribute label being used to distinguish different geological body categories; and generating corresponding closed surface mesh data based on the spatial geometric model and attribute label of each geological body.

3. The tetrahedral classification method based on geological body surface grids according to claim 1, characterized in that, The step of tetrahedralizing the closed surface mesh to obtain a corresponding tetrahedral mesh includes: using a tetrahedral subdivision algorithm to tetrahedralize the closed surface mesh data to obtain a corresponding tetrahedral mesh, wherein the tetrahedral subdivision algorithm generates a tetrahedral mesh suitable for geological simulation requirements based on the boundary and structure of the geological body; adjusting the shape and size of the tetrahedral mesh generated during the subdivision process, wherein the difference between the ratio of the longest side and the shortest side of each tetrahedron and 1 is less than a preset first threshold, and the volume of each tetrahedron is not less than a preset second threshold; and adjusting the control parameters of the tetrahedral subdivision algorithm to adapt to the shape and boundary of different geological bodies and support the automatic subdivision requirements of multiple types of geological body scenes.

4. The tetrahedral classification method based on geological body surface grids according to claim 1, characterized in that, The process of classifying the segmented tetrahedrons based on the spatial geometric information and attribute labels of the geological bodies, using a partitioning strategy and a priority strategy, to obtain tetrahedron classification information includes: classifying each tetrahedron into a corresponding geological body category according to its spatial location based on the spatial geometric information and attribute labels of the geological bodies; classifying tetrahedrons belonging to multiple geological body categories into the category with higher priority based on a preset priority of geological bodies; and obtaining the corresponding tetrahedron classification information based on the geological body category to which the tetrahedron belongs.

5. The tetrahedral classification method based on geological body surface mesh according to claim 1, characterized in that, The process of identifying contact surfaces in the generated tetrahedral mesh and generating a contact surface category matrix based on the boundary relationships between multiple geological bodies includes: identifying contact surfaces in the tetrahedral mesh at the interfaces of multiple geological bodies; classifying each contact surface into a corresponding contact surface category according to its spatial location based on the spatial geometric information and attribute labels of the geological bodies; and constructing a contact surface category matrix based on the contact surface categories of the contact surfaces in the tetrahedral mesh. The contact surface category matrix defines the contact relationship between each pair of geological bodies and the contact surface category to which each contact surface belongs, wherein the contact surface category to which each contact surface belongs is stored in the contact surface category matrix in matrix form.

6. The tetrahedral classification method based on geological body surface mesh according to claim 4 or 5, characterized in that, The method further includes: outputting at least one geological body category to which each tetrahedron belongs, and the contact surface category to which each contact surface belongs, as multi-classification information into the tetrahedral mesh data structure. The multi-classification information is used to identify the boundary relationships between different geological bodies and the geological body category to which each tetrahedron belongs during numerical simulation.

7. The tetrahedral classification method based on geological body surface mesh according to claim 1, characterized in that, The output includes a tetrahedral mesh data structure containing the tetrahedral classification information and the contact surface category matrix, comprising: outputting a tetrahedral mesh data structure containing the tetrahedral mesh, the tetrahedral classification information of the tetrahedral mesh, and the contact surface category matrix in a file format conforming to the input format of the specific numerical simulation software, wherein the specific numerical simulation software directly calls the tetrahedral mesh data in the simulation analysis.

8. The tetrahedral classification method based on geological body surface mesh according to claim 1, characterized in that, The method of classifying the tetrahedrons based on partitioning and priority strategies includes: using OpenCL technology on the GPU to perform parallel computation on the classification of the tetrahedrons based on partitioning and priority strategies, so as to improve the classification efficiency of the tetrahedrons.

9. The tetrahedral classification method based on geological body surface grids according to claim 3, characterized in that, The process of using a tetrahedral partitioning algorithm to tetrahedronize the closed surface mesh data to obtain a corresponding tetrahedral mesh includes: loading closed surface mesh data from multiple geological bodies and repairing duplicate vertices, duplicate faces, zero-area faces, and holes in the closed surface mesh data; calculating the convex hull of all input points in the repaired closed surface mesh data using the BW algorithm; and inserting the boundary faces of the original closed surface mesh data into the convex hull to obtain a corresponding tetrahedral mesh, wherein the tetrahedral mesh contains the geometric information of the original closed surface mesh data.

10. A tetrahedral classification system based on a geological body surface grid, characterized in that, The system includes: an input module for acquiring closed surface mesh data of multiple geological bodies, the closed surface mesh data including spatial geometric information and attribute labels of different geological bodies; a partitioning module for tetrahedral partitioning the closed surface mesh to obtain corresponding tetrahedral meshes, wherein each tetrahedral mesh includes multiple tetrahedra, the difference between the ratio of the longest and shortest sides of each tetrahedron and 1 is less than a preset first threshold, and the volume of each tetrahedron is not less than a preset second threshold; and a classification module for classifying the partitioned tetrahedra according to the spatial geometric information and attribute labels of the geological bodies, based on a partitioning strategy and a priority strategy, to obtain tetrahedral classification information, wherein the classification ensures accurate tetrahedral classification and allows... The tetrahedron has multiple categories; the contact surface identification module is used to identify the contact surfaces in the generated tetrahedral mesh and generate a contact surface category matrix based on the boundary relationship between multiple geological bodies. The contact surface category matrix is ​​used to classify the contact surface categories between different geological bodies to ensure accurate expression of contact surface information under multiple geological body conditions; the output module is used to output a tetrahedral mesh data structure containing the tetrahedral classification information and the contact surface category matrix. The tetrahedral mesh data structure meets the input requirements of specific numerical simulation software, including FLAC3D and 3DEC; wherein, the partitioning strategy includes dividing the entire geological body into multiple sub-regions, and the sub-regions are divided based on different spatial locations.

Images