CA2785793A1 - Method for controlling a hydraulic fracture path in formations containing natural fractures - Google Patents
Method for controlling a hydraulic fracture path in formations containing natural fractures Download PDFInfo
- Publication number
- CA2785793A1 CA2785793A1 CA2785793A CA2785793A CA2785793A1 CA 2785793 A1 CA2785793 A1 CA 2785793A1 CA 2785793 A CA2785793 A CA 2785793A CA 2785793 A CA2785793 A CA 2785793A CA 2785793 A1 CA2785793 A1 CA 2785793A1
- Authority
- CA
- Canada
- Prior art keywords
- fracture
- parameters
- hydraulic
- formation
- path
- 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.)
- Abandoned
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- 230000015572 biosynthetic process Effects 0.000 title claims description 49
- 238000005755 formation reaction Methods 0.000 title claims description 49
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 230000008569 process Effects 0.000 claims abstract description 15
- 239000012530 fluid Substances 0.000 claims description 27
- 238000002347 injection Methods 0.000 claims description 24
- 239000007924 injection Substances 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000002474 experimental method Methods 0.000 claims description 5
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 230000000977 initiatory effect Effects 0.000 abstract 1
- 206010017076 Fracture Diseases 0.000 description 150
- 208000010392 Bone Fractures Diseases 0.000 description 135
- 238000013461 design Methods 0.000 description 13
- 238000011282 treatment Methods 0.000 description 10
- 230000003993 interaction Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 239000011435 rock Substances 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000004886 process control Methods 0.000 description 2
- 238000011269 treatment regimen Methods 0.000 description 2
- 101100491335 Caenorhabditis elegans mat-2 gene Proteins 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 102100040428 Chitobiosyldiphosphodolichol beta-mannosyltransferase Human genes 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- -1 proppants Substances 0.000 description 1
- 230000004936 stimulating effect Effects 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Saccharide Compounds (AREA)
Abstract
The invention relates to methods for controlling, monitoring and optimizing the parameters of a hydraulic fracture in a stratum. The method comprises the following steps: a) firstly producing a matrix of interdependency between a plurality of starting parameters regarding the stratum, pumping-in operation and fracture and a predicted increase in the trajectory of the hydraulic fracture in the stratum, b) initiating the hydraulic fracturing process of the stratum, c) measuring the actual starting parameters of the hydraulic fracture being formed, d) carrying out exploration and production in accordance with the matrix of the predicted increase in the trajectory of the fracture in dependence on the actual starting parameters of the fracture being formed, e) measuring the actual increase in the trajectory of the fracture being formed during the hydraulic fracturing process, f) comparing the actual increase in the trajectory of the fracture being formed with the predicted increase in the trajectory of the fracture, g) and, in the event of a divergence in said values, changing the actual starting parameters of the fracture being formed.
Description
Method for Controlling a Hydraulic Fracture Path in Formations Containing Natural Fractures Field of the Disclosure [0001] This invention relates to methods for controlling and optimizing parameters of a hydraulic fracture created during the hydraulic fracturing of oil-bearing and gas-bearing reservoirs containing an existing network of natural (geological) fractures, and can be used in appropriate oil and gas fields.
Background of the Disclosure [0002] Hydraulic fracturing is a widely used method for stimulating hydrocarbon inflow from a formation into oil wells or gas wells. To ensure that the best economic result is achieved from the hydraulic fracturing treatment, a design model of the treatment is developed. This model is based on mechanical characteristics of the formation (such as formation stresses, moduli of elasticity and plasticity of the formation, cracking resistance, permeability, etc.), as well as on the selection of optimum parameters for injection of a fracturing fluid into the formation, including the selection of a proper fracturing fluid, proppants, injection conditions, etc. The design model of the fracture plays a rather important role which consists in ensuring that the actual parameters of the fracture geometry will be consistent with the predicted values, and that the selected fracturing fluid and proppant, as well as their quantities, injection rates and the proppant schedule are acceptable for successful implementation of the hydraulic fracturing process.
Background of the Disclosure [0002] Hydraulic fracturing is a widely used method for stimulating hydrocarbon inflow from a formation into oil wells or gas wells. To ensure that the best economic result is achieved from the hydraulic fracturing treatment, a design model of the treatment is developed. This model is based on mechanical characteristics of the formation (such as formation stresses, moduli of elasticity and plasticity of the formation, cracking resistance, permeability, etc.), as well as on the selection of optimum parameters for injection of a fracturing fluid into the formation, including the selection of a proper fracturing fluid, proppants, injection conditions, etc. The design model of the fracture plays a rather important role which consists in ensuring that the actual parameters of the fracture geometry will be consistent with the predicted values, and that the selected fracturing fluid and proppant, as well as their quantities, injection rates and the proppant schedule are acceptable for successful implementation of the hydraulic fracturing process.
[0003] Most of the design models offered nowadays for commercial use (STIMPLAN, NSI Technologies; FracProPT, MFRAC, etc.), as well as the data contained, for example, in Thiercelin, M.J. 2009. Hydraulic fracture propagation in discontinuous media. Proc.: International Conference on Rock Joints and Jointed Rock masses, Tucson, 4-10 Jan 2009. pp. 12, Daneshy, A. (2003) Off-balance growth: A new concept in hydraulic fracturing, Journal of Petroleum Technology, 55, 4, April 2003: 78-85, etc.
are based on the assumption that a single hydraulic fracture plane is created in the formation under treatment. The fracture initiates from the wellbore and increases in length and in height over time as the fracturing fluid and proppant are injected. The in-situ stress condition in the reservoir is such that that there is generally a minimum stress in one of the three stress components, and the created hydraulic fracture tends to propagate in the plane normal to the minimum stress. This assumption about a single planar fracture is usually acceptable for the hydraulic fracturing treatment in a formation consisting of horizontally homogeneous layers.
are based on the assumption that a single hydraulic fracture plane is created in the formation under treatment. The fracture initiates from the wellbore and increases in length and in height over time as the fracturing fluid and proppant are injected. The in-situ stress condition in the reservoir is such that that there is generally a minimum stress in one of the three stress components, and the created hydraulic fracture tends to propagate in the plane normal to the minimum stress. This assumption about a single planar fracture is usually acceptable for the hydraulic fracturing treatment in a formation consisting of horizontally homogeneous layers.
[0004] However, the number of hydraulic fracturing treatments in unconventional oil-bearing formations has recently begun to increase. For complex reservoirs (e.g. for shale-gas reservoirs containing a network of natural (geological) fractures), the assumption about the planar geometry of the hydraulic fracture becomes unjustified. In such formations, the fracturing fluid penetrates so-called "connection branches", thus creating a complex network of crisscross fractures.
[0005] It has become obvious that conventional tools intended for development of optimum hydraulic fracturing treatment strategies and based on planar distribution do not produce the required effect when used in the producing layers having a complex three-dimensional configuration. In such cases, new development tools are required for determination of optimum hydraulic fracturing treatment strategies.
[0006] For example, based on the purpose which consists in conducting a hydraulic fracturing process in formations containing natural (geological) fractures, the closest analogue of the claimed invention is the method of fracturing a naturally fractured rock (WO 2008093264, E 21 B 43/26, published on 07.08.2008). The method is used for conducting a hydraulic fracturing process in formations containing geological fractures. The method comprises: a) acquiring subterranean formation layer geomechanical properties, well completion and reservoir data for the subterranean formation, and a natural fracture network description for the subterranean formation; b) simulating a fracture treatment for the formation, the simulation comprising inputting data acquired into a model which simulates propagation of a network of fracture branches by dividing fracture segments into a plurality of elements to form a fracture grid; c) determining and preparing an optimum fracture fluid composition to achieve the fracturing objective; and, d) injecting the fracturing fluid into a wellbore at a pressure sufficient to fracture the subterranean formation.
[0007] However, the known method has disadvantages consisting in the fact that the proposed method is characterized by a great calculation resource intensity, which increases the capital investments required for implementation of this method. In addition, the method does not include any design model correction to be made by comparing the calculated measurements and the actually obtained measurements, which increases the experimental error and results in a low accuracy of the method on the whole.
Summary of the Disclosure [0008] The claimed invention provides higher efficiency and accuracy of the hydraulic fracturing process control (namely, of the hydraulic fracture propagation path control in formations containing natural (geological) fractures), as well as reduced capital investments required for implementation of this invention.
Summary of the Disclosure [0008] The claimed invention provides higher efficiency and accuracy of the hydraulic fracturing process control (namely, of the hydraulic fracture propagation path control in formations containing natural (geological) fractures), as well as reduced capital investments required for implementation of this invention.
[0009] The method comprises the following operations:
a) creating a matrix of relationship between initial formation, injection and fracture parameters and a predicted increment of a hydraulic fracture path;
b) starting a hydraulic fracturing process;
c) measuring actual initial parameters of a hydraulic fracture being created;
d) using the matrix for retrieving and deriving the predicted increment of the fracture path, depending on the actual initial parameters of the fracture being created;
e) measuring an actual increment of a path of the hydraulic fracture being created;
f) comparing the actual increment of path of the fracture being created with the predicted increment of the fracture path;
g) and, in case of a discrepancy between the actual and predicted increments, changing the actual initial parameters of the fracture being created.
Brief Description of the Figures [0010] Figure 1 illustrates a schematic representation of main parameters in case where a hydraulic fracture approaches an interface.
a) creating a matrix of relationship between initial formation, injection and fracture parameters and a predicted increment of a hydraulic fracture path;
b) starting a hydraulic fracturing process;
c) measuring actual initial parameters of a hydraulic fracture being created;
d) using the matrix for retrieving and deriving the predicted increment of the fracture path, depending on the actual initial parameters of the fracture being created;
e) measuring an actual increment of a path of the hydraulic fracture being created;
f) comparing the actual increment of path of the fracture being created with the predicted increment of the fracture path;
g) and, in case of a discrepancy between the actual and predicted increments, changing the actual initial parameters of the fracture being created.
Brief Description of the Figures [0010] Figure 1 illustrates a schematic representation of main parameters in case where a hydraulic fracture approaches an interface.
[0011] Figure 2 is a schematic representation of possible cases of the hydraulic fracture propagation in case where a hydraulic fracture interacts with an interface.
[0012] Figure 3 I a schematic representation of three possible base cases of the hydraulic fracture propagation in case where a hydraulic fracture interacts with a natural fracture.
[0013] Figure 4 shows a fragment of a matrix of relationship between normalized parameters, and characteristics of the offset of the hydraulic fracture at the interface.
Detailed Description [0014] There are many factors that control the complexity of a hydraulic fracture geometry. In case where the hydraulic fracture interacts with a bedding layer or a natural fracture or a fractional fault (all of them being hereinafter referred to as the "interface"), the hydraulic fracture may expand right through the interface, or may end at the interface, or may continue to propagate at a certain offset along the interface (Figure 3).
Detailed Description [0014] There are many factors that control the complexity of a hydraulic fracture geometry. In case where the hydraulic fracture interacts with a bedding layer or a natural fracture or a fractional fault (all of them being hereinafter referred to as the "interface"), the hydraulic fracture may expand right through the interface, or may end at the interface, or may continue to propagate at a certain offset along the interface (Figure 3).
[0015] The key parameters which control the selection of one of possible fracture propagation paths listed above include, but are not limited to, the following: a fracturing fluid injection rate (Q); a fracturing fluid viscosity ( ); a remote stress difference ((Y1-63); interface properties, such as a friction (cp), an adhesion (C), an angle of approach (P) of the hydraulic fracture to the interface, and a gap (d) between the nearest hydraulic fracture tip and the interface (Figure 1).
[0016] Figure 2 shows possible hydraulic fracture propagation paths in case where a hydraulic fracture intersects an interface. Figure 1 and Figure 2 show the interface as a thin horizontal line which cuts each drawing in half;
"Mat 1" and "Mat 2" shows potentially different materials on each side of the interface, and the fracture approaches it from the bottom of the drawing.
For schematic representation of a possible hydraulic fracture propagation path, these drawings are based on the assumption that the fracture has a planar geometry although the process can be expanded to a three-dimensional model.
"Mat 1" and "Mat 2" shows potentially different materials on each side of the interface, and the fracture approaches it from the bottom of the drawing.
For schematic representation of a possible hydraulic fracture propagation path, these drawings are based on the assumption that the fracture has a planar geometry although the process can be expanded to a three-dimensional model.
[0017] By varying independently each of these control parameters through a series of physical experiments or through a numerical simulation, it is possible to determine the conditions under which different scenarios (cases) of the hydraulic fracture intersecting the interface (i.e. the behavior and the direction of the hydraulic fracture), shown in Figure 2, are most likely to happen. This may require a great number of scenarios (series of experiments), depending on the required resolution. But with a numerical simulation made by using state-of-the-art computers, this is quite feasible and industrially applicable.
[0018] The obtained results can be summarized as a matrix (a database, a reference table or other data files) containing key data which determine whether the hydraulic fracture will tend to be more planar during the fracturing fluid injection in one extreme case, or whether it will tend to branch out into many fractures covering a great reservoir volume in the other extreme case, or whether the propagating hydraulic fracture will stop growing at the interface without further propagation. This matrix can be developed through numerical calculations and/or physical experiments.
[0019] In addition, it is also necessary to obtain information on geometrical complexity of the hydraulic fracture propagation through a network of existing natural fractures. This information can be based on the result of the interaction between the hydraulic fracture and existing natural fractures at the interface. The known result of this interaction with a natural fracture can be superimposed to give the resulting path of the hydraulic fracture propagation through a sequence of natural fractures in a fractured formation.
From now on, the hydraulic fracture path change which considers the hydraulic fracture characteristics (such as the behavior, the geometry and the direction of the fracture) will be regarded as the fracture "increment".
[00201 So, the first step of the proposed method is to find a solution to the problem of the interaction between a hydraulic fracture and an existing natural fracture and to keep this solution for further use.
[00211 Since the result of the fracture propagation through the interface depends on such parameters as: 1) formation and natural fracture parameters, 2) fracturing fluid injection parameters, 3) initial fracture path parameters, we consider that, for the purposes of the claimed invention, the above-mentioned parameters form the "initial" parameters of formation, injection and fracture.
[00221 So, the "initial" parameters include, but are not limited to, the following:
- formation and natural fracture parameters (properties): formation mechanical stresses in this point, coefficients of interface adhesion and coefficients of interface friction, a relative angle of the formation and natural fracture slope relative to the hydraulic fracture, location and size of natural fractures, - fracturing fluid injection parameters (properties): a viscosity of a fracturing fluid to be injected and the fracturing fluid injection rate, or an average fracturing fluid pressure in the hydraulic fracture, - initial fracture path (geometry) parameters (properties): a fracture length, a gap (if any) between a hydraulic fracture tip and an interface, etc.
[0023] The results obtained for each parameter under investigation can be kept in a relevant data file (a reference table, a database or a spreadsheet), i.e. in a matrix of relationship, where a relevant predicted path (that is to say, a hydraulic fracture propagation increment) has been derived and kept for each specific set of initial parameters obtained in the course of numerical calculations or in the course of a physical experiment.
[0024] So, the matrix of relationship shall be developed for an independent set of parameters in order to reduce excessive data volume, database creation time and final retrieval time. After the set of parameters has been determined, it is necessary to specify the optimum resolution of parameter values for each independent parameter.
[0025] A hydraulic fracturing process is then implemented and actual measurements of initial parameters of a fracture being created (including formation and natural fracture parameters, fracturing fluid injection parameters, initial fracture path parameters) are taken.
[0026] For example, a number of hydraulic fracture and formation parameters are measured by using measuring instruments (e.g. sensors distributed over relevant zones of the formation, seismic and acoustic measuring instruments, tiltmeters, etc.). For example, the above-mentioned parameters are determined in the real-time mode, which allows to monitor continuously the hydraulic fracturing treatment and to increase the hydraulic fracturing process control efficiency.
[0027] To be able to use the said matrix of relationship during the implementation of the claimed method, it is necessary to use a programming tool capable of performing parameter change operations and logical data processing operations.
[0028] The design tool known from WO 2008093264, E 21 B 43/26 published on 07.08.2008, can be used as such tool. Such a combined tool can be used for developing, monitoring and controlling the hydraulic fracturing treatment in the real-time mode, as well as for evaluating the results of the hydraulic fracturing treatment in a relevant formation containing natural fractures or faults.
[0029] The design tool shall be capable of accessing the matrix of relationship to determine the location of fracture re-initiation or different fracture propagation path (increment).
[0030] So, first of all, when conducting the hydraulic fracturing, it is necessary to provide a design tool with actual initial parameters which characterize a formation and natural fractures, a fracturing fluid injection and an initial hydraulic fracture path, and which are measured by using sensors or by taking dedicated measurements, as described above. The above-mentioned parameters measured during the hydraulic fracturing process shall also be provided in the real-time mode to the design tool.
[0031] These actually measured parameters shall be then converted into independent parameters (normalized values) because the parameters in the matrix of relationship are represented as normalized (non-dimensional) values. The said design tool is therefore capable of converting the parameters actually measured during the hydraulic fracturing process into normalized values (and vice versa).
[0032] In addition, the design tool is capable of accessing the matrix of relationship to retrieve a direct scenario (path) or a scenario (path) obtained by interpolation between two, or more than two, scenarios which represent the predicted increment of the hydraulic fracture propagation path, depending on the actually measured initial parameters of the hydraulic fracture being created.
[0033] The design tool then determines from the matrix of relationship the resulting fracture propagation path (the predicted fracture increment) along with relevant fracture characteristics (parameters), such as the offset of the fracture re-initiation along the interface on the opposite side of the fracture.
And the normalized parameters are then converted into units of measurable physical quantities (for example, the actual distance of the offset).
[0034] While the fracture grows during the hydraulic fracturing process, the actual increment of the resulting fracture path is measured by using sensors or other measurement methods which allow to determine the required characteristics of the fracture. The said data are provided to the design tool.
[0035] The actual increment of the resulting fracture path is compared with the predicted increment of the fracture path. This comparison operation is performed by a control device.
[0036] In fact, the actually measured data on the hydraulic fracture propagation may differ from the results of the calculations made by the control device. To minimize prediction errors and to improve the final result of the hydraulic fracturing treatment, it is possible to use any discrepancies between the predicted data and the measured data to change the actual initial parameters of the resulting fracture in the near-real-time mode (which certainly allows you to optimize the entire method as a whole).
[0037] In such an operation, certain initial parameters are corrected (changed) during the hydraulic fracturing process. For example, it is possible to change the geometry or the pressure, the injection rate and the viscosity of the fracturing fluid in the near-real-time mode, considering the error resulting from the comparison operation.
[0038] Thus, we improve the consistency between the predicted parameters and the measured parameters during the implementation of the method, which, in turn (in combination with the essential features listed above), results in a higher efficiency of the hydraulic fracture path control.
[0039] So, first, it is necessary to create a matrix of relationship between initial parameters characterizing a formation and natural fractures, a fracturing fluid injection and an initial fracture path before the interaction with a natural fracture, on the one hand, and a predicted increment of a hydraulic fracture path after the interaction with the natural fracture, on the other hand. This matrix can be developed, for instance, by finding a numerical solution to the problem of the mechanical interaction between the hydraulic fracture at a constant internal pressure and the natural fracture at the moment of their contact, and can be represented in the form of a table as shown in Figure 4.
[0040] To reduce the number of independent parameters of the problem, it is necessary to normalize all values in the numerical problem, so the solution represents the function of the coefficient of friction on the natural fracture A, the angle of slope of the natural fracture 6, the non-dimensional formation differential stress o = (6, - 63)/(6, + 63) and the non-dimensional hydraulic fracture excess pressure H = 2(p - a3)'(6, +a3)' Considering that, having come across the natural fracture, the hydraulic fracture will continue to propagate if a new tensile crack initiates on the other side of the natural fracture, the offset of the tensile stress peak along the natural fracture and the stress peak itself are then selected as the characteristics of the hydraulic fracture path increment at the natural fracture. As a result of the normalization, the offset will be normalized to the hydraulic fracture length LHF and the stress will be normalized to the average formation stress 6m =(0l+63)/2.
[0041] Then, hydraulic fracturing process is started and actual initial formation, injection and hydraulic fracture parameters are measured. At a certain stage, it is necessary to learn how a hydraulic fracture path will change after a hydraulic fracture has interacted with a natural fracture crossing its path. The measured formation stresses are equal to U, = 6 MPa, 63 =4 MPa, a hydraulic fracture pressure is estimated as p = 5.5 MPa, and a current hydraulic fracture length is equal to LHF = 100 in. The natural fracture has a friction coefficient of 0.5. Let us then assume two possible orientations of the natural fracture with respect to the hydraulic fracture:
a) and b) 40 .
[0042] By using the normalization rules from the table, we obtain: A = 0.5, 0 = 0.2, U = 0.3, and the angles for the first case (a): 6 = 10 , for the second case (b): Q = 40 .
[0043] The matrix of relationship is then used for retrieving and deriving a predicted increment of the fracture path, in particular, an offset value and a stress peak value, depending on the given parameters. Based on the above-mentioned parameters, we retrieve the values of the hydraulic fracture path increment characteristics to find that the offset of the fracture is equal to 0.049 and the stress is equal to -0.0163 in the first case. In the second case, the offset of the fracture is equal to 0.025 and the stress is equal to -0.7929 (the relevant lines of the table shown in Figure 4 are highlighted).
[0044] The normalized values are then converted into dimensional units of relevant characteristics. We obtain for case a): offset = 4.9 in and stress =
0.08 MPa, for case b): offset = 2.5 in and stress = 3.96 MPa.
[0045] Then, based on the comparison of the tensile stress peak with the tensile strength of the rock, a conclusion is made as to whether a secondary fracture will or will not be formed at the natural fracture. Assuming that the tensile strength is equal to 4 MPa, the hydraulic fracture will intersect the natural fracture at a distance of 2.5 m in case a). In case b), no intersection will take place and the hydraulic fracture will stop at this natural fracture.
[0046] Let us then assume that the hydraulic fracture path increment characteristics (namely, the intersection and the offset of the fracture in case a)) have been measured, and these measurements show that the offset is actually equal to 1.3 in.
[0047] The actual increment of the resulting fracture path is then compared with the predicted increment of the fracture path. As these values turned out to be different, the actual initial parameters of the resulting fracture shall be changed. In this case, the fluid pressure value shall be corrected. Having retrieved the required data from the matrix, we find that this value corresponds to a non-dimensional pressure of 0.2 for the given hydraulic fracture length. Conversion into dimensional units gives a new well pressure value of 5 MPa which will be used as a more accurate value in subsequent applications of the method.
[0048] The above embodiments of the invention are made for illustrative purposes only because the invention can be modified and can be practically used in a variety of similar ways which are obvious to persons who are skilled in this art and who can make use of the advantages offered by the concepts described in this application. Moreover, the applicant had no intention to impose any limitations on the details of the designs or developments shown herein, except for those described in the Claims below.
It is therefore obvious that the specific embodiments of the invention, disclosed above, can be changed or modified, and all such modifications shall be considered to be protected in accordance with the scope and the nature of the invention.
From now on, the hydraulic fracture path change which considers the hydraulic fracture characteristics (such as the behavior, the geometry and the direction of the fracture) will be regarded as the fracture "increment".
[00201 So, the first step of the proposed method is to find a solution to the problem of the interaction between a hydraulic fracture and an existing natural fracture and to keep this solution for further use.
[00211 Since the result of the fracture propagation through the interface depends on such parameters as: 1) formation and natural fracture parameters, 2) fracturing fluid injection parameters, 3) initial fracture path parameters, we consider that, for the purposes of the claimed invention, the above-mentioned parameters form the "initial" parameters of formation, injection and fracture.
[00221 So, the "initial" parameters include, but are not limited to, the following:
- formation and natural fracture parameters (properties): formation mechanical stresses in this point, coefficients of interface adhesion and coefficients of interface friction, a relative angle of the formation and natural fracture slope relative to the hydraulic fracture, location and size of natural fractures, - fracturing fluid injection parameters (properties): a viscosity of a fracturing fluid to be injected and the fracturing fluid injection rate, or an average fracturing fluid pressure in the hydraulic fracture, - initial fracture path (geometry) parameters (properties): a fracture length, a gap (if any) between a hydraulic fracture tip and an interface, etc.
[0023] The results obtained for each parameter under investigation can be kept in a relevant data file (a reference table, a database or a spreadsheet), i.e. in a matrix of relationship, where a relevant predicted path (that is to say, a hydraulic fracture propagation increment) has been derived and kept for each specific set of initial parameters obtained in the course of numerical calculations or in the course of a physical experiment.
[0024] So, the matrix of relationship shall be developed for an independent set of parameters in order to reduce excessive data volume, database creation time and final retrieval time. After the set of parameters has been determined, it is necessary to specify the optimum resolution of parameter values for each independent parameter.
[0025] A hydraulic fracturing process is then implemented and actual measurements of initial parameters of a fracture being created (including formation and natural fracture parameters, fracturing fluid injection parameters, initial fracture path parameters) are taken.
[0026] For example, a number of hydraulic fracture and formation parameters are measured by using measuring instruments (e.g. sensors distributed over relevant zones of the formation, seismic and acoustic measuring instruments, tiltmeters, etc.). For example, the above-mentioned parameters are determined in the real-time mode, which allows to monitor continuously the hydraulic fracturing treatment and to increase the hydraulic fracturing process control efficiency.
[0027] To be able to use the said matrix of relationship during the implementation of the claimed method, it is necessary to use a programming tool capable of performing parameter change operations and logical data processing operations.
[0028] The design tool known from WO 2008093264, E 21 B 43/26 published on 07.08.2008, can be used as such tool. Such a combined tool can be used for developing, monitoring and controlling the hydraulic fracturing treatment in the real-time mode, as well as for evaluating the results of the hydraulic fracturing treatment in a relevant formation containing natural fractures or faults.
[0029] The design tool shall be capable of accessing the matrix of relationship to determine the location of fracture re-initiation or different fracture propagation path (increment).
[0030] So, first of all, when conducting the hydraulic fracturing, it is necessary to provide a design tool with actual initial parameters which characterize a formation and natural fractures, a fracturing fluid injection and an initial hydraulic fracture path, and which are measured by using sensors or by taking dedicated measurements, as described above. The above-mentioned parameters measured during the hydraulic fracturing process shall also be provided in the real-time mode to the design tool.
[0031] These actually measured parameters shall be then converted into independent parameters (normalized values) because the parameters in the matrix of relationship are represented as normalized (non-dimensional) values. The said design tool is therefore capable of converting the parameters actually measured during the hydraulic fracturing process into normalized values (and vice versa).
[0032] In addition, the design tool is capable of accessing the matrix of relationship to retrieve a direct scenario (path) or a scenario (path) obtained by interpolation between two, or more than two, scenarios which represent the predicted increment of the hydraulic fracture propagation path, depending on the actually measured initial parameters of the hydraulic fracture being created.
[0033] The design tool then determines from the matrix of relationship the resulting fracture propagation path (the predicted fracture increment) along with relevant fracture characteristics (parameters), such as the offset of the fracture re-initiation along the interface on the opposite side of the fracture.
And the normalized parameters are then converted into units of measurable physical quantities (for example, the actual distance of the offset).
[0034] While the fracture grows during the hydraulic fracturing process, the actual increment of the resulting fracture path is measured by using sensors or other measurement methods which allow to determine the required characteristics of the fracture. The said data are provided to the design tool.
[0035] The actual increment of the resulting fracture path is compared with the predicted increment of the fracture path. This comparison operation is performed by a control device.
[0036] In fact, the actually measured data on the hydraulic fracture propagation may differ from the results of the calculations made by the control device. To minimize prediction errors and to improve the final result of the hydraulic fracturing treatment, it is possible to use any discrepancies between the predicted data and the measured data to change the actual initial parameters of the resulting fracture in the near-real-time mode (which certainly allows you to optimize the entire method as a whole).
[0037] In such an operation, certain initial parameters are corrected (changed) during the hydraulic fracturing process. For example, it is possible to change the geometry or the pressure, the injection rate and the viscosity of the fracturing fluid in the near-real-time mode, considering the error resulting from the comparison operation.
[0038] Thus, we improve the consistency between the predicted parameters and the measured parameters during the implementation of the method, which, in turn (in combination with the essential features listed above), results in a higher efficiency of the hydraulic fracture path control.
[0039] So, first, it is necessary to create a matrix of relationship between initial parameters characterizing a formation and natural fractures, a fracturing fluid injection and an initial fracture path before the interaction with a natural fracture, on the one hand, and a predicted increment of a hydraulic fracture path after the interaction with the natural fracture, on the other hand. This matrix can be developed, for instance, by finding a numerical solution to the problem of the mechanical interaction between the hydraulic fracture at a constant internal pressure and the natural fracture at the moment of their contact, and can be represented in the form of a table as shown in Figure 4.
[0040] To reduce the number of independent parameters of the problem, it is necessary to normalize all values in the numerical problem, so the solution represents the function of the coefficient of friction on the natural fracture A, the angle of slope of the natural fracture 6, the non-dimensional formation differential stress o = (6, - 63)/(6, + 63) and the non-dimensional hydraulic fracture excess pressure H = 2(p - a3)'(6, +a3)' Considering that, having come across the natural fracture, the hydraulic fracture will continue to propagate if a new tensile crack initiates on the other side of the natural fracture, the offset of the tensile stress peak along the natural fracture and the stress peak itself are then selected as the characteristics of the hydraulic fracture path increment at the natural fracture. As a result of the normalization, the offset will be normalized to the hydraulic fracture length LHF and the stress will be normalized to the average formation stress 6m =(0l+63)/2.
[0041] Then, hydraulic fracturing process is started and actual initial formation, injection and hydraulic fracture parameters are measured. At a certain stage, it is necessary to learn how a hydraulic fracture path will change after a hydraulic fracture has interacted with a natural fracture crossing its path. The measured formation stresses are equal to U, = 6 MPa, 63 =4 MPa, a hydraulic fracture pressure is estimated as p = 5.5 MPa, and a current hydraulic fracture length is equal to LHF = 100 in. The natural fracture has a friction coefficient of 0.5. Let us then assume two possible orientations of the natural fracture with respect to the hydraulic fracture:
a) and b) 40 .
[0042] By using the normalization rules from the table, we obtain: A = 0.5, 0 = 0.2, U = 0.3, and the angles for the first case (a): 6 = 10 , for the second case (b): Q = 40 .
[0043] The matrix of relationship is then used for retrieving and deriving a predicted increment of the fracture path, in particular, an offset value and a stress peak value, depending on the given parameters. Based on the above-mentioned parameters, we retrieve the values of the hydraulic fracture path increment characteristics to find that the offset of the fracture is equal to 0.049 and the stress is equal to -0.0163 in the first case. In the second case, the offset of the fracture is equal to 0.025 and the stress is equal to -0.7929 (the relevant lines of the table shown in Figure 4 are highlighted).
[0044] The normalized values are then converted into dimensional units of relevant characteristics. We obtain for case a): offset = 4.9 in and stress =
0.08 MPa, for case b): offset = 2.5 in and stress = 3.96 MPa.
[0045] Then, based on the comparison of the tensile stress peak with the tensile strength of the rock, a conclusion is made as to whether a secondary fracture will or will not be formed at the natural fracture. Assuming that the tensile strength is equal to 4 MPa, the hydraulic fracture will intersect the natural fracture at a distance of 2.5 m in case a). In case b), no intersection will take place and the hydraulic fracture will stop at this natural fracture.
[0046] Let us then assume that the hydraulic fracture path increment characteristics (namely, the intersection and the offset of the fracture in case a)) have been measured, and these measurements show that the offset is actually equal to 1.3 in.
[0047] The actual increment of the resulting fracture path is then compared with the predicted increment of the fracture path. As these values turned out to be different, the actual initial parameters of the resulting fracture shall be changed. In this case, the fluid pressure value shall be corrected. Having retrieved the required data from the matrix, we find that this value corresponds to a non-dimensional pressure of 0.2 for the given hydraulic fracture length. Conversion into dimensional units gives a new well pressure value of 5 MPa which will be used as a more accurate value in subsequent applications of the method.
[0048] The above embodiments of the invention are made for illustrative purposes only because the invention can be modified and can be practically used in a variety of similar ways which are obvious to persons who are skilled in this art and who can make use of the advantages offered by the concepts described in this application. Moreover, the applicant had no intention to impose any limitations on the details of the designs or developments shown herein, except for those described in the Claims below.
It is therefore obvious that the specific embodiments of the invention, disclosed above, can be changed or modified, and all such modifications shall be considered to be protected in accordance with the scope and the nature of the invention.
Claims (20)
1. Method for controlling a hydraulic fracture path in formations containing natural fractures comprising:
a) creating a matrix of relationship between initial formation, injection and fracture parameters and a predicted increment of a hydraulic fracture path, b) starting a hydraulic fracturing process, c) measuring actual initial parameters of a hydraulic fracture being created, d) using the matrix for retrieving and deriving the predicted increment of the fracture path, depending on the actual initial parameters of the fracture being created, e) measuring an actual increment of a path of the hydraulic fracture being created , f) comparing the actual increment of the path of the fracture being created with the predicted increment of the fracture path, and g) in case of a discrepancy between the actual and the predicted increments changing the actual initial parameters of the fracture being created.
a) creating a matrix of relationship between initial formation, injection and fracture parameters and a predicted increment of a hydraulic fracture path, b) starting a hydraulic fracturing process, c) measuring actual initial parameters of a hydraulic fracture being created, d) using the matrix for retrieving and deriving the predicted increment of the fracture path, depending on the actual initial parameters of the fracture being created, e) measuring an actual increment of a path of the hydraulic fracture being created , f) comparing the actual increment of the path of the fracture being created with the predicted increment of the fracture path, and g) in case of a discrepancy between the actual and the predicted increments changing the actual initial parameters of the fracture being created.
2. Method of Claim 1, wherein the matrix of relationship between initial formation, injection and fracture parameters and the predicted increment of the hydraulic fracture path is created through numerical calculations or experiments or both.
3. Method of Claim 1, wherein the initial formation, injection and fracture parameters include formation and natural fracture parameters, fracturing fluid injection parameters and initial fracture path parameters.
4. Method of Claim 3, wherein formation mechanical stresses are used as the formation and natural fracture parameters.
5. Method of Claim 3, wherein coefficients of interface adhesion are used as the formation and natural fracture parameters.
6. Method of Claim 3, wherein coefficients of interface friction are used as the formation and natural fracture parameters.
7. Method of Claim 3, wherein a relative angle between the hydraulic fracture and a natural fracture in a point of their contact is used as the formation and natural fracture parameters.
8. Method of Claim 3, wherein a parameter characterizing the location of natural fractures is used as the formation and natural fracture parameters.
9. Method of Claim 3, wherein a natural fracture size parameter is used as the formation and natural fracture parameters.
10. Method of Claim 3, wherein a viscosity of a fracturing fluid to be injected is used as the fracturing fluid injection parameter.
11. Method of Claim 3, wherein a fracturing fluid injection rate is used as the fracturing fluid injection parameter.
12. Method of Claim 3, wherein an average fracturing fluid pressure in the hydraulic fracture is used as the fracturing fluid injection parameter.
13. Method of Claim 3, wherein a fracture length is used as the initial fracture path parameter.
14. Method of Claim 3, wherein a gap between a hydraulic fracture tip and an interface is used as the initial fracture path parameter.
15. Method of Claim 1, wherein the initial formation, injection and fracture parameters and the predicted increment of the hydraulic fracture path are represented in the matrix as normalized values.
16. Method of Claim 1, wherein operations c), d), e), f), g) are performed in the real-time mode of the hydraulic fracturing process.
17 17. Method of Claim 1, wherein operations c), d), e), f), g) are performed by using a control device.
18. Method of Claim 17, wherein the control device is capable of converting the actually measured values into normalized values, and vice versa.
19. Method of Claim 1, wherein the actual initial parameters of the hydraulic fracture being created are measured by using measuring techniques.
20. Method of Claim 1, wherein the actual increment of the hydraulic fracture being created is measured by using measuring techniques.
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/RU2009/000746 WO2011081544A1 (en) | 2009-12-30 | 2009-12-30 | Method for controlling the trajectory of a hydraulic fracture in strata containing natural fractures |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CA2785793A1 true CA2785793A1 (en) | 2011-07-07 |
Family
ID=44226678
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA2785793A Abandoned CA2785793A1 (en) | 2009-12-30 | 2009-12-30 | Method for controlling a hydraulic fracture path in formations containing natural fractures |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US20130319657A1 (en) |
| CA (1) | CA2785793A1 (en) |
| MX (1) | MX2012007687A (en) |
| RU (1) | RU2505670C1 (en) |
| WO (1) | WO2011081544A1 (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2766748A4 (en) * | 2011-10-11 | 2016-07-13 | Services Petroliers Schlumberger | System and method for performing stimulation operations |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| BR112016028422B1 (en) | 2014-06-05 | 2022-04-19 | Geoquest Systems B.V. | Method for hydraulic fracturing of an underground formation crossed by an exploration well |
| RU2585296C1 (en) * | 2015-03-27 | 2016-05-27 | Открытое акционерное общество "Нефтяная компания "Роснефть" | Method of determining drained hydraulic fracturing crack width and degree of sedimentation of proppant therein |
| CN108386176B (en) * | 2018-01-30 | 2020-07-14 | 中国石油天然气股份有限公司 | Physical model test method for natural crack and artificial crack extension rule |
| CA3037543C (en) | 2018-03-21 | 2023-09-26 | ResFrac Corporation | Systems and methods for hydraulic fracture and reservoir simulation |
| CN110424939A (en) * | 2019-08-12 | 2019-11-08 | 西南石油大学 | A method of increasing gneiss oil-gas reservoir and stitches net volume fracturing effect |
| CN111963149B (en) * | 2020-08-26 | 2024-03-01 | 中国石油天然气集团有限公司 | Post-fracturing stratum pressure solving method taking earth stagnation amount pressurization into consideration |
| CN113341465B (en) * | 2021-06-11 | 2023-05-09 | 中国石油大学(北京) | Directional anisotropic medium ground stress prediction method, device, medium and equipment |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7509245B2 (en) * | 1999-04-29 | 2009-03-24 | Schlumberger Technology Corporation | Method system and program storage device for simulating a multilayer reservoir and partially active elements in a hydraulic fracturing simulator |
| CA2475007A1 (en) * | 2002-02-01 | 2003-08-14 | Regents Of The University Of Minnesota | Interpretation and design of hydraulic fracturing treatments |
| US6991037B2 (en) * | 2003-12-30 | 2006-01-31 | Geosierra Llc | Multiple azimuth control of vertical hydraulic fractures in unconsolidated and weakly cemented sediments |
| RU2327154C2 (en) * | 2004-04-23 | 2008-06-20 | Шлюмберже Текнолоджи Б.В | Method and system for monitoring of cavities filled with liquid in the medium on the basis of boundary waves that are distributed on their surfaces |
| US7386431B2 (en) * | 2005-03-31 | 2008-06-10 | Schlumberger Technology Corporation | Method system and program storage device for simulating interfacial slip in a hydraulic fracturing simulator software |
| US8412500B2 (en) * | 2007-01-29 | 2013-04-02 | Schlumberger Technology Corporation | Simulations for hydraulic fracturing treatments and methods of fracturing naturally fractured formation |
-
2009
- 2009-12-30 US US13/519,572 patent/US20130319657A1/en not_active Abandoned
- 2009-12-30 CA CA2785793A patent/CA2785793A1/en not_active Abandoned
- 2009-12-30 WO PCT/RU2009/000746 patent/WO2011081544A1/en active Application Filing
- 2009-12-30 MX MX2012007687A patent/MX2012007687A/en active IP Right Grant
- 2009-12-30 RU RU2012130906/03A patent/RU2505670C1/en not_active IP Right Cessation
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2766748A4 (en) * | 2011-10-11 | 2016-07-13 | Services Petroliers Schlumberger | System and method for performing stimulation operations |
| US10605060B2 (en) | 2011-10-11 | 2020-03-31 | Schlumberger Technology Corporation | System and method for performing stimulation operations |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2011081544A1 (en) | 2011-07-07 |
| US20130319657A1 (en) | 2013-12-05 |
| RU2505670C1 (en) | 2014-01-27 |
| MX2012007687A (en) | 2012-07-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20130319657A1 (en) | Method for controlling the trajectory of a hydraulic fracture in strata-containing natural fractures | |
| Cipolla et al. | New algorithms and integrated workflow for tight gas and shale completions | |
| CA2873520C (en) | System and method for performing downhole stimulation operations | |
| CA2852044C (en) | System and method for performing stimulation operations | |
| Dohmen et al. | Measurement and analysis of 3D stress shadowing related to the spacing of hydraulic fracturing in unconventional reservoirs | |
| US8775141B2 (en) | System and method for performing oilfield simulation operations | |
| US8794316B2 (en) | Refracture-candidate evaluation and stimulation methods | |
| CA2926788C (en) | Designing wellbore completion intervals | |
| US11753918B2 (en) | Method for multilayer hydraulic fracturing treatment with real-time adjusting | |
| Desroches et al. | Benefits of controlled hydraulic fracture placement: theory and field experiment | |
| Dumitrache et al. | A multidisciplinary approach to production optimization through hydraulic fracturing stimulation and geomechanical modelling in Clair Field | |
| Ramakrishnan et al. | Real-time completion optimization of multiple laterals in gas shale reservoirs: Integration of geology, log, surface seismic, and microseismic information | |
| Gulrajani et al. | Pressure-history inversion for interpretation of fracture treatments | |
| Ejofodomi et al. | Improving well completion via real-time microseismic monitoring: a West Texas case study | |
| Agudelo et al. | Casabe: Water injection optimization and surveillance in a mature field | |
| STRATA-CONTAINING et al. | 12, Patent Application Publication o Pub. No.: US 2013/0319657A1 | |
| Mishra et al. | High Accuracy estimation of hydraulic fracture geometry using crosswell electromagnetics | |
| Malik et al. | Stress Profiling with Microfracturing and Sonic Logs in Antelope Shale | |
| Butula et al. | Re-fracturing considerations of horizontal well multi stage fractured completions in mid permeability formations | |
| Lehman et al. | Calibrating Fracture Models with Direct Diagnostics: A Necessary but Humbling Experience | |
| Pandey et al. | Estimating Fracture Containment in Layered Formations Using 3D Fracture Growth Assumptions | |
| Edwards et al. | New Frac-Pack Redesign Methodology Improves Estimation of Tip Screenout Post-FET | |
| Modeland et al. | Applying petrophysics and systematic completion science to increase potential: A Terryville field tight gas case history | |
| EA042125B1 (en) | MULTI-STAGE HYDRAULIC FRACTURING TREATMENT WITH REAL-TIME ADJUSTMENT | |
| Kundert et al. | Multidiscipline Approach to Real-Time Optimization of Down-Spaced Pay Adds in a Mature Gas Play, San Juan Basin, New Mexico |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FZDE | Discontinued |
Effective date: 20141230 |