Search Documents

By Wilbur H. Somerton

In order to evaluate the quantitative determination of porosity by various neutron logging systems, four test wells were constructed to simulate some of the possible variab1e.r found in field logging conditions. All of the test holes contain sections of limestone with one section of sandstone included in one well. Each well has a representative section of high, medium, and low porosity. Three borehole sizes are represented by three of the test wells, all containing fresh water. One test well contains saturated brine. Statistical and tune constant checks are made for each system evaluated. Runs and repeat runs are made with the sonde centered and eccentric. Factors evaluated are the effect of borehole size, sonde position, salinity, chemical composition of rock, and statistical variation. A semi-quantitative method of evaluation of these effects is used giving the average maximum per cent error possible for each. Drift and calibration procedures are also discussed. The effect of the various factors appears to be largely a function of the phenomena detected as well as the means of detection. Graphic illustrations of some of these effects are presented. It is concluded, that for accurate log interpretation, the hole size and sonde position should be known. INTRODUCTION The use of neutron logs for the derivation of porosity has become an accepted practice in many areas. In recent years, many different neutron logging systems have been offered to the industry. The neutron sources presently available are of the radium-beryllium or polonium -beryllium type. Detector types used are ionization chambers, Geiger-Mueller-type counters, proportional counters, or scintillation detectors used to detect or measure gamma rays, thermal neutrons, epithermal neutrons, fast neutrons or combinations of these. It is well-known that all systems do not respond alike in all environments. As a result, testing facilities and an evaluation method were developed for use in evaluating the porosity discriminating ability of neutron logging systems using conditions simulating those normally encountered in the field. With such an evaluation, the suitability of a logging system can be determined lor selected field conditions. TEST FACILITIES The testing facilities consist of wells containing rocks of various porosities, constructed to include as many parameters as practicable. Three pits were drilled to a depth of 20 ft and the fourth to a depth of 9 ft, all having 48-in. diameters. Rat holes were dug to a depth of 6 ft in the center of each pit and fitted with sections of pipe having appropriate diameters. The top of this pipe was located approximately 8 in. above the pit bottom. A 6-in, layer of cement was then poured on the bottom and while still wet, 42-in. OD corrugated, galvanized culvert pipe was set on bottom and centered in place. This pipe extends approximately 1 ft above the surface of the ground. The annulus around the culvert pipe was then filled with cement to the top. A 1- or 2-in. layer of sand was placed on the cement bottom to provide an adjustable bed on which to place the series of rocks set in each hole. A cross section of a typical test well containing the porous rocks is shown in Fig. 1.

By J. W. Marx, R. H. Langenheim

The injection of heat-bearing fluids may offer a wider application to secondary and tertiary recovery from conventional oil reservoirs than inderground combustion since the process is more easily controlled and the reservoir requirements, in general, are less critical. In all cares, conduction heat losses to overburden and un-derburden will impose an economic limit upon the size of the area which can be swept out from any one injection point for any given set of conditions. The present report describes a method for c.rtin~ating thermal invasion rates, cumulative heated area, and theoretical economic limits for sustained hot-fluid injection at a constant rate into an idealized reservoir. Full allow-ance is made for non-productive reservoir heat losses. Required solutions can be obtained directly from the numerical tables given here, as soon as all operating conditions are specified, withou extensive mathematica1 ~nanipulation on the part of the user. The injection of a heat-bearing fluid into a reservoir is frequently proposed as a means of secondary or tertiary oil production. This form of thermal recovery may offer a wider application to conventional oil reservoirs than underground combustion since the process is more easily controlled and the reservoir requirements, in general, are less critical. Because of its large gross heat capacity, steam appears to be the most efficient heat injection medium, although mixtures of steam and other gases, hot water, hot oil and hot non-condensable gases also have been employed. Whatever the nature of the injection medium, economic evaluation of such thermal recovery processes will depend upon thermal invasion rates, or rates of productive reservoir heating, at any given time after the start of heat injection. In all cases, conduction heat losses to overburden and underburden will impose an economic limit upon the size of the area which can be swept out from any one injection point, for a given set of reservoir conditions, at any given heat injection rate. The present report describes a method for estimating thermal invasion rates, cumulative heated area and theoretical economic limits for sustained heat injection at a constant rate into an idealized reservoir, making full allowance for non-productive reservoir heat losses. Required solutions can be obtained directly from the nu- merical tables given here, as soon as all operating conditions are specified, without extensive mathematical manipulation. It should be noted that ordinary water flooding may be used to move an annular-heated zone outward, even though the economic areal limit for straight heat injection has been reached. Thus well spacing need not be confined to the areal limits as determined here. However, the present calculations do not cover such supplementary heat scavenging programs. The operating element consists of a radial flow system, concentric about the point of injection, with the idealized step function temperature profile illustrated by the dashed line in Fig. 1. The solid curve shown in Fig. 1 gives a qualitative picture of the true radial temperature distribution. It will be shown later that the expressions given below, which are derived from the simplifying step idealization (dashed line in Fig. I), also can be interpreted in terms of the more realistic temperature distribution shown in Fig. 2. For this case, the thermal invasion radius is simply redefined as the distance from injection well to the mid-point of the temperature distribution.

By Harvey T. Kennedy, Erasmo T. Guerrero

In this investigation, the effects of the surface tension of brine and the interfacial tension between oil and brine on the recovery of oil by water flooding of Woodbine sand were studied. Variation of the surface forces was obtained by the addition of small amounts of commercially available chemicals, which did not appreciably change the surface tension of the oil. Fourteen surface-active agents, which were found to be effective in lowering the surface and interfacial tensions, were further investigated to determine the extent of their adsorption by the Woodbine Sand cores. It was found that adsorption was sufficient in all cases to prevent any substantial concentration of chemical at the flood front. The flow system was, therefore, brought to equilibrium with definite concentrations in the flood water by circulation before each test. The data thus have no relation to the results obtainable by using the chemicals in full-scale water floods, but should be helpful in predicting the performance of surface-active chemicals which are not lost by adsorption, if such can be found. Flow tests were made on various concentrations 01 six chemicals selected from the 14 investigated for surface activity and adsorption, with and without the presence of gas in the core. It was found that lowering the interfacial tension had a tendency toward reducing the recovery of oil and that the reduction of the surface tension of the brine had little or no effect. A substantial beneficial effect, however, resulted from the presence of gas in the cores. INTRODUCTION The movement and distribution of crude oil in a porous medium are governed by such factors as surfacc and interfacial tensions, viscosity of the oil, pore size configuration, and wetting characteristics of the medium. The knowledge of the part played by each of these factors is limited. If the oil is assumed to be distributed in droplets whose movement is restricted only by the pore openings, the following equation may be applied: ?P= 2r(1/r1-1/r2) .......(1) In this equation ?P is the pressure drop necessary to force an oil droplet through the construction, r is the oil-brine interfacial tension, r1 is the radius of curvature at the forward and smaller end of the deformed droplet, and r2 is the larger radius. This equation neglects the effects of viscosity, as a static or near static condition is assumed. It would apply equally well for a bubble of gas in liquid if r were taken as the surface tension of the liquid. This relation shows that droplets can be more easily displaced if the interfacial tension is decreased. The displacement of oil from the surface of a sand particle may also be considered from the basis of selective wetting assuming that the oil droplets are smaller than the pore constructions and that the sand surfaces may be considered as planes. A relation of the forces involved can be derived from Dupre's equation to give: Wos = row (1 +Cos?)........(2) where Wos = Work of adhesion between the sand and the oil row = Oil-water interfacial tension 0 = Contact angle ros = Oil-solid interfacial tension This equation would indicate that the displacement of oil from a sand surface is more readily obtained when the water-oil interfacial tension is reduced provided the contact angle remains constant. Evidently. however, the values of Cos 0 may vary along with the water-oil interfacial tension due to changes in the other surfaces of the system. The above theoretical considerations suggest that perhaps even the static behavior of a dispersed oil phase in a porous medium is considerably more complex than may be inferred by the simple systems

Jan 1, 1955

By R. L. Slobod, B. H. Caudle, R. A. Erickson

In predicting the sweepout pattern efficiency to be expected in a .secondary recovery operation, the reservoir engineer is often confronted with a situation in which part of the producing formation lies between the last row of wells and the reservoir boundary. To date, there has been little or no reference in the literalure to sweepout pattern efficiencies in such areas. This paper presents the results, of a laboratory study in which the area contacted around the edge of a reservoir injected in a five-spot pattern was determined. The X-ray shadowgraph technique was used in this study. It was found that surprisingly large portions of the reservoir area lying outside the well network are conracted by the injected fluid before, abndonment conditions are reached. In fact, it can he said that at least 00 per cent of the circa lying outside the last row of wells and within a distance of one well spacing of these wells will ultimately be contacted by the injected fluid in the case of most water floods. Frorn the results given in this paper, it is possible to estitilate not only the extent to which encroachment into these areas will occur but also the injection volume necessary to produce the wells to abnndonment conditions. INTRODUCTION While considerable work has hern done in recent years1.?.'.''.' evaluating the fraction of the reservoir contacted as a function of the mobility ratio for several different well patterns, there has been practically no mention of what happens to the oil in that portion of the reservoir which lies between the last row of wells and the reservoir boundary. In fact, the assumption is generally made that little or no recovery can he expected From this area. This assumption is incorrect and in many cases will lead to seriouq errors in calculating reservoir performance. The fraction of the area outside the regular well pattern contacted in a five-spot injection system has been studied using the X-ray technique5 to follow the progress of actual floods in porous models of reservoir elemental areas. In this paper, the method for carrying out such a study and the results showing the fraction of the area in question swept as a function of mobility ratio and the produced fluids are presented. THE REPRESENTATION OF A RESERVOIR BY ITS ELEMENTAL AREAS She portion of the reservoir lying between the edge of a reservoir and a line connecting the outermost wells may be roughly designated as the part of the reservoir outside the regular well pattern. An example of such an area is shown as the shaded portion in Fig. 1. The case of a square reservoir utilizing a five-spot injection program and with one full well spacing between the outermost wells and the boundary (Fig. 1) was selected as the single example which would provide the most insight into

Jan 1, 1956

By R. P. Alger, C. A. Doh, M. P. Tixier

The principle, the equipment and field operation of sonic logging are described. The tfio-receiver system produces logs independent of hole size and mud. Field experience is given and forms the basis for the interpretation of the log. The derivaiion of porosity values from measured velocities is discussed, according to the type of formations (hard formations, compacted sands, and unconsolidated formations). The time-average equation proposed by M. R. J. Wyllie, A. R. Gregory and L. W. Gardner is uscd as a basis for the computation of porosity in limestones, cemented sandstones and compacted sands. Variations in the lithologic character of limestones do not seem to change the porosity calibration markedly. The compaction of sands is related to the compaction of shale adjacent to them; and, thus, the shale velocity is used for the establishtnent of empirical relations for the computation of porosity in un-consolidated formations. Various forrnu.las are tentatively presented to account for shale and fluid content. Field experience demonstrates that considerable attenuation of sonic energy takes place in unconsolidated formations, particularly when gas bearing, and in fractrured formations. Unusually large attenuation produces skipped cycles, a feature easily recognized. The application of sonic logging to structura1 studies is featured, showing its possible integration with the dipmeter. Correlation and its application to the irrterpretation of seisrmic surveys are reviewed. The paper is illustrated with field examples. Sonic logging is the recording of the time required for a sound wave to traverse a definite length of formation. Sonic travel-times are inversely proportional to the speed of sound in the various formations. The speed of sound in subsurface formations depends upon the elastic properties of the rock matrix, the porosity of the formations and their fluid content and pressure. Below the "weathered" or low-velocity layer extending from 50 to 100 ft or so below the surface, sound velocities may range from about 6,000 ft/sec in shallow shales to as much as 34,000 ft/sec in dolomites. In hard formations (well cemented and/or compacted), the sonic log reflects the amount of fluid in the formations; hence, it correlates well with their porosity. In unconsolidated formations, which are usually of fairly high porosity, the sonic log gives an approach to porosity determination, when its readings are corrected for lack of compaction, shaliness, and fluid content. In such formations the sonic log may also indicate the presence of gas and may distinguish between oil- and water-bearing beds. Field experience with the sonic log in soft formations is more limited than that in hard formations. As a result, interpretation in soft formations is not as well developed at the present time. In all types of formations the sonic log is of considerable value for correlation and for the more detailed and accurate interpretation of seismic data. EQUIPMENT AND PRINCIPLE FIELD OPERATION The basic features of the sonic log" are a two-receiver system and two available short spacings. The log is recorded on film using a standard field unit. An SP curve is run along with the sonic curve to give a more interpretable log and to establish absclute depth control for comparison with other logs run in the same well. If desired, a gamma-ray curve may be recorded simultaneously with the sonic curve. Fig. 1 shows the sonic-gamma ray combination sonde. The sonic portion consists of a sound transmitter and three receivers mounted 3, 4 and 6 ft from the transmitter. When the first and second receivers are used, a sonic log of I -ft spacing (distance between receivers) is obtainea. Use of the first and third receivers provides a 3-ft spacing. The transmitter emits pulses at the rate of l0/sec. After the transmission of a pulse the first arrival of sound energy at each receiver triggers its response system. The differ-

By J. Wade Watkins, E. S. Mardock

The accurate evaluation of reservoir-performance characteristics in the secondary recovery of petroleum by water flooding requires use of a water tracer that may be injected into water-input wells and detected at oil-production wells to supplement data obtained from core analyses, wellhead tests, and subsurface measurements. Radioactive iodine has been used successfully as a water tracer in field tests to determine: (I) relative rates and patterns of flow of injected water between water-input and oil-production wells and (2) zones of excessive water entry into oil-production wells. Laboratory evaluations of potential water tracers, previous tracer studies, the value of using a radioactive tracer, general field procedures, and the use of surface and subsurface instruments for the detection of the emitted gamma radiation, are summarized. Data from the field tests are presented graphically and discussed in detail. It is concluded that the radioactive-tracer method, using radioactive iodine, may be used successfully to measure either the relative rates and patterns of flow or zones of excessive water entry into wells under conditions of comparatively rapid transit time between wells. INTRODUCTION Extensive use is made of data obtained from core analyses, wellhead tests, fluid characteristics, and subsurface measurements in evaluating the sweep efficiency of water injected into oil sands for the recovery of oil. Theoretical flow rates and patterns may be calculated from those data using radial-flow formulas, if it is assumed that the physical conditions in the productive formation are homogeneous. Unfortunately, homogeneous conditions rarely, if ever, exist in oil-productive formations. The use of a tracer that may be injected into an oil sand and detected quantitatively, or even qualitatively, at offsetting oil-production wells provides basic data that may be used in determining more accurately the subsurface rates and patterns of flow of injected water between wells than is possible by theoretical calculations based on assumed conditions. Consideration of the data obtained by using a water tracer assists in the application of remedial measures to water-input wells, such as plugging of channels, or selective plugging of highly permeable zones, thereby effecting a more uniform flood and a greater ultimate oil recovery. A water tracer should have the following characteristics: (1) low adsorption on solid reservoir material; (2) high solubility in water; (3) low or negligible solubility in crude petroleum; (4) a wide range of solubility of compounds that may be formed by chemical reactions with ions present in reservoir rocks or waters; (5) high detectability in low concentrations by portable apparatus; and (6) general availability at low cost. Additional desirable characteristics are that the tracer should be nonhazardous in nature under normal working conditions and quantitatively detectable in the well bore, as well as at the surface. Review of Previous Tracer Studies The need for a water tracer has been recognized for many years. As far back as 1906, Dole1 described the use of fluorescein in flow-rate studies made in France in 1901. In 1921, Ambrose" discussed in detail the use of organic and inorganic dyes, chlorides, nitrates, and other anions, lithium salts, and the Slichter electrical

Jan 1, 1955

By Fred H. Poettmann, A. L. Benham

correlation of waters and identification of mixtures of two and three waters is important in secondary recovery evaluation studies and studies involving underground movement of waters. The importance of such studies suggests the need for a method to correlate a great number of water analyses to accurately identify and improve the understanding of mixtures of waters. Geochemical analyses of waters furnish a quantitative measure of each of the ions present in a water. Table L shows a portion of the data on many oilfield water analysis report forms in use today. Under "parts per million" are listed the quantitative amounts of each ion in milligrams per liter. Following the work of Stabler' and Palmer", some laboratories have reported additional data to indicate the reaction capacity, or reacting value of each ion. Table I lists data under two subheads for the "reacting value". Under the first subhead, the reacting value is reported in gram equivalents per million, hereinafter called equivalents per million. This value is determined by multiplying milligrams per liter of the ion by a factor equal to the reciprocal of the equivalent weight of the ion. Under the second subhead (per cent) is listed the percentage of each ion based on the total equivalents per million. Sodium is not usually an analyzed value in routine analyses, but is determined by the difference between total negative equivalent per million values and the sum of all of positive ion values for which analyses were made. This difference is reported as sodium, hereafter labelled "sodium". On this basis, the subtotals of the equivalent per million and per cent reacting values of the positive and negative ions are equal. CORRELATION OF WATER ANALYSES DATA Previously reported methods for correlating oilfield waters have used part of the data in Table 1, and can be classified into two categories. Correlation by Chemical Character Correlation methods using the chemical character of waters, as represented by the per cent reacting values of the ionic constituents, have been used by Palmer2 and Tickell. These methods have the limitation of neglecting the ionic concentration of the waters. This limitation can be serious in areas where waters have the same chemical character, but differ in concentration. Correlation by Concentration of the Ionic Constituents The Reistle' and the Stiff methods use the ionic concentration in parts per million and equivalents per million. These methods are effective in differentiating waters of similar chemical character and varying concentration, as in the highly saline waters of the Mid-Continent area. These methods have the disadvantage that they correlate waters by comparing many graphs of individual water analyses to determine the similarity to, or the variance from "type" waters. PROPOSED METHOD OF CORRELATING WATER ANALYSES DATA Correlation graphs presented in Fig. 1 represent water analyses data as points to facilitate the correlation of waters. The graph on the left of Fig. 1 is a plot of per cent reacting value of magnesium vs calcium. The right

By P. Hazebroek, C. S. Mathews, E. E. Allen, H. C. Lefkovits

A rigorous study was made of the behavior of reservoirs composed of horizontal layers, unconnected except at the well and filled with a compressible fluid. The report is presented in two parts. Part I deals with the practical implications of the results obtained; Part 11 contains the mathematical derivation of the solution. The main portion of Part I is concerned with application of the results to reservoirs composed of two layers. Relative rates of depletion of the layers are studied, and it is shown that differential depletion between layers exists in the transient stage of the reservoir, which generally is of an order of magnitude longer than the transient stage in a single-layer reservoir. In the transient stage, the more permeable layer is depleted faster than the less permeable. However, the reservoir approaches a steady-state situation where both layers contribute equally to production. Theoretical build-up curves from such reservoirs are presented, and the influence of a skin effect and flow into casing and tubing are also evaluated. The build-up curves are predicted to have, first, the familiar logarithmic straight-line section and, subsequently, a rising and flattened section. The results found are compared with an extension of an approximate theory, and it is found that this simplified theory is applicable over most of the life of the reservoir. Extensions to reservoirs containing more than two layers are indicated, and examples of the performance and build-up curves for three-layered reservoirs are presented. The theoretical results on build-up curves were applied to field examples of such reservoirs, and satisfactory results were obtained. From the build-up curve, the pemeability-thickness product, the wellbore damage and the static pressure are obtained. However, it apparently is not possible to determine the properties of the individual layers from the combined build-up curve. A large number of reservoirs have been found where the producing formation is composed of two or more layers of differing physical characteristics, such as permeability and porosity, and where also the thicknesses of the layers differ. Since in general a differential depletion should exist between the layers of such reservoirs, it would be advantageous to know the amount of the differential depletion under a given production history. It also would be valuable to know the pressure behavior, because it could be undesirable to produce a reservoir in such a manner that the pressure in one of the layers would have fallen below the bubble point of the oil at a time when another layer might still be virtually undepleted. The determination of static pressures in shut-in wells with slow build-up is another feature of importance in layered reservoirs. If the behavior of the reservoir differs appreciably from the predicted behavior of a single-layer reservoir, the type of analysis for build-up curves suggested for single-layer reservoirs is not applicable. Several authors have treated the problem of finding the pressure and production characteristics of multilayer reservoirs. Horner treated the problem of n layers in an infinite field, the layers being connected only at the well. Tempelaar-Lietz2 has presented an approximate treatment of the characteristics of a bounded reservoir composed of two layers of different permeabilities and equal thicknesses. Homer gave no information on differential depletion between layers, and the treatment of Tempelaar-Lietz was not applicable for prediction of pressure build-up behavior. In addition, it was thought that a more rigorous treatment than the one by Tempelaar-Lietz would be valuable in determining the range of validity of this solution. The present study, therefore, was made of a bounded reservoir composed of two or more layers,

By D. R. Parrish, V. J. Berry

Reverse combustion is one thermal method of recovering hydrocarbons from porous undergrortnd formations containing oil or tar. In applying this method, air is introduced via an injection well and the mixture of air and hydrocarbons is ignited in the production well. A combustion zone then recedes toward the injection well, counter-current to the air flow. If the combustion-zorie temperature is sufficiently high, the oil or tar irz place is distilled arid cracked. The hydrocnrhon flows as a vapor to the production well and subsequently is condensed at the surface. Maximum temperature and velocity of morcment arc the two dependent variubles defining the progress of the combustion zone. A theoretical analysis has been made of heat flow in the reverse-comhustion process as.srtrning linear flow in a homogeneous system. The differentirll equatious, which include the oxygen-hydrocar-hon rctrction rate, have been solved numerically. Rescllts indicate that the maximum temperature reached and he cornbustion-zone relocity both increase with an incrense in air-injection rote. Heat loss to surroundings has little effect on the maximum comhustion-zone ten?perature achieved, but it is reflected in a reduced com-bustion-zorzc, velocity. It is also predicted that an increase in the oxygen-hydrocarbon reaction rate results in a reduction in the maximum temperature reached. The calculated results are in agreement with results from reverse-combustion experinments using samples of a tar sand. INTRODUCTION Reverse combustion is one thermal method of recovering hydrocarbons from porous underground forma-tions containing oil or tar.' In applying this method, air is introduced into the underground formation via an in-jection well. In one or more adjacent production wells, the mixture of air and hydrocarbons is ignited at the sand face. The combustion zone thus formed recedes toward the injection well, counter-current to the air flow (Fig. 1). Under certain conditions, the rate of recession of the combustion zone is sufficiently rapid that only a small fraction of the hydrocarbon at any point in the formation is consumed by combustion. If the combustion-zone temperature is sufficiently high, the uncon-sumed oil or tar in place may be partially distilled or cracked. This material then can flow as a vapor with the gaseous combustion products to the production well and subsequently be condensed at the surface. Under other conditions, the combustion-zone temperature may be so low that little, if any, of the hydrocarbon is vaporized. Nevertheless, the viscosity of the hydrocarbon may be reduced enough so that other recovery methods subsequently can be applied to produce oil at attractive rates. The maximum temperature achieved in the combustion zone and the velocity of movement are the two dependent variables defining the progress of the combustion zone. Presumably. knowledge of these factors will permit estimates to be made of oil recovery efficiency, oil and gas producticn rates and average amount of injection air required per barrel of oil recovered. The maximum temperature and combustion-zone vclocity are, in turn, expected to be dependent upon air flow rate, heat losses to surrounding formations, physical properties of the rock and the hydrocarbon. the extent to which the porous formation is saturated with the hydrocarbon, the rate of oxidation of the hydrocarbon and the amount of heat released. It is desirable that the relative importance of each of these factors be understood so that the results of laboratory tests of oil recovery by reverse combustion can be properly interpreted in terms of expected field performance. It is particularly desirable to assess the importance of heat lcsses because these are difficult to control or eliminate

By C. R. Johnson, R. A. Greenkorn, G. W. Widner

This paper describes a short, simple well-testing procedure evolved from unsteady-state concepts to appraise wellbore damage in water-input wells. Input rates are changed at specific intervals over a 30-minute period and the wellhead pressure recorded at the end of each interval. From these data, the pressure drop due to the "skin" can be computed. This pressure drop in terms of psi/B/D of input rate is referred to as the "skin factor". Two of these variable-rate tests are required to determine change of wellbore damage and to evaluate well-treatment techniques. There are several advantages to the variable-rate procedure. The necessary data can be obtained with wellhead gauges and meters. There is little interference with normal operations. The short testing period eliminates many interpretative difficulties that come from interference from adjacent wells. The shut-in period required for standard tests before the plot of pressure vs log of time attains a straight line is often 12 hours or more. Such a time period is incronvenient, especially when a well must be tested repeatedly to evaluate change in wellbore characteristics. This variable-rate procedure has been successfully npplied to water-input wells in three fields. For some of the wells, skin factors were calculated from both the standard and the variable-rate tests. These results show good agreement. THEORY OF THE VARIABLE-RATE PROCEDURE In many waterflood input wells, the pressure in the wellbore differs from the pressure in the reservoir immediately adjacent to the well because of the restriction of flow through the skin of the well. If the true pressure of the reservoir at the well, p,, is that pressure which would exist in the well if no skin were present, and p,' is the actual pressure as measured in the well, this difference, is the skin effect. Fig. 1 shows the change in skin effect with change in flow rate. As the flow is reduced, the skin effect is reduced by an amount equal to the change in wellbore pressure minus the change in reservoir pressure. The equation of pressure change without skin using the point-source solution of the diffusivity equation is (after van Everdingen')

By G. W. Tracy, R. D. Carter

A method based on an approach by Hurst' is developed for calculating water influx behavior. Using this method, superposition calculations may be eliminated. The principal diflerence between this method and that of Hurst is that, over finite intervals of time, constant oil-productton rates are assumed by Hurst whereas constant water influx rates are assumed in the presently described method. By adopting the assumption of constant water influx rate for finite time periods, a method is developed for calculating water influx behavior which lends itself to convenient combination with the material balance equation in the Schilthuis form. Above the bubble point, an equation is given for the explicit step-wise calculation of pressure history for a prescribed oil-production history. Use of this method of calculating water influx behavior should result in a considerable saving of time as compared to superposition calculations when the calculations must be made without using a digital computer. The connection is shown between the method described herein and the method developed by van Ever-dingen and Hursr2 (which requires superposition cal-culations). Also, results obtained using the two methods are compared.

By J. E. Smith, P. B. Crawford, B. L. Laudrum

A method is presented by which the displacement efficiency of a steam-injection process may be determined. The method consists of a determination of the position of the steam front at a particular time and the application of the Buckley-Leverett theory to the hot and cold regions of the reservoir with the appropriate hot or cold fluid properties. A uniform, linear flow system is assumed for this study. The method is applied to a number of reservoir conditions. The primary properties which affect the results of a steam-injection process are the oil viscosities at steam and reservoir temperatures, injection rate, pay thickness and oil saturation. The efficiency of steam injection is compared with that of water injection for a 10-acre, linearly flooded reservoir. In reservoirs where the oil viscosity is about 35 times the water viscosity, the steam process yields recoveries about 20 per cent higher than the water flood. Where the oil viscosity is near 500 cp, the steam process may recover from three to four times as much oil as a water flood. INTRODUCTION In a previous paper,' a method was suggested for evaluating a steam-injection process. Emphasis was placed upon the determination of the steam-front position in the reservoir at various times. In the present paper, the displacement process is studied in more detail. The approach is essentially that employed by Jordan, Rayne and Marshall3 in their study of hot-water injection. The steam-injection process is quite complex, and an exact description of the process has not been found; however, a qualitative description can be made. When a quantity of steam is injected into a relatively cold reservoir, the leading edge of the steam front continaally flows into cool sand, and some steam condenses. As more and more steam condenses, a bank of condensed water will form which is augmented by displaced connate water. The reservoir may be divided in three regions—(1) a region containing displaced fluids at reservoir temperatures, (2) a region contain- ing water at reservoir temperature and (3) a region occupied by steam. The reasons for the lack of a region of hot water or oil of any appreciable size have been discussed.' A method has been presented by which the location of the steam front can be determined even when heat losses to the regions surrounding the reservoir are sizable compared with the amount of heat injected. It remains to show the effectiveness of steam injection in increasing the ultimate oil recovery. To illustrate the technique for determining the increase in displacement efficiency and, hence, the increase in ultimate production, two sets of hypothetical reservoir conditions will be studied in detail.

By B. K. Larkin, H. R. Bailey

A model Of heat flow in an underground combustion process is studied. This model includes convection effects and thus is more general than previous studies which considered conduction as the only mechanistn for heat transfer. Both linear (tube run) and radial (field application) geometries are considered. The effects of ignition heaters, vertical losses and finite source width are considered for the linear case. The results are in the form of equations and are presented in graphical form for a number of cases. Convection effects increase frontal temperatures about 25 per cent over those computed for conductive transfer for typical field operating conditions. This increase in temperature is a result of heat being transported toward the front by the injection gas. Even greater temperature increases are realized as the per cent oxygen in the injection gas decreases. It is well known that compression costs are of considerable importance in estimating the economic feasibility of underground combustion. By assuming an ignition temperature for the combustion fuel, predictions of limiting conditions on fuel density and injection rate necessary to sustain the combustion zone are made. For typical field conditions, at least 0.75 lb/cu ft of fuel are needed with air as the injection gas. If the injection gas is 10 per cent oxygen and 90 per cent nitrogen, this figure is 0.69 lb/cu ft. INTRODUCTION The possibility of increasing oil recovery by underground combustion has been considered for many years. Recent field tests indicate that the underground combustion process is technically feasible. The economic feasibility depends to a large extent on the amount of air which must be injected to sustain combustion. Prediction of the success of employing underground combustion in a particular reservoir must be based on existing field tests, laboratory tube runs"' and solutions of mathematical or analog models. Vogel and Krueger5 devised an electrical analog of the heat transfer problem in an underground combus-ion process. They considered the problem of heat conduction from a cylindrical source with increasing radius assuming no vertical losses. Bailey and Larkin and Ramey solved the corresponding problem including vertical losses by considering a mathematical model of the conduction process. Ref. 9, 10 and 11 also present information related to the subject study. The present paper generalizes these results to include both conduction and convection mechanisms for heat flow. A model of the conduction-convection process is described. The partial differential equations governing this model are written and solved for a number of cases of interest. The formulas are evaluated and the results are presented graphically. A CONDUCTION-CONVECTION MODEL Description of the Model A number of models are employed in this paper, each of which approximate a particular physical situation depending on the geometry, i.e., linear or radial flow, the type of vertical heat losses, the nature of the source, etc. The geometry for both the linear and radial case is shown in Fig. 1. A particular linear model is described in this section, and variations from this model are described in later sections of the paper. Consider a homogeneous, semi-infinite porous prism with lateral sides insulated as shown in Fig. 1. The left end, x = 0, is heated for a time t, with a plane heater supplying a constant heat flux and during this time; gas (nitrogen) is passed through the pores of the prism at a constant mass flux. At time t, the heater is turned off and air is injected in place of the nitrogen. This causes combustion to take place, and a uniform source of width 1 is formed which moves at a constant velocity v, to the right. The frontal velocity v, is related to the gas velocity v,, through stoichiometric considerations

By W. L. Penberthy, H. J. Ramey

Experimental work on the combustion oil recovery process has consisted of both laboratory and field studies. Although field experiments are the ultimate test of any oil recovery process, they are costly, time consuming and difficult to analyze quantitatively. Laboratory combustion tube experiments can be opertzted far more rapidly and cheaply, but are subject to scaling and interpretation problems. This paper points out some important design problems, operational criteria and considerations in2portunt to interpretation of results. An analytical hcut model of movement of a burning front axiully along a cylinder with heat loss through an annular insulation was developed. The result was used to identify steady-state temperatrrre distributions both ahead of and behind the burning front, with and without heat loss. Results indicate potential operating limitations on the minimum burning front rielocity (or air flux) which may be used for any given combustion tube. Results also enable estimating the effective thermal diffusivity and over-all heat loss from experimental data and thickness of the brnning zone. Resu lts of operation of a combustion tube constructed recently verify this preliminary theory in the region immediately ahead of and behind the burning front surprisinRly well. INTRODUCTION Many field and laboratory studies of the forward combustion oil recovery process have been conducted since the early publications of Kuhn and Kochl and Grant and Szasz2 in 1953 and 1954. In view of the complex and costly nature of this type of investigation, it is not surprising that no complete theory of the nature of the forward combustion Drocess is yet available. However, gross effects are understood well enough that reasonable design procedures are available for planning field operations. Nelson and McNeil 3,4 have published two comprehensive papers concerning design procedures. One major consideration in planning field operations is the fuel concentration at the burning front. Fuel concentration controls air requirements — an important cost factor in forward combustion. Although fuel concentration can be estimated from field test results by various methods, results are subject to great uncertainty in view of natural limitations on experimental observations. Nelson and hlcNeil4 recommend that fuel concentration be determined from laboratory combustion tube studies. Fuel concentration is only one of many important factors which can be studied by combustion tube experimentation. An obvious goal of imporrance must be development of a comprehensive theory of the forward combustion process. If a theory of this process can be established which matches controlled laboratory experimentation, it should be possible to apply this theory to field operating conditions with some confidence. Laboratory combustion tube studies have already yielded important information concerning the combustion process. However, details concerning the design, construction and operation of combustion tubes are rare. Combustion tubes used by various investigators vary in size, length and mode of operation. Therefore, one purpose of this paper is to present considerations important to design, construction and operation of a combustion tube. In regard to previous combustion tube studies, attention is called to Refs. 1 through 9. These references describe a wide variety of equipment types and present a great deal of pertinent experimental data. In general, combustion tubes usually consist of thin-walled stainless steel tubes containing an oil sand pack mounted within a pressure jacket. Provisions have often been made to heat the tube ex tern a lly by separately controlled heaters to reduce heat losses. This step usually permits operation at low air fluxes (air rate per square foot burning front surface) similar to those encountered in field operations. Burning is usually conducted from the air inlet end of the tube to the outflow end. The tube orientation used has been vertical or horizontal. For vertical tubes, burning has been conducted vertically downwards. As a result of laboratory studies, many of the

Jan 1, 1967

By A. M. Skov, L. F. Elkins

Inclusion of anisotropic permeability in mathematical analysis of pressure transients observed during development of the huge Spraberry field indicates a major fracture trend which is in good agreement with that observed by fluid-injection tests spread over a 12- by 17-mile area. Delineation of this trend is important in selecting a pattern of injection for the pending large-scale water flooding in this field. Determination of reservoir parameters yielding best agreement between calculated pressures and observed reservoir pressures in newly completed wells was made using an IBM 650 computer. INTRODUCTION The Spraberry field covering 400,000 acres is a tight sand of less than 1-md permeability cut by an extensive system of vertical fractures. Primary recovery dominated by capillary retention of oil in the fractured sand matrix blocks is less than 10 per cent of oil in place. Strong forces of capillary imbibition of water into the sand, coupled with water flow under dynamic pressure gradient, indicate considerable increase in oil recovery can be achieved through water flooding. Best results will occur if the pattern of water injection is selected to force the water flow across the grain of the major fracture system. Existence of an oriented vertical fracture system in the Spraberry, observed first in cores, was highlighted more recently by the 144-fold contrast in permeability along and at right angles to the major fracture trend required to match relative water breakthrough times in Humble Oil & Refining Co.'s waterflood test there. Spraberry operators since have conducted two gas-injection tracer tests for further areal confirmation of the fracture trend. Re-analysis of early reservoir pressure transients for evidence of anisotropic permeability has permitted many more local determinations of major fracture trend without resort to further field tests. This paper is limited to updating analysis of reservoir pressure transients to include anisotropic permeability as a test for orientation of the major fracture trend in the Spraberry. The reader is referred to Refs. 1 and 2 for information about general Spraberry reservoir performance and to Refs. 3 and 4 for information about significance of fracture orientation in selection of the injection-well pattern for water flooding the Spraberry. RESERVOIR PRESSURE DATA—DRIVER AREA During early development of the Spraberry Driver area, Sohio Petroleum Co. made the extra effort to measure the initial pressure in each of the 71 wells in a 5-mile-long area immediately after completion. Progressively greater reductions in pressure ranging up to 400 psi were observed throughout the six-month development period. Detailed data are presented in Ref. 1. Since the reservoir oil was undersaturated some 300 psi initially, early reservoir performance involving 55 new well pressures is subject to analysis as flow of a single compressible fluid in a porous media. Assumption of uniform permeability in all directions yielded good agreement between calculated pressures and observed pressures of these wells in the earlier study,' but subsequent, additional, mathematical development to include anisotropic permeability in the transient pressure considerations and present availability of electronic computers to perform the much more extensive arithmetical calculations now yield even better agreement. The previous analysis, assuming uniform permeability, consisted essentially of calculating pressure reduction expanding circularly around each producing well and summing these effects at the time and location of each newly completed well for comparison with the measured pressure reduction. Permeability, effective fluid and rock compressibility, and permeability X thickness were varied until the best match with measured pressures was obtained. The present analysis, assuming anisotropic

By C. W. Arnold, H. L. Stone, D. L. Luffel

The purpose of this research is to determine (I) the efficiency of small banks of enriched gar driven by methane in displacing oil from a porous medium and (2) the effects of variation in bank size and composition of that efficiency. Most of the experiments were conducted in a sand-packed tube 20-ft long and 1/2-in. in diameter. The hydrocarbon system generally used was methane, butane and decane at 2,500 psia and 160°F. The results of these experiments indicate that, in the regions contacted by the gas, a small bank of an oil-miscible gas driven by methane can displace all of the oil in a piston-like manner. If the enriched gas is of such composition as to remain immiscible with the oil, displacement of oil is less efficient than for the miscible case, and the gas bank travels through the sand with a velocity less than that of the driving gas. These data along with theories discussed imply that smaller banks and less total gas are required when the enriched gas and oil are miscible. INTRODUCTION Widespread application of enriched-gas drive to the recovery of oil rests upon a key factor — the use of limited quantities, or "banks", of enriched gas. At the present time, the value of liquefied petroleum gas or other enriching agents discourages their use in a continuous injection technique, or even in a large bank, except in a few isolated reservoirs. If small banks of enriched gas driven by methane were as effective in displacing oil as is continuous injection, the enriched-gas drive process might be applied to a larger number of reservoirs. Previous research on the mechanics of the enriched-gas drive process reported by Stone and Crurnpl and by Kehn, Pyndus and Gaskell has utilized continuous injection of enriched gas. This work has shown that two types of displacements occur. With gases containing sufficient intermediates. the oil is displaced misciblv and complete recovery is obtained from the regions swept. When gases are used which contain insufficient intermediate hydrocarbon for miscible displacement, oil is displaced immiscibly. In the latter type, selective solution of the intermediate hydrocarbons causes a swelling and reduction in viscosity of the oil and leads to an increased recovery over that obtained by dry-gas (methane) drive. The size of the enriched-gas bank necessary for efficient displacement of oil is determined by those factors which cause deterioration of the bank. A differentiation may be made between those factors which operate on a microscopic scale and those which act on a macroscopic scale. On the smaller scale, the enriched gas mixes in the direction of flow by diffusion and convection with the fluids immediately preceding and following it. On the larger scale, the gas may by-pass the oil by flowing through permeable streaks, by overriding the oil because of density difference, or by fingering because of unfavorable viscosity ratios. In such cases, the enriching material tends to mix with the oil both laterally and in the direction of flow. The increase in effective area available for diffusion and dispersion of the enriching components leads to a faster degradation of the bank and a need for a larger bank than is necessary for those cases in which no by-passing occurs. The effects of such macroscopic factors in the deterioration of enriched-gas banks have been reported in a separate paper by Blackwell, Terry and Rayne. The present study was confined to the factors which operate on the smaller scale, in particular to the behavior of banks of enriched gas in sands uniformly swept by the gas. Experiments were designed to answer the following questions. 1. Can small banks of enriched gas driven by methane be used to secure oil recoveries comparable to those obtained by continuous injection of enriched gas? 2. What is the optimum bank size (the minimum bank size necessary to obtain a recovery comparable to that obtained by continuous injection of the same enriched gas)? 3. How many total pore volumes of gas must be injected to obtain the maximum recovery when the optimum bank size is used? 4. What is the effect of varying the number of enriching components in a gas bank? This report describes the experimental investigations and discusses the results in terms of their significance to reservoir behavior.

By C. van der Poel

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

By D. W. Reed, R. L. Reed, Tracht

Laboratory experiments on the reverse combustion of tar sands in a linear adiabatic system have shown that a highly upgraded oil can be produced from an exceedingly viscous, immobile oil. The dependence on the air-injection rate of peak temperature, combustion-zone velocity, oil recovery, air-oil ratio, residual coke and oil, fuel burned and distribution of product gases is shown graphically. Eflects of initial temperature, oxygen concentration, oil saturation and heat loss are discussed. Experiments bearing on the coking properties of heavy oils are mentioned and some results exhibited. Field application of the process hinges on the existence of adequate air permeability and the rate of reaction under reservoir conditions. INTRODUCTION It has been established that oil can be recovered from underground reservoirs by means of at least two fundamentally distinct methods involving in situ combustion of a certain fraction of the oil. Characteristic of both of these known methods is the production of oil from one or more wells by means of hot gases formed when a high-temperature reaction zone is advanced through the reservoir. In both cases, the reaction zone is created by heating certain of the wells to a sufficiently high temperature prior to the introduction of air, and the zone is maintained and advanced through the reservoir by appropriate control of the air-injection rate. In the first of these methods, which is called "forward combustion",' the combustion zone advances in a direction which is generally the same as that of the air flow; whereas in the second method, "reverse combustion",' the combustion zone moves in a direction generally opposite to that of the air flow. Forward combustion, on the one hand, is an ideal combustion process in the sense that a minimum of the most undesirable fraction of the oil is consumed as fuel in the form of coke, a clean sand is left behind and generated heat is used as efficiently as possible. However, the applicability of forward combustion is limited. Since the products of combustion, vaporized oil and connate water must flow into relatively cold regions of the reservoir, there is an upper limit on the viscosity of oil which can be moved by this process in a practical and economical fashion.' On the other hand, it is characteristic of reverse combustion that the vaporized oil and water together with the products of combustion are produced through sand which is already hot and has had its mobile liquid content eliminated. This means there is no upper limit on oil viscosity; indeed, the oil may be an entirely immobile semi-solid. However, fuel for the process is an intermediate fraction of the original oil, and the most undesirable fraction remains on the sand surface as a substantial deposit of coke. Since this coked material is not burned during reverse combustion, it represents energy which is available for the production of oil but is not used for this purpose. It follows that one can expect economics to be somewhat less attractive with reverse combustion than with forward combustion. Nevertheless, it is a process which is designed for reservoirs where forward combustion is impossible and, as such, has become a subject of experimental and theoretical investigation. In this paper, only experiments made with tar sands are discussed. DESCRIPTION OF THE PROCESS We proceed, then, to consider the process of reverse combustion in greater detail. Fig. 1 illustrates a temperature profile defining a combustion zone which moves from right to left when air flows from left to right. In Zone 1, the temperature is the initial reservoir temperature, and the tar sand is as yet unaltered. This statement must be modified to the extent that physical properties of the oil may be changed by low-temperature oxidation at reservoir temperature. As air passes into Zone 2, which has been warmed by conduction, it assists in vaporization of the very light ends (if there are any), and oxidation occurs at a significant rate. In this region, there is almost no production of carbon monoxide or carbon dioxide because predominantly addition-type reactions take place with the formation of oxygenated compounds such as aldehydes and acids together with water. The hydrocarbon-enriched and slightly oxygen-depleted gas stream enters Zone 3 where

I would like to congratulate the authors on an interesting paper and, in particular, on their justification of the use of superposition in treating gas wells. They state, however, that no theory has previously been developed lor treating the variable rate drawdown case. In a paper published a few years ago by my colleagues and myself, a method is described which allows a constant rate drawdown curve to be derived for any well. regardless of real rate variations. This method has the great advantage that it is not necessary to make any assumptions regarding the form or uniformity of the reservoir, and the derivation suggested by the present authors is simply an idealized case of our more general treatment. In the light of several years' experience in the derivation of constant rate drawdown curves: I would like to point out that the method proposed by Odeh and Jones may lead to confusing answers, particularly if applied over an appreciable time interval. The reason is that, if pressures are measured when a well is flowing, non-linearities in plots like Figs. 1 and 2 can arise due to either (1) changes in the skin factor with time and/or flow rate, or (2) an inhomogeneous reservoir. Separation of these two effects is only possible if some of the pressures used in the analysis are measured when the well flow rate is 0 or, in other words. during brief shut-in periods. REFERENCES 1. Arthur, K. B.. Metunomski, Z. G. and Ruche, Ms "Intetet ct methode d'utilisa tion des mesures de pression [)our Petude des gisemonts" Revue de T'I.F.P. (Nov., 1962) XVII, 1299. 2. Thomere, R. and Arthur, K. R.: "Analysis of the Pressure Bet~avior of a Production Well Subject to Interterence and Its Application to Reservoir Models", Paper SPE 1189 presented at SPE :Amund Fall Mecting Ireld in Denwer, Colo. (Oct. 3-6, 1965).

Jan 1, 1966

By J. P. Heller

The purpose of this discussion is to comment briefly on one of the points raised in the paper. The topic is that referred to in the Conclusion 2, in which the authors note a "disproportionate increase in mixing rate" as they proceed from a displacement performed at viscosity ratio 0.998 to one at 1.002. The recovery curve of this displacement, Fig.2A, shows a quite noticeable difference between the two runs in which displaced and displacing liquids were interchanged. This difference is ascribed by the authors to the change in viscosity ratio between the two experiments. This interpretation of the experimental results must certainly be based on rather precise measurements of the fluid viscosities or, at least, of the ratio of their viscosities. The ordinary Ostwald viscometer, employed in the usual way, is only just able to distinguish between liquids whose viscosities differ by two or three parts per thousand. The experimental refinements which then must have been necessary to measure with any precision a viscosity ratio of 1.002 would be worthy of mention. Indeed, one would have expected the authors to devote at least a paragraph to the display of this technique. For the same reason - to emphasize that the details of such unusually critical experimental work should be reported and discussed in the literature — I would have liked to see some comment by the authors on the requirements for temperature control during both the viscosity measurements and the displacement experiments. The temperature coefficient of viscosity of Soltrol C at room temperature has been measured in this laboratory at about - 1.9%/1°C. How much did the coefficients for the fluid mixtures used in the experiment differ from this and from each other Even if the two fluids were identical in this respect, it would appear that a temperature difference of .l°C between them in the flow system would produce a viscosity con- trast equal to that imposed by their composition. Reference to the experimental solutions of these problems would render Conclusion 2 more credible. Published theoretical work on the stability problem contains no indication of such extreme sensitivity to the viscosity ratio in the immediate neighborhood of unity. On the contrary, perturbation calculations such as those made by Chuoke1 and Perrine2 lead to initial finger growth rates proportional to the first power of the viscosity difference. The later instability growth calculated by Scheidegger,3 using the assumption of zero pressure gradient perpendicular to the average flow, is also simply proportional to the viscosity difference. The conflict regarding the sensitivity of the dependence of instability growth or mixing rate on the viscosity contrast needs explicit recognition in the literature. In particular, it would seem that the experimental work reported here should be supported by a more detailed description of the viscosity controls. Further discussion of possible alternative interpretations of the data would also be in order. For instance, is it possible that the unconsolidstted bead pack used changed its characteristics between the runs (due, perhaps, to the handling and inverting involved in the injection of the front) Or, are the data points of Fig. 2A averages from several runs performed in a series in which the two liquids alternated as displaced and displacing fluids? To summarize, it seems to me that the relationship asserted by Conclusion 2 in the paper needs the support of more detailed experimental evidence. REFERENCES 1. Chouke, R. L., van Mews, P. and van der Poel, C.: "The Instability of Slow Immiscible Viscous Liquid-Liquid Displacements in Porous Media", Trans., AIME (1959) Vol. 216, 188.