Search Documents

By J. T. Raleigh, D. L. Flock

This study investigated the nature of formation plugging with bacteria and attempted to relate its characteristics to physical rock properties. Fifteen core samples of four specific formation types were defined and plugged using a commonly occurring, uniform sized dead bacteria, Bacillus subtilis. Two of the formation types were fairly uniform-grained sandstones and two were heterogeneous carbonates. The injection rate and concentration of solids in the brine were held constant during the test runs. The pore geometry factor G was a significant petro-physical rock characteristic for correlation with depth of plugging, but lacked importance as a parameter for rate of plugging. Although the G factor may not completely describe porous rock geometry for plugging predictions, for reservoirs with multi-pore systems, there is evidence that the pore geometry factor might form a basis for establishing injection water specification in more llornogeneous reservoirs with only primary porosity development. INTRODUCTION GENERAL PURPOSE Injecting water into oil-bearing formations to increase recovery is common practice in the petroleum industry. One of the most serious operational problems involved with water injection wells is the tendency of the porous rock immediately adjacent to the injection wellbore to plug partially. This results in increased injection pressures and/or a decrease in the rate of injection. This plugging in injection wells can result from a number of causes, including bacteria contained in the injected water.1 The purpose of this work was to investigate the nature of injection plugging with bacteria and to attempt to relate its characteristics to physical rock properties. Many variables exist in research of this type, so, to reduce them, injection rate and concentration of bacteria in brine were held constant. The bacteria used for plugging was a uniform-sized, dead bacteria, Bacillus subtilis, a micro-organism commonly found in water, air and soils. Four specific porous rock types were plugged—two consolidated sandstones and two heterogeneous carbonates. The experimental work involved two distinct phases. The first concerned determining the physical rock properties and the second was the experimental work of actually plugging the selected core samples. PORE GEOMETRY FACTOR Capillary pressure is necessary to displace a wetting fluid from a capillary opening with a nonwetting fluid. The factors which govern capillary pressure are the surface or interfacial tension of the fluids involved, the system wettability and the equivalent pore radius. Thus: P,= 2ucososrd..........(1) Thomeer used the idea that the location and shape of the capillary pressure curve reflect characteristics of the pore structure of any porous media sample.' He presented the mathematical description of capillary pressure curves and of indicated differences in pore geometry of samples. where G = pore geometry factor PC = capillary pressure Pi = displacement pressure (V,)p, = fractional bulk volume occupied by the displacing fluid at any capillary pressure (V,), = fractional bulk volume occupied by the displacing fluid at infinite pressure or total interconnected pore volumes. The shape of the curve is defined by the pore geometry factor G. The curve shape is related directly to the sorting and interconnection of the pores of the sample. The location of an infinite number of curves can be defined by the same asymptotes. However, these curves differ in shape, each curve being described by a specific value of G. Thus, the description of a curve is unique when, in addition to the position of its asymptotes, its shape is defined. These curves can be used graphically to determine the characteristic parameters of any arbitrarily located curve. Those parameters that can be approximated from any capillary pressure curve are: (1) total interconnected pore volume, represented by (V,),- ; (2) extrapolated displacement pressure, the abscissa axis; and (3) pore geometry factor, by the curve shape. PROPERTIES OF BRINE AND BACTERIA The fluid used throughout this study was an aqueous solution of sodium chloride, 25,000 ppm by weight. The addition of salt to fresh water was intended to eliminate clay swelling in the sandstone cores and to retard leaching in the limestone cores. The viscosity of the brine solution at 100F, the temp-

Jan 1, 1966

By A. Yoshikawa, T. Sugeno

To detenmine the origin of the sharp fatigue limit in many ferrous metals, S-N curvces were determined in push-pull fatigue at 18.6 kc per sec at room temperature and - 67°C for various kinds of iron. These differed in carbon and nitrogen content and grain size (25 and 2500 grains per sq mm). Sharp knees appeared in the S-N curves of fine-grained specimens with high interstitial content at both temperatures. In coarse-gruined specimens, knees were not sharp .regardless of the interstitial content and the testing temperature. It is concluded from these results that the controlling factov for the appeavance of the sharp fatigue limit is not strain aging but grain size. It has already been mentioned that iron and some iron alloys show fatigue curves characterized by a more or less sharp knee and a fatigue limit.' Explanations of the origin of this characteristic behavior of fatigue in iron appear to fall into two broad classes. One is concerned with interstitial impurities through the aging or locking effect and the other is based on intrinsic plastic properties of iron. Levy and sinclair2 and Lipsit and Horne3 have suggested that the fatigue limit observed generally in iron and steel is due to interstitial solute aging during the course of fatigue. On the other hand, Thompson and wadsworthl have suggested that the fatigue limit may be analogous to the sharp yield point, in the sense that alternating stresses above the fatigue limit are necessary to break down the dislocation locking. However, there are some experimental results which cast doubt on the general validity of these suppositions.4-8 Ferro and Montalenti' have recently suggested that the typical form of the fatigue curve of iron seems to come from an intrinsic plastic characteristic of the iron structure. This suggestion is based on the experimental result that pure iron containing an insufficient amount of interstitial solute to alldw appreciable strain-aging effects still showed a fatigue limit at room temperature. Although systematic investigations on the effect of grain size are scarce, it is thought that this effect is important in relation to the sharp fatigue limit. In fact, sensible differences were observed by Oates and Wilson4 in fatigue behaviors of a low-carbon steel depending on grain size. Ideally, the fatigue limit will be defined as the upper limit of stress amplitudes at which the life is infinite. This ideal fatigue limit can never be determined because of the limited time available for experiments, and moreover scattered data points often leave the construction of S-N curves somewhat arbitrary. And, it is usually stated that a material has a fatigue limit when its statistical S-N curve becomes apparently parallel to the N axis with a more or less sharp knee. The meaning of the fatigue limit and the sharp knee will be discussed in detail in relation to the construction of S-N curves. The origin of this sharp knee was investigated by determining S-N curves of iron for various interstitial contents, testing temperatures, and grain sizes. EXPERTMENTAL PROCEDURE The material as received was a commercial-grade low-carbon steel in the form of 4-mm-diameter rods with the following composition in weight percent: C, 0.09; Si, 0.017; Mn, 0.43; P, 0.006; S, 0.012; Cu, 0.11. Five different groups of specimens were tested. The grain sizes and interstitial contents in each group are given in Table I. The specimens of group C-02 were decarburized by annealing for 100 hr at 700°C, followed by an additional 100 hr at 850°C in wet hydrogen. Specimens of groups N-2 and N-9 nitrided in a mixed atmosphere of ammonia and hydrogen after an initial treatment identical to group C-02. The specimens of group C-9 were prepared in a manner similar to those of group C-02, except heat treatment was carried out in a vacuum. After polishing their surface with emery paper, the specimens of all groups were given a final anneal for 30 min at 930°C and air-cooled in an evacuated quartz tube. For group C*-9 this

Jan 1, 1965

By E. P. Whittemore

Experience with the ASC multipoint gas lift system was obtained in Colonia zone of the West Montalvo field near Oxnard, Calif. The wells in this pool produce from depths varying from 10,500 to 12,000 ft. Oil gravity is generally 14 to 15' API with a few extremes of 12 and 20" API. Some salt water is produced which causes some very viscous emulsions. Viscosities at 150F (which is the approximate wellhead temperature) vary from 5,000 to 100,000 SSU. Most of the production is by gas lift, although a few wells are produced by rod and hydraulic pump. About half of the gas-lift wells are on continuous flow and the remainder are on intermittent lift using large, ported, pilot-operated valves for single-point transfer of gas from casing to tubing. Gas-liquid ratios vary from about 6 to 10 Mcf/bbl of gross fluid lifted. Wells are produced to a 450-psi trap system. The following remarks will be confined to intermittent lift only, since this is the type of lift which has been achieved with the ASC valve system. The maximum gross fluid which has been produced by single-point intermittent lift is about 350 B/D in 3-in. tubing and 200 B/D in 21/2-in. tubing with gas-liquid ratios of approximately 7 to 9 Mcf/bbl. Some design changes could reduce this ratio. The ASC multipoint system has provided production as high as 480 BOPD in 21/2-in. tubing with gas-liquid ratios just under 4 Mcf/bbl. To be able to apply the multipoint system, it is recommended that a detailed explanation be obtained concerning transition-point pressure and stabilizer setting—what its significance is to the string design, how it may work for or against the operation of the well, how it is related to tubing sensitivity and how it affects the unloading operation. The unloading operation may only be of academic interest in a technical paper, but to the production foreman, unloading and setting the valves in operation is a very real problem and should be understood in detail. One item touched lightly in the paper was the unloading valve. This valve controls the maximum pressure at which the well can be operated. When lifting heavy viscous fluids, it is most important to set this valve for the maximum possible realistic operating pressure at the surface. If the well lifts easily, it is simple to set the ASC valves at a lower operating pressure and the unloading valve will remain closed; but if the well happens to be heavier to lift than anticipated, it may be desirable to operate on the unloading valve itself and use all the energy obtainable at the bottom of the hole. In the Colonia pool very heavy wet-gas gradients are experienced due to the viscosity of the liquid and the dense mist which is left behind a slug of fluid. There are many combination strings of 3- and 21/2-in. tubing. This aggravates the wet-gas gradient problem and provides wet-gas gradients of about 50 to 70 psi/1,000. An advantage which multipoint lift has provided is increased slug efficiency through better maintenance of pressure under the slug and decreased fall back as the slug passes up the tubing. By using multipoint injection, wet-gas gradients have been reduced to about 30 psi/1,000. This has reduced bottom-hole operating pressure and given a production increase. The ASC valve is not a simple device. It's operation is difficult to understand, and it must be understood to be used efficiently in gas-lift design. Operating problems are difficult to diagnose—whether they be caused by the fluid lifted, valve malfunction, lift gas rate, or operating pressure. Calculations and reasoning are required to find out what is causing the problem. Inherent in the ASC valve is the inability to create large pressure differentials across a slug. Large differentials may be required to overcome the inertia of very viscous fluid as it is being accelerated in the bottom of the hole. This is tied back to the design of the unloading valve and is one reason for the importance of setting the unloading valve for the highest possible operating pressure. ~u; to the narrow spread the ASC valves provide, it is impossible to cycle slower than about 24 cycles/day on choke control. If small production of 150 BOPD and less is expected, a surface time-cycle controller will be required if the most economical operation is to be achieved. To achieve a satisfactory operation, the operator must keep a record of any changes made in the operating pressure of the ASC valves. Anything which may cause changes in casing pressure in excess of the stabilizer setting will change the valve operating pressure, and if this is not noted from daily inspection of the well casing-tubing pressure recorder charts, the operator will lose control of the well. Significant results can be achieved using ASC valves; however, considerable knowledge is required to design with them, and attention to detail is required for satisfactory field operation.

Jan 1, 1965

By Luis Toron, Salvador Cortes-Obregon

John D. Price (Colorado Fuel and Iron Corp., Pueblo, Colo)—The paper on the coal fields of the Oaxaca district as prepared by engineers Toron and Cortes-Obregon of the staff of the Bank of Mexico bears witness to the thorough and careful way in which the men associated with this organization perform their work. There is little to be added to their paper in way of discussion other than to confirm and amplify some of their statements. Since the only extensive and well-developed field of coking coal lies in the northeastern section of the country adjacent to Sabinas in the state of Coahuila, it follows that blast furnace plants would be located in that same region. Two such plants are now operating at Monterrey and Monclova, using coke produced at the Sabinas district mines. But the nearer of these two plants is 600 miles from Mexico City and even farther from the center of population. Transportation of products from these mills to the market area is therefore expensive, both because of the distance and the difficulty of the terrain over which it must be carried. The development of an integrated steel industry closer to the center of population has therefore long been a goal toward which the Mexican technicians have been striving. While the presence of coal of some grade has been reported in many of the states, and many ideas have been advanced regarding its possible uses in iron and steel production, deposits of anthracite in Sonora and the various coals of the Oaxaca district as reported on in this paper are the only ones that have been explored in a serious manner. The coking coal from the Mix-tepec zone appears to offer promise of producing a coke which could be used in a standard blast furnace. Several problems are indicated, however: 1—The ash in the coal is high as mined, but indications are that it can be washed to an ash content of 15 pct with a recovery of 70 pct of washed coal. 2—Such washing would increase the volatile content from 17.4 pct to about 20 pct, and in a byproduct oven this should give a coke yield of close to 80 pct with an ash content of coke under 20 pct. 3—A free swelling index of 5 appears low for a good coking coal, and below that of the coals from the Sabinas district, which show between 6 and 9. But washing of the coal should result in an improvement in this regard; in the United States coals from Utah with an index even lower than 5 have made a usable coke. 4—A coal with volatile as low as 17.4 in raw coal and 20 in washed coal would come close to being classed as a low-volatile rather than medium-volatile coal, and low-volatile coals are notorious for their high expansion properties. Several plants in the United States are making coke from straight medium-volatile coal of 26 to 28 volatile content, and one at Rosita,, Mexico, from coal of 25 volatile. But no plants to my knowledge are using coal as low as 20 volatile. Since the Rosita coal appears to be a borderline coal from the angle of its expansion properties the coking of one of the straight lower volatile must be approached with caution. 5—There are few coals possessing any degree of coking properties which cannot be used in coke production by careful attention to its preparation and blending. The fact that coals of other types are available in this same region make improvement through blending very possible. 6—There are other workable methods of reducing iron ore other than the conventional coke-blast-furnace method. These will not be discussed here but it is known that their use has been considered. The technicians of not only Mexico but also of the other Latin American countries are keenly aware of their natural resources and their national needs. This paper emphasizes the fact that the Mexican technicians are working on their problem and attempting to speed the day of self-sufficiency for their country. Salvatore Cortes-Obregon (author's reply)—I wish to thank Mr. Price for his kind remarks. The Mixtepec coal as shown in Table II has 30 pct ash and a free swelling index of 5, but when the same coal is washed to 15 pct ash it has a free swelling index of 8 to 9 and the volatiles increased from 17.4 to 20.7 pct. A satisfactory coke has been produced from blends made in the Mexican laboratory using at least 40 pct of the Mixtepec coking coals with the other Oaxaca non-coking coals. Koppers in Germany report good coke obtained from the Oaxaca coal with a blend of 80 pct Mixtepec coal. Consideration is being given the possibility of using methods other than the conventional blast furnace for the reduction of iron ore near the Oaxaca area; electric furnaces appear promising. The non-coking coals could be used to produce cheap electric energy and the coking coals to make metallurgical coke.

Jan 1, 1955

By Z. F. Oszter

This paper discusses the materials handling problems encountered in cement plants in the areas of raw and finished product grinding and cement storage. It illustrates by examples taken from actual installations, how optimum feed rate is maintained in grinding circuits with properly applied bin dischargers (vibrating hoppers, Mexican feed boxes and aerated bin bottoms) for materials with high moisture content and/or difficult handling characteristics. Initial operating problems with the equipment will be discussed and maintenance requirements will be indicated for comparison purposes. In a modern cement plant, properly functioning bins and feeders are as essential as process control, for no process control system, however sophisticated, can perform efficiently if raw materials do not flow into the plant processing facilities at the required time and rate. Most materials used in the manufacture of portland cement - such as silica sand, iron ore fines, crushed limestone, gypsum, clinker or cement itself-have difficult handling characteristics under specific conditions. High moisture content or even dampness due to condensation may create bin hang-ups, arching and ratholing, all of which result in improper functioning of the feeders or proportioning equipment. The occurrence of these problems in existing cement plants spurred a search for the corrective devices discussed in this article. The equipment investigated included bin dischargers in the form of vibrating hoppers for moist raw materials, Mexican feed boxes for clinker and aerated bin bottoms for cement. Operating experience with different types of bin dischargers showed that the actual performance of some did not come up to the expected level under prevailing operating conditions. The aim of this paper is to record the results of this experience and there- by aid designers and/or operators in selecting the most suitable equipment for their requirements. BIN DISCHARGERS (GENERAL) A discharger is a device that causes stored solids to flow out of a hopper without metering the rate of flow as a feeder does. The investigation discussed here encompasses the following types of bin dischargers used in a cement plant: (1) vibrating hoppers (horizontal motion), (2) whirlpool hoppers (horizontal and vertical motion), (3) aerated bin bottoms (air fluidization), and (4) Mexican feed boxes. Most of the bins used in the cement industry are of the plug flow type. Raw materials with difficult handling characteristics such as high moisture content may cause flow problems when fed through small openings (12 x 18 in.), even if two walls of the bin are vertical.' The most appropriate solution to these flow problems would be the installation of separate mass flow bins for the different materials such as silica sand, iron ore fines, crushed limestone and gypsum, with each bin correctly designed from test data obtained by measuring the flow properties of representative samples. However, limited storage capacity and plant space usually preclude this approach. Quantitative design of bins for flow based on measured flow properties of representative samples of the bulk solid to be handled is no longer an "engineering art," thanks to the basic research done on the flow of bulk solids in bins during the last decade.2, Where space limitations exist, however, or where improvements are being sought for equipment already installed, bin dischargers offer a satisfactory solution to flow problems. They can be applied effectively under flat-bottom bins or mounted to existing plug-flow bins with very little modification to these facilities. In both the foregoing cases, plug flow will be maintained, but the flow channel will be considerably increased so that it is larger than the discharger's bin-side diameter, and material flow can be sustained if the discharger is correctly selected. The main purpose of a bin discharger is to draw the bulk solids across the whole area of the bin outlet, to make the opening fully effective. As long as the bin-side opening of the discharger is large enough

Jan 1, 1969

By Keith R. Bock, Norman Parlee, Albert M. Barloga

Coils of pure iron and iron-carbon alloy wire (0.05 to 0.80 pct C) and sufficient sulfur to saturate the solid phase were equilibrated in evacuated or argon filled tubes. After rapid cooling, and removal of the outside nonmetallic layer, the wires were analyzed for carbon and sulfur and the data used to construct an Fe-S binary and isotherms of the Fe-C-S ternary in the range 1149° to 1427°C. THE solid solubility of sulfur in steel is of interest in connection with such phenomena as hot shortness, burning, and so forth. "Burning," the more or less permanent damage that some steels suffer when heated for forging or rolling, has been shown to be related closely to the behavior of sulfur and less closely to carbon and oxygen.' Attempts to interpret burning phenomena in steels fail because of lack of data on the Fe-C-S and the more complex systems in this family. Rosenqvist and Dunicz 2 and Turkdogan, Ignatowicz, and pearson3 have largely elucidated the Fe-S diagram in the region of interest but no information on the Fe-C-S diagram in this region appears to be available in the literature. This paper deals with the elucidation of the Fe-C-S diagram in these interesting ranges. The method employed is different from those used by previous workers2'3 on the Fe-S system. EXPERIMENTAL METHOD Pure iron wires (Ferrovac E) or iron carbon alloy wires of 1 mm in diam were cleaned with acid and acetone, coiled, and placed in silica tubes (7 mm OD and 5 mm ID) previously closed at one end. Enough sulfur was added to assure saturation of the solid iron phase. The filled tubes were either simply evacuated and sealed, or filled with argon at a reduced pressure and sealed. The argon was required at the higher temperatures to prevent collapse of the tubes. The filled and sealed tubes were placed in the middle uniform temperature zone of a Globar tube furnace and equilibrated at temperatures ranging from 2100°F (1149°C) to 2720°F (1493°C). After equilibration the tubes were removed from the furnace and quenched in air or water, the form of quenching being found to have no effect on the results. The tubes were broken open and the coils were placed in a 1 : 1 HC1 solution to remove the sulfide-rich layer. The coils were then cut into small pieces and analyzed for sulfur and carbon. In the early stages of the investigation different equilibration times ranging up to 17 hr were tried and the cores of the wires were analyzed to test for saturation. One hour appeared to be sufficient to reach maximum sulfur content at 1454°C and 2 hr sufficient at 1149°C. The practice adopted was to use at least 3 hr at the higher temperatures and at least 5 hr at the lower temperatures. The iron-carbon alloy wires used were made by carburizing pure iron wire with carbon monoxide gas in a one inch diameter ceramic tube at about 1204°C. Differing carbon contents were obtained by allowing fairly large coils to react with the gas for varying lengths of time at a flow rate of 800 cc per min. The reacting times ranged from 15 min to 4 hr depending upon the amount of carbon desired. The furnaces used were controlled by means of Leeds and Northrup Speedomax Type H instruments with D.A.T. Control attachment. Thermocouples were calibrated against the melting points of gold and copper. The temperatures recorded appear to be accurate within i2.7OC. Starting with run No. 85 and continuing with the

Jan 1, 1962

By C. C. Martin, M. C. Fuerstenau, R. B. Bhappu

Experiments revealed that quartz could not be floated in conductivity water at any pH with a long-chained sulfonate as collector. Various cations, Fe+++, Al+++, Pb++, Mn++, Mg++, Ca++, are shown to function as activators when the pH is such that the cation hydrolyzes. As the solubility product of the various cation-sulfonates was exceeded at the concentrations involved, precipitated cation-hydroxy-sulfonates, e.g. Fe (RSO3)2OH must be functioning as the collector in these systems. A number of interesting and important phenomena were revealed during a study of the response of beryl to sulfonate flotation.' Perhaps the most important observation is the effect of various cations on flotation response. Experiments showed that in the vicinity of pH 3.0, about three ppm Fe+++ increased recovery of leached beryl from about 60 pct to 100 pct, whereas about four ppm Ca+ + reduced the recovery from 60 pct to 2 pct. With the exception of ferric iron, all of the other cations investigated reduced flotation recoveries drastically at this pH. As postulated in the paper Sulfonate Flotation of Beryl,' even with low concentrations of polyvalent cations and sulfonate (on the order of one ppm), cation-sulfonates are precipitated immediately in solution. That is, blueish-white clouds were noted to form after given additions of salts to known solutions of sulfonate. Furthermore, it was also shown that one-hundred-fold more ferric iron than necessary to form the precipitated cloud had to be added to effect complete flotation. In other words, some form of solid ferric sulfonate was found to function as the collector in these systems. The form of the compound was suggested to be Fe(RCO3)2OH and that bonding to the surface occurs through the hydrogen of the hy-droxyl. As the form of this compound was shown to be dependent on the extent of hydrolysis of ferric iron, it is to be expected that ferric iron will be unique in some pH region. Further, if the presence of hydroxyl in the neutral precipitated compound is of primary importance, then other cations should function similarly but at different values of pH, since they will hydrolyze at different concentrations of hydroxyl ion. To determine whether the premise of hydrolysis is correct, the study was extended to quartz systems in the presence of various cations at various values of pH. EXPERIMENTAL MATERIALS AND METHOD Sodium alkyl aryl sulfonate was chosen as the collector for the experimental work. The reagent, sup plied as a courtesy sample by the Shell Chemical Co., has the following physical properties: Physical form Solid — Finely ground Sulfonate content, pct wt 95 to 97 Molecular Weight 450 to 470 Number of carbon atoms in hydrocarbon chain 25 to 30 Melting Point, OF 250 to 260 Solubility Soluble in water with gel formation at concentrations higher than 25 pct by weight sulfonate This particular reagent was chosen because of its high active sulfonate content, long-chain length, and solid form. With the exception of the frother and sulfonate, all other chemicals were of reagent grade quality. Conductivity water was used in the experimental work. This water, made by passing distilled water through an ion exchange column, had an average measured conductivity of one micromho. Experiments were conducted at room temperature in a small glass flotation cell with the following procedure: 1) A redetermined amount of water and salt solution (e.g. CaC12 . 2H2O) were combined and the pH adjusted to a given value. 2) A given amount of sulfonate was added so that final solution volume was 130 cc. 3) One drop of Dow 250 frother was added. 4) Five g of quartz were added and the system conditioned for five minutes. 5) The pH of the system was measured (termed flotation pH).

Jan 1, 1963

By W. N. Miner, R. O. Elliott

Electrical-resistivity )measurments were made be-trueetz room temperatrive and 1.5 oK on five different stocks of cerium metal, and the results were correlated with the types, amounts, and distribution of the impurities present. Magnesium, which appeared to be distributed as a solute throughout the cerium matrix, was found to have a large effect on the residual resistivity, 1 at. pct Mg being sufficient to increase the residual resistivity by about 10 microhm-cm. The formation of ß by repeated thermal cycling between room temperature and liquid-helium temperature is attributed to the deformation and faulting that occurs during the ? D a transformation. In preparation for starting a proposed study of cerium-rich alloys, we examined samples from severa1 stocks of cerium metal in order to choose the purest stock for future use in making the alloys. Low-temperature electrical-resistivity data were used to gain initial information about the purities of the stocks. Additional information was subsequently obtained from chemical analyses, density measurements, optical metallography, and electron-micro-probe analyses. Cerium that has been annealed at high temperature and cooled to room temperature has the fee crystal structure (? phase) with a = 5.16Å. When cooled further, in the region between about 250" and 150°K, some of the fee phase transforms to the double hcp form (ß phase) with a = 3.68 and c = 11.92Å. The remainder of the room-temperature fee form (and perhaps some of the hexagonal phase) transforms electronically in the vicinity of 100° K to the more dense or "collapsed" fee form (a phase), a - 4.85A, by transferring about 0.5 electron per atom from the 4f level to the valence band. Very little work related specifically to the effects of impurities on the physical properties of cerium has been reported in the literature. Gaume-Mahn1, 2 has studied the effects of calcium, magnesium. iron, silicon, and tantalum on the electrical resistivity and magnetic susceptibility of cerium between 60o and 300°K. Gschneidner, Elliott, and McDonald3 have discussed the effects of total impurity content on the hysteresis of the a —? (reversible) transformation and on the formation of the 4 and "a-? intermediate" phases. Smith and Morrice* have reported the effect of total impurity content on the resistivity of cerium at 4 and 293°K. In none of these reports, however, is it apparent that the distributions of the impurities were determined, L.r., whether the impurities were present in inclusions or in solid solution in the matrix is not mentioned; nor does it appear that the phase purity of the cerium was always considered. In the present study, five different stocks of cerium were involved. One was represented by a single sample of electrowon cerium, which was all that had remained from a stock of high-purity material obtained from the U.S. Bureau of Mines, Reno, Nev., in 1959. The other four stocks were recently purchased from commercial sources and were specified to be 99.9 wt pct pure. EXPERIMENTAL Cerium stocks designated as #1, #2, #4, and #5 were obtained from commercial suppliers for the specific purpose of finding a reliable source of high-purity cerium. Although the purity level of the cerium was specified to be 99.9 pet, the specifications did not include limits as to the amounts of various impurities that would be acceptable. Consequently the chemical analyses of the different stocks varied widely. Stock #3, which had been obtained from the U.S. Bureau of Mines and which was the purest on the basis of chemical analyses, consisted of a single specimen of electrowon cerium. It was included in the examination for purposes of comparison. Table I contains the chemical and spectrographic-

Jan 1, 1968

By J. C. Griffiths, D. W. Bennion

A mathematical model, which does not assume a priori that stratification exists, but was designed to test for the stratification was developed. The model segmented the reservoir horizontally into areas of common variance, then divided it vertically into strata (if strata were present). Next, trend surface techniques were used to determine the lateral extent and variation of each stratum. The model was tested on two reservoirs (sandstone and limestone) using porosity as the reservoir rock property. The sandstone reservoir contained approximately 60,000 samples from 2,000 cored wells, while the limestone reservoir had approximately 24,000 samples from 430 cored wells. Within the areas of common variance tested, the model was able to distinguish four separate lithological units or zones in the sandstone reservoir, four in the Marly section and seven in the Vuggy section of the limestone reservoir. INTRODUCTION To accurately predict the movement and production of fluids from a reservoir rock, it is necessary that the reservoir engineer have a knowledge of the rock properties and their distribution throughout the reservoir. For this study, a reservoir rock is defined as a solid containing interconnecting holes or voids which occur relatively frequently and are dispersed within the rock in a regular or random manner. The fraction of bulk volume that these interconnecting voids occupy is called porosity. The rock must also possess sufficient permeability to allow fluids to move through it and it may be composed of one or more strata or lithological units. A stratum or lithological unit as defined in this study is a body or volume of reservoir rock whose properties are so distributed that its lateral and vertical extent can be traced throughout or through a portion of the reservoir. Operationally, the existence of strata are question- able unless the variance within strata is less than among strata. Within the past few years computers have been widely used as a tool for calculating movement and production of fluids from reservoirs. Most prediction models require some type of a mathematical description of certain reservoir rock properties. In the past, most of these models have been relatively simple. Some have assumed vertical variation but no lateral variation, while others have assumed lateral variation but no vertical variation. In general, most of these models have not used functional relationships to predict reservoir rock properties as functions of position. The object of this study was to develop a model which would predict lateral and vertical variations of reservoir rock properties. Since samples are usually available from only a small portion of the total reservoir rock, it seemed logical that if measurements from these samples were to be used to infer the properties of the actual reservoir, the data should be treated statistically as a sample from the total population (reservoir). Therefore, the model developed in this study to predict reservoir rock properties was a stochastic model. The model was tested using porosity (a macroscopic reservoir rock property) measurements from two reservoirs — the first, a sandstone reservoir from which 60,000 samples were obtained; and the second, a limestone reservoir with 24,000 samples. Several investigators have proposed methods to determine the vertical and areal variation of a reservoir rock property. Stiles: Dykstra and Parsons 2 and Suder and Calhoun3 have all developed waterflooding prediction techniques which assume vertical variation but no lateral variation. Stiles developed a method to segment a reservoir arbitrarily using frequency distribution. Law4 suggested that porosity has a normal frequency distribution and that permeability has a log-normal frequency distribution. Elkins and Skov 5 proposed that when a reservoir is segmented vertically, the

Jan 1, 1967

By Lawrence A. Shepard, Richard H. Krock

A series of ductile, two phase W-Ni-Fe composites, sintered in the presence of a liquid phase, were tested in tension. Identical room temperature stress-strain curves were obtained for specimens containing from 80 to 92 wt pct W (58 to 75 vol pct W particles). The composites exhibited a maximum elongation of 29pct at room temperature, and 10.7 pct at 77 °K. The tungsten particles in the composite elongated by the same amount at these temperatures. Single-phase alloy specimens matching the composition of the composite matrix showed about one half the flow stress of the composites. The test results demonstrated that the mechanical properties of W-Ni-Fe composites are determined by the tungsten particles alone and are independent of matrix volume fraction or mean free path over the composition range studied. A general prediction of the plastic behavior of two-phase alloys or composites from a knowledge of the properties, size, shape, and dispersion of the individual components is not, at present, possible. The excellent theoretical and experimental exposition of the dispersion-strengthening problem applies only to cases where the stronger phase is present as fine particles and in volume fractions of a few percent. Materials in which the stronger phase constitutes a significant fraction of the volume, although of broad engineering interest, have received comparatively little analytical attention. This important group of materials is the concern of the present study. The degree of deformation of each phase during the plastic working of a ductile two-phase alloy has been established by Boas and Honeycombe,1 Clare-brough,2 and Clarebrough and berger.3 These workers empirically correlated the recrystallization temperature of the individual phases with the degree of cold work sustained by each. It was shown that the deformations of both phases are equal and therefore the same as that sustained by the composite whole when the volume fraction of the stronger phase exceeds 35 pct. At yielding, however, deformation begins in the weaker phase. A direct correlation of composite strength, as a function of deformation, with the strength of the individual components cannot be made. In the presence of stronger phase particles, the flow stress of the weaker matrix phase is enhanced by an additional hydrostatic component. The development and potential magnitude of matrix strength promotion upon straining has not even been experimentally established. Quantitative work in this field has been confined to ceramic-metal composites. These materials fracture at strains considerably below 1 pct by cracking of the ceramic particles.4-6 An indication of the potential constrained matrix strength capacity is given by the observation by Nishimatsu and Gur-land4 that the maximum tensile strength of a WC-Co composite is almost double the bulk strength of the cobalt matrix. The relationship between composite strength and geometry is ordinarily expressed in terms of the mean-free matrix path (referred to subsequently as MFP) between the stronger particles.* This factor takes into account both the volume fraction of matrix phase and the average size of the particles. unke17 concluded that composite strength varied directly with the decreasing log MFP, following the Gen-samer8 relationship. However, Gurland and ~radzil~ reported a strength maximum at an MFP of the order of 1 p in WC-Co composites. A ductile two phase metallic composite was chosen for the present investigation of composite mechanical properties with the view of avoiding the difficulties and ambiguities noted above. The W-Ni-Fe heavy alloy of Green, Jones, and pitkin9 contains 80 to 94 wt pct W, and nickel and iron in a weight ratio

Jan 1, 1963

By C. Dalton, H. L. McGill

The purpose of this study was to develop a method for predicting the time delay created by metering a fluid through an annular restriction with translational motion of the boundaries. This time delay is utilized to control the opening of valves in various down-hole oilfield tools and forms the basis for the "hydraulic jar". The hydraulic jar is used to give upward impact forces to tubing or other tools that become lodged down-hole. The design of these jars has primarily been by trial and error; varying the clearances, fluid viscosities and dimensions of the jar until a satisfactory time delay was achieved. Obviously, a better understanding of the fluid flow characteristics is required to produce designs that are satisfactory and to eliminate the trial-and-error approach. Previous studies of flow in annular spaces were based on constant length of the annular restriction. Any motion of the boundaries of the annular space was confined to a rotational motion, like that encountered in journal bearing design. Consideration of the effect of translational motion of one boundary with respect to the other was necessary in applying this previous work to a time delay device. A prediction method was developed that takes into account this relative translational motion. Equations were developed to predict the relative velocity, the flow rate, and the time delay in a concentric annular space. The effect of eccentricity also was included. Experimental data were obtained and closely correlated with the prediction method. This correlation is presented graphically in terms of predicted and actual relative velocities for various annular clearances and pressure drops. This extension provides design equations which predict the performance of the hydraulic jar, thereby eliminating trial-and-error solutions and reducing the development time required for satisfactory designs. INTRODUCTION The design of down-hole oilfield tools which use hydraulic delay systems based on flow through annular orifices with relative motion of the boundaries is usually a trial-and-error process. The annular clearances, fluid viscosities and dimensions of the mechanisms were varied until satisfactory delay times were obtained. The development of these tools has been expensive, and a better understanding of the characteristics of annular flow with relative motion of the boundaries could result in superior and less expensive designs. The purpose of this investigation was to provide a method of predicting the performance of down-hole tools which use annular flow and, in particular, the performance of a down-hole tool known as the "hydraulic jar". The hydraulic jar is used to free pipe, packers and other tools that become lodged down-hole. The jar, located in the tubing string between the surface and the object which is stuck, releases the object though the repeated application of upward impact blows. The jar is basically a slip joint that is restrained in a collapsed condition until the desired tubing tension has been developed. The slip joint is then released and accelerates rapidly until it shoulders, delivering an impact blow or jar to the object. The delay obtained by metering the fluid through the annulus formed by the jar valve and housing provides the time delay necessary to stretch and apply tension to the tubing or pipe. Fig. 1A gives a schematic diagram of the hydraulic jar and brief details of its operation. The metering of the hydraulic jar begins when the tension force being applied to the jar by the pipe creates a pressure differential across the valve that is approximately equal to the force divided by the displacement area of the valve. This pressure differential induces flow through the annular restriction formed between the valve and housing. The valve moves with respect to the housing as the

Jan 1, 1969

By Robert Maes, Robert de Strycker

The Fe-Ni-As phase diagram has been established by the study of about a hundred alloys, by microscopic observation, and by thermal analysis, with arsenic contents up to 50 pct. The iron and nickel arsenides present extensive solid-solution fields, owing to the substitution of nickel by iron and vice versa; the extent of the solubility field of each compound has been determined with an accuracy of ±1 pct. In the investigated range of compositions, the solidification reactions were established, and the temperatures of the invariant reaclions detevmined with a preciston of iZ°C in the most favorable and ±5°C in the least favorable cases. The isothermal lines of the liquidus surface have also been drazum, with an accuracy estimated at i5°C. Reactions in the solid state, which take place for the formation or the decomposition of certain phases , were investigated in detail. ThE constitution diagram of the Fe-Ni-As system is fundamental for the understanding of the properties of nickel speiss; these by-products of the extraction of certain nonferrous metals indeed often contain the three elements iron, nickel, and arsenic as main constituents. The Fe-Ni-As system has already been the object of earlier investigations. In 1932, Guertler and Savels-berg' have presented some elements of the ternary phase diagram, for arsenic contents up to 55 pct (all percentages in this paper are given in weight percent), including the vertical section FezAs-Ni&s2. This investigation is, however, incomplete and certain anomalies suggested the necessity to verify the results published by these authors: the section Fe2As-Ni5AsZ, for example, is presented as quasi-binary, with a field of complete miscibility in the solid state, even though these compounds do not have the same structure. Recently, ~useck' established an isothermal section at 800°C in the Fe-Ni-As system by X-ray diffraction and microscopic examination of water-quenched alloys. These techniques are sometimes inaccurate for the determination of the fields where alloys are liquid at the investigated temperature, and they may lead to erroneous conclusions when a compound exists at the investigated temperature but is not stable at room temperature and cannot be maintained by quenching. LIMITING BINARY PHASE DIAGRAMS The Fe-As constitution diagram has been the object of several investigations which have been reviewed by Hansen and Anderko.3 For the Ni-As phase diagram, a recent study has been effected by Yund.4 In this system, Heyding and calvert5 have determined the existence of a compound of unidentified structure at arsenic contents slightly higher than those corresponding to Ni5As2 and at temperatures lower than about 200°C; by analogy with iron and cobalt arsenides, this compound could correspond to the formula Ni2As, as suggested by Kulle-rud; although this is not definitely established. In the region of arsenic contents from 35 to 55 pct, Fried-rich7 detected anomalies in the solidification reactions; an interpretation of these anomalies was given by Hansen ,8 which assumed that the solidification reactions in practice are not in equilibrium, but are meta-stable. The main features of the Fe-Ni constitution diagram are the existence of a complete miscibility field, at least at high temperatures, for the fcc phase (which will be designated by My in the remainder of this work) and of a limited solubility field for the bcc phase (designated by Ma). EXPERIMENTAL METHODS The fundamental technique used in this investigation was microscopic observation, which allowed the determination of the reactions occurring during solidification of the alloys, and possible reactions in the solid state.

Jan 1, 1968

By T. E. Seidel, A. U. MacRae

The dependence of the electrical and crystalline properties of silicon containing ion implanted boron atoms have been studied as a function of the incident dose, substrate temperature, and annealing tempera-ture. Hall effect studies show that the compensation due to defects is negligible after annealing at and above 500°C for a dose of 1013 ions per sq cm. Steps occur in the isochronal annealing curve for the sheet resistance at 300°C and 800°C and the mechanisms responsible for these steps are discussed. The range of the ions at 150 kev is 3500 * 5001. Electron diffraction studies on samples implanted with 150 ker ions show the existence of damaged (polycrystalline) silicon closer to the surface than the range of the ions. Etching and subsequent annealing of etched and unetched specimens suggest that the poly crystalline silicon is closer to the surface than the bulk of the distribution of ions. THE purpose of this study is to investigate some electrical and crystalline properties of ion implanted boron in silicon. Sheet resistance, Hall measurement, depth distribution by chemical staining, thermal probe and electron diffraction have been used to study the properties of samples implanted in a random direction with ions having energy up to 150 kev. The boron-silicon system has been the subject of considerable study'-7 due to its possible application in the fabrication of electronic devices. The implanted ions produce considerable crystalline damage to the substrate crystal before they lose most of their primary energy and come to rest. Fortunately, most, if not all, of this damage can be removed by subsequent annealing or by holding the substrate at an elevated temperature during the implantation process. Sheet resistivity measurements of the p-type layers produced by boron implantation into room-temperature silicon, reveal that the annealing occurs in a two-step process: the first occurring at -300°C and the second at -800°C.' Electron-microscope and low-angle electron diffraction studies of this damaged material indicate that the bombarding ions produce an amorphous silicon region about the trajectory of each ion.' At high doses, these damaged regions overlap and an amorphous surface layer results. Furthermore, the 800°C annealing step, appears to be associated with the elimination of this gross damage. The depth distribution of the as-implanted and also the annealed samples have been studied by an angle lap and staining technique as well as by Hall measurements.3-7 We find the experimentally determined ranges of the boron are less than those calculated by the Lindhard-Scharff-Schiott theorys (hereafter referred to as the LSS theory). Our results are in essential agreement with pre- vious results. In addition, we have made Hall measurements on samples bombarded with a low dose (2 x 1013 ions per sq cm) to minimize impurity banding effects, and find that deep ionized levels influence the carrier concentration and room-temperature hole mobility until the sample is annealed at temperatures above 500°C. Above this temperature, the data corresponds to that obtained when silicon is doped with boron by conventional techniques. We have also made a depth distribution determination of the crystalline damage due to the implanted boron by using glancing angle electron diffraction techniques. It was found that the gross damage depth distribution of heavily implanted samples is closer to the surface than the depth distribution of the electrically active boron. By removing this gross damage, without affecting the distribution of the implanted boron, we have determined that the 800°C annealing stage is not exclusively due to the improvement in the crystalline perfection of the amorphous material but also appears to be associated with the annealing of some point defects. These results and their relation to previous results are discussed in the following sections. IMPLANTATIONS Boron ions were implanted into silicon samples using an accelerator capable of operating from 10 to 150

Jan 1, 1970

By Attilio Ligasacchi

The production of magnetite in copper smelting has been the subject of wide investigation in the past.1-9 It is accepted that most of the magnetite is produced in the converters and returned to the re-verberatory furnace with the converter slag. Part of this magnetite is circulated back to the converters in the matte. As a result, if proper fluxing practice is not followed, and temperatures in the converting process are not strictly controlled, an undesirable circulating load of magnetite may occur. This might ultimately lead to saturation in the reverberatory bath with consequent settling of the magnetite at the furnace bottom. This paper contains the cumulative results of a series of investigations on the formation and control of magnetite presently being conducted by the metallurgical department of the Chino Mines Division of Kennecott Copper Corp. THE METHOD USED FOR THE DETERMINATION OF MAGNETITE AND OTHER MAGNETIC MATERIALS Conventional wet chemical magnetite analyses are neither reproducible nor theoretically sound. These analyses are sensitive to time, temperature, and heating rate-all factors which tend to decrease their accuracy and reproducibility. To overcome the unreliability of the above method, H. W. Mossman1 developed an inductance bridge magnetite "meter" which completely replaces the wet chemical analysis with the advantage of quicker and reproducible results. Fig. 1 is a schematic diagram of the meter in use in the Chino metallurgical laboratory. When a pulverized specimen with magnetic properties is placed in coil A of the inductance bridge circuit, the permeability of the core is increased, causing a voltage drop across the coil. The unbalance is amplified to operate the balancing motor which adjusts the movable iron core in coils C and D to bring the bridge back to balance. An indicator arm connected to the movable iron core indicates directly the percent, or weight, of the magnetic material in the specimen. Coil B acts as a reference coil and a zero adjustment. The meter is standardized with a set of various amounts of -100 mesh pure magnetite crystals (99.8 pet Fe3O4) thoroughly mixed in a test tube with -100 mesh quartz sand. These mixtures are combined with inert transoptic powder and heated to 140° to 160°C (melting point of the powder) and subsequently cooled to room temperature, solidifying the whole mixture, in this manner eliminating any possible preferential orientation of the magnetite grains. In order to determine the effect of the particle size" of the magnetite and other magnetic materials on the meter determinations, five size fractions were obtained by wet screening of the 99.8 pet pure magnetite and several reverberatory matte and slag samples. The following size fractions were investigated: -80 + 100; -100 + 150; -150 + 180; -200 + 250; -250 + 300 and -100 mesh. The deviation of the meter readings for the various sizes of each sample never exceeded ±0.3 pet. This deviation is considered well within the accuracy needed for common smelter operation procedures. Further investigation proved also that grinding slag and matte samples in our laboratory disc grinder normally produced a powder constituted of 95 pet -200 + 300 mesh. For practical purposes a - 100 mesh size was used as standard for our studies. It was also found that for rapid determinations tamping of the sample in the test tube was sufficient and gave comparable results with those samples solidified with transoptic powder. For the purpose of this study we define the magnetite content in copper smelter slags and matte as the equilibrium amount of crystalline Fe3O4 existent in the specimen at a given reference temperature, excluding any secondary magnetite precipitated in the specimen during nonequilibrium conditions. MAGNETITE DETERMINATION IN SLAGS Since it is impossible to ascertain the status of equilibrium of a conventional slag sample, due to

Jan 1, 1965

By L. O. Storey

The McDermitt mine was found as a separate mercury ore-bearing occurrence approximately 305 m (1000 ft) northeasterly and in a different geologic setting from the old Cordero mine, which had been the district's major producer. The following briefly describes the discovery of the new deposit and the Placer Amex lnc. involvement. The Cordero mining property was under lease by Sierra Mineral Management from the owner Fred H. Lenway in February 1971, when the property was brought to the attention of Placer Amex Inc. (known at that time as American Exploration and Mining Co.) by Richard F. Hewlett of Sierra Mineral Management. The property was first examined by Amex geologists in the spring of 1971, who reported favorably on the project. During that year, Sierra Minerals Management had retained Dr. Larry Walters as a consultant. He proceeded to drill test mercury soils anomalies which he had found northeast of the old Cordero mine in an area overlain by recent alluvium and gravel. An area known as the D area was subsequently blocked out and ascribed a tonnage of about 18 1 kt (200,000 st) which occurs in the vicinity of the eastern edge of the currently defined McDermitt ore body. This mineralization was found in flat lying mercury-bearing lakebeds, unlike the steep structurally controlled ore which was mined in the old Cordero mine. It became apparent to Placer Amex that since the chances of much more extensive ore were good in such a geologic environment, additional exploration was warranted. A joint venture was formed in October 1972 by American Exploration and Mining Co. and Sierra Mineral Management to further explore the Cordero property. The project exploration was fully funded and managed by Amex. Sampling began by Placer Amex in the D area where an open trench had exposed the lakebed ore. Visible cinnabar could not account for the high-grade mercury assays taken from the trench, and an unknown mercury mineral was suspected. Subsequent work established the presence of a chlorine-bearing mercury sulfide which may have enriched the lakebeds as a secondary supergene mineral. Early identification of the mercury-bearing mineral corderoite (Hg3Cl2S,) aided in the geologic concepts for the ore occurrence and was helpful in exploring for the expansion of the D area. Corderoite was found occurring with cinnabar in the lakebeds, and in places it is a major contributor. Specimens of the unique ore were given to Dr. Eugene E. Foord and Dr. Pieter Berendsen at Stanford University, who subsequently identified the chemical and crystallographic makeup of the new mineral and described it for the first time as a naturally occurring substance. Metallurgical work on mercury recovery from the newly described mineral was carried out by the Placer Development Ltd. laboratories in Vancouver, and an acceptable recovery system was worked out together with a planned mill design. Step out drilling from the D area by Placer Amex Inc. encountered extensive flat lying lakebed ore. The McDermitt ore body, as subsequently defined by drilling, was found to contain 2.7 Mt (3 million st) of open pittable ore, containing approximately 5 kg/t (10 lb per st) Hg. The ore body gently dips to the north and is found lying mainly in claybeds just above an opalite breccia blanket which overlies the earlier volcanics. Some ore is also found in the opalite breccia. The mine and mill were constructed under an operating joint venture with Placer Amex Inc. and Sterling Mineral Ventures, and production began in June 1975 at the rate of 20,000 flasks per year. The mine is currently in progress.

Jan 1, 1985

By C. van der Poel, R. L. Chuoke P. van Meurs

When an initially planar interface between two im-ttitcihle liquids is displaced at constant rate, U, nor-mat to the front, instability will occur for all rates greater than a critical rate. U, given by provided the Fourier decomposition of the spatial pertmrbcrtiott or deformation of the moving displacetment front contains modes with wavelengths. A, greater than a critical wavelength, given by tit there expressions, p, p and k, with the subscripts I and 2 distinguishing the two 1iquids under consideration, are density, viscosity and effective pemeability coefficients respectively; g is the absolute value of the acceleration due to gravity; and cos (zz') is the direction cosit7e between the vertical carlesian coordinate z' (positive upward) and the z coordinate norma1 to the initially plane macroscopic interface taken positive in the direction from Liquid I to Liquid 2. U is the average volumetric velocity (injected volume of liquid pet. unit time per unit of area normal to Z) and is positive for Liquid I displacing Liquid 2, negative for the reversc displacement and r2 is an effective inter-fucial tension (for displacements between parrallel plates, t0= s the ordinary bulk interfacial tension). further, there is a wavelength of maximum instability given by For natural perturbations, this wavelength will preclotninate and characterize the form of the instability as a quasi-sinusoidal deformation, i.e., viscous fitzgers of peak-to-peak separatiotr, A,,. Comparison of experitient with theory rei'eals tltat these considerations are recsonably valid for displacements in hoth parallel plate chantzels and nnconsoli-clutt'd glass ptnvder packs. These results are. tij particular interest for interpretation of data from laboratory experiments with regard to fielil applica/iotz.s. It is quite possible for A,, although st7zall relative to the large lateral dimension of the reservoir, to rxceed the greatr.st lateral extent of the lahoratory sample or model. If it does, frontal displace- ment occurs in the laboratory experiments, and oil INTRODUCTION Thc puposes of this paper are to present theoretical and experimental evidence for occurrence of macroscopic instabilities in displacement of one viscous fluid by another immiscible with it through a uniform porous medium and to compare available experimental data with some predictions of a theory of instability developed by the first author. The instabilities are referred to as macroscopic in the sense that spatially quasi-sinusoidal, growing fingers of the displacing liquid are formed, the width and peak-to-peak separation (wavelength) of which is large relative to a characteristic length of the particular permeable medium such as grain size. Visual models of two kinds have been used to obtain observations: displacement of oil by water-glycerine solutions through the flow channel formed by closely spaced parallel plates and displacement of oil by water with and without initial interstitial water through unconsolidated glass powder packs, employing the technique of matching indices of refraction. In all cases we have observed macroscopic instabilities or fingers under conditions predicted by the theory to be favorable for their occurrence. The phenomenon discussed here is not the production of streamers due to gross inhomogeneities such as permeability stratification of the porous medium. It is our object to show, on the contrary, that a quantitative prediction of finger spacing is possible in a porous medium known to be macroscopically homogeneous and isotropic throughout. The importance of the phenomenon in its influence on the configuration of oil and water with respect to oil production behavior was noted earlier by Engel-berts and Klinkenberg' who coined the term "viscous fingering". THEORY Necessary and Sufficient Conditions ior Instability anc Initial Kinematics There are several levels of increasing complexity in the theoretical description of instabiIity of fluid displacements in permeable media. Of these, the simplest description, ;tdapted to low permeability systems, is selected for presentation.' More inclusivc tlescriptions are reserved for separate publication.

By D. W. Levinson, L. F. Mondolfo, N. A. Gjostein

MANY investigations are reported in the literature on the age hardening of Al-Mg-Zn alloys but most of them are concerned mainly with mechanical property changes. The present investigation was started to determine the structural changes in an AlMgZn alloy. The alloy composition was chosen in the existing phase diagram" so as to lie on the quasi-binary section Al-(AlZn)Mg Actually, it was found that the diagram as published was inexact at the lower temperatures and that the phase in equilibrium with aluminum at the lower temperatures is not the ternary compound. Technique The alloy used for the investigation was prepared from high purity materials and the resulting composition as determined by chemical and spectro-graphic analysis was: Zn, 6.17 pct; Mg, 2.13 pct; CU, 0.03 pct; Si, 0.009 pct; Fe, 0.005 pct; and the remainder, aluminum. An ingot -in. thick was cast, homogenized, hot and cold-rolled to strip 0.064-in. thick for dilato-metric and powder work, and to foil 0.005-in. thick for single crystal work. Solution treatment temperature for all specimens was 475°±5°C for periods of at least 1 hr, either in air furnaces or lead baths. Aging was done in oil or Pb-Bi baths controlled to within ±3°C. The aging was investigated at the following temperatures: room (30°), 130°, ZOO0, 260°, and 300°C. For the dilatometric study, the length changes were measured continuously, at temperature, using strips of alloy 10 in. long, shaped as channels to eliminate undesirable deflection due to bending. To compensate for thermal expansion and contraction due to the variations in the temperature of the aging baths, the dilatometers were built of the same alloy as the one being tested, with the same shape and size, but with material fully annealed. Gages on which 0.0005 in. could be read and 0.0001 in. estimated were used to measure the dilation. Length changes were also measured at room temperature. Strips approximately 10 in. long were solution-treated, quenched, and measured in a supermicrome-ter in which 0.0001 in. could be read and 0.00001 in. estimated. The specimens were immediately aged for a given time, remeasured, and then the cycle was repeated, including the solution treatment. The error of measurement with this technique was somewhat larger, but the results confirmed those obtained by continuous measurements. The sheet used for the dilatometric work was checked to determine grain orientation. Both X-rays and etching methods revealed a rather coarse grain size but only a negligible amount of preferred orientation. Thus, volume changes calculated from length changes are not invalidated by preferential orientation effects. Pieces of the sheet used for the dilatometric work were also used for the metal -lographic work. Filings for the powder X-ray study were produced from a homogenized piece of sheet. These filings were magnetically purged of iron, screened to 200 mesh, sealed in a Pyrex tube containing an argon atmosphere, and solution-treated. When quenching, the glass tube was broken, so that the powder dropped directly into cold water. After drying with alcohol, the powder was sealed into capillary tubes containing an argon atmosphere and then aged for the required times. In spite of all these precautions the three strongest lines of a( Al2o) appeared in all patterns, and were used as internal standards. By the use of these standards, accurate corrections for film shrinkage were obtained. For the back reflection patterns, annealed nickel powder was mixed with the specimen to act as internal standard. A preliminary investigation was made in small (57.3 mm diam) powder cameras. From these patterns the changes which took place in the solid-

Jan 1, 1957

By Kenneth E. Gray

The turbulent flow drag coefficients, or friction factors, have been experimentally determined for the cut-tings normally encountered in drilling operations. The gas law and average drag coefficients for characteristically flat particles (limestones and Shales) and for angular to sub-rounded particles (sandstones) were used to extend Newton's equation to a more useful form. The resulting equations are expressions for the slip velocity of a particle as a function of the particle size and shape, the bottom-hole injection Pressure, and the Bottom-hole temperature. The velocities necessary to lift actual rock cuttings as observed in the laboratory were compared with velocities predicted from the derived equations. Results indicate that there is no single correct circulating velocity and that the commonly used 3,000 ft/min linear velocity may he sufficient to lift only very mall particles. The advantage, in terms of horsepower requirements, of cycling high pressure air to obtain lower compression ratios was shown. INTRODUCTION A promising departure from standard rotary hydraulic drilling is air or gas drilling. In many instances the rate of penetration and bit life increases resulting from this method have been very substantial. Further the application of air drilling in areas of lost circulation or water sensitive pay zones has been highly successful. Purpose For the drilling contractor, the practical question remains, "How much air pressure and volume should I have?' lt is necessary to have a reasonable knowledge of the air velocities and pressures required in air drilling operations for the proper design and most advantageous use of surface equipment and for adequate removal of cuttings from the borehole. Eased on experience, most operators agree that at least 3,000 ft/min linear velocity is necessary for satisfactory hole cleaning.' It was the purpose of this investigation to evaluate the drag coefficients, or friction factors, for the cuttings normally encountered in drilling operations, i.e., sandstones, limestones and shales, and to utilize the values obtained to develop an expression for the minimum velocities and pressures (hence the volumes) necessary to lift these cuttings. THEORY Several investigators have done extensive work on the ability of drilling mud to lift bit cuttings. They conclude that the ability of a drilling mud to lift cuttings is dependent on the density difference between the rock being drilled and the drilling mud, the cutting size and shape, the mud flow constants, and the flow state (laminar, transitional or turbulent) of the mud. However, drilling muds are not true fluids but have a variation of viscosity with velocity in laminar flow and a constant apparent viscosity in turbulent flow. A survey of the literature revealed that no experimental work had been done to determine the velocities necessary to lift actual rock cuttings using a compressible fluid such as air or gas. In 1953, Nicolson6 stated that in air drilling fluid mechanics the particles could be assumed spherical in shape and that ; constant drag coefficient of 0.5 could be used due to the highly turbulent flow of the circulating air. By equating the turbulent resistance on the particle to the gravitational force on the particle, he obtained the expression, ^ = 2.67^y ,.......(1) where V, is slip velocity of spherical particle, ft/sec; D is particle diameter, in.; p, is particle density, 1b/ft8; and p, is fluid density, lb/ft3. However, as stated earlier, most operators seem to agree that adequate hole cleaning can be obtained by circulating a sufficient volume of air to give a linear velocity of 3,000 ft/rnin in the annular space. Basic Equation The velocity with which a solid particle freely falls through a fluid will increase until the accelerating force, gravity, is equal to the resisting forces, buoyancy due to the fluid displaced by the particle and friction due to the relative motion of the particle through the fluid. When the accelerating and resisting forces become equal,

By N. Mungan

A study was made of the effects of wettability and interfacial tension on the immiscible displacement of a Iiquid by another Iiquid from porous media. The influence of viscosity ratio was also investigated. Porous media used were polytetrafluoroethylene (TFE) cores prepared by compressing TFE powder under different pressures. lt is shown that displacement of a wetting by a nonwetting Iiquid is always less efficient than the displacement of a nonwetting by a wetting fluid, all other things being equal. In the former case, the recovery efficiency can be increased substantially by either reducing the interfacial tension or increasing the viscosity of the displacing fluid. A qualitative discussion is given on the implications of this work to the recovery of crude oil by water flooding. INTRODUCTION The high cost of oil exploration and new recovery schemes makes it imperative that waterflooding be conducted under conditions favoring most efficient oil recovery. To improve oil recovery by waterflooding, it is essential that the role played by interfacial forces in the entrapment of residual oil be studied and understood. Interfacial phenomena in natural rock, connate water and crude oil systems are very complicated because of the complexity of the natural liquids found in petroleum reservoirs, because of our inability to adequately describe the geometrical structure of the porous media and because of a lack of understanding of physical and chemical interactions between the liquids and surface of the pores. The problem becomes further complicated when one tries to elucidate the role of interfacial phenomena in fluid flow. Numerous studies of the displacement of oil by water under different interfacial tension or wet-tability conditions have been made.lm3 These studies have been performed in silica, alundum or sandstone systems using water and paraffin oil and also some surface active material to control the interfacial tension s and the contact angle 0. Unfortunately, the high energies of various interfaces involved favor adsorption and orientation of the surface active material at the interfaces. Also the surface active material concentration at the interfaces exceeds that in the bulk of the liquid phases. Such surface excess may cause the surfactant distribution, the contact angle and the interfacial tension to differ from their measured static equilibrium values and makes interpretation of the displacement experiments difficult. Furthermore, as changes in s also lead to changes in cos 0, the role played individually by one of these parameters in the displacement becomes obscured by the effect of the other. To circumvent these difficulties, a low surface energy solid and true solutions or pure liquids should be used. Use of a low surface energy solid minimizes adsorption and orientation effects at the solid-liquid interfaces. By controlling s and cos 0 through use of selected pairs of pure liquids or true solutions rather than by surfactants, the adsorption effects at liquid-liquid interfaces are eliminated. In the present study TFE cores were used as the porous media. Liquids used were water, sucrose solutions, paraffin oils and benzyl, n-butyl and isobutyl alcohols. The interfacial tension was varied from 40 to 1.1 dynes/cm by suitably choosing the liquid pair. A surface active material was added to the water-oil system only in the case where interfacial tension of 0.5 dyne/ cm was desired. No precise changes of cos 0 were attempted. However, either the displaced or the displacing liquid could be made the one which preferentially wets the TFE surface. Using sucrose solutions and blends of paraffin oils proved to be a convenient way of changing the viscosity ratio between the displaced and displacing liquids. The present investigation examines the effect of interfacial tension, wettability and viscosity ratio on the immiscible liquid-liquid displacement from porous media. Observed data permit a

Jan 1, 1967

By E. C. Barfield, D. G. McCarty

Two-dimensional analyses offer considerable promise in providing the basic information required to effect more precise control of petroleum reservoir performance. This paper describes a method for conducting some engineering analyses of this type using a high-qreed digital computer. The general approach is to (I) develop a representative numerical model of the reservoir and (2) employ a suitable numerical technique to solve the basic equation of flow and to develop the required engineering information for the particular case represented by the model. This technique has been used to calculate pattern performance in connection with several field projects involving water injection into oil reservoirs. This type of analysis involves sequentially (I) calculating the pres-' sure distribution in the reservoir and (2) tracking the progress of the interface between the displaced and displacing fluids in a step-wise manner to provide a depletion history of the operation. The particular analysis presented in this discussion is subject to the restriction that the mobility of the displaced and displacing fluids be equal and assumes that the fluids are incompressible and that gravity and capillary pressure do not affect the shape of the flood pattern. Complete two-dimensional .flexibility is maintained with regard to definition of reservoir rock and fluid characteristics, placement of physical restrains and boundary conditions, investigation of flow characteristics in the reservoir, and movement of the displaced and displacing fluid interface. The results of these studies indicate that the highspeed digital computer is well suited for conducting reservoir performance studies in two-space dimensions and that dependable numerical techniques are available for making such analyses. INTRODUCTION One of the problems facing petroleum engineers is that of providing more precise engineering control over reservoir operations at all stages during the depletion cycle. The most direct approach to this problem is to simulate the reservoir with some type of physical or mathematical model and use the information developed from the behavior of the prototype to predict the performance of the actual reservoir. This principle is, of course, not new, and a variety of simulation techniques have been developed to provide information on the behavior of oil and gas reservoirs. One straightforward refinement of present reservoir analysis techniques is the progress from one-dimensional reservoir models to those where the reservoir can be represented in two- and three-space dimensions. Multidimensional reservoir analyses are desired in order to show not only what is happening in the reservoir but in addition to show where in the reservoir the various phenomena of interest are occurring. Thus far, two-dimensional analyses have been confined predominantly to potentiometric and electrolytic model studies. However, when large numbers of wells are involved and where large-scale models are required, such problems become very tedious and instrumentation problems become acute." It is possible through use of numerical methods and general purpose digital computers to avoid many of the instrumentation difficulties and to relax some of the restrictions of previous techniques. The purpose of this paper is to discuss some recent experiences with the use of high-speed digital computers for predicting pattern performance in fluid injection operations and to describe the numerical procedures which have been used for making such analyses. THE NUMERICAL MODEL The problem of representing the reservoir with an appropriate numerical model is composed of at least two important considerations. The first requires that an adequate description of the reservoir mechanism be formulated in the form of the basic differential equation describing the process. The second consideration is definition of the physical characteristics of the reservoir. The latter requires that representative basic data be obtained, examined, and interpreted to provide sufficient detailed definition of the pertinent physical properties throughout the system. The Differential Equation One of the basic assumptions made for the purpose of performing the two-dimensional studies being dis-