Quantitative Interpretation Of Pulsed Neutron Capture Logs In Thinly-Bedded Formations

Pulsed neutron capture (PNC) logs are commonly used for formation evaluation behind casing and to assess time-lapse variations of hydrocarbon pore volume. Because conventional interpretation methods for sigma logs assume homogeneous formations, errors may arise in thinly-bedded formations when appraising petrophysical properties of hydrocarbon-bearing beds. There exist no quantitative interpretation methods to account for shoulder-bed effects on PNC logs acquired in sandshale laminated reservoirs. Due to diffusion effects between dissimilar beds, sigma logs acquired in such formations do not obey linear mixing laws between the sigma responses of pure-sand and pure-shale end members of the sedimentary sequence. We introduce a new method to rapidly simulate sigma logs in thinly-bedded formations. The method makes use of late-time, thermal-neutron flux sensitivity functions (FSFs) to describe the contribution of multilayer formations to the measured capture cross section. It includes a correction procedure based on diffusion theory that adapts a homogeneous, base-case FSF to simulate the response of vertically heterogeneous formations. Benchmarking exercises indicate that our rapid simulation method yields differences smaller than 2% with respect to PNC logs simulated with rigorous Monte Carlo methods for a wide range of geometrical, petrophysical, and fluid properties. The second part of the paper develops an inversion method to reduce shoulder-bed effects on PNC logs and estimate layer-by-layer capture cross sections based on the previously defined FSFs. We successfully test the estimation method on synthetic examples that include a variety of bed-thickness configurations. Inversion consistently improves the vertical resolution and sigmadefinition of PNC logs across beds thinner than 1.5 ft. Testing on field data confirms the efficiency, reliability, and stability of the inversion procedure. Our fast simulation/interpretation algorithm inverts sigma logs over 100+ ft. intervals in less than 4 minutes of CPU time, and is therefore suitable for joint petrophysical interpretation with other open- and cased-hole logs.

[1]  T. Davis,et al.  Characterization and 3D reservoir modelling of fluvial sandstones of the Williams Fork Formation, Rulison Field, Piceance Basin, Colorado, USA , 2008 .

[2]  W. Hoyle,et al.  Quantitative Interpretation of Thermal Neutron Decay Time Logs: Part I. Fundamentals and Techniques , 1971 .

[3]  W. Hoyle,et al.  QUANTITATIVE INTERPRETATION OF THERMAL NEUTRON DECAY TIME LOGS. PART II. INTERPRETATION EXAMPLE, INTERPRETATION ACCURACY, AND TIME-LAPSE TECHNIQUE. , 1971 .

[4]  H. D. Scott,et al.  Computing thermal neutron decay time environmental effects using Monte Carlo techniques , 1986 .

[5]  A. H. Youmans,et al.  Neutron Lifetime, a New Nuclear Log , 1964 .

[6]  R. A. Haley Pulsed Neutron Capture Log Interpretation In Laminated Formations: A Dual-Exponential-Decay Model , 1995 .

[7]  Per Christian Hansen,et al.  REGULARIZATION TOOLS: A Matlab package for analysis and solution of discrete ill-posed problems , 1994, Numerical Algorithms.

[8]  D. Marquardt An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[9]  C. Torres‐Verdín,et al.  Linear iterative refinement method for the rapid simulation of borehole nuclear measurements: Part I — Vertical wells , 2010 .

[10]  C. Morris,et al.  Applications of Pulsed Neutron Capture Logs in Reservoir Management , 2005 .

[11]  H. D. Scott,et al.  SNUPAR: a nuclear parameter code for nuclear geophysics applications , 1989 .

[12]  S. J. Kimminau,et al.  The Design of Pulsed Neutron Reservoir Monitoring Programs , 1992 .

[13]  M. G. Portal,et al.  Improved pulsed neutron capture logging with slim carbon-oxygen tools: Methodology , 1995 .

[14]  W. Stacey Nuclear Reactor Physics , 2001 .