Borehole compensation during pulsed-neutron porosity logging
Methods, tools, and systems for determining porosity in an earth formation are disclosed. Neutrons are emitted into the formation to induce inelastic scattering gamma rays and thermal capture gamma rays in the formation. The induced gamma rays are detected at a proximal gamma detector and a far gamma detector, which are spaced at different axial distances from the neutron source. A measured proximal-to-far inelastic ratio (a ratio of inelastic scattering gammas detected at the proximal and far detector) and a proximal-to-far thermal capture ratio (a ratio of thermal capture gammas detected at the proximal and far detector) are determined and used to calculate the formation porosity. Techniques are disclosed for removing borehole and casing configuration effects from the measured proximal-to-far thermal capture ratio, leaving only porosity dependence.
1. A method of measuring a porosity of an earth formation traversed by a wellbore, the method comprising:
receiving data generated by a logging tool, wherein the logging tool comprises:
a neutron source configured to emit neutrons into the formation at an energy sufficient to induce inelastic scattering gamma rays and thermal capture gamma rays in the formation,
a proximal gamma detector spaced a first axial distance from the neutron source, and
a far gamma detector spaced a second axial distance from the neutron source, and
wherein the data indicates gamma rays detected at the proximal gamma detector and gamma rays detected at the far gamma detector,
from the data, determining a count of thermal capture gamma rays detected at the proximal gamma detector and a count of thermal capture gamma rays detected at the far gamma detector,
determining a proximal-to-far thermal capture ratio as a ratio of the count of thermal capture gamma rays detected at the proximal gamma detector to the count of thermal capture gamma rays detected at the far gamma detector,
determining a borehole thermal neutron capture cross section, and
using the proximal-to-far thermal capture ratio and the borehole thermal neutron capture cross section to determine the porosity.
2. The method of claim 1 , wherein determining the porosity comprises using the borehole thermal neutron capture cross section to compensate for borehole salinity.
3. The method of claim 1 , wherein the data indicating the gamma rays detected at the proximal gamma detector and the gamma rays detected at the far gamma detector comprises, for each detector, a time spectrum comprising a burst interval indicating gamma ray counts detected while the neutron source is emitting neutrons and a decay interval indicating gamma ray counts detected while the neutron source is not emitting neutrons.
4. The method of claim 3 , wherein determining the count of thermal capture gamma rays detected at the proximal gamma detector and the count of thermal capture gamma rays detected at the far gamma detector comprises, for each detector:
determining a count of thermal capture gamma rays detected during the decay interval at that detector, and
determining a count of thermal capture gamma rays detected during the burst interval at that detector.
5. The method of claim 4 , wherein determining a count of thermal capture gamma rays detected during the decay interval comprises integrating the time spectrum over the decay interval.
6. The method of claim 4 , wherein determining a count of thermal capture gamma rays detected during the burst interval comprises:
fitting a decay function to the decay interval of the time spectrum,
determining a borehole component and a formation component of the decay function over the decay interval,
convolving the borehole component and the formation component of the decay function over the burst interval, and
summing the convolved borehole component and the convolved formation component over the burst interval to determine the count of thermal capture gamma rays detected during the burst interval.
7. The method of claim 1 , wherein determining the borehole thermal neutron capture cross section comprises:
fitting a decay function to the decay interval of the time spectrum, and
determining a borehole component and a formation component of the decay function over the decay interval.
8. The method of claim 7 , wherein the decay function is a dual exponential function (DEF).
9. The method of claim 8 , wherein the dual exponential function (DEF) is expressed as:
DEF= A BH *exp(− t/τ BH )+ A F *exp(− t/τ F ),
where A BH and A F are coefficients, t is time, τ BH is a borehole decay constant, and τ F is a formation decay constant.
10. The method of claim 9 , wherein the borehole thermal neutron capture cross section is determined based on τ BH .
11. A system for measuring a porosity of an earth formation traversed by a wellbore, the system comprising:
a logging tool comprising:
a neutron source configured to emit neutrons into the formation at an energy sufficient to induce inelastic scattering gamma rays and thermal capture gamma rays in the formation,
a proximal gamma detector spaced a first axial distance from the neutron source, and
a far gamma detector spaced a second axial distance from the neutron source, and
a computer configured to:
receive data generated by the logging tool, wherein the data indicates gamma rays detected at the proximal gamma detector and gamma rays detected at the far gamma detector,
from the data, determine a count of thermal capture gamma rays detected at the proximal gamma detector and a count of thermal capture gamma rays detected at the far gamma detector,
determine a proximal-to-far thermal capture ratio as a ratio of the count of thermal capture gamma rays detected at the proximal gamma detector to the count of thermal capture gamma rays detected at the far gamma detector,
determine a borehole thermal neutron capture cross section, and
determine the porosity from the proximal-to-far thermal capture ratio and the borehole thermal neutron capture cross section.
12. The system of claim 11 , wherein determining the porosity comprises using the borehole thermal neutron capture cross section to compensate for borehole salinity.
13. The system of claim 11 , wherein the data indicating the gamma rays detected at the proximal gamma detector and the gamma rays detected at the far gamma detector comprises, for each detector, a time spectrum comprising a burst interval indicating gamma ray counts detected while the neutron source is emitting neutrons and a decay interval indicating gamma ray counts detected while the neutron source is not emitting neutrons.
14. The system of claim 13 , wherein determining the count of thermal capture gamma rays detected at the proximal gamma detector and the count of thermal capture gamma rays detected at the far gamma detector comprises, for each detector:
determining a count of thermal capture gamma rays detected during the decay interval at that detector, and
determining a count of thermal capture gamma rays detected during the burst interval at that detector.
15. The system of claim 14 , wherein determining a count of thermal capture gamma rays detected during the decay interval comprises integrating the time spectrum over the decay interval.
16. The system of claim 14 , wherein determining a count of thermal capture gamma rays detected during the burst interval comprises:
fitting a decay function to the decay interval of the time spectrum,
determining a borehole component and a formation component of the decay function over the decay interval,
convolving the borehole component and the formation component of the decay function over the burst interval, and
summing the convolved borehole component and the convolved formation component over the burst interval to determine the count of thermal capture gamma rays detected during the burst interval.
17. The system of claim 11 , wherein determining the borehole thermal neutron capture cross section comprises:
fitting a decay function to the decay interval of the time spectrum, and
determining a borehole component and a formation component of the decay function over the decay interval.
18. The system of claim 17 , wherein the decay function is a dual exponential function (DEF).
19. The system of claim 18 , wherein the dual exponential function (DEF) is expressed as:
DEF= A BH *exp(− t/τ BH )+ A F *exp(− t/τ F ),
where A BH and A F are coefficients, t is time, τ BH is a borehole decay constant, and τ F is a formation decay constant.
20. The system of claim 19 , wherein the borehole thermal neutron capture cross section is determined based on τ BH .