A method and system for fusing irregular geologic anomalies with three-dimensional geology
By drawing curved surfaces using half-side rules and constructing meshes using octrees, and combining oriented bounding boxes to determine intersecting regular hexahedrons for Boolean operations, the problem of generating and updating irregular anomalies in 3D geology was solved. This enabled the rapid fusion and realistic display of geological anomalies, improving the accuracy and efficiency of advanced geological forecasting.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- SHANDONG UNIV
- Filing Date
- 2023-10-13
- Publication Date
- 2026-06-26
AI Technical Summary
In 3D geological simulation, it is difficult to efficiently generate and update irregular geological anomalies in real time. In particular, it is difficult to set up irregular anomalies in a 3D environment and to display them realistically in real time. Furthermore, Boolean operations are inefficient.
The surface of the 3D point cloud prediction volume is drawn using half-side rules. The 3D geological hexahedral mesh is constructed using an octree. Intersecting hexahedrals are determined by oriented bounding boxes, and polyhedral Boolean operations are performed to quickly fuse irregular geological anomalies with 3D geology.
It enables rapid fusion and realistic display of irregular geological anomalies in three-dimensional geology, improving the accuracy and efficiency of advanced geological forecasting.
Smart Images
Figure CN117574482B_ABST
Abstract
Description
Technical Field
[0001] This invention belongs to the field of advanced geological prediction technology, and relates to a method and system for fusing irregular geological anomalies with three-dimensional geology. Background Technology
[0002] The statements in this section are merely background information related to the present invention and do not necessarily constitute prior art.
[0003] Currently, in simulation experiments for advanced geological prediction data, numerical simulation models are often constructed. These models include three-dimensional geological features and geological anomalies with different wave velocities, such as irregular karst caves, faults, and fracture zones. Two-dimensional stratigraphic and anomaly models are relatively easy to set up; a specific irregular region can be selected on a two-dimensional cross-section. However, three-dimensional environmental simulation geological models are relatively complex, especially in seismic finite difference models. Setting up irregular anomalies within thousands of hexahedral elements is quite difficult because regular three-dimensional anomalies can be obtained using equations for spheres, ellipsoids, etc. Irregular anomalies cannot be generated through equations alone; they must be generated through editing using specialized drawing tools, a complex process, such as creating an irregular three-dimensional fault. Assembling irregular anomalies by selecting individual cubes is extremely difficult.
[0004] In addition, in the actual construction of underground tunnels, the process of displaying the three-dimensional results of the adverse ground bodies ahead, which are measured by the advanced geological prediction method, in real time and realistically in the three-dimensional geology requires the real-time fusion and updating of the three-dimensional anomalies and the three-dimensional geology. This requires the use of boolean operations, which require finding the intersection points by intersecting the three-dimensional anomaly surface with the mesh after the three-dimensional geology is partitioned one by one. This is inefficient and quite difficult. Summary of the Invention
[0005] To address the aforementioned problems, this invention proposes a method and system for fusing irregular geological anomalies with three-dimensional geology. This invention enables the rapid integration of advanced detection results into the preceding three-dimensional geological data during tunnel excavation, providing a more intuitive and realistic representation of geological anomalies within the preceding strata. Furthermore, it allows for the setting of various irregular geological anomalies in numerical simulation experiments for advanced geological prediction, enabling numerical simulation experiments of multiple prediction methods and verifying their accuracy.
[0006] According to some embodiments, the present invention adopts the following technical solution:
[0007] A method for fusing irregular geological anomalies with three-dimensional geology includes the following steps:
[0008] The surface of the predicted volume of 3D point cloud forecast is drawn using half-side rules;
[0009] Using an octree, a three-dimensional geological hexahedral mesh is constructed;
[0010] Find the bounding box of the 3D surface formed by surface rendering, and obtain the regular hexahedron inside the bounding box;
[0011] Determine whether all regular hexahedrons intersect with the 3D curved surface, and identify all intersecting hexahedrons;
[0012] Perform polyhedral Boolean operations on all intersecting hexahedrons and 3D surfaces;
[0013] The results of Boolean operations on polyhedra are integrated into a three-dimensional geological hexahedral mesh.
[0014] As an alternative implementation method, the specific process of constructing a three-dimensional geological hexahedral mesh using an octree includes:
[0015] (1) Set the maximum recursion depth;
[0016] (2) Find the maximum size of the scene and build the first cube based on this size;
[0017] (3) Place the identity elements into the cube that can be contained and has no child nodes in sequence;
[0018] (4) If the maximum recursion depth has not been reached, continue to subdivide into eight equal parts, and then distribute all the unit elements contained in the cube to the eight sub-cubes.
[0019] (5) If it is found that the number of unit elements assigned to the sub-cube is not zero and is the same as that of the parent cube, then the sub-cube will stop subdividing.
[0020] (6) Repeat step (3) until the maximum recursion depth is reached.
[0021] As an alternative implementation method, the specific process of obtaining the oriented bounding box of the three-dimensional surface formed by surface rendering includes obtaining the eigenvectors, i.e. the principal axes of the oriented bounding box, through principal component analysis based on the vertices of the object surface, and obtaining the lowest bounding box by calculating the covariance matrix of the three-dimensional surface, so that the oriented bounding box is along the direction of the eigenvector corresponding to the largest eigenvalue.
[0022] As an alternative implementation, the specific process of obtaining the regular hexahedron inside the orientation bounding box includes: arbitrarily selecting any point P in space, calculating the dot product of the vector DP between one endpoint D of the regular hexahedron inside the orientation bounding box and the vector DE between that endpoint and its horizontally adjacent endpoint E, moving the vector PE to endpoint D, calculating the dot product of vector PE and vector DE, and if the product of the two dot products is greater than zero, then point P is outside the cube; if it is equal to zero, then it is on the cube; if it is less than zero, then it is inside the cube.
[0023] As an alternative implementation, the specific process of determining whether all regular hexahedrons intersect with the three-dimensional curved surface, and identifying all intersecting hexahedrons, includes:
[0024] 1) Establish the positions and index numbers of the six points of all regular hexahedrons inside the oriented bounding box;
[0025] 2) Select one of the regular hexahedrons, and emit a line from a point of the regular hexahedron to positive infinity. Find the intersection point of the ray with the three-dimensional surface of the prediction result. If the number of intersection points is odd, then the point is in the three-dimensional surface of the prediction result.
[0026] 3) Repeat step 2 until all six points of the regular hexahedron have been traversed. If all six points are inside the three-dimensional surface or all six points are outside the three-dimensional surface, then the regular hexahedron does not intersect with the predicted three-dimensional surface; otherwise, they intersect.
[0027] 4) Repeat step 3) to traverse all the regular hexahedrons inside the oriented bounding box and record all the regular hexahedrons that intersect with the 3D surface of the prediction result.
[0028] As an alternative implementation, the specific process of performing polyhedral Boolean operations on all intersecting hexahedrons and three-dimensional surfaces includes:
[0029] 1) In each triangular mesh, extract the set of edge segments from the edges, and extract the set of triangles from the face;
[0030] 2) Filter out the area where the edge segment of one triangular mesh intersects with another triangular mesh;
[0031] 3) When a line segment intersects with the interior of a triangle, calculate the intersection point. For each face related to a line segment, save the relationship between the intersection point and the face of the intersecting triangle.
[0032] 4) When a line segment intersects with a side of a triangle, calculate the intersection point and retain the triangle associated with that side;
[0033] 5) When a line segment intersects a vertex of a triangle, calculate the intersection point and retain the triangle associated with that line segment.
[0034] A system integrating irregular geological anomalies with three-dimensional geology, comprising:
[0035] The surface rendering module is configured to render the predicted volume of the 3D point cloud using half-edge rules.
[0036] The network building module is configured to use an octree to construct a three-dimensional geological hexahedral mesh;
[0037] The acquisition module is configured to obtain the oriented bounding box of the 3D surface formed by surface drawing, and to obtain the regular hexahedron inside the oriented bounding box;
[0038] The intersection extraction module is configured to determine whether all regular hexahedrons intersect with the three-dimensional curved surface, and to identify all intersecting hexahedrons;
[0039] The operation module is configured to perform polyhedral Boolean operations on all intersecting hexahedrons and three-dimensional surfaces;
[0040] The fusion module is configured to fuse the results of polyhedral Boolean operations into a 3D geological hexahedral mesh.
[0041] A computer-readable storage medium storing a plurality of instructions adapted for loading by a processor of a terminal device and executing steps in the method.
[0042] A terminal device includes a processor and a computer-readable storage medium, the processor being configured to implement instructions; the computer-readable storage medium being configured to store a plurality of instructions adapted to be loaded by the processor and executed in accordance with the steps of the method described therein.
[0043] Compared with the prior art, the beneficial effects of the present invention are as follows:
[0044] This invention enables the rapid integration of our advanced detection results into the three-dimensional geological data ahead during tunnel excavation, providing a more intuitive and realistic representation of geological anomalies within the strata. Furthermore, it allows for the setting of various irregular geological anomalies in numerical simulation experiments for advanced geological prediction, enabling numerical simulation experiments of multiple prediction methods and verifying their accuracy.
[0045] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, preferred embodiments are described below in detail with reference to the accompanying drawings. Attached Figure Description
[0046] The accompanying drawings, which form part of this invention, are used to provide a further understanding of the invention. The illustrative embodiments of the invention and their descriptions are used to explain the invention and do not constitute an improper limitation of the invention.
[0047] Figure 1 This is a flowchart of Example 1;
[0048] Figure 2 This is a schematic diagram of the generation of a three-dimensional surface of a karst cave based on half-edge rules in Example 1;
[0049] Figure 3 This is a schematic diagram of the irregular fault three-dimensional anomaly surface drawn based on half-side rules in Example 1;
[0050] Figure 4 This is a schematic diagram of the octree organizational structure in Example 1;
[0051] Figure 5 This is a schematic diagram of the generation of a three-dimensional geological hexahedral mesh based on an octree in Example 1;
[0052] Figure 6 This is a schematic diagram of the OBB calculation principle in Example 1;
[0053] Figure 7 This is an OBB diagram showing the advanced geological prediction results of Example 1;
[0054] Figure 8 This is a flowchart of the process for obtaining all the regular hexahedrons inside the OBB bounding box in Example 1;
[0055] Figure 9 This is a schematic diagram of the process of obtaining all the regular hexahedrons inside the OBB bounding box in Example 1;
[0056] Figure 10 This is a schematic diagram of the process of constructing polygonal line segments by finding intersections in Example 1;
[0057] Figure 11 This is a schematic diagram of the Boolean operation result of the polyhedron in Example 1;
[0058] Figure 12 This is a demonstration of the Boolean operation performed on the three-dimensional tomographic surface and the intersecting regular hexahedron in Example 1;
[0059] Figure 13 This is a schematic diagram of the fault and tunnel fusion result in Example 1;
[0060] Figure 14 This is an actual effect diagram of Example 1;
[0061] Figure 15 This is an actual effect diagram of Example 1;
[0062] Figure 16 This is an actual effect diagram of Example 1;
[0063] Figure 17 This is an actual effect diagram of Example 1. Detailed Implementation
[0064] The present invention will be further described below with reference to the accompanying drawings and embodiments.
[0065] It should be noted that the following detailed description is illustrative and intended to provide further explanation of the invention. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains.
[0066] It should be noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the scope of exemplary embodiments according to the invention. As used herein, the singular form is intended to include the plural form as well, unless the context clearly indicates otherwise. Furthermore, it should be understood that when the terms "comprising" and / or "including" are used in this specification, they indicate the presence of features, steps, operations, devices, components, and / or combinations thereof.
[0067] Example 1
[0068] A method for the rapid fusion of advanced geological prediction of irregular geological anomalies with three-dimensional geology aims to present the three-dimensional geological strata information and geological anomaly information in front of the tunnel as realistically and intuitively as possible.
[0069] In this embodiment, the three-dimensional surface generated based on advanced geological prediction results generated by human experience is re-meshed according to the half-edge rule, and its OBB bounding box is obtained. Using the OBB bounding box rule, the hexahedrons inside the bounding box are obtained. Then, each vertex of the hexahedron inside the bounding box is intersected with the geological anomaly mesh by ray calculation to obtain the hexahedrons that intersect with the geological anomaly mesh. Finally, Boolean operations are performed to achieve rapid integration of the irregular geological anomaly from advanced geological prediction with three-dimensional geology.
[0070] In practical applications, it can provide a basic environmental model for numerical simulation of advanced geological forecasting.
[0071] As a typical embodiment, such as Figure 1 As shown, it includes the following steps:
[0072] 1. Surface rendering of 3D point cloud forecast results based on half-edge rules;
[0073] 2. Construct an octree-based three-dimensional geological hexahedral mesh subdivision structure;
[0074] 3. Obtain the OBB (Oriented Bounding Bix) bounding box of the three-dimensional surface of the forecast result, and obtain the regular hexahedron inside the OBB bounding box;
[0075] 4. Determine whether all the regular hexahedrons inside the OBB bounding box intersect with the predicted 3D surface, and record all intersecting hexahedrons;
[0076] 5. Perform a cgal Boolean operation on all intersecting hexahedrons and the predicted surface;
[0077] 6. Integrate the calculation results into the three-dimensional geological hexahedral mesh model.
[0078] In this embodiment, the predicted volume of advanced geological forecasting is displayed using a three-dimensional surface based on half-edge rules. This process can be achieved using existing technologies, such as... Figure 2 and Figure 3 As shown.
[0079] When constructing an octree-based 3D geological hexahedral mesh organization structure, an octree is a tree-like data structure used to describe 3D space. Each node in an octree represents a volume element of a cube, and each node has eight child nodes. The sum of the volume elements represented by these eight child nodes equals the volume of the parent node. The center point is typically used as the branching center of the node. If an octree is not empty, any node in the tree will have exactly eight or zero child nodes; that is, there will be no child nodes other than 0 and 8.
[0080] The implementation principle of Octree:
[0081] 1) Set the maximum recursion depth;
[0082] 2) Find the maximum size of the scene, and build the first cube using this size;
[0083] 3) Place the identity elements into the cube that can be contained and has no child nodes in sequence;
[0084] 4) If the maximum recursion depth has not been reached, subdivide the cube into eight equal parts, and then distribute all the unit elements contained in the cube to the eight sub-cubes.
[0085] 5) If it is found that the number of unit elements assigned to the sub-cube is not zero and is the same as that of the parent cube, then the sub-cube should stop subdividing. This is because, according to the theory of space partitioning, the space obtained by subdivision must be less. If the number is the same, then no matter how many times it is cut, the number will still be the same, which will result in an infinite number of cuts.
[0086] 6) Repeat step 3 until the maximum recursion depth is reached.
[0087] Specific effects are as follows Figure 5 As shown.
[0088] The process of obtaining the OBB bounding box of the 3D surface of the forecast result and acquiring the regular hexahedrons inside the OBB bounding box is as follows: In simple terms, the OBB generation idea is to obtain the eigenvectors (principal component analysis) from the vertices of the object's surface using PCA (principal component analysis). Principal component analysis is a method that transforms a set of potentially correlated variables into a set of linearly uncorrelated variables (principal components) through orthogonal transformation. The bounding box is obtained by calculating the covariance matrix of the 3D surface. Covariance represents the degree of linear correlation between two variables. The smaller the covariance, the more independent the two variables are, i.e., the weaker the linear correlation.
[0089] cov(X i X j )=E[(X i -μ i (X) j -μ j )]
[0090] The covariance matrix can be obtained using the formula for calculating covariance.
[0091]
[0092] The elements on the main diagonal represent the variance of the variables. Larger elements on the main diagonal indicate stronger signals. The elements off-diagonal represent the covariance between the variables. Larger off-diagonal elements indicate data distortion. To reduce distortion, the linear combination between variables can be redefined by diagonalizing the covariance matrix. The elements of the covariance matrix are real numbers and symmetric. The eigenvectors of the covariance matrix represent the orientation of the OBB bounding box. Larger eigenvalues correspond to larger variances, so the OBB bounding box should be aligned with the direction of the eigenvector corresponding to the largest eigenvalue.
[0093] The results of advanced geological prediction are obtained by OBB calculation as follows: Figure 7 As shown.
[0094] The process of obtaining all the regular hexahedrons inside the OBB bounding box is as follows: Figure 8 As shown, arbitrarily choose a point P in space. Calculate the dot product of vector DP between an endpoint D of the cube inside the bounding box and the point P, and vector DE between that endpoint and its horizontally adjacent endpoint E. Move vector PE to endpoint D, and calculate the dot product of vectors PE and DE. If the product of the two dot products is greater than zero, then point P is outside the cube; if it is equal to zero, then it is on the cube; and if it is less than zero, then it is inside the cube. The specific effect is as follows: Figure 9 As shown.
[0095] The process involves determining whether all the regular hexahedrons inside the OBB bounding box intersect with the predicted 3D surface, and recording all intersecting hexahedrons. The specific implementation process is as follows:
[0096] 1) Establish the positions and index numbers of the six points of all regular hexahedrons inside the OBB bounding box;
[0097] 2) Select one of the regular hexahedrons, and emit a line from a point of the regular hexahedron to positive infinity (set the x, y, z positions of the point at infinity as (100000, 100000, 100000)). Find the intersection point of the ray with the three-dimensional surface of the prediction result. If the number of intersection points is odd, then the point is in the three-dimensional surface of the prediction result.
[0098] 3) Repeat step 2 until all six points of the regular hexahedron have been traversed. If all six points are inside the three-dimensional surface or all six points are outside the three-dimensional surface, then the regular hexahedron does not intersect with the predicted three-dimensional surface; otherwise, they intersect.
[0099] 4) Repeat step 3 to traverse all the regular hexahedrons inside the OBB bounding box and record all the regular hexahedrons that intersect with the predicted 3D surface.
[0100] like Figure 10-12 As shown, the specific process of performing cgal Boolean operations on all intersecting hexahedrons and the predicted surface includes:
[0101] 1) Within each triangular mesh, we extract the segments set from the edges and the triangle set from the faces;
[0102] 2) Use Box_intersection_d to filter out the area where one triangular mesh segment intersects with another triangular mesh;
[0103] 3) Segment and triangle interior intersection: We calculate the intersection point, and for each face associated with the line segment, we save the relationship between the intersection point and the face of the intersecting triangle;
[0104] 4) The segment intersects with the sides of the triangle: Calculate the intersection point and retain the triangle associated with this side;
[0105] 5) When the segment intersects with a triangle vertex: Perform the same operation as in 4).
[0106] 6) Coplanar triangles are filtered out and processed separately.
[0107] Finally, as Figure 13 As shown, the calculation results are integrated into a three-dimensional geological hexahedral mesh model. Figure 14-17 The above embodiments are shown in the actual effect diagram when applied to a project.
[0108] The present invention also provides the following product examples:
[0109] A system integrating irregular geological anomalies with three-dimensional geology, comprising:
[0110] The surface rendering module is configured to render the predicted volume of the 3D point cloud using half-edge rules.
[0111] The network building module is configured to use an octree to construct a three-dimensional geological hexahedral mesh;
[0112] The acquisition module is configured to obtain the oriented bounding box of the 3D surface formed by surface drawing, and to obtain the regular hexahedron inside the oriented bounding box;
[0113] The intersection extraction module is configured to determine whether all regular hexahedrons intersect with the three-dimensional curved surface, and to identify all intersecting hexahedrons;
[0114] The operation module is configured to perform polyhedral Boolean operations on all intersecting hexahedrons and three-dimensional surfaces;
[0115] The fusion module is configured to fuse the results of polyhedral Boolean operations into a 3D geological hexahedral mesh.
[0116] A computer-readable storage medium storing a plurality of instructions adapted for loading by a processor of a terminal device and executing steps in the method.
[0117] A terminal device includes a processor and a computer-readable storage medium, the processor being configured to implement instructions; the computer-readable storage medium being configured to store a plurality of instructions adapted to be loaded by the processor and executed in accordance with the steps of the method described therein.
[0118] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.
[0119] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.
[0120] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.
[0121] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.
[0122] The above description is merely a preferred embodiment of the present invention and is not intended to limit the invention. Various modifications and variations can be made to the present invention by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.
[0123] While the specific embodiments of the present invention have been described above in conjunction with the accompanying drawings, this is not intended to limit the scope of protection of the present invention. Those skilled in the art should understand that various modifications or variations that can be made by those skilled in the art without creative effort based on the technical solutions of the present invention are still within the scope of protection of the present invention.
Claims
1. A method for fusing irregular geological anomalies with three-dimensional geology, characterized in that, Includes the following steps: The surface of the predicted volume of 3D point cloud forecast is drawn using half-side rules; Using an octree, a three-dimensional geological hexahedral mesh is constructed; Find the bounding box of the 3D surface formed by surface rendering, and obtain the regular hexahedron inside the bounding box; Determine whether all regular hexahedrons intersect with the 3D curved surface, and identify all intersecting hexahedrons; Perform polyhedral Boolean operations on all intersecting hexahedrons and 3D surfaces; The results of Boolean operations on polyhedra are integrated into a three-dimensional geological hexahedral mesh.
2. The method for fusing irregular geological anomalies with three-dimensional geology as described in claim 1, characterized in that, The specific process of constructing a three-dimensional geological hexahedral mesh using an octree includes: (1) Set the maximum recursion depth; (2) Find the maximum size of the scene and build the first cube based on this size; (3) Place the identity elements into the cube that can be contained and has no child nodes in sequence; (4) If the maximum recursion depth has not been reached, continue to subdivide into eight equal parts, and then distribute all the unit elements contained in the cube to the eight sub-cubes. (5) If it is found that the number of unit elements assigned to the sub-cube is not zero and is the same as that of the parent cube, then the sub-cube will stop subdividing. (6) Repeat step (3) until the maximum recursion depth is reached.
3. The method for fusing irregular geological anomalies with three-dimensional geology as described in claim 1, characterized in that, The specific process of obtaining the bounding box of the 3D surface formed by surface rendering includes obtaining the eigenvectors, i.e. the principal axes of the bounding box, through principal component analysis based on the vertices of the object surface, and obtaining the lowest bounding box by calculating the covariance matrix of the 3D surface, so that the bounding box is along the direction of the eigenvector corresponding to the largest eigenvalue.
4. The method for fusing irregular geological anomalies with three-dimensional geology as described in claim 1, characterized in that, The specific process of obtaining a regular hexahedron inside an oriented bounding box includes: arbitrarily selecting any point P in space, calculating the dot product of the vector DP between one endpoint D of the regular hexahedron inside the oriented bounding box and the vector DE between that endpoint and its horizontally adjacent endpoint E, moving vector PE to endpoint D, and calculating the dot product of vector PE and vector DE. If the product of the two dot products is greater than zero, then point P is outside the cube; if it is equal to zero, then it is on the cube; and if it is less than zero, then it is inside the cube.
5. The method for fusing irregular geological anomalies with three-dimensional geology as described in claim 1, characterized in that, The specific process for determining whether all regular hexahedrons intersect with a three-dimensional surface, and identifying all intersecting hexahedrons, includes: 1) Establish the positions and index numbers of the six points of all regular hexahedrons inside the oriented bounding box; 2) Select one of the regular hexahedrons, and emit a line from a point of the regular hexahedron to positive infinity. Find the intersection point of the ray with the three-dimensional surface of the prediction result. If the number of intersection points is odd, then the point is in the three-dimensional surface of the prediction result. 3) Repeat step 2 until all six points of the regular hexahedron have been traversed. If all six points are inside the three-dimensional surface or all six points are outside the three-dimensional surface, then the regular hexahedron does not intersect with the predicted three-dimensional surface; otherwise, they intersect. 4) Repeat step 3) to traverse all the regular hexahedrons inside the oriented bounding box and record all the regular hexahedrons that intersect with the 3D surface of the prediction result.
6. The method for fusing irregular geological anomalies with three-dimensional geology as described in claim 1, characterized in that, The specific process of performing polyhedral Boolean operations on all intersecting hexahedrons and three-dimensional surfaces includes: 1) In each triangular mesh, extract the set of edge segments from the edges, and extract the set of triangles from the face; 2) Filter out the area where the edge segment of one triangular mesh intersects with another triangular mesh; 3) When a line segment intersects with the interior of a triangle, calculate the intersection point. For each face related to a line segment, save the relationship between the intersection point and the face of the intersecting triangle. 4) When a line segment intersects with a side of a triangle, calculate the intersection point and retain the triangle associated with that side; 5) When a line segment intersects a vertex of a triangle, calculate the intersection point and retain the triangle associated with that line segment.
7. The method for fusing irregular geological anomalies with three-dimensional geology as described in claim 6, characterized in that, Coplanar triangles are treated separately.
8. A system integrating irregular geological anomalies with three-dimensional geology, characterized in that it includes: The surface rendering module is configured to render the predicted volume of the 3D point cloud using half-edge rules. The network building module is configured to use an octree to construct a three-dimensional geological hexahedral mesh; The acquisition module is configured to obtain the oriented bounding box of the 3D surface formed by surface drawing, and to obtain the regular hexahedron inside the oriented bounding box; The intersection extraction module is configured to determine whether all regular hexahedrons intersect with the three-dimensional curved surface, and to identify all intersecting hexahedrons; The operation module is configured to perform polyhedral Boolean operations on all intersecting hexahedrons and three-dimensional surfaces; The fusion module is configured to fuse the results of polyhedral Boolean operations into a 3D geological hexahedral mesh.
9. A computer-readable storage medium, characterized in that, It stores multiple instructions adapted for loading by the processor of a terminal device and executing the steps of the method according to any one of claims 1-7.
10. A terminal device, characterized in that, It includes a processor and a computer-readable storage medium, the processor being used to implement various instructions; the computer-readable storage medium being used to store a plurality of instructions adapted to be loaded by the processor and executed in the steps of the method of any one of claims 1-7.