CN113431542A - Method for calculating interference strength of horizontal well fracturing fracture - Google Patents

Method for calculating interference strength of horizontal well fracturing fracture Download PDF

Info

Publication number
CN113431542A
CN113431542A CN202010733942.6A CN202010733942A CN113431542A CN 113431542 A CN113431542 A CN 113431542A CN 202010733942 A CN202010733942 A CN 202010733942A CN 113431542 A CN113431542 A CN 113431542A
Authority
CN
China
Prior art keywords
fracturing
fracture
stress
horizontal well
well
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010733942.6A
Other languages
Chinese (zh)
Other versions
CN113431542B (en
Inventor
王增林
张潦源
陈勇
王东英
丁然
杨峰
苏权生
王丽萍
宋李煜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Research Institute of Petroleum Engineering Shengli Co
Original Assignee
China Petroleum and Chemical Corp
Sinopec Research Institute of Petroleum Engineering Shengli Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Research Institute of Petroleum Engineering Shengli Co filed Critical China Petroleum and Chemical Corp
Priority to CN202010733942.6A priority Critical patent/CN113431542B/en
Publication of CN113431542A publication Critical patent/CN113431542A/en
Application granted granted Critical
Publication of CN113431542B publication Critical patent/CN113431542B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/25Methods for stimulating production
    • E21B43/26Methods for stimulating production by forming crevices or fractures
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/06Measuring temperature or pressure

Abstract

The invention provides a method for calculating the interference strength of horizontal well fracturing fractures, which comprises the following steps: step 1, establishing a geomechanical model containing a horizontal well section according to geological data of a development block and well logging data of a development well; 2, carrying out numerical simulation on the fracture initiation and expansion processes of the fracturing fractures of the single-section two adjacent clusters of perforations; step 3, intercepting a vertical section along a well trace of the horizontal well; step 4, adding the near-field stress radiation radii of two adjacent clusters of cracks, and comparing the near-field stress radiation radii with the interval of the fracturing design clusters; step 5, comparing the measured crack steering angle with a design threshold value; and 6, comparing the measured crack development radius with a design threshold value. The method for calculating the interference strength of the horizontal well fracturing fracture can be used for optimizing process parameters such as cluster spacing, construction displacement and construction scale of the horizontal well fracturing process, and provides technical support for improving effective development of oil fields.

Description

Method for calculating interference strength of horizontal well fracturing fracture
Technical Field
The invention relates to the technical field of petroleum and natural gas yield increase, in particular to a method for calculating the interference strength of a horizontal well fracturing fracture.
Background
The 'close cutting, three-dimensional and ultra-long horizontal well' is a new breakthrough of modern volume transformation technology and application, and the core of the method is to further shorten the seepage distance of fluid in a matrix to a crack, greatly reduce driving pressure difference and increase the contact area of the matrix and the crack. The horizontal well dense cutting staged fracturing process is beneficial to increasing the inter-cluster stress interference, improving the complexity of cracks and improving the inter-cluster resource utilization efficiency. The method aims to obtain the maximum reservoir stratum modification volume, overcomes the difference value of two-directional horizontal main stress by using the induced stress formed nearby the fracture, and improves the fracture complexity by deflecting the fracture expansion direction or communicating the natural fracture. Optimizing the cluster spacing and the construction process parameters, enhancing the stress interference degree between cracks, and improving the complexity degree of the cracks and the permeability-increasing transformation degree of a reservoir stratum are the keys of the close-cut fracturing technology.
The method is characterized in that clustering perforation is the key of application of a volume reconstruction technology, when each fracturing section adopts a plurality of clusters of perforation (3 clusters or more), the key of ensuring the opening of each cluster under constant discharge capacity is to limit the number of the perforations in the fracturing section, if the total number of the perforations can ensure that each cluster has enough throttling resistance, all the perforation clusters can be opened, and the clusters which cannot be opened do not need to be opened by adopting a temporary blocking technology in the sections. Due to the influence of factors such as reservoir heterogeneity and perforation hole phase, how to realize the balanced reconstruction of each cluster needs to be analyzed from the aspect of multi-fracture expansion. Whether staged simultaneous fracturing, staged sequential fracturing or alternate fracturing and zipper type synchronous fracturing, the mutual interference among fractures and the control on the fracture morphology are the key points for ideal on-site fracturing effect. Through a large amount of practical cognition on the mutual interference rule existing among a plurality of fracturing clusters in a single fracturing section in the staged multi-cluster hydraulic fracturing process of a horizontal well of a shale oil-gas reservoir, the cluster spacing is considered to be the most main factor influencing the inter-fracture interference, the inter-cluster stress interference is gradually reduced along with the increase of the fracturing cluster spacing, and the expansion of each fracture tends to be uniform; the increase of the formation elastic modulus can promote the longitudinal extension of the cracks of each perforation cluster, but can prevent the transverse opening of the cracks; in addition, the appropriate increase of the discharge capacity of the fracturing fluid in the fracturing process is beneficial to obtaining the long and wide fractures. Although this behavior of interference between fractures is recognized, there is currently no effective, convenient way to quantify the intensity of the interference.
From the geomechanics perspective, the single crack plane always extends along the direction of the maximum horizontal main stress when expanding, and the influence range of the generated crack induced stress field is mainly controlled by the minimum value of the length and the height of the crack; when the single crack is not expanded in a plane, the extending direction of the single crack gradually deviates to the direction of the maximum horizontal main stress, and the crack form is bent as the difference between the initial crack azimuth angle and the main stress is smaller.
In the application No.: 202010035516.5, relates to a method for evaluating a fracture by using an interference well testing theory. The method comprises the steps of surveying the communication condition of a pressure drive well and surrounding adjacent wells, and collecting the well cementation quality of the pressure drive well, a pressure drive target horizon, dynamic production information and a sedimentary phase-to-phase diagram; determining a pressure driving pressure real-time monitoring well; determining the middle depth of an oil reservoir of a pressure real-time monitoring well, the distance between the monitoring well and a pressure drive well, the dynamic production condition, the production pipe column, the well mouth and the well field condition, and determining the testing mode to be pressure gauge real-time monitoring; formulating a real-time monitoring scheme, and determining the testing construction steps and requirements; the pressure of an adjacent well in the pressure flooding process is monitored in real time, and original data are checked and accepted; monitoring a well pressure curve, and performing pressure drive disturbance time lag analysis; and performing pressure drive crack evaluation according to the pressure drive interference time lag interpretation result. The technology can not calculate the fracture interference strength, so that the technological parameters of the horizontal well fracturing technology, such as cluster spacing, construction displacement, construction scale and the like, can not be optimized.
In the application No.: 201711483019.6, relates to a compact reservoir volume fracture network expansion simulation and characterization method, including the following steps: the method comprises the steps that firstly, a compact reservoir multi-fracture stress interference combined ground stress field calculation model considering multi-fracture stress interference is established by utilizing a displacement discontinuity method, mechanical mechanism analysis and a starting and expanding rule, and in the process of solving stress and geometric parameters of each newly added fracture infinitesimal, an additional stress field including a superimposed stress field generated by normal stress and shear stress action must be recalculated at each time step, and the additional stress field must be superimposed on a stress field of the last time step to finally determine new combined stress field distribution in a global coordinate system XOY; step two, simulating the fracture network expansion of the compact reservoir volume fracturing horizontal well, aiming at the problem of stress interference in the volume fracturing multi-fracture expansion process, establishing a flowing pressure drop distribution model of fracturing fluid in the main fracture and the secondary fracture, combining the volume fracturing fracture expansion combined geostress field calculation model established in the step one, forming a compact reservoir volume fracturing horizontal well fracture network expansion theoretical model, and simulating the compact reservoir horizontal well volume fracturing fracture network forming process; and step three, describing and characterizing the space structure of multiple pores of the complex fracture network, comprehensively analyzing the influence of different factors on the structural form of the volume fracture network, and characterizing the structural form and the attribute characteristics of the fracture network by defining a plurality of characteristic parameters. Although the stress interference problem in the multi-fracture expansion process is considered, the fracture interference strength is not defined, and the optimization method of the key parameters of the horizontal well fracturing process is not set forth, so that the process parameters of the horizontal well fracturing process, such as cluster spacing, construction displacement, construction scale and the like, cannot be optimized.
Therefore, a novel method for calculating the interference strength of the horizontal well fracture is invented, and the technical problems are solved.
Disclosure of Invention
The invention aims to provide a method for calculating the interference strength of the horizontal well fracturing fracture, which is simple to implement, low in cost and capable of operating repeatedly from the perspective of joint control of the length and the bending degree of the fracture.
The object of the invention can be achieved by the following technical measures: the method for calculating the interference strength of the horizontal well fracturing fracture comprises the following steps: step 1, establishing a geomechanical model containing a horizontal well section according to geological data of a development block and well logging data of a development well; 2, carrying out numerical simulation on the fracture initiation and expansion processes of the fracturing fractures of the single-section two adjacent clusters of perforations; step 3, intercepting a vertical section along a well trace of the horizontal well; step 4, adding the near-field stress radiation radii of two adjacent clusters of cracks, and comparing the near-field stress radiation radii with the interval of the fracturing design clusters; step 5, comparing the measured crack steering angle with a design threshold value; and 6, comparing the measured crack development radius with a design threshold value.
The object of the invention can also be achieved by the following technical measures:
in step 1, according to geological data of a development block and well logging data of a development well, a geomechanical model containing a horizontal well section is established, wherein the geomechanical model comprises setting geometric parameters, geomechanical parameters and internal and external boundary conditions of the model.
In step 1, when determining the geometric dimensions, the geometric dimensions of the selected model in the direction of the maximum horizontal principal stress and the direction of the minimum horizontal principal stress are 1000m respectively, and the dimensions in the vertical ground stress direction are increased by 500m respectively from top to bottom by taking the horizontal well trajectory as the center.
In the step 1, when determining geomechanical parameters, according to oil deposit geological data, continuous logging data and actual core rock mechanical test data of coring, the elastic modulus, Poisson's ratio, cohesion, internal friction angle, tensile strength, porosity and permeability of a reservoir to be detected and an interlayer are obtained.
In step 1, when determining the boundary condition of the model, the original ground stress data of the actual block is obtained to include the maximum horizontal principal stress sigmaHMinimum horizontal principal stress σhMaximum vertical principal stress σVWill σH、σhAnd σVRespectively applied to six surfaces of the model; and (3) selecting the design position of each cluster of perforation at any section along the well trace of the horizontal well in the model, applying a constant flow boundary condition, and taking the value of the flow as a fracturing construction design parameter.
In the step 2, based on the geomechanical model, according to fracturing process design parameters, carrying out numerical simulation on the initiation and expansion processes of fracturing fractures of two adjacent single-section clusters of perforations, wherein the calculated time length is a fracturing construction design parameter, and the total injection liquid amount is a fracturing construction design parameter; in the simulation process, the simulation can obtain the initiation, the extension process and the stress evolution in the extension process of each cluster of cracks, the tips of the cracks have high stress concentration before the fracture initiation of the fracturing cracks, and once the cracks initiate, the stress is released immediately and transferred to the near-field region of the cracks.
In step 3, after the expansion process and the stress field evolution diagram of the fracturing fracture are obtained through simulation, a vertical section along a horizontal well trajectory is intercepted, the included angles between the two clusters of fracturing fracture trajectories and the normal line of the horizontal well trajectory are measured respectively, and the steering angles a1 and a2 of the non-planar fracture are obtained; and then measuring the respective near-field stress radiation radiuses r1 and r2 of the two adjacent clusters of fractured fractures on the section respectively, wherein the radiation radiuses are judged by comparing the stress value of the fracture near-field region with the initial ground stress value, and the stress values in the fracture near-field region are more than or equal to the initial ground stress value.
In step 4, adding the near-field stress radiation radii of two adjacent clusters of cracks, if the sum of the two near-field stress radiation radii is smaller than the fracturing design cluster spacing, further optimizing the process parameters of the fracturing process, such as the cluster spacing, the construction displacement and the construction scale, and returning the process to step 1; and if the sum of the two is more than or equal to the fracturing design cluster spacing, the flow enters the step 5.
In step 5, comparing the measured fracture steering angle with a design threshold, if the measured fracture steering angle is larger than the design threshold, optimizing process parameters of the fracturing process, such as cluster spacing, construction displacement and construction scale, and returning the process to step 1; if the measured fracture steering angle is less than or equal to the design threshold, the process proceeds to step 6.
In step 6, comparing the measured crack development radius with a design threshold, if the measured crack development radius is smaller than the design threshold, optimizing the process parameters of the fracturing process, such as cluster spacing, construction displacement and construction scale, and returning the process to step 1; and if the measured crack development radius is larger than or equal to the design threshold, ending the process.
The method for calculating the interference strength of the horizontal well fracturing fracture can optimize process parameters such as cluster spacing, construction displacement, construction scale and the like of the horizontal well fracturing process. The method realizes quantitative analysis of the interference strength of the staged fracturing fracture of the horizontal well by jointly controlling three parameters, namely the fracturing fracture steering angle, the fracture development radius and the radiation radius of the stress shadow area, can be used for optimizing process parameters such as cluster spacing, construction displacement, construction scale and the like of the horizontal well fracturing process, and provides technical support for improving effective development of oil fields.
Drawings
FIG. 1 is a schematic representation of a three-dimensional geomechanical model of a single fracture section containing a horizontal well in one embodiment of the present invention;
FIG. 2 is a graph showing respective near-field stress radiation radii of two adjacent fracture clusters in an embodiment of the present invention;
FIG. 3 illustrates the included angle between the two non-planar fracture trajectories and the normal of the horizontal well trajectory and the fracture development radius in an embodiment of the present invention;
FIG. 4 is a result of a quantitative calculation of fracture disturbance intensity based on near-field stress radiation radius, fracture growth radius, and steering angle in an embodiment of the present invention;
fig. 5 is a flowchart of a specific embodiment of the method for calculating the interference strength of the horizontal well fracture.
Detailed Description
It is to be understood that the following detailed description is exemplary and is intended to provide further explanation of the invention as claimed. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the invention. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of the stated features, steps, operations, and/or combinations thereof, unless the context clearly indicates otherwise.
In order to make the technical solutions of the present invention more clearly understood by those skilled in the art, the technical solutions of the present invention will be described in detail below with reference to specific embodiments.
The instruments, reagents, materials and the like used in the following examples are conventional instruments, reagents, materials and the like in the prior art and are commercially available in a normal manner unless otherwise specified. Unless otherwise specified, the experimental methods, detection methods, and the like described in the following examples are conventional experimental methods, detection methods, and the like in the prior art.
If the calculation and analysis of the single factor of the stress shadow zone of the staged fracturing fracture of the horizontal well are only carried out singly, the radiation area of the stress shadow zone hardly reflects the actual mutual interference intensity between two adjacent clusters of fractures accurately; if the calculation of fracture development height (radius) is performed only singly, it is not sufficient to make an accurate judgment on the effectiveness of the fracture initiation extension path. Therefore, the influence and the connection of the stress field of the three-dimensional geomechanical model and the internal cause of the time-space coordinate formed by the crack can be fully considered, and three measurement indexes of the crack deflection angle, the crack development radius and the stress shadow radiation radius are organically combined from different angles and are mutually complemented: the stress shadow region can help to understand and deduce the origin and mode of the observed crack initiation and propagation; the initiation propagation path of the three-dimensional fracture can help explain the law of stress concentration, release and migration of the region.
Based on the thought, quantitative analysis of the interference intensity of the horizontal well staged fracturing fracture is realized through the joint control of three parameters, namely the fracturing fracture steering angle, the fracture development radius and the stress shadow area radiation radius, and as shown in fig. 5, fig. 5 is a flow chart of the method for calculating the interference intensity of the horizontal well fracturing fracture.
Step 101, establishing a geomechanical model containing a horizontal well section according to geological data of a development block and well logging data of a development well, wherein the geomechanical model comprises setting geometric parameters, geomechanical parameters and internal and external boundary conditions of the model;
the geomechanical model of the horizontal well section comprises three aspects of geometric dimension, geomechanical parameters and boundary conditions inside and outside the model. Geometric dimension determination: the geometric dimensions of the selected model in the direction of the maximum horizontal principal stress and the direction of the minimum horizontal principal stress are respectively 1000m, and the dimensions in the vertical ground stress direction (Z direction) are respectively increased by 500m from top to bottom by taking the horizontal well trace as the center. Determining geomechanical parameters: according to the oil deposit geological data, the continuous logging data and the actual coring core rock mechanical test data, the method comprises the following steps: the elastic modulus, Poisson's ratio, cohesion, internal friction angle, tensile strength, porosity and permeability of the reservoir and the interlayer are required. Model boundary conditions: obtaining the data including the maximum horizontal principal stress sigma according to the original ground stress data of the actual blockHMinimum horizontal principal stress σhMaximum vertical principal stress σVWill σH、σhAnd σVApplied to six faces of the model (Z, X, Y three directions), respectively; and (3) selecting the design position of each cluster of perforation at any section along the well trace of the horizontal well in the model, applying a constant flow boundary condition, and taking the value of the flow as a fracturing construction design parameter.
102, performing numerical simulation of the fracture initiation and expansion process of the fracturing fractures of the two adjacent clusters of single-section perforations by using numerical simulation software according to fracturing process design parameters based on the geomechanical model;
the simulation of the fracturing process of the horizontal well geomechanical model can be carried out on any numerical calculation software capable of simulating hydraulic fracturing fracture expansion, the calculation time is fracturing construction design parameters, and the total injection liquid amount is fracturing construction design parameters. In the simulation process, the simulation can obtain the initiation, the extension process and the stress evolution in the extension process of each cluster of cracks, the tips of the cracks have high stress concentration before the fracture initiation of the fracturing cracks, and once the cracks initiate, the stress is released immediately and transferred to the near-field region of the cracks.
Step 103, intercepting a vertical section along a horizontal well trajectory, and firstly respectively measuring the included angle between each of two fracturing fracture trajectories and the normal of the horizontal well trajectory to obtain the steering angle of the non-planar fracture; then respectively measuring the respective average development radius of the two clusters of fracturing fractures; further measuring the respective near-field stress radiation radii of the two adjacent clusters of fracturing fractures on the section;
after the expansion process and the stress field evolution diagram of the fracturing fracture are obtained through simulation, a vertical section along a horizontal well trajectory is intercepted, the included angles between the two fracturing fracture trajectories and the normal line of the horizontal well trajectory are measured respectively, and the steering angles a1 and a2 of the non-planar fracture are obtained; and then measuring the respective near-field stress radiation radiuses r1 and r2 of the two adjacent clusters of fractured fractures on the section respectively, wherein the radiation radiuses are judged by comparing the stress value of the fracture near-field region with the initial ground stress value, and the stress values in the fracture near-field region are more than or equal to the initial ground stress value.
104, adding the near-field stress radiation radii of two adjacent clusters of cracks, and if the sum of the two near-field stress radiation radii is smaller than the fracturing design cluster spacing, further optimizing the process parameters of the fracturing process, such as the cluster spacing, the construction displacement, the construction scale and the like, and repeating the steps;
and adding the measured near-field stress radiation radiuses r1 and r2 of two adjacent clusters of cracks, and if the sum of the two radiuses is smaller than the fracturing design cluster spacing, further optimizing the cluster spacing, construction displacement, construction scale and other process parameters of the fracturing process.
Step 105, comparing the measured crack steering angle with a design threshold, if the measured crack steering angle is larger than the design threshold, optimizing process parameters of the fracturing process, such as cluster spacing, construction displacement, construction scale and the like, and repeating the steps; and comparing the measured crack steering angle a with a design threshold value, wherein the design threshold value is generally 15-30 degrees, and if the steering angle is larger than the set threshold value, optimizing the process parameters of the fracturing process, such as cluster spacing, construction displacement, construction scale and the like.
Step 106, comparing the measured crack development radius with a design threshold, if the measured crack development radius is smaller than the design threshold, optimizing the process parameters of the fracturing process, such as cluster spacing, construction displacement, construction scale and the like, and repeating the steps;
and comparing the measured crack development radius h with a design threshold value, wherein the design threshold value is generally between 100 and 200m, if the design threshold value is smaller than the design threshold value, optimizing the process parameters of the fracturing process, such as cluster spacing, construction displacement, construction scale and the like, and repeating the steps.
In one embodiment of the present invention, the method comprises the following steps:
step one, establishing a geomechanical model containing a horizontal well section, as shown in figure 1, wherein the geometric size of the model is as follows: the geometric dimensions of the selected model in the direction of the maximum horizontal principal stress and the direction of the minimum horizontal principal stress are respectively 1000m, and the dimensions in the vertical ground stress direction (Z direction) are respectively increased by 500m from top to bottom by taking the horizontal well trace as the center. Determining geomechanical parameters: according to the geological data of the oil deposit, the continuous logging data and the actual coring rock core rock mechanical test data, the concrete geomechanical parameters are as follows: the elastic modulus of the reservoir is 25GPa, the Poisson ratio is 0.22, the cohesive force is 19MPa, the internal friction angle is 30 degrees, the tensile strength is 3MPa, the porosity is 12 percent, and the permeability is 1 mD. The elastic modulus of the interlayer is 35GPa, the Poisson ratio is 0.20, the cohesive force is 25MPa, the internal friction angle is 35 degrees, the tensile strength is 4MPa, the porosity is 7 percent, and the permeability is 0.6 mD. Model outer boundary conditions: maximum horizontal principal stress σHMinimum level principal stress sigma 55MPah48MPa, vertical principal stress sigmaV65MPa,. sigmaH、σhAnd σVApplied to six faces of the model (Z, X, Y three directions), respectively; two clusters of perforations at one section are selected as objects along the well trace of the horizontal well in the model, a constant flow boundary condition is applied, the flow value is 8m3/min, and the cluster spacing is 30 m.
And step two, simulating the fracturing process of the horizontal well geomechanical model by using numerical simulation software RFPA, wherein the calculation time is 70min, and the total injection liquid amount is a fracturing construction design parameter. In the simulation process, the simulation can obtain the initiation, the extension process and the stress evolution in the extension process of each cluster of cracks, the tips of the cracks have high stress concentration before the fracture initiation of the fracturing cracks, and once the cracks initiate, the stress is released immediately and transferred to the near-field region of the cracks.
Thirdly, after the expansion process and the stress field evolution diagram of the fracturing fracture are obtained through simulation, a vertical section along the horizontal well trajectory is cut, the included angle between each of the two clusters of fracturing fracture trajectories and the normal line of the horizontal well trajectory is measured, the steering angle a1 of the non-planar fracture is 15 degrees, the steering angle a2 of the non-planar fracture is 32 degrees, and as shown in fig. 2, the brightness of the gray scale represents the stress; measuring the respective development radius h 1-120 m and h 2-100 m of the two clusters of cracks; further, the respective near-field stress radiation radii r1 and r2 of the two adjacent clusters of fracture fractures on the profile are measured respectively to be 10m and 12m, the radiation radii are judged by comparing the stress values of the fracture near-field region with the initial ground stress values, and the stress values in the fracture near-field region are all larger than or equal to the initial ground stress values, as shown in fig. 3.
Step four, adding the measured near-field stress radiation radii r1 and r2 of two adjacent clusters of cracks, if the sum of the two radii is 22m and is smaller than the fracturing design cluster spacing, further optimizing the cluster spacing of the fracturing process, reducing the cluster spacing to 10m, repeating the step one to the step three, enabling the near-field stress radiation radii of the cracks and the fracture steering angle to accord with the design threshold value, and finally determining the reasonable cluster spacing, as shown in fig. 4.
And step five, comparing the measured fracture steering angle a with a design threshold, wherein the value of the design threshold is 15 degrees in the embodiment, and the steering angle is larger than the set threshold, so that the cluster spacing of the fracturing process is optimized, the cluster spacing is increased to 15m, and repeating the steps from the first step to the fourth step.
And step six, comparing the measured crack development radius h with a design threshold, wherein the value of the design threshold is 100m in the embodiment, and the two crack development radii are both larger than the set threshold, so that the cluster spacing of the optimized fracturing process is reasonable to 15 m.
The above embodiments are preferred embodiments of the present invention, but the present invention is not limited to the above embodiments, and any other changes, modifications, substitutions, combinations, and simplifications which do not depart from the spirit and principle of the present invention should be construed as equivalents thereof, and all such changes, modifications, substitutions, combinations, and simplifications are intended to be included in the scope of the present invention.

Claims (10)

1. The method for calculating the interference strength of the horizontal well fracturing fracture is characterized by comprising the following steps:
step 1, establishing a geomechanical model containing a horizontal well section according to geological data of a development block and well logging data of a development well;
2, carrying out numerical simulation on the fracture initiation and expansion processes of the fracturing fractures of the single-section two adjacent clusters of perforations;
step 3, intercepting a vertical section along a well trace of the horizontal well;
step 4, adding the near-field stress radiation radii of two adjacent clusters of cracks, and comparing the near-field stress radiation radii with the interval of the fracturing design clusters;
step 5, comparing the measured crack steering angle with a design threshold value;
and 6, comparing the measured crack development radius with a design threshold value.
2. The method for calculating the interference strength of the horizontal well fracturing fracture according to claim 1, wherein in the step 1, a geomechanical model containing a horizontal well section is established according to geological data of a development block and well logging data of a development well, and the geomechanical model comprises the setting of geometric parameters, geomechanical parameters and internal and external boundary conditions of the model.
3. The method for calculating the interference strength of the horizontal well fracture according to claim 2, wherein in the step 1, when determining the geometric dimensions, the geometric dimensions of the selected model in the direction of the maximum horizontal principal stress and the direction of the minimum horizontal principal stress are respectively 1000m, and the dimensions in the direction of the vertical ground stress are respectively increased by 500m from top to bottom by taking the well trace of the horizontal well as the center.
4. The method for calculating the interference strength of the horizontal well fracturing fracture according to claim 2, wherein in the step 1, when the geomechanical parameters are determined, the elastic modulus, the Poisson's ratio, the cohesion, the internal friction angle, the tensile strength, the porosity and the permeability of the reservoir to be tested and the interlayer are obtained according to the reservoir geological data, the continuous logging data and the actually cored core rock mechanical test data.
5. The method for calculating the disturbance intensity of the horizontal well fracturing fracture according to claim 2, wherein in step 1, when determining the model boundary condition, the original ground stress data of the actual block is obtained to include the maximum horizontal principal stress σHMinimum horizontal principal stress σhMaximum vertical principal stress σVWill σH、σhAnd σVRespectively applied to six surfaces of the model; and (3) selecting the design position of each cluster of perforation at any section along the well trace of the horizontal well in the model, applying a constant flow boundary condition, and taking the value of the flow as a fracturing construction design parameter.
6. The method for calculating the interference strength of the horizontal well fracturing fractures according to claim 1, wherein in the step 2, numerical simulation of the fracture initiation and expansion process of the fracturing fractures of two adjacent single-stage perforating clusters is carried out according to fracturing process design parameters on the basis of the geomechanical model, the calculation time is fracturing construction design parameters, and the total injection liquid amount is fracturing construction design parameters; in the simulation process, the simulation can obtain the initiation, the extension process and the stress evolution in the extension process of each cluster of cracks, the tips of the cracks have high stress concentration before the fracture initiation of the fracturing cracks, and once the cracks initiate, the stress is released immediately and transferred to the near-field region of the cracks.
7. The method for calculating the interference strength of the horizontal well fracturing fractures according to claim 1, wherein in step 3, after the expansion process and the stress field evolution diagram of the fracturing fractures are obtained through simulation, vertical sections along the horizontal well trajectory are intercepted, the included angles between the two clusters of fracturing fracture trajectories and the normal line of the horizontal well trajectory are measured respectively, and the steering angles a1 and a2 of the non-planar fractures are obtained; and then measuring the respective near-field stress radiation radiuses r1 and r2 of the two adjacent clusters of fractured fractures on the section respectively, wherein the radiation radiuses are judged by comparing the stress value of the fracture near-field region with the initial ground stress value, and the stress values in the fracture near-field region are all larger than or equal to the initial ground stress value.
8. The method for calculating the interference strength of the horizontal well fracturing fractures according to claim 1, wherein in step 4, the near-field stress radiation radii of two adjacent clusters of fractures are added, if the sum of the near-field stress radiation radii of the two adjacent clusters of fractures is smaller than the fracturing design cluster spacing, the process parameters of the fracturing process, such as the cluster spacing, the construction displacement and the construction scale, are further optimized, and the process returns to step 1; and if the sum of the two is more than or equal to the fracturing design cluster spacing, the flow enters the step 5.
9. The method for calculating the interference strength of the horizontal well fracturing fracture according to claim 1, wherein in the step 5, the measured fracture steering angle is compared with a design threshold, if the measured fracture steering angle is larger than the design threshold, process parameters such as cluster spacing, construction displacement and construction scale of a fracturing process are optimized, and the process returns to the step 1; if the measured fracture steering angle is less than or equal to the design threshold, the process proceeds to step 6.
10. The method for calculating the interference strength of the horizontal well fracturing fracture according to claim 1, wherein in the step 6, the measured fracture development radius is compared with a design threshold, if the fracture development radius is smaller than the design threshold, process parameters such as cluster spacing, construction displacement and construction scale of a fracturing process are optimized, and the process returns to the step 1; and if the measured crack development radius is larger than or equal to the design threshold, ending the process.
CN202010733942.6A 2020-07-27 2020-07-27 Method for calculating interference strength of horizontal well fracturing fracture Active CN113431542B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010733942.6A CN113431542B (en) 2020-07-27 2020-07-27 Method for calculating interference strength of horizontal well fracturing fracture

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010733942.6A CN113431542B (en) 2020-07-27 2020-07-27 Method for calculating interference strength of horizontal well fracturing fracture

Publications (2)

Publication Number Publication Date
CN113431542A true CN113431542A (en) 2021-09-24
CN113431542B CN113431542B (en) 2023-04-18

Family

ID=77753300

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010733942.6A Active CN113431542B (en) 2020-07-27 2020-07-27 Method for calculating interference strength of horizontal well fracturing fracture

Country Status (1)

Country Link
CN (1) CN113431542B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114707440A (en) * 2022-03-31 2022-07-05 中国石油大学(北京) Multi-crack size calculation method and device and computer equipment

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1994021889A2 (en) * 1993-03-17 1994-09-29 John North Improvements in or relating to drilling and to the extraction of fluids
CN102033242A (en) * 2010-10-22 2011-04-27 中国石油化工股份有限公司 Deep inclined fractured reservoir earthquake amplitude prediction method
CN105201479A (en) * 2015-10-09 2015-12-30 西南石油大学 Optimum design method for staged fracturing perforation cluster parameters of horizontal well of shale reservoir
CN105913155A (en) * 2016-04-21 2016-08-31 中国石油天然气股份有限公司 Tight oil productivity prediction method considering stress interference and fracturing fluid filtration loss and tight oil productivity prediction system thereof
CN105986796A (en) * 2015-02-10 2016-10-05 中国石油大学(北京) Experimental method of fracture interaction of adjacent parallel hydraulic fractures in hydraulic fracturing process
US20170051598A1 (en) * 2015-08-20 2017-02-23 FracGeo, LLC System For Hydraulic Fracturing Design And Optimization In Naturally Fractured Reservoirs
CN106593390A (en) * 2016-12-23 2017-04-26 西南石油大学 Method for simulating horizontal well staged fracturing fracture expanding track
CN106650100A (en) * 2016-12-23 2017-05-10 西南石油大学 Volume alternating fracturing method of horizontal well in experimental shale reservoir
CN107506840A (en) * 2016-06-14 2017-12-22 中国石油化工股份有限公司 A kind of optimization method of shale gas horizontal well perforation cluster spacing
CN107992690A (en) * 2017-12-06 2018-05-04 西南石油大学 A kind of evaluation method of the lower multiple cracking expanded configuration balance degree of induced stress interference
CN109933860A (en) * 2019-02-25 2019-06-25 北京科技大学 Shale gas pressure break horizontal well cluster spacing optimization method
CN111255442A (en) * 2020-01-14 2020-06-09 大庆油田有限责任公司 Method for evaluating fracturing fracture by using interference well testing theory
US10677961B1 (en) * 2019-07-08 2020-06-09 Southwest Petroleum University Method for optimizing perforation parameters to maintain uniform fracture growth in multi-stage hydraulic fracturing of horizontal well

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1994021889A2 (en) * 1993-03-17 1994-09-29 John North Improvements in or relating to drilling and to the extraction of fluids
CN102033242A (en) * 2010-10-22 2011-04-27 中国石油化工股份有限公司 Deep inclined fractured reservoir earthquake amplitude prediction method
CN105986796A (en) * 2015-02-10 2016-10-05 中国石油大学(北京) Experimental method of fracture interaction of adjacent parallel hydraulic fractures in hydraulic fracturing process
US20170051598A1 (en) * 2015-08-20 2017-02-23 FracGeo, LLC System For Hydraulic Fracturing Design And Optimization In Naturally Fractured Reservoirs
CN105201479A (en) * 2015-10-09 2015-12-30 西南石油大学 Optimum design method for staged fracturing perforation cluster parameters of horizontal well of shale reservoir
CN105913155A (en) * 2016-04-21 2016-08-31 中国石油天然气股份有限公司 Tight oil productivity prediction method considering stress interference and fracturing fluid filtration loss and tight oil productivity prediction system thereof
CN107506840A (en) * 2016-06-14 2017-12-22 中国石油化工股份有限公司 A kind of optimization method of shale gas horizontal well perforation cluster spacing
CN106593390A (en) * 2016-12-23 2017-04-26 西南石油大学 Method for simulating horizontal well staged fracturing fracture expanding track
CN106650100A (en) * 2016-12-23 2017-05-10 西南石油大学 Volume alternating fracturing method of horizontal well in experimental shale reservoir
CN107992690A (en) * 2017-12-06 2018-05-04 西南石油大学 A kind of evaluation method of the lower multiple cracking expanded configuration balance degree of induced stress interference
CN109933860A (en) * 2019-02-25 2019-06-25 北京科技大学 Shale gas pressure break horizontal well cluster spacing optimization method
US10677961B1 (en) * 2019-07-08 2020-06-09 Southwest Petroleum University Method for optimizing perforation parameters to maintain uniform fracture growth in multi-stage hydraulic fracturing of horizontal well
CN111255442A (en) * 2020-01-14 2020-06-09 大庆油田有限责任公司 Method for evaluating fracturing fracture by using interference well testing theory

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114707440A (en) * 2022-03-31 2022-07-05 中国石油大学(北京) Multi-crack size calculation method and device and computer equipment

Also Published As

Publication number Publication date
CN113431542B (en) 2023-04-18

Similar Documents

Publication Publication Date Title
CN108009705B (en) Shale reservoir compressibility evaluation method based on support vector machine technology
CN112561144B (en) Method for evaluating main control factor of productivity of tight oil fracturing horizontal well and predicting productivity
CN110644980B (en) Comprehensive classification evaluation method for ultra-low permeability oil reservoir
CN106909758B (en) A kind of method of fine and close oily reservoir-level well multistage sub-clustering perforating site optimization design
CN105735960A (en) Cluster interval optimizing method for segmental multi-cluster fracturing of horizontal well of low-permeability oil and gas reservoir
CN105952427A (en) Method for predicting and evaluating low permeability oil reservoir water-injection induced cracks
CN110359895B (en) Water exploration and drainage method for staged fracturing of horizontal well of heterogeneous huge and thick sandstone
CN115267905B (en) Method for predicting crack lost circulation in drilling engineering of complex structural area
CN116451013B (en) Deep stratum rock in-situ drillability grade value prediction method
CN110439519A (en) A kind of fracturing process and system based on the design of limit current limliting
Liu et al. Natural fractures in tight gas volcanic reservoirs and their influences on production in the Xujiaweizi depression, Songliao Basin, China
CN112377178A (en) Quantitative diagnosis method for boundary water invasion front edge
CN113431542B (en) Method for calculating interference strength of horizontal well fracturing fracture
CN111206921A (en) Description method suitable for favorable reservoir stratum of volcanic overflow phase
CN110671088B (en) New fracture initiation pressure prediction method considering solid-phase plugging main fracture
CN109751018A (en) A kind of construction method for normal pressure shale gas volume fracturing
CN113803044A (en) Method and system for integrally designing unconventional reservoir volume fracturing and well spacing scheme
CN109522579B (en) Method for predicting fracture pressure of horizontal well fracturing construction
CN110764161B (en) Comprehensive method for judging and identifying fractured and broken zone of carbonate rock through drilling data
CN115964836A (en) Method for measuring stress interference intensity among staged multi-cluster fracturing clusters of continental facies shale horizontal well
CN108661629A (en) A kind of engineering dessert quantitative evaluation method for shale formation
CN108979611A (en) A kind of oil-gas reservoir reservoir-level seam horizontal well drilling completion fracturing reform method
CN111335871A (en) Layer-checking hole-repairing technical method based on layered productivity evaluation
CN110872943B (en) Method for determining formation perforation direction
CN113093287B (en) Low-order fault breakpoint identification method

Legal Events

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