Search Documents

By R. G. Grieves, G. Thodos

A method is presented for the accurate calculation of the cricondentherm and cricondenbar pressures of multicomponent hydrocarbon mixtures of known composition. The mixtures may contain six and quite possibly any number of components including paraffins, isoparaffins, olefins, acetylenes, naphthenes and aromatics. The approach is similur to that used for calculating critical pressures and cricondentherm and cricondenbar temperatures. The critical pressure, the normal boiling point and the approximate vapor pressure behavior of each component are all that art required. A stepwise culculation procedure is necessary for mixtures containing more than two components. From an analysis of 22 binary systems and 118 mixtures, the average deviation of calculated cricondentherm pressures from reported values is 2.4 per cent. For nine multicomponent mixtures the average deviation is 1 per cent. Considering 19 binary systems and 108 mixtures, the average deviation of calculated cricondenbar pressures from reported values is 1.7 per cent. For 15 multi-component mixtures, the average deviation is 2.2 per cent. INTRODUCTION A knowledge of the phase behavior in the critical region of multicomponent hydrocarbon mixtures is of value both in industrial processing operations and for the optimum operation of gas condensate reservoirs. Accurate methods of calculating the critical temperatures and critical pressures4,6,7,8,9,23 and cricondentherm and cricondenbar temperatures6p10 of multicomponent hydrocarbon mixtures are available in the literature. If the cricondentherm and cricondenbar pressures could be calculated with equal accuracy, the entire phase diagram of a multicomponent hydrocarbon mixture could be well approximated. The work of Etter and Kay6 is limited to systems containing the normal paraffins and has not been tested on systems containing a heavier component than heptane. In addition, the development of their multicomponent equations is based upon a limited number of mixtures containing from three to six components. Silverman and Thodos have considered systems containing both paraffinic and non-paraffinic hydrocarbons, but their correlation is limited to binary systems and is highly inaccurate for methane systems. Eilerts has done extensive work on the cricondenbar pressure; he has produced an excellent correlation for binary systems. However, his procedure for multicomponent mixtures is chiefly useful for highly complex mixtures requiring a knowledge of the vapor-liquid equilibrium behavior of the mixtures; he has not considered the cricondentherm pressure. The objective of this study was the development of an accurate and rapid method for the calculation of the cricondentherm and cricondenbar pressures of multicomponent mixtures containing all types of hydrocarbons and having a wide volatility range. The approach that was adopted is similar to that used by Grieves and Thodos for cricondentherm and cricondenbar temperaturesl0 and for critical pressures.9 However, methane systems had to be considered separately and a modified stepwise calculation procedure was utilized for the cricondertherm pressure. The correlations were developed in a manner similar to those for critical pressures and cricondentherm and cricondenbar temperatures. Based upon binary data reported in the literature it was observed that the ratios of cricondentherm and cricondenbar pressure to the pseudocritical pressure (molar average), pt/ppc and pp/ppc, respectively, in two-component systems depended upon the mole fraction of the low-boiling component and upon the diversity in properties of the two components. A dimensionless boiling-point parameter Tb,/Tb was chosen to represent the diversity in properties of the components. For a binary system, Ti,'is the molar average of the normal boiling points of the two components involved. For a multicomponent system it is either the molar average of the normal boiling points of all of the components involved or the molar average of the normal boiling point of the pure low-boiling component and of the atmospheric boiling point of the low-boiling-free mixture of

Jan 1, 1965

By M. Prats

The effect of a sand-filled vertical fracture of limited radial extent and finite capacity (fracture capacity is the product of the permeability and width of the fracture) on the flow behavior of a cylindrical reservoir producing an incompressible fluid through a centrally located well has been investigated mathematically. The shape of the lines of equal pressure near the fracture is essentially independent of the size of the reservoir, provided that the field radius is of the order of the fracture length or larger. For reasonable values of production rates and of fluid, reservoir and fracture properties, the total pressure drop between the end of the fracture and the well is gelrerally negligible compared with the pressure drop in the reservoir. With regard to production response, the effect of vertical fractures can be represented by the production response of an equivalent or effective well radius. For a high-capacity fracture, the effective well radius is a quarter of the total fracture length, decreasing with the fracture capacity. When invasion effects are simulated by decreasing the width of the damaged zone with distance from the well, the effect of formation damage around a fracture on the production response is not so seriolls as indicated by the literature. This suggests that frac fluids with a conventional filter-loss response are better than high-spurt-loss frac fluids, provided the effective permeability of the damaged zone is the same. INTRODUCTION This paper considers the effect of the fracture capacity, as well as the formation damage which can result from fracture treatments, on the productive capacity of vertically fractured wells. Other publications, notably those of van Poollen,l,2 consider these same effects. In addition to providing more general results for vertical fractures than are available from the liter- ature, the present paper gives the equivalent well radius of fractures having different lengths and capacities and, also, includes pressure distributions in and around the fractures. The effect of a damaged zone around a fracture on the production response was not found to be so great as that reported by van Poollen. 2 This difference probably stems from the fact that we consider a damaged zone which is widest (but is still small) near the well and thins out toward the extremities of the fracture, whereas van Poollen considers a damaged zone having a uniform width for the entire fracture length. Simple, but adequate, equations which describe the effect of these variables on production response are presented (in Appendixes A and B). Thus, results can easily be extended to values of the variables not specifically considered here. IDEALIZATION AND DESCRIPTION OF THE FRACTURED SYSTEM It is assumed that a horizontal oil-producing layer of constant thickness and of uniform porosity and permeability is bounded above and below by impermeable strata. The liquid is incompressible. It is assumed that, after the well is fractured, the cylindrical outer "boundary" is at a uniform potential, provided that it is not too near the fracture. The fracture system is represented by a single, plane, vertical fracture of limited radial extent, bounded by the impermeable matrix above and below the producing layer (reservoir).3 Fig. 1 indicates the general three-dimensional geometry of the fractured reservoir described above. If gravity effects are neglected, the flaw behavior in the reservoir will be independent of the vertical position ill the oil sand. This means that the flow behavior in the fractured reservoir can be described by the two-dimensional flow behavior shown in the horizontal cross section of the reservoir, Fig. 2. PARAMETERS The effect of a fracture on the pressure distribution in the fracture itself and in an otherwise undamaged (no skin) reservoir indicates that the pressure distribution is a function of three parameters.

By J. J. Gilman

Zinc monocrystals, when compressed nearly parallel to their basal planes (within 20), deform in a new way; and certain deformation markings, called "-bands," are typical of this new mode. Characteristics of these bands are described. Mechanical behavior of low angle zinc monocrystals was investigated in some detail. Recrystallization did not occur upon annealing of banded specimens that were free of deformation twins. -Bands do not seem to be the result of nonbasal slip, but rather, a result of a microbuckling mechanism. u-Bands differ qualitatively and quantitatively from ortho kink bands. A dislocation model, consistent with some of the features of u-bands, is described. THE plastic deformation of zinc monocrystals with slip planes (0001) nearly parallel to the stress axis was studied in this investigation. In a previous investigation, it was shown that some low angle zinc crystals (those with slip planes making angles of more than about 2.5+1° with the compression axis), deform by ortho kinking while others (those with slip planes making angles of less than 2.5+1° with the compression axis) deform in a new way.' This new mode of deformation will be described in some detail here. Some of its striking features are: 1—it operates only for compressive stresses, 2—its temperature dependence is anomalous, and 3—it exhibits a new type of "yield-point" behavior. The mode of deformation to be discussed is not new in the sense that no previous evidence of its existence had been reported. The markings that it creates on the cleavage surfaces of zinc crystals have been described by several investigators.'" However, this is the first time that it has been made to operate under controlled conditions so that its response to mechanical variables, temperature, and aging treatments could be investigated and its geometry completely described. Apparently. Jillson' as well as Holden3 and Pratt" believed that the small bands they observed were simply small versions of ortho kink bands. In fact, two of these authors3,6 "alled the markings "micro kink bands." In spite of the precedence, it is felt that this term is not entirely suitable because these bands do not have the same genesis as kink bands, and they are qualitatively different in structure. As will be indicated by later discussion, the bands seem to be related to the tensile deformation bands observed in aluminum. The present author has chosen to call the bands "micro deformation bands," or simply ",ubands," and the deformation process, ''p-banding." Experimental Methods Crystals were grown from 99.999+ pct Zn obtained from the New Jersey Zinc Co. The details of the technique are given elsewhere." It should suffice to say that the crystals were seeded for orientation control, and precision pyrex tubes were used for molding them. About 40 crystals were used, yielding some 130 specimens. Fig. 1 illustrates the typical form of the specimens and shows how the orientation angles are defined. The crystal specimens were usually 1 3/4 in. long and about 150 mil in diameter. The angle between the slip or cleavage plane (0001) and the specimen axis is designated by xo, while a is the angle between the maximum slope of the (0001) plane and the nearest close-packed direction [1210], vary ing from 0°to 30°. The crystals were cut into specimens with a jeweler's saw and the ends were squared by filing (specimens held in V-block). Crystals of this orientation are stiff toward bending, so this procedure introduced only deformation localized at the ends. Then the specimens were polished either chemically in 50 pct HNO³ or electrolytically in 20 pct chromic acid. The mechanical tests were made with an Instron testing machine. This machine weighs by means of electric strain gages on elastic bars and is quite stiff, since the full-load deflection of its load cells is at most 0.002 in. Most of the testing was done at a

Jan 1, 1956

Salt Lake City was founded July 24, 1847, by Mormons under the leadership of Brigham Young. It had a population of 1 18.1 10, according to the United States census of 1920, and of 15 1,968, according to a civic survey, January 1, 1925. The Salt Lake metropolitan district is estimated to have 180.079 inhabitants. The city's area is 5 1.8 square miles. Its altitude is 4,354 feet. It has 514 miles of streets, of which 91 miles are paved. The streets in the main portion of the city are 132 feet wide and the blocks are 660 feet long. It has 144 miles of street railway and 2,200 electric street lights. It has fifteen parks and playgrounds, forty public schools and six private schools. It has eighty-nine church edifices, of which forty-nine are devoted to denominations other than the Mormon church.

Jan 1, 1925

By B. Trafton, H. L. Hartman

An earlier paper1 introduced a rapid solution for mine ventilation network problems, employing the digital computer. In this paper, refinements and additions to the computer program which expand its scope of application and a solution to a ventilation network having a coal-mine configuration are presented. The revised program-written in FORTRAN for the IBM 7074 computer in Penn State's Computation Center-is applicable to any mine ventilation network employing natural splitting and having its fan located externally. It accomplishes the following: 1) Distributes an assumed total quantity of air flow among all airways of the network, performing iterations until an acceptable level of error is attained and the flow in each airway is known. 2) Computes the pressure drop (head loss) in all airways and the overall network and mine head losses. 3) Determines the operating point for the mine and a given fan, adjusting all air flows and head losses to the actual values. The basis of analysis used is the Hardy Cross method of successive approximation for solution of any fluid network, explained in detail in the earlier article. In brief, corrections to assumed air flows are computed for all the loops of a network; following distribution to the airways of each loop, the process is repeated until acceptable accuracy is attained. When a hand or desk-calculator solution is used, there is advantage to close approximation of the actual flows in the original assumption of air-flow distribution in the network. Since publication of the previous paper, in which assumption of initial quantities was also necessary in the computer solution, the program has been modified to make the first distribution itself, assigning the total quantity to one specified path in a split and then proceeding to correct the flows. The number of iterations required for a high-accuracy solution is little affected, at least in networks of few loops, 10 to 15 being sufficient to reduce the error to 1 pct or less. The principal limitation of the original program was that it contained no provision for handling a fan. Further, it did not determine the overall head loss in the network and was not designed for multiple parallel openings, such as occur in coal mines. The revised program was written to overcome these disadvantages and also to print out the adjusted flows for all airways, corrected for the system operating point. At the present time, the program can handle only an external fan. The procedure employed in solving a ventilation network is now as follows: 1) Input data consisting of airway characteristics and loop designation along with assumed direction of flow, for each branch of the network and any trunk airways, are fed into the computer. 2) Employing the resistance relationship, where R is resistance, K is friction factor, P is perimeter, L is length, A is area, and n is number of airways (where multiple intakes and returns are employed); values are computed for all airways. 3) An arbitrary, continuous route of flow through the circuit—from beginning to end—is selected and the path identified by airway numbers and read in. 4) Based on a total assumed quantity of flow Q (usually 100,000 cfm.) the computer employs the procedure described in the previous article in distributing quantities of flow among all the airways until a satisfactory level of accuracy has been attained (in the current program, the machine is turned off when the cumulative corrective error in cfm, is reduced to five times the number of airways). 5) Airway head losses H are calculated by the relation H = R Q2 and cumulated along the original designated flow route to find the overall head loss for the network. Trunk airways are then included to obtain the mine head.

Jan 1, 1964

By A. S. Yue

Small-angle tilt boundaries have been observed in a lamellar eutectic structure. These boundaries are morphologically similar to a theoretical till boundary originally proposed by Burgers. The production of a small-angle tilt boundary in the lamellar eutectic structure is attributed to the formation of an indentation on a planar solid-liquid interface during growth. Constitutional supercooling due to the accumulation of impurity atoms ahead of the interface is primarily responsible for the generation of an indentation. The distance between two faults of like sign, designated as lamellar fault distance, is calculated by the Burgers equation for small-angle boundaries. It was found that the lamellar fault distance decreases exponentially and the tilt angle increases exponentially with freezing rate. SMALL-ANGLE tilt boundaries have been observed in many crystals in the as-grown states by etching techniques. However, a detailed substantial proof for revealing the presence of dislocations in the form of etch-pits was recently accomplished by Vogel. These boundaries are described as an array of dislocations of like sign. The angle of tilt associated with a boundary is treated theoretically by Burgers3 and is given by 0 = b/D for small tilt-angle where b is the Burgers vector of the dislocations and D is their average spacing. During the course of experimentation on the solidi fication of the lamellar Mg-32 wt pct Al eutectic a microstructure which was strikingly similar to a dislocation tilt boundary was observed. An example of the type of microstructure is given in Fig. 1, which shows a combination of transverse and top sections of the lamellar eutectic structure. The lamellae of the top section on one side of the apparent tilt angle are not parallel with those on the other side. Consequently, some lamellae must be terminated somewhere. Fig. 1 shows that the termination of lamellae is analogous to a dislocation tilt boundary in a crystal lattice originally proposed by Burgers. This work represents a closer look at such structures. The factors governing the production of the small-angle tilt boundary are discussed in terms of the freezing condition and impurity distribution in the melt. The effect of freezing rate on the distance between two faults, designated as lamellar fault distance, and the tilt angle is also described. In addition, the lamellar fault distance is determined by means of the Burgers equation for small-angle tilt boundary. EXPERIMENTAL PROCEDURE The slightly impure Mg-32 wt pct Al eutectic rods, 5/8 in. in diam and 12 in. long, were made from 99.99 wt pct Al and 99.98 wt pct Mg and were cast in a high-purity graphite mold. Some of these rods were purified by zone-refining the eutectic using direct-radiation helical-coil furnaces5 A molten zone about 2 in. long was passed through the cast rods at a speed of 5.3 x 104 cm per sec for ten passes of the molten zone. Before the freezing experiment was begun, 2-1/2 in. of the material from the head and tail ends of the zone-refined eutectic were removed from the rod to eliminate any impurity. The freezing experiments were carried out under controlled conditions,8 using a protective atmosphere of purified argon. At the end of each run, observation of the tilt boundary was made throughout the entire length of the rod with a light microscope on the top, longitudinal, and transverse sections (three mutually perpendicular sections) of the lamellar eutectic solid. The experimental data used for the determination of interlamellar spacing and tilt angle were taken from the middle portion of

Jan 1, 1964

By Y. H. Liu, G. W. Powell

Thin foils made from the gold-rich region of- short (0.15 mm)Au-Ag diffusion couples annealed at 880°C for 1 hr were examined in ax electron microscope and the structures were compared with that of pure gold subjected to identical heat treatments. The couples underwent Plastic deformation (bending) during the diffusional anneal. The foils prepared from the diffusion zone have a higher dislocation density and incidence of subboundaries than pure gold. The plastic bending of the couple can be interpreted in terms of the climb of edge dislocations, which are pavallel to the bend axis, along the direction of diffusion. THE interaction of lattice defects during interdiffu-sion between two metals has been considered by several investigators. In particular, da Silva and Mehl,1 Barnes,' and Bardeen and Herring3 have proposed models based on interactions between point and line defects to account for the Kirkendall effect (marker movement). Each of these models is based on dislocation climb, which can either add or remove lattice planes, induced by an under- or supersaturation of vacancies. A theoretical treatment of the nucleation and growth of voids by the condensation of vacancies has been presented by Seitz.4 Stresses can arise in a binary diffusion couple as a result of the vacancy flux and the gradient in the lattice parameter. rinkmman' has investigated the effects of diffusion-induced stresses on void formation; the results of his analysis suggest that voids can form and grow under the influence of these stresses in the absence of an excess vacancy concentration. The relationships between vacancy supersaturation and diffusion-induced stresses and void formation have been elucidated also by Balluffi and seigle6 who concluded that the stress generated by a vacancy supersaturation aids the growth of voids. Experimental evidence which indirectly shows that stresses are produced by chemical diffusion has been obtained by several researchers.2,7-9 This evidence consists in part of the observation by light microscopy of the formation of subgrains, a process characteristic of the high-temperature deformation of some metals during creep.'0 More striking, however, is the plastic bending of short (0.2 mm) Cu-Ni couples observed by Barnes2 and of boron-doped silicon crystals reported by Queisser.9 Diffusion-induced dislocations in iron and silicon have been investigated by Ainslie et a1.11and Schwutte and Queisser,12 respectively, the former using transmission electron microscopy and the latter an X-ray transmission technique. Matthews and Crawford13 observed the formation of boundaries by the climb of interfacial misfit dislocations in thin, vapor-deposited Au-Pd couples diffused in situ in an electron microscope. In order to obtain additional information on the phenomena associated with defect interaction and diffusion-induced stresses in binary couples, an electron transmission study of foils prepared from the diffusion zone of short Au-Ag couples was undertaken. The initial objective was to investigate the processes occurring at various points of the diffusion zone. Subsequent research demonstrated that foils could be prepared only from the gold-rich region of the diffusion zone using the technique described herein. EXPERIMENTAL PROCEDURE Preparation and Annealing of the Diffusion Couples. The couples were made from annealed, high-purity (99.999 pct). polycrvstalline gold and silver sheets (0.075 mm thick). Shall (0.075 by 5 by 17 mm) elec-tropolished coupons of gold and silver were pressure-bonded by heating the two metals in a graphite device to 770°C and holding at this temperature for approximately 30 min. The quality of the bond was ascertained by metallographic examination at 150 times. Each couple subsequently was annealed at 880°C for 1hr; under these conditions, the extreme portions of the

Jan 1, 1968

By M. E. Wadsworth, C. H. Pitt

Nickel corrosion in sulfuric acid solutions at elevated temperatures and oxygen wer-pressures was imestigated. Weight loss was 1inear with time and varied directly with oxygen concentration. Independent variation of the possible active species in solution indicated that the undissociated sulfuric acid molecule plays an important part in the corrosion mechanism. The rate determining step was considered to be the adsorption of oxygen on a site containing adsorbed H2SO4. THE tendency for a metal to go into solution depends on a number of factors some of which are: 1) the position of the metal in the electrochemical series; 2) the type and concentration of corroding species in solution; 3) the concentration of oxygen in solution; 4) the type of film formed at the liquid-metal interface; and 5) the physical state of the metal. The corrosion of a metal in an aqueous solution may be considered to be a heterogeneous reaction that occurs at the liquid-metal interface. The reaction is generally considered to occur as a result of electrochemical potential differences within the system. Little work has been published in the literature concerning the corrosion of pure nickel. A great deal has been publistled dealing with the experimental and theoretical aspects of the over-all corrosion process. Evans,1 Uhlig,' and recently Todt 3 and Tomashoff,4 have published books which review the results of experiments which have been conducted on corrosion, and have developed some of the basic ideas of the more widely accepted theories of corrosion. Friend5 states that. nickel does not readily discharge hydrogen from any of the common nonoxidiz-ing acids and that a supply of some oxidizing agent is necessary for appreciable rates of corrosion. He further notes that at temperatures between atmospheric and close to boiling, corrosion rates tend to vary with concentration, temperature, and aeration in acid concentrations up to 25 pct. Whitman and Russel6 studied the effect of dissolved oxygen on a number of metals in acid solutions and found that oxygen accelerates the rate of corrosion of nickel in sulfuric acid. Flowers and Kelley7 in studying the action of dissolved CO, on electroplated nickel postulated that nickel dissolves as NiHCO3 in the presence of CO, and air. It was found that oxygen was necessary for corrosion to take place and that the rate increased with increasing amounts of oxygen in solution. Flowers and Kelley attributed the oxygen effect to cathodic depolarization. MacGillavry, Rosenbaum, and wensson' determined metal solution potentials of nickel in various electrolytes. They found that nickel has no tendency to corrode in alkaline solutions in the presence or absence of air. In acid solutions in the presence of air, metal solution potential increases in the order of 1/2 v were observed. This increase was attributed to probable chemisorption of oxygen on the nickel surface and the corrosion mechanism suggested was the interaction of H+ ions with the nickel oxide formed at the surface. Delahay9 observed the reduction of oxygen on various metals. A polarograph was used to measure the course of the reaction. With some metals hydrogen peroxide was produced as an intermediate compound. However, nickel was not considered to be a hdrogen peroxide producer. Van Rysselberghe 10,11 and associates attributed the accelerating effect of carbon dioxide on the corrosion process to the presence of percarbonic acid produced in solution reacting with the hydrogen peroxide. The percarbonic acid can be reduced at less negative potentials than hydrogen peroxide and as a result accelerates the reaction. EXPERIMENTAL Since a complete description of the apparatus used can be obtained elsewhere,12'13 no attempt will be made here to describe the equipment other than to note that it consisted essentially of an autoclave unit from which samples could be withdrawn periodically without disturbing the interior state of the autoclave. A titanium autoclave liner was used to protect the autoclave from the solution.

Jan 1, 1961

By Paul B. Crawford, Melvin B. Burton

A procedure is given which describes how the gelatin model may be-used to study the effect of mobility ratio on areas swept before and after breakthrough in fluid injection programs. Differences in mobility ratio are achieved through the proper selection of electrolytes for the invading and displaced fluids. The progress of the front is followed visually and photographed or traced onto plates (it desirrd positions before and after breakthrough. Example applications to standard flooding pattrrrls, vertically and horizontally fractured reservoirs and the in-situ combustion recovery process ore presentrtl. INTRODUCTION In water flooding, gas cycling and other fluid injection programs, the sweep efficiencies and swept areas after breakthrough are known to depend on the mobility ratio. Previous investigators have described the application of the X-ray technique, stepwise use of the poten-tiometric model and numerical procedures for studying these problems.1,3,4,6 Others have described the application of the electrolytic and gelatin models for studying sweep efficiencies at a mobility ratio near one.2,5,7,8,9 It is a purpose of this paper to describe a procedure for using the gelatin model to study swept areas at mobility ratios other than one. PROCEDURE The technique developed for gelatin model studies is composed of four phases: the construction of the model, the selection and preparation of electrolytes, the running procedure, and the recording of the swept areas. Each phase is discussed independently in order that the over-all procedure will not be confused. Construction of the Model Since some of the advantages of the gelatin model are its economy, and speed, it is desirable to have as small a model as is possihle without creating error through the inability to follow visually the progress of the colored front. For this work, it was found that a distance of 4 to 6 cm between salt bridges is convenient to use. Since the size is relatively small, it follows that the gelatin layer must be thin. For this study, 1/8 in. was used as the gelatin thickness in each case. The model was constructed on a flat 1/2-in. plastic (Lucite) base with boundaries constructed from 1/8-in. plastic strips. The well pattern to be studied was outlined with the 1/8-in. plastic strips on the base with narrow apertures (- - 2 mm) in the boundary to represent the wells. Behind each well opening, a rectangular area of some 8 to 12 sq cm was outlined by additional plastic strips. The plastic strips were fastened to the base by any suitable non-conducting cement so that the boundary system was liquid tight. The electrodes were constructed from sheet copper in a fashion that would allow easy removal. In this work, it was found that by simply forming a large U-shaped clip from a strip of the copper that would slip over the boundary was satisfactory and easily removed. The electrodes should be as wide as one of the sides of the rectangular area behind the well opening. If a fractured pattern is studied, fractures inside the well pattern may be represented by strips of copper of the desired shape mounted in some easily removable fashion. The reason that apertures (or salt bridges) were used for wells instead of metallic electrodes is that the plating of the electroactive ion at the cathode and the introduction of an electroactive ion at the anode has a detrimental effect on the swept area outline if these operations \rere permitted to be carried on within the well pattern itself. The anode may be used within the well pattern

Jan 1, 1957

By Y. Kondo, A. M. Gaudin

A combined method for sizing in the sub-micron range was developed using centrifugation and gamma emission. Its essential features are: 1) Centrifugation of suspended mineral in water carried out at 4°C where water has its maximal density and zero temperature gradient of density 2) Gamma activition of the sample suspension induced by irradiating it in a nuclear reactor prior to centrifugation 3) Measurement of the activity of a beam of suspension as it is sedimenting in the centrifuge, using a scintillation counter provided with a discriminator. Albite, the sodium feldspar, was used as feed to the comminution step, and the induced gamma activity of Na24 was measured. Dissolved sodium ion caused a a non-settling activity in the measurement. Other "noise" is due to scattered gamma radiation from albite particles in the centrifuge outside of the line of sight. These two activities were obtained sep arately and subtracted from the total activity counted. The relative solid concentration in the suspension under the collimating slit is proportional to this net activity. The cumulative weight per cent of undersize particles was calculated and plotted against size on logarithmic scales. The relationship was found to be linear down to about 30 mµ. The specific surfaces calculated from particle size distribution agree with those measured by krypton gas adsorption. The measurement of particle size distribution of comminution products has been of interest in clarifying the mechanism of the process. The Gaudin-Schuhmann equation1,2 relates the cumulative percentage undersize, y, to relative particle size x/k. The particle size, x, and the size modulus, k, together with the distribution modulus, m, characterize the particle size distribution, m is the slope of a plot of log (y) against log (x/k). This equation has been verified for ground particles down to 1 micron in size.3'4 For particle sizes smaller than a few microns, the precision and facility of sedimentation and microscope techniques become questionable. Use of a centrifuge which shortens settling time from weeks to hours is attractive. Many centrifugal methods have been proposed; however, determination of the solid concentration of the suspension at a reference level as a function of time is difficult, and sampling the suspension is likely to introduce errors. Measurement of particle size distribution by X-ray absorption in gravitational5 and in centrifugal6 sedimentation have also been tried. In these methods, the absorption coefficient for the solid is independent of size, since the wave lengths of the X-rays and gamma-rays are orders of magnitude smaller than the particles. Radiometric techniques7 in size distribution measurement have been developed within the last several years. They offer advantages over usual gravimetric and optical methods. PRINCIPLE OF METHOD The combined method of centrifugation and gamma emission for sizing in submicron range described in the present paper has the following essential features. 1) Centrifugation is carried out at 4 ± 0.1°C. Convection currents in the water phase due to density changes caused by temperature variation may disturb completely the settling system. To avoid this difficulty, the unique property of water having zero temperature gradient of density at 4°C is used. 2) Prior to size measurement the sample suspension is irradiated in a nuclear reactor and gamma activity induced. The gamma activity of the particles sedi-mented radially is measured by a scintillation counter provided with a discriminator. By this method, interruption of centrifugation by withdrawing samples is avoided and the counting rate at any time represents the solid concentration in the suspension at the level of the collimating slit of the centrifuge at the particular time at which the measurement is made. 3) By proper selection of the feed material and of the comminution device, gamma activities of foreign

Jan 1, 1964

By J. Lohrenz, R. J. Francis, G. C. Clark

Compositional material balance methods have been extended to solution gas-gas cap-water drive reservoirs. Partial or complete pressure maintenance as well as reservoir repressuring by gas and/or water injection may be considered. Increrments of reservoir withdrawals and injections equivalent to a given time increment are used as the basic iterative step to allow the extension. An example field study shows how the reservoir volumes. extent of the water drive, and relative permeability ratio function may be estimated. Knowing these, the methods can be used to predict reservoir behavior. INTRODUCTION The compositional material balance is similar to conventional material balances in many aspects. The assumption of a pressure-homogenous reservoir with no transients is common to all conventional material balances1-1 although numerical methods are different. The compositional material balance also contains this assumption. Inclusion of phase equilibria directly in reservoir calculations is the key difference and improvement in the compositional material balance. A continuous account of the compositions of the reservoir fluids and effluents is determined. Calculations of fluid phases and properties can be based on these compositions. Excepting fluid phases and properties determined in this manner, the compositional material balance uses procedures exactly analogous to conventional methods. Brinkman and Weinaug5 and Jacoby and Berry6 have given compositional material balances for simple solution gas drive reservoirs; subsequently Jacoby and Berry' extended their procedure to include pressure maintenance by gas injection at constant reservoir pressure. This paper extends compositional material balance techniques to a more general combination drive reservoir with various types of pressure maintenance operations. The procedure for this extended compositional material balance and an example of its field application are presented. PROCEDURE The numerical calculations of this compositional mate- rial balance are performed on a high-speed digital computer. The most comprehensive case the present compositional material balance program can handle is a solution gas-gas cap-water drive reservoir with (1) water production, (2) water influx, (3) separate or commingled production from the oil and gas zones, (4) any surface processing arrangement including LPG and gasoline production, and (5) gas injection into the gas and/or oil zones." This description of the procedure applies to this most comprehensive case. Appropriate deletions in the procedure can be made for more simple mechanisms. ASSUMPTIONS In common with conventional material balances, the compositional material balance makes the following assumptions: (1) the total reservoir is at a uniform pressure at a given time; (2) the volumetric ratio of gas to liquid withdrawal at reservoir conditions is directly proportional to their mobilities; and (3) instead of assuming a pressure-dependent total gas-oil equilibrium, the cornpbsitional material balance makes a sianificantly different and more rigorous assumption as follows. Thermodynamic equilibrium exists for each component in the liquid-gas mixtures in the reservoir and in surface processing. These three assumptions are fundamental in the compositional material balance procedure. Other assumptions involve arbitrary simplifications of the reservoir mechanisms. This compositional material balance uses the following simplifying assumptions: (1) no bulk mixing occurs between the gas and oil zone; (2) total mixing occurs between any injected gas and the contents of the injection zone; (3) no in situ hydrocarbon liquid is produced from the gas zone; (4) no reservoir gas is retained behind any net bottom-water influx; (5) the residual saturation of hydrocarbon liquid at reservoir conditions behind any net bottom-water influx is constant; and (6) fluids retained behind any net bottom-water influx exert no further influence on the reservoir. These six assumptions are not fundamental. Modifications of the present procedure can be made for different assumptions of the reservoir mechanisms. THE TIME INCREMENT Previous compositional material balances have used a pressure increment for simple solution gas drive reservoirs"" and an increment of injected gas for pressure maintenance.' These increments are suitable for these simple cases. They are, however, inadequate for reservoirs with time dependent drive mechanisms including water drive. A time increment removes this inadequacy and is

By M. R. Tek, R. L. Nielsen

The scaling laws as formulated by Rapoport relate dynamically similar flow systems in porous media each involving two immiscible, incompressible fluids. A two-dimensional numerical technique for solving the differential equations describing systems of this type has been employed to assess the practical value of the scaling laws in light of the virtually unscalable nature of relative permeability and capillary pressure curves and boundary conditions. Two hypothetical systems — a gas reservoir subject to water drive and the laboratory scaled model of that reservoir — were investigated with emphasis placed on water coning near a production well. Comparison of the computed behavior of these par ticular systems shows that water coning in the reservoir would be more severe than one would expect from an experimental study of a laboratory model scaled within practical limits to the reservoir system. This paper also presents modifications of the scaling laws which are available for systems that can be described adequately in two-dimensional Cartesian coordinates. INTRODUCTION Present day digital computing equipment and methods of numerical analysis allow realistic and quantitative studies to be carried out for many two-phase flow systems in porous media. Before these tools became available the anticipated behavior of systems of this type could be inferred only from analytical solutions of simplified mathematical models or from experimental studies performed on laboratory models. To reproduce the behavior of a reservoir system on the laboratory scale, certain relationships must be satisfied between physical and geometric properties of the reservoir and laboratory systems. Where the reservoir fluids may be considered as two immiscible and incompressible phases, the necessary relationships have been formulated by Rapoportl and others.2-5 Rapoport's scaling laws follow from inspectional analysis of the differential equation describing phase saturation distribution in such systems. It will be recalled that these scaling laws presuppose three conditions: (1) the relative permeability curves must be identical for the model and prototype; (2) the capillary pressure curve (function of phase saturation) for the model must be linearly related to that of the prototype; and (3) boundary conditions imposed on the model must duplicate those existing at the boundaries of the prototype. These three requirements seldom if ever can be satisfied in scaling an actual reservoir to the laboratory system because: (1) The laboratory medium normally will be unconsolidated (glass beads or sand) while the reservoir usually is consolidated. Relative permeability and capillary pressure curves are usually quite different for consolidated and unconsolidated porous media. (2) The reservoir usually will be surrounded by a large aquifer which could be simulated in the laboratory only to a limited extent. (3) Wells present in the reservoir would scale to microscopic dimensions in the laboratory if geometric similarity is to be maintained. In view of these considerations, rigorous scaling of even a totally defined reservoir probably would never be possible. The purpose of this paper is to assess the practical value of the scaling laws in the light of the unscalable variables. This has been done by carrying out numerical solutions in two dimensions to the differential equations describing the flow of two immiscible, incompressible fluids in porous media for a field scale reservoir and a laboratory model of that reservoir. While both the reservoir and the laboratory model were purely fictional, each has been made as realistic and representative as possible. The field problem selected as the basis for the investigation was an inhomogeneous, layered gas reservoir initially at capillary gravitational equilibrium and subsequently produced in the presence of water drive. The laboratory model of this reservoir was designed to utilize oil and water in a glass bead pack. The numerical treatment employed was similar to that of Douglas, Peaceman and Rachford6 and it included both capillary and gravitational forces as

By E. D. McMurry

INTRODUCTION Before discussing the use of the free piston; it may be well to explain briefly the type of gas lift valves with which the loo] is designed to function. The gas lift valve ii nothing more or less than a small back pressure regulator. It; construction consists of an outer hous ing covering a seamless metallic bellows in such a way as in form a sealed pressure vessel or pressure dome. To the bellow* in attached a stem and to this stem is attached a .rating member. When the valve is run into the hole, a pre-delermined pressure charge is placed in the pressure dome. As shown in Fig. 1. the valve is furnished in four different major assem-blies. to accomplish good results in different types of wells. The valve is mounted on a pup joint (Fig. 2a) : the eoin-plele assembly is entered into the tubing siring at the desired location. Each valve is furnished with a revere flow check valve (Fig. 2-b) to eliminate any back wash through the mechanism. Valve action is as follows: At any time the casing pressure is allowed to exceed the pressure forces within the pressure dome of the valve, the valve simply opens and passes gas from the rasing annulus into the tubing. When the casing pressure forces arc allowed to bleed to a force less than the pressure force within the dome, of the valve. the valve closes. With this type of construction. it is possible to open and close the valve or valves at will from the surface, the medium of the gas pressure in the rasing annnlus. The cycle of operation is governed by a tine cycle motor valve, commonly referred to as a gas lift interutitter or gas life surfacr controller (Fig. 3). In actual practice. after the well is unloaded only one valve is in operation . . . all other valve.: remain in the closed position. On pre-determined time cycles, it is possible to open and close this operating valve at will, and to pass the desired amount of gas through this valve wheneler needed. This method of gas lift is widelv used and has been found In be very efficient. However, it has long been known that even greater efficiency would be obtained if it were possible to place Piston between the injected gas and the fluid head. The free piston which will be discussed here was designed to accomplish this end and has been found to do so with very satisfactory results. DESIGN AND OPERATION OF THE FREE PISTON Basically. the free piston consist- of a telescoping assembly in which are incorporated the various control and operating parts (Fig. A set of sleeve valves which are opened and closed by movement of the telescoping section is contained in the. tool. A cross section of the free piston is shown in Fig. 4-11, 'rile main body cotain- a seamless monel bellows of the same type as that generally used in pressure charged gas lift valves (Fig. 4-c), and a valve core to permit placing a pressure charge in the bellows before the free piston is placed in service. The bellows is the main control unit of the free piston and performs two function: 1. It governs the opening and closing of the sleeve valves. a.: dictated by the producing activities of the well in which the piston is operating. 2. In its normal position, it holds a rubber packer in an elongated position to permit passage of fluid between the packer and the tubing walls. The lower section of the free piston contains a hydraulic shock absorber which cushions the return of the tool to its

Jan 1, 1953

By Stanley J. LeFond, Langtry E. Lynd

Elemental titanium has become famous as a space age metal, because of its high strength/ weight ratio and resistance to corrosion. However, the major use is in the form of titanium dioxide pigment, which because of its whiteness, high refractive index, and resulting light-scattering ability, is unequaled for whitening paints, paper, rubber, plastics, and other materials. A relatively minor use is in welding rod coatings, in the form of the mineral rutile. United States consumption of titanium dioxide in concentrates for these applications in 1982 was as follows: pigments, 630,000 st; titanium metal, 31,000 st; and welding rod coatings, 5,000 st. The only commercially important titanium ore minerals at the present time are ilmenite and its alteration products, and rutile. Titanium was discovered by Gregor in 1790, as a white oxide which he recovered from menaccanite, a variety of ilmenite occurring as a black sand near Falmouth, Cornwall. Barksdale (1966) stated that the fundamental chemical reactions on which the present day titanium industry is based were known before 1800, although it was not until 1918 that these pigments were available commercially on the American market. Pings (1972) outlined the early history of the titanium industry in the United States, referring to the work of Guise (1964) and others: Mining of titanium minerals in the United States began sometime between 1880 and 1900, in Chester County, Pennsylvania. Small quantities of rutile were also produced during that time in North Carolina and Georgia. In 1901 rutile was mined from a deposit near Roseland, VA, and was used in making titanium chemicals and for coloring ceramics. Ilmenite in the deposit was produced as a separate item in 1913, and rutile and ilmenite were obtained from this deposit through 1921. The mining of titanium bearing beach sands began in 1916 near Mineral City and Pablo Beach, FL, for the purpose of making titanium tetrachloride to be used in tracer bullets, flares, and smokescreens. The titanium pigment industry was founded in 1918. By 1928 a large part of the domestic production of rutile, ilmenite, and zircon came from Florida. However production ceased in 1929 in Florida with the mining of newly discovered deposits in Virginia which were mined until 1968. The discovery of new deposits in Florida, and the development of better mining and concentrating methods for low-grade sands, led to a return of activities in this area. The large deposits near Tahawus, NY. were brought into production in 1942 by National Lead Co. (now NL Industries, Inc.), and by 1949 this company was the leading producer of ilmenite in the world. In 1948 the extensive deposits of ilmenite and ilmenite-bearing iron ores were discovered in eastern Quebec. The large Trail Ridge sand deposit in Florida was discovered by E. I. du Pont de Nemours & Co., Inc. geologists in cooperation with the US

Jan 1, 1983

By Frank L. Gaddy

Falls of roof account for over 50% of the fatalities that occur in coal mines in the U.S. Thus, roof control is one of the more important phases of underground mining. In reality, the control of roof influences the system of mining and is a major determinant of the width and spacing of working places, operations at the mining face, ventilation control, and surface subsidence. Frequently, the control of roof is the largest single cost item. Roof control is a never-ending task, not only at the working face and throughout sections where men are working, but along haulways and airways that must be maintained for the life of the mine. This is because roof, even the best of roof, slowly deteriorates, so it must be frequently examined, with corrective support applied when needed. As coal is excavated, stresses are set up in the roof became the previous equilibrium is upset, resulting in pressures that cause fractures and slight, movements that are, frequently, hard to detect. Unless the immediate roof, in the excavated area, is given support by artificial means, it might fall or there might be a succession of falls, varying from a thin scale to several feet, depending on the nature of the top. Exceptions to this are those rare mines, or sections of mines, with hard, strong, nonweathering roof that requires no support. GENERAL CLASSIFICATION OF ROOF There are two broad types of roof as far as support is concerned, the immediate roof above the coal and the main roof. Artificial support for mining purposes is only concerned with the immediate roof as nothing except large blocks of solid coal, or massive concrete, will support the main roof. The immediate roof is generally a few feet thick but can vary from inches to 20 ft or more. There are examples of where there is no immediate roof as the massive main roof lies directly on the coal bed, but this is rare. Impossibility of supporting the main roof and all of the cover

Jan 1, 1973

By Gerald Sherman

MINING asbestos ore by caving at Thetford Mines, Quebec, has been described in the Transactions of the Canadian Institute of Mining and Metallurgy in papers presented by the staffs of The Asbestos Corp. Ltd., in 1934, and Johnson's Co. Ltd., in 1941. The production by caving prior to 1949, in round figures, has been 14,500,000 tons of ore by The Asbestos Corp. and 5,500,000 by Johnson's Co. They are the only caving operations in asbestos at the present time. The Jeffrey mine of the Canadian Johns-Manville Co. at Asbestos, 40 miles southwest of Thetford, is well along in preparations for caving, and The Bell Asbestos Mines Ltd. at Thetford Mines, will follow their examples later. GERALD SHERMAN, Member AIME, is Consulting Engineer, Johnson's Co. and Canadian Johns-Manville Co., New York, N. Y. AIME San Francisco Meeting, Feb. 1949, and New York Meeting, Feb. 1950. TP 2828 A. Discussion (2 copies) may be sent to Transactions AIME before May 31, 1950. Manuscript received Nov. 2, 1949. Nearly all of the asbestos produced by the western hemisphere is mined in those districts of Canada. Since it has been the subject of a number of technical papers, this contribution is offered as a review of progress. Some introductory information seems appropriate, for few members of the American Institute of Mining and Metallurgical Engineers are familiar with the mineral, or how asbestos deposits differ from others in caving characteristics. Asbestos is found in peridotite, a basic intrusive rock irregularly altered to serpentine and erratically fractured, well enough to cave in bulk, but leaving many blocks of solid ore of all sizes scattered through the deposit. Chrysotile is the variety of asbestos produced in Canada. Its fibers are the most valuable commercially, having a good proportion of long fibers, strong, silky, flexible, and suitable for spinning. The fibers occur in all lengths, from the microscopic to 3 in., or more, the shorter are more abundant. Their values depend on length combined with other qualities noted above. They occur as veinlets filling fissures in the peridotite, with crystals compactly arranged in parallel across the cracks. They are sometimes lightly frozen to the walls but the partings are distinct, often marked by particles of magnetite adhering to the fiber ends. Sometimes they break away from the rock in slabs an inch or more in thickness. In wide fissures the fibers are not always continuous, the "broken fibers" being separated by lines of parting with occasional scattered grains of magnetite, which is an injurious mineral to be elirninated when possible. In addition to the chrysotile ''Cross fiber," some "slip fiber" is found in fault fractures in limited quantities of less valuable character. With the exception of the Asbestos Corp, ore, in which the content is said to be more uniform, asbestos is eccentrically distributed through the peridotite. It is the general impression that the higher grades are more often found in hard rock. There are rich pockets and leaner areas which range down through various grades to minor barren sections. Wholesale mining is necessary to make sure of recovering all the high grade for there are no reliable indications of its existence. Ore selection only appears in constructing assay-grade limits for caving. The mill yield of the Thetford ore body in 1934 was 6 to 10 pct of the ore milled, made up of all classes, with occasional runs up to 15 pct or so. Quotations for Canadian asbestos* are as follows: Per Ton No. 1 Crude (% in. and longer) $896 @ $960 No. 2 Crude (Y8 to % in.) 350 @ 545 Spinning (less than % in.) 204 @ 378 Other classes, price diminishing with length. Shorts or "refuse" 24 @ 46 • From Eng. and Min. Jnl., Nov. 1948. Asbestos of even superior quality for electric insulation is found in Arizona limestone. It is white, containing very little combined iron, but occurs only in small quantities.

Jan 1, 1951

By Gerald Sherman

MINING asbestos ore by caving at Thetford Mines, Quebec, has been described in the Transactions of the Canadian Institute of Mining and Metallurgy in papers presented by the staffs of The Asbestos Corp. Ltd., in 1934, and Johnson's Co. Ltd., in 1941. The production by caving prior to 1949, in round figures, has been 14,500,000 tons of ore by The Asbestos Corp. and 5,500,000 by Johnson's Co. They are the only caving operations in asbestos at the present time. The Jeffrey mine of the Canadian Johns-Manville Co. at Asbestos, 40 miles southwest of Thetford, is well along in preparations for caving, and The Bell Asbestos Mines Ltd. at Thetford Mines, will follow their examples later. GERALD SHERMAN, Member AIME, is Consulting Engineer, Johnson's Co. and Canadian Johns-Manville Co., New York, N. Y. AIME San Francisco Meeting, Feb. 1949, and New York Meeting, Feb. 1950. TP 2828 A. Discussion (2 copies) may be sent to Transactions AIME before May 31, 1950. Manuscript received Nov. 2, 1949. Nearly all of the asbestos produced by the western hemisphere is mined in those districts of Canada. Since it has been the subject of a number of technical papers, this contribution is offered as a review of progress. Some introductory information seems appropriate, for few members of the American Institute of Mining and Metallurgical Engineers are familiar with the mineral, or how asbestos deposits differ from others in caving characteristics. Asbestos is found in peridotite, a basic intrusive rock irregularly altered to serpentine and erratically fractured, well enough to cave in bulk, but leaving many blocks of solid ore of all sizes scattered through the deposit. Chrysotile is the variety of asbestos produced in Canada. Its fibers are the most valuable commercially, having a good proportion of long fibers, strong, silky, flexible, and suitable for spinning. The fibers occur in all lengths, from the microscopic to 3 in., or more, the shorter are more abundant. Their values depend on length combined with other qualities noted above. They occur as veinlets filling fissures in the peridotite, with crystals compactly arranged in parallel across the cracks. They are sometimes lightly frozen to the walls but the partings are distinct, often marked by particles of magnetite adhering to the fiber ends. Sometimes they break away from the rock in slabs an inch or more in thickness. In wide fissures the fibers are not always continuous, the "broken fibers" being separated by lines of parting with occasional scattered grains of magnetite, which is an injurious mineral to be elirninated when possible. In addition to the chrysotile ''Cross fiber," some "slip fiber" is found in fault fractures in limited quantities of less valuable character. With the exception of the Asbestos Corp, ore, in which the content is said to be more uniform, asbestos is eccentrically distributed through the peridotite. It is the general impression that the higher grades are more often found in hard rock. There are rich pockets and leaner areas which range down through various grades to minor barren sections. Wholesale mining is necessary to make sure of recovering all the high grade for there are no reliable indications of its existence. Ore selection only appears in constructing assay-grade limits for caving. The mill yield of the Thetford ore body in 1934 was 6 to 10 pct of the ore milled, made up of all classes, with occasional runs up to 15 pct or so. Quotations for Canadian asbestos* are as follows: Per Ton No. 1 Crude (% in. and longer) $896 @ $960 No. 2 Crude (Y8 to % in.) 350 @ 545 Spinning (less than % in.) 204 @ 378 Other classes, price diminishing with length. Shorts or "refuse" 24 @ 46 • From Eng. and Min. Jnl., Nov. 1948. Asbestos of even superior quality for electric insulation is found in Arizona limestone. It is white, containing very little combined iron, but occurs only in small quantities.

Jan 1, 1951

By H. W. Mead, C. E. Birchenall

SELF-DIFFUSION in pure gold has been studied by Sagrubsky,1 McKay,2 Gatos and coworkers,3,4 and Okkerse.5 Sagrubsky used a sectioning technique but reported results far lower than any of the others. The work seems to be carefully done, and it appears that a factor of 60 must have been omitted in his calculations or the units misstated by an equivalent amount. McKay used a surface activity method, and Gatos and Kurtz, following Gatos and Azzam, applied autoradiography to the measurement of activity on an oblique section through the diffusion zone. These two investigations agree well with each other but give results below those reported here. Okkerse employed a very careful sectioning technique. Chemical interdiffusion in Au-Ag alloys has been measured under several concentration conditions. Average values of the diffusion coefficient for the range of 0 to 18 atomic pct Au were obtained by Braune.6 since the diffusion coefficient probably varies markedly over this range, the values may not be too useful. Also, a plot of these data suggests that grain boundary contributions may have been appreciable at the lowest temperature employed. Ebert and Trommsdorf7 give average values for 0 to 9 atomic pct Ag, and Johnson8 reports average values for 45 to 55 atomic pct Au. Johnson also measured silver and gold self-diffusion coefficients in 49.2 atomic pct Au alloys, and self-diffusion in pure silver. The self-diffusion of gold in pure silver has been determined by Jaumot and Sawatsky.10 Other measurements which do not cover a range of temperature or are of low precision are not mentioned in this account since a review by Kubaschewski is available." This contribution reports diffusion coefficients for self-diffusion in gold, self-diffusion of gold in Au-Ag alloys, and chemical diffusion in Au-Ag alloys. Although it had been intended to obtain silver self-diffusion data in the alloys it now appears unlikely that these measurements can be done soon, so a table has been worked out for a self-consistent scheme of diffusion coefficients which should not be too far from the true values. Experimental Details Experiments on pure metals were conducted using fine gold of 99.99+ pct Au and silver of the same purity. Alloys were prepared by melting the two metals together in a graphite crucible in vacuo and allowing the slug to solidify in the crucible. After machining to a cylinder about 0.85 in. diam each slug formed one half of a diffusion couple. Compositions were determined later by analyzing chips taken from the extremities of the diffusion zones. The gold used in the preparation of alloys was sometimes pure gold but was normally either in the form of clean pickled alloy from previous experiments or gold recovered by oxalic acid precipitation from warm aqueous solution. Chips from the various stages of experimentation were dissolved in aqua regia and the solution taken to low bulk, diluted with water, and again taken to low bulk before diluting for the precipitation. This precipitated gold contained a small amount of silver as the only impurity detectable by spectrographic analysis. Disks for use in diffusion couples had opposite plane faces abraded on emery paper as far as zero grit, were annealed for 12 hr at 900°C in argon, and had one face of each disk polished to 000 emery paper. Radioactive gold was plated at about 2 ma per sq cm from a small volume of cyanide solution containing 1 mg of inactive gold as a carrier. Specimens , were freshly abraded on 4/0 emery paper immediately prior to plating. In the case of the active couples disks of a similarly prepared but unplated alloy or pure metal were placed on top and welded to the active disk. With the concentration couples disks of different compositions were welded. Welding was performed after a pressure of 5,000 lb had been applied at

Jan 1, 1958

By J. F. Abel

It can be concluded that rock mechanics instrumentation, geologic mapping, and operations research in combination will produce an accurate estimate of tunnel support requirements for establishing a steel-set load prediction model. This estimate can be used to increase safety and to reduce tunneling costs. Unnecessary supports can be eliminated; overloaded or insufficient supports can be replaced on a more scientific basis than previously possible. Only a limited number of measurements are required to develop a steel set load prediction model. The results of this study can be applied generally; but, every tunnel must be considered as a separate problem. The construction and geologic variables which were found to be significant in the prediction of rock loads in this tunnel were: steel section modulus, set blocking points, percentage of alteration, relative degree of faulting, relative water condition, distance to the nearest fault, and average joint spacing. INTRODUCTION A combined engineering geology and rock mechanics program was performed throughout the Straight Creek tunnel pilot bore in Colorado in order to develop a technique for estimation of steel tunnel support loads. The technique involves the construction of a statistical model from measured steel set loads, geologic and construction factors. The aim of the investigation was to produce a set of charts, based on these factors, which could be used by a trained man equipped with a hand lens and measuring tape, to specify a safe and economical steel size and set spacing. In order to be of practical value the measured parameters had to be rapidly measurable and capable of quantification. It was apparent from the beginning that quantitative values for many of the geologic factors would have to be arbitrary. It was impossible to ascertain in advance which of the readily determi-nable geologic and construction factors significantly affect steel set loading. As the investigation progressed the number of observations was systematically reduced in order to determine the minimum amount of data which would be necessary for a significant mathematical model. This paper presents an analysis of all the variables measured at each instrumented steel set. The statistical significance of the measured variables was determined. On the basis of this determination, a mathematical model was constructed and charts were prepared to select and space the steel tunnel supports at any given point in the tunnel. Once the effect of geologic and construction factors on set loading is sufficiently established, it is possible to exert close control over the selection of steel tunnel support size and spacing. The magnitude of this rock mechanics program can be appreciated when it is known that there were 44 instrumentation stations for which steel set loads were measured. For each of these sets there were three construction variables recorded, and nine geologic parameters determined (Table I). A rough approximation of the cost for both the collection and reduction of the data would be $50,000. This is in addition to the cost of geologic mapping.

Jan 1, 1967

By P. T. Bail, J. E. Bobek

In reservoirs containing soft shales interbedded with oil-producing sands, it is possible to have water movement from the shales into the sands during reservoir depletion. This paper presents the results of scaled flow-model tests performed to investigate the efficiency of oil displacement by this vertical influx of water into oil-producing sands. In these model tests, the effects of rate of influx and rock wettability were considered. Oil recoveries by shale-wafer influx and by water flooding were cornpared after 0.4 hydrocarbon volume of water entered the producing sand. This comparison was made for water-wet and oil-wet systems. In water-wet sand, shale-water influx produces slightly less oil than conventional water flooding; in oil-wet sand, shale-water influx produces about one-half as much oil as conventional water flooding. In another type of comparison, it was found that shale-water influx followed by water flooding recovers less oil than water flooding alone. This is true for both water-wet and oil-wet sands. A comparison of tests run at different rates of shale-water influx shows that oil recovery is insensitive to rate in water-wet sand; in oil-wet sand, oil recovery decreases with decreasing rate of influx. INTRODUCTION Localized subsidence of the earth's surface has accompanied the depletion of some oil reservoirs, particularly in California,1, 2 along the Texas Gulf Coast,3 4 and in the Bolivar Coastal fields of Venezuela. Subsidence occurs when fluid withdrawals reduce the pressure in the reservoir. The reasons for subsidence accompanying a reduction in reservoir pressure are not known with certainty.10° However, there is considerable evidence that compaction of soft shale beds adjacent to producing oil sands can be a major cause of this subsi-dence.3, 6, 10-12 When the shales compact, a limited amount of water moves from the shales into the producing sands. A schematic illustration of this oil-producing mechanism is shown in Fig. 1. As the pressure in the sand declines due to oil production, a predominantly vertical pressure gradient develops at the sand-shale boundary. This causes water to move from the shale into the sand as indicated by the solid arrows. This paper reports the results of scaled flow-model tests to investigate the efficiency of oil displacement by the vertical influx of shale water into producing sands. The model tests were not performed with the aim of duplicating the recovery behavior of any specific reservoir. The primary purpose of the tests was to obtain a better understanding of oil displacement by shale-water influx. A knowledge of the oil- and water-production characteristics of even a simplified reservoir system should be of value to engineers estimating the performance of reservoirs where shale-water influx is occurring in conjunction with other production mechanisms. The following limitations and assumptions apply to the model work. 1. The sand of the prototype reservoir is homogeneous and horizontal. 2. The shale-water influx is uniformly distributed along the top and bottom of the sand member being studied. 3. No effects of edge-water drive or dissolved-gas drive that might accompany shale-water influx are considered. PROTOTYPE RESERVOIRS Scaled-model tests were performed to study oil displacement from three prototype reservoir sections. The reservoir sections modeled consist of single sand layers out of a sequence of -sand and -shale layers, as illustrated in Fig. 1. Each prototype represents a homogeneous, high-permeability, unconsolidated sand layer. In the field, the prototype is oriented in such a manner that one end is at the radius of drainage of a producing well and the other end is pointed toward the producing well (Fig. 2). The sand and fluid properties and the dimensions of the three prototypes are given in Table 1. The capillary pressure relationships for the three prototypes are presented in Figs. 3, 4 and 5. The oil-water per-