Automatic Measurements Of Contact Angle And Curvature In Pore-Space Images

Determining how carbon dioxide (CO2) can be stored securely and safely in underground aquifers, oil and gas flow rate and overall recovery from hydrocarbon reservoirs and shales, the efficiency of gas exchange in leaves in plants and human lungs, the performance of catalysts and fuel cells, how well fabrics resist or soak up water, and the design of water-repellent surfaces, all crucially depend on wettability: how fluid phases interact with solid surfaces within porous structures.

Conventionally, wettability can be inferred through indirect averaged measurements on cm-sized rock samples in the laboratory, such as the Amott and USBM indices, or by studying relative permeability characteristics. At the pore scale, wettability is defined in terms of the contact angle between the fluids residing in the pore space. Moreover, laboratory contact angle measurements are made using smooth mineral surfaces, typically quartz and calcite to represent sandstones and carbonates, respectively.

While wettability has been characterized for external surfaces, we use micron-resolution X-ray three-dimensional imaging to look inside complex porous structures to measure the contact angle and fluid/fluid interface curvature.

In this study, we have developed a novel automated algorithm to measure in situ contact angles and fluid/fluid interface curvature from three-dimensional pore-scale images. The principal idea is to discretize the fluid/fluid and fluid/solid interfaces and define vectors that have a direction perpendicular to these surfaces. The contact angle is found from the dot product of vectors describing the fluid/fluid interface and the solid (rock in this study) surface, where they meet at the three-phase contact line.

We also apply the physical insight that — in equilibrium — the fluid/fluid interfaces should have a uniform curvature. These surfaces are smoothed to increase the accuracy and reliability of the measurements: we apply a volume preserving Gaussian smoothing on all vertices belonging to the generated mesh of the contacted interfaces, which removes the voxelized artifacts from the segmented 3-D image; then, a local uniform curvature smoothing is applied only on the fluid/fluid (oil/brine in our study) interface.

Figure 1. A schematic of wettability and length scales in porous media flows. Right shows images of
oil retained within the pore space of water-wet and oil-wet (or water-repellent) samples. The curvature of the
oil/brine interface is indicated by the colours, while the solid surface in contact with brine is shown in grey. In
a water-wet rock, the oil is trapped in quasi-spherical ganglia with contact angles less than 90◦. For an oil-wet
rock, oil is retained in layers that follow the surface roughness with both positive and negative values of the
curvature κ. This allows the oil to flow to low saturation, facilitating recovery at the km scale (Left). Center
shows how the samples are selected from centimeter-scale cores extracted from a reservoir. Figure republished from PNAS Open Access article:

This novel method has been successfully tested on benchmark synthetic images for which the contact angle and curvature are known analytically and has been applied to a real rock data set. Both a water-wet and a mixed-wet rock have been analyzed. Further, the method has been utilized to study the spatial correlation of the in situ measured distributions of contact angles and oil/brine interface curvatures within reservoir rock samples at subsurface conditions, which originally acquired from a giant oil-producing reservoir in the Middle East.

Measuring in situ contact angles in different types of core samples would allow researchers to define more realistic contact angle distributions, which include natural variability and uncertainty, for reservoir characterization. These distributions could then be used to more accurately predict multiphase fluid flow in pore-scale models. Altered wettability systems that have undergone fluid displacement cycles are likely to have a wide range of contact angles and should not be considered to be uniformly wetting or non-wetting. Furthermore, the curvatures estimated from this algorithm could be used to estimate local capillary pressures.

Figure 2. Contact angle and oil/brine interface curvature measurements. (A) Three-dimensional view
of the raw segmented dataset of sample-1 with a voxel size of 2 × 2 × 2 microns. We cropped out a sub-volume from sample-1 in (B) to illustrate the following automatic contact angle measurement workflow. (C) The oil/brine (green) and oil/solid (red) interfaces are identified. (D) The smoothed mesh M that is found after applying both Gaussian and curvature smoothing. (E) The measured oil/brine interface curvature values of all vertices belonging to the oil/brine interface and contact line (i ∈ {VO B , VC L }). (F) The three-phase contact line. (G) Normal vectors are defined on both the oil/brine and brine/rock interfaces for the contact line set, i ∈ VC L. The contact angle is calculated from the dot product of these two normals. Republished from Wiley Open Access article:

*The codes for measuring contact angle, curvature and roughness are available at:

This work was conducted by Ahmed AlRatrout, Ali Q Raeini, Branko Bijeljic, and Martin J Blunt from Imperial College London.


These findings are described in the articles entitled Automatic measurement of contact angle in pore-space images, recently published in the journal Advances in Water Resources and Spatial Correlation of Contact Angle and Curvature in Pore‐Space Images, recently published in the journal Water Resources Research.

AlRatrout, A., Blunt, M., & Bijeljic, B. (2018b). Wettability in complex porous materials, the mixed-wet state, and its relationship to surface roughness. Proceedings of the National Academy of Sciences, 115 (36), 8901-8906. (USA) doi: 10.1073/pnas.1803734115

About The Author

Ahmed AlRatrout

Ahmed A. AlRatrout is from Imperial College London.

Speak Your Mind!


Meta-Sulfamoyl N-Hydroxybenzamides As HDAC8-Selective Inhibitors

As a family of epigenetic enzymes, histone deacetylases (HDACs) can catalyze the removal of an acetyl group from N-acetyl lysine residues of histones and other proteins. The human HDACs can be grouped into 4 classes based on their homology to yeast prototypes: Class I (HDAC1, 2, 3 and 8), Class II (HDAC 4,5, 6, 7, […]

The Telomerase Gene Therapy For The Treatment Of Age-Related Diseases And Telomere Syndromes Does Not Favor Cancer Development

Telomeres are protective structures at the ends of chromosomes. They consist of tandem repeats of a DNA sequence, which in all vertebrates is TTAGGG, bound by a six-protein complex known as shelterin that exerts the protective function. Telomeres shorten every time a cell divides, and when they become critically short, this induces a persistent DNA […]

“Schrödinger’s Immigrant,” Or Why Immigrants Can Be Both Hardworking And Lazy At The Same Time

We all have heard the argument that immigrants take jobs away from people born in the country. But we also all have heard the argument they are a drain on a country’s resources by taking advantage of social security benefits. As some have pointed out, there seems to be a contradiction between these two arguments. […]

Homo Sapiens May Have Inherited Defenses Against RNA Viruses From Neanderthals

Though it was once thought that the evolution of the human species was a linear progression from ape-like ancestor to modern humans, recent findings have shown that extinct subspecies of early humans existed and that the line from our ancestors to us was substantially more complex. In fact, our Homo sapien ancestors are believed to […]

What Do Catfish Eat?

Catfish eat a wide range of plants and animals but primarily eat insects, water beetles, small fish, plants, algae, crawfish, snails, and just about anything else they can manage to fit in their mouths. Most catfish are not picky eaters and are scavengers, meaning they are happy to eat just about anything. Our diet is […]

Autotroph Vs Heterotroph: The Difference

The difference between an autotroph vs. heterotroph lies in the organism’s ability either produce its own food or need to ingest food. Autotrophs, such as plants, can produce their own food from light via photosynthesis or chemicals via chemosynthesis. Heterotrophs, such as humans, require ingesting food in order to produce energy and thus are not […]

How To Reverse The Unethical Consequences Of Sleep Deprivation

Published by David Welsh Arizona State University, W. P. Carey School of Business These findings are described in the article entitled Overcoming the effects of sleep deprivation on unethical behavior: An extension of integrated self-control theory, recently published in the Journal of Experimental Social Psychology (Journal of Experimental Social Psychology 76 (2018) 142-154). This work […]