Method for dividing three-dimensional rock mechanical layer through well-seismic combination, product and application

Abstract

The invention provides a method for dividing a three-dimensional rock mechanical layer through well-seismic combination, a product and application, and belongs to the technical field of geological information processing. The method for dividing the three-dimensional rock mechanical layer through well-seismic combination comprises the following steps that a rock mechanical index curve is generated according to longitudinal wave time difference and density logging information; dividing a rock mechanics layer according to the rock mechanics index curve; and well seismic calibration is carried out according to the divided rock mechanical layer. According to the method, small, medium and large-scale rock mechanical layer division of a single well and three-dimensional medium and large-scale rock mechanical layer division on an earthquake scale are completed through calculation of a single well rock mechanical index and well earthquake calibration; the problem of division of the three-dimensional rock mechanical layer is solved, and the method can be used for analyzing the control effect of the rock mechanical layer on fracture-crack formation distribution.

Description

Technical Field

[0001] This invention belongs to the field of geological information processing technology, and relates to a method, product and application of well-seismic combined method for dividing three-dimensional rock mechanical layers. Background Technology

[0002] A rock mechanical layer refers to a rock stratum with similar or identical rock mechanical behavior and properties. It is a rock mechanical unit defined by rock mechanical interfaces, possessing similar or identical strength, brittleness, and fracture mechanical properties. It comprises two parts: the rock mechanical unit and the rock mechanical interface. In ultra-deep carbonate rocks, the distribution characteristics of fractures at different scales are key factors determining reservoir formation and high-yield hydrocarbon accumulation. Stress concentration at the tip of secondary fractures controls their extension path and velocity. When encountering strata with lower brittleness or interfaces with significant changes in mechanical properties, the original stress cannot support further fracture propagation, and fracture extension is limited by the rock mechanical interface. Therefore, rock mechanical layers control the formation and distribution of secondary fractures. Delineating rock mechanical layers is of great significance for elucidating the formation and distribution patterns of secondary fracture systems.

[0003] The most commonly used multi-scale rock mechanical layering methods are as follows:

[0004] (1) Petrological methods: The density of fracture development varies greatly among different lithologies. For example, in carbonate strata, dolomite has a higher fracture density than dolomitic limestone, while limestone has the lowest. Fine-grained carbonate strata have a higher fracture density than coarse-grained ones. Therefore, the interface where lithology changes significantly can be used as the basis for dividing rock mechanical layers.

[0005] (2) Sequence stratigraphy method: The unconformity at the top and bottom of the sequence and the different levels of flooding surfaces inside are obvious interfaces of mechanical property change, which can be used as the main basis for dividing the rock mechanical layers.

[0006] (3) Measured rock mechanical parameters method: The rebound hammer can be used to measure and estimate the mechanical parameters such as uniaxial compressive strength, Young's modulus, formation energy or shear modulus (stiffness) of field samples and cores, which are used to calibrate and divide rock mechanical layers.

[0007] (4) Well logging interpretation of mechanical parameters: Various rock mechanical parameters can be directly calculated through well logging curves such as sonic waves and density. By analogy with the well logging method for dividing formations, various rock mechanical indices can be used to divide rock mechanical layers.

[0008] Patent CN114184764A provides a method and system for delineating the rock mechanical layers of tight carbonate reservoirs. The method includes: acquiring static rock mechanical parameters from multiple core samples; acquiring dynamic rock mechanical parameters at the corresponding depth point for each core sample; acquiring corrected static rock mechanical parameters for each depth point based on the static rock mechanical parameters of each core sample, the dynamic rock mechanical parameters at the corresponding depth point, and the dynamic rock mechanical parameters at each depth point; calculating the rock mechanical index for each depth point based on the corrected static rock mechanical parameters and a reservoir rock mechanical index calculation model; and delineating the rock mechanical layers of the reservoir based on the fracture development characteristics of a single well and the rock mechanical index for each depth point. This technology utilizes conventional logging data to calculate rock mechanical indices with high accuracy and computational efficiency, thus enabling rapid, accurate, and reliable delineation of the rock mechanical layers of tight carbonate reservoirs.

[0009] Patent CN116608931A provides a method for delineating rock mechanical layers in deep tight sandstone reservoirs. The method includes the following steps: calculating dynamic rock mechanical parameters (static Young's modulus and dynamic Poisson's ratio) through well logging; converting dynamic and static rock mechanical parameters to obtain static Young's modulus and static Poisson's ratio; and constructing a rock mechanical layer discrimination index model for rock mechanical layer delineation. This invention involves calculating dynamic rock mechanical parameters; establishing a dynamic-to-static conversion model for rock mechanical parameters of the target stratum; constructing a rock mechanical layer discrimination index model; and delineating rock mechanical layers by setting different rock mechanical layer discrimination index thresholds D. This technology proposes a method for delineating rock mechanical layers in deep tight sandstone reservoirs, which has high practical value. The delineation process is simple and clear, and the identification results have practical value for oil and gas exploration and development in oil and gas basins.

[0010] Patent CN112765785A provides a multi-scale well logging method for delineating rock mechanical layers. It constructs a rock mechanical layer discrimination index by calculating dynamic and static rock mechanical parameters and employing an isofrequency transformation method. Through cluster analysis of rock mechanical parameters from well logging points, it delineates rock mechanical layers at a single scale and calculates rock density. By changing the threshold and iteratively performing cluster analysis of rock mechanical parameters from well logging points, it completes the conventional well logging delineation of rock mechanical layers at different scales. This technology has practical value for reservoir geomechanical modeling and stress field numerical simulation.

[0011] Patent CN114428363A provides a method, device, electronic device, and medium for delineating rock mechanical layers in clastic strata. The method for delineating rock mechanical layers in clastic strata includes: dividing a lithological profile; obtaining static rock mechanical parameters; obtaining dynamic rock mechanical parameters through calculation; obtaining the vertical distribution law of rock mechanical properties based on the static and dynamic rock mechanical parameters; and delineating rock mechanical layers according to the lithological profile and the vertical distribution law of rock mechanical properties. This method can be widely applied in strata in my country where capable and incapable layers (such as sandstone and mudstone interbedded) are interbedded in weakly deformable tectonic zones.

[0012] However, the existing rock mechanics layer division methods commonly used above are limited to the division of vertical rock mechanics interfaces and lack a method for dividing three-dimensional rock mechanics layers in space. Summary of the Invention

[0013] In view of this, and in response to the bottleneck problem of the lack of a spatial three-dimensional rock mechanics layer division method in the existing technology, the purpose of this invention is to provide a method, product and application for dividing three-dimensional rock mechanics layers by combining well and seismic analysis, so as to better analyze the distribution law and cause of fractures and cracks.

[0014] To achieve the above-mentioned objectives, this invention provides a method for dividing three-dimensional rock mechanical layers using a combination of well and seismic analysis, comprising the following steps:

[0015] Rock mechanics index curves are generated based on P-wave time difference and density logging data.

[0016] Rock mechanical layers are defined based on rock mechanical index curves;

[0017] Wellbore calibration is performed based on the defined rock mechanical layers.

[0018] Preferably, the step of generating the rock mechanical index curve based on P-wave travel time and density logging data includes the following steps:

[0019] The longitudinal wave time difference is measured, the density value is obtained, the transverse wave time difference is calculated, and the rock mechanical index is calculated based on the longitudinal wave time difference and the transverse wave time difference.

[0020] More preferably, the method for obtaining the density value includes obtaining the density by consulting density logging data and calculating the density based on the P-wave time difference value;

[0021] The method of calculating density based on the P-wave time difference includes the following steps:

[0022] Based on the existing P-wave transit time and well logging density data in this region, the relationship between P-wave transit time and density is obtained. The P-wave transit time in areas with unknown density in this region is determined, and the density is calculated by substituting it into the relationship between P-wave transit time and density.

[0023] More preferably, the relationship between the longitudinal wave time difference and density is as follows:

[0024]

[0025] Where ρ is density, Δt p Let A and B be the longitudinal wave time difference values, where A and B are constants.

[0026] More preferably, when the unit of ρ is g / cm³ 3 And the Δt p When the unit is μs / ft, A is 0.9-0.95 and B is 0.101-0.12.

[0027] More preferably, and as a specific embodiment of the present invention, the unit of ρ is g / cm³. 3 And the Δt p When the unit is μs / ft, A is 0.9044 and B is 0.1101.

[0028] More preferably, the formula for calculating the transverse wave time difference is as follows:

[0029]

[0030] Where, Δt p The P-wave time difference is expressed in μs / ft, ρ b This is the density value, in g / cm³. 3 , Δt s This represents the transverse wave time difference, in μs / ft.

[0031] More preferably, the rock mechanical index is calculated as follows:

[0032]

[0033] Where K is the rock mechanics index, υ is Poisson's ratio, and E is Young's modulus;

[0034] More preferably, the method for calculating the Poisson's ratio is as follows:

[0035]

[0036] Where υ is Poisson's ratio, Δt s The longitudinal wave time difference is expressed in μs / ft, Δt. p This is the transverse wave time difference, in μs / ft;

[0037] More preferably, the method for calculating the Young's modulus is as follows:

[0038]

[0039] Where E is Young's modulus, Δt s The longitudinal wave time difference is expressed in μs / ft, Δt. p The transverse wave time difference is expressed in μs / ft, ρ. b This is the density value, in g / cm³. 3 .

[0040] Preferably, the process of dividing the rock into mechanical layers includes the following steps:

[0041] The rock mechanics layer is divided into three layers: A, B, and C. The spatial relationship between these three layers is as follows:

[0042] Rock mechanics layer A < Rock mechanics layer B < Rock mechanics layer C.

[0043] More preferably, the method for dividing rock mechanical layers A is as follows: the segments with a rock mechanical index greater than or equal to 0.43 and less than 0.5 are classified as Class 2 rock mechanical layers A, and the segments with a rock mechanical index greater than or equal to 0.5 are classified as Class 1 rock mechanical layers A.

[0044] More preferably, the method for dividing rock mechanical layer B is as follows: from the top to the bottom of the target layer, based on the existing surface geological stratification, a curve segment of 60-120m is defined, and a curve segment in which the overall curve value decreases or increases by 10%-20% is defined as a rock mechanical layer B.

[0045] More preferably, the method for dividing the rock mechanical layer C is as follows: from the top to the bottom of the target layer, 2-5 rock mechanical layers B are combined, and the combination of the overall curve value increases or decreases by 30%-45% is used as a rock mechanical layer C.

[0046] More preferably, in the combination of 2-5 rock mechanical layers B, the scale of the combination is 200-400m.

[0047] Preferably, the wellbore calibration includes the following steps:

[0048] The wavelet with the same frequency as the dominant frequency of the well perimeter seismic wave was selected for calibration. The calibration results were verified, and the rock mechanical layers of the single well were calibrated onto the profile. Three-dimensional tracking was performed to complete the division of the three-dimensional rock mechanical layers.

[0049] More preferably, the calibration includes the following steps: during calibration, the interface of the rock mechanical layer B corresponds to the seismic phase axis, the interface with the rock mechanical index increasing downward corresponds to the wave trough reflection, and the interface with the rock mechanical index decreasing downward corresponds to the wave crest reflection.

[0050] More preferably, the test calibration results include the following steps:

[0051] The well-seismic calibration results of the interface of the rock mechanical layer B are verified by the reflection coefficient sequence. The absolute value of the reflection coefficient corresponding to the interface should be the maximum value within 60m above and below it.

[0052] On the other hand, the present invention provides the application of the above method in the analysis of three-dimensional rock mechanical layers.

[0053] Furthermore, the present invention provides an apparatus for combining well-seismic analysis to delineate three-dimensional rock mechanical layers, used to implement the above-mentioned method, comprising:

[0054] Module 1 is used to generate rock mechanical index curves based on P-wave time difference and density logging data;

[0055] Module 2 is used to delineate rock mechanical layers based on rock mechanical index curves;

[0056] Module 3 is used for wellbore calibration based on the defined rock mechanical layers;

[0057] The module 1 transmits the generated rock mechanical index curve to the module 2 to delineate rock mechanical layers;

[0058] The module 2 transmits the obtained rock mechanical layers to the module 3 for well-seismic calibration.

[0059] In another aspect, the present invention provides an electronic device, comprising:

[0060] Memory, which stores executable instructions;

[0061] The processor executes the executable commands in the memory to implement the above-described method of well-seismic combination for dividing three-dimensional rock mechanical layers.

[0062] In another aspect, the present invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for dividing three-dimensional rock mechanical layers using a combination of well and seismic analysis.

[0063] A rock mechanical layer refers to a rock stratum with similar or identical rock mechanical behavior and properties. It is a rock mechanical unit defined by rock mechanical interfaces, possessing similar or identical strength, brittleness, and fracture mechanical properties. It comprises both rock mechanical units and rock mechanical interfaces. In ultra-deep carbonate rocks, the distribution characteristics of fractures at different scales are key factors determining reservoir formation and high-yield oil and gas enrichment. Fracture-crack extension is thus limited by rock mechanical interfaces. Therefore, rock mechanical layers control the formation and distribution of secondary fractures. Currently, commonly used methods for delineating rock mechanical layers are limited to the vertical delineation of rock mechanical interfaces, necessitating a spatial three-dimensional rock mechanical layer delineation method. This invention, through the calculation of single-well rock mechanical indices and well-seismic calibration, completes the delineation of small, medium, and large-scale rock mechanical layers in single wells and the three-dimensional medium and large-scale rock mechanical layers at the seismic scale; it solves the problem of three-dimensional rock mechanical layer delineation in the study area and can be used to analyze the controlling effect of rock mechanical layers on the formation and distribution of fractures. Attached Figure Description

[0064] Figure 1 This is a comparison diagram of the rock mechanical layer division between well SHB84X and well SHB81ZDY in Example 2.

[0065] Figure 2 This is a diagram showing the well seismic calibration results of the mesoscale rock mechanics layer in well SHB84X in Example 2.

[0066] Figure 3 This is a diagram showing the well seismic calibration results of the mesoscale rock mechanics layers in well SHB81ZDY in Example 2. Detailed Implementation

[0067] The present invention will be further described below by way of specific embodiments. The following non-limiting embodiments are intended to enable those skilled in the art to gain a more comprehensive understanding of the present invention, but do not limit the invention in any way. The following content is merely an exemplary description of the scope of protection of the present invention, and those skilled in the art can make various changes and modifications to the invention based on the disclosed content, which should also fall within the scope of protection of the present invention.

[0068] Example 1

[0069] A method for dividing three-dimensional rock mechanical layers by combining well and seismic analysis

[0070] The first step involved statistically analyzing the correlation between P-wave transit time and density curves from 15 wells with available density logging data in the Shunbei area. A formula for calculating density using P-wave transit time was then established using the Gardner formula method.

[0071]

[0072] In the formula, ρ is the density value, in g / cm³.3 Δtp is the longitudinal wave time difference, in μs / ft; A and B are constants, 0.9044 and 0.1101, respectively.

[0073] The transverse wave time difference can then be calculated from these two factors:

[0074]

[0075] In the formula, ρ b Δtp and Δts represent the density value, in g / cm3; Δtp and Δts represent the P-wave and S-wave time differences, respectively, in μs / ft.

[0076] Then, calculate Poisson's ratio (υ, dimensionless) and Young's modulus (E, GPa) using the following formulas:

[0077]

[0078] The strata have now undergone a fracturing process. Strata with high brittleness and high degree of fracturing often exhibit a lower Young's modulus and a higher Poisson's ratio. Therefore, the ratio of 100 times Poisson's ratio to Young's modulus is used as the rock mechanical index K to form the rock mechanical index curve.

[0079]

[0080] The higher the value, the more brittle and fragmented the strata.

[0081] The second step is to divide the rock mechanics layers into small-scale (meter-level), medium-scale (ten-meter-level), and large-scale (hundred-meter-level) layers. The small-scale rock mechanics layers control the development of fracture zones, the medium-scale rock mechanics layers control the development of fracture-cavity bodies, and the large-scale rock mechanics layers control the development of secondary faults and beaded structures.

[0082] First, small-scale rock mechanical layers are delineated. Specifically, layers with a rock mechanical index greater than or equal to 0.43 and less than 0.5 are classified as Class 2 small-scale rock mechanical layers, and layers with a rock mechanical index greater than or equal to 0.5 are classified as Class 1 small-scale rock mechanical layers. Then, medium-scale rock mechanical layers are delineated, analogous to stratigraphic division, using the half-amplitude point of the rock mechanical index curve where there is a large fluctuation as the interface. Finally, large-scale rock mechanical layers are delineated, with the interface being the interface of the higher-level rock mechanical index curve, which is a combination of multiple medium-scale rock mechanical layers with similar rock mechanical indices.

[0083] The third step is to perform fine-grained well-seismic calibration. Rock mechanical interfaces are those where wave impedance changes abruptly, corresponding to the seismic phase axis. A wavelet with the same frequency as the dominant frequency of the well-perimeter seismic profile is selected. During calibration, considering seismic resolution, the mesoscale rock mechanical interface should correspond to the seismic phase axis. Interfaces with decreasing rock mechanical indices indicate decreasing wave impedance below that interface, corresponding to trough reflections; interfaces with decreasing rock mechanical indices indicate increasing wave impedance below that interface, corresponding to crest reflections. Furthermore, the calibration results can be verified using a reflection coefficient sequence. The absolute value of the reflection coefficient corresponding to a rock mechanical layer should be the maximum value in its vicinity. After calibrating the single-well rock mechanical layers onto the seismic profile, these layers can be tracked in three dimensions, completing the three-dimensional delineation of the rock mechanical layers.

[0084] Example 2

[0085] According to the method described in Example 1, taking the Yijianfang Formation and Yingshan Formation of the Ordovician system in the Shunbei area as examples, the Shunbei 84 inclined well (SHB84X) and the Shunbei 81 inclined well straight pilot hole (SHB81ZDY), which have two boreholes penetrating the Yingshan Formation, were selected. Three-dimensional rock mechanics layer division was performed using the method described in the technical solution, combined with 3D seismic data. Seven and twelve small-scale rock mechanics layers were identified in the Shunbei 84 and Shunbei 81 inclined wells, respectively. Based on well correlation, a total of 14 medium-scale rock mechanics layers and 5 large-scale rock mechanics layers (e.g., ...) were identified. Figure 1 (As shown). Among them, the well seismic calibration results of the mesoscale rock mechanics layers in well SHB84X are shown in the figure. Figure 2 The well seismic calibration results of the mesoscale rock mechanics layers in well SHB81ZDY are shown in the figure. Figure 3 It is evident that the interfaces of the large- and medium-scale rock mechanical layers show good calibration results with the seismic phase axis, and can basically correspond to the maximum values ​​of the reflection coefficient (RC). Among them, the rock mechanical indices of the large-scale layers No. 4 and No. 5, located in the lower part of the Yingxia section, are significantly higher than those of the upper layers, indicating a higher degree of fragmentation.

[0086] Example 3

[0087] A device for dividing three-dimensional rock mechanical layers by combining well and seismic analysis can realize the method of dividing three-dimensional rock mechanical layers by combining well and seismic analysis in Examples 1 and 2.

[0088] This device includes the following modules:

[0089] Module 1 is used to generate rock mechanical index curves based on P-wave time difference and density logging data;

[0090] Module 2 is used to delineate rock mechanical layers based on rock mechanical index curves;

[0091] Module 3 is used for wellbore calibration based on the defined rock mechanical layers;

[0092] The module 1 transmits the generated rock mechanical index curve to the module 2 to delineate rock mechanical layers;

[0093] The module 2 transmits the obtained rock mechanical layers to the module 3 for well-seismic calibration.

[0094] Example 4

[0095] An electronic device includes: a memory storing executable instructions; and a processor that executes the executable instructions in the memory to implement the method for dividing three-dimensional rock mechanical layers using a combination of well and seismic analysis as described in this invention.

[0096] In the electronic device described in this embodiment, the memory is used to store non-transitory computer-readable instructions. Specifically, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.

[0097] In the electronic device described in this embodiment, the processor may be a central processing unit (CPU) or other form of processing unit with data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of this disclosure, the processor is used to execute computer-readable instructions stored in the memory.

[0098] Example 5

[0099] A computer-readable storage medium storing a computer program that, when executed by a processor, implements the method for well-seismic combined segmentation of three-dimensional rock mechanical layers as described in this invention.

[0100] The aforementioned computer-readable storage medium stores non-transitory computer-readable instructions. When these non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the methods described in the foregoing embodiments of this disclosure are performed.

[0101] The aforementioned computer-readable media include, but are not limited to:

[0102] Optical storage media, such as CD-ROM and DVD; magneto-optical storage media, such as MO; magnetic storage media, such as magnetic tape or portable hard drives; media with built-in rewritable non-volatile memory, such as memory cards; media with built-in ROM, such as ROM cartridges.

[0103] Finally, it should be noted that the above content is only used to illustrate the technical solution of the present invention, and is not intended to limit the scope of protection of the present invention. Simple modifications or equivalent substitutions made by those skilled in the art to the technical solution of the present invention do not depart from the essence and scope of the technical solution of the present invention.

Claims

1. A method for dividing three-dimensional rock mechanical layers by well-seismic combination, characterized in that, It comprises the following steps: Generating rock mechanics index curve according to the P-wave time difference and density logging data; Dividing rock mechanics layer according to the rock mechanics index curve; Calibrating well-seismic according to the divided rock mechanics layer.

2. The method of claim 1, wherein, The step of generating rock mechanics index curve according to the P-wave time difference and density logging data comprises the following steps: Measuring P-wave time difference value, obtaining density value, calculating S-wave time difference value, and calculating rock mechanics index according to the P-wave time difference value and the S-wave time difference value.

3. The method of claim 2, wherein, The method of obtaining density value comprises consulting density logging data to obtain density and calculating density according to P-wave time difference value; The step of calculating density according to P-wave time difference value comprises the following steps: Fitting P-wave time difference value and density data to obtain P-wave time difference value and density relationship formula, measuring P-wave time difference value of unknown density area in the region, substituting the P-wave time difference value into the P-wave time difference value and density relationship formula, and calculating to obtain density.

4. The method of claim 3, wherein, The P-wave time difference value and density relationship formula is as follows: wherein p is the density in g / cm 3 , At p is the difference in P-wave travel time in μs / ft, and A and B are constants.

5. The method of claim 2, wherein, The formula of calculating S-wave time difference value is as follows: where Δt p is the longitudinal time difference, in μs / ft, ρ b is the density value, in g / cm 3 , Δt s is the shear time difference, in μs / ft.

6. The method of claim 2, wherein, The calculation method of rock mechanics index is as follows: Wherein, K is rock mechanics index, υ is Poisson's ratio, and E is Young's modulus; The calculation method of Poisson's ratio is as follows: where υ is the Poisson's ratio, Δt s is the difference in longitudinal travel time in μs / ft, Δt p is the difference in shear travel time in μs / ft; The calculation method of Young's modulus is as follows: where E is Young's modulus, Δt s is the longitudinal time difference value in μs / ft, Δt p is the shear time difference value in μs / ft, ρ b is the density value in g / cm 3 .

7. The method of claim 1, wherein, The step of dividing rock mechanics layer comprises the following steps: Dividing rock mechanics layer A; dividing rock mechanics layer B; and dividing rock mechanics layer C, and the size relationship of the rock mechanics layer A, the rock mechanics layer B and the rock mechanics layer C in spatial range is as follows: Rock mechanics layer A < rock mechanics layer B < rock mechanics layer C.

8. The method of claim 8, wherein, The method of dividing rock mechanics layer A is that the layer segment with rock mechanics index greater than or equal to 0.43 and less than 0.5 is taken as 2-type rock mechanics layer A, and the layer segment with rock mechanics index greater than or equal to 0.5 is taken as 1-type rock mechanics layer A; The method of dividing rock mechanics layer B is that from the top to the bottom of the target layer, the existing well geological layering is taken as the basis, 60-120m is taken as a curve segment, and the curve segment with overall curve value changing by 10%-20% is taken as a rock mechanics layer B; The method of dividing rock mechanics layer C is that from the top to the bottom of the target layer, 2-5 rock mechanics layers B are taken as a combination, and the combination with overall curve value changing by 30%-45% is taken as a rock mechanics layer C.

9. The method of claim 1, wherein, The well-seismic calibration comprises the following steps: Selecting wavelet with consistent frequency with the seismic main frequency around the well to calibrate, testing the calibration result, calibrating single well rock mechanics layer to profile, performing three-dimensional tracking, and completing three-dimensional rock mechanics layer division.

10. The method of claim 9, wherein, The step of calibrating comprises the following steps: during calibration, the interface of the rock mechanics layer B corresponds to seismic phase axis, the interface with rock mechanics index increasing downward corresponds to wave trough reflection, and the interface with rock mechanics index decreasing downward corresponds to wave peak reflection.

11. The method of claim 9, wherein, The step of testing the calibration result comprises the following steps: The well-seismic calibration result of the interface of the rock mechanics layer B is tested by reflection coefficient sequence, and the absolute value of the reflection coefficient corresponding to the interface should be the maximum value within 60m above and below the interface.

12. The method of any one of claims 1-11, applied to analyzing three-dimensional rock mechanics layers.

13. An apparatus for 3D rock mechanical layer division by combining well and seismic data, for implementing the method of any one of claims 1-11, characterized in that, The method comprises: a module 1 for forming a rock mechanics index curve according to the P-wave time difference and density logging data; a module 2 for dividing rock mechanics layers according to the rock mechanics index curve; a module 3 for well-seismic calibration according to the divided rock mechanics layers; the module 1 transmits the rock mechanics index curve formed to the module 2 for dividing rock mechanics layers; the module 2 transmits the rock mechanics layers divided to the module 3 for well-seismic calibration.

14. An electronic device, comprising: The electronic device comprises: a memory storing executable instructions; a processor running the executable instructions in the memory to implement the method of any one of claims 1-11.

15. A computer-readable storage medium, characterized in that, The computer readable storage medium stores a computer program, which is executed by a processor to implement the method of any one of claims 1-11. The computer readable storage medium stores a computer program, which is executed by a processor to implement the method of any one of claims 1-11.

Images