Search Documents

MELTING in a hearth and in a basket is almost the same thing; charcoal and bellows are used for both. These are made large or small according to the quantity that you wish to melt, and more or less force of the bellows is applied according to the need. The hearth, bowl, or refinery, as the masters call it, is made with bricks and clay like a forge; opposite the middle of the tuyère a round hole is made like a wash- basin, wide at the mouth and narrow at the bottom where there is a hole for emptying it. A little iron plug is put in this so that it may be closed when necessary, or a piece of brick that is cut and pointed. Now, when this has been well covered all over with ashes, and the tuyère has been arranged so that the wind strikes the center, it is first filled with charcoal and baked very well so that the metal may be kept hot once it has been melted. When it has been baked, it is again filled with charcoal, fire is applied, and it is left to inflame gradually by itself. Then put one or

Jan 1, 1942

By H. C., Weed

THE Inspiration Consolidated Copper Co., located in the Globe-Miami district at Inspiration, Ariz., became a producer of copper in 1915. From 1915 until 1948, 116,278,000 tons of ore were produced from underground workings by the block-caving system of mining. In April 1948, open-pit ore production from portions of the ore body was started and earlier plans for block-caving in these areas were abandoned. The change was due to economic conditions. Prior to the war years Inspiration's costs of mining by block-caving methods were comparable to, or even below, the costs of some open-pit operations. By the end of the war these conditions were completely changed. Underground costs had risen steadily for labor and supplies. Difficulty was being experienced in finding competent underground men and underground development was almost hopelessly behind schedule. On the other hand, open-pit costs had been sharply reduced through the development of modern earth-moving machinery and plenty of labor was available for surface work.

Jan 1, 1950

By N. M. Levine, W. M. Fassell

In this laboratory study the problem of aerative conditioning to separate chalcopyrite and pyrite from cobaltite was simply effected with a sulfy-drate collector and pH by proper choice of mixing variables. It was shown that for a given impeller type and geometry, selectivity is governed by a relation between the Froude number and a modified power number multiplied by impeller tank-diameter ratio wherein the power factor is expressed as power intensity or horsepower per unit volume of contained pulp. In applying this approach to other systems, other criteria would be used in place of selectivity. The purpose of this paper is to show how solutions to problems of aerative conditioning of pulps may be enhanced by development of mixing parameters in a given system through application of a fundamental approach to agitation. The term aerative conditioning refers to a combination of processes which include mechanical gas dispersion, gas solution in the liquid phase, transport of dissolved gases and other reactants to the reaction zone, physico-chemico interaction within the liquid phase and/or at the solid surface, and, finally transport of products away from the reaction zone. Important elements of these processes are: 1) Air-Liquid Contact Area: This area limits the rate at which a given gas may dissolve in a given liquid. Thus, an aerative conditioning system should be designed for optimum gas dispersion. 2) Diffusional Factors: Such factors as concentration gradient, diffusion film thickness, and diffusivity may govern the rate at which dissolved constituents are transferred to the zone of interaction at a particle surface and the rate at which reaction products are transported therefrom. Thus, such systems would be designed for maximum turbulence to effect good mixing and to minimize diffusion film thickness. 3) Solids-Turbulent Liquid Contact Time: Sufficient time should be provided for interactions to take place at solid-liquid interfaces. The region of maximum turbulence for mechanical systems being in the impeller zone, it is important that the number of passes through this zone be optimized. This means that a proper balance should be struck be- tween flow and turbulent power to achieve optimum conditioning time. Too high a flow rate and insufficient turbulence may be improper if the reactions are diffusion limited. On the other hand, too much turbulent power and too little time in the reaction zone due to insufficient flow power may also be limiting. 4) Geometry and Impeller Design: The design of the conditioning vessel and the type of impeller employed are important factors in determining the relative distribution between flow and turbulent power and the degree of gas dispersion. Maximization of power input would be to no avail if the distribution between flow and turbulent power were improper or gas dispersion insufficient. 5. Physical Chemistry: The reactions between dissolved gases and solid surfaces may be limiting if gas solubility, activation energies and/or free energies are not favorable. J. H. Rushton, E. W. Costich, and H. J. Everett' have summarized the mechanical elements of the combination of processes listed in a fundamental mixing equation. This equation was derived by dimensional analysis and relates the physical variables for a single impeller, centered in a cylindrical, vertical axis, flat-bottomed tank. It is, where T is tank diameter in feet; H is liquid depth in feet; C is height of impeller off tank bottom in feet; S is pitch of impeller; L is length of impeller blades; W is width of impeller blades; R is number of baffles; D is impeller diameter in feet; Np is the power number, pg/dN3D5; P is power in foot-pounds per second; g is the gravitational constant, feet per second; J is width of baffles; B is number of im-

Jan 1, 1961

By C. R. McEwen

This paper presents a technique for calculating the original amount of hydrocarbon in place in a petroleum reservoir, and for determining the constants characterizing the aquifer performance, based on pressure-production data. A method for doing this based on a least-squares line-fitting computationwus proposed by vanEverdingen,Timmerman andMcMahon in 1953. We found that their method would not work when there is error in the reservoir pressure data — even moderate error. The technique presented here appears to give reasonable answers when pressure data are uncertain to the degree expected in reservoir pressure determinations. The major change introduced in the present analysis is to limit the least-squares line-fitting to yield only one constant— the amount of hydrocarbon in place. The water-influx constant is then taken as proportional to the oil (or gas) in place. The constant of proportionality can be computed from estimates of effective compressibility and reservoir water saturation. It is also pointed out that the commonly used least-squares analysis assumes all of the uncertainty to be in the dependent variable. The material balance should be arranged so that this condition is fulfilled if correct inferences are to be made from statistical calculations. Examples are shown of the application of the new technique to gas reservoirs—both hypothetical and real—and to the oil-reservoir example of van Ever-dingen, Timmerman and McMahon. INTRODUCTION The amount of hydrocarbon originally in place in a petroleum resenroir can be estimated by means of the material-balance calculation. Simultaneous observations of pressure and amounts of produced fluids are required, together with the PVT data applicable to the reservoir fluids. If water encroachment is occurring, it is desirable to try to infer the behavior of the aquifer, as well as the original hydrocarbon in place, from the pressure-production data. This imposes additional demands on the method of calculation, and uncertainty in the data can result in large uncertainty in the answer. In addition, if the size of a gas cap is to be established, the whole problem becomes indeterminate, as pointed out by Woods and Muskat.1 Brownscombe and Collins2 simulated a gas reservoir and its aquifer on a reservoir analyzer and derived quantitative information on the effect of uncertainty in pressure and aquifer permeability on computed gas in place. Among the various techniques of estimating the performance of an aquifer, the method of van Ever-dingen and Hurst? based on compressible flow theory, seems to have been the most generally successful (see Ref. 4, for example). In this paper we shall confine ourselves to their representation of the aquifer. In 1953,' van Everdingen, Timmerman and McMahon5 introduced a statistical technique for deriving the amount of oil originally in place and the parameters which describe the aquifer. (We shall refer to this technique as the "VTM method'', as suggested by Mueller.6) Their example reservoir had no gas cap. It has been our experience that the VTM method gives a reasonable answer when the data are very accurate, but that inaccuracy (particularly in pressure) can cause the method to break down. The effect was first observed in gas reservoirs, but has since been seen in oil reservoirs also. In this paper we present another statistical method which has been successful in achieving a reasonable answer where the VTM method has failed. In the new method, one less parameter is derived from material-balance computations. It is assumed that values can be established for effective compressibility in the aquifer and reservoir water saturation independently of the material-balance calculation. The water-influx constant can then be obtained from these data and the quantity of hydrocarbon in place.

By C. S. Matthews, W. R. Strickler, M. Prats

To gain an insight into the deviations from normal behavior caused by dip. a mathematical study has been made of single-fluid five-spot floods in dipping reservoirs for flooding techniques likely to be employed in such cases. Streamline patterns and sweep eficiencies are computed for a number of different cases, and the results, where possible, are presented in graphical form. An increase in sweep efficiency can be obtained by moving the central injection well updip, and the opti~num position is found for several cases. The applicnhility of the results is discussed. INTRODUCTION Because of the growing importance of water flooding in California and in other locations where large angles of dip are quite common, it is important to know at least some of the effects of dip on the performance characteristics of water floods in such instances. In the past, only reservoirs having relatively mild dips have been wateflooded. Consequently, little is known today as to the effect of dip on the performance characteristics of a flooding operation from the point of view of field experience; and because of the complexity of the problem, little, if anything, has been done on the problem from the point of view of theory. There are two distinct aspects of the role played by gravity in a flooding operation. The first is that of the preferential segregation of the reservoir fluids at the time when the flood is to commence. The second is the role that it may play in affecting the performance of the reservoir during the course of the flooding operation. In this paper the reservoir is assumed to be filled with only one fluid, so that there is no preferential segregation of reservoir fluids at the commencement of the flood. Moreover, the injected and in-place fluids are assumed to have the same characteristics, so that in effect we consider only a single-fluid problem. Single-fluid results have been used as a guide in the past for predicting and interpreting the behavior of regular five-spot floods in horizontal reservoirs. It is thought that single-fluid results might likewise be used for five-spot patterns in dipping reservoirs as a norm to which the effects of the multifluid problems might he referred.* In confining our attention to the single-fluid case, we therefore do not consider the preferential segregation of injected water in the oil bank, or of the oil bank in the unflooded zone. Instead we consider the deviation from the behavior of a horizontal five-spot which is caused by use of flooding techniques (i.e., the flow potential pattern at injection and production wells) likely to be employed during flooding operations in dipping reservoirs. In other words, we consider effects which are caused by the pattern of injection and production well pressures which we choose to apply. Performance characteristics are given for the five-spot spacing pattern as a function of two parameters; one is the ratio of gravity head to potential difference bctwccn the high updip injection well and the low downdip producer of the same five-spot, denoted by Rg; the second is the orientation of the five-spot pat-

Jan 1, 1956

By I. Javandel, P. A. Witherspoon

The finite element method was originally developed in the aircraft industry to handle problems of stress distribution in complex airframe configurations. This paper describes how the method can be extended to problems of transient flow in porous media. In this approach, the continuum is replaced by a system of finite elements. By employing the variational principle, one can obtain time dependent solutions for the potential at every point in the system by minimizing a potential energy functional. The theory of the method is reviewed. To demonstrate its validity, nonsteady-state results obtained by the finite clement method are compared with those of typical boundary value problems for which rigorous analytical solutions are available. To demonstrate the power of this approach, solutions for the more complex problem of transient flow in layered systems with crossflow are also presented. The generality of this approach with respect to arbitrary boundary conditions and changes in rock properties provides a new method of handling problems of fluid flow in complex systems. INTRODUCTION Problems of transient flow in porous media often can be handled by the methods of analytical mathematics as long as the geometry or properties of the flow system do not become too complex. When the analytical approach becomes intractable, it is customary to resort to numerical methods, and a great variety of problems have been handled in this manner. One such method relies on the finite difference approach wherein the system is divided into a network of elements, and a finite difference equation for the flow into and out of each element is developed. The solution of the resulting set of equations usually requires a high speed computer. When heterogeneous systems of arbitrary geometry must be considered, however, this approach is sometimes difficult to apply and may require large amounts of computer time. The finite element method is a new approach that avoids these difficulties. It was developed originally in the aircraft industry to provide a refined solution for stress distributions in extremely complex airframe configurations. 27 Clough has recently reviewed the application of the finite element method in the field of structural mechanics.6 The technique has been applied successfully in the stress analysis of many complex structures.l, 27, 28 Recognition that this procedure can be interpreted in terms of variational procedures involving minimizing a potential energy functional7 leads naturally to its extension to other boundary value problems. In the field of heat flow, there recently have been introduced several approximate methods of solution that are based on variational principles.2-4, 17 By employing the variational principle in conjunction with the finite element idealization, a powerful solution technique is now available for determining the potential distribution within complex bodies of arbitrary geometry. In the finite element approximation of solids, the continuum is replaced by a system of elements. An approximate solution for the potential field within each element is assumed, and flux equilibrium equations are developed at a discrete number of points within the network of finite elements. For the case of steady-state heat flow, the technique is completely described by Zienkiewicz and Cheung.33 Since the flow of fluids in porous media is analogous to the flow of heat, Zienkiewicz et al. have employed the finite element method in obtaining steady-state solutions to heterogeneous and anisotropic seepage problems. 34 Taylor and Brown have used this method to investigate steady-state flow problems involving a free surface.25 The work of Gurtin has been instrumental in laying the groundwork for the application of finite element methods to linear initial-value problems.12 As a result, Wilson and Nickell have recently

Jan 1, 1969

By J. J. Jakosky

The fundamental hydrodynamic principles governing the production of oil from wells have been carefully studied and evaluated by many investigators. These prior studies are quite complete and cover virtually the entire field of oil production and recovery.' One of the most important measurements in production studies is that of the bottom-hole pressure at various rates of production.2 It is the purpose of this paper to outline the technique and economics of fluid-level measurements as applied to pumping wells and to review briefly other methods for obtaining bottom-hole data. The applications of bottom-hole pressure or fluid-level measurements are based on simple hydrodynamic principles. Under equilibrium conditions, pbh = hd + pg + pch [1] when pbk = bottom-hole pressure, Ib. per sq. in. absolute, h = height of fluid in well, ft., d = average fluid density, Ib. per sq. in. per foot, p, = weight of section of gas column 1 in, square from fluid level to casinghead, Pch = casinghead pressure, Ib. per sq. in. absolute. In words, equation 1 states that the bottom-hole pressure is equal to the weight hd of the column of fluid plus the weight of the column of gas between the fluid and the casinghead plus the casinghead pressure. At moderate rates of flow, the rate of production is approximatcly proportional to the difference between the reservoir, or static bottom-hole. pressure and the flowing or producing pressure; i.e.: Q = c(pr - pbh) [2] pr = the reservoir pressure, c = the constant of proportionality. However, at production rates approaching zero and also at high rates the proportionality may break down. By combining equations 1 and 2 the

Jan 1, 1939

By J. J. Jakosky

The fundamental hydrodynamic principles governing the production of oil from wells have been carefully studied and evaluated by many investigators. These prior studies are quite complete and cover virtually the entire field of oil production and recovery.' One of the most important measurements in production studies is that of the bottom-hole pressure at various rates of production.2 It is the purpose of this paper to outline the technique and economics of fluid-level measurements as applied to pumping wells and to review briefly other methods for obtaining bottom-hole data. The applications of bottom-hole pressure or fluid-level measurements are based on simple hydrodynamic principles. Under equilibrium conditions, pbh = hd + pg + pch [1] when pbk = bottom-hole pressure, Ib. per sq. in. absolute, h = height of fluid in well, ft., d = average fluid density, Ib. per sq. in. per foot, p, = weight of section of gas column 1 in, square from fluid level to casinghead, Pch = casinghead pressure, Ib. per sq. in. absolute. In words, equation 1 states that the bottom-hole pressure is equal to the weight hd of the column of fluid plus the weight of the column of gas between the fluid and the casinghead plus the casinghead pressure. At moderate rates of flow, the rate of production is approximatcly proportional to the difference between the reservoir, or static bottom-hole. pressure and the flowing or producing pressure; i.e.: Q = c(pr - pbh) [2] pr = the reservoir pressure, c = the constant of proportionality. However, at production rates approaching zero and also at high rates the proportionality may break down. By combining equations 1 and 2 the

Jan 1, 1939

By Thomas Henderson

Strategically, intermaterials competition is important to companies involved in raw materials extraction, smelting, refining or other processing, semi-finished and finished goods fabrication, and end-use consumers. In other words, just about every company involved in manufacturing. The array of materials covered by other papers in this session is testimony to the subject's breadth. In trying to provide an over- view of such a broad topic, the purpose has been to relate some of the common characteristics of the markets, economics and technology of raw materials, and to describe some of the ways in which Intermaterials competition impacts the markets, economics and technology of raw materials.

Jan 1, 1982

By Frank Woolley, John F. Elliott

The high-temperatzrre solulion calorimeter has been modified and an extension dynamic analysis of i/s transient behavior has provided an improved basis for interpretation of the experimental data. The partial molar heats of mixing, Him, of aluminum, copper, and silicon in the corresponding binary liquid-iron alloys haz,e been measured directly at 1600°C. The resulls in kcal per g-alom are: THE high-temperature1,2 solution calorimeter was planned and designed to obtain directly the partial molar heat of solution of solute elements in a variety of liquid metals that might serve as solvents. Results obtained in the temperature range of 1000° to 1200°C were described in the two earlier publications. However, a long-term goal in the program has been to obtain measurements at 1600°C of the heats of solution of elements in liquid iron. Hence, subsequent to the work last reported,2 and based on the experience of the first two studies, the calorimeter system has been modified and operated at 1600°C to obtain measurements on the Fe-Si, Fe-A1, and Fe-Cu systems. This paper reports the results of the work. THE MODIFIED SYSTEM It was clear from the earlier results that the most precise measurements can be obtained with the high-temperature solution calorimeter when the solution time is less than approximately 30 sec. This can be achieved only when the bath temperature, Tb, is above the melting (liquidus) temperature of the solute material being dropped. A second point of importance is that the refractories used earlier liberated considerable water so that, although the oxygen potential of the atmosphere in the system was satisfactory for work with nickel and copper, it was too high for work with iron, silicon, and so forth. Third, for work at 1550" to 1650°C, an improved temperature-sensing unit is needed. The modified calorimeter in its latest form is shown schematically in Fig. 1. The refractory brick and the Kanthal-Super heating elements which formed the furnace for heating the solute addition have been replaced by a smaller resistance-heated tube furnace which is insulated with bubble-grain alumina. A more elaborate gas purification system was installed to maintain a satisfactorily low oxygen potential in the atmosphere of the calorinieter chamber and of the solute addition furnace. As was the case in the earlier work, at the beginning of a run the calorimeter chamber is evacuated and then filled with purified argon so that a total pressure slightly in excess of 1 atm is maintained during the run. When either of the furnaces within the enclosure attains a red heat, 6 pct H2 is added to the atmosphere. Throughout the run, gas from the calorimeter is circulated through a de-oxidation system at the rate of approximately 5 standard liters per min; the changeover time for the gas within the calorimeter enclosure is approximately 1 hr. The essential stages of the deoxidation system are: a copper furnace at 600°C, a drying tube of anhydrite, and a titanium furnace at 800°C. The oxygen potential of the gas stream in the purification system can be monitored at several points by a stabilized zirconia oxygen cell. This cell is invaluable in detecting the development of leaks in the container vessel and the gas circulating system. The installation of the solute-addition furnace allowed either of two types of temperature-sensing units to be mounted on a water-cooled probe. The first was a small optical device using a silicon diode that was

Jan 1, 1968

By James H. TaVoularis, William J. Wiehagen

This paper discusses the application of commercially available microcomputers for assisting trainers in the delivery and evaluation of new hire, refresher, and occupational training. A brief overview of mine training highlights major issues facing the mine training community. The use of computer-based training has a potential for reducing costs and improving the transfer of conventional training to the job.

Jan 1, 1983

By M. R. Geer, H. F. Yancey

John Grifien (Pittsburgh)—I wish to congratulate the authors on this paper, which, I am sure, will promote a clearer conception of the various criteria which have been advanced as measures of coal-cleaning performance. Their clear descriptions of the more important criteria that have been advanced and explanations of their meaning and application have been sorely needed and, if widely promulgated and followed, should remove much of the confusion which, heretofore, has been evident. I suggest that much of this confusion has also arisen from lack of common agreement as to the exact meaning of the words and terms used in our discussions of this subject. I am sure I have been equally guilty in this respect and suggest it is time that we select a committee to study vocabulary as well as criteria, so that common understanding and agreement can be had. Thus, I suggest, that the title of this paper should have used the words "coal-cleaning performance" and not "coal-washery performance." My understanding is that "washery" refers to a complete plant and it has additional functions to that of coal-cleaning, such as control of size, moisture, etc. Also, I do not agree with the authors that the formulas given in their paper, the Fraser and Yancey formulas, are necessarily the only ones that can be called efficiency formulas. When a dictionary is consulted one realizes that the words "efficiency" and "efficient" have not precise meanings but only very general ones and must be further defined for each particular use, if they are to have any precise meaning. To my mind, the first of the Fraser and Yancey formulas given should never be used as expressing "efficiency" alone but only by the fuller phrase which truly expresses the kind of efficiency it states, which is "coal recovery efficiency." For complete clarity we then need only define "coal" for the particular circumstances under which the formula is used. This would recognize the fact, fully explained by the authors, that the limits in coal qualities imposed by the market have imposed a factor—and a variable factor—which has influenced any answer obtained by this formula. Such a procedure makes it entirely practical and logical to use the values of any of the qualities of coal—ash, sulphur, phosphorus, specific gravity—or a combination of them to define what is coal for that particular use of the formula. Further, it is clear to me that the criteria to be used must be dictated by the objective in view and the use for which the data are obtained. The authors clearly point out that the Fraser and Yancey coal-recovery efficiency formula is influenced by three factors, "one is determined by the coal, one is dictated by market or use conditions, and one is an inherent characteristic of the cleaning unit itself." The formula which uses the inverse of misplaced material as a measure of efficiency is also a coal recovery efficiency formula and is influenced in like manner by these three factors. As the authors point out, both are useless in obtaining data which enables a direct and valid comparison of the performance of a given cleaning unit on two coals where the other two factors vary or of the performance of two or more cleaning units of different types unless the values of the other two factors are identical. Much of the confusion in our discussions of and literature on coal cleaning has come from failure to understand clearly this fact. Many such comparisons have been made from which erroneous conclusions have been drawn. I am, therefore, particularly glad to see the authors' emphasis on 4?-1e criteria which indicate the "sharpness of separation. I feel that these are the essential criteria

Jan 1, 1952

By N. B. Hoffman

Much of the drill steel produced in America is manufactured into hollow rods. After all forging has been finished the entire bar should be thoroughly annealed and heat treated before the point is hardened and tempered, at least where the drill has heavy duty, such as working in deep holes and hard rock. A large number of tools of like character and doing work requiring similar properties are now being made of alloy steels; for instance a steel with a carbon content of 0.75 to 0.85 per cent., manganese 0.25 to 0.35, phosphorus under 0.020, sulfur under 0.030 and vanadium 0.18 to 0.23 has proved very satisfactory for chisels. A steel of this type properly manufactured will be far superior to the ordinary carbon steels, while at the same time is one whose price range should not prove prohibitive. Much of the cutting efficiency of a drill depends on proper heat treatment. This vanadium steel has a far greater safe hardening and forging range than ordinary carbon steels, making it more nearly foolproof. It will harden nicely at 1400" I?. but will also harden well and satisfactorily at 1550" F. retaining its quality of hardness and toughness in either case.

Jan 1, 1922

By K. E. Gray, M. S. Seth

Equations of elasticity for transversely isotropic, axisymmetric, homogeneous, porous media exhibiting pore fluid pressure were formulated. Using an analogy between thermal and porous body stresses, it was shown that the solution for a transversely isotropic porous body may be obtained by incorporating body forces and the stresses due to a boundary load into the corresponding solution for the thermal stress problem. Equations of elasticity ad the thermal analogy method were used to determine transient horizontal, tangential, and vertical stresses and radial displacement in a semi-infinite cylindrical region when either a constant pressure or a constant rate of flow is maintained at the wellbore. The vertical and tangential displacements are zero from the conditions of the problem. A numerical analysis was made of the solutions obtained by using a digital computer to determine the relative influence of each system variable. Considering rock as a porous body with internal fluid pressure generally gives results significantly different than considering the rock to be nonporous; the directional character of rocks leads to significant differences as compared to results based upon the common assumption of isotropy. Stress gradients are high near the wellbore but die out away from the well. Radial stresses are compressive or neutral, whereas tangential stresses are tensile, neutral or compressive, depending upon the boundary conditions and physical properties of the system Vertical stresses are compressive for an unbounded system. For constant wellbore injection rate, the vertical stress is proportional to the rate of fluid injection and decreases with time, whereas the radial and tangential stresses increase with time. At a given location, the radial displacement generally is very dependent upon time. INTRODUCTION A realistic appraisal of the state of stress in subsurface rock formations would be of considerable interest and use to the petroleum industry. For example, knowing the state of stress in proximity to a wellbore would be of fundamental importance in designing a fracturing operation or, more important, of clearly understanding the conditions necessary to produce rock failure of desired dimensions and geometry. Understanding conditions necessary for rock failure at the wellbore would also be of utility in a preventive sense. For example, borehole stability is an important consideration for many rock formations, and knowledge of the stress state at and near the wellbore under conditions of substantial pressure gradient due to fluid flow would be of great value. When fluid flows through a porous body which is initially at some uniform stress level, the following forces generate stresses at any point in the body. (1) Forces due to nonuniform pressure distribution. With increasing pressure the elements of a body are compressed. Such compression cannot proceed freely in a continuum when the pressure is not uniform throughout, and thus, stresses due to flow of fluid are set up. (2) Pore fluid pressure. This gives rise to normal stresses whose value at any point is the product of areal porosity and the fluid pressure. Although the fluid exerts uniform pressure, the stresses it creates in an anisotropic body may not be the same in all directions since the areal porosity in an anisotropic porous body is a direction-dependent quantity. This consideration leads to the concept of directional porosity. (3) Body Forces. The flow of the moving fluid imparts a force to the rock matrix, the magnitude of which depends upon the pressure

Jan 1, 1969

By C. W. Spencer, Che-Yu Li, R. J. Knight

The usual reaction product for the monotectic reaction in the Bi-Se system consists of rod -shaped particles of selenium-rich liquid embedded in a single crystal matrix of At certain crystallo-graphic orientations of the matrix with respect to the reaction interface, more complex shapes are as -sumed by the liquid particles. This reaction has been examined at different growth rates and relative crystallographic orientations of the matrix with respect to the reaction interface. The uhriation in auerage particle size of the liquid phase with reaction rate has been determined at one crystallograpIzic orientation of the matrix. On a plot of logarithm of average particle diameter vs logarithm of reaction rate, data points for crystals grown at rates of 10- to 1.0 cm per hr fall near a curve of slobe -0.4. MONOTECTIC reactions in binary alloy systems are reactions of the eutectic class in which, as the system is cooled through the three-phase reaction isotherm, one liquid phase (L) decomposes to a solid phase and another liquid phase (LII).' A partial phase diagram, as taken from Hansen, showing the monotectic reaction in the Bi-Se system is given in Fig. 1. This diagram specifies the monotectic point at 618°C and 72 at. pct Se, and the monotectic reaction products as the intermetallic Bi Se3 and a liquid (LII) which contains approximately 98 at. pct Se. At a lower temperature, 217°C, this liquid may undergo a eutectic reaction in which the reaction products are BizSe3 and crystalline selenium. However, under conditions of rapid cooling, or quenching, from above 217°C, this selenium-rich liquid may be partially retained to room temperature in an amorphous or glassy state. The intermetallic BizSe3 is rhombohedral, of the space group jm;' however, for convenience it is usually represented by a hexagonal struture. A pronounced plane of cleavage, which lies perpendicular to the main axis of the rhombohedral unit cell or parallel to the basal plane of the hexagonal cell, exists in the structure. The usual monotectic reaction product consists of rod-like particles of LII embedded in a single-crystal matrix of BizSe3, although at certain crystallographic orientations of the matrix BizSe3 other more complex shapes are assumed by LII. A similar rod-like, or cylindrical, morphology for the carbide phase has been reported by westbrook7 in the a Cr-Crz3Ce eutectic structure. In this paper the morphological characteristics of the LII particles, as retained to room temperature, have been examined under conditions of varying growth rate and crystallographic orientation of the matrix BizSe3 with respect to the growth direction. EXPERIMENTAL PROCEDURE Alloy specimens corresponding to the monotectic composition, 72 at. pct, or 50 w pct Se, were prepared from 99.999 pct pure bismuth and selenium obtained from the American Smelting and Refining Co. The two elements, corresponding to a total specimen weight of approximately 30 g, were encapsulated in 13 mm OD Vycor tubing under a vacuum of 5 x X mm of Hg. The capsules were suspended on a fine platinum wire in a Bridgman-type furnace which exhibited a temperature gradient of

Jan 1, 1963

By D. P. Harris

The construction of a model to associate probability of occurrence of some measure of mineral wealth with the geology for each subdivision (cell) of the area is postulated. The questions (1) are the observations made in reconnaissance geology sufficiently related to mineral wealth to distinguish the cells of high value from those of low value, and (2) can this geology be related mathematically to mineral wealth and probability are answered in a case study of 97,000 square miles in Arizona and New Mexico. Measurements on mineral wealth and geologic variables of 243 cells were made. The geologic variables, valued in counts, percentages and lengths, quantified the theoretical concepts of the model. It was found that multiple discriminant analysis and classification analysis by Bayesian statistics and multivariate normal probability function constitute a two-phase probability model that associates probabilities with geologic information and mineral wealth. The cells of explored area must be allotted the value class indicated by their cumulative production and estimate reserves (given well-defined geological relationships in a mineral wealth area), and a discriminant and classification analysis must be performed in which probabilities of the cell belonging to each of the value groups are generated as a function of the geologic variables. This determines whether high value cells can be distinguished from low value cells by their geology as well as indicating those cells misclassified by the model into a higher value group than warranted by their known mineral wealth. The second use of the model is for extrapolation. It is demonstrated how geologic information on an explored area can be utilized in guiding exploration in an unexplored area. Exploration decisions can be predicted on explicitly determined measures of uncertainty (the probabilities). SUMMARY This study postulates the construction of a model to associate probability of occurrence of some measure of mineral wealth with the geology for each subdivision (cell) of the area. Such a model will distinguish those cells deserving further exploration on the basis of their geology: The basic postulates of this study are as follows: v = Q(R,S,F,A) P(V) = G(R,S,F,A,V), where V = a measure of mineral wealth P(V) = the probability of occurrence of V R = age and type of rock F = rock fracturing A =age of igneous activity and contact relationships. There are really two questions to be answered: 1) are the observations made in reconnaissance geology sufficiently related to mineral wealth to distinguish the cells of high value from those of low value, and 2) can this geology be related mathematically to mineral wealth and probability? In order to answer these questions, an area of about 97,000 square miles situated chiefly in Arizona and New Mexico was selected as a case study. This case study area was divided into 243 cells, each 20 miles square, and measurements were made on mineral wealth and geologic variables. Crude geologic variables valued in counts, percentages and lengths are constructed to reflect the above postulated variables. These measurement variables thus provide quantified values defining the theoretical concepts of the model. It was found that multiple discriminant analysis and classification analysis by Bayesian statistics and the multivariate normal probability function constitute a two-phase probability model that associates probabilities with geologic information and mineral wealth. Given the availability of an area that is well explored and upon which the relationships of mineral wealth to geology can be adequately defined, there are two ways in which the analysis outlined in this paper can prove beneficial. First, the cells of the explored area must be allotted the value class indicated by their cumulative production and estimated reserves, and a discriminant and classification analysis must be performed in which probabilities of the cell belonging to each of the value groups are generated as a function of the geologic variables. Such an analysis determines whether cells of high

Jan 1, 1967

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

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

Jan 1, 1967

By Paul A. Farrar, Harold Margolin

The Ti-Al-V system has been delineated from 50 to 100 wt pct Ti and front 600 to 1400°C by X-ray and ntetallographic techniques. Isothermal sections were delineated at 600, 700, 800, 900, 1000, 1100, 1200, 1300, and 1400°C. X-ray and metallographic evidence indicated the presence of two phases and in addition to those reported by previous investigators of the ternary system'-3. T HE titanium-rich region of the Ti-Al-V system as originally proposed by Rausch et UZ, and Jordan and Duwez3 was based on the binary Ti-Al diagram as proposed by Bumps et uZ. As additional phases have been shown to exist5-' in the region previously thought to be all a, the titanium-rich region of the Ti-Al-V system is open to serious question. Therefore, in order to determine the effect of the indicated changes in the Binary Ti-Al system on the ternary Ti-Al-V system the present investigation was initiated. EXPERIMENTAL PROCEDURE The experimental procedures in this investigation are the standard techniques used in the study of titanium-alloy systems which have been described in detail previously.1° Consequently, only details pertinent to the present work are discussed. The alloys employed for the delineation of this system were prepared from titanium obtained from the Bureau of Mines (73 BHN) 99.99 pct Al and 99.7 pct V, by the consumable arc-melting technique. The titanium was premelted into buttons of 25 g because spattering during the first melt caused uncontrollable loss of alloying constituents. Typical analyses of the material used have been reported elsewhere Attempts to hot roll the alloys prior to heat treatment were made on the first series of alloys prepared, but as the area of extensive hot workability was fairly restricted, see Table I, this practice was discontinued. The times of isothermal annealing are given in Table 11. A preliminary heat treatment of 96 hr at 1000°C was used prior to heat treatment at lower temperatures. The standard techniques used for polishing specimens involved belt-grinding, grinding on emery paper, polishing electrolytically, and etching with Remington "A" etch. In the low-alloy region it was found that the oRo etch, developed by Ence and Margolin,o gave better contrast between phases. Where this procedure did not differentiate between phases clearly, the method of stain etching developed by Ence and was used. It has been previously reported15 that filing of titanium alloys is not desirable for producing X-ray powders, since the powder becomes contaminated by the file material. Therefore, for those samples not brittle enough to Crush, the procedure described below was used. The specimen was wrapped in molybdenum sheet and placed in a vacuum system and evacuated to 0.03 p. Palladium-pur if ied hydrogen was then admitted to the system until a partial pressure of 5o Cm of H was reached. The specimen was heated in the system to the temperature of original heat treatment and was slowly cooled over a period of 3 hr to 400°C while the partial pressure of hydrogen was maintained. The sample was removed from the system after it had completely cooled and was crushed to a powder of 230 to 270 mesh. The powder, wrapped in molybdenum sheet, was dehydrogenated at the temperature of original heat treatment until a partial pressure of 0.03 p was maintained in the system for at least 1/2 hr. The tube containing the powders was quenched in iced brine while still connected to the vacuum system. After dehydrogenation, the powder, except for the

Jan 1, 1962

By E. L. Dougherty

A phenomenological theory for a one-dimensional unstable miscible displacement similar in type to the Buckley-Leverett model but including the effects of mixing is proposed An equation giving the fractional flow of pure solvent through an oil phase containing dissolved solvent is derived including the effects of gravity and nonbomogeneities. Numer ical integration of the resulting pair of nonlinear hyperbolic partial differential equations gives volumes of dissolved and undissolved solvent as functions of space and time. Parameters in the model which characterize mixing due to dispersion were determined from experimental data; the parmeters correlated with viscosity ratio. The results indicate that the rate of dispersive mixing is proportional to volumetric flow rate Incorporated into our equations were Koval's heterogeneity factor H. which characterizes mixing due to channeling. This allowed satisfactory predictions of displacement behavior observed in horizontal floods in nonhomo-geneous cores; but in most cases the values of H which we used were considerably larger than those used by Koval. INTRODUCTION It has been shownl,2,3 that for favorable viscosity ratios, the diffusion equation with convection satisfactorily describes the behavior of a miscible displacement in a porous media. Numerous attempts to develop a satisfactory mathematical description for the case of unfavorable viscosity ratio have been reported, but they have met with only partial success. These attempts, which are reviewed by oval! have taken two approaches: (1) development of a flow model akin to the Buckley-Leverett approach for immiscible displacement neglecting the effects of mixing5,6 and (2) simultaneous application of Darcy's Law for flow and the diffusion equation with a convection term for mass transfer.7,8,9 We set out to construct a mathematical model of Type 1 including, though, the effects of mixing. Based upon a set of hypotheses, we derived a pair of nonlinear hyperbolic partial differential equations which describe in a one-dimensional system the combined effects of flow and mixing. To work with these equations it was necessary to develop a fractional flow formula. The combined system was integrated numerically using the method of characteristics. In our equations the mixing process is characterized by four parameters. Three of these, labeled ß, P1 and P2, account for dispersive type mixing. The fourth is the heterogeneity factor H, proposed by Koval to account for channeling due to nonhomogeneities in the porous media. Values of the parameters which would cause agreement between theory and experiment were determined by trial-and-error for miscible floods conducted in horizontal cores of both homogeneous and heterogeneous materials. Calculations were also performed for vertical floods in homogeneous cores. The purpose of this paper is (1) to present the details of the mathematical analysis, (2) to present the results of the calculations, (3) to consider what light the results shed on the mixing process in an unstable miscible displacement, and (4) to provide a firmer foundation for the correlation technique developed by Koval for predicting the behavior of unstable miscible floods. STATEMENT OF PROBLEM The problem is to construct a mathematical model which describes in one dimension the observed behavior in an unstable miscible displacement. The approach is phenomenological in that the equations are based on assumptions which violate certain physical precepts known to apply to the displacement phenomenon. However, the assumptions do allow us to account mathematically for the more essential phenomena. The work is of value if we can quantitatively predict experimental results which heretofore could not be predicted, even though the synthetic nature of the model belies complete explanation of the observations. We assume that the system is comprised of two contiguous flowing phases. One phase, which we call free solvent, has the properties of pure solvent. We call the fraction of the pore volume occupied by this phase the solvent saturation, designated s. The oil phase consists of oil containing dissolved

By R. E. Ogilvie, A. R. Kaufmann, S. H. Gelles

The solid-solubility limits of iron in beryllium were determined between 850o and 1200oC by analysis of differential type multiphase diffusion couples, using an X-ray absorption technique. The maximum value of the solubility limit was found to be 0.92 ± 0.02 at. pct (5.46 wt pet) at the eutectic temperature 1225°C. The solubilities of nickel and beryllium were determined between 900°and 1200°C by the same technique and the maximum solubility was found to be 4.93 + 0.01 at. pct (25.2 wt pet) at the eutectoid temperature, 1065°C. A previously unreported high-temperature phase which decomposes eutectoidally at 1065 °C was found to exist in the beryllium-nickel system at a composition of approximately 8 at. pct Ni (36 wt pet) by diffision-couple analysis. The presence of this phase was confirmed by thermal analysis and metallo-graphic analysis of the structure resulting from the eutectoid decomposition. G. V. Raynor1 has treated the solid solubilities of some of the elements in beryllium on the basis of the "Hume-Rothery" rules2 which have been modified to include ionic size and ionic distortion effects. It was predicted that the solubility of iron and nickel in beryllium should be slightly less than that of copper. The lowering of the solubility, according to Raynor, is due to a more unfavorable relative valency effect and an ionic size effect. Kaufmann and corzine3 have compiled data on the solubilities of elements in beryllium and have discussed them in the light of the Raynor paper. These authors feel that, because the elements having the greatest solubility in beryllium systematically fall in the Group VIII and IB Columns of the periodic table, the electronic structure greatly influences the maximum solid solubility of elements in beryllium. The solubility of iron in beryllium was determined by Teitel and cohen4 as part of the study of the beryllium-iron phase diagram. The determination was carried out by X-ray and thermal analysis and according to the phase diagram presented, the maximum solubility of iron in beryllium is 0.41 at. pct (2.5 wt pct) at 1225oC. However, it is estimated that the uncertainty in the position of the a-beryllium primary solid-solution boundary is about 0.5 at. pct (3wtpct). Losana and Goria3 in studying the beryllium-nickel phase diagram, determined the solid solubility of nickel in beryllium by thermal analysis. They found the maximum solubility to be between 1.65 and 2.65 at. pct (10 to 15 wt pct) at 1240°C. This value decreased rapidly with decreasing temperature. In determining approximate ranges of solubilities for different elements in beryllium, Kaufmann, et al,8 reported a value of between 1.3 and 1.7 at. pct (7.9 to 10.1 wt pct) for the solubility of nickel in beryllium. The value was obtained by metallographic examination of quenched alloys and lattice-parameter measurements. However, the authors also noted a single-phase structure for a 1.7 at. pct Ni alloy (10 wt pct) on cooling from the liquid. This would indicate a higher solubility range than was reported. ~isch,' in his X-ray studies of beryllium-copper, beryllium-nickel, and beryllium-iron intermetallic compounds, reports the disappearance of a second phase (Ni,Be2) in the beryllium primary solid solution at approximately 4 at. pct (20 wt pct). THEORY The analysis of concentration gradients in diffusion couples has proven to be a useful tool in determining phase equilibria.8-14 In this particular study the diffusion couples were chosen to straddle the expected composition range of the phase boundary, then heat treated at a given temperature and the concentration gradient evaluated. The composition of the phase boundary for a given temperature appears at a point of discontinuity of the composition gradient. Examples of typical phase diagrams and the concentration gradients which should be found in such systems are shown in Fig. 1. In the present work, gradients of the form of Fig. l(c) were obtained in diffusion couples made of pure beryllium and two-phase alloys of beryllium with either iron or nickel. The composition at the point where the gradient becomes discontinuous, Cs, corresponds to the solubility limit of either iron or nickel in beryllium. The analysis of the concentration gradients was carried out by an X-ra absorption method developed and applied by Ogilvie and later used by Moll13 and Hilliard.l4 It depends on the fact that the absorption of X-rays by matter is determined by the concentration and type of the various atomic species present. The relationship for the intensity, I, of a monochro-

Jan 1, 1960