Search Documents

By W. S. Ayres

By H. S. Thomsen, W. C. Duesterhoeft, R. E. Hartline

The attenuation and phase shift which formations produce in the electromagnetic field of an induction-type electrical well-logging instrument are great enough to substantially affect the response of the tool to formation conductivity under normal logging conditions. The application of the general equations for the propagation of an electromagnetic field in a conductive medium to the transmitter-receiver coil pair of an induction-logging-tool coil system gives an expression for the response which properly takes the propagation effects into account. A comparison of a calculation of this type with the response computed from the commonly used geometric factor concept, which does not include the propagation effect, leads to a factor Go representing the ratio of the responses computed by the two methods. The value of Go decreases not only with increasing formation conductivity, but also with increasing transmitter-to-detector coil spacing. For a single coil pair with a 40-in. spacing, the value of Go is 0.972 in a 20-millimho/m formation. The value is reduced to 0.915 and 0.812 for conductivities of 200 and 1,000 millimhos/m. As a result, the basic signal generated by the induction-logging coil system is not linearly related to formation conductivity as expected from the geometric factor concept. A uniform conductivity scale can be obtained only by adding a suitable nonlinear element to the recording system. The addition of auxiliary coils, having spacings less than the main-coil span, to achieve focusing results in an even greater departure of the basic signal from linearity with conductivity. The solution of the field equations near the interface between two formations of different conductivity gives the curve shape on an induction logging tool in crossing the interface. The addition of auxiliary coils to achieve focusing can add anomalous character to the curve shape in crossing the interface, which might be mistaken for lithological detail. Preliminary calculations which include the propagation effects in the determination of the conductivity correction for thin beds lead to correction factors which are substantially smaller than those obtained from geometric factor considerations. It is apparent that thin-bed corrections derived from geometric factor calculations are of doubtful validity. INTRODUCTION Although the electromagnetic induction type of electrical well log has been used for a number of years and has gained general acceptance in quantitative well-log interpretation, no technically complete investigation of the mechanism and of its operation has been published. The studies which have been presented are based upon the premise that neither the shape nor the intensity of the induction field of the logging tool is in any way affected by the electrical characteristics of the formation. From the inception of induction logging it has been recognized that under some formation conditions the performance of the instrument could be expected to be substantially different from that predicted by the simplified analysis. This anomalous performance has been called "skin effect" since it is the result of the same phenomenon which causes high-frequency alternating currents to flow only near the surface of metallic conductors. In the range of conductivities encountered in the sediments, the formation actually does not exhibit the sharply defined conducting "skin" shown by metals. The gradual modification of the distribution of the currents in the formation, produced by the induction field of the logging tool, is the direct result of phenomena associated with the propagation of the electromagnetic field through the formation. Therefore, it can be much more appropriately described as the "propagation effect". One purpose of this paper is to review the results of a solution of the equations for the performance of a single transmitter-receiver coil pair induction-logging system which takes into account the propagation effect, showing the manner in which the tool performance established by this analysis differs from that derived from the less comprehensive analyses currently in use. PRINCIPLE OF OPERATION OF THE INDUCTION LOGGING INSTRUMENT The high-frequency alternating current which is maintained at a constant value in the transmitter coil

By B. Wilhelmi, W. H. Somerton

Simultaneous measurement of pore and elastic properties of socks under a wide range of triaxial stress conditions may be made with equipment and methods developed in this paper. Tests were run on three outcrop sandstones at several confining pressures and axial stresses up to approximately 80 percent of failure stress. The greater reduction in porosities and permeabilities occurred in the hydrostatic loading portion of the stress cycle, but substantial reductions in both properties did occur upon application of deviator (triaxial) stress. Triaxial elastic moduli were calculated from strain gauge measurements. Both Young's moduli and shear moduli were found to increase with increasing confining pressure. Poisson's ratios increased with axial stress but showed no consistent change with confining pressure. Only very general relations were observed between pore and elastic properties. The lesser decrease in porosity and permeability at higher stress values is related to increase in rigidity of rocks at higher stresses. INTRODUCTION Many investigators have studied the effects of hydrostatic stress on the properties and behavior of rocks. A few have studied the effects of triaxial stress fields on individual physical properties of rocks. To our knowledge no one has measured simultaneously pore and elastic properties of rocks under a wide range of triaxial stress conditions. This is the purpose of this work. Rocks exist in the subsurface under what Jaegerl refers to as "polyaxial" stress conditions (all principal stresses unequal). When two of the three principal stresses are equal, this is referred to as triaxial loading. The triaxial method of loading is most commonly used in modern testing work because experimental difficulties in measuring physical properties of rocks are greatly reduced over polyaxial testing. Although triaxial loading may not reproduce subsurface conditions precisely, subsurface stress conditions are never known precisely anyway. Triaxial testing would seem to be more realistic than testing under hydrostatic loading (all principal stresses equal) for application of results to subsurface problems. The objectives of this work were threefold: (1) to develop methods of measuring simultaneously pore and elastic properties of rocks under triaxial stress conditions, (2) to determine whether changes in these properties are of sufficient magnitude to be considered in subsurface calculations and (3) to attempt to relate changes in pore properties and elastic properties since they will have been measured at the same time under identical loading conditions. REVIEW OF EARLIER WORK PORE PROPERTIES Several investigators have studied the effects of stress on the physical properties of rocks. Fatt and Davis2 showed the reduction of permeability of rock with hydrostatic loading. McLatchie et al.3 and Knutson and Bohor4 related permeability reduction to compessibility of reservoir rock, again under hydrostatic loading. Wble,5 and later Redmond,6 evaluated the effects of applied radial pressure on conductivity, porosity and permeability of sandstones. In this latter work, no axial stress was applied. Dobrynin7 reviewed much of the previous work done . on the effects of overburden pressure on physical properties of rocks. He concluded that changes in physical properties are controlled to a large extent by pore compressibility of the rock. Several sets of curves developed from empirical relations are presented. These curves permit estimation of change of properties with net overburden pressure, knowing the pore compressibility and porosity of the rock. Net overburden pressure is defined as (Pe - 0.85 Pi ) where Pe is the

By J. A. Putnam, A. D. K. Laird

In laboratory research on the behavior of oil, gas and water in porous materials, no direct method has been devised to measure saturation without disturbing the flow. Indirect methods involving various forms of radiation have been developed for measuring one or two components. The X-ray absorption method for two components has been described previously1-3. The possibility of three component measurement has also been mentioned1. Since two component measurement still requires considerable care to insure accuracy, the fact that three components can be measured satisfactorily is not generally appreciated. The fraction of the energy of an X-ray beam of one wavelength transmitted through a material is given by the transmission factor. in which 1 is the beam's path length through the fluid. p is the fluid density, and 1,. is its mass absorption coefficient. Eq. 1 implies that all parts of the beam pass through equal amounts of the fluid. Practical applications are, therefore, limited to systems in which the fluid is distributed so that I varies negligibly over the cross-section of the beam. Typical variations of the transmission factor with X-ray tube potential arc shown in Fig. 2 from which experimental points have been omitted for clarity. The X-ray absorption of low pressure gases in a core is practically zero so it cannot he distinguished from that a vacuum by instruments that will differentiate between liquids. The transmission factor of a gas in a core is. therefore, unity. Consequently, the saturation of only one gaseous component can be found. If two fluids completely fill a core, their two saturations add to unity and, therefore, only one value of the transmission factor, T, is needed to determine both of them. When two liquids do not completely fill the core, their saturations do not add to unity and transmission factors at two wavelengths must be measured. At one of these wavelengths, 7 for one liquid must be different from T or the other liquid. The change of T for one liquid us the wavelength is changed between the two values must be different from the corresponding change of T in the other liquid. The two values of fluid saturation will also give that of the gas because the three saturations must add to unity. In Fig. 2 it can be seen that there are many large differences of T between various oil solutions and water solutions at any one potential, and between potentials for any one solution. Consequently, conditions for optimum accuracy of measurement of the saturations of a gas and two liquids can be chosen. Once the operating potentials have been selected, however. it is more convenient to cross plot the data from Fig. 2 on transmission factor-solution concentration coordinates for the chosen potentials, as shown in Fig. 3. Fig. 3 was used to plan the experimental procedure for Core N-38-1. Since the voltage control at 45 kv was good, the slopes of the curves in Fig. 2 at this potential would cause little error. At low potentials, however, the voltage control was not good, so 33 kv was chosen because the small slopes of the curves made voltage control less critical. The lowest strength beam that could be measured reasonably at 33 kv corresponded to a minimum T factor of 0.385. Fig. 3 shows that both 4.80 weight per cent of cadmium chloride in water and 24.0 per cent of iodobenzene in crystal oil had this T value, and, consequently, were indistinguishable to the X-rays at 33 kv. Thus, at 33 kv, the liquid saturation was given by the same curve regardless of the relative amounts of oil and water solutions present. This simple calibration curve from which the gas and liquid saturation can be read, is shown in Fig. 4. Fig. 3 also shows that at 45 kv, T for the aqueous solution is 0.610. and that T for the oil solution is 0.220. These values were transferred to Fig. 4 to give the triangular calibration plot from which the brine and oil saturations can be found, after the gas saturation has been determined at 33 kv.

By O. D. Gaither

This paper is an analytical study of the flow of fluids through small vertical conduits. Small conduits are defined as 11/4-in. nominal diameter tubing size and smaller, and approximately twice this area for annular conduits (i.e., 1- X 21/2-in. annulus and smaller). Experimental data are presented for the 1-X2-in. and 11/4- X 2%-in. annuli, and the I-in. and 11/4-in. tubing, since these represent the small conduit sizes and configurations generally encountered in oilfield applications. Data have been gathered for these conduits for single-phase water, single-phase gas and two-phase water-gas mixtures, with particular emphasis on high gas-liquid ratios. Water rates in excess of 2,000 BID and gas rates in excess of 2.5 MMcf/D, and two-phase flow ratios in between these two, represent the scope of the data gathered. Existing equations have been applied to predict flowing pressures and compared with experimental data. New correlations have been developed. INTRODUCTION The increased economic pressure on the domestic oil industry in the United States has constantly required the use of new techniques and equipment designed to reduce the cost of finding and producing oil and gas. Since tangible items are most readily apparent in economic analysis, the advent of lower-cost well completions was inevitable. One of the methods used to reduce costs which has received widespread attention is the slim-hole completion technique where tubing is used as the well casing and in which small conduits are used for tubing if necessary. Small conduits, defined by Kirkpatrick1 as "11/4-in. diameter nominal tubing and smaller for tubing flow and less than twice the 11/4-in. diameter nominal tubing internal flow area for annulus flow", have also found widespread usage as siphon strings for de-watering gas wells and as "kill" strings in deep high-pressure oil and gas wells. The growing use of small-diameter tubing has resulted in an increased need for development of improved methods to measure or predict flowing bottom-hole pressures since the physical dimensions generally preclude the use of subsurface-recording pressure gauges. Even in the cases where small bombs are available, the relatively high velocities encountered at nominal flow rates make it necessary to use excessive weight bars or special hold-down devices. Attempts to use recognized correlations to accurately predict flowing or gas-lift performance in wells equipped with small conduits have been generally unsuccessful. Insufficient field data were available to allow the development of a correlation on this basis, and an experimental approach was applied in an attempt to obtain a workable relation. The experimental approach used to obtain the data presented in this paper was actually a compromise between a field installation and a laboratory study. A test well 1,000 ft in length was used to obtain flow data on single-phase liquid, single-phase gas and two-phase water-gas flowing mixtures. Liquid rates up to 2,200 B/D and gas rates up to 3 MMcf/D were used in the single-phase flow studies. Two-phase flow rates from 100 to 600 B/D with gas-liquid ratios from 500 to 8,000 cu ft/bbl were recorded. Experimental data were obtained for single- and two-phase flow through 1-in and 11/4-in. nominal tubing, and through the annuli between 1- and 2-in. and 11/4- and 2%-in. nominal tubing strings. Experimental results for the two-phase flow are compared to the Poettmann-Carpenter correlation2 which is widely used as a comparative standard for development of multiphase flow predictions in flowing and gas-lift wells. Correlations developed by Tek,3 Baxendell and Thomas" were also investigated. The experimental data recorded herein fell in between the two flow regimes as defined by Ros," and this correlation also failed to yield satisfactory results. The fact that existing correlations failed to confirm the experimental data led to the need for development of a new correlation. Although a two-phase flow study was the primary objective of this investigation, data were also recorded for single-phase flow of water and gas, and constants were developed relating to pipe roughness and equivalent diameters for annular flow. These single-phase studies assisted materially in the development of certain of the two-phase flow results. Considerable previous work has been published which presented relationship of surface measurements to bottom-hole condition. The works of Buthod and Whiteley,6 Jones,' Poettmannb and the Texas Railroad Commission" are classic examples of the successful use of mathematical relationships which allow acceptable predictions of subsurface pressures, when gas is the flowing fluid. Darcy and others have derived relationships which may be used with minor modifications to predict subsurface flowing conditions in injection and water-supply wells. As previously stated, the application of the single-phase flow relationships

By Persifor Frazer

In his " Final Report," Vol. 11, part second, p. 763, Prof. Rogers sums up the metalliferous veins of the mesozoic sandstone by remarking that these are not associated with dykes or trap-rock, but are independent metalliferous injections. But singularly enough, in enumerating the different kinds of ores, no mention is made of

By C. F. Nagel

This paper is written primarily for those who are familiar with the processes mentioned but who desire a further insight into some of the fundamental principles. It does not give a complete description of the various practices employed in the working of aluminum and its alloys, but presents a few illustrations of the essential respects in which aluminum differs from other common metals, and emphasizes the controlling factors of the processes under consideration. Heat Treating Certain aluminum alloys are capable of having their mechanical properties greatly altered by a heat-treatment process. The degree of this irnprovement may be illustrated by some typical properties of dura-lurnin (Table 1). It is essential that the material be in the wrought condition in order to obtain rnaximum heat-treatment effects. In the case of duralumin, the process consists of heating the metal to a prescribed temperature and then quenching it. Immediately after quenching, the material is relatively soft though not as soft as in the annealed state. If allowed to stand at room temperature, the material changes, increasing in hardness and strength without loss of elongation, but at the expense of the plasticity. This action is very rapid at first, gradually diminishing in intensity and to all practical purposes being completed in about four days. This process of hardening is termed

By Oswald J. Heinrich

The great improvement which have been made in late years in the different systems and instruments used to perforate the crust of the earth for purposes of testing and exploring for mineral resources or for other reasons, by means of bore-holes of greater or less diameter and to greater or less depths, are manifold. It may therefore be of interest to make an impartial comparison of them, or of such as are most commonly used. Being well aware of the difficulties which are encountered in such a task, a reasonable allowance most be made in such a comparison. It is a rare instance when even two decidedly different systems or instruments have been used in the same locality, penetrating identically the same strata to the same depth. The diameter of a borehole as well as its greater or less depth must also be considered in regard to its cost per foot, etc. But for practical purposes, for exploration at least, we may assert this rule: Any bore-hole, of whatever diameter, which allows the means of examining the strata most carefully and completely, will be sufficient for this purpose, but that one which will do it most economically and in the least space of time, must be considered the best for it. It will, therefore, partially depend upon circumstances, what system, in spite of other superiorities, may be the best in each case, and the following remarks may aid in deciding this question. To make the attempt here to review all that has been published upon the different systems and means of boring would be out of place. Reference can be had to the able publications treating upon this subject. A difference in this respect exists in regard to the use of the hollow rods in connection with the diamond drill, used for that purpose. The great future, which appears to be evident for this instrument, will be warrant enough for any one, who has had the opportunity to test it, to give his experience to the public. Only in this way can we bring this system to such perfection, as it seems to deserve and

By Alexander L. Holley

IN order to put sufficient work on steel ingots for rails, they must be reduced from about 12 inches square. As this cannot be done at one heat, they are first drawn down to about 7 inches square, and then reheated and rolled into rails. This first reduction or " blooming," is usually done in this country in a 30-inch 3-high rolling-mill, with movable rolls, so as to get several reductions in each set of grooves. The first of these mills has been running at the Troy Steel-works above a year, with great success. Another at the Cambria Iron

By J. Blodget Britton

To the question, " Do the Pennsylvania anthracites contain phosphorus?" asked at the last meeting of the Institute during the discussion on the metallurgical value of Western lignites, I can now give a definite answer. Through the kindness of a gentleman connected with the Philadelphia and Reading Railroad, and another connected with one of the principal ironworks on the Lehigh, there was sent to the Ironmasters' Laboratory a number of samples of the coals now largely shipped from the Schuylkill and Lehigh regions. Seventeen of these I examined for their contents of ash and phosphorus, and obtained results as follows : Ash. Phosphorus. No. 1, . 4.10 .0063 2, . . . . 7.97 .0209 113, . 5.53 .0049 114, . . . . . . . 6.55 .0056 . 5, . 5.39 .0761 . . . . . . . 11 6, 14.15 .0593 7, . 10 43 ,0049 rl 8, . . . . . . . 5.29 .0354 9, . . . . . . . 3.01 .0042 6L10, . . . . . . 3.03 .0042 1 1 . 275 .0014 12, . 2.82 ,0111 l'13, . .3.69 trace L1 14, . . .2.35 trace 15, . . . . . . . 7.22 .0014 16, . . . . . . . 3.49 .0167 17, . 4.63 trace

By L. K. Nemeth, H. T. Kennedy

The investigation presented in this paper was performd to develop a relationship between the dewpoint pressure of a hydrocarbon reservoir fluid and its composition, temperature and characteristics of the heptanes-plus fraction. This is part of an over-all program to develop correlations for the essential properties of reservoir fluids with the composition, temperature and pressure. In all, 579 dewpoint pressures, determined experimentally on 480 different condensate systems along with their compositions and temperatures, were obtained from the literature and from previously unpublished sources. Compositions were expressed in terms of mol fractions of methane through heptanes-plus, nitrogen, carbon dioxide and hydrogen sulfide, and the molecular weight and specific gravity of the heptanes-plus fraction. Distribution of the samples assured a good representation of most of the major oil-producing areas in the world. A multiple regression analysis program written for the IBhl 7094 digital computer was used to compare numerous equations and to calculate their coefficients. The final equation containing 11 coefficients resulted in an average absolute deviation of 7.4 percent for all the samples. Segregation of data on the basis of one or more variables did not materially improve the results. The correlation developed is useful in predicting the dewpoint pressure of condensate reservoir fluids from variables normally measured in the laboratory. INTRODUCTION In any discussion of the dewpoint of hydrocarbon mixtures, it should be noted that there are two kinds of dewpoints which must be distinguished. The first type, or normal dewpoint, usually occurs at low pressures as dry gas is compressed to the point where liquid first forms. The second type, sometimes called the retrograde dewpoint, occurs when a gas mixture containing heavy hydrocarbons in solution is depressured until a liquid is formed. The first type, as applied to condensate mixtures, would normally be below atmospheric pressure and would be of no interest in reservoir performance. The second type is of great significance in interpreting reservoir behavior and is the property referred to throughout this paper. Laboratory measurement of dewpoints is laborious, expensive and subject to many errors, and it is desirable to find a reasonably accurate method to predict dewpoint pressures from the composition and other readily available properties normally measured on the fluids. A possible method of calculating dewpoint pressures would be to employ a set of equilibrium ratios (K-values). 1 However, the K-values themselves depend on composition as well as temperature and pressure; presently their use adds nothing to the accuracy of predicting dewpoints. If future correlations provide a significantly better method of estimating K-values for heavy hydrocarbons, they may be of greater use. REVIEW OF LITERATURE In 1944, Olds, Sage and Lacey 2 studied the behavior of fluids from the Paloma field. Properties of seven mixtures were studied at temperatures ranging from 100 to 250 F and at pressures up to 5,000 psia. From the information obtained, the influence of composition and temperature upon the dewpoint was determined. It was found that removal of the intermediate molecular weight components from the mixture resulted in considerable increases in the dewpoint pressures. Results indicated that the effect of changes in temperature was relatively minor when compared with the effect of modifying the mixture by removing the intermediate components. Thus, removal of propane and butanes resulted in an increase in dewpoint pressure more than 10 times greater than a change in temperature from 100 to 250 F. The authors stated that althouah

By P. A. Witherspoon, R. W. Donovan, T. D. Mueller

The use of petroleum-barren aquifers for underground storage has become extremely important to the natural-gas industry. A critical problem in assessing the feasibility of a specific aquifer for such use is the permeability determination of the caprock over the proposed storage project. The approach used here is to conduct both static and dynamic field tests on the aquifer being analyzed. Valuable information on the possibility of communication between the storage aquifer and any other aquifers above can be obtained by measuring hydrostatic water levels and water analyses. Significant differences in such data give evidence of the lack of communication between the intended storage reservoir and other horizons. The dynamic approach requires that one well be pumped in the storage aquifer, and changes in fluid levels recorded in both the aquifer and its caprock. The interpretation of the data from such pumping tests involves the solution of nonsteady radial flow in an infinite aquifer and the influence on such flow of a leaky caprock. A finite-difference model has been used to investigate this problem, and the transient behavior has been solved numerically with a digital computer. It has been found that the pressure transients in the storage aquifer are not affected significantly by moderate caprock leakage. The pressure behavior of the caprock is a much better indicator of the degree of leakage, and generalized solutions for this behavior are included. Field data are presented to demonstrate both the static and dynamic approach. If is concluded that appropriate investigation of the groundwater hydrology in an aquifer-type gas-storage project can provide much valuable information for determining the effectiveness of the caprock to hold gas. INTRODUCTION Underground storage of natural gas in the United States has been developing at a rapid rate over the past few years. In 1955, the total gas-storage capacity was about 1.6 trillion cu ft; by 1961, this figure was almost 3.2 trillion cu ft, an increase of 100 per cent in six years.' This trend un- doubtedly will continue because the economics favor the development of gas storage, as opposed to the construction of new pipelines, to meet the inherent cyclic demand for fuel in the metropolitan areas of this country.' About 15 per cent of the current underground gas storage has been developed in petroleum-barren aquifers, i.e., geological domes or anticlines in which no commercial quantities of oil or gas had been produced prior to the storage operations. The necessity for using barren aquifers outside many metropolitan areas of this country has been due to the lack of depleted oil or gas fields that were near enough and large enough to meet the demands of such consuming areas. Pipeline companies have developed aquifer storage along their transmission lines to meet the fluctuating needs of their complex systems. Considerable thought has also been given to the problem of storing gas in a structureless aquifer, both in this country' and in the Soviet Union outside the city of Leningrad.'," Conditions such as these have led to the development of aquifer gas-storage projects in many parts of the U. S. Most of these developments have centered in the Mid-Continent area, and the greatest amount of activity has been concentrated in Illinois.6 Thus, the use of petroleum-barren aquifers for gas-storage purposes has become extremely important to the natural-gas industry. There are three basic problems in developing aquifer-type storage: (1) finding an adequate geologic structure, (2) finding a suitable storage reservoir within the structure and (3) determining the tightness of the caprock over the intended storage zone. The first two problems can be solved by applying conventional methods of exploration geology, but once these problems are solved, the question arises as to why no oil or gas is present in an otherwise favorable setting. Two situations are possible: (1) an adequate source bed was never present, or (2) a source bed was present but the petroleum seeped away because of a leaky caprock. Determining the tightness of the caprock is one of the most critical problems in assessing the feasibility of a specific aquifer for storage purposes. In attacking this problem, one usually takes cores of the caprock and subjects them to a rigorous investigation. Such core data are desirable, but they only detail the matrix properties and cannot be expected to reveal the gross characteristics of the caprock. Several gas-storage projects in the U. S. have had considerable leakage where

By C. C. Mattax, J. E. Bobek, M. O. Denekas

ABSTRACT investigations in recent years have shown that rock wettability can exert a profound influence on the displacement of oil by water from oil producing reservoirs. Core analyses frequently show oil recoveries from preferentially water-wet rock to be significantly greater than those from preferentially oil-wet rock. Thus, an accurate prediction of water-drive or waterflood oil recovery is dependent on the evaluation of reservoir rock wettability. Laboratory and field studies have been undertaken to examine and develop methods which may be used to measure rock wettahililiy, to investigate factors which may alter the wettability of reservoir rock, and to determine the wettability of specific reservoirs. It has been found that measurement of the rate and volume of spontaneous imbibition of the wet-ting phase by a rock is a reliable and reproducible test for semi-quantitarive determination of preferential rock wettability. Laboratory tests have shown that the wettability of reservoir rock may depend on both the crude oil composition and the rock type. Field and laboratory tests have indicated that coring flulids and core handling techniques can cause significant changes in the wettability of rock surfaces. However, a few fluids, brine in particular, appear not to affect core wettability and may be used when coring to determine reservoir wettability. Core handling and packing procedures have been developed which preserve core wettability during storage and laboratory testing. The wetfability characteristics of six reservoirs have been evaluated by wettabilify tests performed in the field on fresh cores cut with brine or special water-base muds. The tests indicated that all six reservoirs are water wet; however, the de-gree differed from one reservoir to another. INTRODUCTION Although surface properties of reservoir rock and interfacial characteristics of rock, crude oil, and water systems have been discussed for many years, the full significance of these properties to the production of oil is still in the process of evaluation. It is known, however, that one surface characteristic—preferential wettability of the reservoir rock—-is of paramount importance when the reservoir is being produced by water-flood or water-drive mechanisms. The importance of rock wettability has been pointed out by several authors'.',%* who have found that pref- erentially water-wet cores flood more efficiently than oil-wet cores; that is, more oil is recovered from water-wet cores in the early flooding stages than from oil-wet cores. Other authors'," indicate that waterflood oil recovery from cores of intermediate wettability may be greater than that from either strongly water-wet or strongly oil-wet cores. These recent investigations emphasize that the true wettability of a reservoir must be known before the performance of waterflood or water-drive reservoirs can be accurately evaluated. When reservoir wettability is to be determined by core analysis, precautions must be taken to maintain the original wettability which can be altered readily during coring5 nd testing.' In the coring operation the core may be partially or completely penetrated by the drilling fluid which, if it contains surface active materials, may drastically change the wettability of the core. Also, improper core handling during stcrage and testing may change its wettability because of evaporation of fluids or exposure to oxygen or surface active contaminants. One purpose of this paper is to discuss methods for obtaining, preserving, and testing cores without altering their original wettability. In addition, the results of wettability tests on fresh and preserved cores from six reservoirs are presented where special precautions were taken to obtain samples without altering their wettability during the coring operation.

By P. M. Blair, D. W. Peaceman

The shape and position of the gas-oil transition zone during downdip displacement of oil by gas has been calculated using flow equations which include the effects of gravity, relative permeability, capillary pressure and compressibility of the fluids. The calculations treat the problem in two space dimensions, and results are compared with data from a laboratory model tilted at 30° and 60° from the horizontal on displacements near and above the maximum rate at which gravity segregation prevents channeling of the gas along the top of the stratum. The good agreement between calculated and experimental results demonstrates the validity of the tecbnique as well as that of the flow equations. INTRODUCTION Knowledge of the fluid distribution and movement in an oil reservoir is important in producing operations and estimation of reserves. The history of the oil industry has included steady progress in improving the accuracy of calculations which provide the required knowledge. The earliest method of calculating reservoir performance consisted of material-balance equations based on the assumption that all properties were uniform throughout a reservoir.l For many reservoirs such a simple formulation is still the most useful. However, when large pressure and saturation gradients exist in a reservoir, the assumption of uniform values throughout may lead to significant error. To reduce these errors, Buckley and Leverett2 introduced a displacement eqation which considers pressure and saturation gradients. Methods available at that time permitted solutions to the Buckley-Leverett equation in one space dimension; these solutions have been very useful in solving many problems related to the production of oil. However, the one-dimensional methods are not adequate for systems in which saturations vary in directions other than the direction of flow. An example of such a system is the case of gas displacing oil down a dipping stratum in which the gas-oil contact becomes significantly tilted. Of course, the Buckley-Leverett displacement method cannot predict the tilt of the gas-oil contact. Recent improvements3,4 of the one-dimensional Buckley-Leverett method achieve some success in predicting the tilt of the gas-oil contact at sufficiently low flow rates. However, at rates high enough that the viscous pressye gradient nearly equals or exceeds the gravity gradient, even these improved one-dimen,sional methods incorrectly predict the shape and velocity of the contact. Further progress in estimating such fluid movements in a reservoir appears to require consideration of the problem in more than one space dimension. The recent two-dimensional method of Douglas, Peaceman and Rachford5 appears adaptable to calculate changes with time of the saturation distribution in a vertical cross-section of a reservoir. The movement of saturation contours should represent the moving fluid contacts and include the effects of crossflow due to gravity, as well as variations in the rock and fluid properties. The nonlinear nature of the equations used in the method has prevented proof of the validity of the solutions. Douglas, Peaceman and Rachford made some comparisons with experiment but did not include cases in which gravity was important nor cases involving displacement by the nonwetting phase.For these reasons, a test of the two-dimensional method for a case in which these factors are included would be very desirable. The test selected was a comparison of calculated results with those from a carefully controlled laboratory experiment on a model with measured physical properties. The model selected was one in which gas displaced oil down a tilted, rectangular sand pack. The model can be thought of as representing a vertical cross-section taken parallel to the dip of a reservoir. The displacement thus simulates gas displacing oil downdip that might result from gas-cap expansion or gas injection. This report outlines the calculation method used

By H. W. R. Wardlaw

This paper shows how the cost of the drilling operation per 100-ft drilled can be related directly to bit weight, rotary speed and hydraulic horsepower variations. This enables optimum drilling conditions to be determined in a rapid and concise manner. The paper also shows how data for such calculations cun be obtained quickly and accurately by the "drill-off" method suggested by Lubinski. Examples are given to illustrate the principles and methods used. INTRODUCTION The problem of drilling oil wells at lowest cost and hazard needs little introduction, since it has been with us virtually from the earliest days of our industry. With better understanding of the effect of the various factors of bit weight, rotary speed and hydraulic jet action, however, the problem has recently taken on a more scientific aspect, and engineers or groups of engineers have been assigned to study many of the facets of this problem not immediately apparent to the general observer. Moreover, the availability of better equipment and materials has enabled certain practices to be adopted which were hitherto impractical or uneconomic. This has paved the way for further efforts to reduce footage costs, and so on. The improvements to be made in drilling technology and cost reduction are by no means exhausted, and we still frequently hear of new records being set and new materials and equipment becoming available. Some of the analytical methods used by the author to assist in promoting faster and more economic methods of drilling, particularly in Venezuela, form the basis for this descriptive study. STATEMENT OF THEORY AND DEFINITIONS In a highly competitive industry, it is no longer permissible to regard penetration rate as the sole criterion for drilling efficiency. Every new technique, however spectacular, must nowadays be viewed in the cold light of whether it will eventually lead to lower cost per foot drilled. In the exposition which follows, therefore, actual drilling rates will play only a secondary role, and attention will be focused on the more important aspect of cost per foot drilled. Changes in drilling factors such as weight and rotary speed will be studied in the light of how they directly affect the cost of operations. The advantage of this method is twofold—(1) cost increments are indicated directly; and (2) the relationships are, surprisingly enough, simpler than those relating to drilling rates. To achieve these advantages, it is necessary to introduce a new unit of measurement which, although unfamiliar, is readily understood. The unit in question will be referred to as the "actual drilling cost factor" and will be designated by the symbol h. The unit is defined as the number of rotating hours taken to penetrate 100 ft of formation. The unit can be obtained immediately from the actual drilling rate (r) in feet per hour; thus, This latter transformation is so easily effected that it may be convenient to use h as an international symbol for drilling cost, and to delay conversion to dollars, bolivars, pounds, etc, per foot until actually required (e.g., for fiscal reasons, etc.). In actual operations, h was found to vary in a relatively simple manner with increases in bit weight. The variation in h with rotary speed was found to be somewhat more complex. However, this also could be enormously simplified by introducing a further factor known as L or relative bit life. The factor L is also defined very simply by dividing 1,000 by the rotary speed (N) in revolutions per minute.

By M. Kaitz

A continuing discussion in both the petroleum engineering and economic literature is directed to the difficulties encountered in the use of discounted cash flow rate of return (DCF) as a measure of investment worth. Although useful in most Instances, DCF has been criticized because it is time-consuming in its trial-and-error solution, theoretically invalid, not adaptable to cash flow streams yielding multiple return solutions and not entirely reliable in selecting between mutually exclusive investments. The theoretical invalidity of DCF stems from the reinvestment assumption implicit in the calculation that earnings are reinvested at the DCF rate. The emerging consensus of the economic literature is that the net present value or the present worth of the net cash flow stream (discounting at the average opportunity or cost of capital rate) is more correct and reliable. Other criteria proposed have been ratios of net present value divided by initial investment or by present value of all investments in a project. All of these criteria are simple to determine and explicitly assume a reinvestment rate for [he income generated by a project. This paper develops and discusses a theoretically valid profitability criterion which is simple to compute and retains the appeal of a percent return on investment. It is called "percentage gain on investment" or PGI. It measures the gain an investment is expected to realize over like capital invested in the average opportunity and explicitly considers reinvestment potential. Why add another Concept to the large array of investment criteria now available, any one of which, or perhaps a combination of several, appears to embrace the form's (or individuals) objectives? The answer is that not one of the existing criteria provides both a readily comprehensible and theoretically valid measure of risk coverage that has general application. The proposed PGI does fulfill these requirements. INTRODUCTION An ancient expression warns that "one must yield to the times" — there are better ways of doing things. A review of the petroleum engineering and economic literature on one topic alone, measurement of investment worth, certainly is witness to this truth. In use for a number of years, the DCF has recently received attention, directed mainly to its theoretical invalidity. Several alternatives to DCF have been proposed to provide a valid, simply determined criterion to describe investment worth and to overcome the criticisms previously mentioned. This paper introduces another method called percentage gain on investment (PGI) and is proposed as but one of several yardsticks that should be used in making investment decisions. MEASURES OF INVESTMENT WORTH This paper will consider only those criteria which give weight to the time pattern of future earnings. These criteria are usually compared with an average opportunity rate or cost of capital of the firm to judge the relative worth of the investment. For purposes of demonstration, a 9 percent average opportunity rate will be used throughout this paper. Implicit in the DCF calculation is the assumption that earnings are reinvested at the DCF rate. Some argue, though, that there is no reinvestment assumption,' that the DCF rate is simply that maximum rate of interest one can pay on the investment over the life of the project and break even. The determination of DCF is accomplished by discounting the net cash flow stream at that rate (DCF) which will yield zero. The question is: why should reinvestment potential be explicitly considered in calculating return on investment or other criteria measuring economic worth? Perhaps the answer lies in consistent or equal treatment of future cash flow. It appears entirely illogical to give different present worth value to $1 received, say, 10 years from now, which is the circumstance resulting in comparing projects with different DCF returns. In fact, $1 received in 10 years has the same value regardless of which project generated the income. DCF return thus favors investment projects which are expected to provide early income as compared to those providing long-term income. Not surprisingly, the controversy on reinvestment assumption is an old one. Hoskold' discussed this same problem in 1877. He considered the future income from mineral properties as an annuity or a series of fixed future payments. (As will be demonstrated, his equation can be modified for variable income.) Prior to Hoskold, the value or what one could pay for the mine, was determined with the use of standard, single-interest tables. Here is what Hoskold said with regard to these tables. "This table, and others of its kind, to be found in most works on annuities, is constructed correctly according to

By P. B. Crawford, W. L. Huskey

A study has been made to determine the effects of random vertical fractures existing in a reservoir matrix on the effective matrix properties and producing characteristics of a well. The study indicated that short vertical fractures may distort the isopotentials by 1 percent in the vicinity of the fracture. A correlation was developed which showed the effect of fracture density on the producing capacity of a well. The correlation indicated that fracture shape had little effect on producing capacity, but that total fracture length was closely correlated with the producing capacity. Effective values of reservoir permeability were found to correlate closely with fracture density when measured parallel to the direction of flow. INTRODUCTION In reviewing the petroleum literature, it is found that in most analytical solutions of hydrocarbon reservoir behavior the assumption is made that the rock is uniform throughout although core analysis and well logging surveys indicate that substantially all reservoirs are heterogeneous and some are fractured. Only in recent years have investigators studied the effects of reservoir hetrogeneity. Russell1 discussed the effect of fractures in reservoirs. The literature indicates that fractures greatly affect the exploitation of man reservoirs. This is described by Elkins and Skov2 and Little-field et al. 3 Several types of reservoir heterogeneity have been studied by unsteady-state methods. In 1959, Landrum et al. 4 showed that transient phenomena in hydrocarbon reservoirs could be studied by thermal models. pickering5 studied the flow of heat in stratified linear reservoirs consisting of plates of different metals separated by sheets of rubber which were suddenly exposed to a low-temperature source at one end. Cotman6 studied both laterally and vertically heterogeneous thermal model reservoirs with special reference to cross flow. Another type of reservoir heterogeneity involves the presence of highly permeable, circular lenses around or near the wellbore which was studied by Miesch7 using a thermal and a steady-state model. Many studies have been made of the effects of induced hydraulic fractures on well performance. 8-10 Work regarding vertical fractures that exist initially in the reservoir was reported by Givens, et al.11 Since there appears to be little published work on the importance and influence on reservoir performance of fractures that exist out in the rock matrix, this study was conducted to determine the effects of vertical fractures on pressure distribution, producing capacity and the effective permeability of the reservoir rock. The analogy between the flow of electrical current through a conducting body and the flow of fluid in a porous medium provides a convenient method of studying various aspects of petroleum reservoirs. For electrical radial steady-state flow,12 For steady-state fluid flow, With the above analogy, it is only necessary to establish the ratio of the values of the electrical units in terms of the corresponding flow units to make the analogy quantitative. The ratios are shown below.