Search Documents

By Paul B. Crawford, Bobby L. Landrum

Electrical model studies have been made of the effect of drainhole drilling on production capacity. Drain holes are those lateral holes which are sometimes drilled to assist in increasing the oil drainage into the main wellbore. The studies indicate that if the drainhole length i.r near 0.2 of the drainage radius, the increase in production capacity may be approximately 1.8, 2.08, trnd 2.33 for the case of one, two and four drain holes, respectively. If multiple drain holes are drilled, shorter holes are required to achieve a specific increase in production capacity than if one or two holes are drilled. To achieve increases in production capacities by factors near three it is shown that the drainhole lengths divided by the drainage radii must be approximately 0.64, 0.48, cnd 0.37 for the case of one, two and four drain holes, respectively. IN T RO DUCTION Kecent developments in drilling techniques have permitted the drilling of lateral or horizontal holes to increase the rate of production. The number of these drain holes may vary from 1 to 10, and they have been drilled laterally as much as 100 ft. Drainhole drilling at substantial depths is of recent origin but has been used in the Spraberry field of West Texas, the La Pas field in Western Venezuela, and certain fields in California. The diameter of the drain holes has been approximately 6 in. In drilling drain holes it is desired to know how much increase in production capacity may be expected, what additional benefit may be gained by drilling more than one drain hole, and what increase in production capacity may be brought about by lengthening the existing drain holes. In realization of the importance of lateral drilling, electrical model studies were made to provide quantitative estimates on the effect of drain holes on production capacity. ANALYSIS An electrical model was used in studying this problem. A wooden circular reservoir approximately 30 in. in diameter and 2 in. deep was filled with a dilute salt solution to a depth of 1 in. A copper band formed the periphery of the tank, and a copper wire was inserted in the center of the tank to represent the well. The drainage radius divided by the well radius was equal to 500. A diagram of this type model is shown in the inset of Fig. 1. An electrical potential was applied to the wire well and copper band at the periphery to determine the flow of current. Following this measurement copper wires of various length were soldered to the vertical wire well in such a manner as to simulate lateral wells drilled into the reservoir. Multiple lateral holes of various lengths were studied by varying the number and length of copper wires soldered to the vertical wire well. All multiple lateral wells were located symmetrically about the vertical well. The effect on production capacity by drilling an infinite number of lateral holes was studied by replacing the lateral hole with a copper disk. For each condition the flow of current was measured and compared with the flow of the vertical well only. The model which is thus described provides the steady state increases in production capacity which may be expected in some cases. The data conform to those which would be observed for a uniform homogeneous reservoir. The data may have primary application to those reservoirs producing above the bubble point. It is realized that few lateral holes are drilled in reservoirs of this type, but the data may be of some value in guiding future work or aiding to explain observed phenomena.

Jan 1, 1956

By J. R. Kyte, V. O. Naumann, C. C. Mattax

Imbibition and waterflooding tests at reservoir conditions show that conventional tests at room conditions sometimes lead to recovery estimates that are too low. Core samples which are of intermediate wettability at normal test conditions may be strongly water-wet at reservoir conditions. This was found to be true for samples from two of three reservoirs investigated. For the other reservoir, less difference was observed between the wettability of samples at normal test conditions and at reservoir conditions. Samples which became more water-wet at reservoir temperature and pressure also water flooded more efficiently at reservoir conditions. These results show that core analysis tests on intermediate wettability samples should be performed at reservoir conditions to obtain representative data. INTRODUCTION Throughout the oil industry, core analysis tests normally are performed using refined oils and at room temperature and pressure. This procedure has always been subject to question because the possibility is generally recognized that the reservoir environment may influence displacement behavior. Burkhardt, et al,' have performed waterflood tests at reservoir conditions on East Texas cores. They found little difference between results obtained at reservoir conditions and results of conventional flooding tests at room conditions. As they pointed out, however, this situation might not be general since they obtained data on only one rock and crude-oil system. This paper presents data which demonstrate that it is sometimes necessary to perform core analysis tests in a simulated reservoir environment to properly represent the reservoir wettability and displacement characteristics. The work is specifically concerned with reservoirs which were thought to be of intermediate wettability on the basis of conventional test results at room conditions, i.e., core samples imbibed only small amounts of either oil or water and exhibited no decided preferential wettability characteristics. Tests were performed on cores from three such reservoirs to determine their behavior in the reservoir environment. lmbibition and waterflood tests were used to demonstrate the effect of reservoir environment on reservoir rock wettability and displacement characteristics. Results of tests performed at reservoir conditions were compared with results from conventional tests. Imbibition tests were performed on samples from all three reservoirs. Waterflood tests were performed on samples from only one reservoir. Fluid properties for the three reservoirs studied are shown in Table 1. Tests at reservoir conditions were performed using bottom-hole crude oil at reservoir temperature and pressure. The aqueous phase used in these tests was a synthetic brine approximating the reservoir brine composition. Conventional tests were performed at room conditions with refined oils and synthetic reservoir brine. The oil-to-water viscosity ratios were adjusted to the reservoir values. The cores used in all the tests were taken from the reservoirs using either brine alone or a bentonite-car-boxymethylcellulose mud without additives as coring fluids. Previous tests show that these fluids have no significant effect on the wettability of reservoir rock samples.1 The cores were preserved in brine at the wellsite to prevent wettability changes during shipment and storage. Wellsite imbibition tests have been performed on fresh samples from several reservoirs (including Reservoir A) and compared with later laboratory imbibition tests on samples preserved in brine. The comparison showed that the wettability of the preserved cores was the same as that of fresh, surfaced cores at the wellsite.2 IMBIBITION TESTS Conventional imbibition tests were performed on preserved cores using the apparatus and procedures described by Bobek, et al.2 A sketch of the apparatus for imbibition tests at reservoir conditions is shown in Fig. 1. The preserved samples were first centrifuged to expel excess water and then placed in a stainless-steel container which served

By Jan H. Hommen

TRENDS IN MINE FINANCE Reviewing historical trends may help in dealing with future funding strategies. The first trend is the increase in the use of debt to fund corporate America. Surprisingly, when comparing that with other corporates in the global marketplace, the U.S. does not look too poorly. Japanese, in particular, but also German companies, are much more leveraged than their U.S. counterparts. Note that in Japan and Germany the government and the consumer have shown substantial restraints when it comes to credit. Not so in the U.S. where the savings rate (4%) is much lower than that in Japan (20%) and in Germany (14%). As a result, the cost of leverage is substantially lower in Japan and in Germany than in the U.S. Lower debt costs in these countries allows one to obtain more of it for the same price. The close relationships between financial institutions, corporations and governments in Japan and Germany also make the risk of leverage more bearable and understandable. Also note that interest rate volatility in these countries is more subdued. The investor base in common stock is also changing. Private pension funds have been net sellers of equities during the last few years while public pension funds have been net buyers. Individuals are basically out of the market. It is estimated that over $330 billion of equity was retired in the U.S. through 1986-1988. The reasons are well known: restructurings, buy-backs, takeovers, mergers, and buy-outs. This de-equitization is likely to continue as the cost of equity capital is higher than the cost of debt capital, even on a before-tax basis. Because there has been no market for equity since October 19, 1987, junk bonds have replaced equity capital in the corporate finance structure. The non-investor-grade borrower has access to the public markets through high-yielding debt securities that often look like senior equity. Deregulation has also provided 24-hour markets; securitization can now liquify assets that before were rusting on balance sheets; information processing and communication tools give facts that are instantaneously available, creating the need for products to cover risk and exposures unknown before the start of the information and technology explosion. Thus we now have swaps, options, futures, options on futures, caps, as well as program trading, portfolio insurance, index funds, etc. Globalization will be unstoppable. Today's environment demands it. Corporations are now global competitors, even though they may have never shipped a product off-shore. Boundaries are disappearing and impediments have been removed. Stocks can be issued in a global offering with distribution in all major financial centers of the world. The high leverage of Japanese companies has made common stocks very scarce, which with the high savings rate, makes equity capital inexpensive in Japan. If markets were perfect everyone would issue stock in Tokyo and buy it back in New York. This would certainly stop the bleeding on the balance of payments. With the deficit balance of payments in America, it is no surprise that the cost of money is expensive here. The risk-return relationship on U.S. equity is not favorable enough to appeal to Japanese equity investors. The Tokyo-market, helped by a stronger yen, surpassed New York as the most important equity market in the world. The question is, will the Tokyo market hold up when Japanese domestic consumption increases at the expense of savings? In the U.S. debt markets, the plain vanilla debenture is making a comeback after a period of absence. Commercial paper has replaced short term bank borrowings for working capital funding by larger corporations. The banks have made an interesting comeback with loan sale programs that will continue to attract borrowers. There will be a continuing

Jan 1, 1990

By W. Rostoker, H. Nichols

It has been demonstrated in previous work1'2 that wetting of aluminum alloys by liquid mercury can cause fracture to occur with substantial suppression of prior plastic flow. This has been interpreted as the action of a critical liquid species in reducing the cohesive strength of a stressed solid at sites of potential crack nucleation. In this way the critical transverse stress for fracture nucleation is reached at much lower stresses than is normally the case. The stresses for both fracture nucleation and fracture propagation are depressed by the same mechanism. In so doing the plastic-energy component of fracture energy is reduced to a level which is of the same order of magnitude as the surface, interfacial, or cleavage energy. Recent work2 has illustrated the influence of dispersed precipitates, cold work, and combinations of these on brittle fracture associated with plastic strains of the microyield order of magnitude. This has been rationalized in terms of the contribution of internal microstress fields to the stress magnifications at points of obstruction of active slip bands in a stressed metal. The presently reported studies on an Al-4-1/2 pct Mg alloy amplify this subject by introducing thermal-mechanical conditions which can produce strain aging. As amply demonstrated by Phillips, Swain, and borall, the tensile behavior of aluminum alloys containing magnesium reveals a yield-point elongation and successive discontinuous elongations over the whole plastic-strain range. This gives a stepped or serrated character to the stress-strain diagram. The yield-point elongation is not temperature-dependent. On the other hand, the discontinuous extensions in the plastic-strain range are gradually suppressed by lower temperatures of testing. At -78"C, only the yield-point elongation persists. The yield-point elongation itself can be suppressed by rapid quenching to room temperature from an annealing temperature of 500°C. The argument presented by Phillips et 1. and Cottre114 is that the yield-point elongation derives from surmounting at a critical stress the restraint to slip by grain boundaries which have been reinforced by appreciable magnesium-atom segregation. The serrated character of the plastic-strain range is the result of rapid strain aging during the test. This is made possible by the inevitable proximity of magnesium atoms to any dislocation line. In the Al-4-1/2 pct Mg alloy there is one magnesium atom for every twenty-two aluminum atoms and, with a coordination number of 12, a magnesium atom is on the average only about two atomic distances from any given lattice site. Rapid strain aging is therefore a reasonable expectation. A study of mercury embrittlement of the Al-4-1/2 pct Mg alloy (commercial designation 5083) provides an opportunity to relate the interplay of grain size, grain boundary restraint to slip and strain aging in the brittle-fracture process. One lot of commercial sheet material of 0.125-in. thickness was used for all of the tests. The specimen shape and test procedures used were the same as those described in a previous publication.2 A range of grain sizes between 0.0057 and 0.0125 mm was produced by various combinations of plastic prestrain and annealing. The coarsest grain was somewhat elongated in shape. In this case the measurement was taken in the direction transverse to the axis of tension as more meaningful to a brittle-fracture process. In the annealed condition, specimens wetted with mercury failed with zero permanent elongation at stresses slightly below the normal yield point, i.e., stresses at 0.05 to 0.2 pct strain. Since this occurs before yield-point elongation can begin, fracture derives from slip contained within individual grains and restrained from propagating across grain boundaries. No significant difference in wetted fracture stress was observed between specimens which had been water-quenched or furnace-cooled from the annealing temperatures. This signifies that the degree of magnesium segregation at grain boundaries

Jan 1, 1964

By E. S. Machlin, C. W. Chen

RECENTLY Chen and Machlin' proposed a mechanism for intercrystalline cracking in metals during high-temperature stressing. According to this mechanism the formation of voids at grain boundaries is primarily responsible for intercrystalline cracking. A model has been provided to demonstrate the initiation of submicroscopic cracks and voids at jogs of grain boundaries where grain boundary sliding is inhibited. Submicroscopic cracks can grow as a result of either further grain boundary sliding or localized tensile stress. In the former case, it is believed that grain boundary sliding produces more voids which link together to produce macroscopic cracks and which may lead to the eventual failure of the metal. Concurrently, Nield and Quarrell2 reported their interesting results obtained in a systematic investigation of intercrystalline cracking in two aluminum alloys under various creep conditions. Their observations are generally in agreement with the above mechanism, especially the conclusion that voids are formed as a direct consequence of grain boundary sliding. On the basis of the volume fraction occupied by voids, they could find no confirmation of vacancy condensation as the mechanism for void growth. Furthermore, an important feature of intercrystalline cracking has been likewise pointed out by them that the majority of cracks are not formed at grain corners. Nield and Quarrel1 meanwhile emphasized the importance of grain boundary migration by designating it as one of two main processes operative in the mechanism of crack formation during creep. The designation seems unnecessary because overemphasis of grain boundary migration only complicates the mechanism of void formation. There are two reasons which militate against the inclusion of grain boundary migration in the mechanism of void formation. In the first place, the formation of intercrystalline voids per se does not require grain boundary migration. Secondly, in addition to grain boundary migration, void formation may be minimized by any other process such as plastic deformation which is capable of relieving localized stresses around jogs. Chang and Grant3 have shown that intercrystalline cracking is prevented in aluminum by lattice bending in high stress regions. The influential role played by grain boundary migration in creep has also been realized in our studies. In several cases, voids did not appear at all in grain boundaries whenever grain boundary migration took place prominently during creep. The purpose of this paper is to discuss the effect of grain boundary migration on the formation of intercrystalline voids and other related rupture properties. Experimental evidence indicates that in an unstable grain structure void formation may be markedly modified by structural changes. The abnormally high ductility and the transcrystalline mode of fracture observed in copper at elevated temperatures reveals that effective grain boundary migration minimizes void formation. Grain boundary migration has long been known as a common high-temperature phenomenon in metals with the grain configuration in a nonequilibrium state. Its effect on ductility and mode of fracture has not been recognized until recently.4 Since the introduction of the term "Equicohesive Temperature," the conventional belief has been that below the equicohesive temperature the ductility of metals increases moderately with decreasing temperature until the ductile-to-brittle transition temperature of fracture is reached. In this temperature range, metals fracture in a ductile mode with a "cup-and-cone" fracture surface. Above the equicohesive temperature, ductility remains invariably low and decreases slightly with increasing temperature. The mode of fracture is strictly brittle and intercrystalline. Typical examples have been shown for a brass and copper in Fig. 1 and 9, respectively, in the paper by Greenwood, Miller, and Suiter.5 The observation in the upper temperature range that ductility of metals is low and insensitive to temperature may be explained in terms of void formation at grain boundaries. It is also in accordance with the findings of McLean6 and Fazen, Sherby, and Dorn7

Jan 1, 1961

By Fred H. Poettmann

An equation has been derived for use in calculating the sandface pressure of flowing gas wells in which the variation of the compressibility factor of the gas with pressure is taken into consideration. This variation due to compressibility has been put into both graphical and tabular form. Comparison of calculated results with field measured results were made on 20 dry gas wells from a given field and 11 distillate wells from different fields. The agreement between calculated and observed results is good. In addition to their use in the calculation of sandface pressures of flowing wells, the factors can be used in the direct calculation of the static bottom hole pressure of gas wells, the capacity of gas transmission lines, and in the calculation of the theoretical isothermal horsepower necessary to compress a natural gas. Examples demonstrating the use of various equations are given. INTRODUCTION Many of the equations used to calculate the sandface pressures of flowing gas wells from well head data do not take into consideration rigorously the deviation of the natural gas from ideal gas behavior. For low pressure wells this error is not serious. For high pressure wells flowing at high rates this error call be serious. Sandface pressures are used to determine the producing capacity of gas wells as shown in the U. S. Bureau of Mines' Monograph 7.' An equation has been derived for use in calculating the sandface pressure of flowing gas wells in which the variation of compressibility of the gas with pressure is taken into consideration. This variation due to compressibility has been put into both graphical and tabular form. In addition to their use in the calculation of sandface pressures, these "factors" can also be used to give a direct solution in the calculation of the static bottom hole pressure of gas wells, the calculation of the capacity of gas transmission lines, and the calculation of the isothermal horsepower necessary to compress a gas. CALCULATION OF SANDFACE PRESSURE OF FLOWING GAS WELLS The starting point in the derivation of any specific flow equation is an energy balance on the fluid flowing between any two points in the system. This energy balance based on one pound of flowing fluid is expressed by the well known general flow equation: jp~- vdP+ Ah + -— + wt + r. = o . . (i) P, 2 g. where V = specific volume of flowing fluid in cu ft/lb P = pressure in Ib per sq ft ah = difference in height above datum plane in, ft ------ = kinetic energy change of one lb of flowing fluid, v equals velocity in ft/sec and g. = 32.174 Ws = Work done by fluid while in flow (similar to shaft work in driving a turbine) W, = Znergy losses due to the irreversibilities of the flowing fluid This equation contains no limiting assumptions and can be made the basis of any fluid flow relationship by limitations or substitutions. The energy units are usually expressed in term* of foot-pounds of energy per pound of mass. Considering the vertical flow of gas through tubing or annu-lus, the external work W, done by the gas is zero. Let us also assume that the value of the kinetic energy function is small and can be neglected. Equation (1) is reduced to: /P: VdP + Ah+Wt = 0 .......(2) Pi The energy loszes PI can be expressed in terms of the well known Fanning equatinn: 4/Lf: rf = —— ............(3) 2g,D where L = the integrated average velocity in ft/sec of vertical flow in a uniform flow string -~h = L v = the integrated average velocity in ft/sec / = dimensionless correlating function

Jan 1, 1951

By Fred H. Poettmann

An equation has been derived for use in calculating the sandface pressure of flowing gas wells in which the variation of the compressibility factor of the gas with pressure is taken into consideration. This variation due to compressibility has been put into both graphical and tabular form. Comparison of calculated results with field measured results were made on 20 dry gas wells from a given field and 11 distillate wells from different fields. The agreement between calculated and observed results is good. In addition to their use in the calculation of sandface pressures of flowing wells, the factors can be used in the direct calculation of the static bottom hole pressure of gas wells, the capacity of gas transmission lines, and in the calculation of the theoretical isothermal horsepower necessary to compress a natural gas. Examples demonstrating the use of various equations are given. INTRODUCTION Many of the equations used to calculate the sandface pressures of flowing gas wells from well head data do not take into consideration rigorously the deviation of the natural gas from ideal gas behavior. For low pressure wells this error is not serious. For high pressure wells flowing at high rates this error call be serious. Sandface pressures are used to determine the producing capacity of gas wells as shown in the U. S. Bureau of Mines' Monograph 7.' An equation has been derived for use in calculating the sandface pressure of flowing gas wells in which the variation of compressibility of the gas with pressure is taken into consideration. This variation due to compressibility has been put into both graphical and tabular form. In addition to their use in the calculation of sandface pressures, these "factors" can also be used to give a direct solution in the calculation of the static bottom hole pressure of gas wells, the calculation of the capacity of gas transmission lines, and the calculation of the isothermal horsepower necessary to compress a gas. CALCULATION OF SANDFACE PRESSURE OF FLOWING GAS WELLS The starting point in the derivation of any specific flow equation is an energy balance on the fluid flowing between any two points in the system. This energy balance based on one pound of flowing fluid is expressed by the well known general flow equation: jp~- vdP+ Ah + -— + wt + r. = o . . (i) P, 2 g. where V = specific volume of flowing fluid in cu ft/lb P = pressure in Ib per sq ft ah = difference in height above datum plane in, ft ------ = kinetic energy change of one lb of flowing fluid, v equals velocity in ft/sec and g. = 32.174 Ws = Work done by fluid while in flow (similar to shaft work in driving a turbine) W, = Znergy losses due to the irreversibilities of the flowing fluid This equation contains no limiting assumptions and can be made the basis of any fluid flow relationship by limitations or substitutions. The energy units are usually expressed in term* of foot-pounds of energy per pound of mass. Considering the vertical flow of gas through tubing or annu-lus, the external work W, done by the gas is zero. Let us also assume that the value of the kinetic energy function is small and can be neglected. Equation (1) is reduced to: /P: VdP + Ah+Wt = 0 .......(2) Pi The energy loszes PI can be expressed in terms of the well known Fanning equatinn: 4/Lf: rf = —— ............(3) 2g,D where L = the integrated average velocity in ft/sec of vertical flow in a uniform flow string -~h = L v = the integrated average velocity in ft/sec / = dimensionless correlating function

Jan 1, 1951

By E. Miller, N. Chou, K. Komarek

The solid solubilities of the compounds PbS, PbSe, and PbTe have been determined by equilibrating single crystals with two phase alloys. The solubilities were determined by measuring the Hall coefficient of the equilibrated samples. The observed maximum deviations from stoichiometry on the lead rich and the chalcogen-rich side of the phase fields were: at. pct excess Pb at pct excess S at 850°C; PbSe--4 x 10'" at. pct excess Pb at at. pct excess at 800°C; and PbTe-2.8 x 10 4 at. pct excess Pb at 780°C, 10"1 at. pct excess Te at 720°C. The previously developed theory of retrograde solubility has been expanded to include intrinsic disorder in the compounds, and the experimental results are in agreement with the theory. THE solid solubility limits of intermetallic semiconducting compounds are of considerable interest, since many of these compounds have potential use as thermoelements in thermoelectric devices, as infrared detectors, or in semiconductor devices. The solidus curves of a significant number of these compounds have been reported to be retrograde with the maximum solid solubility occurring between the eu-tectic and the maximum melting temperatures. In a previous paper,' equations were derived which permit calculation of the maximum solid solubility of an AB compound exhibiting retrograde solubility. Derivation of these equations was based on the assumption that the compound is a nondegenerate semiconductor at the temperature under consideration, that the defects were unassociated, and that the vibrational contribution to the entropy of the compound could be evaluated by the Einstein approximation. It was shown that at the temperature of the maximum solid solubility The entropy term ASS is composed of the entropy of formation of the solid stoichiometric compound (AS$)*, and the entropy due to deviations from stoichiometric composition which is the sum of a con-figurational term (aS*), a vibrational term (AS,), and an electronic term (AS,). In the previous paper the positional entropy was evaluated considering the con- N. CHOU, Member AIME, formerly Graduate Student, New York University, N. Y., is with IBM, Yorktown Heights, N. Y. K. L. KOMAREK, Member AIME, formerly with New York University, is Professor, Department of Inorganic Chemistry, University of Vienna, Vienna, Austria. E. MI LLER, Member AIME, formerly with New York University, is Associate Professor of Mechanical Engineering, California State College at Long Beach, Calif. The paper is based on a thesis submitted by N. CHOU to the Graduate Division, School of Engineering and Science, New York University, in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Materials Science. tributions from defects arising from deviations from stoichiometric composition. The derivation is generalized in the appendix to include intrinsic disorder. The electronic entropy term is also recalculated in the appendix. Of the various semiconducting compounds, the solidus curves of Pbs,' and ~b~e,~ have been most extensively studied. The solidus of PbSe has also been determined, but the temperature of maximum solubility has not been accurately established.5129 The present investigation was undertaken to compare the experimental solid solubility data for these systems with the values predicted theoretically. Since the temperature and composition of the maximum solubility of PbSe has not been accurately determined, experimental emphasis was placed on determining the solidus curves for PbSe in the region of maximum solubility, although experiments were also performed on PbS and PbTe at several temperatures. X-ray diffraction studies were also made on carefully prepared powder samples of PbSe to determine the variation of lattice parameter with composition. EXPERIMENTAL PROCEDURE The procedure of Brebrick and Allgaier was employed for determining the solidus data. Monocrystal-

Jan 1, 1970

By John G. Hall

The pumping operations at the Chief mine have been unique in the respect that for many years the entire flow of water into the mine has been disposed of by pumping into natural underground "caverns" or "sinks." Volumes up to 4000 gpm have been pumped into these caverns with no apparent tendency for the caverns to fill with water, and it has never been proven that any of the water returns to the mine workings. Major pumping operations have been carried on during two periods, the first from 1918 to 1927, and the second from 1942 to the present time. Early Pumping Operations CHIEF MINE The Chief Consolidated Mining Co. was organized in 1909, and for the following nine years ore was developed and mined from the 600 level down to the 1800 level, all above the permanent water table, the 1820 level. The ore occurs as extremely-irregular pipes of replacement ore in limestone, and to date a total of approximately 2 1/2 million tons of ore has been produced. The ore was highly oxidized down to the permanent water level, and some highly-oxidized ore runs occur beneath the permanent water table. In 1916, with extensive mining being done on the 1800 level, 20 ft above the water table, there were indications that a zone of secondary enrichment might exist below the water, and in order to follow the ore runs to greater depth, it was realized that arrangements would have to be Made to pump water to the surface. The original pumping installation was designed to pump 350 gpm in five stages against a 1900-ft head to the surface. By January 1919, the water had been lowered to the 1900 level and an excellent grade of silver-lead ore was being produced from the zone of secondary enrichment. At this point. the flow of water equalled the capacity of the pumps, and pumps of an additional 1000 gpln were installed to supplement the original system. By January 1920, the water had been lowered approximately 100 ft, and the total flow of water was 1000 gpm. The water was being pumped from many different places in the mine where ore was being followed downward, with 27 different pumps operating below the 1800 level, mainly Cameron air-driven sinker pumps. These first sinking pumps were cumbersome and inefficient, and each of the larger pumps weighed approximately three tons and required tremeridous volumes of compressed air. During this period, an exploration drift being extended north and west from the Chief shaft on the 1600 level encountered a large natural cavern some 2000 ft from the shaft. This cavern in places measured 60 ft in diameter. Sometime later the disposing of water into this cavern was suggested, and as an experiment to determine the results of such an operation, an attempt was Made to fill the cavern by pumping in a small volume of water, but the water disappeared as fast as it was discharged into the cavern. It was decided then to try the cavern as a drainage outlet for the entire volume of water being pumped to the surface, and the water apparently did not return to the mine workings. Subsequently water up to 3000 gpm was pumped into this cavern. An attempt was made to determine whether any water was being returned from the cavern to the lower workings, and 100 lb of sodium fluorescein salt was added to the water discharged into the cavern, but no trace of it was found in the drainage water in the mine. Results of another such test will be described later.

Jan 1, 1950

By J. L. Moriarty, Irving Johnson, H. M. Feder

The crystallization of ternary alloys has been studied as a coprecipitation phenomenon by measuring the substitutional carrying of trace metals with CeCd11 precipitating from liquid cadmium. The Doerner-Hoskins logarithmic distribution law has been found to apply in all cases; the coprecipitation coefficients obtained were: 1.48 (La), 1.06 (Th), 0.63 (Pr), 0.24 (Gd), 0.17 (Sm), 0.13 (U), 0.10 (Er), 0.09 (Sr). Li, Na, K, Ba, Y, Sc, and Zr were not measurably carried. It has been observed that coprecipitation was most extensive for those tracers which form with cadmium an intermetallic compmdnd isostructural with CeCd11 or closely related thereto. For these tracers an inverse relation between Solubility in cadmium and coprecipitation coefficient was noted. THE partition of impurities between the liquid and solid phases of a metal or a congruently melting intermetallic compound has been extensively studied by the methods of zone melting.' However, when the composition of a liquid alloy and of the solid phase crystallizing from it differ grossly, partition of an impurity cannot be investigated in this manner. Partition of the impurity under equilibrium conditions may be studied by investigation of the appropriate phase diagram and partition under nonequi-librium conditions may be studied by observation of coring in the solidified alloys.2 The nonequilibrium case is of special interest in metallurgical practice because many complex alloys are cooled too rapidly to permit the uniform diffusion of impurities throughout the crystallizing phase(s). With this in mind, the present study was undertaken to obtain quantitative data on the phenomenon of coprecipitation via substitutional solid solution in certain metallic systems. It will be noticed that rather than adopting metallurgical ter- minology we have chosen to refer to the partition of impurities during alloy crystallization as coprecipitation.3 In the terminology of coprecipitation the impurity is referred to as the tracer and the precipitating phase as the carrier. From aqueous solution, coprecipitation may occur as a result of adsorption on the surface of the carrier, by occlusion of the liquid within the crystals, or by solid solution of the tracer in the carrier. By slow cooling of liquid metallic solutions it is relatively easy to form large, regular crystals of the carrier; the effects of adsorption and occlusion may be minimized and the influence of substitutional solid-solution formation may be made to predominate. Previous studies of the crystallization of some intermetallic compounds4 and of copper, silver, and gold5 have shown that the coprecipitation or partition behavior of trace metals is entirely analogous to that observed in aqueous systems when ionic substitution occurs, e.g., coprecipitation of Ra++ with BaSo4. The precipitation of CeCd11 from liquid cadmium was chosen for study. This choice was influenced by the following consideration. If the relative atomic sizes of the tracer atom and the atom for which it substitutes determine the tendency to coprecipi-tate, then a lattice having a substitutional site of well-defined size would be desirable. In CeCd116 (structure type cubic BaHg117) the cerium atom is completely surrounded by a cage of twelve nearly equidistant cadmium atoms. Thus, atoms much larger or much smaller than cerium should be poorly accommodated in this cage. The choice of a lanthanide central atom allowed the testing of the size hypothesis with a number of lanthanide tracers which differed in size but not in valence or electro-negativity. EXPERIMENTAL The general technique of this investigation was to prepare a slightly undersaturated solution of carrier in liquid cadmium to which a tracer was added. The concentration of tracer was always so low that saturation over the entire temperature range investigated was never exceeded. The solution was cooled in a stepwise fashion. After each step a por-

Jan 1, 1964

By L. Burris, C. H. Stockman, I. G. Dillion

Zone melting is a purification process in which separation of impurities is effected by slowly moving a narrow melted zone through a bar of solid material. Equations are presented which 1—predict the concentration profiles after successive passes of the zone and 2—give the limiting distribution of solute along the bar after an infinite number of passes of the molten zone. These equations take into account the effect of normal freezing in the last zone length. These equations were solved numerically using an electronic computer, and the concentration profiles are presented as a function of the solute distribution coefficient and the ratio of molten zone to total bar length. ZONE MELTING is a purification process in which separation of impurities is effected by slowly moving a narrow melted zone through a bar of solid material. Separation occurs if the solid freezing from a solution at any moment has a composition different from that of the liquid from which it is freezing. Solutes which lower the freezing point will tend to remain in solution and will be moved in the direction of zone travel. Solutes which raise the freezing point will preferentially freeze out, and they will be moved in a direction opposite to that of zone travel. After repeated unidirectional passage of a molten zone through a bar of material, impurities will be concentrated in the ends of the bar. A limiting concentration gradient will be reached in the bar after a large number of passes of the molten zone. Both the magnitude of this limiting concentration gradient and the rate of movement of solute in any pass depend on the zone length and on the equilibrium concentrations of solute in the liquid and solid phases. The basic theory of zone melting was thoroughly covered by Pfann.1 The purpose of this paper is to extend the theory of zone melting by presenting 1—equations which relate the concentration profiles in successive passes of the zone and 2—equations which give the limiting distribution of solute along the bar after an infinite number of passes of the molten zone. These equations take into account the effect of normal freezing* in the last zone length concentration of a given solute in the solid to the concentration of the same solute in the liquid. Analytical methods for determining the concentration profiles after any pass have been developed by Lord2 and Reiss. The method developed by Lord does not take into account normal freezing in the last zone and also becomes very complicated when the number of passes becomes large. The method developed by Reiss enables the direct calculation of the complete concentration profiles (including that in the last zone) after any pass for values of k close to unity (0.9 < k < 1.1). Both methods are useful and provide valuable insight into the zone melting process. The equations given in this article are based on the same assumptions used by Pfann. These are: 1—uniform composition in the liquid, 2—negligible diffusion in the solid, 3—constant distribution coefficient, 4—constant length of the liquid zone (except at the ends), and 5—solubility of the solute in the melt is not exceeded at any point. Equations for Concentrations in Successive Passes For a melted zone of length 1 at the left-hand end of a bar of uniform unit cross-section, the amount of solute in the liquid is J0Cn-1(x) dx [1] where the subscript is the number of the pass. Cn-1(x) represents the concentration (expressed as weight per unit volume) obtained in the previous

Jan 1, 1956

By M. R. J. Wyllie, H. W. Patnode

In the quantitative interpretation of electric log data it is essential to know the formation factor, defined as the ratio of the resistivity of the formation 100 per cent saturated with brine to the resistivity of the brine.&apos; This concept implies that the formation factor for any given rock is constant irrespective of the resistivity of the brine or solution with which it is saturated. It is apparent that this can be true only when the rock matrix is wholly non-conductive. As reservoir rocks may contain clay or shale having an appreciable conductance one may expect these constituents to influence the determination of true formation factor. The formation factor may be theoretically derived from log data by dividing the resistivity of the mud filtrate invaded formation by the resistivity of the mud filtrate. Such a method of calculntion has implicit in it the assumption that the fluids in the invaded formation are wholly replaced by mud filtrate. that the solid matrix is non-conducting, and that the formation factor of the invaded zone is identical with the formation factor of the uninvaded zone. Since the validity of the first two basic assumptions regarding the insulating character of porous matrices and the displacement of the interstitial fluids in porous media do not appear to be fully established, the formation factors of a number of cores and clay slurries have been measured under conditions suitable for verifying the assumptions and the extent of any errors inherent in them. CORE STUDIES Apparatus and Procedure Resistivity measurements were made on short cylindrical permeable cores which had been previously mounted in lucite.2 The cores were clamped by means of insulated bolts between brass plates about 5 cm square and 1 cm in thickness. The plates were provided with a ring seal to prevent leakage between them and the core mount, and had a centrally located orifice so that fluids could be flowed through the mounted core. The brass mounting plates were used as electrodes and any current flow between them necessarily passed only through the constituents of the intervening core. Fluid saturation of the core, and fluid flow through the core, were effected by a centrifugal pump connected to the outflow end of the core, thus providing a differential pressure across the core of approximately one atmosphere. The procedure used was first to evacuate each core by means of the pump and then to saturate the core with an aqueous sodium chloride solution of known resistivity. The resistance of the core was measured and flow of the solution through the core continued for sufficient time to establish a constant resistance reading. After equilibrium had been reached with one solution, another sodium chloride solution of a different resistivity was flowed through the core and the flow continued until a constant resistance reading was again obtained. This procedure was continued until resistance measurements of the core had been made with a series of solutions within a salinity range from 1.0N NaCl to distilled water. The data obtained enabled the apparent formation factor of given porous rocks when saturated with fluids of different resistivity to be computed and also provided information regarding the ability of an invading fluid to displace interstitial water from a core. Results Berea Sandstone — The first core examined was a specimen of Berea sandstone." This core had been cut essentially parallel to the bedding for permeability studies. The core was clamped in the core holder and saturated with 1.0N NaCl solution and the resistance determined. It was subsequently flushed and the resistance measured when saturated in turn with 0.1N, 0.025N, 0.01N, .005N sodium chloride solutions and finally with distilled water. In flushing the core with each

Jan 1, 1950

By J. O. Brittain, E. Lautenschlager

Yield point measurements were made on polycrys-talline specimens of a iron deformed various amounts at 770, 1950, or 273°K, aged, and redeformed at 77: 195°, or 273°K. When the prestrain temperature is lower than that subsequent to aging, the yield point return is suppressed relative to specimens deformed at the same temperature both before and after aging. The degree of suppression is a function of both the magnitude of prestrain at the lower temperature and the difference between the deformation temperatures before and after aging. It is proposed that the suppression is caused by the formation of dislocation barriers during the initial low-temperature deformation. The aging kinetics of polycrystalline a-iron wires deformed at 77o or 303oK were studied at 44 " and 64 "C using internal friction measurements. No change was noted as a result of the difference in prestrain temperature. IN an earlier investigation1 it was observed that the return of the yield point was suppressed when a poly-crystalline specimen of a iron was plastically deformed at 77oK, aged at room temperature, and subsequently tested at room temperature. In order to investigate this observation, a research program was undertaken to measure the effect of prestrain temperature on the strain aging of a iron.&apos; According to Thomas and Leak,3 return of the yield point and internal friction are similar measures of strain aging since both are dependent on solute atoms leaving the solid solution. The decrease in internal friction as a function of aging time is indicative of the removal of solute atoms from the solution, while the return of the yield point as a function of aging time gives the rate of arrival of solute atoms at the dislocations. According to Cottrell and Bilby,4 the solute atoms create an "atmosphere" which locks the dislocations and causes a return of the yield point upon restraining. In this experiment, both methods of measuring strain aging were employed in order to have two independent measurements of the strain aging phenomenon. In a previous study,5 the strain aging kinetics of a similar lot of a iron were measured for specimens which were prestrained and tested at the same temperature over the range of 77" to 273°K. Although some differences were found in the correlation of the two methods employed for measuring strain aging, the kinetics of strain aging were unaffected by the change of prestrain and test temperature. In another study,6 dislocation line tension and the density of the solute atmosphere were found to contribute to the return of the yield point when ingot iron specimens were prestrained at 303°K and tested at a lower temperature, namely 77°K. This paper is concerned with changes in strain aging kinetics produced by pre-straining at 77" to 303°K and testing at either the prestrain temperature or at higher temperatures up to 273°K. EXPERIMENTAL PROCEDURE The polycrystalline a iron specimens used in both the tensile and internal friction tests were 0.033-in. diam ingot iron wires, wet hydrogen treated and re-carburized. A more detailed procedure for processing a 3/16 in. ingot iron rod (nominal wt. pet composition: Fe - 0.013; C - 0.0125; Mn - 0.023; Si -0.005; P - 0.0023; S - 0.11; Cu - 0.004 N) to 0.033 in. diam specimens containing approximately 0.020 wt pet C and no detectable nitrogen with a grain diameter of 0.033 to 0.040 mm is given elsewhere.5 Immediately prior to testing, specimens were solution treated for 15 min at 700o ± 2°C in a static vacuum of less than 50 ?. of Hg and rapidly quenched into an ice-water bath while still enclosed in a 3 mm quartz vial. The internal friction specimens had 8-in. gage lengths and the tensile specimens had 2-in. gage lengths. The wires used as tensile specimens were first run as internal

Jan 1, 1962

By R. A. Morse, P. M. Bridges, L. E. Wilsey, Howard N. Hall, P. L. Terwilliger

Theoretical and experimental investigations of a constant pressure gravity drainage system are reported. Experimental data are presented to show that recovery to gas breakthrough by gravity drainage is inversely proportional to rate. The gravity drainage reference rate, which is numerically equal to the so-called "maximum theoretical rate of gravity drainage" is shown to have no particular significance from a recovery standpoint. Before this rate can be used as a basis of comparison for recoveries, it is necessary that the relative permeability and capillary pressure characteristics and displacing fluid viscosities be identical for the systems compared. A method is presented by which accurate prediction of the performance of a gravity drainage system can be made. Close agreement between experimental and calculated drainage performance shows that steady state relative permeability and static capillary pressure data can be used to describe fluid displacement behavior. The very wide range of liquid recoveries before gas breakthrough which result from production rate variation alone demonstrates the importance of this factor in planning depletion of a gravity drainage reservoir. Calculated results are presented which show that little additional recovery can be expected from a high pressure gravity drainage system between the times of gas breakthrough and attainment of such high gas/liquid ratios as to make further pressure maintenance impractical. INTRODUCTION It bas been recognized for some time that gravity forces play an important part in the recovery of oil from some types of reservoirs. Field experience has shown that under certain conditions, gravity drainage can result in very high oil recov- eries. Qualitative reasoning has led most engineers to the general conclusions that: (1) Where gravity drainage is important, the reservoir pressure should be maintained by gas injection at the crest of the structure to prevent shrinkage of the oil in place and to keep a low viscosity so the oil can drain at the fastest possible rate. (2) Recovery by gravity drainage is rate sensitive. A survey of the literature indicates that while considerable work has been done on the effects of gravity in oil production problems, no satisfactory method of calculating the performance of gravity drainage reservoirs has been reported.1,2,3,4,5 In the absence of any proven method of calculating reservoir performance, the level at which pressure should be maintained and the rate of production for most efficient operation has been open to debate. The limits and general character of the recovery us rate relationships are apparent from the following: As rates of production approach zero, oil will drain by virtue of its own weight fast enough that the saturation distribution is always approximately the same as the original static capillary distribution. That is, at the gas-oil contact the liquid saturation goes from 100 per cent to a value of perhaps 10 to 30 per cent in a very short distance, so that as production continues. a relatively sharp gas-oil contact moves down-structure, leaving a very low residual oil saturation in the originally oil saturated zone. At very high rates of production, the oil that drains by virtue of its own weight is negligible and recovery will approach that of a horizontal gas drive, leaving a much higher residual oil saturation. At intermediate rates of production. recovery should be between these extreme*. Despite the lack of quantitative methods, tentative agreement had been reached on the question of optimum producing rate. It has been generally concluded that rather sharp decrease in recovery would occur at rates of production above the "maximum rate of gravity drainage" and hence this rate should not be exceeded. The "maximum rate of gravity drainage" is defined as the rate of production from a 100 per cent

Jan 1, 1951

By L. L. Melton, W. T. Malone

Fluid mechanics research conducted with non-Newtonian fluid systems now permits prediction of the behavior of these fluid systems in both laminar and turbulent modes of flow through circular pipes. Present work concerns non-Newtonian fluid systems currently used in the hydraulic fracturing process. During fracturing treatments, an unsteady-state condition may frequently be encountered arising from. the reaction rate of a chemical additive. This condition must be evaluated in order to predict the actual behavior of a particular fluid during field application. Design and operation of the apparatus used to determine fluid-flow behavior permit obtaining data under such non-equilibrium conditions. This paper shows methods used to obtain rheology measurements, develop hydraulic relationships and evaluate chemical reactions producing unsteady-state conditions. Engineering application of this research is illustrated by employing measured rheological values and developed hydraulic retationships to produce frictional pressure loss (psi/100 ft) us flow rate (bbl/min) charts of common tubing and casing sizes for an example fracturing fluid. How these charts and chemical reaction rate information are then combined to predict actual turbulent hydraulic behavior during unsteady-state field conditions is also discussed INTRODUCTION Many fluids used in hydraulic fracturing contain chemical additives which impart non-Newtonian fluid properties that may drastically alter their hydraulic behavior. Equally drastic alteration in wellhead pressure, injection rate and hydraulic horsepower requirement may result from these fluid properties. Prior research conducted to relate non-Newtonian fluid properties with hydraulic behavior has not as yet produced a universal relationship, particularly for the turbulent flow region. which of the many Possible non-Newtonian fluid properties is responsible has not been conclusively established. A systematic description, suggested by Metzner, of the many possible non-Newtonian fluid properties exhibited by real-fluid behavior, and a current discussion of theoretical and applied aspects of non-Newtonian fluid technoloy can be found in Handbook of Fluid Dynamics. f Little or no research has previously been attempted with fluids exhibiting time - dependent properties. Addition of chemicals during a fracturing treatment is often accomplished by continuously mixing and displacing the fluid. This produces a time-dependent effect or unsteady-state condition while the fluid is progressing through surface and wellbore conductors. This condition is due to solution or chemical reaction of the additive. Considerable departure from conventional methods of obtaining and interpreting data was found necessary to consider these conditions. Therefore methods were developed to obtain hydraulic behavior information under these complex, unsteady-state conditions. Relationships presented in this paper to predict hydraulic behavior of non-Newtonian fluids in circular pipes were obtained by constructing and operating a small pipeline apparatus in the manner of a pilot-plant study. These relationships are suggested as scaleup equations and are not proposed as fundamental rheological parameters. While perhaps deficient from a fundamental research viewpoint, a pilot-plant study does permit the determination and convenient evaluation of variables pertinent to a process. A pilot-plant study can result in a valid engineering application procedure even when fundamental relationships are still undefined. An excellent series of articles by Bowen has appeared in the chemical engineering literature.2-8 These give a thorough review of proposed hydraulic relationships and their limitations for non-Newtonian fluid behavior in circular pipes. A graphical method is presented to scale up data for a fluid exhibiting an anomalous hydrau1ic behavior in the turbulent flow region. Considerable assistance has been obtained from these articles to interpret the anomalous behavior noted during this investigation. These articles also provided assurance that a pilot plant is a practical approach to evaluate the hydraulic behavior of non-Newtonian fracturing fluids. _____________

Jan 1, 1965

By K. H. Coats

This paper presents the development and solution of a mathematical model for aquifer water movement about bottom-water-drive reservoirs. Pressure gradients in the vertical direction due to router flow are taken into account. A vertical permeability equal to a fraction of the horizontal permeability is also included in the model. The solution is given in the form of a dimensionless pressure-drop quantity tabulated as a function of dimensionless time. This quantity can be used in given equations to compute reservoir pressure from a known water-influx rate, to predict water-influx rate (or cumulative amount) from a reservoir-pressure schedule or to predict gas reservoir pressure and pore-volume performance from a given gas-in-place schedule. The model is applied in example problems to gas-storage reservoirs, and the difference between reservoir performances predicted by the thick sand model of this paper and the horizontal, radial-flow model is shown to be appreciable. INTRODUCTION The calculation of aquifer water movement into or out of oil and gas reservoirs situated on aquifers is important in pressure maintenance studies, material-balance and well-flooding calculations. In gas storage operations, a knowledge of the water movemenr is especially important in predicting pressure and pore-volume behavior. Throughout this paper the term "pore volume" denotes volume occupied by the reservoir fluid, while the term "flow model" refers to the idealized or mathematical representation of water flow in the reservoir-aquifer system. The prediction of water movement requires selection of a flow model for the reservoir-aquifer system. A physically reasonable flow model treated in detail to date is the radial-flow model considered by van Everdingen and Hurst.1 In many cases the reservoir is situated on top of the aquifer with a continuous horizontal interface between reservoir fluid and aquifer water and with a significant depth of aquifer underlying the reservoir. In these cases, bottom-water drive will occur, and a three-dimensional model accounting for the pressure gradient and water flow in the vertical direction should be employed. This paper treats such a model in detail — from the description of the model through formulation of the governing partial differential equation to solution of the equation and preparation of tables giving dimensionless pressure drop as a function of dimensionless time. The model rigorously accounts for the practical case of a vertical permeability equal to some fraction of the horizontal permeability. The pressure-drop values can be used in given equations to predict reservoir pressure from a known water-influx rate or to predict water-influx rate (or cumulative amount) when the reservoir pressure is known. The inclusion of gravity in this analysis is actually trivial since gravity has virtually no efFect on the flow of a homogeneous, slightly compressible fluid in a fixed-boundary system subject to the boundary conditions imposed in this study. Thus, if the acceleration of gravity is set equal to zero in the following equations, the final result is unchanged. The pressure distribution is altered by inclusion of gravity in the analysis, but only by the time-constant hydrostatic head. The equations developed are applied in an example case study to predict the pressure and pore-volume behavior of a gas storage reservoir. The prediction of reservoir performance based on the bottom-water-drive model is shown to differ significantly from that based on van Everdingen and Hurst&apos;s horizontal-flow model. DESCRIPTION OF FLOW MODEL The edge-water-drive flow model treated by van Everdingen and Hurst1 is shown in Fig. la. The aquifer thickness b is small in relation to reservoir radius rb. water invades or recedes from the field at the latter&apos;s edges, and only horizontal radial flow is considered as shown in Fig. 1b. The bottom-water-drive reservoir-aquifer system treated herein is sketched in Fig. 2a and 2b. Here the aquifer

By B. W. Levinger

THE system Ti-CI- has been studied by several investigators." &apos; Though titanium and chromium are completely miscible at high temperatures, an intermediate phase of the approximate atomic proportions TiCr, has been identified at various compositions between 60 and 68 wt pct Cr. The structure of this phase was established as the cubic C15 (MgCu2) type with 24 atoms per unit cell and a lattice parameter 6.929 kX.2 During a study of the ternary system, Ti-Mo-Cr,5 the existence of a high temperature modification of the TiCr, phase was discovered. In a series of alloys annealed at 1300°C and containing from 60 to 70 wt pct Cr and 5 pct Mo, a new phase was present which did not correspond to either ß solid solution or TiCr,. By trial, the new phase was identified to be of the hexagonal C14 (MgZn2) structure with 12 atoms per unit cell. Further work showed that the hexagonal structure was. in fact, also present in binary Ti-Cr alloys quenched from high temperatures. To make the identification more certain, the intensities of the various reflections were calculated assuming a composition TiCr2 and C14 structure. Table I lists reflection planes, interplanar distances, calculated intensities, and a comparison of these values with the lines of a binary Ti-Cr alloy (66 pct Cr) annealed for 20 min at 1385°C and water quenched. Since high angle lines were more diffuse, the correspondence between visually observed and calculated intensities was not as good at these angles. Lattice parameters were computed to be a = 4.922 kX, c = 7.945, c/a :- 1.614. The corresponding interatomic distances are: titanium to titanium distance, 3.00:, kx, chromium to chromium distance, 2.43, kX, titanium to chromium distance, 2.88,, kX. These compare to 3.00, 2.45, and 2.87 kX in the cubic phase, respectively.&apos; The atomic positions in the C14 structure are as follows: Titanium 1/3, 2/3, Z; 1/3, 2/3, 1/2: 2/3, 1/3, i; 2/3, 1/3, (1/2-z); Chromium 0, 0, 0: 0, 0, 1/2; 2x, 2, 1/4; 232, x, %; X, 22, 1/4; X, 2x, 3/4; x. x, 1/4; x, x, 3/4; The values of x and z were calculated for the alloy mentioned and found as follows: z — 0.061, x = 0.165. The values assumed in the intensity calculations were those for an ideal close-packed lattice or z = 1/16, x = 1/6. Work in the Ti-Mo-Cr system indicates that y2 the hexagonal modification of TiCr2, is the stable phase above 1300°C, while y1 with the cubic structure is stable below 1000°C. X-ray examination of a number of alloys annealed at intermediate temperatures showed that the hexagonal and cubic phases can coexist. Alloys quenched from 1000°C and annealed at 1100" and 1200°C showed small amounts of y2 at the higher temperature. A considerable amount of y, seemed present, however, in cast alloys annealed at 1100°C. It further appears that there is a small but definite solubility range for both phases in the ternary system. Acknowledgment The results of this paper represent work carried out in conjunction with a project sponsored by Watertown Arsenal, Contract No. DA-11-022-ORD-272. References &apos;M. K. McQuillan: Journal Inst. Metals (1951) 79, p. 379. 2 P. Duwez and J. L. Taylor: Trans. ASM (1952) 44, p. 495. W. J. Van Thyne,. H.. D. Kessler, and M. Hansen: Trans. ASM (1952) 44, p. 974. 1 F. B. Cuff, N. J. Grant, and C. F. Floe: Trans. AIME (1952) 194, p. 848; Journal of Metals (August 1952). Unpublished Research, Armour Research Foundation. "H. J. Wallbaum: Ztsch. Krist. 103, p. 391.

Jan 1, 1954

By H. E. Cline, T. H. Alden

A study zvas made of rate-sensitive plastic deformation and superplasticity in alloys of lead and tin. These alloys (lead, PO-19 pct Sn, Sn-38 pct Pb, Sn-2 pct Pb, tin) were extruded, aged at room temperature, or melted and resolidified to produce a variety of ~nicrostrnctures. Flow stress us strain rate data were obtained frorn tension tests for strain rates between 10-6 and 1 min-1. The curve log o .us log £ was S-shaped jor fine-grained materials and approxinzately linear for coarse-grained nzaterials. At constant grain size, the behavior was relatively insensitive to conposition. High rate sensitivity (d log o/d log k) and superplasticity occurred at intermediate strain rates in fine-grained materials and are believed to be associated with a transition with strain rate in the mode of plastic flow. This transition was shifted to higher strain rates with increasing temperataire or decrensing grain size. The high rate rneclzanisrn is presunzed to be slip, while rnetallographic evidence sug-gests that the low rate lnechanism is pain boundary sliding. PEARSON&apos; discovered that extruded eutectic alloys of Pb-Sn and Bi-Sn could be elongated as much as 2000 pct in tension. The ability for stretching decreased with aging time or annealing to increase the grain size. At high strain rates, or if the material was in the cast condition, normal ductility was observed. Following deformation of the extruded Pb-Sn eutectic at low rates, Pearson observed an equiaxed granular structure on an initially flat surface, while rapid straining produced elongated grains. He attributed the large deformation in the extruded alloys to the persistence of extremely small grains to a temperature sufficiently high for inter granular flow. Avery and backofen2 measured the strain rate sensitivity of the flow stress, m = d In old In 6, in the Pb-Sn eutectic and established a correlation between m and the tensile elongation. In the extruded eutectic alloy the strain rate sensitivity was 0.1 to 0.2 at rapid rates and the elongation was small, while at low rates m increased to 0.4 to 0.6 and superplastic elongation was observed. They also noted a grain size effect. In the case of the cast Pb-Sn eutectic the strain rate sensitivity remained at 0.1 over the range of strain rates of their experiment (3 x 10"3 to 2 min-I). These authors suggested that diffusional (Nabarro-Herring) creep was the mechanism by which deformation occurred at low strain rates, while climb-controlled creep occurred at high rates. Two predictions of this model are that, at low strain rates, the strain rate sensitivity approaches 1, i.e., Newtonian viscosity, and that the strain rate at constant stress varies as D/L2, where D is the coefficient of self-diffusion and L the grain size. More recently, a mechanical and microstructural investigation of the dilute superplastic alloy, Sn-5 pct Bi, has been done.3 The results show that the high strain rate sensitivity is associated with a transition with strain rate of the mechanism of plastic flow. At low strain rates, grain boundary sliding and its associated accommodation processes were evident, while at high strain rates quasi-uniform slip was concluded to be the dominant deformation process. The following experiments were performed to extend our knowledge of the phenomena of superplasticity to a larger range of strain rates, compositions, temperatures, and microstructures in the Pb-Sn alloy system. The question of why a cast structure does not exhibit viscous behavior although it appears to have a fine microstructure2 was investigated.

Jan 1, 1968

By C. Wagner

When an alloy solidifies and the composition of the solid differs from the composition of the liquid, atoms of the alloying elements rejected at the solid-liquid interface have to diffuse toward the bulk liquid. Diffusion may be supported by convection. Theoretical calculations have been made for a liquid without convection, for a liquid involving natural convection, and for solidification of a liquid alloy at the surface of a rotating disk as an example of forced convection. WHEN an alloy solidifies, the equilibrium composition of the originating solid phase differs in general from the composition of the liquid phase. Thus, segregation may take place, i.e., the composition of the solid phase formed at different times differs from the bulk composition. The degree of segregation depends on several factors such as the partition ratio of the alloying elements between the liquid and the solid phase, the rate of solidification, and convection. To clarify, a theoretical analysis has been worked out. The following presuppositions are made: 1—Calculations are confined to the solidification of binary alloys in which the concentration of component A (solvent) is much greater than that of component B (solute). 2—The equilibrium concentration of the alloying element B is assumed to be lower in the solid phase than in the coexisting liquid phase. 3—At the solid-liquid interface, partition equilibrium is assumed. 4—The dependence of the partition ratio of the solute between the liquid and the solid phase on concentration and temperature is disregarded since only small concentrations of the solute are considered. 5—The concentration of the alloying element B in the solid phase is supposed to be lower than the saturation concentration for the formation of a second solid phase. 6—diffusion in the solid phase is disregarded, since diffusion rates in the solid state are comparatively low. When the partition ratio of the solute between the liquid and the solid phase is greater than unity, a part of the solute is rejected at the solid-liquid inter- face and has to diffuse toward the bulk liquid. Thus there is necessarily a concentration difference in the liquid between the solid-liquid interface and the bulk liquid as has been pointed out by Hayes and Chipman,&apos; Rutter and Chalmers,&apos; Tiller, Jackson, Rutter, and Chalmers, and others (see Fig. 1). In view of the concentration gradient in the liquid phase, the ratio of the concentration of the solute in the solid phase to that in the bulk liquid differs from the equilibrium ratio, although at the solid-liquid interface, partition equilibrium is assumed. Under certain limiting conditions, the solid may even have virtually the composition of the bulk liquid as is shown below. General Equations Let u be the linear crystallization velocity in cm per sec, D the diffusion coefficient of component B in the liquid phase in cm2 sec-&apos;, c&apos;,, and c", the concentrations of component B in the liquid and the solid phase, respectively, in mol per cm", and s the distance from the solid-liquid interface toward the bulk liquid at a given time t. The partition ratio is assumed to be equal to a partition constant K. Thus. The difference in the specific volume between the liquid and the solid phase may be disregarded. In view of the discontinuity of the concentration of the solute at the interface according to Eq. 1, u(c&apos;, — c",), , , = tic&apos;, (S = 0)(1 — K) gram-atoms B per unit area per unit time are rejected at the interface and have to diffuse toward the bulk liquid. From Fick&apos;s first law it follows that uc&apos;,,(s = 0) (1 ~K) == - D(dc&apos;,,/?>s). _-„. [2] The concentration profile in the liquid near the interface depends on the interplay of diffusion and convection as is discussed in the following sections. A schematic concentration profile is shown in Fig. 1. The "effective thickness of the diffusion boundary

Jan 1, 1955

By L. E. Wilsey, R. A. Morse, P. L. Terwilliger, Howard N. Hall, P. M. Bridges

Theoretical and experimental investigations of a constant pressure gravity drainage system are reported. Experimental data are presented to show that recovery to gas breakthrough by gravity drainage is inversely proportional to rate. The gravity drainage reference rate, which is numerically equal to the so-called "maximum theoretical rate of gravity drainage" is shown to have no particular significance from a recovery standpoint. Before this rate can be used as a basis of comparison for recoveries, it is necessary that the relative permeability and capillary pressure characteristics and displacing fluid viscosities be identical for the systems compared. A method is presented by which accurate prediction of the performance of a gravity drainage system can be made. Close agreement between experimental and calculated drainage performance shows that steady state relative permeability and static capillary pressure data can be used to describe fluid displacement behavior. The very wide range of liquid recoveries before gas breakthrough which result from production rate variation alone demonstrates the importance of this factor in planning depletion of a gravity drainage reservoir. Calculated results are presented which show that little additional recovery can be expected from a high pressure gravity drainage system between the times of gas breakthrough and attainment of such high gas/liquid ratios as to make further pressure maintenance impractical. INTRODUCTION It bas been recognized for some time that gravity forces play an important part in the recovery of oil from some types of reservoirs. Field experience has shown that under certain conditions, gravity drainage can result in very high oil recov- eries. Qualitative reasoning has led most engineers to the general conclusions that: (1) Where gravity drainage is important, the reservoir pressure should be maintained by gas injection at the crest of the structure to prevent shrinkage of the oil in place and to keep a low viscosity so the oil can drain at the fastest possible rate. (2) Recovery by gravity drainage is rate sensitive. A survey of the literature indicates that while considerable work has been done on the effects of gravity in oil production problems, no satisfactory method of calculating the performance of gravity drainage reservoirs has been reported.1,2,3,4,5 In the absence of any proven method of calculating reservoir performance, the level at which pressure should be maintained and the rate of production for most efficient operation has been open to debate. The limits and general character of the recovery us rate relationships are apparent from the following: As rates of production approach zero, oil will drain by virtue of its own weight fast enough that the saturation distribution is always approximately the same as the original static capillary distribution. That is, at the gas-oil contact the liquid saturation goes from 100 per cent to a value of perhaps 10 to 30 per cent in a very short distance, so that as production continues. a relatively sharp gas-oil contact moves down-structure, leaving a very low residual oil saturation in the originally oil saturated zone. At very high rates of production, the oil that drains by virtue of its own weight is negligible and recovery will approach that of a horizontal gas drive, leaving a much higher residual oil saturation. At intermediate rates of production. recovery should be between these extreme*. Despite the lack of quantitative methods, tentative agreement had been reached on the question of optimum producing rate. It has been generally concluded that rather sharp decrease in recovery would occur at rates of production above the "maximum rate of gravity drainage" and hence this rate should not be exceeded. The "maximum rate of gravity drainage" is defined as the rate of production from a 100 per cent

Jan 1, 1951