Search Documents

By Kenneth A. Grace

Each and every mining operation from a simple sand-and-gravel pit to the largest integrated mine-mill-concentrator-smelter-refining complex has, as its foundation, a quantity of extractable natural resource the reserves. Without reserves, there is no mine, and any concentrating or extracting plant, no matter how sophisticated, is simply so much scrap metal. And, of course, without resources there are no reserves. It is remarkable, therefore, especially in this age of careful planning, high quality engineering, and state-of-the-art financial analysis, that such a fundamental and vital item is so little understood and the subject- of so much confusion. Time and time again it has been the experience of the author, when called in to a mining operation that has run into trouble in reaching projected production targets, to find that the designed-for tons were not reaching the concentrator/extraction plant, and/or the grade of such tonnage was substantially below the forecast and design grade of the plant. The mill blames the mine for delivering below standard tons and/or grade and the mining engineer blames the geologist for overly-optimistic assessment or incorrect interpretation of the deposit. The geologist blames the mining and metallurgy division for wrongly interpreting his data and conclusions, or simply fouling-up a fine orebody.

Jan 1, 1983

By A. L. Barnes, A. M. Rowe

A heat transfer study was made of hot gas injection into oil shale through wells interconnected by vertical fractures. This analysis involved the simultaneous numerical solution of a nonlinear, second-order partial differential equation that describes two-dimensional conduction heat transfer in oil shale and a nonlinear first-order partial differential equation that describes convection heat transfer in the fractures. Three nonlinear, temperature-dependent coefficients were used in this work; they are thermal conductivity, thermal capacity and retorting endothermic heat losses of oil shale. Vertical fractures were considered to be of finite height. Although vertical conduction heat transfer was not considered, an estimate of the error resulting from this limitation was made. The effects of injected gas temperature, injection rate, system geometry, cyclic injection and time upon retorting efficiency were investigated. Results from this study show that the rate of retorting oil shale is a direct function of both injection temperature and rate, and the theoretical producing air-oil ratio:(AOR) is an inverse function of temperature. Retorting rates are constant until "breakthrough" of the 700 F isotherm at the producing well, assuming constant injection parameters. Retorting rates for bounded systems are higher than the analogous unbounded systems and likewise AOR's are less. The use of an alternating injection-soak routine with high injection rates is less efficient than continuous injection at lower rates. These results indicate that injection temperatures on the order of 2000F or greater may give theoretical AOR's in the economic range. INTRODUCTION Over half of the known oil shale reserves are located in the U.S., and most of them lie in the Piceance Creek basin of Western Colorado. The Colorado oil shale outcrops on the edges of the Piceance Creek Basin. At the outcrops the shale beds are relatively thin, from 25 to 50 ft thick. In the center of the basin the oil shale is as great as 2,000 ft thick and is covered with 1,000 ft of overburden. It has been estimated that there are over 1,000 billion bbl of oil in shales having an oil content over 15 gal/ton in this basin. Oil shale does not contain free oil but an organic matter called kerogen. Kerogen yields petroleum hydrocarbons by destructive distillation. It must be heated to approximately 700F, at which temperature it decomposes into shale oil, gases and coke. The U.S. Bureau of Mines and, more recently, oil companies have conducted considerable research on surface retorting methods to economically recover oil from this shale. Another approach to exploit the oil shale deposits, in particular that portion having 1,000 ft of overburden, is to retort the oil shale in place and produce the liquid and gaseous hydrocarbons through wells drilled into the shale. Some research has been done on this approach.l,2 There are several variations to the in situ retorting approach. These variations fall into one of two groups, depending upon the geometry of the system: (1) retorting in a highly fractured or broken up matrix; (2) retorting from single fractures between production and injection wells. The latter is the group studied. Several investigators,3-6 using various assumptions, have studied flow of heat through horizontal systems. The objective of this work was to make a heat transfer study of in situ retorting oil shale by hot gas injection through wells interconnected by single vertical fractures of finite height. The oil shale thermal conductivity, thermal capacity and retorting endothermic heat losses were considered to be functions of temperature. A knowledge of the

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 F. G. Miller, T. D. Mueller, L. B. Davidson

A mathematical model of reservoir behavior during one cycle of the intermittent steam injection process is presented in this paper. A cycle of the process is considered to be composed of three stages: injection, condensation and production. The model concentrates on behavior during the production stage. The influence of several parameters on production improvement is described, using the results obtained with a numerical equivalent of the model. Reservoir system parameters considered are thickness, reservoir thermal conductivity, oil viscosity, solution gas-oil ratio and the change in oil effective permeability with temperature. Operating parameters investigated are heat injection rate, steam temperature and cumulative injection time. In the face of considerable uncertainty, the model appears to give a fair prediction of the behavior observed during an actual test. This makes the model a useful means of selecting projects. It is shown that the model may also be utilized to select an optimum set of operating conditions for a particular project. Results of this study show the extent to which production improvement depends on each of the parameters investigated and give an indication of the range of system conditions within which the process would be effective. INTRODUCTION In recent years the intermittent steam injection or steam soak process of oil recovery has become a popular method of improving recovery in reservoirs bearing high-viscosity hydrocarbons. Production improvement due to the process results from heating the reservoir near the wellbore. The process is applied on an individual well basis and is executed in a series of cycles, each consisting of three stages: injection, condensation (the soaking period) and production. During the injection stage, high-pressure steam is injected into the reservoir surrounding the well-bore. This study assumes that injection is carried out over the entire vertical section of the reservoir rock that can be exposed to fluid injection. The condition of steam entering the reservoir is considered to fall within the following limits: (1) steam pressures less than 1,000 psig, (2) steam temperatures between 300 and 700F and (3) steam quality less than 90 percent. The last restriction is due to equipment limitations. Normally, injection re maintained for less than 3 weeks. After termination of injection, the condensation stage begins. The well is shut in for several days, supposedly to make certain that the injected steam condenses. This condensation would insure the retention of injected heat. The production stage is initiated when the well is opened again. Reservoir response during this stage determines the effectiveness of the process (response is defined as the oil production rates, the water production rates and the wellbore temperatures observed during the entire production stage). The immediate response of the reservoir is sometimes quite good. Sets of case histories by Longl and by Owens and Suter2 give some typical figures for response. Duration of the production stage is equal to the time needed for the oil production rate to decline to some predetermined value. This usually requires 4 to 6 months. Termination of the production stage coincides with the start of the next cycle. The process is continued, cycle by cycle, until it becomes unprofitable to continue. One simple model2 and two more sophisticated models3,4 of the process have been presented. These are based on fundamentally different assumptions than those used in this investigation. Although the three models might be modified to fit the assumptions used in this study, it is felt that greater rigor would result by developing a

By C. S. Matthews, P. Hazebroek, F. Brons

A method has been developed for calculating the average pressure in a bounded reservoir. The reservoir is first divided into the individual drainage volumes of each well, by using the criterion that at steady state each individual drainage volume is proportional to a well's production rate. The average pressure in each drainage volume is then calculated by a method developed in the report. By volumetrically averaging these individual drainage volume pressures, the average pressure in the entire reservoir is obtained. To calculate the average pressure in each drainage volume, a correction is applied to the ordinary extrapolated pressure, i.e., the pressure obtained by extrapolating to infinite time the linear portion of the graph of closed-in pressure versus log[?t/(t + ?t)], where At is the closed-in time and t the production time. The correction, which is a function of the production time, is presented in graphical form for different shapes of the drainage area (horizontal cross section of the drainage volume). INTRODUCTION It is important to be able to find the volumetric average pressure in a reservoir so that the size of the reservoir may be determined from material balance calculations. It is also desirable to be able to find the approximate distribution of pressure within a reservoir for detection of fluid movement. 'The purpose of this paper is to present a method for calculating both the average reservoir pressure and the approximate distribution of pressure within a bounded reservoir — that is, a reservoir with no water drive. In reservoirs where the pressure builds up rapidly after wells are shut in, the determination of average pressure generally poses little problem, for one often need only average the final buildup pressures. It is when pressure buildup is slow that difficulties arise. For practical and economical reasons, the time allowable for closing in wells is limited. If at the maximum allowable closed-in time the pressure has not reached a con- stant value (and this is more often the case than is generally realized), calculation of average pressure presents difficulties. One method which has been used in the past for obtaining an "average" pressure involves the extrapolation to infinite closed-in time of curves of pressure as a function of closed-in time. Contour maps of these extrapolated pressures are used in arriving at an average pressure. In some cases wells are shut in for arbitrary Lengths of time, and contours of these arbitrary shut-in pressures are used for calculation. All these methods are somewhat arbitrary and suffer from lack of fundamental theoretical basis. For example, no precise meaning can be given either to the extrapolated or to the arbitrary shut-in pressure. Furthermore, there seems to be no reason why either set of pressures should be contourable. It was because of the vagueness of these old methods and their lack of a sound theoretical basis that, the present method for determining average reservoir pressure was developed. SYMBOLS A = drainage area of a well, cm' A, = drainage area of well j, cm2 A , = total drainage area of all wells in a reservoir, cm2 c = coefficient of compressibility of fluid at reservoir temperature, atm-1 f = hydrocarbin-filled porosity, fraction h = net thickness of formation, cm k = permeability of formation, darcy p = pressure at a well, atm p0 = initial reservoir pressure, atm p* = pressure obtained by straight-line extrapolation of linear portion of plot of p versus log [At/(t+ At)] p = volumetric average pressure inside drainage volume of a well q = volume flow rate at a well, cc/sec a: prevailing reservoir conditions q, = total volumetric flow rate from a reservoir, cc/sec at prevailing reservoir conditions rw = radius of well bore, cm T = kt/fµcA t = corrected time of production of a well, seconds, defined by

Jan 1, 1955

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 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 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 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 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 M. Prats

Most of the information available on the heat efficiency of hot fluid injection processes, both water and steam, has been obtained from calculated temperature distributions in the pay zone and adjacent formations. The usual approach has been to write the heat balance equations in terms of the temperatures, and then to introduce whatever simplifications are necessary to help obtain an analytical or a numerical solution. In either case, the assumption has often been made that the vertical thermal conductivity in the flooded interval is infinite, so that the temperatures within the flooded interval are then independent of vertical position. Because it was first made by Lauwerier,l and has been used extensively since, this assumption will be referred to as the Lauwerier assumption. Excellent reviews of the literature on the heat efficiency of hot fluid injection processes have been given by Spillette,2 Flock et al.,3 and Ramney.4 Of course, the most significant contributions take into account the effect of a finite vertical thermal conductivity on the vertical temperature profile. The most general analytic expression for the heat efficiency of a hot fluid injection process is that of Antimirov5 who considers injection of a heated incompressible fluid into a reservoir through an arbitrary number of wells. The rate of heat injection into the reservoir, as well as the injection temperature, is an arbitrary function of time. The geometry of the horizontal flow is arbitrary, although the reservoir is considered to be of uniform thickness and properties and to be of infinite areal extent. Heat transfer within the reservoir is by horizontal convection and conduc- tion, and by vertical conduction. In the formations adjacent to the flooded interval, heat transfer is by conduction in any direction. Thus, the available expression for the heat efficiency due to hot water injection has been developed for rather general conditions. This degree of applicability has not been obtained for other thermal recovery processes. By making the Lauwerier assumption we have been able to obtain expressions for the heat efficiency that are essentially independent of the thermal recovery process, be it steam, hot-water, or underground combustion. Clearly, then, the approach followed here is more restrictive than that of Antimirov5 in that it uses the Lauwerier assumption. At the same time less restrictive assumptions are made about the horizontal heat transfer mechanisms, either in the flooded interval or in the formations adjacent to it. The main difference between our approach and that of other investigators is that the heat balance in the pay zone is expressed for the pay zone as a whole rather than for a volume element within it. The Lauwerier assumption is introduced after the problem is developed along these lines as far as possible. The details of the development of the heat efficiency are given in the Appendix, while the main features of the model are discussed in the next section. Results and their implications are discussed later in a separate section. The Model Certain assumptions have been made to determine

Jan 1, 1970

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

By E. H. Herron, J. H. Perry

Mathematical simulation of reservoir behavior may be used to help understand reservoir processes and to predict reservoir behavior, thereby leading to the most economically desirable form of exploitation. In addition, simulation can be used as a tool for reservoir description to learn more about the physical nature of the reservoir and the mode of primary recovery. This use is essential in most reservoir studies and represents one of the more significant applications of simulation. Prediction of reservoir behavior provides information concerning displacement efficiency, optimum well locations, and the comparison of alternative producing processes. Since oil, water and gas all commonly occur in many reservoirs, a simulation that takes into account the simultaneous flow of three immiscible phases is a prerequisite for obtaining such information. Furthermore, a simulation that describes this flow in two dimensions provides the capability of evaluating combination gas- and water-drive reservoirs with either an areal or a cross-sectional description. Many reservoir problems require such a simulation. The objective of the research described here was to develop a mathematical model and associated computer program for accurate, efficient, and economical prediction of the reservoir flow of three phases in two-dimensional geometry. Several authors have discussed multidimensional reservoir simulation of two phases.l-3 Gottfried et al.' presented an analy- sis of three-phase simulation in one dimension. Fagin and Stewart" and Garrett6 have discussed the simulation of three-phase flow in two dimensions, but their solution techniques were different from the method discussed in this paper. In particular, previous models and computer programs have usually involved either fewer capabilities or inefficient and less rigorous descriptions resulting in limited applicability. In the following sections, we present a description of the model, an evaluation of the method, and the application to reservoir problems. Additional and supporting material is presented in the Appendix. The Three-Phase Model General Description The three-phase model consists of the mathematical description of the flow of oil, water and gas in a reservoir. Conceptually this description considers flow in all three space dimensions, but in the subject simulator it is assumed that no flow occurs in the third dimension. The model includes the effects on reservoir behavior of fluid and rock compressibility, viscosity, gravity, capillary pressure, relative permeability and gas solubility. Within the two-dimensional restriction, the reservoir considered may be completely heterogeneous and anisotropic. The differential equations governing multiphase, multidimensional flow in a reservoir are too complex to be solved analytically, so they are approximated by a system of finite difference equations and solved

Jan 1, 1970

By B. P. Nutter, M. Lebourg, J. A. McAlister

The Multi-Flow Evaluator (MFE) is a new system of tools providing an original approach in drill-stem testing. It improves control during the testing operation and gives a more accurate evaluation of the fluid recovered while providing additional pressure information for reservoir analysis. The tools are operated entirely by up and down motion of the drill stem. This up and down operation provides a positive means of control and offers easily observed surface indications of tool operating position. An unlimited number of shut-in and flow periods may be taken with this tool while in the hole. The key to the success of the MFE system of tools is the safety seal packer. Until this development, operation of the tools by up and down motion only had proven unreliable. The success ratio of the MFE now exceeds the ratio achieved by conventional tools. The equipment includes a 2,750-cc chamber in which a representative sample of the flowing formation fluid is trapped at the end of the last flow period and brought to the surface under pressure. The sample can be evaluated at the wellsite or transferred under pressure for laboratory analysis. The sampling feature of the tool allows a sample to be obtained from a reservoir which has suffered minimum influence from production. A representative sample can be obtained for laboratory or empirical analysis by employing a testing technique to minimize drawdown during flow periods. Interpretation methods to take advantage of this additional information are presented and supported by actual field examples. INTRODUCTlON The Multi-Flow Evaluator (MFE) represents the most recent technical advance in formation testing. Drill-stem testing techniques prior to the MFE have been limited in three areas—(1) operaiion, (2) recovery analysis, and (3) reservoir analysis. The tools in the MFE system allow maximum control during testing and provide test information superior to that obtainable with conventional tools. The system of tools and the operation procedures are described to illustrate the mechanical advantages over present tools. Methods of pressure build-up and fluid

Jan 1, 1966

By R. F. Krueger, L. C. Vogel

A design is presented for an electrical analog computer which can solve non-steady state heat transfer problems in an extensive radial forrnation containing a moving cylindrical source. The computer is used to simulate a .simplified thermal oil recovery process in which heat transfer from a moving, constant temperature source takes place radially by conduction only. Temperature distribution curves are .shown for several different assumed modes of travel for the heat source. The data are used to estimate the residual fuel requirements necessary to maintain a self-propagating isotherma1 front for the particular system being studied. Although the computer is designed to represent an unique system. conversion factors can he adjusted to show the effects of changing the assumed values of system constants, such as thermal conductivity of the formation and temperature of the source isotherm. Examples are given to illustrate the relative influence the assumed values of these variab1es will have upon the quantity of fuel required to maintain the source. The maximum effect of heat transfer by the flowing air stream is estimated mathematically, and full requirements for this system are shown to he markedly reduced. Inasmuch as the examples discussrd refer to a simplified process, the results are not directly applicable to a practical field operation. It is believed, however, that the data il1ustr-ate some general trends which are important in thermal recovery processes. Data from laboratory and field experimentation can he used to modify the computer to take into account the influence of heat transfer by various mechanisms, for example, by injected air and by fluids ahead of the hot zone. With .such refinements it would be possible to estimate more accurately the limiting values for the rate of travel of the high temperature front and for the required air- injection rates. Use of this computer should aid in evaluating the economic feasibility of oil recovery by various thermal processes. It is hoped this paper will stimulate further work by others to help accomplish this objective. As the costs of exploration for and development of new oil fields rise. increasing the recovery of oil from established fields becomes more and more important. In fields containing low gravity crude oil, ultimate recoveries by primary means are sometimes less than 10 per cent and often total only 10 to 20 per cent of the oil in place. Conventional secondary recovery methods are not always effective in increasing this total. Methods which have been proposed for increasing the recovery of viscous oils include thermal processes for increasing the mobility of the oil. One method is to provide a moving heat source generated within the formation to heat the oil and a flowing gaseous medium (for example, air) to sweep the oil to the producing wells. The possibility of such a process has intrigued engineers for many years.'. '." but little or no information has been available in the literature relative to the technical or economic feasibility of this method. Recently, however, interest has been revived by publication of laboratory and field test results indicating the feasibility of maintaining

Jan 1, 1956