Search Documents

By H. T. Kennedy, J. E. Little

An empirical equation for the prediction of the viscosity of several pure paraffin hydrocarbons and nitrogen is presented. It involves temperature, pressure and six constants of the material, and it applies reliably to both liquids and gases. The equation is similar in form to van der Waal's equation of state. For the paraffin hydrocarbons methane through r-hexane and nitrogen, an average absolute deviation of 1.9 percent was obtained on 1,006 data points described in the literature by 14 authors. When this equation is extended to complex, liquid hydrocarbon mixtures, a correlation was obtained with an average absolute deviation of 9.9 percent. INTRODUCTION Equations describing the flow of gas and liquid through porous media contain the viscosity coefficient of the fluid. If other pertinent variables remain constant, the volume rate of flow is inversely proportional to this coefficient. In dealing with condensate fluids and volatile oils, however, the compositional effects resulting from changing pressure materially affect the viscosity. The effect of compositional changes also may be significant in certain secondary recovery or pressure maintenance processes, notably miscible displacement or gas injection. Early attempts to describe the performance of reservoirs utilized a volumetric material balance method wherein gas and liquid in the reservoir were identified as produced gas and liquid at the surface. This method of analysis proved adequate for reservoirs at moderate temperature and pressure that contained gas with moderately low amounts of condensable materials. The volumetric material balance procedures for "black oil" reservoirs leave much to be desired when applied to condensate and volatile oil reservoirs because phase behavior and compositional changes are relatively more important in these cases. The alternative is a compositional material balance, which in turn, requires a correlation of properties of the reservoir fluid with composition. This paper supplies this correlation in regard to viscosity, for reservoir crude oils. REVIEW OF LITERATURE The literature contains many empirical equations describing the effects of composition, temperature and pressure on the viscosities of pure liquids and binary liquid mixtures. However, the applicability of a majority of these equations is limited to very low pressures and to a small number of systems. Most of them, when applied to complex hydrocarbon systems, are of little value. The lack of utility of the majority of equations results from the fact that they were developed to show the separate effect of temperature, pressure or composition on viscosity, but not to predict the viscosity as a function of all three variables. And with the few exceptions noted below, they were developed to apply to much simpler systems than hydrocarbon mixtures. Eakin1 developed an equation for the viscosity of pure light hydrocarbons applicable in the gas, liquid and dense fluid regions, and in its most general form is written as µ - µo =A(ebP - e-cpk) The pure component viscosity at atmospheric pressure and temperature T, µo, was calculated using the Sutherland equation: µo = BT3/2/T + S The constants in Eakin's equation permitted the term e-CPk to become quite small at high densities (above 0.215 gm/cc), and the general equation would simplify to

Jan 1, 1969

By G. E. Perry, R. D. Seba

The steam soak process is the most widely applied and most successful thermal supplemental recovery process in use today. This process, which consists of injection of steam in various quantities into a well that is subsequently used as a producer, can be repeatedly applied to the same well. The treatment, if successful, will increase the production rate of the well for a given period. For this reason, the steam soak process should be considered primarily as an acceleration method rather than as a method for producing oil that could not be recovered by primary means. However, application of this process to reservoirs containing viscous crude may produce oil that otherwise would have zero present value or would be produced at such low rates by primary methods that it would not be profitable. A major necessity for the successful application of the steam soak process is a source of natural energy. This source of energy may exist in one or more of the following forms: expansion of fluids, compaction, gravity drainage. This study was undertaken to determine the effect of gravity drainage on steam soak production and recovery and to develop a method for predicting production rates, ultimate recovery, and heat requirements when applying this process to a reservoir in which the primary driving mechanism is gravity drainage. (For additional discussion of some of the factors that affect recovery of viscous oils by the steam soak process, see Ref. 5.) Model The model developed for analyzing the performance of the steam soak process in gravity drainage reservoirs is a modified form of that developed by Matthews and Lefkovits.1,2 The model used in this study consists of two concentric cylindrical volumes of productive formation with cap and base rock of impermeable material. The inner cylinder, extending from the radius of the wellbore to some predetermined radius, is referred to as the hot zone; and the outer cylinder extending from the exterior boundary of the hot zone to the outer drainage boundary will be referred to as the partially heated zone. Fig. 1 is a schematic cross-section of the model, indicating its thermal properties. Horizontally the reservoir is divided into two regions, the upper region having a high gas saturation and residual liquid saturation, and the lower region containing a high liquid saturation. The boundary between these two regions is referred to as the free surface, which drops and changes shape as fluid is produced from the reservoir. The model was used to analyze the steam soak process by applying the following assumptions: 1. The horizontal formation is composed of one or more homogeneous noninterconnected layers containing incompressible fluids. Each layer is treated separately. 2. The hot zone is heated by the initial soak and is maintained at the steam temperature corresponding to atmospheric pressure by resoaking. 3. The height of the free surface at the producing well is at all times zero (producing well is pumped off). 4. The position of the free surface at the start of

Jan 1, 1970

By J. Jones-Parra, R. S. Reytor

The porosities of fractured limestone reservoirs can be divided into two broad types in accordance with their effects on fluid distribution and fluid flow. In the coarse porosity, gravity segregation lakes place freely and the resistance to fluid flow is very small. In the fine porosity there is no segregation and a high resistance to flow, and it has relative permeability characteristics similar to tight sandstones. By analyzing the affect of the two porosities it is concluded that in some cases to recover the maximum amount of oil it is necessary to remove large quantities of gas from the reservoir by producing at high gas-oil ratios. In this manner the fine porosity is drained of its oil and the gas-oil contact drops slowly permitting higher production rates from the oil leg. Since this conclusion is contrary to widely accepted principles of conservation, a mathematical model (Fig. I) was constructed to duplicate the conditions desired. The behavior of the model indicates that under certain conditions it is possible to recover more oil by producing at high gas-oil ratios than by production at low gas-oil ratios, and that the rate of production is affected more by gas shut-offs than by the decrease in pressure (Figs. 2, 3, 4 and 5). INTRODUCTION It has been noted that certain limestone reservoirs behave in such a way as to indicate that there are two distinct types of porous spaces available for fluid storage and fluid flow. Regardless of the nature of these porous spaces the distinguishing characteristic is that in one type of porous space there is definite gravity segregation between the oil and the gas; while in the other, the gas evolved tends to remain distributed equally throughout the reservoir. For the sake of simplicity the porosity in which gravity segregation takes place at a fast enough rate to affect the behavior will be called "coarse po- rosity", while that in which this is not so will be termed "fine porosity". The behavior of a reservoir of this type would be different than that of a sandstone reservoir, especially if most of the oil is in the fine porosity. It has been observed that fluid flow to the wells takes place mostly through the coarse porosity. The large fractures which

By R. H. Adams, A. M. Khan

The Huntington Beach oil field is on the Pacific Coast, southeast of the City of Los Angeles (Fig. 1). It lies both onshore and offshore. Most of the onshore section is operated by Standard Oil Co. of California, WOI. The Tar reservoir is of upper Pliocene age, about 200 ft thick, and consists principally of fine-grained, unconsolidated, somewhat shaly sand, with thin inter-bedded shale stringers throughout (Fig. 2). The reservoir is a broad, highly faulted, elongated anticline, with an areal extent of about 370 acres. Most of the faults bisecting the structure are normal, trend north-south, and are the main cause of reservoir discontinuity (Figs. 3 and 4). From 1964 to Dec., 1967, 259 steam stimulations were performed on a total of 125 Tar sand wells. Eighty-six new wells were drilled and steam stimulated. Of the total, 96 have been stimulated twice, 34 three times and 3 four times. During this same period, a two-injector steam drive pilot has been in progress. The response obtained from this project has been compared with the response to cyclic stimulation. This paper presents an analysis of the field observations for the two methods of thermal recovery. History of Development Drilling to the Tar pool commenced in 1924 with an initial well production of 40 to 50 B/D of 12o API gravity oil. Complete primary development of the pool has been spread over several years. Approxi- mately 40 wells were drilled during the first phase of development, which was completed in 1945. The subsequent primary development started in 1960 with initial well productions of 30 to 40 B/D oil. Thirty-seven wells were completed in the Tar pool during this phase. For the most part, slotted, gravel flow packed liners have been used in the wells in this pool. A cyclic steam stimulation pilot was initiated in March, 1964. Based on the very successful stimulation in the pilot area, the cyclic stimulation project was expanded to include all of the Tar pool wells producing at the end of 1965 (Fig. 5). Also, a development drilling program was begun at this time. To date, 86 new wells have been drilled under this program. The number of wells drilled each year in this secondary phase of development is also shown in Fig. 5. New Well Completion Practices Most of the new wells have been completed to produce from both the Upper and Lower Tar zones. Generally, a casing stub is cemented from the top of Lower Tar to the base of clean Upper Tar sands to exclude intermediate water-bearing sands. The water string is then cemented from the top of Upper Tar to about 1,500 ft with neat cement and 30 percent silica flour. The annulus from 1,500 ft to surface is then filled with a petroleum-based gel through a high-temperature port collar. Once again, a history of sand-free production dictated the use of 60-mesh gravel flow packed liners. The Upper and Lower Tar

Jan 1, 1970

By H. J. de Haan, J. van Lookeren

In view of their size and favorable characteristics, Shell's heavy-oil fields on the eastern coast of Lake Maracaibo, Venezuela — known as the "Bolivar Coast" — are obvious candidates for the application of thermal recovery methods. During the last decade various methods have been tested, such as steam drive,', steam soak and in-situ combustion.3 Steam soak was discovered as a promising production method rather accidentally in 1969, during early steam drive testing in the Mene Grande Tar Sands1 When steam erupted at the surface due to breakdown of the overburden, the injection wells were back-flowed to relieve the reservoir pressure. This resulted in high oil production rates, all the more impressive because the reservoir is unproducible by primary means. It was concluded that injection of limited amounts of steam might be a very effective method for stimulation of heavy-oil wells. To obtain information that would be more widely applicable, testing of the process was continued in the eastern part of the Tia Juana field (Fig. 1). Here even more favorable results were obtained than in Mene Grande. This created sufficient confidence to prompt a large-scale project, which was considered essential for an adequate evaluation of ultimate recovery, steam requirements, mechanical problems and economic prospects. To obtain an early answer on the ultimate recovery a fairly depleted area was selected in the western part of the same field ("A" in Fig. 1). The first steam was injected in March, 1964. This paper describes the performance until the end of 1966. The production mechanism and well inflow performance are discussed in some detail. In view of the favorable results obtained so far, steam soak operations have expanded considerably since the start of this project. Presently, five other projects are under way, of which two are about the same size as the one under discussion, with further extensions being programmed. Description of the Project Reservoir Data Shell's part of the Main Tia Juana reservoir has a proven area of 8,800 acres and an average net oil sand thickness of 130 ft. It contained an initial quantity of 2,750 million bbl of stock tank oil (STOIIP), with an API gravity ranging from 10 to 15 API and a viscosity of 100 to 10,000 cp (at 122F). As a geological unit it is known locally as the Lower Lagun-illas member of the Lagunillas formation, which is early Miocene in age. The sand bodies are believed to be distributary channel sands deposited in a lower coastal plain. Most of the channel sands have lateral and some vertical communication. The field was developed mainly in the period from 1936 to 1950. The ultimate primary recovery is estimated at 18.1 percent of STOIIP, of which 12.5 percent STOIIP had been produced by the end of 1963 from some 540 wells. By then the average pressure had decreased from an initial value of 960 to 550

Jan 1, 1970

By T. K. Perkins, L. R. Kern, R. E. Wyant

Many factors influence the results of hydraulic fracturing to stimulate well productivity. Most of these factors have been studied and their effects discussed in the literature. However, the movement of sand through the crack during fracturing has received little consideration. Exactly where the sand goes in the fracture, how it gets there, and how this might influence results have not been studied, judging from published reports. To get an answer to some of these questions we have studied sand movement through a simulated crack segment in the laboratory. Although the work is preliminary, we have answered some questions about sand movement in a semiquantitative sense. We are reporting our results at this time in the hope of stimulating interest and further study in this important aspect of fracturing. LABORATORY EXPERIMENTS Fig. 1 is a schematic diagram of equipment used in this study. A segment of a vertical crack 0.25-in. wide was formed between two Plexiglas plates, and movement of sand in fluid through the crack was observed under various conditions. Sand and fluid were injected into the crack at constant rates. The top and bottom of the crack were closed so that the fluid flowed horizontally. Because of gravity most of the sand settled to the bottom before moving very far through the crack. If the fluid velocity exceeded a certain critical value, then all the sand injected was washed on through the crack even though it had settled to the bottom. When the fluid velocity was less than this value when sand injection was started, most of the sand which settled in the crack remained in place until enough of the area open to fluid flow was blocked by the settled sand to increase the fluid velocity to the critical value.. At this point all injected sand moved on through the crack even though it settled to the lower portion of the moving fluid section. All sand previously settled remained in place. When fluid velocity was increased at this paint by increasing fluid injection rate (without also increasing sand injection rate), some of the previously settled sand was washed out until enough area was open to flow that the velocity dropped to the critical value. Conversely, when the fluid velocity was decreased by decreasing the fluid injection rate (without changing the rate of sand injection), additional sand would settle until the critical velocity was again attained. In short, the critical velocity is an equilibrium value which the system will automatically seek, if possible. The equilibrium (or critical) velocity was measured under several conditions. The results of these measurements are shown in Pig. 2 where the equilibrium velocity is shown as a function of the amount of sand injected per minute. The points falling in the shaded band shown in Fig. 2 were obtained with water. The higher viscosities were obtained by adding CMC (carboxymethylcellu-lose). Data indicate that the equilibrium velocity is not very sensitive to sand size or transport fluid viscosity, or to sand injection rate (except at very low rates). Points connected by the solid CUrve in Fig. 2 were obtained with water and 20- to 40-mesh steel shot. The equilibrium velocity is apparently fairly sensitive to difference in density between particle and fluid. Only one run was made with a gelled fluid and that was gelled water at zero sand injection rate as shown in Fig. 2. (Measurements of equilibrium velocity at zero sand injection rate were made by filling the crack

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

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

By A. S. Odeh

During the last few years, several authors1-6 have advanced mathematical solutions, both exact and approximate, to the pressure behavior of a well producing from a region bounded by a circular discontinuity. The most recent solutions were reported by Carter, and by Bixel and van Poollen. Varter solved the appropriate equations analytically. Bixel and van Poollen reported a numerical solution and provided correlation charts for the interpretation of flow tests; their method uses an overlay technique and is valid for wide ranges of mobility and storage capacity ratios. The mobility ratio is defined as the ratio between k/p of the two zones, and the capacity ratio is defined as the ratio between $c of the two zones. We have developed independently an analytic solution, programmed for the CDC 1604, that agreed with Carter's.'1 The same differential equation, assumptions, boundary, and initial conditions were used. We ran 30 drawdowns and 30 buildups of a unit capacity ratio and of varied mobility ratios. We also varied the radial distances to the discontinuity, and the radii of drainage. We examined the results and obtained correlations that were simple to apply and useful in flow test interpretation. The results and the correlations reported here in graphical form allow the determination of the mobility of the two zones and the radial distance to the discontinuity. They eliminate the need to fit field data to a curve as is required in the overlay technique.6 Thus, the advantages of the present correlations are that they require less effort and they provide a unique answer. They also provide guidelines for selecting the parts of the flow test plot to be used in the interpretations. Moreover, the results obtained from the programmed analytic solution define the time intervals where the approximate solution reported by Hurstl and discussed by Larkin2 are applicable. Normal buildup and drawdown plots resulted in an early straight line, a long transition zone and a late (second) straight line. However, the second straight line will not occur if the drainage radius is not at least 10 times the radial distance to the discontinuity. The mobility of the imer and outer zones can be calculated from the slopes of the first and second straight lines, respectively, by the normal methods of the flow test analysis. We found that the radial distance to the discontinuity may be determined from the abscissa value for the point of intersection of the extensions of the two straight lines. Results similar to some of these were reported by Adams et aL7 The approximate solution compared favorably with the exact one during the first and second straight-line periods. However, during the transition zone, deviation occurred between the two solutions. Discussion of Figures Figs. 1 through 7 summarize the results of this investigation. On Figs. 3 through 7 there appears the equation of each respective plot. In this paper the mobility ratio M is defined as klp2/kzpl, where subscripts 1 and 2 refer respectively to the inner and outer zones, and dimensionless time = t, = 6.32 X In the investigation, M ranged from 0.25 to 50;

Jan 1, 1970

By L. F. Elkins

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

Jan 1, 1970

By H. S. Price, J. C. Cavendish, R. S. Varga

A numerical formulation of bigh-order accuracy, based on variational methods, is proposed for the solution of multidimensional diffusion-convection-type equations. Accurate solutions are obtained without the difficulties that standard finite difference approximations present. In addition, tests show that accurate solutions of a one-dimensional problem can be obtained in the neighborhood of a sharp front without the need for a large number of calculations for the entire region of interest. Results using these variational methods are compared with several standard finite difference approximations and with a technique based on the method of characteristics. The variational methods are shown to yield higher accuracies in less computer time. Finally, it is indicated how one can use these attractive features of the variational methods for solving miscible displacement problems in two dimensions. INTRODUCTION The problem of finding suitable numerical approximations for equations describing the transport of heat or mass by diffusion and convection simultaneously has been of interest for some time. Equations of this type, which will be called diffusionconvection equations, arise in describing many diverse physical processes. Of particular interest here is the equation describing the process by which one miscible liquid displaces another liquid in a one-dimensional porous medium. The behavior of such a system is described by the following parabolic partial differential equation:

Jan 1, 1969

By J. C. Martin

Equations for three-phase, three-dimensional, compressible flow (including capillarity) are reduced to two-dimensional relations by a partial integration. This reduction allows three-dimensional flow problems to be treated with mathematics for only two spatial dimensions. The results can be used to formulate flow equations for two-dimensional reservoir simulators in which the effects of capillarity and fluid segregation in the third dimension are represented. Such reservoir simulators would retain many of the aduantages of two-dimensional simulators while simulating three-dimensional effects. The principal restriction of the method is that the thickness of the reservoir should be small, compared to the distance across the reservoir. INTRODUCTION In recent years, computers have been used to calculate the performances of many reservoirs. Most of the detailed calculations, however, are based on finite difference solutions of the flow equations, and present day computers are seldom able to handle a sufficient number of cells to produce entirely satisfactory solutions, even for reservoirs represented by two-dimensional arrays of cells. The situation becomes much worse when one wishes to approximate the reservoir by a three-dimensional array. A great saving in computation or a more detailed solution can be obtained for many reservoirs by using the partial integration of the equations of flow presented in this paper. The integration reduces the three-dimensional equations to two-dimensional relations; and, for studies of two-dimensional flows in vertical cross-sections, the equations are reduced to one-dimensional relations. Most reservoir performance calculations currently are based on one- or two-dimensional flow re1ations.l-3 In some cases flow in the third dimension is approximated by assuming a particular type of vertical saturation distribution, such as gravity segregation.l.2. The relations developed in this paper approach those for segregated flow as the capillary pressures approach zero, and they approach those for uniformly distributed saturations as the capillary pressures are increased. For this analysis, the ratio of the reservoir's thickness to the maximum distance across it must be small. The capillary pressures between the oil and water should also be small compared to the maximum pressure difference in the reservoir. This requirement is met by most reservoirs. It is assumed that the capillary-pressure curves are well defined, whether or not hysteresis effects are included. Also, the reservoir must have sufficient vertical permeability to allow the fluids to segregate. The results presented here provide a firm theoretical foundation to Coats' et aL5 assumption of vertical equilibrium and extend the relations to three-phase flow. Coats' assumption of vertical equilibrium, which he verified by calculations and experiment, is developed here mathematically from basic flow equations. DISCUSSION SATURATION AND PRESSURE DISTRIBUTIONS Appendix A presents a mathematical analysis of fluid flow in reservoirs where the ratio of thickness to maximum distance across the reservoir is small. The results indicate: (1) that the fluids along any line perpendicular to such a reservoir's upper surface are in static capillary equilibrium (vertical equilibrium); (2) that, to a first approximation, the fluid pressures and properties are functions of only areal position in the reservoir and time; and (3) that hydrostatic pressure gradients exist along any line perpendicular to the reservoir's upper surface. The results might be expected after studying several physical considerations. First, no flow is allowed normal to the upper and lower reservoir boundaries, which are relatively close together.

Jan 1, 1969

By H. J. de Haan, L. Schenk

Scope for Thermal Recovery in Shell's Heavy Oil Fields in Venezuela The main heavy oil reservoirs on the East coast of Lake Maracaibo (Fig. I), known as "Bolivar Coast", initially contained some 20 billion bbl of oil in place with gravities in the range of 10 to 15° API. The current total production rate is about 400,000 B/D. Since the recovery to date is on an average only 12.5 percent of the initial oil in place, these fields offer a vast potential for secondary recovery methods. The reservoirs are characterized by moderate depth (generally 1,000 to 3,000 ft; maximum 5,000 ft), good formation properties such as net oil sand thickness (50 to 300 ft), high porosity (30 to 40 percent) permeability (1 to 3 darcies) and oil saturation (initially about 80 percent) and high oil viscosity (100 to 10,000 cp in situ), all of which are favorable for the application of thermal recovery processes. Since 1957, steam drive, steam "soak" and underground combustion have been tested or are being tested in this area.1-3 This paper deals with a major steam drive test in the Tia Juana field. The test was commenced in Sept., 1961, on the basis of encouraging results of laboratory experiments and pilot field tests, which will be discussed briefly. Early Laboratory Investigation of the Steam Drive Process In 1956 a series of model experiments was carried out in the Koninklijke/Shell Laboratorium, Amster- dam, to investigate the displacement of heavy oil by steam.' The prototype studied was a horizontal, un-consolidated heavy oil reservoir, subjected to a linear steam drive. The model, representing a slab from this reservoir along the main flow direction, consisted of a steel tube filled with oil-saturated sand, provided with layers of rock at the top and bottom to simulate conductive heat losses to cap and base rock. Lateral heat losses, perpendicular to the main flow direction, were minimized by insulating the sides of the model. Heat flow and viscosity and gravity forces were scaled by taking the steam injection rate and sand permeability inversely proportional to the geometric reduction factor of the model, and using materials (oil, steam, cap rock, etc.) identical with those of the prototype. The main feature of the process was a frontal displacement mechanism that was observed within wide ranges of injection pressure and rate, oil viscosity, initial oil saturation and sand permeability. The high stability of this front can be attributed to the high volumetric steam flow rates, combined with the reduction in oil viscosity. High recoveries were obtained due to the low residual oil saturation of about 15 percent in the zone swept by steam. The latter value pertains to the specific test conditions (heavy oil, uncon-solidated sand) and may be lower for light oil, due to the effects of steam distillation. The applicability of this experimental approach is limited, as it does not take into account the fact that

Jan 1, 1970

By R. W. Parsons

The over-all apparent single-phase permeability of fracture-rock systems was studied using two different two-dimensional models. In a strict sense the results are applicable only to these models, yet may give some insight into real reseruoir behavior. The regular fracture-matrix model can have any number of regular fracture sets. All fractures within each set have the same orientation, width and spacing. The gross behavior is amlytically equivalent to that of an anisotropic permeable medium. In the heterogeneous fracture systems, individual fracture conductivities vary within a two-dinensioml network pattern. A computer program calculates the overall permeability after the conductivities are randomly placed in the network. The calculated permeability depends upon the distribution of these conductivities and the particular fracture pattern into which they are placed. Although detailed information on actual reseruoir fracture systems is not available, this study shows how some of the possible variables may affect obserued gross reservoir characteristics. INTRODUCTION Characteristics of reservoir rocks are determined by direct (core analysis) or indirect (logging and production tests) examination. Both methods, singly and together, are widely used in reservoir evaluations. For fractured reservoirs a problem exists in relating characteristics of cores to the in situ reservoir properties. Knowing fracture characteristics can be of value in determining production mechanisms. Because of the limitations of observation, the intimate details of a fractured rock system will never be known. One approach in analyzing fractures is to study the behavior of some conceptually simple models. The two models analyzed contain features which may be realistic, yet are mathematically tractable. Flow is considered in two dimensions only, with the plane of interest at right angles with the fractures. This concept represents the horizontal flow in a reservoir system with the fractures oriented vertically. Only single-phase, laminar flow is considered. The regular fracture-rock system permeability is calculable from a simple analytical expression. The heterogeneous system requires a computer program for calculating overall permeability.

Jan 1, 1967

By C. F. Weinaug, R. W. Farley, J. F. Wolfe

A rapid, accurate method for predicting the dew points of gas condensate systems and their subsequent normal and retrograde phase behavior with pressure decline has been developed. The method predicts dewpoint pressure, volume percent liquid, condensate content of produced gas, and the liquid and gas compositions at pressures below the dewpoint pressure. Use of a modified form of the BWR equation of state, which was established by regression calculations on experimental data from several gas condensates, provides depletion behavior predictions that are significantly improved over predictions using K-values based on convergence pressures and Standing's correlation for density. INTRODUCTION Prediction of complex hydrocarbon-mixture thermodynamics is fundamental to reservoir and gas process engineering. In particular, accurate descriptions of phase behavior are needed to plan production operations for gas-condensate reservoirs. Although analysis of reservoir behavior during producing and cycling operations can be made for simple reservoir models, extensions of these studies to more complex reservoir systems are still difficult, both with regard to numerical techniques and the phase behavior of the fluids. This report presents a mathematical method for handling the complex phase relationships encountered in retrograde condensates; it represents an advance in reservoir behavior technology in that it will aid the development of closed-form mathematical representations of total production systems. A typical phase diagram for a complex hydrocarbon mixture is shown in Fig. 1. When the reservoir temperature T, is higher than the critical temperature Tc, and lower than the critical-condensation temperature Tct' the mixture is a gas condensate. Decreasing the pressure on a condensate while maintaining nearly constant temperature will result in two phases forming at the dew point. Further decrease in pressure will cause additional liquids to be formed. This phenomenon is called retrograde condensation. As pressure continues to decrease, normal phase behavior regions are encountered (liquid vaporizes as the pressure decreases) until the low-pressure dew point is encountered. When a gas-condensate reservoir is produced, gas is withdrawn and the volume of fluids in the reservoir remains nearly constant as the pressure declines. Heavier components tend to accumulate in the liquid phase as a result of retrograde condensation, while lighter components tend to concentrate in the vapor phase. Thus, as the gas is produced the composition of the gas and the condensate

Jan 1, 1970

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

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

By R. C. Earlougher

Jan 1, 1970

By B. S. Gottfried

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

Jan 1, 1967

By P. J. Closmann

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

Jan 1, 1969

By Robert G. Nisle

The classic theory of pressure build-up in shut-in oil wells as developed by Horner and van Everdingen is based on two-dimensional radial symmetry in the well-reservoir system. Such symmetry does not exist in the case of a well which parfiially penetrates the producing formation. As a result of this lack of symmetry the use of the classic theory in such cases become questionable. In this paper the mathematical theory has been extended to include the case of partially penetrating wells. Numerical solutions illustrating the up-dication of the equations are presented. The effect of partial penetration on pressure build-up is shown by a comparison of synthetic pressure build-up curves derived from the numerical solution of the equation for partially penetrating wells for various degrees of penetration. It is shown that partial penetration is detectable from the characteristic shape of the pressure build-up curve and that formation productivity may be calculated from the pressure buildup data in a manner identical to that described by the classic theory. INTRO DUCTION The use of pressure build-up data on shut-in oil wells is a well-established technique for the measurement of reservoir productivity. The work of Homer' and van Everdingen' is well known. Since the publication of their pioneering work, others"' have elaborated and extended the technique. These authors have dealt exclusively with the case of wells that completely penetrate he producing formation. Two-dimensional radial symmetry is achieved, thereby, and the problem is greatly simplified. This ideal situation is seldom encountered in practice. But, in spite of this, the technique is often employed anyway. The question, therefore, arises: how much confidence can be placed in the results of such calculations when the data are taken on wells that only partially penetrate the producing formation? It is the purpose of this work to answer this question. Only the single-phase fluid case will be considered. van Everdigen' has presented a rigorous solution of the two-dimensional case, but because of mathematical complexity, this solution has not proved to be particularly useful. Horner' utilized Kelvin's point source solution in two dimensions and showed that it is a sufficiently good approximation for oil reservoir studies. The geometry of the well and the producing formation is illustrated in Fig. l. Geometrically the formation may be described as an infinite slab. The infinite slab is defined as the solid (the porous medium in this case) bounded by two parallel planes. Thus, the slab extends to very great distances in both the positive and negative directions of two of the coordinate axes, but is finite in

By H. T. Kennedy, S. M. Avasthi

An equation developed for gaseous hydrocarbon mixtures predicts molal volumes with an average absolute deviation of 0.73 percent when applied to 264 natural gas and condensate systems including 2,043 PVT points. Another equation developed for liquid hydrocarbon mixtures predicts molal volumes with an average absolute deviation of 1.12 percent when applied to 346 crude oil systems including 1,759 PVT points. Both equations require composition of the mixture to be expressed as mole fraction of methane through heptanes-plus, hydrogen sulfide, nitrogen and carbon dioxide, together with the characteristics of the heptanes-plus fraction in addition to the temperature and pressure. The equations cover wide ranges of the variables involved, and their accuracy is considerably better than that of other available methods. The equations were differentiated to allow calculation of the coefficients of isothermal compressibility and isobaric thermal expansion. (In this paper the coefficient of isothermal com-pressibility and the coefficient of isobaric thermal expansion will be expressed as compressibility and thermal expansion coefficient, respectively.) Equations to calculate these quantities are presented. INTRODUCTION Calculations of reservoir performance for petroleum reservoirs require accurate knowledge of the volumetric behavior of hydrocarbon mixtures, both liquid and gaseous. Compressibilities are required in transient fluid flow problems, and thermal expansion coefficients are important in thermal methods of production. An accurate laboratory investigation of the PVT behavior of each reservoir fluid encountered would be costly and time consuming. For this reason various correlations for predicting fluid properties have been developed and recorded in recent literature. Correlations have been presented in the form of graphs, tables and equations. Since an increasing number of studies are being conducted with the aid of electronic computers, recent efforts have been directed toward development of correla-tions suitable for computer programming. Application of computers permits the use of more complex correlations which otherwise are not feasible. Moreover, methods for predicting reservoir performance, particularly those based on the compositional material balance, depend upon the capability of accurately expressing the molal volumes and other fluid properties as functions of pressure, temperature and composition. The coefficient of isothermal compressibility c is defined by

Jan 1, 1969