Hydraulic Fracturing - Well Stimulation PEN715 (Amro Shawky)
Hydraulic Fracturing - Well Stimulation PEN715 (Amro Shawky)
Hydraulic Fracturing - Well Stimulation PEN715 (Amro Shawky)
Presented By:
For
Dr. Eissa Shokir
Well Stimulation
Phd - Summer 2020/2021
CONTENT:
Introduction
Reservoir Candidates for Hydraulic Fracturing
Formation Evaluation
Hydraulic Fracturing Treatment Operation
Fracture Initiation Geomechanics
Fracture Stress Shadowing
Fracture Propagation Models
Hydraulic Fracturing Fluids
Proppants and Fracture Conductivity
Treatment Design and Field Calculations
References
2
INTRODUCTION
WHAT IS HYDRAULIC FRACTURING?
TYPES OF HYDRAULIC FRACTURING
TIP SCREEN OUT FRACTURING
3
WHAT IS HYDRAULIC FRACTURING?
It is the process of injecting pressurized fluid into a formation until the formation rock cracks due to tensile
failure, then extending the fracture by continued fluid injection.
A fracking fluid is used to apply a sufficient force to the formation rock until reaching a point where the rock
breaks in tension and cracks develop outwards guided by a set of internal rock stresses. This increases the contact
between the well and the formation natural fractures and billions of tiny pores within the rock.
A solid proppant is injected into the fractured formation to prevent the fracture from closing. This creates a high
permeability flow path for reservoir fluids to be rapidly produced
Applications of hydraulic fracturing:
Hydraulic fracturing can be categorized into 3 groups according to the process applied to the target formation to
induce the fracture [2] :
1. Hydraulic Fracturing:
• Low rate of pressure loading
• Bidirectional fracture oriented perpendicular to the minimum principal rock stress
• Fracture penetration or reach can be extensive (hundreds of feet)
2. Explosive Fracturing:
• Rapid pressurization of formation
• Results in highly fractured zone around the wellbore in a radial pattern
• Small penetration or reach (approximately 10 feet)
• Peak pressure exceeds both minimum and maximum horizontal in situ stresses
• Near-wellbore stimulation is the primary objective
5
TYPES OF HYDRAULIC FRACTURING (CONT.)
6
TIP SCREEN OUT FRACTURING
Proppant
Concentration 0.5 – 2.0 4 – 12 Figure 3: Comparison between
(lb/ft2) Conventional and TSO hydraulic
fracturing. [2]
8
TIP SCREEN OUT FRACTURING (CONT.)
9
RESERVOIR CANDIDATES FOR HYDRAULIC
FRACTURING
RESERVOIR CANDIDATES
UNCONVENTIONAL RESERVOIR CANDIDATES
CONVENTIONAL RESERVOIR CANDIDATES
10
RESERVOIR CANDIDATES
Hydraulic fracturing reservoir candidates are ones with a substantial volume of OOIP
and/or OGIP that might:
• Suffer from low permeability
• Suffer from reduced flow rates due to near-wellbore damage
• Require a connection between natural fractures and/or cleats in the formation and the producing
well
• Require reduction in pressure drop around the wellbore to reduce sand production
• Require reduction in pressure drop around the wellbore to minimize problems with asphaltenes and/
or paraffin deposition
• Increase the drainage area around the wellbore by increasing contact with the reservoir
11
RESERVOIR CANDIDATES
Reservoirs in need of hydraulic fracturing can be categorized into conventional and unconventional reservoirs.
• Conventional reservoirs are those in which wells can be drilled so that oil and natural gas can be produced at economic flow
rates without large stimulation treatments or any special recovery process. [4]
• On the other hand, unconventional reservoirs are those that cannot be produced at economic flow rates or that does not
produce economic volumes of oil and gas without the assistance of massive stimulation treatments or special recovery
processes and technologies. [4]
14
UNCONVENTIONAL RESERVOIR CANDIDATES (CONT.)
15
CONVENTIONAL RESERVOIR CANDIDATES
16
FORMATION EVALUATION
FORMATION EVALUATION PARAMETERS
WELL LOGGING ANALYSIS
WELL TESTING
CORE SAMPLE ANALYSIS
17
FORMATION EVALUATION PARAMETERS
18
FORMATION EVALUATION PARAMETERS (CONT.)
Petrophysical data required for candidate selection and hydraulic fracture modeling
can be obtained by integrating data from drilling information, well testing, core
sample testing and well logging data analysis. [6]
This data can be divided into two groups.
1. Controllable parameters: mostly related to well design and injection of proppant laden fracture
fluid and they include well landing zone, well orientation, fracture spacing, well completion,
perforation design, fracture fluid and proppant type and volume selection, pump rate, proppant
concentration, etc…
2. Uncontrollable parameters: has the greatest effect on the hydraulic fracturing operation and
includes reservoir geology, in situ stress magnitudes and orientations, formation pressure, natural
fractures, etc…
19
WELL LOGGING DATA ANALYSIS
• True formation porosity is determined using porosity values from both neutron and density logs
20
WELL LOGGING DATA ANALYSIS (CONT.)
• To calculate the water saturation of formations, we can use the total resistivity readings obtained from
resistivity logs
• Where is the formation resistivity factor of clean sand and is the saturation exponent.
21
WELL LOGGING DATA ANALYSIS (CONT.)
Natural
Fracture Characteristics: [1]
• Existence of natural fractures in a reservoir add to the complexity of hydraulic fracture modeling. Natural
fractures affect well log interpretation, the geometry and propagation of hydraulic fractures and fluid leak-off
calculations
1. Estimating fracture departure width using dual laterolog data and mud filtrate resistivity (): [7]
22
WELL LOGGING DATA ANALYSIS (CONT.)
2. Fracture
orientation can be determined using
data from dipmeter and image logs: [8]
23
WELL LOGGING DATA ANALYSIS (CONT.)
4. Fracture porosity can be calculated using integration of data from well testing results, image log processing and
core data analysis: [10]
24
WELL LOGGING DATA ANALYSIS (CONT.)
Dynamic
Young’s Modulus and Poisson’s Ratio:
• To create a continuous profile of the rock elastic properties along the length of the
region of interest, first, we need to calculate the rock dynamic elastic properties
• They can be calculated using correlations that utilize data from acoustic logs like
the correlations shown below: [11]
• The values of dynamic elastic properties are then calibrated by correlating them
with static dynamic elastic properties obtained from core sample analysis to
provide us with accurate and realistic information
25
WELL LOGGING DATA ANALYSIS (CONT.)
• For offshore areas with water depth, you must consider the air gap and water depth, so the overburden stress will be:
26
WELL LOGGING DATA ANALYSIS (CONT.)
• The tectonic strains in the minimum and maximum horizontal stress directions ( and ) are calculated by the following
equations:
Where:
• and are calculated using minimum horizontal stress values from LOT or mini-frac tests from at least two depth intervals
(depth I and depth i+1), then we find the average of the values and .
27
WELL LOGGING DATA ANALYSIS (CONT.)
29
WELL TESTING (CONT.)
Static
Young’s Modulus and Poisson Ratio:
• Static Young’s modulus and Poisson ratio can be obtained by
performing triaxial compression tests on core samples
• A cylindrical core sample with a 2:1 ratio of length to diameter
is loaded axially at a constant confining pressure which presents
the in situ stress conditions
• The axial stress, axial and lateral strains are measured and
recorded during the test and used to calculate the static Young’s
modulus and Poisson ratio
• The static values are then used to calibrate dynamic Young’s
modulus and Poisson ratio obtained from well logging data to
create a continuous rock elastic properties profile along the
length of the formations of interest
Brittleness Index: [12]
• It is calculated using data from triaxial compression tests on
core samples Figure 12: Triaxial compression test apparatus.
[13]
31
CORE SAMPLE ANALYSIS (CONT.)
Fracture
Toughness and Critical Stress Intensity Factor: [1]
• It is a parameter that controls the magnitude of the pressure drop at the tip of a fracture
• Fracture toughness can be determined based on the principle of Linear Elastic Fracture Mechanics (LEFM) which assumes that
when the stresses near the crack tip exceed the material fracture toughness, the crack will grow.
• The crack tip stress field is a function of the location, loading, and geometry
• The fracture toughness or critical stress intensity factor is material specific and can be obtained by experiments performed on core
samples of known crack length
• If the stress intensity factor () is equal to or greater than the critical stress intensity factor (), the crack is reinitiated and the
fracture grows
• As the critical stress intensity factor increases, the pressure required for fracture propagation also increases
Where: is the difference between the pressure inside the fracture and the closure pressure
is a coefficient which depends on the fracture geometry
is the fracture radius
32
CORE SAMPLE ANALYSIS (CONT.)
34
STAGES OF HYDRAULIC FRACTURING TREATMENT OPERATION
35
STAGES OF HYDRAULIC FRACTURING TREATMENT OPERATION (CONT.)
3. Proppant-Laden Stage:
• The proppant solids are mixed into additional fracturing
fluids and the resulting slurry is pumped into the reservoir
to prop open the created fracture. This can be divided into
several sub-stages with each one having a different
proppant concentration.
4. Flush Stage:
• A specific volume of flush fluid is pumped to clear the
tubulars of proppant and the pumps are shut off.
• Well pressure is bled off to allow the fracture to close on
the proppant.
• After the proppants are set in place, thin fluids with vis-
breakers are injected to reduce the viscosity of the pumped
fluids to allow them to flow.
• Recover the injected fluids by flowing or lifting the well
(load recovery) Figure 14: Stages of hydraulic fracturing.
36
FRACTURE INITIATION GEOMECHANICS
OPEN HOLE COMPLETIONS
CASED AND PERFORATED COMPLETIONS
37
FRACTURE INITIATION GEOMECHANICS
A formation fractures when the pressure applied by the mud weight is too high resulting in tensile failure. This
pressure is known as the Fracture Initiation Pressure (FIP).
Tensile failure occurs when minimum effective stress around the wellbore becomes tensile (negative) and when
this tensile stress is greater than or equal to the formation tensile strength.
Horizontal Fractures:
• Horizontal fractures occur in shallow formations
where the overburden stress is the minimum
principal stress
39
OPEN HOLE COMPLETIONS (CONT.)
Vertical Fractures:
• According to the tensile strength criterion, the fracture initiates at the wellbore wall when a principal tensile
stress exceeds the tensile strength of the rock.
• Vertical fractures occur when the minimum principal stress is horizontal and the maximum principal stress is
vertical. Calculating the principal stresses acting on the wellbore wall, first, requires performing a stress
analysis of the induced forces acting on the wellbore wall in the cylindrical (polar) coordinate system.
Figure 17: Transformation of in situ stresses into polar coordinate system. [15]
40
OPEN HOLE COMPLETIONS (CONT.)
Inclined Wellbores:
• To calculate the induced stresses, Kirsch Elastic solution is used which assumes that the rock is homogeneous and exhibits
linear elastic behavior.
• First step is to express the in situ stresses in the Cartesian Coordinate System:
Next,
the induced stresses are expressed in the polar coordinate system:
Finally, we calculate the principal stresses using the calculated induced stresses in the polar coordinate system:
Substituting the values of , and calculated previously using Kirsch equations into the above equation and solving
for , we can obtain the mud weight or FIP required to initiate a fracture for an inclined wellbore. [16]
The general equation of the fracture initiation angle () is:
43
OPEN HOLE COMPLETIONS (CONT.)
44
OPEN HOLE COMPLETIONS (CONT.)
Vertical Wellbores:
• For vertical wellbores, the FIP can be expressed in term of the minimum (most tensile) tangential stress “”
which occurs perpendicular to the minimum horizontal in situ stress direction ( [17]
• Assuming the reservoir is impermeable and no fluid penetration into the formation, the FIP is expressed as:
Where:
Horizontal Wellbores:
• Because the overburden stress is now aligned along the x-axis direction, the general solution for FIP is:
b) For strike slip or reverse faulting stress regimes, the condition gives the minimum FIP required to initiate horizontal
fractures in the maximum horizontal stress direction is:
46
OPEN HOLE COMPLETIONS (CONT.)
48
CASED AND PERFORATED COMPLETIONS
49
CASED AND PERFORATED COMPLETIONS (CONT.)
For vertical wellbores, the wellbore axis aligns with the overburden stress and the two orthogonal
stresses align with the maximum and minimum horizontal in situ stresses.
Assuming the perforation angle is aligned with the minimum in situ stress direction, the tangential
stress can be calculated at (and is called the longitudinal tangential stress ) as well as at (called the
transverse tangential stress ) at the base of the perforation of the wellbore.
Under these conditions, a longitudinal fracture is expected, thus is expected to initiate the anticipated
fracture:
Where:
50
FRACTURE STRESS SHADOWING
51
FRACTURE STRESS SHADOWING
When a well undergoes a hydraulic fracturing stimulation, multiple fractures are created along the length of the
producing zone to maximize the stimulated reservoir volume. Thus, it is important to consider the fracturing-
induced changes to the in situ stresses around the wellbore
The induced changes in the near-wellbore in situ stresses are known as stress shadow effect and they are
considered to be the primary cause of complex interactions between multiple hydraulic fractures and/or
interactions between hydraulic fractures and natural fractures within the reservoir. [21]
Fracture spacing is a key parameter in controlling the stress shadow effect. Narrow spacing between hydraulic
fractures creates larger stress perturbations resulting in increased and values and a higher FIP.
The increase in the minimum horizontal in situ stress ( ) is at a greater scale than that of the maximum horizontal
in situ stress (.
Strong fracture interaction is possible if the stress perturbation () is greater than the horizontal in situ stress
anisotropy ().
52
FRACTURE STRESS SHADOWING (CONT.)
By measuring the stress perturbation () and stress anisotropy, we can construct a model that relates the fracture
spacing with stress anisotropy to determine the spacing required between multiple fractures as shown in the graph
below.
54
2-D FRACTURE MODELS
• Formations with higher Young’s Modulus require a higher net pressure to achieve the same fracture
dimensions then ones with lower Young’s Modulus assuming all other variables remain unchanged
55
2-D FRACTURE MODELS (CONT.)
Later versions of the PKN model were adjusted to account for fluid loss and storage effects.
The storage dominated approximation of the PKN model: [22]
Fracture width:
Fracture half-length:
High leak-off approximation of the PKN model: [22]
Fracture width:
Fracture half-length:
57
2-D FRACTURE MODELS (CONT.)
Radial Model:
• Assumes that a radial or penny-shaped crack
occurs when the fracture is horizontal or when
there are no stress barriers constraining height
growth of a vertical fracture
• Assumes that the length of the fracture is equal
to its height
• Assuming constant pressure, the width of a
static circular fracture of radius “R” is
expressed as:
58
2-D FRACTURE MODELS (CONT.)
Equationsfor the maximum fracture width calculation during field fracture calculation/design
via 2-D models: (2-D Plane Strain Solution) [23]
PKN Model:
KGD Model:
59
3-D FRACTURE MODELS
1. Cell-based Model:
• The fracture length is divided into cells along the length of the fracture
but the fracture shape is not redefined
• The model assumes that the fracture’s vertical extension is slow and the
vertical fluid flow is ignored
• It also assumes that fracture half length is much greater than fracture
height
• The cross-sectional fracture geometry is determined using the net
pressure (e.g. the fracture height at any point of the fracture cross- Figure 29: Fracture shape
section is calculated using the pressure in that cell) according to the Cell-based Model.
[22]
60
3-D FRACTURE MODELS (CONT.)
2. Lumped-Parameter Model:
• Fracture shape is assumed to consist of two half ellipses of
equal horizontal extent but different vertical extent
• The model can calculate the fracture height, length and
width under the condition of slow vertical spreading
• Fluid flow is generally assumed as streamlines from the
perforated interval at the wellbore to the edge of the fracture
• Fracture half-length and the lower and upper fracture heights
are calculated at every time-step based on lumped
parameters
61
3-D FRACTURE MODELS (CONT.)
62
FRACTURE GEOMETRY MODELING DURING PUMPING
The Nolte-Smith pressure analysis technique is a technique to interpret the fracture geometry created during
pumping by analyzing the pressure response.
It is based on the expected pressure response of the 2-D fracture models and then predicts the pressure response
when certain pressure behavior takes place.
For example, the PKN model assumes that the net pressure is directly proportional to the fracture geometry which
is directly proportional to the root of time
Thus, on a log-log plot, the relationship between the net pressure and time is a straight line with a slope of
If the pressure behavior deviates from this ideal case, other types of fracture growth can be identified
63
FRACTURE GEOMETRY MODELING DURING PUMPING (CONT.)
65
FLUID SELECTION CRITERIA
Optimum fracturing fluid design and selection is governed by the following criteria:
1. Viscosity should be high enough for achieving desired fracture dimensions, proppant transportation and fluid-
loss control
2. Cleanliness after flowback to produce the maximum post-fracture conductivity
3. Compatibility with reservoir fluids to avoid creating emulsions or slugs that would plug the formation
4. Surface pump pressure and friction pressure considerations
5. Compatibility with other materials such as resin coated proppants
6. Safety and environmental concerns
7. Cost-effective, stable and easy to produce
8. Compatible with reservoir rock (no swelling of clays or fines migration)
66
TYPES OF FRACKING FLUIDS
67
TYPES OF FRACKING FLUIDS (CONT.)
68
FRACKING FLUID ADDITIVES
69
PROPPANTS AND FRACTURE CONDUCTIVITY
TYPES OF PROPPANTS
FRACTURE CONDUCTIVITY
PROPPANT TRANSPORT MODELING
70
TYPES OF PROPPANTS
An ideal proppant should be strong, resistant to crushing, resistant to corrosion, has a low density and readily available
at low costs.
There are 3 types of proppants that generally fit those criteria:
1. Silica Sand:
• Obtained from sand mining and tested to ensure it has the necessary compressive strength for fracture operations
• Considered the cheapest per pound compared to the other types
• Generally used to prop fractures open in shallow formations
2. Resin Coated:
• Stronger than silica sand and used where higher compressive strength is required to minimize crushing and permeability reduction
• It is a manufactured proppant and not naturally occurring and is more expensive than sand but with a lower effective density
3. Ceramic Proppants:
• Consist of sintered bauxite, intermediate strength proppants (ISP) and lightweight proppants (LWP)
• The strength of this proppant is directly proportional to its density and cost with sintered bauxite being stronger and more expensive than ISP
and LWP
• It is an engineered and manufactured proppant and can be used to stimulate deep wells where large in situ stress values exist
71
TYPES OF PROPPANTS (CONT.)
Silica Sand:
Figure 32: West Texas Dune Sand Figure 33: Canadian Glacial Sand Figure 34: Arkansas River Sand
[1] [1] [1]
72
TYPES OF PROPPANTS (CONT.)
Figure 35: Curable resin coated Figure 36: Curable resin coated
proppant at standard conditions. [1] proppant under reservoir conditions. [1]
73
TYPES OF PROPPANTS (CONT.)
Ceramic Proppants:
Figure 37: High strength ceramic [1] Figure 38: ISP ceramic Figure 39: LWP ceramic
[1] [1]
74
FRACTURE CONDUCTIVITY
The main function of a proppant is to maintain a conductive, post-operation propped fracture to provide an
effective connection between formulation and wellbore
The potential benefit of the propped fracture depends on the capacity of the formation to deliver fluids to the
fracture and the capacity of the fracture to deliver fluids to the well
This relationship can be quantified using the dimensionless fracture conductivity () parameter:
75
FRACTURE CONDUCTIVITY (CONT.)
Fracture conductivity of the propped fracture decreases with time as the well is put on production due to the
increase of the effective stress acting on the proppant
The effective stress acting on the proppant is defined as the difference between the in situ stress and the flowing
pressure in the fracture
As the well continues to produce over time, the flowing bottomhole pressure will decrease which in turn increases
the effective stress
Proppant properties such as density and size (diameter) must be considered carefully when choosing the
appropriate propping agent fro the job, since they play an important role in proppant transportation and settling in
the fracture besides fracturing fluid viscosity
Proppant transport modeling is represented by:
1. Proppant settling: represents the movement of the proppant grains within the fracture fluid system. It can be expressed by
Stoke’s Law for settling velocity for different Reynold’s numbers:
2. Proppant
convection:
• Represents the mechanism of differential
movement of different treatment stages because of
the density differences and fluid viscosity.
• Higher proppant loadings in later slurries result in
higher slurry densities which may displace lighter
slurries already placed in the fracture
79
PRESSURE CALCULATIONS
80
PRESSURE CALCULATIONS (CONT.)
Closure Pressure ( or ):
• It is the minimum required pressure to keep the fracture open
• Can be assumed to be equal to the minimum in situ horizontal stress ()
or:
82
PRESSURE CALCULATIONS (CONT.)
83
PERFORATION DESIGN
In designing the number of perforations per stage for unconventional reservoirs, the practice of limited entry
design is often used
Limited entry involves limiting the number of perforations in a completion design to increase the perforation
friction pressure and create a choking effect
The choking effect encourages the simultaneous entry of the fracturing fluid into multiple zones of varying in situ
stress states to increase perforation efficiency and the success of the hydraulic fracturing stimulation
Limited entry can be achieved through the following steps:
1. Choose the of a single perforation (usually recommended to be chosen between 200 and 300 psi)
2. Once the friction pressure is chosen, calculate the rate per perforation ():
84
PROPPANT CONCENTRATION TREATMENT SCHEDULE
1. Assume a fracture half-length () to achieve and an injection rate () to pump. Calculate the fracture width () using a
selected fracture propagation model (e.g.: KGD model).
2. Based on material balance, solve for the injection fluid volume () from the following equation:
Where:
85
PROPPANT CONCENTRATION TREATMENT SCHEDULE (CONT.)
Since the fluid loss multiplier () depends on the fluid efficiency () which is unknown in the beginning, a numerical
iteration procedure is required.
86
PROPPANT CONCENTRATION TREATMENT SCHEDULE (CONT.)
4. Proppant weight requirement () in ppg and total weight of proppant required () in pounds:
The propped fracture width depends directly on the slurry proppant concentration (.
87
PROPPANT CONCENTRATION TREATMENT SCHEDULE (CONT.)
Figure 43: Table of proppant concentration injection schedule. Figure 44: Graph of proppant concentration schedule
[22]
during a hydraulic fracturing treatment. [22]
88
SLURRY CONCENTRATION CALCULATIONS
89
SLURRY CONCENTRATION CALCULATIONS (CONT.)
90
SLURRY CONCENTRATION CALCULATIONS (CONT.)
Sand per Foot: is the amount of sand per foot calculated on stage and well levels; measured in lb/ft
Water per Foot: is the amount of water used per foot calculated on stage and well levels; measured in bbl/ft
Sand to Water Ratio: is the total sand divided by the total water per stage; measured in ppg
91
SLURRY CONCENTRATION CALCULATIONS (CONT.)
Figure 45: Table of slick water injection schedule before calculations. [22]
92
SLURRY CONCENTRATION CALCULATIONS (CONT.)
Figure 46: Table of slick water injection schedule after calculations. [22]
93
REFERENCES
94
1. Miskimins, J. L. (2019). HYDRAULIC FRACTURING: Fundamentals and Advancements. Society of Petroleum
Engineers.
2. Speight, J. G. (2016). Handbook of Hydraulic Fracturing. Wiley.
3. Walter, G. R., & Thompson, G. M. (1982). A Repeated Pulse Technique for Determining the Hydraulic Properties of
Tight Formations. Ground Water, 20(2), 186–193. https://doi.org/10.1111/j.1745-6584.1982.tb02749.x
4. PetroWiki. (2016, March 22). Unconventional resources of oil and gas from a geologic perspective. PetroWiki.
https://petrowiki.spe.org/Unconventional_resources_of_oil_and_gas_from_a_geologic_perspective#Unconventional_
Resources_Classification_and_Distribution.
5. Donaldson, E. C., Alam, W., & Begum, N. (2014). Hydraulic Fracturing Explained: Evaluation, Implementation, and
Challenges (Gulf Drilling) (1st ed.). Gulf Publishing Company.
6. Wong, S. W. (2019). Hydraulic Fracture Modeling and Design - A Perspective on How things have Changed from
Conventional to Unconventional Reservoirs. Day 2 Wed, March 27, 2019. Published. https://doi.org/10.2523/iptc-
19405-ms
7. Maeso, C., Ponziani, M., Le Nir, I., Kherroubi, J., Quesada, D., Dubourg, I., Luthi, S., Slob, E., Fisher, K., Honeyman,
L., Brown, R., & Zenned, O. (2014). Fracture Aperture Calculations From Wireline and Logging While Drilling
Imaging Tools. All Days. Published. https://doi.org/10.2118/170848-ms
8. Kim, G. Y., & Narantsetseg, B. (2014). Fractures and Breakouts Analysis from Borehole Image Logs: Preliminary
Results for Interpretation of Stress History andIn SituStress State. Marine Georesources & Geotechnology, 33(6),
579–585. https://doi.org/10.1080/1064119x.2014.954179 95
9. Zhang, T., Li, Z., Adenutsi, C. D., & Lai, F. (2017). A new model for calculating permeability of natural fractures in dual-
porosity reservoir. Advances in Geo-Energy Research, 1(2), 86–92. https://doi.org/10.26804/ager.2017.02.03
10. Bahrami, H., Rezaee, R., & Hossain, M. (2012). Characterizing natural fractures productivity in tight gas reservoirs. Journal of
Petroleum Exploration and Production Technology, 2(2), 107–115. https://doi.org/10.1007/s13202-012-0026-x
11. Darvishpour, A., Seifabad, M. C., Wood, D. A., & Ghorbani, H. (2019). Wellbore stability analysis to determine the safe mud
weight window for sandstone layers. Petroleum Exploration and Development, 46(5), 1031-1038. doi:10.1016/s1876-
3804(19)60260-0
12. Rickman, R., Mullen, M. J., Petre, J. E., Grieser, W. V., & Kundert, D. (2008). A Practical Use of Shale Petrophysics for
Stimulation Design Optimization: All Shale Plays Are Not Clones of the Barnett Shale. All Days. Published.
https://doi.org/10.2118/115258-ms
13. Wikipedia contributors. (2020, August 10). Triaxial shear test. Wikipedia. https://en.wikipedia.org/wiki/Triaxial_shear_test
14. Donaldson, E. C., Alam, W., & Begum, N. (2013). Hydraulic Fracturing Explained: Evaluation, Implementation, and
Challenges (Gulf Drilling) (1st ed.). Gulf Publishing Company.
15. Aadnøy, B. S., & Looyeh, M. R. (2011). Petroleum rock mechanics: Drilling operations and well design. Oxford: Gulf
Professional.
16. Zhang, J. J. (2020). Applied Petroleum Geomechanics. Gulf Professional Publishing
17. Hossain, M., Rahman, M., & Rahman, S. (2000). Hydraulic fracture initiation and propagation: roles of wellbore trajectory,
perforation and stress regimes. Journal of Petroleum Science and Engineering, 27(3–4), 129–149.
https://doi.org/10.1016/s0920-4105(00)00056-5
96
18. Yang, L., Cao, J., Tao, Z., Wang, S., Yang, J., & Dong, B. (2020). Fracture Propagation Characteristics in
Horizontal Wells. Chemistry and Technology of Fuels and Oils, 56(3), 405–410. https://doi.org/10.1007/s10553-
020-01151-5
19. Economides, M. J., and T. Martin. 2007. Modern Fracturing - Enhancing Natural Gas Production. Energy
Tribune.
20. Wutherich, Kevin, Lili Xu, Kirby Walker, Walter Sawyer, and Isaac Aso. 2013. “Quantifying the Effect of Drilling
Azimuth in Shale Gas Reservoirs.” SPE Eastern Regional Meeting held in pittsburgh, Pensylvania, USAS, 20-22
August 2013. SPE 165672.
21. Wong, S. W. (2019). Hydraulic Fracture Modeling and Design - A Perspective on How things have Changed from
Conventional to Unconventional Reservoirs. Day 2 Wed, March 27, 2019. Published. https://doi.org/10.2523/iptc-
19405-ms
22. Guo, B., Liu, X., & Tan, X. (2017). Petroleum Production Engineering. Elsevier Gezondheidszorg.
23. Zhang, J. J. (2019). Applied Petroleum Geomechanics (1st ed.). Gulf Professional Publishing.
97