Search Documents

By P. T. Bail, J. E. Bobek

In reservoirs containing soft shales interbedded with oil-producing sands, it is possible to have water movement from the shales into the sands during reservoir depletion. This paper presents the results of scaled flow-model tests performed to investigate the efficiency of oil displacement by this vertical influx of water into oil-producing sands. In these model tests, the effects of rate of influx and rock wettability were considered. Oil recoveries by shale-wafer influx and by water flooding were cornpared after 0.4 hydrocarbon volume of water entered the producing sand. This comparison was made for water-wet and oil-wet systems. In water-wet sand, shale-water influx produces slightly less oil than conventional water flooding; in oil-wet sand, shale-water influx produces about one-half as much oil as conventional water flooding. In another type of comparison, it was found that shale-water influx followed by water flooding recovers less oil than water flooding alone. This is true for both water-wet and oil-wet sands. A comparison of tests run at different rates of shale-water influx shows that oil recovery is insensitive to rate in water-wet sand; in oil-wet sand, oil recovery decreases with decreasing rate of influx. INTRODUCTION Localized subsidence of the earth's surface has accompanied the depletion of some oil reservoirs, particularly in California,1, 2 along the Texas Gulf Coast,3 4 and in the Bolivar Coastal fields of Venezuela. Subsidence occurs when fluid withdrawals reduce the pressure in the reservoir. The reasons for subsidence accompanying a reduction in reservoir pressure are not known with certainty.10° However, there is considerable evidence that compaction of soft shale beds adjacent to producing oil sands can be a major cause of this subsi-dence.3, 6, 10-12 When the shales compact, a limited amount of water moves from the shales into the producing sands. A schematic illustration of this oil-producing mechanism is shown in Fig. 1. As the pressure in the sand declines due to oil production, a predominantly vertical pressure gradient develops at the sand-shale boundary. This causes water to move from the shale into the sand as indicated by the solid arrows. This paper reports the results of scaled flow-model tests to investigate the efficiency of oil displacement by the vertical influx of shale water into producing sands. The model tests were not performed with the aim of duplicating the recovery behavior of any specific reservoir. The primary purpose of the tests was to obtain a better understanding of oil displacement by shale-water influx. A knowledge of the oil- and water-production characteristics of even a simplified reservoir system should be of value to engineers estimating the performance of reservoirs where shale-water influx is occurring in conjunction with other production mechanisms. The following limitations and assumptions apply to the model work. 1. The sand of the prototype reservoir is homogeneous and horizontal. 2. The shale-water influx is uniformly distributed along the top and bottom of the sand member being studied. 3. No effects of edge-water drive or dissolved-gas drive that might accompany shale-water influx are considered. PROTOTYPE RESERVOIRS Scaled-model tests were performed to study oil displacement from three prototype reservoir sections. The reservoir sections modeled consist of single sand layers out of a sequence of -sand and -shale layers, as illustrated in Fig. 1. Each prototype represents a homogeneous, high-permeability, unconsolidated sand layer. In the field, the prototype is oriented in such a manner that one end is at the radius of drainage of a producing well and the other end is pointed toward the producing well (Fig. 2). The sand and fluid properties and the dimensions of the three prototypes are given in Table 1. The capillary pressure relationships for the three prototypes are presented in Figs. 3, 4 and 5. The oil-water per-

By Fred H. Poettmann

An equation has been derived for use in calculating the sandface pressure of flowing gas wells in which the variation of the compressibility factor of the gas with pressure is taken into consideration. This variation due to compressibility has been put into both graphical and tabular form. Comparison of calculated results with field measured results were made on 20 dry gas wells from a given field and 11 distillate wells from different fields. The agreement between calculated and observed results is good. In addition to their use in the calculation of sandface pressures of flowing wells, the factors can be used in the direct calculation of the static bottom hole pressure of gas wells, the capacity of gas transmission lines, and in the calculation of the theoretical isothermal horsepower necessary to compress a natural gas. Examples demonstrating the use of various equations are given. INTRODUCTION Many of the equations used to calculate the sandface pressures of flowing gas wells from well head data do not take into consideration rigorously the deviation of the natural gas from ideal gas behavior. For low pressure wells this error is not serious. For high pressure wells flowing at high rates this error call be serious. Sandface pressures are used to determine the producing capacity of gas wells as shown in the U. S. Bureau of Mines' Monograph 7.' An equation has been derived for use in calculating the sandface pressure of flowing gas wells in which the variation of compressibility of the gas with pressure is taken into consideration. This variation due to compressibility has been put into both graphical and tabular form. In addition to their use in the calculation of sandface pressures, these "factors" can also be used to give a direct solution in the calculation of the static bottom hole pressure of gas wells, the calculation of the capacity of gas transmission lines, and the calculation of the isothermal horsepower necessary to compress a gas. CALCULATION OF SANDFACE PRESSURE OF FLOWING GAS WELLS The starting point in the derivation of any specific flow equation is an energy balance on the fluid flowing between any two points in the system. This energy balance based on one pound of flowing fluid is expressed by the well known general flow equation: jp~- vdP+ Ah + -— + wt + r. = o . . (i) P, 2 g. where V = specific volume of flowing fluid in cu ft/lb P = pressure in Ib per sq ft ah = difference in height above datum plane in, ft ------ = kinetic energy change of one lb of flowing fluid, v equals velocity in ft/sec and g. = 32.174 Ws = Work done by fluid while in flow (similar to shaft work in driving a turbine) W, = Znergy losses due to the irreversibilities of the flowing fluid This equation contains no limiting assumptions and can be made the basis of any fluid flow relationship by limitations or substitutions. The energy units are usually expressed in term* of foot-pounds of energy per pound of mass. Considering the vertical flow of gas through tubing or annu-lus, the external work W, done by the gas is zero. Let us also assume that the value of the kinetic energy function is small and can be neglected. Equation (1) is reduced to: /P: VdP + Ah+Wt = 0 .......(2) Pi The energy loszes PI can be expressed in terms of the well known Fanning equatinn: 4/Lf: rf = —— ............(3) 2g,D where L = the integrated average velocity in ft/sec of vertical flow in a uniform flow string -~h = L v = the integrated average velocity in ft/sec / = dimensionless correlating function

Jan 1, 1951

By V. Kondic, W. M. Rumball

HUNT and Jackson' have recently suggested that binary eutectics might be classified according to the entropy of fusion of the component phases. Such a classification was in fact proposed2 some time ago by the authors and the purpose of this communication is to present some supporting evidence. As pointed out by Hunt and Jackson' typical metals of low entropy of fusion, ?S, tend to form normal eutectics (using Scheil's terminology3). Where one phase is metallic (low ?S) and the other nonmetallic (high AS) then anomalous structures are formed. These anomalous structures may either be regular (e.g., Bi-Sn)or irregular (e.g., A1-Si). Finally Hunt and Jackson suggest that where both phases are of high ?S regular faceted structures are formed. In the proposed classification2 the ratio of the entropies of fusion is taken as the quantitative parameter, i.e., ?Sa/?SB = E, where e > 1.0. An extensive literature survey has been carried out2,4 for reliable photomicrographs of eutectic structures. The structures, numbering 151, have been classified into two classes, viz., normal (i.e., lamellar and rodlike) and anomalous. The latter group has been arbitrarily subdivided into regular and irregular anomalous eutectics. Details and sources of the structures are given elsewhere.2,4 In the absence of better thermodynamic data the value of E has been determined for each of the systems making the following assumptions and approximations: l) The entropy of fusion of solid solutions may be calculated additively;5 ii) Inter metallic compounds are treated similarly but with the addition of the following term for fully ordered structures.5 -4.573(Ni lOg10 Nl + N2 log10 N2) where Nl and N2 are the atomic fractions. It might be argued that the entropies of fusion are meaningless as the eutectic phases form at some temperature (T,°K) below their true melting point (Ta or Tb). As AS, = Lo/To if La the latent heat of fusion of a, is assumed independent of temperature, the actual entropy change AS& is given by La/Te. The ratio of entropies of fusion then becomes Fig. 1 shows E' plotted against cumulative totals for normal and anomalous systems from which it is seen that normal eutectics predominate at low values of E'; e.g., some 90 pct of normal systems lie below E' = 1.5 with 65 pct of anomalous systems above this value. Fig. 1 also shows that the regular anomalous structures occur mostly at low values of c' possibly suggesting growth more akin to normal than anomalous growth. The choice of c as a classification parameter was based on nucleation concepts, it being observed that metals widely separated in Sundquist and Mondolfo's8 nucleation series were anomalous and those close together normal. The basis of these classification views is really outside the intended scope of this communication but it is hoped to explain them fully in the near future.7 Having arrived at this entropy of fusion classification by nucleation considerations it is interesting to note that Hunt and Jackson's ideas are based on growth considerations. Elsewhere one of the authors4 has shown that an even better classification of the metallographic structures reviewed is obtained from purely conventional heterogeneous nucleation concepts. It has been found that systems having a value of the parameter +b less than 0.28 are almost certainly normal and above this anomalous: V, is the gram atomic volume. This parameter is nevertheless closely related to c and c', these being an approximation of p. The parameter is also related to the other classification3,8,9 which have met with any success in the past. All these classifications are to be reviewed shortly by the authorsU7 A significant thing about the classifications based on

Jan 1, 1968

By E. Miller, N. Chou, K. Komarek

The solid solubilities of the compounds PbS, PbSe, and PbTe have been determined by equilibrating single crystals with two phase alloys. The solubilities were determined by measuring the Hall coefficient of the equilibrated samples. The observed maximum deviations from stoichiometry on the lead rich and the chalcogen-rich side of the phase fields were: at. pct excess Pb at pct excess S at 850°C; PbSe--4 x 10'" at. pct excess Pb at at. pct excess at 800°C; and PbTe-2.8 x 10 4 at. pct excess Pb at 780°C, 10"1 at. pct excess Te at 720°C. The previously developed theory of retrograde solubility has been expanded to include intrinsic disorder in the compounds, and the experimental results are in agreement with the theory. THE solid solubility limits of intermetallic semiconducting compounds are of considerable interest, since many of these compounds have potential use as thermoelements in thermoelectric devices, as infrared detectors, or in semiconductor devices. The solidus curves of a significant number of these compounds have been reported to be retrograde with the maximum solid solubility occurring between the eu-tectic and the maximum melting temperatures. In a previous paper,' equations were derived which permit calculation of the maximum solid solubility of an AB compound exhibiting retrograde solubility. Derivation of these equations was based on the assumption that the compound is a nondegenerate semiconductor at the temperature under consideration, that the defects were unassociated, and that the vibrational contribution to the entropy of the compound could be evaluated by the Einstein approximation. It was shown that at the temperature of the maximum solid solubility The entropy term ASS is composed of the entropy of formation of the solid stoichiometric compound (AS$)*, and the entropy due to deviations from stoichiometric composition which is the sum of a con-figurational term (aS*), a vibrational term (AS,), and an electronic term (AS,). In the previous paper the positional entropy was evaluated considering the con- N. CHOU, Member AIME, formerly Graduate Student, New York University, N. Y., is with IBM, Yorktown Heights, N. Y. K. L. KOMAREK, Member AIME, formerly with New York University, is Professor, Department of Inorganic Chemistry, University of Vienna, Vienna, Austria. E. MI LLER, Member AIME, formerly with New York University, is Associate Professor of Mechanical Engineering, California State College at Long Beach, Calif. The paper is based on a thesis submitted by N. CHOU to the Graduate Division, School of Engineering and Science, New York University, in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Materials Science. tributions from defects arising from deviations from stoichiometric composition. The derivation is generalized in the appendix to include intrinsic disorder. The electronic entropy term is also recalculated in the appendix. Of the various semiconducting compounds, the solidus curves of Pbs,' and ~b~e,~ have been most extensively studied. The solidus of PbSe has also been determined, but the temperature of maximum solubility has not been accurately established.5129 The present investigation was undertaken to compare the experimental solid solubility data for these systems with the values predicted theoretically. Since the temperature and composition of the maximum solubility of PbSe has not been accurately determined, experimental emphasis was placed on determining the solidus curves for PbSe in the region of maximum solubility, although experiments were also performed on PbS and PbTe at several temperatures. X-ray diffraction studies were also made on carefully prepared powder samples of PbSe to determine the variation of lattice parameter with composition. EXPERIMENTAL PROCEDURE The procedure of Brebrick and Allgaier was employed for determining the solidus data. Monocrystal-

Jan 1, 1970

By John G. Hall

The pumping operations at the Chief mine have been unique in the respect that for many years the entire flow of water into the mine has been disposed of by pumping into natural underground "caverns" or "sinks." Volumes up to 4000 gpm have been pumped into these caverns with no apparent tendency for the caverns to fill with water, and it has never been proven that any of the water returns to the mine workings. Major pumping operations have been carried on during two periods, the first from 1918 to 1927, and the second from 1942 to the present time. Early Pumping Operations CHIEF MINE The Chief Consolidated Mining Co. was organized in 1909, and for the following nine years ore was developed and mined from the 600 level down to the 1800 level, all above the permanent water table, the 1820 level. The ore occurs as extremely-irregular pipes of replacement ore in limestone, and to date a total of approximately 2 1/2 million tons of ore has been produced. The ore was highly oxidized down to the permanent water level, and some highly-oxidized ore runs occur beneath the permanent water table. In 1916, with extensive mining being done on the 1800 level, 20 ft above the water table, there were indications that a zone of secondary enrichment might exist below the water, and in order to follow the ore runs to greater depth, it was realized that arrangements would have to be Made to pump water to the surface. The original pumping installation was designed to pump 350 gpm in five stages against a 1900-ft head to the surface. By January 1919, the water had been lowered to the 1900 level and an excellent grade of silver-lead ore was being produced from the zone of secondary enrichment. At this point. the flow of water equalled the capacity of the pumps, and pumps of an additional 1000 gpln were installed to supplement the original system. By January 1920, the water had been lowered approximately 100 ft, and the total flow of water was 1000 gpm. The water was being pumped from many different places in the mine where ore was being followed downward, with 27 different pumps operating below the 1800 level, mainly Cameron air-driven sinker pumps. These first sinking pumps were cumbersome and inefficient, and each of the larger pumps weighed approximately three tons and required tremeridous volumes of compressed air. During this period, an exploration drift being extended north and west from the Chief shaft on the 1600 level encountered a large natural cavern some 2000 ft from the shaft. This cavern in places measured 60 ft in diameter. Sometime later the disposing of water into this cavern was suggested, and as an experiment to determine the results of such an operation, an attempt was Made to fill the cavern by pumping in a small volume of water, but the water disappeared as fast as it was discharged into the cavern. It was decided then to try the cavern as a drainage outlet for the entire volume of water being pumped to the surface, and the water apparently did not return to the mine workings. Subsequently water up to 3000 gpm was pumped into this cavern. An attempt was made to determine whether any water was being returned from the cavern to the lower workings, and 100 lb of sodium fluorescein salt was added to the water discharged into the cavern, but no trace of it was found in the drainage water in the mine. Results of another such test will be described later.

Jan 1, 1950

By Klaus Detert, Joseph Ziebs

The progress of primary recrystallization has been studied in alloys of nickel with 0.01, 0.1, and 1.0 at. pct Ta. It can be shown that in one particular alloy containing 0.1 Ta and 0.02 02 a large grain structure (0.5 mm diameter) could he developed by primary recrystallization in thin sheetswhich were cold-worked to 97.5 pct reduction in thickness. Electron transmission microscopy demmstrates that this phenonmena was linked to the occurrence of fine precipitates at subboundaries. A discussion of how fine precipitates can exert under certain circumstances a rather distinct retarding effect on nucleation rate is included. IN a former investigation, the recrystallization and texture formation in a series of binary nickel alloys including pure nickel were studied.' Nickel alloys with some tantalum added displayed peculiarities during primary recrystallization which have been studied in more detail by the present investigation. The alloys were prepared by powder-metallurgy techniques similar to the previous investigations.' Pure metal powders were thoroughly mixed to obtain alloys with 0.01, 0.1, and 1.0 at. pct Ta, referred to as alloy A, B, or C, respectively. The metallic impurities were kept below 0.005 at. pct except iron and cobalt with slightly higher values, by using highly purified nickel and tantalum powders.' The mixtures were pressed into bars and sintered at a temperature of 1300°C. The sintering atmosphere was either high vacuum, dry hydrogen, or wet hydrogen. The sintered rectangular bars were rolled 60 pct to a thickness of 3 mm, then annealed under dry hydrogen at 1000°C and again cold-rolled to about 97.5 pct reduction of thickness. After this heavy deformation, the sheet thickness was 0.07 5 mm. In order to study primary recrystallization which occurred in measurable times between 400" and 600°C, samples of the rolled sheets were annealed in a salt bath at various temperatures for various times in order to study the progress of recrystallization by hardness measurements and microscopic examination. For optical observation, the samples were polished and etched in a solution of 65 cu cm HNO3, and 18 cu cm acetic acid in 20 to 100 cu cm H2O, depending of the particular sample and annealing treatment. For transmission electron microscopic observations, thin foils were prepared by the combined window and Bollmann techniques.273 The electrolyte consisted of 23 pct perchloric acid and 77 pct acetic acid held at 15°C. Polishing was performed at a current density of 0.4 amp per sq cm. The examination of the thin films was accomplished by a Siemens Elmiscope-type electron microscope operating with an accelerating voltage of 100 kv. RESULTS The progress of primary recrystallization was very accurately followed in these alloys in order to establish the time when recrystallization started and ended during an isothermal anneal. Table I gives some of the data obtained by these studies. The time for beginning is a mean value, where somewhat less than 5 pct of the matrix has recrys-tallized. The end of recrystallization is the time when at least 95 pct of the matrix has recrystal-lized. Listed are three different B alloys, with 0.1 Ta distinguished by different sintering atmospheres which influence the rate of primary recrystallization. By the above method, the state of the grain structure when primary recrystallization was completed, but grain coarsening had not started, could be clearly defined. The average grain diameter of the recrystallized matrix is also given in Table I. The grain size did not depend on the temperature of the isothermal anneal. With one exception the grain size of the primary recrystallized matrix was found to be rather fine

Jan 1, 1965

By J. L. Moriarty, Irving Johnson, H. M. Feder

The crystallization of ternary alloys has been studied as a coprecipitation phenomenon by measuring the substitutional carrying of trace metals with CeCd11 precipitating from liquid cadmium. The Doerner-Hoskins logarithmic distribution law has been found to apply in all cases; the coprecipitation coefficients obtained were: 1.48 (La), 1.06 (Th), 0.63 (Pr), 0.24 (Gd), 0.17 (Sm), 0.13 (U), 0.10 (Er), 0.09 (Sr). Li, Na, K, Ba, Y, Sc, and Zr were not measurably carried. It has been observed that coprecipitation was most extensive for those tracers which form with cadmium an intermetallic compmdnd isostructural with CeCd11 or closely related thereto. For these tracers an inverse relation between Solubility in cadmium and coprecipitation coefficient was noted. THE partition of impurities between the liquid and solid phases of a metal or a congruently melting intermetallic compound has been extensively studied by the methods of zone melting.' However, when the composition of a liquid alloy and of the solid phase crystallizing from it differ grossly, partition of an impurity cannot be investigated in this manner. Partition of the impurity under equilibrium conditions may be studied by investigation of the appropriate phase diagram and partition under nonequi-librium conditions may be studied by observation of coring in the solidified alloys.2 The nonequilibrium case is of special interest in metallurgical practice because many complex alloys are cooled too rapidly to permit the uniform diffusion of impurities throughout the crystallizing phase(s). With this in mind, the present study was undertaken to obtain quantitative data on the phenomenon of coprecipitation via substitutional solid solution in certain metallic systems. It will be noticed that rather than adopting metallurgical ter- minology we have chosen to refer to the partition of impurities during alloy crystallization as coprecipitation.3 In the terminology of coprecipitation the impurity is referred to as the tracer and the precipitating phase as the carrier. From aqueous solution, coprecipitation may occur as a result of adsorption on the surface of the carrier, by occlusion of the liquid within the crystals, or by solid solution of the tracer in the carrier. By slow cooling of liquid metallic solutions it is relatively easy to form large, regular crystals of the carrier; the effects of adsorption and occlusion may be minimized and the influence of substitutional solid-solution formation may be made to predominate. Previous studies of the crystallization of some intermetallic compounds4 and of copper, silver, and gold5 have shown that the coprecipitation or partition behavior of trace metals is entirely analogous to that observed in aqueous systems when ionic substitution occurs, e.g., coprecipitation of Ra++ with BaSo4. The precipitation of CeCd11 from liquid cadmium was chosen for study. This choice was influenced by the following consideration. If the relative atomic sizes of the tracer atom and the atom for which it substitutes determine the tendency to coprecipi-tate, then a lattice having a substitutional site of well-defined size would be desirable. In CeCd116 (structure type cubic BaHg117) the cerium atom is completely surrounded by a cage of twelve nearly equidistant cadmium atoms. Thus, atoms much larger or much smaller than cerium should be poorly accommodated in this cage. The choice of a lanthanide central atom allowed the testing of the size hypothesis with a number of lanthanide tracers which differed in size but not in valence or electro-negativity. EXPERIMENTAL The general technique of this investigation was to prepare a slightly undersaturated solution of carrier in liquid cadmium to which a tracer was added. The concentration of tracer was always so low that saturation over the entire temperature range investigated was never exceeded. The solution was cooled in a stepwise fashion. After each step a por-

Jan 1, 1964

By Pablo Hadzeriga

Great Salt Lake brine was subjected, in laboratory scale, to conditions simulating solar evaporation. Solid phases and the variation in composition of the liquid phase throughout the potassium salts crystallization zone are presented. These experimental results are discussed and compare with stable equilibrium data of the quinary system Cl-, SO4, Na+, K+, Mg + + and H2O. The Great Salt Lake brine in Utah is a vast reserve of several potential chemicals which has attracted much attention in the recent past. From the standpoint of an ore reserve, it may be said that the lake contains a low grade ore which has to be concentrated by evaporation of water. Because of the fairly favorable weather conditions in the Great Salt Lake area, it seems reasonable to consider solar evaporation as an economical means of concentration. From the view of phase chemistry, the main constituents of Great Salt Lake brine make it belong to the five component system: C1-, SO4=, Na+, K+, Mg++ and H2O. Because of their low concentration, the presence of other elements may be regarded as having little effect on the equilibrium of the five main components. Equilibrium data for the quinary system of our interest has been available for many years. Using these data, it is possible to calculate how the brine of the lake will concentrate when subjected to solar evaporation. In practicing solar evaporation, the weather factor and supersaturation conditions cannot be ignored, but their introduction into calculations using equilibrium diagrams is almost an impossible task. With the experimental work presented in this paper, an attempt is made to compare the stable equilibrium data of the quinary system as applied to the case of the Great Salt Lake brine and the dynamic equilibrium when that brine is subjected to continuous solar evaporation. Sea water brine belongs to this quinary system when only the major constituents are considered. The study of this system has been carried on since the last decades of the nineteenth century. The reason for these intensive studies was to give a geological explanation to the underground potash deposits. AS a result, many theories were presented trying to relate the equilibrium data with the geological formation, paragenesis and metamorphism of these deposits. Two main approaches were followed for the equilibrium studies of this quinary system. One of them, and the best known to us, was the school started in Germany in the last century by J. H. Van't Hoff. The other one, which started about the same time, was the school of N. S. Kurnakov in Russia. Both approaches were different and as a consequence generated opposite theories about the geological formation of potash deposits. The Van't Hoff school studied the quinary system of sea water under careful laboratory control, resulting in stable phase diagrams. During these studies, Van't Hoff recognized the difficulties in crystallizing some double salts.' ". . . . we unexpectedly found that some compounds whose formation in the solutions under investigation was possible at 25°C, nevertheless totally failed to put in an appearance. These, not to mention the calcium salts, were leonite, MgSO4 . K2SO4 4 H2O; kainite, KC1-MgSO4 . 3 H2O; and kieserite, MgSO4 .H2O. Even exceedingly slow crystallization with the addition of the compounds themselves did not remove the condition of supersaturation in the case of these bodies." However, those salts (leonite, kainite, and kieserite) were included in the stable phase diagram at 25°C as shown in Fig. 1 using D'Ans data.3 The figure also shows the sequence of crystallization of the Great Salt Lake brine according to the stable diagram. In later studies by the Van't Hoff school, especially by D'Ans,3 Karsten,4 Autenrieth,5,6 attention

Jan 1, 1968

By A. Gittens

In two recent papers stevens13,14 has compared the methods used by various investigators to calculate the contribution of grain boundary sliding to the total creep strain. In previous work Gittins15 has shown that the strain due to grain boundary sliding Egb may be estimated from the lateral, d, and vertical, v, displacements of a longitudinal marker line as it crosses a grain boundary. He obtained the formula, where nL is the number of grains per unit length measured on the specimens before creep, and \p and ?2 are, respectively, the angle between the grain boundary and the stress axis in the plane of the surface and normal to the surface (see Ref. 13, Fig. 1). The subscript L indicates that observations are made in a longitudinal direction parallel to the stress axis. It was generally found more convenient to measure the grain size of a specimen before creep; this was necessary to characterize the specimen. Stevens13,14 however prefers to define Egb in terms of NL, the number of grains per unit length after creep. Since nL - NL(l + €1) where cr is the total strain, Eq. 1 11 becomes The quantities d and ?, can be determined directly from the marker lines, while v and ?2 are found by sectioning normal to the surface. In the interior when symmetry considerations suggest that v = d and ?1 - ?2 the strain due to grain boundary sliding may be estimated from a suitable internal longitudinal marker line by the equation To determine cgb at the surface using either Eqs. |1| or [2] requires a large number of observations, particularly if many estimates of cgb are required, and also it is impossible to determine the variation of cgh with creep strain on the same specimen because the specimen must be sectioned to determine ?2. Gittins15 found that for a number of materials the sliding vectors d and v were related to each other so that provided one was measured the other could be estimated. Also during creep, grain boundaries tend to migrate to meet the surface at angles close to 90 deg,16 so that the term nL(v/tan?2)L was comparatively small when 4,- 90 deg. By determining the frequency distributions of ?2 on control specimens this term can be estimated without sectioning the creep specimen. In this way an empirical equation based on Eqs. |1| and |2| enables cgb to be determined with reasonable accuracy with fewer observations of the angles and sliding components. For B brass, the relationship between the sliding vectors indicated that ceb = 0.85nLVL . In Table I of Stevens paper this equation is stated to be incorrectly derived. However it should be pointed out that this is an empirical equation derived from Eq. |1|, based on a knowledge of the relationship between d and v and on the frequency distributions of ?1 and ?2 for a particular material under given experimental conditions. It was used only when these relationships were known, to reduce the number of observations required to estimate cb. Therefore the comment "derivation incorrect does not apply. Gittins did in fact make use of Eqs. |1| and 131. Author's Reply I must apologize to Dr. Gittins for the inadvertent misrepresentation of his work. Other slight errors in Table I have been pointed out to me. Langdon (No. 12 Table I), Gifkins (No. 7, Table I), and possibly other workers did not measure the number of boundaries per unit length after creep (A') but the corresponding quantity bejove creep (n). As Dr. Gittins points out, a completely analogous set of equations to those given in the paper can be derived by the same method: using the number of boundaries per unit length before instead of after creep. The basic equations corresponding to Eqs. [2] and 171 of the paper, from which all the others can be derived, are

Jan 1, 1968

By N. J. Grant, B. C. Giessen, U. Wolff

With the aid of the splat-cooling technique of rapid quenching, metastable alloy phases based on antimony ad bismuth have been prepared. At room temperature, simple cubic phases were found in the Sb + Au, 56 + Yd, and Sb + Ni systems. Another phase was retained in the Au-Bi system at -190°C; in this systent, a microcrystalline phase and a rhombohedrally distorted phase were also observed. The crystal chemistry of these phases is discussed in terms of their valence electron concentration range. A remarkable five-step Phase sequence was found upon heating and equilibration of certain metastable Au-Bi alloys. USING the splat- cooling technique of rapid quenching from the liquid,' a large number of metastable intermediate alloy phases have been prepared in recent years. Thus, Hume-Rothery phases of the hexagonal type based on the valence electron concentration have been produced in systems where they are missing in equilibrium, e.g., Ag-Ge 1 or AU-Sb. 3 In a comprehensive study of metastable alloys based on tin, simple hexagonal metastable y phases were found, raising the total of stable and metastable y phases from 3 to 13. 4 In the In-Bi system, five new metastable intermediate phases were retained.' In most cases, the new metastable phases have some crystal chemical characteristics in common: a) They have rather simple, elementlike structures with few atoms per unit cell, e.g., y phase,4'5 or the metastable primitive cubic phases in the Ag-Te and and Au-Te systems.= b) Their structures are often related to those of the terminal phases, or to elements that are straddled in their position in the periodic table by the constituent elements; this is the case for metastable ß (In-Bi) with the structure of white ß Sn,' or metastable a (Cd-sn). 4 c) They seem to represent the crystal structures appropriate to certain average group number (AGN) ranges in which no or few single-phase structures exist in equilibrium. These average group number ranges can probably be equated in the case of many elements of the Ba(Zn)-Bs(As) groups with the averaged valence electron concentration. In the present work, the splat-cooling technique has been applied to the alloys of antimony and bismuth. In equilibrium, neither element forms intermediate alloy phases with other B elements from groups Bl(Cu) to B4(Si) that contain more than 66.7 at. pct of Sb or Bi, respectively, e.g., AUSb 2. 7 Thus, antimony- or bismuth-rich phases, which might be closely related to antimony or bismuth, and which would reflect struc- tural changes occurring on alloying with small amounts of addition elements, do not exist. Also, the terminal solubility of antimony and bismuth for all B metals is very small. It was therefore of interest to investigate the influence of enforced, nonequilibrium solid solution of addition elements in antimony and bismuth on their crystal structures. It was found that additions of group B 2-B4 metals act distinctly differently than additions of group T 10 and B1 metals. The latter set of additives will be treated here. EXPERIMENTAL TECHNIQUES Alloy preparation and splat-quenching followed closely the procedures described in Refs. 1, 4, and 5; they consisted of the preparation of master alloys from high-purity starting materials (99.99+ pct) and blast-atomizing of 20-mg quantities onto copper sheets. The resulting approximately 1 to 5 u thick foils were investigated by X-ray diffraction. The cooling rates obtained are of the order of 107 to 10 8 C per set,'" causing solidification to proceed at a rat; of approximately 100 cm per sec.' Powder patterns of the resulting fine-crystalline materials were taken with filtered CuK, radiation. For antimony alloys the substrates were held at room temperature and for bismuth alloys at -190°C both during splat foil deposition and X-ray investigation. Annealing treatments at several

Jan 1, 1969

By J. W. Downey, R. A. Morris, M. V. Nevitt

An experimental survey has determined the occurrence of Ti2 Ni-type phases in binary and ternary alloys of Ti, Zr, or Hf with first, second, or third long-period transition metals and also in ternary systems in which 0 is the third component. The existence of eleven previously unreported isomorphs is cited and some approximate phase boundaries are given. Lattice parameters are reported for the new phases. The periodic table relationships involved in the occurrence of the phase are discussed. A previous paper' reported the occurrence of a group of Ti2Ni-type phases (structure-type E93) containing Zr in combination with Rh, Ir, and Pt, and Hf in combination with Rh, Pd, Ir, and Pt, and cited evidence to show that all of the Zr phases and the Hf-Pd phase are not binary phases but occur only when 0 is present. More recently one of the present authors found that in the ternary systems of Ti and 0 with Cr, Mn, Fe, Co, and Ni the occurrence of the Ti2Ni-type phase could be correlated qualitatively with valency electron concentration, and that the role of 0 appeared to be that of an electron acceptor.' This paper describes the results of a broader investigation of the occurrence of Ti2Ni-type phases in binary and ternary systems containing Ti, Zr, or Hf. EXPERIMENTAL PROCEDURE Alloys were prepared by arc melting. The purity of the metals Zr, Rh, Hf, Ir, Pd, and Pt was given in the previous paper.' The additional materials used for this work were metals with a minimum purity of 99.9 pct and reagent grade ZrO2, TiO2, Cr2O3, and NiO. Chemical analysis of several alloys in each system indicated good agreement between intended compositions and cast-button compositions. For the heat treatment precautions were taken to avoid oxygen contamination. The specimens were wrapped in Mo foil and sealed in fused quartz tubes at a room temperature pressure of 1 x mm of Hg. These capsules were then annealed in a vacuum furnace at an external pressure of 5 X mm of Hg. Since the major constituent in most of the alloys studied was the brittle Ti2Ni-type phase, Debye-Scherrer patterns with sharp lines could be obtained with powder crushed from the annealed solid samples without being reannealed. Metallographic, X-ray diffraction, and density determination techniques previously described1, 2 were used. RESULTS AND DISCUSSION Zr-Rh-O, Zr-Ir-Oand Zr-Pt-O Systems-Figs. 1 through 3 show the approximate boundaries of the Ti2Ni-type phase fields at the temperatures indicated in the ternary systems Zr-Rh-O, Zr-Ir-O and Zr-Pt-0. In each isothermal section the Ti2Ni-type phase field is labeled E and the unidentified binary phase found in the previous work' is labeled ?. The dotted lines, which are not phase boundaries, correspond to the composition (Zr2 B)1-xOx. In each system the Ti2Ni-type phase field is displaced to the high-Zr side of this composition line. However, the departures from Zr2B stoichiometry for the metal components are not as high as Zr3B, as was inferred from the previous work. Lattice parameters of the Ti2Ni-type phases in the alloys shown as triangular data points in Figs. 1 through 3 are given in Table IA. The probable compositions of the Ti2Ni-type phases have been estimated from the isothermal sections. Although there are a few exceptions that probably result from uncertainties in phase boundary location and from the effect of variation of oxygen content, the generalization can be

Jan 1, 1961

By L. A. Roe, E. A. Elevatorski

Wollastonite, named after William H. Wollaston, an English chemist, is a calcium metasilicate, CaSiO3; CaO: 48.30%, SiO2: 51.70%. It has a short history as an industrial mineral. The earliest production of wollastonite is reported to be from a deposit near Code Siding, located north of Randsburg, CA. At this locality small tonnages of wollastonite were quarried during 1933-34 and 1938-41, and processed into mineral wool. This operation was largely experimental and virtually no United States production was again reported until the 1950s when a large deposit near Willsboro, NY, was developed by the Cabot Corp. A processing plant was placed onstream in 1953, with nearly continuous production to date. It is currently operated by NYCO, a division of Processed Minerals, Inc. Since 1958, wollastonite deposits in the Little and Big Maria Mountains of Riverside County, and in the Panamint Range of Inyo County, both in California, have operated intermittently for production of both ornamental and commercial wollastonite. During 1980, the United States was the major producing country, furnishing about 75% of the world's output. Current production comes from Finland, Mexico, India, and Kenya. Small amounts have been shipped intermittently from the USSR, New Zealand, Republic of the Sudan, Republic of South Africa, and Namibia (South-West Africa). The principal use of wollastonite is in the manufacture of plastics. Other uses are for paints, ceramics, adhesives, fluxes, glazes, thermal insulation board, and refractory products. Mineralogy Pure wollastonite, CaSiO3, has the composition of 48.3% CaO and 51.7% SiO2. However, it is seldom found in the pure state due to the ease with which it takes into solution the metasilicates of manganese, magnesium, iron, and strontium. Predominantly, wollastonite occurs as a contact metamorphic deposit forming between limestones and igneous rocks. Commonly associated minerals are garnet, diopside, epidote, calcite, and quartz. It has a specific gravity of 2.8 to 3.0, and hardness of 4.5 to 5 on Mohs' scale. When pure, it has a brilliant white color, but with impurities it may be grayish or brownish. Luster is vitreous to pearly. Melting point of wollastonite is about 1540ºC. Wollastonite occurs in coarse-bladed masses, rarely showing good crystal form. It is usually acicular or fibrous, even in the smallest of particles. The most unique property of crushed and ground wollastonite is its cleavage. Fragments of crushed wollastonite tend to be needle-shaped, imparting a high strength, and this property is the basis for many of its uses. The fiber lengths are commonly in the ratio of 7 or 8 to 1, length to diameter. The average diameter of wollastonite is 3.5 µm. Some crystals of wollastonite fluoresce under short-wave or longwave ultraviolet light, or both, with colors ranging from yellow-orange to pink-orange. Specimens may also show phosphorescence. The acicularity of wollastonite is a property of considerable importance to the marketplace. The plastics industry makes utilization of high aspect ratio grades of wollastonite (20:1) for reinforcing thermoplastics and thermoset polymer compounds. The naturally high pH of 9.9 (10% water slurry) is a prime property to the coatings industry. The coatings industry uses milled grades of wollastonite as a pH stabilizer in interior and exterior PVA and acrylic latex systems. Processed wollastonite can have a G.E. brightness of 90 to 93. Chemically, wollastonite is inert and this property makes it useful as a filler and reinforcing agent. There are two polymorphs of calcium silicates: wollastonite, a low temperature form, and pseudowollastonite, a high temperature

Jan 1, 1983

By L. Burris, C. H. Stockman, I. G. Dillion

Zone melting is a purification process in which separation of impurities is effected by slowly moving a narrow melted zone through a bar of solid material. Equations are presented which 1—predict the concentration profiles after successive passes of the zone and 2—give the limiting distribution of solute along the bar after an infinite number of passes of the molten zone. These equations take into account the effect of normal freezing in the last zone length. These equations were solved numerically using an electronic computer, and the concentration profiles are presented as a function of the solute distribution coefficient and the ratio of molten zone to total bar length. ZONE MELTING is a purification process in which separation of impurities is effected by slowly moving a narrow melted zone through a bar of solid material. Separation occurs if the solid freezing from a solution at any moment has a composition different from that of the liquid from which it is freezing. Solutes which lower the freezing point will tend to remain in solution and will be moved in the direction of zone travel. Solutes which raise the freezing point will preferentially freeze out, and they will be moved in a direction opposite to that of zone travel. After repeated unidirectional passage of a molten zone through a bar of material, impurities will be concentrated in the ends of the bar. A limiting concentration gradient will be reached in the bar after a large number of passes of the molten zone. Both the magnitude of this limiting concentration gradient and the rate of movement of solute in any pass depend on the zone length and on the equilibrium concentrations of solute in the liquid and solid phases. The basic theory of zone melting was thoroughly covered by Pfann.1 The purpose of this paper is to extend the theory of zone melting by presenting 1—equations which relate the concentration profiles in successive passes of the zone and 2—equations which give the limiting distribution of solute along the bar after an infinite number of passes of the molten zone. These equations take into account the effect of normal freezing* in the last zone length concentration of a given solute in the solid to the concentration of the same solute in the liquid. Analytical methods for determining the concentration profiles after any pass have been developed by Lord2 and Reiss. The method developed by Lord does not take into account normal freezing in the last zone and also becomes very complicated when the number of passes becomes large. The method developed by Reiss enables the direct calculation of the complete concentration profiles (including that in the last zone) after any pass for values of k close to unity (0.9 < k < 1.1). Both methods are useful and provide valuable insight into the zone melting process. The equations given in this article are based on the same assumptions used by Pfann. These are: 1—uniform composition in the liquid, 2—negligible diffusion in the solid, 3—constant distribution coefficient, 4—constant length of the liquid zone (except at the ends), and 5—solubility of the solute in the melt is not exceeded at any point. Equations for Concentrations in Successive Passes For a melted zone of length 1 at the left-hand end of a bar of uniform unit cross-section, the amount of solute in the liquid is J0Cn-1(x) dx [1] where the subscript is the number of the pass. Cn-1(x) represents the concentration (expressed as weight per unit volume) obtained in the previous

Jan 1, 1956

By M. R. J. Wyllie, H. W. Patnode

In the quantitative interpretation of electric log data it is essential to know the formation factor, defined as the ratio of the resistivity of the formation 100 per cent saturated with brine to the resistivity of the brine.&apos; This concept implies that the formation factor for any given rock is constant irrespective of the resistivity of the brine or solution with which it is saturated. It is apparent that this can be true only when the rock matrix is wholly non-conductive. As reservoir rocks may contain clay or shale having an appreciable conductance one may expect these constituents to influence the determination of true formation factor. The formation factor may be theoretically derived from log data by dividing the resistivity of the mud filtrate invaded formation by the resistivity of the mud filtrate. Such a method of calculntion has implicit in it the assumption that the fluids in the invaded formation are wholly replaced by mud filtrate. that the solid matrix is non-conducting, and that the formation factor of the invaded zone is identical with the formation factor of the uninvaded zone. Since the validity of the first two basic assumptions regarding the insulating character of porous matrices and the displacement of the interstitial fluids in porous media do not appear to be fully established, the formation factors of a number of cores and clay slurries have been measured under conditions suitable for verifying the assumptions and the extent of any errors inherent in them. CORE STUDIES Apparatus and Procedure Resistivity measurements were made on short cylindrical permeable cores which had been previously mounted in lucite.2 The cores were clamped by means of insulated bolts between brass plates about 5 cm square and 1 cm in thickness. The plates were provided with a ring seal to prevent leakage between them and the core mount, and had a centrally located orifice so that fluids could be flowed through the mounted core. The brass mounting plates were used as electrodes and any current flow between them necessarily passed only through the constituents of the intervening core. Fluid saturation of the core, and fluid flow through the core, were effected by a centrifugal pump connected to the outflow end of the core, thus providing a differential pressure across the core of approximately one atmosphere. The procedure used was first to evacuate each core by means of the pump and then to saturate the core with an aqueous sodium chloride solution of known resistivity. The resistance of the core was measured and flow of the solution through the core continued for sufficient time to establish a constant resistance reading. After equilibrium had been reached with one solution, another sodium chloride solution of a different resistivity was flowed through the core and the flow continued until a constant resistance reading was again obtained. This procedure was continued until resistance measurements of the core had been made with a series of solutions within a salinity range from 1.0N NaCl to distilled water. The data obtained enabled the apparent formation factor of given porous rocks when saturated with fluids of different resistivity to be computed and also provided information regarding the ability of an invading fluid to displace interstitial water from a core. Results Berea Sandstone — The first core examined was a specimen of Berea sandstone." This core had been cut essentially parallel to the bedding for permeability studies. The core was clamped in the core holder and saturated with 1.0N NaCl solution and the resistance determined. It was subsequently flushed and the resistance measured when saturated in turn with 0.1N, 0.025N, 0.01N, .005N sodium chloride solutions and finally with distilled water. In flushing the core with each

Jan 1, 1950

By James H. Boettcher, W. Durand Eppler

INTRODUCTION In the second quarter of 1979, Rapid Mining Co., Inc. ("Rapid") inquired whether Crocker National Bank ("CNB") would provide $4 million in term loan financing to help ("Dry Creek") acquire the outstanding shares of the Dry Creek Coal Company ("Dry Creek") in West Virginia. Prior to approaching CNB, Rapid had sought financing with the assistance of a financial advisor from a number of commercial banks without success. A summary of the initial terms requested by Rapid is shown in Table 1 compared with the actual terms on which the transaction closed. It is not difficult to note that the two sets of terms differ significantly. This is often the case when financing requests based on the expectations and vision of a mine developer are juxtaposed with more pragmatic requirements of a commercial bank. The "art" and fund of structuring financing plans for mining ventures (highly leveraged acquisitions or grass roots) comes in working with mine developers to arrive at a final financial structure that meets their objectives and constraints as well as those of the bank or other financial institutions providing much of the financing. Beginning with the terms requested by Rapid and ending with the terms that closed the transaction, this article briefly describes the process we went through with Rapid Mining. The Seller and Buyer The Dry Creek Coal Company was established in 1961, after its owner acquired the coal leases from Republic Steel. Republic has mined the primary coal seam from 1948-1961. Dry Creek was owned as a subchapter S corporation (i.e., for federal income tax purposes, income is taxed at the shareholder level) and had been run on a low key basis. Since the company had virtually no debt and limited cash needs, the owner produced primarily when rail cars were available and income was required. Output was therefore erratic as is shown in the top of Table II, but the result was a profitable little company with an attractive free cash position. Because the owner wished to retire and enjoy his wealth in a more readily spendable form, he had put the mine, with stated recoverable reserves of approximately 13 million tons, up for sale. In contrast, Rapid Coal&apos;s investors consisted of a number of investors of disparate financial situations. They had agreed upon a primary objective of income maximization and planned an annual production rate of 120,000 TPA which they and the current owner believed realizable. The Formulation of a Financing Structure Figure 1 presents a simplified overview of the sequence of steps taken in the financing process which lasted approximately three months from start to finish. Since banks see many requests for funding and considerble time is required to follow up on each, an important step was determining, after a brief review of the materials provided to us, whether the basic "eonomic entity" was viable. In this case, the size and quality of the reserves, the inherent operating cost characteristics, the sales arrangements, the management/labor situation and Rapid&apos;s principles all appeared to be sound enough to warrant investing the time necessary to do a thorough analysis of the project. Once this decision was made, the evaluation and analytic process shown in Figure 2 was initiated. The first and most fundamental concern of any banker is whether the loan will be paid back as scheduled from the cash flow generated by the project. The amount of debt relative to the owner&apos;s equity is of course crucial in such a determination. Before the debt/equity issue was addressed, however, a number of more qualitative factors had to be studied and if

Jan 1, 1985

By J. O. Brittain, E. Lautenschlager

Yield point measurements were made on polycrys-talline specimens of a iron deformed various amounts at 770, 1950, or 273°K, aged, and redeformed at 77: 195°, or 273°K. When the prestrain temperature is lower than that subsequent to aging, the yield point return is suppressed relative to specimens deformed at the same temperature both before and after aging. The degree of suppression is a function of both the magnitude of prestrain at the lower temperature and the difference between the deformation temperatures before and after aging. It is proposed that the suppression is caused by the formation of dislocation barriers during the initial low-temperature deformation. The aging kinetics of polycrystalline a-iron wires deformed at 77o or 303oK were studied at 44 " and 64 "C using internal friction measurements. No change was noted as a result of the difference in prestrain temperature. IN an earlier investigation1 it was observed that the return of the yield point was suppressed when a poly-crystalline specimen of a iron was plastically deformed at 77oK, aged at room temperature, and subsequently tested at room temperature. In order to investigate this observation, a research program was undertaken to measure the effect of prestrain temperature on the strain aging of a iron.&apos; According to Thomas and Leak,3 return of the yield point and internal friction are similar measures of strain aging since both are dependent on solute atoms leaving the solid solution. The decrease in internal friction as a function of aging time is indicative of the removal of solute atoms from the solution, while the return of the yield point as a function of aging time gives the rate of arrival of solute atoms at the dislocations. According to Cottrell and Bilby,4 the solute atoms create an "atmosphere" which locks the dislocations and causes a return of the yield point upon restraining. In this experiment, both methods of measuring strain aging were employed in order to have two independent measurements of the strain aging phenomenon. In a previous study,5 the strain aging kinetics of a similar lot of a iron were measured for specimens which were prestrained and tested at the same temperature over the range of 77" to 273°K. Although some differences were found in the correlation of the two methods employed for measuring strain aging, the kinetics of strain aging were unaffected by the change of prestrain and test temperature. In another study,6 dislocation line tension and the density of the solute atmosphere were found to contribute to the return of the yield point when ingot iron specimens were prestrained at 303°K and tested at a lower temperature, namely 77°K. This paper is concerned with changes in strain aging kinetics produced by pre-straining at 77" to 303°K and testing at either the prestrain temperature or at higher temperatures up to 273°K. EXPERIMENTAL PROCEDURE The polycrystalline a iron specimens used in both the tensile and internal friction tests were 0.033-in. diam ingot iron wires, wet hydrogen treated and re-carburized. A more detailed procedure for processing a 3/16 in. ingot iron rod (nominal wt. pet composition: Fe - 0.013; C - 0.0125; Mn - 0.023; Si -0.005; P - 0.0023; S - 0.11; Cu - 0.004 N) to 0.033 in. diam specimens containing approximately 0.020 wt pet C and no detectable nitrogen with a grain diameter of 0.033 to 0.040 mm is given elsewhere.5 Immediately prior to testing, specimens were solution treated for 15 min at 700o ± 2°C in a static vacuum of less than 50 ?. of Hg and rapidly quenched into an ice-water bath while still enclosed in a 3 mm quartz vial. The internal friction specimens had 8-in. gage lengths and the tensile specimens had 2-in. gage lengths. The wires used as tensile specimens were first run as internal

Jan 1, 1962

By R. A. Morse, P. M. Bridges, L. E. Wilsey, Howard N. Hall, P. L. Terwilliger

Theoretical and experimental investigations of a constant pressure gravity drainage system are reported. Experimental data are presented to show that recovery to gas breakthrough by gravity drainage is inversely proportional to rate. The gravity drainage reference rate, which is numerically equal to the so-called "maximum theoretical rate of gravity drainage" is shown to have no particular significance from a recovery standpoint. Before this rate can be used as a basis of comparison for recoveries, it is necessary that the relative permeability and capillary pressure characteristics and displacing fluid viscosities be identical for the systems compared. A method is presented by which accurate prediction of the performance of a gravity drainage system can be made. Close agreement between experimental and calculated drainage performance shows that steady state relative permeability and static capillary pressure data can be used to describe fluid displacement behavior. The very wide range of liquid recoveries before gas breakthrough which result from production rate variation alone demonstrates the importance of this factor in planning depletion of a gravity drainage reservoir. Calculated results are presented which show that little additional recovery can be expected from a high pressure gravity drainage system between the times of gas breakthrough and attainment of such high gas/liquid ratios as to make further pressure maintenance impractical. INTRODUCTION It bas been recognized for some time that gravity forces play an important part in the recovery of oil from some types of reservoirs. Field experience has shown that under certain conditions, gravity drainage can result in very high oil recov- eries. Qualitative reasoning has led most engineers to the general conclusions that: (1) Where gravity drainage is important, the reservoir pressure should be maintained by gas injection at the crest of the structure to prevent shrinkage of the oil in place and to keep a low viscosity so the oil can drain at the fastest possible rate. (2) Recovery by gravity drainage is rate sensitive. A survey of the literature indicates that while considerable work has been done on the effects of gravity in oil production problems, no satisfactory method of calculating the performance of gravity drainage reservoirs has been reported.1,2,3,4,5 In the absence of any proven method of calculating reservoir performance, the level at which pressure should be maintained and the rate of production for most efficient operation has been open to debate. The limits and general character of the recovery us rate relationships are apparent from the following: As rates of production approach zero, oil will drain by virtue of its own weight fast enough that the saturation distribution is always approximately the same as the original static capillary distribution. That is, at the gas-oil contact the liquid saturation goes from 100 per cent to a value of perhaps 10 to 30 per cent in a very short distance, so that as production continues. a relatively sharp gas-oil contact moves down-structure, leaving a very low residual oil saturation in the originally oil saturated zone. At very high rates of production, the oil that drains by virtue of its own weight is negligible and recovery will approach that of a horizontal gas drive, leaving a much higher residual oil saturation. At intermediate rates of production. recovery should be between these extreme*. Despite the lack of quantitative methods, tentative agreement had been reached on the question of optimum producing rate. It has been generally concluded that rather sharp decrease in recovery would occur at rates of production above the "maximum rate of gravity drainage" and hence this rate should not be exceeded. The "maximum rate of gravity drainage" is defined as the rate of production from a 100 per cent

Jan 1, 1951

By L. L. Melton, W. T. Malone

Fluid mechanics research conducted with non-Newtonian fluid systems now permits prediction of the behavior of these fluid systems in both laminar and turbulent modes of flow through circular pipes. Present work concerns non-Newtonian fluid systems currently used in the hydraulic fracturing process. During fracturing treatments, an unsteady-state condition may frequently be encountered arising from. the reaction rate of a chemical additive. This condition must be evaluated in order to predict the actual behavior of a particular fluid during field application. Design and operation of the apparatus used to determine fluid-flow behavior permit obtaining data under such non-equilibrium conditions. This paper shows methods used to obtain rheology measurements, develop hydraulic relationships and evaluate chemical reactions producing unsteady-state conditions. Engineering application of this research is illustrated by employing measured rheological values and developed hydraulic retationships to produce frictional pressure loss (psi/100 ft) us flow rate (bbl/min) charts of common tubing and casing sizes for an example fracturing fluid. How these charts and chemical reaction rate information are then combined to predict actual turbulent hydraulic behavior during unsteady-state field conditions is also discussed INTRODUCTION Many fluids used in hydraulic fracturing contain chemical additives which impart non-Newtonian fluid properties that may drastically alter their hydraulic behavior. Equally drastic alteration in wellhead pressure, injection rate and hydraulic horsepower requirement may result from these fluid properties. Prior research conducted to relate non-Newtonian fluid properties with hydraulic behavior has not as yet produced a universal relationship, particularly for the turbulent flow region. which of the many Possible non-Newtonian fluid properties is responsible has not been conclusively established. A systematic description, suggested by Metzner, of the many possible non-Newtonian fluid properties exhibited by real-fluid behavior, and a current discussion of theoretical and applied aspects of non-Newtonian fluid technoloy can be found in Handbook of Fluid Dynamics. f Little or no research has previously been attempted with fluids exhibiting time - dependent properties. Addition of chemicals during a fracturing treatment is often accomplished by continuously mixing and displacing the fluid. This produces a time-dependent effect or unsteady-state condition while the fluid is progressing through surface and wellbore conductors. This condition is due to solution or chemical reaction of the additive. Considerable departure from conventional methods of obtaining and interpreting data was found necessary to consider these conditions. Therefore methods were developed to obtain hydraulic behavior information under these complex, unsteady-state conditions. Relationships presented in this paper to predict hydraulic behavior of non-Newtonian fluids in circular pipes were obtained by constructing and operating a small pipeline apparatus in the manner of a pilot-plant study. These relationships are suggested as scaleup equations and are not proposed as fundamental rheological parameters. While perhaps deficient from a fundamental research viewpoint, a pilot-plant study does permit the determination and convenient evaluation of variables pertinent to a process. A pilot-plant study can result in a valid engineering application procedure even when fundamental relationships are still undefined. An excellent series of articles by Bowen has appeared in the chemical engineering literature.2-8 These give a thorough review of proposed hydraulic relationships and their limitations for non-Newtonian fluid behavior in circular pipes. A graphical method is presented to scale up data for a fluid exhibiting an anomalous hydrau1ic behavior in the turbulent flow region. Considerable assistance has been obtained from these articles to interpret the anomalous behavior noted during this investigation. These articles also provided assurance that a pilot plant is a practical approach to evaluate the hydraulic behavior of non-Newtonian fracturing fluids. _____________

Jan 1, 1965

By K. H. Coats

This paper presents the development and solution of a mathematical model for aquifer water movement about bottom-water-drive reservoirs. Pressure gradients in the vertical direction due to router flow are taken into account. A vertical permeability equal to a fraction of the horizontal permeability is also included in the model. The solution is given in the form of a dimensionless pressure-drop quantity tabulated as a function of dimensionless time. This quantity can be used in given equations to compute reservoir pressure from a known water-influx rate, to predict water-influx rate (or cumulative amount) from a reservoir-pressure schedule or to predict gas reservoir pressure and pore-volume performance from a given gas-in-place schedule. The model is applied in example problems to gas-storage reservoirs, and the difference between reservoir performances predicted by the thick sand model of this paper and the horizontal, radial-flow model is shown to be appreciable. INTRODUCTION The calculation of aquifer water movement into or out of oil and gas reservoirs situated on aquifers is important in pressure maintenance studies, material-balance and well-flooding calculations. In gas storage operations, a knowledge of the water movemenr is especially important in predicting pressure and pore-volume behavior. Throughout this paper the term "pore volume" denotes volume occupied by the reservoir fluid, while the term "flow model" refers to the idealized or mathematical representation of water flow in the reservoir-aquifer system. The prediction of water movement requires selection of a flow model for the reservoir-aquifer system. A physically reasonable flow model treated in detail to date is the radial-flow model considered by van Everdingen and Hurst.1 In many cases the reservoir is situated on top of the aquifer with a continuous horizontal interface between reservoir fluid and aquifer water and with a significant depth of aquifer underlying the reservoir. In these cases, bottom-water drive will occur, and a three-dimensional model accounting for the pressure gradient and water flow in the vertical direction should be employed. This paper treats such a model in detail — from the description of the model through formulation of the governing partial differential equation to solution of the equation and preparation of tables giving dimensionless pressure drop as a function of dimensionless time. The model rigorously accounts for the practical case of a vertical permeability equal to some fraction of the horizontal permeability. The pressure-drop values can be used in given equations to predict reservoir pressure from a known water-influx rate or to predict water-influx rate (or cumulative amount) when the reservoir pressure is known. The inclusion of gravity in this analysis is actually trivial since gravity has virtually no efFect on the flow of a homogeneous, slightly compressible fluid in a fixed-boundary system subject to the boundary conditions imposed in this study. Thus, if the acceleration of gravity is set equal to zero in the following equations, the final result is unchanged. The pressure distribution is altered by inclusion of gravity in the analysis, but only by the time-constant hydrostatic head. The equations developed are applied in an example case study to predict the pressure and pore-volume behavior of a gas storage reservoir. The prediction of reservoir performance based on the bottom-water-drive model is shown to differ significantly from that based on van Everdingen and Hurst&apos;s horizontal-flow model. DESCRIPTION OF FLOW MODEL The edge-water-drive flow model treated by van Everdingen and Hurst1 is shown in Fig. la. The aquifer thickness b is small in relation to reservoir radius rb. water invades or recedes from the field at the latter&apos;s edges, and only horizontal radial flow is considered as shown in Fig. 1b. The bottom-water-drive reservoir-aquifer system treated herein is sketched in Fig. 2a and 2b. Here the aquifer

By B. W. Levinger

THE system Ti-CI- has been studied by several investigators." &apos; Though titanium and chromium are completely miscible at high temperatures, an intermediate phase of the approximate atomic proportions TiCr, has been identified at various compositions between 60 and 68 wt pct Cr. The structure of this phase was established as the cubic C15 (MgCu2) type with 24 atoms per unit cell and a lattice parameter 6.929 kX.2 During a study of the ternary system, Ti-Mo-Cr,5 the existence of a high temperature modification of the TiCr, phase was discovered. In a series of alloys annealed at 1300°C and containing from 60 to 70 wt pct Cr and 5 pct Mo, a new phase was present which did not correspond to either ß solid solution or TiCr,. By trial, the new phase was identified to be of the hexagonal C14 (MgZn2) structure with 12 atoms per unit cell. Further work showed that the hexagonal structure was. in fact, also present in binary Ti-Cr alloys quenched from high temperatures. To make the identification more certain, the intensities of the various reflections were calculated assuming a composition TiCr2 and C14 structure. Table I lists reflection planes, interplanar distances, calculated intensities, and a comparison of these values with the lines of a binary Ti-Cr alloy (66 pct Cr) annealed for 20 min at 1385°C and water quenched. Since high angle lines were more diffuse, the correspondence between visually observed and calculated intensities was not as good at these angles. Lattice parameters were computed to be a = 4.922 kX, c = 7.945, c/a :- 1.614. The corresponding interatomic distances are: titanium to titanium distance, 3.00:, kx, chromium to chromium distance, 2.43, kX, titanium to chromium distance, 2.88,, kX. These compare to 3.00, 2.45, and 2.87 kX in the cubic phase, respectively.&apos; The atomic positions in the C14 structure are as follows: Titanium 1/3, 2/3, Z; 1/3, 2/3, 1/2: 2/3, 1/3, i; 2/3, 1/3, (1/2-z); Chromium 0, 0, 0: 0, 0, 1/2; 2x, 2, 1/4; 232, x, %; X, 22, 1/4; X, 2x, 3/4; x. x, 1/4; x, x, 3/4; The values of x and z were calculated for the alloy mentioned and found as follows: z — 0.061, x = 0.165. The values assumed in the intensity calculations were those for an ideal close-packed lattice or z = 1/16, x = 1/6. Work in the Ti-Mo-Cr system indicates that y2 the hexagonal modification of TiCr2, is the stable phase above 1300°C, while y1 with the cubic structure is stable below 1000°C. X-ray examination of a number of alloys annealed at intermediate temperatures showed that the hexagonal and cubic phases can coexist. Alloys quenched from 1000°C and annealed at 1100" and 1200°C showed small amounts of y2 at the higher temperature. A considerable amount of y, seemed present, however, in cast alloys annealed at 1100°C. It further appears that there is a small but definite solubility range for both phases in the ternary system. Acknowledgment The results of this paper represent work carried out in conjunction with a project sponsored by Watertown Arsenal, Contract No. DA-11-022-ORD-272. References &apos;M. K. McQuillan: Journal Inst. Metals (1951) 79, p. 379. 2 P. Duwez and J. L. Taylor: Trans. ASM (1952) 44, p. 495. W. J. Van Thyne,. H.. D. Kessler, and M. Hansen: Trans. ASM (1952) 44, p. 974. 1 F. B. Cuff, N. J. Grant, and C. F. Floe: Trans. AIME (1952) 194, p. 848; Journal of Metals (August 1952). Unpublished Research, Armour Research Foundation. "H. J. Wallbaum: Ztsch. Krist. 103, p. 391.

Jan 1, 1954