Search Documents

By M. Muskat, M. O. Taylor

Theoretical calculations have been made on the performance to be expected of gas-drive reservoirs for various characteristics of the oil and gas and the producing rock. Such performante has been expressed graphically as curves for the reservoir pressure and gas-oil ratios of the production as functions of the cumulative oil recovery. The latter has been expressed in terms of percentage of pore space of the reservoir rock. Such curves automatically give values for the ultimate physical recovery, if the latter is interpreted as the recovery obtained when the reservoir pressure has declined to atmospheric, or to any other pressure chosen as defining the state of practical complete depletion. Calculations of these pressure and gas-oil ratio histories have been made for conditions in which the oil viscosity, the gas solubility and shrinkage, the size of an overlying gas cap, if present, the permeability-saturation characteristics of the rock, and the amount of connate water, have been individually varied, The results serve to show the extent to which the ultimate recoveries are sensitive to the important physical parameters characterizing the oil reservoir. As is to be expected, the ultimate recoveries are found to decrease with increasing oil viscosity. Because of the predominant effect of the oil shrinkage associated with the liberation of the gas in solution, the ultimate recovery will decrease with increasing gas solubility. Increasing gas-cap volumes lead to higher recoveries, although the contribution made by the gas-cap gas is small as compared with the oil expulsion by the equivalent amount of solution gas. As a whole, the oil recovery when expressed in percentage of the pore space is not very sensitive to the details of the permeability-saturation relationship. However, if the rocks possess an equilibrium free-gas saturation, the rise in gas-oil ratio will be retarded, and the 'ultimate recovery somewhat increased' The oil shrinkage associated with the liberation of the dissolved gas also leads to the result that as long as the connate water is immobile and the permeability ratio curve for the rock is a function only of the total liquid saturation, the stock-tank oil recovery will be less for a sand containing no connate water than for One with an Original water content as high as 3' per cent. The space voidage in the former case will be somewhat greater, but the effect of oil shrinkage will lead to smaller values for the equivalent stock-tank recovery. These calculations also give data showing how the productivity indexes for the producing wells will vary during the reservoir history. Because of decreasing Permeability to the oil and increasing oil viscosity, the productivity index will fall continuously as Production proceeds, and may finally reach values as low as 10 per cent of the initial productivity index. Introduction In a previous paper1 was developed the basic theory for the prediction of the production histories of gas-drive reservoirs. The data required for the application of this theory included the characteristics of the petroleum fluids, such as the viscosity of the Oil and gas phases, the solubility of the gas in the oil, and the oil shrinkage— all expressed as functions of the reservoir pressure—and the permeability-saturation characteristics of the producing rock, as

Jan 1, 1946

By T. Vasilos, J. T. Smith

THE mechanism of pressure sintering, or hot pressing, for ceramic materials, has been investigated by several researchers.1-8 Plastic flow has been suggested as the rate-determining mech-anism1,2 while grain boundary sliding3,4 and fragmentation processess have been found controlling in other studies. Vasilos and Spriggs5,6 have reported that pressure sintering occurs by a diffusion-controlled mechanism. For the pressure sintering of metals, investigations with lead7 and with tin6 have not identified the operative mechanism during consolidation of these powders. This note reports the results of a pressure sintering study initiated to determine the mechanism controlling the densification of spherical copper powder. Temperatures of 300°, 400°, 500°, and 600°C at 2500, 5000, 7500, and 10,000 psi pressures were investigated. The copper powder for this study was obtained from Metals Disintegrating Co., Elizabeth, N.J. The minimum purity was reported to be 99.9 pct. The powder size was -325 mesh and a lineal analysis showed the average size to be less than 1 p. An alumina die body and punches were employed. Graphite spacers were placed between the copper powder and the punches before pressing. It was expected that the graphite would form carbon monoxide with any entrapped air and form a reducing atmosphere within the die cavity. This condition was generally realized during pressing although a minor amount of oxidation occurred at 500" and 600°C. The oxide particles were located randomly in the pressed samples and were not considered to influence the densification process. Calibration tests showed the temperature of the copper pellet to be about 20°C less than the die wall temperature, probably resulting from heat losses through the punches. This difference was compensated during the experimental pressing operations. The hot pressing equipment employed for this work has been described by Spriggs et a1.9 and was calibrated with a Baldwin-Lima-Hamilton load cell. Travel of the punches was monitored with a dial indicator showing increments of travel to 0.0005 in. Relative density was calculated at selected time intervals for a known copper weight, punch travel, and die diameter. The accuracy of this calculation was *1 pct. Average grain size after hot pressing was obtained by lineal analysis. Typical densification curves are presented in Fig. 1 and Fig. 2 where relative density is plotted vs log time as a function of applied pressure for the 300" and 600°C pressure sintering temperatures, respectively. The curves for 400" and 500°C were similar to those at 300°C and consequently are not shown. The curves at 300°C as well as 400" and 500°C do not show an end-point or limiting density for a particular applied pressure. At 600°C, a density in excess of 95 pct of theoretical has been reached such that further densification is extremely slow for the 5000, 7500, and 10,000 psi pressures. It is suggested that further densification at these pressure levels may be retarded or prevented by discontinuous grain growth, by entrapped gases, and by other microstructural limitations as described by Cob1e.10 The curve for 2500 psi hot pressing pressure is a straight line over the time interval although a slight change in slope may be occurring at the last time intervals. The plastic flow models,1,2 derived for hot pressing, predict the existence of an end-point (or limiting) density dependent on the yield point of the material at the temperature of interest. The yield stress of copper at a 0.005 strain level has been

Jan 1, 1965

By J. S. Aronofsky

The results of polentiometric model studies and numerical computations are described. The purpose of these studies was to determine the influence of the mobility ratio on flooding efficiencies during water encroachment in petroleum reservoirs, in systems of regular well geometry. The direct line drive well pattern having relative spacing distances of one and one-half to one was investigated for water to oil mobility ratios of 10. 1. and 0.1. The results of these studies indicate that the pattern Sweep efficiency i.; very dependent upon the mobility ratio. INTRODUCTION The almost universal occurrence of water in the immediate neighborhood of oil-hearing sands lends considerable emphasis to the desirability of having a rational analysis of water encroachment into petroleum reservoirs. This problem has 1)become increasingly important with the accelerated use of water injection for secondary recovery operations. The water encroachment problem has been studied from both the theoretical and experimental viewpoints by Muskat.' He has formulated, in a precise manner, differential equations to express flow in a water encroachment system, and rigorous analytical solutions have been developed for the linear, radial, and spherical cases.' However, it becomes exceedingly difficult to solve the encroachment problem rigorously for any general two-dimensional system in which the shape of the two-fluid interface is not immediately evident from the geometry of the system. The difficulty lies in the fart that a rigorous analytical solution requires that the shape of the two-fluid interface be known simultaneously with its instantaneous posititon and with the pressure distribution on both sides of the interface. Although this paper concerns the study of the two-fluid system of oil and water. the methods described in the paper are generally applicable to other systems of incompressible and immiscible fluids in which the displacing fluid sweep-the displaced fluid down to a residual value. In each region it is assumed there is only one mobile fluid and its mobility is constant throughout the region. The Same procedure. with minor alteration, is applicable for gas injection systems provided the assumption can be made that both the displaced and displacing gases can be described by steady state formulas. For the particular two-fluid problem considered in this paper, there is a relative permeability. K, and a viscosity. p. associated with the water zone, and a different permeability and viscosity associated with the oil zone. The K/p (mobility value associated with a particular volume element of the porous medium. at any instant of time. will depend upon whether the element is ahead of. or behind the advancing interface. It is assumed that all volume elements ahead of the interface are characterized by the K/p value for oil, while all the elements behind the interface are characterized by the K/p value for water. That is. the K/p value in the oil zone is constant throughout the oil zone. and the K/p value in the water zone is also constant but different from the value in the oil zone. The symbols are defined at the end of the paper. The well pattern chosen for initial study is the direct line drive; the swept-out area at breakthrough was determined by two different methods which are described in the section on theory. These are: 1. Potentiometric analyzer 2. Numerical computation The use of the potentiometric analyzer for studying water flooding problems is adequately described in the literature, particularly for cases where the mobility in the water zone is assumed equal to the mobility in the oil zone." ' 5

Jan 1, 1952

By C. C. Miller, C. A. Hutchinson, A. B. Dyes

This paper presents a simple and practical method for the estimation of effective permeability and reservoir pressures from the rise in bottom hole pressure when a well is shut in. The equations on which the method is based consider a stabilized well producing from a uniform formation when flow is shut off at the sand face followed by the unsteady state flow of a compressible liquid from regions far from a well to points closer in during the build-up period. An approximation for extending the method for two-phase flow in the reservoir is also presented. In order to compare the effects of extremes in the conditions existing at the external boundary of the reservoir, two general equations defining the build-up characteristics of the reservoir are shown: first, for a reservoir with constant pressure maintained at the radius of drainage, and second, for a reservoir with no influx over the radius of drainage. These equations define the build-up characteristics for the shut-in well as functions of time, the permeability and porosity of the drainage area, the viscosity and compressibility of the fluid flowing, the production rate prior to shut-in, and the radius of drainage of the well. When a well is shut in at the surface, rather than at the sand face, flow will continue into the well bore after shut-in due to the compressibility of the fluids in the bore. The effect upon the build-up characteristics of this flow into the well after shut-in has been studied by means of the electrical analyzer. From the basic relationships and consideration of the limitations due to continued flow after shut-in, methods have been developed for the estimation of the effective permeability of an area from the plot of the build-up characteristics of the well producing that area. Likewise, equations are presented for the estimation of the reservoir pressure from an early portion of the build-up curve. By means of the electrical analyzer, studies have been made concerning the variations in the build-up characteristics resulting from damaged or improved permeability in the immediate vicinity of the well bore. These studies show that the permeability estimated from the build-up curves represents the average effective permeability of the majority of the drainage area and that the permeability in the immediate vicinity of the well does not affect this estimate. Since the major portion of the drawdown is in the immediate vicinity of the well, variations in the permeability of this area greatly affect the productivity of the well as expressed by the productivity index. From pressure build-up and productivity index tests the ratio of the permeability of some arbitrarily specified area immediately surrounding the well and the effective permeability in the majority of the drainage area can be computed. This ratio is a measure of any damage to the well incurred in drilling or reworking and likewise any benefits as a result of remedial measures. In build-up tests and many other tests that are made on a producing well it is at least desirable, if not absolutely necessary, that the well be stabilized prior to testing. A method is presented here for relating the time required to stabilize a well and the reservoir and fluid properties. An equation is also presented for approximation of the length of shut-in time required to obtain sufficient build-up data for the estimation of permeability and reservoir pressure.

Jan 1, 1950

By P. S. Rudman

Crystal-structure occurrences in binary phase diagrams of fifth- and sixth-period elements, zirconium through palladium and hafniun through platinum, are surveyed with the conclusion that average .group number GN is the dominant correlative. This is interpreted in terms of curves of enthalpy differences between competitiue structures that are continuous func -tions of GN. It is further reasoned that the absolute formation energies may also be continuous functions of GN and if so can be simply obtained from the regatives of the heats of sublimation of the elements. It is shown that from such a heat of formation-GN plot the heats of mixing of alloy systems can be quantitatively predicted. For example, for sixth-period alloys very large exothermic heats of mixing are predicted for systems that span GN = 6 (i.e., Hf-Pt, Ta-Os, ...). Thermodynamic data with which to test these predictions do not appear to be available. If in the latter transition (A subgroup) elements we eliminate from consideration the fourth period where magnetic terms contribute to crystal-structure stability, as shown in Table I, there is a very regular sequence of structures. Using the group number (GN) assignments of Table I, which differ from the usual assignments of 8 to the elements listed as 8 through 10, it is perhaps trivial to say that the phase occurrence is a unique function of GN. However, the interesting aspect of alloys from among these elements is that the alloy structures also appear to correlate with GN. The GN for an alloy is given as GN = SixiGNi, where xi and GNi = atomic fraction and GN of component i, respectively. It is difficult to credit anyone in particular with having discovered this correlation, so many people have implicitly or explicitly employed it in the last 10 years. However, Bloom and Grant' in their correlation of s-phase occurrence with total number of d + s electrons were certainly one of the first to appreciate the correlation and Hume -Rothery2 appears to be responsible for the use of average GN as defined above. Only in the last few years have the phase diagrams of enough alloy systems been determined for this correlation to be surveyed with any degree of completeness. We shall here, in an abbreviated format, summarize the available phase-diagram literature and discuss some consequences of the GN correlation. BINARY DIAGRAM SURVEY In Figs. 1 to 3 the phase occurrences are summarized. In Fig. 4 the occurrence limits for the more common structures are summarized. That the GN of an alloy falls within the occurrence limits of a particular structure is a necessary but not a sufficient condition for the appearance of that structure. The preceding figures certainly illustrate that GN is a dominant correlative for structures of these elements. Many of the diagrams have complications, with ordering and other phase changes occurring in intermediate phases, that cannot be clearly accommodated in such a schematic survey and have hence been left out. Such omissions cannot, however, negate the GN correlations discussed here. Relative atomic size is also an important crystal-structure determinant and since atomic size varies quite regularly with GN for these elements it is possible that the apparent GN correlation is basically a size correlation. In Fig. 5 relative size occurrence limits are given for intermediate phases and extended terminal solutions. The relative volume, ?v/(v) - 2(vA - vB)/(vA + uB), has been taken as the measure

Jan 1, 1968

By Richard H. Jahns

Gem materials, comprising those minerals and closely allied natural substances used for personal adornment, for the fashioning of ornamental objects, or for other decorative purposes, have been valued in human commerce for thousands of years. Amber, jet, natural glasses, and at least 15 mineral species were sought and used by primitive man, and most of the materials in common use today were known well before the time of Christ. Their variety is remarkably great, and the traceable relationships among their respective sources, modes of occurrence, properties, and histories of production, preparation, and valuation are extremely complex. Their uses have been affected by numerous shifts in fashion, and from time to time by superstition, by decrees of church and state, by personal interpretations of religious significance, and by various presumptions concerning medicinal value and influence on human behavior (Ball, 1935, 1950; Kunz, 1913; Pogue, 1915; Whitlock, 1946). Thanks to their durability, portability, and relative scarcity, some gem stones have been employed as media of exchange, and the investment value of the more precious stones has long been recognized and exploited, especially in Asiatic and European countries. Demands for gem materials have stimulated trade among various peoples for centuries, and they also were responsible for much early geographic exploration. During more recent times they have prompted numerous geological studies and a considerable amount of research on the growth and properties of crystals. Apart from decorative uses, some naturally occurring gem materials also are employed directly in numerous fields of industry. Substances that are manufactured for gem or industrial uses are termed synthetic gem materials; most of these correspond in composition and crystal structure to certain minerals found in nature. Imitation materials are those whose use is based upon their resemblance to substances of different composition and greater intrinsic value. Gem stone is a designation commonly restricted to those materials suitable for personal adornment, and the complementary term ornamental stone is applied to the remaining materials not used industrially. This distinction is arbitrary and somewhat inconsistent, owing to numerous geographic differences and historical shifts in end uses. Further, there is a complete gradation between the category of ornamental stones and that of dimension stone. Substances that lose their identity in use are excluded from the gem category; thus a transparent crystal of beryl would not be classified as a gem stone if it were processed as a source of beryllium metal. The adjectives precious and semiprecious have been employed as a means for distinguishing different gem materials on the basis of recognized value. Diamond, emerald, ruby, and sapphire consistently have been regarded as precious stones, and in cut form they commonly are referred to as the noble gems. Opal, pearl, turquoise, and chrysoberyl also have been classed as precious at various times. In recent years, however, it has been increasingly recognized that no satisfactory distinction can be based upon assessed value or any group of specific properties, and now all gem stones normally are referred to as precious. Indeed, the tendency is to extend this term to include many of the so-called ornamental stones, and the designation semiprecious has been essentially discontinued. Unlike most mineral commodities of economic value, gem materials are little altered in processing for end uses. Some stones are used in essentially their natural form, and most are reshaped and polished without significant changes in their fundamental properties. Thus the preparation of gem and ornamental stones is intended chiefly to emphasize or enhance

Jan 1, 1983

By Richard H. Jahns

Gem materials, comprising those minerals and closely allied natural substances used for personal adornment, for the fashioning of ornamental objects, or for other decorative purposes, have been valued in human commerce for thousands of years. Amber, jet, natural glasses, and at least 15 mineral species were sought and used by primitive man, and most of the materials in common use today were known well before the time of Christ. Their variety is remarkably great, and the traceable relationships among their respective sources, modes of occurrence, properties, and histories of production, preparation, and valuation are extremely complex. Their uses have been affected by numerous shifts in fashion, and from time to time by superstition, by decrees of church and state, by personal interpretations of religious significance, and by various presumptions concerning medicinal value and influence on human behavior (Ball, 1935, 1950; Kunz, 1913; Pogue, 1915, Whitlock, 1946). Thanks to their durability, portability, and relative scarcity, some gem stones have been employed as media of exchange, and the investment value of the more precious stones has long been recognized and exploited, especially in Asiatic and European countries. Demands for gem materials have stimulated trade among various peoples for centuries, and they also were responsible for much early geographic exploration. During more recent times they have prompted numerous geological studies and a considerable amount of research on the growth and properties of crystals. Apart from decorative uses, some naturally occurring gem materials also are employed directly in numerous fields of industry. Substances that are manufactured for gem or industrial uses are termed synthetic gem materials; most of these correspond in composition and crystal structure to certain minerals found in nature. Imitation materials are those whose use is based upon their resemblance to substances of different composition and greater intrinsic value. Gem stone is a designation commonly restricted to those materials suitable for personal adornment, and the complementary term ornamental stone is applied to the remaining materials not used industrially. This distinction is arbitrary and somewhat inconsistent, owing to numerous geographic differences and historical shifts in end uses. Further, there is a complete gradation between the category of ornamental stones and that of dimension stone. Substances that lose their identity in use are excluded from the gem category; thus a transparent crystal of beryl would not be classified as a gem stone if it were processed as a source of beryllium metal. The adjectives precious and semiprecious have been employed as a means for distinguishing different gem materials on the basis of recognized value. Diamond, emerald, ruby, and sapphire consistently have been regarded as precious stones, and in cut form they commonly are referred to as the noble gems. Opal, pearl, turquoise, and chrysoberyl also have been classed as precious at various times. In recent years, however, it has been increasingly recognized that no satisfactory distinction can be based upon assessed value or any group of specific properties, and now all gem stones normally are referred to as precious. Indeed, the tendency is to extend this term to include many of the so-called ornamental stones, and the designation semi- precious has been essentially discontinued. Unlike most mineral commodities of economic value, gem materials are little altered in processing for end uses. Some stones are used in essentially their natural form, and most are reshaped and polished without significant changes in their fundamental properties. Thus the preparation of gem and ornamental stones is intended chiefly to emphasize or enhance

Jan 1, 1975

By Herbert E. Hawkes

“For his Pioneering leadership in the science and technology of mineral exploration, especially in the development and world-wide application of geochemical methods as major exploration tools, and for his lecture ‘Geothermal Hydrogen’.”

Jan 1, 1981

By E. A. Cisney

The formula for calculating the crystallographic angles of a tetragonal lattice is: C°S =where $ is the angle between planes (HKL) and (hkl). The angles in Table I have been calculated for indium

Jan 1, 1960

By C. C. Templeton, S. S. Rushing

Methods previously developed for the study of air-liquid displacements in microscopic capillaries (inner diameters of 3 to 40 microns) have been used to investigate oil-water displacements in capillaries initially filled with water. Displacement calculations assuming perfect displacement and no capillary pressure hysteresis yielded oil effective viscosities smaller than the macroscopic viscosities. For a given liquid pair, the oil effective viscosity decreased both with decreasing capillary size and with increasing oil-water viscosity ratio. This behavior can be explained by the existence of an annular water film (20 A to 260 A thick) on the capillary wall. When the capillary was first filled with oil, the ratio of the oil effective viscosity to the normal oil vircosity was highest for the first water displacement and decreased with subsequent displacements. Sometimes the oil effective viscosity ratio during the initial water displacement was greater than unity. INTRODUCTION In a previous paper1 a technique was described for studying air-liquid and liquid-liquid displacements in very small capillaries of uniform diameter, in the hope that such microscopic data would further the understanding of the nature of multiphase fluid flow through porous media. That paper contained comprehensive data for air-liquid displacements in Pyrex capillaries, with a few data for oil-water displacements in capillaries initially filled with water. The purpose of this paper is to present more complete results for oil-water displacements in capillaries initially filled with water, and to describe for the first time such observations in capillaries initially filled with oil. In this way the effect of the wetting history of the system upon the displacement process may be studied. METHODS The basic techniques employed were described in the previous paper.1 The measuring procedure and the working equations will be briefly summarized here and a few modilications will be pointed out. For the present study, temperature control within 0.l °C was obtained by placing the microscope and its immediate accessories in a thermostated air bath made of a steel frame covered with plexiglas. The water and oils used in this work were the same as in the previous paper,1 with the addition of "Medium Mineral Oil, U.S.P. supplied by the Harshaw Scientific Co. When the first liquid was introduced into the horizontal capillary, the air-liquid static capillary pressure Poc was measured along the observable length of the capillary. From the relation Poc = 4y cos 0/d, the capillary diameter d could be calculated from Poc if the "microscopic" air-liquid boundary tension, To = y cos On, was taken as equal to its known "macroscopic" value. This calculation involves assuming that cos 6 = 1 and that capillary size or the resulting high interfa-cial curvature has no effect on surface tension, or y cos 6, as was verified by our work, and also by that of Cohan and Meyers on air-liquid-solid systems. For several capillaries these calculated diameter values were compared with values measured visually with a filar micrometer eyepiece. For a given capillary the average values for the two methods agreed within 1 or 2 per cent, but the visual values were less reproducible than the calculated values. Further the calculated values were much less subject to personal bias and took into account any slight ellipticity that might be present. Hence diameters calculated from air-liquid Poc data were used throughout this paper. The second liquid was introduced in such a way that a single spherical interface separated the two liquids. Displacements were always made with a constant pressure, Pt, between the ends of the capillary. The times (t) at which the interface passed certain

Jan 1, 1957

By David Cornell

Back pressure test data on natural gas wells are short time test data of unsteady state nature. Performance curves from which unsteady state effects have been eliminated, called stabilized performance curves, are necessary for understanding the behavior of gas wells and for scheduling gas fields for periods up to 20 years. A calculation procedure has been developed whereby three day, seven day, 30 day, or other duration stabilized performance curves may be calculated from back pressure test data. Such curves are essential for correct gas deliverability calculations and gas field scheduling. I NTRODUCTION The manner in which a back pressure test is made will affect the slope of the back pressure curve because the radius of drainage changes with the length of time that the well is allowed to flow. The reciprocal of the slope of the back pressure curve is defined as n in Equation 1 and varies from 0.5 to 1.0. If the flow is viscous and the radius of drainage constant, n will have a value of 1.0. Because of turbulence in the formation near the wellbore, the pressure drop is usually greater than the viscous flow pressure drop and the exponent n is frequently less than 1.0. In addition, since the radius of drainage moves away from the well-bore during the test, the gas flows a greater and greater distance as the test proceeds, with a resulting increased pressure drop at the high flow rates, and, again, n for the test tends to be less than 1.O. The effect of the changes in the radius of drainage may be calculated and the results used in preparing stabilized performance curves for fixed values of the radius of drainage. Such curves can be made for a defi- nite time period, such as three day, seven day, or 30 day stabilized performance curves. These cukes will show the effect of turbulence near the wellbore and can be used to predict the performance of the gas well at any flow rate or for any production period. In order to discuss the importance of the radius of drainage quantitatively, it should be defined. Muskat" defines the radius of drainage as follows: "... the radius of drainage gives the distance from the well to which the approximately steady state condition has been established ... it being assumed that no gas is removed from any point until the radius of drainage has passed that uoint." An eauivalent definition of the radius of drainage is that radius at which the extrapolated portion of the steady state portion of the pressure gradient gives a pressure equal to the shut-in reservoir pressure. This definition is also implied in the back pressure equation. The radius of drainage is, then, the distance into the reservoir from which the gas appears to originate at any time. It is, however, only a useful concept and does not actually exist. Fig. 1 illustrates the definition of the radius of drainage. At the end of each flow period the radius of drainage may depend on three factors: (1) the duration of each flow period; (2) the rate of gas production; and (3) the rate and duration of the previous flow periods. In general, back pressure tests are conducted with various lengths of production time and at various flow rates. Each individual back pressure test, therefore, requires an individual analysis. The calculations of the change in the radius of drainage have been made, however, for the case of equal flow periods for all four production rates and several sets of increasing flow rates. Fig. 2 illustrates the results of such calculations made by the method of Cornell and Katz2 for the case of 8, = 50,000 and the flow rate increasing in the ratio 1.0:2.0:3.0:4.0 for each flow period. The increase in the radius of drainage during the test is apparent.

Jan 1, 1956

By J. E. Miller, W. E. DeVaney, L. Stroud

During a recent phase study of a natural gas, two stable equilibrium liquid phases were observed at temperatures below —200°F and pressures above 200 psi. This paper reviews the published literature on the occurrence of multiple equilibrium liquid phases and presents analytical data for the vapor and two equilibrium liquid phases of the liquefied natural gas at five experimental conditions. In addition, data for 30 conditions of two-phase equilibria are included. INTRODUCTION The low-temperature phase behavior of gases, most of which contained helium, has been investigated in the laboratories of the Helium Activity for many years. Since 1952, experimental studies of these systems have been continuous as part of the research program at Amarillo, Tex. Because of their value to private industries interested in participating in the Helium Conservation Program, several "Open File" reports containing phase equilibria data for heliurn-bearing natural gases have already been released by the Helium Activity. A paper containing information on the general phase behavior, operating criteria and extensive vapor-liquid data for two helium-containing systems was recently published.2 Additional publications presenting experimental data on the phase relationships of various gas systems are now in process and will be available in the near future. PREVIOUS EXPERIMENTAL WORK Although the formation of multiple liquids has been reported for various systems, to our knowledge this paper is the only substantiated evidence of a vapor-liquid-liquid equilibria in a naturally occurring gas. In 1940, Vink, Ames, and others13 reported the presence of two liquid phases in a hydrocarbon system consisting of mixtures of crude oils, solvents and natural gas. Eilerts and co-workers; published data on the recombined fluids from a gas-condensate well. This condensed gas, containing approximately 76 per cent methane and 24 per cent ethane-plus, exhibited two distinct liquid phases. Weinaug and Bradley 14 observed "unusual" phase behavior in a reservoir mixture. These workers postulated that the anomalous phase behavior was due to the "imminent formation of a second liquid phase". Botkin, Reamer, Sage and Lacey 1 studied two California crude oils that exhibited multiple phases. Kohn and Kuratal0 recently reported two equilibrium liquid phases in the methane-hydrogen sulfide system. Roof and Crawfordl1 and Eakin, et a16 also have reported experiments with binary systems that formed two stable equilibrium liquid phases. APPARATUS AND PROCEDURE A U. S. Bureau of Mines Phase Equilibrium Apparatus was used in conducting this study. The apparatus and procedures employed in its operation have been previously described3 and will not be repeated in detail ill this report. Briefly, the apparatus consists of a windowed cell which can be maintained within ±0.5°F for temperatures between room temperature arid -320°F. Pressure within the cell can be maintained within 0.1 per cent of gauge reading up to 800 psig. Equilibrium vapor and liquid samples are obtained in special containers4 for analysis by a mass spectrometer. Although the accuracy of the analyzer is about ±0.1 mol per cent, the reproducibility of the phase-equilibrium apparatus is considered relatively poor. Values reported from methane and nitrogen are considered accurate to ±l.0 and ±0.6 mol per cent. Data for ethane-plus in the vapor are accurate within 0.2 mol per cent; liquid-phase data for this aggregate compolnent are accurate within 1.5 mol per cent. For helium in the vapor phase, the analytical data are accurate within 0.2 mol per cent; liquid-phase analyses for this component were obtained by the charcoal adsorption method described by Frost8 and are accurate within 0.006 mol per cent. All of these references to the accuracy of reported values are conservative estimates based upon a statistical treatment of reproducibility data obtained with the apparatus.

By Robert G. Nisle

In conducting a pressure build-up test on an oil well, it is often necessary to shut-in the well for a short time prior to initiation of the test. The effect of such a short tern shut-in on the resulting pressure build-up curve is derived. The derivation is based on the Kelvin point-source solution and is a straight forward extension of the fundarnental theory. The magnitude of the effect is evaluatecl and illustrated by an example. It is shown that if the production subsequent to the short term shut-in is at least 10 times the duration of the shrct-in, the effect will he less than 10 per cent. INTRODUCTION It often happens that before a pressure build-up test can be made an additional valve must be mounted on the wellhead, or on arrival at the well, it might be learned that the well had been shut-in several hours the previous day because of a shortage of tankage. In such cases the question is How will such an interruption of production affect the results of a pressure buildup test This work explains how to take account of such a short term shut-in in order that the test can proceed without further delay. The solution presented here makes use of the Kelvin point-source solution used by Horne1 in his treatment of the pressure build-up in wells. Use is also made of the super-position theorem. MATHEMATICAL DERIVATION The well is assumed to have been in production at a constant rate for a relatively long time (several days or more). It is then shut-in for a short time, after which it is again placed in production at the previous constant rate. It continues to produce at this rate until shut-in for a pressure build-up test. The Kelvin point-source solution is expressed by Eq. 1. The sequence of production rates is listed in Table I. TABLE 1 — SEQUENCE OF PRODUCTION RATES Time interval Production Rate At time AT the pressure will be the resultant of the sequence of production rates shown in Table 1. An expression for the super-posed effects may be derived from this table. (Effect of + q for T + ? T) + (Effect of - q for T + ?T-t,) + (Effect of + q for T + ? T - t,) + (Effect of - q for ? T) = Resultant........(2) The desired result is obtained by substituting the effect defined by Eq. 1 for the terms in Eq. 2. The result is Eq. 3. When the logarithmic approximation is made for the exponential integrals and the terms collected the result is Eq. 4. Let T — t2 = T = flowing time following short term shut-in t2 — t1 = 6 = duration of short term shut-in Then T - t1 = r + 6 Substitution of these quantities in Eq. 4 gives p(rw ?T) = P,, -When d<<t - AT then d—&apos;-, T— approaches 1. r + AT In this case the effect of the short term shut-in may be ignored. In order to get some idea of the numerical magnitudes involved take the worst possible case, when AT = 0, that is at shut-in. If the error is made less than some maximum acceptable value at this time, it will be less than the maximum for all subsequent times. d + r Let the difference between the value of y-— and 1 be less than 0.1. Then —6+r — 1 < 0.1 S < O.l T Therefore, if the flowing time following a short term shut-in is at least 10 times the duration of the short term shut-in, the error introduced by ignoring the short term shut-in will not exceed 10 per cent of T+AT A T

Jan 1, 1957

By R. Schuhmann, E. J. Michal

Experimental measurements are reported for the oxygen pressures of iron-silicate slags in equilibrium with solid silica. CO2-CO mixtures were bubbled through the slags in silica crucibles to find equilibrium values of p co2/p co over the entire composition range of silica-saturated slags at 1250°, 1300°, and 1350°C. From the data, the activity vs. composition relationships were calculated for O, FeO, and Fe. SLAGS composed primarily of silica and iron oxides have particular practical importance in copper smelting,&apos; in acid steelmaking, and in accounting for the service behavior of silica refractories. These slags all contain both ferrous and ferric iron, so that, strictly speaking, they are ternary or three-component liquid solutions. For convenience, the components may be taken as the three oxides, FeO, Fe2O3, and SiO2. The constitution of the FeO-Fe,O,-SiO, system at slagmaking temperatures has been studied in a recent investigation described in a separate paper.&apos; This work showed the compositional ranges of iron-silicate melts at temperatures up to 1450°C and also the temperature-composition relationships for melts in equilibrium with various solid phases, including iron, wüstite, magnetite, and silica. Broadly speaking, the compositions of iron-silicate melts at a given temperature should be thought of as varying over a two-dimensional field. In the first place, the silica/total iron oxide relationship varies from maximum silica for silica-saturated melts to minimum silica for wüstite-saturated melts or, at higher temperatures, to zero silica for straight iron oxide melts. In copper smelting, for example, the converter slags approach the low silica limits while the reverberatory-furnace slags have intermediate silica contents, perhaps approaching silica saturation in some circumstances. Liquids formed in the pores of silica refractories presumably are close to silica saturation. The second type of composition variation is in the ferric/ferrous relationship, which depends on how oxidizing or how reducing the process conditions are. Relatively reducing slags, such as acid-steelmaking slags, are lower in ferric oxide and ap- proach saturation with metallic iron. Oxidized slags, such as copper-smelting slags, contain larger percentages of ferric oxides and approach or reach saturation with magnetite. Since these two kinds of compositional variations are closely related to two general properties of interest for all metallurgical slags, basicity, and oxidizing power, it is hoped that the thermodynamic study of iron-silicate slags eventually will throw further light on the significance and nature of these long-recognized but as yet poorly understood slag properties. Measurements of oxygen activity as a function of slag composition, using CO-CO2 gas mixture to fix oxygen activities, afford an experimentally convenient approach to the study of the thermodynamic properties of iron-silicate slags. In a preceding investigation,³ his approach was used in a systematic study of iron-silicate slags in equilibrium with solid iron. The present paper gives the results of a similar study of the iron-silicate slags in equilibrium with solid silica. The experimental work consisted in the determination of CO2: CO ratios in gas mixtures at equilibrium in the system: gas: liquid slag: silica Determinations were made at temperatures from 1250° to 1350°C for slags ranging in composition from iron saturation to magnetite saturation. Also, equilibrium CO,: CO ratios and slag compositions were determined for the univariant systems: gas: liquid slag: silica: iron gas: liquid slag: silica: magnetite Experimental Methods The experimental procedure was essentially the same as that used in the previous study of iron-silicate slags saturated with solid iron," but in order to work with the silica-saturated slags an all-silica crucible assembly was used. A liquid slag was prepared from pure prefused stock materials and brought to equilibrium with the silica crucible at the desired temperature. Then a mixture of CO2 and CO was continuously bubbled through the melt

Jan 1, 1953

By R. V. Smith

The theory for unsteady radial flow of gas, as simplified by Aronofsky and Jenkins, has been reviewed and compared with actual well performance. This study indicated that the theory should be modified by the introduction of an empirical "rate of flow" function "Y(q)". The expansion of the theory on the flow of gas to include a Y(q) term bridges the gap between the theory of unsteady-state flow of gas and actual gas-well behavior. Apparently, the Y(q) term is a function only of the rate of flow for a given well. The completion factor or skin effect is associated with the Y(q) function in such a manner that at least two or more sets of drawdown or pressure build-up tests are needed to separate the Y(q) function and the completion factor. Since the Y(q) function used in this report represents an energy loss dependent on rate of flow near the well-bore of a gas well which is in addition to the loss required by Darcy&apos;s law, the Y(q) function is related to the exponent of the back-pressure curve for a gas well. Techniques and equations are presented which permit the estimation of stabilized and "short-time" deliver-abilities. INTRODUCTION The published theory of radial unsteady-state flow of gas through porous media may be divided into papers1-" which use Darcy&apos;s law as a premise and other papers4-7,16-10 which recognize that Darcy&apos;s law may not hold for flow into gas wells. Miller8 compared theory based on Darcy&apos;s law with actual gas-well behavior and concluded that flow into gas wells deviates from Darcy&apos;s law. Houpeurt4,5 in a thorough investigation of flow into gas wells concluded that the flow rate into a gas well is proportional to the difference between the squares of upstream and downstream pressures elevated to a power between 1.0 and 0.5. Houpeurt also concluded that the deviation from Darcy&apos;s law is caused by irreversible, kinetic energy exchange between the flowing gas and the porous media. Tek9 presented a complete evaluation of the flow problem in gas wells. Several investigators state that the inclusion of gas compressibility into flow equations based on Darcy&apos;s law does not begin to explain the flow problems suggested by gas-well behavior. For many years, natural-gas engineers have realized that gas wells by actual test exhibit behavior which cannot be explained entirely by Darcy&apos;s law. Many tests on gas wells have indicated that a value of 0.85 for the exponent in the back-pressure equation is a fair average for many wells. However, Darcy&apos;s law either for steady or unsteady-state conditions indicates that the exponent should be 1.0. The purpose of this paper has been to examine this deviation from Darcy&apos;s law and to present a reasonable mathematical means of evaluating the deviation. The method presented in this paper combines previous theory with an empirical function to provide a description of gas-well behavior. All equations in this report are for constant terminal-rate-of-flow conditions. While present methods of testing wells usually are not based on this concept, they are believed to approach the constant terminal-rate conditions closely in many tests. However, several gas wells were tested by holding the rates of flow constant in order to approach more closely the constant terminal rate of flow; the results confirm the concepts presented in this paper. REVIEW OF THEORY FOR UNSTEADY-STATE FLOW OF GAS The basic partial differential equation for the flow of natural gas has been given in technical literature; the equation for a cylindrical reservoir with the well in the center is

By T. D. Mueller, F. C. Miller, Jr. Earlougher R. C.

In analyzing pressure buildup tests for field wells producing both oil and gas, the common practice is to use a modification of single-phase flow theory. Validity of such an approximation has been demonstrated for single-well solution-gas-drive systems. This paper indicates that such approximations are also valid for two-well solution-gas&apos;-drive systems. From this it is inferred that the technique can be used for multiple-well systems. A computer was used to simulate the behavior of a two-well solution-gas-drive reservoir to test the validity of the above type of analysis. Simulation results indicate that pressure buildup tests in such a system can be analyzed within engineering accuracy for formation permeability and pressure. A rule of thumb is given for estimating the length of time a well must be shut in for the pressure in a nearby producing well to increase significantly. To observe such an increase, the shut-in well would have to be left shut in much longer than normal. INTRODUCTION An important technique for obtaining data concerning a producing petroleum reservoir is the pressure buildup test. Such a test, when properly conducted and analyzed, provides information on the average reservoir pressure and the permeability in the major drainage area of a well. The theory upon which the analysis of a pressure buildup test is based assumes that the behavior of the fluid in the reservoir is adequately described by the diffusivity equation.1 Use of the diffusivity equation implies that a single fluid of small and constant compressibility is flowing. To analyze a pressure buildup test in situations which involve more than one fluid phase, the single-fluid analysis is extended. This paper verifies the extension of pressure buildup analyses to two-phase, two-well systems. Miller, Dyes and Hutchinson1 have presented a technique for analyzing pressure buildup tests in circular bounded reservoirs. Their method requires that the reservoir be producing at pseudo-steady state (constant pressure gradients) prior to shut-in. An alternate reservoir model, presented by Horner,2 assumes that the reservoir is infinite. The Homer model does not assume a Pseudo-steady state prior to shut-in, but the assumption that the reservoir is infinite implies that the average pressure can be determined only if the total production prior to shut-in is small compared to the total fluid originally present. Limitations imposed by the pseudo-steady state and infinite reservoir assumptions are avoided by the technique proposed by Matthews, Brons and Hazebroek.3 Their method can be used to calculate both the Permeability and the average pressure in bounded single-well systems which are not necessarily at steady state. They also suggested determining the average pressure of a multi-well reservoir producing from pseudo-steady state by volumetrically averaging the static pressures of the individual wells. Matthews and Lefkovits4 used numerical simulation techniques to verify that this method does produce adequate results for single-phase, multiple-well systems. Perrine5 proposed modifications to the single-phase theory so that pressure buildup tests in multiple-phase systems could be analyzed. He suggested replacing the single-phase mobility (k/µ) by the sum of the mobilities of the individual phases, and replacing the fluid compressibility by an average compressibility weighted by the saturations of the separate phases. By using numerical techniques to solve the two-phase flow problem for a single-well system, he concluded that this approximation gives results within engineering accuracy. More recently, Weller6 has shown that, for a gas-oil system with a single well at the center of a circular reservoir, pressure buildup tests can be analyzed adequately by using the modification proposed by Perrine. Feller&apos;s results also indicate that, as the gas saturation increases, this analysis becomes less accurate. To date, no evidence is available to show that the Perrine approximations can be combined with the Matthews-Brons-Hazebroek technique to calculate the average pressure of a multi-well,

By I. S. Levy, J. L. Brimhall, B. Mastel

Specimens of&apos; Hastelloy X-280, a low-cobalt version of the solid-solution- hardened nickel-base alloy Hastelloy X, were given a series of thermomechanical treatttzents. They were then Lensile-tested at room temperature and 1350°F. The Microstructures of the specimens tested at 1350°F were studied by electron Microscopy. Microscopy showed that dislocation tangling, caused by the cold working, mas still in evidence after aging-for 24 hr to temperatures as high as 1550oF. (Cr, Mo)23 C6, (Fe, Mo)6 C&apos;. as well us an unidentifiable phase precipitale due to the combination of the thermomechanical treatments and subsequent tensile strain. The room -temperature and 1350°F tensile properties show a general increase ill strength and decreccse in ductility with increusing cold-work level. This is due to the dislocation taangling-. PerLurbatiot~s to this general trend occur- due to solid-solution weakening at the higher cold-work levels as a direct result of the preciplitation reactions. THE requirements for corrosion resistance and high-temperature strength place nickel-base alloys among the more promising materials for advanced nuclear reactor applications. However, when these alloys are irradiated at service temperatures, they generally suffer severe embrittlement&apos;" and reduction in rupture life.3 Certain experimental treatments resulted in improvements in the tensile and stress-rupture properties of Hastelloy X-280, a solution-hardened nickel-base alloy.&apos; To provide a basis for understanding these improvements, unirradiated tensile specimens were examined by electron microscopy. The results are reported in this paper. EXPERIMENTAL DETAILS Material. Hastelloy X-280, a low-cobalt version of Hastelloy X, is a member of the solution-hardened class of nickel-base alloys. The Hastelloy X-280 was obtained from 4- and 1-in: thick plates in the mill -annealed condition. Table I gives the chemical composition of the material as reported by the National Spectrographic Laboratories. ~x~erimental Treatments. Plates of suitable size, 12 by 3 by 1 in., or by + in., were cut under a water jet from the stock material and examined ultrasonic -ally for homogeneity of grain size and absence of voids. The acceptable plates were solution-treated at 2150°Ffor 1 hr and air-cooled. These plates were then reduced in thickness by cold rolling 0, 10, 15: or 20 pct along the original rolling direction. The plates were machined into buttonhead tensile specimens with a gage diameter of 0.125 in., with a gage length of 1.125 in., and with their tensile axes parallel to the rolling direction. The specimens were examined for quality of finish with a penetrant dye, ultra son ically checked for voids, and finally given a dimensional tolerance check with an optical comparator. The acceptable specimens were aged for 24 hr at either 1000" or 1500 F and air-cooled. These specimens were then held at 540°F in a water loop for approximately 18 days to qualify as controls for similar specimens irradiated under the same conditions of time, temperature, water chemistry, and flow rate. Table U summarizes the treatments. These treatments were investigated to establish whether prior precipitation and stabilization of carbide phases at strain sites within the matrix could reduce the amount of carbides available for precipitation at grain boundaries during service. Tensile-Test Conditions. Specimens were tested to fracture at room temperature or 1350 F in air on an Instron tensile-testing machine at a crosshead speed of 0.05 in. per min. The time for the furnace to reach equilibrium at 1350°F was approximately 30 min; soak time prior to testing was 15 min; and test duration was approximately 5 min.

Jan 1, 1968

By Robert E. Green, Edmund G. Henneke

Expressions are derived for the orientation factors associated with dislocation damping due to dislocation motion in the basal slip plane. The results of the calculation of these factors for zinc crystals aye pvesepzted and the relevance of these results to ultvasonic-attenuation measuvements is discussed. In this paper expressions are derived for the orientation factors associated with damping of sound waves due to dislocation motion in the basal slip plane of hexagonal crystals. For zinc crystals, the orientation factors for dislocation motion in one slip direction only, as well as for simultaneous motion in all three basal slip directions, have been calculated at 1-deg intervals. The results of these calculations are presented in the form of contour diagrams in a sector of a standard polar coordinate graph. The calculations have been carried out for both the traveling- and the standing-wave conditions. The relevance of these results to ultra sonic-attenuation measurements is discussed. Associated with the calculation of the orientation factors, the modes of elastic-wave calculation were also evaluated, according to the procedure previously employed by Hinton and Green.&apos; The results of these calculations are not presented here since similar results have been reported earlier by Mus-grave.2 However, it should be mentioned that, although the qualitative features of Musgrave&apos;s results are preserved, the quantitative values are somewhat altered since different values for the elastic constants of zinc were used. ORIENTATION FACTORS It is convenient to refer to three different sets of coordinate axes, Fig. 1, in the following calculations. First, there are the Miller-Bravais axes which are used in describing the Miller-Bravais indices of planes and directions in a hexagonal crystal. These axes are denoted in the conventional manner by a,, a,, a,, and c. The second set of axes used is a rectangular Cartesian coordinate system located as shown in Fig. 1. Note that the x axis coincides with the close-packed direction al [2110], the z axis coincides with the c axis [0001], and the y axis is orthogonal to these two lying in the [0110 ] direction. Finally, a form of Eulerian angles is used to locate the position of the specimen axis. The angle between the specimen axis and the c axis is denoted by 0 and the angle between the projection of the specimen axis on the basal plane and the a, direction is denoted by ø. The majority of observations on plastic deformation of metal crystals indicates that such deformation is caused by the motion of dislocations on a limited number of crystalline planes and in a close-packed direction. In crystals with hexagonal symmetry the most prominent slip plane is the basal plane (0001) and the slip directiops- are one of the close-packed directions, [2110], [ 12101, or [ 1120]. To properly set the stage for the following work, several points must be clarified. First, we want to consider the contribution to the attenuation of elastic waves by dislocations only. While the Granato-Lucke, model for dislocation damping is presently the one most quoted, it is not necessary to refer the present work to it. The only assumption that need be made is that the motion of dislocations absorbs energy from the elastic wave and therefore attenuates it. It is immaterial here how this energy absorption occurs. However, because dislocations have been observed to move most freely on the (0001) plane for hexagonal crystals undergoing plastic deformation, it will be assumed that it is the dislocation motion on this plane which attributes to the observed attenuation. The basic orientation factor, Ok, is essentially the same as that used by Green and Hinton for fee crystals,4 defined by the equation

Jan 1, 1968

By William C. Hagel

Oxygen-18 exchange between gaseous oxygen, held at a pressure of 125 mm Hg in a PL-IORh chamber, and splzeres of a Fe2O3 containing three or less grains was determined from 9000 to 1250°C. Isotope equilibrum on crystal surfaces appears to be quickly established, and the rate-determining factor is self-difiusion confbvrning to the relation D = 2.04 esp(- 77,900/RT) sq cin sec-&apos;. Changing sphere diameters caused no detectable variation in diffusion coefficients. Oxygen ions diffuse faster than cations at low temperatures, but a crossover may occur at about 1150°C. INVESTIGATIONS by Davies et al.&apos; and Himmel et al2 have suggested that the mechanism of iron oxidation above 600°C involves a predominance of cation migration in the wüstite (FeO) and magnetite (Fe301) layers and oxygen-ion migration in the hematite (a Fe2O3) layer. Evidence for little cation migration during hematite growth consists of: 1) the disappearance of radioactive-silver markers, presumed to vaporize from the oxide-gas interface, and 2) a large discrepancy between parabolic rate constants for magnetite oxidation and those calculated from Lindner&apos;s3 data for cation self-diffusion, D:, using Wagner&apos;s theory.4 Reversing the calculations gave anticipated self-diffusion coefficients, D:, for oxygen ions of 3.6 x 10-9 sq cm sec-1 at 1100°C and 5.9 x 10-10 sq cm sec-1 at 1000°C. Since then Kingery et a1.5 reported three values for oxygen-ion diffusion in polycrystalline, sintered hematite. These were obtained near 1200°C, at an oxygen pressure of approximately 150 mm Hg, and could be represented by the Arrhenius expression, Dg = 1 x 10 11 exp(-146,000/RT) sq cm sec-1. Extrapolation to temperatures below 1100°C infers that oxygen ions are less mobile than cations. This contradicts the aforementioned oxidation evidence and also coble7s study of the initial sintering of hematite. From the latter, apparent diffusion coefficients for the slowest rate-controlling species could be calculated, and fair agreement was found with an extension of Lindner&apos;s data for cation diffusion. Fig. 1 depicts such comparisons graphically. Three points have been included for the determinations of Himmel et al. on cation diffusion in large, natural hematite crystals surrounded by pure oxygen at atmospheric pressure. No dependence on crystal log raphic direction was noted. Since the value of 10-l4 sq cm sec-&apos; at 1000°C actually represents an upper limit, one might interpret that there is a small decrease in D: when single crystals are employed. On the whole, however, cation diffusion appears relatively well-defined,* while the few results tension to lower temperatures. After having gained experience measuring oxygen-ion diffusion in certain glasses7 and in the corundum oxide, a Cr2O3, the purpose of this work was to determine D; values for hematite and thus possibly assist a mechanistic understanding of the kinetics of the oxidation and reduction of solid iron oxides. I) EXPERIMENTAL PROCEDURE A) Material. Hematite powder of moderately high purity was received from the Columbian Carbon Co. An accompanying lot analysis listed 0.03 H2O, 0.008 weight percents, with the remainder Fe2O3. Subsequent treatments were ball milling, drying, prepressing, and hot-pressing for 60 min at 1350°C under 5000 psi. Since graphite dies coated with a slurry of alumina were used for hot-pressing, the surrounding atmosphere at temperature was sufficiently reducing to form some magnetite on the outer surface of each cylinder. After removing any attached alumina, a stabilizing anneal of 50 hr at 1300°C in oxygen was added. The final product was nonmagnetic and exceeded 97 pct of the density (5.24 g cm-3) usually given for FezO3. On polishing cylindrical ends, counts were made of the number of grain boundary intercepts per

Jan 1, 1967

By W. van der Knapp

A theoretical and experimental ana1ysis is given of the change in volume of a porous medium due to changes in external and internal pressures. The result enable one to deduce directly the effect of large increlents in stresses. It is shown that the three-dimensional representation of the volume of a nonlinear porous system is a cylindrical rules surface, the generating lines of Which have a slope according to the solid rock compressibility. Analytical expression for such a ,surface is very simple. It follows that the bulk and pore compressihilities depend only on the difference between pore fluid and external hydrostatic pressure, i.e., the effective rock frame pressure. From measurements on a number of sandstones as well as limestones it can be conduded that for practical values of this effective frame pressure the pore compressibility falls in the range between the compressibility of water and that of undersaturated crude oil. In general, pore compressibility is found to be higher the lower the porosity. For a particular limestone reservoir pore compressibility and porosity could he related by means of a simple empirical formula. Results from the described study find a direct applicxtion in material balance calculations and in prob1ems of liquid flow through porous media. Since the elastic constants of the rock bulk material enter into equatiorls for velocity of acoustic waves, the results are also of importance for calculating rates of wave propagation in these media. It is shown for sandstones with 1.5 to -10 prr cent porosity that laboratory measured porosities may differ from those lender reservoir conditions by about I per cent of their value. The difference for low porosity limestotnes can be of the order of 10 per cent. The subsidence caused by elastic deformation of reservoir rock is shown to he very small. General formulas relating pore and bulk volume to external and internal pressure can be derived, as has been shown by Biot1, Gassmann2 and Geertsma.2 Their theories describe the differential change in volume caused by a variation in pore fluid tension and external hydrostatic tension. It has been known for a long timei that the stress-strain relations of a packing of spheres are essentially noninear. Our experiments show that for almost all reservoir rocks which arc elastically deformable the relation between volume and loading is also a nonlinear one, resembling that of a packing of spheres very closely. The behavior of these media is often described by using the linear stress-strain relations for small incremental stresses and strains in the vicinity of a pre-stressed condition." The present study differs from these considerations in that it enables one to deduce the effect of large increments in the stresses directly. In analogy to the stress tensor for a homogeneous body the stress tensor for a porous body can be written in the conventional way: where the tensor components, o, and t, represent the total forces acting on the faces of a unit cube of porous material. This tensor can be split into two parts. a hydrostatic and a deviatoric part, as follows2,3,1 where the component of the hydrostatic stress system, a. equals The hydrostatic stress causes a change in volume whereas the deviatoric stress system gives rise to a change in shape only. One of the difficulties in describing changes of volume occurring in reservoir rocks is that the relations connecting these volume changes with external and internal stresses are essentially non-linear. A complete description requires the bulk volume, V,,, as well as the pore volume,V,,. to be known as a function of both the changes in pore fluid tension (v) and the hydrostatic tension on the outer boundary.