Search Documents

By E. T. Turkdogan

Jan 1, 1962

By Ling Yang, John Henderson

Self-diffusion coefficients of copper in molten copper have been measured by the capillary reservoir method in the temperature range 1140o to 1260°C. The results can be represented by the equation D = [(1.46 * 0.01) x 10-3) exp [(-9710 * 710)/RT] cm2/sec. The energy of activation agrees within the limit of experimental uncertainty with that for viscous flow. The relationship between the diffusivity and coefficient of viscosity is best described by the equation D = kT/4pnr cm2/sec. MEASUREMENT was made of the self-diffusion coefficient of copper in molten copper in the temperature range 1413º to 1533ºK, using the capillary reservoir technique and CU64 as the radioactive tracer.The procedures for filling the capillaries, making the diffusion run and evaluating Dcu from the counting results were the same as described previously.1,2 All experimental operations involving the melt were carried out in a purified argon atmosphere. A graphite crucible, enclosed in a McDanel tube, was used as a container for the melt because of its satisfactory nature as a refractory material and the low solubility of carbon in copper. Graphite was also used for the capillary material. A McDanel thermocouple sheath, to which graphite radiation shields were affixed, was used to hold the capillaries. By moving this sheath through a rubber seal at the top of the McDanel tube, the capillaries could be lowered into or raised out of the melt. A Globar furnace with a constant temperature zone (± 2°C) of 2 in. was used to heat the assembly. The temperature of the furnace was the equilibrium value corresponding to a given setting of the transformer. Temperature was measured with a Pt-10 pct Rh/Pt thermocouple and varied during a run

Jan 1, 1962

Jan 1, 1964

By R. E. Gilchrist, R. F. Nielsen

When a gas is displaced by mother in a porous medium, and there is a relative immobile liquid present. there is a transition zone in which the gas composition varies from essentially that of the original gas to that of the injected gas. The change in extent of this zone with time may be governed by a mixing process in the direction of flow and by the rate of establishment of gas-liquid equilibrium, if the gases are appreciably soluble in the liquid. The relative importance of these two factors in deterri~ining the growth of the transition zone is discussed, with appropriate equations and experimental data. INTRODUCTION The various solution methods of oil recovery involve a transfer of the components of a mixture, whether within a single phase or between a liquid and gas. The similarity of a natural porous medium to systems used in separation and chromatographic processes makes it seem reasonable that certain concepts used by chemical engineers in "mass transfer" might be applicable to petroleum reservoirs. One approach, the HETP (height equivalent to a theoretical plate), involves a distance or length of time over which equilibrium might be con. sidered as established. Another approach involves direct use of a rate of mixing or rate of establishment of equilibrium. The latter approach has been preferred in relation to production problems. DISCUSSION AND PROCEDURE When a gas or liquid in a porous medium is displaced by one of different composition but miscible with it there will be set up a transition zone in which the composition varies from cssentially that of the original fluid to essentially that of the injected fluid. If the fluid being displaced is a gas, and there is a relatively immobile liquid present, the transition zone may be

By H. T. Kennedy, S. G. Dardaganian, J. L. Clanton

This paper presents equilibrium ratios measured on oil and gas samples taken from the Kelly-Snyder field at the reservoir temperature of 128°F. From these data smoothed curves were constructed for the reservoir fluid up to the bubble-point, 1,874 psia. A material balance calculation, based on the use of equilibrium ratios instead of the usual solubility and shrinkage meayr~rements is made on the field and compared to the con\.entional method. DETERMINATION OF EQUILIBRIUM RATIOS The equilibrium ratios, or K values were obtained by recombining separator samples of oil and gas taken when the reservoir was above its bubble point, and thus represented the reservoir fluid before any change had occurred due to the evolution of gas. The well sampled had a gas-oil ratio of 617 cu ft of high pressure sepalator gas per barrel of separator oil. and a bottom-hole temperature of 128°F. Samples were brought to equilibrium by agitating for 24 hours in a thermostat, maintained at the proper temperature within 0.2ºF, mercury being added when necessary to maintain the pressure constant. Samples of oil and gas were then taken by constant pressure displacement, and analyzed by low temperature frac-tionation. Equilibrium ratios were determined at 128°F and 500, 800, 1,000 and 1,400 psig. At the latter pressure, there were indications that the gas samples were contaminated with small quantities of oil, and only data taken at the first three pressures were included in Fig. I. This figure has been extended to higher and lower pressures by the use of the Fluor charts for 3,000 psi convergence pressure, on which the experimental points fit well. MATERIAL BALANCE EQUATIONS The possible errors involved in conventional material balance equations have been analyzed by Muskat1. Where the errors occur in field measurements, or in obtaining effective averages from variables such as gas-oil ratio and pressure, they will evidently be involved in any attempt to estimate the content of the reservoir. On the other hand, substantial errors may be introduced in conventional material balance calculations by: (1) the neglect of the amount of oil dissolved in gas at high pressures and temperatures in the reservoir, but which condenses under separator conditions; (2) the assumption that the separation of oil and gas in the reservoir may be duplicated by a differential liberation in the

Jan 1, 1969

By P. J. Closmann

Steam zone growth as a function of time has been calculated for the case of constant rate steam injection into a preheated reservoir. To simplify the calculation a linear temperature profile has been assumed in the cap and base rock at the start of steam injection. The results indicate that at early times augmentation of steam zone growth due to preheating should be greatest. At longer times the steam zone development becomes close to that calculated with no preheating. INTRODUCTION With the increasing application of thermal recovery processes to recover viscous oil, cyclic steam injection has become important in many large-scale projects. In this process, the cycle of steam injection followed by oil production is repeated a number of times. At the beginning of the second and later cycles, steam is injected into a reservoir that has already been heated but that has lost part of its heat both in produced fluids and by conductive heat loss away from the injection zone. In such a case, the temperature level of the injection zone and the temperature distribution of the surrounding rock will affect the growth of the steam zone developed during the subsequent steam injection. Knowledge of the size of steam zone developed is important in determining the amount of oil displaced and the extent of heating in the reservoir. It is also useful to be able to compute the size of the steam zone for cases where thief zones take most of the injected steam. This paper presents a fairly straightforward method of estimating the steam zone developed in a preheated formation, based on certain simplifying assumptions. Some cases of steam injection into a reservoir at its original temperature have already been considered elsewhere.1-3 Because it is difficult or impossible to obtain an accurate representation of the temperature distribution in the reservoir some time after initial heating has taken place, in this work, a linear temperature profile in the cap and base rock is assumed (Figs. 1A and 1B). The linear profile has the advantage of yielding heat flow equations that

Jan 1, 1969

By D. A. Elgillani, M. C. Fuerstenau, J. D. Miller

Adsorption of lead and ferric iron on quartz and alumina is presented as a function of pH. Only the hydrolyzed species of these metal ions, FeOH++ and PbOH+, adsorb significantly on each of these minerals. Zeta potentials of quartz were measured as a function of pH in the presence of various additions of aluminum, lead, and magnesium salts. Sign reversal occurs at the pH values where significant concentrations of hydroxy complexes are formed. Models of adsorption are presented in the absence and presence of collector. The mechanism of adsorption of metal and collector species in the activation of nonmetallic minerals is still the center of considerable discussion in spite of the extensive study given this system. Until recently, the view was widely held that hydrated metal ions were responsible for this phenomenon. This fact is surprising, since it had been shown that the pH range in which flotation is effected corresponds to that in which metal ion hydrolysis occurs.' Additional support of this concept was pro vided by the work of Wolstenholme and Schulman in another system.' Later work has shown that hydroxy complexes of metal ions are the active species. These studies have also indicated that basic aqueous complexes of metal-collector may be necessary for flotation.3-7 Depending on the molecular weight of the collector, formation of these complexes may occur before or after precipitation of the metal-collector salt has occurred. Ion pair formation can occur when low molecular weight collectors are involved because relatively high concentrations of metal hydroxy complex and collector can be present in solution before precipitation occurs. This is not the case when high molecular weight collectors are used. A model of aqueous complex adsorption has been presented6,7 in which it has been assumed that adsorption occurs by the formation of water from adsorbed hydroxyl and hydrogen of the hydroxyl contained in the complex. Other authors8 have suggested that hy-droxy complexes are involved in these systems but that collector adsorption occurs by replacement of the hydroxyl of the adsorbed metal hydroxy complex with a fatty acid anion. All of these investigations have involved mineral surfaces that were negatively charged so that it is not possible to state with certainty whether hydrated metal ions or hydroxy complexes are the active species in this system. Use of a positively charged mineral should be helpful, however, since the hydrated metal ion will be repelled from this surface. On the other hand, if the hydroxy complex adsorbs by splitting out water, surface charge should have no effect because of the free energy decrease that the system will experience when water is formed. In this regard, sapphire is well-suited for study, since its surface is positively charged until approximately pH 9. Electrokinetic studies of quartz in the presence of various metal ions should also be useful in establishing the mechanisms of adsorption in these systems. That is, if the hydroxy complexes adsorb by water formation, drastic changes in the zeta potential should occur in the pH range in which hydrolysis occurs. In this view, adsorption of ions of widely different hydrolysis characteristics on both sapphire and quartz was measured. In addition, zeta potentials of quartz were measured as a function of pH in the presence of various metal ions and compared with flotation responses obtained under similar conditions. Experimental Materials and Techniques Reagent-grade aluminum, ferric and magnesium chlorides, and lead nitrate were used in this investigation. Adjustments in pH were made with reagent-grade HCl or KOH. Conductivity water, made by passing distilled water through an ion exchange column, was used in all of the experiments.

Jan 1, 1971

By Frederick Laist

AFTER a series of experiments covering a period of about three years, ending from the spring of 19iJ.2 to the spring of 1915, a 2,000-ton leaching plant for the treatment of the accumulated copper concentrator tailing was built and put into operation. During the experimental period, the first step was that of laboratory experiments or beaker leaches. The results on this small scale were so satisfactory that a small operating plant, capable of handling 5 tons of roasted tailing per day, was installed. Again, the results proved satisfactory and an 80­ton plant was built and operated continuously from August until February, 1913. In this plant, full-sized roasting and leaching units were

Jan 8, 1916

By E. C. Stephens

The eastern Washington— Idaho clay region stretches along the northeast margin of the Columbia basin for 150 mil es. The three better known and developed clay areas are Clayton, Wash., (2 miles north of Spokane); Mica, Wash., (10 miles southeast of Spokane); and Moscow-Deary-Bovill, Latah County, Idaho (about 80 to 90 miles southeast of Spokane) see Fig. 1. Clays from the eastern Washington— Idaho clay area have been mined for more than 50 years. At present Mica, Washington deposits are used by the Gladding—McBean & Co. to make a variety of bricks and ceramic products and Latah County, Idaho clays are used by the A. P. Green Fire Brick Co. to make refractories. The Anaconda Co. made an intensive investigation of the clay deposits in the northeastern Washington-!daho area in 1919 to determine the existence of a supply of raw materials from which alumina could be made. The greater part of the work, which included systematic drilling and shaft sinking, was concentrated at Mica, Wash. This was part of a research project carrted on by the Anaconda Metallurgical Research Dept. under the direction of Frederick Laist and F. Frick and it was terminated in 4920. V. D. Perry, Vice President and Chief Geologist for The Anaconda Co., in 1952 assigned the writer to make a preliminary study of clay and other potential sources of alumina ore in the Pacific Northwest with particular emphasis on those deposits which might be of future importance to The Anaconda Co9s. then new aluminurri plant at Columbia Falls, Mont.

Jan 1, 1961

By B. S. Gottfried

This paper is concerned with possible transport mechanisms which occur during segregated burning (i.e., burning in an oil reservoir in which the oil-bearing formation is overlain by a "clean" porous zone containing only gas). ,Two analytical models are presented. In the first, burning is assumed to occur at the gas-oil interface, and the dominant mechanism is taken to be the rate at which oxygen diffuses to the interface. The second model assumes that combustion occurs within the clean zone, with the dominant mechanism being the mechanism of fuel deposition from the oil-bearing formation to the clean zone. ,Comparison of the results of the analyses with limited experimental data indicates that the latter mechanism is much more plausible. Both models are idealized representations of the actual process and are somewhat conjectural, however, owing to the lack of more definitive laboratory or field data. INTRODUCTION In recent years there has been considerable theoretical and experimental research directed toward understanding the recovery of crude oil by underground combustion. Although two different recovery schemes have been proposed, known as forward and reverse combustion, it is only forward combustion that shows promise of becoming a practical recovery technique.17 Most of the analyses of the forward combustion process have considered it to be pistonlike; i.e., oil is displaced by a vertical combustion zone which, with its ensuing gases, pushes die oil forward as the flame front is propagated through the reservoir.l-8,ll,l2,14-16 However, information gained in field pilot studies indicates that, at least to a certain extent, forward combustion proceeds in a segregated manner in oil reservoirs.10 An ideal segregated system is one in which the oil-bearing formation is overlain by a "clean" porous zone containing only gas. The location of the combustion zone in such a system is not immediately evident, although combustion is believed to occur either at the gas-oil interface or within the clean zone. In either case, the combustion zone may be horizontal and it moves in parallel rather than in series with the produced oil. The purpose of this paper is to explore the likely mechanisms by which segregated burning can occur, and to determine the sensitivity of the resulting models to changes in various physical parameters. Two governing mechanisms are considered: the case of a horizontal combustion zone at the gas-oil interface with the rate of combustion governed by the rate at which oxygen diffuses to the interface, and the case of combustion within the clean zone with the combustion process controlled by the transport of fuel from the oil-bearing formation to the clean zone. Both models are somewhat hypothetical and involve considerable simplification of the actual process. However, the results allow some preliminary understanding of segregated burning which is useful in the absence of more definitive laboratory or field pilot data.

Jan 1, 1967

By J. M. Weehuizen

Complete details are given of the drilling method by which two shafts were sunk at the new Beatm'x mine of Staatsmijnen in the Netherlands. The nature of the rock strata in the Dutch mining district determined choice of the drilling rather than the freezing method. Data for each step and details of the equipment used are given in the illustrations. Staatsmijnen was founded in 1902 with government financial backing, but it is run as a private enterprise, on a commercial basis, with the normal fiscal obligations. The company now has four mines in operation with annual outputs amounting to 7 1/2 million tons of coal. In 1952 it was determined that a new mine—the Beatrix—would be opened. Details of the shaft-sinking operation are described in this paper. SHAFT LOCATION The coal-bearing carboniferous rock in the Dutch mining district is overlain almost completely by soft water-bearing strata, which at the site of the new shafts are about 1600 ft thick. The sinking of the new shafts through this soft water-bearing overburden (Fig. 1) is a costly and time-consuming affair. From the two availaible sinking methods: 1) Poetsch's well known freezing method or 2) the drilling method of Honigmann-de Vooys, the latter was preferred, not only because of its relative cheapness, under the prevailing conditions, but primarily because by using the drilling process, a special type of shaft lining of steel, concrete, and bitumen could be used. The new shafts, located 330 ft apart, have an inner diameter of 18 ft, 4 1/2 in. The boreholes, which have a diameter of 25 ft and which were drilled through nearly 1600 ft of overburden and some 100 ft of carboniferous, each contain more than 29,000 cu yd—by far the largest in the world. DRILLING METHOD The drilling method is based on the use of mud-flush—a thin mixture of clay and water—to keep the borehole filled during the entire drilling period and during insertion of the lining (Fig. 2). The pressure of the mudflush against the borehole wall is greater than that of the water in the surrounding

Jan 1, 1961

By Eugene Stephenson

CRITICAL analyses of the work expended in producing oil by the utilization of gas energy have appeared in the publications of Shaw,1 Pierce and Lewis,2 and Herold,3 authors who have ably discussed the factors which must be taken into consideration in comparing energy consumed per barrel in various wells which are flowed by means of gas-lift. However, two distinct phases of oil production are concerned with the ratio of gas energy consumed to oil produced. One of these phases is the gas-lift, which utilizes the differential gas pressures for the purpose of flowing oil, whether these pressures be natural reservoir pressures or be induced by artificial compression of gas or air. The other phase has to do with prorating oil production in pools which are operated either as units, or are subject to such restrictions as those imposed at Oklahoma City by the Corporation Commission, or, by similar controlling groups in other fields. An abundance of literature is available relative to the problems of flowing wells by gas-lift, but the first suggestion that the energy of the gas within the oil and gas pool might serve as a basis for proration, appeared in the timely article by McWilliams.4 Much of the. dissatisfaction with proration in certain fields is because some operators feel that besides the hazard to the wells which accompanies the test, the method of allocating production on the basis of the "potential" gage is inequitable and gives some operators a distinct advantage over others.5 The most severe penalty may be exacted from the most progressive operators or those who attempt to produce their oil with a minimum quantity of reservoir gas. If :no difficulties ever

Jan 1, 1931

THE 125th meeting of the Institute was held in New York, Feb. 20-23, 1922, inclusive, and was the most successful annual meeting of the Institute ever held; there was a larger registration, there were more sessions of importance, and a greater number of sessions, than in any previous year. The total registration was 1331, comprising 351 New York members, 393 visiting members and 587 guests. At the annual business meeting of the Institute, held on Tuesday, Feb. 21, with President Ludlow presiding, there was a much larger attendance than usual. The report of the tellers was read and the following officers were declared elected: Arthur S. Dwight, Director and President-District 0 J. V. W. Reynders, Director and Vice-president-District 0 A. R. Ledoux, Director and Vice-president-District 0 Charles F. Rand. Director-District 0 George D. Barron, Director-District 0 William H. Bassett, Director-District 1 Thomas B. Stearns, Director-District 8 William Kelly, Director-District 5

Jan 3, 1922

By J. D. Rice, S. C. Pitzer, C. E. Thomas

Interpretation of pressure build-up data obtained in the conventional manner has often been difficult because of the deviation from theoretical behavior. Major causes of this deviation have been attributed to damage and afterflow, and to fluid redistribution in the wellbore which, in extreme cases, can result in a pressure hump in the early portion of the build-up curve. Theoretical investigations show that bottom-hole pressure is definitely influenced by phase redistribution in the tubing column during surface shut-in tests, and that the magnitude of the effect is sensitive to the producing gas-oil ratio and stabilized rate of flow in the well. For field experiments a wire-line tubing packer, which can be run in the tubing against a stabilized flow rate, was developed for bottom-hole shut-in tests. By use of this bottom-hole .shut-in method, pressure humps previously observed in surface shut-in tests were completely eliminated and the effects of afterflow minimized. Analog and digital computer studies have been made to obtain theoretical curves for comparison with field curves, and remarkable agreement between the results of bottom-hole shut-in tests and theoretical curves has been obtained. INTRODUCTION Pressure build-up data from shut-in wells have been used by the petroleum industry to determine the permeability of the formation, to estimate wellbore damage, to evaluate the static reservoir pressure and to speculate reservoir volume. Calculation of these factors is based on methods of analysis developed for theoretical systems whose build-up curves have a charac- teristic shape when the wells are shut-in at the sand faCe.1,2,3,5,8,10 However, field curves obtained by conventional surface shut-in methods do not always exhibit this characteristic. Such factors as stratification, rock heterogeneities and irregular reservoir geometry can cause the character of a build-up curve to deviate from that predicted by theory for a simple system. In addition, all field build-up data are affected by the methods utilized in obtaining the measurements. Conventional surface shut-in techniques allow two effects to occur which contribute directly to the manner in which pressure builds up at the sand face. The first, afterflow, has been recognized for some time and methods are available for estimating the magnitude of its effect and for correcting its influence.' The second, phase segregation in the tubing column after shut-in, has been reported only recently and appears to exert a considerable influence on the character of field build-up curves.!' In extreme cases, phase segregation can produce "pressure humps" in the early portion of the build-up curve. Such curves have been considered anomalous and have defied analysis. It is the purpose of this paper to discuss the effect of phase segregation in the tubing string during build-up and to present a bottom-hole shut-in technique for obtaining field build-up data which minimizes the influence of afterflow and eliminates phase segregation effects. It will be shown that reliable reservoir information can be calculated from build-up measurements obtained using this method, even though conventional surface shut-in tests on the same wells yield anomalous data. DISCUSSION OF PHASE REDISTRIBUTION EFFECTS Some of the effects of phase redistribution on pressure build-up curves have been described previously by Matthews and Stegemeier" who have presented evidence that phase segregation is responsible for the pres-

By L. E. Lynd

Ilmenite concentrates prepared from beach and dune sands of Florida, India, Australia, and other localities deviate to varying degrees from the theoretical composition of pure ilmenite, and generally contain more TiO2 than the formula FeTi03 would indicate, and some Fe2O3. This variation in composition is due to different amounts of altered ilmenite in the concentrates. (') The alteration has been shown to take place by oxidation and leaching of the iron in ilmenite, to yield residual products consisting of amaphous or very finely crystalline titanium dioxide and ferric oxide. Other studies(2) have stressed that ilmenite may alter first to an amorphous titanium-iron oxide, or to an amorphous mixture of titanium and iron oxides, with later alteration to crystalline titanium dioxide. There has been some speculation as to whether alteration of ilmenite is usually caused by weathering or by hydrothermal pr0cesses.(3,4,5 Three main methods of approach were used to study the mechanism and rate of ilmenite weathering. (1) Etch tests were used to learn the relationship between crystallographic orientation and the rate of chemical attack on ilmenite. (2) Altered ilmenites in New Jersey sands were compared with ilmenite in rocks of the area from which the sands were probably derived. (3) Experiments were carried out to determine at what rate, and by what agents, originally fresh ilmenite could be altered in the laboratory under simulated weathering conditions. ETCH FIGURES ON ILMENITE CRYSTAL PLANES A rather large ilmenite crystal from Sannidal, Norway, about 11/2 x 11/2 x 2 inches, bearing a few natural faces with imperfect surfaces, was prepared for etching by grinding and polishing crystal planes whose Bravais-Miller indices (referred to hexagonal axes) had been determined by X-ray methods to be 00.1, 10.0, and 10.2. Metallographic examination showed that the crystal was homogeneous, and uniform in orientation, except for the presence of twin lamellae which were oriented with their broadest surfaces parallel to the 00.1 planes of the large crystal. Since ilmenite is rather inert to standard etch reagents, strong solutions of hot sulfuric acid and cold hydrofluoric acid were used. Preliminary tests established that better etch figures were obtained using HF than with H2So4 solutions. The etch figures produced on 00.1 and 01.0 planes of ilmenite by etching 3 to 5 minutes in cold concentrated HF solution are shown in Figures 1, 2, and 3. The figures on 00.1 are "etch hillocks", rather than pits, and according to Honess(6), this phenomenon

Jan 1, 1961

By J. Douglas, W. T. Ford, G. E. Henderson, I. F. Roebuck

An implicit numerical method is presented for simulating the differential and algebraic relations governing one-dimensional three-phase flow in porous media. The method is based upon compositional representation of the hydrocarbon system. Variable physical properties and water-oil capillary forces are included in the formulation; however, the effect of gravity is ignored. The effects of changing composition and mass transfer are considered through the use of known phase behavior concepts and correlations. The differential and algebraic equations and the numeric approximations to these equations are presented, and the computing algorithm is discussed in detail. This method is used to simulate a single-phase, two hydrocarbon-component gas displacement laboratory experiment. The method is also applied in simulating a solution gas drive and a gas injection problem using a reservoir oil represented as a nine-component mixture. For comparative purposes, the latter two problems also are examined with a one-dimensional volumetric model. The results of the simulation of the laboratory experiment demonstrate the relevance of this method to problems wherein the effects of mass transfer significantly enter the displacement process. The comparison of the results of the volumetric and compositional methods applied to a reservoir system illustrates not only the advantages and benefits derived from examining the hydrocarbon system composition, but also the necessity of this approach in application to a general reservoir problem. INTRODUCTION The advent of large memory, high-speed digital computers enabled reservoir engineers and mathematicians to pool their knowledge in the development of sophisticated techniques for the prediction of oil and gas reservoir performance. These techniques have progressed from volumetric and trend analysis calculations to the development of "mathematical models", or reservoir simulators. These latter models, in general, utilize finite difference approximations to the rather complex partial differential equations that mathematically describe the physics and thermodynamics of fluid flow in porous media. Many forms of this "model" type are available today, and much of the recent literature has been devoted to either improving known methods or developing new techniques for use in these prediction methods. The more advanced techniques in use at this time consist of the solution of differential and algebraic relations that describe a reservoir geometry in one, two or three dimensions and that incorporate reservoir hydrocarbon and water fluid systems that obey specified physical laws, having properties that might easily be measured in the laboratory. This type of model utilizes formation volume factors and other volumetric data, expressed directly as functions of pressure, to represent the reservoir fluid system. Experience gained with this type of simulator, commonly known as a "beta-factor model", has indicated that the approach has many outstanding defects, most of them due to the assumption that those properties measured in the laboratory completely describe the behavior of the fluid system throughout the life of the project under study. The effects of varying hydrocarbon composition and mass transfer between the phases of this hydrocarbon system must, by definition, be ignored. These defects, which may often be ignored in

Jan 1, 1970

By P. van Meurs, C. van der Poel

From observations of the linear displacement of oil by water from a porous medium as visualized in transparent models, new insight into the mechanism of vircous fingering, as occurring in the case of unfavorable oil-water viscosity ratios, has been obtained. A mathematical description of the observed phenomena leads to formulas. expressing both oil production and pressure drop across the formation as a func-tion of cumulative water injection with the oil-water viscosity ratio as a parameter. Theoretical results are compared with those obtained from scaled laboratory experiments in sand-filled rubes in which water was injected at the one end and production was obtained from the other end. In the second part of the paper the same principles are used to calculate a more realistic production scheme in which production of oil and water is obtained from wells regularly distributed over the field. Results show satisfactory agreement with production data of an actual oil field. Application of the theory is extremely simple as only a limited number of readily available parameters have to be known and results are presented in an analytical form. INTRODUCTION When oil is displaced from a porous medium by water having a viscosity lower than the oil. the oil-water interface is essentially unstable and has the tendency to break up into what are called "fingers" or "streamers". This Viscous fingering has been noted before by Engelberts and Klinkenberg' and has recently been demonstrated in transparent models by van Meurs-. It is highly questionable whether it is at all possible to perform relative permeability measurements In the laboratory reproducing the same distribution of water and oil in the cores as Present in practice. Therefore, we based our descriptive theory on an idealized model which was "tailored" to the fingering mechanism as observed in transparent models under scaled conditions. This approach was the more attractive, as a highly simplified model had already yielded an adequate description of the process. Though we are aware our theoretical approach can be extended in several ways, it is given here in its simplest form to explain the fundamental principles and to avoid complicated formulas. SIMPLEST CASE OF A LINEAR DRIVE Observations and Theoretical Model In Fig. 1 a series of photographs shows how water penetrates into a formation saturated with an oil which is considerably more viscous than the water. The technique for obtaining these photographs, using transparent models, has already been described?. The fingers can be clearly observed and are seen to be mainly parallel to the flow direction. Note that the contours of these fingers are not smooth but show protrusions of water enclosing oil pockets.

By L. F. Elkins

In 1946 Littlefield, Gray and Godbold published a thorough geologic description of the West Edmond Hunton Lime reservoir, located in Central Oklahoma, and discussion of its early performance.l They stressed implications of the nature of the reservoir rock on performance to be expected with production both by natural means and by pressure maintenance with gas injection. The essence of their analysis was that presence of an extensive interconnected system of fractures and solution channels constituting about 10 percent of the total reservoir void space in otherwise very low permeability reservoir rock would result in severe channeling of naturally encroaching water or injected gas through the fractures with little or no benefit to ultimate recovery of oil. Considerations were in progress to unitize the reservoir for pressure maintenance by gas injection. Some other engineers involved in the project had considerably different opinions at the time regarding continuity of fractures and regarding effects of such fractures on performance. The field was unitized Oct. 1, 1947. Subsequently, pilot gas injection tests were conducted in four areas of the field; large scale waterflooding was conducted in two parts of the field in the Bois d'Arc or Upper Hunton; and pilot waterflooding was conducted in the Chimney Hill or Lower Hunton. Full scale water-flooding is now in progress in the Chimney Hill section where it is oil productive. The primary purpose of this paper is to divine some of the secrets of this reservoir through analysis of all aspects of its performance, with the hope that it will be of help to engineers and operators in developing and operating other low permeability fractured reservoirs. In part, this considers variance of actual performance from some widely used conventional reservoir engineering assumptions. For brevity, other aspects of unit operation and its benefits to the owners are not discussed. Oil and Gas in Place For unit participation, net pay in the Bois #Arc defined by the gross SP anomaly on electric logs totalled 2,100,000 acre-ft in the 30,000-acre reservoir. Development of West Edmond predated modem logging techniques that now would be used to determine "net" or "productive" pay. In Fig. 1 this log pick is compared with detailed core analysis and core descripfion by Littlefield et a1.l for one well identified as A on the index map in Fig. 2. Not all of the gross Bois #Arc section had oil staining, but all cores with measureable intergranular permeability were stained. (The committee selection of the Bois d'Arc SP interval also includes the Frisco section as identified from cores by Littlefield.) The Bois d'Arc is oil productive throughout the entire 600-ft monoclinal accumulation from the eastern upstructure erosional pinchout to the west-side oil-water contact. (See Littlefield et al.l and McGee and Jenkins for geologic detail.)

Jan 1, 1970

By C. S. Samis, D. P. Seraphim

In the presence of oxygen, galena is oxidized in an aqueous medium containing ammonium acetate in accordance with the following reaction: PbS + 1/2 0, + 2 NH~Ac -» PbAc, + So + 2 NH: + H2O. This reaction was studied in an autoclave under oxygen pressure by polarographic measurement of the concentration of lead ion in solution. Oxidation rate is parabolic. The reaction products are lead acetote and elemental sulfur, the latter forming a film on the galena crystal. Reaction kinetics appear to be determined by the diffusion of lead from the galena phase to the solution interface through the film of sulfur. RECENT applications of pressure leaching on a commercial scale have emphasized the metallurgical importance of the oxidation of sulfide minerals in aqueous media by oxygen under pressure. Research experience has shown that sulfates, thio-nates, and metal oxides are the products of humid oxidation of sulfide minerals in basic solutions. The metal oxides may be solubilized in the presence of a complexing ion.' Certain oxide minerals are also amenable to treatment by this type of oxidation.' Several kinetic studies have been made of reactions of this type. Stenhouse' found that pyrite (FeS?) oxidized in a sodium hydroxide solution to produce a soluble sulfate and an insoluble layer of iron oxide and that the kinetics of the reaction were controlled by oxygen diffusion through the oxide layer. Andersen' studied the oxidation of galena (PbS) in a sodium hydroxide solution under oxygen pressure. He found the lead and sulfur to be simultaneously oxidized to produce plum bite ions and sulfate ions, and that the reaction rate was proportional to the surface area of the galena. The kinetics appeared to be determined by the reaction of oxygen and water with the galena surface. In each of the above reactions the metal and sulfur components of the mineral were oxidized simultaneously and at equal rates. In the present study the alkaline medium has been replaced by an approximately neutral solution of ammonium acetate in which lead sulfate is soluble. Under the conditions no sulfate was produced and no sulfur could be detected in the solution. The metal component of the mineral was solubilized and the liberated sulfur remained as a film of elemental sulfur on the surface of the galena crystal. A kinetic study of this reaction was undertaken to establish the kinetic mechanism and to determine the rate controlling step. Specific reaction rates were determined by oxidizing crystals of galena of measured surface areas in an autoclave maintained at the desired temperature and oxygen pressure. The rate of reaction was followed by measuring the concentration of lead in the ammonium acetate solution with a cathode ray polarograph, employing stationary platinum electrodes incorporated in the autoclave. Experimental Materials and Equipment Crystals of galena obtained from Violamac Mines, Retallack, B. C., were selected for oxidation because

Jan 1, 1957