Search Documents

By C. A. O’Connell

Surface exploration drilling at Mufulira started in 1928, with drillhole No. 1 intersecting ore near the eastern fringe of the deposit. Subsequent drilling along strike intersected ore until the western fringe was fixed some 7500 ft west of drillhole No. 1. Further intensive drilling down dip disclosed a large triangle of ore with its apex 3350 ft deep. It contained approximately 11 6 million tons of ore at 4.41% total copper in three superimposed orebodies. The east and west limits of the orebodies were presumed to have a pronounced rake to the north with a strike length of 8000 ft. Subsequent exploration underground seemed to bear out this presumption to the 1150 ft level. In 1947, attention was directed to the area west of the known orebodies in conjunction with plans for increased production and increased depth of the mine workings.

Jan 12, 1962

[MUNOZ, P.: New Oil Fields in Trinidad, M27, 311 MUNRO, C. H., biography, M84, 354 MUNROE, H. S., biography, M83, 277 Muntz metal (See also 60:40 Brass): corrosion. See Corrosion of Metals. fatigue and corrosion-fatigue, 78, 597 heat-treated: cleavage structure, 83, 498, 505 grain growth, abnormal, 83, 503, 505 machinability, 99, 323 stress-corrosion cracking, 89, 249. 266 MURAKAMI, T.: iron-silicon diagram, 80, 253, 255 MURDOCH, J. A. W., biography, M28, 157 MURDOCK, W. J.: Discussion on Adaptability of Various Coals as Generator Fuel in Manufacture of Water Gas 78, 561, 565, 566, 568 MURPHY. D. W. AND CHIPMAN, J.: Solubility of Nitrogen in Liquid Iron (T.P. 591) 116, 179. Abs., YB35, 33; M35, 114 MURPHY, L. J.: Relative Advantages and Costa of Electric Power in Lease Operations 77, 219]

Jan 1, 1936

By Masatsugu Taniguchi

Onoda Cement Co., Tsukumi is one of the largest limestone producer in Japan. From mines called Suisho and Ganji, it supplies limestone for cement plant, lime manufacture, iron and steel industries and for many purposes. The Suisho mine is highly productive, though it's operations are governed by the utmost regard for its impact on the surrounding environment. However, it became difficult to meet additional expanding demand and to assure continuity of high rate production. Present production rate 1s 6 million tons per year. The Ganji mine, some 2 km south west of the Suisho mine, began first-stage development In mid- 1973 and in january I976 started production. The capacity of major facilities is 9 million tons per year. Full production is predicted in 1979. Mining and milling methods of these mines under strict environmental considerations are discussed.

Jan 1, 1976

By Robin O. Williams

Robin 0. Williams (Oak Ridge National Laboratory)— The authors have questioned the degree to which the coherency strains between the iron-rich and chromium-rich phases are isotropic as proposed in Ref. 5 on the basis of the difference between the elastic properties of the two phases. The relative magnitude of the stresses is determined by the moduli as shown by Eqs. [2], [3], and [4] of Ref. 34. However, the moduli of the two phases have no direct bearing on the uniformity of either the stress or strain within either phase. The idea that the strains are isotropic within each phase (but normally of different magnitude and always of different sign) is based entirely upon the experimental observation that X-ray line broadening has not been detected even when the particles become rather large. It has not proven possible to grow the particles sufficiently large that they lose coherency. Based upon this lack of line broadening one can estimate an upper limit for the nonuniformity of the strains within each phase as follows. It is considered possible to detect line broadening if it is as great as 10 pct of the separation of the K, doublet for the (211) line using chromium radiation. The doublet separation would correspond to a total strain of 0.0017 such that the total variation of lattice parameter relative to the average lattice is now k0.05x0.0017 or something less than ± * For the present case the strain in each phase is roughly 0.002 such that the variation of strain within a phase will not exceed 5 pct. It is stated that the expression derived for strengthening for the hydrostatic straining as observed in this system would substantially overestimate the magnitude due to dislocation flexure. This is contrary to the conclusion reached in the original paper34 for the present range of particle sizes. What is the lowest temperature at which a has been observed to form in this alloy? M. J. Marcinkowski, R. M. Fisher, and A. Szirmae (nutlzors' reply)— -Williams' arguments based on X-ray findings for a chromium-rich precipitate and an iron-rich matrix strained to a common lattice parameter are certainly convincing. This being the case, there are no shear components of strain associated with the precipitate-matrix aggregate to interact with the shear components of the dislocation stress fields, contrary to the opinion expressed by the present authors. On the other hand, the present authors, in spite of this error, did not expect the shear interactions to be significant. The chief objection to Williams' model in the present case is that the various segments of the dislocation line are assumed to pass from one potential valley to the next independently of neighboring segments. This is only true for a highly flexible dislocation line, i.e., one whose radius of curvature is something less than the center to center distance between precipitate particles which amounts to about 90A in the present alloy. In order to maintain this curvature, an externally applied shear stress of at least 230,000 lb per sq in. would be required or about four times the observed stress. It is therefore concluded that the dislocation lines move rather rigidly through the lattice. This being the case, the forces on the dislocation resulting from the hydrostatic interaction between the stress fields of the edge-dislocation components and the precipitate particles should average out to zero; that is particles above the below the slip plane produce forces on the dislocation of opposite sign and therefore will cancel when averaged over the entire length of the dislocation. On the other hand, since the dislocation is not perfectly rigid, Williams' model may lead to some strengthening, but far less than that predicted. A second and equally serious objective to using Williams' strengthening model for the present alloys is that profuse wavy slip due to the motion of screw dislocations played a predominant role not only in the unaged alloys but in the fully aged ones as well. Since the screw dislocation has associated with it only shear components of stress the hydrostatic strengthening model no longer applies. In view of these arguments the present authors must reject Williams' model of strengthening as being pertinent to the present alloy system. The present authors have made no detailed study of the lowest temperature at which a forms in the quenched ferritic alloys. None was ever observed n the alloys aged at 500°C so that forma-tion must occur at temperatures higher than this and was therefore not a factor in the present study.

Jan 1, 1965

By Irving B. Cadoff, Alvin A. Machonis

The resistivity, Hall coefficient, Seebeck coefficient, and thermal conductivity were measured as a function of temperature for cation-rich alloy single crystals covering the composition range across the PbTe-SnTe system. Alloying of PbTe with up to 20 pct SnTe was found to have little effect on the energy gap. Above 20 pct SnTe the alloys were "p" type but below this range the sign could be varied by heat treatment. The lattice thermal resistivity of the compounds SnTe and PbTe is raised by alloying one with the other. Z values in the order the interesting values obtained. THE PbTe-SnTe system has several interesting features. For one, PbTe is a useful thermoelectric material and the possibility of improving its figure of merit by alloying with SnTe, an isomorphous compound, has been suggested since these pseudo-binary solid solutions generally have a more favorable ratio of electrical conductivity to thermal conductivity than either of the components.' Other interesting features relate to the conductivity mechanism, band structure, and stoichiometry of the compounds and their alloys. PbTe is a semiconductor with an energy gap of about 0.29 ev2 at room temperature whose conductivity sign and magnitude can be varied from "n" to "p" by controlling the proportion of lead and tellurium with respect to the stoichiometric ratio.3 Excess lead results in "n"-type conduction. SnTe is found to exist only as a "p"-type material of relatively high conductivity. This behavior is attributed to stoichiometric deviation by Brebrick4 but Sagar and Miller proposed that the behavior of SnTe must be due in part to the presence of an overlapped band. An investigation of alloys of this system, therefore, might give additional information which would permit one to evaluate which of the two proposals is the more appropriate one. Abrikosov et al.' studied the room-temperature electrical properties of these alloys and reported data for Seebeck coefficient and resistivity on poly-crystalline alloys. The present work is a more exhaustive survey of the PbTe-SnTe system. Re- sistivity, Hall coefficient, Seebeck coefficient, and thermal conductivity were measured over a wide temperature range for single crystals at 10-pct intervals of lead/tin ratio across the pseudobinary system. The relative concentration of tellurium was controlled so as to obtain metal-ion excesses in all cases. SAMPLE PREPARATION The crystals were prepared by melting elemental lead, tin, and tellurium in weighed proportions in evacuated Vycor capsules. The lead and tellurium were high-purity grades obtained from American Smelting & Refining Co. The tin was supplied by Comico. The proper calculated proportions of lead, tin, and tellurium were weighed and charged into prepared Vycor capsules prior to evacuation. The capsules were prepared from 15-mm Vycor tubing. A sharp point was worked on one end of the tube. A pyrolytic graphite coating was deposited on the Vycor walls by heating the tip to 800°C in an atmosphere of acetone-saturated argon. An additional coating of graphite was deposited on the pyrolytic coating from an Aquadag suspension. Above the coated tip the tube was reduced in diameter to form a constrictive neck. To avoid scratching the graphite coatings the charge was placed in the tube above the constriction. After a low-temperature bake, the evacuated capsule was sealed. On subsequent heating the charge melted down into the lower portion of the capsule. The crystals were grown by lowering the capsule through a Bridgman-Stockbarger furnace. The lowering rate was 1 in. per 8 hr. The upper portion of the furnace was set for 950°C and the lower portion for 800°C. In general the yield of single crystals was about 25 pct. The mixed compositions were, as expected, the most difficult to grow. The finished crystals were sectioned into 5/8-in. slices. The tip, end, and middle slices from each crystal were analyzed by X-ray fluorescence to determine the lead-to-tin ratio. The resulting values were used to plot a composition vs distance plot for each crystal. Slices were selected from each crystal, with the aid of the composition plots, to cover the complete range of compositions at 10-pct intervals. In general, the slices selected were taken from the seed end of the crystal where the longitudinal segregation (as determined from the X-ray fluorescence analysis) was a minimum. Laue single-crystal analysis and metallographic analysis was used to verify if a slice was single or polycrystal. Any grain boundaries were clearly visible in the as-cut and polished condition. In ad-

Jan 1, 1964

By E. Kendall Cork

The traditional financial objective for a single mine company has been to operate as frugally as possible and to pay out most of the earnings as dividends. If the business is cyclical (as it is for most metals) the dividends might fluctuate quite widely. When the mine is exhausted the company disappears. This is still quite a viable strategy for a single mine company. It is not however a viable strategy for the world as a whole. The mining industry is built by mine development companies who can mobilize the people and capital to bring new mines into production. Their skills must include marketing, engineering, finance and other politics. It is very rare for a property to be brought in without the support of a major company that can provide all these services. The exceptions will usually have some other form of big brother support, for example the U.S. government uranium contracts at guaranteed generous prices. The mine development company will seek as a minimum to perpetuate itself by developing new mines in order to replace those which are running out. The more common and more ambitious objective is to grow -- that is to add to its ore reserves and current production by developing more new mines. The financial objectives for that company are very different. Obviously if all the earnings were paid out in dividends there would be nothing left to work with. The first financial policy then is to spend an appropriate amount on exploration for new properties. The next is to retain enough of the earnings to provide the capital for new projects at least sufficient for the equity. There is no magic formula as to what proportion of earnings should properly be distributed as dividends by a growth-oriented mine development company. As a rough rule of thumb distributing half or more will probably leave too little to work on and 30% or so is probably a good balance. However the circumstances differ widely from company to company. It may be useful to set an objective for the rate of growth of a company's earnings. Some have picked rates such as 15% per annum compounded. Others have set a target in real terms which might appear as 10 or 11% plus inflation. Obviously the arithmetic of compound interest is very attractive; however in practice there is much variation. Indeed current returns from existing operations swing widely with the business cycle and there is no assurance that economic new properties will be found according to someone's arbitrary time schedule. For example, Western Mining Corporation Limited in Australia explored for 30 years with little to show for it, but then found the great Australian nickel deposits and more recently the huge Roxby Downs copper. That long dry spell could not have fitted anyone's arbitrary calendar of growth and yet they would not have found such orebodies without that long period of effort. Should they have abandoned the search? Once a new property has been found or acquired there has to be a threshold rate of return on the new capital to be invested against which to evaluate the property's economics. Conventionally this seems to be 15% after tax, a number common in other heavy industries as well. In some cases it is expressed as a lower number plus allowance for inflation. Discounted cash flow analysis is a very useful tool but it does not make the decision. In the end a "go" decision depends on judgment of many factors some of which are numbers used in the DCF calculation whose credibility must be examined. It is curious how frequently investment proposals come in with the rates of return very close to 15%. The project advocates know that a number much less than 15% will not fly and that a number much more is not necessary. With much higher nominal and real interest rates of recent years, even though before tax, logic suggests that the hurdle rate should also rise. The power of compound interest is so great that 20% is very hard to achieve in any cash flow projection but 18% may be a sensible yard - stick. Once again it is remarkable how many project proposals come in with an 18% return. On the record the mining industry as a whole has not been overly restrictive in choosing its hurdle rates of return. This is shown by the abundance of metals in recent years and the failure of metal prices to keep up with inflation. All of the foregoing is standard text book stuff.

Jan 1, 1985

By E. P. Burtchaell

A method is presented for evaluating the effect of gravity drive upon the reservoir performance of a high relief pool. Conventional forms of reservoir analysis do not consider the alterations in the basic material balance data caused by gravity segregation of reservoir fluids. A procedure is outlined for structurally weighting physical and chemical data for use in the material balance equation. It is demonstrated how actual pool performance data can be utilized to evaluate the future reservoir performance of a gravity drive pool. INTRODUCTION Conventional reservoir engineering. procedure is inadequate for the analysis of an oil pool which has considerable structural relief, steep dips, and good permeability development. In, pools of this type, gravity drainage has an important part in the movement of oil to the wells and the effects of gravity on the overall pool performance should be included in any analysis of reservoir behavior. Many engineers have the opinion that the force of gravity in the movement of oil is not important until the later life of a pool.' Probably the basis for this belief is that gravitational effects may not be readily discernible until a pool is nearing depletion. This would be especially true for pools not having a high degree of structural relief and permeability development. Actually the effects of gravitational forces are at a maximum when the pool pressure is high, for during this period the hydrostatic head of the oil column is at a maximum and the viscosity of the oil is at a minimum. Oil recoveries from pools having favorable gravity drive characteristics may equal or even exceed recoveries which might be expected from water displacement. Field evidence indicates that in some reservoirs gravity drive has resulted in recoveries greater than that which could have been expected from gas expansion or water drive.'.3 Unfortunately, the possible effects of gravity drive on pool performance have been underestimated and other reasons have been sought to explain the high recoveries obtained. There are unquestionably many reservoirs to which the principles of gravity drainage can be effectively applied. It is the purpose of this paper to illustrate one method whereby gravity drive is included in the reservoir analysis of an oil pool. A hypothetical pool, typical of many California reservoirs, is used as an example. As used in this paper, "gravity drive" is defined as the overall effect of gravitational influences on the recovery of petroleum from the reservoir; "gravitational segregation" as the gravity separation of oil and gas within the reservoir; and "gravity drainage" as the downward movement of oil as caused by the force of gravity. SAND VOLUME DATA Fig. 1 presents a structural contour map of the pool under study. Maximum closure is 1950 feet with dips on the south flank approaching 45". The original gas-oil interface was set at -5200 feet. Average thickness of the producing sand was 200 feet. For use in subsequent calculations ill this paper, the pool was subdivided into 100-foot vertical increments and the sand-volume content of each increment was obtained. If the gross sand thickness is small, under 100 feet, the sand-volume content can be obtained by superimposing an isopachous map upon a structural contour map and planimetering the average thickness of each 100-foot increment. For sand thicknesses over 100 feet, one approacli would be to construct a sufficient number of cross-sections of the pool from which the weighted sand-volume of each 100-foot increment could be obtained. Variations in the sand body with depth, as determined by core data, can also be included in the above process. Table I presents a summary of sand-volume calculations, core data, and the original distribution of reservoir hydrocarbons in the pool. Fig. 2 illustrates the structural distribution of the sand-volume content. A total of 171,398 acre-feet is contained within the productive limits of the pool. Assuming an average porosity of 25% and an interstitial water content of 20%, the original hydrocarbon content was computed to be 227,075,000 barrels. DEPTH-PRESSURE DATA The determination of the initial vertical pressure arrangement in the pool is necessary for PVT and material balance calculations. Whenever sufficient data are available, a plot of pressure versus subsea depth of measurement should be made. From this plot a representative fluid pressure gradient can be established. Lacking sufficient initial pressure data, an initial pressure gradient may he estimated or calculated from avail-

Jan 1, 1949

By D. L. Kelley

High-viscosity water injection has been proposed for use in reservoirs containing high-viscosity crude oils. Previous publications have largely ignored the possible effects of the connate water on the proposed process. This paper describes experimental work which indicates that the connate water will be forced ahead of the injected water to form a bank of low-viscosity water. This decreases the oil recovery which would be expected if such a bank were not formed. These effects are shown for a range of fluid mobilities and connate-water saturations for a five-spot injection system. In general, oil recoveries using viscous water are significantly greater than for untreated water even though they are less than would be expected if no connate water bank were formed. INTRODUCTION The effect of mobility ratio on the oil recovery of wa-terfloods has been known for many years. Muskat first pointed out that the fluid mobilities (k/µ) in the oil and water regions would affect the performance of the water-flood, and he estimated the general effect of these variables.' Since this early work, studies of the effect of mobility ratio on secondary recovery have been reported where mathematical,' potentiometric3 and scaled flow models' were used. These studies have shown that a reduction in the mobility ratio between the oil and the displacing fluid would cause additional oil recovery when water-flooding reservoirs containing viscous crude oils. Studies reported by Pye- nd Sandiford 8 have indicated that chemicals to increase injection water viscosity are now available and can be used to reduce the over-all mobility ratio of a waterflood. Where mobility ratios are controlled by the injection of viscous fluids, the connate water of the reservoir can play an important part in the displacement of the reservoir oil. The purpose of this study was to determine the effect of the connate-water saturation in waterfloods where viscous waters are used for injection. DISPLACEMENT OF THE CONNATE WATER Russell, Morgan and Muskat7 were the first to recognize the mobility of connate waters in waterflooding. They conducted waterfloods on oil-saturated cores containing 20 and 35 per cent irreducible water saturations, and found that from 80 to 90 per cent of the "irreducible" water was produced after only one pore volume of water was injected. However, their experiments were conducted at rates of flow significantly higher than those ordinarily occurring in waterfloods. Also, the cores were only from 4.0 to 8.5 cm long. Brown 4 studied a 100-cm linear sand pack which had been prepared to contain connate water and oil. He used 140- and 1.8-cp oils with injection water of essentially the same viscosity as the connate water. He found that all of the connate water was displaced by the injection water in both cases. However, the injection volumes required for complete displacement of the connate water were considerably higher in the case of the more viscous oil. To verify the results of the foregoing experiment, a 10-ft-long linear model was constructed by packing 250-300 mesh sand in a 1/2-in. diameter nylon tube. The model was evacuated, saturated with a brine of 1-cp viscosity, and flooded with a 41-cp mineral oil to the irreducible water saturation of 10.9 per cent. The model was then waterflooded by the injection of a water solution which had an apparent viscosity of 42.6 cp. The solution consisted of 0.5 per cent methylcellulose in distilled water. The viscosities of the oil and connate water were measured with an Ostwald viscosimeter. The viscosity of the polymer solution was calculated by Darcy's law using pressures measured during actual flow conditions. The ratio of the mobility in the oil region to the mobility in the inject ion-water region was approximately 0.32. The mobility ratio of the oil region to the connate-water bank was approximately 14. The mobility ratio between the connate-water bank and the injection water region was 0.024. Approximately 84.5 per cent of the recoverable oil was produced before water breakthrough. Immediately following breakthrough, oil and connate water were produced at an increasing water-oil ratio until the viscous injection water broke through. At viscous-water breakthrough, 96 per cent of the original connate water had been produced. After breakthrough of the viscous water, there was no additional production of connate water or oil. The near-

Jan 1, 1967

By B. J. Kochanowsky

The idea of a portable crusher is not new. Many such crushers are available but they are small and designed for construction work. For many years the author has suggested, both in this country and in Europe, the building of larger portable crushers intended expressly for use in quarries or open pits. Although not applicable under all conditions, there are mining operations where a mobile crusher arrangement could be more profitable than the facilities now used. The primary use of a portable crusher, i.e., a crusher mounted on crawlers or tires, in the rock and mining industries is to reduce costs by permitting the substitution of conveyor belt haulage for truck or track haulage. The usual sequence of operations in surface mining is drilling, blasting, loading, haulage, and crushing. Haulage is normally accomplished by truck or track-mounted cars, the latter method being used for the longer distances. However, by using a portable crusher in the pit, the sequence of operations would be changed so that the crushing stage would occur before haulage (Fig. 1). Such a sequence would permit the use of conveyors to replace the more expensive truck or track haulage methods. Since most quarry and open pit operations normally require a crushing stage, the only additional costs incurred will be due to the investment required to purchase or construct a mobile arrangement for a crusher. But this factor has to be weighed against the advantages to be gained by conveyor haulage. As shown in Fig. 2, transportation of material by belt conveyor over short distances is less expensive than by truck. The inclination of the belt has no effect on belt speed; consequently, the hourly tonnage moved remains the same. Conversely, the output rate of trucks as expressed in tons or ton-miles per shift decreases proportionally to the haulage speed, which is considerably slowed by the steepness of the road (Fig. 3, left). Although maximum possible grades and maximum economic grades of haulage are greater for a belt than for a truck (over the same total lift), the longer haulage distances favor the use of trucks. Although power consumption for hauling on a grade increases for both conveyances, the rate of power consumption increases faster for trucks than for conveyor belts (Fig. 3, right). Since the output rate and related fixed costs are affected by the travel speed, total haulage costs with trucks would increase with the grade more rapidly than the similar costs of conveyor belts (Fig. 4). Travel distance, road grade, speed, size and number of pieces of equipment, efficiency of operation, and many other factors affect such haulage costs. In general terms it can be said that the shorter the distance, the steeper the grade, and the greater the output, the more advantageous the belt becomes in comparison to truck or track haulage. In addition to potential cost savings in haulage procedures, a portable crusher would allow better utilization and performance of shovels. Loading operations would not be interrupted as often by the necessity of waiting for cars or trucks. Unfortunately, the application of belts in open pits for haulage from bench sites is generally not practical under existing conditions because a belt fed directly by a mechanical shovel can be torn, damaged, or worn out quickly by the large rock fragments falling on it during loading. However, by using a mobile crusher this situation can be avoided. As shown in Fig. 1 (b), the shovel feeds rock into the crusher located behind it. The crushed material is initially transported by an extensible and/or movable belt, thence by a longer stationary conveyor to the plant where the material is subjected to further treatment by secondary crushing, screening, etc. The first-mentioned conveyor, needed to bridge the distance between the shovel and the stationary conveyor, is necessarily variable in length owing to the continuous movement of the shovel and the desire to keep the stationary belt at a safe distance from the bench during blasting operations. The remarkable part of mobile crusher operations is the extra-ordinarily high output per man-shift, the low maintenance and power requirements for haulage, and the increased output of the loading shovel. A cement quarry which has been using a portable crusher and conveyor since 1956 requires only three men to operate the shovel and crusher and to transport the crushed rock by belt from the quarry face to the screening plant. If truck haulage

Jan 1, 1961

By Ray W. Carpenter, Ralph Hultgren, Raymond L. Orr

Heats of formation of Au-Zn alloys of compositions Au3Zn and AuZn were rneasured at several temperatures by liquid tin solution calorimetry. The data for Au3Zn show that much smaller heat and entropy effects are associated with the a2- a1 transformation than for the a1 - a transformation. This result is consistent with reported X-ray diffraction studies which indicate that the a1 and (a2 phases have ordered structures which only slightly differ, whereas the a phase is disordered. The heat of formation of AuZn(ß') does not change significantly between 322o and 8000 K, confirming that the ß' alloy remains highly ordered to high temperatures. Not many measurements have been made of the energies of order-disorder transformations in alloy phases. In many cases these energies may be conveniently determined by measuring heats of formation as a function of temperature. Many alloy phases approximately obey Kopp's law of additivity of heat capacities; as a consequence their heats of formation do not change with temperature. In ordered alloys, however, a decrease in degree of order with increasing temperature will contribute an endothermic component to the heat of formation. Thus, for ordered phases where other anomalous contributions to heat capacity in the alloy or its pure components are absent, a measured change in the heat of formation may be attributed to the energy of disordering. A portion of the Au-Zn phase diagram1" is shown in Fig. 1. The composition AuZn(ß') has a well-ordered superlattice structure of the CsCl type at room temperature. The alloy is considered to remain ordered to its melting point,1 although, from X-ray diffraction data and electromotive force measurements of zinc activities, Terpilowski3 concluded that considerable disordering began at -673°K. The Au3Zn composition is a disordered fee solid solution (a) at high temperatures. Between 693o and 543oK, the ailoy orders to form the a1 phase, having a structure based on the L12(AuCu3) type, but with a regular distribution of ordering faults in the c direction, making it tetragonal. Below 543oK, the alloy orders further by means of slight atomic shifts, without diffusion, yielding a complex long-period orthorhom-bic superlattice structure, a2.4-7 Heats of formation of Au3Zn(a2) and AuZn(ß') were measured at 363°K by Biltz et al.8 by aqueous solution calorimetry. Heat capacities of several alloys between 14.9 and 38.4 at. pet Zn in the range from 373" to 723°K were measured by Iwasaki ef al.6 Stoichiometric Au3Zn, however, was not included in the alloys measured, the nearest composition being 26.1 at, pet Zn. Also, the heating rate used for the measurements, 2 deg per min, may not have permitted equilibrium to be reached in the neighborhood of the transformations. Heat contents of solid and liquid AuZn were measured between 488o and 1138°K by Kubaschewski.9 It was decided to determine the ordering energies of Au3Zn and to check on the possible disordering of AuZn by measuring the heats of formation of these alloys at several temperatures. The present paper reports measurements by liquid tin solution calorimetry of the heats of formation of the three phases of Au3Zn at temperatures within their regions of stability, and of AuZn at a low (322°K) and a high (800°K) temperature. EXPERIMENTAL Alloy Preparation and X-Ray Examination. Ten-gram ingots of Au3Zn and AuZn were prepared by melting weighed amounts of gold (99.95 pct Au) and zinc (99.99 pet Zn) together in sealed evacuated Vycor tubes at 1060°K, followed by rapid quenching. The resulting ingots were homogenized at 900°K for 10 days. The Au3Zn alloy lost no weight in preparation, the loss from the AuZn alloy was enough to introduce an uncertainty of only 0.1 at. pet in its composition. Filings for X-ray examination were taken from various parts of each ingot and strain-annealed at 723°K for 20 min in vacuo. The Au3Zn filings were quenched rapidly to retain the disordered phase. Diffraction patterns were determined with an X-ray dif-fractometer using copper Ka radiation.

Jan 1, 1968

By E. P. Burtchaell

A method is presented for evaluating the effect of gravity drive upon the reservoir performance of a high relief pool. Conventional forms of reservoir analysis do not consider the alterations in the basic material balance data caused by gravity segregation of reservoir fluids. A procedure is outlined for structurally weighting physical and chemical data for use in the material balance equation. It is demonstrated how actual pool performance data can be utilized to evaluate the future reservoir performance of a gravity drive pool. INTRODUCTION Conventional reservoir engineering. procedure is inadequate for the analysis of an oil pool which has considerable structural relief, steep dips, and good permeability development. In, pools of this type, gravity drainage has an important part in the movement of oil to the wells and the effects of gravity on the overall pool performance should be included in any analysis of reservoir behavior. Many engineers have the opinion that the force of gravity in the movement of oil is not important until the later life of a pool.' Probably the basis for this belief is that gravitational effects may not be readily discernible until a pool is nearing depletion. This would be especially true for pools not having a high degree of structural relief and permeability development. Actually the effects of gravitational forces are at a maximum when the pool pressure is high, for during this period the hydrostatic head of the oil column is at a maximum and the viscosity of the oil is at a minimum. Oil recoveries from pools having favorable gravity drive characteristics may equal or even exceed recoveries which might be expected from water displacement. Field evidence indicates that in some reservoirs gravity drive has resulted in recoveries greater than that which could have been expected from gas expansion or water drive.'.3 Unfortunately, the possible effects of gravity drive on pool performance have been underestimated and other reasons have been sought to explain the high recoveries obtained. There are unquestionably many reservoirs to which the principles of gravity drainage can be effectively applied. It is the purpose of this paper to illustrate one method whereby gravity drive is included in the reservoir analysis of an oil pool. A hypothetical pool, typical of many California reservoirs, is used as an example. As used in this paper, "gravity drive" is defined as the overall effect of gravitational influences on the recovery of petroleum from the reservoir; "gravitational segregation" as the gravity separation of oil and gas within the reservoir; and "gravity drainage" as the downward movement of oil as caused by the force of gravity. SAND VOLUME DATA Fig. 1 presents a structural contour map of the pool under study. Maximum closure is 1950 feet with dips on the south flank approaching 45". The original gas-oil interface was set at -5200 feet. Average thickness of the producing sand was 200 feet. For use in subsequent calculations ill this paper, the pool was subdivided into 100-foot vertical increments and the sand-volume content of each increment was obtained. If the gross sand thickness is small, under 100 feet, the sand-volume content can be obtained by superimposing an isopachous map upon a structural contour map and planimetering the average thickness of each 100-foot increment. For sand thicknesses over 100 feet, one approacli would be to construct a sufficient number of cross-sections of the pool from which the weighted sand-volume of each 100-foot increment could be obtained. Variations in the sand body with depth, as determined by core data, can also be included in the above process. Table I presents a summary of sand-volume calculations, core data, and the original distribution of reservoir hydrocarbons in the pool. Fig. 2 illustrates the structural distribution of the sand-volume content. A total of 171,398 acre-feet is contained within the productive limits of the pool. Assuming an average porosity of 25% and an interstitial water content of 20%, the original hydrocarbon content was computed to be 227,075,000 barrels. DEPTH-PRESSURE DATA The determination of the initial vertical pressure arrangement in the pool is necessary for PVT and material balance calculations. Whenever sufficient data are available, a plot of pressure versus subsea depth of measurement should be made. From this plot a representative fluid pressure gradient can be established. Lacking sufficient initial pressure data, an initial pressure gradient may he estimated or calculated from avail-

Jan 1, 1949

By C. R. Barrett, N. N. Singh Deo

The high-temperature steady-state creep behavior of a series of dilute copper-silicon alloys was studied to determine the effect of stacking fault energy on the creep-rate. The steady-state creep rate is, when taken at equivalent diffusivities decreases with decreasing stacking fault energy. The stress and temperature dependencies of is suggest that creep is a difusion controlled dislocation climb process. Electron microscopy studies of the creep substructure revealed: 1) the subgrain size is not a function of the stacking fault energy in these alloys, 2) the dislocation density not attributed to the subgrain walls seems to be higher during primary creep and decreases to a lower steady value during steady-state creep, and 3) the dislocation density during steady-state creep decreases with decreasing stacking fault energy. In the past few years numerous investigators have studied the influence of stacking fault energy on high-temperature creep strength. Most of these investigators have confined their attentions to studying the relationship between steady-state creep rate, is, and stacking fault energy, ?, when samples are tested under conditions of comparable stress and temperature. For the case of fcc metals, it was initially shown by Barrett and Sherbyl and since confirmed by many others2"4 that is decreases with decreasing ?, often following an empirical relation of the form i ?m where m is a constant about equal to 3. The application of theory to explain this observation has not been entirely successful. One of the main difficulties has been the almost complete lack of structural information (dislocation density, subgrain size, and so forth) for samples with different stacking fault energies, tested under high-temperature creep conditions. weertman5 has attempted to explain the stacking fault energy dependence of is on the basis of a dislocation climb mechanism. Assuming that both the rate of dislocation core diffusion and the ease of athermal jog formation decreases as ? decreases Weertman has argued that the rate of dislocation climb and hence the creep rate should also decrease as ? decreases. One questionable aspect of Weertman's analysis is the assumption that core diffusion down extended dislocations is slower than core diffusion down unextended dislocations. The only experimental work done in this area, by Birnbaum et al.6 on nickel and Ni-60 Co, has shown the core diffusivity to increase with decreasing ?. Theories of steady-state creep based on the diffusive motion of jogged screw dislocations often seem unable to predict even the qualitative nature of the es- relationship. Assuming that Weertman is correct in his assumption that the dislocation jog density decreases with decreasing ? then the jogged screw theories predict an increasing dislocation velocity with lower ?. It is usually assumed that the increase in dislocation velocity implies a corresponding increase in creep rate. However, two other factors must be considered before such a statement can be made. That is, we must know how both the mobile dislocation density and the effective stress (the difference between applied stress and internal stress) vary with ?. Significant changes in either one of these factors could outweigh any change in dislocation velocity accompanying a change in ?. And with the slower rates of recovery expected in low stacking fault energy materials it seems likely to expect both mobile dislocation density and effective stress to be dependent on ?. Sherby and Burke7 have suggested that stacking fault energy influences the creep rate in an indirect way. These authors cite evidence that the steady-state subgrain size generated during high-temperature creep is a function of ? decreasing with decreasing ?. Assuming the creep rate to be proportional to the area swept out by each expanding dislocation loop and that subgrain boundaries are good barriers to dislocations, then the creep rate should be proportional to subgrain area, hence increasing as ? increases. A critical evaluation of any of the above theories requires more quantitative information concerning the dislocation substructure generated during high-temperature creep. Accordingly this investigation was undertaken with an aim of studying the influence of stacking fault energy on tbe steady-state creep characteristics of a series of dilute copper-silicon alloys. Special emphasis was placed on studying the strain dependence of both the dislocation configuration and density. MATERIALS AND PROCEDURE Dilute copper-silicon alloys of the compositions shown in Table I were tested in tension at constant stress. The relative stacking fault energy of these alloys has been determined and is shown in Table 11. An Andrade-Chalmers lever arm was used to maintain constant stress and testing was carried out in a water

Jan 1, 1970

By J. H. Swisher, E. O. Fuchs

Rate equations were derived to describe the kinetics of internal oxidation of cylinders and spheres. The derived equations for cylinders were checked experimentally by means of sub scale thickness and electrical conductivity measurements on Cu-Cr alloy wires. The properties of the internally oxidized samples were examined with conductivity applications in mind. It was possible to produce uniform dispersions of Cr2O3 in copper with an initial chromium content as high as 3 wt pct. While electrical conductivities only a few pct less than that of OFHC copper were obtained, the Cr2Os particle size and spacing were too large for effective dispersion hardening. T.HE process of internal oxidation has been used widely in basic studies of the permeability of gases in metals. In a review article, Rapp1 has discussed the principles of internal oxidation in considerable detail. From a technological standpoint, internal oxidation is often considered undesirable, since it is a means by which inclusions can be introduced into an otherwise clean material. Another important aspect of internal oxidation is its use as a means of dispersion hardening a material. Broutman and Krock2 discuss this and other methods for making dispersion hardened alloys. The only internally oxidized material known to the authors which is commercially available is a Cu-BeO alloy.3'4 This alloy is made from Cu-Be alloy powder, using a so-called Rhines pack. It has a tensile strength of 80,000 psi and retains its strength at relatively high temperatures. The objectives of the present study were to derive rate equations for the internal oxidation of cylinders and spheres, to check the derived equations for cylinders experimentally, and to examine the structure and properties of internally oxidized Cu-Cr alloys. The Cu-Cr system was chosen for this study because uniform dispersions are obtainable at high alloy contents, which is a desirable characteristic in dispersion hardened materials. RATE EQUATIONS FOR VARIOUS GEOMETRIES A number of authors5--9 have derived equations to describe internal oxidation kinetics. These derivations differ somewhat in mathematical assumptions and approximations, and all except one of the derivations deal exclusively with the internal oxidation of plates. The exception is a brief treatment of cylindrical and spherical geometries given by Meijering and Druy-vesteyn9 as a part of a comprehensive paper on the general subject of internal oxidation. These authors did not obtain rate data to check their derivations, although they did show that the hardness profile across an internally oxidized sample is directly related to the rate of interface movement. For cylindrical and spherical geometries, a quasi-steady-state approximation is needed to circumvent mathematical complications in obtaining a solution to the basic differential equations. In using this approximation, we consider the concentration gradient of dissolved oxygen in the internally oxidized zone or sub-scale to be the same as the gradient which would be present if there were no movement of the subscale interface. The steady-state approximation introduces an error of about 1 pct in computing the rate of internal oxidation of an Fe-1.0 pct Mn alloy plate, if the present method is compared to the more exact method of Wagner.7'10 The details of the derivations of the rate equations for cylinders and spheres are given in the Appendix, and only the results of these derivations are given below. The final equations obtained by Meijering and Druyvesteyn9 can be shown to be equivalent to our Eqs. [1] and [2], although the two approaches are somewhat different. Cylindrical Geometry. [2] where r1 is the outer radius of the cylinder or sphere, cm, r2 is the radius of the unreacted core, cm, see Fig. l(a), D is the diffusion coefficient of oxygen in copper, cm2 per sec, %O is the concentration of dissolved oxygen at the surface of the specimen, wt pct, %Cr is the initial chromium concentration in the alloy, wt pct, and t is the reaction time, sec. Plate Geometry. The analogous rate equation for a plate has been derived previously for internal oxidation of Fe-Al alloys.8'11 For Cu-Cr alloys, we may write the same equation as follows: [3] where r1 is the half-thickness of the plate, cm, and r2 is the distance from the mid-plane to the subscale intherate is An analysis of Eqs. [1], [2], and [3] shows that for a plate the rate is completely parabolic. The initial

Jan 1, 1970

By J. Chipman, G. G. Hatch

One of the important functions of the iron blast furnace is the desulphur-ization of pig iron before it enters the steelmaking furnaces. However, the increasing concentrations of sulphur in the metallurgical coke, source of approximately 90 pct of the sulphur present in the blast furnace charge, and demands for higher rates of production within recent years have increased the need for greater desulphurization within the iron blast furnace. Furnace operators are beginning to look for desulphurizing agents other than blast furnace slag to accomplish the desired degree of desulphurization. A considerable amount of work has been done on desulphurization outside the furnace with soda ash, calcium carbide and various synthetic slags. Whether the desulphurization of pig iron is accomplished wholly inside the furnace or partly inside and the remainder outside, will be determined by the economics involved. Regardless of which is the case, it is believed that it is necessary to have a better understanding of the physical chemistry of desulphurization by blast furnace slags. To this end, it is the object of the present investigation to attempt what is believed to be the first equilibrium study of the distribution of sulphur between liquid pig iron and a wide range of blast furnace slag compositions. Review of Literature There is a considerable amount of information in the literature concerning the desulphurizing power of iron blast furnace slags, the solubility of various sulphides in the slags, and the effect on desulphurization of temperature, of elements dissolved in the liquid iron, and of viscosity. However, there is nothing to indicate that the equilibrium distribution of sulphur between liquid iron saturated with carbon and iron blast furnace slags has been studied experimentally. Wentrupl has made probably the most detailed study of the desulphurization of pig iron to date. He considered that there are three distinct aspects involved, namely: 1. Desulphurization within the blast furnace (by lime and manganese). 2. Subsequent desulphurization by manganese. 3. The effect of subsidiary reactions on the desulphurization by manganese. The experimental work carried out by Wentrup was devoted mainly to obtaining a better understanding of how desulphurization by manganese was accomplished in the mixer and the ladle. Particular attention was given to the part played by carbon, silicon, and phosphorus associated with manganese in the iron, and the effect of temperature on desulphurization. The experimental results indicated that desulphurization by manganese is purely a process of crystallization of manganese sulphide. The addition of silicon to iron melts containing 3.5 pct carbon and less than 0.5 pct manganese had no noticeable effect on desulphurization, but with 1-2 pct manganese the silicon additions improved the desulphurization. Additions of phosphorus also resulted in improved desulphurizati011 by manganese, but the effect was not as marked as in the case of silicon. It was also found that desulphurization by manganese was further improved by lowering the temperature. In order to explain desulphurization inside the blast furnace, Wentrup considered the system iron, sulphur, calcium, oxygen, manganese. (silicon). The distribution of sulphur between the metal and slag was represented by the following equation: (SS) _ (S)Fe + (S)Ca + (S)Mn .... [S] = [s] [1] The parentheses and the brackets represent the equilibrium concentrations in weight per cent of the slag and metal constituents, respectively. Since FeS D (FeS) _ (FeS) (S)Fe LfeS - [FeS] [S] [2] (CaO) + S e (FeO) + (S)Ca _ (FeO)(S)Ca. (S)Ca _ (CaO) Kl = (CaO)[S] [S] ~K1(FeO) [3] Mn + S D (S)Mn (S)Mn (S)Mn K' = [MnpT "1ST = *lIMnJ !4) Substitution of Eq 2, 3, and 4 into Eq 1 resulted in if = L- + *> (Sol + K^ (S) [51 Eq 5 was used to calcu1;lte -f^j and [S] at 1480°C for slags containing 30-50 pct lime, 0.1-2.5 pct iron oxide, 0-26 pct silica, 2 pct sulphur and iron analyzing 1.5 pct manganese. The value for LFaB at 1480°C was found to be equal to 4.5, based on the experimental work of Bardenheuer and Geller.2 The results of the calculations are shown in Table 1. Although the slags are hypothetical and do not represent the range of compositions found in ordinary blast furnace practice, the calculations indicate that lime is effective in controlling desulphurization only if the iron oxide and silica contents of the slag are kept low. Schenck3 did not claim K1 to be a true equilibrium constant, but an empirical value which varied with the silica content of the slag.

Jan 1, 1950

By H. C. Theuerer

GeC14 may be purified by extraction with HCI and chlorine. The process is as effective for the removal of AsCI:, as the more cumbersome distillation methods usually used for this purpose. GERMANIUM for semiconductor use contains impurities at levels no higher than a few parts in ten million. Material of this quality is obtained from highly purified GeC1, by hydrolysis to the oxide and reduction of the oxide in hydrogen. When purifying GeCl,, AsC1, is the most difficult impurity to remove. This is usually accomplished by multiple distillation procedures.1-3 AsC1, may also be removed from GeC1, by extraction with HC1.1-4 Reducing the arsenic to low concentrations is not practicable, however, because of the large number of extractions needed. This paper discusses a new method for the removal of arsenic from GeC1, by extraction with HC1 and chlorine. The method is rapid, leads to little loss of germanium and is at least as efficient as the distillation procedures currently being used. Theory of Extraction Procedures In the simple extraction of GeC1, with HC1, the following reaction occurs ASCl8G8C1 D AsCl3rc1 at equilibrium CA/Cn= K, where K is the distribution coefficient, and C, and C,, are the molar concentrations of AsC13 in HC1 and GeCl,, respectively. The materials balance equation for this reaction is VACA + vncn = VnC,, where V, and Vn are the volumes of HC1 and GeCl4, respectively, and C, is the initial concentration of AsC13 in GeC1,. From this it can be shown that for multiple extractions where C,, is the concentration of AsC13 in GeC14 after n extractions, and r is the ratio of V, to V,,. It is assumed that r is maintained constant, that equilibrium is established during each extraction, and that K is independent of the AsCl3 concentration. By saturating the system with chlorine, the following reaction occurs in the aqueous phase AsCl3 + 4H2O + Cl2 D H5AsO4 + 5HC1 at equilibrium K' = ------------ ai - a4 h2u - aet2 where a is the activity of the various components. The effect of this reaction is to reduce the concentration of the AsC1, in the aqueous layer and, therefore, to promote further extraction of this component from the GeC1, layer. If the arsenic acid remains entirely in the aqueous phase, the net effect of this reaction is to promote the removal of arsenic from the GeC11. The materials balance equation for extraction with HC1 and chlorine with the foregoing reaction is, then, VaCC + VACA + VACn = VnCo where C,. is the molar concentration of H3AsO, in the HC1. With the added assumptions that the activities of AsC13 and H8ASO4 in the aqueous phase are equal to their molar concentrations, it can be shown that for n extractions Cn/Cu = (1/rkK + rK + 1) n where k - K1 a4h2o - acl2/aoncl. It can be seen by comparing Eqs. 1 and 2 that if k is large, the removal of AsC1, by HC1 extraction will be greatly improved by the addition of chlorine. Dilution of the HCI used in the extraction with chlorine would also favor the separation. This, however, would increase the loss of GeCl,, which is undesirable. Experimental Procedure Germanium prepared from oxide of semiconductor purity is n-type with resistivities greater than one ohm-cm. The resistivity is controlled by the donor concentration, which is —lo-: mol pct. Germanium prepared from material with added arsenic will have lower resistivity commensurate with the arsenic concentration. With such material, at arsenic concentrations above 10-1 mol pct the resistivity is controlled by the added arsenic, and effects due to other impurities initially in the oxide are negligible. In this investigation GeO, of semiconductor purity was converted to GeCl,, and -0.01 mol pct As was added. This material was used for the extraction experiments and the purification attained determined by a comparison of the resistivity data for samples of germanium prepared from the initial and purified GeC1,. A method for calculating the arsenic concentration from the resistivity data is discussed later. The details of the experimental procedures used are as follows: Two hundred and thirty cu cm GeC1, were prepared by the solution of GeO, in HC1, followed by

Jan 1, 1957

By G. W. Nabor, A. S. Odeh

The effect of production history of a well on the results of two-rate flow tests, and conventional build-up analyses was investigated. The effect was examined by means of digital computers and an R-C network model, respectively, for wells with infinite and finite radii of drainage. For systems which behave as infinite, it was found that production history and the duration of production at constant rate prior to the initiation of the test have important effects on the results. During build-up time equal to about one-fourth of the stabilized time, correct permeability-thickness product calculations can be made. For wells with finite radii of drainage, the time was determined during which the straight line can be satisfactorily used for permeability-thickness product calculations in case of drawdowns and build-ups. On build-ups, the dimensionless time (based on the external radius) during which the straight line gives reliable results was detertriined to be about 1/12. This is one-fourth as long as that of the drawdown. The investigation was done theoretically, and subsequently was verified by R-C network model runs. General interpretive rules were formulated which, if not followed, could lead to serious errors. Moreover, a recommended testing procedure is reported. INTRODUCTION The method used by most reservoir engineers for estimating formation characteristics in a producing well is the analysis of pressure build-up data. The method originally devised by Horner' makes use of the point source solution to the diffusion equation. This solution is approximated by a logarithmic function and the superposition principle is employed to arrive at the well known pressure build-up equation:where q, the flow rate, is in reservoir B/D; ft is in cp; kh is in md-ft; At is the shut-in time; and t is the producing time. At and t are in any consistent time units. Ey. 1 is applicable to a well of unlimited drainage radius which produces at a constant rate q from zero to time t and is then shut in. Such a constant production rate seldom obtains in practice. Therefore. a correction term must be applied to Eq. I to account for the varying rate. Two theoretically accurate methods are available for treating the variable rate case. The first, originally derived by Horner,' is based on the application of the superposition theorem. It requires knowlege of production history of the well as a function of time and results in lengthy and laborious calculations. The second t*q* niethod is suited for short production tests and requires that the shut-in time be at least one and one-half times the production time. A third method which is not based on any theoretical justification and which was suggested by Horner as a "good working approximation" is the one used by the majority of reservoir analysts. especially when the well has been producing for a long time and the t*q* method is not practicable. The key to this method is in choosing or determining the t that appears in Eq. 1. Horner suggested using a corrected time t, in place of t. t, is calculated by dividing the total cumulative production by the last established rate. Therefore, a normal procedure of pressure build-up testing is to stabilize the well at a constant rate for at least 24 hours before shut in and to use the stabilized rate to calculate t,. The analysis is then made by plotting either or P. and examining the resulting plot for the expected straight line to calculate kh and the original reservoir pressure. Recently, Russell" proposed a method for determining formation characteristics from two-rate flow tests. His method reduces to pressure build-up if the second flow rate is zero. Russell uses the Horner simplified procedure for calculating a corrected t,. His method also requires the stabilization of the well at a constant rate q which is used to calculate t,. Theoretically, the above procedure is valid for a well with an unlimited radius of drainage or with a limited radius as long as the boundary effect has not been felt by the well. Several authors'' derived formulas which allow the estimation of time during which limited reservoirs behave as infinite ones and, thus, can be treated by unsteady-state mechanics. One equation derived by Swift and Kiel' terminates the application of unsteady-state theory when the drainage radius reaches one-half the reservoir radius. Thereafter, steady-state behavior obtains. Another equation derived by Jones'" initiates steady-state

Jan 1, 1967

By L. A. Wilson, P. J. Root

The relative costs of heating a reservoir by steam injection and by combustion have been examined. The comparison was based on a model similar to that proposed by Chu.' The cost of boiler feed water, the price of fuel, pressure and plant capacity were parameters in determin-ing the costs of air compression and steam generation. The analyses compare the cost of heating to the same radius by the two methods. Results suggest that the two primary factors for comparison are the price of fuel and the amount of crude burned during underground combustion. The cost of fuel has a greater effect on the cost of heat from steam than it does on its cost by combustion. As a result, analyses indicate that when the price of fuel is low, steam may be unequivocally cheaper than air. The influence of heat loss is such, however, that as the heated radius increases combustion becomes relatively more competitive depending upon the amount of crude burned. This implies that steam may be cheaper for small stimulation jobs (huff and puff) but combustion may be more economically attractive for heating large areas (flooding). INTRODUCTION Use of thermal methods of recovery is an accepted fact today. After an induction period of several years, processes are being widely used that involve reservoir heating to augment recovery. Of the several techniques, steam injection and forward combustion appear to be destined to dominate the field. Although the objectives of both are the same, the basic differences between generating heat in situ and injecting heat after surface generation influence the cost in different ways. This study compares the cost of heating a reservoir either by steam injection or by forward combustion. There has been no consideration of recovery. Presumably, recovery from the swept region would be high in either case. The sole consideration was the cost of heating to the same radial distance by either process. PROCEDURE THE MODEL The basis for comparison was a mathematical model similar to that used by Chu' for combustion. The model simulates a radial heat wave in two-dimensional cylindricaI coordinates. It includes heat generation, conduction and convection within the reservoir and conduction in the bounding formations. Thus, heat losses from the formation are considered. Three significant modifications were made. 1. Equal logarithmic increments rather than equal increments were used for the mesh spacing in the r direction. By this technique large distances were simulated with relatively few mesh spaces. 2. A backward difference approximation to the convection term was used to avoid troublesome oscillations which result from a central difference approximation when the convection term is large. 3. The radial increments of the combustion zone motion were not necessarily uniquely related to the mesh configuration. The cumbersome step function introduced by the heat of vaporization of steam was circumvented by assuming the enthalpy of the steam to be a linear function of temperature between reservoir temperature and steam temperature. This is equivalent to assuming an average heat capacity numerically equal to the difference between the enthalpy of saturated steam and the enthalpy of water at reservoir temperature divided by the difference between the two temperatures. Heat losses obtained by this model are in essential agreement with those obtained by the analytical solution of Rubenshtein.' A detailed description of the model is presented in the Appendix. Using the model, the times required to heat to particular radial distances were obtained as a function of injection rate and other physical parameters. For the steam case, injected fluid was assumed to be saturated steam at pressures of either 500, 1,000 or 1,500 psia. The corresponding temperatures are 467, 544 and 596F, respectively. Thickness ranged from 10 to 50 ft and injection rate ranged from 100,000 to 1 million Ib/D. Reservoir and overburden temperatures at the injection well were assumed to be that of saturated steam at the injection pressure. The effect of maintaining the overburden temperature at the well at a different temperature (initial reservoir temperature) was examined with no significant change in behavior. The influence of wellbore heat losses for the steam case was determined in the following manner. The rates of heat loss as a function of time were estimated using an approach similar to that suggested by Ramey." he data were based on injection through 2%-in. tubing in 7-in. casing. Integration of these data over the entire iniection period yielded the total heat loss. Total heat losses were then corrected to their equivalents in steam (this number resulted from dividing the total heat loss by the latent heat). This was considered additional steam required to accomplish the reservoir heating and the total cost was increased accordingly.

Jan 1, 1967

By J. E. Flinn, S. A. Duran

The steady-state creep behavior of poly crystalline cad mi inn was studied over a temperature range of (1.56 to 0.94 Tm. Two distinct mechanisms were found to occur over this temperature range. They were described by: where and represerqt the minimum strain rates corresponding to the low- and high-temperature regions, respectirely. The two regions of constant acti11ation energy were connected by a transition region where the strain rate was controlled by both mechanisms acting in parallel. At temperatures below a transition temperature of about 0.7 Tm the agreement between the activation energy value for creep and that for self-diffiision suggests a rate-controlling mechanism of dislocation climb. For cadwzium, steady-state creep at temperatures above 0. 7 Tm appears to be controlled by another mechanism, perhaps involving the behavior of dislocation jogs. FRENKEL et al.1 studied the high-temperature creep of polycrystalline cadmium and reported an activation energy of 21 kcal per mole for the 0.5 Tm < T < 0.8 Tm range. Based on observations of creep rate at only two temperatures, a value of 22.1 kcal per mole was determined by Medbury. These two investigations were for the purpose of showing agreement between the activation energy for creep and that for self-diffusion, reported3 as 18.2 and 19.1 kcal per mole, respectively, for diffusion parallel and perpendicular to the hexagonal axis. Gilman4 investigated prismatic glide in single crystals of cadmium over a higher-temperature range of 0.72 to 0.93 Tm, and found an activation energy of 29 kcal per mole. He also reported5 an activation energy higher than that of self-diffusion for prismatic glide in zinc single crystals deformed at temperatures near the melting point. This value was in good agreement with those found for an equivalent temperature range by Flinn and Munson6 and by Tegart and sherby7 for polycrystalline zinc. These two independent studies also disclosed at lower temperatures another value of activation energy near that for self-diffusion. It would be expected from the creep results on zinc and single-crystal cadmium that creep studies on polycrystalline cadmium, extended to temperatures near the melting point, might yield an activation-energy value higher than the 22 kcal per mole value found in earlier studies. The purpose of this paper is to report the steady-creep behavior of polycrystalline cadmium over a temperature range of approximately 0.5 to 0.9 Tm EXPERIMENTAL METHOD The cadmium used in this study was obtained in the form of as-cast rods, 0.5 in. diam, through the courtesy of the Bunker Hill Mining Co. The material was of 99.995 pct purity, as determined by spectro-chemical analysis. The creep specimens, which were 0.250 in. diam by 0.400 in. long and annealed at 300°C for 45 min to produce a stable average grain diameter 0.25 mm, were tested in compression using an apparatus similar to that described by Sherby.8 The specimen temperature was controlled to within ±0.5°C with the help of appropriate constant-temperature baths. The applied stress was maintained within 1.0 pct of the desired value by the additions of lead shot at fixed strain increments. No barreling was observed over the strains encountered during testing. Isothermal creep tests9 were used in the study with only a few differential temperature tests10 run for comparison purposes. Steady-state creep data were obtained over a temperature range of 60 to 287°C (0.56 to 0.94 Tm) at five stress levels ranging from 28.1 to 140.6 kg per sq cm. RESULTS The minimum or steady-state creep rate may be described by an equation of the following form:" where i is the minimum strain rate, S is the structure factor, F is a stress function, Qc is the energy of activation, T is the absolute temperature, and R is the gas constant. The minimum strain rates obtained in this study for cadmium were recorded on a semilogarithm plot as a function of the reciprocal absolute temperatures for the various stress levels, as shown in Fig. 1. This figure shows a characteristic transitional behavior" with a parallel interaction of two mechanisms. It is obvious that the activation energies corresponding to the individual processes are insensitive to stress because the curves are parallel. The discrete activation-energies values for the low- and high-temperature regions for the various stress levels are reported in Table I, and were determined by the least-mean-square method. For the low-temperature region, an activation energy of 20.7 ± 0.6 kcal per mole was obtained, and for the

Jan 1, 1967

By R. H. Jacoby, V. J. Berry

When dry gas is injected into a reservoir containing a volatile crude oil, a significant amount of the reservoir liquid phase may become vaporized. The resultant rich gas phase, when subsequently produced, con-tributes to tank oil production. This contribution assumes greater importance, the more volatile the oil in-olved. Oil recovery may be sub-stantially greater than that predicted by conventional frontal-drive methods, which do not consider the vaporization equilibrium between the reservoir oil phase and the injected gas. A calculation method has been (developed to account for vaporiza-. tion of the reservoir liquid phase during gas-injectiorz operations, and for the tank oil production which results from this factor. Recovery performance calculations are presented for a reservoir containing a highly volatile oil. Tank oil recovery is calculated to be about twice that predicted by the use of the conven-rional frontal-drive equations. In contrastto usual pressure maintenance performance results, in which the gas-oil ratio rises at an increasing rate after gas breakthrough, the pre-dicated gas-oil ratio rises rapidly to about 12,000 scf/bhl and then rises much less rapidly. During gas inje-tion, most of the reservoir liquid phase contacted is evaporated by the dry injection gar. The gas-oil ratio during this period is dependent upon reservoir pressure. The higher the operating pressure, the lower the gas-oil ratio. The predicted behavior i., in accordance with laboratory PVT tests made to .simulate the vaporization behavior. In addition to recovery performance predictions, results of the calculation procedure include complete wellstrearn composition data of value in the design of gasoline plant facilities often used in con-rrction with gas-injection operations. INTRODUCTION In the cycling of gas-condensate reservoirs, dry gas is injected to maintain reservoir pressure during wet gas production and to thereby eliminate or reduce ultimate loss of liquids due to retrograde condensation within the reservoir. Gas injected into crude oil reservoirs has a dual function. It displaces oil to the producing wells and at the same time serves to partially or fully maintain reservoir pressure. Oil shrinkage which would occur upon pressure reduction is thereby minimized or eliminated. Accepted calculation methods are available&apos;= for predicting recovery performance of either gas-condensate reservoirs or crude oil reservoirs which are being subjected to gas injection. In gas-condensate reservoirs, any retrograde liquid formed does not flow, and it is necessary to account only for the vaporization equilibrium between this liquid and the injected gas. Conversely, in normal crude oil reservoirs, both the oil and gas phases flow, but it has not been considered necessary to account for any vaporization of the reservoir liquid which might occur upon contact with the dry injection gas. Recently, high shrinkage reservoir fluids known as "volatile oils" have been found in increasing amounts. These oils are characterized by tank oil gravities above 4.5" API, solution gas-oil ratios above 1,000 scf/bbl, and reservoir volume factors above two. Special techniques have been devised for predicting depletion performance of reservoirs containing such oils.74.5.1; One of the characteristics of reservoirs producing volatile oils is that the reservoir gas phase carries a significant amount of oil which is recoverable as stock-tank liquid. This unusual vaporization be-havior implies that an appreciable amount of reservoir liquid would be vaporized upon contact of the oil with dry injection gas. In a gas-injection operation, tank oil recovery would he obtained not only through frontal displacement of the reservoir liquid by the injected gas, but also through production of the rich gas phase. This means that not only are improved methods needed to predict recovery of such oils from reservoirs undergoing gas injection, but it would also be expected that high oil recoveries might be obtained by such operations. When relatively dry gas is injected in to a volatile oil reservoir, phase equilibrium between the injected gas and the reservoir oil will tend to bc established. Initially the most volatile components, such as propane, the butanes and the pentanes, will account for most of the material transferred from the oil phase to the gas phase. As the partially stripped oil phase is contacted with additional dry injection gas, the heavier intermediate components, such as the hexanes, the heptanes and the octanes, will gradually transfer to the gas phase in increasing amounts. This is because the supply of lighter components in the oil phase dwindles due to the stripping action of the injected gas. This stripping action will usually continue to be effective down

By G. E. Perry, J. F. Bruskotter, D. G. Russell, J. H. Goodrich

The depletion performance of gas wells has been investigated by mathematical simulation techniques. The gas well model which was studied consisted of a single well located in the center of a bounded, cylindrical, homogeneous reservoir. Dependency of gas compressibility and viscosity on pressure was considered in studies of well performance at both constant mass flow rate and constant flowing pressure conditions. To carry out the investigation, the nonlinear, second-order, partial differential equation which describes Darcy flow of a nonideal gas through porous media was solved numerically. Some of the previous investigations of gas well performance have been of limited general use, because assumptions were introduced to simplify either the gas properties or the basic differential equation. Other studies have been rigorous in these respects but have presented a very limited set of calculated results. The present study was attempted to present a rigorous theoretical model and sufficient numerical results to permit meaningful conclusions to be drawn. It was found that all terms must be retained in the partial differential equation to make accurate predictions. The neglect of higher-order terms, e.g., terms of the order of the "gradient squared", leads to serious material balance errors at large times and to conservative estimates of gas well performance. The higher the gas flow rate and/or the lower the permeability-thickness product of the formation, the more pronounced are these deviations. For example, in a well draining 640 acres in a 25-md-ft formation (8,120 MMcf gas in place) at a constant rate of 993 Mcf/D, the rigorous solution predicts a bottom-hole pressure decline from 4,000 to 1,000 psia in 8.7 years. If higher-order terms are neglected in the differential equation, this decline in pressure is predicted to occur in 5.3 years. With the results of the numerical solution of the differential equation as a basis, simple, easy-to-use approximations for predicting gas well performance for Darcy flow conditions have been developed. These simple approximations are based on the familiar equations for flow of a single, slightly compressible fluid. The approximate methods possess a high degree of accuracy and enable the prediction of long-term gas well performance to be made quickly and accurately without the use of a digital computer. Both transient and stable flow period approximations were developed. INTRODUCTION In recent years income from the sale of natural gas and associated products has represented an ever-incre as ing fraction of the industry&apos;s total revenue from operations. To meet the surge in demand for natural gas, the industry has depended heavily upon established reserves and has actively pursued development of new reserves. The search has progressively led to reservoirs which yesterday were too tight and/or deep to yield the desired return on invested capital. More than ever before, evaluation accuracy is now required to forecast the criteria upon which engineering recommendations and management decisions are based. Considerable effort has been expended by both research and operations personnel on the development and application of methods for analyzing and predicting the performance of gas wells. Fundamentally, the problem is the familiar one of extracting data during the drilling, testing and early production life of a well and applying these data within an accurate simulation model to predict long-term behavior. During the past 30 or more years a voluminous literature dealing specifically with gas field problems has been generated. A recent book&apos; lists a comprehensive bibliography of published material through 1959. Over 1,200 references are cited. Since then 39 additional articles on natural gas technology have been published in Transactions volumes of the Society of Petroleum Engineers of AIME. Most existing theory for predicting gas well performance requires that one or more idealizations (e.g., steady-state flow, ideal gas of constant viscosity, small and constant compressibility and constant-viscosity fluid) be applied. Although existing theory may apply directly or be adapted by various artifices to describe specific gas well and reservoir behavior, no widely applicable method is available, and existing methods appear to be subject to appreciable error unless better limits of applicability are defined. The objectives of this paper are (1) to present numerical solutions to the partial differential equation describing gas flow under conditions of general interest in gas well performance prediction work, (2) to present solutions which possess a high order of accuracy for both transient and stabilized flow periods of a well producing at constant rate or constant pressure, and (3) to develop and present simplified relationships which can be used as high-order approximations to the exact numerical results for fast and accurate predictions of gas well performance at the operating level. Combined, these objectives are designed gen-

Jan 1, 1967