Search Documents

By M. Muskat, M. O. Taylor

Theoretical calculations have been made on the performance to be expected of gas-drive reservoirs for various characteristics of the oil and gas and the producing rock. Such performante has been expressed graphically as curves for the reservoir pressure and gas-oil ratios of the production as functions of the cumulative oil recovery. The latter has been expressed in terms of percentage of pore space of the reservoir rock. Such curves automatically give values for the ultimate physical recovery, if the latter is interpreted as the recovery obtained when the reservoir pressure has declined to atmospheric, or to any other pressure chosen as defining the state of practical complete depletion. Calculations of these pressure and gas-oil ratio histories have been made for conditions in which the oil viscosity, the gas solubility and shrinkage, the size of an overlying gas cap, if present, the permeability-saturation characteristics of the rock, and the amount of connate water, have been individually varied, The results serve to show the extent to which the ultimate recoveries are sensitive to the important physical parameters characterizing the oil reservoir. As is to be expected, the ultimate recoveries are found to decrease with increasing oil viscosity. Because of the predominant effect of the oil shrinkage associated with the liberation of the gas in solution, the ultimate recovery will decrease with increasing gas solubility. Increasing gas-cap volumes lead to higher recoveries, although the contribution made by the gas-cap gas is small as compared with the oil expulsion by the equivalent amount of solution gas. As a whole, the oil recovery when expressed in percentage of the pore space is not very sensitive to the details of the permeability-saturation relationship. However, if the rocks possess an equilibrium free-gas saturation, the rise in gas-oil ratio will be retarded, and the 'ultimate recovery somewhat increased' The oil shrinkage associated with the liberation of the dissolved gas also leads to the result that as long as the connate water is immobile and the permeability ratio curve for the rock is a function only of the total liquid saturation, the stock-tank oil recovery will be less for a sand containing no connate water than for One with an Original water content as high as 3' per cent. The space voidage in the former case will be somewhat greater, but the effect of oil shrinkage will lead to smaller values for the equivalent stock-tank recovery. These calculations also give data showing how the productivity indexes for the producing wells will vary during the reservoir history. Because of decreasing Permeability to the oil and increasing oil viscosity, the productivity index will fall continuously as Production proceeds, and may finally reach values as low as 10 per cent of the initial productivity index. Introduction In a previous paper1 was developed the basic theory for the prediction of the production histories of gas-drive reservoirs. The data required for the application of this theory included the characteristics of the petroleum fluids, such as the viscosity of the Oil and gas phases, the solubility of the gas in the oil, and the oil shrinkage— all expressed as functions of the reservoir pressure—and the permeability-saturation characteristics of the producing rock, as

Jan 1, 1946

By M. Muskat, M. O. Taylor

Theoretical calculations have been made on the performance to be expected of gas-drive reservoirs for various characteristics of the oil and gas and the producing rock. Such performante has been expressed graphically as curves for the reservoir pressure and gas-oil ratios of the production as functions of the cumulative oil recovery. The latter has been expressed in terms of percentage of pore space of the reservoir rock. Such curves automatically give values for the ultimate physical recovery, if the latter is interpreted as the recovery obtained when the reservoir pressure has declined to atmospheric, or to any other pressure chosen as defining the state of practical complete depletion. Calculations of these pressure and gas-oil ratio histories have been made for conditions in which the oil viscosity, the gas solubility and shrinkage, the size of an overlying gas cap, if present, the permeability-saturation characteristics of the rock, and the amount of connate water, have been individually varied, The results serve to show the extent to which the ultimate recoveries are sensitive to the important physical parameters characterizing the oil reservoir. As is to be expected, the ultimate recoveries are found to decrease with increasing oil viscosity. Because of the predominant effect of the oil shrinkage associated with the liberation of the gas in solution, the ultimate recovery will decrease with increasing gas solubility. Increasing gas-cap volumes lead to higher recoveries, although the contribution made by the gas-cap gas is small as compared with the oil expulsion by the equivalent amount of solution gas. As a whole, the oil recovery when expressed in percentage of the pore space is not very sensitive to the details of the permeability-saturation relationship. However, if the rocks possess an equilibrium free-gas saturation, the rise in gas-oil ratio will be retarded, and the 'ultimate recovery somewhat increased' The oil shrinkage associated with the liberation of the dissolved gas also leads to the result that as long as the connate water is immobile and the permeability ratio curve for the rock is a function only of the total liquid saturation, the stock-tank oil recovery will be less for a sand containing no connate water than for One with an Original water content as high as 3' per cent. The space voidage in the former case will be somewhat greater, but the effect of oil shrinkage will lead to smaller values for the equivalent stock-tank recovery. These calculations also give data showing how the productivity indexes for the producing wells will vary during the reservoir history. Because of decreasing Permeability to the oil and increasing oil viscosity, the productivity index will fall continuously as Production proceeds, and may finally reach values as low as 10 per cent of the initial productivity index. Introduction In a previous paper1 was developed the basic theory for the prediction of the production histories of gas-drive reservoirs. The data required for the application of this theory included the characteristics of the petroleum fluids, such as the viscosity of the Oil and gas phases, the solubility of the gas in the oil, and the oil shrinkage— all expressed as functions of the reservoir pressure—and the permeability-saturation characteristics of the producing rock, as

Jan 1, 1946

By Victor Formanek, Paul Chataignon, Maurice Rey

IT is found when floating oxidized lead ores by sulphidization, that the presence of calcium salts in the water, is usually detrimental and lowers the recovery. This effect is particularly marked in dry countries such as North Africa, where the waters often carry large amounts of calcium sulphate and where the ore may even contain gypsum. The effect of calcium salts is readily visible. Whereas in their absence cerussite is quickly stained brown and then black by sodium sulphide, in their presence the mineral remains very light in color. A similar effect is produced when barium sulphide is used as a sulphidizing agent instead of sodium sulphide. Magnesium salts have little or no effect and even tend to reduce the detrimental effect of calcium salts. A study of this phenomenon indicates that it is due to the precipitation of calcium or barium carbonate in contact with the mineral simultaneously with the formation of lead sulphide. The chemical reactions can be interpreted as: PbCO3 + Na2S = PbS + Na2CO3 CaSO4 + Na2CO3 = CaCO3 + Na2 SO4 They might also be written: PbCO3 + S = + Ca++ = PbS + CaCO3 The precipitation of calcium carbonate can be followed by the lowering of the pH with which it is accompanied. Magnesium carbonate is more soluble than calcium carbonate and usually does not precipitate under the conditions prevailing. It is interesting to note that calcium salts have no effect on anglesite (lead sulphate) because calcium sulphate is soluble, but barium salts hinder the sulphidization of anglesite because of the precipitation of barium sulphate. Remedies When calcium sulphate is present in large amounts, the softening of the water with soda ash is usually too expensive to be considered, but the precipitation of the objectionable calcium carbonate can be prevented in two different ways. One is the use of sodium hydrosulphide instead of sodium sulphide. This salt gives a lower pH than sodium sulphide and does not bring about the immediate precipitation of the calcium which remains in solution as calcium bicarbonate. The other procedure is to add ammonium salts such as the sulphate or chloride which have the property of increasing markedly the solubility of calcium carbonate. Ammonium salts have other effects such as cutting down conditioning time and accelerating flotation. They should be added to the flotation cells rather than to the ball mills. Table I. Effect of Sodium Hydrosulphide and/or Ammonium Sulphate on the Flotation of Two Oxidized Lead Ores Mibladen Ore4 La Plagne Ored Con- Tall- Con- Tail-centrate ing centrate ing Reagents Pb, Pct Pb, Pct Pb, Pct Pb, Pot Without addition of CaSO, NaS 56.0 0.68 43.4 0.65 NaSH 56.2 0.60 40.4 0.50 Na2S + (NH4)2SO4d 54.6 0.56 NaSH + (NH4)2SO4d 43.0 0.66 With addition of 40 lb/ton CaSO4 Na2S 39.4 3.85 34.0 3.21 NaSH 49.1 0.70 49.6 0.88 Na2S + (NH4)2SO4d 53.9 0.85 NaSH + (NH4)2SO4d 51.0 0.45 41.1 0.78 In table I are given the results of representative tests on two different oxidized ores. They show that the strongly detrimental effects of calcium sulphate can be offset by the two procedures outlined above. Sodium hydrosulphide is now used regularly on a mill scale on certain ores. Tests are being carried out with ammonium salts. It should be noted that malachite is subject to influences similar to cerussite. One final word of caution—when the ore is rich in primary slime, which is flocculated by the calcium salts, it may be indispensable to remove these by washing or precipitation with sodium carbonate instead of keeping them in solution by the above methods. Acknowledgment We wish to thank Socibtk Miniere & Metallurgique de Penarroya and Minerais et Metaux, for permission to publish these results.

Jan 1, 1951

By Maurice Rey, Victor Formanek, Paul Chataignon

IT is found when floating oxidized lead ores by sulphidization, that the presence of calcium salts in the water, is usually detrimental and lowers the recovery. This effect is particularly marked in dry countries such as North Africa, where the waters often carry large amounts of calcium sulphate and where the ore may even contain gypsum. The effect of calcium salts is readily visible. Whereas in their absence cerussite is quickly stained brown and then black by sodium sulphide, in their presence the mineral remains very light in color. A similar effect is produced when barium sulphide is used as a sulphidizing agent instead of sodium sulphide. Magnesium salts have little or no effect and even tend to reduce the detrimental effect of calcium salts. A study of this phenomenon indicates that it is due to the precipitation of calcium or barium carbonate in contact with the mineral simultaneously with the formation of lead sulphide. The chemical reactions can be interpreted as: PbCO3 + Na2S = PbS + Na2CO3 CaSO4 + Na2CO3 = CaCO3 + Na2 SO4 They might also be written: PbCO3 + S = + Ca++ = PbS + CaCO3 The precipitation of calcium carbonate can be followed by the lowering of the pH with which it is accompanied. Magnesium carbonate is more soluble than calcium carbonate and usually does not precipitate under the conditions prevailing. It is interesting to note that calcium salts have no effect on anglesite (lead sulphate) because calcium sulphate is soluble, but barium salts hinder the sulphidization of anglesite because of the precipitation of barium sulphate. Remedies When calcium sulphate is present in large amounts, the softening of the water with soda ash is usually too expensive to be considered, but the precipitation of the objectionable calcium carbonate can be prevented in two different ways. One is the use of sodium hydrosulphide instead of sodium sulphide. This salt gives a lower pH than sodium sulphide and does not bring about the immediate precipitation of the calcium which remains in solution as calcium bicarbonate. The other procedure is to add ammonium salts such as the sulphate or chloride which have the property of increasing markedly the solubility of calcium carbonate. Ammonium salts have other effects such as cutting down conditioning time and accelerating flotation. They should be added to the flotation cells rather than to the ball mills. Table I. Effect of Sodium Hydrosulphide and/or Ammonium Sulphate on the Flotation of Two Oxidized Lead Ores Mibladen Ore4 La Plagne Ored Con- Tall- Con- Tail-centrate ing centrate ing Reagents Pb, Pct Pb, Pct Pb, Pct Pb, Pot Without addition of CaSO, NaS 56.0 0.68 43.4 0.65 NaSH 56.2 0.60 40.4 0.50 Na2S + (NH4)2SO4d 54.6 0.56 NaSH + (NH4)2SO4d 43.0 0.66 With addition of 40 lb/ton CaSO4 Na2S 39.4 3.85 34.0 3.21 NaSH 49.1 0.70 49.6 0.88 Na2S + (NH4)2SO4d 53.9 0.85 NaSH + (NH4)2SO4d 51.0 0.45 41.1 0.78 In table I are given the results of representative tests on two different oxidized ores. They show that the strongly detrimental effects of calcium sulphate can be offset by the two procedures outlined above. Sodium hydrosulphide is now used regularly on a mill scale on certain ores. Tests are being carried out with ammonium salts. It should be noted that malachite is subject to influences similar to cerussite. One final word of caution—when the ore is rich in primary slime, which is flocculated by the calcium salts, it may be indispensable to remove these by washing or precipitation with sodium carbonate instead of keeping them in solution by the above methods. Acknowledgment We wish to thank Socibtk Miniere & Metallurgique de Penarroya and Minerais et Metaux, for permission to publish these results.

Jan 1, 1951

By I. Tamura, J. O. Brittain, T. Mura

Plain carbon steel in the cold worked or marten-sitic conditions has an internal friction peak at about 250 oC at a frequency of I cps. The influence of substitutional alloying elements on this peak was examined experimentally. The alloys were vacuum melted Co-, Cr-, Mo-, Ni-, and Si-iron (about 3 pct alloying elements), and carburized Si-iron. The internal friction of the alloys in the temperature range 30o to 500 OC was determined by means of a torsion pendulum. The 250 oC peak of the cold-worked alloys was lower in height than that observed in plain C-iron and steel. Carburized Si-iron in the quenched condition had a peak which occurred at a slightly lower temperature than in quenched plain C-steel, but the height of the peak was substantially greater than that found for the same alloy in the cold-worked condition. The activation energy of this peak for martensitic Si-steel was about 35 kcal per mol. ANNEALED a iron has two internal friction peaks due to nitrogen and carbon at 20o and 40o C respectively, at a frequency of 1 cps. If a iron is cold worked, these diffusion peaks become smaller and a new peak appears at about 250o C. The activation energy for the 250o C peak is in the range of 32 to 44 kcal per mol and has an average value of 38 kcal per mo1. It has been experimentally verified that both interstitial solute atoms and dislocations are necessary for the 250o C peak. 3 Wert has also shown that an incubation period precedes the appearance of the 250 C peak.4 With regard to quenched martensitic steel, Ke et al.596 observed two peaks at around 130o and 250o C. They conjectured that the 250o C peak in martensite was due to the same mechanism as the 250°C peak in cold worked a iron, and the 130°C peak was due to the nature of tempered martensite. They did not report the activation energy for these peaks. Chernikova7 found only one peak at about 200° C in tempered martensite. In a previous paperS we reported an activation energy of approximately 40 kcal per mol for the 250° C peak in plain carbon martensitic steels. Ichiyama8 found the activation energy to be about 36 kcal per mol for the 250° C internal friction peak in steels in the martensitic condition. The present authors proposed a theory for the 250°C peak for cold worked and quenched iron and steel.= The theory for the peak is based upon the addition of a strain-energy term to the free energy for the formation of a carbide precipitate due to the presence of a disclocation. The model employed for the relaxation process consists of the growth and solution of a carbide precipitate of critical size as a consequence of an oscillating dislocation. In this report, the influence of substitutional alloy ing elements on the 250° C peak is examined and the experimental results are discussed in terms of our theory for the peak. Most of the experimental work on the alloy steels was confined to the Fe-Si-C alloys.

Jan 1, 1962

By M. R. Tek, K. H. Coats, D. L. Katz

In natural gas storage operations, seasonal pressure fluctuations in the gas reservoir cause the water from the surrounding aquifer to flow into and out of the gas sand. The theory of unsteady-state liquid flow through porous media developed by Van Everdingen and Hurst' has been applied to predict the water movement into and out of the gar bubble for .several .postulated pressure cycles. Those cycles with as many pound-days above the original aquifer pressure as pound-days below, may cause the gas reservoir to slowly grow or shrink rather than hold to a constant volume. Applications to an actual field case study give the predicted gas reservoir monthly pressures and volumes as compared with the observed monthly pressures and volumes. INTRODUCTION A large number of natural gas storage reservoirs are bounded by or adjacent to large water-saturated formations called aquifers. The presence of these aquifers is usually evidenced by the production of water from wells delineating the gas-hearing sands. Some gas storage reservoirs in use today have been purposely developed on aquifers. Cyclic pressure variations in a gas storage reservoir cause water influx and efflux from the surrounding aquifer. This, in turn, results in a varying gas reservoir volume. The prediction of the effect of this aquifer fluid movement on the size and size variation of the gas reservoir can provide information valuable in the study of several reservoir engineering problems. The economics of gas storage operations are directly influenced by the influx of aquifer fluid into the reservoir since a shrinking storage reservoir requires increasing pressures for the storage of the same quantity of gas. Material balance and reserve or recovery calculations also obviously require knowledge of reservoir volume. Various other reservoir engineering calculations such as interpretation of well interference data, evaluation of physical characteristics of porous media, water coning, gas injection and pressure maintenance studies are typical problems where variations of volume due to edge or bottom-water encroachment becomes important. The solution of equations describing influx and efflux of water in the gas sand, application to gas storage reservoir behavior and interpretation of results are the main objectives of this paper. While this work was directed toward gas storage reservoirs, it applies equally well to repressuring of partially depleted oil reservoirs with an active water drive. MATHEMATICAL ANALYSIS As just mentioned, the cyclic pressure variations in an aquifer-surrounded gas storage reservoir cause water movement and varying gas reservoir volume. Prediction of this volume variation from the reservoir pressure cycles or the injection-production schedule requires adoption of a flow model, certain assumptions concerning reservoir and fluid properties, and formulation and solution of the governing differential equation. The flow model used in this study is a gas bubble located on a very large aquifer (Fig. la). Gas bubble radius and aquifer sand thickness are denoted by rr and h, respectively. Fig. lb shows the horizontal, radial two-dimensional porous reservoir attended by the unsteady, compressible motion of the aquifer liquid. Volume variation of the gas bubble is determined from aquifer water flow across the cylindrical boundary, r = r0

By R. E. Shinkosk

This paper outlines the development of a program for increasing the life of woolen bags used for filtering Dwight-Lloyd gases by treating the bags and gases with hydrated lime. Methods and apparatus are described for determining alkalinity of dusts, acidity and breaking strength of bag cloth. Procedure and results, based on several years of operation, are presented. DURING 1939, additional facilities were constructed in the Dwight-Lloyd Blast Furnace and Baghouse departments at the Selby, California, Plant of the American Smelting and Refining Co. In order to handle adequately the increased volume of gases from the resultant increase in production, it was necessary to increase gradually the amount of water used for cooling gases ahead of the sinter machine baghouse. As a result of this increased water cooling, the average bag life dropped from 27 months in 1939 to 14 months in 1941. (Table I). This drop in life meant an increased. bag cost, as well as lower recovery of dust and some curtailment of operation. During 1941, it was found new bags showed as high as 0.3 pct acidity* after two weeks of opera- tion and as much as 2.0 pct acidity after some months of operation. This high acidity was present in spite of the fact that free oxide or relative alkalinity of the unburned dust ran from 5 to 6 pct. In view of these circumstances, a twofold program was started in Nov. 1941.t Part one of this program consisted of vigorously dipping all new bags in a weak lime solution, containing 50 lb of hydrated lime per 50 gal of water. Part two consisted of feeding fine, dry, hydrated lime into the gas stream intake of the sinter baghouse fan. Apparatus for feeding this lime is shown in fig. 1. All baghouse chambers are shaken in rotation about once each hour. On alternate hours, the baghouse operator places 50 lb of hydrated lime (one sack) into the lime feeder, starts feeder and immediately starts the bag shaking machinery. The rate at which lime is fed is set to coincide with the approximate time necessary to shake all sinter bag-house chambers, or about 15 min. It is felt this method of lime addition is most effective for getting lime into the woolen bag fabric. The amount of lime so fed averages about 600 lb per day. The amount of lime fed per day is varied to keep a minimum relative alkalinity of 9 pct in the unburned sinter dust. A daily dust sample is taken for alkalinity by allowing dust to accumulate in a sample pipe over a 24-hr period. This sample pipe, placed in any chamber cellar, is 2 in. in diam, 4 ft long, is sealed on the inner end, and capped on the outer end. It has a 1/2 in. slot cut for 18 in. along the tip end. This slot faces upward and allows the pipe to fill gradually with dust as bags are shaken. Breaking strength of bags has, in most cases, been the deciding factor in bag replacement. Bags that normally test 100 psi breaking strength when new are replaced when they test under 35 lb. The method for determining breaking strength is shown in the description accompanying fig. 2. Since the start of the liming program in 1941, bag life has increased from 14 months to an average of over 23 months, with a consequent material decrease in bag cost per year. Acidity, as per cent sulphuric acid, may be determined by means of a Beckman pH meter as follows: From a piece of bag cloth. which has been thoroughly cleaned of dust, a 5 g sample is weighed on a balance. Cut the sample into fine pieces and place in a 400 cc beaker. Add 100 cc (measured) of distilled water and stir vigorously. Filter on suction funnel, holding cloth pulp in beaker with a stirring rod. Wash cloth sample and filter wash water four additional times, each time with 20 cc distilled water, the last time squeezing cloth pulp over funnel. Discard pulp and rinse funnel and filter paper. Pour wash solution jnto measuring graduate and make up to exactly 300 cc with distilled water. Place into clean 600 cc beaker and measure the pH on meter. The per cent acid in bag cloth is read from the following table:—

Jan 1, 1951

By S. S. Marsden, S. H. Raza

An experimental study of the flow of line-textured. aqueous foams through Pyrex tubes is described. The foams range in quality F (ratio of gas volume to total volume) from 0.70 to 0.96 and behave like pseudoplastic fluids. At lower flow rates they exhibit laminar flow ad have apparent viscosities which increase with quality, and which cover a range of 15 cp to 255 poise for tubes of 0.25- to 1.50-mm radius rj. At higher flow rates a plug-like type of flow is developed, the extent of which increases with both T and ri. When the same foams flow through either open or packed Pyrex tubes, remarkably high streaming potentials ?E are often generated. These can easily reach 50v if nonionic foaming agents are used, but are at least an order of magnitude less for ionic foaming agents. A linear relationship between ?E and the pressure differential ?p is observed; this usually extrapolates to positive values of ?p at ?E of zero. The slope of the line increases with both F and ri. An equation was derived to describe the streaming potential of non- Newtonian fluids in circular tubes and was used to correlate experimental results. The calculated < potentials are of the right order of magnitude. INTRODUCTION Foams are both unusual and inniguing in their physical properties, and have been the subject of many scientific studies. However, present knowledge of foams is still fragmentary, specific and often contradictory. Apparent viscosity of foam is the physical property of greatest interest to both rheo1ogists and engineers. sibreel reported that the apparent viscosity decreased with increasing shear rate in a manner similar to some non-Newtonian fluids. Penny and Blackman2 reported that fire-fighting foams had both a limiting shear stress and a tensile yield stress. There is little doubt that some foams at least behave like non-Newtonian fluids, and have apparent viscosities considerably higher than those of either constituent phase. The high apparent viscosity of foam with its concomittant effect on mobility ratio and sweep efficiency no doubt prompted several attempts by research groups to use foam as a displacing agent in porous media. Based on recent experience, most of these groups probably succeeded in completely blocking fluid flow in the porous media and then abandoned their efforts. Two groups apparently found the successful combination of experimental parameters at about the same time.3,7 Others have recently added to our knowledge of foam flow in porous media and its use as a displacing agent.8-13 An experimental problem encountered by Fried3 was a transient blockage of foam flow in porous media when distilled water was used to prepare the foam-producing solution. Fried surmised that this was due to an electrokinetic effect and he eliminated it by using electrolytes in preparing foaming solutions. He also measured the streaming potential of a number of foams in capillary tubes which he found to be appreciably higher than those obtained when the constituent liquid flowed under comparable conditions. This paper presents results of a more comprehensive study of the streaming potential generated by aqueous foam flowing in both open and packed Pyrex tubes. It also adds to knowledge of the rheology of these foams as deduced from their flow behavior in open tubes. APPARATUS AND PROCEDURE A diagram of the apparatus used is shown in Fig. 1. Details of its construction, testing and use are described elsewhere.14 Careful selection of materials, extreme cleanliness and rather elaborate electrical insulation and shielding were necessary to obtain reproducible results (±5 percent). Both streaming potential and streaming current were measured with an electrometer. The design of the foam generator developed for this work is novel (Fig. 2). It will produce uniform fine-textured foams at pressures sufficiently

By J. J. Pye, J. B. Drew

The surface-diffusion coefficients of Ni63 diffusing on low-index planes of nickel single crystals have been measured over the temperature range from 400° to 1000°C using a precision autoradio-gvaphic technique. The activation energy for sur -face diffusion appears to he the same on the {111}, {110), and (100) planes and amounts to -0.60 ev. This is in good agreement with the value (-0.62 ev) reported for (111) planes obtained by mass-transport techniques. In an earlier paper, the authors have reported self surface-diffusion coefficients on silver.&apos; In that work, a deviation of the diffusion coefficients from linearity in the Arrhenius plot was observed at high temperatures. This effect was attributed to a loss of tracer atoms into the bulk by volume diffusion. The purpose of the present work was to extend the technique developed to the nickel surface-diffusion system up to the point where volume diffusion would come into play. The nickel system especially lends itself to this technique due to the high efficiency of the isotope in darkening the film and the large absorption of the weak radiation (0.063 mev /3- particle of the trapped tracer atoms. The diffusion coefficients and ictivation energy for different orientations on nickel can also be compared with a report published by Blakely and Mykura 2 where a mass-transfer technique was used. This paper reports the self surface diffusion of nickel on the low-index planes over a wide temperature range using a precision autoradiographic technique. EXPERIMENTAL PROCEDURE The experiment is similar to that described in the earlier publication consisting of a radioactive needle as a source resting on a flat, smooth, chemically polished surface annealed in a dynamic dry-hydrogen atmosphere. The needle is held snugly on the surface by a holder fabricated from commercially pure nickel. The point source is made by grinding down a small-diameter nickel rod until a tip of approximately 0.1 mm radius is formed. It is then etched and the needle tip plated with the radioactive isotope Ni63. Nickel single crystals were used whose main impurities were spectrographically determined as Co 0.05, Fe 0.007, C 0.005, and S 0.003 pct. The surfaces, oriented to within 1/2 deg of the low-index planes, were carefully mechanically polished and etched in a 50-30-10-10 hot solution of acetic, nitric, sulfuric, and phosphoric acids giving a smooth polished surface. The crystals were then washed in distilled water and alcohol. Annealing temperatures ranged from 400° to 1000°C with a furnace arrangement in which the crystals could be brought up to temperature in times short compared to diffusing times. The temperature was maintained constant to ±10°C. In the temperature range investigated, the ratio of the surface-diffusion coefficient to the volume-diffusion coefficient (Ds/Du) was 105 or more and the mean bulk-penetration values, calculated from literature values of volume-diffusion coefficients,3 were on the order of 10-5 cm. After the diffusion anneal, the sintered needle is removed, and both the radius of the needle point and spot left on the crystal face are measured. The crystal is then placed on photographic film. Diffusion profiles are obtained by scanning the film with a recording microdensitometer which has a slit width of 10µ. DISCUSSION AND RESULTS The experimental conditions are such that the source strength and needle radius are essentially constant and the time is adjusted such that the mean bulk penetration is small. An exactly analogous heat-flow problem has been considered by Carslaw and Jaeger4 and the solution can be written in integral form as:

Jan 1, 1964

By E. A. Gulbransen, J. W. Hickman

Recent work by the authors1 on the oxide films formed on alloys of titanium and zirconium with nickel and copper indicated that an investigation of the oxides which form on alloys of copper and nickel might permit the relative rates of formation and diffusion of the ions of these metals to be determined. This binary system is of further interest since the alloys which form are of the solid solution type over the entire range. In addition nickel is much more oxidation resistant than copper. A comparison of the free energies of formation of the simple oxides which may occur on these alloys indicates that nickel oxide should be observed over the entire time and temperature ranges of oxidation for all alloy compositions~i In order for nickel oxide to predominate on the surface, it is necessary that nickel ions form and diffuse to the surface. Two factors which may influence the formation of nickel oxide may be the ionic sizes, where cu+ is o.,,6-,4 and ~i++ is o.,oh;, and the charges on the ions. On the basis of size, Ni++ should diffuse more readily than Cu+ while on the basis of charge, Cu+ should form more readily than Ni++. According to Uhlig2 passivity to corrosion in an alloy appears at a critical composition which for copper-nickel alloys is related to the tendency of the d band of electronic energy states to fill with elec-trans. This in in the copper- nickel alloys corresponds to 40 at. pct nickel at which composition Uhlig2 (quoting La Que&apos;s data) states that the alloy system becomes passive. Magnetic and specific heat data indicate that the d band of copper-nickel alloys is flied at 60 at. Pet Or or 4O at. pct nickel. Matt and Jones3 have shown that if, on alloying nickel with copper, one considers a copper atom to be identical with a nickel atom except for an extra electron (nickel has the configuration: 3ds 4s2 and copper: 3 dl0 4s) then at 60 at. pct copper or 40 at. pct nickel enough electrons are contributed by copper to fill the d band of the alloy. At this composition the magnetic susceptibility approximates zero according to Grew.4 Although the electron diffraction data will not determine the rates of the several reactions which may occur during oxidation, yet a determination of the structures of the surface oxides may yield information concerning the relationships, if any, which may exist between the oxide structures which form and that composition where the magnetic susceptibility approximates zero. It is quite possible, of course, that the mechanisms of the oxidation and corrosion reactions may be so widely different that no correlations between them may be possible. Survey OF Literature The metals, copper and nickel, which are face-centered cubic, form solid solution alloys which are face-centered cubic with no superstructure for all compositions. Three oxides of nickel have been investigated structurally by X ray diffraction. Of

Jan 1, 1949

By A. R. Gregory, M. R. J. Wyllie

The literature reveals that scant attention has been paid to the systematic experimental determination of the formation factors of unconsolidated porous media. No experiments appear to have been made on the effect of increasing cementation on the formation factor of an initially unconsolidated porous medium. Measurements have accordingly been made of formation factors as a function of porosity for aggregates of spheres in the porobity range 12-56 per cent and of cubes, cylinders, discs and triangular prisms in the porosity range 30-45 per cent. The results are examined in the light of the theoretical equations of Clerk Maxwell. Fricke and Slawinski. Aggregates of unconsolidated spheres and beach sands have been artificially cemented with silica and the formation factor-porosity relationships determined. A theory is outlined which seeks to explain the results obtained and which postulates that formation factor and porosity data for consolidated sandatones may be used to determine the original porosity of the unconsolidated sands from which the consolidated sandstones were derived. It follows also that for consolidated sandstones F&apos; = C where F is formation factor, fractional porosity and C and k are constants. INTRODUCTION All schemes of electric log interpretation implicitly or explicitly require the use of formation resistivity factors or? as they are now generally called. formation factors. &apos;The originator both of the concept of formation factor and of the use of formation factor data in electric log interpretation was Archie.&apos; More recently it has been suggested that formation factors are of value in describing fluid flow in porous media. A full discussion of this problem together with adequate reference.: has been given by Wyllie and Spangler.&apos; In spite of the theoretical and practical importance of the subject, there appear to be remarkably few basic formation factor data in the petroleum literature. A wider search reveals that this situation is nut confined to the literature of the to industry. at various times over the last century, attempts have been made to compute the theoretical relationship between the formation factor and porosity of simple systems of unconsolidated particle.. However. the formulae derived are not all identical and the experimental evidence adduced in support of each i; not wholly convincitng. Most of the experimental work has been carried out on system. of unconsolidated spheres, although some attention has been given to systems of spheroids and cylinder. Other simple shapes have received neither theoretical nor experimental investigation. The writers feel that a better understanding of the formation factors of complex systems of practical importance requires a firm foundation of reliable experimental data applicable to simple systems. Accordingly, they have reviewed the relevant literature and have determined experimentally the relationship between the formation factor and porosity of unconsolidated aggregates of particles of regular geometrical shapes. Some experiments have also been carried out to investigate the effect of cementation on the formation factor-porosity relationship. PREVIOUS WORK Although the term formation factor was not used. the relationship between the resistivity of a solution containing dielectric solids and the resistivity of a similar solution not containing solids is an old and important problem. The formation factor of a system of dispersed spheres as a function of porosity was considered by Clerk Maxwell. Maxwell arrived theoretically at the following expression: where F is formation factor and fractional porosity. Lord Rayleigh&apos; also examined the problem of the formation factor of systems of spheres and cylinders and concluded that for spheres Clerk Maxwell&apos;.; expression gave a value of F whirl]. at low porosities particularly, was too small. However,

Jan 1, 1953

By M. E. Wadsworth, J. S. Cho

Colloidal silica dissolved in aqueous suspensions of high surface area thoria was permitted to adsorb on the thoria surface. Silica in three forms was identified by means of infrared spectroscopy and characterized by Si-O vibrational modes in the 9.0 to 10.5 region. A band at 9. 1 was attributed to free unad-sorbed silica, while adsorbed silica resulted in bands of greatest absorbency at 9.4 or 10.3p, depending upon the character of the surface bonding. Size considerations, the amount of silica adsorbed, and the character of the infrared spectrum for strongly adsorbed silica has been interpreted as evidence of formation of chain or layer silica linkages on the thoria surface. Infrared spectroscopy has been used successfully in recent years for the identification of minerals and synthetic inorganic substances. Unfortunately, infrared studies of silicate minerals have not proceeded to the point that unambiguous assignments may be made for Si-O vibrational modes as affected by lattice substitutions and structures. Several investigators have presented infrared spectra of silicate minerals, 1,7 and their work has pointed to the potential value of infrared to identify minerals as well as to note the sensitivity of various Si-O modes to structural and substitutional environments. The strong infrared absorption of various Si-O vibrational modes makes the infrared method a very useful tool for the study of silica adsorbed on the surface of solids of high surface area. Spaepen et a18 studied the adsorption of silica on thoria surfaces from aqueous solutions of silicic acid up to 300°C. A general review of the application of the infrared method to the study of surface adsorption has been presented by Pliskin and Eischens.9 EXPERIMENTAL Du Pont colloidal silica, Ludox SM (16.36% SiO2) and Ludox HS (30.25% SiO2) in aqueous suspensions were used as the sources of silica for the adsorption studies. The thoria used for this study was supplied by Union Carbide Nuclear Co.* under the preparation *Chemical Technology Div., Union Carbide Nuclear Co., Oak Ridge, Tenn. . number D-16-650. It was fired from the oxalate at 650°C and had a surface area of 32.4 ± 0.6 m2 per g measured by nitrogen adsorption. The chemical analysis of the D-16 thoria in ppm was: CO3= (5500), NO3- (< 4), so4±PO4 (100), and Si (20). The Ludox silica was placed in contact with the thoria in a two-liter autoclave10 between room temperature and 250°C. Slurries were made by placing 1 g of thoria in one liter of distilled water containing various additions of the silica sol. Samples were removed from the autoclave by opening a small stainless steel valve. Extracted samples were centrifuged for 10 min and excess solution was decanted. The sample was then dried for several hours at 80°C. Infrared samples were prepared by thoroughly mixing 10 mg of the thoria with 400 mg KBr in polystyrene vials on a "Wig-L-Bug" amalgamator. Pellets were pressed under vacuum and had final dimensions of 1.2 cm diam and 1.2 mm thick. Infrared spectra were

Jan 1, 1963

By O. N. Carlson, H. A. Wilhelm, J. M. Dickinson

On the basis of microscopic studies, melting-point observations, and X-ray analyses, a phase diagram is proposed for the Cb-V system. A complete series of solid solutions is formed with a minimum in the solidus at 1810°C near 35 wt pct Cb. No compounds or intermediate phases were found in the system above 650°C. THERE is an ever increasing need for better structural metals and alloys for use in nuclear reactors. In addition to the normal properties of engineering structural materials, such as high temperature strength, resistance to corrosion, and fabric-abil~ty, the nuclear properties of the material must be considered. In a nuclear reactor it is important to conserve neutrons, so a material which removes these neutrons from the reaction excessively is considered to have unfavorable nuclear properties. In nuclear-reactor design the engineer must have nuclear as well as other data available on alloys in order to make a wise selection of materials. Due to the fact that many of the common structural materials have undesirable nuclear properties, it is vital that new alloys of metals having more favorable nuclear properties be investigated. Columbium and vanadium are both high melting metals, both exhibit resistance to chemical attack, and no great difficulty is encountered in fabricating them into desired shapes. With proper treatment both metals can be cold rolled extensively without failure. In addition they have desirable nuclear properties for certain types of reactors. Therefore, the alloys of columbium and vanadium should be of interest in the atomic energy program. Since an alloy-development program is enhanced by a knowledge of the phase equilibria of the components, this investigation was undertaken to establish the phase diagram for the Cb-V system. According to the Hume-Rothery rules for alloying,&apos; the chemical similarity, crystal structure, and atomic-size factor are favorable for a complete series of solid solutions for this system. Both elements are in the same family of group V of the periodic table and thus are quite similar chemically. The crystal structures of columbium and vanadium are compatible for extensive solid solubility, since both have body-centered-cubic structures. The atomic diameters of columbium and vanadium are 2.85 and 2.62Å, respectively. This difference of slightly more than 8 pct is well within the 15 pct maximum difference allowed for extensive solid solubility. Experimental Procedures Source of Materials: Columbium powder and sheet trimmings were obtained from the Fansteel Metallurgical Corp. According to the manufacturer the metal contains less than 1 pct impurity. An analysis of the metal showed approximately 1800 ppm C in the powder while the sheet trimmings contained less than 500 ppm C. Spectrographic analysis showed minor amounts of Ca, Cr, Fe, Mn, Si, Ti, V, and Zr in both forms of the columbium. No commercial source of vanadium having the ductility and purity desired was available to the authors at the beginning of this investigation. As a result, all of the vanadium used in this study was prepared by the bomb reduction of vanadium pen-toxide with calcium employing the method reported by Long.&apos; Yields of massive vanadium normally were about 80 pct. Chemical analysis of the vanadium prepared in this manner showed the presence of 200 to 500 ppm N and 800 to 1000 ppm C. Minor amounts of Ca, Fe, Mn, Si, Zr, Cr, and Cb were detected by spectrographic analysis. This vanadium metal was ductile and was cold rolled into 5 mil sheet. Annealing was not necessary during this rolling and the metal retained its cold-rolling characteristic after are-melting. Preparation of Alloys: The Cb-V alloys were prepared by melting pieces of vanadium sheet togethel-with columbium in the form of sheet or pellets of powder. The melting was carried out under argon in conventional arc-melting equipment employing a tungsten electrode and a water-cooled copper crucible. Each alloy was remelted three or four times, inverting the alloy after each melting in order to assure complete mixing. Alloys normally were obtained as round flat disks, weighed approximately 70 grams, and had roughly the shape of a disk 1 1/2 in. in diameter and 1/4 in. thick. Half of each alloy

Jan 1, 1955

By D. P. Kedzie, R. A. Dodd

The fatigue strength, fatigue hardening, and effect of fatigue deformation on subsequent age hardening of supersaturated Al(Cu) solid solutions have been determined as functions of alloy composition and temperature. The fatigue strength/tensile strength ratios, determined at 150°, 25°, and -195°C, decreased with increase in alloy content for all temperatures, but the F.S./U.T.S. ratios at -195°C decreased much more rapidly than did the ratios for 150" and 25°C. This suggested that strain aging and/or age hardening occurred during tests at higher temperatures. Additionally, the F.S./U.T.S. ratios at 150°C exceeded those at 25°C for all compositions, indicating greater strengthening during fatigue at 150°C. The effect of fatigue or tensile deformation at 150°, 25°, and -195°C on subsequent age hardening showed that the deformatzon increased the rate of precipitation and indicated that mechanically produced vacancies were probably formed during deformation. Fatigue hardening was studied at 150°, 25°, and -195°C, and the effect of room-temperature rests after 10 and 100 cycles was examined. The results confirmed that strain aging occurred at the higher temperatures . DURING the last 15 years various mechanisms of fatigue crack nucleation and growth based on dislocation and vacancy interactions, operating singly or collectively, have been proposed. The probable consensus of present opinion is that the fatigue process in pure metals essentially involves dislocation interactions, and that vacancies formed by such interactions play a minor or inconsequential role. However, there is some evidence that age-hardened alloys tend to overage during fatigue, probably by local vacancy-enhanced diffusion, and strain aging also might be important in selected cases. Furthermore, it has been established that the behavior of quenched-in vacancies in Al(Cu) and other solid solutions is composition-sensitive. Therefore, it seemed worthwhile to investigate various aspects of the fatigue behavior of supersaturated Al(Cu) alloys and to examine the results in terms of vacancy-enhanced effects. EXPERIMENTAL The alloys used in this investigation were prepared from 99.994 wt pet A1 and OFHC copper, the latter containing 0.04 pet 0 as the principal impurity. Six alloys were made, containing, by actual analysis, 0.58, 0.96, 1.96, 2.85, 4.45, and 5.51 wt pet Cu. They were induction-melted in air in graphite crucibles, cast as 7-in. by 7/8-in.-diam rods in graphite molds, and hot-rolled to 5/8 in. diam. The experimental work was a three-part program involving the determination of a) 10&apos; cycle fatigue stresses as a function of alloy composition and temperature; b) the effect of fatigue deformation on subsequent aging of the supersaturated alloys; and c) fatigue hardening as a function of alloy composition and temperature. For determining the 105 cycle fatigue stresses, a portion of the 5/8-in. stock was machined into Krouse rotating cantilever beam fatigue specimens, 2 in. in length by 1/4 in. minimum diameter. These were tested at +150°, 25°, and -195°C (liquid nitrogen) on a Krouse high-speed machine, with special weights to provide lower -than-normal loading ranges, this being necessitated by the small load requirements of those alloys of lower copper content. For the high-temperature tests a small resistance heater was designed to clear the collets and fit in between the chucks, while for the low-temperature tests a hollow nylon cylinder was used, having closed ends drilled to pass a fatigue specimen, and positioned similarly to the heater. A plexiglass container completely enclosed the fatigue machine; rubber gloves fixed to ports in the walls enabled the machine to be operated from outside the box. A tray of magnesium perchlorate dried the air in the container and prevented both atmospheric corrosion fatigue at room and elevated temperatures and troublesome ice build-up at low temperatures. A total of ten to fifteen specimens was used for each combination of alloy composition and temperature. The result of each test was plotted on a standard S-N diagram, and the next stress was selected on the basis of the trends indicated by previous specimens. In this manner a small portion of the S-N diagram was constructed, and the 105 cycle fatigue stress obtained. Small Krouse fatigue specimens were also used to study the effect of cyclic prestrain on subsequent

Jan 1, 1964

By W. R. McMillan, E. A. Gulbransen, K. F. Andrew

Annealed and electrolytically polished pure iron was oxidized between 650° and 850°C at oxygen pressures of 0.1 to 2 microns Hg. Electron optical studies showed that oxidation occurs discontinuously over the surface and orients crystallographically with the substructure of the metal grains. IN previous papers&apos;, " the authors described the formation of oxide films on high purity iron at norrnal pressures and at temperatures below 300°C. Electron optical studies showed that the oxide film consisted of many individual oxide crystallites having nearly random orientation on each metal grain. The size of the oxide crystallites depended upon the temperature and upon the extent of oxidation. In these studies the crystallite size varied from 250 to 1800Å. Recently, Bardolle and Benard3 studied the crystal habit of the oxide crystallites formed on Armco iron between 650" and 850°C in what might be defined as a "poor vacuum" atmosphere in which the pressure could be varied between 10-1 to 10-3 mm of Hg. Their results showed that, at 850°C and a vacuum of l0-2 to 10-3 mm of Hg, a few well oriented nuclei of oxide were formed on the metal grains. At pressures between l0-1 to 10-2 mm of Hg, many oriented nuclei of oxide were formed. At 750°C the metal surface was bright although an oxide was present. Unfortunately, the extent of oxidation, the oxygen pressure, and the influence of carbon in the metal were not determined, and it was, therefore, impossible to give a theoretical interpretation to the very beautiful light micrographs. Bardolle and Benard used the light microscope to determine the crystal habit of the oxide and X-ray diffraction to determine the orientation relationships. It was believed that the use of controlled amounts of oxidation at high temperature is an important method for the study of the initial stages of oxidation of metals at high temperatures. This paper presents the first results of an electron optical study of this problem for iron. Since the electron microscope and electron-diffraction methods using stripped oxide films and surface replicas yielded finer details of the crystal habit and structure of oxide film, these methods were used in this study. Two types of iron having different carbon contents were used: Puron and Armco* iron. A com- parison of the results for these materials made possible a study of the influence of the surface oxide-carbon reaction on the initial stages of oxidation. Before discussing the experimental results, the thermodynamic equilibria of 1—the oxidation reaction, 2—the surface oxide-carbon reaction, and 3—the solubility of oxygen in iron, will be considered. In addition, it is important to consider the kinetic-theory predictions on the rate of collision of oxygen with the iron surface. Thermodynamic and Kinetic-Theory Calculations Only one type of reaction is usually considered in the oxidation of iron at normal pressures at high temperatures: Fe(a) + 1/2O2(g) = FeO(s) 3Fe(a) + 2O2(g) = Fe3O4(s) 2Fe(a) +3/2O2(g) = a-Fe2O3(s). [11 However, two other reactions may occur: FeO(s) + C (solid solution in Fe) + Fe(a) + CO(g) [2] and Fe(a) + O2(g) e Fe(a) (solid solution of 0,). [ 3 ] In addition, iron will react with the higher oxide at certain temperatures to form FeO. For a complete interpretation of the oxidation of iron at low pres-

Jan 1, 1955

By D. W. Lillie

Property changes produced by irradiation of Al-0.4 wt pet Li alloys at 270°C to a burnup of 0.155 pct of all atoms are described. Metallographic evidence of the formation of internal pores and the concentration of gases at grain boundaries is shown. Some conclusions are drawn concerning desirable property characteristics for resistance to radiation swelling. RADIATION induced gases are of great practical significance in nuclear engineering since they affect the behavior of fissionable fuels and control rod materials. In fuels the gases xenon and krypton are the major offenders, producing serious em-brittlement at low irradiation temperatures and gross swelling at high irradiation temperatures. In control rod materials containing boron, the production of helium from the B 10 (n, 0) Li7 reaction produces similar reactions. Since the trend in reactor design is for operation at increasingly high temperatures, the phenomenon of high-temperature swelling by these induced gases has been the subject of considerable recent study. The objectives of such work have been to obtain engineering information and quantitative data on the effects in specific fuels and control rods and to obtain further understanding of the mechanism of gas damage so that more intelligent corrective steps can be made. Analogue systems are often useful for the latter objective and in this case alloys containing Li are suitable in view of the reaction Li6 + n - He4 + He 4. Both products are gases and no serious radiation problem is produced to complicate post irradiation handling. Considerable caution is, of course, necessary in any operations in which tritium release is possible. The present program was initiated to provide a detailed study of the effects of He + H3 generated by irradiation in Li-bearing alloys. The results on Mg-Li solid solution alloys have been reported elsewhere.&apos; This paper presents information on irradiation of an A1-0.4 pct Li alloy in which the lithium was present as highly enriched Li6 in order to maximize the radiation burnup. PROCEDURE The original intent in the selection of the alloy content was to have Li in solid solution at all temperatures. It became apparent, however, that the solubility of Li in A1 was substantially lower than that reported in the original reference consulted.&apos; This was substantiated by the work of Nowak3 who reported solubility of less than 0.01 pct Li at 100°C. The alloy used (0.4 pct Li) was therefore two phase at room temperature. The irradiation temperature was estimated to be about 270° C, and at that temperature all or almost all of the Li would be soluble under equilibrium conditions. The precise state of the alloy under irradiation is not known. Aluminum used was 99.996 pct purity (max. impurities Si 0.002 pct, Mg 0.001 pct, Fe 0.0005 pct, and Cu 0.0005 pct). The lithium was obtained from Oak Ridge National Laboratory; 96.1 + 0.1 pct of the lithium was ~i~, but chemical impurities were 0.25 pct Ca, 0.05 pct Fe, 0.02 pct Na, 0.02 pct Cu, 0.01 pct Sr, and 0.01 pct Ba. The two metals were melted in an alumina crucible in a resistance heated furnace in an atmosphere of helium. A l-in. ingot was cast in an iron mold, machined to %-in. round, hot swaged (350°C) to 0.080-in. round, cold swaged to 0.060-in. round, and finally cold drawn to 0.040-in. round. Four-inch wire lengths were then cut and annealed at 250°C for 4 hr in helium. Irradiations were carried out in a VT hole in the Argonne CP-5 reactor. The wires were held in an aluminum jig, sixteen wires per capsule, and were maintained in a helium atmosphere for maximum heat transfer. Eight wires of pure aluminum and eight of Al-Li alloy were placed in each capsule. Argome Laboratory statements were accepted as establishing flux and exposure levels. These were based on previous experience and flux map-

Jan 1, 1961

By J. E. Kirby, L. B. Schnitz, H. E. Stamm III

The production history of a gas-cap-drive reservoir was reproduced by calculations, and predictions were made for operations under primary depletion, pressure maintenance by gas injection, and pressure maintenance by water injection. Although the general equations used in this analysis have been available to the industry for a number of years, the procedure has not had extensive use because of the unusually large time requirements for developing numerical results. The use of digital computing equipment for processing the calculation work generally reduced the time usually required for an extensive analysis of this type. Calculation techniques were developed that were equally applicable for solving the gar-cap-drive problem, whether digital computing equipment or desk calculators are used. These techniques are outlined and the procedures are illustrated by the field problem. INTRODUCTION The analysis of a gas-cap-drive res-ervoir requires the simultaneous solution of the volumetric balance1 and displacement2 equations. Equations for the analysis of combination drive reservoirs were presented by Wooddy and Moscrip in a previous paper." These equations permit the evaluation of possible operating processes in order that the optimum method of pro- ducing the reservoir may be determined. Although these equations represent the basic technique known for the detailed analysis of gas-capdrive reservoirs, the calculation processes are quite complex and time consuming; hence, the procedure has not had extensive use. Recently, these equations were satisfactorily applied with the aid of a digital computer to investigate the completion history of a gas-cap-drive reservoir under the following operating methods: primary depletion, depletion under full pressure maintenance by gas injection, and depletion under full pressure maintenance by water injection. The use of digital computing equipment substantially reduced the time necessary for making the computations and permitted a more concise investigation of the reservoir&apos;s depletion history than was possible when the desk calculator was used. This more rigorous investigation also permitted a detailed analysis of the calculation procedures necessary in the solution of the gas-cap-drive problem. This paper outlines the procedure of analyzing a gas-cap-drive reservoir and illustrates the process with an actual field problem. CHARACTERISTICS OF THE RESERVOIR ANALYZED Production in the field under consideration is obtained from a Frio formation of Oligocene age at a depth of approximately 5,740 ft sub-sea. The accumulation, which is bounded on the flanks by strand lines, is semirectangular in shape with gentle narrowing toward the base of the structure as shown by the sand limits and completion map in Fig. 1. The reservoir development indicates an oil column of 319 ft that is defined by gas-oil and oil-water contacts. Basic data relative to the characteristics of the formation and the formation fluids are shown in Table 1. The formation exhibits considerable shaliaess and thinning near the oil-water contact indicating that little or no benefit will be received from natural water influx. This is substantiated by the fact that wells completed near the original oil-water contact have not produced water. A sizeable gas cap is present; however, considerable faulting is in evidence, and it is believed that only a limited portion of the gas cap will aid in oil displacement and pressure maintenance. Although the reservoir was in an initial stage of depletion, it was believed that as soon as the equilibrium gas saturation in the oil zone was exceeded and a number of upstructure wells were invaded by the expanding gas cap, the resulting high gas-oil ratio production would rapidly dissipate the natural drive and lead to limited oil recoveries and inefficient reservoir operation. This analysis evaluated the effectiveness of the primary displacement mechanism and

Jan 1, 1958

By Eric J. Jette

In thermodynamic studies of gas mixtures and liquid solutions, the respective problems have been greatly simplified by the use of two general limiting laws; Dalton&apos;s law of partial pressures and Raoult&apos;s law. It is generally recognized that these two laws have a common physical and chemical basis which may be stated as follows: The law may be expected to be strictly valid only for systems in which the different molecular species are so nearly the same in their physical and chemical properties that a molecule is unable to recognize whether its neighbor is one of its own kind or not. Positive and negative deviations from both laws are frequently encountered. The occurrence of negative deviations in both gases and liquids is ascribed to attractive forces which, in the more extreme cases, may be so large in magnitude and so highly specific in character as to lead to the formation of definite chemical compounds. Positive deviations in gases develop only at higher pressures where the intermolecular distances are so decreased that factors such as molecular size and shape, which ultimately become concerned with the repulsive forces, come into play. In liquid mixtures, positive deviations are generally considered as due to differences in the symmetrical (nonpolar) fields of force which surround the various molecular species and which are supposed to give rise to a tendency for each species to associate with its own kind and thus to "squeeze out" the other variety. Fundamentally this also involves the existence of repulsive forces between unlike molecules. If the fields of force are unsymmetrical (polar) a tendency towards compound formation results. Stated in more exact terms, strict adherence to the law means that the "activity" (or any adequate measure of this quantity) of a component in a mixture is a linear function of the molecular (or atomic) fraction of that component between the mol fraction limits 0 and 1. A positive deviation means that the activity is greater than corresponds to the law; a negative deviation that the activity is less. The deviations from Raoult&apos;s law have been fully discussed by many writers1

Jan 1, 1934

By G. L. F. Powell, L. M. Hogan

Large degrees of undercooling have been produced in bulk samples, 400 g, of copper and Cu-O alloys by melting in a slag of commercial soda-lime glass. The maximum degrees of undercooling obtained for copper, hypoeutectic and hypereutectic Cu-O alloy samples were 208°, 218°, and 97°C, respectively. No grain refinement as a function of undercooling was observed for the pure copper samples although in Cu-O alloys there was a marked decrease in grain size at degrees of undercooling greater than 150°C. The grain size change is the result of recrystallization during or immediately following the freezing process. THE first report of a large degree of undercooling in a bulk sample of metal was that by Bardenheuer and Bleckmann1 who produced 258°C undercooling in a 150-g sample of iron by melting in a glass slag. Garbeck2 extended this technique to nickel, cobalt, and copper and reported maximum undercoolings of 220°, 230°, and 60°C, respectively, while Fehling and schei13 undercooled a large number of metals as 2 to 20-g samples slagged in glass. Subsequently, one of the authors4 applied the glass slag technique to bulk samples, -500 g, of silver and obtained a maximum undercooling of 250°C. It was concluded that impurities which normally catalyze nucleation at small degrees of undercooling could be removed by oxidation and solution in a glass slag. Large undercooling of bulk samples of iron, nickel, cobalt, and their alloys has also been reported by other authors.8&apos;10&apos;11 Of the five metals mentioned above, copper was the only one which had not been undercooled by at least 200°C in the form of a large bulk sample. Since copper exhibits high liquid solubility for oxygen, as do iron, nickel, cobalt, and silver, it was considered that a much larger degree of undercooling than that observed by Garbeck for copper should be possible. This was obtained by modification of the technique applied previously to silver and the results of undercooling experiments with 400-g samples of copper form part of this paper. The undercooling behavior of Cu-O alloys was also studied, and the influence of a small oxygen content on the grain structure of undercooled copper was observed metallographically. 1) EXPERIMENTAL The copper was obtained as oxygen-free high-conductivity copper, the major impurities of which were Fe 0.01 pct, Zn 0.015 pct, and Si 0.02 pct. Melting was carried out in an open-ended vertical cylindrical fur- nace wound with Kanthal wire. A schematic diagram of the experimental setup is shown in Fig. 1. The undercooling technique previously used for silver4 involved preoxidation of the samples by melting the silver in contact with the atmosphere. The oxygen content was subsequently reduced by freezing the samples under a glass slag, whereby the oxygen was evolved as gas and the continuous glass cover over the surface prevented re-solution of oxygen when the sample was remelted. Since oxygen is not released as gas when a Cu-O alloy is frozen. the technique had to be modified for use with copper. Also, copper is oxidized during fire refining so that an in situ preoxidation treatment of the samples was not considered to be a prerequisite to large undercooling. Experimentation proved that this supposition was correct. Instead, it was found that care was necessary during melting of the copper under glass to minimize oxygen pickup. Samples of copper weighing 400 g were prepared by adding 50-g pieces of copper to a vitreous silica crucible partly filled with soda lime glass at a temperature of approximately 1000°C. Each piece was quickly immersed in the glass slag and held at this temperature until the oxide coating on the surface decomposed. This change was easily observed and coincided with the formation of gas bubbles in the glass

Jan 1, 1969

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

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