Automated Image Analysis of Mud and Mudrock Microstructure and Characteristics of Hemipelagic Sediments: IODP Expedition 339

Abstract

The microstructural analysis of muds and mudrocks requires very high-resolution measurement. Recent advances in electron microscopy have contributed significantly to the improved characterisation of mudrock microstructures and their consequent petrophysical properties. However, imaging through electron microscopy is limited to small areas of coverage such that upscaling of these properties is a great challenge. In this paper, we develop a new methodology for multiple large-area imaging using scanning electron microscopy through automated acquisition and stitching from polished thin-sections and ion-milled samples. The process is fast, efficient and minimises user-input and bias. It can provide reliable, quantifiable data on sediment grain size, grain orientation, pore size and porosity. Limitations include the time involved for individual runs and manual segmentation, the large amount of computer memory required, and instrument resolution at the nano-scale. This method is applied to selected samples of Quaternary muddy sediments from the Iberian margin at IODP Site 1385. The section comprises finegrained (very fine clayey silts), mixed-composition, biogenic-terrigenous hemipelagites, with a pronounced but non-regular colour cyclicity. There is a multi-tiered and diverse trace fossil assemblage of the deep-water Zoophycos ichnofacies. The sediment microstructures show small-scale heterogeneity in all properties, and an overall random fabric with secondary preferred grain-alignment. These results on the fabric differ, in part, from previous studies of hemipelagic muds. Further work is underway on their comparison with other deep-water sediment facies.

Keywords

mudrocks / microstructure / microporosity / grain-orientation / hemipelagites / trace fossils

References

[1] Aplin A C, Macquaker J H S. Mudstone Diversity: Origin and Implications for Source, Seal, and Reservoir Properties in Petroleum Systems. AAPG Bulletin, 2011, 95(12): 2031-2059.

[2] Bankole S A, Buckman J, Stow D, . Grain-Size Analysis of Mudrocks: A New Semi-Automated Method from SEM Images. Journal of Petroleum Science and Engineering, 2019, 174: 244-256.

[3] Bankole S A, Stow D A V, Lever H, . Microstructure of Mudrock and the Choice of Representative Sample. In: Fifth EAGE Shale Workshop, 2016, Catania, Italy: EAGE

[4] Berens P. CircStat: AMATLAB Toolbox for Circular Statistics. Journal of Statistical Software, 2009, 31(10): 1-21.

[5] Bosl W J, Dvorkin J, Nur A. A Study of Porosity and Permeability Using a Lattice Boltzmann Simulation. Geophysical Research Letters, 1998, 25(9): 1475-1478.

[6] Buckman J. Use of Automated Image Acquisition and Stitching in Scanning Electron Microscopy: Imaging of Large Scale Areas of Materials at High Resolution. Microscopy and Analysis, 2014, 28: 13-15.

[7] Buckman J, Bankole S, Zihms S, . Quantifying Porosity through Automated Image Collection and Batch Image Processing: Case Study of Three Carbonates and an Aragonite Cemented Sandstone. Geosciences, 2017, 7 3 70

[8] Camp W K, Diaz E, Wawak B E. Electron Microscopy of Shale Hydrocarbon Reservoirs, 2013, Tulsa: Association of Petroleum Geologists

[9] Chambers J M, Cleveland W S, Kleiner B, . Graphical Methods for Data Analysis. Journal of the Royal Statistical Society, 1984, 147 3 513.

[10] Curtis M E, Sondergeld C H, Ambrose R J, . Microstructural Investigation of Gas Shales in Two and Three Dimensions Using Nanometer-Scale Resolution Imaging. AAPG Bulletin, 2012, 96(4): 665-677.

[11] Curtis, M. E., Ambrose, R. J., Sondergeld, C. H., et al., 2010. Structural Characterization of Gas Shales on the Micro- and Nano-Scales, Society of Petroleum Engineers—Canadian Unconventional Resources and International Petroleum Conference, Calgary

[12] Davis J C. Statistics and Data Analysis in Geology, 1986, India: Wiley, 656.

[13] Desbois G, Urai J L, Kukla P A. Morphology of the Pore Space in Claystones—Evidence from BIB/FIB Ion Beam Sectioning and Cryo-SEM Observations. eEarth, 2009, 4(1): 15-22.

[14] DeVasto M A, Czeck D M, Bhattacharyya P. Using Image Analysis and ArcGIS® to Improve Automatic Grain Boundary Detection and Quantify Geological Images. Computers & Geosciences, 2012, 49: 38-45.

[15] Fisher N I. Statistical Analysis of Circular Data, 1993, New York: Cambridge University Press, 277

[16] Francus P, Pirard E. Francus P. Testing for Sources of Errors in Quantitative Image Analysis. Image Analysis, Sediments and Paleoenvironments, 2004, Netherlands, Dordrecht: Springer, 87-102.

[17] Grove C, Jerram D A. JPOR: An ImageJ Macro to Quantify Total Optical Porosity from Blue-Stained Thin Sections. Computers & Geosciences, 2011, 37(11): 1850-1859.

[18] Hemes S, Desbois G, Urai J L, . Variations in the Morphology of Porosity in the Boom Clay Formation: Insights from 2D High Resolution BIB-SEM Imaging and Mercury Injection Porosimetry. Netherlands Journal of Geosciences, 2013, 92(4): 275-300.

[19] Hesse R. Turbiditic and Non-Turbiditic Mudstone of Cretaceous Flysch Sections of the East Alps and other Basins. Sedimentology, 1975, 22(3): 387-416.

[20] Hodell D A, Lourens L, Stow D A V, . The “Shackleton Site” (IODP Site U1385) on the Iberian Margin. Scientific Drilling, 2013, 16: 13-19.

[21] Hodell D A, Lourens L, Crowhurst S, . A Reference Time Scale for Site U1385 (Shackleton Site) on the SW Iberian Margin. Global and Planetary Change, 2015, 133: 49-64.

[22] Hoogakker B A A, Rothwell R G, Rohling E J, . Variations in Terrigenous Dilution in Western Mediterranean Sea Pelagic Sediments in Response to Climate Change during the Last Glacial Cycle. Marine Geology, 2004, 211(1/2): 21-43.

[23] Houben M E, Desbois G, Urai J L. Pore Morphology and Distribution in the Shaly Facies of Opalinus Clay (Mont Terri, Switzerland): Insights from Representative 2D BIB-SEM Investigations on mm to nm Scale. Applied Clay Science, 2013, 71: 82-97.

[24] Janssen C, Kanitpanyacharoen W, Wenk H R, . Clay Fabrics in SAFOD Core Samples. Journal of Structural Geology, 2012, 43: 118-127.

[25] Josh M, Esteban L, Delle Piane C, . Laboratory Characterisation of Shale Properties. Journal of Petroleum Science and Engineering, 2012, 88–89: 107-124.

[26] Kameda A, Dvorkin J, Keehm Y, . Permeability-Porosity Transforms from Small Sandstone Fragments. Geophysics, 2006, 71(1): N11-N19.

[27] Keller L M, Schuetz P, Erni R, . Characterization of Multi-Scale Microstructural Features in Opalinus Clay. Microporous and Mesoporous Materials, 2013, 170: 83-94.

[28] Kuila U, Prasad M. Specific Surface Area and Pore-Size Distribution in Clays and Shales. Geophysical Prospecting, 2013, 61(2): 341-362.

[29] Lemmens, H., Richards, D., 2013. Multiscale Imaging of Shale Samples in the Scanning Electron Microscope. American Association of Petroleum Geologists Memoir, 27–35. https://doi.org/10.1306/13391702M1023582

[30] Lonardelli I, Wenk H R, Ren Y. Preferred Orientation and Elastic Anisotropy in Shales. Geophysics, 2007, 72(2): D33-D40.

[31] Loucks R G, Reed R M, Ruppel S C, . Morphology, Genesis, and Distribution of Nanometer-Scale Pores in Siliceous Mudstones of the Mississippian Barnett Shale. Journal of Sedimentary Research, 2009, 79(12): 848-861.

[32] Loucks R G, Reed R M, Ruppel S C, . Spectrum of Pore Types and Networks in Mudrocks and a Descriptive Classification for Matrix-Related Mudrock Pores. AAPG Bulletin, 2012, 96(6): 1071-1098.

[33] Lovie S. Empirical Quantile-Quantile Plots. In: Encyclopedia of Statistics in Behavioral Science, 2005, Chichester: John Wiley and Sons, Ltd., 543-545.

[34] Macquaker J H S, Howell J K. Small-Scale (<5.0 M) Vertical Heterogeneity in Mudstones: Implications for High-Resolution Stratigraphy in Siliciclastic Mudstone Successions. Journal of the Geological Society, 1999, 156(1): 105-112.

[35] Mardia, K. V., Jupp, P. E., 2008. Directional Statistics, Directional Statistics. John Wiley and Sons, Inc. 432. https://doi.org/10.1002/9780470316979

[36] Martínez-Nistal A, Veniale F, Setti M, . A Scanning Electron Microscopy Image Processing Method for Quantifying Fabric Orientation of Clay Geomaterials. Applied Clay Science, 1999, 14(4): 235-243.

[37] Moon C F, Hurst C W. Fabric of Muds and Shales: An Overview. Geological Society, London, Special Publications, 1984, 15(1): 579-593.

[38] Munson E O, Chalmers G R L, Bustin R M, . Utilizing Smear Mounts for X-Ray Diffraction as a Fully Quantitative Approach in Rapidly Characterizing the Mineralogy of Shale Gas Reservoirs. Journal of Unconventional Oil and Gas Resources, 2016, 14: 22-31.

[39] Nishida N. Microstructure of Muddy Contourites from the Gulf of Cádiz. Marine Geology, 2016, 377: 110-117.

[40] Nishida N, Ito M, Inoue A, . Clay Fabric of Fluid-Mud Deposits from Laboratory and Field Observations: Potential Application to the Stratigraphic Record. Marine Geology, 2013, 337: 1-8.

[41] Pal N R, Pal S K. A Review on Image Segmentation Techniques. Pattern Recognition, 1993, 26(9): 1277-1294.

[42] Pickering K T, Hiscott R N. Deep Marine Systems: Processes, Deposits, Environments, Tectonics and Sedimentation, 2015, Chichester: John Wiley and Sons, 657.

[43] Piper D J W. Manual of Sedimentological Techniques. Departments of Geology and Oceanography, 1977, Halifax: Dalhousie University

[44] Rodríguez-Tovar F J, Dorador J. Ichnological Analysis of Pleistocene Sediments from the IODP Site U1385 “Shackleton Site” on the Iberian Margin: Approaching Paleoenvironmental Conditions. Palaeogeography, Palaeoclimatology, Palaeoecology, 2014, 409: 24-32.

[45] Rodríguez-Tovar F J, Dorador J, Grunert P, . Deep-Sea Trace Fossil and Benthic Foraminiferal Assemblages across Glacial Terminations 1, 2 and 4 at the “Shackleton Site” (IODP Expedition 339, Site U1385). Global and Planetary Change, 2015, 133: 359-370.

[46] Saraji S, Piri M. The Representative Sample Size in Shale Oil Rocks and Nano-Scale Characterization of Transport Properties. International Journal of Coal Geology, 2015, 146: 42-54.

[47] Schindelin J, Arganda-Carreras I, Frise E, . Fiji: An Open-Source Platform for Biological-Image Analysis. Nature Methods, 2012, 9: 676-682.

[48] Schindelin J, Rueden C T, Hiner M C, . The ImageJ Ecosystem: An Open Platform for Biomedical Image Analysis. Molecular Reproduction and Development, 2015, 82(7/8): 518-529.

[49] Schneider C A, Rasband W S, Eliceiri K W. NIH Image to ImageJ: 25 Years of Image Analysis. Nature Methods, 2012, 9(7): 671-675.

[50] Sing K S W, Everett D H, Haul R A W, . Reporting Physisorption Data for Gas Solid Systems with Special Reference to the Determination of Surface-Area and Porosity (Recommendations 1984). Pure Applied Chemistry, 1985, 57: 603-619.

[51] Sokolov V N, O’Brien N R. A Fabric Classification of Argillaceous Rocks, Sediments, Soils. Applied Clay Science, 1990, 5(4): 353-360.

[52] Stow D A V. Fine-Grained Sediments in Deep Water: An Overview of Processes and Facies Models. Geo-Marine Letters, 1985, 5(1): 17-23.

[53] Stow D A V, Hernández-Molina F J, Llave E, . The Cadiz Contourite Channel: Sandy Contourites, Bedforms and Dynamic Current Interaction. Marine Geology, 2013, 343: 99-114.

[54] Stow D A V, Huc A Y, Bertrand P. Depositional Processes of Black Shales in Deep Water. Marine and Petroleum Geology, 2001, 18(4): 491-498.

[55] Stow D A V, Tabrez A R. Hemipelagites: Processes, Facies and Model. Geological Society, London, Special Publications, 1998, 129(1): 317-337.

[56] Stow D A V. Sedimentary Rocks in the Field: A Color Guide, 2005, Florida: Taylor and Francis Group, 320

[57] Suttle M D, Genge M J, Russell S S. Shock Fabrics in Fine-Grained Micrometeorites. Meteoritics & Planetary Science, 2017, 52(10): 2258-2274.

[58] Tovey N K, Smart P, Hounslow M W, . Automatic Orientation Mapping of Some Types of Soil Fabric. Geoderma, 1992, 53(3/4): 179-200.

[59] Uchman A, Wetzel A. Deep-Sea Ichnology: The Relationships between Depositional Environment and Endobenthic Organisms. Developments in Sedimentology, 2011, 63: 517-556.

[60] Vanden-Bygaart A J, Protz R. The Representative Elementary Area (REA) in Studies of Quantitative Soil Micromorphology. Geoderma, 1999, 89(3/4): 333-346.

[61] Wang Y, Zhu Y M, . Characteristics of the Nanoscale Pore Structure in Northwestern Hunan Shale Gas Reservoirs Using Field Emission Scanning Electron Microscopy, High-Pressure Mercury Intrusion, and Gas Adsorption. Energy & Fuels, 2014, 28(2): 945-955.

[62] Wenk H R, Houtte P V. Texture and Anisotropy. Reports on Progress in Physics, 2004, 67(8): 1367-1428.

[63] Wenk H R, Lutterotti L, Kaercher P, . Rietveld Texture Analysis from Synchrotron Diffraction Images. II. Complex Multiphase Materials and Diamond Anvil Cell Experiments. Powder Diffraction, 2014, 29(3): 220-232.

[64] Wenk H R, Voltolini M, Martin M, . Preferred Orientations and Anisotropy in Shales: Callovo-Oxfordian Shale (France) and Opalinus Clay (Switzerland). Clays and Clay Minerals, 2008, 56(3): 285-306.

[65] Yang F, Ning Z F, Liu H Q. Fractal Characteristics of Shales from a Shale Gas Reservoir in the Sichuan Basin, China. Fuel, 2014, 115: 378-384.

[66] Zaitoun N M, Aqel M J. Survey on Image Segmentation Techniques. Procedia Computer Science, 2015, 65: 797-806.