  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{\deg}C as X-ray tomography images were recorded. The
other samples were scanned before and after heating at 400{\deg}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.
