A 4D synchrotron X-ray tomography study of the formation of hydrocarbon migration pathways in heated organic-rich shale
Recovery of oil from oil shales and the natural primary migration of hydrocarbons are closely related processes that have received renewed interests in recent years because of the ever tightening supply of conventional hydrocarbons and the growing production of hydrocarbons from low permeability tight rocks. Quantitative models for conversion of kerogen into oil and gas and the timing of hydrocarbon generation have been well documented. However, lack of consensus about the kinetics of hydrocarbon formation in source rocks, expulsion timing and how the resulting hydrocarbons escape from or are retained in the source rocks motivates further investigation. In particular, many mechanisms for the transport of hydrocarbons from the source rocks in which they are generated into adjacent rocks with higher permeabilities and smaller capillary entry pressures have been proposed, and a better understanding of this complex process (primary migration) is needed. To characterize these processes it is imperative to use the latest technological advances. In this study, it is shown how insights into hydrocarbon migration in source rocks can be obtained by using sequential high resolution synchrotron X-ray tomography. Three-dimensional (3D) images of several immature "shale" samples were constructed at resolutions close to 5 micrometers. This is sufficient to resolve the source rock structure down to the grain level, but very fine grained silt particles, clay particles and colloids cannot be resolved. Samples used in this investigation came from the R-8 unit in the upper part of the Green River Shale, which is organic rich, varved, lacustrine marl formed in Eocene Lake Uinta, United States of America. One Green River Shale sample was heated in-situ up to 400°C as X-ray tomography images were recorded. The other samples were scanned before and after heating at 400°C. During the heating phase, the organic matter was decomposed, and gas was released. Gas expulsion from the low permeability shales was coupled with formation of microcracks. The main technical difficulty was numerical extraction of microcracks that have apertures in the 5 to 30 micrometer range (with 5 micrometers being the resolution limit) from a large 3D volume of X-ray attenuation data. The main goal of the work presented here is to develop a methodology to process these 3D data and image the cracks. This methodology is based on several levels of spatial filtering and automatic recognition of connected domains. Supportive petrographic and thermogravimetric data were an important complement to this study. An investigation of the strain field using two-dimensional image correlation analyses was also performed. As one application of the four-dimensional (4D, space + time) microtomography and the developed workflow, we show that fluid generation was accompanied by crack formation. Under different conditions, in the subsurface, this might provide paths for primary migration.