Prediction method and device for maximum structural paleostress

Abstract

The invention relates to the field of petroleum and natural gas, in particular to a method and device for predicting the maximum structural paleostress. The method comprises the following steps: acquiring logging data, selecting a target interval, and determining a to-be-predicted depth point in the target interval according to the logging data and the target interval; obtaining logging calculation effective paleostress and a first included angle of the to-be-predicted depth point; the first included angle is an arc tangent function of the current maximum horizontal stress and the current vertical stress; and determining the maximum structural paleostress of the to-be-predicted depth point based on the first included angle and the logging calculation effective paleostress. According to the method, quantitative calculation is carried out on the maximum structural paleostress of the stratum through conventional logging information, prediction of the maximum structural paleostress is achieved, the unavailability of relying on core testing is reduced, and meanwhile the multiplicity and uncertainty of calculation relying on fracture density, drilling leakage conditions, geophysical data and numerical simulation are filled up.

Description

Technical Field

[0001] The invention relates to the field of petroleum and natural gas, and in particular to a method and a device for predicting maximum tectonic paleostress. Background Art

[0002] Paleotectonic stress refers to the stress of tectonic activity during the geological history, which is determined by the deformation traces left by tectonic activity (macroscopic deformations such as folds, faults, and joints), and can be determined by core-based methods, borehole-based methods, geological methods, geophysical methods, and methods based on underground space. Present-day tectonic stress refers to the tectonic stress currently in action (Mid-Pleistocene of the Quaternary Period), which is an extension or continuation of paleotectonic stress.

[0003] The calculation method of the in-situ maximum paleo-tectonic stress logging data of the formation has always been a hot and difficult issue in oil and gas exploration and development. Due to the complex geological background of the foreland thrust belt and the sensitivity of logging data, it has not been fundamentally solved so far. The physical simulation and mathematical prediction of formation stress that have been carried out by predecessors are mostly to simulate and predict the paleo-stress and present-day stress characteristics of the formation under drilling or outcrop conditions. There is no comprehensive calculation method for the maximum paleo-stress logging of the drilling formation in the strong compression zone of the thrust belt. Summary of the invention

[0004] In order to solve the problem that there are few comprehensive logging calculation methods for the maximum paleostress of drilling formations in strong compression zones of thrust belts, the present invention provides a prediction method and device for the maximum tectonic paleostress, which can achieve quantitative prediction of the maximum tectonic paleostress of in-situ formations based on logging data when the exploration degree is low, drilling core data is lacking and the resolution of seismic data is limited, so as to effectively guide the exploration and development of deep and ultra-deep oil and gas, and fill the gap in the quantitative calculation of tectonic paleostress of formations.

[0005] In order to solve any of the above technical problems, the specific technical solutions of the embodiments of this specification are as follows:

[0006] In one aspect, the present invention provides a method for predicting maximum tectonic paleostress, the method comprising:

[0007] Acquire well logging data and select a target layer segment, and determine a depth point to be predicted in the target layer segment according to the well logging data and the target layer segment;

[0008] Obtaining the logging calculated effective paleostress and the first angle of the depth point to be predicted; the first angle is the inverse tangent function of the current maximum horizontal stress and the current vertical stress;

[0009] The maximum structural paleostress at the depth point to be predicted is determined based on the first angle and the effective paleostress calculated by logging.

[0010] Further, acquiring well logging data and selecting a target layer segment, and determining a depth point to be predicted in the target layer segment according to the well logging data and the target layer segment includes:

[0011] Determine stratigraphic information using seismic depth profile geological interpretation data;

[0012] Obtain well logging data based on regional geological information;

[0013] A target layer section is selected based on the stratigraphic layer information and the logging data, and a depth point to be predicted in the target layer section is determined.

[0014] Furthermore, the depth point ΔGR to be predicted in the target layer segment is obtained according to the following formula:

[0015] ΔGR=(GR-GRmin) / (GRmax-GRmin),

[0016] Among them, GR is the natural gamma logging value, GRmin is the minimum value of natural gamma logging, and GRmax is the maximum value of natural gamma logging.

[0017] Furthermore, the effective paleo-stress S calculated by logging at the depth point to be predicted is obtained according to the following formula: v :

[0018] S v =8.449*LogRT / GR-2.2394*AC+205.447,

[0019] Among them, AC is the logging acoustic wave time difference value, and RT is the logging resistivity value.

[0020] Further, the first angle θ is obtained according to the following formula:

[0021] θ=arctg(P o / σ H )

[0022] Among them, P o is the current vertical stress, σ H The current maximum horizontal stress.

[0023] Furthermore, the current vertical stress P is obtained according to the following formula: o and the current maximum horizontal stress σ H :

[0024]

[0025]

[0026]

[0027] Where TVD is the vertical burial depth, ρ b is the logging density value, p p is the formation pore pressure, α is the Eaton coefficient, which is 1; n is the Eaton index, which is 9, and p pn is the normal compaction pore pressure; p0 is the overlying formation pressure; Δt0 is the mudstone logging acoustic time difference at the depth to be predicted; Δt n is the acoustic time difference of the normal trend line of mudstone corresponding to the calculation point of the depth point to be predicted; ν is the Poisson's ratio; ε H is the correction coefficient of the maximum horizontal principal stress; h is the horizontal minimum principal stress correction coefficient; E is Young's modulus.

[0028] Furthermore, the maximum tectonic paleo-stress S at the depth point to be predicted is determined according to the following formula: max :

[0029] S max =S v *sinθ.

[0030] In another aspect, the present invention provides a device for predicting maximum tectonic paleostress, the device comprising:

[0031] A first acquisition unit is used to acquire well logging data and select a target layer segment, and determine a depth point to be predicted in the target layer segment according to the well logging data and the target layer segment;

[0032] A second acquisition unit is used to acquire the logging calculated effective paleostress and a first angle of the depth point to be predicted; the first angle is an inverse tangent function of the current maximum horizontal stress and the current vertical stress;

[0033] A calculation unit is used to determine the maximum structural paleostress of the depth point to be predicted based on the first angle and the effective paleostress calculated by logging.

[0034] On the other hand, the present invention further provides a computer device, comprising a memory, a processor, and a computer program stored in the memory, wherein the processor implements the above method when executing the computer program.

[0035] On the other hand, the present invention provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and the computer program implements the above method when executed by a processor.

[0036] Finally, the present invention also provides a computer program product, which includes a computer program, and the computer program implements the above method when executed by a processor.

[0037] The maximum tectonic paleostress prediction method of the present invention uses conventional logging data to quantitatively calculate the maximum tectonic paleostress of the formation, and the prediction of the maximum tectonic paleostress of deep and ultra-deep formations is achieved by establishing an effective paleostress model for logging calculation and the conversion of wellbore direction paleostress to vertical stress. This reduces the unavailability of core testing and also fills the multi-solution and uncertainty of calculations based on fracture density, drilling leakage, geophysical data and numerical simulation. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] In order to more clearly illustrate the embodiments of this specification or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the embodiments of this specification. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0039] Figure 1 The figure shows a flow chart of a method for predicting maximum tectonic paleostress in an embodiment of the present specification;

[0040] Figure 2 It is a schematic diagram of the coordinate axes of the current maximum horizontal stress and vertical stress of the wellbore in the embodiment of this specification;

[0041] Figure 3 It is a structural schematic diagram of a shared energy storage optimization configuration device in an embodiment of this specification;

[0042] Figure 4 It is a schematic diagram showing the correlation between the predicted maximum tectonic paleostress and the measured effective paleostress in the embodiments of this specification;

[0043] Figure 5 The figure shows the maximum structural paleostress calculation bar graph of the Jurassic system in Yinan 4 well of Kuche Depression in the embodiment of this specification;

[0044] Figure 6 The figure is a schematic diagram of the structure of a computer device in an embodiment of the present specification.

[0045]

Description of reference numerals

[0046] 201. A first acquiring unit;

[0047] 202. A second acquisition unit;

[0048] 203, computing unit;

[0049] 602. Computer equipment;

[0050] 604, processor;

[0051] 606. Memory;

[0052] 608, driving mechanism;

[0053] 610, input / output module;

[0054] 612. Input devices;

[0055] 614. Output device;

[0056] 616. Presentation equipment;

[0057] 618. Graphical user interface;

[0058] 620, network interface;

[0059] 622, communication link;

[0060] 624. Communication bus. DETAILED DESCRIPTION

[0061] The following will be combined with the drawings in the embodiments of this specification to clearly and completely describe the technical solutions in the embodiments of this specification. Obviously, the described embodiments are only part of the embodiments of this specification, not all of them. Based on the embodiments in the embodiments of this specification, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the embodiments of this specification.

[0062] It should be noted that the terms "first", "second", etc. in the description and claims of the embodiments of this specification and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the embodiments of this specification described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, device, product or equipment that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or equipment.

[0063] It should be noted that the acquisition, storage, use, and processing of data in the technical solution of this application are in compliance with the relevant provisions of relevant laws and regulations. The order of steps listed in the embodiment is only one way of executing the steps among many steps, and does not represent the only execution order. When the actual system or device product is executed, it can be executed in the order or in parallel according to the method shown in the embodiment or the accompanying drawings.

[0064] Mudstone and shale are the least compressive rocks among sedimentary rocks and are the most sensitive to stress. Previous studies have shown that there is basically no montmorillonite with strong cation exchange capacity in the Jurassic mudstone in the western piedmont, and more water will be discharged during tectonic compression. Therefore, the relationship between compression geostress and electrical parameters can be established, and the electrical parameters can be obtained through logging data. However, according to analysis, it can be seen that the response of acoustic waves and resistivity to the maximum paleostress value is poor. Therefore, in the embodiment of the present invention, the maximum paleostress value is not calculated by the single-factor correlation between acoustic waves or resistivity and the maximum paleostress value, but by establishing a maximum effective paleostress calculation model, a model of the angle between the current horizontal maximum principal stress and the vertical stress of the wellbore, and a logging maximum tectonic paleostress model to achieve the calculation of the maximum paleotectonic backlog stress. The specific implementation method is as follows.

[0065] Through the interactive analysis of logging data and measured maximum paleostress values, it is found that acoustic waves and resistivity respond poorly to the maximum paleostress values.

[0066] A specific embodiment of the present invention discloses a method for predicting maximum tectonic paleostress, specifically as follows: Figure 1 As shown, the method includes:

[0067] Step 101, acquiring logging data and selecting a target layer segment, and determining a depth point to be predicted in the target layer segment according to the logging data and the target layer segment. Specifically, the logging data can be acquired from known conventional logging data and the target layer segment can be selected.

[0068] Step 102, obtaining the effective paleostress calculated by logging at the depth point to be predicted and the first angle; the first angle is the inverse tangent function of the current maximum horizontal stress and the current vertical stress. Specifically, Figure 2 As shown, a well logging calculation effective paleostress model and a model of the angle between the current maximum horizontal stress and the current vertical stress are constructed to obtain the well logging calculation effective paleostress and the first angle of the depth point to be predicted. The first angle is the inverse tangent function of the current maximum horizontal stress and the current vertical stress, that is, the inclination angle of the formation, which can be converted into the angle between the well logging calculation effective paleostress and the maximum structural paleostress.

[0069] Step 103, determining the maximum tectonic paleostress at the depth point to be predicted based on the first angle and the effective paleostress calculated by logging. Specifically, after obtaining the effective paleostress calculated by logging and the first angle in step 102, the maximum tectonic paleostress at the depth point to be predicted can be determined based on the angle relationship formed between the effective paleostress calculated by logging and the maximum tectonic paleostress and the inverse tangent function. The maximum tectonic paleostress is also the maximum horizontal tectonic paleostress.

[0070] A method for predicting maximum tectonic paleostress in an embodiment of the present invention is based on new ideas and new methods, uses conventional logging data to quantitatively calculate the maximum tectonic paleostress of the formation, and realizes the prediction of the maximum tectonic paleostress of deep and ultra-deep formations by establishing an effective paleostress model for logging calculation and wellbore direction paleostress-vertical stress conversion. It reduces the unavailability of core testing, and also fills the multi-solution and uncertainty of calculations based on fracture density, drilling leakage, geophysical data and numerical simulation.

[0071] In one embodiment of the present specification, acquiring well logging data and selecting a target layer segment, and determining a depth point to be predicted in the target layer segment according to the well logging data and the target layer segment includes:

[0072] The stratigraphic information is determined by using the geological interpretation data of the seismic depth profile. Specifically, the conventional logging data of key wells is relatively abundant, especially in oil and gas basins. In addition to drilling data, stratigraphic information can also be obtained from the geological interpretation data of the seismic depth profile. On this basis, stratigraphic information can be determined.

[0073] The logging data is obtained based on the regional geological data; the logging data is relatively abundant, and the required logging data can be obtained on demand according to the subsequent calculation requirements. At least it can include: natural gamma logging value, logging acoustic wave time difference value and logging resistivity value.

[0074] Based on the formation layer information and logging data, a target layer section is selected and a depth point to be predicted in the target layer section is determined. After the target layer section is selected, the depth point to be predicted in the target layer section can be determined based on the logging data of the known target layer section obtained above.

[0075] In one embodiment of the present specification, the depth point ΔGR to be predicted in the target layer segment is obtained according to the following formula:

[0076] ΔGR=(GR-GRmin) / (GRmax-GRmin),

[0077] Among them, GR is the natural gamma logging value, GRmin is the minimum value of natural gamma logging, and GRmax is the maximum value of natural gamma logging.

[0078] In one embodiment of the present specification, a well logging calculation effective paleo-stress model is established, and the well logging calculation effective paleo-stress S of the depth point to be predicted is obtained according to the following formula: v :

[0079] S v =8.449*LogRT / GR-2.2394*AC+205.447,

[0080] Among them, AC is the logging acoustic wave time difference value, and RT is the logging resistivity value.

[0081] Among them, in the above embodiment of the present invention, the purpose of introducing the predicted depth point ΔGR=(GR-GRmin) / (GRmax-GRmin) and adopting the calculation method of LogRT / ΔGR is to correct the mud content and fluid of the response of mudstone stress, because the predicted depth point and LogRT / ΔGR have a good correlation with the maximum tectonic paleostress required to be obtained by the present invention.

[0082] In one embodiment of the present specification, the first angle θ is obtained according to the following formula:

[0083] θ=arctg(P o / σ H )

[0084] Among them, P o is the current vertical stress, σ H The current maximum horizontal stress.

[0085] In one embodiment of the present specification, the present vertical stress P is obtained according to the following formula: o and the current maximum horizontal stress σ H :

[0086]

[0087]

[0088]

[0089] Where TVD is the vertical burial depth, ρ b is the logging density value, p p is the formation pore pressure, α is the Eaton coefficient, which is 1; n is the Eaton index, which is 9, and p pn is the normal compaction pore pressure; p0 is the overlying formation pressure; Δt0 is the mudstone logging acoustic time difference at the depth to be predicted; Δt n is the acoustic time difference of the normal trend line of mudstone corresponding to the calculation point of the depth point to be predicted; ν is the Poisson's ratio; ε H is the correction coefficient of the maximum horizontal principal stress; h is the horizontal minimum principal stress correction coefficient; E is Young's modulus; g is the gravitational acceleration.

[0090] In a specific embodiment, since the dynamic rock mechanical physical parameters at a certain depth point are mainly related to the rock density and the P-wave and S-wave time difference at that point, the following formula is introduced to determine the logging density value ρ b and the formation pore pressure p p, and then determine the current horizontal stress σ H :

[0091]

[0092]

[0093]

[0094]

[0095] Among them, Δt0 is the mudstone logging acoustic time difference at the depth point to be predicted, μs / ft; Δt n is the acoustic time difference of the normal trend line of the mudstone corresponding to the depth point to be predicted, μs / ft; n is the Eaton index, which is 9 in this embodiment; p pn is the normal compaction pore pressure, MPa; p0 is the overlying formation pressure, MPa; ρ b is the density of the overlying rock layer, ×10 3 kg / m 3 ρ f is the fluid density, ×10 3 kg / m 3 ; G dyn is the dynamic shear modulus; K dyn is the dynamic bulk modulus; v dyn is the dynamic Poisson's ratio; E dyn is the dynamic Young's modulus; DTSM is the shear wave time difference, μs / ft; DTCO is the longitudinal wave time difference, μs / ft.

[0096] In one embodiment of the present specification, the maximum tectonic paleo-stress S at the depth point to be predicted is determined according to the following formula: max :

[0097] S max =S v *sinθ.

[0098] Based on the same inventive concept, the embodiment of this specification also provides a prediction device for maximum tectonic paleostress, specifically as follows Figure 3 As shown, the device comprises:

[0099] The first acquisition unit 201 is used to acquire well logging data and select a target layer segment, and determine a depth point to be predicted in the target layer segment according to the well logging data and the target layer segment;

[0100] The second acquisition unit 202 is used to acquire the logging calculated effective paleostress and the first angle of the depth point to be predicted; the first angle is the inverse tangent function of the current maximum horizontal stress and the current vertical stress;

[0101] The calculation unit 203 is used to determine the maximum structural paleostress at the depth point to be predicted based on the first angle and the effective paleostress calculated by logging.

[0102] Since the principle of solving the problem by the above device is similar to that of the above method, the implementation of the above device can refer to the implementation of the above method, and the repeated parts will not be repeated.

[0103] In one embodiment of the present specification, the maximum paleo-tectonic stress of the formations of Yinan 4, Dibei 102 and Dixi 1 wells in the Dibei area of ​​the Kuche Depression is calculated. The effective paleo-stress calculated by well logging has a good correlation with the measured paleo-stress. Figure 4 As shown in the figure, it is shown that conventional logging data can be used to calculate the maximum tectonic paleostress of the well. The horizontal axis is the maximum value of the measured paleostress, that is, the maximum tectonic paleostress; the vertical axis is the effective paleostress calculated by logging. R represents the correlation between the effective paleostress calculated by logging and the measured paleostress.

[0104] Figure 5 This is a bar chart of the maximum structural paleostress calculation of the Jurassic in Yinan 4 Well in Kuche Depression. It can be seen that the effective paleostress calculated by logging is highly correlated with the maximum horizontal structural paleostress.

[0105] like Figure 6 As shown, a computer device provided in an embodiment of the present invention, the apparatus in this specification may be a computer device in this embodiment, and the method of the above-mentioned specification is executed. The computer device 602 may include one or more processors 604, such as one or more central processing units (CPUs), and each processing unit may implement one or more hardware threads. The computer device 602 may also include any memory 606, which is used to store any kind of information such as code, settings, data, etc. Non-limiting, for example, the memory 606 may include any one or more combinations of the following: any type of RAM, any type of ROM, flash memory device, hard disk, optical disk, etc. More generally, any memory may use any technology to store information. Further, any memory may provide volatile or non-volatile retention of information. Further, any memory may represent a fixed or removable component of the computer device 602. In one case, when the processor 604 executes an associated instruction stored in any memory or a combination of memories, the computer device 602 may perform any operation of the associated instruction. The computer device 602 also includes one or more drive mechanisms 608 for interacting with any memory, such as a hard disk drive mechanism, an optical disk drive mechanism, etc.

[0106] The computer device 602 may also include an input / output module 610 (I / O) for receiving various inputs (via input devices 612) and for providing various outputs (via output devices 614). A specific output mechanism may include a presentation device 616 and an associated graphical user interface (GUI) 618. In other embodiments, the input / output module 610 (I / O), the input device 612, and the output device 614 may not be included, and the computer device 602 may be used as a computer device in a network. The computer device 602 may also include one or more network interfaces 620 for exchanging data with other devices via one or more communication links 622. One or more communication buses 624 couple the components described above together.

[0107] The communication link 622 may be implemented in any manner, for example, through a local area network, a wide area network (e.g., the Internet), a point-to-point connection, etc., or any combination thereof. The communication link 622 may include any combination of hardwired links, wireless links, routers, gateway functions, name servers, etc. governed by any protocol or combination of protocols.

[0108] The embodiments of the present specification also provide a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and the computer program implements the above method when executed by a processor.

[0109] The embodiments of the present specification also provide a computer-readable instruction, wherein when a processor executes the instruction, the program therein causes the processor to execute the above method.

[0110] It should be understood that in the various embodiments of the present specification, the size of the serial numbers of the above-mentioned processes does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present specification.

[0111] It should also be understood that in the embodiments of this specification, the term "and / or" is only a description of the association relationship of the associated objects, indicating that three relationships can exist. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. In addition, the character " / " in the embodiments of this specification generally indicates that the associated objects before and after are in an "or" relationship.

[0112] Those of ordinary skill in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed in the embodiments of this specification can be implemented with electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to exceed the scope of the embodiments of this specification.

[0113] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0114] In the several embodiments provided in the embodiments of this specification, it should be understood that the disclosed systems, 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 units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices or units, or it can be an electrical, mechanical or other form of connection.

[0115] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the embodiments of this specification.

[0116] In addition, each functional unit in each embodiment of the present specification can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The above integrated unit can be implemented in the form of hardware or in the form of software functional units.

[0117] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the embodiment of this specification is essentially or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the embodiment of this specification. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program code.

[0118] The embodiments of this specification use specific embodiments to illustrate the principles and implementation methods of the embodiments of this specification. The description of the above embodiments is only used to help understand the methods and core ideas of the embodiments of this specification. At the same time, for those skilled in the art, according to the ideas of the embodiments of this specification, there will be changes in the specific implementation methods and application scopes. In summary, the content of this specification should not be understood as a limitation on the embodiments of this specification.

Claims

1. A method for predicting maximum tectonic paleostress, characterized in that: The method comprises: Acquire well logging data and select a target layer segment, and determine a depth point to be predicted in the target layer segment according to the well logging data and the target layer segment; Obtaining the logging calculated effective paleostress and the first angle of the depth point to be predicted; the first angle is the inverse tangent function of the current maximum horizontal stress and the current vertical stress; The maximum structural paleostress at the depth point to be predicted is determined based on the first angle and the effective paleostress calculated by logging.

2. A method for predicting maximum tectonic paleostress according to claim 1, characterized in that: Acquiring well logging data and selecting a target layer segment, and determining a depth point to be predicted in the target layer segment according to the well logging data and the target layer segment includes: Determine stratigraphic information using geological interpretation data from seismic depth profiles; Obtain well logging data based on regional geological information; A target layer section is selected based on the stratigraphic layer information and the logging data, and a depth point to be predicted in the target layer section is determined.

3. A method for predicting maximum tectonic paleostress according to claim 2, characterized in that: The depth point ΔGR to be predicted in the target layer segment is obtained according to the following formula: ΔGR=(GR-GRmin) / (GRmax-GRmin), Among them, GR is the natural gamma logging value, GRmin is the minimum value of natural gamma logging, and GRmax is the maximum value of natural gamma logging.

4. A method for predicting maximum tectonic paleostress according to claim 3, characterized in that: The effective paleo-stress S calculated by logging at the depth point to be predicted is obtained according to the following formula v : S v =8.449*LogRT / GR-2.2394*AC+205.447, Among them, AC is the logging acoustic wave time difference value, and RT is the logging resistivity value.

5. A method for predicting maximum tectonic paleostress according to claim 4, characterized in that: The first angle θ is obtained according to the following formula: θ=arctg(P o / σ H ) Among them, P o is the current vertical stress, σ H The current maximum horizontal stress.

6. A method for predicting maximum tectonic paleostress according to claim 5, characterized in that: The current vertical stress P is obtained according to the following formula: o and the current maximum horizontal stress σ H : Where TVD is the vertical burial depth, ρ b is the logging density value, p p is the formation pore pressure, α is the Eaton coefficient, which is 1; n is the Eaton index, which is 9, and p pn is the normal compaction pore pressure; p0 is the overlying formation pressure; Δt0 is the mudstone logging acoustic time difference at the depth to be predicted; Δt n is the acoustic time difference of the normal trend line of mudstone corresponding to the depth point to be predicted; ν is the Poisson's ratio; ε H is the correction coefficient of the maximum horizontal principal stress; h is the horizontal minimum principal stress correction coefficient; E is Young's modulus.

7. A method for predicting maximum tectonic paleostress according to claim 5, characterized in that: The maximum tectonic paleo-stress S at the depth point to be predicted is determined according to the following formula: max : S max =S v *sinθ。 8. A device for predicting maximum tectonic paleostress, characterized in that: The device comprises: A first acquisition unit is used to acquire well logging data and select a target layer segment, and determine a depth point to be predicted in the target layer segment according to the well logging data and the target layer segment; A second acquisition unit is used to acquire the logging calculated effective paleostress and a first angle of the depth point to be predicted; the first angle is an inverse tangent function of the current maximum horizontal stress and the current vertical stress; A calculation unit is used to determine the maximum structural paleostress of the depth point to be predicted based on the first angle and the effective paleostress calculated by logging.

9. A computer device comprising a memory, a processor, and a computer program stored in the memory, characterized in that: When the processor executes the computer program, the method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images