US20120325462A1 - Method for Determining Spacing of Hydraulic Fractures in a Rock Formation - Google Patents
Method for Determining Spacing of Hydraulic Fractures in a Rock Formation Download PDFInfo
- Publication number
- US20120325462A1 US20120325462A1 US13/531,031 US201213531031A US2012325462A1 US 20120325462 A1 US20120325462 A1 US 20120325462A1 US 201213531031 A US201213531031 A US 201213531031A US 2012325462 A1 US2012325462 A1 US 2012325462A1
- Authority
- US
- United States
- Prior art keywords
- fracture
- fractures
- spacing
- determining
- stress
- 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
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 75
- 239000011435 rock Substances 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000001902 propagating effect Effects 0.000 claims description 25
- 238000011065 in-situ storage Methods 0.000 claims description 15
- 230000000694 effects Effects 0.000 claims description 5
- 230000000644 propagated effect Effects 0.000 claims description 5
- 230000003247 decreasing effect Effects 0.000 claims description 4
- 230000001939 inductive effect Effects 0.000 claims 1
- 206010017076 Fracture Diseases 0.000 description 404
- 208000010392 Bone Fractures Diseases 0.000 description 279
- 238000005755 formation reaction Methods 0.000 description 70
- 229930195733 hydrocarbon Natural products 0.000 description 18
- 150000002430 hydrocarbons Chemical class 0.000 description 18
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 18
- XQCFHQBGMWUEMY-ZPUQHVIOSA-N Nitrovin Chemical compound C=1C=C([N+]([O-])=O)OC=1\C=C\C(=NNC(=N)N)\C=C\C1=CC=C([N+]([O-])=O)O1 XQCFHQBGMWUEMY-ZPUQHVIOSA-N 0.000 description 17
- 239000007789 gas Substances 0.000 description 17
- 238000004519 manufacturing process Methods 0.000 description 10
- 239000004215 Carbon black (E152) Substances 0.000 description 9
- 230000006870 function Effects 0.000 description 8
- 239000003345 natural gas Substances 0.000 description 8
- 238000005553 drilling Methods 0.000 description 6
- 230000000977 initiatory effect Effects 0.000 description 6
- 238000013459 approach Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 239000004576 sand Substances 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 239000011368 organic material Substances 0.000 description 3
- 238000007792 addition Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000011282 treatment Methods 0.000 description 2
- 230000005483 Hooke's law Effects 0.000 description 1
- 235000019738 Limestone Nutrition 0.000 description 1
- 208000006670 Multiple fractures Diseases 0.000 description 1
- 208000002565 Open Fractures Diseases 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 239000004058 oil shale Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000000518 rheometry Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Definitions
- the present disclosure relates in general to well drilling and, more particularly, to a method for determining spacing of hydraulic fractures in a rock formation.
- Hydrocarbon (e.g., oil, natural gas, etc.) reservoirs may be found in geologic formations that have little to no porosity (e.g., shale, sandstone etc.).
- the hydrocarbons may be trapped within fractures and pore spaces of the formation. Additionally, the hydrocarbons may be adsorbed onto organic material of the shale formation.
- the rapid development of extracting hydrocarbons from these unconventional reservoirs can be tied to the combination of horizontal drilling and hydraulic fracturing (“fracing”) of the formations. Horizontal drilling has allowed for drilling along and within hydrocarbon reservoirs of a formation to better capture the hydrocarbons trapped within the reservoirs.
- hydrocarbons may be captured by increasing the number of fractures in the formation and/or increasing the size of already present fractures through fracing.
- the spacing between fractures as well as the ability to stimulate the fractures naturally present in the rock may be major factors in the success of horizontal completions in unconventional hydrocarbon reservoirs.
- a method comprising determining an expected trajectory of induced fractures, analyzing net pressure associated with the induced fractures, and determining at least one of spacing of induced fractures and a property of the induced fractures based on the net pressure.
- Computer-readable medium containing the same are also disclosed.
- a method of optimizing fracture spacing includes propagating an initial fracture, measuring pressure associated with propagating the initial fracture, determining a minimum spacing required to prevent a second fracture from intersecting the initial fracture, and propagating the second fracture at least the minimum spacing distance away from the initial fracture.
- a method of optimizing fracture spacing includes analyzing stresses associated with a first set of at least one fractures associated with a first well, analyzing stresses associated with a second set of at least one fractures associated with a second well, and determining spacing of a fracture associated with a third well such that the fracture associated with the third well does not intersect with the first set and the second set of fractures, the third well running between the first and second wells.
- a method of determining maximum horizontal pressure includes measuring an actual pressure during each stage of fracturing of a rock formation, determining a theoretical expected pressure during each stage of fracturing of the rock formation, and determining a maximum horizontal pressure of the rock formation based at least on a comparison of the theoretical expected pressure and the measured actual pressure.
- FIG. 1 illustrates an example schematic of a gas well configured to extract natural gas from a gas rich shale formation, according to some embodiments of the present disclosure
- FIGS. 2 a and 2 b illustrate examples of the reorientation of stresses in a rock formation due to the placement of a fracture orthogonal to a horizontal well according to some embodiments of the present disclosure
- FIG. 3 illustrates the geometry of a single transverse fracture of a shale formation that includes a pay zone that may include hydrocarbons (e.g., natural gas) and bounding layers that may bound the pay zone, according to some embodiments of the present disclosure
- FIG. 4 illustrates a three dimensional model of multiple transverse fractures in a layered rock formation according to some embodiments of the present disclosure
- FIG. 5 illustrates the calculated propagation of a subsequent fracture (n+1) based on the mechanical stress interference of a previous fracture (n), according to some embodiments of the present disclosure
- FIG. 6 illustrates the results of calculating fracture propagation with each fracture spaced approximately 400 feet apart, according to some embodiments of the present disclosure
- FIG. 7 illustrates the results of calculating fracture propagation with each fracture spaced approximately 300 feet apart, according to some embodiments of the present disclosure
- FIG. 8 illustrates the stress distribution of a rock formation with the stress distribution being influenced by the propagation of a fracture produced during a fourth stage of a fracture treatment, according to some embodiments of the present disclosure
- FIG. 9 illustrates fracture propagation with the fracture spacing reduced to 250 ft., according to some embodiments of the present disclosure
- FIG. 10 illustrates the stress distribution of a rock formation caused by fractures of FIG. 9 , according to some embodiments of the present disclosure
- FIG. 11 illustrates fracture propagation with the fracture spacing reduced to 200 ft., according to some embodiments of the present disclosure
- FIG. 12 illustrates fracture propagation with the fracture spacing reduced to 150 ft., according to some embodiments of the present disclosure
- FIGS. 13 a and 13 b illustrate the impact of fracture spacing on the angle of deviation of the fractures from the orthogonal path, according to some embodiments of the present disclosure
- FIG. 14 illustrates the impact of fracture spacing on the evolution of the net closure stress, according to some embodiments of the present disclosure
- FIGS. 15 a and 15 b illustrate the horizontal stress of a rock formation with a stress reversal region associated with a fracture, according to some embodiments of the present disclosure
- FIGS. 16 a and 16 b illustrate the differences between performing consecutive fracturing and alternate fracturing, according to some embodiments of the present disclosure
- FIGS. 17 a and 17 b illustrate the stress orientation of a rock formation with stress reversal regions associated with a first fracture and a second fracture, according to some embodiments of the present disclosure
- FIG. 18 shows a spacing of fractures for which the stress contrast may be lowest, according to some embodiments of the present disclosure
- FIG. 19 illustrates that the deviatoric stress may approach zero in a near wellbore region in the case of optimum spacing in an alternate fracturing sequence, according to some embodiments of the present disclosure
- FIG. 20 illustrates an example of fracture spacing that may be done with multiple horizontal lateral wells, according to some embodiments of the present disclosure
- FIGS. 21 a and 21 b illustrate the stress distribution between two pairs of fractures propagated from outside lateral wells, according to some embodiments of the present disclosure
- FIG. 22 illustrates the relationship between the length of a “middle fracture” propagating from a center lateral well for different values of the fracture length with respect to the inter-well spacing, according to some embodiments of the present disclosure.
- FIG. 23 illustrates the local stress contrast that may be recorded along the assumed propagation direction of a middle fracture, according to some embodiments of the present disclosure.
- FIG. 1 illustrates an example schematic of a gas well 100 configured to extract natural gas from a gas rich shale formation 102 .
- Well 100 may be drilled using horizontal drilling methods to create a wellbore 104 that runs within and along shale formation 102 .
- wellbore 104 may be drilled such that wellbore 104 runs perpendicular to the maximum horizontal in-situ stresses of shale formation 102 to obtain better production of well 100 .
- Shale formation 102 may produce natural gas that is trapped in fractures and pore spaces of shale formation 102 .
- the natural gas may also be adsorbed in organic material included in the shale of shale formation 102 .
- well bore 104 runs through shale formation 102
- well bore 104 may also run through fractures (not expressly shown) of shale formation 102 .
- the gas in the fractures may enter well bore 104 and may accordingly be retrieved at a drilling rig 106 of well 100 .
- the gas adsorbed on the organic material may be released into the fractures such that the adsorbed gas may also be retrieved.
- the number and/or size of fractures in shale formation 102 may be increased using hydraulic fracturing (“fracing”).
- Fracing may refer to any process used to initiate and propagate a fracture in a rock formation. Additionally fracing may be used to increase existing fractures in a rock formation. Fracing may include forcing a hydraulic fluid in a fracture of a rock formation to increase the size of the fracture and introducing proppant (e.g., sand) in the newly induced fracture to keep the fracture open.
- proppant e.g., sand
- the fracture may be an existing fracture in the formation, or may be initiated using a variety of techniques known in the art.
- the amount of pressure needed to extend and propagate the fracture may be referred to as the “fracturing pressure.”
- producing fractures during fracing may change the stress properties of a rock formation. Accordingly, subsequent transverse fractures initiated from a horizontal well may deviate toward or away from the previous fracture depending on the stress reorientation caused by the fracing.
- the stress reorientation may be a function of mechanical properties of the reservoir rock, fracture spacing, and the orientation of the previous fracture. As described in further detail below, in some instances spacing frac treatments too close together may result in a subsequent fracture intersecting with a previous fracture. Therefore, in such instances, the contribution of the subsequent fracture in hydrocarbon production may be reduced or minimized.
- the spacing of performing fracturing operations for a hydrocarbon well may be determined using net pressure measurements to determine a minimum frac spacing that also reduces the likelihood of subsequent fractures intersecting and interfering with previous fractures.
- net pressure may be determined by surface pressures or down-hole pressures during fracturing. Therefore, the selection of fracture spacing may be such that the number of fractures initiated from a horizontal wellbore may be increased while also reducing the likelihood that the fractures may interfere and/or intersect with each other to allow for better production rates and depletion of hydrocarbons with each induced fracture. Therefore, the economic efficiency of fracturing may be improved and the cost of retrieving hydrocarbons from tight rock formations (e.g., shale formation 102 ) may be reduced.
- FIG. 1 is described as performing fracing with respect to a shale gas formation, the present disclosure may be used to improve frac spacing for any suitable formation (e.g., a tight sand formation, coalbed methane, sandstone, limestone, oil shale).
- any suitable formation e.g., a tight sand formation, coalbed methane, sandstone, limestone, oil shale.
- well 100 is described as being used to extract natural gas, it is understood that the principles described herein may be used to extract any other suitable hydrocarbon.
- FIGS. 2 a and 2 b illustrate examples of the reorientation of stresses in a rock formation due to the placement of a fracture orthogonal (or transverse) to a horizontal well according to some embodiments of the present disclosure.
- the opening of a propped transverse fracture in horizontal wells through hydraulic fracturing may cause a reorientation of stresses in the rock formation surrounding the fracture.
- FIG. 2 a may represent the stress of the rock formation in a horizontal plane (e.g., a plane substantially parallel to the ground. Accordingly, the vertical axis of FIG. 2 a may represent the distance in the x direction from the center of a substantially horizontal wellbore 204 .
- the horizontal axis of FIG. 2 a may represent the distance in the y direction from the center of a fracture 202 opened using a fracturing technique. Fracture 202 in FIG. 2 a may follow the vertical axis in FIG. 2 a , such that fracture 202 may be substantially transverse to horizontal wellbore 204 .
- fracture 202 may extend approximately 500 feet from the center of wellbore 204 .
- FIG. 2 a illustrates that the direction of stress on the rock formation in the area surrounding the fracture is substantially orthogonal to the fracture and also orthogonal to the in-situ direction of maximum horizontal stress of the rock formation.
- the reversal of stress orientation in the area surrounding the fracture may be due to the pressures exerted on the formation from fracturing. This area where the stress orientation is reversed may be referred to as a stress reversal region as shown by stress reversal region 206 of FIG. 2 .
- the degree of reorientation of the stress with respect to the in-situ direction of the maximum horizontal stress may be expressed as an angle from the in-situ direction of stress.
- FIG. 2 b illustrates the stress reorientation with respect to the maximum horizontal stress as caused by creating fracture 202 of FIG. 2 a .
- the orientation of the stress may be substantially orthogonal to the orientation of the in-situ stress such that stress reversal region 206 may be referred to as having a 90° stress reorientation.
- the point where the stress reversal region ends may be referred to as an isotropic point, which may be seen at approximately 140 feet from the center of the fracture along the horizontal axis as shown in FIGS. 2 a and 2 b .
- the distance between fracture 202 and the isotropic point is depicted as s 90° in FIG. 2 b.
- the stresses may still be at an orientation that is not parallel with the maximum horizontal stress of the rock formation.
- the direction of maximum horizontal stress in FIG. 2 a just outside of stress reversal region 206 progressively moves from being parallel with the horizontal axis to being parallel with the vertical axis.
- the orientation angle of the stress adjacent to fracture 204 progressively moves from 90° to 0°. For example, as shown in FIG.
- the rock formation may have a stress reorientation of 10° as shown by s 10° .
- the stress reorientation of rock formation 200 may be 5° at approximately 450 feet from fracture 202 in the positive x-direction and at approximately 400 feet from wellbore 204 in the positive y-direction as shown by s 5° .
- the reorientation of stresses may in turn affect the direction of propagation of subsequent fractures.
- performing fracturing within the stress reversal region of fracture 202 of FIGS. 2 a and 2 b may result in the subsequent fracture propagating parallel to wellbore 204 (longitudinal fracture). In such instances, this phenomenon, often referred to as stress shadowing, may negatively impact the efficiency of the frac stage.
- performing fracturing around s 10° . may cause the subsequent fracture to deviate from a path substantially orthogonal to wellbore 204 by approximately 10°. Therefore, as described in further detail below, by mapping the angle of stress reorientation and the horizontal stress in multiple fractured horizontal wells, the trajectory of each fracture may be estimated. By mapping the trajectory of each induced fracture, the induced fracture spacing may be determined such that it may be minimized without compromising the efficiency of each frac stage.
- FIG. 3 illustrates the geometry of a single transverse fracture 302 of a shale formation 300 that includes a pay zone 304 that may include hydrocarbons (e.g., natural gas) and bounding layers 306 a and 306 b that may bound pay zone 304 , according to some embodiments of the present disclosure.
- Fracture 300 may be modeled based on a variety of properties that may be expressed mathematically. The modeling may be performed by various computer programs, models or combination thereof, configured to simulate and design fracturing operations. The programs and models may include instructions stored on a computer readable medium that are operable to perform, when executed, one or more of the steps described below.
- the computer readable media may include any system, apparatus or device configured to store and retrieve programs or instructions such as a hard disk drive, a compact disc, flash memory or any other suitable device.
- the programs and models may be configured to direct a processor or other suitable unit to retrieve and execute the instructions from the computer readable media.
- the following nomenclature may be used for modeling fracture 302 to describe various properties of fracture 302 :
- the bounding layers 306 a and 306 b may have mechanical properties (E b , v b ) differing from the mechanical properties of pay zone 304 (E p , v p ).
- fracture 302 may be modeled using a numerical model and may have a length in the direction of the x-axis that is equal to 2L f , may have a height in the z-direction that is equal to 2h f and may also have width in the y-direction that is not shown.
- the mechanical behavior of the continuous three-dimensional medium of shale formation 300 may be described mathematically by the equations of equilibrium Eq. (1), the definition of strain Eq. (2) and the constitutive equations Eq. (3).
- the algebraic system of 15 equations for 15 unknowns (6 components of stress ⁇ and strain ⁇ , plus the 3 components of the velocity vector v) may be solved at each node using an explicit, finite difference numerical scheme.
- the Einstein summation convention may apply to indices i, j and k, which take the values 1, 2, 3:
- Pay zone 304 may be homogeneous, isotropic, and purely elastic. Hooke's law relates the components of the strain and stress tensors (constitutive equation):
- poroelastic effects on the stress reorientation around a producing transverse fracture may be ignored in the present example because of the very low permeability of shale and the small amount of fluid leak-off during fracturing.
- the poroelastic effects may be determined and included.
- Shale formation 300 and fracture 302 may also be modeled using a variety of boundary conditions. Displacement along the faces of fracture 302 may be allowed where a constant stress, equal to the net pressure, p net , plus the minimum in-situ horizontal stress ⁇ hmin , is imposed on the faces of fracture 302 to create fracture 302 , or in the present example, to simulate and model the creation of fracture 302 . Therefore, the size (e.g., width, length, height) of fracture 302 may partially be a function of p net .
- fracture 302 may close down on proppant (e.g., sand), which keeps fracture 302 open.
- proppant e.g., sand
- the width of the propped-open fracture may depend on the fractured length and the amount of proppant pumped during the fracturing process.
- the uniform stress boundary condition applied on the fracture face is approximately equal to the pressure required for the proppant to maintain an opening of maximum width w 0 . This pressure value may be smaller than the pressure required to propagate a hydraulic fracture in the same rock.
- the far-field boundaries may be placed at a distance from the fracture equal to at least three times the fracture half-length L f .
- a constant stress boundary condition normal to the “block” faces is applied at outside boundaries. In-situ stresses are initialized prior to the opening of the fracture:
- FIG. 4 illustrates a three dimensional model of multiple transverse fractures in a layered rock formation (e.g., shale formation 300 with pay zone 304 and boundary zones 306 of FIG. 3 ).
- the net pressure required to achieve a specified fracture width may increase with each additional fracture.
- An iterative process may be programmed in order to determine for each fracture, the net pressure corresponding to a given maximum fracture width w 0 .
- the evolution of the net closure stress in the sequential fracturing of a horizontal well is described further below.
- Traditional fracture modeling methods may model fractures perfectly orthogonal to the horizontal wells. However, in order to better quantify the evolution of the direction of propagation of consecutive transverse fractures, it may be advantageous to model subsequent fractures as deviating from the orthogonal path due to the stress reorientation that may be caused by previous fractures.
- Model simplifications may be made in order to tackle this problem.
- multiple inclined fractures are challenging to model on a single numerical mesh.
- the geometry of all fractures may be set from the beginning, which may be very difficult, as the angle of propagation of the subsequent transverse fracture may depend on the mechanical stress perturbation generated by the previous fractures. This may require a complex and time consuming re-meshing after every single fracture stage.
- the net closure stress and the propagation direction may be calculated based on the mechanical stress interference of only the previous fracture.
- FIG. 5 illustrates the calculated propagation of a subsequent fracture (n+1) based on the mechanical stress interference of a previous fracture (n).
- the stress created by the previous propped fracture e.g., fracture (n)
- the net closure stress in the subsequent fracture is equal to the net closure stress of a single transverse fracture (without stress shadow) plus the stresses generated by the previous fracture as shown by Eq. (5) below:
- the trajectory of the subsequent fracture may be approximated by assuming that it will follow the direction of maximum horizontal stress. This may be done by determining the direction of the maximum horizontal stress at one point and having the fracture propagate in that direction. Then, the direction of the maximum horizontal stress may be calculated at another point along the trajectory of the propagation from the previous point and so on to approximate a trajectory for the subsequent fracture as shown for fracture (n+1) of FIG. 5 .
- a determined average angle of deviation may be seen in fracture (n) of FIG. 5 .
- the average angle of deviation may be calculated for fracture (n+1) (e.g., ⁇ f (s f )) from the coordinates of the final position of fracture (n+1). It may be used to model fracture (n+1) in order to calculate the net pressure and trajectory of the subsequent fracture (n+2).
- the propagation direction of subsequent fractures may be a function of the location of the subsequent fracture with respect to areas of the rock formation that have experienced stress reorientation caused by propagating previous fractures. Accordingly, the spacing between a previous fracture and a subsequent fracture may influence the propagation direction of the subsequent fracture.
- FIGS. 6 , 7 , 9 , 11 and 12 described further below illustrate examples of the trajectory of multiple fractures according to various spacing distances between the fractures. The fractures in FIGS. 6 , 7 , 9 , 11 and 12 may be determined using the process described above with respect to FIGS. 3 through 5 and may be done by any suitable computer program. In FIGS. 6 , 7 , 9 , 11 and 12 , the fractures depicted may be induced in separate consecutive stages.
- fracture 1 of FIGS. 6 , 7 , 9 , 11 and 12 may be induced first, fracture 2 may be induced second, etc.
- results of the different simulations with respect to different spacing distances between fractures may be used to determine optimum fracture spacing for a particular well and/or formation.
- TABLE 1 below illustrates the parameters of the rock formation used in the examples of FIGS. 6-23 as taken from a shale gas well in the Barnett shale in Texas.
- FIG. 6 illustrates the results of calculating fracture propagation with each fracture spaced approximately 400 feet apart.
- transverse fractures may propagate away from the previous fracture with a small angle of deviation from the orthogonal path (less than 2°), as expected from the angle of stress reorientation profile shown in FIG. 1 (simulated using the same parameters).
- the spacing is reduced to 300 ft. (as shown in FIG. 7 )
- a closer look at the fracture trajectory shows that after fracture 5 , the fracture may initially propagate toward the previous fracture and then at some distance from the wellbore, the fracture may start propagating away from the previous fracture. Plotting the stress distribution around an oblique fracture reveals the explanation behind this non-trivial trend as shown in FIG. 8 .
- FIG. 8 illustrates the stress distribution of a rock formation 800 with the stress distribution being influenced by the propagation of fracture 4 of FIG. 7 .
- a stress reversal zone 804 a zone where the subsequent fracture (e.g., a fracture 5 of FIG. 7 ) may be attracted by the previous fracture (e.g., an attraction zone 806 ), and another zone where the subsequent fracture may propagate away from the previous fracture (e.g., a repulsion zone 808 ).
- subsequent transverse fractures may propagate in both zones as shown by fracture 5 .
- the size of attraction zone 806 may be function of the net pressure, the in-situ stress contrast and the average angle of deviation from the orthogonal path of fracture 4 .
- the initiation point of fracture 5 may be located within attraction zone 806 caused by the propagation of fracture 4 , thus fracture 5 may initially propagate back toward fracture 4 until it leaves attraction zone 806 .
- FIG. 9 illustrates fracture propagation with the fracture spacing reduced to 250 ft. in accordance with the present example. Because of the closer spacing, the amount of fracture deviation is larger. For instance, fractures 2 , 5 and 8 may propagate away from the previous fracture at an angle ⁇ f >5°. But what mostly stands out is the fact that under a critical value of the fracture spacing, the attraction zone associated with fractures 2 , 5 and 8 may cause fractures 3 , 6 , and 9 , respectively, to intersect fractures 2 , 5 and 8 respectively. The practical consequence of such intersections may be much less efficient drainage of the reservoir, even if the fractures are initiated closer to each other.
- a two-fracture system may be simulated to calculate the stress distribution of the rock formation.
- the stress distribution that may affect the fracture 4 may be modeled based on both the fractures 2 and 3 .
- the stress distribution around the fracture system with respect to the fractures 2 and 3 is shown in FIG. 10 .
- FIGS. 11 and 12 also illustrate fracture propagation as calculated for a 200-ft. and a 150-ft. spacing respectively of the present example simulation.
- the “unsuccessful” fractures e.g., fractures 3 , 5 , 7 and 10 of FIG. 11 and fractures 2 , 4 , 6 , 8 and 10 of FIG. 12
- the unsuccessful fractures may be caused by initiating the fractures within the stress reversal region of the previous fracture, which is located inside the attraction zone associated with the previous fracture (e.g., attraction zone 806 of FIG. 8 ).
- the present example 150-ft. spacing
- only every other fracturing stage effectively stimulates the shale, thus possibly leaving significant portions of the reservoir inadequately drained.
- FIGS. 13 a and 13 b illustrate the impact of fracture spacing on the angle of deviation of the fractures from the orthogonal path.
- the efficiency of fracturing stages may be negatively affected as shown by the large variations in deviation angles with respect to spacings 250, 200 and 150 feet apart. Accordingly, the gain in reservoir drainage at these spacings may be marginal compared to the additional cost represented by an increased number of fracture stages. This result suggests that because of mechanical stress interference, spacing transverse fractures ever closer to each other may not be a desirable completion strategy.
- FIG. 14 illustrates the impact of fracture spacing on the evolution of the net closure stress.
- the net closure stress only increases with each new stage until reaching a plateau.
- the net pressure may have an up and down trend.
- Counting the number of times the net fracturing pressure decreases from one stage to another may indicate the number of unsuccessful fracture stages identified in FIGS. 10 , 12 , 13 a and 13 b .
- the decrease in the fracture closure stress may be a consequence of the smaller mechanical stress interference (stress shadow) generated by the previous fracture when propagating into stimulated regions of the reservoir instead of orthogonal to the well.
- the effective spacing may only equal to 300 ft., as every other fracture may be longitudinal with respect to the wellbore. Accordingly, doubling the number of stages for 150 ft. spacing compared to the 300 ft. spacing may grant very little improvement in well production.
- modeling deviation from the orthogonal path for fractures may reveal a new up-and-down trend in the evolution of the net closure stress.
- This up and down trend may indicate that the spacing between fractures may be too close to generate any improvement in well production. Therefore, the net closure stress at various spacings may be analyzed to determine the closest spacing that may not yield an up and down net pressure such that optimal spacing of fractures may be determined.
- the propagation direction of each fracture may be estimated instead of assuming that the propagation direction is orthogonal to the well as is traditionally done.
- Fracture spacing may also be determined by analyzing the stress reversal region associated with a previous fracture and by initiating the subsequent fracture outside of the stress reversal region.
- FIGS. 15 a and 15 b illustrate the horizontal stress of a rock formation 1500 with a stress reversal region 1502 associated with a fracture (n+1). Fracture (n+1) may run along the vertical axis of FIGS. 15 a and 15 b and a wellbore 1504 may run along the horizontal axis.
- stress reversal region 1502 may extend approximately 230 ft. from fracture 1502 along wellbore 1504 as shown by isotropic point s 90° .
- a subsequent fracture (e.g., a fracture (n+2)) may not be initiated closer than 230 ft. from fracture (n+1) because, as mentioned above, the subsequent fracture may propagate parallel with wellbore 1504 and may not increase hydrocarbon production from wellbore 1504 .
- FIG. 15 b illustrates that in the present example at point s 10° , (e.g., approximately 430 ft. from fracture (n+1)) the stress reorientation of rock formation 1500 may be 10° and at point s 5° (e.g., 600 ft. from fracture (n+1)) the stress reorientation of rock formation 1500 may be 5°.
- the stress reorientations of 10° and 5° may be such that a subsequent fracture (e.g., a fracture (n+2)) initiated between points s 10° and s 5° may not intersect fracture (n+1) although the subsequent fracture may deviate somewhat from an orthogonal path due to the stress reorientation.
- FIGS. 16 a and 16 b illustrate the differences between performing consecutive fracturing and alternate fracturing.
- each fracture starting with fracture “ 1 ” may be initiated one after another in a consecutive order.
- FIG. 16 b it can be seen that two fractures may be initiated consecutively (e.g., fractures “ 1 ” and “ 2 ” of FIG. 16 b ), however the two previous fractures may be sufficiently far apart that a third fracture (e.g., fracture “ 3 ” of FIG. 16 b ) may be initiated between the two previous fractures, such that the fractures alternate.
- a third fracture e.g., fracture “ 3 ” of FIG. 16 b
- FIGS. 17 a and 17 b illustrate the stress orientation of a rock formation 1700 with stress reversal regions 1701 and 1702 associated with a fracture “ 1 ” and fracture “ 2 ” respectively.
- fractures “ 1 ” and “ 2 ” may be placed approximately 650 ft. from each other.
- FIG. 17 b illustrates that the distance between stress reversal regions 1701 and 1702 may be approximately 20 ft. in the present example. Therefore, by initiating fracture “ 3 ” in the middle of fractures “ 1 ” and “ 2 ,” both stress reversal regions 1701 and 1702 may be avoided by a narrow margin of 20 ft.
- FIGS. 17 a and 17 b illustrate that the spacing of fractures “ 1 ” and “ 2 ” may be determined such that stress reversal regions 1701 and 1702 may not intersect, but also such that they are sufficiently far apart to allow for the initiation of a third fracture between them. Accordingly, by analyzing the size of the stress reversal regions of the two “end” fractures (e.g., fractures “ 1 ” and “ 2 ”) in alternate fracturing, the spacing between the two may be more efficiently determined for placement of the “middle” fracture (e.g., fracture “ 3 ”).
- the spacing between the two may be more efficiently determined for placement of the “middle” fracture (e.g., fracture “ 3 ”).
- Hydraulic fracture interaction with pre-existing natural fractures may be a function of a term called the relative net pressure R n . This parameter may be inversely proportional to the local deviatoric stress in which the fracture propagates as shown below in Equation (6).
- High values of the relative net pressure R n may favor fracture path complexity.
- a hydraulic fracture propagating in a region of low stress contrast may create larger networks of interconnected fractures.
- the propensity of the alternate fracturing sequence to generate fracture complexity may be quantified and compared to the more conventional fracturing approach.
- the average value of the stress contrast seen by a propagating middle fracture in the alternate fracturing sequence may be measured for different values of the spacing between the outside fractures ( 2 s f ).
- FIG. 18 shows that the spacing for which the stress contrast is lowest may be equal to the minimum fracture spacing previously calculated (325 ft).
- the minimum fracture spacing in the alternate fracturing sequence may also be the optimum case for creating fracture complexity.
- FIG. 19 illustrates that the deviatoric stress may approach zero in the near wellbore region in the case of the optimum spacing in the alternate fracturing sequence (325 ft).
- the stress contrast may remain lower than 10 psi, which may be equal to 10% of the in-situ stress contrast. It is only in the second half of the fracture propagation that the local stress contrast increases significantly.
- choosing the alternate fracturing sequence may result in high fracture complexity in the near-wellbore region as a result of the propagation of the “middle fracture”.
- FIG. 20 illustrates an example of fracture spacing that may be done with multiple horizontal lateral wells.
- FIG. 20 illustrates three horizontal wells (HW 1 , HW 2 and HW 3 ) that may run substantially parallel to each other through a hydrocarbon reservoir in a horizontal plane that may be substantially parallel with the ground.
- the wells of FIG. 20 may be described by variables representing fracture dimensions (L f , h f ), fracture spacing (s f ) and the inter-well spacing (s w ).
- the middle well (HW 2 ) may be used to propagate a fracture (e.g., fracture “ 3 ”) in between two pairs of fractures previously initiated from the outside wells (e.g., fractures “ 1 ” and “ 2 ” of well HW 1 and fractures “ 1 ′” and “ 2 ′” of well HW 3 ).
- fractures “ 3 ” e.g., fracture “ 3 ”
- Such strategy may allow for benefiting from the propagation of a “middle fracture”, like in alternate fracturing completions, without the need for special downhole tools. Indeed, in each lateral well (e.g., HW 1 , HW 2 and HW 3 ) the fractures may be initiated in a conventional consecutive sequence.
- the spacing between fractures in such multi-lateral sequences may be determined by analyzing the stress distribution (e.g., stress reversal regions) associated with the fractures. For example, the stress distribution between two pairs of fractures (e.g., fractures “ 1 ” and “ 2 ” of well HW 1 and fractures “ 1 ′” and “ 2 ′” of well HW 3 of FIG. 20 ) propagated from the outside laterals HW 1 and HW 3 is shown in FIGS. 21 a and 21 b .
- stress distribution e.g., stress reversal regions
- a fracture spacing s f associated with the fractures of each well e.g., spacing between fractures “ 1 ”, “ 2 ”, “ 4 ”, “ 6 ”, “ 8 ”, etc., of well HW 1 , fractures “ 3 ”, “ 5 ”, “ 7 ”, “ 9 ”, etc., of well HW and fractures “ 1 ′”, “ 2 ′”, “ 4 ′”, “ 6 ′”, “ 8 ′”, etc., of well HW 3
- a fracture spacing s f associated with the fractures of each well e.g., spacing between fractures “ 1 ”, “ 2 ”, “ 4 ”, “ 6 ”, “ 8 ”, etc., of well HW 1 , fractures “ 3 ”, “ 5 ”, “ 7 ”, “ 9 ”, etc., of well HW and fractures “ 1 ′”, “ 2 ′”, “ 4 ′”, “ 6 ′”, “ 8 ′”, etc., of well HW 3
- 600 ft 600 ft.
- a well spacing (s w ) between wells HW 1 , HW 2 and HW 3 may be approximately equal to 500 ft.
- the above spacing may be determined by analyzing the direction of maximum horizontal stress associated with the fractures. For example, the direction of maximum horizontal stress may be reversed everywhere along the outside laterals as shown in FIGS. 21 a and 21 b . Thus, the outside fractures “ 1 ” and “ 2 ” are too closely spaced to allow propagation of a transverse fracture from the outside laterals HW 1 and HW 3 , similarly to the alternate fracturing sequence.
- the direction of maximum horizontal stress may still allow propagation of a transverse fracture.
- the distance of transverse propagation, L transverse , of fracture “ 3 ” of HW 2 may be at a maximum at mid-distance from the previous fractures and may be function of not only the spacing between the outside fractures but also the inter-well spacing (s w ).
- the zone of transverse fracture propagation can also be identified when plotting the angle of stress reorientation as shown in FIG. 21 b.
- FIG. 23 illustrates the local stress contrast that may be recorded along the assumed propagation direction of a middle fracture (e.g., fracture “ 3 ” of FIG. 20 ).
- it may be advantageous to position the horizontal laterals close to each other, but not closer than a distance equal to the fracture half-height. Otherwise, the benefit of propagating long transverse fractures may be lost. Indeed, that may result in fracturing zones of the reservoir that are already stimulated.
- the distance of transverse fracture propagation may be sensitive to the fracture spacing when the inter-well spacing is small.
- Transverse propagation length may be decreased by over 50% as the fracture spacing decreases from 650 ft. to 600 ft., which is only a 50 ft. spacing differential. Therefore, in the present example and similarly to the case of the alternate fracturing sequence in a single well, the spacing between the outside fractures may be at least be equal to 650 ft.
- the middle fracture may be initiated from the middle well (and not from the outside well), where the refracturing interval is wide enough to allow fracture initiation from multiple perforation clusters.
- the spacing of the fractures may be determined to improve production from wells, while also improving the efficiency of each fracturing operation.
- Such stress reorientation analysis may be used for consecutive fracturing, for alternate fracturing and/or for multiple horizontal fracturing operations.
- the in-situ stress contrast which is the difference between the maximum horizontal stress and the minimum horizontal stress, may influence the stress interference created by multiple consecutive fractures, including fracture intersection.
- the evolution of the fracturing pressures during multi-stage fracturing of horizontal wells may be impacted by the in-situ stress contrast, just like it is impacted by the fracture spacing (e.g., as shown in FIG. 14 ).
- the maximum horizontal stress may be more difficult to evaluate in the field. Knowing the value of the maximum horizontal stress may prove useful in modeling multiple engineering problems in the oil and gas industry, including hydraulic fracturing and wellbore stability and sand production issues.
- the proposed method may be used to calculate the evolution of the net closure stress in a given well for different values of the maximum horizontal stress. By comparing the calculated pressure profiles to the field-measured fracturing pressures, the value of the maximum horizontal stress may be determined for the well in question.
- the stress redistribution of a rock formation caused by propagating a fracture may also be a function of the induced fracture length, fracture width, fluid rheology and the injection rates associated with propagating the fracture.
- the propagation of subsequent fractures may be a function of the stress redistribution caused by previous fractures. Therefore, the analysis described above may also be used to determine one or more of the above mentioned properties to better improve fracturing efficiency. For example, in some instances for a particular fracture size, the determined optimal spacing may be too far apart. Accordingly, the spacing may be set at a fixed value and another factor that may affect stress reorientation (e.g., fracture width) may be modified. The stress reorientation, and the propagation and net closure stress of consecutive fractures may be calculated for different values of the fracture width such that an optimum width of the fractures may be determined.
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
Description
- This application claims benefit under 35 U.S.C. §119(e) of U.S. Provisional Application Ser. No. 61/501,003, entitled “METHOD FOR DETERMINING SPACING OF HYDRAULIC FRACTURES IN A ROCK FORMATION,” filed Jun. 24, 2011, the entire content of which is incorporated herein by reference.
- This invention was made with government support under DE-AC26-07NT42677 awarded by The Department of Energy. The United States Government has certain rights in the invention.
- The present disclosure relates in general to well drilling and, more particularly, to a method for determining spacing of hydraulic fractures in a rock formation.
- Hydrocarbon (e.g., oil, natural gas, etc.) reservoirs may be found in geologic formations that have little to no porosity (e.g., shale, sandstone etc.). The hydrocarbons may be trapped within fractures and pore spaces of the formation. Additionally, the hydrocarbons may be adsorbed onto organic material of the shale formation. The rapid development of extracting hydrocarbons from these unconventional reservoirs can be tied to the combination of horizontal drilling and hydraulic fracturing (“fracing”) of the formations. Horizontal drilling has allowed for drilling along and within hydrocarbon reservoirs of a formation to better capture the hydrocarbons trapped within the reservoirs. Additionally, more hydrocarbons may be captured by increasing the number of fractures in the formation and/or increasing the size of already present fractures through fracing. The spacing between fractures as well as the ability to stimulate the fractures naturally present in the rock may be major factors in the success of horizontal completions in unconventional hydrocarbon reservoirs.
- In one embodiment, a method is disclosed comprising determining an expected trajectory of induced fractures, analyzing net pressure associated with the induced fractures, and determining at least one of spacing of induced fractures and a property of the induced fractures based on the net pressure. Computer-readable medium containing the same are also disclosed.
- In an alternative embodiment, a method of optimizing fracture spacing is disclosed. The method includes propagating an initial fracture, measuring pressure associated with propagating the initial fracture, determining a minimum spacing required to prevent a second fracture from intersecting the initial fracture, and propagating the second fracture at least the minimum spacing distance away from the initial fracture.
- In other embodiments, also disclosed is a method of optimizing fracture spacing. The method includes analyzing stresses associated with a first set of at least one fractures associated with a first well, analyzing stresses associated with a second set of at least one fractures associated with a second well, and determining spacing of a fracture associated with a third well such that the fracture associated with the third well does not intersect with the first set and the second set of fractures, the third well running between the first and second wells.
- In still other embodiments, also disclosed is a method of determining maximum horizontal pressure. The method includes measuring an actual pressure during each stage of fracturing of a rock formation, determining a theoretical expected pressure during each stage of fracturing of the rock formation, and determining a maximum horizontal pressure of the rock formation based at least on a comparison of the theoretical expected pressure and the measured actual pressure.
-
FIG. 1 illustrates an example schematic of a gas well configured to extract natural gas from a gas rich shale formation, according to some embodiments of the present disclosure; -
FIGS. 2 a and 2 b illustrate examples of the reorientation of stresses in a rock formation due to the placement of a fracture orthogonal to a horizontal well according to some embodiments of the present disclosure; -
FIG. 3 illustrates the geometry of a single transverse fracture of a shale formation that includes a pay zone that may include hydrocarbons (e.g., natural gas) and bounding layers that may bound the pay zone, according to some embodiments of the present disclosure; -
FIG. 4 illustrates a three dimensional model of multiple transverse fractures in a layered rock formation according to some embodiments of the present disclosure; -
FIG. 5 illustrates the calculated propagation of a subsequent fracture (n+1) based on the mechanical stress interference of a previous fracture (n), according to some embodiments of the present disclosure; -
FIG. 6 illustrates the results of calculating fracture propagation with each fracture spaced approximately 400 feet apart, according to some embodiments of the present disclosure; -
FIG. 7 illustrates the results of calculating fracture propagation with each fracture spaced approximately 300 feet apart, according to some embodiments of the present disclosure; -
FIG. 8 illustrates the stress distribution of a rock formation with the stress distribution being influenced by the propagation of a fracture produced during a fourth stage of a fracture treatment, according to some embodiments of the present disclosure; -
FIG. 9 illustrates fracture propagation with the fracture spacing reduced to 250 ft., according to some embodiments of the present disclosure; -
FIG. 10 illustrates the stress distribution of a rock formation caused by fractures ofFIG. 9 , according to some embodiments of the present disclosure; -
FIG. 11 illustrates fracture propagation with the fracture spacing reduced to 200 ft., according to some embodiments of the present disclosure; -
FIG. 12 illustrates fracture propagation with the fracture spacing reduced to 150 ft., according to some embodiments of the present disclosure; -
FIGS. 13 a and 13 b illustrate the impact of fracture spacing on the angle of deviation of the fractures from the orthogonal path, according to some embodiments of the present disclosure; -
FIG. 14 illustrates the impact of fracture spacing on the evolution of the net closure stress, according to some embodiments of the present disclosure; -
FIGS. 15 a and 15 b illustrate the horizontal stress of a rock formation with a stress reversal region associated with a fracture, according to some embodiments of the present disclosure; -
FIGS. 16 a and 16 b illustrate the differences between performing consecutive fracturing and alternate fracturing, according to some embodiments of the present disclosure; -
FIGS. 17 a and 17 b illustrate the stress orientation of a rock formation with stress reversal regions associated with a first fracture and a second fracture, according to some embodiments of the present disclosure; -
FIG. 18 shows a spacing of fractures for which the stress contrast may be lowest, according to some embodiments of the present disclosure; -
FIG. 19 illustrates that the deviatoric stress may approach zero in a near wellbore region in the case of optimum spacing in an alternate fracturing sequence, according to some embodiments of the present disclosure; -
FIG. 20 illustrates an example of fracture spacing that may be done with multiple horizontal lateral wells, according to some embodiments of the present disclosure; -
FIGS. 21 a and 21 b illustrate the stress distribution between two pairs of fractures propagated from outside lateral wells, according to some embodiments of the present disclosure; -
FIG. 22 illustrates the relationship between the length of a “middle fracture” propagating from a center lateral well for different values of the fracture length with respect to the inter-well spacing, according to some embodiments of the present disclosure; and -
FIG. 23 illustrates the local stress contrast that may be recorded along the assumed propagation direction of a middle fracture, according to some embodiments of the present disclosure. -
FIG. 1 illustrates an example schematic of a gas well 100 configured to extract natural gas from a gasrich shale formation 102. Well 100 may be drilled using horizontal drilling methods to create awellbore 104 that runs within and alongshale formation 102. In the present example,wellbore 104 may be drilled such thatwellbore 104 runs perpendicular to the maximum horizontal in-situ stresses ofshale formation 102 to obtain better production of well 100. -
Shale formation 102 may produce natural gas that is trapped in fractures and pore spaces ofshale formation 102. The natural gas may also be adsorbed in organic material included in the shale ofshale formation 102. As well bore 104 runs throughshale formation 102, well bore 104 may also run through fractures (not expressly shown) ofshale formation 102. The gas in the fractures may enter well bore 104 and may accordingly be retrieved at adrilling rig 106 of well 100. As gas leaves the fractures ofshale formation 102, the gas adsorbed on the organic material may be released into the fractures such that the adsorbed gas may also be retrieved. As the number of fractures ofshale formation 102 that well bore 104 passes through increases, the amount of gas that may be produced by well 100 may also increase. Therefore, increasing the number of fractures inshale formation 102 along well bore 104 may increase the gas production of well 100. - The number and/or size of fractures in
shale formation 102 may be increased using hydraulic fracturing (“fracing”). Fracing may refer to any process used to initiate and propagate a fracture in a rock formation. Additionally fracing may be used to increase existing fractures in a rock formation. Fracing may include forcing a hydraulic fluid in a fracture of a rock formation to increase the size of the fracture and introducing proppant (e.g., sand) in the newly induced fracture to keep the fracture open. The fracture may be an existing fracture in the formation, or may be initiated using a variety of techniques known in the art. The amount of pressure needed to extend and propagate the fracture may be referred to as the “fracturing pressure.” - As shown in further detail with respect to
FIGS. 2 a and 2 b, producing fractures during fracing may change the stress properties of a rock formation. Accordingly, subsequent transverse fractures initiated from a horizontal well may deviate toward or away from the previous fracture depending on the stress reorientation caused by the fracing. The stress reorientation may be a function of mechanical properties of the reservoir rock, fracture spacing, and the orientation of the previous fracture. As described in further detail below, in some instances spacing frac treatments too close together may result in a subsequent fracture intersecting with a previous fracture. Therefore, in such instances, the contribution of the subsequent fracture in hydrocarbon production may be reduced or minimized. - As disclosed in further detail below, the spacing of performing fracturing operations for a hydrocarbon well, such as gas well 100, may be determined using net pressure measurements to determine a minimum frac spacing that also reduces the likelihood of subsequent fractures intersecting and interfering with previous fractures. In some instances, net pressure may be determined by surface pressures or down-hole pressures during fracturing. Therefore, the selection of fracture spacing may be such that the number of fractures initiated from a horizontal wellbore may be increased while also reducing the likelihood that the fractures may interfere and/or intersect with each other to allow for better production rates and depletion of hydrocarbons with each induced fracture. Therefore, the economic efficiency of fracturing may be improved and the cost of retrieving hydrocarbons from tight rock formations (e.g., shale formation 102) may be reduced.
- Modifications, additions or omissions may be made to
FIG. 1 without departing from the present disclosure. For example, althoughFIG. 1 is described as performing fracing with respect to a shale gas formation, the present disclosure may be used to improve frac spacing for any suitable formation (e.g., a tight sand formation, coalbed methane, sandstone, limestone, oil shale). Additionally, although well 100 is described as being used to extract natural gas, it is understood that the principles described herein may be used to extract any other suitable hydrocarbon. -
FIGS. 2 a and 2 b illustrate examples of the reorientation of stresses in a rock formation due to the placement of a fracture orthogonal (or transverse) to a horizontal well according to some embodiments of the present disclosure. The opening of a propped transverse fracture in horizontal wells through hydraulic fracturing may cause a reorientation of stresses in the rock formation surrounding the fracture. -
FIG. 2 a may represent the stress of the rock formation in a horizontal plane (e.g., a plane substantially parallel to the ground. Accordingly, the vertical axis ofFIG. 2 a may represent the distance in the x direction from the center of a substantiallyhorizontal wellbore 204. The horizontal axis ofFIG. 2 a may represent the distance in the y direction from the center of afracture 202 opened using a fracturing technique. Fracture 202 inFIG. 2 a may follow the vertical axis inFIG. 2 a, such thatfracture 202 may be substantially transverse tohorizontal wellbore 204. In the present example,fracture 202 may extend approximately 500 feet from the center ofwellbore 204.FIG. 2 a illustrates that the direction of stress on the rock formation in the area surrounding the fracture is substantially orthogonal to the fracture and also orthogonal to the in-situ direction of maximum horizontal stress of the rock formation. The reversal of stress orientation in the area surrounding the fracture may be due to the pressures exerted on the formation from fracturing. This area where the stress orientation is reversed may be referred to as a stress reversal region as shown by stress reversal region 206 ofFIG. 2 . - The degree of reorientation of the stress with respect to the in-situ direction of the maximum horizontal stress may be expressed as an angle from the in-situ direction of stress.
FIG. 2 b illustrates the stress reorientation with respect to the maximum horizontal stress as caused by creatingfracture 202 ofFIG. 2 a. For example, in stress reversal region 206 ofFIG. 2 a, the orientation of the stress may be substantially orthogonal to the orientation of the in-situ stress such that stress reversal region 206 may be referred to as having a 90° stress reorientation. The point where the stress reversal region ends may be referred to as an isotropic point, which may be seen at approximately 140 feet from the center of the fracture along the horizontal axis as shown inFIGS. 2 a and 2 b. The distance betweenfracture 202 and the isotropic point is depicted as s90° inFIG. 2 b. - Outside of the stress reversal region, the stresses may still be at an orientation that is not parallel with the maximum horizontal stress of the rock formation. For example, the direction of maximum horizontal stress in
FIG. 2 a just outside of stress reversal region 206 progressively moves from being parallel with the horizontal axis to being parallel with the vertical axis. Additionally, it can be seen inFIG. 2 b that the orientation angle of the stress adjacent to fracture 204 progressively moves from 90° to 0°. For example, as shown inFIG. 2 b, in the present example, at approximately 320 feet fromfracture 202 in the positive y-direction and at approximately 300 feet fromwellbore 204 in the positive x-direction the rock formation may have a stress reorientation of 10° as shown by s10°. Similarly, in the present example, the stress reorientation ofrock formation 200 may be 5° at approximately 450 feet fromfracture 202 in the positive x-direction and at approximately 400 feet fromwellbore 204 in the positive y-direction as shown by s5°. - The reorientation of stresses may in turn affect the direction of propagation of subsequent fractures. For example, performing fracturing within the stress reversal region of
fracture 202 ofFIGS. 2 a and 2 b may result in the subsequent fracture propagating parallel to wellbore 204 (longitudinal fracture). In such instances, this phenomenon, often referred to as stress shadowing, may negatively impact the efficiency of the frac stage. As an additional example, performing fracturing around s10°. may cause the subsequent fracture to deviate from a path substantially orthogonal to wellbore 204 by approximately 10°. Therefore, as described in further detail below, by mapping the angle of stress reorientation and the horizontal stress in multiple fractured horizontal wells, the trajectory of each fracture may be estimated. By mapping the trajectory of each induced fracture, the induced fracture spacing may be determined such that it may be minimized without compromising the efficiency of each frac stage. -
FIG. 3 illustrates the geometry of a singletransverse fracture 302 of ashale formation 300 that includes apay zone 304 that may include hydrocarbons (e.g., natural gas) and boundinglayers zone 304, according to some embodiments of the present disclosure. Fracture 300 may be modeled based on a variety of properties that may be expressed mathematically. The modeling may be performed by various computer programs, models or combination thereof, configured to simulate and design fracturing operations. The programs and models may include instructions stored on a computer readable medium that are operable to perform, when executed, one or more of the steps described below. The computer readable media may include any system, apparatus or device configured to store and retrieve programs or instructions such as a hard disk drive, a compact disc, flash memory or any other suitable device. The programs and models may be configured to direct a processor or other suitable unit to retrieve and execute the instructions from the computer readable media. The following nomenclature may be used for modelingfracture 302 to describe various properties of fracture 302: - Ep=Young's modulus of the pay zone, Pa (psi)
Eb=Young's modulus of the bounding layers, Pa (psi)
vp=Poisson's ratio in the pay zone
vb=Poisson's ratio in the bounding layers
K=dry bulk modulus, Pa (psi)
G=shear modulus, Pa (psi)
Lf=fracture half-length, m (ft)
hf=fracture half-height, m (ft)
hp=pay zone half-thickness, m (ft)
w0=maximum fracture width, m (ft)
σv=vertical in-situ stress, Pa (psi)
σhmax=maximum horizontal in-situ stress, Pa (psi)
σhmin=minimum horizontal in-situ stress, Pa (psi) - The bounding layers 306 a and 306 b may have mechanical properties (Eb, vb) differing from the mechanical properties of pay zone 304 (Ep, vp). In the present example,
fracture 302 may be modeled using a numerical model and may have a length in the direction of the x-axis that is equal to 2Lf, may have a height in the z-direction that is equal to 2hf and may also have width in the y-direction that is not shown. - The mechanical behavior of the continuous three-dimensional medium of
shale formation 300 may be described mathematically by the equations of equilibrium Eq. (1), the definition of strain Eq. (2) and the constitutive equations Eq. (3). The algebraic system of 15 equations for 15 unknowns (6 components of stress σ and strain ε, plus the 3 components of the velocity vector v) may be solved at each node using an explicit, finite difference numerical scheme. The Einstein summation convention may apply to indices i, j and k, which take thevalues -
- Pay
zone 304 may be homogeneous, isotropic, and purely elastic. Hooke's law relates the components of the strain and stress tensors (constitutive equation): -
- Where,
-
- and
-
- The impacts of poroelastic effects on the stress reorientation around a producing transverse fracture may be ignored in the present example because of the very low permeability of shale and the small amount of fluid leak-off during fracturing. However, in other models of other rock formations, the poroelastic effects may be determined and included.
-
Shale formation 300 andfracture 302 may also be modeled using a variety of boundary conditions. Displacement along the faces offracture 302 may be allowed where a constant stress, equal to the net pressure, pnet, plus the minimum in-situ horizontal stress σhmin, is imposed on the faces offracture 302 to createfracture 302, or in the present example, to simulate and model the creation offracture 302. Therefore, the size (e.g., width, length, height) offracture 302 may partially be a function of pnet. - At the end of the fracturing process, fracture 302 may close down on proppant (e.g., sand), which keeps
fracture 302 open. The width of the propped-open fracture may depend on the fractured length and the amount of proppant pumped during the fracturing process. The uniform stress boundary condition applied on the fracture face is approximately equal to the pressure required for the proppant to maintain an opening of maximum width w0. This pressure value may be smaller than the pressure required to propagate a hydraulic fracture in the same rock. To simulate a large enough reservoir volume and avoid boundary effects, the far-field boundaries may be placed at a distance from the fracture equal to at least three times the fracture half-length Lf. A constant stress boundary condition normal to the “block” faces is applied at outside boundaries. In-situ stresses are initialized prior to the opening of the fracture: -
- Following modeling of the first fracture, subsequent fractures may also be modeled. After the first fracture is created, its geometry may be represented as being fixed (e.g., no further displacement is allowed). In the present example, it may be assumed that the compression of the proppant placed inside the previous fractures is negligible as subsequent fractures are opened. Subsequent transverse fractures may be modeled using boundary conditions similar to the first fracture as described above.
FIG. 4 illustrates a three dimensional model of multiple transverse fractures in a layered rock formation (e.g.,shale formation 300 withpay zone 304 and boundary zones 306 ofFIG. 3 ). - The net pressure required to achieve a specified fracture width may increase with each additional fracture. An iterative process may be programmed in order to determine for each fracture, the net pressure corresponding to a given maximum fracture width w0. The evolution of the net closure stress in the sequential fracturing of a horizontal well is described further below.
- Traditional fracture modeling methods may model fractures perfectly orthogonal to the horizontal wells. However, in order to better quantify the evolution of the direction of propagation of consecutive transverse fractures, it may be advantageous to model subsequent fractures as deviating from the orthogonal path due to the stress reorientation that may be caused by previous fractures.
- Model simplifications may be made in order to tackle this problem. As opposed to perfectly orthogonal fractures, multiple inclined fractures are challenging to model on a single numerical mesh. In a finite difference model, the geometry of all fractures may be set from the beginning, which may be very difficult, as the angle of propagation of the subsequent transverse fracture may depend on the mechanical stress perturbation generated by the previous fractures. This may require a complex and time consuming re-meshing after every single fracture stage.
- Accordingly, for a more a simplified approach, in the present example, the net closure stress and the propagation direction may be calculated based on the mechanical stress interference of only the previous fracture.
FIG. 5 illustrates the calculated propagation of a subsequent fracture (n+1) based on the mechanical stress interference of a previous fracture (n). For each subsequent fracture (e.g., fracture (n+1)), the stress created by the previous propped fracture (e.g., fracture (n)) in the direction perpendicular to it is computed at some distance from the fracture. The net closure stress in the subsequent fracture is equal to the net closure stress of a single transverse fracture (without stress shadow) plus the stresses generated by the previous fracture as shown by Eq. (5) below: -
p net n+1 +p net 1+Δσyy n(s f) (5) - Based on the stress distribution around a transverse fracture, the trajectory of the subsequent fracture may be approximated by assuming that it will follow the direction of maximum horizontal stress. This may be done by determining the direction of the maximum horizontal stress at one point and having the fracture propagate in that direction. Then, the direction of the maximum horizontal stress may be calculated at another point along the trajectory of the propagation from the previous point and so on to approximate a trajectory for the subsequent fracture as shown for fracture (n+1) of
FIG. 5 . - A determined average angle of deviation may be seen in fracture (n) of
FIG. 5 . The average angle of deviation may be calculated for fracture (n+1) (e.g., θf(sf)) from the coordinates of the final position of fracture (n+1). It may be used to model fracture (n+1) in order to calculate the net pressure and trajectory of the subsequent fracture (n+2). - As mentioned above, the propagation direction of subsequent fractures may be a function of the location of the subsequent fracture with respect to areas of the rock formation that have experienced stress reorientation caused by propagating previous fractures. Accordingly, the spacing between a previous fracture and a subsequent fracture may influence the propagation direction of the subsequent fracture.
FIGS. 6 , 7, 9, 11 and 12 described further below illustrate examples of the trajectory of multiple fractures according to various spacing distances between the fractures. The fractures inFIGS. 6 , 7, 9, 11 and 12 may be determined using the process described above with respect toFIGS. 3 through 5 and may be done by any suitable computer program. InFIGS. 6 , 7, 9, 11 and 12, the fractures depicted may be induced in separate consecutive stages. For example,fracture 1 ofFIGS. 6 , 7, 9, 11 and 12 may be induced first,fracture 2 may be induced second, etc. As described further below, the results of the different simulations with respect to different spacing distances between fractures may be used to determine optimum fracture spacing for a particular well and/or formation. TABLE 1 below illustrates the parameters of the rock formation used in the examples ofFIGS. 6-23 as taken from a shale gas well in the Barnett shale in Texas. -
TABLE 1 Reservoir parameters for Barnett shale gas well Barnett shale gas Pay zone Young's Modulus Ep (psi) 7.3 × 106 Bounding layer Young's Modulus Eb (psi) 3.0 × 106 Poisson's Ratio νp = νb 0.2 σhmax (psi) 6400 σhmin (psi) 6300 Depth (ft) 7000 Pay zone half-thickness hp (ft) 150 Fracture half-height hf (ft) 150 Fracture half-length Lf (ft) 500 Fracture maximum width w0 (mm) 4 -
FIG. 6 illustrates the results of calculating fracture propagation with each fracture spaced approximately 400 feet apart. For a 400-ft. spacing, transverse fractures may propagate away from the previous fracture with a small angle of deviation from the orthogonal path (less than 2°), as expected from the angle of stress reorientation profile shown inFIG. 1 (simulated using the same parameters). When the spacing is reduced to 300 ft. (as shown inFIG. 7 ), the average angle of deviation from the orthogonal path increases to a little over 5° (e.g., afterstage 4, the average angle of deviation converges toward a value θf=5.7°). A closer look at the fracture trajectory shows that afterfracture 5, the fracture may initially propagate toward the previous fracture and then at some distance from the wellbore, the fracture may start propagating away from the previous fracture. Plotting the stress distribution around an oblique fracture reveals the explanation behind this non-trivial trend as shown inFIG. 8 . -
FIG. 8 illustrates the stress distribution of arock formation 800 with the stress distribution being influenced by the propagation offracture 4 ofFIG. 7 . From the stress redistribution caused byfracture 4 it may be possible to draw astress reversal zone 804, a zone where the subsequent fracture (e.g., afracture 5 ofFIG. 7 ) may be attracted by the previous fracture (e.g., an attraction zone 806), and another zone where the subsequent fracture may propagate away from the previous fracture (e.g., a repulsion zone 808). In some cases, subsequent transverse fractures may propagate in both zones as shown byfracture 5. The size of attraction zone 806 may be function of the net pressure, the in-situ stress contrast and the average angle of deviation from the orthogonal path offracture 4. In the present example, for fracture spacings lower than 400 ft., the initiation point offracture 5 may be located within attraction zone 806 caused by the propagation offracture 4, thusfracture 5 may initially propagate back towardfracture 4 until it leaves attraction zone 806. -
FIG. 9 illustrates fracture propagation with the fracture spacing reduced to 250 ft. in accordance with the present example. Because of the closer spacing, the amount of fracture deviation is larger. For instance,fractures fractures fractures fractures - Additionally, it may be noted that to calculate the trajectory of
fractures fracture 3 may intersect withfracture 2, the stress distribution that may affect thefracture 4 may be modeled based on both thefractures fractures FIG. 10 . -
FIGS. 11 and 12 also illustrate fracture propagation as calculated for a 200-ft. and a 150-ft. spacing respectively of the present example simulation. In those examples, the “unsuccessful” fractures (e.g.,fractures FIG. 11 andfractures FIG. 12 ) may not only intersect the previous fracture but may propagate longitudinally to the horizontal well such that increased hydrocarbon production may not be realized through the inducement of these fractures. For such small values of the fracture spacing, the unsuccessful fractures may be caused by initiating the fractures within the stress reversal region of the previous fracture, which is located inside the attraction zone associated with the previous fracture (e.g., attraction zone 806 ofFIG. 8 ). In the present example (150-ft. spacing), only every other fracturing stage effectively stimulates the shale, thus possibly leaving significant portions of the reservoir inadequately drained. -
FIGS. 13 a and 13 b illustrate the impact of fracture spacing on the angle of deviation of the fractures from the orthogonal path. Below a critical value of the fracture spacing, the efficiency of fracturing stages may be negatively affected as shown by the large variations in deviation angles with respect tospacings -
FIG. 14 illustrates the impact of fracture spacing on the evolution of the net closure stress. As shown inFIG. 14 , for fracture spacings of 400 ft. and 300 ft., the net closure stress only increases with each new stage until reaching a plateau. However for the 250, 200 and 150 ft. fracture spacings, the net pressure may have an up and down trend. - Counting the number of times the net fracturing pressure decreases from one stage to another, may indicate the number of unsuccessful fracture stages identified in
FIGS. 10 , 12, 13 a and 13 b. The decrease in the fracture closure stress (from one stage to another) may be a consequence of the smaller mechanical stress interference (stress shadow) generated by the previous fracture when propagating into stimulated regions of the reservoir instead of orthogonal to the well. - Therefore, as an example, in the case of the smallest spacing, while the designed value is 150 ft., the effective spacing may only equal to 300 ft., as every other fracture may be longitudinal with respect to the wellbore. Accordingly, doubling the number of stages for 150 ft. spacing compared to the 300 ft. spacing may grant very little improvement in well production.
- Thus, as shown above, modeling deviation from the orthogonal path for fractures may reveal a new up-and-down trend in the evolution of the net closure stress. This up and down trend may indicate that the spacing between fractures may be too close to generate any improvement in well production. Therefore, the net closure stress at various spacings may be analyzed to determine the closest spacing that may not yield an up and down net pressure such that optimal spacing of fractures may be determined. Additionally, to determine the proper net closure stress, the propagation direction of each fracture may be estimated instead of assuming that the propagation direction is orthogonal to the well as is traditionally done.
- Fracture spacing may also be determined by analyzing the stress reversal region associated with a previous fracture and by initiating the subsequent fracture outside of the stress reversal region. For example,
FIGS. 15 a and 15 b illustrate the horizontal stress of arock formation 1500 with astress reversal region 1502 associated with a fracture (n+1). Fracture (n+1) may run along the vertical axis ofFIGS. 15 a and 15 b and awellbore 1504 may run along the horizontal axis. As shown inFIGS. 15 a and 15 b,stress reversal region 1502 may extend approximately 230 ft. fromfracture 1502 alongwellbore 1504 as shown by isotropic point s90°. Accordingly a subsequent fracture (e.g., a fracture (n+2)) may not be initiated closer than 230 ft. from fracture (n+1) because, as mentioned above, the subsequent fracture may propagate parallel withwellbore 1504 and may not increase hydrocarbon production fromwellbore 1504. - Further,
FIG. 15 b illustrates that in the present example at point s10°, (e.g., approximately 430 ft. from fracture (n+1)) the stress reorientation ofrock formation 1500 may be 10° and at point s5° (e.g., 600 ft. from fracture (n+1)) the stress reorientation ofrock formation 1500 may be 5°. The stress reorientations of 10° and 5° may be such that a subsequent fracture (e.g., a fracture (n+2)) initiated between points s10° and s5° may not intersect fracture (n+1) although the subsequent fracture may deviate somewhat from an orthogonal path due to the stress reorientation. - The above example illustrates how analyzing the size of the stress reversal region may be used to determine the spacing of fractures when the fractures are initiated consecutively, however, the spacing of fractures initiated alternately may also be determined by analyzing the stress reversal region associated with fractures.
FIGS. 16 a and 16 b illustrate the differences between performing consecutive fracturing and alternate fracturing. InFIG. 16 a it can be seen that each fracture starting with fracture “1” may be initiated one after another in a consecutive order. However, inFIG. 16 b it can be seen that two fractures may be initiated consecutively (e.g., fractures “1” and “2” ofFIG. 16 b), however the two previous fractures may be sufficiently far apart that a third fracture (e.g., fracture “3” ofFIG. 16 b) may be initiated between the two previous fractures, such that the fractures alternate. -
FIGS. 17 a and 17 b illustrate the stress orientation of arock formation 1700 withstress reversal regions FIG. 17 b illustrates that the distance betweenstress reversal regions stress reversal regions FIGS. 17 a and 17 b illustrate that the spacing of fractures “1” and “2” may be determined such thatstress reversal regions - Additionally, by analyzing the stress profile of
rock formation 1700 due to fractures “1” and “2,” it can be seen that the stress reorientation caused by fractures “1” and “2” may substantially cancel each other out such that fracture “3” may propagate in a substantially orthogonal path equidistant from fractures “1” and “2”. Accordingly, the advantages of alternate fracturing may be further illustrated and supported by analyzing the stress reversal regions. - The impact of fracture sequencing may also affect fracture complexity. Hydraulic fracture interaction with pre-existing natural fractures may be a function of a term called the relative net pressure Rn. This parameter may be inversely proportional to the local deviatoric stress in which the fracture propagates as shown below in Equation (6).
-
- High values of the relative net pressure Rn may favor fracture path complexity. Thus, a hydraulic fracture propagating in a region of low stress contrast may create larger networks of interconnected fractures. By calculating the local stress contrast experienced by a propagating fracture, the propensity of the alternate fracturing sequence to generate fracture complexity may be quantified and compared to the more conventional fracturing approach. The average value of the stress contrast seen by a propagating middle fracture in the alternate fracturing sequence may be measured for different values of the spacing between the outside fractures (2 s f). In the present example,
FIG. 18 shows that the spacing for which the stress contrast is lowest may be equal to the minimum fracture spacing previously calculated (325 ft). Thus, the minimum fracture spacing in the alternate fracturing sequence may also be the optimum case for creating fracture complexity. - A comparison of the local stress contrast seen by a fracture along its direction of propagation, in the consecutive and alternate fracture sequence, demonstrates improvement in generating fracture complexity using alternate fracturing versus consecutive fracturing.
FIG. 19 illustrates that the deviatoric stress may approach zero in the near wellbore region in the case of the optimum spacing in the alternate fracturing sequence (325 ft). In the present example, along the first half of propagation, the stress contrast may remain lower than 10 psi, which may be equal to 10% of the in-situ stress contrast. It is only in the second half of the fracture propagation that the local stress contrast increases significantly. Thus, choosing the alternate fracturing sequence, may result in high fracture complexity in the near-wellbore region as a result of the propagation of the “middle fracture”. - Analysis of the stress experienced by rock formations may also be used to determine fracture spacing with respect to multiple horizontal lateral wells.
FIG. 20 illustrates an example of fracture spacing that may be done with multiple horizontal lateral wells.FIG. 20 illustrates three horizontal wells (HW1, HW2 and HW3) that may run substantially parallel to each other through a hydrocarbon reservoir in a horizontal plane that may be substantially parallel with the ground. - The wells of
FIG. 20 may be described by variables representing fracture dimensions (Lf, hf), fracture spacing (sf) and the inter-well spacing (sw). In the present example, the middle well (HW2) may be used to propagate a fracture (e.g., fracture “3”) in between two pairs of fractures previously initiated from the outside wells (e.g., fractures “1” and “2” of well HW1 and fractures “1′” and “2′” of well HW3). The same strategy may be adopted in any horizontal completions having an uneven number of laterals (and of course more than just one lateral). Such strategy may allow for benefiting from the propagation of a “middle fracture”, like in alternate fracturing completions, without the need for special downhole tools. Indeed, in each lateral well (e.g., HW1, HW2 and HW3) the fractures may be initiated in a conventional consecutive sequence. - The spacing between fractures in such multi-lateral sequences may be determined by analyzing the stress distribution (e.g., stress reversal regions) associated with the fractures. For example, the stress distribution between two pairs of fractures (e.g., fractures “1” and “2” of well HW1 and fractures “1′” and “2′” of well HW3 of
FIG. 20 ) propagated from the outside laterals HW1 and HW3 is shown inFIGS. 21 a and 21 b. In the present example, for the reservoir properties and fracture geometry of Table 1, it may be determined based on the stress orientation and distribution caused by fractures “1” and “2” of well HW1 and fractures “1′” and “2,′” that a fracture spacing sf associated with the fractures of each well (e.g., spacing between fractures “1”, “2”, “4”, “6”, “8”, etc., of well HW1, fractures “3”, “5”, “7”, “9”, etc., of well HW and fractures “1′”, “2′”, “4′”, “6′”, “8′”, etc., of well HW3) may be equal to 600 ft. and a well spacing (sw) between wells HW1, HW2 and HW3 may be approximately equal to 500 ft. The above spacing may be determined by analyzing the direction of maximum horizontal stress associated with the fractures. For example, the direction of maximum horizontal stress may be reversed everywhere along the outside laterals as shown inFIGS. 21 a and 21 b. Thus, the outside fractures “1” and “2” are too closely spaced to allow propagation of a transverse fracture from the outside laterals HW1 and HW3, similarly to the alternate fracturing sequence. - When considering refracturing the center lateral, the direction of maximum horizontal stress may still allow propagation of a transverse fracture. For example, the distance of transverse propagation, Ltransverse, of fracture “3” of HW2 may be at a maximum at mid-distance from the previous fractures and may be function of not only the spacing between the outside fractures but also the inter-well spacing (sw). The zone of transverse fracture propagation can also be identified when plotting the angle of stress reorientation as shown in
FIG. 21 b. -
FIG. 22 illustrates the relationship between the length of a “middle fracture” propagating from a center lateral well (e.g., fracture “3” of HW2 with a length Ltransverse) for different values of the fracture length with respect to the inter-well spacing. It can be seen that if the wells are spaced too close to each other, the opportunity to propagate a transverse middle fracture from the center lateral well may not exist at all (e.g., sw/Lf=0.1) because the fracture may very quickly intersect with the other lateral wells bordering the center lateral well. For example, if HW1 and HW3 are sufficiently close to HW2 with respect to the length of fracture “3” inFIG. 20 , fracture “3” may quickly intersect with at least one other fracture and/or well. Therefore, the length of transverse fracture propagation may increase with inter-well spacing and may reach its maximum value when the inter-well spacing is at least equal to the fracture length (e.g., sw/Lf=2). - Ltransverse may also increase with the spacing between the outside fractures (sf). Transverse fracture propagation may not be affected if the fracture spacing is at least equal to twice the minimum fracture spacing in the alternate fracturing sequence (2sf=650 ft). In this case, the stress reorientation angle may be equal to zero everywhere along a line equidistant from the outside fractures.
-
FIG. 23 illustrates the local stress contrast that may be recorded along the assumed propagation direction of a middle fracture (e.g., fracture “3” ofFIG. 20 ). This quantity may be minimum for the minimum possible inter-well spacing (sw/Lf=1) and may also be more sensitive to the inter-well spacing than to the fracture spacing. Thus, it may be advantageous to position the horizontal laterals close to each other, but not closer than a distance equal to the fracture half-height. Otherwise, the benefit of propagating long transverse fractures may be lost. Indeed, that may result in fracturing zones of the reservoir that are already stimulated. - Looking back at
FIG. 21 a, the distance of transverse fracture propagation may be sensitive to the fracture spacing when the inter-well spacing is small. Transverse propagation length may be decreased by over 50% as the fracture spacing decreases from 650 ft. to 600 ft., which is only a 50 ft. spacing differential. Therefore, in the present example and similarly to the case of the alternate fracturing sequence in a single well, the spacing between the outside fractures may be at least be equal to 650 ft. - Finally, the optimum multi-lateral completion strategy in the present example of a typical Barnett shale gas well may be summarized below:
-
- The predicted values of the transverse fracture propagation and average stress contrast for the middle fracture are:
-
- We can finally note that while a 650-ft. spacing may not be practical in some alternate fracturing sequence (e.g., when the refracturing interval may only be 20-ft. wide), this spacing may suffice in a multi-lateral completion. In the latter case, the middle fracture may be initiated from the middle well (and not from the outside well), where the refracturing interval is wide enough to allow fracture initiation from multiple perforation clusters.
- Therefore, by analyzing the stress reorientation regions of rock formations due to fracturing operations, the spacing of the fractures may be determined to improve production from wells, while also improving the efficiency of each fracturing operation. Such stress reorientation analysis may be used for consecutive fracturing, for alternate fracturing and/or for multiple horizontal fracturing operations.
- Further, the in-situ stress contrast, which is the difference between the maximum horizontal stress and the minimum horizontal stress, may influence the stress interference created by multiple consecutive fractures, including fracture intersection. As a result, the evolution of the fracturing pressures during multi-stage fracturing of horizontal wells may be impacted by the in-situ stress contrast, just like it is impacted by the fracture spacing (e.g., as shown in
FIG. 14 ). - Although the minimum horizontal stress may be easily obtained from a mini-frac test, the maximum horizontal stress may be more difficult to evaluate in the field. Knowing the value of the maximum horizontal stress may prove useful in modeling multiple engineering problems in the oil and gas industry, including hydraulic fracturing and wellbore stability and sand production issues.
- The proposed method may be used to calculate the evolution of the net closure stress in a given well for different values of the maximum horizontal stress. By comparing the calculated pressure profiles to the field-measured fracturing pressures, the value of the maximum horizontal stress may be determined for the well in question.
- Modifications, additions and omissions may be made to the above FIGURES without departing from the scope of the present disclosure. For example, the above models and FIGURES have been described with respect to specific rock properties and fracture sizes for illustrative purposes only. The principles described above may be used for any other suitable rock formation.
- Additionally, it is also understood that the stress redistribution of a rock formation caused by propagating a fracture may also be a function of the induced fracture length, fracture width, fluid rheology and the injection rates associated with propagating the fracture. As mentioned above, the propagation of subsequent fractures may be a function of the stress redistribution caused by previous fractures. Therefore, the analysis described above may also be used to determine one or more of the above mentioned properties to better improve fracturing efficiency. For example, in some instances for a particular fracture size, the determined optimal spacing may be too far apart. Accordingly, the spacing may be set at a fixed value and another factor that may affect stress reorientation (e.g., fracture width) may be modified. The stress reorientation, and the propagation and net closure stress of consecutive fractures may be calculated for different values of the fracture width such that an optimum width of the fractures may be determined.
Claims (21)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/531,031 US20120325462A1 (en) | 2011-06-24 | 2012-06-22 | Method for Determining Spacing of Hydraulic Fractures in a Rock Formation |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201161501003P | 2011-06-24 | 2011-06-24 | |
US13/531,031 US20120325462A1 (en) | 2011-06-24 | 2012-06-22 | Method for Determining Spacing of Hydraulic Fractures in a Rock Formation |
Publications (1)
Publication Number | Publication Date |
---|---|
US20120325462A1 true US20120325462A1 (en) | 2012-12-27 |
Family
ID=47360733
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/531,031 Abandoned US20120325462A1 (en) | 2011-06-24 | 2012-06-22 | Method for Determining Spacing of Hydraulic Fractures in a Rock Formation |
Country Status (4)
Country | Link |
---|---|
US (1) | US20120325462A1 (en) |
CN (1) | CN103733091A (en) |
CA (1) | CA2852881A1 (en) |
WO (1) | WO2012178026A2 (en) |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014140752A1 (en) * | 2013-03-15 | 2014-09-18 | Ormat Technologies Inc. | Method for selecting the location of a simulating geothermal well |
WO2014158428A1 (en) * | 2013-03-14 | 2014-10-02 | Halliburton Energy Services, Inc. | Determining a target net treating pressure for a subterranean region |
US8967262B2 (en) * | 2011-09-14 | 2015-03-03 | Baker Hughes Incorporated | Method for determining fracture spacing and well fracturing using the method |
US20150075784A1 (en) * | 2013-09-19 | 2015-03-19 | Shell Oil Company | Phased stimulation methods |
US20150212235A1 (en) * | 2012-09-12 | 2015-07-30 | Bp Exploration Operating Company Limited | System and method for determining retained hydrocarbon fluid |
CN105422070A (en) * | 2015-12-07 | 2016-03-23 | 西南石油大学 | Method for optimizing fracture position of ultra-low permeability heterogeneous gas reservoir fractured horizontal well |
US9297250B2 (en) | 2013-03-14 | 2016-03-29 | Halliburton Energy Services, Inc. | Controlling net treating pressure in a subterranean region |
US20160138339A1 (en) * | 2013-08-09 | 2016-05-19 | Landmark Graphics Corporation | Regression relationship approaches |
US20160265331A1 (en) * | 2011-11-04 | 2016-09-15 | Schlumberger Technology Corporation | Modeling of interaction of hydraulic fractures in complex fracture networks |
EP3019693A4 (en) * | 2013-07-08 | 2017-03-01 | Services Pétroliers Schlumberger | Determining geomechanics completion quality |
WO2017172524A1 (en) * | 2016-03-28 | 2017-10-05 | Baker Hughes Incorporated | Completion optimization process based on acoustic logging data in a lateral section of a horizontal well |
US9869170B2 (en) * | 2015-03-17 | 2018-01-16 | Halliburton Energy Services, Inc. | Methods of controlling water production in horizontal wells with multistage fractures |
WO2018102271A1 (en) * | 2016-11-29 | 2018-06-07 | Conocophillips Company | Methods for shut-in pressure escalation analysis |
US10352145B2 (en) | 2011-03-11 | 2019-07-16 | Schlumberger Technology Corporation | Method of calibrating fracture geometry to microseismic events |
US10422208B2 (en) | 2011-11-04 | 2019-09-24 | Schlumberger Technology Corporation | Stacked height growth fracture modeling |
US10544667B2 (en) | 2011-11-04 | 2020-01-28 | Schlumberger Technology Corporation | Modeling of interaction of hydraulic fractures in complex fracture networks |
US10801307B2 (en) | 2016-11-29 | 2020-10-13 | Conocophillips Company | Engineered stress state with multi-well completions |
US10954774B2 (en) | 2013-12-18 | 2021-03-23 | Conocophillips Company | Method for determining hydraulic fracture orientation and dimension |
US20210131250A1 (en) * | 2019-11-04 | 2021-05-06 | Chevron U.S.A. Inc. | Wellbore to fracture connectivity |
US11028679B1 (en) | 2017-01-24 | 2021-06-08 | Devon Energy Corporation | Systems and methods for controlling fracturing operations using monitor well pressure |
CN113530499A (en) * | 2020-04-21 | 2021-10-22 | 中国石油天然气股份有限公司 | Perforating method and perforating device for shale oil-gas horizontal well with small well spacing |
US11209558B2 (en) | 2018-05-09 | 2021-12-28 | Conocophillips Company | Measurement of poroelastic pressure response |
US11365617B1 (en) | 2017-01-24 | 2022-06-21 | Devon Energy Corporation | Systems and methods for controlling fracturing operations using monitor well pressure |
US11859490B2 (en) | 2021-08-19 | 2024-01-02 | Devon Energy Corporation | Systems and methods for monitoring fracturing operations using monitor well flow |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103926620A (en) * | 2014-05-09 | 2014-07-16 | 南京大学 | Hydraulic fracturing monitoring method based on array deconvolution treatment |
CN105178952B (en) * | 2015-09-09 | 2018-04-06 | 中国石油天然气股份有限公司 | Determine the method and device of horizontal well man-made fracture spacing |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5482116A (en) * | 1993-12-10 | 1996-01-09 | Mobil Oil Corporation | Wellbore guided hydraulic fracturing |
US5540093A (en) * | 1994-11-15 | 1996-07-30 | Western Atlas International | Method for optimizing the alignment of a horizontal borehole relative to the strike of rock-layer stress planes |
CA2475007A1 (en) * | 2002-02-01 | 2003-08-14 | Regents Of The University Of Minnesota | Interpretation and design of hydraulic fracturing treatments |
CN1991127A (en) * | 2005-12-30 | 2007-07-04 | 许靖华 | Method for improving crude oil recovery efficiency through three-dimensional water injection between horizontal wells having two human crack surfaces |
AU2008299195B2 (en) * | 2007-09-13 | 2012-08-23 | M-I Llc | Method of using pressure signatures to predict injection well anomalies |
CA2711773C (en) * | 2008-01-31 | 2013-03-19 | Schlumberger Canada Limited | Method of hydraulic fracturing of horizontal wells, resulting in increased production |
US8498852B2 (en) * | 2009-06-05 | 2013-07-30 | Schlumberger Tehcnology Corporation | Method and apparatus for efficient real-time characterization of hydraulic fractures and fracturing optimization based thereon |
-
2012
- 2012-06-22 CN CN201280037328.XA patent/CN103733091A/en active Pending
- 2012-06-22 WO PCT/US2012/043773 patent/WO2012178026A2/en active Application Filing
- 2012-06-22 US US13/531,031 patent/US20120325462A1/en not_active Abandoned
- 2012-06-22 CA CA2852881A patent/CA2852881A1/en not_active Abandoned
Non-Patent Citations (4)
Title |
---|
Addis, M.A. et al., "Estimation of Horizontal Stresses at Depth in Faulted Regions and Their Relationship to Pore Pressure Variations", March 1996, SPE/ISRM Rock Mechanics in Petroleum Engineering Conference, Society of Petroleum Engineers. * |
Li, S. et al., "Maximum Horizontal Stress and Wellbore Stability While Drilling: Modeling and Case Study", December 1 - 3, 2010, SPE Latin American & Caribbean Petroleum Engineering Conference, Society of Petroleum Engineers. * |
Roussel, N.P. et al., "Optimizing Fracture Spacing and Sequencing in Horizontal Well Fracturing", February 10 - 12, 2010, SPE Intenational Symposium and Exhibition on Formation Damage Control, Society of Petroleum Engineers. * |
Sharma, Mukul M. et al., "Improved Reservoir Access Through Refracture Treatments in Tight Gas and Gas Shales", April 7, 2010, Univeristy of Texas at Austin, Petroleum and Geosystems Engineering. * |
Cited By (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10352145B2 (en) | 2011-03-11 | 2019-07-16 | Schlumberger Technology Corporation | Method of calibrating fracture geometry to microseismic events |
US8967262B2 (en) * | 2011-09-14 | 2015-03-03 | Baker Hughes Incorporated | Method for determining fracture spacing and well fracturing using the method |
US10422208B2 (en) | 2011-11-04 | 2019-09-24 | Schlumberger Technology Corporation | Stacked height growth fracture modeling |
US10544667B2 (en) | 2011-11-04 | 2020-01-28 | Schlumberger Technology Corporation | Modeling of interaction of hydraulic fractures in complex fracture networks |
US20160265331A1 (en) * | 2011-11-04 | 2016-09-15 | Schlumberger Technology Corporation | Modeling of interaction of hydraulic fractures in complex fracture networks |
US20150212235A1 (en) * | 2012-09-12 | 2015-07-30 | Bp Exploration Operating Company Limited | System and method for determining retained hydrocarbon fluid |
US10012764B2 (en) * | 2012-09-12 | 2018-07-03 | Bp Exploration Operating Company Limited | System and method for determining retained hydrocarbon fluid |
WO2014158428A1 (en) * | 2013-03-14 | 2014-10-02 | Halliburton Energy Services, Inc. | Determining a target net treating pressure for a subterranean region |
US9217318B2 (en) | 2013-03-14 | 2015-12-22 | Halliburton Energy Services, Inc. | Determining a target net treating pressure for a subterranean region |
US9297250B2 (en) | 2013-03-14 | 2016-03-29 | Halliburton Energy Services, Inc. | Controlling net treating pressure in a subterranean region |
WO2014140752A1 (en) * | 2013-03-15 | 2014-09-18 | Ormat Technologies Inc. | Method for selecting the location of a simulating geothermal well |
US9784090B2 (en) | 2013-03-15 | 2017-10-10 | Ormat Technologies Inc. | Method for selecting the location of a stimulating geothermal well |
EP3019693A4 (en) * | 2013-07-08 | 2017-03-01 | Services Pétroliers Schlumberger | Determining geomechanics completion quality |
US20160138339A1 (en) * | 2013-08-09 | 2016-05-19 | Landmark Graphics Corporation | Regression relationship approaches |
US10344535B2 (en) * | 2013-08-09 | 2019-07-09 | Landmark Graphics Corporation | Regression relationship approaches |
CN104453822A (en) * | 2013-09-19 | 2015-03-25 | 国际壳牌研究有限公司 | Phased stimulation method for subterranean formation having well bore therein |
US20150075784A1 (en) * | 2013-09-19 | 2015-03-19 | Shell Oil Company | Phased stimulation methods |
US10954774B2 (en) | 2013-12-18 | 2021-03-23 | Conocophillips Company | Method for determining hydraulic fracture orientation and dimension |
US11725500B2 (en) | 2013-12-18 | 2023-08-15 | Conocophillips Company | Method for determining hydraulic fracture orientation and dimension |
US11371339B2 (en) | 2013-12-18 | 2022-06-28 | Conocophillips Company | Method for determining hydraulic fracture orientation and dimension |
US9869170B2 (en) * | 2015-03-17 | 2018-01-16 | Halliburton Energy Services, Inc. | Methods of controlling water production in horizontal wells with multistage fractures |
CN105422070A (en) * | 2015-12-07 | 2016-03-23 | 西南石油大学 | Method for optimizing fracture position of ultra-low permeability heterogeneous gas reservoir fractured horizontal well |
US10364672B2 (en) | 2016-03-28 | 2019-07-30 | Baker Hughes, A Ge Company, Llc | Completion optimization process based on acoustic logging data in the lateral section in a horizontal well |
WO2017172524A1 (en) * | 2016-03-28 | 2017-10-05 | Baker Hughes Incorporated | Completion optimization process based on acoustic logging data in a lateral section of a horizontal well |
US11280165B2 (en) | 2016-11-29 | 2022-03-22 | Conocophillips Company | Methods for shut-in pressure escalation analysis |
US10753181B2 (en) | 2016-11-29 | 2020-08-25 | Conocophillips Company | Methods for shut-in pressure escalation analysis |
US10801307B2 (en) | 2016-11-29 | 2020-10-13 | Conocophillips Company | Engineered stress state with multi-well completions |
WO2018102271A1 (en) * | 2016-11-29 | 2018-06-07 | Conocophillips Company | Methods for shut-in pressure escalation analysis |
US11028679B1 (en) | 2017-01-24 | 2021-06-08 | Devon Energy Corporation | Systems and methods for controlling fracturing operations using monitor well pressure |
US11365617B1 (en) | 2017-01-24 | 2022-06-21 | Devon Energy Corporation | Systems and methods for controlling fracturing operations using monitor well pressure |
US11131176B1 (en) | 2017-01-24 | 2021-09-28 | Devon Energy Corporation | Systems and methods for controlling fracturing operations using monitor well pressure |
US11500114B2 (en) | 2018-05-09 | 2022-11-15 | Conocophillips Company | Ubiquitous real-time fracture monitoring |
US11921246B2 (en) | 2018-05-09 | 2024-03-05 | Conocophillips Company | Measurement of poroelastic pressure response |
US11209558B2 (en) | 2018-05-09 | 2021-12-28 | Conocophillips Company | Measurement of poroelastic pressure response |
US11466557B2 (en) | 2019-11-04 | 2022-10-11 | Chevron U.S.A. Inc. | Determining proppant distribution for a plurality of clusters within a fracture stage |
US20210131250A1 (en) * | 2019-11-04 | 2021-05-06 | Chevron U.S.A. Inc. | Wellbore to fracture connectivity |
US11674376B2 (en) * | 2019-11-04 | 2023-06-13 | Chevron U.S.A. Inc. | Wellbore to fracture connectivity |
US11359479B2 (en) | 2019-11-04 | 2022-06-14 | Chevron U.S.A. Inc. | Determining a hydraulic fracture completion configuration for a wellbore |
CN113530499A (en) * | 2020-04-21 | 2021-10-22 | 中国石油天然气股份有限公司 | Perforating method and perforating device for shale oil-gas horizontal well with small well spacing |
US11859490B2 (en) | 2021-08-19 | 2024-01-02 | Devon Energy Corporation | Systems and methods for monitoring fracturing operations using monitor well flow |
Also Published As
Publication number | Publication date |
---|---|
WO2012178026A3 (en) | 2013-05-02 |
WO2012178026A2 (en) | 2012-12-27 |
CN103733091A (en) | 2014-04-16 |
CA2852881A1 (en) | 2012-12-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20120325462A1 (en) | Method for Determining Spacing of Hydraulic Fractures in a Rock Formation | |
US10801307B2 (en) | Engineered stress state with multi-well completions | |
US11727176B2 (en) | Methods for shut-in pressure escalation analysis | |
Roussel et al. | Strategies to minimize frac spacing and stimulate natural fractures in horizontal completions | |
EP3271547B1 (en) | Stacked height growth fracture modeling | |
US20160003020A1 (en) | Methods for time-delayed fracturing in hydrocarbon formations | |
US10267131B2 (en) | Competition between transverse and axial hydraulic fractures in horizontal well | |
Roussel et al. | Optimizing fracture spacing and sequencing in horizontal-well fracturing | |
Saberhosseini et al. | Optimization of the horizontal-well multiple hydraulic fracturing operation in a low-permeability carbonate reservoir using fully coupled XFEM model | |
EP2774066B1 (en) | Modeling of interaction of hydraulic fractures in complex fracture networks | |
Taghichian et al. | Stress shadow size and aperture of hydraulic fractures in unconventional shales | |
US20160265331A1 (en) | Modeling of interaction of hydraulic fractures in complex fracture networks | |
AU2019200654A1 (en) | Modeling of interaction of hydraulic fractures in complex fracture networks | |
Manchanda et al. | Factors influencing fracture trajectories and fracturing pressure data in a horizontal completion | |
Manchanda et al. | Time-delayed fracturing: a new strategy in multi-stage, multi-well pad fracturing | |
Algarhy et al. | Increasing hydrocarbon recovery from shale reservoirs through ballooned hydraulic fracturing | |
Marder et al. | Simple models of the hydrofracture process | |
Manchanda et al. | Simulating the life of hydraulically fractured wells using a fully-coupled poroelastic fracture-reservoir simulator | |
Jabbari et al. | Hydraulic fracturing design for horizontal wells in the Bakken formation | |
Yang et al. | Effect of depletion-induced stress reorientation on infill well fracture propagation | |
Ghassemi | Application of rock failure simulation in design optimization of the hydraulic fracturing | |
Sesetty et al. | Modeling dense-arrays of hydraulic fracture clusters: Fracture complexity, net pressure and model calibration | |
Han et al. | Geomechanical, Geological, and Engineering Controls of Hydraulic Fracturing | |
CN104712299B (en) | It is adapted to gas well control water and increases the design method that air pressure splits | |
Ghassemi | Impact of fracture interactions, rock anisotropy and heterogeneity on hydraulic fracturing: some insights from numerical simulations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: BOARD OF REGENTS, THE UNIVERSITY OF TEXAS SYSTEM, Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:ROUSSEL, NICOLAS P.;SHARMA, MUKUL M.;SIGNING DATES FROM 20120711 TO 20120731;REEL/FRAME:028818/0682 |
|
AS | Assignment |
Owner name: ENERGY, UNITED STATES DEPARTMENT OF, DISTRICT OF C Free format text: CONFIRMATORY LICENSE;ASSIGNOR:UNIVERSITY OF TEXAS AT AUSTIN OFFICE OF TECHNOLOGY COMNMERCIALIZATION, THE;REEL/FRAME:033779/0712 Effective date: 20140123 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |