Script: ex_vector_filter_dip.py

Description

This ExternalAttrib script can be used to apply a vector filter to orientation ( inline and crossline dip) data. The script offers a choice of mean vector, L1 vector median and L2 vector median filters.

Initially the inline and crossline dip data are converted to a normal vector to the local orientation: $ [x_i, y_i, z_i] $.

The Mean Vector Filter averages each of the vector components of the orientation normal vectors in the analysis cube: $$ \Big[x_f, y_f, z_f\Big] = \frac{1}{N} \Big[\sum\limits_{i}^N x_i, \sum\limits_{i}^N y_i, \sum\limits_{i}^N z_i\Big] $$

The L1 vector filter finds the normal vector in the analysis cube whose sum of absolute distance from all the others is a minimum: $$ \Big[x_f, y_f, z_f\Big] = argmin \sum\limits_{i}^N \Big[|x_f-x_i| + |y_f-y_i| + |z_f-z_i|\Big] $$

The L2 vector filter finds the normal vector in the analysis cube whose sum of squared distance from all the others is a minimum: $$ \Big[x_f, y_f, z_f\Big] = argmin \sum\limits_{i}^N \Big[(x_f-x_i)^2 + (y_f-y_i)^2 + (z_f-z_i)^2\Big] $$

The filtered orientation can be output as any of the following:

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.

The script requires the Numba Python package.

Examples


Unfiltered phase dip - crossline dip on an inline
Unfiltered phase dip - crossline dip on an inline


Mean vector filtered phase dip - 3x3x3 (Stepout and ZWindow of 1)
Mean vector filtered phase dip - 3x3x3 (Stepout and ZWindow of 1)


L1 vector median filtered phase dip - 3x3x3 (Stepout and ZWindow of 1)
L1 vector median filtered phase dip - 3x3x3 (Stepout and ZWindow of 1)


L2 vector median filtered phase dip - 3x3x3 (Stepout and ZWindow of 1)
L2 vector median filtered phase dip - 3x3x3 (Stepout and ZWindow of 1)

Input Parameters


ex_vector_filter_dip.py input parameters
ex_vector_filter_dip.py input parameters

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.