Search Documents

By R. L. Parsons, Herman Dykstra

A numerical method for solving partial differential equations in steady state fluid flow is described. This method, known as the "relaxation method," has two advantages over analytical methods: (1) practically any problem can be solved, and (2) a solution can be obtained quickly. A disadvautage is that the solution is not general. The method is applied to core analysis and relative permeability measurement to calculate constriction effects and to calculate the true pressure drop measured by a center tap in a Hassler type relative permeability apparatus. Further applications are suggested. INTRODUCTION Many problems in fluid flow cannot be solved analytically because of the nature of the boundary conditions. For many problems, however. an exact answer is not necessary because boundary conditions are not exactly defined or the parameters describing the porous medium are not accurately known. The relaxation method can be used to obtain an approximate answer easily and quickly for the flow of incompressible fluids in porous media. The method can also be used for other types of problems, such as determining the stress in a shaft under load. or the temperature distribution during steady state heat flow. In this discussion only calculations concerned with the flow of fluids in porous media will be considered. The method was introduced by R. V. Southwell in 1935.' THEORY The treatment given here follows that given by Enimons.2 Consider a porous medium to be replaced entirely by a net of tubes of equal length and uniform cross-sectional area as shown in part in Fig. 1. Assume that the net of tubes behaves exactly like the porous medium which it replaces; that is, the net can be made fine enough to reproduce exactly the porous medium. Assume also that Darcy's Law can be used to calculate the flow from one point to another point through these tubes. The flow from point 1 to point 0 is KA . ------ P-P) .......(11 where a is the distance between points: K is the "permeability" of a tube; A is the cross-sectional area of a tube; is the viscosity of the liquid in the porous medium; and (P1 — P0) is the pressure difference between point 1 and point 0. In like manner the flow can be calculated from points 2, 3, and 4 to point 0. The net flow into point 0 is Qo = KA/µa (P1 + P2 + P33 + P4-4P0) . . (2) MB For an incompressible fluid the net flow into point 0 will be zero or, Q. = 0. This says that at point 0 fluid is neither being accumulated nor depleted. 'Therefore. P1 + P2 + P3 + P4 - 4P0 = 0 .... (3) . If. now. with specified boundary conditions. the pressure i.; known at a finite number of points in a given region, as at the points shown in Fig. 1, Equation (3) will be satisfied at every point. If, on the other hand, the pressure is not known, the pressure can be guessed at these points. Then. unless the guess is perfect. Equation (3) will not be satisfied at all of the points. When Equatiol~ (3,) is not satisfietl. let d = P1 + I?, + P, + P, - If' .,....(4) where 6 is an apparent error and is called the residual at point 0. Equation (4) shows how much the pressure guess is in error at point 0 with respect to the surrounding points. A positive residual means that the pressure is too low, and a negative residual means that the pressure is too high. To bring the residual, 6. to zero in order to satisfy Equation (3). it is necessary to make changes in the pressure guesses. Equation (4) shows that a +1 change in Po will change the residual at point 0 by -4. A +1 change in the pressure at any of the four surrounding points will change the residual at point 0 by +l. Thus it can be seen that a change at any point will affect the residual at that point and the four surrounding points. By changing the pressure from point to point, all of the residuals can eventually be brought nearly to zero and the problem will be solved. This procedure is the essence of relaxation methods and is used to relax the residuals so that Equation (3) is satisfied at every point. The procedure can be most easily explained in detail by solving a simple problem. as Southwell says, "To explain every detail of a practical technique is to risk an appearance of complexity and difficulty which may repel the reader. A

Jan 1, 1951

By V. M. Karve, K. K. Majundar, K. V. Viswanathan, J. Y. Somnay

The effect of calcium, magnesium, iron (both ferrous and ferric) and aluminum ions, which are commonly encountered in a typical beryl ore, was studied in the flotation of pure beryl, soda-feldspar and quartz. The vacuumatic flotation technique was employed. With sodium oleate as collector and in the absence of any activator, beryl floated in a pH range of 3 to 7.5, whereas feldspar and quartz did not float at any pH up to 11.5. The pH range of flotation increased in the presence of the ions studied. With calcium and magnesium ions beryl floated from 3 to 11.5 pH and beyond, soda-feldspar floated beyond pH 6 and quartz floated beyond pH 8. Ferrous ion activation was found to be similar to that of calcium and magnesium. Activation by ferric and aluminium ions was found to be complex and the lower and upper critical pH for all the three minerals was around 2 and 10 respectively. These studies indicated the possibility of separation of beryl from feldspar and quartz even in the presence of calcium, magnesium and ferrous ions between pH 4 and 6. Flotation tests on a mixed feed of pure minerals in a 10 g cell revealed that beryl can be selectively floated from feldspar and quartz if ferric ion is reduced to ferrous state or if it is complexed. Beryl occurs mostly in pegmatites, and hence is associated with feldspar, quartz and micas and small amounts of other minerals such as apatite and tourmaline. The separation of beryl from these minerals is difficult because all the silicates accompanying beryl have more or less the same physical properties. Specific gravities of beryl, feldspar and quartz are 2.70, 2.56 and 2.66 respectively. Electrostatic separation has been suggested but no work has been reported. ' The adsorption of sodium tri-decylate tagged with Cl4 on beryl, feldspar and quartz reveal similarity in surface properties. Much work has been reported on the flotation of beryl from ores, either directly or indirectly as a by-product, but little is known about the fundamental aspects of beryl flotation. Kennedy and O'Meara3 laid emphasis on prior cleaning of the mineral surfaces with HF. Mica is removed first by flotation of beryl with oleic acid, around neutral pH. Runke4 introduced calcium hypochlorite conditioning in a final separation stage for activating beryl in a mixed beryl-feldspar concentrate, and after washing to remove the hypochlorite, floated beryl with petroleum sulphonate. The Snedden and Gibbs5 procedure is somewhat similar to that of Kennedy and O'Meara. Emulsified oleic acid is used as collector. Recently Fuerstenau and Bhappu6 studied the flotation of beryl, feldspar and quartz with petroleum sulfonate in the presence of activators and stressed the importance of iron in the flotation of beryl. From the studies conducted in this laboratory, it was found that feldspar and quartz as such do not float with sodium oleate, but in practice selective flotation of beryl from feldspar and quartz in an ore is found to be impossible with sodium oleate as collector. A glance at the chemical analysis of typical beryl ore indicates the presence of several ions like Ca ++, Mg++, Al + + + and Fe+++ in abundance and Ti++++ and Mn++ in traces. Hence, in an attempt to explain the behaviour of feldspar in the beryl flotation, the effect of Ca++, Mg++, Al+++ and Fe+++, which are known as gangue mineral activators7'8 has been investigated. Materials and Methods: Lumps of beryl ore (hand picked) were boiled with 10% sodium hydroxide and washed with distilled water. They were further boiled many times with 10% hydrochloric acid till no positive test for iron was obtained with ammonium thio cyanate. This was followed by thorough flushing with double distilled water. The lumps were crushed in a porcelain mortar and pestle under water. The minus 65 + 100 mesh fraction was used for testing and was always stored under distilled water. Pure feldspar and quartz were similarly prepared and the minus 65 + 100 mesh fractions collected. Inorganic ions tried as activators were ca++, Mg++ , Fe++, Fe ++ and A1 +++ . Calcium nitrate, magnesium chloride, ferrous ammonium sulfate, ferric ammonium sulfate and aluminum nitrate of G.R.E. Merck grade were used. B.D.H. technical grade sodium oleate was used as a collector. The vacuumatic flotation technique developed by Schuhmann and Prakash was used for studying the effect of pH on flotability. 7 The indications given by this work were confirmed by using 10 g miniature cell.'

Jan 1, 1965

By V. F. Zackay, W. A. Goering

ALLOYS of iron and aluminum up to 35 wt pct aluminum are single-phase solid solutions, and are of potentially wide applicability.1-3 In spite of early and continued interest1-4 little progress has been made until recently in the preparation and evaluation of sound alloys containing more than 6 wt pct aluminum. Vacuum-melting techniques for the production of ductile Fe-A1 alloys have been described recently.1-7 A. procedure for air melting these alloys is presented here. Low-carbon iron is induction melted without a slag in a rammed magnesia crucible. At the beginning of melt-down, aluminum pig (99.95 pct Al), charged in a clay-graphite bottom-pouring crucible is placed in a pot furnace at 1800°F. The primary deoxidation of the molten iron after melt-down is effected by the addition of 0.1 pct aluminum and 0.5 pct manganese. (Hilty and Crafts" have reported a significant increase in the deoxidation efficiency of the aluminum and manganese combination over that of the aluminum alone.) A more drastic deoxidation designed to reduce the oxyen content to the lowest possible level is accomplished by plunging metallic calcium to the bottom of the melt. This is done by wiring small cubes of the metal to a steel rod. A circular shield larger than the diameter of the crucible opening is attached to the rod so that any spa'ttering of the molten metal will not endanger the operator. Since the temperature of the molten metal is above the boiling point of calcium, the bath is vigclrously purged by calcium vapor. It is believed that the calcium-vapor treatment permits a homogeneous distribution of calcium in the melt. Owing to the vigor of the reaction the temperature of the molten metal should be kept below 2900°F prior to the calcium addition. A total of 0.05 pct calcium is added in two stages in this manner. The second calcium deoxidation is made just before charging the molten aluminum into the iron, in order that an excess of calcium be present for the remainder of the melt. The aluminum, which has been removed from the holding furnace, is then hydrogen degassed by bubbling chlorine through a quartz tube immersed in the molten aluminum. The hydrogen-chlorine reaction is an exothermic one preventing the solidification of the aluminum during the 5-min chlorination. Approximately 0.1 pct calcium, based on the amount of aluminum, is then added to the aluminum. A further excess of calcium is introduced into the melt in this manner. The oxide dross is removed, fluorspar is added to the molten iron, and the molten aluminum is poured through the fluorspar slag. The fluorspar should be dried thoroughly prior to its use, as any water present will react with the aluminum. Aluminum oxide formed during the pouring operation reacts with the fluorspar slag to form gaseous aluminum fluoride and calcium oxide. A forced-draft ventilating system is required for this operation as aluminum fluoride is toxic. As soon as the molten aluminum has been added, vigorous manual stirring of the melt is required because the slag-aluminum oxide reaction is highly exothermic and tends to take place near the top of the melt. The combination of high temperature and the slagging action of the fluorspar quickly erodes the crucible at the slag line if the aluminum is not stirred uniformly into the melt. It has been found that at least 4 min of manual stirring combined with induction stirring are necessary to ensure homogeneity. The power is shut off 1 min prior to pouring to allow metal and slag to separate. As much slag as possible is removed from the melt, which is then poured directly into cast-iron molds. A mold wash of aluminum oxide is used to prevent ingot sticking. For slab ingots which are to be rolled into sheet, a carbon-tetrachloride vapor atmosphere or a chlorinated-pitch mold wash is desirable, as the aluminum oxide formed in the pouring operation is subsequently removed by the chlorine in the presence of carbon." As in vacuum melting, a pouring temperature of about 2900°F is recommended. Adequate hot-topping is important as iron-aluminum ingots are subject to very deep piping. Ingots are removed from the molds and buried in vermiculite, where they are allowed to cool slowly to room temperature. The ingots are radiographed,

Jan 1, 1959

By J. S. Bowles

THE martensite transformation in lithium, dis- covered by Barrett,' has been studied extensively by X-ray techniques by Barrett and Trautz,² and Barrett and Clifton.V he present paper reports the results of an investigation into the metallographic characteristics of lithium martensite. Such an investigation has not been carried out before. The spontaneous transformation in lithium consists of a change from a body-centered cubic to a close-packed hexagonal structure with the hexagonal layers in imperfect stacking sequence." As far as is known at present, this transformation can be regarded as being crystallographically equivalent to the body-centered cubic to close-packed hexagonal transformation that occurs in zirconium,5 although stacking errors have not been reported in zirconium. From a study of the orientation relationships in zirconium, Burgers5 as proposed that the martensite transformation, b.c.c. to c.p.h., occurs by a heterogeneous shear on the system (112) [111]. The crystal-lographic principle underlying this proposal is that the configuration of atoms in the (112) plane of a b.c.c. structure is exactly the same as that in the (1010) plane of a close-packed hexagonal structure based on the same atomic radius. The pattern in 2v2 both these planes is a rectangle d X 2v2d where v3 d is the atomic diameter. Thus a close-packed hexagonal structure can be built up from a body-centered cubic structure by displacing the (112) planes relative to each other.* This mechanism leads to orientations that can be described by the relations: (110)b.c.e. // (0001)c,p.h.; [111]b.c.c. // [1120]c.p.h Observations confirm these relations. In zirconium, Burgers' measurements indicated an angle of 0" to 2" between the close-packed directions, while Barrett's measurements on lithium indicated an angle of 3". According to the Burgers' mechanism, the martensite habit plane for this transformation would be expected to be the (112)b.c.c. plane, for this plane would not be distorted by the transformation. One of the purposes of this investigation was to find out whether the observed lithium habit plane agrees with this prediction of the Burgers' mechanism. Experimental Procedure Materials: The lithium was from the same purified ingot used by Barrett and Trautz.² The Bridgman technique was used to produce single crystals. To maintain a temperature gradient in the melt, during the production of these crystals, it was necessary to use a steel mould with a wall thickness of only 0.015 in. Metallographic Techniques: Lithium specimens could be given an excellent metallographic polish by swabbing them gently with cold methyl or ethyl alcohol.? The best results were obtained with methyl alcohol saturated with the reaction product, lithium alcoholate. With higher alcohols the reaction became progressively slower and the attack became an etch pit attack rather than a polish attack. Butyl and amyl alcohols were used for macroetching. After polishing, it was necessary to remove all traces of alcohol from the specimens; otherwise, on subsequent quenching in liquid nitrogen, the alcohol froze to a glassy film. The alcohol was removed with dry benzene. The benzene in turn had to be removed before quenching, but since it does not react with lithium it could be allowed to evaporate. The specimens could then be quickly quenched before they began to tarnish. This operation could be carried out in air on all but excessively humid days when it was advisable to use an atmosphere of dry nitrogen or argon. For examinations at room temperature, the specimens could be transferred directly from the benzene bath into a bath of mineral oil. In mineral oil the specimens oxidized slowly by the diffusion of oxygen through the oil but the structure remained visible for about an hour. Lithium Martensite: Specimens prepared in the manner described above transformed spontaneously to martensite with an audible click when quenched into liquid nitrogen; i.e., M, was above the boiling point of nitrogen (77°K). The disparity between this result and the M, temperature of 71°K, found by Barrett and Trautz, is probably to be attributed to the large grain size and freedom from mechanical deformation of the specimens used in the present work. The relief effects produced by the transformation did not disappear when specimens were quenched from liquid nitrogen into mineral oil at room temperature. This permitted the microstructures to be studied at room temperature where, of course, the martensitic phase was no longer present. Typical micrographs of lithium "martensite" made at room temperature are reproduced in figs. 1, 2, and 3. As anticipated by Barrett and Trautz, the microstruc-

Jan 1, 1952

By D. Havlena, A. S. Odeh

The material balance equation used by reservoir engineers is arranged algebraically, resulting in an equation of a straight line. The straight line method of analysis imposes an additional necessary condition that a successful solution of the material balance equatiott should meet. In addition, this algebraic arrangement attaches a dynamic ineuning to the otherwise static material balance equation. The straight line method requires the plotting of one variable group vs mother variable group. The sequence of the plotted points as well as the general shape of the resulting plot is of utmost importance. Therefore, one cannot progrm the method entirely on a digital computer ar is usually done in the routine solution of the material balance equation. If this method is applied, then plotting and anaIysis are asential. Only the appropriate equations and the method of analysis and interpremtion with comments and discussion are presented in this paper. Illustrative field examples for the various cases treated are deferred to a subsequent writing. INTRODUCTION One of the fundamental principles utilized in engineering work is the law of conservation of matter. The application of this principle to hydrocarbon reservoirs for the purpose of quantitative deductions and prediction is termed "the material balance method of reservoir analysis". While the construction of the material balance equation (MBE) and the computations that go with its application are not difficult tasks, the criteria that a successful solution of the MBE should fulfill have always been a problem facing the reservoir engineer. True and complete criteria should embody necessary and dcient conditions. The criteria which the reservoir engineer uses possess a few necessary but no sufficient conditions. Because of this, the answers obtained from the MBB are always open to question. However, the degree of their acceptability should increase with the increase in the number of the necessary conditions that they should satisfy. Generally, the necessary conditions commonly used are (1) an unspecified consistency of the results and (2) the agreement between the MBE results and those determined volumetrically. This second criterion is usually overemphasized. Actually, the volumetrically determined results are based on geological and petrophysical data of unknown accuracy. In addition, the oil-in-place obtained by the MBE is that oil which contributes to the pressure-production history,' while the volumetrically calculated oil-in-place refers to the total oil, part of which may not contribute to said history. Because of this difference, the disagreement between the two answers might be of paramount importance, and the concordance between them should not be overemphasized as the measure of correctness of either one. In this paper, a third necessary condition of mathematical as well as physical significance is discussed. It is not subject to any geological or petrophysical interpretation, and as such, it is probably the most important necessary condition. It consists essentially of rearranging the MBE to result in an equation of a straight line. This straight line method of the MBE solution has invalidated a few long time accepted concepts. For instance, it has always been advocated that if a water drive exists, but one neglects to take it into account in the MBE, the calculated oil-in-place should increase with time. The straight line method shows that in some cases, depending on the size of the neglected aquifer, the calculated oil-in-place might decrease with time. The straight line method requires the plotting of a variable group vs another variable group, with the variable group selection depending on the mechanism of production under which the reservoir is producing. The most important aspect of this method of solution is that it attaches a significance to the sequence of the plotted points, the direction in which they plot, and to the shape of the resulting plot. Thus, a dynamic meaning has been introduced into the picture in arriving at the final answer. Since the emphasis of this method is placed on the interpretation of the sequence of the points and the shape of the plot, one cannot completely automate the whole sequence to obtain "the best value" as normally done in the routine application of the MBE. If one uses the straight line method, then plotting and analysis are musts. The straight line method was first recognized by van Everdingen, et al,2 but for some reason it was never fully exploited. The advantages and the elegance of this method can be more appreciated after a few cases are carefully treated and worked out by it. SOLUTION OF THE MATERIAL BALANCE EQUATION SATURATED RESERVOIRS The MBE for saturated reservoirs written in AIME symbols is

By H. M. Cyr, T. F. Steele, C. W. Siller

The development of a new metallurgical roasting device is described. It consists of a refractory column into which air is injected at various levels, forming several superimposed fluidized beds with no supporting grates. When pelleted zinc sulphide concentrates are charged, the roasted product needs no further sintering before reduction to metal. WHEN a gas such as air is blown upward with increasing velocities through a loose mass of solid particles, marked changes in the physical behavior of the particles are noted. At first, when the velocity of the gas is insufficient to support any of the solid, the mass constitutes a "fixed bed." As the gas velocity increases until the pressure drop through the bed approaches the effective weight of the bed per unit area, the bed expands until the solid particles are supported by the air rather than by the lower particles. Some vibration of the particles becomes apparent, but little mixing occurs. This condition is called a "quiescent fluid bed." A further increase in gas velocity imparts more separation and more motion to the individual particles until a condition of turbulence is reached. This "turbulent fluid bed" resembles a rapidly boiling liquid with the characteristic highly agitated diffuse surface and many small eruptions of the boiling mass. Different degrees of turbulence can be generated and all produce excellent mixing. The final stage occurs when the gas velocity becomes so great as to create a "dispersed suspension." Here no surface of the mass is defined and the gas carries solid particles out of their original positions. These changing conditions of fluidization have been studied carefully and pertinent nomenclature standardized by a committee of the American Institute of Chemical Engineers.' Many mathematical analyses2-3 have been made of the forces acting in a fluid bed. These analyses are invaluable, especially for the design of column sizes and selection of equipment. However, in a metallurgical process involving solids of many sizes with changing densities, varying temperatures, and changing gas compositions within the bed, calculations based on theory become approximate. Optimum operating conditions then are best determined experimentally. Many applications have been made of the principles of fluid-bed action by mechanical, chemical, and metallurgical engineers. Especially when good con- tact between reacting solids and gases is desired, very effective results are obtained from fluid beds. They permit excellent temperature control and uniformity throughout a mass of solids in fluid action. Heat transfer to walls and any coolers is high, and fast reaction rates are attained because the solid surfaces are continuously swept clean. The main disadvantages of fluid-bed operations are the danger of short-circuiting in a single bed, danger of incipient sintering which stops action, the necessity of avoiding large changes in particle size or density during roasting, and dust losses when particles of the charge are carried out with exit gases. In the metallurgical field the roasting of sulphide ores to form oxides and sulphur dioxide appears to combine several operating conditions which can be achieved to advantage in a fluid bed. Roasting involves a solid-gas reaction where a high reaction rate is necessary for high capacity, where good temperature control is important in order to prevent sintering, where good heat transfer is needed, and where the density of the solids, when changing from sulphides to oxides, is not largely changed. Short-circuiting, however, constitutes a major problem when a single fluid bed is used. Because of the turbulence of the bed, an entering particle may be in the region of the discharge before it is roasted. Hence, to attain a satisfactorily low sulphur in the calcine, a long average residence time with correspondingly low capacity is required. The solution to this difficulty is the use of multiple stages, which in the conventional fluid-bed design requires separate hearths with feed and discharge mechanisms for each stage. A further practical difficulty in fluid-bed roasting of flotation zinc concentrates is their fine particle size which makes a true fluid action without excessive carry-over of dust very difficult to attain, especially when the large air volumes necessary for high capacity are used. A New Design After considerable experimentation in the laboratory and on a semipilot-plant scale, a new method and equipment for roasting were devised which provided a unique solution to these problems. A detailed account of this development appears in the patent literature," and many of the variations of this development reported herein are the subject of

Jan 1, 1955

By G. L. Chierici

Starting from Muskat's theory of water and gas coning, maximum permissible oil production rates without water and/or free-gas production have been determined, in a broad range of reservoir and well parameters, using the potentiometric model technique. The main assumptions made are as follows: (1) the reservoir rock is homogeneeous (either isotropic or anisotropic); (2) the volume of the aquifer underlying the oil zone is very small, so that it does not contribute to reservoir energy; and (3) the gas cap expands at a very low rate, so that it can be assumed to be in quasi-static conditions. The results obtained are presented in the form of diagrams which can be used for solving two types of problems: (1) given the reservoir and fluids characteristics, as well as the position and length of the perforated interval, determine the maximum oil production rate without water and/or free-gas production; and (2) given the reservoir and fluids characteristics only, determine the position and length of the perforated interval which optimize the maximum permissible oil production rate, without water and/or free-gas production. INTRODUCTION In oil reservoirs where the oil-bearing formation is underlain by an aquifer which does not participate in the production mechanism, water-coning is a limiting factor to the flow rates of producing wells. Production rates are usually kept to a value that will prevent the water from entering the wells. The entry of water into a well lowers its productivity by increasing the weigbt of the fluid column; moreover, the separation of water from the effluent, at the surface, may constitute a very difficult problem in cases of heavy viscous oils. A similar situation is encountered in oil reservoirs with a gas cap overlying the oil-saturated zone; here a downward gas cone is induced by the flow of oil towards the producing wells. Production rates must be low enough to prevent the gas from being produced; producing gas from the gas cap would be a waste of energy. Of course, water-coning and gas-coning phenomena can occur at the same time in the same reservoir if the oil-producing formation is both overlain by a gas zone and underlain by a water zone. Due to its relevant practical importance, the mechanism of coning was studied by many people.2,3,5-8 Defining the conditions for getting the maximum water-free and/or gas-free oil production rate is a difficult problem, often encountered under one of the following aspects: 1. Predict the maximum flow rate that can be assigned to a completed well without the simultaneous production of water and/or free-gas. 2. Define the optimum length and position of the interval to be perforated in a well, in order to obtain the maximum water and gas-free production rate. A systematic study of these problems was made by means of the electrical analog technique. The results of this study are presented here, under the form of a set of curves providing solutions for the above stated problems. These curves are valid only for homogeneous forrnations, either isotropic or anisotropic. Should the formation be non-homogeneous (by horizontal or vertical variation of permeability, shale diaphragms, fractures, etc.), a specific potentiometric study would be required for each specific case. Especially when shale diaphragms of some radial extension are present, the critical rates observed are much larger than would be expected from the diagrams. STATEMENT OF THE PROBLEM In the present study the aquifer is supposed to be of such limited volume that it does not contribute to the energy of the reservoir. Moreover, the gas cap is supposed to expand at such a low rate that the potential gradient in the gas cap is negligible. Under static conditions water-oil and gas-oil interfaces (T1 and T2) are both horizontal. When the reservoir production starts, below each well these interfaces take a cone-like shape (Fig. 1) having as an axis the axis of the well. This shape results from the equilibrium between potential gradients in the oil zone and gravitational forces due to density differences between oil and water and between oil and gas. Assuming the oil-bearing formation to be homogeneous and the oil to be incompressible, the analysis of the problem (see Appendix) shows that the oil-water and gas-oil interfaces are stable only if the oil production rate of the well is not higher than the following values.

Jan 1, 1965

By J. S. Rinehart, W. C. Maurer

The mechanics of impact crater formation in rock, particularly sandstone, has been sutdied, the velocity range being approximately that normally associated with oilwell gun perforators. The bullets were small steel spheres having diameters of 3/16, 9/32 and 7/16 in; impact velocities ranged from 300 to 7,000 ft/sec. The craters have two distinct parts — a cylindrical hole (or burrow) with a diameter the same as that of the impacting sphere, and a wide-angle cup comprising most of the volume of the crater. The burrow is fornred as material in front of the projectile is crushed and pushed aside, forming a cylindrical hole surrounded by a high-density zone. The clip forms as fractures are initiated in front of the projectile and propagate along logarithmic spirals, approximaling maximum shear trajectories, to the free surface of the rock. A most significant observation (made for the first time) was that, below the base of the cup in one type of sandstone, there are a group of similar fractures, not extending to the surface, which are spaced uniformly a few millimeters apart. Each fracture follows roughly the contour of the base of the cup and appears to require a certain threshold impulse to initiate it. These fractures comprise a relatively high fraction of the total, newly exposed surface area. The volume of the material removed by crushing varies as the first power of the impact velocity and the volume removed by fracturing, as the second power of the impact velocity. Penetration varies linearly with the impact velocity and is inversely proportional to the specific acoustic resistance of the target material, the proportionality constant being dependent upon the shape of the projectile. INTRODUCTION Yield of oil from a producing well is frequently enhanced by firing bullets and shaped charges through the well casing into the oil-bearing rock, forming craters and fractures from which oil can flow more readily. The purpose of this investigation has been to develop a better understanding of the mechanics of impact crater formation in rock, particularly sandstone, the velocity range being approximately that normally associated with oilwell gun perforators. FORCES OPERATIVE DURING IMPACT When a projectile moving at considerable velocity strikes a- massive target such as oil-bearing sandstone, intense and complex transient stress situations develop within both the projectile and the rock or sandstone against which it is striking. Usually the struck rock fails, the missile or projectile penetrating into the rock to some depth where it comes to rest or is forcibly ejected from its burrow by expansion of a plug of target material compressed in front of it. When the impact velocity is very high, the projectile itself may fail, breaking apart or becoming distorted; this situation is not considered here, the discussion being limited to nondeforming projectiles. Many experimental studies'.' have been carried out to determine the nature of the mechanics of crater formation and the salient features of the forces coming into play, some of the earliest studies being the French Army experiments performed at Metz between 1835 and 1845.' The stratagem in most instances has been to make a post-mortem examination of the crater, measuring volume and depth of penetration and deducing force relationships from these observations rather than performing the more difficult (usually almost impossible) feat of measuring stresses during penetration. In many materials, the force acting during penetration of the projectile is found to be the sum of two components—(1) a constant force, independent of the velocity, representing some inherent strength of the target material; and (2) a component, proportional to the square of the velocity, representing inertial forces. For such materials, the average force per unit area acting on the projectile at any instant while it is in motion and being decelerated may be written F/A = a + bv2 . . . (1) where v is the velocity of the projectile at that instant, A is the cross-sectional area of the penetrating projectile taken normal to its trajectory, and a and b are constants which are dependent upon the target material and the shape of the projectile. It follows that the total penetration s is given by .........(2) where v, is the velocity of the projectile when it just strikes the target. Values of a and b for spherical projectiles impacting in a loose sand-gravel mixture and compacted earth were obtained in the Metz experiments. For sand-gravel, a and b are 620 psi and 0.0115 (psi) (ft/sec)', respectively; and for compacted earthworks, a and b are 432 psi and 0.0008 (psi) (ft/sec)'. Figs 1 and 2

By William H. Breeding

The remarks that follow are based substantially on experience covering 45 years, 80% of which has been in placer work, rather than on a review of available literature. Most commercial placers have been deposited by the action of water. The richer and more- difficult-to-mine placers are those in the headwater areas where gradients are steepest. The most lucrative placers are generally in inter- mediate areas where volumes are greater, fewer boulders are present, and gradients are from 3% to 1-1/2%. The higher volume, lower grade placers are in the lower reaches of river systems where gradients are lower. Where gold-bearing rivers have discharged into the sea, wave action can concentrate values on beaches, past and present. Most of the rich, readily accessible placers were mined by our forefathers. Current opportunities exist: (1) in remote areas where infrastructure has been absent in the past, or development has been prohibited by adverse ownership - political or commercial; (2) in deposits that could not be mined by equipment available to our forefathers; (3) in deposits unidentified by our forefathers; (4) where the-price-of-product/cost ratio is substantially better than in earlier years; or (5) a combination of those factors. When I entered the placer business in the late 1930s, and subsequently, a prevailing opinion believed that glacial deposits should be avoided as irregular in mineral content and composition, and unrewarding to explore and develop; yet an operator has been mining a fluvio-glacial deposit profitably for the past 17 years. Rich buried placer channels, of ten called paleo-channels were worked in the last century, generally by hand methods, and under conditions that would be unacceptable today. Exploration and mining equipment now available make some of these channels attractive targets. Well-known examples are in California and Australia. The formation of a commercial placer requires a source of valuable minerals. Above primary deposits, there may be eluvial deposits formed by the erosion of gangue minerals and the concentration "in situ" of valuable minerals. Down slope from these deposits are the hillside or colluvial deposits, and below them are the alluvial deposits of redeposited material. Most of the great placer fields of the world are the result of several generations of erosion and deposition. Well-known examples are in California and Colombia. Gold is a very resistant and malleable material, and gold placers may extend for 64 or 80 km (40 or 50 miles) along a river system. Platinum is less malleable, but is very resistant to disintegration. Diamonds are extremely hard, and (especially gem diamonds) may be found over great lengths of a river system. Cassiterite is less resistant to disintegration, and tin placers seldom extend over two miles without resupply from an additional source or sources of mineralizaton. Tungsten minerals are generally more friable, and within a few hundred yards of the source disintegrate to the point that they are uneconomical to recover. Rutile, ilmenite and zircon placers generally result from the weathering of massive deposits, and may be encountered over extensive areas; most are fine grained and durable. What does a geologist or mining engineer look for in placer exploration? The old adage to look for a mine near an existing mine is still valid. You need a source of valuable mineral. Then you require conditions for concentration, which means a satisfactory gradient and/or other conditions that will permit heavy minerals to settle. Nicely riffled gravel, often called a shingling of the bars, is conducive to placer formation. Coarser gravel is logically associated with coarser gold. Excessive clay and/or high stream velocities in narrow channels can carry gold far downstream and distribute it uncommercially over a large area. When material is extremely fine, in situ weathering and concentration become more important. Placers frequently occur distant from lode mines, and one must remember that in a larger watershed the exceptional floods that occur once in a hundred or a thousand years can move great quantities of material long distances. The carrying power of water is said to vary with the fifth or sixth power of its velocity. I am not ready to disagree with Waldemar Lindgren and accept that many commercial placers are substantially enriched by the chemical deposition of gold from solutions; however, I have seen crystalline gold in clayey material quite distant from known sources of primary gold that is dif-

Jan 1, 1985

By G. H. Bishop

The kinetics of the inter granular penetration and embrittlement of [100] tilt boundaries in 99.998 pct pure nickel upon exposure to bismuth-rich Ni-Bi liquids have been determined in the temperature range from 700° to 900°C. The kinetics of penetration are parabolic in time at constant temperature over most of the temperature range. In a series of 43-deg bicrystals the rate of penetration is anisotropic with respect to the direction of penetration into the grain boundaries. In lower-angle bicrystals the penetration rate is isotropic. The rate of penetration decreases with tilt angle at 700°C. The activation energy for penetration in the 43-deg bicrystals is 42 kcal per g-atom independent of direction. It is concluded that the intergranular penetration and embrittlement in the presence of the liquid proceeds by a grain boundary diffusion process and not by the intrusion of a liquid film. This was confirmed by a determination that the kinetics of penetration and embrittlement were the same in the 43-deg bicrystals upon exposure to bismuth vapor under conditions such that no bulk liquid phase would be thermodynamically stable. WhEN solid metals are exposed to a corrosive liquid-metal environment, the grain boundaries are sites of preferential attack. Depending on the temperature, the composition of the liquid, and the composition, structure, and state of stress of the solid, a number of modes of attack are possible. This paper reports a study of the kinetics of intergranular penetration and embrittlement of high-purity nickel bicrystals upon exposure to bismuth which, together with an earlier study by Cheney, Hochgraf, and Spencer,&apos; demonstrates that there are at least two modes of intergranular attack possible in the Ni-Bi system. In the study by Cheney et al., columnar-grain specimens of 99.5 pct pure nickel were exposed to liquid bismuth presaturated with nickel in the temperature range 670" to 1050°C. They found that the majority of the boundaries, which were predominantely high-angle boundaries, were penetrated by capillary liquid films, the attack proceeding by a process which will be termed grain boundary wetting. This process occurs in a stress-free solid when twice the liquid-solid surface tension is less than the surface tension of the grain boundary,* i.e., when 2yLs < YGB In this case the penetration of the grain boundary by the liquid occurs at a relatively rapid rate, resulting in the severe embrittlement of a polycrystalline solid. Grain boundary wetting is a common mode of intergranular attack in systems in which the lower melting component is relatively insoluble in the solid, but the solid has an appreciable solubility in the liquid, for example, the Ni-Bi system, Fig. 1. In systems of this type at temperatures above the range of stability of any intermetallic phases, once the liquid is saturated with respect to the solid so that no gross solution occurs, chemical gradients are small, and surface tensions become major driving forces for attack, provided the solid is stress-free. The results of Cheney et al. appear to be typical of those encountered when grain boundary wetting occurs.&apos; Capillary films were observed in the boundaries after quenching from the exposure temperature. The mean depth of penetration increased linearly with time, and the activation energy for the process was found to be 22 kcal per g-atom. In a study of the Cu-Bi system Yukawa and sinott4 found that the depth of penetration of bismuth into high-purity copper bicrystals of orientations from 22 to 63 deg of tilt about (100) at 649°C ranged from 0.05 to 0.25 in. after a 12-hr anneal. This corresponds to a linear rate of 6 to 15 X 10~6 cm per sec. At the same reduced temperature of 0.68 the rate for the Ni-Bi system&apos; was 7 x lo-&apos; cm per sec. In another study of the Cu-Bi system, Scheil and schess15 determined the kinetics of grain boundary wetting in hot-worked commercial rod. While there were several complicating factors present in this study, there is general agreement with the above results. The kinetics of penetration were linear, the activation energy was 20 kcal per g-atom, and at 650°C the rate of wetting was 2 to 5 x 10-6 cm per sec. The rate of wetting in the A1-Ga system6 is somewhat

Jan 1, 1969

By R. O. Williams

Analytical representation of the thermodynamics of solutions is highly desirable from the standpoint of accuracy, compactness, and numerical manipulations. In particular, computer calculations are greatly implemented. Mathematical considerations show that previous expressions have one or more serious defects. This investigation shows a Fourier series to be satisfactory but that it is also possible to derive a new series which fits certain additional conditions. Included examples show the value of analytical expressions in giving a simple characterization of each system using some two to five parameters, the elimination of the Gibbs-Duhem integration, and the es timation of the error for the experimental function as well as derived functions. It is further shown that the present characterization provides easy comparison between systems. IN the past, thermodynamic calculations have depended to a considerable extent on tabular and graphical methods. As the volume and precision of such data increase such methods become less satisfactory. Specifically, the selection of the optimum representation and the estimation of errors require statistical methods which in turn require analytical representation. The utilization of such data require further manipulations which are best done analytically for maximum precision. For example, phase equilibria are determined by common tangents to free-energy curves: a graphical determination is normally of low accuracy. As computers are increasingly used analytical representations become almost mandatory. Insufficient mathematical consideration has been given previously to the selection of empirical expressions. Those expressions having some theoretical justification are generally too inflexible and mathematically unattractive. We consider the problem in some detail and show that a Fourier series can be effectively used. Also a new series is defined which has certain advantages. ANALYSIS We wish to consider the analytical representation of the heat of mixing, AH, the excess free energy, ?Gxs, and the excess entropy, ?sXS, as a function of composition, X, for binary solutions relative to the pure components in the same state. When a distinction is not required, we use W to denote any one of the above functions. One may use a Taylor expansion around X = 0 to generate a power series. As the derivatives are un- known we represent the series as W = A + BX + CX2 + DX3 + EX4 + ... [l] where the constants A , B, C , ..- are to be selected to provide some optimum fit. For the extremes of composition W is necessarily zero so it follows that A = 0 [2a] B +C + D + E +••• = 0 [2b] Nonelectrolytes, which we are considering, appear to satisfy the condition that d3W/dx3 = 0 [3] in the terminal regions. This is the basis of the a, ß, and Q functions used by Hultgren et al.&apos; and others. While this condition does not have a strong theoretical basis it does appear desirable that any analytical relation should satisfy this condition. Darken2 and Turk-dogan and Darken3 have shown that many systems exhibit this behavior over an extended range from each terminal region, departure being restricted to a limited intermediate region. Since we have no a priori knowledge as to where this transition occurs we can require that this condition be satisfied only as a limit at the extreme compositions as a general condition. We will show later how more restricted conditions can be included in specific solutions. Darken2 has called this behavior the quadratic formalism; we call our application the limiting quadratic formalism, LQF. This condition applied to the above power series requires that D = 0 [4a] 4-3-2E +5-4-3_F + 6 • 5 . 4G + ••• =0 [4b] The form of the power series normally used, due to Margules,4 is W=X(1-X)(A + BX + CX2 + DX3 + EX4 + •••) [5] where A, B, C, --. are a new set of constants. (Guggenheim5 has given a variation of this expression in a more desirable form. Since, however, it is contained in the above expression it does not require separate consideration.) This form is precisely what results by incorporating the conditions in Eq. [2] into the power series and regrouping the constants. The LQF requires that B =C [6a] and 4.3.2(D-C) +5-4-3(E-D) + ••• =0 [6b] Thus, the correct form of the Margules expression with two adjustable parameters is w =X(1-X)[A + B +X2-2/3x3)] 171 and the EX4 term must be included before three adjustable parameters are permitted.

Jan 1, 1970

By F. X. Tartaron

This paper deals with basic physical phenomena, which when combined and interpreted, lead to the same mathematical equations that describe comminution phenomena. Thus, a physical model is described that corresponds to the mathematical model presented in the writer&apos;s previous papers. &apos;12 In the mathematical model, the energy consumed in breakage is related to the volume or weight of material broken and the size of particles broken. The equation E=2.303 Ck-a log x1/x2 was derived by multiplying the volume or weight of each size in an ideal Gates-Gaudin-Schumann size distribution by an energy factor. The product of these two factors gives the energy distribution among the different sizes in a single size distribution. The energy of breakage of a specific constant weight of one size distribution to another size distribution is given by the equation E = constant/kn-1. In this case, where the volume or weight is constant, the energy is proportional to the size factor 1/kn-1. In what follows, a physical theory will be presented showing that the energy consumed in comminution is proportional to the volume or weight of the material broken and to the reciprocal of the size of this material raised to a constant exponent. THE VOLUME FACTOR The atomic theory of matter reveals that in solids, atoms or ions are arranged so as to be in equilibrium at specific distances from one another. Although the atoms or ions are oscillating, there is a definite determinable mean distance between them and this distance is a balance between repulsive and attractive electrical forces. It therefore requires force to separate the atoms or ions and when an outside force is applied, it first produces strain in increasing the distance between the atoms or ions. This strain increases to the breakage limit on application of sufficient force. In brittle materials, there is negligible plasticity and when an elastic limit is exceeded, breakage takes place. The work done is the force applied per unit area times the cross sectional area of the ideal particle multiplied by the maximum strain per unit length at right angles to the area times the length of the particle. Thus the work done is proportional to the area times the length, which is equivalent to the volume of the ideal particle. If more than one feed particle is considered broken, each particle must be subjected to sufficient strain so that the breakage limit of its contained atoms or ions is reached in order for the particle to be broken. Thus, the energy of breakage is proportional to the total volume of the particles broken. If the particles are of different sizes, the size factor must be included to get a correct determination of energy of breakage. In the preceding, it has been assumed that there is a constant binding force between the atoms throughout the volume being strained. This, of course, is not true. It is known that there are many irregularities in the structure of matter and the binding force differs markedly in different portions. But the differences are only discernible by examining extremely small subdivisions of matter. In one order of magnitude of volume, cracks can be discerned separately from non-cracked neighboring material. In a smaller subdivision of volume, lattice dislocations can be isolated. When these situations are brought into focus, mechanisms of their behavior can be learned, leading to a fuller understanding of phenomena that occur in larger scale subdivisions of matter. Very often, however, the mechanisms that operate in small scale subdivisions have negligible effect in those of large scale, and there still is a place for deriving a mechanism for large scale conditions. The quantum theory is extremely valuable for use with photons and electrons, but is of negligible use with ordinary atoms and molecules. This paper deals with relatively large scale subdivisions of volume present in comminution phenomena. Hence, the effects of cracks, lattice dislocations, misplaced atoms, etc., are smoothed out in an average, constant for each relatively large subdivision of volume. This attitude is supported by experience. If two 10 cc samples of the same ore were ground identically, the same product would be obtained. However, if two samples, each a cubic micron, were conceived to be broken, then one sample might contain a crack and the other not, hence a different product would be obtained. Experience shows that ordinary samples used in comminution, behave as though no irregularity existed

Jan 1, 1964

By A. S. Grove, C. T. Sah, E. H. Snow, B. E. Deal

A summary of- several recent investigations in to the properties of the metal-oxide-silicon system is presented. A major portion of these studies makes use of the MOS capacitance-z)oltage method of&apos; analysis. The particular areas of investigation which are reported include: 1) a general survey of the electvical properties of thermally oxidized silicon surjbccs; 2) a study of ion migration through silicon dioxide films ; 3) measurements of electron and hole mobilities in surface inversion layers; 4) a study of impurity redistribution due to thermal o.ridatiotz; and 5) measurements of the rates of oxidation oj-heavily doper7. silicon. THE importance of the metal-oxide-semiconductor (MOS) system in the semiconductor industry is well-known. In addition to its importance in the "planar" device technology,&apos; the MOS structure is now also used in the fabrication of active solid-state devices. Consequently, extensive efforts have been made recently to obtain a better understanding of the characteristics of this system. A summary of some studies of the MOS system conducted in our laboratories during the past year is presented. For the most part these studies used silicon as the semiconductor, along with silicon dioxide and aluminum as the other two components of the system. Since the MOS capacitance-voltage method of analysis was used extensively in these studies, we will first briefly describe its nature and consider some of the possible causes of deviation of experimental observations from the simple theory. We will then outline the various related areas of investigation carried out in our laboratories and will briefly indicate some of the results. It should be noted that the purpose of this paper is merely to provide a brief summary of MOS studies. More detailed discussions of the various areas of investigation are given in the references cited. PRINCIPLES OF THE MOS C-V METHOD OF ANALYSIS&apos; A sketch of the MOS structure is shown in the upper portion of Fig. 1. In this case the insulating film is Si02 and the semiconductor p-type silicon. If a large negative bias is applied to the metal field plate, holes are attracted to the silicon surface. The silicon then behaves much like a metal and the capacitance measured is that of the oxide layer alone, Co. If a small positive bias is applied to the aluminum, holes are repelled and a region depleted of majority carriers is formed at the silicon surface. This depletion I-egion adds to the width of the dielectric and the measured capacitance begins to drop. With increasing positive bias, the width of the electrical depletion region increases. At some large positive bias an inzevsion regiotr is formed at the surface and additional charges induced in the silicon appear in the form of electrons in this narrow inversion region. Thus the depletion-region width approaches a maximum value and, consequently, the capacitance reaches a minimum value and then either levels off or rises again depending on the measurement frequency and the rate of equilibration of the minority carriers in the inversion layer.3 Band diagrams, along with the corresponding charge distributions, are shown in Fig. 1 for the above bias conditions. If minority carriers cannot accumulate at the surface to form an inversion region, the depletion-region width continues to increase with increased positive bias and the capacitance drops toward zero as in a reverse biased p-n junction. The effect of a work-function difference $hs between the metal and the silicon, and of surface charges per unit area Qss located at the oxide-silicon interface, is simply to attract charges in the silicon much like the applied bias. It can be shown that this results in a parallel shift of the capacitance-voltage characteristic along the voltage axis by an amount corresponding to AV = -$bIs + Qss/Co. Theoretical curves have been calculated4 giving the capacitance of the MOS structure C normalized to the oxide capacitance Co vs the quantity VG here VG is the voltage applied to the metal field plate. In Fig. 2 such calculations are shown as points for a particular oxide thickness and bulk impurity concentration for a p-type semiconductor. (For an n-type semiconductor the curves would be mirror images of these.) All three cases, i.e., low frequency. high frequency, and depletion, are indicated. Also shown in the figure are recorder tracings of the characteristics of actual devices. These characteristics have been shifted along the voltage axis to compensate the effect of surface charges and work-function difference. It is evident that agreement between experiment and theory is good. The nature of this shift along the voltage axis is

Jan 1, 1965

By Victor deWolfe

INTRODUCTION The Henderson Mine, located near Empire, Colorado, utilizes a continuous panel caving system to extract ore as one of the world&apos;s major producers of molybdenum. Any mine using a caving-by-gravity technique of mining must rely on closely controlled draw of the caved ore. This control is essential to insure proper caving action, to avoid damaging load concentrations of weight and to minimize the dilution of ore with waste material. Henderson is no exception. Draw control is a major factor in all production planning, from long- range plans to short-range and day-to-day ore scheduling. Draw control is reviewed constantly and administered daily in an effort to optimize production efficiency, ore recovery, and cave management. MINING METHOD The cave at Henderson is massive, moving slowly through large panels that are 244 m (800 ft.) wide by 610 m (2,000 ft.) long. Generally two cave areas are drawn at one time. The areas under active draw vary in size but can be as large as 244 m (800 ft.) by 244 m (800 ft. ) containing 400 draw points. Each draw point contains 45,360 mt (50,000 st) on the average and takes about two and one half years to exhaust. A complete panel is worked for seven to ten years. No pillar exists between panels, but rather a buffer zone of broken ore, or "static face," is left in each panel to be drawn with the adjacent, yet-to-be-caved panel in efforts of minimizing dilution of a working area from an exhausted one. (Figure 1) Production drifts are driven on 24.4 m (80 ft.) centers through the ore body. Between the production drifts are funnel-shaped draw bells on 12.2 m (40 ft.) x 24.4 m (80 ft.) centers to receive ore from the cave. Each bell is accessed by two draw points, one from the production drift on either side, thus forming a 12.2 m (40 ft.) x 12.2 m (40 ft.) draw pattern. Extraction of the ore is by rubber-tired, 3.8 m3 (5 yd3) load-haul-dump equipment. The LHDs then tram the ore a maximum of 49 m (160 ft.) to ore passes. Cave initiation and bell development are done from the undercut drifts which are parallel to and 17 m (55 ft.) directly above the production drifts. Longhole rings are drilled and blasted from the undercut drifts to define the bells and establish the undercut for caving. (Figure 2) DRAW CONTROL Since the cave line at Henderson is constantly advancing, it is necessary to be continually initiating new cave at one end while exhausting it at the opposite end. There must exist, therefore, an angle on the ore-waste contact in the broken rock from initiation to exhaustion. The basic concept of draw control is to keep this angle as smooth and even as possible, particularly at the time of exhaustion. If this is achieved, draw points are exhausted more or less in a line, avoiding pockets of remaining ore surrounded by exhausted areas. These pockets would cause spotty ore extraction at the time of exhaustion, increasing the amount of dilution occurring while introducing the potential for significant weight problems in the production area. To arrive at the desired angle on the ore- waste contact, maximum tonnage percentages are assigned to each row of draw points increasing at 10% or 15% increments (depending on cave size and velocity of draw) working away from the cave line. The available tonnage indicated by these percentages is the maximum allowable tonnage to be extracted from each draw point until the available tonnage percent- age is increased. As the cave moves, these percentages increase for each draw point regularly. However, in general the tonnage drawn from each draw point is kept at about 50% of this allowable maximum in order to maintain adequate available tonnage in the cave to sustain production for seven months if cave initiation were to cease. This available tonnage cushion is a safeguard built into the draw control program at Henderson to accommodate fluctuations in the rate of cave advance. When draw points move past the row of 100% tonnage availability, they are drawn past the desired 50% at the same increments per row until exhausted. (Figure 3) To achieve proper draw control, the number of LHD buckets to be taken from each draw point is assigned daily. The actual buckets taken, which may at times deviate from the

Jan 1, 1981

By E. Teghtsoonian, E. M. Schulson

The tensile behavior of bcc ordered P&apos; AuZn single crystals (CsCl structure) has been investigated under varying conditions of temperature, composition, and orientation. Between -0.2 and 0.4 T, multi-stage hardening occurs fm stoichiometric and nonstoichio-metric crystals oriented near the middle of the primary stereographic triangle. At higher and lower temperatures, parabolic type hardening occurs, followed by work - softening at the higher temperatwes. Deviations from stoichiometry give rise to increased flow stresses. Multi-stage hardening was observed for most orientations, except along the [loll-[lll] boundary and near the [001] corner of the stereo -graphic triangle, where parabolic type hardening occurs. Along two slip systems, (hk0)[001] and (, operate simultaneously while in the [001] comer, slip occurs mainly on the system. Electron microscopy of deformed crystals revealed bundles of edge dislocations forming walls approximately Perpendicular to the glide plane. In general the plasticity of 4&apos; AuZn closely resembles the plasticity of bcc crystals. In recent years, considerable interest has arisen concerning the mechanical properties of the CsCl type intermetallic compounds Ag Mg,&apos;- Fe co,&apos; and Ni Al.&apos;-&apos; The compound P&apos;AuZn is structurally similar. It has a low and congruent melting point of 725"~,&apos;" remains ordered up to the melting point,16 and pos-esses a range of solid solubility from 47.5 to 52.0 at. pct Au at room temperature.15 The present paper reports the results of an investigation on the general tensile behavior of material in single crystal form. Some dislocation configurations characteristic of the deformed state are also reported. The results of a detailed study of the slip geometry in AuZn are presented in a separate paper.17 PROCEDURE Alloy preparation, crystal growing techniques, and the procedure followed in selecting specimens of minimum composition variation are reported elsewhere.17 Dumb-bell shaped tensile specimens were prepared by carefully machining single crystals in a jewellers&apos; lathe to a gage length of 0.80 in. and diam of 0.090 in. Back-reflection Laue X-ray patterns and room temperature tensile tests revealed that machining damage could be eliminated by electrochemically polishing 0.005 in. from the machined surface followed by annealing at 300°C for 1 hr. Specimens were polished in fresh 5 pct KCN solution (40°C, 12 v). Experiments were performed by gripping specimens in a self-aligning pin-chuck and threaded collet system, then straining in a floor model Instron tensile machine. All tests were performed in duplicate. Experimental variables included temperature, composition, and orientation. Unless otherwise stated the strain rate was 2.5 x 10"3 per sec. Liquid testing environments included nitrogen (WOK), nitrogen cooled petroleum ether (133" to 293"K), and silicone oil (293" to 488°K). Resolved shear stress-shear strain curves were electronically computed from autographically recorded load-elongation curves. Stress and strain were resolved on the macroscopic noncrystallographic (hkO) [001] system operative under the specific test conditions of temperature, strain rate, and orientation reported earlier.17 RESULTS The temperature dependence of the work-hardening curves is shown in Fig. 1 for gold-rich crystals of 51.0 at. pct Au oriented near the center of the stereo-graphic triangle. Over the range of intermediate temperatures from -200" to 400°K, they are very similar to those classically observed for fcc metals (reviewed by Nabarro et al.).&apos;&apos; The beginning of deformation is characterized by a region of decreasing hardening rate, stage 0, which is followed by a region of low linear hardening, stage I, and then a region of higher linear hardening, stage 11. At the higher temperatures, stage 111 is observed, a region of decreasing hardening rate. Over the intermediate temperature range, the extent of stage 0 and of the slow transition between stages I and I1 decreases with increasing temperature. Total ductility is large, often greater than 300 pct shear. As the temperature is either increased or decreased, the extent of stage I is decreased, giving rise to parabolic type flow and reduced ductility. Similar temperature effects have been reported for bcc ~r~stals.~~-~~ Below -14O°K, hardening is terminated in brittle fracture while above -400°K. initial hardening is followed first by work-softening and then by chisel-edge type ductile fracture. Stoichiometric (50.0 at. pct Au) and Zn-rich (51.0 at. pct Zn) crystals were also tested from 77" to -500°K. The effect of composition on the flow behavior is illus-

Jan 1, 1970

By G. H. F. Gardner

The velocity and attenuation of elastic waves in sandstones were measured as a function of both pressure and fluid saturation. A large change occurs in these quantities if water is added and the rock is not compressed, but the change is small if the rock is subjected to a large overburden pressure. Measurements were made by vibrating cylindrical samples in both the extensional and torsional modes at frequencies up to 30,000 cycles/sec. Formulas were derived which enable the attenuation of dilatational waves in dry rocks to be deduced from the data. Similar experimental methods were used to investigate the properties of unconsolidated sands. Velocities were found to vary with the 1/4 power of the overburden pressure and attenuations to decrease with the 1/6 power. The effects of grain size, amplitude and fluid saturation were studied. Formulas by which the effects produced by a jacket around the sample may be calculated were derived. The practical application of these results to formation valuation is discussed. INTRODUCTION The attenuation of elastic waves in the earth has been of interest to the seismologist and geophysicist for many years, but only recently to the petroleum engineer. Engineering interest has been brought about by the success of velocity logging devices, for it is possible by modification of these instruments to measure the attenuation of sound waves in addition to their velocity and, hence, deduce the mobility of formation fluids as well as the porosities of the rocks which contain them. The main problem is to decide whether field measurements can be made with sufficient accuracy to be of practical use. This problem can only be solved after we know the magnitude of the attenuations which are typical of the earth at various depths. The logarithmic decrement of a fluid-saturated rock is the sum of a "sloshing" decrement and a "jostling" decrement, the former caused by the mobility of the fluid contained within the rock and the latter by the granular framework of the rock. Sloshing decrements can be calculated&apos; using Biot&apos;s theory, but the jostling losses are less well understood. The present paper reports an experimental investigation of jostling losses in consolidated and uncon- solidated sands, particularly with respect to the effect of overburden pressure and fluid saturation. Born&apos; showed that the decrement of a sandstone may increase dramatically when only a few per cent by weight of distilled water is added, and that the additional loss is proportional to the frequency of vibration. His measurements were made with no compressive stress on the framework of the rock. M. Gondouin3 investigated similar phenomena for fluid-saturated plasters but also did not compress the samples. In the present paper it is shown that compression of the framework reduces this effect, so that at depth the jostling decrement of a sandstone may be expected to be almost independent of fluid saturation and frequency. Decrements for many sedimentary rocks have been given by Volarovich,4 but all for the state of zero overburden pressure. Anomalously low velocities have been logged in shallow unconsolidated gas sands. Results of the present investigation confirm that these velocities are not caused by correspondingly high attenuations, because the jostling decrement in a packing of sand grains is small and much less than in a consolidated sandstone at the same depth. Velocities in sands have been measured by Tsareva5 and by Hardin6 as a function of pressure, but the corresponding decrements do not appear to have been measured previously. The widely used "resonant bar method" of measuring velocities and decrements was employed. Comments on variations of this technique have recently been published by McSkimmin.7 The main novelty of the present technique was the application of pressure to the samples. It was found possible to do this by placing the apparatus inside a pressure vessel, provided the conditions leading to large additional losses were avoided. These conditions are discussed below. EXPERIMENTAL TECHNIQUE Cylindrical samples were caused to vibrate in both the extensional and torsional mode of vibration and the amplitude of vibration was measured as a function of frequency in the neighborhood of a resonant frequency. The resonant frequency, fr, is related to the corresponding elastic modulus by the formulas where E and N are Young&apos;s modulus and the modulus of rigidity, p is the density of the sample, and A the wavelength of the vibration.

Jan 1, 1965

By J. F. Radavich

Many heats of steel of low carbon value have been known to produce brittle pieces of steel. The brittleness is believed to be due to the impurities located within the grain boundaries. Such brittle steels have been examined with an optical microscope to ascertain the nature and the amount of the impurities present at the grain boundaries. Due to the relatively low resolving power of the optical microscope, the impurities are not visible in fine detail. The writer obtained some sheet steel and proceeded to determine the location of the impurities and to show the application of the electron microscope to the study of grain boundaries. One sample was known to be capable of becoming embrittled, whereas another sample was believed to be much less susceptible to embrittlement. Treatment of Specimens The specimens were embrittled by annealing above the A3 point under mildly oxidizing conditions. One piece of ingot iron could not withstand a 90" bend, whereas another piece of ingot iron was not affected and could withstand a 90" bend. The brittle piece was then annealed at a high temperature in a hydrogen atmosphere. The annealed ingot iron was termed cured and could withstand a 90" bend very easily. The three specimens examined will be designated as brittle, good. and cured in the discussion that follows. Procedure The sizes of the specimens were as follows: one piece of brittle ingot iron-3/8 by 35 in.; one piece of good ingot iron-96 by 1/8 in.; one piece of cured ingot iron-36 by 54 in. The specimens were too small to be polished by hand and therefore were mounted in bakelite. The polishing procedure was carried out in the conventional manner with the use of 1/0 through 3/0 papers, and the final polish was done with alumina on a billiard cloth. The specimens were then etched in a 4 pct solution of picral in alcohol, and then they were examined through an optical microscope. An area was chosen that showed distinct grain boundaries, and an effort was made to keep near this area when pulling the replicas REPLICA TECHNIQIJE The replica technique used in the preparation of the replicas for examination under the electron microscope is described in Electron Metallography.&apos; It consists essentially of the following steps: 1. Obtaining a suitably etched specimen. 2. Applying a swab of ethylene di-chloride on the surface. 3. Applying a formvar solution on the surface. 4. Placing a screen on any desired spot. 5. Breathing on the fornivar layer. 6. Applying scotch tape on the screen and film. 7. Pulling the film and the screen up with the Scotch tape. 8. Separating the screen from the Scotch tape. This replica technique is very similar to the one described by Harker and Shaefer. However, with the added step, the percentage of replicas removed is very much higher regardless of the length of the time from the etching of the specimen to the actual pulling of the replica. The replicas were then shadow cast with manganese at a filament height to replica distance ratio of 1 1/2:7. This produced a very high contrast replica for use in the electron microscope. One of the dificulties encountered with this study was the restricted area of the specimen. The width of the specimens was the same as that of the 200 mesh nickel supporting screen. In order to increase the effective area, the screens were cut down as shown in Fig 1. The arrow indicates the direction in which the replica was pulled. This operation made it possible to obtain a large percentage of good replicas. Fig 3 shows an electron micrograph of a brittle piece of ingot iron and a grain boundary that was polished mechanically. The surface is very rough probably due to the incomplete removal of the flowed layer by the picral etchant. The grain boundary does show evidence of impurities. It was decided to electropolish the specimens to obtain a much smoother surface than the one obtained by mechanical polishing. ELECTROPOLISHING The specimens were cut in half to expose the metal on the back side. The exposed metal had sufficient area to make good electrical contact and electropolishing was carried out easily. The conditions for electropolishing were 0.9 amp, 35 volts, and 25 sec. in an electrolyte composed of 850 cc of ethyl alcohol, 100 cc distilled water, and 50 cc of perchloric acid. The polished specimens were then etched in the 4 pct picral solution for a shorter time than was necessary for

Jan 1, 1950

By M. Kaitz

A continuing discussion in both the petroleum engineering and economic literature is directed to the difficulties encountered in the use of discounted cash flow rate of return (DCF) as a measure of investment worth. Although useful in most Instances, DCF has been criticized because it is time-consuming in its trial-and-error solution, theoretically invalid, not adaptable to cash flow streams yielding multiple return solutions and not entirely reliable in selecting between mutually exclusive investments. The theoretical invalidity of DCF stems from the reinvestment assumption implicit in the calculation that earnings are reinvested at the DCF rate. The emerging consensus of the economic literature is that the net present value or the present worth of the net cash flow stream (discounting at the average opportunity or cost of capital rate) is more correct and reliable. Other criteria proposed have been ratios of net present value divided by initial investment or by present value of all investments in a project. All of these criteria are simple to determine and explicitly assume a reinvestment rate for [he income generated by a project. This paper develops and discusses a theoretically valid profitability criterion which is simple to compute and retains the appeal of a percent return on investment. It is called "percentage gain on investment" or PGI. It measures the gain an investment is expected to realize over like capital invested in the average opportunity and explicitly considers reinvestment potential. Why add another Concept to the large array of investment criteria now available, any one of which, or perhaps a combination of several, appears to embrace the form&apos;s (or individuals) objectives? The answer is that not one of the existing criteria provides both a readily comprehensible and theoretically valid measure of risk coverage that has general application. The proposed PGI does fulfill these requirements. INTRODUCTION An ancient expression warns that "one must yield to the times" — there are better ways of doing things. A review of the petroleum engineering and economic literature on one topic alone, measurement of investment worth, certainly is witness to this truth. In use for a number of years, the DCF has recently received attention, directed mainly to its theoretical invalidity. Several alternatives to DCF have been proposed to provide a valid, simply determined criterion to describe investment worth and to overcome the criticisms previously mentioned. This paper introduces another method called percentage gain on investment (PGI) and is proposed as but one of several yardsticks that should be used in making investment decisions. MEASURES OF INVESTMENT WORTH This paper will consider only those criteria which give weight to the time pattern of future earnings. These criteria are usually compared with an average opportunity rate or cost of capital of the firm to judge the relative worth of the investment. For purposes of demonstration, a 9 percent average opportunity rate will be used throughout this paper. Implicit in the DCF calculation is the assumption that earnings are reinvested at the DCF rate. Some argue, though, that there is no reinvestment assumption,&apos; that the DCF rate is simply that maximum rate of interest one can pay on the investment over the life of the project and break even. The determination of DCF is accomplished by discounting the net cash flow stream at that rate (DCF) which will yield zero. The question is: why should reinvestment potential be explicitly considered in calculating return on investment or other criteria measuring economic worth? Perhaps the answer lies in consistent or equal treatment of future cash flow. It appears entirely illogical to give different present worth value to $1 received, say, 10 years from now, which is the circumstance resulting in comparing projects with different DCF returns. In fact, $1 received in 10 years has the same value regardless of which project generated the income. DCF return thus favors investment projects which are expected to provide early income as compared to those providing long-term income. Not surprisingly, the controversy on reinvestment assumption is an old one. Hoskold&apos; discussed this same problem in 1877. He considered the future income from mineral properties as an annuity or a series of fixed future payments. (As will be demonstrated, his equation can be modified for variable income.) Prior to Hoskold, the value or what one could pay for the mine, was determined with the use of standard, single-interest tables. Here is what Hoskold said with regard to these tables. "This table, and others of its kind, to be found in most works on annuities, is constructed correctly according to

By J. B. Campbell, V. V. Valleroy, B. T. Willman, L. W. Powers

A steam drive of heavy oil was field tested in a shallow, low oil-saturation formation near Deerfield, Mo. The pilot was conducted in the Warner formation, a sandstone containing an 18&apos; API oil having 1,000-cp viscosity at the 60F origind reservoir temperature. The formation war at a depth of 160 ft. Steam was injected into nine input wells arranged in an array of inverted five-spot patterns. In the completely confined center pattern, 14 temperature observation wells were installed to obtain thermal data and observe test progress. Late in the test, slugs of ammonia were injected to trace the flow paths of injected fluids. From the test area about 7,000 bbl of oil were produced. Data were obtained on areal and vertical temperature distribution, steam front advance, reservoir fluid movement and terminal saturations. This field test of a steam drive (I) demonstrated the feasibility of the method, (2) confirmed that the low residual oil saturations observed in the laboratory are obtained in the steam-swept region in the field and (3) provided recovery and conformance data for one set of field conditions. INTRODUCTION The Deerfield steam drive pilot test was conducted in a shallow sandstone containing 1,000-cp oil. The venture was undertaken cooperatively by the research and production departments of Carter Oil Co., which organizations have since been consolidated into Esso Production Research Co. and Humble Oil & Refining Co.. respectively. The production department was interested in steam injection at Deerfield because it appeared to be the most promising method of commercially producing this heavy oil deposit. The research department was interested in applying the new recovery method and in evaluating its performance in the field. At the time the test was begun, the initial oil saturation was not well known. Subsequent air coring and early pilot results confirmed that there was too little oil in place for profitable commercial exploitation by steam. Pilot termination at that time, however, would have been premature for evaluating field performance of the process, and the tert was continued to obtain additional data on steam injection as a recovery method. The test was located in Vernon County, Mo., about 10 miles north of the town of Deerfield and only a few miles from the Kansas border. The pilot site was selected as typical of the area. The location represented neither the highest nor the lowest oil saturation region in the acreage under lease in 1954. The steam drive was conducted in the Warner sandstone of Lower Pennsylvanian age. At the test site the top of the Warner occurs at about 160 ft subsurface and the formation is a fine- to medium-grained micaceous sandstone that dips gently to the northwest at the rate of 12 to 15 ft/ mile. A cross-section and permeability profile of the test location are shown in Fig. 1. At the pilot location the average total thickness of the Warner formation is about 43 ft, but the effective thickness for steam drive is 26 ft. Two distinct types of hydrocarbon saturation are apparent. The lower portion of the total sand, averaging about 17 ft thick, contains a very heavy asphaltic material that will not flow under the influence of a steam drive. This bottom interval, referred to as a dead oil residue, was not considered as part of the net sand undergoing steam exploitation. The initial formation and fluid properties of the upper 26 ft in the test area are summarized in Table 1, and variation of oil viscosity with temperature is shown in Fig. 2. Imbibition tests on preserved core samples taken at the end of the pilot test showed that the Warner sandstone was then neutral or slightly water-wet. Initially, the reservoir may have been strongly water-wet as indicated by low relative permeability to water during both water injection testing and early steam injection. PRIOR HISTORY Initial production tests of wells at the pilot site produced water with only a faint show of oil. No gas was produced except at Well 7-W in the pilot area and at another well about 1/3 mile northeast of the pilot. Prior to the start of the steam drive, a two-well water injection test and a two-well air injection test were conducted. No oil was produced by either. Water was pumped into Well I-W in the northeast corner of the pilot area with simultaneous production from Well 1 (Fig. 3). The air-injection tat was run at input Well 9-W and its offset, Well 2, in the southwest corner. Air and water injectivities were about the same when corrected for viscosity and pressure differences.

By L. Delisle

This work represents the first step of an attempt to test the applicability of the electron microscope to the study of subgrain structures in copper. Observations on annealed and deformed single crystals and polycrystalline samples of copper are described. IN the course of study of the structure of fine tungsten wires and tungsten rods with the electron microscope, well defined subgrain structures were observed. The size, size distribution, and orientation uniformity of the etch figures varied widely in different samples. Figs. 1 and 2, electron micrographs of a tungsten wire and of a tungsten rod, respectively, are illustrations of the difference in size and size distribution of the etch figures in different samples of the same metal. The observed differences, as pointed out in a previous paper,&apos; appeared to be related to the heat and mechanical treatments of the samples. They were also consistent with the results reported in the literature on the mosaic structure of metals.&apos; For that reason a program of research was initiated in an effort to obtain more systematic evidence of the possible relation of heat and mechanical treatments to the subgrain structure of metals as observed in the electron microscope. The purpose of this paper is to present observations made on the effect of cold work on the subgrain structure of copper. Procedure Starting Materials: Copper was the metal studied because it can be obtained in a high degree of purity, much information is available in the literature on its properties and its response to cold work and heat treatment, it shows no allotropic change, and it is sufficiently hard to be handled without great difficulty. Two groups of specimens were used: 1—single crystals cast from spectroscopically pure copper and 2—polycrystalline samples of oxygen-free high conductivity copper. Single crystals were studied because it was hoped that the elimination of a number of variables, such as grain boundaries, orientation differences, degree of purity, would simplify the problem and perhaps permit a better understanding of the phenomena that would be observed. The polycrystalline samples were designed to give a general picture of the changes considered. The single crystals were made of copper which analyzed spectroscopically to better than 99.999 pct Cu. They were cast in vacuum, by the Bridgman method, in crucibles made of graphite with a maximum ash content of 0.06 pct. The mold design is shown in Fig. 3. It permitted casting crystals of the size and shape required for the experiments, so that the danger of introducing cold work in the original samples by cutting or other machining would be eliminated. The polycrystalline samples were pieces, 3/4 in. long, cut from a rod of oxygen-free high conductivity copper, % in. in diameter. A flat surface, 1/4 in. wide, was milled along the rods, polished, and etched. The samples were then annealed in vacuum at 850°C for 1 hr. Polishing and Etching: Work previously done on tungsten,&apos; polished mechanically and etched chemically," had shown that: 1—the general appearance of the etch figures of a given sample was not altered by repeated polishings and etchings under similar conditions; 2—variations in the time of etching and the concentration of the etchant changed the definition of the etch figures, but did not alter their general size nor orientation distribution within the limits of observation. Further work confirmed the reproducibility of the subgrain structures observed in, 1—single crystals and polycrystalline samples of copper when polishing and etching were repeated under similar conditions, and 2—specimens of tungsten and polycrystalline copper when electrolytic polishing and etching were substituted for mechanical polishing and chemical etching, respectively. On the strength of these observations, it was felt that, if conditions of polishing and etching were kept constant, changes observed in the subgrain structure of a sample upon deformation and annealing would be attributable to such treatments. For that reason the conditions of polishing and etching were kept as constant as possible. The single crystals were polished electrolytically in a bath of orthophosphoric acid in water, in the ratio of 1000 g of acid of density 1.75 g per cc to 1000 cc of solution, under a potential drop of 1.6 to 1.8 V. Electrolytic polishing was selected to prevent the formation of distorted metal in polishing. The same samples were etched by immersion in a 10 pct aque-

Jan 1, 1954