Search Documents

By E. A. Breitenbach

A computer program has been developed to afford rapid and complete quantitative log analysis for exploration and production decisions. The computation consists of automatic selection of tops and bottoms of porous intervals from the digitized data, and then point-by-point calculations within each selected interval. Nearly all log types can be analyzed. This paper presents the calculation techniques found to be appticable to machine evaluation and gives examples of their use. INTRODUCTION As the quantitative interpretation of well logs entails repetitive use of charts and equations, it is natural that digital computer programs would be written. A Universal Log Interpretation Computer Program (ULICP) has been developed to afiord rapid and complete analysis for exploration and production decisions. Application of this program enables the log analyst to: (I) apply rapid quantitative log analysis for exploration and production decisions on either a field, a well or on a single zone basis; (2) apply concepts of interpretation requiring detailed numerical analysis; (3) analyze all porous intervals on each log, rather than a few selected zones, and complete the analysis in less time than previously required. Computations are reported for every digitized point in each zone; (4) use empirical techniques applicable to a given area as an integral part of the computation; (5) experiment with the empirical coefficients and exponents in the interpretation equations to find the best possible solution; and (6) automatically plot both the original and computed data to scales analogous to the field prints. The primary utility of a log interpretation program stems from its ability to do an overwhelming amount of work with very little man-power. If a calculation procedure or thought process can be formalized to the extent that step-by-step logic can be written, a computer program can be developed that follows this logic. The major question is one of economics. A feasibility study of the costs for such a program indicates that digital processing is economical primarily as a means of increasing the productivity of the log analyst. In effect, the probability of missing productive intervals in any well, because of a lack of time to do detailed calculations, is greatly reduced. The nucleus of ULICP is programmed to compute large sections of a suite of logs by selecting zones automatically and then performing all pertinent computations on a data point by data point basis within each zone.'.' An entire suite of logs can be processed in this manner with very little manual intervention. Sufficient programming logic is available so that each log analyst can request computations pertinent to his area. These requests are made by simple additions or deletions of information on the input header cards. Hence, the log analyst is always in complete control of the computation process. The evaluation of a suite of logs requires pre-editing, digitization, computation, presentation of results and interpretation. Work by the log analyst is reqilired only in pre-editing and interpretation. Thus, he is allowed more time for comprehensive interpretation, rather than calculation. For continuity, the discussion of ULICP is organized sequentially: pre-editing, digitization, computation and presentation of results. DISCUSSION PRE-EDITING The extent of pre-editing prior to computation is dependent on the format of the original data. For analog prints, it requires inspection, correlation and editing of the logs, plus the entering of required data an special forms. For digitized data such as magnetic tape field recordings, only the special forms are necessary. The process for analog prints will be given here. The inspection and correlation process involves the selection of sections for digitization and the correlation of the traces to a common depth. To decrease the cost of digitization, traces that exhibit large variations in shale formations can be redrawn to a non-zero baseline. The next step is to enter all pertinent data on the input header cards. The header cards presently used for ULICP are presented as Figs. 1 through 6. Cards 1, 2, 3 and 4 (Fig. 1) are defined as the Main Header Cards. They describe the particular well for output identification and give basic information. Cards 5 through 15 (Figs. 2 through 6) are defined as the Block Header Cards. As such, they define the log types and the interpretation parameters for the block of data immediately following. Cards 1 through 10 are required for every computer run. Cards 11 through 14 pertain only to nuclear logs. Card 11 is required if any nuclear log is supplied, and cards 12 through 14 are required only when gamma-ray, density or neutron log data, respectively, are supplied. Card 15 (Fig. 6) must be supplied when cross-plot calculations are required. Either two- or three-component

Jan 1, 1967

Jan 1, 1969

By P. E. Baker

An improved method of logging formation density has been developed in which the formation is bombarded with a collimuted beam of gamma rays. By means of a scintillation detector and pulse height discriminator, the gamma-ray energy band is accepted and recorded which corresponds to deepest penetration into the formation. Laboratory tests on a field tool revealed no borehole diameter effect for smooth holes and no effect of formation chemistry, except insofar as chemistry affects density. In extensive field tests, the density log has exhibited satisfactory agree-ment with core measurements and has correlated accurately with other logs. INTRODUCTION This paper describes an improved method of logging formation density by means of back-scattered gamma rays to penetrate the formation together with energy discrimination to select only that energy (or frequency) range corresponding to deepest penetration into the forma- tion. This method of density logging was developed at the La Habra laboratory of California Research Corp.* Field instruments have been built and tested by McCullough Tool Co. and commercial service has begun in some areas. All of the data discussed here were obtained for the original field tool. The density log can be a valuable adjunct to existing logs because of its use in determining porosity quantitatively. It is useful also in connection with seismic and gravity surveys where formation density information is of interest. PHYSICAL PRINCIPLES The intensity, I, of a beam of gamma rays after traversing a distance, x, in an absorbing medium is given by I = loe -ax where I, is the intensity of the beam entering the medium and a the attenuation constant or macroscopic cross section of the medium for gamma-ray absorption. The macroscopic cross section depends on the number and types of atoms per unit volume and the gamma-ray energy. In brief, it is the sum of the atomic cross sections of all atoms per cubic centimeter. The density of the material is the sum of the masses of all atoms per cubic centimeter. Thus, the density and the macroscopic cross section for gam-ma-ray absorption vary in propor- tion to the number of atoms per cubic centimeter. The density of the material increases, of course, if the atomic masses increase without a change in the number of atoms per cubic centimeter. When the density increases in this latter way, the macroscopic cross section increases also, but disproportionately. This results from the fact that part of the gamma-ray absorption arises from the photoelectric effect and the pair-production effect**, which are proportional to the sixth and to the second power, respectively, of the atomic number of the absorber. The remainder of the absorption results from the Compton effect, and is directly proportional to the atomic number. The gamma rays that undergo Compton scattering are not immediately absorbed. Since these rays are scattered in all directions, a fraction will emerge from the face of the medium at which they entered. If a gamma-ray shield is placed between the source and detector, most gamma rays reaching the detector will be those scattered out of the medium. The intensity will depend on absorption and scattering in the medium, and so indicates the density of the medium. For gamma rays for which Compton scattering prevails, the intensity

Jan 1, 1958

By A. E. Worthington, R. F. Meldau

An analog computer study was made of the reduction of the magnitude of the SP by formation geometry and resistivity effects. Practical field .situations were studied. The analog data were converted to the form of departure charts. These charts may he used to determine the reduction in the SP due to the above eflects pith good precision Quantitative interpretation of the SP log in terms of formation water resistivity is presently based on the premise that the observed SP is equal to the electrochemical potential developed across a shale between the mud filtrate and the formation water. This procedure inherently assumes that, first, no electro kinetic (streaming) potential is present, and, second, that the SP is not reduced by resistivity or geometry effects in the formation. When, however, such a reduction of the observed SP does occur, one may write, SP = b,E, b, E, (1) where SP is the observed self-potential opposite the center of the bed under consideration. In this case, E, is the electrochemical potential, E, is the electrokinetic potential, b, is the electrochemical reduction factor, and h, is the electrokinetic reduction factor. The first assumption above is equivalent to assuming E, = O in Eq. 1 and the second assumption is equivalent to h, = b, = 1. Examination of electric logs shows that for many formations of interest the observed SP is reduced below its maximum value. The object of this work is to aid quantitative SP interpretation in at least some of these situations by obtaining values for b, and b, in Eq. 1 as functions of realistic combinations of bed thickness, depth of invasion, and formation resistivity. The only previously published work in this field known to the authors is by H. G. Doll' in 1948. He evaluated b, under conditions permitting exact analytical solution. His work clearly established the importance and qualitative nature of formation geometry and resistivity effects upon the SP, but the necessity of obtaining analytical solutions restricted his study to situations of limited practical interest. In order to extend such studies to practical combinations of formation parameters which do not lend themselves to feasible analytical solutions, the results reported in this paper were obtained by use of a special analog computer. THE ANALOG MODEL The analog model used for this study was the Guyed Analog Computer which originally was designed and built for resistivity log departure curve studies. It consists of a two-dimensional mesh of resistors, with the vertical mesh corresponding linearly to vertical distance along the borehole, and the horizontal mesh corresponding to radial distance from the borehole axis. By proper choice of horizontal resistor values, a three-dimensional, axially symmetric, borehole-formation problem can be modeled with this two-dimensional array. These and other details of the model are contained in a paper by H. Guyod'. The geometry used for the SP studies with the analog is depicted in Fig. 1. This configuration has a horizontal plane of symmetry, and, hence, the model, and the following discussion, need only to be concerned with one quadrant such as shown in Fig. 2. A further simplification arises by considering that both the natural potential sources and those used in the model have fixed internal impedance. Therefore, the current flows produced by these sources can be separated and studied independently ( i.e., the superposition theorem can be applied). SP = SP, i SP, = h,E, + b,E, ,(7)

By S. E. Szasz, E. J. Moore, B. F. Whitney

An accurate value of formation water resistivity R, is essential in calculating formation porosity and fluid saturation from electrical well logs. In the cases where R, has not been measured directly, it must be obtained from other data, e.g., the SP curve. This paper deals with another approach: how to calculate R, from the chemical analysis of the formation water. INTRODUCTION It is known that the resistivity of aqueous solutions of pure salts depends on their concentration and on the temperature; the concentrations are given in MPL (mg of solute per liter of solution), or sometimes in ppm (mg of solute per kg of solution): MPL = ppm X specific gravity. Values for different pure salts are available in the literature, but not for solutions of mixtures which are of practical interest. The major component of the dissolved material in almost all formation waters being sodium chloride, it is customary to express the resistivity of formation waters in terms of equivalent sodium chloride concentration, i.e., the concentration of a solution of pure NaCl which has the same resistivity at a given temperature as that of the formation water under consideration. Thus, the problem of calculating R, from the chemical anaylsis can also be stated as how to convert the other constituents of the solute into equivalent NaCl concentration. Salts dissolved in water are at least partly dissociated into ions, and do not conserve their identity. If known amounts of several salts are dissolved in water, the solution does not necessarily contain the same salts in the original proportion, but perhaps some other combination of the ions, along with free ions in solution. This is why the chemical analysis of formation waters is often given in terms of ions, as if all dissolved salts were completely dissociated. Our problem then boils down to how to convert the concentrations of the various ions to equivalent concentrations of Na' and C1-. Dunlap and Hawthorne' have proposed to convert the concentration of all other ions to equivalent Na' and C1-concentrations by means of constant multipliers; e.g., 0.95 for Ca"; 2.0 for Mg"; 0.27 for HCO 3-; 0.5 for SO, -, etc. Their factors were based on measurements made at 68F on 26 formation water samples from the Texas Gulf Coast, ranging in concentration from 1,500 to 75,000 ppm. The Dunlap method is widely used in electric log interpretation, and is often extrapolated beyond its original concentration range. A comparison of R, values calculated by this method and values actually measured on formation water samples has shown large discrepancies, especially at higher concentrations. Therefore, two new methods were developed at Sinclair Oil Corp.'s Tulsa Research Center to calculate equivalent sodium chloride concentration from the chemical analysis of formation water samples. FUNDAMENTAL CONSIDERATIONS The resistivity of a solution, or its reciprocal the conductivity, at a given temperature is determined by the charge, concentration and mobility of the ions actually present. Monovalent ions such as Na' or C1- always carry the same charge. Compounds of polyvalent ions, however. may show incomplete dissociation, e.g., NaSO; + Na' instead of SO,-- + 2Na'. This happens especially in more concentrated solutions. Only very dilute solutions are completely dissociated, as assumed in the chemical analysis report. At higher concentration, the degree of dissociation depends not only on the nature and concentration of the particular salt under consideration but also on the nature and concentrations of the other solutes. Mobility of the ions depends on the viscosity of the solution. It also depends on the degree of hydration of the ions, which in turn is a function of the nature and the charge of the ions and also of the amount of free water available per ion, i.e., the total ionic concentration. The net effect is that the conductivity increases slower than proportional to the concentration, even if a solution contains only one salt such as NaC1, and is different for different salts (Fig. 1). Conductivity can even decline with a further increase in concentration, e.g., if additional salt is little dissociated but ties up some of the free water and/or causes an increase in viscosity. In solutions containing more than one salt, the contribution of one salt to the total conductivity depends not only on the fractional concentration of this same salt, but also on the concentration of all other solutes. A perfect method would give the conductivity or resistivity of a solution as a function of the concentrations of all solutes present. This is so complicated as to be impractical, and a simpler method must be found which is of acceptable accuracy. The Dunlap method, on the other hand. is too simple because it askmes that at any concentration the contrih-

Jan 1, 1967

By J. L. Martin, W. M. Campbell

A new electric logging method, called displacement logging, often gives a direct indication of the presence of mobile hydrocarbons in hydrocarbon-bearing formations. This method is based on the detection of a bank of connate water which is formed between the invading mud filtrate and virgin formation fluids in these cases. The existence of this phenomenon has beer, verified by laboratory studies and by actual field logs taken in known oil sands. INTRODUCTION During the past decade increasing use has been made of conventional electric logs to determine the amount of oil or gas in formations traversed by a borehole. This determination requires knowledge of several factors, such as true resistivity of the formation, formation factor, and water resistivity. Sometimes all of the factors necessary to make a reliable saturation calculation are not available. Also, the computations and corrections which are necessary in such a study often require the services of an experienced log analyst. It has long been the desire of exploration people to have a logging device which will locate hydrocarbon-bearing formations directly without requiring complex interpretation methods. The displacement log, while not completely fulfilling this desire, appears to be a Step in that direction. Basis for the Device The basis for this device stems from a phenomenon which occurs when mud filtrate invades a formation which contains displaceable hydrocarbons. As the mud filtrate invades the formation, most of the oil or gas is pushed out of the invaded zone, leaving residual oil. However, in addition to displacing most of the oil or gas, the connate water is also swept out, since the originally irreducible water is mobile in the presence of a water phase.1,2,3 Tonsequently, this irreducible water piles up as a bank between the virgin oil zone and the invading filtrate. If the mud filtrate is of higher resistivity than the formation water— as is commonly the case—the resulting annular zone, containing the bank of formation water, will be of lower resistivity than either the virgin oil-bearing formation resistivity or the invaded zone resistivity. This resistivity configuration, going from high to low to high again as the distance from the borehole increases, is found only in those zones which contain displaceable hydrocarbons. A device which would reliably detect this low resistivity zone would give a direct indication of mobile hydrocarbons. Illustrations This phenomenon is illustrated in Fig. 1 and 2. Fig. 1 shows a sketch of a zone containing salty connate water plus displaceable hydrocarbons, along with graphs showing fluid content and electrical resistivity vs distance from the borehole, before invasion occurs. Fig. 2 shows this zone, together with fluid content and resistivity graphs. after invasion by a fresh water mud. The resistivity pattern on moving away from the borehole is high: low: high. The absence of the low zone phenomenon in water zones is illustrated in Fig. 3 and 4. Fig. 3 shows a sketch of a zone containing salty formation water and residual oil, along with graphs showing fluid content and electrical resistivity vs distance from the borehole, before invasion occurs. Fig. 4 shows this water zone, together with fluid content and resistivity graphs, after invasion bv a fresh water mud. In this case the resistivity pattern on moving away from the borehole is high: low since there is no change in water saturation at the head ol

Jan 1, 1956

By J. D. Milburn, H. J. Hill

In quantitative interpretation of electrical logs the presence of clay minerals introduces an additional variable which further complicates an already complex problern. Although recognizing the difficulties introduced as a result of the heterogeneity of natural sediments and despite the present incomplete state of knowledge regarding electrochemical behavior of shales, dissemirzated clay minerals and concentrated electrolytes, it was felt that useful empirical correlations might be obtained from experimental investigation. Six typical sandstone formations, having a wide variety of petrophysical properties, were selected for the study. Approximately 45 samples from each formation were selected to satisfactorily represent the range of pore size distribution within the particular formation, As a matter of general interest, four limestone formations were also included in the investigation. Previously proposed equations relating to resistivity, SP and interrelationship of the two phenomena have, where possible, been tested with data obtained in this investigation. These equations do not satisfactorily describe experimental behavior of samples through all degrees of shaliness or throughout the range of brine solution resistivities normally encountered in logging practice. An empirical equation has been developed which quantitatively relates formation resistivity factor to satirrnting solution resistivity, porosity, and "effective clay content." This relation is indicated to be uniformly applicable to clean or shaly reservoir rocks. It is shown that both the SP and resistivity phenomena of shaly samples are related to the sample cation exchange capacity per unit pore volume. The independent chemical determination of this parameter is thus a means of determining the "effective clay content" of samples. Some implications regarding theory and electric log interpretation of shaly sands are discussed. INTRODUCTION The use of electrical resistivity logs as a means for estimating formation porosity is based upon the original work of Archie1. It was shown by empirical methods that the resistivity of a rock, when 100 per cent saturated with brine solution, could be related to rock porosity through the equation, F = Ro/ RW = f-m.......(1) where F, the formation resistivity factor, is defined as the ratio of the resistivity of the 100 per cent water saturated rock, Ro, to the saturating solution resistivity, Rw; f is the ratio of pore volume to bulk volume; and m is the slope of the average line when the data are plotted on log-log paper. This relation was established for relatively shale-free formations containing water having high salt concentrations, and under these conditions F may be considered constant for a given rock sample. Applications of Eq. 1 to interpretation of logs of shale-free formations has been quite successful when satisfactory formation water resistivity data could be obtained. The most readily available means of estimating formation water resistivity is the self-potential log. Conventional methods of SP log interpretation involve correct.

Jan 1, 1957

By G. H. F. Gardner

The velocity and attenuation of elastic waves in sandstones were measured as a function of both pressure and fluid saturation. A large change occurs in these quantities if water is added and the rock is not compressed, but the change is small if the rock is subjected to a large overburden pressure. Measurements were made by vibrating cylindrical samples in both the extensional and torsional modes at frequencies up to 30,000 cycles/sec. Formulas were derived which enable the attenuation of dilatational waves in dry rocks to be deduced from the data. Similar experimental methods were used to investigate the properties of unconsolidated sands. Velocities were found to vary with the 1/4 power of the overburden pressure and attenuations to decrease with the 1/6 power. The effects of grain size, amplitude and fluid saturation were studied. Formulas by which the effects produced by a jacket around the sample may be calculated were derived. The practical application of these results to formation valuation is discussed. INTRODUCTION The attenuation of elastic waves in the earth has been of interest to the seismologist and geophysicist for many years, but only recently to the petroleum engineer. Engineering interest has been brought about by the success of velocity logging devices, for it is possible by modification of these instruments to measure the attenuation of sound waves in addition to their velocity and, hence, deduce the mobility of formation fluids as well as the porosities of the rocks which contain them. The main problem is to decide whether field measurements can be made with sufficient accuracy to be of practical use. This problem can only be solved after we know the magnitude of the attenuations which are typical of the earth at various depths. The logarithmic decrement of a fluid-saturated rock is the sum of a "sloshing" decrement and a "jostling" decrement, the former caused by the mobility of the fluid contained within the rock and the latter by the granular framework of the rock. Sloshing decrements can be calculated' using Biot's theory, but the jostling losses are less well understood. The present paper reports an experimental investigation of jostling losses in consolidated and uncon- solidated sands, particularly with respect to the effect of overburden pressure and fluid saturation. Born' showed that the decrement of a sandstone may increase dramatically when only a few per cent by weight of distilled water is added, and that the additional loss is proportional to the frequency of vibration. His measurements were made with no compressive stress on the framework of the rock. M. Gondouin3 investigated similar phenomena for fluid-saturated plasters but also did not compress the samples. In the present paper it is shown that compression of the framework reduces this effect, so that at depth the jostling decrement of a sandstone may be expected to be almost independent of fluid saturation and frequency. Decrements for many sedimentary rocks have been given by Volarovich,4 but all for the state of zero overburden pressure. Anomalously low velocities have been logged in shallow unconsolidated gas sands. Results of the present investigation confirm that these velocities are not caused by correspondingly high attenuations, because the jostling decrement in a packing of sand grains is small and much less than in a consolidated sandstone at the same depth. Velocities in sands have been measured by Tsareva5 and by Hardin6 as a function of pressure, but the corresponding decrements do not appear to have been measured previously. The widely used "resonant bar method" of measuring velocities and decrements was employed. Comments on variations of this technique have recently been published by McSkimmin.7 The main novelty of the present technique was the application of pressure to the samples. It was found possible to do this by placing the apparatus inside a pressure vessel, provided the conditions leading to large additional losses were avoided. These conditions are discussed below. EXPERIMENTAL TECHNIQUE Cylindrical samples were caused to vibrate in both the extensional and torsional mode of vibration and the amplitude of vibration was measured as a function of frequency in the neighborhood of a resonant frequency. The resonant frequency, fr, is related to the corresponding elastic modulus by the formulas where E and N are Young's modulus and the modulus of rigidity, p is the density of the sample, and A the wavelength of the vibration.

Jan 1, 1965

By M. H. Waxman, L. J. M. Smits

A simple physical model was used to develop an equation that relates the electrical conductivity of a water-saturatedshaly sand to the water conductivity and the cation-exchange capacity per unit pore volume of the rock. This equation fits both the experimental data of Hill and Milburn and data obtained recently on selected shaly sands with a wide range of cation-exchange capacities. This model was extended to cases where both oil and water are present in the shaly sand. This results in an additional expression, relating the resistivity ratio to water saturation, water conductivity and cation-exchange capacity per unit pore volume. The effect of shale content on the resistivity index - water saturation function is demonstrated by several numerical examples. INTRODUCTION A principal aim of well logging is to provide quantitative information concerning porosity and oil saturation of the permeable formations penetrated by the borehole. For clean sands, the relationships between measured physical quantities and porosity or saturation are well known. However, the presence of clay minerals greatly complicates log interpretation, particularly the electrical resistivity and SP logs, and considerably affects evaluation of hydrocarbon-bearing formations. The conductance and electrochemical behavior of shaly sands and their relation to log interpretation have been studied by many workers. wylliel and Lunch2 reviewed this work in some detail. Virtually all laboratory measurements of electrical resistivity and electrochemical potential of shaly sands published to date are the work of Hill and Milburn.3 Their measurements were made on about 300 cores covering a large variety of sedimentary rocks, and a wide range of equilibrating NaCl solution concentrations. Hill and Milburn described their conductivity data by an empirical equation in which the shaly sand conductivity Co was expressed as a function of the solution conductivity C and two parameters b and F01 . The quantity b was shown to be a measure of the effective clay content of the rock, being approximately proportional to the cation-exchange capacity of the rock divided by its pore volume. The latter ratio is designated as 9, in this paper and has the dimensions meq/ml or equiv/liter. Qv is identical with the term representing the concentration of fixed charges in the Meyer-Sieved and Teorel15 theory of permselective membrane behavior. F01 is a formation resistivity factor referred to a hypothetical equilibrating solution resistivity of 0.01 ohm m at 25C,* where clay effects presumably are minimized. F01 was correlated to porosity by an Archie-type equation.6 The Hill-Milburn equation describes their data with a standard deviation of approximately 1 percent and a maximum deviation of + 10 percent. Shaly sands behave as permselective cation-exchange membranes, their electrochemical efficiencies increasing with increasing clay content. The electrochemical potential data of Hill and Milburn were expressed graphically, and demonstrate that the membrane efficiency (or cation transport number) of these sands is a function only of the b value (i.e., Qv) and the respective salt concentrations of the two solutions forming the liquid junction in the sand. The diffusion potentials are not dependent on F01 or any parameter relating to the porosity or pore geometry of the rock. The Hill-Milburn resistivity equation correctly predicts a decreasing sand conductivity Co with decreasing solution conductivity C?. However, at some low value of C,, the calculated Co-C, function goes through a minimum; with further decrease in C?, the predicted sand conductivity increases sharply. As pointed out by Hill and Milburn, an increasing value of Co with decreasing C, is physically meaningless. Since this occurs below the range of practical values of C,, the usefulness of the empirical equation is not affected when

Jan 1, 1969

By B. G. Agnew

In evaluating fracture treatments, the need to answer such questions as "What zone or zones were actually treated?" and "What was the vertical extent of the treatment" is necessary, since determining effectiveness of the fracture treatwrent depends on knowledge of the reservoir portions stimulated and the vertical extent of the fracture system or systems. Injection of hot or cold fluids during fracturing operations will transfer heat to the surrounding formations and fracture faces mainly by heat conduction. Once injection ceases, temperature anomalies will develop opposite the zone fractured because the rate of temperature decay is less than that opposite zones heated or cooled by flow inside the wellbore or along the cement sheath. Thus, temperature surveys can he used to determine where fractures were generated outside the wellbore and to reasonahly estimate the vertical extent of the fractrires. This technique has been used successfully to evaluate over 344 fracture treatments at depths ranging from 1,500 to 15,700 ft. Examples of actual temperature surveys are presented to show how this method of fracture treatment evaluation can be used effectively and beneficially. INTRODUCTION Stimulation of producing or injection wells by fracturing is commonplace in the industry today. In evaluating these treatments, two basic questions which arise are "What zone or zones were actually stimulated?" and "What was the vertical extent of the treatment?". Knowledge of the reservoir portions stimulated and vertical extent of the fracture system or systems is vital to effect efficient and economical well completions and to assure maximum recovery from reservoirs. Not all fracture treatments are successful and many times the desired or anticipated results are not achieved. When this occurs, it is in most cases the engineer's job to determine why. The answer to this problem forms the primary basis for deciding whether to spend more money for development drilling and/or completion attempts or to forego additional expenditures to improve the particular zone. Also, knowledge of the reservoir portion actually stimulated is invaluable in planning future workovers. This is especially true in wells producing from reservoirs containing multiple porosity stringers or zones. For various reasons, all zones within a reservoir may not be perforated on initial completion with the thought of developing remaining zones through future workover operations. Because particular zones are perforated and fracturing fluid pumped down the wellbore, it cannot be assumed that all perforated zones are stimulated, and in turn drained. Occasionally, unperforated zones which are planned for future development are actually stimulated on initial completion. eliminating the need for future recompletion attempts or all zones perforated initially are not stimulated, thus requiring additional completion work to insure that these zones are adequately drained. In most cases, it is costly. if not impossible, to determine with assurance why a fracture treatment is unsuccessful after the job has been completed and the results tested. Further, unpredictable upper and lower fracture barriers. as well as possible borehole communications (cement channels), do not permit it to be taken for granted that the formation opposite the perforations is the only zone receiving the fracture treatment. Thus. evaluation of fracture treatments is necessary for efficient and prudent well completion practices and should be done simultaneously with the fracture job. Several methods have been used to locate fractured zones, all of which are based on detection of radioactive tracer materials added to either the fracturing fluid or the propping agent. These methods have been helpful in increasing knowledge of fracturing operations; however, they are seriously limited because of cost and inability to detect accurately radioactive tracer material more than a few inches away from the borehole. To develop a better diagnostic tool to pinpoint which zone or zones had actually been fractured, temperature surveys were run after fracturing with heated and cooled fluids. Significant temperature anomalies were found to exist opposite zones suspected of receiving the fracture treatment. From the analysis and examples to be shown, it can be seen that temperature surveys run in a stabilized wellbore several hours after a fracture treatment using hot or cold fluid will yield temperature anomalies which show the actual reservoir portion stimulated, the vertical extent of the fracture system or systems and a qualitative indication of the portion of fracture-fluid volume entering a given depth interval. The experience and illustrations used in this paper are based on evaluation of over 344 fracture treatments during the past three years using this diagnostic tool at depths ranging from 1.500 to 15,700 ft and fracture volumes ranging from 2.000 to 110.000 gal. Based on analyses of these temperature surveys, analyses of fracture gradients and evaluation of several open-hole packer impressions in the West Texas area,' vertical fracturing is the predominant mechanism and horizontal fractures occurred frequently, if at all.

Jan 1, 1967

By R. H. Widmyer, G. M. Wood

This investigation was prompted by the absence of a simple, accurate method to predict optimum second stage separation pressure in three stage separation operations. Three stage separation as used here includes two separators plus the stock tank. The first stage pressure is usually fixed by gas sales conditions and the third stage is at atmospheric pressure. Thus, the second stage pressure is the only one that may be varied in-discriminantly during production operations. No variations in temperature were considered in this study, even though it will vary among systems, for a change in temperature has little effect on the optimum pressure. Optimum pressures may be determined from flash calculations, but these are tedious unless a computer is available. An early attempt to develop a simpler method&apos; yielded an equation which gave the optimum pressure only as a function of initial and final pressure. Experience, however, confirms theoretical relationships which show that composition as well as pressure and temperature must be a factor. This indicates that the basic relationship should be an equation of the form, where P, is second stage pressure, P, is first stage pressure, A is a function of P. (stock-tank pressure) and the system composition, and n is a function of the number of stages and the system composition. Fig. 1 shows the basic variation of P2 with P for a fixed system composition. A series of such figures clearly showed the effect of composition. The peak of each curve reflects that second stage pressure at which the stock-tank recovery (gal/MMcf) residue gas is a maximum value. DEVELOPMENT OF NEW CORRELATION Wellstreams having 17 different analyses, ranging from low gravity crude to condensate, were used to determine the effect of composition. Each wellstream was flashed through three stages to determine the optimum second stage pressure for each first stage pressure. Equilibrium vaporization ratios proposed by Latham and Campbell were used throughout. The heptane plus fraction, for simplicity, was assumed to have the properties of octane. The results were subjected to an empirical analysis in an attempt to develop a correlation- that would be (1) simple, (2) accurate, and (3) eliminate the need for any flash vaporization calculations. This analysis yielded two equations, one for wellstreams having a specific gravity greater than one (air = 1.0) and one for those less than one. For crude streams (specific gravFor condensate (specific gravity < 1.0),

By H. Guyod

The response of the Guard Electrode sonde (LL3) and the Laterolog (LL7) can be mathematically computed only for cases that do not represent realistic conditions. A resistance network analogue that simulates a number of typical—although idealized—formations was used instead to obtain, for both devices, a large number of resistivity logs and current distribution maps. From these data, charts were prepared for estimating true resistivities under a number of conditions. Some of these charts and current distribution maps are presented. INTRODUCTION Two main benefits are expected from a resistivity log: (1) obtaining geological information for mapping and other purposes, and (2) obtaining petrophysical data, based on true rock resistivity. The geological information can be derived from any resistivity curve having sufficient vertical detail, but true rock resistivities frequently are difficult to determine with desired accuracy because we do not have sufficient interpretation data. This is particularly true for the records given by the Laterolog tools. The object of this paper is to throw some light on the factors affecting the responses of these tools so that better resistivity data can be obtained from Laterologs. LATEROLOGS Laterolog measuring systems are of two types: the LL3 sonde, frequently called "Guard", and the LL7 sonde. Both measure the resistance between a single exploring electrode placed in the borehole and a faraway return electrode whose influence will be assumed here to be negligible. Although the principles of these tools have already been outlined,1-3 it is desirable to review them briefly in this article. LL3 Fig. 1, left side, represents a long cylindrical, vertical current electrode placed in a medium of uniform resistivity. The current leaving the center section of the electrode flows essentially along horizontal lines, delineating a surface having the shape of a biconcave lens. It is the elec- trical resistance offered to this particular group of current lines that is measured by the Guard tool. This is done in practice by having a long cylindrical electrode consisting of three sections physically and electrically separated by thin insulating washers, as shown in the center of the figure, and keeping the three sections at the same potential. It is obvious, for all practical purposes, that the current and potential distributions about this electrode are the same, as if we had a continuous electrode of same dimensions. Because the three sections are electrically separated, the resistance seen by the center section E can be independently measured. This section will be called the exploring electrbde, or E electrode. Its height generally is from a few inches to 1 ft. The outer sections are the guards G, which are several feet long. The combined length of the three sections is from 10 to 12 ft for the tools currently used in the petroleum industry. Let V denote the potential difference between E and the faraway return electrode, and I be the current flowing between these two electrodes. The resistance (Re) seen by the exploring electrode is VII, and the corresponding apparent resistivity is Ra = k(Re),............(1) where k is the tool calibration factor. The previous description establishes that the Guard tool uses that central portion of the sonde from which, in a uniform medium, the current follows essentially horizontal paths. When the formation being logged is fairly uniform, and when its resistivity is of the same order as

Jan 1, 1965

By R. C. Davis, R. F. Berry, G. W. White, B. T. O’Connell, L. A. Stacha

The formation density logging tool provides a radioactivity measurement that yields formation densities in situ. The relationship between bulk density and porosity is well understood. With knowledge of grain and fluid densities, porosity may be computed from the indicated formation bulk densities in shale-free formations. The porosities thus determined may be used with a resistivity log for water-saturation determinations. One technique utilizes a density vs resistivity plut. A plot of density vs Sonic-log transit time is used for porosity determinations in shaly sands. This method is particularly suited to wells drilled with oil-base or salt muds. In shale-free formations, comparisons of density values nnd neutron-log readings are used to identify lithologr and, thus, to select appropriate values of grain density for porosity computations. Through the use of Formation Density, Sonic and neutron logs, the interpretation problems caused by complex matrix lithologies are simplified. Other applications of the Formation Density log are found in the identification of minerals in evaporite deposits and in yield determinations of "oil shales". INTRODUCTION Widespread field application has demonstrated that the measurement of formation density is a useful and revealing technique for well logging. Responding to the bulk density variations of the formation, the Formation Density log affords valuable data on the porosity, lithology and fluid content of the formations. Although not without its own limitations, the Formation Density tool may provide effective interpretation under conditions that hamper or preclude the use of other porosity devices. When used in conjunction with other devices, such as the Sonic or neutron, complex matrix conditions may be better understood, defined and evaluated. A previous paper&apos; described the valuable contribution of the Formation Density tool in the logging of empty, or gas-filled, holes. The purpose of the present paper is to describe applications of the device in liquid-filled holes. THE FORMATION DENSITY TOOL The Formation Density tool was devised to measure the bulk densities of formations in situ. Actually, the device responds to electron density, but in the overwhelming majority of cases of interest, the electron density is very closely proportional to the bulk density. In other cases appropriate corrections can be made if necessary. Comparisons of log-derived densities with true densities for the more common matrix minerals are shown in Table 1. The tool consists of a gamma-ray source and detector mounted in a skid which, during the logging operation, is in contact with the borehole wall (see Ref. 1 for schematic). Gamma rays are emitted by the source and diffused through the formation. The number of diffused gamma rays reaching the detector at a fixed spacing from the source are counted. An increase in counting rate by the detector indicates a decrease in bulk density in front of the skid. Conversely, a decreased counting rate signifies an increase in bulk density. A measurement of hole diameter, between the face of the skid and that of the eccentering arm, is recorded simultaneously with the measurement of diffused gamma rays. Since the tool is a wall-contact device, the size of the borehole has little influence, provided good contact is maintained with the wall of the hole. Only in extremely large or small boreholes does the actual diameter become significant. However, for most accurate measurements the borehole size should be taken into account. Fig. 1 presents the departure curves for the conversion of counting rate to bulk density for liquid-filled holes of various sizes. It is of utmost importance that the skid maintain good contact with the wall of the hole. If, due to roughness or caving, there is mud between the skid and the formation, the tool will include in its response the relatively low-density mud, and the apparent density reading will be too low. Normally, permeable beds do not cave appreciably, so that reliable readings are usually available for zones of interest. Nevertheless, the caliper reading over the zone of interest should be examined for any wall roughness. When mud cake exists between the skid and the formation, the Formation Density log sees it as part of the formation. When thick mud cakes exist, a correction is required. For best results in wells having thick mud cakes,

By Edward B. Walker, Michel Grosmangin

The geological and economic conditions peculiar to the Greater Oficina area are presented to demonstrate the necessity of a low-cost, well-site method of distinguishing gas-bearing formations. The method of gas detection with the neutron log is briefly discussed and its limitations are explained. A new method is described which consists in running a second neutron log with a larger source-to-detector spacing. Field data indicate that under favorable conditions the larger spacing response is proportionately more influenced by the virgin gas zone than is the shorter spacing response. One important condition is that the gas sands be invaded, with a diameter of invasion equal to about two to three times the hole diameter. The calibration of the longer spacing log is made so that the two logs read the same deflection in shale and in a known water or oil sand. Under such conditions, gas is indicated by a positive separation, i.e., the curve for the longer device is displaced to the right of the curve for the shorter device. Because of the non-linearity of the long spacing response, low porosity formations also show positive separations. In sands with large porosity variations, the gas sands can no longer he located visually. They can be detected on a plot of shorter spacing vs longer spac-ing deflections. INTRODUCTION In the Greater Oficina area of Eastern Venezuela, the rapid and accurate detection of gas from well logs is of primary importance because of a combination of economic, geologic, and operational conditions. Briefly, these conditions are: economic—there is no market for gas at the present time; geologic—a single well may contain as many as 30 or more oil or gas sands for which intervals -to be tested must be selected; opera-tional—first, once opened to the casing through perforations, gas sands are extremely difficult and expensive to squeeze off, and second, gas detection must be rapid because the geologist must decide immediately, at the well, which sand to perforate. During the early period of exploitation in Greater Ofi-cina, when clay base muds were used, sidewall samples, together with electric logs, gave fairly good distinction between oil and gas sands. When exploitation became concentrated on deep fields of poorer sidewall sample recovery, and after the advent of oil emulsion muds, now the standard drilling fluid in Greater Oficina, the regular neutron log became the most reliable indicator of gas. In 1953, R. Norelius devised the dual-spacing neutron technique for gas detection. The dual-spacing neutron log has proved to be highly successful in detecting gas and has been adopted as a regular survey by all of the oil companies of the area. Several hundred dual-spacing neutron logs have been run in Eastern Venezuela and several tests per well are available for positive evaluation of the logs. The findings outlined below are based upon this extensive field experience. GEOLOGICAL CONDITIONS OF THE GREATER OFICINA AREA The Greater Oficina area lies on the south flank of the Eastern Venezuela Basin (Fig. 1). Nearly all accumulations of oil and gas in the area are trapped, in whole or in part, against normal faults. All known oil and gas reservoirs in Greater Oficina occur in the Ofi-cina formation (Miocene-Oligocene) or in the "U" sands. The U sands are the lithologic equivalent of the Merecure group (Otigocene-Eocene), which contains oil in many fields north of Greater Oficina. Most of the petroleum produced in Eastern Venezuela is from the Oficina formation.

Jan 1, 1958

By N. L. Muench, J. S. Osba

Research conducted previously has demonstrated that spectral analysis of gamma radiation from the naturally radioactive elements present in earth materials permits detervination of the relative concentrations of these elements. This research suggested that spectral analysis of gamma radiation resulting from neutron bombardment of elements not nor-mally radioactive might serve to identify these non-radioactive elements, and might permit their quantitative determination. To investigate this possibility, laboratory experiments were conducted in which various earth materials were subjected to bombardment by neutrons from a polonium-beryllium source, and in which the spectra of the resultant gamma radiation were analyzed. The geometry of the lab-ortory system was designed to simulate wellbore conditions. The laboratory work demonstrated that some elements not only can be identified with relative ease, but also can be determined quantitatively. This finding indicates that the use of an intense source of more energetic neutrons may lead to the identification and measurement of additional elements of interest by means of induced nuclear reactions. INTRODUCTION In recent years the applications of nuclear physics to the producing branch of the petroleum industry have become increasingly numerous. One significant development, an outgrowth of nuclear science which is of potential importance in well logging, has been the analysis of gamma-ray spectra of naturally radioactive elements in earth strata. The work of Brannon and Osoba1 has demonstrated that presently available equipment and techniques may be used to determine the relative concentration of thorium, potassium, and uranium in underground strata by means of gamma-ray spectral analysis. Success of the research on spectral analysis of naturally radioactive elements suggested that spectral analysis of induced radioactivity might serve to identify, and perhaps measure quantitatively, those elements which normally are not radioactive, but which become radioactive on neutron bombardment. It was to investigate this possibility that the research described in the following sections was undertaken. THEORY OF INDUCED NUCLEAR REACTIONS Of fundamental importance is the fact that an element, when bom- barded with neutrons, emits gamma rays of energies characteristic of that element. Therefore, in theory, an element may be identified by analysis of the energy spectrum of the neutron-induced gamma radiation. If a neutron impinges upon the nucleus of an atom, any one of several reactions may occur. For identification of elements in earth formations the more important are listed below: 1. The neutron may be elastically scattered, as in a collision between two billiard balls. In this case the total kinetic energy is conserved, and no nuclear radiation is emitted. 2. The neutron may be inelast-ically scattered, as in a collision between two automobiles in which case the total kinetic energy is not conserved. The lost kinetic energy reappears in the form of gamma radiation. 3. The neutron may be captured, in which case the neutron becomes a part of the target nucleus, and a gamma ray is emitted. For example, if the target nucleus is H1, the pro-duced nucleus will be H3, and gamma radiation of a particular energy will be emitted. 4. The neutron may produce nuclear fission, or may enter into miscellaneous reactions which result in the emission of

Jan 1, 1958

By J. L. Dumanoir, Maurice Martin

The induction-electrical log is a combination of a SP curve, a short normal curve, and an induction log; all three logs being recorded simultaneously in a single run. The spacing of the short normal is 16 in. The induction log is designated as 5FF40, which means that the sonde involves five coils and a focusing system, and that the distance between the main coils is 40 in. Examples of induction-electrical logs are shown in Figs. 13 and 15. The right-hand track contains the induction conductivity curve. The scale of this curve is linear with the units of conductivity, i.e., mhos/m or millimhos/m (mmhos/m), increasing from right to left. The same curve can be read in the reverse direction as a resistivity log, with a hyperbolic scale. In the middle track, the solid curve represents the short normal. (An amplified short normal is also shown in Fig. 15). The dotted curve is the resistivity curve of the induction log, i.e., the reciprocal of the conductivity curve, recorded with the same linear scale as the short normal. The SP is in the left-hand track. In this combination, the SP curve and the 16-in. normal are the essential correlation curves, as in the conventional electrical log. The induction log provides a more accurate and detailed record of formation boundaries than the conventional log, especially in hard rock territories. Furthermore, the induction-electrical log makes possible the determination of fluid content, particularly in thin beds—at least qualitatively and very often quantitatively. In such thin beds, because of the adverse effect of the adjacent formations, the responses of the conventional devices are generally doubtful and/or incom-plete. Another advantage of the new combination is that the induction log is little affected by the presence of thin hard streaks interbedded within the permeable sections, whereas such occurrences are the cause of serious distortions of the long normal and lateral curves. The evaluation of saturation from any resistivity log implies the knowledge of the true resistivity R, of the bed and of the resistivity R.. of the flushed zone, or at least of the ratio Rxo/Rc. With the induction-electrical log run in fresh mud, the 16-in. nor-mal has a smaller radius of investigation than the induction log, so that the former reads closer to R.. and the latter closer to Rt. Thus, it has been observed that the ratio of the 16-in. normal to the induction log readings can be used as a basis for the interpretation in terms of saturation, with the help of appropriate charts. The purpose of the present paper is essentially to explain this procedure and to provide the corresponding charts. On the other hand, when a formation is deeply invaded, the response of any log (in particular the induction log) depends upon the radial distribution of formation resistivities behind the wall of the hole, which itself is a function of the arrangement of the fluids saturating the formation pore spaces. This paper, accordingly, will begin with a schematic description of the distribution of fluids and resistivities in the invaded zone. Next, some essential characteristics of the induction log will be reviewed. Finally, the determination of the interpretation charts and their application will be explained. The discussion will be confined to the case of fresh muds, i.e., muds whose resistivity is at least five times as great as that of the formation waters. With the present equipment the induction log is not well adapted for logging in salty mud. The Laterolog and guard electrode log are at their best under conditions of low resistivity mud. Accordingly, when these tools were introduced a few years ago, the use of salty mud was recommended. Now, with the introduction of induction logging, spe-

Jan 1, 1958

By W. B. Nelligam, J. Tittman

Refinements in radiation logging techniques during recent years have involved increasing usage of scintillation ditectors. These detectors produce voltage pulses whose heights are related to the energies of the gamma rays which initiate them. Analysis of the gamma-ray spectrum, as indicated by the pulse heights, yields information about the chemical elements composing the formations surveyed. Refined scintillation counter techniques can furnish chemical information concerning earth formations in situ, from a study of the gamma-ray spectra emitted by the formation either naturally or as a result of neutron bombardment."&apos; Accompanying the rising interest in gamma-ray scintillation spectroscopy, there has been increased activity in the development of accelerator-type neutron sources (in contrast to encapsulated chemical-mixture sources). Such neutron generators are attractive for several reasons: (1) they greatly reduce radiation hazards to personnel; (2) there is a great reduction in contamination danger if they are lost in the hole; (3) they can produce larger neutron intensities than can conveniently available encapsulated sources; and (4) they are capable of being pulsed, thus permitting new techniques in logging. In the past, both accelerator and encapsulated neutron sources have been used by others in conjunction with scintillation-detector pulse-height analysis. The results have not been too encouraging, due to the interference among different gamma-ray spectral lines and to the fact that the gamma-ray peaks were not too clearly distinguishable above the large and ill-defined background "noise".&apos;.&apos; This paper is a status report on laboratory studies of a technique using a borehole accelerator as a neutron source, which gives an improved scintillation spectrum, thus permitting more accurate chemical analyses of the formations penetrated. The Schlumberger-accelerator neutron source&apos; is presented; the origins of inelastic and of thermal-neutron, capture gamma rays are dis- cussed, and results are given for some laboratory measurements performed in borehole geometry. THE NEUTRON GENERATOR-TUBE ACCELERATOR In the attempt to develop an accelerator neutron source which would have the desirable properties outlined in the "Introduction" and which would operate satisfactorily under borehole logging conditions, a small-diameter, cylindrical, neutron generator tube has been developed. This tube utilizes the principle of accelerating deuterons (heavy hydrogen nuclei) by a high voltage so that they bombard a tritium (heavy, heavy hydrogen) target. The resulting reactions produce large numbers of neutrons of 14 Mev energy. Furthermore, the tube is permanently sealed and, thus, the use of pumping techniques in the sonde is avoided. The tube consists of a pressure control, an ion source in which the deuterons are stripped of their electrons, an accelerating gap down which the deuterons are sped by the high voltage, a secondary electron suppressor and a tritium-loaded target. Obviously, the tube requires auxiliary circuitry in a sonde for the control of the ion source, tube pressure, high voltage and other operating variables. Neutron yields, both continuous and pulsed, have been produced under simulated field conditions in the range between 1 and 10 times those conventionally used in neutron logging applications. For the experiments discussed in this paper, neutron-pulse repetition rates in the range between 500 and 5,000 pulses/sec are adequate. Furthermore, by making suitable adjustments in operating conditions, one can vary the pulse width from a minimum of several microseconds up to dc. As will be seen later, the present application to well logging does not require that one have available pulses of neutrons which are appreciably shorter than the time it takes for fast neutrons to slow down to thermal energy in water. As a consequence, much of our interest has been directed towards the operating characteristics of the tube with pulses longer than about 10 microseconds in duration. Our experience to date shows that tubes can be operated in this manner with reasonably constant, average neutron outputs.

Jan 1, 1966

By R. H. Winn

In this paper a heterogeneous carbonate reservoir is considered as a succession of layers of formations of radically different porosities, permeabilities, water saturations, and rock matrix types; i.e., the heterogeneity is vertical rather than radial to the wellbore. Presumably each layer is more or less homogeneous. To determine the "net oil in place" in this reservoir it is desirable to know the porosity and water saturation of each layer. In addition, the section must be subdivided into that which will produce oil, that which will produce water, and that which will produce little or nothing. In numerous wells drilled in these reservoirs radioactivity surveys, guard electrode surveys, and caliper surveys have been made. These logs have been compared with core analysis data, when available, and the analysis of the surveys compared with production data. This study, still in progress, is intended to provide improved analytical techniques if such are possible. This paper is a progress report on the work. DETERMINATION OF POROSITY The neutron-gamma log has long been recognized as a measure of hydrogen concentration.1,2,3 In reservoirs that are free of clay, shale, or other hydrogen containing matrices, it has been accepted as a measure of porosity.4,5 Empirical methods to relate relative neu-tron-gamma intensity to porosity are often employed." To provide a more precise determination of porosity from the neutron-gamma log each instrument is calibrated immediately prior to each survey by a method which refers all tools to a hydrogenous cylinder. To insure stability a calibration check is made after each survey. These calibration records are printed as an integral part of the log that is recorded in standard counts per second. A detailed description of the calibration procedure is available in the literature.&apos; The procedure employed insures that each Survey is numerically reproducible within the limits of statistical variation. To relate the standard neutron counts per second to porosity and hole size, records were made in test wells constructed of blocks of limestone. The porosities of the blocks were accurately measured and the blocks completely saturated with fresh water. Fig. 1 presents some of the results obtained in these test wells.&apos; The lines shown in Fig. 1 were obtained by deriving an empirical equation relating the 11 points shown so that the relationship may be projected to larger and smaller hole sizes. A detailed description of the tool is available in the literature.&apos; Field tests indicate that the tool is adjacent to the wall of the hole most of the time when the surveys are made in uncased holes drilled with saline muds. In this circumstance porosity may be calculated from Fig. 1 if the hole size is known. When fresh muds that form thick mud cakes are employed, the centralizing effect of the mud cake must be considered and the porosity determination is less reliable. In cased holes the combined effect of variations in hole diameter and the unknown eccentricity of the casing in the hole may make porosity determinations extremely unreliable. In any event, an accurate knowledge of hole size is essen-

Jan 1, 1958

By P. E. Baker

The commercial neutron-gamma log measures the intensity of gamma-rays produced by the interaction of neutrons with nuclei in the formation and in the bore hole. This measurement gives information concerning the distribution of neutrons around the neutron source, and in particular gives information on porosity or fluid content. Another type of information is contained in the frequency or energy spectrum of the neutron capture gamma-rays. The capture gamma-rays of each element, in fact of each isotope, have only those energies characteristic of that isotope.&apos;" The neutron capture gamma-ray spectrum, if obtainable, would permit the identification of the elements producing the spectrum. Such information would have considerable value in oilfield exploration work if it could be applied practically in well logging. Experiments have been carried out at the La Habra Laboratory of California Research Cop. to determine the feasibility of a method of well logging in which neutron capture gamma-ray spectra are measured down hole. In these experiments, spectra were obtained for artificial formations under simulated borehole conditions, using actual borehole geometry. This paper will describe the experimental apparatus and methods and will present the spectrum curves for several formations. THEORY Basic Formulas The neutron capture cross section of an element is a measure of its tendency to capture neutrons. The macroscopic cross section, or cross section per cubic centimeter, of a material composed of several elements is the sum of the cross sections of the nuclei making up 1 cc of the material. This is expressed in the equation SB = Sn1s2, ..........(1) where SB is the macroscopic cross section, s, is the cross section of a nucleus of type i, and ni is the number of nuclei of type i/cu cm. If N is the number of neutrons absorbed per second by the composite material, then the number captured by element j is Ni = N nisi/SB = N nisi/Si ni si ......(2) When the nucleus captures a neutron, a new nucleus is formed, which is usually in an excited state and falls to its ground or unexcited state by emitting one or more gamma-ray quanta.* The emission spectrum has been measured for many elements (see Table 1). The intensity of gamma-rays produced by an element is proportional to the number of neutrons captured per unit time by the element, which in turn depends on the capture cross section of that element and the capture cross sections of the other elements present which are competing with the first element for thermal neutrons. In a small sample of material the number of neutrons, N, captured per second is N = ?SB, ..........(3) where $ is the thermal neutron flux and SB is the macro-scopic cross section defined above. The number of neutrons captured by element j is N1 = ?ni si = N ni si/SB . ......(4) The number of gamma-ray quanta of energy, E, pro-duced per unit time by element, j, is I(E)=Nni si/SB P(E) .......(5) where P(E) is the emission probability of element, j, for this gamma-ray energy. Borehole Effects The capture gamma-ray spectra of many elements have been measured accurately in academic and government laboratories. These measurements are made on small samples, and from these measurements the emission spectra of the elements are obtained.

Jan 1, 1958