Search Documents

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 P. E. Hildebrand, S. W. Mao, C. N. Crain

Hydrologic and economic criteria figure in many obvious ways in water resources development, but they are rarely linked quantitatively, and most of the applications are pertinent only to the case in point. In the preparation of a regional plan for the massive program of groundwater development in the Punjab, situations have arisen where the choice among two or more alternatives was made by analyzing the economic implications of critical hydrologic parameters. Two of these situations which are of more or less general interest are described herein. The first involves basic development policy, i.e., the question of whether groundwater development should be based on the demand of the economy, or the conservation concept of safe yield. The analysis demonstrated that mining of groundwater is clearly indicated for areas where alternative supplies are unavailable. The second problem involves integrating hydrologic data into an economic model to determine optimum well design. From an analysis of all pertinent factors, a simple nomograph was derived for determining well specifications on site during construction. The nomograph can be modified for use in any alluvial terrain. In the Northern Zone of the Indus Plains —the region known as the Punjab- Pakistan is proceeding with integrated development of water and land resources on a scale never before contemplated. The heart of the development plan is a massive program of groundwater development to supplement canal irrigation supplies and to provide drainage relief to the irrigated lands. By the close of 1966 about 2500 high-capacity wells were in operation and furnishing the water supply and drainage requirements for nearly two million irrigated acres. The development plan calls for the construction of a total of nearly 30,000 wells with an aggregate capacity of 120,000 cfs by 1980. This is part of an integrated program eventually affecting over 20 million acres of irrigated lands. [This development program is administered by the West Pakistan Water and Power Development Authority (WAPDA) for which the firm of Tipton and Kalm-bach, Inc. serves as regional planning consultant for the Northern Zone. The material presented herein is derived from planning studies which have been made in connection with the preparation of a regional development plan for the Northern zone.] Pakistan is proceeding with this large-scale development of land and water resources at a time when she is at the threshold of an almost calamitous population growth. Furthermore, she is attempting this development at a time when most of her usable land is already occupied and her food demands exceed current production. Here in the Northern Zone is one of the largest continuous areas of irrigated land in the world, easily irrigated, with a warm, dry climate and ample potential supplies of surface water and groundwater. Except for an abundance of land and labor, other primary resources are lacking or in very limited supply. The agriculture sector so dominates the economy that only by means of an effective development of the agricultural resources can other sectors grow and contribute to a viable economy. Because of inadequate irrigation, farming here has not effectively participated in the twentieth century revolution in agriculture, and only by means of massive development of potential supplies can Pakistan catch up. Capital for any kind of mass development must come from the outside. But the overriding restraint to development of the land and labor resources is the availability of water. ECONOMICS OF DEVELOPMENT The decision-making process in planning regional development becomes increasingly complex as the various sectors of the economy become more highly developed. Where, as in this case, one sector dominates the economy and its development can be considered the key to development of other sectors, the decision-making process is simplified and it is easier to identify alternatives and to trace the prob-

Jan 1, 1970

By John Alden

A CONSIDERABLE PROPORTION of the tonnage of aluminum used today has at some stage during its manufacture, been rolled. I. For one, therefore, was particularly pleased to hear that this Society was to hold a symposium on rolling of nonferrous metals. Metallurgists can well afford to direct their close attention to the metallurgy involved in this important operation. To one who has seen today's high production rolling mills for aluminum it is difficult to realize that the first aluminum rolling mill For sheet was installed a scant 55 yr ago, in 1892. This mill was a three-stand train of 2-high mills, driven by a steam engine and was operated for many years successfully rolling aluminum sheet. Mills for the rolling of aluminum rod were installed in 1899 to meet the demand For manufacture of electrical conductor wires and cables. Bars and rods in straight lengths were also made. While these mills were the first built for rolling aluminum products of these types some rod had previously been rolled on brass and copper mills. Not until many years later, in 1929, however, Were facilities provided for the rolling of structural shapes, such as I-beams and channels. Growth of the aluminum sheet business and the necessity for economical manufacture soon brought the realization of the economies resulting From-the rolling of sheet in long lengths and coiling it. A five-stand continuous mill driven by a steam engine, was in operation late in 1907. This took aluminum slabs coating from

Jan 1, 1948

By R. A. Wattenbarger, H. J. Ramey

A systematic study has been made of the application of the real gas pseudo-pressure m(p) to short-time gas well testing. The m(p) function can be used in real gas flow problems to account for the variation of viscosity and z-factor with pressure. A mathematical model was solved numerically to generate solutions of various real gas flow problems. The effects of turbulence, formation damage and wellbore storage were included in the model. The analysis of simulated well tests showed that the interpretation methods normally used for liquid flow are generally accurate when the m(p) is used. For practical rates without turbulence, the solutions are not rate-sensifive, and the flow capacity.. kh and skin effect can be determined accurately from either a buildup or a drawdown test. However, the kh calculated from a drawdown test can be significantly low when turbulence is present. A case was simulated in which this error wa.s 36 percent. Turbulence does not affect the determination of kh from buildup tests. The proper determination of kh from p and pa buildup and drawdown plots is developed by analyzing their relationship to the m(p) method. A simple equation is given that can be used for long-range gas well performance forecasting. This expression is compared with a method pre-sended by Russell et al." Introduction The formation flow capacity and wellbore damage condition can be determined for liquid-producing tests by means of buildup and drawdown tests. These tests make use of short-time pressure transient data rather than stabilized flow tests that are often used for gas well testing. Tracy' presented a method of gas well testing that was based on ideal gas equations and utilized p2 buildup plots. This method was shown to be good for low pressure wells. MatthewsZ suggested plotting p and using an average slope and average gas properties. This method was more successful on high pressure wells. Al-Hussainy, Ramey and Crawfordh howed that variation in gas properties could be simplified by using the real gas pseu do-pressure m(p). Some cases showed that the use of the m(p) function Original manuscript received in Society of Petroleum Engineers office Aug. 8. 1967. Revised manuscript received June 11, 1968. Paper (SPE 1835) was presented at SPE 42nd Annual Fall Meeting held in Houston Tex., Oct. 1-4. 1967, and at SPE 4th Annual Eastern Regional ~eetin$ held in Pittsburgh. Pa., Nov. 2-3. 1967. 0 Copyright 1968 American Institute of Mining, Metallurgical, and Petroleum Engineers, Inc. 'References given at end of paper. provided an accurate method of interpreting buildup and drawdown tests. Because the work of Al-Hussainy et al. included only a few actual cases, there was a need to explore the use of the m(p) method for a greater variety of flow conditions. The occurrence of turbulent flow around the wellbore often is an important factor in gas well testing. Swift and Kiel,' and Carter et a1.' treated turbulence in gas well testing, but did not consider the variation of viscosity and z-factor with pressure. This paper presents the results of an investigation of the application of the m(p) method to buildup and drawdown testing. Vhe investigation included the effects of turbulence, formation damage and wellbore storage. Gas Flow Equations To formulate the mathematical model, many of the assumptions usually used in well testing theory are applied. The system has radial geometry with a closed outer boundary and is composed of a horizontal porous formation that has uniform and isotropic rock properties and uniform thickness. Allowance is made, however, for a radial region of reduced permeability near the wellbore. This region represents formation damage. The geometry of the system is shown in Fig. 1. Darcy's law does not always apply to gas flow. A more general expression is needed for non-Darcy, or "turbulent" REGION OF DAMAGED PERMEABILITY SEALED BOUNDARIES ClRCULAR WELL BORE Fig. 1—Radial flow model. AUGUST, BT7

Jan 1, 1969

By L. S. Darken

A quadratic formalism is developed lor the representation of the excess free energy, and of the activity coefficients of each component of a ternary system in the vicinity of a single component selected as solvent. This formalism, in contrast to the e formalism, is thermodynamically consistent at .finite concentrations. Numerous specific systems are used as illustrations, with cross checks where possible; these systems are Fe-C-X systems in the vicinity of 1600°C and include practically all those for which reliable data are available. It is shown that this formalism is usually quite adequate to rebresent available data up to a solute concentration of 20 to 30 at. pet. It is hoped that the formalism here presented may prove useful for the treatment and tabulation of data on ternary and multi-component solutions. SINCE, as discussed in a prior paper,' we do not as yet have an adequate basic understanding of the nature of binary metallic solutions, we cannot hope to be in better position for ternary and multicomponent systems. However, the accumulation of thermodynamic data on ternary and multicomponent metallic solutions has led to an urgent need for a rational formalism for the presentation and use of such data. Notable attempts to represent the isothermal isobaric thermodynamics of ternary solutions were made by Benedict et 01.' and by Wohl.3 These attempts were directed toward developing an analytic expression adequate to encompass the entire ternary triangle, primarily for organic systems; they generally involve the power-series approach, which has been shown' to be inappropriate, in general, even for binary metallic systems. On the other hand, interest in the thermodynamics of ternary and multicomponent metallic solutions has tended to focus either on general thermodynamic relations4-10 over wide compositional ranges, or on cases where one or more components have been at rather low concentration. This latter approach has been followed by Wagner,= by chipman,11 and by Alcock and Richardson,12 and more recently in many publications by the same and other authors. It has already been shown' that the thermodynamic behavior of binary metallic systems is, in general, relatively simple only in the terminal regions. Hence, it would appear that any approach, at this time, to a simple formalism for ternary and multicomponent systems must, a fortiori, be limited to the vicinity of the major component or solvent. The present treatment is aimed to be useful for a given solvent with moderate solute concentrations ranging up to 20 or 25 at. pet and in exceptional cases to considerably higher concentrations. In order to achieve even this modest objective, it is necessary that any appropriate formalism have built into it 1) a thermodynamic consistency (at least within the desired degree of approximation), and 2) an approach at infinite dilution to Raoult's law for the solvent and to Henry's law for each solute. Let us first formulate the condition of thermodynamic consistency for a ternary system. Under isothermal isobaric conditions, we may write for the molal value G of any extensive function dG - C1 dN1 + G2 dN2 + G3 dN3 [l] where the G's are the corresponding partial molal quantities and the N's are the atom fractions. Taking component 1 as the solvent, this relation may be rewritten in terms of the independent variables N2 and N3

Jan 1, 1968

By A. H. Geisler, J. H. Keeler

Preferred orientations in rolled and annealed titanium sheets were determined by the Geiger counter spectrometer X-ray diffraction technique. Five annealing textures dependent upon the temperature range of annealing were found, and in order of increasing annealing temperature pendent upon the temperature range are: 1—a deformation like texture, 2—a rotated inorder a-recrystallization temperature texture, 3-a retained u-recrysraIlization texture, on annealing at lower temperatures of the ß-region, 4—a transformation texture based on recrystallized a and predicted by the Burgers' relationship, and 5—a ,ß-cube texture. These results are examined in terms of current theories of recrystallization textures. UMEROUS investigators have described the tex- ture obtained by cold rolling the hexagonal metals, titanium, zirconium, and beryllium, which have c/a ratios less than that of ideal packing, 1.633. The basal planes are rotated out of the rolling plane, about the rolling direction, so that the basal poles are tilted toward the transverse direction as shown schematically in Fig. la. In all instances but one,' it was also reported that the [1010] direction was parallel to the rolling direction (see Fig. lb). Hot rolling has been reported as causing a similar tilt of the basal poles in the transverse direction (see Fig. la) and causing the [1010] direction also to be parallel to the rolling direction as shown schematically in Fig. lb. Annealing after deformation does not appreciably change the tilt of the basal poles in the transverse direction." Beryllium2-7 continues to have the [1010] direction in the rolling direction after annealing, and similar observations for titanium and zirconium' . have been reported for annealing at fairly low temperatures, again as in Fig. lb. At higher annealing temperatures, however, the recrystallized grains of titanium" and zirconium have an orientation such that the [1120] direction is approximately in the rolling direction, although the basal poles are still inclined in the transverse direction. Figs. la and lc show the resulting orientations schematically. This change in orientation has been described as a nominally ±30° rotation of the hexagonal crystallites about the basal poles of the cold rolled texture and is apparent from the results which are summarized in Table I for investigations with the X-ray diffraction technique employing film. The angles y, , and ß are indicated in Fig. 2 which represents the stereographic projection of (1070) poles for the mean orientation of a pole figure. Texture determinations for titanium using the Geiger counter spectrometer have provided similar results except that in some instances additional components of the texture were proposed, as shown by the summary of data in the upper half of Table 11. On the other hand, the spectrometer technique, when applied to zirconium,* has revealed a splitting Recently completed studies of the textures of annealed zirconium", show zirconium to possess textures very similar to those reported here for titanium. Therefore, much of this discussion will include zirconium by virtue of its close similarity to titanium in pref erred orientations. of the intense areas of the pole figure for samples annealed at 600°C. This splitting could be described by a 7" rotation of the tilt axis about the normal to the rolling plane. Such a splitting for the annealed texture relative to the cold rolled texture was not observed in other determinations for either zirconium or titanium using the less sensitive film X-ray methoe and makes the relationship between the two types of texture more complex than the simple rotation about the (0001) pole based on film work. The more precise investigations on zirconium permit the descriptions in the lower part of Table 11, which show that the texture depends quantitatively on the temperature of annealing. When zirconium is annealed at temperatures up to 400°C, the texture is similar to the cold rolled texture, while annealing in the range 500" to 900°C produces a texture which is only approximately described as [11%] in the rolling direction. More precisely described results for zirconium show that the two types of splitting ( 1—about an axis in the rolling plane through an angle given in the second column in Table II and 2—about the normal to the rolling plane through an angle given in the third column of Table 11) depend on annealing temperature. The [1120 is the rolling direction only when the annealing temperature is in the vicinity of 900°C

Jan 1, 1957

By L. F. Mondolfo, B. E. Sundquist

The crystallographic orientation relationships resulting when lead is nucleated from the liquid by Ni, Cu, Ag, and Ge were determined. For each nucleating agent several definite orientatioz relationships were found. These relationships seemed to be controlled by good symmetry relations and low crystallographic disregistry between mating planes. For any given nucleating agent the under colling for nucleation was found fairly constant and independent of the orientation relationship and consequent disregistry. It was also found that, upon re melting and refreezing the Pb, the orientation relationship was changed. These findings prove that crystallographic disregistry is not the controlling factor in heterogeneous nucleation from the liquid. The results of this investigation tend to confirm the theory presented in a preceding paper that heterogeneous nucleation starts with the formation of an adsorbed layer of nucleated metal on the nucleat-ing impurity. Evidence is given that cavities in the nucleating agent act as centers of nucleation. IT has long been known' that solid extraneous particles are active in catalyzing phase transformations that occur in a system, particularly condensation and crystallization. It is well established that these heterogeneities act as catalysts by providing surfaces upon which nuclei of the precipitating phase can form with activation energies smaller than those required for homogeneous nucleation. Numerous investigations have shown that in this process of heterogeneous nucleation: a) the nucleus forms with one, or several, definite crystallographic orientation relationships with the nucleating phase2-4 and b) that there is a small range of undercoolings or super saturations characteristic of the nucleation of a given solid on a given Substrate.5-10 Turnbull and vonnegut11 have developed a theory based on theories developed by Volmer12 and Turn-bull and Fisher1= for heterogeneous nucleation from gases and liquids, that relates the super saturation or undercooling required for nucleation to the dis-registry between the lattices of the nucleus and the nucleating agent. This theory predicts that nucleation should occur with the orientation relationship between the nucleus and nucleating agent that minimizes the disregistry. Further, it predicts that the undercooling or super saturation necessary for nucleation should be a function of the disregistry. Numerous investigations have dealt with the orientation relationships resulting from the condensation of vapors onto crystalline solid substrates2,3 and a few with the nucleation of one phase by a second phase in solidification4,14. Others have dealt with the supersaturation8-10 and undercooling5-7 associated with nucleation in condensation and solidification respectively. However, there is virtually no report that gives both of these factors for the same system. In this investigation a study was made of the undercoolings and orientation relationships resulting when Pb is nucleated from the liquid by Ni, Cu, Ag, and Ge. It was the purpose of this investigation to check the Turnbull-Vonnegut theory, i.e., the importance of crystallographic disregistry between nucleating catalyst and nucleated metal. The results indicate that disregistry is not an important factor in nucleation and that the nucleation process is probably somewhat more complex than current theories suggest. EXPERMENTAL PROCEDURE Small single crystals of nickel, copper, silver, and germanium were prepared from materials of four to five nines purity, and the Pb used was also 99.999+ pet pure. Cu and Ag single crystals were prepared by sealing small chips of Cu or Ag in an evacuated quartz capsule and heating the capsule at 2000°F for 1 hr before cooling. Nickel crystals of 200 diam were also prepared in evacuated quartz capsules, but melting was done by heating the capsules in an oxy-acetylene flame for a few minutes. These spheres were invariably polycrystalline so

Jan 1, 1962

By R. J. Tailleur, M. Rosenberg

Oil well drilling muds prepared with lubricating aids such as oil, graphite or mica will not produce a sufficiently strong protective lubricant film for the bit bearing surfaces under high load conditions encountered in drilling operations. Consequently, the damage to the bearing surfaces resulting from the absence of such a film frequently results in a shortening of the drill bit life. Several different types of extreme pressure lubricating additives for use in muds have been developed. These arlrlitives are capable of providing a protective lubricating film for bearing surfaces subjected to high load pressures. A modified version of the Timken lubricant tester, an instrument commonly used for testing the extreme pressure lubricating properties of oils, is used for measuring and controlling the lubricating qualities of the mud. Comparison tests made in the laboratory with 55/8-in. tricone bits drilling on steel plates immersed in extreme pressure lubricating muds and conventional muds showed that the bearing life of those bits tested in the lrrbricating mud was three to four times greater than obtained with conventional muds. Drilling tests made with extreme pressure lubricating muds in West Texas increased over-all bit life three to fourfold, reduced drilling torque and wear on drill pipe and drill collars. INTRODUCTION Roller bit cones contain roller, ball and journal-type bearings that are subjected to very heavy loading during drilling operations in which the only lubricant available for these bearings is either water or drilling mud. These bearings must be able to operate satisfac- torily in order for the cutting structures on the surface of the cones to be fully utilized. Frequently, the bearings wear to the point that the cones freeze and the bit must be pulled before the cutting structures on the cones become worn. In some instances the bearings continue to turn with severe wear but the cones, which are held on the bit by the bearings, drop off; and this causes an expensive fishing job. In other instances the bearings wear to the point that the cones pinch in causing an undergauge hole and delays to ream the hole. The advantage of increasing bit bearing life is self-evident. Because of inability to control the internal wear of bits while drilling, it was felt that improvement of the lubricating qualities of the circulating mud would be helpful. However, the often-discussed lubricating property attributed to muds has been a rather vague and intangible one. The literature has implied that the oil-emulsion muds are more effective lubricants than non-emulsion muds. Actually, there has been no direct method of measuring the lubricating property of muds which would indicate that one mud was better than another; there has been no way to determine whether the lubricating property of any mud was adequate. Furthermore, in comparing lubricating properties of muds, two qualities should be considered—the coefficient of friction and the extreme pressure lubricating qualities. For example, graphite can reduce the coefficient of friction, but does not provide a strong enough lubricating film to protect bit bearings subjected to high loads. The primary object of the investigation reported herein was to find a method for prolonging the bearing life of drill bits; and considering the heavy loads used in drilling, it was felt necessary to find both a method of measuring lubricating properties under high load conditions and methods of establishing effective lubrication under these conditions. Extreme pressure lubrication deals with that realm of lubrication where metal surfaces are in rubbing contact with each other under very high pressures. These pressures may be quite high, as for example,

By Arnulf Muan, Klaus Schwerdtfeger, L. S. Darken

The equilibrium ratios CO/CO, of. a gas pharse coexislitrg with selected oxide phase assemblages of the system "FeO"-MnO-SiO2 and metallic iron have been determined at 1150°C. The data obtained are combined wilh the thermodynamic data on fayalite (Fe2SiOi) and "FeO"-MnO solid solutions to derive the actiuity-coiHpQsiUon relations for the various oxide components of the olivine (Fe2SiO4-,lIn2Si04) and pyroxenoid (standard free energies of formation of ferrosilile (FeSiO,), tephroite , and rhodonite from their oxides have been evaluated. In the Appendicc~s by Darken and Schivcerdtjeger the ideal solution law for solid deri oed. THE quaternary system Fe-Mn-Si-0 is of interest in many important phases of iron- and steelmaking. For instance, reduction in the blast furnace of ores containing Fe-Mn silicates, deoxidation of liquid steel with Mn-Si alloys, formation and subsequent alterations of oxide inclusions in steels, and the oxidation of manganese- and silicon-containing steels during processing are all governed largely by the thermodynamic relations existing in the system Fe-Mn-Si-O. Although the equilibria between liquid iron and Fe-Mn silicates have been studied in detail previously,' experimental data are lacking for equilibria among the solid phases. The present paper provides such data for a selected temperature, 1150°C. I) PREVIOUS WORK The general feature of subsolidus phase relations in the system iron oxide-manganese oxide-silica under strongly reducing conditions is known.' Fig. 1 represents an isothermal section at 1150°C. The phases present, in addition to silica, are three solid-solution series, viz.. the oxide solid-solution series "Fe0"-MnO (manganowiistite), the orthosilicate solid-solution series Fe2SiO4-Mn2SiO4 (olivine), and the metasilicate solid-solution series FeSi03-MnSi03 (pyroxenoid). The "FeO"-MnO solid solution is continuous between the two end members wustite and manganositeV3 as is the orthosilicate series between the two end members Fe2SiO4 and Mn2 SiO4. 4 The metasilicate solid solution, on the other hand, extends only part way from the MnSi0.l end member toward FeSiO because the latter end member is unstable relative to the phase assemblage fayalite + silica. The letters a through m in Fig. 1 identify end members of the individual solid phases dealt with [for instance, c designates fayalite (FeZSiO4) in equilibrium with wustite, e designates fayalite in equilibrium with silica, and so forth]. In addition to the phase relations there are some thermodynamic data available for the end members of the silicate solid solutions and for the "FeO" -MnO solid solution. The free energy of formation of fayalite (Fe2Si04) from its oxide components has been studied by several investigators. Unfortunately, the data obtained show considerable discrepancies. According to recent tabulations by Kelley, the change in standard free energy Affn for the reaction

Jan 1, 1967

By J. K. Elliott, P. H. Kelly

Many engineering functions such as surface metering work and laboratory compressibility check points involve the use of liquid densities of light hydrocarbon mixtures at various pressures and temperatures. However, at the present time, no simple reliable method exists for determining density variation, particularly if the composition of the liquid is unknown. Consequently, a study was undertaken to develop and present a simple and accurate method of predicting density variation of a light hydrocarbon liquid with pressure and temperature, knowing only the density of the liquid at some condition. The experimental liquid compressibility data from API Project 37 by Sage and Lacey' have been considered to be accurate within 0.5 per cent and cover a wide range of pressure (14.7 to 10,000 psia), temperature (100" to 400°F) and molecular weight (up to 150). From these data, a set of liquid density curves, which relate density to pressure, temperature and molecular weight, was developed. These curves make it possible to predict density variation with pressure and temperature. Compared to extensive laboratory compressibility data on a complex, light hydrocarbon liquid, the use of the liquid density curves resulted in an average error of less than 0.5 per cent. Based on the results of this analysis, it is concluded that the set of liquid density curves developed from the data of Sage and Lacey provides an accurate and simple method for predicting the density variation of light hydrocarbon liquids when the density at some condition is known. These curves should be very helpful in many engineering calculations, particularly in the surface metering of light hydrocarbon liquids. INTRODUCTION Many situations arise in field and engineering laboratory work, such as reservoir engineering studies, check of experimentally determined laboratory data and orifice flow-meter formulas, where liquid density factors at various pressure-temperature conditions are required. Also, the need for accurate light hydrocarbon liquid information has become more important with the advent of miscible-type displacements for secondary recovery purposes in oilfield operations. Several reliable methods are available1 - "or determining the density of liquid hydrocarbons if the composition of the liquid is known. However, there is a definite lack of methods for accurately determining the variation of density when the composition of the liquid is unknown. The purpose of this study is to review the various methods for determining hydrocarbon liquid densities and to develop a simple and reliable method of determining variation in density of light hydrocarbon liquids with pressure and temperature when the compositio~n of the liquid is unknown. METHODS FOR DETERMINING DENSITY OF LIQUIDS OF KNOWN COMPOSITION Sage, Lacey and Hicks' have proposed a method to predict the density of light liquid hydrocarbons by using partial molal volumes. Data are available on experimentally developed partial liquid volumes of hydrocarbons over a rather limited range of temperature, pressure and composition. The partial mold volume method has proved satisfactory for determining the density of some hydrocarbon liquids when the composition is known. Within the range covered in the Sage, Lacey and Hicks1 data, the results agree within about 3 per cent of the experimental values. Hanson mentions the limitation of this method to a composition range of approximately 10 per cent by weight of methane, which will not allow this correction to cover most low molecular weight-light hydrocarbon liquids. Standing and Katz2 studied data on light hydrocarbon-liquid systems containing methane and ethane at high temperature and pressure and have presented a method for determining liquid densities, assuming additive volumes for all components less volatile than ethane and using apparent densities for methane and ethane. The compressibility and thermal-expansion curves used by Standing are based on assumptions that compressibility of a hydrocarbon liquid at temperatures below 300°F is a function of the liquid density at 60°F and that thermal expansion of the liquid is affected little by pressure. The information required to use this technique with an example problem is furnished by Standing.' Hanson eports an average error of - 0.5 per cent using the method of apparent densities in calculating liquid densities of several volatile hydrocarbon mixtures. However, as implied, the apparent density method is not applicable for liquid density calculations when the composition of the liquid is unknown. Watson- as presented a method

By O. D. Gonzalez, R. A. Oriani

The migration of hydrogen and of deuterium dissolved in a iron and in nickel induced by an applied electrical potential has been measured over a range of temperature. In all cases the intevstitial solute migrates to the cathode. The deduced charges of transport are positive, independent of or slightly dependent on temperature, and markedly larger for the beavier mass. These results are interpreted as showing predominating momenlum transfer between electron holes and the acticated complex of the elementary diffusion act. THE migration of components of a metallic system under the action of an applied electrical potential is called electromigration or electratransport. This phenomenon has been known since 1861 and is becoming of increasing interest to metal scientists because of the opportunity it offers to explore the connections between mass transport and electron transport. Until recently, reliable data have been very scarce because of the difficulties inherent in the measurements. In addition, theoretical understanding has been either nonexistent or rudimentary and at present this situation is only moderately improved. Whereas the prevalent notion had been (and is still widely held) that electromigration directly reflects the state of ionization of a component in its equilibrium state in the metal, it is now becoming clear that the major, if not the sole, factor is the momentum transfer between the charge carriers, themselves in motion because of the applied field, and the metallic components in their activated state. It becomes important therefore to obtain reliable data over a range of temperature on electrornigration in systems where only one component moves, that is, either in pure metals (vacancy motion) or in interstitial solid solutions. Our choice of systems was dictated by our measurements' of the Soret effect in solid solutions of hydrogen and deuterium in iron and nickel. Indeed, the main reason for making measurements on electromigration at all is that we suspected a mechanistic connection between thermal diffusion and electromigration. In our former paper' we proposed that what is common to both phenomena is the momentum transfer between the charge carriers in the metal and the atom as it jumps from one equilibrium position to another. In Soret diffusion, however, there is in addition the dissipation of the energy of the jumping atom via pho-nons. We have since found that similar ideas for this interconnection have been developed by Fiks.2 In the present paper we will discuss in more detail our understanding of the mechanism of electromigration itself, but we turn first to a presentation of the experiments and results. EXPERIMENTS AND RESULTS Let us consider a metallic system in which only one component, m , can diffuse at the temperature in question, and to which are applied simultaneously an electrical potential gradient, grad cp, and a chemical potential gradient, grad p,. One may write3 for the material flux Jm and the electronic charge flux Je by applying the condition, grad = 0, to Eq. [I] and evaluating grad(um/T) for the isothermal case and for activity coefficient independent of composition over the relevant composition range. ere, c is the concentration and D the diffusivity of the mobile component. The experimental technique that we have used to measure Z* for the isotopes of hydrogen individually dissolved in each of two transition metals has been developed from a method used by Wagner and Heller.4 (A somewhat similar method was also used by Herold.5) A wire or thin rod of the metal under study connects two volumes of gaseous hydrogen at identical pressures; hence grad um = 0. A direct current is made to pass through the specimen, serving to maintain the specimen at the desired temperature and also to produce the electromigration of the hydrogen dissolved in the metal in equilibrium with the gaseous hydrogen. Under these conditions, the flux of hydrogen from one volume to the other is given, from Eq. [I] for grad( um/T) =0by

Jan 1, 1968

By S. K. Tarby, A. E. Rathke

A series of 0.21 pci carbon steel melts was processed under conditions which sinzulated industrial vacuum degassing practices. The results indicated that the efliciency of carbon deoxidation was not significantly increased by vacuunz degassing at chamber pressures below 100 torr. The two common phenomena of vacuu?rz carbon deoxidation, i.e.,an unaccountable carbon loss and nonequilibrium deoxidation, were observed in this study. Oxygen from adsorbed and chenzically combined water of certain refractory structures was probably responsible for the unaccourttable carbon drop. The nonequzlibrium deoxida-tion phenomenon was adequalely correlated by a bubble nucleation )mechanism. Although the major limiting factor in oxygen removal from these melts contained in alumina crucibles could not be eslablished with certainty, there was some evidence that melt-refractory interactions were not the most imp-portant limiting factor. CONSIDERABLE attention has been given to the use of carbon in conjunction with vacuum degassing* as a deoxidation technique. An advantage of utilizing the reaction: for deoxidation purposes is that the reaction product is gaseous and hence can be effectively removed from the melt. The pressure dependency of Eq. [I] is illustrated by the expression for the thermodynamic equilibrium constant, K: where pco is the partial pressure of carbon monoxide, and hc and ho are, respectively, the Henrian activities of carbon and oxygen in molten steel. The thermodynamic aspects of the C-0 reaction in molten Fe-C-0 alloys have been investigated extensively. Theoretical predictions based on Eq. [2] indicate that, at the chamber pressures reached in most vacuum degassing units, carbon is the most effective deoxidizer that can be obtained. ~ornak' has compared industrial values of carbon and oxygen contents with values calculated from theory, and has shown that on a commercial basis vacuum carbon deoxidation has not achieved the degree of success that is predicted from theoretical considerations of the C-0 reaction. The use of industrial units to determine the reason for this deviation is limited by economics and the difficulty in controlling critical process variables. Therefore laboratory-sized experiments in this area should yield valuable information. Investigations on vacuum carbon deoxidation have been made in experiments conducted under vacuum induction melting conditions where pressures in the low micron range, long treatment times, exposure to highly dynamic vacuum conditions, and relatively quiescent bath conditions have been maintained. Inve~ti~ators~~~ utilizing vacuum induction melting for the production of high-purity iron by carbon deoxidation have concluded that, under high vacuum conditions, melt-refractory interactions constantly supply oxygen to the melt. Thus, it has been suggested that refractory dissociation rather than the C-0 reaction establishes the lowest possible oxygen content of a melt. Thomas and Moreau4 and Bennett, Protheroe, and ward5 have studied carbon deoxidation of Fe-C alloys melted in magnesia crucibles under reduced pressures. The results of both studies differed considerably from those predicted by theory, and the limiting factor for lower oxygen contents was attributed to melt-refractory interactions. Because the refractory material used in both investigations was magnesia, no quantitative measure of crucible reaction could be made and the amount of carbon consumed during treatment was utilized as an indirect indication of the extent of crucible reaction. Studies concerned with the determination of the minimum oxygen content attainable in pure iron and Fe-C alloy material vacuum-induction-melted in various crucible materials have also been conducted."' The suitability of a crucible refractory for use under vacuum conditions appears to depend upon the major refractory constituent, the amount and type of impurity oxides, and the composition of the melt. Meadowcroft and Elliott9 have presented a theoretical review of the type and extent of various refractory reactions which should be considered in vacuum refining. ~eadowcroft'~ considered the vacuum degassing process to be dependent upon the "carbon boil" and has proposed that a part of the difference between the experimental and theoretical results for vacuum carbon deoxidation may be explained by the fact that the pco term in Eq. [2] is not the chamber pressure. Instead, the pressure necessary for carbon monoxide

Jan 1, 1970

By P. N. Nikiforuk, M. C. Zoerb

An analog computer model has been developed of a large, walking dragline. This model permits changes in the configuration of the dragline, or changes in its digging cycle, to be readily investigated on a computer. As a result, it is possible to determine ways of improving the dragline performance without having to carry out many expensive experiments on the actual machine. This model can also be used as a simulator for dragline operator training. INTRODUCTION The movement of large volumes of material is always a costly operation and the increasing trend towards the use of large capacity equipment does not necessarily produce the economy desired. This problem is of particular concern to the mining industry in its open-pit operations where large draglines are used to remove the overburden. Since this operation often forms a very large part of the total cost of open-pit mining it is particularly important that it be kept as economical as possible. This requires that the dragline be operated at its maximum possible productivity. In order to optimize the operation of a dragline it is desirable that a suitable model of the dragline be available for study. Once such a model is available the effects of changes in the configuration of the dragline, or changes in its digging cycle, could be readily investigated on a computer and methods of improving its performance could be determined without having to carry out expensive experiments on the actual machine. Such models would be of use not only to dragline users, but also to the manufacturers of draglines for design purposes. Presumably studies of this type could eventually lead to computer-con trolled stripping operations. As far as can be determined very little work1'2 has been done on the derivation of dragline models. Such models must be derived before work of the nature indicated in the previous paragraph can be carried out. As a small contribution to this field this paper describes work that has been carried out by the senior author during the past few years on the derivation of one such model. This model describes the drag, hoist and swing operations of a 35 cu yd, commercially available, walking dragline. This model was derived for use with an analog computer and can be used either for optimization studies, or for the training of dragline operators. DESCRIPTION OF DRAGLINE General: The particular dragline for which an analog computer model was developed is owned by the Great West Coal Company of Estevan, Saskatchewan. It is a walking dragline and its function is to remove 40

Jan 1, 1967

By N. Plint, E. T. Pinkney

Copper segregation developments within the Anglo American Group are described with particular reference to the work of the TORCO Project Team. Some radical differences have been adopted in the mechan-ical approach as compared with methods used by previous workers. Most important among these are the following: (1) The heating unit is a fluidized-bed roaster. (2) A moving bed segregation chamber is used in which a packed bed of heated ore, to which reagents have been added, moves slowly downwards under gravity and is discharged at a controlled rate by means of a suitable valve mechanism. (3) The quantities of reagents used (sodium chloride and some form of carbon) are only a small fraction of those recommended by other investigators. (4) Retention times of only 2 min have been satisfactory for the complete segregation of many types of ore. Project TORCO is the code name given by Anglo American Corp. of South Africa, Ltd., to the investigation undertaken in Zambia to develop a plant capable of treating large tonnages of refractory copper-bearing ores by the segregation process. The process has now been under investigation for six years (1960-66) and during that period work has progressed from bench-scale laboratory tests through small-scale pilot-plant work to a 500-tpd commercial pilot unit. Test work has been carried out mainly on the ores of the Zambian Copperbelt and results have been encouraging. Recoveries on the order of 90% have been obtained on all Copperbelt ores tested with high concentrate grades. PRELIMINARY INVESTIGATION Laboratory scale segregation test work on Nchanga banded sandstone ore and oxidized ore from Kansan-shi was originally conducted by Professor Maurice Rey. SociCtC Miniere et Mktallurgique de Pefiarroya, in 1957. Both ores responded favorably to the process. In November, 1960, the consulting engineers of Anglo American Corp. of South Africa, Ltd., initiated preliminary investigations on the application of the segregation process to Copperbelt ores. In early laboratory work, samples were heated in a closed, externally heated tube furnace for one hour at temperatures ranging from 800 to 830°C. The free moisture content of the ore feed was 4% and reagent additions were 0.5% salt and 1.0% coal. Both salt and coal were pre-crushed to -48 mesh Tyler.* All samples except an ore from the Bancroft mine responded well to the process and recoveries of over 90% were obtained. The presence of dolomite was responsible for the relatively poor recoveries from the Bancroft ore, as tests on one of the other ores with additions of calcite or dolomite confirmed. Recoveries dropped from 94 to 45% Grades of all the rougher flotation concentrates produced were low at 10-15% copper. The laboratory investigation was then followed up on a pilot-plant scale. An attempt was made to carry out the segregation process in a coal-fired 5-ft-diam fluidized-bed roaster. The feed consisted of a banded sandstone from Nchanga to which 0.5% sodium chloride had been added. Reducing conditions were maintained in the roaster by the addition of coal in excess of that required for heating the feed. The excess coal and added salt constituted the segregation reagents and the temperature of the bed was held at approximately 750°C. Recoveries of copper and grades of concentrates were low. The reasons for lack of success were: (1) The reaction is a gaseous one and a fluidized-bed roaster is therefore unsuitable as a segregation unit. (2) Coal in excess of that required for heating was introduced into the fluidized bed and under these conditions some of the copper would be reduced in situ. The term "reduced in situ" is used here to describe any copper mineral reduced to the metallic state in its original site within the ore particles. If this copper is to be recovered by subsequent froth flotation it must first be liberated from the gangue by grinding or some other method which, in practice, is usually not feasible on ores of this type.

Jan 1, 1969

By Hsun Hu, B. B. Ruth

The rates of migration of high-angle boundaries in zone-refined aluminum crystals rolled 20 to 70 pct in the (110)[i12/ orientation were studied. Following a recovery anneal at an appropriate temperature to stabilize the polygonized structure, boundary migration rates of artificially nucleated grains were measwed isothermally at several temperatures. Results indicate that the rate of boundary migration depends strongly on the amount of deformation and on the cell size of the polygonized matrix, and is related to the driving free energy by a power function. The degree of anisotropy in growth 0.f the re crystallized grains nn&apos;th preferred mientation is independent of deformation; the migration rates of the fast-moving and the slow-moping boundary segments of a gowing grain differ by as much as one order of magnitude. The actir\ation energy fm a grain boundary migration, although nearly the same for both the fast-moving and the slow-moving boundaries for a given deformalion, decreases from 45 to 30 kcal per mole with an increase in deformation from 20 to 70 pct reduction. Re crstallization by the growth of the artificially nucleated grains results in preferred orientation. The Percentuge of&apos; grains favorably oriented for growth increases with increasing deformation. None of these grains corresponds to the ideal Kronberg-Wilson orientation relationship. The observed growth aniso-tropy is discussed in terms of boundary structure. The boundary velocity as a function of the cell inter -facial area, or the driving free energy, is discussed in the light of current theories of boundary migration. It is well established that recrystallization with re-orientation occurs by the migration of high-angle boundaries of strain-free grains. The driving force for this process is provided by the free energy stored in the metal during deformation. A quantitative study of the effect of varying driving force on grain boundary migration in deformed metals has not been possible heretofore, primarily because of: 1) concurrent recovery steadily decreasing the available driving free energy for boundary migration, &apos;-3 and 2) in-homogeneity of strain in the deformed metal.4 Aust and Rutter3 studied grain boundary migration in striated single crystals of zone-refined lead. Although the driving free energy in such crystals remains unaltered during annealing, this method does not provide a range of driving free energies over which measurements of grain boundary migration can be made. In the present investigation, the rates of migration of high-angle boundaries in deformed aluminum zone- refined single crystals were studied at various temperatures, after deformation ranging from 20 to 70 pct reduction by rolling at -78°C in the (ll0)[i12] orientation. The boundary migration rates along different crystallographic directions were determined under steady-state conditions, i.e., in the absence of competing recovery processes or impingement of recrystallized grains growing into the deformed single crystal matrix. Simultaneous recovery was eliminated by suitable anneals prior to the boundary migration measurements. The recrystallized grains, which grew a ni so tropically into the homogeneously polygonized matrix, developed flat boundary segments during early stages of growth. These boundary segments subsequently migrated along a direction approximately normal to the boundary plane into the matrix rystal. Increasing deformation over the range employed was estimated to increase the driving free energy for boundary migration by about five times. The kinetics of the boundary migration process, examined under these conditions, indicate that the boundary velocity is greatly affected by a small change of the driving free energy in the matrix crystals. These results were examined in the light of the current theories of grain boundary migration. EXPERIMENTAL PROCEDURES Single crystal strips (9 by 1 by 0.125 in.) of zone-refined aluminum, were seed-grown by the Bridgman method in a high-purity graphite mold (<lo ppm ash) at 1 in. per hr. Precautions were taken to minimize contamination of the metal during crystal preparation and subsequent handling. Spectrographic analysis of the metallic impurities in the grown crystals is Qven in Table I. The crystals were rolled in the (110)[112] orientation at -78°C to various reductions in thickness, ranging from 20 to 70 pct, in 10 pct increments. The desired reduction was achieved by many rolling passes, each being no more than 0.002 in. To minimize surface friction, the crystal was rolled between two thin layers of teflon. For those crystals rolled more than 40 pct, it was necessary to remove the disturbed surface layers by electropolishing at -5" to -10°C at an intermediate stage of rolling. The edges of deformed crystals were removed by a jeweler&apos;s saw while submerged in alcohol at -78° C to obtain samples of about ? by i in. The distorted metal at the cut edges and the surface layers were then removed by electropolishing, with removal of a minimum of 0.004 in. from each surface. The thickness of the crystals prior to rolling was chosen so that the final thickness was 0.025 in. for all samples. These deformed single crystals were each prean-nealed for 1 hr at an appropriate temperature in the range of 130" to 280°C, depending upon the amount of deformation. The purpose of this preannealing was to

Jan 1, 1970

By E. P. Klier, Jane Jellison

The course of the transformations on cooling in a series of Cu-A1 alloys has been followed by means of thermal analysis, X-ray diffraction, and optical metallogrAphy. Specimen size was found to have a pronounced effect on the time of completion of the pearlite transformation due to modification in the unmber of nuclei. The ordeving transformation zvas observed for all cooling rates for all alloys studied, and is considered insuppvessible. At intermediate cooling&apos; rates, a transitional struchire for 011 alloy compositions and based ox the bainite reaction was observed. The alloys of aluminum and copper possess a eu-tectoidal region at from about 9.4 to 16 pct Al. Above the critical temperature, the disordered bcc phase, 0, is stable. At subcritical temperatures, the 0 phase decomposes under equilibrium conditions into the stable phases, a and y,. At the eutectoid composition, these phases develop a lamellar pearlite structure, in appearance much like that found in steels. During rapid quenching from supercritical temperatures, the diffusion processes required for the formation of the equilibrium phases are suppressed, and the 0 transforms to an acicular martensitic phase, designated B or y&apos; depending on the aluminum content, or to a mixture of the two. Early investigators of this system concentrated on determining the equilibrium phase diagram, &apos;&apos;&apos; the structure of the equilibrium phases,3 and later the isothermal decomposition of the 0 phase at subcritical temperatures, or the decomposition of the martensite on tempering. Other work which treated the nonequilibrium transformations concentrated largely on the martensite transformations and the structure of the b and y&apos; phases.15"17 Since the pearlite and martensite transformations in this system resemble those in steels, the one being formed by nucleation and growth of adjacent lamellae of the equilibrium phases, and the other by a diffusionless, shear-type transformation, the present study has been undertaken to investigate the structures formed on continuous cooling, with the objective of constructing cooling transformation diagrams analogous to those derived for alloy steel systems. In particular it was desired to determine whether parallels between the A1-Cu and iron-iron carbide systems on continuous cooling extended to encompass the intermediate phase, namely the bainite structures, which form by a combined martensitic-diffusion mechanism,&apos;&apos; but which are supersaturated with respect to the equilibrium phase upon which their crystal structure is based. To this end the transformations in representative alloys of the Cu-A1 B eutectoid system have been followed, using the methods of thermal analysis, X-ray diffraction, and optical microscopy. EXPERIMENTAL PROCEDURE Materials. The materials tested were those used by Klier and Grmko 5 in their study of the isothermal decomposition of these alloys. The alloys were melted in graphite crucibles and chill-cast into 5 by 2 by 1/4 in. ingots. The ingots were then hot-rolled to sheets 1/8, 1/16, and 1/32 in. thick. The compositions are given in Table I. Specimen Preparation. Specimens for thermal analysis were in the form of small squares approximately 1/4 in. on a side. For the 11.9 and 13.5 Al-Cu alloys, the sheet thickness was 1/16 in. and for the 10.5 A1-Cu alloy the 1/32-in. sheet was used. One lead of a chromel-alumel thermocouple was spot-welded to either face of the square with the specimen itself thus forming the hot junction of the thermocouple. Specimens for X-ray analysis were prepared as filings from martensitic sheet stock and were sealed in glass ampoules under 1/2 atm of argon. Metallographic examination of the powder specimens was carried out on a portion of the heat-treated filings mounted in lucite. All specimens were examined metallographic ally after etching

Jan 1, 1965

By R. J. Burgess, A. R. Ramey, A. R. Adams

During interpretation of pressure buildup tests on gas wells in a tight dolomite gas reservoir, peculiar behavior was noticed. Two straight lines were apparent. Effective permeability to gas taken from either straight line was about the same, and the Miller-Dyes-Hutchinson dimensionless time check for the straight line was proper for both straight lines. Geological data indicated the likelihood of scattered trending fractures in the reservoirs. Since the first straight Iine yielded permeability values close to the geometric mean permeability from core analyses, it was postulated that the reservoir model was that of an acidized well completed in the tight dolomite, but that widely scattered hairline fractures caused the mean permeability of the reservoir distant from the well to be higher than the matrix permeability. Because all other studies of fractured reservoirs to the authors&apos; knowledge assumed that the fracture matrix was dense enough to communicate directly with the well, no interpretative methods were available. The Hurst line-source solution for a radial change in permeability for interference between oil reservoirs was adapted to pressure buildup testing. The result indicated that the first straight line should yield the proper matrix permeability and wellbore skin effect. The second straight line may be extrapolated to obtain static pressure. The time of bend between the straight lines was used to estimate distance to a fracture. Application to field test data is shown. It is believed that the methods developed and the case history presented will add to present tools available for pressure buildup interpretation. Introduction Since the pioneer studies by Miller, Dyes, and Hutchin-son1 and Horner&apos; in 1950 and 1951, well test analysis has become recognized as one of the most powerful tools available to both production and reservoir engineers. Well test analysis serves as a logical basis for well stimulation and completion analysis, and for long-term reservoir engineering. Since the early 19501s, much effort has been placed on the development of well-test analytical methods. Reservoir and well conditions of increasing complexity have been considered systematically to provide the analyst with a catalog of causes and effects. Matthews and Russella state that some 200 papers dealing with this subject have been published in the last 35 years. Developments in well test analysis appear to have originated in one of two ways. Either a physically realistic field condition was anticipated and analytical solutions for the condition achieved, or anomalous field test behavior was recognized and interpretative methods sought for the anomaly. In recent years, it has appeared that the latter has inspired an increasing number of studies. The analyst today finds an increasing number of known cause and effect studies available for well test analysis, the classic of which is that of finding the specific flow problem that generated the answer — the well behavior. Although it may be impossible to achieve this goal uniquely, the analyst often is able to select a useful interpretation that combines all known performance and geologic data — or to show that various logical alternatives would not significantly affect the interpretation. During a recent reservoir study, we observed gas well test behavior that did not appear to fit behavior described previously. Although it cannot be said that we have found a unique interpretation, we shall present in this paper the peculiar behavior observed, and describe the reservoir and interpretative methods developed. Reservoir Description The subject gas reservoir is a 9-mile-long, narrow dolomite reservoir lying within a limestone of Ordovician age. (See Fig. 1.) The dolomitized rock in the field consists of dark brown to buff, dense to coarsely crystalline, vugular dolomite containing numerous hairline fractures, many of which may have been closed in the reservoir and parted when cores were brought to the surface. Larger fractures are also apparent in core, but usually are filled and sealed with euhedral dolomite crystals. Portions of the north flank of the reservoir are known to be cut by a sealing fault downthrown to the north. Gas wells located near the fault have higher open flow potentials than those more distant from the fault. This is believed to be a result of higher permeability near the fault due to more extensive and open fractures. Detailed coring and core analysis have been performed on several of the wells in this reservoir. Fig. 2 presents permeability variation&apos; plots for both horizontal and vertical

Jan 1, 1969

By R. L. Moment

The anisotropic elastic behavior of rolled beryllium sheet has been measured, using a pulse echo technique, and compared with X-ray diffraction data. Calculated elastic stiffness constants compared favorably with published values for beryllium single crystals which were attributed to the strong (0002) rolling plane texture. Variations of Young&apos;s modulus in the yolling plane could be associated with the velative distribution of (0002) planes out of their ideal position in the rollitzg pkule. WHEN a metal is subjected to cold working such as drawing, forming, or rolling, a crystallographic texture develops which can significantly alter its physical properties. One method for detecting this texture is X-ray diffraction, but Alers and Liu&apos; have recently pointed out how the prediction of anisotropic physical properties from pole figures alone is not always accurate due to differences in interpretation. Variations in Young&apos;s modulus with orientation or, more completely, the values of the effective elastic constants of the worked metal, also serve to indicate the presence of a texture. In fact, as Alers and Liu&apos; pointed out, calculated variations in Young&apos;s modulus for assumed orientations, when compared with experimental data, can be used to eliminate some of the uncertainty in interpretation of X-ray pole figures. Thus, elasticity measurements can serve not only to clarify any unusual elastic behavior of worked metal, but also to detect and in part determine the nature of its texture. X-ray determination of the texture of rolled beryllium has been reported by Smigelskas and Barrett,2 who found a strong texture of (0002) in the rolling plane with (1070) planes normal to the rolling direction. In the case of metal rolled at room temperature, they reported that [1010] directions also appeared at positions 60 and 120 deg from the rolling direction in the rolling plane, while in more recent work Keeler3 found these directions were also tilted towards the rolling plane. The texture for beryllium rolled at 80O0C, however, only showed (1010) planes normal to the rolling direction and the spread of (0002) planes out of the rolling plane was less. In looking for elastic anisotropy one might consider unidirectional rolling of a metal as introducing an or-thorhombic symmetry through reorientation of the grains, since the three deformations, compression, extension in the rolling direction, and extension in the cross direction, are orthogonal to each other and unequal in magnitude. Thus the rolled sheet could be treated like an orthorhombic single crystal and the nine stiffness constants of the elasticity tensor used to calculate the anisotropy of Young&apos;s modulus, the shear modulus and Poisson&apos;s ratio. In this case we could write: which is symmetric about its diagonal. Borik and Alers4 have recently used this approach on rolled die steel with very good results. They found, however, that instead of displaying orthorhombic elastic symmetry their specimens could be considered tetragonal in which case Cr1 = c22, c13 = Ca, and c44 =cjj. This conclusion was made solely on the basis of the measured tensor elements, and serves to point out the advantage of this method for studying the anisotropy of rolled metals. Their calculated values for Young&apos;s modulus as a function of angle in the rolling plane also checked very well with direct measurements made on different specimens using the resonance technique. In the present study, cross-rolled beryllium was used which had been unidirectionally rolled about 11 pct for the final reduction. This imparted a slight anisotropy in the rolling plane which was detected both by X-ray techniques and elasticity measurements. For purposes of discussion in this paper, the rolling direction is that direction in which the most reduction passes were made and cross direction is the normal to the rolling direction in the rolling plane. It was also decided to consider the rolled sheet as displaying orthorhombic symmetry for the purpose of obtaining elasticity samples with the direction defined as in Table I. Any change in the final symmetry attributed to the sheet would then be made on the basis of the measured elastic stiffnesses. The final data would then be compared with that expected from the X-ray study and that reported for beryllium single crystals. EXPERIMENTAL PROCEDURE Rolling Schedule. The samples used in this study were taken from a large sheet which, because of its size, had to be unidirectionally rolled for the final reduction. The resulting texture was that of cross-rolled metal with a slight unidirectional texture superimposed. A cast beryllium ingot, 9.500 in. sq by 3.325 in. thick, was cross-rolled to 81 pct reduction followed by unidirectional rolling for an additional 11 pct to give a total reduction of 92 pct. The thickness of the final sheet ranged from 0.265 to 0.280 in. Reduction up to 67 pct was done at 980°C and the final 25 pct at 870°C. Analysis for metallic impurities showed aluminum 0.06 pct, iron 0.19 pct, and silicon 0.11 pct, giving a beryllium purity of 99.64 pct.

Jan 1, 1968

By Charles R. Stewart, William W. Owens

Reservoir performance predictions based on laboratory core test data assume that fluid flow is laminar for the laboratory test. A study has been made to determine the validity of this assumption for laboratory tests on various types of porosity found in producing limestone formation. Data are presented which show that turbulence and slippage can occur during laboratory tests on hetero-geneous-type porosity limestones, thus causing serious errors in measured single-phase permeabilities and two-phase relative permeability characteristics. In single-phase flow tests it is possible to eliminate turbulence and correct for slippage or to eliminate both factors by controlling test conditions. It is not always possible to control test conditions and thereby eliminate turbulence and slippage in two-phase .flow tests. A correction method is presented which can be used to calculate the true two-phase laminar flow relative permeability characteristic even though furbulence and slippage exist. .INTRODUCTION It is customary to make use of Darcy&apos;s law and modifications of this law, together with laboratory data on formation core samples to predict the performance of producing reservoirs. Such predictions are based on an assumption that fluid flow is in the laminar or streamline region for the laboratory test. It was the purpose of this inves- tigation to determine the extent to which turbulent flow may occur in laboratory fluid flow tests on hetero-geneous porosity limestones. Considering that turbulent flow conditions might exist in some laboratory fluid flow tests, additional emphasis was placed on the development of a method to correct for turbulence when laminar flow conditions could not be attained. FLUID FLOW CONCEPTS FOR POROUS MEDIA The Influence of Pore Geometry on Fluid Flow One of the more important factors influencing fluid flow in porous media is the geometry of the pore space which includes such characteristics of the pores as size, shape, distribution, roughness, uniformity, etc. In general, oil- and gas-producing formations can be divided into two broad types on the basis of pore geometry. One has been called sandstone-type porosity media, which is characterized by a small range in pore size, uniformity in shape of the pores, smooth pore surfaces and a regular and uniform distribution of pores. The other type has been called heterogeneous porosity media and is usually limited to the dolomites and limestones. This type is characterized by a wide variation in the size, shape, and distribution of the pores and rough, irregular pore surfaces. It is therefore apparent that conditions are much more favorable for turbulent flow* in heterogeneous-type porosity media than in sandstone-type porosity media. Interrelationship Between Turbulence and Gas Slippage In studying the problem of turbulent flow in laboratory tests on porous media, it is necessary to be aware of the interrelationship between slippage and turbulence for gas flow. As a result of slippage or the Klinken-berg effecta, apparent perrneabilities to gas are greater than the true value because there is no stationary layer of gas in contact with the walls of the flow channels. Gas slippage decreases as the mean free path of the gas molecules decreases. Since the mean free path of any gas decreases with increasing density, increases in static pressure result in lower apparent gas permeabilities. However, a reduction in gas permeability can also be due to turbulence. Therefore, in studying only turbulent flow in porous media, it is necessary to hold gas density, and slippage, constant or to reduce slippage to a negligible value by operating at high static pressures. Presentation of Laminar and Turbulent Flow Data A graphical relationship between permeability and a pseudo-Reynolds number, N,, will be used to show the two types of fluid flow, i.e., laminar and turbulent. The usual graphical method for such a description has been the use of friction factor-Reynolds number charts4. On such a logarithmic diagram, the laminar region appears as a straight line having a slope of 45 degrees. As the friction factor decreases and the Reynolds number increases, the turbulent region is reached and appears as a deviation from the 45-degree slope line. However, in petroleum engineering literature resistance of por-

By Yoichi Ishida, Ching-Yao Cheng, John E. Dorn

Tile creep behavior of a iron was investigated over the range of temperatures from 375° to 1150°K. Apparent activation energies for creep, obtained by the effect of sudden changes in temperature on the creep rate, revealed the presence of- four distinguishable regions. The creep behavior of a iron in Regions II (480° to 774°K) and 111 (775° to 1045°K) was found to correlate well with a model where the creep rate is controlled by the nonconsevualive motion. of jogged screw dislocations. An anomalous increase in the apparent activation energy fur creep in Region III was found to be in harmony with the known decrease in the free activation energy for self-diffusion over the Curie transformation range. Furthermore, the creep rates in Regions II and III were found to increase with stress due not only to the effect of stress on the activation energy but also to an increase in density of mobile dislocations. The evidence suggests that pipe diffusion along moving dislocations could be a significant factor over the lower temperatures of Region II. ALTHOUGH the creep properties of iron and steel have been of interest to metallurgists for some time and an extensive literature is now available, most of the published information has not produced much toward understanding the fundamental dislocation processes controlling the creep behavior of these materials. Whereas the activation energy for high-temperature creep of metals usually agrees well with that for self-diffusion, Sherby, Orr, and Dorn1 deduced from the data of Tapsell -and c1enshaw2 that the activation energy for creep of Armco iron is about 78,000 cal per mole and more recently Lytton and sherby3 reviewed the data of Jenkins and Mellor4 and found an activation energy for creep of a iron of 80,000 cal per mole. Both values are substantially greater than the activation energy for self-diffusion in a iron. Because of the limited data suitable for purposes of identifying the controlling dislocation mechanisms of creep and the unusually high activation energies for creep quoted above, it was considered desirable to reinvestigate in detail the creep of a iron by more reliable techniques which would provide sufficient data for determination of the creep mechanisms. EXPERIMENTAL TECHNIQUE Creep specimens were prepared from 3/8 by 4 in. iron bar stock of the following composition by weight percent: 0.001 C, 0.0120 0, 0.001 N, 0.004 S, and 0.003 P. The as-received bars were cold-rolled to a thickness of 0.100 in., annealed under argon for 30 min at 1113°K, cold-rolled to a final thickness of 0.063 in., and machined into tensile specimens having gage section 0.250 in. wide and 1.70 in. long. Finally they were recrystallized under argon at 1113°K. Specimens to be crept above 1060°K were recrystallized at 1173°K. Both recrystallization treatments gave the same reproducible equiaxed ASTM No. 4 grain size. Creep testing was conducted in machines fitted with Andrade-Chalmers arms that maintained a constant stress to within ±0.2 pct of the reported values. Deformation over the specimen gage section was sensed with a linear differential transformer and recorded auto-graphically as a function of time. The strains deduced from these measurements were sensitive to ±5 x 10-5. During creep the specimens were contained in an argon-filled chamber which was immersed in a temperature-controlled molten tin bath. Creep temperatures, as measured by thermocouples attached to the specimen, were maintained to within ±10.25°K of the reported values. For activation-energy determinations rapid changes in temperature of about 12°K were obtained within 30 sec by direct resistance heating of the creep specimen, and maintained to ±1°K of the reported values. Observations of structural details of specimens before and after creep tests were made by electrolytic polishing and etching in acqueous ammonium persulfate. EXPERIMENTAL RESULTS A typical example of the determination of apparent activation energies, Q, by the effect of small changes in temperature is illustrated in Fig. 1. Q is defined by

Jan 1, 1967