next up previous print clean
Next: New notations and definitions Up: Appendix Previous: Appendix

The wave equation in an acoustic medium

Let xi, where i=1,2,3 be three orthogonal directions and $\vec x$the position vector in a coordinate system associated with the three directions. Let us define $\rho = \rho \left( {\vec x} \right)$ as the mass per volume unit in the acoustic medium, $\vec u = \vec
u\left( {\vec x} \right)$ as the velocity of the acoustic medium and K as the bulk modulus of the acoustic medium. The second law of dynamics states that mass $\times$ acceleration = force = - pressure gradient:

 
 \begin{displaymath}
\rho \frac{{\partial u_i }}{{\partial t}} = - \frac{{\partial
P}}{{\partial x_i }}.\end{displaymath} (15)

Energy can be stored by compression and volume variation. If

 
 \begin{displaymath}
\vec u\left( {\vec x + \delta \vec x} \right) \ne \vec u\left( {\vec
x} \right),\end{displaymath} (16)

we say that the flow diverges (the volume changes). This leads to a pressure variation, proportional to the divergence of the velocity:

 
 \begin{displaymath}
 - \frac{{\partial P}}{{\partial t}} = K\left( {\frac{{\part...
 ...ial x_2 }} +
 \frac{{\partial u_3 }}{{\partial x_3 }}} \right).\end{displaymath} (17)

The wave equation in an acoustic medium can be deduced from (15) and (17) as follows. Derivate (17) with respect to time:

 
 \begin{displaymath}
\frac{{\partial ^2 P}}{{\partial t^2 }} = - K\sum\limits_{i = 1}^3
{\frac{{\partial ^2 u_i }}{{\partial t\partial x_i }}}.\end{displaymath} (18)

Divide (15) by $\rho$ and derivate it with respect to the axis xi:

 
 \begin{displaymath}
\frac{{\partial ^2 u_i }}{{\partial t\partial x_i }} = -
\fr...
 ...tial x_i }}\frac{1}{\rho }\frac{{\partial
P}}{{\partial x_i }}.\end{displaymath} (19)

Plug (19) in (18):

 
 \begin{displaymath}
\frac{{\partial ^2 P}}{{\partial t^2 }} = K\sum\limits_{i = ...
 ...al x_i }}\frac{1}{\rho }\frac{{\partial
P}}{{\partial x_i }}}. \end{displaymath} (20)

Approximation: $\rho$ is a constant that does not depend on the position vector. By denoting the acoustic waves propagation velocity through the medium by v, where

 
 \begin{displaymath}
v^2 = \frac{K}{\rho },\end{displaymath} (21)

we obtain the acoustic wave equation:

 
 \begin{displaymath}
\frac{{\partial ^2 P}}{{\partial t^2 }} = v^2 \sum\limits_{i = 1}^3
{\frac{{\partial ^2 P}}{{\partial x_i ^2 }}}.\end{displaymath} (22)


next up previous print clean
Next: New notations and definitions Up: Appendix Previous: Appendix
Stanford Exploration Project
6/8/2002