The significance of the study begins with the basin itself. The Hetao Basin, cradled along the upper reaches of the Yellow River in Inner Mongolia, is a Cenozoic downwarp filled with thousands of meters of Quaternary lacustrine and fluvial sediments. During the late Quaternary, vast paleolakes and river systems deposited thick successions of organic-rich mud that later became home to thriving methanogenic microbial communities. Those microbes converted buried organic matter into methane at shallow burial depths and low temperatures, generating biogenic gas accumulations in near-surface aquifer systems. Previous research by the same group revealed, through geochemical, genomic, and transcriptomic analyses, the microbial pathways behind this methanogenesis, and linked gas production to climate-driven increases in organic matter burial after the Mid-Pleistocene Transition. But understanding where the gas came from is only half the problem. To predict where it sits, how it migrates, and how it might be produced economically, geologists need to understand the pore structure of the weakly consolidated mudstone that acts as both source and reservoir—a material so friable that it can crumble between fingers, defying the standard laboratory techniques developed for well-cemented shales.

To tame this difficult material, the team turned to a suite of complementary analytical methods. Field-emission scanning electron microscopy allowed them to image pore types directly at nanometer to micron scales, revealing that interparticle pores—voids preserved between loosely packed mineral grains—dominate the pore system of these young sediments. This is a striking contrast with deeply buried, thermally mature shales, where organic matter pores and intraparticle dissolution pores often take center stage. The researchers then employed low-pressure nitrogen and carbon dioxide gas adsorption, the workhorse techniques for probing micropores below roughly 2 nanometers and meso-macropores from 2 to several hundred nanometers, interpreting the isotherms using established frameworks such as the Brunauer-Emmett-Teller surface area method and density functional theory models. The adsorption data yielded detailed pore size distributions spanning the full IUPAC classification, providing the quantitative foundation upon which the multifractal analysis was built.

Multifractal theory is where the study makes its most distinctive technical contribution. Whereas a single fractal dimension reduces a pore network to one number describing overall complexity, a multifractal spectrum captures how complex that network is at every scale and how heterogeneously pore volume is distributed across it. The researchers computed the generalized Rényi dimensions and the singularity spectrum f(α) from the adsorption-derived pore size distributions, following methods pioneered by Chhabra and Jensen and widely applied to soils, coals, and shales. Two parameters emerged as the critical diagnostics: the capacity dimension H, which reflects the connectivity and space-filling character of the pore system, and the width of the singularity spectrum Δα, which measures the degree of heterogeneity across the pore size range. A high H value indicates a well-connected, homogeneous network; a large Δα signals that the pore population is dominated by a few highly concentrated zones while other regions remain sparse.

The results revealed a scale-dependent division of labor within the reservoir. Micropores—the smallest voids, those below about 2 nanometers—exhibit remarkably high and nearly uniform capacity dimensions, with H values ranging from 0.970 to 0.974 and averaging 0.972. Values this close to unity indicate that the micropore network is extremely well connected and space-filling, behaving almost like a homogeneous, smoothly distributed medium. In practical terms, these tiny pores act as the adsorption engine of the reservoir, providing abundant surface area where methane molecules cling to mineral surfaces, and their excellent connectivity suggests that gas stored in them can access migration pathways efficiently. In contrast, the meso- to macropore range, spanning roughly 2 nanometers to tens of nanometers and beyond, displayed substantially greater heterogeneity, with singularity spectrum widths Δα between 0.761 and 0.832 and an average of 0.803. This pronounced heterogeneity means that larger pores are unevenly distributed, clustered in localized zones rather than spread uniformly through the rock. These larger pores, however, are precisely the conduits through which gas migrates and the storage vaults for free gas. The combination—an open, connected micropore network feeding into a heterogeneous network of larger channels—defines the dual storage-and-transport character of these shallow biogenic gas systems.

What controls this heterogeneity? The team’s mineralogical analysis pointed squarely at composition, with clay minerals and quartz emerging as the primary drivers of pore structure variability. Clay-rich intervals, with their platy, high-surface-area grains, generate abundant small pores and complex tortuous networks, while quartz-rich intervals preserve larger interparticle voids that enhance porosity but introduce patchiness. Because weakly consolidated Quaternary mudstone has experienced only shallow burial and minimal diagenetic alteration, mineralogy imposed during deposition remains the dominant imprint on pore architecture, unmasked by the cementation and compaction that homogenize older rocks. The multifractal parameters thus serve as sensitive mineralogical fingerprints, potentially allowing pore heterogeneity to be predicted from routine well logs or X-ray diffraction data in exploration settings where core recovery of unconsolidated material is difficult or impossible.

Perhaps the most intellectually ambitious component of the study is its comparative dimension. By placing the Hetao Basin mudstone alongside shales from a wide spectrum of geological settings—marine shales, continental lacustrine shales, transitional shales, coals, and organic-rich mudrocks spanning maturity ranges from immature to overmature—the researchers distilled the universal controls on pore heterogeneity into four factors: diagenesis, mineral composition, organic matter content, and thermal maturation evolution. In thermally mature shales, organic matter pores created by hydrocarbon generation dominate, and maturation progressively reshapes the pore system; in the Hetao Basin, by contrast, the rocks are so young and thermally inert that organic matter has never generated its own pore network. This natural experiment isolates the roles of the remaining factors with unusual clarity. It demonstrates that the multifractal signature of a reservoir is not a fixed property of mudrocks in general but a moving target shaped by burial history, temperature, and time—knowledge that allows explorationists to anticipate which pore systems to expect before drilling.

The practical implications for energy exploration are considerable. Shallow biogenic gas is regarded as a clean unconventional resource precisely because it is nearly pure methane, produced by microbes rather than thermal cracking, and it has already proven commercially viable in China’s Qaidam Basin, where the Sebei gas field demonstrates that Quaternary mudstone-hosted biogenic gas can sustain large-scale production. The Hetao Basin, with its thick organic-rich Quaternary fills and documented methanogenic activity in aquifer systems, represents one of China’s most promising frontier plays for this resource. Accurate reservoir evaluation in such settings demands knowledge of which intervals store gas, which transmit it, and where seals might be weak—questions that hinge directly on pore structure. The multifractal parameters established in this study offer quantitative targets that can be inverted from geophysical data, refined with machine learning approaches, and used to rank drilling prospects. Better prediction of migration pathways also carries environmental dividends, helping to avoid unintended gas leakage into shallow groundwater or the atmosphere during development.

The study also contributes to a growing appreciation that fractal geometry is the natural language of porous geomaterials. From soils and coals to carbonates and shales, researchers have repeatedly found that pore systems obey scale-invariant statistics, and multifractal analysis has proven especially powerful for bridging nano- and micron-scale observations obtained from gas adsorption, mercury intrusion, nuclear magnetic resonance, and computed tomography. The Hetao Basin work extends this toolkit to one of its most challenging applications: rocks that are barely consolidated, where sample preparation itself is delicate and where the very concept of a “pore” blurs between sedimentary void and soft-sediment fabric. That the method yields clean, interpretable spectra even in this difficult material suggests it can be deployed broadly across young sedimentary basins worldwide, from the Dutch North Sea to Hangzhou Bay, wherever microbial gas accumulations are being evaluated.

Looking forward, the researchers suggest that their findings will aid prediction of gas migration and storage in unconsolidated mudrock reservoirs and support exploration and development of shallow biogenic gas systems. The work was funded by the Fundamental and Interdisciplinary Disciplines Breakthrough Plan of the Ministry of Education of China, the Geological Exploration Fund of Inner Mongolia Autonomous Region, and the Natural Science Foundation of Hubei Province. As global energy systems transition toward lower-carbon sources, resources like shallow biogenic gas—abundant, shallow, cheap to produce, and composed almost entirely of methane—may occupy an increasingly important niche as both a transitional fuel and a subject of fundamental geoscience. With this study, the invisible pore networks of the Hetao Basin’s youngest rocks have finally been rendered in mathematical detail, turning a crumbly, overlooked mudstone into a quantitatively characterized reservoir and demonstrating that even the most mundane sediments can yield profound insights when examined through the lens of fractal geometry.