Chapter 2

Physics of Acoustic Wave

Propagation
To expedite the data processing and simplify the theory explorationists
often invoke the acoustic approximation, i.e., shear effects in the data
are negligible. This is an acceptable approximation for somewhat layered media, near-offset traces recorded by vertical component phones,
and surface-wave filtered data. To deepen our understanding of this
acoustic approximation we now present an overview of the physics of
the acoustic wave equation.

2.1

Acoustic Media and Acoustic Waves

Assume a compressible, non-viscous (i.e., no attenuation) fluid with no

shear strength and in equilibrium (i.e., no inertial forces). Small localized displacements of the fluid will propagate as an acoustic wave, also
known as a compressional wave. Due to the lack of shear strength, localized deformations of the medium do not result in shear deformations
but instead cause changes in the fluid elements volume, as shown in
Figure 2.1.
The equilibrium force/unit area on the face of a volume element
will be called the time-independent equilibrium pressure Peq (r), while
the change in pressure due to a localized compressional wave will be
denoted as P (r, t). For example, the atmospheric pressure decreases
1

2CHAPTER 2. PHYSICS OF ACOUSTIC WAVE PROPAGATION

Equilibrium

Rarefraction

Compression

Figure 2.1: Cube of air in (top) equilibrium and (bottom) disturbed

from equilibrium. The lower left diagram depicts a rarefraction where
the surrounding medium pulls the cube into a larger volume so that
the air density is less than the surrounding medium. The lower right
diagram is similar, except the surrounding medium compresses the cube
into a smaller volume resulting in denser air.

with elevation increases and can be considered to be independent of

time. If I begin talking, however, I excite transient acoustic waves
P (r, t) that disturb the equilibrium pressure.
Snapshots of the particle distribution for a condensation wave and
a rarefraction wave are shown in Figure 2.1. In the compressional
case the element volume is filled with denser (shaded) air while in the
rarefraction case the element volume has lighter (unshaded) air than
the surrounding medium. We can physically create the condensation
wave by injecting air from our lungs into the medium (HELLLL!), and
the rarefraction wave by sucking air into our lungs (UHHHHH!). Using
a spring-mass model, rarefractions are created by pulling a spring (i.e.,
tension) and condensations by compressing a spring (compressions).

2.2. ACOUSTIC HOOKES LAW: P = U

2.2

Acoustic Hookes Law: P = u

Hookes law for an acoustic medium says stress is linearly proportional

to strain for small enough strains. A simple 1-D example will first
be given to demonstrate Hookes law, and then we will apply it to the
case of an acoustic medium.
1. 1-D Spring: The force on a mass connected to a spring disturbed
from equilibrium (see Figure 2.2) is given by

F = k(du/l)k,

(2.1)

where du is the displacement from equilibrium, k is the spring

constant, l is the length of the spring in equilibrium, and k is
the downward point unit vector in Figure 2.2. The ratio du/l is
also known as the compressional strain of the spring model. Note
that in the equilibrium position there is no motion because the
gravitational force balances the elastic force. When disturbed,
the net force on the mass is given by the above equation.
2. 3-D Acoustic Springy Cube: Hookes Law says that the surrounding disturbed medium exerts a pressure on the face of the
cube that is linearly proportional to the corresponding volume
change of the cube. That is,
P = dV /V,

(2.2)

where dV is the change in volume after disturbance, V is the

undisturbed cube volume, and is the bulk modulus. Note that a
volume change is the 3-D equivalent of a 1-D displacement change.
Figure 2.2 shows that the relative volume change dV /V is given
by:
dV
V

= [dxdydz (dx + u)(dy + v)(d

udydz + vdxdz + wdxdy

+ O(2nd order terms)
dxdydz

u v w
+
+
,
x y
z

4CHAPTER 2. PHYSICS OF ACOUSTIC WAVE PROPAGATION

1-D Spring

du

disturbed state

3-D Acoustic Cube

dv

dy

du

disturbance

dz
dx

equilibrium

du dydz + dv dxdz + dw dxdy

dx dy dz
du
dx

dV =

equilibrium

+ dv + dw
dy

dz

du dv dw

dw {

disturbed state

Figure 2.2: (Top) Spring in (left) equilibrium and (right) disturbed from
equilibrium. (Bottom) Elemental cube of air in (left) equilibrium and
(right) disturbed from equilibrium. In this case the cube has expanded
so net tensional forces of surrounding medium must be expanding cube.
Note that the 1st-order volume change is outlined by the solid heavy
lines, where the higher-order terms are associated with the corner parts.

2.3. NEWTONS LAW: P/X = (X, Y, Z) 2 U/T 2

where V = dxdydz is the cubes volume before the disturbance.

Substituting equation 2.3 into equation 2.2 yields:
P = (x, y, z) u + S(x, y, z, t),

(2.4)

where u = (u, v, w) are the cartesian components of displacement

along the x, y, and z axes and S(x, y, z, t) is a time dependent
source term that is independent of the displacement field. For
example, S(x, y, z, t) can be a source that injects material into
the medium such as an air gun used for marine seismic surveys.
Often the spatial and temporal variables are suppressed for the
P and u field variables, but the spatial coordinates are explicitly
expressed for the bulk modulus to remind the reader that the
physical properties of the medium can vary with location.
Note:
Neglecting second-order terms is the small displacement approximation, valid for V /V < 104 or sound quieter than a jet engine
(Kinsler and Frey, 1961, Fundamentals of Acoustics).
Sign convention: P is the force/area that the surrounding medium
exerts on the face of the elemental cube, where tensions (force
pointing away from cube into external medium) are negative and
compressions (force pointing into cube) are positive. For example,
the divergence of u is positive if the volume expands by tensional
forces, which is consistent with the sign in equation 2.4.

2.3

Newtons Law: P/x = (x, y, z) 2u/t2

The external force on an acoustic cube is given in Figure 2.3. These

forces have a non-zero gradient along the x-axis, and so there is a net
elastic force imposed upon the cube by the external medium. This net
force must be balanced by an inertial force (i.e., it accelerates) so that
Newtons law says:
[P (x + dx, y, z, t) P (x, y, z, t)]dydz [(x, y, z)dxdydz]
u(x, y, z, t),
(2.5)

6CHAPTER 2. PHYSICS OF ACOUSTIC WAVE PROPAGATION

TENSIONS ARE PULLS (-)

COMPRESSIONS ARE PUSHES (+)

..

..

P(x, y, z, t)

P(x+dx, y, z, t)

Figure 2.3: External (top) tension and (bottom) compressional forces

on an elemental cube. In the lower case there exists a spatial gradient
of the disturbed pressure field along the x-axis.
where the double dot corresponds to two time derivatives. Expanding
the LHS in a Taylor series about the point (x, y, z) we get
[P (x, y, z, t) P (x, y, z, t) + P/xdx + High Order T erms]dydz
[dxdydz]
u(x, y, z, t), (2.6)
and dividing by dxdydz and neglecting higher-order terms we get Newtons Law:
P (x, y, z, t)/x = (x, y, z)
u(x, y, z, t).

(2.7)

For an arbitrary force distribution the general form of Newtons law is:
P = (x, y, z)
u(x, y, z, t),
(x, y, z, t) is the particle acceleration vector.
where u
Note:

(2.8)

2.4. ACOUSTIC WAVE EQUATION

The minus sign is used so that we are consistent with the notation
for pressure. If P (x + dx, y, z, t) is positive and greater than
P (x, y, z, t) in Figure 2.3 then the cube should accelerate to the
left, which it will according to the above form of Newtons Law.

2.4

Acoustic Wave Equation

Applying 2 /t2 to equation 2.4 and applying to equation 2.8 after

dividing by /(x, y, z) gives:
y, z, t),
+ S(x,
P = (x, y, z) u

(2.9)

.
[1/(x, y, z)P ] = u

(2.10)

The above two equations are the first-order equations of motion in an

acoustic medium. Substituting equation 2.9 into equation 2.10 yields
the 2nd-order acoustic wave equation:
y, z, t)/(x, y, z).(2.11)
(1/(x, y, z)P ) 1/(x, y, z)P = S(x,
This equation is valid for arbitrary velocity and density distributions.
Assuming negligible density gradients the above equation reduces
to the oft-used wave equation:
y, z, t)/(x, y, z),
2 P c2 2 P/t2 = (x, y, z)S(x,

(2.12)

where
c=

(x, y, z)/(x, y, z) ,

(2.13)

and c is the compressional wave propagation velocity.

The above equation is known as the inhomogeneous acoustic wave
equation for negligible density variations, expressed as
2 P c2 2 P/t2 = F ,

(2.14)

y, z, t)/(x, y, z) is the inhomogeneous source term

where F = S(x,
that specifies the location and time history of the source. For example,
F specifies the strength of the hammer blow and its location in the
Liberty Park experiment.

8CHAPTER 2. PHYSICS OF ACOUSTIC WAVE PROPAGATION

2.5

Solutions of the Wave Equation

The physics of wave propagation will now be explained using some special solutions to the wave equation. A harmonic plane wave propagating
in a homogeneous medium will be first examined, and then the case of
a 2-layered medium will be studied.

2.5.1

1-D Wave Propagation in an Homogeneous

Medium

A harmonic wave oscillates with period T and has a temporal dependence usually given by eit , where = 2/T is the angular frequency
inversely proportional to the period T . A plane wave is one in which
the wavefronts line up along straight lines, and can be described by the
following function:
P (x, t) = A0 ei(kxt)

(2.15)

which also solves the homogeneous wave equation 2.12. This can be
shown by plugging equation 2.15 into equation 2.12 to get
(k 2 (/c)2 )P (x, z, t) = 0,

(2.16)

which admits non-trivial solutions if

k 2 = (/c)2 .

(2.17)

This equation is known as the dispersion equation and imposes a constraint on the temporal and spatial variables in the Fourier domain.
It says that the wavefront will move a distance of one wavelength
during one period T of elapsedqtime, and the propagation velocity of
this movement is equal to c = /. See Figure 2.40 for an illusration
of a plane wave propagating plotted as Offset vs Time.
The real part of equation 2.15, i.e., cos(kxt), plots out as slanted
lines in x t space, and these lines move to the right as t increases. In
an x y z volume, these lines (or wavefronts) become planes that are
perpendicular to the x-axis. Therefore this function represents a plane
wave propagating along the x-axis. The shortest distance between 2

2.5. SOLUTIONS OF THE WAVE EQUATION

Wavelength
0

0.02

0.5

0.04
Period

Time (s)

0.06

0.5
0.08

1
0.2

0.1
0.15

0.2
0.15

0.1
0.1
0.05
Time (s)

0.12

0.05
0

Offset (m)

0.05

0.1
Offset (m)

0.15

Figure 2.4: 2-D wavefront propagating along the Offset coordinate.

adjacent peaks of the wavefront is defined to be the wavelength

and its reciprocal is given by k = 2/, where k is known as the
wavenumber. Using this definition of k and that for = 2/T the
dispersion relation can be re-expressed as:

c = /T.

(2.18)

The constant
q c is also defined in terms of the material properties of the
rock as c = /, otherwise known as the compressional wave velocity.
Note: The function ei(kxt) describes a right-going plane wave
and ei(kx+t) describes a left-going wave. To see this, note that we
follow a wavefront if the phase = (kx t) stays constant for an
increase in both t and simultaneous increase in x (such that x/t = /k).
Because x increases this means that the wave is moving to right with the
compressional velocity given in equation 2.18. Conversely, = (kx+t)
is a constant if t increases and x decreases; thus, the wave is moving to
left.

10CHAPTER 2. PHYSICS OF ACOUSTIC WAVE PROPAGATION

2.5.2

2-D Wave Propagation in an Homogeneous

Medium

The following function

P (x, z, t) = A0 ei(krt) ,

(2.19)

solves the wave equation 2.12, where the wavenumber vector is given by
k = (kx , kz ) and the observation vector is given by r = (x, z). Similar to
the 1-D dispersion equation, the 2-D dispersion equation can be derived
by plugging 2.19 into equation 2.12 to get
kx2 + kz2 (/c)2 = 0.

(2.20)

The real part of equation 2.19, i.e., cos(kx x + kz z t), plots out as
straight lines perpendicular to the wavenumber vector k, and these lines
propagate in a direction parallel to k as t increases. This is easy to prove
because the general equation for a straight line is given by k r = cnst,
where k is a fixed vector perpendicular to the straight line. The locus of
points that satisfy this equation defines the wavefront where the phase
(i.e., = kx x + kz z t) is a constant. Thus as the time increases, i.e.
as cnst increases, the straight line also moves such that the direction
of movement is parallel to the fixed k vector, as shown in Figure 2.5.
Therefore equation 2.19 represents a harmonic plane wave propagating along the direction parallel to k. Similar to the discussion for a 1-D
plane wave, the shortest distance between two adjacent peaks of the
wavefront isqdefined to be the wavelength and is given by = 2/k,
where k = kx2 + kz2 is known as the wavenumber. Using this definition of k and that for the 2-D dispersion relation takes the same form
as equation 2.18. A good illustration of the relationship between the
wavenumber vector and the direction of wave propagation is given in
the movie. Here the left figure is a plot of points in (kx,kz) space, and
the right figure is the corresponding snapshot of the harmonic plane
wave. Note that as the length of the wavenumber vector increases the
wavelength decreases, and as the wavenumber direction changes so does
the direction of the propagating wave.

2.5. SOLUTIONS OF THE WAVE EQUATION

11

APPARENT WAVELENGTHS
T

x = /cos( )

z = /sin( )

Figure 2.5: 2-D wavefront propagating along the direction parallel to

the wavenumber k vector.

0.02

0.5

0.04

Time (s)

0.06

0.5
0.08

1
0.2

0.1
0.15

0.2
0.15

0.1
0.1
0.05
Time (s)

0.12

0.05
0

Offset (m)

0.05

0.1
Offset (m)

0.15

Figure 2.6: (Left) 3D seismograms and (right) 2D seismograms.

12CHAPTER 2. PHYSICS OF ACOUSTIC WAVE PROPAGATION

Incident Wave

Reflected Wave

ik z

-ik z
R e

x
z

T e

ikz

Transmitted Wave

Figure 2.7: Plane wavefront normally incident on a flat interface that

separates two homogeneous media. The unprimed medium indicates
that of the incident wave.

2.5.3

Plane Wave Propagation in a 2-Half Space

Medium

Figure 2.7 depicts a plane harmonic wave normally incident on an interface separating two half-spaces of unequal stiffness. The functions
are those for the up and downgoing solutions of the wave equation,
but it is understood that the geophones record the sum of the up- and
down-going wavefields, the total wavefield. That is, the total pressure
fields in the upper (+) and lower (-) media are expressed as
P + (z) = eikz + Rp eikz ,

(2.21)

P (z) = Tp eikz ,

(2.22)

where Rp and Tp denote the pressure reflection and transmission coefficients, respectively. The harmonic function eit has been harmlessly
dropped because it cancels out in the final expressions for Pp and Tp .
The two unknowns in these linear equations, Rp and Tp , can be
determined by imposing two equations of constraints at the interface

2.5. SOLUTIONS OF THE WAVE EQUATION

13

at z = 0: continuity of pressure across the interface

P + (z = 0) = P (z = 0) [eikz + Rp eikz ]|z=0 = Tp eikz |z=0 ,
or
1 + R p = Tp ,
(2.23)
and continuity of particle velocity (recall Newtons Law in equation 2.8)
across the interface
1/P + /z|z=0 = 1/0 P /z|z=0 1//z(eikz + Rp eikz )|z=0
= 1/0 /z(Tp eikz )|z=0 ,
or
(k/)(1 Rp ) = (k 0 /0 )Tp ,
(2.24)
Setting k = /c and k 0 = /c0 , and solving for Rp in equations 2.23
and 2.24 yields the pressure reflection and transmission coefficients for
a normally incident plane wave on a flat interface:
Rp = (0 c0 c)/(0 c0 + c),

(2.25)

Tp = 20 c0 /(0 c0 + c).

(2.26)

Here c is known as the impedance of the medium, and roughly indicates the stiffness of a medium. For example, a plane harmonic plane
wave in a homogeneous medium exerts a pressure denoted by P =
eikxit and has a particle velocity denoted by u = k/()P = 1/(c)P .
Therefore, the ratio P/u becomes
P/u = c.

(2.27)

This says that decreasing impedances imply larger particle velocities for
a fixed elastic pressure on a cubes face. This is exactly what one would
expect in a really soft medium: larger displacements for springier-softer
rocks, which is one of the reasons that earthquakes shake sediment-filled
valleys more than the surrounding bedrock. Conversely, stiffer media
lead to smaller displacements for a given elastic pressure.
Note:

14CHAPTER 2. PHYSICS OF ACOUSTIC WAVE PROPAGATION

P = 0 at the Free Surface

-ikz

Upgoing

x
z

- e

ikz

Downgoing

-ikz

- e

ikz

Total Field
Node

Figure 2.8: The pressure at the free surface is always zero because the
air has no stiffness to resist motion. Mathematically, the downgoing
wave has equal and opposite amplitude to the upgoing wave at the free
surface.
The pressure reflection coefficient is negative if the impedance of
the incident layer is greater than that of the refracting layer, i.e.,
c > 0 c0 . Thus gas sands (which typically have lower velocity
than the overlying brine sand) have negative polarity reflections.
The free-surface reflection coefficient due to an upcoming wave
is Rp = 1 because the rock impedance of the incident layer is
greater than that of air (impedance of air 0 so c > 0 c0 = 0).
Equation 2.21 says that the total pressure field value on the free
surface is P = 1 + R= 1 1 = 0! See Figure 2.8.
In a land experiment geophones record particle velocity of the
ground while a marine experiment records pressure with hydrophones.
If P = 0 on the free surface then we must lower the hydrophones
several feet beneath the water surface, otherwise we record nothing.
The pressure transmission coefficient Tp is larger than 1 if the incident medium has a lower impedance than the refracting medium,

2.5. SOLUTIONS OF THE WAVE EQUATION

15

i.e. c < 0 c0 . For a given displacement, a really stiff medium

exerts a greater elastic pressure than a soft medium. Check out
the a href=fd1mov.mpgmovie/a of a 1-D propagating wave
and see if the transmission amplitude changes and the polarity of
the reflection changes according to the reflection and transmission
formulae derived above.

2.5.4

Reflection Coefficient for Particle Velocity

Marine experiments measure the pressure field, so this is why the hydrophones must be sufficiently below the sea surface in order to measure
a non-zero pressure. On the other hand, land experiments use geophones that measure the particle velocity. Typically, only the vertical
component of particle velocity is measured. The reflection coefficient
for particle velocity has a different form than that for pressure. To see
this, assume that the up and downgoing vertical particle-velocity fields
in the top and bottom layers are given by
w + (z) = eikz + Rw eikz ,
0
w (z) = Tw eik z .

(2.28)

The boundary conditions at the interface are continuity of vertical /z|z=0 ,

particle velocity w + |z=0 = w |z=0 and pressure w + /z|z=0 = 0 w
2
where Hookes law is used w = P/z = c P/z.
These two continuity conditions yield the following boundary conditions:
1 + Rw = Tw ,
c(1 Rw ) = (1 + Rw )0 c0 ,

(2.29)

which can be solved for the particle velocity reflection and transmission
coefficients:
Rw = (c 0 c0 )/(0 c0 + c),

(2.30)

Tw = 2c/(0 c0 + c),

(2.31)

16CHAPTER 2. PHYSICS OF ACOUSTIC WAVE PROPAGATION

where the unprimed variables again refer to the medium of the incident
wave, and the subscript w denotes vertical-particle velocity. Note that
the reflection coefficient above will have opposite polarity compared
to the pressure reflection coefficients. Also, note that in some cases
the transmitted amplitude can be greater than the amplitude of the
incident wave! Does this violate conservation of energy? No, energy is
the squared modulus of amplitude scaled by the impedance (see later
section). Thus, a weaker medium with weak rocks (small impedance)
can have transmit larger amplitude waves than the incident waves in a
much stronger (larger impedance) medium. It takes much more energy
to rapidly displace strong rock 1 mm than it does in a weak rock.
Free-Surface Reflection Coefficent. The particle velocity reflection
coefficient is equal to +1 at the free surface, so the total particle velocity
field at the free surface is 1 + Rw = 2. Thus the free surface, because
it straddles a zero stiffness medium, can oscillate with great vigor and
has the largest amplitude compared to the underlying rock motion in
a half-space. When an earthquake hits, dig a hole, jump in, and cover
yourself with dirt! See Figure 2.9.

2.6

Oblique Incidence Angles, Reflection

Coefficients and Snells Law

If the incidence angle of the plane wave is non-zero (measured w/r to

the normal) then the normal derivatives in equation 2.24 bring down
kz rather than k. Therefore, equation 2.24 becomes
(kz /)(1 Rp ) = (kz0 /0 )Tp ,
(2.32)
or kz = 2cos()/ and kz0 = 2cos0 /0
(cos/())(1 Rp ) = (cos0 /(0 0 ))Tp ,

(2.33)

and because = c/f we have

(cos/(c))(1 Rp ) = (cos0 /(c0 0 ))Tp .

(2.34)

2.6. OBLIQUE INCIDENCE ANGLES, REFLECTION COEFFICIENTS AND SNELLS LAW17

Solving for the reflection coefficient yields
Ru = (cos0 c cos0 c0 )/(cos0 c0 + cos0 c),

(2.35)

If the horizontal wavenumbers in the upper (kx = sin/c) and

lower medium (kx0 = sin/c0 ) are equated, kx0 = kx , then this implies
Snells law:
sin /c = sin0 /c0 .

(2.36)

This means that transmission rays bend across an interface, bending towards the vertical when entering a slower velocity medium and bending
towards the horizontal when entering a faster medium (see Figure 2.10).
At the critical incidence angle crit the refraction angle of the transmitted ray 0 is 90 degrees so that Snells law says crit = arcsin(c/c0 ) if
c0 > c. This gives rise to refraction head waves that propagate parallel
to the interface at the velocity c0 of the underlying medium.
A consequence of Snells law is that a medium with a velocity that
increases linearly with depth always turns a downgoing ray back towards the surface, as shown in Figure 2.10. This can easily be shown
by approximating the linear velocity gradient medium with a stack of
thinly-spaced layers, each with a homogeneous velocity that slightly increases with depth. The velocity increase is the same across each layer.
Applying Snells law to a downgoing ray shows that each ray transmitted across an interface bends a little bit more towards the horizontl
until it goes back up. As the thickness of eash layer decreases, the ray
trajectory will be the arc of a circle.
Why did we equate the upper and lower medium wavenumbers? We
did this because the apparent horizontal velocities along the interface
must be equal to one another (i.e., vx = /kx = vx0 = /kx0 ). The horizontal velocities of the incident and transmitted waves must be equal
at the interface otherwise they would break away from one another and
not be coupled at the interface.
How does the formula for transmission coefficient in equation 2.31
change for an oblique incidence angle?

18CHAPTER 2. PHYSICS OF ACOUSTIC WAVE PROPAGATION

Vertical Displacement at the Free Surface

-ikz

+e

Upgoing

ikz

Downgoing

-ikz

+ e

ikz

Total Field

Anti-Node

Figure 2.9: The vertical displacement at the free surface is maximum

because no elastic resistance at z = 0 means vigorous ground shaking. Mathematically, the downgoing wave has equal amplitude to the
upgoing wave at the free surface.

Linear Increasing Velocity V(z) Medium

V(z+dz) > V(z)
dz

V1
V2
V3

V4

sin = sin
V3

V4

Figure 2.10: Downgoing rays bend across an interface towards the

horizontal if the velocity increases with depth. For a medium where
c(z) = a + bz, all downgoing rays eventually bend back towards the
surface.

2.7. ENERGY OF PROPAGATING ACOUSTIC WAVES

2.7

19

Energy of Propagating Acoustic Waves

Elastic energy is stored in a cube of acoustic material as it is deformed

from equilibrium. That is, squeeze a cube of acoustic material, release
it, and then the cube moves to perform work on the medium. In the
deformed state the potential to perform work takes the form of elastic
potential energy. Figure 2.11 shows that the work (i.e., area force/area
distance cube moves) performed byR the surrounding medium on a
cube along the z axis is given by (P dxdy)dz, where the limits of
integration are from the undeformed volume to the deformed volume
at some given time. This figure shows that, using equation 2.3 and 2.2
and = c2 , the expression for instantaneous potential energy density
is given by
PE

P 2 /(2c2 ).

(2.37)

Therefore the total instantaneous energy density of an acoustic wave

propagating along the x-axis in a homogeneous medium is given by
= E/V

= 1/2[u 2 + P 2 /(c2 )],

(2.38)

where the kinetic energy is given by the second to the last term and
the potential energy is given by the last term. For a harmonic plane
wave c = P/u this equation becomes:
= u 2 .

(2.39)

As one might expect, it takes more energy to move denser rock with
the same particle velocity as moving lighter rock.

2.8

Elastic Wave Equation

If the medium has non-zero shear strength then there can be shear
strains supported by the rock. This means that the shape of a cube
can be distorted into a, e.g., trapezoidal-like shape after application
of a shear stress on the cube. Although most of our treatment of
exploration seismology will assume the acoustic approximation, we will

20CHAPTER 2. PHYSICS OF ACOUSTIC WAVE PROPAGATION

Instantaneous Potential Energy Density

deformed

[ P dx dy ]dz = - P dV = V P dP = 1/2 P V
undeformed
c2
c2
Hookes
Law
-

dy
dz

dx

undeformed
deformed

dV = - V dP
c2

Figure 2.11: The potential energy of a cube of deformed material along

the z-axis.
need to acknowledge the underling physics of elastic wave propagation
when dealing with surface waves and shear waves in our data.
The notation for the shear strain is ij = 1/2(ui /xj + uj /xi ),
which says that a shear strain exists if there is a non-zero gradient
of displacement that is perpendicular to the direction of displacement.
We can also use Einstein index notation so that ij = 1/2(ui,j + uj,i ),
where the index follwoing a comma indicates a partial derivative with
respect to that indiexs coordinate. Figure ?? illustrates the idea of a
shear strain that distorts the shape of a cube.
Hookes Law for elastic strain in an isotropic medium Aki and Rcihards, 1980) is given by
ij = ij kk + 2ij ,

(2.40)

where ij is the stress tensor thgat denotes the force/area imposed by

the outside medium in the ith direction along the face with normal j.
Also, is the shear modulus and is Lames constant. A consequence
of a medium supporting shear strain is that the shear stresses can be
non-zero.
Newtons law is given by
ui = i1 /x1 + i2 /x2 + i3 /x3 = ij,j ,

(2.41)

2.8. ELASTIC WAVE EQUATION

21

where Eintein notation says that repeated indicies indicate summation

over all three components. Inserting equation 2.40 into equation 2.41
gives the elastic wave equation in terms of components only:
ui = ( + )jj,i + ui,jj

(2.42)

or in vector notation:
= c2p ( u) c2s ( u),
(2.43)
u

where the compressional velocity is given by cp = + 2 and the

shear velocity is given by cs = .

For a homogeneous medium we can decompose u into a sum of
compressional and shear potentials u = + , and plug it
into equation 2.43 to get the wave equation for a compressional wave
propagating at velocity cp :
= c2p 2

(2.44)

and a shear wave propagating at velocity cs :

= c2s 2

(2.45)

Here we use the identities = 0, phi = 0, and 2 =

( ) ( ).

2.8.1

Reflection Coefficients

Figure 2.12 depicts the plane wave impinging upon a horizontal elastic
interface. Here, there are 5 different wave types to consider because
a shear wave can be generated at the interface. The shear wave has
particle motion that is perpedicular to the propagation direction while
the compressional components is parallel to the direction of propagation. Imposing continuity of vertical and horizontal particle velocity
and normal and shear stress tensors provide four equations of constraint. Similar to the acoustic case, we can solve these four equations
for the unknown amplitudes P P, P S, P S 0 , andP P 0 . The P P reflection
coeficcient at the free surface is given in Aki and Richards (1980):
PP =

4 4 p2 cos(i)cos(j)/() (1 2 2 p2 )2
4 4 p2 cos(i)cos(j)/() + (1 2 2 p2 )2

(2.46)

22CHAPTER 2. PHYSICS OF ACOUSTIC WAVE PROPAGATION

Plane Wave Reflection Coefficients (4 unknowns)
P

PS
PP

PP

PS

Figure 2.12: Plane wave impinging upon a horizontal interface.

where , , and indicate P-wave velocity, S-wave velocity and density
respectively.

2.8.2

Small Angle Approximation to Reflection Coefficients

For a two layered medium separated by a horizontal interface the PP

reflection coefficient for small jumps in the medium parameters is given
by (Aki and Richards, 1980, p. 153):

4 2 p2

,
P P = .5(1 4 p )/ +
2 cos2 i

2 2

(2.47)

where = 2 1 , = 2 1 , and = 2 1 . here the

1 subscript indicates the upper medium where the plane P-wave is
incident and 2 indicates the lower medium.
Further manipulation of this equation by defining the Poisson factors
= 1 2 ; 2 = 2 (1 2)/(2[1 ]),
= (1 + 2 )/2 ; = (1 + 2 )/2;

(2.48)

2.8. ELASTIC WAVE EQUATION

23

gives
R() = A + B sin2 + C(tan2 sin2 )

(2.49)

where
A = R0 ; R0 1/2(/ + /); B = B0 +

B0 = C 2(1 + C)

;
(1 2 )R0

1 2
/
; C=
;
1
/ + /

(2.50)

(2.51)

The first angle-dependent term in equation 2.49 significantly contrubutes

for 0 < < 30 degrees, while the second starts to significantly contribute for > 30 degrees.
For < 30 degrees, geophysicists will use the small angle approximation to equation 2.49:
R() A + B sin2 ,

(2.52)

and plot up crossplot curves (Foster et al., 1997) to assess geology. For
example, Figure 2.13 depicts the crossplot of A and B pairs taken from
a well log (i.e., estimate density and P- and S-wave velocities from sonic
log). It shows a linear trend, and the idea is that any deviations from
this trend represent a significant change of geology such as oil or gas
bearing rocks. The departures can be estimated by finding A and B
pairs from the R() vs curves estimated from the seismic reflection
amplitudes along a horizon of interest.
Another measure of hydrocarbon anomalies can be estimated by
finding the linearization of the fluid factor (Shuey, 1985):
Rf () = R0 + 1/2/tan2 ,

(2.53)

where tan2 sin2 for small angles. This factor is used in the case
history described in the next chapter.
Details for implementing this AVO (i.e., Amplitude vs Offset) procedure are non-trivial because much data processing must be performed
before the A and B pairs can be picked. Nevertheless, significant oil
and gas deposits have been discovered by the AVO method.

24CHAPTER 2. PHYSICS OF ACOUSTIC WAVE PROPAGATION

Figure 2.13: Crossplot of slope (B) and intercept (A) pairs using well
log data from the North Sea. The dashed line corresponds to the fluid
line for / = 1.9. The A and B pairs tend to lie on a trend that is
consistent with the fluid line. Note, the fluid line assumes a constant
density, the actual density profile is used to compute A and B.

Bibliography
[1] Aki, K., and Richards, P., 1980, Quantitative seismology: Theory
and Methods: Freeman Co., NY, NY.
[2] Foster, D., Keys, R., and Reilly, J., 1997, Another perspective on
AVO crossplotting: The Leading Edge, September issue.
[3] Kinsler, L. and Frey, A., 1961, Fundamentals of Acoustics, . ley and
Sons, NY, NY.
[4] Shuey, R.T., 1985, Application of the Zoeppritz equations: Geophysics, 609-614.

2.9

Problems

1. Identify the direct arrival, air wave, surface waves, refraction arrivals, and reflection arrivals in the CSG shown in Figure 2.14.
Estimate the apparent velocity in the x-direction Vx and the associated period for each event. From these calculations determine
the wavelengths. Show work.
2. Which arrivals have the same apparent velocity as the actual
propagation velocity of that event? Why?
3. The 1-D SH wave equation is the same form as the 1-D acoustic
wave equation, except c becomes the shear wave velocity, P becomes the y-component of displacement v, c=sqrt(mu/rho) where
mu is the shear modulus, and the SH wave equation is
1/c2 2 v/t2 2 v/z 2 = 0
25

(2.54)

26

BIBLIOGRAPHY

Figure 2.14: Shot gather from salt lake valley. The trace interval is 5
feet along the horizontal axis and the time units along the vertical axis
are seconds.
SH (or shear horizontal) refers to the fact that the shear wave
particle motion is perpendicular to the direction of particle motion, and is along the horizontal direction (in and out of plane
of paper). The SH continuity conditions at the interface at z=0
are a). Continuity of y-displacement v + = v ., b). Continuity
of shear traction: v/z + = v/z , where is the shear
modulus.
Derive the y-displacement reflection and transmission coefficients
for a plane SH wave normally incident on a planar interface in an
elastic medium.