Continental facies shale interlayer bedding seam activation capability evaluation method and device and medium

Abstract

The invention provides a continental facies shale interlayer bedding seam activation capability evaluation method and device and a medium, and belongs to the technical field of oil and gas reservoir development. The method comprises the following steps: establishing a continental facies shale three-dimensional cohesion unit method multi-cluster fracturing model, performing multi-factor simulation analysis by using the model to obtain fracturing crack propagation forms under different factors and activation conditions of a shale interlayer bedding interface, and drawing an interlayer bedding interface activation chart. According to the method, the three-dimensional cohesion unit method multi-cluster fracturing model is adopted to describe the three-dimensional multi-cluster fracturing fracture form of the continental facies shale, the fracture expansion form and the bedding interface opening condition can be accurately evaluated, and the model can consider the mutual disturbance effect among the multi-cluster fractures and the stress interference effect of the adjacent bedding interfaces; the activation condition of each bedding interface under the conditions of different clusters, crack initiation positions, displacements and liquid viscosities is obtained through variable analysis, and help is provided for the fracturing construction process, segment cluster division optimization and construction parameter optimization.

Description

Technical Field

[0001] This invention belongs to the field of oil and gas reservoir development technology, specifically relating to a method, device, and medium for evaluating the activation capacity of interlayer bedding fractures in continental shale. Background Technology

[0002] In the development of shale oil and gas, hydraulic fracturing technology is typically required to modify the reservoir, create volumetric fractures, form a network of high-conductivity fractures, and thus construct flow channels between oil and gas within the shale reservoir and the wellbore. Continental shale typically has well-developed bedding, which can be activated during hydraulic fracturing, thereby affecting the fracture propagation morphology.

[0003] To accurately obtain the mechanical characteristics of continental shale bedding and its influence on fracture morphology, Chinese patent publication CN108254262A proposes a method and device for predicting shear parameters of rock bedding fractures. This involves compressing core samples with different bedding fracture dip angles that meet preset conditions to obtain the compression strength; based on the compression strength, the shear parameters of the bedding fracture surface are determined. Chinese patent publication CN115203807A proposes a method for calculating the geometric parameters of shale bedding fracture propagation. This involves establishing a geometric, mechanical, and fracture propagation model of shale bedding, dynamically simulating fracture propagation during shale gas volumetric fracturing, and obtaining the calculation results of the geometric parameters of shale bedding fracture propagation. Chinese patent publication CN110348032A proposes a numerical simulation method for the height of hydraulic fractures in bedding-developed shale formations. This method establishes a fluid-structure interaction shale model based on geomechanical parameters and determines whether bedding is activated and the fracture height during fracture propagation through simulation. Chinese patent publication CN106545324A proposes a method for opening multiple horizontal bedding fractures in shale gas. The method involves activating the bedding fractures through acid treatment and low-volume, low-viscosity slickwater, and then using proppant to form complex fracture support, thereby providing flow conductivity.

[0004] The existing technologies mentioned above focus more on the strength of bedding fracture opening, how to open more bedding fractures, the morphology of the opened fractures, and their support and conductivity. However, for continental shale, there are significant vertical lithological differences, bedding development, and large differences in interlayer rock mechanical properties and in-situ stress. Whether bedding fractures can be opened during fracturing is affected by a combination of factors, including interface strength, in-situ in-situ stress, and stress interference from already opened fractures. Therefore, it is necessary to establish an evaluation method for bedding interface activation to accurately guide the optimization of fracturing technology and parameters for continental shale and improve the fracturing effect. Summary of the Invention

[0005] The purpose of this invention is to solve the problems existing in the prior art and to provide a method, device and medium for evaluating the activation capacity of interlayer bedding fractures in continental shale, which can help with fracturing construction technology, segment cluster division and optimization and construction parameter optimization.

[0006] This invention is achieved through the following technical solution:

[0007] The first aspect of this invention provides a method for evaluating the activation capacity of interlayer bedding fractures in continental shale. A three-dimensional cohesive unit method multi-cluster fracturing model of continental shale is established. The model is used to conduct multi-factor simulation analysis to obtain the fracturing fracture propagation morphology and the activation status of interlayer bedding interfaces under different factors, and an interlayer bedding interface activation chart is drawn.

[0008] A further improvement of the present invention is that:

[0009] The method includes the following steps:

[0010] Step 1: Obtain logging data for the target shale reservoir;

[0011] Step 2: Conduct core testing to obtain core mechanical parameters and strength parameters of bedding interfaces;

[0012] Step 3: Correct the logging data using the core mechanical parameters obtained from core testing;

[0013] Step 4: Establish a three-dimensional cohesive element method multi-cluster fracturing model for continental shale;

[0014] Step 5: Use the established model to conduct multi-factor simulation analysis to obtain the crack propagation morphology and bedding interface opening under different factors.

[0015] Step 6: Draw the activation template for the interlayer layering interface.

[0016] A further improvement of the present invention is that:

[0017] The logging data in step 1 includes: elastic modulus, Poisson's ratio, compressive strength, tensile strength, permeability, porosity, formation pressure, and triaxial stress.

[0018] A further improvement of the present invention is that:

[0019] Step 2 involves obtaining the core mechanical parameters, specifically including:

[0020] First, core samples were obtained from each sub-layer of the target shale reservoir, excluding interlayer bedding. Then, the core samples from each sub-layer were tested to obtain the rock mechanical parameters under formation pressure conditions.

[0021] The tests include: triaxial compression test, tensile strength test, shear test, porosity and permeability test, and Kaiser stress test.

[0022] A further improvement of the present invention is that:

[0023] Step 2 obtains the strength parameters of the bedding interface. Specific operations include:

[0024] First, core samples containing interlayer bedding between two adjacent layers in the target shale reservoir are obtained. The shear strength and tensile strength of the bedding interface under formation pressure are then tested to obtain the shear strength and tensile strength of the bedding interface, which are the strength parameters of the bedding interface.

[0025] A further improvement of the present invention is that:

[0026] Step 4 involves establishing a three-dimensional cohesive element method multi-cluster fracturing model for continental shale. Specific operations include:

[0027] First, based on the sub-layer division of the target shale reservoir section, a fracturing geometric model containing bedding interfaces is established;

[0028] Then, the corrected logging data and bedding interface strength parameters are input into the fracturing geometry model containing bedding interfaces to obtain the continental shale fracturing fracture propagation model, which is the three-dimensional cohesive element method multi-cluster fracturing model.

[0029] A further improvement of the present invention is that:

[0030] Step 5 utilizes the established model to conduct multi-factor simulation analysis, obtaining crack propagation morphology and bedding interface opening under different factors. Specific operations include:

[0031] Using at least one factor as a variable, the variable is input into the three-dimensional cohesive unit method multi-cluster fracturing model of continental shale established in step 4 for numerical simulation, so as to obtain the fracture propagation morphology and bedding interface opening under different factors.

[0032] A further improvement of the present invention is that:

[0033] The factors include the initiation layer, the number and location of initiation points, the viscosity of the liquid, and the displacement.

[0034] A second aspect of the present invention provides a device for evaluating the activation capacity of interlayer bedding fractures in continental shale, comprising:

[0035] The acquisition unit is used to acquire logging data of the target shale reservoir;

[0036] The testing unit is used to perform core testing to obtain core mechanical parameters and strength parameters of bedding interfaces;

[0037] The correction unit is used to correct well logging data using core mechanical parameters obtained from core testing;

[0038] Model building unit, used to establish a three-dimensional cohesive element method multi-cluster fracturing model of continental shale;

[0039] The simulation analysis unit is used to perform multi-factor simulation analysis using the established model to obtain the crack propagation morphology and bedding interface opening under different factors.

[0040] The drawing unit is used to draw the activation drawing of the interlayer layer interface.

[0041] A third aspect of the present invention provides a computer-readable storage medium storing at least one computer-executable program, which, when executed by the computer, causes the computer to perform the steps in the method for evaluating the activation capacity of interlayer bedding fractures in continental shale.

[0042] Compared with the prior art, the beneficial effects of the present invention are:

[0043] This invention corrects well logging data based on core testing results and obtains true interface parameters based on interface testing results. It employs a three-dimensional cohesive element method multi-cluster fracturing model to describe the morphology of three-dimensional multi-cluster fracturing fractures in continental shale, accurately evaluating fracture propagation morphology and bedding interface opening. The model considers the mutual disturbance between multiple fracture clusters and the stress interference effect of adjacent bedding interfaces. Through variable analysis, it obtains the activation status of each bedding interface under different clusters, fracture initiation locations, flow rates, and fluid viscosities, and plots bedding interface activation charts, providing assistance for fracturing construction technology, optimal segment cluster division, and optimized construction parameters.

[0044] This invention conducts experimental tests on continental shale rocks and bedding interfaces under formation pressure conditions to obtain actual parameters. When establishing a fracturing model, parameter values ​​are assigned based on the corrected well logging data, avoiding the uniform assignment of shale parameters and bedding interfaces in previous models, thus improving the accuracy of the analysis results.

[0045] This invention employs a three-dimensional cohesive element method for multi-cluster fracturing model analysis of fracture propagation morphology and bedding activation. It can take into account the mutual disturbance of multi-cluster fracture propagation in shale, as well as the stress interaction between adjacent bedding interfaces, greatly improving the accuracy of bedding interface activation analysis. Attached Figure Description

[0046] Figure 1 This is a flowchart of a method for evaluating the activation capacity of interlayer bedding fractures in continental shale according to the present invention;

[0047] Figure 2 This is a flowchart illustrating the establishment of a multi-cluster fracturing model for continental shale using the three-dimensional cohesive element method in an embodiment of the present invention.

[0048] Figure 3 This describes the propagation of hydraulic fracturing fractures and the activation of bedding in continental shale, as illustrated in this embodiment of the invention.

[0049] Figure 4a This is the activation diagram of the layer interface of layer ② and layer ③ in this embodiment of the invention;

[0050] Figure 4b This is the activation diagram of the layer interface of layer ③ and layer ④ in the embodiment of the present invention. Detailed Implementation

[0051] The present invention will now be described in further detail with reference to the accompanying drawings:

[0052] In fracturing operations, bedding interfaces in continental shale are easily activated, significantly impacting fracture height. Activation of bedding interfaces can lead to complex fracture morphologies and enhance fracturing effectiveness. Therefore, accurately assessing the activation status of bedding interfaces between different layers is crucial. Based on this, this invention proposes a method for evaluating the activation capacity of bedding fractures between continental shale layers. The method involves correcting logging data based on core testing results and obtaining true interface parameters based on interface testing results. A fracturing geometric model incorporating bedding interfaces is established. Corrected logging interpretation data is used to assign values ​​to shale matrix material parameters, and bedding interface strength parameter test results are used to assign values ​​to bedding interface parameters. A three-dimensional cohesive element method multi-cluster fracturing model of continental shale is established. The method analyzes the fracture propagation morphology and activation status of individual shale bedding interfaces under different wellbore trajectories, cluster divisions, fluid displacements, and performance conditions. Activation charts for each bedding interface are then generated, providing guidance for wellbore trajectory selection, fracturing techniques, optimal cluster division, and optimization of construction parameters.

[0053]

Example 1

[0054] This invention provides a method for evaluating the activation capacity of inter-layer bedding fractures in continental shale, such as... Figure 1 As shown, the specific steps include:

[0055] Step 1: Obtain logging data for the target shale reservoir;

[0056] Step 2: Conduct core testing to obtain core mechanical parameters and strength parameters of bedding interfaces;

[0057] Step 3: Correct the logging data using the core mechanical parameters obtained from core testing;

[0058] Step 4: Establish a three-dimensional cohesive element method multi-cluster fracturing model for continental shale;

[0059] Step 5: Use the established model to conduct multi-factor simulation analysis to obtain the crack propagation morphology and bedding interface opening under different factors.

[0060] Step 6: Draw the activation template for the interlayer layering interface.

[0061]

Example 2

[0062] Step 1 involves obtaining logging data from the target shale reservoir. Specific operations include:

[0063] Well logging data includes: rock mineral composition, elastic modulus, Poisson's ratio, compressive strength, tensile strength, permeability, porosity, formation pressure, and triaxial stress.

[0064] The above logging data can all be obtained using existing technical methods, and will not be described in detail here.

[0065]

Example 3

[0066] Step 2 involves core testing to obtain core mechanical parameters and strength parameters of bedding interfaces. Specific operations include:

[0067] During core testing, the core samples must meet the following requirements:

[0068] Select a complete full-diameter core sample from this well or a pilot well on the same platform, which should include all layers of the target stratum in the vertical direction.

[0069] The operations for obtaining core mechanical parameters include:

[0070] First, core samples were obtained from each sub-layer of the target shale reservoir, excluding interlayer bedding. Then, the core samples from each sub-layer were tested to obtain the rock mechanical parameters under formation pressure conditions.

[0071] The tests include triaxial compression testing, tensile strength testing, shear testing, porosity and permeability testing, and Kaiser stress testing, all of which are existing technical testing methods and will not be described in detail here.

[0072] The operations for obtaining the strength parameters of the bedding interface include:

[0073] First, core samples containing interlayer bedding between two adjacent layers in the target shale reservoir were obtained. The shear strength and tensile strength of the bedding interface under formation pressure were then tested to obtain the shear strength and tensile strength of the bedding interface, i.e., the strength parameters of the bedding interface.

[0074]

Example 4

[0075] Step 3 involves using the core mechanical parameters obtained from core testing to correct the logging data. Specific operations include:

[0076] By using core test parameters from the same depth to correct the logging data, the accuracy of the analysis results was improved.

[0077]

Example 5

[0078] In step 4, a three-dimensional cohesive element method multi-cluster fracturing model of continental shale is established, such as... Figure 2 As shown, the specific operations include:

[0079] First, based on the division of each sub-layer in the target shale reservoir section (specifically referring to the height of each sub-layer), a fracturing geometric model containing the bedding interface is established.

[0080] Then, material parameters are assigned to the shale matrix based on the corrected logging data, and parameters are assigned to the bedding interfaces based on the strength parameters of the bedding interfaces. In other words, the corrected logging data and the strength parameters of the bedding interfaces are respectively input into the fracturing geometric model containing the bedding interfaces to obtain the continental shale fracturing fracture propagation model, which is the three-dimensional cohesive element method multi-cluster fracturing model.

[0081] In this embodiment of the invention, numerical simulation software such as ABAQUS, MATLAB, and COMSOL are used to establish a three-dimensional cohesive element method multi-cluster fracturing model.

[0082]

Example 6

[0083] Step 5 utilizes the established model to conduct multi-factor simulation analysis, obtaining crack propagation morphology and bedding interface opening under different factors. Specific operations include:

[0084] Using at least one factor as a variable, the variable is input into the three-dimensional cohesive element method multi-cluster fracturing model of continental shale established in step 4 for numerical simulation to obtain the fracture propagation morphology and bedding interface opening under different factors. The fracture propagation morphology includes fracture length and fracture height.

[0085] The factors mentioned include the initiation layer, the number and location of initiation points, the viscosity of the liquid, and the discharge rate.

[0086] For example, the liquid viscosity parameters are set to 3, 6, 10, 15, and 20 mPa·s, respectively, and then input into the three-dimensional cohesive unit method multi-cluster fracturing model of continental shale established in step 4 for numerical simulation to obtain the fracture propagation morphology and bedding interface opening under different liquid viscosities.

[0087] Set the liquid displacement parameters to 5, 10, 15, and 20 m respectively. 3 The value is / min, which is input into the three-dimensional cohesive unit method multi-cluster fracturing model of continental shale established in step 4 for numerical simulation to obtain the fracture propagation morphology and bedding interface opening under different fluid displacements.

[0088]

Example 7

[0089] Step 6 involves drawing the interlayer layering interface activation template. Specific operations include:

[0090] Based on the numerical simulation results in step 6, the opening status of each bedding interface under different factors is obtained, and the activation chart of each bedding interface is drawn to obtain its activation capability during the fracturing operation.

[0091]

Example 8

[0092] The method of the present invention will be further explained and illustrated below through a specific embodiment.

[0093] A horizontal well was sidetracked in a continental shale reservoir with a vertical depth of 2600m. The well was completed with casing, and the horizontal section was 1800m long. The target shale layer was divided into five sub-layers with a total height of 28m. Sub-layers ① to ⑤ had heights of 7m, 5m, 8m, 5m, and 3m, respectively. There were certain differences in lithology and mechanical parameters between the sub-layers, and the bedding interfaces were clearly developed. A full-diameter shale core of the reservoir section was obtained through coring in the pilot well. The activation capacity of the bedding interfaces in this continental shale reservoir was evaluated using the method of this invention. The specific implementation steps are as follows:

[0094] (1) Based on the geological parameters, logging data, and core test results of this well, the geological and mechanical characteristics of this well were interpreted. The stress distribution (vertical stress / maximum horizontal stress / minimum horizontal stress) is approximately 63 / 62 / 57 MPa, the formation pressure coefficient is 1.32, the elastic modulus is 24–28 GPa, the Poisson's ratio is 0.22–0.24, and the rock tensile strength is 5.1–6.7 MPa. From top to bottom, there are four bedding interfaces from layer ① to layer ⑤, with compressive strengths of 3 MPa, 2.4 MPa, 2.8 MPa, and 3.5 MPa, and shear strengths of 5 MPa, 4.3 MPa, 4.6 MPa, and 5.3 MPa, respectively.

[0095] (2) According to the fracturing design of this well, the experienced trajectory of this well traverses 100% of the ③ sub-layer, with a total of 23 fracturing operations. Each section has 3-5 perforation clusters, and each cluster has 16-18 spiral perforations with an included angle of 60° and a diameter of 9.5mm. The main fracturing flow rate is 18-21m³ / h. 3 / min, single-stage liquid volume approximately 2220~2300m³ 3 The overall sand ratio is 7.3%. The viscosity of low-viscosity slickwater is 3-5 mPa·s, the viscosity of high-viscosity slickwater is 9-12 mPa·s, and the viscosity of ultra-high-viscosity slickwater is 25-30 mPa·s.

[0096] (3) Taking the fracturing design of the 8th stage as an example, the perforation in this stage consists of 5 clusters with a cluster spacing of 10m, and the scale of low-viscosity slickwater used is 1750m. 3 The high-viscosity slickwater has a scale of 300m. 3 The ultra-high viscosity slippery water volume is 200m. 3 The 70 / 140 mesh proppant has a particle size of 90m. 340 / 70 mesh proppant has a size of 60m 3 30 / 50 mesh proppant has a size of 20m 3 .

[0097] (4) Establish a three-dimensional cohesive element method multi-cluster fracturing model for the continental shale in this well. This embodiment is based on the COMSOL platform, establishing a three-dimensional finite element multi-cluster fracturing model. Rock mechanical parameters, geostress parameters, and interface strength parameters are input into the model, and calculations are performed to obtain the fracturing fracture morphology and bedding interface opening morphology, such as... Figure 3 As shown.

[0098] (5) Conduct multi-factor variable simulation analysis to determine the opening status of the bedding interface under different fracturing parameter conditions. Based on the fracturing design parameters, this embodiment mainly analyzes the variables of fracturing flow rate and fracturing fluid viscosity. For example, according to the simulation results, for the bedding interface between layer ② and layer ③, when the fracturing fluid viscosity is 3 mPa·s, when the fracturing flow rate reaches 8 mPa·s... 3 At a flow rate of 10 m³ / min, this bedding interface can be activated; for the bedding interface between layers ③ and ④, when the fracturing flow rate reaches 10 m³ / min, the interface can be activated. 3 The layer interface can be activated when / min.

[0099] (6) Draw the stratification interface activation diagram. Simulations were conducted with displacement and viscosity as variables. The results show that the stratification interface between layers ② and ③ is the easiest to open, while the stratification interface between layers ③ and ④ requires a larger displacement to open. The stratification interfaces between layers ① and ②, and between layers ④ and ⑤, are more difficult to open. The stratification interface activation diagram is shown below. Figure 4a and Figure 4b As shown.

[0100] (7) Optimization of fracturing construction parameters. According to the simulation results, when using fracturing fluid with a viscosity of 3 mPa·s, the flow rate is gradually increased from small to large. When the fracturing fracture propagates, the bedding interface between layer ② and layer ③ is opened first, and then the bedding interface between layer ③ and layer ④ is opened. Because the bedding interface is activated first, the fracture propagates along the bedding, but the fracture fails to penetrate the layers in the longitudinal direction, and the overall height of the fracture is only 8m. When using fracturing fluid with a viscosity of 15 mPa·s, the flow rate is gradually increased from small to large. When the fracturing fracture propagates, the bedding interface between layer ② and layer ③ and between layer ③ and layer ④ is penetrated first, and then the bedding interface is opened, and the overall height of the fracture reaches 18m. When the viscosity reaches 25 mPa·s, the overall height of the fracture reaches 25m.

[0101] (8) Optimize construction parameters based on simulation results. The simulation results show that using fracturing fluid with a viscosity of 25 mPa·s can enable the fracture to penetrate the bedding interface first and achieve longitudinal expansion of the fracture. Then, using fracturing fluid with a viscosity of 3 mPa·s can activate the bedding interface between layer ② and layer ③ and between layer ③ and layer ④, thus increasing the complexity of the fracture to a limited extent.

[0102] (9) Repeat steps (3) to (8) for other fracturing sections to simulate and optimize the fracture parameters of all fracturing sections until a more ideal fracture morphology is obtained.

[0103]

Example 9

[0104] This invention provides a device for evaluating the activation capacity of inter-layer bedding fractures in continental shale, the device comprising:

[0105] The acquisition unit is used to acquire logging data of the target shale reservoir;

[0106] The testing unit is used to perform core testing to obtain core mechanical parameters and strength parameters of bedding interfaces;

[0107] The correction unit is used to correct well logging data using core mechanical parameters obtained from core testing;

[0108] Model building unit, used to establish a three-dimensional cohesive element method multi-cluster fracturing model of continental shale;

[0109] The simulation analysis unit is used to perform multi-factor simulation analysis using the established model to obtain the crack propagation morphology and bedding interface opening under different factors.

[0110] The drawing unit is used to draw the activation drawing of the interlayer layer interface.

[0111] This invention conducts experimental tests on continental shale rocks and bedding interfaces under formation pressure conditions to obtain actual parameters. When establishing a fracturing model, parameter values ​​are assigned based on the corrected well logging data, avoiding the uniform assignment of shale parameters and bedding interfaces in previous models, thus improving the accuracy of the analysis results.

[0112] This invention employs a three-dimensional cohesive element method for multi-cluster fracturing model analysis of fracture propagation morphology and bedding activation. It can take into account the mutual disturbance of multi-cluster fracture propagation in shale, as well as the stress interaction between adjacent bedding interfaces, greatly improving the accuracy of bedding interface activation analysis.

[0113]

Example 10

[0114] This invention also provides a computer-readable storage medium storing at least one computer-executable program, which, when executed by the computer, causes the computer to perform the steps in the method for evaluating the activation capacity of interlayer bedding fractures in continental shale.

[0115] Those skilled in the art will understand that all or part of the processes in the above embodiments can be implemented by a computer program instructing related hardware. The program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments described above. Any references to memory, storage, databases, or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), RAMbus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and RAMbus dynamic RAM (RDRAM), etc.

[0116] The above technical solution is only one embodiment of the present invention. For those skilled in the art, based on the principles disclosed in the present invention, it is easy to make various types of improvements or modifications, and not limited to the technical solutions described in the specific embodiments of the present invention. Therefore, the foregoing description is only a preferred option and is not restrictive.

Claims

1. A method of evaluating interlayered bed fracture activation capacity in a continental shale, characterized by, A three-dimensional cohesive zone element method multi-cluster fracturing model of continental shale is established, and the fracture propagation morphology and the activation of interlayer interface under different factors are obtained by using the model for multi-factor simulation analysis, and an interlayer interface activation chart is drawn.

2. The method of claim 1, wherein, The method comprises the following steps: Step 1, obtaining the logging data of the target shale reservoir; Step 2, performing core testing to obtain core mechanical parameters and interlayer interface strength parameters; Step 3, correcting the logging data by using the core mechanical parameters obtained by core testing; Step 4, establishing a three-dimensional cohesive zone element method multi-cluster fracturing model of continental shale; Step 5, using the established model to perform multi-factor simulation analysis to obtain the fracture propagation morphology and the opening of the interlayer interface under different factors; Step 6, drawing an interlayer interface activation chart.

3. The method of claim 2, wherein, The logging data in step 1 includes elastic modulus, Poisson's ratio, compressive strength, tensile strength, permeability, porosity, formation pressure and three-dimensional ground stress.

4. The method of claim 2, wherein, In step 2, the core mechanical parameters are obtained, and the specific operation includes: First, obtain the cores in each small layer of the target shale reservoir without interlayer, and then test each small layer core to obtain the rock mechanical parameters under the condition of formation pressure; The test content includes triaxial compression test, tensile strength test, shear test, porosity, permeability test and Kaiser stress test.

5. The method of claim 2, wherein, In step 2, the strength parameters of the interlayer interface are obtained, and the specific operation includes: First, obtain the cores containing interlayer between adjacent two layers in the target shale reservoir, and test the shear strength and tensile strength of the interlayer interface under the condition of formation pressure to obtain the shear strength and tensile strength of the interlayer interface, which are the strength parameters of the interlayer interface.

6. The method of claim 2, wherein, In step 4, the three-dimensional cohesive zone element method multi-cluster fracturing model of continental shale is established, and the specific operation includes: First, according to the small layer division of the target shale reservoir section, a fracturing geometric model containing interlayer interface is established; Then, the corrected logging data and the interlayer interface strength parameters are input into the fracturing geometric model containing the interlayer interface to obtain the fracturing fracture propagation model of continental shale, which is the three-dimensional cohesive zone element method multi-cluster fracturing model.

7. The method of claim 2, wherein, In step 5, the established model is used for multi-factor simulation analysis to obtain the fracture propagation morphology and the opening of the interlayer interface under different factors, and the specific operation includes: Taking at least one factor as a variable, inputting the variable into the three-dimensional cohesive zone element method multi-cluster fracturing model of continental shale established in step 4 to perform numerical simulation to obtain the fracture propagation morphology and the opening of the interlayer interface under different factors.

8. The method of claim 7, wherein, The factors include the fracturing layer, the number and position of the fracturing points, the liquid viscosity and the displacement.

9. A device for evaluating the capacity of interlayered bedded seams in a continental shale to be activated, characterized in that, It comprises: An acquisition unit is configured to acquire logging data of a target shale reservoir; A testing unit is configured to perform core testing to obtain core mechanical parameters and interlayer interface strength parameters; A correction unit is configured to correct the logging data by using the core mechanical parameters obtained by core testing; A model construction unit is configured to establish a three-dimensional cohesive zone element method multi-cluster fracturing model of continental shale; A model construction unit is configured to establish a three-dimensional cohesive zone element method multi-cluster fracturing model of continental shale; The simulation analysis unit is configured to perform multi-factor simulation analysis by using the established model to obtain the crack propagation pattern and the bedding interface opening condition under different factors. The chart drawing unit is configured to draw an interlayer bedding interface activation chart.

10. A computer-readable storage medium, characterized in that, The computer readable storage medium stores at least one program executable by the computer, and the at least one program, when executed by the computer, causes the computer to perform the steps in the continental shale interlayer bedding fracture activation capacity evaluation method according to any one of claims 1-8.

Images