A comparison of clustering algorithms applied to fluid characterization using NMR T1-T2 maps of shale

Abstract Nuclear magnetic resonance T1-T2 maps are popular for characterizing fluids in shale. The characterization, however, is often done manually, which is difficult for shale due to its complicated nature. This work investigates the clustering approach for fluid characterization on T1-T2 maps by comparing 6 different algorithms: k-means, Gaussian mixture model, spectral clustering, and 3 hierarchical methods. T1-T2 maps are collected on shale samples at as-received and dried conditions. We propose two cluster validity indices to select the optimal cluster number and best algorithm. Our results validate the capability of the two indices. Gaussian mixture model is found to be the best algorithm in most of the cases, as its fluid partitioning shows the highest consistency with theoretical fluid boundaries. In addition, 5 fluid components are identified from Gaussian mixture model, and their values are qualitatively validated by comparing with those in literature. Results also indicate that clustering is sensitive to the fluid distribution. Drying the sample producing better clustering by revealing the footprint of organic matter. This work provides a practical guide for applying cluster analysis in fluid characterization in Nuclear magnetic resonance T1-T2 maps.

[1]  Hao Li,et al.  Prediction of Subsurface NMR T2 Distributions in a Shale Petroleum System Using Variational Autoencoder-Based Neural Networks , 2017, IEEE Geoscience and Remote Sensing Letters.

[2]  Hao Li,et al.  Long Short-Term Memory and Variational Autoencoder With Convolutional Neural Networks for Generating NMR T2 Distributions , 2019, IEEE Geoscience and Remote Sensing Letters.

[3]  Geoffrey E. Hinton,et al.  A View of the Em Algorithm that Justifies Incremental, Sparse, and other Variants , 1998, Learning in Graphical Models.

[4]  Peter Adams,et al.  The EMMIX software for the fitting of mixtures of normal and t-components , 1999 .

[5]  Marc Fleury,et al.  Characterization of shales using T1–T2 NMR maps , 2016 .

[6]  Carlos Torres-Verdín,et al.  New Methods for Interpretation of NMR Measurements Acquired in Spatially Heterogeneous Rocks , 2015 .

[7]  Mukul M. Sharma,et al.  Characterization of Free and Bound Fluids in Hydrocarbon Bearing Shales Using NMR and Py GC-MS , 2014 .

[8]  E. Purcell,et al.  Relaxation Effects in Nuclear Magnetic Resonance Absorption , 1948 .

[9]  Han Jiang,et al.  Fracture capture of organic pores in shales , 2017 .

[10]  George J. Hirasaki,et al.  Fluid Typing and Pore Size in Organic Shale Using 2D NMR in Saturated Kerogen Isolates , 2016 .

[11]  J. Birdwell,et al.  Updated methodology for nuclear magnetic resonance characterization of shales. , 2013, Journal of magnetic resonance.

[12]  Clemens Reimann,et al.  Factor analysis applied to regional geochemical data: problems and possibilities , 2002 .

[13]  Hugh Daigle,et al.  Altering nuclear magnetic resonance surface relaxation on nanoparticles by Adsorption of Fe(III) , 2018 .

[14]  Fionn Murtagh,et al.  A Survey of Recent Advances in Hierarchical Clustering Algorithms , 1983, Comput. J..

[15]  Joydeep Ghosh,et al.  Data Clustering Algorithms And Applications , 2013 .

[16]  Carl H. Sondergeld,et al.  NMR Response of Brine, Oil and Methane in Organic Rich Shales , 2014 .

[17]  Anil K. Jain,et al.  Data clustering: a review , 1999, CSUR.

[18]  Benjamin King Step-Wise Clustering Procedures , 1967 .

[19]  Mohamed Mehana,et al.  Shale characteristics impact on Nuclear Magnetic Resonance (NMR) fluid typing methods and correlations , 2016 .

[20]  Tian Zhang,et al.  BIRCH: an efficient data clustering method for very large databases , 1996, SIGMOD '96.

[21]  Gaël Varoquaux,et al.  Scikit-learn: Machine Learning in Python , 2011, J. Mach. Learn. Res..

[22]  Carl H. Sondergeld,et al.  NMR Study of Shale Wettability , 2011 .

[23]  J. MacQueen Some methods for classification and analysis of multivariate observations , 1967 .

[24]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[25]  Han Jiang,et al.  Characterization of Thermal Evolution of Pores and Fluids in Shales Using NMR 2D Measurement , 2018 .

[26]  Paolo Ferraris,et al.  Characterization of Underlying Pore and Fluid Structure Using Factor Analysis on NMR Data , 2013 .

[27]  J. H. Ward Hierarchical Grouping to Optimize an Objective Function , 1963 .

[28]  Carlos Torres-Verdín,et al.  Multidimensional NMR inversion without Kronecker products: Multilinear inversion. , 2016, Journal of magnetic resonance.

[29]  Ralf Littke,et al.  Geochemical and petrophysical source rock characterization of the Vaca Muerta Formation, Argentina: Implications for unconventional petroleum resource estimations , 2017 .

[30]  J.T.H. Ng,et al.  An Improved Temperature-Viscosity Correlation For Crude Oil Systems , 1983 .

[31]  J. Mullen Petrophysical Characterization of the Eagle Ford Shale in South Texas , 2010 .

[32]  Justin E. Birdwell,et al.  Multivariate Analysis Relating Oil Shale Geochemical Properties to NMR Relaxometry , 2015 .

[33]  Hugh Daigle,et al.  Direct Determination of Surface Relaxivity in Isolated Kerogen by Pulsed-Field Gradient NMR , 2016 .

[34]  Gerald Hamon,et al.  Estimating Saturations in Organic Shales Using 2D NMR , 2016 .

[35]  Anil K. Jain Data clustering: 50 years beyond K-means , 2008, Pattern Recognit. Lett..

[36]  Jitendra Malik,et al.  Normalized cuts and image segmentation , 1997, Proceedings of IEEE Computer Society Conference on Computer Vision and Pattern Recognition.

[37]  N. E. Day Estimating the components of a mixture of normal distributions , 1969 .

[38]  Lalitha Venkataramanan,et al.  An Unsupervised Learning Algorithm to Compute Fluid Volumes From NMR T 1 -T 2 Logs in Unconventional Reservoirs , 2018 .

[39]  E. Purcell,et al.  Effects of Diffusion on Free Precession in Nuclear Magnetic Resonance Experiments , 1954 .

[40]  Kathryn E Washburn,et al.  Detection of intermolecular homonuclear dipolar coupling in organic rich shale by transverse relaxation exchange. , 2017, Journal of magnetic resonance.

[41]  Han Jiang,et al.  Multiscale pore networks and their effect on deformation and transport property alteration associated with hydraulic fracturing , 2017 .

[42]  Lizhi Xiao,et al.  Advanced fluid-typing methods for NMR logging , 2011 .

[43]  Mukul M. Sharma,et al.  Porosity Evaluation of Shales Using NMR Secular Relaxation , 2014 .

[44]  Hugh Daigle,et al.  Porosity-Deformation Relationships in Organic-Rich Shale , 2019, Memoir 120: Mudstone Diagenesis: Research Perspectives for Shale Hydrocarbon Reservoirs, Seals, and Source Rocks.

[45]  S. Meiboom,et al.  Modified Spin‐Echo Method for Measuring Nuclear Relaxation Times , 1958 .