Search Documents

By W. R. Aufricht, E. H. Koepf

The theory, measurement, and application of capillary forces in reservoir rock have been treated ex-tensively in the various petroleum publications.1,2,3,4 These forces are of prime importance in oil production, and data on capillary properties are among the most important tools used in the study of reservoir mechanics. The principal application of measured capillary pressure data by the production and reservoir engineer is in the determination of the initial distribution and mobility of fluids in the reservoir. The distribution refers to the magnitude of the oil and interstitial water saturations under capillary equilibrium conditions and to their variation throughout the reservoir as a result of changes in rock properties and nearness to free water. The magnitude of the water (or oil) saturation, as well as the porosity, is required in volumetric calculations of oil in place. In low permeability reservoirs, where long transition zones are common, the magnitude of the oil content, even below the depth interval where commercial rates of flow of oil can be obtained, must be evaluated since this oil contains energy available to the production of the remainder of the reservoir. The mobility of the water phase is of particular importance in water-flood calculations since mobile water may condemn what would otherwise be a successful water flood. In reservoirs with predominantly small pores, high immobile water contents are possible. Conversely, reservoirs with predominantly large pores may contain mobile water even though the interstitial water contents are relatively low. This presentation reviews the distribution and mobility of interstitial water contents as determined from capillary pressure data and shows typical data for various types of carbonate reservoirs. A discussion of the variations found in different types of carbonate reservoirs is given. DISCUSSION Interpretation of Laboratory Data Fig. 1 is a typical laboratory capillary Pressure curve where water saturation as Per cent of Pore space is plotted vs capillary Pressure (psig) for an air-water system. This par- ticular curve was obtained by measurement on a sample of San Andres dolomite with intergranular-type porosity. The sample studied had a permeability of 10.0 md and a porosity of 18.8 per cent. The capillary pressure-saturation relationship illustrated is an inherent property of the rock (a function of the size and distribution of the pores) and represents the ability of this water-wet rock to retain water against the pressure differential exerted by a non-wetting fluid (oil or gas) to remove it. In the laboratory a known pressure differential is exerted from a pressure source; in the reservoir the pressure differential is controlled by the density difference between the hvdrocarbons and the water. The laboratory capillary pressure values may be converted into terms of equivalent height of water column using known values of reservoir oil and water density and measured or average values of interfacial tension

Jan 1, 1958

M. AlbERtson, Chairman. Production Engineer, Shell Petroleum Corporation, Houston, Texas. R. P. McLaughlin, Associate Chairman. General Manager, Burnham Exploration Co., Los Angeles, Calif. Eugene A. Stephenson, Secretary-Treasurer. Professor of Petroleum Engineermg, Missouri School of Mines and Metallurgy, Rolla, Mo. Ezecutive Commillee M. Albertson, Chairman. C. E. Beecher, Empire Oil & Refining Co., Bartlesville, Okla. L. L. Foley, Petroleum Engineer and Geologist, Tulsa, Okla. R. B. Kelly, Division Manager, Pure Oil Co., Fort Worth, Texas. Hallan N. Marsh, Production Engineer, General Petroleum Corporation of California, Los Angeles, Calif. R. E. Somers, Geologist, The Gulf Companies, Pittsburgh, Pa. Production Engineering T. V. Moore, Vice-chairman. Charge Production Research Division, Technical Service and Development Dept., Humble Oil & Refining Co., Houston, Texas. Joseph Chalmers. De Bataafsche Pet. Maatschappij, The Hague, Netherlands. M. L. Haider, Petroleum Engineer, Carter Oil Co., Tulsa, Okla. Emile HUGUENIN, Deputy State Oil and Gas Supervisor, Los Angeles, Calif. Carl E. Reistle, Jr., Petroleum Engineering Dept. Humble Oil & Refining CO., Houston. Texas. Produclion J. Terry DUCE, Vice-chairman. Geologist, The Texas Company, New York, N. Y. Basil B. Zavoico, Associate Vice-chairman. Consulting Geologist, Houston, Texas. Engineering Research R. L. Huntington, Vice-chairman. Associate Professor, Petroleum Engineering, University of Oklahoma, Norman Okla. Economics Sidney A. Swensrud, Vice-chairman. Assistant to the President, Standard Oil CO., Cleveland, Ohio. Refinery Engineering Walter Miller, Vice-chairman. Vice-president, Continental Oil Co., Ponca City, Okla.

Jan 1, 1937

By Ernest E. Thurlow

The western coal industry, long overshadowed by the petroleum and metallic minerals industries, has been revitalized by several new coal mining operations, with the promise that this is only the beginning of a tremendous surge in production of coal that will carry into the next century. The content of this paper is confined to the coal fields and recent developments in eastern Montana, northeastern Wyoming and western North Dakota. The area is commonly referred to as the Powder River-Fort Union region although the Powder River Basin and the Fort Union formation are separate entities and are not coextensive. Coincidentally or otherwise, the extent of the subbituminous part of the coal region coincides approximately with the Powder River Basin. Environmental considerations together with the low sulfur content of much of the western coal have triggered demand since this coal could not compete at present with eastern and midwestern coal solely on a delivered-cost-per-million-btu basis. However, increasing costs of underground mining and the comparatively thin seams of the eastern and midwestern deposits may soon bring a balance to these costs regardless of the sulfur content. Another factor which has made western coal competitive at this time is the use of the unit-trains to significantly reduce the cost of transportation. Fast loading and unloading facilities and trains of 100 cars or more in turn-around service are the keys to success in this effort.

Jan 5, 1974

By L. W. Casteel

St. Joseph Lead Company's Southeast Missouri mines have been successfully converted to the use of ammonium nitrate and fuel oil as a blasting agent. This is significant to the company because it means a savings of a quarter of a million dollars per year. It is significant to the mining industry because of the unusual conditions involved and the methods and materials used to overcome them. St. Joe has been mining the low grade lead deposits of southeast Missouri since 1864. This ore occurs as disseminated galena in widely scattered pockets within the Bonneterre dolomite. The daily production of 20,000 tons is obtained from several hundred small, room and pillar stopes covering an area 10 miles long and five miles wide. Distribution of explosives and other materials is a real problem. Drilling in the stopes is done with air legs, column- mounted machines or jumbos, depending on the size and location of the stope. Drill hole sizes vary from 1 3/8-in. diam down to 1 1/8-in. Each driller charges his own holes and shoots them at the end of the shift. Numerous shafts provide natural ventilation which is supplemented in isolated areas by churn drill holes equipped with blowers. Water, in most stopes, is not a serious problem, although the mines do produce about 20,000 gpm.

Jan 5, 1962

The resolution drafted by J. W. Richards, J. V. N. Dorr, and Allen H. Rogers on the death of J. E. Johnson, Jr., and adopted by the Board of Directors, is as follows: "The Boards of Directors of the American Institute of Mining and Metallurgical Engineers wishes to place on record its great personal loss in the sudden death of Mr. J. E. Johnson, Jr., our fellow Director of the activities of the Institute, a member our Executive Committee, Chairman of the Institute's Committee on Admissions, and Vice-chairman of its Iron and Steel Committee. "A mere recital of these positions of trust and responsibility, which he filled with conscientious efficiency and marked ability, shows how greatly the Institute has lost by his death. No one else of the Institute's officers could have been so illy spared. We, his colleagues, take up his tasks with a sense of our inability to discharge them as well as he, and this with a doubled sense of our loss. "Others have written and spoken of the blow which his death has been to the furtherance of the metallurgy of iron and steel. His pre-eminence in this line is one of the facts of which the Institute is justly proud. Here his personality, originality, courageousness in advocating new principles, clearness of writing and forcible presentation of facts, are irreplaceable. His two great treatises on the blast furnace are his literary and scientific memorial.

Jan 7, 1919

By R. G. Hawthorne

This report is concerned with fluid displacement in porous media, in those cases where viscous and gravitational forces control the displacement. Such a system would usually be found in a sand body of large physical dimensions such as an oil reservoir, although it is possible to create such a system in the laboratory. It is shown that the position of the fluid interface can be predicted by numerical calculations using a basic idea presented by Dietz. Fluid flow is considered in a vertical plane in a homogeneous, porous medium of sufficient thickness that the capillary transition zone is small in comparison with the total reservoir. A theory developed by Dietz' is used to make numerical calculations of the position of the fluid interface. The results for several conditions are compared with scaled model experiments. The results show that, for gas drive in a reservoir of steep dip, a relatively low flow rate can displace large volumes of oil before gas breakthrough. On the other hand, water injection at favorable mobility ratio and low dip may show best performance at high rates. Water tends to underride the oil and, given sufficient rime, will break through without much oil displacement. For certain conditions, which include relatively low flow rate, the interface is a straight line and its behavior is simple to calculate. At higher flow rates, the interface is unstable, and a numerical solution was programed for an automatic computer. In general, good agreement is shown between the fluid model and the computed results so long as gravitational forces have control. For a water drive at very unfavorable mobility ratio, many small water fingers appear. These viscous fingers are not controlled by the relatively small gravitational forces. When viscous fingering becomes the controlling factor, the mathematical model is oversimplified, and results do not check the fluid flow model. INTRODUCTION Present methods of reservoir analysis depend upon certain simplifying assumptions to obtain mathematical descriptions of practical use. Material-balance methods (Muskat2 or Tarner2) assume uniform fluid saturations in the entire reservoir, or in a few subdivisions of the reservoir. An unsteady-state flow calculation by West, er at considered pressure and saturation changes in flow to a well during solution gas drive, and neglected gravity effects. Results showed only a 4 per cent difference for ultimate oil recovery by the Muskat method, even though the case chosen for study was one in which unsteady-state effects should be high. The Buckley-Leverett5 method commonly assumes a one-dimensional flow system. It is applicable at high flow rates where viscous forces predominate over gravity forces. Simultaneous, parallel flow of the two fluids is assumed, and the concept of a fluid interface is not introduced. Permeabilities to each fluid for a given saturation must be known. The method is not applicable for a two-dimensional system where cross flow becomes possible. Less well known is the displacement equation derived by Dietz. This method is designed for two-dimensional flow systems and assumes a definable fluid interface within the porous medium. Dietz showed that, for a range of low flow rates, the interface would be stable. straight and at an angle of inclination which could be simply calculated. At a certain critical flow rate, the calculated interface tilt would equal the formation dip. For higher flow rates, a finger of displacing fluid would invade the displaced fluid. Dietz indicated that his method applied only to macroscopic reservoir behavior, while the Buckley-Leverett method applied to the small transition zone at the fluid interface. The examples worked out in this report are based on the fluid-displacement theory of Dietz. It is shown that the Dietz theory may be used to derive equations analogous to the Buckley-Leverett equations. In contrast to the Buckley-Leverett method, flow is considered in a plane rather than being limited to a line. Rather than a frontal advance, the movement of a fluid interface is followed. For flow rates substantially exceeding the critical rate and for high viscosity ratio, many fingers of invading fluid occur-—rather than the single finger assumed by Dietz. On the other hand, so long as some gravitational influence remains, the flow is not entirely parallel to the bedding planes as assumed by Buckley and Leverett; therefore, both methods fail to give an adequate descrip-

By H. H. Rachford

Numerical solutions of immiscible flow problems in which dispersive effects of capillarity are dominated by convection require excessively fine grid spacing with attendant high computing costs. The use of coarser spacing reduces cost but often produces oscillation or undue dispersion associated with displacement fronts, A numerical formulation is proposed here which should be applicable to two-dimensional flow problems. it is in part analogous to an approach previously tested for miscible systems. The convective transport is approximated using a change of variables to yield a coordinate system moving approximately with the local characteristic velocity. The capillarity-induced dispersive terms in the differential system describing the process are approximated with respect to a fixed coordinate system by the usual implicit formulation. One-dimensional tests of the procedure yielded results in which the saturation profiles tended smoothly to the zero-capillary pressure solution as the ratio of viscous to capillary forces was successively increased in a sequence of calculations. This contrasted favorably with solutions by other numerical procedures which would require attendant grid refinements to approach the zero capillary pressure results. INTRODUCTION Numerical solution of displacement problems has until recently relied on applying methods developed primarily for transient heat-flow problems. Such problems are classified as parabolic in type, and where the heat transport is purely by diffusion their solutions are characterized by a high degree of smoothness. It is not surprising, therefore, that for approximating these solutions available finite difference methods are quite adequate. In flow problems the transport is partly by diffusion, partly by convection or flow. Although the problem remains of parabolic type because the dispersive effects of capillary forces or diffusion play some role in every displacement, at high flow rates the problem is dominated by convection, and solutions tend toward those of equations of the hyperbolic type. Solutions of hyperbolic problems are characterized by the translation of fronts, or discontinuities, that may progressively increase in sharpness. Numerical methods for treating parabolic problems become less and less satisfactory as displacement rates increase and the role of dispersion due to concentration or capillary pressure gradients becomes small relative to transport due to flow. In computation the difficulty manifests itself as an error associated with the grid size chosen. 1-6 In summary, if the heat-flow type approximations are to include the terms arising due to convection, one of several choices may be made: (1) an upstream (to the direction of flow) approximation for the convection terms may be used; (2) a centered-in-distance (CID) approximation may be used; or (3) a recently developed approximation based on the theory of oscillation matrices may be chosen.6 The last appears to have significant promise for one-dimensional flow problems; its extendibility to two or three dimensions is an open question. In either of the first two approaches, a suitably small ratio of v&/D must be maintained, where v is the velocity, & is the grid spacing and D the effective dispersivity in the direction of flow. In the first choice, the approximation of the convective part is only first-order correct and errors introduced appear as a numerically induced dispersivity of magnitude proportional to v?x. In the CID choice, the approximation can be second-order correct, but the difference formulation fails to satisfy the maximum principle unless a condition on v?x/D is met. Practically, this means that for high flow rates oscillatory solutions may result in the neighborhood of a front unless exceedingly small grid intervals are taken. While the procedure proposed by Stone and Brian4 permits a less severe limitation to be placed on this ratio, ultimately the flow rates increase relative to the dispersivity the oscillation obtains. Further, extensions of their approach to higher dimensional systems may be attended by considerable

Jan 1, 1967

By J. L. Huitt, S. R. Darin

The use of a partial monolayer of propping agent to obtain a high flow capacity for a hydraulically induced fracture is discussed. From the results of laboratory work it was shown that a modified form of the Kozeny-Carman relation could be used to describe the flow in the partial monolayer propped fracture. With equations presented in the paper, the density pattern of the propping agent (number of particles per unit of fracture surface) that results in the maximum flow capacity for the fracture can be determined. The maxitnum flow capacity obtained with a partial monolayer is often an order of magnitude greater than the flow capacity obtained in greater width fractures containing multilayers of the propping agent. INTRODUCTION One of the predominant factors controlling the success of a hydraulic fracturing operation is the propping of the fracture. The trend in hydraulic fracturing recently has been an increase in the ratio of propping agent-to-fluid. The primary purpose of this increased ratio is to sustain a propped fracture of greater width. It was shown in a recent study' that for some formations low concentrations of propping agent would result in "closing" or "healing" of the fracture. In some cases, in an effort to insure sufficient propping agent concentration, the fracturing operation is designed to obtain a pack of the propping agent in the fracture. The placing of a pack of propping agent in the fracture provides a fracture of maximum width; however, the width alone does not control the flow capacity of the fracture. The flow capacity is dependent on the permeability of the fracture, as well as the fracture width. Thus, increasing the permeability to obtain larger flow capacities is as important as obtaining a greater width fracture. By increasing the concentration and particle diameter of the propping agent, fractures of sufficient flow capacities are obtainable for effective well stimulation in a number of formations. However, in some of the recovery techniques, the fracture capacity obtained by commonly used practices in fracturing is inadequate. One such recovery technique is gravity drainage via a horizontal fracture placed near the base of the productive zone. For this recovery technique, the fracture capacity desired is usually one to two orders of magnitude greater than the fracture capacity normally obtained in fractures propped with the size of sand commonly used (20-40 mesh). By the use of particles of larger diameter, the permeability of a pack2 is increased; however, even a multilayer pack results in a fracture flow capacity lower than that desired for some formations. For example, a fracture of 1/4-in. width, packed with 8-12 mesh sand of about 2,000 darcies permeability, would result in a fracture flow capacity (permeability X fracture width) of only 40,000 md-ft. If the fracture flow capacity needed is an order of magnitude greater than that cited in the example, an approach different from the formation of a pack in the fracture is needed. One approach is a partial monolayer of large-diameter propping agent. In working with this type of fracture system, it was found that the flow problems were unique to the system; therefore, the study discussed in this paper was made. THEORY AND RELATED STUDIES Previous studies, although related to that described in this paper, have not pertained directly to the flow in a fracture propped with a partial monolayer of propping agent. However, studies such as flow between parallel plates,' flow in packed columns,' fracture flow capacity' and flow in simulated fractures reduce the scope of this flow problem considerably. A fracture propped with a partial monolayer of dropping agent may be considered to be between two extremes — between an open fracture and a bed of particles packed In a conduit of fracture geometry. Thus, it seems appropriate to review the extreme cases.

By P. F. Gnirk, J. B. Cheatham

Single bit-tooth penetration experiments under static load were conducted on six rocks at confining pressures of O to 5,000 psi using sharp wedge-shaped teeth with included angles ranging from 30 to 120° In general, the force-displacement curves for all rocks exhibit an increasingly nonlinear and discontinuous behavior with decreasing confining pressure. The confining pressure at which a rock exhibits a macroscopic transition from predominantly ductile to predominantly brittle behavior during penetration vmies from about 500 to 1,000 psi for the limestones to greater than 5,000 psi for dolomite. The correlation between calculated values of force per unit penetration based on plasticity theory and experimental values is quite encouraging, even at confining pressures as low as 1,000 psi A qualitative correlation between volume of fragmented rock per unit energy input for a single bit-tooth and drilling rate for microbits appears to exist over a confining pressure range of 0 to 5,000 psi. INTRODUCTION Laboratory experiments utilizing a small-scale drilling apparatus have demonstrated that penetration rates are reduced considerably as a result of increasing the confining pressure from atmospheric to a few thousand psi?-3 This undesirable situation can, in general, be attributed to a combination of decreased efficiency of chip removal at the bottom of the borehole, increased rock-failure strength, and a possible change in the mechanism of chip generation and rock fragmentation with increasing confining pressure. To more fully understand the principles underlying the last circumstance, it is the purpose of this investigation to experimentally study the mechanism of single bit-tooth penetration into dry rock at low confining pressures and, in particular, to establish the confining pressure at which the penetration mechanism may undergo a brittle to ductile transition for various rock types commonly encountered in drilling. Confining pressure as used here refers to the differential pressure between the borehole fluid pressure and the formation-pore fluid pressure. EXPERIMENTAL PROCEDURE Using an experimental apparatus previously described: a single, sharp wedge-shaped tool was forced under a "statically" applied load into an effectively semi-in finite dry rock sample subjected to a prescribed confining pressure. To prevent the invasion of the con fining-pressure fluid into the pores of the rock sample during penetration, the exposed surface of the rock was jacketed with a layer of silicon putty.* Electrical instrumentation incorporated into the apparatus yielded a graphical plot of force on the tool as a function of penetration or displacement of the tool into the rock during an experiment. During the course of the experimentation the following conditions were maintained constant: (1) pore pressure — atmospheric (i.e., the rock was dry); (2) temperature — 75F; (3) rate of loading — essentially static (approximately 0.002 in./sec); (4) bit tooth — a sharp wedge-shaped tool loaded normal to the rock surface; (5) rock surface — smooth and flat; (6) drilling fluid — hydraulic oil; and (7) maximum depth of penetration — approximateiy 0.1 in. In addition, each experiment was performed on a different rock sample so the rock surface is free of a layer of cuttings and of any previous indentation craters. The influence of the corners of a borehole was neglected, since each rock sample was cemented inta a section of aLuminum tubing to simulate a semi-infinite body. Indentation experiments were made on a sandstone, two limestones, a marble, a dolomite and a schist with a 60' bit-tooth over a range of confining pressures from atmospheric to 5,000 psi and at a constant confining pressure of 1,000 psi for a variety of bit-teeth with angles ranging from 30 to

Jan 1, 1966

By P. A. Robbins

FAR NORTH, on the bare Arctic tundra, 11 mi. above the mouth of the Keewalik River where the latter discharges into Kotzebue Sound. several ., Eskimos garbed in parkies and muck lucks mingle with a small group of similarly clad white men. All of them shield their eyes from the white glare of the arc as a gasoline-driven arc-welding set rains molecules of metal on the break in a bull wheel, the large gear wheel which swings the boom and cab of a dragline shovel.

Jan 1, 1937

By P. Chiotti, J. T. Mason

Thermal, X-ray, metallographic, and resistivity data were obtained to establish the phase diagram for the Yb-Zn system. Six intermetallic compounds and five eutectics have been identified. The compounds position at 695°C. No evidence for the compound YbZn,, was observed. Whether this compound is stabilized or inhibited by small amounts of impurities is uncertain. SEVERAL aspects of the Yb-Zn system have been reported in the literature. Much of this work deals with the crystallographic identification of some of the compounds. The YbZn compound has been reported by Iandelli and palenzonal and the YbZn, structure and thermal data for the composition range between YbZn and Yb3Znll have been reported by Michel et al. A compound of uncertain stoichiometry in the range 60 to 68 wt pct Zn with a cubic structure has been reported by Kuz'ma et aL4 The YbzZn17, YbZn,,, and YbZn13 compounds have been reported by both Kuz'ma et a1 and Iandelli and Palenzona.' YbnZn17 has also been identified by Laube. This investigation was undertaken to establish the phase diagram of this system. Ytterbium in many respects behaves as an alkaline earth rather than a rare earth metal6 and it is of interest to note the resemblance between the Yb-Zn system, the Ca-Zn,' and typical rare earth-zinc systems.g'10 MATERIALS AND EXPERIMENTAL PROCEDURES The materials used in this investigation were Ames Laboratory ytterbium and Bunker-Hill slab zinc, 99.99 pct pure. Analyses of the ytterbium showed 38 ppm Ta, 82 ppm 0, and 15 ppm N with Ca, Cu, C, Mg, Al, Si, Cr, Fe, Ni, Sc, Y, Er, Tm, and Lu all less than 50 ppm. The ytterbium was dehydrogenated by heating in tantalum under vacuum for 8 hr at 900°C. The metal was kept under argon except when being cut with shears and when being weighed. Tantalum was found to be a suitable container for alloys in the composition range between ytterbium and Yb2Zn17. With higher zinc concentrations the tantalum was penetrated by the alloy and failed after several days at temperatures of 600" to 700°C. Approximately 20 g of zinc and ytterbium metal were sealed under argon in *-in.-diam tantalum cruci- bles equipped with thermocouple wells. The crucibles were in turn sealed in stainless-steel jackets. The alloys were equilibrated in an oscillating furnace and given various heat treatments in a soaking furnace. The steel jackets were removed and differential thermal analysis data were obtained under an atmosphere of argon. The thermal analysis apparatus was calibrated against the melting points of pure zinc, aluminum, silver, and copper and found to give thermal arrests within 1°C of the accepted melting temperatures for these elements. The thermal analysis apparatus has been described in an earlier paper.g Alloys were further heat-treated as necessary and subjected to metallographic examination and X-ray diffraction analysis. Most alloys were etched with Nital (0.5 to 10 pct nitric acid in alcohol). EXPERIMENTAL RESULTS The results of the thermal analysis are summarized by the points on the phase diagram, Fig. 1. Six compounds and five eutectics were observed. The eutectic temperature and compositions in weight percent zinc are 491°C and 12 pct, 645°C and 32 pct, 687°C and 54 pct, 742°C and 67 pct, and 750°C and 77 pct. The compounds YbZn, YbZn2, Yb2Zng, Yb2Zn17, and YbZnll melt congruently at 650°, 751°, 75Z°, 754", and 755"C, respectively. The compound Yb3Znll undergoes a peritectic decomposition at 695°C. The transformation and melting temperatures for the ytterbium were observed to be 797" and 817°C. re-

Jan 1, 1969

By J. Geertsma

The relationship between porosity and the speed of propagation of acoustic waves in fluid-saturated porous rocks as measured by the Sonic log and by ultrasonic tecbniques is analyzed. Biot's continuum theory 1 is used to explain the difference in acoustic wave propagation between a dry and a liquid-saturated porous material. The porosity is a variable in this theory. However, the acoustic wave propagation in the dry rock depends too on porosity, and this dependence is not predicted by the theory. Frequently in dry sandstones, a nearly linear relationship between reciprocal acoustic wave velocity and porosity is observed in the low-porosity range. The physics behind this behavior is outlined. An empirical relationship of the form, 1/v = A + F, applies accordingly for many porous dry rocks, provided the porosity is the only variable. The presence of a liquid in the pores changes the value of B, and this change is found to be in agreement with the Biot theory. The time-average relation introduced some years ago2 results in an equation of the same type — 1/V = F/Vf + (1 - F)/Vr — but is not based on a sound physical picture. Still, this relation sometimes predicts approximately correct A and B values. Carbonate rocks with their complicated pore structures sometimes show a different relationship between wave velocity and porosity, unfavorable for log interpretation. Examples are presented. The simultaneous presence of calcite, dolomite and anhydrite, with their different grain densities and matrix compressibilities, complicates acoustic-log interpretation in carbonate rocks still further. Other complicating effects of acoustic-log interpretation are discussed. They are related to the influence of shale streaks and natural fractures on the average wave velocity observed by the logging tool and to the effect of adsorption phenomena on wave propagation in unstressed rocks particularly in sandstones. INTRODUCTION The velocity of propagation of sound waves in porous sediments as a function of depth is a quantity frequently measured. Velocity loggers of various d-sign may be used, for which velocity is defined as the. distance between a wave generator and a wave detector (one-receiver system) or the distance between two wave detectors (two-receiver system), divided by the shortest time required for a vibrational pulse to cross this distance. Velocity loggers ale used to assist in seismic prospecting, to differentiate between the various types of sedimentary rock layers, and also to determine rock porosity and pore fluid content. This paper is especially concerned with the relation between velocity and porosity and the effect of the pore fluid. Most velocity-log interpretation in terms of porosity is based on a formula advocated by Wyllie, et al,2 suggested earlier by Hughes and Jones,3 and known as the "time-average relation". This formula app1ies for a model consisting of a layered system of parallel alternating slices of a solid and a liquid, crossed by the wave path perpendicular to the interfacrbs. This must be an unsuitable model for the derivation of wave propagation properties of liquid-saturated porous media. It suggests that only rock matrix and fluid properties influence the wave velocity. The surprising fact, nevertheless, is that applications of this formula to clean sandstones under specific conditions (sufficient effective stress) results in porosity predictions which are close to reality. Other authors, such as Gassmann,4 White and Sengbush,5 Brandt,6 and Hicks and Berry,7 make use of more realistic models to arrive at a wave velocity formula for porous sediments. Their models consist of various liquid-saturated packings of frictionless spherical grains and are, therefore, also highly idealized. A number of experimental results obtained in the laboratory can be better explained with these model theories than by the time-average relation. A study of a general character, not depending a priori on a special model of the porous structure, was still lacking.

By J. Szekely, G. P. Martins

A formulation is given and a general analytical solution is presented for the problem of spherical phase growth in multicomponent systems. The results are valid for describing the rate of bubble growth in the absence of translational motion and, therefore, should be relevant to heterogeneous nucleation processes. Typical worked examples are given, which include bubble growth due to the simultaneous diffusion of hydrogen and nitrogen in molten steel, and the growth of CO bubbles in steel, due to the reaction: C + 0 = CO. The calculations predict that the presence of additional components would enhance the growth rate and also increase the absolute rate at which the individual components are transferred into the bubble. These findings should be relevant to degassing processes. THE problem of spherical bubble growth in the absence of translational motion is appropriate to the description of the heterogeneous nucleation and subsequent growth of gas bubbles while attached to a solid surface. Such problems are of considerable practical relevance to various vacuum degassing and steelmaking processes where bubble nucleation and subsequent growth play an important part in the overall reaction mechanism. A definitive formulation of spherical phase growth has been given in the now classical paper by Scriven1 and the analogous problem for solid phase growth has been discussed by Horvay,2 furthermore a thorough review of the relevant literature has also become available.3 All this previous work was concerned with one component or binary systems. The purpose of this paper is to show that the technique developed by Scriven may be generalized to alow the consideration of multicomponent systems. The results of this generalization would be particularly relevant, as gas bubble nucleation in the metallurgical systems previously mentioned frequently involves the simultaneous diffusion of several components. In the following a general formulation is given and typical worked examples are presented. FORMULATION Consider the growth of a spherical gas bubble due to diffusion in a quiescent liquid of infinite extent, as sketched in Fig. 1. The formulation of such problems requires the statement of the following conservation equations relating to the liquid phase: Equation of continuity Equation of motion Conservation of the diffusing species—Diffusion Equation. The following assumptions are made: a) The fluid is incompressible, quiescent, and of infinite extent b) The bubble is stationary, i.e., it undergoes no translational motion c) Spherical symmetry is observed throughout d) The mass diffusivity is constant e) Chemical equilibrium exists at the surface of the bubble, and any reactions occurring there are fast compared to the diffusion process f) Inertia and surface tension effects are neglected. Equations derived on the basis of these assumptions by previous investigators were found to provide a reasonable representation of bubble growth while attached to a surface in organic and aqueous systems4,5 For gas bubble growth in molten metal systems assumptions a) to d) chould hold reasonably well. The exclusion of chemical kinetics as a possible rate limiting step may not be universally valid, but should hold for many systems operating at high temperatures. The effect of hydrodynamic inertia may be neglected for moderate growth rates; finally, the assumption neglecting the effect of interfacial tension on the bubble growth was made for the sake of mathematical convenience; its inclusion would have required the use of quite cumbersome numerical techniques. While surface tension effects could be appreciable on the initial values of the overall growth rate, especially at low absolute pressures and in case of small bubbles, they are not thought to modify the relative effects present in multicomponent systems, the study of which is the objective of this paper. Within the framework of the above assumptions

Jan 1, 1970

Albert Sauveur, Cambridge, Mass. (written discussion).— In his introductory paragraph, the author states that he has come to bury Beilby, not to praise him. The paper, however, gives the impression that while he may have dug the grave, Beilby is still at large. I am not one of those whom the author calls enthusiastic supporters of Beilby's amorphous theory and in support of this claim quote from a book I am writing: "Attractive as the amorphous theory is, and notwithstanding the ability with which it has been presented and the renown of its sponsors and supporters, it should not be accepted without reservations as a demonstrated fact.... So long as it explains more satisfactorily than other theories a greater number of phenomena in which it may play a part, we should use it in preference to others as the best available. Let us avoid the danger, however, into which too many fall, of accepting as indisputable truths, statements and theories emanating from those whom we hold in great esteem. It is not what those men themselves would wish us to do. I am entitled, I think, to occupy a position in the middle ground. Neither should Howe be classified with the enthusiastic supporters for he has written that the theory should be classed rather among the precious working theories than among those firmly established," a very wise and prudent statement. Jeffries, on the other hand, I would not hesitate to describe as an enthusiastic supporter. The author reproaches the advocates of the amorphous theory with having taken whatever properties they needed and with having assigned them arbitrarily to the hypothetical amorphous phase. He then applics their method to small crystals. If the existence of an amorphous phase is assumed, its greater tenacity and greater hardness is satisfactorily explained in the words of Howc: '(This greater strength is natural, because the amorphous state lacks the cleavages and slip planes which weaken crystalline substances." Can the greater strength and hardness claimed for small crystals be as satisfactorily explained? Nor do I think that the cold deformation of metals through successive slips within the crystalline grains can he accounted for as acceptably by the "small crystals theory."

Jan 1, 1923

By Will L. Clark

He who enters a mining claim within an Indian reservation of the United States of America acquires no rights thereby, because of the fact that the lands within such Indian reservation are not a part of the public domain, and cannot be located or acquired as mineral land. However, where such an entry is authorized as to a particular reservation by an Act of Congress, such claims as may be ontered under such authorization are thereby segregated from the reservation, and the Indian title is oxtinguished. Under the Choctaw Constitution, any citizen of that nation, in other words any Choctaw Indian, who discovered a coal mine, acquired an exclusive right to all coal within a radius of 1 mile. The laws of the Chickasaw nation provided for the formation of corporations to develop coal and other mines, with authority to contract with others to develop and work the mines. Under these laws, leases were made; but Congress subsequently abrogated all existing leases, prohibited all persons from receiving royaltics from such mines, and provided that all coal within the nation should remain the common property of the tribes. , Such leases are now expressly prohibited by an Act of Congress (32 United States Statutes at Large, 655). All leases of mineral lands upon Indian reservations must now be made under regulations promulgated by the Secretary of the Interior, and the royalties must' be paid into the United States Treasury for the benefit of the tribes. In Canada the Indians have no rights to the royal mines and minerals, and the government can make no stipulation with the Indians which would affect the rights of the Province in which the lands are situated. Therefore, while it is true that under our present law no rights can be acquired to a mining claim within an Indian reservation, it seems to be fully within the authority of Congress to make laws, giving such rights to mineral locators as it may deem necessary or expedient. There is a demand for a general law that will permit mining upon Indian reservations, and I am informed that most of the members of both houses of Congress from the Western States are in favor of Congress pass-

Jan 1, 1915

Albert Sauveur, Cambridge, Mass. (written discussion).— In his introductory paragraph, the author states that he has come to bury Beilby, not to praise him. The paper, however, gives the impression that while he may have dug the grave, Beilby is still at large. I am not one of those whom the author calls enthusiastic supporters of Beilby's amorphous theory and in support of this claim quote from a book I am writing: "Attractive as the amorphous theory is, and notwithstanding the ability with which it has been presented and the renown of its sponsors and supporters, it should not be accepted without reservations as a demonstrated fact.... So long as it explains more satisfactorily than other theories a greater number of phenomena in which it may play a part, we should use it in preference to others as the best available. Let us avoid the danger, however, into which too many fall, of accepting as indisputable truths, statements and theories emanating from those whom we hold in great esteem. It is not what those men themselves would wish us to do. I am entitled, I think, to occupy a position in the middle ground. Neither should Howe be classified with the enthusiastic supporters for he has written that the theory should be classed rather among the precious working theories than among those firmly established," a very wise and prudent statement. Jeffries, on the other hand, I would not hesitate to describe as an enthusiastic supporter. The author reproaches the advocates of the amorphous theory with having taken whatever properties they needed and with having assigned them arbitrarily to the hypothetical amorphous phase. He then applics their method to small crystals. If the existence of an amorphous phase is assumed, its greater tenacity and greater hardness is satisfactorily explained in the words of Howc: '(This greater strength is natural, because the amorphous state lacks the cleavages and slip planes which weaken crystalline substances." Can the greater strength and hardness claimed for small crystals be as satisfactorily explained? Nor do I think that the cold deformation of metals through successive slips within the crystalline grains can he accounted for as acceptably by the "small crystals theory."

Jan 1, 1923

By J. L. Lytton, T. E. Tietz

The recovery of flour stress of recrylstallized are-cast molyhdenum after 10 pet prestrain at 150°C was studied for recovery temperatures of 725o 801o and 901oC. The yield phenomenon observed after recovery was such that the yield drop diminished and disappeared as recovery proceeded. A fractional flow-stress recovery parameter fr, was developed and it was shown that, independent of recovery ternperatzwe, equal f, values of about 0.16 produced superposition of stress-strain curves except for the initial yield phenomenon. The activation energy for low degrees oj recovery was found to he 32,600 cal per mole, in close agreement with a reported valzie for difjzdsion of carbon in molybdenum. For values of fr greater than about 0.15 where recovery was thought to predominate over strain aging, the activation energy varied from 43,600 to 82,500 cal per mole and increased with increasing recovery and recovery temperature. A 10 pet pre-strain shifted the ductile-brittle transition temperature, as indicated by tensile elongation, from 10°C. to about -22"C. A 100-hr recovery anneal at 801oC was shown to shift the transition temperature to about —5oC, the major shift occurring during the first hour of recovery. ThE recovery of cold-worked metals is accompanied by changes in the arrangement of the dislocations generated during cold work. In recent years, the use of transmission electron microscopy has allowed direct observation of some of the dislocation configurations formed, such as "tangles", networks, cells, subboundaries, and loops. These observations have recently been reviewed by Keh and Weissmann.' Accordingly, some understanding has been gained concerning the nature of the dislocation rearrangements involved. The relationship between recovery and the decrease of mechanical strength, however, has only recently been explored. Previous efforts of the authors were directed toward investigating the decrease of tensile flow stress of aluminum and Al-1 pet Mg during recovery2 and the effects of concurrent elastic and plastic strain on the flow-stress recovery of aluminum.3 In addition, an investigation of the effect of four different solute additions on the tensile flow-stress recovery of aluminum has been made.4 Many important metals, however, have the bee crystal structure, e.g., iron, molybdenum, tungsten, vanadium, columbium, tantalum, and chromium. These metals might be expected to behave differently during recovery from those with other crystal structures for two reasons. First, the previous and present investigations of aluminum and its alloys have shown that tensile flow-stress recovery is strongly affected by foreign atoms in the lattice. In the bee metals, interstitial impurities are usually present and might be expected to have various effects upon both flow behavior and recovery. For example, heterogeneous yielding might occur during prestraining. The diffusion of these interstitials can also relock dislocations during recovery so that a yield point might be obtained upon restraining (strain aging). Second, some of the bee metals (e.g., iron, molybdenum, and tungsten) undergo the familiar structure-sensitive ductility transition and exhibit brittle fracture at suitably low temperatures. It is of considerable interest to determine the effects of recovery on the transition temperature of such metals. The purpose of the present investigation was to extend our knowledge of flow-stress recovery to the bee metals and to study the effects of recovery on the ductile-brittle transition temperature. Unalloyed recrystallized are-cast molybdenum was chosen as a test material, since it exhibits a well-defined ductility transition near room temperature. As in the previous studies, the degree of recovery was measured in terms of the decrease of tensile flow stress. Prestraining and restraining were carried out at 150°C to avoid the possibility of brittle fracture. EXPERIMENTAL PROCEDURE The are-cast molybdenum sheet used in this investigation was obtained from the Climax Molybdenum Co. of Michigan. It was 0.060 in. thick by 10

Jan 1, 1964

By P. M. Blair

This paper presents numerical solutions of the equations describing the imbibition of water and the countercurrent flow of oil in porous rocks. The imbibition process is of practical importance in recovering oil from heterogeneous formations and has been studied principally by experimental means. Calculations were made for imbibition of water into both linear and radial systems. Imbibition in the linear systems was allowed to take place through one open, or permeable, face of the porous medium studied. ln the radial system, water was imbibed inward from the outer radius. The effects on rate of imbibition of varying the capillary pressure and relative permeability curves, oil viscosity and the initial water saturation were computed. For each case studied, the rate of water imbibition and the saturation and pressure profiles were calculated as functions of time. The results of these calculations indicate that, for the porous medium studied, the time required to imbibe a fixed volume of water of a certain viscosity is approximately proportional to the square root of the viscosity of the reservoir oil whenever the oil viscosity is greater than the water viscosity. Results are also presented illustrating the effects on rate of imbibition of the other variables studied. INTRODUCTION The process of imbibition, or spontaneous flow of fluids in porous media under the influence of capillary pressure gradients, occurs wherever there exist in permeable rock capillary pressure gradients which are not exactly balanced by opposing pressure gradients (such as those resulting from the influence of gravity). The importance of such capillary movement in the displacement of oil by water or gas was recognized in early investigations and described by Leverett, Lewis and True in 1942.1 Methods advanced by these authors for studying the process using dynamically scaled models were rendered more general and flexible by the research of later workers.2"4 The influence of capillary forces in laboratory water floods has also been discussed by several authors.6 While imbibition plays a very important role in the recovery of oil from normal reservoirs, Brownscombe and Dyes7 pointed out that imbibition might be the dominant displacement process in water flooding reservoirs characterized by drastic variations in permeability, such as in fractured-matrix reservoirs. In water-wet, fractured-matrix reservoirs, water will be imbibed from fractures into the matrix with a countercurrent expulsion of oil into the fractures. If the imbibition occurs at a sufficiently rapid rate, a very successful water flood can result; if the imbibition proceeds slowly the project might not be economically attractive. Scaled-model studies have demonstrated the vital importance of imbibition in secondary recovery in fractured reservoirs. It is therefore important in the evaluation of waterflooding prospects to develop a thorough understanding of the quantitative relationships of the factors which control the rapidity of capillary imbibition. The imbibition process serves reservoir engineers in still another important way by providing a technique for studying the wettability of reservoir core samples. Such experiments are usually conducted by observing the rate of expulsion of oil or water from core samples submerged in the appropriate fluid. Several papers have been 'published on the experimental techniques involved?? Although Handy has recently published a method for calculating capillary pressures from experiments with gas-saturated cores, it has not yet been possible to deduce quantitative information regarding water-oil relative permeability and capillary pressure characteristics of the rock from the experimental results. Thus a technique is needed for studying the quantitative dependence of imbibition rate on oil and water viscosity, initial water saturation, relative permeability-saturation, and capillary pressure-saturation relations. The development of such information, including saturation and pressure profiles by laboratory experiments, would be very difficult. However, recent advances in the applica-

Jan 1, 1965

By H. E. Cline

A criterion for stability of a lamellar structure in two-phase alloys ulas derived froflr perturbation analysis of the steady-state solutions for eutectic growth. Both the mass and thermal diffusion were considered with boundary conditions for coupled lamellar growth. The concept of an effective liquidus slope was introduced. Excellent agreement was found with the experiments of Mollard and Flemings and the theory predicted a range of compositions on one side of the eulectic composition in which stable growth may occur with zero temperature gradient. The lamellar structure was found to be more stable than predicted by simple constitutional supercooling. The effect of impurity additions was also considered. DIRECTIONALLY solidified eutectic alloys may have a structure consisting of parallel plates or rods. It is important to know the range of composition and growth conditions that produce this ideal composite structure. Breakdown of the planar solid-liquid interface during growth may lead to a dendritic or cellular structure as in single-phase alloys. In the case of single-phase alloys, the concept of constitutional supercooling1 may explain breakdown, but in two-phase alloys the concept is no longer valid because the melting point is not given by the liquidus of one phase. Both phases must be considered in the mass and thermal diffusion problem. In the following theory a boundary condition at the solid-liquid interface was determined by using the steady-state equations of Jackson and ~unt' to relate the average temperature and composition of the liquid at the interface. Then a sinusoidal perturbation in the shape of the interface was introduced to test for stability as in the analysis of Mullins and ~ekerka.~ The wavelength of the important perturbation was assumed to be much larger than the lamellar spacing, h, and the perturbation was assumed not to change the lamellar spacing, A. Both the mass and thermal diffusion equation were used with averaged quantities at the boundary to determine whether the perturbation grows or shrinks in time. A stability criterion was determined by considering the case in which the perturbation remains stationary. An example of Pb-Sn alloys was used to illustrate the validity of this theory. I) THEORY A) Temperature at the Interface. The temperature of the interface in coupled two-phase growth is close to the eutectic temperature TE even for alloys off eutectic composition. It is assumed that the interface is in local equilibrium and one expects the temperature at the intersection of the liquid phase and the solid phases a and 0 to be near the eutectic temperature. The interface is undercooled because of curvature and the deviation of the composition in the liquid from the eutectic composition. The average interface temperature Ti for steady-state growth of a lamellar eutectic with a plane interface was determined by Jackson and Hunt' as given in their Eq. [17c]. Instead of using the ratio of the lamellar widths hp/h, as a parameter, the volume fraction of the 0 phase, +, will be used as it simplifies the algebra. The average interface temperature using this notation becomes: The average temperature of the interface depends on the volume fraction of the 0 phase, +, which may be related to the composition of the solid or the bulk composition of the liquid Co by: where C, and Cp are the average compositions of the solid phases. In steady-state growth, the composition of the solid is equal to the composition in the liquid far from the interface. The variation of the average interface temperature with composition is illustrated in Fig. 1 for the case of Pb-Sn alloys with a 10"4 cm lamellar spacing as calculated from Eqs. [I] and [2]. The interface temperature goes through a maximum at a composition near but not exactly at the eutectic composition. There is no reason why the composition corresponding to the maximum temperature should be equal to the eutectic

Jan 1, 1969

By B. B. Wood, M. J. Rzasa, E. W. Hough

An apparatus for the determination of fluid-fluid interfacial tension by the pendant drop method has been constructed. The apparatus is refined beyond those previously described in that the samples are introduced into an evacuated windowed cell without contacting air or mercury. Interfacial tension data are reported for the water-methane system between 74 and 280°F and 15 and 15,000 psi. The variation of wettability of the stainless-steel dropper tip with pressure and temperature was observed incidentally to the measurements. INTRODUCTION The most refined previous apparatus --for the determination of fluid-fluid interfacial tension at reservoir conditions was built by Hauser and Michaels.' These authors used the pendant drop method, previously investigated by Andreas, Hauser, and Tucker.' Recently, Fordham3 has made the pendant drop method an absolute one by numerical integration of the differential equation of drop shape. It is necessary to know the difference in density of the two fluids and the acceleration due to gravity in order to evaluate the interfacial tension from measurements of drop dimensions. Contact angle is not involved in the measurements by this method. Some precise tests of the method have been made at atmospheric pressure and 77°F by Douglas," who found agreement with established values to within about 0.1 per cent for water and air, cyclohexane and air, benzene and air, and benzene and water. No values of interfacial tension from the recent pendant drop apparatusJ have been published by Hauser and Michaels. However, values for the water-benzene and water-decane systems are given in Michaels' thesis5 or Dressures from 15 .to 10,000 psi and temperatures from 73 to 268ºF. Some determinations on gas-oil and gas-water systems in the range of 15 to 4,000 psi and 78 to 178°F have been made by Hozott,' by the drop-volume method. The capillary rise method has been used to investigate the interfacial tension in gas-oil systems at 88°F to several thousand psi by Swartz.7 Other determinations in gas-oil systems at pressures under 1,000 psia have been made by Jones: and Beecher and Parkhurst.9 The scarcity of data for pressures above 5.000 psi is apparent. APPARATUS An assembly drawing of the working section of the interfacial tenqion apparatus is given in Fig. 2. 'The endanant drop (A) is housed in a pressure chamber having glass windows (B). An agitator (6) is provided. The metal parts of the apparatus coming into contact with the sample are made of AISI: Type 316 steel. Drain grooves are provided around the whdows so that the sample in contact with the O-rings is not returned to the working section when the pressure is reduced. A schematic arrangement of the apparatus is given in Fig. 3. The working section containing the pendant drop (A) is filled by samples from bombs (B) with the aid of the vacuum system (C). The O-ring sealed pistons (D) used to separate the sample from the hydraulic fluid are shown together with the pressure production and measurement system (E). A front view of the apparatus is shown in Fig. 1. Most of the hydraulic system is at the top of the panel. The pump and displacement piston are at the lower right. The Variacs for the air bath heaters. the camera. and the vacuum gauge are a1so visible on the panel. The working section together with a window closure and a drop-forming tip closure are shown in Fig. 9. The carbon arc light source with filter is visible in the upper part of the picture. Two sides of the air bath with the impeller and a part of the convection shield can also be seen. The assembled air bath is shown in Fig. 10.

Jan 1, 1951