Geomechanical Modeling of In-Situ Stresses Around a Borehole

Total Page:16

File Type:pdf, Size:1020Kb

Geomechanical Modeling of In-Situ Stresses Around a Borehole Geomechanical Modeling of In-Situ Stresses Around a Borehole Samantha Grandi, Rama Rao, and M. Nafi Toks¨oz Earth Resources Laboratory Dept. of Earth, Atmospheric, and Planetary Sciences Massachusetts Institute of Technology Cambridge, MA 02139 Abstract In this paper, we present a modelling of the in-situ stress state associated with the severe hole enlargement of a wellbore. Geomechanical information is relevant to assure wellbore stability, i.e., to prevent damages in the formation and later on, the casing. Many of the drilling parameters, as mud weight or the optimal orientation of the borehole, require some knowledge of the mechanical behaviour of the rock. The lack of these kind of data in exploratory areas, where there are usually insufficient constraints for the geological model, increases even more the risk, hence the costs. The present model uses the concepts of poroelasticity theory to compute the stationary 2D, brittle response of the formation around a borehole that is submitted to effective compressive horizontal stresses. The numerical solution is obtained using a finite element approximation. The initial stress state at the far field was estimated combining a frictional-failure theory with the observations of dipmeter caliper in a particular borehole that presents elongations in a preferential di- rection. The direction and relative extension of the observed breakouts at a particular depth are modelled successfully using formation realistic parameters and dimensions, although the exact shape of the borehole (at all angles) was unknown. For the particular case study, the orientation of the breakout is NE-SW, at about 82 degrees azimuth. Therefore, the maximum horizontal stress lies at approximately 350 degrees azimuth. The ratios of horizontal principal stresses to vertical stress that best honor the observations are SHmax = 2.3Sv and Shmin = 1.7Sv. The compressive strength necessary for the rock to fail, as indicated by the caliper data under this stress field, is about 140 MPa. 1 Introduction In this paper, we present a numerical modelling of the geomechanical conditions of a wellbore where an heterogeneous stress field is responsible for numerous breakouts. The governing equations are the equations of motion (Navier’s equations) for perfectly elastic materials, therefore, only a brittle behavior is assumed and no plasticity is included. A numerical solution computed using the finite element method is presented below. The modelling carried out is a mechanical one and thermal stresses are disregarded. The displacements are assumed invariables with time. To have the ability of predicting the geomechanical behaviour of the column rock that is to be drillled is a dream-like scenario for the oil industry. In exploration, it is specially relevant, because the uncertainties of the area difficult the selection of drilling parameters as mud weigth, liners depths and so forth. The issue of wellbore stability has gained particular importance in the last 10 years as a response to the increasing of exploration of complex areas which represent major engineering challenges in drilling and production (Charlez and Onaisi, 2001). Severe hole enlargement has a long list of consequences, from reducing the drill bit life to causing cleaning hole problems, stuck pipe, poor cementing, logging problems, and often the need to sidetrack (Last, 2001). There are several problems in engineering to which wellbore stability refers. The one addressed in this project, is related to drilling performance where instability is mainly driven by stress concentrations around the cavity. 1 During drawdown, depletion of fluid could alter pressures significantly. Sometimes, and very much associated with the kind of lithology supporting the reservoir, the casing could suffer important damage and even the total loss of the wellbore (Frederick et al., 1998). A third problem considered as wellbore instability, is sanding production in the completion stage. During drilling, there are other sources which can also contribute to the development of stability problems. In particular, the orientation of the well has great relevance. Selecting a drilling orientation that does not agree with the natural weak orientation of the formation could introduce difficulties that translate in time and money. On the other hand, there are some laboratory experiments that indicate the damaging effect of some mud components reactions (Charlez and Onaisi, 2001), however, these could only be responsible for washouts of uniform-radius occurrence and not for breakout areas with preferential direction. There are well established methods to derive the in-situ stress state in a borehole. Among them, breakout analysis, hydraulic fracturing experiments and earthquake focal plane mechanisms, are of wide used (Zoback, 1989). Barton and Zoback (1988) described another stress measurement technique based on stress-induced polarization of tube waves (Stoneley). Sinha et al. (1994) introduced a method that uses the flexural mode recorded in cross-dipole data to distinguish induced stresses in the borehole. Finally, the resolution improvement of borehole imaging logs, have made the in-situ stress derivation another application of this tool, overcoming the limitation of four-arms caliper data to provide all-azimuth measurements (Brudy and Kjorholt, 2000). Each method has its own limitations. In all of these techniques (with the exception of focal mechanisms), the stress state inverted includes the effects of any drill-induced stress in the wellbore superimposed to the tectonic field. If the formation is anisotropic, the stress anisotropy would amplify or alter by the presence of the bore (Zoback, 1989), being these effects sometimes larger than the in-situ stresses existing before drilling the well. Knowledge of the ambient stress field is important in elucidating tectonic processes. To understand the behaviour of a wellbore in a tectonically pre-stressed region, forward modelling is a necessary procedure. In this project we are willing to quantify the tectonic stress effect in the formation around a well that does not present any significant anisotropy, at least at the depth of study. However, the heterogeneties of the far field stresses are responsible for important geomechanical damages. The paper is structured in the following way: first, a background theory is presented to facilitate the description of the numerical modelling. Second, the problem is described along with the derivation of the model parameters. Finally, the results are presented and discussed. In brief, the formation parameters are computed from the well data available and the initial range of values of the principal stresses is obtained by a standard analysis of breakouts inferred from the dipmeter caliper data. 2 Theoretical Background 2.1 Structural Mechanics When a layer of rock, buried at certain depth, is submitted to a stress field, the displacements and displace- ment gradients are sufficiently small for the theory of linear elasticity to be valid. The equations of motion, or equilibrium, can be considered to be satisfied in the undeformed reference configuration. We chose the Lagrangian description of motion to express the global conservation laws. Because of the fixed material viewpoint, the conservation of mass is automatically satisfied. Using conservation of momentum, the Cauchy’s equations of motion are obtained (Malvern, 1969). In vectorial notation, ∂u ∇T + f = ρ , (1) ∂t where u is the velocity, ρ is the density of the material, f represents the body forces per unit volume and T is the stress tensor defined as: σx τxy τzx T = τxy σy τyz (2) τzx τyz σz σi represent the normal stresses and τij the shear stresses. 2 ∂u In the special case of static equilibrium, the acceleration ∂t is zero and equation (1) reduces to: ∇T + f = 0 (3) In general, Tij = Tji, and there are only six independent unknown stress components instead of nine. However, the 3 equations (1) are still not enough to solve for the stresses, and constitutive equations have to be used. The stress-strain relations describe the ideal macroscopic behavior of the material under consideration. For an ideal elastic solid, the constitutive equations are known as Hooke’s law : T = C, (4) where C is a fourth-rank tensor with 36 independent coefficients (in the most general case), and is the strain tensor x γxy γzx = γxy y γyz (5) γzx γyz z By definition, the strains are related to the displacements,u ¯ = (u,v,w), as: γ 1 ∂u ∂v = ∂u ; = xy = + x ∂x xy 2 2 ∂y ∂x γ 1 ∂v ∂w = ∂v ; = yz = + (6) y ∂y yz 2 2 ∂z ∂y γ 1 ∂w ∂u = ∂w ; = zx = + z ∂z zx 2 2 ∂x ∂z For an elastically isotropic material, the coefficients of C can be expressed in terms of only two elastic modulii. Usually, these are either the Lame constants, λ and µ, the Young’s modulus (E) and Poisson ratio ν, or the bulk modulus K and the shear modulus (G). The bulk modulus is a measure of the fractional volume change when the pressure is increased by δP . It is obtained through isotropic compression experiments. Most frequently, uniaxial compression is simpler to perform, in which case, it is the Young’s modulus that it is measured. By definition, σ E = zz (7) zz On the other hand, the deformations that take place in the transversal direction are described by the Poisson ratio: ν = − yy (8) zz (Gueguen, 1994). Using the Young’s modulus (E), and the Poisson’s ration (ν), Hooke’s law reduces to: E ν σx = 1+ν x + 1−2ν (x + y + z) ; τxy = Gγxy E ν σy = 1+ν y + 1−2ν (x + y + z) ; τyz = Gγyz (9) E ν σz = 1+ν z + 1−2ν (x + y + z) ; τzx = Gγzx where G is defined as: E G = (10) 2(1 + ν) Finally, substituting equation (6) in (9) and equation (9) in (3), Navier’s displacement equations of motion are: −∇C∇u = f, (11) 3 where, 2G + µ 0 0 0 µ 0 0 0 µ 0 G 0 G 0 0 0 0 0 0 0 G 0 0 0 G 0 0 0 G 0 G 0 0 0 0 0 C = µ 0 0 0 2G + µ 0 0 0 µ (12) 0 0 0 0 0 G 0 G 0 0 0 G 0 0 0 G 0 0 0 0 0 0 0 G 0 G 0 µ 0 0 0 µ 0 0 0 2G + µ Eν and µ = (1+ν)(1−2ν) Once the derivatives of the displacements have been obtained, all strain and stress components can be derived.
Recommended publications
  • Mapping Thixo-Visco-Elasto-Plastic Behavior
    Manuscript to appear in Rheologica Acta, Special Issue: Eugene C. Bingham paper centennial anniversary Version: February 7, 2017 Mapping Thixo-Visco-Elasto-Plastic Behavior Randy H. Ewoldt Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA Gareth H. McKinley Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA, USA Abstract A century ago, and more than a decade before the term rheology was formally coined, Bingham introduced the concept of plastic flow above a critical stress to describe steady flow curves observed in English china clay dispersion. However, in many complex fluids and soft solids the manifestation of a yield stress is also accompanied by other complex rheological phenomena such as thixotropy and viscoelastic transient responses, both above and below the critical stress. In this perspective article we discuss efforts to map out the different limiting forms of the general rheological response of such materials by considering higher dimensional extensions of the familiar Pipkin map. Based on transient and nonlinear concepts, the maps first help organize the conditions of canonical flow protocols. These conditions can then be normalized with relevant material properties to form dimensionless groups that define a three-dimensional state space to represent the spectrum of Thixotropic Elastoviscoplastic (TEVP) material responses. *Corresponding authors: [email protected], [email protected] 1 Introduction One hundred years ago, Eugene Bingham published his first paper investigating “the laws of plastic flow”, and noted that “the demarcation of viscous flow from plastic flow has not been sharply made” (Bingham 1916). Many still struggle with unambiguously separating such concepts today, and indeed in many real materials the distinction is not black and white but more grey or gradual in its characteristics.
    [Show full text]
  • Using Effective Stress Theory to Characterize the Behaviour of Backfill
    Paper 27 Minefill Using effective stress theory to characterize the behaviour of backfill A.B. Fourie, Australian Centre for Geomechanics, University of Western Australia, Perth, Western Australia, M. Fahey and M. Helinski, University of Western Australia, Perth, Western Australia ABSTRACT The importance of understanding the behaviour of cemented paste backfill (CPB) within the framework of the effective stress concept is highlighted in this paper. Only through understand- ing built on this concept can problems such as barricade loads, the rate of consolidation, and arch- ing be fully understood. The critical importance of the development of stiffness during the hydration process and its impact on factors such as barricade loads is used to illustrate why conven- tional approaches that implicitly ignore the effective stress principle are incapable of capturing the essential components of CPB behaviour. KEYWORDS Cemented backfill, Effective stress, Barricade loads, Consolidation INTRODUCTION closure (i.e. the stiffness or modulus of compressibil- ity) that is of primary interest. As the filled stopes The primary reason for backfilling will vary from compress under essentially confined, one-dimen- site to site, with the main advantages being one or a sional (also known as K0) conditions, the stiffness combination of the following: occupying the mining increases with displacement. Well-graded hydaulic void, limiting the amount of wall convergence, use as fills have been shown to generate the highest settled a replacement roof, minimizing areas of high stress density and this, combined with the exponential con- concentration or providing a self-supporting wall to fined compression behaviour of soil (Fourie, Gür- the void created by extraction of adjacent stopes.
    [Show full text]
  • Landslide Triggering Mechanisms
    kChapter 4 GERALD F. WIECZOREK LANDSLIDE TRIGGERING MECHANISMS 1. INTRODUCTION 2.INTENSE RAINFALL andslides can have several causes, including Storms that produce intense rainfall for periods as L geological, morphological, physical, and hu- short as several hours or have a more moderate in- man (Alexander 1992; Cruden and Vames, Chap. tensity lasting several days have triggered abun- 3 in this report, p. 70), but only one trigger (Varnes dant landslides in many regions, for example, 1978, 26). By definition a trigger is an external California (Figures 4-1, 4-2, and 4-3). Well- stimulus such as intense rainfall, earthquake shak- documented studies that have revealed a close ing, volcanic eruption, storm waves, or rapid stream relationship between rainfall intensity and acti- erosion that causes a near-immediate response in vation of landslides include those from California the form of a landslide by rapidly increasing the (Campbell 1975; Ellen et al. 1988), North stresses or by reducing the strength of slope mate- Carolina (Gryta and Bartholomew 1983; Neary rials. In some cases landslides may occur without an and Swift 1987), Virginia (Kochel 1987; Gryta apparent attributable trigger because of a variety or and Bartholomew 1989; Jacobson et al. 1989), combination of causes, such as chemical or physi- Puerto Rico (Jibson 1989; Simon et al. 1990; cal weathering of materials, that gradually bring the Larsen and Torres Sanchez 1992)., and Hawaii slope to failure. The requisite short time frame of (Wilson et al. 1992; Ellen et al. 1993). cause and effect is the critical element in the iden- These studies show that shallow landslides in tification of a landslide trigger.
    [Show full text]
  • Slope Stability
    Slope stability Causes of instability Mechanics of slopes Analysis of translational slip Analysis of rotational slip Site investigation Remedial measures Soil or rock masses with sloping surfaces, either natural or constructed, are subject to forces associated with gravity and seepage which cause instability. Resistance to failure is derived mainly from a combination of slope geometry and the shear strength of the soil or rock itself. The different types of instability can be characterised by spatial considerations, particle size and speed of movement. One of the simplest methods of classification is that proposed by Varnes in 1978: I. Falls II. Topples III. Slides rotational and translational IV. Lateral spreads V. Flows in Bedrock and in Soils VI. Complex Falls In which the mass in motion travels most of the distance through the air. Falls include: free fall, movement by leaps and bounds, and rolling of fragments of bedrock or soil. Topples Toppling occurs as movement due to forces that cause an over-turning moment about a pivot point below the centre of gravity of the unit. If unchecked it will result in a fall or slide. The potential for toppling can be identified using the graphical construction on a stereonet. The stereonet allows the spatial distribution of discontinuities to be presented alongside the slope surface. On a stereoplot toppling is indicated by a concentration of poles "in front" of the slope's great circle and within ± 30º of the direction of true dip. Lateral Spreads Lateral spreads are disturbed lateral extension movements in a fractured mass. Two subgroups are identified: A.
    [Show full text]
  • Soil Mechanics Lectures Third Year Students
    2016 -2017 Soil Mechanics Lectures Third Year Students Includes: Stresses within the soil, consolidation theory, settlement and degree of consolidation, shear strength of soil, earth pressure on retaining structure.: Soil Mechanics Lectures /Coarse 2-----------------------------2016-2017-------------------------------------------Third year Student 2 Soil Mechanics Lectures /Coarse 2-----------------------------2016-2017-------------------------------------------Third year Student 3 Soil Mechanics Lectures /Coarse 2-----------------------------2016-2017-------------------------------------------Third year Student Stresses within the soil Stresses within the soil: Types of stresses: 1- Geostatic stress: Sub Surface Stresses cause by mass of soil a- Vertical stress = b- Horizontal Stress 1 ∑ ℎ = ͤͅ 1 Note : Geostatic stresses increased lineraly with depth. 2- Stresses due to surface loading : a- Infintly loaded area (filling) b- Point load(concentrated load) c- Circular loaded area. d- Rectangular loaded area. Introduction: At a point within a soil mass, stresses will be developed as a result of the soil lying above the point (Geostatic stress) and by any structure or other loading imposed into that soil mass. 1- stresses due Geostatic soil mass (Geostatic stress) 1 = ℎ , where : is the coefficient of earth pressure at # = ͤ͟ 1 ͤ͟ rest. 4 Soil Mechanics Lectures /Coarse 2-----------------------------2016-2017-------------------------------------------Third year Student EFFECTIVESTRESS CONCEPT: In saturated soils, the normal stress ( σ) at any point within the soil mass is shared by the soil grains and the water held within the pores. The component of the normal stress acting on the soil grains, is called effective stressor intergranular stress, and is generally denoted by σ'. The remainder, the normal stress acting on the pore water, is knows as pore water pressure or neutral stress, and is denoted by u.
    [Show full text]
  • 3. Stresses in the Ground
    3-1 3. STRESSES IN THE GROUND 3.1 STRESSES IN DRY SOIL Let A in Fig. 3.1 (a) represent a small cubical shaped element of soil at a depth z in an extensive uniform soil deposit in which the ground surface is horizontal and which has been formed by the gradual accretion of material on the ground surface. Because the soil deposit is extensive (by comparison with distance z) in the horizontal direction the stresses on element A will be identical with the stresses on an adjacent element at the same depth below the ground surface. This means that there cannot be any shear stresses existing on the vertical or horizontal planes which bound element A. In other words the vertical stress ( σv) and horizontal stress ( σH) are principal stresses. The vertical stress on element A can be determined simply from the mass of the overlying material. If ρd represents the density of the soil, the vertical stress is σv = ρd gz (3.1) The horizontal stress is customarily expressed as a proportion of the vertical stress σH = K' o σv = K' o ρd gz (3.2) where K' o = coefficient of earth pressure at rest in terms of effective stresses (see equation (3.7)). This coefficient contains the words “at rest” since the soil was deposited under conditions of zero horizontal strain. In other words, because of the large lateral extent of the soil deposit, the vertical planes on any soil element A do not experience any lateral movement as the stresses increase as a consequence of the accretion of material on the ground surface.
    [Show full text]
  • Critical Borehole Orientations – Rock Mechanics Aspects
    DRILLING Critical Borehole Orientations – Rock Mechanics Aspects By R. BRAUN* Abstract again and again during the development of 3 Significant “Rock Mechanics” This article discusses rock mechanics as- the inclination, and these can often only be Factors pects of the relationship between borehole controlled to a limited extent by classical 3.1 In situ stresses stability and borehole orientation. Two countermeasures (variation of mud pressure The magnitudes and orientations of in situ kinds of instability are considered. One is and its composition). It is therefore worth- stresses have a decisive influence on the failures resulting from normal stresses in while in every case to take a more careful borehole loadings. In order to make a rea- planes perpendicular to the borehole, or look at the rock which is to be penetrated and sonable assessment of these in situ loadings along its axis. These may be tension or shear at the possible interactions of this with the it is essential to distinguish between the vari- failures in the direction of the minimum borehole. These interactions may be geolog- ous different components. and/or the maximum normal stress compo- ical/mineralogical or geochemical, but the First the external loading (also called the to- nent related to the borehole, or peeling/loos- ones discussed here concern rock mechan- tal loading) is to be considered. The vertical ening around the complete borehole perime- ics. The significant rock mechanics factors component of this is the total overburden ter. These can be avoided, or at least mini- and the way in which they can influence pressure and it can therefore be determined mized, by the selection of an appropriate borehole stability are discussed in the directly (for example from the average den- borehole pressure.
    [Show full text]
  • Consideration on Landslide Mechanism Based on Pore Water Pressure Loading Test Examan Du Mécanisme De Glissement De Terrains Par Essai De Suppression Hydrostatique
    INTERNATIONAL SOCIETY FOR SOIL MECHANICS AND GEOTECHNICAL ENGINEERING This paper was downloaded from the Online Library of the International Society for Soil Mechanics and Geotechnical Engineering (ISSMGE). The library is available here: https://www.issmge.org/publications/online-library This is an open-access database that archives thousands of papers published under the Auspices of the ISSMGE and maintained by the Innovation and Development Committee of ISSMGE. Consideration on landslide mechanism based on pore water pressure loading test Examan du mécanisme de glissement de terrains par essai de suppression hydrostatique S-Ohtsuka & YMiyata — Nagaoka University of Technology, Nagaoka, Niigata, Japan ABSTRACT: Pore water pressure loading test of clay was conducted under the constant deviator stress. From the test results, the yield and failure stresses were newly defined and the yield stress was found located on the critical state line. With the increase in pore water pressure over the yield stress, the shear deformation developed with long time. The progressive failure with water migration was pointed out as the reason for the slow movement of landslide. Numerical simulation on shear strength of clay against repetitive pore water pressure loading clarified the process of deterioration in shear strength of clay caused by repeats of landslide movement. RESUME: L’essai du chargement de la pression de l’eau interstitielle de argile a été procédé sous la contrainte constante de cisaillement. Selon les résultats, la contrainte d’écoulement et celle de rupture ont été définies et la première a été trouvée aux environs de la ligne de l’état critique. En augmentant la pression de l’eau interstitielle jusqu’au-dessus de la contrainte d’écoulement, la déformation de cisaillement s’est formée et s’est développée en longues heures.
    [Show full text]
  • Lecture 1 : Macroscopic Characterisation of Soils
    Lecture 1 : Macroscopic Characterisation of Soils Kenichi Soga, University of Cambridge 1.1 Introduction All aspects of soil stability—bearing capacity, slope stability, the supporting capacity of deep foundations, and penetration resistance, to name a few—depend on soil strength. The stress-deformation and stress-deformation-time behaviour of soils are important in any problem where ground movements are of interest. Most relationships for the characterization of the stress-deformation and strength properties of soils are empirical and based on phenomenological descriptions of soil behaviour. The Mohr–Coulomb equation is by far the most widely used for strength. It states that τ ff = c + σ ff tanφ (1.1) τ ff = c′ + σ ′ff tanφ′ (1.2) where τff is shear stress at failure on the failure plane, c is a cohesion intercept, σff is the normal stress on the failure plane, and φ is a friction angle. Equation (1.1) applies for σff defined as a total stress, and c and φ are referred to as total stress parameters. Equation (1.2) applies for σ′ff defined as an effective stress, and c′ and φ′ are effective stress parameters. As the shear resistance of soil originates mainly from actions at interparticle contacts, the second equation is the more fundamental. In reality, the shearing resistance of a soil depends on many factors, and a complete equation might be of the form Shearing resistance = F(e,c',φ',σ ',C, H,T,ε,ε, S) (1.3) in which e is the void ratio, C is the composition, H is the stress history, T is the temperature, ε is the strain, ε is the strain rate, and S is the structure.
    [Show full text]
  • Cyclic Response of a Sand with Thixotropic Pore Fluid
    Cyclic Response of a Sand with Thixotropic Pore Fluid C.S. El Mohtar1, J. Clarke1, A. Bobet1, M. Santagata1, V. Drnevich1 and C. Johnston2 1School of Civil Engineering, Purdue University, West Lafayette, IN 2Department of Agronomy, Purdue University, West Lafayette, IN ABSTRACT Saturated specimens of Ottawa sand prepared with 0%, 3% and 5% bentonite by dry mass of sand are tested under cyclic loading to investigate the effects of bentonite on the cyclic response. For the same skeleton relative density and cyclic stress ratio (CSR), the cyclic tests on the sand-bentonite mixtures show a significant increase of the number of cycles required for liquefaction compared to the clean sand. This is caused , as observed in resonant column tests, by an increase of the elastic threshold due to the presence of bentonite, which delays the generation of excess pore pressure. Such behavior can be explained by the rheological properties of the pore fluid. Oscillatory tests conducted with a rheometer on bentonite slurries show that for shear strains as large as 1% these materials exhibit elastic behavior with a constant shear modulus. Moreover, due to the thixotropic nature of the bentonite slurries, their storage modulus shows a marked increase with time. This observation is consistent with the increase in the liquefaction resistance of the sand-bentonite mixtures with time also observed in cyclic triaxial experiments. INTRODUCTION Liquefaction is an important cause of damage to civil infrastructures during earthquakes. Notable examples are: the collapse
    [Show full text]
  • Tensile Strength, Shear Strength, and Effective Stress for Unsaturated Sand
    TENSILE STRENGTH, SHEAR STRENGTH, AND EFFECTIVE STRESS FOR UNSATURATED SAND A Dissertation Presented to the Faculty of the Graduate School University of Missouri – Columbia In Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy By RAFAEL BALTODANO GOULDING Dr. William J. Likos, Dissertation Supervisor MAY 2006 Abstract It is generally accepted in geotechnical engineering that non-cohesive materials such as sands exhibit no or negligible tensile strength. However, there is significant evidence that interparticle forces arising from capillary and other pore-scale force mechanisms increase both the shear and tensile strength of soils. The general behavior of these pore-scale forces, their role in macroscopic stress, strength, and deformation behavior, and the changes that occur in the field under natural or imposed changes in water content remain largely uncertain. The primary objective of this research was to experimentally examine the manifestation of capillary-induced interparticle forces in partially saturated sands to macroscopic shear strength, tensile strength, and deformation behavior. This was accomplished by conducting a large suite of direct shear and direct tension tests using three gradations of Ottawa sand prepared to relatively “loose” and relatively “dense” conditions over a range of degrees saturation. Results were compared with previous experimental results from similar tests, existing theoretical formulations to define effective stress in unsaturated soil, and a hypothesis proposed to define a direct relationship between tensile strength and effective stress. The major conclusions obtained from this research include: Theoretical models tended to underpredict measured tensile strength. Analysis of results indicates that shear strength may be reasonably predicted using the sum of tensile strength and total normal stress as an equivalent effective stress (σ’ = σt + σn).
    [Show full text]
  • 11. the Stability of Slopes
    11-1 11. THE STABILITY OF SLOPES 11.1 INTRODUCTION The quantitative determination of the stability of slopes is necessary in a number of engineering activities, such as: (a) the design of earth dams and embankments, (b) the analysis of stability of natural slopes, (c) analysis of the stability of excavated slopes, (d) analysis of deepseated failure of foundations and retaining walls. Quite a number of techniques are available for these analyses and in this chapter the more widely used techniques are discussed. Extensive reviews of stability analyses have been provided by Chowdhury (1978) and by Schuster and Krizek (1978). In order to provide some basic understanding of the nature of the calculations involved in slope stability analyses the case of stability of an infinitely long slope is initially introduced. 11.2 FACTORS OF SAFETY The factor of safety is commonly thought of as the ratio of the maximum load or stress that a soil can sustain to the actual load or stress that is applied. Referring to Fig. 11.1 the factor of safety F, with respect to strength, may be expressed as follows: τ F = ff (11.1) τ where τ ff is the maximum shear stress that the soil can sustain at the value of normal stress of σn, τ is the actual shear stress applied to the soil. Equation 11.1 may be expressed in a slightly different form as follows: c σ tan φ = n τ F + F (11.2) Two other factors of safety which are occasionally used are the factor of safety with respect to cohesion, F c, and the factor of safety with respect to friction, F φ.
    [Show full text]