Dip and Azimuth

various scripts to estimate orientation

Description

These Python External Attribute scripts implement various algorithms to estimate orientation, ie dip or dip azimuth.

All scripts will estimate at least the following attributes:

OUTPUT DESCRIPTION
Inline Dip Event dip observed on a crossline in microseconds per metre for time surveys and millimetres per metre for depth surveys. Output can be positive or negative with the convention that events dipping towards larger inline numbers producing positive dips.
Crossline Dip Event dip observed on an inline in microseconds per metre for time surveys and millimetres per metre for depth surveys. Output can be positive or negative with the convention that events dipping towards larger crossline numbers producing positive dips.
True Dip Event dip in microseconds per metre for time surveys and millimetres per metre for depth surveys. Output is always positive.
Dip Azimuth Azimuth of the True Dip direction relative to the survey orientation. Output ranges from -180 to 180 degrees. Positive azimuth is defined from the inline in the direction of increasing crossline numbers. Azimuth = 0 indicates that the dip is dipping in the direction of increasing crossline numbers. Azimuth = 90 indicates that the dip is dipping in the direction of increasing inline numbers.

Some scripts may offer additional outputs such as a measure of event coherency or planarity.

All of the scripts require the numba Python package.

Orientation from gradients

Script: ex_gradient_dip.py

Calculates orientation from inline, crossline and Z gradients. No filtering is applied

Unfiltered gradient dip - crossline dip on an inline

Unfiltered gradient dip - crossline dip on an inline

Orientation from vector filtered gradients

Script: ex_vf_gradient3_dip.py

Uses Kroon’s (2009) 3 point derivative filter to estimate data gradients. Next gradient normal unit vectors are determined and smoothed using a vector filter.

NAME DESCRIPTION
Output What to calculate - choice of inline dip, crossline dip, true dip or dip azimuth.
Z window (+/-samples) Specifies the extent of the analysis cube in the Z direction. Number of Z samples in cube will be \((2*Zwindow+1)\).
Stepout Specifies the inline and crossline extent of the analysis cube. Number of samples in each direction will be \((2*Stepout+1)\).
Filter Choice of Mean Dip, L1 Vector Median or L2 Vector Median.

The aperture of the vector filtering is \((2*ZStepOut-1)\) \(Z\) samples and \((2*Stepout-1)\) samples in the inline and crossline direction. For example for a 5x5x5 analysis cube \((Zwindow=2, Stepout=2)\) the gradients and associated normal unit vectors are generated on a 3x3x3 cube and vector filtered. The outer samples are only used in the gradient calculation.

Vector filtered gradient dip - crossline dip on an inline - 5x5x5 input

Vector filtered gradient dip - crossline dip on an inline - 5x5x5 input

Orientation by the gradient structure tensor

Scripts: ex_gradient3_st_dip.py & ex_gradient5_st_dip.py

Uses either Kroon’s (2009) 3 point or the Farid and Simoncelli (2004) 5 point derivative filter to estimate data gradients which are then used to form the gradient structure tensor.

NAME DESCRIPTION
Output What to calculate - choice of inline dip, crossline dip, true dip or dip azimuth.
Z window (+/-samples) Specifies the extent of the analysis cube in the Z direction. Number of Z samples in cube will be \((2*Zwindow+1)\).
Stepout Specifies the inline and crossline extent of the analysis cube. Number of samples in each direction will be \((2*Stepout+1)\).

For the ex_gradient3 script the structure tensor is formed from an aperture of \((2*ZStepOut-1)\) \(Z\) samples and \((2*Stepout-1)\) samples in the inline and crossline direction.

For the ex_gradient5 script the structure tensor is formed from an aperture of \((2*ZStepOut-2)\) \(Z\) samples and \((2*Stepout-1)\) samples in the inline and crossline direction.

Gradient3 structure tensor dip - crossline dip on an inline - 5x5x5 input

Gradient3 structure tensor dip - crossline dip on an inline - 5x5x5 input

Orientation from the 3D complex trace phase

Script: ex_phase3_dip.py

Calculates orientation from the 3D complex trace phase gradients as per Barnes (2007). Kroon’s (2009) 3 point derivative filter is used to compute gradients.

NAME DESCRIPTION
Output What to calculate - choice of inline dip, crossline dip, true dip or dip azimuth.
Z window (+/-samples) Specifies the length \((2*Zwindow+1)\) of the time domain operator used to generate the complex analytic signal (recommend >= 15)
Band Specifies the proportion of the frequency band to include when generating the complex analytic signal (recommend 0.9).

Unfiltered phase dip - crossline dip on an inline

Unfiltered phase dip - crossline dip on an inline

Orientation from vector filtered 3D complex trace phase

Script: ex_vf_phase3_dip.py

Calculate orientation unit normal vectors using the 3D complex trace phase gradient and apply a vector filter. Kroon’s (2009) 3 point filter used to compute gradients.

NAME DESCRIPTION
Output What to calculate - choice of inline dip, crossline dip, true dip or dip azimuth.
Z window (+/-samples) Specifies the length \((2*Zwindow+1)\) of the time domain operator used to generate the complex analytic signal (recommend >= 15)
Stepout Specifies the inline and crossline extent of the analysis cube. Number of samples in each direction will be \((2*Stepout+1)\).
Filter Choice of Mean Dip, L1 Vector Median or L2 Vector Median.
Vector Filter ZStepOut Specifies the extent of the analysis cube for vector filtering in the Z direction. Number of Z samples in cube will be \((2*ZStepOut+1)\).
Band Specifies the proportion of the frequency band to include when generating the complex analytic signal (recommend 0.9).

The aperture of the vector filter is \((2*ZStepOut+1)\) \(Z\) samples and \((2*Stepout-1)\) samples in the inline and crossline direction.

Mean Vector Filtered phase dip - crossline dip on an inline - 3x3x3

Mean Vector Filtered phase dip - crossline dip on an inline - 3x3x3

Orientation using the envelope weighted 3D complex trace phase structure tensor

Script: ex_weighted_phase3_st_dip.py

Forms a structure tensor from the 3D complex trace phase gradients. Tensor elements are weighted by the trace envelope as per Luo etal (2006). Kroon’s (2009) 3 point filter is used to compute gradients.

NAME DESCRIPTION
Output What to calculate - choice of inline dip, crossline dip, true dip or dip azimuth.
Z window (+/-samples) Specifies the length \((2*Zwindow+1)\) of the time domain operator used to generate the complex analytic signal (recommend >= 15)
Stepout Specifies the inline and crossline extent of the analysis cube. Number of samples in each direction will be \((2*Stepout+1)\).
Tensor ZStepOut Specifies the extent of the analysis cube for vector filtering in the Z direction. Number of Z samples in cube will be \((2*ZStepOut+1)\).
Band Specifies the proportion of the frequency band to include when generating the complex analytic signal (recommend 0.9).

The aperture of the structure tensor is \((2*ZStepOut+1)\) \(Z\) samples and \((2*Stepout-1)\) samples in the inline and crossline direction.

Structure tensor phase dip - crossline dip on an inline - 3x3x3 tensor input

Structure tensor phase dip - crossline dip on an inline - 3x3x3 tensor input


Last modified April 17, 2022: First update (afe101d)