Fault slippage risk identification method and device, electronic equipment, storage medium and program product

Abstract

The embodiment of the invention provides a fault slippage risk identification method and device, electronic equipment, a storage medium and a program product. The method comprises the steps that a three-dimensional space grid model of a fault of a to-be-recognized target series of strata is established, the grid point occurrence of a grid on the fault of the to-be-recognized target series of strata is acquired, and the grid point occurrence comprises the dip angle and the trend of each grid point; obtaining a pore fluid pressure value of the grid surface of the grid; determining principal stress data on the grid surface of the grid according to a preset three-dimensional structure geologic model; determining a slip index value of the grid according to the occurrence of the grid, the pore fluid pressure value of the grid surface of the grid and the principal stress data of the grid; wherein the slip index value represents whether the grid surface of the grid slips or not; and determining the slippage risk coefficient of the fault of the to-be-identified target series of strata according to the number of the slippage grids. The method is used for achieving the effect of accurately determining the risk coefficient of slippage of the fault of the target series of strata.

Description

Technical Field

[0001] This application relates to the field of oil extraction technology, and in particular to a method, apparatus, electronic device, storage medium and program product for identifying fault slip risk. Background Technology

[0002] In the process of exploration and development of deep shale gas, it is necessary to determine the risk coefficient of fault slippage in the target strata to provide a reference for the design, construction and subsequent maintenance of deep shale gas wells. For example, the target strata are unconventional oil and gas reservoirs such as shale.

[0003] In existing technologies, the faults of the target stratigraphic system are subjected to overall data processing and analysis in order to determine the risk coefficient of the faults of the target stratigraphic system being slipped.

[0004] However, the above methods are not accurate in determining the risk factor of the fault slip of the target stratum. Therefore, there is an urgent need for a scheme that can accurately determine the risk factor of the fault slip of the target stratum. Summary of the Invention

[0005] This application provides a method, device, electronic device, storage medium, and program product for identifying fault slip risk. By integrating multi-source data such as geological, seismic, and well logging data, and combining three-dimensional stress field numerical simulation and three-dimensional spatial fault layer stress analysis, a quantitative method for identifying fault slip risk is established.

[0006] In a first aspect, embodiments of this application provide a method for identifying the slip risk of a fault, comprising:

[0007] A three-dimensional spatial mesh model of the fault of the target stratigraphic system to be identified is established. The three-dimensional spatial mesh model includes the attitude of the mesh points on the fault of the target stratigraphic system. The attitude of the mesh points includes the dip angle and orientation of the mesh points on the fault. The dip angle represents the dip angle of the mesh points on the mesh, and the orientation of the mesh points represents the orientation of the mesh points on the mesh. The pore fluid pressure value of the mesh surface is obtained. The pore fluid pressure value of the mesh surface represents the pressure borne by the fluid in the pores on the mesh surface.

[0008] Based on a preset three-dimensional structural geological model, the principal stress data on the grid surface are determined; wherein, the preset three-dimensional structural geological model represents the stratigraphic and fault information of the target strata to be identified, and the principal stress data represents the principal stress borne by the grid surface, including the maximum horizontal principal stress, the minimum horizontal principal stress, and the vertical principal stress of the grid.

[0009] Well logging is performed on the grid to obtain the directions of the horizontal maximum principal stress, the horizontal minimum principal stress, and the vertical principal stress of the grid.

[0010] The slip index value of the mesh is determined based on the orientation of the mesh points, the pore fluid pressure value of the mesh surface, the principal stress data of the mesh, the direction of the maximum horizontal principal stress, the direction of the minimum horizontal principal stress, and the direction of the vertical principal stress of the mesh; wherein, the slip index value characterizes whether the mesh surface will slip.

[0011] The slip risk value of the fault of the target stratigraphic system to be identified is determined based on the number of grids in which slip occurs; wherein the slip risk value characterizes the risk coefficient of slip of the fault of the target stratigraphic system to be identified.

[0012] In one possible implementation, the orientation of the grid points on the fault is the direction of the line connecting the grid point to the first adjacent grid point. The first adjacent grid point is the next grid point located on the same grid line as the grid point in the first direction, which is the horizontal direction in the grid division direction of the three-dimensional spatial grid model.

[0013] The slip index of the mesh is determined based on the mesh point orientation, the pore fluid pressure value of the mesh surface, the principal stress data of the mesh, the direction of the maximum horizontal principal stress, the direction of the minimum horizontal principal stress, and the direction of the vertical principal stress, including:

[0014] The angle between the orientation of the grid points and the due north direction is determined as the angle value of the grid.

[0015] Based on the grid point tilt angle and the included angle of the grid, determine the cosine value of the first included angle and the cosine value of the second included angle of the grid, and based on the grid point tilt angle, determine the cosine value of the third included angle of the grid; wherein, the first included angle of the grid is the angle between the normal to the grid surface and the direction of the maximum horizontal principal stress of the grid, the second included angle of the grid is the angle between the normal to the grid surface and the direction of the minimum horizontal principal stress of the grid, and the third included angle of the grid is the angle between the normal to the grid surface and the direction of the vertical principal stress of the grid.

[0016] The effective normal stress value of the mesh is determined based on the pore fluid pressure value of the mesh surface, the maximum horizontal principal stress of the mesh, the minimum horizontal principal stress of the mesh, the vertical principal stress of the mesh, the cosine value of the first included angle of the mesh, the cosine value of the second included angle of the mesh, and the cosine value of the third included angle of the mesh; wherein, the effective normal stress value of the mesh characterizes the effective normal stress borne by the mesh surface.

[0017] The shear stress value of the mesh is determined based on the maximum horizontal principal stress, the minimum horizontal principal stress, the vertical principal stress, the cosine of the first included angle, the cosine of the second included angle, and the cosine of the third included angle; wherein, the shear stress of the mesh characterizes the shear stress borne by the mesh surface.

[0018] The slip index of the mesh is determined based on the effective normal stress value and the shear stress of the mesh.

[0019] In one possible implementation, the cosine of the first included angle of the grid is cosβ1 = sinαsinγ. The cosine of the second included angle of the grid is cosβ2 = sinαcosγ; the cosine of the third included angle of the grid is cosβ1 = sinαsinγ; where α is the grid point tilt angle of the grid, and γ is the included angle value of the grid.

[0020] In one possible implementation, the effective normal stress value of the mesh is σ = (S i -P)cos 2 β1+(S h -P)cos 2 β2+(S v -P)cos 2 β3.

[0021] The shear stress of the mesh is

[0022]

[0023] Among them, S i S is the maximum horizontal principal stress of the mesh. h S is the minimum principal stress at the horizontal level of the mesh. v Let P be the vertical principal stress of the mesh, and P be the pore fluid pressure at the mesh surface. The slip index of the mesh is H = τ / (μσ + C); where σ is the effective normal stress of the mesh, τ is the shear stress of the mesh, μ is the preset friction coefficient of the mesh surface, and C is a preset constant.

[0024] In one possible implementation, a three-dimensional spatial mesh model of the fault of the target stratigraphic system to be identified is established, including:

[0025] Based on the three-dimensional seismic data volume, a three-dimensional spatial distribution feature is established; wherein, the three-dimensional seismic data volume includes the fault data of the target stratum to be identified; the three-dimensional spatial distribution feature characterizes the characteristics of the faults of the target stratum to be identified in three-dimensional space.

[0026] The three-dimensional spatial distribution features are divided into grids to obtain the grid of the fault of the target layer to be identified; and the connecting lines of the grids are determined; wherein, the connecting lines of the grids are the lines connecting a grid point to the second adjacent grid point adjacent to the grid point, and the second adjacent grid point adjacent to the grid point is the next grid point located on the same grid line as the grid point in a second direction, and the second direction is the perpendicular direction in the grid division direction in the three-dimensional spatial grid model.

[0027] The angle between the connecting line of the mesh and the physical horizontal line is determined as the tilt angle of the mesh point to obtain the three-dimensional spatial mesh model.

[0028] In one possible implementation, obtaining the pore fluid pressure value of the mesh surface includes:

[0029] Acquire well logging data, including formation velocities; and establish a relationship model based on preset measured pore fluid pressure values ​​and formation velocities in the well logging data, the model representing the correspondence between pore fluid pressure values ​​and formation velocities. Acquire three-dimensional seismic material data, and use the well logging data to perform seismic inversion calculations on the three-dimensional seismic material data to obtain the seismic layer velocities of the grid surfaces of the faults of the target strata to be identified.

[0030] The seismic layer velocity of the grid surface is used as the formation velocity and input into the relational model to obtain the pore fluid pressure value of the grid surface.

[0031] In one possible implementation, based on a preset three-dimensional geological model, the principal stress data on the grid surface are determined, including:

[0032] Acquire three-dimensional seismic material data of the target stratigraphic system to be identified, wherein the three-dimensional seismic material data includes stratigraphic interpretation data and fault interpretation data of the target stratigraphic system.

[0033] Based on the three-dimensional seismic material data, the preset three-dimensional structural geological model is established.

[0034] The mechanical parameters of the target layer to be identified are detected; wherein the mechanical parameters characterize the mechanical properties of the target layer.

[0035] The preset three-dimensional structural geological model is assigned the parameters from the mechanical parameters of the target strata to be identified, thus obtaining a three-dimensional mechanical model.

[0036] The three-dimensional mechanical model is meshed to obtain a mesh model for finite element analysis; wherein, the mesh model for finite element analysis is a three-dimensional mechanical model with a mesh.

[0037] The stress field numerical simulation is performed on the mesh model of the finite element analysis to obtain the initial principal stress data on the mesh surface.

[0038] Based on the pre-acquired core measured stress values ​​of the target strata to be identified, the initial principal stress data on the grid surface of the grid is corrected to obtain the principal stress data on the grid surface; wherein, the core measured stress value characterizes the principal stress borne by the core of the target strata.

[0039] Well logging is performed on the grid to obtain the directions of the horizontal maximum principal stress, the horizontal minimum principal stress, and the vertical principal stress of the grid, including:

[0040] Based on the pre-acquired imaging logging technology, the drilling-induced fracture phenomenon and wellbore collapse phenomenon on the fault of the target stratum to be identified are identified, and the directions of the horizontal maximum principal stress, the horizontal minimum principal stress and the vertical principal stress of the grid are obtained.

[0041] In one possible implementation, the method further includes:

[0042] If the slip risk value of the fault of the target stratum to be identified indicates that the risk coefficient of the fault of the target stratum to be identified is a high risk coefficient or a medium risk coefficient, a warning message is issued to the user; wherein, the warning message indicates that the fault of the target stratum to be identified will slip.

[0043] Secondly, embodiments of this application provide a device for identifying fault slip risk, comprising:

[0044] A module is established to create a three-dimensional spatial mesh model of the fault of the target stratigraphic system to be identified. The three-dimensional spatial mesh model includes the attitude of the mesh points on the fault of the target stratigraphic system to be identified. The attitude of the mesh points includes the dip angle of the mesh points on the fault and the orientation of the mesh points. The dip angle of the mesh points represents the dip angle of the mesh points on the mesh, and the orientation of the mesh points represents the orientation of the mesh points on the mesh.

[0045] The acquisition module is used to acquire the pore fluid pressure value of the mesh surface of the mesh, wherein the pore fluid pressure value of the mesh surface represents the pressure borne by the fluid in the pores on the mesh surface.

[0046] The first determining module is used to determine the principal stress data on the grid surface of the grid according to a preset three-dimensional structural geological model; wherein, the preset three-dimensional structural geological model represents the stratigraphic information and fault information of the target strata to be identified, and the principal stress data represents the principal stress borne by the grid surface, and the principal stress data includes the maximum horizontal principal stress of the grid, the minimum horizontal principal stress of the grid, and the vertical principal stress of the grid.

[0047] The logging module is used to log the grid to obtain the direction of the horizontal maximum principal stress, the direction of the horizontal minimum principal stress, and the direction of the vertical principal stress of the grid.

[0048] The second determining module is used to determine the slip index value of the grid based on the grid point orientation, the pore fluid pressure value of the grid surface, the principal stress data of the grid, the direction of the horizontal maximum principal stress, the direction of the horizontal minimum principal stress, and the direction of the vertical principal stress of the grid; wherein, the slip index value characterizes whether the grid surface will slip.

[0049] The third determining module is used to determine the slip risk value of the fault of the target stratum to be identified based on the number of grids that have slipped; wherein the slip risk value characterizes the risk coefficient of the fault of the target stratum to be identified slipping.

[0050] In one possible implementation, the orientation of the grid points on the fault is the direction of the line connecting the grid point to the first adjacent grid point. The first adjacent grid point is the next grid point located on the same grid line as the grid point in the first direction, which is the horizontal direction in the grid division direction of the three-dimensional spatial grid model.

[0051] The second determining module includes:

[0052] The first determining submodule is used to determine the angle between the orientation of the grid points and the due north direction, which is the angle value of the grid.

[0053] The second determining submodule is used to determine the cosine value of the first included angle and the cosine value of the second included angle of the mesh based on the inclination angle of the mesh points and the included angle value of the mesh, and to determine the cosine value of the third included angle of the mesh based on the inclination angle of the mesh points; wherein, the first included angle of the mesh is the angle between the normal of the mesh surface and the direction of the maximum horizontal principal stress of the mesh, the second included angle of the mesh is the angle between the normal of the mesh surface and the direction of the minimum horizontal principal stress of the mesh, and the third included angle of the mesh is the angle between the normal of the mesh surface and the direction of the vertical principal stress of the mesh.

[0054] The third determining submodule is used to determine the effective normal stress value of the mesh based on the pore fluid pressure value of the mesh surface, the maximum horizontal principal stress of the mesh, the minimum horizontal principal stress of the mesh, the vertical principal stress of the mesh, the cosine value of the first included angle of the mesh, the cosine value of the second included angle of the mesh, and the cosine value of the third included angle of the mesh; wherein, the effective normal stress value of the mesh characterizes the effective normal stress borne by the mesh surface of the mesh.

[0055] The fourth determination submodule is used to determine the shear stress value of the mesh based on the maximum horizontal principal stress, the minimum horizontal principal stress, the vertical principal stress, the cosine value of the first included angle, the cosine value of the second included angle, and the cosine value of the third included angle; wherein the shear stress of the mesh characterizes the shear stress borne by the mesh surface.

[0056] The fifth determining submodule is used to determine the slip index value of the mesh based on the effective normal stress value and the shear stress of the mesh.

[0057] In one possible implementation, the cosine of the first included angle of the grid is cosβ1 =

[0058] sinαsinγ; the cosine of the second included angle of the grid is cosβ2=sinαcosγ; the cosine of the third included angle of the grid is cosβ1=sinαsinγ; where α is the grid point tilt angle of the grid, and γ is the included angle value of the grid.

[0059] In one possible implementation, the effective normal stress value of the mesh is σ = (S i -P)cos 2 β1+(S h -P)cos 2 β2+(S v -P)cos 2 β3.

[0060] The shear stress of the mesh is

[0061] τ=√(S i -S h ) 2 cos 2 β1cos 2 β2+(S h -S v ) 2 cos 2 β2cos 2 β3+(S v -S i ) 2cos 2 β1cos 2 β3. Wherein, S i S is the maximum horizontal principal stress of the mesh. h S is the minimum principal stress at the horizontal level of the mesh. v Let P be the vertical principal stress of the mesh, and P be the pore fluid pressure at the mesh surface. The slip index of the mesh is H = τ / (μσ + C); where σ is the effective normal stress of the mesh, τ is the shear stress of the mesh, μ is the preset friction coefficient of the mesh surface, and C is a preset constant.

[0062] In one possible implementation, the establishment module includes:

[0063] The first submodule is used to establish three-dimensional spatial distribution features based on the three-dimensional seismic data volume; wherein, the three-dimensional seismic data volume includes data on faults of the target stratigraphic system to be identified; and the three-dimensional spatial distribution features characterize the features of the faults of the target stratigraphic system to be identified in three-dimensional space.

[0064] The sub-module is used to divide the three-dimensional spatial distribution features into a grid to obtain the grid of the fault of the target layer to be identified; and to determine the connecting lines of the grid; wherein, the connecting line of the grid is the line connecting the grid point to the second adjacent grid point adjacent to the grid point, and the second adjacent grid point adjacent to the grid point is the next grid point located on the same grid line as the grid point in the second direction, and the second direction is the perpendicular direction in the grid division direction in the three-dimensional spatial grid model.

[0065] The first generation submodule is used to determine the angle between the connecting line of the mesh and the physical horizontal line, which is the tilt angle of the mesh point, so as to obtain the three-dimensional spatial mesh model.

[0066] In one possible implementation, the acquisition module includes:

[0067] The acquisition submodule is used to acquire well logging data, which includes formation velocity; and to establish a relationship model based on the preset measured pore fluid pressure value and the formation velocity in the well logging data, wherein the relationship model is the correspondence between the pore fluid pressure value and the formation velocity.

[0068] The first processing submodule acquires three-dimensional seismic material data, uses the well logging data to perform seismic inversion calculations on the three-dimensional seismic material data, and obtains the seismic layer velocity of the grid surface of the fault of the target stratum to be identified.

[0069] The second processing submodule is used to input the seismic layer velocity of the grid surface as the formation velocity into the relational model to obtain the pore fluid pressure value of the grid surface.

[0070] In one possible implementation, the first determining module includes:

[0071] The second submodule is used to acquire three-dimensional seismic material data of the target stratigraphic system to be identified, wherein the three-dimensional seismic material data includes stratigraphic interpretation data and fault interpretation data of the target stratigraphic system; and to establish the preset three-dimensional structural geological model based on the three-dimensional seismic material data.

[0072] The detection submodule is used to detect the mechanical parameters of the target layer to be identified; wherein the mechanical parameters characterize the mechanical properties of the target layer.

[0073] The second generation submodule is used to assign the parameters of the mechanical parameters of the target strata to be identified to the preset three-dimensional structural geological model to obtain a three-dimensional mechanical model; and to perform meshing on the three-dimensional mechanical model to obtain a mesh model for finite element analysis; wherein, the mesh model for finite element analysis is a three-dimensional mechanical model with mesh.

[0074] The third processing submodule is used to perform stress field numerical simulation processing on the mesh model of the finite element analysis to obtain the initial principal stress data on the mesh surface.

[0075] The correction submodule is used to correct the initial principal stress data on the grid surface of the grid based on the pre-acquired core measured stress values ​​of the target stratum to be identified, so as to obtain the principal stress data on the grid surface; wherein, the core measured stress value represents the principal stress borne by the core of the target stratum.

[0076] The logging module is specifically used to: identify drilling-induced fracture phenomena and wellbore collapse phenomena on the faults of the target strata to be identified based on pre-acquired imaging logging technology, and obtain the direction of the horizontal maximum principal stress, the direction of the horizontal minimum principal stress, and the direction of the vertical principal stress of the grid.

[0077] In one possible implementation, the device further includes:

[0078] The warning module is used to issue a warning message to the user if the slip risk value of the fault of the target stratum to be identified indicates that the risk coefficient of the fault of the target stratum to be identified is a high risk coefficient or a medium risk coefficient; wherein the warning message indicates that the fault of the target stratum to be identified will slip.

[0079] Thirdly, embodiments of this application provide an electronic device, including: a memory and a processor;

[0080] The memory stores computer-executed instructions.

[0081] The processor executes computer execution instructions stored in the memory, causing the processor to perform the method provided as in the first aspect and any embodiment of the first aspect.

[0082] Fourthly, embodiments of this application provide a computer-readable storage medium storing computer-executable instructions, which, when executed by a processor, are used to implement the first aspect and / or various possible implementations of the first aspect.

[0083] Fifthly, embodiments of this application provide a computer program product, the computer program product including a computer program that, when executed by a processor, implements the methods provided in the first aspect and any embodiment of the first aspect. Attached Figure Description

[0084] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with this application and, together with the description, serve to explain the principles of this application.

[0085] Figure 1 A flowchart illustrating a method for identifying fault slip risk provided in this application. Figure 1 ;

[0086] Figure 2a A schematic diagram of a three-dimensional spatial grid structure of a fault in a target stratum to be identified, provided for the purposes of this application;

[0087] Figure 2b Left view of a schematic diagram of a three-dimensional spatial grid structure of a fault in a target stratigraphic system to be identified, provided in this application;

[0088] Figure 2c A front view of a three-dimensional spatial grid structure diagram of a fault of a target stratigraphic system to be identified, provided in this application;

[0089] Figure 3 A flowchart illustrating a method for identifying fault slip risk provided in this application (Figure 2);

[0090] Figure 4 A schematic diagram of the forces acting on a grid point of a grid in a three-dimensional spatial grid model of a fault of the target stratigraphy to be identified in this application;

[0091] Figure 5 A schematic diagram of the structure of a fault slip risk identification device provided in this application Figure 1;

[0092] Figure 6 Schematic diagram 2 of a fault slip risk identification device provided in this application;

[0093] Figure 7 A schematic diagram of an electronic device structure for identifying the slip risk of a fault, as provided in this application.

[0094] The accompanying drawings illustrate specific embodiments of this application, which will be described in more detail below. These drawings and descriptions are not intended to limit the scope of the concept in any way, but rather to illustrate the concept of this application to those skilled in the art through reference to particular embodiments. Detailed Implementation

[0095] Exemplary embodiments will now be described in detail, examples of which are illustrated in the accompanying drawings. When the following description relates to the drawings, unless otherwise indicated, the same numbers in different drawings denote the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with this application. Rather, they are merely examples of apparatuses and methods consistent with some aspects of this application as detailed in the appended claims.

[0096] In the development of deep shale gas, due to complex geological conditions, especially the interaction of factors such as geostress field, fracture occurrence and pore pressure, the problem of casing deformation in horizontal wells occurs frequently, which seriously affects the extraction efficiency and safety of shale oil and gas.

[0097] Currently, the likelihood of fault slip is primarily assessed by studying the influence of factors such as geostress field, fault orientation, and pore pressure on fault slip mechanisms. However, most of these studies remain at the qualitative analysis level, lacking quantitative identification methods, especially in identifying fault slip risks in three-dimensional space.

[0098] Therefore, the method, apparatus, equipment, storage medium, and program product for identifying fault slip risk provided in this application can accurately identify fault slip risk, provide a scientific basis for the design, construction, and subsequent maintenance of deep shale gas wells, and ensure the efficient and safe development of shale gas.

[0099] The technical solution of this application and how the technical solution of this application solves the above-mentioned technical problems are described in detail below with specific embodiments. These specific embodiments can be combined with each other, and the same or similar concepts or processes may not be described again in some embodiments. The embodiments of this application will now be described with reference to the accompanying drawings.

[0100] Figure 1A flowchart illustrating a method for identifying fault slip risk provided in this application. Figure 1 ,like Figure 1 As shown, the method includes:

[0101] S101. Establish a three-dimensional spatial mesh model of the fault of the target stratigraphy to be identified. The three-dimensional spatial mesh model includes the attitude of the mesh points on the fault of the target stratigraphy to be identified. The attitude of the mesh points includes the dip angle of the mesh points on the fault and the orientation of the mesh points. The dip angle of the mesh points represents the dip angle of the mesh points on the mesh, and the orientation of the mesh points represents the orientation of the mesh points on the mesh.

[0102] In one example, the target layer is, for example, a shale layer.

[0103] For example, establishing a three-dimensional spatial mesh model of the faults in the shale layer to be identified requires collecting geological data about the shale layer and its faults, including geological exploration reports, seismic interpretation data, drilling data, etc. This data will serve as the basis for establishing the three-dimensional spatial mesh model. Subsequently, the data undergoes preprocessing, including noise reduction, interpolation, and format conversion, to ensure data accuracy and consistency. In the preprocessed data, geological interpretation software or manual identification methods are used to determine the specific location, morphology, and strike of the shale layer faults. Based on the fault identification results and the spatial distribution characteristics of the shale layer, a three-dimensional mesh is created. The mesh division should consider the fault morphology, strike, and the thickness and dip angle of the shale layer to ensure that the mesh accurately reflects the spatial structure of the shale layer and its faults. The three-dimensional spatial mesh model includes the attitude of the grid points on the faults of the shale layer to be identified, where the attitude of the grid points mainly includes the dip angle and strike of the grid points on the fault. The grid point dip angle refers to the angle of inclination of each grid point relative to the physical horizontal line. It represents the dip angle of the rock strata at the location of the grid point. The dip angles of these grid points are assigned values ​​based on geological data and interpretation results. The grid point strike refers to the direction of extension of the fault plane projected onto the Earth's surface at a certain grid point on the fault.

[0104] S102. Obtain the pore fluid pressure value of the mesh surface. The pore fluid pressure value of the mesh surface represents the pressure borne by the fluid in the pores on the mesh surface.

[0105] For example, 3D seismic data, sonic logging data, and measured pore fluid pressure data are collected. Formation velocities at different depths are obtained through sonic logging data, which reflects the variation of formation velocity with depth. The raw data is processed, including outlier removal and filtering, to ensure accuracy and reliability.

[0106] Using sonic logging data, the relationship between formation porosity and parameters such as depth and sonic transit time can be established, thereby predicting formation pore fluid pressure. Sonic transit time, the time required for a sound wave to travel a certain distance through rock formations, is closely related to the porosity of the rock layer. Higher porosity results in slower sound wave propagation and a larger sonic transit time. Under normal compaction conditions, formation porosity gradually decreases with increasing formation depth, and sonic transit time decreases accordingly. However, in abnormally high-pressure formations, the increased pore fluid pressure leads to a relative increase in formation porosity, causing the sonic transit time to deviate from the normal trend. Formation velocity typically refers to the speed at which sound waves propagate through the formation; its reciprocal is the sonic transit time. Increased pore pressure alters the effective stress between rock particles, thus affecting the rock's elastic modulus and density, ultimately leading to changes in formation velocity. Generally speaking, an increase in pore pressure will reduce formation velocity because it reduces the effective stress between rock particles, making the rock more "soft" and slowing down the propagation speed of sound waves.

[0107] Detailed geological stratigraphic interpretation is performed on 3D seismic data volumes to clarify the spatial location and morphology of each stratigraphic layer. Well logging data is used as constraints to perform inversion processing on the seismic data, obtaining high-precision seismic layer velocity data. Seismic layer velocity data reflects the physical properties of rocks in the formation and is closely related to pore fluid pressure.

[0108] The seismic layer velocity data obtained from seismic inversion are substituted into the previously established pore fluid pressure-formation velocity relationship model to calculate the formation pore fluid pressure distribution in three-dimensional space. Based on the fault grid division results, the calculated three-dimensional pore fluid pressure data are distributed to each grid surface to obtain the pore fluid pressure value on each grid surface.

[0109] S103. Based on the preset three-dimensional structural geological model, determine the principal stress data on the grid surface. The preset three-dimensional structural geological model represents the stratigraphic and fault information of the target strata to be identified. The principal stress data represents the principal stress borne by the grid surface, including the maximum horizontal principal stress, the minimum horizontal principal stress, and the vertical principal stress of the grid.

[0110] For example, 3D seismic data, sonic logging data, and geological exploration reports are acquired to construct a 3D geological model. The 3D seismic data is then interpreted in detail to identify the spatial distribution characteristics of faults in the shale layer to be identified, including the fault's strike, dip angle, and variations in 3D space. Geological modeling software (such as Petrel and Gocad) is then used to integrate the seismic interpretation results and sonic logging data into the 3D geological model, constructing a refined 3D structural geological model. This model should accurately reflect the geometric morphology and positional relationships of the target strata and faults.

[0111] Core samples of the target shale layer were collected for rock mechanics experiments to obtain mechanical parameters such as Young's modulus and Poisson's ratio. Based on the differences in rock mechanical properties between the fault and the surrounding shale, mechanical parameters were assigned to the fault. Typically, the Young's modulus of the fault is taken as 40%–70% of that of the shale, and the Poisson's ratio is 0.02–0.1 higher than that of the shale, and can be further adjusted according to the fault grade. Based on the three-dimensional geological model, the fault was meshed. The mesh size should be determined according to the computational accuracy and resource requirements to ensure accurate reflection of stress changes along the fault plane.

[0112] The three-dimensional geological model, after being assigned mechanical parameters, is converted into a three-dimensional mechanical model for stress field numerical simulation. Stress boundary conditions in the study area are clarified using single-well experimental test data, including the direction and magnitude of the maximum and minimum horizontal principal stresses, and vertical stress. The mechanical model is meshed to establish a finite element analysis mesh model. The finite element method is used for stress field simulation and prediction; by solving the mechanical equations, the stress field distribution in three-dimensional space is calculated. The stress field simulation results are corrected using core-measured stress values ​​to ensure accuracy. The correction process may require multiple iterations until the error between the simulation results and measured values ​​is controlled within 10%.

[0113] Stress data were extracted from each grid surface from the three-dimensional stress field numerical simulation results. Based on stress tensor theory, the maximum horizontal principal stress, minimum horizontal principal stress, and vertical principal stress on each grid surface were calculated. These principal stress data characterize the main stress states experienced by the grid surface.

[0114] S104. Perform well logging on the grid to obtain the directions of the horizontal maximum principal stress, the horizontal minimum principal stress, and the vertical principal stress.

[0115] In one example, pre-acquired imaging logging data is used to identify drilling-induced fractures and wellbore collapse phenomena on shale faults through image analysis techniques. These phenomena are often related to the direction of the geostress field, particularly the orientation of drilling-induced fractures and the morphology of wellbore collapse, which can indirectly reflect the direction of the maximum horizontal principal stress. Drilling-induced fractures refer to fractures formed during drilling due to pressure changes within the wellbore or the intrusion of drilling fluid, causing the wellbore rock to break. The direction of these fractures is often closely related to the stress state around the wellbore, especially the direction of the maximum principal stress. Wellbore collapse refers to the collapse of the wellbore rock during drilling due to stress imbalance. The direction and pattern of wellbore collapse can also reflect the stress state around the wellbore, particularly the directions of the maximum and minimum horizontal principal stresses. The direction of the minimum principal stress can be inferred by being perpendicular to the direction of the maximum principal stress; the direction of the vertical principal stress is perpendicular to the plane determined by the directions of the maximum and minimum horizontal principal stresses, and the starting point of this vertical principal stress direction is located at the intersection of these orthogonal planes.

[0116] S105. Based on the orientation of the grid points, the pore fluid pressure value of the grid surface, the principal stress data of the grid, the direction of the maximum horizontal principal stress, the direction of the minimum horizontal principal stress, and the direction of the vertical principal stress, determine the slip index value of the grid. The slip index value characterizes whether the grid surface will slip.

[0117] For example, Figure 2a This is a schematic diagram of the three-dimensional spatial grid structure of the fault of the target layer to be identified; Figure 2b This is a left view of a schematic diagram of the three-dimensional spatial grid structure of the fault of the target layer to be identified; Figure 2c It is a front view of a schematic diagram of the three-dimensional spatial grid structure of the fault of the target layer to be identified.

[0118] like Figure 2aAs shown in the schematic diagram of the three-dimensional spatial grid structure of the fault of the target layer to be identified, the horizontal direction in the grid division direction of the three-dimensional spatial grid model is the first direction, that is, the direction of the line connecting grid points A, B, C, D, and E in the figure; the direction of the line connecting grid points K, P, G, S, and W in the figure; the direction of the line connecting grid points L, O, H, T, and X in the figure; and the line direction consistent with the direction of this connection. The vertical direction in the grid division direction of the three-dimensional spatial grid model is the second direction, that is, the direction of the line connecting grid points F, G, C, H, and I in the figure; the direction of the line connecting grid points Q, P, B, O, and N in the figure; the direction of the line connecting grid points R, S, D, T, and U in the figure; and the line direction consistent with the direction of this connection.

[0119] In one example, such as Figure 2b As shown in the left view of the three-dimensional spatial grid structure diagram of the fault of the target layer to be identified, point A is taken as the grid point of the study grid. The direction of grid point A is the direction of the line connecting grid point A to the next grid point (grid point B) in the first direction. The included angle of grid point A is the angle between the direction of grid point A and the due north direction, denoted by θ.

[0120] In one example, such as Figure 2c As shown in the front view of the three-dimensional spatial grid structure diagram of the fault of the target layer to be identified, point I is taken as the grid point of the study grid. The inclination angle of grid point I is the angle between the line connecting grid point I and the next grid point (grid point H) in the second direction and the physical horizontal line, denoted by α. The physical horizontal line refers to the line perpendicular to the direction of Earth's gravity, i.e., the horizontal line.

[0121] The pore fluid pressure value of the mesh surface is the pore fluid pressure value obtained in step S102, denoted by P. The principal stress data of the mesh, namely the maximum horizontal principal stress, the minimum horizontal principal stress, and the vertical principal stress, are the principal stress data obtained in step S103, wherein the maximum horizontal principal stress of the mesh is denoted by S. i The minimum horizontal principal stress is represented by S. h The vertical principal stress is represented by S. v The directions of the horizontal maximum principal stress, the horizontal minimum principal stress, and the vertical principal stress of the mesh are those obtained in step S104.

[0122] The cosine of the first included angle of the mesh is cosβ1 = sinαsinγ; the cosine of the second included angle of the mesh is cosβ2 = sinαcosγ; the cosine of the third included angle of the mesh is cosβ1 = sinαsinγ; where α is the tilt angle of the mesh points, and γ is the included angle value of the mesh; the effective normal stress of the mesh is σ = (S i -P)cos 2 β1+(S h -P)cos 2 β2+(S v -P)cos 2 β3; the shear stress of the mesh is The mesh slip index is H = τ / (μσ + C).

[0123] S106. Determine the slip risk value of the fault of the target stratum to be identified based on the number of grids that have slipped; wherein, the slip risk value represents the risk coefficient of the fault of the target stratum to be identified slipping.

[0124] For example, a slip risk value H is calculated for each grid on the fault of the shale layer to be identified. The H value reflects the probability of slippage on the grid surface. If H > 0, it indicates that the grid surface has a slip risk. The number of slip risk values ​​H for all grids on the fault of the shale layer to be identified is counted. If more than 75% of the grids on the fault of the shale layer to be identified have a slip risk, the fault of the shale layer to be identified is defined as a high slip risk fault; if the number of grids with slip risk on the fault of the shale layer to be identified is between 50% and 75%, the fault of the shale layer to be identified is defined as a medium slip risk fault; if the risk coefficient of slippage on the fault of the shale layer to be identified is a high risk coefficient or a medium risk coefficient, a warning message is issued to the user; the warning message indicates that the fault of the shale layer to be identified will slip.

[0125] This application provides a method for identifying the slip risk of faults. This method involves establishing a three-dimensional spatial grid model of the fault in the target strata, obtaining the included angle and dip angle of the grid points, acquiring the pore fluid pressure values ​​on the grid surfaces, and determining the principal stress data and directions on the grid surfaces. The slip index of each grid is calculated to obtain the slip risk coefficient of the fault in the target strata. If the slip risk coefficient of the fault in the target strata is high or medium, a warning message is issued to the user. This method can accurately identify the slip risk of faults in the target strata, providing a scientific basis for drilling operations and geological disaster prevention.

[0126] Figure 3 A flowchart illustrating a method for identifying fault slip risk provided in this application is shown in Figure 2. Figure 3 As shown, the method includes:

[0127] S301. Based on the three-dimensional seismic data volume, establish three-dimensional spatial distribution characteristics; wherein, the three-dimensional seismic data volume includes the fault data of the target strata to be identified; the three-dimensional spatial distribution characteristics characterize the fault features of the target strata to be identified in three-dimensional space.

[0128] For example, three-dimensional seismic data volumes are collected. These data volumes, acquired through seismic exploration techniques, contain reflected wave information from underground rock structures (including faults in shale layers). The collected three-dimensional seismic data volumes undergo preprocessing, including denoising, filtering, and amplitude recovery, to improve data quality, reduce interference signals, and make the characteristics of underground structures clearer. After data preprocessing, seismic data interpretation is performed. By comparing the characteristics of seismic waveforms (such as amplitude, frequency, and phase) and combining them with geological background knowledge, the location and morphology of the shale layers to be identified, as well as the faults within them, are identified. Fault identification is based on features such as the dislocation and distortion of seismic reflection phase axes. Based on the results of seismic data interpretation, geological modeling software or three-dimensional visualization technology is used to reconstruct and display the faults of the shale layers to be identified in three-dimensional space. By constructing a three-dimensional mesh model, the spatial morphology, positional relationships, and scale of the shale layers and their faults are described, forming a three-dimensional spatial distribution feature.

[0129] S302. The three-dimensional spatial distribution features are divided into grids to obtain the grid of the fault of the target layer to be identified; and the grid connection lines are determined; wherein, the grid connection line is the line connecting the grid point to the second adjacent grid point adjacent to the grid point, and the second adjacent grid point adjacent to the grid point is the next grid point located on the same grid line as the grid point in the second direction, and the second direction is the perpendicular direction in the grid division direction in the three-dimensional spatial grid model.

[0130] For example, a mesh is generated for the pre-constructed 3D spatial distribution features (i.e., the 3D morphology of the faults in the shale layer to be identified). The purpose of this step is to discretize the continuous 3D space into a series of regular or irregular mesh units for numerical analysis or visualization. After mesh generation, each mesh unit is defined by a mesh point, and a mesh point has six directly adjacent mesh points. The mesh connection line is the line connecting the mesh point to the next mesh point in the horizontal direction of the mesh division. The next mesh point is the mesh point to the left of the current mesh point in the horizontal direction. The horizontal direction is the second direction in the mesh division direction of the 3D spatial mesh model.

[0131] In one example, such as Figure 2a As shown, the line connecting grid point H is the line connecting grid point H and grid point C; the line connecting grid point C is the line connecting grid point C and grid point G; and the line connecting grid point G is the line connecting grid point G and grid point F.

[0132] S303. Determine the angle between the mesh connection line and the physical horizontal line, which is the mesh point tilt angle, to obtain a three-dimensional spatial mesh model.

[0133] For example, the tilt angle of a grid point is the angle between the line connecting the grid point and the physical horizontal line, where the physical horizontal line refers to a horizontal line perpendicular to the direction of Earth's gravity.

[0134] In one example, such as Figure 2c As shown, the inclination angle of grid point I is the angle between the connecting line IH of grid point I and the physical horizontal line of grid point I, i.e., α1; the inclination angle of grid point G is the angle between the connecting line FG of grid point G and the physical horizontal line of grid point G, i.e., α2.

[0135] S304. Acquire well logging data, which includes formation velocity; and establish a relationship model based on the preset measured pore fluid pressure value and the formation velocity in the well logging data. The relationship model is the correspondence between the pore fluid pressure value and the formation velocity.

[0136] For example, well logging data containing various physical properties of underground rock formations is acquired; among these, the main data includes formation velocities and the propagation speed of acoustic waves within the formation. Measured pore fluid pressure values ​​are obtained, which refer to the pressure of fluids (usually water, oil, or gas) in the pores of underground rock. The acquired data is processed and cleaned to ensure its quality and integrity. Then, necessary corrections are performed to eliminate the influence of instrument errors or environmental factors. A model relating the measured pore fluid pressure values ​​to formation velocities is established, using methods such as the Bowers method and the Fillippone method. The accuracy and reliability of the established model are cross-validated, and the model is adjusted and optimized based on the validation results.

[0137] S305. Acquire three-dimensional seismic material data, use well logging data to perform seismic inversion calculations on the three-dimensional seismic material data, and obtain the seismic layer velocity of the grid surface of the fault of the target stratum to be identified; use the seismic layer velocity of the grid surface as the formation velocity, input it into the relational model, and obtain the pore fluid pressure value of the grid surface.

[0138] For example, a refined geological stratigraphic interpretation is performed on the 3D seismic material data. "Refined" refers to using the same mesh generation method as the 3D spatial mesh model in step S101 to mesh the 3D seismic material data. Formation velocity information from well logging data is combined with the 3D seismic material data to provide additional constraints for seismic inversion. After completing the seismic inversion calculation, the seismic layer velocities of the fault mesh faces are obtained. Seismic inversion can be performed using pre-stack inversion methods such as seismic wave travel time inversion and tomographic inversion. The seismic layer velocities of the mesh faces are input as formation velocities into the previously established correspondence model between pore fluid pressure values ​​and formation velocities. This model is used to calculate the pore fluid pressure values ​​of the mesh faces, ultimately obtaining the pore fluid pressure values ​​of the fault mesh faces of the shale layer.

[0139] S306. Based on the preset three-dimensional structural geological model, determine the principal stress data on the grid surface; wherein, the preset three-dimensional structural geological model represents the stratigraphic and fault information of the target strata to be identified, and the principal stress data represents the principal stress borne by the grid surface, including the maximum horizontal principal stress, the minimum horizontal principal stress, and the vertical principal stress of the grid.

[0140] In one example, step S306 includes the following process:

[0141] Acquire three-dimensional seismic material data of the target stratigraphic system to be identified, including stratigraphic interpretation data and fault interpretation data of the target stratigraphic system; establish a preset three-dimensional structural geological model based on the three-dimensional seismic material data.

[0142] The mechanical parameters of the target layer to be identified are detected; where the mechanical parameters characterize the mechanical properties of the target layer.

[0143] The parameters of the mechanical parameters of the target strata to be identified are assigned to the preset three-dimensional structural geological model to obtain the three-dimensional mechanical model.

[0144] The three-dimensional mechanical model is meshed to obtain the mesh model for finite element analysis; the mesh model for finite element analysis is a three-dimensional mechanical model with a mesh.

[0145] The stress field numerical simulation of the mesh model from the finite element analysis is performed to obtain the initial principal stress data on the mesh surface.

[0146] Based on the pre-obtained core measured stress values ​​of the target strata to be identified, the initial principal stress data on the grid surface are corrected to obtain the principal stress data on the grid surface; among them, the core measured stress value characterizes the principal stress borne by the core of the target strata.

[0147] For example, three-dimensional seismic material data of the study area are acquired and organized, and then imported into professional geophysical interpretation software such as Petrel. In the software, the top and bottom interfaces of the shale layers to be identified are used to determine the stratigraphic interpretation by utilizing the amplitude, frequency, and other attributes of the seismic data. Similarly, by analyzing discontinuities and anomalies in the seismic data, faults are identified and interpreted, including their strike, dip, and displacement, thus completing fault interpretation. Based on the results of stratigraphic and fault interpretation, a three-dimensional structural geological model is constructed.

[0148] Shale cores were collected from the shale layer to be identified and processed into 50mm × 25mm plunger samples according to standards. Triaxial loading tests were performed on the plunger samples using the GCTS rock mechanics testing system to determine their Young's modulus and Poisson's ratio. The experimentally obtained mechanical parameters were assigned to the shale portion of the geological model. Based on the difference in mechanical properties between the fault and the surrounding shale, a relatively low Young's modulus and a relatively high Poisson's ratio were assigned to the fault, considering the influence of fault grade on the parameters.

[0149] The mechanical model is meshed using finite element analysis software (such as ANSYS, ABAQUS, etc.), ensuring that the mesh quality meets the analysis requirements, such as mesh size, shape, and distribution. This yields a mesh model suitable for finite element analysis.

[0150] Stress boundary conditions in the study area were determined using single-well experimental test data (such as well logging data and geostress testing). Boundary conditions were set in finite element analysis software, and a simulation was run to calculate the stress field distribution in three-dimensional space. The simulation results were output, including the horizontal maximum principal stress, horizontal minimum principal stress, and vertical principal stress data.

[0151] Measured stress values ​​from rock cores were obtained using methods such as rock acoustic emission and differential strain analysis. These measured stress values ​​were then compared with numerical simulation results. By adjusting parameters such as boundary loads or boundary displacements, the simulation results were made to match the measured values, with the error controlled within 10%, thus obtaining an accurate three-dimensional stress field distribution characteristic model.

[0152] S307. Perform well logging on the grid to obtain the directions of the maximum horizontal principal stress, the minimum horizontal principal stress, and the vertical principal stress of the grid.

[0153] In one example, step S307 includes the following process: based on the pre-acquired imaging logging technology, identify drilling-induced fracture phenomena and wellbore collapse phenomena on the fault of the target stratum to be identified, and obtain the direction of the horizontal maximum principal stress, the direction of the horizontal minimum principal stress, and the direction of the vertical principal stress of the grid.

[0154] For example, this step can refer to step S104 above, and will not be repeated here.

[0155] In one example, the directions of the horizontal maximum principal stress, the horizontal minimum principal stress, and the vertical principal stress of the grid are obtained through paleomagnetic and wave velocity anisotropy experiments. The paleomagnetic experiment, by analyzing the directions of ancient magnetic fields preserved in the rock and combining this with geological history, can infer the stress state and direction at the time of rock formation. The wave velocity anisotropy experiment, by measuring the differences in the propagation speed of sound waves in different directions within the rock, infers the stress distribution and principal stress directions within the rock. The direction of the minimum principal stress can be inferred by being perpendicular to the direction of the maximum principal stress; the direction of the vertical principal stress is perpendicular to the plane determined by the directions of the horizontal maximum and minimum principal stresses, and the starting point of this vertical principal stress direction is located at the intersection of these orthogonal planes.

[0156] In one example, by measuring the faults in the shale layer to be identified using multipole acoustic logging, the directions of the horizontal maximum principal stress, the horizontal minimum principal stress, and the vertical principal stress of the grid are obtained. Multipole acoustic logging is an advanced logging technique that analyzes the mechanical properties of the rock and the direction of the maximum principal stress by transmitting acoustic signals in multiple directions in the well and receiving the reflection and transmission of these signals in the rock. The direction of the minimum principal stress can be inferred by being perpendicular to the direction of the maximum principal stress; the direction of the vertical principal stress is perpendicular to the plane determined by the directions of the horizontal maximum and minimum principal stresses, and the starting point of the vertical principal stress direction is located at the intersection of these orthogonal planes.

[0157] S308. Based on the orientation of the grid points, the pore fluid pressure value of the grid surface, the principal stress data of the grid, the direction of the maximum horizontal principal stress, the direction of the minimum horizontal principal stress, and the direction of the vertical principal stress, determine the slip index value of the grid. The slip index value characterizes whether the grid surface will slip.

[0158] In one example, the grid point orientation includes the grid point orientation on the fault. The grid point orientation is the direction of the line connecting the grid point to the first adjacent grid point. The first adjacent grid point is the next grid point located on the same grid line as the grid point in the first direction, which is the horizontal direction in the grid division direction of the 3D spatial grid model. Step S208 includes the following process:

[0159] Determine the angle between the orientation of the grid points and the due north direction, which is the angle value of the grid.

[0160] Based on the grid point tilt angle and the grid's included angle, determine the cosine value of the first included angle and the cosine value of the second included angle of the grid, and based on the grid point tilt angle, determine the cosine value of the third included angle of the grid; where the first included angle of the grid is the angle between the normal to the grid surface and the direction of the maximum horizontal principal stress of the grid, the second included angle of the grid is the angle between the normal to the grid surface and the direction of the minimum horizontal principal stress of the grid, and the third included angle of the grid is the angle between the normal to the grid surface and the direction of the vertical principal stress of the grid;

[0161] The effective normal stress value of the mesh is determined based on the pore fluid pressure value of the mesh surface, the maximum horizontal principal stress of the mesh, the minimum horizontal principal stress of the mesh, the vertical principal stress of the mesh, the cosine value of the first included angle of the mesh, the cosine value of the second included angle of the mesh, and the cosine value of the third included angle of the mesh; whereby the effective normal stress value of the mesh characterizes the effective normal stress borne by the mesh surface.

[0162] The shear stress value of the mesh is determined based on the maximum horizontal principal stress, the minimum horizontal principal stress, the vertical principal stress, the cosine of the first included angle, the cosine of the second included angle, and the cosine of the third included angle; where the shear stress of the mesh characterizes the shear stress borne by the mesh surface.

[0163] The slip index of the mesh is determined based on the effective normal stress and shear stress of the mesh.

[0164] For example, the normal direction of the mesh surface is determined. In a three-dimensional spatial mesh model, each mesh surface has a unique normal direction, which is perpendicular to the mesh surface. The first angle is the angle between the normal of the mesh surface and the direction of the horizontal maximum principal stress, the second angle is the angle between the normal of the mesh surface and the direction of the horizontal minimum principal stress, and the third angle is the angle between the normal of the mesh surface and the direction perpendicular to the principal stress.

[0165] In one example, Figure 4 This is a schematic diagram of the forces acting on a grid point of a grid in a three-dimensional spatial mesh model of the fault of the target stratigraphic system to be identified, as provided in this application. Figure 4 As shown, line segment The direction of the normal to the mesh surface; line segment The direction of maximum principal stress; line segment The direction of minimum principal stress; line segment The direction is perpendicular to the stress direction. The first included angle β1 of the mesh is ∠DOB; the second included angle β2 of the mesh is ∠DOA; the second included angle β3 of the mesh is ∠DOC. (Line segment) The force in the direction is the finite principal stress of the mesh surface; line segment The force in the direction is the maximum principal stress value minus the pore fluid pressure value, that is, S.i -P; line segment The force in the direction is the minimum principal stress value minus the pore fluid pressure value, that is, S. h -P; line segment The force in the direction is the vertical principal stress value minus the pore fluid pressure value, that is, S. v -P.

[0166] In one example, the cosine of the first included angle of the grid is cosβ1 = sinαsinγ; the cosine of the second included angle of the grid is cosβ2 = sinαcosγ; and the cosine of the third included angle of the grid is cosβ1 =

[0167] sinαsinγ; where α is the tilt angle of the grid points, and γ is the included angle of the grid. The effective normal stress of the grid is σ=(S i -P)cos 2 β1+(S h -P)cos 2 β2+(S v -P)cos 2 β3.

[0168] The shear stress of the mesh is

[0169]

[0170] Among them, S i S represents the maximum horizontal principal stress of the mesh. h S is the minimum principal stress of the horizontal grid. v Let P be the vertical principal stress of the mesh, and P be the pore fluid pressure value of the mesh surface.

[0171] In one example, the slip index of the mesh is H = τ / (μσ + C); where σ is the effective normal stress of the mesh, τ is the shear stress of the mesh, μ is the preset friction coefficient of the mesh surface, and C is a preset constant.

[0172] S309. Determine the slip risk value of the fault of the target stratum to be identified based on the number of grids that have slipped; wherein, the slip risk value represents the risk coefficient of the fault of the target stratum to be identified slipping.

[0173] For example, this step is the same as step S106 above, and will not be repeated here.

[0174] S310. If the slip risk value of the fault of the target stratum to be identified indicates that the risk coefficient of the fault of the target stratum to be identified is a high risk coefficient or a medium risk coefficient, then a warning message is issued to the user; wherein, the warning message indicates that the fault of the target stratum to be identified will slip.

[0175] For example, this step is the same as step S106 above, and will not be repeated here.

[0176] This application provides a method for identifying fault slip risk. It establishes a model of the relationship between pore fluid pressure and formation velocity using methods such as the Bowers method and the Fillippone method. High-precision seismic layer velocities are obtained by constraining seismic inversion with well logging data. Three-dimensional spatial formation pore pressure is calculated using the well logging-established model of formation velocity and pore pressure relationship, along with the seismically inverted layer velocities, allowing for more accurate pore fluid pressure values ​​for each grid cell. Stress field numerical simulation analysis is performed using finite element method software to obtain three-dimensional stress field data. Then, measured stress values ​​from rock cores are obtained using methods such as rock acoustic emission and differential strain methods. The measured and simulated values ​​are corrected by adjusting boundary loads or boundary displacements to obtain an accurate three-dimensional stress field distribution characteristic model. All obtained values ​​are combined to calculate the slip risk value for each grid cell. Based on the slip risk value of each grid cell, the slip risk coefficient of the fault to be developed is obtained, providing a scientific basis for deep shale gas development.

[0177] Figure 5 A schematic diagram of the structure of a fault slip risk identification device provided in this application Figure 1 ,like Figure 5 As shown, the fault slip risk identification device 50 provided in this embodiment includes:

[0178] Module 501 is used to establish a three-dimensional spatial mesh model of the fault of the target stratigraphy to be identified. The three-dimensional spatial mesh model includes the attitude of the mesh points on the fault of the target stratigraphy to be identified. The attitude of the mesh points includes the dip angle of the mesh points on the fault and the orientation of the mesh points. The dip angle of the mesh points represents the dip angle of the mesh points on the mesh, and the orientation of the mesh points represents the orientation of the mesh points on the mesh.

[0179] The acquisition module 502 is used to acquire the pore fluid pressure value of the mesh surface. The pore fluid pressure value of the mesh surface represents the pressure borne by the fluid in the pores on the mesh surface.

[0180] The first determining module 503 is used to determine the principal stress data on the grid surface of the grid according to the preset three-dimensional structural geological model. The preset three-dimensional structural geological model represents the stratigraphic information and fault information of the target strata to be identified, and the principal stress data represents the principal stress borne by the grid surface. The principal stress data includes the maximum horizontal principal stress of the grid, the minimum horizontal principal stress of the grid, and the vertical principal stress of the grid.

[0181] The logging module 504 is used to log the grid to obtain the direction of the horizontal maximum principal stress, the direction of the horizontal minimum principal stress, and the direction of the vertical principal stress.

[0182] The second determining module 505 is used to determine the slip index value of the mesh based on the mesh point orientation, the pore fluid pressure value of the mesh surface, the principal stress data of the mesh, the direction of the maximum horizontal principal stress, the direction of the minimum horizontal principal stress, and the direction of the vertical principal stress of the mesh; wherein, the slip index value characterizes whether the mesh surface will slip.

[0183] The third determining module 506 is used to determine the slip risk value of the fault of the target stratum to be identified based on the number of grids that have slipped; wherein the slip risk value characterizes the risk coefficient of the fault of the target stratum to be identified slipping.

[0184] Figure 6 A schematic diagram of the structure of a fault slip risk identification device provided in this application is shown in Figure 2. Figure 6 As shown, the fault slip risk identification device 60 provided in this embodiment includes:

[0185] Module 601 is used to establish a three-dimensional spatial mesh model of the fault of the target stratigraphy to be identified. The three-dimensional spatial mesh model includes the attitude of the mesh points on the fault of the target stratigraphy to be identified. The attitude of the mesh points includes the dip angle of the mesh points on the fault and the orientation of the mesh points. The dip angle of the mesh points represents the dip angle of the mesh points on the mesh, and the orientation of the mesh points represents the orientation of the mesh points on the mesh.

[0186] The acquisition module 602 is used to acquire the pore fluid pressure value of the mesh surface, which represents the pressure borne by the fluid in the pores on the mesh surface.

[0187] The first determining module 603 is used to determine the principal stress data on the grid surface of the grid according to the preset three-dimensional structural geological model. The preset three-dimensional structural geological model represents the stratigraphic information and fault information of the target strata to be identified, and the principal stress data represents the principal stress borne by the grid surface. The principal stress data includes the maximum horizontal principal stress of the grid, the minimum horizontal principal stress of the grid, and the vertical principal stress of the grid.

[0188] The logging module 604 is used to log the grid to obtain the direction of the horizontal maximum principal stress, the direction of the horizontal minimum principal stress, and the direction of the vertical principal stress.

[0189] The second determining module 605 is used to determine the slip index value of the mesh based on the mesh point orientation, the pore fluid pressure value of the mesh surface, the principal stress data of the mesh, the direction of the maximum horizontal principal stress, the direction of the minimum horizontal principal stress, and the direction of the vertical principal stress of the mesh; wherein, the slip index value characterizes whether the mesh surface will slip.

[0190] The third determining module 606 is used to determine the slip risk value of the fault of the target stratum to be identified based on the number of grids that have slipped; wherein, the slip risk value represents the risk coefficient of the fault of the target stratum to be identified slipping.

[0191] In one example, the orientation of the grid points on the fault is the direction of the line connecting the grid point to the first adjacent grid point. The first adjacent grid point is the next grid point that is on the same grid line as the grid point in the first direction, which is the horizontal direction in the grid division direction in the three-dimensional space grid model.

[0192] The second determining module 605 includes: a first determining submodule 6051, used to determine the angle between the orientation of the grid points and the due north direction, which is the angle value of the grid.

[0193] The second determining submodule 6052 is used to determine the cosine value of the first included angle and the cosine value of the second included angle of the mesh based on the mesh point tilt angle and the included angle value of the mesh, and to determine the cosine value of the third included angle of the mesh based on the mesh point tilt angle; wherein, the first included angle of the mesh is the angle between the normal of the mesh surface and the direction of the maximum horizontal principal stress of the mesh, the second included angle of the mesh is the angle between the normal of the mesh surface and the direction of the minimum horizontal principal stress of the mesh, and the third included angle of the mesh is the angle between the normal of the mesh surface and the direction of the vertical principal stress of the mesh.

[0194] The third determining submodule 6053 is used to determine the effective normal stress value of the mesh based on the pore fluid pressure value of the mesh surface, the maximum horizontal principal stress of the mesh, the minimum horizontal principal stress of the mesh, the vertical principal stress of the mesh, the cosine value of the first included angle of the mesh, the cosine value of the second included angle of the mesh, and the cosine value of the third included angle of the mesh; wherein, the effective normal stress value of the mesh characterizes the effective normal stress borne by the mesh surface of the mesh.

[0195] The fourth determination submodule 6054 is used to determine the shear stress value of the mesh based on the maximum horizontal principal stress, the minimum horizontal principal stress, the vertical principal stress, the cosine value of the first included angle, the cosine value of the second included angle, and the cosine value of the third included angle; wherein, the shear stress of the mesh characterizes the shear stress borne by the mesh surface.

[0196] The fifth determination submodule 6055 is used to determine the slip index value of the mesh based on the effective normal stress value and the shear stress of the mesh.

[0197] In one example, the cosine of the first included angle of the grid is cosβ1 = sinαsinγ; the cosine of the second included angle of the grid is cosβ2 = sinαcosγ; and the cosine of the third included angle of the grid is cosβ1 =

[0198] sinαsinγ; where α is the tilt angle of the grid points and γ is the included angle of the grid.

[0199] In one example, the effective normal stress value of the mesh is σ = (S i -P)cos 2 β1+(S h -P)cos 2 β2+(S v -P)cos 2 β3.

[0200] The shear stress of the mesh is

[0201]

[0202] Among them, S i S represents the maximum horizontal principal stress of the mesh. h S is the minimum principal stress of the horizontal grid. v Let σ be the vertical principal stress of the mesh, and P be the pore fluid pressure at the mesh surface. The mesh slip index is H = τ / (μσ + C); where σ is the effective normal stress of the mesh, τ is the shear stress of the mesh, μ is the preset friction coefficient of the mesh surface, and C is a preset constant.

[0203] In one example, module 601 is created, including:

[0204] The first submodule 6011 is used to establish three-dimensional spatial distribution features based on the three-dimensional seismic data volume; wherein, the three-dimensional seismic data volume includes the fault data of the target strata to be identified; the three-dimensional spatial distribution features characterize the characteristics of the faults of the target strata to be identified in three-dimensional space.

[0205] The submodule 6012 is used to divide the three-dimensional spatial distribution features into a grid to obtain the grid of the fault of the target layer to be identified; and to determine the connection line of the grid; wherein, the connection line of the grid is the line connecting the grid point to the second adjacent grid point adjacent to the grid point, and the second adjacent grid point adjacent to the grid point is the next grid point located on the same grid line as the grid point in the second direction, and the second direction is the perpendicular direction in the grid division direction in the three-dimensional spatial grid model.

[0206] The first generation submodule 6013 is used to determine the angle between the mesh connection line and the physical horizontal line, which is the mesh point tilt angle, in order to obtain a three-dimensional spatial mesh model.

[0207] In one example, module 602 is retrieved, including:

[0208] The acquisition submodule 6021 is used to acquire well logging data, which includes formation velocity; and to establish a relationship model based on the preset measured pore fluid pressure value and the formation velocity in the well logging data. The relationship model is the correspondence between the pore fluid pressure value and the formation velocity.

[0209] The first processing submodule 6022 acquires three-dimensional seismic material data, uses well logging data, performs seismic inversion calculations on the three-dimensional seismic material data, and obtains the seismic layer velocity of the grid surface of the fault of the target stratum to be identified.

[0210] The second processing submodule 6023 is used to input the seismic layer velocity of the grid surface as the formation velocity into the relational model to obtain the pore fluid pressure value of the grid surface.

[0211] In one example, the first determining module 603 includes:

[0212] The second submodule 6035 is used to acquire three-dimensional seismic material data of the target stratigraphic system to be identified. The three-dimensional seismic material data includes stratigraphic interpretation data and fault interpretation data of the target stratigraphic system. Based on the three-dimensional seismic material data, a preset three-dimensional structural geological model is established.

[0213] The detection submodule 6031 is used to detect the mechanical parameters of the target layer to be identified; wherein, the mechanical parameters characterize the mechanical parameters of the target layer.

[0214] The second generation submodule 6032 is used to assign the parameters of the mechanical parameters of the target strata to be identified to the preset three-dimensional structural geological model to obtain a three-dimensional mechanical model; to perform meshing on the three-dimensional mechanical model to obtain a mesh model for finite element analysis; wherein, the mesh model for finite element analysis is a three-dimensional mechanical model with mesh.

[0215] The third processing submodule 6033 is used to perform stress field numerical simulation processing on the mesh model of finite element analysis to obtain the initial principal stress data on the mesh surface.

[0216] The correction submodule 6034 is used to correct the initial principal stress data on the grid surface of the grid based on the pre-acquired core measured stress values ​​of the target stratum to be identified, so as to obtain the principal stress data on the grid surface; wherein, the core measured stress value represents the principal stress borne by the core of the target stratum.

[0217] Logging module 604 is specifically used for:

[0218] Based on pre-acquired imaging logging technology, drilling-induced fractures and wellbore collapse phenomena on faults of the target strata are identified, and the directions of the horizontal maximum principal stress, the horizontal minimum principal stress, and the vertical principal stress of the grid are obtained.

[0219] In one example, the apparatus provided in this application embodiment further includes:

[0220] The warning module 607 is used to issue a warning message to the user if the slip risk value of the fault of the target stratigraphy to be identified indicates that the risk coefficient of the fault of the target stratigraphy to be identified is a high risk coefficient or a medium risk coefficient; wherein the warning message indicates that the fault of the target stratigraphy to be identified will slip.

[0221] This embodiment provides a fault slip risk identification device that can execute the method provided in the above method embodiment. Its implementation principle and technical effect are similar, and will not be described in detail here.

[0222] Figure 7 A schematic diagram of an electronic device structure for identifying fault slip risk provided in this application. Figure 7 As shown, the electronic device 70 provided in this embodiment includes at least one processor 701 and a memory 702. Optionally, the device 70 further includes a communication component 703. The processor 701, memory 702, and communication component 703 are connected via a bus 704.

[0223] In a specific implementation, at least one processor 701 executes computer execution instructions stored in memory 702, causing at least one processor 701 to perform the above-described method.

[0224] The specific implementation process of processor 701 can be found in the above method embodiments, and its implementation principle and technical effect are similar. It will not be repeated here.

[0225] In the above embodiments, it should be understood that the processor can be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), etc. The general-purpose processor can be a microprocessor or any conventional processor. The steps of the method disclosed in this invention can be directly implemented by a hardware processor, or implemented by a combination of hardware and software modules within the processor.

[0226] The memory may include random access memory (RAM) and may also include non-volatile memory (NVM), such as at least one disk storage device.

[0227] The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus, etc. Buses can be categorized as address buses, data buses, control buses, etc. For ease of illustration, the buses shown in the accompanying drawings are not limited to a single bus or a single type of bus.

[0228] This application also provides a computer program product, including a computer program that, when executed by a processor, implements the above-described method.

[0229] This application also provides a computer-readable storage medium storing computer-executable instructions, which, when executed by a processor, implement the above-described method.

[0230] The aforementioned readable storage medium can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic storage, flash memory, magnetic disk, or optical disk. The readable storage medium can be any available medium accessible to a general-purpose or special-purpose computer.

[0231] An exemplary readable storage medium is coupled to a processor, enabling the processor to read information from and write information to the readable storage medium. Of course, the readable storage medium can also be a component of the processor. The processor and the readable storage medium can reside in an Application Specific Integrated Circuit (ASIC). Alternatively, the processor and the readable storage medium can exist as discrete components in the device.

[0232] The division of units is merely a logical functional division; in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be indirect coupling or communication connection through some interfaces, devices, or units, and may be electrical, mechanical, or other forms.

[0233] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0234] In addition, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0235] If a function is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this invention, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of the various embodiments of this invention. 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.

[0236] Those skilled in the art will understand that all or part of the steps of the above-described method embodiments can be implemented by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When executed, the program performs the steps of the above-described method embodiments; and the aforementioned storage medium includes various media capable of storing program code, such as ROM, RAM, magnetic disks, or optical disks.

[0237] Finally, it should be noted that other embodiments of the invention will readily occur to those skilled in the art upon consideration of the specification and practice of the invention disclosed herein. This invention is intended to cover any variations, uses, or adaptations of the invention that follow the general principles of the invention and include common knowledge or customary techniques in the art not disclosed herein, and is not limited to the precise structures described above and shown in the accompanying drawings, and various modifications and changes can be made without departing from its scope. The scope of the invention is limited only by the appended claims.

Claims

1. A method for identifying the slip risk of a fault, characterized in that, include: A three-dimensional spatial mesh model of the fault of the target stratigraphic system to be identified is established. The three-dimensional spatial mesh model includes the attitude of the mesh points on the fault of the target stratigraphic system. The attitude of the mesh points includes the dip angle of the mesh points on the fault and the orientation of the mesh points. The dip angle represents the dip angle of the mesh points on the mesh, and the orientation of the mesh points represents the orientation of the mesh points on the mesh. The pore fluid pressure value of the mesh surface is obtained. The pore fluid pressure value of the mesh surface represents the pressure borne by the fluid in the pores on the mesh surface. Based on a preset three-dimensional structural geological model, the principal stress data on the grid surface are determined; wherein, the preset three-dimensional structural geological model represents the stratigraphic and fault information of the target strata to be identified, and the principal stress data represents the principal stress borne by the grid surface, including the maximum horizontal principal stress, the minimum horizontal principal stress, and the vertical principal stress of the grid. Well logging was performed on the grid to obtain the directions of the horizontal maximum principal stress, the horizontal minimum principal stress, and the vertical principal stress of the grid; The slip index value of the mesh is determined based on the orientation of the mesh points, the pore fluid pressure value of the mesh surface, the principal stress data of the mesh, the direction of the maximum horizontal principal stress, the direction of the minimum horizontal principal stress, and the direction of the vertical principal stress; wherein, the slip index value characterizes whether the mesh surface will slip. The slip risk value of the fault of the target stratigraphic system to be identified is determined based on the number of grids in which slip occurs; wherein the slip risk value characterizes the risk coefficient of slip of the fault of the target stratigraphic system to be identified.

2. The method according to claim 1, characterized in that, The orientation of the grid points on the fault is the direction of the line connecting the grid point to the first adjacent grid point. The first adjacent grid point is the next grid point located on the same grid line as the grid point in the first direction. The first direction is the horizontal direction in the grid division direction in the three-dimensional space grid model. The slip index of the mesh is determined based on the mesh point orientation, the pore fluid pressure value of the mesh surface, the principal stress data of the mesh, the direction of the maximum horizontal principal stress, the direction of the minimum horizontal principal stress, and the direction of the vertical principal stress, including: The angle between the orientation of the grid points and the due north direction is determined as the angle value of the grid. Based on the grid point tilt angle and the included angle of the grid, determine the cosine value of the first included angle and the cosine value of the second included angle of the grid, and based on the grid point tilt angle, determine the cosine value of the third included angle of the grid; wherein, the first included angle of the grid is the angle between the normal to the grid surface and the direction of the maximum horizontal principal stress of the grid, the second included angle of the grid is the angle between the normal to the grid surface and the direction of the minimum horizontal principal stress of the grid, and the third included angle of the grid is the angle between the normal to the grid surface and the direction of the vertical principal stress of the grid; The effective normal stress value of the mesh is determined based on the pore fluid pressure value of the mesh surface, the maximum horizontal principal stress of the mesh, the minimum horizontal principal stress of the mesh, the vertical principal stress of the mesh, the cosine value of the first included angle of the mesh, the cosine value of the second included angle of the mesh, and the cosine value of the third included angle of the mesh; wherein, the effective normal stress value of the mesh characterizes the effective normal stress borne by the mesh surface. The shear stress value of the mesh is determined based on the maximum horizontal principal stress, the minimum horizontal principal stress, the vertical principal stress, the cosine of the first included angle, the cosine of the second included angle, and the cosine of the third included angle; wherein, the shear stress of the mesh characterizes the shear stress borne by the mesh surface. The slip index of the mesh is determined based on the effective normal stress value and the shear stress of the mesh.

3. The method according to claim 2, characterized in that, The cosine of the first included angle of the grid is cosβ1 = sinαsinγ; the cosine of the second included angle of the grid is cosβ2 = sinαcosγ; the cosine of the third included angle of the grid is cosβ1 = sinαsinγ; where α is the tilt angle of the grid point and γ is the included angle value of the grid.

4. The method according to claim 2, characterized in that, The effective normal stress value of the grid is σ = (S i -P)cos 2 β1+(S h -P)cos 2 β2+(S v -P)cos 2 β3; The shear stress of the mesh is Among them, S i S is the maximum horizontal principal stress of the mesh. h S is the minimum principal stress at the horizontal level of the mesh. v P is the vertical principal stress of the mesh, and P is the pore fluid pressure value of the mesh surface; the slip index value of the mesh is H=τ / (μσ+C); where σ is the effective normal stress value of the mesh, τ is the shear stress of the mesh, μ is the preset friction coefficient of the mesh surface, and C is a preset constant.

5. The method according to claim 1, characterized in that, Establish a three-dimensional spatial mesh model of the faults of the target stratigraphic system to be identified, including: Based on the three-dimensional seismic data volume, a three-dimensional spatial distribution feature is established; wherein, the three-dimensional seismic data volume includes the fault data of the target stratum to be identified; the three-dimensional spatial distribution feature characterizes the characteristics of the fault of the target stratum to be identified in three-dimensional space; The three-dimensional spatial distribution features are divided into grids to obtain the grid of the fault of the target layer to be identified; and the connecting lines of the grids are determined; wherein, the connecting lines of the grids are the lines connecting a grid point to the second adjacent grid point adjacent to the grid point, and the second adjacent grid point adjacent to the grid point is the next grid point located on the same grid line as the grid point in a second direction, and the second direction is the perpendicular direction in the grid division direction in the three-dimensional spatial grid model; The angle between the connecting line of the mesh and the physical horizontal line is determined as the tilt angle of the mesh point to obtain the three-dimensional spatial mesh model.

6. The method according to claim 1, characterized in that, Obtaining the pore fluid pressure value at the mesh surface of the mesh includes: Acquire well logging data, which includes formation velocity; and establish a relationship model based on the preset measured pore fluid pressure value and the formation velocity in the well logging data, wherein the relationship model is the correspondence between the pore fluid pressure value and the formation velocity. Acquire three-dimensional seismic material data, and use the well logging data to perform seismic inversion calculations on the three-dimensional seismic material data to obtain the seismic layer velocity of the grid surface of the fault of the target stratum to be identified; The seismic layer velocity of the grid surface is used as the formation velocity and input into the relational model to obtain the pore fluid pressure value of the grid surface.

7. The method according to claim 1, characterized in that, Based on a pre-defined three-dimensional geological model, the principal stress data on the grid surface are determined, including: Acquire three-dimensional seismic material data of the target stratigraphic system to be identified, wherein the three-dimensional seismic material data includes stratigraphic interpretation data and fault interpretation data of the target stratigraphic system; Based on the three-dimensional seismic material data, the preset three-dimensional structural geological model is established; The mechanical parameters of the target layer to be identified are detected; wherein the mechanical parameters characterize the mechanical properties of the target layer. The preset three-dimensional structural geological model is assigned the parameters from the mechanical parameters of the target strata to be identified, thereby obtaining a three-dimensional mechanical model; The three-dimensional mechanical model is meshed to obtain a mesh model for finite element analysis; wherein, the mesh model for finite element analysis is a three-dimensional mechanical model with a mesh. The stress field numerical simulation is performed on the mesh model of the finite element analysis to obtain the initial principal stress data on the mesh surface; Based on the pre-acquired core measured stress values ​​of the target strata to be identified, the initial principal stress data on the grid surface of the grid is corrected to obtain the principal stress data on the grid surface; wherein, the core measured stress value characterizes the principal stress borne by the core of the target strata; Well logging is performed on the grid to obtain the directions of the horizontal maximum principal stress, the horizontal minimum principal stress, and the vertical principal stress of the grid, including: Based on the pre-acquired imaging logging technology, the drilling-induced fracture phenomenon and wellbore collapse phenomenon on the fault of the target stratum to be identified are identified, and the directions of the horizontal maximum principal stress, the horizontal minimum principal stress and the vertical principal stress of the grid are obtained.

8. The method according to any one of claims 1-7, characterized in that, The method further includes: If the slip risk value of the fault of the target stratum to be identified indicates that the risk coefficient of the fault of the target stratum to be identified is a high risk coefficient or a medium risk coefficient, a warning message is issued to the user; wherein, the warning message indicates that the fault of the target stratum to be identified will slip.

9. A device for identifying the slip risk of a fault, characterized in that, The apparatus includes a module for performing the method as described in any one of claims 1-8.

10. An electronic device, characterized in that, include: Memory, processor; The memory stores computer-executed instructions; The processor executes computer execution instructions stored in the memory, causing the processor to perform the method as described in any one of claims 1-8.

11. A computer-readable storage medium / computer program product, characterized in that, The computer-readable storage medium stores computer-executable instructions, which, when executed by a processor, are used to implement the method as described in any one of claims 1-8; and The computer program product includes a computer program that, when executed by a processor, implements the method described in any one of claims 1-8.

Images