Search Documents

By S. W. Churchill

A digital computer was used to obtain an exact numerical solution for the transient behavior of the insulation and earth adjacent to an isothermal, submerged flat surface for a single set of parametric values. Comparison of the computed results with analytical solutions for limiting conditions revealed that a complete and general solution for all parametric values could be constructed from these limiting solutions. Complete and general solutions for insulated spheres and cylinders can also be constructed from limiting solutions. The procedure is illustrated in part for an insulated spherical tank. The heat flow was found to fall rapidly to a pseudo-steady state; after some time, it then decreased slowly to zero for a flat plate and cylinder, and to a low steady-state rate for a sphere. The accumulative heat flow during the initial falling-rate period may be a significant fraction of the heat flow during the entire first year. INTRODUCTION Underground storage of liquefied natural gas and liquefied petroleum gas has received considerable recent attention. The rate of heat flow from the earth to the storage cavity and the resulting temperature field in the earth are important factors in the technical and economic evaluation of such storage facilities. The objective of this paper is to indicate how complete solutions can be developed for the transient flux and temperature field for various geometrical configurations. The representative properties and dimensions, and the resulting parameters utilized in the illustrative results, are indicated in Table 1. The results presented herein are for dry earth. The effect of the latent heat of solidification of water in the soiI will be described in a subsequent paper. The temperature field in the insulation and earth is determined by energy balances, boundary conditions and initial conditions. The physical problem can be described mathematically with virtually no idealizations insofar as physical properties are known. It appears possible to solve the equations by numerical integration on a high-speed digital computer for any geometrical configuration and conditions. For complex situations, however, the computations are expensive and the results are highly specific. Analytical solutions have been developed for a few simple but important geometrical configurations and conditions, including one-dimensional heat transfer from (or to) earth at an initially uniform temperature to (or from) flat plates, spheres and circular cylinders. The solution for an insulated flat plate has also been derived but is in the form of a slowly converging series involving taklated functions. It was planned to use a computing machine to evaluate this series for a number of representative conditions. However, upon examination of the results of preliminary computations it was discovered that a complete, general and accurate solution could be developed by interpolation between several much simpler solutions for limiting conditions. This technique then was used to develop a a complete solution for an insulated sphere. The equations presented herein are derived in Carslaw and Jaegerl and other books on applied mathematics, or they are simple extensions of these previous results. Hence, all derivations are omitted. Emphasis is instead placed on the results and their

By N. Marusov, D. K. Lowe, J. C. Karp

This paper considers, from an engineering viewpoint, several factors involved in creating, designing and locating horizontal barriers for controlling water coning. This is an effort to consolidate new concepts with previous information so that a reasonable selection of barrier materials, dimensions and vertical position can be made. Coning theories previously developed are briefly reviewed and an effort i,s made to reduce the results of coning-theory calculations to a point where routine calculations can be made with a desk calculator. It is expected that these simplified calculations will give usable predictions of the amounts of improvement attainable with barriers of various dimensions. Apparatus and procedures used for testing the suitability of commercial cements are described and test results are presented. INTRODUCTION Just about as soon as the phenomenon of water coning was recognized as an oilfield problem, there were suggestions that the production of water coning could be controlled or completely suppressed by means of horizontal barriers. It was also suggested that natural barriers such as shale streaks were helpful in restricting bottom-water production. The implication was that wells not having the benefit of continuous shale streaks to suppress bottom waters should, in some way, be supplied with an artificial barrier. With regard to the selection of barrier materials, there are many references which give the merits of liquid barriers.',' These liquids include surface-active agents," recipitates' and emu1sions. There is no real need to review the merits of these materials again. This paper deals with the use of solid barriers which could be made from various cements. The two major restrictions placed on these barrier cements are that (1) they must be commercially available in large quantities, and (2) they must be applicable by means of standard tools and equipment with techniques that have been fully developed and are in general use at the present time. A schematic drawing of a horizontal barrier in a reservoir is shown in Fig. 1. Two maximum stable water cones are illustrated-—-one for the wellbore diameter and the other for a barrier. A maximum stable cone is one that has attained its maximum volume at the critical production rate and is on the threshold of producing water. The volume displaced by each water cone and the rise in water table are indications of the amount of water-free oil produced. Corrections for porosity, connate water and residual oil, of course, must be applied. The water-oil interface which defines a maximum stable cone is called a free surface. It is a surface over which the pressure is constant. This free surface has the shape of a limiting streamline at and below which there is no movement or transport. The shape of the water cone describes the pressure drop in the reservoir. The effect of a completely impermeable barrier on the cone shape is essentially the same as extending the wellbore out to the barrier radius. The placement of a horizontal barrier requires that a horizontal fracture be created using a single-point entry technique. Two methods may be employed for filling this fracture. One method consists of propping the fracture before filling with cement. The other consists of fracturing with fracture fluid, then immediately following with cement

By E. B. Brauer, B. F. Pennington, E. W. Hough, G. L. Stegemeier

Interfacial tension divided by the difference in density between the liquid and the vapor phases was determined experimentally by the pendant drop method on several isotherms in the two phase region below the critical point for the methane-normal decane system. The density difference data of Sage and Lacey was used in the calculation of interfacial tension. Both interfacial tension and and interfacial tension divided by density difference were found to vanish at the critical point. Interfacial tensions of less than one dyne/centimeter were observed as far as 1,000 pounds per square inch below the critical pressure. EXPERIMENTAL PROCEDURE The interfacial tension divided by the density difference for the methane-normal decane system was determined at the 100°, 130°, 160' and 190°F isotherms, from pressures of about 1,000 psi to the critical pressure, which is more than 5,000 psi for these isotherms. Particular emphasis was placed upon the investigation at pressures slightly below the critical pressure where the interfacial tension is less than 0.5 dyne/cm. Volumetric properties in the two-phase region, including the critical pressures and temperatures, were taken from the work of Sage and Lacey.1 The experimental pendant-drop technique used for the determination of interfacial tensions at high pressures incorporated the ideas of Michaels and Hauser,2 Hough, et a1,3 Walker4 and Heuer.5 In addition, the technique for determination of extremely small interfacial tension by Jennings6 was utilized in the region near the critical points. A detailed description of the apparatus is given in a dissertation by one of the authors.7 Cleaning operations on the stainless-steel sample system included successive washings with chromic acid, tap water and, finally, distilled water. Subsequent cleanings were performed with re-distilled normal pentane, which had an extremely low residue upon evaporation. Specific composition requirements necessitated a fairly precise sample introduction — although, fora two-phase, two-component system, the composition of each phase is completely determined if pressure and temperature are controlled. The normal decane was delivered into the evacuated sample system as a liquid from a burette. The methane was then introduced into the system from a calibrated isothermal container, so that pressure differentials could be used to determine the amount introduced. High pressures were obtained by compressing the sample with a mercury injection pump until the critical pressure was reached for the particular isotherm being studied. Experimental data were then obtained for specific pressures by first decreasing the pressure slightly so that two phases would appear, and then photographing a drop at that pressure. Subsequent photographs were made at increments throughout the pressure range. Calculation of interfacial tension divided by density difference was made from measurements

By L. D. Mullins, W. R. Bartlett, R. H. Barham, W. L. Wahl

This paper describes an electrical model and its applicution to the analysis of four reservoirs in Saudi Arabia. The model has 2,501 mesh points and represents 35,000 sq miles of the Arab-D member. Details of modeling such as mesh size, control problems and standards of performance in matching reservoir history are discussed. The particular performance match achieved for the Arab-D rnember is presented. Details such as permeability barriers, aquifer depletion and iwerference between oil fields are given. The performance match realized in the Abqaiq pool is presented in detail. INTRODUCTION The resistor-capacitor network and associated control equipment described in this paper comprise an electrical analog of a reservoir system. Similar equipment has been used to study the transient response of reservoirs for many years. The unique feature of the model and application to be described is the extremely large size of the model and reservoir system, and the detail observed in simulating the reservoir with the model. The Arabian American Oil Co. first became interested in analog computers for simulation of oil reservoirs in 1949. Since that time, several models have been developed, each more elaborate and refined so that the reservoir system might be more closely simulated. The current model is the latest in a series designed, built and operated by the Field Research Laboratory of Socony Mobil Oil Co. in collaboration with Aramco. It has been and continues to be used to study the regional performance of the Arab-D member limestone reservoir. The Arab-D member is one of the Middle East's most prolific producing horizons.. THE MODEL The theory of simulating a reservoir system with an electrical system has been presented in the literature.1-5 Therefore, this paper will not discuss the theoretical aspect of the problem except to point out the correspondence between the fluid system and electrical system, as shown in Table 1. In general, the complete model is made up of input devices, output devices, central control and a resistance-capacitance (RC) network. At times, the RC network alone is referred to as the "model". However, it should be evident from the text which meaning is attached to the word "model". A discussion of the equipment follows. THE RESISTANCE-CAPACITANCE NETWORK The RC network consists of 2,501 capacitance decades interconnected through 4,900 resistance decades. The components are arranged to form a rectangular network of 2,501 mesh points in a 41- X 61-mesh array. Imposing the mesh grid system on the continuous reservoir system divides the reservoir into discrete areal segments. These discrete segments may be of various sizes. More precisely, the mesh size need not be uniform throughout the model. The RC network is fabricated in two sections which are connected at the top. An inside view of the "tunnel" formed by the two sections is shown in Fig. 1. The height and width of the tunnel are shown in the figure. Numerals appear along the bottom and along the back opening of the tunnel. These numbers denote the x and y co-ordinate positions of mesh points. Fig. 2 presents a rear view of one-half the model. The length dimensions of the model, as well as a rear view of the capacitor decades, are shown in this figure. The control dials used in adjusting the resistance and capacitance values on the model can be seen in the enlarged portion of the model shown in Fig. 3. The electrical capacity at any mesh point can range from 0 to 1.0 microfarads set to the nearest tenth of a microfarad. The electric resistance connecting any two mesh points can range from 0 to 9,990,000 ohms set to the nearest 1,000 ohms. External capacitors may be added to any or all mesh points if the need arises. The values of electrical resistance and capacitance are adjusted manually by manipulating the two types of decade units. INPUT EQUIPMENT A considerable quantity of equipment is used to control the input to the RC network. These input devices are

By M. R. Tek

A new method for correlating the data on multiphase flow through vertical pipe is presented. The correlation is based on a "two-phase f factor" concept which was developed and successfully applied to horizontal multiphase flow by Bertuzzi, Tek and Poettmann. Field data previously published on several flowing and gas-lift wells have been used as the basis for the developed correlation. The application of the method to actual two-phase flow problems indicate that this method is capable of predicting the pressure distribution in vertical multiphase flow strings well within the accuracy range usually desired for common engineering and design calculations. A new working chart developed for calculation of two-phase pressure gradients and a graphical step-by-step procedure for the computation of pressure distribution are presented along with an example problem. INTRODUCTION Multiphase flow through vertical pipe is encountered in many engineering installations. In petroleum, chemical process, nuclear engineering and many other industries, problems associated with simultaneous flow of two or more phases through vertical pipe have been of interest for a long time. This interest has increased considerably during recent years due to applications to new processes in petroleum production and refining and to problems of steam generation and heat removal from nuclear reactors. One prominent example of vertical two-phase flow is provided by the gas-lift process where oil, water and gas flow simultaneously. If the pressure profile in a gas-lift well can be predicted within reasonable accuracy, it would be possible to get good estimates of the power required to lift the oil, the optimum depth, the pressure and the rate at which to inject gas. Furthermore, the effect of production rate and tubing sizes on these quantities can be evaluated before any design decision is made on the installation and operation of the flow string. The majority of experimental work available in the literature deals with two-component systems where individual phase flow rates in and out of the pipe remain constant. The general problem of prediction of pressure drops in multiphase flow systems is very complicated. The co-existence of numerous flow patterns of widely different geometry and mechanism, conditions of surface instability and the nature of force fields acting upon the system are among the major difficulties commonly encountered. The classical approach of fluid dynamics which would be based upon the formulation and solution of Navier-Stokes equations has been found by many investigators completely devoid of any hope of success-—not only because of inherent nonlinearities but also because of insurmountable analytical difficulties standirng in the way of setting up the boundary conditions. The presence and effect of interfacial forces on mu1tiphase flow systems further complicate the theoretical approach. For these reasons, many investigators choose to adopt semi-empirical if not purely empirical approaches in order to obtain solutions of engineering utility. A prominent practical solution has been given by Poettmann and Carpenter' in the form of an empirical correlation. In their paper, total flowing densities of fluids and solubility effects of gas in oil have been taken into account for the correlation of field data covering a wide range of operating conditions for oil wells. They treated the gas. oil and water as a single phase of combined propcrties and correlated the multiphase friction factor as a function of the product (pvd), i.e., density X velocity X diameter of the flow string. However, because the product (pvd) is dimensional while the frictior factor is dimensionless, the generality of their result's is somewhat restricted. It seems that the omission 01 viscosity effect may be one of the reasons for the scatter of data as shown in their correlation chart. The concepts of the "two-phase f factor" and the "two-phase Reynolds' number function" were recently developed and successfully applied to correlate horizontal multphase flow by Bertuzzi, Tek and Poettmann.' Recently, two new methods of correlation by Tek and Chan' have been presented on simultaneous flow of liquid and gas through vertical pipe. These two methods of correlation and the working equations in field units necessary for their application are included in this paper. The extension of the conccpt of two-phasc Reynolds' number function successfully developed for horizontal flow into vertical multiphase flow systems, development and evaluation of working charts permitting calculation of two-phase pressure gradients and

By H. D. Attra

The effects of phase behavior on reservoir economics is a function of both fluid composition and flow properties of the reservoir rock. In some operations such as gas injection into volatile crude-oil resertoirs, the stock-tank recovery derived from vaporization of the reservoir oil by the displacing gas can approach the recovery derived from the displacement mechanism. In other cases of gas-cycling reservoirs, the dry injected gas can effectively vaporize a retrograde liquid phase, allowing a substantial reduction in cycling pressure with no appreciable loss in recoverable liquids. Because of these effects, it is essential that phase behavior be considered in developing the economics of any secondary recovery program in which gas is the injected fluid. The inclusion of phase behavior in conventional calculation procedures has been retarded because of the complexity of the calculation and the large volume of laboratory data required. However, by applying high-speed computers and new laboratory techniques, the phase-behavior displacement calculation has been developed to economically include this basic fundamental of gas displacement. The calculation is performed by reducing the reservoir to a series of one-dimensional segments. Phase changes, saturation distributions and changing fluid properties are evaluated for each segment as a function of time. The relative volumes of equilibrium liquid and gas displaced from each segment are determined by a trial-and-error procedure to meet the requirements of phase equilibrium and two-phase flow. The produced fluids are flashed through conditions of surface separation approximating those in the field to obtain actual values for stock-tank liquid recovery and producing gas-oil ratios. Based on the results of successful field studies, the technique appears to have application to reservoirs producing under a nonmiscible gas-displacement mechanism and should substantially increase the reliability of this type of analysis when phase behavior is an important factor. Predicting reservoir performance under gas-injection operations involves three basic principles of reservoir engineering. First, a material balance exists between the displacing gas and the displaced reservoir fluid. Second, the oil and gas flow in relationship with the two-phase flow property of the reservoir rock; and third, phase equilibrium exists between the displacing gas and the native reservoir fluid. To illustrate these three principles, consider the one-dimensional reservoir shown in Fig. 1. The reservoir is divided into a series of segments representing a one-dimensional cross section between an injection area and a producing area. So represents the oil saturation, Cw represents the connate-water saturation and Sg represents the displacing gas saturation. The gas is moving into Segment 1 and the reservoir oil is being produced out of Segment N. Consider the case when gas has just moved into Segment 1 and the displacing gas-oil contact is at position a-a The volume of oil pushed out of Segment 1 is a function only of mobility and the two-phase flow properties of the reservoir rock, that is, Kg/Ko and µo/µg. Behind the front if the injected gas is foreign to the reservoir oil (a nonequilibrium gas), there will either be vaporization of the oil into the gas or condensation of the gas into the oil. In either case the net result is a mass transfer and a re-distribution of composition and liquid and gas saturation. If the displacement is at constant pressure and is continued, the gas moves into new segments, proceeding succesively into Segment 2, 3, etc. The total volume of oil produced from the reservoir

By L. W. Holm, A. K. Csaszar

A series of laboratory experiments was conducted in which oil was displaced from a porous medium by water-driven slugs of alcohols or similar solvents. The solvents used were soluble to some degree in both oil and water and covered the range of solubilities from complete solubility in oil to complete solubility in water. Displacement experiments were conducted on 2-and 3%-in. in diameter consolidated cores and 1-and 2-in. in diameter unconsolidated sand packs. The cores and sand packs ranged in length from 1 to 30 ft, and they were saturated with brine and crude or refined oils. The solvents used included ethyl alcohol, isopropyl alcohol (IPA), tertiary butyl alcohol (TBAI, secondary butyl alcohol (SBA), n-amyl alcohol (NAA), methyl ethyl ketone, acetone and others. It was found that all of the oil present in a porous medium could be miscibly displaced by injecting a slug of mutual solvent and driving it with water. The oil-recovery efficiency was dependent upon (1) the relative solubilities of the solvent in oil and water, and (2) the distance traversed by the flood. For complete oil recovery from cores, a smaller amount of a preferentially oil-soluble solvent was required, compared to the amount of preferentially water-soluble solvent needed. The size of the solvent slug required varied inversely with the linear flooding-path length raised to the 0.65 power. Water-driven dual solvent combinations (an oil-soluble solvent slug followed by a water-soluble solvent slug) were found to effect complete oil recovery with less total solvent than any single solvent used. In these dual-solvent displacement experiments, the slug size required varied inversely with length raised to the 0.55 power. Based upon the experimental results, a theory was developed to describe the displacement of oil and water by mutual solvents, and equations are presented to predict the production history in a linear system. These equations take into account the properties of the solvents and the porous medium. INTRODUCTION Oil-recovery processes which utilize displacing fluids that are miscible with the reservoir fluids have been studied extensively in recent years. Because of the poor contact efficiency and high pressure requirements of the LPG-gas displacement process there has been considerable interest in the alcohol-water process, and a number of studies have been made on the recovery of oil through the use of solvents which are mutually soluble in oil and water. An investigation by Sievert, Dew and Conleyl indicated that the use of mutual solvents would be limited because the presence of water in a porous medium would cause a phase break in the leading edge of the displacing solvent. Their study also showed that, in consolidated cores containing. only oil, the displacement of oil by a water-driven mutual-solvent slug of tertiary butyl alcohol (TBA) was affected by the viscosity ratios of the fluids involved. Gatlin and slobod3 concluded that an isopropyl alcohol (IPA) slug acts as a miscible piston, completely displacing both oil and water until the alcohol content of the mixing zone falls below the concentration necessary to maintain miscibility. Their study was conducted on uniform unconsolidated sand packs. They concluded further that IPA could be used effectively to recover oil from a watered-out sand. In a paper by Taber, Gamath and Reed4 relating an investigation on sandstone cores, it was stated that the displacement of oil by mutual solvents, particularly IPA, was not a miscible displacement and that no improvement in efficiency could be expected with increase in flooding-path length. However, their analytical analysis of the displacement mechanism using TBA is, in fact, one which indicates that the displacement is controlled by miscible mixing. They suggested that the lack of improvement in efficiency with flooding - path

By M. Prats, D. G. Russell

The performance of a well in a bounded, layered reservoir with interlayer crossflow has been investigated mathematically. The system studied comprises a centrally located well in a bounded cylindrical reservoir composed of two layers of contrasting physical properties. The reservoir contains a single fluid of small and constant compressibility. In the case of a well producing at constant pressure, the production rate of the crossflow system declines exponentially after the early transient behavior dies out. For wells which produce with constant rate, the pressure-time relationship is linear following the transient. For very early production time, the production rate of the two-layer crossflow system behaves as that of the two-layer stratified system (the case in which the layers are in communication at the well only). Except for early time, in most practical cases the performance of a two-layer reservoir with cross-flow can be duplicated by that of a single-layer reservoir. The equivalent single-layer reservoir must have the same pore volume and the same drainage and wellbore radii as in the two-layer case. In addition, the equivalent reservoir must have a "kh" (permeability x thickness) product equal to the sum of the "kh" products of the respective layers in the crossflow system. If the porosity and permeability distributions are known, the production and pressure behavior of a well with and without crossflow can be predicted. From a comparison of the observed pressure or production performance of a well and the theoretical behavior, the occurrence of interlayer crossflow can be detected. INTRODUCTION When flow is possible across the bedding plane from a low-permeability layer into an adjacent layer of moderate or high permeability, the earlier depletion of the more permeable layer will cause such flow to occur. If this situation exists, the oil recovery from tight layers can be appreciable; thus, the economic ramifications of interlayer crossflow are apparent. The possible importance and extent of flow between adjacent and connected layers — whether resulting from capillary, gravitational, or viscous forces — have been the subject of considerable speculation for many years. Jacquard1 presented a rigorous analysis of the crossflow due to single-fluid expansion occurring between two adjacent and connected layers when the well produces at constant rate. His presentation was entirely mathematical, however, and no numerical results from evaluation of the mathematical solutions were presented. Our objectives in the preparation and presentation of the present paper are essentially twofold— (1) to analyze the behavior of a well which produces a single compressible fluid from a bounded, layered reservoir with interlayer crossflow, and (2) to present practical numerical results and simplified formulas to facilitate engineering analysis of reservoirs with crossflow. The constant-terminal-pressure case was analyzed rigorously, whereas the constant-terminal-rate case was treated only for semisteady-state conditions. DESCRIPTION OF THE IDEALIZED TWO-LAYER RESERVOIR The reservoir to be studied is assumed to be bounded and cylindrical in shape, with a completely penetrating well at its center. The reservoir is composed of two adjacent horizontal layers between which fluid is free to flow. Each layer is assumed to be homogeneous, isotropic and completely filled with a fluid of small and constant compressibility. The total thickness of the reservoir is small enough

By J. B. Agan, R. J. Fernandes

This paper utilizes the layered-system approach, modified to include areal sweep efficiency, to determine the miscible-slug size required and to predict the perfort7zance of a miscible-slug process in a highly stratified reservoir having a permeability variation of 0.76. The slug size in each layer was determined from theoretical calculations and laboratory data which showed that, to maintain miscibility in the reservoir, the slug size had to be 8.2 per cent of the floodable hydrocarbon pore volume. Economic considerations, however, such as the value of the miscible product injected, the cost of injection and the cost of reprocessing the miscible product, limited the amount of miscible product which could be injected into the lower-permeability layers. This resulted in subjecting only 60 per cent of the reservoir to a slug size sufficient to maintain miscibility. Utilizing the relationship of producing gas-oil ratio and areal sweep efficiency as a function of pore volume injected (as obtained from data published by Dyes, Caudle and Erickson2), the theoretical performance was determined. This prediction is compared with actual performance to date. Actual performance reflects an increasing percentage of tniscible product in the produced fluid to a particular level, which remains static for a short period and then increases again to a higher level. This performance suggests that the reservoir has performed as a layered system. It is concluded that the layered-system approach can be used to predict performance in a highly stratified reservoir and that economic consideration is a limiting factor in determining the amount of the reservoir which can be miscibly swept. INTRODUCTION The Paloma zone in the Paloma oil field, Kern County, Calif., is comprised of several sand members separated by siltstone, shale and/or cherty sections. Rapid facies changes within the zone indicate that the sediments were deposited in a "near shore" environment. As a result, the Paloma sands have a high permeability variation, and the Paloma zone is comprised of several sand reservoirs. Since one of the sand reservoirs containing crude oil was not receiving the benefits of gas injection, consideration was given to recompleting certain wells to isolate the sand in order to perform a miscible-phase injection project and, thereby, to increase its ultimate oil recovery. Before proceeding with the project, a review of the literature was made to determine the miscible-slug size that would be required to maintain miscibility throughout the displacement process. Although the review showed that considerable information was available for homogeneous systems, insufficient information was available to determine slug-size requirements and to predict performances for a highly stratified reservoir. This paper has been prepared to describe the method used to determine the miscible-slug-size requirement and to predict the performance of the miscible-phase injection project initiated in the Paloma oil field. Although the Paloma project has only been in operation since Aug., 1958, sufficient production history has been obtained so that a comparison can be made with the predicted performance. DESCRIPTION OF RESERVOIR The reservoir in which the miscible-phase project was initiated in the Paloma field is shown on Figs. 1 and 2. Reservoir properties are as follows: average sand thickness, 43 ft; porosity, 19.1 per cent; water saturation, 47.2 per cent: oil viscosity, 0.328 cp; reservoir temperature, 250°F; reservoir pressure, 3,500 psi: gas saturation, 21.2 per cent at time of miscible product injection; and total net hydrocarbon pore volume, 19,993,000 reservoir bbl. LIMITATIONS The method described herein utilizes a layer-system approach, where injection into and production from each layer is assumed to be proportional to the layer's fractional capacity. The effects of increased fluid mobility on injection and production rates, in layers which have had breakthrough, are not reflected in this method. There has been considerable controversy over the slug size necessary to maintain miscibility throughout the displacement process. In this method, the slug-size

By L. W. Emery

Undergrorlnd combustion operations were initiated in a 60-acre Bartlesville sand "shoe-string" reservoir in Allen Connty, Kans., in 1956. Tests in separate patterns were conducted using various co~nbinations of air and recycle gas to propagate combustion fronts from the injection toward the producing wells. These patterns were made up of 6 injection and 20 prodrrcing wells Gas and liquid prorluctiorz from each pattern was measured on an individual-well basis, and comparisons were made between the three patterns to ascertain the relative effects of injected gas composition on production behavior. Breakthrough of the combustion front at a well was characterized by an increase in water production from the well followed by an increase in bottomhole temperatrrre to approximately 250" F. After burning fronts had broken through at five producing wells, operations were terminated in 1960. From the total project approximately 79,000 bbl of oil were produced during thermal operations at a cumulative produced GOR of 23 Mcf/bbl. No appreciable change in the character of the produced crude was observed. Combustion in the reservoir was maintained with injected gas compositions ranging fronz 6 per cent oxygen in recycle gas to 100 per cent air. lnjectiotz of large quantities of recycle gas resulted in higher producing GOR's from offset wells than were measured from a pattern into ~vhich straight air ~vas injected. The air required to move the combustion front through I acre-ft of reservoir was computed to be 20 MMscf. This valrre was found to be relatively independent of the quantities of recycle gas injected. The recovery efficiency from the swept area was esti~izated to be about 59 per cent. Areas swept were similar in shape to tlzose obtained with a laboratory potentiometric model. Samples of sund taken from behind the burning front by coring indicated almost total oil removal from the sand. Petrographic analysis of the core samples indicated that the sand had been heated to peuk temperature of rlbout 1,200" F. No rignificant difference in peak temperature was forrnd in two areas where compositions of injected gas were quite different. Compression costs for thermal recovery were estimated to be $1,20/bhl of produced oil. INTRODUCTION The use of the "forward combustion" process as an oil recovery method has received a great deal of attention. This method involves ignition of the formation in an injection well, followed by propagation of a combustion front through the reservoir. Combustion is maintained by the injection of an oxygen-containing gas to react with reservoir hydrocarbons. As the flame front progresses through the reservoir, oil and formation water are vaporized, driven forward in the gaseous phase and recondensed in the cooler part of the formation. In turn, the condensed fluids push oil into the producing wellbores. Completed field tests of the process were first reported by Kuhn and Koch,' and by Grant and Szasz.' Results from other tests have since been reported by Walter,3 by Moss, White and McNeil,' and by Gates and Ramey." ach of these tests essentially utilized a single injection well surrounded by four or more producing wells. Sinclair Research, Inc., elected to do field experimental work using a number of test patterns in a single field in order that comparisons between various operating schemes could be made. The site selected and purchased in 1955 for this experimental work was a 60-acre Bartlesville sand reservoir located in Allen County, Kans. Combustion operations were initiated in mid-1956. Between that time and termination of the project in mid-1960, combustion fronts were propagated from injection wells to producers in three separate well patterns, using different mixtures of air and recycle gas. The test was terminated before sweep of the three patterns was complete so that information about the effect of combustion on the swept areas could be obtained by coring. Results from the test in the form of injection and producing well performance have been carefully recorded, and these form the general basis for this paper. DESCRIPTION OF RESERVOIR The reservoir in which the combustion tests were conducted is a Bartlesville sand "shoe-string", typical of a number of small reservoirs in Southeastern Kansas. Average reservoir characteristics are shown in Table 1. Fig. 1 is an isopachous map of the producing sand showing the reservoir to be approximately 400-ft wide and 2,500-ft long. Maximum net productive sand thickness is 21 ft. Fig. 2 shows a typical core analysis obtained by coring with water-base mud. The reservoir has no appreciable dip and is closed on the sides by degradation of sand into shale. The main body of sand is heavily laminated with shale stringers, which are not continuous between wells. The main reservoir is overlain by 30 to 40 ft of laminated low-permeability sand and shale streaks. No information is available on the original properties of

By P. E. Baker

Approximate solutions are presented for the heat-flow equations in a loss-free linear system with a moving source and with heat transfer by convection and conduction, representing in situ combustion in an oil reservoir. The rate of heat generation and the velocity vf of the source (flame front) are assumed constant. A cold front or convection front caused by the cooling effect of in-jected air develops behind the flame front and moves with velocity vc, which is determined by the thermal properties of the air and rock. Because of thermal conduction the fronts are not sharp; i.e., they do not represent temperature discontinuities. The following solutions apply after a sufficiently long time. where Q is heat generated per unit volume of formation, k is thermal conductivity, CF is volumetric heat capacity, a is thermal diffusivity, and subscripts 1 and 2 refer to values behind and ahead of the flame front, respectively. The parameter CF2 is contained in v . Te is initial reservoir temper- . C2 e ature and air-injection temperature. Other mathematically possible solutions are presented. The different solutions result from the use of different values for various physical parameters. The solution just presented is thought to be most commonly applicable. INTRODUCTION In this paper, approximate analytical expressions are derived for the temperature profile developed in the underground combustion method of oil recovery. A qualitative description of the process may be found in a paper by Tadema.1 The flame front is regarded as a moving source with a constant rate of heat generation. Heat transfer is by conduction and convection. Previous publications, of which a few are cited here, 2-5 have also attacked this problem and have presented solutions pertaining to different methods of approximating or idealizing the process studied. The present paper also treats an idealized case: linear flow is assumed for heat and fluids, and vaporization and condensation effects are neglected. The results obtained are similar to those derived by Bailey and Larkin,3 but are obtained by an entirely different procedure. The method used here should give considerably greater insight into the physics of the in situ combustion process. Four types of cases are considered, distinguished by thermal and fluid-flow properties of the system. A constant rate of heat generation and constant rate of flame-front advance are assumed. Flame-front temperature is determined for each case using a heat balance which involves the entire temperature profile. The applicability of the different solutions to real systems is discussed after they are all presented. Again, this is done for the purpose of understanding better the physics of the process. Consideration of the nonapplicable solutions and of the reasons for their not applying is significant. In a paper by Cooperman,4 the only case solved is one considered here as not applicable. Fluid flow, as such (pressures, saturations, oil recovery, etc.), is not treated. This subject has been covered in a paper by Wilson, et al.6 Also, we are not concerned with requirements for continuance of the process; we only attempt to determine the temperature history for cases where the process has continued for some time. PHYSICAL SYSTEM AND APPROACH TO PROBLEM A linear flow (heat and fluid) system is assumed and the x-axis is taken parallel to the flow direction. The combustion zone (flame front) is vertical,

By G. Thodos, R. B. Grieves

A method has been developed For predicting the critical temperatures and critical pressures of binary mixtures of carbon dioxide, hydrogen sulfide, nitrogen, hydrogen, carbon monoxide and oxygen with the normal paraffin hydrocarbons. For carbon-dioxide and hydrogen-sulfide systems, relations are presented that take into account the peculiar behavior of mixtures with closely boiling components, such as carbon dioxide-ethane and hydrogen sulfide-propane mixtures which exhibit minimum critical temperature points. For hydrogen, nitrogen and carbon-monoxide systems, the extreme critical behavior caused by wide differences in pure component properties is established In addition, those fixed gas-paraffin systems which resemble paraffin-paraffin systems are also accounted for. For a mixture of known composition, the pure component critical temperatures, critical pressures and normal boiling points are all that are required to determine its critical point. Graphical relations are presented relating Tc and PC , of the mixture to the pure component properties. From the treatment of 12 carbon-dioxide and hydrogen-sulfide systems reported in the literature (74 mixtures), the expected error lor the critical temperature is approximately 1.5 per cent, and for the critical pressure, approximately 2 per cent. From the treatment of six hydrogen, nitrogen and carbon-monoxide systems reported in the literature (30 mixtures), the expected error for both the critical temperature and critical pressure is approximately 2.5 per cent. The relationships, which have been developed with only normal paraffins as the hydrocarbon components, may be extended to those isoparaffins and olefins which fall within the allowable volatility ranges. INTRODUCTION Many of the fixed gases — carbon dioxide, hydrogen sulfide, nitrogen, hydrogen, carbon monoxide and oxygen — occur in natural mixtures with hydrocarbons. Carbon dioxide and hydrogen sulfide are frequent components of the fluids produced from underground petroleum reservoirs. Nitrogen, carbon dioxide and hydrogen sulfide are present in varying quantities in most natural gases and gas-condensate well effluents. Hydrogen mixtures are of considerable interest in many phases of refining processes of petroleum. The determination of the critical temperatures and critical pressures of such mixtures is of value in vapor-liquid equilibrium studies, for the prediction of the characteristics of underground reservoirs, and for reduced-state correlations of PVT, transport and thermodynamic properties. The accurate estimation of the critical point for binary mixtures is an important initial step toward a complete analysis for the establishment of the critical temperatures and pressures of multi-component mixtures. Methods for predicting the critical temperatures and critical pressures of binary hydrocarbon systems have already been presented in the literature.8 It is possible to apply these existing methods to fixed gas-paraffin mixtures, but due to their unusual critical behavior, values calculated deviate considerably from experimental values. For systems containing trace quantities of the fixed gases, these methods are acceptable; however, for systems containing more than 5 mol per cent of the fixed gases, these utterly fail to produce reasonable critical values. Consequently, in this study a method has been developed for handling such binary mixtures over the entire composition range. CARBON-DIOXIDE AND HYDROGEN-SULFIDE SYSTEMS The critical behavior and the vapor pressure behavior of mixtures of carbon dioxide and hydrogen sulfide with paraffinic hydrocarbons may be quite similar or quite dissimilar to that of paraffin-paraffin mixtures, depending on the volatilities of the components involved. The critical temperature and normal boiling point of carbon dioxide are very close to the corresponding values for ethane, while its critical pressure is considerably higher than that of ethane. The normal boiling point of

By R. L. Perrine

A stability theory is developed for miscible liquid-liquid displacement within a porous medium. In the usual case considered, a high-density high-viscosity "oil" is displaced downdip by a low-density low-viscosity "solvent'! Perturbation methods are used to find the conditions under which the spreading mechanism changes from the stable dispersion process to unstable viscous fingering. We find that instability is conditional, that there is a dependence on the shape of a disturbance leading to a"diameter" effect and that very difficult experimental scaling problems may result. A useful consequence is the definition of a minimum "slug size" for stable miscible displacement. This should make possible optimum use of the solvent process for oil recovery. The results may apply to many situations in which one fluid displaces another of somewhat different fluid properties within a porous medium. INTRODUCTION The process of miscible liquid-liquid displacement looked very favorable when it first received serious consideration as a means to increase petroleum recovery. Some early laboratory results were interpreted to mean that only a small "slug" of solvent was needed, perhaps 2 to 3 per cent of a hydrocarbon pore volume. This "slug" could recover the oil in an entire reservoir provided the fluids remained miscible.1,2 Only a slow growth of the "mixing zone", or gradual transition from oil to solvent, should occur. This would follow naturally from mixing by a dispersion mechanism such as that described by Scheidegger.3 And, indeed, laboratory cores have shown this kind of behavior provided solvent viscosity was at least as great as that of the oi1.2,4-6 The same kind of behavior has also been observed with solvent viscosity lower than that of the oil after the distance traversed became large. Not all laboratory results have been this favorable, however.7-9 Mixing zone growth may become proportional to the distance traveled rather than the square root of this quantity. Such behavior accompanies high rates in short systems of large diameter, provided the viscosity ratio is adverse. If this more rapid spreading were to persist in a reservoir, the solvent requirement for miscible "slug" displacement could exceed 30 per cent of a pore volume. There is considerable economic importance in the difference between 3 and 30 per cent solvent. A 1ogical conclusion is that two different kinds of flow behavior are possible, with each kind leading to a different result. One kind gives efficient displacement, the other does not. In the efficient case, the solvent bank spreads out only by the mechanism we have termed dispersion. Under these conditions, displacement is as near piston-like as possible. In the second kind of flow, found only with adverse viscosity ratios, viscous fingering occurs. That is, permeability variations cause a small finger of low-viscosity solvent-rich fluid to move ahead of its average position within the mixing zone. This creates a path of low resistance to flow, and an even greater amount of solvent-rich fluid follows. Thus, the process is autocatalytic. Once started, the fingering mechanism rapidly becomes dominant. The question to be answered is this. Under what conditions will each of these two different kinds of flow occur? In particular, what conditions denote the transition from dispersion to viscous fingering as the solvent spreading mechanism? The answer may tell us whether or not the desirable, near piston-like miscible displacement process is practical. An answer to this problem can be obtained from theory by the use of perturbation methods. The procedure is as follows. We first formulate a mathematical representation of the system. Then a small disturbance (or perturbation) in solvent concentration profile is introduced and observed to see what happens. Unstable flow is indicated when a small disturbance will grow larger. This will lead to eventual viscous fingering. If on the other hand the disturbance dies out, the displacement is stable. In this latter case, with dispersion as the spreading mechanism, near piston-like displacement is possible. Thus, this paper presents the development of a stability theory for miscible liquid-liquid displace-

By O. G. Kiel, G. W. Swift

The concept of "a continuous succession of steady-states", which has been applied successfully by Aromfsky and Jenkins to obtain a solution for the nonlinear partial differential equation describing the transient Darcy flow of gas through porous media, is demonstrated to be equally valid for transient non-Darcy flow. A mathematical model, which numerically solves the partial differential equation, is used to check the validity of the succession of steady-states solution. Comparison of sand-face pressure histories compelted by the two methods shows excellent agreement. The utility of the succession of steady-states solution in predicting performance of gas wells rests in the fact that no special computation equipment is required. The development of the succession of steady-states solution leads also to a practical method for determining and analyzing field test data. A method for taking gas-well test data under constant-rate conditions is presented. Experimental data obtained in the field by employing the constant-rate method are presented and analyzed in accordance with the succession of steady-states solution. Analysis of data in this fashion is demonstrated to give direct "in situ" information for reservoir permeability, porosity and turbulence or non-Darcy coefficient. INTRODUCTION The economics of gas production are dependent upon the transient behavior of flow within the reservoir. For production from a finite reservoir, the transient flow behavior can be subdivided into two parts. At first, the transient caused by the movement of the pressure "wave" into the reservoir is of importance. Later in the production history, the pressure-wave movement ceases and the second transient stage of material depletion becomes controlling. For reservoirs of relatively high permeability, it can be shown that the pressure wave moves into the reservoir and stabilizes quite rapidly. In the case of relatively impermeable reservoirs, quite the opposite is true. Although it is theoretically possible to compute the production capability of a well from the properties of the reservoir as determined by static tests and core analyses, more reliable information is obtained by conducting flow tests on the well and thereby obtaining some measure of "in situ" formation properties. For gas wells, there are two basic types of tests in existence: the flow-after-flow method,' and the isochronal method.' Both of these techniques are tailored to obtain data that can be analyzed in accordance with the empirical performance equation: In addition, the isochronal method makes provision for the sluggish nature of pressure-wave movement in "tight" formations by requiring pressure build-up between flows and by stipulating that data obtained on successive flows be analyzed at equal elapsed flow times. It can be demonstrated that either test is valid for reservoirs of high permeability. Further, since it has been pointed out that the pressure wave stabilizes rapidly for reservoirs of this type, tests of relatively short duration will give stabilized information on the performance of a well. Further decline of sand-face pressure and/or production rate may be determined by employing material-balance techniques. Cullender' points out that for relatively impermeable reservoirs the flow-after-flow method gives invalid results. (See Appendix C.) If the isochronal method of testing is used, there are two alternatives: (1) the tests must be conducted for a sufficient length of time to obtain stabilized information (which may require months to accomplish); or (2) some method for extrapolating the results of short-term isochronal tests must be employed. The first alternative is impracticable because of manpower, conservation and economic considerations. Recourse to the second alternative requires some assurance regarding the reliability of the extrapolation technique. Poettmann and SchilsonJ present an empirical method for predicting stabilized performance. The present investigation was originally initiated to determine the reliability of this technique. To do this, a mathematical model was developed to simulate the Darcy and non-Darcy flow of gas through porous media. The model consisted of a finite-difference approximation of the nonlinear partial differential equation which was solved on an IBM 7090 computer. Long-term production histories were simulated by the model and compared against predictions obtained from the Poettmann-Schilson method. As the work progressed, it became apparent that a straightforward predictive equation could be developed by utilizing the concept of a succession of steady-states. As a result, the emphasis of the work was redirected to exploit the advantages of the new method.

By G. Thodos, E. M. Phillips

A reduced density correlation for the gaseous and liquid regions of ethane has been developed from all the experimental PVT data available in the literature for this substance. Saturated vapor and liquid densities reported for methane, propane and n-butane were utilized to establish the rectilinear-diameter line for ethane. This line was then used with the saturated-vapor data for ethane to produce the saturated envelope for this substance. Most of the saturated liquid data reported for ethane were found to be inconsistent with the saturated-liquid curve resulting from this study. The correlation covers the regions 0.65 < TR <6.0 and 0 < PR < 15. For regions for which no experimental data have been reported, the isobars of the correlation were extended by the use of the theorem of corresponding states. This correlation is presented both on rectilinear and log-log coordinates. INTRODUCTION At the present time, there is considerable interest in the exact PVT behavior of pure hydrocarbons. Matschke and Thodos12 recently have developed a comprehensive compilation of experimental vapor and liquid density data for methane. Similar studies have been conducted for other pure substances.4, 7 In the present study the density data available in the literature for gaseous and liquid ethane have been critically reviewed and utilized to develop a reduced density correlation for this substance. SATURATED LIQUID AND VAPOR DENSITIES The rectilinear-diameter rule proposed by Cailletet and Mathias5 states that the sum of the saturated liquid and vapor densities of a pure substance is a linear function of temperature which passes through the critical point. This law, although empirical, has been well established by continued verification since its conception in 1886. This rule has been applied to the saturated density data of methane,6 propane 17 and n-butane8 to produce the straight lines presented in Fig. 1 which pass through the critical point TR = 1.00, Pr = 2.00. The rectilinear-diameter relationship for ethane was established from the lines for the other substances in Fig. 1 by plotting the quantity p~, + PR against zc for different reduced temperatures. A critical analysis of the saturated vapor and liquid densities reported for ethane by Kuenen,10 Maass and McIntosh,11 and Sage, Webster and Lacey18 revealed that the values of PR1 + pr resulting from these studies were inconsistent with the rectilinear-diameter line obtained from the reliable density data of the other substances, as shown in Fig. 1. In the analysis of these data, reduced quantities were calculated using the critical con-

By J. H. Moran, E. E. Finklea

The pressure build-up technique is a recognized method of determining permeability from conventional drillstem tests. In this paper an effort is made to extend such techniques to the interpretation of data obtained from the wireline formation tester. Such a study is necessary because of the differences, for this case, in the magnitude of the flow parameters (rate of flow, amount of recovered fluids) and in the flow geometry (flow through a perforation vs flow across the face of the wellbore, etc.) involved in the solution of the equations of flow for compressible fluids. The perforation is replaced by a spherical hole, and the effect of the borehole is neglected, so that the flow can be considered to be radial in a spherical co-ordinate system. Arguments are presented to justify this idealization. Assuming single-phase flow, general relations between pressure and flow rate are developed for a homogeneous medium. The study is then extended to permeable beds of finite thickness. It is shown that the early stages of pressure build-up tend towards spherical flow, while the later stages tend towards cylindrical flow. The thinner the bed, the more quickly flow approaches the cylindrical model. The prevalence of thin beds in practical work makes this analysis quite important. Cases involving permeability anisotropy are treated. INTRODUCTION From wireline formation tester operation, two types of data are obtained: (1) the nature and amount of recovered fluids, and (2) the pressure history recorded during the test. A number of papers have been written dealing with the interpretation of formation production on the basis of the recovered fluids.&apos;.&apos; In general, the methods described have been quite accurate for both high- and low-permeability formations. The present paper will deal with an analysis of the pressures observed. An analysis of the pressure build-up curves obtained in hard-rock country has already been attempted on the basis of the formula proposed by Hor-ner. Although this approach has met with success in many instances, some questions have been raised as to its validity. It is the aim of the present study to place the analysis of pressure build-up in the formation tester on a firmer basis, from which more detailed methods of interpretation can evolve. Because of the great differences between the operation of the wireline formation tester and the conventional drillstem test, modifications are necessary in the interpretation. The major difference relates to the flow geometry. Once the flow geometry has been established other features such as multiphase flow, skin effect, afterflow, etc., well described in the literature, can be introduced. It will be assumed that the mechanical operation of the formation tester is already known to the reader.6 t will suffice here merely to state that the tester provides the means for taking a relatively small sample of the fluid immediately adjacent to the borehole, and for recording the subsequent pressure response. In comparison with conventional drillstem tests, the time required for a satisfactory pressure build-up response is much shorter, because of the relatively small quantity of fluid withdrawn by the wireline tester. This feature is highly desirable in the case of low-permeability formations. For an analysis of the pressure response within the formation, three simple flow geometries are considered— linear, cylindrical and spherical. The spherical and cylindrical flow geometries are most pertinent to the formation tester; therefore, they will receive the major emphasis. Since the configuration of the borehole and the perforation made by the tester complicate the flow geometry, it is necessary to allow for them in the drawdown response. However, because of the volume of formations contributing to the pressure-response, the details of the perforation shape are unimportant in the build-up period. Since relatively small amounts of fluid are withdrawn from the formation, in contrast to a conventional drill-stem test, a study of the "depth of investigation" and the significance of drawdown as well as build-up data will be included. Because the "depth of investigation" will be shown to be rather large, the effect on the build-up curves of the finite thickness of the permeable bed is considered. It is this consideration that leads to the importance of cylindrical flow geometry. Also included is a discussion of permeability anisotropy and its effect on the interpretation of the tester results. The pressure curves recorded by the formation tester will follow two general patterns, depending upon whether the formation is of high or low permeability. Fig. I (a and b) schematically illustrates these two responses. In Fig. 1(a), the high pressure recorded during fill-up of the tool is essentially the pressure differential across the choke in the system. In Fig. l(b), the flow rate is

By H. D. Outmans

In steady vertical flow, the interface of an immiscible liquid-liquid displacement is horizontal for any flow rate below the critical. In nonvertical flow, however, the shape of the interface in the steady state does depend on the flow rate, and the purpose of this paper is to calculate the unsteady interfaces during the transition of one steady state of flow to another. A knowledge of these transient interfaces is of considerable importance in reservoir engineering where the calculation of breakthrough recovery depends on the instant the interface reaches the producing wells and on the shape of the interface at that time. Although the emphasis is put on transient interlaces, which eventually approach stable equilibrium, it is shown that if the displacement exceeds a critical rate no equilibrium is possible. The interface is then unstable and viscous fingers are formed during the displacement. The critical rate and the shape of the transient and equilibrium interfaces are affected by the effective interfacial tension; but since this effective interfacial tension appears in the calculations only in combination with the inverse square of the thickness of the medium, its effect in the reservoir would appear to be negligible compared to its significance in model experiments. INTRODUCTION Stability criteria and the early growth of interfacial disturbances in a plane parallel to the boundaries of a dipping formation in which oil is displaced by an incompressible fluid were described in a previous paper.l This type of instability is significant in thin reservoirs. However, if the reservoir has appreciable thickness, then interfacial stability in vertical planes, normal to the upper and lower boundaries, also becomes important (the displacement is supposed to be parallel to these vertical planes). The difference between the two stability problems is that, in the first case, the intersections of the interface with planes parallel to the boundaries are normal to the direction of the displacement; in the second case, the intersections, this time with vertical planes, are not normal to the displacement. Instead, they are tilted at an angle which depends on the displacement rate. The tilt of steady interfaces was calculated by Dietz2 who also determined the critical rate of displacement for stability in the vertical plane by assuming that this rate would coincide with an interfacial tilt equal to the dip of the formation. The critical rate thus calculated is the same as has been found for thin reservoirs (see Eq. 1.1 of Ref. 1 and of the present paper). Dietz&apos;s calculation of the stable tilt was verified by laboratory experiments and the agreement was found to be fairly good.3 It is doubtful, however, that stable tilts actually exist in the reservoir because a change in production rate is not followed by an instantaneous adjustment of the interface to the new rate but, rather, by a transition period during which the interface changes from one equilibrium tilt to the other. The principal objective of this paper has been to describe these transient interfaces without putting any restrictions on the flow conditions or the shape of the interface, as had been done previously. The second objective was to compute the critical velocity, taking into account capillary effects, and the third was to evaluate, at least qualitatively, the shape of the front at rates above the critical, again without making the simplifying assumptions introduced by previous investigators.2,3 In the following sections two examples are given of the calculation of interfacial motion. The first describes this motion for an initially horizontal interface in a dipping layer, and the second for a vertical interface in a horizontal layer. The mathematical formulation of the problem is nonlinear in the boundary conditions, and this prohibits its solution in closed form. Instead, the solution is obtained in the form of higher-order approximations. 1 Before proceeding to a description of the mathematical model, however, we define two quantities

By D. R. Parrish, K. W. Beaver, H. W. Wood, R. W. Rausch

A pilot test of forward combustion in the Shannon pool, Salt Creek field, Wyo., is described. The Shannon sand, 950-ft deep, contains a heavy (25" API), viscous (76 cp) oil. Natural reservoir energy is limited. Primary production, intermittent since 1889, recovered only about 2 per cent of the oil in place. The field is operated by Pan American Petroleum Corp. for the Midwest Oil Corp., the owner. The original pilot was a 1.32-acre five-spot. The expanded pilot has eight producing wells surrounding a roughly triangular area of about five acres with the injection well near the center. A control or comparison well was also recompleted in another part of the field. Operation of the pilot has been little different front an ordinary gas drive. Little special equipment was found to be absolutely necessary. Except for some use of a temperature-resistant cement, all wells were conventionally completed. In spite of poor oxygen consumption, the over-all performance of the pilot has been good. Total oil recovery to June 1, 1961, was 73,971 bbl. The wells of the original pilot alone had produced about 24,000 bbl, equivalent to 50 per cent of the oil in place, when fire breakthrough at the first well occurred. These wells have now produced oil equivalent to more than 74 per cent of the oil in place in the original pilot area and production is continuing. It appears that ultimate recovery will approach theoretical maximums before the wells must be abandoned. Performance of the pilot has been encouraging, and expansion to a fieldwide combustion operation is being investigated. INTRODUCTION The results of both laboratory investigations and field tests of underground combustion have been reported previously.&apos; However, most of the field tests were primarily experimental. More information and experience are needed before forward-combustion operations can be engineered with confidence. The purpose of this paper is to present the results of a successful pilot test of forward combustion. These re- sults should increase confidence in forward combustion as a practical method for commercial oil recovery. HISTORY OF THE SHANNON POOL The Shannon pool is located on the north end of the Salt Creek field in Natrona County, Wyo. It is approximately 50-miles north of the city of Casper. The Shannon pool&apos;s3 discovery well was completed in 1889, making this one of the oldest oil fields in the Rocky Mountain region. Three more wells were drilled in 1890. First production was hauled in wooden barrels by horse and wagon to the railroad in Casper. In 1894 a small refinery, the first in Wyoming, was built in Casper to process the Shannon crude. In the following years the field changed hands several times. It appears that each new owner did some development drilling as several wells were completed in each of the years 1895, 1902, 1905 and 1912, with negligible development in the intervening years. Forty-eight wells were drilled in this period, but many were later abandoned. Discovery of the more prolific Salt Creek field proper ultimately forced suspension of operations at the Shannon pool. After 1915 there was only sporadic production, mostly to supply cheap boiler fuel to drilling rigs in the Salt Creek field. But even this was discontinued in 1931. Since then the pool had been dormant until the recent operations, which are the subject of this paper. The Shannon pool is now owned by the Midwest Oil Corp. Field operations are conducted for Midwest by Pan American Petroleum Corp. THE RESERVOIR Fig. 1 is a map showing subsurface contours of the Shannon pool. The reservoir is on a nose of the Salt Creek anticline dipping to the north at about 500 ft/mile. The trap is provided by a shallow fault on the updip side of the productive area. The downdip limits are bounded by water, but this water has not provided an effective source of reservoir energy. The Shannon sand outcrops at many places in the immediate vicinity, providing good surface indications of the Salt Creek anticline. At the Shannon pool the sand has been lowered by faulting and is overlain by about 900 ft of shale and other sands. The Shannon sand consists of two members. The upper, a water sand, is about 40-ft thick and is separated from the lower member by several feet of sandy shale. The oil pay occurs in the lower mem-

By H. J. Ramey

As fluids move through a wellbore, there is transfer of heat between fluids and the earth due to the diflerence between fluid and geothermal temperatures. This type of heat transmission is involved in drilling and in all producing operations. In certain cases, quantitative knowledge of wellbore heat transmission is very important. This paper presents an approximate solution to the wellbore heat-transmission problem involved in injection of hot or cold fluids. The solution permits estimation of the temperature of fluids, tubing and casing as a function of depth and time. The result is expressed in simple algebraic form suitable for slide-rule calculation. The solution assumes that heat transfer in the wellbore is steady-state, while heat transfer to the earth will be unsteady radial conduction. Allowance is made for heat resistances in the wellbore. The method used may be applied to derivation of other heat problems such as flow through multiple strings in a wellbore. Comparisons of computed and field results are presented to establish the usefulness of the solution. INTRODUCTION During the past few years, considerable interest has been generated in hot-fluid-injection oil-recovery methods. These methods depend upon application of heat to a reservoir by means of a heat-transfer medium heated at the surface. Clearly, heat losses between the surface and the injection interval could be extremely important to this process. Not quite so obvious is the fact that every injection and production operation is accompanied by transmission of heal between wellbore fluids and the earth. Previously, the interpretation of temperature logs&apos;,&apos; has been the main purpose of wellbore heat studies. The only papers dealing specifically with long-time injection operations are those of Moss and White3 and Lesem, et al.&apos; The purpose of the present study is to investigate wellbore heat transmission to provide engineering methods useful in both production and injection operations, and basic techniques useful in all wellbore heat-transmission problems. The approach is similar to that of Moss and White:&apos; DEVELOPMENT The transient heat-transmission problem under consideration is as follows. Let us consider the injection of a fluid down the tubing in a well which is cased to the top of the injection interval. Assuming fluid is injected at known rates and surface temperatures, determine the temperature of the injected fluid as a function of depth anti time. Consideration of the heat transferred from the injected fluid to the formation leads to the following equations. For liquid, Eqs. 1, 1A and 2 are developed in the Appendix. These equations were developed under the assumption that physical and thermal properties of the earth and wellbore fluids do not vary with temperature, that heat will transfer radially in the earth and that heat transmission in the wellbore is rapid compared to heat flow in the formation and. thus, can be represented by steady-state solutions. Special cases of this development have been presented by Nowakl and Moss and White.3 Both references are recommended for excellent background material. Nowak&apos; presents very useful information concerning the effect of a shut-in period on subsequent temperatures. Since one purpose of this paper is to present methods which may be used to derive approximate solutions for heat-transmission problems associated to those specifically considered here, a brief discussion of associated heat problems is also presented in the Appendix. Analysis of the derivation presented in the Appendix will indicate that many terms can be re-defined to modify the solution for application to other problems. Before Eqs. 1, 1A and 2 can be used, it is necessary to consider the significance of the over-all heat-transfer coefficient U and the time function f(t). Thorough discussions of the concept of the over-all heat-transfer coefficient may be found in many references on heat transmission. See McAdams5 or Jakob," for example. Briefly, the over-all coefficient U considers the net- resistance to heat flow offered by fluid inside the tubing, the tubing wall, fluids or solids in the annulus, and the casing wall. The effect of radiant heat transfer from the tubing to the casing and resistance to heat flow caused by scale or wax on the tubing or casing may also be included in the over-all coefficient. According to McAdams, on page 136 of Ref. 5>

By O. R. Wagner, R. O. Leach, H. W. Wood, C. F. Harpke

PAN AMERICAN PETROLEUM CORP.TULSA, OKLA. CASPER, WYO A field test has been made in which additional oil recovery was obtained from a previously-waterflooded "oil-wet" sandstone reservoir. This recovery improvement was accomplished by adjusting the reservoir wettability through chemical treatment of the flood water. The test was made in the Muddy sand of the West Harrisburg Unit, Neb. The chemical used, sodium hydroxide, was injected as a slug of dilute caustic solution through a water-injection well. The natural wettability of the reservoir and the chemical requirements for reversing wettability were determined in the laboratory from contact-angle studies using field water and oil samples. Laboratory flood tests with a synthetic system had shown that reversing the wettability of an oil-wet consolidated core would lead to improved oil recovery. The field performance indicates that the mechanism by which increased oil retovery is obtained in the field is the same as that observed in the lahoratory. Laboratory studies indicate that higher ultimate recoveries and decreascd water-injection requirements result when the adjusting agent is added early in the life of the flood. However, a previously waterflooded area was intentiorzally chosen for the field test so that unambiguous conclusions could be made about the effects of chemical injection on wettability and the extent of oil-recovery improvernent afforded by wettability reversal. Forces existing at the fluid-solid and fluid-fluid interfaces in a porous medium have an important effect on oil recovery during a water flood. Modifying these interfacial forces in the reservoir to improve oil recovery has been the object of much research. Several papers&apos; &apos; have discussed the use of surfactants to lower the water-oil interfacial tension. This paper describes work concerned with improving oil recovery by modifying the forces at the fluid-solid interfaces— that is, by changing the preferential wettability. An earlier paper: presented evidence to show that some reservoir systems could be changed from preferentially oil-wet to preferentially water-wet by the action of simple chemicals added to the water to increase waterflood oil recovery. The process was attractive for further study because it used only inexpensive chemicals and because it called for a gross "wettability reversal" rather than for a precise adjustment. The present paper is an extension of the earlier work. Laboratory flow tests were made using treated water and oil in a consolidated core to determine the amount of additional oil recovery and the producing performance which might be expected from reversing the wettability of an oil-wet reservoir during a flood. A contact-anglc study using water and oil from the Muddy "J" sand of the Harrisburg field, Banner County, Neb., indicated the reservoir to be preferentially oil-wet and susceptible to wettability reversal through chemical injection. Based on these studies, a field trial of wettability-reversal water flooding was initiated. in the Harrisburg field. The primary purpose of the trial was to determine if wettability-reversal flooding would improve oil recovery in an actual field situation. The field trial was also an important test of the laboratory tecnique used to determine the effect on oil recovery of certain wettability manipulation, as well as a test of the contract-angle method of determining reservoir wetra-bility. The paper is presented in two parts. The first part covers the laboratory experiments leading up to the field trial, and the second part covers the field trial. PART I—LABORATORY STUDY The laboratory experiments consisted of displacement tests using an idealized fluid system and of a contact-angle study with crude and water from the Muddy "J" sand of the West Harrisburg Unit, Banner County, Neb. The displacement tests were designed to determine the effects of a particular wettability manipulation, oil-wet to water-wet, on oil displacement in a water flood. The contact-angle measurements were made to determine the natural wettability of the reservoir and to determine if beneficial wetting changes could be brought about by a chemical addition to the flood water. PROCEDURE DISPLACEMENT TESTS The displacement tests were perfornled in a consolidated sandstone core using a refined oil and water. The wetting properties of this system could