Automatic Data Reduction and Quantification of X-Ray Computed Tomography Images of Sedimentary Cores: Method and Illustration

This paper presents a procedure from which information contained in 3-Dimensional single energy X-ray computed tomography (XR-CT) images of sedimentary rocks is converted into sub-mm scale resolution core scalar and core image logs. This new data provide a quantitative and compact (data volume reduction of ~90%) description of the XR-CT images. Density-related outputs are calibrated through automatic integration with continuous digital visual core description (VCD) and discrete moisture and density (MAD) property index measurements of selected samples. After lithology-based calibration of the X-ray attenuation coefficients into density values, quantitative displays include: 1) histogram of the distribution of density values and its related statistical parameters, 2) radial and angular distributions of core density values, 3) volume, average density and mass contributions of three core fractions defined by density thresholds corresponding to voids or vugs (VV, density ≤ ~1 g∙cm), and a break in the histogram of distribution of the density values showing the limit between the damaged (DM) and non-damaged (ND) fractions of the core material, and so, 4) providing a sub-mm scale bulk density core log free of any drilling disturbance. The procedure is illustrated on data from the 365 m deep Hole C9001C drilled off-shore Shimokita (northeast coast of Honshu, Japan). Usage of the outputs include: 1) derivation of sub-mm scale porosity core log, 2) correction of volume sensitive measurements in case of poor core quality and partially filled core liner, and 3) seismic modeling and well ties.


Introduction
Computer Tomography (CT) is a nondestructive, radiation recording technique, which maps the numerical values of the linear attenuation coefficient of X-or gamma-(photon) ray through a section of the investigated object. It requires the acquisition of two-dimensional radiographs for different positions during step-wise rotation around a central axis, whereby either the source and detector or the sample is moved. This is followed by the reconstruction of two-dimensional cross-sections perpendicular to the axis of rotation [1] [2]. First developed for applications in medicine [3], CT technology has been applied to various domains of geosciences [4] [5]. For more recent references, [6] reviews recent advances in computed tomography of geomaterials.
For work conducted on sedimentary rocks, CT investigation has both a qualitative and a quantitative approach. Qualitative investigations reveal various details of the core structure as bioturbation [7], heterogeneous discontinuities [8], alternate beddings in sandstones or partially filled fractures, vein selvage and vuggy porosity in limestone [9], biogenetic structures [10], invasion of the mud induced by the capillarity of the small scale porosity [11], facies succession composed of fine and thick rhythmites [12], and rotational fractures in air dried oil shales [13]. By assembling together more tomographic cross-sections, it is possible to obtain genuine 3-Dimensional (3D) representations of the investigated structures [14]. Because sedimentary structures reflect the environmental conditions during their accumulation, CT has been also used in description of facies as a result of transport and climatic processes [12] [15]. The quantitative approach in single energy CT scanner is related to the determination of the sediment bulk density. Recent quantitative use of dual energy CT scanner is not discussed in this paper. Based on core-to-core calibration, single energy CT has been used with good results to determine both vertical and transverse profiles of the bulk density of marine cores [12] [16]. The same technique has shown, for the oil shales of Cassuarina Basin (Australia), the presence of a strong positive correlation between attenuation coefficients and bulk density of air dried shales together with a moderately strong inverse correlation between the attenuation coefficients and modified Fisher assay oil yield [13]. The same kind of linear calibration curves have been experimentally confirmed for a number of synthesized cores composed by sand, silt and clay, saturated either with marine or fresh water [17]. Reference [18] quantitatively evaluated the attenuation contrast of various minerals using a medical X-ray CT scanner. More recently, [19] investigated the relationship between longitudinal (topograms) and transversal (tomograms) intensity scales to extract bulk density from its. Results showed excellent correlations between converted topogram values and Hounsfield Unit (HU) values obtained on the tomograms, as well as between relative bulk densities de-rived from the converted topogram values, and the absolute bulk density of various minerals.
In the framework of the Integrated Ocean Drilling Program (http://www.iodp.org/), non-invasive XR-CT scanning of cores is now routinely available on-board the Japanese drilling vessel, Chikyu, operated by the Japan Agency for Marine-Earth Science and Technology (JAMSTEC). As CT generates large data volumes, it requires considerable computer resources for visualization, analysis and archiving. As an example, medical-type X-ray Computed Tomography (XR-CT) imaging of a 1m-long core provides 1600 core transverse slices (i.e. one 512 × 512 pixel image every 0.625 mm) for a total data volume of about 200 Megabytes. A standard two months IODP expedition usually recovers at least 2 km of core resulting in 0.4 Terabytes of data and more than 3 million slices. The objective of this paper is to present a procedure for automatic data reduction (by about 90%) and automatic quantification of CT images facilitating integration of such 3D data set with standard on-board data and post-cruise results. The proposed method uses medical type CT images of cores to automatically 1) identify voids/vugs; 2) identify the damaged core fractions of the core and convert it into a semi-quantitative core quality indicator; 3) identify local variations of core linear attenuation coefficient in terms of radial and angular distributions; and once properly calibrated with lithology-based discrete moisture and density (MAD) property index measurements 4) derive a sub mm-resolution bulk density core log to further compare with routine multi-sensor core logger (MSCL) continuous gamma ray attenuation (GRA) density core log. The procedure is described step-by-step and illustrated with diatomaceous silty clays core data collected in 2007 at Hole C9001C (offshore Shimokita, northeast coast of Honshu, Japan) [20]. Result of NantroSEIZE Stage 1 (IODP Expeditions 314, 315 and 316) are available on-line at http://www.pgaillot.com/Work. Available core materials and on-board main core analysis techniques are presented. The principles, assumptions, acquisition and systematic errors of the MAD methodology are first introduced as of prime importance during calibration of the XR-CT derived density. Next, the principles and acquisition settings of GRA and XR-CT data collected on-board are presented. Detailed steps of the procedure for reduction of stored data volume and quantification of XR-CT images are described. MAD-density, GRA-density and XR-CT derived density are compared and discussed in term of data quality. Finally, usages of the outputs of the data reduction of the 3D XR-CT images for 1) derivation of sub-mm scale porosity core log, 2) correction of volume sensitive measurements in case of poor core quality and partially filled core liner, and 3) seismic modeling (synthetic seismogram) and well ties, are presented.

Materials
The procedure for automatic data reduction of core CT images and production P. Gaillot [20].
Following a standard on-board core processing workflow, the recovered cores (core barrel ~ 9.5 m long) were cut into ~1.5 m long sections and scanned with a medical-type XR-CT scanner (GE Medical System Light Speed Ultra16; [21]).
The X-ray tube voltage and current were set to 120 kV and 100 mA, respectively.
Measurement interval of the helical scan has been set to 0.625 mm under the 16 arrays of detectors mode. In this scan mode, three minutes are required to scan a 1.5 m section, five minutes for mathematical reconstruction using Image Works software (http://www.imageworkscorporation.com/productdrivers), and five additional minutes for data format conversion and transfer.
After XR-CT scanning, selected whole core segment (typically 10 -30 cm long) were taken for further investigation of ephemeral properties (interstitial water and microbiology analyses). The remaining whole core sections were scanned using multi-sensor core logger at 4-cm measurement intervals for gamma ray attenuation (GRA) density, magnetic susceptibility (MS), natural gamma radiation (NGR), resistivity (RES) and P-wave (PWV) velocity. The core sections were subsequently split, visually described in term of lithology and sedimentary features (visual core description, VCD), and finally sampled for other physical and chemical analyses [20]. For the purpose of this study, and especially validation of the XR-CT derived bulk density (ρ XRCT ), the most relevant among all these discrete analyses are MAD property index measurements resulting in high-quality (i.e. drilling disturbance free) core bulk density ( ρ MAD Bulk ), grain density ( ρ MAD SOLID ) and porosity ( φ MAD ) data. Before presenting the results of data reduction of XR-CT images and their applications, the principles, acquisition, and processing of the MAD, GRA and XR-CT scan data to the two-phase (minerals and interstitial water) system model for sediments and rocks are first summarized

Moisture and Density of Discrete Sample
In a technique known as MAD, moisture (porosity) and density of discrete samples can be estimated from mass and volume measurements of the sample before and after removal of interstitial pore fluid through drying via a convection oven.
The principle of this technique relies on phase-relation equations (extracted water by drying is assumed to be interstitial pore fluid) and correction for the mass and volume of salt precipitated during the drying process [22].
Practically, MAD measurements of the current data set were routinely meas- Open Journal of Geology ured using ~10 cm 3 samples from split cores collected at a frequency of two per core (i.e., ~ every 5 m), taking care to sample undisturbed parts of the core and avoid drilling slurry and "biscuits" (broken pieces of cores). Sampling frequency was increased as needed to characterize all significant changes of lithology. Immediately after samples were collected, wet sediment mass was measured. Samples were then placed in a convection oven for 24 hours at a temperature of 105˚C ± 5˚C. After drying, dry sediment mass and dry sediment volume were measured. Sample mass was determined on board to a precision of ± 0.01 g using an electric motion compensated balance system (error less than 1% with 0.1 g of sample). Volumes were determined using a helium five-chambered pycnometer with an approximate precision of ± 0.02 cm 3 .

Gamma Ray Attenuation Density
Bulk density of sediments and rocks can also be estimated from the measurement of gamma ray attenuation (GRA) of medium-energy gamma rays (100 -1000 keV) [23]. The continuous measurement of whole core section density relies on a 137 Cs source (used transmits gamma rays at 660 keV) and a standard NaI scintillation detector coupled with a photo-multiplier. In the on-board equipment setting, the collimator used is 5 mm in diameter and the spatial resolution of the sensor is about 0.5 cm. Measurements with the present system have typically count rates of 20,000 cps (soft mud). If measured for 4s, as in the current core flow, the statistical error is therefore less than 40,000 ± 200, or 0.5%.
However, accuracy is limited by the assumption that the measured material has the same attenuation coefficient as the calibration standards used. For general sediment water mixtures, this should be the case and errors should be less than 5%.

X-Ray Computed Tomography Images
The operating energy range of the present-day XR-CT scanners is between 30 and 200 keV, and the linear attenuation coefficient μ is a function of both bulk density and effective atomic number of the sample.
Image quality optimization requires a reduction of artifacts, which are inherent to the XR-CT technique. Several common types of artifacts can already be reduced to some extent by optimizing acquisition conditions, but a further reduction during reconstruction is often needed [24]. The most important artifacts in XR-CT images are related to "beam hardening", which refers to the preferential attenuation of low-energy X-rays when a polychromatic X-ray beam, typical of medical XR-CT scanners equipped with X-ray tubes, passes through an object.
For beam-hardening artifact reduction during reconstruction, linearization procedures are highly recommended for monomineralic systems [24]. For polymineralic systems, a dual energy approach allows a dramatic improvement in image quality [25].
For inhomogeneous samples, the attenuation coefficient has different values where w µ and a µ are the linear attenuation coefficients of the water and air, respectively.

Data Reduction of XR-CT Images
The result of the mathematical reconstruction and reformatting of the 3D XR-CT images is a set of transversal core slices (tomograms) with a separating distance of 0.625 mm between slices and a pixel size of 0.1875 mm × 0.1875 mm within the slice. A total of ~ 580,000 slices were recorded for this core dataset, and, thus the need for development of a procedure for automatic data reduction.
As the 3D images of cores is a set of co-axial (aligned) core transversal 2D slices, the procedure for automatic data reduction of XR-CT images of cores is auto-

Core Extraction (Step 1)
XR-CT slices are 512 pixels × 512 pixels in size DICOM (Digital Imaging and Communication in Medicine) format images in HU. The first step of the procedure is to extract the core image by detecting the plastic core liner and so, isolate the core material representing the only portion of the image of geological interest (i.e. the core itself). At the resolution of the image, the position, diameter, thickness and shape of the core liner can vary from slice to slice. However, the core liner can be detected by the contrast in attenuation coefficient of air (negative HU) and plastic (positive HU; see Table 1). This a-priori simple step is sometimes complicated by localized mathematical artifacts and/or partial volume effects inherent to CT reconstruction. Example of artifacts and impact on core liner detection are illustrated in Figure

Calibration of X-Ray Attenuation Coefficient into Bulk Density (Step 2)
An empirical relationship between XR-CT attenuation coefficient and bulk density, ρ, can be established by scanning materials of known bulk density. To avoid  where the background is automatically set to black, the core-liner is set to gray and the internal core liner fraction is set to white. Note that in this set of images, the core liners are not initially displayed as a continuous thick gray zone as it could be expected in correctly reconstructed images where the core liner detection is successful (b).
any complications related to XR-CT scanner geometry, mathematical reconstruction algorithm and physical photoelectric and Compton effects, a normal practice is to use pellets of controlled bulk density made of the same rock material as the core analyzed [26]. In the current workflow, discrete MAD samples combined with digital VCD have been used as lithology-controlled natural pellets. The range of calibration has been expended by using air, water, glass and aluminum standards. The resulting empirical relationship (Equation (2)) relating ρ MAD BULK to XR-CT HU values is slightly non-linear (second order polynom, R 2 ~ 0.75; Figure 2).   Table 1). The size of the symbols is proportional to the standard deviation of the HU and density values of the calibration pieces.

Distribution of Bulk Density Values (Step 3)
Once the XR-CT attenuation coefficients of each core slice are converted into density values (Figure 3 Table 2.

Thresholding and Phase Identification (Step 4)
From the histogram of distribution of density values, one can define two thresholds (Figure 3(b)). The first threshold (Tvv) corresponding to the standard pore water density ( ρ STANDARD   Two XR-CT-derived bulk density measurements are of primary importance while comparing those density estimates with available MAD-and GRA-density data for this core dataset: 1) the first measurement is the average bulk density of the slice, ρ XRCT TOTAL , that equally integrates the VV, DM and ND fractions of the core material, 2) the second measurement is the average density of the isolated ND fraction of the core, bulk density of the slice, ρ XRCT ND , thought to be free of any drilling disturbance as voids/vugs and core damaged effects are not included in this measurement. In the example of Figure 3, where the VV and DM fractions respectively equal 10.54% and 11.21% of the slice volume, ρ XRCT TOTAL equals 1.31 g•cm −3 while ρ XRCT ND = 1.42 g•cm −3 , a 7.8 % difference in bulk density value.
Percentage of VV and DM core fraction can further be used to automatically provide a semi-quantitative core quality estimator color-coded according to user define rules. The rules used in this work are given in Table 3 and resulting distribution of core quality estimates is given in Figure 5(f) and Figure 6.

Angular and Radial Distribution (Step 5)
After having computed the center of the core, the angular and radial distributions of the original XR-CT attenuation coefficients and/or bulk density values after calibration (see Step 2) are computed. The angular XR-CT attenuation coefficients (resp. density values) distribution is computed by averaging the XR-CT coefficients (resp. density values) contained in different rotating angular sectors of the core (Figure 4(a)). The radial distribution of XR-CT attenuation coefficients (resp. density values) is computed by averaging the coefficient (resp. density values) contained in running averaged annuli of constant thickness and increasing diameter (Figure 4(b)). The size of the angular sector and the thickness of the annulus used in moving average computations are user-defined parameters of the analysis. For the current work, angular sector is set to 10 degrees and annulus thickness is set to 1.5 mm (8 pixels), as a good compromise between computation time and resolution of the outputs.

Visualization-Plotting (Step 6)
Combining all the quantitative information automatically extracted on all the slices, synthesis depth-plots are automatically produced using the Borehole Geology and Geophysics package [27] relying on the graphical capabilities of the Generic Mapping Tools [28]. An example of summary figure for the section C9001C-11H4 (96.82 -98.32 m, curated core depth below seafloor) is given in Figure 5.  (d) Two color-coded measures of dispersion of the density values based on the distance between XR-CT-derived mean density ( ρ XRCT TOTAL ) and 50%-tile density ( ρ XRCT 50% ) and the width at half height maximum (WHH) of the histogram of distribution of density values (step 3; Table 2). Based on this example, ρ ρ − XRCT XRCT TOTAL 50% appears to be a better indicator of the degree of homogeneity of the core than WHH.
(e) Volume fraction of the voids/vugs (VV), damaged (DM) and non-damaged (ND) core fractions (step 4) clearly identifying voids related to clay expansion/contraction and core biscuit; (f) Color-coded core-quality estimator automatically derived from the abundance of VV and DM fractions (step 4, Table 3); in particular, this core-quality indicator constitutes a reliable and automatic source of information for assessing the impact of drilling/coring disturbance and partly filled core liner on vo- Open Journal of Geology lume-sensitive core logging data such as gamma ray attenuation density ( ρ GRA ), magnetic susceptibility (MS) and natural gamma radiation (NGR).
(g) Very-high vertical resolution (mm-scale) XR-CT derived average density of the non-damaged core fraction ( ρ XRCT ND ; step 4), thought to be free of all drilling/coring disturbance and, thus, best to estimate core bulk density ( ρ CORE BULK ); (h), (i) Angular and radial distributions of the density values (step 5 after calibration of attenuation values into density, step 2); When all slices of a given core section are combined, the unwrapped image of the angular distribution is similar to borehole wall images or a 360-degrees digital scan of the core surface where beddings and structural features can be recognized as sine waves ( Figure 5(h)).
This output can further be integrated with geographically oriented borehole wall images, such as wireline electrical micro-imaging or logging while drilling density images, if available, to back orient the cores [29]. Indeed, when non-ambiguous common features are recognized on both sets of images, cores can be geographically oriented and so provide unique source of information for detailed structural and paleomagnetic studies. Because such data were not available for hole C9001C, this possible use of the outputs is not discussed in the Application section. The radial distribution plot ( Figure 5(i)) provides information on radial variation that can help in diagnosing zones of bias in CT reconstruction (beam hardening characterized by radial attenuation), and/or zone of core expansion (low CT attenuation coefficient/low density in the center vs. high attenuation coefficient/density at the core periphery) or damage due to friction at the core liner interface (low CT attenuation coefficient/density at the core liner interface).
To take into account depth overlap between adjacent expanded cores (recovered core length > drilled interval), drilling (bit depth, length of cored interval) and curatorial information (length of recovered core) are integrated to produce a new depth scale, known as scaled core depth below seafloor. In this scale depth, a curated core length longer than its corresponding drilled interval is linearly compressed to fit the drilled interval. In case the curated core length is shorter than the drilling interval (incomplete core recovery or no core expansion) the curated length is not modified and equals the scaled curated core depth (IODP depth scale terminology, IODP-MI internal report, not published). Providing that affine transforms between depth scales are available, both scalar and array (image) outputs can be displayed in any reference depth.
Array (image) logs, such as 1) the 2D vertical and horizontal longitudinal sections of the 3D XR-CT images, 2) the histogram of distribution of the density values, and 3) the angular and radial distributions of the density values usually carry too much resolution/details to be displayed at the borehole scale (assuming standard A4 or letter display size). Scalar logs, possibly filtered, and color-coded semi-quantitative indicators of core quality and density measure dispersion are more suitable for such overview plots.

Applications
By themselves quantitative and semi-quantitative outputs of the previously described procedure for data reduction and quantification of the XR-CT images are valuable parts of the full data integration framework. In particular the above described outputs and their 2D graphical representations provide a valuable screening tool to further guide core sampling and/or more detailed analysis of the 3D XR-CT images with specialized software. Based on the data available for Hole C9001C, this section quickly presents four examples of application of the outputs of the presented procedure. Examples include 1) validation of sub-mm scale XR-CT derived core bulk density, 2) derivation of sub-mm scale porosity, 3) correction of volume sensitive MSCL measurements (NGR, MS and possibly ρ GRA ), and 4) core-seismic integration in the absence of well log (and acoustic data).  The shallow depth interval (Cores C9001C-1H to 9H, ~ 0 -76 m curated core depth below seafloor) where the formation is mostly unconsolidated is characterized by a continuous and significant amount of VV and DM (average VV + DM ~ 13%). Down to Core C99901C-34H (~ 302 m curated core depth below seafloor) where the hydraulic piston coring system was switched to the extended shoe coring system, core quality is overall good, with reduced core damaged and more localized core voids related to core decompression processes in more consolidated formation (VV + DM usually < 6%). Below ~ 333 m curated core depth below seafloor, overall core quality degrades and VV by itself mostly exceed 20% (partly filled core liner) due to erosion/surface friction and fracturing In any case, they provide a lower bound estimate of the core bulk density. In terms of methodology, even if ρ GRA suffers from localized artificial gaps due to early sampling of cores, the match between ρ GRA and ρ XRCT TOTAL provides an independent validation of the calibration of ρ XRCT N data by taking advantage of the more continuous and the higher vertical resolution nature of the ρ GRA data (~6800 measurements) in respect to the ρ MAD BULK data (77 measurements). In summary, ρ XRCT ND provides a ρ CORE BULK estimate of quality comparable to ρ MAD BULK , with a vertical resolution much higher than ρ GRA which value is highly dependent of core quality and which can possibly suffer from artificial sampling gaps. Note that even if ρ XRCT

X-Ray CT Derived Porosity
After having validated the sub-mm scale resolution ρ XRCT ND as best ρ CORE BULK estimate, a sub-mm scale porosity core log ( φ XRCT ) is derived by combining ρ XRCT ND with the interpolated grain (solid) density core log, ρ MADint SOLID , according to Equation (3): where ρ CORE FLUID is the averaged density of the fluid (pore) fraction of the core, assumed here to be 1.025 g•cm −3 (seawater). Computation of φ XRCT is technically simplified if 1) ρ MADint  Figure 9(a). ρ MAD SOLID and its spline interpolated version, ρ MADint SOLID , are plotted in Figure 9(b). Assuming the interpolation is valid, the resulting very-high vertical resolution φ XRCT is compared with φ MAD in Figure 9 ), the match between φ XRCT and φ MAD is generally within 2%. However, use of discrete MAD data both as input ( ρ MAD SOLID ) and output ( φ MAD ) does not provide a validation of the interpolated values. A complete independent validation would have required a separate data set, such as a neutron porosity logs, that were unfortunately not available for this wellbore. In any case, the gain in vertical resolution on one hand, and the simplicity of the procedure on the other hand makes this XR-CT derived porosity estimate valuable for other investigations, in particular to assess magnitude of "porosity rebound" on density estimates [30].

Volume Correction of Core Physical Properties Measurement
After proper calibration of the MSCL sensors, if the core liner is full, the MSCL measurement, MSCL M P , reflects the core property, CORE P , in a one-to-one relationship ( Figure 10(a)) as Equation (4): This simple relationship breaks down when the core liner is not completely full ( Figure 10(b)). In that case, MSCL M P include both contributions of the VV (ν) and core material fraction (CORE) according to the following general mixture relationship shown as Equation (5): where ν is the integrated volume fraction of VV over the MSCL spatial tool response, and ν P is the value of the given property, P, for the VV fraction.
This latter Equation (5) is obviously reduced into Equation (4)  response and tabulated value of ν P used in this work are given in Table 4.    Figure 11(d) and Figure 11(e)) core logs in case of partly filled core liner.
As an additional test, the volume correction procedure has been applied to ρ GRA assuming VV filled with air ( Figure 11(f)). Figure 11(f) confirms the good match between the raw ρ GRA (black) and ρ XRCT TOTAL (red, 0.1 m boxcar filtered) previously discussed. In general volume corrected ρ GRA (Figure 11  drilling disturbance free values, ρ CORE BULK , with a much higher vertical resolution (sub-mm for XR-CT vs. cm scale for MSCL), such development has not been conducted.

Seismic Modeling -Well Ties
Correlations between 1D physical property logs and 2D or 3D seismic reflection data are useful to provide geological information on the nature of seismic reflections and the lateral extent of geological layers [29]. The different resolution of core and seismic data is one major hurdles seismic interpreters face when tying both data sets, core data having a resolution that is at least two orders of magnitude greater [31] [32]. In the absence of sonic log data, [33] used the XR-CT derived density log ( ρ XRCT ND ) to compare three different approaches (binomial filtering, seismic wavelet filtering, and Discrete Wavelet Transform filtering (DWT)) to upscale the density log to the seismic scale. For Hole C9001C, best P. Gaillot et al. Open Journal of Geology ties with the seismic data were obtained when the density log was processed using the DWT, as this method best preserves the physical property changes representing the burial history of geological layers following their deposition.
The high-quality and high-resolution nature of the density log on one hand and the multi-scale characteristics embedded in a wavelet transform [34] making it sensitive to both cyclic (e.g. eustatic cycles) and non-cyclic components (e.g. local sediment inputs) of the signal on the other hand, made the seismic modeling successful, even in the absence of sonic log data ( Figure 12).

Conclusions
The presented methodological development for automatic data reduction and quantification of single energy XR-CT scanner images of cores allows 1) identification and quantification of voids/vugs (VV), damaged (DM) and non-damaged (ND) core fractions in terms of volume fraction, average density and mass contribution, 2) estimation of core quality based on abundance of VV and DM fractions, 3) characterization of local variation of core density in terms of radial and angular density distributions, and ultimately 4) derivation of a very-high resolution bulk density core log. Following standard approaches for a quantitative use of XR-CT images for estimating sediment bulk density determination [12] [16], the proposed procedure relies on core-to-core calibration, and also takes into account lithology-controlled mineralogical effects. Additional strength of the developed procedure is that it relies on the full integration of various data sets In addition to simple HU-conversion, the procedure automatically produces a histogram of distribution of the density values from which the ND fraction of the core can be isolated. As a result, a core bulk density estimate, ρ ρ CORE XRCT BULK ND , free of any drilling disturbance is provided. For slurry cores or cores damaged by the rotary core system, ρ XRCT TOTAL and ρ GRA should be disregarded in terms of absolute density estimates but can still be used for independent validation of the XR-CT derived densities ( ρ XRCT TOTAL and ρ XRCT ND ). As illustrated in this paper, the simplest and most direct usage of ρ XRCT ND is the derivation of a sub-mm scale resolution porosity, φ XRCT provided that ρ CORE SOLID and ρ CORE FLUID are known. In addition to provide semi-quantitative core quality indicators, quantitative estimates of VV and DM augmented by a knowledge of MSCL sensor responses can be used to apply correction to volume sensitive measurements (NGR, MS and possibly ρ GRA ) when the core liner is not entirely filled by cored material. Coupled or not with sonic data, the third use of the continuous ρ XRCT ND is the generation of seismic impedances/reflection coefficient to generate reflection coefficient series that are later used to build synthetic seismograms to assess seismic-to-well ties and for the seismic upscaling of physical property logs [33].
Finally, the azimuthal distribution of density/HU values can be used to reorient core if orientated borehole wall images are available and if non-ambiguous features can be recognized on both data sets.
In conclusion, the presented procedure and algorithms for automatic data reduction and quantification of single energy XR-CT images of cores are fast and efficient. They reduce data volume by ~ 90% and provide valuable outputs in a compact and directly usable form. The procedure provides an integrated framework maximizing the value of drilling/curation, core (VCD, MAD, MSCL), borehole (wireline, logging while drilling scalar and image logs) and seismic data, thus opening avenues for core-log-seismic integration and upscaling of physical properties log. Routine application of this workflow to large scale projects with XR-CT data, in particular IODP expeditions, is recommended.