geost.analysis.find_top_sand#

geost.analysis.find_top_sand(lith: ndarray, top: ndarray, bottom: ndarray, min_sand_frac: float, min_sand_thickness: float) → float[source]#

Find the top of sand depth in a borehole. The top of sand is defined by the first layer of a specified thickness that contains a minimum percentage of sand. By default: when the first layer of sand is detected, the next 1 meter is scanned. Within this meter, if more than 50% of the lenght has a main lithology of sand, the initially detected layer of sand is regarded as the top of sand. If not, continue downward until the next layer of sand is detected and repeat.

Parameters:

  • lith (ndarray) – Numpy array containing the lithology of the borehole.

  • top (ndarray) – Numpy array containing the top depth of the layers of the borehole.

  • bottom (ndarray) – Numpy array containing the bottom depth of the layers of the borehole.

  • min_sand_frac (float) – Minimum percentage required to be sand.

  • min_sand_thickness (float) – Minimum thickness of the sand to search for.

Returns:

top_sand – Top depth of the sand layer that meets the requirements.

Return type:

float