Combining inelastic and capture gamma ray spectroscopy for determining formation elemental

Patent metadata

Field	Value
Publication number	US-10429540-B2
Application number	US-201213656287-A
Country	US
Kind code	B2
Filing date	Oct 19, 2012
Priority date	Dec 15, 2011
Publication date	Oct 1, 2019
Grant date	Oct 1, 2019

Abstract

A method for determining weight fractions of a plurality of elements in a subsurface formation penetrated by a wellbore includes determining weight fractions of a first plurality of elements using measurements of capture gamma rays made in the wellbore. The capture gamma rays result from bombardment of the formations with high energy neutrons. A weight fraction of a second plurality of elements is determined using measurements of inelastic gamma rays made in the wellbore resulting from bombardment of the formations with high energy neutrons. The weight fraction for the second plurality is determined by setting a transformation factor for those elements common to both the first and second plurality to result in a statistically equivalent weight fraction for the common elements as determined for the first plurality of elements. Weight fractions of elements in the second plurality not common to the first plurality of elements is determined using the set transformation factor.

First claim

What is claimed is:

1. A method for determining weight fractions of a plurality of elements in a subsurface formation penetrated by a wellbore, comprising:
in a computer, determining weight fractions of individual elements in a first plurality of elements using measurements of capture gamma rays made in the wellbore, the capture gamma rays resulting from bombardment of the formations with neutrons, and a capture transformation factor Fc;
in the computer, determining weight fractions of individual elements in a second plurality of elements using measurements of inelastic gamma rays made in the wellbore, the inelastic gamma rays resulting from bombardment of the formations with neutrons, and an inelastic transformation factor Fi;
in the computer, optimizing the inelastic transformation factor Fi by minimizing weighted average difference between inelastic and capture weight fractions for elements present in both the first plurality of elements and the second plurality of elements;
in the computer, determining weight fractions of individual elements that are present in the second plurality but that are not present in the first plurality of elements using the optimized inelastic transformation factor.

2. The method of claim 1 wherein the determining weight fractions of individual elements in the first plurality of elements comprises minimization of a closure function.

3. The method of claim 1 wherein the bombardment of the formation comprises emitting controlled duration bursts of neutrons from a pulsed neutron generator.

4. The method of claim 1 wherein the detecting gamma rays comprises activation thereby of a scintillation detector.

5. The method of claim 1 wherein the optimization of the inelastic transformation factor Fi is performed according to equation:
Σ j {( FiYi,j/Si,j−Wc,j ) 2 /(( FiσYi,j,/Si,j ) 2 +(σ Wc,j ) 2 )}

wherein j designates an individual element, Y designates gamma-ray yield, W designates weight fraction, S designates elemental sensitivity, a designates standard deviation, i designates inelastic gamma rays, and c designates capture gamma rays.

6. The method of claim 1 further comprising, in the computer, determining an enhanced transformation factor Fic by taking into account of measurements of both inelastic gamma rays and capture gamma rays.

7. The method of claim 2 wherein determining weight fractions of individual elements in the first plurality of elements comprises determining enhanced weight fractions thereof and recalculating a transformation factor for the first plurality of elements until differences between the determined weight fractions of the first plurality of elements and the enhanced weight fractions is at least one of minimized and fall below a selected threshold.

8. The method of claim 3 wherein the inelastic gamma rays are detected in a first selected length time interval occurring during and/or after at least one of the neutron bursts.

9. The method of claim 3 wherein the capture gamma rays are detected in a second selected length time interval after an end of at least one of the neutron bursts.

10. The method of claim 4 wherein the determining weight fractions of individual elements comprises analyzing numbers of gamma rays having particular energy levels by measuring an amplitude of light flashes emanating from the scintillation detector upon detection gamma rays.

11. The method of claim 6 , wherein Fic is determined by solving equation:
Fic {(Σ nXn ( Yc,n/Sc,n )}+Σ mXmWic,m= 1

wherein m designates a subset of individual elements measured by both inelastic gamma rays and capture gamma rays, n designates another subset of individual elements measured by capture gamma rays but are excluded in subset m, Y designates gamma-ray yield, W designates weight fraction, S designates elemental sensitivity, i designates inelastic gamma rays, and c designates capture gamma rays.

12. The method of claim 11 further comprising, in the computer, determining an enhanced weight fraction Wic by solving equation:
Wic,n=FicYc,n/Sc,n.

13. A method for well logging comprising:
moving a well logging instrument along a wellbore drilled through subsurface formations, the well logging instrument including a source of neutrons and at least one gamma ray detector, the gamma ray detector coupled to a spectral analyzer;
in a computer, determining weight fractions of individual elements in a first plurality of elements using measurements of capture gamma rays made in the wellbore, the capture gamma rays resulting from bombardment of the formations with neutrons, and a capture transformation factor Fc;
in the computer, determining weight fractions of individual elements in a second plurality of elements using measurements of inelastic gamma rays made in the wellbore, the inelastic gamma rays resulting from bombardment of the formations with neutrons, and an inelastic transformation factor Fi;
in the computer, optimizing the inelastic transformation factor Fi by minimizing weighted average difference between inelastic and capture weight fractions for elements present in both the first plurality of elements and the second plurality of elements;
in the computer, determining weight fractions of individual elements that are present in the second plurality but that are not present in the first plurality of elements using the optimized inelastic transformation factor.

14. The method of claim 13 wherein the determining weight fractions of individual elements in the first plurality of elements comprises minimization of a closure function.

15. The method of claim 13 wherein the moving the well logging instrument along the wellbore drilled through subsurface formations comprises extending and/or retracting an electrical cable having the well logging instrument coupled to an end thereof.

16. The method of claim 13 wherein the moving a well logging instrument along a wellbore drilled through subsurface formations comprises moving a drill string having the well logging instrument forming a part thereof along the wellbore.

17. The method of claim 13 wherein the bombardment of the formation comprises emitting controlled duration bursts of neutrons from a pulsed neutron generator.

18. The method of claim 13 wherein the detecting gamma rays comprises activation thereby of a scintillation detector.

19. The method of claim 14 wherein determining weight fractions of individual elements in the first plurality of elements comprises determining enhanced weight fractions thereof and recalculating a transformation factor for the first plurality of elements until differences between the determined weight fractions of the first plurality of elements and the enhanced weight fractions is at least one of minimized and fall below a selected threshold.

20. The method of claim 17 wherein the inelastic gamma rays are detected in a first selected length time interval occurring during and/or after at least one of the neutron bursts.

21. The method of claim 17 wherein the capture gamma rays are detected in a second selected length time interval occurring after an end of at least one of the neutron bursts.

22. The method of claim 18 wherein the determining weight fractions of elements comprises analyzing numbers of gamma rays having particular energy levels by measuring an amplitude of light flashes emanating from the scintillation detector upon detection gamma rays.

Assignees

• Schlumberger Technology Corp

Inventors

• Grau James A

Classifications

• G01V5/101 Primary

  and detecting the secondary Y-rays produced in the surrounding layers of the bore hole · CPC title