Search Documents

By M. R. J. Wyllie

Eight laboratory-prepared aqueous base drilling muds representing common mud types, and 15 aqueous base drilling muds sampled in the field, have been used in an experimental investigation of the relationship between the streaming potential developed across a mud filter cake and the differential pressure causing filtrate flow. It has been established that the relationship is of the form, ES = kPy, where Es is the streaming potential, P the differential pressure and k and y are constants for any particular mud system at constant temperature. Beyond a certain minimum thickness this relationship is independent of filter cake thickness. Experiment has shown that both k and y are temperature dependent parameters and that in certain mud systems there is a critical value of P above which the simple relationship Es = kPY ceases to hold. The constant k tends to be a function of the mud resistivity at 77ºF, but the exponent y has not been found to be related to any commonly measured mud property. Certain field experiments are described which serve to demonstrate that the total SP (self potential) developed opposite a formation may be quantitatively accounted for as the algebraic sum of an electrochemical potential and an electrokinetic potential. A chart to permit the simple use of mean solution activities in the interpretation of electrochemical SP data is given and its use described. Various field methods for computing the electrokinetic SP component of the total SP are considered and rejected as being unduly complicated and inaccurate. It is shown that the employment of low resistivity drilling muds is essential if optimum accuracy is desired in the use of the S.P. curve for the purpose of computing connate water salinities by the electrochemical SP method. The use of such low resistivity drilling muds is accordingly recommended. INTRODUCTION In a previous paper' a quantitative analysis of the electrochemical compovent of the SP (self potential) curve was made and it was shown that the electrochemical* emf, E, could be expressed as E = ZtRT/F In a,/a,, where t is the transference number of the chloride ion, R is the gas constant, T the absolute temperature of the formation water and mud, a, the mean activity of the connate water (assumed to be sodium chloride) and a, the mean activity of the mud (as sodium chloride). The essential accuracy of this expression was confirmed by an investigation 2 of the SP kicks opposite sands of known water salinity. The logs examined had been routine runs in oil wells in a number of Mid-Continent states. Further confirmation of the original analysis of the electrochemical SP is implicit in the paper by Tixier3 wherein a relationship between the SP and the ratio of connate water and mud resistivities, based only on field data, is given. Allowing for the fact that mud and connate water resistivities were used instead of activities, and that electrokinetic effects were ignored, the constant reported by Tixier is essentially 2tRT/F. The error involved in the universal substitution of resistivities for activities will receive further reference below. Possible reasons for the fact that the electrochemical properties of shales do not appear to be significantly dependent upon either the nature or the quantity of cation exchange materials, such as clays, they contain has been considered at some length by Wyllie and Patnode.4 In the same paper, the probable effects of shale thickness and compaction on electrochemical efficiency, and hence the inherent difficulty of obtaining reliable laboratory-scale data, were also pointed out. However, it is

Jan 1, 1951

By Milton Williams

A comparison is made between the amount of interstitial water found by analysis in cores and that estimated from the resistivity curve of the electric log for the corresponding strata of producing reservoirs. It is shown that while the determined and the estimated water saturations are in good agreement in some cases, this agreement may be due to a fortuitous compensation of errors. Estimates of water saturations made from the electric 1og by use of existing relations should, therefore, be employed with caution. Also presented are data on the variation of salinity of interstitial waters in reservoir sands, and on the effect of variation of salinity and porosity on the validity of estimates of water content made from the electric log. INTRODUCTION Throughout the history of electric logging there have been recurrent cycles of optimism and pessimism as to the value of the electric log in the quantitative estimation of reservoir fluid content. The very early days of electric logging were an optimistic period, in which there was widespread faith in the ability of the resistivity curve to reflect fluid content, and in the ability of the potential curve to indicate textural characteristics. In time the accumulation of exceptions to this idealistic behavior brought about a general disillusionment, and caused the more pessimistic of those dealing with electric logs to believe that the sole importance of the log lay in stratigraphic correlation. Within the last several years optimism as to the quantitative use of the electric log has been revived, primarily as the result of the stimulating work of C. E. Archie' on the estimation of interstitial water in sands from resistivity data. Even though the general trend of thought on the quantitative value of the electric log reflects optimism, it is inevitable that there are some widely divergent opinions of the validity of quantitative estimates based on the electric log. A valid method of estimation of fluid content by means of the electric log and a general knowledge of the textural character of the reservoir rock would be of signal importance, not only in that it would save time and money in securing information on reservoirs now being drilled, but also because it would permit the evaluation of reserves in many fields drilled without adequate core data subsequent to the advent of electric logging. Because of this importance, Humble's Production Research Division, in cooperation with various other company divisions, initiated several years ago series of tests designed to evaluate the reliability of estimates of fluid content made from the electric log. The investigation was limited at the outset to sand reservoirs, since it was felt that the inclusion of limestone reservoirs would entail needless complications. In devising these tests it was immediately apparent that the chief difficulty in evaluating the worth of any method of estimation lay in establishing the true value of fluid content. The results of direct analyses of cores cut with ordinary mud were not sufficient, because of the probability of contamination of the fluids by filtrate from the drilling mud; the results of restored state tests on cores were deemed unsuitable because these results can hardly be said to furnish a primary standard. The method adopted was the direct analysis of cores cut with an oil-base drilling fluid. Such cores, when properly handled, are believed to afford the most precise determination of water content available at the present, and, in addition, lend themselves readily to the determination of various other properties which are pertinent to the estimation of fluid content from the electric log. Since it was desired to compare the fluid contents estimated from the electric log run in ordinary mud with cores cut in the same hole, the use of oil-base fluid necessitated the changing over, after all cores were cut, to a conventional mud. This procedure presented several advantages, one of which was that the use of oil-base fluid permitted

Jan 1, 1950

By M. R. J. Wyllie, H. W. Patnode

In the quantitative interpretation of electric log data it is essential to know the formation factor, defined as the ratio of the resistivity of the formation 100 per cent saturated with brine to the resistivity of the brine.' This concept implies that the formation factor for any given rock is constant irrespective of the resistivity of the brine or solution with which it is saturated. It is apparent that this can be true only when the rock matrix is wholly non-conductive. As reservoir rocks may contain clay or shale having an appreciable conductance one may expect these constituents to influence the determination of true formation factor. The formation factor may be theoretically derived from log data by dividing the resistivity of the mud filtrate invaded formation by the resistivity of the mud filtrate. Such a method of calculntion has implicit in it the assumption that the fluids in the invaded formation are wholly replaced by mud filtrate. that the solid matrix is non-conducting, and that the formation factor of the invaded zone is identical with the formation factor of the uninvaded zone. Since the validity of the first two basic assumptions regarding the insulating character of porous matrices and the displacement of the interstitial fluids in porous media do not appear to be fully established, the formation factors of a number of cores and clay slurries have been measured under conditions suitable for verifying the assumptions and the extent of any errors inherent in them. CORE STUDIES Apparatus and Procedure Resistivity measurements were made on short cylindrical permeable cores which had been previously mounted in lucite.2 The cores were clamped by means of insulated bolts between brass plates about 5 cm square and 1 cm in thickness. The plates were provided with a ring seal to prevent leakage between them and the core mount, and had a centrally located orifice so that fluids could be flowed through the mounted core. The brass mounting plates were used as electrodes and any current flow between them necessarily passed only through the constituents of the intervening core. Fluid saturation of the core, and fluid flow through the core, were effected by a centrifugal pump connected to the outflow end of the core, thus providing a differential pressure across the core of approximately one atmosphere. The procedure used was first to evacuate each core by means of the pump and then to saturate the core with an aqueous sodium chloride solution of known resistivity. The resistance of the core was measured and flow of the solution through the core continued for sufficient time to establish a constant resistance reading. After equilibrium had been reached with one solution, another sodium chloride solution of a different resistivity was flowed through the core and the flow continued until a constant resistance reading was again obtained. This procedure was continued until resistance measurements of the core had been made with a series of solutions within a salinity range from 1.0N NaCl to distilled water. The data obtained enabled the apparent formation factor of given porous rocks when saturated with fluids of different resistivity to be computed and also provided information regarding the ability of an invading fluid to displace interstitial water from a core. Results Berea Sandstone — The first core examined was a specimen of Berea sandstone." This core had been cut essentially parallel to the bedding for permeability studies. The core was clamped in the core holder and saturated with 1.0N NaCl solution and the resistance determined. It was subsequently flushed and the resistance measured when saturated in turn with 0.1N, 0.025N, 0.01N, .005N sodium chloride solutions and finally with distilled water. In flushing the core with each

Jan 1, 1950

By E. S. Mardock, Robert E. Bush

The object of this paper is to report recent developments in the quantitative interpretation of radioactivity logs. The use of reference lines is described in the application of the new zero radioactivity reference line for the determination of neutron derived porosities is discussed. Illustrated examples are shown of logs from the reef limestone of the North Snyder Field and from the Lansing-Kansas City limestone of central Kansas. A method of using an interrelation factor calculated from the neutron and gamma ray curves, called the neutron Productivity number, for the estimation of productivity from sandstone formations is discussed. When used in a relative manner the neutron productivity number has aided in Predicting whether a given zone will produce gas or liquid or whether it is too shaly to produce. Examples of logs from wells in the Freites District, Venezuela, and in East Texas are shown to illustrate the application of the neutron productivitv number. INTRODUCTION The petroleum industry was introduced to some of the preliminary results of investigations into the quantitative interpretation of radioactivity logs in October, 1949." It was stated at that time that the methods presented were entirely relative because they involved the use of reference points obtained from the radioactivity log, i.e., shale reference lines, dense zone reference lines, or a combination of both shale and dense zone reference lines. These relative methods offered many advantages to the industry in interpreting radioactivity logs in terms of porosity with particular application to limestone and dolomite reservoirs. However, these relative methods did not offer the wide range of applicability desired because the esta1)lishment of reliable reference points from shale in some areas were difficult. It was recognized that a more reliable mean; of establishing reference lines was vitally needed. Such a method is now available through the introduction of an instrumeutal zero device. It is the purpose of this paper to discuss the application of this new zero reference to the quantitative interpretation of radioactivity logs. Previously all applications of quantitative interpretations of radioactivity logs have been made in limestone reservoirs with little or no emphasis placed on sandstone reservoirs. A brief resumé of investigations into the quantitative interpretation of radioactivity logs in sandstone reservoirs is presented here. This discussion deals primarily with the interrelation of the gamma ray and neutron curves and the response of the latter to the density of the fluids filling the pore space in relatively homogeneous sandstones of uniform porosity. The authors believe that the methods presented herein are only applicable locally to specific fields and reservoirs, and that they probably do not offer a wide range of applicability. However, tile authors believe that the wider use of the new instrumental zero references in sand areas will reveal more information about the specific response of the gamma ray and neutron curves to changes in lithology and saturation. Two applications of quantitative interpretation using the new instrumental zero references are presented. Examples are given which show the comparison of logs and porosity in the limestones and dolomites of West Texas and Kansas. Logs of wells in 'Venezuela and East Texas are compared with the drill stem test data, production data, etc.. to illurtrate the interrelation of the gamma ray and neutron curve: and the rise of methods involving these interrelations to determine certain saturation variables. THE QUANTITATIVE APPLICATION OF RADIOACTIVITY LOGS TO LIMESTONE AND DOLOMITE RESERVOIRS Tile theoretical basis for the gamma ray and neutron curves has been discussed elsewhere.9 he method of recording the logs and apparatus involved has also been thoroughly covered and will not be discussed here.'."".'" It has been shown that the neutron curve responds exponentially to liquid filled porosity over the 1 1/2 to 35 per cent

Jan 1, 1951

By Robert E. Bush, E. S. Mardock

The object of this paper is to report recent developments in the quantitative interpretation of radioactivity logs. The use of reference lines is described in the application of the new zero radioactivity reference line for the determination of neutron derived porosities is discussed. Illustrated examples are shown of logs from the reef limestone of the North Snyder Field and from the Lansing-Kansas City limestone of central Kansas. A method of using an interrelation factor calculated from the neutron and gamma ray curves, called the neutron Productivity number, for the estimation of productivity from sandstone formations is discussed. When used in a relative manner the neutron productivity number has aided in Predicting whether a given zone will produce gas or liquid or whether it is too shaly to produce. Examples of logs from wells in the Freites District, Venezuela, and in East Texas are shown to illustrate the application of the neutron productivitv number. INTRODUCTION The petroleum industry was introduced to some of the preliminary results of investigations into the quantitative interpretation of radioactivity logs in October, 1949." It was stated at that time that the methods presented were entirely relative because they involved the use of reference points obtained from the radioactivity log, i.e., shale reference lines, dense zone reference lines, or a combination of both shale and dense zone reference lines. These relative methods offered many advantages to the industry in interpreting radioactivity logs in terms of porosity with particular application to limestone and dolomite reservoirs. However, these relative methods did not offer the wide range of applicability desired because the esta1)lishment of reliable reference points from shale in some areas were difficult. It was recognized that a more reliable mean; of establishing reference lines was vitally needed. Such a method is now available through the introduction of an instrumeutal zero device. It is the purpose of this paper to discuss the application of this new zero reference to the quantitative interpretation of radioactivity logs. Previously all applications of quantitative interpretations of radioactivity logs have been made in limestone reservoirs with little or no emphasis placed on sandstone reservoirs. A brief resumé of investigations into the quantitative interpretation of radioactivity logs in sandstone reservoirs is presented here. This discussion deals primarily with the interrelation of the gamma ray and neutron curves and the response of the latter to the density of the fluids filling the pore space in relatively homogeneous sandstones of uniform porosity. The authors believe that the methods presented herein are only applicable locally to specific fields and reservoirs, and that they probably do not offer a wide range of applicability. However, tile authors believe that the wider use of the new instrumental zero references in sand areas will reveal more information about the specific response of the gamma ray and neutron curves to changes in lithology and saturation. Two applications of quantitative interpretation using the new instrumental zero references are presented. Examples are given which show the comparison of logs and porosity in the limestones and dolomites of West Texas and Kansas. Logs of wells in 'Venezuela and East Texas are compared with the drill stem test data, production data, etc.. to illurtrate the interrelation of the gamma ray and neutron curve: and the rise of methods involving these interrelations to determine certain saturation variables. THE QUANTITATIVE APPLICATION OF RADIOACTIVITY LOGS TO LIMESTONE AND DOLOMITE RESERVOIRS Tile theoretical basis for the gamma ray and neutron curves has been discussed elsewhere.9 he method of recording the logs and apparatus involved has also been thoroughly covered and will not be discussed here.'."".'" It has been shown that the neutron curve responds exponentially to liquid filled porosity over the 1 1/2 to 35 per cent

Jan 1, 1951

By M. H. Waxman, M. Gondouin, H. J. Hill

Streaming-potential experiments were conducted within the Muddy- and Dakota-sandstone interval of a Denver basin well. Analysis of the data shows that, for this case, streaming potentials opposite sands were higher than those opposite shales when the mud was fresh. Streaming potentials opposite shales, however, were of sufficient magniiIlde to be important in SP interpretation. They were linearly dependent upon pressure differential, and they increased with mud filtrate resistivity. The magnitude of streaming potentials opposite sands was influenced by the characteristics of the original mud cake formed at the time of drilling. Use of mud samples taken at the time of logging may not always be satisfactory for determining streaming-potential corrections. Both the streaming and electrochemical potentials were found to be affected by the previous history of the borehole system. Following a change in mud characteristics, several days were required for SP stabilization. This was, in part, due to the mud cake acting as an imperfect shale membrane. INTRODUCTION Laboratory investigations by Wyllie,1 Salisch,2 Moore," Sen Gupta and Bannerjee,0 and Hill and Anderson' have shown that streaming potentials developed across mud cakes deposited on permeable filter beds are of sufficient magnitude to be critically important in SP interpretation. Streaming potentials across shales have been reported by Schenck,' Gondouin and Scala; and Hill and Anderson." The magnitude of the shale streaming potentials obtained by Schenck and by Gondouin et a1 led these authors to conclude that in many cases the shale streaming potential would approximately equal the mud-cake streaming potential. Hill and Anderson, however, obtained somewhat lower magnitudes of shale streaming potential. The field experiment described in the present paper was a cooperative project carried out by the respective operations and research departments of Schlumberger Well Surveying Corp., Shell Development Co. and Shell Oil Co. The primary objective of the experiment was to determine in the borehole the streaming-potential behavior of both shales and mud cakes under conditions approaching, as nearly as possible, those normally encountered in logging practice. METHODS AND EQUIPMENT WELL A Muddy- and Dakota-sandstone vie11 in the Denver basin was chosen for the experiment. Both sandstones were water-bearing and, since reservoir pressures in these sandstones are abnormally low (gradient of approximately 0.25 psi/ft), it was expected that studies of streaming potential could be conducted from low-pressure differentials (approximately 200 psi) to about 27000 psi by bailing mud from the hole, filling up and then applying wellhead pressure. The selected well encountered the Muddy sand-

By R. H. Winn

A method for the selection of the most suitable corrosion inhibitor for a particular system is given. The method involves the evaluation of surface passivity by means of copper ion displacement after the metal surface has been subjected to the corrosive atmosphere. The application of the method to the selection of suitable inhibitors for use in producing wells is described. Correlation with recognized laboratory tests has beet1 found to be good within the limitations of the test. The use of the test to evaluate inhibitor persistencc is described. INTRODUCTION Corrosion as it is associated with the industry shows wide variations in its qualitative as well as its quantitative aspects. This complicates the work of the corrosion engineer. In the past few years the efforts of the corrosion engineer have been further taxed by the large numbers of products being marketed as corrosion inhibitors by various manufacturers. These products are submitted to him together with a limited amount of data to support their use as inhibitors. The data are usually authentic and reliable yet all too often when these products are applied to an actual corrosive system the results are disappointing due to the ineffectiveness of the inhibitor and the corrosive system. The corrosion engineer must select a comparatively limited number of inhibitors for field trial from the large number available. The selection is usually made on the basis of a standardized laboratory test such as that adopted by the National Association of Corrosion Engineers or by field experience in an area; however, testing of this type is discouraging on large numbers of compounds, particularly in view of the specificity of environment. It has been an object of this work to devise a method which will allow the corrosion engineer to select the best inhibitor for a specific environment under field conditions. The test would of necessity be rapid and inexpensive. It would preferably eliminate the need for expensive and delicate instruments. STATEMENT OF THEORY The test as developed makes use of the fact that metal ions will be displaced from solution by metals having a more positive electrode potential. Thus, if an active iron surface is exposed to a solution of a copper salt, free copper will plate spontaneously on the iron surface as it is displaced from the salt solution by iron ions. If the iron surface is passified, as by coating with an insulating-type organic corrosion inhibitor, it is no longer capable of displacing the copper ions due to its own limited ionization. The result will be a change in the amount of copper plated upon exposure to a neutral copper salt solution. In practice it has been found possible to measure the magnitude of the alteration of the cell potential by determining the copper deposited on a standard iron surface under controlled conditions. This means that, as the film builds up, the resulting shift in the electrode potential tends to cause the iron surface to act more like a noble metal, resulting in a reduction in the cell potential to a point that copper will no longer plate spontaneously. Obviously, the salts of other metals such as silver, gold or platinum might be used for this operation; however, copper was selected because of the color difference of the plate, the comparatively low cost of copper salts and the ease of evaluation of the plated copper. Furthermore, its high cell potential provides a substantial driving energy. The determination of the copper is readily done colorimetrically by dissolving it in 5 per cent ammonium hydroxide solution containing 2 per cent potassium per-sulfate. The blue color of the copper ammonic complex ion so developed may be compared with standards prepared by dissolving known quantities of copper powder in the same reagent. In routine screening it is seldom desirable to determine the amount of copper plated. It is usually more important to determine the concentration of inhibitor required to prevent the plating of copper.

By P. Majani, F. P. Kokesh, M. Grosmangin

Determination of the quality of cementation of casing in oil wells in the past has involved inflow and circulation tests to insure that the producing zones were adequately sealed off from the adjacent zones. Existing logging methods, such as temperature and radioactivity surveys, may detect the presence of cement behind the casing. However, the qualities of the cement (i.e., its hardness and particularly its bond to the casing) are not indicated. The new logging method described in this paper operates on the principle that the attenuation of a sonic pulse transmitted by a casing is greatly increased when that casing is bonded to an outer annulus of hard material (such as set cement) which has an appreciably smaller sonic-wave velocity than that of the casing. The down-hole tool contains a source of recurrent sound pulses which are detected by a receiver spaced a few feet from the source The amplitude of the detected casing-borne pulse is measured, and the resulting signal is transmitted to the surface where it is recorded vs depth. Since amplitude is a function of attenuation, the log is readily interpreted. Laboratory studies have shown straightforward relationship between attenuation and such variables as source-detector spacing and per cent of circumference bonded. It is shown that cement not set or not bonded to the casing has compara- tively little attenuating effect. Field examples show not only the cement top, but also the variation in cementation quality below the top. Further, the increase of bonding with time and after squeeze cementation is depicted. The detection of poor cement jobs is confirmed by production tests and by formation-test results. It is anticipated that the method will have wide application in evaluating cementation quality prior to formation testing in completions and recompletions. The analysis it affords may aid in further improving cementation techniques. INTRODUCTION The main purpose of oilwell cementing is to isolate a production zone from other undesirable zones. To investigate whether this purpose has been accomplished, several logging methods have been used such as temperature logs, radioactive tracer logs, etc. While these logs all respond to the presence of cement behind the casing, they do not indicate the degree of bonding of the cement to the casing. Early in the application of sonic logging, it was noticed that considerable attenuation of sound signals takes place in cemented pipe and is often made evident on the standard Sonic log by cycle skipping.' The development of a circuit capable of continuously recording the amplitude of the casing-borne sound signal has made possible an extensive series of laboratory and field tests, which gave the following results. The amplitude of a sound signal after it has traveled in a firmly cemented pipe is only a small fraction of that recorded by the same device in free pipe. This provides a wide spectrum of energy levels; for given local conditions, empirical values of the amplitude can be correlated with the quality of cementation. Interpretations made in this manner generally have been confirmed by production tests, circulation tests and squeeze cementation. The purposes of this paper are to give a general description of the new logging method and to present some laboratory and field results. CYCLE SKIPPING Early attempts to study the quality of the cement behind the casing were performed with a standard Sonic log, which measures the transit time At, and were based on the well known phenomenon of cycle skipping.' Cycle skipping normally is interpreted as being a manifestation of weak signals at the receivers. The log of Fig. 1 was run with the recording instrumentation adjusted to enhance cycle skipping. A mirror-image presentation was used for better visual interpretation. The transit time of sound in steel is about 58 microseconds/ft (corresponding to a velocity of 17,000 ft/ sec), and the portions of the log where this value is recorded are interpreted as zones with no cement bond. Where cycle skipping produces a higher value of At, weak signals (or a high rate of attenuation) are indicated, and a good cement bond may be present. The occurrence of cycle skipping, however, depends too much on instrument adjustment to give a uni-

By Leendert de Witte

A review is given of the principles on which recently proposed methods of electric log interpretation in shaly sands are based and of the evidence brought up in support of the theoretical derivations. It is. pointed out that the qualitative techniques suggested to date have a tendency to be too pessimistic in the prediction of the presence of commercial hydrocarbon accumulations. In the quantitative treatment, a new concept is introduced, namely that of strongly reduced activity of the double layer counter ions which are present near the negatively charged rock surfaces. This lemmu when applied to the calculation of the equilibrium concentratiotls of interstitial waters yields a set of very simple relations. The resulting expressions for electrochemical potentials across rock samples appear in satisfactory agreement with laboratory experiments. The formulae obtained for the electrical conductivity of shaly formatioris are of the same form as those arrived at empirirally by previous workers. Combination of the expres.riorls for potentia1s and conductivities gives a direct proof of the Tixier relation which states that for shaly water sands the product of apparent formation water resistivity and apparent formation factor equals the resistivity of the sand 100 per cent .saturated with formalion water and which was verified by the work of Wyllie and Perkins and their co-workers. The relations for the resistivities and spontaneous porer11inI.s have been extended to the case of hydrocarbon hearing shaly formations, thus laying a formal basis for the quantitative interpretation of electric logs in shaly oil sands. INTRODUCTION The presence of disseminated clays in porous rocks saturated with electrolytic solutions has a strong influence on the transference of ions taking place under electrical or chemical potential gradients. Both the electrochemical emf's in evidence on the SP curves of electric logs and the electrical conductivity measured on resistivity logs arc directly dependent on the ionic transference in the interstitial waters of the formations traversed by boreholes. A large number of researches have been carried out on the nature and magnitude of the effects of clays on ionic transference. The most recent of these have been described in papers by McCardell and Winsauer,1 2 Wyllie and Southwick,3 poupon, Loy, and Tixier,4 and by Perkins et al.5 Most of these authors consider the effects of the clays as due to the fixed negative charges on the clay lattice. These charges cause the formation of diffuse surface layers with a high concentration of positive ions. The internal solution of the shaly reservoir rock has therefore a different ionic composition and concentration than the electrolyte with which it is in equilibrium (connate water produced from the formation or mud filtrate by which the formation is invaded). To calculate the equilibrium concentration of the internal solution, use is made of the Donnan6 relation, which states that the mean ionic activities of the internal solution and of the external electrolyte must be equal. From this and the

Jan 1, 1956

By G. R. Pickett

Examples are presented which show that the velocity~ amplitude, attenuation and apparent frequency of several acoustic waves can be recorded in the borehole. Examination of such recordings, termed "character" logs, indicates that the wave types observed include a refracted compressional wave and a wave which travels with formation shear velocity. Laboratory data are used to show that compressional and shear wave velocities are dependent on porosity, effective stress and lithology; but that the change in reciprocal velocity per unit change in porosity is larger for shear waves than for compressional waves. We, therefore, conclude that. the accuracy of porosity determinations can sometimes be improved by use of shear wave velocities, provided that the shear wave amplitudes are large enough to delineate the shear arrival from the preceding compressional arrival on the character log. Borehole data are presented which show that the difference between shear wave and compressional wave reciprocal velocities can be used to predict porosities. This is a refinement which may allow the prediction of porosities from single-receiver acoustic logs without introduction of errors from borehole fluid traveltimes. Laboratory and field data are presented to show that the relationship between compressional and shear wave velocities can be used to indicate lithology. An example is presented to show that fractures usually cause a greater reduction in borehole shear wave amplitudes than in compressional wave amplitudes, an effect which may offer a more reliable means of detecting fractures. The complexity of the borehole acoustic wave train can rake presently available cement bond logs highly sensitive to the gate and bias settings used. The character log offers a means to circumvent possible misinterpretations by recording all amplitudes, from which the interpreter can select the appropriate data for evaluating the cement bond. Character logs may also be used as a quality control for open-hole transit-time logs when existence of small compressional wave amplitudes interferes with the proper functioning of bias-controlled timing devices. Evaluation of the potential uses of character log data is not complete; but a character log presented in a form convenient for routine use would be a desirable addition to currently available logs. To summarize, possible applications for such a log in formation evaluation include the following (1) quality control of transit-time logs, (2) refinement of porosity predictions, (3) determination of lithology, (4) improvement of fracture detection and (5) improvement of cement bond evaluation. Suggestions are made regarding the requirements for a suficient but practical character log for routine use. INTRODUCTION Acoustic logs have become a widely used porosity tool in formation evaluation. In addition, there is a growing application of acoustic logs in cement bond evaluation and fracture detection. These applications have mainly involved the use of logs of first-arrival transit times and amplitudes and have not included detailed studies of the complete signal. The purpose of this paper is to show that significant benefits in formation evaluation can be gained by a more complete use of the acoustic wave train generated in the borehole by an acoustic logging tool. We hope that this discussion will also stimulate further development of logs suitable for routine use so that these benefits may be realized. Examples of acoustic wave train logs, termed "character" logs, are presented to show that several identifiable acoustic waves are present in the borehole. The measurable characteristics of these acoustic waves and some of their relations to formation properties of interest are also discussed. The more obvious potential uses of character logs are listed, and some suggestions are made regarding the requirements for a sufficient but practical character log for routine use. CHARACTER LOGS Some 10 years ago, Vogel' and Summers and Broding' noted that the signals received uphole from an acoustic logging tool located in a borehole had a number of interesting characteristics. The logging tool consisted of two or more pressure transducers spaced on an acoustically insulated body (Fig. la). One of the pressure transducers was used as a transmitter to generate pressure waves in the borehole fluid. The other transducer served as a receiver to detect any pressure waves reaching it in the borehole. The receiver then converted these pressure waves to electrical signals which were transmitted to the surface and displayed on an oscilloscope as a record of time vs receiver-signal amplitude. Fig. lb is a schematic representation of a typical record. The interesting characteristics seen in the earlier' and subsequent experiments were (1)

By M. Felsenthal

Engineers are frequently faced with the problem of having to predict oil recovery from a solution gas drive reservoir in the early life of a field. This is often the time when actual laboratory or field-derived gas-oil permeability (kg/k0) data are not yet available. Although excellent collections of k0/k0 data for various types of rocks may be found in the literature,1,2,3 a need exists for a more detailed correlation of k0/k0 data with core properties such as porosity, permeability, clay content or pore-size distribution. It is the purpose of this article to present this type of correlation for sandstone material. The correlation is based on laboratory tests performed on over 300 core samples from 19 reservoirs in the United States and Canada. EXPERIMENTAL Cores were analyzed by routine procedures, using samples of conventional size (7/8-in. diameter, 7/8-in long) cut parallel to the bedding planes. Air permeability was measured at 1.67 atm mean pressure in a Hassler-type apparatus. Porosity was determined with the aid of a Ruska bulk volume meter and a Boyle's Law rock volume meter, An evaporation method' was used to evaluate irreducible water saturations. Core samples were saturated with 25,000 ppm NaCl brine. The brine was reduced to the irreducible saturation, and the remainder of the pore space was filled with kerosene. This was followed by multiphase flow tests using the so-called "dispersed feed" technique.5 n this technique, gas and oil (kerosene and kerosene-saturated air) were dispersed with the aid of a porous Alun-durn feed head, the face of which contained a waffle-shaped system of Lucite-lined grooves. Oil entered the back side of the feed head and emerged at the raised areas of the front face. Gas was introduced directly into the system of grooves at the front face. The gas-oil mixture then entered the test core which was kept in capillary contact with the feed head. Permeabilities were determined from rate and pressure measurements and saturations from weights. All determinations were performed at equilibrium or steady-state conditions. Pore-size distributions were determined on selected samples by a mercury injection technique.6 GENERAL CORRELATIVE TRENDS When examining data from a given reservoir, it was noted that k0/k0 vs saturation curves became generally less steep as specific air permeability increased, a development that has a favorable effect on oil displacement efficiency by solution gas drive. This trend has also been reported by others in the literature.' An effect of porosity on k0/k0 data was also noted. This effect was not generally discernible in a study of relative permeability data for a given reservoir but became apparent when data for sandstone reservoirs of similar lithology but differing average porosity were compared. For instance, in a comparison of argillaceous and/or calcareous sandstones from 11 reservoirs ranging in average porosities from 14 to 28 per cent, a definite trend was noted, indicating that for a given permeability k0/k0 curves became less favorable (i.e., steeper) as porosity increased. A similar trend was observed for a group of comparatively clean sandstones from five reservoirs ranging in average porosities from 15 to 21 per cent. For a given permeability and porosity, comparatively clean sandstones gave more favorable k0/k0 curves than argillaceous and/or calcareous sandstone or chert reservoirs. The laboratory studies also showed that least favorable k0/k0 curves were exhibited by shaly sandstone, conglomerate, and sandstone containing carbon inclusions. Table 1 lists the various source reservoirs classified by lithology types. To facilitate a correlation of k0/k0 data with core characteristcis, it was found expedient to work with

By I. Fatt

A recently developed system for measuring electrical resistivity of liquids without use of electrodes offers some interesting possibilities in electric logging technology. The equipment as supplied by the manufacturer is satisfactory for continuous mud logging on a drilling rig or for measuring mud or filtrate resistivity in the laboratory. A simple modification of the commercially available instrument makes it suitable for measuring resistivity of core samples in the laboratory. The continuous measurement of mud resistivity on a drilling rig is a convenient means for detecting mixing of formation water with the drilling mud. Such information is useful to the geologist, the mud engineer and the logging engineer. However, continuous mud resistivity logging by conventional electrode-type resistivity cells is beset with difficulties. The mud, sand and rock chips abrade the electrodes, thereby changing the cell constant and eventually destroying the cell. Also, additives and crude oil in the mud may poison the electrodes by coating them with a nonconductive material. An electrode-type resistivity cell. therefore, may give erroneous readings under certain conditions. Electric logging companies circumvent the electrode poisoning problem by using a four-electrode resistivity cell for measurement of mud resistivity. In this cell, change in electrode area does not change the cell constant. However, the four-electrode cell is difficult to adapt for continuous reading and does not solve completely the problem of electrode abrasion by the sand and cuttings in the mud. The measurement of electric logging parameters on core samples in the laboratory encounters some of the same problems discussed in connection with mud logging. Ideally, the electrical resistivity of a core sample should be measured by placing platinum black electrodes in direct contact with the plane ends of a cylindrical or rectangular sample. Platinum black electrodes however, are much too fragile and easily abraded to be brought in contact with a rock sample. Also, oil or other constituents in the fluid contained in the sample will poison platinum black. In practice, gold-plated brass electrodes, in an AC bridge circuit operating at about 1,000 cps, are used for routine core analysis. For more precise work in research studies, a four-electrode scheme is used.',' Preparation of the samples for the four-electrode method is much too involved for routine core analysis. An apparatus for measuring resistivity of liquids without use of electrodes was described by Guthrie and Boys3 in 1879. They suspended a beaker, containing the electrolyte, by a torsion wire and rotated a set of permanent bar magnets around the vessel. The eddy currents induced in the electrolyte reacted against the rotating magnetic field to develop a torque, which was measured as a deflection of the torsion wire. In 1879 this method could not be made precise or convenient because of the lack of strong permanent magnets. The writer described a very greatly improved apparatus similar to that of Guthrie and Boys, but it was not suitable for continuous measurements or core samples.' Many electrodeless resistivity devices using radio frequency current are described in the literature.5, 6 These generally are suitable only for noting the end-point in a chemical titration. They do not measure resistivity, instead measuring a complex quantity which includes the dielectric constant and the magnetic permeability. The first description of the apparatus to be discussed in this paper was given by Relis.7 Improvements and modifications are described by Fielden,s Gupta and Hills,> and Eichholz and Bettens.10 DESCRIPTION OF APPARATUS The apparatus used in this study is based on the principle that the solution under measurement can form a loop coupling two transformer coils, as shown in Fig. 1. For a fixed AC voltage applied across Coil A, the voltage appearing across Coil B is a function of the resistance of the liquid-filled loop. The details of the voltage generating and measuring circuits are given in Refs. 7, 8, 9 and 10. A block diagram of the equipment is given in Fig. 7. Special features worth mentioning are the operating frequency of 18,000 cps and the automatic temperature compensation which results in the given resistivity readings being automatically correlated to 25°C. The liquid loop supplied by the manufacturer, shown in Figs. 1 and 2, was modified for use in core analysis (Fig. 3). The core sample under test is substituted for a section of the original loop. As shown in Fig. 3, the unit accepts only plastic-mounted cylindrical core specimens. A Hassler-type sleeve easily can be designed for the unit if unmounted cores are to be measured. EXPERIMENTAL PROCEDURE MUD LOGGING A simulated mud line was set up in the laboratory.

By J. P. Timmons, J. H. Moran, G. K. Miller, M. A. Coufleau

A prototype equipment has been designed and built for the digital recording of well logs on magnetic tape at the same time that the regular film recording is made. The format of the digital tape produced is such that it can be used directly at the input of the ZBM 704, 7090 or other models of ZBM computers which accept digital magnetic tape. This apparatus has been used for the experimental field recording of dipmeter tape logs which were subsequently computed by means of an ZBM 704 or 7090. In this paper the equipment and the digital tape are described briefly, and their application to the computer-interpretation of dipmeter data is discussed. A principal element in the interpretation of the dipmeter log is the correlation of the three microresirtivity dipmeter curves to determine the depth displacements between them. Several correlation methods for computer use are considered, with particular attention to their sensitivity to error and their consumption of computer time. The tape data were used to compute information content of the dipmeter microresistivity curves in terms of their frequency spectra. The results show that the sampling rate used in recording the digital information is quite adequate and illustrate a use of the digital tape in evaluating the characteristics of new tools. Some examples of field results are shown. It can be foreseen that, when digital tape recording becomes available for general field use, a whole new realm of possibilities will be opened up for the processing of other well logs through computations, which hitherto were not feasible because they were too laborious and time-con.sunzing. INTRODUCTION The last few years have seen a revolution in the design and production of data-processing equipment. Stored-pro-gram digital computers have progressed from a research curiosity to the basis of a major industry. There are now hundreds of such machines in daily use in the United States. With the acceptance of a technique that was, in fact, already clearly described by John von Neumann in 1945, the last decade has seen great strides in the development'of components, reliability, programming systems and, most spectacularly, in the sheer number of machines built and in use. In 1957 there were enough digital computers available to the oil industry to justify the suggestion that it would be worthwhile to investigate the possibility of using these machines in processing well log data.' The first result of this investigation was the appearance of what may be referred to as the input-output bottleneck. Well logs are customarily recorded on film. To get these data into a machine required then (and still does): a time-consuming semi-automatic reading of the film; conversion of the log data to digital form; and recording these digital data in some medium acceptable for computer input, such as cards, magnetic tape, or punched paper tape. However, the recording, reading, and re-recording could only result in deterioration of the data. Therefore, it was concluded that the fist step should be the development of a new, more direct recording technique supplemental to the film recording, which would provide easy access to the digital computer. There are many solutions to the problem of recording log data in an easily recoverable form. After careful consideration it was decided to adopt the boldest solution which, it was felt, was also the most elegant. It was decided to record well logs directly, in the field, on magnetic tape in such a way that this tape could be used without further modification as an input to the IBM 704 or 7090 computer. To realize practical field recording of magnetic tape logs, it became necessary to develop in a rather small package, an analog-to-digital converter, a tape recorder, and the necessary multiplexing and control circuits to allow the simultaneous recording of a multiplicity of logging signals. The magnetic tape recording was to be made simultaneously with the conventional logging operation in such a way as not to interfere with it. Along with the development of hardware, it was necessary to begin development of interpretation techniques and machine programs that would exploit the power of the digital computer. Here, again, there is a long list of possible applications. After much consideration it was decided to concentrate on the interpretation of the dipmeter log as a first application. It is the object of this paper to describe in some detail the developments sketched in the last three paragraphs.

By C. S. Matthews, H. M. Girner, C. D. Williams, M. J. Edwardson, H. R. Parkison

Quantitative interpretation of electric logs requires knowledge of formation temperature. In this paper, methods are developed for computing changes in formation temperature caused by circulation of mud during drilling operations. The basis of the method is the mathematical solution of the differential equation of heat conduction. The solution of this equation is presented in a series of graphs. These graphs are used to determine formation temperature disturbance at various radii for arbitrary mud circulation histories. Example comparisons with field results show reasonable agreement. It is concluded that, in general, the temperature disturbances caused by circulating mud are small beyond 10 ft from the wellbore but are quite significant near the wellbore. INTRODUCTION Quantitative interpretation of electric logs requires knowledge of the formation temperature in order to establish the resistivity of the formation water with accuracy. To determine the formation temperature, the temperature disturbances produced by circulating drilling mud must be evaluated. -The objective of this investigation was to develop a method for the numerical determination of these temperature disturbances at any distance from the wellbore as a function of time. The first step was to solve the differential equation describing the temperature behavior in the formation during and after mud circulation. This solution was used to calculate a series of curves relating temperature disturbance to shut-in time for various circulating periods. An "exact" method for computing formation temperatures with the use of these temperature-disturbance plots is described, and the results of an application of this method to a well in Montana are presented. Procedures for approximating formation temperatures are also discussed. BASIC EQUATIONS AND ASSUMPTIONS With the assumptions that (1) cylindrical symmetry exists, with the borehole as the axis, (2) heat flow is 1. The formation can be treated as though it is radially infinite and homogeneous in extent, as regards heat flow. 2. The presence of mud cake can be disregarded. 3. The effect of heat generated by bit action is negligible. 4. After mud circulation ceases, the rate of radial heat flow at the wellbore is negligible and is assumed zero for purposes of calculating temperature change with time. These assumptions are believed to be reasonable because of the agreement achieved between calculated and measured borehole temperature in a well, to be discussed later. Further, the agreement at the wellbore indicates that temperatures in the formation were also approximately correct. since it is the formation temperatures that determine wellbore temperature after ceasing circulations SOLUTION OF THE DIFFERENTIAL EQUATION AND ITS APPLICATION The solution of Eq. 2 for the temperature behavior during and after a single period of mud circulation is

By G. H. Pardue, L. H. Gollwitzer, J. H. Moran, R. L. Morris

Experimental and theoretical studies of casing-borne signals have provided data on the attenuation efect, compressive strength and shear modulus for different types of cement. These data resulted in the establishment of a definite relationship between the cement compressive strength and the Cement Bond log reading. The relationship is influenced by the casing dimensions, but is independent of the type of cement and the curing conditions. The theoretical analysis provides a mathematical statement of the phenomenon and shows that the attenuation rate is basically a function of the shear modulus of the cement and the thickness of the casing. The lab tests provided a connection between the shear modulus and the compressive strength, and permitted a comparison of the experimental and theoretical work. The results of the experimental and theoretical study have provided means of quantitative interpretation of the Cement Bond log. This interpretation is made more convenient with the use of a nomograph. INTRODUCTION Past interpretation of Cement Bond logs has been essenlially based upon a qualitative comparison between the log response in the zone of interest and the responses in the "fully cemented" and the "uncemented" intervals. Most of the variables influencing the magnitude of the signals were recognized in an earlier paper'; however, insufficient laboratory data were available at that time for quantitative evaluation of these variables. Early attempts to use the numerical amplitudes of the signals without full knowledge of their significance resulted in erroneous comparisons between various types of cement. An extensive program was undertaken to evaluate, both experimentally and theoretically, the attenuation phenomenon and to study the effects of the more common variables. A major objective of the experimental program was to obtain data on the relationship between cement-bond-tool response and physical property standards used by the industry so that practical interpretation methods could be devised. The studies included a comparison of attenuation rate with these properties as a function of time, temperature and pressure, using common sizes of casing. The main result was that a definite relationship was established between the compressive strength of the cement in inti- mate contact with the casing and the response of the cement bond tool. For a given casing thickness, regardless of the type of cement or the curing conditions used, the attenuation rate is a function of the compressive strength. A theoretical analysis showed that the attenuation rate for casing-borne signals should be a function of (among other things) the shear modulus of the cement and the thickness of the casing. The relationship between the compressive strength and shear modulus for curing times similar to those encountered in the field was determined from laboratory tests. Using this relationship, the computed attenuation rates showed reasonable agreement with the measured attenuation rates. The purpose of this paper is to report the results of the experimental and theoretical investigations and to present the first practical method for a quantitative interpretation of Cement Bond logs based on known variables. This method is applicable where the casing-borne signal is not affected by high-velocity formation arrivals. DESCRIPTION OF LOGGING TOOL The Cement Bond logging tool shown schematically in Fig. 1 has been described previously.' The tool consists of two sections: acoustic and electronic. The acoustic section contains a transmitter and a receiver, usually spaced 3-ft apart. Sonic energy emittted by the transmitter travels through the casing and is detected at the receiver. The electronic section measures the amplitude of the

Jan 1, 1969

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

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

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