Search Documents

By K. C. Strebig, C. W. Schultz, A. A. Selim

The effect of some organic additives in diamond drilling of quartzite was investigated in the laboratory. The drill was designed to measure the rate of penetration, the thrust, and the torque and to record each as a function of the distance drilled. A statistical plan was followed throughout the investigation and an equation representing the bit performance was derived from the mathematical theory of reliability. The drilling torque was found to increase with the use of additives while the sliding coefficient of friction and the wear coefficient decreased. The energy per unit volume consumed was found to increase with the use of additives. A hypothesis explaining the effects achieved with the additives is given. INTRODUCTION For some time researchers have attempted to increase drilling efficiency by adding certain agents to the flushing medium. Much of this interest is the result of the observations of Rehbinder, Shreiner and Zhigach1 They postulated that "...in processes for the mechanical destruction of solids, a region of increased crack formation is created in the deformed layers adjoining that surface of destruction, forming the predestmction zone. The surrounding medium (for example, a liquid) which wets a given solid will penetrate the embryo microcracks of the predestruc-tion zone. Increasing the molecular affinity of this liquid for the surface of the solid considerably facilitates the deformation and destruction by softening the solid in the zone of increased crack formation. In this manner, the hardness or strength of the deformed body is decreased by the influence of the penetration of the liquid into the predestraction zone." This is known as the "Rehbinder effect". Rehbinder's work has prompted many investigators to examine the effect of additives on the mechanical destruction of rocks. Engelmann, Terichow and Selim reported that the maximum penetration of the rock by a diamond point in a solution environment was at the isoelectric point of the solution.* Robinson reported that his results did not appear to support the theory of Rehbinder.3 Joris and .McLaren suggested that additives with surface tension lower than water will be more efficient in cooling the diamonds.4 Their effectiveness in removing the fine cuttings from the face of the bit will also improve the efficiency of the system. They reported an improvement in impregnated bit life and penetration rates with the use of soluble oil additives. In contrast, Long and Agnew, drilling with surface set bits, reported that lubricating agents increase bit temperature and wear.5 A significant change was reported6 in the mechanical properties of rock with polar additives. It was concluded that this change in the mechanical properties of the rock in the presence of the additive is due to the rock hydrophilic properties. Ethylene glycol, glycerine and anionic detergent were selected for this investigation because their surface tension is lower than that of water. The surface tension of these additives for the concentrations studied was 42 to 75 percent of the surface tension of water. A low surface tension additive would wet the surface of the diamonds more easily and hence would be more effective in cooling the diamonds. Another characteristic of these low surface tension additives is their effectiveness in removing the cuttings and preventing the flocculation of the formed particles in the drilling process. This article summarizes the results of an investigation performed at the laboratories of USBM, Twin Cities Mining Research Center. The objective of the investigation was to study the over-all effects of some organic additives on the performance of diamond impregnated bits in atmospheric drilling conditions only and not in downhole pressured conditions.

Jan 1, 1970

By M. Sheffield

This paper presents a technique lor predicting the flow of oil, gas and water through a petroleum reservoir. Gravitational, viscous arid capillary lorces are considered, and all fluids are considered to be slightly compressible. Some theoretical work concerning the fluid flow in one-, two- and three-space dimensions is given along with example performance predictions in one- and two-space dimensions. INTRODUCTION Since the introduction of high speed computing equipment, one of the goals of reservoir engineering research has been to develop more accurate methods of describing fluid movement through underground reservoirs. Various mathematical methods have been developed or used by reservoir engineers to predict reservoir performance, 193-5,7,8 The work reported in this paper extends previously published work on three-phase fluid flow (1) by including a rigorous treatment of capillary forces and (2) by showing certain theoretical mathematical results proving that these equations can be approximated by certain numerical techniques and that a unique solution exists. DISCUSSION The method of predicting three-phase compressible fluid flow. in a reservoir can be summarized briefly by the following steps. 1. The reservoir, or a section of a reservoir, is characterized by a series of mesh points with varying rock and fluid properties simulated at each mesh point. 2. Three partial differential equations are written to describe the movement at any point in the reservoir of each of the three compressible fluids. All forces influencing movement are considered in the equations. 3. At each mesh point, the partial differential equations are replaced by a system of analogous difference equations. 4. A numerical technique is used to solve the resulting system of difference equations. Capillary forces have been included in two-phase flow calculations.1,4,8 The literature, however, does not contain examples of prediction techniques for three-phase flow that include capillary forces. Where capillary forces are considered, each of the three partial differential equations previously discussed has a different dependent variable, namely pressure in one of the three fluid phases. Therefore, three difference equations must be solved at each point in the reservoir. Where large systems of equations must be solved simultaneously, an engineer might question whether a unique solution to this system of equations actually exists and, if so, what numerical techniques may be used to obtain a good approximation to the solution. It is shown in the Appendix that a unique solution to the three-phase flow problem, as formulated, always exists. It is also shown that several methods may be used to obtain a good approximation to the solution. The partial differential equations and differenceequations used are shown in the Appendix. Matrix notation has been used in developing the mathematical results. Two sample problems were solved on a CDC 1604. They illustrate the type of problems that can be solved using a three-phase prediction technique. SAMPLE ONE-DIMENSIONAL RESERVOIR PERFORMANCE PREDICTION A hypothetical reservoir was studied to provide an example of a one-dimensional problem that can be solved. Of the several techniques available, the direct method of solution as shown in the Appendix was used. The reservoir section studied was a truncated, wedge-shaped section, 2,400 ft long, with a 6' dip. (A schematic is shown in Fig. 1.) This section was represented by 49 mesh points, uniformly spaced at 50-ft intervals. The upper end of the wedge was 2 ft wide, and the lower end was 6 ft wide. The section

Jan 1, 1970

By R. R. Goddard

By use of the frequency response method with a radioactive tracer, it was possible to determine fluid dispersion and distribution in a natural consolidated and an unconsolidated medium. Measurements were made in a linear flow system at oleic saturations of 69 per cent in the consolidated medium, and 100 per cent in both media. Dispersion and distribution were obtained by measuring the amplitude attenuation and the phase velocity of sinusoidal waves with a dual monitor apparatus. The gamma ray emissions permitted in situ measurements at any distance along the porous samples. One result of importance was that the effective diffusivity increased as the wave length increased. As a consequence, a dispersion coefficient appropriate for the injection of large slugs might exceed the value measured by use of small slugs. Since flow models based solely on fluid velocity and an effective diffusivity coefficient imply that the diffusivity should be independent of frequency, such representations were not adequate for the data of this study. A comparison was made with a capacitance model of porous media with deadend PVs, but even this model was not completely adequate. By using attenuation and phase velocity data, fluid dispersion can be predicted without postulating a differential equation satisfied by the tracer concentration, thereby eliminating the need of a complicated model to represent dispersion. INTRODUCTION The flow of similar miscible fluids through a porous medium can be fairly adequately described by two parameters: the average fluid velocity and the effective diffusivity.1-3 It has been pointed out recently, however, that significant discrepancies exist between this representation and the experimental data.4-7 An improved agreement can be obtained by introducing additional parameters based on the concept of dead-end pores. The purpose of the present investigation was to find out whether the frequency response method could be used to measure the relevant parameters. The method was used in the following form. A stream of fluid was flowed at a constant rate through a sample of porous material and the concentration of a radioactive tracer in the fluid was varied sinusoidally at a fixed frequency. The effects of flow through a porous medium are a decrease in the amplitude of the concentration wave and an increase in the velocity of the peaks of the waves above the average velocity. Attenuation and phase veIocity of the waves were measured as a function of frequency and fluid velocity. The simple two-parameter model implies that the diffusivity should be independent of frequency. Data reported in this paper show that the diffusivity decreases as the frequency increases. Hence, as shown also by many others, the two-parameter model is not completely adequate. Coats and Smith 5 used two additional parameters in their model: the volume of the dead-end pores and the rate of mass transfer between dead-end pores and the flowing stream. Their capacitance model of a porous medium containing some stagnant fluid, to which transfer occurs by molecular diffusion, did not explain the dispersion results of either the present study or of theirs. Instead, the capacitance effect can be better described as the result of extreme velocity variations within the pores of the medium, with transfer between the velocity zones by convection. EXPERIMENTAL PREPARATION OF POROUS MEDIA Data reported here were obtained with a natural Berea sandstone, and with two packings of 20/25 mesh sand. Pertinent information about these media is listed in Table 1. The sandstone core was first sealed with a non-

Jan 1, 1967

By Mark J. Klein

An internal-friction profile induced by nitrogen in Cr-35 at. pet Re was studied as a function of nitrogen concentration and heat treatment. From these studies, the solubility of nitrogen in this alloy was deduced. The equilibrium solubility of nitrogen in chromium does not appear to be substantially changed by the addition of 35 at. pet Re to chromium. However the precipitation kinetics of nitrogen from solution are significantly decreased. IT has been found that rhenium has a beneficial effect on the low-temperature ductility of Group VI-A metals. As part of a study to investigate the mechanism of this increased ductility, it was of interest to compare the concentration of interstitials retained in solution in one of these ductile bee alloys with the concentration retained in solution in its base metal. In this regard, a study of the solubility of nitrogen in ductile Cr-35 at. pet Re was selected for this phase of the investigation because the solubility of nitrogen in the base metal, chromium, has been determined over a fairly large temperature span.1-3 In addition, nitrogen is known to have an embrittling effect on chromium,4 but does not appear to have a comparable embrittling effect on Cr-35 at. pet Re.5 The experimental method selected for this study was the measurement of internal friction. This technique allows interstitial-concentration measurements to be made in bee metals that are sensitive to small changes in the concentration of a particular interstitial, even when other interstitials are present. When subjected to a cyclic stress, interstitial atoms in solid solution give rise to an internal-friction peak whose origin is the stress-induced ordering of interstitials into preferred lat- tice sites. These internal-friction peaks signify a dissipation of elastic energy which is greatest when the period of oscillation is comparable with the relaxation time for the attainment of an equilibrium interstitial distribution. For dilute interstitial binary alloys, the height of the internal-friction peak is proportional to the concentration of interstitial atoms in solid solutions. snoek8 has shown that the internal friction, Q-1, may be expressed by the equation where A is the relaxation strength, ? is the angular frequency, and 7 is the relaxation time. The relaxation strength, A, is proportional to the concentration of interstitials taking part in the relaxation if no interaction occurs between solute atoms. Thus values of peak height, which can be determined as a function of temperature, give a measure of the relative concentration of interstitial atoms in solution at various temperatures. If equilibrium has been attained, the peak heights will be proportional to the interstitial solubilities. The constant that relates A to the interstitial concentration can be determined experimentally from peak-height measurements as a function of concentration for specimens of constant crystallographic orientation. When this has been done, Q-1 measurements can be used to determine solubilities according to the equation C - KQ-1 where C is the soluble interstitial concentration and K is a constant that can be experimentally determined. It is possible, however, that lattice imperfections may bind interstitials to certain lattice sites. Internal friction does not measure these interstitial~, but does measure interstitials which are free to move in the normal lattice sites when subjected to an applied oscillating stress. On alloying the pure metal with a substitutional

Jan 1, 1965

By M. J. Pool, with Appendix by Larry Kaufman, R. W. Krenzer

Thermodynmic properties of the Cu-Mn system have been deterrrzined in the temperature range 973" to 1273°K by measuring the vapor pressure of manganese in equilibrium with alloys of compositions varying from 0.13 to 0.91 rnole fraction manganese using the Knudsen effusion technique. Activities of manganese calculated from the measured vapor pressures show large positive deviations. Calculated heats and entropies of mixing are positive with a maximum occurring at approximately 0. 7 mole fraction manganese. ThE investigation of thermodynamic properties of binary systems involving transition metals is of considerable interest at the present time. The incomplete subvalency electron levels of transition metals lead to unusual behavior in some alloys, and considerable information can be obtained by determining the thermodynamic characteristics of these systems and relating the thermodynamic behavior to observed physical properties. The Cu-Mn system, shown in Fig. 1,' provides an example of an alloy system in which the determination of the equilibrium phase diagram and various physical properties has received considerable attention, but very little work has been done on measurement of the thermodynamic properties. Elemental copper and manganese have comparable electronegatives, atomic sizes, and cohesive energies, but differ in crystal structure. They form a continuous series of solid solutions only at high temperatures and the structure varies in a continuous manner from fcc at the copper-rich end to fct at the manganese-rich end.' The purpose of this investigation was to measure the vapor pressure of manganese over a series of Cu-Mn alloys and from these values determine the thermodynamic properties of the binary system. Vapor pressure measurements were made in the temperature range 973" to 1273"K, remaining within the solid-solution region for each alloy studied. EXPERIMENTAL PROCEDURE The relative thermodynamic properties of solid Cu-Mn alloys were determined from vapor pressure measurements by means of the Knudsen effusion tech- nique. Several discussions of the theory can be found in the literature. In the case of copper and manganese, it is assumed that manganese is the only component subliming, so the vapor pressure of manganese over the various alloys can be determined directly by measuring the weight loss of the effusion cell. This is a valid assumption since the vapor pressure of manganese is about four orders of magnitude higher than the vapor pressure of copper in the temperature region of interest. The vapor pressure of manganese over the alloys at any given temperature is then related to the weight loss of the effusion cell per unit area per unit time. An upper limit of 1.5 x 10"s has been measured for the ratio of polyatomic to monoatomic molecules for manganese vapor in equilibrium with solid manganese at 1180°K, so it was assumed that manganese vapor is monoatomic. The effusion cells, fabricated from high-purity alumina, are 26 mm high, 18 mm OD, and have a wall thickness of 1 mm. Three different sized holes of 2.03, 2.54, and 3.18 mm diam were drilled in the lids from the inside face to minimize chipping, and the lids were then slip-cast onto the cells to provide a gas-tight seal. Alumina was selected for this particular application both because of its resistance to thermal shock and reaction with manganese vapor. In the ideal case, the orifice of the effusion cell should have a knife edge so that effusion through the orifice

Jan 1, 1970

By Milton C. Head, Aulton E. Roland, William E. Hawes

INTRODUCTION Uranium mining in the Grants Mineral Belt has undergone changes from methods pioneered during the 1950's and has now become fairly standardized throughout the district. The orebodies vary in size and shape and have been most successfully mined by using modified room and pillar methods with pillar recovery. This paper will expand on methods and equipment that are being used and modifications of these that will increase productivity and efficiency. GEOLOGICAL BACKGROUND The petrology of the Grants Mineral Belt consists of alternating layers of sandstone and shales (mudstone), with occasional thin layers of coal. The sandstones are aquifers which are generally artesian to some extent. The shales are bentonitic and decompose into extremely sloppy mud when exposed to water. Figure 1 shows a typical geologic section. Most of the uranium produced in the Grants Mineral Belt comes from the Westwater Canyon member of the Morrison Formation and occurs as scattered ore bodies throughout the formation. This formation consists of layers of bentonitic shale and sandstone of varying degrees' of cementing and grain size. It is also a very significant artesian aquifer. The unconfined compressive strength of the various formations is not high, as shown in Table I. It is important to note that these represent average values and a great range may exist within a single formation, especially In the Westwater which is the ore formation. [ ]

Jan 1, 1982

By C. D. Wiseman

THE concept of alloying to increase the strength of metals originated during the bronze age. However, at the present time there is no single theory capable of explaining all of the observed strengthening effects of alloying elements on the mechanical properties. In general, it may be said that solid solutions are stronger than the pure metal and that the strengthening effect increases with increasing alloy content. Various theories for solution hardening have been proposed, the first of which was presented by Rosenhain who suggested that the alloying elements roughened the slip planes, thereby requiring a higher stress to initiate slip. Recent work by French and Hibbard2 confirmed Norbury's3 conclusion that those alloying elements which cause a greater change in lattice parameter have larger effects on hardness. However, a given change in lattice parameter produced by different alloying elements does not cause the same change in hardness. A further refinement was suggested by Dorn, Pietrokowski, and Tietz4 who considered differences in valence between the solute and solvent with the change in lattice parameter to explain the plastic properties of aluminum alloys. Allen, Schofield, and Tate5 further investigated the valency contribution to solution strengthening by testing some copper-base solid-solution alloys having constant electron densities. Tensile tests performed on these alloys gave identical stress-strain curves even though the atomic size effect, which varied by a factor of two, was neglected. The advent of dislocation theory6-8 provided a new and powerful tool capable of explaining many phenomena of plastic deformation and indicated new approaches to solution hardening. The basic con- cept of recent theories is that solution hardening is not due to stresses associated with randomly dispersed solute atoms, but rather to a nonuniform distribution of the solute. This conclusion was reached by consideration of the effective tension present on each dislocation line and the ability of a dislocation to pass an obstacle. cottrell suggested the possibility that the internal stresses due to randomly distributed solute atoms would offer essentially no resistance to the passage of dislocations. He reasoned that an edge dislocation may be considered as a rigid line discontinuity. As the dislocation moved through a crystal, the random stress fields would be equally distributed along its length. Since half of the stress fields would be acting to move the dislocation forward and the other half trying to move it back, the net effect on the dislocation would be zero. Thus the randomly distributed solute atoms would not contribute appreciably to solution strengthening. Mott and Nabarro10 considered the dislocation as a flexible, rather than a rigid, discontinuity. They obtained a simple model for a dislocation moving in a lattice having precipated particles, where the dislocation would actually bend around the obstacle until it would loop back onto itself to rejoin and then continue onward. Their work indicated that the minimum radius of curvature of the dislocation line, under the action of an external stress and in the region of an internal stress, would be such that the dislocation would be unable to pass through a matrix containing small closely spaced precipates by the looping process. The dislocations would tend to remain rather straight, approaching the condition of Cottrell's rigid dislocation. Hence, the strengthening must not be due to random solute atoms; however, it might be due to the concentration of atoms, either as clusters of like atoms or as short-range ordered regions, or to the concentration of solute atoms around dislocations. To date, four mechanisms of solution strength-

Jan 1, 1959

By G. V. Jere, V. Krishnan, C. C. Patel

Standard free energy and standard enthalpy changes as a function of temperature have been calculated for the chlorination reactions of different oxides constituting columbite and tantalite. The tallies of standard lee-energy change (F) indicate the possibility of pvefeential chlorination of different oxides but the ease 0.f chlorination of columbium pentoxide at about 250°C and the chlorinating tendency of columbium pentachloride in promoting the chlorination of tantalum and other metal oxides, p7-evetzt the separation of columbiu?n and tantalum from other oxides and pom each othel.. It is likely that low temperatuves of chlorination and continuous 1err2oval of columbium pentachloride from the ..eaction zone may help partial separation of columbium from tantalum as well as from titanium. The standavd enthaply change (AH?) values shotr that 1ou.-temperature chlorination around 300" to 500°C is likely to sustain the chlorination , without external heating. COLUMBIUM (niobium) and tantalum metals, their alloys and interstitial compounds are finding extensive uses in nuclear reactors and high-temperature equipment.la-3 These metals also find uses in chemical process plants, electrical rectifiers and capacitors and as "getters" in electronic tubes.la With the development of Kroll process of reduction of volatilized titanium tetrachloride by magnesium for the production of titanium metal, chlorine metallurgy is receiving greater attention. Chlorine metallurgy has several advantages, the foremost amongst them are: the ease of chloride formation, separation of the different constituents of minerals by selective chlorination, and the reduction of the chlorides in vapor phase to the respective metals. Recently, Si-bert, Kolk, and Steinberg4 have considered a number of methods for the preparation of columbium metal and the most promising is found to be the halide reduction. Since the minerals columbite and tantalite are comparatively abundant in nature, the future of the metallurgical processes depends on the utilization of these minerals for the production of columbium and tantalum chlorides by chlorination. In this paper, an attempt has been made to interpret thermodynamically the chlorination of individual oxides constituting columbite and tantalite, with the help of the recent thermodynamic data.5-8 For this purpose, standard free energy and standard enthalpy changes of the chlorination reactions of the oxides with chlorine, both in absence and presence of carbon as the reductant, have been calculated as a function of temperature and represented graphically. Similar calculations for the reactions employing carbon tetrachloride and carbonyl chloride as chlorinating agents have also been incorporated in this paper. REPORTED WORK ON CHLORINATION Cuvelliezg has patented a process for the separation of columbium and tantalum from their concentrate, wherein 75 pct Cb05 is claimed to be chlorinated at 1050°C in the course of 7 hr using chlorine or chlorine containing nonreducing gases, the products of reaction being CbC1, or CbOC1,. At this temperature, Ta,05 is claimed to remain unaffected. It is further claimed1° that at higher temperatures, the proportion of Ta205 chlorinated increases while that of Cb205 diminishes. During the studies on the action of chlorine on individual metal oxides, Kangro and Jahn" have observed that Ta,05 is not chlorinated even at 1200" in the stream of chlorine, while 75 pct of Cb,05 is converted into CbOC1, around 900 to 1100°C. The chlorination of Cb205 by chlorine in the presence of carbon, investigated by Lind and Ingle,' and

Jan 1, 1962

By L. Yarborough, L. R. Smith

This paper presents results of a laboratory study of retrograde condensate recovery by revaporization into dry injection gas. Flow tests were performed in 10.6-ft long sand packs at 100F and 1,500 psi. In three runs methane revaporized the liquid from a n-pentane-methane mixture in the presence of immobile water. TWO of these tests were water-wet, and the third was totally oil-wet. In the three runs n-pentane recovery was complete after 2.5 hydrocarbon PV of injection. There was no significant performance difference between the two wettability extremes. In a fourth experiment, a methane-hydrogen sulfide mixture revaporized a synthetic light, sour condensate. No water saturation was present. Equilibrium compositions and volumetric data were obtained for the four-component condensate. The heavy component, n-beptane, was removed after 6 PV production. Comparison of the effluent fluid compositions with known equilibrium data shows that the flowing fluid was equilibrium vapor and that the mixing zone between equilibrium vapor and dry injection gas was short. Data indicated that complete recovery of retrograde liquid occurred after it was contacted by a sufficient quantity of dry gas. INTRODUCTION When pressure declines below the fluid dew point in a gas condensate reservoir, a liquid phase forms. In this process, referred to as retrograde condensation, the quantity of liquid formed is frequently small enough that the liquid is not a flowing phase. To prevent loss of valuable retrograde liquids, the process of dry gas cycling has been employed for several years as a more or less standard practice. In this procedure the reservoir pressure is maintained above the fluid dew point so that the liquid components may be produced as vapor and then separated at the surface. Although full pressure maintenance by gas cycling seems ideal in terms of preventing liquids loss, several factors can reduce the atuactiveness of such an operation. From a study of a condensate reservoir in Alberta, Canada, Havlena et a1.1 con-cluded that cycling under conditions of declining pressure leads to economic advantages and to a high recovery of hydrocarbon liquids. This study considered effects of volumetric sweep efficiency, retrograde behavior of the original wet gas and revaporization characteristics of the retrograde liquid when contacted by dry gas. The first major work concerning revaporization of liquid in a gas condensate system is that of Standing et a1.2 Calculations based upon the PVT behavior of a recombined gas condensate fluid indicated that all retrograde liquid can be recovered if it is contacted by a sufficient quantity of dry gas. The paper considered the effect of variable permeability upon the recovery of retrograde liquid. Standing et al. concluded that recovery of heavier components in the retrograde liquid is greatest if reservoir pressure is allowed to decline below the dew point prior to dry gas injection. Since the work of Standing et al., several laboratory studies have been reported which show that recovery of hydrocarbon liquids by vaporization into dry injected gas can contribute to increased recovery above that obtained by ordinary production practices. Vaporization from retrograde condensate,3 conventional oil4-7 and volatile oil8 reservoirs has been considered. There is little work that deals with revaporization recovery from condensate reservoirs. This study demonstrates that vapor-liquid equilibrium exists when a dry gas is injected into a porous medium containing retrograde condensate; investigates the effect of wettability on revaporization of retrograde condensate; and provides an indication of the amount of dry gas necessary to

Jan 1, 1969

By C. K. Eilerts

Finite difference equations were programmed and used to integrate the second-order, second-degree, partial differential equation with variable coefficients that represents the transient linear flow of gas-condensate fluids. Effect was given to the change with pressure of the compressibility factor, the viscosity, and the effective permeability and to change of the absolute permeability with distance. Integrations used as illustrations include recovery of fluid from a reservoir at a constant production rate followed by recovery at a declining rate calculated to maintain a constant pressure at the producing boundary. The time required to attain such a limiting pressure and the fraction of the reserve recovered in that time vary markedly with properties of the fluid represented by the coefficients. Fluid also is returned to the reservoir at a constant rate, until initial formation pressure is attained at the input boundary, and then at a calculated rate that will maintain but not exceed the limiting pressure. The computing programs were used to calculate the results that would be obtained in a series of back-pressure tests made at selected intervals of reservoir depletion. If effect is given to the variations in properties of the fluid that actually occur, then a series of back- pressure curves — one for each stage of reserve depletion — is required to indicate open-flow capacity and related flow characteristics dependably. The isochronal performnce method for determining flow characteristics of a well was simulated by computation. INTRODUCTION The back-pressure test procedurel is based on a derivation of the equation for steady-state radial flow of a gas, the properties of which are of necessity assumed to remain unchanged in applying the test results. The properties of most natural gases being recovered from reservoirs change as the reserve is depleted and pressures decline, and the results of an early back-pressure test may not be dependable for estimating the future delivery capacity of a well. The compressibility factor of a fluid under an initial pressure of 10,000 psia can change 45 per cent and the viscosity can change 70 per cent during the productive life of the reservoir. There are indications that the effective permeability to the flowing fluid can become 50 per cent of the original absolute permeability before enough liquid collects in the structure about a well as pressure declines to permit flow of liquid into the well. The advent of programmed electronic computing made practicable the integration of nonlinear, second-order, partial differentia1 equations2.3 pertaining to flow in reservoirs. Aronofsky and Porter4 represented the compressibility factor and the viscosity by a linear relationship, and integrated the equation for radial flow of gas for pressures up to 1,200 psi. Bruce, Peaceman, Rachford and Rice5 integrated the partial differential equations for both linear and radial unsteady-state flow of ideal gas in porous media. Their published results were a substantial guide in this study of integration of the partial differential equation of linear flow with coefficients of the equation variable. The computing program was developed to treat effective permeability as being both distance-dependent and pressure-dependent. In this study all examples of effective permeability are assumptions designed primarily to aid in developing programs for giving effect to this and other variable coefficients. The accumulation of data for expressing the pressure dependency of the effective permeability is the objective of a concurrent investigation. FUNDAMENTAL DIFFERENTIAL EQUATIONS DIFFERENTIAL EQUATION OF FLOW Linear transient fluid flow is represented by the partial differential equation,

Jan 1, 1965

By W. R. Kube, W. H. Oppelt

Five low-rank coals, including two lignites, a steam-dried lignite, and two subbituminous coals, were carbonized at 940°F, in a bench-scale carbon-ize~ with a nitrogen and hydrogen atmosphere, or both, using pressures from atmospheric to 1000 psig. Coking characteristics were observed, and yields, composition, and distribution of products from thermal degradation were determined. The influence of pressure on yields from coal carbonization, as well as on modification of coke structure, has long been of interest. One of the first attempts at pressure carbonization was that of Capps and Hulett1 who in 1917 investigated the effects of pressure ranging from 1 to 20 atm on yields at 600° C for bituminous and subbituminous coals. High-pressure carbonization has not been used commercially for coke production. More recently, carbonization at relatively high pressure in the presence of catalysts and hydrogen has been proposed for the production of metallurgical-grade coke from low-rank fuels.' Another approach involves carbonization at high pressure in a hydrogen atmosphere to produce heating gas and a solid residue that is gasified with oxygen and steam to produce the hydrogen required for hydro-genation 3-4 or, alternatively, is burned to raise steam. Recently it has been proposed to carbonize coal under pressure, burn the gas in a gas turbine to produce power, and use the solid residue as fuel in a conventional steam-power plant.5 Our primary interest in pressure carbonization stems from an attempt to determine the influence of pressure on the coking characteristics of lignite to design a high-pressure slagging gasifier using oxygen and steam. At atmospheric pressure, North Dakota lignite is generally considered to be non-coking. However, lignite and subbituminous coals have been reported to show increased agglomeration properties in steam, hydrogen, or carbon-dioxide atmospheres at elevated pressures, although no strong coke was formed at pressures up to 3000 psig.6 If lignite would agglomerate seriously under gasification conditions, proper flow of charge through the gasifier could not be maintained without installing a mechanical stirring device or destroying the caking properties of the coal before gasification. Data on the caking properties of lignite or other gasification fuel, as well as on yields of products under conditions simulating the carbonization process in the pregasification zone, were required in the preliminary design stage of the gasifier. No information was available concerning the influence of pressure and type of atmosphere on North Dakota lignite. To provide this basic data, the U. S. Bureau of Mines investigated the effect of pressure and a hydrogen or nitrogen-rich atmosphere on the yield of products and on coke formation in a bench-scale pressure carbon-izer. SCOPE OF INVESTIGATION Five coals were tested at 940°F in nitrogen and hydrogen atmospheres, at pressures from atmospheric to 1000 psig. Yields and composition of solid, liquid, and gaseous products were determined as a function of the process pressure, and attention was given to the physical appearance of the solid residue. Percentage of carbon converted, sulfur retained in char, and thermal balances were calculated. In this paper, primary consideration is given to carbonization in a nitrogen atmosphere and to data concerning the solid phase. The effect of pressure in hydrogen is presented in detail for only one lignite. More complete information on the response of the other coals in a hydrogen atmosphere will be published in a USBM Report of Investigations. DESCRIPTION OF EQUIPMENT Details of the bench-scale pressure carbonizer are given in Fig. 1. The body of the carbonizer was

Jan 1, 1961

By V. R. Garfias, C. O. Isakson

The production of petroleum in Mexico during 1930 was approximately 39,600,000 bbl., or 5,188,000 bbl. less than in 1929. The decline which has continued for the last eight years, will place Mexico in fifth or sixth rank as an oil-producing country. The production of the northern and southern fields near Tampico declined over 5,250,000 bbl., while that of the fields in the Isthmus of Tehuantepec moderately increased. The wildcat operations which have been conducted for some years near the international border along a 290-mile belt, extending roughly 170 miles northwest and 120 miles southeast from Laredo, Texas, have so far met with little encouragement, although the oil possibilities of the region warrant the belief that commercially productive fields will be discovered eventually. An interesting development during the year was the starting of construction of the oil pipe line from Palma Sola in the South Fields near sea level—over Sierra Madre at 8660 ft.—to Atzcapotzalco, a suburb of Mexico City at an elevation of 7350 ft. This line, about 146 miles long, has a diameter of 10 in. and an estimated daily capacity of 10,000 bbl., and with its terminal refinery at Mexico City will form a transporting and refining unit which should materially affect the distribution and marketing of refined products in the Republic. Development During the past year 148 wells were drilled, as compared with 214 in 1929 and 371 in 1928, the percentage of producers being 49 as against 47 in 1929. With the exception of one well in San Isidro, at the extreme southern end of the Golden Lane, the Jardln Paso Real area was the only sector of the old fields where new development could be considered as successful, all other producers brought in during the year being drilled in nearly exhausted pools to be operated as ('strippers." Northern Fields A total of 71 wells were drilled in the northern fields, of which 22 were producers; one well in Ebano was estimated at 5000 bbl. per day, while the general average of all producers was less than 300 bbl. Results

Jan 1, 1931

By Robert E. Green

Experiments were performed in order to investigate the influence of magnitude of driving force, recouery, and previous heat treatment on grain boundary migration in deformed aluminum crystals. The frequency of occurrence of strain-free grains possessing Preferred crystallographic orientations was studied both as a function of heat treatment and as a function of amount of deformation of the matrix crystal. Aluminum of such composition as to exhibit precipitation-hardening behavior was used. CONSIDERABLE work has been devoted to the factors influencing grain boundary migration in metals.1'52 Aust and Rutterl-9 have studied grain boundary migration into matrix crystals containing a striation substructure, which provided a thermally stable driving force during the course of the migration experiments. These authors state that, while variations in the striation pattern were sometimes observed, it was possible to obtain reasonably uniform and reproducible striation structure by careful control of the solidification conditions. Admittedly, this source of driving force is a unique one and has served well in numerous experiments. However, although the striation substructure provides a thermally stable driving force, the subject of its reproducibility from crystal to crystal is questionable. In fact, Aust and Rutter3 attribute the scatter shown in some of their experimental points and several anomalies in boundary migration rates to variations in the substructure driving force along a single specimen. A further limitation on striation-induced driving force is that the driving force cannot be varied over wide ranges in a controlled manner. More often used as a source of driving force is the strain energy stored in a crystal during plastic deformation. This type of driving force has the disadvantage that it is not thermally stable due to recovery taking place simultaneously with the grain boundary mobility measurements. On the other hand, this strain-induced driving force is very reproducible from crystal to crystal since it is very easy to grow single crystals possessing the same purity and crystallographic orientation and then to strain them identically. In order to measure grain boundary migration rates two different methods are primarily used. Most often the technique employed is the heat-cool-etch method1' which involves thermal cycling, while the other technique, developed by Graham and Cahn11 using an X-ray diffractometer and goniometer furnace, involves continuous heating only. Thus the influence of thermal cycling vs continuous heating is a factor to be considered in grain boundary migration rate measurements, particularly in cases where the source of driving force is the strain energy stored in the deformed matrix, and recovery is likely to occur. Numerous recent investigations have shown that extremely small impurity additions to zone-refined metals exert a marked effect on grain boundary migration rates1-5, 9, 12-27 and a few have considered the effect of the distribution of impurities.28-35 The concentration of impurities and their distribution as influenced by heat treatment must also be taken into account in migration-rate studies. In the present article consideration is given to the influence of the magnitude of driving force, recovery, and previous heat treatment on grain boundary migration rates. In order to vary the driving force over a large range it was necessary to use as a source of driving force the energy stored in the matrix crystal during plastic deformation. A majority of the investigations were performed on aluminum specimens which were of such compositions as to exhibit precipitation-hardening behavior.

Jan 1, 1965

By M. R. Tek, M. L. Katz

The flow of fluids in stratified porous systems, in which adjacent layers of the system are in communication over part or all of their common interface, in many instances involves crossflow between layers. A theoretical study has been made of the unsteady-state flow of slightly compressible fluids during depletion of bounded stratified porous systems. Consideration has been restricted to single-phase flow, with negligible capillary and gravity effects. Analytical expressions have been derived for pressure distribution and for (luxes at the producing face in two-layer systems with complete fluid communication between layers. These expressions occur in the form of double Fourier and Fourier-Bessel series. The method of extending the solution to any number of layers is demonstrated by deriving the equation for pressure distribution in a three-layer system. On the basis of the exact results for two-layer systems, certain generalizations regarding the behavior of stratified systems in general are presented. The performance of stratified systems, in terms of cumulative flux as a function of time, is shown to lie at all times between bounds established from single-layer theory. The upper bound is given by treating the system as a single layer with arithmetic average physical properties. The corresponding lower bound is the summation of the fluxes from each layer as treated individually. Variation in performance of stratified systems between the upper and lower bounds is determined by the extent to which crossflow occurs in the system. The principal parameters affecting the amount of crossflow are the vertical permeability of the system, the ratio of system thickness to "radius of drainage" and, for radial systems, the ratio of system thickness to interior radius. At very early times, the system performance is near the lower bound. As time and the radius of drainage increase, the per- formance tends toward the upper bound. INTRODUCTION One of the principal problems confronting the reservoir engineer today is the incorporation of the effect of reservoir heterogeneity into predictions of reservoir performance. The first part of this problem is the characterization of the actual reservoir heterogeneity in a mathematically useful manner.1, 2 It is then necessary to calculate the effect of this heterogeneity on the system performance. It is frequently possible to estimate this effect by comparison with the known behavior of an idealized system. One such idealized system is the stratified or "layer-cake" model. Many studies have been made of the behavior of stratified systems, from both a theoretical and experimental standpoint. The great majority of the theoretical studies have been concerned with systems of separated layers, with no fluid crossflow. Among the more noteworthy of such theoretical papers is that recently published by Lefkovits, et al,3 in which a comprehensive analysis was presented for the behavior of bounded stratified systems which are unconnected except at the wellbore. Three recent references have dealt with stratified systems including the effect of crossflow. Jacquard4 has presented mathematical derivations pertaining to depletion of two layer systems producing at constant rate. A doctoral dissertation by Katz5 contains analytical solutions and tables of results for pressure distribution and fluid flux during constant-terminal-pressure depletion of two-layer and three-layer systems with crossflow, together with tables of eigenvalues necessary in the computations. The dissertation also presents results from experimental studies of stratified systems using a heat-transfer analog model. This paper is largely based on the material in this dissertation. A third paper has been presented by Russell and Prats6 since completion of the work resented herein, containing mathematical derivations for and discussion of the behavior

By R. N. Thurston, E. A. Nesbitt, G. Y. Chin

A general relalionship between the amount of glide shear (due to slip) and the macroscopic shape change has been developed. Since the deformation can be large, finite strain analysis is employed. In this treatnzent, the shape change is expressed by a deformation gradient matrix, F - [axi/aXj1, where X and x refer to the initial and final positions of a particle. This matrix is readily evaluated in terms of the amounl of glide shear, a , and the direction cosines of the slip-plane normal and the slip direction for a single active slip system. In the case of deformation from slip on several systems, the product F(a) of the several deformation gradient matrices is first calculated. Then, by assuming that the final configuration is reached by a long series of small shears of magnitude a, occurring more or less alternately in the several slip sys- IN problems of crystal plasticity, it is often necessary to relate the amount of glide in the operating slip systems to the macroscopic strain components. This relationship is relatively simple when only small strains are considered. In this case the separate strain contributions from several slip systems are additive. Such a procedure is incorrect in the case of large plastic deformation. To the authors' knowledge, no general treatment of the latter problem has appeared in the literature. In specific instances, Mark, Polanyi, and schmid3 have derived the relationship between glide and axial elongation during tensile pulling of a single crystal when only one slip system is active. A similar relationship for duplex slip was worked out by v. Goler and Sachs.4 Taylor and Elam5 have also studied in detail problems of large plastic deformation. Due to uncertainties as to the slip systems at the time, terns, the final deformation gradient matrix can be obtained. Mathematically, the problem is reduced to finding the limit of F(dN as N — m while a — 0 in such a way that the Product Na = 0, a ftnite conslant designatzng the accunzulated amount of shear. It turns out that this limit is simply , where F, is the matrix whose elements are the coefficients of a in the F(a) matix. Application is made of the present treatment to a fcc crystal of Perlrzalloy compressed on the (110) plane and contrained to elongate zn the [112] directzon. In addition, it is shown that one rnay readily obtain from the general analyszs the well-known formulas relating elongation, amount of glide shear, and amount of lattice rotation for crystals deforming by single and double slip under tension. however, they were mainly concerned with proving that slip occurs on {111}(110) systems in fcc metals. The present general treatment arose out of recent investigations of magnetic anisotropy in cold-worked Fe-Ni alloys6 as well as some related problems in strength anisotropy of single crystals.' Detailed application of the general analysis will be provided in the case of a Permalloy single crystal compressed on the (110) plane and constrained to elongate only in the [112] direction. It will be seen that, as expected, the small strain approximation leads to significant errors after moderate straining. It will also be shown that the present general treatment yields the formulas of Mark et al. and of v. Goler and Sachs in those specific cases. GENERAL CONSIDERATIONS Consider a homogeneous deformation in which a material point initially at moves to both positions being referred to the same set of Cartesian axes. Such a deformation can be

Jan 1, 1967

By W. J. Hillenbrand

The soda-ash method for eliminating arsenic from a molten copper bath has been described previously in some detail.' Briefly, the process consists in blowing powdered sodium carbonate under the surface of the partially oxidized melt and in removing the slag formed, which will contain a certain proportion of the arsenic formerly in the bath. Sodium carbonate melts readily at the temperature of molten copper, and when molten it dissolves arsenic oxide rather rapidly. In the process, sodium arsenate is formed which is removed along with the excess of molten soda ash. In order that arsenic may be eliminated in this way, it must first be oxidized to the pentavalent condition. This, however, is easily accomplished by blowing air into the molten copper until a fair excess of cuprous oxide is present. The manner in which soda ash is introduced into the furnace has a considerable effect upon the efficiency of the process. In the earliest attempts at arsenic removal, soda ash—often in combination with slaked lime, for sodium carbonate is not the only solvent that will remove arsenic—was thrown upon the surface of the bath, was allowed to melt, and was removed after it had dissolved a sufficient amount of arsenic. At present powdered soda ash is generally blown in under the surface of the melt with the aid of compressed air. Melting immediately as it comes into contact with the molten copper, it can react with considerably more of the bath than by the former method, and thus can remove arsenic more rapidly and more completely. Since the resulting slag is removed at once, its very destructive effect upon the furnace refractories is minimized. The effectiveness of the soda-ash method may be gathered from the results of the following test. A 900-lb. charge of copper was melted in a small, oil-fired reverberatory furnace, and sufficient metallic arsenic was added to the bath to bring the concentration of that element to about 0.5 per cent. Air was blown into the melt until its oxygen content was just under 1 per cent by weight. The bath was then sampled and soda ash was blown in under the surface in 5-lb. batches. After each addition

Jan 1, 1934

By R. H. Hiltz, R. W. Douglass

Large-grain specimens of iodide titanium prepared metal-lographically were stain etched using the technique of New York University as modified by Watertown Arsenal Laboratories. Orientations of grains showing an abnormal response to stain etching were determined by back-reflection X-ray techniques. X-ray results showed that a-titanium grains, oriented such that their basal planes were parallel to the polished surface, have a slower reactivity to stain etching than any other orientation. Grains oriented such that their basal planes made small angles with the polished surface exhibited reaction rates between that of the true basal orientation and the normal rate. As the angle with the surface increases the reactivity increases up to a maximum angle of 14 deg. Basal planes making an angle with the surface of greater than 14 deg show normal reactivity to stain etching. These data, plus work on polycrystalline specimens, indicate that the oxidation rate of titanium might be reduced by correct alignment of crystal planes. ThE earliest observed phenomenon in the micrography of titanium was the occurrence of a grains which exhibited no birefringency. Even though a titanium possesses a hexagonal-close-packed structure and should be optically anisotropic, it can appear optically isotropic, when the crystal is so oriented that its intersection with the polished surface is perpendicular to the axis of symmetry. During the develoment of color staining techniques for titanium: it was noted that the a phase did not stain uniformly. Staining of specimens of commercially pure titanium resulted in a micro-structure of blue-colored grains, among which were randomly oriented a few grains exhibiting colors between the initial yellow and the final blue color. Under polarized light, these off-colored grains appeared to be optically isotropic. More intense study with polarized light, however, showed only the yellow grains to be isotropic. The intermediate colored grains exhibited weak birefringency, which at first glance is not detectable due to the strong birefringency of the surrounding structure. Aminoff2 and Straumanis3 have previously shown that in zinc the basal orientation is less reactive than other orientations. Leidheiser4 showed the same thing to be true for the hexagonal-close -packed configuration of cobalt. Thus, it would appear that the basal-plane orientation is responsible for the abnormality in stained titanium. Later work by Leidheiser5 on cobalt has shown that other planes besides the basal plane exhibited reaction rates slower than the normal rate. Reported planes were the (100), the (1011) and several planes making small angles with the basal plane. The last observance appears not to be due to the reactivity of the reported plane but due to the reactivity of the basal plane which is very close to being parallel to the surface. There are, then, two possible explanations for the range of colors observed in stain-etched a titanium. Since extensive research on the effect of orientation on surface reaction has failed to provide any definite trend, these two explanations are not adequate in themselves. Planes showing the retarded reactivity differ not only for different crystal structures but for metals of the same crystal structure and even for different reactions of the same metal. Tragert,6 Leidheiser,7 Tammann,6 and wagner9 have shown copper to exhibit retarded reactivity on the (111) plane. However, Leidheiser10 has shown that under other conditions the retarded reactivity of copper is associated with the (100)

Jan 1, 1960

By A. M. Gaudin

It is only in recent years that substantial progress has been made in developing a physico-chemical picture of the solid state. The molecular concept, so useful in dealing with gases, was naturally carried over to the field of crystalline solids; it has resulted in the firm implanting of erroneous concepts that only now are being eradicated successfully. Two simply written and effective books on the subject are Crystal Chemistry by Stillwell' and An Introduction to Crystal Chemistry by Evans. A more complete treatment of the subject is available in a book by Wells. The monograph by Hucke14 also touches on a subject most easily perceived by reference to the brief introduction by Bragg; that master of the field. These texts make it clear that while some solids are built of molecules, many are built of ions, others of atoms connected in the form of giant molecules, and still others of cations and electrons. The various lattices in which these entities are systematically arranged, the joint use of several of these entities, and the transitional character of some of the bonds make for the great variety of crystalline solids. The purpose of this paper is to point out a most important corollary of present views on the structure of solids. In this connection the classification of solids in Table I is helpful.

Jun 1, 1955

By B. L. Averbach, Morris Cohen, S. A. Kulin

The martensitic transformation can be initiated by elastic stresses at temperatures above M. in a steel containing 20 pct Ni and 0.5 pct C. Shear strains and normal tensile strains acting on a potential habit plane promote the transformation but compressive strains oppose it. THE strain sensitivity of the martensitic transformation has long been recognized. It is now well known that martensite formation can be induced by plastic deformation at temperatures below, and not too far above, the M, point.' However, the exact role of elastic vs. plastic strains and of the concomitant stresses has proved to be a very elusive matter. Many investigators2-' of the steel hardening reaction have made use of the stresses resulting from volume changes to account for such observations as the existence of retained austenite and its variation across the section of a bar, the change in the amount of retained austenite with quenching rate, and the self-stopping nature of the transformation when the cooling is stopped. The effect of applied strain as an independent variable has been studied in 70 pct Fe-30 pct Ni alloys by Scheil7 and McReynolds.8 Some data are also available for lithium and Li-Mg alloys,9 Cu-Zn,10 Cu-Sn and Cu-A1 alloys," and aus-tenitic stainless steels."-" It has been found that martensite formation can be induced isothermally, even at temperatures above M,, by mechanical deformation of the parent phase. The competitive nature of plastic yielding by slip and by the martensitic reaction has been depicted by Scheil.7 Basically, Scheil postulated that: 1—at temperatures not too far above M, (where austenite is less stable than martensite), a critical resolved shear stress (within the elastic range) is required to promote the transformation to martensite; 2—at and below M,, the austenite lattice becomes mechanically, as well as thermodynamically, unstable and shears over spontaneously into martensite without the application of external stress; and 3—at temperatures sufficiently far above M8, the critical resolved shear stress for martensite formation increases to a level above that required for slip, and plastic flow then supersedes the transformation when external stress is applied. Scheil's concepts were tested by McReyno1ds8 ho found no change in the elastic moduli (measured both statically and dynamically) of 71 pct Fe-29 pct Ni alloys in the vicinity of the M, temperature. Since the moduli did not approach zero or become negative on cooling to M,, it was concluded that the austenite lattice does not become mechanically unstable at M2, as postulated by Scheil. McReynolds also reported that M, is not raised by elastic stresses. Accordingly, plastic deformation was considered to stimulate the transformation by generating martensite nuclei in the distorted regions of the parent phase. The above issues with regard to the role of applied stress require clarification, if a basic understanding of the martensitic transformation is to be obtained. For example, in the reaction-path theory,'"," it is postulated that strain embryos exist in the austenite that provide part of the activation energy for the nucleation process. These embryos are visualized to be regions of localized strain in which the atoms in the austenite lattice are displaced part way along the path to their ultimate positions in the martensite. Consequently, it would

Jan 1, 1953

By F. I. Stalkup

Vapor-liquid phase equilibrium experiments have been conducted in a static equilibrium cell on mixtures of a light, 45 API stock- tank gravity reservoir fluid and a rich hydrocarbon gas containing approximately 55 mole per cent of intermediate hydrocarbons. Both apressure-vs-composition study of the gas and a simulated reservoir fluid, and a multiple-batch contact simulation of the condensing-gas-drive oil recovery process were performed. In the latter experiments equilibrium gas and liquid compositions were analyzed. Also, conventional, "condensing-gas-drive", long-tube displacement experiments of the reservoir fluid and gases of various richness were performed. The results of these experiments could not be satisfactorily interpreted by the conventional pseudo-ternary-diagram representation of multicom-ponent phase behavior. The results seem to be explained better by considering a bubble-point surface and a dew - point surface joined in a plait-point locus. Portions of the plait-point locus cannot be "seen" directly by the rich hydrocarbon gas because of curvature of the dew-point surface. In such a system, continuous injection of the rich gas over stationary reservoir fluid might form a zone of contiguously miscible compositions from pure rich gas to pure reservoir fluid by: (1) saturating the reservoir fluid with injected gas to the bubble-point surface; (2) creating by mass transfer with fresh injected gas a path of contiguously miscible compositions along the bubble - point surface to the plait-point locus; and (3) creating by mass transfer with additional injected gas a path of gas compositions along the dew-point surface up to the point where direct miscibility results between dew-point fluid and the injected rich gas. INTRODUCTION The use of the pseudo-ternary-phase diagram to illustrate miscible displacement phase behavior has been discussed by several authors.5,7 Such a representation of phase behavior is not rigorous, but the ternary diagram nevertheless gives a qualitative picture of what actually occurs in a miscible displacement process. Fig. 1 is a typical illustration of miscible displacement phase behavior by a ternary diagram. The multicomponent hydrocarbon system is divided into three pseudo-components: a light fraction containing methane and nitrogen, an intermediate fraction containing ethane through hexanes plus carbon dioxide, and a heavy fraction containing heptane and heavier components. A two-phase region is bounded by a dew-point curve and a bubble-point curve, which are joined at the critical point. The concept deduced from such a representation for miscible displacement by a condensing-gas-drive process is as follows: a rich gas G, which lies to the right of the limiting tie line through the critical point C, is injected into the reservoir and contacts reservoir fluid L, saturating the reservoir fluid to give bubble-point fluid L1 and equilibrium dew-point gas GI. Continued injection of rich gas changes the composition of the saturated liquid L1 through a series of liquid compositions lying along the bubble-point curve, until the critical composition C is reached, at which point direct miscibility with the rich gas is achieved. Some equilibrium gas with compositions lying along the dew-point curve from G1 to C is also formed in this process. Thus, continuous injection of rich gas forms a

Jan 1, 1966