Dynamic visualization method, system and device for three-dimensional geological model, and storage medium

Abstract

Disclosed in the present invention are a dynamic visualization method, system and device for a three-dimensional geological model, and a storage medium. The method comprises: performing voxelization processing on a geological body model, so as to obtain cell position attribute information corresponding to the geological body model; determining whether the cell position attribute information meets a preset model position condition; if the cell position attribute information meets the preset model position condition, obtaining cell value attribute information on the basis of sample point data and by means of a three-dimensional interpolation algorithm, and obtaining a regular hexahedral grid model on the basis of layer number information, a cell index number, the cell position attribute information and the cell value attribute information; performing triangular-face intersection on the geological body model and the regular hexahedral grid model, so as to obtain a corner-point grid model; and performing dynamic visualization on the basis of the regular hexahedral grid model and the corner-point grid model. In the present invention, the regular hexahedral grid model and the corner-point grid model are obtained by means of a voxelization algorithm, the three-dimensional interpolation algorithm, and an intersection algorithm for a geological body and a regular hexahedron, thereby realizing a dynamic visualization function of a three-dimensional geological model.

Description

Three-dimensional geological model dynamic visualization method, system, equipment and storage medium Technical Field

[0001] The present invention relates to the field of geographic information technology, and in particular to a three-dimensional geological model dynamic visualization method, system, equipment and storage medium. Background Art

[0002] Most current 3D geological models use irregular triangulated networks (ITMs) for data organization and management. By constructing the surface of a 3D object, the entire object achieves a three-dimensional spatial effect, making it easier to visualize. 3D geological models constructed based on surface models are pseudo-3D, not true 3D geological models. They cannot express internal heterogeneous geological phenomena and are difficult to handle for local spatial analysis and queries. This is particularly true for coal mine geological model applications, where intelligent mining requires efficient processing such as single-point queries and multi-attribute analysis. Therefore, how to quickly implement dynamic visualization capabilities for regular hexahedral and corner point mesh models has become an urgent issue.

[0003] The above content is only used to assist in understanding the technical solution of the present invention and does not constitute an admission that the above content is prior art.

[0004] Summary of the Invention

[0005] The main purpose of the present invention is to provide a three-dimensional geological model dynamic visualization method, system, equipment and storage medium, aiming to solve the technical problem of how to quickly realize the dynamic visualization function of regular hexahedral grid models and corner point grid models.

[0006] To achieve the above-mentioned object, the present invention provides a three-dimensional geological model dynamic visualization method, which includes:

[0007] voxelizing the geological model to obtain grid position attribute information corresponding to the geological model;

[0008] Determining whether the grid position attribute information meets a preset model position condition;

[0009] If satisfied, grid numerical attribute information is obtained through a three-dimensional interpolation algorithm according to the sample point data, and a regular hexahedral grid model is obtained according to the layer number information, the grid index number, the grid position attribute information and the grid numerical attribute information;

[0010] Performing triangular intersection on the geological body model and the regular hexahedral grid model to obtain a corner point grid model;

[0011] The three-dimensional geological model is dynamically visualized based on the regular hexahedral grid model and the corner point grid model.

[0012] Optionally, before the step of obtaining a regular hexahedral grid model according to the layer number information, the grid index number, the grid position attribute information, and the grid value attribute information, the step further includes:

[0013] Determining layer number information according to the geological body model, and determining grid three-dimensional coordinate information according to the grid position attribute information;

[0014] A grid index number is determined according to the grid three-dimensional coordinate information.

[0015] Optionally, the step of performing triangular intersection of the geological body model and the regular hexahedral grid model to obtain a corner grid model includes:

[0016] Performing triangulation intersection on the geological body model and the regular hexahedral grid model to obtain a triangulation intersection result;

[0017] The triangulated mesh intersection result is processed to obtain a corner point mesh model.

[0018] Optionally, after the step of dynamically visualizing the three-dimensional geological model based on the regular hexahedral grid model and the corner point grid model, the method further includes:

[0019] Determining the grid model data formats of the regular hexahedral grid model and the corner point grid model respectively;

[0020] The regular hexahedral grid model and the corner point grid model are stored according to the grid model data format.

[0021] Optionally, the step of storing the regular hexahedral grid model and the corner grid model according to the grid model data format includes:

[0022] The regular hexahedral grid model and the corner point grid model are respectively divided into blocks by a hierarchical block organization strategy to obtain a divided regular hexahedral grid model and a divided corner point grid model;

[0023] The divided regular hexahedral grid model and the divided corner grid model are stored according to the grid model data format.

[0024] In addition, to achieve the above-mentioned purpose, the present invention further proposes a three-dimensional geological model dynamic visualization system, which includes:

[0025] A processing module, configured to voxelize the geological model to obtain grid position attribute information corresponding to the geological model;

[0026] A judgment module, used to judge whether the grid position attribute information meets the preset model position condition;

[0027] A generation module, configured to obtain grid numerical attribute information based on the sample point data through a three-dimensional interpolation algorithm if the conditions are met, and obtain a regular hexahedral grid model based on the layer number information, the grid index number, the grid position attribute information, and the grid numerical attribute information;

[0028] The generating module is further configured to perform triangular intersection between the geological body model and the regular hexahedral grid model to obtain a corner point grid model;

[0029] A display module is used to dynamically visualize the three-dimensional geological model based on the regular hexahedral grid model and the corner point grid model.

[0030] In addition, to achieve the above-mentioned purpose, the present invention also proposes a three-dimensional geological model dynamic visualization device, which includes: a memory, a processor, and a three-dimensional geological model dynamic visualization program stored in the memory and executable on the processor, wherein the three-dimensional geological model dynamic visualization program is configured to implement the steps of the three-dimensional geological model dynamic visualization method described above.

[0031] In addition, to achieve the above-mentioned purpose, the present invention also proposes a storage medium, on which a three-dimensional geological model dynamic visualization program is stored. When the three-dimensional geological model dynamic visualization program is executed by a processor, the steps of the three-dimensional geological model dynamic visualization method described above are implemented.

[0032] The present invention first voxelizes a geological body model to obtain grid position attribute information corresponding to the geological body model. It then determines whether the grid position attribute information meets a preset model position condition. If so, a three-dimensional interpolation algorithm is used based on sample point data to obtain grid numerical attribute information. A regular hexahedron grid model is then obtained based on layer number information, grid index number, grid position attribute information, and grid numerical attribute information. The geological body model and the regular hexahedron grid model are then intersected by triangular faces to obtain a corner grid model. Finally, the three-dimensional geological model is dynamically visualized based on the regular hexahedron grid model and the corner grid model. In the present invention, the regular hexahedron grid model and the corner grid model are obtained through the voxelization algorithm, the three-dimensional interpolation algorithm, and the geological body and regular hexahedron intersection algorithm, thereby realizing a model dynamic visualization function based on the regular hexahedron grid model and the corner grid model. BRIEF DESCRIPTION OF THE DRAWINGS

[0033] FIG1 is a schematic diagram of the structure of a three-dimensional geological model dynamic visualization device in a hardware operating environment according to an embodiment of the present invention;

[0034] FIG2 is a schematic flow chart of a first embodiment of a method for dynamic visualization of a three-dimensional geological model according to the present invention;

[0035] FIG3 is a schematic diagram of a geological body model according to a first embodiment of a method for dynamic visualization of a three-dimensional geological model of the present invention;

[0036] FIG4 is a voxelization diagram of a first embodiment of a method for dynamic visualization of a three-dimensional geological model according to the present invention;

[0037] FIG5 is a schematic diagram of sample point data of the first embodiment of the dynamic visualization method of a three-dimensional geological model according to the present invention;

[0038] FIG6 is a schematic diagram of interpolation results of the first embodiment of the method for dynamic visualization of a three-dimensional geological model according to the present invention;

[0039] FIG7 is a schematic diagram of a regular hexahedral grid model of a first embodiment of a method for dynamic visualization of a three-dimensional geological model according to the present invention;

[0040] FIG8 is a schematic diagram of a corner point grid model of a first embodiment of a method for dynamic visualization of a three-dimensional geological model according to the present invention;

[0041] FIG9 is a schematic diagram of the hierarchical and block organization of a grid model according to the first embodiment of the method for dynamic visualization of a three-dimensional geological model of the present invention;

[0042] FIG10 is a schematic diagram of a grid model data format according to a first embodiment of a method for dynamic visualization of a three-dimensional geological model of the present invention;

[0043] FIG11 is a structural block diagram of the first embodiment of the three-dimensional geological model dynamic visualization system of the present invention.

[0044] The purpose, features and advantages of the present invention will be further described with reference to the accompanying drawings and in conjunction with the embodiments. DETAILED DESCRIPTION

[0045] It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0046] 1 , which is a schematic diagram of the structure of a three-dimensional geological model dynamic visualization device in a hardware operating environment according to an embodiment of the present invention.

[0047] As shown in Figure 1, the 3D geological model dynamic visualization device may include: a processor 1001, such as a central processing unit (CPU), a communication bus 1002, a user interface 1003, a network interface 1004, and a memory 1005. The communication bus 1002 is used to enable communication between these components. The user interface 1003 may include a display and an input unit, such as a keyboard. Optionally, the user interface 1003 may also include a standard wired interface or a wireless interface. The network interface 1004 may optionally include a standard wired interface or a wireless interface (such as a Wireless Fidelity (Wi-Fi) interface). The memory 1005 may be a high-speed random access memory (RAM) or a stable non-volatile memory (NVM), such as a disk drive. Optionally, the memory 1005 may be a storage system independent of the processor 1001.

[0048] Those skilled in the art will appreciate that the structure shown in FIG1 does not limit the three-dimensional geological model dynamic visualization device, and may include more or fewer components than shown, or combine certain components, or arrange the components differently.

[0049] As shown in FIG1 , the memory 1005 as a storage medium may include an operating system, a network communication module, a user interface module, and a three-dimensional geological model dynamic visualization program.

[0050] In the 3D geological model dynamic visualization device shown in FIG1 , the network interface 1004 is primarily used for data communication with a network server; the user interface 1003 is primarily used for data interaction with a user; the processor 1001 and the memory 1005 in the 3D geological model dynamic visualization device of the present invention can be provided in the 3D geological model dynamic visualization device. The 3D geological model dynamic visualization device calls the 3D geological model dynamic visualization program stored in the memory 1005 via the processor 1001 and executes the 3D geological model dynamic visualization method provided in an embodiment of the present invention.

[0051] An embodiment of the present invention provides a method for dynamic visualization of a three-dimensional geological model. Referring to FIG. 2 , FIG. 2 is a flow chart of a first embodiment of the method for dynamic visualization of a three-dimensional geological model according to the present invention.

[0052] In this embodiment, the three-dimensional geological model dynamic visualization method includes the following steps:

[0053] Step S10: voxelizing the geological body model to obtain grid position attribute information corresponding to the geological body model.

[0054] It is easy to understand that the execution subject of this embodiment can be a three-dimensional geological model dynamic visualization system with functions such as data processing, network communication and program running, or other computer equipment with similar functions, etc., and this embodiment is not limited thereto.

[0055] It should also be noted that voxelization processing (i.e., voxelization algorithm) includes surface voxelization processing and internal voxelization processing. The voxelization algorithm determines the internal and external relationship between the grid point and the geological body by calculating the number of intersection points between the ray at the grid point location and the geological body.

[0056] In this embodiment, refer to Figure 3, which is a schematic diagram of the geological body model of the first embodiment of the dynamic visualization method of a three-dimensional geological model of the present invention. The geological body is then divided into layers, and the surface grid and the internal grid are obtained by voxelizing the surface and the interior of the geological body model triangulation network. Refer to Figure 4, which is a voxelization schematic diagram of the first embodiment of the dynamic visualization method of a three-dimensional geological model of the present invention.

[0057] It should be understood that after surface voxelization and interior voxelization are performed to obtain the surface grid and the interior grid, the position attributes of each grid need to be marked.

[0058] Step S20: Determine whether the grid position attribute information meets the preset model position condition.

[0059] It should also be understood that the preset model position condition is that the grid is on and within the geological body model.

[0060] Step S30: If satisfied, grid numerical attribute information is obtained through a three-dimensional interpolation algorithm based on the sample point data, and a regular hexahedral grid model is obtained based on the layer number information, the grid index number, the grid position attribute information and the grid numerical attribute information.

[0061] It should also be noted that the layer number information is determined based on the geological model, the grid three-dimensional coordinate information is determined based on the grid position attribute information, and the grid index number is determined based on the grid three-dimensional coordinate information. Among them, the grid unique index number (i.e., the grid index number) is obtained by calculating the three-dimensional coordinate value of the grid.

[0062] In this embodiment, the sample point data refers to Figure 5, which is a schematic diagram of the sample point data of the first embodiment of the three-dimensional geological model dynamic visualization method of the present invention. The sample point data is within the scope of the geological body model, and then the numerical attribute information of each grid is obtained through the three-dimensional interpolation algorithm. Refer to Figure 6, which is a schematic diagram of the interpolation result of the first embodiment of the three-dimensional geological model dynamic visualization method of the present invention.

[0063] The 3D interpolation algorithm supports the 3D inverse distance weighted algorithm and the 3D Kriging interpolation algorithm.

[0064] It should also be noted that the grid numerical attribute value refers to the attribute value type is a numeric type, such as 0.52, 2.6, not a character attribute value such as "big" or "small".

[0065] In a specific implementation, the layer number information, grid unique index number, grid position attribute and grid numerical attribute value are stored in a database to obtain a regular hexahedral grid model. Referring to FIG7 , FIG7 is a schematic diagram of a regular hexahedral grid model of the first embodiment of the dynamic visualization method of a three-dimensional geological model of the present invention.

[0066] Step S40: performing triangular intersection on the geological body model and the regular hexahedral grid model to obtain a corner point grid model.

[0067] Furthermore, the geological body model and the regular hexahedral mesh model are intersected by triangulation to obtain the corner point mesh model. The processing method is to intersect the geological body model and the regular hexahedral mesh model by triangulation to obtain the triangulated mesh intersection result; and the triangulated mesh intersection result is processed to obtain the corner point mesh model.

[0068] It should also be understood that the algorithm for intersecting the triangular face of the geological model with the regular hexahedron is to intersect each triangle in the triangular face of the geological model with the grid information in the regular hexahedron grid model taken out from the database to obtain the resulting triangular face (i.e., the triangulated mesh intersection result), process the intersection result and record the attribute values ​​to obtain the corner point grid model and the corresponding corner point attribute values. Refer to Figure 8, which is a schematic diagram of the corner point grid model of the first embodiment of the three-dimensional geological model dynamic visualization method of the present invention.

[0069] Step S50: dynamically visualizing the three-dimensional geological model based on the regular hexahedral grid model and the corner point grid model.

[0070] After the step of dynamically visualizing the three-dimensional geological model based on the regular hexahedron grid model and the corner point grid model, the grid model data formats of the regular hexahedron grid model and the corner point grid model are determined respectively; and the regular hexahedron grid model and the corner point grid model are stored according to the grid model data formats.

[0071] In this embodiment, the regular hexahedral grid model and the corner grid model are respectively divided into blocks using a hierarchical block organization strategy to obtain a divided regular hexahedral grid model and a divided corner grid model. Referring to FIG9 , FIG9 is a schematic diagram of the hierarchical block organization of the grid model of the first embodiment of the dynamic visualization method of a three-dimensional geological model of the present invention. The divided regular hexahedral grid model and the divided corner grid model are then stored according to the grid model data format. Referring to FIG10 , FIG10 is a schematic diagram of the grid model data format of the first embodiment of the dynamic visualization method of a three-dimensional geological model of the present invention. Refer to the following table, which is a table of the data format of the regular hexahedral model and the data format of the corner grid model, respectively:

[0072] Regular hexahedron model data format table

[0073] Corner point grid model data format table

[0074] In practice, data storage for 3D geological grid models (i.e., regular hexahedral grid models and corner grid models) is implemented using the Grid Storage class. This class supports data storage in MongoDB, HDB, and PostgreSQL databases, and supports fast block and batch storage, as well as multi-attribute storage.

[0075] The database stores the geometric parameters (including the starting point of the geological model, grid step size, number of grids, and Lod level) of the grid model (i.e., regular hexahedral grid model and corner point grid model), grid sequence number (IJK sequence set) and multi-attribute set, supports numerical attribute values ​​and flow attribute values, and does not store grid geometric coordinate information.

[0076] A hierarchical and block-based organization strategy is employed, with the grid model organized hierarchically based on the level of detail (LoD) technique. The hierarchical grid model is divided into blocks using a tree structure such as a quadtree, octree, or KD tree. Referring to Figure 9, which illustrates the hierarchical and block-based organization of the grid model in the first embodiment of the dynamic visualization method for a 3D geological model of the present invention, the center point coordinates and attribute values ​​of each regular hexahedral cell in the bounding box are stored. Internal points are assigned attribute values, external points are assigned invalid values ​​(NULL) or marked values, and boundary cells are assigned marked values. The generated regular hexahedral grid model is stored in a database, and the triangulation information of the corner grid model (external boundary point coordinates and external boundary triangulation) is stored in a local cache.

[0077] It should also be understood that the PostgreSQL database table is used for storage, the Mpdh_item metadata table mainly describes the metadata information of the grid model, and the GridName_LODX grid data table is used to store the hierarchical information.

[0078] Grid models of different LODs are stored as separate collections, named in the format GridName_LODX. The collections contain four fields: gn (grid number), layer (geological layer to which the grid belongs), valid (whether the grid is valid), and att1 (grid attribute value).

[0079] The storage interface design supports field expression updates, eliminating the need to fetch all attribute data to the client for calculation and then update back to the database. It supports queries based on multiple condition combinations, and allows for specifying query content.

[0080] Furthermore, the grid cache class generates different grid caches for different grid types. For corner grid types, the hexahedral grid is used to intersect with the original geological model triangulation network, and the result of the intersection is processed. The triangular faces within the hexahedron are assigned the attribute values ​​of the hexahedron grid; for regular hexahedrons, the corresponding grid and attributes are directly taken from the database to generate a hexahedral grid and generate a cache.

[0081] The mesh cache is optimized for storage, using the center coordinates of the mesh as the actual display coordinates. Shaders are used for triangulation on the display. The cache records relevant mesh index information, allowing for quick query, highlighting, and filtering of a single mesh.

[0082] The grid cache class also includes the processing logic of multi-level LOD grids. For each grid to be displayed, it is subdivided according to the LOD level, the internal and external relationships of the subdivided grid are determined, and the interpolation results are calculated and stored in the corresponding database table.

[0083] In this embodiment, a geological body model is first voxelized to obtain grid position attribute information corresponding to the geological body model. A determination is then made as to whether the grid position attribute information satisfies a preset model position condition. If so, a three-dimensional interpolation algorithm is used based on the sample point data to obtain grid numerical attribute information. A regular hexahedral grid model is then obtained based on the layer number information, grid index number, grid position attribute information, and grid numerical attribute information. The geological body model and the regular hexahedral grid model are then intersected by triangular faces to obtain a corner grid model. Finally, dynamic visualization of the three-dimensional geological model is performed based on the regular hexahedral grid model and the corner grid model. In this embodiment, the regular hexahedral grid model and the corner grid model are obtained through the voxelization algorithm, the three-dimensional interpolation algorithm, and the intersection algorithm between the geological body and the regular hexahedron, thereby achieving a dynamic model visualization function based on the regular hexahedral grid model and the corner grid model.

[0084] 11 , which is a structural block diagram of a first embodiment of a three-dimensional geological model dynamic visualization system according to the present invention.

[0085] As shown in FIG11 , the three-dimensional geological model dynamic visualization system proposed in an embodiment of the present invention includes:

[0086] The processing module 1101 is used to perform voxel processing on the geological body model to obtain grid position attribute information corresponding to the geological body model.

[0087] It should also be noted that voxelization processing (i.e., voxelization algorithm) includes surface voxelization processing and internal voxelization processing. The voxelization algorithm determines the internal and external relationship between the grid point and the geological body by calculating the number of intersection points between the ray at the grid point location and the geological body.

[0088] In this embodiment, refer to Figure 3, which is a schematic diagram of the geological body model of the first embodiment of the dynamic visualization method of a three-dimensional geological model of the present invention. The geological body is then divided into layers, and the surface grid and the internal grid are obtained by voxelizing the surface and the interior of the geological body model triangulation network. Refer to Figure 4, which is a voxelization schematic diagram of the first embodiment of the dynamic visualization method of a three-dimensional geological model of the present invention.

[0089] It should be understood that after surface voxelization and interior voxelization are performed to obtain the surface grid and the interior grid, the position attributes of each grid need to be marked.

[0090] The judgment module 1102 is used to judge whether the grid position attribute information meets the preset model position condition.

[0091] It should also be understood that the preset model position condition is that the grid is on and within the geological body model.

[0092] The generation module 1103 is used to obtain grid numerical attribute information through a three-dimensional interpolation algorithm based on the sample point data if the conditions are met, and to obtain a regular hexahedral grid model based on the layer number information, the grid index number, the grid position attribute information and the grid numerical attribute information.

[0093] It should also be noted that the layer number information is determined based on the geological model, the grid three-dimensional coordinate information is determined based on the grid position attribute information, and the grid index number is determined based on the grid three-dimensional coordinate information. Among them, the grid unique index number (i.e., the grid index number) is obtained by calculating the three-dimensional coordinate value of the grid.

[0094] In this embodiment, the sample point data refers to Figure 5, which is a schematic diagram of the sample point data of the first embodiment of the three-dimensional geological model dynamic visualization method of the present invention. The sample point data is within the scope of the geological body model, and then the numerical attribute information of each grid is obtained through the three-dimensional interpolation algorithm. Refer to Figure 6, which is a schematic diagram of the interpolation result of the first embodiment of the three-dimensional geological model dynamic visualization method of the present invention.

[0095] The 3D interpolation algorithm supports the 3D inverse distance weighted algorithm and the 3D Kriging interpolation algorithm.

[0096] It should also be noted that the grid numerical attribute value refers to the attribute value type is a numeric type, such as 0.52, 2.6, not a character attribute value such as "big" or "small".

[0097] In a specific implementation, the layer number information, grid unique index number, grid position attribute and grid numerical attribute value are stored in a database to obtain a regular hexahedral grid model. Referring to FIG7 , FIG7 is a schematic diagram of a regular hexahedral grid model of the first embodiment of the dynamic visualization method of a three-dimensional geological model of the present invention.

[0098] The generating module 1103 is further configured to perform triangular intersection between the geological body model and the regular hexahedral grid model to obtain a corner grid model.

[0099] Furthermore, the geological body model and the regular hexahedral mesh model are intersected by triangulation to obtain the corner point mesh model. The processing method is to intersect the geological body model and the regular hexahedral mesh model by triangulation to obtain the triangulated mesh intersection result; and the triangulated mesh intersection result is processed to obtain the corner point mesh model.

[0100] It should also be understood that the algorithm for intersecting the geological body model triangular face with the regular hexahedron is to intersect each triangle in the geological body model triangular face with the grid information in the regular hexahedron grid model retrieved from the database to obtain the resulting triangular face (i.e., the triangulated mesh intersection result), process the intersection result and record the attribute values ​​to obtain the corner point grid model and the corresponding corner point attribute values. Refer to Figure 8, which is a schematic diagram of the corner point grid model of the first embodiment of the three-dimensional geological model dynamic visualization method of the present invention.

[0101] The display module 1104 is configured to dynamically visualize the three-dimensional geological model based on the regular hexahedral grid model and the corner point grid model.

[0102] After the step of dynamically visualizing the three-dimensional geological model based on the regular hexahedron grid model and the corner point grid model, the grid model data formats of the regular hexahedron grid model and the corner point grid model are determined respectively; and the regular hexahedron grid model and the corner point grid model are stored according to the grid model data formats.

[0103] In this embodiment, the regular hexahedral grid model and the corner grid model are respectively divided into blocks using a hierarchical block organization strategy to obtain a divided regular hexahedral grid model and a divided corner grid model. Referring to FIG9 , FIG9 is a schematic diagram of the hierarchical block organization of the grid model of the first embodiment of the dynamic visualization method of a three-dimensional geological model of the present invention. The divided regular hexahedral grid model and the divided corner grid model are then stored according to the grid model data format. Referring to FIG10 , FIG10 is a schematic diagram of the grid model data format of the first embodiment of the dynamic visualization method of a three-dimensional geological model of the present invention. Refer to the following table, which is a table of the data format of the regular hexahedral model and the data format of the corner grid model, respectively:

[0104] Regular hexahedron model data format table

[0105] Corner point grid model data format table

[0106] In practice, data storage for 3D geological grid models (i.e., regular hexahedral grid models and corner grid models) is implemented using the Grid Storage class. This class supports data storage in MongoDB, HDB, and PostgreSQL databases, and supports fast block and batch storage, as well as multi-attribute storage.

[0107] The database stores the geometric parameters (including the starting point of the geological model, grid step size, number of grids, and Lod level) of the grid model (i.e., regular hexahedral grid model and corner point grid model), grid sequence number (IJK sequence set) and multi-attribute set, supports numerical attribute values ​​and flow attribute values, and does not store grid geometric coordinate information.

[0108] A hierarchical and block-based organization strategy is employed, with the grid model organized hierarchically based on the level of detail (LoD) technique. The hierarchical grid model is divided into blocks using a tree structure such as a quadtree, octree, or KD tree. Referring to Figure 9, which illustrates the hierarchical and block-based organization of the grid model in the first embodiment of the dynamic visualization method for a 3D geological model of the present invention, the center point coordinates and attribute values ​​of each regular hexahedral cell in the bounding box are stored. Internal points are assigned attribute values, external points are assigned invalid values ​​(NULL) or marked values, and boundary cells are assigned marked values. The generated regular hexahedral grid model is stored in a database, and the triangulation information of the corner grid model (external boundary point coordinates and external boundary triangulation) is stored in a local cache.

[0109] It should also be understood that the PostgreSQL database table is used for storage, the Mpdh_item metadata table mainly describes the metadata information of the grid model, and the GridName_LODX grid data table is used to store the hierarchical information.

[0110] Grid models of different LODs are stored as separate collections, named in the format GridName_LODX. The collections contain four fields: gn (grid number), layer (geological layer to which the grid belongs), valid (whether the grid is valid), and att1 (grid attribute value).

[0111] The storage interface design supports field expression updates, eliminating the need to fetch all attribute data to the client for calculation and then update back to the database. It supports queries based on multiple condition combinations, and allows for specifying query content.

[0112] Furthermore, the grid cache class generates different grid caches for different grid types. For corner grid types, the hexahedral grid is used to intersect with the original geological model triangulation network, and the result of the intersection is processed. The triangular faces within the hexahedron are assigned the attribute values ​​of the hexahedron grid; for regular hexahedrons, the corresponding grid and attributes are directly taken from the database to generate a hexahedral grid and generate a cache.

[0113] The mesh cache is optimized for storage, using the center coordinates of the mesh as the actual display coordinates. Shaders are used for triangulation on the display. The cache records relevant mesh index information, allowing for quick query, highlighting, and filtering of a single mesh.

[0114] The grid cache class also includes the processing logic of multi-level LOD grids. For each grid to be displayed, it is subdivided according to the LOD level, the internal and external relationships of the subdivided grid are determined, and the interpolation results are calculated and stored in the corresponding database table.

[0115] In this embodiment, a geological body model is first voxelized to obtain grid position attribute information corresponding to the geological body model. A determination is then made as to whether the grid position attribute information satisfies a preset model position condition. If so, a three-dimensional interpolation algorithm is used based on the sample point data to obtain grid numerical attribute information. A regular hexahedral grid model is then obtained based on the layer number information, grid index number, grid position attribute information, and grid numerical attribute information. The geological body model and the regular hexahedral grid model are then intersected by triangular faces to obtain a corner grid model. Finally, dynamic visualization of the three-dimensional geological model is performed based on the regular hexahedral grid model and the corner grid model. In this embodiment, the regular hexahedral grid model and the corner grid model are obtained through the voxelization algorithm, the three-dimensional interpolation algorithm, and the intersection algorithm between the geological body and the regular hexahedron, thereby achieving a dynamic model visualization function based on the regular hexahedral grid model and the corner grid model.

[0116] Other embodiments or specific implementations of the three-dimensional geological model dynamic visualization system of the present invention can refer to the above-mentioned method embodiments and will not be described in detail here.

[0117] It should be noted that, in this document, the terms "comprises," "includes," or any other variations thereof are intended to encompass non-exclusive inclusion, such that a process, method, article, or system comprising a series of elements includes not only those elements but also other elements not explicitly listed, or elements inherent to such process, method, article, or system. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of other identical elements in the process, method, article, or system comprising the element.

[0118] The serial numbers of the above embodiments of the present invention are for description only and do not represent the advantages or disadvantages of the embodiments.

[0119] Through the description of the above embodiments, those skilled in the art can clearly understand that the above embodiment methods can be implemented by means of software plus the necessary general hardware platform, and of course can also be implemented by hardware, but in many cases the former is a better embodiment. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art can be embodied in the form of a software product, which is stored in a storage medium (such as read-only memory / random access memory, magnetic disk, optical disk), and includes a number of instructions for enabling a terminal device (which can be a mobile phone, computer, server, or network device, etc.) to execute the methods described in each embodiment of the present invention.

[0120] The above are only preferred embodiments of the present invention and are not intended to limit the patent scope of the present invention. Any equivalent structure or equivalent process transformation made using the contents of the present invention description and drawings, or directly or indirectly applied in other related technical fields, are also included in the patent protection scope of the present invention.

Claims

1. A three-dimensional geological model dynamic visualization method, characterized in that: The three-dimensional geological model dynamic visualization method comprises the following steps: voxelize the geological body model to obtain grid position attribute information corresponding to the geological body model; Determining whether the grid position attribute information satisfies a preset model position condition; If satisfied, then the grid numerical attribute information is obtained through a three-dimensional interpolation algorithm according to the sample point data, and a regular hexahedral grid model is obtained according to the layer number information, the grid index number, the grid position attribute information and the grid numerical attribute information; Perform triangular surface intersection of the geological body model and the regular hexahedral grid model to obtain a corner point grid model; The three-dimensional geological model is dynamically visualized based on the regular hexahedral grid model and the corner point grid model.

2. The method according to claim 1, characterized in that Before the step of obtaining a regular hexahedral grid model according to the layer number information, the grid index number, the grid position attribute information and the grid value attribute information, the step further includes: Determining layer number information according to the geological body model, and determining grid three-dimensional coordinate information according to the grid position attribute information; A grid index number is determined according to the grid three-dimensional coordinate information.

3. The method according to claim 2, characterized in that The step of performing triangular surface intersection of the geological body model and the regular hexahedral mesh model to obtain a corner point mesh model comprises: The geological body model and the regular hexahedral mesh model are intersected by triangular faces to obtain a triangular Angular grid intersection results; The triangulated network intersection result is processed to obtain a corner point mesh model.

4. The method according to any one of claims 1 to 3, characterized in that: After the step of dynamically visualizing the three-dimensional geological model based on the regular hexahedral grid model and the corner point grid model, the method further includes: Determining the mesh model data formats of the regular hexahedral mesh model and the corner point mesh model respectively; The regular hexahedral mesh model and the corner point mesh model are stored according to the mesh model data format.

5. The method according to claim 4, characterized in that The step of storing the regular hexahedral mesh model and the corner point mesh model according to the mesh model data format comprises: The regular hexahedral mesh model and the corner point mesh model are respectively divided into blocks by a hierarchical block organization strategy to obtain a divided regular hexahedral mesh model and a divided corner point mesh model; The divided regular hexahedral mesh model and the divided corner point mesh model are stored according to the mesh model data format.

6. A three-dimensional geological model dynamic visualization system, characterized in that: The three-dimensional geological model dynamic visualization system includes: A processing module, used for voxelizing the geological body model to obtain grid position attribute information corresponding to the geological body model; A judgment module, used to judge whether the grid position attribute information meets the preset model position condition; The generation module is used to obtain the grid through the three-dimensional interpolation algorithm based on the sample point data if it meets Numerical attribute information, and obtaining a regular hexahedral grid model according to the layer number information, the grid index number, the grid position attribute information and the grid numerical attribute information; The generation module is further used to perform triangular surface intersection between the geological body model and the regular hexahedral grid model to obtain a corner point grid model; A display module is used to dynamically visualize the three-dimensional geological model based on the regular hexahedral grid model and the corner point grid model.

7. A three-dimensional geological model dynamic visualization device, characterized in that: The device comprises: a memory, a processor, and a three-dimensional geological model dynamic visualization program stored in the memory and executable on the processor, wherein the three-dimensional geological model dynamic visualization program is configured to implement the steps of the three-dimensional geological model dynamic visualization method according to any one of claims 1 to 5.

8. A storage medium, characterized in that: The storage medium stores a three-dimensional geological model dynamic visualization program, which, when executed by a processor, implements the steps of the three-dimensional geological model dynamic visualization method according to any one of claims 1 to 5.

Images