_morph_drain_config

_morph_drain_config(image, radius)

Estimate the multiphase configuration of the pore space using the morphological drainage method following Hilper & Miller (2001): :rtype: Tuple[ndarray, float]

1. Perform erosion for the pore space with specified radius

2. Label the eroded pore space, and leave only the pore space that is still connected with the non-wetting phase.

3. Perform the dilation for the labelled pore space with specified radius

This method was adapted from LBPM (OPM/LBPM)

Parameters:

• image – Binary 3D image with values 0 and 1 representing grain space and pore space, respectively

• radius – radius of the structuring element used for erosion and dilation

Returns:

A 3D image with estimated multiphase configuration and the wetting saturation at the specified invasion radius.

Return type:

Tuple[np.ndarray, float]