Search Documents

By Orville R. Lyons

THE relative performance of coal washing equipment, or the effectiveness with which any type or make of equipment removes impurities from coal, has been most difficult to evaluate in the past. The most widely used yardstick is the Frazer and Yancey efficiency formula developed in 1922,1 but Yancey in a later article states that "washers treating coals of different density composition or operating at different densities of separation cannot be compared directly on the basis of this criterion."' Prior to and since 1922, a variety of other methods has been used for comparison purposes, including the distribution curve, the error area, and the "ecart probable" or probable error. Yancey and Geer in discussing these methods conclude, "Performance can be evaluated in a number of different ways, with the choice of the proper method to use being dictated by the objectives of the investigation and the data available."' It is true that performance can be evaluated in a variety of ways, but if the equipment is to be evaluated on an effectiveness basis, there should be only one universal comparison method. Varying methods have been used because one universal comparison method has not been found or developed. In the article previously quoted, Yancey and Geer state in clear terms the primary concept for a universal comparison method: "One of the simplest, and certainly one of the most obvious evaluations of washery performance is the quantity of sink material in the washed coal and the float material in the refuse. If the washery products are tested at the density at which the washing unit is operated, the sink in the washed coal and the float in the refuse represent material that has been misplaced." The quantity of misplaced material was used as a criterion of washery performance by Lincoln in 1913,8 by the United States Bureau of Mines in 1938,4 by Hancock in 1947,3 and by the national French research agency Cerchar in recent years.2 In 1950 Anderson' proposed the use of this criterion as an efficiency value to replace the Frazer and Yancey formula. However, none of the above-mentioned investigators used the misplaced material concept in a manner that would provide universal coal-cleaning equipment comparisons. The Correlation Theory The ideal coal cleaning process would treat all sizes and would make a perfect separation at any given specific gravity. All material lower in density than the desired value would report in the coal product and all material higher in density would report in the refuse product. Unfortunately, no known cleaning process achieves this goal and there seems little likelihood that any process yet to be invented will do more than approach it. When coal is treated in volume under operating conditions, it is impossible to avoid mechanical entrapment, fluctuations in throughput and effective gravity of separation, and the creation of turbulent currents, even when a true heavy-liquid bath is used and the feed is closely sized and contains little intermediate gravity material. This being so, it is possible to appreciate the difficulties inherent in trying to obtain a perfect separation when treating a wide range of sizes and a feed containing high percentages of intermediate material, using turbulent currents to help create the effective separation gravity, under operating conditions which normally tend to be on the overload side. When coal is separated from refuse in any coal cleaning equipment, some refuse always reports to the coal and some coal to the refuse; the writer therefore assumed that there should be a relationship between the total amount of misplaced material produced by any given piece of equipment and the difficulty of separation as represented by the percentage of near gravity material in the feed. With small amounts of near gravity or ±0.1 material in the feed there should be less misplacement of material than would occur with large amounts of near

Jan 1, 1952

By N. A. Gjostein

The formal descnption of the decay of an isolated scratch can be written in terms of an appropriate Fourier integral. With the application of certain approximations, this description leads to the second-moment analysis developed by King and Mullins. It is shown that the second-moment analysis does not adequately describe the behavior of a decaying, isolated scratch, and that it is difficult to obtain internally consistent values of the surface self-diffusion coefficient, D,,from this analysis. This inconsistency arises from contributions of higher-order moments to the decay process. Two new analytical methods, based on the Fourier integval representation of a scratch profile, have been developed for determining values of D, from the decay kinetics of isolated scratches. In the profile-matching method, a given initial profile is resolved into its Fourier components, the components are allowed to decay exponentially over an annealing time interval, At, and finally they are synthesized to form a decayed profile. The predicted profile is adjusted to fit the experimental profile by varying only the parameter, D,At; the contribution from volume diffusion is directly taken into account and the integrals involved in the Fourier analysis are evaluated by numerical integration on an electronic computer. In the kinetic method, the widths, W, of a series of predicted profiles are determined for various volume-diffusion contributions, and the results are plotted in the form W vs DSAt. For a fixed ratio of surface to volume diffusivities, R, these plots are linear; the slope of the lines increases as R decreases. A simi-lar analysis can be performed on the depths of the predicted profiles. In this case, plots of (l/df3 vs D,At are found to be nonlinear. The proper value of Ds is found by superimposing the data for the experimental widths and depths on Plots of this kind and adjusting Ds until a fit is obtained with the line having the appropriate value oj R. An analysis of the slopes of these plots for R —= (surface diffusion only) reveals that they do not give the correct value, 2276, for width tf second-moment standard profile. This finding allows one to modify the second-moment analysis in such a way as to obtain internally consistent values of Ds. ThERE are several processes, involving a flow of matter along a surface, due to gradients in curvature of the surface, that may be used to determine the surface self-diffusion coefficient, Ds. These processes are commonly known as 1) thermal grooving, 2) single and 3) multiple scratch decay, 4) blunting of field-emitter tips, and 5) sintering. In this series of papers, the first three of these processes are examined in detail in an attempt to subject the phenomenological theories governing these processes to a critical test and, furthermore, to make a critical comparison of the values of D, obtained from these techniques. Although the thermal grooving technique has been established'' as a reliable method for determining D, , this has not been the case for the process of scratch flattening. Mullins and coworkers374 have proposed phenomenological laws governing the growth of the width and the decay of the amplitude of an isolated scratch, but only in one study5 has the width-growth law been supported by experimental evidence. In paper I of this series detailed experimental evidence is presented that points up the deficiencies of the second-moment analysis; the origins of these deficiencies are also discussed. Following this development, some new analytical procedures, based on representing the scratch profile by Fourier integrals, are presented as a means of obtaining reliable values of D, from the kinetics of the scratch-flattening process. Further applications of the Fourier integral approach are reported in papers 11' and 111,~ where the method has been useful in interpreting the anisotropy of flattening rates of isolated scratches (paper 11) and in describing the multiple scratch decay process (paper 111). THEORY OF SCRATCH DECAY Following ~ullins" approach we consider an arbitrarily shaped scratch profile, Fig. 6, at some time, t. The profile is designated by a function z(x,t) which can be expressed in terms of a Fourier integral as follows:

Jan 1, 1967

By John W. Weeton, Donald W. Petrasek

Relatively few metal-metal systems exist that would permit the creation of fiber-metal composites consisting of mutually insoluble constituents. It is anticipated that most high strength-to-weight ratio and high -temperature fiber composites ultimately to be produced will utilize high-strength fibers embedded in a highly alloyed matrix. An investigation was conducted to determine the effect of alloying on the tensile Properties and micro-structure of tungsten-fibey-reinforced composites. Composites were made of tungsten fibers infiltrated with copper binary alloys that contain elements of varying solubility in tungsten. Room-temperature tensile tests were made on the composites, and a metallographic study of the microstructure of the fiber-metal-matrix interfaces was conducted. It was shown that, as the depth of penetration of some alloying elements into the fiber increased, a decrease in both tensile strength and ductility of the composite results. A notch effect, due to the formation of a brittle alloy zone with the tungsten fiber, was observed and gave rise to greater reductions in tensile properties and ductility behavior than would be Predicted by a law-of-mixture relation alone. In recent years, investigators have considered combining fibrous materials with relatively weak binder materials. The interest in fiber-reinforced composites results from the fact that fibers or wires may be exceedingly strong and exhibit mechanical properties superior to those of the bulk materials from which they are derived. For example, tungsten wires drawn to less than 1 mil in diameter may have strengths of over 600,000 psi,' while steel wires ranging from 10 to 3 mils in diameter may have strengths of the order of 300,000 to over 600,000 psi.2 Both metallic and ceramic whiskers are known to have very high strengths. For example, iron whiskers have been shown to have strengths of 1.9 million psi3 and sapphire whiskers have been reported to have strengths of the order of 1.7 million psi.4 If all or part of the strengths of such fibers could be retained, subsequent to the incorporation of the fibers in an engineering material, a superior material could result. In previous work done at the Lewis Research Center, McDanels, Jech, and weeton5 demonstrated the feasibility of fabricating metal-fiber—metal-matrix composites that utilized the full strengths of the fiber and matrix materials. Tungsten fibers were combined with a copper matrix by a liquid-infiltration method such that the tungsten fibers were uniaxially oriented parallel to the tensile axis of the composite. A linear relation was established that related the tensile strength of the composite with the volume percent of the fibers. Mutually insoluble constituents were selected for the study so that the strength of the composite could be quantitatively related to the strength of the individual constituents. The fact that there are only a few metals that are mutually insoluble with tungsten, e.g., copper, silver, gold, and zinc, and that these either have low melting points, high vapor pressures, or poor oxidation resistance, makes it necessary to utilize other metals that are soluble in tungsten, where high-temperature applications are ultimately contemplated. It was anticipated that, for most practical fiber-reinforced composite materials, it would be necessary to use a fiber and metallic matrix that would have varying degrees of solubility for each other. A logical extension of the earlier work was thus felt to be a study that would determine the effects of elements that had mutual solubility with the tungsten. Numerous investigators have found that tungsten fibers may be damaged by thermal treatments or contamination of the surface of the fiber by alloying elements.6-10 On the other hand, surface treatments may improve the strength of such fibers."-" Relatively little work with composites containing tungsten fibers have been reported. The work of parikh14 has shown that kinked, short-length tungsten, steel, and molybdenum fibers could str'engthen such matrices as silver, gold, and copper. Some of the combinations of fibers and matrices studied by

Jan 1, 1964

By J. E. Steiner, J. M. Hodge, M. A. Orehoski

Ingots of four steels (1045, 1080, Ni-Mo-V, and Ni-Cr-Mo-V) were cast at pressures varying from about 1 to 760 mm of mercury, so as to obtain a range of hydrogen contents in each steel. The susceptibility of these steels to flaking was evaluated by cooling hot-rolled blooms at a range of cooling rates, varying from water quenching to slag cooling, and studying the occurrence of flaking as a function of cooling rate. The steels exhibited differences in susceptibility to flaking, and a tlzreshold value for hydrogen content for each steel was found, below which the susceptibility was very low and above which the susceptibility increased markedly with increasing hydrogen content. Mechanical property evaluations, metallogvaplzic exantinations, and X-my diffraction determinations o,t retained austenite were made on samples from the blooms. A mechanism of flaking which is consistent with the observed behaviors, the effects of hvdrogen content, and the differences in susceptibility among the steels investigated is postulated, based on the results of mechanical property evaluations and metallographic and X-ray diffraction studies. ThE important role of hydrogen as a factor in the occurrence of "flakes" or "hair-line cracks" in steel has now been established by many investigations. These have included investigations by steel producers and forging manufacturers, to whom the prevention of flaking is a very real and serious problem, and many research studies aimed primarily at establishing the mechanism of this behavior Several investigators3 have reported that steels cast at reduced pressures (absolute pressures of about 1 mm) not only are low in hydrogen content but are essentially insensitive to the formation of flakes. These investigations have shown that a general relationship exists between hydrogen content and susceptibility to flaking. The difficulty of obtaining suitable samples of steels of a given composition with the desired range of hydrogen contents, particularly in the relatively large sections in which flaking most commonly occurs, has been the principal obstacle to quantitative evaluations of these relationships. Steels with intermediate hydrogen contents between 1 ppm, and the 4ppm or higher of air-cast steels are generally not available in commercially produced steels. Laboratory procedures, in which the hydrogen is introduced into solid steel samples by heating in a hydrogen atmosphere under pressure, generally involve a limitation in the size of the sample which can be treated, and thus are not representative of the behaviors in large sections. The availability of a facility at the South Chicago Works of U.S. Steel, in which several ingots could be cast from a single heat of steel at total air pressures ranging from 0.25 to 760 mm of Hg, has made it possible to overcome this obstacle and to obtain ingots of a given steel covering the intermediate range of hydrogen contents desired. These ingots, furthermore, are of a size suitable for the production of product with section sizes large enough to be representative of the products in which flaking most commonly occurs in commercial production. In these experiments, this facility was used to study the relationships between hydrogen content and the susceptibility to flaking of four steels, two plain carbon steels at 0.45 and 0.80 pct C and two low-carbon alloy steels representing a medium and a rather high hardenability level. These flaking studies were supplemented by evaluations of mechanical properties and by metallographic studies on samples from the blooms, in order to obtain pertinent information on the factors other than hydrogen content which influence the susceptibilities of these steels to flaking. Based on the results of these tests, a mechanism for the initiation and propagation of the flaking cracks, which is consistent with the observed behaviors, is postulated. MATERIALS AND EXPERIMENTAL WORK Steelmaking and Conversion. To provide material for the investigation, five 29- by 35-in., hot-topped, big-end-up ingots were cast from each of four 80-ton, basic electric-furnace heats of four different steels (a 1045 steel, a 1080 steel, a Ni-Mo-V steel, and a Ni-Cr-Mo-V steel). A double-slag practice was used on all steels and the two plain carbon steels were Si-A1 deoxidized with an addition of 2 lb of A1 per ton, while the Ni-Mo-V and Ni-Cr-Mo-V steels were deoxidized only with silicon and

Jan 1, 1964

By D. E. Mikkola, G. E. Poquette

X-ray diffraction was used to study the growth of antiphase domains in quenched (or "disordered&apos;? Cu3Au annealed in the range 300" to 385°C. Measurements of the long-range order parameter indicate that a high degree of order is established rapidly at these annealing temperatures. Average domain sizes obtained by Fourier analysis of the superlattice peak profiles show that antiphase domain growth is analogous to classical metallurgical grain growth (i.e.,fol-lows a D2 US kt relation). The activation energy for the growth process, 44 kcal per mole, indicates diffusion control. The anisotropy in the observed effective domain sizes from seven different superlattice peaks can be best accounted for by a configuration involving antiphase domain boundaries on {100 } and (111 ) planes, where the (100) boundary is the low-energy type, {lOO} Type I. The disappearance of the (111) antiphase domain boundary during domain growth is dif-fusion-controlled with an activation energy of 48.5 kcal per mole. On the other hand, the behavior of the parameter representing the amount of (100) Type I boundary indicates that there may be a conversion of the higher-energy boundary types to (100) Type I boundary during the growth process. THE order-disorder transformation in the alloy Cu3Au has been the subject of numerous experimental investigations over the past 40 years, most of which have involved measurements of convenient physical properties such as resistivity and hardness as a function of heat treatment. The interpretation of these data is difficult because measurements of this type have the disadvantage that they yield total changes resulting from several different atomic processes which may be occurring simultaneously. The present study was undertaken to establish by means of X-ray diffraction the kinetic behavior of one of these atomic processes, namely antiphase domain growth. In the completely ordered state Cu3Au has an L12 structure which can be described as an fcc type arrangement of atoms in which the gold atoms occupy the corner sites and copper atoms the face-centered sites. Three other equivalent unit cells can be formed by allowing two gold atoms to occupy opposite face-centered sites and copper atoms to occupy the remaining corner and face-centered sites. These four different unit cells are crystallographically related by a shift of 1/2<110>. Quenching from above the critical temperature (392°C) to room temperature retains the "disordered" structure-a "random" arrangement of copper and gold atoms having the fcc structure. Annealing the disordered alloy below the critical temperature causes ordering to occur through the nucleation and growth of ordered regions at different points throughout the crystal. These small ordered regions, each based on one of the four types of unit cell, grow by consuming the disordered material until they impinge to form a domain structure in which the domains are antiphase by a shift of 4(110). The boundaries between domains are called antiphase domain boundaries (APDB) and they have been found to lie primarily on &apos; {100} planes in highly ordered Cu3Au. Two types of APDB can occur on ( 100) planes. An example of an (001)1/2[110] shear type APDB is shown in Fig. I. Because of the layer of copper atoms separating the domains in this case, there are no near-neighbor violations across the boundary and it has a low surface ener This type of APDB will be referred to as {lOO} Type I. Fig. 2 shows an (001): [011] climb type APDB, (100) Type 11. This type of APDB can be described as resulting from removal of a (200) plane followed by a shift of 1/2<110> to bring the separated regions together. Near-neighbor violations across (100) Type I1 APDB make it a high-energy boundary. There are numerous other possible APDB configurations; however, it is only worth noting that shear type APDB can be produced on the (111) plane by motion of a normal fcc dislocation with Burgers vector 1/2<110> through the ordered structure. A (111); [110] shear type APDB is shown in Fig. 3. Again, because of the near-neighbor violations this type of APDB has a higher energy than the (100) Type I APDB.

Jan 1, 1970

By A. K. Mukherjee, J. E. Dorn, J. D. Mole

Investigations were carried out on the effect of polyslip on the strain hardening of aluminum single crystals. The orientations investigated were those lor which the tensile axis was in the [001], [111], [112], and [012] directions plus another for which the Schmid angles for {111}(110) slip were 1 deg. The experimental data were analyzed on a model based on the intersection of dislocations with particular emphasis on the effect of polyslip on the activation volume for inter-section. It is shown that the rate of strain hardening inreases for those orientations wherein attractive dislocation intersections occur and that those orientations which produce the greater number of such intersections exhibit the greater strain hardening. Good correlation of the data is obtained with the concept that attractive junctions, as proposed by Saada, Play an important role in accounting for the rate of strain hardening. EXISTING concepts on the nature and cause of strain hardening in fcc metals have been deduced principally from experiments on the deformation of single crystals under single slip. The effect of crystal orientation on the shapes of the stress-strain curves for single slip have been summarized by seegerl and more recently by Clarebrough and Hargreaves.2 Tensile specimens whose axes fall near the center of the [001]-[011] line of the standard triangle of the stereographic projection exhibit the longest range of easy glide (Stage I) and the lowest rates of linear hardening (Stage 11) whereas specimens whose axes lie near fie [001]-[111] line of the standard triangle have limited or no easy-glide range and exhibit somewhat higher linear strain-hardening rates. Specimens whose axes lie near the [001] or the [ill] poles do not exhibit easy glide and have the highest rates of linear hardening. Kocks3 has shown that the highest rates of linear hardening Occur under Polyslip when the tensile axis coincides with the [111] or the [ 001] pole. Since the rate of strain hardening is sensitive to specimen orientation and the incidence of polyslip, these relationships might help to discriminate between various dislocation models for strain hardening in fcc metals. Previous attempts to analyze the effect of orientation on strain hardening,1"3 however, did not provide a unique answer to this problem. Consequently the present investigation was undertaken wherein additional data, particularly that for the effect of polyslip on the activation volume for intersection, was also determined in order to provide more complete information on the details of strain hardening. Whereas analyses of these data reveal that several recommended models for strain hardening are at variance with the facts, good correlation of the data is obtained with the concept that attractive junctions4 play an important role in accounting for the rate of strain hardening. I) EXPERIMENTAL APPROACH seegerl demonstrated that slip in fcc crystals at low temperatures is dependent on thermally activated intersection of glide dislocations with forest dislocations. This has been confirmed by tests on single crystals of aluminum by Mitra, Osborne and Dorn5 and on polycrystalline aluminum by Mitra and Darn.&apos; Thus in accord with Seeger&apos;s theory, the shear strain rate, ?, below a critical temperature, T, is where .V is the number of points of contact per unit volume between forest dislocations and glide dislocations, A is the area swept out per successful intersection, b is the Burgers vector, v is the frequency of vibration of the segment of the glide dislocation undertaking intersection, U is the activation energy for intersection, k is Boltzmann&apos;s constant, and T is the absolute temperature. For aluminum, which has an extremely high staeking-fault energy, the constriction energy is negligibly small and therefore the activation energy decreases practically linearly with the stress according to as will be recomfirmed later, where Uo is jog energy at the absolute zero, G and Go are the shear moduli at the test temperature T and O°K, respectively, L is the spacing of the forest dislocations, t is the applied shear stress for slip, and tGo is the stress field that must be surmounted athermally.

Jan 1, 1965

By Evans B. Mayo

This study examines the structural framework of the Southwest for evidence of the four principal trends of lineament tectonics. It attempts to classify their intersections and compares the positions of those trends that appear most favorable with the positions of the presently known mining districts. This is a controversial topic. The author, in presenting his analysis, is aware that the study of lineament tectonics and relation of ore districts to regional structure is complicated by insufficient data and, unavoidably, by personal bias. The development of lineament tectonics has been summarized by Umbgrove.&apos; Early attempts to fit ore districts into the Cordilleran framework were made by Billingsley and Locke, who do not refer to lineament tectonics, although their appoach is similar; they observe that heat and fluids, including the ore-depositing fluids, are most likely to rise at or near intersections of major structures where the crust is fractured, or weakened, to great depth. As a result of studies distributed over the earth-including ocean basins as well as continents— some tectonists recognize four dominant structural trends: 1) northwest; 2) northeast; 3) nearly east-west, or equatorial; and 4) nearly north-south, or meridional. Baker&apos; proposed a theory to account for these trends and Sonder" called their world-wide arrangement the regmatic shear pattern. Moody and Hill proposed a much more complicated shear network which, although fascinating and perhaps ultimately useful, will not be followed here. In a recent review of deformation within the Cordillera, Wisser mentioned the four fundamental directions. The fact that many geologists deny the existence of the regmatic shear pattern implies that the fundamental structures are far from obvious. It may mean, also, that some geologists are not accustomed to examine regional and world maps analytically. The maps require much study, and certain features should be isolated on overlays. Even so, with the present limited knowledge, uncertainties remain. The following analysis is a qualitative experiment, subject to change as information accumulates, and should be supplemented by the western sheet of the Tectonic Map of the United States.&apos; Many have probably gained the impression that the Cordillera of the West is oriented northwest-southeast and from this, unless experience rules otherwise, it is natural to assume that the structure likewise trends northwest-southeast. To an important extent this is true, yet anyone tracing off the western sheet of the Tectonic Map all the recorded north west-southeast structures may be surprised to see what a small part of the entire area they occupy. In southwestern U. S. (Fig. 1) the northwest-southeast structures are mostly restricted to the eastern, southern, and western parts. The eastern margin of the Cordillera from northern Colorado to the International Boundary near El Paso, Tex., is obviously determined by some structure other than northwesterly ones. In eastern Nevada and far southward toward the Gulf of California many mountain ranges, valleys, and faults are meridional. In the Rocky Mts. of Colorado, and at many places in southern Arizona, the crystalline Pre-Cambrian is foliated northeast-southwest. The Uinta Mts. of Utah trend approximately east-west, in much the same way as a broad belt of transverse, west-northwest structures—the Texas lineament of Hill",&apos; and Ransome10 in southern Arizona, southwestern New Mexico, and southern California. It seems, then that there are four regmatic shear directions in the Southwest, but at many places they are discontinuous, and their projections must be inferred. To clarify these trends the four sets have been isolated into two systems: 1) northwest-northeast and 2) east-west-north-south (Figs. 2 and 3). Northwest-Northeast System: A number of prominent northwest-trending zones of structure are easily recognized. They are designated by circled Roman numerals, the northeast-trending structures by circled capital letters. Perhaps no two geologists would agree completely on the positions of all these belts. Names given below are for convenience only, and may be discarded where other names have priority. (I) The Sierra Nevada-Lower California belt contains the Jurassic-Cretaceous granitic massifs of the Sierra Nevada and Lower California. These

Jan 1, 1959

By B. H. Caudle, L. H. Silberberg

Reservoir depletion by natural water drive is typified by the movement of water from an aquifer into the adjacent oil-bearing formation. Prior studies of this type 01 water movement have generally neglected the resistance to flow in the aquifer and its effect on the movement of water into the oil bearing zone. A method for designing and operating scaled models of such reservoir systems is presented. Experimental data on a model of an edge-water-drive reservoir are shown and discussed. INTRODUCTION Most petroleum reservoirs derive at least a part of their productive capacity from water influx. This water may be injected from the surface or it may come from an aquifer adjacent to the oil zone. In either case, the reservoir engineer must be able to estimate the advance of the water as a function of either the elapsed time or the fluid produced. This paper describes a type of fluid flow model which includes the effect of viscous and gravitational forces in the reservoir and the surrounding aquifer. Natural water influx can be divided into three general (but widely overlapping) categories according to the direction of flow in the aquifer. These are shown schematically in Fig. 1. Fig. l(a) illustrates an edge-water-drive mechanism. In this case, water advances updip along the stratum but the movement of water is essentially horizontal and very little of the oil is actually underlain by water. The bottom-water drive (Fig. l(b)) is characterized by a thick aquifer underlying the oil zone. The water movement is generally vertical in the aquifer. The third category is illustrated in l(c). Referred to as the "thin oil column", this type of oil accumulation consists of oil over water in a relatively thin stratum. Fluid movement is horizontal in both the oil and water zones, except close to the producing well. One characteristic which all three of these types have in common is that part of the water flow takes place in the water invaded region of the oil reservoir while the remainder of the water flow occurs in the aquifer. The resistance to the flow of water usually will not be the same in these two regions. This occurs because the porous rock normally contains only water in the aquifer, while the water-invaded region always contains microscopic globules of bypassed oil (residual oil) which interfere with the flow of water through the rock. The production history of any type of a water-drive system is a function of two phenomena — where the water goes, and how it displaces the oil in the area invaded by the water. It is usually impractical to study these two phenomena at the same time. In most cases one type of analysis is used to predict the gross movement of water in the reservoir, while a different type of analysis is employed to determine the amount of oil to be recovered from the gross volume contacted by the water. The portion of the oil reservoir which is invaded by water is mainly a function of the resistances to fluid flow in the several parts of the flow system. This gross water influx in the three types of natural water drive typically results in water cusping to the production wells either from the side or from below (water coning). In 1747, Muskat described a greatly simplified mathematical model for predicting the area invaded by water in a bottom-water-drive reservoir.1 This model used the Laplace equation and suitable boundary conditions to describe the isopotentials and streamlines in the flow system. The use of this model assumes that (1) the reservoir rock is homogeneous in nature, (2) the oil and water mobilities are equal, (3) there is no oil flowing in the water invaded region, and (4) external forces — such as gravity — do not affect the flow. This method has been used for both natural water drives and water injection projects. In this, as in most models, the method of images is used to reduce the size of the model necessary to describe the reservoir. The assumption of equal viscosities and permea-

Jan 1, 1966

By K. F. Veith, R. L. Marovelli, T. S. Chen

An analytical study is made of thermal stress distribution in a thin circular disc subjected to a peripheral thermal shock at various rates of heat transfer. The problem is of importance in predicting the thermal shock response of a rock body of finite size. The theoretical analysis is based on radial heat flow by conduction in the disc and heat exchange by convection between the disc and the surroundings. The case of constant properties and plane stress is treated. Solutions of the stress distribution are presented for both cooling and heating shocks and an average stress theory is formulated. Preliminary experimental verification was obtained from the results of shock tests on thin rock discs insulated on both flat end faces so that heat exchange took place through the exposed peripheral surface. Physical properties vital to the analysis are Young&apos;s modulus, tensile strength, coefficient of linear thermal expansion, thermal diffusivity and thermal conductivity. Plots of these properties are presented. Historically, thermal energy has been used for primary rock removal throughout the ages. Although it is not widely known, the fire-setting method of thermal rock fragmentation was still in use 100 years ago in systematic underground mining operations in Scandinavia.&apos; There, the availability of cheap labor and abundant wood fuel made the method competitive with the early blasting powders. By 1880, as high explosives replaced powder, use of the method declined and it was finally abandoned about 1885. Only 40 years ago, some experimental fire-setting methods were used for rock removal underground at Pribram and Zinnwald in Central Europe.2 Compressed air-oil burners replaced the earlier wood fuels. About 15 years ago, gaseous oxygen-oil burners were introduced in the United States and gained acceptance for difficult blasthole drilling and some quarry opera- tions. At present, there is renewed interest in the compressed air-oil burner designs. The literature included information on thermal secondary processes for rock weakening or particle liberation3,4,5,6,7,8 on thermal spallation of rock 9,10,11,12,13 and on the thermal fracture of materials other than rock. The 1955 Symposium on Thermal Fracture14 and the later work of Hassel-man, cover recent developments in thermal shock investigations on brittle ceramics. Although the number of published analytical and experimental investigations conducted on ceramic and refractory materials is large, the literature reveals that there is a lack of information on both the theoretical and actual response of rock to thermal shock. Many of the thermal shock studies on ceramic and refractory bodies are based on the case of heat transfer at the solid&apos;s surface by convection in either liquid or air. A similar heat transfer situation was adopted for a theoretical analysis of the thermal stresses in a thin circular disc subjected to a peripheral thermal shock. Experimental thermal shock tests were then conducted on rock specimens to check the analytical results. Specially insulated quartzite, basalt and taconite discs were shocked in either liquid or still air bath at various rates of heat transfer. The work was done in order to resolve the question as to whether theoretically determined tensile stresses are actually realized in a finite rock body undergoing thermal shock. Specifically, if theory predicts that the tensile strength of a rock should be exceeded in a particular type of controlled thermal shock, will fracture occur? This information is needed by Bureau researchers working on thermal and electrical rock fragmentation methods. The object of this paper is to present the main features of the theory and procedure being used in the authors&apos; approach to quantitative prediction of rock response to a thermal shock. THEORETICAL ANALYSIS In the analysis to follow, simple geometry of a thin circular disc is considered. The disc is assumed to be isotropic, homogeneous and perfectly elastic, and to have constant thermal and mechanical properties. The effect of radiation is neglected so that heat is transferred by conduction inside the disc and by convection with its surrounding. In addition, the end

Jan 1, 1967

By Rex B. McLellan

A simple model for dilute interslitial solid solutions is set up which enables a solubility equation for the equilibrium between the solid solution and the gaseous solute to be deduced. This solubility equation is used to estimate the partial thermody -namic functions of Ag-O and Fe-N dilute interstitial solid solutions from the measured equilibrium between these solutions and the gas phase. It is shown that the partial excess entropy of the solute atom in solution is a vibrational entropy which is due Principally to the oscillation of the interstitial atom and does not arise from electronic or elastic interaction between the solute atoms and solvent matrix. In this work measurements of the equilibrium between dilute solid interstitial solutions and a gaseous phase are used to deduce the important thermody-namic functions of the solid solution with the aid of a theoretical solubility equation formulated from a knowledge of the thermodynamic functions of the gas and a simple model for the solid solution. The quantity of principal interest which is deduced from the equilibrium data is the change in vibrational entropy occurring when a solute atom (u) is placed in an interstial site in the solvent (u) lattice. This paper forms part of a more general study in which the properties of condensed phases are deduced from measurements of the equilibrium between the condensed phase and a gas. Previous investigations have dealt with the equilibrium between krypton and liquid metals,&apos; a pure metal and its vapor,&apos; and dilute substitutional alloys and the vapor,3 and the subject has been popularly reviewed by McLellan.4 In the latter investigation3 the vapor pressures of very dilute Cu-Ag and Cu-Au alloys were measured with a Knudsen cell technique over a large range of temperatures and the analysis of the vapor-pressure data showed that the relative partial excess entropies were large and positive [ds= -Sudk = 30 for Au-Ag alloys and 150 for Cu-Au alloys; these values should be compared to unity, the value assumed for these exponentials in the quasi-chemical theory of solutions]. For many interstitial solutions experimental data is available for the equilibrium between the solution and a gaseous phase, and this can be used to deduce the vibrational entropies of the solute atoms in solution. The general method employed is to calculate the chemical potential of the solute atoms in the gas and equate this to pi, the chemical potential of solute atoms in the solid, which is deduced from a simple model set up for the dilute interstitial solution (in which the solute atoms act as bound oscillators in interstitial sites). This enables a theoretical solubility equation to be set up. In the case of Ag-O and Fe-N solid solutions the solubility c, is found to vary with temperature according to a relation of the type, where PU2 is the pressure of the O2 or Np molecules. Essentially, the theoretical solubility equation enables the parameter A(T), which gives the vibration entropy, and the parameter AH, which gives the energy of solution of a solute atom, to be estimated from the measured variation of solubility with temperature. The solutions treated are dilute enough to enable solute-solute interactions to be neglected so that the partial energy and excess entropy of solution E, and GS are not functions of composition and the configurational entropy is that of a random solution. Large positive values have been found for sxSu this quantity has not been given sufficient consideration in recent work on interstitial solutions and it is felt that the partial thermodynamic quantities of very dilute solutions should be understood before the thermodynamics of concentrated solutions can be fully appreciated. EXPERIMENTAL DATA 1) The Silver-Oxygen System. The measurements of Steacie and Johnson,= using silver foils, indicate that there is a pronounced minimum in the solubility vs temperature curve for O-Ag solutions. Eichenauer and Miiller6 have recently redetermined the solubility of oxygen in silver using a method developed in the investigation of metal-hydrogen systems.7 The solubility determinations were carried out on compact specimens of differing shapes and sizes and the results indicate a much lower solubility than was obtained by Steacie and Johnson with no minimum at 400°C. Subsidiary experiments showed that the erroneous results of Steacie and

Jan 1, 1964

By J. B. Cohen, O. F. Kimball

Aging of specimens of Au-20 at. pct Ni and Au-40 at. pct Ni has been followed with electron microscopy, resistance, magnetic measurements, and X-ray diffraction. The periodic structures observed by electron microscopy in Au-Ni alloys by Fisher and Embury are shown to be due to nickel-rich prepre-cipitates. Assuming this alloy system undergoes spinodal decomposition, our observations and available data on thermodynamic properties and elastic constants indicate that the gradient energy coefficient in this system is of the order of 10-5 ergs cm-1 . Shifts of X-ray peaks in this system have been found to be due primarily to residual stresses, not faulting. A large effect of applied stress during aging was detected by X-rays and electron microscopy for the Au-20 at. pct Ni alloy. The presence of a preprecipitation stage during aging of gold-rich Au-Ni alloys for short times below 300°C has been indicated by several previous investigations. Changes in electrical resistivity,1"3 elastic modulus,4 superparamagnetisrn,5 positions of X-ray diffraction peaks,6 and hardness7 have been studied. Electron diffraction and microscopic studies by Fukano8 and Fisher and Embury9 respectively showed that the precipitation process resulted in the formation of periodic composition fluctuations—perhaps spinodal decomposition. MOSS10 detected a satellite with X-rays in quenched specimens, but the wavelength was much smallerthan that found by Fukano. Recently Woodilla and Averbac11 detected a modulation similar to that found by Moss, but with electron diffraction. They suggested that the microstructures observed by Fisher and Embury may be simply due to a Moir&apos;e effect from a back deposit of gold during electropolishing. Satellites in electron diffraction can be enhanced by small bends or tilts of a foil. Even if this were not the case, asymmetric satellites, as are found in Au-Ni are difficult to analyze. They are principally influenced by the large differences in atomic size in this system, and obtaining information on the associated composition fluctuations depends on assumptions as to whether or not there are isolated clusters or composition "waves", and in the latter case whether or not the entire specimen has decomposed. Assuming that decomposition was complete in alloys aged at low temperatures, Woodilla and Averbach suggested that the satellites indicated a composition which varied about the average value by only± 5 pct, much less than that predicted from the phase limits, Fig. 1; they attributed this small variation to a cessation of the reaction caused by the loss of mobile vacancies. In Ref. 12 it was shown that vacancies are trapped by the preprecipitate in this system. Furthermore, after reversion, aging proceeds at a much slower rate than immediately after a quench from the solution temperature,5&apos;11,13 indicating the importance of vacancies in the kinetics. In some senses Au-Ni is a favorable system for studying preprecipitation and reversion mechanisms because of the simplicity of the equilibrium diagram14 and the availability of good thermodynamic15 and diffusion data.16,l7 The calculations necessary to relate the decomposition process to classical nucleation and growth or to spinodal decomposition, as well as those necessary for understanding reversion mechanisms, can, therefore, be made. As an example, the position of the chemical spinodal, calculated from the thermo-dynamic data of Sellars and Maak15 with the method described by Rundman and Hilliard,18 is included in Fig. 1. It is compared to the chemical spinodal calculated from the equilibrium miscibility gap by Cook and Hilliard&apos;s method.l9 Included is the coherent spinodal suggested by Golding and MOSS," using measured elastic constants for their calculations and the same ther-modynamic data. The spinodal estimated by Woodilla and Averbach,11 by determining the temperature limits for the observation of the satellites discussed above, is also given. Actually this last curve is open to question on the following grounds: As a phase boundary is approached the spacing of particles or the wavelength

Jan 1, 1970

By F. Bonaccorso, G. Venturello, C. Antonione

A study of the effect of small amounts of interstitial impurities on recovery and re crystallization in high-purity iron (99.995 pct) has been undertaken. This paper gives results on the effect of carbon, introduced in small-dosed amounts @om 0.0005 to 0.0086 wt pct) by heating iron specimens of high purity m a static atmosphere of CO + ,. The materials prepared in this way, cold-worked 80 pct and subjected to a series of isochronal and isothermal annealings, were submitted to examinations by X-rays, micrographs, and hardness tests. It was observed that effect of carbon is remarkable in the sense of blocking recovery of mechanical properties up to the temperature at which re crystallization begins. On the contrary, carbon has a negligible effect on primary re crystallization temperature, when compared with the known effect of substitutional impurities. This is in agreement with the high mobility of interstitials. In effect, only a slight decrease, -20 pct, of the grain-boundary motion rate was noted, due to the interaction between the grain boundary and the carbon atoms. On the other hand, in the samples in which the carbon content is above the solubility limit at the temperature at which re crystallization occurs, a slight increase of nucleation frequency is noted due to the presence of precipitated carbides. ThE effect of purity on recovery and recrystalli-zation phenomena has been known for a long time; however, recently, new attention has been given to the problem since new methods for obtaining metals with extremely low impurity contents have become available. Most of this research work essentially concerns the effects of impurities which cause precipitates or give rise to substitutional solid solutions. The works of Bolling and winegard1 and Aust and utter&apos; on lead, and Vandermeer and orddon&apos; on aluminum, 01sen4 on nickel, and Abrahamson and Blakeney on iron should in particular be referred to. On the effect of this type of impurities a quantitative theory has been formulated by Detert and Lucke6 and has lately been discussed critically by Cahn7 and Gordon and vandermeer.&apos; Interstitial solid solutions in iron have not as yet been studied. As the study of the effect of interstitials is of great interest both from a theoretical and practical standpoint, it was deemed useful to examine the effects of carbon, nitrogen, and boron on recovery and primary recrystallization of iron. There already is some work by Chaudron et al.,1-" on the effect of interstitials on iron; their work, however, mainly concerns secondary recrystallization. The present paper refers in particular to the effect of carbon. EXPERIMENTAL PART Preparation of Materials. The pure iron used for this research was obtained from FeCIS recrystal-lized and purified by extraction with isopropylic ether. From the ferric chloride purified in this manner, hydroxide was precipitated with a solution of very pure ammonia, and, by calcination in pure sintered alumina crucibles, oxide was obtained. Reduction of the oxide to iron sponge was performed in sintered alumina tubes with very pure hydrogen at 650°C; at the end of the operation temperature was increased to 900°C. Specimens for the experiments were obtained by sintering the sponge at 1480°C in pure Hz after a pre sinter ing treatment at 900° C. It is important to note that the treatment at 1480°C in Ha produces a further purification from more volatile elements such as zinc, cadmium, arsenic, lead, and tin. Details on the preparation and characteristics of this type of very pure iron are given in a previous work.&apos;&apos; Only the complete analysis performed by neutron-activation methodsz3 is given here, Table I. Some of the specimens prepared in this way were carburized in the 0 region with very low amounts of carbon by treating them at 700°C in a static atmosphere of Ha containing a definite amount of CO. The set-up used is described in Fig. 1. A gas-tight quartz tube containing the specimen to be carburized and an internal friction control specimen, after being evacuated, was filled with Hz

Jan 1, 1963

By C. S. Matthews, P. Hazebroek, H. Rainbow

It ha been suggested that lormation fractures created by well stimulation treatments will adversely affect sweep-out efficrency in injection operations. Fluid-flow model studies involving vertical fractures of various lengths and fluid systems of various mobility ratios have been carried out to study this subject. In addition, limited data have been obtained on one model containing a horizontal fracture. It was found that relatively long and highly conductive fractures (not generally obtained fracturing operations) were required to affect the sweep-out efficiency substantially. In a given case in the field an approximate distinction can be made between the presence of long conductive fractures and shorter or less conductive ones. This is done with pressure build-up analyses along with data on the relationship of fracture length and conductivity to well productivity. This type analysis shows that usually the fractures induced are either short or of limited conductivity and therefore do not damage sweep-out efficiency. INTRODUCTION Improved well productivity and in-jectivity can frequently be exploited in injection operations. Higher total throughput can yield improved economics. On occasion, the achievement of an increased productivity or injectivity in specific wells can bring about a more uniform sweep of the reservoir. Higher rates can be exploited particularly in water floods of "depleted" reservoirs where a rapid "fill-up" is desired and where low pressures contribute to low well productivity. The creation of fractures local to the wellbore is an excellent means for achieving these objectives. Even though fracturing has been employed in some floods with success,1"3 there still seems to be some reluctance to employ this tool for fear of undue damage to the flood pattern and ultimate recovery. We therefore need to examine the influence which fractures of varying length may have on flood performance and then determine the length of fractures which obtain in the field with conventional fracture treatments. A substantial influence of fractures on the recovery obtained at breakthrough of the injected fluid has been presented by Crawford and Collins for equal fluid mobilities for the line-drive pattern.4,5 In addition, the influence which a fracture has on the production performance after breakthrough and on the ultimate recovery warrants consideration in reaching a conclusion concerning the use of induced fractures in flooding operations. This report presents this type of data for the five-spot injection pattern for several fluid mobilities. In arriving at some conclusion concerning the fracture lengths obtained in field operations we must examine the performance characteristic most affected by the fracture. This is the change in the flow system as reflected in the change in productivity and pressure build-up behavior. A study of these changes is also presented in this report. RESERVOIR ANALOGS The X-ray shadowgraph technique, employing miscible displacement in porous models, has been used in the study of the influence of fractures on pattern sweep-out efficiency. The X-ray shadowgraph procedure is described in detail in an earlier report8. Fractures were represented by leaving the proper portion of the model surface exposed to either injection or production. This assumes the fracture resistance to be negligible compared to that of the formation. Actually the flow resistance in propped fractures obtained in the field is sometimes not negligible so that the results with this model indicate the maximum influence of the fracture. Two types of models were necessary to represent vertical fractures in a five-spot flood. These pattern elements are illustrated in Fig. 1. In studying the influence of the horizontal fracture, only one well spacing length to thickness ratio of 50 was used. The pool unit of a Carter Electric Analyzer was used in studying the influence of fractures on productivity and build-up behavior. The square drainage system of a well was represented by a network of 576 elements of equal volume and resistance. One-fourth of this drainage area was studied as a model unit using a network of 144 condensers and resistors. Vertical fractures were represented by a shunt directed from the well perpendicular to the drainage boundary. Horizontal fractures were represented by a circular shunt. Resistance of a shunt was varied in representing different conductivities for the fractures. FRACTURE DIRECTION AND LENGTH The direction at which a vertical fracture extends into the formation and its length and conductivity influence the sweep behavior. The two

By S. S. Marsden, S. A. Khan

Externally generated foam was injected continuously into short porous media. Both flow rate and pressure drop were measured. Liquid saturation was determined by electrical conductivity. Foam yuality I, expressed as the ratio of gas volume to total volume, was varied from 0.70 to 0.96. As measured with a modified fann VG Meter, apparent viscosity of this foam µa decreases with increasing shear rate but usually falls within the range of 50 to 500 cp. At a given shear rate, µa increases almost linearly with I-. When measured with a Bendix Ultraviscoson, kinematic pa is independent of V but absolute pa increases with T from about 3 to 8 cp. The effective permeability -apparent viscosity ratio ke/µa decreases almost linearly with V for porous media of high permeability, but the rate of decrease becomes less for tighter ones. The relative permeability-apparent viscosity ratio kr/µa us V data does not fall on a single line. The kr/µa ratio increases with liquid saturation in the porous medium and with surfactant concentration. Estimates of pa for foam in porous media vary from 30 to 100 cp. INTRODUCTION Although research on the development of a foam-drive, oil recovery process has been going on for almost a decade, most of the significant publications have appeared within the last several years. This illustrates well the rate at which interest in this process is accelerating. Bond and Holbrook 1 were the first to describe the use of foam to improve oil recovery in their patent of 1958. They proposed that an aqueous foaming agent slug be injected into the formation and that this be followed by gas to produce a foam in situ. Fried 2 studied the injection of foam into porous media which had already been subjected to conventional gas or water drives and found that gas could be used to drive a foam bank which would, in turn, displace additional oil in the form of an oil bank. He attributed the increased oil recovery to the high effective viscosity of foam flowing in porous media. His microscopic observations showed the importance of foam generation and regeneration within the porous medium. By injecting both air and aqueous surfactant solution, Bernard3 generated foams within the porous medium in which oil displacement was being studied. In a separate empirical test, he also measured the dynamic foaming characteristics of the same surfactants in water and/or oil. With some exceptions and for the seyen surfactants studied, there seems to be a qualitative relationship between the efficiency of liquid rlisplacrment and the dynamic foaming testtsed. This relationship was not consistent enough to eliminate the necessity of actual foam flood tests in porous media for surfactant selection. In a study basic to gas storage in aquifers, Bennett4 described the displacement of brine by foam in consolidated porous media. Among other things, he stated that the ability of a surfactant solution to foam is more important than the stability of its foam. The presence of a foam bank between the displacing air and the displaced brine improved both breakthrough and ultimate recovery. In a continuation of this work Kolb5 attributed the great reduction in surfactant solution production rate as displacement by air progressed to a decrease in relative permeability to gas. These several effects reported by both Bennett and Kolb can all be attributed to the high apparent viscosity of foam which was obviously flowing in the porous media. Bernard and Holm6 observed that gas permeability in the presence of foam was much less than in the absence of foam when both were measured at the same gas saturations. This reduction in permeability was proportionately greater for more permeable sands and sandstones than for the less permeable. Viscosities of the surfactant solution and air were used for their calculations rather than apparent foam viscosities. The foam stability increased with decreasing permeability, which was probably due to the natural regenerative properties of porous media, and the permeability reduction increased

Jan 1, 1967

By S. K. Tarby, A. E. Rathke

A series of 0.21 pci carbon steel melts was processed under conditions which sinzulated industrial vacuum degassing practices. The results indicated that the efliciency of carbon deoxidation was not significantly increased by vacuunz degassing at chamber pressures below 100 torr. The two common phenomena of vacuu?rz carbon deoxidation, i.e.,an unaccountable carbon loss and nonequilibrium deoxidation, were observed in this study. Oxygen from adsorbed and chenzically combined water of certain refractory structures was probably responsible for the unaccourttable carbon drop. The nonequzlibrium deoxida-tion phenomenon was adequalely correlated by a bubble nucleation )mechanism. Although the major limiting factor in oxygen removal from these melts contained in alumina crucibles could not be eslablished with certainty, there was some evidence that melt-refractory interactions were not the most imp-portant limiting factor. CONSIDERABLE attention has been given to the use of carbon in conjunction with vacuum degassing* as a deoxidation technique. An advantage of utilizing the reaction: for deoxidation purposes is that the reaction product is gaseous and hence can be effectively removed from the melt. The pressure dependency of Eq. [I] is illustrated by the expression for the thermodynamic equilibrium constant, K: where pco is the partial pressure of carbon monoxide, and hc and ho are, respectively, the Henrian activities of carbon and oxygen in molten steel. The thermodynamic aspects of the C-0 reaction in molten Fe-C-0 alloys have been investigated extensively. Theoretical predictions based on Eq. [2] indicate that, at the chamber pressures reached in most vacuum degassing units, carbon is the most effective deoxidizer that can be obtained. ~ornak&apos; has compared industrial values of carbon and oxygen contents with values calculated from theory, and has shown that on a commercial basis vacuum carbon deoxidation has not achieved the degree of success that is predicted from theoretical considerations of the C-0 reaction. The use of industrial units to determine the reason for this deviation is limited by economics and the difficulty in controlling critical process variables. Therefore laboratory-sized experiments in this area should yield valuable information. Investigations on vacuum carbon deoxidation have been made in experiments conducted under vacuum induction melting conditions where pressures in the low micron range, long treatment times, exposure to highly dynamic vacuum conditions, and relatively quiescent bath conditions have been maintained. Inve~ti~ators~~~ utilizing vacuum induction melting for the production of high-purity iron by carbon deoxidation have concluded that, under high vacuum conditions, melt-refractory interactions constantly supply oxygen to the melt. Thus, it has been suggested that refractory dissociation rather than the C-0 reaction establishes the lowest possible oxygen content of a melt. Thomas and Moreau4 and Bennett, Protheroe, and ward5 have studied carbon deoxidation of Fe-C alloys melted in magnesia crucibles under reduced pressures. The results of both studies differed considerably from those predicted by theory, and the limiting factor for lower oxygen contents was attributed to melt-refractory interactions. Because the refractory material used in both investigations was magnesia, no quantitative measure of crucible reaction could be made and the amount of carbon consumed during treatment was utilized as an indirect indication of the extent of crucible reaction. Studies concerned with the determination of the minimum oxygen content attainable in pure iron and Fe-C alloy material vacuum-induction-melted in various crucible materials have also been conducted."&apos; The suitability of a crucible refractory for use under vacuum conditions appears to depend upon the major refractory constituent, the amount and type of impurity oxides, and the composition of the melt. Meadowcroft and Elliott9 have presented a theoretical review of the type and extent of various refractory reactions which should be considered in vacuum refining. ~eadowcroft&apos;~ considered the vacuum degassing process to be dependent upon the "carbon boil" and has proposed that a part of the difference between the experimental and theoretical results for vacuum carbon deoxidation may be explained by the fact that the pco term in Eq. [2] is not the chamber pressure. Instead, the pressure necessary for carbon monoxide

Jan 1, 1970

By S. C. Jones, J. A. Davis

This paper describes displacement mechanisms of micel-lur solution slugs, displuced by a thickened water "mobility buffer", in a glass micromodel and in consolidated Berea sandstone cores. Colored motion pictures, taken through a microscope, recorded interactions on a pore-to-pore bash in the micromodel. Production and differential pressure data from consolidated cores previously waterflooded to residua1 oil saturation revealed displacement mechani.sm.s macroscopically. Oil is displaced miscibly, while interstitia1 water is displaced immiscibly by the injected micellar. solution slug. Part of this water is dispersed and carried near the lead ing edge of the .slug, and some may remain as a residual saturation. Thickened water displaces any residual water. It also emulsifies with the. slug material, thereby displacing it. The displacement sequence forms banks of oil and water with relatively constant fractional flows, saturations, and mobilities. The mobility of oil and water flowing ahead of the .slug defines .slug viscosity required for stable displacement. Equations are presented that describe. several aspecis of the displacement flow behavior: oil breakthrough, oil cut during production of the stabilized bank, breakthrough of banked water and the resulting decreased fractional flow of oil, fluid mobilities, etc. Breakthroughs, calculated from the equations, are compared with experimental data. Slugs of micellar sollution, ranging from 1 to 5 percent PV, produced unusually high oil recoveries (up to 100 percent) in previously waterflooded consolidated cores. Investigation of the effect of mobility buffer size on oil recovery indicated that recovery is reduced in 3-in. X 4-ft cores for mobility buffer sizes .smaller than 50 percent of the pore volume. Introduction A new oil recovery process,* conceived in 1962, utilizes a slug of micellar solution&apos; that has sufficient viscosity to avoid an unfavorable mobility ratio with respect to the oil and water it displaces. These micellar solutions are composed of hydrocarbons, water, and surfactants, and may contain electrolytes and co-surfactants. The micellar solution is displaced by a "mobility buffer". The mobility buffer is a bank of reduced-mobility fluid, preferably an aqueous polymer solution. The entire system is displaced by unthickened water. This paper describes some of the microscopic and macroscopic displacement mechanisms observed during laboratory development of the fluids for a potential field test. Microscopic phenomena were observed photographically in a specially constructed glass flow model and macroscopic data were taken from constant rate, tertiary displacements conducted in Berea sandstone cores. Qualitative Description of Micellar Solution Flooding This section qualitatively describes displacement mechanisms that occur when a slug of micellar solution is injected into a previously waterflooded core. The slug completely removes oil from the portion of core it contacts, and banks oil in the process. The oil saturation in this bank and the flowing WOR reach definite, fixed values that are determined by water and oil mobilities and slug characteristics. This zone of constant WOR and oil saturation will be called "stabilized oil bank" in the following discussion. When this stabilized bank reaches the outlet end of the core, the oil cut jumps from zero to a constant value equal to the fractional flow of oil. Saline water has little or no solubility in the slug. Because of this. the slug displaces in-place brine immiscibly. During slug injection, most of the brine is displaced ahead of the slug, some is carried along in the leading portion of the slug, and part of the brine may become immobile. The immobile, residual brine is later miscibly displaced by thickened water. Displacement of all the oil and only part of the water during slug injection results in a higher fractional flow of oil ahead of the slug than occurs when all the water is displaced after injection of thickened water is begun. The slug spontaneously emulsifies fresh (low-salt content) water. By emulsification, the slug can become infinitely diluted, and hence totally displaced, by thickened water. Fluid saturations in the core are shown schematically in Fig. 1. Parts A and B show conditions before and after the start of the thickened-water injection, respectively. Since thickened water displaces any residual brine by-

Jan 1, 1969

By R. D. West

Miscible displacenzent recovers all oil in the area contacted by the injected .fluid, whereas water or immiscible gas drives usually leave substantial amounts of oil as residual. However, the Door mobility ratios associated with a gas-driven miscible displacement cause the sweep pattern efficiency to be much lower than that obtained with water flooding. One way in which the sweep eficiency in a miscible displacement process can be increased is by decreasing the mobility behind the flooding front. This can be achieved by injecting water along with the gas which drives the miscible slug. This water reduces the relative permeability to gas in this area and thus lowers the total mobility. The main operating conditions for the simultaneous injection -vrocess are that a zone of gas exists between the miscible slug and the leading edge of the water and that a su,@cient amount of gas be injected with the water to form the pas volume which is being left in the water zone. Laboratory model studies have shown that the ultimate sweep pattern efficiency can be as high as 90 per cent for a five-spot flooding system. If gar alone is used as the driving medium an ultimate sweep-out efficiency of about 60 per cent would be obtained in the same system. I INTRODUCTION The miscible displacement processes are a step towards total oil recovery. Conventional gas or water drives usually leave 25 to 5.0 per cent of the oil as residual in the swept portion of the reservoir. This residual can be eliminated if the oil is driven by a fluid with which it is miscible. At some reservoir conditions natural gas will become miscible with the oil. This is the "high pressure gas process".&apos; More often, the oil does not contain enough light hydrocarbons to cause the gas to become miscible with the oil at reasonable pressures. In these cases a small band of fluid which is miscible both with the oil and gas must be kept between them2. Less than 2 per cent of the reservoir volume of the slug material is needed to keep the displacement miscible. Both processes work in the same manner, recovering all of the oil in the portion of the reservoir contacted by the injected fluids. The only difference is the manner in which the miscibility between the oil and the injected gas is obtained. Previous publications have contained detailed descriptions of these processes.1,2,3,4 However, total displacement of the oil in the swept region does not guarantee an efficient recovery process. The amount of oil to be recovered is also determined by the fraction of the reservoir contacted by the flood. This fraction is largely determined by the mobilities of the fluids. (The fluid mobility is the permeability of the rock to that fluid divided by the fluid&apos;s viscosity, k/p). This. dependence of the fraction swept on the mobility ratio has been shown in previous studies. Fig. 1 shows the ultimate fraction swept in a five-spot system as a function of the mobility ratio. The small drawings show the location of the areas left unswept for two different mobility ratios. The ultimate fraction of the reservoir swept is here considered to be attained when the producing stream contains less than 5 per cent oil at reservoir conditions. THE GAS-DRIVEN MISCIBLE DISPLACEMENT Since there is no oil left in the swept region after miscible displacement the mobility in this region is very high. It is often 50 times the mobility in the unswept regions. This means that the fraction of the reservoir contacted by the injected fluid will be less for a gas-driven miscible displacement than for a conventional water or gas drive. For a five-spot injection system, water would contact the entire reservoir volume, and the low pressure gas would contact about 90 per cent of this volume, while a gas-driven miscible displacement would only contact about 65 per cent of the reservoir. This poor sweep efficiency often offsets the benefits obtained through miscible displacement. Fig. 2 shows what the recovery curves for the three processes might look like for a five-spot system. The curves show the fraction of the in-place oil recovered as a function of reser-

By N. Brown, R. Kossowsky

Plain carbon steels with 0.48 and 0.95 pct C were quenched and tempered at 705°C to produce carbide dispersions with spacings on the order of 1 p. The morphology of the structure consisted of a carbide-dislocation network. The strengthening due to the dispersion was found to vary linearly with M-½ where M is the mean free-ferrite path determined by the entire network. No preyield microstrain preceded the upper yield point. After prestraining, the lowest stress at which dislocation movement could be detected was 104 psi; this frictional stress was independent of the dispersion and prestrains up to 6.5 pct. The Cottrell-Petch equation for grain-size strengthening was used to discuss the dispersion strengthening in this investigation. The results support an impurity mechanism for the upper yield point rather than one based on the Johns ton- Gilman theory. IT was pointed out by Gensamerl that the mean free path in the matrix is the important variable which controls the degree of dispersion strengthening. Gensamer&apos;s data on steels showed a linear relationship between the yield point and the logarithm of the mean free-ferrite path. The first theory was by Orowan,2 who suggested the mechanism of dislocations bowing between particles, with the resulting relationship that where a is the yield point, P is the distance between particles, and G is the shear modulus. The Orowan equation applies to that range of dispersion where P is large compared to the particle size and the particles are not coherent, so that the matrix is essentially free of internal stresses. As was pointed out by Orowan, there are different degrees of dispersion which, in turn, will influence the strength-controlling mechanism. However, in this investigation, we wish to confine ourselves to the region of coarse dispersion, where the Orowan mechanism should, intuitively, be applicable. It was soon evident that the data, which existed at about the time that the Orowan theory was proposed, did not agree with the Orowan equation in that the actual strengthening was always greater than the theoretical predictions. Thus, Fisher, Hart, and pry3 modified the Orowan theory by suggesting that the Orowan process took place in the microstrain region preceding macroscopic yielding and the subsequent, rapid work hardening in the form of residual dislocation loops around the particles largely determined the observed macroscopic yield point. The F-H-P theory stated that the increment of strengthening due to the work-hardening mechanism was proportional to f3/2 where f is the volume fraction of the precipitate. F-H-P used data by Shaw, Shepard, Starr, and Dorn4 on A1 3-5 pct Cu to support their theory. Roberts, Carruthers, and Averbach5 were the first to make a microstrain study of dispersion strengthening, and they found that for steel the Gensamer relationship was obeyed. Hayman and Nutting6 suggested that in the case of tempered steel 1) the ferrite grain boundary was the primary obstacle and 2) the carbide particles simply formed part of the grain boundary obstacles. They found that the strength varied as G"" where G is the ferrite grain size. Turkalo and LOW&apos; determined the structure of quenched and tempered plain carbon steel using a replica technique; they concluded that the carbide particles did not necessarily lie in the ferrite grain boundaries. When they defined the mean free-ferrite path as being determined by both particles and the ferrite path boundaries, their data obeyed the Gensamer relationship. Meiklejohn and skoda8 showed that the strengthening from iron particles in mercury was a function of the particle size and the distance between particles. Dew-Hughes9 explained the Meiklejohn and Skoda results by the following theory: 1) grown-in dislocations which surrounded the particles were produced by the thermal stresses during the time the material was cooled to the testing temperature, and 2) the observed strengthening was associated with the cutting of the grown-in dislocations by the glide dislocations. Ansell and Lenel10 proposed a theory in which the glide dislocations pile-up or surround the particle until the number of dislocations in the pile-up is sufficient to yield or fracture the particle. The resulting theory says that the yield point varies as p-1/2, where P is the distance between particles. The Ansell-Lenel theory is almost identical to the

Jan 1, 1965

By C. W. Arnold, H. L. Stone, D. L. Luffel

The purpose of this research is to determine (I) the efficiency of small banks of enriched gar driven by methane in displacing oil from a porous medium and (2) the effects of variation in bank size and composition of that efficiency. Most of the experiments were conducted in a sand-packed tube 20-ft long and 1/2-in. in diameter. The hydrocarbon system generally used was methane, butane and decane at 2,500 psia and 160°F. The results of these experiments indicate that, in the regions contacted by the gas, a small bank of an oil-miscible gas driven by methane can displace all of the oil in a piston-like manner. If the enriched gas is of such composition as to remain immiscible with the oil, displacement of oil is less efficient than for the miscible case, and the gas bank travels through the sand with a velocity less than that of the driving gas. These data along with theories discussed imply that smaller banks and less total gas are required when the enriched gas and oil are miscible. INTRODUCTION Widespread application of enriched-gas drive to the recovery of oil rests upon a key factor — the use of limited quantities, or "banks", of enriched gas. At the present time, the value of liquefied petroleum gas or other enriching agents discourages their use in a continuous injection technique, or even in a large bank, except in a few isolated reservoirs. If small banks of enriched gas driven by methane were as effective in displacing oil as is continuous injection, the enriched-gas drive process might be applied to a larger number of reservoirs. Previous research on the mechanics of the enriched-gas drive process reported by Stone and Crurnpl and by Kehn, Pyndus and Gaskell has utilized continuous injection of enriched gas. This work has shown that two types of displacements occur. With gases containing sufficient intermediates. the oil is displaced misciblv and complete recovery is obtained from the regions swept. When gases are used which contain insufficient intermediate hydrocarbon for miscible displacement, oil is displaced immiscibly. In the latter type, selective solution of the intermediate hydrocarbons causes a swelling and reduction in viscosity of the oil and leads to an increased recovery over that obtained by dry-gas (methane) drive. The size of the enriched-gas bank necessary for efficient displacement of oil is determined by those factors which cause deterioration of the bank. A differentiation may be made between those factors which operate on a microscopic scale and those which act on a macroscopic scale. On the smaller scale, the enriched gas mixes in the direction of flow by diffusion and convection with the fluids immediately preceding and following it. On the larger scale, the gas may by-pass the oil by flowing through permeable streaks, by overriding the oil because of density difference, or by fingering because of unfavorable viscosity ratios. In such cases, the enriching material tends to mix with the oil both laterally and in the direction of flow. The increase in effective area available for diffusion and dispersion of the enriching components leads to a faster degradation of the bank and a need for a larger bank than is necessary for those cases in which no by-passing occurs. The effects of such macroscopic factors in the deterioration of enriched-gas banks have been reported in a separate paper by Blackwell, Terry and Rayne. The present study was confined to the factors which operate on the smaller scale, in particular to the behavior of banks of enriched gas in sands uniformly swept by the gas. Experiments were designed to answer the following questions. 1. Can small banks of enriched gas driven by methane be used to secure oil recoveries comparable to those obtained by continuous injection of enriched gas? 2. What is the optimum bank size (the minimum bank size necessary to obtain a recovery comparable to that obtained by continuous injection of the same enriched gas)? 3. How many total pore volumes of gas must be injected to obtain the maximum recovery when the optimum bank size is used? 4. What is the effect of varying the number of enriching components in a gas bank? This report describes the experimental investigations and discusses the results in terms of their significance to reservoir behavior.

By K. I. Savage, S. C. Sun

This paper describes a process developed to recover coal, clays and pyrite from coal wastes. The process consists of fine grinding followed by coal and pyrite flotation which leaves the clays in the flotation pulp. A bituminous coal refuse containing 10% sulfur and 30% carbonaceous material was treated by this method to yield a coal product containing 4% sulfur, 10% ash; a pyrite product containing 45% sulfur (84% FeS2), 1% carbonaceous material; and a clay product containing 2% sulfur (3.5% FeS2). The coal yield was about 89%. The pyrite yield was about 77%. The process steps may be entirely flotation, or gravity separation (hydrocycloning) may be used to increase the pyrite : coal ratio in the flotation feed. Cost estimates for the process show a profit of $2.28 per ton of low pyrite grade refuse, but these do not include labor, maintenance, overhead and plant depreciation. The development of this process consisted of three parts: (1) exploratory tests, (2) op-timazation tests and (3) confirmatory tests. The objectionable qualities that sulfur imparts to coal have been commented upon from early times, and they have become more objectionable as the uses of coal have grown. In whatever form coal takes — raw, carbonized or gasified —the sulfur content remains objectionable, and therefore its compounds are removed as completely as possible. It has long been known that sulfur occurs in coal in different forms. In 1861, sulfur was said to exist in the state of sulfuric acid in combination with a base; in combination with iron as iron pyrites; as bisulfides of iron; and in combination with the organic elements of coal.&apos; Pyritic sulfur is a term loosely used to cover the sulfur associated with iron in its various forms. The mineralogically recognized forms are pyrite (FeS2), pyrrhotite (Fe 1-x S) and marcasite (FeS2). The particle sizes of pyrites vary widely. Isolated grains of marcasite smaller than 15 microns have been found disseminated through coal.2 Hair-like "veins" of pyrite filling cracks in vitrain have been found. At the other extreme, lumps or nodules of pyrites large enough for removal by hand picking have been encountered. Organic sulfur, unlike pyritic sulfur, does not exist as discrete particles, but is instead intimately associated with the coal structure and thus it is impossible to remove it or reduce its concentration by physical or mechanical means.2 In the preparation of coal for its various markets, the pyrite minerals (pyrites) are separated from the raw coal feed. This separation process concentrates the pyrites in the tailings or other waste products. It would be desirable to recover these pyrites for three reasons: (1) Pyrites are potential sources of sulfur and iron. In 1967, for the fifth consecutive year, Free World consumption of sulfur exceeded production.3 Propelled by the shortage, the domestic price of sulfur has risen from $24 to $38 per long ton (bright). (2) The refuse, when placed in piles, becomes ignited. The pyrites (FeS2) bum, giving off SO,. Thus, the pyrites are a cause of air pollution. (3) Also, the pyrites undergo chemical reaction when exposed to air. The refuse is leached by waters which result in stream pollution due to the water-soluble iron and acidic reaction products. Coal refuse also contains coal minerals and clay minerals. Therefore, any process for recovering the pyrites must successfully separate them from the coal and clay minerals. In the study discussed here, ten different bituminous coal refuse samples were successfully upgraded in pyrite content. These samples represented a wide variety of coal waste materials from Pennsylvania and other states. The variabilities of the sulfur and coal contents are shown in Table I. The extremes are Sample J, a high sulfur-low coal, and Sample B, a low sulfur-high coal. Definitions of some of the symbols or terms used in this report are given below: Mesh size—Tyler standard mesh screen sieves with an opening based on the square root of two. Fe —Indicates iron, as determined by a stannous reduction-dichromatic oxidation method. S—Indicates sulfur, as determined by the ASTM "Eschka" method for sulfur in coal.

Jan 1, 1969