Scientific Code

modelr.wavelet — Module Containing wavelet functions.

TODO: long doc here

modelr.wavelet.klauder(points, f_lower, f_upper)

The autocorrelation of a linearly swept frequency-modulated sinuusoidal signal used in Vibroseis. It is zero-phase, and symmetrical about time t = 0 ms. Requires f_upper and f_lower, and a duration, which we can interpret from points.

modelr.wavelet.ormsby(points, f1, f2, f3, f4)

The Ormsby wavelet requires four frequencies: f1 = low-cut frequency f2 = low-pass frequency f3 = high-pass frequency f4 = hi-cut frequency

Together, the frequencies define a trapezoid shape in the spectrum. The Ormsby wavelet has several sidelobes, unlike Ricker wavelets which only have two, one either side.

modelr.wavelet.ricker(points, a)

Also known as the “mexican hat wavelet”, models the function:

modelr.wavelet.ricker_freq(points, f)

Also known as the “mexican hat wavelet”, models the function:

A = (1-2 pi^2 f^2 t^2) e^{-pi^2 f^2 t^2}

The amplitude ‘’A’’ of the [[Ricker wavelet]] with peak frequency ‘’f’’ at time ‘’t’’ is computed like so:

Params points:adsf
Params f:asdf
Returns:adf

modelr.reflectivity – TODO short docstring here

TODO: long doc here

modelr.reflectivity.create_theta(pad, thickness, prop0, prop1, theta, f, points, reflectivity_method)

Create a 2D array where the first dimension is time and the second is angle.

Parameters:
  • pad – pad the array top and bottom in ms
  • thickness – the distanc between the two interfaces
  • prop0 – rock properties 1
  • prop1 – rock properties 1
  • theta – array of angles
  • f – the frequency for the wavelet
modelr.reflectivity.create_wedge(ntraces, pad, max_thickness, prop0, prop1, theta, f)

Create a wedge model.

Parameters:
  • ntraces – number of traces
  • pad – pad the array top and bottom in ms
  • max_thickness – The thickest part of the wedge
  • prop0 – rock properties 1
  • prop1 – rock properties 1
  • theta – angle of incidence
  • f – the frequency for the wavelet

modelr.rock_properties – TODO short docstring here

TODO: long doc here

class modelr.rock_properties.RockProperties(vp, rho, vs=None)

Class to store rock properties.

Parameters:
  • vp – pressure wave coef
  • vs – shear wave coef
  • rho – TODO
modelr.rock_properties.zoeppritz(Rp0, Rp1, theta1)

Wrapper around long zoeppritz funcion.

Parameters:
  • Rp0
  • Rp1
  • theta1
modelr.zoep_calc.zoeppritz(vp1, vs1, rho1, vp0, vs0, rho0, theta1)

Documentation goes here.