Modeling method and analysis method of geologic body finite element model containing fault three-dimensional fine structure

Abstract

The invention discloses a modeling method and an analysis method of a geologic body finite element model containing a fault three-dimensional fine structure, and relates to the field of earth science. According to the method, the three-dimensional fault model is properly smoothed, so that the three-dimensional geometrical characteristics of the fault are reserved, and subsequent grid division is facilitated. And meanwhile, during grid division, even if a relatively complex geometric form is adopted, hexahedron units are still adopted, so that more accurate computational simulation is facilitated, and finally, a geologic body finite element model containing a fault three-dimensional fine structure is obtained. According to the invention, flexible adjustment can be carried out according to different geological backgrounds and scientific problems, including but not limited to addition of corresponding material media, friction coefficients, physical equations, fault contact and other parameters, and geologic bodies and faults in a research area can be simulated more accurately and finely. Meanwhile, the geologic body finite element model mesh created by the method can reach million level, and more refined research, simulation and processing work are facilitated.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of earth science, in particular to a modeling method and an analysis method of a finite element model of a geological body containing a three-dimensional fine structure of a fault. BACKGROUND

[0002] The fault geometry in the three-dimensional finite element geological body model widely created in the current research is often relatively simple, and the grid division is relatively sparse, the number of grids is small, and the size of the grid unit is large, which greatly simplifies the real geological body condition. Thus, it is difficult to carry out further high-precision research based on the three-dimensional fine structure of the fault.

[0003] In the process of creating a three-dimensional finite element model, the complex boundary geometry division between blocks, the embedding of fault structure and the grid division are all key points and difficulties. Especially, the geometry division of blocks and the grid division based on the three-dimensional fine structure of the fault model are even more difficult. SUMMARY

[0004] In view of the above problems in the prior art, the modeling method and the analysis method of the finite element model of the geological body containing the three-dimensional fine structure of the fault provided by the present application solve the problems that the existing finite element model of the geological body is difficult to retain the three-dimensional fine structure of the fault, and it is difficult to carry out further high-precision stress and strain evolution analysis for the complex geometric structure fault.

[0005] In order to achieve the above-mentioned application purposes, the technical scheme adopted by the present application is as follows: A modeling method of a finite element model of a geological body containing a three-dimensional fine structure of a fault is provided, which comprises the following steps: Obtaining geological and geophysical data, and constructing a three-dimensional fault model through a three-dimensional geological modeling software; Importing the original three-dimensional surface of the fault in the three-dimensional fault model into Rhino software, and performing smoothing processing on the original three-dimensional surface of the fault through the Rhino software to obtain a smoothed three-dimensional surface of the fault; Importing the smoothed three-dimensional surface of the fault into HyperMesh software to perform geometry division of blocks, and obtaining a geometry model containing a plurality of geometric bodies; Performing grid division on the geometry model in the HyperMesh software to obtain a three-dimensional grid model; Performing block arrangement on the three-dimensional grid model according to the geometry model in the HyperMesh software, and obtaining a block of the three-dimensional grid model arranged and divided; Enter the contact module, create the contact surface pair, select the three-dimensional grid model block, select the grid surface according to the position and depth of the fault surface, and obtain the finite element model fault surface; wherein each finite element model fault selects two surfaces on both sides of the corresponding block to form a contact pair; the fault interface on both sides retains a fault element zone with a set thickness; All three-dimensional grid model blocks are combined into a geological body finite element model, and different contact relationships are set for the contact pairs to simulate the relative movement of the two plates of the fault in the geological body finite element model.

[0006] A geological body finite element model analysis method based on a three-dimensional fine structure containing a fault, comprising the following steps: Import the geological body finite element model.inp file from the hypermesh software into the abaqus software; In the abaqus software, according to the actual research situation, different layers and blocks are given corresponding material medium and elastic parameters, and the contact condition and boundary condition are set; Create a time step, select a static or dynamic time step as needed, customize the step size and name of each time step, set the output field according to the requirement, submit the analysis calculation, and obtain the generated odb file; Open the generated odb file, create a display group, extract the fault surface, visualize the results on the fault surface, and show the dynamic changes of the physical field on the fault surface over time.

[0007] The beneficial effects of the present application are: The present application can be adjusted flexibly according to different geological backgrounds and scientific problems, including but not limited to adding corresponding material medium, friction coefficient, physical equation and fault contact parameters, and can more accurately simulate the geological body and three-dimensional fault of the research area. At the same time, the grid of the geological body finite element model created by the present application can reach millions, which is beneficial to more accurate research, simulation and processing work.

[0008] The present application can be adjusted flexibly according to different geological backgrounds and scientific problems, including but not limited to adding corresponding material medium, friction coefficient, physical equation and fault contact parameters, and can more accurately simulate the geological body and three-dimensional fault of the research area. At the same time, the grid of the geological body finite element model created by the present application can reach millions, which is beneficial to more accurate research, simulation and processing work. BRIEF DESCRIPTION OF DRAWINGS

[0009] Figure 1 It is a flowchart of the present three-dimensional fine structure containing a fault geological body finite element modeling method; Figure 2 Operation diagram for data interpretation and information extraction in Example 1; Figure 3 Operation diagram for exporting the constructed three-dimensional fault model into a DXF format data file in Example 1; Figure 4 Operation diagram for selecting entity editing in Example 1; Figure 5 Operation diagram for entering the contact module in step S6 in Example 1; Figure 6 Operation diagram for creating contact surface pairs in step S6 in Example 1; Figure 7 Diagram of the contact surface pairs formed in step S6 in Example 1; Figure 8 Operation diagram for assigning corresponding material media and elastic parameters to different layers and blocks in Example 2; Figure 9 Operation diagram for setting contact conditions in Example 2; Figure 10 Operation diagram for setting boundary conditions in Example 2; Figure 11 Operation diagram for creating time steps in Example 2; Figure 12 Operation diagram for setting output fields in Example 2; Figure 13 Diagram of the strain rate results of the simulation calculation in Example 2; Figure 14 Diagram of the display group setting in Example 2; Figure 15 Diagram of the simulation calculation results on the fault surface in Example 2; Figure 16 Vector diagram of the simulation calculation results on the fault surface in Example 2; Figure 17 Diagram of the three-dimensional fine structure model of the Longmenshan fault in Example 1. DETAILED DESCRIPTION

[0010] The specific embodiments of the present application are described below to facilitate the understanding of the present application for those skilled in the art, but it should be clear that the present application is not limited to the scope of the specific embodiments, and for those skilled in the art, any changes that are obvious within the spirit and scope of the present application as defined and determined by the appended claims are within the scope of the present application.

[0011] Example 1: As Figure 1 shown, the finite element modeling method of the geological body containing a three-dimensional fine structure of a fault includes the following steps: S1, obtaining geological and geophysical data, and constructing a three-dimensional fault model through a three-dimensional geological modeling software; S2, importing the original fault three-dimensional surface in the three-dimensional fault model into Rhino software, and performing smoothing processing on the original fault three-dimensional surface through the Rhino software to obtain a smoothed fault three-dimensional surface; S3, importing the smoothed fault three-dimensional surface into HyperMesh software to perform geometric division of blocks, and obtaining a geometric model containing a plurality of geometric bodies; S4, performing mesh division on the geometric model in the HyperMesh software to obtain a three-dimensional mesh model; S5, performing block arrangement on the three-dimensional mesh model according to the geometric model in the HyperMesh software, and obtaining a three-dimensional mesh model block arranged and divided; S6, entering the contact module, creating a contact surface pair, selecting the three-dimensional mesh model block arranged and divided, selecting a mesh surface according to the position and depth of the fault surface, and obtaining a finite element model fault surface; wherein each finite element model fault selects two surfaces on both sides of the corresponding block to form a contact surface pair; a fault element zone with a set thickness is reserved on both sides of the fault boundary surface; S7, assembling all the three-dimensional mesh model blocks into a geological body finite element model, and simulating the relative movement of the two plates of the fault in the geological body finite element model by setting different contact relationships for the contact surface pairs.

[0012] In this embodiment, the acquisition of geological and geophysical data needs to collect geological and seismic data and geophysical data. Data integration and loading are also needed, that is, various types of detection data, observation results and maps at different depths from the ground to the underground are unified in the same coordinate system and loaded into a three-dimensional modeling software or platform, such as SKUA-GoCad, Petrel, GOCAD, etc. The following takes SKUA-GoCad as an example: Data interpretation and information extraction: interpreting the loaded various types of data, identifying and extracting the spatial position (point, line, surface form) and attribute information of the active fault. For example Figure 2As shown, select import, according to the data format such as point for csv format, then can select column-based file for import, line and surface are generally shp file format, select ArcViewShp for file import. Then select surface function, according to the imported point (small earthquake data) to select PointSet to generate surface from point, Curve to generate surface from multiple fault profile lines, and the fault surface is preliminarily built.

[0013] The point elements mainly include field survey measured fault points, trench excavation revealed fault points, borehole revealed fault points, or interpreted fault points by borehole joint profile, and artificial earthquake exploration profile interpreted upper fault points; the line elements mainly include active fault mapping obtained fault lines, borehole joint profile interpreted fault lines, artificial earthquake exploration profile interpreted fault lines, based on strong earthquake aftershock relocation results interpreted fault lines, InSAR inversion obtained fault surface traces, and non-seismic geophysical exploration interpreted fault lines; and the surface elements mainly include three-dimensional seismic exploration interpreted fault surfaces, strong earthquake source mechanism solution obtained fault planes, coseismic deformation observation data inversion obtained seismogenic fault surfaces, and three-dimensional geophysical method joint inversion obtained fault surfaces.

[0014] Based on the above established preliminary fault surface, based on geological rules and spatial topological relations, appropriate interpolation algorithm is selected to perform spatial interpolation on the interpreted fault point, line and surface data, to generate continuous and overall smooth fault surface, and to construct three-dimensional fault model. As shown in the following figure, Figure 3 As shown in the following figure, the finally constructed three-dimensional fault model is exported as a DXF format data file.

[0015] In this embodiment, the specific method of step S2 includes the following steps: (1) Import the original fault three-dimensional surface in the three-dimensional fault model into Rhino, use RhinoResurf to disperse the surface into point set, and export the point set for Z axis inversion, because in general case, the fault surface is positive with Z axis downward, while in Rhino, it is positive with Z axis upward; (2) Import the inverted point set into Rhino, import the adjusted fault points into Rhino, adjust the fault points, generate surface according to the point cloud, and adjust the control points and order to make the generated surface smoother; (3) Generate mesh from point, extract mesh frame, cut the surface by frame through Trim instruction, and then generate smooth fault three-dimensional surface.

[0016] In this embodiment, the specific method of step S3 includes the following steps: (1) Create the plane of the study area in Rhino, and create different planes according to the depth of the fault, export the plane and the fault plane in.igs format from Rhino, and import them into HyperMesh; (2) In HyperMesh, select the plane dragging into a body at Z=0, and select the depth according to the study; (3) For the geological body which needs to be cut according to the depth of the fault, select solid edit and trim with plane / surf, select the entity to be cut and the plane in turn, cut the geological body according to the maximum depth of the plane, and cut the geological body again according to the fault plane of the depth, complete a complete cutting operation; select the next deep geometric plane of the last complete cutting operation, repeat the complete cutting operation process, and obtain several geometric bodies; (4) If the cutting according to the depth of the fault is not needed, the shallow layer (such as 20 km) can be cut first, and then trim with lines is selected to cut the remaining entity with the cutting trace on the 20 km plane cut by the last step; (5) Trim with plane / surf is selected to cut with planes of different depths; (6) At this point, the construction of the geometric body is completed, and the horizontal block and vertical layer are divided, and the geometric body model containing several geometric bodies is obtained.

[0017] In the embodiment, the specific method of step S4 includes the following steps: (1) In hypermesh, the geological body geometric model constructed is first subjected to two-dimensional meshing, the ground plane is selected, and 2d automesh is performed, wherein the mesh size and the mesh type (herein, the quadrilateral is taken as an example, and the size is about 5 km) are adjusted according to the requirement; (2) The surface of the geometric body model is selected, the mesh size (element size) and the mesh type (meshtype) are input, and mesh is clicked to complete the operation; (3) The number of nodes is adjusted by left-clicking the node number addition and right-clicking the reduction, and then the mesh is adjusted. After the meshing is completed, check is performed, and the operation is completed if no problem is found. Then, the next step is performed. (4) According to the two-dimensional grid to divide the three-dimensional grid, the three-dimensional grid is divided by transverse block and longitudinal layer. Select solid map in 3D, select the first block in the first layer, drag according to the two-dimensional grid on the surface, select one volume function, select the number of layers to be divided into at the density option, adjust the grid density in the longitudinal direction, and try to keep consistent with the size of the two-dimensional grid. Then the same operation is performed on other blocks; (5) Similarly, complete the three-dimensional grid division. During the three-dimensional grid division, since the division is performed by block and layer, the grid of each small block is relatively independent, so the problem can be found and adjusted. After the three-dimensional grid division is completed, the two-dimensional grid can be hidden or deleted.

[0018] The embodiment first divides the two-dimensional grid according to the ground surface and adjusts it, and then performs three-dimensional grid dragging according to the requirements of grid size and layer density to complete the three-dimensional grid division. This method has two advantages: ①Divide the two-dimensional grid first, which is simpler, faster and easier to adjust, so as to prevent the direct division of the three-dimensional grid from being time-consuming and inconvenient to adjust, and the problem needs to be re-divided, which is time-consuming.

[0019] ②The three-dimensional grid is divided by dragging, which has the advantage that the layer can be dragged layer by layer according to the block, the layer density of the grid can be adjusted according to the needs, and since it is layer-by-layer dragging, errors or problems in the division process can be found in time and corrected layer by layer without the need to re-divide.

[0020] When the three-dimensional grid model is divided into blocks, the three-dimensional grid model is divided horizontally according to the horizontal cutting of the geometric model; the three-dimensional grid model is divided vertically according to the structure of the lithosphere. The lithospheric structure of different regional environments will be different, and in this embodiment, the following division can be made: the first 20 km is the upper crust, 20 km~30 km is the middle crust, 30 km~63 km is the lower crust, and 63 km~100 km is the upper mantle.

[0021] In this embodiment, the operation schematic diagram of step S6 is as shown in Figure 5 The contact surface formed is as shown in Figure 6 After creating all the faults, display them in the grid model, and export the contact surface and the grid model to abaqus together. Select the folder to be exported to the folder.

[0022] Example 2: This embodiment is a further improvement based on embodiment one. In this embodiment, a fault-containing three-dimensional fine structure geological body finite element model analysis method based on a fault-containing three-dimensional fine structure geological body finite element model modeling method is provided, which comprises the following steps: A1, import the geological body finite element model.inp file from the hypermesh software into the abaqus software; A2, as shown in Figure 8 , Figure 9 and Figure 10 , in the abaqus software, according to the actual research situation, different layers and different blocks are given corresponding material medium and elastic parameters, and the contact condition and boundary condition (such as velocity boundary) are set; it can also be selected to present fault activity by setting contact relationship, or to give low strength medium properties to the fault zone elements on both sides of the fault to represent the fault; A3, create time step, select static or dynamic time step according to need, customize the step length and name of each time step, set output field according to demand, submit analysis calculation, and obtain generated odb file; A4, open the generated odb file, create displaygroup, extract the fault surface, visualize the results on the fault surface and show the dynamic changes of the physical field on the fault surface over time.

[0023] In this embodiment, as shown in Figure 13 , Figure 14 , Figure 15 and Figure 16 , create displaygroup, you can directly select the contact surface created in hypermesh from surface, you can choose single side or double side, so as to extract the fault surface, you can directly visualize the results on the fault surface, such as stress vector, strain rate, etc., at the same time, you can see the dynamic changes of strain rate or other physical field on the fault surface over time, you can make video or dynamic graph to show, which is convenient for subsequent in-depth analysis.

[0024] This embodiment can directly use python program to extract model node data and output field data, and perform more subsequent calculation and processing or mapping processing according to need. But sometimes the model grid is too much, and it is difficult to extract all, so the accurate data extraction can be carried out according to need. Select file, Run script, select the code script you need to run.

[0025] When the data is extracted accurately for subsequent processing, taking the fault surface as an example, the required fault node data can be extracted from the surface in the cae file or inp file and saved as a csv file through a python program, a python program is written to extract the required time step and output field data from the odb file according to the node csv file, and finally the data is processed into a graph through the program.

[0026] In addition to the fault surface, different display groups can also be created according to the needs, and the required time step, node data and output field data are extracted, so that more complex calculation and processing or mapping processing can be carried out through the python program.

[0027] At this point, the creation and post-processing analysis process of the entire three-dimensional fine structure of the geological body finite element model containing the fault are completed.

[0028] In the specific implementation process of the present application, the Longmenshan area is taken as an example, all published literatures and active fault mapping reports related to the fault are collected, and information on the target fault surface geometric distribution and segmentation, motion properties, underground cutting depth, tendency and inclination is found. At the same time, according to the surface fault trace, different types of earthquake catalog, crustal velocity structure model, coseismic rupture, earthquake relocation and focal mechanism solution, seismic reflection profile and velocity structure data, the three-dimensional fine structure modeling of the fault is carried out on SKUA-GOCAD, which covers the characteristics of the spatial extension, tendency and inclination of the fault, and finally the three-dimensional fine structure of the Longmenshan fault is obtained by the method as shown in Figure 17 .

[0029] As described above, the three-dimensional fault model is appropriately smoothed in the present application, not only the three-dimensional geometric characteristics of the fault are retained, but also the subsequent mesh division is facilitated. At the same time, even if the geometric shape is relatively complex, hexahedral elements are still used in the mesh division, which is conducive to more accurate calculation simulation, so that the three-dimensional fine structure of the geological body finite element model containing the fault is finally obtained.

Claims

1. A method for modeling a finite element model of a geological body containing a three-dimensional fine structure of a fault, characterized by, The method comprises the following steps: Obtaining geological and geophysical data, and constructing a three-dimensional fault model through a three-dimensional geological modeling software; Importing the original three-dimensional fault surface in the three-dimensional fault model into Rhino software, and performing smoothing processing on the original three-dimensional fault surface through the Rhino software to obtain a smoothed three-dimensional fault surface; Importing the smoothed three-dimensional fault surface into HyperMesh software to perform geometric division of blocks, and obtaining a geometric model comprising a plurality of geometric bodies; Performing mesh division on the geometric model in the HyperMesh software to obtain a three-dimensional mesh model; Performing block arrangement on the three-dimensional mesh model according to the geometric model in the HyperMesh software, and obtaining a three-dimensional mesh model block arranged and divided; Entering a contact module, creating a contact surface pair, selecting the three-dimensional mesh model block arranged and divided, selecting mesh surfaces according to the positions and depths of the fault surfaces, and obtaining a finite element model fault surface; wherein each finite element model fault selects two surfaces on the two sides of the corresponding block to form a contact surface pair; a fault element zone with a set thickness is reserved on the two sides of the fault boundary surface; Assembling all the three-dimensional mesh model blocks into a geological body finite element model, and simulating the relative movement of the two plates of the fault in the geological body finite element model by setting different contact relationships for the contact surface pairs.

2. The faulted three-dimensional fine structure containing geological body finite element model modeling method according to claim 1, characterized in that, The specific method for obtaining geological and geophysical data and constructing a three-dimensional fault model through a three-dimensional geological modeling software comprises the following steps: Collecting geological, seismic and geophysical data; Unifying various types of exploration data, observation results and drawings at different depths from the ground to the underground into the same coordinate system, and loading them into a three-dimensional modeling software or platform; wherein the three-dimensional modeling software or platform is SKUA-GoCad, Petrel or GOCAD; Interpreting the loaded various types of data, identifying, extracting and establishing the spatial position and attribute information of active faults; wherein the spatial position of the active faults comprises point data, line data and surface data; Forming a surface by PointSet function through point data through a surface function; forming a surface through Curve function through line data; and obtaining each surface as a preliminary fault surface; Based on the preliminary fault surface, generating a continuous and overall smooth fault surface by spatial interpolation of the interpreted fault point data, line data and surface data through an interpolation algorithm based on geological rules and spatial topological relations, thereby obtaining a three-dimensional fault model.

3. The faulted three-dimensional fine structure containing geological body finite element modeling method according to claim 2, characterized in that, The point data comprises field investigation and measurement fault points, trench excavation revealed fault points, drilling revealed fault points, joint profile interpreted fault points through drilling, and artificial earthquake exploration profile interpreted fault points; The line data comprises active fault mapping obtained fault lines, joint profile interpreted fault lines through drilling, artificial earthquake exploration profile interpreted fault lines, fault lines interpreted based on strong earthquake aftershock relocation results, InSAR inversion obtained fault surface traces, and non-seismic type geophysical exploration interpreted fault lines; The fault surface data include a fault surface interpreted from three-dimensional seismic exploration, a fault surface obtained from a strong earthquake focal mechanism solution, a fault surface obtained from coseismic deformation observation data inversion, and a fault surface obtained from three-dimensional geophysical method joint inversion.

4. The faulted three-dimensional fine structure containing geological body finite element modeling method according to claim 1, characterized in that, The specific method for smoothing the original three-dimensional fault surface by using Rhino software includes the following steps: The original three-dimensional fault surface in the three-dimensional fault model is discretized into a point set by using RhinoResurf function, and the point set is exported for Z-axis inversion; The point set after Z-axis inversion is imported into Rhino software for adjustment, and a surface is formed according to the adjusted point set; wherein the adjustment method includes deleting points with a value exceeding a threshold on the fault surface; A grid is generated from the adjusted point set, the grid frame is extracted, the surface formed by the point set is trimmed by the grid frame through the Trim instruction, and the smoothed three-dimensional fault surface is obtained.

5. The faulted three-dimensional fine structure containing geological body finite element modeling method according to claim 1, characterized in that, The specific method for importing the smoothed three-dimensional fault surface into HyperMesh software for geometric division of blocks includes the following steps: A plane in the research area range is created in Rhino software, different planes are created according to the fault under-cutting depth, and the created plane and the smoothed three-dimensional fault surface are exported to HyperMesh software; In HyperMesh, the plane Z=0 is selected to be dragged into a body; wherein the depth is selected as the required depth during the plane dragging process, and the geological body to be cut is obtained; Geometric body model containing several geometric bodies is obtained by cutting the geological body: For the geological body that needs to be cut according to the fault depth, select entity editing and select the trim with plane / surf function, and then select the entity to be cut and the surface in sequence. The geometric plane with the maximum fault depth is used to cut the geological body in depth, and the geological body is re-cut by the fault surface with the depth, and a complete cutting operation is completed. The next deep geometric plane of the last complete cutting operation is selected, and the complete cutting operation process is repeated to obtain several geometric bodies. For the geological body that does not need to be cut according to the fault depth, the shallow geological body is cut first to obtain the cutting trace on the shallow surface; then the trim with lines function is selected, and the cutting trace on the shallow surface is used to cut the remaining geological body; Or, for the geological body that does not need to be cut according to the fault depth, the trim with plane / surf function is selected, and the geological body is cut by planes with different depths to obtain several geometric bodies.

6. The faulted three-dimensional fine structure containing geological body finite element modeling method according to claim 1, characterized in that, The specific method for meshing the geometric body model in HyperMesh software includes: Two-dimensional meshing of the geometric body model in HyperMesh software: select the surface of the geometric body module, input the mesh size and mesh type, and click the generate button to adjust the number of control nodes for quality optimization of two-dimensional meshing; For the geometric body model that has been two-dimensionally meshed, select the surface of the geometric body module, input the mesh size and mesh type, and click the generate button to adjust the number of control nodes for quality optimization of two-dimensional meshing; Carrying out single-layer three-dimensional grid division: selecting the solid map function in 3D, selecting the first block in a certain layer in the longitudinal direction, and carrying out longitudinal dragging according to the two-dimensional grid on the surface, so as to adjust the longitudinal grid density, and the two-dimensional grid size is consistent, thereby obtaining the corresponding three-dimensional grid of the layer; Repeating the single-layer three-dimensional grid division operation until the corresponding three-dimensional grid of each layer is obtained, thereby obtaining the three-dimensional grid model.

7. The faulted three-dimensional fine structure containing geological body finite element modeling method according to claim 6, characterized in that, In the process of grid division of the geometric model in the HyperMesh software, the density option in the one volume function is selected to select the grid layer to be divided.

8. The faulted three-dimensional fine structure containing geological body finite element modeling method according to claim 1, characterized in that, The specific method for block arrangement of the three-dimensional grid model according to the geometric model in the HyperMesh software includes: The three-dimensional grid model is divided horizontally with the fault cut horizontally in the geometric model as the boundary; and the three-dimensional grid model is divided longitudinally according to the lithospheric structure.

9. A method of analyzing a finite element model of a geological body having a three-dimensional fine structure including a fault, based on the method of modeling a finite element model of a geological body having a three-dimensional fine structure including a fault according to any one of claims 1 to 8, characterized by, The method includes the following steps: The finite element model.inp file of the geological body is imported from the HyperMesh software into the Abaqus software; In the Abaqus software, different layers and blocks are given corresponding material media and elastic parameters according to the actual research situation, and the contact conditions and boundary conditions are set; The time step is created, the static or dynamic time step is selected according to the needs, the step length and name of each time step are customized, the output field is set according to the needs, the analysis calculation is submitted, and the generated odb file is obtained; The generated odb file is opened, the display group is created, the fault surface is extracted, the results on the fault surface are visualized, and the dynamic changes of the physical field on the fault surface over time are shown.

10. The faulted three-dimensional fine structure containing geological body finite element model analysis method according to claim 9, characterized by, When creating the display group and extracting the fault surface, the contact surface pair created in the HyperMesh software is selected from the surface function, and the selection mode is to select one side or both sides; the results on the fault surface include stress vector, strain rate and displacement; the dynamic changes of the physical field on the fault surface over time include the changes of strain rate, stress and displacement field on the fault surface over time.

Images