Search Documents

By C. K. Eilerts

Finite difference equations were programmed and used to integrate the second-order, second-degree, partial differential equation with variable coefficients that represents the transient linear flow of gas-condensate fluids. Effect was given to the change with pressure of the compressibility factor, the viscosity, and the effective permeability and to change of the absolute permeability with distance. Integrations used as illustrations include recovery of fluid from a reservoir at a constant production rate followed by recovery at a declining rate calculated to maintain a constant pressure at the producing boundary. The time required to attain such a limiting pressure and the fraction of the reserve recovered in that time vary markedly with properties of the fluid represented by the coefficients. Fluid also is returned to the reservoir at a constant rate, until initial formation pressure is attained at the input boundary, and then at a calculated rate that will maintain but not exceed the limiting pressure. The computing programs were used to calculate the results that would be obtained in a series of back-pressure tests made at selected intervals of reservoir depletion. If effect is given to the variations in properties of the fluid that actually occur, then a series of back- pressure curves — one for each stage of reserve depletion — is required to indicate open-flow capacity and related flow characteristics dependably. The isochronal performnce method for determining flow characteristics of a well was simulated by computation. INTRODUCTION The back-pressure test procedurel is based on a derivation of the equation for steady-state radial flow of a gas, the properties of which are of necessity assumed to remain unchanged in applying the test results. The properties of most natural gases being recovered from reservoirs change as the reserve is depleted and pressures decline, and the results of an early back-pressure test may not be dependable for estimating the future delivery capacity of a well. The compressibility factor of a fluid under an initial pressure of 10,000 psia can change 45 per cent and the viscosity can change 70 per cent during the productive life of the reservoir. There are indications that the effective permeability to the flowing fluid can become 50 per cent of the original absolute permeability before enough liquid collects in the structure about a well as pressure declines to permit flow of liquid into the well. The advent of programmed electronic computing made practicable the integration of nonlinear, second-order, partial differentia1 equations2.3 pertaining to flow in reservoirs. Aronofsky and Porter4 represented the compressibility factor and the viscosity by a linear relationship, and integrated the equation for radial flow of gas for pressures up to 1,200 psi. Bruce, Peaceman, Rachford and Rice5 integrated the partial differential equations for both linear and radial unsteady-state flow of ideal gas in porous media. Their published results were a substantial guide in this study of integration of the partial differential equation of linear flow with coefficients of the equation variable. The computing program was developed to treat effective permeability as being both distance-dependent and pressure-dependent. In this study all examples of effective permeability are assumptions designed primarily to aid in developing programs for giving effect to this and other variable coefficients. The accumulation of data for expressing the pressure dependency of the effective permeability is the objective of a concurrent investigation. FUNDAMENTAL DIFFERENTIAL EQUATIONS DIFFERENTIAL EQUATION OF FLOW Linear transient fluid flow is represented by the partial differential equation,

Jan 1, 1965

By W. R. Kube, W. H. Oppelt

Five low-rank coals, including two lignites, a steam-dried lignite, and two subbituminous coals, were carbonized at 940°F, in a bench-scale carbon-ize~ with a nitrogen and hydrogen atmosphere, or both, using pressures from atmospheric to 1000 psig. Coking characteristics were observed, and yields, composition, and distribution of products from thermal degradation were determined. The influence of pressure on yields from coal carbonization, as well as on modification of coke structure, has long been of interest. One of the first attempts at pressure carbonization was that of Capps and Hulett1 who in 1917 investigated the effects of pressure ranging from 1 to 20 atm on yields at 600° C for bituminous and subbituminous coals. High-pressure carbonization has not been used commercially for coke production. More recently, carbonization at relatively high pressure in the presence of catalysts and hydrogen has been proposed for the production of metallurgical-grade coke from low-rank fuels.' Another approach involves carbonization at high pressure in a hydrogen atmosphere to produce heating gas and a solid residue that is gasified with oxygen and steam to produce the hydrogen required for hydro-genation 3-4 or, alternatively, is burned to raise steam. Recently it has been proposed to carbonize coal under pressure, burn the gas in a gas turbine to produce power, and use the solid residue as fuel in a conventional steam-power plant.5 Our primary interest in pressure carbonization stems from an attempt to determine the influence of pressure on the coking characteristics of lignite to design a high-pressure slagging gasifier using oxygen and steam. At atmospheric pressure, North Dakota lignite is generally considered to be non-coking. However, lignite and subbituminous coals have been reported to show increased agglomeration properties in steam, hydrogen, or carbon-dioxide atmospheres at elevated pressures, although no strong coke was formed at pressures up to 3000 psig.6 If lignite would agglomerate seriously under gasification conditions, proper flow of charge through the gasifier could not be maintained without installing a mechanical stirring device or destroying the caking properties of the coal before gasification. Data on the caking properties of lignite or other gasification fuel, as well as on yields of products under conditions simulating the carbonization process in the pregasification zone, were required in the preliminary design stage of the gasifier. No information was available concerning the influence of pressure and type of atmosphere on North Dakota lignite. To provide this basic data, the U. S. Bureau of Mines investigated the effect of pressure and a hydrogen or nitrogen-rich atmosphere on the yield of products and on coke formation in a bench-scale pressure carbon-izer. SCOPE OF INVESTIGATION Five coals were tested at 940°F in nitrogen and hydrogen atmospheres, at pressures from atmospheric to 1000 psig. Yields and composition of solid, liquid, and gaseous products were determined as a function of the process pressure, and attention was given to the physical appearance of the solid residue. Percentage of carbon converted, sulfur retained in char, and thermal balances were calculated. In this paper, primary consideration is given to carbonization in a nitrogen atmosphere and to data concerning the solid phase. The effect of pressure in hydrogen is presented in detail for only one lignite. More complete information on the response of the other coals in a hydrogen atmosphere will be published in a USBM Report of Investigations. DESCRIPTION OF EQUIPMENT Details of the bench-scale pressure carbonizer are given in Fig. 1. The body of the carbonizer was

Jan 1, 1961

By M. R. Tek, M. L. Katz

The flow of fluids in stratified porous systems, in which adjacent layers of the system are in communication over part or all of their common interface, in many instances involves crossflow between layers. A theoretical study has been made of the unsteady-state flow of slightly compressible fluids during depletion of bounded stratified porous systems. Consideration has been restricted to single-phase flow, with negligible capillary and gravity effects. Analytical expressions have been derived for pressure distribution and for (luxes at the producing face in two-layer systems with complete fluid communication between layers. These expressions occur in the form of double Fourier and Fourier-Bessel series. The method of extending the solution to any number of layers is demonstrated by deriving the equation for pressure distribution in a three-layer system. On the basis of the exact results for two-layer systems, certain generalizations regarding the behavior of stratified systems in general are presented. The performance of stratified systems, in terms of cumulative flux as a function of time, is shown to lie at all times between bounds established from single-layer theory. The upper bound is given by treating the system as a single layer with arithmetic average physical properties. The corresponding lower bound is the summation of the fluxes from each layer as treated individually. Variation in performance of stratified systems between the upper and lower bounds is determined by the extent to which crossflow occurs in the system. The principal parameters affecting the amount of crossflow are the vertical permeability of the system, the ratio of system thickness to "radius of drainage" and, for radial systems, the ratio of system thickness to interior radius. At very early times, the system performance is near the lower bound. As time and the radius of drainage increase, the per- formance tends toward the upper bound. INTRODUCTION One of the principal problems confronting the reservoir engineer today is the incorporation of the effect of reservoir heterogeneity into predictions of reservoir performance. The first part of this problem is the characterization of the actual reservoir heterogeneity in a mathematically useful manner.1, 2 It is then necessary to calculate the effect of this heterogeneity on the system performance. It is frequently possible to estimate this effect by comparison with the known behavior of an idealized system. One such idealized system is the stratified or "layer-cake" model. Many studies have been made of the behavior of stratified systems, from both a theoretical and experimental standpoint. The great majority of the theoretical studies have been concerned with systems of separated layers, with no fluid crossflow. Among the more noteworthy of such theoretical papers is that recently published by Lefkovits, et al,3 in which a comprehensive analysis was presented for the behavior of bounded stratified systems which are unconnected except at the wellbore. Three recent references have dealt with stratified systems including the effect of crossflow. Jacquard4 has presented mathematical derivations pertaining to depletion of two layer systems producing at constant rate. A doctoral dissertation by Katz5 contains analytical solutions and tables of results for pressure distribution and fluid flux during constant-terminal-pressure depletion of two-layer and three-layer systems with crossflow, together with tables of eigenvalues necessary in the computations. The dissertation also presents results from experimental studies of stratified systems using a heat-transfer analog model. This paper is largely based on the material in this dissertation. A third paper has been presented by Russell and Prats6 since completion of the work resented herein, containing mathematical derivations for and discussion of the behavior

By R. N. Thurston, E. A. Nesbitt, G. Y. Chin

A general relalionship between the amount of glide shear (due to slip) and the macroscopic shape change has been developed. Since the deformation can be large, finite strain analysis is employed. In this treatnzent, the shape change is expressed by a deformation gradient matrix, F - [axi/aXj1, where X and x refer to the initial and final positions of a particle. This matrix is readily evaluated in terms of the amounl of glide shear, a , and the direction cosines of the slip-plane normal and the slip direction for a single active slip system. In the case of deformation from slip on several systems, the product F(a) of the several deformation gradient matrices is first calculated. Then, by assuming that the final configuration is reached by a long series of small shears of magnitude a, occurring more or less alternately in the several slip sys- IN problems of crystal plasticity, it is often necessary to relate the amount of glide in the operating slip systems to the macroscopic strain components. This relationship is relatively simple when only small strains are considered. In this case the separate strain contributions from several slip systems are additive. Such a procedure is incorrect in the case of large plastic deformation. To the authors' knowledge, no general treatment of the latter problem has appeared in the literature. In specific instances, Mark, Polanyi, and schmid3 have derived the relationship between glide and axial elongation during tensile pulling of a single crystal when only one slip system is active. A similar relationship for duplex slip was worked out by v. Goler and Sachs.4 Taylor and Elam5 have also studied in detail problems of large plastic deformation. Due to uncertainties as to the slip systems at the time, terns, the final deformation gradient matrix can be obtained. Mathematically, the problem is reduced to finding the limit of F(dN as N — m while a — 0 in such a way that the Product Na = 0, a ftnite conslant designatzng the accunzulated amount of shear. It turns out that this limit is simply , where F, is the matrix whose elements are the coefficients of a in the F(a) matix. Application is made of the present treatment to a fcc crystal of Perlrzalloy compressed on the (110) plane and contrained to elongate zn the [112] directzon. In addition, it is shown that one rnay readily obtain from the general analyszs the well-known formulas relating elongation, amount of glide shear, and amount of lattice rotation for crystals deforming by single and double slip under tension. however, they were mainly concerned with proving that slip occurs on {111}(110) systems in fcc metals. The present general treatment arose out of recent investigations of magnetic anisotropy in cold-worked Fe-Ni alloys6 as well as some related problems in strength anisotropy of single crystals.' Detailed application of the general analysis will be provided in the case of a Permalloy single crystal compressed on the (110) plane and constrained to elongate only in the [112] direction. It will be seen that, as expected, the small strain approximation leads to significant errors after moderate straining. It will also be shown that the present general treatment yields the formulas of Mark et al. and of v. Goler and Sachs in those specific cases. GENERAL CONSIDERATIONS Consider a homogeneous deformation in which a material point initially at moves to both positions being referred to the same set of Cartesian axes. Such a deformation can be

Jan 1, 1967

By R. H. Hiltz, R. W. Douglass

Large-grain specimens of iodide titanium prepared metal-lographically were stain etched using the technique of New York University as modified by Watertown Arsenal Laboratories. Orientations of grains showing an abnormal response to stain etching were determined by back-reflection X-ray techniques. X-ray results showed that a-titanium grains, oriented such that their basal planes were parallel to the polished surface, have a slower reactivity to stain etching than any other orientation. Grains oriented such that their basal planes made small angles with the polished surface exhibited reaction rates between that of the true basal orientation and the normal rate. As the angle with the surface increases the reactivity increases up to a maximum angle of 14 deg. Basal planes making an angle with the surface of greater than 14 deg show normal reactivity to stain etching. These data, plus work on polycrystalline specimens, indicate that the oxidation rate of titanium might be reduced by correct alignment of crystal planes. ThE earliest observed phenomenon in the micrography of titanium was the occurrence of a grains which exhibited no birefringency. Even though a titanium possesses a hexagonal-close-packed structure and should be optically anisotropic, it can appear optically isotropic, when the crystal is so oriented that its intersection with the polished surface is perpendicular to the axis of symmetry. During the develoment of color staining techniques for titanium: it was noted that the a phase did not stain uniformly. Staining of specimens of commercially pure titanium resulted in a micro-structure of blue-colored grains, among which were randomly oriented a few grains exhibiting colors between the initial yellow and the final blue color. Under polarized light, these off-colored grains appeared to be optically isotropic. More intense study with polarized light, however, showed only the yellow grains to be isotropic. The intermediate colored grains exhibited weak birefringency, which at first glance is not detectable due to the strong birefringency of the surrounding structure. Aminoff2 and Straumanis3 have previously shown that in zinc the basal orientation is less reactive than other orientations. Leidheiser4 showed the same thing to be true for the hexagonal-close -packed configuration of cobalt. Thus, it would appear that the basal-plane orientation is responsible for the abnormality in stained titanium. Later work by Leidheiser5 on cobalt has shown that other planes besides the basal plane exhibited reaction rates slower than the normal rate. Reported planes were the (100), the (1011) and several planes making small angles with the basal plane. The last observance appears not to be due to the reactivity of the reported plane but due to the reactivity of the basal plane which is very close to being parallel to the surface. There are, then, two possible explanations for the range of colors observed in stain-etched a titanium. Since extensive research on the effect of orientation on surface reaction has failed to provide any definite trend, these two explanations are not adequate in themselves. Planes showing the retarded reactivity differ not only for different crystal structures but for metals of the same crystal structure and even for different reactions of the same metal. Tragert,6 Leidheiser,7 Tammann,6 and wagner9 have shown copper to exhibit retarded reactivity on the (111) plane. However, Leidheiser10 has shown that under other conditions the retarded reactivity of copper is associated with the (100)

Jan 1, 1960

By B. L. Averbach, Morris Cohen, S. A. Kulin

The martensitic transformation can be initiated by elastic stresses at temperatures above M. in a steel containing 20 pct Ni and 0.5 pct C. Shear strains and normal tensile strains acting on a potential habit plane promote the transformation but compressive strains oppose it. THE strain sensitivity of the martensitic transformation has long been recognized. It is now well known that martensite formation can be induced by plastic deformation at temperatures below, and not too far above, the M, point.' However, the exact role of elastic vs. plastic strains and of the concomitant stresses has proved to be a very elusive matter. Many investigators2-' of the steel hardening reaction have made use of the stresses resulting from volume changes to account for such observations as the existence of retained austenite and its variation across the section of a bar, the change in the amount of retained austenite with quenching rate, and the self-stopping nature of the transformation when the cooling is stopped. The effect of applied strain as an independent variable has been studied in 70 pct Fe-30 pct Ni alloys by Scheil7 and McReynolds.8 Some data are also available for lithium and Li-Mg alloys,9 Cu-Zn,10 Cu-Sn and Cu-A1 alloys," and aus-tenitic stainless steels."-" It has been found that martensite formation can be induced isothermally, even at temperatures above M,, by mechanical deformation of the parent phase. The competitive nature of plastic yielding by slip and by the martensitic reaction has been depicted by Scheil.7 Basically, Scheil postulated that: 1—at temperatures not too far above M, (where austenite is less stable than martensite), a critical resolved shear stress (within the elastic range) is required to promote the transformation to martensite; 2—at and below M,, the austenite lattice becomes mechanically, as well as thermodynamically, unstable and shears over spontaneously into martensite without the application of external stress; and 3—at temperatures sufficiently far above M8, the critical resolved shear stress for martensite formation increases to a level above that required for slip, and plastic flow then supersedes the transformation when external stress is applied. Scheil's concepts were tested by McReyno1ds8 ho found no change in the elastic moduli (measured both statically and dynamically) of 71 pct Fe-29 pct Ni alloys in the vicinity of the M, temperature. Since the moduli did not approach zero or become negative on cooling to M,, it was concluded that the austenite lattice does not become mechanically unstable at M2, as postulated by Scheil. McReynolds also reported that M, is not raised by elastic stresses. Accordingly, plastic deformation was considered to stimulate the transformation by generating martensite nuclei in the distorted regions of the parent phase. The above issues with regard to the role of applied stress require clarification, if a basic understanding of the martensitic transformation is to be obtained. For example, in the reaction-path theory,'"," it is postulated that strain embryos exist in the austenite that provide part of the activation energy for the nucleation process. These embryos are visualized to be regions of localized strain in which the atoms in the austenite lattice are displaced part way along the path to their ultimate positions in the martensite. Consequently, it would

Jan 1, 1953

By Freeman Bishop

President Allende says he will expropriate all American investments in Chile and he's already well on the way to accomplishing this objective. What the average observer doesn't realize is that this has not only created big lossess for the copper companies. It could also take a chunk out of the U.S. taxpayers' pockets. Premier Castro expropriated the U.S. taxpayers' investment of over $100 million in the Nicaro (Cuba) nickel and cobalt operations-and U.S. Government protest over the loss was never heard above a whisper. Where the taxpayer comes in today is through the Overseas Private Investment Corporation, which was set up to guarantee U.S. companies' foreign investments against losses due to political changes. As of now, metals firms hold more than $300 million worth of obligations against this U.S. Government agency. The agency itself may not be able to raise more than $85 million, and a supplemental appropriation may be needed to meet its guarantees.

Jan 1, 1971

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 John W. Weeton, Max Quatinetz, Robert J. Schafer

Nickel powders of 1-, 0.4-, and 0.2- average particle size were combined with 0.05- p MgO powder. Oxide was added in quantities of 4, 12, and 20 vol pct. The mixtures were vacuum hot pressed and extruded at 2100°F. In spite of a general tendency to growth of oxide particles during processing, dispersions of fineness comparable to those of Sintered Aluminum Powder, SAP, were achieved. A small amount of oxide considerably improved 1800°F stress-rupture strength, but further oxide addition resulted in a decrease of strength. RECENTLY there has been much speculation concerning the requisites for strong and stable dispersion-strengthened alloys. The nature of the dispersed phase, the stability of the alloys, and variables influencing the mechanical properties of alloys have been given consideration. Such topics have been reviewed by Bunshah and Goetzel. Of the many parameters related to the mechanical properties of dispersion alloys, mean free path between dispersed particles (a measure of the fineness of the dispersion), particle size, and vol pct of the dispersed phase have been shown to be very important. Generally, strength increases with decreasing mean free path, and there is an optimum vol pct of dispersed phase. Just how fine a dispersion ought to be is somewhat open to debate, but Irmann, Lenel, Backensto, and Rose, Gregory and Grant,4 and others have shown in the original SAP-type alloys, aluminum plus aluminum oxide alloys, that better strength is generally obtained with finer dispersions and the best strengths were obtained when the mean free path was less than . Further, Cremens, and rant' showed that in nickel plus aluminum oxide alloys, stress-rupture strength at 1500°F was definitely increasing with decreasing mean free path. In their investigation, however, the finest dispersions that were produced were about 2.0 mean free path. To determine the potential of the nickel plus oxide type of product, it would be neces- sary to obtain dispersions analogous to those of SAP alloys; namely, ones where the mean free path between oxides is less than . Unfortunately, however, it is very difficult to produce extremely fine dispersions of oxides in metals other than those which naturally form refractory oxide skins. There are several approaches to this problem including such procedures as mechanical mixing, internal oxidation, reduction of mixed oxides, and metal deposition from solution. Work at the Lewis Research Center has shown that metal powders 0.1p or finer could be produced by ball milling, when selected grinding aids were utilized.' Because these unusually fine powders were available, it was an initial objective of this study to determine if by mixing of these fine Ni powders with various amounts of MgO (0.05) a fine dispersion of MgO in Ni comparable to those in Al-Al203 might be achieved. A second objective was then to attempt to relate the parameters of structure to 1800°F stress-rupture strength of the alloys. Nickel powders of 1-, 0.4-, and 0.2-p average particle size were combined with 0.05- average particle size MgO powder. This oxide was added in quantities of 4, 12, and 20 vol pct. The mixtures were hot pressed and extruded at 2100°F.

Jan 1, 1962

By F. I. Stalkup

Vapor-liquid phase equilibrium experiments have been conducted in a static equilibrium cell on mixtures of a light, 45 API stock- tank gravity reservoir fluid and a rich hydrocarbon gas containing approximately 55 mole per cent of intermediate hydrocarbons. Both apressure-vs-composition study of the gas and a simulated reservoir fluid, and a multiple-batch contact simulation of the condensing-gas-drive oil recovery process were performed. In the latter experiments equilibrium gas and liquid compositions were analyzed. Also, conventional, "condensing-gas-drive", long-tube displacement experiments of the reservoir fluid and gases of various richness were performed. The results of these experiments could not be satisfactorily interpreted by the conventional pseudo-ternary-diagram representation of multicom-ponent phase behavior. The results seem to be explained better by considering a bubble-point surface and a dew - point surface joined in a plait-point locus. Portions of the plait-point locus cannot be "seen" directly by the rich hydrocarbon gas because of curvature of the dew-point surface. In such a system, continuous injection of the rich gas over stationary reservoir fluid might form a zone of contiguously miscible compositions from pure rich gas to pure reservoir fluid by: (1) saturating the reservoir fluid with injected gas to the bubble-point surface; (2) creating by mass transfer with fresh injected gas a path of contiguously miscible compositions along the bubble - point surface to the plait-point locus; and (3) creating by mass transfer with additional injected gas a path of gas compositions along the dew-point surface up to the point where direct miscibility results between dew-point fluid and the injected rich gas. INTRODUCTION The use of the pseudo-ternary-phase diagram to illustrate miscible displacement phase behavior has been discussed by several authors.5,7 Such a representation of phase behavior is not rigorous, but the ternary diagram nevertheless gives a qualitative picture of what actually occurs in a miscible displacement process. Fig. 1 is a typical illustration of miscible displacement phase behavior by a ternary diagram. The multicomponent hydrocarbon system is divided into three pseudo-components: a light fraction containing methane and nitrogen, an intermediate fraction containing ethane through hexanes plus carbon dioxide, and a heavy fraction containing heptane and heavier components. A two-phase region is bounded by a dew-point curve and a bubble-point curve, which are joined at the critical point. The concept deduced from such a representation for miscible displacement by a condensing-gas-drive process is as follows: a rich gas G, which lies to the right of the limiting tie line through the critical point C, is injected into the reservoir and contacts reservoir fluid L, saturating the reservoir fluid to give bubble-point fluid L1 and equilibrium dew-point gas GI. Continued injection of rich gas changes the composition of the saturated liquid L1 through a series of liquid compositions lying along the bubble-point curve, until the critical composition C is reached, at which point direct miscibility with the rich gas is achieved. Some equilibrium gas with compositions lying along the dew-point curve from G1 to C is also formed in this process. Thus, continuous injection of rich gas forms a

Jan 1, 1966

By Paul A. Beck, M. N. Parthasarathi

By determining the (220) pole figure for OFHC copper rolled to 96 pct R. .A., the occurrence of four texture components of the type (135) [211] was confirmed. It was found that the total volume fraction of material close to the four symmetrical orientations of this type is larger than the volume fraction of the other components put together. THE rolling texture of copper was described by early investigators1'" as a mixture of components of the (110)[l12] and (112) [11l] types. Other studies, still using the qualitative photographic method of texture determination arrived, instead, at orientations of the type (135) [533] or (135) [211] for copper3'4 and at the closely related "Z-texture" for a fee nickel-iron alloy.5 The prevalence of this type of orientation, near (123) [412], in the rolling texture of aluminum and of copper was later confirmed by means of quantitative diffractometer methods.' However, in a recent paper7 Jones and Fell stated that the (110) [112] and (112) [111] double texture was in better agreement with the measurements of Young's modulus for cold-rolled copper.' These authors also expressed their belief that the "roll texture" can be determined by means of two photographic quasi-fiber texture patterns taken with sheet specimens rotated around the rolling direction and the transverse direction, respectively. Jones and Fell found7 that the results obtained from such patterns for cold-rolled nickel did not agree with any of the orientations mentioned above, but they suggested that the observed diffraction spots may have been brought about by the overlapping of reflections due to the (110) [112] and (112) [11l] double texture components. Jones and Fell concluded that the (123) 14121 type orientation is inconsistent with their patterns. They also stated, without offering any supporting evidence, that this orientation is inconsistent with the (220) pole figure. "Ideal orientations" are often used for a very rough description of the principal orientations encountered in a texture, It is, however, well to remember the fact that, no matter how "accurate", these ideal orientations are only schematic representatives of a continuous distribution9, which often (as in the case of rolled copper) includes extensive orientation spreads. The attempt by Jones and Fell7 to use elastic modulus data in deciding between various ideal orientations in the rolling texture of copper neglects this fact. The inadequacy of two-dimensional "axis density figures" in describing orientation distributions in sheet textures was previously discussed.9 It was also pointed out9 that there is no reliable evidence to prove that orientation distributions derived by the "quasi-fiber method" from rotated sheet specimens are fully equivalent to the corresponding information obtainable from quantitative pole figures. The view that a sheet texture could be adequately determined by means of two photographic quasi-fiber texture patterns is obviously erroneous9. The question as to the consistency with the (220) pole figure of ideal orientations derived from the (111) and (200) pole figures is, on the other hand, certainly very relevant. Since no quantitative (220) pole figure for cold-rolled copper appears to be available in the literature, the following work was undertaken. EXPERIMENTAL METHOD An extruded rod of OFHC copper, 1 in. in diam, was given alternate rolling (30 to 35 pct R.A.) and

Jan 1, 1962

By W. J. Joslin

The frontal-advance equation can determine how the fluid withdrawal rate and subsurface operating pressure influence oil recovery from pressure-maintained reservoirs having characteristics favorable for efficient vertical segregation. For effective vertical segregation in formations having low dip angles, the requirements are: thick, massive sands, low viscosity oil and high vertical permeability with few discontinuous shale members. This approach has been successfully applied to the LL-370, which is a massive Eocene sandstone reservoir containing crude with an average gravity of 26" APT, a dip of only 3" and vertical permeability of ½ darcy. The average gas cap saturation was calculated for various rates of vertical gas-oil contact movement corresponding to a wide range of reservoir producing rates. This method also evaluated the effect of different pressure maintenance levels on oil recovery. Applicability of this procedure for calculating the present oil saturation in the LL-370 gas cap has been substantiated from field performance, where the oil saturation was computed by comparing oil recovery with the gas cap volume vacated. INTRODUCTION During the natural depletion history of the LL-370 reservoir, located in Lake Maracaibo, Venezuela, excellent vertical segregation resulted in the formation of a thin gas zone over a substantial portion of the oil-productive area. Subsequently, gas injection forced this thin gas-saturated zone slowly downward, reducing the oil saturation in the gas cap to a low value. The oil recovery has been raised considerably above earlier predictions where vertical segregation effects were under-estimated. In order to determine the maximum efficient producing rate and optimum operating pressure for this reservoir, a vertical gas frontal-advance analysis has been made through application of the fractional flow equation. LL-370 RESERVOIR CHARACTERISTICS The reservoir consists of a truncated monocline of massive. highly-permeable Eocene sandstone, called the B-b formation, covering an area of more than 13,500 acres. The original oil in place was 2.17 billion STB. Internal non-sealing faults subdivide the structure into four blocks called the B-6-x.6, 10, 11 and 13, which are effectively bounded by major faulting on the flanks, an unconformity updip to the west, and an immobile aquifer to the southeast. The producing formation occurs at an average depth of 5,250 ft subsea and has a closure of more than 1,500 ft. A typical producing section has a net oil sand thickness of 170 ft. Scattered throughout the section are a few randomly-positioned shale lenses which increase in thickness and areal extent downdip. The reservoir rock and fluid properties vary with depth. The permeability and produced oil gravity are 1600 md and 28" API at the top, compared to 300 md and 18° API at the water-oil contact. The productivity of updip wells is, therefore, much higher than the reservoir average. The reservoir oil was originally saturated at the crest, grading to a slight degree of undersaturation at the water-oil contact. There was no original gas cap; and material balance and the field performance indicate no water influx. Except for the northern B-6-x.6 block, the fluid and rock properties have already been published.' Average LL-370 reservoir properties are shown in Table 1. RESERVOIR PERFORMANCE UNDER NATURAL DEPLETION Fig. 1 shows production performance from 1939 to date. Early development was gradual, with 30 wells pro-

Jan 1, 1965

By J. Lemish

Considerable research has been conducted in Iowa in recent years on carbonate aggregate reactions in concrete. Groups outside the state have become active in the past few years after recognizing that certain carbonate aggregates are reactive in concrete. The paper reviews the nature and extent of the problem and a summary is given of research by the author and his associates. In recent years considerable research has been conducted in Iowa on carbonate aggregate reactions in concrete. The research is part of a general study on the behavior of carbonate rocks as concrete aggregates sponsored by the Iowa State Highway Commission. Research on carbonate aggregates has not been extensive in the past. Laughlin1 published one of the first papers on carbonate aggregates in 1928 and was followed by others, notably H. S. sweet,' B. Mather et al.,3 and K. Mather4 who emphasized the petrographic and physical aspects of carbonate aggregates. The recognition that certain carbonate aggregates are reactive in concrete has led to increased participation in research in the past few years by groups outside of Iowa. Significant contributions have been made by Swenson and Gillot,5 Hadley,6 and Chaiken and Halstead7 in work which emphasizes the expansion of carbonate rocks in alkaline environments and the reaction of carbonate minerals or dedolomitization with alkalies. It is the purpose of this paper to review briefly the nature and extent of the problem and present a summary of the research by the author and his associates. EXTENT AND NATURE OF THE PROBLEM The problem is the deleterious behavior of certain carbonate aggregates in concrete. It was first recognized in Iowa over 20 years ago when poor service records of certain highways were correlated with aggregates produced from the Glory and LeGrand quarries. Deleterious behavior of carbonate aggregates have been reported in Canada and, more recently, in New York and Virginia. As more knowledge is gained concerning the problem, distress in concrete related to carbonate aggregates may be recognized over a widespread area. The problem may potentially exist in the Midwest, the Northeast, and other parts of the country underlain by thick sections of carbonate rocks readily available for utilization as aggregates. The nature of distress in Iowa highways is a general weakening and deterioration of concrete characterized by progressive cracking and spalling away from joints and edges (Fig. 1). The distress generally occurs after 15 years of service and becomes progressively worse. Specimens of such concrete lack the toughness of sound concrete and break easily under the hammer in the manner of a soft brick. The matrix is generally highly carbonated and dark borders, or reaction shells or rims outline the argillaceous carbonate aggregate fragments (Fig. 2). Cracks observed in cores cut across the aggregates and concrete matrix. Although localized volume change (expansion) is associated with spalling, there is no definite field evidence in Iowa that excessive growth of concrete slabs has occurred. The experience with carbonate aggregates in Iowa indicates that the concrete deterioration is not the same as that induced by the alkali-aggregate reaction associated with siliceous aggregates. The principal difference is that in Iowa the distressed concrete is weakened throughout the slab. In concrete affected by the alkali aggregate reaction, the

Jan 1, 1961

By D. R. Parrish, W. E. Lamoreaux, W. W. Owens

Several methods are now being used by the industry for determining the gas-oil flow characteristics of reservoir rock samples. Most of the laboratory experirnerltal rlzethods can be classified either as steady-state, unsteady-state, or Hassler-type tests. Although the validity of the results obtained by any specific method is not questiorzed, in general, the complexity of these methods leaves much to be desired. This investigation dealt primarily with the development of a simplified unsteady-state (gas drive) laboratory test method for measuring kg/ko re1ationships of small core samples. A procedure for selecting uniform and representative sandstone rock samples was developetl since it was observed that the reliability of the test resu1ts was dependent upon having a sample which was uniform in porosity and permeability. The importance of factors known to influence the results of this type test, such as oil viscosity, gas expansion, and pressure differential were eva/uated. The method of calculating the test results was found to give valid data over a wide range of gas expansion. The limits of pressure differential and oil viscosity which would insure reliable results were defined. In cornparisorz with other presently used testing methods, the gar drive method requires considerably less apparatus, is less time consuming, and is more readily adaptable to a routine type core test. INTRODUCTION The detcrmination of the gas-oil flow characteristics of reservoir rocks has been given concentrated study by the oil industry for the past decade. Both analytical and experimental methods have been considered. The analytical methods have generally relied on a simplified characterization of porous materials in order that the fundamental laws of fluid flow and capillarity could be utilized. More direct experimental methods have been proposed by numerous production research laboratories with the primary goal being a simple and inexpensive, yet reliable method for measuring the gas-oil flow characteristics of reservoir rock samples. The experimental methods which have been reported in the literature generally can be grouped into one of the following classifications: 1. Steady state1,2,3,4 in which a fixed ratio of flowing fluids is forced through the test sample until equilibrium conditions of saturation and pressure gradient are established. 2. Hassler or stationary liquid method5,6,7 in which a uniform liquid saturation is attained, usually through capillary desaturation, then gas or oil is flowed through the sample at a pressure drop insufficient to disturb the pre-established saturation. :. Non-steady state method1,8,9,10 In which only the gas phase is injected into the sample, causing displacement of the oil phase and resulting in a system of continually changing average saturation and saturation gradient. The steady-state and Hassler methods have been most popular because, no doubt, of the simplicity of the actual relative permeability calculations. The problems involved in these methods were limited to apparatus design and measurement techniques. In the steady-state method, for instance, a constant flowing fluid ratio is established which, at equilibrium, results in an essentially uniform saturation over the entire length of the test sample. Therefore, the basic flow characteristics could be obtained by correlating the average saturation of the test sample with the kg/ko determined at various flowing gas-oil ratios.

Jan 1, 1957

By P. A. Beck, Hsun Hu

As described by means of quantitative pole figures, the annealing texture of highly rolled aluminum consists of the four retained components of the rolling texture near (123) [121], rather more sharply developed, and of a cube texture component. Local reorientation corresponds fairly well to 40 rotation around a [111] axis. In copper strip rolled 96 pct, the annealing texture is mainly the cube texture, with the four twin orientations as minor components. The annealing texture of highly rolled brass strip consists of four components of the (225) [734] type. REVIEWS by Wassermann,' Barrett,' Dunn,3 and Richards' give detailed accounts of the numerous publications dealing with annealing textures in rolled polycrystalline metals. Although a considerable amount of information of an empirical nature has been available for many years, even relatively recently the state of recrystallization-texture theory was considered uncertain. The reviews referred to above deal mainly with the description of experimental results, without developing generally applicable principles, suitable for explaining the observed annealing textures. The investigations of Burgers and Louwerse5 and Barrett6 showed that in deformed single crystals of aluminum the orientation of the recrystallized grains deviates drastically from that of the deformed matrix. According to Barrett and to Beck and Hu,7 the relationship between the two orientations can be best described as a rotation of about 40° around a [111] axis. Recent investigations led to similar orientation relationships for various face-centered cubic metals with a simple texture, regardless of whether the annealing process was recrystallization (strain-free grains growing in a highly deformed large crystal of aluminum' or copper: or in deformed cube texture copper") or coarsening (strain-free grains growing at the expense of an essentially strain-free, fine grained, polycrystalline matrix of single orientation texture copper,'" aluminum,' or nickel-iron). Barrett" pointed out the possibility that this orientation relationship may be a result of the orientation dependence of the rate of grain boundary migration. More recently, considerable evidence became available in favor of this view.12,13 In apparent contradiction to the above results, it has been known for a long time that the recrystallization texture of polycrystalline metals is often similar to the deformation texture. This was found, in particular, for the fiber textures obtained by wire drawing, compression, or compression rolling.' Following a suggestion by Barrett,2 Kronberg and Wilson and Beck and Hu7 described a mechanism by which a [Ill] fiber texture, resulting from plastic deformation, may be retained upon recrystallization, even though the orientation of individual recrystallized grains differs from that of the local area of the matrix in which they grow. This mechanism assumes that the reorientation corresponds to a rotation around the [1ll] axis, which is also the fiber axis. The present investigation was undertaken to examine in detail the orientation relationship between deformation texture and recrystallization texture in rolled face-centered cubic metals, using the quantitative methods of texture determination that have recently become available.", '" Quantitative pole figures for the rolling textures of such metals were reported recently." Quantitative pole figures for annealing textures are given in the present paper. It was planned to investigate whether or not the orientation relationship previously found for the recrystallization of deformed single crystals, and for coarsening of single-orientation texture material, is capable of accounting for the observed, often complex, annealing textures in rolled polycrystalline metals. The question whether local reorientation is predominant, even in the early phases of the annealing process, was to be studied with high purity aluminum. by making use of an oxide film and sensitive tint illum-

Jan 1, 1953

By E. L. Dougherty, J. W. Sheldon

Using the numerical techniques shown in this paper it is possible to compute the simultaneous dynamical behavior of multiple fluid-fluid interfaces in two dimensions. Hence, fluid-fluid inter face models of several oil recovery processes can be constructed which allow prediction of the effect of fluid flow in two dimensions on the behavior of these processes while taking into account various physical effects, such as saturation gradients, phase changes and thermal stimulation. Techniques for constructing fluid-fluid interface approximations for several oil recovery processes are reuiewed. The validity of these techniques is established by comparison to experimental results. The availability of a computer program for computing the behavior of multiple fluid-fluid interface problems makes the use of potentiometric models uneconomic. The results indicate the areal sweep efficiencies obtained from potentiometric model studies at high mobility ratios may be somewhat optimistic. The results also indicate that the areal sweep efficiency curves obtained at high viscosity ratios 2—4 using miscible fluids in porous media are misleading and that considerable care must be taken in defining the mobility ratio for a displacement process, especially for miscible fluids. INTRODUCTION Many of the methods used to estimate areal sweep efficiency assume that the oil recovery process can be represented by one or more fluid-fluid interfaces. When the mobilities of the displaced and displacing phases are assumed equal, this quantity can be computed using analytical mathematical methods.l.2 For non-unit mobility ratio, experimental methods based upon the fluid-fluid interface concept include potentiometric and electrolytic blotter models. Scaled models of porous media have also been used to study this problem. Fay and Prats3 and Aronofsky4 used numerical techniques to compute breakthrough sweep efficiency, but neither of these studies were extended beyond breakthrough. Numerical techniques which have been developed for solving a general fluid-fluid interface problem are presented in a companion paper.5 These techniques have been incorporated into a program for a large-scale digital computer which can treat simultaneously the dynamical behavior of as many as six fluid-fluid interfaces. Thus it is possible to simulate the behavior of several oil recovery processes taking into account such effects as saturation gradients. Techniques for representing several oil recovery processes by fluid-fluid interfaces are presented. Processes considered are conventional water flooding with and without a mobile gas phase, hot water flooding, miscible flooding and enriched gas drive. Other processes, whose effects could be studied by this technique, but are not considered in detail here, are underground combustion and steam injection processes and various slug processes such as gas-water injection and alcohol-solvent processes. Theoretical results obtained with single interfaces for mobility ratios of 0.25, 1, 2, 4, 8, 16, 32 and 64 in a repeated five-spot are compared to experimental results for this geometrical configuration obtained using miscible fluids and conductive analogs. Predicted recovery curves for water flooding in a repeated five-spot are compared to experimental results reported in the literature. Example calculations are shown for a conventional water flood in the presence of a free gas phase and for a hot-water flood. MATHEMATICAL PRELIMINARIES We assume that the rock matrix is homogeneous with constant thickness, absolute permeability and porosity. The fluids are incompressible. We also neglect the effects of gravity.

Jan 1, 1965

By M. A. Selim, W. H. Somerton

Thermal expansions and heats of reaction of three typical sandstones were measured in the temperature range of 25°to 1,000°C. The significance of these data in subsurface heat-transfer calculations is discussed. Linear thermal expansions were measured both parallel and perpendicular to the bedding planes. Volume expansions are reported as the perpendicular expansion plus two times the parallel expansion. Expansion behavior of the sandstones was found to be controlled by the expansion characteristics of the quartz content. Differential expansion of the quartz grains and other minerals included in the sandstones caused permanent deformation of the heated samples after they were cooled to room temperature. Structural damage resulting from heating is probably an important cause of the reduction of thermal conductivity of heated samples. Measurements of heats of reaction were based on the known heat required for a- quartz inversion. Thermal reactions, which probably include deby-droxylation of clay minerals and decomposition of carbonate minerals contained within the samples, were found to require more than one-fourth of the amount of heat necessary to raise the temperature of the rock alone. In shales and limestones, the reaction heat could be substantially greater than that required from specific heat considerations alone. INTRODUCTION In earlier work, l-3 methods of measuring thermal properties were developed and thermal data for several sedimentary-rock types were reported. Data included thermal conductivity, specific heat and thermal diffusivity. In addition, calculations of heat transfer in subsurface formations may require data on the thermal-expansion characteristics of rock and the heats of reaction of mineral constituents. Thermal expansion of rocks is relatively small in magnitucle and, from the standpoint of change in bulk density, has only minor effects on heat-transfer cliaracteristics. Thermal-expansion behavior, however, may have significant effects on the structure of rocks. Differential expansion of different minerals and along different crystal-lographic axes may result in structural damage which could effectively reduce the thermal conductivity of the rock. Many mineral constituents of sedimentary rocks undergo phase or related changes when heated to sufficiently high temperatures. Although some of these reactions (such as the a-ß quartz inversion) are reversible and the heat absorbed is returned to the system upon cooling, many of the major reactions are irreversible. Through certain temperature ranges, the additional heat energy required to complete the reactions may be nearly as great as that necessary to raise the temperature of the rock in the absence of thermal reactions. If the heat requirements for these reactions are not considered, serious errors in heat-transfer calculations may result. Thermal-tzxpansion and thermal-reaction data have been obtained for three typical sandstones. Methods of measurement and results obtained are reported and discussed in the following. THERMAL EXPANSION OF SANDSTONES The linear expansions of three outcrop sandstones (Bandera, Berea and Boise) have been measured in the temptbratue range of 25° to 1,000° C. Expansion measurements were made in the directions parallel and perpendicular to the bedding. The differential expansion apparatus used is that described by Mitoff and Pask.4 The test sample is heated in an electric furnace with a temperature rise of 6°C/minute. The lengthening of the samples upon heating is compared with the small and known expansion of a fused silica rod. The change in length is transmitted to an X-Y recorder by means of a Statham transducer. A maximum error of ± 1.5 per cent has been obtained with this apparatus

By T. L. Berger, M. C. lnman, M. F. Rose

AFTER shock loading, mechanical twinning has been observed by several authors1-4 for a variety of fcc metals and alloys. It has been shown2,3 that low stacking fault energy materials deform primarily by the formation of mechanical twins and stacking faults during dynamic loading. Further, Inman, Murr, and Rose2 have shown that a correlation exists between the average hardness and the volume percent of mechanically twinned material in stainless steel. It has been observed in nickel' that mechanical twins are produced during shock deformation if the peak pressure associated with the shock wave is greater than 350kbar. It is the purpose of this work to discuss shock-hardening mechanisms in nickel and to determine to what extent the initial condition influences the final response. The "flying plate technique" was used to generate planar shock waves with peak pressures as high as 710 kbar. The details of the technique5" have been discussed elsewhere and will not be repeated here. The material used in this study was polycrystalline nickel of 99.98 pct purity and had an initial grain size of 0.18 mm. Prior to any deformation, all samples were given a stress-relief anneal of 6 hr at 600°C in a vacuum furnace. After annealing, the initial condition of the samples was changed by cold rolling individual samples from 0 pct to a maximum of 80 pct reduction in thickness. This provided samples with controlled initial histories for shock deformation. Immediately after shock loading, the specimens were slowed and quenched by suitable arrangement of buffer material and water bath. This technique minimized specimen damage and the effects of shock heating on the samples. Since deformation twins in high stacking fault energy fcc metals are only produced at high stress levels, the initial defect configuration prior to mechanical twinning should greatly influence the number of twins formed. In fact, Zukas and Fowler7 have found that prestraining in iron and steel greatly suppresses the formation of mechanical twins during the shock deformation process. Following the methods of Warren and warekois8 and Cohen and Wagner,9 X-ray diffraction profile analysis was performed to determine the incidence of stacking faults and mechanical twins in the deformed nickel. Analysis was made on samples in the annealed, cold-rolled, shock-loaded, and cold-rolled—shock-loaded states. The estimated accuracy of the peak position measurements is *0.01 deg in 2?. Therefore the accuracy in the measurement of the relative peak shifts between the (400) and (222) diffraction line profiles is on the order of *0.04 deg in 2?. None of our measurements indicated peak shifts greater than 0.04 in 28. This would correspond to a fault density of -1 pct and sets the limit on the sensitivity of our instrument. Similarly, analysis of the shift in the center of gravity of the diffraction line profiles yields the minimum measurable value of the twin fault probabilitya to be -0.002. Fig. 1 shows mechanical twinning in nickel shock-loaded to 400 kbar. In this photograph, fault intersections are clearly visible. The thickness of the individual fault lamellae was determined from the image widths of their lines of intersections by the

Jan 1, 1968

By A. M. Skov, L. F. Elkins

Inclusion of anisotropic permeability in mathematical analysis of pressure transients observed during development of the huge Spraberry field indicates a major fracture trend which is in good agreement with that observed by fluid-injection tests spread over a 12- by 17-mile area. Delineation of this trend is important in selecting a pattern of injection for the pending large-scale water flooding in this field. Determination of reservoir parameters yielding best agreement between calculated pressures and observed reservoir pressures in newly completed wells was made using an IBM 650 computer. INTRODUCTION The Spraberry field covering 400,000 acres is a tight sand of less than 1-md permeability cut by an extensive system of vertical fractures. Primary recovery dominated by capillary retention of oil in the fractured sand matrix blocks is less than 10 per cent of oil in place. Strong forces of capillary imbibition of water into the sand, coupled with water flow under dynamic pressure gradient, indicate considerable increase in oil recovery can be achieved through water flooding. Best results will occur if the pattern of water injection is selected to force the water flow across the grain of the major fracture system. Existence of an oriented vertical fracture system in the Spraberry, observed first in cores, was highlighted more recently by the 144-fold contrast in permeability along and at right angles to the major fracture trend required to match relative water breakthrough times in Humble Oil & Refining Co.'s waterflood test there. Spraberry operators since have conducted two gas-injection tracer tests for further areal confirmation of the fracture trend. Re-analysis of early reservoir pressure transients for evidence of anisotropic permeability has permitted many more local determinations of major fracture trend without resort to further field tests. This paper is limited to updating analysis of reservoir pressure transients to include anisotropic permeability as a test for orientation of the major fracture trend in the Spraberry. The reader is referred to Refs. 1 and 2 for information about general Spraberry reservoir performance and to Refs. 3 and 4 for information about significance of fracture orientation in selection of the injection-well pattern for water flooding the Spraberry. RESERVOIR PRESSURE DATA—DRIVER AREA During early development of the Spraberry Driver area, Sohio Petroleum Co. made the extra effort to measure the initial pressure in each of the 71 wells in a 5-mile-long area immediately after completion. Progressively greater reductions in pressure ranging up to 400 psi were observed throughout the six-month development period. Detailed data are presented in Ref. 1. Since the reservoir oil was undersaturated some 300 psi initially, early reservoir performance involving 55 new well pressures is subject to analysis as flow of a single compressible fluid in a porous media. Assumption of uniform permeability in all directions yielded good agreement between calculated pressures and observed pressures of these wells in the earlier study,' but subsequent, additional, mathematical development to include anisotropic permeability in the transient pressure considerations and present availability of electronic computers to perform the much more extensive arithmetical calculations now yield even better agreement. The previous analysis, assuming uniform permeability, consisted essentially of calculating pressure reduction expanding circularly around each producing well and summing these effects at the time and location of each newly completed well for comparison with the measured pressure reduction. Permeability, effective fluid and rock compressibility, and permeability X thickness were varied until the best match with measured pressures was obtained. The present analysis, assuming anisotropic

By I. Fatt, D. N. Saraf

A method is described for measuring two- and three-phase relative permeabilities in sandstones or sand packs using a nuclear magnetic resonance (NMR) technique to determine fluid saturations. Two- and three-phase relative permeabilities have been determined on Boise sandstone using the NMR technique of saturation measurement. Three-phase relative permeability to water was found to depend only on the water saturation, whereas three-phase permeability to oil depended on both the water and oil saturations. Relative permeability to gas in three-phase flow was found to depend only on the total liquid saturation. INTRODUCTION Three-phase relative permeabilities are extremely useful in calculating field performance for reservoirs being produced by simultaneous water and gas drives. Three-phase relative permeability data are also needed for analyzing solution gas-drive reservoirs which are partially depleted and are being produced by water drive. Some thermal recovery processes involve three-phase flow which require three-phase relative permeability data for predicting reservoir behavior. Unfortunately, three-phase relative permeability measurements have rarely been made. Also, because of the scarcity of three-phase data, it has not been possible to date to relate other measured rock characteristics to the relative permeabilities with a great certainty. Leverett and Lewis,l Reid2 and Snell3 have reported three-phase relative permeability data on unconsolidated sands. Leverett and Lewis used ring electrodes spaced along the length of the sand sample to measure the resistivity of the sample which was assumed to be monotonically related to brine saturation. Gas saturation was determined from pressure-volume measurements. Oil saturation was obtained by material balance on the cell containing the sand sample. This method is involved and time consuming. Another difficulty arises from the fact that the resistivity of the sand is a function not only of saturation of brine but also of the distribution and saturation history of the brine in the pore spaces. Reid2 used a gamma ray absorption technique for measuring liquid saturation. This method has the disadvantage that total liquid saturation rather than oil or brine saturation is all that can be measured and still another method is required to determine the saturations of individual components. snel13 used a neutron bombardment method which also required a separate determination of the individual component saturations. Caudle et al.4 measured three-phase relative permeability on consolidated sandstones using vacuum distillation for determining fluid saturations. Distillation after each reading makes this technique very lengthy and time consuming. Corey et a1. 5 and Naar and Wygal 6 measured three-phase relative permeability on sandstones by the capillary-pressure method. Semipermeable diaphragm assemblies were used at each end of the core specimen to keep the water phase in the core. Gravimetric methods were used to determine fluid saturations. Sarem7 recently reported an unsteady-state method for measuring three-phase relative permeability on sandstones. This method is an extension of Welge's method8 for measuring two-phase relative permeability. Although Sarem's method is simple and comparatively fast, the assumptions involved may oversimplify the problem. Sarem's assumption, that in all rocks relative permeability to each fluid will depend only on the saturation of that fluid, seems to be rather unrealistic. Neglecting capillary effects at the end of the core is also a weak assumption. Donaldson and dean9 measured three-phase relative permeability using a method similar to Sarem's.