Ugwu Thesis
Ugwu Thesis
Ugwu Thesis
A Thesis
by
MASTER OF SCIENCE
December 2008
A Thesis
by
MASTER OF SCIENCE
Approved by:
December 2008
ABSTRACT
Cement Fatigue and HPHT Well Integrity with Application to Life of Well Prediction.
(December 2008)
In order to keep up with the world’s energy demands, oil and gas producing companies
have taken the initiative to explore offshore reserves or drill deeper into previously
existing wells. The consequence of this, however, has to deal with the high temperatures
For an oil well to maintain its integrity and be produced effectively and
completion. This complete zonal isolation can be compromised due to factors that come
into play when oil well cement experiences cyclic loading conditions which can lead to
cement material depending on stress levels and number of cycles. There have been a lot
concrete structures but the fatigue behavior of oil well cement is still relatively unknown
to engineers.
iv
Research in the area of oil well cement design has led to improved cement
designs and cementing practices but yet many cement integrity problems persist and this
casing cement bond under HPHT well conditions that can lead to best practices and a
model to predict well life. An analytical model, which can be used to evaluate stresses in
the cement sheath based on actual wellbore parameters, was developed and combined
effectively with finite element models to evaluate the fatigue and static loading behavior
of a well.
casing, cement and formation as well loading conditions play a very big role in the static
Finally, recommendations for future work on this subject were also presented in
order to understand all tenets of cement fatigue and to develop governing equations.
v
DEDICATION
To God who gave me the strength and intellect to venture into this work, this work and
encouragement, even in the face of difficulty, to me, this work is especially dedicated to
you.
To everyone else who helped along the way, you are a great part of this work.
vi
ACKNOWLEDGMENTS
I would love to express my profound gratitude to my committee chairs, Dr. Xin-Lin Gao
and Dr. Catalin Teodoriu, for their untiring efforts in seeing to the success of this work.
Your willingness to devote your time to answer my questions whenever I called and your
Further appreciation is expressed to Dr. Palazzolo and Dr. Muliana, for making time
I would also not forget to acknowledge the minerals management services (MMS)
To all my friends here at Texas A&M, especially Jerome Rajnauth, your warmth
and care made my transition through this great school much easier and enjoyable.
vii
NOMENCLATURE
LS Low Strength
HS High Strength
SF Safety Factor
E Young’s Modulus
v Poisson Ratio
σa Stress Amplitude
σm Mean Stress
viii
TABLE OF CONTENTS
Page
ABSTRACT…………………………………………………………………………… iii
DEDICATION…………………………………………………………………………. v
ACKNOWLEDGMENTS……………………………………………………………… vi
NOMENCLATURE…………………………………………………………………… vii
LIST OF FIGURES…………………………………………………………………….. x
CHAPTER
1.1 Background……………………………………………........................... 2
1.1.1 HPHT Well Integrity……………………………………………... 2
1.1.2 Failure Criteria……………………………………………………. 3
1.1.3 Fatigue…………………………………………………………….. 4
1.1.4 Fatigue in Cement/Concrete………………………………………. 6
1.2 Research Objectives………………………………………....................... 11
1.3 Approach...……………………………………………………………… 11
II OILFIELD CEMENTING…………………………………………………………… 12
CHAPTER Page
5.1 Conclusions……………………………………………………..…....... 88
5.2 Recommendations for Future Work…………………………………… 90
REFERENCES…………………………………………………………………………. 91
VITA………………………………………………………………………..................... 94
x
LIST OF FIGURES
FIGURE Page
1.2 Comparison of Fatigue Strain for High and Low Strength Cements ………….. 7
1.3 Comparison of Max. Stress Levels to Number of Cycles for Different Cement
Strengths………………………………………………………..………………. 8
2.3 Tangential Stress for Hard and Soft Cement Systems……… ………………… 31
3.4 Failure Envelope for Triaxial Compression and Tensile Stress State …………. 38
3.8 Equivalent Stress for Casing –Cement –Formation Model with Meshing ……. 48
3.9 Equivalent Stress Comparison for Analytical and FEA Models ………………. 48
FIGURE Page
3.13 Zero Based Loading, Goodman’s Diagram and Fatigue Options ……………... 53
3.14 Fully Reversed Loading, Goodman’s Diagram and Fatigue Options ………….. 54
6
4.1 Equivalent Stress with Formation E=3 x 10 psi and v=0.42 for Scenario 1………. 57
6
4.2 Equivalent Stress with Formation E=1 x 10 psi and v=0.3 for Scenario 1………... 58
4.3 Equivalent Stress in Three Cement Systems with pi=15000 psi and pf=1000 psi 59
4.4 Tangential Stress in Three Cement Systems with pi=15000 psi and pf=1000 psi 59
4.5 Radial Stress in Three Cement Systems with pi=15000 psi and pf=1000 psi....... 60
6
4.6 Von Mises Stress with Formation E=1 x 10 psi and v=0.3 for Loading case 2…… 61
4.7 Equivalent Stress in Three Cement Systems with pi=15000 psi and pf=0 psi...... 62
4.8 Tangential Stress in Three Cement Systems with pi=15000 psi and pf=0 psi..... 62
4.9 Radial Stress in Three Cement Systems with pi=15000 psi and pf=0 psi………. 63
4.10 Equivalent Stress in Three Cement Systems with pi=4000 psi and pf=10000 psi 64
4.11 Tangential Stress in Three Cement Systems with pi=4000 psi and pf=10000 psi 64
4.12 Radial Stress in Three Cement Systems with pi=4000 psi and pf=10000 psi…... 65
4.13 Equivalent Stress with pi=15000 psi, pf=1000 psi and ΔT =150oF ……………. 66
4.14 Tangential Stress with pi=15000 psi, pf=1000 psi and ΔT =150oF ……………. 66
4.15 Radial Stress with pi=15000 psi, pf=1000 psi and ΔT =150oF ………................ 67
4.16 Equivalent Stress with pi=15000 psi, pf= 0 psi and ΔT =150oF ……………….. 68
4.17 Tangential Stress with pi=15000 psi, pf= 0 psi and ΔT =150oF ………………. 68
4.18 Radial Stress with pi=15000 psi, pf= 0 psi and ΔT =150oF …………………… 69
4.19 Equivalent Stress with pi= 4000 psi, pf=10000 psi and ΔT =150oF.................... 70
xii
FIGURE Page
4.20 Tangential Stress with pi= 4000 psi, pf=10000 psi and ΔT =150oF ……….…. 71
4.21 Radial Stress with pi= 4000 psi, pf=10000 psi and ΔT =150oF ……..………... 71
6
4.22 Equivalent Stress for Static Loading with Formation E=1 x 10 psi and v=0.42..... 72
4.25 Fatigue Sensitivity to Life Plot for Zero Based Loading …………….................. 75
4.26 Fatigue Sensitivity to Safety Factor Plot for Zero Based Loading …………....... 75
4.29 Alternating Stress for Zero Based Loading with Formation Property:
6
E=1 x 10 , v=0.3………………………………………………………………... 77
4.30 Alternating Stress for Fully Reversed Loading with Formation Property:
6
E=1 x 10 , v=0.3………………………………………………………………… 78
6
4.31 Life Cycle for Fully Reversed Loading with Formation: E=1 x 10 , v=0.3……… 79
4.33 Alternating Stress for Cement System 2 under a Zero Based Cyclic Loading … 80
4.34 Life Cycle for Cement System 2 under a Zero Based Cyclic Loading………… 81
4.35 Fatigue Sensitivity to Life for a Zero Based Cyclic Loading with Cement 2…. 82
4.36 Fatigue Sensitivity to Damage Plot for Zero Based Loading with Cement 2…. 82
4.37 Fatigue Sensitivity to Safety Factor Plot for a Zero Based Cyclic Loading with
Cement 2….......................................................................................................... 83
xiii
FIGURE Page
4.38 Equivalent Stress for Cement System 1 under Static Loading ……………....... 84
4.39 Alternating Stress for Cement System 1 under Zero Based Cyclic Loading....... 84
4.40 Life Cycle for Cement System 1 under Zero Based Cyclic Loading………….. 85
4.41 Fatigue Sensitivity to Life for Zero Based Cyclic Loading with Cement 1……. 86
4.42 Fatigue Sensitivity to Safety Factor for Zero Based Cyclic Loading with
Cement 1………………………………………………………………………… 86
4.43 Fatigue Sensitivity to Damage for Zero Based Cyclic Loading with Cement 1… 87
xiv
LIST OF TABLES
TABLE Page
3.1 Stress-Life Data for Cement System 1 (Strength; 26 MPa/ 3771 psi)………… 49
3.2 Stress-Life Data for Cement System 2 (Strength; 84 Mpa /12,183 psi)……….. 49
1
CHAPTER I
For a well, whether oil or gas, to maintain its integrity and produce effectively and
completion. This complete zonal isolation, however, can be compromised due to factors
that come into play during the operative life of the completed well. Such factors may
come in the form of thermal or pressure loads generally regarded as HPHT (high
temperature-high pressure) loads which can manifest itself as a static/cyclic load or both,
the number of cycles and even the mechanical properties of the well cement, cyclic
cement.
failure of structures like buildings and bridges but the fatigue behavior of well cement is
still relatively unknown to engineers. Research has led to improved cement designs and
cementing practices, yet many cement integrity problems persist and this further
strengthens the need to understand the mechanism of cement fatigue. Even though most
structural failures are as a result of fatigue rather than static loading, insights on the role
of both static and fatigue loading conditions on the failure of cement sheath would
1.1 Background
A well can be said to have maintained its integrity if it effectively achieves zonal
isolation over its productive life. However, maintaining integrity is not always the case
in real life oilfield practice as case histories abound where the integrity of the well was
compromised due to failure of cement sheath and leading to loss of money and
production.
In order to keep up with the world’s energy demands, oil and gas producing
companies have taken the initiative to explore offshore reserves or drill deeper into
previously existing wells. The consequence of this, however, is that they have to deal
with the high temperatures and pressures encountered at increasing depths. The industry
acknowledges the threshold for high temperature and high pressure conditions as 3000F
and 10,000 psi respectively. For temperatures and pressures above these values, only a
rigorous design would ensure the integrity of the well. High pressure /high temperature
(HPHT) scenarios can be seen in the case of Tuscaloosa trend in Louisiana drilled to a
depth of 23,000 ft and with bottom hole temperature of 4000F and pressures between
17000 to 20000 psi or even in case of the Shearwater field in the East Central Graben
area of the North sea with a depth of 16000 ft, temperature of 3600F and pressure of
sheath as a result of casing pressure testing and formation integrity tests set in motion
events which could compromise the long term integrity of the cement sheath due to
fatigue. Knowledge of the mechanism of fatigue in cement and factors that affect it such
as the magnitude of the load, strength and composition of the cement, mechanical
properties of the cement and pattern of load cycles are important to achieve a realistic
design of a cement system that will be subjected to fatigue loading. Such a design will
go a long way to ensure the long term integrity of a well operating under HPHT
conditions.
• debonding
• radial cracking
These are not new failure modes but just a petroleum engineer’s term for the usual failure
modes in mechanics of materials. Debonding can also be regarded as shear failure and can
exist in two forms - debonding from casing and debonding from formation. It is however
important to note that debonding can also occur as a result of cement shrinkage which in
this case cannot be regarded as a failure due to shear. Radial cracking is a failure mode by
4
fracture as a result of the gradual growth of cracks when the cement is subjected to fatigue
loading. Usually, the surface exhibits no sign of deformation and would finally fail under a
gradually increasing load perpendicular to the loading axis in tension and inclined to the
Cook and Young1 discussed different classical theories of failure for brittle
materials including the maximum normal stress and Mohr-Coulomb’s criteria which may
partly define some of the failure modes described above. These failure criteria are used to
predict if a given material, in this case cement, will fail under a given stress condition.
Concrete under triaxial stresses fail in a unique manner and the Mohr-Coulomb’s criterion
can be used to approximately predict failure when concrete is under compressive and
tensile stresses. Neither this criterion nor the maximum normal stress criterion will suffice
in the case of triaxial compressive state. These failure criteria will be discussed in detail in
Chapter III. It is important to note that these failure criteria work by comparing a given test
condition such as a stress state with an already known parameter of the material like
ultimate strength.
1.1.3 Fatigue
subjected to repeated loading. Fatigue failure due to cyclic loading conditions occurs as
a result of frequently repeated stresses applied at stress levels lower than the ultimate
5
stress under static conditions. The magnitude of this stress decreases as the number of
Since the third stage is very unstable, the fatigue life may be related to the transition
Results of static studies are used as the basis in defining a fatigue study. Fatigue
parameters include
• Fatigue Strength – the greatest amount of stress that can lead to the failure of a
Fig. 1.1 shows the tensile and compressive stress fluctuations for a body under cyclic
loading.
6
Stress
Tension
σa
0
Cycles
σa
Compression
The nature of fatigue in well cement is generally unknown and only a few studies exist on
the fatigue of construction cement. The differences between oil well and cement used in the
construction industry will be discussed in the next chapter. The fatigue strength of
properties of the cement, environmental and loading conditions and water-cement ratio of
concrete. A number of studies have been done on the fatigue of construction cement and it
was found that due to the heterogeneous nature of cement, experimental results show a
large scattering in the concrete behavior due to the cyclic loading and few data sets may not
Studies were conducted by Kim and Kim2 on the fatigue behavior of high strength
concrete using a type I Portland cement to which Elkem micro silica (powder) was added.
A constant minimum stress level of 25 percent of the static uniaxial compressive strength
7
was maintained while the cyclic tests were conducted at maximum stress levels of 75, 80,
85 and 95% of the static strength. The first cycle of loading was loaded at a standard rate,
and the other cycles were loaded in the frequency of 1 Hz. The test results (Fig.1.3)
indicated that under the same stress levels, fatigue life decreases as the concrete strength
increases, and then the fatigue resistance of high strength concrete seems to be inferior to
that of low strength concrete. Figure 1.2 below shows the relationships between the number
of loading cycles and the fatigue strains of concrete low strength (LS) and high strength
(HS). Although the fatigue strain of HS concrete is smaller than that of LS concrete, the
slope of strain increment curve of HS concrete is steeper than that of LS concrete, i.e., the
rate of strain increment increases with the strength of concrete. Therefore high strength
concrete is more brittle than low strength concrete under fatigue loading.
Fig 1.2: Comparison of Fatigue Strain for High and Low Strength Cements [2]
8
Fig 1.3: Comparison of Max Stress Levels to Number of Cycles for Different Cement
Strengths [2]
paste using 2 specimens; one with a high-water cement ratio of 0.7 and another with a
low water-cement ratio of 0.45. A high stress level of 80% (percentage of the ultimate
static compressive strength of the cement) was used in conducting this investigation. It
was observed that the degree to which shrinkage stresses were present in the specimens
were proportional to the water content. This lead to the 0.7 mixture undergoing more
shrinkage due to a more extensive capillary pore system. It was also observed that at
equivalent percentages of the compressive strength, the 0.7 water-cement mixture was
capable of withstanding more cycles to failure than the 0.4 water-cement mixture. From
these results, the author suggested that shrinkage stresses play a greater role in fatigue
strength because they serve to restrain crack propagation3. Crack propagation was slower
in the open capillary structure cement (0.7 mixture) than in the dense structure cement
9
(0.4 mixture) because the high water cement ratio paste is less brittle and can re-adjust
Breitenbucher et al. 4 noticed in their investigations that cyclic loading could lead
to the reduction of stiffness of concrete and that fatigue strain plays a role in the
degradation of the mechanical properties of concrete. The level of damage due to fatigue
loading can be ascertained from the degradation of stiffness at certain number of cycles.
It was also observed that as the longitudinal strain increases at 60% stress level, the
concrete properties (Young’s modulus, fracture energy) decreased faster up to the first
2.0 millions of cycles thereafter, whereas the compressive strength almost remained
constant. There was no observed failure due to fatigue for at least 25.5 million cycles.
Similar observations were made at 70% & 75% stress levels. This shows that the effects
of number of cycles appear to be negligible and therefore the damage is governed only
by the evolution of fatigue strain. These results were also corroborated by the findings of
Breitenbucher and Ibuk5, who in addition noticed that small differences in the upper load
inclusion of steel fibers in ordinary cement paste helped to improve its damage
resistance. From Fig. 1.4, it is observed that at stress levels of 80%, failure occurred at
2.7 x 104 cycles and at a reduced stress level of 72%, the specimens did not fail even
after 2 x 106 cycles. A sharp decline in Young’s modulus during the first 104 cycles was
observed with no appreciable decay in the composite strength. These were also in line
Fig 1.4: Damage Development Depending on Different Cyclic Load Levels [5]
From the existing literatures reviewed, it can be observed that the fatigue of cement can
It should be noted that these studies were conducted using construction cement/concrete
but they should be an indicative of the fatigue of well cement as they both contain the
The objective of this research is to develop a better understanding of the performance of the
casing – cement bond under HPHT well conditions, leading to a model to predict well life.
• to identify the factors that affect the casing-cement integrity under HPHT
conditions
Based on the knowledge acquired from completing these tasks, it will be possible to
analyze the mechanics of casing-cement systems under HPHT conditions for the long term
1.3 Approach
The approach of this research would involve a combination of analytical and finite element
The analytical model of well bore stresses will be developed to provide insights into
the casing –cement –formation interactions under static loading conditions and also
The finite element model will seek to utilize the results from the analytical model to
develop an understanding of fatigue life for cement sheath under HPHT conditions.
12
CHAPTER II
OILFIELD CEMENTING
Well cements are specially formulated for the exploratory drilling of oil and gas wells.
• to enable zonal isolation by preventing the movement of fluids through the annular
space outside the casing, stopping the movement of fluids into fractured formation
or simply by plugging off an abandoned portion of the well. This is the primary
The choice of a particular cement for zonal isolation centers largely on down hole and
formation conditions. But in general, oil well cements should have some features8 which
are necessary for a successful completion job and which would qualify them to be used
• The cement should be able to maintain its integrity in terms of durability and
being free from strength retrogression during the operational life of the well at
• There should be an optimal setting time for the cement .Too reactive a slurry will
result in a short setting time and an insufficiently reactive slurry may take too
long to set
Well cement and construction cement have one thing in common- they are both
Portland cements. Their difference lies in the fact that well cement, in addition to its
Portland cement base, is also mixed up with additives in order to tailor it to a particular
of about 2,600°F to 3,000°F. There are four principal compounds in Portland cement,
which are
result of reaction between water and these constituting components of the cement. The
rate of hydration depends on temperature, size of cement particle and the percentage of
each component present, with C3A hydrating most rapidly followed by C3S and then by
C4AF and finally by C2S. This hydration reaction results in reduction of volume which
makes Portland cement to shrink when set. Expansive cement which is a modified
Portland cement is used to compensate for volume decrease due to shrinkage and to
Pure Portland cement looses its compressive strength and increases its
from the breakdown of its crystalline structure at such temperatures. This would render
the Portland cement unusable for high temperature applications. The strength
When Portland cement is mixed with water, tricalcium silicate (C3S) and
dicalcium silicate (C2S) hydrate to form calcium silicate hydrate (C-S-H) gel and
hydrated lime (Ca(OH)2)2. At temperatures higher than 230°F, C-S-H gel converts to α-
dicalcium silicate hydrate (α-C2SH). Conversion to the α-C2SH phase results in the loss
The American Society for Testing and Materials (ASTM) Specification C-150
provides for eight types of Portland cement I-VIII with type I cement being the normal,
general-purpose cement used for construction purposes10. More than 92% of Portland
cement produced in the United States are type I and II (or Type I/II). Type III accounts
for about 3.5% of cement production. Type IV cement is only available on special
request, and type V may also be difficult to obtain (less than 0.5% of production).
The American Petroleum Institute ,on the other hand, has defined Specifications
for materials and testing for well cements (API Specification 10A) , which includes
requirements for eight classes of oil-well cements (classes A through H) and three grades
(Grades O - ordinary, MSR - moderate sulphate resistant, and HSR - high sulphate
resistant). Each class is applicable for use at a certain range of well depths, temperatures,
15
used in the United States for well cementing. The petroleum industry also uses
conventional types of Portland cement with suitable cement-modifiers. Tables 2.1 and
2.2 shows the ASTM and API cement classifications and their uses.
Special situations call for innovative actions and so, the need to drill deeper and produce
oil and gas under HPHT environments has motivated drilling engineers to come up with
conventional cements that are modified so as to improve its long term sealing integrity in
HPHT wells. In some cases, these designer cements are non-Portland based. The non-
Portland based cements may include Pozzolanic cements, Epoxy Resin11, geo-polymers,
graphite and fibers (glass, steel) which polymerizes at suitable temperatures and/or
cements are not actually cements but react at ordinary temperature with calcium
• Foamed cement,
conditions as the cement for the design of one well may not be appropriate for the design
cement systems in securing a long term integrity for HPHT wells. These have been
conducted using new cement systems which includes a Portland based expansive cement
18
experimental results, the expansive cement exhibited good qualities which make it a
good candidate for a HPHT scenario. When set, it generates an internal compressive
systems. The tensile stress generated by the pressure within the wellbore annulus serves
first to reduce the compressive pre-stress present in the cement before the material
realizes a net tensile stress. As a result, the effective compressive strength of the cement
is increased by the compressive preload applied. This is shown in Fig. 2.1 below:
cement, retarded, with 35% silica) showed the expansive cement system providing a
very good seal in a hard formation. Based on the number of cycles at each pressure
applied (up to 10,000 psi with increments of 1000 psi and temperature of 2000F), the
expansive Cement absorbed 5 times the energy before failure than the conventional
cement system.
From the mechanical shear bond testing, the conventional cement was about 140 psi
times.
The resistance to various materials to allowing water to flow through or past a plug
bond of 3800 psi as opposed to 6000 psi for the expansive cement system.
against confinement and the higher the level of confinement, the better its performance
since expansive cements are generally strong in compression and weak in tension. As a
epoxy resin cement systems but with more inconclusive than positive results as
compared to the expansive cement system. They observed that for epoxy resin, the
20
pumping time is relatively insensitive to pressure and weighting materials, which makes
its design simpler than conventional cement systems. Different hardeners and diluents
can be added to push the temperature higher. The HPHT annular seal testing failed
laboratory investigations but was successful in field trials, and this led the investigators
to conclude that more tests and observation are necessary to understand the mechanism
of sealing using epoxy resin. The mechanical shear bond, however, was both high and
constant, and seems to be due to the material folding up and mechanically resisting the
Foamed cement has also been employed to solve difficult HPHT well integrity
problems. This can be seen in the case of Shearwater field12,13. The presence of higher
than anticipated B annulus pressures in wells drilled at the Shearwater field in the East
central Graben area of the north sea resulted in serious concerns about the long term
integrity of the wells. The field is regarded as a HPHT well with initial reservoir
temperature and pressure of 360oF and 15200 psi respectively at 17900 MD. The
The shearwater field team proposed the use of a foamed cement system citing the
following advantages:
removal properties, has a lesser overall fluid loss and compensates for shrinkage
• Improved Ductility
Compared to conventional systems, foamed cements are more flexible and possess the
ability to withstand both high temperature and high pressure cycling induced stresses.
Foamed cement high tensile strength would make it more resistant to tensile cracking
Under a HPHT scenario, foamed cement provides a cost effective life cycle design even
though the initial cost may be higher than that of conventional cement, and it also
The use of foamed cement system proved very effective in dealing with the
shearwater field problems. However, it was pointed out that a comprehensive analysis is
required to assess the risk of damage to the cement sheath due to downhole well events.
It was also suggested that cement systems should be pre-tested in a laboratory to ensure
Depending on downhole conditions, certain qualities may be required of the cement used
in completing the well. Additives when added to the Portland cement base could be used
to achieve the desired qualities. They could also be used to extend the properties of the
22
base cement. For instance, with additives, Portland cement may be modified to sustain
• Accelerators
These are cement additives that generally tend to reduce the thickening time of cement
slurry and increase the rate of development of compressive strength. Since the hydration
process, which results in the setting of cement, occurs at a faster rate at higher
temperatures, the setting of cement might be a problem while cementing wells drilled in
areas of low temperatures and also result in a lot of waiting time .To counter such,
accelerators like CaCl2, NaCl, sodium silicate, sea water etc are used to speed up the
thickening time.
• Retarders
Retarders are the opposite of accelerators as the name suggests. They act to increase the
thickening time of cement slurry. They are mostly lignosulfonate which are polymers
derived from wood pulp. Examples include calcium, sodium and chemically modified
lignosulfonates.
• Weighting Agents
These are added to cement to increase the density of the final cement mix. They are very
important when designing wells with high temperature and high pressure conditions in
order to give the base cement more strength to sustain high pressures and also to prevent
These set of additives prevents from phase separation under downhole temperature and
pressure conditions. Such a separation would result in fluid being lost to the formation.
• Extenders
These set of additives help to lower the density of the cement mix. Examples include
The drilling and completion of a well is a capital project that runs into millions of dollars
and hence, it is necessary to have a comprehensive design of the cement used for
completion of a particular well and also to avoid remedial cement work which would
add extra cost to the project. Cement design is usually streamlined to a particular well
Ravi and Xenakis 14 discussed a three step approach to cement design. Step one
would involve a detailed engineering analysis. It would require identifying the nature of
the formation- is it a hard or a loose formation? It would require identifying all forces
that would come into play as the well is being produced- are there high temperatures,
high pressures or both? Is it normally or abnormally pressured? Step one also includes
static and fatigue loading analysis to determine if the cement sheath would sustain the
24
series of cyclic loads it would encounter during its lifetime. The answers to step one
questions would lead to step two which would involve designing the cement slurry based
on factors identified in step one. Here properties of the cement like tensile strength,
hydration, and post-cement slurry hydration are chosen so as to effectively match the
The data from the laboratory tests and the analysis of step one are then analyzed
together to evaluate performance. Step three would involve adhering to best drilling and
cementing practices such as centering of casing and effectively cleaning out hole of all
mud so as not to undermine the performance of the designed slurry. It also involves
monitoring during the life of the well. Fig.2.2 below summarizes the design process;
25
The Knowledge of the stress and deformation fields in the vicinity of a HPHT well is very
and high pressure variations exerts excessive load between the protection casing strings and
ultimately on the cement sheath. Studies have been conducted on the effects of high
temperature or high pressure or both on well integrity. Some studies describing such efforts
• Case Study 1
Very high temperatures change the crystalline structure of cement. Stiles15 conducted
cement. Five different cement formulations were exposed to a high temperature of 645oF
and the variations of the Young’s modulus, tensile strength and Poisson ratio of these
cement systems at this temperature (and pressure of 2,133 psi) were observed over a 2
System 1- Conventional (class G) cement system with 40% BWOC silica flour, 4.3%
BWOC gypsum, fluid loss control additive and CaCl2 accelerator. It was mixed at 15.4
lb/gal.
System 2- A low density thixotropic cement system with 92% class A cement with 40%
BWOC silica flour, 8% gypsum, fluid loss control additive and CaCl2 accelerator and
System 3 – A foamed cement system of density of 11.65 lb/gal with class G cement (with
40% BWOC silica flour) as base slurry, a surfactant additive with 26% gas production
System 4 – High strength low density cement of density of 11.65 lb/gal, containing
hollow ceramic microspheres designed with maximized solid volume fraction technique
with silica content in excess of 40% BWOC and with added fluid loss control additive
System 5- A flexible and expanding low density cement system of density of 12.5 lb/gal
containing solid flexible and an MgO base expanding agent designed with maximized
solid volume fraction technique similar to system 4.The flexible particles were added at
50% BVOB.
27
From the analysis of the experimental data, it was observed that the conventional
and foamed cement systems exhibited brittleness after curing at 6450F while other
systems possessed mechanical parameters (low Young’s modulus and high tensile
strength) that are as good as or better than the flexible blends. From this study it can also
be deduced that the ratio of tensile strength to young modulus gives an indication of
resistance to failure under tension. The thixotropic cement and the flexible cement had
The combined presence of high temperatures downhole with high pressure loads lead to
excessive pressure loading in annuli of the casing strings. The pressure load may be as a
result of an increase in pressure around the wellbore region due to pressure integrity
tests, increase of mud weight, casing perforation, stimulation, gas production etc and
these can result in considerable damage to the mechanical properties of cement sheath.
Godwin and Crook16 observed that failure in cement sheath due to excessive pressure
would normally occur in the bottom one half to three-quarters of the casing string and
• Case Study 2
Godwin and Crook16 investigated the effects of high temperatures and excessive annulus
pressure using laboratory experiments and field trials. They circulated hot oil at a
temperature of 350oF through the annulus of the test specimen while gradually
28
increasing the pressure up to 10,000 psi with 2,000 psi increments using the following
cement systems;
System 1- Cement/siliceous material mixture system with 30% BWOW latex, 1.25%
fluid loss control additive and 0.5% gelling agent. It was mixed at 12.1 lbm/gal with
10.81 gal mix water/sack. The yield was 2.49ft3/sack. Compressive strength, Young’s
modulus and Poisson ratio was 1000psi, 0.69 x 106 and 0.42 respectively.
System 2- Cement/Pozzolan mixture system with 30% BWOW latex, 2gal/sack fluid loss
control additive, 10 lbm silica additive and 0.25% gelling agent. It was mixed at 13.1
lbm/gal with 6.48 gal mix water/sack. The yield was 1.76 ft3/sack. Compressive strength,
Young’s modulus and Poisson ratio was 2500psi, 0.8 x 106 and 0.32 respectively.
System 3 is the same as system 2 but without the latex and was mixed at 13.1 lbm/gal.
Compressive strength, Young’s modulus and Poisson ratio was 2000psi, 0.9 x 106 and
0.3 respectively.
System 4- Class H cement system with 35% BWOC silica flour mixed at 18 lbm/gal.
Compressive strength, Young’s modulus and Poisson ratio was 9600psi, 2.4 x 106 and
0.11 respectively.
8000 psi and no appreciable change in the permeability of systems 2 through 4 with
increasing casing pressures. All the cement systems withstood pressures of 2000 psi and
fractures where observed at pressures of 4000 psi with catastrophic effects after 6000 psi
stress cycles. A microanulus occurred in cement 2 with internal pressures up to 6000 psi.
System 1 exhibited elastic properties and provided full casing support until its failure at
29
6000 psi. Field trials were also conducted with systems 2, 3 & 4. After being subjected
to a 12,150 psi internal casing pressure (casing pressure plus casing fluid hydrostatic
pressure) system 4 failed while system 3 remained intact. No evaluation logs were run
for system 2.
As the awareness of the need to analyze the structural behavior of the cement is
currently increasing, so does the need for guidelines and quantitative results. This has led
to the use of simple models or the more complex finite element models to analyze the
casing-cement-formation system and to base the cement selection on the results of this
analysis. Fleckenstein et al.17 proposes finite element analysis as the best way to analyze
the casing-cement-formation properties during the design phase. This proposal has been
corroborated by others18,19,20,21.
and solving the problem using FEA. With FEA, stress variations at different points along
the cement sheath radius can be effectively modeled. Several studies 18,19,20,21 show how
finite element analysis was effectively used to model HPHT wells in order to understand
and overcome cement sheath failure problem. Ravi et al22 sounds a note that cement
which may be suitable under one set of conditions may not be suitable under a different
30
set of conditions. Thus, a rigorous design procedure using finite element analysis should
be used to select the cement system that will satisfy each specification.
Fleckenstein et al.17 used finite element analysis to investigate the role which mechanical
properties of the cement plays in withstanding the internal casing pressure by modeling a
cemented wellbore with two cement systems- a soft and hard cement, the soft (ductile)
cement being one with a high Poisson’s ratio and low modulus of elasticity and the hard
(brittle) cement being one with a low Poisson’s ratio and relatively high modulus of
elasticity.
The hard cement was Class H cement mixed with 35% silica flour at 18 ppg,
with the following mechanical properties: compressive strength: 9,500 psi, Young’s
Modulus: 2,400,000 Psi, Poisson’s Ratio: 0.11 while the soft cement system was a
cement/siliceous material mixed with 30% latex at 12.1 ppg with Compressive strength:
1,000 Psi, Young’s Modulus: 690,000 psi, Poisson’s Ratio: 0.42. These cement systems
were also studied by Godwin and Crook16 .The results indicated that there is little
confining stress outside the cement sheath would increase the burst resistance of the
casing.
31
Fig 2.3: Tangential Stress for Hard and Soft Cement Systems [17]
The confining stress on the cement sheath may come from the formation and its effect on
the cement sheath was also modeled. It was run with a 10,000 psi internal burst pressure
acting upon the inner surface of the 5-½” casing and a 3000 psi confining stress from the
formation.
The results show reductions in von Mises stress from 6,099 psi to 4,292 psi and
in tangential stress from 894 Psi (tension) at the casing outer diameter to -283 Psi
(compression) at the borehole wall. The investigators noticed a difference in the von
Mises stress that is generated by the hard and soft cement systems. At 1000 psi, hard
cement generate twice the von Mises stress when compared with the soft cement but this
reduces as the confining stress increases with a reduction of less than 20% at 5000psi
confining stress.
32
From these investigations they deduced that radial cracking is less likely to occur
with soft cement systems because hard cement systems are likely to generate significant
tangential stresses (Fig. 2.3), which increases the likelihood of forming radial cracks in
It should also be noted23 that with the use of flexible and expanding cement
systems, the stiffness of the formation plays an important role. If the stiffness of the
formation is low (low Young’s modulus) compared to the cement system, debonding and
CHAPTER III
3.1.1 Background
infinitesimal element within the cement sheath, these stresses act in a three dimensional
fashion and the cement sheath can therefore be regarded as being under a triaxial stress
state, as shown in Figure 3.2. The third component (not shown in the diagram) is the
axial stress component which is perpendicular to the two stresses shown in Figure 3.2.
The radial stress is always compressive in nature while the tangential stress could be
tensile or compressive depending on the loading conditions. Figure 3.1 shows the nature
and profile of tangential and radial stresses under different loading conditions.
cylinder with three concentric cylinders and a perfect bonding is assumed to exist
between the cement and casing and the cement and formation. The pressure and/or
temperature changes induce stress concentrations near the casing – cement and the
cement-formation boundaries.
34
+
σr σr σθ
- σθ - -
(a)Internal pressure
(b)External pressure
Fig 3.1: Radial and Hoop Stress Profile Due to Loading Conditions [1]
Formation
Cement
Casing
3.1.2 Assumptions
The following assumptions were made to facilitate the analytical modeling of the
wellbore:
• The composite cylinder undergoes plane strain deformation. This implies that the
discontinuities. This implies that the radial displacements and radial stresses are
• The cement sheath and formation is treated as a thick walled pressure vessel
The failure criteria employed in predicting the failure of cement sheath in this model
includes
• Mohr-Coulomb’s criterion
• Experimental Investigation
The maximum normal stress criterion predicts that an isotropic material will fail when the
largest principal stress reaches a limiting value. This implies that failure would occur when
36
σ1
≥1 (3.1)
σf
where σ 1 is the maximum principal stress and σ f is the limiting stress. If σ 1 is tensile,
then σ f is the limiting tensile stress and the other two smaller principal stresses σ 2 & σ 3
becomes
σ3
≥1 (3.2)
σf
where σ 3 is the magnitude of the minimum principal stress. It is should be noted that
this criterion would be inaccurate if all three principal stresses are compressive.
In the Mohr-Coulombs’ criterion, σ 2 does not to play any role and failure is predicted to
occur when
σ1 σ3
− ≥1 (3.3)
σ tensile σ compressive
where σ tensile and σ compressive are the tensile and compressive strengths respectively.
approximation when the tangential, radial and axial stresses are all compressive. In such
cases, different failure theories based on experiments would be required to determine the
onset of failure. Avram et al.24 discussed concrete fracture under triaxial stresses and
proposed a new failure criterion given in Eq. (3.4) below which complies with Mohr-
Coulombs’ criterion.
37
0.86
σ1 ⎛σ ⎞
= 1 + 3.7⎜⎜ 3 ⎟⎟
fc ⎝ fc ⎠ (3.4)
where f c is the compressive strength of the cement, σ 1 is the major principal compressive
stress at failure and σ 3 is the minor principal compressive stress with σ 1 〉 σ 2 〉 σ 3 . Fig 3.3
and Fig 3.4 shows the experimental results and failure envelope respectively for concrete
Fig 3.3: Concrete Failure Criterion under Triaxial Compressive Stresses [24]
38
Fig 3.4: Failure Envelope for Triaxial Compression and Tensile Stress State [24]
It should be noted that this experimental failure criterion was proposed for plain concrete
In a composite cylinder model under consideration, the internal pressure pi acting on the
inner surface of the casing in conjunction with temperature increase will expand the
casing radially , while the cement sheath will resist the expansion. As a result, a contact
pressure (pc1) will develop at the interface between the casing and the cement.
39
Considering the casing-cement interface as shown in Fig. 3.5, pc1 is the contact
pc2
b pc1
a pi
pc1
1
εθ = [σ θ − ν (σ z + σ r )] + αΔT (3.5)
E
1
εz = [σ z −ν (σ θ + σ r )] + αΔT (3.6)
E
But since the axial strain is negligible considering the large depth, then εz≈0 (i.e. plane
σ z = ν [σ r +σ θ ]− α E Δ T (3.7)
40
εθ =
1
E
[ ( ) ( )
σ θ 1 − ν 2 − ν + ν 2 σ r + (1 + ν )αEΔT ] (3.8)
δr =
r
E
[ ( ) ( )
σ θ 1 − ν 2 − ν + ν 2 σ r + (1 + ν )αEΔT ] (3.9)
Let the radii a, b and c be represented by ra , rb & rc respectively. Considering the casing
prm
σ r = − p and σ θ = (3.10)
ts
where p = pi − pc1 , r m is the mean radius of the casing and ts is the thickness of the
a (Pi − Pc1 ) ⎡ rm ⎤
δ r − ca sin g = { ⎢ ( 1 −ν s
2
) + (ν s + ν s2 )⎥ } + [ (1 + ν s )aα s ΔT ] (3.11)
Es ⎣ ts ⎦
∂ (ΔT )
Considering the cement sheath as a thick walled cylinder and assuming = 0 , the
∂r
pc1b 2 ⎛ c 2 ⎞ pc 2c 2 ⎛ b 2 ⎞
σr = ⎜1 − ⎟ − ⎜1 − ⎟ (3.12)
c 2 − b 2 ⎜⎝ r 2 ⎟⎠ c 2 − b 2 ⎜⎝ r 2 ⎟⎠
pc1b 2 ⎛ c 2 ⎞ pc 2c 2 ⎛ b 2 ⎞
σθ = ⎜1 + ⎟ − ⎜1 + ⎟ (3.13)
c 2 − b 2 ⎜⎝ r 2 ⎟⎠ c 2 − b 2 ⎜⎝ r 2 ⎟⎠
σ r = − pc1 (3.14)
41
⎛ c2 + b2 ⎞ ⎛ 2c 2 ⎞
= ⎜ ⎟
σ θ pc1 ⎜ 2 2 ⎟ pc 2 ⎜⎜ 2 2 ⎟⎟
− (3.15)
⎝c −b ⎠ ⎝c −b ⎠
This gives the radial expansion in the cement sheath at r = b when Eq. (3.14) and Eq.
b ⎡ ⎡ ⎛ b2 + c2 ⎞ ⎛ 2c 2 ⎞⎤ ⎤
δ r − cement = ⎢( 1 − ν c )⎢ pc1 ⎜ 2
2
⎜ ⎟
2 ⎟
− p ⎜
2⎜ 2
⎟
2 ⎟ ⎥ + p (
ν + ν 2
)⎥ + [ (1 + ν c )bα c ΔT ]
⎣ ⎝c −b ⎠ ⎝ c − b ⎠⎦
c c 1 c c
Ec ⎣⎢ ⎦⎥
(3.16)
Since both radial expansions are equal, it follows from Eq. (3.11) & Eq. (3.16) that
⎧⎪ b
( )
2 ⎡b + c ⎤
( a ⎡ rm
) ( 2 ⎤⎫
) (⎡b
) ⎛ 2c 2 ⎞
( ⎤
)
2 2
⎪
pc1 ⎨ [ 1 − vc ⎢ 2 ]
2
⎥ + v + v c + ⎢ 1 − v
2
+ v + v s ⎥ ⎬ − pc 2 ⎢ ⎜⎜ 2 ⎟ 1 − vc 2 ⎥ =
2 ⎟
⎣c − b ⎦ ⎝c −b ⎠
2
⎪⎩ Ec ⎦ ⎪⎭
c s s
Es ⎣ t s ⎣ Ec ⎦
pi a ⎡ rm
Es ⎣ t s
(
2
) (
2 ⎤
)
⎢ 1 − vs + vs + v s ⎥ + [ (1 + ν s )aα s ΔT ] − [ (1 + ν c )bα c ΔT ]
⎦
(3.17)
where
⎧⎪ b
A = ⎨ [ 1 − vc (
2
)
⎡ b2 + c2 ⎤
⎢ 2 2⎥
+ vc + v ( 2
c ) ] + Ea ⎡⎢ rt (1 − v ) + (v + v )⎤⎥⎫⎪⎬⎪ (3.19)
m 2 2
s
⎪⎩ Ec ⎣c − b ⎦
s s
s ⎣ s ⎦⎭
⎡ b ⎛ 2c 2 ⎞ ⎤
B = - ⎢ ⎜⎜ 2 ⎟ 1 − vc 2 ⎥
2 ⎟
( ) (3.20)
⎣ Ec ⎝ c − b ⎠ ⎦
C= (
pi a ⎡ rm 2
) (
⎢ 1 − v s + vs + v
Es ⎣ t s
2
s )⎤⎥ + [ (1 +ν )aα ΔT
s s ] − [ (1 + ν c )bα c ΔT ] (3.21)
⎦
42
Pf
pc2
c
Pc2
b
pc1
∂ (ΔT )
Considering the cement sheath and assuming = 0 , at r = c,
∂r
σ r = − pc 2 (3.22)
⎛ 2b 2 ⎞ ⎛ c2 + b2 ⎞
σ θ = pc1 ⎜⎜ 2 2 ⎟⎟ − pc 2 ⎜⎜ 2 2 ⎟⎟ (3.23)
⎝c −b ⎠ ⎝c −b ⎠
43
This gives the radial expansion in the cement sheath at r = c when Eq. (3.22) and Eq.
c ⎡ ⎡ ⎛ 2b 2 ⎞ ⎛ c 2 + b 2 ⎞⎤ ⎤
δ r − cement = ⎢( 1 − ν c )⎢ pc1 ⎜⎜ 2
2
⎟
2 ⎟
− p ⎜
2⎜ 2
⎟
2 ⎟ ⎥ + p 2
ν + ν 2
(
⎥ + [ (1 + ν c )cα c ΔT ) ]
⎣ ⎝c −b ⎠ ⎝ c − b ⎠⎦
c c c c
Ec ⎣⎢ ⎥⎦
(3.24)
Considering the formation as a thick walled pressure vessel with a finite radius d into the
∂ (ΔT )
formation and also assuming = 0 , at r = c
∂r
σ r = − pc 2 (3.25)
⎛ c2 + d 2 ⎞ ⎛ 2d 2 ⎞
σ θ = pc 2 ⎜⎜ 2 2 ⎟⎟ − p f ⎜⎜ 2 2 ⎟⎟ (3.26)
⎝d −c ⎠ ⎝d −c ⎠
when Eq. 3.24 and Eq. 3.25 are substituted into Eq. 3.9 , it follows that
c ⎡ ⎡ ⎛ c2 + d 2 ⎞ ⎛ 2d 2 ⎞⎤ ⎤
δ r − formation = ⎢( 1 − ν f )⎢ pc 2 ⎜⎜ 2
2
⎟
2 ⎟
− p ⎜
f⎜ 2 2 ⎟⎥
(
⎟ + pc 2 ν f + ν 2f ⎥ + [ (1 + ν f )cα f ΔT ) ]
E f ⎣⎢ ⎣ ⎝d −c ⎠ ⎝ d − c ⎠⎦ ⎥⎦
(3.27)
Since both radial expansions are equal, it follows from Eq. 3.24 & Eq. 3.27 that
[ (1 − v f 2 )⎢ d 2 + c 2 ⎥ + (v f ) ] + Ec [1 − v ]⎛⎜⎜ bc ( )⎫⎪⎬⎪ − p
⎧⎪ c ⎡c ⎤
⎡ ⎤ + c2 ⎞ ⎛ 2b 2 ⎞
( )
2 2 2
2 2
pc 2 ⎨ +v
2
⎟ − vc + v c1 ⎢
⎜⎜ 2 ⎟ 1 − vc 2 ⎥ =
− b 2 ⎟⎠ 2 ⎟
f c
⎣d − c ⎦ ⎝c −b ⎠
c 2
⎪⎩ E f c ⎝ ⎭ ⎣ Ec ⎦
⎡ c
pf ⎢
⎛ 2d 2 ⎞
⎜⎜ 2 (
⎟ 1 − vf 2
2 ⎟
)⎤⎥ − [ (1 + ν )cα ΔT ] + [ (1 + ν c )cα c ΔT ]
⎝d −c ⎠
f f
⎢⎣ E f ⎥⎦
(3.28)
⎡ c ⎛ 2b 2 ⎞ 2 ⎤
where D = - ⎢ ⎜⎜ 2 ⎟
2 ⎟
1 − v (⎥ ) (3.30)
⎣ Ec ⎝ c − b ⎠
c
⎦
⎡ p c ⎛ 2d 2 ⎞
F = ⎢ f ⎜⎜ 2 ⎟ 1− vf 2
2 ⎟
( )⎤⎥ − [ (1 + ν )cα ΔT ] + [ (1 + ν c )cα c ΔT ] (3.32)
⎢⎣ E f ⎝ d − c ⎠
f f
⎥⎦
From a simultaneous solution of Eq.(3.18) and Eq. (3.29), the contact pressures pc1 and
FB − KC
pc1 = (3.33)
DB − AK
⎡ FB − KC ⎤
C−⎢ A
⎣ DB − AK ⎥⎦
pc 2 = (3.34)
B
From the analysis presented above, the circumferential, radial and axial stresses present
pc2
pc1
The radial, tangential and axial stresses in the cement sheath are then calculated using
b2 ⎡ c2 ⎤ c2 ⎡ b2 ⎤
σ r −cement = Pc1 ⎢1 − 2 ⎥
− Pc2 2 ⎢1 − 2 ⎥ (3.35)
c2 − b2 ⎣ r ⎦ c − b2 ⎣ r ⎦
b2 ⎡ c2 ⎤ c2 ⎡ b2 ⎤
σ θ −cement = Pc1 ⎢1 + 2 ⎥
− Pc2 ⎢1 + 2 ⎥ (3.36)
c2 − b2 ⎣ r ⎦ c2 − b2 ⎣ r ⎦
(Pc1 − Pc 2 )b 2 c 2
τ max = (3.38)
(c 2
− b 2 )r 2
At r = b,
(Pc1 − Pc 2 ) c2
τ max = (3.39)
(c 2
− b2 )
At r = c
τ max =
(Pc1 − Pc 2 ) b 2 (3.40)
(c 2
− b2 )
The response of the cement sheath to different static and fatigue loading conditions will be
• Cement system 1 - Ductile cement system with compressive strength of 3000 psi,
tensile strength of 1000 psi, a young’s modulus of 0.69 x 106 psi and a Poisson ratio
of 0.4
46
• Cement system 2 - Brittle cement system with compressive strength of 9500 psi,
tensile strength of 3000 psi, a young’s modulus of 2.4 x 106 psi and a Poisson ratio
of 0.1.
• Cement system 3- A low young’s modulus and a low Poisson ratio cement system
with compressive strength of 2500 psi, tensile strength of 1000 psi a young’s
The finite element analysis was done with ANSYS workbench 11.0 and since the
casing –cement –formation is axisymetric, a quadrant was used for the 2D modeling.
47
element simulations were compared. An internal pressure of 15,000 psi was applied inside
the casing with no formation pressure. The meshing was done with 6648 elements and
20,607 node density. Plane strain state was assumed. The boundary conditions applied
include
• U=0, τxy = 0 on X = 0
• V=0 , τyx= 0 on Y= 0
• σ rr = pi on r = a
• σ rr = p f on r = d
The casing –cement –formation model was also assumed fully bonded with no separation at
the interfaces.
Figure 3.8 shows the distribution of equivalent von Mises stress under the loading
condition described above and Figure 3.9 below compares the analytical values for von
Mises equivalent stress with those from the finite element analysis. The analytical and
Fig 3.8: Equivalent Stress for Casing –Cement –Formation Model with Meshing
9000
8500
8000
Equivalent stress
7500
FEA
7000
Analytical
6500
6000
5500
5000
4.5 4.7 4.9 5.1 5.3 5.5 5.7 5.9
Radius(inches)
Fig 3.9: Equivalent Stress Comparison For Analytical and FEA Models
49
To study the fatigue of the cement sheath, finite element analysis was utilized. The fatigue
properties of concrete were used for this analysis. Such properties include S-N curves,
strain-life curves and stress-strain curve. Tables 3.1 and 3.2 below gives the fatigue
material data used in the finite element modeling for cement systems 1 and 2 based on data
extracted from Fig.1.3 using a special software. The fatigue life for the cement systems will
be predicted from these curves using ANSYS. The S-N curve is shown in Fig.3.10.
Table 3.1: Stress-Life Data for Cement System 1 (Strength; 26 MPa/ 3771 psi)
No of Cycles Stress(psi)
0.99692 4332.322
4.65426 4031.543
57.8064 3590.166
683.674 3203.529
8084.64 2852.871
100398 2535.519
1.00E+06 2276.105
Table 3.2: Stress-Life Data for Cement System 2 (Strength; 84 Mpa /12,183 psi)
No of Cycles Stress(psi)
1.01527 13376.842
1.87358 13009.757
19.6824 11594.051
165.699 10353.456
1740.34 9208.586
16164.6 8206.784
1.46E+05 7314.037
50
Newman and Choo25 gave a relationship between the number of cycles to failure
and load ratio derived through experimental studies for concrete with densities greater
f c. max
= 1 − 0.0685(1 − R) log10 N 3.26 (3.41)
fc
f c. max
R is the stress ratio and
f c. min
N is the number of cycles to failure.
The fatigue properties of the cement sheath were also studied under two loading conditions;
These loading conditions shown in Fig. 3.11 below, presents a cyclic loading condition due
to shear at the coupling and that due to the internal and formation pressure on areas where
Casing
τ
Pi
Casing Po
Coupling
Cement
For both cases, a constant amplitude loading condition is assumed .A fully reversed
loading is also assumed while modeling the fatigue behavior of the cement under shear
cyclic loading while a zero based (compression) loading is assumed for the internal and
external pressures. The Goodman diagram for fatigue analysis was chosen as it gives a
good description of the fatigue behavior of brittle materials. Figure 3.12 shows the
Goodman diagram. Figs. 3.13 and 3.14 present the different loading conditions and options
Fig 3.12: Goodman Diagram for Brittle and Ductile Materials [1]
σ σm 1
a
+ = (3.42)
σ fs σu SF
53
where σm =
(σ max + σ min ) (3.43)
2
σa =
(σ max − σ min ) (3.44)
2
Fig 3.13: Zero Based Loading, Goodman’s Diagram and Fatigue Options
54
Fig 3.14: Fully Reversed Loading, Goodman’s Diagram and Fatigue Options
• Fatigue life - The fatigue life plot gives an indication of number of cycles to failure
of a particular material;
• Damage - The fatigue damage plot gives the fatigue damage at a given design life;
55
• Factor of safety – This gives the factor of safety with respect to fatigue failure at a
• Fatigue sensitivity – This plot gives an indication of how loading conditions affect
CHAPTER IV
To understand the effect of static loading on the integrity of cement sheath, the analytical
model and finite element analysis were used to examine the responses of different cement
described in the previous chapter. The temperature change was neglected and the analytical
model was used to show the trends in the von Mises equivalent, tangential and axial
With an inner pressure of 15,000 psi and a formation pressure of 1000 psi, the finite
element model gives us an idea of the amount of stress generated within the steel casing
and the actual amount transferred to the cement sheath through the casing-cement interface.
The same can also be said of the transfer of stresses from the formation to the cement
sheath through the cement-formation interface. Figures 4.1 and 4.2 shows the trend for von
Mises stress using cement system 2, which has the same mechanical properties for steel
casing but with different formation properties. From the diagrams, it is evident the amount
of support the steel casing provides the cement sheath. In Fig. 4.1, the equivalent stress
within the casing decreases non-linearly from 84,178 psi to 66,784 psi. At the casing-
57
cement interface, there is a large decrease from 66,784 to 8,384 psi. Due to the fact that the
properties of the cement (E=2.4 x 106 psi, v=0.4) and the formation (E=3 x 106 psi, v=0.42)
are similar, there is a little change in the von Mises equivalent stress across the cement-
formation interface (5,806 to 5736 psi). Fig. 4.2 shows somewhat a similar trend but with
higher stresses imposed on the steel casing (93,575 to 74,325 psi) and a lower stress on the
cement sheath (8,160 to 5,560 psi).As opposed to Fig. 4.1, there a significant difference in
Fig 4.1: Equivalent Stress with Formation E=3 x 106 psi and v=0.42 for Scenario 1
58
Fig 4.2: Equivalent Stress with Formation E=1 x 106 psi and v=0.3 for Scenario 1
Keeping the formation and casing properties unchanged, the analytical model was
used to study the response to a high internal pressure and low formation pressures on the
casing-cement-formation system as shown in Figs. 4.3 to 4.5. The results show that cement
system 1 generates compressive (negative) radial and tangential stresses. The tangential
stress fluctuates between a maximum value of -1,161 psi to a minimum of -590 psi. In
contrast, cement systems 2 and 3 generate tensile tangential stresses. Cement system 3
posses more of a ductile quality and generates less von Mises stress than system 2 but a
10000
9000
8000
Equivalent Stress(psi)
7000
6000
cement system 1
5000 cement system 2
4000 cement system 3
3000
2000
1000
0
4.5000 5.0000 5.5000 6.0000
Radius(inches)
Fig 4.3: Equivalent Stress in Three Cement Systems with pi=15000 psi and pf=1000 psi
6000
5000
4000
Tangential stress(psi)
3000
cement system 1
2000 cement system 2
cement system 3
1000
0
4.5000 5.0000 5.5000 6.0000
-1000
-2000
Radius(inches)
Fig 4.4: Tangential Stress in Three Cement Systems with pi=15000 psi and pf=1000 psi
60
600
-400
4.5000 5.0000 5.5000 6.0000
-1400
Radial stress(psi)
-4400
-5400
-6400
Radius(inches)
Fig 4.5: Radial Stress in Three Cement Systems with pi =15000 psi and pf =1000 psi
The trend in the radial and tangential stresses in the three cement systems subjected to a
high inner pressure of 15,000 psi and a zero formation pressure is strikingly similar to the
trend in scenario 1 with values for radial and tangential stresses being a little higher as
shown in Figs. 4.7 to 4.9.The maximum equivalent stress in this case is 8,692 psi for
cement system 2 as opposed to 9,000 psi for the same cement system under scenario 1.The
same can be said about the values for the radial and tangential stresses for these two
loading cases.
61
Fig 4.6: Von Mises Stress with Formation E=1 x 106 psi and v=0.3 for Loading Case 2
The finite element modeling of case 2 with cement system 2 (Fig. 4.6) shows a
higher amount of equivalent stress on the casing (91,151 to 73,398 psi) in comparison with
case 1 with the same formation properties (84,178 to 66,784 psi).The transitions at the
boundary show a different trend from that of scenario 1.At the casing-cement interface, the
pressure reduces from 73,398 psi in the casing to 8,210 psi in the cement .The trend is the
opposite at the cement-formation boundary with 5,586 psi at the cement side of the
10000
9000
7000
6000
cement system 1
5000 cement system 2
4000 cement system 3
3000
2000
1000
0
4.5000 5.0000 5.5000 6.0000
Radius(inches)
Fig 4.7: Equivalent Stress in Three Cement Systems with pi=15000 psi and pf=0 psi
6000
5000
4000
Tangential stress(psi)
1000
0
4.5000 5.0000 5.5000 6.0000
-1000
Radius(inches)
Fig 4.8: Tangential Stress in Three Cement Systems with pi=15000 psi and pf=0 psi
63
600
4.5000
-400 5.0000 5.5000 6.0000
-1400
Radial stress(psi)
cement system 1
-2400 cement system 2
cement system 3
-3400
-4400
-5400
Radius(inches)
Fig 4.9: Radial Stress in Three Cement Systems with pi=15000 psi and pf=0 psi
Fig 4.10-4.12 below presents the case when a casing pressure of 4000 psi and a high
formation pressure of 10,000 psi is applied to the casing –cement-formation model. The
trend here is significantly different t from that seen in the first two cases. Here, both the
tangential and radial stresses are compressive in nature. Cement systems 2 and 3 generates
almost equal amounts of tangential stresses and are much lower than those of cement
system 1.However, cement system 1 has the lowest value of radial stress (-5297 psi) of the
7000
6000
5000
Equivalent Stress(psi)
4000
cement system 1
cement system 2
3000 cement system 3
2000
1000
0
4.5000 5.0000 5.5000 6.0000
Radius(inches)
Fig 4.10: Equivalent Stress in Three Cement Systems with pi=4000 psi and pf=10000 psi
0
4.5000 5.0000 5.5000 6.0000
-1000
-2000
Tangential stress(psi)
-5000
-6000
-7000
Radius(inches)
Fig 4.11: Tangential Stress in Three Cement Systems with pi=4000 psi and pf=10000 psi
65
500
300
100
-100
4.5000 5.0000 5.5000 6.0000
Radial stress(psi)
-700
-900
-1100
-1300
Radius(inches)
Fig 4.12: Radial Stress in Three Cement Systems with Pi=4000 psi and Pf=10000 psi
It is no doubt that temperature change plays an important role in adding to the stresses that
would ultimately lead to the failure of the cement sheath. As expected, the combined effect
of temperature and pressure shows a trend that is a bit different, in terms of the magnitude
of the stress, from that when the effects of pressure are considered alone. Figs. 4.13 to 4.15
give us a view of that trend when an inner pressure of 15,000 psi, formation pressure of
1000 psi and temperature change (ΔT) of 150oF respectively are applied to the casing-
cement-formation model.
66
12000
8000
cement system 1
6000 cement system 2
cement system 3
4000
2000
0
4.5000 5.0000 5.5000 6.0000
Radius(inches)
Fig 4.13: Equivalent Stress with pi=15000 psi, pf=1000 psi and ΔT =150oF
6000
5000
4000
Tangential stress(psi)
3000
cement system 1
2000 cement system 2
cement system 3
1000
0
4.5000 5.0000 5.5000 6.0000
-1000
-2000
Radius(inches)
Fig 4.14: Tangential Stress with pi=15000 psi, pf=1000 psi and ΔT =150oF
67
600
-400
4.5000 5.0000 5.5000 6.0000
-1400
Radial stress(psi)
-4400
-5400
-6400
Radius(inches)
Fig 4.15: Radial Stress with pi=15000 psi, pf=1000 psi and ΔT =150oF
With the effect of temperature change, the trend in the tangential stress distribution
across the cement sheath remains pretty much the same although the stress value is higher.
The highest von Mises stress with ΔT=0 is about 9,000 Psi (cement system 2) compared
with 10,000 psi with the same cement system with ΔT=150. The tangential stress
12000
10000
Equivalent Stress(psi)
8000
cement system 1
6000 cement system 2
cement system 3
4000
2000
0
4.5000 5.0000 5.5000 6.0000
Radius(inches)
Fig 4.16: Equivalent Stress with pi=15000 psi, pf=0 psi and ΔT =150oF
7000
6000
5000
Tangential stress(psi)
4000
cement system 1
3000 cement system 2
cement system 3
2000
1000
0
4.5000 5.0000 5.5000 6.0000
-1000
Radius(inches)
Fig 4.17: Tangential Stress with pi=15000 psi, pf= 0 psi and ΔT =150oF
69
600
4.5000
-400 5.0000 5.5000 6.0000
-1400
Radial stress(psi)
cement system 1
-2400 cement system 2
cement system 3
-3400
-4400
-5400
Radius(inches)
Fig 4.18: Radial Stress with pi=15000 psi, pf= 0 psi and ΔT =150oF
In this case, the radial stress distribution profile is similar to that of Fig. 4.9 with
almost the same magnitude of stress but the tangential stress profile is different. While the
tangential stress profile in Fig. 4.8 was all negative (compressive) stress for cement system
1, it fluctuates from 85 psi (tensile) to -468 psi (compressive) when the effect of
temperature change is considered.. The equivalent stresses decreases non -linearly from the
inner surface of the cement ( 2928 psi for cement 1, 9,652 psi for cement 2, 4,687 psi for
cement 3) to 1,977 psi for cement 1, 6,802 psi for cement 2, 3,164 psi for cement 3 at the
This case generates the higher values for tangential and radial stresses in the three cement
systems as compared to the case without temperature change. The radial stress in the three
cement systems also increases non-linearly from the inner surface of the cement to the
outer surface where the opposite is the case without temperature change. The results are as
6000
5000
Equivalent Stress(psi)
4000
cement system 1
3000 cement system 2
cement system 3
2000
1000
0
4.5000 5.0000 5.5000 6.0000
Radius(inches)
Fig 4.19: Equivalent Stress with pi=4000 psi, pf=10000 psi and ΔT =150oF
71
0
4.5000 5.0000 5.5000 6.0000
-1000
-2000
Tangential stress(psi)
-5000
-6000
-7000
Radius(inches)
Fig 4.20: Tangential Stress with pi=4000 psi, pf=10000 psi and ΔT =150oF
600
-400
4.5000 5.0000 5.5000 6.0000
-1400
-2400
Radial stress(psi)
-3400
cement system 1
-4400 cement system 2
cement system 3
-5400
-6400
-7400
-8400
-9400
Radius(inches)
Fig 4.21: Radial Stress with pi=4000 psi, pf=10000 psi and ΔT =150oF
72
The aim of the fatigue study was to determine if a cement system capable of sustaining a
static load would be able to sustain a similar load under cyclic loading conditions. A
fatigue analysis together with a related static analysis was done on a casing – cement -
formation setup with ANSYS using cement systems 2 with an inner pressure of 7,000 psi
and formation pressure of 2,000 psi. As was done in the static study, the properties of the
casing remained the same while the fatigue behavior of cement was studied with two
different formation properties; E=3 x 106 psi, v=0.42 and E=1 x 106 psi, v=0.3. Fully
reversed and zero based loading conditions were used in the fatigue analysis.
Fig 4.22: Equivalent Stress for Static Loading with Formation E=1 x 106 psi and v=0.42
73
With 7,000 psi internal pressure and 2000 psi formation pressure, the finite element
model in Fig. 4.22 gives an indication that cement system 2 would withstand such pressure
loads under static loading conditions. The effective von Mises stress changes from 22,760
psi at the cement-casing boundary on the casing side to 4330 psi on the cement sheath side
and from 2994 psi on the cement side of cement-formation boundary to 1,842 psi on the
formation side of the same boundary. The equivalent alternating stresses for a zero based
cyclic loading condition with the formation properties E=1 x 106 psi, v=0.42, however, is a
The zero based fatigue loading results in a lower equivalent alternating stress in the
casing but higher values in the cement and the formation, as shown in Fig. 4.23.The
equivalent alternating stress decreases non-linearly from 18,199 psi to 13,719 psi at the
casing-cement boundary with a life of 5.86 x 105 cycles (Fig. 4.24). The cement part of the
casing –cement boundary has a life of 5.6 x 104 cycles (Fig.4.24) and an alternating stress
of 7,783 psi (Fig 4.23) psi as opposed to 4,330 psi (fig 4.22) under static loading
conditions. The progression from the cement to the formation at the cement –formation
boundary however sees a large increase in the equivalent alternating stress from 2,992 psi
to 1 x 1032 psi.
Figures 4.25 to 4.27 show the fatigue sensitivity plot for life, damage and safety
factor but does not explicitly show the fatigue performance of the different regions of the
Fig 4.25: Fatigue Sensitivity to Life Plot for Zero Based Loading
Fig 4.26: Fatigue Sensitivity to Safety Factor Plot for Zero Based Loading
76
Fig 4.27: Fatigue Sensitivity to Damage Plot for Zero Based Loading
With the formation properties changed to E=1 x 106 psi and v=0.3, while the casing
and cement properties remain the same, an insight is gained on how mechanical properties
like Young’s modulus and Poisson ratio affects stress distribution and ability to withstand
fatigue loading. Figure 4.27 shows the response to fatigue with zero based load and
formation properties E=1 x 106 psi and v=0.3.Under the different formation properties, both
show similar trends in the static equivalent stresses and alternating stresses but as shown in
Figs. 4.28 and 4.29, the formation with more ductile properties (E=1 x 106 psi and v=0.3)
has slightly higher stress values for the static equivalent stresses and almost the same value
alternating stress.
77
Fig 4.28: Equivalent Stress for Static Loading with Formation Property: E=1 x 106, v=0.3
Fig 4.29: Alternating Stress for Zero Based Loading with Formation Property: E=1 x 106,
v=0.3
78
FEA modeling of the fully reversed cyclic shear loading of cement sheath at the
couplings (Figs. 4.30 and 4.31) also showed that the casing-cement boundary, which is the
Fig 4.30: Alternating Stress for Fully Reversed Loading with Formation Property:
Fig 4.31: Life Cycle for Fully Reversed Loading with Formation Property: E=3 x 106,
v=0.3
The fatigue sensitivity plots for life, damage and safety factor provides us with the
number of life cycles to failure, the amount of damage done and the factor of safety for the
casing-cement-formation system, but they do not clearly express the effects of fatigue
loading on the cement sheath. In order to achieve this, the cement sheath was isolated and
the contact pressures obtained from the analytical model were applied. A temperature
change of 1500F was also considered. The response of cement systems 1 and 2 were
Fig 4.32: Equivalent Stress for Cement System 2 under Static loading
Fig 4.33: Alternating Stress for Cement System 2 under a Zero Based Cyclic Loading
81
Fig 4.34: Life Cycle for Cement System 2 under a Zero Based Cyclic Loading
The fatigue sensitivity plots are as shown in Figs. 4.35 to 4.37 below.
82
Fig 4.35: Fatigue Sensitivity to Life for a Zero Based Cyclic Loading with Cement 2
Fig 4.36: Fatigue Sensitivity to Damage plot for a Zero Based Cyclic Loading with Cement
2
83
Fig 4.37: Fatigue Sensitivity to Safety Factor plot for a Zero Based Cyclic Loading with
Cement 2
From the fatigue sensitivity plots, it can be seen that the maximum load ratio that
can be supported by cement system 2 before failure is about 0.75. The same analysis for the
ductile cement (cement system 1), and the results are shown in Figs.4.38 to 4.40 below.
84
Fig 4.38: Equivalent Stress for Cement System 1 under Static Loading
Fig 4.39: Alternating Stress for Cement System 1 under Zero Based Cyclic Loading
85
Fig 4.40: Life Cycle for Cement System 1 under Zero Based Cyclic Loading
For cement system 1, the magnitude of the alternating stress appears to be smaller
than the equivalent stress under static loading (Figs. 4.38 & 4.39).This trend is quite
different from what was obtained with cement system 2 where the minimum alternating
stress under a zero based cyclic loading is greater than the maximum equivalent stress
under static loading (Figs. 4.32 and 4.33).The fatigue sensitivity plots (Figs.4.41 to 4.43)
show that the maximum load ratio that can be supported by cement system 1 before failure
is about 1.1.
86
Fig 4.41: Fatigue Sensitivity to Life for Zero Based Cyclic Loading with Cement 1
Fig 4.42: Fatigue Sensitivity to Safety Factor for Zero Based Cyclic Loading with Cement
1
87
Fig 4.43: Fatigue Sensitivity to Damage for Zero Based Cyclic Loading with Cement 1
88
CHAPTER V
5.1 Conclusions
The drilling and completion of a well is a capital project that needs to be executed
forces that may affect the integrity of a well throughout its life span.
cement design. Designs based solely on static loading conditions may or may not be
enough to ensure long term integrity depending on prevailing downhole conditions thus
the need to take the analysis further by also examining the effect of fatigue.
Life prediction is possible through models that define the operating stress
principles in a well. Such are incorporated into the analytical model presented in this
thesis, which can be used to evaluate stresses in the cement sheath based on actual
wellbore parameters and in synergy with finite element models, can be combined
effectively to evaluate the fatigue and static loading behavior of the well thereby helping
to predict the life of the well. These tools will also help in optimizing the design with
regards to the material properties of both the cement and casing which would also help
Insights were also thrown into some fatigue and static behavior of well cements.
Fatigue failure in cement occurs when microscopic damage within the microstructure of
the cement caused by initial cyclic loading turns into macroscopic cracks under
89
gradually increasing loads. Cyclic loading impacts initial damage and if loading is
continued at load ratios above the critical ratio for a particular cement mix, failure is
imminent but may undergo many cycles when loaded below this ratio. Loading
conditions may affect the fatigue property of cement only when the mechanical
The mechanical properties of cement play a very important role in the static and
fatigue performance of cement. Ductile cement s systems – cements with low Young’s
modulus and a high Poisson ratio generally perform better under static and cyclic
loading conditions as compared to brittle cement systems i.e. cement systems with a high
Young’s modulus and low Poisson ratio. Ductile cement systems generates a
significantly lower value of tangential and radial stresses while brittle cements are more
likely to generate higher tensile and radial stresses within its microstructure under a
The magnitude of confining stress and the mechanical properties of the formation
also play an important role in the static and fatigue behavior of both the cement and
casing. A large far field stress (formation pressure) act to increase the performance of the
casing and counteracts high internal pressures ensuring a minimal transfer to the cement
sheath. Also the more brittle the formation (in terms of Young’s modulus and Poisson
ratio), the more stresses that will be transmitted to the casing and cement sheath.
90
This work includes report on the effect of both static and fatigue behavior of well cement
based on analytical and finite element models. A lot of experimental work is required in
Developing equations specific to well cement. Equations 3.3, 3.4 and 3.26 can be
derived specifically for well cement from experimental data and linear regression
analysis. New failure mechanism, crack initiation and propagation and failure theories
can also be developed from these data. This would help expand the analytical model to
• The findings reported in this work are centered mainly on the mechanical
properties of the cement and on loading conditions. The effect of other factors
• The effect of additives on the static and fatigue properties of well cement.
• Performance of new cement system with special properties like foam and
expansive cements should also be studied and data generated for them
91
REFERENCES
1. Cook, R.D. and Young, W.C.: Advanced Mechanics of Materials, 2nd ed. Prentice Hall,
Upper Saddle River, NJ. 1999.
2. Kim, J. and Kim, Y.: “Experimental Study of the Fatigue Behavior of High Strength
Concrete,” Cement and Concrete Research, Vol. 26, No. 10, 1996. 1513-1523.
3. Antrim, J.D.: “The Mechanism of Fatigue in Cement Paste and Plain Concrete”,
Highway Research Record, 210 no 6, 1965.
4. Breitenbucher, R., Alawieh, H., and Ibuk, H.: “Influence of Cyclic Loading on the
Degradation of Mechanical Concrete Properties,” Advances in Construction Materials,
(August 2007) 317–324.
6. Nayeb, H., Hashemi I., Cohen, M. D., and Erturk, T.: “Evaluation of Fatigue Damage
on the Mechanical Properties of Fiber Reinforced Cement Pastes,” Cement and Concrete
Research,(May 1985)15, 879-888.
7. Bourgoyne, A.T., Millheim, K.K., Chenevert, M.E., and Young, F.S.: Applied Drilling
Engineering, SPE Textbook Series, Vol. 2, 1986.
8. Nelson, E.B.: “Well Cementing,” Development in Petroleum Science, Vol. 28, 1990.
11. Edgley, K., D., Sabins, F., L. and Watters, L., T.: “Supercement for Annular Seal and
Long Term Integrity in Deep, Hot Wells “ ‘Deep Trek’,” Phase II Annual Report, CSI
Technologies, Houston. August 2005.
12. Pine, M., Hunter, L., Mutch, J., Adam, J., and Griffith, E.J.: “Selection of Foamed
Cement for HPHT Gas Wells Proves Effective for Zonal Isolation-Case History,” paper
IADC/SPE 79909, IADC/SPE Drilling Conference, Amsterdam, Netherlands, February
2003.
92
13. Ravi, K., Bosma, M., and Hunter, L.: “Optimizing the Cement Sheath Design in
HPHT Shearwater Field,” paper IADC/SPE 79905, IADC/SPE Drilling Conference,
Amsterdam, Netherlands, February 2003.
14. Ravi, K and Xenakis, H.: “Cementing Process Optimized to Achieve Zonal
Isolation”. Halliburton, January 2007.
16. Godwin, K.J., and Crook R.J.: “Cement Sheath Stress Failure,” SPE Drilling
Engineering, pp. 291-296, December 1992.
17. Fleckenstein, W.W., Eustes, A.W., and Miller, M. G.: “Burst Induced Stresses in
Cemented Well Bores,” paper SPE 62596, SPE/AAPG Western Regional Meeting, Long
Beach, California, June 2000.
18. Gray, K.E., Podnos, E., and Becker, E.: “Finite Element Studies of Near-Wellbore
Region during Cementing Operations: Part 1,” Paper SPE 106998, SPE Production and
Operations Symposium, Oklahoma City, Oklahoma, April 2007.
19. Heathman, J., and Beck, F.E.: “Finite Element Analysis Couples Casing and Cement
Design for HP/HT Wells in East Texas,” paper SPE 98896, IADC/SPE Drilling
Conference, Miami, Florida, February 2006.
20. Rodriguez, W. J., Fleckenstein W.W., and Eustes A.W.: “Simulation of Collapsed
Load on Cemented Casing Using Finite Element Analysis,” paper SPE 84566, SPE
Annual Technical Conference, Denver, Colorado, October 2003.
21. Schubert, J.J., Shahri, M.A., and Amani, M.: “Detecting and Modeling Cement
Failure in High Temperature/High Pressure (HP/HT) Wells Using Finite Element
Methods,” paper IPTC 10961, IPTC Conference, Doha, Qatar, November 2005.
22. Ravi, K., Bosma, M., and Gastebled, O.: “Improve the Economics of Oil and Gas
Wells by Reducing the Risk of Cement Failure,” paper SPE 74497, IADC/SPE Drilling
Conference, Dallas, Texas, February 2002.
23. Krusche, K., Johnson, C.R., Braud, N.Y., and Ghazi, H.B.: “Application of Cement
Engineered Solutions to Solve Long Term Cement Integrity Issues in Tunisia,” paper
SPE 100390, SPE Annual Technical Conference, San Antonio, Texas, September 2006.
93
24. Avram, C., Facaoaru, I., Mirsu, O., Filimon, I., and Tertea, I.: “Concrete Strength and
Strains,” Developments in Civil Engineering, Vol. 3, Elsevier Scientific Publishing
Company, New York, 1981.
25. Newman, J., and Choo, B.S.: “Concrete Properties”, Advanced Concrete Technology,
Elsevier Butterworth-Heinemann, Boston, MA. 2003.
2. Bocca, P. and Crotti, M.: “Variations in the Mechanical Properties and Temperature
of Concrete Subjected to Cyclic Loads, Including High Loads,” Materials and
Structures, (February 2003)36, 40-45.
3. Morris, A. D., and Garrett, G. G.: “A Comparative Study of the Static and Fatigue
Behavior of Plain and Steel Fiber Reinforced Mortar in Compression and Direct
Tension,” The International Journal of Cement Composites and Lightweight Concrete,
(May 1981) 3, no 2, 73-91.
5. Alliche, A., and Francois, D.: “Fatigue Behavior of Hardened Cement Paste,” Cement
and Concrete Research, (January 1986)16, 199-206.
6. LeRoy-Delage, S., Baumgarte, C., Thiercelin M., and Vidick B.: “New Cement
Systems for Durable Zonal Isolation,” paper SPE 59132, IADC/SPE Drilling
Conference, New Orleans, Louisiana, February 2000.
7. O’Brien, B.T.: “Case against Cementing Casing - Casing Annuli,” paper SPE 35106,
IADC/SPE Drilling Conference, New Orleans, Louisiana, March 1996.
10. Thiercelin, M.J., Dargaud, B., Baret, J.F., Rodriguez, W. J.: “Cement Design Based on
Cement Mechanical Response,” SPE Annual Technical Conference, San Antonio, Texas,
337–348, 1997.
94
VITA