Search Documents

By A.J. PHILLIPS, Wm. B. Price

ADMIRALTY nickel is a new corrosion-resisting and heat-resisting white metal alloy composed of 70 per cent. copper, 29 per cent. nickel and 1 per cent. tin. It has been given the trade name "Adnic." In 1924, Win. B. Price1 described this alloy and its successful use as a diaphragm metal. The physical properties of Adnic at ordinary and elevated temperatures were given in Mining & Metallury.2 A patent3 was granted on Adnic, April 28, 1925. Work was first started on this alloy during December, 1921. An investigation of the patent situation in this country and abroad, and also a search of the literature, did not reveal any work on the ternary diagram, copper-nickel-tin; and the only wrought copper-nickel-tin alloy found was that known as "Newloy," composed of copper, 64 per cent.; nickel, 35 per cent.; tin, 1 per cent. This was being manufactured and sold by Wm.. Gallimore & Sons, Ltd., Sheffield, England, for unplated spoon and fork stock, and first came to our attention in January, 1922. The preliminary work in connection with the development of Adnic was started by casting four, binary copper-nickel alloys of the compositions shown in Table 1. Various amounts of tin were added to each of these mixtures, to produce the series of ternary copper-nickel-tin alloys also given in Table 1. A 200-gm. sample of each of these was cast by melting in a clay graphite crucible and allowing the melt to slowly cool with the furnace. A longitudinal section from each of these castings was polished and etched for-examination. The castings were not homogeneous solid solutions, but showed the usual cored structure associated with cast alloys. However, from an examination of these sections we roughly determined

Jan 1, 1928

By D. E. Baldwin

A method has been devetoped for predicting the produced concentration profile for a miscible slug in a five-spot pattern. The technique consists of dividing the five-spot into radial etements and applying an approximate radial solution of the dispersion equation to each element. To test the method, laboratory experiments were conducted in five-spot models. The models were 2-ft X 2-ft X 2-in. Berea sandstone slabs simulating one quarter of a five-spot. An aqueous system was used with tritiated water as the tracer fluid. Effiuent concentration profiles were obtained for slug sizes ranging from 0.7 to 10.2 per cent of the pore volume. Agreement between predicted and experimental profiles was excellent. A prediction was also made for a field tracer test. For this case the predictive technique was modified to account for stratification by use of a layer cake model. INTRODUCTION Dispersion in porous media is of growing interest to the petroleum industry because of the increasing importance of secondary recovery operations. One of the controlling factors in the recovery of oil by miscible displacement or by the use of waterflood additives is the degree of mixing between the fluids of interest. Consequently, one of the prerequisites of any recovery prediction is an adequate method of predicting dispersion. Also, any quantitative interpretation of waterflood tracer tests is dependent upon the ability to determine the dispersion of the tracer as it flows through the reservoir. During the past few years there have been reported in the literature the results of a large number of investigations, both theoretical and experimental, concerning dispersion in porous media. (An excellent review of these results has been presented by Perkins and Johnston." Only a few of these investigations, however, deal with non-lin-ear flow systems. Raimondi et al.9 have presented an approximate solution of the dispersion equation for radial, diverging flow. Lau et al. 7 and Bentsen 1 have presented laboratory data for radial systems and have shown that the dispersion can be adequately predicted by Raimondi's solution. Bentsen has also presented a small amount of data for a radial, converging system. Brigham and Smith' have applied Raimondi's solution, with approximating assumptions, to a five-spot pattern to predict the behavior of a field tracer test. The afore-mentioned investigations have been concerned almost entirely with radial diverging flow and have neglected the effect of converging flow. For this reason the present investigation was undertaken to develop a method of predicting dispersion in a five-spot pattern including the effects of both diverging and converging flow. The method which was developed consisted of dividing the five-spot into diverging-converging radial elements and applying an approximate solution of the disperson equation to each element. DEVELOPMENT OF PREDICTION METHOD The development of the prediction method consisted of two phases: construction of a radial model of a five-spot pattern and derivation of equations to describe the dispersion in the radial model. FIVE-SPOT MODEL A radial model which simulated the flow characteristics of a five-spot pattern was constructed by dividing the five-spot into diverging-converging radial elements (Fig. 1). Because of symmetry the model represented only an octant of a complete five-spot. Each element had an included angle of l° (making a total of 45 elements) and a radius r, such that the sum of the individual element areas was equal to the area of the five-spot octant. The normal five-spot fractional flow behavior was determined by averaging the experimental five-spot data of Caudle and Witte,' Dyes, Caudle and Erickson 5 and Fay and Prats. 4 The resulting curve is shown as a solid line in Fig. 2. To match this curve with the radial model it was

Jan 1, 1967

By Richard H. Siller, Rex B. McLellan

A model for interstitial solid solutions has been considered in which a repulsive Potential exists between interstitial atoms in the solvent lattice. It has been shown that this model is consistent with the thermodynamic Properties of certain interstitial solid solutions. The model is then applied to the calculation of the variation of the diffusivity of carbon in au-tenite with composition. Although the energy levels at the saddle point, and thus the absolute values of dif-fusivity, cannot be calculated, it is shown that by taking the repulsive interactions between the diffusing atoms into account the functional form of the dependence of diffusivity on solute concentration can be derived. ThE composition dependence of the diffusivity, D, of carbon in austenite at a given temperature is well-known.'-3 smith3 tried to reconcile the observed composition dependence of D at 1000°C with the previously observed4 composition dependence of carbon activity using equations for D derived by arken5 and Fisher, Hollomon, and Turnbull,6 but was unsuccessful. In this paper we shall modify the approach of Fisher et al. by considering the "blocking" model7 for interstitial solid solutions, based on concept due to Speiser and Spretnak,' in which a repulsive potential exists between interstitial solute atoms in the solvent lattice, causing a certain number of interstitial sites to be "blocked" from occupancy by further solute atoms. We shall thus be able to describe the intimate relationship between the composition dependence of carbon diffusivity, a kinetic phenomenon, and the composition dependence of carbon activity in austenite. THEORY Consider for convenience one-dimensional interstitial diffusion in an fcc lattice in the [loo] direction. Fig. 1. Assume that solute atoms occupy octahedral sites only.7,9'10 These octahedral sites are located on {lOO} planes a/2Å apart. A solute atom, A, occupying an octahedral site on a {100 } plane denoted by 2, has twelve nearest-neighbor octahedral sites: four on plane 1, four on plane 2, and four on plane 3. Diffusion occurs when there is a net flux of solute atoms, through an imaginary plane shown with dashed lines, from plane 2 to plane 3. Fisher et al6. use the concept of activated complexes, as given in the theory of absolute reaction rates." They show that the number of solute atoms per cubic centimeter per second which pc — Fig. 1-Diagram showing the arrangement of octahedral sites (x) in the fcc lattice. pass in the direction of diffusion from activated to new interstitial position is: ?where C is the concentration of the solute in atoms per cubic centimeter, l/m + is the fraction of new interstitial positions corresponding to motion in the forward diffusion direction (transmission coefficient), ?F* is the free-energy difference between the activated complex and the reactants when each is in its standard state at the temperature of the reaction (independent of temperature and composition in the case of carbon in austenite), y is the activity coefficient of the solute species, and, finally, ym is the activity coefficient of the activated complex (assumed constant). The flux of atoms moving in the forward direction through an area of 1 sq cm perpendicular to the direction of diffusion is therefore: m+ Ym l ' atoms per second where is the distance between planes. The blocking model states that a nearest-neighbor octahedral site of a given solute atom has a probability, f, of being inaccessible to other solute atoms. The number of available sites for diffusion on a given plane will, therefore, depend on f, which is constant, and Nu, the number of solute atoms on that plane. The effect is to have more available sites in the diffusion direction, thereby increasing the net flux of solute atoms. This is reflected in the transmission coefficient, l/m. To calculate the transmission coefficient for solute

Jan 1, 1970

By J. Weertman

It is shown that conditions suitable for the conversion of straight dislocations into helices are common in crystals hardened either through long-range dislocation interaction or by jog formation on dislocation lines. For dislocations other than screw dislocations or dislocations nearly screw in orientation, helices can be formed with the aid of core difhsion at temperatures low relative to the melting point of the material. It is assumed that, once formed, helices will degenerate into dislocation tangles. THE dislocation tangle is one major dislocation arrangement observed by the transmission electron microscope which had not been predicted previously from dislocation theory. Since first being observed, the tangle itself has become important in applied dislocation theory. The Cambridge school (P. Hirsch) bases part of its theory of work hardening upon this phenomenon.' In spite of the fact that the tangle obviously is a structure of prime importance to a cold-worked crystal, until recently relatively little attention has been paid to the question of exactly why or how tangles arise. It was generally assumed at first that cross slip and perhaps jog formation will lead to tangles. However Wilsdorf and schmitz2 concluded that cross slip cannot markedly contribute to the tangles. The Wilsdorfs and Maddin3'* have proposed that tangles arise by means of prismatic dislocation loops and superjogs which are produced by a precipitation of point defects. Dislocations with superjogs or which have joined with prismatic loops are themselves dislocation sources and they can produce new dislocations, a process they call "mushrooming", which then tangle with each other. There is no question that the Wilsdorfs et al. mechanism as well as the cross slip one can pro- duce tangles. However these may not be the only paths leading to tangle formation. In this paper we wish to consider another way tangles can form. Our work will be based on the theory of dislocation screw We wish to show that whenever a crystal is cold-worked there always is present a "driving force" leading to helix formation and thus to tangle formation. We do not expect to observe well-developed, regular helices in cold-worked metals. Rather we would expect that after a helix is formed random interactions with nearby dislocations will cause a well-developed helical dislocation to go over into an irregular form. A dislocation tangle is, of course, a grouping of dislocations with irregular forms. Therefore we feel that, in accounting for the formation of dislocation tangles, it is sufficient to show that many dislocations in a cold-worked crystal will take on helical form. The crucial step is helix formation. In order to have a straight dislocation assume a helical form, diffusion of vacancies, or interstitials, must occur. Our discussion of tangle formation therefore will be limited to crystals which have been deformed at temperatures ranging from about one fifth to half their absolute melting point. In this tem per a ture range appreciable core diffusion is possible. (There is no problem in explaining tangle formation at temperatures greater than half the melting point. At such temperatures dislocation climb through a bulk diffusion of point defects gives an extra degree of freedom to dislocation motion. This extra freedom of motion makes tangle formation easy.) At temperatures much below one fifth the absolute melting temperature diffusion processes are sufficiently slow that helix formation does not occur. Of course helix or tangle formation could occur in material deformed at very low tempera-tures which was subsequently warmed in order to make the thin specimen necessary for electron-microscope examination. Since within the temperature range of approximately one fifth to half the melting temperature the diffusion of point defects in the core of the dislocation is rapid, we can reasonably assume that within

Jan 1, 1963

By R. W. Fountain, J. F. Libsch

ONSIDERABLE time and effort have been ex- pended recently in research designed to provide a better understanding of the solid state transformations leading to the permanent magnet qualities of many commercial alloys. The knowledge derived from these investigations has resulted in the development of new magnetic materials as well as improvements in some of the older materials. The permanent magnet qualities of the newer alloys having high residual inductions and high coercive force values, have been attributed to the precipitation of a second phase or an order-disorder transformation in a solid solution. The ternary Co-Fe-V alloys provide an interesting group of magnetic materials having a wide range of magnetic properties. The low vanadium alloys' possess soft magnetic properties similar to the binary Co-Fe alloys, while at higher vanadium contents (about 10 pct) interesting permanent magnet properties are obtained.2,3 Previous work by Nesbitt' on the ternary Co-Fe-V permanent magnet materials, resulted in the conclusion that the precipitation of a second phase, ?, was responsible for the permanent magnet qualities of these alloys. However, unlike the conventional precipitation hardening alloys, the ? phase, which is stable at high temperatures, is precipitated in the low temperature a phase upon aging. Recent work by Geisler and Martin"." has established that the order-disorder reaction inherent in the binary Co-Fe alloys is also present in the ternary Co-Fe-V alloys. On the basis of their investigation, it was proposed that ordering of the a phase is also a possible explanation of the mechanical and magnetic hardening in these alloys. Thus, it appears that the Co-Fe-V alloys are rather unique in that two solid state transformations are taking place which may account for the properties of these alloys. In view of these facts, the present investigation was undertaken to study the mechanism of mechanical and magnetic hardening in the 10 pct V-Co-Fe alloy with particular reference to the role of the ? phase; and to correlate the results of this study with other magnetic properties. In addition, the influence of increasing vanadium contents on the transformation of Co-Fe alloys was also undertaken to supplement the main objective of the investigation. Experimental Procedure The Co-Fe-V alloys used for study consisted of: 1—a commercial lot of Vicalloy I, a 10 pct V alloy, and 2—alloys prepared by powder metallurgy and containing 52 pct Co and 0 to 14 pct V, balance iron. The commercial alloy was received in the cast and forged condition and was used primarily for the detailed study of hardness, electrical resistivity, and magnetic properties of the 10 pct V alloy. The powder metallurgy alloys were prepared from electrolytic iron, commercial cobalt, and ferro-vanadium powders. The processing cycle was as follows: 1—Elemental powders were weighed and mechanically mixed for 24 hr. 2—Mixed powders were pressed at 60,000 psi into the form of bars 5/8 x 5/8 x 5 in. 3—Bars were sintered in an atmosphere of purified dry hydrogen for 24 hr at temperatures ranging from 1380" to 1340°C, depending upon the vanadium content. The as-sintered bars were cut in half; and one half was machined to 7/16 in. diam for thermal analysis. The remainder of the bar was hot forged to approximately 0.050 in. for X-ray diffraction analysis. Additional details of sample preparation are presented subsequently according to the type of test performed. Thermal Analysis Measurements One of the simplest and most widely used tools for determining phase relations as they change with temperature and composition is to study the rate of change of temperature of a material as heat is supplied or extracted at a constant rate. The thermal analysis method used in this investigation is the inverse-rate method described by Smith.'

Jan 1, 1954

By Charles B. Dudley

THEO. N. ELY, Esq., Superintendent of Motive Power. DEAR SIR: It is now nearly three years since my first report to you on the subject of steel rails was written. That report, as you will remember, dealt principally with the question of the relation between the chemical composition and physical properties of steel rails and their power to resist crushing and fracture in actual service. Other matters were referred to or touched upon in that report, but the main question was, Why do some rails crush or break in service while others do not? You will doubtless remember that the principal conclusion arrived at was, that the softer rails are less liable to crush or break in service than the hard ones. Or, in other words, so far as conclusions could be drawn from the chemical analysis and physical test of 25 samples of steel rails which had actually been in service, these conclusions were that those rails which have the smaller amounts of carbon, phosphorus, silicon, and manganese are less liable to crush or break in service than those which have larger amounts of these elements. Or, again, looked at in the light * At the Lake George meeting of the Institute in October, 1878, I had the honor of presenting to the Institute, through the kind permission of the officers of. the Pennsylvania Railroad Company, the results of a study of twenty-five samples of steel rails, which had all been in actual service. Considerable discussion followed the publication of that paper, and there seemed to be a strong disinclination, especially on the part of the steel-rail manufacturers, to accept the conclusions presented. One of the principal objections urged against the conclusions drawn was, that they were based on too few samples; in other words, that no conclusions safe to act upon could be drawn from the examination of twenty-five rails. In view of this criticism it was decided to repeat the investigation with a larger number of samples, and with the aid of the experience gained in the first investigation. The results of this second study of steel rails are, by the permission of the officers of 'the Pennsylvania. Railroad Company, herewith presented to the Institute, with the sincere desire that they may aid in adding to our knowledge of this most important product. Like the previous paper, this is in the form of a report to one of the officers of the company, which will account for the style, and the manner of presenting the data

Jan 1, 1881

By F. V. Lenel, G. S. Ansell

The creep behavior of an aluminum -aluminum oxide alloy, A T 400, fabricated by compacting an atomized aluminum powder, extruding the compact, cold working, and recrystallizing the extrusion, was investigated. This alloy has a dispersion spacing of 1 . Previously Ansell and Weertman determined the steady-state creep rate of a recrystallized aluminum-aluminum oxide SAP-type alloy, MD 2100, with a 0.48-p dispersion. This alloy was found to have a steady-state creep rate of less than 10-6 min-1, too small to be measured accurately. By employing a more sensitive method of measuring creep rate, it was possible to determine the steady-state creep rate of AT 400 alloy, which is somewhat faster than that of MD 2100, at a series of temperature and applied stresses. Above a stress of Burger's vector, A = dispersion spacing, u = is a shear modulus) the rate follms an equation of the type: Steady-state creep rate K = A exp. where A is a constant, Q is the activation energy for sey-diffusion, and is the applied stress. At stresses lower than the steady-state creep rate drops off sharply and is no longer in agreement with this equation. This behavior is in good agreement with that predicted by the model for dispersion strengthened alloys in the stress and temperature range investigated. The absolute value of the creep rate is, however, four orders of magnitude lower than that predicted if the normal number of active dislocation sources are assumed to be present. It is concluded that in this SAP-type alloy the usual three-dimensional dislocation network is not present. Instead, short dislocation segments extending from one dispersed particle and terminating at a neighboring particle may act as dislocation sources. With this provision the creep model proposed by Ansell and Weertman reasonably accounts for the creep behavior which was observed. RECENTLY, Ansell and Weertman1 proposed a dislocation model from which they derived creep equations to describe the steady-state creep behavior of dispersion - strengthened alloys. Following Schoeck,2 they assumed that the rate controlling process for steady-state creep is the climb of dislocations over the dispersed particles. In order to verify their creep model, Ansell and Weertman determined the steady-state creep behavior of an aluminum-aluminum oxide SAP-type alloy, MD-2100, in both the fine-grained as-extruded, and coarse-grained recrystallized conditions. Since the grain size was the order of the dispersion-spacing, their model was not applicable for the as-ex- truded alloy. It should, however, be applicable to the coarse-grained recrystallized alloy. The results they obtained were rather unexpected. No meas-ureable steady-state creep was observed for the recrystallized MD-2100 alloy. If their model is correct, and if the second-phase particles act solely to hinder dislocation motion, then some measureable steady-state creep would have been expected. On this basis, they postulated that, in this SAP-type alloy, the main effect of the fine dispersion was to inactivate dislocation sources rather than to hinder the movement of dislocations. In order to understand more fully this unusual creep behavior, and to determine the validity of the model proposed,' the steady-state creep behavior of an aluminum-aluminum oxide recrystallized SAP-type alloy, with a somewhat coarser dispersion spacing, was investigated. The results of this study are presented in this paper.

Jan 1, 1962

By John Johnston

WHEN I was honored by being invited to give the Howe Memorial Lecture, I decided to read Howe's book, "The Metallography of Steel and Cast Iron," published in 1916-that is, about 25 years ago-in search of a text. I found the book written in fine, clear, attractive, yet precise style, as I expected, knowing that his father Dr. Samuel Gridley Howe, his mother Julia Ward Howe, and all of his four brothers and sisters had lived in a New England literary atmosphere and had written books. This leads me to suggest that many authors of metallurgical papers and books could well study Howe's writings with a view to the improvement of their own. I failed to find a specific text, but did find a subject, for in reading his book I was struck by the few instances in which time is considered explicitly as a factor in the phenomena examined, and by the fact that many of the hottest arguments would have been settled had time been taken properly into account. Indeed, Howe hardly mentions time except in a rather incidental way, as in a few references to lag; for instance: The lowering of temperature at which the transformation occurs by rapid cooling . . . is a phenomenon of undercooling or lag. It is in part explained by the fact that the transformation itself is a time-consuming reaction which, setting in during a rapid fall of temperature, is naturally protracted till the temperature has sunk far below that at which it is due. He does discuss, quite fully, the iron-carbon equilibrium diagram, as known at that time, introducing it by stating: As he is a reasonable man who is not deterred from using a hand' truck for moving his trunk by the consideration that he has thus to move a trunk plus a truck, so am I reasonable in putting into your hands this useful, indeed indispensable, tool for mastering the ABC of iron metallurgy. Somewhat later, in introducing a discussion of the phase rule, which "it would be hardly proper to pass ... by without an attempt to outline its meaning," for "a clear exposition of its applicability or jurisdiction is needed to restrain the half¬initiated from misusing it and the uninitiated from being misled by these misuses," he writes-and I quote his statement because we cannot properly discuss the influence of time on a process until we know the state of equilibrium which would ultimately be reached by the system under consideration: It is a most remarkable and valuable generalization, its conceptions help greatly toward getting a broad outlook on metallography, but its misconception has brought out a flood of obscuring writings. It tells us about the constitution towards which alloys tend, that which they reach when equilibrium is complete, when all tendencies have corrected themselves and have been complied with completely. Its application, in absence of true equilibrium, is like determining density with a telescope or even with a tuning fork; it is like the mad tea party at which the butter put into the works of the watch was the best butter. If the reader learns no more than to beware of attempts to decipher the conditions of inequilibrium by a law which touches only the conditions of exact

Jan 1, 1942

By G. A. Mancini, V. Bharucha, J. W. Spretnak, G. W. Powell

The characteristics of the motion of the a interface during the "down" transformation was studied in zone-refined iron and dilute binary alloys containing nickel and molybdenum by means of the thermionic emission microscope and high speed photography. The frontal movements were found in all cases to be definitely discontinuous with time and an appreciable number of incremental jumps were found to proceed in the reverse direction (a to b). Molybdenum was found to be particularly effective in increasing the frequency of the increments of transformation from the stable to the metastable phase. A knowledge of the mechanism by which alloying elements affect the hardenability of a steel is desirable from both the practical and scientific viewpoints. It is well established that microstructures consisting of tempered martensite and lower bainite exhibit the best performance in terms of toughness and resistance to crack propagation. The presence of some upper transformation products, namely ferrite and pearlite, is in general undesirable. Various researches have indicated that the ferrite matrix grain structure is a very important difference between the upper and lower transformation products, namely equi-axed vs acicular. Any comprehensive theory of the decomposition of austenite should explain the dependency of the morphology of the matrix on undercooling. The contribution of various alloying elements to hardenability has been determined experimentally and catalogued in considerable detail. A rationalization of the basic mechanism or mechanisms by which individual elements affect the kinetics of the decomposition of austenite to products other than the phase "martensite" has not been evolved as yet. Interest in this problem was stimulated by the effects of minute quantities of boron in hardenable steels. This paper deals with one aspect of this problem, namely the characteristics of the ?-a transformation in high-purity iron and in two carbon-free iron binary alloys. A detailed historical account of the scientific developments in the study of the decomposition of austenite would be a prodigious undertaking and will not be attempted here. There are several excellent reviews in the literature on the formation of pearlite, martensitic reactions, and kinetics of solid-state reactions. The thinking on reaction mechanisms has been influenced to a large degree by the "classical" Tammann nucleation and growth model, in which an embryo springs into being by a fluctuation, and grows to a critical size through atom by atom deposition. Such a process is expected, accordingly, to be diffusion controlled. Considerable research attempting to relate the contribution of alloying elements to hardenability on the basis of their effect on the diffusivity of carbon in iron has not been productive from a mechanistic viewpoint. The formation of pearlite and ferrite has been treated in terms of nucleation and growth processes in which presumably the lattice transformation also occurs by a diffusion process. The other prominently considered mechanism is the "martensitic" or cooperative shearing mechanism. It seems generally, but not entirely, agreed that the phase "martensite" forms by this mechanism, and that it is involved to some extent in

Jan 1, 1962

By R. R. Vandervoort

The creep behavior of Nb(Cb), Ta, Mo, and W was determined under conditions of constant atomic dif-fzisivity, constant stress to elastic modulus ratio, and nearly equivalent grain size, and the steady-state creep rates obtained from these tests were correlated with calculated stacking fault energies for the metals. These results, in conjunction with similar data for several fccMetals,13 suggest that stacking fault energy may influence the creep strength ofbcc metals. The interrelationship between steady-state creep rate, subgrain size, and stacking fault energy was examined. It was found that the subgrain size for a given creep stress, increased as stacking fault energy increased, but that this relationship did not cormpletely account for the effect of stacking fault energy on creep rate. The crystallography and energetics of stacking fault formation in bcc metals has been discussed by a num-ber of authors,1-5 and impurity stabilized stacking faults on (112) planes have been observed in Nb,6,7 w,8,9 Fe,] and V" by transmission electron microscopy. However, a crucial question is whether or not stack-ing faults influence the mechanical strength of bcc metals. Potentially, stacking faults could increase strength by reducing the mobility of the partial dis-locations bounding the fault, by acting as barriers to slip dislocations, and by retarding the climb of dislo-cations during high-temperature deformation. The objective of this study was to seek a correlation be-tween creep strength and stacking fault energy for several bcc metals; namely, Nb, Ta, Mo, and W. The creep behavior of most polycrystalline metals and alloys at high temperatures and moderate stresses can be described by the following relation:11,12 im=Af(s) where i, = minimum creep rate, A = constant, j(s) = a function involving metallurgical structure, a = applied stress, E = average elastic modulus at the test tempera-ture, w = constant (equal to 5 for most pure metals), D = diffusion coefficient. One factor in the structure function F(s) which sig- R. R. VANDERVOORT, Member AlME is Research Metallurgist, Process and Materials Development Division, Chemistry Department, Lawrence Radiation Laboratory, University of California, Livermore, Calif. Manuscript submitted February 28, 1969. IMD nificantly affects the creep resistance of fcc metals is stacking fault energy, and creep rate has been shown to vary directly with stacking fault energy to the 3.5 power." In the latter investigation, four fcc metals of widely different stacking fault energies (Ag, Cu, Ni, and Al) were creep tested at a constant stress to modulus ratio of 1.21 x 10-4, at a constant diffusivity of 2.7 x 10-12 sq cm per sec, and at nearly equivalent grain sizes of about 0.7 mm. The creep data were then correlated with stacking fault energies. In the present study, a similar procedure was followed. All materials used in this work were consolidated by powder metallurgy techniques. Impurity contents in the as-received materials are listed in Table I. Chemical analyses showed that no measurable contamination of the test specimens occurred during pretest annealing treatments or creep testing. Specimens with a gage section 0.75 by 0.125 by 0.050 in. were creep tested in tension in a vacuum of less than 10-9 torr. Deformation at temperature was measured by tracking fiducial marks on the gage section of the specimen with an optical comparator. Optical deformation measurements also permitted observation of the macroscopic characteristics of the deformation Table I. Typical Specimen Impurity Content, ppm Nb Ta Mo W C 45 10 155 6 O 185 30 4 10 N 30 6 3 2 H 5 I 1 <1 als 3 10 2 15 Ca <5 I3 5 Cr 5 <3 10 <5 Cu 10 50 2 15 Fc 10 10 150 35 Ni 2 150 20 <5 Si <I0 1 3 <10 Ta 100 Ti 10 8 1 Zi 15 50 1 3 Table II. Test Conditions for Constant Stress-Modulus Ratio of 6 X 10.&apos; and Constant Diffusivity of 2.7 X 10-12 sq cm per see, and Grain Size Values for the Given Pretest Annealing Treatments Literature references Pretest Annealing for E and D Treatment Stress, Temperature, ___"&apos;Values__ Grain Tempera-Metal psi "C E D Size, mm ture, .C Time hr Nb 745 1525 14 15 to 17 0.83 1650 I Ta 1220 1770 18 19.20 O.91 1800 I Mo 1975 1630 18 21 0.77 2200 I W 2140 2265 18 22 040 2400 5

Jan 1, 1970

By Craig S. Hartley

Expressions are obtained for the energy changes associated with the reaction of (a& (111) slip dislocations on intersecting (110)planes in anisotropic bcc metals. An energy criterion for assessing the likelihood of dissociation of the products of such reactions is also presented. It is found that the "burrier reactions" which form a(100) dislocations at the intersection of two active {110) slip planes are more energetically favorable in metals which exhibit a high value of Zener&apos;s anisotropy factor, A, than those which have a low value. The results are presented in a form which permits the stacking fault energy to be obtained from a measurement of the separation between par-tials in a dissociated configuration. However, until accurate calculations or measurements of the stacking fault energies involved are available, it is not possible to assess the physical importance of dissociated dislocations. In a recent paper,&apos; the energy changes associated with several types of reactions between two slip dislocations, (a/2)(111){110), in bcc structures were calculated.* Isotropic elasticity and the approxima- tion v = -3- were employed. The purpose of this work is to present calculations of the energy changes for many of the same reactions using anisotropic elasticity. The problem of dissociation of a(100) and a(110) dislocations is also considered, and maximum fault energies for which dissociation will be energetically favorable are calculated for several bcc metals. Two general types of reactions are considered; those for which the reactant (a/2)(111) dislocations have long-range attractive forces and those for which the reverse is true. An example of the former is: (a/2)[lll] + (a/2)[lll]-a[l00] while the latter are typified by: (a/2)[lll] + (a/2)[111] -a[011] Only reactants lying in different slip planes are considered; therefore, the products must lie along (111) or (100) directions, which are the intersection of two {llO} planes. It will be assumed that the reactants and products are infinitely long parallel dislocations, since in this case the energy change associated with the reactions is a maximum.&apos; THEORY The self-energy per unit length of a straight mixed dislocation in an anisotropic medium can be written? where b is the Burgers vector, K is an appropriate combination of the single-crystal elastic constants, and R and ro are, respectively, outer and inner cut-off radii of the elastic solution. The energy given by Eq. [I] does not account for any variation of the core energy with orientation. This could be manifested by an orientation dependence of the core radius or, equivalently, the Peierls width, of the dislocation. However, the energy contribution due to this source is expected to be small, and current models of the dislocation core are not sufficiently accurate to justify such a refinement. It has already been shown that for the isotropic case the energy contributions due to nonzero tractions across the cores of the reactants and products exactly cancel one another in the reaction.&apos; Accordingly, it will be assumed that this contribution to the total energy change in the anisotropic case is small. In the subsequent discussion it is also assumed that the core radii of the reactant and product dislocation are the same and that, where stacking faults are formed, the faulted region is bounded by the centers of the partials. Consequently only changes in elastic energy due to the reactions will be considered. When the dislocation is parallel to either the (111) or the (100) directions, K may be written:375 K = (Ke sin2 a + Ks cos2 a) [2] where K, and Ks are the combination of elastic constants corresponding to an edge and screw dislocation lying along the same direction as the mixed dislocation, and a is the angle between the direction tangent to the dislocation line and the Burgers vector. Eq. [2] should not be confused with the isotropic approximation to the variation in energy with line Orientation.6 It should be noted that the essentially isotropic expression for K is a result of the characteristic symmetry of the (111) and (100) directions and is not, in general, valid for other dislocation directions in anisotropic cubic metals. The energy* change for a reaction in which the re- actant and product dislocations are parallel perfect dislocations can be written: where Ep and E, refer to the self-energies of the products and reactants, respectively. For dislocations parallel to (100) and (111) directions, Eq. [3] becomes:

Jan 1, 1969

By James K. Hallenburg

INTRODUCTION The delayed fission neutron, or DFN technique for uranium ore body analysis uses the first down-hole method for detecting uranium in place quantitatively. This technique detects the presence of and measures the amount of uranium in the formation. DFN TECHNIQUE DESCRIPTION The DFN technique depends upon inducing a fission reaction in the formation uranium with neutrons, resulting in an anomalous and quantitative return of neutrons from the uranium. Since there are no free, natural neutrons in formation, a good, low noise assessment may be made. There are several methods available for determining uranium quantity in situ. The method used by Century uses an electrical source of neutrons. This is a linear accelerator which bombards a tritium target with high velocity deuterium ions. The resulting reaction emits high energy neutrons which diffuse into the surrounding formation. They lose most of their energy until they come to thermal equilibrium with the formation. Upon encountering a fissile material, such as uranium, these thermal neutrons will react with the material. These reactions produce additional neutrons, the number of which is a function of the number of original neutrons and the amount of fissile material exposed. The particular source used, the linear accelerator, has several distinct advantages over other types of sources: 1. It can be turned off. Thus, it does not constitute a radioactive hazard when it is not in use. 2. It can be gated on in short bursts (6 to 8 microseconds). This results in measurements free of a high background of primary neutrons. 3. The output can be controlled. Thus, the neutron output can be made the same in a number of tools, easily and automatically. There are several interesting reactions which take place during the lifetime of the neutrons around the source. During the slowing down or moderating process the neutron can react with several elements. One of these is oxygen 17. This results in a background level of neutrons in any of the measurements which must be accounted for in any interpretation technique. These elements are usually uninteresting economically. The high energy neutrons will also react with uranium 238. However, the proportions of uranium 235 and 238 are nearly constant. Therefore, this reaction aids detection of uranium mineral and need not be seperated out. Upon reaching thermal energy the neutrons will react with any fissile material, uranium 235, uranium 234, and thorium 232. At present, we do not have good techniques for seperating out the reaction products of uranium 234 and thorium 232. However, uranium 234 is a small (.0055%) percentage of the uranium mineral and thorium 232 is usually not present in sedimentary deposits. When the uranium 235 reacts with thermal neutrons it breaks into two or more fragments and some neutrons. This occurs within a few microseconds after the primary neutrons have moderated and is the prompt reaction. One system uses this; the PFN or prompt fission neutron technique. We don&apos;t use this method because the neutron population is low and, therefore, the signal is small and difficult to work with, accurately. Within a few microseconds to several seconds the fission fragments also decay with the emmission of additional neutrons. Now, with a long time period available and a large neutron population we gate off the generator and measure the delayed fission neutrons after a waiting period. These neutrons can be a measure of the amount of uranium present around the probe. Thermal neutrons are detected with the DFN technique instead of capture gamma rays to avoid some of the returns from other elements than uranium. LOGGING TECHNIQUE The exact logging technique will depend, to some extent, upon the purpose of the measurement. However, the general technique is to first run the standard logs. These will include: 1. The gamma ray log for initial evaluation of the mineral body and for determining the position of the borehole within the mineral body, 2. The resistance or resistivity log for determining the formation quality, lithology, and porosity. 3. The S. P. curve for estimating the redox state and shale content, and measuring formation water salinity, 4. The hole deviation for locating the position, depth, and thickness of the mineral (and other formations), and 5. The neutron porosity curve. The neutron porosity curve is most important to the interpretation of the DFN readings. The neutrons from this tool are affected in the same way by bore hole and formation fluids as the DFN neutrons are. Therefore, we can use this curve to determine effect of the oxygen 17 in the water. Of course, this curve can be used to determine formation porosity. It can also be used to calculate formation density.

Jan 1, 1979

By Michael Epstein, Daniel E. Rosner

Fume nucleation sufficiently close to vaporizing suvfaces can augment net vaporization rates into cooler environments. Environmental conditions favoring large vaporization rate enhancements are briefly discussed and a previous theoretical treatment of this nucleation phenomenon is generalized to account for the self-regulating effect of condensalion-heat release within the boundary layer. Despite kinetic limitations on homogeneous nucleation, and latent heat release, non-diffusive condensate removal processes appear to make possible large enhancements in steady-state vaporization rates, provided surface temperatures are well below the boiling point. When condensed phases vaporize (or dissolve) into cooler media, the diffusion-limited mass loss rate can be strongly influenced by the process of nucleation/con-densation (or precipitation) within the thermal boundary layer. This condensation process (which typically leads to mists or fumes in the case of evaporation into cooler gases) has the effect of steepening the vapor pressure profiles near the evaporating surface, since the condensation zone acts as a vapor sink. However. the resulting enhancement in the diffusion-limited evaporation rate can be estimated (as first done by Turkdogan1 for the case of molten iron/nickel alloys evaporating into helium) only if one has independent knowledge of the critical supersaturation, sCrit(T), required to homogeneously nucleate the vapor.* In a recent reformulation and generalization of the theoretical model of Ref. 1 it has been shown that, when log sCrit is approximately linear in reciprocal temperature, rather simple expressions can be derived4 for the ratio of the actual rate of vaporization j" to either the minimum (no condensation) rate j"min, or the maximum (equilibrium condensation) rate j"max In the present communication we wish to briefly report on further developments and implications of the formulation of Ref. 4, with emphasis on i) environmental conditions favoring large enhancements in vaporization rate, and ii) the self-regulatory influence of condensation heat release (neglected in Refs. 1 to 4) on predicted vaporization rates. Additionally, we take this opportunity to correct several misprints appear- ing in Ref. 4, and comment on Elenbaas&apos;s recent criticism5 of Ref. 1. MAXIMUM POSSIBLE VAPORIZATION RATE IN PRESENCE OF CONDENSATION A nonequilibrium theory is of interest because of the very large difference between the minimum (no condensation) and maximum (equilibrium condensation) vaporization rate. The magnitude of this maximum possible enhancement can be shown quite clearly by combining a result of Refs. 3 and 4 with the fact that for most liquids there is a simple relationship between the molar heat of evaporation, A, and its boiling point, i.e., A/(RTBp) = C, where the constant C, often called the Trouton ratio, takes on values not very different from 11.* More generally, for any substance (including The Trouton ratio (which for water is 13, for methane, 10, and so forth) will be recognized as the ratio of the molar entropy change upon vaporization (at TBP or Ttransf) to the unlversal gas constant R. Its near constancy reflects the fact that the change in atomic order upon vaporization depends only weakly on the kinds of molecules involved. those that sublime under ordinary conditions) we can define a characteristic transformation temperature. Ttransf, by a relation of the form Ttransf =A/(CR), and then examine the maximum possible evaporation rates as a function of how far removed from Ttransf are the surface temperature, Tw, and ambient temperature, T. Subject to the assumptions: 1) equilibrium vapor pressure, pv,eq, everywhere small compared to prevailing total pressure, p, and 2) negligible effect of condensation heat on temperature profile, the maximum enhancement ratio was found (Eq. [17], Ref. 4) to be: where, for most vapors, Nu/NuD (the ratio of heat transfer coefficient to mass transfer coefficient for the same configuration) is a number near unity.* Ex- *An alternative derivation of the Nu = NuD special case of this equation. revealing its validity for arbitrary velocity/temperature profiles in a laminar boundary layer, is given in Ref. 3. amining this result for a "Trouton substance", one obtains the results shown in Fig. 1, constructed for C = 11. Since we are concerned with vaporization enhancements (j&apos;max/J"min > 1) at surface temperatures below Ttransf, this area of interest is shown unshaded. One notes that at a fixed ambient temperature (hence, T/TtranSf) there is a unique surface temperature, 2T , at which j"max/j"min attains its peak value; moreover, the peak enhancement ratio, see dashed locus. Fig. 1, is: (NuA/NuD)(C/4)(Ttransf/T,). Hence, if Nu = NuD, when the ambient temperature is less than 1/4 of TtranSf the peak enhancement exceeds the Trouton

Jan 1, 1969

By F. Selig, A. S. Odeh

A second-order approximation to the exact solution of the diffusivity equation corresponding to the pressure build-up of a well producing at a variable rate is derived. This approximation is applicable when the well&apos;s shut-in time is larger than the total time elapsed since the well was first produced. The resulting equations are compact in form and easy to use. Thus, the need for Horner&apos;s&apos; theoretically precise but rather laborious solution to the above problem is eliminated. In addition, these equations apply where the use of Horner&apos;s widely known approximate method is questionable. From a practical point of view, the reported method is best suited for analysis of drill-stem tests and short production tests conducted on new wells. INTRODUCTION The utility of drill-stem and short production tests in reservoir studies has long been recognized by the reservoir engineer. If interpreted correctly they could lead to a wealth of information upon which may depend the success or failure of reservoirs&apos; analyses. Initial reservoir pressure and the average flow capacity are two quantities that are normally sought from a drill-stem and/or a short production test analysis. Pressures are the most valuable and useful data in reservoir engineering. Directly or indirectly, they enter into all phases of reservoir engineering calculations. Therefore, their accurate determination is of utmost importance. The flow capacity kh of the reservoir is indicative of its commercial capability. In addition, it can indicate the presence of a darnaged zone around the wellbore and, thus, the necessity for remedial measures. Of the several methods used to analyze drill-stem and short production tests, Horner&apos;s&apos; method is by and large the most common. It applies to an infinite reservoir and or a limited reservoir where the effect of production has not been felt by the boundary. Horner&apos;s method makes use of the so-called "point-source" solution of the diffusivity equation. The point-source solution is approximated by a logarithmic function and the superposition theorem is utilized to give the familiar pressure build-up equation where is the shut-in time, q is in reservoir barrels per day and the rest of the symbols conform with AIME nomenclature. Eq. 1 was derived for a well which produced at a constant rate q from time zero to time t and was then shut in. In actuality, such a constant rate of production does not normally obtain. Therefore, a correction must be applied to Eq. 1 to account for the varying rates of production. Horner suggested two methods. The first, which results in a theoretically accurate solution, is rather lengthy and laborious and, thus, it is not suited for routine analysis. The second which has been termed a "good working approximation" is the one used by the majority of the reservoir engineers. In the second method, Eq. 1 is modified by simply introducing a corrected time t, and writing where q is the last established production rate prior to shut-in, and t, is obtained by dividing the total cumulative production by the last established rate. Horner&apos;s original paper does not give any indication that this method of correction is based on any theoretical justification. In addition, there is a question as to what constitutes the last established rate. In case of a drill-stem test some engineers use the average rate obtained by dividing the total fluid produced by the total flow time, while others calculate the average rate by dividing the total fluid produced by the last flow-period time. Obviously, different results obtain for the different flow rates used. Because of this, a simple method to the varying-rate case was developed which is theoretically sound and which defines clearly the flow rate and its associated time to be used in the calculations. The final equation arrived at is where q* and t* are a modified rate and time, respectively, and can be easily calculated. In addition, it is shown theoretically that Horner&apos;s approximate method, if used for a variable-rate case, gives the correct pressure but would not be expected to give the correct flow capacity. MATHEMATICAL ANALYSIS The general equation governing the flow of slightly compressible fluid in porous media may be written as The elementary solution to Eq. 4, representing an instantaneous withdrawal of Q units volume of fluid at the origin at t = 0, is known as the instantaneous sink

By R. Wild, F. A. Wright

When a blast furnace has a certain defined burden and is operated under fixed conditions of blast temperature, etc., the fuel efficiency is determined by the extent to which the reducing gases can remove oxygen from the burden in the furnace stack. This is determined by two distinct factors: 1) The uniformity of gas-solid contact, and 2) The ease with which oxygen can be removed from individual pieces of burden. This latter is often called burden reducibility. When burdens were poorly prepared the first factor was by far the most important and a study of the reducibility of individual lumps was of rather academic interest. In recent years good burden preparation with emphasis on uniformly sized material has led to greatly improved gas distribution in the stack, and thus the second factor has become much more important and there has been a marked increase in interest in methods of measuring reducibility. This paper explores the Linder method of measuring such reducibility. The measurement of reducibility of burden materials must be carried out under conditions duplicating, as nearly as possible, those of the blast furnace stack. This is very difficult since the blast furnace process is a counter-current one, and thus the initial conditions encountered by the solid (gas temperature, composition, etc.) are the result of heat and mass transfer occurring lower down the stack. Any method of burden testing which does not take this into account is, at least to some extent, based on arbitrary assumptions. In an attempt to study blast furnace reactions under non-arbitrary conditions BISRA adopted the SCICE technique as a method of investigation. This technique has been used with a measure of success.&apos; The SCICE technique, however, was found to be too slow for use as a routine test for burden materials and it was decided to construct additional equipment for burden testing. A test was required which would: 1) Be as realistic as possible. 2) Be quick and easy to operate. 3) Give some indication of the breakdown likely to take place during reduction in addition to a reducibility index. After a critical assessment of the reducibility tests which have been proposed it was decided to adopt the Linder test equipment and procedure as a basis for burden testing. THE LINDER TEST APPARATUS AND PROCEDURE The apparatus which was constructed (Fig. 1) was the same as that described by Linder2 except for minor changes in design. Linder also laid down a test procedure which he had derived from the results of investigations on Swedish blast furnaces. The variations of temperature and gas composition during the test were defined; these are shown diagram-matically in Fig. 2. BISRA&apos;s intention was to use the standard temperature and gas composition programmes for testing a variety of burden materials and also to investigate the influence of different programmes on standard burden materials, making use of information from the SCICE apparatus wherever this is possible. Up to the present, effort has been concentrated on the first part of the programme, and work on the second part has only just commenced. For each test 200 g of coke and 500 g of burden material were used. Linder had recommended that the coke and burden material should be between 1 and 1 1/2 in. and this was adhered to in early experiments on ores and sinters. Since the eventual aim of this work was to relate the test results to blast furnace operation, it was decided to carry out subsequent experiments using burden material in the size range used in the blast furnace, as far as this was possible. If the main interest was, for example, a comparison of the products resulting from different methods of agglomeration, then there would be advantages in using burden materials as close as possible to a standard size. After the charge had been placed in the reaction tube and this had been connected to the gas supply, rotation of the reaction tube at 30 rpm was started and the reduction programme was commenced, the gas temperature and composition being manually controlled according to the programme shown in Fig. 2. After the reduction test the charge was cooled in a nitrogen atmosphere. It was then removed from the reaction tube, the coke and burden separated, and the extent of burden breakdown assessed by screening it at 10 and 30 mesh. The extent of reduction was then

Jan 1, 1964

By R. G. Treuting, W. C. Ellis

The twinned state is characterized by a lattice of coincidence sites. Imperfections are required at stable lateral twin interfaces. Twinned regions can occur with relative ease in the diamond cubic IN recent contributions1,2 on the origin and growth of cubic annealing twins, attention has been directed to the orientation relations between such twinned components and their parent matrix. There are some aspects of twinning which may be illuminated by a more detailed consideration of the twinned state" alone. As an extreme example, the dense twinning in cast ingots of germanium,&apos; as contrasted with the rarity of twins in cast face-centered cubic metals, is yet to be accounted for. It has been this that has led us to the present work, which, it will be noted, uses methods and constructions in many respects similar to those of Kronberg and Wilson.&apos; In the cubic systems, a 70" 32&apos; rotation about a <110> axis is angularly equivalent, as to twinning, to the more usually considered 180" rotation about a <111> axis. Figs. 1 and 2 show a (110) projection of a twinned face-centered cubic lattice and a twinned diamond cubic lattice. In both figures, the two adjacent planes A and B, shown by the larger and smaller circles, are sufficient to represent the entire array. In each case a section of lattice, the original atom sites of which are shown by open circles, has been rotated as indicated through 70" 32&apos; to bring an original. [112] direction into coincidence with the [112] diiection. The latter is the intercept on the (110) projection of the (ill) plane normal thereto, the twinning plane. In the face-centered cubic case the rotation can be performed about an axis passing through an atom-site; the mirror plane then is also a composition plane containing atoms common to both twinned and untwinned lattices. The diamond cubic lattice may be construed as two interpenetrated face-centered lattices. Its (111) planes recur in a sequence of alternately short and long interspacings. Consequently a mirror plane for twinning cannot be a composition plane, but must be the bisector of one of the spacings. When the longer spacing is selected, the closest distance of approach across the mirror plane in the [ill] direc- tion is identical with that in the untwinned structure. In each case periodically recurring (ill) planes (parallel with the twinning plane) are found, on which there is coincidence of atom sites of the pre-twinned and twinned orientations; these are indicated by the cross-hatched circles. In the face-centered lattice there is such coincidence every third (ill) plane; in the diamond cubic lattice, on two adjacent planes in every six. At the twinning interface in the latter, there is on each side of the mirror plane a (ill) plane of atoms common to both twin components. Conceivably, there is little influence on a plane of atoms about to be adhered to such a pair of coincidence planes, whether it be laid down in a normal or in a twinned position with respect to the previously formed structure. Slawson% as attributed the high incidence of twinning in diamond to this boundary state. Further examination shows that the motion of intermediate planes can consist of various pairs of equal and opposite translations, for example of (ill) planes in the [l&apos;;i2) direction, the familiar twinning shear, indicated in the small schematics in the figures. Since the translations form a system of shears of alternating sign between coincidence planes, twinning could take place by such a mechanism over an extended region without extensive shear; in fact, in this case any atom moves but the distance in the [1i2] direction. One alternative construction for the face-centered cubic lattice leading to the same end result is illustrated in Fig. 3. The plane (711) with respect to the pretwinning orientation (the twinning plane of Fig. 1) is given, the twinned region arbitrarily bounded by <110> and <112> directions. The coupled shear is identical to that of Fig. 1. The "rotational" movement about coincidence sites generating the same twinned position could consist as shown of the translation a,/d% for each atom of a group of three in the B layer in a different one of the three <112> directions, and a similar translation of the underlying three atoms in the C layer in either the same or the opposite sense. This is not dissimilar to Kronberg and Wilson&apos;s construction for their 22" rotation of three adjacent (111) planes.

Jan 1, 1952

By M. D. Arnold, H. J. Gonzalez, P. B. Crawford

A method is presented for estimating the effective directional permeability ratio and the direction of maximum and minimum permeabilities in anisotropic oil reservoirs. The method is based on the principle that production from a well in an anisotropic reservoir results in elliptical isopo-tentials about the well, rather than circular. Bottom-hole pressure data from three observation wells surrounding a producing well are required to apply the method. The method involves fitting field pressure data to a set of general charts of isopotentials and making a few simple calculations until a solution is found. The method is based on a steady-state equation for homogeneorrs fluid pow. In addition to the method, a brief discussion of the theory underlying it is presented. INTRODUCTION The existence of a different permeability in one direction than another in oil reservoirs has been mentioned in several papers. Hutchinson&apos; reported laboratory tests on 10 limestone cores and pointed out that one-half of them showed significant, preferential, directional permeability ratios, the average being about 16:1. Johnson and Hughesz reported a permeability trend in the Bradford field in the northeast-southwest direction with flow being 25 to 30 per cent greater in that direction. Barfield, Jordan and Moore -eported an effective permeability ratio of 144:1 in the Spraberry. Crawford and Landrum4 showed that sweep efficiencies could often vary by a factor of two to four, and sometimes considerably more, due to variations in flooding direction and patterns in anisotropic media. These findings indicate that the poss&apos;bility of anisotropy may be worthy of consideration in the development of an oil field. In considering this, it should first be determined if anisotropy exists. If it does, the direction of the maximum and minimum permeabilities and the ratio of their magnitudes are quantities which can be of value in planning the most efficient well-spacing patterns. Past methods of determining these quantities have included analysis of oriented cores and analysis of flooding performance of pilot injection patterns. In recent work, Elkins and Skov5 resented an analysis of the pressure behavior in the Spraberry which accounted for anisotropic permeability. This work was based on the transient pres- sure distribution in a porous and permeable medium, with the solution expressed as an exponential integral function involving rock and fluid properties. The purpose of this study is to provide a method, based on steady-state equations, of estimating the direction and relative magnitude of permeabilities in an oil reservoir from field pressure data and well locations only. The method presented is based on work by Muskat6 which shows that Laplace&apos;s equation represents the steady-state pressure distribution for homogeneous fluid flow in homogeneous, anisotropic media if the co-ordinates of the system are shrunk or expanded by replacing x with it is desirable that data be obtained early in the history of a field because knowledge of an anisotropic condition would allow new wells to be spaced in such a manner that reservoir development and subsequent secondary recovery programs could be planned more efficiently. THEORETICAL CONSIDERATIONS A brief discussion of the theoretical basis on which the graphical solution was developed is presented in this section. Muskat&apos;s two-dimensional6 olution for the pressure distribution in an homogeneous, anisotropic medium with an homogeneous fluid flowing can be algebraically manipulated to show that the isobaric lines are perfect ellipses. The ratio of the major axis to the minor axis, a/b, is related to the permeability ratio, k,/k,, as follows. alb = dk,/k,--...........(1) It can also be shown that the pressure varies linearly with the logarithm of the radial distance from the producing well. However, the gradient along any ray is a function of the orientation of that ray, and a ..xiable is present when anisotropy exists which cancels out for a radial (isotropic) system. For a system such as that described, a dimensionless pressure-drop ratio was developed which is completely independent of the actual magnitude of the pressures. This was done by arranging Muskat&apos;s solution in such a way that aIl variables cancelled out except k,/k, and well positions. However, this solution depends on having a co-ordinate system with axes coinciding with the major and minor axes of the elliptical isobars. Thus, it was necessary to introduce a co-ordinate system rotation factor. The two unknown variables are then k,/k. and 0, and the two measured dimensionless pressure-drop ratios are related to the unknown variables as follows.

By Keijiro Nishio

At a time of great antiquity when our Yomato tribe had not yet found its way throughout the country, there lived in Japan barbarous tribes of the stone age, whose dwellings were vertical caves covered with roof&apos;s of twigs and weeds. In a cave lately discovered at the village of Morita, 10 miles west from the port of Aomori, there were found pieces of micaceous iron-ore and psilomelane, with oval holes pierced from both sides. Apparently these objects were regarded as ornaments or curios, without being applied to any practical use. Among the traditions given in the Nihonshoki, an ancient history of Japan, there is one relating to the mining industry. According to this tradition, in the mythological age Izanami-no-Mikoto gave birth to Kanayamahiko-no-Mikoto, god of mines, and Amaterasu-Okami (the sun-goddess) was armed with swords. It came to pass at one time that she concealed herself in Ama-no-Iwato ("Heavenly Cave"), and several gods made contrivances to entice her out. Among them, Ishikori-dome, a smith, made a mirror with copper obtained from the Ama-no-Kagoyama, by means of a pair of bellows known as the Ama-no-habuki, which was made of a deer-skin. It is recorded also that Susanb-no-Mikoto, her younger brother, observed the abundance of gold and silver in Korea, where his descendants were to reign, and urged them to provide themselves with boats to convey them across the Sea of Japan. In the course of time he accomplished his object, proceeding to Korea with his son, Isotakeru-no-Mikoto. These accounts

Jan 1, 1913

By J. H. Henderson, T. E. Morrisson

An irvestigation has been made of the production of oil through horizontal fractures of high capacity and large radius placed at the base of producing formations. The specific aitns of the study were to (I) determine by model studies the reservoir performance when gravity drainage is the only producing mechanism, and (2) develop a method of predicting such perfortnance from commonly measured reservoir proper-ties. Although the theoretical approach tailed to develop a rigorous analytic sorution to the transient problem, it (lid lead to a practical method of predicting performance. Predictions (Ire in good agreement with rates and recoveries observed in scaled-model experinzents. Capillary effects and the influence of relative permzenbil-ity are taken into account. The experimental program was mainly concerned with the influenre of fracture radius, fracture capacily and capillary hold-up on the total recovery and the rate of recovery. Much of the experimental work was done on a triangular model which represented one-eighth of a square well pattern, so constritcted that horizontal fractures of varying radii and flow capacity could be studied. An electrolytic model was used to relate the initial production rates to fracture radius, formation thickness and the drainage radius of the well pattern. The conclusions reached as a resrtlt of the work are as follows. I. The production of oil under gravity force alone into a well-propped horizontal fracture, whose radius is of the order of hundreds of feet, is a recovery method which combines high recoveries with very acceptable production rates. 2. The mechanism is effective in "dead oil" reservoirs as well as in reservoirs which possess the energy of solution gas. 3. A method of predicting performance has been developed which gives good agreement with experimental behavior in cases of practicnl interest. Even for very thick beds and large fracture radii, where actual conditions do not conform to the assumptions made, the method gives order of magnitude predictions which would not otherwise be possible. INTRODUCTION Gravity drainage of oil from a formation which has continuous ver- tical permeability over a reasonable depth interval is a surprisingly effective method of recovery. However, for a conventional well the production rates which result from gravity drainage are so low that the economics are seldom attractive. The same formation may produce rapidly and profitably by gravity alone through a high capacity, extensive (i.e., > 100-ft radius) horizontal fracture placed near the bottom of the pay section. Fig. 1 shows a schematic diagram of the plan. Energy in the reservoir gas is not a necessary condition for this recovery process. With a sufficiently large fracture, attractive production rates will result from gravity flow alone, and the producing GOR need only be the solution GOR. It is, of course, necessary to keep the fluid level in

The Sub-committee on Coal Mining, since its appointment, has been collecting data about subsidences in coal mines in this country, and thinks that almost all of the available data about occurrences of this sort in bituminous mines has been collected. Similar data is being compiled for the anthracite regions, which it is hoped will be ready for the 1927 meeting. As most of the bituminous seams are nearly level, the records obtained give the problem in its simplest form, without complications due to the pitch of the seam. Copies of the questionnaire, Appendix A, were mailed to all of the coal operators&apos; associations, to most of the large operating companies and to many engineers, with the request that any data on the subject be returned to the committee, or that advice of any subsidences known be given: the replies to most of the questionnaires were general, although some specific cases were given. Many of the companies were unable to reply to all questions, as they lacked the necessary data, but in nearly every case where accuratc data was available, it was cheerfully furnished, and any additional information requested was secured. All specific cases of subsidence were followed up by correspondence, and in several cases by personal inspection. Figs. 1 to 31 were made from the information furnished. Table 1 gives measurements of elevations taken in one case at intervals from October, 1913, to May, 1918. With one exception, where the evidence cannot be released on account of a lawsuit, it is believed that the cases submitted with this report include practically all that have happened in bituminous coal mining in this country of which actual measurements were secured. Some of these have been published by Young and Stoek,l and some have been secured from other publications. In collecting these data, the subcommittee had the Engineering Societies Library extend its bibliography to date; this list is given in Appendix B. To assist in the study of the problem, all of the subsidence literature in English was examined for cases in which spccific measurements were given, and the essential data for each case are shown in Appendix C.

Jan 1, 1927