Search Documents

By P. G. Meikle, C. T. Holland

The effect of the condition of friction upon the ultimate unit compressive strength as determined between the compressed faces of precisely machined coal specimens and specially prepared testing platens has been examined. Reductions of 30 to 40% in the ultimate unit compressive strength of these model coal pillars have been verified when lubricants areused between the faces as compared to non-lubricated faces. Increase in unit compressive strength due to confinement at each of three levels of friction appears to follow generally the widely accepted "1/2 power law", i.e. S - C (L/T)1/2 where S is the unit ultimate strength, C is a constant dependent on the coal being tested and the coefficient of friction acting at the compression surfaces and L and T are the least width and thickness of the specimen, respectively. Results suggest that the surface condition of both the specimen and platens should be mentioned in detail in conjunction with compression testing, for although the term (L/T)1/2 is adequate for all practical purposes, the coefficient C is directly dependent on the condition of interfacial friction, all else remaining constant. One phase of the investigation of the confining effect of friction on the ultimate unit compressive strength of small coal specimens has recently been completed at the School of Mines of West Virginia University. The efficient control of ground is of fundamental importance to successful mining operations and the design of support pillars has always presented the mining engineer with a difficult and intriguing engineering problem. Since about the turn of the century, many investigators have performed extensive laboratory tests on small coal specimens in efforts to solve some of the perplexing problems of strata control and pillar action in coal mines. A detailed review of this work would serve no purpose here, but a brief mention of the more significant contributions is in order. Muller16 was perhaps the first to illustrate the effect of confining pressure on the strength of coal specimens in compression. His work has recently been more carefully expanded by the experiments of Murrell17 in Great Britain. Rice and Enzian l9 observed that as the size of cubical test specimens increased, the unit compressive strength decreased. These observations were later supplemented by Lawall and Holland, 13 by Greenwald4,5 et a1 at the U. S. Bureau of Mines, by Gaddy3 as well as by Steart,22 Millard15 et a1 and Evans2 et al. Gaddy3 seems to have been first to quantify this relationship with his equation: S = K/vd , where S - ultimate unit compressive strength (psi), K = a constant dependent on the coal, and D =edge dimension (in inches). Holland7 concluded that the strength of coal as affected by lateral dimensions and thickness could be expressed as where K = a coefficient dependent upon the coal, L = least width of pillar (in.) and T = thickness of pillar (in.). Data supplied by this formula have been compared to the estimated stresses in pillars that had been standing for long periods in coal mines and in the cases so far studied the results indicate that the formula predictions are sufficiently accurate for engineering purposes (Holland11 ). This equation therefore is believed to be the simplest and most accurate tool available at the present time for approximating the strength of a coal pillar (over a specified range of L/T) in situ. However, Holland7,9 has indicated that friction between the coal pillars and the top and bottom rock

Jan 1, 1965

By Guy Borden, Michael J. Rzasa

Laboratory data on bubble point pressures and reservoir volume factors have been correlated as functions of solution gas-oil ratio, calculated gas gravity of the pentanes-and-lighter fraction of the entire fluid, differential residual oil gravity, and reservoir temperatures. INTRODUCTION Several correlations of crude oil properties have appeared in the literature. D. L. Katz4 in 1942 presented five methods of predicting oil shrinkage. these being of decreasing accuracy for decreasing amounts of information available. M. B. Standing- in 1947 published three correlations of laboratory flash vaporization data of California crudes. From values of GOR (gas-oil ratio), gas gravity, liquid gravity, and temperature, his correlations will predict bubble point pressure, formation volumes of bubble point liquids, and two-phase formation volumes. Curtis and Brinkley2 in 1949 presented several correlations. From the gas-oil ratio, an approximation of reservoir volume factor and barrels of condensate recoverable per barrel of reservoir space may be obtained; along with liquid gravity and reservoir temperature, the GOR will allow prediction of bubble point pressure. These last cor- relations seem to be more qualitative than quantitative. Generally, laboratory bottom hole sample tests furnish information on solution gas-oil ratios, residual oil gravities. bubble point pressures, viscosities of oils, liquid shrinkages, and occasionally gas gravities. Each of these data has its own applications and use in reservoir engineering calculations. The particular uses of correlated bottom hole sample data are found in (1 Providing a basis for obtaining estimates of formation crude properties in fields where bottom hole sampling is impractical or impossible. (2) Greatly reducing the time in obtaining the desired information. (3) Determining the applicability of the results from various bottom hole samples to particular field problems. (4) Avoiding, in many cases, the uncertainties of sampling by replacing it with an element over which greater control can be exercised. (5) Permitting use of preliminary field data in application of production procedures before a bottom hole sample can be obtained and analyzed in the laboratory. (6) Serving as a check on data which may appear out of line. (7) Estimating for a particular type crude the appropriate equilibrium constants by working backward from the bubble point pressure. (8) Estimating original or other past history properties of reservoirs that were not sampled in the past. PROCEDURE Application of the published correlations4,6 to Stanolind laboratory data indicated that the general scheme presented by Standing" could give dedrable results if changes were made in parameter positions and scales. The correlation curves were drawn with all the variables having consistent gradations except the temperature increments which were drawn in to best fit the data. The variables from available Stanolind laboratory data are defined below: (1) Gas-oil Ratio: Gas is liberated at reservoir temperature by differential vaporization (or rather by a series of flashes, approaching differential vaporization) and measured at atmospheric pressure and temperature, at which the compressibility factor is assumed to be unity. The oil is the residual liquid remaining after the pressure has been reduced to atmospheric. For the gas-oil ratio both volumes are corrected to standard conditions of 14.7 psia and 60°F. (2) Gas Gravity: It was decided to arbitrarily divide the hydrocarbons of the entire bottom hole sample into pentanes-and-lighter and hex-anes-and-heavier, and use a calculated gas gravity of the pentanes-and-lighter for a correlating variable. (Sample calculation is shown in Table 111.) (3) Liquid Gravity: This is the API gravity of the residual liquid from the differential vaporization. The gravity is measured at room temperature and corrected to 60°F.

Jan 1, 1950

By Granville S. Borden

MANY complaints have been filed by mine prospectors, explorers and producers about inequities in mining taxation. They have insisted upon improvement, and in the past 35 years there has been some progress. Many of the best taxation provisions made to date are those contained in the Revenue Act of 1951. Mining economics presents special problems, a fact even recognized by many textbooks. Certain provisions have been made to take care of mining taxation rights, but these have been made by Congressional act and do not represent basic, constitutional, guarantees. Despite special provision, in the form of depletion allowances, there remain injustices. The industry faces the problem of gaining adequate income tax allowances for exploration and development expenditures, operating losses, and bonus and subsidy payments. Problems in mine taxation arose soon after the adoption of the 16th Amendment. This amendment gave Congress the power to impose income taxes. Supreme Court decisions against a number of mining companies held that Congress, under its power to tax income, could tax net operating revenue before depletion allowance was deducted. The court also ruled that the right to deduct depletion is a matter of Congressional permission. There is no statute provision permitting mining stockholders to recover their investment out of dividends while the mine is being operated. The Supreme Court did not recognize that part of mining proceeds are actually amortization of the cost of mineral rights. Both these omissions stem from failure to recognize that by its very nature mining consumes its mineral rights and therefore that a part of gross proceeds are not income but recovery of capital. Suppose a Congress unfriendly to mining interests repeals the provisions that allow deductions for depletion. Under such a situation assume that a corporation buys a mineral deposit for $300,000. They mine the deposit and earn an operating profit of $500,000. At the current rates they would pay a Federal Income tax of in excess of $250,000. This venture results in a loss of $50,000 as shown in the table below. For the example given the absence of depletion provisions in income tax law would turn ore into waste rock, because the word "ore" connotes commercially mineable rock. [ ] In retrospect the conclusion is clear: If the Supreme Court had recognized that a part of mining proceeds represent in reality sales proceeds of mineral rights, then profits derived from such sales would have been taxed at the favorable long-term capital rates, instead of being taxed at the high normal, surtax and excess profits rates. It is a question whether percentage depletion allowances offset this tax situation. History The history of depletion allowances has been one of general improvement. The original act did not include depletion, but an allowance at the rate of 5 pct of gross proceeds was provided in 1913. The Acts of 1916 and 1917 allowed depletion upon cost of property, and in 1918 the law was further changed to permit depletion upon cost, or market value as of date of discovery. In 1932 metal, coal and sulphur properties were permitted to elect depletion in percentages of gross income, limited to 50 pct of net income, in lieu of cost depletion. Shortly thereafter the rule was changed and percentage depletion is now allowable in any case where it is higher than cost depletion. Since 1942 percentage depletion allowances have been extended to cover many nonmetals, and in 1951 the rate for coal mines was increased from 5 to 10 pct of gross income. It is interesting to note that these allowances have been retained, and even

Jan 1, 1952

By T. R. Ingraham

When aluminum sulfate was decomposed in an evacuated apparatus and the sample cooled, Q alumina was found. The thermodynamics of the decomposition reaction were resolved by measurements of the total pressure of SO3, SO2, and O2 generated over the sample at various temperatures. The results were correlated with other published data to calculate the rmges of thermodynamic stability of the various compounds involved. The thermal decomposition of aluminum sulfate was studied by Wijhler, Pluddeman, and Wöhler1 and by Warner and Ingraham.2 The results of both studies were critically reviewed by Kellogg,3 who observed that "each set of data shows considerable scatter and strong curvature on a plot of log pso3 vs 1/T".3 He agreed with Warner and Ingraham2 that both the discrepancy in the two sets of results and their curvature may arise from the indefinite properties of the particular form of A12O3 which is formed during the decomposition. While studying the thermal stability of other sulfate systems, we had noted the presence of nonequilibrium oxide properties in both the cobalt4 and manganese5 systems, but neither had shown any appreciable curvature on a logPSO, vs 1/T relationship, and the data from both systems agreed reasonably well with prior work. It therefore seemed reasonable to reexamine the technique which had been used in our earlier work on aluminum sulfate.2 The principal difference was that a tightly compacted pellet had been used for the earlier aluminum sulfate work, whereas a loose powder, contained in a platinum boat, was used for our work on other sulfate systems. In addition, a better technique had been developed for dehydrating the samples. From work done on the nickel system,6 we had become aware of the highly protective nature of some sulfate layers on their oxjdes, and hence of the necessity of providing a large surface area of P3e sample when equilibrium was to be attained in a reasonable time. That observation suggested that perhaps a steady state, rather than equilibrium, had been reached in the earlier work on aluminum sulfate. This paper reports a repetition of the experimental work with large, finely divided samples of anhydrous aluminum sulfate. APPARATUS AND PROCEDURE The apparatus used in the experimental work was described in earlier papers.2,7 The source material was Anachemia reagent-grade Al2(sO4)3 - 18H2O which was first partly dehydrated by heating at 400°C for 48 hr and then completely dehydrated by repeated baking of an aluminum sulfate-oleum slush to dryness at 450°C. ? Alas was prepared by decomposing some of the sulfate in a muffle furnace at 800° C overnight. Approximately equal amounts of the 7) alumina and aluminum sulfate were used in the experiments. The finely divided samples were handled under dry box conditions, and in addition some of the aluminum sulfate was decomposed in the apparatus while it was under vacuum, to ensure its complete dehydration prior to sealing the apparatus. The samples were heated to a series of constant temperatures on both increasing- and decreasing-tem-perature programs. Two to three days were required to obtain steady readings of the total decomposition pressure on an increasing-temperature program, and 5 to 6 days on a decreasing-temperature program. Three samples were examined on both programs. Within experimental error, the results for each sample and each program could not be distinguished.

Jan 1, 1969

By Paul F. Sikora, John W. Weeton, Philip A. Clarkin

Effects of variations in carbon content and micro-structure on the stress-rupture properties of a liquid-phase-sirztered, Izeut-resistunt alloy were studied. Using the powder rnetallurgy technique, specimens could be produced having stress-rupture profierties equivalent to or better than those of the alloy in wrought form. In addition, sintered specimens were easily hot worked even when carbon contents were as high as I pct. Correlations between stress-rupture properties and micro structure were obtained in some cases. IT is well known that many alloys used in high-temperature service are quite complex both from the standpoint of chemical composition and microstruc-ture. During the development of these alloys their complexity has in some cases led to treatment schedules which, while practical, often do not fully utilize available alloying elements. For example, in some alloys, massive carbides that form during ingot practice may never be dissolved by reasonable solution treatments and therefore are unavailable for subsequent precipitation hardening, Ref. 1. This investigation was undertaken to determine whether the powder metallurgy method of consolidation might provide distinct advantages over the conventional melting method. First, segregation of alloying elements might be eliminated or reduced and thus a more efficient use of these elements obtained. Second, the structures produced in powder specimens might differ considerably from those normally seen in castings, which would perhaps result in better response to subsequent mechanical or thermal treatments. If these results could be achieved, it was felt that both forgeability and stress-rupture properties might be improved. The alloy used in this investigation was S-816 produced by powder metallurgy methods. The influence of carbon content over the range of 0.1 to 1.0 wt pct was studied for the alloy in the as-sintered condition and after three additional treatments involving heat treatment, hot swaging, and combinations thereof. EXPERIMENTAL Standard S-816 has the composition shown in Table I. It is a cobalt-base alloy having a nominal carbon content of 0.4 pct and is classified as a precipitation-hardening alloy. This alloy was chosen because it is representative of a class of heat-resistant alloys and because adequate information was available concerning its response to various treatments and its subsequent properties. For this study, -100 mesh alloy powder having essentially the same composition as shown in Table I but without the 0.4 pct C was used. Desired carbon variations, achieved by wet mixing the alloy powder with graphite, ranged from 0.1 to 1.0 wt pct. Specimens were hydrostatically pressed and then sintered in vacuum by slowly heating until incipient melting was observed. Sintering this way gave a liquid phase in the interstices of the powder particles and achieved reasonably high densities. The average density of as-sintered specimens was about 95 pct of wrought alloy density. In order to obtain various microstructures, some specimens were heat treated (solution treated for 24 hr at 2150°F, water quenched; aged 16 hr at 1400°F, air cooled); others were hot swaged at 2150 o F (50 pct reduction in area), while still others were hot swaged (50 pct reduction in area) and then given the above heat treatment. During hot swaging it was observed that forgeability was quite good even when specimens contained as much as 1.0 pct C. Specimens were evaluated by stress-rupture test-

Jan 1, 1962

S. Tandanand, Twin Cities Mining Research Center, U. S. Bureau of Mines As an addition to the fine reviews on drilling presented in the foregoing papers, we feel it would be appropriate to mention briefly some work in this area by the U. S. Bureau of Mines. The Bureau has had under way for several years at its Twin Cities Mining Research Center studies directed toward a better understanding of basic drilling mechanisms. The work includes drillability investigations with diamond rotary, and with percussive drilling systems, as well as fundamental studies of rock behavior under drill bit types of loading. Results of our work to date have contributed the following to the state of the art. MECHANISM OF ROCK FAILURE-Rock fragmentation through use of mechanical energy systems, including drilling systems, consists of two major components-machine and rock. The two components are related to each other by their functional behavior and their resultant motion which the components mutually influence. Forces in the systems are transmitted directly to the bit which in turn induces stresses in the rock. When the induced stresses exceed the strength of the rock, the rock fails resulting in penetration. Rock strength is a function of the forces applied to the rock. Stresses under a drill bit, regardless of bit configurations, are essentially compressive, resulting in a condition where the rock has the highest resistance to failure. The ratio of principal stresses induced under drill bits is an indication of the mode of failure in the rocks. Initially, compressive stresses under the drill bit cause rock to undergo inelastic deformation. The impression observed at the first stage of crater formation under the drill bit shows that rock is no longer elastic. Associated subsurface failure at this stage precedes total rupture which subsequently appears at the rock surface. Thus, explanation of crater formation, by applying theory of elasticity alone, is likely to be over-simplified. ENERGY VERSUS VOLUME OF ROCK REMOVED-Experimental results show that energy per unit volume of rock removed is constant for a particular rock and bit type. Rock offers a definite resistance to penetration, depending on the method of the energy application; this may be defined as the drilling strength of a rock. Energy consumption in a drilling system may be divided into two components: 1) Energy required to move the machine parts; and 2) energy consumption in breaking rock. Component 1 is dependent on the type of the machine. Component 2 is dependent on the strength properties of rock and the efficiency of the bit in rock cutting. The sum of the energy required in both components determines the total energy required in breaking rock. Our work in this area has included consideration of both diamond and percussive drilling systems. PENETRATION RATE-For a drilling rate expression to be useful in predicting field performance, it should include the operating characteristics of the drilling machine. Eq. 37, for example, is simply a statement of the fact that the drilling rate is determined by the power output of the machine (2pNT) and the power requirements to break the rock in the hole (AER). In order to predict the performance of a particular drill and the effect, for example, of the rotary speed, on the drilling rate, R, Eq. 37 must be modified to include the

Jan 1, 1967

By Meguru Nagamori

DURING the smelting of the sulfide ores of copper, nickel, and certain other metals, a matte phase is formed. For a better understanding of the physico-chemical properties of mattes, binary systems such as Cu-S, Ni-S, Fe-S, and so forth, should be studied, especially in the vicinity of liquid compounds such as Cu2S, Ni3S2, FeS, and so forth. These compounds are relatively stable in the solid state although they can exhibit appreciable deviations from stoichiometry. With small deviations from stoichiometry there are large changes in the activity of both sulfur and metal this is a characteristic of both solid and liquid sul-fides near the stoichiometric compositions. In the vicinityof molten Cu2S,1 for instance, the activity of sulfur increases about 100 times for a 0.012 increase in the mole fraction of sulfur and, correspondingly, the activity of copper decreases greatly. In this sens the minor deviations from stoichiometry should be carefully taken into account in studies of such metal-sulfur melts. Metal-sulfur melts may be considered to have properties intermediate between metallic solutions and molten salts, i.e., neither entirely metallic nor entirely ionic. The radius of the ion sZ- is 1.84 as compared with 1.04? for a single bond covalent s.2 It is therefore inferred that when the sulfur is more ionic its partial molar volume is larger, although metal-metal, metal-sulfur, and sulfur-sulfur interactions may bring about further complexities. Thus, it is to be expected that the molar volumes, which are derived from density-composition relations, will help in resolving the structure of the melts. In this study the maximum bubble pressure method with bubbling gases of H2, H2S, and Ar, .was employed to measure the density. EXPERIMENTAL Materials. Silver of 99.95 pct purity was melted in an alumina crucible and hydrogen gas was then bubbled through the liquid for the deoxidation. In preparation of Cu2S, electrolytic copper powder (>99.9 pct Cu) was mixed with precipitated sulfur and heated to 1200°C in an alumina crucible under a hydrogen atmosphere. Mixtures having the approximate compositions FeS and Ni3S2 were prepared by heating reagent grade electrolytic iron or reduced nickel powder with precipitated sulfur under a hydrogen atmosphere. The H2S gas was of 98.5 pct purity, the impurities being mostly CO2 and CS2 which hardly contaminate the melts. Apparatus. The crucible system is shown schematically in Fig. 1. To minimize the variation of melt composition during a density measurement, the apparatus was designed in such a way that either Ar, H2, H2S, or mixtures of these gases could be used as the bubbling gas. When the melt had a high partial pressure of sulfur, H2S gas was used. When the partial pressure of sulfur was so low that oxygen might dissolve in the melt, the use of H2 was preferred. In intermediate cases, argon or mixtures of H2-H2S were used. Density of Dibutylphtalate. Dibutylphtalate was used as the manometer fluid. Relative to distilled water, the density of the dibutylphtalate, pdib g per cc, was measured and found to be: pdib (25°C) = 1.046 i 0.002. With the use of the temperature coefficient reported by Kemppinen and Gokcen,3 the following value was

Jan 1, 1970

By L. W. Sink, B. H. Kear, S. M. Copley, M. F. Hornbecker

IT is well known that under appropriate conditions of chill casting, many metals and alloys develop a columnar grain structure having a high degree of preferred orientation.' In cubic materials, this principle has been employed to produce shaped castings with a unidirectional cube oriented columnar grain structure as well as shaped single crystals with a cube axis parallel to the growth ais.' "he purpose of this note is to describe how this same principle can be enlployed to produce wires, rods, or more complex crystal Shapes with a variety of conlpositions in cube or random orientations by a single casting operation. The procedure for directional solidification of in Fig. 2. Examination under crossed polarizers indicated that the striations were a second phase (presumably the B19 phase) which frequently forms initially at primary grain boundaries and then grows toward the interior as thin ribbons. X-ray studies on filings of this alloy reannealed at 1300°C for 15 min indicated a strong pattern of the fcc phase together with a substantial pattern of the B19 phase. The un-annealed filings, however, showed only the pattern of a cold-worked, fcc structure. In these two-phase alloys it appears, therefore, that the B19 phase is transformed into a disordered fcc structure by cold-working. The mechanical properties of these alloys would be of some interest since a fcc structure might be expected to possess greater ductility than the or-thorhombic B19 structure from which it forms. The author wishes to thank Mr. Harvey Yakowitz and Mr. Richard Manuszewski for obtaining the photomicrograph in Fig. 2. This research was partially supported by Research Grant DE-02455 from the National Institute of Dental Research. columnar grained ingots involves pouring a melt into a preheated mold having a water-cooled copper chill at its base? and controlling the temperature gradient in the mold so as to maintain columnar grain growth from end to end of the ingot.4 This is illustrated schematically in Fig. l(a). It should be noted that the longitudinal temperature gradient promoting unidirectional columnar grain growth is maintained by heat loss primarily through the solidifying ingot to the chill, so that both the temperature gradient and the growth rate diminish with increasing distance from the chill. As the distance from the chill increases? the columnar structure coarsens, and the perfection in preferred orientation tends to increase progressively. Using this same procedure of gradient solidification, Fig. l(b) shows a simple scheme for producing a batch of crystals having the same cube orientation rather than an assemblage of columnar grains. Wires and rods are produced by cementing ceramic tubes to the wall of the mold. The steep temperature gradient produced in the central chamber is controlled by the mold heater to produce favorable conditions for crystal growth. Clearly the success of the operation lies in the selection of individual growing columnar grains by some kind of orifice at the base of each channel in the mold. Using simple constrictions, experience with nickel-base alloys has shown that 4 -in.-diam holes at a distance of -1 in. from the chill yield crystals within 7 deg of (100); closer control of orientation is possible by placing the crystal selector farther away from the chill.

Jan 1, 1970

By Frederik W. de Klerk

INTRODUCTION During the past 40 years, there has been in the United States of America and other industrialised Western countries a distinct displacement in the procurement of minerals from domestic sources towards an increased dependence on imported minerals. This was caused in part by one or each of the following: a lack of local supply of particular minerals, domestic source depletion or exhaustion and the relative ease and inexpensiveness of obtaining minerals from comparatively unexploited sources in underdeveloped and some developing countries, notably in the Middle East, South America, Southern Africa, and Australasia. Governments of the industrialised nations, concerned about the reliance of their economies on imported raw minerals examined both the degree of reliability and the extent of mineral resources in the countries and regions supplying the essential minerals. This brought a painful awareness of the vulnerability of industrialised countries as a result of the possible interruption of their mineral imports. In this regard, the United States of America may be quoted as a case in point: a report by a House of Representatives Subcommittee on Mines and Mining (Santini Subcommittee)1 revealed that the United States in 1980 was dependent on foreign sources for more than 50 per cent for 24 of the 32 minerals essential for national survival. The Subcommittee also acknowledged that although the United States may develop its own alternative energy sources, there are in many cases no substitutes for the strategic + minerals imported from foreign sources, from countries which are often unstable at best, hostile at worst. This resulted in various policy responses as well as contingency measures, amongst other, renewed emphasis on the utilisation of domestic mineral resources, an intensified search for substitutes of critical minerals and increased stockpiling programmes. Increased attention was devoted to analyses of the political and economic conditions prevailing in the regions from which strategic and other minerals are exported to the industrialised countries, and the formulation of policy measures vis-a-vis those regions. One region which attracted a considerable degree of attention has been Southern Africa, and in particular South Africa.2 The Santini Subcommittee affirmed that while the bulk of the minerals imported by the United States is supplied by a number of source countries, a host of critical[t] minerals vital to the well-being of that country (and this includes all the Western industrialised countries) originates from Southern Africa, a region, deservedly referred to as the "Persian Gulf of Minerals". Apart from the Republic of South Africa the mineral-rich countries concerned are Angola (diamonds, iron ore, oil), Botswana (diamonds, coal), South West Africa/Namibia (diamonds, uranium), Zambia (copper, cobalt), Zaire (industrial diamonds, copper, cobalt) and Zimbabwe (gold, nickel, chromium, asbestos). The Subcommittee however noted that, with the exception of cobalt, the greatest portion of the reserves of the critical minerals listed, lies within the borders of South Africa. SOUTH AFRICA A RELIABLE PURVEYOR OF MINERALS In view of the continual instability and political risk, along with the nationalisation of mineral industries in third world countries, importing industrialised nations have in recent years become more and more concerned about security of long-term supply, and have in consequence begun to diversify their mineral enquiries to areas where the following conditions prevail: • Long-standing record of consistency of supply; • An economic system of private enterprise; • Political stability; • A well integrated infrastructure, with a maintenance of transportation and development costs at realistic levels, and • Attractive investment opportunities, with security of the return of profits.

Jan 1, 1982

By W. E. Showalter

This paper discusses some of the results of combstion-drive tests which were made in a test cell using a sand bed 10 in. in diameter x 10-ft long. The test method is illustrated and described. The relationship between the API gravity of the in situ oil and the amount of air required for combustion drive is discussed in detail. Other things constant, the air requirement for combustion drive increases as the APl gravity of the in situ oil decreases. If the test results apply to actual reservoirs, the lowest-priced oils may cost the most to recover by this method. Information is shown which indicates that the effect of pressure on the amount of hydrocarbon burned is not large. A method of predicting air requirements from the API gravity of the in situ oil is presented. INTRODUCTION Combustion drive is the term used to identify the process of interstitial or in situ burning as an oil recovery method. Part of the in situ oil is burned to generate the energy needed to produce the remainder of the oil. Combustion drive as an oil recovery mechanism remains an economic uncertainty in spite of all the work that has been done by the industry in both laboratory and field?-6 This paper will show some of the results of tests which were made in a test cell for the purpose of studying the nature of the combustion-drive process. It will present data which indicate that the API gravity of the in situ oil is a significant indicator of the amount of air required to drive a burning front through oil sand. Air requirement varies inversely as the API gravity of the in situ oil. EXPERIMENTAL The tests were performed in a cell which utilized a cylindrical sand section 10 in. in diameter x 10-ft long. The thin-walled metal pipe which held the sand was wound with twenty external electrical resistance heaters which, by means of an automatic controller, maintained adjacent sections of the wall of the pipe at temperatures equal to the temperatures of the contained sand. Each heater covered 6 linear in. of the pipe. By this means lateral heat loss from the sand section was minimized, thereby causing the sand section to simulate more closely a horizontal increment of a combustion-drive reservoir. Fig. 1 shows a schematic diagram of the test assembly. Thermocouples to measure the temperature in the sand were located every 6 in. along the length of the sand section. The pipe containing the sand was enclosed in a cell designed for an operating pressure of 500 psig. The inlet air pressure was controlled at the inlet, and the gas flow rate was controlled and measured at the outlet of the cell. The oil sand used for the tests was prepared by mixing first water and then oil with the noncon-solidated sand using a closed mixer similar to a cement mixer. Table 1 shows a screen analysis of the sand. Ninety percent of the sand was 100 mesh or finer. This sand was a mixture of 80 per cent No. 120 Nevada White Sand and 20 per cent Tennessee Hi-Fusion Moulding Sand No. 3. The Nevada sand was a clean silica sand. The Tennessee sand con-

By J. F. Gormley

Drill pipe, and the connections used to join the pipe together, have had a long history of development and improvement. With the growing use of high tensile-strength pipe, the rotary shouldered connections used to join the lengths of drill pipe have become more critical. Fatigue-type failures of rotary-shouldered connections, or tool joints, have become a cause of much conern to operators and manufacturers. Most fatigue tests on tool joints have been made on the basis of bending moment applied vs number of cycles, rather than on the basis of stress vs cycles. A desire to learn more about the nature of stress in a tool joint due to bending caused a rather intensive experimental study to determine what actually happens to cause fatigue failure. It was known that a tool joint with insufficient make-up torque was subject to early fatigue failure. Therefore, the establishment of the relationship between make-up torque and allowable bending moment would enable a testing engineer to better correlate the results of various tests. The problem of calculating the stress in a tool joint due to an applied bending moment was complicated by: the effect of unknown stress concentrations; the applicability of the well known formula S = Mc/I; and the question of whether to use the combined cross section of the pin and box in calculating the moment of inertia I. TEST EQUIPMENT The test equipment consisted of the following. 1. A lathe modified to act as a fatigue machine for full-size tool joint or drill-collar specimens. 2. An adapter plate bolted to the headstock of the lathe and into which a short adapter was screwed. 3. Tool-joint pins on sections of drill collars bucked into the adapter. 4. A hydraulic load cylinder which rested on the bed of the lathe and applied bending moments to the joints. 5. Paper-backed strain gauges (FAP-12-12, 1/8-in. gauge length and A-5-1, 1/2-in. gauge length). 6. Eastman 910 contact cement for bonding gauges to joint and Epoxylite 222 waterproofing. 7. A 12-position switchbox, SR-4 bridge balancing unit connected to an SR-4 strain indicator to read the gauge output. 8. A four-channel Sanborn recorder, Model 150, was used as the readout instrument for part of the testing. The fatigue machine was well suited as a loading device. Because of the long lever arm available, the change of moment with distance was small; consequently, the vertical shear stress was small, less than 300 psi. Some difficulty was experienced with the hydraulic load cylinder at high bending moments due to leakage. The strain gauge equipment was simple and quite reliable. A 41/2-in. extra hole (E. H.) joint with "as machined" threads and shoulders and a 2 7/8-in. internal flush (I. F.) joint with copper-plated threads and shoulders were selected for testing because of their common usage in the field. In addition, the 4 1/2-in. E. H. was suitable because it had an inside diameter large enough to permit the insertion and accurate positioning of strain gauges without using special jigs, while the 2 7/8-in. I. F. joint required much lower bending loads to cause separation of the shoulders of the joint when made up to the recommended torque. The joints were made to standard API specifications, except for very small bevels on the box and pin shoulders and cylindrical stress-relief grooves added to the pins. Larger bevels would have

By E. R. Buchanan, R. E. Reed-Hill, F. W. Caldwell

Poly crystalline zirconium tensile specimens containing a sizable fraction of grains unfavorably oriented for slip were deformed at 77°K to strains as high as 9 pct. The contributions of the various twinning modes to the observed strain were deter~nlned quantitatively using measured twin volume fractions and Schrnid orientation factors. The results indicate that at very small strains twinning can account for almost all of the observed strain. At higher strains slip becomes the dominant means of deformation. The most significant form of twinning was {1121}. PLASTIC deformation in polycrystalline aggregates is poorly understood even today. This is primarily because slip cannot be easily observed or studied in the interior of a metal. On the other hand, twinning can be discerned inside a specimen after deformation. Twinning studies can therefore often yield quantitative information about polycrystalline deformation not otherwise obtainable. This investigation uses the above aspect of twinning in determining the contributions of twinning to the strain during tensile deformation at a cryogenic temperature (77°K). Polycrystalline zirconium was used as the object of study. There are several reasons why this information should be of value. In contrast to many metals, both the strength and ductility of zirconium increase with a corresponding decrease in temperature. This may be related to the fact that twinning becomes more profuse and complex with decreasing temperature. Secondly, it has been discovered that certain room-temperature properties of zirconium can be greatly improved by about 1 pct prestrain at 77°K. These include an increase in tensile elongation in critical directions by over 50 pctl and a large mechanical hysteresis effect.2 Both effects are believed to be related to twins formed by prestraining: plastic twin-boundary movements accounting for the increased ductility, and anelastic twin-boundary movements for the internal friction. Thirdly, hexagonal metals which twin primarily on only the (1012) mode (such as beryllium, zinc, and magnesium), unlike zirconium, are very limited in ductility at cryogenic temperatures. Therefore, the significance of the various types of twinning in zirconium at cryogenic temperatures is of more than academic interest. Three twinning mechanism-s have beenobserved in quantity at WOK, {1012), {1121), and {1122). Each must be considered as a unique deformation mechanism because of their markedly different twinning shears and basically different angles between twinning and basal planes, as shown in Table I. It is therefore necessary to separately consider the role of each of these twins during plastic deformation at 77°K. Fortunately, each twin possesses unique morphological features which makes its identification possible on a properly prepared metallographic section. The use of the polarized microscope3 can also greatly aid in identification of the various forms of twins. EXPERIMENTAL PROCEDURE Arc-melted sponge zirconium plate, previously describedY4 was machined into 1/8 by 1 1/4 in. threaded-end, round tensile test specimens. These specimens were deformed to the desired strains in a liquid-nitrogen cryostat mounted in an Instron testing machine. The tensile axis was the transverse plate direction. This specimen direction resulted in large numbers of grains in orientations unfavorable for deformation by prism or basal slip. Deformed specimens were sectioned, mounted in epoxy cold-mount resin, ground on metallographic papers, electropolished, and anodized to provide surfaces sensitive to polarized light. Standard quantitative metallographic procedures5

Jan 1, 1965

By L. T. Sanchez, W. T. Rogers

This paper presents the results of a production experiment evaluating the effect of the use of oxygen in the bessemer converter with respect to its relation to blowing time, amount of steel scrap or cold pig iron which can be melted and teeming temperatures of the molten steel. THE simplest and probably the most efficient and economical means for producing steel from molten pig iron is the bessemer converter. This process is essentially a mechanism for bringing oxygen into contact with the molten metal thereby oxidizing the elements silicon, manganese, and carbon. The source of oxygen is that contained in the air and is limited by that amount. In theory, there is an intermediate chemical reaction due to the mass action of iron and air which is represented by the equation: 2Fe + O2- = 2Fe0 The FeO then becomes the oxidizing agent and reacts with the silicon, manganese, and carbon to produce silicon dioxide (SiO,), manganese oxide (MnO), and carbon monoxide (CO), in the order mentioned. The balancing of the various chemical reactions shows that in the ordinary bessemer process more than sufficient oxygen is delivered to satisfy these reactions, which indicates that a further increase in the oxygen content of the air would serve no useful purpose. In spite of this reasoning, however, it will be demonstrated that when melting stelel scrap, commercially pure industrial oxygen introduced into the normal air blast will permit the melting of a larger proportion of steel scrap than is possible with air alone and also will increase the production rate by a substantial reduction in blowing time. It will also be shown that when melting cold pig iron, it is possible to melt up to twice as much pig iron with the use of oxygen as can be melted without added oxygen with a nominal increase in blowing time. At the plant where this study was made, oxygen is purchased from a nearby supplier and piped directly into the plant. It is supplied to the bessemer converter by being introduced directly into the blast line through an automatic control system which may be adjusted by the bessemer operator or "blower." Normally, iron is supplied by five blast furnaces to an open hearth shop containing twelve 165-ton capacity stationary furnaces and a bessemer converting shop having three 30-ton converters. The molten metal from the blast furnaces is allocated to one shop or the other in a manner to best provide for maximum production. At various times, it is desirable to direct the maximum amount of available blast furnace production to the open hearth shop. This has the effect of reducing the production of the bessemer shop. When this condition exists, it is desirable to supplement the reduced iron charge at the bessemer with steel scrap or cold pig iron to the maximum possible extent in order to maintain high total ingot production. This condition becomes aggravated when one of the blast furnaces is taken out of production for any repair job and particularly for extended periods such as when relining is necessary. The purpose of the investigation and experimentation with oxygen was to determine whether or not it would contribute to increased flexibility in the use of the raw materials, molten iron, steel scrap, and cold pig iron. In addition to the effect of oxygen on flexibility of operation, it was also necessary to determine the relation of its use to production rate or "blowing time," final finishing temperature of the steel or "teeming temperature," and steel quality.

Jan 1, 1953

By D. J. Pye

Certain high molecular weight synthetic polymers in very dilute solutions decrease water mobility in porous media 5 to 20 times more than would be expected from the solution viscosity. This indicates that increasing water viscosity to reduce adverse mobility ratios in water floods is economically feasible in many situations. Laboratory results are consistent with accepted theory with respect to improved areal sweep, permeability distribution and displacement efficiencies as would be predicted from the solution mobility measured in the core. The laboratory results have been confirmed by field pilot floods which recovered 80 per cent and 100 per cent more oil than comparative water-flood pilots in fields containing 16 cp and 130 cp oil respectively. Economic analyses based on pilot data indicate the method to be profitable. INTRODUCTION The desirability of improving the water-oil mobility ratio in waterflood-ing operations has long been appreciated, especially for high viscosity crudes. Methods of reducing the oil viscosity by application of heat, gas resaturation and miscible drives have met with some success. Conversely, it has been generally conceded that increasing the water viscosity would also be effective, but the amount of thickening agent required to effect an appreciable increase in viscosity is discouragingly high. The use of such materials as glycerin, sugar, or glycols is completely out of the question economically, so that considerable attention has been given to the use of the much more efficient synthetic water-soluble polymers as viscosity improvers; but even these materials are not economically attractive. If increased water viscosity could be economically realized, marked improvements in areal sweep efficiency, as discussed by Caudle and Wittel and others, would be realized in field operations. Similarly, the improved mobility ratio would bring about increased recoveries by correcting permeability distribution problems as discussed by Stiles2 and Dykstra and Parsons,' and improvements in displacement efficiency would be obtained as recognized by Buckley and Leverett4 and Welge.5 LABORATORY STUDIES THE RESISTANCE EFFECT In the course of research to develop a more efficient water-soluble polymer viscosity builder, we discovered that a very few types of water-soluble polymers exhibited a most unusual and very interesting property not previously observed. Among the earliest materials showing such activity were polymers containing acry-lamide. Fig. 1 shows a typical plot of the solution viscosity vs concentration for one of these polymer solutions as determined in an Ostwald viscometer. As usual, such a plot is an approximate straight line on semi-logarithmic paper. Now, an inspection of Darcy's law, shows that a viscosity value can also be determined in a formation sample if the pressure, flow rate and permeability are known. If the core permeability is previously determined with wa- ter or brine, then the viscosity of another aqueous fluid may be determined from the equation. The viscosities so determined for glycerin (glycol or polyvinyl alcohol solutions, for example) agree with the viscometer determinations. However, in the case of the water-soluble polymer solutions in which we are interested, the viscosities measured in the formation sample depart very markedly from the viscometer values. A typical relationship for a polymer solution is also shown in Fig. 1. This unusual departure from the expected response is referred to as a resistance property of the polymer and is quantified as "resistance factor". Since the resistance factor R can best be determined in a core under study, the term has been defined on the basis of the ratio of the brine mobility to the polymer solution mobility at residual oil saturation so that:

Jan 1, 1965

By K. P. Gupta

The Cb-rich boundary of the (Cb,Al) a phase field at 1250OC is near 41 pct Al. The Al atoms tend to occupy the C. N. 12 sites in this structure. A homologous (Ta,Al)a phase was identified. No a phase was found in the Mo-Al system at 1200°C. The first binary o-phase with a non-transition element as a major component was found1,2 in the Cb-Al system. The composition of this phase was reported2 as approximately 34 at. pct Al, based on X-ray diffraction patterns taken of alloys in the as-cast condition. Further study of the composition of this o-phase, therefore, appeared desirable. Since Si was previously found3,4 to occupy preferentially the lowest C.N. sites (I and IV) in the a phase, and preliminary resultsa indicated that Al may show similar ordering, one of the aims of the present work was to obtain more complete information on ordering in (Cb, Al) o. Among the binary systems of Al with Ti-, V-, and Cr-group elements, the Ti-Al, Zr-Al, V-Al, and Cr-Al systems are known to comprise no o-phase. In the present work a Ta-Al and a Mo-Al alloy were studied in an effort to find additional o phases with Al as a major component. EXPERIMENTAL PROCEDURES AND RESULTS An alloy with the intended composition of Ta + 40 at. pct Al was prepared by arc melting in He atmosphere. In both the as-cast condition and after annealing for 10 days at 1250°C and quenching in cold water this alloy consisted of two phases, as seen in the X-ray pattern and the microstructure. One of these phases was identified as the o-phase; the observed and calculated d-spacings are shown in Table I. The lattice parameters of this o-phase are a = 9.972 A, c = 5.214 A and c/a =0.5228, as determined from a Debye-Scherrer pattern using Fe K a radiation. The d-spacings for the unknown second phase are given in Table 11. The alloy contained approximately 25 pct of the unknown phase after annealing, as estimated on the basis of microscopic examination. A Mo-Al alloy was prepared by arc melting, with an intended Al-content of 35 at. pct. The alloy button was annealed for 7 days at 1200°C and water quenched. The X-ray pattern of the alloy in both the as-cast and the annealed state showed the reflections of Mo3Al (Cr30 type structure) and an unknown phase different from a. Six Cb-Al alloys were prepared in the composition range 30 to 50 pct Al. All alloys were arc melted in He atmosphere. After annealing for three days at 1175oC in evacuated silica capsules and quenching in water, the alloys with 30 to 40 pct Al contained three phases, namely the o-phase, Cb3Al (Cr3O-type structure) and an unknown phase. The alloys with more than 40 pct Al contained two phases, namely, what appeared metallographically, the o-phase and the same unknown phase. In the alloys containing less that 38 pct Al, the unknown phase was not present after annealing for 4 days at 1250°c,

Jan 1, 1962

By M. Bevis, A. G. Crocker

A. G. Crocker and M. Bevis (Battersea College of Technology)— Richman and Richman and conard17 have recently published the very interesting results of their studies on the plastic deformation of Fe-Ni-C martensites and of an Fe-5 wt pct Be alloy. Other work on the range of Fe-Be alloys containing between 16 and 27 at. pct Be has also been reported by Richman." One of the fascinating features of this work is that the materials investigated deform primarily by means of twinning and that several twinning modes are operative in addition to the well-established (112 ) mode of bcc materials The composition planes of these unusual twins were determined using a single-surface analysis and reported as (0131, "{045)" or "{089)", and "{2,3,10)" or "(127)". (It is now thought1' that the "{0,1,13)" habit reported for the Fe-Ni-C alloys'' is more likely to be near "{127)". Although the full twinning elements of these unusual twinning modes could not be determined directly from the experimental results, the authors attempt to deduce the elements by considering twin intersectionslQ and by using the projection technique of Jaswon and Dove.' The basic assumption of their approach is that the (013) habits are truly rational planes but that the other habits are irrational and merely approximated to by the indices quoted above. They thus deduce that the twinning mode describing the (013) habit is either Type I or compound and that the modes describing the other habits are Type 11. This deduction must however be false as all twinning modes in cubic structures obeying the classical orientation relationships of deformation twinning must be compound. Type-I and Type-I1 modes can occur in all of the other crystal systems (indeed in the triclinic system compound modes are not allowed), but in the cubic system only compound modes can arise. This fact, which does not appear to be well-known, fol- lows immediately from the relations between the four twinning elements defining a twinning mode. Thus for a Type-I mode which is uniquely defined by the rational twinning elements K1 = hi and tj2 = vi (i = 1, 2, 3), the Kz qd 771 elements are given by: where the summation convention of the tensor calculus has been adopted. For cubic materials the metrics cij Iand cij reduce to the unit tensors 6ij and 6ij so that we have It is clear from these equations that both ki and ui must take on integral values; i.e., Kz, and ql must be rational, if hi and vi are integers, i.e., if K1 and 77~ are rational. Thus for cubic materials a Type-I mode must degenerate into a compound mode and the same result also arises for Type-II modes. The apparently irrational habits observed by the authors can thus not be explained in terms of Type-I1 twinning if the classical laws of deformation twinning are to be obeyed. An alternative explanation of the irrational habits is that they may arise from twinning modes with four irrational twinning elements. These modes, which do not obey the classical orientation relationships associated with deformation twinning, have been shown to arise as the result of double twinning processes in magnesium,22 but may also describe a simple shear on an irrational plane in an irrational direction. Indeed double twinning processes in bcc

Jan 1, 1963

By P. Hazebroek, M. Prats, E. E. Allen

Well locations in water floods frequently do not coincide with the ideal well locations associated with the flood pattern which is used. A question of importance is. "How much oil is lost because of this?". This is the problem treated in this paper. It is assumed that each pattern in a five-spot flood has the proper number of production and injection wells but that the wells may be off-pattern. The distance by which a well is off-patern is given in terms of the pattern spacing, defined as the length of one side of a five-spot. The influence of off-pattern wells on the performance of a five-spot flood has been investigated mathematically for unit mobility ratio. All production and injection rates are assumed to remain equal. The loss in recovery resulting from pattern irregularitie: has been calculated at a cumulative water injection corresponding to a 98 per cent water-cut in a regular five-spot. Displacing a single production or injection well by about one-third the pattern distance reduces the total field recovery by 4.4/N per cent of the recovery obtainable with regular spacing, where N is the number of five-spots in the flood (assumed large). When off-pattern wells are separated by at least one regular five-spot and are displaced by less than one-third the pattern distance, the total loss in recovery in a large-scale flood is, at most, 1 per cent of the recovery obtainable from that field with a regular five-spot pattern. However, much higher losses may occur if the off-pattern wells are not separated by a regular five-spot, if the rates are imbalanced, or if the mobility of the injected water is higher than that of the crude. Although the loss in oil recovery due to pattern irregularity often may be low, the off-pattern case will always require injection of more water to reach the same final water-cut. The volume of water injected per unit volume of oil displaced has been calculated for each five-spot adjacent to a five-spot containing a single 08-pattern well. The cumulative water injection at the time of 98 per cent water-cut may be increased by as much as 37 per cent when a production well is displaced one-third the pattern distance toward an injection well. The existence of differences between actual subsurface locations and the idealized well patterns used in secondary recovery operations is well recognized. However, the methods commonly used for predicting waterflood behavior ignore such differences, even when subsurface locations are known. The reason for this is simple. The idealized well patterns are highly symmetrical and, as such, the production history of the pattern (and even of the field) can be determined from that of a few production wells—usually from only one. Consideration of the many irregular well spacings found in a specific field would entail a large amount of work, and then the results would be applicable only to that specific field. In an attempt to provide an answer to the more general problem, this paper considers a water flood having an idealized well pattern, except for one well which is not at its regular location. When several off-pattern wells are present, their contributions are assumed to be additive unless the off-pattern wells are close to each other. Only four positions of the off-center well (at x, y. of Fig. 1) are considered in this paper. These are shown

By F. W. Jessen, R. W. Durie

The dissolution of salt in the development of salt cavities is controlled by free-convection boundary layer flow along the salt surfaces. It is the purpose of this paper to expand upon results published previously, and primarily to illustrate the influence of surface features on the boundary layer and hence the dissolution process in the development of underground salt cavities. General dissolution equations are presented that apply to a surface inclined to the vertical. A salt surface inclined so as to be overlain by water is dissolved by boundary layer flow, while an inverted surface is dissolved by a combination of boundary layer and cellular flow. Detailed investigation of the boundary layer was limited to laminar free-convection flow. Truly laminar free-convection boundary layer flow may be encountered only in laboratory solution cavities, but the boundary layer will control dissolution even where cavity dimensions are such that free-convection turbulence will occur. The factors found to significantly affect the laminar boundary layer will affect dissolution with a turbulent boundary layer. INTRODUCTION Underground salt cavities have been used extensively for the storage of LPG products, and the use of these cavities is expanding each year. The increasing requirements for such cavities suggest the need for a better understanding of the mechanisms involved, for the purpose of the control of cavity geometry and the most effective use of the injected water. A previous paper outlined a quantitative approach to the determination of salt dissolution rates in the laminar free-convection dissolution of salt. It is the purpose of this paper to expand upon the results previously published and to point out significant factors in the application of these results to the dissolution of field cavities. The dissolution of salt in laboratory scale cavities is controlled by the free-convection boundary layer, which is the region of flow along the salt surface caused by the difference in buoyancy forces resulting from the salt dissolution. The transfer of salt into solution at the salt face causes an increase in salinity, and hence a greater density, in the water at the salt surface. The region of higher density tends to fall through the lighter, less saline water and a flowing boundary layer is the result. The water remains stationary at the contact with the salt surface,l so the resulting velocity profile is such that the influence of gravity is opposed by a viscous shear force at the salt surface. For salt to be transferred into solution from a solid surface to a boundary layer of fluid, it must diffuse through a stationary film of water at the salt contact. According to Fick's first law, the rate of mass transfer by diffusion is the product of the diffusion coefficient and the concentration gradient, with diffusion transfer in the direction of decreasing salinity. Once the salt has passed into the boundary layer, it will be transferred by a combination of molecular diffusion and convection. The convection transfer will in turn affect the rate of dissolution by influencing the gradient of the salt concentration at the surface. In the previous investigation by the authors to determine the basic mechanism of the salt dissolution process, an equation was developed that agreed very well with the experimental rates of free-convection salt dissolution from a smooth vertical surface. The technique was taken from the work of Eckert2 in the analogous heat transfer system and will be referred to as the integral method. The rate of salt transfer from a vertical section of height H was found to be:

Jan 1, 1965

By H. D. Yoo, M. R. Tek, D. L. Katz

The behavior of gas-storage reservoirs subject to water drive is investigated through analog simulation on an electronic differential analyzer. The simulation technique developed on an LM-10 computer permits the prediction of reservoir volume or pressure resulting from the movement of water in the surrounding aquifer. The method developed on the analog computer consists of setting up an appropriate transfer-function circuit and feeding the arbitrary time-varying boundary conditions as an input signal. The input may be specified as gas reservoir pressure, pore volume or the water flux. Several cases studied include an isolated gas reservoir on a limited aquifer, interference among three reservoirs adjacent to a common aquifer and the growth of gas-storage volume on an aquifer. It is concluded that the method developed on an electronic differential analyzer provides an excellent technique to simulate and investigate the behavior of gas reservoirs subject to water drive. The agreement between the reservoir performance as predicted from the simulation technique and as measured from actual field data is found to be better than the range usually encountered in predicting water-drive behavior. INTRODUCTION It is generally known that some gas-storage reservoirs are located on top of blanket sands of large extent, saturated with brine called aquifers. Because the volume of the body of water associated with aquifers is usually very large and water is compressible, the cyclic pressure variations encountered in normal storage service inevitably cause unsteady, compressible flow conditions in the adjacent aquifers. The solution to radial diffusivity equation for a limited aquifer for constant terminal conditions has been know111 since the early 1930's. The solution for the constant terminal conditions for an infinite aquifer was puhlished in 1949.2 The proldem of handling the unsteady flow of water through a porous medium with time-varying boundary conditions has been studied by many investigators through the use of graphical,3 digital4 computing techniques or on an R-C type of electrical analog called a reservoir analyzer.5 The digital computer performs a series of arithmetic calculations to superimpose the solutions of constant terminal conditions, while the R-C type of analyzer permits the simulation of the continuous medium by means of "lumped-cell" electrical resistance and capacitance elements based on finite-difference approximations. The graphical method is based on the convolution integral solution of the partial differential equalion. The purpose of the present study was to simulate the behavior of an aquifer on a general-purpose, electronic differential analyzer to obtain the influx of water or the change in pressure at the reservoir continuous1r and quickly. The use of the electronic differential analyzer was based on well known solutions for the constant terminal conditions already available in the literature.2,6-8 THEORETICAL DEVELOPMENT MODEL. The gas reservoir (R) is described in relation to an aquifer in Fig. 1. The permeability and porosity of an aquifer sand, the viscosity of the brine, and the compressibility of the brine and the formation are assumed to be constant. It further is assumed that the blanket sand is not thick enough to give appreciable vertical pressure distribution. The model in consideration is circular with radial geometry. The aquifer radius can be of limited or infinite extent. The gas-water interface is assumed to be fully segregated, and the displacement of the gas by

By N. S. Stoloff, R. G. Davies

The annealing behavior of a heavily cold-worked FeCo-V alloy has been studied at temperatures both above and below Tc, the critical temperature for ordering. It was found that re crystallization and grain growth can take place below Tc in the ordered state and that, in agreement with results for fcc alloys, order retards grain growth. The activation energy for pain growth is -55 kcal per mole in the disordered condition; it was not possible to assign a value to the activation energy for grain growth below Tc because the degree of order is changing rapidly over the temperature interval investigated. Precipitation of an fcc y phase was observed upon annealing the cold-worked alloy at temperatures of 675°C and below. SINCE it is known that diffusion constants are greatly affected by long-range order,' it is not unreasonable to anticipate that order will influence the diffusion-controlled phenomena of recrystallization and grain growth. In the course of a study of strain aging of ordering alloys around the composition Ni3Fe Vidoz, Lazarevic, and cahn2 noted that a disordered cold-worked alloy annealed at a temperature T1 just below Tc will order but not recrystallize; Tc is the critical temperature for order. The same alloy annealed at a temperature T2 just above Tc, or an alloy of different (nonordering) composition annealed at either T1 or T2, will recrystallize after equivalent amounts of cold work. Initially Selisskii and coworkers,3,4 who examined a series of cold-worked and annealed bcc Fe- Co alloys, reported that alloys forming the B2 FeCo super lattice could not be recrystallized below Tc. However in a later paper5 it was shown that FeCo ordered alloys could be recrystallized below Tc. Thus the first question to be answered by the present investigation was whether or not a heavily cold-worked Fe-Co-V alloy could be recrystallized below Tc. When it was found that the alloy would indeed re-crystallize below Tc, the influence of long-range order upon subsequent grain growth was studied. EXPERIMENTAL PROCEDURE The grain-growth measurements were made upon FeCo-2 pct V wire of 0.031 in. diam; the alloy was the same as that used in an earlier study of the mechanical properties of FeCO-v.6 The wire was produced by cold swageing and drawing a rod 0.125 in. diam that had been quenched into iced brine from 850°C to disorder the alloy; the initial grain size was 0.05 mm. Rates of grain browth were measured after isother-mally annealing the wire in the temperature range 675° to 825°C; T, is -720°C. Grain diameters were measured on the transverse cross section of the wires. The line-intercept method of determining the grain diameter was used; conventional light photomicrographs were taken when the grain size was greater than 10 . Carbon replicas were produced and examined in the electron microscope if the grain size was less than Dislocation configurations after recrystallization anneals were studied by electron transmission microscopy of foils produced by cold rolling a 0.125-in.-thick strip of the quenched material to 0.012 in. The foils, after heat treatments similar to the wire samples, were thinned by electropolishing in a chromic-

Jan 1, 1967