Search Documents

By S. C. Miller, R. D. Carter, H. G. Riley

A merhod based on theoretical considerations is presented for estimating srabilized performance of gas wells from data obtained from short-term flow tests. Calculated results are presented for II gas wells evaluated with four-to five-hour flow tests. For seven wells on which comparison war available, satisfactory agreement between measured and calculated performance was obfained. Although the final results are favorable, diflerences in flow capacity k,h obtained from flow tests conducted at different rates were observed. This inconsistency did not seriously affect the calculated results in the test wells. INTRODUCTION In order to negotiate the sale of gas from newly drilled gas wells and to plan otherwise for the orderly exploitation of gas reserves, it is necessary to predict the long-term withdrawal rates physically possible for each well. Since low permeability gas wells do not stabilize except after a lengthy period of flow, the operator is unable, without wasting large amounts of gas, to measure directly the stabilized performance prior to contracting for sale of the gas and connecting the well to a pipeline. This problem is not new, but is one which has recently received considerable attention because of the rising importance of natural gas in the economy and new emphasis on the development of gas reserves. Thus, a method is needed for estimating stabilized performance from information obtained from short-flow tests in which stabilized flow is not attained. Such short tests would be suitable for wells not connected to a pipeline since the tests require only small amounts of time, labor and production of gas. Other requirements of such a method are short calculation time and use of a model of the gas well and its reservoir drainage area which is as consistent as possible with present theoretical, experimental and field knowledge. The method presented herein was developed in an effort to meet these requirements. It has been named the "two-Row method". In addition to the section on the evelopment and description of the method, a description of field application of the method and a discussion of well test results are in- cluded. Finally, there is a section on an additional consideration prompted by some of the test results. The Appendix contains certain details of the derivation of the two-flow method and a sample calculation for Test Well 1. MATHEMATICAL DEVELOPMENT AND DESCRIPTION OF THE METHOD The method is based upon the following assumptions: (1) The system can be described as a finite cylindrical gas reservoir produced by a single cylindrical well located at the center of the reservoir. Flow is strictly radial. (2) Flow obeys Darcy's law except very near the well where a quadratic flow law is assumed. (3) The system is also characterized by a laminar flow skin factor. The use of a quadratic or turbulent flow term is more fully justified in Refs. 1 and 2. When such a system is opened at the wellbore to constant rate flow from an initial condition of uniform pressure distribution, the radius or drainage quickly passes through the non-Darcy flow and "skin effect" region imrnediately surrounding the well. After the drainage radius or "transient" has passed beyond this region and before it reaches the external boundary, drawdown can be expressed by'." After the drainage radius reaches the external boundary, stabilized flow exists and drawdown is given by'," Eqs. 1 and 2 form the basis for the method. The validity of these equations will now be examined in the light of some well test data previously available. It has been found by actual measurements on gas wells that the relationship between drawdown and flow rate, that is the performance curve, is closely approximated by' Eq. 3 is a straight line on log-log graph paper. Eq. 2 and, for any lixed value of time, Eq. 1 plot as a curve which is slightly concave upward on log-log graph

By S. R. B. Cooke, Y. S. Kim, I. Iwasaki

lron oxides or the gangue minerals of iron ores can be made to float, depending upon the choice of collector. The selectivity of the separation is controlled by the accompanying chemical and operating conditions. In this report the findings on the general flotation characteristics of magnetite are summarized and the results are compared with those obtained for other iron oxides. Depending upon the choice of collector, either the iron oxides or the gangue minerals of iron ores can be made to float, the selectivity of the separation being controlled by the accompanying chemical and operating conditions. The mechanism by which a collector exerts preferential affinity for either group of minerals has been a subject of extended discussion. Previous investigations on the electrokinetic properties and flotation characteristics of hematite,' goethite,2 and quartz3 have demonstrated the electrical nature of the interaction of the collector with these oxide surfaces; anionic collectors are effective on positively charged surfaces and cationic collectors on negatively charged surfaces. Thus, from the differences in the electrokinetic properties of the above-named oxides, a theoretical basis has been provided for their separation by some typical collectors. Examination of the results reported for hematite and goethite shows that their electrokinetic properties and flotation characteristics are remarkably similar. It is of considerable interest, therefore, to examine similar information for magnetite, another common iron oxide, in order to consolidate the flotation characteristics of the three principal minerals of oxide iron ores. Such information is practically significant in the up-grading of both natural and artificial mag- netite produced from taconite, and in the treatment of partially oxidized magnetite ores. Information on maghemite and lepidocrocite, polymorphs of hematite, and goethite, respectively, is also of some theoretical interest in ascertaining the effects of crystal structure and of hydration in the iron oxide group. This article reports the findings on the general flotation characttxistics of magnetite, and compares the results with those obtained for other iron oxides. Experimentally, tlhe surface properties of magnetite were investigated by electrophoresis, and by measurement of the aclsorption of some typical flotation collectors. Its flotation behavior was tested by contact angle measurements and by simple flotation tests using a modified Hallimond tube. The basic information thus gathered was shown to be directly applicable to the interpretation of the results of batch flotation tests made on an artificial mixture of magnetite and quartz. EXPERIMENTAL MATERIALS Magnetite: Magnetite concentrate from Port Henry, N. Y., was further purified by fine crushing and magnetic separaiion, and was used throughout this investigation. For flotation tests the 150/200 mesh fraction was iso1ated and was cleaned by leaching briefly with warn 0.1 N HC1 followed by repeated washing with demineralized water. The final sample analyzed 71.71 pct Fe, 22.37 pct Fe" and 0.47 pct insoluble, and was stored in demineralized water. For e1ectrophoresis and adsorption measurements, the —400 mesh fraction of the magnetite concentrate was ground dry in a pebble mill. The surface area of the ground sample, as determined by the gas permeability method, was 70,000 sq cm per g. Since artificial magnetite has been reported to possess somewhat different physical properties than those of natural magnetite, it also was briefly examined. Preliminary flotation tests on a sample of artificial magnelite, prepared by reducing Brazilian hematite in a C0/CO2 mixture and analyzing 71.94 pct Fe and 23.78 pct Fe'', gave, however, nearly identical results with those using natural magnetite. It was concluded that artificial magnetite behaves similarly to natural magnetite. Samples which had been magnetically separated showed a tendency to cluster together due to residual

Jan 1, 1962

By Milton F. Goodrich

The idea of using water-only cyclones to scalp clean coal from the feed to other devices has recently been gaining in popularity. 1,2,6 An indication of this popularity is that water-only cyclone scalping is already employed in a few commercial plants. Water-only cyclones scalp the feed to a heavy media vessel in the Warwick plant,3 to Deister tables in the Rowland plant,4 and to a heavy media cyclone circuit in the Amherst Coal Company's McGregor plant.5 Despite this trend, we found surprisingly little information in the literature on the performance of water-only cyclone scalping. Therefore, we undertook a study of its technical merits, limiting our investigation to a water-only cyclone scalper in conjunction with heavy-media cyclones. We felt this to be a reasonable combination of equipment, since both processes operate well on the same size coal. First, we developed simple mathematical relationships revealing general characteristics of the system. Next, we compared the clean coal yield and quality obtainable by the system on four different coals with that obtainable by heavymedia cyclone circuit alone and water-only cyclones middlings recirculation circuit alone. This was done by using computer program to simulate the operation of the three systems. Finally, we considered the effect that implementation of the system would have on a typical preparation plant circuit. The System The system chosen for study [(shown schematically in Fig. 1)] comprises two stages: a primary stage (PI) consisting of water-only cyclones, and a secondary stage (P2) consisting of heavy media cyclones. The water-only cyclones separate raw coal at a low specific gravity, producing a relatively small amount of high quality clean coal and a refuse containing much good coal. This refuse is then recleaned in the secondary stage, where the sharper separation of heavy media cyclones, operated at a higher specific gravity, reclaims additional clean coal. Finally, the clean coal from both stages is combined as the system's clean coal, while the refuse from the secondary stage becomes the system's refuse output. Theoretical Analysis Two assumptions lie behind the following theoretical analysis. The first assumption is that partition curves can represent not only the performance of each stage but also the performance of the system as a whole. These curves show partition number as a function of specific gravity, where partition number is either the percent or fraction of a raw coal with a given specific gravity reporting to clean coal. The second assumption is that these curves are independent of the raw coal characteristics. At any specific gravity the system's partition number is solely a function of the partition numbers of each stage. For convenience, we use the following variables in describing this relationship: g = specific gravity, PS(g) = the partition number of the system at specific gravity g, P1(g) = the partition number of the primary stage at specific gravity g, P2(g) = the partition number of the secondary stage at a specific gravity g. For this system, we can show algebraically that the relationship is: [Ps (g) = PI (g) + (1- P1(g)) P2 (9)(1)] With this equation the system's partition curve can be determined if the partition curves of each stage are known. By definition, the partition numbers P1(g) and P2(g) must each be less than, or equal to, one. Placing these constraints on Eq. (1), we arrive at the relations: [Ps(g) % P1(g) and Ps(g) % P2(g)(2)] Thus, at any specific gravity the system's partition number is always greater than that of both stages. [Figure 2] shows typical partition curves. Curve I represents the primary stage, which consists of water-only cyclones operating at a separating gravity of 1.3. (The separating gravity or the

Jan 1, 1979

By J. B. Cheatham

Applications of the mathematical theory of plasticity promises to lead to the solution of many drilling and rock mechanics problems. ,Because of mathematical considerations, the inelastic behavior of rock has frequently been represented by a perfectly plastic model in conjunction with a yield criteria of the Coulomb or Mohr type. The totality of all stress states for which a solid ceases to behave elastically can be represented as a limit surface in stress space. ,Probing of such limit surfaces indicates details of strain hardening which are not provided by the standard triaxial testing procedure. Probing tests of the limit surfaces have been performed on Cardova Cream limestone to provide data for extending plasticity theory to cover situations in which consolidation and strain hardening are present. Test results indicate that this highly porous limestone undergoes a permanent volume decrease when it is subjected to hydrostatic pressures in excess of 3,500 psi. A virgin sample tested under a confining pressure of 1,500 psi has a yield strength of 1,700 psi; however, if the sample is subjected to a consolidation pressure of 5,000 psi, before testing at 1,500 psi, the yield strength is raised to 2,300 psi. Thus, both consolidation and strain hardening are important considerations in describing the mechanical behavior of this limestone. Tests conducted with the axis of the core having different orientations indicate that this rock is also anisotropic. Portions of the initial and subsequent limit surfaces are determined for samples loaded either perpendicular or parallel to the bedding planes. INTRODUCTION Previous experimental work in rock mechanics indicates that no mathematically tractable constitutive theory is inclusive enough to describe the mechanical behavior exhibited by all types of rocks under all conditions of stress and temperature. Indeed, the type of deformation encountered in a single type of rock is known to depend upon the stress and temperature conditions in the rock during deformation.l-5 Certain rocks and minerals, notably those minerals composed of ionic salts, have been shown to exhibit plastic deformation when tested under conditions of high confining pressure.2 Since the mathematical theory of plasticity provides simplifications over the theory of elasticity in certain types of problems, such as those in which limit analysis can be applied, it is of interest to know under what conditions plasticity theory may be applied to rock mechanics problems. The following factors determine the nature of the deformation a particular specimen will undergo: (1) the microscopic structure of the rock, i.e., the structure visible under an optical microscope, including number of phases, porosity, distribution of phases, (2) the mineralogical structure of the solid phases, (3) the conditions of stress and the rate of change of stress and (4) the temperature. Extensive experiments on Yule marble, Carthage marble, Solenhofen limestone and other calcerous rocks indicate that these relatively nonporous rocks deform plastically under certain conditions of loading, and creep under other conditions of loading.l-3 This study is concerned with the behavior of Cordova Cream limestone (Austin chalk) which is also composed almost entirely of calcite and thus has the same mineralogical composition but, because of a rather large porosity, it possesses a different microscopic structure. This investigation was undertaken to learn if Cordova Cream limestone deforms plastically despite the embrittling effect of the pore spaces, and to provide data which can be used to determine whether the mathematical theory of plasticity can describe the mechanical behavior of Cordova Cream limestone. A brief discussion is given of some of the underlying principles of plasticity theory concerning stress space and limit surfaces. Strain hardening is defined and a discussion is given of how the mechanical behavior of a material can depend on the previous history of the loading. This is followed by a description of the experimental procedures

By M. C. Flemings, A. A. Tzavaras

Analytical and experimental procedures are descrihed for quantitative determination of volume-fraction microporosity in metals by mimorarlio-graphg. The analysis assumes spherical pores but allows for the observed fact that pores or portions of pores smaller than a given thickness do not project an image on the film. Volume-fraction microporosity was measuured in cast low-alloy steel and shown to vary from 0.6 x 10-4 to 10 x 10 depending on mode of solidification and on whether the material was air- or vacuum-cast. The aim of work described herein was to develop a method of quantitative measurement of the volume fraction of microporosity in cast metals. Microporosity, present in nearly all castings and ingots, is always interdendritic in nature and may be large enough to be seen on a polished surface with the naked eye. More generally, it is smaller in size, and its detection more difficult. Three techniques might be considered for measurement of microporosity: 1) lineal analysis of a polished metallographic surface, 2) density measurement, or 3) quantitative microradiography. As will be shown later, volume fractions of porosity in directionally solidified castings may be the order of 10-4 to 10-3; such small amounts of porosity are difficult to identify and measure on a polished surface. Density measurement is also not a suitable technique (for alloys) when the volume fraction of porosity is low. This is because the density of the base metal is often not known within the necessary accuracy, due to microsegregation or partial phase transformations. In general, it is not possible to completely eliminate microsegregation without also altering the porosity. Because of the limitations of the foregoing techniques, quantitative microradiography was employed to determine amount of porosity. In recent work, microradiography has been found to be useful for qualitative measurement of porosity, including determination of size and shape of individual pores;&apos;)&apos; however, quantitative determination of amount of porosity has not been attempted. I) CALCULATION OF VOLUME-FRACTION POROSITY The following analysis employs procedures similar to those of Cahn and Nutting3 and Hilliard.4 Major simplifications are that pores are assumed to be spherical and to be present in sufficiently small quantities that no overlapping occurs. A basic difference from the previous treatments is that the limited "sensitivity" of microradiography is considered. Assume first 1) cavities are spherical and of uniform diameter, 2) cavities are randomly distributed, 3) volume-fraction porosity and specimen thickness to pore radius ratio are sufficiently small that overlapping of the projected images of cavities is negligible, 4) sensitivity is limited (i.e., holes or portions of holes of less than a given thickness do not project an image on the film). The sensitivity, is defined as: where 5 = minimum thickness of void from which an image will be projected on the film and t = specimen thickness. Now, from Fig. 1, it is possible to distinguish between two different types of holes. These are 1) holes whose centers lie within the region t/2 > X > 6/2, where X is distance from one surface of the specimen (positive into the specimen), and 2) holes whose centers lie within the region (6 — r) < X < 6/2, Holes whose centers lie outside X = (6 -r) project no image on the film. These two types of pores are considered separately below. Type 1 Pores. When 6 is greater than 0, the image of a spherical cavity of radius r appears on the film as a circle of radius a which is smaller than r. From geometry, Fig. 1: a2=r2-52/4 [2] where a = radius of image of pore on film and r = radius of pore. Since the pores are randomly distributed and spherical V =4/3pr2Nv [3] where V = volume fraction of voids and Nv = number of pores per unit volume and the number of images of Type 1 pores per unit area of film, N&apos;A, is

Jan 1, 1965

By L. A. Roe, E. A. Elevatorski

Wollastonite, named after William H. Wollaston, an English chemist, is a calcium metasilicate, CaSiO3; CaO: 48.30%, SiO2: 51.70%. It has a short history as an industrial mineral. The earliest production of wollastonite is reported to be from a deposit near Code Siding, located north of Randsburg, CA. At this locality small tonnages of wollastonite were quarried during 1933-34 and 1938-41, and processed into mineral wool. This operation was largely experimental and virtually no United States production was again reported until the 1950s when a large deposit near Willsboro, NY, was developed by the Cabot Corp. A processing plant was placed onstream in 1953, with nearly continuous production to date. It is currently operated by NYCO, a division of Processed Minerals, Inc. Since 1958, wollastonite deposits in the Little and Big Maria Mountains of Riverside County, and in the Panamint Range of Inyo County, both in California, have operated intermittently for production of both ornamental and commercial wollastonite. During 1980, the United States was the major producing country, furnishing about 75% of the world&apos;s output. Current production comes from Finland, Mexico, India, and Kenya. Small amounts have been shipped intermittently from the USSR, New Zealand, Republic of the Sudan, Republic of South Africa, and Namibia (South-West Africa). The principal use of wollastonite is in the manufacture of plastics. Other uses are for paints, ceramics, adhesives, fluxes, glazes, thermal insulation board, and refractory products. Mineralogy Pure wollastonite, CaSiO3, has the composition of 48.3% CaO and 51.7% SiO2. However, it is seldom found in the pure state due to the ease with which it takes into solution the metasilicates of manganese, magnesium, iron, and strontium. Predominantly, wollastonite occurs as a contact metamorphic deposit forming between limestones and igneous rocks. Commonly associated minerals are garnet, diopside, epidote, calcite, and quartz. It has a specific gravity of 2.8 to 3.0, and hardness of 4.5 to 5 on Mohs&apos; scale. When pure, it has a brilliant white color, but with impurities it may be grayish or brownish. Luster is vitreous to pearly. Melting point of wollastonite is about 1540ºC. Wollastonite occurs in coarse-bladed masses, rarely showing good crystal form. It is usually acicular or fibrous, even in the smallest of particles. The most unique property of crushed and ground wollastonite is its cleavage. Fragments of crushed wollastonite tend to be needle-shaped, imparting a high strength, and this property is the basis for many of its uses. The fiber lengths are commonly in the ratio of 7 or 8 to 1, length to diameter. The average diameter of wollastonite is 3.5 µm. Some crystals of wollastonite fluoresce under short-wave or longwave ultraviolet light, or both, with colors ranging from yellow-orange to pink-orange. Specimens may also show phosphorescence. The acicularity of wollastonite is a property of considerable importance to the marketplace. The plastics industry makes utilization of high aspect ratio grades of wollastonite (20:1) for reinforcing thermoplastics and thermoset polymer compounds. The naturally high pH of 9.9 (10% water slurry) is a prime property to the coatings industry. The coatings industry uses milled grades of wollastonite as a pH stabilizer in interior and exterior PVA and acrylic latex systems. Processed wollastonite can have a G.E. brightness of 90 to 93. Chemically, wollastonite is inert and this property makes it useful as a filler and reinforcing agent. There are two polymorphs of calcium silicates: wollastonite, a low temperature form, and pseudowollastonite, a high temperature

Jan 1, 1983

By Larry J. Tucker

INTRODUCTION The orebody at Pea Ridge strikes NE-SW and dips at about 800 to the SE. The ore and surrounding rocks are Pre-Cambrian in age. These rocks are usually very competent. Locally strong jointing has caused problems and toward the hangingwall dikes have caused particular problems during or after the open stoping phase of production. About 396 m of sediments lie over the ore. The shape and mineral content of the orebody has been fairly consistent so far with production areas lying from the 1375 to the 2475 levels. The deeper limits of the ore- body haven&apos;t been found yet. The ore is high grade magnetite with a strike length of over 600 m and a maximum width of over 180 m on some levels at the center of the orebody. Production at Pea Ridge for the first three years of operation came from two types of open stopes. The original stopes were recovered with horizontal production rings being drilled and blasted from vertical raises. Later a more typical sublevel stoping method using vertical production rings from one or more drilling sub- levels was used. There were a number of reasons for the selection of open stoping as the initial production method at Pea Ridge. First, it was judged that the ore was not of a block cave quality. After fifteen years of production, this assumption can be verified. There was some concern during the planning days for the potential for high flow rates of water into the mine from aquifers in the sediments lying over the orebody. Dewatering wells from surface and underground were put in as a precaution and nothing came of this potential problem. Open stoping and its resulting pillar structure put off this problem by supporting the sediments a few years longer. For a moderate investment in level development some production from open stoping was quickly realized. The first open stopes were brought in on what was then the bottom mine production level (2275), near the shaft pillar and situated on the footwall. This production practice of stoping on the heels of development from shaft pillar westward and from footwall to hangingwall continued through the open stope phases of production at Pea Ridge for many years. From the standpoint of mine efficiency and pillar stability, this would not seem to be a recommended mining practice. 74 570 w electric slushers, using 1.37 m dippers, pulled the open stope ore out into 18 ton bottom dump cars. Ordinarily five of these cars made up a train, which was pulled to the dump area by a trolley powered electric locomotive. Open stoping had the disadvantage that 100% of the recoverable portions of the orebody had to be developed in order to recover 33% of the ore. 67% of the ore was left behind in a large pillar system. For many years we had no planned method for recovery of these pillars. Mining methods are of course, limited by the existing technology. Advancing technology, coupled with the limited production capabilities of open stoping, low efficiencies and high relative cost early in the history of the operation recommended that other mining methods be looked at. In Sweden, considerable practical experience, coupled with a large amount of theoretical work had become available on sublevel caving. To Pea Ridge it offered the apparent benefits of: 1) The flexibility of diesel powered trackless haulage. 2) The efficiency of two machines to one operator on longhole drilling. 3) A multilevel production picture with a large number of drawpoints for ore recovery. 4) 100% recovery of the mineable portion of the orebody in one production pass, with no pillar recovery problems left behind. This paper outlines the experiments and experiences which Pea Ridge has gone through with sublevel caving from 1967 to this date.

Jan 1, 1981

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 R. W. Heckel, R. A. Tanzilli

Treatment of the finite, diffusion-controlled, two-phase, moving-interface homogenization problem has been carried out by numerical methods for a wide range of input parameters. Solutions have been obtained for planar, cylindrical, and spherical diffusion geometries. The effects of parameters such as the interdiffusion coefficients (concentration-independent) of- both phases, the concentrations (time-independent) of each phase at the interlace, mean composition, initial size of the dissolving- phase, and initial composition of- the two phases are considered. The results are displayed to show both the position of the interjace and the degree of homogenization (expressed in terms of concentration at the matrix phase symmetry boundary) as a function of diffusion time. Typical exaw~ples to which the results may be applied are discussed. A limited number of exact, closed-form solutions to the multiphase, diffusion-controlled, moving-interface problem exists in the literature, e.g., Refs. 1 to 7. The exact, closed-form solutions generally consider the outermost phases to be infinite in extent and also assume an interfacial velocity function. It is through use of such boundary conditions that closed-form solutions are readily obtained. An extensive literature also exists on numerical solutions to the mathematically analogous solidification problem, e.g., Refs. 8 to 12. In addition, a numerical solution involving mass transfer with spherical interfaces and infinite boundaries has recently been carried out.13 A practically important class of diffusion problems exists, however, in which the boundaries of each phase away from the interface are finite. For example the finite boundary condition arises when one sets up a model for analyzing the homogenization behavior of two-component powder compacts in two-phase systems or for studying the solution behavior of precipitates. Although exact, closed-form solutions have as yet not been found to this class of problems, solutions are possible using numerical methods. The purpose of this paper is twofold: 1) to describe a numerical method for treating the two-phase, diffusion-controlled, moving-interface problem involving finite boundaries with planar, cylindrical, and spherical interfaces; 2) to present the results of a parametric study which identifies the role played by each variable appearing in the analytical model. R. A. TANZlLLI,Student MemberAIME,onleaveofabsence from the General Electric Co., Re-entry Systems Department and A. W. HECKEL, Member AIME, are Graduate Student in Materials Engineering, and Professor, Department of Metallurgical Engineering, respectively, Drexel Institute of Technology, Philadelphia, Pa. Manuscript submitted May 6. 1968. IMD MATHEMATICAL MODEL The three geometries analyzed are shown in Fig. 1. For each configuration a symmetry element is shown, across whose boundaries the mass transfer is zero. Mean composition, C, and the characteristic solute-rich phase half-thickness (or radius), 1/2, are the independent variables while the characteristic symmetry element half-thickness, L/2 (measured from solute-rich phase center to the matrix phase plane of zero mass transfer), is a dependent variable. The symmetry element for a planar interface accurately represents a real diffusion situation involving the stacking of symmetry elements of constant size. Because cylindrical and spherical symmetry elements of constant size cannot be stacked to fill space, the meaning of the dimension, L/2, for these situations may be considered in terms of the following approximations: and where NA is the number of parallel cylindrical elements per unit area on a surface normal to the cylinder axes and Nv is the number of spherical elements per unit volume. The model consists of a solute-rich phase matrix phase composite of limited mutual solubility. The concentrations at the interface are assumed constant with time; no assumption of local equilibrium at the interface is necessary except for the case where the

Jan 1, 1969

By J. A. Musselman, P. F. Gnirk

A study was made of indexed penetrations by a single dull bit tooth under statically applied loads into rock subjected to confining pressures from atmospheric to 5,000 p.si and atmospheric pore pressure. Experimental results obtained with a 45" wedge tooth over the above range of pressures are presented for two limestones and a sandstone, a variely of indexing distances and two degrees of tooth dullness. The optimum distance between successive bit-tooth penetrations required for maximum rock damage and chip formation decreases substantially with increa.sing confining pressure above the brittle-to-ductile transition pressure of a particular rock. However, the distance remains approximately constant for a variation in confining pressure below the transition pressure. At a given confining pressure, the bit-tooth force required lor chip formation is constant for indexing distances greater than optimum, but generally decreases linearly with decreasing indexing distance for distances less than optimum. The chip-formation force data obtained at confining pressures for which the chip-generation mechanism is macroscopic-ally of a pseudoplastic nature compare favorably with previous theoretical results for indexed dull hit-tooth penetrntion into an idealized, rigid, perfectly plastic rock. INTRODUCTION To study the cutting action of a roller-cone bit, the interaction of a penetrating bit tooth with craters formed in a rock surfacc by the passage of previous bit teeth must be considered. Current knowledge of the basic mechanics of bit tooth-rock intcract:on under simulated borehole environmental conditions has been extensively reviewed.&apos;,&apos; Particularly of current interest are the effects of bit-tooth shape, distance between successive penetrations and differential pressure at the rock surface on the extent of the interaction. Pertinent to this paper are results of an experimental investigationb f the interaction between successive penetrations by sharp, wedge-shaped bit teeth. It was demonstrated that both bit-tooth angle and differential pressure influence the extent of rock damage between successive bit-tooth penetrations. Specifically, the optimum or minimum distance between successive penetrations required for maximum interact&apos;on or chip generation tends to decrease with decreasing bit-tooth angle and increasing differential pressure. In addition, at differential pressures on the order of 2,000 psi for four va- rieties of limestones and sandstones, the macroscopic mechanism of chip formation exhibits a transition from brittle to ductile. In this paper, experimental consideration is given to successive or indexed penetrat ons by wedge-shaped bit teeth. Since pore pressure in a rock sample is maintained at atmospheric pressure, the differential pressure at the fluid-rock interface is equal to the confining pressure. The degrees of bit-tooth dullness include teeth with flat and round apexes (Fig. 1). Of primary concern in this paper are the bit-tooth forces required for chip formation, the macroscopic mechanism of chip formation and the minimum distance between successive penetrations, i.e., optimum indexing distance required for maximum rock damage, as functions of differential pressure. Of further interest is a comparison of actual bit-tooth chip formation forces at elevated differential pressures with calculated forces from previous theoretical results4 for indexed dull bit tooth penetration into a rigid-plastic rock. EXPERIMENTAL APPARATUS AND PROCEDURE The experimental apparatus consists of a pressure vessel equipped with a piston through the side of the vessel (Fig. 2). A single bit tooth inserted in the lower end of the piston is forced at a desired rate into the rock when the ram pressure chamber is pressurized. Electrical instrumentation incorporated into the apparatus yields a graphical plot of force on the piston as a function of bit-tooth penetration or displacement into the rock. Since the piston and ram assembly are in force equilibrium for a constant confining pressure in the main vessel, the fluid volume in the vessel remains constant during bit-tooth penetration. Hence, the force resisting penetration of the tooth is independent of the force exerted on either end of the as-

By T. D. Mueller

Many investigators have used the response of the "dimensionless aquifer" to a unit pressure drop or a unit fluid-withdrawal volume to calculate the performance of an aquifer in supplying water influx to an oil reservoir. In the past, these response functions have been calculated with the aid of the Laplace transform. With the advent of ultra highspeed digital computers, it becomes practical to solve for the response functions with finite-difference techniques. The computer method also permits extension of the dimensionless-aquifer concept to include the nonhomogeneous aquifer wherein the permeability and other properties vary as a function of the space co-ordinates. This paper gives results of calculating the response functions for a series of nonhomogeneous aquifers. Response functions are presented for both linear and radial aquifers whose thickness, permeability-viscosity ratio and porosity-compressibility vary. These functions are new and should prove useful to the petroleum engineer in analyzing the behavior of nonbomogeneous aquifers. Results are presented in the form of charts that can be easily used by the field engineer. INTRODUCTION Aquifers which surround many oil and gas reservoirs have the ability to supply water influx to such reservoirs as oil and gas are withdrawn. This water influx, called natural water drive, provides one of the most effective driving mechanisms for the production of oil and gas. In producing a reservoir, therefore, it behooves one to make the maximum use of natural water drive. To achieve the maximum use, the reservoir engineer must be able to predict the performance of an aquifer under a variety of production schemes that may be proposed for the reservoir. Unfortunately, the physical properties which dictate aquifer behavior often are known only within limits. Seldom do wells penetrate the porous strata of the aquifer. Even when they do, quantitative information regarding porosity, permeability and water compressibility is seldom available. It is known, however, that the water efflux from most aquifer systems is governed by a single, relatively simple, linear, partial differential equation. Also, the general physical location of the aquifer boundaries often are known. A technique originall proposed by Hurst1 and van Everdingen and Hurst2 has been found useful in analyzing reservoirs in this situation. The idea was later expanded by van Everdingen, Timmerman and McMahan3 to include the mathematical technique of least-squares fitting. This latter method will be referred to as the VTM method. The basic assumptions of the VTM method include the following. 1. The location of the physical boundaries of the aquifer are known. 2. The flow conditions at these physical bound-daries are also known. 3. The aquifer is homogeneous; e.g., thickness, permeability, water compressibility and porosity are constant throughout. In the VTM method, a material balance is made on the fluids entering and leaving the reservoir. In the balance equations, the water-influx term is represented as the product of the water influx from an arbitrarily-selected, dimensionless aquifer system times an unknown conversion number. This balance can be formed as many times as there are data points in the history of the reservoir. Each time, the conversion number can be evaluated. If the reservoir engineer has picked the correct dimensionless aquifer to represent the real aquifer, the conversion number remains constant for all balances that have been made over the history period. If such a situation occurs, the reservoir engineer can then predict the performance of the reservoir for any type of production scheme by using the function associated with that particular dimensionless system and the derived conversion number. These functions will be referred to as the response function of the aquifer. If the aquifer is nonhomogeneous (e.g., if the porosity, permeability, thickness, porosity, or

By J. A. Elias, A. J. Heckler

A simplified method is described for the prepauatior~ of comnposite samples for determination of one quadrant of an X-ray pole figure of a sheet material by the Schulz reflection method. Other special conzposite sampling methods were used to investigate specific areas of the pole figures. Plotting nets for these sampling plans were drawn by a computer. Step-scanning techniques were used to obtain detailed information about specific regions of a pole figure. Satisfactory pole figures were obtained for large-grained (0.5 to 1 in. diam) sheet materials. COMPOSITE sampling techniques for preparing one quadrant of an X-ray pole figure have been reported by a number of authors.&apos;-5 All of these methods are similar in that they use a composite sample prepared by stacking and holding together coupons cut from a sheet material. The schulze reflection technique is applicable to the composite sample so that only minor corrections are needed for defocusing and absorption effects. Also, a relatively large area of a sheet is sampled through the use of multiple cross sectional pieces. In general, sampling methods necessary to obtain a composite-type specimen have been formidable enough to discourage their use in routine pole figure determinations. This paper presents a versatile sampling technique which is simpler than those previously described along with an accurate method for obtaining stereographic plotting nets for any orientation of a composite specimen surface. A step-scanning technique is used to obtain more quantitative determinations of specific regions in pole figures, and the composite technique is satisfactorily applied to determine pole figures of very large grained sheet materials. REVIEW OF PREVIOUS TECHNIQUES Mueller and ~nott&apos; first proposed a technique for obtaining one complete quadrant of an X-ray pole figure by using a sheet sample and six different composite samples in which the individual coupons were aligned and held together by two types of jigs. Data from the seven different samples were overlapped to obtain a single quadrant. ~eieran~ and Lopata and ~ula~ obtained a complete quadrant of a pole figure from a single sample prepared by forming a cube from coupons cut from a sheet. These coupons were bonded together using an epoxy cement, and a cut was made at 54.7 deg to the normal, transverse, and rolling directions of the sheet, Fig. I. The normal direction of the cut face projects into one quadrant of a stereographic projection at equal angles to the three principal orthogonal reference directions mentioned above. X-ray reflection data obtained along a spiral path in physical space could then be plotted on a stereographic net offset 54.7 deg from the normal, as shown in Fig. 2. A rotation about an axis representing the surface normal is designated as a, and a rotation about an axis lying in the plane of the surface is 4. From Fig. 2 it is apparent that a tilting angle of no more than 55 deg @ is necessary to obtain the complete quadrant. ~eber~ devised an improved sampling technique by cutting coupons 45 deg to the rolling direction of a sheet, stacking and bonding these with an epoxy cement at an angle of 54.7 deg to a reference plane, and finally grinding and lapping a flat specimen surface parallel to this plane. An important consideration of this technique and that of Lopata and ~ula~ is variation in the diffracting power of the specimen surface due to the bonding medium. While any composite-type sample will have some such variation due to boundaries between coupons, the problem can be minimized by elimination of the bonding medium and by avoiding prolonged etching times which could lead to significant inter lame llar attack of the coupon faces. EXPERIMENTAL TECHNIQUE We were working concurrently though independently on a technique similar to ~eber&apos;s,&apos; but which does not

Jan 1, 1968

By Sidney A. Swensrud

MUCH discussion has occurred as to the use or uses which might be made after the war of either or both .of the large diameter (24 inch and 20 inch) War Emergency pipe lines built by y the United States Government during the war, to move crude petroleum and refined products such as gasoline and light distillate fuels. Indeed, the disposition and future use of these lines constitutes, in the minds of many people, one of the most important problems of postwar readjustment in the Oil Industry. The two lines, running from points in Texas to the New York-New Jersey-Philadelphia area, were constructed only because of the inability during the war to continue normal ocean tanker movements from the Gulf Coast to the East Coast. Physically they could continue after the war to move crude or refined products, and it appears that lines of such large diameter as these, if operated at or near capacity, could move these liquids at costs per barrel which might be considered in the same general cost range as for tanker movements.1 If, therefore, we were to look forward at the end of the war to a shortage of tankers, it is probable that the continued use of these pipe lines for petroleum movements to the East Coast would be worth considering, as against building new tankers to take over the movement, provided there were no other economic or national defense reasons for restoring tanker fleets and no other better use for the pipe lines. Because of the extensive tanker construction program which has been and still is in progress, and because of our Navy&apos;s marked success in combating the submarine problem, it appears to be quite possible, however, that instead of a shortage of tankers at the end of the war, we may have a surplus instead. If this should turn out to be the case, then there would seem to be considerable doubt that these lines would continue to be used for moving petroleum or products to the East Coast at the expense of tying up tankers. It would hardly seem probable that the owner of any such tankers would ship via one of the large lines at the expense of tying up one of its own tankers, unless the pipe line tariff rate charged was lower than the out-of-pocket cost of operating a tanker as against the out-of-pocket cost of keeping it tied up. And considering that a substantial part of the cost of operating a tanker goes on even when it is tied up, it is very difficult to figure out how any pipeline tariff could be established that would be low enough to compete with individually owned tankers that otherwise would be idle. Certainly a pipe line tariff rate high enough to cover direct costs plus even a low rate of depreciation and return on the investment would be far above the out-of-pocket tanker operating costs. In view of these considerations, it has seemed desirable to examine with some care some of the various alternatives that have been suggested as to other possible uses to which these lines might be put. One of these is for transportation of natural gas to the New York-New Jersey Philadelphia area. This Eastern area, located at the destination of these lines, contains the heaviest concentration of population in the United States, and is not now served with natural gas. At the other end of the lines lie the greatest discovered natural gas reserves in the United States, if not in the world, and many

Jan 1, 1944

By Ralph W. M. Lai, D. W. Fuerstenau

This paper presents the results of a study of the principles which control the distribution of ultrafine particles between an oil phase and an aqueous phase. Alkyl sulfonates were used to control the wettability of fine alumina particles and to stabilize the water-isooctane emulsion system in the extraction process. The effects of sulfonate chain length, sulfonate adsorption, pH, and oil-water ratio on the recovery of solids in the oil phase were investigated. Higher extraction was found under conditions which yield higher hydrophobicity and large interfacial area of the oil-water emulsion. Today there is considerable interest in developing methods for the recovery of fine particles from disseminated ores. One such method might be the extraction of fine particulate solids from aqueous suspensions into an oil phase, and one example of a successful industrial process based on oil-water distribution phenomena is the emulsion flotation process used in concentrating manganese ore.&apos; The oil agglomeration process of Farnand, Smith, and Puddington has been suggested for concentrating tin ore. Patents 3-4 for oil flotation processes extend back to that of William Haynes in 1860. The basis of all these processes is directly related to the stabilization of emulsions with solid particles. The first extensive study of oil-water emulsions stabilized by finely divided solids was made by Pickering. Later, Briggs,6 weston, 7 Chessman and King,8 and Leja and Schulman9 experimentally demonstrated how solid particles control water-oil and oil-water systems. Hildebrand10 suggested that fine particles can stabilize an emulsion because the fine particles are partially wettable by both oil and water, that is, that the contact angle must be finite. Takakuwa and Takamori 11 recently applied emulsion inversion to the study of collector adsorption in flotation systems. At the Royal School of Mines, Shergold and Mellgren12 are currently carrying out a comprehensive study of the extraction of fine particles into an organic liquid phase. In the present study of the liquid extraction or oil flotation of ultrafine particles, the recovery of hydro-philic alumina particles of 0.1-µ average diameter was determined after separating an isooctane-water mixture (emulsion) into its two phases. Alkyl sulfonates of different chain lengths were used to change the wettability of the alumina and to effect emulsifi-cation of the oil-water mixture. The effects of the oil-water volume ratio, sulfonate adsorption, and pH on the recovery of the fine particles were investigated. Since the surface charge on an oxide such as alumina is determined by the adsorption or interaction of H+ and OH- with the surface, pH control is extremely important. Furthermore, since the zero point of charge (zpc) of alumina13-14 occurs at pH 9 to 9.4, major adsorption effects with anionic collectors in these systems will occur below pH 9. The results of the foregoing investigation are discussed in terms of the recovery of the fine particles, the oil-water ratio, contact angles, surface tension, the adsorption of sulfonate onto the solid surface, and the emulsifica-tion of oil. MATERIALS AND EXPERIMENTAL METHODS Materials: The particulate solid used in all the experiments was Linde "A" alumina, a high-purity, readily available form of a-alumina. The surface area of the alumina, as evaluated from measurements of stearic acid adsorption from benzene, was 14.5 sq m per gm, based on a value of 20.5 sq Å for the cross-sectional area of stearic acid.13 All aqueous solutions were prepared from triply distilled conductivity water. Other materials included high-purity crystals of sodium alkyl sulfonates, isooctane, and, for pH adjustments, reagent grade hydrochloric acid and sodium hydroxide-Because sulfonic acids are strong acids, they are

Jan 1, 1969

By C. van der Poel, R. L. Chuoke P. van Meurs

When an initially planar interface between two im-ttitcihle liquids is displaced at constant rate, U, nor-mat to the front, instability will occur for all rates greater than a critical rate. U, given by provided the Fourier decomposition of the spatial pertmrbcrtiott or deformation of the moving displacetment front contains modes with wavelengths. A, greater than a critical wavelength, given by tit there expressions, p, p and k, with the subscripts I and 2 distinguishing the two 1iquids under consideration, are density, viscosity and effective pemeability coefficients respectively; g is the absolute value of the acceleration due to gravity; and cos (zz&apos;) is the direction cosit7e between the vertical carlesian coordinate z&apos; (positive upward) and the z coordinate norma1 to the initially plane macroscopic interface taken positive in the direction from Liquid I to Liquid 2. U is the average volumetric velocity (injected volume of liquid pet. unit time per unit of area normal to Z) and is positive for Liquid I displacing Liquid 2, negative for the reversc displacement and r2 is an effective inter-fucial tension (for displacements between parrallel plates, t0= s the ordinary bulk interfacial tension). further, there is a wavelength of maximum instability given by For natural perturbations, this wavelength will preclotninate and characterize the form of the instability as a quasi-sinusoidal deformation, i.e., viscous fitzgers of peak-to-peak separatiotr, A,,. Comparison of experitient with theory rei&apos;eals tltat these considerations are recsonably valid for displacements in hoth parallel plate chantzels and nnconsoli-clutt&apos;d glass ptnvder packs. These results are. tij particular interest for interpretation of data from laboratory experiments with regard to fielil applica/iotz.s. It is quite possible for A,, although st7zall relative to the large lateral dimension of the reservoir, to rxceed the greatr.st lateral extent of the lahoratory sample or model. If it does, frontal displace- ment occurs in the laboratory experiments, and oil INTRODUCTION Thc puposes of this paper are to present theoretical and experimental evidence for occurrence of macroscopic instabilities in displacement of one viscous fluid by another immiscible with it through a uniform porous medium and to compare available experimental data with some predictions of a theory of instability developed by the first author. The instabilities are referred to as macroscopic in the sense that spatially quasi-sinusoidal, growing fingers of the displacing liquid are formed, the width and peak-to-peak separation (wavelength) of which is large relative to a characteristic length of the particular permeable medium such as grain size. Visual models of two kinds have been used to obtain observations: displacement of oil by water-glycerine solutions through the flow channel formed by closely spaced parallel plates and displacement of oil by water with and without initial interstitial water through unconsolidated glass powder packs, employing the technique of matching indices of refraction. In all cases we have observed macroscopic instabilities or fingers under conditions predicted by the theory to be favorable for their occurrence. The phenomenon discussed here is not the production of streamers due to gross inhomogeneities such as permeability stratification of the porous medium. It is our object to show, on the contrary, that a quantitative prediction of finger spacing is possible in a porous medium known to be macroscopically homogeneous and isotropic throughout. The importance of the phenomenon in its influence on the configuration of oil and water with respect to oil production behavior was noted earlier by Engel-berts and Klinkenberg&apos; who coined the term "viscous fingering". THEORY Necessary and Sufficient Conditions ior Instability anc Initial Kinematics There are several levels of increasing complexity in the theoretical description of instabiIity of fluid displacements in permeable media. Of these, the simplest description, ;tdapted to low permeability systems, is selected for presentation.&apos; More inclusivc tlescriptions are reserved for separate publication.

By N. H. Harrison, J. K. Rodgers, S. Regier

This paper presents comparisons of data obtained from a laboratory reservoir study and from a calculated behavior prediction with the actual production history of a condensate reservoir. A small non-commercial discovery was depleted under closely controlled conditions and the well fluids were sampled at frequent intervals. Data on the reservoir and production variables were accumulated on a fixed schedule. A laboratory reservoir study war made using the initial well fluid samples as charging stock. The production procedures and operating conditions were held constant throughout the study wherever possible and in general paralleled the field work. The well fluid compositions and the cumulative recoveries ar a function of the reservoir pressure were also calculated using conventional flash vaporization procedures and equilibrium constants. Comparisons based on the composition of the well fluid show good agreement, the laboratory study agreeing within experimental accuracy with the field work and the calculated data comparing equally well. The gas-oil ratios are also in good agreement, but with somewhat greater deviations at the higher pressures. In the overall picture, it is believed that a model st,~tiy can predict within experimental accuracy the production history of a condensate reservoir. Better equilibrium constants for the heavier hydrocarbons are needed in order to attain improved composition accuracy by calculation. INTRODUCTION In Aug., 1955, a gas condensate well was completed in San Juan County, Utah, that was initially thought capable of good commercial production. These conclusions were derived principally from core data and electric logs, which indicated good permeability, porosity and gas content. However, after the usual series of potential tests it was found that the reservoir pressure had declined some 22 per cent, and it was obvious that the zone tapped was but a small pocket or trap. It became apparent that, with a controlled depletion of a small reservoir, a unique opportunity was available to compare laboratory and calculated studies with an actual field depletion and to further the present knowledge of condensate reservoirs. FIELD WORK The Coalbed Canyon Well No. 1 was conventionally completed in the Paradox limestone formation to a total depth of 5,912 ft. The producing zone from 5,762 to 5,806 ft was perforated with four jet shots per ft. The wellhead and field equipment were also conventional, the major items consisting of a two-pass indirect fired line heater, a high- and a low-pressure separator with the necessary controls and accessories, gas meters, back-pressure regulators, flare stacks and condensate stock tanks. The initial testing of this well con-sisted of a series of flow potential and pressure build-up tests during which some 30 MMcf of gas was produced. The reservoir pressure declined from an estimated 2,300 to 1,782 psig during this period, from which it was concluded that the reservoir was very small. In order to approach steady-state conditions in the reservoir and so provide optimum conditions for making comparisons, the field depletion was programmed to approach, if possible, constant production conditions. Bi-hourly readings were taken of the tubing pressure, the pressure and temperature of the separators, oil and gas rates, and other pertinent operating data. The gas rate, as indicated by the orifice meter, was held constant by the adjustment of the choke in the line heater. The temperature of the first stage separator was held constant by adjustment of the line heater jacket temperature. Practical considerations of production made the maintenance oi a constant gas rate impossible. The test started with a gas rate of 4 MMcf/D and a separator pressure of 250 psig. This rate was maintained until the choke was fully opened. The gas rate then declined with the falling tubing pressure and production was continued until the rate was about 2

By S. W. Kennedy, M. Cohen, L. D. Calvert

A high-temperature X-ray diffraction method has been used to study the composition and the kinetics of formation of oxide scales at 800 °C on iron and pure iron-nickel alloys containing 25.6, 75, and 84 pct Ni. This technique has been combined with chemical analysis, and metallographic examination, and determination of weight-gain time curves. The 75 and 84 pct Ni alloys form a three-layer scale consisting of porous NiO next to the metal, a spinel layer, and a -Fe203 on the outside. The 25 pct Ni alloy produces a scale of spinel and Fe,O,, in which, however, traces of FeO and NiO are observed near to the metal. On all the nickel alloys the Fe20, is orieqted with the (107) plane parallel to the surfaces. The exposed surface of the Fe20, is in the form of needles whose axis is perpendicular to the specimen surface. Although the alloys show nonuniform oxidation rates as determined from weight-gain measurements, the X-ray data indicate that the total scale and individml scale layers grow uniformly and parabolically. Deep subscale formation is observed in the 25 pct Ni alloy only. This appeared to be associated with cracking of the scale and penetration of the oxygen along the grain boundaries. ALTHOUGH film formation has been investigated by high-temperature electron-diffraction techniques, for example by Gulbransenl and Hickman,&apos; there have apparently been no investigations of scaling by methods permitting study of the scale during formation at the reaction temperature. Methods involving quenching of scales are based on the working assumption that no significant changes take place on cooling. In general, however, compounds may disproportionate, and the different thermal expansions of constituents may cause disruption of the scale on cooling. In particular, some scales cannot be quenched without spalling. These difficulties are avoided by X-ray diffraction examination at the temperature of formation. High-temperature diffraction procedures also provide a method of determining the kinetics of formation of each phase separately in a multicom-ponent scale. The scaling of binary iron-nickel alloys Containing up to 30 pct Ni at a temperature of 800°C has been studied by Benard and Moreau using cooled specimen~,~ whilst Foley, Druck, and Fryxel14 made an electron-diffraction study of cooled oxide films formed on an iron 42 pci Ni alloy. In the present work the scaling of three iron-nickel alloys containing 25.6, 75, and 84 pct by weight of nickel, at 800°C in oxygen at atmospheric pressure was investigated by X-ray diffraction techniques using a Norelco Geiger-counter diffrac-tometer. The oxidation of iron was investigated for comparison. In addition normal metallographic examinations and chemical and X-ray analyses were made, and weight-gain time curves obtained. EXPERIMENTAL Materials Used—The analyses of the alloys were as follows: Percent by Weight Nickel 25.6 pct 75.4 pct 83.2 pct Aluminum 0.004 None found 0.002 Cobalt — 0.05 0.05 Chromium — — — Silicon — 0.01 None found Carbon All less than 0.05 pct The chemical analysis of the iron was: Al— 0.002 i 0.001 pct, C—0.007 pct, Cr less than 0.001 pct, Cu—0.002 ± 0.001 pct, Mn—0.003 ± 0.001 pct, N—0.0002 pct, 0-0.013 pct, P—less than 0.0002 pct, Pb—not detected, S—0.0017 pct, Si—0.0095 pct, Sn—none detected. This iron was the same as that used to make the alloys and was obtained from the same source-Vacuum Metals Corp. One specimen of pure iron from The British Iron and Steel Research Association was used for preliminary investigation. Specimens were machined to shape from cast rod, polished (the final polish being with 600 Sic paper), degreased in trichlorethylene, washed, and etched. The etching and neutralizing procedure adopted by

Jan 1, 1960

By M. R. J. Wyllie, Michael A. Torcaso

With the continued deep drilling of today, increasing numbers of high pressure and high temperature gas-condensate reservoirs are being discovered. Correspondingly, the ranges of properties of gas-condensate reservoir fluids are also being increased. The literature1*&apos;." contains several excellent reviews of the property ranges in presently known condensate reservoirs. When one examines these property ranges, it becomes apparent that the need for reliable data is great if the engineer is to evaluate properly or appraise a gas-condensate reservoir. It is the purpose of this paper to present phase-behavior data that will point out some of the important variables that influence recovery from gas-condensate reservoirs produced by pressure depletion. All these variables may be resolved into two major categories: the hydroc:arbon composition of the reservoir fluid and the physical properties of the reservoir. But the variables which we will consider are those which are usually determined from well tests shortly after discovery and are, consequently, readily available to the engineer. They are as follows: (1) the producing separator gas-stock-tank liquid ratio, (2) stock-tank liquid gravity, and (3) reservoir temperature. Laboratory depletion studies have been performed on several reservoir fluids to illustrate the effect of changes in the cronsidered variables upon ultimate recoveries. The deviation factors, Z, of a series of reservoir fluids were determined experimentally to illustrate their order of magnitude and their effect upon the calculation of in-place reserves existing at high pressurcs and temperatures. The procedures followed in conjunction with the performance of these studies and their associated calculations are illustrated in Fig. 1. For the purpose of this paper, the depletion recoveries are expressed only in terms of pcr cent of initial stock-tank liquid, with the abandonment pressure assumed to be 500 psig. VARIATION OF RECOVERY WITH INITIAL PRODUCING SEPARATOR GAS-STOCK-TANK LIQUID RATIO The producing separator gas-stock-tank liquid ratio has long been one of the principal criteria in assumptions and generalized correlations&apos; used for the estimation of ultimate recovery from gas-condensate reservoirs. Fig. 2 presents a plot of the cumulative recovery of stock-tank liquid at abandonment as a function of the initial producing separator gas-stock-tank liquid ratio. The graphical illustration represents the results of 44 laboratory depletion studies performed on actual gas-condensate fluids. These data are presented not only to illustrate the general trend of ultimate recovery increase with initial producing ratio increase but to emphasize the considerable variation in possible recovery at any given producing ratio. Table 1 presents the properties of two actual reservoir fluids chosen to illustrate the approximate extremities of producing ratios for gas-condensate systems. The fluids designated as A and B exhibited initial producing ratios of 3,400 and 114,000 scf/bbl, respectively. Fig. 3 presents a graphical comparison of the cumulative stock-tank liquid recoveries, calculated from the well-stream compositions and produced volumes for each fluid during depletion. It is to be noted that Fluid A exhibited a recovery of only 16.5 per cent when produced to an assumed abandonment pressure of 500 psig while the corresponding recovery for Fluid B was calculated to be 62.2 per cent. This difference in recovery, which exemplifies the general trend of recovery increase with initial producing ratio increase, is readily explained by consid-

By John P. Huger

me phase diagram for the system Ag-Si has been established by thermal analysis. The diagram is a simple eutectic type with a eutectic at 3.0 pct Si and 840°C. From the liquidus curve and with the aid of regular solution theory, the thermodynamic properties of the liquid Ag-Si system at 1420 °C have been determined. Silicon exhibits positive deviations from Raoultain behavior at all compositions, while silver exhibits both positive and negative deviations. MEASUREMENTS of silicon activities in liquid iron alloys by means of distribution studies between the immiscible liquids silver and iron have been reported by Chipman, Fulton, Gokcen, and caskey.&apos; The applicability of this technique is dependent, however, on the accuracy of the silicon activity values in the Ag-Si alloys, particularly in the dilute silicon region. Darken and Gurry2 have shown that log ?si/(1-xSi)2 is approximately linear with composition, but uncertainty exists in the extrapolation to dilute solutions. Their calculations are based on the old Ag-Si phase diagram of Arrivaut3 and any variation in the position of the liquidus curves will strongly influence the extrapolation. The purpose of this study is to redetermine the Ag-Si phase diagram and calculate the associated thermodynamic properties. EXPERIMENTAL The Ag-Si phase diagram was determined by thermal analysis of liquid alloys. The alloys were contained in a 35 mm alumina crucible under a protective stream of purified argon. The crucible was placed inside a 64 mm alumina tube fitted with a gas-tight, water-cooled, copper head. The tube and assembly were heated in a manually controlled globar furnace. The temperature was measured with a Pt, Pt + 10 pct Rh thermocouple in a high-temperature porcelain protection tube. The thermocouple calibration4 was checked against the melting point of pure silver at the start and during the course of the experiments, and against both gold and silver at the completion of the experiments. The alloys were prepared from granular silver analyzing 99.95 pct Ag and powdered silicon analyzing 99.88 pct Si and 0.026 pct Fe. Approximately 150 g of alloy were used for each run. The composition of the alloy was determined by withdrawing samples of the liquid alloy with a vycor tube-aspirator bulb device. Analysis of solidified samples agreed, for alloys containing less than 23 pct Si, within ±0.3 pct Si of the calculated composition. Duplicate samples showed a variance within 50.05 pct Si of the mean analysis. The alloy was heated to 100&apos; to 20pC above the estimated liquidus temperature and then the furnace was adjusted to give an average cooling rate of 2" to 7°C per min. A cooling curve was obtained by taking potentiometer readings at 15-sec intervals. Several analyses were made by recording the thermocouple output on a high-speed, Stip-chart recorder, with the exception of run NO. 11, a minimum of two thermal analyses were made for each alloy and the cooling curves then obtained duplicated the liquidus and eutectic temperatures to about 1°C. RESULTS The results of the thermal analysis are given in Table 1, and are represented as the Ag-Si phase dlagram, Fig. 1. The eutectic composition was found

Jan 1, 1963

By J. W. Marx, R. H. Langenheim

The injection of heat-bearing fluids may offer a wider application to secondary and tertiary recovery from conventional oil reservoirs than inderground combustion since the process is more easily controlled and the reservoir requirements, in general, are less critical. In all cares, conduction heat losses to overburden and un-derburden will impose an economic limit upon the size of the area which can be swept out from any one injection point for any given set of conditions. The present report describes a method for c.rtin~ating thermal invasion rates, cumulative heated area, and theoretical economic limits for sustained hot-fluid injection at a constant rate into an idealized reservoir. Full allow-ance is made for non-productive reservoir heat losses. Required solutions can be obtained directly from the numerical tables given here, as soon as all operating conditions are specified, withou extensive mathematica1 ~nanipulation on the part of the user. The injection of a heat-bearing fluid into a reservoir is frequently proposed as a means of secondary or tertiary oil production. This form of thermal recovery may offer a wider application to conventional oil reservoirs than underground combustion since the process is more easily controlled and the reservoir requirements, in general, are less critical. Because of its large gross heat capacity, steam appears to be the most efficient heat injection medium, although mixtures of steam and other gases, hot water, hot oil and hot non-condensable gases also have been employed. Whatever the nature of the injection medium, economic evaluation of such thermal recovery processes will depend upon thermal invasion rates, or rates of productive reservoir heating, at any given time after the start of heat injection. In all cases, conduction heat losses to overburden and underburden will impose an economic limit upon the size of the area which can be swept out from any one injection point, for a given set of reservoir conditions, at any given heat injection rate. The present report describes a method for estimating thermal invasion rates, cumulative heated area and theoretical economic limits for sustained heat injection at a constant rate into an idealized reservoir, making full allowance for non-productive reservoir heat losses. Required solutions can be obtained directly from the nu- merical tables given here, as soon as all operating conditions are specified, without extensive mathematical manipulation. It should be noted that ordinary water flooding may be used to move an annular-heated zone outward, even though the economic areal limit for straight heat injection has been reached. Thus well spacing need not be confined to the areal limits as determined here. However, the present calculations do not cover such supplementary heat scavenging programs. The operating element consists of a radial flow system, concentric about the point of injection, with the idealized step function temperature profile illustrated by the dashed line in Fig. 1. The solid curve shown in Fig. 1 gives a qualitative picture of the true radial temperature distribution. It will be shown later that the expressions given below, which are derived from the simplifying step idealization (dashed line in Fig. I), also can be interpreted in terms of the more realistic temperature distribution shown in Fig. 2. For this case, the thermal invasion radius is simply redefined as the distance from injection well to the mid-point of the temperature distribution.