[go: up one dir, main page]

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

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

Info

Publication number
CN118653803B
CN118653803B CN202410710377.XA CN202410710377A CN118653803B CN 118653803 B CN118653803 B CN 118653803B CN 202410710377 A CN202410710377 A CN 202410710377A CN 118653803 B CN118653803 B CN 118653803B
Authority
CN
China
Prior art keywords
fault
stress
hydrate
main stress
different depths
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202410710377.XA
Other languages
Chinese (zh)
Other versions
CN118653803A (en
Inventor
任金锋
李廷微
徐梦婕
王笑雪
何玉林
梁劲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangzhou Marine Geological Survey
Original Assignee
Guangzhou Marine Geological Survey
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangzhou Marine Geological Survey filed Critical Guangzhou Marine Geological Survey
Priority to CN202410710377.XA priority Critical patent/CN118653803B/en
Publication of CN118653803A publication Critical patent/CN118653803A/en
Application granted granted Critical
Publication of CN118653803B publication Critical patent/CN118653803B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/01—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells specially adapted for obtaining from underwater installations
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells
    • E21B47/001—Survey of boreholes or wells for underwater installation

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

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

Method, device, equipment and medium for determining maximum bearable gas column height of fault
Technical Field
The application relates to the technical field of geological exploration and development of oil and gas reservoirs, in particular to a method, a device, equipment and a medium for determining the height of a fault maximum bearable gas column.
Background
Hydrates are often present in land permanent frozen soil layers 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 faults of hydrate accumulation areas, and once the height of the accumulated gas column in the faults exceeds the bearable range of the faults, the faults slide, possibly causing stratum instability and inducing geological disasters such as landslide, ground subsidence and even submarine debris flow. This not only poses a threat to the mining activity, but may also have a serious impact on the surrounding environment.
The prior art has the problem of low accuracy when determining the maximum bearable gas column height of faults.
Therefore, a method for determining the maximum sustainable air column height of a fault is needed to improve the accuracy of determining the maximum sustainable air column height of the fault.
Disclosure of Invention
The application provides a method, a device, equipment and a medium for determining the maximum bearable gas column height of a fault, which are used for improving the accuracy of determining the maximum bearable gas column height of the fault.
In a first aspect, an embodiment of the present application provides a method for determining a maximum sustainable air column height of a fault, including:
The method comprises the steps of obtaining vertical main stress, horizontal maximum main stress and horizontal minimum main stress corresponding to faults of a hydrate aggregation area at different depths, wherein the hydrate aggregation area is a shallow hydrate aggregation area in an extension stress state;
Determining a sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the hydrate aggregation zone fault at different depths, wherein the sliding pressure minimum value is used for indicating a minimum pressure value required for enabling the hydrate aggregation zone fault to slide;
And determining the maximum bearable gas column height corresponding to the fault of the hydrate aggregation area according to the minimum value of the sliding pressure.
Optionally, the vertical principal stress, the horizontal maximum principal stress and the horizontal minimum principal stress corresponding to the fault of the hydrate aggregation area at different depths all include a magnitude and a direction, and determining the sliding pressure minimum corresponding to the fault of the hydrate aggregation area according to the vertical principal stress, the horizontal maximum principal stress and the horizontal minimum principal stress corresponding to the fault of the hydrate aggregation area at different depths includes:
Acquiring a three-dimensional model of the fault of the hydrate aggregation area, and determining corresponding dip angles of the fault of the hydrate aggregation area at different depths according to the three-dimensional model;
determining effective normal stress and shear stress corresponding to the fault of the hydrate aggregation area at different depths according to inclination angles, vertical main stress, horizontal maximum main stress and horizontal minimum main stress corresponding to the fault of the hydrate aggregation area at different depths;
And determining a sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to the effective normal stress and the shearing stress corresponding to the hydrate aggregation zone fault at different depths.
Optionally, determining the effective normal stress and the shear stress corresponding to the hydrate aggregation zone fault at different depths according to the inclination angle, the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the hydrate aggregation zone fault at different depths includes:
determining the direction of the normal line of the unit surface corresponding to the fault of the hydrate aggregation area at different depths according to the corresponding inclination angles of the fault of the hydrate aggregation area at different depths;
Determining a first included angle between the vertical main stress and the normal line of the unit surface, a second included angle between the horizontal maximum main stress and the normal line of the unit surface and a third included angle between the horizontal minimum main stress and the normal line of the unit surface, which correspond to the hydrate aggregation area fault at different depths, according to the direction of the normal line of the unit surface, the direction of the vertical main stress, the direction of the horizontal maximum main stress and the direction of the horizontal minimum main stress, which correspond to the hydrate aggregation area fault at different depths;
and determining effective positive stress and shear stress corresponding to the fault of the hydrate aggregation area at different depths according to the vertical main stress, the horizontal maximum main stress, the horizontal minimum main stress, the first included angle, the second included angle and the third included angle corresponding to the fault of the hydrate aggregation area at different depths.
Optionally, determining, according to the magnitude of the vertical main stress, the magnitude of the horizontal maximum main stress, the magnitude of the horizontal minimum main stress, the first included angle, the second included angle and the third included angle corresponding to the fault of the hydrate aggregation area at different depths, the effective normal stress and the shear stress corresponding to the fault of the hydrate aggregation area at different depths includes:
Calculating effective positive stress and shear stress corresponding to the fault of the hydrate aggregation zone at different depths by using the following formula;
σn＝σvcos2αn+σHcos2βn+σhcos2γn
Wherein σ n is effective normal stress, τ is shear stress, σ v is vertical main stress, σ H is horizontal maximum main stress, σ h is horizontal minimum main stress, α 1 is a first included angle, α 2 is a second included angle, and α 3 is a third included angle.
Optionally, determining the sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to the effective normal stress and the shear stress corresponding to the hydrate aggregation zone fault at different depths includes:
obtaining the friction coefficient corresponding to the fault of the hydrate aggregation area, and according to a formula Determining sliding pressure values corresponding to the hydrate aggregation zone faults at different depths, wherein P s is the sliding pressure value, sigma n is effective positive stress, tau is shear stress, mu f is friction coefficient;
And determining the minimum sliding pressure corresponding to the fault of the hydrate aggregation zone according to the sliding pressure values corresponding to the fault of the hydrate aggregation zone at different depths.
Optionally, determining a maximum bearable gas column height corresponding to the fault of the hydrate aggregation area according to the sliding pressure minimum value includes:
Calculating the maximum bearable gas column height corresponding to the fault of the hydrate aggregation zone by using the following formula:
Wherein h g is the maximum bearable gas column height, ρ w is the sea water density of the hydrate aggregation area, ρ g is the gas density aggregated by the hydrate aggregation area fault, and P g is the sliding pressure minimum value corresponding to the hydrate aggregation area fault.
In a second aspect, an embodiment of the present application provides a fault maximum sustainable air column height determining apparatus, including:
the acquisition module is used for acquiring vertical main stress, horizontal maximum main stress and horizontal minimum main stress corresponding to faults of the hydrate aggregation area at different depths, wherein the hydrate aggregation area is a shallow hydrate aggregation area in an extension stress state;
The first determining module is used for determining a sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the hydrate aggregation zone fault at different depths, wherein the sliding pressure minimum value is used for indicating a minimum pressure value required for enabling the hydrate aggregation zone fault to slide;
And the second determining module is used for determining the maximum bearable gas column height corresponding to the fault of the hydrate aggregation area according to the minimum value of the sliding pressure.
In a third aspect, an embodiment of the present application provides an electronic device, including a processor, and a memory communicatively coupled to the processor;
the memory stores computer-executable instructions;
the processor executes computer-executable instructions stored in the memory to implement the method as described in the first aspect.
In a fourth aspect, embodiments of the present application provide a computer-readable storage medium having stored therein computer-executable instructions for performing the method of the first aspect when executed by a processor.
In a fifth aspect, embodiments of the present application provide a computer program product comprising a computer program which, when executed by a processor, implements a method as described in the first aspect.
The method comprises the steps of obtaining vertical main stress, horizontal maximum main stress and horizontal minimum main stress corresponding to faults of a hydrate aggregation area at different depths, enabling the hydrate aggregation area to be a shallow hydrate aggregation area in an extension stress state, determining sliding pressure minimum values corresponding to the faults of the hydrate aggregation area according to the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the faults of the hydrate aggregation area at different depths, wherein the sliding pressure minimum values are used for indicating minimum pressure values required for enabling the faults of the hydrate aggregation area to slide, determining the maximum bearable air column height corresponding to the faults of the hydrate aggregation area according to the sliding pressure minimum values, improving accuracy of determining the maximum bearable air column height of the faults, improving safety in a hydrate exploitation process and reducing geological disasters caused by stratum instability.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the application and together with the description, serve to explain the principles of the application.
FIG. 1 is an application scenario diagram provided in an embodiment of the present application;
FIG. 2 is a schematic flow chart of a method for determining a maximum sustainable air column height of a fault according to an embodiment of the present application;
FIG. 3 is a schematic diagram of a three-dimensional model of a fault of a hydrate aggregation area according to an embodiment of the present application;
FIG. 4 is a schematic structural diagram of a fault maximum sustainable air column height determining device according to an embodiment of the present application;
fig. 5 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Specific embodiments of the present application have been shown by way of the above drawings and will be described in more detail below. The drawings and the written description are not intended to limit the scope of the inventive concepts in any way, but rather to illustrate the inventive concepts to those skilled in the art by reference to the specific embodiments.
Detailed Description
Reference will now be made in detail to exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, the same numbers in different drawings refer to the same or similar elements, unless otherwise indicated. The implementations described in the following exemplary examples do not represent all implementations consistent with the application. Rather, they are merely examples of apparatus and methods consistent with aspects of the application as detailed in the accompanying claims.
Hydrates are often present in land permanent frozen soil layers 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 faults of hydrate accumulation areas, and once the height of the accumulated gas column in the faults exceeds the bearable range of the faults, the faults slide, possibly causing stratum instability and inducing geological disasters such as landslide, ground subsidence and even submarine debris flow. This not only poses a threat to the mining activity, but may also have a serious impact on the surrounding environment.
In addition, in the geological sealing process of the carbon dioxide in the sea area, if the gas injection pressure is too high and exceeds the pressure-bearing upper limit of the fault, the risk of carbon dioxide leakage can be greatly increased.
The prior art has the problem of low accuracy when determining the maximum bearable gas column height of faults.
In view of the above, the application provides a method for determining the maximum bearable gas column height of a fault, which is used for determining the maximum bearable gas column height corresponding to the fault of a shallow hydrate aggregation zone, firstly, obtaining the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the fault of the hydrate aggregation zone at different depths, determining the sliding pressure minimum corresponding to the fault of the hydrate aggregation zone according to the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the fault of the hydrate aggregation zone at different depths, and finally determining the maximum bearable gas column height corresponding to the fault of the hydrate aggregation zone according to the sliding pressure minimum. According to the method, the sliding pressure minimum value corresponding to the fault is determined according to the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the fault at different depths, and the maximum bearable gas column height corresponding to the fault in the hydrate aggregation area is determined according to the sliding pressure minimum value, so that the accuracy of determining the maximum bearable gas column height of the fault can be improved, the safety in the hydrate exploitation process is further improved, and the probability of geological disasters caused by stratum instability is reduced.
Fig. 1 is an application scenario diagram provided by the embodiment of the present application, as shown in fig. 1, a user inputs vertical main stress, horizontal maximum main stress and horizontal minimum main stress corresponding to a fault at different depths through a user device, the input form may be a functional form, and by way of example, the user may input a function of vertical main stress changing with depth through the user device, the function of horizontal maximum main stress changing with depth and the function of horizontal minimum main stress changing with depth are determined, after the input is completed, click a determination button, the user device sends the function of vertical main stress changing with depth, the function of horizontal maximum main stress changing with depth and the function of horizontal minimum main stress changing with depth to a server, after the server receives the function of vertical main stress changing with depth, the function of horizontal maximum main stress changing with depth and the function of horizontal minimum main stress changing with depth, and determines a sliding pressure minimum value corresponding to a hydrate aggregation area fault, and sends a corresponding maximum bearable air column height corresponding to the hydrate aggregation area to the user device, so that the maximum bearable air column height can be displayed on the user device.
The following describes the technical scheme of the present application and how the technical scheme of the present application solves the above technical problems in detail with specific embodiments. The following embodiments may be combined with each other, and the same or similar concepts or processes may not be described in detail in some embodiments. Embodiments of the present application will be described below with reference to the accompanying drawings.
Fig. 2 is a schematic flow chart of a method for determining a maximum sustainable air column height of a fault according to an embodiment of the present application. The execution body of the embodiment may be any device having a data processing function, and the present application is specifically described with reference to a server as an execution body. As shown in fig. 2, a method for determining a maximum sustainable air column height of a fault provided by an embodiment of the present application may include:
Step 201, vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to faults of a hydrate aggregation area at different depths are obtained, and the hydrate aggregation area is a shallow hydrate aggregation area in an extension stress state.
Specifically, the vertical principal stress, the horizontal maximum principal stress and the horizontal minimum principal stress are respectively compressive stresses to which the rock is subjected in three-dimensional space. In the present application, the fracture of the hydrate aggregation zone is subjected to compressive stress in three-dimensional space.
The direction of the vertical main stress is along the vertical direction, and the vertical main stress is the main stress with the largest value among the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress in the stretching stress state.
The horizontal maximum principal stress and the horizontal minimum principal stress are along the horizontal direction, and the horizontal maximum principal stress and the horizontal minimum principal stress are mutually perpendicular. The horizontal maximum principal stress is the principal stress with the numerical value centered in the vertical principal stress, the horizontal maximum principal stress and the horizontal minimum principal stress in the stretching stress state, and the horizontal minimum principal stress is the principal stress with the numerical value minimum in the vertical principal stress, the horizontal maximum principal stress and the horizontal minimum principal stress.
Wherein for a point of depth in the fault of the hydrate accumulation zone, the vertical principal stress to which the point is subjected is formed by the pressure of sea water above the sea floor together with the gravity of all formations above the point.
The functional relation of the vertical main stress corresponding to the fault of the hydrate aggregation area along with the change of the depth is as follows:
The method comprises the steps of obtaining a hydrate aggregation area, wherein sigma v is vertical main stress, the unit is MPa, h is a depth value corresponding to the depth to be calculated of the hydrate aggregation area fault, the unit is m, rho (h) is a function of the change of the target stratum density along with the depth, the function of the change of the target stratum density along with the depth is obtained by density logging, the target stratum comprises the hydrate aggregation area fault and a stratum above the fault, g is gravity acceleration, 9.80m/s 2;ρw can be taken as the sea water density above the hydrate aggregation area, the unit is g/cm 3;h0, and the unit is m.
And aiming at the fault of the hydrate aggregation area, taking the rock core in the fault area to perform rock mechanical property test to obtain rock mechanical parameters, and substituting the rock mechanical parameters and vertical main stress into a Huang model to obtain horizontal minimum main stress and horizontal maximum main stress.
Wherein sigma v is vertical main stress, the unit is MPa, P w is pore fluid pressure, the unit is MPa, mu is static Poisson ratio, alpha is Biot coefficient, beta 1 and beta 2 are structural stress coefficients, and the static Poisson ratio is an important parameter in material mechanics and is used for describing the relation between transverse strain and longitudinal strain of a material under a stress state. The Biot coefficient is a parameter used in rock mechanics and geomechanics to describe the relationship between effective stress and total stress of a pore medium (e.g., rock or soil) under external stress. The structural stress coefficient is an important parameter used for describing the stress state of the crust in geomechanics and rock mechanics, reflects the characteristics of the structural stress field of the crust, the static poisson ratio and the Biot coefficient are obtained through rock mechanical property tests of core samples in faults of a hydrate aggregation area, and the structural stress coefficients beta 1 and beta 2 can be empirical values selected according to the actual conditions of the hydrate aggregation area.
The rock mechanical property test refers to the test and analysis of mechanical properties of a rock core so as to know the mechanical behaviors and characteristics of the rock under different stress conditions. The rock mechanical parameters can include compressive strength, tensile strength, shear strength, elastic modulus, plastic modulus, poisson's ratio, friction angle, failure criterion parameters, and the like.
Specifically, a user can input vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to the fault of the hydrate aggregation area at different depths through the user equipment, the input form can be a function form, and the user can input a function of vertical principal stress changing along with the depth, a function of horizontal maximum principal stress changing along with the depth and a function of horizontal minimum principal stress changing along with the depth on an input interface of the user equipment, after the input is completed, click a confirmation button, the user equipment sends the function of vertical principal stress changing along with the depth, the function of horizontal maximum principal stress changing along with the depth, the function of horizontal minimum principal stress changing along with the depth to a server, and the server obtains the function of vertical principal stress changing along with the depth, the function of horizontal maximum principal stress changing along with the depth and the function of horizontal minimum principal stress changing along with the depth.
And 202, determining a sliding pressure minimum value corresponding to the fault of the hydrate aggregation zone according to the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the fault of the hydrate aggregation zone at different depths, wherein the sliding pressure minimum value is used for indicating a minimum pressure value required for enabling the fault of the hydrate aggregation zone to slide.
Specifically, the server determines a sliding pressure minimum value corresponding to the fault of the hydrate aggregation area according to the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the fault of the hydrate aggregation area at different depths, namely a minimum pressure value required for enabling the fault of the hydrate aggregation area to slide.
Optionally, the vertical principal stress, the horizontal maximum principal stress and the horizontal minimum principal stress corresponding to the fault of the hydrate aggregation area at different depths all include a magnitude and a direction, and determining the sliding pressure minimum corresponding to the fault of the hydrate aggregation area according to the vertical principal stress, the horizontal maximum principal stress and the horizontal minimum principal stress corresponding to the fault of the hydrate aggregation area at different depths includes:
Acquiring a three-dimensional model of the fault of the hydrate aggregation area, and determining corresponding dip angles of the fault of the hydrate aggregation area at different depths according to the three-dimensional model;
determining effective normal stress and shear stress corresponding to the fault of the hydrate aggregation area at different depths according to inclination angles, vertical main stress, horizontal maximum main stress and horizontal minimum main stress corresponding to the fault of the hydrate aggregation area at different depths;
And determining a sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to the effective normal stress and the shearing stress corresponding to the hydrate aggregation zone fault at different depths.
The three-dimensional model of the fault of the hydrate accumulation zone can be obtained by explaining the seismic data corresponding to the fault of the hydrate accumulation zone.
For each point on the fault of the hydrate accumulation zone, the corresponding tangent plane of the point is tangent to the fault, the effective normal stress corresponding to the point is the normal stress perpendicular to the corresponding tangent plane of the point, and the shear stress corresponding to the point is the stress parallel to the corresponding tangent plane of the point.
The method comprises the steps of obtaining a three-dimensional model of a fault of a hydrate aggregation area, determining corresponding inclination angles of the fault of the hydrate aggregation area at different depths according to the three-dimensional model, determining corresponding effective normal stress and shear stress of the fault of the hydrate aggregation area at different depths according to the corresponding inclination angles of the fault of the hydrate aggregation area at different depths, the magnitude and direction of vertical main stress, the magnitude and direction of horizontal maximum main stress and the magnitude and direction of horizontal minimum main stress, and finally determining corresponding sliding pressure minimum value of the fault of the hydrate aggregation area according to the corresponding effective normal stress and shear stress of the fault of the hydrate aggregation area at different depths.
Therefore, according to the inclination angle, the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the faults of the hydrate aggregation area at different depths, the effective normal stress and the shearing stress corresponding to the faults of the hydrate aggregation area at different depths are determined, the accuracy of the determined effective normal stress and shearing stress can be improved, and the accuracy of the sliding pressure minimum value corresponding to the determined faults is further improved.
Optionally, determining the effective normal stress and the shear stress corresponding to the hydrate aggregation zone fault at different depths according to the inclination angle, the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the hydrate aggregation zone fault at different depths includes:
determining the direction of the normal line of the unit surface corresponding to the fault of the hydrate aggregation area at different depths according to the corresponding inclination angles of the fault of the hydrate aggregation area at different depths;
Determining a first included angle between the vertical main stress and the normal line of the unit surface, a second included angle between the horizontal maximum main stress and the normal line of the unit surface and a third included angle between the horizontal minimum main stress and the normal line of the unit surface, which correspond to the hydrate aggregation area fault at different depths, according to the direction of the normal line of the unit surface, the direction of the vertical main stress, the direction of the horizontal maximum main stress and the direction of the horizontal minimum main stress, which correspond to the hydrate aggregation area fault at different depths;
and determining effective positive stress and shear stress corresponding to the fault of the hydrate aggregation area at different depths according to the vertical main stress, the horizontal maximum main stress, the horizontal minimum main stress, the first included angle, the second included angle and the third included angle corresponding to the fault of the hydrate aggregation area at different depths.
For each point of the hydrate aggregation zone fault at each depth, the tangent plane of the fault at the point is a unit surface, and the normal line of the unit surface is perpendicular to the unit surface, that is to say, the inclination angle is added by 90 degrees, namely, the included angle between the normal line of the unit surface and the horizontal plane.
Specifically, the server determines the direction of the normal line of the unit surface corresponding to the fault of the hydrate aggregation area at different depths according to the corresponding inclination angles of the fault of the hydrate aggregation area at different depths, then determines the included angle between the vertical main stress corresponding to the fault of the hydrate aggregation area at different depths and the normal line of the unit surface, namely a first included angle, determines the included angle between the horizontal maximum main stress corresponding to the fault of the hydrate aggregation area at different depths and the normal line of the unit surface, namely a second included angle, and determines the included angle between the horizontal minimum main stress corresponding to the fault of the hydrate aggregation area at different depths and the normal line of the unit surface, namely a third included angle.
And finally, determining effective normal stress and shear stress corresponding to the fault of the hydrate aggregation area at different depths according to the vertical main stress, the horizontal maximum main stress, the horizontal minimum main stress, the first included angle, the second included angle and the third included angle corresponding to the fault of the hydrate aggregation area at different depths.
In this way, according to the inclination angles corresponding to faults at different depths, the direction of the normal line of the unit surface is determined, then three included angles are determined according to the direction of the normal line of the unit surface, and finally the effective normal stress and the shear stress are determined according to the magnitudes of the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress and the three included angles, so that the accuracy of the determined effective normal stress and the determined shear stress can be improved.
Optionally, determining, according to the magnitude of the vertical main stress, the magnitude of the horizontal maximum main stress, the magnitude of the horizontal minimum main stress, the first included angle, the second included angle and the third included angle corresponding to the fault of the hydrate aggregation area at different depths, the effective normal stress and the shear stress corresponding to the fault of the hydrate aggregation area at different depths includes:
Calculating effective positive stress and shear stress corresponding to the fault of the hydrate aggregation zone at different depths by using the following formula;
σn＝σvcos2αn+σHcos2βn+σhcos2γn
wherein, sigma n is effective normal stress, tau is shear stress, sigma v is vertical main stress, sigma H is horizontal maximum main stress, sigma h is horizontal minimum main stress, alpha 1 is a first included angle, alpha 2 is a second included angle, alpha 3 is a third included angle, and alpha is a third included angle.
Specifically, the server calculates the effective normal stress corresponding to the fault of the hydrate aggregation area at different depths by using a formula σn＝σvcos2αn+σHcos2βn+σhcos2γn, and uses the formulaThe shear stress corresponding to the hydrate accumulation zone fault at different depths was calculated.
Thus, the accuracy of effective normal stress and shear stress corresponding to the determined fault of the hydrate aggregation zone at different depths can be improved according to the formula.
Optionally, determining the sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to the effective normal stress and the shear stress corresponding to the hydrate aggregation zone fault at different depths includes:
obtaining the friction coefficient corresponding to the fault of the hydrate aggregation area, and according to a formula Determining sliding pressure values corresponding to the hydrate aggregation zone faults at different depths, wherein P s is the sliding pressure value, sigma n is effective positive stress, tau is shear stress, mu f is friction coefficient;
And determining the minimum sliding pressure corresponding to the fault of the hydrate aggregation zone according to the sliding pressure values corresponding to the fault of the hydrate aggregation zone at different depths.
Specifically, the server obtains the friction coefficient corresponding to the fault of the hydrate aggregation area input by the user through the user equipment, and according to the formulaAnd determining sliding pressure values corresponding to faults of the hydrate aggregation area at different depths, wherein P s is the sliding pressure value, the unit is MPa, sigma n is effective positive stress, the unit is MPa, tau is shear stress, the unit is MPa, and mu f is friction coefficient.
The sliding pressure value is a function changing along with the depth, and the minimum value corresponding to the function is calculated, namely the sliding pressure minimum value corresponding to the fault of the hydrate aggregation area.
Therefore, the functional relation of the sliding pressure value along with the change of the depth is determined, and then the minimum value of the function is calculated to determine the minimum value of the sliding pressure corresponding to the fault of the hydrate aggregation area, so that the efficiency and the accuracy of determining the minimum value of the sliding pressure can be improved.
And 203, determining the maximum bearable gas column height corresponding to the fault of the hydrate aggregation area according to the minimum value of the sliding pressure.
Specifically, the server determines the maximum bearable air column height corresponding to the fault of the hydrate aggregation area according to the minimum value of the sliding pressure, and finally determines the fault stability according to the maximum bearable air column height, wherein the greater the maximum bearable air column height is, the stronger the fault stability is, the smaller the maximum bearable air column height is, and the weaker the fault stability is.
Optionally, determining a maximum bearable gas column height corresponding to the fault of the hydrate aggregation area according to the sliding pressure minimum value includes:
Calculating the maximum bearable gas column height corresponding to the fault of the hydrate aggregation zone by using the following formula:
Wherein h g is the maximum bearable gas column height, ρ w is the sea water density of the hydrate aggregation area, ρ g is the gas density aggregated by the hydrate aggregation area fault, and P g is the sliding pressure minimum value corresponding to the hydrate aggregation area fault.
Specifically, the maximum sustainable air column height as a function of sliding pressure minimum is Wherein h g is the maximum acceptable air column height, the unit is m, ρ w is the sea water density of the hydrate aggregation area, the unit is g/cm 3;ρg is the gas density aggregated by the fault of the hydrate aggregation area, the unit is g/cm 3;Pg is the sliding pressure minimum value corresponding to the fault of the hydrate aggregation area, and the unit is MPa. The server can determine the maximum bearable gas column height corresponding to the fault of the hydrate aggregation area according to the functional relation.
Thus, the accuracy of the determined maximum sustainable air column height of the fault can be improved according to the formula.
The method for determining the maximum bearable gas column height of the fault provided by the application can obtain the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the fault of the hydrate aggregation area at different depths, wherein the hydrate aggregation area is a shallow hydrate aggregation area in an extension stress state, the sliding pressure minimum value corresponding to the fault of the hydrate aggregation area is determined according to the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the fault of the hydrate aggregation area at different depths, the sliding pressure minimum value is used for indicating the minimum pressure value required for enabling the fault of the hydrate aggregation area to slide, the maximum bearable gas column height corresponding to the fault of the hydrate aggregation area is determined according to the sliding pressure minimum value, the accuracy of determining the maximum bearable gas column height of the fault can be improved, the safety of the hydrate in the mining process is further improved, and the probability of geological disasters caused by stratum instability is reduced.
The method of the application is adopted to determine the maximum bearable gas column height corresponding to the fault of the shallow hydrate aggregation area under a certain stretching stress state.
The gas source of the hydrate accumulation zone mainly comes from a deep gas reservoir, and the accumulation forms are hydrate, low-saturation natural gas and medium-high-saturation natural gas reservoirs. The relatively high pore sandstone reservoir developed in the area is a main aggregation layer system of the hydrate and the natural gas, and the sand on the plane is widely and continuously distributed, so that good reservoir conditions are provided for the aggregation of the natural gas and the hydrate. In addition, a plurality of positive faults develop in the region, and the positive faults control the distribution rule of natural gas and hydrate.
The application has the basic condition that the research area has better three-dimensional seismic data and logging interpretation data so as to provide comprehensive basic data.
The implementation process comprises the following steps:
1. Calculating vertical main stress, namely, the depth of a water body above the sea bottom of a hydrate gathering area is 1700m, fitting a regression curve of density along with the depth by using the density log of 6 wells of the hydrate gathering area, thereby obtaining the pressure of an overlying stratum, and finally, calculating the size of the vertical main stress according to a vertical main stress calculation formula to obtain the relation between the vertical main stress sigma v and the depth h as sigma v =0.0175 h-13.305.
2. Calculating the horizontal maximum principal stress and the horizontal minimum principal stress by substituting rock mechanical parameters corresponding to a rock core in a fault of a hydrate aggregation zone, such as static poisson ratio mu=0.27, a structural stress coefficient beta 1 of 0.475 and a structural stress coefficient beta 2 of 0.235 into a yellow model, wherein the change relation of the horizontal maximum principal stress with depth is sigma H = 0.0161h-10.691, and the change relation of the horizontal minimum principal stress with depth is
σh=0.0143h-7.4975。
3. The coefficient of friction of the hydrate aggregation zone rock is obtained. Rock mechanical property test is carried out on the rock core in the fault of the hydrate aggregation zone to obtain the friction coefficient mu f =0.5.
4. And constructing a fault three-dimensional model of the hydrate aggregation area through seismic interpretation, determining corresponding dip angles of faults at different depths according to the fault three-dimensional model, determining the direction of the normal line of the unit surface according to the dip angles, and finally determining the included angles of the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress with the normal line of the unit surface respectively.
Fig. 3 is a schematic diagram of a fault three-dimensional model of a hydrate aggregation area according to an embodiment of the present application.
5. And calculating the effective normal stress sigma n and the shearing stress tau corresponding to different depths of the fault by using a formula according to the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress and three included angles.
6. Using the formulaAnd determining the sliding pressure minimum value corresponding to the fault and the position corresponding to the sliding pressure minimum value, wherein the sliding pressure minimum value corresponding to the fault is 0.29MPa.
7. Using the formulaCalculating the maximum bearable gas column height corresponding to the fault, wherein the sea water density is 1g/cm 3, the methane average gas density is 0.155g/cm 3, g is the gravity acceleration, 9.8m/s 2, and the maximum bearable gas column height corresponding to the fault is 35.020m.
Table 1 is a schematic representation of calculated data corresponding to each depth in a fault three-dimensional model according to an embodiment of the present application, and corresponds to the three-dimensional model in fig. 3.
TABLE 1
Corresponding to the method for determining the maximum sustainable air column height of the fault, the embodiment of the application also provides a device for determining the maximum sustainable air column height of the fault. Fig. 4 is a schematic structural diagram of a fault maximum bearable air column height determining device according to an embodiment of the present application. As shown in fig. 4, the apparatus includes:
The acquisition module 401 is configured to acquire vertical principal stress, horizontal maximum principal stress and horizontal minimum principal stress corresponding to faults of the hydrate aggregation area at different depths, where the hydrate aggregation area is a shallow hydrate aggregation area under an extended stress state;
A first determining module 402, configured to determine a sliding pressure minimum value corresponding to a fault of a hydrate aggregation zone according to a vertical principal stress, a horizontal maximum principal stress, and a horizontal minimum principal stress corresponding to the fault of the hydrate aggregation zone at different depths, where the sliding pressure minimum value is used to indicate a minimum pressure value required for sliding the fault of the hydrate aggregation zone;
And a second determining module 403, configured to determine a maximum bearable gas column height corresponding to the fault of the hydrate accumulation zone according to the sliding pressure minimum value.
Optionally, the vertical principal stress, the horizontal maximum principal stress, and the horizontal minimum principal stress corresponding to the fault of the hydrate collection zone at different depths each include a magnitude and a direction, and the first determining module 402 is specifically configured to:
Acquiring a three-dimensional model of the fault of the hydrate aggregation area, and determining corresponding dip angles of the fault of the hydrate aggregation area at different depths according to the three-dimensional model;
determining effective normal stress and shear stress corresponding to the fault of the hydrate aggregation area at different depths according to inclination angles, vertical main stress, horizontal maximum main stress and horizontal minimum main stress corresponding to the fault of the hydrate aggregation area at different depths;
And determining a sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to the effective normal stress and the shearing stress corresponding to the hydrate aggregation zone fault at different depths.
Optionally, the first determining module 402 is specifically configured to, when determining the effective normal stress and the shear stress corresponding to the hydrate collection zone fault at different depths according to the inclination angle, the vertical principal stress, the horizontal maximum principal stress, and the horizontal minimum principal stress corresponding to the hydrate collection zone fault at different depths:
determining the direction of the normal line of the unit surface corresponding to the fault of the hydrate aggregation area at different depths according to the corresponding inclination angles of the fault of the hydrate aggregation area at different depths;
Determining a first included angle between the vertical main stress and the normal line of the unit surface, a second included angle between the horizontal maximum main stress and the normal line of the unit surface and a third included angle between the horizontal minimum main stress and the normal line of the unit surface, which correspond to the hydrate aggregation area fault at different depths, according to the direction of the normal line of the unit surface, the direction of the vertical main stress, the direction of the horizontal maximum main stress and the direction of the horizontal minimum main stress, which correspond to the hydrate aggregation area fault at different depths;
and determining effective positive stress and shear stress corresponding to the fault of the hydrate aggregation area at different depths according to the vertical main stress, the horizontal maximum main stress, the horizontal minimum main stress, the first included angle, the second included angle and the third included angle corresponding to the fault of the hydrate aggregation area at different depths.
Optionally, the first determining module 402 is specifically configured to, when determining the effective normal stress and the shear stress corresponding to the fault of the hydrate collection area at different depths according to the magnitude of the vertical main stress, the magnitude of the horizontal maximum main stress, the magnitude of the horizontal minimum main stress, the first included angle, the second included angle, and the third included angle corresponding to the fault of the hydrate collection area at different depths:
Calculating effective positive stress and shear stress corresponding to the fault of the hydrate aggregation zone at different depths by using the following formula;
σn＝σvcos2αn+σHcos2βn+σhcos2γn
wherein σ n is effective normal stress, τ is shear stress, σ 1 is vertical main stress, σ 2 is horizontal maximum main stress, σ 3 is horizontal minimum main stress, α 1 is a first included angle, α 2 is a second included angle, and α 3 is a third included angle.
Optionally, the first determining module 402 is specifically configured to, when determining the sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to the effective normal stress and the shear stress corresponding to the hydrate aggregation zone fault at different depths:
obtaining the friction coefficient corresponding to the fault of the hydrate aggregation area, and according to a formula Determining sliding pressure values corresponding to the hydrate aggregation zone faults at different depths, wherein P s is the sliding pressure value, sigma n is effective positive stress, tau is shear stress, mu f is friction coefficient;
And determining the minimum sliding pressure corresponding to the fault of the hydrate aggregation zone according to the sliding pressure values corresponding to the fault of the hydrate aggregation zone at different depths.
Optionally, the second determining module 403 is specifically configured to:
Calculating the maximum bearable gas column height corresponding to the fault of the hydrate aggregation zone by using the following formula:
Wherein h g is the maximum bearable gas column height, ρ w is the sea water density of the hydrate aggregation area, ρ g is the gas density aggregated by the hydrate aggregation area fault, and P g is the sliding pressure minimum value corresponding to the hydrate aggregation area fault.
The specific implementation principle and effect of the fault maximum bearable air column height determining device provided by the embodiment of the application can be referred to the foregoing embodiment, and will not be repeated here.
Fig. 5 is a schematic structural diagram of an electronic device according to an embodiment of the present application. As shown in fig. 5, the electronic device of the present embodiment may include:
at least one processor 501, and
A memory 502 communicatively coupled to the at least one processor;
Wherein the memory 502 stores instructions executable by the at least one processor 501 to cause the electronic device to perform the method as described in any of the embodiments above.
Alternatively, the memory 502 may be separate or integrated with the processor 501.
The implementation principle and technical effects of the electronic device provided in this embodiment may be referred to the foregoing embodiments, and will not be described herein again.
The embodiment of the application also provides a computer readable storage medium, wherein computer executable instructions are stored in the computer readable storage medium, and when a processor executes the computer executable instructions, the method of any of the previous embodiments is realized.
Embodiments of the present application also provide a computer program product comprising a computer program which, when executed by a processor, implements a method as described in any of the preceding embodiments.
In the several embodiments provided by the present application, it should be understood that the disclosed apparatus and method may be implemented in other ways. For example, the above-described embodiments of the apparatus are merely illustrative, and for example, the division of the modules is merely a logical function division, and there may be additional divisions when actually implemented, for example, multiple modules may be combined or integrated into another system, or some features may be omitted or not performed.
The integrated modules, which are implemented in the form of software functional modules, may be stored in a computer readable storage medium. The software functional modules described above are stored in a storage medium and include instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) or processor to perform some of the steps of the methods described in the various embodiments of the application.
It should be appreciated that the Processor may be a processing unit (Central Processing Unit, CPU), but may also be other general purpose processors, digital signal processors (DIGITAL SIGNAL Processor, DSP), application SPECIFIC INTEGRATED Circuit (ASIC), etc. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of a method disclosed in connection with the present application may be embodied directly in a hardware processor for execution, or in a combination of hardware and software modules in a processor for execution. The Memory may include a high-speed random access Memory (Random Access Memory, RAM), and may further include a Non-Volatile Memory (NVM), such as at least one magnetic disk Memory, and may also be a U-disk, a removable hard disk, a read-only Memory, a magnetic disk, or an optical disk.
The storage medium may be implemented by any type of volatile or non-volatile Memory device or 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 (Programmable Read-Only Memory, PROM), read-Only Memory (ROM), magnetic Memory, flash Memory, magnetic disk or optical disk. A storage media may be any available media that can be accessed by a general purpose or special purpose computer.
An exemplary storage medium is coupled to the processor such the processor can read information from, and write information to, the storage medium. In the alternative, the storage medium may be integral to the processor. The processor and the storage medium may reside in an Application SPECIFIC INTEGRATED Circuits (ASIC). It is also possible that the processor and the storage medium reside as discrete components in an electronic device or a master device.
It should be noted that, in this document, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
The foregoing embodiment numbers of the present application are merely for the purpose of description, and do not represent the advantages or disadvantages of the embodiments.
From the above description of the embodiments, it will be clear to those skilled in the art that the above-described embodiment method may be implemented by means of software plus a necessary general hardware platform, but of course may also be implemented by means of hardware, but in many cases the former is a preferred embodiment. Based on such understanding, the technical solution of the present application may be embodied essentially or in a part contributing to the prior art in the form of a software product stored in a storage medium (e.g. ROM/RAM, magnetic disk, optical disk) comprising instructions for causing a terminal device (which may be a mobile phone, a computer, a server, an air conditioner, or a network device, etc.) to perform the method according to the embodiments of the present application.
The foregoing description is only of the preferred embodiments of the present application, and is not intended to limit the scope of the application, but rather is intended to cover any equivalents of the structures or equivalent processes disclosed herein or in the alternative, which may be employed directly or indirectly in other related arts.

Claims (5)

1. A method for determining the maximum sustainable air column height of a fault, comprising:
The method comprises the steps of obtaining vertical main stress, horizontal maximum main stress and horizontal minimum main stress corresponding to faults of a hydrate aggregation area at different depths, wherein the hydrate aggregation area is a shallow hydrate aggregation area in an extension stress state;
Determining a sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the hydrate aggregation zone fault at different depths, wherein the sliding pressure minimum value is used for indicating a minimum pressure value required for enabling the hydrate aggregation zone fault to slide;
Determining the maximum bearable gas column height corresponding to the fault of the hydrate aggregation area according to the minimum value of the sliding pressure;
The method for determining the sliding pressure minimum value corresponding to the fault of the hydrate aggregation zone according to the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the fault of the hydrate aggregation zone at different depths comprises the following steps:
Acquiring a three-dimensional model of the fault of the hydrate aggregation area, and determining corresponding dip angles of the fault of the hydrate aggregation area at different depths according to the three-dimensional model;
determining effective normal stress and shear stress corresponding to the fault of the hydrate aggregation area at different depths according to inclination angles, vertical main stress, horizontal maximum main stress and horizontal minimum main stress corresponding to the fault of the hydrate aggregation area at different depths;
determining a sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to the effective normal stress and the shearing stress corresponding to the hydrate aggregation zone fault at different depths;
The method for determining the effective normal stress and the shear stress of the hydrate aggregation zone fault at different depths according to the inclination angle, the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the hydrate aggregation zone fault at different depths comprises the following steps:
determining the direction of the normal line of the unit surface corresponding to the fault of the hydrate aggregation area at different depths according to the corresponding inclination angles of the fault of the hydrate aggregation area at different depths;
Determining a first included angle between the vertical main stress and the normal line of the unit surface, a second included angle between the horizontal maximum main stress and the normal line of the unit surface and a third included angle between the horizontal minimum main stress and the normal line of the unit surface, which correspond to the hydrate aggregation area fault at different depths, according to the direction of the normal line of the unit surface, the direction of the vertical main stress, the direction of the horizontal maximum main stress and the direction of the horizontal minimum main stress, which correspond to the hydrate aggregation area fault at different depths;
Determining effective positive stress and shear stress corresponding to the fault of the hydrate aggregation area at different depths according to the vertical main stress, the horizontal maximum main stress, the horizontal minimum main stress, the first included angle, the second included angle and the third included angle corresponding to the fault of the hydrate aggregation area at different depths;
According to the vertical main stress, the horizontal maximum main stress, the horizontal minimum main stress, the first included angle, the second included angle and the third included angle corresponding to the fault of the hydrate aggregation area at different depths, the effective normal stress and the shear stress corresponding to the fault of the hydrate aggregation area at different depths are determined, and the method comprises the following steps:
Calculating effective positive stress and shear stress corresponding to the fault of the hydrate aggregation zone at different depths by using the following formula;
σn＝σvcos2αn+σHcos2βn+σhcos2γn
Wherein, sigma n is effective normal stress, tau is shear stress, sigma v is vertical main stress, sigma H is horizontal maximum main stress, sigma h is horizontal minimum main stress, alpha 1 is a first included angle, alpha 2 is a second included angle, and alpha 3 is a third included angle;
the method for determining the sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to the effective normal stress and the shearing stress corresponding to the hydrate aggregation zone fault at different depths comprises the following steps:
obtaining the friction coefficient corresponding to the fault of the hydrate aggregation area, and according to a formula Determining sliding pressure values corresponding to the hydrate aggregation zone faults at different depths, wherein P s is the sliding pressure value, sigma n is effective positive stress, tau is shear stress, mu f is friction coefficient;
Determining a sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to sliding pressure values corresponding to the hydrate aggregation zone fault at different depths;
and determining the maximum bearable gas column height corresponding to the fault of the hydrate aggregation area according to the minimum value of the sliding pressure, wherein the method comprises the following steps:
Calculating the maximum bearable gas column height corresponding to the fault of the hydrate aggregation zone by using the following formula:
Wherein h g is the maximum bearable gas column height, ρ w is the sea water density of the hydrate aggregation area, ρ g is the gas density aggregated by the hydrate aggregation area fault, and P g is the sliding pressure minimum value corresponding to the hydrate aggregation area fault.
2. A fault maximum sustainable air column height determining apparatus, comprising:
the acquisition module is used for acquiring vertical main stress, horizontal maximum main stress and horizontal minimum main stress corresponding to faults of the hydrate aggregation area at different depths, wherein the hydrate aggregation area is a shallow hydrate aggregation area in an extension stress state;
The first determining module is used for determining a sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the hydrate aggregation zone fault at different depths, wherein the sliding pressure minimum value is used for indicating a minimum pressure value required for enabling the hydrate aggregation zone fault to slide;
The second determining module is used for determining the maximum bearable gas column height corresponding to the fault of the hydrate aggregation area according to the minimum value of the sliding pressure;
the hydrate aggregation zone fault comprises a vertical main stress, a horizontal maximum main stress and a horizontal minimum main stress which correspond to different depths, wherein the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress comprise the size and the direction, and the first determination module is specifically used for:
Acquiring a three-dimensional model of the fault of the hydrate aggregation area, and determining corresponding dip angles of the fault of the hydrate aggregation area at different depths according to the three-dimensional model;
determining effective normal stress and shear stress corresponding to the fault of the hydrate aggregation area at different depths according to inclination angles, vertical main stress, horizontal maximum main stress and horizontal minimum main stress corresponding to the fault of the hydrate aggregation area at different depths;
determining a sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to the effective normal stress and the shearing stress corresponding to the hydrate aggregation zone fault at different depths;
the first determining module is specifically configured to, when determining the effective normal stress and the shear stress corresponding to the fault of the hydrate aggregation area at different depths according to the inclination angle, the vertical main stress, the horizontal maximum main stress and the horizontal minimum main stress corresponding to the fault of the hydrate aggregation area at different depths:
determining the direction of the normal line of the unit surface corresponding to the fault of the hydrate aggregation area at different depths according to the corresponding inclination angles of the fault of the hydrate aggregation area at different depths;
Determining a first included angle between the vertical main stress and the normal line of the unit surface, a second included angle between the horizontal maximum main stress and the normal line of the unit surface and a third included angle between the horizontal minimum main stress and the normal line of the unit surface, which correspond to the hydrate aggregation area fault at different depths, according to the direction of the normal line of the unit surface, the direction of the vertical main stress, the direction of the horizontal maximum main stress and the direction of the horizontal minimum main stress, which correspond to the hydrate aggregation area fault at different depths;
Determining effective positive stress and shear stress corresponding to the fault of the hydrate aggregation area at different depths according to the vertical main stress, the horizontal maximum main stress, the horizontal minimum main stress, the first included angle, the second included angle and the third included angle corresponding to the fault of the hydrate aggregation area at different depths;
The first determining module is specifically configured to determine, when effective normal stress and shear stress corresponding to the fault of the hydrate aggregation area at different depths according to the magnitude of vertical main stress, the magnitude of horizontal maximum main stress, the magnitude of horizontal minimum main stress, the first included angle, the second included angle and the third included angle corresponding to the fault of the hydrate aggregation area at different depths:
Calculating effective positive stress and shear stress corresponding to the fault of the hydrate aggregation zone at different depths by using the following formula;
σn＝σvcos2αn+σHcos2βn+σhcos2γn
Wherein, sigma n is effective normal stress, tau is shear stress, sigma v is vertical main stress, sigma H is horizontal maximum main stress, sigma h is horizontal minimum main stress, alpha 1 is a first included angle, alpha 2 is a second included angle, and alpha 3 is a third included angle;
The first determining module is specifically configured to, when determining a sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to effective normal stress and shear stress corresponding to the hydrate aggregation zone fault at different depths:
obtaining the friction coefficient corresponding to the fault of the hydrate aggregation area, and according to a formula Determining sliding pressure values corresponding to the hydrate aggregation zone faults at different depths, wherein P s is the sliding pressure value, sigma n is effective positive stress, tau is shear stress, mu f is friction coefficient;
Determining a sliding pressure minimum value corresponding to the hydrate aggregation zone fault according to sliding pressure values corresponding to the hydrate aggregation zone fault at different depths;
the second determining module is specifically configured to, when determining a maximum sustainable air column height corresponding to the fault of the hydrate aggregation area according to the minimum sliding pressure value:
Calculating the maximum bearable gas column height corresponding to the fault of the hydrate aggregation zone by using the following formula:
Wherein h g is the maximum bearable gas column height, ρ w is the sea water density of the hydrate aggregation area, ρ g is the gas density aggregated by the hydrate aggregation area fault, and P g is the sliding pressure minimum value corresponding to the hydrate aggregation area fault.
3. An electronic device comprising a processor and a memory communicatively coupled to the processor;
the memory stores computer-executable instructions;
the processor executes computer-executable instructions stored in the memory to implement the method of claim 1.
4. A computer readable storage medium having stored therein computer executable instructions which when executed by a processor are adapted to implement the method of claim 1.
5. A computer program product comprising a computer program which, when executed by a processor, implements the method of claim 1.
CN202410710377.XA 2024-06-03 2024-06-03 Method, device, equipment and medium for determining the maximum tolerable gas column height of a fault Active CN118653803B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410710377.XA CN118653803B (en) 2024-06-03 2024-06-03 Method, device, equipment and medium for determining the maximum tolerable gas column height of a fault

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410710377.XA CN118653803B (en) 2024-06-03 2024-06-03 Method, device, equipment and medium for determining the maximum tolerable gas column height of a fault

Publications (2)

Publication Number Publication Date
CN118653803A CN118653803A (en) 2024-09-17
CN118653803B true CN118653803B (en) 2025-02-28

Family

ID=92707140

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410710377.XA Active CN118653803B (en) 2024-06-03 2024-06-03 Method, device, equipment and medium for determining the maximum tolerable gas column height of a fault

Country Status (1)

Country Link
CN (1) CN118653803B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109057757A (en) * 2018-08-24 2018-12-21 广州海洋地质调查局 A kind of gas hydrate mining methods and device
CN113656744A (en) * 2021-08-16 2021-11-16 中国科学院广州能源研究所 A method for estimation of pit depth and free gas layer evaluation in marine hydrate area

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8074720B2 (en) * 2004-09-28 2011-12-13 Vetco Gray Inc. Riser lifecycle management system, program product, and related methods
CN117993327B (en) * 2024-03-27 2024-06-11 广州海洋地质调查局三亚南海地质研究所 Reservoir simulation method, device, equipment and medium for seepage type natural gas hydrate

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109057757A (en) * 2018-08-24 2018-12-21 广州海洋地质调查局 A kind of gas hydrate mining methods and device
CN113656744A (en) * 2021-08-16 2021-11-16 中国科学院广州能源研究所 A method for estimation of pit depth and free gas layer evaluation in marine hydrate area

Also Published As

Publication number Publication date
CN118653803A (en) 2024-09-17

Similar Documents

Publication Publication Date Title
Robertson Cone penetration test (CPT)-based soil behaviour type (SBT) classification system—an update
CA2690992C (en) Method for predicting well reliability by computer simulation
US8423337B2 (en) Method for multi-scale geomechanical model analysis by computer simulation
US8548782B2 (en) Method for modeling deformation in subsurface strata
US8868390B2 (en) Method of developing a petroleum reservoir from a flow model calibrated through pressure and saturation map scaling
US20160349389A1 (en) Method for developing a geomechanical model based on seismic data, well logs and sem analysis of horizontal and vertical drill cuttings
WO2012071103A1 (en) System and method for estimating fluid distribution in a subterranean reservoir
CN109509111B (en) Method and system for predicting formation pressure of exploratory well
Sturm On the liquefaction potential of gravelly soils: Characterization, triggering and performance
Bunn et al. Geologic trends in shear strength properties inferred through three‐dimensional back analysis of landslide inventories
US8718992B2 (en) Method for history matching of a geological model comprising a sub-seismic fault network
Falcon‐Suarez et al. CO2‐brine substitution effects on ultrasonic wave propagation through sandstone with oblique fractures
US20150219792A1 (en) Simulation of insular karstification
Chakraborty et al. Seismic slope stability analysis of a hydraulic fill dam
CN118653803B (en) Method, device, equipment and medium for determining the maximum tolerable gas column height of a fault
CN118779559B (en) Method, device, equipment and medium for determining maximum bearable gas column height of cover layer
Steer et al. 2D and 3D modeling of rock mechanical properties of Khasib formation in East Baghdad oil field
Dudley et al. Predicting accelerating subsidence above the highly compacting Luconia carbonate reservoirs, offshore Sarawak Malaysia
Younessi et al. Evaluation of reservoir compaction and geomechanical related risks associated with reservoir production using a two-way coupled 4D geomechanical model for a deep-water field in South East Asia
Ringrose et al. The property model
US20220342105A1 (en) Method for validating rock formations compaction parameters using geomechanical modeling
Lescoffit et al. Quantifying the impact of fault modeling parameters on production forecasting for clastic reservoirs
Damico et al. Challenging geostatistical methods to represent heterogeneity in CO2 reservoirs under residual trapping
CN113514890A (en) Method, device and equipment for predicting formation pressure coefficient by using seismic data
Petrie et al. Modeling strain across mechanical sedimentary lithologic interfaces: geomechanical models derived from outcrop analysis

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant