Search Documents

By J. S. Levine, M. Prats

A homogeneous and uniform cylindrical reservoir containing oil and gas is fractured vertically on completion and is produced at a constant bottom-hole pressure. The fracture has an infinite flow capacity, is of limited lateral extent and is bounded above and below by the impermeable strata defining the vertical extent of the reservoir. Results show that such a fractured reservoir can be represented by a reservoir of circular symmetry having very nearly the same production history. The well radius of this circular reservoir is about 1/4 the fracture length and is essentially the same as that obtained previously for a single fluid of constant compressibility. At the same value of cumulative oil production, gas-oil ratios of fractured reservoirs producing at constant terzinal pressure are larger than those of reservoirs having no fractures. This leads to more inefficient use of the reservoir energy in fractured wells and results in lower reservoir pressures for the same cumulative oil production. The reduction in operating life due to fracturing a reservoir is not as great as that for a slightly compressible fluid. This diflerence can be accounted for by the lower reservoir pressure in the fractured reservoir and its adverse effect on the average mobility and compressibility of the oil. As anticipated, the reduction in operating life increases czs the reservoir permeability decreases. The type of results presented in this report can be used to determine the economic attractiveness of fracture treatments per se, to setect the initial spacing to be used in developing a field, and to compare the relative merits of fracturing available wells and infill drilling. INTRODUCTION The effect of vertical fractures on a reservoir producing either an incompressible or a compressible liquid has already been discussed in the 1iterature.l,2 Those results indicate that the production history of such a reservoir is essentially the same as that of a circular reservoir having an effective well radius of approximately one-fourth the fracture length. The present work reports on the effect of a vertical fracture on a reservoir producing two compressible fluids —oil and gas—by solution gas drive. Because of the empirical nature of the PVT and relative permeability data used to obtain the performance of such reservoirs, results can only be obtained numerically and with the aid of high-speed computers. Since reservoirs lose their radial symmetry when fractured vertically, pressure and saturation can no longer be given only in terms of distance from the well. Two coordinates (such as x and y) must now be used to describe the pressure and saturation within the reservoir, and, since we are dealing with compressible fluids, time is also a variable. Thus the solution of a vertically fractured reservoir requires finding two unknowns (pressure and saturation) in two space variables (say x and y) and in time (t). Since no means are readily and generally available for solving such problems at the present time, we have used the results of previous work1,2 to approximate the effect of a vertical fracture on a reservoir producing both oil and gas by depletion. The purpose of the present wmk, then, is to investigate the possibility of using available numerical techniques (limited at the moment to one space variable) to study the two-space-variable flow behavior resulting from a vertical fracture. Results obtained in the course of this investigation are also reported and discussed. Input and output data of the numerical methods used are given in practical units: BOPD, feet, psi, cp, and md. Results are discussed fist in terms of specific reservoir and crude properties and geometries. Later, dimensionless parameters are introduced in order to extend results to different values of some of the reservoir and fracture properties. 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 reservoir has an impermeable, circular, cylindrical outer boundary of radius r,. 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). It is assumed that there is no pressure drop in the fracture due to fluid flow. 1 indicates the general three-dimensional geometry of the fractured reservoir. Gravity effects and the effects of differential depletion resulting from variations in hydrostatic head (pressure) will be neglected. Thus, the flow behavior in the fractured reservoir is described by the

By M. E. Shank, J. Wulff

At the present time, the designer of dies for metal powder pressing is handicapped by relative ignorance of stress distribution and frictional effects at the interior surface of the die. Unckell was the first to develop a method for the study of wall friction. He used three Brinell balls on which the die rested during pressing. The total frictional wall force was determined by the size of impression these balls left on a soft metal plate. Since the method does not give radial pressures, or distribution of such pressures, coefficients of friction could not be determined. Although Unckel measured density distribution, he was not able to determine radial or shear stresses. Shaler2 has proposed theoretical expressions for the stress and density distribution within cylindrical compacts during pressing, in accordance with the experimental results of Kamm, Steinberg, and Wulff.3 By application of Siebel's method,4 Kamm et a13 plotted stress trajectories for two compacts. From the stress trajectories they calculated coefficients of friction from point to point along the die wall. As pointed out by Shaler in the discussion of Ref. 3, these values are based on progressive point-to-point calculations on finite size grid squares across the compact. In the region of the die wall such calculated values may therefore have considerable cumulative error. The purpose of the present paper is to develop an experimental method by which the nonhydrostatic pressures and shears acting on the interior wall of a cylindrical die can be measured. Such measurements can then he correlated with existing data to aid in the explanation of the pressing process. The method used is based on the elastic: properties of the thick-walled tube used as the die. The principle of super-position of force systems on an elastic body is assumed to hold. Electric strain gauges were mounted in adjacent positions on the exterior die wall in order to get an exact measurement of the variation of tangential strain over the length of the die during pressing. While in this paper, measurements in terms of only tangential strains are considered, it is well to note that similar calculations may be set up for axial strains. The latter are not preferred, since they tend to be smaller than the tangential strains and therefore permit less sensitive measurements. Discussion in this work is restricted to compacts pressed from both ends, since the elastic deformation of the die is then more amenable to analysis. Before choosing the electric strain gauge method, a more direct line of attack was considered and discarded. The discarded idea was the insertion of a pressure gauge through a hole in the die wall.* The gauge would have been in the form of a small piston. If pressure were exerted against such a gauge, it would move outward along a radius of the die. One disadvantage of the scheme is its inability to measure shears along the die wall. Another more serious disadvantage is the disturbance caused by the device itself. It would serve to change the forces it was designed to measure. No matter how small the movement of the gauge, when pressure is applied a discontinuity would exist in the wall surface at that point. Due to the stress concentration caused by the hole, abnormal deflections of the die wall would occur around the gauge. During pressing, powder would be forced into the resulting depression. The depression would then become larger with increasing compacting pressure. Powder, not being a fluid, is capable of supporting shear. The ease with which it would flow into the die wall depression to further move the piston is an indication, not of the radial pressure at that point, but of the state of shear retarding the movement. Thus the "pressure" gauge is really a criterion of flowability, and of the capability of the powder to support shear. For these reasons, it was decided that the electric strain method, herein employed, was more reliable, if more indirect. The gauges and lead wires, mounted on the external die wall do not in any way affect the behavior of the metal powder or the die during pressing. Theory of the Method THE EFFECT OF RADIAL PRESSURE ON THE DIE WALL Effect of a Single Small Band of Hydrostatic Pressure Consider a die which is a thick-walled cylinder of outer radius R. and inner radius Ri. If over a small finite length L there is a normal pressure P, a tangential strain distribution at the outer wall results. This is shown schematically in Fig 1. The exact shape of the curve may he predicted by an extension of the theory of a semi-infinite beam on an elastic foundation.6 This

Jan 1, 1950

By Shiou-Chuan Sun

THE presence of an excessive amount of tough froth in the flotation of minerals, particularly coals, may create trouble in dewatering, filtering, and handling. Froth is also a nuisance in many chemical industries.' This paper presents a study on the destruction of extremely tough froths with intense high-frequency sound. The data indicate that sound waves can be employed for continuous atandsoundwavescan instantaneous defrothing. A powerful high-frequency siren was used in obtaining the data. Also tested was an ultrasonorator of the crystal type with a frequency range of 400, 700, 1000, and 1500 kc per sec and a maximum power output from its amplifier of 198 w. The results, not presented, indicate that as now designed this machine is not suitable for defrothing. Although the sound generators of the magnetostriction type2,3 and of the electromagnetic type'.' were not available, it is beelectromagneticlieved they are capable of producing the required sound intensity for defrothing. The use of ultrasonics for defrothing was suggested by Ross and McBain1 in 1944. Ramsey8 reported in 1948 that E. H. Rose mentioned a supersonic device that broke down flotation froth but with low capacity. The writer has not been able to find any published literature containing practical experiments. Theoretical Considerations The mechanism of defrothing by sound is attributed to the periodically collapsing force of the propagated sound waves and the induced resonant vibration of the bubbles. The collapse of froth is further facilitated by the sonic wind and the heat of the siren. Sound waves can exert a radiation pressure'," against any obstacle upon which they impinge. When a froth surface is subjected to the periodic puncturing of sound waves, the bubbles are broken. According to Rayleigh9 and Bergmann,12 the radiation pressure of sound, P, in dynes per sq cm is given as: P = 1/2 (r+1)i/v where r is the ratio of the specific heats of the medium through which sound is traveling and is equal to 1 on the basis of Boyle's law; i is the sound intensity in ergs per sec per sq cm, and v is the sound velocity in cm per sec. In this case, the accuracy of the formula is only approximate, because a perfect reflection can hardly result from a column of froth. In addition to the radiation pressure, the propagated sound waves cause the bubbles of the froth to have a resonant vibration.'" he vibratory motion of the bubbles causes collision and coalescence, thereby weakening if not breaking the bubble walls. Sonic wind and heat were also generated." The sonic wind can exert pressure on the froth surface, and the heat can evaporate the moisture content of the bubble walls as well as expand the enclosed air. Apparatus The defrothing apparatus, shown in Figs. 1 and 2, consists of a powerful high-frequency siren, a glass or stainless steel beaker of 2-liter capacity with 12.4 cm diam and 17.1 cm height, and a metal reflector. The beaker was placed 2 in. above the top point of the siren. The metal reflector was adjusted to reflect and focus the generated sound waves into the central part of the beaker. Fig. 2 shows the crystal probe microphone used to measure the acoustic intensity and the mandler bacteriological filter employed to introduce compressed air into the beaker for frothing. The apparatus was enclosed in a soundproof cabinet equipped with a glass window. The siren, shown in Fig. 3, consists of a rotor that interrupts the flow of air through the orifices in a stator. The rotor, a 6-in. diam disk with 100 equally spaced slots, is driven by a 2/3 hp, Dumore W2 motor at 133 rps. The frequency of the siren can be varied from 3 to 34 kc. The maximum chamber pressure is about 2 atm, yielding acoustic outputs of approximately 2 kw at an efficiency of about 20 pct. The siren itself is relatively small and can be operated in any orientation. A detailed description of the siren has been given by Allen and Rudnick.11 Collapse of Froth To study the sequence of the collapse of froth, the glass beaker was partially filled with 920 cc water, 100 g of —150 mesh bituminous coal, 0.3 cc petroleum light oil, 0.2 cc pine oil and 1.54 cc Pyrene foam compound. This mineral pulp was agitated for 5 min and then aerated through a mandler filter until the empty space of the beaker, approximately 9 cm high, was filled completely with min-

Jan 1, 1952

By Robin O. Williams

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

Jan 1, 1965

By Neils Hansen

The microstructure of dispersion-strengthened aluminum aluminum-oxide products containing from 0.2 to 4.7 wt pct of aluminum oxide has been examined by optical and transmission electron microscopy, and the flow stress has been determined at room temperature and at 400C by tensile testing. Products were examined as recrystallized and as high-temperature extruded, and the microstructures consisted of a fine dispersion of oxide particles in a matrix divided by respectively recrystallized grain boundaries and subgrain boundaries. The flow stress (0.2 pct offset) at room temperature of recrystallized dispersion strengthened aluminum aluminum-oxide products is the superposition of dispersion strengthening and grain boundary strengthening. This superposition has been found to be linear. The flow stress (a) can be related to the grain size (t) by the Petch equation: ing content of oxide and k is a constant independent of the oxide content. For extruded products a similar relation has been found by replacing the grain size by the subgrain size. The k-value is of the same order for the two types of structure, which shows that the subgrain boundaries are as effective slip barriers as grain boundaries. Tensile testing at 400C of re-crystallized and extruded products shows that oxide dispersion strengthening is very effective, whereas the strengthening effect of grain boundaries and subgrain boundaries is small. THE microstructure of dispersion-strengthened products consists of hard particles finely distributed in a metal matrix. The strengthening effect of the dispersed phase has been fairly well established,1 and it has been found that the size and volume fractions of the dispersed particles are important structural parameters. However, in many dispersion-strengthened products which have been worked and heat-treated during manufacture the matrix is divided into well-defined grains or sub-grains, which may also have a strengthening effect. A model of the matrix strengthening in dispersed products worked during manufacture has been proposed,2 introducing the energy of the structure as a strengthening factor, especially at low temperatures. A difficulty in this model is, however, to relate this (stored) energy to the structural parameters directly observable as for instance grain size. The strengthening effect of the matrix grain size after recrystallization has been in- vestigated for nickel-thoria (TD-Nickel) products3 and for copper aluminum-oxide products. Conclusive results were, however, not obtained as the grain size of TD-nickel was constant. 5 to II , after recrystallization at temperatures from 700 to 1200°C and as the copper products containing 5 to 1 wt pct of aluminum oxide could not be recrystallized even after severe cold reduction and heat treatment at 1050C. For aluminum aluminum-oxide products containing from 1 to 5 wt pct of aluminum oxide it has been shown that the tensile strength at room temperature decreases when an extruded product is cold-worked and recrystallized. The matrix in the extruded products is divided into well-defined subgrains of micron size, and as the grain size of the recrystallized products is about two orders of magnitude higher, it is obvious that grain boundary strengthening occurs. Preliminary results8 have indicated that the flow stress containing no grain boundaries, A is a constant and t is the subgrain size. At elevated temperatures the effect of boundaries is more complex; it has been shown11 that recrystallized products having an oxide content of about 3 wt pct are more creep resistant than extruded material in the temperature range 400° to 600°C, whereas on application of a higher strain rate the tensile flow stress (0.2 pct offset) is higher in extruded than in recrystallized aluminum—5 wt pct aluminum oxide products at temperatures from room temperature to 427°C (800), Finally it has been shown12 that the Brinell hardness at 350°C of extruded products having about the same content of aluminum oxide increases with decreasing grain size, determined by X-ray line-width measurements. The present study was undertaken to obtain a quantitative relationship between the tensile strength and the grain size of aluminum aluminum-oxide products in the recrystallized as well as in the extruded state. The tensile testing was performed at room temperature and at 400uC. The grain size of the recrystallized products was varied by changing the degree of cold-work preceding the recrystallization heat treatment. In extruded products grain (or subgrain) size variations were obtained by high-temperature heat treatment after extrusion. EXPERIMENTAL a) Materials. Aluminum aluminum-oxide products have been manufactured by consolidation of aluminum powder covered with a layer of aluminum oxide formed during powder manufacturing. The products examined were manufactured from atomized powder containing

Jan 1, 1970

By F. M. d’Heurle

Films of molybdenum have been prepared by sputtering onto oxidized silicon substrates. The resistivity. lattice parameter, orientation, and grain size were studied as a function of substrate temperature and substrate bias. Under normal sputtering conditions, the resistivity of the films was found to be quite high (600 x 10 ohm-crn). However, with the use of the negative substrate bias of 100 v and a substrate temperature of 350°C, films weve produced with a resistivity of ahout twice that of bulk molybdenum. The lattice parameters measured in a direction nornzal to the surface of the films weve found to be gveatev than the bulk value. This was interpreted as being at least partly due to the presence of compressive stresses. The effects of annealing in an Ar-H atmosphere were studied in terms of diffraction line width, lattice parameter, and resistivity. BECAUSE of its relatively low bulk resistivity (5.6 x 106 ohm-cm)' molybdenum is potentially interesting as a thin-film conductor in integrated circuits. An additional feature which makes it attractive for this purpose is its low coefficient of expansion (5.6 x KT6 per "c),' which is fairly well matched to that of silicon (3.2 x 10 per c). It is possible to deposit molybdenum films by evaporation but generally films produced in this manner have a high resistivity. In order to achieve resistivities close to bulk value, Holmwood and Glang found it necessary to operate in a vacuum of about 107 Torr and to maintain the substrates at 600 C during film deposition. Sputtered molybdenum films have been examined by Belser et a1.7 and, recently, by Glang et al.' This paper describes the results of an attempt to extend some of that work and examine the effects of annealing and getter sputtering on the physical and structural properties of the films produced. SPUTTERING APPARATUS AND PROCEDURE The apparatus used for most of the film sputtering work described here consisted of two "fingers" serving as anode and cathode, respectively, which were mounted within an 18-in.-diam glass chamber. A liquid nitrogen-trapped 6-in. diffusion-pump system was used to achieve a vacuum of about 1 x 107 Torr within the chamber prior to sputtering. The essential features of the equipment are shown in Fig. 1. Cathode and anode fingers are stainless-steel tubes isolated from the top and bottom plates by Teflon collars. In order to limit the discharge to the space between anode and cathode, each finger is surrounded by an aluminum hield, at ground potential, having an internal diameter 18 in. larger than the outside diameter of the finger. The cathode and anode fingers are 6 and 4 in. in diam, respectively. A 116-in.-thick sheet of molybdenum is brazed with a 10 pct Pd, 58 pct Ag, 32 pct Cu alloy to a copper disc which is mounted by means of screws and a large 0 ring onto the lower end of the cathode finger. The disc is cooled during sputtering by water circulation inside the finger. The use of several feet of plastic tubing for the water input and outputg reduces leakage to ground to less than 1 ma when the cathode potential is raised to 5 kv. The upper end of the anode finger is terminated by a brazed-on copper block. A variety of specimen holders can be easily mounted on the upper face of this block. Substrate heating or cooling is achieved by use of an appropriate unit attached to the lower face of the same block. Heating is achieved by means of cartridge-type heaters and cooling by copper coils fed with forming gas under pressure. The inner chamber of the specimen finger constitutes a small vacuum chamber of its own which is evacuated by an auxiliary mechanical pump in order to limit heating element oxidation and heat transfer by convection currents. An advantage of the finger arrangement is the absence of cooling and heating coils and wires within the main chamber. The stain less-steel shutter is useful to establish a discharge for cleaning the cathode at the beginning of each sputtering run. Water cooling of the shutter reduces heating and the out-gassing of impurities which might condense on the nearby substrates. Unless otherwise specified, the substrates used in these experiments were 1-in.-diam oxidized silicon wafe:s, 0.007 in. thick, having an oxide thickness of 6000A. The substrate holders were large copper discs onto the surface of which a number of molybdenum discs, 116 in. thick and 78 in. in diam, were brazed. The wafers were clamped to the molybdenum discs

Jan 1, 1967

By A. W. Adamson

A model for transport processes in liquid mixtures is discussed which supposes that the elementary act involves a position exchange between two species and that the exchange is so confined by the solvent cage as to occur nearly isosterically. The rate-determining step, thus, is likened to a bi-molecular reaction and is so treated, using absolute rate theory. The cage model has been applied to diffusion, thermal diffusion, sedimentation and viscosity, but only the first and last of these phenomena are emphasized in the present paper. The model leads to semi-empirical relationships between the absolute value for a digusion coefficient and the activation energy for diffusion, between mutual and self-diffusion coefficients and for the variation of the viscosity of a binary mixture with composition. These are discussed in relation to experimental data for various systems, including hydrocarbon mixtures. It is shown that the proposed viscosity equation and seven other commonly used ones all may be regarded as special cases of a single general relationship; a brief critical analysis is made of the basis of selection of one or the other for data fitting or interpolation. INTRODUCTION AND GENERAL THEORY The present paper covers a brief discussion of a cage model for transport processes in liquid mixtures and how this model may be useful in treating the diffusional behavior and the viscosity of such systems. Since diffusion requires the more detailed treatment, it will be taken up first, and the model then applied to viscosity. There are two types of diffusion coefficients that may be measured experimentally, apart from thermal diffusion quantities. The first is the mutual or binary diffusion coefficient, D which may be defined in terms of Fick's first law. This states that the permeation, or flux P, is proportional to the concentration gradient. In the usual experiment, P is measured relative to a frame of reference fixed with respect to the medium (e.g., the diaphragm in a diffusion cell); as a consequence, the same value of D is obtained regardless of whether P and C refer to Component 1 or to Component 2; i.e., there is only one independent mutual diffusion coefficient for a binary system. In addition to D there will be various self-diffusion coefficients. defined in terms of the gradient in labelled species i and its permeation in an otherwise uniform medium. The thermodynamic approach to mutual diffusion supposes that the actual driving force is the gradient of the chemical potential, i.e., that In the case of a dilute solution of solute, Eqs. 1 and 3 lead to the Einstein equation, If the solution is ideal and the friction coefficient is taken to be then the familiar Stokes- Einstein equation results. Mutual and self-diffusion coefficients can not be related on general thermodynamic grounds; it is necessary to invoke some additional assumptions, i.e., a model; several such have been proposed. Hartley and Crank' supposed the existence of separate, intrinsic diffusion coefficients (Dl and D2) for each component, essentially corresponding to the two self-diffusion coefficients. The two flows can not be independent, however, but must be coupled through the usual restriction that there be no net volume flow. For an ideal solution. one then obtains' Glasstone, et al' treated diffusion in terms of absolute rate theory, but their approach otherwise resembled the previously mentioned one in that each species was considered to move with respect to the general medium in a manner determined by its individual jump distance and specific rate constant. For other than dilute solutions, a coupling of flows leading to an equation such as Eq. 6 would again be present. However, as required by Eq. 6, one does expect that the self-diffusion coefficient for the solute and the mutual-diffusion coefficient for the system become identical at infinite dilution. Lamm4 recognized that there should be three distinctive interactions in a two-component system-1-1, 1-2 and 2-2 — and, therefore, proposed three rather than two fundamental friction coefficients. Mutual diffusion resulted from 1-2 interactions only, and self-diffusion resulted from 1-2 plus either 1-1 or 2-2 interactions. Again, a collective coupling between all motions was imposed to meet the condition of no net volume flow. Laity' has shown how to convert the Onsager equations to a form very similar to Lamm's. Cage Model For Diffusion Work in this laboratory on diffusion in aqueous sucrose solutions made it apparent that three, rather than two, interactions were indeed needed," but considera-

By Hollis L. Caswell

The present status of integrated circuits utilizing. superconductive switching. elements is reviewed with special attention given to fabrication techniques, methods for interconnecting completed circuits, and refrigeration requirements. Cryoelectronics has been largely an "inte- grated-circuit" technology since its conception because the switching speed of superconductive devices is attractive only when these devices are fabricated with thin-film techniques. It is true that cryotron circuits can be constructed from wires of appropriate materials (as indeed was done by Dudley Buck 1 in his early investigations) but these circuits will switch in times characteristic of milliseconds whereas similar circuits fabricated by thin-film methods have potential switching times of nanoseconds. Furthermore, cryo-electronic devices such as the cryotron lend themselves readily to fabrication by thin-film techniques since these components may be made from polycrys-talline thin films and are relatively insensitive to the presence of impurities (as measured by semiconductor standards). Therefore, during the past decade considerable effort has been devoted to developing techniques for batch fabricating circuit arrays containing superconductive switching elements. Technology had developed to the point several years ago that fabrication of cryoelectronic arrays containing up to one hundred devices was rather straightforward. However, larger arrays containing between lo4 and 106 components which are required for commercial development of cryoelectronics still pose very severe yield problems. Thus in a sense cryoelectronics found itself in 1962 at the point semiconductor technology finds itself today; namely, individual devices and small groups of integrated devices could be fabricated with acceptable yield and the outlook for building larger integrated-circuit arrays was bright. Unfortunately, problems associated largely with yield have made fabrication of these larger arrays difficult. Unlike semiconductor technology, cryoelectronics had to solve the problems of large-scale integration before it could become economically attractive. This has proven to be a sizable burden to bear. Since several reviews exist on superconductivity,2 superconductive devices,3 and cryoelectronic technology, no attempt will be made in this paper to summarize these areas. Instead a few specific topics will be dealt with in more detail. First, a brief description is given of selected superconducting switching and storage devices with special attention to several metallurgical techniques which improve the performance of these devices. Second, techniques used to fabricate cryoelectronic devices are described with emphasis on problems affecting yield. Third, techniques for interconnecting a number of cryoelectronic planes are described. And last, refrigeration of cryoelectronic components is discussed briefly since the low operating temperature of superconductive devices is an important consideration in this technology. SUPERCONDUCTING STORAGE AND SWITCHING DEVICES The basic superconductive switching device is the thin-film cryotron. The geometry of this device is attractively simple, since it involves only the intersection of two lines that are electrically insulated from each other. The switching element (gate) and control element (control) of a crossed-film cryotron are arranged as illustrated in Fig. 1. The material for the gate is selected to permit the gate to be switched from the superconducting to the normal (resistive) state by the application of a control current. Tin, which has a critical temperature (T,) of 3.7°K, is commonly used for the gate and the cryotron is operated at a temperature just below T, (for example, 3.5°K). The control material (normally lead, with T, = 7.2°K) is chosen so that the control is never driven normal during circuit operation. To improve cryotron operation, a ground plane, also of lead, is placed under all of the circuitry to act as a diamagnetic shield and improve the current-density uniformity across the width of various thin-film elements. Normally, line widths vary from 0.005 to ^ 0.020 in. and film thicknesses from 300 to 10,000A, although new fabrication techniques make narrower lines feasible. In fabricating cryotrons it is important that the edges of the gate elements be geometrically sharp to avoid undesirable switching characteristics associated with a thinner edge region, Fig. 2. One technique which has been used extensively to form patterns consists of placing a physical mask containing the film pattern between the evaporation source and the substrate and depositing through the mask. Film strips formed in this manner possess a penumbra at the film edges due to shadowing of the evapor-ant under the mask. Several techniques have been proposed for minimizing effects due to this penumbra. One of the more promising metallurgical techniques

Jan 1, 1967

By E. Eichen, G. S. Ansell. L. R. Sefton

A Technique for the Preparation of Thin Films of Two-Phase Alloys Suitable for Use in Transmission Electron Microscopy In order to obtain foil sufficiently thin to permit transmission electron microscopy of two-phase alloys, particularly alloys in which the second phase is present as a very finely dispersed second phase, a technique has been developed to permit the thinning of these materials without either the large degree of pitting present in the electrochemical or chemical thinning process or the slicing of the second phase particles accompanying the mi-crotoming process. The metal or alloy to be thinned is first ground to initial sheet thickness of approximately 0.002 in. by hand grinding on 4/0 emery paper. This starting sheet of material is then reduced in cross section by an air-abrasion technique using an S. S. White dental abrasive unit. This unit consists of a supply of compressed CO, which provides a stream of high velocity gas impinging on the metal surface. Incorporated in the machine is a hopper-type of device in which fine abrasive particles of alumina are introduced into this high velocity gas stream. The abrasive particles are irregular in shape, approximately 30 across. If any particles are introduced into the sample during thinning, they should be readily observable. Impingement of this gas abrasive mixture on the metal surface causes rapid removal of the metal in a rather uniform manner. The rates of metal removal may be adjusted by both regulating the gas pressure and amount of abrasive material introduced into the gas stream. By moving the gas stream abrasive mixture slowly about the specimen or sheet surface, a large area of the sheet may be uniformly reduced to a thickness which is then suitable for use in transmission electron microscopy. Attendant with this metal removal is a result in cold working of the metal surface due to the impingement of the abrasive particles. The resultant foils, therefore, if made from materials which are desirable to view in the annealed condition, have to be subsequently heat treated. If the alloy from which these foils have been made has been initially heavily worked or deformed during manufacture, then no subsequent heat treatment is, necessary as in the case of the SAP-type alloys. Following this initial abrasive thinning technique, the material is then given a final electropolish which is just sufficient to remove any oxide which has formed on the metal surface during the abrasion technique. Care has to be exercised in this technique to make sure that none of the attendant problems associated with the chemical or electrochemical removal process are encountered. An electrolyte solution which has been found suitable for aluminum-aluminum oxide SAP materials is one containing 78 ml of perchloric acid, 120 ml of distilled water, 700 ml of ethanol, and 100 ml of butylcellosolve. The current density used is 1 amp per sq cm and a polishing time of 3 sec on each side of the foil is used. Extreme care must be used in choosing the current density and time to ensure that just the oxide film on the specimen is removed. Times or current density which are greater than this will lead to perforation of the film due to the removal of the A1,0, particles or a complete loss of the film due to complete solution of the matrix. Once the oxide film is removed the specimens must be kept in alcohol until placed in the microscope to prevent any further oxidation. Fig. 1 shows -an electron micrograph at X28,000 of a thin foil of an aluminum-aluminum oxide SAP-type alloy prepared using this technique. The alloy consists of a matrix of commercial purity aluminum containing a very finely dispersed second phase of aluminum oxide flakes. In the electron micrograph, the dark patches are the aluminum oxide particles. The lighter portions of the micrograph are the aluminum matrix. Within the aluminum matrix, the dislocation structure of the alloy is clearly visible. This structure consists of intersecting twin boundaries and single dislocations. No abrasive particles were observed in the samples. In addition, thin films of commercial purity aluminum were prepared using this method and observed in the electron microscope. Here also, there was no evidence of abrasive particles introduced in the samples. The technique presented provides a satisfactory method for the production of thin films of aluminum-aluminum oxide SAP-type alloys for use in transmission electron microscopy where none of the usual thinning techniques appears to be satisfactory.

Jan 1, 1962

By W. Hurst, R. E. Leeser, W. C. Goodson

Three aspects of gas deliverability are presented in this paper. The first treats with the gas deliverability or availability of a normal depletion-type dry gas field. Such encompasses not only the period of stabilized constant rate, but more so, the "tailings" when a fixed abandonment pressure is reached and the rate by necessity must decline. A comprehensive work plot is offered, developed from mathematics herein included, that removes the triai-and-errnr computations that attended such undertakings in the past. The second part treats with the discount factor of the open flow potential constant from what is observed initially in testing a gas well to what is evidenced when stabilization is reached. This prevails in tight formations, such as the Kansas Hugoton field which is offered as the example. The means of establishing this factor are pressure build-up curves which, as sustained by analytical deductions, reproduce this entire period of transient flow under conditions of a constant rate influx. Finally, what is offered in this paper is the deliverability performance of an exceedingly rich gas condensate field producing from a tight formation. The example shown is the Knox Bromide field in Oklahoma, producing from the Bromide formations. The results are ominous, showing early reduction in permeability to gas pow, due to the retrograde condensate forming in the pore space, with the attending early logging-up of these wells. The analytics of lowered permeability are incorporated in the gas deliverability formula along with the PVT data that gives the increased condensate liquid saturation as the gas flows to the well bore. This paper would not be complete without a critique oflered at the end. With the many gas wells now in production and those that have completed their life, there has been no factual information collected by any source as to what constitutes that permeability range where a gas well would be unimpaired in its gas deliverability by the presence of rich condensate content, and the lowered range where such would be harmful. This question confronts all producers. INTRODUCTION Various aspects of gas deliverability are presented in this paper that includes depletion-type reservoirs, deteriora- tion factor of the gas deliverability constant, and the performance of a rich gas condensate reservoir producing from a tight sand. With respect to the presentation of gas deliverability and its tailings for depletion-type gas reservoirs, one notes that this is essentially the information offered by every gas transmission company and producer appearing before the Federal Power Commission for Letters of Conveyance in the dedication of reserves. In the ordinary procedure, as many engage upon this study, trial-and-error calculations are included, particularly as apply to the tailings. For many years one of the writers has employed mathematical analyses to perform this step and avoid the complexities so associated. In the preparation of this paper these analyses have been amplified to include any slope n for the open flow potential relationship for which the tailings can be determined from Fig. 1. With reference to the deterioration or discount factor of the open flow potential constant as such occurs in the gas deliverability formula, this for the most part has been an unexplored subject. Although the issue first appeared in the Kansas Hugoton field, where such was surmised but only recently resolved, this situation of a deterioration of the gas deliverability constant can occur wherever dry gas production from a tight sand is encountered. The first concerted attacks upon this problem were the presentations of Hurst' and Goodson' before the Kansas Corporation Commission to show that transient fluid flow and unsteady-state flow formulas prevailed. This was amplified later before the Federal Power Commission3 to show that this deterioration factor could be identified from pressure build-up curves. This has been reported by McMahon.4 Its importance to the industry merits the review of these essential features in completing the program on the aspects of gas deliverability. Finally, as illustrated here, for a low permeability formation such as the Knox Bromide field where the gas is rich, representing some 165 bbl of condensate per MMcf of effluent gas, the gas deliverability can be of limited extent in the life of the field, leaving substantial amounts of condensate and gas unrecovered. In cases such as this, gas cycling is mandatory. This is particularly revealed by the fluid mechanics introduced here, employing factual field as well as laboratory data, to show this-restriction upon gas deliverability. PRESSURE DEPLETION What will now be offered is the study of gas deliver-

By J. U. Messenger

A technique is described whereby the resistance of an emudian to breaking can be quantitatively determined. Produced ailfield emulsions are usually the water-in-oil type and, accordingly, do not conduct an electrical current. However, there is a threshold of A-C voltage pressure above which an emulsion will break and current will flow. The more stable an emulsion, the higher the required voltage. A Fann Emulsion Tester, modified so that low voltages (0 to 10 v) can be accurately measured, is suitable. This technique has application in evaluating the effect of a demuksifier on the stability of an emulsion. Emulsions can, in essence, be titrated with demulsifiers by adding a quuntity of demulsifier, stirring, and measuring the voltage required to cause current to flow. Any synergistic effect of two or more materials added simultaneously can be followed accurately. A demulsifier that significantly lowers the threshold voltage (from 100 to 400 v to 0 to 10 v for the emulsions in this study) is effective and can cause the enlulsion to break. A demulsifier that will bring about this drop in the threshold voltage at low concentration ir very desirable. The technique is also well adapted for rapidly screening demulsifiers. INTRODUCTION Stable emulsions in produced reservoir fluids resulting from certain well stimulation and completion procedures are common problems. The use of suitable demulsifiers can often mitigate these difficulties. At the present time, a rapid and efficient method for selecting satisfactory demulsifiers is not available. It is badly needed. Reliance is now placed primarily on trial-and-error procedures. A new test method has been developed which permits a more rapid and precise selection of demulsifiers. It involves measuring the electrical stability potential of an emulsion before and after a demulsifier has been added. This paper describes this method and shows where it should have application in field emulsion problems. NATURE OF OILFIELD EMULSIONS Two immiscible components must be present for an emuhion to form; we are concerned here with crude oil and water. An emulsifier must be present for tin emulsion to be stable. J Emulsifiers can be substances which are soluble in oil and /or mter and which lower interfacial tension. They can be colloidal solids such as bentonite, carbon, graphite, or asphalt which collect at the interface and are preferentially wet by one of these phases. Unrefined crude oils can contain both types of emulsifiers, A popular theory is that, of the two phases in an emulsion, the dispersed phase will be the one contributing most to the interfacial tension.' Usually this phase contains the least amount of emulsifier. The stability of a water-in-oil emulsion is affected by the fol1owing: l) viscosity; (2) particle or droplet size; (3) interfacial tension between the phases; (4) phase-volume ratios; and (5) the difference in density between the phases. A stable emulsion is usually characterized by high-viscosity, small droplets, low interfacial tensions, small differences in density between its phases, and slow separatian of the phases. It also has low conductivity (high electrical stability potential). Water-in-oil and oil-in-water emulsions"' are both common; however, oil field emulsions are predominantly water-in-oil emulsions. The emulsions which commonly occur during oompletion and stimulation operations contain a combination of several of the following: acids, fracturing fluids (oil, water, acid), and formation water and oil. Produced emulsions usually contain formation water and oil. Emulsions form in oil wells because oil and water are mixed together at a high rate of shear in the presence of a naturally occurring or unavoidably produced emulsifier. During the completion and stimulation of productive zones, and while formation fluids are being produced, oil and water are very often commingled. These mixtures are formed into emulsions by agitation which occurs when the fluids are pumped from the surface into the matrix of the formation or produced through the formation to the surface. Restrictions to flow (such as perforations, pumps, and chokes)".'" increase the level of agitation; tight emulsions are more likely to form under these conditions. Often an emulsified droplet is an emulsion itself.'" Therefore, emulsion-breaking problems can be quite complex. The complexity can be even greater if a third phase (gas) is included. Demulsifiers operate by tending to reverse the form of the emulsion. During this process, droplets of water become bigger, viscosity is lowered, color becomes darker, separation of the phases faster and electrical stability potential approaches zero. Any of these effects could be followed as a means of determining emulsion stability. However, electrical stability potential is the most reproducible and most easily measured parameter for following the stability of a water-in-oil emulsion.

Jan 1, 1966

By J. Chipman, N. J. Grant, B. M. Shields

The vacuum tin-fusion method of analysis for hydrogen, developed by Carney, Chipman, and Grant, has been modified to permit the analysis of the evolved gases for hydrogen by means of a thermal conductivity cell. A properly prepared sample can be analyzed in 10 min with a probable error of ±0.12 ppm. A study of various methods for storage of hydrogen samples shows that samples can be safely held in a dry ice-acetone bath as long as six days. Storage in liquid nitrogen is necessary for samples to be held one week or more. HE vacuum tin-fusion method, as developed by I- Carney, Chipman and Grant,' is the only analytical procedure which has shown promise of being fast enough for use in the control of hydrogen during steelmaking. It was felt that further simplification and faster speed of operation could be effected by the use of thermal conductivity measurements for analysis of the gases evolved in the tin-fusion method. The application of conductivity measurements to the tin-fusion method is possible because: 1—the evolved gas is essentially a mixture of hydrogen, nitrogen and carbon monoxide with a hydrogen content usually over 50 pct, 2—the evolved gas is collected at a relatively low pressure, and 3— the thermal conductivities of CO and N2 are practically identical while that of hydrogen is very much greater. The major part of this research program was devoted to the construction and calibration of a vacuum tin-fusion apparatus which analyzes the evolved gases for hydrogen by means of a thermal conductivity cell. The second phase of the problem was associated with the development of a procedure for storage of samples prior to analysis. With the rapid quenching method for hydrogen sampling,' which seems to be the most practical for steel mill use, it is necessary that the samples be stored safely during the interval between sampling and analysis if the hydrogen content of the molten metal is to be maintained in the supersaturated solid samples. The thermal conductivity bridge has been used for a number of years in the analysis of certain gas mixtures. An elementary discussion of the theory and practice of gas analysis by thermal conductivity measurements is given by Minter.3 A more comprehensive discussion of the theory and of the various measuring circuits is presented by Daynes.' A complete knowledge of the theory and properties of the thermal conductivity of gases and gaseous mixtures can be gained by a study of the standard textbooks on the kinetic theory of gases."' The existing data on the thermal conductivity of single gases are reviewed by Hawkins: that for a number of binary gas mixtures by Daynes' and Lindsay." The thermal conductivity method may be applied to the determination of the composition of a binary mixture if: 1—the thermal conductivity of the mixture varies monotonically with composition, and 2— the two gases have measurably different thermal conductivities. The greater the difference between the two gases, the greater the sensitivity of the method.10 he method is applicable to the analysis of multicomponent mixtures when all of the gases in the mixture except one have nearly the same thermal conductivity. Fortunately, the mixture of hydrogen, nitrogen, and carbon monoxide evolved by the tin-fusion analysis' falls in this latter classification. The thermal conductivities of nitrogen and carbon monoxide are practically equal; and the thermal conductivity of hydrogen is approximately seven times that of the other two. Therefore, the thermal conductivity of a gaseous mixture of hydrogen, nitrogen, and carbon monoxide at known temperature and pressure can be related directly to the percentage of hydrogen in the mixture by suitable calibration. Usually the thermal conductivity of a mixture of gases is measured at atmospheric pressure where the thermal conductivity is independent of pressure over a wide pressure range. At very low pressures (below 1 mm Hg), the thermal conductivity of gases varies with the pressure. This phenomenon has been utilized in the Pirani vacuum gage for the measurement of pressures in the range of 10" to 10-0 mm of mercury.= Very little has been published concerning the variation of thermal conductivity with pressure at intermediate pressures between 1 mm Hg and 1 atm. However, preliminary measurements indicated that the thermal conductivities did vary with pressure over the range of pressures (up to 10 mm Hg) at which gases are delivered from the vacuum pump. Therefore, the calibration of the thermal conductivity cell had to be planned to include the effects of both gas composition and pressure. Such a calibration chart is shown in Fig. 4. Most industrial applications of the thermal conductivity method of gas analysis have used a compensated Wheatstone bridge circuit containing two

Jan 1, 1954

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

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

Jan 1, 1970

By W. R. Purcell

An apparatus is described whereby capillary pressure curves for porous media may be determined by a technique that involves forcing mercury under pressure into the evacuated pores of solids. The data so obtained are compared with capillary pressure curves determined by the porous diaphragm method, and the advantages of the mercury injection method are stated. Based upon a simplified working hypothesis, an equation is derived to show the relationship of the permeability of a porous medium to its porosity and capillary pressure curve, and experimental data are presented to support its validity. A procedure is outlined whereby an estimate of the permeability of drill cuttings may be made with sufficient acuracy to meet most engineering requirements. INTRODUCTION The nature of capillary pressures and the role they play in reservoir behavior have been lucidly discussed by Lev-rett', Hassler, Brunner, and Deah12, and others. As a result of these publications the value of determining capillary pressure curves for cores has come to be generally recognized within the oil industry. While considerable attention has been directed toward the subject in an effort to provide a reliable method of estimating percentages of connate water, it has been recognized that capillary pressure data may prove of value in other equally important applications. This paper describes a method and procedure for determining capillary pressure curves for porous media wherein mercury is forced under pressure into the evacuated pores of the solids. The pressure-volume relationships ob- tained are reasonably similar to capillary pressure curves determined by the generally accepted porous diaphragm method. The advantages of the method lie in the rapidity with which the experimental data can be obtained and in the fact that small, irregularly shaped samples, e.g., drill cuttings, can be handled in the same manner as larger pieces of regular shape such as cores or permeability plugs. Based upon a simplified working hypothesis, a theoretical equation will be derived which relates the capillary pressure curve to the porosity and permeability of a porous solid, and experimental data will be presented to support its validity. This relationship aplied to capillary pressure data obtained for drill cuttings by the procedure described provides a means for predicting the permeability of drill cuttings. METHODS FOR DETERMINING CAPILLARY PRESSURES Several techniques have so far been employed in determining capillary pressure curves and these fall into two principal categories: (1) Liquid is removed from, or imbibed by, the core through the medium of a high displacement pressure porous diaphragm (2) Liquid is removed from the core which is subjected to high centrifugal forces in a centrifuge4,'. There are? however, certain limitations inherent in both methods. The greatest capillary pressure which can be observed by method (I), above, is determined by the maximum displacement pressure procurable in a permeable diaphragm which at the present time appears to be less than 100 psi. An even more serious limitation of the diaphragm method is imposed hy the fact that several days may be required to reach saturation equilibrium at a given pressure; hence, the time re- quired to obtain a well-defined curve may be measured in terms of weeks. Furthermore, to date, no suitable technique for handling relatively small, irregularly shaped pieces of rock, such as drill cuttings, has been reported and, therefore, measurements must be made, in general, on cores, or portions thereof. The centrifuge method offers the distinct advantage over the porous diaphragm method of arriving at saturation equilibrium in a relatively short time by virtue of the elimination of the transfer medium for the liquid. The calculation of capillary pressures from centrifuge speeds is somewhat tediousa, however, and the equipment required is fairly elaborate. While there exists the possibility that this method might be adaptable to the determination of the capillary pressures of cuttings, this particular ramification has not been investigated, as far as is known. In view of the limitations of the two principal methods for determining capillary pressures, the apparatus described in the following sections has been devised in order that difficulties previously encountered might be circumvented. MERCURY INJECTION METHOD FOR DETERMINING CAPILLARY PRESSURES Theory The methods described above for determining capillary pressures are characterized by the fact that one of the fluids present within the pore spaces of the solid is a liquid which "wets" the solid, i.e., the contact angle which the liquid forms against the solid is less than 90" as measured through that phase. For these "wetting" liquids the action of surface forces is such that the fluid spontaneously fills the voids within the solid. These forces likewise oppose the withdrawal of the fluid from the pores of the solid.

Jan 1, 1949

By W. R. Purcell

An apparatus is described whereby capillary pressure curves for porous media may be determined by a technique that involves forcing mercury under pressure into the evacuated pores of solids. The data so obtained are compared with capillary pressure curves determined by the porous diaphragm method, and the advantages of the mercury injection method are stated. Based upon a simplified working hypothesis, an equation is derived to show the relationship of the permeability of a porous medium to its porosity and capillary pressure curve, and experimental data are presented to support its validity. A procedure is outlined whereby an estimate of the permeability of drill cuttings may be made with sufficient acuracy to meet most engineering requirements. INTRODUCTION The nature of capillary pressures and the role they play in reservoir behavior have been lucidly discussed by Lev-rett', Hassler, Brunner, and Deah12, and others. As a result of these publications the value of determining capillary pressure curves for cores has come to be generally recognized within the oil industry. While considerable attention has been directed toward the subject in an effort to provide a reliable method of estimating percentages of connate water, it has been recognized that capillary pressure data may prove of value in other equally important applications. This paper describes a method and procedure for determining capillary pressure curves for porous media wherein mercury is forced under pressure into the evacuated pores of the solids. The pressure-volume relationships ob- tained are reasonably similar to capillary pressure curves determined by the generally accepted porous diaphragm method. The advantages of the method lie in the rapidity with which the experimental data can be obtained and in the fact that small, irregularly shaped samples, e.g., drill cuttings, can be handled in the same manner as larger pieces of regular shape such as cores or permeability plugs. Based upon a simplified working hypothesis, a theoretical equation will be derived which relates the capillary pressure curve to the porosity and permeability of a porous solid, and experimental data will be presented to support its validity. This relationship aplied to capillary pressure data obtained for drill cuttings by the procedure described provides a means for predicting the permeability of drill cuttings. METHODS FOR DETERMINING CAPILLARY PRESSURES Several techniques have so far been employed in determining capillary pressure curves and these fall into two principal categories: (1) Liquid is removed from, or imbibed by, the core through the medium of a high displacement pressure porous diaphragm (2) Liquid is removed from the core which is subjected to high centrifugal forces in a centrifuge4,'. There are? however, certain limitations inherent in both methods. The greatest capillary pressure which can be observed by method (I), above, is determined by the maximum displacement pressure procurable in a permeable diaphragm which at the present time appears to be less than 100 psi. An even more serious limitation of the diaphragm method is imposed hy the fact that several days may be required to reach saturation equilibrium at a given pressure; hence, the time re- quired to obtain a well-defined curve may be measured in terms of weeks. Furthermore, to date, no suitable technique for handling relatively small, irregularly shaped pieces of rock, such as drill cuttings, has been reported and, therefore, measurements must be made, in general, on cores, or portions thereof. The centrifuge method offers the distinct advantage over the porous diaphragm method of arriving at saturation equilibrium in a relatively short time by virtue of the elimination of the transfer medium for the liquid. The calculation of capillary pressures from centrifuge speeds is somewhat tediousa, however, and the equipment required is fairly elaborate. While there exists the possibility that this method might be adaptable to the determination of the capillary pressures of cuttings, this particular ramification has not been investigated, as far as is known. In view of the limitations of the two principal methods for determining capillary pressures, the apparatus described in the following sections has been devised in order that difficulties previously encountered might be circumvented. MERCURY INJECTION METHOD FOR DETERMINING CAPILLARY PRESSURES Theory The methods described above for determining capillary pressures are characterized by the fact that one of the fluids present within the pore spaces of the solid is a liquid which "wets" the solid, i.e., the contact angle which the liquid forms against the solid is less than 90" as measured through that phase. For these "wetting" liquids the action of surface forces is such that the fluid spontaneously fills the voids within the solid. These forces likewise oppose the withdrawal of the fluid from the pores of the solid.

Jan 1, 1949

By J. B. Cohen, O. F. Kimball

Aging of specimens of Au-20 at. pct Ni and Au-40 at. pct Ni has been followed with electron microscopy, resistance, magnetic measurements, and X-ray diffraction. The periodic structures observed by electron microscopy in Au-Ni alloys by Fisher and Embury are shown to be due to nickel-rich prepre-cipitates. Assuming this alloy system undergoes spinodal decomposition, our observations and available data on thermodynamic properties and elastic constants indicate that the gradient energy coefficient in this system is of the order of 10-5 ergs cm-1 . Shifts of X-ray peaks in this system have been found to be due primarily to residual stresses, not faulting. A large effect of applied stress during aging was detected by X-rays and electron microscopy for the Au-20 at. pct Ni alloy. The presence of a preprecipitation stage during aging of gold-rich Au-Ni alloys for short times below 300°C has been indicated by several previous investigations. Changes in electrical resistivity,1"3 elastic modulus,4 superparamagnetisrn,5 positions of X-ray diffraction peaks,6 and hardness7 have been studied. Electron diffraction and microscopic studies by Fukano8 and Fisher and Embury9 respectively showed that the precipitation process resulted in the formation of periodic composition fluctuations—perhaps spinodal decomposition. MOSS10 detected a satellite with X-rays in quenched specimens, but the wavelength was much smallerthan that found by Fukano. Recently Woodilla and Averbac11 detected a modulation similar to that found by Moss, but with electron diffraction. They suggested that the microstructures observed by Fisher and Embury may be simply due to a Moir'e effect from a back deposit of gold during electropolishing. Satellites in electron diffraction can be enhanced by small bends or tilts of a foil. Even if this were not the case, asymmetric satellites, as are found in Au-Ni are difficult to analyze. They are principally influenced by the large differences in atomic size in this system, and obtaining information on the associated composition fluctuations depends on assumptions as to whether or not there are isolated clusters or composition "waves", and in the latter case whether or not the entire specimen has decomposed. Assuming that decomposition was complete in alloys aged at low temperatures, Woodilla and Averbach suggested that the satellites indicated a composition which varied about the average value by only± 5 pct, much less than that predicted from the phase limits, Fig. 1; they attributed this small variation to a cessation of the reaction caused by the loss of mobile vacancies. In Ref. 12 it was shown that vacancies are trapped by the preprecipitate in this system. Furthermore, after reversion, aging proceeds at a much slower rate than immediately after a quench from the solution temperature,5'11,13 indicating the importance of vacancies in the kinetics. In some senses Au-Ni is a favorable system for studying preprecipitation and reversion mechanisms because of the simplicity of the equilibrium diagram14 and the availability of good thermodynamic15 and diffusion data.16,l7 The calculations necessary to relate the decomposition process to classical nucleation and growth or to spinodal decomposition, as well as those necessary for understanding reversion mechanisms, can, therefore, be made. As an example, the position of the chemical spinodal, calculated from the thermo-dynamic data of Sellars and Maak15 with the method described by Rundman and Hilliard,18 is included in Fig. 1. It is compared to the chemical spinodal calculated from the equilibrium miscibility gap by Cook and Hilliard's method.l9 Included is the coherent spinodal suggested by Golding and MOSS," using measured elastic constants for their calculations and the same ther-modynamic data. The spinodal estimated by Woodilla and Averbach,11 by determining the temperature limits for the observation of the satellites discussed above, is also given. Actually this last curve is open to question on the following grounds: As a phase boundary is approached the spacing of particles or the wavelength

Jan 1, 1970

By D. H. Yardley

GEOCHEMICAL investigation of trace elements in surface materials was begun near Ely, Minn., in 1953 along the basal contact of Duluth gabbro with Giants Range granite (Fig. 1). This article presents data on the distribution of copper and nickel in till and in stream sediments in the area and proposes an explanation for the types of distribution found. The Duluth gabbro, one of the world's largest basic intrusives, intrudes rocks which range in age from Keewatin to middle Keweenawan. Within the test area the gabbro is in contact with granite except for short sections where it is in contact with remnants of iron formation. Sulfide mineralization occurs within the gabbro, near and parallel to the basal contact for a distance of several miles. Schwartz and Davidson' have described the geologic setting of the mineralization. The sulfides, believed to be syngenetic, include chalcopyrite, cubanite, pentlandite, pyrrhotite, and minor amounts of bornite. They occur disseminated in the silicates and as small interstitial masses. The ratio of copper to nickel is about 3.8:1, based on 66 chemical analyses of rock samples from various outcrops (Ref. 1, p. 702, and Ref. 2). Test Procedures: With specified exceptions, all nickel and copper tests were made by the chromo-graph method,:' which measures the intensity of a colored spot formed by a reaction between the metal being determined and special reagent paper. The intensity is then compared to the intensity of spots prepared from samples of known metal content. Details of the test procedure are outlined in another article (Ref. 4, pp. 77,78). All soil samples tested in this investigation to date have been weighed on an analytical balance. However, a volumetric scoop designed to provide about 0.1 g of soil adds to the speed and ease of testing and has been found to give satisfactory results (Ref. 5, p. 531, and Ref. 19). The size of the samples used for the tests was 0.1 g. Whenever such small samples are used there is some question as to whether they are representative of the several grams in the field sample. Many repeat tests of the samples used in this investigation demonstrated that results can be reproduced within the limits of accuracy of the method without formal mixing beyond that inherent in screening the soil fractions. Furthermore, the 0.1 g is probably as representative of the field sample as the field sample is of its area of influence. Hawkes and Lakin (Ref. 6, p. 291), who considered the general problem, compared ground and quartered bulk samples of 500 g with 5-g grab samples. They concluded that "there is no significant loss in accuracy of data by substituting grab samples for bulk samples." The term soil implies somewhat different things to the geologist, engineer, and soil scientist.' For convenience the term as used in this article refers to unconsolidated material (the mantle) overlying bed rock. Sampling Procedure: Samples were taken at 100-ft intervals along north-south traverse lines across the gabbro-granite contact. The soil (till) samples were taken at an average depth of 1 ft, which was below the high-humus surface layer and into clean till. Samples taken at 1-ft intervals down to ledge showed as high a metal content at 1 ft below the air-surface as at greater depths and in two instances were slightly higher. The till at 1-ft depth did not appear to differ from material at greater depths. Total depth to bedrock has been tested at only a few points and where measured varied from 1 to 10 ft. Aerial Distribution Contours and Profiles: Plotting of copper, nickel, and cobalt content in contour form (Fig. 2) shows that anomalous amounts of these metal ions occur in till over and closely adjacent to mineralized areas of the gabbro. Contouring nickel content alone, or the copper content, outlines the same target area. Contours of the copper content provide a more distinct anomaly than nickel because of the higher copper concentration. The traverses are rather widely spread for interpolation; however, drilling has confirmed the target area essentially as shown.

Jan 1, 1959

By Paul J. Fopiano

SOME relationships between the flow characteristics of iron single crystals of 99.9 pct purity and the behavior of imperfections have been investigated. X-ray rocking-curve measurements and etch-pit counts were made as a function of plastic strain, and compared to the stress-strain curve obtained on a modified Polyani tensile machine. Crystals grown from rolled strips of vacuum-melted iron by the strain-anneal method1 had a high preference for a (110) longitudinal direction and a (211) face normal. The tensile specimens were prepared from 2 by i by 0.040 in. single crystals having a gage area of 3 by \ in. Rocking-curve measurements were carried out with a highly perfect germanium monochromating crystal in which the dazz spacing was matched to that of the dZl1 in the ir0n.l Well-collimated CuKal radiation was used throughout. These procedures practically eliminated errors due to geometrical and wavelength resolution. Inasmuch as the rocking-curve half breadth may vary markedly from point to point in the specimen being irradiated, the crystals were strained in place by mounting a hydraulic loading device on the double-crystal spectrometer. The rocking curves were taken after each increment of strain in the unloaded condition, since no observable difference was found in the rocking-curves between the loaded and unloaded states. The rocking-curve half breadths of the as-grown specimens were in the range 90 to 120 sec of arc when the beam irradiated an area of about -£ by -& in. on the specimen. DeMarco and weiss3 have shown that, for a well-colli- mated X-ray beam, irradiating about 10"! sq in. of the very same material, half breadths within 10 pct of the Darwin natural half breadth were observed. Since the rocking-curve specimens were stressed by the load-unload technique, the strain achieved at any given stress depended on the time of holding because of low-temperature creep. Fig. 1 shows the rocking-curve half breadth (also area/peak height) as a function of plastic strain for a relatively short holding time (2 to 5 min) at each stress level. For strains less than 0.1 pct the rocking-curve breadth is essentially constant; it is only for larger strains that there occurs a significant increase in this breadth. Where the holding times at each stress level were longer (by well over an order of magnitude) there occurs a significant increase in the rocking-curve breadth only after plastic strains of the order of 0.6 pct had been introduced into the specimen. This observation is related to the time dependence of creep phenomena and emphasizes the difficulty in comparing data obtained by two such different straining methods. Etch-pit results were obtained using a 2 pct nital etch on specimens strained in the range of 0 to 1 pct. Prior to etching, all specimens were annealed for 3 hr at 150°C, the carbon content being sufficient to decorate the dislocations for strains of at least 1 p~t.~ The data points were all taken from parts of the same single crystal which had been strained with short holding times at stress in increments of strain of the order of tenths of 1 pct. The (211) plane is particularly difficult to etch-pit in vacuum-melted iron and therefore it is felt that these values are as much as an order of magnitude low. Fig. 2 shows the etch-pit density as a function of plastic strain. The smooth curve passing through the data points is not meant to infer a quantitative correlation with the rocking-curve data. What is of interest, however, is the change in etch-pit density in the region of 0.2 pct plastic strain. The first three increments in strain (points 2,3, and 4) did not produce a measurable change in the etch-pit density while subsequent increments did produce a measurable change. While the absolute values of these results do not appear to be cor-

Jan 1, 1967

By W. A. Griffith

Fractional replication of factorial design, a general method for planning experimentation and for analysis of data obtained, is described as applied to a flotation investigation. This method leads to determination of optimum process conditions with minimum experimental effort. Its advantages over simple factorial design are demonstrated. A METHOD for planning experimentation and for analyzing the data secured will be demonstrated. This method, termed fractional replication of factorial design, employs a relatively small number of individual experiments to determine which of a large number of process variables are controlling, to determine which combination of levels of these variables is most likely to produce optimum results, and also to predict what results will be obtained with a particular combination of conditions not already tested. Although the general method is not new, having been developed by Finney in 1945,' the extent to which it can improve the effectiveness of applied research generally has not been recognized by metallurgists. The fractional replication procedure is particularly useful in flotation experimentation and an example from such an investigation will be used in the paper. However, it has equal value in any investigation in which similar experimental difficulties are encountered. In developing a flotation process for a particular mineral separation, the investigator is inevitably confronted with the following difficulties: 1—There are a large number of potentially important process variables. 2—Results of individual experiments are not highly reproducible, due in part to the difficulty in precisely controlling all the variables. 3—Considerable effort is expended in conducting individual experiments. 4—There are practical limits on the number of individual experiments which can be made. In situations of just such a type, modern statistical methods of planning experimentation and analyzing data have their greatest value. Applications of one such technique, called factorial design, to problems of this type have been described by Dorenfeld and others.'-' The simple factorial design is an efficient procedure when the investigator hopes to provide a comprehensive understanding of the effects and interrelationships of a small number of variables over a limited range. In applied research, this is seldom the main objective. Rather, the investigator usually wishes to determine which of the many potentially important variables are in fact controlling, which levels of the controlling variables will provide opti- mum metallurgical results, and what these results will be at optimum conditions. Interest in detailed trends is limited to the controlling variables and to levels in the region of optimum conditions. Simple factorial design has serious deficiencies for such objectives and is not the most efficient method of experimental design. Deficiencies of Factorial Design In a simple factorial design, an experiment must be made at every possible combination of each level of every variable, once these have been chosen and the levels of each to be included have been decided upon. As the number of variables or levels of each increases, the experimental program quickly reaches prohibitive size. For example, an investigation of only four variables, each at four levels, requires 256 individual experiments. Often upon completion of such an extensive program, it is found that one or more of the variables has metallurgically unimportant effects or that a poor estimate has been made as to the appropriate range of levels to be investigated. The result is that only a small proportion of the data obtained falls in the range of real metallurgical interest. Indeed, it frequently can be anticipated that certain combinations of levels of variables will not produce results of interest, but they still must be included so that the essential balance, or orthogonality, of the design will be retained. It may be true that factorial design will provide the greatest amount of information from a given number of experiments, but it does not necessarily follow that it will lead to the desired information with a minimum number of experiments. Much of the information provided may be of little value. Advantages of Fractional Replication The disadvantages of simple factorial design are overcome to a great extent by a modification known as fractional replication. This is a technique for sampling systematically the potential data of a full factorial experiment, that is, the data which would have been obtained if the complete factorially designed program had been completed. Only a fraction of the total array of experiments is made, but the experiments are chosen in such a way that the important advantages of factorial design and the accompanying analysis of variance are retained. The data obtained from the first group of experiments are used to determine which of several variables are controlling and which levels of these variables are most likely to produce the desired result. Unimportant variables and levels of variables then may be

Jan 1, 1956

By D. W. Fuerstenau, P. C. Kapur

The kinetics of green pelletization in a laboratory balling drum have been studied, using pulverized limestone as a model system. The growth characteristics of green pellets were found to be extremely sensitive to the moisture content of the material. Empirical kinetic equations, which incorporate a function of specific surface of the pellets as the criterion for growth potential, have been found to describe growth in a nucleation region and in a ball growth region. The rate constants in the kinetic equations are strongly dependent on the moisture content of the material being pelletized. Size distributions of the balls at different stages of pelletiz-ing are also discussed. In many industrial chemical processes, particulate matter can only be utilized if it is in an agglomerated form, such as pellets. Pelletizing is now widely used in iron ore technology1, and it has also been applied to a number of diverse fields such as the production of cement-kiln feed2, fertilizers3, and fluorspar4. Recently, it has been proposed to pelletize dispersion-type ceramic nuclear fuel elements5. In iron ore technology, for example, the production of agglomerates by pelletizing involves two major steps: 1) the preparation of green balls by rolling particles in a suitable balling device and 2) the firing of the green balls to form compact, strong bodies upon sintering. The critical step in a successful iron ore pelletizing operation is generally considered to be the balling operation1. In this paper, which is not concerned with the sintering of green pellets, the words green pelletizing and balling will be used interchangeably. Green pelletizing, or balling, will be defined as the process of forming larger bodies by rolling fine particles on a surface without the application of direct pressure. Two recent literature surveys6" indicate that in spite of the considerable amount of industrial pelletizing, very little is known about the fundamental principles of balling and its kinetics. The first reported research on the kinetics of pelletizing is the work of Newitt and Conway-Jones8. Using silica sands of different sizes in a batch laboratory balling drum, they found that the average green pellet diameter increased linearly with time at constant drum speed, and qualitatively the growth rate increased with moisture content. Generalized conclusions cannot be drawn from their research since the materials which they pelletized were sand and sand-silt mixtures rather than comminuted materials. Moreover, Newitt and Conway-Jones used testing sieves to estimate the size distribution of the green pellets, and this technique limited the range over which they could study the growth kinetics of the pellets. Bhrany and co-workers9 investigated the kinetics of balling iron ore fines on disk pelletizers ranging in diam from 1 to 18 ft. In their investigation, balling was carried out as a continuous operation, and growth kinetics were studied in terms of retention time of the material on the disk. Although the feed material in their study was quite coarse (the maximum size being about 1/2 in.), they also found qualitative relationships between pellet growth and water content, and feed size. In the present investigation, a number of innovations were introduced that refined the experimental measurements and established the reproducibility of balling experimentation. This enabled extension of the range of measurements to include study of agglomerate nucleation phenomena in the fractional mm size. This paper presents a detailed analysis of the nucleation and growth of green pellets in a laboratory balling drum. MATERIALS AND METHOD Pulverized limestone of specific gravity 2.72 was used as a model system in these studies. It has already been established that the balling characteristics of limestone and silica are similar to those of iron ore concentrates 2,8,10, depending on physical, rather than chemical, properties of the particles. The size distribution of the limestone was determined by a wet-dry sieving technique in the sieve range and by a sedimentation balance in the sub-sieve range. Fig. 1 presents the size distribution of the limestone used in this research. This figure shows that the material is finer than 200p (65 mesh) and that 25% of it is finer than 12. The specific surface area of this powder, as measured by BET gas adsorption methods,

Jan 1, 1964