Method, device, equipment and medium for determining the maximum tolerable gas column height of a fault

Abstract

The present application provides a method, device, equipment, and medium for determining the maximum tolerable gas column height of a fault, wherein the method comprises: obtaining the vertical principal stress, the maximum horizontal principal stress, and the minimum horizontal principal stress corresponding to the fault in a hydrate accumulation zone at different depths; the hydrate accumulation zone is a shallow unconsolidated seepage hydrate accumulation zone; according to the vertical principal stress, the maximum horizontal principal stress, and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths, determining the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone; wherein the minimum sliding pressure is used to indicate the minimum pressure value required to make the fault in the hydrate accumulation zone slide; according to the minimum sliding pressure, determining the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation zone, can improve the accuracy of determining the maximum tolerable gas column height of the fault, thereby improving the safety during hydrate exploitation and reducing the probability of geological disasters caused by formation instability.

Description

Technical Field

[0001] The present application relates to the technical field of geological exploration and development of oil and gas reservoirs, and in particular to a method, device, equipment, and medium for determining the maximum tolerable gas column height of a fault. Background Art

[0002] Hydrates are usually found in permafrost on land or shallow marine sediments and are considered to be one of the important clean energy sources in the future. These areas have highly sensitive and unstable geological structures, especially the faults in the hydrate accumulation area. Once the height of the gas column accumulated in the fault exceeds the tolerance of the fault, the fault will slide, which may cause stratum instability and trigger geological disasters such as landslides, ground subsidence and even submarine mudslides. This not only poses a threat to mining activities, but may also have a serious impact on the surrounding environment.

[0003] The existing technology has the problem of low accuracy when determining the maximum tolerable gas column height of a fault.

[0004] Therefore, there is an urgent need for a method for determining the maximum tolerable air column height of a fault to improve the accuracy of determining the maximum tolerable air column height of a fault. Summary of the invention

[0005] The present application provides a method, device, equipment, and medium for determining the maximum tolerable air column height of a fault, so as to improve the accuracy of determining the maximum tolerable air column height of a fault.

[0006] In a first aspect, an embodiment of the present application provides a method for determining a maximum tolerable gas column height of a fault, comprising:

[0007] Obtaining the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths; the hydrate accumulation zone is a shallow hydrate accumulation zone under a tensile stress state;

[0008] Determine the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone according to the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths; wherein the minimum sliding pressure is used to indicate the minimum pressure value required to make the fault in the hydrate accumulation zone slide;

[0009] According to the minimum sliding pressure, the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation zone is determined.

[0010] Optionally, the vertical principal stress, the maximum horizontal principal stress, and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths all include magnitudes and directions. According to the vertical principal stress, the maximum horizontal principal stress, and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths, determining the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone includes:

[0011] Acquire a three-dimensional model of the fault in the hydrate accumulation area, and determine the corresponding dip angles of the fault in the hydrate accumulation area at different depths according to the three-dimensional model;

[0012] Determine the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths according to the inclination angle, vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths;

[0013] According to the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths, the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone is determined.

[0014] Optionally, determining the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths according to the dip angle, vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths includes:

[0015] Determining the direction of the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths according to the inclination angle corresponding to the fault in the hydrate accumulation zone at different depths;

[0016] According to the direction of the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths, the direction of the vertical principal stress, the direction of the horizontal maximum principal stress and the direction of the horizontal minimum principal stress, determine the first angle between the vertical principal stress and the unit surface normal, the second angle between the horizontal maximum principal stress and the unit surface normal, and the third angle between the horizontal minimum principal stress and the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths;

[0017] According to the magnitude of the vertical principal stress, the magnitude of the horizontal maximum principal stress, the magnitude of the horizontal minimum principal stress, the first angle, the second angle and the third angle corresponding to the fault in the hydrate accumulation zone at different depths, the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths are determined.

[0018] Optionally, determining the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths according to the magnitude of the vertical principal stress, the magnitude of the horizontal maximum principal stress, the magnitude of the horizontal minimum principal stress, the first angle, the second angle, and the third angle corresponding to the fault in the hydrate accumulation zone at different depths includes:

[0019] The effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths are calculated using the following formula;

[0020] σ n =σ v cos 2 α n +σ H cos 2 β n +σ h cos 2 γ n

[0021]

[0022] Among them, σ n is the effective normal stress; τ is the shear stress; σ v is the vertical principal stress; σ H is the maximum horizontal principal stress; σ h is the horizontal minimum principal stress; α1 is the first angle; α2 is the second angle; α3 is the third angle.

[0023] Optionally, determining the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone according to the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths includes:

[0024] Obtain the friction coefficient corresponding to the fault in the hydrate accumulation area, and use the formula Determine the sliding pressure values ​​corresponding to the fault in the hydrate accumulation area at different depths; wherein, P s is the sliding pressure value, σ n is the effective normal stress; τ is the shear stress, μ f is the friction coefficient;

[0025] According to the sliding pressure values ​​corresponding to the fault in the hydrate accumulation zone at different depths, the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone is determined.

[0026] Optionally, determining the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation zone according to the minimum sliding pressure includes:

[0027] The maximum tolerable gas column height corresponding to the fault in the hydrate accumulation area is calculated using the following formula:

[0028]

[0029] In the formula, h g is the maximum tolerable air column height; ρ w is the seawater density in the hydrate accumulation area; ρ gP is the gas density accumulated in the fault of the hydrate accumulation area; g is the minimum sliding pressure corresponding to the fault in the hydrate accumulation area.

[0030] In a second aspect, an embodiment of the present application provides a device for determining a maximum tolerable gas column height of a fault, comprising:

[0031] An acquisition module is used to obtain the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation area at different depths; the hydrate accumulation area is a shallow hydrate accumulation area under a tensile stress state;

[0032] A first determination module is used to determine the minimum sliding pressure corresponding to the fault in the hydrate accumulation area according to the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation area at different depths; wherein the minimum sliding pressure is used to indicate the minimum pressure value required to make the fault in the hydrate accumulation area slide;

[0033] The second determination module is used to determine the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation area according to the minimum sliding pressure.

[0034] In a third aspect, an embodiment of the present application provides an electronic device, comprising: a processor, and a memory communicatively connected to the processor;

[0035] The memory stores computer-executable instructions;

[0036] The processor executes the computer-executable instructions stored in the memory to implement the method according to the first aspect.

[0037] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium, wherein the computer-readable storage medium stores computer-executable instructions, and when the computer-executable instructions are executed by a processor, they are used to implement the method described in the first aspect.

[0038] In a fifth aspect, an embodiment of the present application provides a computer program product, including a computer program, which implements the method described in the first aspect when executed by a processor.

[0039] The present application provides a method, device, equipment, and medium for determining the maximum tolerable gas column height of a fault, the method comprising: obtaining the vertical principal stress, the maximum horizontal principal stress, and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths; the hydrate accumulation zone is a shallow hydrate accumulation zone under a tensile stress state; according to the vertical principal stress, the maximum horizontal principal stress, and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths, determining the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone; wherein the minimum sliding pressure is used to indicate the minimum pressure value required to make the fault in the hydrate accumulation zone slide; according to the minimum sliding pressure, determining the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation zone can improve the accuracy of determining the maximum tolerable gas column height of the fault, thereby improving the safety during hydrate exploitation and reducing the probability of geological disasters caused by formation instability. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0041] Figure 1 An application scenario diagram provided for an embodiment of the present application;

[0042] Figure 2 A schematic flow chart of a method for determining the maximum tolerable gas column height of a fault provided in an embodiment of the present application;

[0043] Figure 3 A schematic diagram of a three-dimensional model of a hydrate accumulation zone fault provided in an embodiment of the present application;

[0044] Figure 4 A schematic diagram of the structure of a device for determining the maximum tolerable gas column height of a fault provided in an embodiment of the present application;

[0045] Figure 5 A schematic diagram of the structure of an electronic device provided in an embodiment of the present application.

[0046] The above drawings have shown clear embodiments of the present application, which will be described in more detail later. These drawings and text descriptions are not intended to limit the scope of the present application in any way, but to illustrate the concept of the present application to those skilled in the art by referring to specific embodiments. DETAILED DESCRIPTION

[0047] Exemplary embodiments will be described in detail herein, examples of which are shown in the accompanying drawings. When the following description refers to the drawings, the same numbers in different drawings represent the same or similar elements unless otherwise indicated. The implementations described in the following exemplary embodiments do not represent all implementations consistent with the present application. Instead, they are merely examples of devices and methods consistent with some aspects of the present application as detailed in the appended claims.

[0048] Hydrates are usually found in permafrost on land or shallow marine sediments and are considered to be one of the important clean energy sources in the future. These areas have highly sensitive and unstable geological structures, especially the faults in the hydrate accumulation area. Once the height of the gas column accumulated in the fault exceeds the tolerance of the fault, the fault will slide, which may cause stratum instability and trigger geological disasters such as landslides, ground subsidence and even submarine mudslides. This not only poses a threat to mining activities, but may also have a serious impact on the surrounding environment.

[0049] In addition, during the geological storage of carbon dioxide in the sea, if the injection pressure is too high and exceeds the upper pressure limit of the fault, the risk of carbon dioxide leakage will be greatly increased.

[0050] The existing technology has the problem of low accuracy when determining the maximum tolerable gas column height of a fault.

[0051] In view of this, the present application provides a method for determining the maximum tolerable gas column height of a fault, which is used to determine the maximum tolerable gas column height corresponding to the fault in the shallow hydrate accumulation zone. First, the vertical principal stress, the horizontal maximum principal stress and the horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths are obtained; according to the vertical principal stress, the horizontal maximum principal stress and the horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths, the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone is determined; finally, according to the minimum sliding pressure, the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation zone is determined. The present application determines the minimum sliding pressure corresponding to the fault according to the vertical principal stress, the horizontal maximum principal stress and the horizontal minimum principal stress corresponding to the fault at different depths, and then determines the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation zone according to the minimum sliding pressure. This can improve the accuracy of determining the maximum tolerable gas column height of the fault, thereby improving the safety during hydrate exploitation and reducing the probability of geological disasters caused by stratum instability.

[0052] Figure 1 An application scenario diagram provided for an embodiment of the present application, such as Figure 1As shown, the user inputs the vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault at different depths through the user device. The input form can be in the form of a function. For example, the user can input the function of vertical principal stress varying with depth, the function of horizontal maximum principal stress varying with depth and the function of horizontal minimum principal stress varying with depth through the user device. After the input is completed, click the OK button. The user device sends the function of vertical principal stress varying with depth, the function of horizontal maximum principal stress varying with depth and the function of horizontal minimum principal stress varying with depth input by the user to the server. After receiving the function, the server determines the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone according to the function of vertical principal stress varying with depth, the function of horizontal maximum principal stress varying with depth and the function of horizontal minimum principal stress varying with depth, and determines the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation zone according to the minimum sliding pressure, and sends the maximum tolerable gas column height corresponding to the fault to the user device, so that the user device displays the maximum tolerable gas column height corresponding to the fault on the display interface.

[0053] The technical solution of the present application and how the technical solution of the present application solves the above-mentioned technical problems are 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 may not be repeated in some embodiments. The embodiments of the present application will be described below in conjunction with the accompanying drawings.

[0054] Figure 2 The flowchart of a method for determining the maximum tolerable gas column height of a fault provided in an embodiment of the present application is shown in FIG. The execution subject of this embodiment can be any device with data processing function. This application uses a server as the execution subject for specific description. Figure 2 As shown, a method for determining the maximum tolerable gas column height of a fault provided by an embodiment of the present application may include:

[0055] Step 201, obtaining the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths; the hydrate accumulation zone is a shallow hydrate accumulation zone under a tensile stress state.

[0056] Specifically, the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress are respectively the compressive stresses that the rock is subjected to in three-dimensional space, and in this application, the compressive stress that the fault in the hydrate accumulation zone is subjected to in three-dimensional space.

[0057] The direction of the vertical principal stress is along the vertical direction. Under the tensile stress state, the vertical principal stress is the principal stress with the largest value among the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress.

[0058] The maximum horizontal principal stress and the minimum horizontal principal stress are both in the horizontal direction, and the directions of the maximum horizontal principal stress and the minimum horizontal principal stress are perpendicular to each other. Under the tensile stress state, the maximum horizontal principal stress is the principal stress with the middle value among the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress, and the minimum horizontal principal stress is the principal stress with the smallest value among the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress.

[0059] Among them, for a point at a certain depth in the fault of the hydrate accumulation area, the vertical principal stress at this point is formed by the sea water pressure above the seabed and the gravity of all the strata above this point.

[0060] The functional relationship between the vertical principal stress corresponding to the fault in the hydrate accumulation area and the depth is:

[0061]

[0062] Among them, σ v is the vertical principal stress, in MPa; h is the depth value corresponding to the depth to be calculated of the fault in the hydrate accumulation zone, in m; ρ(h) is the function of the target formation density varying with depth, which is obtained by density logging; the target formation includes the fault in the hydrate accumulation zone and the formation above the fault, and g is the gravitational acceleration, which can be taken as 9.80 m / s 2 ρ w is the seawater density above the hydrate accumulation area, in g / cm 3 ; h0 is the seawater depth above the hydrate accumulation area, in meters.

[0063] For the fault in the hydrate accumulation area, the cores in the fault were taken to test the rock mechanical properties and obtain the rock mechanical parameters. The rock mechanical parameters and vertical principal stress were substituted into the Huang model to obtain the horizontal minimum principal stress and the horizontal maximum principal stress.

[0064]

[0065] Among them, σ v is the vertical principal stress, in MPa; P wis the pore fluid pressure, in MPa; μ is the static Poisson's ratio; α is the Biot coefficient; β1 and β2 are tectonic stress coefficients. The static Poisson's ratio is an important parameter in material mechanics, which is used to describe the relationship between the lateral strain and the longitudinal strain of the material under stress. The Biot coefficient is a parameter used in rock mechanics and geomechanics to describe the relationship between the effective stress and the total stress of porous media (such as rock or soil) under external stress. The tectonic stress coefficient is an important parameter used in geomechanics and rock mechanics to describe the stress state of the crust, which reflects the characteristics of the tectonic stress field of the crust; the static Poisson's ratio and Biot coefficient are obtained by testing the lithology mechanical properties of core samples in the fault of the hydrate accumulation area; the tectonic stress coefficients β1 and β2 can be selected as empirical values ​​based on the actual situation of the hydrate accumulation area.

[0066] Among them, rock mechanical property testing refers to the testing and analysis of the mechanical properties of rock cores to understand the mechanical behavior and characteristics of rocks under different stress conditions. Among them, rock mechanical parameters can include: compressive strength, tensile strength, shear strength, elastic modulus, plastic modulus, Poisson's ratio, friction angle and failure criterion parameters.

[0067] Specifically, the user can input the vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths through the user device. The input can be in the form of a function. For example, the user can input the function of vertical principal stress varying with depth, the function of horizontal maximum principal stress varying with depth, and the function of horizontal minimum principal stress varying with depth on the input interface of the user device. After the input is completed, click the confirmation button, and the user device sends the function of vertical principal stress varying with depth, the function of horizontal maximum principal stress varying with depth, and the function of horizontal minimum principal stress varying with depth to the server, and the server obtains the function of vertical principal stress varying with depth, the function of horizontal maximum principal stress varying with depth, and the function of horizontal minimum principal stress varying with depth.

[0068] Step 202: determine the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone according to the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths; wherein the minimum sliding pressure is used to indicate the minimum pressure value required to make the fault in the hydrate accumulation zone slide.

[0069] Specifically, the server determines the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone based on the vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths; that is, the minimum pressure value required to cause the fault in the hydrate accumulation zone to slide.

[0070] Optionally, the vertical principal stress, the maximum horizontal principal stress, and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths all include magnitudes and directions. According to the vertical principal stress, the maximum horizontal principal stress, and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths, determining the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone includes:

[0071] Acquire a three-dimensional model of the fault in the hydrate accumulation area, and determine the corresponding dip angles of the fault in the hydrate accumulation area at different depths according to the three-dimensional model;

[0072] Determine the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths according to the inclination angle, vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths;

[0073] According to the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths, the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone is determined.

[0074] Among them, the three-dimensional model of the fault in the hydrate accumulation area can be obtained by interpreting the seismic data corresponding to the fault in the hydrate accumulation area.

[0075] For each point on the fault in the hydrate accumulation zone, the corresponding section of the point is the surface tangent to the fault. The effective normal stress corresponding to the point refers to the normal stress perpendicular to the corresponding section of the point, and the shear stress corresponding to the point refers to the stress parallel to the corresponding section of the point.

[0076] Specifically, a three-dimensional model of the fault in the hydrate accumulation zone is obtained, and the inclination angle corresponding to the fault in the hydrate accumulation zone at different depths is determined according to the three-dimensional model; the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths are determined according to the inclination angle corresponding to the fault in the hydrate accumulation zone at different depths, the magnitude and direction of the vertical principal stress, the magnitude and direction of the maximum horizontal principal stress, and the magnitude and direction of the minimum horizontal principal stress; finally, the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone is determined according to the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths.

[0077] In this way, the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths can be determined according to the inclination angle, vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths, which can improve the accuracy of the determined effective normal stress and shear stress, and then improve the accuracy of the minimum sliding pressure corresponding to the determined fault.

[0078] Optionally, determining the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths according to the dip angle, vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths includes:

[0079] Determining the direction of the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths according to the inclination angle corresponding to the fault in the hydrate accumulation zone at different depths;

[0080] According to the direction of the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths, the direction of the vertical principal stress, the direction of the horizontal maximum principal stress and the direction of the horizontal minimum principal stress, determine the first angle between the vertical principal stress and the unit surface normal, the second angle between the horizontal maximum principal stress and the unit surface normal, and the third angle between the horizontal minimum principal stress and the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths;

[0081] According to the magnitude of the vertical principal stress, the magnitude of the horizontal maximum principal stress, the magnitude of the horizontal minimum principal stress, the first angle, the second angle and the third angle corresponding to the fault in the hydrate accumulation zone at different depths, the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths are determined.

[0082] Among them, for each depth point in the fault of the hydrate accumulation area, the section of the fault at this point is the unit surface, and the normal of the unit surface is perpendicular to the unit surface, that is, the inclination plus 90 degrees is the angle between the normal of the unit surface and the horizontal plane.

[0083] Specifically, the server determines the direction of the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths according to the inclination angle corresponding to the fault in the hydrate accumulation zone at different depths, and then determines the angle between the vertical principal stress corresponding to the fault in the hydrate accumulation zone at different depths and the unit surface normal, i.e., the first angle, according to the direction of the unit surface normal, the direction of the vertical principal stress, the direction of the horizontal maximum principal stress, and the direction of the horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths; determines the angle between the horizontal maximum principal stress corresponding to the fault in the hydrate accumulation zone at different depths and the unit surface normal, i.e., the second angle; determines the angle between the horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths and the unit surface normal, i.e., the third angle.

[0084] Finally, the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths are determined according to the magnitude of the vertical principal stress, the magnitude of the horizontal maximum principal stress, the magnitude of the horizontal minimum principal stress, the first angle, the second angle and the third angle corresponding to the fault in the hydrate accumulation zone at different depths.

[0085] In this way, the direction of the unit surface normal is determined according to the corresponding inclination angles of the fault at different depths, and then the three angles are determined according to the direction of the unit surface normal. Finally, the effective normal stress and shear stress are determined according to the magnitude of the vertical principal stress, the horizontal maximum principal stress and the horizontal minimum principal stress and the three angles. This can improve the accuracy of the determined effective normal stress and shear stress.

[0086] Optionally, determining the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths according to the magnitude of the vertical principal stress, the magnitude of the horizontal maximum principal stress, the magnitude of the horizontal minimum principal stress, the first angle, the second angle, and the third angle corresponding to the fault in the hydrate accumulation zone at different depths includes:

[0087] The effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths are calculated using the following formula;

[0088] σ n =σ v cos 2 α n +σ H cos 2 β n +σ h cos 2 γ n

[0089]

[0090] Among them, σ n is the effective normal stress, in MPa; τ is the shear stress, in MPa; σ v is the vertical principal stress, in MPa; σ H is the maximum horizontal principal stress, in MPa; σ h is the horizontal minimum principal stress, in MPa; α1 is the first angle, in °; α2 is the second angle, in °; α3 is the third angle, in °.

[0091] Specifically, the server uses the formula σ n =σ v cos 2 α n +σ H cos 2 β n +σ h cos 2 γ n Calculate the effective normal stress of the fault in the hydrate accumulation area at different depths using the formula Calculate the shear stress corresponding to the fault in the hydrate accumulation area at different depths.

[0092] In this way, according to the above formula, the accuracy of the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths can be improved.

[0093] Optionally, determining the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone according to the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths includes:

[0094] Obtain the friction coefficient corresponding to the fault in the hydrate accumulation area, and use the formula Determine the sliding pressure values ​​corresponding to the fault in the hydrate accumulation area at different depths; wherein, P s is the sliding pressure value, σ n is the effective normal stress; τ is the shear stress, μ f is the friction coefficient;

[0095] According to the sliding pressure values ​​corresponding to the fault in the hydrate accumulation zone at different depths, the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone is determined.

[0096] Specifically, the server obtains the friction coefficient corresponding to the fault in the hydrate accumulation area input by the user through the user device, and calculates it according to the formula Determine the sliding pressure values ​​corresponding to the faults in the hydrate accumulation area at different depths; where P s is the sliding pressure value, in MPa, σ n is the effective normal stress, in MPa; τ is the shear stress, in MPa, μ f is the friction coefficient.

[0097] The sliding pressure value is a function that changes with depth. The minimum value corresponding to the function is calculated, which is the minimum sliding pressure corresponding to the fault in the hydrate accumulation area.

[0098] In this way, the functional relationship between the sliding pressure value and the depth is first determined, and then the minimum value of the function is calculated to determine the minimum sliding pressure corresponding to the fault in the hydrate accumulation area, which can improve the efficiency and accuracy of determining the minimum sliding pressure.

[0099] Step 203: Determine the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation zone according to the minimum sliding pressure.

[0100] Specifically, the server determines the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation area based on the minimum sliding pressure, and finally determines the stability of the fault based on the size of the maximum tolerable gas column height. The larger the maximum tolerable gas column height, the stronger the stability of the fault, and the smaller the maximum tolerable gas column height, the weaker the stability of the fault.

[0101] Optionally, determining the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation zone according to the minimum sliding pressure includes:

[0102] The maximum tolerable gas column height corresponding to the fault in the hydrate accumulation area is calculated using the following formula:

[0103]

[0104] In the formula, h g is the maximum tolerable air column height; ρ w is the seawater density in the hydrate accumulation area; ρ g P is the gas density accumulated in the fault of the hydrate accumulation area; g is the minimum sliding pressure corresponding to the fault in the hydrate accumulation area.

[0105] Specifically, the functional relationship between the maximum tolerable air column height and the minimum sliding pressure is: In the formula, h g is the maximum tolerable air column height, in m; ρ w is the seawater density in the hydrate accumulation area, in g / cm 3 ρ g is the gas density accumulated in the fault of the hydrate accumulation area, in g / cm 3 ;P g is the minimum sliding pressure corresponding to the fault in the hydrate accumulation area, in MPa. The server can determine the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation area according to the functional relationship.

[0106] In this way, the accuracy of the maximum tolerable gas column height of the fault can be improved according to the above formula.

[0107] The method for determining the maximum tolerable gas column height of a fault provided in the present application can obtain the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths; the hydrate accumulation zone is a shallow hydrate accumulation zone under a tensile stress state; according to the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths, the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone is determined; wherein the minimum sliding pressure is used to indicate the minimum pressure value required to make the fault in the hydrate accumulation zone slide; according to the minimum sliding pressure, the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation zone is determined, which can improve the accuracy of determining the maximum tolerable gas column height of the fault, thereby improving the safety during hydrate exploitation and reducing the probability of geological disasters caused by formation instability.

[0108] The method of the present application is used to determine the maximum tolerable gas column height corresponding to the fault in the shallow hydrate accumulation zone under a certain tensile stress state.

[0109] The gas source of this hydrate accumulation area mainly comes from deep gas reservoirs, and the accumulation forms are hydrates, low-saturation natural gas and medium-to-high-saturation natural gas reservoirs. Relatively high-porosity and permeability sandstone reservoirs are developed in the area, which are the main accumulation layers of hydrates and natural gas. The sand bodies are widely distributed and continuous on the plane, providing good storage conditions for the accumulation of natural gas and hydrates. In addition, there are many normal faults in the area, which control the distribution of natural gas and hydrates.

[0110] Basic conditions for implementing this application: The study area has good three-dimensional seismic data and well logging interpretation data to provide comprehensive basic data.

[0111] Implementation process:

[0112] 1. Calculation of vertical principal stress: The water depth above the seabed in the hydrate accumulation area is 1700m. The density logging of 6 wells in the hydrate accumulation area is used to fit the regression curve of density with depth, so as to obtain the overlying formation pressure. Finally, the magnitude of the vertical principal stress is calculated according to the vertical principal stress calculation formula to obtain the vertical principal stress σ v The relationship with depth h is σ v =0.0175h-13.305.

[0113] 2. Calculate the maximum horizontal principal stress and the minimum horizontal principal stress: Substitute the rock mechanical parameters corresponding to the core in the fault of the hydrate accumulation area, such as the static Poisson's ratio μ = 0.27, the tectonic stress coefficient β1 is 0.475, the tectonic stress coefficient β2 is 0.235, and the vertical principal stress corresponding to the points at different depths into the Huang model to calculate the relationship between the maximum horizontal principal stress and the depth σ H =0.0161h-10.691, the relationship between the horizontal minimum principal stress and depth is:

[0114] σ h =0.0143h-7.4975.

[0115] 3. Obtain the friction coefficient of the rock in the hydrate accumulation area. Test the rock mechanical properties of the core in the fault of the hydrate accumulation area to obtain the friction coefficient μ f =0.5.

[0116] 4. Construct a three-dimensional fault model of the hydrate accumulation area through seismic interpretation, and determine the corresponding inclination of the fault at different depths based on the three-dimensional fault model, determine the direction of the unit surface normal based on the inclination, and finally determine the angles between the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress and the unit surface normal.

[0117] Figure 3 A schematic diagram of a three-dimensional model of a hydrate accumulation zone fault provided in an embodiment of the present application.

[0118] 5. According to the magnitude and three angles of vertical principal stress, maximum horizontal principal stress and minimum horizontal principal stress, use the formula to calculate the effective normal stress σ corresponding to different depths of the fault n , shear stress τ.

[0119] 6. Use formulas The minimum sliding pressure corresponding to the fault and the position corresponding to the minimum sliding pressure are determined, wherein the minimum sliding pressure corresponding to the fault is 0.29 MPa.

[0120] 7. Use formulas Calculate the maximum tolerable air column height corresponding to the fault, where the seawater density is 1g / cm 3 The average gas density of methane is 0.155g / cm 3 , g is the acceleration due to gravity, take 9.8m / s 2 , and the maximum tolerable air column height corresponding to the fault is 35.020m.

[0121] Table 1 is a schematic table of calculation data corresponding to each depth in a three-dimensional fault model provided in an embodiment of the present application, corresponding to Figure 3 Schematic diagram of the three-dimensional model.

[0122] Table 1

[0123]

[0124]

[0125] Corresponding to the above-mentioned method for determining the maximum tolerable air column height of a fault, an embodiment of the present application also provides a device for determining the maximum tolerable air column height of a fault. Figure 4 This is a schematic diagram of the structure of a device for determining the maximum tolerable gas column height of a fault provided in an embodiment of the present application. Figure 4 As shown, the device comprises:

[0126] An acquisition module 401 is used to obtain the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation area at different depths; the hydrate accumulation area is a shallow hydrate accumulation area under a tensile stress state;

[0127] The first determination module 402 is used to determine the minimum sliding pressure corresponding to the fault in the hydrate accumulation area according to the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation area at different depths; wherein the minimum sliding pressure is used to indicate the minimum pressure value required to make the fault in the hydrate accumulation area slide;

[0128] The second determination module 403 is used to determine the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation area according to the minimum sliding pressure.

[0129] Optionally, the vertical principal stress, the maximum horizontal principal stress, and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths all include magnitudes and directions. The first determination module 402 is specifically used to:

[0130] Acquire a three-dimensional model of the fault in the hydrate accumulation area, and determine the corresponding dip angles of the fault in the hydrate accumulation area at different depths according to the three-dimensional model;

[0131] Determine the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths according to the inclination angle, vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths;

[0132] According to the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths, the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone is determined.

[0133] Optionally, when determining the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths according to the dip angle, vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths, the first determining module 402 is specifically used to:

[0134] Determining the direction of the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths according to the inclination angle corresponding to the fault in the hydrate accumulation zone at different depths;

[0135] According to the direction of the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths, the direction of the vertical principal stress, the direction of the horizontal maximum principal stress and the direction of the horizontal minimum principal stress, determine the first angle between the vertical principal stress and the unit surface normal, the second angle between the horizontal maximum principal stress and the unit surface normal, and the third angle between the horizontal minimum principal stress and the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths;

[0136] According to the magnitude of the vertical principal stress, the magnitude of the horizontal maximum principal stress, the magnitude of the horizontal minimum principal stress, the first angle, the second angle and the third angle corresponding to the fault in the hydrate accumulation zone at different depths, the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths are determined.

[0137] Optionally, when the first determination module 402 determines the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths according to the magnitude of the vertical principal stress, the magnitude of the horizontal maximum principal stress, the magnitude of the horizontal minimum principal stress, the first angle, the second angle, and the third angle corresponding to the fault in the hydrate accumulation zone at different depths, it is specifically used to:

[0138] The effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths are calculated using the following formula;

[0139] σ n =σ v cos 2 α n +σ H cos 2 β n +σ h cos 2 γ n

[0140]

[0141] Among them, σ n is the effective normal stress; τ is the shear stress; σ1 is the vertical principal stress; σ2 is the maximum horizontal principal stress; σ3 is the minimum horizontal principal stress; α1 is the first angle; α2 is the second angle; α3 is the third angle.

[0142] Optionally, when determining the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone according to the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths, the first determining module 402 is specifically used to:

[0143] Obtain the friction coefficient corresponding to the fault in the hydrate accumulation area, and use the formula Determine the sliding pressure values ​​corresponding to the fault in the hydrate accumulation area at different depths; wherein, P s is the sliding pressure value, σ n is the effective normal stress; τ is the shear stress, μ f is the friction coefficient;

[0144] According to the sliding pressure values ​​corresponding to the fault in the hydrate accumulation zone at different depths, the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone is determined.

[0145] Optionally, the second determining module 403 is specifically configured to:

[0146] The maximum tolerable gas column height corresponding to the fault in the hydrate accumulation area is calculated using the following formula:

[0147]

[0148] In the formula, h g is the maximum tolerable air column height; ρ w is the seawater density in the hydrate accumulation area; ρ g P is the gas density accumulated in the fault of the hydrate accumulation area; g is the minimum sliding pressure corresponding to the fault in the hydrate accumulation area.

[0149] The specific implementation principle and effect of the device for determining the maximum tolerable gas column height of a fault provided in the embodiment of the present application can be found in the aforementioned embodiment and will not be repeated here.

[0150] Figure 5 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application. Figure 5 As shown, the electronic device of this embodiment may include:

[0151] at least one processor 501; and

[0152] A memory 502 communicatively connected to the at least one processor;

[0153] The memory 502 stores instructions that can be executed by the at least one processor 501, and the instructions are executed by the at least one processor 501 to enable the electronic device to execute the method described in any of the above embodiments.

[0154] Optionally, the memory 502 may be independent or integrated with the processor 501 .

[0155] The implementation principle and technical effects of the electronic device provided in this embodiment can be found in the aforementioned embodiments and will not be described in detail here.

[0156] An embodiment of the present application further provides a computer-readable storage medium, wherein the computer-readable storage medium stores computer-executable instructions. When a processor executes the computer-executable instructions, the method described in any of the above embodiments is implemented.

[0157] An embodiment of the present application further provides a computer program product, including a computer program, which implements the method described in any of the above embodiments when executed by a processor.

[0158] In the several embodiments provided in this application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic, for example, the division of the modules is only a logical function division, and there may be other division methods in actual implementation, such as multiple modules can be combined or integrated into another system, or some features can be ignored or not executed.

[0159] The above-mentioned integrated module implemented in the form of a software function module can be stored in a computer-readable storage medium. The above-mentioned software function module is stored in a storage medium, including a number of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) or a processor to perform some steps of the method described in each embodiment of the present application.

[0160] It should be understood that the above-mentioned processor can be a processing unit (Central Processing Unit, CPU), or other general-purpose processors, digital signal processors (Digital Signal Processor, DSP), application-specific integrated circuits (Application Specific Integrated Circuit, ASIC), etc. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor, etc. The steps of the method disclosed in the application can be directly embodied as being executed by a hardware processor, or executed by a combination of hardware and software modules in the processor. The memory may include a high-speed random access memory (Random Access Memory, RAM), and may also include non-volatile memory (Non-Volatile Memory, NVM), such as at least one disk storage, and can also be a USB flash drive, a mobile hard disk, a read-only memory, a disk or an optical disk, etc.

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

[0162] An exemplary storage medium is coupled to a processor so that the processor can read information from the storage medium and write information to the storage medium. Of course, the storage medium can also be a component of the processor. The processor and the storage medium can be located in an application specific integrated circuit (ASIC). Of course, the processor and the storage medium can also exist as discrete components in an electronic device or a main control device.

[0163] It should be noted that, in this article, the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the sentence "comprises a ..." does not exclude the existence of other identical elements in the process, method, article or device including the element.

[0164] The serial numbers of the above-mentioned embodiments of the present application are for description only and do not represent the advantages or disadvantages of the embodiments.

[0165] Through the description of the above implementation methods, those skilled in the art can clearly understand that the above-mentioned embodiment methods can be implemented by means of software plus a necessary general hardware platform, and of course by hardware, but in many cases the former is a better implementation method. Based on such an understanding, the technical solution of the present application, or the part that contributes to the prior art, can be embodied in the form of a software product, which is stored in a storage medium (such as ROM / RAM, magnetic disk, optical disk), and includes a number of instructions for a terminal device (which can be a mobile phone, computer, server, air conditioner, or network device, etc.) to execute the methods described in each embodiment of the present application.

[0166] The above are only preferred embodiments of the present application, and are not intended to limit the patent scope of the present application. Any equivalent structure or equivalent process transformation made using the contents of the present application specification and drawings, or directly or indirectly applied in other related technical fields, are also included in the patent protection scope of the present application.

Claims

1. A method for determining the maximum tolerable gas column height of a fault, characterized in that: include: Obtain the vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault in the hydrate accumulation area at different depths; The hydrate accumulation zone is a shallow hydrate accumulation zone under a tensile stress state; Determine the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone according to the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths; wherein the minimum sliding pressure is used to indicate the minimum pressure value required to make the fault in the hydrate accumulation zone slide; Determining the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation zone according to the minimum sliding pressure; Wherein, the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths all include magnitude and direction. According to the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation zone at different depths, the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone is determined, including: Acquire a three-dimensional model of the fault in the hydrate accumulation area, and determine the corresponding dip angles of the fault in the hydrate accumulation area at different depths according to the three-dimensional model; Determine the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths according to the inclination angle, vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths; Determine the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone according to the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths; Wherein, according to the dip angle, vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths, the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths are determined, including: Determining the direction of the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths according to the inclination angle corresponding to the fault in the hydrate accumulation zone at different depths; According to the direction of the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths, the direction of the vertical principal stress, the direction of the horizontal maximum principal stress and the direction of the horizontal minimum principal stress, determine the first angle between the vertical principal stress and the unit surface normal, the second angle between the horizontal maximum principal stress and the unit surface normal, and the third angle between the horizontal minimum principal stress and the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths; Determine the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths according to the magnitude of the vertical principal stress, the magnitude of the horizontal maximum principal stress, the magnitude of the horizontal minimum principal stress, the first angle, the second angle and the third angle corresponding to the fault in the hydrate accumulation zone at different depths; Among them, according to the magnitude of the vertical principal stress, the magnitude of the horizontal maximum principal stress, the magnitude of the horizontal minimum principal stress, the first angle, the second angle and the third angle corresponding to the fault in the hydrate accumulation zone at different depths, the effective normal stress and the shear stress corresponding to the fault in the hydrate accumulation zone at different depths are determined, including: The effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths are calculated using the following formula; s n =s v cos 2 a n +s H cos 2 b n +s h cos 2 c n Among them, σ n is the effective normal stress; τ is the shear stress; σ v is the vertical principal stress; σ H is the maximum horizontal principal stress; σ h is the horizontal minimum principal stress; α1 is the first angle; α2 is the second angle; α3 is the third angle; Wherein, according to the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths, determining the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone includes: Obtain the friction coefficient corresponding to the fault in the hydrate accumulation area, and use the formula Determine the sliding pressure values ​​corresponding to the fault in the hydrate accumulation area at different depths; wherein, P s is the sliding pressure value, σ n is the effective normal stress; τ is the shear stress, μ f is the friction coefficient; Determining the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone according to the sliding pressure values ​​corresponding to the fault in the hydrate accumulation zone at different depths; Wherein, determining the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation zone according to the minimum sliding pressure includes: The maximum tolerable gas column height corresponding to the fault in the hydrate accumulation area is calculated using the following formula: In the formula, h g is the maximum tolerable air column height; ρ w is the seawater density in the hydrate accumulation area; ρ g P is the gas density accumulated in the fault of the hydrate accumulation area; g is the minimum sliding pressure corresponding to the fault in the hydrate accumulation area.

2. A device for determining the maximum tolerable air column height of a fault, characterized in that: include: An acquisition module is used to obtain the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation area at different depths; The hydrate accumulation zone is a shallow hydrate accumulation zone under a tensile stress state; A first determination module is used to determine the minimum sliding pressure corresponding to the fault in the hydrate accumulation area according to the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation area at different depths; wherein the minimum sliding pressure is used to indicate the minimum pressure value required to make the fault in the hydrate accumulation area slide; A second determination module is used to determine the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation area according to the minimum sliding pressure; Among them, the vertical principal stress, the maximum horizontal principal stress and the minimum horizontal principal stress corresponding to the fault in the hydrate accumulation area at different depths all include magnitude and direction. The first determination module is specifically used for: Acquire a three-dimensional model of the fault in the hydrate accumulation area, and determine the corresponding dip angles of the fault in the hydrate accumulation area at different depths according to the three-dimensional model; Determine the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths according to the inclination angle, vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault in the hydrate accumulation zone at different depths; Determine the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone according to the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths; Among them, the first determination module, when determining the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation area at different depths according to the dip angle, vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault in the hydrate accumulation area at different depths, is specifically used to: Determining the direction of the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths according to the inclination angle corresponding to the fault in the hydrate accumulation zone at different depths; According to the direction of the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths, the direction of the vertical principal stress, the direction of the horizontal maximum principal stress and the direction of the horizontal minimum principal stress, determine the first angle between the vertical principal stress and the unit surface normal, the second angle between the horizontal maximum principal stress and the unit surface normal, and the third angle between the horizontal minimum principal stress and the unit surface normal corresponding to the fault in the hydrate accumulation zone at different depths; Determine the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths according to the magnitude of the vertical principal stress, the magnitude of the horizontal maximum principal stress, the magnitude of the horizontal minimum principal stress, the first angle, the second angle and the third angle corresponding to the fault in the hydrate accumulation zone at different depths; Among them, the first determination module is specifically used to determine the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation area at different depths according to the magnitude of the vertical principal stress, the magnitude of the horizontal maximum principal stress, the magnitude of the horizontal minimum principal stress, the first angle, the second angle and the third angle corresponding to the fault in the hydrate accumulation area at different depths: The effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths are calculated using the following formula; s n =s v cos 2 a n +s H cos 2 b n +s h cos 2 c n Among them, σ n is the effective normal stress; τ is the shear stress; σ v is the vertical principal stress; σ H is the maximum horizontal principal stress; σ h is the horizontal minimum principal stress; α1 is the first angle; α2 is the second angle; α3 is the third angle; Wherein, when the first determination module determines the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone according to the effective normal stress and shear stress corresponding to the fault in the hydrate accumulation zone at different depths, it is specifically used to: Obtain the friction coefficient corresponding to the fault in the hydrate accumulation area, and use the formula Determine the sliding pressure values ​​corresponding to the fault in the hydrate accumulation area at different depths; wherein, P s is the sliding pressure value, σ n is the effective normal stress; τ is the shear stress, μ f is the friction coefficient; Determining the minimum sliding pressure corresponding to the fault in the hydrate accumulation zone according to the sliding pressure values ​​corresponding to the fault in the hydrate accumulation zone at different depths; Wherein, when the second determination module determines the maximum tolerable gas column height corresponding to the fault in the hydrate accumulation zone according to the minimum sliding pressure, it is specifically used to: The maximum tolerable gas column height corresponding to the fault in the hydrate accumulation area is calculated using the following formula: In the formula, h g is the maximum tolerable air column height; ρ w is the seawater density in the hydrate accumulation area; ρ g P is the gas density accumulated in the fault of the hydrate accumulation area; g is the minimum sliding pressure corresponding to the fault in the hydrate accumulation area.

3. An electronic device, characterized in that: include: A processor, and a memory communicatively connected to the processor; The memory stores computer-executable instructions; The processor executes the computer-executable instructions stored in the memory to implement the method according to claim 1.

4. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer-executable instructions, which are used to implement the method according to claim 1 when executed by a processor.

5. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the method according to claim 1 is implemented.

Images