Search Documents

By Ronald J. McGregor, Anton R. Hendriksz

INTRODUCTION The NIMCIX contactor was developed in the late sixties at the National Institute for Metallurgy in South Africa, for the main purpose of extracting uraniuy from relatively low-grade uranium solutions. ' The concept was sufficiently developed by 1974 to permit the operation of a 1500 litre per minute column to treat unclarified urgnium-bearing solution under industrial conditions. The plant operated so successfully that several large-scale units have subsequent18 teen constructed and are operating efficiently. ' Total resin consumption figures as210w as 3 per cent per annum have been reported, and the extraction efficiencies of the installations have been exceptionally good. All the solutions treated have been acidic in nature and in every case have been of such a nature as to significantly depress the maximum uranium loading on the resin beads. The chemical constituents of the leach solution fed to the loading column determines the ease with which the soluble uranium can be absorbed by an ion exchange resin. The concentration of competing ions and various resin foulants define the nature of the continuous ion exchange equipment required for efficient extraction performance. The basic NIMCIX concept is retained in these various alternate designs, but batch or continuous elution methods with fully automated or manual control systems may be selected for the different process solution conditions. BRIEF DESCRIPTION OF THE NIMCIX CONTACTOR The NIMCIX ion exchange column is a multi-stage fluidized bed contactor which is designed to permit a counter-current flow of solution and resin in order to maximize extraction efficiency. The flow o resin is periodic, based on the Cloete and Streat 5 principle of reverse flow. The basic operating cycle can be divided into five periods, as illustrated in Figure 1. For the purposes of describing the cycle, the loading column is considered; how- ever, a similar sequence may be applied to the operation of a NIMCIX continuous elution column. S'I'EP 1: During the forward flow period, Solution flows up through the column and fluidizes the resin contained in each stage. The uranium depleted solution overflows the upper flared-section of the column as the tail solution. Step 2: The flow of solution to the column is temporarily stopped and the resin is allowed to Settle onto the plates in the column. STEP 3: During the reverse flow period, a fixed volume of solution is drained from the column to the transfer vessel. This flow causes resin transfer to occur from stage to stage and also out of the bottom stage. At the end of this period, which expires when the low level switch (LS) is actuated, some of the resin is left in the conical section of the column and the pipe interconnecting the column and the vessel. STEP 4: The cone flush period ensures that all the resin left in the cone of the column and the pipe is transferred to the transfer vessel before forward flow commences. STEP 5: The loaded resin in the transfer vessel, which has been isolated from the column, can either be hydraulically transferred to the top of the elution column or eluted in the vessel before being transferred to the top of the loading column. The major advantages associated with the NIMCIX concept are: (i) Simple method of operation (ii) Contactor can readily be scaled up to large diameters. (iii) Because the resin does not pass through pumps or airlifts, and valves cannot be closed on resin beads at any time, the resin losses are exceptionally low. (iv) Operating costs are low. (v) Relatively few valves in the system. (vi) Stable operation is achieved over a wide range of feed rates and concentrations.

Jan 1, 1980

By Forrest F. Craig

Both early and more recent1 laboratory measurements of gas-oil relative permeabilities were made by subjecting oil-saturated cores to external gas drives. In these runs it was generally assumed that the flowing gas/oil ratio was constant at every point throughout the core. Using this assumption, the resulting equation for the gas-oil relative permeability ratio becomes A where Q is the flow rate at the producing end of the system. This relative permeability was then specified for the average saturation of the core. The assumption and resulting equation fail to recognize that in an external drive, the oil flow rate varies from zero at the injection end to a maximum at the producing end. The assumption also fails to consider the existence of a saturation gradient during an external drive. This saturation gradient is due solely to non-steady-state flow. Some authors2,3,4 have used this erroneous method of calculating relative permeability ratios from external drives. Other authors5,6,7 have used relative permeability ratios in calculating expected field performance by external drive, using Equation (1) and making the incorrect assumption that the saturation is uniform within the driven portion of the reservoir. To find the correct relation between the relative permeability. ratio at the average saturation and the produced gas/oil ratio, the insertion in Equation (1) of a factor, F, to adjust it for oil velocity gradient has been considered. The equation then takes the form Assuming that the oil velocity is directly proportional to the distance from the injection point leads to the conclusion that, for a linear system, the average oil velocity is one-half that at the producing end. Thus, the average gas/oil ratio is twice the produced gas/oil ratio, and Equation (2) is transformed into However, here again the existence of a saturation gradient is not recognized. The produced gas/oil ratio is related to the relative permeability ratio at the saturation of the outlet end of the system. It can be shown quite readily that, if the true log k,/k, saturation relationship is linear over small ranges of saturation, then d log kjko Where aaxv and u2 are the average and outlet saturations, respectively. Examination of Equation (3) shows that the factor F cannot have a single value over the whole saturation range. Another method of dealing with external drives is that of Buckley and Leverett.' Their method, which assumes no gas expansion and outlet end effects, considers fully the saturation gradient present due solely to the displacement mechanism during an external drive, and applies Equation (1) to each point along this gradient. The calculation of production history using relative permeability ratio is straightforward. However, the calculation of relative permeability ratios from production history by this method is impractical. In a recent paper, Welge,1 using an approach similar to that of Buckley and Leverett, presents a simplified method for determining the saturation at the producing end of the system. At thi; saturation, the relative permeability ratio calculated by Equation (1) applies exactly, after breakthrough of the stabilizej front. By this method, relative permeability ratios or production history may be calculated with equal ease. Use of this method as a laboratory tool in determination of relative permeability ratios by external drive requires the knowledge of the total injected fluid volume; in addition to the average saturation and producing ratio. As a preliminary test of this method of calculation, laborntory external gas drives were run on a 3% in. diameter, 11 in. long sample of Nellie Bly consolidated sandstone. This core sample had a permeability of 824 md and a porosity of 28.1 per cent. The driving fluid was air. Three different oils were used, having viscosities of 1.4, 9.8, and 125 cp, respectively. The runs were made under identical pressure differentials,

Jan 1, 1952

By G. S. Coe, G. F. Bolling

The columnar-to-equiaxed transition (CET) has often been studied as a function of concentration, C, temperature gradient in the melt, G, and rate of solidification, R. Although it may not be immediately obvious it is sure that the CET, expressed as some ratio G/R necessary to suppress equiaxed nucleation for a given C,' cannot be a single-valued function of C. This restriction which is ignored in theoretical presentations is most trivially illustrated by the following simple reasoning. An ordinary casting is an unconstrained system where the thermal conditions change continuously with progressive solidification, and thus where the ratio G/R decreases from the outset. A CET occurring early in a casting (much equiaxed material) must then mean a higher value of G/R at transition, in comparison with the late CET which implies a more columnar ingot. Therefore, when we choose to experiment under constant superheat for all alloys across the simplest binary, and we find an inversion in the amount of equiaxed material, we shall find two values of C satisfying the G/R condition. However, in addition to this simple instance, the early experiments of ~orthcott' seem to show that the amount of equiaxed material in cast ingots may be a multivalued function of composition in a given binary system. In the present investigation, the CET has been examined in Pb-Sb, Pb-Ag, Pb-Sn, and Sn-Pb under several types of casting conditions. First, we have used the "bidirectional" ingot section Simulation;~ and we have used .'i~z situ" ingot solidification where a hot mold already filled with liquid is chilled so that there is no disturbance provoked by pouring. Differences in method are trivial to our purpose except for the fact that certain variations are more simply observed by one technique than by the other. Second, we have varied the mode of superheating from the simple extreme of casting all alloys from the same temperature regardless of composition to the careful casting at a standard superheat, calculated with respect to the alloy liquidus. Third, we have made provisions for decanting and interface observation, so that some "castings" were not completed. The results obtained with Pb-Sb are exemplary; some of these are shown in Fig. 1 along with explanatory captions. At 9 to 10 at. pct Sb, and thus above the solid-solubility limit Cs = 5.8 at. pct Sb and yet below the eutectic concentration Ceu = 17.5 at. pct Sb, there is a persistent anomaly in all the variations. Measurement was easy enough and standard enough so that we may forgo discussion of either method or possible error, in view of the forceful self-consistency shown by experiment. Solute distributions were measured along the bidirectional ingots in order to determine if a solute buildup had occurred at the CET. No anomalous changes in concentration were observed along the bars at either side of the solute anomaly of 10 at. pct Sb. The microstructures of the ingots were also examined across the range of composition investigated. At concentrations less than 9 at. pct Sb the micro- and macrostructures were typical of the cellular-dendritic growth usually reported for the columnar zone. Above the anomaly, both the columnar grain size and the amount of primary phase decreased. This microstruc-tural observation can be correlated with observations made of decanted interfaces. Below 9 at. pct Sb, the interface reveals the "normal" columnar dendritic structure; however above the anomalous concentration range, preferential growth of certain dendrites appears to form a "super dendrite" pattern. In Fig. 2 a microstructural view of the etched '(super dendrites'' shows their constitution both in terms of relative phase proportion and in terms of the ordinary smaller dendrites.

Jan 1, 1969

By D. Turnbull, R. E. Cech

THE supercooling behavior of pure liquid metal droplets has been described.' The solidification behavior of small droplets of Cu-Ni alloys as a function of composition is described herein. The Cu-Ni system is of interest for such an investigation because the components are miscible in all proportions in both the liquid and solid state and the phase diagram of the system, Fig. 1, is a type often encountered in metal systems. It was noted by Van Riemsdyk2 and others that the nature of the "blick" of gold droplets during assaying is not perceptibly affected by the content of the platinum group metals that are soluble in both liquid and solid gold. With the exception of these very qualitative observations there appear to have been no prior investigations of the supercooling of alloy droplets. A series of nine alloys varying in composition by 10 pct (nominal weight) steps was made using certified OFHC copper, 99.996 pct purity, obtained from the American Brass Co., and electrolytic nickel pellets, 99.92 pct purity, obtained from the International Nickel Co. The alloys were melted and cast in vacuo in a high frequency induction furnace. Segregation was minimized by chill casting into 1-in. diam steel molds. Each alloy rod was homogenized at a temperature 50°C lower than the solidus for 100 hr in dry hydrogen. Samples for solidification experiments and chemical analysis were machined from the same portion of each rod and at a depth of in. from the rod surface. A Car-boloy tool was used in taking turnings to minimize contamination of sample by the cutting tool. As a check, a number of samples for microscopic observation were chipped from the ingot with a fragment of pyrex glass. The resultant data in all cases agreed, within experimental error, with the data obtained from machined turnings. A number of small particles, 30 to 50 micron diam, of a Cu-Ni alloy were placed in a quartz vial 1/32 in. diam x % in. long. This was inserted into a platinum ribbon heater and mounted in a high temperature microscope stage. The general procedure followed in inserting samples into the stage and measuring the temperature of melting and solidification was the same as described in an earlier paper.' A chromel-alumel thermocouple was used to study alloys of 88.93 pct Cu, 80.62 pct Cu, and 72.27 pct Cu. A platinum-10 pct rhodium thermocouple was used for the remainder of the alloys. Procedure The assembled high temperature microscope stage was first evacuated and the sample given a 2 min in situ treatment at 900 °C in dry hydrogen to remove surface oxides which might otherwise catalyze crystal nucleation. After this treatment the atmosphere was changed to purified helium and melting and solidification were observed microscopically. It was not possible to determine the solidus with certainty. However, the liquidus temperature was measured very accurately since the surface of the particle changed in appearance from coarse and uneven to smooth, mirror-like within a few degrees. The visually determined liquidus temperature was checked by finding the highest temperature to which the particle could be heated such that it did not supercool upon lowering the temperature. This temperature which checked within experimental error with that found by observing surface change was assumed to be characteristic of the thermody-namic liquidus. In order to minimize the effect of selective evaporation of copper from the alloy, samples were changed frequently and the droplets were heated no more than 50 °C in excess of the liquidus temperature T,. When the latter conditions were fulfilled, T, values were reproducible. After melting, the temperature was decreased slowly until the particles solidified. The sudden release of heat of fusion on solidification caused the

Jan 1, 1952

By C. F. Baes, H. H. Kellogg

The effect of dissolved sulphur on the surface tension of liquid copper has been measured by a sessile-drop method at 1120°C. Sulphur is highly "surface active" in liquid copper; the surface tension lowering is about 650 dynes per cm for a sulphur content of 0.6 pct. The degree of adsorption of sulphur on the liquid copper surface has been calculated, and various models of the surface layer have been tested against the results. THE science and technology of surface chemistry and physics at room temperature are relatively advanced, and the engineer possesses a background of information that allows him to predict and control many phenomena—such as frothing, emulsifica-tion, wetting of solids by liquids, flow through capillaries, adsorption, mineral flotation, penetration of porous bodies by liquids—which are related to the behavior at interfaces. Surface phenomena in high temperature systems are probably of equal importance to the metallurgist and ceramist, but at present are incompletely understood. Frothing of liquid slags, penetration of refractories by liquid metals and slags, infiltration of powder metallurgy "skeletons" by liquid metals, bonding of metals to ceramic oxides in the manufacture of "cermets," wetting of one metal by another for purposes of soldering and welding, are but a few of the many phenomena of metallurgy that are controlled by surface behavior. Concerning the most elementary of surface properties—the gas-liquid surface tension—the existing data on liquid metals are inadequate. For the surface tension of "pure" metals the results of different investigators differ by 10 or 20 pct. The role of surface-active agents in liquid metals is relatively unknown, and it is not unlikely that the discrepancy in the values for the surface tension of pure metals is due to the presence of traces of surface-tension-lowering impurities in the samples. Based on analogy to room-temperature systems, the highly surface-active impurities in liquid metals should be elements or compounds of limited solubility in the liquid metal and possessing considerably weaker in-termolecular bonding forces than the metal itself. The nonmetals of the V, VI, and VII groups of the periodic table might be expected to act in this way. The experimental work in this paper represents a study of the surface tension of pure liquid copper in hydrogen at 1120°C, and the effect of dissolved sulphur (which is a strongly surface-active agent in liquid copper) on the surface tension of copper. The Gibbs equation has been applied to the results and an attempt made to correlate the surface tension behavior with two models of the interfacial layer. The sessile-drop method for surface-tension measurement at high temperatures which was used in this investigation has been described by several investigators.1,2 It involves the measurements of the dimensions of a stationary liquid drop on a horizontal surface. From these dimensions and the liquid density, the surface tension of the liquid may be calculated either by means of various approximate formulas, or more accurately by means of the Bash-forth and Adams tables.- or liquid drops which exhibit a contact angle greater than 90°, by far the simplest dimensions to measure are those indicated in Fig. 1, i.e., the diameter of the maximum horizontal cross section (2X) and the maximum height (Z) above this plane. In order to calculate the contact angle and the volume of the drop, the diameter at the base of the drop (2a) and the total height (h) also must be measured. Upon investigation of the literature, the authors failed to discover a concise statement of the accuracy of the method when used in conjunction with the Bashforth and Adams tables, hence a careful examination of the tables was made and the curve shown in Fig. 2 was constructed. The "error factor," plotted as the ordinate, is the factor by which the

Jan 1, 1954

By Heribert Wiedemeier, Ali Khan

The tetnperature-conposition phase diagram of the syste,n MnS-CdS has been investigated over a range fro)n 600° to 1000° C. The small degree of solid solution on the MnS rich slde cornpared to the large CdS rich single-plzase region could be esplained by the difference in strzrctures of the cotnponents and the relative sizes of the metal ions. The change in lattice parameter as a function of cot)zposition is linear atzd jollous Vegavd's rule within experittzental limits of error. Mns and CdS have gained considerable attention in recent years because of their interesting solid state properties. In order to investigate such properties and thermodynamic quantities of mixed compositions of transition metal chalcogenides, it is necessary to establish their phase diagrams. Two earlier publications are concerned with the MnS-CdS system. Schnaase' found complete miscibility when MnS and CdS were co-precipitated in aqueous solutions at room temperature. However, there is evidence that the mixed crystals were not stable over the entire range of composition. roger's' investigation of the ternary system ZnS-CdS-MnS at 900" C indicates a solution of 48 mole pct MnS in CdS. During the course of this work which covers a temperature range from 600" to 100O° C, emphasis has been placed on the use of high-purity materials and well-defined modifications of starting products. EXPERIMENTAL Starting Materials. The stable green modification of MnS with NaCl structure was prepared by direct fusion of the elements and by chemical transport reaction. Stoichiometric amounts of manganese of 99.99+ pct purity and of sulfur of 99.999 pct purity from Gallard-Schlesinger Corp. were sealed in an evacuated quartz tube at a pressure of 10"5 mm Hg and annealed for 48 to 60 hr at 800°C. The quartz tubes used throughout this work were outgassed under high vacuum at a temperature of 90OU to 1000° C for 12 to 18 hr. For the synthesis of MnS by chemical transport3 stoichiometric ratios of elemental manganese and sulfur and small amounts of iodine were sealed under vacuum in a quartz tube, which was placed in the temperature gradient, 900/800-C, of a two zone furnace. Single crystals of MnS in form of platelets and octahedra were deposited at the cooler end of the tube. The manganese content was determined by EDTA titratin. A typical result of several analyses was 63.05 pct Mn which agreed well with the computed value, 63.14 pct. Commercially available CdS of 99.9+ pct purity, a mixture of the cubic and hexagonal forms, was converted into the hexagonal high-temperature modification by annealing CdS under vacuum in quartz tubes at 700° C for about 36 hr. As a second method, cadmium sulfide was transported in a temperature gradient, 70O° /650" C, with in vacuo sublimed iodine from the hot to the cooler zone of the tube. Large sized crystals in form of hexagonal columns and some platelets were obtained. An average result of the cadmium analysis by EDTA titration4 was 77.93 pct in good agreement with the computed value, 77.81 pct. Only the NaCl modification of MnS and the wurtzite modification of CdS were used for the solid solutions. Preparation of Solid Solutions. Finely powdered MnS and CdS samples were thoroughly mixed in molar ratios and compressed into pellets of about 6 mm diam and 2 mm in height at a pressure of approximately 13,500 atm.

Jan 1, 1969

By G. C. Wallick, J. S. Aronofsky

This paper presents an experimental verification of numerical solutions of the differential equation describing the transient flow of an ideal gas through a porous material including gas slip effects. Several experiments previously proposed by one of the authors are described together with the correlation of the experimental results with the numerical solutions. Excellent agreement is shown between experiment and theory, thus demonstrating the validity of both the basic assumption made in the derivation of the differential equation and of the numerical solutions of the differential equation. INTRODUCTION In a recent paper,' one of the authors has discussed a numerical solution of the differential equation describing the unsteady flow of an ideal gas through a porous material. This paper was primarily concerned with the effect of gas slip on the flow behavior. In steady state measurements, Klinkenberg2 and others have observed that the apparent gas permeability Kn of a porous material may be expressed as a linear function of the reciprocal average pressure 1/P of the form, ka = b + m/P........(1) where the intercept b is a constant for a given porous material and the slope m is a constant for a given porous material-gas combination. Aronofsky1 assumed that Equation (1) could be used to describe transient flow if P was replaced by the average pressure in an infinitesimal volume element. A number of numerical solutions of the resultant differential equation were presented and specific flow experiments were proposed as a means of testing the numerical solutions. The purpose of the present paper is to describe the results of several of these experiments and the correlation of the experimental data with the results of the numerical analysis. Reference should be made to the previous paper' for details of the mathematical and numerical analysis. EXPERIMENTAL PROCEDURE A simplified schematic drawing of the experimental apparatus is shown in Fig. 1. An 8.89 cm diameter limestone core of length L = 14.82 cm was mounted in a conventional rubber sleeve permeameter cell. The external sleeve pressure in this cell was maintained equal to 500 psig during all experimental tests, thus permitting only longitudinal flow through the core. One end of the core, X = L, was connected to an adjustable cavity of volume V as well as a pressure gauge or manometer. The other end of the core, X = 0, was connected to a Nullmatic pressure regulator which was in turn connected to a nitrogen cylinder or vacuum pump. In each experiment, the core and adjacent dead volume were charged with nitrogen to a uniform gas pressure P1. The pressure at X = 0 was then suddenly raised or lowered to a constant pressure P, and the resulting pressure buildup or pressure decline at X — L was recorded as a function of the time. The porosity + and the permeability constants m and b were determined using the transient method which was suggested by the results of the numerical analysis and which was described in the earlier paper.' The following average values for these constants were obtained as a result of 19 separate evaluations: 171 -0.259 2 0.005 md-atm; b = 0.374 5 0.008 md; and F = 0.098 ± 0.005*. It should be noted that in order that transient phenomena will be slow enough to be observed conveniently, it is desirable to use a porous material of very low permeability, e.g., less than I md. Moreover, it is in this low permeability range that the gas slip effects are most significant. Thus, such a low permeability sample is very suitable for use in an experimental verification of the numerical solutions previously described.

Jan 1, 1955

By U. N. Bhrany, N. Arbiter

F. C. Bond (Consulting Engineer, Processing Machin-Dept., AllisChalmers Manufacturing Co., Milwaukee) This is a very comprehensive paper. It deals with 1) size distribution functions, 2) energy-particle size relationships, and 3) criteria for grinding efficiency. The method of analysis used consists of tracing the increase and decrease of successive Tyler scale screen fractions as grinding progresses. The treatment is based upon the assumption that the work done in grinding is concentrated entirely upon the coarse size fractions in the mill. It is postulated that when a feed consisting entirely of particles of one size fraction is batch ground that finer fractions increase at a constant rate as grinding progresses, until finally each in turn begins to be ground and to decrease in amount. It is shown that under these conditions the empirical Gates-Gaudin-Schuhmann log-log size distribution plot can be logically explained. It is difficult to imagine any condition of ball mill operation, or any probability analysis, which would result in no grinding contacts on the finer particles present. However, the authors' postulate must be accepted if their conclusions are to be sustained. They conclude that the equation expressing the for mation rate of each size fraction can be resolved into the G.-G.-S. size distribution equation, and that the curvilinear portion of the log-logsize distribution line represents the only particle sizes which are being ground as well as produced. The author's statement that the G.-G.-S. Eq. 4 applies only to the finer portion of the size distribution while the Rosin-Rammler Eq. 6 applies over the whole size distribution range certainly supports the view that an exponential equation of the Rosin-Rammler type corresponds more closely to actual size distribution than does the power law G.-G.-S. equation. However, the G.-G.-S. equation is used to derive Eq. 8 and Eq. 9, which state that the energy required to grind a ton through a particular mesh varies inversely as the particle size to the exponent -a , whic h is the slope of the extended G.-G.-S straight line. This is stated to apply to both batch and continuous grinding. The G.-G.-S. slope exponent thus becomes very important in all crushing and grinding analyses by this method. The authors include some data which can be used to test their conclusions, as described below: We will consider grinding a feed composed of 14 mesh (1190u) to -200 mesh (74u). The average of 17 Allis-Chalmers grindability tests on quartz gives a work index (Wi) of 12.77. The -200 mesh product should have 80 pct passing 41.6u or P = 41.6. The equivalent 80 pct passing size of the 14 mesh feed particles is considered to be the particle diameter, as is true for spherical particles. Thus F is 1190u. In the work index equation 10 wi 10 wi 127.7 127.7 P f 41.6 1190 = 16.1 kw-hrperton. In Table I the ball mill exponent for quartz is 1.0 and from Eqs. 9 and 10: w=kw-hr/ton =Constan/(do/di) = 16.1 ton so that (1190) Constant/ Constant and the Constant equals 259. The tons per horsepower-hour value is 0.0464, in agreement with Table IV. The exponents for other minerals are listed in Tabl e 11. In Table IIA, the kilowatt hours per ton required to grind these minerals between the same sizes as the quartz, both by the authors' method and the work index equation, are listed. According to the authors kw-hrperton = 259/(16.1)".

Jan 1, 1961

By A. D. Mutch

IT has long been known that there are wide variations in the thermo-electric potential of pyrite. In the last few years a study of this variation and other properties of pyrite was carried out at the University of Toronto by F. G. Smith on specimens of pyrite from many loyalities. As a result of these studies the pyrite geothermometer was developed by Professor Smith. The calibration of the pyrite geothermometer is based on the assumption that "crystals of any one electronically conducting mineral species deposited at a high temperature have a more positive thermoelectric potential than crystals deposited at a low temperature, and the thermo-electric potential varies continuously between any given limits of the temperature of deposition." The cause of the variation was postulated to be crystal discontinuities, such as lineage boundaries. As yet this postulated relationship has not been established as fact. On theoretical grounds it is indicated that the thermo-electric potential of pyrite or any other partial conductor is dependent upon a number of factors as well as the temperature of formation. For this and other reasons there is ground for doubt as to the validity of the relationship postulated above. To date a large number of readings have been made with the pyrite geothermometer which have been recorded in terms of the calibration of the instrument in degrees centigrade. So far most of the pyrite tested by this technique has been in association with gold ore. The values for the temperature of formation obtained in these tests appear to be reasonable and, in relation to geological problems, to vary in the right sense. Determinations on some specimens not in association with gold ore have given readings that appear to be much too high, taking into account the probable maximum temperature of the host rock. To sum up, although many determinations made with the pyrite geothermometer appear to give reasonable values for the temperature of formation of the pyrite, there are at the same time other readings which seem anomalous. On the basis of an examination of available results of work with the pyrite geothermometer, two problems must be considered: 1—the validity of the calibration of the pyrite geothermometer; and 2— the nature of the distribution of the thermo-electric potentials of pyrite in relation to the specific problems that have been examined. In this paper the academic problem relating to the validity of the calibration of the pyrite geothermometer will not be discussed. Instead, the variations of thermo-electric potential of pyrite in association with gold veins will be discussed with regard to specific examinations. The first of these two examinations deals with the distribution of thermo-electric potentials of pyrite in association with a gold vein and was carried out in the McIntyre Porcupine Mine, Schumacher, Ont., in the summer of 1949 as part of a study with this technique. The second problem deals with variations of thermo-electric potential of pyrite in gold ore specimens from all over the world. The specimens used in this study were part of the private collection of Frank Ebbutt, geologist for Howe Sound Exploration Co. Results of Detailed and Reconnaissance Studies In discussion of the results of determinations made in these two studies the values of thermo-electric potential in millivolts will be used. It will be noted that in some of the diagrams the determinations are recorded in terms of pyrite geothermometer temperatures. Fig. 1 is a calibration chart relating pyrite geothermometer temperatures and thermoelectric potentials. A full discussion of the technique is given in the original paper by F. G. Smith.' Briefly, the determinations are made by placing two stainless steel probes, one hot and the other cold, against the surface of the pyrite crystal being tested. The thermoelectric potential that is developed is read with a null-point potentiometer that is calibrated in what

Jan 1, 1953

By C. van der Poel

When oil is displaced from a horizontal formation by another fluid of lower density, the latter tends to override the former in the shape of a tongue owing to gravity segregation. This gravity tongue has an adverse influence on the oil recovery. If the fluids are miscible, diffusion (mixing) takes place at the interfacial boundary of the gravity tongue. This mixing should have a favorable effect on oil recovery. The report describes a laboratory study of the magnitude of the mixing zone under various conditions, so as to assess the effect of diffusion on oil recovery both in laboratory experiments and under actual field conditions. The technique used enables visual observation and measurement of the size of the mixing zone in transparent glass-powder packs. The results show that in experiments in small models and cores the width of the mixing zone may well be of the same order of magnitude as the height of the model. In such cases oil recovery is favorably affected by mixing. It can further be concluded that, under conditions prevailing in the field, mixing of the injected fluid with reservoir oil is equal to that caused by molecular diffusion alone, eddy-mixing not taking place to any appreciable extent. A simple calculation, then, shows that molecular diffusion is too small for a beneficial effect to be expected from the injection of miscible fluids in horizontal or nearly horizontal reservoirs unless pay zones are thin. INTRODUCTION This paper gives results of experiments made on the mixing which occurs when a miscible displacement is carried out in a horizontal reservoir. This mixing takes place at the interfacial boundary of the gravity tongue formed when the lighter injected fluid overrides the oil present in the reservoir. The object of the experiments was to simulate the field case where, for example, propane with a viscosity of 0.075 cp under reservoir conditions displaces an oil of viscosity 0.6 cp. In the experiments the lighter fluid had a viscosity of about 1 cp and the heavier one a viscosity of about 8 cp, so as to obtain the same viscosity ratio. In order to enable the results to be compared with those published in the literature, a set of experiments with viscosity ratio equal to one was also performed. EXPERIMENTAL TECHNIQUE The technique developed for the purpose enables the width of the mixing zone to be studied as a function of time and place. A glass-powder pack is saturated with a water-glycerine mixture of suitable viscosity, which represents the reservoir fluid. The pack is then rendered transparent by dissolving sufficient ammonium thiocyanate in the mixture to obtain a solution with a refractive index matching that of the glass powder. Alkaline water to which phenol-phthalein has been added is used as the lighter and less-viscous displacing fluid. In those places where mixing or diffusion of the two liquids occurs, the slightly acidic ammonium-thiocyanate solution neutralizes the alkali in the water, and the glass pack shows up white against the red-colored invading water. If a black cloth is hung over the back of the apparatus, the transparent part of the glass pack appears black. In this way the position of the two phases and the transition zone between them is clearly visible, as shown in Fig. 1 (where the red-colored invading water is the grey-shaded zone lying uppermost). In all experiments the amount of alkali in the water was chosen such that the upper contour corresponded to a concentration of 5 per cent of the dense liquid. The lower contour is determined by the size of the glass grains and the thickness of the pack and, consequently, varies from experiment to experiment. Corresponding concentrations of the dense liquid range from about 97 per cent for packs of small grain size (permeability = 1 darcy) to about 95 per cent for those with large grains

By A. D. Mutch

IT has long been known that there are wide variations in the thermo-electric potential of pyrite. In the last few years a study of this variation and other properties of pyrite was carried out at the University of Toronto by F. G. Smith on specimens of pyrite from many localities.1-3 As a result of these studies the pyrite geothermometer was developed by Professor Smith. The calibration of the pyrite geothermometer is based on the assumption that "crystals of any one electronically conducting mineral species deposited at a high temperature have a more positive thermo-electric potential than crystals deposited at a low temperature, and the thermo-electric potential varies continuously between any given limits of the temperature of deposition." The cause of the variation was postulated to be crystal discontinuities, such as lineage boundaries. As yet this postulated relationship has not been established as fact. On theoretical grounds it is indicated that the thermo-electric potential of pyrite or any other partial conductor is dependent upon a number of factors as well as the temperature of formation. For this and other reasons there is ground for doubt as to the validity of the relationship postulated above. To date a large number of readings have been made with the pyrite geothermometer which have been recorded in terms of the calibration of the instrument in degrees centigrade. So far most of the pyrite tested by this technique has been in association with gold ore. The values for the temperature of formation obtained in these tests appear to be reasonable and, in relation to geological problems, to vary in the right sense. Determinations on some specimens not in association with gold ore have given readings that appear to be much too high, taking into account the probable maximum temperature of the host rock. To sum up, although many determinations made with the. pyrite geothermometer appear to give reasonable values for the temperature of formation of the pyrite, there are at the same time other readings which seem anomalous. On the basis of an examination of available results of work with the pyrite geothermometer, two problems must be considered: 1-the validity of the calibration of the pyrite geothermometer; and 2-the nature of the distribution of the thermo-electric potentials of pyrite in relation to the specific problems that have been examined. In this paper the academic problem relating to the validity of the calibration of the pyrite geothermometer will not be discussed. Instead, the variations of thermo-electric potential of pyrite in association with gold veins will be discussed with regard to specific examinations. The first of these two examinations deals with the distribution of thermo-electric potentials of pyrite in association with a gold vein and was carried out in the McIntyre Porcupine Mine, Schumacher, Ont., in the summer of 1949 as part of a study with this technique. The second problem deals with variations of thermo-electric potential of pyrite in gold ore specimens from all over the world. The specimens used in this study were part of the private collection of Frank Ebbutt, geologist for Howe Sound Exploration Co. Results of Detailed and Reconnaissance Studies In discussion of the results of determinations made in these two studies the values of thermo-electric potential in millivolts will be used. It will be noted that in some of the diagrams the determinations are recorded in terms of pyrite geothermometer temperatures. Fig. 1 is a calibration chart relating, pyrite geothermometer temperatures and thermo-electric potentials. A full discussion of the technique is given in the original paper by F. G. Smith.' Briefly, the determinations are made by placing two stainless steel probes, one hot and the other cold, against the surface of the pyrite crystal being tested. The thermo-electric potential that is developed is read with a null-point potentiometer that is calibrated in what

Jan 1, 1952

By C. H. Li, J. A. Sarteli, W. Bonfield

The stress to cause a permanent micros train of 2 x 10-6 in. per in. (defined as the microscopic yield stress) in beryllium is found to be very sensitive to surface condition. The initiation of plastic flow in as-machined specimens, which contain a high density of twins and large residual compressive stresses to a depth of about 0.010 in. from the surface, occurs by the nucleation of slip from high stress fields around twin tips in the surface layers. Removal of the damaged surface layers by chemical polishing results in an appreciable increase in the microscopic yield stress, which is attributed to the remoual of stress raising twins rather than to the release of residual stresses. IT is well known that the macroscopic mechanical properties of beryllium are very sensitive to the nature of its surface finish. Matthews and coworkers1 have shown that the tensile strength of rolled beryllium sheet can be considerably enhanced by chemical polishing, presumably due to the removal of the high concentration of notches, twins, and large residual compressive stresses normally found in the surface layers of ground specimens. There was no similar effect of surface condition on the macroscopic yield stress (the stress to produce 0.2 pct strain), a point to which we shall return. To date there have been no extensive studies on the effect of surface condition on the microstrain region of beryllium. Recent measurements2 of the stress to cause a permanent microstrain of the order of 10-6 in. per in. on machined hot pressed QMV beryllium, gave values ranging from about 600 to 3500 psi. Reasons for the considerable scatter in the results were not discussed. In this paper, we wish to show that the microstrain region of polycrystalline beryllium is, in fact, very sensitive to surface condition. The stress to cause a permanent strain of 2 x 10-6 in. per in. which for convenience we will refer to as the microscopic yield stress, can be increased by almost a factor of 2 by minimizing the amount of surface damage. The variation of microscopic yield stress is attributed to the presence of twins in the surface layers, rather than to residual compressive stresses. EXPERIMENTAL PROCEDURE The starting material for this investigation was commercial purity hot pressed QMV beryllium, with an average grain size of about 20 . The effects of surface treatment were assessed by comparing the microstructure with residual stress measurements and microstrain properties. Specimens were prepared for metallographic examination by wet grinding with 240, 400, and 600 papers, followed by polishing with a suspension of fine alumina in warm 20 pct oxalic acid. After this treatment, the surface could be examined, either in the unetched condition with polarized light, or after etching in 2 pct hydrofluoric or boiling 20 pct oxalic acid with normal light. Thin films of beryllium suitable for electron transmission study were prepared from bulk specimens by chemical polishing with hydrofluoric and phosphoric acids, and were examined in an RCA-EMU-3 electron microscope. Residual stress measurements were made by a procedure similar to that of Treuting and Read.3 Specimens 4.0 in. by 0.625 in. by 0.080 in. were prepared by three different machining procedures, namely: turning, grinding, and sawing, in order to introduce different levels of residual stresses. An electrical resistance strain gage was attached to one side of the specimen, and that surface was then protected by an acid-resistant enamel, while the opposite surface was removed by etching approximately one thousandth of an inch at a time. The change in specimen curvature in the longitudinal direction was noted and the residual stress was calculated by the following formula:

Jan 1, 1963

By Orville R. Lyons

CENTRIFUGAL force has been utilized for the dewatering of fine coal for over 50 years by means of machines commonly called centrifugal dryers. In any centrifuge the coal and water are subjected to a spinning action which usually inctedtocreases in intensity as the coal progresses through the machine. This spinning action, or the centrifugal force that it induces, tears water away from coal particles and produces a dewatered coal. The effectiveness of the dewatering action for any particular machine is governed by the size-consist of the coal and the centrifugal force imparted to the water on the coal. Four makes of centrifuge are currently being used to dewater fine bituminous coal in the United States. These include the vertical, basket-screen C.M.I. (Centrifugal and Mechanical Industries) Carpenter, Reineveld, and the horizontal solid-bowl Bird. The C.M.I. and Reineveld centrifuges are similar in construction and are variations of the Elmore centrifuge, varying principally in the slope of the basket, operating speed, size of perforations in the basket, and diameter of the basket. The same operating description suffices for both. These centrifuges are made up of two rotating elements, an outside conical screen frame and an inside solid cone which carries spiral pusher blades. The screen frame is supported by a hollow shaft, and the solid cone by a spindle shaft passing through the hollow shaft. Both the screen frame and the solid cone rotate at high speeds; the solid cone, with its scraping flights, rotating somewhat more slowly than the screen. The wet coal enters the machine at the top, falls on the cone, and is thrown against the screen by centrifugal force. The coal slides down the screen until it meets the upper end of the flights, where it is moved slowly along the flights until discharged at the bottom. The elfluent passes through the sieve basket or screen and discharges into a launder. The Carpenter centrifuge consists of a conical rotating element with a vertical axis, built up of three rows of perforated stainless-steel plates or stainless-steel wedge-wire or round-wire sections cut and rolled to conform to the surface of the cone. Wet material is delivered to a conical hopper inside the centrifuge casing and is distributed in a thin, even layer over the inner surface of the cone by a disk mounted on the rotor shaft. Impacting of the coal against the screen and the centrifugal force exerted as the material moves downwards creates a dewatering effect. The moisture is forced through the coal bed and then through the perforated plates of the basket, collecting in a circular trough built into the centrifuge. The effluent is conducted from the trough through the centrifuge casing by outlet pipes. The dewatered coal is discharged over the inside bottom edge of the rotor. The Bird centrifuge consists of a tank or truncated conical shell, which is revolved at the desired speed by means of a drive sheave. A screw conveyor rotates inside the cone or bowl at a slightly lower speed in the same direction of rotation. The feed entrance, in the center of the large end of the truncated cone, is high enough to allow formation of a pool of slurry. Adjustable effluent-discharge ports are located in the large end of the bowl so the level of the liquid can be regulated. The solids are moved forward by the screw conveyor as fast as deposited and carried above the level of the pool for an interval before leaving the bowl. In the past these four makes of centrifuge could be compared only by indirect means. The author published a rough comparison method in 1949' but

Jan 1, 1952

By J. E. Hilliard

AS is well-known, the length per unit volume, Lv, of any lineal feature (such as grain edges in a polycrys-talline specimen) can be estimated from a count of the number of point intersections with a random plane of polish. Specifically:' Lv = 2NA in which NA is the expected number of intersections per unit area of the plane. As with all other measurements in quantitative metallography, the estimation of Lv is subject to a sampling error, the magnitude of which will depend on the number of observations. The analysis is usually performed by counting the number of intersections falling within a closed test figure of known area applied to a succession of different areas on the plane of polish. The test figure can be drawn directly on the screen of a projection microscope or else superimposed on the image of the structure by means of a reticule in a focusing eyepiece of an ordinary microscope. The use of micrographs is very inefficient unless, as in the present study, these are also needed for some other purpose. Consider an analysis in which n applications are made of a test figure within which the expected number of intersections is N. The expected total number of intersections counted is then given by NT = nN If (N) is the standard deviation in N for one application of the figure and if the latter is applied to non-overlapping regions so that the N's are independent in the statistical sense, then the standard deviation in Nt is given by: a(NT) = a(N)vn Eliminating n from Eqs. [2] and [3] we obtain: Thus, k will be a constant for a given specimen and test figure. From Eqs. [ I] and [4] we find As is to be expected, therefore, the fractional error in the estimation of Lv is inversely proportional to the square root of the total number of intersections counted. It is very desirable to have at least a rough estimate of k before starting an analysis, since it will then be possible to decide in advance by use of Eq. [5J how many counts should be made in order to achieve the desired accuracy in Lv. Without such information the investigator, unless lucky, will either do needless (and tedious) work or else obtain a result with too large an error. It is very easy to calculate k for the limiting case of N — 0, i.e., when the test figure is vanishingly small. In this case there is a negligible probability of two or more intersections falling within the figure. There are therefore only two possible outcomes: a count of one with a probability of, say, p and a zero count with a probability of (I-p). Thus and Therefore: We might expect k to depend on N and therefore on the size of the test figure. Such a dependence was observed2 in the analogous case involving the estimation of grain-boundary area per unit volume by a count of the number of intersections of a test circle with traces of the boundaries on the plane of polish. In this case the theoretical value of k at the limit N — 0 was 1.41. It was also predicted that k would display an initial sharp decrease with increasing N. Experimentally, k was observed to decrease to a value of 0.62 at N - 6 and thereafter remain approximately constant. Furthermore, the value of k for N > 6 was found to be insensitive to the form of the structure provided no pronounced anisotropy or grain-size variation was present. In order to determinewhether k would exhibit a similar dependence on N in the estimation of grain-

Jan 1, 1967

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

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

Jan 1, 1965

By M. H. Gaskell, D. C. Lindley

Small, steeply inclined reservoirs without natural water drives often are found associated with salt domes or other highly faulted structures. Frequently, only one well may be economically justified in these reservoirs, from which oil is produced by an expanding gas-cap drive or by a dissolved-gas drive. Although a large fraction of the oil upstructure from the well may be recovered during primary depletion, only a small fraction of the downdip oil can be produced because of gravity effects. This downdip oil has been defined as cellar oil991 in contradistinction to "attic oil" which is the oil above the structurally highest well in reservoirs where a natural water drive is present. In this study, two methods of producing cellar oil by water flooding a one-well reservoir were investigated by means of a laboratory model. This model was scaled to represent a prototype reservoir having an area of 30 acres, a permeability of 500 md, and a thickness and dip of 30 ft and 30°, respectively. One method considered, the simultaneous injection of water and production of oil from a single well, would require that the lower portion of the productive interval of a well be completed for water injection and the upper portion for oil production. The other method considered, intermittent injection, would involve alternate injection of water into and production of oil from a single well. Injection of water into the formation would repressure the reservoir, causing gas to go back into solution. Since some of the injected water would flow updip from the well, water injection usually would be followed by a waiting period during which the water could drain from above the wellbore. Production of oil by gravity drainage and gas-cap drive could be initiated when water drainage became sufficiently complete. The amount of water produced along with the oil would naturally depend on how soon after water injection ceased that production was started. This process could be repeated in cycles until oil could no longer be produced economically. Of these two methods, the intermittent injection technique would be simpler to put into operation since it does not require a dual completion. Also, it is probably more generally applicable than the simultaneous injection-production method. It is to be expected that in some reservoirs the rock and fluid properties would be such that practical oil-production rates could not be achieved by the simultaneous technique. However, in reservoirs where the simultaneous method will result in practical rates of oil production, it does have the advantage of permitting immediate, continuous oil production; whereas in the intermittent method, oil production is interrupted by periods of water injection and drainage. It was the objective of this study to determine whether these two methods of cellar oil production are practical for the representative prototype investigated and to evaluate their relative merits. DESIGN OF LABORATORY MODEL SCALING REQUIREMENTS The prototype reservoir and fluids were selected to have the dimensions and physical properties shown in Table 1. The laboratory model was then designed in accordance with well known scaling criteria, as discussed by Croes, et al,2 except for one modification. According to Croes, et al, one requirement of scaling a laboratory model to a reservoir prototype is that both the relative permeabilities and dimensionless capillary pressure be the same functions of water saturation for the model and the prototype. Since unconsolidated sand packs are used in the laboratory model to represent the consolidated porous media of the prototype, it is virtually impossible to satisfy this requirement. Much of the difficulty stems from the fact that unconsolidated sands usually exhibit lower residual oil and connate-water saturations than do reservoir rocks. To satisfy the scaling requirement as nearly as possible, Perkins and Collins&apos; have suggested that the relative permeabilities be re-defined as

By Brown, R. W., G. H. Neill, R. G. Loper

All facets of ball sealer behavior must be known and understood to design for their optimum use in well treatments. The down-hole factors including the inertial forces, drag forces and holding forces have been discussed. Equations have been developed and presented for these forces. The use of these equations in this process will contribute to the most efficient use of an economical and proven ball sealer process. The urocedure presented has been used to successfully analyze field results indicating problems involving the use of ball sealers. These problems have been isolated and corrected in future treatments by using the procedure presented to revise the perforating program, the injection rate, the design of ball sealer stages and/or the type and size of ball sealers. The procedure presented also indicates that dislodgement of ball sealers occurs primarily in ultra slim-hole completions. INTRODUCTION The perforation sealing process has been proven highly successful and economical since it was introduced to the oil and gas industry in early 1956.&apos; Since it has proved so successful in multistage fracture treatments at a much lower cost than conventional packers, the ball sealer process has become a byword in well completions.&apos; Not only has it changed the design of fracture treatments, it has changed the concepts of selectively perforating and has greatly aided the success of the single-point and limited-entry techniques." Yet with all the advantages, there are still facets of ball sealer behavior that are not widely recognized. These must be known and understood to better design and execute the optimum ball sealer treatment. The down-hole behavior of ball sealers is influenced by a number of factors including an "inertial force", a "drag force" and a "holding force". These are considered in two phases—the initial seating of the ball and the forces tending to unseat the ball after contact has been made (see Fig. 1). This paper presents a theoretical study of these factors and formulas for predicting their influence. The authors acknowledge the fact that the flow equations presented herein describe the behavior of Newtonian fluids and do not cover the so-called power-law fluids. The assumptions made, however, are reasonable and well within practical limitations of the conditions encountered during field application. For example, it is assumed that the perforations are round with a discharge coefficient of 0.82. While this may not be exactly true, it has been verified within practical field limitations.&apos; The discharge coefficient would not vary significantly even if the perforation size or configuration varied pronouncedly. The approach of assuming uniform perforation size is not as strained as it might first appear. Normally, a perforator contains charges equally spaced with the same number in each vertical (or horizontal) plane. If the perforating gun is not centralized, the holes in one plane will be smaller than average, while the holes in opposing planes will tend to be larger than average. The net result is that the two effects tend to balance. FACTORS AFFECTING CONTACT WITH PERFORATlON Ball sealers are injected into well treating fluids for the purpose of contacting and sealing those perforations which are accepting the fluid flow. The efficiency of the sealers is primarily influenced by: (1) the velocity of the balls down the pipe, and (2) the fluid velocity through the perforations. To divert the sealer to the perforation, the inertial force of the ball must be overcome by the drag force created by the fluid velocity through the perforation. BALL VELOCITY The final, or stable, velocity of the ball is the sum of the actual fall velocity and the fluid velocity. This will be

By B. Haimson, C. Fairhurst

A criterion is Proposed for the initiation of vertical hydraulic fracturing taking into consideration the three stress fields around the wellbore. These fields arise from (1) nonhydrostatic regional stresses in earth, (2) the difference between the fluid pressure in the wellbore and the formation fluid pressure and (3) the radial fluid flow through porous rock from the wellbore into the formation due to this pressure difference. The wellbore fluid pressure required to initiate a fracture (assuming elastic rock and a smooth wellbore wall) is a function of the porous elastic constants of the rock, the two unequal horizontal principal regional stresses, the tensile strength of the rock and the formation fluid pressure. A constant injection rate will extend the fracture to a point where equilibrium is reached and then, to keep the fracture open, the pressure required is a function of the porous elastic constants of the rock, the component of the regional stress normal to the plane of the fracture, the formation fluid pressure and the dimensions of the crack. The same expression may also be used to estimate the vertical fracture width, provided all other variables are known. The derived equations for the initiation and extension pressures in vertical fracturing may be employed to solve for the two horizontal, regional, principal stresses in the rock. INTRODUCTION Well stimulation by hydraulic fracturing is a common practice today in the petroleum industry. However, this stimulation process is not a guaranteed success; hence, the deep interest shown by the petroleum companies in better understanding the mechanism that brings about rock fracturing, fracture extension and productivity increase. Geologists and mining people became interested in hydraulic fracturing from a different point of view: the method may possibly be employed to determine the magnitude and direction of the principal stresses of great depth. Numerous articles in past years have dealt with the theory of hydraulic fracturing, but they all seem to underestimate the effect of stresses around the wellbore due to penetration of some of the injected fluid into the porous formation.1-4 Excellent papers on stresses in porous materials due to fluid flow have been published5,6 but no real attempt has been made to show the effect of these stresses in the form of a more complete criterion for vertical hydraulic fracturing initiation and extension. This paper is such an attempt. ASSUMPTIONS It is assumed that rock in the oil-bearing formation is elastic, porous, isotropic and homogeneous. The formation is under a nonhydrostatic state of regional stress with one of the principal regional stresses acting parallel to the vertical axis of the wellbore. This assumption is justified in areas where rock formations do not dip at steep angles and where the surface of the earth is relatively flat. This vertical principal regional stress equals the pressure of the overlying rock, i.e., S33 = —pD, where S33 is the total vertical principal stress (positive for tension), p is average density of the overlying material and D is the depth of the point where S33 is calculated. The wellbore wall in the formation is considered to be smooth and circular in cross-section. The fluid flow through the porous elastic rock obeys Darcy&apos;s law. The whole medium is looked upon as an infinitely long cylinder with its axis along the axis of the wellbore. The radius of the cylinder is also very large. Over the range of depth at which the oil-bearing formation occurs, it will be assumed that any horizontal cross-section of the cylinder is subjected to the same stress distribution, and likewise that it will deform in the same manner. STATE OF STRESS PRIOR TO DRILLING OF WELLBORE At a depth D from the earth surface within the oil-bearing strata, it is assumed that the following are 1References given at end of paper.

By R. S. Wagner, H. Brown

Criteria for nucleation with the aid of twin planes during crystal growth from the melt are discussed. It is shown experimentally that bismuth crystals can grow from a subercooled melt in a twinned bismdtic form. A mechanism of growth based on the presence of at least me twin plane is proposed. It is well known that faceted growth can be observed in crystals grown from the vapor phase, dilute solutions, or from the melt. Faceted growth, in which the crystal is bounded by its slowest growing faces, will occur if the growth rate is strongly orientation dependent. The addition of atom layers on a facet plane can only occur by two-dimensional surface nucleation or with the aid of an imperfection. Nucleation with the aid of dislocations was proposed by rank&apos; and observed to be of particular importance in crystal growth from the vapor phase. Nucleation with the aid of two-dimensional imperfections, notably twin planes, was suggested by stranski2,3 and Frank.4 A twin plane, terminating at the growth interface, can form a reentrant step at which nucleation can occur at a much smaller super-saturation (or supercooling) than required for two-dimensional surface nucleation. Nucleation by a twin plane usually results in a fast growth rate along a crystallographic direction parallel to the twin plane. Such a nucleation mechanism, based on the presence of a twin plane, has been used by Dawson5 to explain the growth and morphology of sheet-like crystals of paraffin n-hectane. In recent years more evidence was found that nucleation by twin planes occurs frequently in crystal growth from the vapor phase and from dilute solutions, for example, in cadmium,&apos; barium titanate,7 silver,8 silicon9 and gallium arsenide.10 In crystal growth from a supercooled melt, however, nucleation by twin planes has so far only been observed in materials with the diamond or zinc blende structure. Bennett and Longini11 were the first to recognize the important role a single twin plane plays in the growth of germanium dendrites. Subsequently, it was found that a prerequisite for continued propagation of germanium dendrites is the presence of two parallel twin planes which form a self-perpetuating system of nucleation sites.12 Sili- con and some III-V compounds such as InSb also grow dendritically by twin plane nucleation. It appears that the following criteria must be satisfied in order that nucleation with the aid of twin planes can occur: A) The growth form of the material must be faceted." B) The formation of twin boundaries, by a fault in layer sequence during growth, by plastic deformation, or by an accidental meeting of properly oriented crystals, must be possible. C) Twinning must form a reentrant step where the twin boundary terminates at the faceted growth interface. D) These three conditions are necessary, but in certain cases not sufficient for continued crystal growth with the aid of a twin plane. Depending on interface morphology, a combination of more twin planes is sometimes required to form a self-perpetuating system of nucleation sites. This investigation was carried out with the idea in mind that nucleation with the aid of twin planes is a more general phenomenon and should occur in those materials, grown from a supercooled melt, which satisfy the criteria outlined above. It has been found that crystals of high-purity bismuth grown from a supercooled melt can grow by the aid of a twin plane. In the first two sections the faceted growth form and the twinning elements of bismuth are discussed. In the last section prismatic crystals of bismuth, which appear to grow with the aid of a deformation twin are discussed.

Jan 1, 1962

By Y. Kondo, A. M. Gaudin

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

Jan 1, 1964