STABView Software
STABView Software
Adv
ew
ancedWellPl
anningSof
t
TM
war
e
t
oAnalyze:
Wel
lbor
eSt
abi
li
ty
SandPr
oduct
ion
Los
tCi
rcul
ati
on&Fr
act
uri
ng
Faul
torBeddi
ngPl
aneSl
ip
STABView™ Version 3.8
Wellbore Stability, Lost Circulation,
Fracturing, Sand Production and
Discontinuity Slip Analysis Software
USER’S MANUAL
November, 2008
Weatherford Advanced Geotechnology, Suite 1200, 333 – 5th Ave. S.W., Calgary, Alberta, CANADA T2P 3B6
Phone: 1(403)693-7530 Fax: 1(403)693-7541 E-Mail: software@advgeotech.com Web: www.advgeotech.com
STABView 3.8 User’s Manual ii November, 2008
TERMS OF USE
TABLE OF CONTENTS
TERMS OF USE ..................................................................................................................... ii
1. INTRODUCTION ........................................................................................................... 1
1.1 Background............................................................................................................. 1
1.2 Problems You Can Analyze with STABView ...................................................... 1
1.2.1 Drilling Applications.................................................................................... 1
1.2.2 Well Control Applications........................................................................... 1
1.2.3 Formation Evaluation Applications ........................................................... 1
1.2.4 Completions Applications............................................................................ 1
1.2.5 Stimulation Applications ............................................................................. 2
1.2.6 Production Engineering Applications ........................................................ 2
1.2.7 New and Novel Applications ....................................................................... 2
1.3 Copyright and Legal Information......................................................................... 2
4. THEORETICAL BACKGROUND............................................................................. 17
4.1 Rock Mechanics Principles and Applications.................................................... 17
4.1.1 Fundamentals ............................................................................................. 17
4.1.2 Stress, Strain and Sign Conventions ........................................................ 17
4.1.3 Rock Yielding Mechanisms ....................................................................... 18
4.1.4 Rock Strength and Mechanical Properties.............................................. 19
4.1.5 Apparent Capillary Cohesion ................................................................... 24
4.1.6 In-situ Stresses and Stress-Depletion Response ...................................... 27
4.2 Prediction of Rock Yielding and Failure Around Boreholes ........................... 31
4.2.1 Introduction................................................................................................ 31
4.2.2 2D Elastic Models....................................................................................... 32
4.2.3 3D Elastic Models....................................................................................... 36
4.2.4 Plane of Weakness Effects on Stability .................................................... 39
4.2.5 Thermal Effects on Stability ..................................................................... 40
4.2.6 Effect of Swab Pressures on Stability....................................................... 41
4.2.7 Calibration of Elastic Stability Models .................................................... 42
4.3 Elasto-plastic Models for Boreholes and Cylindrical Perforations ................. 43
4.3.1 Borehole and Perforation Yielded Zone Models..................................... 43
4.3.2 Elasto-plastic Models for Hemi-spherical Perforations ......................... 48
4.3.3 Kinematic Considerations ......................................................................... 49
STABView 3.8 User’s Manual v November, 2008
LIST OF APPENDICES
INDEX
STABView 3.8 User’s Manual vii November, 2008
TECHNICAL SUPPORT
Weatherford Advanced Geotechnology targets to provide a response to all enquires related to its
software products within 12 hours of the contact, during our normal business work week
(Monday to Friday, 8 am to 5 pm Mountain Standard Time). Normally an AG specialist can be
reached by telephone for rapid response during our normal work hours. On weekends and during
holidays the target response time is within 3 days.
Depending upon the nature of the software problem, our response will be email or telephone. In
some cases it may be possible to discuss the problem on the telephone using an internet
accessible version of the software using WebEx.
For software problems of a more serious nature, that may involve code changes or programming,
we attempt to service all clients with current support agreements as quickly as possible, by
sending a patch or software upgrade, if necessary. The typical response time for such
improvements is between 1 and 5 days, but may take more time depending on the complexity of
the issue and amount of work required.
STABView Technical Support
Weatherford Advanced Geotechnology
1200, 333 – 5th Ave. S.W.
Calgary, Alberta, CANADA T2P 3B6
Tel: 1-403-693-7530
Fax: 1-403-693-7541
Email: software@advgeotech.com
Web: www.advgeotech.com
STABView 3.8 User’s Manual viii November, 2008
STABView 3.8 User’s Manual 1 November, 2008
1. INTRODUCTION
1.1 Background
Welcome to the User’s Manual for STABView - a software program for wellbore stability, sand
production, lost circulation and discontinuity slip risk analysis. STABView was developed for
technical drilling, completion, production and exploration personnel to conduct rapid parametric
analyses for a wide range of near-wellbore integrity problems. Efficient, rapidly converging,
analytical and semi-analytical algorithms are used to solve the fundamental equations that
describe the stresses, pore pressures, and temperature distributions for near-wellbore problems.
Additional information regarding these models is included in Section 4 of this manual.
Information in this document is subject to change without notice. No part of this document or
STABView software package may be reproduced or transmitted in any form or by any means,
electronic or mechanical, for any purpose, other than as permitted in the license agreement,
without the expressed written permission of AG.
© Weatherford International Ltd., 2008. All rights reserved. Printed in Canada.
See the STABView License Agreement in Appendix A of this manual for further legal and
copyright details.
STABView 3.8 User’s Manual 4 November, 2008
STABView 3.8 User’s Manual 5 November, 2008
2. INSTALLING STABVIEW
2.1 Introduction
STABView requires a minimum of approximately 250 megabytes of hard drive space in order to
install the software and required components. Table 2.1 below outlines the minimum and optimal
system requirements in order to run STABView.
2.3 Installation
Follow the procedures listed below to install STABView. It is recommended that you close all
other applications before starting the installation procedure. For some systems, the installation
process may be adversely affected if anti-virus software is running in the background. In some
cases you may be required to restart your system after or during installation.
STABView 3.8 User’s Manual 6 November, 2008
Use the following instructions if you have been provided with a download link, username,
website password and installation password from Advanced Geotechnology.
1. Log on as the primary STABView user, this user should also have local Admin privileges
on the PC.
2. Remove any USB security keys used to run Advanced Geotechnology software from the
PC.
3. Click “Start” -> Control Panel -> Add/Remove Programs.
4. Click “Remove” (uninstall) all previous versions of STABView.
5. Close the Add/Remove Programs dialog.
6. Copy and paste into a web browser the download link provided to you by Advanced
Geotechnology. If Advanced Geotechnology has not provided you with a download link
then you will not be able to install the software in this manner.
7. When prompted for credentials from the website enter the username and password also
provided by Advanced Geotechnology.
8. You will be prompted to Open or Save the file. Save the setup.exe file to a well known
location on your hard drive.
9. After the file has completed downloading double click to launch the setup.exe file.
10. Enter the installation password, again this should have been provided to you by
Advanced Geotechnology.
11. Follow the instructions provided by the setup utility to install STABView.
12. If you do encounter an error during installation (common on Windows 2000) just click
OK to continue.
13. Copy and paste into a web browser or click on the following link to download the version
of DirectX from Microsoft.
http://advgeotech.com/STABViewPrivate/directx_nov2007_redist.exe
(If you are prompted to login use the same credentials as above.)
14. You will be prompted to Open or Save the file. Save the directx_nov2007_redist.exe file
to a well known location on your hard drive.
15. After the file has completed downloading double click to launch the file and follow the
instructions provided by the setup utility to install DirectX.
16. Execute the newly installed version of STABView 3.8 from the shortcut placed on your
desktop or through the “Start” menu.
17. STABView has been developed to run under Windows with the English language. You
may encounter problems if your system is running a different language or has Unicode
support for languages other than English.
18. STABView will prompt you for a new key code. Please copy and send the serial number
and lock code from the registration dialog to software@advgeotech.com or 1(403) 693-
7530. Software support will then issue you a new key code. A single use license permits
you to install STABView on up to three PCs, but you will still be required to contact us
for a new key code for each installation.
Enter the key code we provide you and click “Unlock”. You will then be prompted by an
activation dialog to make sure your PC is connected to the internet. This is required to authorize
and confirm this installation of STABView. Make sure your PC is connected and click
STABView 3.8 User’s Manual 7 November, 2008
“Activate”. If you are unable to activate the software you will prompted with a dialog requesting
additional information. Contact software support for further instructions.
Use the following instructions if you have been provided with an installation CD from Advanced
Geotechnology.
1. Log on as the primary STABView user, this user should also have local Admin privileges
on the PC.
2. Remove any USB security keys used to run Advanced Geotechnology software from the
PC.
3. Click “Start” -> Control Panel -> Add/Remove Programs.
4. Click “Remove” (uninstall) all previous versions of STABView.
5. Close the Add/Remove Programs dialog.
6. Insert the STABView installation CD into a drive.
7. This will launch the STABView installation and you will next be prompted with a
password dialog. If the application does not launch right away then use Windows
Explorer to open the CD drive and execute setup.exe. Enter your assigned installation
password that should have been included with the CD. If you are unable to locate the
install password contact software support.
8. Follow the instructions provided by the setup utility to install STABView.
9. If you do encounter an error during installation (common on Windows 2000) just click
OK to continue.
10. When the installation has completed insert the USB security key and launch STABView
11. Execute the newly installed version of STABView 3.8 from the shortcut placed on your
desktop or through the “Start” menu.
12. STABView has been developed to run under Windows with the English language. You
may encounter problems if your system is running a different language or has Unicode
support for languages other than English.
13. STABView will prompt you for a new key code. Please copy and send the serial number
and lock code from the registration dialog to software@advgeotech.com or 1(403) 693-
7530. Software support will then issue you a new key code. A single use license permits
you to install STABView on up to three PCs, but you will still be required to contact us
for a new key code for each installation.
14. Enter the key code we provide you and click “Unlock”. You will then be prompted by an
activation dialog to make sure your PC is connected to the internet. This is required to
authorize and confirm this installation of STABView. Make sure your PC is connected
and click “Activate”. If you are unable to activate the software you will prompted with a
dialog requesting additional information. Contact software support for further
instructions.
Assistance with the installation of a licensed copy of STABView on a local area network can be
provided to licensees as part of their purchase agreement. Please, contact STABView Technical
Support at 1(403)693-7530 or e-mail software@advgeotech.com if you require information or
assistance.
After successful installation of the software, a new folder with a default name of "STABView
3.8" will be created in the "Program" folder accessed from the Windows "Start" button. This
folder will contain the following shortcuts:
1. STABView 3.8 Manual (PDF)
2. STABView 3.8
3. Uninstall STABView 3.8
Click on "STABView 3.8" to run the program, or select "STABView 3.8 Manual" to access the
STABView documentation. Note that the help file is also accessible from the "Help" menu once
the STABView 3.8 program has been activated.
If your installation of STABView came shipped with a USB hardware key, please make sure to
always have this key inserted into an available USB port when running STABView. When
inserting the USB key please wait a few seconds for your PC to recognize the new hardware
before launching STABView. In the event your USB key becomes lost or damaged please
contact STABView Technical Support at 1(403) 693-7530 or e-mail: software@advgeotech.com
to order a replacement key. A replacement fee will be charged in order to replace any lost or
stolen keys.
If your installation of STABView came shipped with a network USB hardware key, please see
the documentation that accompanied the network installation package. In most cases you will
need to enlist the aid of your network admin or IT support staff in order to use the software in
this manner.
STABView 3.8 User’s Manual 9 November, 2008
A software security interface will appear when STABView is initially run. The software will not
run until it has been unlocked by entering the required "key" code. This key code can be
obtained by registered users of the software only, by contacting STABView Technical Support.
The STABView software key code is only valid for one personal computer. Additional codes are
required for each computer that STABView is installed on. Single user licensees for STABView
may install the software on up to three computers although it can only be run on one computer at
any time with its USB hardware key. Installations of the software on more than three PCs can be
arranged for an additional fee.
2.5.4 Activation
Some installations of STABView will require that you activate the software over the internet. If
you are prompted to activate the software after unlocking, then please make sure that your PC is
connected to the internet and then click “Activate”. If for some reason you are unable to activate
your installation of STABView please contact your network administrator first to obtain access
to the internet. If there are security or access problems then you should contact STABView
Technical Support at 1.403.693.7541 or software@advgeotech.com.
STABView 3.8 User’s Manual 10 November, 2008
STABView 3.8 User’s Manual 11 November, 2008
3. STABVIEW FEATURES
Note: Technical features marked with an asterisk (*) are currently in development. Contact
Advanced Geotechnology if you require further details.
• Evaluate instability risks as function of ECD, EMW, BHP, SG, mud density, over-
or underbalance pressures
• 3D linear elastic models
• 2D elastoplastic models
• 3D passive shear failure initiation (borehole ballooning)
• Swab and surge pressure effects
• Underbalanced drilling features
• Calculate a mud density, EMW, ECD or bottomhole pressure required to achieve a
tolerable instability risk, e.g., hole enlargement
• Calculate wellbore instability risks for a specified mud density, ECD, EMW or
bottomhole pressure
• Calculate and display the drilling “mud weight window”
• Thermo-elastic effects due to steady-state conductive heat transfer
• Predict a profile of enlarged hole sizes for hydraulics optimization
• 2D elastoplastic models
• 3D linear elastic models
• Extent of horizontal well collapse (rubble fill percentage)
• Near-wellbore skin can be used as an input parameter
• Perforated, openhole, slotted liner, screen and expandable completions can be analyzed
• Cylindrical or hemi-spherical perforation cavities
• Input pressure drop across liners, screens and expandables
• First-order estimate of isotropic loading on liners, screens or expandables
• First-order estimate of borehole wall deformations
• Account for the effects of ubiquitous, weak discontinuities on the risk of borehole
collapse
• Mohr-Coulomb failure criterion
• Shear failure and slip tendency analysis for weak bedding planes, faults, natural fractures,
cleats
• Re-opening pressure analysis for weak bedding planes, faults, natural fractures, cleats
* Display log-derived bedding planes faults, natural fractures and cleats on polar plots*
• Steady-state flow with pore pressure gradients into or out of the wellbore
• Capillary threshold pressure for OBMs and pseudo-OBMs
• Apparent capillary strength for weak, partially saturated sands
• Filter cake efficiency model for permeable sandstones
• Wall coating efficiency model for shales
• Effect of an instantaneous BHP change
• Fluid viscosity and permeability effects
• Formation damage and skin effects
• Steady-state non-Darcy flow effects for high rate gas wells and perforations
• Compressible fluid effects in the near-wellbore area
• Skin damage option for underbalanced drilling
• Biot or Terzaghi effective stress
STABView 3.8 User’s Manual 13 November, 2008
2D Elastoplastic strain-weakening
• Mohr-Coulomb Model
3D Linear Elastic Models
• Mohr-Coulomb
• Modified Lade
* Mogi-Coulomb
• Non-linear Hoek-Brown
• Tensile fracture criteria
• Passive shear failure initiation
Other Options
∗ Time-dependent loss or gain of rock strength and elastic properties
• Osmotic pressure model for reactive clay inhibition effects based on shale and mud
activities (API RP13B specification)
• Handles many common oil-based and water-based drilling fluids
• Database of published shale water activities and membrane efficiencies for a variety of
shales and fluids
• Drilling fluid activity calculator
STABView 3.8 User’s Manual 14 November, 2008
• Fix an acceptable BHP, drawdown pressure or mud density based on the performance of
an offset well
• User-defined borehole breakout angle criterion
• Calibration options for borehole collapse and fracture breakdown
STABView 3.8 User’s Manual 15 November, 2008
3.17 Units
3.18 Functionality
4. THEORETICAL BACKGROUND
4.1 Rock Mechanics Principles and Applications
4.1.1 Fundamentals
Rock mechanical analysis in petroleum engineering involves the assessment of stresses induced
in rock, due to well drilling, completion, fluid injection, temperature changes and the associated
deformations and yielding that may occur. The application of rock mechanics to engineering
design involves the evaluation of these stresses, deformations and yielding events against
established performance criteria. If the predicted rock response fails to meet these criteria, a
more conservative design is often desirable or required. This section provides brief explanations
of the fundamental terms and principles that underlie most rock mechanical models used in
STABView. Readers who are unfamiliar with rock or soil mechanics should consult textbooks
on either of these subjects for more detailed background information, e.g., Fjaer et al, 1992.
Stress is a measure of the internal distribution of force per unit area within a body that balances
and reacts to the loads applied to it. Stress is a tensor quantity with nine terms, but can be
described fully by six terms due to symmetry: 3 normal stress components (denoted σii, or
sometimes simply σi) and 3 shear stress components (denoted τij). Stress can occur in liquids,
gases, and solids. Liquids and gases support normal stress (pressure), but flow under shear stress.
Solids support both shear and normal stress.
At a given point in the subsurface it is possible to locate three orthogonal planes on which the
shear stress vanishes. These planes are called the principal planes, while the normal stresses on
these planes are the principal stresses.
As illustrated in Figure 4.1, a portion of the total mechanical stress applied to a saturated porous
medium is carried by the pore fluid, while the remainder of the stress is carried by the skeleton of
mineral grains. The latter stress is referred to as the effective stress. Effective stress is defined as
follows:
σ ii' = σ ii − α p (4.1.1)
where:
σ'ii = effective normal stress
σii = total normal stress
α = Biot’s coefficient
p = pore pressure
The value of Biot’s coefficient is often assumed to be very close to 1, which is true for most
weak, porous rocks. Because fluids cannot sustain shear stresses, shear stress components acting
on mineral grains within a porous rock are unaffected by fluid pressure.
For porous media containing pressurized pore fluids, it is the effective stresses that govern most
rock deformation, yielding, and failure conditions.
STABView 3.8 User’s Manual 18 November, 2008
σii
Figure 4.1: Effective stress concept: fluid pressure within a porous media acts against an
externally applied normal stress, σii, reducing the stress that acts at grain contacts.
Strain is the geometrical expression of deformation caused by the action of stress on a physical
body. Strain is a tensor quantity with nine terms, but can be described fully by six terms due to
symmetry: 3 compressive strain components (denoted εii, or sometimes simply εi) which
represent the change in length of a line, and 3 shear strain components (denoted εij) which
represent the change in angle between two lines.
The geomechanics sign convention, which will be used throughout this manual, is that
compressive stresses are positive and contractile strains are positive.
Rock around a borehole or perforation can yield or hydraulically fracture in different ways
resulting in wellbore instability, or total or partial collapse sand production, depending upon the
type of rock, its properties and other factors. See Figure 4.2 for a summary of mechanisms that
can result in sand production. When effective normal stresses become negative, they are said to
be tensile. Tensile stresses act to force a solid body apart. Tensile stresses can develop in
producing wells under drawdown conditions when a steep pressure gradient exists near the
borehole or perforation wall. If they exceed the tensile strength of the rock (which is often very
small), these tensile stresses will result in tensile yielding and concentric shells of rock will
detach from the borehole or perforation wall. Sand production or spalling resulting from tensile
yielding tends to be transient and short-lived in nature, occurring most often immediately after an
increase in drawdown pressure.
In a compressive stress regime, if the largest principal stress component (denoted σ1) is
significantly greater than the smallest principal stress (denoted σ3), large shear stresses will exist
on planes oriented obliquely to these principal stresses. If the shear stresses resulting from this
anisotropic compression exceed the rock’s shear strength, compressive shear yielding will occur.
The type of stress regime required to induce compressive shear yielding exists in most near-well
STABView 3.8 User’s Manual 19 November, 2008
Compressive
τ Shear Failure
Tensile Erosion
Failure
Drawdown
σr ’ < 0 Q
Induced
Transient flow induced
1991)
environments; i.e., the radial stress acting normal to the borehole or perforation surface (σr) tends
to be small, while the “hoop or tangential stress acting tangent to the borehole or perforation
surface (σθ) tends to be large. Excessive amounts of compressive shear yielding can result in
severe borehole instability sand or solids productions problems, and potentially catastrophic
collapse of boreholes and perforations.
At high flow rates (e.g., high-rate gas wells), the flowing fluid may induce a sand erosion
mechanism, in which individual grains detach from the rock. In some cases, or if this type of
failure persists for long periods of time, the resultant sand production may be problematic.
Erosion of sand grains is often more of a concern in cases where compressive shear yielding has
occurred, in which case the rock’s susceptibility to erosion is enhanced. As such, the yielding
mechanism that poses the greatest risk for sand production is compressive shear yielding.
Consequently, this yielding mechanism will be emphasized in this section and subsequent parts
of this manual. Sand erosion models are not currently part of STABView.
4.1.4 Rock Strength and Mechanical Properties
Figure 4.3 shows typical axial and radial stress-strain responses measured on a cylindrical rock
sample during a triaxial compression test. Before reaching the peak stress level, the rock behaves
more-or-less as a linear elastic solid. This behavior is described by Young’s modulus (which is a
measure of the axial stiffness of the rock) and Poisson’s ratio (which is a measure of the lateral
expansion of the rock in response to an axial load). Definitions for these parameters are given in
Figure 4.4 and procedures for interpreting them from a triaxial compression test are shown in
Figure 4.5.
STABView 3.8 User’s Manual 20 November, 2008
Peak strength
σa
ε a = ( L − L*) / L
(positive)
ε r = (r − r*) / r
(negative)
Poisson’s ratio:
L L**
−εr
υ=
εa
Young’s’Modulus
r
σa
E=
r* εa
Figure 4.4: Elastic properties of rock: Young’s modulus and Poisson’s ratio
εa denotes axial strain, and εr denotes radial strain.
STABView 3.8 User’s Manual 21 November, 2008
120
Peak Strength, σp
100
Deviatoric Stress, MPa
80
20
0
-0.2 -0.1 0 0.1 0.2 0.3
Figure 4.5: Basic rock elastic properties determined from triaxial compression tests.
8 peak strength 3 15
5 20
Shear Stress (MPa)
7 c = 0.88 MPa
6 phi = 31°
5
4
3
2
1
0
0 5 10 15 20
Normal
Effective Stress
Normal (MPa)(MPa)
Stress
Figure 4.6: Mohr circles representing the stress states at which compressive shear yielding
occurred in triaxial compression tests conducted at four different confining stresses.
STABView 3.8 User’s Manual 22 November, 2008
The peak stress level shown in Figure 4.5 occurs when compressive shear yield occurs. For
frictional materials such as rock and soil, this peak stress level tends to increase with increasing
confining stress, as illustrated by the data plotted in Figure 4.6. The shear yield criterion for a
rock is a mathematical function that defines the triaxial stress state at which compressive shear
yielding occurs.
The most commonly used yield criterion for rocks is the linear Mohr-Coulomb criterion, which
is illustrated in Figure 4.7. This criterion is expressed by the following equation:
where:
τp = shear stress at compressive yield (peak shear stress)
c = cohesion
σ’n = effective normal stress
φ = friction angle
Another commonly used yield condition in rock mechanics is the Hoek and Brown criterion
(Hoek, 2000). This non-linear criterion is illustrated in Figure 4.8, and is defined as follows:
a
⎛ σ' ⎞
σ = σ + UCS ⎜⎜ m 3 + s ⎟⎟
'
1
'
3 (4.1.3)
⎝ UCS ⎠
where:
UCS = unconfined compressive strength
m = frictional strength parameter
s = 1 for intact rock; < 1 for broken rock masses
a = 0.5 for intact rock; < 0.5 for broken rock masses
Several other failure criteria have been used in rock mechanics, some of which account for the
effect of the intermediate principal stress (σ2). An example of such a criterion is the modified
Lade criterion, which has been used for sand production prediction models by Ewy (2001).
According to this criterion, yield occurs when the following conditions are met:
where:
( ) ( ) (
I 1 ' ' = σ 1' + S1 + σ 2' + S1 + σ 3' + S1 ) (4.1.5)
STABView 3.8 User’s Manual 23 November, 2008
10
9 τ = c + σ n' tan φ
Slope = tan(φ)
8
Shear Stress (MPa)
7
6
eld
5
4 Yi
d
3 y iel
2 No
1
0
0 5 10 15 20
Cohesion (c) Normal
Effective Stress
Normal (MPa)
Stress (MPa)
Figure 4.7: Linear Mohr-Coulomb yield criterion for the data plotted in Figure 4.6.
25 0.5
⎛ σ3' ⎞
σ 1 ' = σ 3 '+2.79⎜14.6 + 1⎟
xyz sandstone ⎝ 2.79 ⎠
Maximum Effective Stress (MPa)
20 peak strength
15
0
0 1 2 3 4 5 6
Minimum Effective Stress (MPa)
Figure 4.8: Hoek and Brown yield criterion for the data plotted in Figure 4.6.
STABView 3.8 User’s Manual 24 November, 2008
( )( )(
I 3 = σ 1' + S1 σ 2' + S1 σ 3' + S1 ) (4.1.6)
Only two rock strength parameters are required for this criterion, which are analogous (and direct
functions of) cohesion and friction angle, as follows:
c
S1 = (4.1.7)
tan φ
4 tan 2 φ (9 − 7 sin φ )
η= (4.1.8)
(1 − sin φ )
As the stress state in a rock approaches its peak level, microcracks begin to develop. At peak
strength, compressive shear yielding occurs as a consequence of these microcracks coalescing to
form a macroscopic shearing surface. Several types of post-peak behavior may occur, depending
on factors such as rock lithology, mean effective stress and temperature. Figure 4.9 illustrates
three idealized post-peak responses, where a more-or-less constant residual stress level is often
reached, due largely to the frictional resistance to sliding along the newly developed shear
surface or surfaces. This residual strength can play an important role in the stabilization of
yielded rocks. Like peak strength, residual strength can be characterized using various yield
criteria. An example of linear Mohr-Coulomb yield criteria for peak and residual strength in a
sandstone is shown in Figure 4.10.
4.1.5 Apparent Capillary Cohesion
Capillary pressures exist within porous rocks containing two or more immiscible fluids. In the
case of a water-wet, primarily gas-saturated, sandstone reservoir, for example, these pressures
exist because the water phase tends to adhere to the mineral grains. These capillary pressures act
to “pull” the mineral grains together, contributing an additional component of strength to the
rock. Although the magnitude of this apparent capillary cohesion will be negligible for many
rocks, it may represent a significant proportion of the cohesive strength of weak, poorly-
cemented, fine grained sandstones and siltstones.
It is reasonable to expect that the magnitude of apparent capillary cohesion should change as
fluid saturations change. Consequently, apparent capillary cohesion has been identified as one
possible mechanism to account for the significant increase in sand production commonly
experienced with increasing water cut. Bruno et al. (1998) demonstrated this effect with
laboratory experiments, and with numerical modeling of a partially-saturated assemblage of disk-
shaped particles using the discrete element method. In spite of the apparent importance of this
effect, practical models are not usually applied for predicting the magnitude and saturation-
sensitivity of apparent capillary cohesion for sandstones.
STABView 3.8 User’s Manual 25 November, 2008
Residual strength
Residual
strength
Stress
Stress
Stress
-plastic
Residual strength
10
9 τ p = c p + σ n' tan φ p
Slope = tan(φp)
8
Shear Stress (MPa)
7
6
5
Slope = tan(φr)
4
3
2 τ r = cr + σ n' tan φr
cp 1
cr 0
0 5 10 15 20
Normal
Effective Stress
Normal (M Pa)
Stress (MPa)
Figure 4.10: Characterization of peak and residual strength for rock using linear Mohr-
Coulomb yield criteria. Subscript p denotes peak, and r denotes residual.
STABView 3.8 User’s Manual 26 November, 2008
In the discipline of unsaturated soil mechanics, the effects of water-air saturation on the shear
strength of soil has been studied extensively (Vanapalli et al., 1996). Borrowing from this body
of work, and converting to parameters used in petroleum engineering, the following expression is
proposed for predicting apparent capillary cohesion (cc) as a function of water saturation (Sw):
S w − S wi
cc = pc tan φ (4.1.9)
1 − S gi − S wi
where:
Swi = irreducible water saturation (fraction)
Sgi = irreducible gas saturation (fraction)
pc = capillary pressure (see Section 4.6.2)
Figure 4.11 shows a sample calculation of apparent capillary cohesion using equation 4.1.9,
using a weak sand friction angle of 40°. Apparent capillary cohesion is observed to increase
steadily as water saturation begins to drop below its maximum value, reach a relatively constant
value of roughly 3 psi (20 kPa) over the intermediate-low range of Sw, then decreases abruptly as
Sw approaches its minimum value. Given that air-water and methane-water interfacial tensions
are approximately equal, air-water capillary pressure data of the type shown in this example
could be used to calculate apparent capillary cohesion in gas reservoirs.
300
Capillary Pressure, Capillary Cohesion (kPa)
200
150
100
Swi 1-Sgi
50
0
0 0.2 0.4 0.6 0.8 1
(a) (b) Water Saturation
Figure 4.11: (a) Distribution of water and gas in a sandstone reservoir, at low water
saturation. (b) Apparent capillary cohesion calculated from air-water capillary pressure
data for a poorly-sorted, medium fine-grained sandstone.
STABView 3.8 User’s Manual 27 November, 2008
According to the Laplace equation (Bear, 1972) as described in Section 4.6.2 of this manual, the
magnitude of the capillary pressure is inversely proportional to the square of pore throat radius.
As such, the magnitude of this effect is greatest for fine-grained rocks. Further, it is important to
note that these capillary effects are subject to change if the mineral grain structure of a rock is
compressed (e.g., due to elastic compression or consolidation induced by stress or pressure
changes) or loosened (e.g., due to yielding), hence changing pore throat sizes.
4.1.6 In-situ Stresses and Stress-Depletion Response
The native stresses that exist in a reservoir prior to hydrocarbon extraction are referred to as in-
situ stresses. Characterization of the complete in-situ stress tensor requires the determination of
the magnitudes and orientations of three principal stress components, denoted from largest to
smallest, as σ1, σ2 and σ3. In settings with gentle surface relief, a common assumption is that the
ground surface is flat and horizontal. In such a case, one of the principal stress orientations is
approximately vertical. This vertical in-situ stress, denoted by σV, results from the weight of the
overlying rocks, sediments and the fluids contained within them. Due to the orthogonal nature of
principal stresses, the other two principal stresses will be in the horizontal plane. These principal
stresses are denoted by σHmin and σHmax for the minimum and maximum horizontal in-situ
stresses, respectively (see Figure 4.12). Commonly, especially in extensional sedimentary
basins, σV is larger than both σHmin and σHmax. In convergent or compressive basins, however,
lateral compressive forces may result in stress regimes wherein σV is either the intermediate or
the least principal stress.
fa ce
Sur
n d
r ou
G
σ Hmax
σV
σHmin
The pore pressure changes associated with hydrocarbon production or injection alter the stress
state in and around reservoirs. If a reservoir was a free body, production- or injection-induced
effective stress changes would simply result in reservoir contraction or expansion, without any
change in total in-situ stresses. However, the reservoir is “attached” to the surrounding rock,
within which there is no driving force for displacement because pore pressures are not changing.
Due to this competition between internal driving forces and external constraints, anisotropic
changes in total stress may be induced, as illustrated in Figure 4.13. The relationship between
pore pressure changes and total stress magnitudes is often referred to as the reservoir’s stress-
depletion response. It is important to understand the stress-depletion response of a reservoir for
several practical petroleum engineering problems. For example, a reservoir may initially be
produced sand-free, but at a later time, depending on the magnitude of the depletion-induced
stress changes, it can cross a threshold and begin producing sand.
Normalized stress arching ratios, defined as follow, are useful parameters for characterizing the
stress-depletion response of a reservoir:
∆σ H ∆σ V
γH = , γV = (4.1.10, 4.1.11)
∆p ∆p
Ground
surface
Stress change in
Reservoir rocks surrounding
the reservoir
Stress change
within the
reservoir
where:
γH = horizontal normalized stress arching ratio
γV = vertical normalized stress arching ratio
∆σH = induced change in horizontal stress magnitude
∆σV = induced change in vertical stress magnitude
∆p = reservoir (pore) pressure change (negative for depletion)
For the special case of a laterally infinite reservoir, poroelastic theory predicts the following
expressions for stress arching ratios:
1 − 2ν
γH =α , γV = 0 (4.1.12, 4.1.13)
1 −ν
Adopting the usual assumption that Biot’s coefficient (α) is close to 1 for porous reservoir rocks,
it is evident that production- or injection-induced horizontal stress changes are controlled by
Poisson’s ratio (υ). Horizontal stress arching ratios predicted with equation 4.1.11 have proven
to be consistent with field observations in several sedimentary basins. Specifically, γH tends to
fall in the range from 0.5 to 0.8 (Addis et al., 1996). In published accounts of stress-depletion
response, it has commonly been assumed that the vertical stress arching ratio is zero. Although, it
is reasonable to expect this to hold true in reservoirs with large lateral dimensions, it is not
expected to be true for reservoirs with thicknesses that are comparable to their lateral
dimensions.
Using poroelastic theory, Soltanzadeh and Hawkes (2007) have shown that the normalized stress
arching ratios for a plane strain, elliptical reservoir of width (w) and thickness (h) are defined as
follows:
1 − 2ν 1 1 − 2ν e
γH =α , γV = α (4.1.14, 4.1.15)
1 −ν 1 + e 1 −ν 1 + e
where:
e = reservoir aspect ratio
=h/w
It is also reasonable to expect that contrasts in the mechanical properties of the reservoir and its
surrounding rocks should have an effect on reservoir stress-depletion response. Soltanzadeh et al.
(2007) have recently derived the following equations for stress arching ratios in such conditions,
for a plane strain reservoir that is elliptical in cross-section:
C1 C2
γH = , γV = (4.1.16, 4.1.17)
C3 C3
where:
C1 = α (1 − 2υ * )[ RG [e(1 − 2υ ) + 2(1 − υ )] + e] (4.1.18)
STABView 3.8 User’s Manual 30 November, 2008
Sensitivity analyses conducted with equation 4.1.16 and 4.1.17 have lead to the following
observations regarding stress changes within a reservoir:
• Induced stress changes are sensitive to the reservoir aspect ratio (height/width), but not its
absolute size;
• Induced vertical stress changes increase as the reservoir becomes more equi-dimensional
in cross-section;
• Induced horizontal stress changes decrease as the reservoir becomes more equi-
dimensional in cross-section, except for cases where the shear modulus of the reservoir is
less than the “caprock” (i.e., surrounding rock);
• Vertical and horizontal induced stress changes both decrease as the reservoir:caprock
shear modulus ratio increases; and
• Vertical and horizontal induced stress changes increase as the Poisson’s ratio of the
reservoir decreases
• Poisson’s ratio of the caprock has limited effect on induced stresses.
If pressure depletion occurs in a reservoir located within such a setting, poroelastic effects would
tend to decrease the magnitude of the minimum horizontal stress. However, a limit is placed on
the magnitude of this reduction, as the yield criterion in equation 4.20 must be satisfied. In such a
case, assuming that fault plane cohesion is negligible, the reservoir stress-depletion response is
given by:
2α sinφ
γH = (4.1.22)
(1 + sinφ )
STABView 3.8 User’s Manual 31 November, 2008
STABView currently gives the user the choice of using either a “poroelastic” horizontal stress
correction governed by equations 4.1.10 and 4.1.11, or a passive basin, limit equilibrium-based
horizontal stress correction governed by equation 4.1.21. The dialog for “stress depletion”,
which is accessed under the “Advanced Features” tab in “Model Properties”, is shown below in
Figure 4.14. New, general purpose, poro-elastic stress corrections for depletion or injection
processes, based on Soltanzadeh and Hawkes (2007) are in development and will be
implemented in a future version of STABView.
4.2.1 Introduction
In this section, a number of models for predicting the onset and extent of compressive shear
yielding will be discussed. Many of these models were originally developed to predict yielding
around boreholes. However, none of these models is scale-dependent, so they can all be applied
to any cylindrical cavity in the subsurface. As such, these models can also be used for
perforations that are long and thin, as shown in Figure 4.15. The only differences between a
borehole and a perforation are the radius and the orientation, e.g., perforations are horizontal in a
vertical well. The one exception to be presented here pertains to a sand production model
available in STABView, for hemi-spherical cavities, which can be used for perforations that are
short and fat.
Figure 4.14: STABView dialog showing options available for poro-elastic and limit
equilibrium-based stress corrections to account for depletion effects.
STABView 3.8 User’s Manual 32 November, 2008
casing
ro
Fluid
cement flow Lp Lp dp
Lp = Perforation length h
dw
φ
D Lp= Perforation diameter
p
p
Dwd=w Wellbore diameter
ro
R0φ = Outer radius
φ = Phase angle
In-situ stresses act within all sub-surface rock formations and it is usually assumed that they do
not exceed the rock strength. However, the act of creating a borehole or perforation induces a
change in the stress state within the surrounding rock. The focus of the following section will be
the analysis of these induced stresses, and whether or not they exceed the rock strength near a
borehole or perforation.
The simplest type of rock mechanical model for assessing the stability of rock around a borehole
or cylindrical perforation is derived assuming that the formation is a homogeneous, isotropic,
linear-elastic continuum. For cases where the borehole or perforation axis is parallel to one of the
principal in-situ stresses, the effective stresses at any point on the borehole or perforation wall
can be found as follows:
σ r' = p w − αp a (4.2.1)
τ rθ = τ rz = τ θz = 0 (4.2.4)
where:
σ r' = effective radial stress
STABView 3.8 User’s Manual 33 November, 2008
=α
(1 − 2υ )
(1 − υ )
υ = Poisson’s ratio
θ = angle (measured counter-clockwise) from the σMax direction to the point of
interest on the borehole (or perforation) wall
An example of the radial, tangential and axial stresses around a borehole is given in Figure 4.16
for the following input parameters, which are representative of an overbalanced drilling case:
would correspond to the minimum horizontal stress. Alternatively, to analyze stresses around a
cylindrical perforation tunnel oriented parallel to the maximum horizontal stress, σMax would
correspond to the greater of the vertical stress and the minimum horizontal stress, and σMin would
correspond to the lesser of these two in-situ stresses.
Figure 4.17 graphically depicts the orientations of radial, tangential and axial stress components
at a point on the borehole wall, and illustrates a typical profile of these stresses with increasing
radial distance.
Compressive shear yielding is most likely to occur on the borehole or perforation wall at a point
where the maximum principal stress (which is commonly σθ) is large and the minimum principal
stress (which is commonly σr) is small. As demonstrated by equations 4.2.1 through 4.2.4, and
illustrated in Figure 4.16, this will occur at points that are rotated 90° and 270° from σMax. At
such a point, the effective stresses can be found as follows:
σ r' = p w − αp a (4.2.5)
7000
σHmin radial stress
tangential stress
6000
axial stress
σHmax θ
Effective Stress (psi)
5000
4000
3000
2000
1000
0
0 30 60 90 120 150 180 210 240 270 300 330 360
Angle from σHmax (°)
po
where:
UCS = unconfined compressive strength
φp = peak friction angle
Np = Mohr-Coulomb peak strength parameter
1 + sin φ p
=
1 − sin φ p
Substituting equations 4.2.5 and 4.2.6 into equation 4.2.8, the following equation for the
bottomhole pressure at which compressive shear yielding initiates is obtained:
Different yield criteria, such as Hoek and Brown (Hoek, 2000) or modified Lade (Ewy, 2001),
can also be used to assess this critical value of bottomhole pressure. However, the basic
principle remains the same. The first step is to calculate the induced stresses, and the second step
is to assess how theses stresses compare to the rock’s strength.
The equations presented in the preceding section are valid for boreholes or perforations that are
oriented parallel to one of the principal in-situ stress components. In cases when the borehole is
not in a principal stress direction (or within a tolerable range), it is still possible to evaluate the
stresses at the borehole or perforation wall analytically, but it is not possible to derive closed
formed solutions for the critical bottomhole pressure at which shear yielding initiates.
A useful pair of coordinate systems for inclined well problems is illustrated in Figure 4.18. A
Cartesian coordinate system that is aligned with the principal in-situ stress directions is used as a
frame of reference. The azimuth (aw) of the borehole or perforation is measured counter-
clockwise from the σHmax direction. The inclination (iw) of the axis is measured with respect to
vertical. All calculations of induced stresses must be conducted in a coordinate system that is
aligned with the borehole or perforation axis. In Figure 4.18, such a coordinate system has been
drawn with the z-axis parallel to the borehole or perforation axis, the x-axis oriented parallel to a
radial line passing through the highest point on the borehole or perforation wall, and the y-axis
normal to the plane define by the x- and z-axes. For such a coordinate system, the in-situ stress
tensor can be evaluated as follows:
σx
σHmax
aw
σy
σHmin
iw
σz
σV
Figure 4.18: Coordinate systems used for in-situ stress tensor transformations
in inclined well problems.
STABView 3.8 User’s Manual 37 November, 2008
τxy = Lxx Lyx σHmax + Lxy Lyy σHmin + Lxz Lyz σV (4.2.13)
τyz = Lyx Lzx σHmax + Lyy Lzy σHmin + Lyz Lzz σV (4.2.14)
τxz = Lzx Lxx σHmax + Lzy Lxy σHmin + Lzz Lxz σV (4.2.15)
where:
Lxx = cos(aw) cos(iw) ; Lyx = -sin(aw) ; Lzx = cos(aw) sin(iw)
Subsequently, the induced stresses at a given point on the borehole or perforation wall can be
calculated as follows (after Aadnoy, 1987):
σ r = pw (4.2.16)
For this inclined well geometry, θ represents the angle counter-clockwise from the x-axis.
Yield criteria are often most conveniently expressed in terms of principal stresses. As such, for a
given bottomhole pressure, in order to evaluate if compressive shear yield occurs at a selected
point on the borehole or perforation wall, it is necessary to determine the principal stress
components at that point. Given that the borehole or perforation wall is bounded by fluid, the
radial stress is, by definition, a principal stress (see Figure 4.19). In the plane normal to the radial
stress (i.e., tangent to the borehole surface), the maximum and minimum principal stresses (σtmax
and σtmin, respectively) can be found as follows:
1
σθ + σ z ⎡ 2 ⎛ σ θ − σ z ⎞2 ⎤ 2
σ t max = + ⎢τ θz + ⎜ ⎟ ⎥ (4.2.20)
2 ⎢⎣ ⎝ 2 ⎠ ⎥⎦
STABView 3.8 User’s Manual 38 November, 2008
1
σθ + σ z ⎡ 2 ⎛ σ θ − σ z ⎞2 ⎤ 2
σ t min = − ⎢τ θz + ⎜ ⎟ ⎥ (4.2.21)
2 ⎢⎣ ⎝ 2 ⎠ ⎥⎦
Before using these stresses to evaluate the yield condition, it is necessary to determine the
relative values of these principal stresses; i.e.,
Thus, for a given bottomhole pressure, yielding at a point on a borehole or perforation wall is
predicted if:
(σ 1 − αp a ) ≥ UCS + N p (σ 3 − αp a ) (4.2.24)
Principal stresses, as determined using equations 4.2.22 and 4.2.23, can similarly be used in other
yield criteria such as Hoek-Brown and modified Lade. Regardless of the yield criterion, the
determination of the critical bottomhole pressure at which compressive shear yield initiates is
solved in STABView in an iterative process. Starting at a low value of bottomhole pressure (pw),
equations 4.2.16 through 4.2.23 are evaluated for 0 ≤ θ ≤ 180°. If yield is predicted at more than
one point, this indicates that pw falls below the critical value then pw is increased by a small
increment, and the process is repeated until a condition is reached where yield is predicted at
only one point.
x
y
Plane tangent to the
borehole surface
θ
σtmin
σr
σtmax
Shales, by definition, possess fissile, weak bedding planes. If these bedding planes are closely
spaced relative to the diameter of a borehole (as illustrated in Figure 4.20), shear yielding on
these weak planes can have a severe, negative impact on borehole stability. The effect of
pervasive planes of weakness on borehole stability is accounted for in STABView’s 3D elastic
models are as follows:
1. For a given point on the borehole wall, the complete stress tensor in the borehole
coordinate system (i.e., a Cartesian coordinate system with the z-axis aligned parallel to
the borehole axis) is calculated using equations 4.2.16 through 4.2.19.
2. The critical pressure for the initiation of compressive shear yielding, based on intact rock
strength, is calculated as described in equation 4.2.24. The value of this pressure is stored
as a parameter called pintact.
3. Using 3D stress transformation algorithms available in most rock mechanics textbooks
(e.g., Jaeger and Cook, 1979), the stress tensor is transformed from the borehole
coordinate system into a second coordinate system (x′, y′, z′) in which the z′-axis is
normal to the plane of weakness.
4. The critical bottomhole pressure for the initiation of shear yielding on the plane of
weakness (ppow) is calculated, using the following equations:
• σz′ in the new coordinate system represents the normal stress acting on the plane
of weakness (termed σnpow).
• The maximum shear stress acting on the plane of weakness (τpow) is:
• Assuming a linear Mohr-Coulomb yield criterion for the plane of weakness, with
cohesion cpow and friction angle φpow, the borehole pressure at which the near-well
pore pressure pa satisfies the following condition is found:
σnpow
τpow
When heated or cooled, rock has a tendency to expand or contract. Consequently, in situations
where the fluid temperature in a borehole is greater than the in-situ temperature of the
surrounding rock, an increase in compressive stress will occur in any plane tangent to the
borehole surface, resulting in an increased risk of compressive shear yielding. Conversely, if the
fluid temperature in a borehole is cooler than the surrounding rock, a stress reduction will be
induced and the risk of compressive shear yielding will decrease. As shown in Figure 4.21 both
heating and cooling scenarios may occur during drilling operations. In such cases, when
calculating the critical bottomhole pressure at which compressive shear yield initiates, the
equations presented in section 4.2.3 must be corrected as follows:
(σ θ )corr = σ θ + AT ∆T (4.2.27)
(σ z )corr = σ z + AT ∆T (4.2.28)
where:
STABView 3.8 User’s Manual 41 November, 2008
(a) (b)
Temperature Stress
Hot
T of mud tangential stress
r
Depth
Mud Temperatures Thermal Borehole Stresses
(after Dusseault, 1994)
Figure 4.21 Typical temperature profile while drilling a well and the resultant thermally
induced stress changes around a borehole.
AT = thermoelastic constant
Eα T
=
1−υ
αT = linear thermal expansion coefficient
∆T = bottomhole temperature – formation temperature
Note: Thermal effects are currently only accounted for in STABView for stress calculations
conducted at the borehole wall. As such, when calculating collapse pressures using the default
3D linear elastic model, the thermal option is available. However, when calculating yielding
risks (usually some parameter related to yielded zone size) for a specified bottomhole pressure,
the thermal option is not available.
During tripping, drill pipe run out of the borehole too quickly can generate a sudden reduction in
fluid pressure referred to as a swab pressure. This swabbing effect reduces the magnitude of the
support pressure that the drilling mud applies to the borehole wall, hence increasing the risk of
compressive shear yielding. For the elastic tensile fracturing models in STABView, it is assumed
that the decrease in pressure (∆pswab) occurs instantaneously, then dissipates rapidly. The
mechanical effect of this swab pressure on near-well stresses is assumed to be instantaneous.
However, it is assumed that there is insufficient time for the short-lived, transient pressure pulse
STABView 3.8 User’s Manual 42 November, 2008
to affect fluid pressures within the surrounding rock. Therefore, the near-well pore pressure
regime is assumed to be unchanged.
To account for a swab pressure in calculating the bottomhole pressure at which compressive
shear yield initiates, the stresses given in section 4.2.3 must be corrected as follows:
Predictions made with linear elastic models have a tendency to be too conservative for many
practical well engineering problems. The principal reason for this is the fact that, in many cases,
the initiation of compressive shear yielding alone does not necessarily lead to hole collapse that
constitutes a drilling problem or to sand production that is a serious issue. In order to retain the
relative ease of use of elastic models, while obtaining more realistic predictions, the following
approaches may be used in STABView.
Critical Breakout Angle Criterion: The procedures described in Sections 4.2.2 and 4.2.3
predict the critical bottomhole pressure at which compressive shear yielding initiates at
two diametrically opposed points on the borehole or perforation wall. An alternative
approach is to calculate the bottomhole pressure at which compressive shear yielding
occurs over a specified portion of the wall. A parameter commonly used to quantify this
type of yielding distribution is the breakout angle (θbrkt), which is illustrated in Figure
4.22. The rationale behind the critical breakout angle criterion is that the borehole or
perforation will be stable as long as there is sufficient unyielded rock on the wall to
provide a stable support structure that prevents catastrophic collapse. Some researchers
have gone as far as to recommend allowing a certain amount of breakout to purposely
restrict the development of a large yielded zone (see Zhang, 1998).
Calibration: If the sand production or borehole collapse pressure for a given depth and
well trajectory is known from previous experience, output from linear elastic models can
be calibrated to provide predictions that are consistent with this experience. One method
is to use a linear multiplier, in which case the ratio of the known pressure to the
theoretical pressure calculated with the selected model is determined. Subsequently, all
of the calculated output is multiplied by this ratio. Alternatively, a linear offset method
may be used, in which the difference between the known pressure and the theoretical
pressure calculated with the model is determined. Subsequently, this difference is
subtracted from all of the calculated output.
STABView 3.8 User’s Manual 43 November, 2008
Borehole or perforation
Yielded rock
θbrkt
Users of STABView have the option of selecting a critical breakout angle which is input as the
percentage of yielding on the borehole circumference, or the linear or offset calibration
procedure to “tune” a conservative 3D elastic model prediction to realistic field performance.
The linear elastic models presented in the preceding sections fail to account for an important
factor in determining the stability of a borehole or perforation. After rock has yielded, it may
weaken, yet retain some amount of strength, known as the residual strength. Stresses are
redistributed as yielding occurs, resulting in a low-stress zone of damaged rock near the borehole
or perforation (see Figure 4.23). The load supported in the yielded zone can play an important
role in stabilizing a borehole or perforation.
A method that has proven to provide more realistic wellbore stability and sand production
predictions than linear elastic models involves the use of elasto-plastic model. These models
enable the prediction of the bottomhole pressure at which some rock yielding occurs, yet the
yielded zone size is small enough that borehole or perforation stability is maintained (McLellan
et al., 2000).
For idealized conditions, closed-form solutions are available for predicting yielded zone size
around boreholes or perforations. Figure 4.24 shows an idealized elastic-brittle-plastic material
behavior model (which represents a special class of elasto-plastic behavior), which is a
reasonable approximation of material behavior in many rocks such as poorly cemented
sandstones, stiff shales, coal and microfractured carbonates. For this type of material behavior,
the most important segment of the stress-strain curve is the post-peak portion. As such, linear
elastic behavior is assumed until peak strength is reached. Then, instantaneous strain weakening
is assumed, followed by “plastic” behavior at a constant residual stress level.
STABView 3.8 User’s Manual 44 November, 2008
σ θ
STRESS
σ r
RADIAL DISTANCE
yielded rock
elastic rock
rw
R
Several authors, including Wang and Dusseault (1991), McLellan and Wang (1994), and Hawkes
and McLellan (1996) have published closed-form solutions for elastic-brittle-plastic materials,
with pore pressure effects included. Additional assumptions required by their solutions include:
• the borehole or perforation axis is parallel to one of the principal in situ stresses;
• the principal stresses normal to this orientation are equal;
• the peak and residual strengths of the rock can be represented by linear Mohr-Coulomb
criteria;
• plane strain conditions prevail;
• strains are small;
• the material is a homogeneous, isotropic continuum; and
• the rock and pore fluid are at constant temperature.
The general solution for such conditions was obtained by Hawkes and McLellan (1996) by
evaluating the following equation for yielded zone radius, R:
STABView 3.8 User’s Manual 45 November, 2008
peak strength
Deviator Stress
residual strength
Actual
Idealized
Axial Strain
1 ⎡ ⎤
{σ Ave − α p( R, t )}(1 − N p ) + S p + Ap ⎧⎨ o ( N p − 1) + p( R, t )⎫⎬⎥
p
⎢
N p +1 ⎣ ⎩2 ⎭⎦
N r −1
Sr ⎡ S r ⎤⎛ R ⎞ R
+ − ⎢ pw + ⎥⎜ ⎟ − α p (1 − N r ) R N r −1 ∫ r − N r p (r , t ) dr (4.3.1)
1− Nr ⎣ 1 − N r ⎦⎜⎝ rw ⎟⎠ rw
Ap po
+ σ Ave − =0
2
where Nr and Sr are residual strength parameters defined as follows:
1 + sin(φ r )
Nr = (4.3.2)
1 − sin(φ r )
− 2c r cos(φ r )
Sr = (4.3.3)
1 − sin(φ r )
and:
σAve = in-situ stress in the plane normal to the borehole or perforation axis,
αp = poroplastic coefficient, which is an effective stress parameter for yielded rocks
that is analogous to Biot’s coefficient for elastic rocks. There is limited theoretical
information or laboratory data pertaining to αp, and it is generally assumed to
have a value of 1.0.
rw = borehole radius or perforation radius, (rperf) depending on the completion method
for the well being analyzed.
STABView 3.8 User’s Manual 46 November, 2008
For complex near-well pore pressure regimes (e.g., transient flow conditions; complex near-well
permeability distributions; compressible fluid flow such as natural gas), the integral in equation
4.3.1 must be evaluated numerically. However, for steady-state radial flow of an incompressible
fluid, equation 4.3.1 simplifies to the following form:
⎡ ⎛ ⎞ ⎤
N r −1 (1 − N r )⎢2σ Ave + 2Y ⎜⎜ p a − A ln R ⎟⎟ − Ap po − α BG − UCS ⎥ − (1 + N p )(α p A − S r )
⎛R⎞ ⎣ ⎝ rw ⎠ 2G + λ ⎦
⎜⎜ ⎟⎟ =
⎝ rw ⎠ [ ]
S r − α p A + (1 − N r )( p w − α p p a ) (1 + N p )
.................................................................................................................................... (4.3.4)
where the following parameters describe rock mechanical properties:
α (N p − 1) + A p − α p (N p + 1)
Y= (4.3.5)
2
G = shear modulus
E
= (4.3.6)
2(1 + υ )
λ = Lamé’s elastic constant
Eν
= (4.3.7)
(1 + ν )(1 − 2ν )
and the following parameters describe near-well pore pressures:
pa − po
A= (4.3.8)
⎛r ⎞ ⎛r ⎞
k r ln⎜ d ⎟ − ln⎜ w ⎟
⎝R⎠ ⎝R⎠
pa = pore pressure at the borehole (or perforation) wall
po = native pore pressure, formation pressure, or reservoir pressure
B = kr A
ky = yielded rock permeability
ke = elastic (intact) rock permeability
kr = permeability ratio (ky / ke)
The parameters defining near-well pore pressures are shown in Figure 4.25. As illustrated, the
solution given in equation 4.3.4 is capable of handling situations where yielded rock
permeability is different from elastic rock permeability (e.g., permeability increase due to
microcracking during yielding), and where there is a step change in pore pressure at the borehole
or perforation wall (i.e., pw ≠ pa). The latter capability is well-suited to addressing situations
STABView 3.8 User’s Manual 47 November, 2008
where residual filter cake, or other types of wellbore skin mechanisms, result in permeability
impairment at the borehole or perforation wall.
Alternatively, the step-change pressure drop feature may be used as a simple means of analyzing
the complex problem of sand production during transient flow conditions. Specifically, in cases
where bottomhole pressure is reduced very suddenly and/or reservoir permeability is relatively
low, very steep pore pressure gradients will initially develop. Numerical or semi-analytical
techniques are generally required to rigorously model this transient stage of fluid flow, during
which there is an increased risk of sand production (albeit this sand production will likely be
short-lived). A simple means of addressing the worst-case scenario for transient flow conditions
is to assume that the bottomhole pressure drop will be applied instantaneously, resulting in a step
change in pressure at the borehole wall. Such a condition can be simulated using equation 4.3.4,
by setting pa = po. This end member of the transient case is available in STABView as an
advanced option.
Strictly speaking, the solutions presented in equations 4.3.1 and 4.3.4 are only valid for cases
where the in-situ stress state in the plane normal to the borehole or perforation axis is isotropic.
However, it has been demonstrated that these equations can be used to calculate the average
radius of the yielded zone for cases where the stress state is anisotropic. In such cases, the
average of the two principal stresses in the plane normal to the borehole or perforation axis can
used. The average yielded zone radius (as illustrated in Figure 4.26) is often adequate as a sand
production risk parameter. Equations for calculating the predicted length and width of the elasto-
plastic yielded zone can be found in Detournay and St. John (1988). An extension to these
equations is used in STABView for estimating elliptical hole sizes and shapes.
⎛r ⎞ ⎛r ⎞
p = pa − A ln⎜⎜ ⎟⎟ p = po − B ln⎜⎜ ⎟⎟
⎝ rw ⎠ ⎝ rd ⎠
p
po
ky ke
pa
pw
r
rw R rd
The cylindrical geometry assumed for perforations in the preceding sections is not well-suited
for modeling short, large-diameter perforations. Bratli and Risnes, R. (1981) have derived an
elasto-plastic model wherein perforations are assumed to have a hemi-spherical geometry, as
illustrated in Figure 4.27. Their model assumes elastic perfectly-plastic behavior (i.e., no strain
weakening) and an isotropic in-situ stress state (i.e., σV = σHmax = σHmin = σAve). The radius of the
yielded zone (R) can be evaluated for this model using the following expression:
T +1
⎛ R ⎞
T
3 (σ Ave − po )
⎜ ⎟ = T +3 (4.3.9)
⎜r ⎟ T +1 µQ
⎝ perf ⎠ 4c p tan γ −
T 4πkrperf
where:
rperf = perforation radius
γ = failure angle
π φp
= +
4 2
T = 2{tan 2 (γ ) − 1}
µ = fluid viscosity
rw
R
a
Figure 4.26: Elliptical yielded zone predicted for an anisotropic stress regime,
and the equivalent average yielded zone radius, R.
STABView 3.8 User’s Manual 49 November, 2008
σAve
σ
σr σθ
r rperf σAve
σAv σθ
σr
y
σAve
θ
pw
po
x
po pw
σAv rperf rd
R
r
Figure 4.27: Geometry and typical stress distribution for elasto-plastic modeling of
hemi-spherical perforations (after Bratli and Risnes, 1981).
STABView users have the ability to analyze hemi-spherical perforations using this model by
selecting the perforation stability option under the class of sand production models. Figure 4.28
shows the additional information the user may provide to more accurately characterize the
perforating design. A future version of STABView will display the 3D nature of multiple
perforation tunnels with varying orientations, phasing, size, depth and yielding.
Horizontal wells which develop an elliptical yielded zone under drilling or production conditions
can result in the detachment of this yielded or failed material due to gravity. This “rubble” will
accumulate on the bottomside of the hole which may or may not have a tubular present such as
drill pipe or a slotted liner, for instance. It is possible to provide a first-order calculation of the
cross-sectional area of yielded rock which may detach and then bulk up by a factor βr . Using
theory developed for the flow of particulates in bins, (Dreschler, 1991) the wedge angle (αw) can
be estimated as a function of the yielded rock’s residual friction angle (φr) and a roughness
angles (φi) which reflects the uneven detachment surface.
STABView 3.8 User’s Manual 50 November, 2008
π (φ r + φ i )
αw = − (4.3.10)
2 2
Figure 4.29 shows a wellbore cross-section where a slotted liner is depicted on the bottomside of
a borehole, which exhibits a large elliptical yielded area under drawdown conditions. The
yielded shape is predicted with STABView’s 2D elastoplastic model with peak and residual
failure parameters. STABView presently allows the user to calculate the Rubble Fill Percentage
(RFP) for this problem, without accounting of the presence of the liner. RFP is defined as
follows:
A1
RFP = 100β r (4.3.11)
πrw 2
where:
ro
Elastoplastic
In Section 4.3.1, a closed-form solution for predicting yielded zone size around boreholes or
perforations for idealized elasto-plastic material behavior was described. This idealized material
model involves linear elastic behavior until peak strength is reached, followed by instantaneous
strain weakening, then “plastic” behavior at a constant residual stress level. If a constant angle of
dilation (ψd) is assumed during this plastic behavior (see Figure 4.30), it is possible to calculate
the strain and displacement magnitudes at the borehole or perforation wall.
The equations presented below are modified after Wang and Dusseault (1991a) and Kooijman et
al. (1996). Additional assumptions required by their solutions include:
• the borehole or perforation axis is parallel to one of the principal in situ stresses;
STABView 3.8 User’s Manual 52 November, 2008
peak strength
Deviator Stress
residual strength
Axial Strain
Volumetric Strain
Expansion
⎛ 2 sinψ d ⎞
arctan ⎜⎜ ⎟⎟
⎝ 1 − sinψ d ⎠
Axial Strain
Contraction
• the principal stresses normal to this orientation are equal (i.e., radial symmetry);
• the peak and residual strengths of the rock can be represented by linear Mohr-Coulomb
criteria;
• plane strain conditions prevail;
• strains are small;
• the material is a homogeneous, isotropic continuum; and
• the rock and pore fluid are at constant temperature.
For radially symmetric borehole geometries, radial displacement (u) of any point on the borehole
wall can be calculated as follows:
u = ε θT rw (4.3.12)
where:
ε θT = total tangential strain
rw = borehole (or perforation) radius
( ) ( )
Total tangential strain is the sum of elastic ε θe and plastic ε θp tangential strains; i.e.,
ε θT = ε θe + ε θp (4.3.13)
STABView 3.8 User’s Manual 53 November, 2008
The elastic component of tangential strain can be calculated from stresses at the borehole or
perforation wall using Hooke’s law, as follows:
1+υ
ε θe =
E
[ {
(1 − υ ) σ θ' (rw ) − υσ r' (rw ) }] (4.3.14)
where:
υ = Poisson’s ratio
E = Young’s modulus
σ θ' (rw ) = effective tangential stress at rw
σ r' (rw ) = effective radial stress at rw
The plastic component of tangential strain is derived from the compatibility equation and the
plastic flow rule, providing the following result:
ε θp = C1 + C 2 + C3 (4.3.15)
where:
⎡ 2
⎤
(1 − υ ){sin (ψ d ) − 1} [
N r −1+
⎛ ⎞
(N r − 1)σ r' (rw ) − S r + A]⎢⎢1 − ⎜⎜ R ⎟⎟
1− sin (ψ d )
C1 = ⎥ (4.3.16)
{1 − sin (φ r )sin (ψ d )}2G ⎝ rw ⎠ ⎥
⎣⎢ ⎦⎥
⎡ 2
⎤
−A ⎛ A p ⎞⎢ ⎛ R ⎞ 1−sin (ψ d ) ⎥
C2 = (1 − υ ){sin (ψ d ) − 1}⎜⎜1 − ⎟⎟⎢1 − ⎜⎜ ⎟⎟
⎥
(4.3.17)
2G 2 ⎠ r
⎝ ⎢⎣ ⎝ w ⎠ ⎥⎦
2
C3 =
(1 − υ ) (S ⎛ R ⎞ 1−sin (ψ d )
− S p )⎜⎜ ⎟⎟ (4.3.18)
r
2G ⎝ rw ⎠
and:
ψd = dilation angle
G = shear modulus
E
=
2(1 + υ )
φr = peak friction angle
φr = residual friction angle
cp = peak cohesion
STABView 3.8 User’s Manual 54 November, 2008
cr = residual cohesion
1 + sin(φ r )
Nr =
1 − sin(φ r )
− 2c p cos(φ p )
Sp =
1 − sin(φ p )
− 2c r cos(φ r )
Sr =
1 − sin(φ r )
Ap = poroelastic constant
=α
(1 − 2υ )
(1 − υ )
pa − po
A=
⎛r ⎞ ⎛r ⎞
k r ln⎜ d ⎟ − ln⎜ w ⎟
⎝R⎠ ⎝R⎠
pa = pore pressure at the borehole (or perforation) wall
po = native pore pressure, formation pressure, or reservoir pressure
B = kr A
ky = yielded rock permeability
ke = elastic (intact) rock permeability
kr = permeability ratio (ky / ke)
R = yielded zone radius
rd = drainage radius
An algorithm for calculating the stress on a slotted liner or screen is given by Kooijman et al.
(1996). This solution is valid for openhole completions in which elastoplastic deformation of the
borehole wall has closed the gap between the sandface and the slotted liner or screen. The model
makes use of a number of simplifying assumptions, most notably that the slotted liner or screen
is rigid, and that the in-situ stress state in the plane normal to the borehole axis is isotropic
(hence the deformation of the borehole wall is axisymmetric).
In principal, the liner load algorithm is the same as the solutions given above for yielded zone
radius and near-well stresses (Section 4.2.2) and borehole wall deformations (this section). The
only difference is that, when there is contact between the liner and the sandface, the effective
radial stress at the borehole wall is exerted by the liner rather than the borehole fluid(s).
The liner load solution algorithm works as follows:
1. Generate an initial guess for the yielded zone radius (R).
STABView 3.8 User’s Manual 55 November, 2008
( )
2. Calculate the effective radial stress at the borehole wall ( σ r rw ).
'
( )
3. Calculate the total tangential strain at the borehole wall ( ε θ rw ).
T
5. If u is equal to the annular gap between the borehole wall and the slotted liner or screen,
the effective radial stress calculated in step 2 represents the magnitude of the stress acting
on the liner.
6. If u is NOT equal to the annular gap between the borehole wall and the slotted liner or
screen, adjust the guessed value for R and repeat steps 2 through 5 until a solution is
found.
The elastic models presented in Section 4.2 can also be used to predict tensile yielding or
fracture initiation, which may occur while overbalanced drilling, completing or stimulating a
well, or during water injection. An example of the radial, tangential and axial stresses around a
vertical borehole is given in Figure 4.31. These effective stresses were calculated using equations
6000
σHmin
radial stress
5000 tangential stress
4000
3000
2000
1000
Compression
Tension
-1000
0 30 60 90 120 150 180 210 240 270 300 330 360
Figure 4.31 Example showing the distribution of effective stresses around the
circumference of a borehole in which effective tensile stresses have developed.
STABView 3.8 User’s Manual 56 November, 2008
4.2.1 through 4.2.3 for the following input parameters, which are representative of a drilling case
with a weighted mud system:
3σ Min − σ Max − Ap po + σ T
pb = (4.4.3)
(1 + α − A )
p
In these equations, σT is the rock tensile strength. For partially-penetrating fluids, the magnitude
of the pressure drop (∆p) at the borehole wall can be related to the formation of a filter-cake, or
wall coating as described in Section 4.6 of this manual.
For boreholes that are not drilled parallel to an in-situ stress direction, the equations for borehole
stresses presented in Section 4.2.3 may be used to evaluate fracture breakdown pressure. Using
equation 4.2.21 for σtmin (i.e., the smallest stress in the plane tangent to the borehole surface), the
criterion for fracture breakdown at a given point on the borehole wall is:
STABView 3.8 User’s Manual 57 November, 2008
1
σθ + σ z ⎡ 2 ⎛ σ θ − σ z ⎞2 ⎤ 2
σ t min = − ⎢τ θz + ⎜ ⎟ ⎥ = −σ T (4.4.4)
2 ⎢⎣ ⎝ 2 ⎠ ⎥⎦
− C B − C B2 − 4C A C C
pb = (4.4.5)
2C A
where:
CA =
1
4
[
(2 A p ε fp − 1 − 2αε fp )2 − 1 ]
C B = C 4 (2 A p ε fp − 1 − 2αε fp ) + C 5
C C = C 42 − C 32 − C 52
C1 − C 2
C5 =
2
C1 + C 2 − 2 A p ε fp p o
C4 = − αp o (1 − ε fp ) + σ T
2
C 3 = 2τ yz cos θ − 2τ xz sin θ
[
C 2 = σ zz − υ 2(σ xx − σ yy ) cos 2θ + 4τ xy sin 2θ ]
C1 = σ xx + σ yy − 2(σ xx − σ yy ) cos 2θ − 4τ xy sin 2θ
and εfp denotes the fluid penetration coefficient, which will be discussed further in Section 4.6.
This parameter ranges from 0.0 for a non-penetrating fluid to 1.0 for a fully-penetrating fluid.
To determine the fracture breakdown pressure at a given depth in a well, equation 4.4.14 is
evaluated over the range 0 ≤ θ ≤ 180°, and the smallest value found is retained as the solution.
The deep penetrating fluid model is similar to the default linear elastic fracture breakdown model
described above. The fundamental difference is that the deep penetrating fluid model is
applicable for cases where the extent of fluid penetration during injection is very large compared
to the thickness of the reservoir, resulting in a poro-elastic response different from the default
model.
STABView 3.8 User’s Manual 58 November, 2008
The solution for fracture breakdown pressure using the deep fluid penetration option is given as
follows:
− C B − C B2 − 4C A C C
pb = (4.4.6)
2C A
where:
CA =
1
4
[
(2 A p ε fp − 1 − 2αε fp )2 − C 6 ]
C B = C 4 (2 A p ε fp − 1 − 2αε fp ) + C 5 C 6
C C = C 42 − C 32 − C 52
C 6 = −2 A p ε fp + 1
C1 − C 2 − 2 A p ε fp p o
C5 =
2
and the definitions of constants C1 through C4 are identical to those given for the default fracture
breakdown model.
The modified poro-elastic effect for deep fluid penetration was taken from "Waterflooding
Under Fracturing Conditions" by E.J.L. Koning (Ph.D. Thesis, Delft University of Technology,
1988). Note that Koning's equations were derived for vertical wells in horizontal reservoirs. The
same equations have been applied to the general 3D problem in STABView. As such, this model
is most applicable for cases where fluid flow is radially symmetrical about the borehole, and
occurs in a plane normal to the borehole axis.
This model computes the fracture breakdown pressure for an inclined borehole in a setting with
uniform in-situ horizontal stresses. This model is based on the semi-analytical fracture
mechanics model derived by Morita et al. (1996). This model is useful for borehole fluids that
contain fracture plugging additives, such that fracture breakdown does not occur until a pre-
existing or induced tensile crack at the borehole wall reaches a critical aperture, at which point
whole drilling mud is able to penetrate the crack.
The geometric parameters relevant to this model are shown in Figure 4.32. The solution process
is as follows:
1. An initial guess for the critical crack length (A) is made.
2. For the specified crack length, the fracture breakdown pressure (pb) is calculated.
STABView 3.8 User’s Manual 59 November, 2008
3. For the specified crack length and fracture breakdown pressure, the aperture of the crack
(W) is calculated.
4. If the calculated crack aperture corresponds to the critical aperture at which whole mud
can penetrate the crack, the fracture breakdown pressure calculated in step 2 is the correct
value and the calculation process is complete.
5. If the calculated crack aperture does not correspond to the critical aperture at which
whole mud can penetrate the crack, a revised guess for crack length is made and steps 2
through 4 are repeated as required.
For relatively short crack lengths, the equation required for step 2 is:
Kc
p b = 3σ Min − σ Max + (4.4.7)
1.215 π ( A − rw )
Where Kc is the fracture toughness of the rock. For longer crack lengths, the equations for pb can
be found in Table 1 of Morita et al. (1996).
The equations required for step 3 are:
W = W1 + W 2 + W3 (4.4.8)
where:
(
W1 = 4 1 − υ 2 ) AG
E
1
σ Min
(
W2 = 4 1 − υ 2 ) AGE (σ
2
Max − σ Min )
(
W3 = 8 1 − υ 2 ) πrE ⎡⎢G
w
3
⎛ A ⎞⎤
+ ln⎜⎜ ⎟⎟⎥ p b
⎢⎣ ⎝ rw ⎠⎥⎦
and the values of the fitting functions G1 through G3 can be found in Table 1 of Morita et al.
(1996).
The version of this model implemented in STABView assumes that no mud pressure penetration
occurs, and that isothermal conditions exist. Variations on this model that account for these
effects are also included in Morita et al. (1996).
Note: For cases where the horizontal stresses are unequal, STABView will calculate fracture
breakdown pressures for this model using the average of σHmax and σHmin.
STABView 3.8 User’s Manual 60 November, 2008
σMax
σMin
W
pb
po
rw
This model computes the stresses around a borehole and calculates the bottomhole pressure at
which a tensile fracture will initiate. These calculations are based on the model for a transversely
isotropic medium, as derived in "Stresses Around Horizontal Boreholes Drilled in Sedimentary
Rocks" by B.S. Aadnoy (SPE paper 17119, July, 1987). This model assumes that Young's
modulus is constant in the medium's symmetry plane (e.g., bedding), but different from the value
normal to the symmetry plane. Poisson's ratio is assumed to be isotropic. The model is valid for
wells that are drilled within the symmetry plane. As such, it is suitable for the case of a
horizontal well in flat or dipping strata.
If the temperature in a well is different from the in-situ temperature of the surrounding rock,
thermally-induced stress changes will result. For example, it is common during water injection
for the bottomhole temperature to be lower than the reservoir temperature. As such, the rock
around the borehole will tend to contract, resulting in tensile stress changes. To accurately
predict the fracture breakdown pressure for such cases, it is necessary to account for this thermal
effect.
For the default elastic tensile fracturing model, it is assumed that the temperature change is
confined to the near-well area. This scenario is representative, for example, of a drilling case in a
low-permeability formation, or drilling in a high-permeability formation in which limited leakoff
STABView 3.8 User’s Manual 61 November, 2008
occurs due to the development of an effective filter cake. For the default tensile fracturing model,
the tangential and axial stresses at a point on the borehole wall are given as follows:
where:
AT = thermoelastic constant
Eα T
=
1−υ
αT = linear thermal expansion coefficient
∆T = bottomhole temperature – formation temperature
Radial and shear stresses are unaffected by temperature effects, and can be calculated using
equations 4.2.16 and 4.2.19, respectively. The procedure for calculating fracture breakdown
pressure is unchanged from Section 4.4.1, except the constant C4 must be corrected as follows:
(C 4 )corr = C 4 + AT ∆T (4.4.11)
For the linear elastic stress model with deep fluid penetration, it is assumed that the wellbore
fluid flows deep into the surrounding rock, such that the rock temperature becomes equal to the
bottomhole temperature. For this fracturing model, the tangential stress at a point on the borehole
wall is given as follows
Radial, axial and shear stresses are unaffected by temperature effects for this model, and can be
calculated using equations 4.2.16, 4.2.18 and 4.2.19, respectively. The procedure for calculating
fracture breakdown pressure is unchanged from Section 4.4.2, except the constants C4 and C5
must be corrected as follows:
(C 4 )corr = C 4 + AT ∆T (4.4.13)
(C 5 )corr = C 5 + AT ∆T (4.4.14)
During tripping, drill pipe run rapidly into the borehole can generate large surge pressures in the
borehole-drillpipe annulus which can lead to tensile fracturing of the formation. For the elastic
tensile fracturing models in STABView, it is assumed that the increment in pressure (∆psurge)
occurs instantaneously, then dissipates rapidly. The mechanical effect of this surge pressure on
near-well stresses is assumed to be instantaneous. However, it is assumed that there is
STABView 3.8 User’s Manual 62 November, 2008
insufficient time for the short-lived, transient pressure pulse to affect fluid pressures within the
surrounding rock. Therefore, the near-well pore pressure regime is assumed to be unchanged.
To account for a surge pressure in calculating the fracture breakdown pressure using the default
elastic fracturing model, constants C4 and C5 given in Section 4.4.1 must be corrected as follows:
∆p surge
(C 4 )corr = C4 − (4.4.15)
2
∆p surge
(C 5 )corr = C5 − (4.4.16)
2
Similarly, to account for a surge pressure in calculating the fracture breakdown pressure using
the deep fluid penetration fracturing model, constants C4 and C5 given in Section 4.4.2 must be
corrected as follows:
∆p surge
(C 4 )corr = C4 − (4.4.17)
2
∆p surge
(C 5 )corr = C5 − (4.4.18)
2
This model computes the stresses around a borehole drilled in an arbitrary orientation, and uses
these results to assess the risk of passive shear failure initiation at the borehole wall. In some
cases, especially for very weak, unconsolidated formations, passive shear failure can occur at
bottomhole pressures lower than the tensile fracture initiation pressure. The creation of a zone of
rock that has yielded in shear can affect tensile fracture initiation. The assumptions and
limitations of this model are similar to the tensile fracture initiation model.
The term passive shear yielding refers to a failure mode in which radial stress is large, and the
minimum stress in the plane tangent to the borehole surface is small, such that compressive shear
yielding occurs. The solution for passive shear yielding is found using equations 4.2.16, 4.2.21
and 4.2.24, substituting σr for σ1 and σtmin for σ3. At a given point on the borehole wall, the
critical pressure for the initiation of passive shear yielding can be found as follows:
− C B − C B2 − 4C A C C
pw = (4.5.1)
2C A
where
Np
C A = C5 −
4
STABView 3.8 User’s Manual 63 November, 2008
C B = 2C 4 C 5 + N p2 C 6
(
C C = C 42 − N p2 C 32 − C 62 )
C1 − C 2 − ∆p surge
C6 =
2
2 − 2αε fp + N p (1 − 2 A p ε fp + 2αε fp )
C5 =
2
⎡ C1 + C 2 − ∆p surge − 2 A p ε fp p o ⎤
C 4 = αp o (ε fp − 1) + ∆p surge + S p − N p ⎢ + αp o (ε fp − 1)⎥
⎣ 2 ⎦
and the definitions of constants C1 through C3 are identical to those given for the default fracture
breakdown model (Section 4.4.1).
To determine the critical pressure for passive shear yielding at a given depth in a well, equation
4.5.1 is evaluated over the range 0 ≤ θ ≤ 180°, and the smallest value found is retained as the
solution.
Note: The solution given above accounts for surge pressure effects. Thermal effects are not
currently accounted for in the passive shear yielding model implemented in STABView.
Near-well pore pressures have a strong influence on stresses, deformations and yielding in the
surrounding rock. A convenient parameter for characterizing near-well pore pressures is the filter
cake or wall coating efficiency (εfc). Specifically, this parameter characterizes the extent to
which the pore pressure at the borehole wall increases in response to an overbalance pressure.
Mathematically, the filter cake or wall coating efficiency is defined as follows:
pw − pa
ε fc = (4.6.1)
p w − po
A filter cake or wall coating efficiency of 1.0 indicates that no mud pressure penetration will
occur, whereas a value of 0.0 indicates that no filter cake or wall coating is present. Figure 4.33
illustrates a typical near-well pore pressure profile, for an overbalanced drilling case in which the
filter cake or wall coating efficiency has a value of roughly 0.5.
Equation 4.6.1 can be re-worked to solve for other pressure parameters that have been used in
various models described this chapter. Notable variables of this type include the pressure right at
the borehole wall:
STABView 3.8 User’s Manual 64 November, 2008
p a = p w − ε fc ( p w − p o ) (4.6.2)
∆p = p a − p o = ( p w − p o )(1 − ε fc ) (4.6.3)
The fluid penetration coefficient (εfp), which is specified for fracture or lost circulation problems,
is defined as follows:
pa − po
ε fp = (4.6.4)
p w − po
A fluid penetration coefficient of 1.0 indicates a fully-penetrating fluid, whereas a value of 0.0
indicates a non-penetrating fluid. Numerically, the fluid penetration coefficient is equal to one
minus the filter cake or wall coating efficiency:
ε fp = 1 − ε fc (4.6.5)
εfc or εfp are key inputs for most of the algorithms used in STABView for borehole collapse and
fracture breakdown calculations. They influence the 3D elastic models that assess conditions at
the borehole wall, and they control the boundary conditions for the steady state flow equations
that are solved within the elastic and elasto-plastic models used to predict yielded zone size.
Filter cake
pw
Support
pa
po
rw r
Figure 4.33: Pressure drop across a filter cake.
STABView 3.8 User’s Manual 65 November, 2008
The capillary threshold pressure (pc) is the overbalance pressure that must be exceeded before a
non-wetting drilling fluid (typically oil- or synthetic-based) will penetrate a water-wet formation.
As described by the Laplace law, the magnitude of this threshold pressure is inversely
proportional to the radius of pore throats in the rock:
2γ cos ϕ
p cap = p non − wet − p wet = (4.6.6)
r pore
where:
pnon-wet = pressure in the non-wetting fluid phase (e.g., OBM)
pwet = pressure in the wetting fluid phase (e.g., native pore fluid)
γT = interfacial tension
ϕ = contact angle of the fluid interface where it meets the mineral grain surface
rpore = pore throat radius
Large capillary threshold pressures can exist for intact shales with small pore throat radii. If the
overbalance pressure is smaller than the capillary threshold pressure, no mud pressure
penetration occurs (Figure 4.34). Mathematically, this scenario is equivalent to drilling with a
mud that has a filter cake or wall coating efficiency of 1.0. In all of STABView models that use
this parameter, it is assumed that mud pressure penetration occurs (i.e., εfc = 0) once the capillary
threshold pressure has been exceeded.
Figure 4.34: Prevention of mud pressure penetration due to a high capillary threshold
pressure in shale (Hawkes et al., 2000).
STABView 3.8 User’s Manual 66 November, 2008
The near-wellbore fluid pressure parameter required for geomechanical modeling of sand
production and openhole stability is the pore pressure at the borehole wall or the “sandface”
pressure. For wells completed with slotted liners, screens or pre-packed screens, a pressure drop
will occur across the sand control device as illustrated in Figure 4.35. The magnitude of this
pressure drop may be significant, depending on the dimensions and total cross-sectional area
openings in the liner, screen or pre-pack, and the viscosity and flow rate of the produced fluid.
For analyses that are conducted for a specified bottomhole pressure, it is necessary to add the
liner/screen pressure drop to the bottomhole pressure in order to determine the appropriate
sandface pressure to use in the calculations.
Similarly, when modeling near-well stresses and yield conditions for specified fluid production
rates rather than specified bottomhole pressures, the sandface pressure required to achieve the
target flow rate must be calculated first. Subsequently, the corresponding bottomhole pressure
may be calculated by subtracting the liner/screen pressure drop.
STABView users can specify a value for the additional pressure drop from the formation “sand-
face” to the interior of the wellbore.
Screen Formation
p
po
pa
∆p
pw
r
OD rw
Figure 4.35: Pressure drop across a slotted liner, screen or pre-packed screen gives rise
to a sandface pressure that is different from the wellbore pressure inside the tubular.
The flow equations that are used by default within the elastic and elasto-plastic models for
yielded zone size calculations assume steady-state, radial flow of an incompressible fluid. This
assumption works well for water and oil, for which density variation with pressure is relatively
STABView 3.8 User’s Manual 67 November, 2008
small. However, given the highly compressible nature of natural gas, a different flow model is
required.
The radial flow equation for compressible fluids can be linearized approximately by replacing all
pressure terms with the so-called gas “pseudo-pressure”, m(p), which is defined as follows:
p
p
m( p ) = 2 ∫ dp (4.6.7)
p ref
µZ
where:
p = the pressure of interest
pref = a reference pressure (which can be chosen arbitrarily; often standard conditions)
µ = fluid viscosity
Z = gas compressibility factor
Steady-state solutions for radial inflow conditions, cast in terms of gas pseudo-pressure (mp),
have been presented by several investigators (e.g., Butler, 1994):
2
⎛ 2p ⎞⎛ Q sc T ⎞ r 2βMW ⎛ Q sc Tp sc ⎞ ⎛ 1 1⎞
m( p ) = m( p w ) + ⎜⎜ sc ⎟⎟⎜ ⎟ ln + ⎜⎜ ⎟⎟ ⎜⎜ − ⎟⎟ (4.6.8)
⎝ Tsc ⎠⎝ 2πkh ⎠ rw RT µ ⎝ 2πhTsc ⎠ ⎝ rw r ⎠
where:
m(pw) = gas pseudo-pressure at bottomhole pressure (pw)
psc = pressure at standard conditions
Tsc = temperature at standard conditions
T = reservoir temperature
Qsc = volumetric flow rate at standard conditions
T = reservoir temperature
k = reservoir permeability
h = reservoir thickness
r = radial distance
rw = borehole radius
β = turbulence factor
MW = gas molecular weight
STABView 3.8 User’s Manual 68 November, 2008
For a specified bottomhole pressure, pore pressure at any point around a well is evaluated by
solving equation 4.6.8 at the point of interest, then using an inverse form of equation 4.6.7 to
convert gas pseudo-pressure to gas pressure. To calculate yielded zone size for gas wells using
the 2D elastoplastic model, the integral in equation 4.3.1 is solved using Gaussian quadrature
(Gerald and Wheatley, 1989), wherein gas pressures are calculated at nodal points using the
procedure presented above.
As the rate of fluid flow through a porous medium increases, inertial and/or turbulence effects
become significant. In such cases, Darcy’s law alone does not adequately describe the
relationship between pressure gradient and fluid flow rate. The second term on the right hand
side of equation 4.6.8 accounts for non-Darcy flow effects, which are common for high-rate gas
wells. The relevant parameter in this term is the turbulence factor (β). This rock property
depends on factors such as porosity, pore size distribution, pore shape and tortuosity. Given that
these factors also control permeability, it is not surprising that empirical correlations exist for
estimating β as a function of k. An historical example based on data for several sandstones and
carbonates is given in Butler (1994) is:
ln β = 28.83 − 1.201 ln k (4.6.9)
The gas pseudo-pressure solution given in Section 4.6.4 requires the characterization of gas
viscosity and compressibility factor over a finite pressure range. Although they can be measured
in laboratory tests, these properties are commonly estimated using empirical correlations. The
only input parameter required of STABView users is the gas gravity (γg), which is the ratio of the
molecular weight of natural gas to that of air. The molecular weight of air is normally taken as
28.9625 g/gmol.
Pseudo critical temperature (Tpc) and pseudo critical pressure (ppc) of the gas are estimated from
the user-specified gas gravity using the correlation of Sutton (2005):
γg = gas gravity
STABView 3.8 User’s Manual 69 November, 2008
Note: Sutton’s equations are applicable for 0.55 < γg < 2.1 with < 2 mol% CO2, < 2 mol% N2
and < 2 mol% H2S. If non-hydrocarbon components are present, the gas gravity should be
based on the hydrocarbon components only.
The gas compressibility factor (Z) is then estimated using the correlation of Dranchuk and Abou-
Kassem (1975), as presented in Londono et al. (2002):
⎡ A A A A ⎤ ⎡ A A ⎤ 2
Z = 1 + ⎢ A1 + 2 + 33 + 44 + 55 ⎥ ρ pr + ⎢ A6 + 7 + 28 ⎥ ρ pr
⎣⎢ Tpr Tpr Tpr Tpr ⎥⎦ ⎣⎢ Tpr Tpr ⎦⎥
(4.6.l2)
⎡ A7 A8 ⎤ 5 ρ pr2 − A11 ρ 2pr
− A9 ⎢ + 2 ⎥ ρ pr + A10 (1 + A11ρ pr ) 3 e
2
where:
Z = gas compressibility factor
ppr = pseudo reduced pressure = p/ppc
Tpr = pseudo reduced temperature = T/Tpc
p pr
ρ pc = 0.27 (4.6.13)
Z ⋅ Tpr
An iterative procedure is used to solve for Z and ρpc that satisfy both equations 4.6.l and 4.6.13.
Finally, gas viscosity is estimated using the correlation of Lee at al. (1966) as reported in Sutton
(2005):
µ = Y1 exp(Y5)/10,000 (4.6.14)
where:
Y1 = (22.65 + 0.03873 MW)T 1.5
(209.2 + 19.26 MW + 1.8 T )
Y2 = 3.448 + 548 + 0.01009 MW
T
Y3 = 2.447 – 0.2224 Y2
Y4 = 0.4222 MW
10,000 Bg
Y5 = Y2 (Y4Y3)
The gas formation volume factor required for the above correlation is calculated as follows:
Bg = 0.351 Z T (4.6.15)
p
The flow of fluid (e.g., pore water or filtrate) and solute (e.g., dissolved ions or molecules) into
or out of shales can have a profound influence on near-wellbore pore pressures, stresses, and
deformations. The models presented in this section account for hydraulically-driven flow; i.e.,
STABView 3.8 User’s Manual 71 November, 2008
flow which is driven by fluid pressure gradients. Another mechanism which drives fluid flow,
especially for many shales, is referred to as “chemical osmosis” or simply “osmosis.” The basic
principle of osmosis revolves around the concept of a “semi-permeable membrane,” i.e., an
interface which is permeable to water molecules but impermeable to solute ions or molecules.
As illustrated in Figure 4.36, if a semi-permeable membrane separates solutions having different
solute concentrations, an “osmotic pressure potential” exists. Water will flow through the
membrane from the low concentration solution to the high concentration solution. More
specifically, water will flow from the solution having a high water activity to the solution having
a low water activity. As a consequence, fluid pressure in the concentrated solution will increase
until the counter-flow induced by the newly created hydraulic head balances the flow induced by
the concentration difference. At equilibrium, the latter hydraulic pressure difference is referred to
as the “osmotic pressure” of the fluid membrane system.
The recognition of chemical osmosis in shales has led to the common practice of attempting to
control water flow out of (or into) shales by adjusting the solute concentration of drilling muds.
Traditionally it had been believed that the osmotic process was only applicable for oil-based
emulsion drilling muds. In such cases, the film of emulsifier at the interface between water (or
brine) droplets is believed to act as the semi-permeable membrane. More recently, it has become
apparent that some shales may act to some extent as semi-permeable membranes when contacted
with water-based muds. The exact mechanism by which the shale acts as a membrane is not fully
understood, although it is likely related to the reduced mobility of ions or molecules of a given
size, charge or charge distribution within small clay pores bounded by clay particles having
negative surface charges. The efficiency of the solute exclusion process is not perfect in such
cases, so these shales are referred to as “non-ideal” or “leaky” membranes.
The following equation can be used to estimate the change in pore pressure at the borehole wall
due to chemical osmosis:
RT ⎛a ⎞
∆p osm = M eff ln⎜⎜ mud ⎟⎟
Vw ⎝ a shale ⎠ (4.7.1)
where:
To assess the effect of osmosis on near-well stresses and yielding, the following correction can
be applied to terms used throughout this chapter to characterize pore pressure conditions at the
borehole wall:
While rock mechanical factors can lead to rock yielding and the eventual detachment of failed
material into a wellbore, there is another significant mechanism for hole enlargement. Borehole
STABView 3.8 User’s Manual 73 November, 2008
Permafrost
Interval
T < 0°C
dh
db
Figure 4.37: Schematic showing the mechanism of hole enlargement while drilling
permafrost based on the model of Kutasov and Caruthers (1988).
erosion, due to the action of fluids circulating at sufficient velocity in an annulus, can also cause
intact or yielded rock to detach. The treatment of the hydraulic forces that will result in erosion
is a complicated subject and beyond the scope of this manual. STABView does, however, allow
the user to predict the potential for hole enlargement in the special case of a vertical borehole
through permafrost.
As shown in Figure 4.37 the erosive action of a warm drilling fluid circulating down drillpipe
and back along the annulus between a permafrost-rich formation and the drillpipe can cause the
borehole to “washout” to a new diameter with time. This new diameter is a function of many
factors including: the fluid’s temperature and thermal conductivity, the latent heat of fusion for
ice, the ice saturation and porosity of the formation, the circulation rate of the fluid, and it’s
rheological characteristics.
Appendix D of this manual describes a simple model for permafrost hole enlargement that is
available in STABView. This model is based on the published work of Kutasov and Caruthers
(1988) and Kutasov (1999).
If elevated fluid pressures in a borehole (e.g., due to drilling with a weighted mud system or due
to fluid injection) migrate into natural fractures or faults, the effective normal stress acting on
these discontinuities is reduced. Such a reduction may result in shear failure or slippage on these
discontinuities, which can negatively impact drilling, completion and production operations,
and/or compromise the hydraulic integrity of a reservoir.
STABView 3.8 User’s Manual 74 November, 2008
The following 2D stress transform can be used to calculate the normal stress (σnpow) and shear
stress (τpow) acting on a discontinuity inclined at an angle θ from the x-axis, when σx and σy are
principal stresses, and the strike of the discontinuity is parallel with the third principal stress (σz),
as illustrated in Figure 4.38:
Planar discontinuity
y
σnpow
σx τpow
θ
x
σy
Figure 4.38 Resolution of stresses acting on a planar discontinuity, for the 2D case.
STABView 3.8 User’s Manual 75 November, 2008
In a manner very similar to the induced slip risks described in Section 4.9, it is also possible to
re-open natural fractures and faults if elevated borehole pressures exceed the magnitude of the
stress component acting normal to the discontinuity plane. Depending on the discontinuity
orientation and the in-situ stress regime, it is possible for this to occur at pressures lower than the
critical slip pressure (pslip).
For the 2D case illustrated in Figure 4.38, the normal stress acting on a weak discontinuity can
be calculated using equation 4.9.1. The critical pressure required to re-open the discontinuity
(popen) is then given by:
p open = σ npow
(4.10.1)
STABView includes a model for fracture re-opening analysis. The methodology used for these
calculations follows directly from equations 4.9.1 and 4.9.4. The only difference is that a general,
3D stress transformation is used to resolve the normal stress acting on the discontinuity.
Equations for such stress transformations can be found in most rock mechanics textbooks (e.g.,
Jaeger and Cook, 1979).
Sections 4.2, 4.4, 4.5, 4.9 and 4.10 present the equations and methodologies for calculating the
critical bottomhole pressures at which various types of yielding or failure events initiate (e.g.,
active or passive shear yielding; tensile fracturing; re-opening or slip on weak discontinuities).
An alternative approach for analyzing these yielding and failure events is to calculate the
relevant strength factor for a prescribed bottomhole pressure or mud density. In a general sense,
strength factor is defined as the ratio of rock strength to the stress(es) that drive yielding or
failure. As such, a strength factor less than 1.0 indicates that yielding or failure has initiated.
Various types of strength factors are used in STABView. Their specific definitions are given as
follows:
For compressive shear failure (both active and passive):
σ1f − σ 3
SF = (4.11.1)
σ1 − σ 3
where:
SF = strength factor
Note: Given that there is limited practical value in reporting strength factors that are excessively
large, the output generated by STABView is truncated at a maximum value of 10.
STABView 3.8 User’s Manual 77 November, 2008
The analysis type is the primary input required to define the type of problem the user is asking
the software to solve. Several graphs and options will only be available in a specific analysis
type and/or mode. Tables 5.1 and 5.2 provides a detailed overview of the graphs available for
each analysis type and mode. Table 5.3 briefly summarizes types of analysis available in
STABView.
This input allows the user to specify the type of analysis to be conducted, and whether the
objective is to calculate risk factors for specified wellbore pressures (or mud densities) as
opposed to calculating the wellbore pressure required to achieve specified risk parameter values.
A user may switch between analysis modes but when the user switches modes they will be
presented with a warning dialog. Figure 5.1 shows the Switch Analysis Mode warning dialog
that is displayed when a user switches from “Pressure” to “Risk” analysis mode. A similar dialog
is presented when the user switches from “Risk” to “Pressure” mode. A user should take note of
the warnings presented and proceed with caution if they choose to continue. The user should
also note that not all STABView graphs are available in both “Risk” and “Pressure” modes.
After switching analysis mode or type all open graphs will continue to function in the mode and
type they were originally opened in.
STABView 3.8 User’s Manual 78 November, 2008
Hole Collapse While Use this analysis type to calculate the collapse Controlled Pressure Drilling
Drilling risks associated with a specified wellbore Underbalanced Drilling
pressure (or mud density) or to calculate a Overbalanced Drilling
wellbore pressure (or mud density) that prevents
hole collapse.
Fracture Breakdown / Use this analysis type to calculate the fracture Overbalanced Drilling
Lost Circulation initiation or breakdown risks for a specified Completions
wellbore pressure or to calculate a wellbore Workovers
pressure that prevents fracture breakdown. This Stimulation
option can also be used to investigate the risk of
fracture re-opening or fault reactivation.
Hole Collapse and Use this analysis type to calculate the collapse Overbalance Drilling
Fracture Breakdown / and fracture breakdown risks associated with a Controlled Pressure Drilling
Lost Circulation specified wellbore pressure or to calculate a
wellbore pressure that prevents hole collapse and
fracture breakdown. This is where a “mud
weight window” can be generated in pressure
mode.
Sand Production / Open Use this analysis type to calculate the sanding Sand Production Prediction
Hole Stability Under risks associated with a specified drawdown Openhole completion
Drawdown pressure or to calculate a drawdown pressure that Cavity completions
prevents or mitigates sanding. Openhole, slotted
liner, screen and perforated completion options
available.
Figure 5.1: Switching Analysis Mode dialog shown when a user switches
from “Pressure” to “Risk” mode.
STABView 3.8 User’s Manual 79 November, 2008
Table 5.2: Available Graphs in STABView by Analysis Mode, Type and Model - Drilling Applications
Combined - Hole Collapse and Fracture Breakdown / Lost Circulation
Analysis Type
Hole Collapse While Drilling Fracture Breakdown / Lost Circulation
Model Type Borehole Collapse Models Fracture Models
2D Permafrost 2D Elastoplastic 3D Linear Elastic 3D Linear Elastic 3D Linear Elastic 3D Anisotropic 3D Slip Analysis Fracture Re- Fracture Mechanics
Hole with Mohr- with Mohr- Model for Tensile Model for the Elastic Properties for a Weak opening of a Weak Based Model
Enlargement Coulomb Failure Coulomb , Hoek Fracture Initiation of Passive Model for Tensile Discontinuity Discontinuity Tensile Fracture
Calculation
Criterion Brown or Lade Breakdown with or Shear Yielding Fracture Breakdown
Models Failure Criteria without Deep Fluid ("Borehole
Penetration Ballooning")
Pressure & Flow Borehole Collapse Borehole Collapse Pressure & Flow Pressure & Flow Pressure & Flow Pressure & Flow Pressure & Flow Pressure & Flow
Properties Options Options Properties Properties Properties Properties Properties Properties
Stress Depletion Pressure & Flow Pressure & Flow Fracture Calibration Fracture Calibration Fracture Calibration Fracture Calibration Fracture Calibration Fracture Calibration
Properties Properties
Osmotic Effects Collapse Calibration Collapse Calibration Stress Depletion Stress Depletion Stress Depletion Stress Depletion Stress Depletion Stress Depletion
Available
Model Thermal Effects Stress Depletion Stress Depletion Osmotic Effects Osmotic Effects
Options Temperature Osmotic Effects Osmotic Effects Thermal Effects
Permafrost Thermal Effects Thermal Effects Temperature
Temperature Temperature
Strains &
Deformations
Depth Profile - Collapse Risk//Pressure Depth Profile - with Strength Factor/Fracture Breakdown Pressure
Radial Distance - Pressure Well Profile - Plan with Strength Factor/Fracture Breakdown Pressure
Radial Distance - Stress and Pressure Well Profile - VS Plane with Strength Factor/Fracture Breakdown Pressure
Cross‐section
Strength Factor
Well Profile 3D with Pressure Window
Permafrost - Hole Azimuth - Collapse Azimuth - Fracture Breakdown Pitch - Fracture Dip Angle - Slip Polar Plot - Re- Inclination - Fracture
Enlargement Risk/Pressure Risk/Pressure Breakdown Pressure Risk/Pressure opening Breakdown Pressure
Risk/Pressure
Table 5.3: Available Graphs in STABView by Analysis Mode, Type and Model - Sand
Production Application
Analysis Type Sand Production / Open Hole Stability Under Drawdown
Model Type Borehole or Perforation Collapse Models
Calculation 2D Elastoplastic with Mohr-Coulomb Failure 3D Linear Elastic with Mohr-Coulomb, Hoek Brown
Models Criterion or Lade Failure Criterion
Borehole Collapse Options Borehole or Perforation Collapse Options
Pressure & Flow Properties Pressure & Flow Properties
Collapse Calibration Collapse Calibration
Available Stress Depletion Stress Depletion
Model Options
Completion Type Completion Type
Strains & Deformations Gas Properties
Gas Properties
Depth Profile - Collapse Pressure / Risk**
Depletion - Critical Drawdown
Pressure - Collapse Risk
Pressure - Deformation Risk
Tornado/Spider Plot - Collapse Pressure / Risk
Standard Radial Distance - Pressure
Output* Radial Distance - Stress and Pressure
Well Profile - 3D with Collapse Pressure / Risk**
Well Profile - Plan with Collapse Pressure / Risk**
Well Profile - VS Plane with Collapse Pressure / Risk**
Cross-section Yielded Zone
Cross-section Strength Factor
Azimuth - Collapse Pressure / Risk
Inclination - Collapse Pressure Risk
Additional
Output for Polar Plot - Collapse Pressure / Risk
Boreholes and Radial Distance - Normal Stresses & Pressure
Perforations Radial Distance - Principal Stresses & Pressure
Radial Distance - Shear Stresses & Pressure
Model Properties dialog showing the Analysis Type and Case Properties tabs.
This input form allows the user to specify background information for the current analysis. The
data entered on this form are not used by the STABView calculation algorithms, but they are
listed in the text-based output reports generated by the program.
The selection of “Single Depth Analysis” or “Multiple Depth Analysis” for the “Data Entry
Mode” option affects the configuration of subsequent “Model Properties” input forms. The
remainder of this section assumes that the user has selected the “Multiple Depth Analysis”
option, and describes the procedures required to operate in this mode. In cases where the “Single
Depth Analysis” mode is selected, a limited number of model properties are required before
proceeding to the “Input Rock Mechanical and Pressure Data” form.
The selection of an “Onshore” or “Offshore” well location will affect how vertical stress
gradients are calculated (due to the absence or presence of a column of sea-water).
For rapid analyses, the user may specify a simple well trajectory in the form a vertical well, an
inclined well with constant trajectory from surface to total depth, or a horizontal well segment
occurring at a specified vertical depth such as a well branch. For more rigorous simulations of
openhole fracturing risks or more complex well designs, the user can elect to build or specify a
directional well survey. Required inputs for this survey are measured depth, well inclination and
well azimuth. The user can enter these values manually, or paste data generated in or provided
by another application into the input form.
STABView 3.8 User’s Manual 82 November, 2008
From the three columns of well survey input data STABView calculates by the Minimum
Curvature Method the true vertical depth (TVD), north-south and east-west offset distances and
dogleg severity (DLS).
Well Data and Well Survey tabs from Model Properties Dialog.
The user has the option to specify a uniform borehole diameter for the entire well, or to enter
values in a table that lists the borehole diameters used for different depth intervals of the well.
The borehole diameter data for selected analysis depths are listed in the text-based output reports
generated by STABView.
5.1.7 Tubulars
This input data form allows the user to specify the internal and external dimensions and
measured depth interval for casing strings, liners, screens or production tubulars in the well being
analyzed. These data are not currently used by STABView's calculation algorithms, but they are
In this dialog the user first enters names for the stratigraphic
units to be analyzed for a given well with STABView. The
lithology of each unit can be manually entered, or selected
from the drop-down list of pre-defined lithologies. Example
patterns for various lithologies are shown in Figure 5.2. By
clicking on the colored button beside the lithology descriptor,
the user selects a color and pattern for each stratigraphic unit.
Figure 5.2: Example of
The top and bottom depth for each unit can be entered as a lithology patterns and colors
measured depth or a true vertical depth. It is NOT required that available in STABView.
stratigraphic units be specified from the ground surface level
(or sea-floor for offshore wells) to total depth, nor is it required that the bottom depth of a given
stratigraphic unit is equal to the top depth of the underlying unit. STABView will only calculate
borehole collapse and fracture pressures or risks in well intervals for which stratigraphic units
have been defined.
The user is required to enter in-situ stresses and formation pressures for each stratigraphic unit.
These values can be specified as average gradients or as absolute values. When the gradient
option is selected, the user may choose to enter a constant gradient for each stratigraphic unit by
selecting the “Use Identical Gradient Values at Top and Bottom of Each Unit” check box. If this
STABView 3.8 User’s Manual 84 November, 2008
check box is not selected, or if the “Absolute Values” option is selected, the user must enter
stress and pressure values at the top and at the bottom of each unit.
When a formation pressure or in-situ stress is entered as an average gradient, the absolute value
of this parameter is calculated assuming that the gradient specified is with respect to vertical
depth from the kelly bushing (KB) or rotary table. For example, the vertical stress at a given
vertical depth (TVD) is calculated from a vertical stress gradient as follows:
σv = σv gradient x TVD (KB) (5.1)
Note that the same procedure is used in STABView to calculate vertical and horizontal stresses
or formation pressures for offshore and onshore wells. While this calculation will in fact
introduce an error for offshore wells, that depends upon the water depth and the depth of the
target interval below the seabed, the procedure is consistent with the typical practice used by
most service and operating companies. If the user wants a more accurate profile of stresses and
gradients, for example, a deep water case where the initial casing string is the concern, then the
user should enter absolute stress and pressure data into STABView, not gradients.
The user is required to specify values for all of the properties listed for each stratigraphic unit.
The default Mohr-Coulomb shear strength parameters (peak cohesion and peak friction angle),
static elastic properties (Poisson’s ratio and Young’s modulus), tensile strength and Biot’s
coefficient are well-known rock mechanical properties. Depending on the model or model
options a user has selected other rock mechanical properties may need to be entered on other
dialogs.
STABView 3.8 User’s Manual 85 November, 2008
The user is required to specify values for all of the properties listed for each stratigraphic unit.
The capillary threshold pressure is the overbalance pressure that must be exceeded before a non-
wetting drilling fluid (typically oil- or synthetic-based) will penetrate a water-wet formation.
Once the threshold pressure has been exceeded, it is assumed that mud pressure penetration
occurs.
The fluid penetration coefficient (εfp), which is specified for fracture or lost circulation problems,
and was defined in Section 4.6.1, is a measure of how well the pressure penetrates the near well
area. A fluid penetration coefficient of 1.0 indicates a fully-penetrating fluid, whereas a value of
0.0 indicates a non-penetrating fluid.
Mathematically, the filter cake or wall coating efficiency is equal to one minus the fluid
penetration coefficient.
The linear elastic algorithms used for 3D fracture breakdown and borehole collapse calculations
at the borehole wall in STABView are sensitive to the pore pressure within the formation right at
the wall. Pore pressures are actually calculated as a function of radius for steady-state inflow or
outflow conditions. The latter pressures are used for elastic and elastoplastic yielded zone size
calculations.
For any given stratigraphic unit, the user may activate any of a number of modeling options. The
first options pertain to the selection of the type of fracture breakdown or borehole collapse
model. These model types are described in Section 4. The additional input parameters required
for these model types are displayed on the "Model Options" tab once the appropriate model type
has been selected. For certain model types, input parameters specified on the “Rock and Fluid
Properties” input tab are implicitly modified or not used.
The availability of additional model options depends on the selected model type. Once the
model type is selected, only the relevant options are displayed. Tables 5.2 and 5.3 summarizes
the available model options available for each model type. A brief description of all of these
options follows.
Borehole Collapse Options: A number of additional options are available when using the
borehole collapse models. There is an option to specify the amount of maximum tolerable
yielding around the borehole. By default the collapse models used in STABView predict the
initiation of borehole yielding, where in reality some amount of yielding may generally be
tolerated. Other options are available for the calculation of effective stress and analyzing the
effects of a plane of weakness on borehole collapse risks.
Calibration: If the fracture breakdown or borehole collapse pressure for a given depth and well
trajectory is known from previous experience, two options are provided to calibrate
STABView 3.8 User’s Manual 87 November, 2008
STABView’s linear elastic models to provide output that is consistent with these observations.
For the Linear Multiplier option, the ratio of the known pressure to the theoretical pressure
calculated with the selected model is determined. Subsequently, all of the calculated output is
multiplied by this ratio. For the Linear Offset option, the difference between the known pressure
and the theoretical pressure calculated with the selected model is determined. Subsequently, this
difference of offset is subtracted from all of the calculated output.
Stress Depletion or Re-pressurization: Formation pressure changes in many reservoirs are
known to alter the total in-situ stress magnitudes within the reservoir. The user may select
whether to estimate these changes using a simple uniaxial strain, poro-elastic model suitable for
passive basins, or using simple plane-strain solutions suitable for normal or thrust-faulted basins.
Formation re-pressurization due to injection may also be accounted for in this option. See Addis
(1997) for more technical details on these options.
Osmotic Pressure Effects: Many shale, mudstone and claystone formations can act as partially
effective osmotic (i.e., semi-permeable) membranes. For such shales, near-well pore pressure
changes can occur as a result of water flows driven by chemical potential gradients which exist
when the chemical activity of water in the borehole is different from the chemical activity of the
formation water. For these conditions, STABView can calculate the osmotically-induced pore
pressure change (which may be positive or negative) at the borehole wall, which is additive to
the pore pressure change caused by fluid leakoff. For a complete review of the subject of
physico-chemical interaction in shales, see the Gas Research Institute report prepared by
Advanced Geotechnology (Hawkes et al., 2000).
Permafrost: The 2D Permafrost Hole Enlargement model requires several additional input
parameters that describe the drilling mud properties, rheological model and hydraulic
parameters. A more thorough discussion of the input parameters and hole enlargement model in
permafrost used in STABView can be found in Appendix D.
Strains and Deformations: If the 2D Elastoplastic model has been selected the user has the
option of analyzing the affects of borehole deformation or liner loading. Further details about
these calculations can be found in Section 4.3.4 Elasto-plastic Modelling of Borehole
Deformations and Liner Loading.
Thermal Effects: The introduction of a borehole fluid with a temperature that is different from
the initial formation temperature results in the generation of thermal stresses in the formation.
When this option is selected, the effect of these additional stresses on borehole collapse and
fracture breakdown pressures is included in the calculated output.
Temperature: When either of the “Osmotic Effects” or “Thermal Effects” advanced options are
selected, wellbore and formation temperatures are required input parameters.
STABView 3.8 User’s Manual 88 November, 2008
Edit the rock mechanical and pressure data Switch analysis mode from pressure to risk
Create a new graph or report Switch analysis mode from risk to pressure
Switch to tornado plot from spider plot Switch to spider plot from tornado plot
Set pressure or risk parameter type Tile the open graph windows
“Copy to Clipboard” places a copy of the currently selected graphic or text into the clipboard, so
it can be pasted into another application.
“Edit the Model Properties” opens the “Model Properties” input data form, allowing the user
to change any of the entries. Once the “OK” button is clicked on this form, any existing output
data windows are re-calculated and refreshed. This capability can also be accessed using the
toolbar button for editing the model properties on the main toolbar.
“Edit the Rock Mechanical and Pressure Data” opens the “Input Rock Mechanical and
Pressure Data” input data form, allowing the user to change any of the existing data entries or
advanced options. Model output is re-calculated and refreshed when the "OK" button is clicked.
This capability can also be accessed using the toolbar button for editing the rock mechanical and
pressure data on the main toolbar.
“Create a New Graph or Report” allows the user to create a new output window that contains
an output graph or text-based report. This capability can also be accessed using the toolbar
button for creating a new graph or report on the main toolbar. See Section 5.4 Adding Graphical
or Text Based Output Windows for more details.
STABView 3.8 User’s Manual 90 November, 2008
“Run Case” allows a user to recalculate the case after changes have been made to model
properties and input data. This menu is enabled by default and allows the user to make several
modifications to case parameters without forcing recalculation. A user may disable the “Run”
menu from the “Tools – Set Run Case Preferences” menu option
“Run Graph” provides the user the option to just generate the output for selected graph. This
can be useful when a user has a number of graphs open and would like to make a change on just
one graph without running the entire case again.
“Define A Custom Units Set” opens a utility that lets the user assign SI Metric or US Oilfield
units for each data type used by the program.
“Set Sensitivity” allows the user to select a
sensitivity parameter or collapse model for
comparison on the current output graph. Graph
sensitivity properties can also be accessed using
the toolbar button for setting sensitivity on the
active graph window toolbar. After selecting the
menu or toolbar option a select sensitivity group
dialog will open as shown in Figure 5.5. This Figure 5.5 Sensitivity Group Selection
dialog allows the user to select from the various Dialog.
sensitivity options available for display on the
selected graph. Figure 5.6 shows the Input
Sensitivity dialog for selecting general input
parameter sensitivities which allows the user to
specify input ranges as absolute values or as a
percentage of the base value.
The current version of STABView provides
support for showing general input parameter
sensitivities and collapse model comparisons.
Future versions of STABView will expand on
this capability and allow the user to add
additional parameters such as calibration, stress
depletion models, fracture models and wellbore
pressures.
Figure 5.6 Sensitivity Input Parameter
“Set Pressure or Risk Parameter Type” Dialog.
allows the user to specify whether the output
should be expressed as a bottomhole pressure, an
average pressure gradient, or an equivalent circulating density (ECD). This feature can also be
accessed using the toolbar button for setting pressure or risk parameter type on the active graph
window toolbar.
STABView 3.8 User’s Manual 91 November, 2008
“Set Scaling” allows the user to specify the maximum and minimum value for the pressure axis
on the current graph. For polar plots, no data will be plotted for areas where the data values fall
outside the specified range. This feature can also be accessed using the toolbar button for setting
scaling on the active graph window toolbar.
“Set Run Case Preferences” allows the user to disable or enable the “Run” button and menu
“Set Analysis Mode” allows the user to switch between analysis modes quickly. The user
should still be aware of the potential risks as outlined in section 5.1.2.
“Change Well Profile Settings” allows the user to set certain properties for all well profile
plots. The user can adjust the graphics quality and calculation speed of these plots.
“Export Datasets to File” allows the user
to export the raw data for the current graph
window to a text file.
“Export Graphic Image to File” allows
the user to export the current graph
window as an extended metafile (.emf) or
bitmap (.bmp) format graphical image.
“View Datasets in Notepad” dumps the
raw data for the current graph into the
Notepad text editor tool.
“Preferences” displays a dialog (Figure
5.7) that allows the user to customize
global application settings. A detailed
description of the available settings and
options is provided in Table 5.5.
Figure 5.7: Example of Preferences dialog box.
“Tile Output” and “Cascade Output” provide two different methods for organizing and
viewing the output windows. These features can also be accessed using the toolbar button for
tiling the open graphs and the toolbar button for cascading the open graph windows on the main
toolbar.
STABView 3.8 User’s Manual 92 November, 2008
General Preferences
Run Case Set case to run automatically after changes are Yes Yes
made on the Model or Input dialogs. No
Default Units Defines the unit set used when generating a new Metric (Modified SI) Metric
case but does not have any effect when opening an US Oilfield (Modified SI)
existing case.
Metric High Precision
US Oilfield-High Precision
User defined units
Well Profile Resolution Sets the resolution all graphs showing a well 1-4 1
profile (i.e., depth profile plot) will be displayed 1 = Best speed
at.
4 = Best graphics
Cross-Sections
Display tubulars and rock Display tubulars and rock deformations if Yes Yes
deformations? available. No
Display complete set of input Display an additional set of parameters on the Yes Yes
parameters? right side of a cross-section. No
Display breakout angle? Displays the breakout angle and draws the angle Yes Yes
on the borehole cross-section for the 3D elastic No
model. The 2D elastoplastic model does not
generate a breakout angle.
Polar Plots
Fracture Breakdown Pressure Use minimum in-situ stress as lower limit for Yes Yes
breakdown pressure? No
Color Display all polar plots in color or gray scale. Plot in color Plot in color
Plot in grayscale
Resolution Sets the resolution a polar plot will be displayed 1-4 1
at. 1 = Best speed
4 = Best graphics
Angles (Apply to wells and perforation tunnels)
Vertical Well Angle Tolerance Near Vertical Wells: For well inclinations greater 0-45° 15°
than this value, cross-sections of the wellbore are
referenced to the top side of the hole. (TOH).
Horizontal Well Angle Near Horizontal Wells: When the difference 0-45° 15°
Tolerance between the well inclination and a horizontal
trajectory is less than this value the top side of the
hole is labeled as the vertical stress direction.
Angular Tolerance for the 2D Maximum inclination and azimuth angle at which 0-45 25°
Elastoplastic Model a well trajectory can deviate from a principal stress
direction and still use the 2D elastoplastic model.
Wells with an inclination angle greater than this
value cannot use this model.
Horizontal Stress Ratio Maximum allowable ratio of the two horizontal in- 5-100% 5%
Tolerance for the 2D situ principal stresses , that will allow a user to use
Elastoplastic Model the 2D elastoplastic model, for a horizontal well
section oriented within the angular tolerance of a
horizontal principal stress
STABView 3.8 User’s Manual 93 November, 2008
This menu entry provides access to the help file, as well as basic information about the version of
STABView installed on the user’s computer.
STABView offers considerable flexibility in terms of the quantity and character of output data
displayed. The user can generate practically any number of output windows, each of which may
contain any of a number of graphical or text-based displays. Below is a list of the currently
available output windows in STABView. Note that the availability of these graphs depends on
the analysis mode/type and modeling options that have been selected. Appendix F provides a
detailed overview of the graphs available for each analysis type/mode and options.
STABView Output Graphs and Reports
General Output (all analysis types)
• Text Report – Multiple Depth Input Report
• Text Report – Multiple Depth Output Report
• Text Report – Single Depth Analysis Report
• Well Profile – 3D View
• Well Profile – Plan View
• Well Profile – Vertical Section (VS) Plane View
Hole Collapse While Drilling
General Output (both analysis modes)
• Cross-Section – Yielded Zone
• Cross-Section – Strength Factor (not 2D Elastoplastic Model)
• Pressure – Collapse Risk
• Pressure – Deformation Risk
• Radial Distance – Normal Stresses and Pore Pressure
• Radial Distance – Shear Stresses and Pore Pressure
STABView 3.8 User’s Manual 94 November, 2008
Risk Mode
All of the same graphs available for Hole Collapse While Drilling and Fracture Breakdown /
Lost Circulation
STABView 3.8 User’s Manual 96 November, 2008
When the user clicks on the depth profile properties button a graph properties dialog will open as
shown in Figure 5.11. The Axis tab on the graph properties dialog allows the user to change the
depth range, the depth type, change the depth and X-Axis scaling, and change the position of the
X-Axis. Any change the user makes will immediately be displayed on the relevant graph.
Clicking the Ok button will accept any changes the user has made and likewise clicking the
Cancel button will cancel any changes.
The Data tab of the graph properties dialog as shown in Figure 5.12 displays all the data series as
generated for the graph type and any sensitivities the user has added. A user may modify the
name, window bounds, visibility and color of any of the available data series.
STABView 3.8 User’s Manual 98 November, 2008
Figure 5.10: Example of a Depth Profile graph showing a Mud Weight Window.
STABView 3.8 User’s Manual 99 November, 2008
STABView provides the user with a wide selection of standard graph types for output as detailed
in Section 5.3. Figure 5.22 shows an example of two such graphs. Typical of the standard
graphs in STABView is the depth bar on the left side of the graph. This depth bar provides a
high level view of the stratigraphic units and their lithologies along the well survey. The arrow
indicator as seen in Figure 5.22 represents the depth at which the graph is rendered. At any time
the user may change the point of calculation and create a new graph by clicking above or below
the point indicated by the arrow. A user may also change the point of calculation by entering a
Depth value in Calc Depth text box and press the return key. The depth bar may be changed to
MD or TVD by clicking the appropriate radio button. Hovering the mouse over the depth bar
will provide a small popup window with information about the current depth and the collapse
and/or fracture models being used in the zone.
The toolbar at the top of all standard graphs displays a collection of buttons that will allow the
user to customize and/or change the default graph. Table 5.6 summarizes the available icons
and their function.
STABView 3.8 User’s Manual 105 November, 2008
Allows the user to change the pressure or risk parameter being graphed.
Table 5.7 summarizes the STABView demonstration cases which are normally shipped with the
program. These example files should be examined and re-run to evaluate and test the capabilities
of the software. For some STABView licensees a different set of demonstration cases have been
installed, depending upon the user requirements. Please contact us at software@advgeotech.com
if you have any problems with these cases.
STABView 3.8 User’s Manual 106 November, 2008
Table 5.7: STABView well planning sample cases typically included with STABView as demonstration files
No. STABView Type of No. Zones Analysis Analysis Well Options/Features Selected Comments
Sample Case Model Type Mode Trajectory
Name
1 3D Drilling 3D - elastic single Hole Pressure inclined Effect of a plane of weakness on Evaluation of the effects of a fissile shale on wellbore stability
Collapse Case collapse wellbore collapse. Collapse with a calibrated 3D elastic model.
with Weak while calibration.
Bedding drilling
2 3D Hole 3D - elastic multiple Hole Risk inclined Failure criteria sensitivity. Depth Pseudo-yielded zone calculation for an inclined well based on
Enlargement collapse profile. the extent of the failed zone determined from the default 3D
Model while linear elastic model. It is interesting to compare the results
drilling obtained when Modified Lade and Hoek-Brown failure criteria
are used.
3 Horizontal Well 2D - elasto- single Sand Risk horizontal Steady - state pore pressure gradient Horizontal sand production case to decide if a slotted liner,
2D Sand plastic production effects. Depletion induced stresses. screen or gravel pack are required. Yielded zone size and
Production Risk Barefoot openhole completion. deformations predicted for various sensitivities, e.g. dilation
Calculation of the rubble fill angle, peak cohesion, residual cohesion.
percentage.
4 Horizontal Well 2D - elasto- multiple Open hole Risk horizontal Liner, screen or pre-packed screen Horizontal sand production case to assess if hole closure results
with Liner plastic stability/ completion. Calculation of in contact between the formation and a stiff slotted liner or
sand formation stresses on a stiff screen screen used to complete this well.
production or liner due to hole closure.
Calculation of plastic strain at the
wellbore wall.
5 Extended Reach 3D - elastic multiple Hole Pressure build plus Collapse calibration. Shale-mud Determine safe range of EMW and ECD to prevent borehole
Well with collapse and horizontal interaction and osmotic pressure collapse and fracture breakdown in the build section of the well
Horizontal Leg fracture (survey) effects on wellbore collapse. Mud and within a depleted reservoir section of a horizontal well.
breakdown weight window.
6 Inclined Well 3D - elastic multiple Fracture Pressure survey Steady - state thermal effects on lost Effect of cool fluid injection on fracture breakdown pressure
with Cold Fluid breakdown circulation and fracturing risks. for an open hole completion in a directional well.
Injection
7 Injection Well 3D - elastic single Lost Pressure vertical 3D slip risk analysis on a weak Analysis to determine if elevated pore pressures during water
Slip Pressure Circulation discontinuity, e.g., natural fractures, injection will cause shear failure on weak natural fractures or
Analysis faults or bedding. bedding planes.
8 Permafrost Hole 2D - single Hole Risk vertical 2D permafrost hole enlargement Example of a vertical well drilled through permafrost
Enlargement permafrost collapse model.
erosion while
model drilling
9 Permafrost Hole 2D - multiple Hole Risk vertical 2D permafrost hole enlargement Example of a vertical well drilled through rich, weak erodable
Enlargement permafrost collapse model. permafrost intervals in the Mackenzie Delta, N.W.T., Canada.
Multi - Zone erosion while
model drilling
10 Underbalanced 3D – elastic multiple Hole Pressure survey Steady - state pore pressure gradient Multi-zone wellbore stability analysis for an offshore, extended
Drilling in Fissile and 2D collapse effects. Depletion induced stresses. reach underbalanced well penetrating a fissile shale at angle.
Shale elasto-plastic while Effects of a plane of weakness on Case takes several minutes to run 2D elasto-plastic and 3D
drilling wellbore collapse. Collapse elastic models compared for high angle penetrations.
calibration. Mud weight window.
STABView 3.8 User’s Manual 107 November, 2008
Table 5.7: STABView well planning sample cases typically included with STABView as demonstration files (Cont.)
No. STABView Type of No. Zones Analysis Analysis Well Options/Features Selected Comments
Sample Case Model Type Mode Trajectory
Name
11 Underbalanced 2D - elasto- multiple Hole Risk survey Steady - state pore pressure Multi-zone wellbore stability analysis for an offshore, extended
Drilling plastic collapse gradient effects. Stress depletion reach underbalanced well. Bottomhole pressure profile used as
while model for changing reservoir input.
drilling pressures. Yielded zone
predictions.
12 Vertical Gas Well 2D - elasto- multiple Open hole Risk vertical Barefoot openhole completion. Vertical gas well with barefoot open hole completion. The user
Open hole plastic stability Effects of gas properties, has input a gas flow rate, from which STABView calculates the
Completion under temperature and non-Darcy flow bottomhole pressure for its use in a borehole collapse model.
drawdown rates on openhole stability.
13 Vertical Oil Well 2D - elasto- multiple Open hole Risk vertical Perforated liner completion with Vertical well with perforated completion. User has specified a
plastic stability cylindrical perforation cavities. liquid production rate, from which STABView calculates the
under Calculation of the perforation bottomhole pressure based on the steady-state radial flow
drawdown stability for the most critical equation for an incompressible fluid.
perforation orientation.
14 Vertical Well with 2D - elasto- multiple Hole Risk vertical Steady - state pore pressure 2D elasto-plastic analysis of the yielded zone size for a vertical
Collapse Risk plastic collapse gradient effects (two-zone mobility well drilled overbalanced through multiple zones.
while model).
drilling
15 Inclined Well 3D - elastic single Sand Risk/ inclined Effect of well trajectory and Calibrated 3D analysis model for cylindrical and hemi-spherical
with Perforations production Pressure perforation orientation on fracture perforation cavities. Assess the value of oriented perforating.
breakdown pressures and sand
production risk
16 Horizontal Well 2D - elasto- multiple Sand Pressure survey Effect of reservoir depletion on the 2D elastoplastic analysis model used for examining the effects
with Effects of plastic Production sand production critical bottom of reservoir depletion for a range of critical bottom hole
Depletion 2D hole pressure or drawdown pressures.
17 Complex Multi- 2D multiple Hole Pressure survey Demonstrates using a combination A complex drilling case showing multiple zones and models
Zone Offshore elastoplastic & collapse and of 2D/3D analysis models in a used to calculate wellbore stability collapse and fracture
Well 3D elastic fracture drilling problem in order to pressures for an offshore well.
breakdown determine a safe operating mud
weight window.
STABView 3.8 User’s Manual 108 November, 2008
STABView 3.8 User’s Manual 109 November, 2008
Q: How does STABView differ from other commercial wellbore stability programs
available on the market?
A: Drawing on a combined total of more than 30 years of experience in solving industry
geomechanics problems, Advanced Geotechnology has developed several unique approaches to
wellbore stability, lost circulation and sand production risk assessment. STABView was
primarily designed for technical drilling, completion and production engineers. We have
evaluated and tested many approaches to integrating the theoretical and practical aspects of well
construction, drilling fluids engineering, production optimization and oilfield operations.
STABView combines the best of these approaches into a user-friendly interface with rapid
analysis capabilities.
The following are some of the unique features and capabilities in STABView:
major program additions will be available to current licensees for an upgrade price. Companies
purchasing the support and maintenance option will also receive a discount on in-house or
external short courses on borehole stability and sand production offered by AG. If you are
interested in more details, contact STABView Technical Support at software@advgeotech.com
to receive a current list of software purchase, maintenance and training options.
Q: How does STABView evaluate the risk of borehole collapse while drilling?
A: STABView has two primary classes of models for evaluating borehole collapse risk. Both of
these model classes assess the mechanical integrity of the rock around the circumference of the
borehole. STABView's 2D model evaluates the size of the zone of yielded rock around the
borehole using the theory of elastoplasticity, including pore pressure effects during over- or
underbalanced drilling. This model is useful for vertical wells, or for horizontal wells that are
approximately parallel to either the maximum or minimum horizontal in-situ stress. STABView's
3D elastic model evaluates the bottomhole pressure, static mud density or equivalent circulating
density (ECD) at which shear yielding or tensile hydraulic fracturing initiates on the borehole
wall using linear elastic theory. When calibrated based on drilling experience for offset wells,
this model has been shown to reliably predict the bottomhole pressure mud density or ECD at
which catastrophic hole collapse will occur.
Q: How can STABView help determine drilling parameters that will result in a more stable
well?
A: STABView can evaluate the sensitivity of borehole instability risks to a large number of
controlling parameters. The most notable of these parameters are the drilling mud's equivalent
mud weight or circulating density, or its bottomhole pressure, the effectiveness of filter cake or
wall coating additives for preventing mud pressure penetration into the formation, and the
trajectory of the wellbore relative to the principal in-situ stresses and weak bedding or fracture
planes. Other controllable factors which can be assessed with STABView include capillary
threshold pressure effects for oil-based or pseudo oil-based drilling muds, the effects of high
salinity muds on shale inhibition, and the consequences of sudden decreases in bottomhole
pressures due to swabbing effects. Casing setting depths may also be optimized with
STABView using the software’s “Mud Weight Window” graph.
Q: Can I use STABView to assess borehole collapse risks for drilling wells underbalanced?
A: Yes. Both the 2D elasto-plastic and 3D elastic models in STABView are capable of assessing
borehole collapse risks for conditions where the bottomhole pressure is less than the pore
pressure of the formation being drilled. The effects of steady-state drawdown pressures in the
near-borehole area on the stability of the rock are accounted for in STABView. It is also
possible to include the effects of an external or internal filter cake at the borehole wall which has
not been removed or destroyed. Transient pressures in the near-borehole area, which may be
relevant in low permeability rocks are not currently accounted for in STABView. It is however
possible to investigate the most severe transient case at initial drawdown or instantaneous
drawdown conditions, i.e., where there is a step-change from the borehole pressure to the
formation pressure, and there is no reduction in the formation pressure near the well.
STABView 3.8 User’s Manual 111 November, 2008
Q: Can I use STABView to assess wellbore instability risks for directional and horizontal
wells?
A: Yes. The 3D linear elastic model in STABView is capable of assessing wellbore instability
risks to hole collapse or lost circulation for all possible well trajectories. The 2D elasto-plastic
model is capable of assessing borehole instability risks for wells that are within a specified angle
of a principal in-situ stress orientation. As such, slightly deviated wells and horizontal wells that
are sub-parallel to the horizontal in-situ stresses can be simulated using the 2D model.
Q: What features are available for assessing slip or re-opening risks on weak planes or
discontinuities such as faults, bedding planes, natural fractures or cleats?
A: Two models are available for assessing these risks. One model predicts the bottomhole
injection pressure or equivalent circulating density at which shear failure will occur on weak
planes. The second model predicts the pressure at which a discontinuity will re-open, assuming
no cohesion remains on the plane. These models are useful for identifying the maximum
injection pressure that can be used for waste liquids or solids disposal in the subsurface,
waterflooding, sand or slop waste injection or gas storage wells.
STABView 3.8 User’s Manual 113 November, 2008
7. SELECTED BIBLIOGRAPHY
7.1 Cited References
Aadnoy, B.S., Stresses Around Horizontal Boreholes Drilled in Sedimentary Rocks, SPE
paper 17119, unsolicited manuscript submitted July, 1987.
Aadnoy, B.S. and Chenevert, M.E., Stability of Highly Inclined Boreholes, SPE Drilling
Engineering, December, 1987, pp. 364-374.
Aadnoy, B.S., Stresses Around Horizontal Boreholes Drilled in Sedimentary Rocks, Journal
of Petroleum Science and Engineering, Vol. 2, pp. 349-360, 1989.
Addis, M.A., The Stress-Depletion Response of Reservoirs, SPE 38720, presented at the SPE
Annual Technical Conference and Exhibition, San Antonio, Texas, October 5-8, 1997.
Advanced Geotechnology Inc., ROCKSBank Rock Mechanical and Petrophysical Properties
Database, User’s Manual, Version 2.0, May 2006.
Bear, J., Dynamics of Fluids in Porous Media, Dover Publications Inc., New York, 1972,
764 p.
Bradley, W.B., Failure of Inclined Boreholes, Journal of Energy Resources Technology /
Transactions of the ASME, 1979, pp. 232-239.
Bratli, R.K. and Risnes, R., Stability and Failure of Sand Arches, SPE Journal, April, 1981,
pp. 236-248.
Butler, R., Horizontal Wells for the Recovery of Oil, Gas and Bitumen: Monograph Number
2, CIM Petroleum Society, Calgary, 1994, 228 p.
Detournay, E. and Fairhurst, C., Two-dimensional Elasto-plastic Analysis of a Long,
Cylindrical Cavity Under Non-hydrostatic Loading, International Journal of Rock Mechanics,
Mining Science, and Geomechanics Abstracts, Vol. 24, No. 4, 1987, pp. 197-211.
Detournay, E. and St. John, C.M., Design Charts for a Deep Circular Tunnel Under Non-
uniform Loading, Rock Mechanics and Rock Engineering, 21, 1988, pp. 119-137.
Dranchuk, P.M. and Abou-Kassem, J.H., Calculation of Z-Factors for Natural Gases Using
Equations of State, Journal of Canadian Petroleum Technology, vol. 14, July-Sept. 1975, pp.
34-36.
Drescher, A., Analytical Methods in Bin-Load Analysis, Elsevier, 1991.
Ewy, R.T., Wellbore Stability Predictions Using a Modified Lade Criterion, SPE Paper
47251, Presented at SPE/ISRM Eurock’98, Trondheim, Norway, July 8-10, 1998.
Fjaer, E., Holt, R., Horsrud, P., Raaen, A. and Risnes, R., Petroleum Related Rock Mechanics,
Elsevier, New York, 1992, 337 p.
STABView 3.8 User’s Manual 114 November, 2008
Gerald, C.G. and Wheatley, P.O., Applied Numerical Analysis, Addison Wesley, Don Mills,
Ontario, 1989, 679 p.
Hawkes, C.D., McLellan, P.J., Ruan, C.G. and Maurer, W.C., Wellbore Instability in Shales: A
Review of Fundamental Principles, Physico-Chemical Mechanisms in Mud-Shale
Interaction and GRI-Funded Research, CD-ROM, Part 1: Wellbore Instability in Shales: A
Review of Fundamental Principles and GRI-Funded Research, Part 2: Review of Physico-
Chemical Mechanisms in Mud-Shale Interaction and their Effects on Wellbore Stability, Gas
Research Institute Report No. GRI-99/0025.3, December, 2000.
Hoek, E., Kaiser, P.K., and Bawden, W.F., Support of Underground Excavations in Hard
Rock, Balkema, 1995, 215 p.
Hoek, E. and Brown, E.T., Underground Excavations in Rock, Institute of Mining and
Metallurgy, London, 1980, 527 p.
Jaeger, J.C. and Cook, N.G.W, Fundamentals of Rock Mechanics, Chapman and Hall, London,
1979, 593 p.
Koning, E.J.L., Waterflooding under Fracturing Conditions, Ph.D. thesis, Delft University of
Technology, The Netherlands, 1988.
Kutasov, I.M. and Caruthers, R.M., Hole Enlargement Control During the Arctic Drilling,
SPE 17442, presented at the SPE California Regional Meeting, Long Beach, California, 1988.
Kutasov, I.M., Applied Geothermics for Petroleum Engineers, Elsevier, 1999.
Lee, A.L., Gonzalez, M.H. and Eakin, B.E., The Viscosity of Natural Gases, Journal of
Petroleum Technology, August 1966, pp. 997-1000.
Londono, F.E., Archer, R.A. and Blasingame, T.A., Simplified Correlations for Hydrocarbon
Gas Viscosity and Gas Density – Validation and Correlation of Behavior Using a Large-
Scale Database, SPE 75721, Presented at the SPE Gas Technology Symposium, Calgary, AB,
April 30 – May 2, 2002.
Mody, F.K. and Hale, A.H., Borehole-Stability Model to Couple the Mechanics and
Chemistry of Drilling-Fluid/Shale Interaction, Journal of Petroleum Technology, November,
1993, pp.1093-1101.
Morita, N., Fuh, G.F. and Black, A.D., Borehole Breakdown Pressure with Drilling Fluids -
II. Semi-analytical Solution to Predict Borehole Breakdown Pressure, International Journal
of Rock Mechanics, Mining Sciences, and Geomechanics Abstracts, Vol. 33, No. 1, 1996,
pp. 53-69.
Risnes, R., Bratli, R.K. and Horsrud, P., Sand Stresses Around a Wellbore, SPE Journal,
December 1982, pp. 883-898.
Soltanzadeh, H. and Hawkes, C.D. Semi-Analytical Models for Stress Change and Fault
Reactivation Induced by Reservoir Production and Injection, manuscript submitted to the
Journal of Petroleum Science and Engineering, publication anticipated in 2007.
STABView 3.8 User’s Manual 115 November, 2008
Sutton, R.P., Fundamental PVT Calculations for Associated and Gas/Condensate Natural
Gas Systems, SPE 97099, Presented at the SPE Annual Technical Conference and Exhibition,
Dallas, TX, Oct. 9-12, 2005.
Wang, Y. and Dusseault, M.B., Borehole Yield and Hydraulic Fracture Initiation in Poorly
Consolidated Rock Strata - Part I. Impermeable Media, International Journal of Rock
Mechanics, Mining Sciences, and Geomechanics Abstracts, Vol. 28, No. 4, 1991, pp. 235-246.
Wang, Y. and Dusseault, M.B., Borehole Yield and Hydraulic Fracture Initiation in Poorly
Consolidated Rock Strata - Part II. Permeable Media, International Journal of Rock
Mechanics, Mining Sciences, and Geomechanics Abstracts, Vol. 28, No. 4, 1991, pp. 247-260.
McLellan, P.J., Hawkes, C.D., Comparison of Methods for Modeling Borehole Instability
and Its Consequences, SPE 78178. Presented at the SPE 2002 Oil Rock Conference, Irvin,
Texas, October 20-22, 2002.
McLellan, P.J., Hawkes, C.D., Smith, S., Coupled Modeling of Borehole Instability and
Multiphase Flow for Underbalanced Drilling, SPE 74447 presented at the IADC
Underbalanced Technology Conference, Aberdeen, Scotland, Nov. 27-28, 2001 and at the
SPE/IADC Drilling Conference, Dallas, Texas, Feb 26-28, 2002.
McLellan, P.J., and Hawkes, C.D., Practical Risk Assessment Techniques for Evaluating
Borehole Instability, Presented at the CADE/CAODC Drilling Conference, Calgary, Alberta,
October 23-24, 2001.
McLellan, P.J., and Hawkes, C.D., Sand Production and Control in Horizontal Wells for Gas
Storage Reservoirs, SPE 65510, presented at the 2000 SPE/CIM International Conference on
Horizontal Well Technology, Calgary, Alberta, November 6-8, 2000.
Hawkes, C.D., McLellan, P.J., Ruan, C.G. and Maurer, W.C., Wellbore Instability in Shales: A
Review of Fundamental Principles, Physico-Chemical Mechanisms in Mud-Shale
Interaction and GRI-Funded Research, CD-ROM, Part 1: Wellbore Instability in Shales: A
Review of Fundamental Principles and GRI-Funded Research, Part 2: Review of Physico-
Chemical Mechanisms in Mud-Shale Interaction and their Effects on Wellbore Stability, Gas
Research Institute Report No. GRI-99/0025.3, December, 2000.
Hawkes, C.D. and McLellan, P.J., A New Model for Predicting Time-dependent Failure of
Shales: Theory and Application, Journal of Canadian Petroleum Technology, Vol. 38, No. 12,
49-55, December 1999.
McLellan, P.J., Hawkes, C.D., and Yuan, Y., Minimizing Borehole Instability Risks in Build
Sections Through Shales, Presented at the 7th SPE/CIM One Day Conference on Horizontal
Well Technology, Calgary, Alberta, November 3, 1999.
Hawkes, C. and McLellan, P.J. Borehole Stability Analysis for Underbalanced Drilling, Paper
99-07, Presented at the Joint CSPG/Petroleum Society Conference, Calgary, Alberta, June 14-18,
1999.
STABView 3.8 User’s Manual 116 November, 2008
McLellan, P.J. and Hawkes, C.D., User-Friendly Borehole Stability Software for Designing
Horizontal and Deviated Wells, Paper 99-101, Presented at the CADE/CAODC Spring Drilling
Conference, Calgary, Alberta, April 7-8, 1999.
McLellan, P.J. and Hawkes, C.D., Applications of Probabilistic Techniques for Assessing
Sand Production an Borehole Instability Risks, SPE/ISRM Paper 47334, Presented at
SPE/ISRM Eurock ’98, Trondheim, Norway, July 8-10, 1998.
McLellan, P.J., Assessing the Risk of Wellbore Instability in Inclined and Horizontal Wells,
Journal of Canadian Petroleum Technology, Vol. 35, No. 5, 21-32, May 1996.
Hawkes, C. and McLellan, P.J., Transient Poro-elasto-plastic Modeling of Yielded Zone
Enlargement Around an Unstable Wellbore, Presented at the 2nd North American Rock
Mechanics Symposium - NARMS ‘96, Montreal, Quebec, June 1996.
McLellan, P.J. and Cormier, K., Borehole Stability in Dipping Fissile Shales, Northeastern
British Columbia, SPE 35634, SPE Gas Technology Symposium, Calgary, Alberta, April, 1996.
Hawkes, C.D., Effects of Pore Pressure Penetration on the Extent of the Yielded Zone
Around a Wellbore, Ph.D. thesis, University of New Brunswick, July 1996.
Warren, B.W., McLellan, P.J., and Pratt, K., Wellbore Stability and Drilling Fluids Design for
Horizontal Wells in Peace River Oil Sands, Alberta, SPE/IADC 29426, Amsterdam,
Netherlands, Feb. 1995.
McLellan, P.J., Wang, Y., Predicting the Effects of Pore Pressure Penetration on the Extent
of Wellbore Instability: Application of a Versatile Poro-Elasto-plastic Model, SPE/ISRM
28053, EUROCK 94 - Rock Mechanics in Petroleum Engineering, A Joint SPE/ISRM Meeting,
The Hague, Netherlands, August 1994.
Bell, J.S., Price, P.R., and McLellan, P.J., In-Situ Stress in the Western Canada Sedimentary
Basin, Chapter 29 in the Geological Atlas of the Western Canada Sedimentary Basin,
ISPG/ARC/CSPG, 1994.
STABView 3.8 User’s Manual A-1 November, 2008
APPENDIX A
BY DOWNLOADING OR INSTALLING THIS SOFTWARE, CLICKING “I ACCEPT”, AND/OR SIGNING THIS AGREEMENT, YOU
ACKNOWLEDGE THAT YOU HAVE READ ALL THE TERMS OF THIS AGREEMENT AND YOU ARE AGREEING TO BE BOUND
BY ALL OF THE TERMS OF THIS AGREEMENT. IF YOU DO NOT AGREE TO ALL OF THE TERMS DESCRIBED HERE, YOU
MUST NOT DOWNLOAD OR INSTALL THE SOFTWARE.
This is a legal agreement between the individual, corporation, or organization downloading or installing this software (the
"Licensee") and Advanced Geotechnology (A Division of Hycal Energy Research Laboratories Ltd.) and its successors ("AG"). This
License Agreement states the terms and conditions upon which AG is providing the STABView software together with all related
documentation and accompanying items including, but not limited to, executable programs, test versions of the software,
demonstration versions of the software, drivers, libraries, printed materials, graphics, images, animations and data files associated
with such programs (collectively, the "Software").
1. Grant of License. In consideration of downloading or installing the Software and the fees paid to AG hereunder, AG
hereby grants to the Licensee a non-exclusive, irrevocable and royalty-free license to use the Software subject to the terms of this
License Agreement. This License Agreement is the Licensee's proof of license to exercise the rights granted herein and must be
retained by the Licensee. As between the Licensee and AG (and, to the extent applicable, its licensors), AG retains all title to and
ownership of the Software and reserves all rights not expressly granted to the Licensee. AG grants to the Licensee the right to use
all or a portion of this Software provided that (a) the Software is not distributed for profit to other parties; (b) the Software may NOT
be modified, copied or decompiled; (c) all copyright notices, proprietary notices, legends and logos are maintained in the Software;
and (d) the Licensee, and its contractors, consultants and its end users agree to be bound by the terms of this License Agreement.
The specific options applicable to this license and the Licensee shall be set out in the Invoice or letter rendered by AG to the
Licensee which accompanies and forms a part of this License Agreement.
2. For Demonstration Use (Demo Option). If the Licensee has licensed the demonstration version of the Software, it may
be used only on a single computer by a single user at any time, for internal demonstration purposes only, for the specified limited
time set out in the attached letter. Notwithstanding anything to the contrary herein, the Licensee may not make any back-up copy of
the Software and may not use the Software for any commercial purpose, other than to determine the suitability of the Software for
the Licensee.
3. For Use on a Single Computer (Single User Option). If the Licensee has licensed the single user version of the
Software it may be used only on a single computer by a single user at any time. The Licensee may transfer the machine-readable
portion of the Software from one computer to another computer, provided that (a) the Software (including any portion or copy
thereof) is erased from the first computer, and (b) there is no possibility that the Software will be used on more than one computer at
a time. The Licensee agrees to report the location of all copies of the Software to AG at AG's request. This version of the Software
cannot be installed or used on an Application Service Provider (ASP) nor on a network where multiple users can access it. These
limitations also apply to a Licensee who possesses two or more stand alone PC licenses to the Software.
4. For Use on a Network (Network Option). If the Licensee has licensed the network version of the Software with a defined
number of seats it may only be used on a network by the same defined number of users at any time equal to the number of seats
licensed, as set out in the Invoice. The maximum authorized number of potential Software users for a network seat is 10. The
Licensee may only transfer the machine-readable portion of the Software from one installation to another installation provided that
(a) the Software (including any portion or copy thereof) is erased from the first network, and (b) there is no possibility that the
Software will be used on more than one network at a time. The Licensee agrees to report the location of all copies of the Software to
AG at AG's request. This version of the Software cannot be installed or used on an ASP.
5. For Use on an ASP (ASP Option). If the Licensee has licensed the ASP version of the Software, with a defined number
of seats, the Software may only be hosted with AG or an ASP approved by AG, and may only be accessed by the same defined
number of users at any time equal to the number of seats licensed, as set out in the Invoice. The Licensee may only transfer the
machine-readable portion of the Software from one host to another host provided that (a) the Software (including any portion or copy
thereof) is erased from the network of the first host, and (b) there is no possibility that the Software will be used with more than one
ASP at a time. The Licensee agrees to report the location of all copies of the Software to AG at AG's request.
6. Stand-Alone Basis. If the Licensee has licensed the stand-alone version of the Software, the Licensee may use the
Software only on a stand-alone basis, such that the Software and the functions it provides are accessible only to persons who are
physically present at the location of the computer on which the Software is installed. The Licensee may not allow the Software or its
functions to be accessed remotely, or transmit all or any portion of the Software through any network or communication line.
Licenses for multiple stand-alone users, local area network (LAN) use, or ASP use can be purchased from AG and are subject to a
separate or amended license agreement.
7. Copyright. The Software is owned by AG and is protected by Canadian and United States copyright laws and
international treaty provisions. The Licensee may not remove the copyright notice from any copy of the Software or any manuals or
written materials accompanying the Software.
8. No Merger or Integration. The Licensee may not merge any portion of the Software into, or integrate any portion of the
Software with any other program, except to the extent expressly permitted by the laws of the jurisdiction where the Licensee is
located or if permission to do so has been obtained from AG. Any portion of the Software merged into or integrated with another
program will continue to be subject to the terms and conditions of this License Agreement, and the Licensee must reproduce on the
merged or integrated portion all copyright and other proprietary rights notices included on the originals of the Software.
9. Attribution. The Licensee agrees to provide attribution to AG in any and all public-domain publications or presentations in
which results obtained from the Software are described, printed or presented.
STABView 3.8 User’s Manual A-3 November, 2008
10. Transfer of License. The Licensee may not sell, or transfer the license of the Software to a third party without AG's
written permission. The recipient must agree to the terms of this License Agreement.
11. Restrictions on Use.
a. Except to the extent expressly permitted by this License Agreement, or by the laws of the jurisdiction where the Licensee
acquired the Software, or for back-up purposes, the Licensee may not copy or modify the Software or make any derivative works
from the Software.
b. The Licensee may not sell, transfer, sublicense or distribute the Software or any part thereof to other parties;
c. The Licensee acknowledges that the Software contains trade secrets and other proprietary information of AG and its licensors.
Except to the extent expressly permitted by this License Agreement or by the laws of the jurisdiction where the Licensee is located,
the Licensee may not decompile, disassemble or otherwise reverse engineer the Software, or engage in any other activities to
obtain underlying information that is not visible to the user in connection with normal use of the Software. In particular, the Licensee
agrees, unless otherwise stipulated in writing from AG, not for any purpose to (i) transmit the Software or display the Software's
object code on any computer screen or to make any hard copy memory dumps of the Software's object code; and (ii) transmit or
display the Software to any of the following direct software competitors of AG: Advantek International, Geomechanics International,
Inc., International Technology Consultants, GeoScience Limited, Schlumberger, Halliburton, Baker Hughes, Maurer Technology,
and any other firms who develop or market a competing product to this Software.
12. Request for Information. If the Licensee believes it requires information related to the inter-operability of the Software
with other programs, the Licensee shall not decompile or disassemble the Software to obtain such information, and the Licensee
agrees to request such information from AG. Upon receiving such a request, AG shall determine, in its sole judgment, whether the
Licensee requires such information for a legitimate purpose and, if so, AG will provide such information to the Licensee within a
reasonable time and on reasonable conditions. In any event, the Licensee will notify AG of any information derived from reverse
engineering or such other activities, and the results thereof will constitute the confidential information of AG that may be used only in
connection with the Software.
13. Maintenance and Support. If the Licensee has purchased the maintenance and support option from AG as set out in the
Invoice, then AG will provide the maintenance and support services set out in Schedule “A”. Except to the extent specified in
Schedule “A” (if applicable), AG shall not be obliged to provide the Licensee with major new software releases, updates,
modifications, or enhancements to the Software.
14. Verification Rights. AG reserves the right, with reasonable notice and during regular business hours, to audit the
facilities, systems and records of the Licensee to audit and verify compliance with the terms of this license. In the event any such
audit or review discloses any breach of the terms of the License, AG may, at its option, terminate the license, and/or seek any other
remedy available under this Agreement or at law.
15. Termination. The Licensee may terminate this License Agreement at any time by destroying the installed version of the
Software (including any backup copies and/or portions thereof) then currently in the Licensee's possession or control. The Licensee
must provide a letter to AG stating that it has destroyed the Software and any copies, should this agreement be terminated. The
license will also terminate automatically upon notice from AG in the event the Licensee fails to comply with any term or condition of
this License Agreement. The Licensee agrees that upon any such termination, it will destroy the Software (including any backup
copies and/or portions thereof). Upon termination, AG may also enforce any and all rights provided by law. The provisions of this
License Agreement that protect the proprietary rights of AG will continue in force after termination.
16. Limited Warranty and Remedies. (a) AG warrants that AG holds title to the Software and has the right to grant this
license (b) AG represents and warrants that to the best of its knowledge, the use of the Software, in compliance with the terms and
conditions of this License Agreement, will not infringe any Canadian, United States or European Union patent or copyright of any
third party. (c) AG warrants the media on which the Software is furnished to be free from defects for a period of ninety days from the
date of receipt. (d) The Licensee assumes full responsibility for the selection of the Software to achieve its intended results, and for
the downloading, use and results obtained from the Software. (e) AG does not warrant that the functions contained in the Software
will meet the Licensee's requirements or that the operation of the Software will be uninterrupted or error-free. (f) AG does not
warrant that the operation of the Software will not be interrupted by reason of any defect or that AG can or will correct all errors or
render all improvements requested. (g) The Software is provided “AS IS” and AG makes no other representation or warranty,
express or implied, and specifically excludes the implied warranties and conditions of merchantability and fitness for a particular
purpose.
17. No Liability for Damages, Including Without Limitation Consequential Damages. Except as otherwise provided
herein, AG and its licensors shall not be liable for any damages whatsoever (including, without limitation, incidental, direct, indirect,
or consequential damages, damages for loss of business profits, business interruption, loss of business information, or other
pecuniary loss) arising out of the use or inability to use this Software, even if AG or its licensors have been advised of the possibility
of such damages In no event will AG’s liability for any claim made or arising during the first 12 months of the term of this Agreement,
whether for breach of contract or warranty, negligence or otherwise, exceed the fee paid by Licensee for the Software that is the
subject of such claim. After the first 12 months, Licensee’s only recourse or remedy for any such claim, damage or loss will be the
termination of this license.
18. Indemnification. The Licensee hereby agrees to indemnify, hold harmless and defend AG from and against any and all
claims or lawsuits, including attorney's fees and costs that arise, resulting from or are connected with the use or distribution of the
Software in violation of this License Agreement.
19. Confidentiality. AG agrees that in the event that it becomes privy to proprietary information or data that the Licensee
considers confidential and Licensee advises AG in writing of the confidential nature of such information, it will not disclose any such
STABView 3.8 User’s Manual A-4 November, 2008
information or data to any third party without the Licensee's written consent, unless: 1. The information is already publicly available;
2. The information was already in AG’s possession prior to the effective date of this Agreement, 3. The information has already been
lawfully and legitimately disclosed to AG by a another source. Confidential information supplied by AG to the Licensee with respect
to the Software’s algorithms, code structure, security and all other information are to be treated as trade secrets of AG, and cannot
be disclosed to other third parties.
20. Security. The Licensee will take all reasonable steps necessary to prevent unauthorized access to and the use of the
Software and its data from any points of access that are under the Licensee’s control.
21. Activation. The license right granted under this Agreement requires activation over the internet when the Software is first
installed and used, or in the event of computer hardware modification, reinstallation, or running multiple instances of the Software. In
the event that the Software is unable to be activated, the Software may fail to execute and the Licensee is required to contact AG
directly.
22. Payment of Applicable Fees. In consideration of AG providing this license and (if applicable) the support and
maintenance services for an agreed upon period of time, the Licensee will pay AG the agreed upon license fee and other fees
according to the terms specified in the Invoice sent to the Licensee. Invoices for this Software must be paid in full before the
Software can be properly activated on an end-user’s computer. AG reserves the right to uninstall the Software by physical removal
or by placing a time limit on the license, should an invoice not be paid on time. Terms of payment are included in the Invoice.
23. General. This License Agreement is binding on the Licensee, its employees, consultants, contractors and agents, and on
any successors, partners and assignees. Neither the Software nor any information derived from it from may be exported except in
accordance with the laws of Canada, the U.S. or other applicable provisions. This License Agreement is governed by the laws of the
Province of Alberta without regard to its conflict of laws rules, and the federal laws applicable therein. The Licensee agrees that AG
will not have any liability for any untrue statement or representation made by its agents or anyone else (whether innocently or
negligently) upon which the Licensee relied upon entering this License Agreement. If any provision of this License Agreement is
deemed invalid or unenforceable by any country or government agency having jurisdiction, that particular provision will be deemed
modified to the extent necessary to make the provision valid and enforceable, and the remaining provisions will remain in full force
and effect. For questions concerning this Agreement, or for product or technical matters, contact AG software support.
24. Compliance with laws. Licensee agrees to comply with all applicable laws, including export controls and economic
sanctions laws imposed by the Governments of Canada and/or the United States. Without limiting the foregoing, licensee agrees
that it shall not export or re-export this Software or any product incorporating this Software without first obtaining all necessary
licenses or approvals. Licensee acknowledges that shipments of this software are subject to the export control laws of the United
States, and accepts that this product may not be exported or re-exported to Cuba, Iran, Syria, Sudan or other countries subject to
U.S. economic sanctions.
APPENDIX B
APPENDIX B
1 + sin φ r
=
1 − sin φ r
OBM oil-based mud
p fluid pressure (within wellbore, perforation, tubular, etc.); pore pressure
pa pore pressure at the borehole or perforation wall
pc capillary threshold pressure
pintact critical bottomhole pressure at which compressive shear yield initiates in intact
rock
pnon-wet pressure in the non-wetting fluid phase of a two-phase system
po native pore pressure, formation pressure, or reservoir pressure
popen critical pressure to re-open a discontinuity
ppc pseudo critical pressure
ppow critical bottomhole pressure at which shear yield initiates on a weak plane
ppr pseudo reduced pressure = p/ppc
pref reference pressure, used for gas pseudo-pressure calculations
psc pressure at standard conditions (i.e., standard atmospheric pressure)
pslip critical pressure for slip on a weak plane or discontinuity
pw bottomhole pressure
pwet pressure in the wetting fluid phase of a two-phase system
∆p difference between pore pressure at the borehole (or perforation) wall and native
pore pressure (or formation pressure or reservoir pressure)
= pa - po
∆posm near-well pressure increase induced by chemical osmosis
∆psurge surge pressure; increase in pressure while running drill pipe into a well.
∆psawb swab pressure; decrease in pressure while running drill pipe out of a well.
Q volumetric flow rate
Qsc volumetric flow rate measured at standard pressure and temperature
R average radius of yielded zone
R universal gas constant (0.0822 litre·atm/mol·K; 8.314 J/mol·K)
rd drainage radius
rperf radius of perforation cavity
rpore pore throat radius
rw borehole radius
s Hoek-Brown criterion strength parameter
SB seabed
SF strength factor
T shear strength parameter for hemispherical cavity sanding model
{ }
= 2 tan 2 (γ ) − 1
T temperature
Tsc temperature at standard conditions (15.5°C)
Tpc pseudo critical temperature
Tpc pseudo reduced temperature T/Tpc
TVD true vertical depth
UCS unconfined compressive strength
u radial displacement
STABView 3.8 User’s Manual B-4 November, 2008
Greek Symbols:
α Biot’s coefficient
=1 – (Kbulk / Kgrain)
αp poroplastic coefficient
αw wedge angle
β turbulence factor (for non-Darcy flow conditions)
Br bulking factor (1.0 to 4.0)
εfc filter cake efficiency
εfp fluid penetration coefficient
ε θT total tangential strain
εθe elastic tangential strain
ε θP plastic tangential strain
ψd dilation angle
γ failure angle
π φp
= +
4 2
γg gas gravity
γH horizontal normalized stress arching ratio
γT interfacial tension
γV vertical normalized stress arching ratio
ϕ contact angle of the fluid interface where it meets the mineral grain surface
µ fluid viscosity
µ average fluid viscosity
υ Poisson’s ratio
ρ fluid density
ρpc pseudo critical density
ρpr pseudo reduced density=ρ/ρpc
Ò1 maximum principal stress
σ2 intermediate principal stress
σ3 minimum principal stress
σ1f value of σ1 at which compressive shear failure occurs
εfc filter-cake or wall-coating efficiency
STABView 3.8 User’s Manual B-5 November, 2008
APPENDIX C
GLOSSARY OF TERMS
STABView 3.8 User’s Manual C-2 November, 2008
APPENDIX C
GLOSSARY OF TERMS
advection: The process by which solute ions or molecules are transported by the motion of
flowing pore water.
apparent capillary cohesion: The small increment of strength in weakly cemented geomaterials
that is attributed to capillary pressures when two fluid phases are present, e.g., water and gas.
axial stress: The normal stress component that acts parallel to the axis of a wellbore or a
cylindrical rock sample.
balanced activity mud: Drilling mud with a water activity equal to the water activity of the
shale being drilled.
Biot coefficient: A constant which determines the magnitude of the total stress component that is
supported by the hydraulic pressure of fluid(s) within the pores of a rock or soil. The magnitude
of this coefficient is a function of the bulk stiffness of a rock or soil relative to the stiffness of its
constituent mineral grains or particles.
Biot effective stress: Total stress minus (pore pressure times the Biot coefficient). This is the
default option for borehole collapse calculations in STABView.
borehole ballooning: Radial expansion of a borehole, typically experienced in shales and
mudrocks in response to high bottomhole pressures that result in yielding or fracturing.
breakout angle: The angle from a borehole centre that subtends the limits of shear failure on
half of the borehole wall.
capillary threshold pressure: The overbalance pressure that must be exceeded in order for a
non-wetting fluid such as oil to penetrate a water-wet rock or soil.
chemical diffusion: The process by which dissolved ions and molecules move from areas of
higher concentration to areas of lower concentration when there is no barrier or membrane.
chemical osmosis: The movement of water through a semi-permeable membrane which
separates solutions of unequal water activities.
claystone: A fine-grained sedimentary rock consisting predominantly of clay minerals, which
displays no lamination or fissility.
cloud point: The critical temperature (which is affected by the salinity of the solution) at which
an emulsion of insoluble glycol droplets in water develops. The resulting emulsion has a cloudy
appearance, as light is scattered by the fine droplets.
compressibility: The change in volume per unit volume per change in pressure for a fluid, rock
or soil (inverse of bulk modulus).
compressive stress: A normal stress tending to shorten a body in the direction in which it acts.
STABView 3.8 User’s Manual C-3 November, 2008
contact angle: The angle of intersection of the interface between two fluids at a solid surface.
Used in the estimation of capillary threshold pressure.
convection: The process by which fluid moves from areas of higher hydraulic pressure to areas
of lower hydraulic pressure (also see hydraulic diffusion).
dielectric constant: The factor by which the force between two charged particles is decreased
due to the presence of a solid, liquid or gas in the space between the particles.
diffusion osmosis: The movement of solute and associated water through a non-ideal semi-
permeable membrane which separates solutions of unequal solute concentrations.
discontinuity: A collective term for all fractures or planar features in a rock mass such as weak
bedding planes, faults, natural fractures and joints having little or relatively low tensile strength.
drilling fluid: The fluid phase of a drilling mud, including any emulsified liquid phases and/or
dissolved solutes (e.g., salts).
drilling mud: A suspension, usually in water but sometimes in oil or synthetic non-polar liquids,
typically circulated down the drill pipe, through the drill bit, and back to surface through the
annulus between the drill pipe and the borehole wall. Common additives include bentonitic
clays, dissolved ions such as salts, weighting agents such as barite, and polymers.
drilling mud filtrate: The drilling mud constituents that are able to permeate formations
surrounding a wellbore.
Drucker-Prager failure criterion: A three dimensional failure criterion for frictional materials
that is popular in materials engineering. Its rupture surface is a cone in principal stress space.
ductile: A property of a material that can sustain permanent deformation without losing the
ability to sustain load.
effective salinity of a shale pore fluid: The salt concentration in a pure aqueous solution that
would result in a water activity equal to the true value for the pore fluid. The latter value can be
affected not only by dissolved ions, but also by clays, non-ionic solutes, and any other
components present in the shale.
effective stress: The average normal force per unit area transmitted from grain to grain in a rock
or soil mass. Effective stress is equal to the total stress less the pore pressure in most soils and
many rocks.
elasticity: Property of a material that returns to its original form or condition after the applied
force is removed.
electro-osmosis: The movement of liquid through a semi-permeable membrane under the
influence of an applied electric field.
equivalent circulating density (ECD): The density of a fluid which includes a component that
accounts for friction losses due to circulation, and swab or surge pressures.
STABView 3.8 User’s Manual C-4 November, 2008
equivalent mud weight (EMW): In this software EMW is defined as the static mud density plus
or minus a swab or surge pressure. It does not apply to dynamic situations where the fluid is
being circulated.
fabric: The spatial and geometrical configuration of the grains, particles, crystals and cement of
which a sedimentary rock is composed.
factor of safety: See safety factor.
failure: The condition of a material or structure that can no longer adequately support the forces
applied to it or otherwise perform its engineering function.
failure criterion: Specification of the mechanical condition under which materials fail by
fracturing or by deforming beyond some specified limit. This specification may be in terms of
stresses, strains, stress rates, strain rates, or some combination of these quantities in the
materials.
filter cake efficiency: This parameter characterizes the extent to which pore pressure at the
borehole wall penetrates into the formation through a filtercake in response to an overbalance
pressure.
fissility: The tendency to split easily along closely spaced, parallel planes. Shales, by definition,
possess fissility but claystones or mudstones do not.
fluid penetration coefficient: Parameter used to describe the ability of a fluid to penetrate into
the surrounding formation.
formation pressure: The hydraulic pressure of fluid(s) within the pores of a rock formation.
fracture breakdown pressure: The hydraulic pressure at which a tensile fracture initiates at the
wellbore wall. Sometimes called the fracture initiation pressure.
fracture closure pressure (FCP): The fluid pressure in a hydraulic fracture just at the point
where the fracture closes. This pressure is equal to and counteracts the minimum principal stress
in the rock perpendicular to the fracture plane.
fracture propagation pressure (FPP): The fluid pressure required to extend an existing
hydraulic fracture. It is greater than the fracture closure pressure and depends on the size and
shape of the fracture, the roughness of the fracture, the fluid rheology and other factors.
friction angle: The rate of increase in the shear strength of a geomaterial with increasing normal
stress.
fugacity: A thermodynamic function defined by the equation dG = RT d ln f , where G is the
Gibb’s free energy, R is the universal gas constant, T is absolute temperature and f is the
fugacity. Fugacity is expressed in units of pressure.
Gibb’s free energy: A thermodynamic state function related to the enthalpy, temperature and
entropy of a system.
STABView 3.8 User’s Manual C-5 November, 2008
glycol: The common name for 1,2-ethanediol, also known as ethylene glycol, which is an
alcohol containing two OH groups. Glycol drilling muds are often based on more complex
polyglycol structures.
hoop stress: Same as tangential stress.
Hoek-Brown failure criterion: A non-linear empirical rock strength criterion to describe the
brittle failure of intact rock and/or fractured rock masses.
hydration: The movement of water into the interlayer sites within clay particles.
hydrational stress: The mechanical stress that develops as a consequence of hydration.
hydraulic diffusion: The process by which fluid moves from areas of higher hydraulic pressure
to areas of lower hydraulic pressure (also see convection).
hydraulic flow: See hydraulic diffusion.
hydraulic pressure: The force per unit area exerted by a fluid.
hydrostatic pressure: The pressure exerted by a fluid at rest at a given depth.
hydrostatic stress: A state of stress in a material in which the normal stresses acting on any
plane are equal and where shearing stresses do not exist.
inhibition: The arresting or slowing of the hydration, swelling and disintegration of reactive
clays in claystones, mudstones and shales.
in-situ stress: The state of stress in a rock formation prior to disturbance by drilling or
production activities.
interfacial tension: The surface tension at the interface between two fluids.
interlayer water: The water molecules located between clay platelets. Interlayer water that is
able to move due to hydraulic gradients is referred to as “free” water, and water that is adsorbed
onto clay platelet surfaces and is restricted from flow in normal conditions is referred to as
“bound” water.
instantaneous shut-in pressure (ISIP): The fluid pressure in a hydraulic fracture immediately
after shut-in; i.e., after pumping has ceased.
ion exchange: A process by which an ion in the crystal structure of a mineral is replaced by
another ion that was present in an aqueous solution.
linear elasticity: A specific type of elasticity in which the strain in a material is linearly
proportional to the applied stress.
lost circulation: The loss of a drilling mud, cement or completion fluid into natural or induced
fractures, vugs or high-porosity zones, resulting in the inability to circulate fluids to surface.
maximum principal stress: The largest of the three principal stress components.
STABView 3.8 User’s Manual C-6 November, 2008
maximum horizontal in-situ stress: The larger principal stress component acting in a horizontal
plane within the earth’s crust.
membrane efficiency: A parameter related to the relative mobility of solvent and solute
particles through a semi-permeable membrane. It is defined numerically as the ratio of the
observed osmotic pressure which develops across a semi-permeable membrane separating
solutions of unequal water activities to the osmotic pressure predicted for an ideal semi-
permeable membrane.
methyl glucoside: A polyol monomer which has a compact, two-tiered cyclic structure that is a
component in some low activity drilling fluid systems.
minimum principal stress: The smaller of the three principal stress components.
minimum horizontal in-situ stress: The smaller principal stress component acting in a
horizontal plane within the earth’s crust.
Mohr-Coulomb failure criterion: A linear rock failure criteria characterized by a friction angle
and a cohesion intercept. The rupture surface forms an irregular pyramid in principal stress
space.
moisture content: The weight of water contained in the pore space of a rock divided by the
weight of solid material.
mud pressure penetration: A process by which high hydraulic pressures in a wellbore diffuse
into the rock formations surrounding the wellbore, resulting in a near-well increase in pore
pressure.
mudstone: A fine-grained sedimentary rock which displays no lamination or fissility and
contains approximately equal proportions of clay and silt.
non-linear failure criterion: A failure criterion where the maximum principal stress upon
yielding is non-linearly related to the minimum principal stress.
normal strain: The change in length per unit length normal to a given plane.
normal stress: The stress component normal to a given plane.
oil-based emulsion: A fluid consisting of a continuous oil phase with emulsified water or brine
droplets.
oil-based mud (OBM): A drilling mud in which the external phase is an oil-related product,
such as diesel oil or mineral oil.
osmosis: The movement of water from one aqueous system to another through a semi-permeable
membrane driven by activity differences between the two systems.
osmotic membrane efficiency: See membrane efficiency.
osmotic pressure: The hydraulic pressure difference across a semi-permeable membrane
resulting from a difference in water activities of the solutions on either side of the membrane.
STABView 3.8 User’s Manual C-7 November, 2008
osmotic pressure coefficient: The ratio of the observed osmotic pressure which develops across
a semi-permeable membrane separating solutions of unequal water activities to the osmotic
pressure predicted for an ideal semi-permeable membrane (also see membrane efficiency and
reflection coefficient).
overbalanced drilling: Drilling when the pressure in the wellbore exceeds the pressure of fluids
in the formation. This excess pressure prevents reservoir fluids (oil, gas or water) from entering
the wellbore
overbalance pressure: The difference between the hydraulic pressure in the wellbore and the
formation pore pressure, for cases where the former is greater than the latter.
overburden stress or pressure: (See vertical in-situ stress).
partial molar volume: The volume of a given component in 1 mole (6.02×1023 molecules) of a
solution.
peak strength: The maximum shear stress that a rock can sustain under a given set of
conditions.
physico-chemical interaction: A process or reaction in which matter and/or chemical energy are
exchanged between substances.
plastic: The property of a material whose deformation is not recoverable when the applied load
is removed.
Poisson’s ratio: An elastic property that is a measure of the deformability of a material
perpendicular to the applied stress. It is expressed as the ratio of the lateral to the longitudinal
strain in a test sample.
pore: The volume between mineral grains and clay particles in a rock or soil mass.
pore pressure: The hydraulic pressure at a given point of fluid(s) within the pores of a rock or
soil mass.
pore throat: A narrow channel, where mineral grains or particles are in close contact,
connecting larger interstitial pore volumes.
principal stress: A stress acting normal to one of three mutually perpendicular planes
intersecting at a point in a body, on which the shearing stress is zero.
radial stress: The normal stress component that acts in a direction pointing radially outward
from the centre of a wellbore or a cylindrical rock sample.
reflection coefficient: (See membrane efficiency).
relative humidity or relative water vapor: The vapor pressure of water in a solution, drilling
mud or shale divided by the vapor pressure of pure water at the same temperature. This ratio is
usually expressed as a percentage.
STABView 3.8 User’s Manual C-8 November, 2008
residual strength: The minimum shear strength that a rock drops to after the peak strength has
been exceeded.
rock failure: The process where the maximum strength of a rock, or the stress or strain
requirement of a specific design, is exceeded.
rock yielding: The process where rock undergoes deformations which are not recovered when
the applied stress is removed.
safety factor: The ratio of the strength of a material to the stresses acting at a point in that
material. Many specific forms of safety factors exist, depending on the type of failure criterion
used, both for shear and tensile yielding.
semi-permeable membrane: An interface which is permeable to water (or some other solvent)
but impermeable to solute ions or molecules.
shale: A fine-grained detrital sedimentary rock which has a finely laminated structure and
exhibits fissility.
shear strain: A measure of the amount by which parallel lines have been sheared past one
another by deformation. The ordinary or engineering definition of this parameter is the tangent of
the change in angle between initially perpendicular lines. The infinitesimal-strain theory
definition is half the tangent of the change in angle between initially perpendicular lines.
shear stress: The component of stress acting tangential to a given plane.
solute: The minor component of a solution, or the component that is solid in its pure state.
solution: A homogenous mixture of two or more substances.
solvent: The major component of a solution, or the component that is liquid in its pure state.
stiffness: The deformation of a rock specimen divided by the normal or shear stress that caused
the deformation.
stress: The force per unit area acting within a rock mass.
stress-depletion: The pore pressure changes associated with hydrocarbon production or injection
which alters the stress state in the reservoirs.
support pressure: The hydraulic pressure drop across the wellbore wall, which acts as a radially
oriented normal stress that mechanically supports the rock around the wellbore.
surface tension: The force per unit length that acts across any line in a surface, tending to pull
the surface open.
surge pressure: A wellbore pressure increase caused by a downward pipe movement in a well.
swab pressure: A wellbore pressure decrease caused by an upward pipe movement in a well.
STABView 3.8 User’s Manual C-9 November, 2008
swelling: A volume increase resulting from the hydration of clay or other water-absorbing
minerals.
swelling pressure: The stress required to prevent swelling when a water-sensitive rock or soil is
contacted with a fluid that has a different chemical composition than its pore fluid.
tangential stress: The normal stress component that acts in a direction tangential to the
circumference of a wellbore at a given point.
tectonically stressed: A condition where the large-scale deformations occurring in the upper
crust of the earth have resulted in elevated in-situ stresses in a region, especially the maximum
horizontal in-situ stress.
tensile stress: A normal stress tending to lengthen a body in the direction in which it acts.
tensor (second order): A mathematical entity expressed as a 3×3 matrix which describes the
relationships between two vector entities. In the case of stresses, one set of vectors describes the
three orthogonal planes for a chosen reference orientation, and the other set of vectors describes
the magnitudes of the stresses which act in normal and shear directions on these planes.
Terzaghi effective stress: Total stress minus pore pressure, where the Biot coefficient is equal
to 1.0.
total stress: The total force per unit area acting within a mass of rock, which includes both the
intergranular or effective stress and the pore pressure.
triaxial strength test: A test in which a cylindrical specimen of rock encased in an impervious
membrane is subjected to a confining pressure and then loaded axially to failure.
unconfined compressive strength (UCS): The peak strength of an axially loaded rock or soil
specimen with no lateral confinement.
underbalanced drilling (UBD): Drilling when the pressure in the formation exceeds the
pressure of fluids in the wellbore.
underbalance pressure: The magnitude of the difference between the fluid pressure in a
wellbore and the formation pore pressure, for cases where the former is less than the latter.
vapor pressure: The pressure exerted by a vapor in equilibrium with its liquid (or solid) phase.
vertical in-situ stress: The principal stress component oriented vertically within the earth’s
crust, generated by the weight of the overburden.
water activity: The ratio of the vapor pressure of water in a solution, drilling mud or shale pore
to the vapor pressure of pure water at the same temperature. [Note: The strict definition of
activity is in terms of fugacity, but for practical purposes vapor pressures are sufficiently
accurate.
water-based mud (WBM): A drilling fluid in which water or brine is the major liquid phase and
the wetting (external) phase.
STABView 3.8 User’s Manual C-10 November, 2008
water-wet shale: A shale in which the wettability of mineral grain and particle surfaces is
greater for water than for other liquids such as oil.
wellbore failure: A condition where a wellbore is no longer able to perform its intended
operational purpose.
wellbore instability: A condition where rock around a wellbore is undergoing yield and/or
failure.
wettability: The preference of a solid to contact one liquid or gas, known as the wetting phase,
rather than another.
yield: The condition where a solid material begins to undergo irreversible deformation.
yield criterion: Specification of the mechanical condition under which solid materials begin to
undergo irreversible deformation. This specification may be in terms of stresses, strains, stress
rates, strain rates, or some combination of these quantities in the materials.
Young’s modulus: An elastic property that is the ratio of longitudinal stress to longitudinal
strain.
________________________________________________________________________
Sources Consulted for this Glossary
American Association for Testing and Materials, Standard Definitions of Terms and Symbols
Relating to Soil and Rock Mechanics, ASTM D 653 – 80, 1981.
American Geological Institute, Glossary of Geology, R.L. Bates and J.A. Jackson (eds.),
American Geological Institute, Falls Church, Virginia, 1980, 751 p.
American Geological Institute, Dictionary of Geological Terms, R.L. Bates and J.A. Jackson
(eds.), Anchor Press/Doubleday, New York, 1984, 571 p.
Fetter, C.W., Applied Hydrogeology, 3rd edition, Maxwell Macmillan, Toronto, 1994, 691 p.
Giancoli, D.C., Physics: Principles with Applications, 2nd edition, Prentice Hall, Toronto, 1985,
811 p.
Gillespie, R.J., Humphreys, D.A., Baird, N.C. and Robinson, E.A., Chemistry, Allyn and Bacon,
Toronto, 1986, 891 p.
STABView 3.8 User’s Manual D-1 November, 2008
APPENDIX D
APPENDIX D
ANALYSIS OF BOREHOLE ENLARGEMENT IN PERMAFROST
D.1 Model Description
When drilling through permafrost, a problem arises when ice contained in the pores of the
formation melts as a consequence of circulating drilling fluids with a higher temperature. To
estimate the maximum borehole enlargement due to permafrost melting, the model developed by
Kutasov and Caruthers (1988) has been implemented in Advanced Geotechnology’s STABView
software. Additional background on this class of problems can be found in Kutasov (1999) and
Kudryashov and Yakovlev (1983). The following equations were programmed as shown on the
flowchart, Figure D.1, to calculate K, the washout coefficient.
Laminar Flow
2
1 ⎛1 ⎞ Tm k m t
K= + ⎜ − 1 ⎟ + 57.6 × 10 6 (D.1)
D ⎝D ⎠ L S i φ f ρ i d 2b
Turbulent Flow
B Tm k m t
57 .6 × 10 6 = F ( KD ) − F ( D ) (D.2)
L S i φ f ρ i d p2
where:
5 10
F (x) = ( x + 1) 2 . 8 − ( x + 1) 1 .8 (D.3)
14 9
0 .8 0 .4
⎛ 4 × 10 9 q ⎞ ⎛ υK ρ C ⎞
B = 0 . 046 ⎜ ⎟ ⎜⎜ ⎟⎟ (D.4)
⎜ πυ d ⎟
⎝ K p ⎠ ⎝ 1000 k m ⎠
Input Parameters
Y
NRe>2100
Y
N Calculate washout.
Laminar Flow (Eqn. D.1)
ti < t
End
Figure D.1: Flowchart showing the permafrost hole enlargement routine in STABView.
The criterion used to define the flow regime is based on Reynolds Number (NRe). Using
Kutasov’s approach the flow is considered turbulent when NRe ≥ 2100 and laminar for NRe
<2100. Reynolds Number is calculated as follows:
v(d h − d p )
N = (D.5)
1000 υ
Re
K
vd = υ K / ρ (D.6)
1. Ice is considered the only cementing material within the permafrost interval. Once the ice
melts, the soil is assumed to washout and become part of the mud. This is a conservative
assumption for any materials with appreciable cohesion or cementation. Furthermore, in
permeable materials drilled overbalanced, a filter cake will build which also provides
support to the wellbore wall, preventing collapse and/or washout. To assess this mode of
failure an elastoplastic borehole stability analysis could be run in STABView to predict
the size of the potential yielded zone and the related deformations.
2. The initial formation temperature is assumed to be 0°C. This is clearly a limitation of the
Kutasov approach for colder permafrost.
3. While drilling, thawed material is entrained in and removed from the borehole by the
drilling fluid. Hole cleaning calculations are not accounted for since the model assumes
that all of the dislodged material can be lifted from the hole despite the low annular
velocities that result in the ever increasing hole size.
4. These equations were developed for a vertical well.
5. Additional turbulence developed by pipe rotation has not been accounted for.
6. Scour due to hydraulic jet impact forces at the bit has not been included.
7. Erosion due to shearing at the fluid - wall contact is assumed to be 100% efficient, i.e., all
of the material is removed or melted and becomes part of the drilling fluid.
8. In this implementation of STABView we have assumed that the permafrost interval is
100% saturated with ice, i.e., no appreciable unfrozen water or gas is present.
The following improvements would make the model more suitable for practical applications:
1. Include additional drilling fluid rheological properties, i.e., Power Law and Hershel
Buckley models.
2. Include variable permafrost temperatures below 0 ºC.
3. Include turbulence effects due to pipe rotation by considering the Taylor number.
STABView 3.8 User’s Manual D-5 November, 2008
4. Develop a new option for back-analyses that will use a time history of various input
parameters such as the drilling fluid temperature, rheological properties, fluid thermal
properties, and pumping rate.
Dallimore, S.R., Uchida, T., and Collett, T.S., Scientific Results from JAPEX / JNOC / GSC
Mallik 2L-38 Gas Hydrate Research Well, Mackenzie Delta, Northwest Territories,
Canada. Geological Survey of Canada Bulletin 544, 1999.
Dallimore, S.R., Laframboise, R.R., Fotiou, M.J., and Medioli, B.E., JAPEX / JNOC / GSC
Mallik 2L-38 Gas Hydrate Research Well, Mackenzie Delta, Northwest Territories,
Canada: Interactive Data Viewer, GSC Open File D3726, CD-ROM, November,
1999.
Kudryashov, B.B., and Yakovlev, A.M., Drilling in Permafrost (translated from Russian),
Nedra Publishers, Moscow, 318p, 1983.
Kutasov, I.M., Applied Geothermics for Petroleum Engineers, Elsevier, 1999.
Kutasov, I.M. and Caruthers R.M., Hole Enlargement Control During Arctic Drilling, SPE
17442, SPE California Regional Meeting, Long Beach, California, 1988.
STABView 3.8 User’s Manual D-6 November, 2008
Table D.1: Example input and output from STABView 3.8 for
a hole enlargement case in permafrost
CASE DATA
Case Title Hole Enlargement in the Permafrost Section at 100 mKB
Well Name Research Well Mallik 2L-38, Mackenzie Delta, NWT
Company Advanced Geotechnology Inc.
Client Geological Survey of Canada
User name GS/PM
Unit Type Metric (Modified SI)
WELL DATA
Offshore/Onshore Onshore
KB Height (relative to ground level) 7.1 m
Well Survey Vertical Well
FORMATION DATA
Stratigraphic Unit Iperk Sequence
Lithology unconsolidated sediment
SINGLE DEPTH ANALYSIS DATA
Figure D3: Predicted Reynold’s number with time in the permafrost interval
at 100 mKB for the example problem.
STABView 3.8 User’s Manual D-8 November, 2008
STABView 3.8 User’s Manual E-1 November, 2008
APPENDIX E
Table E.1: STABView well planning sample cases typically included with STABView as demonstration files
No. STABView Type of No. Zones Analysis Analysis Well Options/Features Selected Comments
Sample Case Model Type Mode Trajectory
Name
1 3D Drilling 3D - elastic single Hole Pressure inclined Effect of a plane of weakness on Evaluation of the effects of a fissile shale on wellbore stability
Collapse Case collapse wellbore collapse. Collapse with a calibrated 3D elastic model.
with Weak while calibration.
Bedding drilling
2 3D Hole 3D - elastic multiple Hole Risk inclined Failure criteria sensitivity. Depth Pseudo-yielded zone calculation for an inclined well based on
Enlargement collapse profile. the extent of the failed zone determined from the default 3D
Model while linear elastic model. It is interesting to compare the results
drilling obtained when Modified Lade and Hoek-Brown failure criteria
are used.
3 Horizontal Well 2D - elasto- single Sand Risk horizontal Steady - state pore pressure gradient Horizontal sand production case to decide if a slotted liner,
2D Sand plastic production effects. Depletion induced stresses. screen or gravel pack are required. Yielded zone size and
Production Risk Barefoot openhole completion. deformations predicted for various sensitivities, e.g. dilation
Calculation of the rubble fill angle, peak cohesion, residual cohesion.
percentage.
4 Horizontal Well 2D - elasto- multiple Open hole Risk horizontal Liner, screen or pre-packed screen Horizontal sand production case to assess if hole closure results
with Liner plastic stability/ completion. Calculation of in contact between the formation and a stiff slotted liner or
sand formation stresses on a stiff screen screen used to complete this well.
production or liner due to hole closure.
Calculation of plastic strain at the
wellbore wall.
5 Extended Reach 3D - elastic multiple Hole Pressure build plus Collapse calibration. Shale-mud Determine safe range of EMW and ECD to prevent borehole
Well with collapse and horizontal interaction and osmotic pressure collapse and fracture breakdown in the build section of the well
Horizontal Leg fracture (survey) effects on wellbore collapse. Mud and within a depleted reservoir section of a horizontal well.
breakdown weight window.
6 Inclined Well 3D - elastic multiple Fracture Pressure survey Steady - state thermal effects on lost Effect of cool fluid injection on fracture breakdown pressure
with Cold Fluid breakdown circulation and fracturing risks. for an open hole completion in a directional well.
Injection
7 Injection Well 3D - elastic single Lost Pressure vertical 3D slip risk analysis on a weak Analysis to determine if elevated pore pressures during water
Slip Pressure Circulation discontinuity, e.g., natural fractures, injection will cause shear failure on weak natural fractures or
Analysis faults or bedding. bedding planes.
8 Permafrost Hole 2D - single Hole Risk vertical 2D permafrost hole enlargement Example of a vertical well drilled through permafrost
Enlargement permafrost collapse model.
erosion while
model drilling
9 Permafrost Hole 2D - multiple Hole Risk vertical 2D permafrost hole enlargement Example of a vertical well drilled through rich, weak erodable
Enlargement permafrost collapse model. permafrost intervals in the Mackenzie Delta, N.W.T., Canada.
Multi - Zone erosion while
model drilling
10 Underbalanced 3D – elastic multiple Hole Pressure survey Steady - state pore pressure gradient Multi-zone wellbore stability analysis for an offshore, extended
Drilling in Fissile and 2D collapse effects. Depletion induced stresses. reach underbalanced well penetrating a fissile shale at angle.
Shale elasto-plastic while Effects of a plane of weakness on Case takes several minutes to run 2D elasto-plastic and 3D
drilling wellbore collapse. Collapse elastic models compared for high angle penetrations.
calibration. Mud weight window.
STABView 3.8 User’s Manual E-3 November, 2008
Table E.1: STABView well planning sample cases typically included with STABView as demonstration files (Cont.)
No. STABView Type of No. Zones Analysis Analysis Well Options/Features Selected Comments
Sample Case Model Type Mode Trajectory
Name
11 Underbalanced 2D - elasto- multiple Hole Risk survey Steady - state pore pressure Multi-zone wellbore stability analysis for an offshore, extended
Drilling plastic collapse gradient effects. Stress depletion reach underbalanced well. Bottomhole pressure profile used as
while model for changing reservoir input.
drilling pressures. Yielded zone
predictions.
12 Vertical Gas Well 2D - elasto- multiple Open hole Risk vertical Barefoot openhole completion. Vertical gas well with barefoot open hole completion. The user
Open hole plastic stability Effects of gas properties, has input a gas flow rate, from which STABView calculates the
Completion under temperature and non-Darcy flow bottomhole pressure for its use in a borehole collapse model.
drawdown rates on openhole stability.
13 Vertical Oil Well 2D - elasto- multiple Open hole Risk vertical Perforated liner completion with Vertical well with perforated completion. User has specified a
plastic stability cylindrical perforation cavities. liquid production rate, from which STABView calculates the
under Calculation of the perforation bottomhole pressure based on the steady-state radial flow
drawdown stability for the most critical equation for an incompressible fluid.
perforation orientation.
14 Vertical Well with 2D - elasto- multiple Hole Risk vertical Steady - state pore pressure 2D elasto-plastic analysis of the yielded zone size for a vertical
Collapse Risk plastic collapse gradient effects (two-zone mobility well drilled overbalanced through multiple zones.
while model).
drilling
15 Inclined Well 3D - elastic single Sand Risk/ inclined Effect of well trajectory and Calibrated 3D analysis model for cylindrical and hemi-spherical
with Perforations production Pressure perforation orientation on fracture perforation cavities. Assess the value of oriented perforating.
breakdown pressures and sand
production risk
16 Horizontal Well 2D - elasto- multiple Sand Pressure survey Effect of reservoir depletion on the 2D elastoplastic analysis model used for examining the effects
with Effects of plastic Production sand production critical bottom of reservoir depletion for a range of critical bottom hole
Depletion 2D hole pressure or drawdown pressures.
17 Complex Multi- 2D multiple Hole Pressure survey Demonstrates using a combination A complex drilling case showing multiple zones and models
Zone Offshore elastoplastic & collapse and of 2D/3D analysis models in a used to calculate wellbore stability collapse and fracture
Well 3D elastic fracture drilling problem in order to pressures for an offshore well.
breakdown determine a safe operating mud
weight window.
STABView 3.8 User’s Manual E-4 November, 2008
Case Title: 3D Drilling Collapse Case - Weak Bedding 1 of 4 Andes Exploration S.A.
CASE PROFILE
Analysis Type Hole Collapse While Drilling
Analysis Mode Calculate required bottomhole pressures or mud densities to achieve tolerable borehole
instability risks.
Borehole Collapse Model Type Default (3D Linear Elastic Model with Mohr-Coulomb Failure Criterion)
Date and Time 26-Nov-2008, 14:18
Filename: Case 1 - 3D Drilling Collapse Case with Weak Bedding Effects.bhs 26-Nov-2008 14:18
Case Title: 3D Drilling Collapse Case - Weak Bedding 2 of 4 Andes Exploration S.A.
Filename: Case 1 - 3D Drilling Collapse Case with Weak Bedding Effects.bhs 26-Nov-2008 14:18
Case Title: 3D Drilling Collapse Case - Weak Bedding 3 of 4 Andes Exploration S.A.
Äçä = 1.00
cñ = 290 psi
Åñ = 35 °
E = 0.70 E+06 psi
10.5 Í = 0.30
À = 1.00
Plane of Weakness
Dip = 35 °
Dip Direction = 300 °
cç = 29 psi
Åç = 20 °
10.0
0 30 60 90 120 150 180 210 240 270 300 330 360
Formation: Shale
12.0 Lithology: shale
Äçä = 1.00
cñ = 290 psi
10.5 Åñ = 35 °
E = 0.70 E+06 psi
Í = 0.30
À = 1.00
10.0 Plane of Weakness
Dip = 35 °
Dip Direction = 300 °
cç = 29 psi
Åç = 20 °
9.5
0 10 20 30 40 50 60 70 80 90
Filename: Case 1 - 3D Drilling Collapse Case with Weak Bedding Effects.bhs 26-Nov-2008 14:18
Case Title: 3D Drilling Collapse Case - Weak Bedding 4 of 4 Andes Exploration S.A.
N
Analysis Mode: Pressure 12.07
Analysis Type: Drilling Collapse
Äçä = 1.00
cñ = 290 psi
Åñ = 35 ° 10.14
E = 0.70 E+06 psi
Í = 0.30 ppg
À = 1.00
Bedding Friction 15 20 25 °
Angle
Bedding Dip
26 35 44 °
Filename: Case 1 - 3D Drilling Collapse Case with Weak Bedding Effects.bhs 26-Nov-2008 14:18
Case Title: 3D Elastic Yielded Zone Analysis for an Inclined Well 1 of 7 Hornet Exploration Ltd.
CASE PROFILE
Analysis Type Hole Collapse While Drilling
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Date and Time 26-Nov-2008, 14:23
BOREHOLE PROPERTIES
Bottom Borehole
No. Depth Diameter
ft KB MD in
1 4005.12 8.66
TUBULAR PROPERTIES
* No tubular data were entered for this case.
STRATIGRAPHIC UNITS
No. Stratigraphic Unit Lithology Pattern Top Depth Btm Depth Top Depth Btm Depth
ft KB MD ft KB MD ft KB TVD ft KB TVD
1 Interval 1 dolomite 3684.66 3905.01 3018.30 3198.80
CASE PROFILE
Analysis Type Hole Collapse While Drilling
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Date and Time 26-Nov-2008, 14:23
Unit No. 1 - Stratigraphic Unit: Interval 1, Lithology: dolomite, Model Type: Default (3D Linear Elastic Model with Mohr-Coulomb Failure
Criterion)
Lithology: dolomite
14
Model: 3D Linear Elastic
Failure Criteria: Mohr Coulomb
Äçä = 1.00
10 cñ = 290 psi
Åñ = 35 °
E = 1.16 E+06 psi
Í = 0.30
À = 1.00
8
0 30 60 90 120 150 180 210 240 270 300 330 360
Azimuth of Min
90 135 180 °
Horz Stress
Base Case R/rw = 1.10
CASE PROFILE
Analysis Type Sand Production / Openhole Stability Under Drawdown
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Borehole Collapse Model Type 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
Date and Time 26-Nov-2008, 14:24
cñ = 800 psi
Åñ = 40 °
1.2 có = 220 psi
Åó = 33 °
E = 0.20 E+06 psi
Í = 0.30
À = 1.00
1.1
500 600 700 800 900
Model: 2D Elastoplastic
0.22 Failure Criteria: Mohr Coulomb
Òüîêï Azi = 45 °
cñ = 800 psi
0.16 Åñ = 40 °
có = 220 psi
Åó = 33 °
E = 0.20 E+06 psi
0.14 Í = 0.30
À = 1.00
0.12
500 600 700 800 900
700
0 5 10 15
Model: 2D Elastoplastic
Failure Criteria: Mohr Coulomb
4000 Ò÷ grad = 0.97 psi/ft
Òüîâù grad = 0.93 psi/ft
Òüîêï grad = 0.80 psi/ft
psi
Òüîêï Azi = 45 °
3000 Pó grad = 0.40 psi/ft
Pó = 1600 psi
Pó grad(new) = 0.22 psi/ft
Pó(new) = 880 psi
Pø = 700 psi
2000 ¿P = -180 psi
cñ = 800 psi
Åñ = 40 °
có = 220 psi
1000 Åó = 33 °
E = 0.20 E+06 psi
Í = 0.30
À = 1.00
0
0 5 10 15
Lithology: sandstone
Model: 2D Elastoplastic
50 Failure Criteria: Mohr Coulomb
cñ = 800 psi
35 Åñ = 40 °
có = 220 psi
Åó = 33 °
E = 0.20 E+06 psi
30 Í = 0.30
À = 1.00
25
500 600 700 800 900
Peak Cohesion
600 800 1000 psi
Residual Cohesion
165 220 275 psi
Pressure Drop
5 20 100 psi
Across Liner
CASE PROFILE
Analysis Type Sand Production / Openhole Stability Under Drawdown
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Date and Time 26-Nov-2008, 14:26
BOREHOLE PROPERTIES
Bottom Borehole
No. Depth Diameter
ft KB MD in
1 4015.01 12.00
2 5600.00 8.50
TUBULAR PROPERTIES
No. Tubular Type Top Depth Btm Depth ID OD
ft KB MD ft KB MD in in
1 casing 0.00 4015.01 9.00 9.60
2 slotted liner 4015.01 5600.00 6.60 7.00
STRATIGRAPHIC UNITS
No. Stratigraphic Unit Lithology Pattern Top Depth Btm Depth Top Depth Btm Depth
ft KB MD ft KB MD ft KB TVD ft KB TVD
1 Eocene Sand sandstone 4000.00 5600.00 3930.06 3974.14
INITIAL FORMATION STRESS AND PRESSURE GRADIENT DATA
No. Stratigraphic Unit Top Top SHmin Sv Sv SHmax SHmax SHmin SHmin Pr Pr
Depth Depth Azi. Grad Grad Grad Grad Grad Grad Grad Grad
ft KB MD ft KB TVD (Top) (Btm) (Top) (Btm) (Top) (Btm) (Top) (Btm)
psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft
1 Eocene Sand 4000.00 3930.06 45 0.95 0.95 0.95 0.95 0.95 0.95 0.40 0.40
WELLBORE PRESSURE OR PRODUCTION RATE DATA
Top Top Wellbore Wellbore
No. Stratigraphic Unit Lithology Depth Depth Pressure (Top) Pressure (Btm)
ft KB MD ft KB TVD psi psi
1 Eocene Sand sandstone 4000.00 3930.06 950 950
BASIC ROCK MECHANICAL PROPERTIES
No. Stratigraphic Unit Lithology Top Top Peak Peak Apparent Poisson's Young's Biot's
Depth Depth Cohesion Friction Capillary Ratio Modulus Coefficient
ft KB MD ft KB psi Angle Cohesion E+06 psi
TVD psi
1 Eocene Sand sandstone 4000.00 3930.06 650 40 0.00 0.30 0.15 1.00
FLUID PENETRATION PROPERTIES
No. Stratigraphic Unit Lithology Top Top Capillary Filter Cake
Depth Depth Threshold or Wall
ft KB MD ft KB TVD Pressure Coating
psi Efficiency
CASE PROFILE
Analysis Type Sand Production / Openhole Stability Under Drawdown
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Date and Time 26-Nov-2008, 14:26
Unit No. 1 - Stratigraphic Unit: Eocene Sand, Lithology: sandstone, Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
Model: 2D Elastoplastic
Failure Criteria: Mohr Coulomb
cñ = 650 psi
Åñ = 40 °
có = 145 psi
Åó = 35 °
E = 0.15 E+06 psi
0.4 Í = 0.30
À = 1.00
0.3
500 600 700 800 900 1000
500
1000
TVD (ft)
1500
2000
2500
3000
3500
4000
4500
-1500 -1000 -500 0 500 1000 1500 2000 2500 3000 3500
Well Inclination = 86 °
Well Azimuth = 315 °
1.8 Formation: Eocene Sand
Lithology: sandstone
Model: 2D Elastoplastic
Failure Criteria: Mohr Coulomb
1.6
Ò÷ grad = 0.95 psi/ft
Òüîâù grad = 0.95 psi/ft
Òüîêï grad = 0.95 psi/ft
Òüîêï Azi = 45 °
1.4
Pó grad = 0.40 psi/ft
Pó = 1587 psi
cñ = 650 psi
1.2 Åñ = 40 °
có = 145 psi
Åó = 35 °
E = 0.15 E+06 psi
Í = 0.30
1.0 À = 1.00
0.8
700 800 900 1000 1100 1200
Model: 2D Elastoplastic
Failure Criteria: Mohr Coulomb
1.8
Ò÷ grad = 0.95 psi/ft
Òüîâù grad = 0.95 psi/ft
Òüîêï grad = 0.95 psi/ft
Òüîêï Azi = 45 °
1.6
Pó grad = 0.40 psi/ft
Pó = 1587 psi
cñ = 650 psi
1.4 Åñ = 40 °
có = 145 psi
Åó = 35 °
E = 0.15 E+06 psi
Í = 0.30
1.2 À = 1.00
1.0
700 800 900 1000 1100 1200
Model: 2D Elastoplastic
Failure Criteria: Mohr Coulomb
400
Ò÷ grad = 0.95 psi/ft
Òüîâù grad = 0.95 psi/ft
Òüîêï grad = 0.95 psi/ft
Òüîêï Azi = 45 °
300
Pó grad = 0.40 psi/ft
Pó = 1579 psi
cñ = 650 psi
200 Åñ = 40 °
có = 145 psi
Åó = 35 °
E = 0.15 E+06 psi
Í = 0.30
100 À = 1.00
0
0 500 1000 1500
Residual Friction 25 35 30 °
Angle
Apparent Capillary
10 0 30 psi
Cohesion
Base Case R/rw = 1.58
CASE PROFILE
Analysis Type Hole Collapse and Fracture Breakdown / Lost Circulation
Analysis Mode Calculate required bottomhole pressures or mud densities to achieve tolerable borehole
instability risks.
Date and Time 26-Nov-2008, 14:49
Filename: Case 5 - Build and Horizontal Well EMW Window.bhs 26-Nov-2008 14:49
Case Title: Extended Reach Horizontal Well 2 of 9 ABC Oil & Gas
Depth Incl. Azimuth Depth N/S E/W DLS Depth Incl. Azimuth Depth N/S E/W DLS
ft KB MD ° ° ft KB TVD Distance Distance deg/100 ft KB MD ° ° ft KB TVD Distance Distance deg/100
ft ft ft ft ft ft
7545.84 40.00 22.00 6523.07 2905.15 1173.76 0.0 15747.84 90.00 80.00 10023.57 6325.01 6179.20 0.0
8202.00 40.00 22.00 7025.72 3296.21 1331.76 0.0 15944.69 90.00 80.00 10023.57 6359.19 6373.06 0.0
8858.16 40.00 22.00 7528.37 3687.27 1489.75 0.0
9514.32 40.00 22.00 8031.02 4078.33 1647.75 0.0
BOREHOLE PROPERTIES
Bottom Borehole
No. Depth Diameter
ft KB MD in
1 820.00 36.00
2 1970.00 16.00
3 9840.00 12.25
4 15944.69 8.50
TUBULAR PROPERTIES
No. Tubular Type Top Depth Btm Depth ID OD
ft KB MD ft KB MD in in
1 surface casing 0.00 804.00 28.75 30.00
2 casing 0.00 1952.00 12.40 13.40
3 casing 0.00 9826.00 8.65 9.60
STRATIGRAPHIC UNITS
No. Stratigraphic Unit Lithology Pattern Top Depth Btm Depth Top Depth Btm Depth
ft KB MD ft KB MD ft KB TVD ft KB TVD
1 Tertiary Mudstone mudstone 3280.80 7500.10 3216.81 6488.03
2 Cretaceous Shale shale 7500.10 11318.75 6488.03 9387.62
3 Reservoir Unit A sandstone 11318.75 11909.30 9387.62 9725.90
4 Reservoir Unit B limestone 11909.30 12303.00 9725.90 9887.48
5 Reservoir Unit C sandstone 12303.00 15944.60 9887.48 10023.57
INITIAL FORMATION STRESS AND PRESSURE GRADIENT DATA
No. Stratigraphic Unit Top Top SHmin Sv Sv SHmax SHmax SHmin SHmin Pr Pr
Depth Depth Azi. Grad Grad Grad Grad Grad Grad Grad Grad
ft KB MD ft KB TVD (Top) (Btm) (Top) (Btm) (Top) (Btm) (Top) (Btm)
psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft
1 Tertiary Mudstone 3280.80 3216.81 40 0.93 0.99 0.77 0.84 0.69 0.75 0.44 0.44
2 Cretaceous Shale 7500.10 6488.03 40 0.99 1.02 0.84 0.87 0.75 0.79 0.44 0.44
3 Reservoir Unit A 11318.75 9387.62 40 1.02 1.02 0.87 0.88 0.79 0.80 0.44 0.44
4 Reservoir Unit B 11909.30 9725.90 40 1.02 1.02 0.90 0.90 0.81 0.81 0.44 0.44
5 Reservoir Unit C 12303.00 9887.48 40 1.02 1.02 0.88 0.88 0.80 0.80 0.44 0.44
BASIC ROCK MECHANICAL PROPERTIES
No. Stratigraphic Unit Lithology Top Top Peak Peak Poisson's Young's Tensile Biot's
Depth Depth Cohesion Friction Ratio Modulus Strength Coefficient
ft KB MD ft KB psi Angle E+06 psi psi
TVD
1 Tertiary Mudstone mudstone 3280.80 3216.81 350 30 0.35 0.29 15 1.00
2 Cretaceous Shale shale 7500.10 6488.03 400 35 0.33 0.40 20 1.00
3 Reservoir Unit A sandstone 11318.75 9387.62 580 30 0.30 0.73 45 1.00
4 Reservoir Unit B limestone 11909.30 9725.90 1500 40 0.20 0.58 100 0.80
5 Reservoir Unit C sandstone 12303.00 9887.48 725 35 0.25 0.44 75 0.90
FLUID PENETRATION PROPERTIES
No. Stratigraphic Unit Lithology Top Top Capillary Filter Cake Fluid
Depth Depth Threshold or Wall Penetration
ft KB MD ft KB TVD Pressure Coating Coefficient
psi Efficiency
Filename: Case 5 - Build and Horizontal Well EMW Window.bhs 26-Nov-2008 14:49
Case Title: Extended Reach Horizontal Well 3 of 9 ABC Oil & Gas
Borehole Collapse Model Type: Default (3D Linear Elastic Model with Mohr-Coulomb Failure Criterion)
Calibrate Borehole Collapse Pressure Based on the following Well:
Inclination 0 °
Azimuth 0 °
Vertical Depth 6400.00 ft
Calibration Method: Linear Multiplier
Borehole Collapse Type: EMW at Borehole Collapse
Calculated Calibration Ratio 0.73
Unit No. 2 - Stratigraphic Unit: Cretaceous Shale, Lithology: shale (Top MD = 7500.10 ft ; Top TVD = 6488.03 ft)
Fracture Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
Borehole Collapse Model Type: Default (3D Linear Elastic Model with Mohr-Coulomb Failure Criterion)
Calibrate Borehole Collapse Pressure Based on the following Well:
Inclination 0 °
Azimuth 0 °
Vertical Depth 9000.00 ft
Calibration Method: Linear Multiplier
Borehole Collapse Type: EMW at Borehole Collapse
Calculated Calibration Ratio 0.72
Unit No. 3 - Stratigraphic Unit: Reservoir Unit A, Lithology: sandstone (Top MD = 11318.75 ft ; Top TVD = 9387.62 ft)
Fracture Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
Borehole Collapse Model Type: Default (3D Linear Elastic Model with Mohr-Coulomb Failure Criterion)
Calibrate Borehole Collapse Pressure Based on the following Well:
Inclination 0 °
Azimuth 0 °
Vertical Depth 9515.00 ft
Calibration Method: Linear Multiplier
Borehole Collapse Type: EMW at Borehole Collapse
Calculated Calibration Ratio 0.73
Unit No. 4 - Stratigraphic Unit: Reservoir Unit B, Lithology: limestone (Top MD = 11909.30 ft ; Top TVD = 9725.90 ft)
Fracture Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
Borehole Collapse Model Type: Default (3D Linear Elastic Model with Mohr-Coulomb Failure Criterion)
* No Model Options selected.
Unit No. 5 - Stratigraphic Unit: Reservoir Unit C, Lithology: sandstone (Top MD = 12303.00 ft ; Top TVD = 9887.48 ft)
Fracture Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
Borehole Collapse Model Type: Default (3D Linear Elastic Model with Mohr-Coulomb Failure Criterion)
In-situ Stress Regime: Passive Basin
Depleted Formation Pressure Gradient 0.37 psi/ft
Depleted Formation Pressure 3658 psi
In-situ stresses after depletion:
Max Horz Stress Gradient 0.84 psi/ft
Max Horz Stress 8286 psi
Min Horz Stress Gradient 0.76 psi/ft
Min Horz Stress 7495 psi
Filename: Case 5 - Build and Horizontal Well EMW Window.bhs 26-Nov-2008 14:49
Case Title: Extended Reach Horizontal Well 4 of 9 ABC Oil & Gas
CASE PROFILE
Analysis Type Hole Collapse and Fracture Breakdown / Lost Circulation
Analysis Mode Calculate required bottomhole pressures or mud densities to achieve tolerable borehole
instability risks.
Date and Time 26-Nov-2008, 14:49
3280.83 3216.84 2443 0.76 14.61 5691.83 5102.82 4246 0.83 16.01
3582.18 3475.83 2698 0.78 14.94 5993.21 5333.69 4471 0.84 16.13
3883.56 3717.42 2951 0.79 15.28 6294.59 5564.56 4700 0.84 16.25
4184.94 3948.47 3161 0.80 15.40 6595.96 5795.42 4932 0.85 16.38
4486.31 4179.34 3372 0.81 15.52 6897.34 6026.29 5167 0.86 16.50
4787.69 4410.21 3586 0.81 15.65 7198.72 6257.16 5404 0.86 16.62
5089.07 4641.08 3803 0.82 15.77 7500.07 6488.01 5645 0.87 16.74
5390.45 4871.95 4023 0.83 15.89
Unit No. 2 - Stratigraphic Unit: Cretaceous Shale, Lithology: shale, Model Type: Default (3D Linear Elastic Model for Tensile Fracture
Breakdown)
7500.13 6488.06 5743 0.89 17.03 9682.19 8159.61 7554 0.93 17.81
7772.86 6696.98 5962 0.89 17.13 9954.95 8368.56 7790 0.93 17.91
8045.62 6905.93 6183 0.90 17.23 10227.71 8577.50 8027 0.94 18.01
8318.38 7114.87 6407 0.90 17.33 10500.47 8786.45 8267 0.94 18.11
8591.14 7323.82 6632 0.91 17.42 10773.23 8995.40 8509 0.95 18.20
8863.90 7532.77 6859 0.91 17.52 11045.99 9198.22 9094 0.99 19.02
9136.67 7741.71 7089 0.92 17.62 11318.72 9387.60 9718 1.04 19.92
9409.43 7950.66 7320 0.92 17.72
Unit No. 3 - Stratigraphic Unit: Reservoir Unit A, Lithology: sandstone, Model Type: Default (3D Linear Elastic Model for Tensile Fracture
Breakdown)
Filename: Case 5 - Build and Horizontal Well EMW Window.bhs 26-Nov-2008 14:49
Case Title: Extended Reach Horizontal Well 5 of 9 ABC Oil & Gas
Depth Depth BHP Pressure ECD Depth Depth BHP Pressure ECD
ft KB MD ft KB TVD psi Gradient ppg ft KB MD ft KB TVD psi Gradient ppg
psi/ft psi/ft
11318.78 9387.65 11462 1.22 23.49 11656.21 9592.22 11426 1.19 22.92
11360.93 9414.18 11468 1.22 23.44 11698.39 9616.42 11417 1.19 22.85
11403.11 9440.73 11468 1.21 23.38 11740.57 9638.97 11404 1.18 22.77
11445.30 9467.29 11465 1.21 23.30 11782.75 9660.70 11389 1.18 22.68
11487.48 9493.85 11459 1.21 23.23 11824.94 9682.43 11373 1.17 22.60
11529.66 9519.62 11451 1.20 23.15 11867.12 9704.17 11357 1.17 22.52
11571.84 9543.82 11444 1.20 23.07 11909.27 9725.88 11342 1.17 22.44
11614.03 9568.02 11435 1.20 23.00
Unit No. 4 - Stratigraphic Unit: Reservoir Unit B, Lithology: limestone, Model Type: Default (3D Linear Elastic Model for Tensile Fracture
Breakdown)
11909.33 9725.91 12735 1.31 25.20 12134.27 9824.27 12587 1.28 24.65
11937.42 9738.44 12718 1.31 25.13 12162.39 9834.80 12568 1.28 24.59
11965.54 9750.99 12699 1.30 25.06 12190.51 9845.34 12549 1.27 24.53
11993.66 9763.54 12680 1.30 24.99 12218.64 9855.88 12530 1.27 24.46
12021.79 9776.09 12662 1.30 24.92 12246.76 9866.41 12511 1.27 24.40
12049.91 9788.64 12643 1.29 24.85 12274.88 9876.95 12492 1.26 24.34
12078.03 9801.18 12625 1.29 24.79 12302.97 9887.47 12473 1.26 24.27
12106.15 9813.73 12606 1.28 24.72
Unit No. 5 - Stratigraphic Unit: Reservoir Unit C, Lithology: sandstone, Model Type: Default (3D Linear Elastic Model for Tensile Fracture
Breakdown)
12303.03 9887.49 11050 1.12 21.51 14383.91 10023.57 10353 1.03 19.87
12563.11 9960.37 10794 1.08 20.85 14644.03 10023.57 10353 1.03 19.87
12823.23 10005.80 10573 1.06 20.33 14904.14 10023.57 10353 1.03 19.87
13083.34 10023.26 10358 1.03 19.89 15164.26 10023.57 10353 1.03 19.87
13343.46 10023.57 10353 1.03 19.87 15424.37 10023.57 10353 1.03 19.87
13603.57 10023.57 10353 1.03 19.87 15684.49 10023.57 10353 1.03 19.87
13863.69 10023.57 10353 1.03 19.87 15944.57 10023.57 10353 1.03 19.87
14123.80 10023.57 10353 1.03 19.87
3280.83 3216.84 1439 0.45 8.61 5691.83 5102.82 2865 0.56 10.81
3582.18 3475.83 1664 0.48 9.21 5993.21 5333.69 3032 0.57 10.94
3883.56 3717.42 1889 0.51 9.78 6294.59 5564.56 3200 0.58 11.07
4184.94 3948.47 2050 0.52 9.99 6595.96 5795.42 3369 0.58 11.19
4486.31 4179.34 2211 0.53 10.18 6897.34 6026.29 3540 0.59 11.30
4787.69 4410.21 2372 0.54 10.35 7198.72 6257.16 3711 0.59 11.41
5089.07 4641.08 2536 0.55 10.51 7500.07 6488.01 3884 0.60 11.52
5390.45 4871.95 2700 0.55 10.66
Unit No. 2 - Stratigraphic Unit: Cretaceous Shale, Lithology: shale, Model Type: Default (3D Linear Elastic Model with Mohr-Coulomb Failure
Criterion)
7500.13 6488.06 3640 0.56 10.80 9682.19 8159.61 4802 0.59 11.33
7772.86 6696.98 3783 0.56 10.87 9954.95 8368.56 4950 0.59 11.38
8045.62 6905.93 3927 0.57 10.94 10227.71 8577.50 5098 0.59 11.44
8318.38 7114.87 4072 0.57 11.01 10500.47 8786.45 5247 0.60 11.49
8591.14 7323.82 4217 0.58 11.08 10773.23 8995.40 5396 0.60 11.54
8863.90 7532.77 4363 0.58 11.14 11045.99 9198.22 5587 0.61 11.69
9136.67 7741.71 4509 0.58 11.21 11318.72 9387.60 5818 0.62 11.93
9409.43 7950.66 4655 0.59 11.27
Unit No. 3 - Stratigraphic Unit: Reservoir Unit A, Lithology: sandstone, Model Type: Default (3D Linear Elastic Model with Mohr-Coulomb
Failure Criterion)
11318.78 9387.65 4730 0.50 9.70 11656.21 9592.22 4941 0.52 9.91
11360.93 9414.18 4757 0.51 9.72 11698.39 9616.42 4966 0.52 9.94
11403.11 9440.73 4784 0.51 9.75 11740.57 9638.97 4990 0.52 9.96
Filename: Case 5 - Build and Horizontal Well EMW Window.bhs 26-Nov-2008 14:49
Case Title: Extended Reach Horizontal Well 6 of 9 ABC Oil & Gas
Depth Depth BHP Pressure EMW Depth Depth BHP Pressure EMW
ft KB MD ft KB TVD psi Gradient ppg ft KB MD ft KB TVD psi Gradient ppg
psi/ft psi/ft
11445.30 9467.29 4811 0.51 9.78 11782.75 9660.70 5013 0.52 9.99
11487.48 9493.85 4838 0.51 9.81 11824.94 9682.43 5036 0.52 10.01
11529.66 9519.62 4865 0.51 9.83 11867.12 9704.17 5059 0.52 10.03
11571.84 9543.82 4890 0.51 9.86 11909.27 9725.88 5082 0.52 10.05
11614.03 9568.02 4916 0.51 9.89
Unit No. 4 - Stratigraphic Unit: Reservoir Unit B, Lithology: limestone, Model Type: Default (3D Linear Elastic Model with Mohr-Coulomb
Failure Criterion)
11909.33 9725.91 4694 0.48 9.29 12134.27 9824.27 4832 0.49 9.46
11937.42 9738.44 4712 0.48 9.31 12162.39 9834.80 4846 0.49 9.48
11965.54 9750.99 4730 0.49 9.34 12190.51 9845.34 4860 0.49 9.50
11993.66 9763.54 4748 0.49 9.36 12218.64 9855.88 4874 0.49 9.52
12021.79 9776.09 4766 0.49 9.38 12246.76 9866.41 4888 0.50 9.53
12049.91 9788.64 4784 0.49 9.40 12274.88 9876.95 4901 0.50 9.55
12078.03 9801.18 4801 0.49 9.43 12302.97 9887.47 4915 0.50 9.56
12106.15 9813.73 4818 0.49 9.45
Unit No. 5 - Stratigraphic Unit: Reservoir Unit C, Lithology: sandstone, Model Type: Default (3D Linear Elastic Model with Mohr-Coulomb
Failure Criterion)
12303.03 9887.49 6050 0.61 11.77 14383.91 10023.57 6295 0.63 12.09
12563.11 9960.37 6174 0.62 11.93 14644.03 10023.57 6295 0.63 12.09
12823.23 10005.80 6256 0.63 12.03 14904.14 10023.57 6295 0.63 12.09
13083.34 10023.26 6295 0.63 12.08 15164.26 10023.57 6295 0.63 12.09
13343.46 10023.57 6295 0.63 12.09 15424.37 10023.57 6295 0.63 12.09
13603.57 10023.57 6295 0.63 12.09 15684.49 10023.57 6295 0.63 12.09
13863.69 10023.57 6295 0.63 12.09 15944.57 10023.57 6295 0.63 12.09
14123.80 10023.57 6295 0.63 12.09
Filename: Case 5 - Build and Horizontal Well EMW Window.bhs 26-Nov-2008 14:49
Case Title: Extended Reach Horizontal Well 7 of 9 ABC Oil & Gas
Äçä = 0.30
11 cñ = 400 psi
Åñ = 35 °
E = 0.40 E+06 psi
Í = 0.33
À = 1.00
10
9
0 10 20 30 40 50 60 70 80 90
N
Analysis Mode: Pressure 13.02
Analysis Type: Drilling Collapse
Äçä = 0.30
cñ = 400 psi 10.56
Åñ = 35 °
E = 0.40 E+06 psi ppg
Í = 0.33
À = 1.00
¤ Base Case
S Òüîêï
Filename: Case 5 - Build and Horizontal Well EMW Window.bhs 26-Nov-2008 14:49
Case Title: Extended Reach Horizontal Well ABC Oil & Gas
8 of 9
2000
3000 11.22
4000
10.35
5000
TVD (ft)
6000 9.48
7000
8.61
8000 ppg
9000
10000
11000
12000
0 2000 4000 6000 8000 10000
-1000 1000 3000 5000 7000 9000
2000
3000
Tertiary Mudstone
Measured Depth (ft KB)
4000
5000
6000
7000
9000
10000
11000
Reservoir Unit A
12000 Reservoir Unit B
Reservoir Unit C
13000
14000
15000
Filename: Case 5 - Build and Horizontal Well EMW Window.bhs 26-Nov-2008 14:49
Case Title: Extended Reach Horizontal Well 9 of 9 ABC Oil & Gas
11.22
5000 N
North/South Axis (ft)
10.35
4000 N
9.48
3000 N
2000 N 8.61
ppg
1000 N
1000 S
0 2000 E 4000 E 6000 E
1000 W 1000 E 3000 E 5000 E 7000 E
Filename: Case 5 - Build and Horizontal Well EMW Window.bhs 26-Nov-2008 14:49
Case Title: Effects of Cool Fluid Injection on Fracture Breakdown 1 of 6 CRX Oil and Gas
CASE PROFILE
Analysis Type Fracture Breakdown / Lost Circulation
Analysis Mode Calculate required bottomhole pressures or mud densities to achieve tolerable borehole
instability risks.
Date and Time 26-Nov-2008, 14:27
BOREHOLE PROPERTIES
Bottom Borehole
No. Depth Diameter
ft KB MD in
1 1968.50 16.00
2 8038.00 12.25
3 9317.47 8.50
Filename: Case 6 - Inclined Well with Cool Fluid Injection.bhs 26-Nov-2008 14:27
Case Title: Effects of Cool Fluid Injection on Fracture Breakdown 2 of 6 CRX Oil and Gas
TUBULAR PROPERTIES
No. Tubular Type Top Depth Btm Depth ID OD
ft KB MD ft KB MD in in
1 casing 0.00 1903.00 12.40 13.35
2 casing 0.00 7972.00 8.50 9.40
STRATIGRAPHIC UNITS
No. Stratigraphic Unit Lithology Pattern Top Depth Btm Depth Top Depth Btm Depth
ft KB MD ft KB MD ft KB TVD ft KB TVD
1 Pliocene - Oligocene mudstone 1640.00 4921.00 1639.99 4910.96
2 Eocene shale 4921.00 8005.00 4910.96 7851.68
3 Forties Formation sandstone 8005.00 9317.50 7851.68 8882.77
INITIAL FORMATION STRESS AND PRESSURE GRADIENT DATA
No. Stratigraphic Unit Top Top SHmin Sv Sv SHmax SHmax SHmin SHmin Pr Pr
Depth Depth Azi. Grad Grad Grad Grad Grad Grad Grad Grad
ft KB MD ft KB TVD (Top) (Btm) (Top) (Btm) (Top) (Btm) (Top) (Btm)
psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft
1 Pliocene - Oligocene 1640.00 1639.99 10 0.97 0.97 0.84 0.84 0.71 0.71 0.44 0.44
2 Eocene 4921.00 4910.96 13 0.99 0.99 0.88 0.88 0.75 0.75 0.44 0.44
3 Forties Formation 8005.00 7851.68 15 0.99 0.99 0.91 0.91 0.80 0.80 0.44 0.44
BASIC ROCK MECHANICAL PROPERTIES
No. Stratigraphic Unit Lithology Top Top Poisson's Young's Tensile Biot's
Depth Depth Ratio Modulus Strength Coefficient
ft KB MD ft KB E+06 psi psi
TVD
1 Pliocene - Oligocene mudstone 1640.00 1639.99 0.35 0.22 30 1.00
2 Eocene shale 4921.00 4910.96 0.30 0.44 75 1.00
3 Forties Formation sandstone 8005.00 7851.68 0.25 0.70 100 1.00
FLUID PENETRATION PROPERTIES
No. Stratigraphic Unit Lithology Top Top Fluid
Depth Depth Penetration
ft KB MD ft KB TVD Coefficient
Unit No. 2 - Stratigraphic Unit: Eocene, Lithology: shale (Top MD = 4921.00 ft ; Top TVD = 4910.96 ft)
Fracture Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
* No Model Options selected.
Unit No. 3 - Stratigraphic Unit: Forties Formation, Lithology: sandstone (Top MD = 8005.00 ft ; Top TVD = 7851.68 ft)
Fracture Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
Steady-state Thermal Effects:
Thermal Expansion Coefficient 20.0 E-06/°F
Temperature Data:
Surface Temperature 41 °F
Wellbore Temperature (Top) 110 °F
Wellbore Temperature (Btm) 110 °F
Formation Temperature (Top) 140 °F
Formation Temperature (Btm) 140 °F
Filename: Case 6 - Inclined Well with Cool Fluid Injection.bhs 26-Nov-2008 14:27
Case Title: Effects of Cool Fluid Injection on Fracture Breakdown 3 of 6 CRX Oil and Gas
CASE PROFILE
Analysis Type Fracture Breakdown / Lost Circulation
Analysis Mode Calculate required bottomhole pressures or mud densities to achieve tolerable borehole
instability risks.
Date and Time 26-Nov-2008, 14:27
1640.03 1640.02 1356 0.83 15.91 3417.21 3413.63 2790 0.82 15.73
1776.71 1776.68 1467 0.83 15.88 3553.92 3549.82 2900 0.82 15.72
1913.42 1913.34 1577 0.82 15.86 3690.63 3686.01 3011 0.82 15.72
2050.13 2049.97 1688 0.82 15.84 3827.33 3822.19 3121 0.82 15.71
2186.83 2186.55 1798 0.82 15.83 3964.04 3958.38 3231 0.82 15.71
2323.54 2323.10 1909 0.82 15.81 4100.75 4094.57 3341 0.82 15.70
2460.25 2459.55 2019 0.82 15.80 4237.46 4230.76 3452 0.82 15.70
2596.96 2595.99 2129 0.82 15.78 4374.17 4366.93 3562 0.82 15.69
2733.67 2732.37 2240 0.82 15.77 4510.88 4503.01 3671 0.82 15.69
2870.38 2868.74 2350 0.82 15.76 4647.58 4639.08 3781 0.82 15.68
3007.08 3005.02 2460 0.82 15.75 4784.29 4775.16 3891 0.81 15.68
3143.79 3141.25 2570 0.82 15.74 4920.97 4910.93 3999 0.81 15.67
3280.50 3277.44 2680 0.82 15.73
Unit No. 2 - Stratigraphic Unit: Eocene, Lithology: shale, Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
4921.03 4910.99 4354 0.89 17.06 6591.50 6540.12 5731 0.88 16.86
5049.50 5038.50 4464 0.89 17.05 6720.00 6662.89 5834 0.88 16.85
5178.00 5166.03 4575 0.89 17.04 6848.50 6784.01 5932 0.87 16.83
5306.50 5293.45 4684 0.88 17.03 6977.00 6905.13 6030 0.87 16.81
5435.00 5419.76 4791 0.88 17.01 7105.50 7026.25 6129 0.87 16.78
5563.50 5546.07 4897 0.88 16.99 7234.00 7146.55 6225 0.87 16.76
5692.00 5672.38 5003 0.88 16.97 7362.50 7265.24 6319 0.87 16.74
5820.50 5797.90 5107 0.88 16.95 7491.00 7383.52 6414 0.87 16.72
Filename: Case 6 - Inclined Well with Cool Fluid Injection.bhs 26-Nov-2008 14:27
Case Title: Effects of Cool Fluid Injection on Fracture Breakdown 4 of 6 CRX Oil and Gas
Depth Depth BHP Pressure ECD Depth Depth BHP Pressure ECD
ft KB MD ft KB TVD psi Gradient ppg ft KB MD ft KB TVD psi Gradient ppg
psi/ft psi/ft
5949.00 5922.29 5212 0.88 16.94 7619.50 7501.80 6509 0.87 16.70
6077.50 6046.69 5316 0.88 16.92 7748.00 7618.77 6604 0.87 16.68
6206.00 6171.08 5420 0.88 16.90 7876.50 7735.23 6699 0.87 16.66
6334.50 6294.36 5524 0.88 16.89 8004.97 7851.65 6793 0.87 16.65
6463.00 6417.24 5628 0.88 16.88
Unit No. 3 - Stratigraphic Unit: Forties Formation, Lithology: sandstone, Model Type: Default (3D Linear Elastic Model for Tensile Fracture
Breakdown)
8005.03 7851.71 6786 0.86 16.63 8715.94 8455.72 7193 0.85 16.37
8059.69 7901.24 6830 0.86 16.63 8770.63 8495.64 7217 0.85 16.35
8114.38 7950.80 6873 0.86 16.63 8825.31 8534.73 7248 0.85 16.34
8169.06 7999.70 6915 0.86 16.63 8880.00 8573.41 7290 0.85 16.36
8223.75 8048.43 6956 0.86 16.63 8934.69 8612.08 7325 0.85 16.37
8278.44 8097.15 6997 0.86 16.63 8989.38 8650.75 7360 0.85 16.37
8333.13 8145.88 7038 0.86 16.63 9044.06 8689.42 7394 0.85 16.38
8387.81 8191.62 7045 0.86 16.55 9098.75 8728.09 7430 0.85 16.38
8442.50 8237.74 7081 0.86 16.54 9153.44 8766.77 7467 0.85 16.39
8497.19 8283.86 7117 0.86 16.53 9208.13 8805.44 7505 0.85 16.40
8551.88 8328.94 7139 0.86 16.49 9262.81 8844.11 7544 0.85 16.41
8606.56 8372.44 7154 0.85 16.44 9317.47 8882.77 7581 0.85 16.42
8661.25 8414.77 7173 0.85 16.40
Filename: Case 6 - Inclined Well with Cool Fluid Injection.bhs 26-Nov-2008 14:27
Case Title: Effects of Cool Fluid Injection on Fracture Breakdown 5 of 6 CRX Oil and Gas
0.82
0.80
0 30 60 90 120 150 180 210 240 270 300 330 360
N
Analysis Mode: Pressure 9750
Analysis Type: Fracture Breakdown
Äçñ = 0.80
cñ = 0 psi 7192
Åñ = 0 °
E = 0.70 E+06 psi psi
Í = 0.25
À = 1.00
¤ Base Case
S Òüîêï
Filename: Case 6 - Inclined Well with Cool Fluid Injection.bhs 26-Nov-2008 14:27
Case Title: Effects of Cool Fluid Injection on Fracture Breakdown 6 of 6 CRX Oil and Gas
1000
0.87
2000
3000 0.85
TVD (ft)
4000
0.83
5000
0.81
6000
psi/f t
7000
8000
9000
10000
-4000 -3000 -2000 -1000 0 1000 2000 3000 4000 5000 6000
0.8
0 30 60 90 120 150 180 210 240 270 300 330 360
Filename: Case 6 - Inclined Well with Cool Fluid Injection.bhs 26-Nov-2008 14:27
Case Title: Shear Slip Analysis for Weak Natural Fracture Planes 1 of 3 Black Oil International
CASE PROFILE
Analysis Type Fracture Breakdown / Lost Circulation
Analysis Mode Calculate required bottomhole pressures or mud densities to achieve tolerable borehole
instability risks.
Fracture Model Type 3D Slip Analysis on a Weak Discontinuity
Date and Time 26-Nov-2008, 14:28
Case Title Shear Slip Analysis for Weak Natural Fracture Planes
Well Name Water Disposal Well Ls-2
Company AGI
Client Black Oil International
User Name B. Smith
Unit Type US Oilfield
Comments Analysis to determine if elevated pore pressures during water injection will cause shear failure
on weak natural fractures or bedding planes.
CASE PREFERENCES
Well Profile Resolution Low
Cross-Sections
Display Tubulars & Rock Deformations Yes
Display Input Parameters Yes
Display Breakout Angle Yes
Polar Plots
Fracture Breakdown Pressure Limit Set No
Color/Grayscale Color
Graphic Resolution Medium
Tolerance Angles
2D Elastoplastic Model Inclination Tolerance 25°
Vertical Well Inclination Tolerance 15°
Horizontal Well Inclination Tolerance 15°
WELL DATA
Offshore/Onshore Onshore
KB Height (relative to ground level) 16.40 ft
Well Survey Vertical Well
FORMATION DATA
Stratigraphic Unit Nisku
Lithology dolomite
SINGLE DEPTH ANALYSIS DATA
Calculation Depth (TVD) 9900.00 ft KB
Borehole Diameter 8.50 in
Azimuth 0.00 °
Inclination 0.00 °
IN-SITU STRESSES AND PRESSURES
Azimuth of Minimum Horizontal Stress 135 °
Vertical Stress Gradient 1.15 psi/ft
Vertical Stress 11366 psi
Max Horz Stress Gradient 1.17 psi/ft
Max Horz Stress 11583 psi
Min Horz Stress Gradient 0.95 psi/ft
Min Horz Stress 9405 psi
Formation Pressure Gradient 0.36 psi/ft
Formation Pressure 3564 psi
ROCK MECHANICAL PROPERTIES
Peak Cohesion 0 psi
Peak Friction Angle 0 °
Poisson's Ratio 0.22
9000
8500
0 10 20 30 40 50 60 70 80 90
N
Analysis Mode: Pressure 22.70
Analysis Type: Fracture Breakdown
16.92
Strength of Planes of Weakness
Base Case ¤ ppg
cç = 75 psi
Åç = 35 °
Dip = 70 ° ¤ Base Case
Dip Dir. = 10 °
S Òüîêï
Lower Hemisphere Plot Shows Poles to Planes of Weakness
CASE PROFILE
Analysis Type Hole Collapse While Drilling
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Borehole Collapse Model Type
Date and Time 26-Nov-2008, 14:29
2.2
2.0
Normalized Washout Zone
1.8
Radius
1.6
1.4
1.2
1.0
0 50 100 150
Time (hours)
5000
4500
4000
Reynolds Number
3500
3000
2500
2000
1500
1000
0 50 100 150
Time (hours)
CASE PROFILE
Analysis Type Hole Collapse While Drilling
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Date and Time 26-Nov-2008, 14:30
BOREHOLE PROPERTIES
Bottom Borehole
No. Depth Diameter
m KB MD mm
1 150.00 438.0
TUBULAR PROPERTIES
* No tubular data were entered for this case.
STRATIGRAPHIC UNITS
No. Stratigraphic Unit Lithology Pattern Top Depth Btm Depth Top Depth Btm Depth
m KB MD m KB MD m KB TVD m KB TVD
1 Holocene unconsolidated 100.00 110.00 100.00 110.00
sediment
2 Iperk Sequence unconsolidated 110.00 150.00 110.00 150.00
sediment
unconsolidated
1
Holocene sediment 100.00 100.00 0.0 0.30
unconsolidated
2
Iperk Sequence sediment 110.00 110.00 0.0 0.30
MODEL OPTIONS & PROPERTIES
Unit No. 1 - Stratigraphic Unit: Holocene, Lithology: unconsolidated sediment (Top MD = 100.00 m ; Top TVD = 100.00 m)
Borehole Collapse Model Type:
Drilling Mud Properties for 2D Permafrost Model:
Mud Thermal Conductivity 0.565 W/m·°C
Mud Specific Heat 4.120 kJ/kg·°C
Mud Temperature 1 °C
Mud Density 1000 kg/m³
Circulation Rate 0.002 m³ / s
Drill Pipe Outer Diameter 127.0 mm
Circulation Time 140.00 hours
Formation Porosity 50 %
Drill Bit Diameter 438.0 mm
Newtonian Fluid Rheological Model:
Dynamic Viscosity 1.5 mPa·s
Unit No. 2 - Stratigraphic Unit: Iperk Sequence, Lithology: unconsolidated sediment (Top MD = 110.00 m ; Top TVD = 110.00 m)
Borehole Collapse Model Type:
Drilling Mud Properties for 2D Permafrost Model:
Mud Thermal Conductivity 0.565 W/m·°C
Mud Specific Heat 4.187 kJ/kg·°C
Mud Temperature 2 °C
Mud Density 1000 kg/m³
Circulation Rate 0.002 m³ / s
Drill Pipe Outer Diameter 127.0 mm
Circulation Time 110.00 hours
Formation Porosity 44 %
Drill Bit Diameter 438.0 mm
Newtonian Fluid Rheological Model:
Dynamic Viscosity 1.5 mPa·s
CASE PROFILE
Analysis Type Hole Collapse While Drilling
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Date and Time 26-Nov-2008, 14:30
Unit No. 2 - Stratigraphic Unit: Iperk Sequence, Lithology: unconsolidated sediment, Model Type: 2D Permafrost Hole Enlargement
2.2
2.0
Normalized Washout Zone
1.8
Radius
1.6
1.4
1.2
1.0
0 50 100 150
Time (hours)
5000
4500
4000
Reynolds Number
3500
3000
2500
2000
1500
1000
0 50 100 150
Time (hours)
1.8
1.6
Radius
1.4
1.2
1.0
0 50 100 150
Time (hours)
CASE PROFILE
Analysis Type Hole Collapse While Drilling
Analysis Mode Calculate required bottomhole pressures or mud densities to achieve tolerable borehole
instability risks.
Date and Time 26-Nov-2008, 14:35
BOREHOLE PROPERTIES
Bottom Borehole
No. Depth Diameter
ft KB MD in
1 650.00 16.00
2 2500.00 12.25
3 12900.00 8.50
4 15771.43 6.00
TUBULAR PROPERTIES
No. Tubular Type Top Depth Btm Depth ID OD
ft KB MD ft KB MD in in
1 casing 0.00 650.00 11.80 13.00
2 casing 0.00 2500.00 8.50 9.50
3 liner 2450.00 12900.00 6.00 7.00
STRATIGRAPHIC UNITS
No. Stratigraphic Unit Lithology Pattern Top Depth Btm Depth Top Depth Btm Depth
ft KB MD ft KB MD ft KB TVD ft KB TVD
1 Rotliegendes 1 sandstone 13000.00 14100.00 7631.42 7809.56
2 Lower Cretaceous shale w/ siltstone 14100.00 14400.00 7809.56 7858.14
3 Rotliegendes 2 sandstone 14400.00 15771.00 7858.14 8080.17
INITIAL FORMATION STRESS AND PRESSURE GRADIENT DATA
No. Stratigraphic Unit Top Top SHmin Sv Sv SHmax SHmax SHmin SHmin Pr Pr
Depth Depth Azi. Grad Grad Grad Grad Grad Grad Grad Grad
ft KB MD ft KB TVD (Top) (Btm) (Top) (Btm) (Top) (Btm) (Top) (Btm)
psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft
1 Rotliegendes 1 13000.00 7631.42 162 0.90 0.90 0.72 0.72 0.72 0.72 0.52 0.52
2 Lower Cretaceous 14100.00 7809.56 162 0.90 0.90 0.72 0.72 0.72 0.72 0.52 0.52
3 Rotliegendes 2 14400.00 7858.14 162 0.90 0.90 0.72 0.72 0.72 0.72 0.52 0.52
BASIC ROCK MECHANICAL PROPERTIES
No. Stratigraphic Unit Lithology Top Top Peak Peak Poisson's Young's Biot's
Depth Depth Cohesion Friction Ratio Modulus Coefficient
ft KB MD ft KB psi Angle E+06 psi
TVD
1 Rotliegendes 1 sandstone 13000.00 7631.42 400 45 0.32 0.20 1.00
shale w/
2
Lower Cretaceous siltstone 14100.00 7809.56 300 35 0.35 0.15 1.00
3 Rotliegendes 2 sandstone 14400.00 7858.14 800 45 0.30 0.40 1.00
FLUID PENETRATION PROPERTIES
No. Stratigraphic Unit Lithology Top Top Capillary Filter Cake
Depth Depth Threshold or Wall
ft KB MD ft KB TVD Pressure Coating
psi Efficiency
Unit No. 2 - Stratigraphic Unit: Lower Cretaceous, Lithology: shale w/ siltstone (Top MD = 14100.00 ft ; Top TVD = 7809.56 ft)
Borehole Collapse Model Type: Default (3D Linear Elastic Model with Mohr-Coulomb Failure Criterion)
Maximum Tolerable Wellbore Instability Failure Criteria: Percentage of the Hole Circumference
Tolerable Yield Circumf., % 60.00
Bedding Dip 20 °
Bedding Dip Direction 340 °
Bedding Cohesion 75 psi
Bedding Friction Angle 30 °
Calibrate Borehole Collapse Pressure Based on the following Well:
Inclination 0 °
Azimuth 0 °
Vertical Depth 7850.00 ft
Calibration Method: Linear Multiplier
Borehole Collapse Type: Bottomhole Borehole Collapse Pressure
Calculated Calibration Ratio 0.79
Unit No. 3 - Stratigraphic Unit: Rotliegendes 2, Lithology: sandstone (Top MD = 14400.00 ft ; Top TVD = 7858.14 ft)
Borehole Collapse Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
Residual Cohesion 300 psi
Residual Friction Angle 40 °
CASE PROFILE
Analysis Type Hole Collapse While Drilling
Analysis Mode Calculate required bottomhole pressures or mud densities to achieve tolerable borehole
instability risks.
Date and Time 26-Nov-2008, 14:35
13000.03 7631.42 3951 0.52 9.96 13595.83 7727.91 4004 0.52 9.97
13045.83 7638.84 3955 0.52 9.96 13641.67 7735.33 4008 0.52 9.97
13091.67 7646.26 3959 0.52 9.96 13687.50 7742.75 4012 0.52 9.97
13137.50 7653.68 3963 0.52 9.96 13733.33 7750.18 4016 0.52 9.97
13183.33 7661.11 3967 0.52 9.96 13779.17 7757.60 4021 0.52 9.97
13229.17 7668.53 3971 0.52 9.96 13825.00 7765.02 4025 0.52 9.97
13275.00 7675.95 3975 0.52 9.97 13870.83 7772.44 4029 0.52 9.97
13320.83 7683.37 3979 0.52 9.97 13916.67 7779.87 4033 0.52 9.98
13366.67 7690.80 3984 0.52 9.97 13962.50 7787.29 4037 0.52 9.98
13412.50 7698.22 3988 0.52 9.97 14008.33 7794.71 4041 0.52 9.98
13458.33 7705.64 3992 0.52 9.97 14054.17 7802.14 4045 0.52 9.98
13504.17 7713.06 3996 0.52 9.97 14099.97 7809.55 4049 0.52 9.98
13550.00 7720.49 4000 0.52 9.97
Unit No. 2 - Stratigraphic Unit: Lower Cretaceous, Lithology: shale w/ siltstone, Model Type: Default (3D Linear Elastic Model with Mohr-
Coulomb Failure Criterion)
14100.03 7809.56 4016 0.51 9.89 14262.50 7835.87 4030 0.51 9.90
14112.50 7811.58 4017 0.51 9.89 14275.00 7837.90 4031 0.51 9.90
14125.00 7813.61 4018 0.51 9.89 14287.50 7839.92 4032 0.51 9.90
14137.50 7815.63 4019 0.51 9.90 14300.00 7841.95 4034 0.51 9.90
14150.00 7817.66 4020 0.51 9.90 14312.50 7843.97 4035 0.51 9.90
Depth Depth BHP Pressure EMW Depth Depth BHP Pressure EMW
ft KB MD ft KB TVD psi Gradient ppg ft KB MD ft KB TVD psi Gradient ppg
psi/ft psi/ft
14162.50 7819.68 4021 0.51 9.90 14325.00 7846.00 4036 0.51 9.90
14175.00 7821.70 4022 0.51 9.90 14337.50 7848.02 4037 0.51 9.90
14187.50 7823.73 4023 0.51 9.90 14350.00 7850.04 4038 0.51 9.90
14200.00 7825.75 4025 0.51 9.90 14362.50 7852.07 4039 0.51 9.90
14212.50 7827.78 4026 0.51 9.90 14375.00 7854.09 4040 0.51 9.90
14225.00 7829.80 4027 0.51 9.90 14387.50 7856.12 4041 0.51 9.90
14237.50 7831.83 4028 0.51 9.90 14399.97 7858.14 4042 0.51 9.90
14250.00 7833.85 4029 0.51 9.90
Unit No. 3 - Stratigraphic Unit: Rotliegendes 2, Lithology: sandstone, Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
14400.03 7858.15 4202 0.53 10.29 15142.63 7978.41 4286 0.54 10.34
14457.12 7867.39 4208 0.53 10.29 15199.75 7987.66 4292 0.54 10.34
14514.25 7876.65 4213 0.53 10.29 15256.88 7996.91 4297 0.54 10.34
14571.38 7885.90 4219 0.54 10.30 15314.00 8006.16 4303 0.54 10.34
14628.50 7895.15 4225 0.54 10.30 15371.13 8015.41 4308 0.54 10.34
14685.63 7904.40 4241 0.54 10.32 15428.25 8024.66 4314 0.54 10.35
14742.75 7913.65 4246 0.54 10.33 15485.38 8033.92 4320 0.54 10.35
14799.88 7922.90 4252 0.54 10.33 15542.50 8043.17 4325 0.54 10.35
14857.00 7932.15 4258 0.54 10.33 15599.63 8052.42 4331 0.54 10.35
14914.13 7941.40 4263 0.54 10.33 15656.75 8061.67 4337 0.54 10.35
14971.25 7950.66 4269 0.54 10.33 15713.87 8070.92 4342 0.54 10.35
15028.37 7959.91 4275 0.54 10.33 15770.97 8080.17 4348 0.54 10.35
15085.50 7969.16 4280 0.54 10.34
1000 S
North/South Axis (ft)
2000 S
3000 S
4000 S
5000 S
6000 S
7000 S
8000 S
9000 S
10000 S
2000 W 0 2000 E 4000 E 6000 E 8000 E
1000 W 1000 E 3000 E 5000 E 7000 E 9000 E
Well Inclination = 81 °
Well Azimuth = 167 °
Formation: Lower Cretaceous
Lithology: shale w/ siltstone
10.5
Model: 3D Linear Elastic
Failure Criteria: Mohr Coulomb
Allowed Yield
Circumference = 60.0 %
Äçä = 0.40
9.5 cñ = 300 psi
Åñ = 35 °
E = 0.15 E+06 psi
Í = 0.35
À = 1.00
Plane of Weakness
Dip = 20 °
9.0 Dip Direction = 340 °
0 10 20 30 40 50 60 70 80 90cç = 75 psi
Åç = 30 °
N
Analysis Mode: Pressure 10.02
Analysis Type: Drilling Collapse
Bedding Friction 25 30 35 °
Angle
Peak Cohesion
225 300 375 psi
CASE PROFILE
Analysis Type Hole Collapse While Drilling
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Date and Time 26-Nov-2008, 14:37
BOREHOLE PROPERTIES
Bottom Borehole
No. Depth Diameter
ft KB MD in
1 650.00 23.00
2 2500.00 16.00
3 7500.00 12.25
4 12900.00 8.50
5 15771.43 6.00
TUBULAR PROPERTIES
No. Tubular Type Top Depth Btm Depth ID OD
ft KB MD ft KB MD in in
1 casing 0.00 650.00 19.00 20.00
2 casing 0.00 2500.00 12.42 13.38
3 casing 0.00 7500.00 8.68 9.63
4 liner 7470.00 12900.00 6.18 7.00
STRATIGRAPHIC UNITS
No. Stratigraphic Unit Lithology Pattern Top Depth Btm Depth Top Depth Btm Depth
ft KB MD ft KB MD ft KB TVD ft KB TVD
1 Reservoir 1 sandstone 13000.00 14100.00 7631.42 7809.56
2 Cretaceous shale w/ siltstone 14100.00 14400.00 7809.56 7858.14
3 Reservoir 2 sandstone 14400.00 15771.00 7858.14 8080.17
INITIAL FORMATION STRESS AND PRESSURE GRADIENT DATA
No. Stratigraphic Unit Top Top SHmin Sv Sv SHmax SHmax SHmin SHmin Pr Pr
Depth Depth Azi. Grad Grad Grad Grad Grad Grad Grad Grad
ft KB MD ft KB TVD (Top) (Btm) (Top) (Btm) (Top) (Btm) (Top) (Btm)
psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft
1 Reservoir 1 13000.00 7631.42 77 0.90 0.90 0.80 0.80 0.72 0.72 0.52 0.52
2 Cretaceous 14100.00 7809.56 77 0.90 0.90 0.80 0.80 0.72 0.72 0.52 0.52
3 Reservoir 2 14400.00 7858.14 77 0.90 0.90 0.80 0.80 0.72 0.72 0.52 0.52
WELLBORE PRESSURE OR PRODUCTION RATE DATA
Top Top Underbalance Underbalance
No. Stratigraphic Unit Lithology Depth Depth Pressure (Top) Pressure (Btm)
ft KB MD ft KB TVD psi psi
1 Reservoir 1 sandstone 13000.00 7631.42 500 400
2 Cretaceous shale w/ siltstone 14100.00 7809.56 400 350
3 Reservoir 2 sandstone 14400.00 7858.14 350 200
BASIC ROCK MECHANICAL PROPERTIES
No. Stratigraphic Unit Lithology Top Top Peak Peak Poisson's Young's Biot's
Depth Depth Cohesion Friction Ratio Modulus Coefficient
ft KB MD ft KB psi Angle E+06 psi
TVD
1 Reservoir 1 sandstone 13000.00 7631.42 400 45 0.32 0.20 1.00
shale w/
2
Cretaceous siltstone 14100.00 7809.56 300 35 0.35 0.15 1.00
3 Reservoir 2 sandstone 14400.00 7858.14 800 45 0.30 0.40 1.00
FLUID PENETRATION PROPERTIES
No. Stratigraphic Unit Lithology Top Top Capillary Filter Cake
Depth Depth Threshold or Wall
ft KB MD ft KB TVD Pressure Coating
psi Efficiency
Unit No. 2 - Stratigraphic Unit: Cretaceous, Lithology: shale w/ siltstone (Top MD = 14100.00 ft ; Top TVD = 7809.56 ft)
Borehole Collapse Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
Residual Cohesion 150 psi
Residual Friction Angle 35 °
Unit No. 3 - Stratigraphic Unit: Reservoir 2, Lithology: sandstone (Top MD = 14400.00 ft ; Top TVD = 7858.14 ft)
Borehole Collapse Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
Residual Cohesion 300 psi
Residual Friction Angle 40 °
CASE PROFILE
Analysis Type Hole Collapse While Drilling
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Date and Time 26-Nov-2008, 14:37
Unit No. 1 - Stratigraphic Unit: Reservoir 1, Lithology: sandstone, Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
Unit No. 2 - Stratigraphic Unit: Cretaceous, Lithology: shale w/ siltstone, Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure
Criteria
Unit No. 3 - Stratigraphic Unit: Reservoir 2, Lithology: sandstone, Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
13300
13500
Measured Depth (ft KB)
13700
13900
14100
Cretaceous
14300
Reservoir 2
14500
14700
14900
15100
15300
15500
15700
1.0
7 8 9 10 11 12 13
Residual Friction 25 35 40 °
Angle
CASE PROFILE
Analysis Type Sand Production / Openhole Stability Under Drawdown
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Date and Time 26-Nov-2008, 14:38
BOREHOLE PROPERTIES
Bottom Borehole
No. Depth Diameter
m KB MD mm
1 1000.00 216.0
TUBULAR PROPERTIES
No. Tubular Type Top Depth Btm Depth ID OD
m KB MD m KB MD mm mm
1 casing 0.00 900.00 140.0 152.0
STRATIGRAPHIC UNITS
No. Stratigraphic Unit Lithology Pattern Top Depth Btm Depth Top Depth Btm Depth
m KB MD m KB MD m KB TVD m KB TVD
1 Interval 1 sandstone 900.00 1000.00 900.00 1000.00
CASE PROFILE
Analysis Type Sand Production / Openhole Stability Under Drawdown
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Date and Time 26-Nov-2008, 14:38
/ Original Borehole
/ Yielded Zone
(Drawn to Scale)
0.00
6 7 8 9 10 11
CASE PROFILE
Analysis Type Sand Production / Openhole Stability Under Drawdown
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Date and Time 26-Nov-2008, 14:39
BOREHOLE PROPERTIES
Bottom Borehole
No. Depth Diameter
m KB MD mm
1 800.00 220.0
TUBULAR PROPERTIES
No. Tubular Type Top Depth Btm Depth ID OD
m KB MD m KB MD mm mm
1 casing 0.00 700.00 140.0 152.0
STRATIGRAPHIC UNITS
No. Stratigraphic Unit Lithology Pattern Top Depth Btm Depth Top Depth Btm Depth
m KB MD m KB MD m KB TVD m KB TVD
1 Interval 1 sandstone 700.00 800.00 700.00 800.00
CASE PROFILE
Analysis Type Sand Production / Openhole Stability Under Drawdown
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Date and Time 26-Nov-2008, 14:39
Model: 2D Elastoplastic
Radius, R/rperf
1.1
0 5 10 15
Well Inclination = 0 °
Well Azimuth = 0 °
Formation: Interval 1
Lithology: sandstone
20
Model: 2D Elastoplastic
Failure Criteria: Mohr Coulomb
cñ = 4.0 MPa
10 Åñ = 35 °
có = 1.0 MPa
Åó = 30 °
E = 8.00 GPa
Í = 0.30
À = 1.00
5
0 5 10 15
Lithology: sandstone
1.4
Model: 2D Elastoplastic
Failure Criteria: Mohr Coulomb
Radius, R/rperf
1.0
4 5 6 7 8
CASE PROFILE
Analysis Type Hole Collapse While Drilling
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Date and Time 26-Nov-2008, 14:43
BOREHOLE PROPERTIES
Bottom Borehole
No. Depth Diameter
ft KB MD in
1 600.00 12.25
2 3280.80 8.67
TUBULAR PROPERTIES
* No tubular data were entered for this case.
STRATIGRAPHIC UNITS
No. Stratigraphic Unit Lithology Pattern Top Depth Btm Depth Top Depth Btm Depth
ft KB MD ft KB MD ft KB TVD ft KB TVD
1 Paleocene Mudstone mudstone 656.10 1640.40 656.10 1640.40
2 Cretaceous Shale shale 1640.40 2788.60 1640.40 2788.60
3 Jurassic Reservoir sandstone 2788.60 3280.80 2788.60 3280.80
Unit No. 2 - Stratigraphic Unit: Cretaceous Shale, Lithology: shale (Top MD = 1640.40 ft ; Top TVD = 1640.40 ft)
Borehole Collapse Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
Residual Cohesion 220 psi
Residual Friction Angle 20 °
Pore Pressure Effects:
Formation Fluid Viscosity 0.5 cp
Fluid Viscosity in Yielded Zone 0.5 cp
Elastic Rock Permeability 0.001 md
Yielded Rock Permeability 0.005 md
Maximum Swab Pressure 20 psi
Unit No. 3 - Stratigraphic Unit: Jurassic Reservoir, Lithology: sandstone (Top MD = 2788.60 ft ; Top TVD = 2788.60 ft)
Borehole Collapse Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
Residual Cohesion 220 psi
Residual Friction Angle 28 °
Fracture Roughness Angle 45 °
Rubble Bulking Factor 2.00
Pore Pressure Effects:
Formation Fluid Viscosity 0.5 cp
Fluid Viscosity in Yielded Zone 0.5 cp
Elastic Rock Permeability 100 md
Yielded Rock Permeability 200 md
CASE PROFILE
Analysis Type Hole Collapse While Drilling
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Date and Time 26-Nov-2008, 14:43
Unit No. 1 - Stratigraphic Unit: Paleocene Mudstone, Lithology: mudstone, Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure
Criteria
Unit No. 2 - Stratigraphic Unit: Cretaceous Shale, Lithology: shale, Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
Unit No. 3 - Stratigraphic Unit: Jurassic Reservoir, Lithology: sandstone, Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure
Criteria
900
1100
Measured Depth (ft KB)
1300
1500
1900
2100
2300
2500
2700
Jurassic Reservoir
2900
3100
Residual Friction 13 15 20 °
Angle Base Case R/rw = 1.44
Lithology: sandstone
1.3 Model: 2D Elastoplastic
Failure Criteria: Mohr Coulomb
Radius, R/rw
Äçä = 0.80
1.1 cñ = 580 psi
Åñ = 35 °
có = 220 psi
Åó = 28 °
E = 2.00 E+06 psi
Í = 0.30
À = 1.00
1.0
7 8 9 10 11
CASE PROFILE
Analysis Type Sand Production / Openhole Stability Under Drawdown
Analysis Mode Calculate borehole instability risks for specified bottomhole pressures or mud densities.
Borehole Collapse Model Type Default (3D Linear Elastic Model with Mohr-Coulomb Failure Criterion)
Date and Time 26-Nov-2008, 14:45
1.10
1 3 5 7 9 11
Drawdown, MPa
cñ = 3.0 MPa
Åñ = 35 °
1.23 E = 1.00 GPa
Í = 0.30
À = 1.00
1.22
1.21
2.0 2.5 3.0 3.5 4.0
Drawdown, MPa
CASE PROFILE
Analysis Type Sand Production / Openhole Stability Under Drawdown
Analysis Mode Calculate required bottomhole pressures or mud densities to achieve tolerable borehole
instability risks.
Borehole Collapse Model Type 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
Date and Time 26-Nov-2008, 14:46
Filename: Case 16 -Horizontal Well with Effects of Depletion 2D.bhs 26-Nov-2008 14:46
Case Title: Extended Reach Horizontal Well 2D with Depletion 2 of 3 ABC Oil & Gas
Filename: Case 16 -Horizontal Well with Effects of Depletion 2D.bhs 26-Nov-2008 14:46
Case Title: Extended Reach Horizontal Well 2D with Depletion 3 of 3 ABC Oil & Gas
Model: 2D Elastoplastic
25 Failure Criteria: Mohr Coulomb
20
Òüîêï Azi = 80 °
cñ = 9.0 MPa
10 Åñ = 40 °
có = 3.0 MPa
Åó = 35 °
E = 3.03 GPa
5 Í = 0.25
À = 0.90
0
0 10 20 30 40
Model: 2D Elastoplastic
Failure Criteria: Mohr Coulomb
Radius, R/rw
cñ = 9.0 MPa
Åñ = 40 °
có = 3.0 MPa
1.5 Åó = 35 °
E = 3.03 GPa
Í = 0.25
À = 0.90
1.4
5 7 9 11 13 15
Drawdown, MPa
Filename: Case 16 -Horizontal Well with Effects of Depletion 2D.bhs 26-Nov-2008 14:46
Case Title: Complex Multi-zone Offshore Well 1 of 9 ABC Oil Ltd.
CASE PROFILE
Analysis Type Hole Collapse and Fracture Breakdown / Lost Circulation
Analysis Mode Calculate required bottomhole pressures or mud densities to achieve tolerable borehole
instability risks.
Date and Time 26-Nov-2008, 14:47
Depth Incl. Azimuth Depth N/S E/W DLS Depth Incl. Azimuth Depth N/S E/W DLS
ft KB MD ° ° ft KB TVD Distance Distance deg/100 ft KB MD ° ° ft KB TVD Distance Distance deg/100
ft ft ft ft ft ft
16955.70 61.39 249.08 16512.04 -426.51 -1115.75 2.7 20300.00 87.96 219.60 17166.80 -2499.53 -3505.59 0.0
17000.00 61.39 249.08 16533.26 -440.39 -1152.07 0.0 20400.00 87.96 219.60 17170.36 -2576.53 -3569.30 0.0
17100.00 61.39 249.08 16581.14 -471.74 -1234.08 0.0 20503.17 87.96 219.60 17174.03 -2655.97 -3635.02 0.0
17200.00 61.39 249.08 16629.03 -503.09 -1316.08 0.0 20600.00 87.96 219.60 17177.48 -2730.54 -3696.70 0.0
17300.00 61.39 249.08 16676.91 -534.43 -1398.08 0.0 20700.00 87.96 219.60 17181.04 -2807.54 -3760.40 0.0
17355.34 61.39 249.08 16703.41 -551.78 -1443.46 0.0 20800.00 87.96 219.60 17184.60 -2884.54 -3824.10 0.0
17400.00 62.08 248.08 16724.56 -566.15 -1480.08 2.7 20900.00 87.96 219.60 17188.15 -2961.54 -3887.81 0.0
17500.00 63.66 245.90 16770.16 -600.94 -1561.98 2.7 21003.17 87.96 219.60 17191.83 -3040.99 -3953.53 0.0
17600.00 65.27 243.78 16813.27 -639.31 -1643.64 2.7
BOREHOLE PROPERTIES
Bottom Borehole
No. Depth Diameter
ft KB MD in
1 9633.00 36.00
2 14250.00 16.00
3 17355.00 12.25
4 21003.17 8.50
TUBULAR PROPERTIES
No. Tubular Type Top Depth Btm Depth ID OD
ft KB MD ft KB MD in in
1 surface casing 0.00 9833.00 35.25 36.00
2 casing 0.00 14250.00 15.50 16.00
3 casing 0.00 17355.00 9.00 9.65
4 casing 0.00 18902.00 8.00 8.50
STRATIGRAPHIC UNITS
No. Stratigraphic Unit Lithology Pattern Top Depth Btm Depth Top Depth Btm Depth
ft KB MD ft KB MD ft KB TVD ft KB TVD
1 Mud mudstone 10950.00 12650.00 10950.00 12650.00
2 Sand sandstone 12650.00 14121.00 12650.00 14121.00
3 Mudstone mudstone 14121.00 18142.12 14121.00 17000.00
4 Sand sandstone 18142.12 18600.00 17000.00 17091.97
5 Sand sandstone 18600.00 19125.00 17091.97 17124.97
6 Shale shale 19125.79 21003.17 17125.00 17191.83
INITIAL FORMATION STRESS AND PRESSURE GRADIENT DATA
No. Stratigraphic Unit Top Top SHmin Sv Sv SHmax SHmax SHmin SHmin Pr Pr
Depth Depth Azi. Grad Grad Grad Grad Grad Grad Grad Grad
ft KB MD ft KB TVD (Top) (Btm) (Top) (Btm) (Top) (Btm) (Top) (Btm)
psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft psi/ft
1 Mud 10950.00 10950.00 110 0.49 0.54 0.48 0.54 0.48 0.53 0.44 0.48
2 Sand 12650.00 12650.00 110 0.54 0.58 0.54 0.57 0.53 0.55 0.48 0.49
3 Mudstone 14121.00 14121.00 110 0.58 0.65 0.57 0.63 0.55 0.61 0.49 0.55
4 Sand 18142.12 17000.00 110 0.65 0.65 0.63 0.63 0.61 0.61 0.55 0.55
5 Sand 18600.00 17091.97 110 0.65 0.65 0.63 0.63 0.61 0.60 0.55 0.53
6 Shale 19125.79 17125.00 110 0.65 0.65 0.63 0.62 0.60 0.59 0.53 0.50
BASIC ROCK MECHANICAL PROPERTIES
No. Stratigraphic Unit Lithology Top Top Peak Peak Poisson's Young's Tensile Biot's
Depth Depth Cohesion Friction Ratio Modulus Strength Coefficient
ft KB MD ft KB psi Angle E+06 psi psi
TVD
1 Mud mudstone 10950.00 10950.00 262 21 0.43 0.05 0 1.00
2 Sand sandstone 12650.00 12650.00 525 19 0.42 0.16 0 1.00
3 Mudstone mudstone 14121.00 14121.00 550 22 0.40 0.20 0 1.00
4 Sand sandstone 18142.12 17000.00 700 22 0.38 0.27 0 1.00
5 Sand sandstone 18600.00 17091.97 700 22 0.38 0.27 0 1.00
6 Shale shale 19125.79 17125.00 1030 19 0.36 0.27 0 1.00
FLUID PENETRATION PROPERTIES
No. Stratigraphic Unit Lithology Top Top Capillary Filter Cake Fluid
Depth Depth Threshold or Wall Penetration
ft KB MD ft KB TVD Pressure Coating Coefficient
psi Efficiency
No. Stratigraphic Unit Lithology Top Top Capillary Filter Cake Fluid
Depth Depth Threshold or Wall Penetration
ft KB MD ft KB TVD Pressure Coating Coefficient
psi Efficiency
Unit No. 2 - Stratigraphic Unit: Sand, Lithology: sandstone (Top MD = 12650.00 ft ; Top TVD = 12650.00 ft)
Fracture Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
Borehole Collapse Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
Residual Cohesion 150 psi
Residual Friction Angle 16 °
Maximum Tolerable Wellbore Instability Failure Criteria: Average Yielded Zone Radius (R/rw)
Tolerable Yield Zone, R/rw 1.50
Unit No. 3 - Stratigraphic Unit: Mudstone, Lithology: mudstone (Top MD = 14121.00 ft ; Top TVD = 14121.00 ft)
Fracture Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
Borehole Collapse Model Type: 3D Linear Elastic Model with Modified Lade Failure Criterion
Modified Lade Cohesive Strength Parameter 1361 psi
Modified Lade Frictional Strength Parameter 6.66
Maximum Tolerable Wellbore Instability Failure Criteria: Percentage of the Hole Circumference
Tolerable Yield Circumf., % 50.00
Unit No. 4 - Stratigraphic Unit: Sand, Lithology: sandstone (Top MD = 18142.12 ft ; Top TVD = 17000.00 ft)
Fracture Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
Borehole Collapse Model Type: Default (3D Linear Elastic Model with Mohr-Coulomb Failure Criterion)
Maximum Tolerable Wellbore Instability Failure Criteria: Percentage of the Hole Circumference
Tolerable Yield Circumf., % 50.00
Unit No. 5 - Stratigraphic Unit: Sand, Lithology: sandstone (Top MD = 18600.00 ft ; Top TVD = 17091.97 ft)
Fracture Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
Borehole Collapse Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
Residual Cohesion 100 psi
Residual Friction Angle 22 °
Unit No. 6 - Stratigraphic Unit: Shale, Lithology: shale (Top MD = 19125.79 ft ; Top TVD = 17125.00 ft)
Fracture Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
Borehole Collapse Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
Residual Cohesion 100 psi
Residual Friction Angle 22 °
Maximum Tolerable Wellbore Instability Failure Criteria: Average Yielded Zone Radius (R/rw)
Tolerable Yield Zone, R/rw 1.50
CASE PROFILE
Analysis Type Hole Collapse and Fracture Breakdown / Lost Circulation
Analysis Mode Calculate required bottomhole pressures or mud densities to achieve tolerable borehole
instability risks.
Date and Time 26-Nov-2008, 14:47
10950.03 10950.03 5391 0.49 9.47 11870.83 11870.83 6144 0.52 9.96
11020.83 11020.83 5448 0.49 9.51 11941.67 11941.67 6204 0.52 10.00
11091.67 11091.67 5504 0.50 9.55 12012.50 12012.50 6264 0.52 10.03
11162.50 11162.50 5561 0.50 9.59 12083.33 12083.33 6324 0.52 10.07
11233.33 11233.33 5618 0.50 9.62 12154.17 12154.17 6385 0.53 10.11
11304.17 11304.17 5675 0.50 9.66 12225.00 12225.00 6446 0.53 10.15
11375.00 11375.00 5733 0.50 9.70 12295.83 12295.83 6507 0.53 10.18
11445.83 11445.83 5791 0.51 9.74 12366.67 12366.67 6569 0.53 10.22
11516.67 11516.67 5849 0.51 9.77 12437.50 12437.50 6630 0.53 10.26
11587.50 11587.50 5907 0.51 9.81 12508.33 12508.33 6692 0.54 10.30
11658.33 11658.33 5966 0.51 9.85 12579.17 12579.17 6755 0.54 10.33
11729.17 11729.17 6025 0.51 9.88 12649.97 12649.97 6817 0.54 10.37
11800.00 11800.00 6084 0.52 9.92
Unit No. 2 - Stratigraphic Unit: Sand, Lithology: sandstone, Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
12650.03 12650.03 7134 0.56 10.85 13446.79 13446.79 7724 0.57 11.05
12711.29 12711.29 7179 0.56 10.87 13508.08 13508.08 7770 0.58 11.07
12772.58 12772.58 7223 0.57 10.88 13569.38 13569.38 7816 0.58 11.08
12833.88 12833.88 7268 0.57 10.90 13630.67 13630.67 7862 0.58 11.10
12895.17 12895.17 7314 0.57 10.91 13691.96 13691.96 7909 0.58 11.12
12956.46 12956.46 7359 0.57 10.93 13753.25 13753.25 7955 0.58 11.13
13017.75 13017.75 7404 0.57 10.94 13814.54 13814.54 8002 0.58 11.15
13079.04 13079.04 7449 0.57 10.96 13875.83 13875.83 8048 0.58 11.16
13140.33 13140.33 7495 0.57 10.98 13937.13 13937.13 8095 0.58 11.18
Depth Depth BHP Pressure ECD Depth Depth BHP Pressure ECD
ft KB MD ft KB TVD psi Gradient ppg ft KB MD ft KB TVD psi Gradient ppg
psi/ft psi/ft
13201.63 13201.63 7540 0.57 10.99 13998.42 13998.42 8142 0.58 11.19
13262.92 13262.92 7586 0.57 11.01 14059.71 14059.71 8189 0.58 11.21
13324.21 13324.21 7632 0.57 11.02 14120.97 14120.97 8236 0.58 11.22
13385.50 13385.50 7678 0.57 11.04
Unit No. 3 - Stratigraphic Unit: Mudstone, Lithology: mudstone, Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
14121.03 14121.03 8165 0.58 11.13 16299.11 16119.92 10081 0.63 12.03
14288.55 14288.55 8312 0.58 11.19 16466.65 16233.66 10196 0.63 12.09
14456.09 14456.09 8460 0.59 11.26 16634.20 16338.45 10300 0.63 12.13
14623.64 14623.55 8610 0.59 11.33 16801.75 16433.82 10391 0.63 12.17
14791.19 14790.38 8762 0.59 11.40 16969.29 16518.55 10469 0.63 12.20
14958.73 14955.57 8915 0.60 11.47 17136.84 16598.78 10545 0.64 12.22
15126.28 15118.40 9070 0.60 11.54 17304.39 16679.01 10620 0.64 12.25
15293.83 15278.00 9224 0.60 11.62 17471.93 16757.36 10656 0.64 12.24
15461.37 15433.21 9377 0.61 11.69 17639.48 16829.28 10672 0.63 12.20
15628.92 15583.60 9528 0.61 11.77 17807.03 16894.17 10684 0.63 12.17
15796.47 15728.38 9676 0.62 11.84 17974.57 16951.02 10694 0.63 12.14
15964.01 15866.16 9819 0.62 11.91 18142.09 16999.99 10701 0.63 12.11
16131.56 15996.89 9954 0.62 11.97
Unit No. 4 - Stratigraphic Unit: Sand, Lithology: sandstone, Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
18142.15 17000.01 10813 0.64 12.24 18390.14 17057.73 10795 0.63 12.18
18161.20 17005.08 10812 0.64 12.23 18409.22 17061.40 10794 0.63 12.17
18180.28 17010.16 10810 0.64 12.23 18428.29 17064.76 10792 0.63 12.17
18199.35 17015.24 10809 0.64 12.22 18447.37 17068.13 10791 0.63 12.17
18218.43 17019.78 10808 0.64 12.22 18466.45 17071.49 10790 0.63 12.16
18237.51 17024.29 10807 0.63 12.21 18485.53 17074.86 10789 0.63 12.16
18256.59 17028.81 10805 0.63 12.21 18504.61 17078.08 10787 0.63 12.15
18275.67 17033.33 10804 0.63 12.21 18523.69 17080.86 10786 0.63 12.15
18294.75 17037.85 10803 0.63 12.20 18542.77 17083.64 10784 0.63 12.15
18313.82 17041.95 10801 0.63 12.20 18561.84 17086.42 10782 0.63 12.14
18332.90 17045.89 10799 0.63 12.19 18580.92 17089.20 10781 0.63 12.14
18351.98 17049.84 10798 0.63 12.19 18599.97 17091.97 10780 0.63 12.14
18371.06 17053.79 10796 0.63 12.18
Unit No. 5 - Stratigraphic Unit: Sand, Lithology: sandstone, Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
18600.03 17091.98 10780 0.63 12.14 18884.38 17116.13 10632 0.62 11.95
18621.88 17094.48 10766 0.63 12.12 18906.25 17117.19 10625 0.62 11.94
18643.75 17096.99 10751 0.63 12.10 18928.13 17117.96 10621 0.62 11.94
18665.63 17099.49 10737 0.63 12.08 18950.00 17118.74 10618 0.62 11.94
18687.50 17102.00 10723 0.63 12.07 18971.88 17119.52 10614 0.62 11.93
18709.38 17104.22 10710 0.63 12.05 18993.75 17120.30 10611 0.62 11.93
18731.25 17106.04 10699 0.63 12.04 19015.63 17121.08 10607 0.62 11.92
18753.13 17107.86 10688 0.62 12.02 19037.50 17121.86 10603 0.62 11.92
18775.00 17109.68 10676 0.62 12.01 19059.38 17122.63 10600 0.62 11.91
18796.88 17111.51 10665 0.62 11.99 19081.25 17123.41 10596 0.62 11.91
18818.75 17112.73 10657 0.62 11.98 19103.13 17124.19 10592 0.62 11.90
18840.63 17113.87 10649 0.62 11.97 19124.97 17124.97 10589 0.62 11.90
18862.50 17115.00 10640 0.62 11.96
Unit No. 6 - Stratigraphic Unit: Shale, Lithology: shale, Model Type: Default (3D Linear Elastic Model for Tensile Fracture Breakdown)
19125.82 17125.00 10189 0.59 11.45 20142.71 17161.20 10139 0.59 11.37
19204.01 17127.78 10186 0.59 11.44 20220.93 17163.98 10135 0.59 11.36
19282.24 17130.57 10182 0.59 11.44 20299.15 17166.77 10131 0.59 11.36
19360.46 17133.35 10178 0.59 11.43 20377.38 17169.55 10127 0.59 11.35
19438.69 17136.14 10174 0.59 11.42 20455.60 17172.34 10123 0.59 11.34
19516.91 17138.92 10170 0.59 11.42 20533.82 17175.12 10119 0.59 11.34
19595.13 17141.71 10166 0.59 11.41 20612.05 17177.91 10115 0.59 11.33
19673.36 17144.49 10162 0.59 11.41 20690.27 17180.69 10111 0.59 11.32
19751.58 17147.28 10158 0.59 11.40 20768.50 17183.47 10107 0.59 11.32
19829.81 17150.06 10155 0.59 11.39 20846.72 17186.26 10103 0.59 11.31
19908.03 17152.84 10151 0.59 11.39 20924.95 17189.04 10099 0.59 11.31
19986.26 17155.63 10147 0.59 11.38 21003.14 17191.83 10095 0.59 11.30
20064.48 17158.41 10143 0.59 11.37
Depth Depth BHP Pressure ECD Depth Depth BHP Pressure ECD
ft KB MD ft KB TVD psi Gradient ppg ft KB MD ft KB TVD psi Gradient ppg
psi/ft psi/ft
10950.03 10950.03 4681 0.43 8.23 11870.83 11870.83 5483 0.46 8.89
11020.83 11020.83 4740 0.43 8.28 11941.67 11941.67 5546 0.46 8.94
11091.67 11091.67 4799 0.43 8.33 12012.50 12012.50 5597 0.47 8.97
11162.50 11162.50 4858 0.44 8.37 12083.33 12083.33 5661 0.47 9.01
11233.33 11233.33 4918 0.44 8.42 12154.17 12154.17 5725 0.47 9.06
11304.17 11304.17 4978 0.44 8.47 12225.00 12225.00 5789 0.47 9.11
11375.00 11375.00 5038 0.44 8.52 12295.83 12295.83 5854 0.48 9.16
11445.83 11445.83 5099 0.45 8.57 12366.67 12366.67 5919 0.48 9.21
11516.67 11516.67 5172 0.45 8.64 12437.50 12437.50 5971 0.48 9.24
11587.50 11587.50 5233 0.45 8.69 12508.33 12508.33 6037 0.48 9.29
11658.33 11658.33 5295 0.45 8.74 12579.17 12579.17 6103 0.49 9.34
11729.17 11729.17 5357 0.46 8.79 12649.97 12649.97 6169 0.49 9.38
11800.00 11800.00 5420 0.46 8.84
Unit No. 2 - Stratigraphic Unit: Sand, Lithology: sandstone, Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
12650.03 12650.03 5740 0.45 8.73 13446.79 13446.79 6347 0.47 9.08
12711.29 12711.29 5783 0.45 8.76 13508.08 13508.08 6407 0.47 9.13
12772.58 12772.58 5827 0.46 8.78 13569.38 13569.38 6453 0.48 9.15
12833.88 12833.88 5871 0.46 8.80 13630.67 13630.67 6499 0.48 9.17
12895.17 12895.17 5929 0.46 8.85 13691.96 13691.96 6545 0.48 9.20
12956.46 12956.46 5973 0.46 8.87 13753.25 13753.25 6591 0.48 9.22
13017.75 13017.75 6018 0.46 8.90 13814.54 13814.54 6638 0.48 9.25
13079.04 13079.04 6062 0.46 8.92 13875.83 13875.83 6700 0.48 9.29
13140.33 13140.33 6107 0.46 8.94 13937.13 13937.13 6746 0.48 9.31
13201.63 13201.63 6166 0.47 8.99 13998.42 13998.42 6793 0.49 9.34
13262.92 13262.92 6211 0.47 9.01 14059.71 14059.71 6840 0.49 9.36
13324.21 13324.21 6256 0.47 9.04 14120.97 14120.97 6887 0.49 9.39
13385.50 13385.50 6302 0.47 9.06
Unit No. 3 - Stratigraphic Unit: Mudstone, Lithology: mudstone, Model Type: 3D Linear Elastic Model with Modified Lade Failure Criterion
14121.03 14121.03 6814 0.48 9.29 16299.11 16119.92 8761 0.54 10.46
14288.55 14288.55 6962 0.49 9.38 16466.65 16233.66 8886 0.55 10.53
14456.09 14456.09 7111 0.49 9.47 16634.20 16338.45 9002 0.55 10.60
14623.64 14623.55 7262 0.50 9.56 16801.75 16433.82 9107 0.55 10.66
14791.19 14790.38 7418 0.50 9.65 16969.29 16518.55 9201 0.56 10.72
14958.73 14955.57 7576 0.51 9.75 17136.84 16598.78 9281 0.56 10.76
15126.28 15118.40 7731 0.51 9.84 17304.39 16679.01 9361 0.56 10.80
15293.83 15278.00 7885 0.52 9.93 17471.93 16757.36 9434 0.56 10.83
15461.37 15433.21 8038 0.52 10.02 17639.48 16829.28 9501 0.56 10.86
15628.92 15583.60 8190 0.53 10.11 17807.03 16894.17 9558 0.57 10.89
15796.47 15728.38 8340 0.53 10.20 17974.57 16951.02 9606 0.57 10.90
15964.01 15866.16 8486 0.53 10.29 18142.09 16999.99 9650 0.57 10.92
16131.56 15996.89 8626 0.54 10.38
Unit No. 4 - Stratigraphic Unit: Sand, Lithology: sandstone, Model Type: Default (3D Linear Elastic Model with Mohr-Coulomb Failure Criterion)
18142.15 17000.01 9601 0.56 10.87 18390.14 17057.73 9629 0.56 10.86
18161.20 17005.08 9602 0.56 10.87 18409.22 17061.40 9630 0.56 10.86
18180.28 17010.16 9604 0.56 10.86 18428.29 17064.76 9633 0.56 10.86
18199.35 17015.24 9607 0.56 10.87 18447.37 17068.13 9636 0.56 10.86
18218.43 17019.78 9611 0.56 10.87 18466.45 17071.49 9634 0.56 10.86
18237.51 17024.29 9615 0.56 10.87 18485.53 17074.86 9637 0.56 10.86
18256.59 17028.81 9613 0.56 10.86 18504.61 17078.08 9638 0.56 10.86
18275.67 17033.33 9619 0.56 10.87 18523.69 17080.86 9641 0.56 10.86
18294.75 17037.85 9620 0.56 10.87 18542.77 17083.64 9638 0.56 10.86
18313.82 17041.95 9621 0.56 10.86 18561.84 17086.42 9641 0.56 10.86
18332.90 17045.89 9625 0.56 10.87 18580.92 17089.20 9644 0.56 10.86
18351.98 17049.84 9628 0.56 10.87 18599.97 17091.97 9644 0.56 10.86
18371.06 17053.79 9629 0.56 10.86
Depth Depth BHP Pressure EMW Depth Depth BHP Pressure EMW
ft KB MD ft KB TVD psi Gradient ppg ft KB MD ft KB TVD psi Gradient ppg
psi/ft psi/ft
Unit No. 5 - Stratigraphic Unit: Sand, Lithology: sandstone, Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
18600.03 17091.98 9608 0.56 10.82 18884.38 17116.13 9509 0.56 10.69
18621.88 17094.48 9604 0.56 10.81 18906.25 17117.19 9507 0.56 10.69
18643.75 17096.99 9601 0.56 10.81 18928.13 17117.96 9506 0.56 10.69
18665.63 17099.49 9575 0.56 10.78 18950.00 17118.74 9483 0.55 10.66
18687.50 17102.00 9572 0.56 10.77 18971.88 17119.52 9482 0.55 10.66
18709.38 17104.22 9569 0.56 10.77 18993.75 17120.30 9481 0.55 10.66
18731.25 17106.04 9544 0.56 10.74 19015.63 17121.08 9480 0.55 10.65
18753.13 17107.86 9542 0.56 10.73 19037.50 17121.86 9479 0.55 10.65
18775.00 17109.68 9539 0.56 10.73 19059.38 17122.63 9478 0.55 10.65
18796.88 17111.51 9537 0.56 10.72 19081.25 17123.41 9477 0.55 10.65
18818.75 17112.73 9513 0.56 10.70 19103.13 17124.19 9476 0.55 10.65
18840.63 17113.87 9512 0.56 10.69 19124.97 17124.97 9453 0.55 10.62
18862.50 17115.00 9510 0.56 10.69
Unit No. 6 - Stratigraphic Unit: Shale, Lithology: shale, Model Type: 2D Elastoplastic Model with Mohr-Coulomb Failure Criteria
19125.82 17125.00 9029 0.53 10.15 20142.71 17161.20 8959 0.52 10.05
19204.01 17127.78 9028 0.53 10.14 20220.93 17163.98 8959 0.52 10.04
19282.24 17130.57 9028 0.53 10.14 20299.15 17166.77 8936 0.52 10.02
19360.46 17133.35 9028 0.53 10.14 20377.38 17169.55 8936 0.52 10.02
19438.69 17136.14 9005 0.53 10.11 20455.60 17172.34 8936 0.52 10.01
19516.91 17138.92 9005 0.53 10.11 20533.82 17175.12 8914 0.52 9.99
19595.13 17141.71 9005 0.53 10.11 20612.05 17177.91 8913 0.52 9.98
19673.36 17144.49 9004 0.53 10.11 20690.27 17180.69 8913 0.52 9.98
19751.58 17147.28 8982 0.52 10.08 20768.50 17183.47 8913 0.52 9.98
19829.81 17150.06 8982 0.52 10.08 20846.72 17186.26 8891 0.52 9.95
19908.03 17152.84 8982 0.52 10.08 20924.95 17189.04 8890 0.52 9.95
19986.26 17155.63 8960 0.52 10.05 21003.14 17191.83 8890 0.52 9.95
20064.48 17158.41 8959 0.52 10.05
11000 Mud
12000
Measured Depth (ft KB)
Sand
13000
14000
Mudstone
15000
16000
17000
18000
Sand
Sand
19000
Shale
20000
21000
10000
10.25
11000
9.57
12000
TVD (ft)
13000 8.90
14000 8.23
ppg
15000
16000
17000
18000
-2000 -1000 0 1000 2000 3000 4000 5000 6000 7000
0 10.25
North/South Axis (ft)
500 S
9.57
1000 S
8.90
1500 S
2000 S 8.23
ppg
2500 S
3000 S
3500 S
4000 S
4000 W 3000 W 2000 W 1000 W 0
4500 W 3500 W 2500 W 1500 W 500 W 500 E
APPENDIX F
* To determine the complete STABView version and build number, please check under Help Æ About STABView 3.8, e.g., Version
3.8 (Build 201). This last build number is very important.
Please attach a sample STABView input and/or output file or printed output that produces the problem.
AG Use Only
INDEX
Discontinuity, 73, 75
Aadnoy, B., 37 Discontinuity slip, 112
Addis, T., 29 Drilling applications, 1
Advanced Geotechnology, vii Drilling instability, 11
Analysis mode, 77 Dusseault, M., 44
Analysis type, 77
Anisotropic elastic properties, 60 Effective axial stress, 33
Anisotropic stress regime, 48 Effective radial stress, 32
Apparent capillary cohesion, 24, 26 Effective stress, 17, 18
Average yielded zone radius, 48 Effective tangential stress, 33
Axial strain, 20 Elastic models, 13, 32, 36
Axial stress, 34 Elastic properties, 20
Elastic-brittle-plastic material, 43, 45
Barefoot openhole completion, 106, 107, 2, 3 Elasto-plastic model, 13, 43, 48
Bear, J., 27 Elasto-plastic modeling, 49
Bedding planes, 39 Elasto-plastic yielded zone, 47
Biot’s coefficient, 17, 33, 84 Emulsifier, 71
Borehole deformations and strains, 13 Equivalent circulating density, 110
Borehole radius, 45 Ewy, R., 22, 36
Borehole sizes, 82 Extended reach well, 107, 3
Borehole stresses and pressures, 13 External filter cake, 110
Bottomhole pressure, 33, 35
Bratli, R., 48 Failure angle, 48
Breakout angle, 42, 43 Fault(s), 12, 75
Bruno, M., 24 Filter cake, 47, 63
Butler, R., 67 Filter cake efficiency, 85
Buttons, 88 Fissile shale, 106, 2
Fluid penetration coefficient, 85
Calibration, 14, 42, 86 Fluid viscosity, 48
Capillarity, 12 Formation evaluation applications, 1
Capillary pressures, 24 Formation fluid penetration properties, 85
Capillary threshold pressure, 65 Formation pressures, 83
Cartesian coordinate system, 36 Fracture breakdown, 106, 111, 2
Case properties, 81 Fracture initiation, 55, 56
Chemical effects, 13 Fracture re-opening tendency, 75
Chemical osmosis, 71 Fracturing, 11
Cohesion, 22 Friction angle, 22
Collapse calibration, 106, 2 Fully-penetrating fluid, 56
Completions, 1
Compressive shear yielding, 21, 31, 34 Gas Compressibility, 66
Coordinate systems, 36 Gas compressibility factor, 67, 69
Copyright, 2 Gas formation volume factor, 70
Critical breakout angle, 42 Gas gravity, 68
Cross section, 103 Gas property correlations, 68
Cylindrical perforation, 107, 3 Gas Research Institute, 72
Gas viscosity, 70
Database(s), 14
Deep fluid penetration, 57 Hawkes, C., 29, 31, 44
Demonstration cases, 105 Help & Documentation, 15
Depletion induced stresses, 28, 106, 2 Hemi-spherical cavities, 31
Depth profile, 97 Hemi-spherical perforations, 48, 107, 3
Detournay, E., 47 Hoek and Brown yield criterion, 22, 23, 36, 106, 2
Directional wells, 111 Hoek, E., 36
STABView 3.8 User’s Manual I-4 November, 2008