In a preprint currently under review for the data-focused journal Earth System Science Data, Heather Anne Sheldon of CSIRO Mineral Resources in Canberra and her colleagues across CSIRO’s energy, environment and information management divisions describe the creation of a Global Petrophysical Database, or GPD, that consolidates more than 112,000 individual rock property measurements drawn from over 600 publications. The compilation spans data from 70 countries and incorporates ten existing databases that had previously been maintained in isolation. Crucially, the database includes both in-situ measurements, taken where the rock actually sits in the subsurface, and laboratory measurements performed on samples brought to the surface, allowing users to compare how properties behave under natural conditions and under controlled experimental ones.

The technical scope of the resource is what sets it apart from earlier attempts. Previous petrophysical databases have tended to be limited in geographic extent, or restricted in the range of properties, rock types or measurement techniques they cover. The GPD deliberately cuts across those boundaries, encompassing igneous, metamorphic and sedimentary rocks alike, and focusing primarily on the properties that control subsurface fluid flow, mass transport and heat transport. Each measurement is accompanied by associated metadata, including rock type, geographic location and depth, the contextual information that transforms a bare number into something a modeller can actually use. Without that context, a permeability value is nearly meaningless, because the same rock type can behave radically differently at two kilometres depth than it does in a surface outcrop.

The practical motivation behind the project is easy to understand for anyone who has attempted a literature synthesis in this field. As the authors note, collating petrophysical property data from the literature is time consuming and is frequently hampered by inconsistent or incomplete metadata. A hydrogeologist modelling aquifer recharge might spend weeks tracking down porosity measurements for a particular sandstone formation, only to find that half the relevant papers omit the depth at which samples were collected or fail to specify the measurement technique. Those omissions matter enormously. Permeability measured on a plug in a laboratory at ambient conditions can differ by orders of magnitude from the effective permeability of the same formation in the ground, where confining stress, temperature and natural fracturing all play their part.

The applications the authors identify span some of the most consequential geoscience challenges of the coming decades. Groundwater management depends on accurate representations of how water moves through aquifers and the low-permeability layers that confine them. Carbon geo-sequestration, the injection of carbon dioxide into deep geological formations as a climate mitigation strategy, requires confident knowledge of reservoir porosity and caprock integrity. Resource exploration, whether for minerals, hydrocarbons or geothermal energy, relies on petrophysical properties to interpret geophysical surveys and to build predictive models of ore-forming fluid flow. And subsurface waste disposal, including the long-term storage of hazardous materials, demands an exceptionally rigorous understanding of how slowly, or quickly, fluids and dissolved species can migrate through host rocks over thousands of years.

What makes the database genuinely powerful for the modern era is its accessibility and its machine-readability. The GPD is publicly available through a persistent digital object identifier hosted by CSIRO, and it is supported by a graphical web interface called the Global Petrophysical Database Explorer, which is designed to facilitate rapid exploration and visualisation of the data. A researcher can query the collection, filter by rock type or property, and visualise distributions without writing a single line of code, while those building numerical models can ingest the underlying data directly. In an age when machine learning is increasingly applied to geoscience problems, a large, well-documented, consistently structured dataset of this kind is exactly the raw material such methods require, and its absence has been a persistent bottleneck.

The scale of the underlying effort should not be understated. Harmonising data from more than 600 publications means confronting a babel of units, measurement standards, rock classification schemes and reporting conventions that have evolved independently across subdisciplines and decades. A permeability reported in millidarcies by a petroleum engineer, one reported in square metres by a hydrogeologist, and a third reported only qualitatively in an old mining report must all be reconciled into a coherent framework. The decision to fold in ten pre-existing databases, rather than starting from scratch, reflects a pragmatic philosophy: the community has already generated an enormous volume of valuable measurements, and the highest-value contribution is to make them findable, comparable and reusable, in line with the open-data principles that journals such as Earth System Science Data were created to champion.

The timing is significant. As governments accelerate investment in carbon capture and storage, geothermal energy and managed aquifer recharge, demand for reliable subsurface property data is rising sharply, and the cost of poor data is measured not just in wasted research time but in flawed decisions about where to inject carbon dioxide or how to site critical infrastructure. A single, globally scoped reference dataset allows modellers to benchmark their assumptions against real measurements from analogous rock types and settings worldwide, rather than defaulting to textbook values that may bear little resemblance to the formations under study. It also exposes gaps: where the database is thin, it points directly to where new field and laboratory campaigns would deliver the greatest scientific return.

The work is currently at the preprint stage, published as a discussion paper with peer review open and ongoing, so the details may evolve as referees and community members weigh in. But the core contribution is already tangible and freely available. For the researchers who built it, the goal was straightforward: to save others the time they themselves had spent hunting for rock property values, and to give the geoscience community a shared foundation for modelling and decision making. If the database achieves the uptake its designers hope, the era of rebuilding the same scattered literature review for every new subsurface project may finally be drawing to a close, replaced by a single, searchable, global record of what we know about the rocks beneath our feet.

Subject of Research: A global petrophysical database compiling rock property measurements for subsurface modelling

Article Title: A global petrophysical database for igneous, metamorphic and sedimentary rocks

Article References: Sheldon, H. A., Dewhurst, D. N., Bekele, E., Raiber, M., Mallants, D., Turnadge, C., & Williams, G. (2026). A global petrophysical database for igneous, metamorphic and sedimentary rocks. https://doi.org/10.5194/essd-2026-556

Image Credits: AI Generated

DOI: 10.5194/essd-2026-556

Keywords: petrophysics, rock properties, porosity, permeability, thermal conductivity, database, groundwater, carbon sequestration, geothermal energy, subsurface modelling, CSIRO, open data

The research, conducted by a team from GEOPIC, the research arm of Oil and Natural Gas Corporation in Dehradun, along with collaborators at Pandit Deendayal Energy University and Maharashtra Institute of Technology World Peace University, focused on the gas-bearing tight sandstones of the Jhuran Formation. These rocks, deposited during the syn-rift and early post-rift phases that followed the fragmentation of Gondwana, are anything but uniform. Permeability values range from a near-impermeable 0.01 millidarcies to a still-modest 10 millidarcies, and the formation is sliced into parallel fault blocks by northwest-southeast trending structures that complicate any attempt to extrapolate reservoir quality between wells. For an emerging hydrocarbon province like the Kachchh Offshore Basin, that heterogeneity is a formidable barrier to exploration.

The conventional toolkit offers little comfort. Classic permeability estimation leans on empirical porosity-permeability transforms calibrated to laboratory core measurements. In homogeneous reservoirs such relationships perform adequately, but in tight, clay-rich sandstones the link between porosity and permeability weakens dramatically, because flow is governed by pore geometry, grain packing, clay distribution and cementation rather than by pore volume alone. Compounding the problem, core plugs are expensive and sparse; in this study, direct core measurements existed for just one well, and only fifteen of them. The team therefore turned to machine learning, training four supervised algorithms on wireline log data from four wells: Random Forest, Support Vector Regression, a Multi-Layer Perceptron neural network, and Extreme Gradient Boosting.

The input features were the workhorses of petrophysics: gamma ray, deep resistivity, bulk density, neutron porosity and sonic transit time, supplemented by engineered variables such as volume of clay, total and effective porosity, and a density-neutron-resistivity combination parameter. Nuclear magnetic resonance logging provided an independent permeability benchmark in two wells, derived through the Schlumberger-Doll Research equation, which relates permeability to porosity and the logarithmic mean of the NMR T2 relaxation distribution. Crucially, the NMR-derived permeability served only as the target variable, never as an input, a design choice that prevents data leakage and ensures the models can be applied to wells lacking NMR coverage. Spearman rank correlation analysis confirmed that no input pair approached the severe multicollinearity threshold, supporting their joint inclusion.

Validation was deliberately stringent. Wells A, C and D supplied the training data, split 70-30 into training and testing subsets, with five-fold cross-validation used for hyperparameter tuning. Well B, which contained the fifteen core measurements, was held out entirely and used only once, after model finalisation, as a blind validation well. On that single held-out well, Random Forest achieved a coefficient of determination of 0.77 against the core and NMR benchmarks, the best cross-well generalisation of the four algorithms, and posted the lowest prediction variance, with a residual standard deviation of 1.41. Support Vector Regression, by contrast, fit its training data beautifully, reaching a training R-squared of 0.92, but its validation performance collapsed, the classic signature of overfitting.

The contrast between algorithms carried a practical lesson. The Multi-Layer Perceptron and XGBoost both achieved near-perfect alignment on training data yet generalised poorly, with XGBoost showing the largest gap between training and blind-validation performance and the highest prediction uncertainty in Well A, where residual standard deviations reached 7.81. Random Forest, averaging hundreds of decision trees each trained on random subsets of the data, proved the most stable performer across all wells, with residual standard deviations as low as 0.38 in Well C. The authors suggest a heuristic for similarly data-starved settings: Random Forest as the default choice, Support Vector Regression where smooth low-variance surfaces are needed, and neural networks or gradient boosting only with strong regularisation and blind-well checks, or once more wells become available.

What elevates the study beyond a routine algorithm bake-off is its use of SHAP, or Shapley Additive Explanations, a technique borrowed from game theory that decomposes each prediction into the contributions of individual input features. The SHAP analysis revealed a clear hierarchy of geological control. Effective porosity exerted the strongest first-order influence on permeability, followed by sonic transit time, which acts as a proxy for compaction. Most strikingly, shale volume displayed strongly negative SHAP values beyond a threshold of roughly 35 to 40 percent, marking a permeability collapse in clay-rich intervals. In other words, permeability in the Jhuran sandstones is governed by nonlinear interactions among shale content, porosity and compaction, not by porosity alone, precisely the kind of behaviour that linear empirical transforms cannot capture.

The team was candid about the limits of its evidence base. Four wells with uneven log and core availability is a small dataset, and the NMR calibration constants, drawn from standard literature values for water-wet clastic systems, were checked only qualitatively against the fifteen core samples from Well B rather than formally recalibrated for the Jhuran Formation. The reported R-squared of 0.77 reflects a single blind well, and the authors explicitly caution that basin-wide generalisation requires confirmation with additional held-out wells and more core control. Well D, lacking both core and NMR data, relied on a facies-derived pseudo-benchmark, an independent consistency check but not a direct measurement. These caveats frame the work honestly as a proof of concept rather than a finished predictive product.

Even so, the implications reach well beyond one basin. Frontier exploration provinces worldwide face the same predicament: geologically complex reservoirs, a handful of wells, and sparse core data, all while drilling decisions carry enormous financial risk. The study demonstrates that a transparent, reproducible machine learning workflow, built on open-source libraries and validated against independent physical measurements, can extract usable permeability predictions from ordinary wireline logs under exactly those constraints. By pairing ensemble learning with explainability, the researchers have shown that machine learning need not be a black box; it can double as an analytical instrument that tells geoscientists which geological processes matter most. As exploration accelerates along India’s northwestern continental margin, that combination of prediction and insight may prove as valuable as the predictions themselves, offering a transferable template for characterising tight reservoirs wherever data is scarce and heterogeneity is the rule.

Subject of Research: Interpretable machine learning prediction of permeability in tight Jhuran Formation sandstones of the Kachchh Offshore Basin, India

Article Title: Data driven permeability modelling in tight Jhuran formation sandstones of the Kutch offshore Basin India

Article References: Jyotsna, G. D. D., Anees, M., Mavani, S. H., Thakur, S., Saxena, A., Singh, S. K., Desai, D. B. S. G., & Singh, D. A. (2026). Data driven permeability modelling in tight Jhuran formation sandstones of the Kutch offshore Basin India. Discover Geoscience, 4(1), Article 366. https://doi.org/10.1007/s44288-026-00724-x

Image Credits: AI Generated

DOI: 10.1007/s44288-026-00724-x

Keywords: permeability prediction, machine learning, tight sandstone, Jhuran Formation, Kachchh Offshore Basin, Random Forest, SHAP analysis, well logging, NMR logging, reservoir characterization, petrophysics, explainable AI

The central problem the researchers set out to address is a mismatch of scales and definitions. Laboratory porosimetry, and especially helium porosimetry on core plugs, provides accurate reference values for pore parameters, but these are discrete point measurements that do not adequately capture the spatial variability of pore fabrics across a reservoir. Digital image analysis, by contrast, offers a non-destructive way to quantify the visible pore space in thin-section images of rock. The catch is that image analysis yields a two-dimensional geometric porosity value that says nothing about whether the pores it counts are actually connected, accessible to fluids, or partially blocked by clay minerals that coat and choke pore throats. A rock can look porous in a photograph yet perform poorly as a reservoir.

The new workflow, which the authors call DIA–FIS, was tested on sandstone samples from the Petro-Masilah reservoir. The first stage relies on thin sections impregnated with blue epoxy, a standard petrographic preparation in which the resin fills the pore space and renders it a distinctive color under the microscope. The researchers applied pre-fixed thresholds in the hue–saturation–value color space, a representation of digital images that separates chromatic content from brightness and is often more robust to uneven illumination than working in red–green–blue channels alone. Segmentation in this space allows pore space to be isolated from grains and cements, and the resulting image-derived porosity can be calibrated against measured values.

A crucial step in the study was an explicit test of how sensitive this segmentation is to the choice of threshold boundaries. Because any threshold-based segmentation involves a judgment about where pore pixels end and grain pixels begin, the researchers conducted stability tests by perturbing their hue–saturation–value thresholds by plus or minus five percent and plus or minus ten percent. This kind of perturbation analysis quantifies the uncertainty that is usually left implicit when an analyst simply picks a threshold by eye. It also addresses a known weakness in the broader literature: image processing variability, including the choice between conventional thresholding and machine learning approaches, can meaningfully change pore-structural results.

The second stage of the workflow takes the calibrated image-derived porosity together with the measured clay volume ratio and feeds both into a Mamdani-type fuzzy inference system. Mamdani systems, one of the classic architectures of fuzzy logic, operate on linguistic variables and human-readable IF–THEN rules, in contrast to black-box machine learning models whose internal reasoning cannot easily be inspected. In this framework, statements such as rules linking higher image porosity and lower clay content to better reservoir quality are encoded explicitly, and the degree to which each rule fires is computed through membership functions that allow partial, gradual membership in categories such as low, medium, and high.

The output of the fuzzy stage is a defuzzified score the authors denote S_FIS, a dimensionless reservoir-quality score derived from the combined interpretation of image-derived porosity and clay volume. The authors are careful to define what this score is not: it is a rule-based assessment of reservoir quality rather than a physical porosity measurement, and it is not equivalent to effective, connected, non-clay, or helium-measured porosity. This distinction matters, because conflating a fuzzy quality index with a measured porosity would be a serious methodological error. The fuzzy output instead captures the gradual transitions between pore types and reservoir-quality classes, the very vagueness that hard classification schemes handle poorly.

The performance of the digital image analysis stage was evaluated against helium porosimetry performed on the same core plugs, one of the most reliable reference measurements available in petrophysics. The comparative analysis demonstrated good correlation for the image analysis stage, with a coefficient of determination greater than 0.95 and a mean absolute deviation below 3.5 percent. Those figures indicate that, once calibrated, simple threshold-based image segmentation tracks laboratory measurements closely for the visible pore fabric captured in two dimensions, even before any fuzzy reasoning is applied.

The fuzzy stage then did what it was designed to do: it penalized clay. Samples classified as clay-affected received lower reservoir-quality scores and lower quality classes from the inference system, while cleaner sandstones received more favorable interpretations. This behavior reflects genuine petrophysical reasoning. Clays, whether dispersed between grains or lining pore walls, reduce effective pore throats, impair permeability, and make a rock’s total porosity a misleading indicator of its fluid-storage and fluid-delivery capacity. By explicitly including clay volume as an input to the rule base, the model forces the interpretation to account for a factor that a purely geometric porosity number would ignore.

Interpretability is the distinguishing selling point of the approach. Deep learning methods have made remarkable inroads into digital rock physics, from segmenting digital rock images to predicting porosity from well logs, but their opacity can be a barrier in an industry where decisions worth millions ride on petrophysical interpretations and where regulators and partners may demand to know why a model produced a given answer. A Mamdani fuzzy system, by contrast, can be audited rule by rule: a geologist can read the IF–THEN statements, check whether they accord with established reservoir engineering understanding, and adjust membership functions if local rock properties demand it. The workflow thus positions itself as a transparent complement to, rather than a replacement for, conventional laboratory characterization.

The authors describe the DIA–FIS workflow as a reproducible, non-destructive complement to traditional reservoir characterization, and the framing is apt. The image analysis can be repeated on archived thin sections without consuming additional core material, the threshold perturbation tests make the segmentation uncertainty explicit rather than hidden, and the fuzzy rules encode petrophysical reasoning in a form that domain experts can verify and refine. For heterogeneous, clay-influenced sandstone reservoirs such as those of the Sayun-Masilah region from which the study samples derive, the approach offers a way to extract more interpretable information from routine petrographic preparations. The dataset underpinning the study comprises petrographic and chemical characterization of sandstone samples from the Petro-Masilah reservoir; because the core material is proprietary, the data are curated by the corresponding author and available upon reasonable request for academic and non-commercial research purposes. As machine learning continues to spread through the geosciences, studies like this one make the case that sometimes the most valuable model is not the most powerful one, but the one whose reasoning a reservoir engineer can actually read.

Subject of Research: Hybrid digital image analysis and fuzzy logic for porosity and reservoir-quality estimation in clay-influenced sandstone reservoirs

Article Title: Interpretable hybrid digital image analysis model with fuzzy logic system for porosity estimation in clay-influenced sandstone reservoirs

Article References: Interpretable hybrid digital image analysis model with fuzzy logic system for porosity estimation in clay-influenced sandstone reservoirs. (n.d.). https://doi.org/10.1007/s12145-026-02230-7

Image Credits: AI Generated

DOI: 10.1007/s12145-026-02230-7

Keywords: porosity estimation, digital image analysis, fuzzy logic, fuzzy inference system, sandstone reservoirs, thin sections, HSV segmentation, clay content, reservoir quality, helium porosimetry, petrophysics, interpretability

Researchers led by Bomi Kim of the Korea Institute of Geoscience and Mineral Resources, together with colleagues at Kongju National University, turned their attention to the Myeonsan titanium deposit in the Taebaeksan Basin of eastern South Korea. The deposit, hosted by the Cambrian Myeonsan Formation, is interpreted as a paleo-placer: titanium-bearing minerals such as rutile, ilmenite, and hematite were eroded from Precambrian basement rocks, transported by ancient rivers, and concentrated in conglomeratic sandstone during sedimentary deposition. Because the ore horizons are surrounded by a diverse suite of lithologies, including limestone, clastic sedimentary rocks, intrusive dykes, and granite, the site provides an ideal natural laboratory for testing which physical measurements can reliably separate ore-bearing rock from barren host rock.

The team selected 39 representative drill-core samples from borehole KDMS 2023-10, which penetrates the mineralized interval and its surroundings. The samples were divided into five groups: nine Myeonsan ore specimens, twelve limestones, four dyke rocks, eleven clastic sedimentary rocks, and three granites. Each specimen was prepared as a polished cylinder, a geometry that promotes uniform current distribution and simplifies the calculations needed for electrical measurements. Where core material was scarce, half-core semi-cylindrical specimens were prepared for the titanium ores, ensuring that every sample could pass through the full battery of tests.

That battery was deliberately broad. The researchers measured density and porosity using the Archimedes buoyancy method, following International Society for Rock Mechanics procedures: specimens were dried at 110 degrees Celsius for 24 hours, vacuum-saturated with water, and weighed in both saturated and submerged states. Magnetic susceptibility was recorded with a Bartington MS-3 meter at six positions on each sample to average out mineralogical anisotropy. Direct-current electrical resistivity, time-domain induced polarization (TDIP), and electrochemical impedance spectroscopy (EIS), known in geophysics as spectral induced polarization or SIP, were measured on fully saturated specimens using Solartron instruments with copper-mesh and filter-paper electrodes. Portable X-ray fluorescence (p-XRF) provided elemental compositions, and petrographic thin sections revealed the mineralogical texture behind the numbers.

The results deliver a clear hierarchy of diagnostic power. Density proved the most straightforward discriminator: the titanium ores averaged 3.17 grams per cubic centimeter, well above limestone at 2.74, dyke rocks at 2.73, clastic sediments at 2.84, and granite at 2.68. The elevated density reflects the abundance of heavy oxide minerals such as rutile and ilmenite packed into the ore horizons. The larger standard deviation among ore samples also revealed heterogeneous mineral enrichment within the mineralized interval, a useful warning that ore grades can vary significantly over short distances.

Chargeability, the capacity of a rock to store and slowly release electrical charge, produced an even more dramatic contrast. The ore samples averaged 91.57 millivolts per volt, roughly sixteen times higher than limestone at 5.64, dyke rocks at 3.41, clastic sediments at 4.62, and granite at 6.56. This signal arises because semiconductive ore minerals accumulate charge at mineral-fluid interfaces when an electric current is applied, then release it gradually after the current stops. In practical terms, an induced-polarization survey over this deposit should light up the ore zones while leaving most host rocks nearly invisible, making TDIP one of the most valuable exploration tools for sedimentary titanium deposits.

Not every measurement proved useful, and the negative results are as informative as the positive ones. Porosity values, ranging from 0.36 to 1.19 percent across groups, showed no meaningful separation between ore and host rocks, reflecting instead the compaction and cementation history of the sedimentary sequence. Magnetic susceptibility failed as a discriminator because the ore assemblage itself is magnetically quiet: rutile is only weakly magnetic and hematite is antiferromagnetic, so the intrusive dyke rocks, averaging 10.00 millisiemens per meter, outmagnetized the ores at 2.82. Electrical resistivity was similarly ambiguous, controlled mainly by pore structure and water saturation rather than ore content, with dyke and granite samples showing lower resistivities than the ores themselves.

The frequency-domain measurements added a powerful new dimension. In EIS spectra, several ore specimens, including MO1, MO2, MO5, and MO9, displayed strong phase responses at low frequencies below about 100 hertz, a signature of interfacial polarization between semiconductive ore minerals and pore fluids. These same samples showed the highest TDIP chargeability and elevated titanium contents above 3.11 weight percent and iron above 8.82 weight percent. Other ore samples and most host rocks instead showed phase values that increased toward higher frequencies, a pattern associated with dielectric polarization in insulating silicate and carbonate minerals. The two spectral fingerprints give geophysicists a way to distinguish genuine ore polarization from background dielectric effects, sharpening the interpretation of field survey data.

Petrography tied the physics back to mineralogy. Under the microscope, the ore specimens were dominated by opaque rutile, ilmenite, and hematite with minor quartz, while limestones were composed of calcite, dyke rocks of quartz, feldspar, and amphibole, clastic sediments of alternating coarse and fine silicate grains, and granites of quartz, feldspar, mica, and tourmaline. Sericitization of feldspar in the dykes and granites pointed to hydrothermal alteration. These observations directly explain the measured contrasts: dense, polarizable oxide minerals produce high density, high chargeability, and strong low-frequency phase responses, whereas insulating silicate and carbonate host minerals generate only weak, high-frequency dielectric behavior.

The study also exposed practical limits. One ore specimen, MO3, contained the highest titanium (9.85 weight percent) and iron (24.27 weight percent) of all samples yet showed weaker polarization than its peers, likely because p-XRF was measured at a single spot on a heterogeneous specimen. Because the core samples had to be preserved intact for the mining project, destructive analyses such as X-ray diffraction and ICP-MS were not possible. The authors conclude that spot geochemistry is best used to flag the presence of metallic minerals, while volume-integrating measurements such as density, TDIP, and EIS provide the more robust bulk characterization. Together, the integrated dataset establishes quantitative criteria for distinguishing titanium-bearing rocks from their hosts and demonstrates how combined petrophysical, geochemical, and mineralogical measurements can sharpen the search for a critical mineral whose secure supply increasingly depends on finding it faster and smarter.

Subject of Research: Integrated physical property measurements of titanium-bearing ores and host rocks at the Myeonsan deposit for mineral exploration

Article Title: Integrated physical property measurements of titanium-bearing and host rocks for mineral exploration

Article References: Kim, B., Shin, S., Shin, D., & Lee, G. (2026). Integrated physical property measurements of titanium-bearing and host rocks for mineral exploration. Environmental Earth Sciences, 85(16), Article 405. https://doi.org/10.1007/s12665-026-13140-8

Image Credits: AI Generated

DOI: 10.1007/s12665-026-13140-8

Keywords: titanium, mineral exploration, Myeonsan deposit, induced polarization, spectral induced polarization, density, chargeability, paleo-placer deposit, petrophysics, critical minerals, South Korea, rutile