Method, device and equipment for analyzing causes of complex strike-slip fault system and storage medium

Abstract

This application relates to the field of petroleum geology and discloses a method, apparatus, equipment, and storage medium for analyzing the genesis of complex strike-slip fracture systems. The method includes dividing the area where the target strike-slip fracture system is located into multiple tectonic blocks, with fracture distribution characteristics within the same tectonic block tending to be consistent; constructing a regional tectonic kinematic model based on the boundary types and internal deformation characteristics of each tectonic block to characterize the quantitative geometric relationships and kinematic processes of interactions between the blocks, thereby revealing the genetic mechanism of the fracture system; constructing a discrete element numerical model based on the regional tectonic kinematic model to simulate the interaction process between the blocks, obtaining simulated geological characteristics of the target strike-slip fracture system; comparing the simulated geological characteristics with actual geological characteristics, and if the consistency of the comparison results reaches a preset standard, the genetic mechanism revealed by the regional tectonic kinematic model is confirmed to be correct. This method more accurately reveals the genetic mechanism of complex strike-slip fracture systems.

Description

Technical Field

[0001] The embodiments of the present invention relate to the field of petroleum geology technology, and in particular to a method, apparatus, equipment and storage medium for analyzing the genesis of complex strike-slip fracture systems. Background Technology

[0002] Deep strike-slip fault-controlled oil and gas reservoirs are an important target type for oil and gas exploration in the Tarim Basin. Their reservoir development and hydrocarbon accumulation are controlled by strike-slip fault systems. Compared with conventional oil and gas reservoirs, these reservoirs are characterized by greater burial depth, stronger heterogeneity, complex fault systems, and uneven reservoir distribution. Specifically, the spatial combination, activity phases, and reactivation processes of faults at multiple scales directly control fracture development, reservoir space, and hydrocarbon migration pathways. Therefore, accurately analyzing the formation mechanism of multi-system tract strike-slip fault systems is fundamental to the effective evaluation and favorable area prediction of deep fault-controlled oil and gas reservoirs.

[0003] However, existing research on strike-slip fault systems mainly uses faults as the unit of analysis, which makes it difficult to reveal the formation mechanism of complex strike-slip fault systems and restricts a deeper understanding of the development law and hydrocarbon accumulation mechanism of fault-controlled reservoirs. Summary of the Invention

[0004] The purpose of this invention is to provide at least one method, apparatus, device, and storage medium for analyzing the formation of complex strike-slip fracture systems. This invention can at least solve the technical problem that existing methods mainly use fracture as the analysis unit, making it difficult to reveal the formation mechanism of complex strike-slip fracture systems. It can at least achieve a more accurate analysis of the formation mechanism of complex strike-slip fracture systems.

[0005] To address the aforementioned technical problems, at least one embodiment of this application provides a method for analyzing the causes of complex strike-slip fracture systems, including: Based on the distribution characteristics of the target strike-slip fault system identified by the 3D seismic data, the area where the target strike-slip fault system is located is divided into multiple tectonic blocks, wherein the fault distribution characteristics within the same tectonic block tend to be consistent. Based on the boundary types and internal deformation characteristics of each of the structural blocks, a regional tectonic kinematic model is constructed. The regional tectonic kinematic model is used to characterize the quantitative geometric relationship and kinematic process of the interaction between each of the structural blocks, so as to reveal the causal mechanism of the target strike-slip fracture system. Based on the regional tectonic kinematic model, a discrete element numerical model is constructed, and the interaction process between the tectonic blocks is simulated using the discrete element numerical model to obtain the simulated geological characteristics of the target strike-slip fault system. The simulated geological features are compared with the actual geological features. If the consistency between the simulated geological features and the actual geological features reaches a preset standard, then the causal mechanism revealed by the regional tectonic kinematic model is confirmed to be correct.

[0006] At least one embodiment of this application also provides a device for analyzing the causes of complex strike-slip fracture systems, comprising: The partitioning module is used to divide the area where the target strike-slip fault system is located into multiple structural blocks based on the distribution characteristics of the target strike-slip fault system as determined by the three-dimensional seismic data, wherein the fault distribution characteristics within the same structural block tend to be consistent. The first construction module is used to construct a regional tectonic kinematic model based on the boundary type and internal deformation characteristics of each of the tectonic blocks. The regional tectonic kinematic model is used to characterize the quantitative geometric relationship and kinematic process of the interaction between each of the tectonic blocks, so as to reveal the causal mechanism of the target strike-slip fracture system. The second construction module is used to construct a discrete element numerical model based on the regional tectonic kinematic model, and to use the discrete element numerical model to simulate the interaction process between the tectonic blocks, so as to obtain the simulated geological characteristics of the target strike-slip fault system. The confirmation module is used to compare the simulated geological features with the actual geological features. If the consistency between the simulated geological features and the actual geological features reaches a preset standard, the causal mechanism revealed by the regional tectonic kinematic model is confirmed to be correct.

[0007] At least one embodiment of this application also provides an electronic device, including: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, the instructions being executed by the at least one processor to enable the at least one processor to perform the above-described method for analyzing the causes of complex strike-slip fracture systems.

[0008] At least one embodiment of this application also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for analyzing the causes of complex strike-slip fracture systems.

[0009] The method for analyzing the genesis of complex strike-slip fracture systems provided in this application, based on the division of fracture system domains, introduces tectonic blocks, establishes a regional tectonic kinematic model through block boundary identification and tectonic kinematic modeling, and constructs a numerical simulation model under geological constraints on this basis. This achieves a unified analysis of fracture geometry, tectonic evolution process, and dynamic mechanism, thereby systematically revealing the differential genesis mechanisms of multiple fracture system domains. This solves the technical problem that existing methods, which mainly use fractures as the analysis unit, are unable to reveal the formation mechanism of complex strike-slip fracture systems.

[0010] In some optional embodiments, the division of the area where the target strike-slip fault system is located into multiple tectonic blocks based on the distribution characteristics of the target strike-slip fault system identified by 3D seismic data includes: Based on the fracture geometry and activity phases obtained from the three-dimensional seismic data, the target strike-slip fracture system is divided into multiple fracture system domains; Based on the division results of the fracture system domain, and combined with the continuity and spatial separation relationship of the fracture, the area where the target strike-slip fracture system is located is divided into multiple structural blocks. Establish the correspondence between the fracture system domain and the structural block, wherein each fracture system domain corresponds to one or more structural blocks with similar deformation characteristics, and adjacent structural blocks are separated by fracture zones or structural transition zones.

[0011] In this embodiment, after clarifying the distribution characteristics of the target strike-slip fault system based on 3D seismic data, the tectonic blocks are divided. This allows for the precise division of the fault system domain, the reasonable division of the area into tectonic blocks, and the establishment of corresponding relationships. This helps to clearly understand the geological structure and provides a scientific basis for subsequent resource exploration, geological research, and other tasks.

[0012] In some optional embodiments, constructing a regional tectonic kinematic model based on the boundary types and internal deformation characteristics of each of the tectonic blocks includes: Based on at least one of the fracture orientation, fracture combination pattern and fracture density in the internal deformation characteristics of each structural block, determine the displacement direction and relative displacement between adjacent structural blocks. Based on the boundary type between each of the structural blocks, the constraint method corresponding to different boundaries is determined, wherein the constraint method is used to limit the displacement of adjacent structural blocks in the direction perpendicular to the boundary and / or parallel to the boundary. The displacement direction and the relative displacement amount are used as the motion boundary conditions of the corresponding structural blocks, and the constraint method is applied to the corresponding boundary to establish the kinematic relationship between adjacent structural blocks. By integrating the kinematic relationships between adjacent structural blocks, the kinematic model of the region is obtained.

[0013] In this embodiment, a regional tectonic kinematic model is constructed based on the boundary type and internal deformation characteristics of the tectonic blocks. This model can accurately determine the displacement direction and relative displacement of adjacent blocks, reasonably set boundary constraints, and effectively establish and integrate the kinematic relationships between blocks. This helps to deepen the understanding of the laws of regional tectonic movement and provides reliable model support for geological disaster prediction, resource development, and other purposes.

[0014] In some optional embodiments, constructing a discrete element numerical model based on the kinematic model of the region includes: Based on the kinematic model of the region and the boundary conditions of the motion of each of the structural blocks, an initial discrete element numerical model is constructed. Rock mechanics parameters and contact parameters are set in the initial discrete element numerical model. The rock mechanics parameters are used to characterize the mechanical properties of rock materials in different parts of the target strike-slip fracture system, and the contact parameters are used to characterize the mechanical response characteristics of the internal structure and the boundary fracture zone of the structural block. Based on the kinematic characteristics of each of the structural blocks, at least one of displacement boundaries, velocity boundaries, or stress boundaries is set on the boundary of the initial discrete element numerical model to drive the model to deform and move, thereby obtaining the discrete element numerical model.

[0015] In this embodiment, a discrete element numerical model is constructed based on a regional tectonic kinematic model. This model can accurately combine the motion boundary conditions to build an initial model, reasonably set the rock and contact parameters to truly reflect the geological characteristics, and scientifically set the boundary to drive the deformation motion of the model. This helps to more accurately simulate the dynamic process of the target strike-slip fault system and provides reliable support for in-depth exploration of geological tectonic evolution.

[0016] In some optional embodiments, comparing the simulated geological features with actual geological features includes: At least one of the following characteristics in the simulated geological features—fracture geometry, fracture distribution, fracture combination pattern, and fracture density distribution—is compared with the actual geological features.

[0017] In this embodiment, comparing multiple key elements in the simulated geological features with actual geological features can accurately verify the reliability of the simulation results, promptly identify differences between the simulation and reality, and help correct simulation parameters or models to make the simulation closer to the real geological conditions.

[0018] In some optional embodiments, the method further includes: If the consistency between the simulated geological features and the actual geological features does not meet the preset standard, the regional tectonic kinematic model is adjusted, the discrete element numerical model is adjusted based on the adjusted regional tectonic kinematic model, and simulation is performed again based on the adjusted discrete element numerical model until the consistency meets the preset standard, thereby obtaining the formation mechanism of the target strike-slip fault system.

[0019] In this embodiment, by continuously adjusting the regional tectonic kinematic model and discrete element numerical model and repeatedly simulating when the consistency between the simulation and the actual geological characteristics does not meet the preset standard, it is possible to gradually approximate the real geological conditions, effectively eliminate error interference, and accurately reveal the complex causal mechanism of the target strike-slip fault system. Attached Figure Description

[0020] One or more embodiments are illustrated by way of example with reference to the accompanying drawings, and these illustrative descriptions do not constitute a limitation on the embodiments.

[0021] Figure 1 This is a flowchart of a method for analyzing the causes of complex strike-slip fracture systems provided in one embodiment of this application. Figure 1 ; Figure 2 This is a flowchart of a method for analyzing the causes of complex strike-slip fracture systems provided in one embodiment of this application. Figure 2 ; Figure 3 This is a schematic diagram of the fracture system distribution characteristics and system domain division results in the study area provided by an embodiment of this application; Figure 4 This is a schematic diagram of the block division results of the study area provided in one embodiment of this application; Figure 5 This is a schematic diagram of the boundary type of the research area construction block provided in one embodiment of this application; Figure 6 This is a schematic diagram of a kinematic model of the study area tectonic blocks under the action of multiple phases of tectonic activity, provided in one embodiment of this application. Figure 7 This is a schematic diagram of a discrete element numerical simulation model of the study area provided in one embodiment of this application; Figure 8 This is a schematic diagram of the fracture evolution and development in discrete element numerical simulation provided in one embodiment of this application; Figure 9 This is a schematic diagram illustrating the formation mechanism of a multi-fracture domain system provided in one embodiment of this application; Figure 10 This is a schematic diagram of a complex strike-slip fracture system cause analysis device provided in another embodiment of this application. Detailed Implementation

[0022] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the various embodiments of this application will be described in detail below with reference to the accompanying drawings. However, those skilled in the art will understand that many technical details have been provided in the various embodiments of this application to help readers better understand this application. However, the technical solutions claimed in this application can be implemented even without these technical details and various changes and modifications based on the following embodiments. The division of the various embodiments below is for the convenience of description and should not constitute any limitation on the specific implementation of this application. The various embodiments can be combined with and referenced by each other without contradiction.

[0023] To facilitate understanding of the embodiments of this application, relevant content regarding strike-slip fracture systems will be introduced first.

[0024] Existing research on strike-slip fault systems mostly focuses on single fault zones or local tectonic units, assuming that the fault system in the study area has a relatively consistent tectonic background and evolution process. Analysis is conducted within a unified stress field or a single tectonic stage, centering on the fault itself. Related studies primarily focus on the geometric and kinematic analyses of the faults, emphasizing the description of their geometric morphology, combination patterns, and formation mechanisms.

[0025] However, for complex strike-slip fault systems, especially in regions where multiple fault system domains (defined as one or more fault systems trending within a certain area, with differences in fault distribution characteristics and domain boundaries) coexist, different fault systems often exhibit significant differences in formation age, tectonic setting, and evolutionary processes. Existing research methods lack an analytical approach that considers the overall regional tectonic structure and integrates fault system domains with tectonic blocks, making it difficult to characterize the tectonic pattern of "block segmentation-boundary control-deformation coordination" and to explain the spatial differences in fault development from a dynamic perspective.

[0026] Regarding the formation mechanism and evolution analysis, existing studies reconstruct the periods of fault activity through seismic interpretation and stratigraphic correlation. However, an analytical framework that combines multi-stage tectonic activity with block kinematic processes is still lacking, making it difficult to systematically reveal the evolution process and dynamic mechanism of strike-slip fault systems under the superposition of multiple tectonic activities. While existing studies have introduced numerical simulation methods, they mostly employ idealized models for forward modeling, with insufficient consideration given to boundary conditions and block kinematic characteristics during model construction. Boundary condition settings often rely on empirical assumptions and lack constraints based on geological analytical results. Furthermore, the lack of effective connection between fault geometry, block division, and numerical simulation leads to discrepancies between simulation results and actual geological phenomena, making it difficult to quantitatively verify the genesis mechanism of faults.

[0027] Taking the northern and southern parts of the Tarim Basin as examples, the fault systems within different strike-slip fault system domains exhibit significant differences in formation age, development scale, geometric characteristics, tectonic setting, and evolutionary process. Current technologies still have shortcomings in unified analysis of multiple fault system domains, block boundary identification, kinematic modeling, numerical simulation, and comprehensive analysis of multi-stage tectonic evolution. There is an urgent need to establish a comprehensive analytical method for the genetic mechanisms of multi-system domain strike-slip fault systems, thereby overcoming the difficulty of existing methods in explaining the differential development and genetic mechanisms of faults within complex strike-slip fault zones.

[0028] To address the technical problem that existing methods, which primarily use fracture as the analysis unit, struggle to reveal the formation mechanism of complex strike-slip fracture systems, this invention proposes a method for analyzing the genesis of complex strike-slip fracture systems. The implementation details of this embodiment's method are described below. The following content is merely for ease of understanding and is not essential for implementing this solution.

[0029] Example 1: The method for analyzing the causes of complex strike-slip fracture systems in this embodiment can be applied to electronic devices with communication, computing, and data storage capabilities. The specific process can be as follows: Figure 1 As shown, it includes: Step 110: Based on the distribution characteristics of the target strike-slip fault system identified by the three-dimensional seismic data, the area where the target strike-slip fault system is located is divided into multiple structural blocks, wherein the fault distribution characteristics within the same structural block tend to be consistent. Specifically, the core of this step lies in using high-precision 3D seismic data to perform a detailed analysis of the geometric characteristics of the target strike-slip fault system, thereby providing an objective basis for the division of tectonic blocks. In practice, the planar distribution morphology, scale, and internal fault combination patterns of the strike-slip fault zone must first be identified based on information such as the reflection wave phase axis characteristics, seismic phase data, and coherence properties from the 3D seismic data.

[0030] Based on this, regions are divided according to the differences in fault distribution. Areas with statistically similar fault density, strike, combination morphology, and activity intensity are grouped into the same tectonic block, while areas where fault characteristics change abruptly are defined as boundaries between different blocks. For example, if within a region, faults mainly trend northeast and have similar fault combination patterns (such as en echelon), and these faults were active at similar times in geological history, then the areas affected by these faults can be classified as a single tectonic block.

[0031] This division breaks down the originally complex strike-slip fracture system into several independent units with relatively uniform internal deformation and clear boundary features, laying a refined geometric framework for revealing the interaction relationships between the blocks.

[0032] Step 120: Based on the boundary type and internal deformation characteristics of each of the structural blocks, a regional tectonic kinematic model is constructed. The regional tectonic kinematic model is used to characterize the quantitative geometric relationship and kinematic process of the interaction between each of the structural blocks, so as to reveal the causal mechanism of the target strike-slip fracture system. Specifically, after completing the division of tectonic blocks, a quantitative regional tectonic kinematic model needs to be constructed from the perspective of dynamic evolution to explain the phenomenon of "how these blocks move and lead to the current fault pattern".

[0033] The construction of the regional tectonic kinematic model requires the boundary types of each tectonic block (such as strike-slip fracture control boundary, thrust tectonic control boundary, or tectonic transition boundary, etc.) and the deformation characteristics inside the block (such as strike-slip deformation, compression deformation, extension deformation, or combined deformation, etc.) determined in step 110, and then the kinematic parameters such as the relative motion direction, displacement, and rotation angle between each block are set.

[0034] The kinematic model of a region's tectonic structure is not merely a simple geometric diagram, but a theoretical framework that quantitatively describes how interactions (such as compression, tension, or shear) between tectonic blocks regulate the total strain of the region. For example, it can be represented by the relative sliding of rigid blocks along strike-slip boundaries and the escape or rotation of adjacent blocks. By constructing this model, quantitative hypotheses about the formation mechanism of a target strike-slip fracture system can be proposed kinematically, providing clear boundary conditions and driving mechanisms for subsequent numerical simulations.

[0035] Step 130: Construct a discrete element numerical model based on the regional tectonic kinematic model, and use the discrete element numerical model to simulate the interaction process between the tectonic blocks to obtain the simulated geological characteristics of the target strike-slip fault system. Specifically, the Discrete Element Method (DEM) is a numerical simulation method based on discrete media mechanics theory. It treats geological bodies as a series of discrete particles or blocks that interact with each other through contact forces. The DEM can better simulate discontinuities in geological bodies, such as fractures and joints. When constructing a DEM, each structural block needs to be discretized into a certain number of particles or blocks, and contact parameters between the particles or blocks, such as contact stiffness and friction coefficient, need to be set.

[0036] When constructing a discrete element numerical model based on a regional tectonic kinematic model, it is necessary to accurately transform the geometry, boundary conditions, and kinematic characteristics of the tectonic blocks into parameters within the discrete element model. For example, based on the boundary type of the tectonic blocks, corresponding boundary constraints are set to restrict the motion of particles or blocks; based on the kinematic characteristics of the tectonic blocks, initial velocity or displacement boundary conditions are set to drive the model to generate deformation and motion. Simultaneously, it is also necessary to reasonably set rock mechanics parameters, such as elastic modulus and Poisson's ratio, to characterize the mechanical properties of rock materials in different parts of the target strike-slip fracture system.

[0037] By simulating the interaction between various tectonic blocks using a discrete element method (DEM) numerical model, the simulated geological characteristics of the target strike-slip fault system can be obtained, such as the fault geometry, distribution characteristics, combination patterns, and density distribution. Analysis of the simulation results provides a direct understanding of the formation process and evolution mechanism of the strike-slip fault system. For example, it allows observation of how faults gradually extend and connect under stress, and how the movement between different tectonic blocks affects the morphology and distribution of the faults.

[0038] Step 140: Compare the simulated geological features with the actual geological features. If the consistency between the simulated geological features and the actual geological features reaches a preset standard, then the causal mechanism revealed by the regional tectonic kinematic model is confirmed to be correct.

[0039] Specifically, when comparing simulated geological features with actual geological features, it is necessary to select appropriate comparison indicators and methods. For example, key indicators such as fault geometry, fault distribution characteristics, fault combination patterns, and fault density distribution characteristics can be selected for comparison. The comparison method can combine qualitative description and quantitative analysis. For instance, by drawing fault distribution maps, the location, strike, and extension length of simulated and actual faults can be visually compared; by calculating parameters such as fault density and dominant fault strike, the consistency between simulated and actual geological features can be quantitatively evaluated.

[0040] Preset criteria serve as the basis for judging the consistency between simulated and actual geological features. These criteria can be reasonably set according to the research objectives and the complexity of the actual geological conditions. For example, maximum permissible error for fault location and similarity threshold for fault strike can be set. When the consistency between simulated and actual geological features reaches the preset criteria, it indicates that the constructed regional tectonic kinematic model can reflect the actual geological conditions well, and the genetic mechanisms it reveals have high reliability.

[0041] In summary, the genetic analysis method for complex strike-slip fracture systems provided in this embodiment, based on the division of fracture system domains, introduces tectonic blocks, establishes a regional tectonic kinematic model through block boundary identification and tectonic kinematic modeling, and constructs a numerical simulation model under geological constraints on this basis. This achieves a unified analysis of fracture geometry, tectonic evolution process, and dynamic mechanism, thereby systematically revealing the differential genetic mechanisms of multiple fracture system domains. This addresses the technical problem of existing methods, which primarily use fractures as the analysis unit, making it difficult to reveal the formation mechanism of complex strike-slip fracture systems.

[0042] In an optional embodiment, the step of dividing the area where the target strike-slip fault system is located into multiple tectonic blocks based on the distribution characteristics of the target strike-slip fault system determined by the 3D seismic data includes: dividing the target strike-slip fault system into multiple fault system domains according to the fault geometric features and activity periods obtained from the 3D seismic data; dividing the area where the target strike-slip fault system is located into multiple tectonic blocks based on the division results of the fault system domains and the continuity and spatial separation of the faults; and establishing a correspondence between the fault system domains and the tectonic blocks, wherein each fault system domain corresponds to one or more tectonic blocks with similar deformation characteristics, and adjacent tectonic blocks are separated by fault zones or tectonic transition zones.

[0043] Specifically, traditional tectonic zoning typically relies solely on the planar distribution of faults (such as strike and density). However, this embodiment must also consider the active phases of faults (i.e., the chronological order of their formation or the stages of their activity). For example, in 3D seismic data, the formation and activity of faults can be determined by analyzing fault growth indices, variations in stratum thickness on either side of a fault, or the age of strata crossed by the fault. Faults that formed within the same tectonic period, possess similar geometric characteristics, and have a generative connection are grouped into the same "fault system domain." For instance, faults with similar strikes and the same active phase are grouped into a strike-slip fault system domain; faults with compressional characteristics and similar active ages are grouped into a thrust fault system domain.

[0044] The continuity of a fault reflects its spatial extension and integrity. In 3D seismic data, the continuity of a fault can be determined by tracing its reflection characteristics. The spatial separation relationship between faults is an important basis for dividing tectonic blocks. Faults in different fault system domains may intersect, truncate, or restrict each other, thus dividing the subsurface space into different regions. Based on the division of fault system domains, further consideration is given to the spatial continuity and spatial separation relationship of faults. Accordingly, the area where the target strike-slip fault system is located is divided into several "tectonic blocks".

[0045] The activity of each fault system domain induces specific deformation in the surrounding rocks, thus forming tectonic blocks with similar deformation characteristics. For example, the activity of a strike-slip fault system domain causes shear deformation in the surrounding rocks, forming structures such as shear zones and lenses; the activity of a thrust fault system domain causes compression deformation in the rocks, forming structures such as folds and thrust faults. By analyzing the deformation characteristics of tectonic blocks, their correspondence with fault system domains can be determined. Based on the results of deformation characteristic similarity analysis, the correspondence between fault system domains and tectonic blocks is established. Typically, each fault system domain corresponds to one or more tectonic blocks with similar deformation characteristics. For example, a strike-slip fault system domain may correspond to two adjacent tectonic blocks that have undergone similar shear deformation under the action of strike-slip faults.

[0046] It is worth noting that this application did not adopt the conventional approach of directly dividing tectonic blocks based on fault zones. Instead, it employed a process of "first dividing the fault system domains, then dividing the tectonic blocks, and finally establishing the corresponding relationships." The reason for this is that a "tectonic block" may have undergone multiple stages of deformation throughout its geological history. Its interior may contain "old fault system domains" that were formed early and later dissected by faults, as well as "new fault system domains" that formed concurrently with the block's boundary activity. If tectonic blocks are divided solely based on the current location of fault zones, it is easy to confuse all faults of different ages and origins within the block, making it impossible to distinguish which are inherent deformations within the block and which are products derived from the movement of the block's boundaries.

[0047] The structural block partitioning process in this embodiment aims to decouple structural events of different phases and scales, thereby more accurately defining the connotation of "structural block". By introducing a "time dimension" (fracture system domain), this embodiment solves the "overlap" problem in pure geometric partitioning, enabling subsequent kinematic models to distinguish between the active action of block boundaries and the passive response inside the block, thus more realistically simulating the cause of fracture systems.

[0048] In this embodiment, after clarifying the distribution characteristics of the target strike-slip fault system based on 3D seismic data, the tectonic blocks are divided. This allows for the precise division of the fault system domain, the reasonable division of the area into tectonic blocks, and the establishment of corresponding relationships. This helps to clearly understand the geological structure and provides a scientific basis for subsequent resource exploration, geological research, and other tasks.

[0049] In some optional embodiments, constructing a regional tectonic kinematic model based on the boundary type and internal deformation characteristics of each of the structural blocks includes: determining the displacement direction and relative displacement between adjacent structural blocks based on at least one of the fracture dominance orientation, fracture combination pattern, and fracture density in the internal deformation characteristics of each structural block; determining the constraint method corresponding to different boundaries based on the boundary type between each structural block, wherein the constraint method is used to limit the displacement of adjacent structural blocks in the direction perpendicular to the boundary and / or parallel to the boundary; using the displacement direction and the relative displacement as the kinematic boundary conditions of the corresponding structural blocks, and applying the constraint method to the corresponding boundaries to establish the kinematic relationship between adjacent structural blocks; and integrating the kinematic relationship between each adjacent structural block to obtain the regional tectonic kinematic model.

[0050] Specifically, the dominant fracture orientation within a structural block reflects the predominant direction of tectonic stress on that block. When adjacent structural blocks exhibit different dominant fracture orientations, their relative motion direction can be inferred. For example, if one block has a dominant fracture orientation of NE and the other NE, strike-slip movement may exist between them, and the angle between the strike-slip direction and the fracture orientation can help determine the specific displacement direction. For instance, if structural block A has a dominant fracture orientation of NE 45° and structural block B has a dominant fracture orientation of NE 30°, analysis suggests that strike-slip movement may exist between the two blocks, with the strike-slip direction roughly NE-SW.

[0051] Different fracture assemblages correspond to different tectonic deformation mechanisms. For example, en echelon fracture assemblages are usually associated with strike-slip motion. By analyzing the arrangement direction and length of en echelon fractures, the strike-slip displacement direction and relative displacement between adjacent blocks can be inferred. Broom-shaped fracture assemblages, on the other hand, may be related to rotational motion. Based on the convergence or divergence direction of the broom-shaped fractures, the rotation center and rotation angle can be determined, thus yielding the displacement direction and relative displacement. For instance, if en echelon fractures develop between tectonic blocks A and B, with a length of 5 km and an arrangement direction consistent with the inferred strike-slip direction, the relative strike-slip displacement between them can be estimated to be 3 km based on the characteristics of en echelon fractures.

[0052] Fracture density reflects the intensity of deformation within a structural block. Generally, areas with higher fracture density experience more intense deformation, and the relative displacement between adjacent blocks may be greater. By comparing the fracture densities of adjacent blocks and combining this with other deformation characteristics, the relative displacement between adjacent structural blocks can be determined more accurately. For example, structural block A has a fracture density of 5 fractures per square kilometer, while structural block B has a fracture density of 3 fractures per square kilometer. Combined with the previous analysis, this further confirms that the relative displacement between the two structural blocks is relatively large, consistent with the estimated 3 km.

[0053] The constraint method is determined based on the boundary type of the structural blocks. For example, the boundary between structural blocks A and B is a strike-slip fracture control boundary. For strike-slip fracture control boundaries, the displacement perpendicular to the boundary direction is set to zero, and only displacement parallel to the boundary direction is allowed. Simultaneously, based on the regional tectonic background and analogy with similar structures, the maximum strike-slip displacement is set to 5 km to limit the movement range of the blocks.

[0054] The determined displacement direction (northeast to southwest) and relative displacement (3 km) are used as the boundary conditions for the motion of structural block A and structural block B. The set constraints (zero displacement perpendicular to the boundary direction, maximum strike-slip displacement of 5 km) are applied to the strike-slip boundary. A kinematic model between adjacent structural blocks is established using numerical simulation software to simulate the relative motion of structural block A and structural block B under tectonic stress.

[0055] Furthermore, by combining the kinematic relationships between all adjacent blocks (for example, block A moves relative to block B at a speed and direction XX, with a strike-slip boundary and a normal displacement limit of 0; block A remains stationary relative to block C, etc.), a complete, closed, and self-consistent kinematic network is formed. This kinematic network is the final "regional tectonic kinematic model," which can quantitatively answer: If block X moves to the right, how will block Y move? What kind of internal strain will it produce? In this embodiment, a regional tectonic kinematic model is constructed based on the boundary type and internal deformation characteristics of the tectonic blocks. This model can accurately determine the displacement direction and relative displacement of adjacent blocks, reasonably set boundary constraints, and effectively establish and integrate the kinematic relationships between blocks. This helps to deepen the understanding of the laws of regional tectonic movement and provides reliable model support for geological disaster prediction, resource development, and other purposes.

[0056] In some optional embodiments, the construction of a discrete element numerical model based on the regional tectonic kinematic model includes: constructing an initial discrete element numerical model based on the regional tectonic kinematic model and the motion boundary conditions of each tectonic block; setting rock mechanics parameters and contact parameters in the initial discrete element numerical model, wherein the rock mechanics parameters are used to characterize the mechanical properties of rock materials in different parts of the target strike-slip fracture system, and the contact parameters are used to characterize the mechanical response characteristics of the interior and boundary fracture zones of the tectonic block; and setting at least one of displacement boundaries, velocity boundaries, or stress boundaries on the boundary of the initial discrete element numerical model according to the kinematic characteristics of each tectonic block to drive the model to deform and move, thereby obtaining the discrete element numerical model.

[0057] Specifically, firstly, based on the geometric shape, spatial distribution, and boundary types and kinematic relationships between the structural blocks determined by the regional tectonic kinematic model, an initial discrete element numerical model is established. This initial model typically consists of a large number of randomly or regularly arranged particle units, and the initial contact relationships between the particles must reflect the spatial separation characteristics of each structural block.

[0058] Based on this, two types of key parameters are further set in the model: one type is rock mechanics parameters, including the stiffness, bond strength, internal friction angle and density of particles inside different structural blocks, which are used to characterize the mechanical anisotropy and strength differences of rock materials in different parts of the target strike-slip fracture system (such as the core area of ​​rigid blocks, ductile deformation zones, fracture zones, etc.); the other type is contact parameters, mainly including the friction coefficient between particles and between particles and boundaries, normal and tangential stiffness, which are used to characterize the interaction between particles inside structural blocks and the mechanical response characteristics of boundary fracture zones under stress. For example, strike-slip boundaries can be set to a low friction coefficient to allow easy sliding.

[0059] Finally, based on the kinematic characteristics of each tectonic block determined by the regional tectonic kinematic model, such as displacement direction, relative displacement, and motion rate, appropriate boundary conditions are applied to the corresponding boundaries of the initial discrete element numerical model to drive the model to deform and move. These boundary conditions can be displacement boundaries that directly control the displacement, velocity boundaries that control the motion rate, or stress boundaries that simulate far-field tectonic stress. If necessary, multiple boundary conditions can be combined. For example, a constant velocity boundary can be applied to the active tectonic block to drive strike-slip motion, while a stress boundary can be applied to the passive tectonic block to simulate the resistance it experiences.

[0060] The discrete element numerical model established through the above steps can be used to simulate the interaction process between various structural blocks and reproduce the dynamic evolution of the strike-slip fracture system.

[0061] In this embodiment, a discrete element numerical model is constructed based on a regional tectonic kinematic model. This model can accurately combine the motion boundary conditions to build an initial model, reasonably set the rock and contact parameters to truly reflect the geological characteristics, and scientifically set the boundary to drive the deformation motion of the model. This helps to more accurately simulate the dynamic process of the target strike-slip fault system and provides reliable support for in-depth exploration of geological tectonic evolution.

[0062] In some optional embodiments, the method further includes: if the consistency between the simulated geological features and the actual geological features does not reach a preset standard, then adjusting the regional tectonic kinematic model, adjusting the discrete element numerical model based on the adjusted regional tectonic kinematic model, and performing simulation again based on the adjusted discrete element numerical model until the consistency reaches the preset standard, thereby obtaining the causal mechanism of the target strike-slip fault system.

[0063] Specifically, similarity thresholds and corresponding preset standards can be set for different geological features. For example, for morphological parameters such as the strike, dip, and dip angle of faults, higher similarity thresholds can be set, such as a strike similarity threshold of 90%, a dip similarity threshold of 85%, and a dip angle similarity threshold of 80%. The simulated fault morphological parameters must deviate from the actual measured values ​​within the corresponding threshold range. If the simulated fault strike deviates from the actual strike by more than 10% (i.e., failing to reach 90% similarity), the dip deviation exceeds 15%, and the dip angle deviation exceeds 20%, then the simulation of the fault morphological features is considered to have failed to meet the preset standards.

[0064] If the simulated fault strike deviates significantly from the actual strike, the analysis suggests that the assumptions regarding the direction of the regional tectonic stress field may be inaccurate. In this case, adjusting the directional parameters of the tectonic stress field in the regional tectonic kinematic model—for example, changing the originally assumed near-east-west stress field to a northeast-southwest direction—and recalculating the direction and manner of movement of the tectonic blocks will affect the fault strike.

[0065] Furthermore, based on the adjusted regional tectonic kinematic model, tectonic blocks and generated particle sets are redefined. Particle distribution and contact parameters near fault boundaries are optimized, for example, by increasing particle density at fault boundaries and adjusting contact stiffness and friction coefficient, to better simulate fault morphological development.

[0066] In this embodiment, by continuously adjusting the regional tectonic kinematic model and discrete element numerical model and repeatedly simulating when the consistency between the simulation and the actual geological characteristics does not meet the preset standard, it is possible to gradually approximate the real geological conditions, effectively eliminate error interference, and accurately reveal the complex causal mechanism of the target strike-slip fault system.

[0067] Example 2: Based on the above embodiments, this embodiment provides an application example of a method for analyzing the genesis of complex strike-slip fracture systems. For example... Figure 2 The diagram shown illustrates the flowchart of the genetic analysis method for complex strike-slip fracture systems in this embodiment. This method is based on the division of fracture system domains, introduces tectonic blocks, and establishes a regional tectonic kinematic model through block boundary identification and tectonic kinematic modeling. Based on this, a numerical simulation model under geological constraints is constructed to achieve a unified analysis of fracture geometry, tectonic evolution, and dynamic mechanisms, thereby systematically revealing the differential genetic mechanisms of multiple fracture system domains. The specific steps include the following: Step 1: Fracture System Analysis and System Domain Delineation ① Fracture analysis and geometric analysis Based on 3D seismic data, fracture analysis is carried out to identify the geometric information and distribution characteristics of major fractures, obtain parameters such as their strike, extension length, branching pattern and combination relationship, and analyze the spatial combination characteristics and distribution patterns between fractures. ② Determination of key tectonic phases of the fault Based on 3D seismic data, key strata were selected using a combination of planar and cross-sectional analysis to conduct fault analysis and stratigraphic interface identification. Based on the fault-stratigraphic contact relationship, the formation and main active phases of the fault system were analyzed and determined. ③ Division of fracture system domains and tectonic blocks The study area is divided into zones based on the geometric combination and activity characteristics of the fracture system. Regions with similar fracture combination patterns and activity characteristics are classified into the same fracture system domain. Based on the division results of the fracture system domain, combined with fracture continuity, combination characteristics and spatial separation relationship, the study area is divided into multiple tectonic blocks, clarifying the correspondence between the fracture system domain and the tectonic blocks, and realizing the transformation of research from "fracture" to "block".

[0068] Each fracture system domain corresponds to one or more tectonic blocks with relatively consistent deformation characteristics; adjacent fracture system domains are separated by fracture zones or tectonic transition zones. Step 2: Tectonic kinematic modeling of multi-fracture system domains ① Determining the structural properties of building blocks Based on fracture geometry and kinematic properties, the boundaries between different structural blocks are classified into different types: (1) Strike-slip fault control boundary; (2) Thrust structure control boundary; (3) Tectonic transition boundary (weak deformation zone).

[0069] Based on kinematic parameters, the internal fracture characteristics and overall deformation patterns of each structural block are analyzed. These parameters include: (1) Main deformation direction (fracture dominance and combination characteristics); (2) Deformation type (strike slip, compression, extension or composite); (3) Deformation strength (fracture density).

[0070] ② Constructing boundary conditions and establishing a kinematic model Based on the boundary types and internal deformation characteristics of tectonic blocks, a regional tectonic kinematic model is constructed, including: (1) Boundary conditions for the movement of each structural block (displacement direction and relative displacement magnitude); (2) Constraint methods corresponding to different boundary types (slip-slip control or compression control); (3) Overall deformation mode of multiple structural phases with the coordinated action of multiple blocks.

[0071] Step 3: Numerical simulation of the genetic mechanism of multi-fracture system tracts ① Construction of numerical simulation model Based on the established kinematic model and the boundary conditions of the tectonic blocks, a discrete element numerical simulation model is constructed to reasonably characterize the fracture system domain, tectonic blocks, and boundary conditions. Corresponding rock mechanics parameters and contact parameters are set to reflect the mechanical differences in different parts. According to the kinematic characteristics of the tectonic blocks, the model loading method (displacement boundary, velocity boundary, and stress boundary) is set. Different tectonic stages are simulated to realize the simulation of composite stress fields such as strike-slip and compression. ② Comparative analysis of numerical model results and geological features Discrete element numerical simulations are conducted to simulate the initiation, propagation, segmental connection, and penetration of fractures. The simulation results of fracture geometry, distribution characteristics, combination forms, and density distribution are compared and analyzed with the actual fracture development results to evaluate the consistency between the two. If the consistency is high, combine the construction analysis and numerical simulation of Step 2 to clarify the cause mechanism; if the consistency is insufficient, go back to adjust the construction kinematic model of Step 2, and carry out numerical simulation again based on the new model. ③ Verification of structural evolution model and determination of causal mechanism Based on the simulation-actual comparison results, the rationality of the structural kinematic model is verified, and the formation mechanism of the multi-system domain strike-slip fracture system is clarified, including the multilateral coaxial compression and the control of the tectonic regulation system on fracture development.

[0072] This embodiment uses the method to perform detailed seismic-geological analysis, tectonic block boundary identification and kinematic modeling, genetic mechanism analysis, and discrete element numerical simulation verification on the complex strike-slip fault system in the northern and southern parts of the Tarim Basin. Details are as follows: ① Based on 3D seismic data, fault analysis was conducted at the T74 interface in the northern and southern parts of the Tarim Basin, clarifying the geometric information and planar distribution characteristics of the main strike-slip fault system. Based on key seismic profiles, fault analysis and identification of key stratigraphic interfaces were carried out, analyzing fault-stratigraphic matching relationships and identifying the main active period of fault development and evolution as the Late Ordovician. Based on the geometric combination and activity characteristics of the faults, the fault system was divided into zones, with faults exhibiting similar distribution characteristics and combination patterns grouped into the same fault system domain. For example... Figure 3 The diagram shows the distribution characteristics of the fault system and the division of system domains in the study area. The study area is the Tarim Basin's northern and southern regions, which is also the area where the target strike-slip fault system is located in the above embodiments.

[0073] ② Based on the division results of three fault system domains in the northern and southern parts of the Tarim Basin, the study area is divided into three corresponding tectonic blocks. For example... Figure 4 The diagram shows a schematic representation of the tectonic block division in the study area. The types of tectonic block boundaries are identified, including thrust fault boundaries, strike-slip fault boundaries, and tectonic transition boundaries. Figure 5 The figure shown is a schematic diagram of the boundary types of the tectonic blocks in the study area.

[0074] ③ Based on the fracture characteristics and deformation modes within the block, determine its stress mechanism, construct a regional tectonic kinematic model, clarify the tectonic block to which each fracture system belongs, its boundary conditions and deformation characteristics, and deduce the kinematic connections between different fracture systems.

[0075] like Figure 6 The diagram shows a kinematic model of the tectonic blocks in the study area under multiple phases of tectonic activity in this embodiment. The diagram illustrates multiple tectonic blocks and their stress and displacement. Different blocks undergo displacement under compressive stress, with the direction of displacement indicated by arrows. The diagram also marks "Act 1," "Act 3," and other information indicating the tectonic activity phases. "Act 1" represents the earliest tectonic activity in the region. During this stage, under a specific stress field, the tectonic blocks began to displace, forming some early faults or folds, laying the foundation for subsequent tectonic evolution. "Act 3" represents a relatively later tectonic activity phase. During this period, the regional stress state may have changed, and the tectonic blocks continued to move under the new stress drive, modifying early structures or forming new tectonic features, making the regional geological structure more complex.

[0076] ④ Based on the Late Ordovician tectonic setting, tectonic block boundary conditions, and kinematic models, a paleogeomechanical numerical simulation model was constructed. The discrete element method was used to reproduce the development process of multi-system tract strike-slip faults, revealing their formation mechanism. For example... Figure 7 The figure shown is a schematic diagram of the discrete element numerical simulation model for the study area.

[0077] ⑤ Discrete element numerical simulation was conducted to simulate the entire process of fracture initiation, propagation, segmental connection, and penetration. The simulation results of fracture geometry, distribution characteristics, combination forms, and density distribution were compared with the actual fracture analytical results. The assessment concluded that the two were highly consistent. Figure 8 The figure shown is a diagram of fracture evolution and development in discrete element numerical simulation.

[0078] ⑥ For example Figure 9The diagram illustrates the genetic mechanism of the multi-fault domain system. It reveals the genetic mechanisms of three fault domains: Fault Domain 1 (TP39-TP29 in the north of Tarim), Fault Domain 2 (SB4-SB18 in the Shuntogor Low Uplift), and Fault Domain 3 (SB7-SB11 in the Shuntogor Low Uplift). Specifically, it shows the multi-stage activation mechanism dominated by conjugate shearing in the TP39-TP29 fault domain; the left-lateral strike-slip mechanism driven by thrust along the southern margin of the SB4-SB18 fault domain; and the right-lateral strike-slip mechanism dominated by thrust along the southwestern margin of the SB7-SB11 fault domain.

[0079] In summary, the method described in this embodiment was applied to conduct a systematic analysis of the complex strike-slip fault system developed within the Tarim Basin craton, achieving an integrated study of multi-fault system tract feature identification, tectonic evolution reconstruction, and genetic mechanism determination. Based on the fine-grained 3D seismic analysis of fault geometry and kinematic characteristics, three strike-slip fault system tracts were identified: TP39-TP29 in the northern Tarim Basin, SB4-SB18 and SB7-SB11 in the Shuntogole Low Uplift. Combined with fault-stratigraphic relationships, the multiple tectonic activity processes, including the Late Ordovician, were clarified. Furthermore, through boundary identification and kinematic modeling of tectonic blocks, a regional tectonic kinematic model was established. A discrete element numerical simulation model under geological constraints was further constructed, successfully reproducing the development and evolution process and spatial distribution characteristics of the multi-system tract strike-slip faults. The simulation results showed high consistency with the actual fault development characteristics, verifying the rationality of the tectonic evolution model.

[0080] Compared with the prior art, the present invention has the following advantages: 1. It achieves unified analysis of multiple fracture system domains, breaking through the limitations of traditional single fracture system research, and can effectively reveal the differences and co-evolutionary relationships between different fracture system domains; 2. By introducing structural block analysis units and boundary control mechanisms, an analysis framework of "block segmentation-boundary control-deformation coordination" is established to explain the spatial differences in fracture development from a dynamic perspective; 3. A numerical simulation method under geological constraints is constructed to achieve coupled analysis of fracture geometry, tectonic evolution process and genetic mechanism, overcoming the problems of high idealization and disconnection from actual geological conditions in existing numerical simulation models; 4. A quantitative analysis of the genetic mechanism of multiple fracture system domains was achieved, clarifying the formation mechanism of different fracture system domains under multiple tectonic processes.

[0081] The method described in this embodiment breaks through the traditional understanding of weak deformation within cratons, improves the accuracy and reliability of multi-system tract structural analysis of strike-slip faults, provides important technical support for the prediction of favorable areas of deep carbonate fault-controlled reservoirs and the deployment of oil and gas exploration, and has reference value for the study of strike-slip tectonic evolution in ancient craton basins around the world.

[0082] Example 3: Another embodiment of this application relates to a device for analyzing the genesis of complex strike-slip fracture systems. The implementation details of this embodiment's device are described below. The following details are provided for ease of understanding and are not essential for implementing this solution. A schematic diagram of this embodiment's device for analyzing the genesis of complex strike-slip fracture systems can be seen as follows: Figure 10 As shown, it includes a partitioning module 1010, a first construction module 1020, a second construction module 1030, and a confirmation module 1040.

[0083] The partitioning module 1010 is used to divide the area where the target strike-slip fault system is located into multiple structural blocks based on the distribution characteristics of the target strike-slip fault system as determined by the three-dimensional seismic data, wherein the fault distribution characteristics within the same structural block tend to be consistent. The first construction module 1020 is used to construct a regional tectonic kinematic model based on the boundary type and internal deformation characteristics of each of the tectonic blocks. The regional tectonic kinematic model is used to characterize the quantitative geometric relationship and kinematic process of the interaction between each of the tectonic blocks, so as to reveal the causal mechanism of the target strike-slip fracture system. The second construction module 1030 is used to construct a discrete element numerical model based on the regional tectonic kinematic model, and use the discrete element numerical model to simulate the interaction process between the tectonic blocks to obtain the simulated geological characteristics of the target strike-slip fault system. The confirmation module 1040 is used to compare the simulated geological features with the actual geological features. If the consistency between the simulated geological features and the actual geological features reaches a preset standard, then the causal mechanism revealed by the regional tectonic kinematic model is confirmed to be correct.

[0084] It is worth mentioning that all modules involved in this embodiment are logical modules. In practical applications, a logical unit can be a physical unit, a part of a physical unit, or a combination of multiple physical units. Furthermore, to highlight the innovative aspects of this application, this embodiment does not introduce units that are not closely related to solving the technical problems proposed in this application; however, this does not mean that other units are absent in this embodiment.

[0085] In some optional embodiments, the partitioning module includes: The first division unit is used to divide the target strike-slip fault system into multiple fault system domains based on the fault geometric features and activity periods obtained from the three-dimensional seismic data. The second division unit is used to divide the area where the target strike-slip fracture system is located into multiple structural blocks based on the division results of the fracture system domain and in combination with the continuity and spatial separation relationship of the fracture. The correspondence establishment unit is used to establish the correspondence between the fracture system domain and the structural block, wherein each fracture system domain corresponds to one or more structural blocks with similar deformation characteristics, and adjacent structural blocks are separated by fracture zones or structural transition zones.

[0086] In some optional embodiments, the first building module includes: The first analysis unit is used to determine the displacement direction and relative displacement between adjacent structural blocks based on at least one of the fracture dominance orientation, fracture combination pattern and fracture density in the internal deformation characteristics of each structural block. The second analysis unit is used to determine the constraint method corresponding to different boundaries based on the boundary type between each of the structural blocks, wherein the constraint method is used to limit the displacement of adjacent structural blocks in the direction perpendicular to the boundary or parallel to the boundary. A unit is established to use the displacement direction and the relative displacement as the motion boundary conditions of the corresponding structural block, and to apply the constraint method to the corresponding boundary to establish the kinematic relationship between adjacent structural blocks. An integration unit is used to integrate the kinematic relationships between adjacent structural blocks to obtain the kinematic model of the region.

[0087] In some alternative embodiments, the second building module includes: An initial model building unit is used to build an initial discrete element numerical model based on the kinematic model of the region and the motion boundary conditions of each of the building blocks. The parameter setting unit is used to set rock mechanics parameters and contact parameters in the initial discrete element numerical model. The rock mechanics parameters are used to characterize the mechanical properties of rock materials in different parts of the target strike-slip fracture system, and the contact parameters are used to characterize the mechanical response characteristics of the interior and boundary fracture zones of the structural block. The driving mode setting unit is used to set at least one of displacement boundaries, velocity boundaries, or stress boundaries on the boundaries of the initial discrete element numerical model according to the kinematic characteristics of each of the structural blocks, so as to drive the model to deform and move, thereby obtaining the discrete element numerical model. In some alternative embodiments, the apparatus further includes: An adjustment unit is used to adjust the regional tectonic kinematic model if the consistency between the simulated geological features and the actual geological features does not meet a preset standard, adjust the discrete element numerical model based on the adjusted regional tectonic kinematic model, and perform simulation again based on the adjusted discrete element numerical model until the consistency meets the preset standard, thereby obtaining the formation mechanism of the target strike-slip fault system.

[0088] Example 4: Another embodiment of this application relates to an electronic device, including: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, the instructions being executed by the at least one processor to enable the at least one processor to perform the complex strike-slip fracture system causal analysis method in the above embodiments.

[0089] In this embodiment, the memory and processor are connected via a bus, which can include any number of interconnected buses and bridges, connecting various circuits of one or more processors and the memory together. The bus can also connect various other circuits, such as peripheral devices, voltage regulators, and power management circuits, which are well known in the art and therefore will not be further described in this embodiment. A bus interface provides an interface between the bus and the transceiver. The transceiver can be a single element or multiple elements, such as multiple receivers and transmitters, providing a unit for communicating with various other devices over a transmission medium. Data processed by the processor is transmitted over the wireless medium via an antenna, which further receives data and transmits it to the processor.

[0090] The processor manages the bus and general processing, and also provides various functions, including timing, peripheral interfaces, voltage regulation, power management, and other control functions. Memory is used to store data used by the processor during operation.

[0091] Example 5: Another embodiment of this application relates to a computer-readable storage medium storing a computer program. When executed by a processor, the computer program implements the method embodiments described above.

[0092] That is, those skilled in the art will understand that all or part of the steps in the methods of the above embodiments can be implemented by a program instructing related hardware. This program is stored in a storage medium and includes several instructions to cause a device (which may be a microcontroller, chip, etc.) or processor to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0093] Those skilled in the art will understand that the above embodiments are specific embodiments for implementing this application, and in practical applications, various changes can be made to them in form and detail without departing from the spirit and scope of this application.

Claims

1. A method for analyzing the genesis of complex strike-slip fracture systems, characterized in that, include: Based on the distribution characteristics of the target strike-slip fault system identified by the 3D seismic data, the area where the target strike-slip fault system is located is divided into multiple tectonic blocks, wherein the fault distribution characteristics within the same tectonic block tend to be consistent. Based on the boundary types and internal deformation characteristics of each of the structural blocks, a regional tectonic kinematic model is constructed. The regional tectonic kinematic model is used to characterize the quantitative geometric relationship and kinematic process of the interaction between each of the structural blocks, so as to reveal the causal mechanism of the target strike-slip fracture system. Based on the regional tectonic kinematic model, a discrete element numerical model is constructed, and the interaction process between the tectonic blocks is simulated using the discrete element numerical model to obtain the simulated geological characteristics of the target strike-slip fault system. The simulated geological features are compared with the actual geological features. If the consistency between the simulated geological features and the actual geological features reaches a preset standard, then the causal mechanism revealed by the regional tectonic kinematic model is confirmed to be correct.

2. The method for analyzing the genesis of complex strike-slip fracture systems according to claim 1, characterized in that, Based on the distribution characteristics of the target strike-slip fault system identified by 3D seismic data, the area where the target strike-slip fault system is located is divided into multiple tectonic blocks, including: Based on the fracture geometry and activity phases obtained from the three-dimensional seismic data, the target strike-slip fracture system is divided into multiple fracture system domains; Based on the division results of the fracture system domain, and combined with the continuity and spatial separation relationship of the fracture, the area where the target strike-slip fracture system is located is divided into multiple structural blocks; Establish the correspondence between the fracture system domain and the structural block, wherein each fracture system domain corresponds to one or more structural blocks with similar deformation characteristics, and adjacent structural blocks are separated by fracture zones or structural transition zones.

3. The method for analyzing the genesis of complex strike-slip fracture systems according to claim 1, characterized in that, The construction of a regional tectonic kinematic model based on the boundary types and internal deformation characteristics of each of the aforementioned structural blocks includes: Based on at least one of the fracture orientation, fracture combination pattern and fracture density in the internal deformation characteristics of each structural block, determine the displacement direction and relative displacement between adjacent structural blocks. Based on the boundary type between each of the structural blocks, the constraint method corresponding to different boundaries is determined, wherein the constraint method is used to limit the displacement of adjacent structural blocks in the direction perpendicular to the boundary and / or parallel to the boundary. The displacement direction and the relative displacement amount are used as the motion boundary conditions of the corresponding structural blocks, and the constraint method is applied to the corresponding boundary to establish the kinematic relationship between adjacent structural blocks. By integrating the kinematic relationships between adjacent structural blocks, the kinematic model of the region's structure is obtained.

4. The method for analyzing the genesis of complex strike-slip fracture systems according to claim 3, characterized in that, The construction of a discrete-element numerical model based on the kinematic model of the region includes: Based on the kinematic model of the region and the boundary conditions of the motion of each of the structural blocks, an initial discrete element numerical model is constructed. Rock mechanics parameters and contact parameters are set in the initial discrete element numerical model. The rock mechanics parameters are used to characterize the mechanical properties of rock materials in different parts of the target strike-slip fracture system, and the contact parameters are used to characterize the mechanical response characteristics of the internal structure and the boundary fracture zone of the structural block. Based on the kinematic characteristics of each of the structural blocks, at least one of displacement boundaries, velocity boundaries, or stress boundaries is set on the boundary of the initial discrete element numerical model to drive the model to deform and move, thereby obtaining the discrete element numerical model.

5. The method for analyzing the genesis of complex strike-slip fracture systems according to claim 1, characterized in that, The comparison of the simulated geological features with the actual geological features includes: At least one of the following characteristics in the simulated geological features—fracture geometry, fracture distribution, fracture combination pattern, and fracture density distribution—is compared with the actual geological features.

6. The method for analyzing the genesis of complex strike-slip fracture systems according to claim 1, characterized in that, The method further includes: If the consistency between the simulated geological features and the actual geological features does not meet the preset standard, the regional tectonic kinematic model is adjusted, the discrete element numerical model is adjusted based on the adjusted regional tectonic kinematic model, and simulation is performed again based on the adjusted discrete element numerical model until the consistency meets the preset standard, thereby obtaining the formation mechanism of the target strike-slip fault system.

7. A device for analyzing the genesis of complex strike-slip fracture systems, characterized in that, include: The partitioning module is used to divide the area where the target strike-slip fault system is located into multiple structural blocks based on the distribution characteristics of the target strike-slip fault system as determined by the three-dimensional seismic data, wherein the fault distribution characteristics within the same structural block tend to be consistent. The first construction module is used to construct a regional tectonic kinematic model based on the boundary type and internal deformation characteristics of each of the tectonic blocks. The regional tectonic kinematic model is used to characterize the quantitative geometric relationship and kinematic process of the interaction between each of the tectonic blocks, so as to reveal the causal mechanism of the target strike-slip fracture system. The second construction module is used to construct a discrete element numerical model based on the regional tectonic kinematic model, and to use the discrete element numerical model to simulate the interaction process between the tectonic blocks, so as to obtain the simulated geological characteristics of the target strike-slip fault system. The confirmation module is used to compare the simulated geological features with the actual geological features. If the consistency between the simulated geological features and the actual geological features reaches a preset standard, the causal mechanism revealed by the regional tectonic kinematic model is confirmed to be correct.

8. The apparatus for analyzing the genesis of complex strike-slip fracture systems according to claim 7, characterized in that, The first building module includes: The first analysis unit is used to determine the displacement direction and relative displacement between adjacent structural blocks based on at least one of the fracture dominance orientation, fracture combination pattern and fracture density in the internal deformation characteristics of each structural block. The second analysis unit is used to determine the constraint method corresponding to different boundaries based on the boundary type between each of the structural blocks, wherein the constraint method is used to limit the displacement of adjacent structural blocks in the direction perpendicular to the boundary or parallel to the boundary. A unit is established to use the displacement direction and the relative displacement as the motion boundary conditions of the corresponding structural block, and to apply the constraint method to the corresponding boundary to establish the kinematic relationship between adjacent structural blocks. An integration unit is used to integrate the kinematic relationships between adjacent structural blocks to obtain the kinematic model of the region.

9. An electronic device, characterized in that, include: At least one processor; as well as, A memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor to enable the at least one processor to perform the genetic analysis method for complex strike-slip fracture systems as described in any one of claims 1 to 6.

10. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the method for analyzing the causes of complex strike-slip fracture systems as described in any one of claims 1 to 6.

Images