Search Documents

By A. L. Lee, M. H. Gonzalez, R. F. Bukacek

Experimental viscosity data for methane are presented for temperatures from 100 to 340F and pressures from 200 to 8,000 psia. A summary is given of the available data for methane, and a comparison is presented for that data common to the experimental range reported in this paper. Correlation of the data is presented, and predicted values are given for temperatures up to 460F and pressures up to 10,000 psia. INTRODUCTION The increasing ranges of temperature and pressure at which fluids are produced, transferred and processed in the petroleum and chemical industries stress the need for accurate information on physical properties, both for engineering calculation and for improving the methods used to estimate physical properties. A survey of the literature reveals that disagreements between published data on the viscosity of methane are common and that most investigations have been conducted over restricted temperature and pressure ranges. This situation made it desirable to undertake an investigation with a scope common to most of the literature data available. APPARATUS AND MATERIALS The apparatus used in this investigation was the capillary tube viscometer'described by Dolan,10 with modifications introduced in the general design and operation of the instrument. It has an effective pressure range from 14.7 to 10,000 psia and a temperature range from room temperature to 400F. The design of the viscometer is based on the establishment of a manometric head between two vessels containing the test fluid and a volume of mercury. The reservoirs are connected by a capillary tube through which the test fluid flows and a tube through which mercury flows. A pressure gradient is established by elevating one of the vessels above the other, the resulting flow of mercury displacing the test fluid through the capillary. The schematic diagram of the system (Fig. 1) shows the arrangement of the equipment auxiliary to the viscometer. The density cells assembly (E, Fig. 1) was not used in this investigation since reliable data are available on methane. 23 The assembly consists of a bank of eight 316 stainless steel pycnometers. Experimental values of the density of a natural gas19 and i-butane14 samples have been obtained for pressures up to 8,000 psia and temperatures up to 340F. The methane used was obtained from the Southern California Gas Co. Mass spectrometric analysis showed a composition of 99.8 percent CH4, 0.1 percent C2H6 and 0.1 percent C3H8. EXPERIMENTAL VALUES AND ACCURACY OF DATA Determinations were carried out at several temperatures. With the temperature in the cabinet controlled at a given level, a series of runs were made at various pressures. To establish pressures lower than that obtained from the methane cylinder, it was necessary to vent the gas slowly to the atmosphere until the desired pressure was reached. The consistency of the test run time for a particular condition can be used as a criterion for the accuracy of the data. A further check is obtained by comparison of the viscosity values obtained using different driving heads. To establish the instrument's reproducibility, a large number of runs were made at three completely different conditions of temperature and pressure. Data obtained were analyzed statistically to determine the most reliable value for the variance of the experimental data. It was found that the variances of the data taken at different experimental conditions were not statistically equivalent, but they were small. The largest variance was found at 8,000 psia and 100F, where the standard deviation was 0.385 micropoise. Therefore, 99.9 percent of the data can be expected to fall within 1.15 micropoises of the population mean, except for the intervention of systematic errors in operation.

By P. J. Root, J. E. Warren

An idealized model has been developed for the purpose of studying the characteristic behavior of a permeable medium which contains regions which contribute significantly to the pore volume of the system but contribute negligibly to the flow capacity; e.g., a naturally fractured or vugular reservoir. Unsteady-state flow in this model reservoir has been investigated analytically. The pressure build-up performance has been examined in some detail; and, a technique for analyzing the build-up data to evaluate the desired parameters has been suggested. The use of this approach in the interpretation of field data has been discussed. As a result of this study, the following general conclusions can be drawn: 1. Two parameters are sufficient to characterize the deviation of the behavior of a medium with "double porositym from that of a homogeneously porous medium. 2. These parameters can be evaluated by the proper analysis of pressure build-up data obtained from adequately designed tests. 3. Since the build-up curve associated with this type of porous system is similar to that obtained from a stratified reservoir, an unambiguous interpretation is not possible without additional information. 4. Differencing methods which utilize pressure data from the final stages of a build-up test should be used with extreme caution. INTRODUCTION In order to plan a sound exploitation program or a successful secondary-recovery project, sufficient reliable information concerning the nature of the reservoir-fluid system must be available. Since it is evident that an adequate description of the reservoir rock is necessary if this condition is to be fulfilled, the present investigation was undertaken for the purpose of improving the fluid-flow characterization, based on normally available data, of a particular porous medium. DISCUSSION OF THE PROBLEM For many years it was widely assumed that, for the purpose of making engineering studies, two param- eters were sufficient to describe the single-phase flow properties of a producing formation, i.e., the absolute permeability and the effective porosity. It later became evident that the concept of directional permeability was of more than academic interest; consequently, the degree of permeability anisotropy and the orientation of the principal axes of permeability were accepted as basic parameters governing reservoir performance. 1,2 More recently, 3-6 it was recognized that at least one additional parameter was required to depict the behavior of a porous system containing regions which contributed significantly to the pore volume but contributed negligibly to the flow capacity. Microscopically, these regions could be "dead-end" or "storage" pores or, macroscopi-cally, they could be discrete volumes of low-permeability matrix rock combined with natural fissures in a reservoir. It is obvious that some provision for the inclusion of all the indicated parameters, as well as their spatial variations, must be made if a truly useful, conceptual model of a reservoir is to be developed. A dichotomy of the internal voids of reservoir rocks has been suggested.7,8 These two classes of porosity can be described as follows: a. Primary porosity is intergranular and controlled by deposition and lithification. It is highly interconnected and usually can be correlated with permeability since it is largely dependent on the geometry, size distribution and spatial distribution of the grains. The void systems of sands, sandstones and oolitic limestones are typical of this type. b. Secondary porosity is foramenular and is controlled by fracturing, jointing and/or solution in circulating water although it may be modified by infilling as a result of precipitation. It is not highly interconnected and usually cannot be correlated with permeability. Solution channels or vugular voids developed during weathering or burial in sedimentary basins are indigenous to carbonate rocks such as limestones or dolomites. Joints or fissures which occur in massive, extensive formations composed of shale, siltstone, schist, limestone or dolomite are generally vertical, and they are ascribed to tensional failure during mechanical deformation (the permeability associated with this type of void system is often anisotropic). Shrinkage cracks are the result

By M. T. Hepworth, W. B. Sample

HYDROGEN solubility and linear contraction measurements were made at constant temperature on an alloy of titanium of constant oxygen-to-titanium molal ratio but variable hydrogen content. A cylindrical specimen of titanium (0.25 by 1.6 in.) containing 0.75 at. pct 0 and less than 0.03 at. pct N was cast from a master button of iodide process titanium which had been alloyed with titanium dioxide. Four vacuum fusion analyses, which were made along the length of the cylinder after the hydrogen solubility studies were made, revealed a variation of oxygen content of less than 0.02 at. pct. The specimen was introduced into a quartz furnace tube, annealed and outgassedat 1050"~ for about 1 hr, and then held at 800" * 0.5" for about 30 hr, during which latter time experimental data were taken. The leak rate to the furnace chamber was not sufficient to contaminate the specimen during the run as determined by leak-rate studies and vacuum fusion analyses for oxygen and nitrogen before and after the runs. The furnace chamber was attached to a modified Sieverts apparatus used to introduce measured quantities of purified hydrogen to the specimen. A quartz pushrod transmitted changes in length of the specimen to a Chevenard dilatometer. The push-rod passed through a vacuum grease seal to the dilatometer which amplified motion of the rod by a factor of 319. During hydrogen addition, the pushrod did not contact the specimen. Measured quantities of hydrogen were added to the specimen and equilibrium hydrogen pressures recorded as a function of hydrogen content. The solubility data, which appear in Fig. 1, agree well with a separate investigation made with a different apparatus. The solubility limit of the hcp a phase appears as an abrupt discontinuity in slope at 3.4 at. pct H. The solubility plot has been extended as a dashed line at this discontinuity in order to emphasize the discontinuity. The hydrogen content at which bcc P appears is not clearly defined. An arrow indicates solubility data were taken upon hydrogen addition. Upon reaching the arbitrary maximum hydrogen content of the run, the quartz pushrod was pressed against the specimen in order to record the change in length with hydrogen removal. A pressure of about 100 g per sq cm was required to overcome the friction of the vacuum seal in order to obtain reproducible linear contraction data. The contraction ratio ~2/1", change in length divided by original length, was set equal to zero for the maximum hydrogen content. The relative contraction of the specimen was observed upon removal of hydrogen in the following manner. The system was opened to the pumps for a few seconds and an unknown quantity of hydrogen was removed. The system was then isolated from the pumps by closing a stopcock and allowed to return to equilibrium. When the temperature returned to 800°C, the pressure stabilized, and contraction ceased, measurements of contraction and pressure were made. The hydrogen content of the specimen was then deduced by measuring the equilibrium PIydrogen pressure and referring to the solubility plot. Linear contraction was then plotted against interpolated values of atomic percent hydrogen. The justification for this procedure is the established reversibility of the hydrogen equilibrium.' Uniaxial compression was assumed to have no effect on hydrogen pressure. The isothermal linear contraction data are plotted on Fig. 1. An arrow indicates data were taken upon removal of hydrogen. The solubility limits of the two phase a i /3 region can be identified by discontinuities ill slope of the dilatometric plot at 9.1 and 3.4 at. pct H, agreeing well with other studies.' Dashed extensions of the plots emphasize the discontinuities. Isothermal studies of expansion made upon hydrogen addition gave curves of much smaller slope and less pronounced discontinuities. A uniaxial compressive stress tends to accentuate volume changes in going from the more open bcc phase to the hcp phase.2 Incomplete data on ex-

Jan 1, 1962

By B. A. Wilcox, A. H. Clauer

The tensile and creep deformation characteristics of fiber reinforced composites have been studied, primarily at room temperature. The matrix was an alloy of Mg-14 wt pct Li-1 wt pct Al (LA141A alloy) and the reinforcing fiber consisted of continuous wires of a high strength precipitation hardening stainless steel (AFC-77 alloy). The ultimate tensile strength and strength-to-density ratio of composites were much improved over that of the mmix, both at room ternperatwe and up to 200°C. Creep studies at room temperature revealed logarithmic weep behavior for wires and composites and steady state creep for the nonreinforced matrix. It was concluded that creep of the fibers controlled composite creep. The instantaneous logarithmic weep rate of the composite, 6, was related to the average tensile stress on the fibers, where n = 3.4. Mg-Li alloys are the lightest structural metals commercially available, e.g., the alloy LA141A (Mg-14 wt pct Li-1 wt pct Al) is 27 pct less dense than beryllium Although Mg-Li alloys have moderate tensile strengths at normal strain rates and ambient temperature, their creep strength is relatively poor, even at room temperature. Reinforcement by fibers offers a means of increasing the tensile strength as well as overcoming the limitations of low creep strength. With a proper choice of reinforcing filament, composites could be designed with higher strength-to-density and modulus-to-density ratios than the matrix alloy. This communication reports a study on the tensile and creep characteristics of alloy LA141A reinforced with continuous wires of a high-strength precipitation hardening stainless steel (AFC-77 alloy). The matrix has a bcc structure, with possible precipitates of AlLi and MgLi2A1,3 and is ductile at room temperature Stainless steel wire was chosen as the -reinforcing filament since previous work at Battelle4 showed that a good bond between wire and matrix could be produced in this system. For all volume fractions of reinforcing wire, the theoretical composite strength-to-density ratio is greater than that of the matrix and the modulus-to-density ratio is somewhat lower than that of the matrix. EXPERIMENTAL PROCEDURES The 0.004-in. diam AFC-77 wire was obtained from two suppliers: National-Standard Co.* and Crucible 0.13 Si, 0.18 V, 0.011 P, 0.012 S, balance Fe. In the as-received condition, wires from both sources had room-temperature ultimate tensile strengths of 585,000 to 590,000 psi with approximately 2 pct total elongation and 40 pct reduction in area. Composite preparation was accomplished by the vacuum infiltration technique. This involved placing the continuous wires in a mild steel tube (1/8-in. ID), capping one end of the tube with a sheet of the matrix alloy, and evacuating the other end of the tube. When the capped end was placed in a crucible of the molten matrix alloy and held for 31/2 min at 70O°C, the liquid metal was drawn up the tube approximately 7 to 8 in. A cross-sectional view of a typical composite is shown in Fig. 1. This preparation technique promoted excellent bonding between the matrix and wires, and porosity was rarely observed. Test specimens for both creep and tension studies were prepared by machining off the steel tube. Tensile tests on the wires and composite specimens were made in an Instron using a cross-head speed of 0.01 in. per min. The composite specimens had no reduced section and were held in self-tightening grooved serrated grips. Most of the specimens broke in the effective 1-in. gage section. Even though several tests fractured in the grip section, their ultimate strengths were essentially the same as those of specimens which fractured in the gage section. Creep specimens of the composite were gripped in a similar manner, and

Jan 1, 1970

By I. H. Silberberg, R. L. Reed, O. K. Kimbler

interfacial films have frequently been observed at interfaces between certain crude oils and water. Several investigators have postulated that the presence of these films should influence the efficiency of oil recovery in water drive or waterflood operations. They may also influence the stability of emulsions which are sometimes a problem in petroleum production, and may be a factor in the fonation of paraffin deposits in oil well tubing and flow lines. This paper presents a technique with which a modified Langmuir film balance may be used to study the compressibility and collapse pressure of these natural interfacial films. Experimental data are presented for several nude oil-water systems. Data developed are used to infer the phase state of the film as a function of such variables as rate of reduction of interfacial area, ionic composition of the subtrate and pH of the subtrate. A film of known physical characteristics is shown to have a sinnificant effect on oil recovery from an unconsolidated sand pack. Possible applications of these results to petroleum production are discussed. INTRODUCTION The use of water to displace petroleum from reservoir rocks is of major importance both as a primary and a secondary recovery process. As water invades the rock, oil is completely displaced from some pores and left as a discontinuous phase in other pores. The manner in which water moves from pore to pore is strongly influenced by capillary forces. In view of the complexity of reservoir fluid systems, there can be little doubt that complicated interactions take place at both the liquid-solid and oil-water interfaces. One of the more interesting, and least understood, of the phenomena which take place at the oil-water interface is the formation of interfacial films. These films are believed to result from the adsorption of high molecular weight polar molecules at the interface.l.2 Presence of such molecules may cause a striking alteration in interfacial tension. When the oil-water interfacial area of certain crudes is rapidly reduced, a thin region (film) about the interface assumes the appearance of a solid membrane, and striations, wrinkles and gross distortions may occur. If such a film is solid, it should greatly alter the interfacial tension normally assumed to exist between the oil and water phases. If the membrane is continuous, a solid phase would separate the oil and water. Interfacial films between crude oil and water were observed in 1949 by Bartell and Niederhausers who commented upon the apparent rigidity of the films and their possible importance in the petroleum industry. Morrell and Egloff4 had earlier attributed the extreme stability of emulsions of sea water in fuel oil to very stable asphaltic films. Numerous investigators have observed rigid films in the course of crude oil-water interfacial tension determinations by the pendent drop method. Several investigators5,6,2 have separated inter-facially active materials from crudes and attempted to characterize them chemically. Reisberg and Doscher,2 using Ventura crude, showed the interfacial tension against water (as measured by the pendent drop method) to be affected by aging, contraction and expansion of the interface, and the pH of the water. These investigators attributed the adhesion of oil to a water-wetted surface and the distortions of flow paths in glass capillaries to the presence of rigid films. Dodd7 has studied the interfacial viscosity of adsorbed films and found them to be non-Newtonian in behavior. Craighead and Harvey8 reported a series of displacements in tubes packed with 60 mesh glass beads. They interpreted the results as indicating an effect of stearic acid films on waterflood recovery and imply that natural films may produce similar results.

Jan 1, 1967

By G. S. Cole, G. F. Bolling

Recent experitnents on crystal growth have indicated that thermal convection in the liquid influences the conditions of solidification and subsequent crystal characteristics. However, the important columnar-equiaxed transition in castings has not been examined as a function of this free convection. An experiment is presented here to demonstrate a change in the transition as convection is suppressed. Lead- and tin-based alloys were cast in stainless-steel molds either without special procedures or with thin grids placed horizontally across the tnold diameter in portions of the castings. While the grids allowed some diffusion of heat and mass transport to occur during solidification, massive fluid flow was hindered in the vegiotz of the grids. A larger percentage of columnar grains appeared in the grid region than in other parts of the special castings or than in all regions of normal castings. This structural difference can be explained in terms of thermal convection. THERE is little information about the possible effects of natural convection which may occur during the solidification of a casting. Ruddlel-3 has stated that the movement of liquid metal produced by convection might assist in the transfer of heat, but in the absence of confirmatory evidence he suggested that this idea be treated with caution. chalmers4 considers, but does not assess, the importance of fluid motion in casting, while Tiller5 suggests that, in liquid metals which have extremely high thermal diffusivities, the effects of thermal convection may be negligible. In addition, throughout the literature, there are statements indicating the possible presence of thermal convection but without proofs or evaluations of its influence. Natural convection may be divided into two categories: thermal convection, where liquid flow is due to density differences caused by temperature gradients in the liquid (e.g., liquid near a cold wall falling with respect to hotter liquid further away from the wall); and solute convection, where flow is due to density differences as a result of concentration gradients, see Wagner.6 In most instances thermal convection predominates over solute convection. Nevertheless, even thermal convection has not been studied in any detail and is difficult to handle, since it can vary quantitatively from one experiment to another depending upon the exact positions of heat sources and sinks or the attitude of the container with respect to gravity. In spite of these difficulties it is possible to realize that a higher rate of heat transfer must result with convection than in any situation where conduction acts alone. Thus, the degree of increased heat transfer is the important undecided question in considerations of casting. It will be demonstrated here that a reduction in convective flow leads to a suppression of the columnar-equiaxed transition in a casting and that this suppression is related to a decrease in the degree of heat transfer. Fortunately it is possible to do this by comparative rather than quantitative measures. The insertion of horizontal plates or grids into a simple mold does not appreciably influence heat transfer per se, but in reducing the allowable height through which colder liquid may fall freely it reduces natural thermal convection. This leads to a reduced fluid velocity, a different flow pattern, and an altered heat transfer which can be measured. EXPERIMENTS The alloys were prepared from elements of 99.9 pct nominal purity and cast into an uncoated 2-in.-diameter stainless-steel mold insulated at the bottom with a suitable refractory so as to ensure predominantly radial heat flow during freezing. Casting was performed either in empty molds or in molds containing grids or plates in various configurations . The grids were constructed of aluminum screen (1/32 in. thick) which was cut to the desired shape, cleaned in aqua regia, and coated with a mixture of 2 parts aqua dag graphite -1 part Sauereisen cement mixture, thus supplying a hard, protective, and insulating layer approximately 0.006 in. thick. The

Jan 1, 1965

By M. M. Singh, H. L. Hartman

This paper reports on the results of a fundamental investigation of the behavior of rock under impact of a chisel-shaped bit. Drop-tester studies were conducted with strain gauges attached to the rock and striking bit-edge. It was established that rock behaves viscoelastically. Little correlation was found to exist between the strain waveforms detected on the rock surface and in the bit; thus, strain-time curves obtained with transducers mounted on the bit cannot be used to explain the sequence of events in rock breakage. It is estimated that nearly 1 per cent of the total energy of the blow is expended in permanent deformation of the rock. The strain- energy distribution on a rock surface as a result of chisel impact, which follows a simple law, is presented graphically. Experiments also were conducted with gauges embedded in the rock vertically below the line of impact and at an angle of approximately 45' to the vertical. No essential difference was observed in the type of wave obtained with embedded gauges and with gauges on the surface, as may have been anticipated from theory. INTRODUCTION Through the ages great strides have been taken towards improvement of the various techniques applied in the breakage of rock. Indeed, numerous applications of the rock-failure process are an integral part of present-day industry, including oil well drilling. Yet, astonishingly, little is known regarding the behavior of rock under dynamic loading or the mechanism of rock failure. The value of impact loading studies from the point of view of drilling is evident, however, and requires little elaboration because impact of a chisel edge on rock forms the basic action constituting most common methods of drilling. RELATED STUDIES Impact studies of communition and fragmentation conducted by Fahrenwald, et a1,8 Poncelet,l7 Axelson, et a1,2 Piret,16 Hartman and pfleider,ll Charles and deBruyn,5 and others have been of considerable help in understanding rock-breakage phenomena. Further insight into the basic action of drilling and penetration has been obtained by the work of Irving,13 Cherkasov and Shlyapin,6 BOrcherdt,3 Shepherd,19 Shreiner and Pavlova,20 Drilling Research Inc.7 Hartman,10 Appl and Gatley1 and several other investigators. However, knowledge concerning the stresses and strains induced in the rock and bit during impact remains scant. Recent photoelastic studies by Tandanand and Hartman (1961. 1962)22, 23 have shed light on stress distributions and Fracturing produced in plastics by chisel-shaped bits. In the field of transient strain Pulses and their role in rock drilling, the only extensive work has been performed by Drilling Research Inc. (DRI). They attached strain gauges to the shank of a bit and conducted drop-tester studies, using a die-shaped bit to obtain force-time diagrams. Their work led to the conclusion that, for a given rock type, the first peak of force occurs at a fixed level regardless of the total energy in the system, provided it exceeds a certain minimum threshold level. Further, on this basis an explanation was offered for the sequence of events in crater formation as the result of a single blow on rock. PURPOSE OF PRESENT STUDY The investigations reported in this paper were conducted in a manner similar to the DRI approach, but with the objective of determining the strain distribution in the rock around the area of bit impact and, to a limited extent, the strain in the chisel-shaped bit itself. It was expected that a knowledge of the disposition of these strains, along with the rate of attenuation of strain energy in the rock

By D. de Fontaine

A convolution procedure is described for obtaining an analytical expression for the composition profiles of solute concentration:; in binary solid solutions. The free energy of solutions containing such compositional inhomogeneities can then be evaluated and the stability of spherically symmetric clusters (G.P. zones, for short) us a function of their radii and central compositions can be studied accordingly. Two characteristic zone radii are defined: the minimum radius (RM) corresponding to the smallest zone for which growth is energetically favored over dissolution, and the critical radius RK) corresponding to zones associated with the minimum work of formation. The values of RM and RK coincide at the two classical limits of nucleation and growth and spinodal decomposition. As the spinodal is approached, however, RM and RK diverge increasingly and it is proposed that the initial] decomposition morphology is governed more closely by the minimum, rather than by the critical radius since the former possesses no singularity within the m eta stable phase boundary. Reversion is briefly discussed in the light of the foregoing and it is shown that a set of reversion temperatures located between the coherent miscibility gap and the coherent spinodal must correspond to a given distribution of low temperature zones. SINCE the early works of Guinier1 and preston,2 the kinetics and morphology of GP zones have been studied extensively, see Ref. 3, for example, for a recent review. Nevertheless, a completely satisfying theoretical treatment of the decomposition of meta-stable or unstable binary solid solutions is still lacking due to the almost unsurmountable mathematical difficulties that one encounters in attempting to solve the pertinent diffusion equation which, in principle, describes the whole process. This equation, first derived by Cahn for the initial stages of spinodal decomposition,4 is a nonlinear partial differential equation of the fourth order in anisotropic space, a general analytical solution of which is completely out of the question. From this equation, some prediction can be made, however, concerning the later stages of spinodal decomposition5 and a recent numerical solution6 of the equation in the case of one- and two-dimensional composition modulations of a solid solution inside the spinodal has led to the following considerations: the sinusoidal morphology7 is to be expected only for solutions of average composition c, located near the center of the coherent miscibility gap.8 For solutions with co closer to the spinodal decomposition, depleted shell zones9 should constitute the first recognizable decomposition product, later evolving to iso- lated, enriched preprecipitates in a uniformly depleted matrix as this must, in all cases, constitute the final metastable morphology, barring, of course, loss of coherency and change in crystal structure. In agreement with the qualitative description of Bonfiglioli and Guinier,l0 it thus appears that the various models discussed by previous authors are, in fact, idealized successive stages of the general decomposition process. Unfortunately, corresponding calculations for a solid solution outside the spinodal are not yet available because of two additional difficulties: it is essential that the diffusion equation be solved in three-dimensional space, and furthermore, one must include in the equation an appropriate "source term" describing the statistics of the creation of composition fluctuations required to trigger the nucleation events. At present, a general treatment of this problem does not appear to be tractable. Nevertheless, by using a modelistic approach, i.e., by selecting a priori certain composition variation profiles and by evaluating their associated free energies, one can draw certain conclusions relative to the stability of these isolated composition modulations or "zones", both outside and inside the spinodal. This is the object of the present investigation. THEORY We shall now describe a standard four-parameter spherically symmetric GP zone, evaluate its free energy according to the Cahn and Hilliard11 formula (Eq. [5], below), and plot the corresponding energy surfaces as a function of zone amplitude and radius as

Jan 1, 1970

By H. Biloni, B. Chalmers

Metallographic studies of the segregation of solute in and near the surfaces of chill-cast alloys have shown that, under sufficiently severe conditions, nucleation on the cold mold surface is followed by the formation of a disc of solid of the same composition as the liquid. This disc is surrounded by two concentric rings in which the concentration first decreases and then increases. The discs and rings are surrounded by a region in which arms of low solute concentration separate regions of high concentration. These arms start at the periphery of the outer ring and are initially radial; they make a gradual curved transition from the radial direction to a crystallographic direction. A normal dendritic substrcture is established beyond the position at which the arms have become crystallographically aligned. When the cooling is less severe, the disc becomes less irregular in shape and has a concentration minimum at the center. The term "Predendritic" is proposed for these transitions. A qualitative explanation is outlined for the phenomena observed at the surface and in the equiaxed region of the ingots. It is well-established1-' that the normal mode of solidification for alloys under most conditions is dendritic. A characteristic of dendritic solidification is the growth of rods or plates in specific crys-tallographic directions. The remaining liquid solidifies later, usually with a higher concentration of solute than the original liquid. Very little is known, however, about the transition from a crystal nucleus until the dendritic pattern is established. It is clear that there must be a transition because the nucleus is much smaller than the characteristic spacings found in the "steady-state" structure. And it also seems clear that the initial conditions must be important since the undercooling existing during the nucleation of the grains has a remarkable influence upon the grain size4 and the segregation pattern. The present paper reports a study of the micro-segregation in regions where crystals have nucleated, and the conclusions that can be drawn regarding the early stages of growth, in connection with the details of the solidification process. Aluminum of 99.993 purity, A1-Cu alloys containing 0.5 to 20 wt pct Cu, and Pb-Sn alloys with 5 wt pct Sn were cast and examined metallographically at the as-cast surface, close to the surface, and in the interior after casting. The alloys were cast from controlled temperatures against surfaces at different temperatures ranging from room temperature to approximately the liquidus temperature of the alloy. The surfaces were, in various experiments, graphite, stainless steel, copper, and aluminum. The methods used for metallographic examination were as follows. 1) The Nomarski phase contrast interferometer6 was used to obtain information on the topography of the as-cast surface. 2) The distribution of any solute present in A1 99.993 pct and copper in the A1-Cu alloys was studied by the anodic film technique originated by Lacombe and Mouflard7 and previously applied in the study of the solidification of aluminum by Biloni.8 The electropolished surface of the specimen is anodized under controlled conditions so that the thickness of the oxide film and therefore the interference colors under white light are sensitive to the local composition. This technique was applied to the as-cast surfaces and to metallographically prepared surfaces from the interior. 3) If in Al-Cu alloys, the surface is anodically oxidized for several minutes, instead of seconds as in the color-film technique, the copper segregation is outlined by a different level with respect to the matrix. This method is less sensitive but complementary to that for the color films. 4) The distribution of copper revealed and mapped by the metallographic techniques was confirmed by electron probe. 5) In samples with high concentration of copper, evidence for the copper distribution could be seen optically after electropolishing.

Jan 1, 1965

By D. H. Shaffer, W. F. Harris, W. M. Hickam, M. H. Loeffler

Experiments have been performed on the quantitative addition to niobium of enriched isotopes 018 and 017 at the 0.1 - 40 micro-gram level and their subsequent recovery. A measured quantity of molecular oxygen of known isotopic concentration is added to niobium metal by oxidation. Recovery of the enriched isotopes as carbon monoxide is by the vacuum fusion technique. Quantitative determinations of the enriched isotopes are made in standard volumes at constant temperature using accurately calibrated micromanometers and the mass spectrometer as pressure measuring instruments. THE presence of low concentrations of nonmetallic impurities in metals can often adversely affect the physical properties of the metal. The increase in hardness and the embrittlement of refractory metals due to oxygen is well known. In metallurgical studies dealing with the diffusion rates or the effect of oxygen on the physical properties of refractory metals, it is often necessary to add small quantities of oxygen to the metal and later recover the total quantity present. It is desirable, therefore, to establish the accuracy with which microgram quantities of oxygen can be added to niobium and recovered by vacuum fusion. The successful addition of microgram quantities of oxygen to metals by the usual manometric techniques and the subsequent recovery by vacuum fusion is difficult because of the uncertainty in the accurate determination of the quantity of oxygen originally present in the metal and because of the chance of contamination by atmospheric oxygen. In order to minimize these uncertainties earlier workers have added quantities of O16 at levels much higher than the oxygen concentrations that are of interest.'-3 In order to eliminate the uncertainty associated with the O16 additions to metals at low microgram levels, experiments have been performed on the quantitative addition of oxygen highly enriched in O17 and O18 at the 0.1 to 40 pg level to niobium and the subsequent quantitative recovery of the isotopes by vacuum fusion. Since the natural abundance of O17 is 0.04 pct and of 0" is 0.2 pct this enriched oxygen provides a uniqueness from interference not achieved in earlier work.'-= A measured quantity of molecular oxygen of known isotopic concentration is diffused into the metal. Recovery of the isotopes in the form of carbon monoxide is by vacuum fusion. The total quantity of gas is determined in standard volumes at constant tempera- ture using accurately calibrated Consolidated Electrodynamics micromanometers and the mass spectrometer for pressure measurements. The molecular and isotopic composition is determined mass spectrometrically. APPARATUS A schematic of the oxygen addition apparatus is shown in Fig. 1. Take-offs are provided for the attachment of the isotope container and for sample bulbs which are used to transfer samples to the mass spectrometer. Metal valves are used throughout the system including the sample bulbs in order to achieve a vacuum of 10"6 mm Hg. Evacuation is accomplished with a 3-stage oil-diffusion pump backed by a mechanical pump. A liquid nitrogen trap is used to remove condensables and to prevent oil vapor from entering the system. A micromano-meter is used for pressure measurements. The furnace assembly is shown in Fig. 2. It consists of a quartz tube joined to the apparatus through a 40/W joint with Apiezon W low-pressure wax. The furnace cap, also a 40/W pyrex joint, is provided with a pyrex hook. The sample is suspended from this hook with a platinum wire. The total volume of the furnace section is 435 ml. All volumes were determined by expanding a known volume of gas at a known pressure into the evacuated system and re-measuring the pressure. Both the micromanometer used on the oxygen addition system and the micromanometer used on the mass spectrometer were compared to a precision

Jan 1, 1961

By R. L. Pober, L. Kaufman, E. V. Clougherty

Composite structure of hafnium, zirconium, and titanium diboride with additions of metals and/or compound phases were prepared by reactive high-pressure hot pvessing and evaluated in air and in mixtures of oxygen and helium at temperatures up to 2400°K and at linear flow rates of 0.1 to 0.2 ft per sec stp. Oxidation characteristics, principally depth of conversion of diboride to dioxide, were determined by a gas analysis method and by postoxidation metallo-graphic analyses. The best oxidation resistance was observed for hafnium diboride-silicon carbide composite which exhibited a 6-mil conversion to dioxide for a 1-hv oxidation in air at 2200°K compared to a 20-mil recession for hafnium diboride with no additive. Analogous improvements in oxidation resistance were observed for the addition of silicon carbide to zirconium. EARLIER investigationof the oxidation behavior of the diborides of the transition elements of Groups IVA and VA in flowing gas mixtures have produced a ranking of these compounds according to which HfB2 is the most oxidation-resistant followed in turn by ZrB2 > TiB2 > TaBz > NbB2. Other ~tudies~-~ have examined the physical, thermal, and thermodynamic properties of these diborides to provide a sound basis for the more detailed and specialized investigations which are required to generate the information needed to assess the usefulness of such materials in high-temperature oxidizing environments. Recent oxidation studies7 have examined the effect of variations of metal and boron content on the oxidation characteristics of zirconium and hafnium diboride in flowing gas mixtures including air and helium plus oxygen. Results obtained in the latter investigation confirmed the superior oxidation resistance of HfB2 over ZrB2. Moreover, metal-rich compositions of hafnium diboride showed improved oxidation resistance over boron-rich compositions at temperatures up to 2000°K; at higher temperatures the differences are not distinguishable. In particular HfBI.7 exhibited a diboride to dioxide conversion depth of 4 mils in 1 hr at 2000°K; HfB2. 2 exhibited a conversion depth of 16 mils for the same exposure.~ Specimens for the study of metal and boron stoichiometry on the oxidation characteristics of HfB2 and ZrB2 were prepared by high-pressure hot pressing, a specialized type of fabrication which produced dense crack-free bodies of the desired compositions with no significant impurity introductions. The purpose of the present investigation was to explore the possibility of further improving the oxidation resistance of the metal-rich ZrB, and HfB2. The plan to accomplish this objective was based in the formulation of several diboride compositions containing other compounds and phases which are generally known to have good oxidation resistance in one or more types of test environment. Ideally, such additives would also enhance or at least not diminish other boride property values principally mechanical strength and those parameters related to thermal stress resistance. The additions selected for HfBz included zirconium, aluminum, and (Ta + TaB2) designed to produce ternary metal diborides, hafnium to form a discrete HfB phase, chromium to provide a metal skeletal phase, (Hf + Si) to form a discrete hafnium silicide phase, and Sic to provide a discrete second phase. Additions to ZrB2 included A1 + B designed to produce a ternary metal diboride, Zr + Si to form a discrete zirconium silicide phase, and Sic to provide a discrete second phase. All the specimens were prepared by high-pressure hot pressing which was used to advantage to fabricate billets of the variety of compositions desired and suitable for oxidation testing. The present investigation also included the evaluation of high-pressure hot-pressed TiB2and TiB2-20 vol pct Sic, conventionally hot-pressed Boride 2,' and KT silicon carbide. Boride

Jan 1, 1969

By C. A. Fothergill

An improved blotter model is described which makes use of a plastic sheet as base for the blotter field and pins to represent the input wells. The model is simple to construct and adapt and is particularly useful for demonstration purposes. Operational details are given. INTRODUCTION The blotter model was first applied to the investigation of oilfield problems by Wyckoff, Botset and Mus-kat,' primarily in order to demonstrate the flood patterns which could be expected from various arrangements of input and output wells in waterflooding projects. In recent years the model has been largely superseded for analog studies by the potentiometric and electric analog models which are more refined and versatile, since they can simulate the form, and, in certain cases, permeability variations which it is not possible to do with the blotter model. However, the latter is still used for analog studies where the reservoir can be represented as being uniform in thickness and permeability, and where a visual record of the flood advance is required. Ramey andNabor,2 for example, have applied the model to the investigation of areal sweeps in oil recovery by in-situ combustion. Because of its limitations, the blotter model can be used in general for only very simple arrangements of input and output wells, but within these limitations it has a definite advantage over all other types, in that it provides a picture of the advancing flood using the simplest and most easily constructed apparatus, by contrast with the gelatin model, for example, which requires more elaborate equipment. Moreover, the experiments can be quickly repeated and modified if necessary, whereas considerable time is required to prepare a gel field, compared with that required for the preparation of sheets of conducting blotting paper. For these reasons the blotter model can be very useful for demonstration and teaching purposes, enabling the student to experiment for himself with different arrangements of input and output wells. A model has been constructed for this purpose which is simple to operate and adapt. The model (Fig. 1) consists of a base platform or box, approximately 12-in. square, on top of which are screwed three lathes of 1/2 -in. square section to form three sides of a square 6 in. in length. This forms the support for a sheet of Cobex, transparent plastic, thickness .003 in., on which is carried the blotting paper representing the field. The plastic is flexible and relatively soft, and can be attached to the support by cellulose tape so that it can be easily removed when required. The input wells are represented by brass drawing pins pressed through the plastic sheet from the underside until the points are just projecting through the top; these are connected in the circuit as shown in Fig. 2. The sheet of blotting paper of required shape and size is saturated with a solution of 0.2 n-potassium sulfate to which has been added a saturated alcohol solution of phenol-phthalein at a strength of 10 cc/liter of stock solution. The surplus water is removed and the blotting paper placed in position on the plastic sheet making contact with the drawing pins and the sheet (which is an insulator). The anode consists of a second sheet of blotting paper cut as shown and saturated with the same solution. This is suspended from a laboratory stand, being connected to the conducting wire by means of a crocodile clip and making contact with the horizontal blotting paper at the projecting point. This is the method used by Ramey and Nabor and which permits the hydroxyl ions or colored front to reach the anode. As shown in Fig. 2, the electrical voltage used is 220-v DC and a shunt is placed across the terminals to step down the voltage as required. Further adjustment can be made as necessary using a lamp board in series with the shunt to obtain the

Jan 1, 1958

By F. D. Rosi

The slip and twinning behavior in extended titanium crystals were studied in some detail. The formation and appearance of coarse kink bands are discussed. Their crystallographic geometry was determined by X-ray analysis. A phenomenological interpretation of the complexities in kink band development is also presented. The critical resolved shear stress, coefficient of shear hardening, and plane of fracture were determined for several crystals extended at room temperature. THE slip and twinning elements observed in the room-temperature deformation of titanium were enumerated in a previous paper1 in which considerations were advanced regarding the nature and selection of these elements and their effect on the known mechanical properties of this metal. The present study concerns the crystallographic, microscopic, and mechanical aspects of flow in relation to the slip and twinning elements, and includes a prediction of slip systems, nature of slip and twin markings, inhomogeneities of plastic flow, and stress-strain characteristics. Unless otherwise noted, arc-melted titanium sponge (99.77 pct) was used in these experiments. The method of production of crystals, their dimensions, and the surface preparation for micrographic examination have been reported.&apos; For obtaining stress-strain characteristics, only those crystals which traversed the entire width of the specimen and were at least 8 mm in length were used. Tensile deformation of the crystals was performed with conventional grips for sheet specimens and a constant-stress loading beam, designed after the method of Andrade and Chalmers.&apos; The specimens were loaded by allowing sand to flow from a reservoir into a bucket suspended from the longer end of a balanced 6:1 lever arm at a rate controlled to load the specimen approximately 2 kg per min. Strain measurements were made using the Baldwin SR-4, bonded, resistance-wire strain gage, Type A-8, which permitted a reading accuracy of 2 microinches per inch. The formulas used in the evaluation of shear stress and shear strain, as in deriving the coefficient of shear hardening, are given by Schmid and Boasv n terms of the original orientation and change in length of the crystal. For calculating the critical resolved shear stress, the standard equation was used. The crystallographic nature of the unpredictable slip observed in a number of specimens was determined by the single-surface X-ray method of analysis as described in ref. 1. Experimental Results Prediction of Slip System: It was reported&apos; that room temperature slip in titanium takes place predominantly on {10i0} and in a <1120> direction, giving three potential slip systems. Hence, it should be possible to predict the operative primary system in the manner used by Taylor and Elam&apos; for alu- minum (i.e., the one with the greatest component of shear stress in the direction of slip). The stereo-graphic construction in Fig. 1 shows that this is true. In all cases, slip was found to occur initially on the (0110) plane and, in a number of cases, in the [21f0] direction. (The direction of slip was not determined for all orientations.) It follows that duplex slip can be expected when two slip systems are geometrically equally favorable for slip, as demonstrated in Fig. 2. In both crystals slip took place simultaneously on the (1700) and (olio) prismatic planes, as indicated by Fig. 2a and b. In the case of crystal B the movement of the specimen axis with increasing extension was toward the (10i0) direction, which represents the stable end orientation resulting from alternate slip on the [2ii0] and [1120] directions. This is analogous to the duplex slip process in face-centered cubic metals." Unpredictable Slip: In a number of crystals whose orientations were well within the operation of a single slip system, secondary slip occurred on planes not predicted by the criterion of maximum shear stress. Examples are shown in Fig. 3. In addition to

Jan 1, 1955

By L. W. Reynolds, R. F. Campbell, K. G. Carroll, S. H. Ballard

Based on metallographic and X-ray data probable equilibrium conditions from 1340" to 2400°F are presented for the composition range investigated. These are correlated with investigations of Takei and Kuo. Provision is made for existence of Mo, C and (Fe, Mo)2,C, carbides heretofore omitted in diagrams in this composition range. These carbides and (Fe, Mo),C, in equilibrium with ferrite and austenite, restrict the austenite and austenite + ferrite fields to lower carbon contents and higher temperatures than in the Fe-C system. OK groups of steel of different molybdenum content, but with approximately the same carbon variation in each group, Table I, were studied. The alloys were melted in an acid-lined 30-lb high-frequency induction furnace, and the ingots forged to l-in. rounds, which were then ground to 0.75-in. diam. To minimize chemical segregations a short length of each bar was given a high-temperature heat treatment as listed in Table I. Specimens 0.062 to 0.125-in. thick were thereafter cut from the homogenized bars for heat treatment to determine phase changes in the alloy. Heat treatments were performed in either a lead bath or muffle-type electric furnace, depending upon the temperature required. Maximum temperature variation was * 2" F in the lead bath and * 5"F up to 1900°F and +8"F above 1900°F in the electric furnace. Adequate protection of the specimen from the oxidizing atmosphere in the electric furnace was obtained by sealing the specimen in either a pyrex or a silica capsule which was evacuated or not as circumstances required. For austenitizing followed by quenching to an isothermal temperature, the specimen was placed in a closed-end refractory tube the mouth of which, extending beyond the furnace, was connected to a vacuum pump. The specimen could be withdrawn from this tube and quenched as desired. IDENTIFICATION OF PHASES The principal method of identification of phases was metallographic examination conducted on surfaces well removed from exposure to lead or atmosphere effects during heat treatment. Hardness measurements (Dph/lOkg) were also made on these surfaces. The etching reagents used and their effects are listed in Table II. Martensite and ferrite were readily recognized by their customary appearance after etching with saturated picral or 1 pct nital. The presence of liquid at high temperature was indicated by curved interfaces characteristic of boundaries between solid and melt, by intergranular pools of dendritic areas in a carbide matrix, and by a black network coincident with austenite grain boundaries. Examples of these structures are shown in Figs. l(a) and (0). The types of carbide which were attacked by the electrolytic chromic acid etch and the alkaline KMnO, etch, Table 11, were identified by X-ray analyses of extracted carbides of five alloys. Carbide residue was obtained from a metallographic specimen by etching the specimen for 24 hr in hydrochloric-picric acid in alcohol (grain-size reagent), after which the carbides, standing in relief, were scraped from the surface of the specimen and mounted on a glass slide. Diffraction patterns were obtained with a Philips spectrometer using iron-filtered cobalt radiation. Figs. 2, 3, 4, and 5 illustrate typical etching behaviors with chromic acid and alkaline KMnO,. Chromic acid attacked Mo,C as shown in Figs. 2(a) and (b)-but not (Fe, Mo),C, Fig. 3(b). Alkaline KMnO, attacked (Fe, Mo),C, Fig. 3(c), and also attacked M0,C. In accordance with these behaviors, a typical differentiation of coexisting Mo,C and (Fe, Mo),~ is shown in Figs. 4(a) and (b). The carbide (Fe, Mo),,C, could not be obtained free of other coexisting carbides for etching experiments. The behavior. however. of coexisting carbides identified by x-ray diffraction as Mo,C ind (Fe, Mo),,C, indicate that (Fe, Mo)23C6 is not etched by chromic

Jan 1, 1961

By William Hurst, A. F. van Everdingen

For several years the authors have felt the need for a source from which reservoir engineers could obtain fundamental theory and data on the flow of fluids through permeable media in the unsteady state. The data on the unsteady state flow are composed of solutions of the equation Two sets of solutions of this equation are developed, namely. for "the constant terminal pressure case" and "the constant terminal rate case." In the constant terminal pressure case the pressure at the terminal boundary is lowered by unity at zero time, kept constant thereafter, and the cumulative amount of fluid flowing across the boundary is computed, as a function of the time. In the constant terminal rate case a unit rate of production is made to flow across the terminal boundary (from time zero onward) and the ensuing pressure drop is computed as a function of the time. Considerable effort has been made to compile complete tables from which curves can be constructed for the constant terminal pressure and constant terminal rate cases, both for finite and infinite reservoirs. These curves can be employed to reproduce the effect of any pressure or rate history encountered in practice. Most of the information is obtained by the help of the Laplace transformations, which proved to be extremely helpful for analyzing the problems encountered in fluid flow. The application of this method simplifies the more tedious mathematical analyses employed in the past. With the help of Laplace transformations some original developments were obtained (and presented) which could not have been easily foreseen by the earlier methods. INTRODUCTION This paper represents a compilation of the work done over the past few years on the flow of fluid in porous media. It concerns itself primarily with the transient conditions prevailing in oil reservoirs during the time they are produced. The study is limited to conditions where the flow of fluid obeys the diffusivity equation. Multiple-phase fluid flow has not been considered. A previous publication by Hurst&apos; shows that when the pressure history of a reservoir is known, this information can be used to calculate the water influx, an essential term in the material balance equation. An example is offered in the literature by Old² in the study of the Jones Sand, Schuler Field, Arkansas. The present paper contains extensive tabulated data (from which work curves can be constructed), which data are derived by a more rigorous treatment of the subject matter than available in an earlier publication.&apos; The application of this information will enable those concerned with the analysis of the behavior of a reservoir to obtain quantitatively correct expressions for the amount of water that has flowed into the reservoirs, thereby satisfying all the terms that appear in the material balance equation. This work is likewise applicable to the flow of fluid to a well whenever the flow conditions are such that the diffusivity equation is obeyed. DIFFUSITY EQUATION The most commonly encountered flow system is radial flow toward the well bore or field. The volume of fluid which flows per unit of time through each unit area of sand is expressed by Darcy&apos;s equation as where K is the permeability, the viscosity and aP/ar the pressure gradient at the radial distance r. A material balance on a concentric element AB, expresses the net fluid traversing the surfaces A and B, which must equal the fluid lost from within the element. Thus, if the density of the fluid is expressed by p, then the weight of fluid per unit time and per unit sand thickness, flowing past Surface A, the surface nearest the well bore, is given as The weight of fluid flowing past Surface B, an infinitesimal distance 6r, removed from Surface A, is expressed as

Jan 1, 1949

By A. F. van Everdingen, William Hurst

For several years the authors have felt the need for a source from which reservoir engineers could obtain fundamental theory and data on the flow of fluids through permeable media in the unsteady state. The data on the unsteady state flow are composed of solutions of the equation Two sets of solutions of this equation are developed, namely. for "the constant terminal pressure case" and "the constant terminal rate case." In the constant terminal pressure case the pressure at the terminal boundary is lowered by unity at zero time, kept constant thereafter, and the cumulative amount of fluid flowing across the boundary is computed, as a function of the time. In the constant terminal rate case a unit rate of production is made to flow across the terminal boundary (from time zero onward) and the ensuing pressure drop is computed as a function of the time. Considerable effort has been made to compile complete tables from which curves can be constructed for the constant terminal pressure and constant terminal rate cases, both for finite and infinite reservoirs. These curves can be employed to reproduce the effect of any pressure or rate history encountered in practice. Most of the information is obtained by the help of the Laplace transformations, which proved to be extremely helpful for analyzing the problems encountered in fluid flow. The application of this method simplifies the more tedious mathematical analyses employed in the past. With the help of Laplace transformations some original developments were obtained (and presented) which could not have been easily foreseen by the earlier methods. INTRODUCTION This paper represents a compilation of the work done over the past few years on the flow of fluid in porous media. It concerns itself primarily with the transient conditions prevailing in oil reservoirs during the time they are produced. The study is limited to conditions where the flow of fluid obeys the diffusivity equation. Multiple-phase fluid flow has not been considered. A previous publication by Hurst&apos; shows that when the pressure history of a reservoir is known, this information can be used to calculate the water influx, an essential term in the material balance equation. An example is offered in the literature by Old² in the study of the Jones Sand, Schuler Field, Arkansas. The present paper contains extensive tabulated data (from which work curves can be constructed), which data are derived by a more rigorous treatment of the subject matter than available in an earlier publication.&apos; The application of this information will enable those concerned with the analysis of the behavior of a reservoir to obtain quantitatively correct expressions for the amount of water that has flowed into the reservoirs, thereby satisfying all the terms that appear in the material balance equation. This work is likewise applicable to the flow of fluid to a well whenever the flow conditions are such that the diffusivity equation is obeyed. DIFFUSITY EQUATION The most commonly encountered flow system is radial flow toward the well bore or field. The volume of fluid which flows per unit of time through each unit area of sand is expressed by Darcy&apos;s equation as where K is the permeability, the viscosity and aP/ar the pressure gradient at the radial distance r. A material balance on a concentric element AB, expresses the net fluid traversing the surfaces A and B, which must equal the fluid lost from within the element. Thus, if the density of the fluid is expressed by p, then the weight of fluid per unit time and per unit sand thickness, flowing past Surface A, the surface nearest the well bore, is given as The weight of fluid flowing past Surface B, an infinitesimal distance 6r, removed from Surface A, is expressed as

Jan 1, 1949

By W. M. Rutherford

A knowledge of the limits of miscibility between reservoir oil and possible injection fluids is required for selection of the optimum miscible-injection fluid. Limits of miscibility can be estimated from the results of equilibrium phase-behavior experiments. They can also be determined by means of displacement experiments conducted in a high-pressure sandpack. This paper describes the equipment and procedure which have been developed for determining miscibility conditions by stable displacement. A systematic series of displacements of a West Texas reservoir oil was carried out. The results indicate that, at constant pressure, miscibility is a function only of the pseudocritical temperature of the injection gas. This fact, together with improved experimental methods, makes the displacement technique a rapid, reliable means for determining miscibility conditions. In conjunction with the displacement experiments, phase diagrams were constructed for the oil with dry gas and propane and with dry gas and ethane. Phase behavior of the methane-ethane-propane system was determined at 110°F. The experimental work demonstrates the feasibility of using ethane-rich gases to reduce cost and pressure requirements for miscible displacement. INTRODUCTION In recent years, interest in the miscible displacement of oil by light hydrocarbon mixtures has been high. Many pilot and a few field scale projects have been started.1 These projects have made use of various methods for achieving miscibility: (1) the LPG-slug process, (2) the enriched-gas-drive process and (3) the high-pressure gas-drive process. Some field projects have been successful; the results of others are debatable. In general, projects which have performed best have involved the injection of an appreciable fraction of a pore volume of miscible material. Economical application of miscible displacement depends strongly on the cost of the miscible-injection fluid. If an appreciable fraction of a pore volume of material is required for successful application of these methods, a precise knowledge of the minimum requirements for miscibility in terms of composition of injection fluid is essential. Therefore, reliable experimental methods for determining miscibility conditions are important, and a procedure for estimating these conditions from the composition of the reservoir fluid is highly desirable. The subject of this paper is the problem of determining conditions which result in miscible displacement of oil by light hydrocarbon mixtures. Miscibility conditions can be estimated by means of equilibrium experiments conducted in a PVT cell, or they can be determined by means of high-pressure displacement experiments. This paper describes the equipment and procedure which have been developed for the determination of miscibility by high-pressure displacement experiments. These methods have been applied to the displacement of a West Texas reservoir oil with mixtures of dry gas, ethane and propane. In conjunction with the displacement experiments, triangular phase diagrams have been constructed for mixtures of the oil with dry gas and propane and with dry gas and ethane. The effect of injection-gas composition on conditions for miscible displacement in high-pressure sandpacks and cores has been the subject of several published papers.2-5 The experimental methods used in these investigations were such that displacements were unstable, and the effects of fingering and/or gravity layover are clearly evident in the results. Miscibility conditions were probably correct in spite of the instability phenomena, but the experiments evidently were time-consuming, and limited data were reported. Systematic high-pressure flow studies which would support a correlation of miscibility conditions have not been reported; however, Wilson5 has proposed the use of pseudocritical temperature of the injection gas as a parameter and Benham, et al,6 have based a miscibility correlation on observed and calculated equilibrium data.

By R. W. Heckel, R. A. Tanzilli

Treatment of the finite, diffusion-controlled, two-phase, moving-interface homogenization problem has been carried out by numerical methods for a wide range of input parameters. Solutions have been obtained for planar, cylindrical, and spherical diffusion geometries. The effects of parameters such as the interdiffusion coefficients (concentration-independent) of- both phases, the concentrations (time-independent) of each phase at the interlace, mean composition, initial size of the dissolving- phase, and initial composition of- the two phases are considered. The results are displayed to show both the position of the interjace and the degree of homogenization (expressed in terms of concentration at the matrix phase symmetry boundary) as a function of diffusion time. Typical exaw~ples to which the results may be applied are discussed. A limited number of exact, closed-form solutions to the multiphase, diffusion-controlled, moving-interface problem exists in the literature, e.g., Refs. 1 to 7. The exact, closed-form solutions generally consider the outermost phases to be infinite in extent and also assume an interfacial velocity function. It is through use of such boundary conditions that closed-form solutions are readily obtained. An extensive literature also exists on numerical solutions to the mathematically analogous solidification problem, e.g., Refs. 8 to 12. In addition, a numerical solution involving mass transfer with spherical interfaces and infinite boundaries has recently been carried out.13 A practically important class of diffusion problems exists, however, in which the boundaries of each phase away from the interface are finite. For example the finite boundary condition arises when one sets up a model for analyzing the homogenization behavior of two-component powder compacts in two-phase systems or for studying the solution behavior of precipitates. Although exact, closed-form solutions have as yet not been found to this class of problems, solutions are possible using numerical methods. The purpose of this paper is twofold: 1) to describe a numerical method for treating the two-phase, diffusion-controlled, moving-interface problem involving finite boundaries with planar, cylindrical, and spherical interfaces; 2) to present the results of a parametric study which identifies the role played by each variable appearing in the analytical model. R. A. TANZlLLI,Student MemberAIME,onleaveofabsence from the General Electric Co., Re-entry Systems Department and A. W. HECKEL, Member AIME, are Graduate Student in Materials Engineering, and Professor, Department of Metallurgical Engineering, respectively, Drexel Institute of Technology, Philadelphia, Pa. Manuscript submitted May 6. 1968. IMD MATHEMATICAL MODEL The three geometries analyzed are shown in Fig. 1. For each configuration a symmetry element is shown, across whose boundaries the mass transfer is zero. Mean composition, C, and the characteristic solute-rich phase half-thickness (or radius), 1/2, are the independent variables while the characteristic symmetry element half-thickness, L/2 (measured from solute-rich phase center to the matrix phase plane of zero mass transfer), is a dependent variable. The symmetry element for a planar interface accurately represents a real diffusion situation involving the stacking of symmetry elements of constant size. Because cylindrical and spherical symmetry elements of constant size cannot be stacked to fill space, the meaning of the dimension, L/2, for these situations may be considered in terms of the following approximations: and where NA is the number of parallel cylindrical elements per unit area on a surface normal to the cylinder axes and Nv is the number of spherical elements per unit volume. The model consists of a solute-rich phase matrix phase composite of limited mutual solubility. The concentrations at the interface are assumed constant with time; no assumption of local equilibrium at the interface is necessary except for the case where the

Jan 1, 1969

By Henry Leidheiser, Melvin C. Jr. Hobson

The rate of formation of the intermetallic compound, indium antimonide, at the interface between iudium and antimony at 100°C is greatly increased when a composite electrode of indium electrode -posited on antimony is made the cathode in an electrolytic cell. The amount of indium antimonide formed during electrolytic treatment at constant potential appeared to he proportional to the square root of time, but the scatter of the data, particularly at long times, makes it impossible to state this conclzlsi1:ely. The rate constant for the electrochemical synthesis was estimated to he 1.5 x 10-12 sq cm per sec. A similar- set of experiments was conducted in the absence of an applied potential but under otherwise similar conditions. The rate constant lor the thermal synthesis was less than 2 x 10-14 sq cm per, sec, at least two orders of magnitude less than obtained during electrolytic treatment. In a series of electrochemical studies it was observed that intermetallic compounds formed on the surface of polarized electrodes.1&apos;2 The possibility of using this technique for the electrochemical preparation of intermetallic compounds was investigated and subsequently applied to the synthesis of the semiconducting III-V compounds. Gallium antimonide, indium arsenide. indium antimonide, and indium bismuthide were successfully formed by electrodepositing the Group III element from a boiling acid solution onto an electrode made of the Group V element. The following kinetic studies were undertaken on indium antimonide to obtain quantitative information on the rate of formation of this representative III-V compound under a potential gradient. EXPERIMENTAL The electrolytic cell for these experiments was a modified, three-neck, round-bottom, 500-ml distillation flask. Sealed into the side of the flask was a Luggin probe which monitored the surface of the working electrode and extended out of the flask to a small reservoir containing a saturated calomel electrode. The counter electrode was platinum. The three electrodes were connected to a Wenking po-tentiostat which maintained a preselected constant potential near the surface of the working electrode with reference to a saturated calomel electrode. All potentials noted in this paper are with reference to a saturated calomel electrode at 45° ± 5°C. The potential and the current were both monitored using Sargent Model MR recorders. The antimony electrodes were cast in graphite molds to form spheres about 15 mm in diameter with a shank 6 mm in diameter and 10 mm long. A heavy copper wire (12B & S gage) was soldered to this shank. Collodion was applied for protection to the shank and the connecting part of the copper lead. A snug-fitting O-ring was slipped over the shank and the electrode was then fitted into a teflon holder such that only the spherical portion was exposed. The electrode surfaces were cleaned and chemically polished prior to assembly and introduction into the electrolyte by etching in a modified CP-4 solution (50 parts glacial acetic acid, 10 parts concentrated nitric acid, 2 parts 48 pct hydrofluoric acid). The electrolyte was 10-2 M In+3 in 1 MH2So4 and its temperature was held constant by refluxing. The reflux condenser was open to the atmosphere and no temperature correction was made for changes in barometric pressure. The indium concentration in the electrolyte was checked from time to time by the analytical method outlined below for the determination of the amount of indium antimonide formed electrochemically. A schematic of the recorded potential is shown in Fig. 1. At point A the antimony electrode was introduced into the cell and the open-circuit potential recorded. At point B the potentiostat, preset at the desired potential, was switched on. It was switched off at point C and the potential fell to the open-circuit value of the electrodeposited indium. The

Jan 1, 1965

By M. H. Waxman, L. J. M. Smits

A simple physical model was used to develop an equation that relates the electrical conductivity of a water-saturatedshaly sand to the water conductivity and the cation-exchange capacity per unit pore volume of the rock. This equation fits both the experimental data of Hill and Milburn and data obtained recently on selected shaly sands with a wide range of cation-exchange capacities. This model was extended to cases where both oil and water are present in the shaly sand. This results in an additional expression, relating the resistivity ratio to water saturation, water conductivity and cation-exchange capacity per unit pore volume. The effect of shale content on the resistivity index - water saturation function is demonstrated by several numerical examples. INTRODUCTION A principal aim of well logging is to provide quantitative information concerning porosity and oil saturation of the permeable formations penetrated by the borehole. For clean sands, the relationships between measured physical quantities and porosity or saturation are well known. However, the presence of clay minerals greatly complicates log interpretation, particularly the electrical resistivity and SP logs, and considerably affects evaluation of hydrocarbon-bearing formations. The conductance and electrochemical behavior of shaly sands and their relation to log interpretation have been studied by many workers. wylliel and Lunch2 reviewed this work in some detail. Virtually all laboratory measurements of electrical resistivity and electrochemical potential of shaly sands published to date are the work of Hill and Milburn.3 Their measurements were made on about 300 cores covering a large variety of sedimentary rocks, and a wide range of equilibrating NaCl solution concentrations. Hill and Milburn described their conductivity data by an empirical equation in which the shaly sand conductivity Co was expressed as a function of the solution conductivity C and two parameters b and F01 . The quantity b was shown to be a measure of the effective clay content of the rock, being approximately proportional to the cation-exchange capacity of the rock divided by its pore volume. The latter ratio is designated as 9, in this paper and has the dimensions meq/ml or equiv/liter. Qv is identical with the term representing the concentration of fixed charges in the Meyer-Sieved and Teorel15 theory of permselective membrane behavior. F01 is a formation resistivity factor referred to a hypothetical equilibrating solution resistivity of 0.01 ohm m at 25C,* where clay effects presumably are minimized. F01 was correlated to porosity by an Archie-type equation.6 The Hill-Milburn equation describes their data with a standard deviation of approximately 1 percent and a maximum deviation of + 10 percent. Shaly sands behave as permselective cation-exchange membranes, their electrochemical efficiencies increasing with increasing clay content. The electrochemical potential data of Hill and Milburn were expressed graphically, and demonstrate that the membrane efficiency (or cation transport number) of these sands is a function only of the b value (i.e., Qv) and the respective salt concentrations of the two solutions forming the liquid junction in the sand. The diffusion potentials are not dependent on F01 or any parameter relating to the porosity or pore geometry of the rock. The Hill-Milburn resistivity equation correctly predicts a decreasing sand conductivity Co with decreasing solution conductivity C?. However, at some low value of C,, the calculated Co-C, function goes through a minimum; with further decrease in C?, the predicted sand conductivity increases sharply. As pointed out by Hill and Milburn, an increasing value of Co with decreasing C, is physically meaningless. Since this occurs below the range of practical values of C,, the usefulness of the empirical equation is not affected when

Jan 1, 1969