Wavelets

A collection of seismic data processing wavelet

Berlage

# Seismic.BerlageFunction.

Berlage(; <keyword arguments>)

Create a Berlage wavelet.

Arguments

Keyword arguments

  • dt::Real=0.002: sampling interval in secs.
  • f0::Real=20.0: central frequency in Hz.
  • m::Real=2: exponential parameter of Berlage wavelet.
  • alpha::Real=180.0: alpha parameter of Berlage wavelet in rad/secs.
  • phi0::Real: phase rotation in radians.

Example

julia> w = Berlage(); plot(w);

Reference

  • Aldridge, David F., 1990, The berlage wavelet: GEOPHYSICS, 55, 1508–1511.

source

Ormsby

# Seismic.OrmsbyFunction.

Ormsby(; <keyword arguments>)

Create a Ormsby wavelet sampled every dt seconds with corner frequencies defined by the vector f = [f1, f2, f3, f4]. The final wavelet is multiplied by a Hamming window.

Arguments

Keyword arguments

  • dt::Real=0.002: sampling interval in secs.
  • f::Vector{Real}=[2.0, 10.0, 40.0, 60.0]: corner frequencies in Hz.

    ^

    1 | *** | * * | * * | * * | * * ––––––––––––––-> f f1 f2 f3 f4

Example

julia> w = Ormsby(); plot(w);

source

Ricker

# Seismic.RickerFunction.

Ricker(; <keyword arguments>)

Create a Ricker wavelet.

Keyword arguments

  • dt::Real=0.002: sampling interval in secs.
  • f0::Real=20.0: central frequency in Hz.

Examples

julia> w = Ricker(); plot(w);
julia> w = Ricker(dt=0.004, f0=20); plot(w);

Reference

Sheriff, Robert, 2002, Encyclopedic Dictionary of Applied Geophysics, fourth ed.: Society of Exploration Geophysicists. Geophysical Reference Series No. 13.

source