Subsections

The lithostatic stress gradient is the variation of total vertical stress with vertical depth (usually referred as true depth in petroleum engineering). The following subsections review the fundamental concepts of stress, stress equilibrium, and effective stress.

## 2.1.1 Stress concept and equilibrium 1D

Consider the solid bar in Figure 2.1. Two forces and pull the bar in opposite directions. Underlining indicates the variable is a vector in (three dimensions). Equilibrium requires summation of (vectorial) forces to be zero, hence, . If the solid bar has a weight, then equilibrium requires , the difference is the weight .

Stress is defined as force over area (perpendicular to the force) such that

 (2.1)

The units of stress are [Force]/[Area]: MPa, psi, etc., but it does not necessarily mean stress is a pressure! Stress depends on the direction in which is measured, while pressure is the same in all directions. The directionality of stress is a result of the solid capacity to resist shear stresses.

## 2.1.2 Rock vertical stress gradient

In a 3D porous solid with volume (Figure 2.2), equilibrium requires the summation of (vectorial) forces in all directions to be zero . Equilibrium of forces in the vertical direction (gravity direction) requires (forces in vertical direction). Hence,

 (2.2) (2.3) (2.4)

Considering infinitesimal variations yields the following equation

 (2.5)

The term within the integral in the right-hand-side is called the vertical total stress gradient, or sometimes, simply as lithostatic stress gradient. In a semi-infinite medium (e.g. approximation of the Earth's surface) .

If , then vertical stress as a function of depth is

 (2.6)

The total vertical stress is a compressive stress and by convention in geomechanics we assign it a positive sign.

EXAMPLE 2.1: Assume a rock made of 100% quartz (mass density ) with 20% porosity filled with water (mass density ). What is the lithostatic stress gradient?

SOLUTION
The bulk (volume average) rock mass density depends on porosity , volume fractions of mineral phases, and volume fractions of fluid phases. For a water-saturated rock:

Let's assume typical values kg/m and kg/m (check out https://nature.berkeley.edu/classes/eps2/wisc/glossary2.html, look for specific gravity SG). Then,

 kg/m    kg/m    kg/m

This bulk mass density is also equal to 144.8 lb/ft and 19.4 ppg (pounds per gallon), where 1 kg = 2.2 lb, 1 ft = 0.305 m, and 1 ft = 7.48 US gallon. The stress gradient results

 kg/m    m/s    Pa/m    MPa/km

Note that 1 N = 1 kg m/s, 1 Pa = 1 N/m, 1 MPa = 10 Pa and 1 Km = 10 m.

Let us calculate the stress gradient in petroleum field units (psi/ft). First, we need to convert kg to lb (1 kg = 2.205 lb) and m to ft (1 m = 3.281 ft), thus     kg/m    lb/ft. Second, we will pass pounds mass (lb) to pounds force (lbf). This takes into account multiplying for , so you should not multiply any numerical value going from lb to lbf. Third, we will separate the denominator in area times length, and convert ft to in:

 lbf/ft    lbf/ft    1/ft    lbf/(12 in)    1/ft    psi/ft

Typical vertical stress gradient are around psi/ft MPa/km for porosity . Hydrostatic pore pressure gradient is psi/ft MPa/km. You may use fluid saturations if the rock has two or more fluids in the pore space to accurately calculate the vertical stress gradient. You may also use the corresponding mineral volume fractions if the rock is comprised of two or more minerals, e.g., a dolomite-rich shale.

Now that we have the stress gradient, we can calculate total stress at a given depth. Let's consider two cases: onshore and offshore.

### 2.1.2.1 Onshore total vertical stress

First, consider an onshore scenario in which the surface coincides with the water table (in most practical applications this is a reasonable assumption). Assume a constant fluid density and constant bulk rock mass density (remember that it includes fluids within pores). Then, the pore pressure gradient as a function of depth will be

 (2.7)

as long as there is a connected pore network path from the surface to the depth , such that the fluid is in hydrostatic equilibrium (it may not happen sometimes). The total vertical stress gradient as a function of depth is calculated with Equation 2.6 because in onshore conditions the assumption of constant bulk mass density with depth is acceptable. Both (hydrostatic) and increase linearly with depth with constant mass densities (Figure 2.3).

The difference between total stress and pore pressure is called “effective stress”; hence, the effective vertical stress is

 (2.8)

Effective stresses in the subsurface are mostly compressive. Compressive effective stress holds rock together and compacts them. Notice that must hold in order to have a compressive effective vertical stress. We will see later that effective stress is a very important quantity and dictates rock deformation and failure. Figure 2.4 shows an example of effective stress making cohesionless ground coffee strong and stiff like a brick. You will see this example in class how I stand on a “brick” of ground coffee thanks to Terzaghi's effective stress and how effective stress strengthens the coffee pack due to friction forces.

Similar to the membrane effect conceptualized in Figure 2.4, the mudcake that forms around a wellbore during drilling provides an effective stress to the surrounding rock. This is a result of the sharp pressure gradient between the mud pressure in the well and the pore pressure in the formation (Figure 6.3).

EXAMPLE 2.2: Compute the pore pressure and vertical stress at 4000 ft of depth (True Vertical Depth: TVD) underneath an onshore drilling rig. The pore pressure gradient is hydrostatic with brine mass density of 1.04 g/cm and the average rock bulk mass density is about 2.35 g/cm. Make a plot of pressure and vertical stress versus depth. Show results in MPa and psi.

SOLUTION
The corresponding hydrostatic and lithostatic gradients are (keep in mind that 1 g/cm = 1 g/cc = 1000 kg/m):

 kg/m    m/s    MPa/km

and

 kg/m    m/s    MPa/km

Hence, the pore pressure and hydrostatic gradient at 4000 ft (1219.2 m) of depth are

 MPa/km    km    MPa    psi

and

 MPa/km    km    MPa    psi

Note: 1 MPa = 145 psi.

### 2.1.2.2 Offshore total vertical stress

In an offshore case we have to take into account the presence of overlying water weight proportional to thickness . Again, let us assume constant mass densities for simplicity. In the water, you cannot define a stress because there is no solid phase. Within the water zone, fluid pressure will increase according to Eq. 2.7. This fluid pressure will continue increasing within the seafloor sediments as long as there is a connected pore network from the seafloor to depth . Stresses start to develop at the seafloor beyond depth (Figure 2.6), such that the total vertical stress will be

 for (2.9)

Notice that the total stress does not start at 0 at the seafloor . However, the important and “physically meaningful” quantity is effective stress and this one does start at zero at the seafloor.

 for (2.10)

This makes sense because otherwise an octopus would not be able to dig in sand at the seafloor (take a look at this video https://www.youtube.com/watch?v=LSeYVTo7DFs). We will see later that stress is a “tensor” (with magnitude that depends on the orientation of the plane in which it is measured) and that pore pressure affects effective stresses in the normal directions.

EXAMPLE 2.3: Compute the pore pressure and vertical stress at 9000 ft of depth (TVD) underneath an offshore drilling rig with water depth of 2000 ft. The pore pressure gradient is 0.44 psi/ft and the lithostatic gradient is 1 psi/ft. The hydrostatic pressure gradient above the seafloor is 0.44 psi/ft as well. Make a plot of pressure and vertical stress versus depth. Calculate effective stress at 9000 ft of depth. Show results in MPa and psi.

SOLUTION
The pore pressure and hydrostatic gradient at = 9000 ft of total depth with water depth = 2000 ft are

 psi/ft    ft    psi/ft    ft    psi

and

 psi/ft    ft    psi/ft    ft    psi

hence,

 psi    psi    psi

The values in SI units are     MPa,     MPa, and     MPa.

For an interactive example of calculation of vertical stress with depth, check my Jupyter notebook at https://mybinder.org/v2/gh/dnicolasespinoza/GeomechanicsJupyter/master?filepath=PorePressureVerticalStress_Widget.ipynb

### 2.1.2.3 General calculation of total vertical stress

In the general case rock bulk mass density varies with depth. For example, rock lithology will vary with depth. Rocks usually have lower porosity as depth increases. Brine mass density also changes with depth. In general, vertical stress is calculated from integration of the equation

 (2.11)

where is obtained from (gamma ray) density well logs and is obtained from the true vertical depth (TVD) calculated from the well measured depth (MD) and the well deviation survey (Figure 2.7). Accurate calculations should account for the difference between the rotary table or kelly bushing (from where measured depth is often obtained) and the actual ground level. This difference can be significant in offshore cases. In deviated wellbores, you should also take into account wellbore deviation and compute TVD from MD and well trajectory.

Well logs contain “discrete” data, hence you have to do numerical integration. Equation 2.11 can be approximated with

 (2.12)

where is the i-th interval at depth . You can code this equation as the “trapezoidal rule" in a for loop or in a spreadsheet. For example, the vertical stress at the -th layer is

 (2.13)

where is the average height of the trapezoid, is the base of the trapezoid, and is the total stress at the previous depth data point calculated with the same rule, but one level above. The result is the addition of the weight of all layers above the point under consideration.

EXAMPLE 2.4: Go to https://github.com/dnicolasespinoza/GeomechanicsJupyter and download the files 1_14-1_Composite.las and 1_14-1_deviation_mod.dev. The first one is a well logging file (.LAS). You will find here measured depth (DEPTH - Track 1) and bulk mass density (RHOB - Track 8). Track 3 also shows bulk density correction (DRHO). Add RHOB to DRHO to obtain the corrected bulk mass density. The second file has the deviation survey of the well. Use this file to calculate true vertical depth subsea (TVDSS) as a function of measured depth (MD) in the well logging file. You may assume an average bulk mass density of 2 g/cc between the seafloor and the beginning of the density data. Apply Equation 2.12. Summations with discrete data sets can be easily done through a for loop.

SOLUTION
See Figure 2.8. You may use the library https://pypi.org/project/lasio/ to open LAS files with python or write your own code to read the text files .LAS and .DEV.