Search Documents

By Paul B. Crawford, Melvin B. Burton

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

Jan 1, 1957

By J. Szekely, G. P. Martins

A formulation is given and a general analytical solution is presented for the problem of spherical phase growth in multicomponent systems. The results are valid for describing the rate of bubble growth in the absence of translational motion and, therefore, should be relevant to heterogeneous nucleation processes. Typical worked examples are given, which include bubble growth due to the simultaneous diffusion of hydrogen and nitrogen in molten steel, and the growth of CO bubbles in steel, due to the reaction: C + 0 = CO. The calculations predict that the presence of additional components would enhance the growth rate and also increase the absolute rate at which the individual components are transferred into the bubble. These findings should be relevant to degassing processes. THE problem of spherical bubble growth in the absence of translational motion is appropriate to the description of the heterogeneous nucleation and subsequent growth of gas bubbles while attached to a solid surface. Such problems are of considerable practical relevance to various vacuum degassing and steelmaking processes where bubble nucleation and subsequent growth play an important part in the overall reaction mechanism. A definitive formulation of spherical phase growth has been given in the now classical paper by Scriven1 and the analogous problem for solid phase growth has been discussed by Horvay,2 furthermore a thorough review of the relevant literature has also become available.3 All this previous work was concerned with one component or binary systems. The purpose of this paper is to show that the technique developed by Scriven may be generalized to alow the consideration of multicomponent systems. The results of this generalization would be particularly relevant, as gas bubble nucleation in the metallurgical systems previously mentioned frequently involves the simultaneous diffusion of several components. In the following a general formulation is given and typical worked examples are presented. FORMULATION Consider the growth of a spherical gas bubble due to diffusion in a quiescent liquid of infinite extent, as sketched in Fig. 1. The formulation of such problems requires the statement of the following conservation equations relating to the liquid phase: Equation of continuity Equation of motion Conservation of the diffusing species—Diffusion Equation. The following assumptions are made: a) The fluid is incompressible, quiescent, and of infinite extent b) The bubble is stationary, i.e., it undergoes no translational motion c) Spherical symmetry is observed throughout d) The mass diffusivity is constant e) Chemical equilibrium exists at the surface of the bubble, and any reactions occurring there are fast compared to the diffusion process f) Inertia and surface tension effects are neglected. Equations derived on the basis of these assumptions by previous investigators were found to provide a reasonable representation of bubble growth while attached to a surface in organic and aqueous systems4,5 For gas bubble growth in molten metal systems assumptions a) to d) chould hold reasonably well. The exclusion of chemical kinetics as a possible rate limiting step may not be universally valid, but should hold for many systems operating at high temperatures. The effect of hydrodynamic inertia may be neglected for moderate growth rates; finally, the assumption neglecting the effect of interfacial tension on the bubble growth was made for the sake of mathematical convenience; its inclusion would have required the use of quite cumbersome numerical techniques. While surface tension effects could be appreciable on the initial values of the overall growth rate, especially at low absolute pressures and in case of small bubbles, they are not thought to modify the relative effects present in multicomponent systems, the study of which is the objective of this paper. Within the framework of the above assumptions

Jan 1, 1970

By J. L. Moriarty, Irving Johnson, H. M. Feder

The crystallization of ternary alloys has been studied as a coprecipitation phenomenon by measuring the substitutional carrying of trace metals with CeCd11 precipitating from liquid cadmium. The Doerner-Hoskins logarithmic distribution law has been found to apply in all cases; the coprecipitation coefficients obtained were: 1.48 (La), 1.06 (Th), 0.63 (Pr), 0.24 (Gd), 0.17 (Sm), 0.13 (U), 0.10 (Er), 0.09 (Sr). Li, Na, K, Ba, Y, Sc, and Zr were not measurably carried. It has been observed that coprecipitation was most extensive for those tracers which form with cadmium an intermetallic compmdnd isostructural with CeCd11 or closely related thereto. For these tracers an inverse relation between Solubility in cadmium and coprecipitation coefficient was noted. THE partition of impurities between the liquid and solid phases of a metal or a congruently melting intermetallic compound has been extensively studied by the methods of zone melting.' However, when the composition of a liquid alloy and of the solid phase crystallizing from it differ grossly, partition of an impurity cannot be investigated in this manner. Partition of the impurity under equilibrium conditions may be studied by investigation of the appropriate phase diagram and partition under nonequi-librium conditions may be studied by observation of coring in the solidified alloys.2 The nonequilibrium case is of special interest in metallurgical practice because many complex alloys are cooled too rapidly to permit the uniform diffusion of impurities throughout the crystallizing phase(s). With this in mind, the present study was undertaken to obtain quantitative data on the phenomenon of coprecipitation via substitutional solid solution in certain metallic systems. It will be noticed that rather than adopting metallurgical ter- minology we have chosen to refer to the partition of impurities during alloy crystallization as coprecipitation.3 In the terminology of coprecipitation the impurity is referred to as the tracer and the precipitating phase as the carrier. From aqueous solution, coprecipitation may occur as a result of adsorption on the surface of the carrier, by occlusion of the liquid within the crystals, or by solid solution of the tracer in the carrier. By slow cooling of liquid metallic solutions it is relatively easy to form large, regular crystals of the carrier; the effects of adsorption and occlusion may be minimized and the influence of substitutional solid-solution formation may be made to predominate. Previous studies of the crystallization of some intermetallic compounds4 and of copper, silver, and gold5 have shown that the coprecipitation or partition behavior of trace metals is entirely analogous to that observed in aqueous systems when ionic substitution occurs, e.g., coprecipitation of Ra++ with BaSo4. The precipitation of CeCd11 from liquid cadmium was chosen for study. This choice was influenced by the following consideration. If the relative atomic sizes of the tracer atom and the atom for which it substitutes determine the tendency to coprecipi-tate, then a lattice having a substitutional site of well-defined size would be desirable. In CeCd116 (structure type cubic BaHg117) the cerium atom is completely surrounded by a cage of twelve nearly equidistant cadmium atoms. Thus, atoms much larger or much smaller than cerium should be poorly accommodated in this cage. The choice of a lanthanide central atom allowed the testing of the size hypothesis with a number of lanthanide tracers which differed in size but not in valence or electro-negativity. EXPERIMENTAL The general technique of this investigation was to prepare a slightly undersaturated solution of carrier in liquid cadmium to which a tracer was added. The concentration of tracer was always so low that saturation over the entire temperature range investigated was never exceeded. The solution was cooled in a stepwise fashion. After each step a por-

Jan 1, 1964

By H. Kurushima, S. Tsunoda

An investigation of the hydrometallurgy of Cu-Zn concentrates by The Dowa Mining Co. of Japan is described and a detailed flowsheet of the Kosaka hydro-metallurgical plant is given. Sufficient operating experience has shown that the Dorrco FluoSolids system has proven successful for recovery of both zinc and copper by roasting in a fluidized bed. At present, 2400 metric tons of concentrates per month are being treated and recoveries of 93 pct Cu and 65 pct Zn are being effected. A SURVEY of the Japanese chemical industry in 1949 revealed, among other things, that sulphur-bearing blast furnace gas from the Kosaka smelter of The Dowa Mining Co. Ltd., Kosaka, Japan, was being exhausted to the atmosphere as waste. In addition, it was costing that company an average of ten million yen ($28,000) annually for damages to the surrounding property caused by the sulphur dioxide in the gas. It was recommended by this survey that the gas be utilized for the production of sulphuric acid by the contact process. Mr. Kurushima visited the United States in May 1950 to investigate the engineering, construction, and operation of a contact acid plant. Any acid produced at Kosaka would necessarily have to be 98 pct H,SO, to prevent corrosion in railroad tank cars transporting it to the point of use. However, the production of 98 pct acid requires constant volume and SO, content of the gas. As the converter gas at Kosaka fluctuates in both volume and quality, it was decided to evaluate the smelting process with relation to other methods of recovery. This evaluation revealed a relatively new type of roaster, the Dorrco FluoSolids reactor, which was reported to be capable of roasting the copper sulphide concentrate in an excess of air with the production of an SO1 gas strong enough for acid production and a water-soluble copper sulphate calcine. The copper sulphate then could be leached from the roasted ore and electrolytic copper recovered from the leach solution in much the same manner as practiced by the Chile Exploration CO. for treating naturally occuring copper sulphate ores. Use of a reactor of this type would eliminate the need for a blast furnace and converter and, in addition, would effect savings in labor and fuel. The Kosaka ore is a microscopically fine mixture of copper, zinc, and lead sulphides. The lead and a small amount of the barite present can be separated by flotation, but any separation of copper and zinc proved extremely difficult. In blast furnace operation at Kosaka, zinc had been slagged off until a stockpile of about 12,000,000 tons containing 5 pct Zn had been accumulated. Mr. Kurushima felt that if any of the zinc in the concentrate was sulphat-ized and leached with the copper, it also could be recovered electrolytically. The Dorr Co., sponsors of fluidization technique in the metallurgical field, suggested that sulphate roasting be tested out on a laboratory scale and that, if successful, a commercial Dorrco FluoSolids system be constructed. They indicated that the results from such laboratory tests could be transposed to a commercial scale without pilot planting the roasting operation. It was therefore decided to send samples of flotation concentrate to the West-port Mill, laboratories of The Dorr Co., for testing. The objective of Mr. Kurushima's visit to the United States thus had changed from investigation of sulphuric acid plants to a fundamental examination of the hydrometallurgy of copper and zinc. Testing Program In February 1951, Mr. Tsunoda went to the United States to witness the application of the FluoSolids system to sulphate roasting at The Dorr

Jan 1, 1956

By Pablo Hadzeriga

Great Salt Lake brine was subjected, in laboratory scale, to conditions simulating solar evaporation. Solid phases and the variation in composition of the liquid phase throughout the potassium salts crystallization zone are presented. These experimental results are discussed and compare with stable equilibrium data of the quinary system Cl-, SO4, Na+, K+, Mg + + and H2O. The Great Salt Lake brine in Utah is a vast reserve of several potential chemicals which has attracted much attention in the recent past. From the standpoint of an ore reserve, it may be said that the lake contains a low grade ore which has to be concentrated by evaporation of water. Because of the fairly favorable weather conditions in the Great Salt Lake area, it seems reasonable to consider solar evaporation as an economical means of concentration. From the view of phase chemistry, the main constituents of Great Salt Lake brine make it belong to the five component system: C1-, SO4=, Na+, K+, Mg++ and H2O. Because of their low concentration, the presence of other elements may be regarded as having little effect on the equilibrium of the five main components. Equilibrium data for the quinary system of our interest has been available for many years. Using these data, it is possible to calculate how the brine of the lake will concentrate when subjected to solar evaporation. In practicing solar evaporation, the weather factor and supersaturation conditions cannot be ignored, but their introduction into calculations using equilibrium diagrams is almost an impossible task. With the experimental work presented in this paper, an attempt is made to compare the stable equilibrium data of the quinary system as applied to the case of the Great Salt Lake brine and the dynamic equilibrium when that brine is subjected to continuous solar evaporation. Sea water brine belongs to this quinary system when only the major constituents are considered. The study of this system has been carried on since the last decades of the nineteenth century. The reason for these intensive studies was to give a geological explanation to the underground potash deposits. AS a result, many theories were presented trying to relate the equilibrium data with the geological formation, paragenesis and metamorphism of these deposits. Two main approaches were followed for the equilibrium studies of this quinary system. One of them, and the best known to us, was the school started in Germany in the last century by J. H. Van't Hoff. The other one, which started about the same time, was the school of N. S. Kurnakov in Russia. Both approaches were different and as a consequence generated opposite theories about the geological formation of potash deposits. The Van't Hoff school studied the quinary system of sea water under careful laboratory control, resulting in stable phase diagrams. During these studies, Van't Hoff recognized the difficulties in crystallizing some double salts.' ". . . . we unexpectedly found that some compounds whose formation in the solutions under investigation was possible at 25°C, nevertheless totally failed to put in an appearance. These, not to mention the calcium salts, were leonite, MgSO4 . K2SO4 4 H2O; kainite, KC1-MgSO4 . 3 H2O; and kieserite, MgSO4 .H2O. Even exceedingly slow crystallization with the addition of the compounds themselves did not remove the condition of supersaturation in the case of these bodies." However, those salts (leonite, kainite, and kieserite) were included in the stable phase diagram at 25°C as shown in Fig. 1 using D'Ans data.3 The figure also shows the sequence of crystallization of the Great Salt Lake brine according to the stable diagram. In later studies by the Van't Hoff school, especially by D'Ans,3 Karsten,4 Autenrieth,5,6 attention

Jan 1, 1968

By M. R. J. Wyllie, H. W. Patnode

In the quantitative interpretation of electric log data it is essential to know the formation factor, defined as the ratio of the resistivity of the formation 100 per cent saturated with brine to the resistivity of the brine.' This concept implies that the formation factor for any given rock is constant irrespective of the resistivity of the brine or solution with which it is saturated. It is apparent that this can be true only when the rock matrix is wholly non-conductive. As reservoir rocks may contain clay or shale having an appreciable conductance one may expect these constituents to influence the determination of true formation factor. The formation factor may be theoretically derived from log data by dividing the resistivity of the mud filtrate invaded formation by the resistivity of the mud filtrate. Such a method of calculntion has implicit in it the assumption that the fluids in the invaded formation are wholly replaced by mud filtrate. that the solid matrix is non-conducting, and that the formation factor of the invaded zone is identical with the formation factor of the uninvaded zone. Since the validity of the first two basic assumptions regarding the insulating character of porous matrices and the displacement of the interstitial fluids in porous media do not appear to be fully established, the formation factors of a number of cores and clay slurries have been measured under conditions suitable for verifying the assumptions and the extent of any errors inherent in them. CORE STUDIES Apparatus and Procedure Resistivity measurements were made on short cylindrical permeable cores which had been previously mounted in lucite.2 The cores were clamped by means of insulated bolts between brass plates about 5 cm square and 1 cm in thickness. The plates were provided with a ring seal to prevent leakage between them and the core mount, and had a centrally located orifice so that fluids could be flowed through the mounted core. The brass mounting plates were used as electrodes and any current flow between them necessarily passed only through the constituents of the intervening core. Fluid saturation of the core, and fluid flow through the core, were effected by a centrifugal pump connected to the outflow end of the core, thus providing a differential pressure across the core of approximately one atmosphere. The procedure used was first to evacuate each core by means of the pump and then to saturate the core with an aqueous sodium chloride solution of known resistivity. The resistance of the core was measured and flow of the solution through the core continued for sufficient time to establish a constant resistance reading. After equilibrium had been reached with one solution, another sodium chloride solution of a different resistivity was flowed through the core and the flow continued until a constant resistance reading was again obtained. This procedure was continued until resistance measurements of the core had been made with a series of solutions within a salinity range from 1.0N NaCl to distilled water. The data obtained enabled the apparent formation factor of given porous rocks when saturated with fluids of different resistivity to be computed and also provided information regarding the ability of an invading fluid to displace interstitial water from a core. Results Berea Sandstone — The first core examined was a specimen of Berea sandstone." This core had been cut essentially parallel to the bedding for permeability studies. The core was clamped in the core holder and saturated with 1.0N NaCl solution and the resistance determined. It was subsequently flushed and the resistance measured when saturated in turn with 0.1N, 0.025N, 0.01N, .005N sodium chloride solutions and finally with distilled water. In flushing the core with each

Jan 1, 1950

By B. H. Caudle, R. L. Slobod, E. R. Brownscombe

An earlier publication has discussed three methods lor obtaining relative permeability data on small core sample-: and the apparatus and technique for the capillary pressure displacement method. This paper describes the apparatus and technique for the solution gas and the dynamic displacement methods and presents a routine procedure for obtaining oil-gas and water-oil relative permeability data. Theoretical and experimental considerations are presented to show that the end effect commonly associated with the dynamic flow mechanism is extremely small where constant rates are employed in the flowing phase. An effect of flow rate on relative permeabilities obtained by the dynamic system is found only when gas is one of the flowing phases and this effect is ascribed to a form of channeling in the capillary system. The apparatus and procedure* used to obtain relative permeability data with three phases flowing are described and some preliminary results of the use of this method are shown. INTRODUCTION In a previous publication from this laboratory there appeared a brief discussion of the concepts behind three basic methods for obtaining relative permeability data.' These three methods were called the capillary pressure displacement method. the solution gas displacement method, and the dy. namic displacement method — the names being suggestive of the type of process used for obtaining the desired saturation prior to making the permeability measurements. In the same publication, the apparatus and technique for the capillary pressure displacement method were described and some typical results obtained by this method were presented. This paper will present: 1. The routine procedure used to obtain complete relative permeabilitv data (In small core samples with two phasec flowing; 2. The experimental techniques for the solution gas and dynamic displacement methods for obtaining relative permeability data; 3. The results of some studies on the mechanism of fluid flow through consolidated porous media; and, 4. The preliminary results on the determination of relative permeability for a system in which three phases are flowing. TECHNIQUE AND APPARATUS FOR THE DETERMINATION OF RELATIVE PERMEABILITY WITH TWO PHASES FLOWING In order to obtain relative permeability data on a small core sample which will describe the various phases of reservoir depletion. it is necessary that the methods used duplicate, as closely as possible, the various mechanisms involved in rescbrvoir depletion. It is important that. in any type of relative permedbility determination. the core be originally at virgin reservoir conditions with respect to the fluids within the pore *paces. (for example. oil used to displace water from core until connate water saturation estahlished). and that the subsequent saturation changes in the core be analagous to tile reservoir depletion mechanism for which the data are to he used. PROCEDURE OF DETERMINING RELATIVE PERMEABILITY DATA WITH TWO PHASES FLOWING A standardized procedure has been set up for obtaining relative permeability data and other information on small core samples. This procedure essentially consists of 11 steps as shown below: 1. Determination of tile capillary pressure curve.'! This information is used to help select samples for relative permeability determinations and to charactrrize. in general. the pore size distribution of the sample. 2. Determination of the specific permeabilities to air, oil. and water. These values are those on which the relative pern~eahility calculations are based. 3. Reduction of the core to connate water by an oil drive using the high speed centrifuge. 1. Determination of the critical gas saturation value and the associated relative permeability to oil by the solution gas drive method. 5. Determination of the relative permeabilities to oil and gab in the region where both phases are flowing by the capillary presqure method. 6. Reduction of the sample to residual oil in the centrifuge by an air drive. and the determination of the relative permeability to gas at this point.

Jan 1, 1951

By J. L. Lytton, T. E. Tietz

The recovery of flour stress of recrylstallized are-cast molyhdenum after 10 pet prestrain at 150°C was studied for recovery temperatures of 725o 801o and 901oC. The yield phenomenon observed after recovery was such that the yield drop diminished and disappeared as recovery proceeded. A fractional flow-stress recovery parameter fr, was developed and it was shown that, independent of recovery ternperatzwe, equal f, values of about 0.16 produced superposition of stress-strain curves except for the initial yield phenomenon. The activation energy for low degrees oj recovery was found to he 32,600 cal per mole, in close agreement with a reported valzie for difjzdsion of carbon in molybdenum. For values of fr greater than about 0.15 where recovery was thought to predominate over strain aging, the activation energy varied from 43,600 to 82,500 cal per mole and increased with increasing recovery and recovery temperature. A 10 pet pre-strain shifted the ductile-brittle transition temperature, as indicated by tensile elongation, from 10°C. to about -22"C. A 100-hr recovery anneal at 801oC was shown to shift the transition temperature to about —5oC, the major shift occurring during the first hour of recovery. ThE recovery of cold-worked metals is accompanied by changes in the arrangement of the dislocations generated during cold work. In recent years, the use of transmission electron microscopy has allowed direct observation of some of the dislocation configurations formed, such as "tangles", networks, cells, subboundaries, and loops. These observations have recently been reviewed by Keh and Weissmann.' Accordingly, some understanding has been gained concerning the nature of the dislocation rearrangements involved. The relationship between recovery and the decrease of mechanical strength, however, has only recently been explored. Previous efforts of the authors were directed toward investigating the decrease of tensile flow stress of aluminum and Al-1 pet Mg during recovery2 and the effects of concurrent elastic and plastic strain on the flow-stress recovery of aluminum.3 In addition, an investigation of the effect of four different solute additions on the tensile flow-stress recovery of aluminum has been made.4 Many important metals, however, have the bee crystal structure, e.g., iron, molybdenum, tungsten, vanadium, columbium, tantalum, and chromium. These metals might be expected to behave differently during recovery from those with other crystal structures for two reasons. First, the previous and present investigations of aluminum and its alloys have shown that tensile flow-stress recovery is strongly affected by foreign atoms in the lattice. In the bee metals, interstitial impurities are usually present and might be expected to have various effects upon both flow behavior and recovery. For example, heterogeneous yielding might occur during prestraining. The diffusion of these interstitials can also relock dislocations during recovery so that a yield point might be obtained upon restraining (strain aging). Second, some of the bee metals (e.g., iron, molybdenum, and tungsten) undergo the familiar structure-sensitive ductility transition and exhibit brittle fracture at suitably low temperatures. It is of considerable interest to determine the effects of recovery on the transition temperature of such metals. The purpose of the present investigation was to extend our knowledge of flow-stress recovery to the bee metals and to study the effects of recovery on the ductile-brittle transition temperature. Unalloyed recrystallized are-cast molybdenum was chosen as a test material, since it exhibits a well-defined ductility transition near room temperature. As in the previous studies, the degree of recovery was measured in terms of the decrease of tensile flow stress. Prestraining and restraining were carried out at 150°C to avoid the possibility of brittle fracture. EXPERIMENTAL PROCEDURE The are-cast molybdenum sheet used in this investigation was obtained from the Climax Molybdenum Co. of Michigan. It was 0.060 in. thick by 10

Jan 1, 1964

By F. J. Fayers, J. W. Sheldon

A theoretical discussion is presented of the behavior of the interface between two fluids of different physical properties when displacement is occurring along a thin tilted bed. An approximate equation of motion is deduced for the motion of the interface, suitable boundary conditions are derived and the conditions necessary for the approximation to be valid are determined. Steady-state behavior and the circumstances under which a steady-state interface can develop are discussed. Finite-difference solutions of the equation of motion for two characteristic water-drive problems are presenled. The theory of displacement in a thin bed may be regarded as being the porous-medium-equivalent of the classical "shallow-water wave" theory of hydrodynamics. INTRODUCTION The solution of moving interface problems are important in predicting the behavior of oil reservoirs under various conditions of production. Kidder4,5 has obtained analytic solutions for the fingering problem and for another particular case of the motion of the interface between two immiscible fluids of unequal density in a porous medium. In the present analysis, we analyze the problem of the motion of an interface along a slightly dipping bed whose length is long compared with its thickness. The requirement of a thin bed is necessary in order to give a simple differential equation for the motion. The two fluids are of unequal density and mobility. In the derivation of the equation of motion, it is assumed that the lighter fluid is being injected at constant rate at the high end of the porous layer, while heavy fluid is produced at the low end of the layer. The analysis applies equally well, however, to the reverse procedure in which the heavier fluid is injected at the low end of the tilted layer. It will be shown that the behavior of the system is critically dependent on which fluid has the higher mobility and that the interface can in certain circumstances achieve a simple progressive steady-state shape. The fluids may be considered to be immiscible or miscible. It is assumed, however, that the thickness of the interface region is negligible in comparison with other dimensions of the system, and that the strength of the saturation or concentration discontinuity remains constant during the motion. In conjunction with Buckley-Leverett theory and areal sweep-out calculations, the tilted interface problem is one of the fundamental problems of reservoir engineering; it defines the very important underpassing and overpassing phenomena associated with water-drives, gas-drives and secondary recovery operations. The tilted interface problem has been dealt with previously by Dietz.2 In the present paper, we endeavor to derive the pertinent mathematical results more rigorously, and in doing so we obtain modifications in both the differential equation and boundary conditions for those given by Dietz. We also attempt to show the connection between these results and some of the recent work on the stability of interfaces. Finally, we attempt to demonstrate the practical importance of the thin bed theory by showing the difference in behavior of two characteristic waterflood problems. In the favorable mobility case recovery is improved by increasing production rate, while in the unfavorable case undercutting or overcutting is minimized by reducing the displacement rate. PROPERTIES OF THE EQUATION OF MOTION Fig. 1 illustrates a porous layer of constant thickness H and tilt 8 with Fluid 1 of mobility ?l and density pl below Fluid 2 of mobility ?2 and density p2. The fluids are incompressible and flowing with volumetric flux q along the layer. Both the slope a of the interface with respect to the bed and the curvature of the interface are assumed to be small. What is meant by "small" is indicated in the derivation given in Appendix A of an approximate equation of motion for the interface. Other symbols are defined in the "Nomen-

By G. Sachs, E. J. Riping

All ferrous alloys can be made brittle by straining at sufficiently low temperatures. However, the changes in mechanical properties for different ferrous materials with decreasing testing temperature do not appear to follow any universal law. In particular, complex effects of testing temperature have been observed if cold worked steels were subjected to tensile tests at low temperatures. EFFECTS OF PRESTRAINING AT A TEMPERATURE ABOVE THE TESTING TEMPERATURE ON THE FRACTURE CHARACTERISTICS‡ A considerable amount of data has been presented by a number of investigators on the fracture characteristics of various cold worked steels.' These data usually relate to a two-step procedure first used by Davidenkov,² consisting of stretching by tension to a certain strain or " prestraining" at room temperature and then completing the tensile test, or "testing" at the temperature of liquid air. On the basis of data of this type made available to date, ferrous alloys can be classified into two groups with respect to their fracture characteristics observed on testing at a low temperature after prestraining at room temperature. Very complex phenomena have been observed for all annealed steels. On the other hand, some observations indicate that steels in certain other conditions represent a fundamentally simpler relation. In considering the effect of prestraining at one temperature on the properties obtained at some lower temperature, it should be expected that any cold work reduces the ductility retained in proportion to the magnitude of the cold work. However, the test data presented to date indicate that only heat treated (quenched and tempered) steels appear to conform to this expectation, according to the very limited test data in Fig 1, presented by McAdam, Geil and Mebs. This series of tests, Fig 1, shows that the retained ductility became smaller the larger the prestrain. However, this decrease is less than the amount of prestrain or, in other words, less than the decrease in ductility at the prestraining temperature. These effects of prestraining are fundamentally identical with those observed by Bridgman3 when a steel was prestrained (in tension) under hydrostatic pressure and then subjected to a regular tension test. It appears, therefore, that the basic effect of prestraining under conditions which yield a higher ductility than the (subsequent) testing procedure, consists of gradually increasing the total ductility of the metal from the initial low value of testing to the higher value of prestraining. On the other hand, all annealed steels, if cold worked by tension at room temperature and then tested at a lower temperature, provided that some ductility is retained, showed a far more complex behavior than that discussed above. This relationship is exemplified by Fig 2, for three pearlitic steels.4,5 These particular steels were selected because of the various shapes of their ductility-prestrain curves. All three steels, and any other investigated so far, suffered a rapid decrease in retained (and consequently also in total) ductility when subjected to small pre-strains. Then, after exceeding a certain prestrain it was generally observed that the ductility recovered.* This rather complex behavior of annealed steels may be tentatively correlated with the presence of stretcher strains at small restrains. Heat treated steels do not exhibit stretcher strains. This fact might possibly explain their simpler prestrain-retained ductility relationship. The destructive effect of these stretcher strains may possibly be associated with the presence of triaxiality, which may occur in any highly non-uniform stress and strain state (and which is retained as residual stress after unloading). Even a small degree of triaxiality may then cause a reduction in ductility, or even embrittlement at a sufficiently low testing temperature. The discussion presented above deals with the effects of cold working at a high temperature on the fracturing

Jan 1, 1950

By R. L. Slobod, B. H. Caudle, E. R. Brownscombe

An earlier publication has discussed three methods lor obtaining relative permeability data on small core sample-: and the apparatus and technique for the capillary pressure displacement method. This paper describes the apparatus and technique for the solution gas and the dynamic displacement methods and presents a routine procedure for obtaining oil-gas and water-oil relative permeability data. Theoretical and experimental considerations are presented to show that the end effect commonly associated with the dynamic flow mechanism is extremely small where constant rates are employed in the flowing phase. An effect of flow rate on relative permeabilities obtained by the dynamic system is found only when gas is one of the flowing phases and this effect is ascribed to a form of channeling in the capillary system. The apparatus and procedure* used to obtain relative permeability data with three phases flowing are described and some preliminary results of the use of this method are shown. INTRODUCTION In a previous publication from this laboratory there appeared a brief discussion of the concepts behind three basic methods for obtaining relative permeability data.' These three methods were called the capillary pressure displacement method. the solution gas displacement method, and the dy. namic displacement method — the names being suggestive of the type of process used for obtaining the desired saturation prior to making the permeability measurements. In the same publication, the apparatus and technique for the capillary pressure displacement method were described and some typical results obtained by this method were presented. This paper will present: 1. The routine procedure used to obtain complete relative permeabilitv data (In small core samples with two phasec flowing; 2. The experimental techniques for the solution gas and dynamic displacement methods for obtaining relative permeability data; 3. The results of some studies on the mechanism of fluid flow through consolidated porous media; and, 4. The preliminary results on the determination of relative permeability for a system in which three phases are flowing. TECHNIQUE AND APPARATUS FOR THE DETERMINATION OF RELATIVE PERMEABILITY WITH TWO PHASES FLOWING In order to obtain relative permeability data on a small core sample which will describe the various phases of reservoir depletion. it is necessary that the methods used duplicate, as closely as possible, the various mechanisms involved in rescbrvoir depletion. It is important that. in any type of relative permedbility determination. the core be originally at virgin reservoir conditions with respect to the fluids within the pore *paces. (for example. oil used to displace water from core until connate water saturation estahlished). and that the subsequent saturation changes in the core be analagous to tile reservoir depletion mechanism for which the data are to he used. PROCEDURE OF DETERMINING RELATIVE PERMEABILITY DATA WITH TWO PHASES FLOWING A standardized procedure has been set up for obtaining relative permeability data and other information on small core samples. This procedure essentially consists of 11 steps as shown below: 1. Determination of tile capillary pressure curve.'! This information is used to help select samples for relative permeability determinations and to charactrrize. in general. the pore size distribution of the sample. 2. Determination of the specific permeabilities to air, oil. and water. These values are those on which the relative pern~eahility calculations are based. 3. Reduction of the core to connate water by an oil drive using the high speed centrifuge. 1. Determination of the critical gas saturation value and the associated relative permeability to oil by the solution gas drive method. 5. Determination of the relative permeabilities to oil and gab in the region where both phases are flowing by the capillary presqure method. 6. Reduction of the sample to residual oil in the centrifuge by an air drive. and the determination of the relative permeability to gas at this point.

Jan 1, 1951

By S. S. Marsden, S. H. Raza

An experimental study of the flow of line-textured. aqueous foams through Pyrex tubes is described. The foams range in quality F (ratio of gas volume to total volume) from 0.70 to 0.96 and behave like pseudoplastic fluids. At lower flow rates they exhibit laminar flow ad have apparent viscosities which increase with quality, and which cover a range of 15 cp to 255 poise for tubes of 0.25- to 1.50-mm radius rj. At higher flow rates a plug-like type of flow is developed, the extent of which increases with both T and ri. When the same foams flow through either open or packed Pyrex tubes, remarkably high streaming potentials ?E are often generated. These can easily reach 50v if nonionic foaming agents are used, but are at least an order of magnitude less for ionic foaming agents. A linear relationship between ?E and the pressure differential ?p is observed; this usually extrapolates to positive values of ?p at ?E of zero. The slope of the line increases with both F and ri. An equation was derived to describe the streaming potential of non- Newtonian fluids in circular tubes and was used to correlate experimental results. The calculated < potentials are of the right order of magnitude. INTRODUCTION Foams are both unusual and inniguing in their physical properties, and have been the subject of many scientific studies. However, present knowledge of foams is still fragmentary, specific and often contradictory. Apparent viscosity of foam is the physical property of greatest interest to both rheo1ogists and engineers. sibreel reported that the apparent viscosity decreased with increasing shear rate in a manner similar to some non-Newtonian fluids. Penny and Blackman2 reported that fire-fighting foams had both a limiting shear stress and a tensile yield stress. There is little doubt that some foams at least behave like non-Newtonian fluids, and have apparent viscosities considerably higher than those of either constituent phase. The high apparent viscosity of foam with its concomittant effect on mobility ratio and sweep efficiency no doubt prompted several attempts by research groups to use foam as a displacing agent in porous media. Based on recent experience, most of these groups probably succeeded in completely blocking fluid flow in the porous media and then abandoned their efforts. Two groups apparently found the successful combination of experimental parameters at about the same time.3,7 Others have recently added to our knowledge of foam flow in porous media and its use as a displacing agent.8-13 An experimental problem encountered by Fried3 was a transient blockage of foam flow in porous media when distilled water was used to prepare the foam-producing solution. Fried surmised that this was due to an electrokinetic effect and he eliminated it by using electrolytes in preparing foaming solutions. He also measured the streaming potential of a number of foams in capillary tubes which he found to be appreciably higher than those obtained when the constituent liquid flowed under comparable conditions. This paper presents results of a more comprehensive study of the streaming potential generated by aqueous foam flowing in both open and packed Pyrex tubes. It also adds to knowledge of the rheology of these foams as deduced from their flow behavior in open tubes. APPARATUS AND PROCEDURE A diagram of the apparatus used is shown in Fig. 1. Details of its construction, testing and use are described elsewhere.14 Careful selection of materials, extreme cleanliness and rather elaborate electrical insulation and shielding were necessary to obtain reproducible results (±5 percent). Both streaming potential and streaming current were measured with an electrometer. The design of the foam generator developed for this work is novel (Fig. 2). It will produce uniform fine-textured foams at pressures sufficiently

By W. C. Knopf, G. Samsel

For several years International Minerals & Chem-ical Corp. has used a radiation method to assay potassium content of products from potash and feldspar beneficiation. The procedure is rapid, accurate, and inexpensive and is easily handled, if necessary, by nontechnical personnel. Data are presented to show the accuracies that may be expected if the technique is used for control of pilot plants or large-scale plants. Natural potassium is composed of three isotopes —K, K and K. Five other isotopes, artificially produced, have been detected. In natural potassium only the K isotope, with a half life of 1.4x10 years, is radioactive. Its modes of decay are beta and gamma rays, positrons, and x-rays. Only the more abundant beta rays, however, are important in radiation measurements using geiger tubes. Barnes and Salley&apos; measured radiation of potassium salts in solution in 1943. Gaudin and Parnell2 first reported the use of solid salts arranged in a vessel surrounding the ß counter in 1948. Later, Scheel used the method to assay potassium content of various crude and finished potash salts. Experimental: Dip-type Geiger-Muller tubes (TGC-5) supplied by Tracerlab Co. were used in all radiation counting measurements. The tube was encompassed by a glass container about M in. larger in diameter than the tube. This container, which served as a sample holder, was closed at the bottom with a rubber stopper to facilitate removing the samples. Samples were not weighed, but about 40 g of sample were needed for each test. Electrical impulses from the tube were counted on either a Tracerlab scaler or a rate meter. To obviate the need for constant attendance, preset features were incorporated into the scalers, so that at a predetermined time or count the units were automatically deactivated. When rate meters were used, counting rates were read directly from instrument or from a recorder used in conjunction with the rate meter. Each tube was calibrated, according to the manufacturer&apos;s instructions, to determine the correct operating voltage to insure reproducible results. A linear calibration plot of percent K2O vs counting rate was also obtained for each tube by recording the counting rate plus background count of several standards of different K2O content chemically analyzed by the chloroplatinate method.&apos; These standards were of approximately the same particle size distribution as samples to be assayed. Specific background count caused by cosmic rays and other surrounding radiation was checked periodically by counting a sample of pure NaCl. Fig. 1 is a representative calibration plot. Assays were determined directly from the standard plot or calculated by subtracting the background count from the counting rate, then multiplying by the slope of the line. Chemical assays of the test samples were made by the Sherril method," a flame photometric technique, or the chloroplatinate method. Data: The radiation assay method has been used for three years by International Minerals & Chemical Corp. in the laboratory and in pilot plants processing ores containing potassium. These dry concentration plants process sylvinite and langbeinite ores at Carlsbad, N. M., and feldspar ores at Top-sham, Maine. Sylvinite ore is predominantly a mixture of syl-vite (impure KCl) and halite (impure NaCl). Contaminants include sulfates, such as langbeinite (2MgSO - K2SO4), kieserite (MgSO H2O), poly-halite (MgSO4.K2SO, 2CaSO, . 2H2O), and slimes. Langbeinite ore is predominantly langbeinite and halite with small amounts of sylvite, sulfates other than langbeinite, and slimes. Several variables, in addition to correct calibration and standardization of tubes, were found to be important in obtaining reproducible results. Particle size distribution of standard samples should be

Jan 1, 1959

By G. H. Enzian

THE effects of phosphorus and nitrogen on the properties and behavior of low carbon steels are important considerations to both the manufacturer and the user of such material. For one thing, these two elements are the principal points of distinction between the chemical compositions of Bessemer and open-hearth grades. Further, open-hearth steels themselves may vary considerably in nitrogen content depending on the raw materials and ..practice used. In addition, appreciable amounts of rephosphorized open-hearth steels are being used at the present time in some low alloy, high strength grades and in certain other cases. For these and other reasons, it was felt desirable to obtain more complete knowledge of the effects of simultaneous small variations in phosphorus and nitrogen content on such properties of low carbon steel as tensile strength, impact toughness, sensitivity to cold work, and strain aging. A long range investigational program designed to furnish such information, was carried out, and the present paper summarizes some of the data currently obtained. General Discussion The strengthening effect of relatively small amounts of phosphorus in steel has been recognized, and to some extent utilized, for many years. A number of investigators have pointed out, from time to time, that under proper conditions of control, phosphorus can be a useful alloying element in low carbon steels.1-3 Nevertheless, in spite of the volume of data on the subject, the use of phosphorus for this purpose has always been handicapped by widespread feeling that its presence makes the steel brittle, difficult to cold work, and unreliable at low temperatures. This belief on the part of some steel users is not entirely justified in many cases, and often stems from early experience with acid Bessemer steels which have appreciably higher phosphorus (and nitrogen) contents than basic open-hearth grades. Nitrogen in steel, and its effect on properties, has also received considerable attention in the literature.4 In 1906, Braune5 ascribed the brittleness of Bessemer rails to the nitrogen content, and proposed a schedule of nitrogen limits for steels for various uses. This was one of the first papers to recommend directly the use of nitrogen content as a quality control measure, but it did not meet with much favor or support; as a matter of fact, the Jernkon-toret was rather sharply criticized for publishing the paper. In the intervening years a number of papers have been published on the subject of nitrogen in steel. As a result, it is now generally recognized that nitrogen can have an important effect on the properties and behavior of steel and that the presence of this element should be considered in any systematic study of factors affecting steel quality. In 1932, Graham6 introduced the concept that the quality of various steels could conveniently be evaluated and compared in terms of "sensitivity", which was defined as "the characteristic reaction rate of a piece of steel toward brittleness or hardness." This concept formed the basis for a broad research program into factors causing variations in the quality of commercial steels. At an early stage of this work, it became apparent that the observed effect of either phosphorus or nitrogen on the sensitivity of steel depended, in no small measure, on how much of the other of these two elements was also present. Indeed, the combined effect of phosphorus and nitrogen in steel may often be more important than the specific effect of either. With a few exceptions, it is not possible to predict accurately the effect on steel properties of a change in the amount of either phosphorus or nitrogen in the steel without also specifying the approximate amount of the other element. It seems desirable, at this point, to discuss briefly the distinction between strain-sensitivity and strain aging as used in this paper. In the early stages of its development, the work-brittleness test7 (which will be described in more detail later) was used

Jan 1, 1951

By D. E. Baldwin

A method has been devetoped for predicting the produced concentration profile for a miscible slug in a five-spot pattern. The technique consists of dividing the five-spot into radial etements and applying an approximate radial solution of the dispersion equation to each element. To test the method, laboratory experiments were conducted in five-spot models. The models were 2-ft X 2-ft X 2-in. Berea sandstone slabs simulating one quarter of a five-spot. An aqueous system was used with tritiated water as the tracer fluid. Effiuent concentration profiles were obtained for slug sizes ranging from 0.7 to 10.2 per cent of the pore volume. Agreement between predicted and experimental profiles was excellent. A prediction was also made for a field tracer test. For this case the predictive technique was modified to account for stratification by use of a layer cake model. INTRODUCTION Dispersion in porous media is of growing interest to the petroleum industry because of the increasing importance of secondary recovery operations. One of the controlling factors in the recovery of oil by miscible displacement or by the use of waterflood additives is the degree of mixing between the fluids of interest. Consequently, one of the prerequisites of any recovery prediction is an adequate method of predicting dispersion. Also, any quantitative interpretation of waterflood tracer tests is dependent upon the ability to determine the dispersion of the tracer as it flows through the reservoir. During the past few years there have been reported in the literature the results of a large number of investigations, both theoretical and experimental, concerning dispersion in porous media. (An excellent review of these results has been presented by Perkins and Johnston." Only a few of these investigations, however, deal with non-lin-ear flow systems. Raimondi et al.9 have presented an approximate solution of the dispersion equation for radial, diverging flow. Lau et al. 7 and Bentsen 1 have presented laboratory data for radial systems and have shown that the dispersion can be adequately predicted by Raimondi&apos;s solution. Bentsen has also presented a small amount of data for a radial, converging system. Brigham and Smith&apos; have applied Raimondi&apos;s solution, with approximating assumptions, to a five-spot pattern to predict the behavior of a field tracer test. The afore-mentioned investigations have been concerned almost entirely with radial diverging flow and have neglected the effect of converging flow. For this reason the present investigation was undertaken to develop a method of predicting dispersion in a five-spot pattern including the effects of both diverging and converging flow. The method which was developed consisted of dividing the five-spot into diverging-converging radial elements and applying an approximate solution of the disperson equation to each element. DEVELOPMENT OF PREDICTION METHOD The development of the prediction method consisted of two phases: construction of a radial model of a five-spot pattern and derivation of equations to describe the dispersion in the radial model. FIVE-SPOT MODEL A radial model which simulated the flow characteristics of a five-spot pattern was constructed by dividing the five-spot into diverging-converging radial elements (Fig. 1). Because of symmetry the model represented only an octant of a complete five-spot. Each element had an included angle of l° (making a total of 45 elements) and a radius r, such that the sum of the individual element areas was equal to the area of the five-spot octant. The normal five-spot fractional flow behavior was determined by averaging the experimental five-spot data of Caudle and Witte,&apos; Dyes, Caudle and Erickson 5 and Fay and Prats. 4 The resulting curve is shown as a solid line in Fig. 2. To match this curve with the radial model it was

Jan 1, 1967

By K. Mukherjee

Instantaneous resistance change during martensitic transformation in Fe-30 at. pct Ni alloy has been studied using a low-noise electronic system. The smallest signals recorded in this investigation were of the order of 10 pv. These small signals did not show the initial increase reported by earlier investigators to support a double shear mechanism.&apos; A variety of pulse shapes and sizes was observed during "burst" transformations, and these signals were neither a linear function of the size of martensite plates nor of the volume fraction transformed. Even the smallest signals, presumably from single plates, showed complications resulting from elastic and/or plastic deformation, shock wave propagation, and so forth, which are inherent in the growth process of a martensite plate. The time of formation and interface velocity determined froni this experiment agree very well with the pioneering work of Bunshalz and Mehl.&apos; AN Fe-30 at. pct Ni (c0.005 pct C) alloy was used to study the time of formation of martensite plates. Martensitic transformations, unlike the nucleation and growth processes, are usually characterized by an extremely fast reaction rate. During the transformation, the crystal structure changes and this change is propagated through a finite volume of the parent crystal lattice within a very short time interval. The interface between the two phases probably consists of moving dislocations, and the time of formation of the plate depends on how fast this dislocation interface can move. Knapp and Dehlin~er or Cottrell and Bilby mechanisms, or similar dislocation models, can be used to calculate the time of formation of a martensite plate of a given size. In the past, conventional photographic technique475 as well as electrical resistivity techniques"~-" have been used to measure the time of formation of a martensite plate. The specific resistance of the martensitic phase in Fe-30 at. pct Ni alloy is lower than that of the austenite phase. Therefore if a current is passed through a small austenite sample in series with a large external resistance, then one would observe a voltage change across the sample as martensite is formed. Such dc measurements are shown in Fig. 1. The initial large voltage drop in Fig. 1, at the M, temperature, is due to a "burst" transformation where a large number of plates are formed in a cataclysmic manner within a very short interval of time. Even during the later stages of transformation as shown in Fig. 1, one can observe discontinuous voltage changes due to small bursts. In order to resolve the time scale of such discontinuous stages, the voltage signals are amplified and fed into a highspeed oscilloscope. In this investigation a low-noise amplifier was used for the first stage of amplification. The signal to noise ratio was further improved by the design of a pulse transformer. This transformer was capable of amplifying the signal by a factor of six. Voltage signals were studied over a wide time scale ranging from 2 sec to 0.1 psec. I) EXPERIMENTAL PROCEDURE A) Specimen Preparation. Specimens were cut from induction-melted Fe-30 at. pct Ni (<0.005C) alloy and cold-rolled to 0.2 by 1 mm cross section and an average length of 5 cm. During mechanical shaping, the specimens were strain-annealed at 950°C for 24 hr in vacuum. Finally the finished samples were vacuum-annealed at 1000°C for 6 hr and water-quenched to room temperature. The average grain diameter was 0.02 cm. B) Cooling System. The M, temperature of the alloy was -55°C. Fig. 2 shows the arrangement which was used to cool the specimen at a slow rate. The specimens were mounted inside a i -in.-diam hole in the 2-in.-long section of the copper specimen holder. This cavity was filled with silica gel. A dummy specimen was used under identical conditions for

Jan 1, 1969

By L. V. Pirela, S. M. Farouq Ali

Some interest has been expressed recently in the application of solvents in conjunction with a thermal drive, such as a steamflood. At least one field project of this type has been reported. This paper makes the first attempt to provide some of the basic information necessary for choosing a solvent for such an application. Results of displacements conducted at elevated temperatures are also discussed. This paper shows that, in the case of alcohol-hydrocarbon-water or brine ternary systems, the system miscibility may increase, decrease, or both increase and decrease within the same system when the temperature increases. The distribution coefficient was consistently found to increase in favor of the oleic phase with an increase in temperature, which is advantageous from the standpoint of oil recovery. Results of solvent displacements in a sandstone core showed that even a relatively small temperature effect in the favorable direction (increase in miscibility and distribution coefficient) can lead to a 50 percent increase in oil recovery. On the other hand, if the system miscibility decreases with the temperature, even though the distribution coefficient increases, the oil recovery at the higher temperature may not increase. INTRODUCTION Miscible displacement has been the subject of many investigations. Miscible-phase solvent flooding, employing solvents such as alcohols that are miscible with both oil and water, has received considerable attention. Investigations conducted by Gatlin and Slobod,l Taber, Kamath and Reed,2 Holm and Csaszar,3 and Farouq Ali and Stah14 have helped in understanding the mechanism of solvent flooding. While the use of alcohols and similar solvents as oil recovery agents is economically unattractive, these solvents nevertheless provide a means of understanding the mechanistic aspects of important miscible-phase oil recovery processes such as the Maraflood process (discussed by Gogarty5). As shown by Taber et a1.,2 the phase behavior of a solvent-oil-water system is of .utmost importance in determining the efficiency of the oil displacement. This paper makes the first attempt to obtain the ternary phase behavior data for nine alcohol-hydrocarbon-water/brine systems at elevated temperatures. Results of core tests at elevated temperatures are also presented which corroborate the effect of temperature on oil recovery, as judged on the basis of phase-behavior data. The use of solvents at elevated temperatures has been suggested earlier by Farouq Ali.6 At least one field test of this type has been reportedO7 With the widespread use of thermal recovery techniques, it is possible in some situations (especially in steamflooding) that the use of a suitable solvent in conjunction with the heat carrier may be economically feasible. This paper attempts to provide some of the basic information needed for judging such feasibility. LOCATION OF THE PLAIT POINT The phase behavior of a typical ternary system consisting of a solvent, hydrocarbon and brine (or water) can be represented by a triangular diagram such as the one shown in Fig. 1 (for temperature TI). Point Pon the binodal curve represents the plait point, being the limit of tie lines such as Y1 Y2 The compositions of the coexisting oleic and aqueous phases are given by points Y2 and Y1, respectively, for any mixture composition along Y1Y2. The characteristics of ternary diagrams have been discussed by Findlay8 and others. Taber et aL3 have shown that the position of the plait point in a particular ternary system plays the decisive role in determining whether or not theoleic phase would be displaced in a continuous manner during an alcohol displacement. The area above the binodal curve also is important, since it represents

Jan 1, 1969

By C. C. Martin, M. C. Fuerstenau, R. B. Bhappu

Experiments revealed that quartz could not be floated in conductivity water at any pH with a long-chained sulfonate as collector. Various cations, Fe+++, Al+++, Pb++, Mn++, Mg++, Ca++, are shown to function as activators when the pH is such that the cation hydrolyzes. As the solubility product of the various cation-sulfonates was exceeded at the concentrations involved, precipitated cation-hydroxy-sulfonates, e.g. Fe (RSO3)2OH must be functioning as the collector in these systems. A number of interesting and important phenomena were revealed during a study of the response of beryl to sulfonate flotation.&apos; Perhaps the most important observation is the effect of various cations on flotation response. Experiments showed that in the vicinity of pH 3.0, about three ppm Fe+++ increased recovery of leached beryl from about 60 pct to 100 pct, whereas about four ppm Ca+ + reduced the recovery from 60 pct to 2 pct. With the exception of ferric iron, all of the other cations investigated reduced flotation recoveries drastically at this pH. As postulated in the paper Sulfonate Flotation of Beryl,&apos; even with low concentrations of polyvalent cations and sulfonate (on the order of one ppm), cation-sulfonates are precipitated immediately in solution. That is, blueish-white clouds were noted to form after given additions of salts to known solutions of sulfonate. Furthermore, it was also shown that one-hundred-fold more ferric iron than necessary to form the precipitated cloud had to be added to effect complete flotation. In other words, some form of solid ferric sulfonate was found to function as the collector in these systems. The form of the compound was suggested to be Fe(RCO3)2OH and that bonding to the surface occurs through the hydrogen of the hy-droxyl. As the form of this compound was shown to be dependent on the extent of hydrolysis of ferric iron, it is to be expected that ferric iron will be unique in some pH region. Further, if the presence of hydroxyl in the neutral precipitated compound is of primary importance, then other cations should function similarly but at different values of pH, since they will hydrolyze at different concentrations of hydroxyl ion. To determine whether the premise of hydrolysis is correct, the study was extended to quartz systems in the presence of various cations at various values of pH. EXPERIMENTAL MATERIALS AND METHOD Sodium alkyl aryl sulfonate was chosen as the collector for the experimental work. The reagent, sup plied as a courtesy sample by the Shell Chemical Co., has the following physical properties: Physical form Solid — Finely ground Sulfonate content, pct wt 95 to 97 Molecular Weight 450 to 470 Number of carbon atoms in hydrocarbon chain 25 to 30 Melting Point, OF 250 to 260 Solubility Soluble in water with gel formation at concentrations higher than 25 pct by weight sulfonate This particular reagent was chosen because of its high active sulfonate content, long-chain length, and solid form. With the exception of the frother and sulfonate, all other chemicals were of reagent grade quality. Conductivity water was used in the experimental work. This water, made by passing distilled water through an ion exchange column, had an average measured conductivity of one micromho. Experiments were conducted at room temperature in a small glass flotation cell with the following procedure: 1) A redetermined amount of water and salt solution (e.g. CaC12 . 2H2O) were combined and the pH adjusted to a given value. 2) A given amount of sulfonate was added so that final solution volume was 130 cc. 3) One drop of Dow 250 frother was added. 4) Five g of quartz were added and the system conditioned for five minutes. 5) The pH of the system was measured (termed flotation pH).

Jan 1, 1963

By W. van der Knapp

A theoretical and experimental ana1ysis is given of the change in volume of a porous medium due to changes in external and internal pressures. The result enable one to deduce directly the effect of large increlents in stresses. It is shown that the three-dimensional representation of the volume of a nonlinear porous system is a cylindrical rules surface, the generating lines of Which have a slope according to the solid rock compressibility. Analytical expression for such a ,surface is very simple. It follows that the bulk and pore compressihilities depend only on the difference between pore fluid and external hydrostatic pressure, i.e., the effective rock frame pressure. From measurements on a number of sandstones as well as limestones it can be conduded that for practical values of this effective frame pressure the pore compressibility falls in the range between the compressibility of water and that of undersaturated crude oil. In general, pore compressibility is found to be higher the lower the porosity. For a particular limestone reservoir pore compressibility and porosity could he related by means of a simple empirical formula. Results from the described study find a direct applicxtion in material balance calculations and in prob1ems of liquid flow through porous media. Since the elastic constants of the rock bulk material enter into equatiorls for velocity of acoustic waves, the results are also of importance for calculating rates of wave propagation in these media. It is shown for sandstones with 1.5 to -10 prr cent porosity that laboratory measured porosities may differ from those lender reservoir conditions by about I per cent of their value. The difference for low porosity limestotnes can be of the order of 10 per cent. The subsidence caused by elastic deformation of reservoir rock is shown to he very small. General formulas relating pore and bulk volume to external and internal pressure can be derived, as has been shown by Biot1, Gassmann2 and Geertsma.2 Their theories describe the differential change in volume caused by a variation in pore fluid tension and external hydrostatic tension. It has been known for a long timei that the stress-strain relations of a packing of spheres are essentially noninear. Our experiments show that for almost all reservoir rocks which arc elastically deformable the relation between volume and loading is also a nonlinear one, resembling that of a packing of spheres very closely. The behavior of these media is often described by using the linear stress-strain relations for small incremental stresses and strains in the vicinity of a pre-stressed condition." The present study differs from these considerations in that it enables one to deduce the effect of large increments in the stresses directly. In analogy to the stress tensor for a homogeneous body the stress tensor for a porous body can be written in the conventional way: where the tensor components, o, and t, represent the total forces acting on the faces of a unit cube of porous material. This tensor can be split into two parts. a hydrostatic and a deviatoric part, as follows2,3,1 where the component of the hydrostatic stress system, a. equals The hydrostatic stress causes a change in volume whereas the deviatoric stress system gives rise to a change in shape only. One of the difficulties in describing changes of volume occurring in reservoir rocks is that the relations connecting these volume changes with external and internal stresses are essentially non-linear. A complete description requires the bulk volume, V,,, as well as the pore volume,V,,. to be known as a function of both the changes in pore fluid tension (v) and the hydrostatic tension on the outer boundary.

By D. S. Rowley, F. C. Appl

Assuming that rock behavior, during cutting with a single diamond, may be approximated by that of a rigid, Coulomb, plastic material, a theory of single diamond cutting action has been developed. Using this theory, the stresses on the diamond cutting surface and the components of the cutting force have been determined. Theoretical results agree reasonably well with available experimental data. The theory of cutting of single diamonds may be used as the basis for subsequent study of the perfonnance of surface set diamond drill bits and other surface set diamond cutting tools. INTRODUCTION For more than twenty years, surface set diamond bits have been used successfully in the mining and petroleum industries to drill and core rock formations that range from medium hard to the hardest known formations. Most extremely hard rock formations cannot be economically drilled or cored by any other means. In view of the growing use and importance of surface set diamond drill bits, it is reasonable to study the way in which these tools cut, and the factors that influence wear of the diamonds. Such knowledge will contribute to better bit design and will improve operating procedures. To analyze the cutting performance of these tools it is appropriate to consider the way in which an individual diamond cuts. Once this is established, the over-all performance of the bits can be determined by considering the total effect of the individual diamonds. In recent years there has been a growing interest in rock mechanics and there now exist several publications dealing with the deformation and failure of rocks. Notable experimental work has been reported by Bredthauer,1 Handin and Hager,2, 3 and Gnirk and Cheatham.4 An extensive collection of data pertaining to rock properties at low confining pressures has been given by Wuerker.5 Cheatham and Gnirk have recently presented an excellent review of the state of knowledge regarding analytical and experimental rock deformation and failure.6 Generally rocks are extremely weak in uniaxial tension, and weak or only moderately strong in uniaxial compression. Both tensile failure and compressive failure are generally brittle in nature. This means that there is little or no deformation or "flow" of the material before shear or tensile fracture occurs. Thus, until fracture occurs, rock behaves very nearly like an elastic material in uniaxial stress. Significant changes in behavior of rock occur under conditions of triaxial compressive stress. There is generally a marked increase in compressive strength and in ductility. It has been demonstrated that the failure becomes plastic (in a macroscopic sense) at sufficiently high values of confining pressure. Sharp wedge chisel indentation studies conducted by Gnirk and Cheatham indicate that many rocks behave in a plastic manner at relatively low triaxial compressive stress states.4 If "failure" under triaxial compressive conditions is considered to be either complete rupture or macroscopic plastic flow, then the general failure characteristics of a rock can be reasonably well represented by the commonly known Mohr envelope. This certainly is not a perfect failure criterion (since the effect of the intermediate principal stress is neglected) but it has been found to be adequate for most engineering purposes. Mohr envelopes usually are constructed graphically using experimental triaxial test data. For purposes of analysis it is convenient to represent the envelope in equation form. Among the several approximate forms appearing in the literature, a parabolic form is suggested by Cheatham,7 and an exponential form is suggested by Paone and Tandanand.8 The exponential form has been adopted for this study, and based on published data, appropriate parameters have been established for 19 rock types. A typical plot of the data compared with the approximate envelope, and a summary of the rock properties and

Jan 1, 1969