Induced Gamma Ray Spectroscopy

Key Takeaways

• Neutron-capture versus inelastic scattering gamma ray spectra provide complementary elemental information at different formation depths and neutron energy ranges: inelastic scattering gamma rays (produced when fast neutrons interact with nuclei before slowing to thermal energies) carry information about carbon (TOC estimation), oxygen, and silicon content in the near-wellbore formation (within approximately 10 to 20 centimeters of the borehole wall) at a high neutron flux that allows fast measurement even in short logging intervals; capture gamma rays (produced when thermalized neutrons are captured by nuclei after the neutrons have slowed to near-room-temperature thermal velocities) penetrate deeper into the formation (20 to 30 centimeters) and include the full suite of capture cross-section elements, with hydrogen and chlorine captures being the strongest and providing information about formation water salinity and fluid content alongside the lithology elements; modern spectroscopy tools measure both spectra simultaneously, using the capture spectrum for stable elemental ratios that form the basis of mineralogy calculations and the inelastic spectrum for carbon and oxygen ratios that estimate hydrocarbon versus water in the pore space and TOC in organic-rich shales.
• Mineralogy quantification from spectroscopy is the primary application driving formation evaluation in complex lithologies including organic-rich shale gas and tight oil reservoirs, where the conventional triple-combo log suite (gamma ray, density, neutron porosity) cannot reliably distinguish between illite, chlorite, and kaolinite clays (which have different cation exchange capacities and thus different effects on resistivity interpretation), between dolomite and calcite (which have similar density and neutron responses but very different mechanical and geochemical properties), or between pyrite and clay (both of which are heavy minerals that reduce bulk density and affect neutron porosity in opposite directions): the spectroscopy-derived mineralogy resolves these ambiguities by providing direct element-based discrimination between mineral types, with silicon from quartz and feldspar, calcium from calcite and dolomite (with sulfur distinguishing calcium sulfate from calcium carbonate), iron from pyrite, siderite, and iron-rich clays, and titanium from rutile and ilmenite (accessory minerals that correlate with detrital clay input in clastic reservoirs); the mineralogy volumes are used in the total porosity and effective porosity calculation (by subtracting the clay-bound and structural water associated with each clay mineral from the total water volume) and in geomechanical property estimation (by calculating the rock brittleness index from quartz and carbonate content, which determines the hydraulic fracturability of the shale).
• Total organic carbon (TOC) estimation from induced gamma ray spectroscopy uses the carbon-oxygen (C/O) ratio from the inelastic scattering measurement to estimate the fraction of carbon in the formation that is associated with organic matter rather than carbonate minerals: in a carbonate-free formation, all carbon is organic (from kerogen and bitumen in the rock matrix), and the spectroscopy C/O ratio provides a direct organic carbon measure after correction for the carbon contribution from CO2 in the formation fluid; in formations with both carbonate and organic carbon, the calcium content (from the capture spectrum) is used to estimate the carbonate carbon fraction (CaCO3 contains one carbon per three oxygens), and the remaining carbon after subtracting the carbonate contribution is assigned to organic matter; spectroscopy-based TOC estimates are calibrated against core TOC values (from Rock-Eval pyrolysis or LECO combustion analysis) in the same well and applied to spectroscopy-only wells in the basin where core data is not available, providing a continuous TOC log that guides organic shale resource mapping and the identification of the highest-TOC intervals for horizontal well landing zone selection.
• Environmental corrections for induced gamma ray spectroscopy are simpler than for density and neutron porosity in some respects (the spectroscopy measurement is less sensitive to borehole diameter and fluid type) but more demanding in others (the neutron source activation of borehole and formation materials creates background gamma radiation that must be subtracted from the detected spectrum before elemental analysis): in large or irregular boreholes (washouts), the volume of borehole fluid between the tool and the formation increases the contribution of borehole elements (chlorine from NaCl in brine, silicon from silicate muds, barium from barite-weighted mud) to the detected spectrum and must be corrected using a borehole sigma (neutron capture cross-section of the borehole fluid) input determined from mud weight and composition; the gadolinium (Gd) correction is particularly important in wells drilled with muds containing gadolinium compounds as neutron absorbers (used in reservoir salt sections to limit neutron tool activation), because gadolinium's extremely high neutron capture cross-section creates a borehole contribution to the Gd yield that overwhelms the formation Gd signal; formation Gd content from spectroscopy (where not overwhelmed by borehole Gd) is used as a shale indicator in some reservoirs because Gd concentrates in clay minerals through its affinity for clay exchange sites.
• Simultaneous log interpretation using spectroscopy-derived mineralogy integrates the elemental composition from spectroscopy with the density, neutron, and resistivity measurements from the triple combo to provide a fully consistent multi-mineral petrophysical model that simultaneously satisfies all measured log responses using the mineral volumes as the fitting parameters: in a complex lithology with three to five mineral components (quartz, calcite, dolomite, illite, and pyrite, for example), the five unknown mineral volumes cannot be uniquely determined from three log responses (density, neutron, PEF) without additional constraints; the spectroscopy elemental ratios provide the additional constraints needed to over-determine the mineral system (more equations than unknowns), enabling a least-squares inversion that finds the mineral volumes most consistent with all log measurements simultaneously; the resulting mineral volumes are self-consistent by construction (they honor all logs simultaneously rather than reconciling logs in a sequential workflow) and are the basis for a petrophysical model that can be confidently used for reserve calculation, simulation model property population, and completion design in heterogeneous formations where single-measurement-based porosity and saturation calculations are unreliable.