Simulation method and device for strike-slip fracture layering deformation process based on three-dimensional discrete element

Abstract

The invention provides a strike-slip fracture layering deformation process simulation method and device based on a three-dimensional discrete element, and relates to the technical field of geological structure deformation discrete element numerical simulation. The method comprises the steps that plane distribution characteristics of different strike-slip fractures in a research area are analyzed, drilling data and logging data are combined, and a stratigraphic mechanics layering pattern is obtained through division; according to the stratigraphic mechanics layering pattern, different mesoscopic parameters are given to a discrete element model, and a three-dimensional discrete element model with layering characteristics is constructed; according to the three-dimensional discrete element model, different extrusion conditions are set for iterative calculation, and a simulation result to be analyzed is determined; and comparing the simulation result with an analysis result of three-dimensional seismic interpretation, and determining the influence of different mechanical stratification on strike-slip fracture deformation characteristics.

Description

Technical Field

[0001] The present invention relates to the technical field of discrete element numerical simulation of geological structure deformation, and in particular to a method and device for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements. Background Art

[0002] This section is intended to provide a background or context for embodiments of the present invention. No description herein is admitted to be prior art by virtue of its inclusion in this section.

[0003] Strike-slip faults are characterized by deep burial depth, weak fault activity, and large differences in lithologic properties between the strata they penetrate, resulting in strong stratified deformation differences. The so-called stratified activity of strike-slip faults refers to the formation of multiple sets of strata with different capabilities within a certain area under the influence of multi-stage tectonic-sedimentary evolution. Influenced by differences in stress backgrounds at different times, strike-slip faults produce vertical stratified structural deformation. Large strike-slip faults are typically important areas for oil and gas enrichment, and the study of strike-slip faults plays an important role in controlling oil and gas migration, accumulation, and differential enrichment. Existing technologies do not offer a good solution for accurately analyzing the stratified activity and deformation process of strike-slip faults.

[0004] For example, existing discrete element numerical simulation methods mostly focus on simulating the evolution process of strike-slip fault structures, usually simplifying the strata into a homogeneous model and ignoring the impact of differences in mechanical properties between different strata on fault deformation.

[0005] In summary, existing technologies have significant deficiencies in analyzing the factors controlling stratified deformation of strike-slip faults. Existing methods limit our understanding of the stratified deformation characteristics of strike-slip faults and their tectonic deformation mechanisms.

[0006] Therefore, there is an urgent need for a technical solution that can overcome the above-mentioned defects and realize the numerical simulation of strike-slip faults based on mechanical stratification characteristics to enhance the understanding of the dynamic mechanism of longitudinal stratification deformation. Summary of the Invention

[0007] To address the challenges of existing technologies, this paper proposes a method and device for simulating the layered deformation process of strike-slip faults using three-dimensional discrete elements. By designing and analyzing three-dimensional discrete element numerical simulation experiments of strike-slip faults and comparing the simulation results with seismic profiles of strike-slip faults, this paper reveals the influence of factors controlling layered deformation on the deformation characteristics of strike-slip faults.

[0008] The overall solution solves the problem that traditional homogeneous models cannot reveal the layered deformation mechanism through the technical path of layered modeling, numerical simulation, and data verification. It provides an analysis tool that is closer to actual geological conditions for strike-slip fault research, and has both scientific value and engineering application potential.

[0009] In a first aspect of an embodiment of the present invention, a method for simulating the layered deformation process of a strike-slip fault based on a three-dimensional discrete element is proposed. The method includes:

[0010] The planar distribution characteristics of different strike-slip faults in the study area were analyzed, and the formation mechanical stratification pattern was obtained by combining drilling data and logging data.

[0011] According to the stratum mechanics stratification pattern, different microscopic parameters are assigned to the discrete element model to construct a three-dimensional discrete element model with stratification characteristics;

[0012] According to the three-dimensional discrete element model, different extrusion conditions are set to perform iterative calculations to determine the simulation results to be analyzed;

[0013] The simulation results are compared with the analysis results of 3D seismic interpretation to determine the influence of different mechanical layers on the deformation characteristics of strike-slip faults.

[0014] In a second aspect of an embodiment of the present invention, a device for simulating the stratified deformation process of a strike-slip fault based on a three-dimensional discrete element is proposed, the device comprising:

[0015] The mechanical layer analysis module is used to analyze the planar distribution characteristics of different strike-slip faults in the study area and, combined with drilling data and logging data, to obtain the formation mechanical layering pattern;

[0016] A three-dimensional discrete element model construction module is used to assign different microscopic parameters to the discrete element model according to the formation mechanics layering pattern, so as to construct a three-dimensional discrete element model with layered characteristics;

[0017] The discrete element numerical simulation module is used to set different extrusion conditions according to the three-dimensional discrete element model for iterative calculation to determine the simulation results to be analyzed;

[0018] The impact analysis module is used to compare the simulation results with the analysis results of the three-dimensional seismic interpretation to determine the impact of different mechanical layers on the deformation characteristics of the strike-slip fault.

[0019] In a third aspect of an embodiment of the present invention, a computer device is proposed, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, a method for simulating the stratified deformation process of a strike-slip fault based on three-dimensional discrete elements is implemented.

[0020] In a fourth aspect of an embodiment of the present invention, a computer-readable storage medium is proposed, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, a method for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements is implemented.

[0021] In a fifth aspect of an embodiment of the present invention, a computer program product is proposed, which includes a computer program. When the computer program is executed by a processor, a method for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements is implemented.

[0022] The method and device for simulating the layered deformation process of strike-slip faults based on three-dimensional discrete elements proposed in the present invention combine seismic interpretation profiles and drilling and logging data to accurately identify and describe the layered characteristics of different strata, thereby constructing a three-dimensional discrete element model with clear layered characteristics. By applying different extrusion speeds and extrusion distances to the discrete element model, the strike-slip shear process under different mechanical layering is simulated. This method can not only analyze the influence of each mechanical layering on the strike-slip fault deformation process, but also reveal the dominant role of layering characteristics in fault behavior, providing a new method for understanding the influence of different mechanical differences on strike-slip fault behavior. Comparing the simulation result profile with geophysical data further verified the key role of layering characteristics in the strike-slip fault process, and provided a strong experimental basis and theoretical reference for oil and gas exploration controlled by deep strike-slip faults. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following is a brief introduction to the drawings required for use in the description of the embodiments. Obviously, the drawings described below are some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0024] Figure 1 It is a flow chart of a method for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements according to an embodiment of the present invention.

[0025] Figure 2 It is a schematic diagram of a three-dimensional discrete element initial model of different mechanical layering characteristics according to a specific embodiment of the present invention.

[0026] Figure 3 It is a schematic diagram of the strike-slip fracture shear process of a specific embodiment of the present invention.

[0027] Figure 4 It is a schematic diagram of three-dimensional strike-slip fracture shear strain in a specific embodiment of the present invention.

[0028] Figure 5 It is a schematic cross-sectional view of a strike-slip fracture shear strain according to a specific embodiment of the present invention.

[0029] Figure 6 It is a schematic diagram of the architecture of a device for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements according to an embodiment of the present invention.

[0030] Figure 7It is a schematic diagram of the structure of a computer device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0031] The principles and spirit of the present invention will be described below with reference to several exemplary embodiments. It should be understood that these embodiments are provided solely to enable those skilled in the art to better understand and implement the present invention, and are not intended to limit the scope of the present invention in any way. Rather, these embodiments are provided to make this disclosure more thorough and complete, and to fully convey the scope of the present disclosure to those skilled in the art.

[0032] Those skilled in the art will appreciate that the embodiments of the present invention may be implemented as a system, apparatus, device, method, or computer program product. Therefore, the present disclosure may be implemented in the following forms: entirely in hardware, entirely in software (including firmware, resident software, microcode, etc.), or in a combination of hardware and software.

[0033] According to an embodiment of the present invention, a method and device for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements are proposed, which relates to the technical field of discrete element numerical simulation of geological structural deformation. The overall scheme of the present invention can break through the limitations of homogeneous models. Existing methods often ignore the mechanical differences of the strata. The present invention, through layered modeling, incorporates the mechanical properties of different strata (such as strength and friction coefficient) into discrete element simulation for the first time, revealing the main control mechanism of layered deformation. The present invention proposes multi-data fusion simulation, combining seismic data with drilling and logging data, to achieve accurate docking of geological characteristics with numerical simulation, and improve the authenticity and reliability of the model. By comparing the simulation results under different layers, the overall scheme clarifies the influence of mechanical stratification on the morphology and strain distribution of strike-slip faults, providing a new perspective for analyzing the dynamics of fault evolution. The simulation results can be combined with geophysical data to assist in identifying oil and gas enrichment areas controlled by deep strike-slip faults, providing experimental basis and theoretical reference for exploration deployment, and are suitable for the study of strike-slip faults under different geological conditions, expanding new scenarios for the application of discrete element methods in the field of structural deformation.

[0034] The principles and spirit of the present invention are explained in detail below with reference to several representative embodiments of the present invention.

[0035] Figure 1 FIG1 is a flow chart of a method for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements according to an embodiment of the present invention. Figure 1 As shown, the method includes:

[0036] S101: Analyze the planar distribution characteristics of different strike-slip faults in the study area and, combined with drilling and logging data, obtain the stratigraphic stratification pattern.

[0037] S102, assigning different mesoscopic parameters to the discrete element model according to the formation mechanics layering pattern, and constructing a three-dimensional discrete element model with layered characteristics;

[0038] S103, setting different extrusion conditions according to the three-dimensional discrete element model to perform iterative calculations to determine the simulation results to be analyzed;

[0039] S104, comparing the simulation results with the analysis results of the 3D seismic interpretation to determine the influence of different mechanical layers on the deformation characteristics of the strike-slip fault.

[0040] The method and device for simulating the layered deformation process of strike-slip faults based on three-dimensional discrete elements proposed in the present invention combine seismic interpretation profiles and drilling and logging data to accurately identify and describe the layered characteristics of different strata, thereby constructing a three-dimensional discrete element model with clear layered characteristics. By applying different extrusion speeds and extrusion distances to the discrete element model, the strike-slip shear process under different mechanical layering is simulated. This method can not only analyze the influence of each mechanical layering on the strike-slip fault deformation process, but also reveal the dominant role of layering characteristics in fault behavior, providing a new method for understanding the influence of different mechanical differences on strike-slip fault behavior. Comparing the simulation result profile with geophysical data further verified the key role of layering characteristics in the strike-slip fault process, and provided a strong experimental basis and theoretical reference for oil and gas exploration controlled by deep strike-slip faults.

[0041] In order to explain more clearly the above-mentioned method for simulating the layered deformation process of strike-slip faults based on three-dimensional discrete elements, each step is explained in detail below.

[0042] In one embodiment, S101 analyzes the planar distribution characteristics of different strike-slip faults in the study area, combines drilling and logging data, and divides the formation mechanical stratification pattern. The specific process includes:

[0043] The planar distribution characteristics of different strike-slip faults in the study area were analyzed using seismic interpretation profiles, and the planar and cross-sectional deformation characteristics of deep faults were obtained through analysis.

[0044] Acquire drilling data and logging data, and analyze vertical changes, lithology, thickness, and physical and mechanical properties of the formation based on the planar and cross-sectional deformation characteristics and in combination with the drilling data and logging data to determine mechanical differences in the formation;

[0045] According to the mechanical differences of the strata, the stratigraphic mechanical stratification pattern in the study area was divided.

[0046] Specifically, the team used seismic interpretation profiles to analyze the planar distribution of different strike-slip faults within the study area, and further analyzed the planar and cross-sectional deformation characteristics of deep faults. In combination with drilling and logging data, the team divided the stratigraphic stratification pattern and analyzed the evolution of structural deformation and its vertical development pattern.

[0047] In one embodiment, S102 , different microscopic parameters are assigned to a discrete element model according to the formation mechanics layering pattern, so as to construct a three-dimensional discrete element model with layered characteristics.

[0048] Specifically, a three-dimensional discrete element model is constructed according to the size, base and lateral boundaries of a preset three-dimensional simulation box; wherein the rectangular area enclosed by the base and the lateral boundaries corresponding to the multiple lateral baffles is used to limit the movement range of the particles.

[0049] The method for constructing the three-dimensional discrete element model includes:

[0050] The base is a movable rigid base composed of spheres; the material of the lateral baffle is set to be a non-bonded material, and the particles are given the properties of a bonded material;

[0051] The particle radius, Young's modulus, Poisson's ratio, density, friction coefficient, tensile strength, shear strength, gravitational acceleration, time step safety factor, and local damping coefficient are set, and the particles are continuously deposited under the action of gravity to form a discrete element initial model;

[0052] According to the lithologic characteristics of different strata in the study area, the mechanical layering information of different thicknesses is converted into discrete element information in combination with the stratum histogram;

[0053] For the discrete element initial model, rock mechanical parameters of different strengths of the corresponding strata are assigned, wherein the rock mechanical parameters at least include friction coefficient, tensile strength and shear strength, and a three-dimensional discrete element model with layered characteristics is constructed.

[0054] In one embodiment, S103 , different extrusion conditions are set according to the three-dimensional discrete element model to perform iterative calculations to determine simulation results to be analyzed.

[0055] Specifically, different extrusion rates and extrusion distances are set for discrete element numerical simulation. For the three-dimensional discrete element model, the particle walls on the left and right sides are divided into independent areas along the central axis, and different horizontal extrusion rates are applied on both sides of the model; horizontal extrusion rates are applied to the left and right areas respectively, and the horizontal extrusion rates on the left and right sides are parallel and in opposite directions; by setting the extrusion distance, the degree and range of shear deformation are controlled, and after iterative calculations, the three-dimensional strike-slip fault simulation results to be analyzed are completed.

[0056] In one embodiment, S104, the simulation results are compared with the analysis results of the 3D seismic interpretation to determine the influence of different mechanical layers on the deformation characteristics of the strike-slip fault. The specific process includes:

[0057] Slice the simulation results of the 3D discrete element model, select the preset section for segmentation, and obtain the strike-slip deformation and fracture results of different strata;

[0058] The strike-slip deformation and fracture results were compared with the analysis results of 3D seismic interpretation to determine the control effect of stratification characteristics on strike-slip fracture behavior and obtain the influence of different mechanical stratification characteristics on the formation and evolution of strike-slip fractures.

[0059] The following is an explanation of the method for simulating the layered deformation process of strike-slip faults based on three-dimensional discrete elements proposed in the present invention with reference to an exemplary scenario.

[0060] In step 1, the morphological characteristics of strike-slip faults are described in detail by combining seismic interpretation profiles within the study area. In combination with drilling and logging data, the vertical variations, lithology, thickness, and physical and mechanical properties of the strata are analyzed to further determine the mechanical differences between the strata. Based on the mechanical properties of the strata, a stratum mechanical stratification pattern is established within the study area, providing a modeling foundation for subsequent 3D discrete element simulations.

[0061] Step 2) Based on the divided formation mechanics stratification pattern, a discrete element model with clear stratification characteristics is constructed.

[0062] Specifically, the dimensions of the 3D simulation chamber were set to 120m in length, 260m in width, and 60m in height. Furthermore, a rectangular area was constructed, consisting of a base and multiple lateral baffles, to limit the range of particle movement.

[0063] In terms of material settings, the wall material is defined as an unbound material, while the particles are assigned the properties of a bounded material. Initial particles are then randomly and uniformly generated within the rectangular area and continuously deposited under the action of gravity to form the initial discrete element model. The particle information in the initial model is exported to the sample.txt file.

[0064] Create a simulation box and its boundary with the same size as the initial model, import the external particle data file sample.txt; set the microscopic parameters of the particle material, where the microscopic parameters include at least: particle radius, Young's modulus, Poisson's ratio, density, friction coefficient, tensile strength, shear strength, gravitational acceleration, time step safety factor, and local damping coefficient.

[0065] The microscopic parameters of the particles are shown in Table 1.

[0066] Table 1 Microscopic parameters of particles

[0067]

[0068] To investigate the impact of different mechanical stratifications on strike-slip fault deformation, different strata were set up in the discrete element model. In this example, gypsum and carbonate rocks were distinguished based on their mechanical properties. The rock mechanical parameters for the different strata are shown in Table 2. Specifically, gypsum and carbonate rocks exhibit significant differences in tensile and shear strength.

[0069] Table 2 Rock mechanical parameters of different layers

[0070] Strata Friction coefficient μ Tensile strength T(MPa) Shear strength C (MPa) Gypsum rock (red) 0.0 0 0 Carbonate rock (gray) 0.1 4.5 5

[0071] The specific steps of the layered method are:

[0072] First, define the material parameters and set the microscopic parameters corresponding to the mud salt and carbonate rock layers.

[0073] Secondly, based on the mechanical stratification pattern obtained in step 1), the height ranges of different strata are defined based on the total height of the model; wherein the upper and lower limits of the strata can be determined by the proportional coefficient.

[0074] Finally, all particles are traversed and their positions in space are used to determine whether they fall within the height range and horizontal range of a certain layer. If the conditions are met, the material properties of the particle are set to the material of the corresponding layer, and particles of different layers are distinguished by color. Figure 2 , which is a schematic diagram of a three-dimensional discrete element initial model of different mechanical layering characteristics according to a specific embodiment of the present invention.

[0075] Step 3), set different extrusion conditions to perform iterative calculations to complete the simulation results to be analyzed.

[0076] Specifically, along the central axis, the particle walls on the left and right sides are divided into independent areas, and different horizontal extrusion speeds are applied on both sides of the model.

[0077] The horizontal extrusion speed V1 = (0.0, 0.5, 0.0) m·s is applied to the left area -1 The right side area is subjected to a horizontal extrusion velocity V2=(0.0,-0.5,0.0)m·s, which is parallel to but in the opposite direction. -1 At the same time, the extrusion distance is set to 2.6m to accurately control the degree and range of shear deformation. Then, continuous iterative calculations are performed to complete the three-dimensional strike-slip fault simulation results to be analyzed. Figure 3 FIG. 1 is a schematic diagram of a strike-slip fracture shear process according to a specific embodiment of the present invention; Figure 4 FIG. 1 is a schematic diagram of three-dimensional strike-slip fracture shear strain according to a specific embodiment of the present invention.

[0078] Step 4) Slice the discrete element model results to obtain the strike-slip deformation and fracture results of different strata.

[0079] Specifically, the 3D discrete element model is sliced ​​and segmented by selecting specific sections to analyze the strike-slip deformation and fault behavior in different strata. Figure 5 The figure shows a schematic diagram of a shear strain cross section of a strike-slip fault according to a specific embodiment of the present invention. The cross section is compared with the seismic cross section to reveal the influence of mechanical layering characteristics on the formation and evolution of strike-slip faults.

[0080] The proposed method for simulating the stratified deformation process of strike-slip faults based on three-dimensional discrete elements fully considers the influence of different mechanical stratifications on the deformation results of strike-slip faults. Based on seismic interpretation profiles and stratigraphic histograms, a discrete element model with clear stratification characteristics was constructed. By setting different microscopic parameters and model boundary motion parameters, the strike-slip shear process under different mechanical stratifications in the study area was simulated. The simulation results were compared with the analysis results of three-dimensional seismic interpretation to reveal the dominant role of stratification characteristics on strike-slip fault behavior, providing experimental basis and theoretical reference for oil and gas exploration controlled by deep strike-slip faults.

[0081] It should be noted that although the operations of the method of the present invention are described in a specific order in the above embodiments and drawings, this does not require or imply that these operations must be performed in this specific order, or that all illustrated operations must be performed to achieve the desired results. Additionally or alternatively, certain steps may be omitted, multiple steps may be combined into one step, and / or one step may be broken down into multiple steps.

[0082] After introducing the method of the exemplary embodiment of the present invention, next, reference is made to Figure 6 A device for simulating the delamination deformation process of a strike-slip fault based on three-dimensional discrete elements according to an exemplary embodiment of the present invention is introduced.

[0083] The implementation of the device for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements can be referred to in the implementation of the above-mentioned method, and the repeated parts will not be repeated here. The terms "module" or "unit" used below can refer to a combination of software and / or hardware that implements a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, implementation using hardware, or a combination of software and hardware, is also possible and contemplated.

[0084] Based on the same inventive concept, the present invention also proposes a device for simulating the layered deformation process of strike-slip faults based on three-dimensional discrete elements, such as Figure 6 As shown, the device includes:

[0085] Mechanical layer analysis module 610 is used to analyze the planar distribution characteristics of different strike-slip faults in the study area and, by combining drilling data and logging data, to obtain the formation mechanical layer pattern;

[0086] A three-dimensional discrete element model construction module 620 is used to assign different microscopic parameters to the discrete element model according to the formation mechanics layering pattern, so as to construct a three-dimensional discrete element model with layered characteristics;

[0087] The discrete element numerical simulation module 630 is used to set different extrusion conditions according to the three-dimensional discrete element model to perform iterative calculations and determine the simulation results to be analyzed;

[0088] The impact analysis module 640 is used to compare the simulation results with the analysis results of the 3D seismic interpretation to determine the impact of different mechanical layers on the deformation characteristics of the strike-slip fault.

[0089] In one embodiment, the mechanical delamination analysis module 610 is specifically used to:

[0090] The planar distribution characteristics of different strike-slip faults in the study area were analyzed using seismic interpretation profiles, and the planar and cross-sectional deformation characteristics of deep faults were obtained through analysis.

[0091] Acquire drilling data and logging data, and analyze vertical changes, lithology, thickness, and physical and mechanical properties of the formation based on the planar and cross-sectional deformation characteristics and in combination with the drilling data and logging data to determine mechanical differences in the formation;

[0092] According to the mechanical differences of the strata, the stratigraphic mechanical stratification pattern in the study area was divided.

[0093] In one embodiment, the three-dimensional discrete element model building module 620 is specifically used to:

[0094] A three-dimensional discrete element model is constructed according to the size, base and lateral boundaries of the preset three-dimensional simulation box; wherein the rectangular area enclosed by the base and the lateral boundaries corresponding to the multiple lateral baffles is used to limit the movement range of the particles.

[0095] In one embodiment, the three-dimensional discrete element model building module 620 is further specifically configured to:

[0096] The base is a movable rigid base composed of spheres; the material of the lateral baffle is set to be a non-bonded material, and the particles are given the properties of a bonded material;

[0097] The particle radius, Young's modulus, Poisson's ratio, density, friction coefficient, tensile strength, shear strength, gravitational acceleration, time step safety factor, and local damping coefficient are set, and the particles are continuously deposited under the action of gravity to form a discrete element initial model;

[0098] According to the lithologic characteristics of different strata in the study area, the mechanical layering information of different thicknesses is converted into discrete element information in combination with the stratum histogram;

[0099] For the discrete element initial model, rock mechanical parameters of different strengths of the corresponding strata are assigned, wherein the rock mechanical parameters at least include friction coefficient, tensile strength and shear strength, and a three-dimensional discrete element model with layered characteristics is constructed.

[0100] In one embodiment, the discrete element numerical simulation module 630 is specifically used to:

[0101] Different extrusion rates and extrusion distances were set for discrete element numerical simulation. For the three-dimensional discrete element model, the particle walls on the left and right sides were divided into independent regions along the central axis, and different horizontal extrusion rates were applied on both sides of the model. Horizontal extrusion rates were applied to the left and right regions, respectively, and the horizontal extrusion rates on the left and right sides were parallel and in opposite directions. By setting the extrusion distance, the degree and range of shear deformation were controlled. After iterative calculations, the three-dimensional strike-slip fault simulation results to be analyzed were completed.

[0102] In one embodiment, the impact analysis module 640 is specifically configured to:

[0103] Slice the simulation results of the 3D discrete element model, select the preset section for segmentation, and obtain the strike-slip deformation and fracture results of different strata;

[0104] The strike-slip deformation and fracture results were compared with the analysis results of 3D seismic interpretation to determine the control effect of stratification characteristics on strike-slip fracture behavior and obtain the influence of different mechanical stratification characteristics on the formation and evolution of strike-slip fractures.

[0105] It should be noted that while the detailed description above mentions several modules of the 3D discrete element-based strike-slip fault layered deformation simulation device, this division is merely exemplary and not mandatory. In practice, depending on the embodiments of the present invention, the features and functions of two or more modules described above can be embodied in a single module. Conversely, the features and functions of a single module described above can be further divided and embodied by multiple modules.

[0106] Based on the above invention concept, Figure 7 As shown, the present invention also proposes a computer device 700, including a memory 710, a processor 720, and a computer program 730 stored in the memory 710 and executable on the processor 720. When the processor 720 executes the computer program 730, the aforementioned method for simulating the layered deformation process of strike-slip faults based on three-dimensional discrete elements is implemented.

[0107] Based on the aforementioned inventive concept, the present invention proposes a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, it implements the aforementioned method for simulating the layered deformation process of strike-slip faults based on three-dimensional discrete elements.

[0108] Based on the aforementioned inventive concept, the present invention proposes a computer program product, which includes a computer program. When the computer program is executed by a processor, it implements a method for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements.

[0109] The method and device for simulating the layered deformation process of strike-slip faults based on three-dimensional discrete elements proposed in the present invention combine seismic interpretation profiles and drilling and logging data to accurately identify and describe the layered characteristics of different strata, thereby constructing a three-dimensional discrete element model with clear layered characteristics. By applying different extrusion speeds and extrusion distances to the discrete element model, the strike-slip shear process under different mechanical layering is simulated. This method can not only analyze the influence of each mechanical layering on the strike-slip fault deformation process, but also reveal the dominant role of layering characteristics in fault behavior, providing a new method for understanding the influence of different mechanical differences on strike-slip fault behavior. Comparing the simulation result profile with geophysical data further verified the key role of layering characteristics in the strike-slip fault process, and provided a strong experimental basis and theoretical reference for oil and gas exploration controlled by deep strike-slip faults.

[0110] The acquisition, storage, use, and processing of data in the technical solution of this application comply with relevant laws and regulations.

[0111] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, apparatus, or computer program products. Thus, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0112] The present invention is described with reference to flowcharts and / or block diagrams of methods and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowcharts and / or block diagrams, as well as combinations of processes and / or blocks in the flowcharts and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowcharts and / or block diagrams. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0113] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0114] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0115] Finally, it should be noted that the above-described embodiments are only specific implementation methods of the present invention, which are used to illustrate the technical solutions of the present invention, rather than to limit them. The scope of protection of the present invention is not limited thereto. Although the present invention has been described in detail with reference to the above-described embodiments, those skilled in the art should understand that any person skilled in the art can modify or easily conceive of changes to the technical solutions described in the above-described embodiments within the technical scope disclosed by the present invention, or replace some of the technical features therein with equivalents. Such modifications, changes, or replacements do not deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should be included in the scope of protection of the present invention. Therefore, the scope of protection of the present invention shall be subject to the scope of protection of the claims.

Claims

1. A method for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements, characterized in that: The method includes: The planar distribution characteristics of different strike-slip faults in the study area were analyzed, and the formation mechanical stratification pattern was obtained by combining drilling data and logging data. According to the stratum mechanics stratification pattern, different microscopic parameters are assigned to the discrete element model to construct a three-dimensional discrete element model with stratification characteristics; According to the three-dimensional discrete element model, different extrusion conditions are set to perform iterative calculations to determine the simulation results to be analyzed; The simulation results are compared with the analysis results of 3D seismic interpretation to determine the influence of different mechanical layers on the deformation characteristics of strike-slip faults.

2. The method for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements according to claim 1, characterized in that: By analyzing the planar distribution characteristics of different strike-slip faults in the study area and combining drilling and logging data, the formation mechanics stratification pattern was obtained, including: The planar distribution characteristics of different strike-slip faults in the study area were analyzed using seismic interpretation profiles, and the planar and cross-sectional deformation characteristics of deep faults were obtained through analysis. Acquire drilling data and logging data, and analyze vertical changes, lithology, thickness, and physical and mechanical properties of the formation based on the planar and cross-sectional deformation characteristics and in combination with the drilling data and logging data to determine mechanical differences in the formation; According to the mechanical differences of the strata, the stratigraphic mechanical stratification pattern in the study area was divided.

3. The method for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements according to claim 1, characterized in that: According to the formation mechanics layering pattern, different microscopic parameters are assigned to the discrete element model to construct a three-dimensional discrete element model with layered characteristics, including: A three-dimensional discrete element model is constructed according to the size, base and lateral boundaries of the preset three-dimensional simulation box; wherein the rectangular area enclosed by the base and the lateral boundaries corresponding to the multiple lateral baffles is used to limit the movement range of the particles.

4. The method for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements according to claim 1, characterized in that: The method for constructing the three-dimensional discrete element model includes: The base is a movable rigid base composed of spheres; the material of the lateral baffle is set to be a non-bonded material, and the particles are given the properties of a bonded material; The particle radius, Young's modulus, Poisson's ratio, density, friction coefficient, tensile strength, shear strength, gravitational acceleration, time step safety factor, and local damping coefficient are set, and the particles are continuously deposited under the action of gravity to form a discrete element initial model; According to the lithologic characteristics of different strata in the study area, the mechanical layering information of different thicknesses is converted into discrete element information in combination with the stratum histogram; For the discrete element initial model, rock mechanical parameters of different strengths of the corresponding strata are assigned, wherein the rock mechanical parameters at least include friction coefficient, tensile strength and shear strength, and a three-dimensional discrete element model with layered characteristics is constructed.

5. The method for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements according to claim 1, characterized in that: According to the 3D discrete element model, different extrusion conditions are set for iterative calculation to determine the simulation results to be analyzed, including: Different extrusion rates and extrusion distances were set for discrete element numerical simulation. For the three-dimensional discrete element model, the particle walls on the left and right sides were divided into independent regions along the central axis, and different horizontal extrusion rates were applied on both sides of the model. Horizontal extrusion rates were applied to the left and right regions, respectively, and the horizontal extrusion rates on the left and right sides were parallel and in opposite directions. By setting the extrusion distance, the degree and range of shear deformation were controlled. After iterative calculations, the three-dimensional strike-slip fault simulation results to be analyzed were completed.

6. The method for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements according to claim 1, characterized in that: The simulation results are compared with the analytical results of 3D seismic interpretation to determine the influence of different mechanical layers on the deformation characteristics of strike-slip faults, including: Slice the simulation results of the 3D discrete element model, select the preset section for segmentation, and obtain the strike-slip deformation and fracture results of different strata; The strike-slip deformation and fracture results were compared with the analysis results of 3D seismic interpretation to determine the control effect of stratification characteristics on strike-slip fracture behavior and obtain the influence of different mechanical stratification characteristics on the formation and evolution of strike-slip fractures.

7. A device for simulating the layered deformation process of a strike-slip fault based on three-dimensional discrete elements, characterized in that: The device includes: The mechanical layer analysis module is used to analyze the planar distribution characteristics of different strike-slip faults in the study area and, combined with drilling data and logging data, to obtain the formation mechanical layering pattern; A three-dimensional discrete element model construction module is used to assign different microscopic parameters to the discrete element model according to the formation mechanics layering pattern, so as to construct a three-dimensional discrete element model with layered characteristics; The discrete element numerical simulation module is used to set different extrusion conditions according to the three-dimensional discrete element model for iterative calculation to determine the simulation results to be analyzed; The impact analysis module is used to compare the simulation results with the analysis results of the three-dimensional seismic interpretation to determine the impact of different mechanical layers on the deformation characteristics of the strike-slip fault.

8. The device for simulating the stratified deformation process of a strike-slip fault based on three-dimensional discrete elements according to claim 7, characterized in that: The mechanical delamination analysis module is specifically used for: The planar distribution characteristics of different strike-slip faults in the study area were analyzed using seismic interpretation profiles, and the planar and cross-sectional deformation characteristics of deep faults were obtained through analysis. Acquire drilling data and logging data, and analyze vertical changes, lithology, thickness, and physical and mechanical properties of the formation based on the planar and cross-sectional deformation characteristics and in combination with the drilling data and logging data to determine mechanical differences in the formation; According to the mechanical differences of the strata, the stratigraphic mechanical stratification pattern in the study area was divided.

9. The device for simulating the stratified deformation process of a strike-slip fault based on three-dimensional discrete elements according to claim 7, characterized in that: The three-dimensional discrete element model building module is specifically used for: A three-dimensional discrete element model is constructed according to the size, base and lateral boundaries of the preset three-dimensional simulation box; wherein the rectangular area enclosed by the base and the lateral boundaries corresponding to the multiple lateral baffles is used to limit the movement range of the particles.

10. The device for simulating the stratified deformation process of strike-slip faults based on three-dimensional discrete elements according to claim 7, characterized in that: The three-dimensional discrete element model building module is further specifically used for: The base is a movable rigid base composed of spheres; the material of the lateral baffle is set to be a non-bonded material, and the particles are given the properties of a bonded material; The particle radius, Young's modulus, Poisson's ratio, density, friction coefficient, tensile strength, shear strength, gravitational acceleration, time step safety factor, and local damping coefficient are set, and the particles are continuously deposited under the action of gravity to form a discrete element initial model; According to the lithologic characteristics of different strata in the study area, the mechanical layering information of different thicknesses is converted into discrete element information in combination with the stratum histogram; For the discrete element initial model, rock mechanical parameters of different strengths of the corresponding strata are assigned, wherein the rock mechanical parameters at least include friction coefficient, tensile strength and shear strength, and a three-dimensional discrete element model with layered characteristics is constructed.

11. The device for simulating the stratified deformation process of a strike-slip fault based on three-dimensional discrete elements according to claim 10, characterized in that: The discrete element numerical simulation module is specifically used for: Different extrusion rates and extrusion distances were set for discrete element numerical simulation. For the three-dimensional discrete element model, the particle walls on the left and right sides were divided into independent regions along the central axis, and different horizontal extrusion rates were applied on both sides of the model. Horizontal extrusion rates were applied to the left and right regions, respectively, and the horizontal extrusion rates on the left and right sides were parallel and in opposite directions. By setting the extrusion distance, the degree and range of shear deformation were controlled. After iterative calculations, the three-dimensional strike-slip fault simulation results to be analyzed were completed.

12. The device for simulating the stratified deformation process of a strike-slip fault based on three-dimensional discrete elements according to claim 7, characterized in that: The impact analysis module is specifically used to: Slice the simulation results of the 3D discrete element model, select the preset section for segmentation, and obtain the strike-slip deformation and fracture results of different strata; The strike-slip deformation and fracture results were compared with the analysis results of 3D seismic interpretation to determine the control effect of stratification characteristics on strike-slip fracture behavior and obtain the influence of different mechanical stratification characteristics on the formation and evolution of strike-slip fractures.

13. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the method according to any one of claims 1 to 6 is implemented.

14. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 6 is implemented.

15. A computer program product, characterized in that The computer program product comprises a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 6 is implemented.

Images