Multi-cluster fracturing crack propagation simulation method considering natural crack interference

Abstract

The invention discloses a multi-cluster fracturing fracture propagation simulation method considering natural fracture interference. The method comprises the steps that reservoir geology, physical property and fracturing construction parameters are collected; constructing a two-dimensional geologic model representing the spatial distribution characteristics of the natural fractures; based on a finite element and global cohesion unit technology, establishing a fluid-solid coupling model of segmented multi-cluster fracturing; a pipe flow unit is introduced, non-uniform distribution of inter-cluster fluid in the simultaneous crack expansion process is simulated, and a dynamic flow distribution model containing shaft friction resistance and perforation hole friction resistance is established; parameter sensitivity analysis is carried out, and the influence of perforation parameters, ground stress difference, hydraulic load and other factors on crack propagation behaviors is studied; according to the method, key factors such as natural fracture interference, inter-fracture stress coupling and flow competition are fully considered, the fracture expansion process and complex fluid distribution characteristics can be effectively reflected, and theoretical support and technical reference are provided for fracture design and optimization; the invention further comprises a system, equipment and a storage medium for implementing the method.

Description

Technical Field

[0001] This invention relates to the field of unconventional oil and gas reservoir development technology, and in particular to a method for simulating the propagation of multi-cluster fracturing fractures that takes into account the interference of natural fractures. Background Technology

[0002] With the continuous development of unconventional oil and gas resources, deep gas reservoirs have become an important target for exploration and development due to their abundant reserves and huge development potential. For these reservoirs, horizontal well segmented multi-cluster fracturing technology has become a core engineering means to improve single-well production and effectively modify complex reservoir structures.

[0003] However, reservoirs generally exhibit characteristics such as well-developed natural fractures, large differences in geostress, and strong heterogeneity. This makes hydraulic fractures susceptible to interference from natural fractures and stress coupling between fractures during propagation, leading to complex fracture morphology evolution and uneven flow distribution among clusters. Consequently, this affects the fracturing effect and final production capacity. Traditional fracturing design methods are mostly based on ideal condition assumptions, neglecting the mutual interference and flow competition between fractures in multi-cluster fracturing, and thus failing to accurately reflect the dynamic propagation behavior of fractures.

[0004] Furthermore, existing numerical modeling studies on the simultaneous propagation of multiple fracture clusters generally lack a systematic integration of key factors such as the effects of natural fractures, stress disturbance mechanisms, and non-uniform flow distribution. For example, Tang Xuanhe's paper "A Complex Fracture Cross-Propagation Model for Hydraulic Fracturing in Shale Gas Reservoirs Based on FEM-DFN," published in *Natural Gas Industry*, 2023, 43(1):162-176, analyzes the propagation law of multiple hydraulic fracture clusters in a natural fracture network, but does not consider the influence of flow competition between clusters. Such methods are insufficient for conducting sensitivity analysis of fracture behavior and forming regular understandings that can be used for engineering decision-making. Therefore, it is urgent to establish a fracture propagation simulation method that can comprehensively consider the above-mentioned multiple influencing factors, in order to deeply explore the mechanism of action of fracturing parameters on fracture propagation behavior, and provide theoretical support and engineering basis for unconventional reservoir fracturing construction. Summary of the Invention

[0005] In order to overcome the shortcomings of the prior art, the present invention aims to provide a multi-cluster fracturing fracture propagation simulation method that considers the interference of natural fractures. This method can comprehensively consider the influence of natural fracture distribution, inter-fracture stress interference, and inter-cluster flow competition, and is used to study the influence of fracturing construction parameters on fracture propagation behavior.

[0006] To achieve the above objectives, the technical solution of the present invention is as follows: A method for simulating the propagation of multi-cluster hydraulic fracturing fractures that takes into account the interference of natural fractures includes the following steps: Step 1: Obtain the geological parameters, construction parameters, and completion parameters of the reservoir in the fracturing block; Step 2: Based on the geological parameters, reservoir structural characteristics, and geological statistical laws from Step 1, establish a reservoir natural fracture model that reflects the spatial distribution characteristics of natural fractures. Step 3: Based on the reservoir natural fracture model in Step 3, the finite element method and the method of inserting global cohesive elements are used to establish the rock constitutive equation, fluid-structure interaction equation, fluid flow equation within the fracture, fracture propagation criterion equation, cohesive element damage equation, and stiffness degradation criterion equation for the stress field and flow field of the segmented multi-cluster fracturing in a two-dimensional horizontal well. These equations are used to describe the evolution process of the fracture and the rock damage behavior. Step 4: Based on Step 3, a pipe flow unit is introduced to simulate the non-uniform distribution behavior of fluid flow between clusters under the condition of simultaneous expansion of multiple fractures. A dynamic flow distribution model considering fracture competition effect and inter-cluster interference is established. The dynamic flow distribution model includes the energy conservation equation of the fluid system, the wellbore friction pressure drop equation, and the perforation hole friction pressure drop equation. Step 5: Based on the model constructed in Step 4, set multiple sets of parameter conditions to obtain the influence of different perforation parameters, ground stress difference, and hydraulic load on the crack propagation behavior and complexity. Step 6: Visualize the results, which include fracture propagation morphology, fluid distribution patterns, and the influence of various fracturing parameters on fracture behavior.

[0007] The geological parameters mentioned in step one include natural fracture angle, natural fracture length, Young's modulus of rock, Poisson's ratio of rock, shear modulus of rock, and filtration coefficient; Construction parameters include fracturing fluid density, fracturing fluid viscosity, and injection rate. Completion parameters include the number of perforations and the perforation diameter.

[0008] The establishment of the reservoir natural fracture model in step two is specifically as follows: The starting point of the crack ( , The fractures are randomly distributed within the reservoir to ensure minimal distance between them. To avoid overlap, the starting position of the crack is randomly generated using the following formula: in, The length of the reservoir's edge; the Random() function is used to generate random numbers within a specified range; the length of the fracture. and angle By specifying the minimum length and maximum length and minimum angle and maximum angle To generate randomly: The other end of the crack ( , Based on the starting point ( , ), length of the crack and angle Perform the calculation: Step three specifically involves: 3.1 Rock constitutive equation: The linear elastic constitutive equation of rock is shown below: in, It is Poisson's ratio (dimensionless). It is the shear modulus. ; , For Young's modulus, , where is the Biot coefficient (dimensionless). Pore ​​pressure, Constitutive equations based on stress and solid strain express; 3.2. Establish fluid-structure interaction equations, including the coupling control equations for rock solid skeleton deformation and fluid flow, the mass conservation equation for fluid seepage, and the fluid flow velocity equation within the rock, expressed as follows: in, For effective stress, It is the identity matrix; For virtual strain rate, It is a surface force vector per unit area. , It is the volume force per unit volume. ; It is a virtual velocity field. ; per unit volume ; Represents the area of ​​the unit subjected to surface forces. ; It represents the volume change of fluid in a porous medium (dimensionless). For fluid density, ; The porosity is related to the porosity of porous media; For spatial vectors, ; The seepage velocity of the fluid. ; For time; It is the acceleration due to gravity. , The penetration matrix, ; 3.3 Establish the fluid flow equations within the fracture, including the tangential fluid flow equation, the fluid mass conservation equation, and the fracturing fluid loss equation, as specifically expressed below: In the above formula, This is tangential flow within the crack. , The pressure gradient is along the length of the element. , The viscosity of the fracturing fluid. , The thickness of the crack propagation. , , The volumetric flow rate per unit time is the volumetric flow rate of the upper and lower surfaces of the crack. ; , The filtration coefficients are the upper and lower surface filtration coefficients. , The pore pressure at the upper and lower surfaces of the crack. ; The fluid pressure of the unit. ; 3.4. Establish the crack propagation criterion equation, the cohesive element damage equation, and the stiffness degradation criterion equation for the damage evolution of the cohesive element, as detailed below: In the above formula, , , These represent the critical normal stress and the critical tangential stresses in two directions at the point of cohesive unit failure. , , It is the normal stress component before the damage. It is the tangential stress component before the damage; This indicates the overall degree of damage to the cohesive unit, which increases linearly from 0 to 1 after the element damage begins. This represents the maximum displacement amplitude reached by the element during the loading process; This represents the displacement amplitude when the element is completely destroyed. This represents the displacement amplitude at the initial damage point of the element.

[0009] In step four, the established flow dynamic distribution model specifically includes the fluid energy conservation equation, the wellbore friction pressure drop equation, and the perforation orifice friction pressure drop equation, as shown below: In the above formula, It is the injection rate. Is entering the cluster The flow rate, It is the number of clusters. and These represent the pressure drop as fluid passes through the wellbore element and the vertical distance between the nodes at both ends of the wellbore element, respectively. and These represent the flow velocity and density of the fracturing fluid in the wellbore, respectively. and These represent the acceleration due to gravity and the coefficient of friction of the wellbore, respectively. This refers to the direction-dependent drag coefficient in the wellbore. and These are the pressures of the fracturing fluid before and after it flows through the perforation element. This indicates the number of perforations in each cluster; The diameter of the perforation hole; It is a dimensionless flow coefficient used to inversely map the shape changes of orifices caused by wear.

[0010] The present invention also includes: A system comprising a parameter acquisition module, a modeling and simulation module, and a display module, capable of simulating the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference, wherein the parameter acquisition module implements step one of the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference; the modeling and simulation module implements steps two through five of the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference; and the display module implements step six of the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference.

[0011] An apparatus comprising: Memory: Used to store the computer program for the simulation method of multi-cluster fracturing fracture propagation considering the interference of natural fractures; Processor: Used to implement the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference when executing the computer program.

[0012] A computer-readable storage medium storing a computer program that, when executed by a processor, implements the aforementioned method for simulating multi-cluster fracturing fracture propagation considering interference from natural fractures.

[0013] Compared with the prior art, the beneficial effects of the present invention are: 1. Step two of this invention is based on the geological parameters of step one, and introduces the distribution characteristics of natural fractures and reservoir structural characteristics. It uses a random generation method to simulate the spatial distribution of natural fractures to ensure that it can reflect the distribution of fractures in the actual reservoir.

[0014] 2. Step four of this invention establishes a more reasonable dynamic distribution model of flow by considering the influence of flow competition between clusters. This model can truly reflect the mutual influence between fractures in multi-cluster fracturing, making the simulation results more consistent with actual working conditions.

[0015] 3. The combination of steps four and five in this invention enables the analysis of the influence of key factors such as different perforation parameters and geostress differences on fracture propagation under multi-cluster fracturing conditions, and has the advantages of flow distribution prediction and fracture propagation prediction.

[0016] 4. Step six of this invention, based on steps one to five, visualizes the crack propagation morphology, fluid distribution patterns, and the influence of fracturing parameters on crack behavior, providing an intuitive and effective basis for fracturing optimization.

[0017] In summary, this invention, based on the finite element method and fluid-structure interaction model, combined with the dynamic fluid distribution characteristics under multi-cluster fracturing, can simulate the non-uniform flow and competition effect of fluid between fractures, providing theoretical support for reservoir fracturing optimization. This method can analyze the influence of different fracturing construction parameters on fracture propagation behavior, which helps to optimize construction parameters and thus improve fracturing effect. Attached Figure Description

[0018] Figure 1 These are the implementation steps of the present invention.

[0019] Figure 2 This is a model of natural fractures in the reservoir.

[0020] Figure 3 A dynamic flow distribution model considering natural cracks and inter-cluster interference is proposed.

[0021] Figure 4 The flow distribution diagram for the three fracture clusters is shown under different injection rates.

[0022] Figure 5 The flow distribution diagram for the three clusters of fractures under different perforation diameters is shown.

[0023] Figure 6 This is a diagram showing the crack morphology. Detailed Implementation

[0024] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the embodiments and accompanying drawings.

[0025] Reference Figure 1A method for simulating the propagation of multi-cluster hydraulic fracturing fractures that takes into account the interference of natural fractures includes the following steps: Step 1: Obtain the geological parameters, construction parameters, and completion parameters of the reservoir in the fracturing block. The geological parameters include natural fracture angle, natural fracture length, Young's modulus of the rock, Poisson's ratio of the rock, shear modulus of the rock, and filtration coefficient. Construction parameters include fracturing fluid density, fracturing fluid viscosity, and injection rate. Completion parameters include the number of perforations and the perforation diameter.

[0026] As shown in Table 1; Table 1 Step 2: Using the natural fracture angles and lengths obtained in Step 1, generate natural fractures in the reservoir using formulas (1) to (3), and establish a natural fracture model of the reservoir, such as... Figure 2 As shown, the influence of natural cracks on crack propagation was considered, providing a basis for subsequent crack propagation simulation.

[0027] The starting point of the crack ( , The fractures are randomly distributed within the reservoir to ensure minimal distance between them. To avoid overlap, the starting position of the crack is randomly generated using the following formula: in, It is the side length of the reservoir. The Random() function is used to generate random numbers within a specified range.

[0028] Length of the crack and angle By specifying the minimum length and maximum length and minimum angle and maximum angle To generate randomly: The other end of the crack ( , Based on the starting point ( , ), length of the crack and angle Perform the calculation: Step 3: Based on Step 3, a pipe flow unit is introduced to simulate the non-uniform distribution of fluid flow between clusters under the condition of simultaneous propagation of multiple fractures. Based on formulas (4) to (7), fluid-structure interaction equations of stress field and flow field in two-dimensional horizontal well segmented multi-cluster fracturing, fluid flow equations within fractures, fracture propagation criterion equations, rock constitutive equations, cohesive unit damage equations, and stiffness degradation criterion equations for cohesive unit damage evolution are established to describe the evolution process of fractures and rock damage behavior. 3.1 Rock constitutive equations: Using a linear elastic constitutive model as the constitutive relation for rocks can satisfy the computational accuracy requirements. The linear elastic constitutive equation for rocks is shown below: in, It is Poisson's ratio (dimensionless). It is the shear modulus. ; , For Young's modulus, , where is the Biot coefficient (dimensionless). Pore ​​pressure, Constitutive equations based on stress and solid strain express.

[0029] 3.2. Establish fluid-structure interaction equations, including the coupling control equations for rock solid skeleton deformation and fluid flow, the mass conservation equation for fluid seepage, and the fluid flow velocity equation within the rock, expressed as follows: in, For effective stress, It is the identity matrix; For virtual strain rate, It is a surface force vector per unit area. , It is the volume force per unit volume. ; It is a virtual velocity field. ; per unit volume ; Represents the area of ​​the unit subjected to surface forces. ; It represents the volume change of fluid in a porous medium (dimensionless). For fluid density, ; The porosity is related to the porosity of porous media; For spatial vectors, ; The seepage velocity of the fluid. ; For time; It is the acceleration due to gravity. , The penetration matrix, .

[0030] 3.3 Establish the fluid flow equations within the fracture, including the tangential fluid flow equation, the fluid mass conservation equation, and the fracturing fluid loss equation, as specifically expressed below: In the above formula, This is tangential flow within the crack. , The pressure gradient is along the length of the element. , The viscosity of the fracturing fluid. , The thickness of the crack propagation. , , The volumetric flow rate per unit time is the volumetric flow rate of the upper and lower surfaces of the crack. ; , The filtration coefficients are the upper and lower surface filtration coefficients. , The pore pressure at the upper and lower surfaces of the crack. ; The fluid pressure of the unit. .

[0031] 3.4. Establish the crack propagation criterion equation, the cohesive element damage equation, and the stiffness degradation criterion equation for the damage evolution of the cohesive element, as detailed below: In the above formula, , , These represent the critical normal stress and the critical tangential stresses in two directions at the point of cohesive unit failure. , , It is the normal stress component before the damage. It is the tangential stress component before the damage. This indicates the overall degree of damage to the cohesive unit, which increases linearly from 0 to 1 after the element damage begins. This represents the maximum displacement amplitude reached by the element during the loading process; This represents the displacement amplitude when the element is completely destroyed. This represents the displacement amplitude at the initial damage point of the element.

[0032] Step 4: Based on formula (8), introduce a pipe flow unit to simulate the non-uniform distribution behavior of fluid flow between clusters under the condition of simultaneous propagation of multiple fractures. Establish a dynamic flow distribution model considering fracture competition effect and inter-cluster interference. The dynamic flow distribution model includes the energy conservation equation of the fluid system, the wellbore friction pressure drop equation, and the perforation hole friction pressure drop equation; such as Figure 3 As shown, this illustrates that the non-uniform distribution of fluid flow during multi-crack propagation is caused by the mutual influence and competition between cracks.

[0033] The fluid energy conservation equation, the wellbore friction pressure drop equation, and the perforation hole friction pressure drop equation are shown below: In the above formula, It is the injection rate. Is entering the cluster The flow rate, It is the number of clusters. and These represent the pressure drop as fluid passes through the wellbore element and the vertical distance between the nodes at both ends of the wellbore element, respectively. and These represent the flow velocity and density of the fracturing fluid in the wellbore, respectively. and These represent the acceleration due to gravity and the coefficient of friction of the wellbore, respectively. This refers to the direction-dependent drag coefficient in the wellbore. and These are the pressures of the fracturing fluid before and after it flows through the perforation element. This indicates the number of perforations in each cluster; The diameter of the perforation hole. It is a dimensionless flow coefficient used to inversely map the shape changes of orifices caused by wear.

[0034] Step 5: Based on the model constructed above, set multiple sets of parameter conditions and conduct fracturing parameter sensitivity analysis to obtain the influence of key factors such as different perforation parameters, ground stress difference, and hydraulic load on fracture propagation behavior and complexity. Figure 4 The flow distribution diagram for the three fracture clusters under different injection rates is shown. Figure 5 The diagram shows the flow distribution of three fracture clusters under different perforation diameters, illustrating the flow distribution patterns during fracture propagation under different injection rates and different perforation diameters.

[0035] Step Six: Results Visualization. Results include fracture propagation morphology, fluid distribution patterns, and the impact of various fracturing parameters on fracture behavior, providing a basis for optimizing fracturing operations. Figure 6 This is a diagram of crack propagation morphology based on Table 1.

[0036] The present invention also includes: A system comprising a parameter acquisition module, a modeling and simulation module, and a display module, capable of simulating the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference, wherein the parameter acquisition module implements step one of the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference; the modeling and simulation module implements steps two through five of the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference; and the display module implements step six of the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference.

[0037] An apparatus comprising: Memory: Used to store the computer program for the simulation method of multi-cluster fracturing fracture propagation considering the interference of natural fractures; Processor: Used to implement the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference when executing the computer program.

[0038] A computer-readable storage medium storing a computer program that, when executed by a processor, implements the aforementioned method for simulating multi-cluster fracturing fracture propagation considering interference from natural fractures.

[0039] The above description is merely a preferred embodiment of the present invention and is not intended to limit the invention in any way. Any simple modifications, alterations, and equivalent changes made to the above embodiments based on the inventive essence shall fall within the protection scope of the present invention.

Claims

1. A method for simulating the propagation of multi-cluster hydraulic fracturing fractures considering the interference of natural fractures, characterized in that, Includes the following steps: Step 1: Obtain the geological parameters, construction parameters, and completion parameters of the reservoir in the fracturing block; Step 2: Based on some geological parameters, reservoir structural characteristics and geological statistical laws from Step 1, establish a reservoir natural fracture model that reflects the spatial distribution characteristics of natural fractures. Step 3: Based on the reservoir natural fracture model in Step 3, the finite element method and the method of inserting global cohesive elements are used to establish the rock constitutive equation, fluid-structure interaction equation, fluid flow equation within the fracture, fracture propagation criterion equation, cohesive element damage equation, and stiffness degradation criterion equation for the stress field and flow field of the segmented multi-cluster fracturing in a two-dimensional horizontal well. These equations are used to describe the evolution process of the fracture and the rock damage behavior. Step 4: Based on Step 3, introduce a pipe flow unit to simulate the non-uniform distribution behavior of fluid flow between clusters under the condition of simultaneous expansion of multiple fractures, and establish a dynamic flow distribution model that considers fracture competition effect and inter-cluster interference; a dynamic fluid distribution model between fractures; and the energy conservation equation, wellbore friction pressure drop equation and perforation hole friction pressure drop equation in the fluid system. Step 5: Based on the model constructed in Step 4, set up multiple sets of parameter conditions to explore the influence of different perforation parameters, ground stress difference, and hydraulic load on the crack propagation behavior and complexity. Step 6: Visualize the results, including fracture propagation morphology, fluid distribution patterns, and the influence of various fracturing parameters on fracture behavior.

2. The method for simulating the propagation of multi-cluster fracturing fractures considering the interference of natural fractures according to claim 1, characterized in that, The geological parameters mentioned in step one include natural fracture angle, natural fracture length, Young's modulus of rock, Poisson's ratio of rock, shear modulus of rock, and filtration coefficient; Construction parameters include fracturing fluid density, fracturing fluid viscosity, and injection rate. Completion parameters include the number of perforations and the perforation diameter.

3. The method for simulating the propagation of multi-cluster fracturing fractures considering the interference of natural fractures according to claim 1, characterized in that, The establishment of the reservoir natural fracture model in step two is specifically as follows: The starting point of the crack ( , The fractures are randomly distributed within the reservoir to ensure minimal distance between them. To avoid overlap, the starting position of the crack is randomly generated using the following formula: in, This refers to the reservoir's edge length; the `Random()` function generates random numbers within a specified range. The fracture length is also mentioned. and angle By specifying the minimum length and maximum length and minimum angle and maximum angle To generate randomly: The other end of the crack ( , Based on the starting point ( , ), length of the crack and angle Perform the calculation: （3） 。 4. The method for simulating the propagation of multi-cluster fracturing fractures considering the interference of natural fractures according to claim 1, characterized in that, Step three specifically involves: 3.1 Rock constitutive equation: The linear elastic constitutive equation of rock is shown below: in, It is Poisson's ratio (dimensionless); It is the shear modulus. ; , For Young's modulus, , where is the Biot coefficient (dimensionless). Pore ​​pressure, Constitutive equations based on stress and solid strain express; 3.

2. Establish fluid-structure interaction equations, including the coupling control equations for rock solid skeleton deformation and fluid flow, the mass conservation equation for fluid seepage, and the fluid flow velocity equation within the rock, expressed as follows: in, For effective stress, It is the identity matrix; For virtual strain rate, It is the surface force vector per unit area. , It is the volume force per unit volume. ; It is a virtual velocity field. ; per unit volume ; Represents the area of ​​the unit subjected to surface forces. ; It represents the volume change of fluid in a porous medium (dimensionless). For fluid density, ; The porosity is related to the porosity of porous media; For spatial vectors, ; The seepage velocity of the fluid. ; For time; It is the acceleration due to gravity. , The penetration matrix, ; 3.3 Establish the fluid flow equations within the fracture, including the tangential fluid flow equation, the fluid mass conservation equation, and the fracturing fluid loss equation, as specifically expressed below: In the above formula, This refers to tangential flow within the crack. , The pressure gradient is along the length of the element. , The viscosity of the fracturing fluid. , The thickness of the crack propagation. , , The volumetric flow rate per unit time is the volumetric flow rate of the upper and lower surfaces of the crack. ; , The filtration coefficients are the upper and lower surface filtration coefficients. , The pore pressure at the upper and lower surfaces of the crack. ; The fluid pressure of the unit. ; 3.

4. Establish the crack propagation criterion equation, the cohesive element damage equation, and the stiffness degradation criterion equation for the damage evolution of the cohesive element, as detailed below: In the above formula, , , These represent the critical normal stress and the critical tangential stresses in two directions at the point of cohesive unit failure. , , It is the normal stress component before the damage. It is the tangential stress component before the damage; This indicates the overall degree of damage to the cohesive unit, which increases linearly from 0 to 1 after the element damage begins. This represents the maximum displacement amplitude reached by the element during the loading process; This represents the displacement amplitude when the element is completely destroyed. This represents the displacement amplitude at the initial damage stage of the element.

5. The method for simulating the propagation of multi-cluster fracturing fractures considering the interference of natural fractures according to claim 1, characterized in that, In step four, the dynamic fluid distribution equations among multiple fracture clusters are established, including the fluid energy conservation equation, the wellbore friction pressure drop equation, and the perforation hole friction pressure drop equation, as detailed below: In the above formula, It is the injection rate. Is entering the cluster The flow rate, It is the number of clusters. and These represent the pressure drop as fluid passes through the wellbore element and the vertical distance between the nodes at both ends of the wellbore element, respectively. and These represent the flow velocity and density of the fracturing fluid in the wellbore, respectively. and These represent the acceleration due to gravity and the coefficient of friction of the wellbore, respectively. This refers to the direction-dependent drag coefficient in the wellbore. and These are the pressures of the fracturing fluid before and after it flows through the perforation element. This indicates the number of perforations in each cluster; The diameter of the perforation hole; It is a dimensionless flow coefficient used to inversely map the shape changes of orifices caused by wear.

6. A system, characterized in that, The system includes a parameter acquisition module, a modeling and simulation module, and a display module, and is capable of simulating the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference as described in any of claims 1-5. The parameter acquisition module implements step one of the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference as described in any of claims 1-5; the modeling and simulation module implements steps two through five of the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference as described in any of claims 1-5; and the display module implements step six of the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference as described in any of claims 1-5.

7. A device, characterized in that, include: Memory: for storing the computer program of the multi-cluster fracturing fracture propagation simulation method considering the interference of natural fractures as described in any one of claims 1-5; Processor: Used to implement the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference when executing any of the computer programs described in claims 1-5.

8. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the multi-cluster fracturing fracture propagation simulation method considering natural fracture interference as described in any one of claims 1-5.

Images