pc_curve#
- pc_curve(im, sizes=None, pc=None, seq=None, sigma=0.072, theta=180, voxel_size=1)[source]#
Produces a Pc-Snwp curve given a map of meniscus radii or capillary pressures at which each voxel was invaded
- Parameters:
im (ndarray) – The voxel image of the porous media with
True
values indicating the void spacesizes (ndarray, optional) – An image containing the sphere radii at which each voxel was invaded during an invasion experiment.
pc (ndarray, optional) – An image containing the capillary pressures at which each voxel was invaded during an invasion experiment.
seq (ndarray, optional) – An image containing invasion sequence values, such as that returned from the
ibip
function.sigma (float, optional) – The surface tension of the fluid-fluid system of interest. This argument is ignored if
pc
are specified, otherwise it is used in the Washburn equation to convertsizes
to capillary pc.theta (float) – The contact angle measured through the invading phase in degrees. This argument is ignored if
pc
are specified, otherwise it is used in the Washburn equation to convertsizes
to capillary pressures.voxel_size (float) – The voxel resolution of the image. This argument is ignored if
pc
are specified, otherwise it is used in the Washburn equation to convertsizes
to capillary pressures.
- Returns:
pc_curve – A custom object with the following data added as named attributes:
Attribute
Description
pc
The capillary pressure, either as given in
pc
or computed fromsizes
(see Notes).snwp
The fraction of void space filled by non-wetting phase at each pressure in
pc
- Return type:
Results object
Notes
If
sizes
is provided, then the Washburn equation is used to convert the radii to capillary pressures, using the givensigma
andtheta
values, along with thevoxel_size
if the values are in voxel radii. For more control over how capillary pressure model, it can be computed by hand, for example:$$ pc = frac{-2*0.072*np.cos(np.deg2rad(180))}{sizes cdot voxel_size} $$
then passed in as the
pc
argument.Examples
Click here to view online example.