Search Documents

By R. H. Barham, G. W. Nabor

Linear aquifers, either limited or essentially infinite, may be encountered in reservoir engineering practice. In areas where faulting fixes reservoir boundaries, the fault block reservoir may have an aquifer of limited extent whose geometry is best approximated as linear. An infinite linear aquifer can occur as a regional feature whenever water movement through the aquifer member is constrained to one direction. Such constraints can arise from major faults, facies changes or pinchout of the member. Miller* pointed out that linear aquifers have received only meager attention in the past. He analyzed the performance of finite and infinite aquifers, developed working equations and curves, and presented examples. While Miller's curves may be used fairly easily, a separate one is required for each size of aquifer. In this paper, Miller's equations have been used as a starting point. By modifying them slightly, they can be reduced to a form which yields a single working curve, applicable to any size of aquifer. Thus. interpolation between curves is eliminated and accuracy is improved. Miller's results for finite aquifers covered only the boundary condition of no flow across the outer aquifer boundary. This paper also includes the case of constant pressure at the outer aquifer boundary.

Jan 1, 1965

By B. H. Caudle, B. M. Hickman, I. H. Silberberg

Secondary recovery projects often are not started in oil reservoirs until dictated by rising GOR's or declining oil production. Such circumstances require a well dispersed injection pattern to prevent the oil bank from bypassing the production wells. Many times this purpose is served by using the five-spot or line-drive patterns, in which half of the wells are made injection wells. A regular injection pattern is proposed that is based on the same well spacing as the five-spot but in which only a third of the wells need be converted to injection wells. Fluid flow model studies on this proposed pattern show that production histories will be as good as those for a five-spot pattern at mobility ratios of three and above. The model data are presented in a form that may be used to predict recoveries from both homogeneous and stratified formations. Introduction The injection of fluid into a petroleum reservoir usually is prompted by an actual, or predicted, loss in the productivity of the reservoir. Often, an oil reservoir is produced by pressure depletion until high GOR's dictate a change in policy. Fluid injection also causes increased ultimate recovery, but this usually is of less importance to the operator than is the maintenance of the production rate. It is a fact that the ultimate recovery in any particular fluid injection project is relatively independent of the number and spacing of the injection wells as long as the reservoir voidage through the production wells is balanced by the total injection rate. However, the response of a production well is directly affected by the distance it lies from an injection well. In places where the well productivities are declining, the injection wells should be spaced as evenly as possible among the producing wells ("dispersed" injection pattern). More particularly, if the production response depends upon the formation (and subsequent production) of an oil bank, each production well should be surrounded by injection wells. Otherwise, the oil productivity will increase very little as the oil bank moves past the well.'.' An evenly spaced injection pattern has the best chance of generating an oil kick in the production wells. If reservoir voidage and injection rates are to be equal, the ratio of injection wells to production wells must be equal to the ratio of the average injectivity to the average productivity (in reservoir volumes). This ratio is in strumental in choosing the injection pattern. Table 1 shows the production well to injection well ratio for the dispersed injection patterns that have been described in the literature.'-' All of these patterns will have comparable recoveries and response times (for the same basic distance between wells). The choice of pattern there fore depends upon the relative injectivity and productivity of the wells and upon the type of drilling pattern employed. Another important factor may be the ability of the proposed production wells to make up the lost productivity from those wells converted to injectors. These factors become more pronounced at the higher oil viscosities (low gravity oils). All other things being equal, the probable ratio of single well injectivity to single well productivity will be equal to the water-oil mobility ratio of the proposed flood. Thus, high mobility ratios (due to low gravity oil) should call for fewer injection wells — as long as these injection wells are "dispersed" to obtain more even pro duction response for the remaining wells. Secondary re covery projects in oil reservoirs seldom, if ever, will have mobility ratios requiring more injection wells than production wells. Thus, only the last five of the patterns listed in Table 1 will be of significant interest in oil recovery. Almost every development drilling program calls for the wells to be drilled along equally spaced rows and

Jan 1, 1969

By J. E. Berry

AII evaluation is made of the acoustic velocity log for measurement of formation porosity. Plots of field-observer1 velocities vs core-measured porosities of sandstones and limestotnes with inter intergranular porosity show that the velocity log gives a useful measure of porotity, in agreement with published data. A new method of plotting electrical and velocity log data makes it easier to recognize hydrocarbon-water coritnct~ from log darn and to make serni-quantitative estinlates of hydrocarbon saturation. Two examples of this interpretation technique are given. This paper also discusses the various factors that affect acoustic velocity and shows how corrections can be made for some of these factors. INTRODUCTION Although the continuous velocity log was developed primarily as an aid to seismic interpretation, it is finding widespread use for measurement of porosity. The purpose of this paper is to evaluate this use of the log. Results are presented for the velocity-porosity relationship in many sandstones and one limestone. Practical applications of porosity values derived from the velocity log to improve interpretation of electrical surveys are presented. Various parameters which affect velocity in porous media are considered. Porosity, composition, cementa-tion, pressure difference (overburden pressure minus fluid pressure), fluid saturation, and wellability are discussed. Some of the ideas are opinions based on general considerations of well logging problems while others are supported by published data and our own observations. Major conclusions from these studies arc given in the body of the paper and detailed discussions of some of the various factors affecting porosity are given in the Appendix. POROSITY Sandstones The primary factor affecting the velocity of sound in porous media is porosity. Wyllie, et al,1 have reported a time average equation for the relationship between velocity and porosity.* This has found considerable acceptance in the industry. Other relationships between velocity and porosity have been proposed, none of which are entirely satisfactory. It is doubtful that any one relationship will completely describe all situations. We have tried an empirical relationship based on correlating log velocities with core-measured porosities. Measured porosities were averaged over a depth interval equal to the spacing used in recording the velocity logs. The depth range for the information used was from 2,000 to 12,000 ft. Velocity measurements were normalized for the effects of depth and overburden pressure. A detailed discussion of the basis for this normalization is given in the Appendix. The change is usually insignificant and when working with individual formations it can be neglected. All velocity measurements were normalized to a depth of 10,000 ft. The velocity-porosity data for sandstones in Fig. 1 came from 20 wells scattered throughout the Mid-Continent and Gulf Coast areas. The scatter of the data points is attributed to three causes: 1. There is no theoretical basis for expecting an exact velocity vs porosity relationship and it is not possible to correct precisely for the effects of such factors as impurities, grain shape, manner of cementation, etc. Some of these factors are discussed in the Appendix. 2. Porosities measured by core analysis are not necessarily representative of what the velocity log "sees". 3. Inaccuracies appear in the velocity log reading. This is particularly true of velocities from one-receiver logging tools.*" A linear-time curve (time-average equation) is plotted in Fig. 1 along with the best straight line through the points. Velocities used for the linear-time curves were 19,500 ft/sec for the matrix velocity (zero porosity) and 5,500 ft/sec for fluid velocity. For clean sandstones with velocities greater than 10,000 ft/sec, there is little difference between the velocity-porosity and the time-porosity curves. Both curves are well within the scatter of data points. The data points in Fig. 1 were restricted primarily to clean sandstones, using the SP and/or gamma-ray logs as clues to the lithology. Where a lithologic de-

By G. E. Warner

This paper presents a review and analysis of a highly stratified Burbank sand waterflooding project in Osage County, Okla. Permeability values in this reservoir range from less than 0.1 md to nearly 3 darcys, and 90 percent of the permeability capacity is contained in approximately 26 percent of the reservoir volume. To date, the project has recovered 103 bbl/flood-acre-ft since waterflood initiation, which represents 50 percent of the total preflood cumulative and 12 percent of the original stock-tank oil in place (OSTOIP). The outstanding feature of this project has been the economical recovery of a sizable quantity of oil at water-oil ratios higher than those normally considered prohibitive. Nearly 70 percent of the recovery to date has been achieved at water-cut values greater than 95 percent. An ultimate waterflood recovery of 138 bbl/flood-acre-ft, or 67 percent of the preflood recovery and 16 percent of the OSTOIP, is predicted. The ultimate primary and secondary recovery will be 302 bbl/acre-ft and approximutely 40 percent of the OSTOIP. The actual performance is compared with the predicted performance by both the Stiles and the Dykstra-Parsons methods, and applicabiliry of these methods to this particular reservoir is evaluated. The performance of this project indicafes highly stratified reservoirs can be successfully and economically waterflooded when either reservoir or economic conditions preclude the use of mechanical and chemical means of mdbility control. Introduction The concept of the stratified reservoir has been used by reservoir engineers for many years to describe the vertical variations in permeability within a reservoir. Numerous techniques have been advanced to predict the expected waterflood performance of the stratified reservoir, and in most cases they prove reasonably accurate. The advisability of undertaking a prospective waterflooding project is determined, based on these predictions. This paper presents a review of a waterflood project in a highly stratified reservoir that normally would be considered a marginal prospect, at best, when evaluated with any prediction technique properly adjusted for reservoir conditions. This kind of project can be successful if the predictions are used to design a waterflooding system capable of economically recovering an oil cut of 1 to 2 percent. Experience gained in the early stages of the project's development was profitably applied in the later development of the entire reservoir and could prove beneficial to other projects with similar reservoir stratification or fracture characteristics. The East Burbank pool waterflooding operations have been developed in three separate stages, each of which has been a separate project: Mid-Burbank Unit, Flood A, and Stanley Waterflood. When referred to collectively, these projects will be called the Composite Project. The Skelly Oil Co.'s Barber lease — the only portion of the pool not included in the three projects — is operated as a separate facility with a cooperative injection-well agreement. It contains less than 80 productive acres and would not significantly affect the over-all performance picture.

Jan 1, 1969

By O. L. Culberson, J. J. McKetta

INTRODUCTION The equilibrium constants for methane and for water, and for ethane and water have been calculated from experimental data for the two binary systems.2,3,11,12 These constants are for the two-phase systems for the temperature range of 100" to 340°F and the pressure range of 200 to 10,000 psia. The constants for ethane are greater than those for methane. Equilibrium constants for a natural gas from a natural gas-water system6 are about the same as those for methane. Equilibrium constants for water in all three systems are very nearly the same over the greater part of the range of temperatures and pressures studied. In ideal solutions, the vapor-liquid equilibrium constants for a component should be the same in all such solutions at a given temperature and pressure.5 Few solutions are ideal, and consequently the use of generalized constants often leads to error in situations where they do not actually apply. In decidedly non-ideal systems, it is desirable to have equilibrium constant data for the specific system. Vaporization equilibrium constants may be calculated either by thermodynamic methods or from experimental vapor-liquid equilibrium data. In the thermodynamic calculation. the equilibrium constant K may be shown to be" K = fP,/fl'..........(1) where fPi is the fugacity of the pure component i at the given temperature and at the vapor pressure of the pure component, and fp is the fugacity of pure i in the gaseous state at the temperature and total pressure of the system. At best the thermodynamic calculation results in a K applicable to ideal solutions. Equilibrium constants calculated from Equation (1) have been found to have only limited agreement with experimental values, even though the solutions for which they were calculated approached ideality."'" If the more volatile component is at a temperature above its critical, the accuracy of the equilibrium constant for the component is made more doubtful by the use of the vapor pressure term under such conditions. In order to make the thermodynamic calculation on a component above its critical temperature, one must resort to one of the empirical methods of extrapolation of the vapor pressure to conditions above its critical point.13 Application of the thermodynamic equilibrium constant then to the methane-water or the ethane-water systems would be risky since the systems are non-ideal, and the methane and ethane are above their critical in most petroleum operations. ETHANE-WATER SYSTEM In the temperature range of from 100° to 340°F and for pressures up to 10,000 psia, the vapor-liquid equilibrium data 'References given at end of paper. Manuscript received in the office of the Petroleum Branch March 16, 1951.

Jan 1, 1951

By J. J. McKetta, O. L. Culberson

INTRODUCTION The equilibrium constants for methane and for water, and for ethane and water have been calculated from experimental data for the two binary systems.2,3,11,12 These constants are for the two-phase systems for the temperature range of 100" to 340°F and the pressure range of 200 to 10,000 psia. The constants for ethane are greater than those for methane. Equilibrium constants for a natural gas from a natural gas-water system6 are about the same as those for methane. Equilibrium constants for water in all three systems are very nearly the same over the greater part of the range of temperatures and pressures studied. In ideal solutions, the vapor-liquid equilibrium constants for a component should be the same in all such solutions at a given temperature and pressure.5 Few solutions are ideal, and consequently the use of generalized constants often leads to error in situations where they do not actually apply. In decidedly non-ideal systems, it is desirable to have equilibrium constant data for the specific system. Vaporization equilibrium constants may be calculated either by thermodynamic methods or from experimental vapor-liquid equilibrium data. In the thermodynamic calculation. the equilibrium constant K may be shown to be" K = fP,/fl'..........(1) where fPi is the fugacity of the pure component i at the given temperature and at the vapor pressure of the pure component, and fp is the fugacity of pure i in the gaseous state at the temperature and total pressure of the system. At best the thermodynamic calculation results in a K applicable to ideal solutions. Equilibrium constants calculated from Equation (1) have been found to have only limited agreement with experimental values, even though the solutions for which they were calculated approached ideality."'" If the more volatile component is at a temperature above its critical, the accuracy of the equilibrium constant for the component is made more doubtful by the use of the vapor pressure term under such conditions. In order to make the thermodynamic calculation on a component above its critical temperature, one must resort to one of the empirical methods of extrapolation of the vapor pressure to conditions above its critical point.13 Application of the thermodynamic equilibrium constant then to the methane-water or the ethane-water systems would be risky since the systems are non-ideal, and the methane and ethane are above their critical in most petroleum operations. ETHANE-WATER SYSTEM In the temperature range of from 100° to 340°F and for pressures up to 10,000 psia, the vapor-liquid equilibrium data 'References given at end of paper. Manuscript received in the office of the Petroleum Branch March 16, 1951.

Jan 1, 1951

By Harold B. Janzen, M. A. Nobles

Campbell and Martin1 described a new exploratory logging tool which was developed as the result of a new concept of the phenomenon that takes place when a porous formation containing mobile hydrocarbons and irreducible interstitial water is invaded by mud filtrate. They postulate that invasion will create a bank of interstitial water between the flushed zone and the undisturbed formation. This postulate was confirmed by Wyllie.' Thus, if R,, is greater than R, and mobile hydrocarbons are present, the resistivity profile with increasing radii from the borehole is high:low:high, as shown in Fig. 1. If mobile hydrocarbons are absent, the bank of interstitial water is indistinguishable from the undisturbed formation and the resistivity profile, when R is greater than R,, is high:low, as shown in Fig. 2. The zone containing the bank of interstitial water and residual oil was named the "low zone". Campbell and Martin conducted extensive electrolytic wedge and resistivity analog model studies to select a logging system that would reflect the presence or absence of the low zone. The system they selected was that of a family of lateral curves. The model studies indicated that the presence of a low zone could be detected in approximately 60 to 70 per cent of the producing sand reservoirs with a family of lateral curves provided that the sonde effectively occupies 90 per cent of the diameter of the hole. These studies were conducted with D, = 5d. The presence or absence of the low zone is best determined by a plot of the apparent resistivity vs spacing as shown in Figs. 3 and 4. A dip or inflection in the profile that repeats at successive depth intervals is considered diagnostic. To reduce the shunting effect of the hole and mud, a sonde of insulating material that is surrounded by a series of flexible rubber bumpers was designed. The bumper efficiency is considered to be the ratio of the bumper diameter to the hole diameter. Field tests conducted by Campbell and Martin on 10 wells with a tool that provided for station readings indicated that dips or inflections could be detected in 19 out of 25 known producing zones. INSTRUMENT DEVELOPMENT Halliburton developed a logging system to simultaneously and continuously record eight lateral curves. The instrumentation is described in the literature.' The original logs made

By A. Lubinski, C. R. Stewart, K. A. Blenkarn

Storage porosity has been considered one of the important pore geometry characteristics of heterogeneous-porosity limestones. Storage pores are only containers for fluids, in contrast to flow channel pores which both contain fluids and provide continuity for fluid flow. The concept of another geometry characteristic, "poro-striction", is presented as a second pertinent variable in describing limestone pore space. In simple terms, poro-striction is a measure of the flow resistance between storage and flow channel pores. Alternating-flow core-testing theory provides flow relationships which can be used to divide the pore space of heterogeneous-porosity media into flow channel and storage pores and to measure the "porostriction" of the latter. Experimental application of this theory to naturally occurring heterogeneous limestones shows that "porostriction" and the ratio of storage to flow channel pores can be estimated. Porostriction and porosity ratio are microscopic characteristics which should influence oil-recovery efficiencies during certain types of displacement processes. A knowledge of pore geometry should be valuable in designing or selecting the most effective oil-recovery process for heterogeneous limestones containing storage porosity. INTRODUCTION The primary and secondary recovery of oil from limestone (carbonate) reservoirs has been and will continue to be a major source of the world's supply of petroleum. Oil recovery from this type of formation is strongly influenced by the size, shape and arrangement of the pore spaces which hold the formation fluids. Therefore, it is expected that better knowledge of pore geometry would lead to design of more effective recovery processes. The concept of storage and flow channel pores has been presented as an important characteristic of limestone porosity.1 In this concept, storage pores are only containers for fluids, whereas flow channel pores are those pores which both contain and transmit the fluids. The work reported herein involves the concept of an additional pore characteristic called "porostriction". This characteristic can be visualized as the flow resistance between the storage and flow channel pores. This resistance is an important pore geometry characteristic because it determines how easily gas and oil can replenish fluids removed from the flow channel pores. Fatt and co-workers at the U. of California have studied the effect of storage porosity and a factor similar to porostriction on transient-flow response of models of porous media during a pressure build-up or drawdown test.55 These models consisted of porous sandstone cores with artificially simulated storage pores. It is the purpose of this paper to describe tests which permit the storage porosity and porostriction of actual reservoir core samples to be estimated. The tests used for this purpose involved not transient tests but, rather, an alternating-flow technique in which pressure and flow are varied according to a sine-wave relationship. This technique may offer advantages over transient methods. THE POROSTRICTION CONCEPT Before discussing alternating-flow (A.F.) tests in core samples, it is necessary to define the porostriction parameter. For this purpose attention is focused upon an elemental volume v representative of a porous medium having both flow channel pores and restricted-access storage pores. For such an elemental volume, it is considered that the storage porosity is combined to form a single pore and that the access into the storage porosity may be represented by a single capillary tube connecting the single pore to the flow channels of the elemental volume. Let the radius of the capillary tube be designated as Rc and its length as L,. When fluid flows into or out of the pore, the following equation (Poiseuille's law) relates flow rate q to pressure drop p, and the physical characteristics of the tube and the fluid of viscosity jx.

By O. C. Baptist

The Umiat oil field is in Naval Petroleum Reserve No. 4 between the Brooks Range and Arctic Ocean in far-northern Alaska. The Umiat anticline has been tested by 11 wells, six of which produced oil ; however, [lie productive capacity and recoverable reserves of the field are subject to considerable speculation because of unusual reservoir conditions and because several wells appear to have been .seriously damaged during drilling and completion. Oil is produced at depths of 275 to 1,100 ft; the depth to the bottom of the permanently frozen zone varies from about 800 to 1,100 ft, .so that most of the oil reserves are in the permafrost Reservoir pressures are estimated to range from 50 to 350 psi, increasing with depth, and the small amount of gas dissolved in the oil is the major source of energy for production. Laboratory tests were made on cores under simulated permafrost conditions to estimate oil recoverable by solution-gas expansion from low saturation pressures. The cores were also tested for clay content and susceptibility to productivity impaiment by swelling clays and increased water. content if exposed to fresh water. The results indicate that oil can be produced fronz reservoir rocks in the permafrost and that substantial amounts of oil can be produced from depletion-drive reservoirs by a pre.s.r~lrr drop of as little as 100 psi below the saturation pressure. Freezing of formation water reduces oil productivity much more than that due to increased oil viscosity: Failure of we1ls drilled with rtuter-base mud to produce is attributed to freezing of water in the urea immediately surrounding the wellbore. Swelling clays apparently contributed very little to the plugging of the wells. INTRODUCTION Naval Petroleum Reserve No. 4 lies between the Brooks Range and the Arctic Ocean in northern Alaska. The Umiat oil field is located in the southeastern part of the Reserve and is about 180 miles southeast of Point Barrow (the only permanent settlement in the Reserve and the primary supply point for drilling of the wells at Umiat). Eleven wells were drilled for the U. S. Department of the Navy, Office of Naval Petroleum and Oil Shale Reserves, between 1944 and 1953 to test the oil and gas possibilities of the Umiat anticline. Six of these wells produced oil in varying quantities and the best one pumped about 400 B/D.' Estimates of recoverable oil range from 30 to 100 million bbl. The main oil-producing zones are two marine sandstone beds in the Grandstand formation of Cretaceous age: these are referred to as the upper and lower sands. Good oil shows were found throughout the sand settions in the first three wells drilled on the structure, but the highest rate of oil production obtained on any 01 the many tests was about 24 BOPD. These first wells were drilled with conventional rotary methods using water-base mud; later wells were drilled either with cablc tools using brine or rotary tools using oil or oil-base mud. These experiments were successful as is shown by comparing the oil production from Well No. 2 with that from No. 5. These two wells are only 200 ft apart and are located at about the same elevation on the structure. Well No. 2. drilled with a rotary rig using water-base mud, was abandoned as a dry hole after all formation tests were negative. Well No. 5. drilled with cable tools and reamed with a rotary using oil, pumped 400 BOPD which was the maximum capacity of the pump and less than the capacity of the well. These field results indicated that the producing sands were extremely "water sensitive" and it was assumed that the cause of this sensitivity was the presence of swelling clays in the sands. Because of the very unusual reservoir conditions and the difficulties encountered in completing oil wells in the permafrost. the Navy asked the U. S. Bureau of Mines to make laboratory studies under simulated permafrost conditions to assist them in estimating the production potential of the field and the recoverable reserves. These tests were designed to determine the cause of the plugging of wells in the permafrost and to test oil recovery from frozen sand by solution-gas expansion with the oil gas-saturated at very low pressures. EXPERIMENTAL METHODS AND PROCEDURES Samples Analyzed Core samples were analyzed that represent the lower sand in Umiat Well No. 7, the upper sand in No. 3. and both the upper and lower sands in No. 9. These sands should be productive in all of the wells because of their location on the structure. Core samples from

By H. H. Rachford, D. S. Howard

In the withdrawal of fluids from tilted aquifers it is of value to be able to predict pressure patterns during the course of the pressure decline. As an example of this, in the displacement of fluids from oil reservoirs being produced by water-drive, prediction of the reservoir pressure is an important factor in fixing production rate. The prediction of pressure patterns in the depletion of an inclined aquifer has conventionally been made by solving the well known flow equations for the pressure Patterns in a horizontal aquifer, and correcting these patterm for the pressure difference due to the hydrostatic head associated with the vertical distance between a point in the horizontal aquifer and the corresponding point in the tilted aquifer. Inasmuch as the equation describing flow in tilted strata is different from that for horizontal flow. there was question as to the validity of the simple corrective procedure current1y in use. To examine the magnitude of error made using this procedure, an inclined aquifer was chosen for which the error was expected to he large, and the equation for flow in this tilted stratum was solved numerically. The results of the conventional procedure of correcting the solution of the equation for horizontal flow was compared with the solution of the correct equation by studying the difference between predicted pressures by the two methods. Three forms of the hydrostatic head correction were studied. The results for all three forms of the head correction are similar. The error in the approximate solution is quite dependent on the position of the horizontal aquifer with respect to the inclined one. The error is greatest at the producing end ad increases with amount of fluid withdrawn. The magnitude of the errors is not large if a good choice is made for the position of the horizontal datum. The most satisfactory position for the datum in all cases is near or slightly about the mid-point of the inclined stratum. Use of the best approximate method studied with the best datum tested yielded predicted pressures in error no greater than 7 psi. INTRODUCTION The movement of water by expansion from an aquifer to displace fluids from a contiguous oil reservoir was first recognized over 20 years ago.9 Since that time much work has been done to relate this mechanism to reservoir behavior mathematically. The result of this effort on the part of many workers has led to an understanding of underground fluid movement which has enabled us to follow and predict the behavior of the water-drive reservoir. One reason that good progress has been made in the application of mathematics to the study of pressure patterns and water movements in aquifers has been that the differential equation which describes this process is of a relatively simple form. The equation for the isothermal, horizontal, viscous flow of water through homogeneous sandstone is identical with that for the diffusion of heat through homogeneous solids, and this equation and means of obtaining solutions for it with many types of boundary conditions have been studied extensively. That this work has been highly successful is clearly emphasized by the large number of solutions developed for problems of interest both by methods of formal mathemat-

Jan 1, 1957

By J. J. Arps

This paper reviews the methods currently in use for estimating primary oil reserves and discusses the principles on which these methods are based. Particular emphasis is placed on how these methods change with the type of information available during the life cycle of an oil property. This paper contains various novel estimating methods and shortcuts heretofore unpublished. INTRODUCTION Estimating oil reserves is one of the most important phases of the work of a petroleum engineer since the solutions to the problems he deals with usually depend on a comparison of the estimated cost in terms of dollars, with the anticipated result in terms of barrels of oil. His recommendations to management regarding the best course of action are therefore normally based on the most favorable balance between these two. Specific engineering problems which require such a knowledge of recoverable oil reserves and a projection of future rates are: [a.] the exploitation and development of an oil reservoir; [b.] the construction of gasoline plants, pipelines and refineries; [c.] the division of ownership in unitized projects; [d.] the price to be paid in case of a sale or purchase of an oil property, and the magnitude of the loan which it will support; [e.] the proper depreciation rate for the investment in oil properties; and [f.] evaluation of the results of an exploration program. This discussion will be confined to the various methods and tools which are currently in use for estimating oil reserves to be obtained during the primary phase of an oil-producing reservoir and for a projection of the future production rates. Reserves which may be obtained by secondary recovery methods or fluid injection programs and gas and gas condensate reserves will not be discussed in this paper. Unfortunately, reliable oil reserve figures are most urgently needed during the early stages when only a minimum amount of information is available. Management's interest in the oil recovery from a property—aside from its use for accounting purposes—usually declines when the property approaches its economic limit, just at the time when the reliability of the estimates is at its best. To illustrate this general idea a chart is presented (Fig. 1) showing the three periods in the life of an imaginary oil property. Time is shown on the horizontal axis, while the cumulative production and estimated ultimate recovery are plotted vertically. No particular units are used and this schematical chart is not to scale. During the first period, before any wells are drilled on a property, any estimates will of necessity be of a very general nature, based on experience from similar pools or wells in the same area, and usually expressed in barrels per acre. This will

Jan 1, 1957

By A. Timur

Fluid flow properties of porous media have long been of interest in such varied disciplines as geology, geophysics, soil mechanics, and chemical, civil, and mechanical engineering. This interest has resulted in numerous models of porous media that have been proposed, tested, and found to be useful under, at best, only special circumstances. A critical review of most of these models is given by Scheidegger. 1 The application of each new tool to the study of porous media has helped to test some of the existing theories and to form new ideas based on the new parameters being measured. In the application of nuclear magnetic resonance (NMR) techniques to the study of the properties of fluids in porous media, the theoretical studies of Korringa et al.2 resulted in a model for the relaxation of spin polarization of protons in a hydrogenous fluid in a pore of a solid. Seevers3,4 verified this model, established a method for measuring the surface-to-volume-ratio distribution in a porous medium, and proposed a technique for determining specific permeability in sandstones. In the present investigation, the uses of NMR methods are extended and a simple model of porous media is developed from an analysis of pulsed NMR measurements. In this model, the pore spaces of a porous medium are divided into three groups on the basis of their surface-to-volume ratios. To evaluate predictions of the NMR model, laboratory measurements of spin-lattice relaxation time, porosity, permeability, and residual (irreducible) water saturation were conducted on more than 150 sandstone samples obtained from four oil fields in North America. Applications are discussed for estimating the volume of movable fluid, ?p, and the specific permeability, k, of sandstone samples. The former parameter, which can be considered as producible porosity, is expressed by ?p = [1 -- (Swr/100)]? , .... (1) where + is porosity in percent of bulk volume and Swr is the residual (irreducible) water saturation in percent of pore volume at a capillary pressure of 50 Psi. Fast, nondestructive laboratory methods are described for determining porosity, movable fluid, and permeability of sandstone samples from measurements with a pulsed NMR (spin-echo) apparatus. Procedures are given for adapting these methods to samples of irregular shape and small size. Theory General In the Korringa, Seevers, and Torrey2 (KST) model the observed decrease in the spin-lattice relaxation time, T1, of protons of a hydrogenous liquid contained in the pore spaces of a porous solid is attributed to an increase in the correlation time for the random motion of the water molecules and to the presence of paramagnetic centers at the liquid-solid interface. This KST model predicts the observed relaxation time of a hydrogenous liquid in a pore of a solid through

Jan 1, 1970

By D. O. Wyble

Five water-driver fields, for which pressure predictions had been made within a few years after discovery, were re-studied when four to five years of additional pressure and production data became available. It was fortnd that forecasts of the pressure in water-dri~fe reservoirs may be subject to considerable unccrtainty. Predictions car be erroneous when there is insufficient presswe history to define the cllaracteristics of the aquifer, and when there are errors in the pressure, production or geologic data. INTRODUCTION The efficient operation of an oil or gas reservoir requires that the important characteristics of the reservoir be determined soon after its discovery. It is especially desirable to evaluate the pressure-production relationship quickly, for this relationship provides a basis for deducing the type of production mechanism that is operative, and for making pressure predictions. The pressure predictions, in turn, provide a basis for calculating efficient rates of production, determining whether water or gas should be injected, etc. The importance of pressure prediction is such that forecasts should be made as accurate as possible, and the limits of error well defined. Historically, pressure forecasts for some fields have proved very accurate, while for others the forecasts were so erroneous that they resulted in the adoption of wrong, or at best, unnecessary production practices. Pressure predictions made for water-drive fields are especially subject to error. Accurate forecasting for these fields requires that much be known not only about the reservoir proper but about the adjacent water-filled rock as well. The water-containing rock, or aquifer, may extend for many miles, and often there is only meager data regarding its thickness, permeability and porosity. Much of the information must be inferred from the previous pressure behavior of the reservoir. A brief pressure history may provide little information about the characteristics of the aquifer. Hence, pressure predictions made early in the life of water-drive fields can be incorrect. The work described in this paper was done to determine the degree of reliability which can be expected of pressure predictions made within a few years after water-drive fields are discovered. Five fields, for which the original pressure predictions had been based on less than seven years' history, were re-examined with the aid of four to five years of additional pressure data. Where the predictions were inaccurate, the reasons for failure were ascertained. The studies were accomplished by use of the electric reservoir analyzer.' The analyzer was first adjusted to repeat the original pressure-production match using the same values for the various parameters as were used in the original study. Then, the actual production rates experienced during the time elapsed since the initial study were set into the analyzer, and the pressure behavior indicated by the machine for this period was compared with that which actually OCcurred. Where the pressure profile traced out by the analyzer deviated from the observed field behavior, the analyzer was readjusted to obtain a better over-all match. The old and the new solutions were compared to determine the reasons for any differences. COMPARISON OF PREDICTED RESERVOIR PRESSURE BEHAVIOR WITH ACTUAL PERFORMANCE Field A The original pressure-production match is entirely adequate as evidenced by the pressure match of the field's history subsequent to Jan., 1952 (see Fig. 1)'.

By J. R. Kyte, L. A. Rapoport

The material balance equation for partiai or complete 1:trter-drivc reser1,oir.s been re-arranged to include a pressure irrferferencr ternr. This pressure interfernce term was ohtained from the theory established hy Morinda. A .successful progranl deissed to solve the resrclritlg c,qrration in a few hours on a nleelitr .sizeti, digital compuier. No prior knowdelge of the properties of the eqrrifer is required in the solution. 7' ri~rl calculations areatle using various assumed con!binuiions of value of the constant expretsing the qquifer properties correspoinds to the best fit ohtained between observed and calculated field history. Details of the method of calculaion are illus trated by an example. INTKCDUCTICN A method 01' calculating the material balance equation for a partial water-drive reservoir completed in an aquifer has been described in the literature." Pressure interference i'rc,m other pools iocated in the same nquiucr is not taken into account in this solution. When present, the pres- sure interference effect must be included in order to make a complete performance analysis. A correction for this effect has been developed using the theory cstahlished by Mor-tada.' Manual calculation of the material balance for water-drive reservoirs is a formidable task. Such manual ca1culation when interference eflects are includcd would be impossible to completc in any reasonable time period. However, a medium-sized stored program, digltal conlputer can he used satisfactorily. The method is described in the following sections. THECKETICAL CONSIDEKATIOhS The material balance equation for a reservoir producing by partial or total water drive, as proposed by Van Everdingen, This relationship is linear for Z/a from which thc slopc, Z, and the intercept, N, can be obtained by the least-squares method from the data points. It is not necessary lo know the aquifer properties to uktain oplimum values tor 2nd IV. Various assumptions are made tor At,,: where 11 is production time, seconds (oilfield between the equal successive time intervals into which the history is divided. t,, h It, times the number of equal time periods from i.e.. the total dimen-sionless time interval from It is noted that the dinlensionless time interval, At,,, involves all the relevant aquifer properties. Accordingly, any assunlption of is an assumption of the combined aquifer properties which pertain to the solution of the equation. Deviations of the data points from a straight line arc calculated for each assumption of and the value of A, correspond ing to the minimunl deviation is chosen as the optimum solution. This method is described in Ref. 1 and 2. Thc material balance solution is usually made with an initial assumption that the reservoir being studied is located in an infinite aquifer from which there are no other withdrawals. Deviation of the points by amounts greater than can be attributed to inaccuracies in field data may be the result of: (1) boundary effects of a limited aquifer, (2) pressure interference from othcr pools completed in the .sarne aquifer, and (3) a combination of these effects. Differentiation of these effects can be made only on the basis of geologic data. Pressure interference from other fields will result in a pressure drop at the original oil-water boundary having an interference component, The measured instantaneous pressure drops at the oil-water

By B. L. Landrum, P. B. Crawford

The rise of a new laboratory model for studying tran-sient fluid flow problems, is described. The theory of he model is based on the analogy between the equa-ions which describe the flow of compressible fluids in porous media and the conduction of heal m solids. The "thermal model" is used to study the Iransient Pressure distributions during fluid displacement in five Standard well patterns. The studies indicate quantita-ively how die pressures change with lane at any poit in the reservoir. It was found that, in general, the equipressure contours are more radial about the wells early in the program than after steady-state conditions are established. In addition two studies of fractured five-spot patterns are presented to illustrate the opplicability of the model for a wide range of reservoir conditions. INTRODUCTION In reviewing the oil recovery techniques involving fluid displacement it was found that a very considerable quantity, if not all of the fluid within the pattern, mav be displaced (during a time-dependent or transient period rather than the usually applied steady-state conditions. Since the oil is displaced during a transient period, it was believed desirable to study the nature of the transients for several of the commonly used fluid displacement patterns. A literature review indicated fairly complete studies of the streamlines and isopotentials for steady-state: yet few data were available on the transient displacement process. Several simulated steady-state multiwell displacement programs were carefully mapped by Foster and Lodge on an electrical conducting sheet as early as 1875, hut the more recent work of Muskat provides the best source for steady-state phenomena having direct petroleum production applications. In order to study the transient phenomena a new lab-oratory model was developed which is based on the similarity of the equations that describe the flow of fluids and the conduction of heal. The equations describing heat flow were developed early in the nineteenth century. Regorous solutions tor most transient problems are in gcneral complex and usually limited to fairly restricted houndary conditions.2 the thermal model may be used to obtain approximate solutions for many transient problems without some of the restrictions required for analytical solution. The application of the analogy between heat conduc-. tion in solids. electrical flow and fluid flow, in a porous media has long been recognised. Techniques originally developed for the solution of heat conduction prol?lems have in many instances proved useful in ohtaining solutions to fluid flow prohlcms. Ir is not surpri5ing then that a heat conduction or thermal model might bf usefull in sludying transient fluid flow. As an example ol the problems suitable for the study with the thermal model, transient pressure distributions during fluid displacement programs are studied for the five-spot, line-drive, staggered line-diibe. seven-spot 211 two fractured patterns. A model was constructed of an iron plate of geometry similar to the reservoir patterns under study. The principle of mirror images as suggested by Bryan was employed in shaping the model. A temperature differential was obtained by soldering two sections of copper tubing to the plate at points analogous to the well location in the field. Hot and cold water was passed through the tubing. The temperature distribution, analogous to the pressure distribution in the field. was obtained a a function of time from thermocouplcs distributed over the plate.

By N. Mungan, F. W. Smith, J. L. Thompson

Adsorption of polymers and transport, rheology and oil recovery efficiency of their solutions were studied in the laboratory to evaluate the use of polymers in waterflood-ing. While a tenfold mobility reduction was obtained with polymer concentrations as low as 0.05 per cent by weight, the mobility reduction depended on the type of polymer, molecular weight, salinity and pH of water, crude oil and capillary properties of the porous media. Choice of a suitable polymer and a workable concentration will have to be tailored for each application. Little reduction in the residual oil saturation can be expected from polymer flooding. Improvement in the volumetric sweep efficiency is possible hut the extent of the improvement can best be evaluated by properly designed field testing. Some aspects of the field use of polymer floods are discussed. INTRODUCTION Waterflooding is a simple, inexpensive secondary recovery method and is being used widely. Innumerable laboratory studies have been made to unravel the fundamentals of the displacement of oil by water and to find the ways of most efficient oil recovery. These studies and a great many field case histories have revealed that the prime cause of poor oil recovery is the inefficient and incomplete sweep of reservoir volume by the injected water. Sweep efficiency is affected by many factors of which the mobility ratio is an important one. Mobility ratio M is defined here as the ratio of water to oil mobilities: M = (k»/y,r)/(k,JJJJ........(1) In Eq. I, the permeabilities are the effective permeabilities and depend on fluid saturations and, hence, change during the different depletion stages in a flood. A wide practice is to use the effective water permeability at residual oil saturation and the effective oil permeability at interstitial water saturation in Eq. 1. If the mobility ratio is greater than one, the mobility ratio is unfavorable and water, being more mobile than oil, would finger through the oil zone resulting in poor oil recovery efficiency. If the mobility ratio is favorable (one or less) the displacement of oil by water occurs more or less in a pistonlike fashion. In some waterfloods. the mobility ratio is unfavorable and any additives by which the mobility of water can be decreased would favor more efficient oil recovery. The thing to bear in mind, however, is whether or not the improvement in oil recovery is sufficient to more than pay out the cost of the additives needed. For example, materials like sugars, alcohols and glycerine reduce water mobility by thickening the water, but the cost of material requirement precludes any field application. For an additive to be useful in water-flooding, it must bring about a large reduction in water mobility at low concentrations; it must be adsorbed only negligibly; and it must not completely plug up the formation. Some synthetic organic polymers have shown promise of meeting these requirements and have been used in the field.'-W owever, no in-depth studies of the rheological, adsorption and oil displacement characteristics of polymer solutions have been reported. The present work is a study of these properties. EXPERIMENTAL In this work, concentrations are given on a weight per volume basis; 0.5 per cent concentration means 0.5 gm of polymer is dissolved in enough water (or NaCl solution) to make 100 ml. A bactericide, usually 0.1 to 0.2 per cent by volume of 38 per cent formaldehyde solution, was used in the polymer solutions. The NaCl solution was 30,000 ppm. Some properties of the polymers studied are given in Table I. Physical properties of all cores used are in Table 2. Flow behavior of polymer solutions was studied by three consecutive flow tests in cores. First, water (or brine) was injected at constant rate of about 1 ft/D to obtain the water mobility. Then, filtered polymer solution (prepared in water or brine) was flowed through the core. Since the rate was constant, increase in the pressure drop across the core reflects decrease in the mobility. Finally, the core was flooded with water (or brine) to study recovery of mobility. The Alundum cores which were used in

Jan 1, 1967

By H. A. Deans

A three-parameter mathematical model for one-dimensional flow in porous media is developed. The objective of the model is to predict accurately the longitudinal dispersion associated with the flow of either gases or liquids over a wide range of Reynolds number. A qualitative analysis of the model parameters is given. Published experimental pulse-response curves are compared with results predicted by the model. Several general types of problems are given for which the model can be used. INTRODUCTION The problem of multicomponent single-phase flow through porous media is encountered in the study of petroleum reservoirs, gas chromatographic and ion-exchange columns, industrial fixed-bed con-tacters and elsewhere. These particular examples span the wide range of flow conditions possible; Reynolds numbers of less than l0-4 are not unusual in oil-production problems, while values in excess of 103 are common in large fixed-bed operations. The study of flow-dependent transport phenomena is complicated both by the changes in the character of the flow over this range, and by the irregularity of the flow boundaries inherent to porous media. Dispersion is one of the important phenomena known to depend fundamentally on flow conditions as well as on fluid and medium properties. As used in this paper, the term "dispersion" refers to the observed mixing of fluid elements of different composition which occurs in flow systems. The actual mechanism may be one or more of a number listed below. Only dispersion in the direction of the mean flow (referred to as axial or longitudinal dispersion, or mixing) is considered here, although lateral dispersion arises as part of certain coupled mechanisms. MECHANISMS FOR LONGITUDINAL DISPERSION A number of distinct mechanisms are known3,5,7,8 to contribute to the phenomenon of longitudinal dispersion. The more important of these are as follows. 1. Molecular diffusion in the flow direction. 2. Turbulent (cell) mixing. 3. Lateral transport processes coupled with velocity and/or residence time distributions, including: (a) "Taylor" diffusion caused by the interaction of velocity profiles in individual voids with lateral molecular diffusion; (b) separation and remixing or interdiffusion of streams having different velocities around particles; and (c) the couplipg of gross velocity profiles, caused by viscous instability or inhomogeneous porosity, with lateral dispersion. 4. Finite mass-transfer rate between a porous matrix and the flowing phase, and finite diffusion rate inside elements of the porous matrix. It is a well documented fact that these mechanisms, whether acting individually or in combinations, produce essentially similar integral effects. The observed dispersion can be described approximately by solutions of the diffusion equation with a properly chosen dispersion coefficient. Necessary conditions for equivalence between dispersion and diffusion are discussed by Klinkenberg and Sjenitzer.7 Apparently these conditions are not satisfied exactly in many real systems. "Best fit" solutions of the diffusion equation often depart appreciably from measured breakthrough curves for the usual pulse- or step-forced experimental systems. In particular, pulse response curves usually show some degree of asymmetry and "tailing" which cannot be reproduced by the one-parameter diffusion model. The three-parameter model proposed herein attempts to produce these effects in a realistic manner. Asymmetry and tailing are predicted for certain sets of parameter values, while for other sets the new model agrees with the diffusion model. VELOCITY DEPENDENCE OF DISPERSION Molecular diffusion apparently controls at sufficiently low Reynolds number in both gas and liquid systems.*l3 At Re > 20, turbulent cell mixing is dominant for gases flowing at low pressure; the data of McHenry and Wilhelm8 show that the mixing cell is approximately 1 particle diameter long in random sphere packs. on the other hand, the data

By L. T. Stanley

A method is described whereby the behavior of the front is approximated when gas is injected into a thick reservoir having reasonably homogeneous properties. The method is applied to the Abqaiq field of Saudi Arabia. The producing member is the Arab-D, a high-permeability calcarenite which has no continuous betlding planes. Field performance indicates the injected gas to be partially miscible with the reservoir oil because gas volume factors experienced are lower than those predicted by the gas law. Eqlrilibril~m calculations confirnz ihis and also show that the volume factor may vary appreciably, depending on the relative quantities of oil and gas in contact in the reservoir; consequently, volutnetric balances are made for the purpose of determining the behavior of the gas volume factor during injection history. The effect of gravity segregation on position and shape of the gas front is approximated by dividing the producing section into three zones vertically and tracing the movement of the average interface in each zone, with counterflow of oil and gas taking place betweell zones behind the front. This results in a gas front that advances more rapidly along the top of [he formation than along the base, producing an "umbrella" effect in cross section that becomes more pronounced as the front progresses. The calculations cover the gas-injection period of Abqaiq field history through 1958, and the position of the gas front is plotted at one-year intervals. Finally, some comparisotls between the calcrlated fronial position and the gas-oil contact as Measured by neutron logs are made. INTRODUCTION Fig. 1 is a structure contour map of the Abqaiq field. For analytical purpcses, the field has been divided into two areas designated A and B, as shown on the map. Gas injection into Area A was begun in early 1954 into two wells located near the crest of the structure on the north-south axis of the anticline. During the period covered by this paper, the average injection rate has been 122 MMscf/D with a maximum sustained rate of 160 MMscf/D. The reservoir pressure has been maintained at 2,450 psig since start of gas injection. The reservoir originally was undersaturated but, prior to pressure maintenance, the pressure was drawn below the bubble point in most of the Area A section, creating an in-place free-gas saturation of 8 per cent at the crest of the structure. The average critical gas saturation is 13 per cent established by laboratcry core tests. That the entire reservoir remained below critical gas saturation prior to gas injection is borne out by the fact that none of the crestal wells have exhibited high gas-oil ratios during production history. The Area A reservoir is a carbonate section about 200-ft thick having average porosity and permeability of the order of 22 per cent and 500 md, respectively. The productive horizon is the Arab-D member of Jurassic age and is encountered at an average depth of

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

First response to large-scale water flooding in the fractured very low permeability Spraberry sand has led to a new unique cyclic operation. Capacity water injection is used to restore reservoir pressure. This is followed by marly months production without water irzjection and the cycle repeated. Expansion of the oil, rock and water during pressure decline expels part of the fluids but capillary forces hold much of the injected water in the rock. At least with reservoir pressure restored and with partial water flood development, field performance has proved this cyclic operation is capable of producing oil from the nzatrix rock at least 50 per cent faster and with lower water percentage than is imbibition of water at stable reservoir pressure. INTRODUCTION The Spraberry Field of West Texas presents unusual problems for both primary production and water flooding. Extensive interconnected vertical fractures in the fractional-md sandstone permitted recovery of oil on 160-acre well spacing, but they made capillary end effects dominant. Primary recovery by solution gas drive is less than 10 per cent of oil in place. The concept of displacement of oil from the sand matrix by capillary irnbibition of water has led to field techniques which promise greatly increased oil recovery. Free exchange of laboratory research, reservoir information and results of field pilot tests among the various companies has been very important in development of this technology. Five units covering a total of 170,000 acres have been formed for water flooding, and 10 other areas covering an additional 175,500 acres are in various stages of unitization. Part of the Driver Unit reaching fillup first has demonstrated very unusual waterflood behavior and indicated numerous operating problems that will develop within and among the various units. SPRABERRY ROCK AND PRIMARY PERFORMANCE The Spraberry, discovered in February, 1949, is a 1,000-ft section of sandstones, shales and limestones with two main oil productive members: a 10-15 ft sand near the top and a 10-15 ft sand near the base. In part of the field some thinner intermediate sands are oil productive, and others are water bearing. All sands have permeabilities of 1 md or less and porosities of 8-15 per cent. Ordinary core analysis and electric and radiation logs are ineffective in differentiating between oil productive and nonprcductive sands. Sands capable of containing producible oil are best identified by mercury injection capillary pressure measurement and, in some cases, by core water saturation. About 3,500 wells have been drilled in the 500,000-acre trend. Vertical fractures were observed in practically all Spraberry cores. Continuity and interconnection of fractures were confirmed by pressure interference among wells during early development.' Major fractures trend northeast-southwest as indicated by oriented cores and confirmed by five fluid injection tests, by analysis of the pressure transients observed during development,''' and by three interference tests in the Driver Unit Water Flood reported herein. Fracture spac- ing probably averages inches to a few feet. Spraberry wells typically produced 100-400 BOPD initially after hydrauLic fracture treatments. By 1962 oil production had declined to an average of 12 bbl/well/day, near the economic Limits of operation. Reservoir pressure had declined from 2,300 psi initially in the Upper Spraberry and 2,500 psi in the Lower Spraberry to 500-1,000 psi. Partial closing of the fractures with declining reservoir pressure is believed to be the cause of such low oil production rates at these relatively high reservoir pressures. Cumulative recovery of 208 million bbl of oil is 80 to 90 per cent of that recoverable by primary means. Performance of the entire reservoir is summarized in Fig. 1. IMBIBITION WATER FLOODING By 1952 reservoir performance indicated low primary recoveries. Most engineers, expecting serious channeling of injected fluids through the fractures, held little hope for secondary recovery. With its extensive background of research on the fundamentals of fluid flow within reservoir rocks, Atlantic's Research and Development Division on short notice in 1952 conceived that displacement of oil by capillary imbibition of water into the rock might significantly increase Spraberry recovery. Laboratory data reported by Brownscombe and Dyes scaled to probable reservoir conditions showed potential waterflood recovery equal to or greater than primary recovery with a 10-15 year flood life.= A pilot test using three 40-acre injection wells, one central producing well and 18 surrounding observation wells demonstrated technical feasibility of the process. Injection of 1.5 million bbl of water from November 1952 to August 1955 proved water entered the rock and displaced oil

By J. N. Sicking, F. H. Brinkman

A new method for obtaining equilibrium vaporization ratios (K-values) for reservoir fluids has been developed and tested. By application of the method, complex experimental measurements of liquid and vapor phase compositions are eliminated. This simplified technique reduces the cost of experimental equilibrium-ratio data for reservoir studies of condensates and volatile crude-oil systems. The method is designed for systems of constant composition and, therefore, is best suited for depletion studies where compositional changes at high pressures are minor. The basic data required, in addition to the composition of the initial reservoir fluid, are the relative vapor-liquid volumes and densities at reservoir temperature and variom reservoir pressures. Tests demonstrated that the method predicts equilibrium ratios accurately for condensates. A single test on a crude oil was not conclusive; further testing will be necessary before the accuracy of the method can be determined for crude-oil systems. In addition to determining equilibrium ratios, the calculation method provides information on the physical properties of the "plus" component in the vapor and liquid phases. The "plus" component is that mixture of components heavier than the least volatile fraction analyzed. This information is useful in studies of both natural depletion and cycling operations for condensate reservoirs where the heptanes-plus component in the gas phase is produced from the reservoir. INTRODUCTION As more volatile oil and condensate reservoirs are found, the use of phase behavior techniques to predict their performance is increasing in importance. These techniques have long been used for condensate fields and have more recently been applied to crude-oil fields containing oils of medium-to-high volatility. In these phase behavior methods, equilibrium ratios (K-values) are used to predict compositional changes in the reservoir fluids—thereby accounting for the recoverable oil that exists in the gas phase. The reliability of the predictions depends to a large extent on the equilibrium ratios used. These values must be obtained for each component for the entire pressure range being investigated. Unfortunately, because of the complex nature of hydrocarbon mixtures, accurate K-values are hard to obtain. The equilibrium ratios for a particular component will vary not only with the temperature and pressure, but also with over-all composition of the system. The importance of composition is quite critical at elevated pressures, but becomes negligible at pressures below about 300 psia. Therefore, because most phase behavior problems involve the high-pressure region, each fluid system becomes a special case. Experimental programs to determine characteristic K-values are quite difficult and time-consuming. Thus, it is often necessary to resort to approximations of the K-value data. Charts giving K-values for various mixtures and classes of mixtures are available in the literature. However, there are two major difficulties in using them: (1) the K-values of the "plus" component (that mixture of components heavier than the last one analyzed) must be obtained by extrapolation from the K-values of the other components; and (2) the K-values obtained must finally be adjusted by trial and error to agree with observed volumetric data. To eliminate these difficulties, a new method of determining equilibrium ratios was developed. Briefly, after the composition of the system as a whole has been analyzed, the method uses empirical correlations and the gross fluid properties of the system (relative vapor-liquid volumes and densities) to calculate K-values. Because the calculative procedure is long, it is best solved on a digital computer. About one hour of machine time on an IBM 650 computer is required to develop a K-chart for the fluid being examined. DEVELOPMENT OF THE METHOD OF OBTAINING K-VALUES Equilibrium ratios are defined as the ratio of the mole fraction of a component in the vapor phase to its mole fraction in the liquid phase. This statement is expressed in Eq. 1. A typical plot of equilibrium ratios for a particular system is shown in Fig. 1. It should be noted that, at pressures near the saturation pressure, the K-values appear to converge to a common point. This apparent convergence point is called the convergence pressure and is a characteristic of the system involved. Various empirical correlations of K-values have been noted. It has been observed that an isothermal plot of