Modeling method and analysis method of finite element model of geological body with three-dimensional fine structure of fault

Abstract

The application discloses a kind of geological body finite element model modeling method and analysis method containing three-dimensional fine structure of fault, and relates to the field of earth science.The application is handled by being appropriately smoothed to three-dimensional fault model, not only retains the three-dimensional geometric feature of fault, but also is favorable to subsequent grid division.Meanwhile in grid division, even if it is more complex geometric form, still adopts hexahedron unit, is favorable to more accurate calculation simulation, and finally obtains geological body finite element model containing three-dimensional fine structure of fault.The application can be flexibly adjusted according to different geological background and scientific problem, including but not limited to adding corresponding material medium, friction coefficient, physical equation and fault contact and the like parameter, and geological body and fault of research area can be more accurately simulated.Meanwhile, the grid of geological body finite element model created by the application can reach million level, is favorable to more fine research and simulation and processing work.

Description

Technical Field

[0001] This invention relates to the field of Earth sciences, specifically to a finite element modeling method and analysis method for geological bodies with three-dimensional fine structures containing faults. Background Technology

[0002] The fault geometries in the three-dimensional finite element geological models widely created in current research are often quite simple, with sparse meshes, a small number of meshes, and large mesh element sizes, significantly simplifying the actual geological conditions. Therefore, it is difficult to conduct further high-precision studies based on the fine three-dimensional geometric structure of faults.

[0003] In the process of creating a 3D finite element model, the complex boundary geometry between blocks, the embedding of fault structures, and mesh generation are all key and challenging aspects. Particularly challenging is the geometric division of blocks and the subsequent mesh generation for a detailed 3D fault geometry model. Summary of the Invention

[0004] To address the aforementioned shortcomings in existing technologies, this invention provides a finite element modeling and analysis method for geological bodies containing three-dimensional fine fault structures. This method solves the problems that existing finite element models of geological bodies struggle to preserve the three-dimensional fine fault structures and to conduct further high-precision stress-strain evolution analysis on faults with complex geometries.

[0005] To achieve the above-mentioned objectives, the technical solution adopted by this invention is as follows:

[0006] A finite element modeling method for geological bodies with detailed three-dimensional fault structures is provided, which includes the following steps:

[0007] Acquire geological and geophysical data, and construct three-dimensional fault models using geological three-dimensional modeling software;

[0008] Import the original three-dimensional fault surface from the three-dimensional fault model into Rhino software, and smooth the original three-dimensional fault surface using Rhino software to obtain the smoothed three-dimensional fault surface.

[0009] The smoothed fault 3D surface is imported into HyperMesh software for geometric partitioning of blocks, resulting in a geometric model containing several geometric bodies.

[0010] In HyperMesh software, the geometric model is meshed to obtain a three-dimensional mesh model;

[0011] In HyperMesh software, the 3D mesh model is divided into blocks based on the geometric model to obtain the divided 3D mesh model blocks.

[0012] Enter the contact module, create contact surfaces, select the segmented 3D mesh model block, and select the mesh surface according to the location and depth of the fault plane to obtain the finite element model fault plane; for each finite element model fault, select two surfaces on both sides of the corresponding block to form contact surfaces; retain fault element bands of a set thickness on both sides of the fault interface.

[0013] All three-dimensional mesh model blocks are combined into a geological finite element model, and the relative motion of the two disks of the fault layer in the geological finite element model is simulated by setting different contact relationships on the contact surfaces.

[0014] A finite element modeling method for geological bodies with three-dimensional fine-grained fault structures, comprising the following steps:

[0015] Import the geological finite element model .inp file from Hypermesh software into Abaqus software;

[0016] In the Abaqus software, based on the actual research situation, different layers and blocks are assigned corresponding material medium and elastic parameters, and contact conditions and boundary conditions are set.

[0017] Create time steps, select static or dynamic time steps as needed, customize the step size and name of each time step, set the output field according to requirements, submit analysis and calculation, and obtain the generated odb file;

[0018] Open the generated odb file, create a displaygroup, extract the fault plane, visualize the results on the fault plane, and show the dynamic changes of the physical field on the fault plane over time.

[0019] The beneficial effects of this invention are as follows:

[0020] This invention, by appropriately smoothing the three-dimensional fault model, not only preserves the three-dimensional geometric features of the fault but also facilitates subsequent mesh generation. Furthermore, even with relatively complex geometries, hexahedral elements are used during mesh generation, which is beneficial for more accurate computational simulation, ultimately resulting in a finite element model of a geological body containing the fine three-dimensional structure of the fault.

[0021] This invention can be flexibly adjusted according to different geological backgrounds and scientific problems, including but not limited to adding appropriate material media, friction coefficients, physical equations, and fault contact parameters, enabling more accurate and detailed simulation of geological bodies and three-dimensional faults in the study area. Furthermore, the finite element model mesh of the geological body created by this invention can reach the million-level, facilitating more precise research, simulation, and processing. Attached Figure Description

[0022] Figure 1 A flowchart illustrating the finite element modeling method for this geological body with a detailed three-dimensional fault structure;

[0023] Figure 2 This is a schematic diagram illustrating the data interpretation and information extraction operations in Example 1;

[0024] Figure 3 This is a schematic diagram illustrating the process of exporting the constructed three-dimensional fault model as a DXF format data file in Example 1.

[0025] Figure 4 This is a schematic diagram illustrating the operation of selecting entity editing in Example 1;

[0026] Figure 5 This is a schematic diagram of the operation of entering the contact module in step S6 of Example 1;

[0027] Figure 6 This is a schematic diagram of the operation of creating contact surfaces in step S6 of Example 1;

[0028] Figure 7 This is a schematic diagram of the contact surface formed in step S6 of Embodiment 1;

[0029] Figure 8 This is a schematic diagram illustrating the operation of assigning corresponding material mediums and elastic parameters to different layers and blocks in Example 2;

[0030] Figure 9 This is a schematic diagram illustrating the operation of setting contact conditions in Example 2;

[0031] Figure 10 This is a schematic diagram illustrating the operation of setting boundary conditions in Example 2;

[0032] Figure 11 This is a schematic diagram illustrating the operation of creating a time step in Example 2;

[0033] Figure 12 This is a schematic diagram illustrating the operation of setting the output field in Example 2;

[0034] Figure 13 This is a schematic diagram of the strain rate simulation results in Example 2;

[0035] Figure 14 This is a schematic diagram of setting up a displaygroup in Example 2;

[0036] Figure 15 This is a schematic diagram of the simulation calculation results at the interruption level in Example 2;

[0037] Figure 16This is a vector diagram showing the simulation results on the interrupted plane in Example 2;

[0038] Figure 17 This is a schematic diagram of the three-dimensional fine structure model of the Longmenshan Fault in the embodiment. Detailed Implementation

[0039] The specific embodiments of the present invention are described below to enable those skilled in the art to understand the present invention. However, it should be understood that the present invention is not limited to the scope of the specific embodiments. For those skilled in the art, various changes are obvious as long as they are within the spirit and scope of the present invention as defined and determined by the appended claims. All inventions utilizing the concept of the present invention are protected.

[0040] Example 1:

[0041] like Figure 1 As shown, the finite element modeling method for this geological body with a detailed three-dimensional fault structure includes the following steps:

[0042] S1. Obtain geological and geophysical data, and construct a three-dimensional fault model using geological three-dimensional modeling software;

[0043] S2. Import the original three-dimensional fault surface from the three-dimensional fault model into Rhino software, and smooth the original three-dimensional fault surface using Rhino software to obtain the smoothed three-dimensional fault surface.

[0044] S3. Import the smoothed fault 3D surface into HyperMesh software to perform geometric division of blocks, and obtain a geometric model containing several geometric bodies.

[0045] S4. Mesh the geometric model in HyperMesh software to obtain a three-dimensional mesh model;

[0046] S5. In HyperMesh software, divide the 3D mesh model into blocks based on the geometric model and obtain the divided 3D mesh model blocks.

[0047] S6. Enter the contact module, create contact surfaces, select the divided 3D mesh model blocks, and select the mesh surfaces according to the location and depth of the fault plane to obtain the finite element model fault plane; for each finite element model fault, select two surfaces on both sides of the corresponding block to form contact surfaces; the fault interface retains fault element bands of a set thickness on both sides.

[0048] S7. Assemble all the three-dimensional mesh model blocks into a geological finite element model, and simulate the relative motion of the two disks of the fault layer in the geological finite element model by setting different contact relationships on the contact surfaces.

[0049] In this embodiment, acquiring geological and geophysical data requires the collection of geological and seismic data as well as geophysical data. Data integration and loading are also necessary, unifying various types of exploration data, observation results, and maps from different depths from the surface to the subsurface into the same coordinate system and loading them into 3D modeling software or platforms, such as SKUA-GoCad, Petrel, GOCAD, etc. The following example uses SKUA-GoCad:

[0050] Data interpretation and information extraction: This involves interpreting various types of loaded data to identify and extract the spatial location (point, line, or surface form) and attribute information of active faults. For example... Figure 2 As shown, select "import". Depending on the data format, if points are in CSV format, you can choose "column-based file" to import. Lines and surfaces are generally in SHP file format, so select "ArcViewShp" to import the file. Then select the "surface" function. Depending on the imported points (small earthquake data), you can select "PointSet" to create a surface from points, or "Curve" to create a surface from multiple imported fault profile lines. This will allow you to initially construct the fault plane.

[0051] Point features mainly include fault points measured in field surveys, fault points exposed by trench excavation, fault points exposed by boreholes, or fault points interpreted through borehole combined profiles, and upper fault points interpreted by artificial seismic exploration profiles; linear features mainly include fault lines obtained from active fault mapping, fault lines interpreted through borehole combined profiles, fault lines interpreted through artificial seismic exploration profiles, fault lines interpreted based on the relocation results of strong earthquake aftershocks, fault surface traces obtained by InSAR inversion, and fault lines interpreted by non-seismic geophysical exploration; areal features mainly include fault planes interpreted by 3D seismic exploration, fault nodules obtained from strong earthquake focal mechanism solutions, seismogenic fault planes obtained by inversion from coseismic deformation observation data, and fault planes obtained by joint inversion using 3D geophysical methods.

[0052] Based on the preliminary fault planes established above, and based on geological patterns and spatial topological relationships, appropriate interpolation algorithms are selected to spatially interpolate the interpreted fault point, line, and surface data, generating continuous and overall smooth fault planes, thus constructing a three-dimensional fault model. For example... Figure 3 As shown, the final constructed 3D fault model is exported as a DXF format data file.

[0053] In this embodiment, the specific method of step S2 includes the following steps:

[0054] (1) Import the original three-dimensional fault surface in the three-dimensional fault model into Rhino, use RhinoResurf to discretize the surface into a set of points, export the set of points and reverse the Z-axis, because in general the fault surface is positive with the Z-axis downward, while in Rhino the Z-axis is positive with the Z-axis upward.

[0055] (2) Import the inverted point set into Rhino, import the adjusted fault points into Rhino, adjust the fault points, form a surface based on the point cloud, and adjust the control points and order to make the surface smoother.

[0056] (3) Generate a mesh from points, extract the mesh border, and use the Trim command to trim the surface with the border, thereby generating a smooth three-dimensional fault surface.

[0057] In this embodiment, the specific method of step S3 includes the following steps:

[0058] (1) Create a plane of the study area in Rhino and create different planes according to the fault incision depth. Export the plane and fault plane in .igs format from Rhino and import them into HyperMesh.

[0059] (2) In HyperMesh, select the plane with Z=0, drag it into a solid, and select the depth required for the study.

[0060] (3) For geological bodies that need to be cut according to the fault depth, select solid edit and select the trim with plane / surf function. Select the entity and surface to be cut in sequence. Use the geometric plane with the maximum fault depth to cut the geological body at a depth. Cut the geological body again through the fault plane at this depth to complete a complete cutting operation. Select the second deepest geometric plane of the last complete cutting operation and repeat the complete cutting operation process to obtain several geometric bodies.

[0061] (4) If it is not necessary to cut according to the fault depth, you can cut the shallow layer (e.g., 20km) first, and then select trim withlines to cut the remaining solid using the cutting traces on the 20km surface after the previous step.

[0062] (5) Select trim with plane / surf to cut with planes of different depths;

[0063] (6) At this point, the construction of the geometry is completed, and the horizontal block division and vertical layer division are achieved, resulting in a geometry model containing several geometry objects.

[0064] In this embodiment, the specific method of step S4 includes the following steps:

[0065] (1) In hypermesh, the constructed geological geometry model is first divided into two-dimensional meshes. Select the surface plane and perform 2D automesh. Adjust the mesh size and mesh type as needed (taking a quadrilateral as an example, with a size of about 5km).

[0066] (2) Select the surface of the geometry model, enter the mesh size and mesh type, and click mesh;

[0067] (3) Adjust the number of nodes by left-clicking to add nodes and right-clicking to decrease nodes, thereby adjusting the grid. After dividing the grid, check it. If there are no problems, the two-dimensional grid division of the ground surface is complete, and proceed to the next step.

[0068] (4) Divide the 3D mesh according to the 2D mesh. The 3D mesh is divided into horizontal blocks and vertical layers. Select solid map in 3D, select the first block of the first vertical layer, and drag vertically according to the 2D mesh on this surface. Select the one volume function, and select the number of mesh layers to divide it into in the density option. Adjust the vertical mesh density to keep it as consistent as possible with the size of the 2D mesh. Then drag. Then perform the same operation on other blocks.

[0069] (5) Continue this process to complete the 3D mesh generation. During the 3D mesh generation process, because the mesh is generated in blocks and layers, each small mesh block is relatively independent, so any problems can be detected and adjusted promptly. After the 3D mesh generation is complete, the 2D mesh can be hidden or deleted.

[0070] This embodiment first divides and adjusts a two-dimensional grid based on the ground surface, and then drags and drops a three-dimensional grid according to the user's requirements for grid size and layer density to complete the three-dimensional grid division. This method has two advantages:

[0071] ① First, perform two-dimensional mesh generation, which is simpler, faster, and easier to adjust. This avoids the time-consuming and inconvenient process of directly generating a three-dimensional mesh, which may require re-generating the mesh if problems are found, thus wasting time.

[0072] ② Drag and drop the 3D mesh to complete the 3D mesh division. The advantage is that you can drag and drop layers vertically according to the blocks. The layer density of the mesh can be adjusted according to your needs. At the same time, since it is dragged layer by layer, you can find errors or problems in the division process in time and make layer corrections without having to re-divide the whole thing.

[0073] When organizing the 3D mesh model into blocks, the 3D mesh model is divided horizontally by the faults that cut laterally in the geometric model; and vertically by the lithospheric structure. The lithospheric structure varies in different regional environments. In this embodiment, the following division can be made: the first 20km is the upper crust, 20km~30km is the middle crust, 30km~63km is the lower crust, and 63km~100km is the upper mantle.

[0074] In this embodiment, the operation diagram of step S6 is as follows: Figure 5 As shown, the formed contact surface is as follows Figure 6 As shown. After creating all the faults, display them in the mesh model. Export the contact surfaces and the mesh model together to Abaqus, selecting the desired export folder.

[0075] Example 2:

[0076] This second embodiment is a further improvement upon the first embodiment. In this embodiment, a finite element model analysis method for geological bodies with three-dimensional fine-grained fault structures is provided, which includes the following steps:

[0077] A1. Import the geological finite element model .inp file from Hypermesh software into Abaqus software;

[0078] A2, such as Figure 8 , Figure 9 and Figure 10 As shown, in the Abaqus software, based on the actual research situation, different layers and blocks are assigned corresponding material medium and elastic parameters, and contact conditions and boundary conditions (such as velocity boundaries) are set. Alternatively, fault activity can be represented by setting contact relationships, or fault zone units on both sides of the fault plane can be assigned low-strength medium properties to represent the fault.

[0079] A3. Create time steps. Select either static or dynamic time steps as needed, customize the step size and name of each time step, set the output field according to requirements, submit the analysis and calculation, and obtain the generated odb file.

[0080] A4. Open the generated odb file, create a displaygroup, extract the fault plane, visualize the results on the fault plane, and show the dynamic changes of the physical field on the fault plane over time.

[0081] In this embodiment, as Figure 13 , Figure 14 , Figure 15 and Figure 16As shown, by creating a displaygroup, you can directly select the contact surfaces created in Hypermesh from the surface. You can choose to select one side or both sides to extract the fracture surface. You can directly visualize the results on the fracture surface, such as stress vectors and strain rates. At the same time, you can see the dynamic changes of strain rate or other physical fields on the fracture surface over time. You can make videos or dynamic graphs to display them, which will facilitate further in-depth analysis.

[0082] This embodiment allows direct extraction of model node data and output field data using a Python program, followed by further computational processing or mapping as needed. However, sometimes the model mesh is too large, making complete extraction difficult. Therefore, precise data extraction can be performed on a case-by-case basis. Select "file," then "Run script," and choose the desired script to run.

[0083] When accurately extracting data for subsequent processing, taking fault planes as an example, the required fault node data can be extracted from the CAE or INP file using a Python program and saved as a CSV file. Then, a Python program can be written to extract the required time steps and output field data from the ODB file based on the node CSV file. Finally, the program can be used for plotting.

[0084] In addition to the cross-sectional area, different display groups can be created according to your needs to extract the required time steps, node data and output field data, so that more complex calculations or plotting can be performed using Python programs.

[0085] At this point, the creation and post-processing analysis of the finite element model of the geological body with a detailed three-dimensional fault structure are complete.

[0086] In the specific implementation of this invention, taking the Longmenshan area as an example, all published literature and active fault mapping reports related to the faults were collected to find information on the surface geometric distribution and segmentation, motion properties, subsurface cutting depth, dip direction, and dip angle of the target fault. Simultaneously, based on surface fault traces, earthquake catalogs of different types, crustal velocity structure models, coseismic rupture, earthquake relocation and focal mechanism solutions, earthquake reflection profiles, and velocity structures, a detailed three-dimensional fault structure model was performed on SKUA-GOCAD, covering the spatial extension, dip direction, and dip angle variation characteristics of the fault. Finally, the detailed three-dimensional structure of the Longmenshan fault was obtained using this method. Figure 17 As shown.

[0087] In summary, this invention, by appropriately smoothing the three-dimensional fault model, not only preserves the three-dimensional geometric features of the fault but also facilitates subsequent mesh generation. Furthermore, even with relatively complex geometries, hexahedral elements are used during mesh generation, which promotes more accurate computational simulation, ultimately resulting in a finite element model of a geological body containing a fine three-dimensional fault structure.

Claims

1. A finite element modeling method for geological bodies with fine three-dimensional fault structures, characterized in that, Includes the following steps: Acquire geological and geophysical data, and construct three-dimensional fault models using geological three-dimensional modeling software; Import the original three-dimensional fault surface from the three-dimensional fault model into Rhino software, and smooth the original three-dimensional fault surface using Rhino software to obtain the smoothed three-dimensional fault surface. The smoothed fault 3D surface is imported into HyperMesh software for geometric partitioning of blocks, resulting in a geometric model containing several geometric bodies. In HyperMesh software, the geometric model is meshed to obtain a three-dimensional mesh model; In HyperMesh software, the 3D mesh model is divided into blocks based on the geometric model to obtain the divided 3D mesh model blocks. Enter the contact module, create contact surfaces, select the segmented 3D mesh model block, and select the mesh surface according to the location and depth of the fault plane to obtain the finite element model fault plane; for each finite element model fault, select two surfaces on both sides of the corresponding block to form contact surfaces; retain fault element bands of a set thickness on both sides of the fault interface. All three-dimensional mesh model blocks are combined into a geological finite element model, and the relative motion of the two disks of the fault layer in the geological finite element model is simulated by setting different contact relationships on the contact surfaces. The specific method for smoothing the original three-dimensional fault surface using Rhino software includes the following steps: The RhinoResurf function discretizes the original three-dimensional fault surface in the three-dimensional fault model into a point set, and then exports the point set for Z-axis inversion. Import the point set after Z-axis inversion into Rhino software for adjustment, and then form a surface based on the adjusted point set; the adjustment method includes deleting points on the cross-section whose protrusion value exceeds the threshold. A mesh is generated from the adjusted point set, the mesh border is extracted, and the face formed by the point set is trimmed with the mesh border using the Trim command to obtain the smoothed fault 3D surface. Specific methods for meshing geometric models in HyperMesh software include: In HyperMesh software, perform 2D mesh generation on the geometric model: Select the ground plane and perform 2dautomesh processing, where the mesh size and mesh type are adjusted as needed; For a geometric model that has been meshed in two dimensions, select the surface of the geometric module, input the mesh size and mesh type, click the generate button, and adjust the number of control nodes to optimize the quality of the two-dimensional mesh; To perform single-layer 3D mesh generation: Select the solid map function in 3D, select the first block of a certain vertical layer, and drag it vertically according to the 2D mesh on that surface. Adjust the vertical mesh density to be consistent with the 2D mesh size by dragging to obtain the corresponding 3D mesh for that layer. Repeat the single-layer 3D mesh generation operation until the corresponding 3D mesh for each layer is obtained, thus obtaining the 3D mesh model.

2. The finite element modeling method for geological bodies with three-dimensional fine structures containing faults according to claim 1, characterized in that, The specific methods for acquiring geological and geophysical data and constructing three-dimensional fault models using geological three-dimensional modeling software include the following steps: Collect geological, seismic, and geophysical data; Unify various types of detection data, observation results, and maps from different depths from the surface to underground into the same coordinate system and load them into 3D modeling software or platform; the 3D modeling software or platform is SKUA-GoCad, Petrel, or GOCAD; The various types of loaded data are interpreted to identify, extract, and establish the spatial location and attribute information of active faults; the spatial location of active faults includes point data, line data, and area data. The surface function is used to form surfaces from point data using the PointSet function; the Curve function is used to form surfaces from line data; and the resulting surfaces are used as preliminary cross-sectional surfaces. Based on the preliminary fault plane, and taking into account geological patterns and spatial topological relationships, the fault point data, line data, and surface data obtained from the interpretation are spatially interpolated using an interpolation algorithm to generate a continuous and overall smooth fault plane, thus obtaining a three-dimensional fault model.

3. The finite element modeling method for geological bodies with three-dimensional fine structures containing faults according to claim 2, characterized in that, Point data includes fault points measured in the field, fault points exposed by trench excavation, fault points exposed by boreholes, fault points interpreted by combined borehole profiles, and fault points interpreted by artificial seismic exploration profiles. Line data includes fault lines obtained from active fault mapping, fault lines interpreted from borehole joint profiles, fault lines interpreted from artificial seismic exploration profiles, fault lines interpreted from strong earthquake aftershock relocation results, fault surface traces obtained from InSAR inversion, and fault lines interpreted from non-seismic geophysical exploration. The surface data includes fault planes interpreted from 3D seismic exploration, fault nodal planes obtained from strong earthquake focal mechanism solutions, seismogenic fault planes obtained from coseismic deformation observation data inversion, and fault planes obtained from joint inversion using 3D geophysical methods.

4. The finite element modeling method for geological bodies with three-dimensional fine structures containing faults according to claim 1, characterized in that, The specific method for importing the smoothed 3D fault surface into HyperMesh software for geometric partitioning of blocks includes the following steps: In Rhino software, create a plane representing the study area and create different planes based on the fault incision depth. Export the created planes and the smoothed 3D fault surface to HyperMesh software. In HyperMesh, for a plane with Z=0, select the plane and drag it to form a solid; during the plane dragging process, select the depth required for the study to obtain the geological body to be cut; Perform geological body cutting to obtain a geometric model containing several geometric shapes: For geological bodies that need to be cut according to fault depth, select solid editing and select the trim with plane / surf function. Select the entity and face to be cut in sequence. Use the geometric plane with the maximum fault depth to cut the geological body at a depth. Then cut the geological body again through the fault plane at that depth to complete one complete cutting operation. Select the second deepest geometric plane of the last complete cutting operation and repeat the complete cutting operation process to obtain several geometric bodies. For geological bodies that do not require cutting based on fault depth, first cut the shallow geological body to obtain cutting traces on the shallow surface; then select the trim with lines function to use the cutting traces on the shallow surface to cut the remaining geological body. Alternatively, for geological bodies that do not require cutting based on fault depth, select the trim with plane / surf function to cut the geological body with planes of different depths to obtain several geometric shapes.

5. The finite element modeling method for geological bodies with three-dimensional fine structures containing faults according to claim 1, characterized in that, In HyperMesh software, when meshing a geometric model, select the number of mesh layers to be created using the density option in the one volume function.

6. The finite element modeling method for geological bodies with three-dimensional fine structures containing faults according to claim 1, characterized in that, In HyperMesh software, specific methods for dividing and organizing a 3D mesh model into blocks based on the geometric model include: The three-dimensional mesh model is divided horizontally by faults that cut laterally in the geometric model; and vertically by the lithosphere structure.

7. A finite element model analysis method for geological bodies with three-dimensional fine structures containing faults, based on the finite element modeling method for geological bodies with three-dimensional fine structures containing faults as described in any one of claims 1 to 6, characterized in that, Includes the following steps: Import the geological finite element model .inp file from Hypermesh software into Abaqus software; In the Abaqus software, based on the actual research situation, different layers and blocks are assigned corresponding material medium and elastic parameters, and contact conditions and boundary conditions are set. Create time steps, select static or dynamic time steps as needed, customize the step size and name of each time step, set the output field according to requirements, submit analysis and calculation, and obtain the generated odb file; Open the generated odb file, create a displaygroup, extract the fault plane, visualize the results on the fault plane, and show the dynamic changes of the physical field on the fault plane over time.

8. The finite element model analysis method for geological bodies with three-dimensional fine structures containing faults according to claim 7, characterized in that, When creating a displaygroup and extracting the fracture surface, select the contact surface created in the HyperMesh software from the surface function. The selection method can be either selecting one side or both sides. The results on the fracture surface include stress vector, strain rate, and displacement. The dynamic changes of the physical field on the fracture surface over time include the changes in strain rate, stress, and displacement field on the fracture surface over time.

Images