Method and system for calculating fault plane segmentation rock stratum and electronic equipment

Abstract

The invention provides a calculation method and system for a fault plane segmented rock stratum and electronic equipment, and relates to the technical field of geological modeling, and the method comprises the steps: obtaining a rock stratum model; obtaining a fault plane model; performing interactive analysis on the rock stratum model and the fault plane model, calculating an intersection point or an intersection line of the rock stratum and the fault plane, and determining a rock stratum model after the fault plane is segmented based on the intersection point or the intersection line of the rock stratum and the fault plane; according to the geometrical characteristics of the fault plane, calculating the offset of the rock stratum model after the fault plane is segmented; and outputting the rock stratum model after the fault plane segmentation and the offset of the rock stratum model for mine pressure analysis, oil and gas simulation and disaster early warning. The method can adapt to a more complex curved surface fault plane and an irregular rock stratum structure, and high precision of a segmentation result is ensured. In a complex geological scene, efficient processing of the relation between a complex fault plane and a rock stratum is achieved, and efficient and accurate geological modeling requirements are met.

Description

Technical Field

[0001] The present invention relates to the technical field of geological modeling, and particularly to a calculation method, system and electronic device for dividing rock formations by fault planes. Background Art

[0002] With the development of modern geological exploration technologies, in transparent geological projects (such as coal mining, oil and gas exploration, etc.), simulating the tunneling process of the coal seam working face and monitoring the rock pressure and oil and gas state of the rock formations have become important links. These simulation processes can effectively predict potential disaster risks and achieve timely early warnings.

[0003] However, related technologies often rely on simple grid division and two-dimensional calculation methods, and are unable to effectively handle the division and optimization of three-dimensional complex geological bodies. Especially when the rock formation interacts with the fault plane, the difficulty of intersection point calculation and geometric optimization increases significantly. In the face of complex fault planes, it is impossible to quickly and effectively recalculate and obtain the rock formation model divided by the fault plane. Especially when dealing with the case of a curved fault plane, the calculation efficiency is low and it is difficult to meet the real-time requirements.

[0004] Therefore, there is an urgent need for a geometric algorithm that can improve the calculation efficiency of dividing rock formations by fault planes, especially in complex geological structures, to meet the needs of efficient geological modeling. Summary of the Invention

[0005] The present invention provides a calculation method, system and electronic device for dividing rock formations by fault planes, aiming to solve the problems existing in the above background art.

[0006] To solve the above technical problems, the present invention is implemented as follows: In a first aspect, the present invention provides a calculation method for dividing rock formations by fault planes, the method comprising: Obtaining a rock formation model, which is used to describe the geometric characteristics of the rock formation; Obtaining a fault plane model, which is used to describe the geometric characteristics of the fault plane, and the fault plane includes a curved fault plane; Performing interactive analysis on the rock formation model and the fault plane model, calculating the intersection points or intersection lines between the rock formation and the fault plane, and determining the rock formation model divided by the fault plane based on the intersection points or intersection lines between the rock formation and the fault plane; Calculating the offset of the rock formation model divided by the fault plane according to the geometric characteristics of the fault plane; Outputting the rock formation model divided by the fault plane and its offset for use in rock pressure analysis, oil and gas simulation and disaster early warning.

[0007] Optionally, the performing interactive analysis on the rock formation model and the fault plane model, calculating the intersection points or intersection lines between the rock formation and the fault plane, includes: Perform mesh generation on the rock stratum model and the fault plane model to obtain the first mesh elements of the rock stratum model and the fault plane model respectively; Based on the intersection points between the first mesh elements of the rock stratum model and the fault plane model respectively, determine the intersection points or intersection lines of the rock stratum and the fault plane by means of a three-dimensional geometric intersection algorithm.

[0008] Optionally, the determining the rock stratum model after the fault plane divides, based on the intersection points or intersection lines of the rock stratum and the fault plane, includes: Determine the dividing surface according to the intersection points or intersection lines of the rock stratum and the fault plane, and use a geometric cutting algorithm to divide the rock stratum model to obtain multiple independent rock stratum blocks; Perform geometric optimization on the rock stratum blocks located on the dividing surface to obtain the rock stratum model after the fault plane divides.

[0009] Optionally, the performing geometric optimization on the rock stratum blocks located on the dividing surface to obtain the rock stratum model after the fault plane divides, includes: For the dividing points on the dividing surface, reconstruct the triangular mesh to obtain the second mesh elements corresponding to the rock stratum blocks located on the dividing surface; Adjust the second mesh elements corresponding to the rock stratum blocks located on the dividing surface by means of an adaptive mesh refinement technique to obtain the rock stratum model after the fault plane divides.

[0010] Optionally, the calculating the offset of the rock stratum model after the fault plane divides according to the geometric characteristics of the fault plane, includes: According to a preset given value or the geometric characteristics of the fault plane, perform block-by-block calculation on the rock stratum model after the fault plane divides by means of a geometric transformation algorithm to obtain the offset of each rock stratum block after the fault plane divides.

[0011] Optionally, the method further includes: Use a BVH tree to perform spatial partitioning on the first mesh elements of the rock stratum model and the fault plane model respectively to obtain a plurality of bounding volume structures, and each bounding volume structure includes one or more first mesh elements; Based on the plurality of bounding volumes, perform interactive analysis on the rock stratum model and the fault plane model, and calculate the intersection points or intersection lines of the rock stratum and the fault plane, so as to improve the calculation speed and accuracy of the fault plane cutting the rock stratum.

[0012] Optionally, the obtaining the rock stratum model includes: Obtain geological exploration data, and the geological exploration data includes the geometric characteristics of each rock stratum; According to the geological exploration data, use data interpolation technology and mesh generation technology to digitally process the geometric shape, physical properties and their spatial distribution of the rock stratum to obtain the rock stratum model; The obtaining of the fault plane model includes: Obtaining geological exploration data, where the geological exploration data includes the geometric characteristics of the fault plane; According to the geological exploration data, using data interpolation technology and grid division technology, digitize the geometric shape and position of the fault plane to obtain the fault plane model.

[0013] Optionally, the method further includes: Through multi-threading technology or GPU acceleration technology, allocate computing tasks to multiple processing units for parallel computing to support large-scale parallel computing and meet the real-time simulation requirements of complex geological models.

[0014] In a second aspect, the present invention provides a computing system for fault plane-segmented rock formations, and the system includes: A first input module for obtaining a rock formation model, where the rock formation model is used to describe the geometric characteristics of the rock formation; A second input module for obtaining a fault plane model, where the fault plane model is used to describe the geometric characteristics of the fault plane; An intersection module for performing interactive analysis on the rock formation model and the fault plane model, calculating the intersection points or intersection lines of the rock formation and the fault plane, and determining the rock formation model after fault plane segmentation based on the intersection points or intersection lines of the rock formation and the fault plane; An offset calculation module for calculating the offset of the rock formation model after fault plane segmentation according to the geometric characteristics of the fault plane; An output module for outputting the rock formation model after fault plane segmentation and its offset for use in mine pressure analysis, oil and gas simulation, and disaster warning.

[0015] Optionally, the intersection module includes: A first grid division module for performing grid division on the rock formation model and the fault plane model to obtain the respective first grid cells of the rock formation model and the fault plane model; An intersection sub-module for determining the intersection points or intersection lines of the rock formation and the fault plane based on the intersection points between the respective first grid cells of the rock formation model and the fault plane model through a three-dimensional geometric intersection algorithm.

[0016] Optionally, the intersection module includes: A segmentation plane determination sub-module for determining a segmentation plane according to the intersection points or intersection lines of the rock formation and the fault plane, and using a geometric cutting algorithm to segment the rock formation model to obtain multiple independent rock formation blocks; A geometric optimization sub-module for performing geometric optimization on the rock formation blocks located on the segmentation plane to obtain the rock formation model after fault plane segmentation.

[0017] Optionally, the geometric optimization sub-module includes: A construction unit, configured to reconstruct a triangular mesh for the segmentation points on the segmentation surface, and obtain a second mesh unit corresponding to the rock formation block located on the segmentation surface; An adjustment unit, configured to adjust the second mesh unit corresponding to the rock formation block located on the segmentation surface through an adaptive mesh refinement technique, and obtain a rock formation model after fault plane segmentation.

[0018] Optionally, the offset calculation module includes: A block-by-block calculation sub-module, configured to perform block-by-block calculation on the rock formation model after fault plane segmentation through a geometric transformation algorithm according to a preset given value or geometric characteristics of the fault plane, and obtain the offset of each rock formation block after fault plane segmentation.

[0019] Optionally, the system further includes: A BVH tree application module, configured to use a BVH tree to perform spatial partitioning on the first mesh units of the rock formation model and the fault plane model respectively, and obtain a plurality of bounding volume structures, where each bounding volume structure includes one or more first mesh units; An interaction analysis module, configured to perform interaction analysis on the rock formation model and the fault plane model based on the plurality of bounding volumes, calculate the intersection points or intersection lines between the rock formation and the fault plane, so as to improve the calculation speed and accuracy of the fault plane cutting the rock formation.

[0020] Optionally, the input module includes: A first acquisition sub-module, configured to acquire geological exploration data, where the geological exploration data includes geometric characteristics of each rock formation; A first processing sub-module, configured to digitally process the geometric shape, physical properties, and spatial distribution of the rock formation according to the geological exploration data by using data interpolation technology and mesh partitioning technology, and obtain the rock formation model; The input module includes: A second acquisition sub-module, configured to acquire geological exploration data, where the geological exploration data includes geometric characteristics of the fault plane; A second processing sub-module, configured to digitally process the geometric shape and position of the fault plane according to the geological exploration data by using data interpolation technology and mesh partitioning technology, and obtain the fault plane model.

[0021] Optionally, the system further includes: A parallel computing module, configured to allocate computing tasks to multiple processing units for parallel computing through multi-threading technology or GPU acceleration technology, so as to support large-scale parallel computing and meet the real-time simulation requirements of complex geological models.

[0022] In a third aspect, an embodiment of the present invention provides an electronic device, including: a processor, a memory, and a computer program stored on the memory and capable of running on the processor. When the computer program is executed by the processor, it implements the steps of a calculation method for segmenting a rock formation by a fault plane.

[0023] The technical solution provided by the present invention at least brings the following beneficial effects: By obtaining a rock formation model and a fault plane model and performing interactive analysis on them, the present invention accurately calculates the intersection points or intersection lines between the rock formation and the fault plane based on this, thereby providing a scientific basis for subsequent rock formation segmentation. Different from traditional two-dimensional segmentation techniques, the three-dimensional intersection point calculation method of the present invention can adapt to more complex curved fault planes and irregular rock formation structures, ensuring high-precision segmentation results. In practical applications, the present invention calculates the offset amount of the rock formation after being segmented by the fault plane by using the physical properties of the rock formation, improving the accuracy of the calculation. In complex geological scenarios, the present invention realizes the efficient processing of the relationship between complex fault planes and rock formations, meeting the requirements of efficient and accurate geological modeling. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or in the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can be obtained based on these drawings.

[0025] Figure 1 It is a schematic diagram of the steps of a calculation method for segmenting a rock formation by a fault plane provided by an embodiment of the present invention; Figure 2 It is a calculation flow chart of the intersection line between a rock formation and a fault plane in an embodiment of the present invention; Figure 3 It is a schematic diagram of the triangular mesh reconstruction method in an embodiment of the present invention; Figure 4 It is a structural block diagram of a calculation system for segmenting a rock formation by a fault plane provided by an embodiment of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0026] The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the drawings in the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art without creative efforts based on the embodiments of the present invention belong to the scope of protection of the present invention.

[0027] With the continuous progress of modern geological exploration technologies, geological modeling is increasingly widely applied in fields such as coal mining, oil and gas exploration, and engineering geology. In these fields, especially for the study of complex geological structures, the precise segmentation and simulation analysis of fault planes and rock formations have become key links. When current technologies deal with the segmentation of rock formations by fault planes, several important problems are often faced: First, when facing complex curved fault planes and irregular rock formations, existing calculation methods are difficult to quickly and accurately complete the segmentation task; Second, due to the dependence of traditional methods on two-dimensional models or simplified assumptions, their processing efficiency is low and cannot meet the requirements of real-time simulation of large-scale geological models and complex scenarios; Finally, after the rock formation segmentation by related technologies, the physical offset of the rock formation under the action of the fault cannot be accurately calculated, resulting in insufficiently accurate prediction results for subsequent applications such as mine pressure analysis, oil and gas exploration, and disaster warning. To address the above problems, the present invention proposes a calculation method for segmenting rock formations by fault planes. By adopting a three-dimensional geometric model, an accurate intersection point calculation method, and a mesh division technology, the present invention can accurately solve the intersection line between the fault plane and the rock formation, and on this basis, perform accurate rock formation segmentation, thereby realizing efficient and accurate calculation of segmenting rock formations by fault planes.

[0028] Figure 1 is a schematic diagram of the steps of a calculation method for segmenting rock formations by a fault plane provided by an embodiment of the present invention, as Figure 1 shown, the method includes: Step S101, obtaining a rock formation model, which is used to describe the geometric characteristics of the rock formation.

[0029] The rock formation model is three-dimensional grid data obtained from geological exploration data. The rock formation model can be created by modeling software (such as CAD, Blender, etc.) and exported in a standard format (such as obj, vtu, etc.). Each grid unit of the rock formation model contains the geometric characteristics, physical characteristics, and other related attributes of the rock formation. Specifically, the geometric characteristics of the rock formation refer to information such as the spatial shape, thickness, and inclination angle of each layer of the rock formation. The physical characteristics of the rock formation include physical parameters of the rock formation, such as density, elastic modulus, and strength. The rock formation model can be three-dimensionally visualized on a computer to help technicians intuitively understand the structure and distribution of the rock formation.

[0030] Step S102, obtaining a fault plane model, which is used to describe the geometric characteristics of the fault plane, and the fault plane includes a curved fault plane.

[0031] The fault plane model represents the spatial position, geometric shape, and its geometric characteristics of the fault plane. Similarly, the fault plane model can also be described by three-dimensional grid data obtained from geological exploration data, and can represent a fault plane with a complex shape. The spatial position of the fault plane is used to determine the position and shape of the fault in three-dimensional space, including information such as the depth and azimuth of the fault plane underground. The geometric shape of the fault plane refers to the shape of the fault plane in three-dimensional space. The geometric characteristics of the fault plane include the inclination angle and displacement direction of the fault plane. In particular, the present invention is applicable to a curved fault plane.

[0032] When dealing with the fault plane, related technologies usually rely on simple plane models or assume that the fault plane is a regular geometric shape, which has great limitations for complex curved fault planes. Traditional methods are difficult to accurately describe the complex geometric characteristics on the fault plane, resulting in the inability to accurately capture the interaction relationship between the curved fault plane and the rock formation when dealing with actual geological situations. Especially when performing rock formation segmentation, the calculation accuracy is low and the efficiency is not high. The subsequent method of the present invention adopts a three-dimensional geometric model and intersection point calculation method, which can accurately describe and process a curved fault plane with a complex shape, uses a geometric cutting algorithm to achieve efficient segmentation, and combines grid division and optimization technologies to ensure high accuracy and high efficiency even in complex fault situations. Therefore, the present invention can overcome the limitations of traditional methods and provide precise processing support for curved fault planes.

[0033] Step S103: Perform an interactive analysis on the rock formation model and the fault plane model, calculate the intersection points or intersection lines of the rock formation and the fault plane, and determine the rock formation model after the fault plane is segmented based on the intersection points or intersection lines of the rock formation and the fault plane.

[0034] Interactive analysis refers to the analysis of the geometric intersection of the rock formation model and the fault plane model. The rock formation model and the fault plane model can be represented by three-dimensional grid data. Through interactive analysis, the intersection relationship between the rock formation and the fault plane is calculated, and the intersection points or intersection lines in three-dimensional space are found. Intersection point calculation is to determine the intersection part between the rock formation model and the fault plane model through three-dimensional geometric methods. In three-dimensional space, the rock formation model and the fault plane may intersect, and the intersection points or intersection lines refer to the specific positions where they intersect. These intersection points or intersection lines constitute the boundary line where the rock formation is cut by the fault plane. In this step, a three-dimensional geometric intersection point calculation method (such as a three-dimensional intersection point solving algorithm) is used to calculate the specific coordinates of these intersection points or intersection lines. By finding the intersection of the fault plane and the rock formation grid model, the boundary line where the fault plane and the rock formation intersect is solved, and this boundary line defines the boundary where the rock formation is segmented by the fault, which will be described in detail later.

[0035] Step S104: Calculate the offset of the rock formation model after the fault plane is segmented according to the geometric characteristics of the fault plane.

[0036] The offset refers to the displacement amount or offset that the rock formation model undergoes under the action of the fault plane, indicating the movement of the rock formation in the direction perpendicular to or parallel to the fault plane. The geometric characteristics of the fault plane (such as the inclination angle of the fault plane and the displacement direction) determine the response of the rock formation under the action of the fault plane. Softer or thinner rock formations may experience larger displacements, while harder or thicker rock formations may have smaller displacements. In the embodiments of the present invention, according to the geometric characteristics of the fault plane, the local offset is calculated, that is, the displacement amount of the rock formation located on the interface under the action of the fault plane is determined. For the entire rock formation model, the individual local offsets are combined to obtain the overall offset of the rock formation.

[0037] Step S105: Output the rock formation model and its offset after the fault plane segmentation for use in strata pressure analysis, oil and gas simulation, and disaster warning.

[0038] The offset obtained in step S104 is output as an independent data field and saved together with the data of the rock formation model to accurately describe the deformation characteristics of the rock formation. For the convenience of subsequent analysis and use, the segmented rock formation model and offset are saved in a standard file format (such as.obj,.vtu), supporting subsequent visualization, simulation, and computational analysis. These results can be used for strata pressure analysis, oil and gas resource simulation, and disaster prediction, etc.

[0039] The present invention obtains the rock formation model and the fault plane model, and conducts interactive analysis on them. Based on this, the intersection points or intersection lines of the rock formation and the fault plane are accurately calculated, thereby providing a scientific basis for subsequent rock formation segmentation. Different from traditional two-dimensional segmentation techniques, the three-dimensional intersection point calculation method of the present invention can adapt to more complex curved fault planes and irregular rock formation structures, ensuring high-precision segmentation results. In practical applications, the present invention calculates the offset of the rock formation after the fault plane segmentation by using the physical properties of the rock formation, improving the calculation accuracy. In complex geological scenarios, the present invention realizes the efficient processing of the relationship between complex fault planes and rock formations, meeting the requirements of efficient and accurate geological modeling.

[0040] In an alternative embodiment, the interactive analysis of the rock formation model and the fault plane model to calculate the intersection points or intersection lines of the rock formation and the fault plane includes: Step S1031: Perform mesh division on the rock formation model and the fault plane model to obtain the respective first mesh units of the rock formation model and the fault plane model.

[0041] The rock stratum model and the fault plane model are meshed and transformed into a discretized format that is more suitable for calculation. Considering that rock strata and fault planes are usually composed of complex geometric shapes and curved surfaces, it is difficult to directly calculate the intersection points or intersection lines. Through meshing, these continuous geometric objects can be transformed into a discrete model composed of a number of relatively simple first grid cells, which is convenient for calculation and analysis. The first grid cell is a small part representing the rock stratum or fault plane model, and it is a grid cell with geometric coordinates and morphological characteristics. After meshing, when calculating the intersection points or intersection lines, it can be deduced through the intersection relationship between grid cells, avoiding directly dealing with complex geometric shapes and improving the calculation speed and efficiency.

[0042] For the rock stratum model, the meshing method can be determined according to the morphological complexity of the geological body and the calculation accuracy requirements. For relatively complex areas, finer grids can be used to improve the calculation accuracy; for relatively simple areas, larger grids can be used to improve the calculation efficiency. Similarly, for the fault plane model, fine meshing can be carried out according to the shape of the fault plane, so as to accurately capture the interaction characteristics between the fault plane and the rock stratum in the intersection point calculation.

[0043] Step S1032, through a three-dimensional geometric intersection algorithm, based on the intersection points between the respective first grid cells of the rock stratum model and the fault plane model, determine the intersection points or intersection lines between the rock stratum and the fault plane.

[0044] As mentioned above, the intersection point is the intersection position in three-dimensional space of the grid cell of the rock stratum model and the grid cell of the fault plane model. The intersection line is the boundary line along a certain path between the rock stratum and the fault plane, and the intersection line describes the boundary where the rock stratum is cut by the fault. The three-dimensional geometric intersection algorithm adopted in the embodiments of the present invention is used to calculate the intersection between two three-dimensional geometric bodies, that is, to calculate the intersection points or intersection lines between the grid cells of the rock stratum and the grid cells of the fault plane. Its core idea is to deduce the position of the intersection point based on a mathematical model by considering the spatial position, shape and size of the two geometric bodies. Specifically, in the rock stratum model and the fault plane model, the calculation of the intersection point is based on the geometric shape of the grid cell and the position relationship of the fault plane. By detecting whether the grid cell of the rock stratum intersects with the grid cell of the fault plane and calculating the precise coordinates of the intersection point. If a continuous path (i.e., the intersection line) is formed between the intersection points, then further calculate the connection method between the intersection points to determine the geometric shape of the intersection line.

[0045] The intersection points or intersection lines are calculated through the three-dimensional geometric intersection algorithm, and these results are used to determine the interaction relationship between the rock stratum and the fault plane and help the subsequent steps to accurately simulate the behavior of the rock stratum under the action of the fault. Figure 2 It is the calculation flow chart of the intersection line between the rock stratum and the fault plane in an embodiment of the present invention. Please refer to Figure 2 , Figure 2The main steps of the shown process include the determination of intersection points, the generation of intersection lines, the dissection and recombination of rock formations. As previously mentioned, through the three-dimensional geometric intersection point algorithm, first, according to the intersection between the rock formation grid cells and the fault plane grid cells, the positions of the intersection points are calculated. If a continuous path is formed between the intersection points, an intersection line will be further formed, which describes the boundary line between the rock formation and the fault plane. After the rock formation is cut by the fault plane, the next step is to process and recombine these divided parts. As Figure 2 shown, when the fault plane intersects the rock formation, the intersection line divides the rock formation into two parts. At this time, the intersection line is equivalent to the "division boundary" of the rock formation. Under the action of the fault plane, the rock formation is cut into two independent blocks. Next, the three-dimensional geometric intersection point algorithm recombines these two parts. Specifically, when two surfaces are recombined, they will form a new rock formation body. Through the action of the fault plane, these two divided surfaces are reconnected to form a new geometric body. For the dissected rock formation, the fault plane is covered on the surface of the original rock formation. Figure 2 The "recombination" shown in it is equivalent to "cutting" the original rock formation, generating two divided parts at the intersection line position. Then, by covering the curved surface of the fault plane on the original rock formation, the fault plane defines a new boundary. At this time, the boundaries of the two divided rock formation blocks will change accordingly according to the shape and position of the fault plane. After the rock formation is cut in half, the middle part will become empty. Because the fault plane cuts the rock formation, the area between the intersection lines no longer has an actual rock formation surface. This part of the void needs to be filled. According to the geometric shape of the fault plane and the geometric relationship between the cut rock formation blocks, the filling method for these voids is calculated. This process can be regarded as filling the "gap", sewing the remaining cut fault part (i.e., the "seam") back onto the rock formation surface. The gap filling is achieved by sewing the remaining fault "seam" onto the surface of the rock formation to fill the void between the fault plane and the rock formation, thereby forming a complete rock formation body. Through this process, the final rock formation body will contain complete geometric surfaces, avoiding the appearance of "hollow" areas after being cut by the fault plane. Through the above steps of intersection point calculation, intersection line generation, dissection, surface recombination, and gap filling, a new rock formation body is finally generated.

[0046] By combining grid generation and geometric intersection point algorithms, complex geological bodies can be effectively and accurately modeled.

[0047] In an alternative embodiment, the determination of the rock formation model after the fault plane division based on the intersection points or intersection lines between the rock formation and the fault plane includes: Step S1033, based on the intersection points or intersection lines between the rock formation and the fault plane, determine the division surface, and use the geometric cutting algorithm to divide the rock formation model to obtain multiple independent rock formation blocks.

[0048] Based on the intersection points or intersection lines, the splitting plane is determined as the plane passing through the intersection points or intersection lines. The splitting plane divides the contact area between the rock formation model and the fault plane model. The splitting plane can be a complex curved surface, usually a region formed by multiple intersection points or intersection lines in three-dimensional space.

[0049] The core idea of the geometric cutting algorithm is to cut the rock formation model according to the intersection points or intersection lines based on the geometric shape of the splitting plane. The geometric cutting algorithm checks each first grid cell of the rock formation model one by one to determine which grid cells have an intersection with the splitting plane and divides the first grid cells into different blocks. Generally, the cutting algorithm uses the intersection points or intersection lines as the cutting boundaries and "cuts" the rock formation model along these boundaries. According to the intersection points of the splitting plane and these grid cells, it is calculated which grid cells are on one side of the splitting plane and which are on the other side. Through the cutting operation, the grid cells that originally belonged to the same rock formation are divided into different rock formation blocks.

[0050] Through the geometric cutting algorithm, the rock formation model is divided into multiple independent rock formation blocks. The rock formation blocks represent the rock formation units under the action of the fault plane, and the position and shape of each rock formation block in space have changed according to the position and shape of the fault plane. The rock formation blocks can be analyzed subsequently through their independent geometric information and provide basic data support for the mine pressure analysis, oil and gas simulation, etc. of the rock formation.

[0051] Step S1034, perform geometric optimization on the rock formation blocks located on the splitting plane to obtain the rock formation model after being split by the fault plane.

[0052] The rock formation model is divided into multiple independent rock formation blocks, but the geometric shapes of these divided rock formation blocks may have irregular, non-smooth boundaries or corners, and may even have distorted or interlaced shapes. These problems may affect the accuracy and stability of subsequent analysis. Therefore, in the embodiments of the present invention, geometric optimization is performed on these rock formation blocks. The optimization goal is to make the geometric shape of each rock formation block more regular and smooth, and reduce unnecessary geometric distortion or intersection. The optimized model will be more suitable for subsequent analysis such as mine pressure, oil and gas simulation, and disaster warning.

[0053] In an optional implementation manner, the performing geometric optimization on the rock formation blocks located on the splitting plane to obtain the rock formation model after being split by the fault plane includes: reconstructing a triangular grid for the splitting points on the splitting plane to obtain a second grid cell corresponding to the rock formation blocks located on the splitting plane; adjusting the second grid cell corresponding to the rock formation blocks located on the splitting plane through an adaptive grid refinement technique to obtain the rock formation model after being split by the fault plane.

[0054] Around the segmentation points on the segmentation surface, the grid is reconstructed to improve the accuracy. That is, the local accuracy of the grid is enhanced in the area where the fault plane intersects with the rock formation. In the embodiments of the present invention, triangular grids are used, which can better adapt to the curved surface morphology. The purpose of reconstructing the grid is to make the connection between the segmentation points and the edges and nodes of the adjacent grids closer by refining and optimizing the grid distribution, so as to more accurately describe the intersection of the rock mass and the fault plane.

[0055] Figure 3 It is a schematic diagram of the reconstruction method of triangular grids in an embodiment of the present invention. Figure 3 It shows that the triangular grid division demonstrates different cutting situations, and the redivision boundaries are marked by dotted lines. Specifically, Figure 3 In it, the black triangles represent the original triangular grid units, the white triangles represent the adjacent triangular grid units, the dark lines represent the cutting paths, and the white dotted lines represent the redivision cutting paths. If there are no adjacent triangles intersecting around the black triangles (that is, the grids near the segmentation points are not involved), these triangles do not need to be redivided. Therefore, only the triangular grids that intersect with the segmentation points of the adjacent triangles need to be reconstructed. When reconstructing the triangular grids, there are six situations as Figure 3 shown. In the first situation (the first one in the first row), the cutting path passes through the interior of the black triangle. At this time, the fault plane just passes through the middle of the black triangle, and the black triangle is directly cut along the white dotted line into five new triangular grids. The second situation (the second one in the first row) and the third situation (the third one in the first row) are the same as the first situation. In the fourth situation (the first one in the second row), the cutting path passes through the vertex of the black triangle and extends into the interior of the black triangle, but does not completely pass through the entire interior of the triangle. At this time, the fault plane only passes through one vertex of the triangle, and the cutting does not completely pass through the black triangle, so the black triangle is redivided into three new triangular grids. In the fifth situation (the second one in the second row), the cutting path only coincides with the upper half of the side of the black triangle. At this time, the fault plane does not completely cut the triangle, so the adjacent white triangle is redivided along the white dotted line. In the sixth situation (the third one in the second row), the cutting path coincides with the side of the black triangle. In this case, the fault plane cuts along the side of the triangle, and no processing needs to be done to the adjacent triangles.

[0056] After reconstructing the triangular grids, the corresponding second grid units are obtained. The second grid units are more refined and better adapt to the position and morphology of the segmentation surface than the original first grid units. It should be noted that the first grid units and the second grid units only exist as a distinction form in different grid division stages, and do not represent different types of grid units. Essentially, they are all obtained by dividing the model into grids.

[0057] Through the above process, the second grid cells can provide a more accurate geometric approximation locally, enabling the geometric relationship at the junction of the fault plane and the rock formation to be represented more precisely. This is the first step in geometric optimization, which helps create a more suitable grid structure.

[0058] Adaptive Mesh Refinement (AMR) is a technique that dynamically adjusts the grid density according to the needs of specific regions. Different from the traditional uniform grid division method, the adaptive mesh refinement technique refines in regions that require high precision (such as near the segmentation plane or the fault plane region), while maintaining a coarser grid in other regions to improve computational efficiency. In the rock formation blocks after fault plane segmentation, especially at the boundaries in contact with the fault plane, a higher resolution is required to capture the complex geological morphology. In the embodiments of the present invention, the adaptive mesh refinement technique is used to adjust the size of the second grid cells according to the needs of these specific regions, making the second grid cells near the fault plane more detailed, thereby improving the local geometric accuracy. The refinement technique dynamically adjusts the grid density according to the deformation degree of the grid cells, geometric complexity, or numerical accuracy requirements of the solution. For example, in the region near the junction of the fault plane and the rock formation, the second grid cells are denser, so that the local deformation of the fault and the changes in the rock formation can be better simulated.

[0059] Specifically, through the adaptive mesh refinement technique, the regions that require higher precision are first identified, and then smaller second grid cells are generated within these regions. The number of nodes and edges of these newly generated second grid cells will increase significantly, thereby enhancing the accuracy of the local geometric morphology. Through this process, the accuracy of the rock formation model is improved, while still maintaining the efficiency of the overall grid, avoiding unnecessary waste of computing resources.

[0060] After geometric optimization, including reconstructing the grid and adaptive mesh refinement, the finally obtained rock formation model will have a high-precision local structure. In particular, the geometric shape of the rock formation blocks near the segmentation plane will be more refined. The optimized rock formation model will accurately reflect the role of the fault plane, showing the interaction between the rock formation and the fault plane, ensuring that accurate results can be obtained in subsequent geological analyses (such as oil and gas simulation, mine pressure analysis, etc.).

[0061] In an optional embodiment, calculating the offset of the rock formation model after fault plane segmentation according to the geometric characteristics of the fault plane includes: calculating each rock formation block after fault plane segmentation one by one through a geometric transformation algorithm according to a preset given value or the geometric characteristics of the fault plane, to obtain the offset of each rock formation block after fault plane segmentation.

[0062] In the embodiments of the present invention, according to a given value (the given value can be a specific displacement or angle, which is derived from existing assumptions in previous geological surveys or engineering designs), or according to the physical properties of the fault plane (such as the inclination angle and fault plane displacement), the offset of the segmented rock layer model is calculated. The geometric transformation algorithm adopted in the embodiments of the present invention realizes the precise offset of the rock layer through geometric transformation. Under the action of the fault, the rock layer will be offset to varying degrees according to its physical properties. That is to say, the geometric properties of the fault plane determine the offset of each segmented rock layer block. Through the geometric transformation algorithm, according to the geometric properties of the fault plane, the displacement that occurs during the deformation of the rock layer is calculated, and the displacement can be a combination of deformation methods such as translation, rotation, and shear.

[0063] Specifically, each grid unit of the rock layer block is regarded as an independent geometric unit through the geometric transformation algorithm, and according to its position and relative coordinate system, the initial state of these grid units in three-dimensional space is determined. Taking the physical properties of the fault plane as input data, using the geometric transformation algorithm, the offset of each grid unit of each rock layer block under the action of the fault is calculated block by block.

[0064] The geometric transformation algorithm can accurately handle the interaction relationship between complex fault planes and irregular rock layers, has high calculation accuracy and efficiency, and is applicable to fields such as oil and gas exploration, geological modeling, and underground engineering design. Through the geometric transformation algorithm, the accuracy of fault segmentation can be effectively improved, and the reliability of geological modeling and resource assessment can be enhanced.

[0065] In an alternative embodiment, the method further includes: Step S201, using a BVH tree to perform spatial partitioning on the first grid units of the rock layer model and the fault plane model respectively, obtaining a plurality of bounding volume structures, and each bounding volume structure includes one or more first grid units.

[0066] The present invention proposes that when dealing with the interaction between the fault plane model and the rock layer model, a BVH tree data structure is used for spatial partitioning to reduce the complexity and scope of calculation, thereby improving efficiency. The BVH tree (Bounding Volume Hierarchy) is a tree-shaped data structure used to accelerate collision detection, interaction analysis, and geometric calculation in space. The principle of the BVH tree is to enclose the geometric bodies in the scene in simple geometric shapes (such as boxes, spheres, or cylinders), and these bounding volumes are used to represent more complex geometric shapes. By dividing the space of the object into a tree structure, the number of geometric units that need to be calculated in detail can be significantly reduced, thereby improving the calculation efficiency.

[0067] In an embodiment of the present invention, the role of the BVH tree is to perform spatial partitioning on the first grid cells of the rock formation model and the fault plane model. By assigning one or more bounding volumes to each grid cell, the interactive analysis between the models is accelerated, especially when finding the intersection between the fault plane and the rock formation model.

[0068] To accelerate the interactive analysis, first, the first grid cells of the rock formation model and the fault plane model are enclosed within larger geometric shapes (bounding volumes). The bounding volume is a simple geometric shape, such as a rectangular box (AABB, axis-aligned bounding box) or a sphere, etc. The role of the bounding volume is to reduce the number of first grid cells that need to be calculated by converting complex geometric bodies into simple geometric shapes.

[0069] Step S202: Based on the multiple bounding volumes, perform an interactive analysis on the rock formation model and the fault plane model, and calculate the intersection points or intersection lines between the rock formation and the fault plane to improve the calculation speed and accuracy of the fault plane cutting the rock formation.

[0070] When the fault plane interacts with the rock formation model, first check whether these two bounding volumes intersect. If they do not intersect, then these first grid cells do not need to perform further precise interactive calculations. Only when the two bounding volumes intersect will more detailed collision or intersection point calculations be performed. In this way, by assigning bounding volumes to each first grid cell, a large amount of computational effort can be reduced from the detailed level to a higher level. If the bounding volumes intersect, it means that there may be intersection points or intersection lines in these two regions. On this basis, further refine the calculation to check whether there are actual intersection points or intersection lines for these grid cells. The detailed calculation is performed on the grid cells within the intersecting bounding volumes, and the calculation steps for further interactive analysis are as described above.

[0071] In an alternative embodiment, the obtaining of the rock formation model includes: obtaining geological exploration data, where the geological exploration data includes the geometric characteristics of each rock formation; according to the geological exploration data, using data interpolation technology and grid division technology, digitize the geometric shape, physical properties, and spatial distribution of the rock formation to obtain the rock formation model.

[0072] Geological exploration data is the basis for constructing the rock formation model. It can come from geological surveys, drilling, seismic exploration, etc., and contains the geometric characteristics of the underground rock formations, specifically including information such as the spatial form, thickness, and inclination angle of each rock formation layer. Using data interpolation technology to convert these geological exploration data into a digital rock formation model, the data interpolation technology is used to process and fill in the gaps in the geological exploration data to generate a complete model.

[0073] Data interpolation technology refers to inferring the values in unknown areas based on known data points. For example, if the physical properties or geometric shape data of certain rock formations are known, data interpolation technology can help estimate the values of these data in other areas. The goal of data interpolation technology is to infer relevant information at other unmeasured positions in space through existing discrete data. Interpolation methods can include common interpolation methods such as Lagrange interpolation method, spline interpolation method, and Kriging interpolation method, which can infer the changes of rock formations in the entire space based on geological exploration data points.

[0074] Mesh generation technology discretizes the data in space into a finite number of mesh elements to facilitate subsequent calculations and analyses. This has been described in detail previously and will not be elaborated here.

[0075] The obtaining of the fault plane model includes: obtaining geological exploration data, where the geological exploration data includes the geometric characteristics of the fault plane; based on the geological exploration data, using data interpolation technology and mesh generation technology, digitizing the geometric shape and position of the fault plane to obtain the fault plane model.

[0076] Similarly, geological exploration data is the basis for constructing the fault plane model, which can be obtained through geological surveys, drilling, seismic exploration, etc., and contains the geometric characteristics of the fault plane, specifically including the dip angle and displacement direction of the fault plane. Using data interpolation technology, these geological exploration data are transformed into a digital fault plane model. Data interpolation technology is used to process and fill in the blanks in geological exploration data to generate a complete fault plane model.

[0077] Mesh generation technology discretizes the data in space into a finite number of mesh elements to facilitate subsequent calculations and analyses. This has been described in detail previously and will not be elaborated here.

[0078] Generally speaking, through the processing of these geological exploration data, using data interpolation technology and mesh generation technology, discrete exploration data can be transformed into a complete and continuous digital model of rock formations and fault planes.

[0079] In an alternative embodiment, the method further includes: through multi-threading technology or GPU acceleration technology, distributing computing tasks to multiple processing units for parallel computing to support large-scale parallel computing and meet the real-time simulation requirements of complex geological models.

[0080] Multithreading technology refers to dividing computational tasks into multiple independent subtasks, which can run simultaneously on multiple processor cores, leveraging the advantages of multi-core processors to accelerate the computational process. In the case of processing a large amount of geological data and performing complex calculations in the embodiments of the present invention (such as the interactive analysis of fault planes and rock formations, the segmentation and offset calculation of models, etc.), by dividing these computational tasks into multiple parallel subtasks, the computational time can be significantly reduced. For example, for the calculation of a three-dimensional geological model, traditional single-threaded calculation needs to process each data point step by step. With multithreading technology, these data points can be simultaneously allocated to multiple threads for parallel processing, greatly enhancing the computational speed and efficiency.

[0081] GPU (Graphics Processing Unit) acceleration technology utilizes a large number of parallel processing units in the graphics card to accelerate computational tasks. Compared with traditional CPUs, GPUs have more computational cores and can parallelly process a large number of simple computational tasks. This makes GPUs more advantageous than traditional CPUs when dealing with tasks such as large-scale computing, graphics rendering, and scientific computing.

[0082] In specific implementation, for subtasks with small tasks and less computational volume, they can be allocated to multiple CPU cores for parallel execution through multithreading technology. For tasks with large computational volume and high parallelizability, leveraging the powerful parallel computing ability of GPUs can greatly accelerate the computational process. For example, when performing rock formation segmentation and fault plane offset calculation, multiple computational tasks can run simultaneously in multiple processing units of the GPU.

[0083] By combining multithreading technology and GPU acceleration technology, the computational tasks of complex geological models can be parallelized, making full use of the computational power of modern hardware. It can achieve fast geological model simulation in a large-scale parallel computing environment, significantly improve computational efficiency, meet real-time simulation requirements, and provide strong support for complex geological analysis and decision-making.

[0084] Figure 4 It is a structural block diagram of a computational system for segmenting rock formations by a fault plane provided by an embodiment of the present invention, as Figure 4 shown, the system includes: A first input module 301, configured to obtain a rock formation model, where the rock formation model is used to describe the geometric characteristics of the rock formation; A second input module 302, configured to obtain a fault plane model, where the fault plane model is used to describe the geometric characteristics of the fault plane; An intersection module 303, configured to perform interactive analysis on the rock formation model and the fault plane model, calculate the intersection points or intersection lines of the rock formation and the fault plane, and determine the rock formation model after segmentation by the fault plane based on the intersection points or intersection lines of the rock formation and the fault plane; An offset calculation module 304 is configured to calculate the offset of the rock formation model after the fault plane is segmented according to the geometric characteristics of the fault plane; An output module 305 is configured to output the rock formation model after the fault plane is segmented and its offset for use in mine pressure analysis, oil and gas simulation, and disaster warning.

[0085] In an optional implementation manner, the intersection module includes: A first mesh division module is configured to perform mesh division on the rock formation model and the fault plane model to obtain respective first mesh units of the rock formation model and the fault plane model; An intersection sub-module is configured to determine the intersection points or intersection lines between the rock formation and the fault plane based on the intersection points between the respective first mesh units of the rock formation model and the fault plane model through a three-dimensional geometric intersection algorithm.

[0086] In an optional implementation manner, the intersection module includes: A dividing surface determination sub-module is configured to determine a dividing surface according to the intersection points or intersection lines between the rock formation and the fault plane, and use a geometric cutting algorithm to segment the rock formation model to obtain multiple independent rock formation blocks; A geometric optimization sub-module is configured to perform geometric optimization on the rock formation blocks located on the dividing surface to obtain the rock formation model after the fault plane is segmented.

[0087] In an optional implementation manner, the geometric optimization sub-module includes: A construction unit is configured to reconstruct a triangular mesh for the segmentation points on the dividing surface to obtain respective second mesh units corresponding to the rock formation blocks located on the dividing surface; An adjustment unit is configured to adjust the respective second mesh units corresponding to the rock formation blocks located on the dividing surface through an adaptive mesh refinement technique to obtain the rock formation model after the fault plane is segmented.

[0088] In an optional implementation manner, the offset calculation module includes: A block-by-block calculation sub-module is configured to perform block-by-block calculation on the rock formation model after the fault plane is segmented according to a preset given value or the geometric characteristics of the fault plane through a geometric transformation algorithm to obtain the offset of each rock formation block after the fault plane is segmented.

[0089] In an optional implementation manner, the system further includes: A BVH tree application module is configured to perform spatial division on the respective first mesh units of the rock formation model and the fault plane model by using a BVH tree to obtain a plurality of bounding volume structures, and each bounding volume structure includes one or more first mesh units; An interaction analysis module, configured to perform interaction analysis on the rock formation model and the fault plane model based on the multiple bounding volumes, and calculate the intersection points or intersection lines between the rock formation and the fault plane, so as to improve the calculation speed and accuracy of the fault plane cutting the rock formation.

[0090] In an alternative embodiment, the input module includes: A first acquisition sub-module, configured to acquire geological exploration data, where the geological exploration data includes the geometric characteristics of each rock formation; A first processing sub-module, configured to digitally process the geometric shape, physical properties and their spatial distribution of the rock formation according to the geological exploration data by using data interpolation technology and grid division technology, so as to obtain the rock formation model; The input module includes: A second acquisition sub-module, configured to acquire geological exploration data, where the geological exploration data includes the geometric characteristics of the fault plane; A second processing sub-module, configured to digitally process the geometric shape and position of the fault plane according to the geological exploration data by using data interpolation technology and grid division technology, so as to obtain the fault plane model.

[0091] In an alternative embodiment, the system further includes: A parallel computing module, configured to distribute computing tasks to multiple processing units for parallel computing through multi-threading technology or GPU acceleration technology, so as to support large-scale parallel computing and meet the real-time simulation requirements of complex geological models.

[0092] The embodiments of the present disclosure also provide an electronic device, including a processor, a memory, and a computer program stored on the memory and capable of running on the processor. When the computer program is executed by the processor, it implements each process of the above-mentioned embodiment of the calculation method for a fault plane to divide a rock formation, and can achieve the same technical effect. To avoid repetition, it will not be elaborated here.

[0093] Those skilled in the art should understand that the embodiments of the present invention can be provided as methods, devices, electronic devices, and storage media. Therefore, the embodiments of the present invention can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the embodiments of the present invention can take the form of a computer program product implemented on one or more computer-readable storage media (including but not limited to disk memories, CD-ROMs, optical memories, etc.) containing computer-usable program codes.

[0094] Embodiments of the present invention are described with reference to the flowcharts and / or block diagrams of methods and apparatuses according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing terminal devices to generate a machine, such that the instructions executed by the processor of the computer or other programmable data processing terminal devices generate a device for implementing the functions specified in one process Figure 1 one process or multiple processes and / or blocks Figure 1 one block or multiple blocks. These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing terminal device to work in a specific manner, such that the instructions stored in the computer-readable memory generate a manufactured article including an instruction device, and the instruction device implements the functions specified in one process Figure 1 one process or multiple processes and / or blocks Figure 1 one block or multiple blocks. These computer program instructions can also be loaded onto a computer or other programmable data processing terminal device, such that a series of operation steps are executed on the computer or other programmable terminal device to generate a computer-implemented process. Thus, the instructions executed on the computer or other programmable terminal device provide steps for implementing the functions specified in one process Figure 1 one process or multiple processes and / or blocks Figure 1 one block or multiple blocks.

[0095] Although the preferred embodiments of the embodiments of the present invention have been described, those skilled in the art can make additional changes and modifications once they learn the basic creative concepts. Therefore, the appended claims are intended to be construed to include the preferred embodiments and all changes and modifications falling within the scope of the embodiments of the present invention.

[0096] Finally, it should also be noted that in this text, relational terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the term "including" or any other variant thereof is intended to cover non-exclusive inclusion, such that a process, method, article, or terminal device including a series of elements not only includes those elements, but also includes other elements not explicitly listed, or elements inherent to such process, method, article, or terminal device. Without further limitation, the elements defined by the statement "including..." do not exclude the existence of additional identical elements in the process, method, article, or terminal device including the said elements.

[0097] The above has introduced in detail a calculation method, system and electronic device for dividing rock formations by fault planes. In this article, specific examples are used to elaborate on the principle and implementation manner of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea; at the same time, for those of ordinary skill in the art, according to the idea of the present invention, there will be changes in the specific implementation manner and application scope. In summary, the content of this specification should not be construed as a limitation to the present invention.

Claims

1. A calculation method for dividing rock formations by fault planes, characterized in that The method includes: Obtaining a rock formation model for describing the geometric characteristics of the rock formation; Obtaining a fault plane model for describing the geometric characteristics of the fault plane, where the fault plane includes a curved fault plane; Performing interactive analysis on the rock formation model and the fault plane model, calculating the intersection points or intersection lines between the rock formation and the fault plane, and determining the rock formation model after being segmented by the fault plane based on the intersection points or intersection lines between the rock formation and the fault plane; Calculating the offset of the rock formation model after being segmented by the fault plane according to the geometric characteristics of the fault plane; Outputting the rock formation model after being segmented by the fault plane and its offset for use in mine pressure analysis, oil and gas simulation, and disaster warning.

2. The method according to claim 1, wherein The performing interactive analysis on the rock formation model and the fault plane model, calculating the intersection points or intersection lines between the rock formation and the fault plane, includes: Performing mesh division on the rock formation model and the fault plane model to obtain the respective first mesh units of the rock formation model and the fault plane model; Based on the intersection points between the respective first mesh units of the rock formation model and the fault plane model, determining the intersection points or intersection lines between the rock formation and the fault plane through a three-dimensional geometric intersection algorithm.

3. The method according to claim 1, characterized in that, The determining the rock formation model after being segmented by the fault plane based on the intersection points or intersection lines between the rock formation and the fault plane includes: Determining a segmentation plane according to the intersection points or intersection lines between the rock formation and the fault plane, and using a geometric cutting algorithm to segment the rock formation model to obtain multiple independent rock formation blocks; Performing geometric optimization on the rock formation blocks located on the segmentation plane to obtain the rock formation model after being segmented by the fault plane.

4. The method according to claim 3, wherein The performing geometric optimization on the rock formation blocks located on the segmentation plane to obtain the rock formation model after being segmented by the fault plane includes: For the segmentation points on the segmentation plane, reconstructing a triangular mesh to obtain the corresponding second mesh units of the rock formation blocks located on the segmentation plane; Adjusting the corresponding second mesh units of the rock formation blocks located on the segmentation plane through an adaptive mesh refinement technique to obtain the rock formation model after being segmented by the fault plane.

5. The method according to claim 3, wherein The calculating the offset of the rock formation model after being segmented by the fault plane according to the geometric characteristics of the fault plane includes: According to a preset given value or the geometric characteristics of the fault plane, performing block-by-block calculation on the rock formation model after being segmented by the fault plane through a geometric transformation algorithm to obtain the offset of each rock formation block after being segmented by the fault plane.

6. The method according to claim 2, wherein The method further includes: Using a BVH tree to perform spatial division on the respective first mesh units of the rock formation model and the fault plane model to obtain a plurality of bounding volume structures, where each bounding volume structure includes one or more first mesh units; Based on the plurality of bounding volumes, performing interactive analysis on the rock formation model and the fault plane model, calculating the intersection points or intersection lines between the rock formation and the fault plane to improve the calculation speed and accuracy of the fault plane cutting the rock formation.

7. The method according to claim 1, characterized in that The obtaining the rock formation model includes: Obtaining geological exploration data, where the geological exploration data includes the geometric characteristics of each rock formation; According to the geological exploration data, using data interpolation technology and mesh division technology to digitally process the geometric shape, physical properties, and spatial distribution of the rock formation to obtain the rock formation model; The obtaining the fault plane model includes: Obtain geological exploration data, where the geological exploration data includes the geometric characteristics of the fault plane; According to the geological exploration data, use data interpolation technology and grid division technology to digitize the geometric shape and position of the fault plane to obtain the fault plane model.

8. The method according to any one of claims 1-7, characterized in that, The method further includes: Through multi-threading technology or GPU acceleration technology, allocate computing tasks to multiple processing units for parallel computing to support large-scale parallel computing and meet the real-time simulation requirements of complex geological models.

9. A calculation system for dividing rock formations by fault planes, characterized in that, The system includes: A first input module for obtaining a rock formation model, where the rock formation model is used to describe the geometric characteristics of the rock formation; A second input module for obtaining a fault plane model, where the fault plane model is used to describe the geometric characteristics of the fault plane; An intersection module for performing interactive analysis on the rock formation model and the fault plane model, calculating the intersection points or intersection lines of the rock formation and the fault plane, and determining the rock formation model after the fault plane segmentation based on the intersection points or intersection lines of the rock formation and the fault plane; An offset calculation module for calculating the offset of the rock formation model after the fault plane segmentation according to the geometric characteristics of the fault plane; An output module for outputting the rock formation model after the fault plane segmentation and its offset for use in mine pressure analysis, oil and gas simulation, and disaster warning.

10. An electronic device, characterized in that, Includes: A processor, a memory, and a computer program stored on the memory and capable of running on the processor. When the computer program is executed by the processor, it implements the steps of the method according to any one of claims 1-8.

Images