Search Documents

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 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 S. R. Bardsley, S. T. Algermissen

Induction, nuclear and sonic well-logging methods were employed in a Green River formation oil-shale analysis and evaluation study conducted in northeastern Uintah County, Utah. The physical and chemical properties of an oil-shale section, previously cored and assayed by conventional methods, were used to evaluate the response of the various logs. The logging program was designed to measure variable properties of oil shale which relate to oil-yield potential in order that a relationship between assay-oil yield and log-determined properties could be identified, thereby permitting a direct determination of yield potential from logs alone. All of the major oil-shale zones and section ttlarkers are recognizable on the logs used in the study. The relationship between the response of the Density and Sonic logs and the assay-oil yield in gallons per ton was suficient to perinit the derivation of equations expressing the relationships. These equations can be used to determine the potential oil yield in gallons per ton of an oil-shale zone or section. The Neutron log response distinguishes the rich oil-shale intervals from the lean intervals, but does not appear to permit establishment of a quantitative relationship. The Garnma-Ray and Induction logs indicate only a qualitative relationship to oil-shale yields. Logging oil shale by Gamma-Gamrna Density and/ or Sonic logging methods will permit a fast, economical and accurate means of evaluating the potential yield of oil-shale deposits. Considerable sums of money have been allocated to Green River formation oil-shale evaluation during the past years by both private enterprise and Federal and state governments. The result of the work accomplished to date is commendable, but the task of fully exploring and evaluating one of the world's greatest reserves of potential energy for economical exploitation is enormous and much information is still needed. Present methods of evaluation consist primarily of sampling and assaying the oil shale for potential yield. The data received are excellent, but are both costly (inasmuch as it requires funds specifically allocated for the evaluation purpose) and time-consuming (as each representative sample must be assayed in the laboratory by specially trained personnel). At the present time, the Green River formation oil-shale province is one of the most active areas in the Rocky Mountain region for exploratory oil and gas drilling. This activity could play a twofold role and serve the function of spear-heading many oil-shale evaluation programs, as it could provide both cutting samples for assaying, geophysical data and modern well logs. This paper deals with a recent study conducted on Green River formation oil shale in which the application of induction, sonic and nuclear well logs to oil-shale evaluation was tested and proved. THEORY AND DEFINITIONS Oil shale may be described as siliceous marlstone, rich in solid organic matter called "kerogen". Kerogen is only slightly soluble in organic solvents, but it will decompose and yield oil vapors and gas when heated to destructive distillation temperatures at about 800 F. The mineral constituents of Green River formation oil shale are found in essentially uniform proportions with one another. The dominant types are, in the general order of abundance, dolomite, calcite, feldspars and quartz. Kerogen is also essentially uniform in its composition, and is about 80 per cent carbon and 10 per cent hydrogen by weight. The ratio of kerogen to the mineral constituents determines the richness or potential yield of the oil shale. This ratio also determines the total physical and chemical properties of the rock and provides a basis by which oil-shale richness may be determined from well logs. Because Green River formation oil shale was deposited in a lacustrine environment, the lithology is laterally consistent over wide areas. The vertical section consists predominantly of thin bands of alternating rich and lean oil shale which have been described by Bradley' as varves. The varves are occasionally broken by thin bands of volcanic ash and tuff and by zones of lean vugular oil shale. The effect that the thin non-oil shale beds and vugular zones would have on the well log analysis was a prime concern in the study. To cope with this possible problem, the logs used in the study were of two types: (1) logs believed to respond favorably to those properties of oil shale which are dependent upon the kero-gen-to-mineral ratio, these logs being the Gamma-Gamma Density log, the Sonic log and the Neutron log; and (2) logs that would respond to litho-logies and properties within the oil-shale section not associated with yield that might result in anomalous responses by the yield-measuring logs. Logs of the second type were the Gamma-Kay log, the Induction log and the Caliper log. THE PRESENT STUDY The present logging study was con-

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

By C. Scala, M. Gondouin

Published laboratory data have established that very significant streaming potentials can exist across mud cakes subjected to pressure differentials such as exist between a mud column and formation fluids. Since shale is similar in kind and in properties to a well-formed mud cake, it appeared possible that streaming potentials might also arise across shales. A practical consequence would be that the streaming potential contribution to the SP would be smaller than might otherwise be assumed since the SP is measured with respect to the "shale bare line". The resultant contribution would be the difference between the streaming potential across the shale and that across the mud cake. In order to investigate this possibility, a laboratory apparatus was devised for measuring streaming potentials under pressure conditions normally encountered in wells. Reported data on streaming potentids across mud cakes were verified. More important, appreciable streaming potentials were found across shale samples. As in the case of mud cakes, the shale potentials varied with the concentration and nature of the impregnating electrolyte. Dependence on pressure was In a different rela-tion than that for mud cakes, possibly because of the lower compressibility of shales. The differences in streaming potentials between the shales md mud cakes were generally small, although they might vary considerably in practice depending on the relative electro-kinetic properties of the particular shale-mud cuke combination. A "poorly eflective" shale with a "very eflec-tive mud cake" would produce the largest negative addition to the electrochemical component of the SP. Results of field experiments with pressurized wells were considered in the light of the laboratory results. INTRODUCTION Field experiments wherein the pressure of the mud column in a borehole was varied have often indicated the presence of a streaming potential, but its magnitude was not easily predictable and its sign sometimes changed from one formation to another in the same well.' In the interpretation of SP logs from standard surveys, nevertheless, the correct value of R,, is very frequently obtained, within the limits of field accuracy, using only the simple electrochemical expression of the SP provided that proper corrections are made for bed thickness, ion activity, etc.' Accordingly, it was apparent that the streaming potential component of the SP in a borehole was considerably more complex than commonly assumed and clarification was indeed necessary The close similarity in nature and constitution of a shale to a well-compacted mud cake suggested that a streaming potential might also exist across shales as pointed out by M. P. Tixier in 1953 and as recently emphasized by Wyllie." e practical consequence would be that the streaming potential across the mud cake would be canceled, at least in part, on electrical logs since the SP is recorded with respect to the shale potential as a base line. This appeared to be a plausible explanation for many of the effects observed but no supporting data were available. In 1955 preliminary experiments at low pressures were conducted in our laboratories which supported the hypothesis of the shale streaming potential and formed the basis for later studies at borehole pressures.* EARLY LOW PRESSURE EXPERIMENTS The early experiments were conducted primarily to determine the possibility of a streaming potential across shales. The shale used was a compact type designated as No. 7 from Reagan County, Tex. (Permian). Cut samples were saturated for over one month in .2 to 2.0-N NaCl solutions. In order to eliminate all electrode polarization effects, a dynamic method was used. A transient pressure of about 1 atm was applied across

By R. P. Alger, J. S. Wahl, C. W. Johnstone, J. Tittman

To overcome the problems of mudcake and hole irregularities, the new compensated formation density logging device employs two detectors spaced at different intervals from the source. The detector at the shorter spacing is particularly sensitive to the density of material immediately adjacent to the face of the pad. The contribution of this material, which includes mudcake and minor wall irregularities, affects the response of each detector to a different degree. The signals from both detectors are combined to give automatically a density correction which is added to the uncompensated density information from the detector at the longer spacing. Both the compensated density measurement and the amount of compensation are then recorded on the log. With the unwanted borehole effects removed, the measurement is recorded directly in terms of bulk density on a linear scale. INTRODUCTION Formation density is a very useful and diagnostic parameter for formation evaluation. When the matrix lithology is known, porosity is accurately and readily computed from density data.'-" hen the density measurements are used with other porosity-responsive measurements, such as Sonic and Neutron logs, both lithology and porosity may be defined with good accuracy."'" In the past, however, application of density logs was complicated by the need for corrections to obtain true formation densities from the log values. Mudcakes on permeable formations and roughness of the borehole walls have been particularly troublesome since, in each case, the device was prevented from direct contact with the formation. ~ttempts have been made to correct for such conditions by making a simultaneous caliper log, estimating pad standoff, and assuming knowledge of the composition of the intervening materials. While such corrections are frequently accurate, they are not always reliable because of uncertainties in the pad standoff and material composition. A further complication was the necessity for manual conversion of log counting rates to density. The newly introduced compensated formation density logging device (FDC) automatically corrects for mudcake and minor wall irregularities, and provides a log that is scaled directly in bulk density. The purpose of this paper is to discuss briefly the equipment, to explain the theoretical considerations involved in the compensation, and to present field examples of this new log. COMPENSATED FORMATION DENSITY EQUIPMENT The FDC tool incorporates two major advances over previous density logging devices: a second detector, spaced closer to the gamma ray source to provide increased sensitivity to material close to the pad; and the automatic derivation of a corrected formation density from the two counting rates. The arrangement of the source and detectors is shown in Fig. 1. The longer spacing detector is at the same spacing from the source and gives the same response as the single detector in the previous uncompensated device. By appropriate combination of the outputs of the two detectors, the FDC tool is able to compensate in large measure for effects of mudcake and hole rugosity. This is achieved, as will be described later, without specific knowledge of either the thickness or the nature of material between the pad and formation. To obtain corrected formation density, the two count rates are automatically processed by an analog computer in the surface panel. More specifically, the computer derives a value of uncompensated formation density from the long-spacing counting rate, computes a density correction p using both the long and the short spacing count-

Jan 1, 1965

By W. A. Nagel, R. P. Alger, H. Sherman, J. Tittmann

A sidewall epithermal neutron tool has been developed to substantially reduce environmental effects that have previously complicated neutron log interpretation. Designed for operation in uncased wells, the device provides increased accuracy in both liquid-filled and empty holes. A brief discussion of neutron moderation, diffusion and capture shows that logs using epithermal neutron detection depend on a smaller number of formation-characterizing parameters than those using thermal neutron or capture gamma ray detection. Thus, they come closer to providing an unambiguorls determination of hydrogen content. In this new device a directionally sensitive epithertnal neutron detection system has been incorporated in a side-wall source-detector skid to minimize borehole effects. The effects of variations in borehole size and shape, mud type, temperature and salinity are greatly reduced. Small residual borehole effects are then computationally accounted for in the surface control panel to provide a borehole-corrected neutron log. The log presents a direct recording of neutron-derived porosity on a linear scale. With corrections for borehole effects already applied, this direct recording of porosity simplifies log interpretation. Furthermore, comparison with a linear porosity presentation of a formation density log permits sandstones, limestones and dolomites to be readily identified. Thus, in complex or variable lithology, porosity is determined with greater accuracy and reliability than heretofore. Laboratory data and field results demonstrate the improvements in neutron logging afforded by this new side-wall epithermal neutron logging device. INTRODUCTION The new Sidewall Neutron Porosity (SNP) logging system, designed for use in uncased wells, provides reliability and accuracy never before achieved with neutron logs. The effects of variations in borehole diameter and shape, fluid salinity, mud weight and temperature—parameters that have long complicated neutron log interpretation— are suppressed or corrected for by this sidewall epithermal neutron detection system. Furthermore, to simplify interpretation the SNP log presents a direct recording of computed porosity on a linear scale. The performance improvements achieved by this new system arise primarily from the combination of two im- portant design features. First. though not unique to this new tool, an epithermal neutron detection system is used,',' Epithermal neutron detection substantially reduces the perturbing influences of the thermal neutron absorption properties of rock matrices and water salinity. Second, the neutron detection system is mounted in a directionally sensitive sidewall skid to greatly minimize borehole effects. The purpose of this paper is to explain briefy the advantages offered by the detection of epithermal neutrons, to describe the SNP equipment and the log, to present calibration and correction data and to give examples of interpretation methods made especially convenient by this new system. ADVANTAGES OF EPITHERMAL NEUTRON DETECTION Epithermal neutron detection, of all the neutron methods in current commercial use, provides the simplest determination of formation hydrogen content. Advantages of epithermal neutron detectors, as compared with thermal neutron and gamma ray detectors, are probably best explained by a brief review of the life of a source-emitted neutron. A fast neutron from the source will eventually be cap-tured by the nucleus of an atom. However. before capture is likely to occur, the fast neutron (energy of 100,000 electron volts (ev) or more) will be slowed down until it is in thermal equilibrium with its surroundings. The neutron is then considered to be a "thermal" neutron (average energy of 0.025 ev at 25C). The overwhelming majority of the slowing down thus occurs in reaching epithermal energies (energies just above thermal—0.5 to several ev). After reaching thermal energy. the neutron will "diffuse" through the formation with, on the average, no further energy change until capture. Upon the capture of the thermal neutron, relatively high energy gamma rays arc usually emitted. Hydrogen plays a very important role in the process of slowing down fast neutrons. On the other hand, the absorption effects of water salinity and low concentrations of other strong thermal neutron absorbers—relatively unimportant in the moderation process-—are very important in thermal diffusion, capture and the production of gamma rays. Thus, by detecting only epithermal neutrons with the SNP, spurious effects due to these thermal neutron absorbers are greatly minimized as compared to the effects on tools using thermal neutron or gamma ray detectors. The epithermal neutron detection system affords another important advantage. Log response in various litho-

Jan 1, 1967

By A. A. Stripling

The Cretaceous limestone wells of the Mara/Maru-caibo Dist. of Venezuela are extremely prolific producers. To maintain production on cessation of natural flow, large scale gas-lift operations were commenced involving high production rates on caring Flow. For the design of these gas-lift installations it was essential to have some knowledge of the pressure gradients involved in casing flow, so that the required injection pressures, optimum gas injection raws, etc., could be forecast. This paper presents a method of calculating these annular pressure gradients. Basically the method makes use of an energy balance equation coupled with an empirical energy loss factor derived from field data. The calculations have been put in a form suitable for "punch card-type" calculating machines. A set of gradient curves for Mara/Maracaibo conditions covering flow rates up to 12,000 B/D, and GOR's of 500 to 2,000 cu fd/bbl is presented. From these curves the flowing BHP can be predicted from surface data and the vertical pow performance under varying conditions of gas injection can be estimated. The method has been applied to La Paz and Mam gas-lift operations for over three years and has given consistently accurate results. INTRODUCTION The Cretaceous limestone fields of the Mara/Maracaibo area in Venezuela are characterized by the presence of widely distributed fissure systems. Under these conditions individual wells frequently have extremely high potentials and such wells are generally flowed on the casing-tubing annulus in order to reduce the pressure loss in vertical flow. To maintain production on cessation of natural flow, an extensive gas-lift system was planned. As a first approach to this gas-lift design problem, an attempt was made to construct a set of empirical flowing gradient curves from the available subsurface pressure data. Such a method has been used successfully for tubing-flow conditions in California by Gilbert'. Unfortunately, the empirical construction of casing-flow gradient curves presents additional difficulties, as pres- sure bombs cannot bc run in the casing-tubing annulus against high rates of flow. Thus, the Full flowing gradients cannot be defined by actual measurement and the best that can be done is to obtain spot pressures at the tubing shoe. Under these conditions a vast accumulation of BHP data is required to construct a set of gradients. To the La Paz and Mara fields the majority of the available flowing pressure measurements fall into a relatively narrow GOR range, in the neighborhood of the solution ratio. While there appeared to be reasonable prospects of constructing gradient curves for this limited range, measurements at higher GOR's, which were of major interest for gas-lift operations, were relatively few and the prospects of constructing empirical curves cor~spondingly remote. As an alternative approach, the possibilities of a mathematical analysis of vertical flow conditions were investigated. May and Laird2fV n 1934 presented an excellent analysis of the vertical-flow conditions in the Anglo-Iranian Oil Co.'s fields. Their analysis was based on an energy-balance equation and calculated the losses due to kinetic energy, friction and slippage separately. By combining these losses with the total calculated available energy, the required depth-pressure traverse could be deduced. On applying the method to La Paz conditions, good results were obtained for tubing flow at high tubing-head pressures. However, for lower wellhead pressures (300 psi and under) the discrepancy between calculated and measured values increased appreciably. It was considered that this was probably due to energy losses caused by slippage or liquid hang-up in a continuous gas phase, for which it was difficult to visualize any analytical approach. In 1952 Poettmann and Carpenter4 presented a method of pressure gradient calculation which, while based on an energy-balance equation similar to that employed by May and Laird, made no attempt to evaluate the various components making up the total energy loss resulting from irreversibilities of flowing fluids. Instead, they proposed a method of analysis which assumed that all significant losses for multiphase flow could be correlated in a form similar to that of the Fanning equation for frictional losses in single-phase flow. They then derived an empirical relationship linking measurable field data with a factor which, when applied to the standard form of the Fanning equation, would enable the energy losses to be determined. The total available

By W. A. Hoyer

Neutron activation analysis for mineral identification has attracted the attention of numerous investigators, judging from the publications on this subject in the past few years. 1-9 The appealing aspect of this method of mineral analysis is its proven sensitivity in the laboratory. In principle, activation analysis is one of the few petro-physical methods whereby a direct analysis for a mineral itself is possible in a borehole. The successful application of activation analysis requires the use of particular neutron reactions, measurement of selected gamma-ray energy intervals, and optimum bombardment and detection times. For copper analysis, Labanov et al., 10,11 using a low-energy PuBe source, reported an acceptable logging means. On the other hand, Pradzynski12, in analyzing ores and process tailings, used high-energy neutrons and photons for surface measurements of process samples. The work of Labanov et al. was based on measuring the rate of decay of the activated copper and hence was not as definitive as might be desired. Equipment is commercially available to carry out activation studies in the field 13,14 Gamma-ray spectra can be successfully transmitted over ordinary logging cables and, with borehole high-energy neutron accelerators and lower-energy chemical neutron sources, a variety of neutron activation methods can be studied. The purpose of the study discussed in this paper was to determine the feasibility and practicability of logging for copper using neutron activation techniques and commercially available equipment. The laboratory study showed that Cu could be detected in rock samples in the laboratory by high-energy neutron activation of 63 Cu. In the field, however, due in part to the large volume irradiated, the direct detection of copper using high-energy neutrons was not as promising. Detection of copper by measuring both the direct gamma-ray emission and the excessively strong low-energy gamma-ray absorption appears promising. Copper was not detected either in the laboratory or

Jan 1, 1972

By J. A. Franklin, E. Broch, G. Walton

To make the best use of rock core, the mechanical character of the rock should be evaluated concurrently with conventional geological observations. The advent of simple techniques, such as the point load 'nutcracker' test, makes field logging of mechanical properties competitive with sampling and testing in the laboratory

Dec 1, 1971

A combination of the needs of business, employee desires and safety and health considerations is stimulating increasing attention to the time spent by people at work and the timing of that work. Particularly we are seeing: - innovative work schedules for 24-hour mining businesses; and - response from industry to accommodate the needs, even the desires, of employees in return for systems which are more productive and more capital efficient. There appears to be an evolutionary trend for business management to become more sophisticated and responsive to the need for achieving close to best cost work schedules to sustain continuous improvement. The underlying motivating forces driving this trend for improved techniques in organising attendance at, and time spent at work, is the seemingly inseparable combination of technological advance and ever-increasing accumulation of capital. The mining industry has not remained untouched by this thrust for ever-increasing need for best use of labour, to adopt advances in technology and to provide attractive opportunities for putting accumulated capital to work. An essential component of this on-going thrust is a substantial effort by the mining industry to improve the performance of labour by providing most-productive work schedules.

Jan 1, 1996