Search Documents

By S. S. Arentz

PLANNED health and safety programs have become an essential part of American industry because such programs lead to increased operating efficiency, improved labor relations, better public relations, and to substantial savings in compensation insurance. Those of you who have had the unpleasant duty of informing the wife or widow of one of your men of his serious injury or death while on the job, know that all the benefits of a successful safety program do not show on the balance S. S. ARENTZ, Member AIME, is General Superintendent, Nevada Operations, Combined Metals Reduction Co., Pioche, Nevada. AIME San Francisco Meeting, February 1949. TP 2741 A. Discussion of this paper (2 copies) may be sent to Transactions AIME before March 31, 1950. Manuscript received Jan. 6, 1949. sheet. These programs are of particular importance to the mining ,industry because mining's reputation as an unusually hazardous industry and the commonly isolated location of mining operations tend to focus attention on these problems. Description of Operations: Before proceeding with a discussion of our health and safety programs at Pioche, it may be proper to give a brief description of Pioche and of our operations there. Pioche is one of the early Nevada mining camps. It was founded shortly after the discovery of high grade silver ore in 1863 and mining has continued with more or less regularity to the present day. In an era of lawlessness, Pioche was notorious. The story persists that 75 men died with their boots on before one died a natural death, and old payroll records show that nearly as many gunmen were employed to stand off claim jumpers as there were miners working the mine. That was probably as close to a safety program as the times permitted. Pioche is situated in southeastern Nevada on the main highway between Ely and Las Vegas. The camp is on the flank of "Treasure Hill," near the original silver discovery, at an elevation of about 6000 ft. The present day population of about 2000 is primarily dependent upon the mines of the area, although Pioche also serves as the county seat of Lincoln Couqty and as the center of the surrounding livestock industry. The camp is served by a branch of the Union Pacific Railroad and receives power from the generators at Hoover Dam. The Pioche operations of the Combined Metals Reduction Co. were started in 1923 when the first complex lead-zinc ore was shipped to the company's mill at Bauer, Utah. The modern mill at Pioche was completed in 1941. The operations are medium sized in the nonferrous field, employing an average of 350 men in the mine, mill, and related works. The complex lead-zinc ore is mined from replacement deposits in a comparatively flat, extensively faulted, limestone horizon. Mining methods vary from stull-supported open stopes to filled square-set stopes. The thin bedded limestone and shale overlying the ore is allowed to cave as areas are mined out and caving frequently follows closely upon ore extraction. The relatively heavy ground and the numerous faults add to the problems of safe mining. The mine is well mechanized and the mill and surface plant are modern and well equipped. Labor is organized in a C.I.O. union and labor-management relations have been unusually harmonious. During most of the period since 1923 a competent supervisory staff worked to reduce safety hazards but the primary responsibility for safety rested on the individual workman. Accidents happened and all too frequently they were regarded by all concerned as unavoidable. In October 1939, the late Robert L. Dean became superintendent at Pioche. Most of his previous experience had been in the fields of iron and coal mining and from that experience he brought the concept that no accident is unavoidable. Many of the features of our present health and safety programs were initiated by Mr. Dean during his term as superintendent. Health Program: Our health program centers in Dr. Q. E. Fortier and his new, well-equipped, and well-staffed, modern hospital in Pioche. The program starts with a thorough pre-employment physical examination and is followed by yearly re-examinations at the expense of the company. The Pioche Mutual Benefit Association, to which all Pioche mine operators and employees belong, pays benefits covering hospitalization and surgery expense incurred by employee members and their families. The Association is governed by a board of directors elected by its members. The mine operators of the district donated the original capital and pay the monthly dues of the employee members. The employees pay the dues covering members of their families. Though not strictly a part of the

Jan 1, 1951

By Herman Dykstra, R. L. Parsons

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 John Johnston

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

Jan 1, 1942

By T. A. Pollard, C. R. Sandberg, E. B. Elfrink

This paper presents an evaluation of compressibility factor data and a discussion of their application to the estimation of gas reserves. A correlation is presented which provides compressibility factors for use in both two-phase and single-phase hydrocarbon systems. Accuracies comparable to those obtained previously for single-phase systems only can be expected. A simple means of predicting the presence or absence of a liquid phase in a condensate system of known composition is illustrated. The correlation is based on 1,030 compressibility determinations from 21 hydrocarbon samples taken from eight oil fields. Of the data used, 75 per cent were from California, 15 per cent were from the Mid-Continent area, and 10 per cent were from South America. The average numerical deviation of the experimental data from this compressibility chart is 1.22 per cent. Charts and tables are included and discussed which illustrate the errors involved through the misuse or nonuse of compressibility factors in estimates of gas-in-place. INTRODUCTION The compressibility factor is a coefficient which expresses the deviation of a gas of given composition from the Perfect Gas Laws. A factor such as this is necessary in any calculation involving volumes of gaseous mixtures, and finds extensive application in the estimation of natural gas and condensate reserves. It is used in the decline curve, the volumetric, and the material balance methods of gas reserves estimation. The behavior of gases at high pressures has been investigated extensively in recent years and considerable data have been assembled on the compressibility characteristics of a number of gases. Notable among the investigators of high pressure gas relationships are Kvalnes and Caddya, Kay" Sage and Laceys, and Brown, Standing" and Katz4. Other contributors are Smith and Watson, Roland and Kaveler', and Stevens and Vance". The data on compressibility factors assembled by these investigators have been presented and correlated in numerous ways; particularly noteworthy is the method first presented by Kay5 which relates compressibility factor "Z" to pseudo-reduced temperature and Pseudo-reduced pressure. Other investigators have used this empirical relationship successfully and numerous charts covering a large number of gases have been constructed on this basis. It is the purpose of this paper to evaluate compressibility factor data and to illustrate the importance of their proper use in the estimation of gas-in-place for both gas and condensate reservoirs. EVALUATION OF COMPRESSIBILITY FACTOR DATA The compressibility chart presented by Standing, Katz, Brown and Holcomb in 1941 is probably the most recent and reliable chart for natural gases. The reliability of the Katz chart has been investigated by the writers for a fairly large number of gases (single-phase) in the pressure range from 500 to 10,000 psia and temperature range from 80 to 400°F. A summary of this comparison is presented in Table I. The numerical average deviation of the chart compressibility factors from the experimental compressibility factors is 1.05 per cent, and the algebraic deviation is -0.28 per cent. The Katz correlation appears to be quite accurate for single-phase hydrocarbon systems. It was based on single-phase systems and does not necessarily provide a means of accurately predicting the behavior of two-phase systems which might be encountered in condensate work. Data published recently by Sage and Oldsl' have enlarged the knowledge of phase behavior in condensate systems by covering extensively the Compressibilities of systems in the two-phase as well as single-phase regions for a wide range of temperatures, pressures and gas-oil ratios. An evaluation of these data indicates that conditions can arise, especially in condensate systems having fairly low gas-oil ratios, in which two phases occur and the usual gaseous compressibility factors do not accurately apply. The measured PVT data of Sage and Olds covering two-phase systems have been correlated with additional such data from the Mid-Continent area and South America; these are presented in Fig. 1. This chart is the result of 1,030 compressibility* determinations from 21 hydrocarbon samples taken from eight oil and condensate fields. Of the data used, 74.75 per cent were from California, 15.45 per cent were from the Mid-Continent area, and 9.80 per cent were from South America. The chart construction is based on the original

Jan 1, 1949

By W. F. Kieschnick, C. C. Miller, E. R. Brownscombe

The possibility has been mentioned that large pressure gradients in a solution gas driven field caused by high production rates might lead to a reduction in the ultimate recovery obtainable compared to that which would be obtained by a very slow rate of production. In the present study the reservoir conditions accompanying a high rate of production and the corresponding ultimate recovery were compared with those obtained if the reservoir were produced at some marginal rate throughout its entire life. The method of calculation involved the application of fluid flow - material balance analysis to a series of successive steady state conditions in the reservoir. In the case studied very little difference in recovery was obtained at the abandonment pressure. An analysis of the basic factors involved indicates that the sane results would hold for considerable variation in properties of reservoir fluids, permeability of the formation or well spacing. Even if the reservoir has a water drive it appears that no harm will be done by- open flow production if the rate is cut back before any appreciable free gas is produced. However, any condition which leads to disproportionate withdrawal rates and which causes large pressure differences over large areas might result in substantial loss in ultimate recovery. INTRODUCTION When the pressure declines in a reservoir, gas is evolved from solution causing a shrinkage in the reservoir oil so that a stock tank barrel of oil with its residual dissolved gas occupies a smaller volume than originally. For a given reservoir oil saturation, therefore. stock tank oil content will be higher at the lower pressure. At the same time the evolved gas which remains in the reservoir reduces the reservoir oil saturation, and this would be expected perhaps to cause a reduction in stock tank oil content. However, around the well the shrinkage occurs so rapidly that stock tank oil may actually accumulate in this region in spite of the simultaneous reduction in reservoir oil saturation. If this accumulation should continue, spreading into a larger zone, the loss might be appreciable and because larger pressure gradients accompany higher rates the accumulation might possibly be large if wells are produced at too high a rate. It is the purpose of this paper to determine the effect of rate of flow on ultimate recovery by solution gas drive. The problem under consideration involves the unsteady state flow of oil in a closed radial system from which oil is produced by the evolution of gas as a finite drawdown is estalllished and the reservoir pressure declines. A somewhat similar unsteady state production problem dealing, however, with a linear system has been considered by Muskat and Meres, but these authors do not discuss the effects which rate may have on recovery. Other methods have been reported for estimaling recovery as a function of pressure decline when producing at an infinitesimal rate, and Muskat and Taylor' have shown the effect of fluid and reservoir characteristics on recovery under these conditions. Assuming steady state to exist, Babson has given a procedure for computing the decline in well productivity as the field is depleted, and Moyer has combined these methods to estimate abandonment conditions and the resulting ultimate recovery. These latter approaches to the problem either do not consider the effect of rate or assume steady state can be established. In detail this paper presents the behavior of the reservoir producing under unsteady state conditions and considers the effect of production rate during this unsteady state flow on the recovery attained. The analysis also shows that there will exist at abandonment a type of steady state, and on this basis a method has been clevelopetl for rapid calculation of abandonment conditions.

Jan 1, 1949

By P. B. Baxendell, R. Thomas

Work on the calculation of vertical two-phase flow gradients by Cia. Shell de Venezuela has been based mainly on the "energy-loss" method proposed by Poett-mann and Carpenter in 1952. The "energy-loss-factor" correlation proposed by Poettmann and Carpenter was based on relatively low-rate flow data. This correlation proved inapplicable to high-rate flow conditions. In an attempt to establish a satisfactory correlation for high rates, a series of experiments was carried out at rates up to 5,000 BID in Cia. Shell de Venezuela's La Paz field in Venezuela, using tubing strings fitted with electronic surface-recording pressure elements. As a result of these experiments a correlation between energy-loss factor and mass flow rate was established which is believed to be applicable to a wide range of conduit sizes and crude types at high flow rates (e.g., above 900 BID for 27/8-in. OD tubing). It is anticipated that the resulting gradient calculations will have an accuracy of the order of % 5 per cent. At lower flow rates the energy-loss factor cannot be considered as constant for any mass rate of flow, but varies with the free gas in place and the mixture velocity. No satisfactory correlating parameter was obtained. As a practical compromise for low flow rates, a modification of the curve proposed by Poettmann and Carpenter was used. In practice this was found to give gradient accuracies of approxirnately ± 10 per cent clown to flow rates as low as 300 B/D in 27/8-in. tubing. INTRODUCTION Production operations in Cia. Shell de Venezuela's light- and medium-crude fields are principally concerned with high-rate flowing or gas-lift wells. Under these conditions the analysis of well performance, the selection of production strings and the design of gas-lift installations are vitally dependent on an accurate knowledge of the pressure gradients involved in vertical two-phase flow. Initially, attempts were made to establish the gradients empirically as done by Gilbert,' but the results were not reliable due to scarcity of data over a full range of rates and gas-oil ratios. Several methods of calculation based on energy-balance considerations were attempted, but the computations were cumbersome and the results cliscouraging. In 1952 a paper was published by Poettmann and Carpenter' which proposed a new approach. Their method was also based on an energy-balance equation. but it was original in that no attempt was made to evaluate the various components making up the total energy losses. Instead, they proposed a form of analysis which assumed that all the significant energy losses for mutiphase flow could be correlated in a form similar to that of the Fanning equation for frictional 1osses in single-phase flow. They then derived an empirical relationship linking measurable field data with a factor which, when applied to the standard form of the Fanning equation, would enable the energy losses to be determined. The basic method was applied in Venezuela to the problem of annular flow gradients in the La Paz and Mara fields" This involved establishing a new energy-loss-factor correlation to cover high flow rates and, also, some adaptation of the method to permit mechanized calculation using punch-card machines. The final result was 1 set of gradient curves for La Paz and Mara conditions which proved to be surprisingly accurate. With the encouraging results of the annular flow calculations, several attempts were made to obtain a corresponding set of curves for tubing flow. Here, unfortunately, little progress could be made. The original correlation of Poettmann and Carpenter was based on rather 1imited data derived from low-rate observations in 23/8- and 27/8-in. OD tubing. It did not cover the higher range of production rates, and extrapolation proved unsuccessful. A new correlation covering high flow rates was required, but this proved to be extremely difficult to establish since tubing flow pressure measurements at high rates did not exist—due to the difficulty of running pressure bombs against high-velocity flow. The necessity for reliable tubing flow data increased with the development of the new concessions in Lake Maracaibo, where high-rate tubing flow from depths of 10,500 ft became routine. Thus. it was decided to set up a full-scale test to establish a reliable energy-loss factor for tubing flow conditions. A La. Paz field light-oil producer with a potential of approximately 12,000 B/D on annular flow was chosen. To obtain full pressure gradients, a special tubing string was installed which was equipped with electronic surface-recording pressure measuring devices,

By F. C. Bond, J. T. Wang

H. J. Kamack (E. I. du Pont de Nemours & Co., Inc., Wilmington, Del.)—Rittinger's law usually is stated to the following effect: "The work (or energy) consumed in particle size reduction is proportional to the new surface area produced." The law has been stated substantially in this way by Taggart20 Berry21 Dalla Valle22 Coghill and DeVaney,23 Richards and Locke," Gross,25 and many others, and, according to Gaudin,26 was originally expressed by Rittinger in the same form. Consequently there can be little doubt that this is the understanding of the law among most workers in the field of particle size reduction. Bond and Wang, however, express the law in the form "the useful work accomplished . . . ." (italics theirs). The distinction is critical, for in the form used by Bond and Wang the law becomes, as they themselves remark "merely a more or less arbitrary definition of useful work," while in its usual sense the law expresses a physical hypothesis which has been verified experimentally within certain limitations. Considering the way the law has been used, it might be stated more explicitly as follows: "In a given machine operating under a given set of constant operating conditions, the work consumed in the particle size reduction of a given material is proportional to the new surface area produced." Or, as Coghill and deVaney have said,27 "the (Rittinger) law holds only when the tests being compared are made under analogous conditions." There are occasions when the law transcends these limitations; for example, the surface area produced per unit energy consumption for a given material in a ball mill does not vary much over a fairly wide range of operating conditions. But by and large, the surface area production per unit energy consumption will vary with the operating conditions, the type of machine, and the material. The essence of Rittinger's law is that the surface area production per unit energy consumption is independent of the particle size, and this has been verified experimentally by numerous workers for numerous materials, within certain limits. An important limitation is that when one grinds to very small particle sizes, agglomerative forces may tend to interfere with size reduction so that the surface area increases less rapidly than Rittinger's law would predict. Within such limitations, Rittinger's law can be regarded as empirically established. The law has, however, certain theoretical implications, and it seems to be chiefly against these that Bond and Wang direct their criticism. Solids are believed to possess surface energy which is proportional to surface area. Thus, Rittinger's law implies a proportionality between surface energy produced and mechanical energy expended (for a particular material in a particular machine). It does not imply that all or most of the mechanical energy is transformed into surface energy; in fact it is known that most of the mechanical energy is transformed into heat. Bond and Wang assert that most of this heat arises from the damping of elastic vibrations of stressed particles. This may possibly be true for crushing (with which they are chiefly concerned), although in grinding it is probable that much of the heat arises from friction between particles. However, the fact that the surface energy is small compared to the heat energy does not invalidate Rittinger's law, which implies merely that they are proportional. The authors also criticize Rittinger's law on the grounds that "this theory cannot be justified mathematically, since work is the product of force times distance, and the distance factor is ignored," and "energy input must be the product of force times distance, and the Rittinger theory completely ignores large variations in the distance (strain dimension or deformation) throughout which a force must act to produce breakage of different materials." However, the quantities force and distance as such are irrelevant to Rittinger's law, which considers energy input as an integral quantity. The fact that different materials require different forces and strains (hence, different amounts of energy) to break them is incontrovertible but again is irrelevant to Rittinger's law, which, as mentioned before, applies -only to a single material. Even in theory, it would not be expected that the surface area produced per unit energy consumption would be the same for different materials, because the surface energy per unit area is presumably specific to each material. Bond and Wang advance another argument of a

Jan 1, 1952

By D. A. Dukelow, K. Li, G. C. Smith

Decarburization in the Fe-C system by oxygen top blowing has been studied in laboratory -scale experiments. It is shown that equilibrium models fail to explain or predict either the course of refining or endpoint conditions, giving results which either are incompatible with the chemistry of the system or do not satisfy material balance requirements. Also the path of decarburization was found to vary even for heats made under apparently identica1 conditions. A promising approach to analyzing the decarburization results is to relate oxygen efficiency fm carbon removal to bath carbon content. This relationship for Fe-C heats shows the same range of oxygen efficiencies as is obtained in pilot-plant and commercial heats using hot metal-scrap charges. This implies that oxygen transfer is primarily controlled by the decarburization reaction itself, independent of other refining reactions. Therefore, it should be possible to study separately decarburization and slag-metal reactions. DECARBURIZATION is probably the most important reaction in steelmaking. Not only is it a main reaction in the refining of iron to steel but it also provides the stirring action in the bath necessary for the diffusion processes to proceed at reasonable rates so as to make a steelmaking process practical. Kinetics of decarburization in the open-hearth process has been a subject of investigation for many years.&apos;-B It is generally accepted that at steelmaking temperatures the rate of homogeneous C-0 reaction is extremely high and cannot constitute a rate-controlling step. Diffusion of oxygen through a boundary film in the metal phase has been suggested by arken&apos; as rate-determining. Recently, Larsen and sordah16 concluded from experiments in a laboratory furnace that, with oxygen supplied from air or combustion gases, the rate of "steady-state" carbon boil is controlled essentially by a diffusion process of O2, Co2, or H2O through a film of nitrogen above the slag surface. Displacing this diffusion film by a stream of nearly pure oxygen produced a ten-fold increase in the rate of carbon boil with the rates of slag-metal oxygen transfer, bubble nucle-ation, and other steps all apparently able to keep pace. In the top-blown basic oxygen process, however, the transport of oxygen takes a more direct route. and the state of bath agitation is much more turbulent than in the open-hearth process. In addition, direct contact of the gas with the metal phase provides opportunity for direct oxidation of carbon. It is likely that the rate-limiting factor for the decarburization reaction will be different. However, only a few descriptive discussions of the subject have been reported in the literature.10-l2 Studies of the decarburization kinetics based on plant or pilot-plant data are necessarily complicated and are influenced by other refining reactions which occur simultaneously. In order to investigate the mechanism of decarburization, experiments have been conducted in which carbon-saturated iron melts were top-blown with pure oxygen over a range of conditions. It is hoped that this study will form a foundation on which a more basic understanding of this important reaction may be built. EXPERIMENTS One group of blowing experiments was made in a standard 200-lb induction furnace and another group in a 500-lb induction furnace. The furnaces were modified to the general shape of a basic oxygen furnace by adding a rammed refractory cone section to the regular crucible body. Crucible and cone were of high MgO (95 pct) material. A water-cooled lance, 1/2 in. in diam and threaded at one end to take a nozzle, was used for blowing oxygen. The lance with its water and oxygen lines was supported on a cantilever arrangement so that it could be moved up, down, or sideways. Oxygen of 99.5 pct purity was supplied from a cylinder and metered through a rotameter equipped with pressure and temperature gages. Another pressure gage was located at the top of the lance. A schematic diagram of the assembly is shown in Fig. 1. Before each experiment, a weighed amount of ingot iron, containing 0.02 pct C, < 0.01 pct Si, 0.10 pct Mn, 0.019 pct P, and 0.015 pct S, was charged in the furnace and melted down by induction heating. Graphite was then added to the molten charge until it became saturated. When the temperature of the charge reached the desired level, the lance was lowered to a predetermined height above the bath

Jan 1, 1964

By R. Taggart, D. H. Polonis, W. C. Gallaugher

In the Ti-Cu system, the a&apos; phase can be produced over a wide range of alloy composition witJwut the retention of measurable amounts of the ß or ? phases. This paper reports on the decomposition of this hexagonal martensite phase front the standpoint of solute corzcentratiorz, tempering temperature, and coherency strain condition. The mechanism of precipitation during tempering comprises localized precipitation which is followed by discontinuous precipitation if sufficient coherency strains are present. The localized precipitation process in hypoeutectoid alloys is described by the generalized Johnson-hiehl equation and an activation energy of 46,000 cal per mole. In hypereutectoid alloys the corresponcEing activation energy is 51,000 cal per mole. The rate-controlling process has been proposed as the difision of copper along a&apos; platelet interfnces. SEVERAL investigations have been reported concerning the kinetics of decomposition of the mar-tensitic a&apos; phase in titanium binary alloys.&apos; In the Ti-Mo system,&apos; two modes of a&apos; decomposition have been observed in the presence of /3 phase; one of these involves the precipitation of fine a from a&apos;, and the other involves diffusion across the a&apos;-ß interface. In studies of the Ti-Cr system, Rostoker2 did not detect the formation of TiCr2 during tempering. In Ti-Ni alloys,3 the intermediate phase Ti2Ni has been found to precipitate directly from a&apos; along the interfaces between plates. With the exception of the study of Ti-Ni alloys the previous investigations of tempering phenomena in substitutional martensites are mainly qualitative and do not present a detailed description of the precipitation processes. The following limitations restrict the correlation of the microstruc-tural processes and reaction kinetics during tempering in most binary-alloy systems. 1) A mixture of at least two phases characterizes the constitution of quenched alloys. 2) Difficulties have been encountered in obtaining uniform structures throughout quenched samples. 3) The reaction products, and in particular their morphology. have not been clearly resolved. 4) The martensitic a&apos; phase forms over only a limited composition range for most titanium-base binary alloys. Gomez and polonis4 showed that the Ti-Cu system provides an excellent basis for investigating the tempering of a&apos; over a range of composition without the complication of retained P or the transition w phase. In the present work the effects of solute concentration, tempering temperature, and coherency strain conditions are considered with reference to the over-all precipitation process in quenched alloys of the Ti-Cu system. Both microstructural observations and kinetic data are correlated to define the rate-controlling processes which govern the observed localized and discontinuous precipitation reactions. An earlier paper5 was devoted to a discussion of the modes of heterogeneous nucleation of Ti2Cu from the martensitic a&apos; structure. The present paper will therefore emphasize the progress of the tempering reaction beyond the initial stages. EXPERIMENTAL METHODS The Ti-Cu alloys for this study were prepared by conventional arc-melting procedures. Chemical analysis revealed that the alloy buttons contained approximately 0.02 wt pct 0, 0.04 wt pct N, and within 0.1 wt pct of the intended copper concentration. The progress of the tempering reaction was followed by means of electrical-resistance measurements utilizing a Kelvin double bridge, microhard-ness readings with a 400-g load, and X-ray dif-fractometry patterns to reveal line-broadening effects. Metallographic specimens examined at magnifications greater than X1000 were shadowed with germanium to reveal fine structural details. Direct carbon replicas were prepared for the electron-microscopy studies. EXPERIMENTAL RESULTS Property Changes. The changes of microhardness that accompany the precipitation of Ti2Cu from a&apos; are shown in Figs. 1 and 2. As the solute concentration increases the peak hardness, for a given tempering temperature, increases. In alloys of given composition the time to reach the peak hardness agreed with the time to attain maximum X-ray diffraction peak breadth, as shown in Fig. 3. The maximum hardness and the maximum line breadth increased with lower tempering temperatures, and the time to reach the maximum also increased. The a&apos; peaks broadened during the initial stages of tem-

Jan 1, 1965

By J. M. W. Mackenzie, C. J. Wood

In attempting to describe the size classification performance of a hydrocyclone, most workers have elected to use either an equilibrium orbit theory or an non-equilibrium orbit theory. The equilibrium orbit theory has been used by the majority of workers including Lilge,&apos; Bradley; and Yoshioka and Hotta. In applying this theory, it is argued that particles in the body of a hydrocyclone attain an equilibrium radial position where the drag force on the particle resulting from the inward radial fluid velocity is balanced by the outward centrifugal force caused by the tangential component of fluid flow. When considered over the full height of the hydrocy-clone, attainment of this radial equilibrium orbit results in the particle following a conical equilibrium envelope. It is then argued that if this envelope lies outside the envelope of "zero vertical velocity," the particle will report to the underflow, while if the equilibrium envelope lies inside the envelope of "zero vertical velocity," the particle will report to the overflow or vortex finder product. The d50-sized particle which reports in equal quantities to the underflow and overflow is assumed to correspond to particles whose equilibrium envelope is coincident with the envelope of "zero vertical velocity." In considering the equilibrium orbit theory, it is apparent that the horizontal position of the particles in the feed inlet pipe should have no effect on their ultimate destination on the hydrocyclone. Each particle should attain an equilibrium position which depends on the density, size, and shape of the particle; the density and viscosity of the fluid; and the flow patterns within the hydrocy-clone. The nonequilibrium orbit or unsteady state theory has been largely developed by Rietema4 and Mizrahi.6 Mizrahi has listed four main objections to the equilibrium orbit theory. These objections center on the short residence time in the hydrocyclone, the fact that the experimental classification curve is much less sharp than is theoretically predicted, and the absence of negative efficiency conditions in hydro cyclones operating on a feed material which is much finer in size than d50. Proponents of the nonequilibrium orbit theory argue that for a particle to discharge with the underflow it must have sufficient outward radial velocity to reach the downward-flowing region close to the hydrocyclone wall in which the flow lines are parallel to the wall and the ratio of vertical to radial velocity is constant. It is then postulated that a d50 particle entering the cyclone at the center of the feed inlet will just reach this downward-flowing region as it reaches the apex. Thus for uniform distribution of particles across the feed inlet, half the d50 particles—that is, those injected in the half of the inlet area nearest the cyclone wall —will report to the underflow while those injected in the other half will not reach the downward-flowing region and will be carried inward to the center of the cyclone and thus report in the overflow. The exact thickness of the down-ward-flowing region of fluid adjacent to the outer wall of the hydrocyclone is uncertain but Mizrahi considers it to be equal to the apex radius minus the air core radius. Particles larger than d50 have a greater outward centrifugal force acting on them than the d50 particles and may reach the wall even if injected at a distance from the wall greater than Di/2 (Di is inlet diameter). Conversely, particles smaller than d50 may not reach the wall even if injected at a distance less than Di/2 from the cyclone wall. Since the equations put forward by the proponents of both theories yield approximately the same values of d50, it is not possible to decide between these theories by measurement of d50. It should be possible however to examine the theories by injecting a small stream of solids into the feed inlet of a hydrocyclone running on clear water. If the efficiency or classification curve is measured for various horizontal injection positions, then the curves should be coincident if the equilibrium orbit theory holds. If, however, the unsteady state theory describes the cyclone operation, then the classification curves should show finer d50 sizes for particles injected close to the cyclone wall. Experimental A 6-in.-diam hydrocyclone with geometry as in Figs. 1 and 2 was used. Quartz particles were injected as a 50% by wt pulp via an 1/8-in. steel probe. For each in-

Jan 1, 1971

By John D. Verhoeven

A technique has been developed for carrying out normal freezing experiments with a current density of 2000 amp per sq cut passing through the solid-liquid interface. The equation relating the effective distribution coefficient to the equilibrium distribution coefficient in electric field-aided solidification, originally developed by Huckc et al.,1 has been modified for the case of concentrated solutions. Preliminary experiments on the Sn-Bi system give qualitative agreement with the equation. The data are analyzed by a slightly novel use of the normal freeze equation which allows one to determine the effective distribution coefficient more easily. Very extensive mixing in the liquid was found at these high current densities and it is postulated that the mixing results from a vertical component of the magnetic Lorentz force generated by the electric current. In the search for techniques of obtaining ultrahigh-purity metals the inefficient but very effective technique of electrotransport has received little attention. Electrotransport is most effective in the liquid state and a natural application, therefore, is to apply an electric field across the liquid zone of a zone-melting experiment. The present investigation was undertaken to study the effect of an electric field upon solidification of metals, so that the usefulness of electrotransport in such solidification experiments as zone melting could be determined. In zone-melting and normal-freezing experiments it is difficult to achieve complete mixing in the liquid in the immediate vicinity of the solidifying interface. Consequently a solute build-up will occur at the interface in the portion of the liquid where complete mixing does not occur (an equilibrium distribution coefficient, ko, less than one, and unidirectional atom motion will be implied throughout). This local solute build-up produces a corresponding rise in the solute concentration in the solid so that the ratio of the solute concentration between the solid and the bulk liquid is larger than the equilibrium distribution coefficient. This ratio is defined as the effective distribution coefficient, k,. The differential equation describing the solidification process may be derived by applying the continuity equation to an expression for the net solute flux at the interface. The solution to this differential equation then allows one to determine the solute distribution in the liquid and the relationship between k0 and ke. One of the most useful solutions to this equation was first derived by Burton, Prim, and Slichter,&apos; in which they assumed that a) the equilibrium distribution of solute was maintained on the plane of the interface, 11) the solute build-up ahead of the interface in the liquid disappeared at a distance 6 from the interface, and c) the solute distribution in the liquid was invariant with time. The following well-known relation between ko and ke was then obtained, where R is the rate of solidification and D the diffusion coefficient of the solute in the liquid. This equation appears to correlate the data from a number of different types of solidification experiments very well. Application of an electric field across the solid-liquid interface can produce an additional flow of solute atoms as a result of the electrotransport. When the polarity of the field is such as to direct the electrotransport flux away from the interface the solute build-up may be diminished, even to the point of producing a solute depletion and a consequent ke smaller than ko. The quantitative description of this process and the resulting form of Eq. [I] was first given by Hucke et al.1* and then inde- where E is the electric-field intensity and U is the differential mobility, i.e., the velocity of the solute atoms with respect to the solvent atoms per unit electric field. Both authors follow the method of Burton, Prim, and Slichter in their derivation, the only difference being the additional electrotransport term in the flux equation. It has been pointed out1,3 that Eq. [2] predicts the possibility of a noticeable increase in the purification of materials by solidification in an electric field. The validity of Eq. 121 has not been checked experimentally and it is possible that other factors&apos; arising from the presence of an electric field across

Jan 1, 1965

By James H. Boettcher

INTRODUCTION A complex interaction of worldwide economic and political forces is increasingly requiring 3 primary participants for the successful development of large natural resource projects in developing countries. They are the host government of the developing country (the "host government"), one or a consortium of international natural resource companies, (the "Resource Company") and one, or more likely, a syndicate of international lenders. The host government typically controls the rights to its country&apos;s natural resources and establishes the economic ground rules for their development. However, because of worldwide inflation and the depletion of easily accessible, high-grade resources, the nominal and real capital development and operating costs of projects have been spiraling. Moreover, the depletion of high-grade and mineralogically uncomplicated resources also raises the technical and operating risks associated with resource recovery. Thus, despite increasing political self-consciousness and economic nationalism in developing countries, international resource companies are still often necessary project sponsors since they can provide, among other things, the required technical and management expertise and equity capital. The presence of these factors, and often the addition of worldwide market access in a project, is often critical in determining whether or not the project will be considered creditworthy by international commercial lenders. The international commercial lender is often the critical third party in project development, due in a large part to the ever escalating-costs of project developments relative to the financial resources of many companies and even a growing number of countries. In addition, lenders have other reasons for taking more than a passive interest in all aspects of a project&apos;s development. These reasons include the following: 1. Projects are often financially structured such that the primary security for a loan, once completion and performance standards are met, is the cashflow of the project rather than recourse to its sponsors. 2. Lenders, in addition to the "project risk" just described, are frequently asked to assume various elements of "political" risk which nay put them in a sensitive economic position between the host government and international resource firms (This is discussed later). 3. International banks today have large loan portfolios of developing country debt with final maturities ranging from one to fifteen years. The various current projects undertaken are the building blocks of most developing country economies that will hopefully contribute to the long run economic development of a country as well as generate the foreign exchange necessary to repay foreign loans in the future. In this regard it is useful for lenders to realize whether a given project "makes sense" in terms of a country&apos;s natural resource endowments and competitive advantages over other countries. The type of framework discussed in this article would help lenders to better make such assessments. In order to facilitate the continued international flow of capital and technology to developing countries and of readily available supplies of raw materials from them, it is important for each of these three parties to increase its awareness of the goals and objectives of the other parties as well as the bargaining process that results in the legal, financial, credit and fiscal structures by which a joint venture project ultimately proceeds. Each of the parties faces increasingly complicated accept/reject decisions when they are faced with difficult choices among. a wide range of alternative structural combinations each of which has different implications for risk and reward to them. Therefore, by increasing the mutual, 3-way awareness of the decision processes of the other parties, it should be possible to: • provide better information and results in the formation of more rational expectations by all parties going into negotiations, • facilitate the project formulation process by allowing parties to focus more clearly on the areas of common interest and conflict, • result more frequently in final structures and agreements that better reflect the economic and political realities associated with a project and thus • result in more creative, innovative agreements that should reduce the need for changes to them and thus promote the kind of long run political stability conducive to attracting capital and further investment.

Jan 1, 1982

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.&apos; 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&apos;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. &apos;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&apos; .,....(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 C. C. Miller, W. F. Kieschnick, E. R. Brownscombe

The possibility has been mentioned that large pressure gradients in a solution gas driven field caused by high production rates might lead to a reduction in the ultimate recovery obtainable compared to that which would be obtained by a very slow rate of production. In the present study the reservoir conditions accompanying a high rate of production and the corresponding ultimate recovery were compared with those obtained if the reservoir were produced at some marginal rate throughout its entire life. The method of calculation involved the application of fluid flow - material balance analysis to a series of successive steady state conditions in the reservoir. In the case studied very little difference in recovery was obtained at the abandonment pressure. An analysis of the basic factors involved indicates that the sane results would hold for considerable variation in properties of reservoir fluids, permeability of the formation or well spacing. Even if the reservoir has a water drive it appears that no harm will be done by- open flow production if the rate is cut back before any appreciable free gas is produced. However, any condition which leads to disproportionate withdrawal rates and which causes large pressure differences over large areas might result in substantial loss in ultimate recovery. INTRODUCTION When the pressure declines in a reservoir, gas is evolved from solution causing a shrinkage in the reservoir oil so that a stock tank barrel of oil with its residual dissolved gas occupies a smaller volume than originally. For a given reservoir oil saturation, therefore. stock tank oil content will be higher at the lower pressure. At the same time the evolved gas which remains in the reservoir reduces the reservoir oil saturation, and this would be expected perhaps to cause a reduction in stock tank oil content. However, around the well the shrinkage occurs so rapidly that stock tank oil may actually accumulate in this region in spite of the simultaneous reduction in reservoir oil saturation. If this accumulation should continue, spreading into a larger zone, the loss might be appreciable and because larger pressure gradients accompany higher rates the accumulation might possibly be large if wells are produced at too high a rate. It is the purpose of this paper to determine the effect of rate of flow on ultimate recovery by solution gas drive. The problem under consideration involves the unsteady state flow of oil in a closed radial system from which oil is produced by the evolution of gas as a finite drawdown is estalllished and the reservoir pressure declines. A somewhat similar unsteady state production problem dealing, however, with a linear system has been considered by Muskat and Meres, but these authors do not discuss the effects which rate may have on recovery. Other methods have been reported for estimaling recovery as a function of pressure decline when producing at an infinitesimal rate, and Muskat and Taylor&apos; have shown the effect of fluid and reservoir characteristics on recovery under these conditions. Assuming steady state to exist, Babson has given a procedure for computing the decline in well productivity as the field is depleted, and Moyer has combined these methods to estimate abandonment conditions and the resulting ultimate recovery. These latter approaches to the problem either do not consider the effect of rate or assume steady state can be established. In detail this paper presents the behavior of the reservoir producing under unsteady state conditions and considers the effect of production rate during this unsteady state flow on the recovery attained. The analysis also shows that there will exist at abandonment a type of steady state, and on this basis a method has been clevelopetl for rapid calculation of abandonment conditions.

Jan 1, 1949

By N. De Nevers

A calculation method has been developed for carbonated water flooding. This method takes into account the effects of oil viscosity reduction and oil swelling, due to carbon dioxide transferred to the reservoir oil from the carbonated flood water. It shows the effect of changes in carbonation pressure and carbonated water slug size. The method is based on a Buckley-Leverett-type linear flow model. The mathematical approach is similar to that developed by Welge, et al, for enriched gas drives. However, it contains a very important advance over the previous papers in its solution to the problem of injection of a "slug" of carbonated water rather than continuous injection of fluid of constant composition. The method also shows the chromatographic nature of the transport of Co 2 in a carbonated water flood. Sample calculations indicate that viscosity reduction is the most important source of incremental oil recovery. Swelling contributes to a lesser extent. Raising the carbonation pressure increases recovery, but this effect diminishes as pressure is increased. increasing the slug size also increases recovery, but again this effect diminishes as slug size is increased. INTRODUCTION In recent years there has been considerable interest in carbonated water flooding as a secondary recovery method. Several carbonated waterflood field tests are now in progress.1&apos;3 Yet there is no published method for predicting the results of a carbonated water flood. This paper presents a method for predicting the oil recovery for ; carbonated water flood. In a carbonated water flood. a carbonated water slug is injected into the reservoir. Then this slug is driven forward by injecting plain water behind it. The field operator must choose the carbonated water slug size and the CO2:water ratio in the slug. The calculation method presented here shows the effects of both slug size and CO2:water ratio on oil recovery. The calculation assumes that the addition of carbon dioxide to flood water increases the oil recovery by two mechanisms: 1. Carbon dioxide transfers from the flood water to the oil, thereby reducing the viscosity of the oil. At the lower viscosity, the oil flows more readily and is more easily swept out by flood water. 2. The carbon dioxide which transfers to the oil phase causes the oil phase to swell. Assume that the residual oil saturation (volume per cent of pore space) is the same whether the oil is carbonated or not. Then the same volume will contain less stock tank oil if it is carbonated, due to oil swelling. The plain water driving the carbonated water slug forward will strip the CO2 out of this residual oil., causing this oil to shrink. This swelling-followed-by-shrinking leaves behind a lower residual oil saturation than would have been left behind by a plain water flood. Since some oils swell to almost twice their original volume at reasonable CO2 pressures, this effect may be large. Two additional benefits sometimes claimed for carbonated water flooding are: 1. Adding CO2 to flood water greatly increases the water injectivity in some cases.1 (Injectivity = BPD of water injected/ft of sand x differential bottom-hole pressure, psi.) CO2 forms a weak acid in water; injecting carbonated water is, therefore, somewhat like acid treating. If a formation responds to acid treating, it will probably have a higher injectivity for carbonated water than for plain water. Although the increased injectivity increases the profitability of a carbonated water flood over a plain water flood, it does not change the production curve based on total fluid injected. Therefore, the method presented here does not show any effect of increased injectivity. 2. It is claimed in one patent4 that carbon dioxide will react with some constituents of the oil to form detergents. These detergents are expected to increase the oil recovery. There is little data available to support this idea. The calculation method presented here takes no credit for extra oil recovery due to the formation of any such detergents.

Jan 1, 1965

By E. I. Organick, B. H. Golding

A simple correlation is presented for the prediction of saturation pressures in condensate-gas and volatile-oil mixtures. Saturation pressure is related directly to the composition of the mixture with the aid of two generalized composition characteristics, B, the molal average boiling point, and W, the modified weight average equivalent molecular weight. The need for hydrocarbon equilibrium constants is entirely eliminated. Only the values of the two conlposition parameters need be obtained, and these are calculated in a simple and reproducible manner from the stream analysis of the mixture and the ASTM distillation for the heavy ends. The correlation is given .in the form of a set of 14 working charts in which saturation pressures, either retrograde dew points or bubble points, are plotted cs temperature, forming partial phase envelopes for mixtures having discrete values of the composition parameters, B and Wul. A locus of critical states on each chart intersects the partial phase envelopes, distinguishing retrograde dew point pressures from bubble pressures. A comparison Between experimental and predicted saturation pressures reveals that: the probable error of the predicted values, when considering the data as a whole (214 experimental points), is five per cent; for high gasjoil ratio condensate fluids (approximately 40.000 or more cu ft/bbl) the probable error is about 8 per cent. In either case, however, the probable error is comparable with the overall uncertainty of the data from which the correlation is derived. INTRODUCTION There is a present need for a simplified, yet reasonably accurate method for predicting upper phase boundaries of reservoir fluids. In the case of high gas/oil ratio fluids, namely the so-called highly volatile oils and the condensate gases, no simple method, accurate within 5 to 10 per cent, is now available by which bubble point and retrograde dew point pressures may be predicted. In order to calculate the upper phase boundary of a volatile reservoir fluid, using correlations now available, it is necessary to find and employ the proper K-values (equilibrium volatility ratios of the hydrocarbons in that fluid). This process of finding the Proper K-values usually involves a double trial and error computation. Since the K-values are in themselves dependent upon the composition of the fluid, their use in the calculation of a saturation pressure actually represents an unnecessary complication. The correlation described in this paper relates saturation pressure directly to the stream composition, entirely eliminating the need for K-values. The only data required to calculate the Stream composition are obtained from the customary low temperature fractional analyses of separator gas and liquidt and from the ASTM distillation of the c,+ fraction of the separator liquid. A considerable simplification in the procedure for predicting saturation pressures is thereby effected without significant loss in accuracy. BASIS FOR A CORRELATION Empirical relationships between the saturation pressure and the composition of the fluid in question have already been partially explored by one of the authors.V n particular, a relationship between the critical pressure for a hydrocarbon mixture and its composition was developed. The critical pressure is in reality a saturation pressure, as it repregents a limiting case either of a dew point or of a buhhle point pressure for any given mixture. It was found that the critical pressure may be correlated as a function of two generalized composition characteristics or variables. One characteristic expressed a molal average prop-erty, the other a weight average property of the mixture. Specifically, the molal average boiling point (MABP, or B as it

Jan 1, 1952

By R. W. Snyder, H. J. Ramey

This paper presents the results of applying the Buckley-Leverett&apos; displacement theory to petroleum reservoirs consisting of a finite number of layers. The layers are assumed to communicate only in the wellbores, and the reservoir may be represented as a linear system. Most previous investigations of this nature were limited by assumptions and by inconsistent calculation techniques. This study improves on previous work by applying the Buckley-Lev-erett displacement theory to a noncommunicating layered system where permeability, porosity, initial saturation, residual saturation and relative permeability vary from layer to layer in a logical and consistent manner. Gravity and capillary-pressure effects are neglected. A modification of the Higgins-Leighton calculation method was used in this study. Waterflood predictions were made with all properties varying, and then with only permeability varying using several inability ratios. These results were compared with the Stiles and Dykstra-Parsons predictions. It is shown that the latter methods generally give poor values for the breakthrough recovery and pessimistic predictions for the performance after breakthrough. Similar results were obtained for a gas-displacement case. lNTRODUCTION Field experience with immiscible displacement usually shows constant producing conditions until breakthrough of the displacing fluid. Then oil production continues at increasing displacing-to-displaced fluid ratios until the economic limit is reached. Three different ideal mechanisms are known that will produce this behavior: (1) relative permeability effects as described by Buckley-Leverett frontal advance theory,&apos; (2) vertical stratification as considered by Stiles,2 Dykstra and Parsons5 and others and (3) different path lengths involved in areal (two-dimensional) flow between wells as described by Dyes et al.4 Without question, a combination of these factors modified by formation heterogeneity and other known and unknown factors actually does control the behavior of real systems. This paper presents results of an investigation of certain factors that should affect performance but which have received little attention to date. In 1944, Law5 demonstrated that porosity and perme- ability are often found to have normal and logarithmic-normal distributions, respectively. throughout cored intervals in natural formations. This led to the concept of the noncommunicating, multilayered reservoir model for immiscible displacement. This model assumes that the reservoir is composed of a number of layers that communicate only at the wellbores. Each layer is individually homogeneous, but may be different from every other layer. Stiles&apos; presented one of the earliest applications of this model to waterflood performance. In addition, Stiles assumed that the initial saturations and relative permeabilities were the same for each layer, porosity was the same. displacement was piston-like, fluids were incompressible and injection into each layer was proportional to that layer&apos;s permeability capacity (permeability-thickness product). The last assumption would be true if the mobility ratio for the displacement were unity.21 Dykstra and Parsons" used the same model as Stiles, but rigorously included mobility ratios other than unity for piston-like displacement. Dykstra and Parsons used their general result to produce charts for log-normal permeability distributions between layers. Similarly, Muskat6 Pub1ished analytical solutions for linear and exponential permeability distributions. In 1959, Roberts&apos; described a scheme for calculating water-drive performance for the noncommunicating, layered reservoir model which considered two-phase flow in the displaced region. Roberts used the same model and assumed that the injection rate into a layer was proportional to that layer&apos;s permeability capacity, but that flood front locations could be evaluated from the Dykstra-Parsons results. These assumptions are inconsistent, and a material balance cannot be maintained except for a mobility ratio of unity. At the same time, Kufus and Lynch8 coupled Buckley-Leverett displacement theory with the layered model to provide an improvement of the Dykstra-Parsons method that was consistent. In 1960, Higgins and Leighton9 resented a numerical method for calculating waterflood performance also considering two-phase flow in the displaced region. The result was used to investigate variation in absolute permeability and oil viscosity. An excellent, detailed history of using the noncomrnunicating, layered reservoir model was presented by Nielsen.&apos;" The preceding techniques (and many related ones) were similar in that differences in initial saturations, residual saturations and relative permeabilities from layer to layer were neglected. It is well known that the irreducible water saturation is an important function of absolute permeability. Calhoun11 showed that the irreducible water saturation

By M. R. Tek, R. L. Nielsen

The scaling laws as formulated by Rapoport relate dynamically similar flow systems in porous media each involving two immiscible, incompressible fluids. A two-dimensional numerical technique for solving the differential equations describing systems of this type has been employed to assess the practical value of the scaling laws in light of the virtually unscalable nature of relative permeability and capillary pressure curves and boundary conditions. Two hypothetical systems — a gas reservoir subject to water drive and the laboratory scaled model of that reservoir — were investigated with emphasis placed on water coning near a production well. Comparison of the computed behavior of these par ticular systems shows that water coning in the reservoir would be more severe than one would expect from an experimental study of a laboratory model scaled within practical limits to the reservoir system. This paper also presents modifications of the scaling laws which are available for systems that can be described adequately in two-dimensional Cartesian coordinates. INTRODUCTION Present day digital computing equipment and methods of numerical analysis allow realistic and quantitative studies to be carried out for many two-phase flow systems in porous media. Before these tools became available the anticipated behavior of systems of this type could be inferred only from analytical solutions of simplified mathematical models or from experimental studies performed on laboratory models. To reproduce the behavior of a reservoir system on the laboratory scale, certain relationships must be satisfied between physical and geometric properties of the reservoir and laboratory systems. Where the reservoir fluids may be considered as two immiscible and incompressible phases, the necessary relationships have been formulated by Rapoportl and others.2-5 Rapoport&apos;s scaling laws follow from inspectional analysis of the differential equation describing phase saturation distribution in such systems. It will be recalled that these scaling laws presuppose three conditions: (1) the relative permeability curves must be identical for the model and prototype; (2) the capillary pressure curve (function of phase saturation) for the model must be linearly related to that of the prototype; and (3) boundary conditions imposed on the model must duplicate those existing at the boundaries of the prototype. These three requirements seldom if ever can be satisfied in scaling an actual reservoir to the laboratory system because: (1) The laboratory medium normally will be unconsolidated (glass beads or sand) while the reservoir usually is consolidated. Relative permeability and capillary pressure curves are usually quite different for consolidated and unconsolidated porous media. (2) The reservoir usually will be surrounded by a large aquifer which could be simulated in the laboratory only to a limited extent. (3) Wells present in the reservoir would scale to microscopic dimensions in the laboratory if geometric similarity is to be maintained. In view of these considerations, rigorous scaling of even a totally defined reservoir probably would never be possible. The purpose of this paper is to assess the practical value of the scaling laws in the light of the unscalable variables. This has been done by carrying out numerical solutions in two dimensions to the differential equations describing the flow of two immiscible, incompressible fluids in porous media for a field scale reservoir and a laboratory model of that reservoir. While both the reservoir and the laboratory model were purely fictional, each has been made as realistic and representative as possible. The field problem selected as the basis for the investigation was an inhomogeneous, layered gas reservoir initially at capillary gravitational equilibrium and subsequently produced in the presence of water drive. The laboratory model of this reservoir was designed to utilize oil and water in a glass bead pack. The numerical treatment employed was similar to that of Douglas, Peaceman and Rachford6 and it included both capillary and gravitational forces as

By W. A. Tiller

The unidirectional freezing of a semiinfinite liquid from one end has been treated by calculating the solute and temperature distributions in the liquid and solid phases with time, when the solid-liquid interface advances at a rate proportional to (time)lt2. The nucleation probability for new grain formation is calculated as a function of the constitutional and thermal properties of the alloy, the pouring conditions and the catalytic efficiency of the nucleating agents. The transition from columnar crystal formation to equiaxed crystal formation is assumed to occur at a critical nucleation frequency for a particular alloy and the quantitative dependence of the ratio (columnar zone length/equiaxed zone length) upon the freezing conditions has been predicted. The effects of (i) fluid flow in the liquid, (ii) heat transfer through a mold wall, and (iii) a finite ingot length, have been explicitly tveated as departures from ideality. These departures from ideality are found to have very important consequences on the ingot stmcture. In Part I&apos; of this series, a qualitative description of the general ingot structure was developed. The factors that determine i) the ratio of the equiaxed zone length to columnar zone length and ii) the grain sizes in these zones, were discussed. In the present paper, topic (i) will be treated quantitatively, and it is intended that topic (ii) will be treated quantitatively in Part 111. The complete treatment has been split in this manner for ease of development and coherency of presentation. The treatments are quite mathematical in nature and, although some of the conclusions to Part 11 depend upon the treatment presented in Part 111, it was felt that the intended audience would find the material more readable and more useful in its present form. To this same end, it seems worthwhile to reiterate some of the qualitative aspects of grain nucleation during ingot solidification before embarking on the mathematical description. QUALITATIVE ANALYSIS In the conventional method of casting, a liquid alloy with a certain degree of superheat is poured into a mold. The outer rim of liquid metal cools rapidly to the liquidus temperature of the alloy and begins to undercool. As the degree of undercooling increases, nuclei of solid begin to form in this chill layer, this nucleation being catalyzed by certain heterogeneous particles in the liquid. The grains in the chill layer which are generally of random orientation, begin to grow inwards because of heat conduction through the mold. The grains usually develop a columnar shape since the frequency of nucleation of new grains in the liquid ahead of the advancing interface is generally insufficient to seriously impede the growth of the original grains. Due to the slow diffusion rate of solute in the liquid, a solute distribution is formed in the liquid adjacent to the interface which is either an excess or a depletion if the partition coefficient, ko, of the solute is < 1 or > 1, respectively. Thus, the liquidus temperature of the liquid adjacent to the interface is a function of position. The actual temperature in the liquid is also a function of position and the temperature gradient at the interface decreases as the interface moves away from the mold wall. Since the rate of advance of the interface decreases with time, the equilibrium temperature distributions and the actual temperature distributions as a function of interface position, X, can be illustrated as in Fig. 1. Thus, a region of liquid may exist in which the actual temperature is below the liquidus temperature. This liquid is contitutionally supercooled and therefore unstable with respect to the nucleation of new grains. Now, since the supercooling in Fig. 1 increases as the interface moves further away from the mold wall, the nucleation rate of new grains will increase. When the nuclea-

Jan 1, 1962

By K. E. Barrett

Exxon owns a uranium mill and holds two mining leases in Live Oak County, Texas, about halfway between San Antonio and Corpus Christi. The properties make up the Felder Uranium Operations which was reopened earlier this year. Exxon held an oil, gas, and other minerals lease on the J. C. Felder tract, which was adjacent to a relatively shallow uranium discovery by Susquehanna-Western, Inc. on the Marrs-McLean lease immediately south of the Felder property. Drilling in 1967 and 1968 confirmed the presence of reduced uranium mineralization in the basal sand unit of the Oakville formation on the Felder tract, which formed the major part of the roll-front deposit. In 1969 Exxon and Susquehanna-Western, Inc. entered into a sale and purchase agreement which provided for Susquehanna to mine and process Felder ore and purchase recovered uranium. Susquehanna moved an alkaline-leach mill from Wyoming, erected it on the Ray Point property, and placed it into operation late in 1970. Susquehanna mined and processed ore from the Felder and McLean properties through March 1973. Susquehanna ceased operations in March 1973. Exxon then acquired the mill and mill property. Exxon also purchased the mineral rights to the McLean lease, re-negotiated a mining lease for that property, and carried out shut-down programs for the mining and mill areas in the fall of 1973. The project was put on a standby basis until late 1973, when Exxon initiated mine feasibility studies for the project. MINE PLANNING EVALUATION The feasibility study for reopening the Felder Project began in late 1975 and was not completed until late 1976. I will discuss several areas of the feasibility study that required additional work prior to making the decision to renew operations. Ore Reserves Preparations for estimating the ore reserves began with the re-evaluation of more than 1500 natural radioactivity logs from exploration and pre-development drilling that had been completed on the property. These gamma ray logs of non-core rotary drill holes were the principal source of data used in making the estimate. Chemical assays of cores from the deposit were also used in the reserve determination. Electric resistivity and self-potential logs were made along with the gamma ray log. In December 1975 an additional core drilling project was undertaken to confirm the in-place density and radiometric equilibrium characteristics of the ore deposits. Comparison of chemical assays of cores with the U308 values calculated from the logs showed that the unoxidized ores were in radiometric equilibrium. In contrast, cores from anomalies occurring in near surface, weathered, and oxidized zones were in radiometric disequilibrium. Several important decisions were made in developing a mine plan or schedule of production from the Felder and McLean ore bodies. Disposal of Produced Mine Water: The ore bodies of the Felder Uranium Project occur at a point below the ground water table. The ore zones to be mined must first be dewatered to allow removal of mineralized material. In the open pit operations, this is accomplished by maintaining a perimeter ditch around the periphery of the open pit, allowing the interior of the pit to drain and collect into a sump and be pumped from the mine. In addition to anticipated water production from future mining operations, approximately 200M gallons of water was contained in three open pits left from prior mining operations. In two of these existing pits, the water was to be removed and disposed to allow for planned backfilling of waste material into these pits. The third pit would also have to be drained to allow continued mining of an area left from the prior operations. Essentially no ground water information was available for this area. The only data available was water production history from Susquehanna&apos;s mining operation. Two water wells were drilled early in 1976 on the Felder lease for use in obtaining hydrological data. A long term draw-down test was performed by pumping one water well and measuring water level drawdown in both the pumped well and the observation well. From these data, values for permeability and storage coefficient were calculated. These data were then used in modeling the aquifer to allow calculation of water influx into the mining area versus time. Once a schedule of water production, including the stored volume in the existing pits was calculated, alternate solutions for disposal were evaluated. The first system evaluated was a series of deep injection wells. The wells were designed to inject at a depth of approximately 3500 feet. Again very little information concerning reservoir characteristics of the receiving sand units was known. Using assumed values for reservoir permeability and storage coefficients, an injection well system was designed to allow for disposal of produced mine water. The biggest

Jan 1, 1979