Skip to main content
Log in

Estimating Organic-Rich Shale Fractal Dimensions from Gas Adsorption Isotherms: Combining Different Methods Leads to More Reliable Values and Insight

  • Original Paper
  • Published:
Natural Resources Research Aims and scope Submit manuscript

Abstract

Fabric and roughness of the pore-size distributions in organic-rich shales determine their fluid flow and storage capabilities. Accurate estimation of the pore-surface fractal dimension (D) provides valuable insight to these qualities in shales. Low-pressure gas adsorption isotherms are widely used for determining D, typically applying the Frenkel–Halsey–Hill (FHH) method. Other D estimation methods, proposed by Neimark (NM) and Wang and Li (WL), are theoretically consistent and mathematically related to the FHH model but yield distinctive D values for many shales. This study evaluates the mathematical relationships between the FHH, NM and WL fractal determination methods, and with the aid of twenty-six published adsorption isotherms from shales around the world, compares their similarities and differences. Uncertainties exist in establishing best-fit lines to curved data trends in the FHH and NM methods, and in fitting power curves to data trends in the NM and WL methods. The FHH and WL D values are found to be more consistent for whole isotherm and isotherm segment analysis than the NM D values, which are systematically higher. The reasons for this are explained in terms of their graphical relationships. This leads to a novel 10-step protocol for a more thorough determination of shale D values that incorporates all three methods and involves graphical analysis that clearly exposes the uncertainties associated with the values determined. Applying this protocol should derive reliable D values to compare with key shale properties such as surface area, surface volume, thermal maturity and organic richness in future research.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7

Similar content being viewed by others

References

  • Alturki, A. A., Maini, B. B., & Gates, I. D. (2014). The effect of wall roughness on two-phase flow in a rough-walled Hele-Shaw cell. Journal of Petroleum Exploration and Production Technology, 4(4), 397–426.

    Article  Google Scholar 

  • Avnir, D., Farin, D., & Pfeifer, P. (1984). Molecular fractal surfaces. Nature, 308(5956), 261–263.

    Article  Google Scholar 

  • Avnir, D., Farin, D., & Pfeifer, P. (1992). A discussion of some aspects of surface fractality and of its determination. New Journal of Chemistry, 1, 439–449.

    Google Scholar 

  • Cai, J., Zhang, L., & Wei, W. (Eds.), (2020). Modelling of flow and transport in fractal porous media, Elsevier, 255 https://doi.org/10.1016/C2018-0-02631-6

  • Chalmers, G. R. L., & Bustin, R. M. (2007). The organic matter distribution and methane capacity of the Lower Cretaceous strata of Northeastern British Columbia, Canada. International Journal of Coal Geology, 70, 223–239.

    Article  Google Scholar 

  • Clarkson, C. R., Solano, N., Bustin, R. M., Bustin, A. M. M., Chalmers, G. R. L., He, L., & Blach, T. P. (2013). Pore structure characterization of North American shale gas reservoirs using USANS/SANS, gas adsorption, and mercury intrusion. Fuel, 103, 606–616.

    Article  Google Scholar 

  • Gao, Z., Fan, Y., Xuan, Q., & Zheng, G. (2020). A review of shale pore structure evolution characteristics with increasing thermal maturities. Advances in Geo-Energy Research, 4(3), 247–259.

    Article  Google Scholar 

  • Gregg, S. J., & Sing, K. S. W. (1982). Adsorption, surface area, and porosity (2nd ed.). Academic Press.

    Google Scholar 

  • Gu, Y., Ding, W., Yin, M., Jiao, B., Shi, S., Li, A., Xiao, Z., & Wang, Z. (2019). Nanoscale pore characteristics and fractal characteristics of organic-rich shale: An example from the lower Cambrian Niutitang Formation in the Fenggang block in northern Guizhou Province, South China. Energy Exploration & Exploitation, 37(1), 273–295.

    Article  Google Scholar 

  • Halsey, G. (1948). Physical adsorption on non-uniform surfaces. The Journal of Chemical Physics, 16, 931.

    Article  Google Scholar 

  • Hazra, B., Wood, D. A., Kumar, S., Saha, S., Dutta, S., Kumari, P., & Singh, A. K. (2018). Fractal disposition, porosity characterization and relationships to thermal maturity for the Lower Permian Raniganj basin shales, India. Journal of Natural Gas Science and Engineering, 59, 452–465.

    Article  Google Scholar 

  • Hill, T. L. (1952). Theory of physical adsorption. Advanced catalysis IV (pp. 211–257).

  • Jaroniec, M. (1995). Evaluation of the fractal dimension from a single adsorption isotherm. Langmuir, 11, 2316.

    Article  Google Scholar 

  • Kiselev, A. V., & Pavlova, L. F. (1965). Use of general equations of isotherms for the adsorption of benzene-n-hexane solutions on the adsorbents of various characters. Russian Chemical Bulletin, 14, 15–23.

    Article  Google Scholar 

  • Laaksonen, A., Malila, J., Nenes, A., Hung, H.-M., & Chen, J.-P. (2016). Surface fractal dimension, water adsorption efficiency, and cloud nucleation activity of insoluble aerosol. Science and Reports, 6, 25504.

    Article  Google Scholar 

  • Li, A., Ding, W., He, J., Dai, P., Yin, S., & Xie, F. (2016). Investigation of pore structure and fractal characteristics of organic-rich shale reservoir: A case study of Lower Cambrian Qiongzhusi formation in Malong block of eastern Yunnan Province, South China. Marine and Petroleum Geology, 70, 46–57.

    Article  Google Scholar 

  • Liu, K., Ostadhassan, M., Jang, H. W., Zakharova, N. V., & Shokouhimehr, M. (2021). Comparison of fractal dimensions from nitrogen adsorption data in shale via different models. RSC Advances, 11, 2298–2306.

    Article  Google Scholar 

  • Lu, Y., Zhou, J., Li, H., Chen, X., & Tang, J. (2020). Different effect mechanisms of supercritical CO2 on the shale microscopic structure. ACS Omega, 5, 22568–22577.

    Article  Google Scholar 

  • Mahamud, M. M., & Novo, M. F. (2008). The use of fractal analysis in the textural characterization of coals. Fuel, 87, 222–231.

    Article  Google Scholar 

  • Mandelbrot, B. B. (1975). Les objects fractals: Forme, hasard et dimension. Flammarion.

    Google Scholar 

  • Neimark, A. (1992). New approach to the determination of the surface fractal dimension of porous solids. Physica A: Statistical Mechanics and Its Applications, 191(1–4), 258–262.

    Article  Google Scholar 

  • Neimark, A., & Unger, K. K. (1993). Method of discrimination of surface fractality. Journal of Colloid and Interface Science, 158, 412–419.

    Article  Google Scholar 

  • Pfeifer, P. (1983). Chemistry in noninteger dimensions between two and three. I. Fractal theory of heterogeneous surfaces. The Journal of Chemical Physics, 79, 3558. https://doi.org/10.1063/1.446210

    Article  Google Scholar 

  • Pfeifer, P. (1984). Fractal dimension as working tool for surface roughness problems. Applied Surface Science, 18(1–2), 146–164.

    Article  Google Scholar 

  • Pfeifer, P., & Cole, M. W. (1990). Fractals in surface science: Scattering and thermodynamics of adsorbed films II. New Journal of Chemistry, 14, 221.

    Google Scholar 

  • Pfeifer, P., Wu, Y. J., Cole, M. W., & Krim, J. (1989). Multilayer adsorption on a fractally rough surface. Physical Review Letters, 62(17), 1997.

    Article  Google Scholar 

  • Powles, J. G. (1985). On the validity of the Kelvin equation. Journal of Physics A: Mathematical and General, 18(9), 1551.

    Article  Google Scholar 

  • Rigby, S. P., Jahan, H., Stevens, L., Clement Uguna, C., Snape, C., Macnaughton, B., Large, D. J., & Fletcher, R. S. (2020). Pore structural evolution of shale following thermochemical treatment. Marine and Petroleum Geology, 112, 104058.

    Article  Google Scholar 

  • Sahouli, B., Blacher, S., & Brouers, F. (1996). Fractal surface analysis by using nitrogen adsorption data: The case of the capillary condensation regime. Langmuir, 12(11), 2872–2874.

    Article  Google Scholar 

  • Sakhaee-Pour, A., & Li, W. (2016). Fractal dimensions of shale. Journal of Natural Gas Science and Engineering, 30, 578–582.

    Article  Google Scholar 

  • Shao, X., Pang, X., Li, Q., Wang, P., Chen, D., Shen, W., & Zhao, Z. (2017). Pore structure and fractal characteristics of organic-rich shales: A case study of the lower Silurian Longmaxi shales in the Sichuan Basin, SW China. Marine and Petroleum Geology, 80, 192–202.

    Article  Google Scholar 

  • Sing, K. S. W., Everett, D. H., Haul, R. A. W., Moscou, L., Pierotti, R. A., Rouquerol, J., Rouquerol, F., & Siemieniewskat, T. (1985). Reporting physisorption data for gas/solid systems with special reference to the determination of surface area and porosity (Recommendations 1984). Pure and Applied Chemistry, 57, 603–619.

    Article  Google Scholar 

  • Sinha, A.S. (2017). Surface area study in organic-rich shales using nitrogen adsorption. MSc Thesis University of Oklahoma (USA) 127 pages. https://shareok.org/handle/11244/52243 [Accessed 11 Mar 2021]

  • Sun, M., Yu, B., Hu, Q., Chen, S., Xia, W., & Ye, R. (2016). Nanoscale pore characteristics of the Lower Cambrian Niutitang Formation Shale: A case study from Well Yuke #1 in the Southeast of Chongqing, China. International Journal of Coal Geology, 154–155, 16–29.

    Article  Google Scholar 

  • Tian, Z., Wei, W., Zhou, S., Wood, D. A., & Cai, J. (2021). Experimental and fractal characterization of the microstructure of shales from Sichuan Basin. China. Energy & Fuels, 35(5), 3899–3914.

    Article  Google Scholar 

  • Wang, F., & Li, S. (1997). Determination of the surface fractal dimension for porous media by capillary condensation. Industrial and Engineering Chemistry Research, 36(5), 1598–1602.

    Article  Google Scholar 

  • Wang, Y., Du, B., Dou, X., Liu, J., Shi, B., Wang, D., & Tang, H. (2007). Study on the pore surface fractal dimension and surface acid–base properties of natural particles around Guanting reservoir. Colloids and Surfaces A: Physicochemical and Engineering Aspects, 307, 16–27.

    Article  Google Scholar 

  • Wood, D. A., & Hazra, B. (2017). Characterization of organic-rich shales for petroleum exploration & exploitation: A review- part 1: Bulk properties, multi-scale geometry and gas adsorption. Journal of Earth Science, 28(5), 739–757.

    Article  Google Scholar 

  • Yang, F., Ning, Z., & Liu, H. (2014). Fractal characteristics of shales from a shale gas reservoir in the Sichuan Basin, China. Fuel, 115, 378–384.

    Article  Google Scholar 

  • Yang, R., He, S., Yi, J., & Hu, Q. (2016). Nano-scale pore structure and fractal dimension of organic-rich Wufeng-Longmaxi shale from Jiaoshiba area, Sichuan Basin: Investigations using FE-SEM, gas adsorption and helium pycnometry. Marine and Petroleum Geology, 70, 27–45.

    Article  Google Scholar 

  • Zhao, H., Lai, Z., & Firoozabadi, A. (2017). Sorption hysteresis of light hydrocarbons and carbon dioxide in shale and kerogen. Scientific Reports, 7, 16209.

    Article  Google Scholar 

Download references

Acknowledgements

Many thanks are due to Kouqi Liu for clarifying aspects of the calculation methods used in Liu et al. (2021) and Jianchao Cai and his associates (Tian et al., 2021) for providing some useful additional insights with respects of fractal dimensions and their analysis.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to David A. Wood.

Appendix: details of 26 published shale isotherms evaluated

Appendix: details of 26 published shale isotherms evaluated

The organic-rich shale adsorption isotherms evaluated with the fractal-dimension calculation methods have been compiled from published data from shale basins distributed around the world. These adsorption isotherms are displayed in Figs. 8, 9, 10 and 11. Details of the samples from which these isotherms are derived are listed below and available in the cited publications.

Figure 8
figure 8

Adsorption isotherms displayed for samples #1 to #7. N refers to the volume of adsorbate adsorbed by the sample

Figure 9
figure 9

Adsorption isotherms displayed for samples #8 to #14. N refers to the volume of adsorbate adsorbed by the sample

Figure 10
figure 10

Adsorption isotherms displayed for samples #15 to #21. N refers to the volume of adsorbate adsorbed by the sample

Figure 11
figure 11

Adsorption isotherms displayed for samples #22 to #26. N refers to the volume of adsorbate adsorbed by the sample

#1 Duvernay: shale (Canada) from Fig. 2a in Clarkson et al. (2013)

#2 Barnett: shale (USA) from Fig. 2a in Clarkson et al. (2013)

#3 Haynesville: shale (USA) from Fig. 2a in Clarkson et al. (2013)

#4 Woodford: shale (USA) from Fig. 2a in Clarkson et al. (2013)

#5 Eagle Ford: shale (USA) from Fig. 2a in Clarkson et al. (2013)

#6 Mont2 (Milk River): shale (Canada) from Fig. 2a in Clarkson et al. (2013)

#7 Muska1 (Horn River): shale (Canada) from Fig. 2a in Clarkson et al. (2013)

#8 Muska2 (Horn River): shale (Canada) from Fig. 2a in Clarkson et al. (2013)

#9 Marcellus 1: shale (USA) from Fig. 2a in Clarkson et al. (2013)

#10 Marcellus2: shale (USA) from Fig. 28d in Sinha (2017)

#11 Wolfcamp: shale (USA) (surface area 25 m2/g) from Fig. 28c in Sinha (2017)

#12 Bakken Immature: shale (USA) (surface area 16.44 m2/g) from Fig. 28a in Sinha (2017)

#13 Bakken1: sample #1 shale (USA) from Fig. 1 in Liu et al. (2021)

#14 Bakken 2: sample #16 (with low micropore volume) from Fig. 1 in Liu et al. (2021)

#15 Bakken 3: sample #18 (with high micropore volume) from Fig. 1 in Liu et al. (2021)

#16 Kimmeridge: from Fig. 1 (Kimmeridge Blackstone, UK) Zhao et al. (2017)

#17 Bowland (untreated): Upper Bowland Shale (UK) from Fig. 1 in Rigby et al. (2020)

#18 Bowland (treated): Thermally treated #17 samples from Fig. 1 in Rigby et al. (2020)

#19 Raniganj3: Sample CG1003 (Raniganj Basin, India) from Fig. 9 in Hazra et al. (2018)

#20 Raniganj13: Sample CG1013 (Raniganj Basin, India) from Fig. 10 in Hazra et al. (2018)

#21 Raniganj19: Sample CG1019 (Raniganj Basin, India) from Fig. 10 in Hazra et al. (2018)

#22 Longmaxi1: Sample W1-1 (Sichuan Basin, China) from Fig. 3a in Tian et al. (2021)

#23 Longmaxi5: Sample W1-5 (Sichuan Basin, China) from Fig. 3a in Tian et al. (2021)

#24 Fuling1: untreated sample from Fig. 12a Lu et al. (2020)

#25 Fuling2: sample treated with Argon at 4 MPa pressure from Fig. 12a Lu et al. (2020)

#26 Fuling3: sample treated with CO2 at 16 MPa pressure from Fig. 12d Lu et al. (2020)

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wood, D.A. Estimating Organic-Rich Shale Fractal Dimensions from Gas Adsorption Isotherms: Combining Different Methods Leads to More Reliable Values and Insight. Nat Resour Res 30, 3551–3574 (2021). https://doi.org/10.1007/s11053-021-09909-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s11053-021-09909-2

Keywords

Navigation