A one-way wave equation is a first-order partial differential equation describing one wave traveling in a direction defined by the vector wave velocity. It contrasts with the second-order two-way wave equation describing a standing wavefield resulting from superposition of two waves in opposite directions (using the squared scalar wave velocity).[1] [2] [3] In the one-dimensional case it is also known as a transport equation,[4] and it allows wave propagation to be calculated without the mathematical complication of solving a 2nd order differential equation. Due to the fact that in the last decades no general solution to the 3D one-way wave equation could be found, numerous approximation methods based on the 1D one-way wave equation are used for 3D seismic and other geophysical calculations, see also the section .[5] [6]
The scalar second-order (two-way) wave equation describing a standing wavefield can be written as:where
x
t
s=s(x,t)
c
Due to the ambiguity in the direction of the wave velocity,
c2=(+c)2=(-c)2
+x
-x
where
s+
s-
+c
-c
When a one-way wave problem is formulated, the wave propagation direction has to be (manually) selected by keeping one of the two terms in the general solution.
Factoring the operator on the left side of the equation yields a pair of one-way wave equations, one with solutions that propagate forwards and the other with solutions that propagate backwards.[7]
The backward- and forward-travelling waves are described respectively (for
c>0
The one-way wave equations can also be physically derived directly from specific acoustic impedance.
In a longitudinal plane wave, the specific impedance determines the local proportionality of pressure
p=p(x,t)
v=v(x,t)
with
\rho
The conversion of the impedance equation leads to:[3]
A longitudinal plane wave of angular frequency
\omega
s=s(x,t)
The pressure
p
v
s
E
\sigma
\sigma=E\varepsilon
\varepsilon=
\DeltaL | |
L |
These relations inserted into the equation above yield:
With the local wave velocity definition (speed of sound):
directly(!) follows the 1st-order partial differential equation of the one-way wave equation:
The wave velocity
c
+c
-c
For wave propagation in the direction of
+c
and for wave propagation in the
-c
There also exists a spherical one-way wave equation describing the wave propagation of a monopole sound source in spherical coordinates, i.e., in radial direction. By a modification of the radial nabla operator an inconsistency between spherical divergence and Laplace operators is solved and the resulting solution does not show Bessel functions (in contrast to the known solution of the conventional two-way approach).[12]
The one-way equation and solution in the three-dimensional case was assumed to be similar way as for the one-dimensional case by a mathematical decomposition (factorization) of a 2nd order differential equation.[13] In fact, the 3D One-way wave equation can be derived from first principles: a) derivation from impedance theorem and b) derivation from a tensorial impulse flow equilibrium in a field point. It is also possible to derive the vectorial two-way wave operator from synthesis of two one-way wave operators (using a combined field variable). This approach shows that the two-way wave equation or two-way wave operator can be used for the specific condition ∇c=0, i.e. for homogeneous and anisotropic medium, whereas the one-way wave equation resp. one-way wave operator is also valid in inhomogeneous media.[14]
For inhomogeneous media with location-dependent elasticity module
E(x)
\rho(x)
c(x)
The method of PDE factorization can also be transferred to other 2nd or 4th order wave equations, e.g. transversal, and string, Moens/Korteweg, bending, and electromagnetic wave equations and electromagnetic waves.