Method and apparatus for identifying slip risk of fault, and program product

Abstract

Provided in the embodiments of the present application are a method and apparatus for identifying the slip risk of a fault, and a program product. The method comprises: establishing a three-dimensional spatial grid model for a fault of a target formation series to be identified, and acquiring a grid point attitude of grids on the fault of said target formation series, wherein the grid point attitude comprises the dip angle and strike of each grid point; acquiring pore fluid pressure values of the grid surfaces of the grids; on the basis of a preset three-dimensional structural geological model, determining principal stress data on the grid surfaces of the grids; on the basis of the grid point attitude, the pore fluid pressure values of the grid surfaces of the grids, and the principal stress data of the grids, determining slip index values of the grids, wherein the slip index values indicate whether the grid surfaces of the grids will undergo slip; and on the basis of the number of grids that undergo slip, determining a slip risk coefficient of the fault of said target formation series. The method is used to achieve the effect of accurately determining a slip risk coefficient of a fault of a target formation series.

Description

Method, device and program product for identifying slip risk of fault

[0001] The present application claims priority to the Chinese patent application No. 202411337332.9, filed on September 24, 2024, and entitled "Method, device, electronic device, storage medium and program product for identifying slip risk of fault", the whole content of which is incorporated herein by reference. TECHNICAL FIELD

[0002] The present application relates to the technical field of oil exploitation, and more particularly, to a method, device and program product for identifying slip risk of fault. BACKGROUND

[0003] In the process of deep shale gas exploration and development, it is necessary to determine the risk coefficient of slip of the fault of the target layer system, so as to provide a reference for the design, construction and later maintenance of deep shale gas wells, for example, the target layer system is a shale or other unconventional oil and gas reservoir.

[0004] In the prior art, the data of the fault of the target layer system is analyzed as a whole to determine the risk coefficient of slip of the fault of the target layer system.

[0005] However, the risk coefficient of slip of the fault of the target layer system determined in the above manner is not accurate, and thus there is an urgent need for a scheme that can accurately determine the risk coefficient of slip of the fault of the target layer system. SUMMARY

[0006] The embodiments of the present application provide a method, device and program product for identifying slip risk of fault, which integrates geological, seismic, logging and other multi-source data, combines three-dimensional stress field numerical simulation and three-dimensional space fault plane stress analysis, and establishes a quantitative fault slip risk identification method.

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

[0008] establishing a three-dimensional space grid model of the fault of the target layer system to be identified, wherein the three-dimensional space grid model comprises a grid point occurrence of a grid on the fault of the target layer system to be identified, the grid point occurrence comprises a grid point dip angle and a grid point strike of the grid on the fault, the grid point dip angle represents the dip angle of the grid point on the grid, and the grid point strike represents the strike of the grid point on the grid; and obtaining a pore fluid pressure value of a grid face of the grid, the pore fluid pressure value of the grid face of the grid represents the pressure borne by the fluid in the pore on the grid face.

[0009] According to a preset three-dimensional tectonic geological model, main stress data on a grid surface of the grid is determined; wherein the preset three-dimensional tectonic geological model represents horizon information and fault information of a target layer series to be identified, the main stress data represents a main stress borne by the grid surface, and the main stress data includes a horizontal maximum main stress of the grid, a horizontal minimum main stress of the grid, and a vertical main stress of the grid.

[0010] Logging is performed on the grid to obtain a direction of the horizontal maximum main stress of the grid, a direction of the horizontal minimum main stress of the grid, and a direction of the vertical main stress of the grid.

[0011] According to the grid point occurrence of the grid, the pore fluid pressure value of the grid surface of the grid, the main stress data of the grid, the direction of the horizontal maximum main stress of the grid, the direction of the horizontal minimum main stress of the grid, and the direction of the vertical main stress of the grid, a slip index value of the grid is determined; wherein the slip index value represents whether the grid surface of the grid will slip.

[0012] According to the number of grids that slip, a slip risk value of the fault of the target layer series to be identified is determined; wherein the slip risk value represents a risk coefficient of the slip of the fault of the target layer series to be identified.

[0013] In a possible implementation, the grid point strike of the grid on the fault is a direction of a line connecting the grid point to a first adjacent grid point adjacent to the grid point, and the first adjacent grid point adjacent to the grid point is a next grid point located on a same grid line as the grid point in a horizontal direction of a division direction of the grid in a three-dimensional space grid model.

[0014] According to the grid point occurrence of the grid, the pore fluid pressure value of the grid surface of the grid, the main stress data of the grid, the direction of the horizontal maximum main stress of the grid, the direction of the horizontal minimum main stress of the grid, and the direction of the vertical main stress of the grid, a slip index value of the grid is determined, including:

[0015] An angle between the grid point strike of the grid and the north direction is determined as an angle value of the grid.

[0016] According to the grid point dip angle of the grid and the angle value of the grid, a cosine value of a first angle of the grid and a cosine value of a second angle of the grid are determined, and according to the grid point dip angle of the grid, a cosine value of a third angle of the grid is determined; wherein the first angle of the grid is an angle between a normal of the grid surface and the direction of the horizontal maximum main stress of the grid, the second angle of the grid is an angle between the normal of the grid surface and the direction of the horizontal minimum main stress of the grid, and the third angle of the grid is an angle between the normal of the grid surface and the direction of the vertical main stress of the grid.

[0017] determining an effective normal stress value of the grid according to the pore fluid pressure value of the grid face of the grid, the horizontal maximum principal stress of the grid, the horizontal minimum principal stress of the grid, the vertical principal stress of the grid, the cosine value of the first included angle of the grid, the cosine value of the second included angle of the grid and the cosine value of the third included angle of the grid; wherein the effective normal stress value of the grid represents the effective normal stress borne by the grid face of the grid.

[0018] determining a shear stress value of the grid according to the horizontal maximum principal stress of the grid, the horizontal minimum principal stress of the grid, the vertical principal stress of the grid, the cosine value of the first included angle of the grid, the cosine value of the second included angle of the grid and the cosine value of the third included angle of the grid; wherein the shear stress of the grid represents the shear stress borne by the grid face of the grid.

[0019] determining a slip index value of the grid according to the effective normal stress value of the grid and the shear stress of the grid.

[0020] In a possible implementation, the cosine value of the first included angle of the grid is cosβ1=sinαsinγ. The cosine value of the second included angle of the grid is cosβ2=sinαcosγ; and the cosine value of the third included angle of the grid is cosβ1=sinαsinγ; wherein α is the grid point inclination of the grid, and γ is the included angle value of the grid.

[0021] In a possible implementation, 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.

[0022] The shear stress of the grid is

[0023] wherein S i is the horizontal maximum principal stress of the grid, S h is the horizontal minimum principal stress of the grid, S v is the vertical principal stress of the grid, and P is the pore fluid pressure value of the grid face of the grid. The slip index value of the grid is H=τ / (μσ+C); wherein σ is the effective normal stress value of the grid, τ is the shear stress of the grid, μ is the preset friction coefficient of the grid face of the grid, and C is a preset constant.

[0024] In a possible implementation, a three-dimensional spatial grid model of faults of a target layer system to be identified is established, including:

[0025] According to a three-dimensional seismic data body, a three-dimensional spatial distribution feature is established; wherein the three-dimensional seismic data body includes data of faults of the target layer system to be identified; and the three-dimensional spatial distribution feature represents a feature of the faults of the target layer system to be identified in three-dimensional space.

[0026] The three-dimensional spatial distribution feature is grid-divided to obtain a grid of the faults of the target layer system to be identified; and a connection line of the grid is determined; wherein the connection line of the grid is a line connecting a grid point to a second adjacent grid point adjacent to the grid point, and the second adjacent grid point adjacent to the grid point is a next grid point located on a same grid line as the grid point in a second direction, and the second direction is a vertical direction among division directions of the grid in the three-dimensional spatial grid model.

[0027] An angle between the connection line of the grid and a physical horizontal line is determined as a dip angle of the grid point of the grid, so as to obtain the three-dimensional spatial grid model.

[0028] In a possible implementation, a pore fluid pressure value of a grid face of the grid is obtained, including:

[0029] Logging data is obtained, and the logging data includes a formation velocity; and a relationship model is established according to a preset measured pore fluid pressure value and the formation velocity in the logging data, and the relationship model is a corresponding relationship between the pore fluid pressure value and the formation velocity. Three-dimensional seismic material data is obtained, and the logging data is used to perform seismic inversion calculation and processing on the three-dimensional seismic material data, so as to obtain a seismic interval velocity of a grid face of a grid of the faults of the target layer system to be identified.

[0030] The seismic interval velocity of the grid face of the grid is input as the formation velocity into the relationship model, so as to obtain the pore fluid pressure value of the grid face of the grid.

[0031] In a possible implementation, according to a preset three-dimensional structural geology model, main stress data on the grid face of the grid is determined, including:

[0032] The three-dimensional seismic material data of the target layer system to be identified is obtained, wherein the three-dimensional seismic material data includes horizon interpretation data and fault interpretation data of the target layer system.

[0033] The preset three-dimensional structural geology model is established according to the three-dimensional seismic material data.

[0034] Mechanical parameters of the target layer system to be identified are detected; wherein the mechanical parameters represent parameters that the target layer system has in mechanics.

[0035] Assign a parameter in mechanical parameters of the preset three-dimensional structural geological model to a target layer system to be identified to obtain a three-dimensional mechanical model.

[0036] Grid division is performed on the three-dimensional mechanical model to obtain a grid model for finite element analysis; wherein the grid model for finite element analysis is a three-dimensional mechanical model with grids.

[0037] Stress field numerical simulation processing is performed on the grid model for finite element analysis to obtain initial principal stress data on a grid surface of the grid.

[0038] Based on a core measured stress value of the target layer system to be identified obtained in advance, the initial principal stress data on the grid surface of the grid is corrected to obtain principal stress data on the grid surface of the grid; wherein the core measured stress value represents the principal stress borne by the core of the target layer system.

[0039] Logging is performed on 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, including:

[0040] Based on imaging logging technology obtained in advance, drilling-induced fracture phenomenon and wellbore collapse phenomenon on the fault of the target layer system to be identified are identified 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.

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

[0042] If the slip risk value of the fault of the target layer system to be identified represents a high risk coefficient or a medium risk coefficient of the risk of the fault of the target layer system to be identified to slip, a warning information is sent to a user; wherein the warning information represents that the fault of the target layer system to be identified will slip.

[0043] In a second aspect, an embodiment of the present application provides a device for identifying the slip risk of a fault, including:

[0044] A establishing module is configured to establish a three-dimensional space grid model of a fault of a target layer system to be identified, wherein the three-dimensional space grid model includes grid point occurrences of grids on the fault of the target layer system to be identified, the grid point occurrences include grid point dips and grid point strikes of the grids on the fault, the grid point dip represents the dip of the grid point on the grid, and the grid point strike represents the strike of the grid point on the grid.

[0045] The acquisition module is configured to acquire a pore fluid pressure value of a grid face of the grid, the pore fluid pressure value of the grid face of the grid representing a pressure borne by fluid in a pore on the grid face.

[0046] The first determination module is configured to determine main stress data on the grid face of the grid according to a preset three-dimensional tectonic geological model, wherein the preset three-dimensional tectonic geological model represents stratigraphic information and fault information of a target layer system to be identified, the main stress data represents a main stress borne by the grid face, and the main stress data includes a horizontal maximum main stress of the grid, a horizontal minimum main stress of the grid, and a vertical main stress of the grid.

[0047] The logging module is configured to perform logging on the grid to obtain a direction of the horizontal maximum main stress of the grid, a direction of the horizontal minimum main stress of the grid, and a direction of the vertical main stress of the grid.

[0048] The second determination module is configured to determine a slip index value of the grid according to the occurrence of the grid point of the grid, the pore fluid pressure value of the grid face of the grid, the main stress data of the grid, the direction of the horizontal maximum main stress of the grid, the direction of the horizontal minimum main stress of the grid, and the direction of the vertical main stress of the grid, wherein the slip index value represents whether the grid face of the grid will slip.

[0049] The third determination module is configured to determine a slip risk value of a fault of the target layer system to be identified according to the number of grids that slip, wherein the slip risk value represents a risk coefficient of the fault of the target layer system to be identified slipping.

[0050] In a possible implementation, the grid point strike of the grid on the fault is a direction of a line connecting the grid point to a first adjacent grid point adjacent to the grid point, the first adjacent grid point adjacent to the grid point is a next grid point on a same grid line as the grid point in a horizontal direction of a division direction of the grid in a three-dimensional space grid model, and the first direction is the horizontal direction.

[0051] The second determination module includes:

[0052] The first determination submodule is configured to determine an included angle between the grid point strike of the grid and a due north direction as the included angle value of the grid.

[0053] The second determining submodule is configured to determine a cosine value of a first included angle of the grid and a cosine value of a second included angle of the grid according to the grid point dip angle of the grid and the included angle value of the grid, and determine a cosine value of a third included angle of the grid according to the grid point dip angle of the grid; wherein the first included angle of the grid is an included angle between a normal of a grid face and a direction of a horizontal maximum principal stress of the grid, the second included angle of the grid is an included angle between the normal of the grid face and a direction of a horizontal minimum principal stress of the grid, and the third included angle of the grid is an included angle between the normal of the grid face and a direction of a vertical principal stress of the grid.

[0054] The third determining submodule is configured to determine an effective normal stress value of the grid according to the pore fluid pressure value of the grid face of the grid, the horizontal maximum principal stress of the grid, the horizontal minimum principal stress of the grid, the vertical principal stress of the grid, the cosine value of the first included angle of the grid, the cosine value of the second included angle of the grid, and the cosine value of the third included angle of the grid; wherein the effective normal stress value of the grid represents an effective normal stress borne by the grid face of the grid.

[0055] The fourth determining submodule is configured to determine a shear stress value of the grid according to the horizontal maximum principal stress of the grid, the horizontal minimum principal stress of the grid, the vertical principal stress of the grid, the cosine value of the first included angle of the grid, the cosine value of the second included angle of the grid, and the cosine value of the third included angle of the grid; wherein the shear stress of the grid represents a shear stress borne by the grid face of the grid.

[0056] The fifth determining submodule is configured to determine a slip index value of the grid according to the effective normal stress value of the grid and the shear stress of the grid.

[0057] In a possible implementation, the cosine value of the first included angle of the grid is cosβ1=sinαsinγ; the cosine value of the second included angle of the grid is cosβ2=sinαcosγ; and the cosine value of the third included angle of the grid is cosβ1=sinαsinγ; wherein α is the grid point dip angle of the grid, and γ is the included angle value of the grid.

[0058] In a possible implementation, 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.

[0059] The shear stress of the grid is

[0060] wherein S i is a horizontal maximum principal stress of the grid, S h is a horizontal minimum principal stress of the grid, S v is a vertical principal stress of the grid, P is a pore fluid pressure value of a grid face of the grid. A slip index value of the grid is H = τ / (μσ+C); wherein σ is an effective normal stress value of the grid, τ is a shear stress of the grid, μ is a preset friction coefficient of a grid face of the grid, and C is a preset constant.

[0061] In a possible implementation, the establishing module comprises:

[0062] A first establishing submodule is configured to establish a three-dimensional spatial distribution feature according to a three-dimensional seismic data volume; wherein the three-dimensional seismic data volume comprises data of a fault of a target layer system to be identified; and the three-dimensional spatial distribution feature represents a feature of the fault of the target layer system to be identified in a three-dimensional space.

[0063] A dividing submodule is configured to divide the three-dimensional spatial distribution feature into grids to obtain a grid of the fault of the target layer system to be identified; and determine a connecting line of the grid; wherein the connecting line of the grid is a line connecting a grid point to a second adjacent grid point adjacent to the grid point, and the second adjacent grid point adjacent to the grid point is a next grid point located on a same grid line as the grid point in a second direction, and the second direction is a vertical direction in a division direction of a grid in a three-dimensional spatial grid model.

[0064] A first generating submodule is configured to determine an included angle between the connecting line of the grid and a physical horizontal line as a grid point dip angle of the grid to obtain the three-dimensional spatial grid model.

[0065] In a possible implementation, the obtaining module comprises:

[0066] An obtaining submodule is configured to obtain logging data, wherein the logging data comprises a formation velocity; and establish a relationship model according to a preset measured pore fluid pressure value and the formation velocity in the logging data, wherein the relationship model represents a corresponding relationship between a pore fluid pressure value and a formation velocity.

[0067] A first processing submodule is configured to obtain three-dimensional seismic material data, use the logging data to perform seismic inversion calculation processing on the three-dimensional seismic material data, and obtain a seismic interval velocity of a grid face of a grid of a fault of a target layer system to be identified.

[0068] A second processing submodule is configured to input the seismic interval velocity of the grid face of the grid of the fault of the target layer system to be identified as a formation velocity to the relationship model to obtain a pore fluid pressure value of the grid face of the grid.

[0069] In a possible implementation, the first determining module comprises:

[0070] The second establishing sub-module is configured to acquire three-dimensional seismic material data of the target layer series to be identified, wherein the three-dimensional seismic material data comprises horizon interpretation data and fault interpretation data of the target layer series; and establish the preset three-dimensional structural geological model according to the three-dimensional seismic material data.

[0071] The detecting sub-module is configured to detect a mechanical parameter of the target layer series to be identified, wherein the mechanical parameter represents a parameter that the target layer series has in mechanics.

[0072] The second generating sub-module is configured to assign a parameter in the mechanical parameter of the target layer series to be identified to the preset three-dimensional structural geological model to obtain a three-dimensional mechanical model; and perform grid division on the three-dimensional mechanical model to obtain a grid model for finite element analysis, wherein the grid model for finite element analysis is a three-dimensional mechanical model with grids.

[0073] The third processing sub-module is configured to perform stress field numerical simulation processing on the grid model for finite element analysis to obtain initial principal stress data on a grid surface of the grid.

[0074] The correction sub-module is configured to correct the initial principal stress data on the grid surface of the grid based on a core measured stress value of the target layer series to be identified, to obtain principal stress data on the grid surface of the grid, wherein the core measured stress value represents a principal stress that a core of the target layer series bears.

[0075] The logging module is specifically configured to identify, based on pre-acquired imaging logging technology, a drilling-induced fracture phenomenon and a borehole wall collapse phenomenon on a fault of the target layer series to be identified, to obtain a direction of a horizontal maximum principal stress, a direction of a horizontal minimum principal stress, and a direction of a vertical principal stress of the grid.

[0076] In a possible implementation, the device further comprises:

[0077] The warning module is configured to issue a warning information to a user if a slip risk value of the fault of the target layer series to be identified represents a high risk coefficient or a medium risk coefficient of a risk of the fault of the target layer series to be identified to slip, wherein the warning information represents that the fault of the target layer series to be identified will slip.

[0078] In a third aspect, an embodiment of the present application provides an electronic device, comprising: a memory, a processor;

[0079] The memory stores computer execution instructions.

[0080] The processor executes computer-executed instructions stored in the memory, so that the processor executes the method provided by the first aspect and any implementation of the first aspect.

[0081] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium, wherein computer-executed instructions are stored in the computer-readable storage medium, and the computer-executed instructions are executed by a processor to implement the method in the first aspect and / or various possible implementation manners of the first aspect.

[0082] In a fifth aspect, an embodiment of the present application provides a computer program product, which comprises a computer program, and the computer program is executed by a processor to implement the method in the first aspect and any implementation of the first aspect.

[0083] The method, device and program product for identifying the slip risk of a fault provided by the embodiments of the present application can establish a three-dimensional space grid model of a fault of a target layer system to be identified, obtain an included angle and a dip angle of a grid point, obtain a pore fluid pressure value of a grid surface of the grid, determine principal stress data and a direction on the grid surface of the grid, calculate a slip index of each grid, and thus obtain a slip risk coefficient of the fault of the target layer system to be identified. When the slip risk coefficient of the fault of the target layer system to be identified is a high risk coefficient or a medium risk coefficient, a warning information is sent to a user. The method can accurately identify the slip risk of the fault of the target layer system, and provides a scientific basis for drilling operation and geological disaster prevention. BRIEF DESCRIPTION OF DRAWINGS

[0084] FIG. 1 is a flowchart of a method for identifying the slip risk of a fault according to an embodiment of the present application;

[0085] FIG. 2a is a schematic diagram of a three-dimensional space grid structure of a fault of a target layer system to be identified according to an embodiment of the present application;

[0086] FIG. 2b is a left view of a three-dimensional space grid structure of a fault of a target layer system to be identified according to an embodiment of the present application;

[0087] FIG. 2c is a front view of a three-dimensional space grid structure of a fault of a target layer system to be identified according to an embodiment of the present application;

[0088] FIG. 3 is a flowchart of a method for identifying the slip risk of a fault according to another embodiment of the present application;

[0089] FIG. 4 is a schematic diagram of forces borne by a grid point of a grid of a three-dimensional space grid model of a fault of a target layer system to be identified according to an embodiment of the present application;

[0090] FIG. 5 is a schematic diagram of a structure of a device for identifying the slip risk of a fault according to an embodiment of the present application;

[0091] Fig. 6 is a structural schematic diagram of a device for identifying the slip risk of a fault according to an embodiment of the present application;

[0092] Fig. 7 is a structural schematic diagram of an electronic device for identifying the slip risk of a fault according to an embodiment of the present application. DETAILED DESCRIPTION

[0093] The exemplary embodiments will be described in detail herein with reference to the attached drawings. In the following description, the same numbers are used to indicate the same or similar elements, unless otherwise represented. The embodiments described in the following exemplary embodiments do not represent all the embodiments consistent with the present application. Rather, they are merely examples of apparatuses and methods consistent with some aspects of the present application as detailed in the appended claims.

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

[0095] At present, the possibility of fault slip is preliminarily judged by studying the influence of factors such as geostress field, fracture occurrence, pore pressure and the like on the fault slip mechanism. However, most of these studies are limited to qualitative analysis, lacking quantitative identification means, especially in the identification of fault slip risk in three-dimensional space.

[0096] Therefore, the method, device, equipment, storage medium and program product for identifying the slip risk of a fault provided by the present application can accurately identify the fault slip risk, provide a scientific basis for the design, construction and post-maintenance of deep shale gas wells, and ensure the efficient and safe development of shale gas.

[0097] The technical solutions of the present application and how the technical solutions of the present application solve the above technical problems will be described in detail below with specific embodiments. The following specific embodiments can be combined with each other, and the same or similar concepts or processes can not be described again in some embodiments. The embodiments of the present application will be described below with reference to the drawings.

[0098] Fig. 1 is a flowchart of a method for identifying the slip risk of a fault according to an embodiment of the present application, as shown in Fig. 1, the method comprises:

[0099] S101, establish a three-dimensional space grid model of the fault of the target layer system to be identified, wherein the three-dimensional space grid model comprises grid point occurrences of the grid on the fault of the target layer system to be identified, and the grid point occurrences comprise a grid point dip angle and a grid point strike of the grid on the fault, the grid point dip angle represents a dip angle of the grid point on the grid, and the grid point strike represents a strike of the grid point on the grid.

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

[0101] Illustratively, to establish a three-dimensional space grid model of the fault of the shale layer to be identified, geological data about the shale layer to be identified and its fault are collected, including geological exploration reports, seismic interpretation data, drilling data, etc. These data will serve as the basis for establishing the three-dimensional space grid model. Subsequently, the data are preprocessed, including denoising, interpolation, format conversion, etc., to ensure the accuracy and consistency of the data. In the preprocessed data, the specific location, shape and strike of the fault in the shale layer are determined by using geological interpretation software or manual identification methods. According to the fault identification results and the spatial distribution characteristics of the shale layer, the grid is divided in the three-dimensional space. The division of the grid should take into account the shape, strike of the fault and the thickness, dip angle of the shale layer, etc., to ensure that the grid can accurately reflect the spatial structure of the shale layer and its fault. The three-dimensional space grid model comprises grid point occurrences of the grid on the fault of the shale layer to be identified, wherein the grid point occurrences mainly comprise a grid point dip angle and a grid point strike of the grid on the fault. The grid point dip angle refers to the inclination angle of each grid point on the grid relative to the physical horizontal line, which represents the dip angle of the rock layer at the location of the grid point, and the grid point dip angle is assigned according to the geological data and interpretation results; the grid point strike refers to the extension direction of the fault plane on the surface at a certain grid point on the fault.

[0102] S102, obtain a pore fluid pressure value of a grid face of the grid, wherein the pore fluid pressure value of the grid face of the grid represents the pressure borne by the fluid in the pore on the grid face.

[0103] Illustratively, three-dimensional seismic data, acoustic logging data, measured pore fluid pressure data, etc. are collected, and the formation velocity at different depths is obtained through the acoustic logging data, which reflects the variation law of the formation velocity with depth. The original data are processed, including removing outliers, filtering, etc., to ensure the accuracy and reliability of the data.

[0104] The relationship between the formation porosity and the depth, the acoustic travel time and other parameters can be established by using the acoustic logging data, and then the formation pore fluid pressure can be predicted. The acoustic travel time is the time required for the acoustic wave to pass through a certain rock formation, and its size is closely related to the porosity of the rock formation. The greater the porosity, the slower the speed of acoustic wave propagation, and the greater the acoustic travel time. Under normal compaction conditions, as the burial depth of the formation increases, the formation porosity gradually decreases, and the acoustic travel time also decreases accordingly. However, in the abnormal high pressure formation, due to the increase of the pore fluid pressure, the formation porosity increases relatively, and the acoustic travel time deviates from the normal trend line. The formation velocity generally refers to the speed of acoustic wave propagation in the formation, and its reciprocal is the acoustic travel time. The increase of pore pressure will change the effective stress between rock particles, and then affect the elastic modulus and density of the rock, ultimately leading to the change of the formation velocity. Generally speaking, the increase of pore pressure will make the formation velocity decrease, because the increase of pore pressure will reduce the effective stress between rock particles, making the rock become more "soft" and the speed of acoustic wave propagation slower.

[0105] Fine geological horizon interpretation is performed on the three-dimensional seismic data volume to determine the spatial position and shape of each horizon. Logging data is used as a constraint condition to perform inversion processing on the seismic data to obtain high-precision seismic interval velocity data. The seismic interval velocity data reflects the physical properties of the rock in the formation and is closely related to the pore fluid pressure.

[0106] The seismic interval velocity data obtained by seismic inversion is substituted into the previously established pore fluid pressure and formation velocity relationship model to calculate the pore fluid pressure distribution in three-dimensional space. According to the results of fault grid division, the calculated three-dimensional pore fluid pressure data is distributed to each grid surface to obtain the pore fluid pressure value on each grid surface.

[0107] S103, according to the preset three-dimensional structural geological model, the principal stress data on the grid surface of the grid is determined. The preset three-dimensional structural geological model represents the horizon information and fault information of the target layer system to be identified, and the principal stress data represents the principal stress borne by the grid surface. The principal stress data includes the horizontal maximum principal stress of the grid, the horizontal minimum principal stress of the grid, and the vertical principal stress of the grid.

[0108] Illustratively, three-dimensional seismic data, acoustic logging data, geological exploration reports, etc. are obtained to construct a three-dimensional geological model. The three-dimensional seismic data is finely interpreted to identify the spatial distribution characteristics of the faults of the shale layer to be identified, including the strike, dip angle and changes in three-dimensional space. The seismic interpretation results, acoustic logging data, etc. are integrated into the three-dimensional geological model using geological modeling software (such as Petrel, Gocad, etc.) to construct a fine three-dimensional structural geological model. The model should accurately reflect the geometric shape and positional relationship of the target layer system and faults.

[0109] Shale cores of the target layer are collected, and rock mechanics experiments are performed to obtain the Young's modulus and Poisson's ratio of the shale. According to the difference in rock mechanics properties between the fault and the surrounding shale, the mechanical parameters of the fault are assigned. Generally, the Young's modulus of the fault is 40% to 70% of that of the shale, and the Poisson's ratio is 0.02 to 0.1 greater than that of the shale, and can be further adjusted according to the fracture level. On the basis of the three-dimensional geological model, the fault is divided into grids. The size of the grid should be determined according to the calculation accuracy and the demand for computing resources to ensure that the stress change of the fault plane can be accurately reflected.

[0110] The three-dimensional geological model with assigned mechanical parameters is converted into a three-dimensional mechanical model for stress field numerical simulation. The stress boundary conditions of the research area are determined through single-well experimental test data, including the directions and sizes of the horizontal maximum and minimum principal stresses, and the vertical stress. The mechanical model is divided into grids to establish a grid model for finite element analysis. The finite element method is used for stress field simulation and prediction, and the stress field distribution in the three-dimensional space is calculated by solving the mechanical equation. The stress field simulation results are corrected using the measured stress values of the cores to ensure the accuracy of the simulation results. The correction process may need to be iterated several times until the error between the simulation results and the measured values is controlled within 10%.

[0111] The stress data on each grid plane are extracted from the three-dimensional stress field numerical simulation results. According to the stress tensor theory, the horizontal maximum principal stress, the horizontal minimum principal stress, and the vertical principal stress on each grid plane are calculated. These principal stress data represent the main stress state of the grid plane.

[0112] S104, 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.

[0113] In one example, using pre-acquired imaging logging data, drilling-induced fracture and borehole collapse phenomena on faults in shale formations are identified by image analysis techniques. These phenomena are often related to the orientation of the stress field, especially the orientation of drilling-induced fractures and the pattern of borehole collapse, which can indirectly reflect the orientation of the horizontal maximum principal stress. Drilling-induced fractures are fractures formed in the borehole wall due to changes in borehole pressure or the invasion of drilling fluid during drilling. The direction of these fractures is often closely related to the stress state around the borehole, especially the direction of the maximum principal stress. Borehole collapse refers to the collapse of the borehole wall due to stress imbalance during drilling. The direction and pattern of borehole collapse can also reflect the stress state around the borehole, especially the directions of the horizontal maximum and minimum 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 horizontal maximum and minimum principal stresses, and the starting point of the vertical principal stress direction is at the intersection of the orthogonal plane.

[0114] S105, determining a slip index value of the grid according to the grid point occurrence of the grid, the pore fluid pressure value of the grid face of the grid, the principal stress data of the grid, the direction of the horizontal maximum principal stress of the grid, the direction of the horizontal minimum principal stress, and the direction of the vertical principal stress, wherein the slip index value represents whether the grid face of the grid will slip.

[0115] For example, FIG. 2a is a schematic diagram of a three-dimensional spatial grid structure of a fault of a target layer system to be identified; FIG. 2b is a left view of the schematic diagram of the three-dimensional spatial grid structure of the fault of the target layer system to be identified; and FIG. 2c is a front view of the schematic diagram of the three-dimensional spatial grid structure of the fault of the target layer system to be identified.

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

[0117] In one example, as shown in FIG. 2b, in the left view of the three-dimensional spatial grid structure diagram of the fault of the target layer series to be identified, point A is the grid point of the study grid, the strike of the A grid point is the direction of the connecting line from the A grid point to the next grid point (B grid point) of the A grid point in the first direction, and the included angle of the A grid point is the included angle between the strike of the A grid point and the north direction, which is denoted by θ.

[0118] In one example, as shown in FIG. 2c, in the front view of the three-dimensional spatial grid structure diagram of the fault of the target layer series to be identified, point I is the grid point of the study grid, the dip angle of the I grid point is the included angle between the connecting line of the I grid point and the next grid point (H grid point) of the I grid point in the second direction and the physical horizontal line, which is denoted by α, wherein the physical horizontal line refers to the line perpendicular to the direction of the gravity of the earth, i.e. the horizontal line.

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

[0120] The cosine value of the first included angle of the grid is cosβ1=sinαsinγ; the cosine value of the second included angle of the grid is cosβ2=sinαcosγ; the cosine value of the third included angle of the grid is cosβ1=sinαsinγ; wherein α is the dip angle of the grid point of the grid, and γ is the included angle value of the grid; 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; and the shear stress of the grid is The slip index value of the grid is H=τ / (μσ+C).

[0121] In step S106, according to the number of grids that slip, the slip risk value of the fault of the target layer series to be identified is determined; wherein the slip risk value represents the risk coefficient of the slip of the fault of the target layer series to be identified.

[0122] Exemplarily, the slip risk value H is calculated for each grid on the fault of the shale layer to be identified, the H value reflecting the possibility of the grid surface slipping, and if H>0, it indicates that the grid surface has a slip risk. The number of slip risk values H of all grids on the fault of the shale layer to be identified is counted, and 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 a 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 the fault of the shale layer to be identified slipping is a high-risk coefficient or a medium-risk coefficient, a warning information is sent to the user; wherein the warning information indicates that the fault of the shale layer to be identified will slip.

[0123] The method for identifying the slip risk of a fault provided by the embodiments of the present application can accurately identify the slip risk of the fault of the target layer system, and provide a scientific basis for drilling operations and geological disaster prevention.

[0124] FIG. 3 is a flowchart of a method for identifying the slip risk of a fault provided by the present application, as shown in FIG. 3, the method comprises:

[0125] S301, according to the three-dimensional seismic data body, a three-dimensional spatial distribution characteristic is established; wherein the three-dimensional seismic data body comprises data of the fault of the target layer system to be identified; the three-dimensional spatial distribution characteristic represents the characteristics of the fault of the target layer system to be identified in the three-dimensional space.

[0126] Exemplarily, a three-dimensional seismic data volume is collected. These data volumes are obtained through seismic exploration techniques, containing reflection wave information of underground rock structures (including faults of shale layers). The collected three-dimensional seismic data volume is preprocessed, including steps such as denoising, filtering, amplitude restoration, etc., to improve data quality, reduce interference signals, and make the features of underground structures clearer. After data preprocessing, seismic data interpretation is performed. By comparing the features of seismic waveforms (such as amplitude, frequency, phase, etc.), combined with geological background knowledge, the location and shape of the to-be-identified shale layer and the faults therein are identified. Fault identification is based on the features such as breakage and distortion of seismic reflection events. Based on the results of seismic data interpretation, using geological modeling software or three-dimensional visualization technology, the faults of the to-be-identified shale layer are reconstructed and displayed in three-dimensional space. By constructing a three-dimensional grid model, the spatial form, positional relationship, size, etc. of the shale layer and its faults are described, forming three-dimensional spatial distribution characteristics.

[0127] S302, grid division is performed on the three-dimensional spatial distribution characteristics to obtain a grid of the faults of the to-be-identified target layer system; and a connection line of the grid is determined; wherein the connection line of the grid is a line connecting a grid point to a second adjacent grid point adjacent to the grid point, and the second adjacent grid point adjacent to the grid point is a next grid point located on the same grid line as the grid point in a second direction, and the second direction is a vertical direction among the division directions of the grid in the three-dimensional spatial grid model.

[0128] Exemplarily, grid division is performed on the already constructed three-dimensional spatial distribution characteristics (i.e. the three-dimensional form of the faults of the to-be-identified shale layer). The purpose of this step is to discretize the continuous three-dimensional space into a series of regular or irregular grid cells for numerical analysis or visualization display. After grid division is completed, each grid cell is defined by a grid point, and a grid point has 6 directly adjacent grid points. The connection line of the grid is the connection line of the grid point and the next grid point in the horizontal direction among the division directions of the grid, wherein the next grid point is the grid point on the left of the grid point in the horizontal direction. Among them, the horizontal direction is the second direction among the grid division directions in the three-dimensional spatial grid model.

[0129] In one example, as shown in FIG. 2a, the connection line of the H grid point is the line connecting the H grid point and the C grid point; the connection line of the C grid point is the line connecting the C grid point and the G grid point; and the connection line of the G grid point is the line connecting the G grid point and the F grid point.

[0130] S303, an included angle between the connection line of the grid and the physical horizontal line is determined as the dip angle of the grid point of the grid, to obtain a three-dimensional spatial grid model.

[0131] Exemplarily, the inclination of the grid point is the angle between the connecting line of the grid point and the physical horizontal line, wherein the physical horizontal line refers to a horizontal line perpendicular to the direction of the earth's gravity.

[0132] In one example, as shown in FIG. 2c, the inclination of the grid point I is the angle between the connecting line IH of the grid point I and the physical horizontal line of the grid point I, that is, α1; and the inclination of the grid point G is the angle between the connecting line FG of the grid point G and the physical horizontal line of the grid point G, that is, α2.

[0133] S304, obtaining logging data, wherein the logging data comprises formation velocity; and establishing a relationship model according to the measured pore fluid pressure value and the formation velocity in the logging data, wherein the relationship model is a corresponding relationship between the pore fluid pressure value and the formation velocity.

[0134] Exemplarily, the logging data containing various physical properties of the underground rock formation is obtained; wherein the logging data mainly includes formation velocity; the formation velocity is the propagation speed of acoustic waves in the formation. The measured pore fluid pressure value is obtained, wherein the measured pore fluid pressure value refers to the pressure of the fluid (usually water, oil or gas) in the underground rock pores; the obtained data is sorted and cleaned to ensure the quality and integrity of the data, and then the data is corrected as necessary to eliminate instrument errors or environmental factors. The relationship model between the measured pore fluid pressure value and the formation velocity is established, and the specific modeling methods include Bowers method, Fillippone method, etc. Whether the established model is accurate and reliable is verified by cross verification, and the model is adjusted and optimized according to the verification result.

[0135] S305, obtaining three-dimensional seismic material data, using the logging data to perform seismic inversion calculation and processing on the three-dimensional seismic material data to obtain the seismic interval velocity of the grid surface of the grid of the fault of the target layer system to be identified; inputting the seismic interval velocity of the grid surface of the grid as the formation velocity into the relationship model to obtain the pore fluid pressure value of the grid surface of the grid.

[0136] Exemplarily, the three-dimensional seismic material data is interpreted in a fine geological horizon, where the fine refers to that the three-dimensional seismic material data is also divided into grids in the form of the grid division of the three-dimensional spatial grid model in the S101 step. The formation velocity information in the logging data is combined with the three-dimensional seismic material data to provide an additional constraint condition for seismic inversion. After the seismic inversion calculation, the seismic interval velocity of the grid surface of the fault grid is obtained. The seismic inversion can select a pre-stack inversion method such as seismic travel time inversion and tomographic inversion. The seismic interval velocity of the grid surface of the grid is input as the formation velocity into the previously established corresponding relationship model between the pore fluid pressure value and the formation velocity, and the pore fluid pressure value of the grid surface of the grid is calculated by using the relationship model, and finally the pore fluid pressure value of the grid surface of the fault grid of the shale layer is obtained.

[0137] S306, determining the principal stress data on the grid surface of the grid according to the preset three-dimensional structural geological model; wherein, the preset three-dimensional structural geological model represents the horizon information and fault information of the target layer system to be identified, the principal stress data represents the principal stress borne by the grid surface, and the principal stress data includes the horizontal maximum principal stress of the grid, the horizontal minimum principal stress of the grid and the vertical principal stress of the grid.

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

[0139] Obtaining three-dimensional seismic material data of the target layer system to be identified, wherein the three-dimensional seismic material data includes horizon interpretation data and fault interpretation data of the target layer system; and establishing a preset three-dimensional structural geological model according to the three-dimensional seismic material data.

[0140] Detecting the mechanical parameters of the target layer system to be identified; wherein the mechanical parameters represent the parameters possessed by the target layer system in mechanics.

[0141] Assigning the parameters in the mechanical parameters of the target layer system to be identified to the preset three-dimensional structural geological model to obtain a three-dimensional mechanical model.

[0142] Grid dividing the three-dimensional mechanical model to obtain a grid model for finite element analysis; wherein the grid model for finite element analysis is a three-dimensional mechanical model with grids.

[0143] Stress field numerical simulation processing is performed on the grid model for finite element analysis to obtain the initial principal stress data on the grid surface of the grid.

[0144] Based on the previously obtained core measured stress value of the target layer system 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 of the grid; wherein the core measured stress value represents the principal stress borne by the core of the target layer system.

[0145] Exemplarily, three-dimensional seismic material data of the research area is acquired and organized, and the three-dimensional seismic material data is imported into a professional geology interpretation software such as Petrel. In the software, the top and bottom interfaces of the shale layer to be identified are identified and tracked by using the amplitude, frequency and other attributes of the seismic data, and the horizon interpretation is completed. In the software, the faults are identified and interpreted by analyzing the discontinuity and abnormal characteristics in the seismic data, including the strike, dip and displacement of the faults, and the fault interpretation is completed. Based on the results of the horizon interpretation and the fault interpretation, a three-dimensional tectonic geological model is constructed.

[0146] Shale cores are collected from the shale layer to be identified, and are processed into 50mm*25mm plunger samples according to the standard. The GCTS rock mechanics test system is used to conduct triaxial loading experiments on the plunger samples to measure the Young's modulus and Poisson's ratio. The mechanical parameters obtained by the experiment are assigned to the shale part of the geological model. Based on the difference in mechanical properties between the faults and the surrounding shale, relatively low Young's modulus and relatively high Poisson's ratio are set for the faults, considering the influence of fault level on the parameters.

[0147] The finite element analysis software (such as ANSYS, ABAQUS, etc.) is used to divide the grid of the mechanical model, to ensure that the grid quality meets the analysis requirements, such as grid size, shape and distribution. The grid model for finite element analysis is obtained.

[0148] The stress boundary conditions of the research area are determined by using single-well experimental test data (such as logging data, in-situ stress testing, etc.). The boundary conditions are set in the finite element analysis software, and the simulation is run to calculate the stress field distribution in the three-dimensional space. The simulation results, including the horizontal maximum principal stress, the horizontal minimum principal stress and the vertical principal stress data, are output.

[0149] The measured stress values of the cores are obtained by using rock acoustic emission method, differential strain method and other methods. The measured stress values are compared with the numerical simulation results. By adjusting the boundary load or boundary displacement parameters, the simulation results are made consistent with the measured values, with an error controlled within 10%, so as to obtain an accurate three-dimensional stress field distribution characteristic model.

[0150] S307, logging is performed on 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.

[0151] In one example, step S307 includes the following process: based on the pre-acquired imaging logging technology, the drilling-induced fracture phenomenon and the wellbore collapse phenomenon on the faults of the target layer system to be identified are identified, and 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 are obtained.

[0152] Exemplarily, this step can refer to step S104 described above, and will not be described again.

[0153] In one example, 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 are obtained by a paleomagnetic experiment method and a wave velocity anisotropy experiment method. The paleomagnetic experiment method analyzes the direction of the ancient magnetic field preserved in the rock, and in combination with the geological history, the stress state and direction of the rock when it is formed can be inferred. The wave velocity anisotropy experiment method measures the difference in the propagation speed of sound waves in the rock in different directions to infer the stress distribution and the principal stress direction inside 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 direction of the horizontal maximum principal stress and the direction of the horizontal minimum principal stress, and the starting point of the direction of the vertical principal stress is located at the intersection of the orthogonal plane.

[0154] S308, determining a slip index value of the grid according to the grid point occurrence of the grid, the pore fluid pressure value of the grid face of the grid, the principal stress data of the grid, the direction of the horizontal maximum principal stress of the grid, the direction of the horizontal minimum principal stress of the grid, and the direction of the vertical principal stress of the grid, wherein the slip index value represents whether the grid face of the grid will slip.

[0155] In one example, the grid point occurrence of the grid includes the grid point strike of the grid on the fault, and the grid point strike is the direction of the connecting line from the grid point to the first adjacent grid point adjacent to the grid point, and the first adjacent grid point is the next grid point on the same grid line as the grid point in the first direction, and the first direction is the horizontal direction in the division direction of the grid in the three-dimensional space grid model. Step S208 includes the following processes:

[0156] determining the included angle between the grid point strike of the grid and the north direction as the included angle value of the grid;

[0157] determining the cosine value of the first included angle of the grid and the cosine value of the second included angle of the grid according to the grid point dip angle of the grid and the included angle value of the grid, and determining the cosine value of the third included angle of the grid according to the grid point dip angle of the grid; wherein the first included angle of the grid is the included angle between the normal of the grid face and the direction of the horizontal maximum principal stress of the grid, the second included angle of the grid is the included angle between the normal of the grid face and the direction of the horizontal minimum principal stress of the grid, and the third included angle of the grid is the included angle between the normal of the grid face and the direction of the vertical principal stress of the grid;

[0158] determining the effective normal stress value of the grid according to the pore fluid pressure value of the grid face of the grid, the horizontal maximum principal stress of the grid, the horizontal minimum principal stress of the grid, the vertical principal stress of the grid, the cosine value of the first included angle of the grid, the cosine value of the second included angle of the grid, and the cosine value of the third included angle of the grid; wherein the effective normal stress value of the grid represents the effective normal stress borne by the grid face of the grid;

[0159] The shear stress value of the grid is determined according to the horizontal maximum principal stress of the grid, the horizontal minimum principal stress of the grid, the vertical principal stress of the grid, the cosine value of the first included angle of the grid, the cosine value of the second included angle of the grid, and the cosine value of the third included angle of the grid; wherein the shear stress of the grid represents the shear stress borne by the grid surface of the grid;

[0160] The slip index value of the grid is determined according to the effective normal stress value of the grid and the shear stress of the grid.

[0161] Exemplarily, the normal direction of the grid surface is determined, in the three-dimensional space grid model, each grid surface has a unique normal direction, which is perpendicular to the grid surface. The first included angle refers to the included angle between the normal of the grid surface and the horizontal maximum principal stress direction, the second included angle refers to the included angle between the normal of the grid surface and the horizontal minimum principal stress direction, and the third included angle refers to the included angle between the normal of the grid surface and the vertical principal stress direction.

[0162] In one example, FIG. 4 is a schematic diagram of forces borne by grid points of a grid of a three-dimensional space grid model provided by the present application to identify a target fault layer system, as shown in FIG. 4, line segment is the normal direction of the grid surface; line segment

[0163] is the maximum principal stress direction; line segment is the minimum principal stress direction; line segment is the vertical stress direction. The first included angle β1 of the grid is ∠DOB; the second included angle β2 of the grid is ∠DOA; and the second included angle β3 of the grid is ∠DOC. Line segment is the force in the direction of the finite principal stress of the grid surface; line segment is the force in the direction of the maximum principal stress value minus the pore fluid pressure value, that is, S i -P; line segment is the force in the direction of the minimum principal stress value minus the pore fluid pressure value, that is, S h -P; line segment is the force in the direction of the vertical principal stress value minus the pore fluid pressure value, that is, S v -P.

[0164] In one example, the cosine value of the first included angle of the grid is cosβ1=sinαsinγ; the cosine value of the second included angle of the grid is cosβ2=sinαcosγ; and the cosine value of the third included angle of the grid is cosβ1=sinαsinγ; wherein α is the grid point inclination angle of the grid, and γ is the included angle value of the grid. 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.

[0165] The shear stress of the grid is

[0166] Wherein, S i is the horizontal maximum principal stress of the grid, S h is the horizontal minimum principal stress of the grid, S v is the vertical principal stress of the grid, and P is the pore fluid pressure value of the grid face of the grid.

[0167] In an example, the slip index value of the grid is H = τ / (μσ+C); wherein σ is the effective normal stress value of the grid, τ is the shear stress of the grid, μ is the preset friction coefficient of the grid face of the grid, and C is a preset constant.

[0168] S309, according to the number of grids that slip, determining the slip risk value of the fault of the target layer system to be identified; wherein the slip risk value represents the risk coefficient of the slip of the fault of the target layer system to be identified.

[0169] Exemplarily, the present step refers to the above step S106, and will not be repeated here.

[0170] S310, if the slip risk value of the fault of the target layer system to be identified represents that the risk coefficient of the slip of the fault of the target layer system to be identified is a high risk coefficient or a medium risk coefficient, a warning information is sent to the user; wherein the warning information represents that the fault of the target layer system to be identified will slip.

[0171] Exemplarily, the present step refers to the above step S106, and will not be repeated here.

[0172] The method for identifying the slip risk of the fault provided by the embodiments of the present application establishes the pore fluid pressure and formation velocity relationship model through the Bowers method, the Fillippone method and the like; obtains the higher-precision seismic layer velocity by using the logging data to constrain the seismic inversion; performs the three-dimensional space formation pore pressure calculation by using the formation pore pressure relationship model established by the logging and the seismic inversion layer velocity, so that the pore fluid pressure value of each grid can be obtained more accurately. Through the finite element software analysis, the stress field numerical simulation analysis is performed, the three-dimensional stress field data are obtained, the measured stress value of the core is obtained through the rock acoustic emission method, the differential strain method and the like, the measured value and the simulated value are corrected by using the method of adjusting the boundary load or the boundary displacement, so that the accurate three-dimensional stress field distribution characteristic model is obtained; all the obtained numerical values are integrated, the slip risk value of each grid is calculated. According to the slip risk value of each grid, the slip risk coefficient of the fault to be identified is obtained, so as to provide a scientific basis for the deep shale gas development.

[0173] Figure 5 is a structural schematic diagram of a fault slip risk identification device provided by the present application. As shown in Figure 5, the fault slip risk identification device 50 provided by the present embodiment comprises:

[0174] The establishing module 501 is configured to establish a three-dimensional space grid model of the fault of the target layer series to be identified, wherein the three-dimensional space grid model comprises grid point occurrences of grids on the fault of the target layer series to be identified, and the grid point occurrences comprise grid point dips and grid point strikes of the grids on the fault, the grid point dip represents a dip angle of a grid point on a grid, and the grid point strike represents a strike of the grid point on the grid.

[0175] The obtaining module 502 is configured to obtain a pore fluid pressure value of a grid face of the grid, and the pore fluid pressure value of the grid face of the grid represents a pressure borne by a fluid in a pore on the grid face.

[0176] The first determining module 503 is configured to determine main stress data on the grid face of the grid according to a preset three-dimensional tectonic geological model, wherein the preset three-dimensional tectonic geological model represents horizon information and fault information of the target layer series to be identified, the main stress data represents a main stress borne by the grid face, and the main stress data comprises a horizontal maximum principal stress of the grid, a horizontal minimum principal stress of the grid, and a vertical principal stress of the grid.

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

[0178] The second determining module 505 is configured to determine a slip index value of the grid according to the grid point occurrences of the grid, the pore fluid pressure value of the grid face of the grid, the main stress data of the grid, the direction of the horizontal maximum principal stress of the grid, the direction of the horizontal minimum principal stress of the grid, and the direction of the vertical principal stress of the grid, wherein the slip index value represents whether the grid face of the grid will slip.

[0179] The third determining module 506 is configured to determine a slip risk value of the fault of the target layer series to be identified according to a number of grids that have slipped, wherein the slip risk value represents a risk coefficient of the fault of the target layer series to be identified slipping.

[0180] Figure 6 is a structural schematic diagram of a fault slip risk identification device provided by the present application. As shown in Figure 6, the fault slip risk identification device 60 provided by the present embodiment comprises:

[0181] The establishing module 601 is configured to establish a three-dimensional space grid model of faults of a target layer system to be identified, wherein the three-dimensional space grid model comprises grid point occurrences of grids on the faults of the target layer system to be identified, and the grid point occurrences comprise grid point dips and grid point strikes of the grids on the faults, the grid point dip represents a dip angle of a grid point on a grid, and the grid point strike represents a strike of the grid point on the grid.

[0182] The acquiring module 602 is configured to acquire a pore fluid pressure value of a grid face of the grid, and the pore fluid pressure value of the grid face of the grid represents a pressure borne by a fluid in a pore on the grid face.

[0183] The first determining module 603 is configured to determine main stress data on the grid face of the grid according to a preset three-dimensional tectonic geological model, wherein the preset three-dimensional tectonic geological model represents horizon information and fault information of the target layer system to be identified, the main stress data represents a main stress borne by the grid face, and the main stress data comprises a horizontal maximum main stress of the grid, a horizontal minimum main stress of the grid, and a vertical main stress of the grid.

[0184] The logging module 604 is configured to log the grid to obtain a direction of the horizontal maximum main stress of the grid, a direction of the horizontal minimum main stress of the grid, and a direction of the vertical main stress of the grid.

[0185] The second determining module 605 is configured to determine a slip index value of the grid according to the grid point occurrences of the grid, the pore fluid pressure value of the grid face of the grid, the main stress data of the grid, the direction of the horizontal maximum main stress of the grid, the direction of the horizontal minimum main stress of the grid, and the direction of the vertical main stress of the grid, wherein the slip index value represents whether the grid face of the grid will slip.

[0186] The third determining module 606 is configured to determine a slip risk value of the faults of the target layer system to be identified according to a number of grids that slip, wherein the slip risk value represents a risk coefficient of the faults of the target layer system to be identified slipping.

[0187] In one example, the grid point strike of the grid on the fault is a direction of a line connecting the grid point to a first adjacent grid point adjacent to the grid point, the first adjacent grid point adjacent to the grid point is a next grid point on a same grid line as the grid point in a first direction, and the first direction is a horizontal direction in a division direction of the grid in the three-dimensional space grid model.

[0188] The second determining module 605 comprises a first determining sub-module 6051 configured to determine an included angle between the grid point strike of the grid and a due north direction as the included angle value of the grid.

[0189] The second determining sub-module 6052 is configured to determine a cosine value of a first included angle of the grid and a cosine value of a second included angle of the grid according to the grid point dip angle of the grid and the included angle value of the grid, and determine a cosine value of a third included angle of the grid according to the grid point dip angle of the grid; wherein the first included angle of the grid is an included angle between a normal of a grid face and a direction of a horizontal maximum principal stress of the grid, the second included angle of the grid is an included angle between the normal of the grid face and a direction of a horizontal minimum principal stress of the grid, and the third included angle of the grid is an included angle between the normal of the grid face and a direction of a vertical principal stress of the grid.

[0190] The third determining sub-module 6053 is configured to determine an effective normal stress value of the grid according to the pore fluid pressure value of the grid face of the grid, the horizontal maximum principal stress of the grid, the horizontal minimum principal stress of the grid, the vertical principal stress of the grid, the cosine value of the first included angle of the grid, the cosine value of the second included angle of the grid, and the cosine value of the third included angle of the grid; wherein the effective normal stress value of the grid represents an effective normal stress borne by the grid face of the grid.

[0191] The fourth determining sub-module 6054 is configured to determine a shear stress value of the grid according to the horizontal maximum principal stress of the grid, the horizontal minimum principal stress of the grid, the vertical principal stress of the grid, the cosine value of the first included angle of the grid, the cosine value of the second included angle of the grid, and the cosine value of the third included angle of the grid; wherein the shear stress of the grid represents a shear stress borne by the grid face of the grid.

[0192] The fifth determining sub-module 6055 is configured to determine a slip index value of the grid according to the effective normal stress value of the grid and the shear stress of the grid.

[0193] In one example, the cosine value of the first included angle of the grid is cosβ1=sinαsinγ; the cosine value of the second included angle of the grid is cosβ2=sinαcosγ; and the cosine value of the third included angle of the grid is cosβ1=sinαsinγ; wherein α is the grid point dip angle of the grid, and γ is the included angle value of the grid.

[0194] In one example, 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.

[0195] The shear stress of the grid is

[0196] wherein S i is the horizontal maximum principal stress of the grid, S h is the horizontal minimum principal stress of the grid, and Sv P is the pore fluid pressure value of the grid face of the grid. The slip index value of the grid is H = τ / (μσ + C); wherein, σ is the effective normal stress value of the grid, τ is the shear stress of the grid, μ is the preset friction coefficient of the grid face of the grid, and C is a preset constant.

[0197] In one example, the establishing module 601 comprises:

[0198] The first establishing submodule 6011 is configured to establish a three-dimensional spatial distribution feature according to a three-dimensional seismic data volume; wherein, the three-dimensional seismic data volume comprises data of a fault of a target layer system to be identified; and the three-dimensional spatial distribution feature represents a feature of the fault of the target layer system to be identified in a three-dimensional space.

[0199] The dividing submodule 6012 is configured to divide the three-dimensional spatial distribution feature into grids to obtain a grid of the fault of the target layer system to be identified; and determine a connecting line of the grid; wherein, the connecting line of the grid is a line connecting a grid point to a second adjacent grid point adjacent to the grid point, and the second adjacent grid point adjacent to the grid point is a next grid point located on a same grid line as the grid point in a second direction, and the second direction is a vertical direction in a division direction of the grid in the three-dimensional spatial grid model.

[0200] The first generating submodule 6013 is configured to determine an included angle between the connecting line of the grid and a physical horizontal line as a grid point dip angle of the grid to obtain the three-dimensional spatial grid model.

[0201] In one example, the obtaining module 602 comprises:

[0202] The obtaining submodule 6021 is configured to obtain logging data, wherein the logging data comprises a formation velocity; and establish a relationship model according to a preset measured pore fluid pressure value and the formation velocity in the logging data, wherein the relationship model is a corresponding relationship between the pore fluid pressure value and the formation velocity.

[0203] The first processing submodule 6022 is configured to obtain three-dimensional seismic material data, and perform seismic inversion calculation processing on the three-dimensional seismic material data using the logging data to obtain a seismic interval velocity of a grid face of a grid of the fault of the target layer system to be identified.

[0204] The second processing submodule 6023 is configured to input the seismic interval velocity of the grid face of the grid as the formation velocity into the relationship model to obtain a pore fluid pressure value of the grid face of the grid.

[0205] In one example, the first determining module 603 comprises:

[0206] The second establishing submodule 6035 is configured to acquire three-dimensional seismic material data of a target layer system to be identified, wherein the three-dimensional seismic material data comprises horizon interpretation data and fault interpretation data of the target layer system; and a preset three-dimensional structural geological model is established according to the three-dimensional seismic material data.

[0207] The detecting submodule 6031 is configured to detect a mechanical parameter of the target layer system to be identified, wherein the mechanical parameter represents a parameter that the target layer system has in mechanics.

[0208] The second generating submodule 6032 is configured to assign a parameter in the mechanical parameter of the target layer system to be identified to the preset three-dimensional structural geological model to obtain a three-dimensional mechanical model; and 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 meshes.

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

[0210] The correction submodule 6034 is configured to correct the initial principal stress data on the mesh surface of the mesh based on a core measured stress value of the target layer system to be identified, to obtain principal stress data on the mesh surface of the mesh, wherein the core measured stress value represents a principal stress that a core of the target layer system bears.

[0211] The logging module 604 is specifically configured to:

[0212] Based on the imaging logging technology acquired in advance, the drilling-induced fracture phenomenon and the wellbore collapse phenomenon on the fault of the target layer system to be identified are identified to obtain a direction of the horizontal maximum principal stress, a direction of the horizontal minimum principal stress, and a direction of the vertical principal stress of the mesh.

[0213] In one example, the apparatus provided by the embodiment of the present application further includes:

[0214] The warning module 607 is configured to issue a warning information to a user if the slip risk value of the fault of the target layer system to be identified represents a high risk coefficient or a medium risk coefficient of the risk of the fault of the target layer system to be identified to slip, wherein the warning information represents that the fault of the target layer system to be identified will slip.

[0215] The apparatus for identifying the slip risk of the fault provided by the embodiment of the present application can execute the method provided by the method embodiment, and has similar implementation principles and technical effects, which will not be described here in detail.

[0216] Fig. 7 is a structural schematic diagram of an electronic device for identifying a slip risk of a fault provided in the present application. As shown in Fig. 7, the electronic device 70 provided in the embodiment includes at least one processor 701 and a memory 702. Optionally, the device 70 further includes a communication component 703. The processor 701, the memory 702 and the communication component 703 are connected through a bus 704.

[0217] In the implementation process, the at least one processor 701 executes the computer-executable instructions stored in the memory 702, so that the at least one processor 701 performs the above-mentioned method.

[0218] The specific implementation process of the processor 701 can refer to the method embodiments described above, which have similar implementation principles and technical effects, and will not be described here in detail.

[0219] In the above embodiments, it should be understood that the processor can be a central processing unit (CPU), and can also be other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), etc. The general-purpose processor can be a microprocessor or any conventional processor, etc. The steps of the method disclosed in the application can be directly embodied as execution completed by a hardware processor, or executed by a combination of hardware and software modules in the processor.

[0220] The memory can include a random access memory (RAM), and can also include a non-volatile memory (NVM), such as at least one disk memory.

[0221] The bus can be an industry standard architecture (ISA) bus, a peripheral component (PCI) bus or an extended industry standard architecture (EISA) bus, etc. The bus can be divided into an address bus, a data bus, a control bus, etc. For the convenience of representation, the bus in the drawings of the present application does not limit to only one bus or one type of bus.

[0222] The present application also provides a computer program product, including a computer program, which is executed by a processor to implement the above-mentioned method.

[0223] The application further provides a computer readable storage medium, wherein computer execution instructions are stored in the computer readable storage medium, and when a processor executes the computer execution instructions, the method described above is realized.

[0224] The readable storage medium described above can be realized 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 memory, flash memory, magnetic disk or optical disk. The readable storage medium can be any available medium that can be accessed by a general-purpose or special-purpose computer.

[0225] An exemplary readable storage medium is coupled to the processor, so that the processor can read information from the readable storage medium and write information to the readable storage medium. Of course, the readable storage medium can also be an integral part of the processor. The processor and the readable storage medium can be located in an application specific integrated circuit (ASIC). Of course, the processor and the readable storage medium can also exist as discrete components in the device.

[0226] The division of units is only a logical function division, and in actual implementation, there can be another division mode, for example, a plurality of units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the coupling or direct coupling or communication connection between the units or components shown or discussed can be indirect coupling or communication connection through some interfaces, devices or units, which can be electrical, mechanical or other forms.

[0227] The units described as separate components can or can not be physically separated, and the components shown as units can or can not be physical units, that is, they can be located in one place, or they can be distributed on a plurality of network units. According to actual needs, some or all of the units can be selected to achieve the purpose of the embodiment.

[0228] In addition, the functional units in each embodiment of the application can be integrated in one processing unit, or each unit can exist physically, or two or more units can be integrated in one unit.

[0229] If the functions are implemented in the form of software function units and sold or used as independent products, they can be stored in a computer readable storage medium. Based on this understanding, the technical solutions of the present application or the parts of the technical solutions that essentially contribute to the prior art can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes a number of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the embodiments of the method of the present application. The aforementioned storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various media that can store program codes.

[0230] It can be understood by those skilled in the art that all or part of the steps of the above-mentioned method embodiments can be completed by program instruction related hardware. The aforementioned program can be stored in a computer readable storage medium. When the program is executed, the steps of the above-mentioned method embodiments are executed; and the aforementioned storage medium includes: ROM, RAM, magnetic disk or optical disk, and various media that can store program codes.

[0231] Finally, it should be noted that: those skilled in the art will easily think of other embodiments of the present application after considering the specification and practicing the application disclosed herein. The present application is intended to cover any variations, uses or adaptations of the present application that follow the general principles of the present application and include common knowledge or conventional technical means in the art that are not disclosed in the present application, and is not limited to the precise structure described above and shown in the drawings, and various modifications and changes can be made without departing from the scope thereof. The scope of the present application is only limited by the appended claims.

Claims

1. A method of identifying a slip risk of a fault, characterized by, The method comprises the following steps: establishing a three-dimensional space grid model of faults of a target layer system to be identified, wherein the three-dimensional space grid model comprises grid point occurrences of grids on the faults of the target layer system to be identified, the grid point occurrences comprise a grid point dip angle of the grids on the faults and a grid point strike of the grids, the grid point dip angle represents a dip angle of a grid point on the grid, and the grid point strike represents a strike of the grid point on the grid; and obtaining a pore fluid pressure value of a grid face of the grid, which represents a pressure borne by fluid in a pore on the grid face; determining main stress data on the grid face of the grid according to a preset three-dimensional tectonic geological model, wherein the preset three-dimensional tectonic geological model represents horizon information and fault information of the target layer system to be identified, the main stress data represents a main stress borne by the grid face, and the main stress data comprises a horizontal maximum main stress of the grid, a horizontal minimum main stress of the grid, and a vertical main stress of the grid; logging the grid to obtain a direction of the horizontal maximum main stress of the grid, a direction of the horizontal minimum main stress of the grid, and a direction of the vertical main stress of the grid; determining a slip index value of the grid according to the grid point occurrences of the grid, the pore fluid pressure value of the grid face of the grid, the main stress data of the grid, the direction of the horizontal maximum main stress of the grid, the direction of the horizontal minimum main stress of the grid, and the direction of the vertical main stress of the grid, wherein the slip index value represents whether the grid face of the grid will slip; determining a slip risk value of the faults of the target layer system to be identified according to a number of grids that slip, wherein the slip risk value represents a risk coefficient of the faults of the target layer system to be identified slipping.

2. The method of claim 1, wherein, The grid point strike of the grid on the fault is a direction of a line connecting the grid point to a first adjacent grid point adjacent to the grid point, and the first adjacent grid point adjacent to the grid point is a next grid point on a same grid line as the grid point in a horizontal direction of a division direction of the grid in the three-dimensional space grid model; determining the slip index value of the grid according to the grid point occurrences of the grid, the pore fluid pressure value of the grid face of the grid, the main stress data of the grid, the direction of the horizontal maximum main stress of the grid, the direction of the horizontal minimum main stress of the grid, and the direction of the vertical main stress of the grid comprises: determining an included angle between the grid point strike of the grid and a north direction as an included angle value of the grid; determining a cosine value of a first included angle of the grid and a cosine value of a second included angle of the grid according to the grid point dip angle of the grid and the included angle value of the grid, and determining a cosine value of a third included angle of the grid according to the grid point dip angle of the grid, wherein the first included angle of the grid is an included angle between a normal line of the grid face and the direction of the horizontal maximum main stress of the grid, the second included angle of the grid is an included angle between the normal line of the grid face and the direction of the horizontal minimum main stress of the grid, and the third included angle of the grid is an included angle between the normal line of the grid face and the direction of the vertical main stress of the grid; Determine the effective normal stress value of the grid according to the pore fluid pressure value of the grid face of the grid, the horizontal maximum principal stress of the grid, the horizontal minimum principal stress of the grid, the vertical principal stress of the grid, the cosine value of the first included angle of the grid, the cosine value of the second included angle of the grid and the cosine value of the third included angle of the grid; wherein the effective normal stress value of the grid represents the effective normal stress borne by the grid face of the grid; Determine the shear stress value of the grid according to the horizontal maximum principal stress of the grid, the horizontal minimum principal stress of the grid, the vertical principal stress of the grid, the cosine value of the first included angle of the grid, the cosine value of the second included angle of the grid and the cosine value of the third included angle of the grid; wherein the shear stress of the grid represents the shear stress borne by the grid face of the grid; Determine the slip index value of the grid according to the effective normal stress value of the grid and the shear stress of the grid.

3. The method of claim 2, wherein, The cosine value of the first included angle of the grid is cosβ1=sinαsinγ; the cosine value of the second included angle of the grid is cosβ2=sinαcosγ; and the cosine value of the third included angle of the grid is cosβ1=sinαsinγ; wherein α is the grid point inclination of the grid, and γ is the included angle value of the grid.

4. The method of claim 2, wherein, 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 wherein S i is the horizontal maximum principal stress of the grid, S h is the horizontal minimum principal stress of the grid, S v is the vertical principal stress of the grid, P is the pore fluid pressure value of the grid face of the grid; the slip index value of the grid is H = τ / (μσ + C); wherein σ is the effective normal stress value of the grid, τ is the shear stress of the grid, μ is the preset friction coefficient of the grid face of the grid, and C is a preset constant.

5. The method according to any one of claims 1 to 4, characterized in that, Establish a three-dimensional space grid model of the fault of the target layer system to be identified, including: According to the three-dimensional seismic data body, establish a three-dimensional space distribution characteristic; wherein the three-dimensional seismic data body includes the data of the fault of the target layer system to be identified; and the three-dimensional space distribution characteristic represents the characteristics of the fault of the target layer system to be identified in the three-dimensional space; Grid division is performed on the three-dimensional space distribution characteristic to obtain the grid of the fault of the target layer system to be identified, and the connection line of the grid is determined; wherein the connection line of the grid is a 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 vertical direction in the division direction of the grid in the three-dimensional space grid model; An included angle between the connection line of the grid and the physical horizontal line is determined as the grid point inclination of the grid to obtain the three-dimensional space grid model.

6. The method according to any one of claims 1 to 5, characterized in that, Obtaining the pore fluid pressure value of the grid face of the grid includes: Obtaining logging data, wherein the logging data includes formation velocity; and establishing a relationship model according to the measured pore fluid pressure value and the formation velocity in the logging data, wherein the relationship model is the corresponding relationship between the pore fluid pressure value and the formation velocity; Obtaining three-dimensional seismic material data, performing seismic inversion calculation and processing on the three-dimensional seismic material data using the logging data to obtain the seismic interval velocity of the grid face of the grid of the fault of the target layer system to be identified; Inputting the seismic interval velocity of the grid face of the grid as the formation velocity into the relationship model to obtain the pore fluid pressure value of the grid face of the grid.

7. The method according to any one of claims 1 to 6, characterized in that, According to a preset three-dimensional structural geological model, determine the principal stress data on the grid face of the grid, including: acquire three-dimensional seismic material data of the target layer system to be identified, wherein the three-dimensional seismic material data comprises horizon interpretation data and fault interpretation data of the target layer system; establish a preset three-dimensional structural geological model according to the three-dimensional seismic material data; detect a mechanical parameter of the target layer system to be identified, wherein the mechanical parameter represents a parameter of the target layer system in mechanics; assign a parameter in the mechanical parameter of the target layer system to be identified to the preset three-dimensional structural geological model to obtain a three-dimensional mechanical model; perform grid division on the three-dimensional mechanical model to obtain a grid model for finite element analysis, wherein the grid model for finite element analysis is a three-dimensional mechanical model with grids; perform stress field numerical simulation processing on the grid model for finite element analysis to obtain initial principal stress data on a grid surface of the grid; correct the initial principal stress data on the grid surface of the grid based on a core measured stress value of the target layer system to be identified, to obtain principal stress data on the grid surface of the grid, wherein the core measured stress value represents a principal stress borne by a core of the target layer system; perform logging on the grid to obtain a direction of a horizontal maximum principal stress, a direction of a horizontal minimum principal stress, and a direction of a vertical principal stress of the grid, including: identify a drilling-induced fracture phenomenon and a borehole wall collapse phenomenon on a fault of the target layer system to be identified based on pre-acquired imaging logging data 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.

8. The method according to any one of claims 1 to 7, characterized in that, The method further includes: if a slip risk value of the fault of the target layer system to be identified represents a high risk coefficient or a medium risk coefficient of a risk of slip of the fault of the target layer system to be identified, issuing a warning information to a user, wherein the warning information represents that the fault of the target layer system to be identified will slip.

9. An apparatus for identifying a slip risk of a fault, characterized by The apparatus includes modules for performing the method of any one of claims 1-8.

10. An electronic device, comprising: including: a memory, a processor; the memory stores computer execution instructions; the processor executes the computer execution instructions stored in the memory, so that the processor performs the method of 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 execution instructions, and the computer execution instructions are executed by the processor to implement the method of any one of claims 1-8; and The computer program product includes a computer program, and the computer program is executed by the processor to implement the method of any one of claims 1-8.

Images