File Exchange

image thumbnail

2D/3D specific surface of porous media using 2D images

version 2.0.0.0 (11.2 KB) by Arash Rabbani
This code calculates 2D/3D specific surface of porous media using 2D images

7 Downloads

Updated 09 Feb 2019

View Version History

View License

2D/3D specific surface of porous media using 2D images

Specific surface of the porous media has two different definitions:
1- 3-D Specific surface of the void/porous space, which is equal to the surface
are of pores divided by their volume. In 2-D the definition is: perimeter of pores
divided by their area.
2- 3-D Specific surface of the solid/grain space, which is equal to the surface
are of grains divided by their volume. In 2-D the definition is: perimeter of grains
divided by their area.
This code computes all of these parameters and visualizes the grain/pore boundaries.
The unit is: micron ^-1. The resolution of the image is set to be 5 micron per pixel.
Consider that you can change the input resolution in beginning lines of the code.
the Rabbani et al. (2014) have founded that the 3-D specific surface is approximately 1.35 times
greater than its 2-D value. In addition this code computes the porosity of the porous media.
If porosity obtained irregular you should invert the color of your input image.
If you want to use this code in your researches you can cite these articles:

1- Rabbani, A., Jamshidi, S., & Salehi, S. (2014). Determination of specific surface of
rock grains by 2D imaging. Journal of Geological Research, 2014.

2- Rabbani, A., & Jamshidi, S. (2014). Specific surface and porosity relationship for
sandstones for prediction of permeability. International Journal of Rock Mechanics and
Mining Sciences, 71, 25-32.

By Arash Rabbani
arashrabbani.com
rabarash@yahoo.com

Cite As

1- Rabbani, A., Jamshidi, S., & Salehi, S. (2014). Determination of specific surface of rock grains by 2D imaging. Journal of Geological Research, 2014. 2- Rabbani, A., & Jamshidi, S. (2014). Specific surface and porosity relationship for sandstones for prediction of permeability. International Journal of Rock Mechanics and Mining Sciences, 71, 25-32.

Comments and Ratings (6)

Arash Rabbani

Hi Maria
It supposed to be used like this
P=bwperim(A,Conn);
It is the pixel connectivity criteria when checking the perimeter pixels. You can check the MATLAB help for this.

Maria Huguett

Dear Arash Rabbani,
could you tell me please, what is the meaning of Conn=8 in the code? and why it takes the value 8?
Thank You

Ashane Fernando

Arash Rabbani

Yes, soon I will upload such code to find the throat sizes.

Muhammad RSMY

dear Arash Rabbani

could you please tell me if possible to calculate the pore-throat??

thanks

MATLAB Release Compatibility
Created with R2015a
Compatible with any release
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!