Search Documents

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' 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 G. H. F. Gardner

The boundary conditions F = 1.0 and 4 = 1.0 (100 per cent) are dictated by the definition of F in the equation 1/F = aFp. These conditions have physical significance as the "limiting conditions" for the observed behavior in systems of dispersed particles. The other boundary conditions, cited by Gardner in his Eq. 3 as applying to this system, have no significance in the systems we studied. Maxwell's equation (Eq. 2) cited in Gardner's comments was developed for and limited to a regular, rather than random, dispersion of spheres which are sufficiently removed from each other so that they do not interfere with current flow about their neighbors. These restrictions also apply to Rayleigh's analysis. The title of Rayleigh's paper cited by Gardner is "On the Influence of Obstacles Arrayed in Rectangular Order upon the Properties of a Medium". He is not discussing randomly dispersed spheres. It might be inferred from Gardner's comments that Rayleigh specifically considered the differential equation (Gardner's Eq. 3). This is not the case. Mathematically, Gardner's Eq. 3 is exactly true for the system to which it applies, however, as both Rayleigh and Maxwell note that their equation for F in terms of F is not for randomly dispersed spheres. At the time we wrote our paper, we were aware of the work of Wyllie and Gregory as cited by Gardner. We cited their work in Ref. 5 of our paper. Our comment on page 287 of the paper about the validity of an m value of 1.3 for spheres with a porosity of 38 per cent remains as our answer to Gardner's question on this point. Our data, at 38 per cent porosity, show that the condition a = 1 and m = 1.3 holds for spheres near 40 per cent porosity.

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 W. A. Hoyer

A pulsed high-energy neturon-induced spectral logging tool has been built and field-tested. The reaction of deuterium on tritium is used to generate pulses of 14-Mev neutrons. By detecting only the prompt gamma rays produced by neutron inelastic reactions in the formation, the presence and relative abundance of carbon, oxygen, calcium, silicon and other important elements may be ascertained from a gamma-ray spectrum. Gamma-ray spectra obtained in a shallow test well and in experimental field use show that it is possible to identify formations and their contained fluids. INTRODUCTION The penetrating gamma rays from naturally occurring radioactive elements in subsurface formations have been used for a number of years in well logging as a means of characterizing and distinguishing strata. Still another nuclear logging method which has been employed for some years consists of the bombardment of strata with neutrons and the measurement of the number of gamma rays produced by neutron capture reactions involving elements in subsurface strata.' Since one of the principal elements entering into capture reactions is hydrogen, the latter procedure essentially results in a hydrogen log, whether the hydrogen be in combined form in either oil or water; thus, the log gives an indirect measure of porosity. These methods, with various refinements, have been developed to such an extent that they have become routine procedures in formation evaluation. Neither gamma-ray logging nor conventional neutron logging, however, yields sufficient information to permit unequivocal identification of the mineralogic composition of formations, and neither method gives information which may be used for the positive identification of hydrocarbons in strata. Accordingly, a number of efforts have been made in recent years to gain additional information from nuclear logs. Brannon and Osoba have shown that it is possible to identify naturally occurring radioactive elements in subsurface formations by spectral analysis of gamma rays emanating from these elements. Such an identification is of value in the characterization of strata. A simplified form of spectral analysis of gamma radiation resulting from neutron capture reactions between elements in earth materials and bombarding neutrons has been used with some suc- cess under favorable conditions to differentiate between petroleum and water. This method relies upon the relatively high energy of gamma radiation from neutron capture by chlorine and, in effect, furnishes a chlorine log. In areas in which interstitial water is of sufficiently high salinity, this log can give valuable information on water saturation and, thus, indirectly on hydrocarbon saturation. Still another approach to obtaining more information by nuclear techniques is "activation logging", in which certain elements yield short-lived radioactive isotopes on neutron activation and, thus, can be detected by gamma-ray spectral analysis. From the standpoint of determining by nuclear logging methods the mineralogic composition of strata and the presence or absence of hydrocarbons, it is essential that information be obtained on the presence and relative amounts of several elements such as carbon, oxygen, hydrogen, calcium, silica and others. Some elements, such as hydrogen, chlorine and sulfur, can be determined by spectral analysis of gamma radiation resulting from neutron capture reactions."" Others — carbon and oxygen, for example — do not enter readily into capture reactions but do yield gamma rays of characteristic energies from inelastic scattering reactions with high-energy neutrons. Accordingly, the latter reactions are of particular interest as a means of identifying hydrocarbons in subsurface strata. Several years ago it was recognized that information requisite to the identification of subsurface strata and of contained fluids could be obtained by pulsed operation of a subsurface neutron generator and associated gamma-ray detector. It was contemplated that pulsed operation would be effective in discriminating in time that gamma radiation which results predominantly from inelastic scattering reactions from that gamma radiation arising from neutron capture reactions. Accordingly, research was continued and construction of such a device was initiated. Although many problems were encountered, they were solved successfully. The following sections describe the tool which has evolved, its performance and the results which have been obtained with it to date. THEORY OF OPERATION OF NEUTRON-INDUCED GAMMA-RAY SPECTRAL LOGGING TOOL THEORY OF INDUCED NUCLEAR REACTIONS Of fundamental importance is the fact that an element, when bombarded with neutrons, emits gamma radiation of energy characteristic of that element. In theory, therefore, an element may be identified by

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."' 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".'.' 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' 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.' 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.' 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.' 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.'" 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

By E. C. Hopkinson, A. H. Youmans, R. A. Bergan, H. I. Oshry

A new log has been developed for quantitative formation evaluation which is based on a measurement of the length of time slow neutrons survive before they are captured in the rocks and fluids. The logging instrument employs a cyclically pulsed neutron generator and a gated scintillation counter which is synchronized with the source. The source emits short, intense bursts of 14 mev neutrons once every 1,000 microsec and is quiescent between bursts. During the period the source is quiescent, the detector is electronically actuated for two independent preselected intervals. A comparison of the counting rates during these two intervals gives a measure of the rate of decay of the slow neutrons and of the associated gamma radiation. The average neutron lifetime in most earth formations is in the range from 50 to 500 microsec. It can be measured during a continuous logging operation at conventional logging speeds. The design of the logging instrument is described and the results of tests are compared with theoretical predictiom. Formulas are developed which give the relationship between log response and formation properties. It is shown that the method is particularly sensitive to formation fluid salinity, and that salt water saturation can be measured accurately in either cased or open hole. The measurement can be made independent of borehole size, fluid type, casing and tool position in the hole by properly selecting the intervals during which the measurements are made. The results of tests with a prototype logging tool are given. INTRODUCTION A new nuclear logging system has been developed which employs the Accelatron,* an accelerator-type neutron source, and accurately measures formation brine saturation in an entirely new way. It has produced a type of formation log with better sensitivity, greater sampling depth and simpler quantitative interpretation than any other nuclear log thus far suggested. The new Neutron Lifetime Log* employs a pulsed electromechanical neutron source and a synchronously gated radiation detector. A prototype instrument has been field tested during recent months to demonstrate the operability of the apparatus and the feasibility of the method. Tests in wells and simulated boreholes have confirmed theoretical predictions and have shown that formation param ters can be measured independent of casing and other borehole parameters. Preliminary results of field tests have indicated that the log may have important and widespread applications. BASIC PRINCIPLE OF NEUTRON LIFETIME LOG The Neutron Lifetime Log is based on the fact that neutrons emitted by a source in a well are rapidly but not instantly captured by the material around the source. Their capture is a matter of statistical probability; the greater the number of capturing nuclei and the greater the "capture cross section", the greater is the probability that a neutron will be captured quickly. The average life of a thermal neutron in vacuum is about 13 minutes, but in common earth materials, the average neutron life ranges between extremes of about 5 rnicrosec for rock salt and perhaps 900 microsec for quartzite. The Neutron Lifetime Log responds to variations in this average neutron life. The theoretical basis for a log of this general type has been well understood by nuclear logging experts in many laboratories both in America and in Russia, and develop mental work along these lines has been in progress for many years. The Russian literature has reported both theoretical and experimental work1,2 but in this country there have been no published reports of progress toward a practical logging instrument. The logging instrument is designed to measure radiation produced by slow neutrons during selected intervals when no neutrons are being emitted by the source. The source is arranged to emit neutrons in bursts or pulses. During the quiescent interval between the pulses, it is possible to observe the exponential "decay" of the neutrons and the neutron-induced radiation as the individual neutrons progressively disappear due to capture by atoms in the formation or the borehole. When a short pulse of 14 mev neutrons is emitted by a source in a borehole, the individual neutrons are slowed to thermal energy within a few microsec. Thus, a cloud of "slow" neutrons is formed around the source within 10 to 50 microsec after the pulse. This cloud is most dense within a few inches of the source, and is progressively less dense out to a radius of about 3 ft, where radiation from the source is practically undetectable.

Jan 1, 1965

By J. T. Dewan

An experimental setup is described whereby the responses of Neutron logging instruments have been determined opposite formations of different porosities and under various borehole conditions. The results of these tests are presented in several charts which show how the Neutron log is affected by the following factors: hole diameter, casing, cement, salinity and weight of drilling mud, and position of logging tool with respect to the wall of the hole. The responses in the case of empty holes are also shown. A chart is furthermore presented which shows the shapes of the Neutron log, as a function of logging speed and time constant, at the boundaries between two formations with different porosities. These were derived mathematically using as a basis laboratory determination of the shape of the curve under static conditions, i.e. at very low logging speed. Corrections for the readings made opposite thin beds are given by these computations. The correction is found negligible under usual logging conditions for beds thicker than about 4 ft. INTRODUCTION The Neutron log is used in borehole surveys both for the delineation of porous formations and for the quantitative estimation of their porosities.1,2,3,4,5 It is common practice to calibrate the response of the Neutron log in terms of porosity values for each given geological unit by means of empirical comparisons between the deflections of the neutron curve and the porosity values given by the core analysis in a sufficient number of wells. Usually, the plots of the porosity values vs Neutron log deflections show a definite trend. It is therefore possible to trace an average line across the plotted points. This average line constitutes a calibration curve which can be used to evaluate the porosities from the Neutron logs in the other wells drilled in the same region. Futhermore, when the porosity values are plotted on a logarithmic scale and the Neutron log deflections on a linear scale, the calibration curve is practically a straight line. Accordingly the calibration curves are often determined by simply tracing a straight line between two points corresponding to two widely different porosities. Different choices are made for these two points, depending on local conditions. For example, one point will be defined by plotting the deflection of the Neutron curve opposite a formation known to be very tight, for which an average porosity value of 1 or 2 per cent is assumed. The other point may correspond to shale sections, to which a porosity of around 40 per cent is ascribed. The information provided by the eletrical logs — in particular the Laterolog, Limestone log or Micro-Latero-log — is also very often a help in the calibration of the Neutron logs. Empirical approaches such as these make possible the derivation of fairly reliable porosity values from the Neutron logs in many cases, at least in open holes and in essentially clean formations. There are, however, several causes of uncertainty which may sometimes entail appreciable errors. On the one hand, the Neutron log readings are influenced by the borehole conditions: hole diameter, presence or absence of casing and cement, presence or absence of drilling mud, characteristics of the mud, and position of the logging instrument with respect to the wall of the hole. The amplitudes of the curves opposite thin beds are furthermore atfected by the bed thickness,

Jan 1, 1957

By R. P. Alger, M. P. Tixier, D. R. Tanguy

In the combination induction-electrical log used at present in the field, the induction logging tool is appropriate for the investigation of moderately invaded formations. A new induction sonde with a radius and investigation about twice as large has been developed recently for the case of deep invasion. It has very nearly the same vertical resolution as the present sonde so that thin beds are defined as accurately as before. The characteristics of the new tool are described, the corresponding interpretation charts are given and field examples are discussed. The design of the sonic logging tool has been modified to improve the calibration and the reliability. The fact that porosity can be accurately recorded by means of the sonic log has prompted new interpretation procedures for saturation estimation, wherein the data concerning the various permeable beds in a given well are correlated. One approach consists of plotting transit time vs true resistivity, with an appropriate scale. With this approach, saturations can be estimated conveniently even in cases where formation water resistivity is not well known. In another approach, a comparison is made of the values of the formation waters computed from the re- sistivity and sonic logs. Using the concept of continuity, this procedure makes possible a quick determination of zones of saturation in shaly sands and/or in case of appreciable variations of formation salinities with depth. It has been found that the comparison of porosity from the sonic log with the apparent porosity computed from a short-investigation resistivity log may reveal, in many cases, the presence of residual oil and thus detect potentially productive formations; this procedure is valuable when the true formation resistivity and the resistivity of the formation water are in doubt. INTRODUCTION During the past year, the efficiency of log interpretation has been vastly improved. The improvements have largely resulted from the introduction of a deep-investigation induction device and from the application of new interpretation techniques that utilize sonic vs resistivity readings. Since the new interpretation techniques depend, in part, upon good values of true formation resistivity, the new induction log will be discussed under Part I. The sonic interpretation techniques will be studied under Part 11. Early in 1959, the 6FF40 induction equipment was introduced in the field. This device was designed for a better approach to true formation resistivities in deeply invaded zones. The greatly improved radial investigation of the 6FF40 equipment has been achieved without sacrificing vertical resolution. The first combination induction-electrical log, the 5FF40, was introduced as a standard tool in 1956 for the logging of wells drilled with fresh muds. The tool has received wide industry acceptance in the United States. 'The 5FF40 induction log has a radial investigation sufficient to overcome average depths of mud filtrate invasion. At 5d invasion, for example, the 5FF40 induction log will read about 1.4 R. in a water sand where R., = 10R,. At 10d invasion, such induction log would read 2.45 R. in the same water sand. In either case, the effects of invasion would not be sufficiently great to cause a water sand to be mistaken for a shale-free oil- or gas-producing zone. Some formations, however, invade deeply — in excess of 10d. Such water zones could be mistaken for oil-or gas-saturated sands unless the porosity balance' can clearly make the distinction. It is for these deeply invaded formations that the 6FF40 was developed. CHARACTERISTICS OF THE 6FF40 Radial Investigation Characteristics To describe the comparative responses of the 5FF40

By R. J. S. Brown, B. W. Ganison

A new logging method has been developed, based on measurement of the nuclear magnetism of formation fluids. The nuclear magnetism log (NML) is the only log that responds solely to formation fluids. It operates equally well in both oil-base and water-base muds and in empty holes, and can be used in all kinds of formations except strongly magnetic ones. Two separate NML measurements can be made, one of which provides a continuous formation fluid curve. This fluid curve is called the free fluid log (FFL) and is believed to indicate a minimum effective porosity in most formations. The FFL not only delineates fluid-containing zones, but provides an excellent correlation curve that can be obtained under conditions where conventional correlation logs are ineffective. Preliminary tests indicate that the second kind of NML measurement may help distinguish oil and water zones and provide information concerning permeability and wettability. (The FFL itself appears to provide some information on permeability.) The second kind of NML measurement requires stopping the logging tool for a short time opposite a zone of interest and taking more extensive NML data that can be displayed as nuclear magnetic relaxation curves. In some instances, oil and water saturations for the region immediately adjacent to the borehole can be read from these relaxation curves. INTRODUCTION In 1946, Bloch, Hansen and Packard and Purcell, Torrey and Pound3 independently announced the successful demonstration of the phenomenon of nuclear magnetic resonance. During the past 13 years, there have been many applications of nuclear magnetic resonance, including applications to the study of chemical structure and to the measurement of magnetic field strengths. Preliminary experiments on the feasibility of using nuclear magnetism measurements in well logging were made independently by California Research Corp. and Varian Assoc., the Varian work being sponsored by the Byron Jackson Tools, Inc. Since then a cooperative research program on nuclear magnetism logging has been carried out by the Byron Jackson Div. and Research Center of Borg-Warner Corp., and California Research Corp., subsidiary of Standard Oil Co. of California. The use of nuclear magnetism in well logging is of special interest because it offers a way of making direct measurements on the hydrogen in the formation fluids and not on the rock matrix. Within the past 1 1/2 years, successful measurements have been made with a research model logging tool in wells in California, Tex as, Utah, Louisiana and Wyoming. NUCLEAR MAGNETISM SIGNALS Polarization, Relaxation and Precession Many atomic nuclei possess magnetic moments and spins; that is, they are similar in some respects to bar-magnet and gyroscope combinations. Molecules and their nuclei are subject to thermal motion, which has a scrambling effect, tending to leave as many nuclear spins oriented in any one direction as in any other. However, if a magnetic field is applied, the magnetic nuclei tend to align in the direction of the field. The scrambling and aligning forces compete with each other, with the result that a few more spins are oriented parallel to the field than in other directions. This gives a net magnetization, or polarization, which is directly proportional to the strength of the applied magnetic field (aligning influence) and inversely proportional to the absolute temperature (scrambling influence). When the magnetic field in, or temperature of, a liquid sample containing protons is changed, the new equilibrium value of proton polarization is not established immediately but requires an amount of time which depends on the nature of the hydrogen-containing materials. The process of approaching the equilibrium value of polarization is called relaxation.'.' Polarization is a vector quantity, and the components parallel to and perpendicular to the magnetic field must be considered separately. Relaxation of the component parallel to the field is called "thermal relaxation", or "longitudinal relaxation", and the corresponding time for this component of non-equilibrium polarization to decay by a factor of e (natural log base) is denoted T. The relaxation of the perpendicular component is called "transverse relaxation", and the corresponding relaxation time is denoted T,. The potential energy of a magnet in a uniform field depends on the angle the magnet makes with the field; therefore, a change of the component of net polarization parallel to the magnetic field involves an exchange of energy between the spin system and the thermal motion of the molecules, leading to the term thermal relaxation for the relaxation of this component. Suppose we subject a sample to a strong magnetic field at right angles to the earth's field for a time greater than TI. A polarization, thus, is established at right an-

By J. E. 214-000-000-010 Warren, M. R. J. Wyllie, T. Meidav, L. Scharon, R. Uhley, A. J. deWitte

By considering the stoichiometry of the underground combustion process, an equrztion was derived relating the point velocity of the combustion front as a function of the air flux, fuel content, efficiency' of oxygen uitilization, hydrogen:carbon ratio, and ratio of carbon dioxide In carbon monoxide produced. In addition, an equation for cnlcrllating the amount of injected air necessary 10 produce the oil from I Cu. ft of reservoir space was derived. This equation has a direct bearing on evaluating compression costs for the process. Values calculated from these equations are in good agreement with experimentnlly obscrvcd values obtained from laborator) tube runs reported in the literature. The equations are also u\eful for atlcrlyzinf: experinlental data for inconsistencie.r.. I NTRODUCTION The injection of energy in the form of heat into an oil-bearing formation as a means of increasing the recovery of oil has been the subject of much discussion and experimentation as far back as 1920. The ideas have ranged from the injection of steam to the direct use of electricity as a means of getting thermal energy into a formation. However, most laboratory and field studies have been confined to the process whereby a portion of the crude in the reservoir is burned to generate the necessary heat to increase the mobility of the oil by reducing its viscosity. A recent paper by W. L. Martin, et al.', described the results of a number of laboratory experiments in which the combustion process was studied in a sandpack. Studies such as this indicate that with the temperatures involved, the oil is partially coked in this process and it is this coke which provides the fuel for combustion. The combustion of this residual coke can be analyzed theoretically by applying simple stoich-iometric principles to the combustion process. It is the purpose of this paper to compare the results reported in these experiments with those obtained from theoretical considerations. The variables studied by Martin, et al, among others, include combustion front velocity, air flux, oxygen utilization, fuel concentration, ratio of carbon dioxide and carbon monoxide in combustion gases, and the hydrogen:carbon ratio of the residual fuel. All these variables come into consideration when the stoichiometry of the combustion process is studied. The experimental conditions and results are tabulated in Tables 3 and 4 of Ref. 1. For purposes of comparison with calculated results, it was necessary to change the units of a number of the variables. UNDERGROUND COMBUSTION STOICHIOMETRY The combustion or oxidation of the carbonaceous residue or "coke" on the sand can be described in terms of the equation: where n is atomic ratio of hydrogen to carbon, m is ratio of moles of CO2 to CO produced, and CH. is taken as the unit molecular weight, realizing that actually the coke has a large range of high molecular weights. Assuming that all the fuel laid down in the sand is burned, the ratio of moles of oxygen consumed to moles of coke can be written as follows. Moles 0MolesMolescoke In terms of moles O2 per pound of coke, where 12 + n is the molecular weight in pounds per mole. If the 0, consumed is measured in standard cubic feet (60°F and 1 atm) and the coke in pounds per cubic foot of formation, 2, then:

By M. P. Tixier

In several log interpretation methods, water saturation is evaluated by use of the ratio of the readings of a short spacing resistivity device and a long spacing resistivity device plus information on the mud filtrate and connate water resistivities, which is often derivable from the SP. These methods are valid over a certain range of conditions, usually specified by moderate invasion and R.. greater than R,. Since these methods involve no explicit evaluation of the formation factor, F, the saturafiion so found may be used in the standard Archie equation to derive a computed formation factor which, as a check, may be compared with the formation factor known from other independent measurements. Discrepancies in the two values of formation factor generally indicate that the ratio method is not within its range of applicability. If the computed F is too low, the corresponding S, is too low; if the computed F is too great, the corresponding S, is too great. By means of this "porosity balance" check, and some knowledge of the probable conditions of invasion, the interpretation can often be improved. The poqosity balance check is discussed for the cases of the Induction-Electrical Log Interpretation Method, the Rocky Mountain Method, and the R.,/R, Method. For the first method, a discussion is also given for the case of shaly sands. INTRODUCTION Al though many logging tools are available to obtain a practical value of porosity, the only logging method reflecting the water saturation is electrical logging. Water saturation is commonly computed from electrical log data by means of the basic Archie formula, Where F is the formation factor, R, the resistivity of the formation water, and R, the true resistivity of the formation. One difficulty in this procedure is a possible error in the evaluation of R, from the logs— for example, in the case of very deep invasion. There is no possibility, inherent in the procedure, which would provide a check of the accuracy of R, and, hence, of the computed saturation. Other approaches have been proposed which are based essentially on the ratio of the readings made with two dzerent devices: one device, with a short radius of investigation, provides a value close to the resistivity of the invaded zone; the second device, with a long radius of investigation, gives a value close to the true resistivity of the bed. By these methods, the determination of saturation is made independently of any direct knowledge of the formation factor. Therefore, if the formation factor (or porosity) is known from other sources (core analysis, neutron log, MicroLog . . .), it becomes possible to take advantage of this additional information in order to check the value of saturation derived from these latter methods. Inasmuch as this proposed check procedure relies on a knowledge of the value of formation factor or porosity, it has been deemed appropriate to call this verification a porosity balance. The purpose of this paper is to discuss the principle and the application of the porosity balance, depending on the method of interpretation being used. Basis of the Method There are in existence three interpretation methods by means of which the water saturation can bc found without a knowledge of the formation factor (or porosity). As just stated, they are based on the fact that the ratio of the two resistivities (one measuring the resistivity close to the hole, the other measuring the resistivity far from the hole) is independent of the formation factor since it affects the two resistivity measurements in the same way. Chronologically, the resistivity ratio methods are: (1) the Rocky Mountain Method, which has been used since 1943 and was explained in a paper in 1949;' (2) the R../R, Method or Shaly Sand Method, presented in 1954, which permits the calculation of the water saturation in both clean and shaly sands2; and (3) the Induction-Electrical Log Method, presented in 1957 for the interpretation of the induction-electrical log combination." Each method requires certain definite conditions in order for the interpretation to be reliable. The porosity balance method will let us know whether these conditions are fulfilled. In clean formations we know that the fundamental Archie relation (Eq. 1) must always be verified. This relation can be rewritten:

By A. A. Venghiattis

A rational approach to the selection of the appropriate perforator to use in each specific zone of an oil well is presented. The criteria presently in use for this choice bear little resemblance with actual down-hole condilions. These environmental conditions affect the elastic properties of rocks. One of these elastic properties, acoustic velocity, is suggested as the leading parameter to adopt for the choice of a perforator because, being currently measured in the natural location of the formation, it takes into account all of the effects of compaction, saturation, temperature, etc., which are overlooked in the laboratory. Equations and curves in relation with this suggestion are given to allow the prediction of the depth of perforation of bullets and shaped charges when an acoustic log has been run in the zone to be perforated. INTRODUCTION When an oil company has to decide on the perforator to choose for a completion job, I wonder if it is really understood that, to date, there is no rational way of selecting the right perforator on the basis of what it will do down-hole. This situation stems from the fact that the many varieties of existing perforators, bullets or shaped charges, are promoted on the basis of their performance in the laboratory, but very little is said on how this performance will be affected by subsurface conditions such as the combination of high overburden pressure and high temperature, for example. The purpose of this paper is to show the limitations of the existing ways of evaluating the performance of perforators, to show that performances obtained in laboratories cannot be extended to down-hole conditions because the elastic properties of rocks are affected by these conditions and, finally, to suggest and justify the use of the acoustic velocity of rocks, as the parameter to utilize for the anticipation of the performance of a perforator in true down-hole environment. EVALUATING THE PERFORMANCE OF A PERFORATOR It is natural, of course, to judge the performance of a perforator from the size of the hole it makes in a predetermined target. Considering that the ultimate target for an oilwell perforator is the oil-bearing formation preceded in most cases by a layer of cement and by the wall of a steel casing, the difficulties begin with the choice of an adequate experimental target material. For obvious reasons of convenience, the first choice that came to the mind of perforator designers was mild steel. This is a reasonable choice for the comparison of two perforators in first approximation. Mild steel is commercially available in a rather consistent state and quality, and is comparatively inexpensive. The trouble with mild steel is that it represents a yardstick very much contracted; minute variations in depth of penetration or hole diameter and shape may be significant though difficult to measure. The penetration of projectiles in steel being a function of the Brinell hardness of the steel (Gabeaud, O'Neill, Grun-wood, Poboril, et al), it is often difficult to decide whether to attribute a small difference in penetration to a variation on the target hardness or to an actual variation on the efficiency of the projectile. Another target material which has been widely used for testing the efficiency of bullets or shaped charges in an effort to represent a formation—a mineral target as opposed to an all-steel target—is cement cast in steel containers. This type of target, although offering a larger scale for measuring penetrations, proved so unreliable because of its poor repeatability that it had to be abandoned by most designers. The drawbacks of these target materials, and particularly their complete lack of similarity with an oil-bearing formation, became so evident that a more realistic target arrangement was sought until a tacit agreement was reached between customers and designers of oilwell perforators on a testing target of the type shown on Fig. 1. This became almost a necessity about seven years ago because of the introduction of a new parameter in the evaluation of the efficiency of a perforator, the well flow index (WFI). The WFI is the ratio (under predetermined and constant conditions of ambiance, pressure and temperature) of the permeability to a ceitain grade of kerosene of the target core (usually Berea sandstone) after verforation. to its vermeabilitv before perforation. The value of this index ;or the present state if the perforation technique varies from 0 to 2.5, the good perforators presently available rating somewhere around 2.0 and the poor ones around 0.8, There is no doubt that, to date, the WFI type of test is by far the most significant one for comparing perforators. It is obvious that a demonstration of a perforator

By J. P. Myers, A. H. Flagg, E. S. Mardock, J. L. P. Campbell, J. M. Terry

The development process for the use of radioactive tracers as a means of locating zones of permeability is discussed. The general techniques for the safe handling of radioactive materials is given as developed by the Lane-Wells Co. and Well Surveys, Inc. The prob1em and. successes with tracers in waterflood systetns, oil and gas injection profiles, fracture sand tracing, squeeze (ement tracing. lost circulution, and cement top localion are discussed and illustrated. INTRODUCTION The exact location of the permeable zones lying within the productive horizon has been a major problem faced by the petroleum engineer for decades, and it has been a problem which has seldom been solved with certainty. Reservoir engineers usually have to incorporate a question mark in their estimates of future reservoir performance because of blanks in the core data. Although the use of radioactivity or electrical logs for the estimation of porosity has been of great assistance in furnishing a record for the reservoir engineer, this application of logging data has not completely solved the problem. In short, the industry needs an in-situ well surveying method which can locate and estimate the permeability of the zones existing within a given pay section. The initial development work with radioactive tracers has been towards filling the need for permeability information. Progress to date permits the location of permeable zones with a qualitative estimate of their relative importance with respect to permeability. The calibration of the method is still incomplete. but present development work on tracer preparations, borehole effect correction factors, and instrument limitations give a promise of quantitative "permeability profiles" in the near future. The development work on the radioactive tracer method discussed here began several years ago in the laboratories of Well Surveys, Inc., in Tulsa. At the outset of this work, three approaches to the problem were visualized as follows: 1. A fluid carrying a soluble, chemically inert, gamma ray emitting isotope is pumped down the borehole and out into the permeable zones of the formation. The greater gamma activity of the permeable zones thus permits them to be differentiated from the impermeable zones which contain no active tracer. It was found difficult to obtain a tracer material which was chemically inert. Absorption to some degree was noted in all experiments. Also, it proved impossible to hold the tracer near the horehole where it could be logged. 2. A fluid carrying a soluble, gamma ray emitting isotope is pumped down the borehole into the permeable formations. In contrast with proposal 1, above. the tracer is chemically active with the fluids and materials in the pore spaces of the permeable zones. The permeable zones are thus indicated by their higher radioactivity level. The main difficulty in this proposal was finding a tracer material which would not be absorbed by shale as well as in the permeable porosity. Very little field work was done with this system. 3. Radioactive particles suspended within the injection fluid are pumped down the well. As the fluid enters the permeable formations, the radioactive particles are filtered out on borehole wall. The permeable zones are differentiated from the impermeable by the gamma ray activity deposited on the formation wall.

Jan 1, 1956

By J. W. Graham

For many years the author has been cognizant of the difficulty encountered by some in treating with the water influx formulas for unsteady-state fluid flow as pertain to the material balance equation. This has particularly applied in establishing reservoir performance and identifying reservoir pressure, which to the practicing engineer has entailed a trial-and-error procedure, and for others has necessitated resorting to computing devices and reiteration processes. In retrospect this difficulty stems from the fact that reservoir pressure in the material balance formulas, as well as associated with the water influx equations, is an inexplicit term, and the work reported in the past is irrefutable. However, what will be presented in this paper is another approach to the problem, whereby the entire material balance equation will be treated by the Laplace transformation, and reservoir pressure which hereto has been inexplicit, can now be isolated by mathematical procedure to relate that parameter with all the factors contributing to its change. This is the simplification entailed, that treats first with an undersaturated oil reservoir as an integrated effect from the inception of production. The second phase pertains to saturated oil reservoirs that encompass a survey traverse. Although both methods of approach are necessarily different in aspect, the most interesting fact is that the mathematics so deduced are identical. Both the linear and radial water-drive systems are incorporated. for which an illustrated factual example is offered for the latter, treating with a saturated oil reservoir. INTRODUCTIO N What is performed in this work is the simplification of an involved computation by advanced analysis. Although such may be construed as a contradiction when one treats with higher mathematics; nevertheless, when direction is given to such an undertaking the results car. be most revealing. Likewise, it is to be mentioned that the bases for these mathematics have been developed on the expediency of the occasion. This is not to be inferred as a qualification of this work, but rather the demands frequently placed upon the author in his private prac- tice in meeting a time limit. A situation, instead of being fraught with hazards, often has given emphasis to creative thought. What will be entailed in this work is the simplification of the material balance formulas by the Laplace Transformation., Although this reveals entirely new horizons that will be given expression in a forthcoming tract, it suffices in the present instance to limit our attention to this phase of the development that treats both with an undersaturated and saturated oil reservoir. To orient the reader's thoughts as to what is involved in this simplification is the recognition that reservoir pressure, as such, is an inexplicit term in the material balance equation. This is the independent parameter that defines the total history of performance in the author's' unsteady-state water influx formulas, as well as the basis for the physical dependency of fluid behavior within the formation as prescribed in the Schil-thuis' material balance equation. Therefore, to isolate reservoir pressure, which is the most essential factor in any reservoir study, is rather a cumbersome procedure entailing either a trial-and-error calculation for the engineer; or as some prefer, a reiteration process performed on a computing device. However, once such an equation can be transcribed as a Laplace transformation, this inexplicitness so expressed can be alleviated to identify reservoir pressure as an explicit function of all the factors contributing to its change. This is the simplification encompassed, that will treat first with an undersaturated oil reservoir as an integrated effect from the inception of production, and secondly, with a saturated oil reservoir as a survey traverse. Although the two approaches are necessarily different because of the uhvsics involved. it is an interesting commentary that the mathematics are identical, showing the interdependency of the two methods. In order to acquaint the reader with this development, the simplest case will be treated first; namely, an under-saturated oil reservoir subject to a linear water drive. However, what may be construed for this example as an idealistic case is actually a most practical application in certain parts of the world, where the size of the fields are so large that radial water-drive approaches the configuration of a linear drive. Further, to avoid the repetition of much symbolism, frequent references will be made to the work of the author and an associate on Laplace Transformations3,

By M. R. Tek, K. H. Coats, D. L. Katz

A large number of boundary value problems encountcred in unsteady-state heat transfer, fluid flow through porous media, neutron diffusion and mass transfer involve the solution of a linear, parabolic partial differential equation commonly refcrred to as the diffusivity equation, where U is the dependent potential variable. K is the diffusivity (hydraulic, thermal, neutron, etc.) and 1 is the time variablc. Solutions to Eq. I are available in the literature for a wide variety of initial and houndary conditions. The great majority of these solutions are ohtained for geometric boundaries corresponding to linear. cylindrical or spherical flow models. A typical engineering application where the solution to Eq. 1 is required is the calculation of underground water encroachment across the boundaries of oil or natural gas reservoirs, In this particular area of application, the reservoir boundary is invariably approximated by circular geometry. However. the areal shape of many reservoirs can be better approximated by elliptic rather than circular boundaries. Thus the need for a general method of solving the diffusivity equation in elliptic coordinates arises in this problem as well as in other engineering applications involving elliptic boundaries. The solution to the diffusivity equation usually in volves the Error Function for the linear flow model, Bessel functions for the radial flow model and trigonometric or Legendre functions for the spherical flow model. It is well known that the general solution to the diffusivity equation in elliptic coordinates involves Mathieu functions. The significance of Mathieu functions in the analytical treatment of the diffusivity equation in elliptic coordinates is discussed in the literature However, these references do not providc analytical solutions useful in practical engineering problems. The objectives of this papcr are the development of the equations describing the unsteady-state liquid flow through a porous medium with an elliptic inner boun- dary, the development of a numerical method of solving these equations and, finally, a comparison of the water encroachment quantities calculated from the elliptic flow equation with those calculated from the radial flow equation. While the specific problem treated in this paper relates to unsteady-state liquid flow through a porous medium, the basic equations and computational techniques devcloped will apply equally well to problems occurring in the other areas of engincering interest mentioned previously. The solution given here is limited to a single case in which thc outer boundary encloses an area 100 times that of the inner boundary. DESCRIPTION OF THE FLOW MODEL Fig. 1 shows the flow model upon which the calcu lations presented in this paper are based. The inncr and outer houndaries of the flow model are represented by two confocal ellipses with major and minor axes, respectively, equal to 2a1, 2b1, and 2a, and 26,. The height of the elliptic cylinder flow model is denoted by 17. The following assumptions are employed in the development 01 the equations governing the unsteady-state liquid flow through the described flow model. 1. Uniform porosity and permeability throughout the flow model